地震勘探全波形反演的应用与发展分析
收稿日期: 2013-08-22
修回日期: 2014-03-03
网络出版日期: 2014-05-10
基金资助
国家自然科学基金项目（41374061、41204048、41240027）；博士后基金项目（2012M510533）。
Application and Development of FullWaveform Inversion Research in the Seismic Exploration
Received date: 2013-08-22
Revised date: 2014-03-03
Online published: 2014-05-10
本文首先对20世纪80年代发展起来的全波形反演应用及其在勘探地球物理领域的发展进行了分析；其次，面对定量化、精细化的地震勘探要求，提出了将地震勘探全波形反演与其他数据处理环节或处理技术相结合的研究设想，并展望了全波形反演的发展趋势；最后，论述了全波形反演研究中地震波场数值模拟、反演初始速度模型获取、目标函数形式选择、寻优算法启用及各向异性介质中的应用等关键问题，并总结了通过Laplace域的全波形反演获取反演初始速度模型、结合射线追踪并充分发挥并行计算之于波动方程方法来模拟地震波场的巨大优势，及灵活选用反演目标函数形式和寻优算法更新速度模型参数来加快全波形反演方法的实用化进程。
黄金, 高星, 王伟 . 地震勘探全波形反演的应用与发展分析[J]. 地球信息科学学报, 2014 , 16(3) : 396 -401 . DOI: 10.3724/SP.J.1047.2014.00396
Firstly, this paper gives a simple review of the research of full waveform inversion (FWI) which emerged in the 80's of the 20th century. Aimed at the increasingly quantitative and precision requirement of the seismic exploration, a viewpoint is proposed in this paper that FWI should combine with other data processing steps or techniques in the future. In terms of the data processing of seismic exploration or from the standpoint of the final geological interpretation, a combination with other data process steps or techniques would be the development trend of FWI. Finally, author holds a discussion about the development of FWI, especially focusing on the key problems of FWI research, including the seismic wavefield modeling, inversion initial velocity model picking, object function style selecting, optimization algorithm choosing and its application in the anisotropic medium. By using Laplace FWI to obtain inversion-based initial velocity model, making full use of the advantages of combining the wave equation method which benefited from parallel computing and ray tracing to simulate the seismic wavefield, and selecting the form of inversion objective function and optimization algorithm flexibly to update velocity model parameters, it could be concluded that the practical use of FWI will be accelerated greatly, especially its application in the anisotropic medium.
[1] Tarantola A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics, 1984,49(8):1259-1266.
[2] Tarantola A. A strategy for nonlinear elastic inverse of seismic reflection data[J]. Geophysics,1986,51(10):1893-1903.
[3] Pica A, Diet J P, Tarantola A. Nonlinear inversion of seismic reflection data in laterally invariant medium[J]. Geophysics, 1990,55(3):284-292.
[4] Crase E, Wideman C, Noble M, et al. Nonlinear elastic waveform inversion of land seismic reflection data[J]. Journal of Geophysical Research, 1992,97(B4):4685-4703.
[5] Pratt R G. Frequency domain elastic wave modeling by finite difference: A tool for crosshole seismic imaging[J]. Geophysics, 1990,55(5):626-632.
[6] Pratt R G, Shin C. Gauss-Newton and full Newton methods in frequency space seismic waveform inversion[J]. Geophysical Journal International, 1998,133(2):341-362.
[7] Pratt R G, Shipp R M. Seismic waveform inversion in the frequency domain, part 2: fault delineation in sediments using cross-hole data[J]. Geophysics, 1999,64(3):902-914.
[8] Shin C, Cha Y H. Waveform inversion in the Laplace domain[J]. Geophysical Journal International, 2008,173(3): 922-931.
[9] Pecher I A, Minshull T A, Singh S C, et al. Velocity structure of a bottom simulating reflector offshore Peru: Results from full waveform inversion[J]. Earth and Planetary Science Letters, 1996,139(3):459-469.
[10] 宋海斌, Osamu M, Shin' ichi K.天然气水合物似海底反射层的全波形反演[J].地球物理学报,2003,46(1):42-46.
[11] 霍元媛, 张明. 南海北部似海底反射层速度结构全波形反演[J]. 物探化探计算技术, 2010,32(2):182-187.
[12] Ma Y, Hale D, Gong B, et al. Image-guided sparse-model full waveform inversion[J]. Geophysics,2012, 77(4):189-198.
[13] 师学明,王家映.遗传算法[J].工程地球物理学报,2008,5 (2):129-140.
[14] 张世俊, 杨慧珠, 董渊, 等. 遗传算法反演HTI 介质各向异性参数[J]. 石油地球物理勘探, 2002,37(1):24-28.
[15] 师学明, 王家映. 模拟退火法[J]. 工程地球物理学报, 2007,4(3):165-174.
[16] 任义庆, 徐仲达, 马在田. 应用模拟退火法反演横波速度[J]. 石油地球物理勘探, 1996,31(5):677-684.
[17] 刘巧霞, 田仁飞, 钟宙灿. 模拟退火地震参数反演研究[J]. 长江大学学报(自然科学版), 2009,6(2):186-188.
[18] Wang B L, Gao J. Fast full waveform inversion of multishot seismic data[C]. Expanded Abstract of 80th Annual International SEG Meeting, 2010,1055-1058.
[19] Berkhout A J. Combining full wavefield migration and full waveform inversion[C]. Expanded Abstract of 81th Annual International SEG Meeting, 2011,3321-3325.
[20] Berkhout A J. Combining full wavefield migration and full waveform inversion, a glance into the future of seismic imaging[J]. Geophysics, 2012,77(2):43-50.
[21] 李信富, 李小凡, 张美根. 地震波数值模拟方法研究综述[J]. 防灾减灾工程学报, 2007,27(2):241-248.
[22] 朱多林, 白超英. 基于波动方程理论的地震波场数值模拟方法综述[J]. 地球物理学进展, 2012,26(5):1588-1599.
[23] 张文生, 何樵登. 二维横向各向同性介质的伪谱法正演模拟[J]. 石油地球物理勘探, 1998,33(3):310-319.
[24] 汪文帅, 李小凡, 鲁明文, 等. 基于多辛结构谱元法的保结构地震波场模拟[J]. 地球物理学报, 2012,55(10): 3427-3439.
[25] 张永刚. 地震波场数值模拟方法[J].石油物探, 2003,42 (2):143-148.
[26] 杨顶辉, 滕吉文. 各向异性介质中三分量地震记录的 FCT 有限差分模拟[J]. 石油地球物理勘探, 1997,32(2): 181-190.
[27] 邓飞, 王美平, 周杲, 等. 高斯射线束法地震记录合成系统的研究与开发[J]. 大庆石油地质与开发. 2006,25(6): 93-97.
[28] 邓飞, 刘超颖, 赵波, 等. 高斯射线束正演与偏移[J]. 石油地球物理勘探, 2009,44(3):265-269,281.
[29] Beydoun W B, Keho T H. The paraxial ray method[J]. Geophysics, 1987,52(12):1639-1653.
[30] 陈国金, 高志凌, 吴永栓. 井间地震层析成像中自动生成初始速度模型的方法研究[J]. 石油物探,2005,44(4): 339-342.
[31] Chung W, Shin C, Pyun S. 2-D elastic waveform inversion in the Laplace domain[C]. Expanded Abstract of 80th Annual International SEG Meeting, 2010,1059-1064.
[32] Brossier R, Operto S, Virieux J. Which data residual norm for robust elastic frequency-domain full waveform inversion[J]. Geophysics, 2010,75(3):37-46.
[33] Bube K P, Langan R T. Hybrid L1/L2 minimization fitting with applications to tomography[J]. Geophysics, 1997,62(4):1183-1195.
[34] Bube K P, Nemeth T. Fast line searches for the robust solution of linear systems in the hybrid L1/L2 and Huber norms[J]. Geophysics, 2007,72(2):A13-A17.
[35] 卞爱飞, 於文辉, 周华伟. 频率域全波形反演方法研究进展[J]. 地球物理学进展, 2010,25(3):982-993.
[36] 程金星, 董敏煜, 秦顺亭, 等. 三种算法联合迭代反演求取最佳剩余静校正量[J]. 地球物理学报, 1996,39(3): 416-423.
[37] Ma Y, Hale D. A projected Hessian matrix for full waveform inversion[C]. Expanded Abstract of 81th Annual International SEG Meeting,2011,2401-2404.
[38] 朱培民, 王家映. 共轭梯度法[J].工程地球物理学报, 2008,5(4):381-386.
[39] Wang B L, Gao J, Zhang H, et al. CUDA-based acceleration of full waveform inversion on GPU[C]. Expanded Abstract of 81th Annual International SEG Meeting, 2011, 2528-2531.
[40] Podgornova O, Charara M. Multiscale time-domain fullwaveform inversion for anisotropic elastic media[C]. Expanded Abstract of 81th Annual International SEG Meeting, 2011,2459-2464.
[41] Samec P, Blangy J P. Viscoelastic attenuation, anisotropy, and AVO[J]. Geophysics, 1992,57(3):441-450.
[42] Carcione J M, Helle H B, Zhao T. Effects of attenuation and anisotropy on reflection amplitude versus offset[J]. Geophysics, 1998,63(5):1652-1658.
[43] Alkhalifah T, Tsvankin I. Velocity analysis for transversely isotropic media[J]. Geophysics, 1995,60(5):1550-1566.
[44] Alkhalifah T. An acoustic wave equation for anisotropic media[J]. Geophysics, 2000,65(4):1239-1250.
[45] Gholami Y, Ribodetti A, Brossier R, et al. Sensitivity analysis of full waveform inversion in VTI media[C]. Extended Abstract of the 72nd EAGE Conference & Exhibition, 2010.
[46] Plessix R E, Rynja H. VTI full waveform inversion: a parameterization study with a narrow azimuth streamer data example[C]. Expanded Abstract of 80th Annual International SEG Meeting, 2010,962-966.
[47] Lee H Y, Koo J M, Min D J, et al. Frequency-domain elastic full waveform inversion for VTI media[J]. Geophysical Journal International, 2010,183(2):884-904.
[48] Jeong W, Min D, Lee G, et al. 2D frequency domain elastic full waveform inversion using finite-element method for VTI media[C]. SEG Technical Program Expanded Abstracts 2011,2654-2658.
[49] Kim Y, Cha H Y, Shin C S, et al. Improved logarithmic waveform inversion considering the power-spectrum of the wavefield[C]. Expanded Abstract of 78th Annual International SEG Meeting, 2008,1966-1969.
[50] Resnick J R. Stratigraphic filtering[J]. Pure and Applied Geophysics, 1990,132(1-2):49-65.
[51] Yang Y, Li Y, Liu T. 1D viscoelastic waveform inversion for Q structures from the surface seismic and zero-offset VSP data[J]. Geophysics, 2009,74(6):141-148.
[52] 杨午阳, 王西文, 雍学善, 等. 地震全波形反演方法研究综述[J]. 地球物理学进展,2013,28(2):766-776.
/
| 〈 |
|
〉 |