Seismic wave modeling based on modified symplectic scheme and quasi-particles method
Su Bo1,3, Tuo Xianguo2, Liu Zhigui3
1. Graduate School, CAEP, Mianyang, Sichuan 621900, China; 2. Sichuan University of Science & Engineering, Zigong, Sichuan 643002, China; 3. School of Computer Science and Technology, Southwest University of Science and Technology, Mianyang, Sichuan 621010, China
Abstract:Based on Hamilton mechanics,we develop a quasi-particles system to discretize seismic wave equation from the viewpoint of molecular dynamics.Each particle in the system only interacts with the particles which located at the upper,the lower,the left,and the right,and the four diagonal particles.The force and the relative displacement between the particles can be considered as approximately linear.In this paper,the interaction coefficients of the particles are derived.We get a modified symplectic scheme with third-order on the basis of two-order symplectic schemes to deal with temporal discretization.Theoretical analysis shows the new scheme possesses weaker numerical dispersion and larger stability than those of the conventional symplectic schemes.The modified symplectic scheme is suitable for long-term computational owing to all positive symplectic coefficients which consistent with iterative computation.In order to test the accuracy of the symplectic scheme and the spatial quasi-particles system,we adopt Lamb's problem to investigate the accuracy and efficiency of the numerical simulation of elastic wave.In addition,the Sigsbee2B velocity model is selected to test the stability of the proposed method.Numerical experiments demonstrate improvements in the numerical dispersion suppression and numerical stability of the modified symplectic scheme.
曹禹,杨孔庆.对声波和弹性波传播模拟的Hamilton系统方法.物理学报,2003,52(8):1984-1992.Cao Yu,Yang Kongqing.Hamiltonian system approach for simulation of acoustic and elastic wave propagation.Acta Physica Sinica,2003,52(8):1984-1992.
[2]
Alford R,Kelly K,Boore D M.Accuracy of finite-diffe-rence.Geophysics,1974,39(6):834-842.
[3]
Virieux J.P-SV wave propagation in heterogeneous media:velocity-stress finite-difference method.Geophysics,1986,51(4):889-901.
[4]
Lysmer J,Drake L A.A finite element method for seismology.Methods in Computational Physics Advances in Research & Applications,1972,11:181-216.
[5]
Moser F,Jacobs L J,Qu J.Modeling elastic wave propagation in waveguides with the finite element method.Ndt & E International,1999,32(4):225-234.
[6]
Liu S,Li X,Wang W et al.A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modelling.Journal of Geophysics and Engineering,2014,11(5):055009.
[7]
Fornberg B.High-order finite difference and the pseudospectral method on staggered grids.Siam Journal on Numerical Analysis,1990,27(4):904-918.
[8]
马德堂,朱光明.有限元法与伪谱法混合求解弹性波动方程.地球科学与环境学报,2004,26(1):61-64.Ma Detang,Zhu Guangming.Hybrid method combining finite difference and pseudospectral method for solving the elastic wave equation.Journal of Earth Sciences and Environment,2004,26(1):61-64.
[9]
唐怀谷,何兵寿.一阶声波方程时间四阶精度差分格式的伪谱法求解.石油地球物理勘探,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.
[10]
Komatitsch D,Tromp J.The spectral-element me-thod,Beowulf computing,and global seismology.Science,2002,298(5599):1737-1742.
[11]
汪文帅.基于辛-谱元方法的地震波场模拟研究[学位论文].北京:中国科学院大学,2013.
[12]
Yang D H,Teng J,Zhang Z et al.A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media.BSSA,2003,93(2):882-890.
[13]
Chen Y,Song G,Xue Z et al.An improved optimal nearly analytical discretized method for 2D scalar wave equation in heterogeneous media based on the modified nearly analytical discrete operator.Geophy-sics,2014,79(6):T349-T362.
[14]
汪勇,段焱文,王婷等.优化近似解析离散化方法的二维弹性波波场分离模拟.石油地球物理勘探,2017,52(3):458-467.Wang Yong,Duan Yanwen,Wang Ting et al.Numerical simulation of elactic wave separation in 2D isotropic medium with the optimal nearly-analytic discretization.OGP,2017,52(3):458-467.
[15]
Takekawa J,Madariage R,Mikada H et al.Numerical simulation of seismic wave propagation produced by earthquake by using a particle method.Geophysical Journal International,2012,191(3):1305-1316.
[16]
Suzuki Y,Koshizuka S,Oka Y.Hamiltonian moving-particle semi-implicit (HMPS) method for incompressible fluid flows.Computer Methods in Applied Mechanics and Engineering,2007,196(29-30):2876-2894.
[17]
Liu S,Yang D,Ma J.A modified symplectic PRK scheme for seismic wave modeling.Computers & Geosciences,2017,99:28-36.
[18]
罗明秋,刘洪,李幼铭.地震波传播的哈密顿表述及辛几何算法.地球物理学报,2001,44(1):120-128.Luo Mingqiu,Liu Hong,Li Youming.Hamiltonian description and symplectic method of seismic wave propagation.Chinese Journal of Geophysics,2001,44(1):120-128.
[19]
Chen J B.High-order time discretizations in seismic modeling.Geophysics,2007,72(5):SM115-SM122.
[20]
Liu S,Li X,Wang W et al.A new kind of optimal se-cond-order symplectic scheme for seismic wave simulations.Science in China:Earth Sciences,2014,57(4):751-758.
[21]
Liu S,Yang D,Lang C et al.Modified symplectic schemes with nearly-analytic discrete operators for acoustic wave simulations.Computer Physics Communications,2017,213:52-63.
[22]
Liu S,Li X,Wang W et al.A modified symplectic scheme for seismic wave modeling.Journal of Applied Geophysics,2015,116(5):110-120.
[23]
冯康.哈密尔顿系统的辛几何算法.浙江杭州:浙江科技出版社,2003.
[24]
马毅恒,井西利,邢小宁等.弹性波的三维Hamilton系统方法模拟.石油地球物理勘探,2012,47(4):524-531.Ma Yiheng,Jing Xili,Xing Xiaoning et al.3D Hamiltonian system approach simulation for elastic wave.OGP,2012,47(4):524-531.
[25]
刘少林,李小凡,汪文帅等.求解弹性波方程的辛RKN格式.地球物理学报,2015,58(4):1355-1366.Liu Shaolin,Li Xiaofan,Wang Wenshuai.A symplectic RKN scheme for solving elastic wave equations.Chinese Journal of Geophysics,2015,58(4):1355-1366.
[26]
Li X,Li Y,Zhang M et al.Scalar seismic-wave equation modeling by a multisymplectic discrete singular convolution differentiator method.BSSA,2011,101(4):1710-1718.
[27]
Basabe J D D,Sen M K.Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations.Geophysics,2007,72(6):T81-T95.
[28]
梁文全,王彦飞,杨长春.声波方程数值模拟矩形网格有限差分系数确定法.石油地球物理勘探,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.OGP,2017,52(1):56-62.
[29]
De Hoop A.A modification of Cagniard's method for solving seismic pulse problems.Applied Scientific Research,Section B,1960,8(1):349-356.