频率域全波形反演方法研究进展
|
摘要
全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力.根据研究需要,全波形反演既可在时间域也可在频率域实现.频率域相对于时间域反演具有计算高效、数据选择灵活等优势.近十几年来频率域全波形反演理论在波场模拟方法、反演频率选择策略、目标函数设置方式、震源子波处理方式、梯度预处理方法等方面取得了进展.目标函数存在大量局部极值的特性是影响反射地震全波形反演效果的重要内在因素之一.如果将Laplace域波形反演、频率域阻尼波场反演、频率域波形反演三种方法有机结合,可以降低反演的非线性程度.
The full waveform inversion (FWI) method utilizes the kinematic and dynamic information of prestack seismic data to rebuild underground velocity structure.It has the potential of revealing detailed structure and lithology characteristics under complex geological background.FWI can be carried out in the either frequency or time domain. Frequency-domain FWI has the advantages of high computation efficiency and data selection flexibility over its time-domain counterpart.Recent advances in frequency-domain FWI including waveform modeling methods,frequency selection strategies,object function configuration styles,source wavelet processing methods,and gradient precondition methods make it attractive for realistic problems.The multiple local minimal nature of the object function js one of the major obstacles to good reflection waveform inversion results.The frequency-domain damped wavefield inversion method fills the gap between the Laplace-domain and frequency-domain FWI methods.The combination of the three methods can mitigate the strong nonlinearity of the waveform inversion problem.
The full waveform inversion (FWI) method utilizes the kinematic and dynamic information of prestack seismic data to rebuild underground velocity structure.It has the potential of revealing detailed structure and lithology characteristics under complex geological background.FWI can be carried out in the either frequency or time domain. Frequency-domain FWI has the advantages of high computation efficiency and data selection flexibility over its time-domain counterpart.Recent advances in frequency-domain FWI including waveform modeling methods,frequency selection strategies,object function configuration styles,source wavelet processing methods,and gradient precondition methods make it attractive for realistic problems.The multiple local minimal nature of the object function js one of the major obstacles to good reflection waveform inversion results.The frequency-domain damped wavefield inversion method fills the gap between the Laplace-domain and frequency-domain FWI methods.The combination of the three methods can mitigate the strong nonlinearity of the waveform inversion problem.
引文
[1] 宋海斌,张关泉.层状介质弹性参数反演问题研究综述[J].地球物理学进展,1996,13(4) :67～78． Song H B,Zhang G Q.A summary of inverse problem research on elastic parameters of stratified media[J].Progress in Geophys.(in Chinese),1996,13(4) : 67～78．
[2] 苑书金.叠前地震反演技术的进展及其在岩性油气藏勘探中的应用[J].地球物理学进展,2007,22(3) :879～885． Yuan S J.Progress of pre-stack inversion and application in exploration of the lithological reservoirs [ J ].Progress in Geophys.(in Chinese),2007,22(3) : 879～885．
[3] 杨文采.非线性波动方程地震反演的方法原理及问题[J].地球物理学进展,1992,7(1) :9～19． Yang W C.Methods and problems in seismic inversion of nonlinear wave equations [ J ].Progress in Geophys.(in Chinese),1992,7(1) :9～19．
[4] 杨文采.非线性地球物理反演方法:回顾与展望[J].地球物理学进展,2002,17(2) :255～261． Yang W C.Non-linear geophysical inversion methods: review and perspective[J].Progress in Geophys.(in Chinese),2002,17(2) : 255～261．
[5] Williamson P R.A guide to the limits of resolution imposed by scattering in ray tomography[J].Geophysics,1991,56(2) : 202～207．
[6] Williamson P R,Worthington M H.Resolution limits in ray tomography due to wave behavior: Numerical experiments[J].Geophysics,1993,58(5) :723～735．
[7] Wu R S,Toksoz M N.Diffraction tomography and multisource holography applied to seismic imaging.Geophysics,1987,52(1) :11～25．
[8] Pratt R G.Seismic waveform inversion in the frequency domain,Part 1: Theory and verification in a physical scale model[J].Geophysics,1999,64(3) :888～901．
[9] Pratt R G,Shipp R M.Seismic waveform inversion in the frequency domain,Part 2: Fault delineation in sediments using crosshole data[J].Geophysics,1999,64(3) :902～914．
[10] 曹辉,陈国金,郭建,等.频率域2． 5-D井间波形层析成像及其实际应用[J].地球物理学报,2008,51(3) :91 5～923． Cao H,Chen G J,Guo J,et al.The 2． 5-D cross-well waveform tomography in the frequency domain and its practical application[J].Chinese J.Geophys.(in Chinese),2008,51(3) :915～923．
[11] Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49 (8) : 1259 ～ 1266．
[12] Tarantola A.A strategy for nonlinear elastic inversion of seismic reflection data[J].Geophysics,1986,51 (10) : 1893 ～ 1903．
[13] Pratt R G.Frequency-domain elastic wave modeling by finite differences: A tool for crosshole seismic imaging [J].Geophysics,1990,55(5) :626～632．
[14] Pratt R G,Worthington M H.Inverse theory applied to multi-source cross-hole tomography,Part 1 : Acoustic waveequation method[J].Geophysical Prospecting,1990,38:287～ 310．
[15] Pratt R G.Inverse theory applied to multi source cross-hole tomography,Part 2: Elastic wave-equation method [ J ].Geophysical Prospecting,1990,38:311～329．
[16] Pratt R G,Shin C,Hicks G J.Gauss-Newton and full Newton methods in frequency domain seismic waveform inversion[J].Geophys.J.Int.,1998,133:341～362．
[17] Shin C,Cha Y H.Waveform inversion in the Laplace domain [J].Geophys.J.Int.,2008,173:922～931．
[18] Operto S,Virieux J,Dessa J X,et al.Crustal seismic imaging from multifold ocean bottom seismometer data by frequency domain full waveform tomography: Application to the eastern Nankai trough[J].J.Geophys.Res.,2006,111,B09306,doi: 10． 1029/2005JB003835．
[19] Bleibinhaus F,Hole J A,Ryberg T,et al.Structure of the California Coast Ranges and San Andreas Fault at SAFOD from seismic waveform inversion and reflection imaging[J].J.Geophys.Res.,2007,112,B06315,doi: 10． 1029/2006JB004611．
[20] Malinowski M,Operto S.Quantitative imaging of the Permo-Mesozoic complex and its basement by frequency domain waveform tomography of wide-aperture seismic data from the Polish Basin[J].Geophysical Prospecting,2008,56: 805～825．
[21] Gao F C,Levander A,Pratt R G,et al.Waveform tomography at a groundwater contamination site: VSP-surface data set[J].Geophysics,2007,71(1) :H1～H11．
[22] Gao F C,Levander A,Pratt R G,et al.Waveform tomography at a groundwater contamination site: Surface reflection data[J].Geophysics,2007,72(5) :G45～G55．
[23] Ben-Hadj-Ali H,Operto S,Virieux J.Velocity model building by 3D frequency-domain,full-waveform inversion of wide-aperture seismic data [J].Geophysics,2008,73 (5) :VE101～VE117．
[24] Hicks G J,Pratt R G.Reflection waveform inversion using local descent methods Estimating attenuation and velocity over a gas-sand deposit[J].Geophysics,2001,66(2) : 598～612．
[25] 许琨,王妙月.声波方程频率域有限元参数反演[J].地球物理学报,2001,44(6) ,852～864． Xu K,Wang M Y.Finite element inversion of the coefficients of acoustic equation in frequency domain [J].Chinese J .Geophvs.(in Chinese),2001,44(6) ,852～864．
[26] 许琨,王妙月.利用地质规则块体建模方法的频率域有限元弹性波速度反演[J].地球物理学报,2004,47(4) :708～717． Xu K,Wang M Y.Frequency-domain finite-element inversion of elastic-wave velocity using the geological regularblocky-model method[J].Chinese J .Geophys.(in Chinese),2004 ,47 (4) :708～717．
[27] Xu K,Greenhalgh S A,Wang M Y.Comparison of sourceindependent methods of elastic waveform inversion [ J ].Geophysics,2006,71(6) :R91～R100．
[28] 吴国忱,梁锴.VTI介质频率-空间域准P波正演模拟[J]. 石油地球物理勘探,2005,40(5) :535～545． Wu G C,Liang K.Quasi P-wave forward modeling in frequency-space domain in VTI media[J].OGP,2005,40 (5) : 535～545．
[29] 吴国忱,梁锴.VTI介质准P波频率空间域组合边界条件研究[J].石油物探,2005,44(4) ,301～307． Wu G C,Liang K.Combined boundary conditions of quasi-P wave within frequency-space domain in VTI media[J].GPP,2005,44(4) : 301～307．
[30] 殷文,印兴耀,吴国忱,等.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报,2006,49(2) :561～568． YinW,Yin X Y,Wu G C ,et al.The method of finite difference of high precision elastic wave equations in the frequency domain and wavefield simulation[J].Chinese J .Geophys.(in Chinese) ,2006 ,49 (2) : 561～568．
[31] 梁锴,吴国忱,印兴耀.TTI介质qP波方程频率-空间域加权平均有限差分算子[J].石油地球物理勘探,2007,42(5) :516～525． Liang K,Wu G C,Yin X Y.Weighted mean finite-difference operator of qP wave equation in frequency-space domain for TTImedia[J].OGP,2007,42(5) :516～525．
[32] 吴国忱,罗彩明,梁锴.TTI介质弹性波频率-空间域有限差分数值模拟[J].吉林大学学报(地球科学版),2007,37(5) :1023～1033． Wu G C,Luo C M,Liang K.Frequency-space domain finite difference numerical simulation of elastic wave in TTI media [J].Journal of Jilin University (Earth Science Edition),2007,37(5) : 37(5) :1023～1033．
[33] 吴国忱,梁锴.VTI介质qP波方程高精度有限差分算子[J].地球物理学进展,2007,22(3) :896～904． Wu G C,Liang K.High precision finite difference operators for qP wave equation in VTI media[J].Progress in Geophys.(in Chinese),2007,22(3) :896～904．
[34] 吴永栓,曹辉,陈国金,等.井间地震波形反演与应用[J].勘探地球物理进展,2006,29(5) :337～340． Wu Y S,Cao H,Chen G J,et al.Crosswell seismic waveform inversion and its application [ J ].Progress in Exploration Geophysics,2006,29(5) : 337～340．
[35] 张荣忠,袁诚.地震波形反演技术的应用[J].天然气工业,2009,29(2) :52～54． Zhang R Z,Yuan C.Application of the seismic waveform inversion technique[J].Natural Gas Industry,2009,29(2) :52～54．
[36] 龙桂华,李小凡,张美根,等.频率域粘弹性声波透射波形速度反演.地震学报,2009,31(1) :32～41． Long G H,Li X F,Zhang M G,et al.Visco-acoustic transmission waveform inversion for velocity structure in space-frequency domain[J].Acta Seismologica Sinica,2009,31 (1) : 32～41．
[37] Kosloff D,Baysalt E.Forward modeling by a Fourier method [J].Geophysics,1982,47(10) :1402～1412．
[38] Kosloff D,Reshef M,Loewenthal D.Elastic wave calculations by the Fourier method [J].Bull.Seis.Soc.Am.,1984,74(3) :875～891．
[39] Virieux J.SH-wave propagation in heterogeneous media[J].Geophysics,1984,49(11) :1933～1957．
[40] Virieux J.P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method [ J ].Geophysics,1986,51(4) :889～901．
[41] Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms: Numerical results [J].Geophysics,1986,51(7) :1387～1403．
[42] Mora P.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,1987,52(9) : 1211 ～1228．
[43] Mora P.Elastic wave-field inversion of reflection and transmission data[J].Geophysics,1988,53(6) : 750～759．
[44] Luo Y,Schuster G T.Wave-Equation Traveltime + WaveformInversion[C].60th Ann.Internat.Mtg.,Sot.Expl.Geophys.,Expanded abstracts,1990: 1223～1225．
[45] Zhou C X,Cai W Y,Luo Y,et al.Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data [J].Geophysics,1995,60(3) : 765～773．
[46] Zhou C X,Schuster G T,Hassanzadeh S,et al.Elastic wave-equation traveltime and waveform inversion of crosswell seismic data[J].Geophysics,1997,62(3) : 853～868．
[47] 丁继才,常旭,刘伊克,等.基于声波方程的井间地震数据快速WTW反演方法[J].地球物理学报,2007,50(5) :1527 ～1533 Ding J C ,Chang X,Liu Y K,et al.Rapid method for acoustic wave-equation WTW inversion of crosshole seismic data[J].Chinese J .Geophys.(in Chinese),2007 ,50 (5) : 1527～1533
[48] Bunks C,Saleck F M,Zaleski S,et al.Multiscale seismic waveform inversion [J].Geophysics,1995,60 (5) : 1457 ～ 1473．
[49] Yoon K,Shin C,Marfurt K J.Waveform inversion using time-windowed back propagation[C].73th Ann.Internat.Mtg.,Sot.Expl.Geophys.,Expanded abstracts,2003: 690～693．
[50] Marfurt K J.Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations [J ].Geophysics,1984,49(5) : 533～549．
[51] Davis T A.Algorithm 832: UMFPACK,an unsymmetricpattern multifrontal method [J ].ACM Transactions on Mathematical Software,2004,30(2) : 196～199．
[52] Amestoy P R,Duff I S,L'Excellent J Y.Multifrontal parallel distributed symmetric and unsymmetric solvers[J].Comput.Methods Appl.Mech.Engrg.,2000,184: 501～520．
[53] Van der Vorst H A.BI-CGSTAB: A fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems[J].SIAM J.Sci.Statist.Comput.,1992,13: 631～644．
[54] Saad Y,Schultz M.GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems[J].SIAM J.Sci.Statist.Comput.,1986,7: 856～869．
[55] Riyanti C D,Kononov A,Erlangga Y A,et al.A parallel Multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation [ J ].Journal of Computational Physics,2007,224 : 431～448．
[56] Plessix R E.A Helmholtz iterative solver for 3D seismicimaging problems[J].Geophysics,2007,72 (5) : SM185 ～ SM194．
[57] Song Z M,Willianson P R.Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅰ: 2． 5-D modeling method[J].Geophysics,1995,60(3) :784～795．
[58] Song Z M,Willianson P R,Pratt R G.Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅱ: Inversion method,synthetic experiments and real-data results[J].Geophysics,1995,60(3) :796～809．
[59] Jo C,Shin C,Suh J H.An optimal 9-point,finite-difference,frequency-space 2d scalar wave extrapolator[J].Geophysics,1996,61(2) :529～537．
[60] Shin C,Sohnz H.A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator [J].Geophysics,1998,63(1) :289～296．
[61] Stekl I,Pratt R G.Accurate visco-elastic modeling by frequency-domain finite differences using rotated operators [J].Geophysics,1998,63(5) :1779～1794．
[62] Min D J,Shin C,Kwon B D,et al.Improved frequencydomain elastic wave modeling using weighted averaging difference operators[J].Geophysics,2000,65 (3) : 884～ 895．
[63] Min D J,Yoo H S,Shin C,et al.Weighted averaging finiteelement method for scalar wave equation in the frequency domain[J].J.Seism.Explor.,2002,11:197～222．
[64] Min D J,Shin C,Pratt R G,et al.Weighted-Averaging Finite-Element Method for 2D Elastic Wave Equations in the Frequency Domain[J].Bull.Seis.Soc.Am.,2003,93(2) : 904～921．
[65] Luo Y,Schuster G.Parsimonious staggered grid finitedifferencing of the wave equation[J].Geophys.Res.Lett.,1990,17(2) :155～158．
[66] Hustedt B,Operto S,Virieux J.Mixed grid and staggered grid finite-difference methods for frequency-domain acoustic wavemodelling[J].Geophys.J.Int.,2004,157:1269～ 1296．
[67] Brossier R,Virieux J,Operto S.Parsimonious finite-volume frequency-domain method for 2-D P-SV-wave modelling[J].Geophys.J.Int.,2008,175:541～559．
[68] Operto S,Virieux J,Amestoy P,et al.3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study[J].Geophysics,2007,72 (5) :SM195～ SM211．
[69] Cao S,Greenhalgh S.2． 5-D modeling of seismic wave propagation: Boundary condition,stability criterion,and efficiency[J].Geophysics,1998,63(6) :2082～2090．
[70] Zhou B,Greenhalgh S A.Composite boundary valued solution of the 2． 5-D Green's function for arbitrary acoustic media seismograms[J].Geophysics,1998,63 (5) :1813～ 1823．
[71] Zhou B,Greenhalgh S.An adaptive wavenumber sampling strategy for 2． 5D seismic-wave modeling in the frequency domain[J].PureAppl.Geophys.,2006,163: 1399～1416．
[72] Shin C.Sponge boundary condition for frequency-domain modeling[J].Geophysics,1995,60(6) :1870～1874．
[73] Liao Q B,McMechan G A.Multifrequency viscoacoustic modeling and inversion[J].Geophysics,1996,61(5) : 1371～ 1378．
[74] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].J.Comput.Phys.,1994,114: 185～200．
[75] Min D J,Shin C,Yoo H S.Surface boundary condition in finite-difference elastic wave modeling[J].Bull.Seis.Soc.Am.,2004,94(1) :237～250．
[76] Hustedt B,Operto S,Virieux J.A multi-level directiterative solver for seismic wave propagation modeling space and wavelet approaches[J].Geophys.J.Int.,2003,157: 1269～1296．
[77] Warner M,Stekl I,Umpleby A.Full wavefield seismic tomography-Iterative forward modelling in 3D[C].69th Annual Conference and Exhibition,EAGE,ExtendedAbstracts,2007,C025．
[78] Sirgue L,Pratt R G.Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies[J].Geophysics,2004,69(1) :231～248．
[79] Wang Y H,Rao Y.Reflection seismic waveform tomography [J].J.Geophys.Res.,2009,114,B03304,doi:10． 1029/ 2008JB005916．
[80] Shin C,Min D J.Waveform inversion using a logarithmic wavefield[J].Geophysics,2006,71(3) :R31～R42．
[81] Shin C,Pyun S,Bednar J B.Comparison of waveform inversion,part 1: Conventional wavefield vs.logarithmic wavefield[J].Geophysical Prospecting,2007,55:449～464．
[82] Bednar J B,Shin C,Pyun S.Comparison of waveform inversion,part 2: phase approach [ J ].Geophysical Prospecting,2007,55 :465～475．
[83] Pyun S,Shin C,Bednar J B.Comparison of waveform inversion,part 3: amplitude approach [J].Geophysical Prospecting,2007,55:477～485．
[84] Clarbout J,Muir F.Robust modeling with erratic data [J].Geophysics,1973,38(5) :826～844．
[85] Guitton A,Symes W W.Robust inversion of seismic data using the Huber norm[J].Geophysics,2003,68(4) :1310～ 1319．
[86] Chung W,Ha T,Ha W,et al.Robust seismic waveform inversion using back-propagation algorithm[C].77th Annual International Meeting,SEG,Expanded Abstracts,2007: 1780～1784．
[87] Shin C,Ha W.2008． A comparison between the behavior of objective functions for waveform inversion in the frequency and Laplace domains [J ].Geophysics,73 (5) : VE119 ～ VE133．
[88] Zhou B,Greenhalgh S A.Crosshole seismic inversion with normalized full-waveform amplitude data [J ].Geophysics,2003,68(4) :1320～1330．
[89] Lee K H,Kim H J.Source independent full-waveform inversion of seismic data [J].Geophysics,2003,68 (6) : 2010～2015．
[90] Choi Y,Shin C,Min D J,et al.Efficient calculation of the steepest descent direction for source-independent seismic waveform inversion: An amplitude approach[J].J.Comput.Phys.,2005,208:455～468．
[91] Greenhalgh S A,Zhou B.Surface seismic imaging by multifrequency amplitude inversion[J].Exploration Geophysics,2003,34(4) : 217～224．
[92] Shin C,Yoon K G,Marfurt K J,et al.Efficient calculation of a partial derivative wavefield using reciprocity for seismic imaging and inversion[J].Geophysics,2001,66(6) :1856～ 1863．
[93] Shin C,Jang S and Min D J.Improved amplitude preservation for prestack depth migration by inverse scattering theory [J].Geophysical prospecting,2001,49: 592～606．
[94] Jang U,Min D J,Shin C.Comparison of scaling methods for waveform inversion[J].Geophysical Prospecting,2009,57:49～59．
[95] Choi Y,Min D J,Shin C.Frequency-domain elastic full waveform inversion using the new pseudo-Hessian matrix: Experience of elastic Marmousi-2 synthetic data [J].Bull.Seis.Soc.Am.,2008,98(5) :2402～2415．
[96] Vigh D,Starr E W.Comparisons for Waveform Inversion,Time Domain or Frequency Domain? [C].78th Annual International Meeting,SEG,Expanded Abstracts,2008:1890～1894．
[97] Shi Y M,Zhao W Z,Cao H.Nonlinear process control of wave-equation inversion and its application in the detection of gas[J].Geophysics,2007,72(1) : 9～18．
[2] 苑书金.叠前地震反演技术的进展及其在岩性油气藏勘探中的应用[J].地球物理学进展,2007,22(3) :879～885． Yuan S J.Progress of pre-stack inversion and application in exploration of the lithological reservoirs [ J ].Progress in Geophys.(in Chinese),2007,22(3) : 879～885．
[3] 杨文采.非线性波动方程地震反演的方法原理及问题[J].地球物理学进展,1992,7(1) :9～19． Yang W C.Methods and problems in seismic inversion of nonlinear wave equations [ J ].Progress in Geophys.(in Chinese),1992,7(1) :9～19．
[4] 杨文采.非线性地球物理反演方法:回顾与展望[J].地球物理学进展,2002,17(2) :255～261． Yang W C.Non-linear geophysical inversion methods: review and perspective[J].Progress in Geophys.(in Chinese),2002,17(2) : 255～261．
[5] Williamson P R.A guide to the limits of resolution imposed by scattering in ray tomography[J].Geophysics,1991,56(2) : 202～207．
[6] Williamson P R,Worthington M H.Resolution limits in ray tomography due to wave behavior: Numerical experiments[J].Geophysics,1993,58(5) :723～735．
[7] Wu R S,Toksoz M N.Diffraction tomography and multisource holography applied to seismic imaging.Geophysics,1987,52(1) :11～25．
[8] Pratt R G.Seismic waveform inversion in the frequency domain,Part 1: Theory and verification in a physical scale model[J].Geophysics,1999,64(3) :888～901．
[9] Pratt R G,Shipp R M.Seismic waveform inversion in the frequency domain,Part 2: Fault delineation in sediments using crosshole data[J].Geophysics,1999,64(3) :902～914．
[10] 曹辉,陈国金,郭建,等.频率域2． 5-D井间波形层析成像及其实际应用[J].地球物理学报,2008,51(3) :91 5～923． Cao H,Chen G J,Guo J,et al.The 2． 5-D cross-well waveform tomography in the frequency domain and its practical application[J].Chinese J.Geophys.(in Chinese),2008,51(3) :915～923．
[11] Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49 (8) : 1259 ～ 1266．
[12] Tarantola A.A strategy for nonlinear elastic inversion of seismic reflection data[J].Geophysics,1986,51 (10) : 1893 ～ 1903．
[13] Pratt R G.Frequency-domain elastic wave modeling by finite differences: A tool for crosshole seismic imaging [J].Geophysics,1990,55(5) :626～632．
[14] Pratt R G,Worthington M H.Inverse theory applied to multi-source cross-hole tomography,Part 1 : Acoustic waveequation method[J].Geophysical Prospecting,1990,38:287～ 310．
[15] Pratt R G.Inverse theory applied to multi source cross-hole tomography,Part 2: Elastic wave-equation method [ J ].Geophysical Prospecting,1990,38:311～329．
[16] Pratt R G,Shin C,Hicks G J.Gauss-Newton and full Newton methods in frequency domain seismic waveform inversion[J].Geophys.J.Int.,1998,133:341～362．
[17] Shin C,Cha Y H.Waveform inversion in the Laplace domain [J].Geophys.J.Int.,2008,173:922～931．
[18] Operto S,Virieux J,Dessa J X,et al.Crustal seismic imaging from multifold ocean bottom seismometer data by frequency domain full waveform tomography: Application to the eastern Nankai trough[J].J.Geophys.Res.,2006,111,B09306,doi: 10． 1029/2005JB003835．
[19] Bleibinhaus F,Hole J A,Ryberg T,et al.Structure of the California Coast Ranges and San Andreas Fault at SAFOD from seismic waveform inversion and reflection imaging[J].J.Geophys.Res.,2007,112,B06315,doi: 10． 1029/2006JB004611．
[20] Malinowski M,Operto S.Quantitative imaging of the Permo-Mesozoic complex and its basement by frequency domain waveform tomography of wide-aperture seismic data from the Polish Basin[J].Geophysical Prospecting,2008,56: 805～825．
[21] Gao F C,Levander A,Pratt R G,et al.Waveform tomography at a groundwater contamination site: VSP-surface data set[J].Geophysics,2007,71(1) :H1～H11．
[22] Gao F C,Levander A,Pratt R G,et al.Waveform tomography at a groundwater contamination site: Surface reflection data[J].Geophysics,2007,72(5) :G45～G55．
[23] Ben-Hadj-Ali H,Operto S,Virieux J.Velocity model building by 3D frequency-domain,full-waveform inversion of wide-aperture seismic data [J].Geophysics,2008,73 (5) :VE101～VE117．
[24] Hicks G J,Pratt R G.Reflection waveform inversion using local descent methods Estimating attenuation and velocity over a gas-sand deposit[J].Geophysics,2001,66(2) : 598～612．
[25] 许琨,王妙月.声波方程频率域有限元参数反演[J].地球物理学报,2001,44(6) ,852～864． Xu K,Wang M Y.Finite element inversion of the coefficients of acoustic equation in frequency domain [J].Chinese J .Geophvs.(in Chinese),2001,44(6) ,852～864．
[26] 许琨,王妙月.利用地质规则块体建模方法的频率域有限元弹性波速度反演[J].地球物理学报,2004,47(4) :708～717． Xu K,Wang M Y.Frequency-domain finite-element inversion of elastic-wave velocity using the geological regularblocky-model method[J].Chinese J .Geophys.(in Chinese),2004 ,47 (4) :708～717．
[27] Xu K,Greenhalgh S A,Wang M Y.Comparison of sourceindependent methods of elastic waveform inversion [ J ].Geophysics,2006,71(6) :R91～R100．
[28] 吴国忱,梁锴.VTI介质频率-空间域准P波正演模拟[J]. 石油地球物理勘探,2005,40(5) :535～545． Wu G C,Liang K.Quasi P-wave forward modeling in frequency-space domain in VTI media[J].OGP,2005,40 (5) : 535～545．
[29] 吴国忱,梁锴.VTI介质准P波频率空间域组合边界条件研究[J].石油物探,2005,44(4) ,301～307． Wu G C,Liang K.Combined boundary conditions of quasi-P wave within frequency-space domain in VTI media[J].GPP,2005,44(4) : 301～307．
[30] 殷文,印兴耀,吴国忱,等.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报,2006,49(2) :561～568． YinW,Yin X Y,Wu G C ,et al.The method of finite difference of high precision elastic wave equations in the frequency domain and wavefield simulation[J].Chinese J .Geophys.(in Chinese) ,2006 ,49 (2) : 561～568．
[31] 梁锴,吴国忱,印兴耀.TTI介质qP波方程频率-空间域加权平均有限差分算子[J].石油地球物理勘探,2007,42(5) :516～525． Liang K,Wu G C,Yin X Y.Weighted mean finite-difference operator of qP wave equation in frequency-space domain for TTImedia[J].OGP,2007,42(5) :516～525．
[32] 吴国忱,罗彩明,梁锴.TTI介质弹性波频率-空间域有限差分数值模拟[J].吉林大学学报(地球科学版),2007,37(5) :1023～1033． Wu G C,Luo C M,Liang K.Frequency-space domain finite difference numerical simulation of elastic wave in TTI media [J].Journal of Jilin University (Earth Science Edition),2007,37(5) : 37(5) :1023～1033．
[33] 吴国忱,梁锴.VTI介质qP波方程高精度有限差分算子[J].地球物理学进展,2007,22(3) :896～904． Wu G C,Liang K.High precision finite difference operators for qP wave equation in VTI media[J].Progress in Geophys.(in Chinese),2007,22(3) :896～904．
[34] 吴永栓,曹辉,陈国金,等.井间地震波形反演与应用[J].勘探地球物理进展,2006,29(5) :337～340． Wu Y S,Cao H,Chen G J,et al.Crosswell seismic waveform inversion and its application [ J ].Progress in Exploration Geophysics,2006,29(5) : 337～340．
[35] 张荣忠,袁诚.地震波形反演技术的应用[J].天然气工业,2009,29(2) :52～54． Zhang R Z,Yuan C.Application of the seismic waveform inversion technique[J].Natural Gas Industry,2009,29(2) :52～54．
[36] 龙桂华,李小凡,张美根,等.频率域粘弹性声波透射波形速度反演.地震学报,2009,31(1) :32～41． Long G H,Li X F,Zhang M G,et al.Visco-acoustic transmission waveform inversion for velocity structure in space-frequency domain[J].Acta Seismologica Sinica,2009,31 (1) : 32～41．
[37] Kosloff D,Baysalt E.Forward modeling by a Fourier method [J].Geophysics,1982,47(10) :1402～1412．
[38] Kosloff D,Reshef M,Loewenthal D.Elastic wave calculations by the Fourier method [J].Bull.Seis.Soc.Am.,1984,74(3) :875～891．
[39] Virieux J.SH-wave propagation in heterogeneous media[J].Geophysics,1984,49(11) :1933～1957．
[40] Virieux J.P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method [ J ].Geophysics,1986,51(4) :889～901．
[41] Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms: Numerical results [J].Geophysics,1986,51(7) :1387～1403．
[42] Mora P.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,1987,52(9) : 1211 ～1228．
[43] Mora P.Elastic wave-field inversion of reflection and transmission data[J].Geophysics,1988,53(6) : 750～759．
[44] Luo Y,Schuster G T.Wave-Equation Traveltime + WaveformInversion[C].60th Ann.Internat.Mtg.,Sot.Expl.Geophys.,Expanded abstracts,1990: 1223～1225．
[45] Zhou C X,Cai W Y,Luo Y,et al.Acoustic wave-equation traveltime and waveform inversion of crosshole seismic data [J].Geophysics,1995,60(3) : 765～773．
[46] Zhou C X,Schuster G T,Hassanzadeh S,et al.Elastic wave-equation traveltime and waveform inversion of crosswell seismic data[J].Geophysics,1997,62(3) : 853～868．
[47] 丁继才,常旭,刘伊克,等.基于声波方程的井间地震数据快速WTW反演方法[J].地球物理学报,2007,50(5) :1527 ～1533 Ding J C ,Chang X,Liu Y K,et al.Rapid method for acoustic wave-equation WTW inversion of crosshole seismic data[J].Chinese J .Geophys.(in Chinese),2007 ,50 (5) : 1527～1533
[48] Bunks C,Saleck F M,Zaleski S,et al.Multiscale seismic waveform inversion [J].Geophysics,1995,60 (5) : 1457 ～ 1473．
[49] Yoon K,Shin C,Marfurt K J.Waveform inversion using time-windowed back propagation[C].73th Ann.Internat.Mtg.,Sot.Expl.Geophys.,Expanded abstracts,2003: 690～693．
[50] Marfurt K J.Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations [J ].Geophysics,1984,49(5) : 533～549．
[51] Davis T A.Algorithm 832: UMFPACK,an unsymmetricpattern multifrontal method [J ].ACM Transactions on Mathematical Software,2004,30(2) : 196～199．
[52] Amestoy P R,Duff I S,L'Excellent J Y.Multifrontal parallel distributed symmetric and unsymmetric solvers[J].Comput.Methods Appl.Mech.Engrg.,2000,184: 501～520．
[53] Van der Vorst H A.BI-CGSTAB: A fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems[J].SIAM J.Sci.Statist.Comput.,1992,13: 631～644．
[54] Saad Y,Schultz M.GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems[J].SIAM J.Sci.Statist.Comput.,1986,7: 856～869．
[55] Riyanti C D,Kononov A,Erlangga Y A,et al.A parallel Multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation [ J ].Journal of Computational Physics,2007,224 : 431～448．
[56] Plessix R E.A Helmholtz iterative solver for 3D seismicimaging problems[J].Geophysics,2007,72 (5) : SM185 ～ SM194．
[57] Song Z M,Willianson P R.Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅰ: 2． 5-D modeling method[J].Geophysics,1995,60(3) :784～795．
[58] Song Z M,Willianson P R,Pratt R G.Frequency-domain acoustic-wave modeling and inversion of crosshole data: Part Ⅱ: Inversion method,synthetic experiments and real-data results[J].Geophysics,1995,60(3) :796～809．
[59] Jo C,Shin C,Suh J H.An optimal 9-point,finite-difference,frequency-space 2d scalar wave extrapolator[J].Geophysics,1996,61(2) :529～537．
[60] Shin C,Sohnz H.A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator [J].Geophysics,1998,63(1) :289～296．
[61] Stekl I,Pratt R G.Accurate visco-elastic modeling by frequency-domain finite differences using rotated operators [J].Geophysics,1998,63(5) :1779～1794．
[62] Min D J,Shin C,Kwon B D,et al.Improved frequencydomain elastic wave modeling using weighted averaging difference operators[J].Geophysics,2000,65 (3) : 884～ 895．
[63] Min D J,Yoo H S,Shin C,et al.Weighted averaging finiteelement method for scalar wave equation in the frequency domain[J].J.Seism.Explor.,2002,11:197～222．
[64] Min D J,Shin C,Pratt R G,et al.Weighted-Averaging Finite-Element Method for 2D Elastic Wave Equations in the Frequency Domain[J].Bull.Seis.Soc.Am.,2003,93(2) : 904～921．
[65] Luo Y,Schuster G.Parsimonious staggered grid finitedifferencing of the wave equation[J].Geophys.Res.Lett.,1990,17(2) :155～158．
[66] Hustedt B,Operto S,Virieux J.Mixed grid and staggered grid finite-difference methods for frequency-domain acoustic wavemodelling[J].Geophys.J.Int.,2004,157:1269～ 1296．
[67] Brossier R,Virieux J,Operto S.Parsimonious finite-volume frequency-domain method for 2-D P-SV-wave modelling[J].Geophys.J.Int.,2008,175:541～559．
[68] Operto S,Virieux J,Amestoy P,et al.3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study[J].Geophysics,2007,72 (5) :SM195～ SM211．
[69] Cao S,Greenhalgh S.2． 5-D modeling of seismic wave propagation: Boundary condition,stability criterion,and efficiency[J].Geophysics,1998,63(6) :2082～2090．
[70] Zhou B,Greenhalgh S A.Composite boundary valued solution of the 2． 5-D Green's function for arbitrary acoustic media seismograms[J].Geophysics,1998,63 (5) :1813～ 1823．
[71] Zhou B,Greenhalgh S.An adaptive wavenumber sampling strategy for 2． 5D seismic-wave modeling in the frequency domain[J].PureAppl.Geophys.,2006,163: 1399～1416．
[72] Shin C.Sponge boundary condition for frequency-domain modeling[J].Geophysics,1995,60(6) :1870～1874．
[73] Liao Q B,McMechan G A.Multifrequency viscoacoustic modeling and inversion[J].Geophysics,1996,61(5) : 1371～ 1378．
[74] Berenger J P.A perfectly matched layer for the absorption of electromagnetic waves[J].J.Comput.Phys.,1994,114: 185～200．
[75] Min D J,Shin C,Yoo H S.Surface boundary condition in finite-difference elastic wave modeling[J].Bull.Seis.Soc.Am.,2004,94(1) :237～250．
[76] Hustedt B,Operto S,Virieux J.A multi-level directiterative solver for seismic wave propagation modeling space and wavelet approaches[J].Geophys.J.Int.,2003,157: 1269～1296．
[77] Warner M,Stekl I,Umpleby A.Full wavefield seismic tomography-Iterative forward modelling in 3D[C].69th Annual Conference and Exhibition,EAGE,ExtendedAbstracts,2007,C025．
[78] Sirgue L,Pratt R G.Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies[J].Geophysics,2004,69(1) :231～248．
[79] Wang Y H,Rao Y.Reflection seismic waveform tomography [J].J.Geophys.Res.,2009,114,B03304,doi:10． 1029/ 2008JB005916．
[80] Shin C,Min D J.Waveform inversion using a logarithmic wavefield[J].Geophysics,2006,71(3) :R31～R42．
[81] Shin C,Pyun S,Bednar J B.Comparison of waveform inversion,part 1: Conventional wavefield vs.logarithmic wavefield[J].Geophysical Prospecting,2007,55:449～464．
[82] Bednar J B,Shin C,Pyun S.Comparison of waveform inversion,part 2: phase approach [ J ].Geophysical Prospecting,2007,55 :465～475．
[83] Pyun S,Shin C,Bednar J B.Comparison of waveform inversion,part 3: amplitude approach [J].Geophysical Prospecting,2007,55:477～485．
[84] Clarbout J,Muir F.Robust modeling with erratic data [J].Geophysics,1973,38(5) :826～844．
[85] Guitton A,Symes W W.Robust inversion of seismic data using the Huber norm[J].Geophysics,2003,68(4) :1310～ 1319．
[86] Chung W,Ha T,Ha W,et al.Robust seismic waveform inversion using back-propagation algorithm[C].77th Annual International Meeting,SEG,Expanded Abstracts,2007: 1780～1784．
[87] Shin C,Ha W.2008． A comparison between the behavior of objective functions for waveform inversion in the frequency and Laplace domains [J ].Geophysics,73 (5) : VE119 ～ VE133．
[88] Zhou B,Greenhalgh S A.Crosshole seismic inversion with normalized full-waveform amplitude data [J ].Geophysics,2003,68(4) :1320～1330．
[89] Lee K H,Kim H J.Source independent full-waveform inversion of seismic data [J].Geophysics,2003,68 (6) : 2010～2015．
[90] Choi Y,Shin C,Min D J,et al.Efficient calculation of the steepest descent direction for source-independent seismic waveform inversion: An amplitude approach[J].J.Comput.Phys.,2005,208:455～468．
[91] Greenhalgh S A,Zhou B.Surface seismic imaging by multifrequency amplitude inversion[J].Exploration Geophysics,2003,34(4) : 217～224．
[92] Shin C,Yoon K G,Marfurt K J,et al.Efficient calculation of a partial derivative wavefield using reciprocity for seismic imaging and inversion[J].Geophysics,2001,66(6) :1856～ 1863．
[93] Shin C,Jang S and Min D J.Improved amplitude preservation for prestack depth migration by inverse scattering theory [J].Geophysical prospecting,2001,49: 592～606．
[94] Jang U,Min D J,Shin C.Comparison of scaling methods for waveform inversion[J].Geophysical Prospecting,2009,57:49～59．
[95] Choi Y,Min D J,Shin C.Frequency-domain elastic full waveform inversion using the new pseudo-Hessian matrix: Experience of elastic Marmousi-2 synthetic data [J].Bull.Seis.Soc.Am.,2008,98(5) :2402～2415．
[96] Vigh D,Starr E W.Comparisons for Waveform Inversion,Time Domain or Frequency Domain? [C].78th Annual International Meeting,SEG,Expanded Abstracts,2008:1890～1894．
[97] Shi Y M,Zhao W Z,Cao H.Nonlinear process control of wave-equation inversion and its application in the detection of gas[J].Geophysics,2007,72(1) : 9～18．
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心