李安昱,尤加春,刘学伟,韩文功,张光德,杨德宽

(1.中国地质大学(北京)地球物理与信息技术学院,北京100083;2.成都理工大学地球物理学院,四川成都610059;3.中国石化石油工程地球物理有限公司,北京100728;4.中国石化石油工程地球物理有限公司胜利分公司,山东东营257086)

近年来,为满足岩性勘探的需要,保幅偏移技术成为地震资料处理中成像技术的热点之一。目前主流保幅偏移方法大致可分为射线类保幅偏移方法、单程波真振幅偏移方法和逆时偏移方法三大类。传统射线类保幅偏移算法对射线积分算子加权实现保幅,该类算法对于复杂构造成像的能力有限[1-2]。单程波真振幅偏移算法对常规单程波算子加入一定的振幅修正项以实现保幅,但单程波本身是对全声波方程的近似,导致横向变速较大时对大角度构造的成像效果较差[3-5]。逆时偏移算法理论上能处理波在任意速度变化介质中的传播问题,但存在低频噪声、多次波虚假成像和计算量大的缺点[6-7]。目前的保幅偏移方法在解决岩性勘探的难题时依然受诸多限制,射线类和单程波保幅方法虽然一定程度上提高了振幅保真能力,但该类方法解决复杂构造保幅成像的能力偏弱,工业应用的性价比较低。逆时偏移对计算机硬件性能要求较高,尽管发展了多种加速方法,仍然不能广泛应用于工业生产[8-9]。从偏微分方程数值求解的角度分析,声波方程是一个关于时间和空间变量的二阶偏导数微分方程,理论上在深度域或时间域求解需要两个边界条件方可准确求解。……