Expanded approximate analytic discretization and elastic wave numerical simulation
Wang Yong1,2, Duan Yanwen1,2, An Yifan1,2, Wang Xiaocong3, Gui Zhixian1,2
1. Geophysics and Oil Resource Institute, Yangtze University, Wuhan, Hubei 430100, China; 2. Key Laboratory of Exploration Technologies for Oil and Gas Resources, the Ministry of Education, Yangtze University, Wuhan, Hubei 430100, China; 3. BGP International, BGP Inc., CNPC, Zhuozhou, Hebei 072751, China
Abstract:The seismic wave numerical simulation is the foundation of the seismic wave theory,migration imaging and seismic inversion.The accuracy of numerical simulation has a great significance for seismic exploration.This paper discusses an expanded approximate analytic discretization method for seismic wave numerical simulation,and theoretically analyses its accuracy,dispersion relations,stability,and computational efficiency.This method is one order higher in time difference than before,thereby the calculation accuracy is improved.Simulation results demonstrate that the maximum error rate is 88% lower than before.As other approximate analytic discretization methods,the proposed method has advantages of shorter operator radius and better adaption to too-coarse grids.In the smallest main wavelength it needs only 5.9 grid points.Compared with the dispersion curves and the simulation results with the fourth-order LWC and the staggered-grid finite difference scheme,the proposed method can better suppress numerical dispersion.Finally with the proposed method we numerically simulate elastic waves in the isotropic homogeneous medium model and horizontal layered medium model.Very clear and accurate seismic characteristics demonstrate the practicability of the proposed method.
Cerveny V.Seismic rays and ray intensities in inhomogeneous anisotropic media.Geophyscial Journal International,1972,29(1):1-33.
[2]
Kosloff D D,Baysal E.Forward modeling by a Fourier method.Geophysics,1982,47(10):1402-1412.
[3]
唐怀谷,何兵寿.一阶声波方程时间四阶精度差分格式的伪谱法求解.石油地球物理勘探,2017,52(1):71-80.Tang Huaigu,He Bingshou.Pseudo spectrum method of first-order acoustic wave equation finite-difference schemes with fourth-order time difference accuracy.OGP,2017,52(1):71-80.
[4]
Drake L A.Rayleigh waves at a continental boundary by the finite element method.Bulletin of the Seismological Society of America,1972,62(5):1259-1268.
[5]
Bouchon M,Coutant O.Calculation of synthetic seismograms in a laterally varying medium by the bounda-ry element-discrete wavenumber method.Bulletin of the Seismological Society of America,1994,84(6):1869-1881.
[6]
Komatissch D,Barnes C,Tromp J.Simulation of anisotropic wave propagation based upon a spectral element method.Geophysics,2000,65(4):1251-1260.
[7]
Alterman Z,Karal F C.Propagation of seismic wave in layered media by finite difference methods.Bulletin of the Seismological Society of America,1968,58(1):367-398.
[8]
Kelly K R,Ward R W,Treitel S et al.Synthetic seismograms:a finite-difference approach.Geophysics,1976,41(1):2-27.
[9]
Virieux J.P-SV wave propagation in heterogeneous media:Velocity-stress finite-difference method.Geo-physics,1986,51(4):889-901.
[10]
Graves R W.Simulation seismic wave propagation in 3D elastic media using staggered-grid finite differences.Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
[11]
董良国,马在田,曹景忠等.一阶弹性波方程交错网格高阶差分解法.地球物理学报,2000,43(3):411-419.Dong Liangguo,Ma Zaitian,Cao Jingzhong et al.The staggered-grid high-order difference method of first-order elastic equation.Chinese Journal of Geophysics,2000,43(3):411-419.
[12]
董良国,马在田,曹景忠.一阶弹性波方程交错网格高阶差分解法稳定性研究.地球物理学报,2000,43(6):856-864.Dong Liangguo,Ma Zaitian,Cao Jingzhong.The stability study of the staggered-grid high-order difference method of first-order elastic equation.Chinese Journal of Geophysics,2000,43(6):856-864.
[13]
董良国.弹性波数值模拟中的吸收边界条件.石油地球物理勘探,1999,34(1):46-56.Dong Liangguo.Absorptive boundary condition in elastic-wave numerical modeling.OGP,1999,34(1):46-56.
[14]
王书强,杨顶辉,杨宽德.弹性波方程的紧致差分方法.清华大学学报(自然科学版),2002,42(8):1128-1131.Wang Shuqiang,Yang Dinghui,Yang Kuande.Com-pact finite difference scheme for elastic equations.Journal of Tsinghua University(Science and Techno-logy),2002,42(8):1128-1131.
[15]
Saenger E H,Gold N,Shapiro S A.Modeling the propa-gation of elastic waves using a modified finite-difference grid.Wave Motion,2000,31(1):77-92.
[16]
徐文才,杨国权,李振春等.横向各向同性介质拟声波一阶速度-应力方程.石油地球物理勘探,2016,51(1):87-96.Xu Wencai,Yang Guoquan,Li Zhenchun et al.First order velocity-stress equation in TI media.OGP,2016,51(1):87-96.
[17]
Kondoh Y,Hosaka Y,Ishii K.Kerrel optimum nearly analytical discrimination algorithm applied to parabo-lic and hyperbolic equations.Computers & Mathema-tics with Applications,1994,27(3):59-90.
[18]
杨顶辉,滕吉文,张中杰.三分量地震波场的近似解析离散模拟技术.地球物理学报,1996,39(增刊):283-291.Yang Dinghui,Teng Jiwen,Zhang Zhongjie.Nearly-analytical discretization modeling technique of 3-component seismic wave-fields.Chinese Journal of Geophysics,1996,39(S):283-291.
[19]
Yang D H,Teng J W,Zhang Z J et al.A nearly-analytic discrete method for acoustic and elastic wave equation.Bulletin of the Seismological Society of America,2003,93(2):882-890.
[20]
Yang D H,Lu M,Wu R S et al.An optimal nearly-ana-lytic discrete method for 2D acoustic and elastic wave equations.Bulletin of the Seismological Society of America,2004,94(5):1982-1991.
[21]
Yang D H,Song G J,Lu M.Optimally accurate nearly analytic discrete scheme for wave-field simulation in 3D anisotropic media.Bulletin of the Seismological Society of America,2007,97(5):1557-1569.
[22]
Yang D H,Song G J,Chen S et al.An improved nearly analytical discrete method:an efficient tool to simulate the seismic response of 2-D porous structures.Journal of Geophysics & Engineering,2007,4(1):40-52.
Tong P,Yang D H,Hua B L.High accuracy wave simulation-revised derivation,numerical analysis and testing of a nearly analytic integration discrete method for solving acoustic equation.International Journal of Solid and Structures,2011,48(1):56-70.
[26]
Lukacova M M,Warnecke G.Lax-Wendroff type se-cond order evolution Galerkin methods for multidimensional hyperbolic systems.East-West Journal of Numerical Mathematics,2000,8(2):127-152.
[27]
董良国,李培明.地震波传播数值模拟中的频散问题.天然气工业,2004,24(6):53-56.Dong Liangguo,Li Peiming.Dispersive problem in seismic wave propagation numerical modeling.Natural Gas Industry,2004,24(6):53-56.
[28]
张文生.科学计算中的偏微分方程有限差分法.北京:高等教育出版社,2006.
[29]
高正辉,孙建国,孙章庆等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟.石油地球物理勘探,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 perfectly matched layer.OGP,2016,51(6):1128-1133.
[30]
汪勇,段焱文,王婷等.优化近似解析离散化方法的二维弹性波波场分离模拟.石油地球物理勘探,2017,52(3):458-467.Wang Yong,Duan Yanwen,Wang Ting et al.Numerical simulation of elastic wave separation in 2D isotropic medium with the optimal nearly-analytic discretization.OGP,2017,52(3):458-467.