工程井间地震波正演模拟及层析成像技术研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
井间地震由于具有多波接收、高精度、高分辨率、能量传播距离短、接近探测目标、避开低速带等一系列优点,在工程中得到了越来越广泛的应用。但是由于井间地震波的多样性和复杂性,为了更有效了解不同波场分量在均匀半空间的分布情况,分析井间波场中的反射、透射、绕射、散射以及能量的衰减等运动学和动力学的各种细节特征,更好地指导地质解释,本论文开展了基于典型工程地质模型的井间地震波的交错网格正演模拟;并在作者多年从事水利水电工程地震层析成像技术研究和应用的基础上,总结了工程应用中的一些重要技术措施,可为工程井间地震层析成像技术的推广应用,提供了科学依据,具有指导意义。本论文就工程井间地震波正演和层析成像技术两大内容,开展了以下工作:
     (1)在水利水电工程勘察中进行井间地震层析成像时,不同于油气田的井间地震勘探,其井深度较小,大多在100米深度内,当震源、检波器接近地表或井中水位面顶部时,地表面和井中水位面的影响不可忽略。为此,作者在本文中,结合水利水电工程的特点,建立了包含地表和井水位面的井间地震模型,并对其进行了正演模拟,得到了模拟的井间地震记录和水平粒子分量波场快照。
     (2)作者系统地总结了前人的研究成果及当前主要的地震波场正演算法的基础上,推导了二、三维的一阶速度-应力弹性波方程及其任意阶时间和空间精度的交错网格差分格式。在正演计算过程中,为节约内存空间和模拟计算时间,作者采用了非规则的可变有限差分网格和变化的时间步长理论,在保证计算精度的前提下,较大地提高了计算速度。
     (3)详细地探讨了交错网格差分法正演计算时的稳定性条件,得出了偶阶精度差分的时间和空间采样步长须满足的稳定性条件;分析了地震波数值频散产生的原因,阐述了通量校正传输(FCT)和提高差分精度两种频散压制方法;归纳了几种经典吸收边界条件的分类方法及这些吸收边界条件的优缺点,着重介绍了完全匹配层(PML)吸收边界条件,并导出了二、三维各向同性介质一阶应力-速度弹性波方程的相应的PML边界条件的差分格式,为跨孔地震波场的差分法正演模拟及其相应的程序编制打下基础。
     (4)结合Visual C++可视化程序设计语言,编制了井间交错网格的差分法地震波正演程序,把该自制程序应用于典型工程地质模型的正演计算中,得到了相应模型的波场快照和地震记录,记录中可以清晰观测到初至波、反射波、透射波及多次波的运动学和动力学特征。通过分析这些正演记录和波场快照,可更深刻地理解地震波场的传播规律,可加深对井间地震波场快照的认识,提高井间地震资料处理的解释精度。此外,对不同力源(水平方向单力源与水平方向力偶源)在快、慢地层产生的地震波,进行了全波场正演模拟,得到了压力、水平应力、水平粒子速度、垂直粒子速度的地震记录及波场快照,通过这些跨孔地震全波场的模拟,可全面把握不同场源、不同粒子速度的特征。
     (5)详细分析了走时(速度)层析成像技术的原理,总结了不同层析成像射线追踪技术的优缺点,并应用ART算法和SIRT通过对数字化模型分别进行重构,比较了两种算法的优劣。在此基础上编制了层析成像正反演处理系统软件,通过对四侧和两侧观测系统下的正方形缺陷模型和十字裂隙模型进行正演模拟和反演成像,深入地研究了观测系统对成像精度的影响。综合利用地质地球物理资料、速度测井、孔内摄像和地震记录,提出了建立初始速度模型和反演边界约束条件的综合方法,耿得了更准确和清晰的反演成像效果。然后,深入探讨了地震波吸收层析成像的机理,得知地震吸收层析成像方法,它可为地震波速层析成像提供有益补充。针对这两种层析成像技术的优缺点,提出了利用多参数(波速、吸收系数)综合地震层析成像的想法,使波速层析成像与吸收层析成像两者能相互印证、互为补充。作者把它们应用于工程实践中,有效地提高地震层析成像技术在工程应用中的解释精度,为地质解释与工程诊断提供更充分的依据。
