摘要
波动方程有限差分法是一种使用广泛的地震波数值模拟方法。但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率。为了克服常规有限差分算子的数值频散,该文针对VTI介质地震波数值模拟问题,构造了频率‐空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯‐牛顿法确定了高精度有限差分算子的优化系数。利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明:有限差分优化算子具有较高的渡场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率‐空间域qP波正演模拟奠定了基础。