碳酸盐岩储层地震相控非线性反演技术及应用

高 君，黄捍东，季 敏，苑书金

（1.中国石油化工股份有限公司石油勘探开发研究院，北京100083；2.中国石油大学（北京）地球物理学院，北京102249）

碳酸盐岩的矿物组成、结构以及成岩演化过程非常复杂[1]。碳酸盐岩储集空间类型及其孔缝尺度的多样性导致该类储层具有强烈的非均质特征，其不同类型岩相的波阻抗等地球物理参数呈现交叉叠置的现象。因此直接进行目标层系的地震波阻抗反演难以有效预测和描述储层特征[2]。长期以来，碳酸盐岩储层的地震精细预测一直困扰着地球物理工作者。

碳酸盐岩的发育及其储层特征总体上受沉积环境、成岩作用及其改造作用的影响，尤其是不同沉积相带的储集性能和岩石物理特征均存在差异[3-4]，因此，刻画沉积微相及其展布对碳酸盐岩储层的预测具有重要的指示意义[5-6]。为进一步提高复杂储层的地震预测精度，充分考虑沉积相带的影响，在地震相模式约束下将非线性随机模拟技术与地震反演相结合，使地震反演储层预测结果更符合实际情况。HUANG 等[7]、KORNEEV 等[8]、ZHANG 等[9]和ZHOU 等[10]基于非线性理论及统计学方法，探索并提出了高精度非线性混沌反演算法，有效降低了地震反演的多解性，从本质上提高了地震预测的效果。慎国强等[11]通过叠前叠后的联合反演进一步提高了地震预测精度。黄捍东等[12]、张志伟等[13]、郭恺等[14]、张显文等[15]、杨子川等[16]和李斌等[17]针对塔河、阿姆河等不同地区的碳酸盐岩储层，在地震地质特征分析的基础上，建立岩相与地震反射特征的关系，而后采用相带约束的方法在等时格架的基础上开展地震属性特征分析及相控反演研究，改进了地震反演方法，从而进一步提高了储层预测的精度。关达等[18]针对大牛地气田D 井区储层非均质性强的特点，采用相控非线性地震反演技术，一定程度上克服了运算结果的多解性问题，该技术在华北油田复杂储层地震预测中发挥了重要作用，取得了良好的应用效果。

本文研究区的盐下碳酸盐岩储层埋深大、非均质性强，储层与非储层波阻抗值域范围严重重叠。为提高复杂碳酸盐岩储层的地震预测精度，首先在岩石物理分析及地震相特征研究的基础上，根据实际地震资料建立包含目标层段的宏观相控模型，然后在测井资料的约束下采用非线性随机算法逐道、逐次地进行外推反演，再进一步对反演数据体进行地震精细解释，最终结合神经网络技术实现储层的定量预测。

1 方法原理

针对复杂的碳酸盐岩储集体，采用地震与沉积相相结合的技术思路，将地质研究成果与地震反演有机结合，使地震反演结果最大限度地排除边界效应，减少储层地震预测的多解性。地震相控逐道外推反演方法是从井旁道出发进行反演，在精确的层序地层格架控制下实现反演，并将每次得到的反演结果作为下一道反演的初始模型，以此类推。

由于地震道与波阻抗的关系是非线性的，所以基于地震道的逐道外推反演也可称为非线性随机反演。实际操作时可依据井的波阻抗和地震相划分结果建立地下波阻抗模型，然后利用井点波阻抗和井旁地震道求出控制参数，即可进行井约束、相控非线性随机反演。在实际资料反演过程中，该反演方法可以有效地提高地震资料的垂向分辨率，并充分考虑目标地层层序的地质特征，使反演结果更符合实际情况。由此，基于地震道非线性最优化反演思想，假设岩石密度为常数，将波阻抗反演转化为速度反演，建立地震道与波阻抗关系的目标函数为：

式中：V为速度；SΔi为模型响应，即速度预测结果对应的合成地震记录；Di为实际地震记录，i为地震记录的采样点序号，范围为[0，n－1]。

利用Taylor 公式将（SΔi－Di）在初始模型响应Si处展开，为方便求解将二次项以上的高次项省略，可得：

式中：Si代表速度初始模型对应的合成地震记录；ΔVk，ΔVj代表不同模型参数的摄动量；k，j分别代表模型数和参变量。

对（2）式求导，可得：

对（1）式的ΔV求导并展开，可得：

图1 L油田盐下地震反射剖面及地层层序划分

图2 BVE100层段岩石物理多参数交会分析结果

表1 BVE100层段不同类型储层的弹性参数统计结果

2.2.3 相控地震反演