Cross-hole seismic has been applied in project extensively because of many advantages, such as multi-wave receiving, high accuracy, high resolution, short of energy propagation length, approach detecting object, avoiding low speed zones. But because of the complexity and diversity of cross-hole seismic wave, for understanding the distributing instance of different wave field component in homogeneous half space effectively, analysis diversified detailed characteristics of kinematics and dynamics of reflection, transmission, diffraction, scattering and power attenuation in cross-hole seismic, and for better geological explanations, this text develops the forward simulation based on the cross-hole seismic wave equation and the seismic tomography technology. Based on the long time research and application of seismic tomography technology in water resources and hydropower project, the author sums up some important technology measures in project application, which could provide scientific basis and guiding meaning for the development and application of tomography technology of project cross-hole seismic. For the forward simulation and tomography technology of project cross-hole seismic wave, this text develops the below works:
     (1)When doing the cross-hole seismic tomography in water resources and hydropower project, the depth of the hole is small, and most of them are less than 100 meters. when the seismic source, wave detector are close to the ground or the top of water level surface in the hole, the influence of ground and the water level surface could not be ignored. So, in this text, combining with the characteristic of water resources and hydropower project, the author establishes cross-hole seismic model including ground and water level surface in the hole, does the forward simulation to it, calculate and completely analysis the propagation rule and characteristic of wave field in different field source and different particle velocity.
     (2)Based on systematically summing up the previous author's research and the mainly seismic wave field forward calculation method, the author derives first-order coupled elastic equations of 2D and 3D particle velocity and stress and the any order accurate time and space formulation of staggered-grid, finite-difference scheme. During the forward simulation calculation, the author takes irregular variable finite grids and variable time step theory. This would improve the calculation speed while guarantee the accuracy.
     (3)The author discusses detail the stability condition of staggered-grid, finite-difference method, and gets the stability condition which the time and space sampling step of high order accuracy difference must be satisfied. Analysis the reason of seismic wave numerical dispersion, and discusses the two dispersion compacting methods which are FCT and improving finite-difference accuracy. Conclude the classify method of some classical absorbing boundary condition and the advantages and disadvantages of them, then mainly introduces the basic theory used in forward simulation of PML absorbing boundary condition, and derives 2D or 3D PML absorbing boundary condition of first-order coupled elastic equations of particle velocity and stress in isotropic media. This could set the foundation for the cross-hole seismic wave field and the according program compiling.
     (4)Combining with Visual C++ visual program design language, compiles difference seismic wave forward program of cross-hole staggered-grid. Then apply this program into the forward simulation calculation of typical project geological models, and get the wave field snapshot and seismic record of the according model. In these results, it could clearly observe the kinematics and dynamics characteristics of primary wave, reflection wave, transmission wave and multiple wave of reflection. By analysis these forward record and wave field snapshot, it could understand the propagation regulation of seismic wave field more profundity, improve the comprehension to cross-hole seismic wave field snapshot and provide the technology support to the inversion interpretation. Moreover, do the whole wave field forward simulation to the seismic wave produced in quick, slow formation of different force source (level direction single force source and level direction double force source). It could get the seismic record and wave field snapshot of pressure, level stress, level particle velocity and vertical particle velocity. By the whole wave field simulation of this cross-hole seismic, it could grasp the characteristic of different field source and particle velocity completely.
     (5)analysis in detail the theory of travel-time tomography, sum up the advantages and disadvantages of different tomography ray tracing technology, and compare the advantages and disadvantages of ART and SIRT by applying them to rebuilding the numerical model separately. Based on this, the author compiles the inversion tomography process systematic software. By the forward simulation and inverse imaging of square defect model and cross fracture model in four sides and two sides observing system, it could research the influence of observing system to the imaging accuracy deeply. Comprehensively utilizing the geologic geophysical data, velocity logging, hole imaginary and seismic record, put forward the comprehensive method of setting up the initial velocity model, and gets the more accuracy and clearly inversion imaging effect.
     Then, deeply discuss the mechanism of seismic wave absorbing tomography; get the seismic absorbing tomography method, which could provide the useful supplement for the seismic wave absorbing tomography. According to the advantages and disadvantages of these two tomography technology method, the author puts forward utilizing multi-parameter to integrate the seismic tomography, which could make the wave velocity tomography and absorption tomography to makeup and confirm each other. The author applies these to project practice, which could improve the interpretations accuracy of seismic tomography technology in project, and provide more abundant basis to geological interpretations and project diagnose.
