摘要
间断有限元(Discontinuous Galerkin:DG)方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D′Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D′Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D′Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.
The Discontinuous Galerkin （DG） method has great advantages in suppressing numerical dispersion and dealing with complex structures. Therefore, in this paper, we apply a new DG method to numerical simulations in two-phase and viscoelastic media, and suggest a DG method to solve both Biot elastic wave equations and the D＇Alembert wave equations. For this, we first transform the Biot equations and the D＇Alembert wave equations into a system of first-order equations with respect to time-space by introducing auxiliary variables. Then we transform the first-order equations into a semi-discrete ordinary differential equation （ODE） system using the DG method.Finally, we use a weighted Runge-Kutta method show that the DG method works very well for Alembert wave equations, and can effectively accurate information on the wave-field. to solve the ODE system. The numerical results solving the Biot elastic wave equations and Dr suppress the numerical dispersion and provide
作者
张金波
杨顶辉
贺茜君
马啸
ZHANG JinBo1 , YANG DingHui1, HE XiJun2 , MA Xiao3(Department of Mathematical Sciences, Tsinghua University, Beij ing 100084, China Department of Mathematics, Colledge of Information Science and Technology, Hainan University, Haikou, 570228, China Department of Applied Mathematics, Northwestern Polytechnical University, Xi＇an 710072, Chin)
出处
《地球物理学报》
SCIE
EI
CAS
CSCD
北大核心
2018年第3期926-937,共12页
Chinese Journal of Geophysics
基金
国家自然科学基金重大项目(41390452)的资助
关键词
数值模拟
间断有限元
双相介质
黏弹性介质
Numerical simulation
Discontinuous Galerkin method
Two-phase medium
Viscoelasticmedium
