G 147井区新近系为河流相沉积, 岩性是砂岩与泥岩的互层组合, 储层发育, 分布广, 厚度大。其中馆陶组属辫状河沉积, 纵向上形成“ 砂包泥” 沉积剖面, 碎屑层尤为发育, 一般大于剖面厚度的40%, 砂岩单层厚度平均大于20 m; 明化镇组属于曲流河沉积, 纵向上形成“ 泥包砂” 沉积剖面, 碎屑层一般占剖面厚度的30%左右, 单砂层厚度平均10 m左右, 河道砂体发育, 为油气富集提供了良好的储集空间。

G 147井区明化镇组受控于X 1-13西断层和G 196东断层, 构造幅度达50 m, 埋深1 030~1 080 m。圈闭含油面积0.46 km²。G 147井区南部靠近港西断层, 北部为港西北坡的一部分, 构造整体上为南高北低。该井区内部被南北向断层所切割, 分为X 34-14井区和X 36-12井区。其中X 34-14井区NmⅡ -8小层顶界高点埋深1 040 m, 闭合幅度20 m, 北部受控于岩性下倾尖灭, 圈闭面积0.15 km², 圈闭落实可靠。G 147井区NmⅡ -8小层顶界高点埋深1 030 m, 闭合幅度50 m, 圈闭面积0.31 km², 圈闭落实可靠（图1）。

图1

	Figure Option ViewDownloadNew Window
	图1 G 147井区NmⅡ -8小层含油面积图

G 147井区油藏类型为构造-岩性油藏, 曲流河是河流成熟期发育的地貌特征, 以侧向侵蚀和加积作用为主, 其成因类型多、砂体横向变化快及空间展布复杂^[1]。该区块构造幅度为50 m, 产状较陡, 邻井之间的高低关系变化大, 造成地层对比难。

该井区NmⅡ -8小层有效厚度为1.5~5 m, 而地震波长为30~35 m, 油层薄造成波阻抗反射界面系数较小, 地震响应特征表现为弱反射结构。常用的地震单一属性只对某些地质特征敏感, 利用单一属性来预测有效储层在精度上存在着不足; RGB属性融合技术对图像进行融合, 更全面、准确地说明了目标信息, 但也相对增加了一些储层干扰信息, 在“ 好” 属性中容易出现假储层。两者在精度以及多解性的问题上都存在不足, 给储层预测带来多解性和局限性, 已不适用于该井区。

精细层位标定是把地震波同相轴赋予具体地质意义的关键方法, 是油气藏综合解释过程中的关键步骤, 也是确定各目的层及出油点地震反射特征, 并进行储层沿层属性提取和横向预测的基础^{[2, 3, 4]}。例如：以G 147井NmⅡ -8小层日产油1.13 t的含油砂体为依据, 在合成地震记录上进行标定（图2）。在明化镇组中, 地震波在泥岩中的波阻抗大于砂岩中的波阻抗, 造成反射系数小于0, 储层的地震波波形响应为波谷形态, 同时, 通过卡住地震解释与该井钻遇层位NmⅠ 底, 最终确定NmⅡ -8小层在地震剖面中的位置。以G 147井NmⅡ -8小层为依据, 夹在上下两套泥岩之间的砂体含油后, 地震波的传播速度降低, 由地震剖面上表现为上下对称的两个相对较弱的反射波峰, 从而建立储层响应特征模型。由地震剖面得知, 地震波在NmⅡ -8小层的传播时间为1 016 ms, 深度为1 043 m, 因此可以计算出地震波在储层中的传播速度υ =2× 1 043/1 016=2.05 m/ms。

图2

	Figure Option ViewDownloadNew Window
	图2 G 147井地震合成记录标定

在港北地区完钻层位较浅, 浅层各层段保留地层较全, 基于该区这一构造特征开展井震结合确定各地震层位反射特征及构造格架, 由单井标定到多井联合标定进行地层精细划分和对比, 统一地质地震层位, 即建立井之间的横向特征对比桥梁模型, 能够解决较小幅度构造、砂泥岩薄互层等复杂的地质问题^{[5, 6, 7, 8]}。如图3所示, 橙色区域为目的层NmⅡ -8, 通过井震结合精细地层对比, 精准地完成了G 147井-X 36-11-4井的对比, 解决了由于构造幅度较大造成地层对比难的问题。通过开展井震结合精细地层对比, 统一地质地震层位, 将目的层段拉平, 剔除构造因素的影响, 对沉积现象进行逐层的解释, 进而勾绘出各种地质体和岩性变化。

