初至波波形反演方法及其数值模拟试验
|
摘要
近地表速度模型的反演精度是影响复杂地表条件下地下地质构造成像效果的主要因素之一。初至波波形反演方法将速度建模问题转化为寻求地震数据信息和理论模型的最佳拟合,即转化为求解目标函数极小化问题,利用观测记录和模型正演记录之间的最佳匹配关系建立地下模型。梯度类波形反演方法的关键因素是目标函数、梯度和步长,梯度问题可以通过确定合适的目标函数来解决,选取合适的步长可以有效地提高波形反演的计算效率和反演精度。初至波波形反演方法类似于对地表模型同时做偏移和层析,能够同时分辨出高波数和低波数成分,并且同时恢复速度值和速度界面。初至波波形反演结果与初至波射线走时层析反演结果的对比分析表明,对于速度变化比较缓慢的起伏地表,前者的反演效果明显优于后者。
The inversion accuracy of near surface velocity model is one of the main factors for influencing the imaging of subsurface geologic structure with complex surface conditions.First arrival waveform inversion transforms velocity modeling problem into finding best fitting between seismic data information and theoretical model, which is to solve minimization problem of objective function,and use the best matching relationship between observation record and model forwarding record to establish subsurface model.The key factors for gradient-type waveform inversion method are objective function,gradient and step length.Gradient problem can be solved by defining proper objective function.Selecting proper step length can effectively improve the computation efficiency and inversion accuracy. First arrival waveform inversion method is similar to simultaneous migration and tomography on surface model,which can simultaneously recognize high wave number and low wave number component.and recover velocity value and interface.The first arrival waveform inversion results are compared to first arrival ray travel time tomography inversion results.The comparison results show that the former inversion is obviously superior to the later one at the rough surface with slow velocity change.
The inversion accuracy of near surface velocity model is one of the main factors for influencing the imaging of subsurface geologic structure with complex surface conditions.First arrival waveform inversion transforms velocity modeling problem into finding best fitting between seismic data information and theoretical model, which is to solve minimization problem of objective function,and use the best matching relationship between observation record and model forwarding record to establish subsurface model.The key factors for gradient-type waveform inversion method are objective function,gradient and step length.Gradient problem can be solved by defining proper objective function.Selecting proper step length can effectively improve the computation efficiency and inversion accuracy. First arrival waveform inversion method is similar to simultaneous migration and tomography on surface model,which can simultaneously recognize high wave number and low wave number component.and recover velocity value and interface.The first arrival waveform inversion results are compared to first arrival ray travel time tomography inversion results.The comparison results show that the former inversion is obviously superior to the later one at the rough surface with slow velocity change.
引文
1 Pica A,Diet J P,Tarantola A.Nonlinear inversion of seismic reflection data in a laterally medium[J].Geophysics, 1990,55(3):284-292
2 Tarantola A.Inversion of seismic data in acoustic approximation [J].Geophysics,1984,49(8):1259- ??1266
3Zhou C X,Gerard T.Waveform inversion of subwell velocity structure[J].Expanded Abstracts of 63~(rd) Annual Internat SEG Mtg,1993,106-109
4 Woodward M J.Wave-equation tomography[J].Geophysics, 1992,57(1):15-26
5 Devanvey A J.Geophysical diffraction tomography[J]. IEEE Transactions on Geoscience and Remote Sensing, 1984,GE-22(1):3-13
6Wu R S,Toksoz M N.Diffraction tomography and multisource holography applied to seismic imaging[J]. Geophysics,1987,52(1):11-25
7Pratt R G,Goulty N R.Combining wave-equation imaging with traveltime tomography to form high-resolution images from crosshole data[J].Geophysics,1991, 56(2):208-224
8 Schuster G T,Aksel Q B.Wavepath eikonal traveltime inversion:Theory[J].Geophysics,1993,58(9): 1314-1323
9 Pratt R G,Shin C S,Hicks G J.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International,1998, 133(2):341-362
10 Mora P.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,1987,52(9): 1211-1228
11吴永栓,曹辉,陈国金,等.井间地震波形反演与应用[J].勘探地球物理进展,2006,29(5):337-340
12Luo Y,Schuster G T.Wave-equation traveltime inversion [J].Geophysics,1991,56(5):645-653
13Sheng J M,Leeds A,Buddensiek M,et al.Early arrival waveform tomography on near-surface refraction data [J].Geophysics,2006,71(4):U47-U57
14Sheng J M.High resolution seismic tomography with the generalized radon transform and early arrival waveform inversion[D].Salt Lake City:University of Uath, 2004.20-54
15 Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms:Numerical results[J].Geophysics,1986,51(7):1387-1403
16 Tarantola A.Inverion problem theroy and methods for model parameter estimation.philadelphia[M].Philadelphia: Society for Industrial and Applied Mathematics, 2005.203-214
17 Mora P.Inversion=migration+tomography[J].Geophysics, 1989,54(12):1575-1586
18 Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8): 1259-1266
2 Tarantola A.Inversion of seismic data in acoustic approximation [J].Geophysics,1984,49(8):1259- ??1266
3Zhou C X,Gerard T.Waveform inversion of subwell velocity structure[J].Expanded Abstracts of 63~(rd) Annual Internat SEG Mtg,1993,106-109
4 Woodward M J.Wave-equation tomography[J].Geophysics, 1992,57(1):15-26
5 Devanvey A J.Geophysical diffraction tomography[J]. IEEE Transactions on Geoscience and Remote Sensing, 1984,GE-22(1):3-13
6Wu R S,Toksoz M N.Diffraction tomography and multisource holography applied to seismic imaging[J]. Geophysics,1987,52(1):11-25
7Pratt R G,Goulty N R.Combining wave-equation imaging with traveltime tomography to form high-resolution images from crosshole data[J].Geophysics,1991, 56(2):208-224
8 Schuster G T,Aksel Q B.Wavepath eikonal traveltime inversion:Theory[J].Geophysics,1993,58(9): 1314-1323
9 Pratt R G,Shin C S,Hicks G J.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International,1998, 133(2):341-362
10 Mora P.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,1987,52(9): 1211-1228
11吴永栓,曹辉,陈国金,等.井间地震波形反演与应用[J].勘探地球物理进展,2006,29(5):337-340
12Luo Y,Schuster G T.Wave-equation traveltime inversion [J].Geophysics,1991,56(5):645-653
13Sheng J M,Leeds A,Buddensiek M,et al.Early arrival waveform tomography on near-surface refraction data [J].Geophysics,2006,71(4):U47-U57
14Sheng J M.High resolution seismic tomography with the generalized radon transform and early arrival waveform inversion[D].Salt Lake City:University of Uath, 2004.20-54
15 Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms:Numerical results[J].Geophysics,1986,51(7):1387-1403
16 Tarantola A.Inverion problem theroy and methods for model parameter estimation.philadelphia[M].Philadelphia: Society for Industrial and Applied Mathematics, 2005.203-214
17 Mora P.Inversion=migration+tomography[J].Geophysics, 1989,54(12):1575-1586
18 Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8): 1259-1266
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心