波动方程有限差分数值模拟是研究地震波在地下介质中的波场特征和传播机理的重要手段。对于常规有限差分技术,当采用大网格对计算空间进行差分离散时会出现严重的数值频散问题,降低了计算精度。通量校正(Flux-corrected transport method,FCT)技术能够有效压制粗网格情况下有限差分的数值频散。本文研究了具有垂直对称轴横向各向同性(Vertical Transverse Isotropy,VTI)介质的交错网格FCT有限差分技术。首先从一阶速度—应力弹性波方程出发,在交错网格空间中给出了该方程的高阶有限差分法格式及稳定性条件,在此基础上研究了波动方程正演过程中的数值频散FCT压制技术,二者结合实现了该方程的高精度有限差分数值模拟。同常规算法相比,本文算法不额外增加内存需求,少量增加计算量,但可有效压制VTI介质中弹性波动方程正演的数值频散现象。当采用大网格进行数值模拟时,本文方法明显提高了波场模拟精度。
