Elastic wavefield forward modeling in heterogeneous media based on the quasi-regular grid high-order finite difference
LI Qingyang1,2, WU Guochen1,2, DUAN Peiran1
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266071, China
Abstract:In this paper, we propose a quasi-regular grid high-order finite-difference seismic forward mo-deling to accurately describe elastic wave propagation in heterogeneous media.The quasi-regular grid strategy is to stagger the displacement component of the elastic wave equation and make a new center difference operator for displacement, while the si-mulation accuracy and stability are the same as those of the staggered grid.But the memory usage of the quasi-regular grid is reduced by 60% in 2D and 66.7% in 3D compared to the staggered grid.Then the equivalence between the quasi-regular grid and staggered grid is mathematically demonstrated, and the source loading, boundary conditions and stability are analyzed.A numerical test is carried out on a layered model to validate the accuracy of the quasi-regular grid method by comparing with staggered grid and regular grid methods.Finally, a test on Marmousi-2 model proves the applicability and stability of the proposed method.
Pestana R C, Stoffa P L.Time evolution of the wave equation using rapid expansion method[J].Geophy-sics, 2010, 75(4):T121-T131.
[4]
Tessmer E.Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration[J].Geophysics, 2011, 76(4):S177-S185.
[5]
Fomel S, Ying L, Song X.Seismic wave extrapolation using low rank symbol approximation[J].Geophysical Prospecting, 2013, 61(3):526-536.
[6]
Wu Z, Alkhalifah T.The optimized expansion based low-rank method for wavefield extrapolation[J].Geophysics, 2014, 79(2):T51-T60.
[7]
Fang G, Fomel S, Du Q, et al.Low rank seismic-wave extrapolation on a staggered grid[J].Geophysics, 2014, 79(3):T157-T168.
[8]
Kelly K R, Ward R W, Treitel S, et al.Synthetic seismograms:A finite-difference approach[J].Geophy-sics, 1976, 41(1):2-27.
[9]
Robertsson J O A, Blanch J O, Symes W W, et al.Galerkin-wavelet modeling of wave propagation:Optimal finite-difference stencil design[J].Mathematical and Computer Modelling, 1994, 19(1):31-38.
[10]
Robertsson J O A, Blanch J O, Symes W W.Visco-elastic finite-difference modeling[J].Geophysics, 1994, 59(9):1444-1456.
[11]
Alterman Z, Karal Jr F C.Propagation of elastic waves in layered media by finite difference methods[J].Bu-lletin of the Seismological Society of America, 1968, 58(1):367-398.
[12]
Boore D M.Finite difference methods for seismic wave propagation in heterogeneous materials[J].Methods in Computational Physics, 1972, 11:1-37.
[13]
Madariaga R.Dynamics of an expanding circular fault[J].Bulletin of the Seismological Society of America, 1976, 66(3):639-666.
[14]
Virieux J.SH-wave propagation in heterogeneous media:Velocity-stress finite-difference method[J].Geophysics, 1984, 49(11):1933-1942.
Levander A R.Fourth-order finite-difference P-SV seismograms[J].Geophysics, 1988, 53(11):1425-1436.
[17]
Graves R W.Simulating seismic wave propagation in 3D elastic media using staggered-grid finite diffe-rences[J].Bulletin of the Seismological Society of America, 1996, 86(4):1091-1106.
[18]
Moczo P, Kristek J, Halada L.3D fourth-order staggered-grid finite-difference schemes:Stability and grid dispersion[J].Bulletin of the Seismological Society of America, 2000, 90(3):587-603.
[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.
[20]
刘洋, 李承楚, 牟永光.任意偶数阶精度有限差分法数值模拟[J].石油地球物理勘探, 1998, 33(1):1-10.LIU Yang, LI Chengchu, MOU Yongguang.Finite-difference numerical modeling of any even-order accuracy[J].Oil Geophysical Prospecting, 1998, 33(1):1-10.
[21]
董良国, 马在田, 曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报, 2000, 43(6):856-864.DONG Liangguo, MA Zaitian, CAO Jingzhong.A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation[J].Chinese Journal of Geophysics, 2000, 43(6):856-864.
[22]
殷文, 印兴耀, 吴国忱, 等.高精度频率域弹性波方程有限差分方法及波场模拟[J].地球物理学报, 2006, 49(2):561-568.YIN Wen, YIN Xingyao, WU Guochen, et al.The method of finite difference of high precision on elastic wave equations in the frequency domain and wavefield simulation[J].Chinese Journal of Geophysics, 2006, 49(2):561-568.
[23]
杨庆节, 刘财, 耿美霞, 等.交错网格任意阶导数有限差分格式及差分系数推导[J].吉林大学学报(地球科学版), 2014, 44(1):375-385.YANG Qingjie, LIU Cai, GENG Meixia, et al.Staggered grid finite difference scheme and coefficients deduction of any number of derivatives[J].Journal of Jilin University(Earth Science Edition), 2014, 44(1):375-385.
[24]
杜启振, 郭成锋, 公绪飞.VTI介质纯P波混合法正演模拟及稳定性分析[J].地球物理学报, 2015, 58(4):1290-1304.DU Qizhen, GUO Chengfeng, GONG Xufei.Hybrid PS/FD numerical simulation and stability analysis of pure P-wave propagation in VTI media[J].Chinese Journal of Geophysics, 2015, 58(4):1290-1304.
[25]
黄金强, 李振春, 黄建平, 等.TTI介质Lebedev网格高阶有限差分正演模拟及波型分离[J].石油地球物理勘探, 2017, 52(5):915-927.HUANG Jinqiang, LI Zhenchun, HUANG Jianping, et al.Lebedev grid high-order finite-difference mode-ling and elastic wave-mode separation for TTI media[J].Oil Geophysical Prospecting, 2017, 52(5):915-927.
[26]
高正辉, 孙建国, 孙章庆, 等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟[J].石油地球物理勘探, 2016, 51(6):1128-1133.GAO Zhenghui, SUN Jianguo, SUN Zhangqing, et al.2.5D acoustic approximate equation numerical simulation with a new construction method of the perfect matched layer[J].Oil Geophysical Prospecting, 2016, 51(6):1128-1133.
[27]
梁文全, 王彦飞, 杨长春.声波方程数值模拟矩形网格有限差分系数确定法[J].石油地球物理勘探, 2017, 52(1):56-62.LIANG Wenquan, WANG Yanfei, YANG Changchun.Acoustic wave equation modeling with rectangle grid finite difference operator and its linear time space domain solution[J].Oil Geophysical Prospecting, 2017, 52(1):56-62.
[28]
Di Bartolo L, Dors C, Mansur W J.A new family offinite-difference schemes to solve the heterogeneous acoustic wave equation[J].Geophysics, 2012, 77(5):T187-T199.
[29]
Dablain M A.The application of high-order differencing to the scalar wave equation[J].Geophysics, 1986, 51(1):54-66.
[30]
李青阳, 吴国忱, 梁展源.基于PML边界的时间高阶伪谱法弹性波场模拟[J].地球物理学进展, 2018, 33(1):228-235.LI Qingyang, WU Guochen, LIANG Zhanyuan.Time domain high-order pseudo spectral method based on PML boundary for elastic wavefield simulation[J].Progress in Geophysics, 2018, 33(1):228-235.
[31]
吴建鲁, 吴国忱, 王伟, 等.基于等效交错网格的流固耦合介质地震波模拟[J].石油地球物理勘探, 2018, 53(2):272-279.WU Jianlu, WU Guochen, WANG Wei, et al.Seismic forward modeling in fluid-solid media based on equi-valent staggered grid scheme[J].Oil Geophysical Prospecting, 2018, 53(2):272-279.
[32]
Reynolds A C.Boundary conditions for the numerical solution of wave propagation problems[J].Geophy-sics, 2012, 43(6):155-165.
[33]
Clayton R, Engquist B.Absorbing boundary condi-tions for acoustic and elastic wave equations[J].Bulletin of the Seismological Society of America, 1977, 67(6):1529-1540.
[34]
Cerjan C, Dan K, Kosloff R, et al.A nonreflecting boundary condition for discrete acoustic and elastic wave equations[J].Geophysics, 1985, 50(4):705-708.
[35]
Berenger J P.A perfectly matched layer for the ab-sorption of electromagnetic waves[J].Physics of Plasmas, 1994, 114(2):185-200.