周东红　李景叶　陈　莉

(①中海石油(中国)有限公司天津分公司渤海石油研究院，天津塘沽 300459； ②中国石油大学(北京)，北京 102249)

1　引言

地震反演在油气勘探与开发中发挥着十分重要的作用[1-4]。叠后地震反演方法是地震勘探界应用较广泛的反演方法之一，在实际勘探开发中也取得了较好的效果。叠后反演方法可分为递推反演、模型反演、地震属性反演、地震统计学反演和测井曲线反演等[5]。与叠前反演相比，叠后反演具有更高的计算效率和更强的抗噪性，一直以来都被认为是储层识别与表征的核心技术。然而目前的叠后反演大多基于褶积模型，在各个界面单独计算反射系数，而忽略了波的传播效应以及地层厚度的影响，并且一般只能求取地下介质的弹性阻抗或者是声阻抗，难以直接获取储层的弹性参数，如纵波速度等。

众所周知，地震波在地下介质中的传播会受到传播效应的影响，例如层间多次波、表层多次波、透射损失和转换波等[6]，且在入射角大于30°时这种影响尤为明显。弹性波动方程能够进行全波场模拟，基于弹性波动方程的反演精度更高，但受计算效率、反演稳定性和地震资料品质等因素的制约，目前弹性波动方程反演还难以在实际生产中推广应用。

反射率法(The Reflectivity Method)是一种向量化的递推计算方法，是基于一维层状地球模型假设下进行全波场模拟的方法，该方法模拟的波场包括反射波、透射波、转换波和多次波，考虑了地震波相位变化、地层厚度以及透射损失的影响，计算精度高于有限差分法和射线法，而计算成本介于两者之间[7,8]。……