弹性波传播的低数值频散波场模拟方法
作者: 张朝元
关键词: 弹性波方程 Runge-Kutta 方法 近似解析离散化算子 通量校正传输技术 数值频散 波场模拟
摘要：压制数值频散以提高计算精度是检验地震波数值模拟方法的一个重要标准。基于弹性波传播方程，建立了低数值频散波场模拟的八阶 FNRK 方法。该方法以 Runge-Kutta 方法对时间导数进行三阶离散，以近似解析离散算子替代差分算子对空间偏导数进行八阶离散，结合通量校正传输技术消除离散后的数值频散。弹性波场模拟结果表明，与高阶有限差分方法相比，该方法能在压制数值频散方面具有明显的优势，计算精度提高，且适应于地震波在大规模复杂介质中传播的波场模拟。