ISSN 1000-1239 CN 11-1777/TP

计算机研究与发展 ›› 2016, Vol. 53 ›› Issue (5): 1147-1155.doi: 10.7544/issn1000-1239.2016.20148440

• 系统结构 • 上一篇    下一篇

一种求解地震波方程的高效并行谱元格式

林灯,崔涛,冷伟,张林波   

  1. (中国科学院数学与系统科学研究院 北京 100190) (科学与工程计算国家重点实验室(中国科学院数学与系统科学研究院) 北京 100190) (lind@lsec.cc.ac.cn)
  • 出版日期: 2016-05-01
  • 基金资助: 
    国家“九七三”重点基础研究发展计划基金项目(2011CB309703);国家“八六三”高技术研究发展计划基金项目(2012AA01A309);国家自然科学基金项目(11171334,11321061,11101417);中国科学院国家数学与交叉科学中心资助课题

An Efficient Parallel Spectral Element Scheme for Solving Seismic Wave Equation

Lin Deng, Cui Tao, Leng Wei, Zhang Linbo   

  1. (Academy of Mathematics and System Sciences, Chinese Academy of Sciences, Beijing 100190) (State Key Laboratory of Scientific and Engineering Computing (Academy of Mathematics and System Sciences, Chinese Academy of Sciences), Beijing 100190)
  • Online: 2016-05-01

摘要: 地震波数值模拟在地震学和地震勘探中扮演着非常重要角色.在已有工作的基础上,提出1种高效并行的地震波PML方程谱元格式.PML被引入地震波方程以吸收外向波进而模拟无界区域.进一步,为了适应复杂地形同时允许时间显式推进,谱元方法被用来离散地震波PML方程.由此得到地震波PML方程谱元格式.在此基础上,阐述了单元刚度矩阵分解性质,并说明了利用单元刚度矩阵分解可以大幅减少刚度矩阵存储量同时显著加速刚度矩阵与向量乘积,进而显著减少格式的计算量和存储量.此外,算法复杂性分析表明格式无论在计算量上还是在存储量上都优于几种已知的1阶地震波PML方程谱元格式.结合并行技术,给出了高效并行的地震波PML方程谱元格式.数值实验验证了格式的正确性、良好的强弱并行可扩展性以及对复杂地形的适应性.

关键词: 地震波方程, 数值模拟, 完美匹配层, 谱元方法, 单元刚度矩阵分解, 并行计算

Abstract: Numerical simulation of seismic waves plays an essential role in seismology and seismic exploration. We propose here an efficient parallel spectral element scheme for seismic wave equation with perfectly matched layer (PML). PML is integrated into the seismic wave equation to absorb out-going waves and mimic unbounded domain. Ulteriorly, to enable adapting complex topography and explicit time stepping, the spectral element method (SEM) is used to discretize seismic wave equation with PML, which results in a spectral element scheme. In addition, we demonstrate that element stiffness matrices can be decomposed, which can be used to greatly reduce the storage of stiffness matrix and accelerate stiffness matrix-vector multiplication and thus remarkably speed up the scheme and cut down memory cost. Furthermore, we study several spectral element schemes known and show that our scheme is superior to others in both calculation and storage. Combined with parallel technique, an efficient parallel spectral element solver for seismic wave equation with PML is present. Numerical experiments show that our scheme is correct, well stronglyweakly scalable and of good adaptation to complex topography.

Key words: seismic wave equation, numerical simulation, perfectly matched layer (PML), spectral element method (SEM), element stiffness matrices decomposition, parallel computing

中图分类号: