叠前地震反演已成为获取储层流体敏感参数的重要途径并得到了广泛应用.叠前地震反演方法大体可分为三类：基于波动方程的叠前反演、基于Zoeppritz方程或其近似的AVO/AVA反演和叠前弹性阻抗反演(Mallick，1995).弹性阻抗反演方法简洁、方便、效率高，在抗噪能力和稳定性方面比其 他两种方法具有优势(Cambois，2000)，其实际应用 效果也更好(Karimi et al., 2013; Westeng et al., 2014).

弹性阻抗的概念最初由Connolly(1999)提出，由于其对非零炮检距数据的反演具有较高的稳定性，随后得到了广泛的认可.Whitcombe等(2002)提出了扩展的弹性阻抗的概念，扩大了弹性阻抗的定义域；Verwest等(2000)根据射线传播理论提出了射线阻抗的概念，将弹性阻抗理论扩展到射线域；倪逸(2003)结合相邻层速度和密度的不同变化规律，提出了基于范数动态可调的弹性阻抗计算方法；马劲风(2003)提出了适用于叠后资料的广义弹性阻抗；Wang 等(2006)提出了基于Gray近似方程的弹性波阻抗反演方法，通过反演和参数提取可以直接得到拉梅常数；印兴耀等(2010)基于Russell提出的多孔流体饱和弹性介质的近似方程，提出了包含流体因子f的弹性波阻抗方程及反演方法；Morozov(2010)将反射波、透射波和转换波的弹性阻抗都表示为矩阵形式，无需近似从而提高了计算精度；Zhang等(2012)根据弹性阻抗和射线阻抗的概念提出了新的射线阻抗概念，不仅提高了反射系数的精度而且提高了解释能力；Zong等(2013)推导出包含泊松比和杨氏模量的弹性阻抗方程，建立了相应的直接反演算法.随着地质统计学反演理论的快速发展，一些学者提出利用地质统计学反演方法来反演弹性阻抗，取得了良好的效果(Ruiz et al., 2012；Azevedo et al., 2013; Francis，2013).

目前，广泛使用的弹性阻抗方程大多是由Aki近似方程及其他基于Aki近似方程得到的AVO近似方程经过推导得到的.方程通常是由三项组成的，而且包含密度项.在缺乏大角度地震资料或者地震资料数据信噪比较低的情况下，大角度反演的结果品质较差(李爱山等，2009)，提取参数的过程中，参数矩阵的稳定性较差，造成最终提取的弹性参数误差较大.另外，由于密度项的参数权值在角度有限情况下较小，反演的难度更加大.为此，本文基于入射角AVO近似方程，结合弹性参数间的岩石物理关系，在保证方程精度的前提下推导得到了包含Russell流体项的两项AVO近似方程，并进一步推导出两项流体弹性阻抗方程，建立了相应的概率化叠前反演方法.该方法由于未知项只有两项，对叠前资料的角度要求较小，可以在缺少大角度地震资料的情况下使用，并得到较为合理的反演结果. 2 基于入射角的流体弹性阻抗 2.1 包含流体项的两项近似方程

Aki近似方程是常用的AVO近似方程，在实际应用中为方便使用经常将入射角作为角度参数.由于Aki近似方程是基于入射角和透射角的平均角度进行推导得到的，所以在界面两侧物性变化较大或入射角较大的情况下，该近似方程的计算结果会产生较大的误差.张广智等(2012)提出基于入射角的两项AVO近似方程来解决误差问题，笔者将该理论依据应用于Aki三项AVO近似公式，可以得到如下的近似方程：

式中：v_P是纵波速度 v_S是横波速度；ρ 是密度；是入射角； R(θ_i) 表示纵波反射系数.

基于公式(1)推导得到的弹性阻抗反演方法能够得到比常规弹性阻抗反演方法更加精确的结果(Li et al., 2012)，因此我们在此基础上进行推导.

根据Gardner公式(Gardner et al., 1974)以及纵横波阻抗的表达式可以得到：

式中，

