高效高保真地模拟地震波传播对深入研究地球内部精细结构十分重要。然而现有的地震波模拟算法存在着不同的瑕疵在非均匀分界面精度不足、粗网格条件下会产生数值频散、长时计算后会发生振幅改变、相位漂移等现象。因此发展具有高精度、低数值频散且能保持真振幅的数值模拟方法不仅是地震学研究领域的基础性问题之一,也是该领域亟待解决的重要问题之一。本项目拟通过研究复杂介质中的地震波传播方程,结合模拟地震波场的近似解析离散化方法、模拟大气海洋环流的平方守恒格式、并匹配以合适的保辛算法以实现时间层推进计算,得到具有高精度、低数值频散,长时计算保持真振幅的地震波场模拟算法,并研究算法在不同集群架构上的并行算法及其并行性质。本项目的完成将对研究地球内部精细结构、理解地震波观测资料细节及地震波反演产生重要的影响;同时也将为研制和使用固体地球物理专用的并行集群提供经验支持,有利于从软硬件两个方面推进地球物理学研究的发展。
英文主题词seismic wave field simulation;MNAD;WONAD;NSPRK;parallel algorithm