地震波数值模拟是地震勘探和地震学研究的重要部分。尽管现有的如有限差分、有限元等技术基本能满足小尺度问题的模拟需要,然而对于大尺度,尤其是全球尺度地震波传播的计算而言,由于涉及到长时程追踪,且前述方法大多是非高精度保结构方法,长时程计算时不可避免地产生高积累误差,导致数值频散严重,很难精确处理全球尺度地震波传播这一长时程追踪问题。因此,发展针对于全球尺度地震波传播数值模拟的高精度保结构算法成为关键。本项目针对球坐标下非均匀介质中弹性波方程拟构造谱元法结合新推的力梯度辛算法的离散方法,通过对其保结构性、精度及效率等问题的研究,发展适用于全球尺度地震波传播的高精度保结构数值模拟方法,使长时程模拟时的积累误差大幅降低,数值频散得以压制;并用间断算子处理地球内部的非均匀层间断面,显著提高模拟精度。保结构新方法的推出,将为高精度处理类似全球尺度三维非均匀介质中地震波传播的长时程计算问题提供新的途径。
英文主题词seismic wave propagation;structure-preserving algorithm;spectral element method;Earth free oscillations;