引文
[1]汪兴旺,杨勤海,孙党生,等.岩溶探测中井间地震波层析成像的应用[J].物探与化探,2008,32(1):105-108
    [2]黄金莉,刘家琦,欧阳代俊,等.工程物探中地震层析成像的研究[J].地球物理学报,1994,37(5):682-687
    [3]孙党生,李洪涛,任晨虹,等.井间地震波CT技术在水库渗漏勘察中的应用[J].勘察科学技术,2000,128(5):57-59
    [4]M C Williams.井间地震成像:这种技术的全盛时期来到了吗?邵祝华译[J].石油物探译从,1997,(5):39-46
    [5]陈世军,刘洪,周建宇,等.井间地震技术的现状与展望[J].地球物理学进展,2003,18(3):524-529
    [6]张文波.井间地震交错网格高阶差分数值模拟及逆时偏移成像研究:[博士学位论文].长安:长安大学,2005
    [7]J W Rector.Cross-well methods:Where are we,where are we gong?Geophysics,1995,60(3):627-630
    [8]孔庆丰.井间地震波场数值模拟技术研究与应用[J].勘探地球物理进展,2006,29(05):333-336
    [9]薛云峰.混凝土防渗墒质量控制及检测技术研究:[博士学位论文].长沙:中南大学,2007
    [10]牟永光.三维复杂介质地震物理模拟[M].山东:石油大学出版社,2003
    [11]窦玉坛.基于正演方法的井间地震响应特征研究:[硕士学位论文].青岛:中国科学院海洋研究所,2007
    [12]周建宇.井间地震研究与应用:[博士学位论文].兰州:中国科学院兰州地质所,2002
    [13]裴正林.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟[J].石油物探,2005,44(4):308-315
    [14]韩文功,印兴耀,王兴谋,等.地震技术新进展[M].山东:石油大学出版社,2006
    [15]佘德平.数值模拟技术的发展现状[J].勘探地球物理进展,2003,26(5-6):407-412+422
    [16]Crampin S.Evaluation of anisotropy by shear-wave spilitting[J].Geophysics,1985,50(1):142-152
    [17]赵鸿儒,唐文榜,郭铁栓.超卢地震模型试验技术及应用[M].北京:石油工业出版社,1986
    [18]赵群,郝守玲,宋遐龄,等.井间定向接收的物理模型实验[J].勘察科学技术,2000,128(5):57-59
    [19]曹均,贺振华,黄德济,等.储层孔(裂)隙的物理模拟与超声波实验研究[J].地球物理学进展,2004,19(2):386-391
    [20]黎霞.浅层地震模型试验研究[J].矿冶工程,2006,26(2):14-17
    [21]曹旭东,狄帮让,许洪波.地震物理模拟系统设计技术[J].微计算机信息,2006,22(1-2):160-162
    [22]邵志刚,傅容珊,黄建华,等.线性粘弹介质中地震波场数值模拟[J].地球物理学进展,2007,22(4):1135-1141
    [23]孙建国.复杂地表条件下地球物理场数值模拟方法评述[J].世界地质,2007,26(3):345-362
    [24]褚春雷.非规则网格有限差分法声波方程地震正演模拟及逆时偏移:[硕士学位论文].青岛:中国海洋大学,2003
    [25]张永刚.地震波场数值模拟方法[J].石油物探,2003,42(2):143-148
    [26]傅旦丹,朱宏彰,何汉漪.复杂介质多波地震正演模拟的动力学射线追踪法[J].中国海上油气(地质),1999,13(5):328-333.
    [27]曹俊兴,严忠琼.地表波跨孔旅行时层析成像分辨率的估计[J].成都理工学院学报,1995,22(4):290-293
    [28]何樵登,熊维纲.应用地球物理教程[M].北京:地质出版社,1991
    [29]程冰洁,李小凡,徐大吉.复杂非均匀介质伪谱法波场数值模拟[J].石油物探,2007,46(1):1-19
    [30]张建中,陈世军,余大祥.最短路径射线追踪方法及其改进[J].地球物理学进展,2003,18(1):146-150
    [31]陈景波,秦孟兆.辛几何算法在射线追踪中的应用[J].数值计算与计算机应用,2000,43(4):254-265
    [32]马在田.计算地球物理学概论[M].上海:同济大学出版社,1997
    [33]宋常瑜,裴正林.井间地震粘弹性波场特征的数值模拟研究[J].石油物探,2006,45(5):513
    [34]Zou Zhihui,Yu Wenhui.Wave field forward modeling and theoretical analysis of weakness in discrete media[J].Applied Geophysics,2006,3(2):75-81
    [35]Wu C L,Harris J M.Borehole seismic modeling with inclusion of tube wave and other tube-wave-related arrivals.73th Ann Intemat Mtg Soc Expl Geophys,2003:2239-2242
    [36]董良国.复杂地表条件下地震波传播数值模拟[J].勘探地球物理进展,2005,28(3):187-194
    [37]Yang D H,Lu M,Wu R S,et al.An optimal nearly analytic discrete method for 2D acoustic and elastic wave equations[J].Bull Seis Soc Am,2004,94:1982-1991
    [38]Fomberg B.The pseudospectral mthod:Comparisons with finite difference for the elastic wave equation[J].Geophysics,1987,52(4):483-501
    [39]李信富,李小凡,张美根.地震波数值模拟方法研究综述[J].防灾减灾工程学报,2007,27(2):241-248
    [40]郝奇.三维TTI介质地震波场正演模拟:[硕士学位论文].吉林:吉林大学,2006
    [41]杜启振,刘莲莲,孙晶波.各向异性黏弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报,2004,56(12):6143-6148
    [42]Carcione J M.Modeling an elastic singular surface wave in the earth.Geophysics,1992,57: 781-792
    [43]王秀明,GSERIANI,林伟军.利用谱元法计算弹性波场的若干理论问题[J].中国科学(G辑:理学力学天文学),2007,37(1):41-59
    [44]王润秋,耿伟峰,王尚旭.地震波正演模拟的任意差分精细积分方法[J].石油地球物理勘探,2003,38(03):252-257
    [45]李景叶,陈小宏.横向各向同性介质地震波场数值模拟研究[J].地球物理学进展,2006,21(3):700-705
    [46]刘大利,陈磊,桂冰.二维交错网格高分辨格式的并行实现[J].南京林业大学学报(自然科学版),2006,30(5):71-75
    [47]陈建忠,封建湖,史忠科,等.求解双曲型守恒律的半离散中心差分格式[J].应用力学学报,2005,22(4):536-540
    [48]奚先,姚姚.二维弹性随机介质中的波场特征[J].地球物理学进展,2005,20(1):147-154
    [49]孙卫涛,杨慧珠.双相各向异性介质弹性波场有限差分正演模拟[J].固体力学学报,2004,25(1):21-28
    [50]杨宽德,杨顶辉,王书强.基于BISQ高频极限方程的交错网格法数值模拟[J].石油地球物理勘探,2002,37(5):463-468
    [51]张文生,宋海斌.三维正交各向异性介质三分量高精度有限差分正演模拟[J].石油地球物理勘探,2001,36(4):422-432
    [52]Falk J,Tessmer E,Gajewski D.Tube wave modeling by the finite-difference method with varying grid spacings.Pageoph,1996,148:77-93
    [53]Pilant W L.Elastic Waves in the Earth.Amsterdam:Elsevier Press,1979.59-63
    [54]White J E,Lessenger M A.Caliper effect on borehole coupling.Expl Geophys,1988,19:201-205
    [55]Mo L W,Harris J M.Analysis and attenuation of tube waves in crosswell seismic survey.65th Ann Internat Mtg Soc Expl Geophys,1995:438-441
    [56]陶果,张友生,张洪娥,等.用于声波测井的大型三维有限差分模拟程序[J].测井技术,2001,25(4):273-277
    [57]Tam C K W,Webb J C.Dispersion-relation-preserving finite difference schemes for computational acoustics.J Comput Phys,1993,107:262-281
    [58]Park C B,Miller R D,Xia J.Multichannel analysis of surface waves[J].Geophysics,1999,64:800-808
    [59]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,(3):411-419
    [60]裴正林.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟[J].中国石油大学学报(自然科学版),2006,30(2):16-20
    [61]裴正林.双相各向异性介质弹性波传播交错网格高阶有限差分法模拟[J].石油地球物理勘探,2006,41(2):137-143
    [62]裴正林,王尚旭.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟[J].地震学报,2005,27(4):441-451
    [63]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5):23-29
    [64]裴正林.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟[J].石油地球物理勘探,2004,39(6):629-634
    [65]裴正林,牟永光.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),2003,27(6):17-21
    [66]Xia Fan,Dong Liangguo,Ma Zaitian.The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid,Finite Difference Scheme[J].Applied Geophysics,2004,1(1):38-41
    [67]Moczo P.Finite-difference technique for SH-wave in 2D media using irregular grids-application to the seismic response problem.Geophys J Int,1989,99:321-329
    [68]Robert W Graves.Simulating Seismic Wave Propagation in 3D Elastic Media Using Staggered-Grid Finite Differences[J].Bulletin of Seismological Society of America,1996,86(4):1091-1106
    [69]Rune Mittet.Free-surface boundary conditions for elastic staggered-grid modeling schemes [J].Geophysics,2002,67:1616-1623
    [70]Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress f'mite-difference method[J].Geophysics,1986,51(4):889-901
    [71]Carcione J M,G C Herman,A P E ten Kroode.Seismic modeling[J].Geophysics,2002,67:1304-1325
    [72]Graves R W.Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences[J].Bull.Seism.Soc.Am,1996,86:1091-1106
    [73]Lindman E L.Free space boundary conditions for the time dependent wave equation[J].J.Comp.Phys.1975,18:66-78
    [74]Robertsson J O A.A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J].Geophsics,1996,61:1921-1934
    [75]Francis Collino,Chrysoula Tsoka.Applicition of the perfectly matched absorbing layer model to the liner elastodynamic problem in anistropic heterogeneous media[J].Geophysics,2001,66:294-307
    [76]常君勇.交错网格有限差分法地震数值模拟:[硕士学位论文].北京:中国地质大学,2007
    [77]董良国,李培明.地震波传播数值模拟中的频散问题[J].天然气工业,2004,24(6):53-56
    [78]李国平.组分k型黏弹性介质及其波场特征研究:[硕士学位论文].北京:中国地质大学,2007
    [79]张文忠.Biot介质的交错网格差分法波场模拟研究:[硕士学位论文].北京:中国地质大学,2007
    [80]Falk J,Tessmer E,Gajewski D.Tube wave modeling by the finite-difference method with varying grid spacings[J].Pageoph,1996,148:77-93
    [81]Ivo Oprsal,Jirl Zahradnik.Elastic finitic-difference method for irregular grids[J]. Geophysics,1999,64(1),240-250
    [82]王春镛,李幼铭.地震层析成像及其应用[M].北京学术书刊出版社,1989
    [83]杨国宪,张可霓,李丕武.层析成像在黄河小浪底工程地下厂房区域的应用研究[J].地球物理学报,1993,36(1):106-112
    [84]杨文采.应用于地质勘探和矿产开发的地表层析成像[J].中国地质,1991,(7):22-26
    [85]杨文采,杜剑渊.煤矿陷落柱地震层析成像试验[J].物探与化探,1992,16(3):174-180
    [86]杨文采,杜剑渊.层析成像新算法及其在工程检测中的应用[J].地球物理学报,1994,37(1):239-244
    [87]王旭,陈志强,熊华,等.联合代数重建算法中基于像素的投影计算方法[J].核电子学与探测技术,2005,25(6):784-788
    [88]杨文采.地球物理反演的理论与方法[M].北京:地质出版社,1997
    [89]魏亦文.基于复杂介质最短路径射线追踪层析成像正演研究:[硕士学位论文].长沙:中南大学,2007
    [90]刘永华,田宗勇,喻振华,等.工程VSP与地震CT联合探测方法及其在岩土工程的应用[J].地球物理学进展,2005,20(1):267-272
    [91]李志鹏,丛鹏,邬海峰.代数迭代算法进行CT图像重建的研究[J].核电子学与探测技术,2005,25(2):184-186
    [92]张顺利,张定华,熬波,等.ART算法高质量重建研究[J].核电子学与探测技术,2007,27(4):719-723
    [93]张云妹,顾汉明,师学明.基于MOSER曲线追踪的SIRT声波层析成像[J].工程地球物理学报,2005,2(3):167-176
    [94]宿淑春,王守东,吴律.井间层析成像的平滑SIRT算法[J].石油大学学报(自然科学版),2001,25(6):29-31
    [95]李尊营,宋一中等.复杂流场的代数迭代重建算法[J].光电子.激光,2006,17(10):1259-1263
    [96]刘劲松,刘福田,刘俊等.地震层析成像LSQR算法的并行化[J].地球物理学报,2006,49(2):540-545
    [97]于更新.有限差分模拟走时波场的地震层析成像方法与应用研究:[硕士学位论文].北京:中国地质大学,2006
    [98]张赛民.跨孔地震层析成像研究:[硕士学位论文].长沙:中南大学,2003
    [99]杨文采.地球物理反演和地震层析成像[M]北京:地质出版社,1989
    [100]赵永贵,李勤,郭鸿,等.频率域地震吸收CT方法及其工程应用[J].中国科学,2002,30(1):107-112
    [101]杜世通.地震波动力学[M].山东:石油大学出版社,2003
    [102]孙成禹.地震波理论与方法[M].山东:石油大学出版社,2003
    [103]王春燕.高阶交错网格有限差分地震波场计算:[硕士学位论文].成都:成都理工大学,2007
    [104]Levander A R.Fourth-order finite-difference P-SV seismograms[J].Geophsics,1988,53:1425-1436
    [105]Groenenboom J,Falk J.Scattering by hydraulic fractures:Finite-difference modeling and laboratory data.Geophysics,2000,65:612-622
    [106]Moczo P,Bystricky E,Kristek J,et al.Hybrid modeling of P-SV seismic motion at inhomogeneous viselastic topographic structures.Bull Seis Soc Am,1997,87:1305-1323
    [107]Yang D H,Peng J M,Lu M,et al.Optimal nearly analytic discrete approximation to the scalar wave equation.Bull Seis Soc Am,2006,96:1114-1130
    [108]周竹生,刘喜亮,熊孝雨.弹性介质中瑞雷面波有限差分法止演模拟[J].地球物理学报,2007,50(2):567-573
    [109]长文忠.Biot介质的交错网格差分法波场模拟研究:[博士学位论文].北京:中国地质大学,2007
    [110]Emmerich,H.and M.korn,Incorporation of attenuation into time-domain computations of seismic wave fields[J].Geophysics,1987,52:1252-1264
    [111]Forbriger T.Inversion of shallow-seismic wavefields:1.Wavefield transformation [J].Geophysics J Int,2003,153:719-734
    [112]刘喜亮.基于波动方程的瑞利面波有限差分法正演模拟:[硕士学位论文].长沙:中南大学,2006
    [113]刘洋,李承楚,牟永光.任意偶数阶精度有限差分法数值模拟[J].石油地球物理勘探,1998,33(01):1-10
    [114]熊章强.复杂介质中瑞雷面波的正演模拟及传播特征研究:[博士学位论文].武汉:中国地质大学,2006
    [115]董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶筹分解法稳定性研究[J].地球物理学报,2000,43(6):856-864
    [116]轩义华.倾斜横向各向同性介质_TTI_参数反演方法研究:[博士学位论文].吉林:吉林大学,2007
    [117]Day S M,J B Minster.Numerical simulation of wave fields using a Pade approximant method[J].Geophysics,1984,78:105-118
    [118]Joachim Falk.Efficient finite-difference modeling of seismic waves using locally adjustable time steps[J].Geophysical prospecting,1998,46:603-616
    [119]Dablain M A.The application of high-order differencing to the scalar wave equation[J]..Geophysics,1986,51:54-66
    [120]Fornberg B.Generation of finite difference formulas on arbitrary spaced grids[J].Math.Computation,1988,51:699-706
    [121]Oprsal I,Zahradnik J.Elastic finite-difference method for irregular grids.Geophysics,1999,64:240-250
    [122]Klimes L.Accuracy of elastic finite differences in smooth media[J].Pageoph,1996,148,39-76
    [123]Bohlen T,Saenger E H.Accuracy of heterogeneous staggered-grid finite-difference modeling of Rayleigh waves[J].Geophysics,2006,71(4):109-115
    [124]Saenger,E.H.and T.Bohlen.Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid[A].Geophysics 69,2004,583-591.
    [125]Moczo P,J Kristek,V Vavrycuk,et al.3D heterogeneous staggered-grid finite-difference modeling of seismic motion with volume harmonic and arithmetic averaging of elastic moduli and densities[J].Bull.Seism.Soc.Am.,2002,92:3042-3066.
    [126]Graves R W.Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences[J].Bull Seism Soc Am,1996,86:1091-1106
    [127]Saenger E H,Gold N,Shapiro S A.Modeling the propagation of elastic waves using a modified finite-difference grid[J].Wave Motion,2000,31:77-92
    [128]奚先,姚姚.二维粘弹性随机介质中的波场特征分析[J].地球物理学进展,2004,19(3):608-615
    [129]孙卫涛,杨慧珠.各向异性介质弹性波传播的三维不规则网格有限差分方法[J].地球物理学报,2004,47(2):332-337
    [130]王秀明,张海澜,王东.利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播[J].地球物理学报,2003,46(6):842-849
    [131]孟庆生,何樵登,朱建伟,王德利.基于BISQ模型双相各向同性介质中地震波数值模拟[J].吉林大学学报(地球科学版),2003,33(2):227-221
    [132]赵福云,汤广发,周安伟.基于交错网格的高效离散格式[J].湖南大学学报(自然科学版),2002,29(3):147-150
    [133]杨宽德,杨顶辉,王书强.基于横向各向同性BISQ方程的弹性波传播数值模拟[J].地震学报,2002,24(6):599-606
    [134]胡勇,钟兵,杨雅和.用交错网格建立气水两相数值试井模型[J].油气井测试,2001,10(3):10-14
    [135]邱建贤,戴嘉尊.一类求解Hamilton-Jacobi方程的交错网格差分格式[J].南京航空航天大学学报,2000,32(5):573-578
    [136]李本文,赫冀成,刘日新.不同差分格式在同位网格系统中的计算效果比较[J].东北大学学报(自然科学版),2000,21(4):412-415
    [137]杜启振,刘莲莲,孙晶波.各向异性粘弹性孔隙介质地震波场伪谱法正演模拟[J].物理学报,2007,56(10):6143-6149
    [138]王德利,雍运动,韩立国,等.三维粘弹介质地震波场有限差分并行模拟[J].西北地震学报,2007,29(1):30-34
    [139]程冰洁,李小凡,徐天吉.复杂非均匀介质伪谱法波场数值模拟[J].石油物探,2007,29(1):30-34
    [140]轩义华,何樵登,林炎,等.双相各向异性介质正演研究[J].石油天然气学报(江汉石油学院学报),2006,28(6):69-72
    [141]熊晓军,贺振华,黄德济.复杂地表条件下的地震波场模拟方法研究[J].西安石油大学学报(自然科学版),2006,21(6):5-7
    [142]轩义华,何樵登,孟庆生,等.基于BISQ机制的双相EDA介质的波场分析[J].石油地球物理勘探,2006,41(5):550-556
    [143]熊煜,李录明,罗省贤.各向异性介质弹性波场正演及偏移[J].成都理工大学学报(自 然科学版),2006,33(3):310-316
    [144]熊晓军,贺振华,黄德济.三维波动方程正演及模型应用研究[J].石油物探,2005,44(6):554-556
    [145]杜正聪,贺振华,黄德济.基于混合方程的裂隙介质地震波场数值模拟[J].石油天然气学报,2005,27(2):197-200
    [146]郑海山,张中杰.横向各向同性(VTI)介质中非线性地震波场模拟[J].地球物理学报,2005,48(3):660-671
    [147]杜正聪,贺振华,黄德济,等.基于弹性方程与声学方程的裂缝模型正演与偏移[J].物探化探计算技术,2005,27(1):1-4
    [148]何凯涛,ZHANG Hai-jiang.小波多分辨分析方法在波动方程正演模拟中的应用[J].地球物理学进展,2001,16(2):73-80
    [149]阮爱国,李清河.非均匀各向异性弹性波场正演[J].西北地震学报,2001,23(4):318-329+338
    [150]孟凡顺,郭海燕,和转,王渝明.复杂地质体粘滞弹性波止演模拟的有限差分法[J].青岛海洋大学学报(自然科学版),2000,30(2):315-320
    [151]杜正聪,贺振华,黄德济,等.基于弹性方程与声学方程的裂缝模型正演与偏移[J].物探化探计算技术,2005,27(1):1-5
    [152]邵治龙,贺振华,黄德济.正交各向异性介质中地震记录的有限差分数值模拟[J].物探化探计算技术,1998,20(4):300-305
    [153]牟永光,裴正林.三维复杂介质地震数值模拟[M],北京:石油工业出版社,2005
    [154]李振春,张华,刘庆敏,等.弹性波交错网格高阶有限差分法波场分离数值模拟[J].石油地球物理勘探,2007,42(5):510-515
    [155]程冰洁,李小凡,徐天吉.非均匀介质中交错网格高阶有限著分数值模拟[J].物探化探计算技术,2006,28(4):294-298
    [156]孙卫涛,杨慧珠,舒继武.非均匀介质弹性波动方程的不规则网格有限差分方法[J].计算力学学报,2004,21(2):135-142
    [157]李胜军,孙成禹,倪长宽,等.声波方程有限差分数值模拟的变网格步长算法[J].工程地球物理学报,2007,4(3):207-212
    [158]王一博,杨慧珠.方位各向异性介质的多尺度有限差分法波场模拟[J].清华大学学报(自然科学版),2006,46(2):293-296
    [159]Wu C L,Harris J M.Single-well seismic modeling in viscoelastic media using a variable-grid finite-difference method.71th Ann Intemat Mtg Soc Expl Geophys,2001:1155-1158.
    [160]Aoi S,Fujiwara H.3D finite-difference method using discontinuous grids.Bull Seis Soc Am,1999,89:918-930
    [161]赵海波,王秀明.一种优化的交错变网格有限差分法及其在井间卢波中的应用[J].科学通报,2007,52(12):1387-1395
    [162]Wang Y,Xu J,Schuster G T.Viscoelastic wave simulation in basins by a variable-grid finite-difference method.Bull Seis Soc Am,2001,91:1741-1749
    [163]Mikumo T,Miyatake T.Dynamic rupture processes on a dipping fault,and estimates of stress drop and strength excess from the results of waveform inversion.Geophys J Int,1993,112:481-496
    [164]Pitarka A.3D elastic finite-difference modeling of seismic motion using staggered grids with nonuniform spacing.Bull Seis Soc Am,1999,89:54-68
    [165]Wu C L,Harris J M,Nihei K T,et al.Two-dimensional finite-difference seismic modeling of an open fluid-filled fracture:Comparison of thin-layer and linear-slip models.Geophys,2005,70:T57-T62
    [166]Hayashi K,Burns D,Toks?z M N.Discontinuous-grid finite-difference seismic modeling including surface topography.Bull Seis Soc Am,2001,91:1750-1764
    [167]Jastram C,Tessmer E.Elastic modeling on a grid with vertically varying spacing.Geophys Prosp,1994,42:357-370
    [168]Fornberg B.Generation of finite difference formulas on arbitrarily spaced grids.Math Comput,1988,51(184):699-706
    [169]Adriano S,Oliveira M.A fourth-order finite-difference method for the acoustic wave equation on irregular grids.Geophysics,2003,68:672-676
    [170]Ye F,Chu C L.Dispersion-relation-preserving finite difference operators:derivation and application.75th Ann Internat Mtg Soc Expl Geophys,2005:1783-1786
    [171]Charles Cerjian.A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50:705-708
    [172]Reynolds A C.Boundary conditions for the numerical solution of wave propagation problems[J].Geophysics,1978,42(6):109-1110
    [173]夏凡,董良国,马在田.三维弹性波数值模拟中的吸收边界条件[J].地球物理学报,2004,47(1):132-136
    [174]董良国.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探,1999,34(1):45-56
    [175]Liao Z P,Wong H L,Yang B P,et al.Atransmitting boundary for gansient wave analyses [J].Scientia Sinica(Sefies A),1984,27(10):1063-1076
    [176]赵海波,王秀明,王东,等.完全匹配层吸收边界在孔隙介质弹性波模拟中的应用[J].地球物理学报,2007,50(2):581-591
    [177]Wang T,Tang X M.Finite-difference modeling of elastic wave propagation:A nonsplitting perfectly matched layer approach[J].Geophysics,2003,68:1749-1755
    [178]Festa G,Nielsen S.PML absorbing boundaries[J].Bull Seism Soc Am,2003,93:891-903
    [179]Cerjan C,Kosloff D,Kosloff R,et al.A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50:705-708
    [180]Zeng Y Q,He J Q,Liu Q H.The application of the perfectly matched layer in numerical modeling of wave propagation in poroelastic media[J].Geophysics,2001,66:1258-1266
    [181]Robertsson J O A.A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J].Geophysics,1996,61:1921-1934.
    [182]Liao Z,H L Wong,Y Baipo,et al.A transmitting boundary for transient wave analysis[J].Scienfia Sinica,1984,27:1063-1076
    [183]Higdon R L.Absorbing boundary conditions for elastic waves[J].Geophysics,1991,56:231-241.
    [184]Clayton R,B Engquist.Absorbing boundary conditions for acoustic and elastic wave equations[J].Bull.Seism.Soc.Am.1977,67:1529-1540
    [185]Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].J.Comp.Phys,1994,114:185-200
    [186]Komatitsch D,J Tromp.A Perfectly Matched Layer(PML)absorbing condition for the second-order elastic wave equation[J].Geophys.J.Int,2003,154:146-153
    [187]Kosloff R,D Kosloff.Absorbing boundaries for wave propagation problems[J].J.Comp.Phys,1986,63:363-376
    [188]Chew W C,Q H Liu.Perfectly matched layers for elastodynamics:a new absorbing boundary condition[J].J.Comp.Acoustics,1996,4:341-359
    [189]Randall C J.Absorbing boundary condition for the elastic wave equation:velocity-stress formulation[J].Geophysics,1989,54:1141-1152
    [190]Abarbanel S,D.Gottlieb.A mathematical analysis of the PML method[J].Comput Phys,1997,134,357-363
    [191]Zhao H B,Wang X M,Zhang H L.Studies on effective and stable absorbing boundary conditions in ultrasonic wave modeling.Proc IEEE Ultrason Syrup,2005,3:1472-1475
    [192]Robertsson,.J.O.A.,A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J].Geophsics,1996,61,1921-1934.
    [193]Charles Cerjian,Dan Kosloff,Ronnie Kosloff,etc.A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics,1985,50,705-708
    [194]陈浩,王秀明,赵海波.旋转交错网格有限差分及其完全匹配层吸收边界条件[J].科学通报,2006,51(17):1985-1994
    [195]李景叶,陈小宏.TI介质地震波场数值模拟边界条件处理[J].西安石油大学学报(自然科学版),2006,21(4):20-23
    [196]Song Ruolong,Ma Jun,Wang Kexie.The Application of the Nonsplitting Perfectly Matched Layer in Numerical Modeling of Wave Propagation in Poroelastic Media[J].Applied Geophysics,2005,2(4):216-222
    [197]Collino F,Tsogka C.Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2001,66:294-307
    [198]杨三女,陈永清,陈建国.伪谱法求解波动方程算法的VC++实现[J].华尔地质学院学报,2003,26(2):183-186
    [199]李信富,李小凡.伪谱法地震波场数值模拟[J].桂林工学院学报,2007,27(5):178-181
    [200]成谷,马在田,耿建华等.地震层析成像发展回顾[J].勘探地球物理学进展,2002,25(3):6-12
    [201]裴正林.井间地震层析成像的现状与进展[J].地球物理学进展,2001,16(3):91-97
    [202]王辉,黄鼎成.地震层析成像方法及其在岩体结构研究中的应用[J].工程地质学报,2000,8(1):109-117
    [203]王运生.弯曲射线地表波透射层析成像的一种实现方法[J].河海大学学报,1993(40):227-230
    [204]左德堃等.地震衰减层析成像理论基础及其实用性探讨[J].煤炭工程师,1994,2:44-49+4
    [205]黄政宇,黄靓,缪仑.混凝土超声波幅值衰减层析成像研究[J].铁道科学与工程学报,2004,1(2):35-38.
    [206]何卫红.地震吸收系数层析成像方法及在工程中的应用[J].西部探矿工程,2004,96(5):4-6.
    [207]严又生,宜明理,魏新,等.井间地震速度和Q值联合层析成像及应用[J].石油地球物理勘探,2001,36(1):9-17.
    [208]赵永贵,李勤,郭鸿,等.频率域地震吸收CT方法及其工程应用[J].中国科学D辑,2000,30(1):107-112.
    [209]邢庆祝.地质CT成像技术中正演模拟技术的研究与应用:[硕士学位论文].福州:福州大学,2005
    [210]吴律.层析基础及其在井间地震中的应用[M].北京:石油工业出版社,1997
    [211]杨文采.应用地震层析成像[M]北京:地质出版社,1993
    [212]马昭军,刘洋.地震波衰减反演研究综述[J].地球物理学进展,2005,20(4):1074-1082
    [213]杨文采,李幼铭.应用地震层析成像[M].北京:地质出版社,1993
    [214]杨文采.地球物理反演的理论与方法[M].北京:地质出版社,1993
    [215]陈仲侯,王兴泰,杜世汉.工程与环境物探教程[M].北京:地质出版社,1993
    [216]Guust Nolet.地震层析技术[M].北京:地质出版社,1991
    [217]G.T.赫尔曼,H.k图伊,K.J.兰振伯格,等.层析成像和反演问题的基本方法[M].北京:石油工业出版社,1997
    [218]骆循.ART算法及其在地球物理反演中的作用[J].物探化探计算技术,1987,9(1):36-38
    [219]梁国钱,苏全.层析成像技术在拱坝质量检测中的应用[J].水利发电,2002,(5):20-22.
    [220]李盛东,李月华.地震走时层析成像算法与比较[J].中国矿业大学学报,2000,29(2):213-214
    [221]曹小林,严忠琼,朱介寿.井间地震层析成像技术及其在岩溶勘察中的应用[J].地质灾害与环境保护,1997,8(2):53-58.
    [222]孙跃军,俞国柱,石桂.井间地震层析成像技术及其在岩溶勘察中的应用[J].CT理论与应用研究,2001,10(4):10-13
    [223]王超凡,赵永贵,郭鸿等.地震吸收CT方法及其在采空区探测中的应用[J].地球物理学报,1998,41(增刊):367-375
    [224]苏全.井间层析成像技术在水利结构工程检测中的应用:[硕士学位论文].淅江:淅江大学,2007
    [225]喻振华,冯德山,汤井田.地震波速与吸收系数综合层析成像[J].地球物理学进展, 2006,21(4):1338-1341
    [226]赵连锋.井间地震波速与衰减联合层析成像方法研究:[博士学位论文].成都:成都理工大学,2002
    [227]陈培善.地震波衰减研究在我国的进展[J].地球物理学报,1994,37(S1):231-240
    [228]王辉,常旭,高峰等.井间地震波衰减成像的几种方法[J].地球物理学进展,2001,16(1):104-109
    [229]Quan Y L,Jerry M H.Seismic attenuation tomography using the frequency shift method[J].Geophysic,1997,62(3):895-905
    [230]朱介寿,严忠琼,曹俊兴.井间地球物理层析成像软件系统研究[J].物探化探计算技术,1994,16(4):310-320
    [231]王家映.地球物理反演理论[M].北京:高等教育出版社,2002
    [232]杜剑渊,杨文采.跨孔地震层析成像的级联方法[J].地球物理学报,1991,31(6):441-443
    [233]杨辉,戴世坤,宋海斌等.综合地球物理联合反演综述[J].地球物理学进展,2002,17(2):262~271
    [234]金亚荣,陈国平.层析技术在岩溶地区岩土工程勘察中的应用[J].中山大学学报,2000,39(2):114-118

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700