修正辛—准粒子法地震波数值模拟

摘要

基于Hamilton力学,从分子动力学角度建立了针对地震波传播的空间离散的准粒子体系,在该体系中各个粒子仅与上、下、左、右和4个对角粒子发生相互作用,且认为其作用力和相对位移呈近似线性关系,给出了粒子间的相互作用系数,并将该体系变换至Hamilton系统.在时间离散上,通过在二阶辛格式的基础上构建了一个新的修正辛格式,该格式具有三阶时间精度,并在理论上分析了其数值频散和稳定性.修正辛格式所有的系数均为正系数,符合时间演进方向,且具有长时间计算能力.为了测试修正辛格式和空间准粒子体系的准确性,利用拉梅问题验证对弹性波数值模拟的计算精度和效率,为了测试稳定性,选取Sigsbee 2B速度模型进行验证.结果表明,修正辛格式在数值频散压制和数值稳定性提升等方面具有明显的优势.