应用球型差分模板的低秩有限差分法纯qP波逆时偏移

摘要

各向异性介质拟声波方程正演模拟及逆时偏移存在伪横波干扰及数值不稳定现象,成像质量不佳;同时计算效率也是各向异性逆时偏移的重要影响因素。为此,提出一种兼具成像质量和计算效率的纯qP波逆时偏移成像方法。首先,从波动方程伪解析解和各向异性精确频散关系出发,分别构建了VTI和TTI介质纯qP波伪解析解波场延拓公式,避免了推导显式波动方程,无需对频散关系做平方处理,由此消除了波场模拟存在的伪横波干扰;然后,通过设计三维球型或二维圆形差分模板,并借助低秩有限差分法求解与模型相适应的差分系数,推导了基于低秩有限差分的纯qP波波场延拓公式,在时间方向延拓时,无需进行多次Fourier变换,降低了计算复杂度;最后,在此基础上,利用GPU并行计算正、反向地震波场,提高了逆时偏移的成像效率。典型二维和三维模型试算结果表明:该方法能够消除伪横波干扰,保证计算过程稳定,在二维和三维各向异性介质中都具有较强的适应性和较高的计算效率。