简介：频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播，但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性，本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验，验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理，基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像，说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力，我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理，对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力，具有一定实际应用价值。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。