基于逆时间频散变换和三阶辛积分的长时程数值模拟

摘要

地震波场模拟是全波形反演的重要基础.基于有限差分的直接解法是地震勘探领域应用最为广泛的数值模拟方法.为了降低三维大规模问题的计算量和内存需求,往往需要采用较大的空间步长和时间步长.然而,大步长通常会导致有限差分方法(空间变量和时间变量)在离散过程中会出现数值频散,严重降低长时程波场的精度.严重的数值频散会导致全波形反演失败,尤其是对深层目标危害很大.本文针对时间频散的压制展开了研究,希望能够采用接近于稳定性上限的时间步长以大幅度提升长时程数值模拟的计算效率.