Russell等(2003)对饱含流体多孔介质的AVO理论进行了研究，提出了流体项F(F=fρ).根据纵横波速度、密度与模量的关系可以得到：

与普通弹性阻抗方程推导过程类似，反射系数可以表示为：

由公式(6)可得：

类似于常规的弹性阻抗公式，公式(9)也存在求取的弹性阻抗值随着角度的变化在量纲尺度上有很大变化的问题，这不利于进行不同角度的弹性阻抗值之间的对比以及与波阻抗值的对比.在综合分析波阻抗与弹性阻抗时，首先要将弹性阻抗变换到波阻抗的量纲尺度上，这给实际工作带来了不便.为了解决这个问题，消除入射角变化对量纲尺度的影响，要对推导出的弹性阻抗公式进行标准化处理(Whitcombe，2002).引入了2个参考常数： F₀、μρ₀，并把弹性阻抗方程修改为

如果这些常数值被定义为 F、μρ 曲线的平均值，这样求得的弹性阻抗值就会在单位1附近变化，这一修改使得函数不再依赖尺度.为了使得弹性阻抗的尺度与波阻抗一致，利用θ_i=0°时的弹性阻抗值进行加权处理，即

为了验证新弹性阻抗方程的准确性，使用经典的Ostrander模型进行精确度定量分析，模型数据如表 1所示.

表 1(Table 1)

表 1 Ostrander的三层含气砂岩与页岩模型 Table 1 Model for gas-bearing s and and shale designed by Ostrander VP/(m·s-1) VS/(m·s-1) ρ/(g·cm-3) σ VP/VS 页岩 3048 1244 2.4 0.4 2.45 含烃砂岩 2438 1625 2.14 0.1 1.50 页岩 3048 1244 2.4 0.4 2.45

表 1 Ostrander的三层含气砂岩与页岩模型 Table 1 Model for gas-bearing sand and shale designed by Ostrander

入射介质波阻抗小于透射介质波阻抗的界面称为正波阻抗界面，反之称为负波阻抗界面.上述模型中，上覆页岩和下伏含气砂岩的反射界面就是负波阻抗界面，而上覆砂岩和下伏页岩的反射界面为正波阻抗界面.以上述模型为基础分别用精确的Zoeppritz方程、公式(11)(记为TEI)、Connolly弹性阻抗方程(记为CEI)、基于入射角的弹性阻抗方程(记为IEI)计算了界面处的反射系数，如图 1a，1b所示.

图 1

Fig. 1

图 1 精确Zoeppritz方程(红色实线)、TEI(蓝色点线)、CEI(绿色点线)、IEI(粉色点线)计算的反射系数 (a)负波阻抗界面反射系数对比；(b)正波阻抗界面反射系数对比. Fig. 1 AVO curve computed with the exact Zoeppritz equation(red line)，TEI(blue line)，CEI(green line) and IEI(purple line) (a)Comparison of coefficients at positive interface;(b)Comparison of coefficients at negative interface.

从图 1a中可以看出，在角度较小(小于15°)时，几种公式计算结果接近.随着角度逐渐增大，IEI计算结果与精确Zoeppritz方程的计算结果吻合度比较高，这符合IEI的性质；TEI作为IEI的两项近似也有较高的精确度，在角度增加到40°时与精确值的误差仍然非常小；而CEI在角度大于20°以后误差逐渐增大，精度已经不满足弹性阻抗反演的需求.在图 1b中可以得到类似的结论，因此可以认为新的弹性阻抗方程的精度符合弹性阻抗反演的要求，而且高于常规弹性阻抗方程. 3 基于贝叶斯理论的弹性阻抗反演

基于入射角的弹性阻抗反演需要综合地质、测井和地震的资料，它以包含丰富地下信息的地震反射资料为主要资料，以地质和测井资料作为约束，来揭示地下储层的属性及其含流体特征.主要包括以下步骤：地震资料处理、测井资料处理、子波的提取与标定、弹性阻抗反演.为了更好地利用待反演参数的先验信息和提高反演结果的合理性，本文采用基于贝叶斯理论的弹性阻抗反演方法.

