苏波^1,2,,
李怀良³,
刘少林^4,,
杨顶辉⁴
1. 西南科技大学计算机科学与技术学院, 四川绵阳 621010
2. 中国工程物理研究院研究生院, 四川绵阳 621900
3. 核废物与环境安全国防重点学科实验室, 四川绵阳 621010
4. 清华大学数学科学系, 北京 100084

基金项目: 国家自然科学基金（41604034，41774118），国家重点研发计划项目课题（2017YFC150031）资助


详细信息

作者简介: 苏波, 男, 1977年生, 西南科技大学计算机科学与技术学院教师, 主要从事地震波正反演的研究工作.E-mail:bosu@foxmail.com
刘少林, 男, 1988年生, 新加坡南洋理工大学、清华大学数学系博士后, 从事地震波传播与成像研究.E-mail:shaolinliu88@gamil.com

中图分类号: P315

收稿日期:2018-09-03
修回日期:2018-12-02
上线日期:2019-04-05

Modified symplectic scheme with finite element method for seismic wavefield modeling

SU Bo^1,2,,
LI HuaiLiang³,
LIU ShaoLin^4,,
YANG DingHui⁴
1. School of Computer Science and Technology, Southwest University of Science and Technology, Mianyang Sichuan 621010, China
2. Graduate School, China Academy of Engineering Physics, Mianyang Sichuan 621900, China
3. Fundamental Science on Nuclear Wastes and Environmental Safety Laboratory, Mianyang Sichuan 621010, China
4. Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China


MSC: P315

--> Received Date: 03 September 2018
Revised Date: 02 December 2018
Available Online: 05 April 2019


摘要
摘要:三角网格有限元法具有网格剖分的灵活性，能有效模拟地震波在复杂介质中的传播.但传统有限元法用于地震波场模拟时计算效率较低，消耗较大计算资源.本文采用改进的核矩阵存储（IKMS）策略以提高有限元法的计算效率，该方法不用组合总体刚度矩阵，且相比于常规有限元法节省成倍的内存.对于时间离散，将有限元离散后的地震波运动方程变换至Hamilton体系，在显式二阶辛Runge-Kutta-Nystr m（RKN）格式的基础之上加入额外空间离散算子构造修正辛差分格式，通过Taylor展开式得到具有四阶时间精度时间格式，且辛系数全为正数.本文从理论上分析了时空改进方法相比传统辛-有限元方法在频散压制、稳定性提升等方面的优势.数值算例进一步证实本方法具有内存消耗少、稳定性强和数值频散弱等优点.
关键词: 有限元法/
修正辛算法/
波场模拟/
数值频散/
稳定性

Abstract:Finite element method (FEM) with triangular mesh has the properties of flexibility and adaptability in complex medium. However, the traditional finite element method is inefficient in seismic wavefield modeling because this method commonly requires a large amount of computation resources. Here we propose a so-called improved kernel matrix storage (IKMS) strategy to improve the computation efficiency of FEM. The improved FEM does not need to assemble global stiffness matrix and the memory requirement is many times less than that of the conventional FEM. In terms of temporal discretization, the elastic wave equation after FEM discretization is first transformed into a Hamilton system. Then, an additional spatial discretization term is added to the second-order explicit Runge-Kutta-Nystr?m (RKN) scheme. Based on Taylor series expansion, we obtain a fourth-order symplectic scheme with all positive symplectic coefficients. Theoretical analysis shows that the temporal-spatial modified numerical scheme is superior over the traditional symplectic FEM in suppressing numerical dispersion and increasing numerical stability. Furthermore, numerical results verify that this method requires less computer resource and achieves higher numerical accuracy.
Key words:Finite element method/
Modified symplectic scheme/
Wavefield simulation/
Numerical dispersion/
Stability



PDF全文下载地址:

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