张金波^1,,
杨顶辉^1,,,
贺茜君²,
马啸³
1. 清华大学数学科学系, 北京 100084
2. 海南大学信息学院数学系, 海口 570228
3. 西北工业大学应用数学系, 西安 710072

基金项目: 国家自然科学基金重大项目（41390452）的资助


详细信息

作者简介: 张金波, 男, 1991年生, 清华大学数学科学系硕士研究生, 主要研究方向为地震波动方程的数值方法及波场模拟.E-mail:zhangjb14@mails.tsinghua.edu.cn

通讯作者: 杨顶辉, 教授, 主要从事计算地球物理、孔隙介质波传播理论、地震层析成像等研究.E-mail:dhyang@math.tsinghua.edu.cn

中图分类号: P315;P631

收稿日期:2017-02-13
修回日期:2017-06-19
上线日期:2018-03-05

Discontinuous Galerkin method for solving wave equations in two-phase and viscoelastic media

ZHANG JinBo^1,,
YANG DingHui^1,,,
HE XiJun²,
MA Xiao³
1. Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China
2. Department of Mathematics, Colledge of Information Science and Technology, Hainan University, Haikou, 570228, China
3. Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072, China


More Information

Corresponding author: YANG DingHui,E-mail:dhyang@math.tsinghua.edu.cn

MSC: P315;P631

--> Received Date: 13 February 2017
Revised Date: 19 June 2017
Available Online: 05 March 2018


摘要
摘要:间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此，本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟，发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组，然后对该方程组进行DG空间离散，得到半离散化的常微分方程组.最后，对此常微分方程组，应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明，DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程，并能很好地压制因离散求解波动方程而产生的数值频散，获得清晰的各种地震波震相.
关键词: 数值模拟/
间断有限元/
双相介质/
黏弹性介质

Abstract: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 to solve the ODE system. The numerical results show that the DG method works very well for solving the Biot elastic wave equations and D'Alembert wave equations, and can effectively suppress the numerical dispersion and provide accurate information on the wave-field.
Key words:Numerical simulation/
Discontinuous Galerkin method/
Two-phase medium/
Viscoelastic medium



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14416