基于贝叶斯理论得到的反射系数 r 的后验概率密度函数可表示为：

假设地震资料背景噪声服从高斯分布，地震数据和反射系数之间的关系可以用服从Gauss分布的似然函数来描述.假设先验信息服从Cauchy分 布，将先验分布函数和似然函数代入贝叶斯公式可得：

将(13)式代入边缘化公式，取对数后得到最大化后验概率分布目标函数为

因此，式(15)即为最终的目标函数，对目标函数求梯度并令其等于0，得到如下反演方程：

从式(16)所示的反演方程可以看出，该方程有一定的弱非线性，因此采用迭代重加权最小二乘算法进行求解即可得到最终的结果.

从反演得到的弹性阻抗数据中提取流体因子等弹性参数是叠前地震反演的重要步骤，在此考虑直接提取Russell流体项F，一方面可避免间接计算带来的累计误差；另一方面两个角度的弹性阻抗建立的方程稳定性更高. 4 模型试算

为了验证基于入射角的流体弹性阻抗反演方法的可行性，本文选取Marmousi2模型(Martin et al., 2006)的部分数据建立二维模型(图 2a，2b)，在时间0.9～1.0 s、CDP300～500处有一含气储层，时间1.45～1.55 s、CDP250～850处有一含油储层.从图 2a中可以很清楚地看到气层和油层，说明Russell流体项对于区分油气层效果非常好，而在图 2b中油气层都不清楚，说明 μρ 本身的指示能力较差.

图 2

Fig. 2

图 2 模型数据 (a)Russell流体项；(b)μρ. Fig. 2 The model (a)Russell fluid term;(b)μρ.

利用图 2数据进行正演得到1°到30°角度地震数据，得到中心角分别为5°、10°、15°的三个角度叠加道集.为了得到公式(11)中的常数b值，从模型中抽取几道作为伪井，对纵波速度和密度取对数进行交会，如图 3所示.经过线性拟合可以得到二者的线性关系表达式，其中斜率参数0.26即为b的取值.利用本文方法对5°和10°数据进行叠前弹性阻抗反演，并提取弹性参数直接得到Russell流体项和μρ，如图 4a，4b所示.分别对比图 2a，2b和图 4a，4b，可 以看出通过本方法估算得到的Russell流体项和 μρ 与模型值基本一致，精确度较高.为了检验方法的抗噪性，在叠前道集数据中加入信噪比为4 : 1的随机噪声，再进行如前的叠加和反演.估算得到的Russell流体项如图 5a所示.对比图 5a和图 2a可以看出，虽然受到噪声的影响，但是估算得到的Russell流体项精确度依然非常高.

图 3

Fig. 3

图 3 纵波速度与密度做对数后的交会 Fig. 3 Crossplotting of the logarithm of P-velocity and density

图 4

Fig. 4

图 4 无噪声情况反演结果 (a)Russell流体项；(b)μρ. Fig. 4 Inversion result without noise (a)Russell fluid term;(b)μρ.

为了体现直接提取Russell流体项的优势，采取基于Connolly弹性阻抗方程的常规弹性阻抗反演方法，利用三个角度道集进行反演得到纵横波速度、密度后间接计算得到Russell流体项，如图 5b所示.对比图 5a和5b可以看出，由于经过了间接计算，反演过程中的误差得到了进一步的累积，因此图 5b的结果受噪声影响更大，Russell流体项预测效果差于图 5a的结果.为了更加直观地比较两种结果的差别，抽取CDP400的同时经过气层和油层的单 道数据进行对比，如图 6.从图中可以看出，直接提取的结果与模型精确值更加接近，间接方法得到的结果误差较大.通过以上比较可以看出，在缺少大角度地震资料的情况下，本文的反演方法依然可以得到准确的反演结果，而且由于系数矩阵条件数较小，稳定性得到了保证.

图 5

Fig. 5

