韩磊，刘俊州，杨睿，张广智，周游
（1.中国石化石油勘探开发研究院，北京 100083；
2.中国石油大学（华东）地球科学与技术学院，山东青岛 266580）
当前，随着非常规裂缝性油气藏勘探开发的深入进行，利用地震反演方法预测储层裂缝成为一种有效的技术手段。对于致密砂岩储层，地下裂缝系统是重要的油气储集空间和流体渗流通道，大大提高储层有效孔隙度和整体渗透率，对油气开采和油气田产量有着重要影响。因此，分析地震波在致密砂岩储层中传播时的地震响应特征，基于方位叠前地震数据开展弹性阻抗反演方法，准确估算储层弹性参数和表征裂缝发育强度的各向异性参数，进而预测储层位置，为致密砂岩储层的勘探和开发提供可靠的地球物理依据[1−3]。
地震各向异性理论是叠前地震反演预测裂缝的基础[4]。地震各向异性的强弱能够有效表征地下裂缝的发育程度，进而预测储层的位置，最终指导致密砂岩储层的勘探和开发。THOMSEN[5]通过研究地震波在各向异性介质中的传播，基于弱各向异性假设构建了表征介质各向异性强弱的3 个参数。RÜGER[6]推导了多种各向异性介质的反射系数近似公式，为裂缝性储层叠前地震反演奠定了理论基础。针对储集空间复杂且埋藏较深的致密砂岩油气藏，地震资料信噪比很低，限制了叠前AVO反演的精度。CONNOLLY[7]首次提出弹性阻抗反演方法，相比以往的叠前AVO反演方法，具有更好的稳定性和鲁棒性；
WHITCOMBE 等[8−9]在其基础上提出了扩展弹性阻抗（EEI）反演，用于储层的岩性检测和流体识别；
WANG 等[10]推导了以拉梅参数表示的弹性阻抗方程，并在实际工区中进行了应用。陈天胜等[11]推导了简化的HTI 介质弹性阻抗近似公式。陈怀震等[12]利用基于方位弹性阻抗反演方法提取了储层弹性参数和各向异性强度参数。
对于低孔低渗的致密砂岩储层，研究致密砂岩中裂缝的发育程度对储层的预测与评估具有重要意义。国内外许多专家学者在描述裂缝性储层中地震响应与裂缝特征的关系上做了许多探索。PEREZ等[13]基于方位AVO 技术利用P 波地震数据实现实际工区的裂缝检测。MALLICK 等[14]首先利用海洋地震资料的方位纵波AVO 技术预测垂直裂缝。GRAY等[15]利用方位3D 地震数据，开展叠前AVAZ 反演方法研究，预测实际工区的裂缝特征（裂缝方位和密度）。印兴耀等[16]提出一种无约束的稳定方位AVO梯度反演方法，用来预测储层的裂缝发育情况。宗兆云等[17]充分利用叠前地震资料，开展基于贝叶斯理论的弹性阻抗反演方法，实现了碳酸盐岩和碎屑岩的裂缝含气性检测。撒利明等[18]详细讨论了缝洞型储层的地震响应特征与识别方法，并通过模型和实际数据分析实现了地下溶洞特征描述。对于页岩储层预测及特征描述的研究，张广智等[19]基于叠前纵横波联合反演方法开展了页岩气储层的识别及流体预测的研究。陈怀震[20]基于岩石物理模型，针对裂缝性储层实际特点，开展了弹性参数和裂缝特征的地震反演预测方法研究。对于致密砂岩裂缝性储层的预测方法研究，国内罗群、马如辉、黄绪德等做了一些工作[21−23]，但基于叠前地震反演方法实现致密砂岩储层裂缝特征描述的相关文献资料较少。
主要针对川西坳陷新场地区须家河组致密砂岩裂缝性储层，针对工区储层特征及预测中存在的问题，从等效VTI 介质理论出发，利用RÜGER 推导的VTI介质近似公式，开展基于各向异性弹性阻抗反演的裂缝参数提取方法研究。利用工区的实际数据进行反演方法的试算，验证了本方法的可行性。
1.1 工区概况
须家河组发育三角洲平原—前缘沉积砂体，主要砂岩储层的发育层段为须二段和须四段。储层类型复杂，分为4 种类型：裂缝型，裂缝−孔隙型、孔隙−裂缝型、孔隙型。新场须二构造是由5个局部高点组成的NEE 向复式背斜，构造南陡北缓。断裂走向复杂，主要为近SN、NE、EW 向，规模较大的断层主要为SN 向。纵向上划分为3 个亚段，其中主产层为上亚段T2x2和T2x3，中亚段T2x4、T2x5以及下亚段T2x7。纵向物性变化：新场构造带纵向上T2x4、T2x5砂体物性最好，其次为T2x2、T2x1砂体，均属低孔低渗致密储层。须家河组致密砂岩裂缝性储层有机质含量丰富，但地质条件复杂、储层物性条件差、非均质性强、流体关系复杂。目前仍存在几个关键问题亟待解决：
1）缺少表征裂缝发育特点的各向异性参数曲线。致密碎屑岩储层地质特点复杂。裂缝往往与优质储层相关联，裂缝发育导致各向异性特征强烈，需要根据储层特点，建立合适的岩石物理模型，进而估测各向异性参数曲线，确定优质储层位置。
2）缺少针对性的各向异性参数反演方法。须家河组气层纵横向分布非均质性极强，主要是受复杂的地质条件影响，导致了常规地震反演储层识别难度很大，需要针对储层发育特征，选择合适的各向异性反演方法，提高各向异性叠前地震反演的稳定性。
1.2 针对性技术手段
针对川西新场工区地震道集的特点,形成了一套工区地震资料的处理流程，处理后信噪比与分辨率明显提升，为开展储层裂缝预测提供了数据基础。
综合考虑工区岩石矿物、裂缝特征，根据等效介质理论构建了低角度裂缝发育的VTI介质等效模型，实现了裂缝参数的定量表征，为后期开展储层各向异性反演提供井中先验信息。
常规波阻抗反演对致密砂岩储层中流体和特殊岩性的预测不够敏感，并且角道集资料的信噪比又是制约其广泛应用的重要因素。基于叠前地震资料的弹性阻抗反演是部分角度叠加数据反演的结果，通常远角度弹性阻抗比近角度阻抗对流体含量以及裂缝储层更加敏感，可以更好地适应生产；
针对每个角道集数据提取子波，消除通常出现在远炮检距道集中的低频转移现象，从而更加稳定地反演储层参数以及裂缝参数。
须家河组气层纵横向分布非均质性极强的特点，导致常规地震反演储层识别难度加大。针对这一问题，根据储层发育特征，研究形成了一套合适的各向异性反演方法流程，提高了各向异性叠前地震反演的稳定性。同时，完成了川西新场工区地震资料各向异性参数反演，裂缝预测效果与钻井产能结果较为吻合，为开展裂缝分布规律研究提供了地球物理支持。
2.1 等效VTI介质模型
沉积地层中最常见的各向异性模型就是横向各向同性介质（TI）。须家河组致密砂岩储层低角度缝、平缝普遍发育，岩心照片如图1所示。根据等效介质理论可以将其等效为具有垂直裂缝的横向各向同性（VTI）介质。VTI介质主要用于描述形成泥页岩过程中各种构造运动及外界作用下形成的成层状和片状的岩石特征，图2为VTI介质模型的示意图。
图1 岩心照片Fig.1 Core sample
图2 VTI介质模型Fig.2 VTI medium model
2.2 VTI介质中的弹性阻抗方程
根据RUGER 推导的反射透射近似公式，利用THOMSEN 各向异性参数来确定各向异性强弱，用以衡量地震波的反射透射系数的变化规律，并以各向同性介质为背景，推导出了基于各向异性介质反射透射系数的线性化近似表达[6]：
RÜGER 将上述关系式分解为各向同性部分和各向异性扰动部分：
其中，
选取能够适用于纵波阻抗和横波阻抗参数反演的裂缝性储层反射系数各向同性部分（）Fatti 近似公式[32]：
其中，a(θ)=1+tan2θ，b(θ)=−8gsin2θ，c(θ)=4gsin2θ−tan2θ，d(θ)=sin2θ，e(θ)=sin2θtan2θ。
式（1）—式（5）中：θ为入射角，°；
Vp为纵波速度，km/s；
Vs为横波速度，km/s；
ρ为密度，kg/m3；
IP为纵波阻抗，kg·m−2·s−1；
IS为横波阻抗，kg·m−2·s−1；
δ、ε为无量纲的Thomsen各向异性参数。
式中：EI为弹性阻抗，其随入射角的变化而变化；
K为与纵横波速度比的平方相关的常量。
对上述RPP反射公式两边同时积分再取对数可得到VTI介质中的弹性阻抗表达式：
方程（7）是开展VTI 介质弹性阻抗反演的基础，这里通过反演的不同入射角（小、中、大）对应的弹性阻抗，进一步提取储层弹性参数和THOMSEN 各向异性参数，进而有效指导目标区域的成层状地质特征预测。
2.3 基于各向异性弹性阻抗的裂缝参数提取流程
根据褶积公式，结合已知的井数据资料，合成不同入射角、不同的地震记录。
式中：S(θ)为不同的地震记录；
W(θ)为不同的地震子波（此处采用雷克子波）。
当确定弹性阻抗后，提取裂缝性储层弹性参数（IP、IS）和各向异性参数（Δδ、Δε）的方法如下：
从前面推导出的各向异性弹性阻抗公式可知，弹性参数（IP、IS）和各向异性参数（δ、ε）与弹性阻抗之间是非线性指数关系，直接求解比较复杂。因此，需对各向异性弹性阻抗公式两侧同时取对数，将其线性化：
从以上公式可以发现，将具有n个入射角的地震道集代入求解，则不同入射角的弹性阻抗的对数方程组为：
使用求解混定问题所使用的马奎特反演方法（又称阻尼最小二乘法）求解式（10）中的纵横波阻抗以及裂缝各向异性参数。
综合上述讨论分析，提出了利用各向异性弹性阻抗进而提取弹性参数和各向异性参数的流程，具体流程如下：
1）对已经确定的叠前地震道集做角度叠加，得到不同角度的道集数据；
2）提取相应的不同角度的地震子波；
3）将岩石物理模型和测井数据作为初始模型约束；
4）对不同角度、不同方位的地震数据进行稀疏脉冲反演得到弹性阻抗数据体，并取对数；
5）基于获取的弹性阻抗数据体和初始模型约束提取弹性参数和各向异性参数。
选取某裂缝性储层工区的C 井建模，如图3所示，利用该模型制作合成地震记录并添加不同信噪比的高斯噪声，采用上面提出的方法流程进行基于方位弹性阻抗反演的致密砂岩储层弹性参数和各向异性参数提取方法试算，如图4、图5所示。反演结果表明，在无噪声情况下，该方法能够获取与真实值基本吻合的纵横波阻抗及裂缝各向异性参数，验证了方法的可行性。在适当噪声的情况下，纵横波阻抗及裂缝各向异性参数的反演结果与真实值仍有较高的吻合度。
图3 三层VTI介质模型Fig.3 Three-layer VTI medium model
图4 信噪比为5时测井数据模型参数反演结果与真实值对比Fig.4 Comparison between estimated results and real data when S/N is 5
图5 信噪比为1时测井数据模型参数反演结果与真实值对比Fig.5 Comparison between estimated results and real data when S/N is 1
选取的实际工区位于四川盆地新场工区须家河组，本区目标储层是典型的致密砂岩气藏。工区主要发育高、低角度裂缝，储集类型为孔隙−裂缝型，以裂缝性储层为主。首先对地震资料进行品质分析，工区的部分CMP 道集及其频谱分析如图6、图7所示，可以看出，叠前CMP 道集信噪比较低，随机高频噪音严重，能量分布不均匀。需要做进一步处理，提高地震资料品质。经过处理后的叠前CMP 道集如图8所示。利用处理后的叠前CMP 道集进行该实际工区反演试算，共有11口井参与模型约束，提取储层弹性参数以及各向异性参数，如图9、图10所示。实际工区反演结果表明，通过地震数据反演得到的岩石弹性参数IP、IS与井上信息匹配效果良好，井符合率在80%以上，能够很好地指示砂体位置。反演估算得到的裂缝各向异性参数Δδ、Δε呈现出明显的高值特征，能够很好地指示砂体裂缝发育区域，与测井岩心照片匹配一致。
图6 原始叠前CMP道集Fig.6 Original pre-stack CMP seismic profiles
图7 叠前CMP道集频谱Fig.7 Spectrum of pre-stack CMP seismic profiles
图8 处理后叠前CMP道集Fig.8 Processed pre-stack CMP seismic profiles
图9 实际工区弹性参数IP、IS反演结果Fig.9 Estimated inversion results of elastic parameters IP and IS in real workflow
图10 实际工区裂缝岩石物理参数Δδ、Δε反演结果Fig.10 Estimated inversion results of fracture petrophysical parameters Δδ and Δε in real workflow
针对非常规致密砂岩储层，分析裂缝诱导的储层各向异性特征，推导相应的反射系数近似公式，进一步推导各向异性弹性阻抗公式，建立基于各向异性弹性阻抗的致密砂岩储层各向异性反演方法流程，并最终实现一维井数据和实际工区的储层弹性参数以及各向异性参数的提取。一维模型测试结果表明，该方法具有较好的抗噪性和可靠性；
实际工区的反演结果表明，相比于常规波阻抗反演方法对流体和裂缝性储层预测并不敏感的问题，该方法可较为准确地估算裂缝各向异性参数，进而指导储层中的砂体位置以及裂缝发育区域，对裂缝性储层的勘探和开发具有较好的指导作用。