图3

	Figure Option ViewDownloadNew Window
	图3 G 147井-X 36-11-4井地层对比

人工-智能解释技术是以样本点特征波形划分以及样本点层位分类标签概率化进行深度神经网络建模, 利用深度神经网络模型从样本点出发对其邻域点进行层位追踪概率预测以及样本点与追踪点层位时间差质控, 通过深度神经网络自动输出追踪遍历所有追踪点, 实现了复杂地质条件下的高精度解释, 确保后续沉积、储层、流体及成藏研究成果精度, 以降低井位实施风险。实现流程见图4a。

图4

	Figure Option ViewDownloadNew Window
	图4 人工-智能解释技术实现流程和目标区成果图

在给出电脑第一道解释剖面的解释样本后, 电脑会根据人工的层位解释方案、地震资料微构造的变化规律等情况, 自动推测出下一道范围内相对合理的解释方案。在这个过程中需要人工监控, 如果认为电脑解释结果不合理, 实时进行人工修正, 电脑会根据新的解释方案自动闭合确保可控性, 直至完成目标区域的层位解释。目标区成果图见图4b。

地震属性交会技术是应用识别储层的多种敏感属性求“ 交集” , 预测储层“ 甜点” 展布图。利用地震属性交会技术将岩性分成不同的组, 提取每一组的地震属性对这些属性做分类统计, 优选出最优属性^{[9, 10, 11, 12]}。以目标井区为例, 从上百种属性中依次选出5种属性进行交会, 每一种属性可以分为3个值域, 两两值域形成一种组合, 每种组合又对应1种颜色, 共出现243种颜色, 在这些颜色中优选出与储层展布最相关的属性, 通过分类统计可得到相应的敏感属性组合为：均方根振幅+弧长。

基于多元线性回归的地震属性融合技术, 既结合了单一属性对某种油层的识别优势, 又利用属性融合技术相互弥补各自属性局限性的特点, 剔除无关属性的干扰, 使有效信息最大化, 从而提高储层预测的精度^{[13, 14, 15]}。

3.5.1 时窗的选取

时窗大小的确定应根据地震反射波的视周期和研究区地质体的规模而定。当计算时窗过小时, 看不到一个完整的波峰或者波谷, 提取或计算出的某些地震属性并不真实; 当计算时窗过大时, 多个地震反射同相轴同时出现, 降低了对目的层的分辨率, 将使提取的地震属性具有多解性, 或者因包含的地质体类型过多导致所代表的地质含义难以确定。为准确地提取出目标储层的地球物理信息, 本文以层序为边界, 在等时地层格架控制下进行地层对比, 对层序界面进行追踪解释, 确定时窗, 开展针对目的层小时窗的地震属性分析研究, 时窗顶底定量的确定能够较好地约束地质体, 以便反映地下储层的真实情况。

3.5.2 属性融合计算

属性融合是指将两种属性融合于一体, 并对这两种属性进行数学运算变换, 同时结合这两种属性对储层的影响因素, 利用井位计算出两种地震属性的融合比重, 从而预测出有利储层的发育位置, 实现属性融合, 得到最优化的结果, 最终输出地震属性图。

以G 147井区为例, 目的层段地震资料主频为22 Hz, 有效频带为5~40 Hz, 据此可在有效频带内, 分别提取该砂体的均方根振幅（A）属性和弧长（B）属性, 通过已知井的情况得出比例系数A为2/3, B为1/3, 进行融合, 形成新的属性图。与常规方法相比, 均方根与弧长融合的属性储层预测精度更高（图5）。

图5

	Figure Option ViewDownloadNew Window
	图5 G 147井区NmⅡ -8小层常规属性与融合属性对比

3.5.3 储层预测精度验证

应用储层预测成果图中所钻遇储层的完钻新井井数占目标区中新井总井数的比重来验证储层预测精度。利用含油储层高精度预测技术, 实现了G 147井区滚动扩边, 并且针对明化镇组砂体在港北地区部署10口井位, 应用均方根振幅与弧长属性预测G 147井区储层情况与10口新井钻遇储层情况匹配关系良好, 储层预测精度达到100%。新井投产后, 初期日产量达到120 t, 井区目前累产油2.2× 10⁴ t。