Full-waveform inversion of elastic waves in variable grids based on variable difference coefficients
ZHANG Hongjing1,2, HE Huili3, SUN Wenbo4,5, SUN Pengyuan3, LI Honghui3, ZHOU Hui6,7,8
1. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, Jiangxi 330013, China; 2. School of Geophysics and Measurement-Control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China; 3. Geophysical Exploration Technology Research Center, BGP Inc, CNPC, Zhuozhou, Hebei 072751, China; 4. Exploration and Development Research Institute of CNOOC Research Institute Co., Ltd, Beijing 100028, China; 5. National Engineering Research Center of Offshore Oil and Gas Exploration, Beijing 100028, China; 6. National Key Laboratory of Petroleum Resources and Engineering, China University of Petroleum (Beijing), Beijing 102249, China; 7. CNPC Key Laboratory of Geophysical Exploration, China University of Petroleum (Beijing), Beijing 102249, China; 8. College of Geophysics, China University of Petroleum (Beijing), Beijing 102249, China
Abstract:Full-waveform inversion fully utilizes all the information of seismic wave propagation, including amplitude, phase, and travel time, which can obtain inversion results with higher resolution and accuracy compared with travel time tomography. When the velocity in the shallow layer of the medium is low, fine grids are usually used to sample the low-velocity layer, so as to ensure the accuracy of forward modeling. However, subdividing the entire model into finer grids will increase computational complexity and storage requirements, and oversampling will occur in the high-velocity regions of the model. These problems are further amplified in the process of full-waveform inversion. To address these issues, this study introduced a finite difference forward modeling method of the elastic wave equation in variable grids based on variable difference coefficients. Firstly, the difference coefficient of variable grid wave field simulation was derived based on Taylor expansion, and the variable grid wave field simulation was realized. Secondly, the forward modeling method of variable difference coefficients was applied to the forward modeling, residual backpropagation, and wave field reconstruction of the full-waveform inversion, and the full-waveform inversion of elastic waves in variable grids based on the variable difference coefficient was realized. In the full-waveform inversion, the multi-scale inversion strategy and the conventional conjugate gradient method were used for the iterative solution. The fine grids were used to divide shallow low-velocity layers, and coarser grids were used to divide the middle and deep layers with high velocity, which ensured the accuracy of inversion for the shallow layers while avoiding oversampling in the middle and deep layers. The inversion results of the model data demonstrate that the full-waveform inversion in variable grids based on variable difference coefficients can effectively characterize the low-velocity anomalies compared with full-waveform inversion in uniform coarse grids. The test of noisy data shows that the proposed full-waveform inversion method has strong noise resistance.
QIN Ning, LIANG Hongxian, GUO Zhenbo, et al. A parallel implementing strategy for full waveform inversion of 3D elastic waves based on domain decomposition[J]. Oil Geophysical Prospecting,2023, 58(2):351-357.秦宁, 梁鸿贤, 郭振波, 等. 基于区域分解的3D弹性波全波形反演并行实现策略[J]. 石油地球物理勘探, 2023, 58(2):351-357.
[2]
YANG Ruidong, HUANG Jianping, YANG Zhenjie, et al. Multi-scale full waveform inversion based on double diagonal flux correction transport[J]. Oil Geophysical Prospecting, 2022, 57(5):1120-1128.杨瑞冬, 黄建平, 杨振杰, 等. 基于双对角通量校正的多尺度全波形反演[J]. 石油地球物理勘探, 2022, 57(5):1120-1128.
[3]
LI Qingyang, WU Guochen, WANG Yumei, et al. Elastic full-waveform inversion of land single-component seismic data based on optimal transport theory[J].Oil Geophysical Prospecting, 2021, 56(5):1060-1073.李青阳, 吴国忱, 王玉梅, 等. 基于最优输运原理的陆上单分量资料弹性波全波形反演[J]. 石油地球物理勘探, 2021, 56(5):1060-1073.
[4]
CHEN Liang, HUANG Jianping, WANG Ziying, et al. Seismic numerical simulation and reverse time migration with adaptive variable-grid and compact difference method[J].Oil Geophysical Prospecting, 2023,58(3):641-650, 712.陈亮, 黄建平, 王自颖, 等. 应用自适应变网格紧致差分的地震波数值模拟及逆时偏移[J]. 石油地球物理勘探, 2023, 58(3):641-650, 712.
[5]
MOCZO P. Finite difference technique for SH waves in 2D media using irregular grid:Application to the seismic response problem[J]. Geophysical Journal International, 1989, 99(2):321-329.
[6]
FALK J, TESSMER E, GAJEWSKI D. Tube wave modeling by the finite-difference method with varying grid spacing[J]. Pure and Applied Geophysics, 1996, 148(1):77-93.
[7]
FALK J, TESSMER E, GAJEWSKI D. Efficient finite-difference modelling of seismic waves using locally adjustable time steps[J]. Geophysical Prospec-ting, 1998, 46(6):603-616.
[8]
MUFTI I R, PITA J A, HUNTLEY R W. Finite-difference depth migration of exploration-scale 3-D seismic data[J]. Geophysics, 1996, 61(3):776-794.
[9]
WANG Y, XU J, SCHUSTER G T. Viscoelastic wave simulation in basins by a variable-grid finite-difference method[J]. Bulletin of the Seismological Society of America, 2001, 91(6):1741-1749.
[10]
JASTRAM C, BEHLE A. Acoustic modeling on a grid of velocity varying spacing[J]. Geophysical Prospecting, 1992, 40(2):157-169.
[11]
JASTRAM C, TESSMER E. Elastic modelling on a grid with vertically varying spacing[J]. Geophysical Prospecting, 1994, 42(4):357-370.
[12]
SUN Chengyu, DING Yucai. Accuracy analysis of wave equation finite difference with dual-variable grid algorithm[J]. Oil Geophysical Prospecting, 2012, 47(4):545-551.孙成禹, 丁玉才. 波动方程有限差分双变网格算法的精度分析[J]. 石油地球物理勘探, 2012, 47(4):545-551.
[13]
TESSMER E. Seismic finite-difference modeling with spatially varying time steps[J]. Geophysics, 2000, 65(4):1290-1293.
[14]
ADRIANO S, OLIVEIRA M. A fourth-order finite-difference method for the acoustic wave equation on irregular grids[J]. Geophysics, 2003, 68(2):672-676.
[15]
KANG T S, BAAG C E. An efficient finite-difference method for simulating 3D seismic response of localized basin structures[J]. Bulletin of the Seismological Society of America, 2004, 94(5):1690-1705.
[16]
KANG T S, BAAG C E. Finite-difference seismic simulation combining discontinuous grids with locally variable time steps[J]. Bulletin of the Seismological Society of America, 2004, 94(1):207-219.
[17]
SUN Weitao, YANG Huizhu. A 3-D finite difference method using irregular grids for elastic wave propagation in anisotropic media[J]. Chinese Journal of Geophysics, 2004, 47(2):332-337.孙卫涛,杨慧珠. 各向异性介质弹性波传播的三维不规则网格有限差分方法[J]. 地球物理学报, 2004, 47(2):332-337.
[18]
NIE S, WANG Y, OLSEN K, et al. Fourth-order staggered-grid finite-difference seismic wavefield estimation using a discontinuous mesh interface[J]. Bulletin of the Seismological Society of America, 2017, 107(5):2183-2193.
[19]
KRISTEK J, MOCZO P, GALIS M. Stable discontinuous staggered grid in the finite-difference modelling of seismic motion[J]. Geophysical Journal International, 2010, 183(3):1401-1407.
GAO L, BROSSIER R, VIRIEUX J. Using time filtering to control the long-time instability in seismic wave simulation[J]. Geophysical Journal International. 2016, 204(3):1443-1461.
[22]
ZHANG Z, ZHANG W, LI H, et al. Stable discontinuous grid implementation for collocated-grid finite-difference seismic wave modelling[J]. Geophysical Journal International, 2012, 192(3):1179-1188.
[23]
ZHANG Jianfeng. Non-orthogonal grid finite-difference method for numerical simulation of elastic wave propagation[J]. Chinese Journal of Geophysics, 1998, 41(S1):357-366.张剑锋. 弹性波数值模拟的非规则网格差分法[J]. 地球物理学报, 1998, 41(增刊1):357-366.
[24]
WU C, 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[J]. Geophysics, 2005, 70(4):T57-T62.
[25]
ZHU Shengwang, WEI Xiucheng. Differential forward modeling of wave-equation having irregular grid and any-order precision[J]. Oil Geophysical Prospecting, 2005, 40(2):149-153.朱生旺, 魏修成. 波动方程非规则网格任意阶精度差分法正演[J]. 石油地球物理勘探, 2005, 40(2):149-153.
[26]
ZHANG W, ZHANG J. Full waveform tomography with consideration for large topography variations[C]. SEG Technical Program Expanded Abstracts, 2011, 30:2539-2542.
[27]
HUANG Chao, DONG Liangguo. Staggered-grid high-order finite-difference method in elastic wave simulation with variable grids and local time-steps[J]. Chinese Journal of Geophysics, 2009, 52(11):2870-2878.黄超, 董良国. 可变网格与局部时间步长的交错网格高阶差分弹性波模拟[J]. 地球物理学报, 2009, 52(11):2870-2878.
[28]
ZHANG Hui, LI Zhenchun. Forward modeling of carbonate fracture reservoir based on time-space dual variable grid algorithm[J]. Journal of China University of Petroleum, 2011, 35(3):51-57.张慧, 李振春. 基于时空双变网格算法的碳酸盐岩裂缝型储层正演模拟[J]. 中国石油大学学报, 2011, 35(3):51-57.
[29]
PITARKA A. 3D elastic finite-difference modeling of seismic motion using staggered grids with nonuniform spacing[J]. Bulletin of Seismological Society of Ame-rica, 1999, 89(1):54-68.