| 基于多辛结构谱元法的保结构地震波场模拟 |
| |
| 引用本文： | 汪文帅,李小凡,鲁明文,张美根.基于多辛结构谱元法的保结构地震波场模拟[J].地球物理学报,2012,55(10):3427-3439. |
| |
| 作者姓名： | 汪文帅 李小凡 鲁明文 张美根 |
| |
| 作者单位： | 1. 中国科学院地球深部重点实验室, 中国科学院地质与地球物理研究所, 北京 100029;2. 宁夏大学数学计算机学院, 银川 750021 |
| |
| 基金项目： | 国家自然科学基金项目，国家重点基础研究发展计划“973”计划 |
| |
| 摘 要： | 近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一. 本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法. 本文给出的多组数值试验对比结果表明,本算法无论在内存消耗、稳定性及计算耗时,还是长时程跟踪能力方面都有上佳的表现; 另外,本文给出的起伏地表多层介质模型的数值算例验证了该算法处理复杂几何形状和复杂介质时的有效性. 该多辛结构谱元法的发展将为长时程地震波传播的计算及模拟提供更为广泛而有效的选择.
|
| 关 键 词： | 谱元法 辛算法 地震波模拟 |
| 收稿时间： | 2012-03-19 |
Structure-preserving modeling for seismic wavefields based upon a multisymplectic spectral element method |
| |
| WANG Wen-Shuai
,
LI Xiao-Fan
,
LU Ming-Wen
,
ZHANG Mei-Gen.Structure-preserving modeling for seismic wavefields based upon a multisymplectic spectral element method[J].Chinese Journal of Geophysics,2012,55(10):3427-3439. |
| |
| Authors: | WANG Wen-Shuai LI Xiao-Fan LU Ming-Wen ZHANG Mei-Gen |
| |
| Institution: | 1. Key Laboratory of the Earth's Deep Interior, CAS, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;2. School of Mathematics and Computer Science, Ningxia University, Yinchuan 750021, China |
| |
| Abstract: | In recent years, structure-preserving algorithms effective for treating problems of high-precision and long-time tracing have emerged as one attractive topic of the numerical simulation of seismic waves. In this paper, a new kind of time-space structure-preserving algorithm, which combines the spectral element method (SEM) in spatial discretization with a new three-stage third-order symplectical algorithm (NTSTO) in temporal discretization, is designed to simulate seismic wave propagation. The comparisons of numerical experiments show that the NTSTO-SEM algorithm is much superior to the other methods in terms of efficiency, long-time tracing ability, stability, and low storage requirement. Finally, the numerical experiments of a multilayered model with irregular topography are given to illustrate the effectiveness of NTSTO-SEM for handling complex geometry and complex media. These appealing features of the multisymplectic algorithm would make it effective to model the large-scale and long-time seismic wave propagation. |
| |
| Keywords: | Spectral element method Symplectical algorithm Seismic modeling |
| 本文献已被 万方数据 等数据库收录！ |
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
