文献类型：期刊文章

作　　者：张金波[1] 杨顶辉[1] 贺茜君[2] 马啸[3]

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)

机构地区：[1]清华大学数学科学系,北京100084 [2]海南大学信息学院数学系,海口570228 [3]西北工业大学应用数学系,西安710072

出　　处：《地球物理学报》

基　　金：国家自然科学基金重大项目(41390452)的资助

年　　份：2018

卷　　号：61

期　　号：3

起止页码：926-937

语　　种：中文

收录情况：AJ、BDHX、BDHX2017、CAS、CSA、CSA-PROQEUST、CSCD、CSCD2017_2018、EI、GEOBASE、GEOREFPREVIEWDATABASE、IC、JST、PA、RCCSE、SCIE、SCOPUS、WOS、ZGKJHX、核心刊

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

关 键 词：数值模拟  间断有限元 双相介质 黏弹性介质

分 类 号：P315] P631]