研究区三维地震资料面积240 km2，选定其中的7口井参与分析和运算，确定地震反演的3个目标层段，即4 个 约 束 层 段 界 面（BVE100、BVE200、BVE300及ITP层段顶界面），以地层层序界面和前期获得的地震相模型为约束展开地震相控非线性反演。钻井揭示，BVE100层段发育叠层石灰岩及颗粒灰岩组合，与上覆厚层盐岩相比，具有明显的低波阻抗特征；BVE200层段是大套的泥灰岩夹薄层灰岩，与上、下围岩的地震反射界面清晰；BVE300层段与下伏地层存在明显的角度不整合特征，主要发育大套泥灰岩组合，其波阻抗值通常较大。相控地震反演过程如下：首先通过井震综合标定，建立基于层序特征的地震相控低频模型，并将其作为反演的基础条件；然后以模型参数作为区域化变量，采用上述非线性逐道外推方法展开随机模拟；再调整变差函数以确立模型参数点之间的统计相关关系，基于目标函数最小化原则将模型参数变量作为符合高斯分布的随机变量；最后通过高精度随机模拟实现约束反演结果与地震道的最佳匹配，进而获得研究区的地震反演速度、密度等参数。图3为地震相控非线性反演的波阻抗与反演得到的地震剖面的叠合显示（红线分别为4个约束层段顶界面）。反演得到的波阻抗不仅与地震波组的振幅、频率、相位的细节变化高度一致，而且反演得到的地震剖面中地震反射的中、弱信号更加清晰，目标层内幕细节更加突出和丰富。

图3 地震相控非线性反演的波阻抗与反演得到的地震剖面叠合显示

2.2.4 应用效果

2.2.4.1 岩性岩相和微生物礁特征

图4为过8-LL-9井与过1-R-682A 井的地震反演波阻抗剖面与井点处目标层段岩性岩相特征，井震特征具有良好的一致性。图4a中的灰岩储层和致密的泥灰岩、图4b顶部的高速石膏和下部较高速度的白云岩波阻抗均存在明显差异，波阻抗反演结果进一步提高了目标层段的纵、横向分辨率。由图5a可知，经地震相控非线性反演后，微生物礁体显示得更清晰，但受地震分辨率的限制，从图5b中无法辨识微生物礁体的细节，只能看到一些较弱的内部反射特征。图5a不仅更清晰地反映了微生物礁体的叠瓦状轮廓，也进一步揭示了其内幕结构特征。

2.2.4.2 储层厚度分布预测

基于地震相控非线性反演结果，利用纵波速度、波阻抗等参数与井点储层厚度进行多参数拟合分析，最终预测出研究区储层厚度的分布情况，进而编制出目标层段一类和二类储层以及三类储层的预测厚度平面展布（图6）。从图6a可以看出，BVE100层段一类和二类优质储层集中发育在两个构造隆起区，最厚的储层发育在南部，厚度达70 m，其次是北部的东西相邻井区，储层最大厚度分别为50 m 和25 m，平均厚度约为20 m。从图6b可以看出，BVE100层段局部发育三类储层，厚度较小，主要分布于研究区的中部，呈条带状南北向局限分布，储层平均厚度约10 m。

对研究区14口井依次统计3个层段的地震反演预测厚度与测井解释厚度吻合情况，将二者相对误差大于20%视为不吻合，统计结果如表2所示。储层厚度预测的总吻合率为90.9%，表明地震相控非线性反演技术能够有效实现复杂碳酸盐岩储层的精细预测。

图4 过8-LL-9井（a）与过1-R-682A 井（b）地震反演波阻抗剖面与井点处目标层段岩性岩相特征

图5 地震反演得到的波阻抗剖面（a）与对应的地震剖面（b）

图6 BVE100层段储层预测厚度平面展布

表2 一类和二类优质储层预测厚度及误差统计结果

3 结论认识

1）相控地震反演技术的核心思想是将地震相时空分布模型与地震反演运算有机结合，使地震反演最大限度地排除边界效应。理想的地震反演需要相控模型的良好匹配和复杂运算的多次迭代。其中的基于层序特征的地震相模型控制了反演的低频趋势，为地震反演提供了必要条件。

2）相控低频模型是实现复杂碳酸盐岩储层精细预测的重要基础，基于相界面和测井资料建立的低频速度模型理想地控制了岩性岩相的宏观变化，从根本上减少了地震反演的多解性。

3）实际资料应用统计结果表明，后验井符合率大于90%，储层厚度与测井资料的相对误差小于20%。该技术具有运算效率高、反演结果随机性小的特点，可以更好地体现相控的思想，并有效实现复杂储层的定量预测。