地震偏移的最优可分近似算法实现

针对地震偏移算法中单程波算子的特征函数近似,采用最优可分表示法将该特征函数展开为空间变量(g)和水平波数(k)的可分表达式,以此可分近似表达式为基础,运用正反Fourier变换重新构造单程波算子.为了克服特征函数最优可分近似计算在奇点及其领域产生的数值振荡,引入等价黏性技巧以增强算法的数值计算稳定性,并采用分频最优可分及对空间变量(g)线性插值的方法,不仅提高了计算精度,也节约计算机时.文中具体研究了单程波算子的脉冲响应和二维叠前深度偏移.结果表明,在不甚大步长情况下,本文构造的算子具有适应横向强变速的能力,运用Marmousi模型验证了本文方法适合复杂构造的成像.

单程波最优化可分近似偏移方法精度分析(英文)

单程波算子的可分近似形式是关于空间变量和波数变量的函数,这使得快速傅里叶变换的利用成为可能,从而极大提高了计算效率。从函数近似的角度来看,最优可分近似方法与其他几种可分近似方法一样都具有单程波算子可分表达的形式,但最优可分近似方法是其中唯一的函数整体近似的方法,这就造成了该方法在相位误差曲线、脉冲响应、模型偏移结果上表现出与其他可分近似方法不一样的特征:该方法精度较高,且随着阶数的增大,对速度变化的敏感度降低。