图 5 信噪比为4 : 1时的反演结果 (a)本文方法得到的Russell流体项;(b)间接方法得到的Russell流体项. Fig. 5 Inversion result when S/N=4 : 1 (a)Russell fluid term inverted by the new method;(b)Russell fluid term calculated by indirect method.

图 6

Fig. 6

图 6 信噪比为4:1时的反演结果抽取过CDP400的单道 Fig. 6 CDP 400 of the inversion result when S/N=4 : 1

实际资料来自中国东部某油田，由于目标储层埋藏较深，无法获取高质量的大角度地震数据，因此应用两项弹性阻抗公式进行叠前反演和储层预测.

进行叠前反演之前，需要对地震数据、测井数据进行预处理.地震数据预处理主要是保幅处理，测井数据的预处理主要是井数据的拼接和校正.

通过岩石物理分析，可以得到符合工区地质特点的 γ²_dry、γ²_sat、b 等参数值.弹性参数的统计和交会可以得到工区储层对于弹性参数的敏感性，如图 7是 Russell流体项 F和μρ 的交会图，图中红点表示含油砂岩，绿点表示含水砂岩，蓝点表示泥岩.从图中可以看出，Russell流体项能够很好地区分含油砂岩和含水砂岩，因此可以作为本工区的流体指示参数.

图 7

Fig. 7

图 7 Russell流体项和 μρ 的交会图 Fig. 7 Crossplotting of Russell fluid term and μρ

经过叠前处理，将共反射点道集由偏移距域转换到角度域，并得到中心角为8°(3°~13°)和18°(13°~23°)的部分角度叠加道集数据，如图 8a，8b所示.

图 8

Fig. 8

图 8 不同角度的部分角度叠加道集 (a)小角度叠加道集；(b)中角度叠加道集. Fig. 8 Seismic profiles with different angle range (a)Seismic profile of small angle-stack；(b)Seismic profile of middle angle-stack.

图 8a，8b中CDP405处有一口井，井上的柱状曲线是储层流体解释曲线，红色部分是含油砂岩，绿色部分是含水砂岩，蓝色部分是非储层.利用本文方法进行弹性阻抗反演，得到的Russell流体项和 μρ 如图 9a，9b所示.从图 9a中可以看到，井上所示2.63~2.66 s处的油层对应较低的Russell流体项值，而2.39~2.41 s处和2.48~2.51 s处的水层也对应较低的异常值，但是明显高于含油层的值，储层流体得到了很好的识别.该结果与岩石物理交会分析结果以及实际钻井资料匹配较好.

图 9

Fig. 9

图 9 反演后提取的弹性参数 (a)Russell流体项；(b)μρ. Fig. 9 Elastic parameters extracted with inversion results (a)Russell fluid term;(b)μρ.

上述工区取得较好的应用效果后，在另外一个工区进行实际应用.该工区目的层埋深与前述工区相比较浅，但是由于数据采集条件的限制无法得到高质量的大角度数据，因此使用相同方法进行弹性阻抗反演和储层流体识别.图 10是最终得到的弹性参数剖面.与图 9类似，Russell流体项剖面可以较好地指示含油储层，应用效果较好.

图 10

Fig. 10

图 10 另一工区反演得到的弹性参数 (a)Russell流体项；(b)μρ. Fig. 10 Elastic parameters extracted in another project (a)Russell fluid term；(b)μρ.

本文针对地震数据角度范围较小的情况，基于孔隙弹性介质理论和入射角AVO理论推导出了包含Russell流体项的两项弹性阻抗方程，并建立了基于贝叶斯框架的弹性阻抗反演方法，能够直接提取Russell流体项，不仅使解的稀疏性更高，而且避免了累积误差，提高了计算结果的精度.经过经典模型的试验和实际资料的应用，证明了基于入射角的两项流体阻抗反演方法有较高的准确性和较强的实用性.该方法由于未知项只有两项，对叠前资料的角度要求较小，可以在缺少大角度地震资料的情况下使用，并得到较为合理的反演结果.