断裂对构造圈闭的形成、油气运聚分析、储层物性改善等都具有重要意义, 控制着油气富集^[1]。随着中浅层油气勘探程度的不断提高, 深层已成为当今油气勘探的重点领域之一, 对深层地质认识、断裂预测精度的要求越来越高。深层地震资料信噪比较低, 背景噪声大, 对断裂及储层裂缝预测带来很大挑战。

地震解释人员对叠后地震数据断裂预测技术进行了诸多探索, 主要是基于相干^[1-2]、蚂蚁体^[3-4]、曲率类属性^[5]和边缘检测^[6]计算过程的技术优化, 通过地震数据滤波、平滑、分频融合、多属性等提高地震资料品质, 减小背景噪声, 改善断裂或储层裂缝的识别精度^[1-9]。蚂蚁体压制噪声能力较差, 改善断裂识别效果有限, 主要用于裂缝预测^{[3, 10]}。最大似然属性一般用于小尺度断裂、裂缝的预测^[11-13], 但是, 在构造导向控制下, 能够较好地预测出中、小级别的断裂, 有效指导断裂解释, 提高断裂解释的准确性和效率^[14-18]。

似然属性目前在油气和煤勘探领域的裂缝预测、潜山储层预测、隐形走滑断裂预测等得到初步应用^[11-19], 其断裂、裂缝预测效果相比常规相干、方差、蚂蚁体、基于方差体的蚂蚁体断裂预测技术的效果和精度明显得到提高, 能够有效预测不同尺度断裂, 但深层地震资料品质较差时, 断裂预测连续性较差, 整体较为破碎。一般情况下, 断裂与裂缝预测采用同一种技术方法、参数时, 不能有效区分主干断裂及微、小断裂, 两者之间存在相互干扰。本文针对珠江口盆地深水区白云凹陷深层地震资料信噪比较低、背景噪声严重导致断裂预测困难, 断裂组合关系不明, 储层裂缝预测精度较差的问题, 首次提出基于构造导向高清似然属性方法并探索了相关技术流程。通过对比分析相干、方差、蚂蚁体、基于方差体的蚂蚁体和最大似然属性断裂预测技术的优缺点, 进一步优化滤波技术组合、参数组合以降低噪声干扰、提高信噪比; 探索利用边缘增强迭代技术, 通过主干断裂和微、小断裂倾角、走向数据的差异性控制高清似然属性的计算, 达到区分主干断裂与裂缝预测的目的, 得到分别适用于主干断裂和裂缝预测的两个不同的高清似然属性数据, 突出需要的单一信息, 减小相互干扰, 进而有效地指导断裂解释、预测裂缝性储层发育区。

断裂复杂发育区, 各种级别断裂和裂缝广泛发育, 地层倾角、方位角变化较大, 中深层地震资料品质差, 随机噪声严重, 直接使用叠后资料预测断裂或者构造裂缝发育区, 背景噪声干扰严重, 不能有效识别主要断裂和裂缝发育区。

倾角和方位角是地质构造精准描述需要包含的两个方向属性, 也是制约地震数据体处理、断裂识别刻画精度的重要因素。倾角构造导向技术将倾角、方位角控制和滤波处理、属性计算相结合, 是一种在时窗范围内依据采样点倾角和方位角来追踪叠后地震数据同相轴进行滤波和多道属性处理的技术^[20-22], 可以在保留地震同相轴有效信息的基础上压制随机噪声、增强同相轴的连续性, 突出边缘特征, 提高地震同相轴追踪的准确性和连续性, 有利于构造的精细刻画和解释。利用倾角构造导向技术进行地震数据体处理或多道属性计算时, 首先利用地震数据体生成表示地质体构造形态(倾向和走向)空间展布特征的导向体, 确定每一个采样点的倾角和方位角信息, 然后在导向数据体的约束下确定滤波处理或者多道属性计算方向, 进行纯数据驱动的构造处理和多道属性计算^[23-25], 达到保护断裂信息消减随机噪声的目的, 并且能够有效提高滤波的效果及似然属性计算精度, 减小背景噪声干扰, 突出断裂、裂缝的预测效果。

目前相干、方差断裂识别方法主要通过对比时窗范围内的地震反射同相轴的波形、振幅和相位等相似性描述地层、岩性的横向不连续性以检测断裂; 蚂蚁体属性依据蚁群觅食方式形成的仿生算法, 可以实现自动优化搜索断裂信息。在深层地震信噪比较低、构造复杂时, 相干、方差及蚂蚁体等属性容易受背景噪声干扰, 无法准确识别断裂信息。

断裂似然属性(Likelihood)是一种以突出断裂信息为导向, 基于多道相似性第二代相干算法(C2)改进的断裂属性预测技术, 通过放大相邻样本点间相似性的对比关系, 有效增强断裂和非断裂属性响应的差异性, 压制背景噪声, 突出断裂成像效果^[18-19]。相邻采样点的相似性表示为:

式中: S为地震数据相邻采样点的相似属性(数值限定为0~1);g为地震数据体; 〈·〉_s表示地震数据体计算时窗内采样点构造导向平均值; 〈·〉_f表示沿断裂倾向、走向(倾角和方位角)平滑滤波, 增加相似性属性计算的稳定性。S值越小, 越能突出断裂造成的局部不连续强反射特征, 可以更加准确地预测断裂。

断裂似然属性计算公式为:

断裂似然属性(F_l, 表示断裂发育可能性大小)似然值范围为0~1。公式(2)中指数8扩大了相似性属性高值与低值之间的差异性, 使断裂特征更加突出, 背景噪声更小^[26-27]。

高清似然属性计算公式为:

式中: 高清似然属性(F_h)值范围为0~1;ω是M_l(最大似然属性)属性背景噪声最大值, 可以结合似然属性和实际地震剖面综合确定; F是迭代函数; n是迭代次数; D_E表示边缘检测增强技术参数。在不同倾角和走向的共同控制下, 通过边缘检测增强技术多次迭代可以增强断裂或者裂缝特征, 弱化噪声干扰, 突出主干断裂或者裂缝平面展布规律。

构造导向高清似然属性计算思路就是在对断裂进行增强性滤波处理、构造导向滤波基础上, 在断裂倾角、方位角(倾向和走向)约束下, 沿断裂面对时窗范围内扫描的倾向和倾角采样点分别计算其最大似然属性, 在背景噪声值约束下用Clip Attribute处理最大似然属性体, 消除背景噪声, 在不同倾角和走向约束下, 通过边缘检测增强迭代处理得到符合实际需求的高清似然属性, 从而更加准确地突出主干断裂信息, 增加断裂垂向和横向的延展性, 改善断裂的预测效果。基于构造导向高清似然属性断裂预测技术流程如图 1所示, 其关键步骤包括: ①增强断裂显示效果的构造滤波处理, 增强地震同相轴连续性, 整体压制地震资料噪声, 提高信噪比; ②基于初步处理的地震数据体计算倾角和方位角导向体, 利用构造导向体约束对地震数据进行滤波处理, 突出断裂位置, 形成构造导向滤波数据体; ③基于构造导向滤波数据体再次计算倾角和方位角导向体, 在倾角和方位角导向体约束下计算断裂似然属性及最大似然属性, 利用Tracking体自动统计分析不同规模断裂走向、倾角; ④结合原始地震剖面特征, 分析断裂似然属性体, 综合确定背景噪声值分布范围, 采用Clip Attribute对最大似然属性数据体进行处理, 在倾角和方位角控制下, 通过边界增强迭代处理得到符合实际需求的高清似然属性, 指导构造精细解释及微断裂发育区的预测。

图 1(Fig. 1)

图 1 基于构造导向高清似然属性断裂预测技术流程

白云凹陷位于珠江口盆地深水区, 整体呈NW-NE走向, 北邻番禺低隆起, 东靠东沙隆起, 西接云开低凸起(图 2), 水深为300~2 000 m, 自下至上发育始新统文昌组、恩平组, 渐新统珠海组、珠江组、韩江组、粤海组地层。目前浅层珠江组、韩江组已发现多个中型气田。近年针对深层珠海及恩平组钻探B1井和B2井, 测井解释出近100 m气层, 展现出深层极大的勘探潜力, 是珠江口盆地深水区未来的勘探方向, 但是井区储层渗透性较差, 需要寻找裂缝发育区进一步证实该区油气潜力。

图 2(Fig. 2)

图 2 白云凹陷研究区位置

研究区位于白云凹陷断裂叠置转换带, 发育张扭性质断裂, 北东向断裂与北西向断裂相互交织, 深层的珠海组、恩平组地震同相轴连续性差、反射较弱、复波较多, 同时受一统暗沙走滑断裂、火山底辟影响, 断裂垂向断距较小, 深层断裂不易准确预测(图 3)。从常规的相干、方差及其改进技术断裂预测切片结果(图 4a, 图 4b)可以看出, 背景噪声很大, 几乎不能有效识别出目的层段的断裂, 更无法指导构造裂缝发育区的优选。蚂蚁体及基于方差体的蚂蚁技术断裂预测的效果(图 4c, 图 4d)也不佳, 出现杂乱现象, 有效信息几乎全部被背景噪声淹没, 难以准确刻画断裂平面展布及相互关系, 不能有效反映裂缝实际发育情况。

图 3(Fig. 3)

图 3 白云凹陷深层断裂地震反射特征(位置见图 2)

图 4(Fig. 4)

图 4 研究区4 900 m深度不同技术方法的断裂预测切片结果 a 基于构造导向的相干; b 倾角校正方差体; c 蚂蚁体; d 基于方差体的蚂蚁体

首先对试验区叠后地震资料进行预处理, 为了突出研究区断裂刻画效果, 需要通过增强断裂构造滤波改善断裂显示, 平滑地震同相轴, 增强连续性, 减小地震数据因构造、地层快速变化产生的随机噪声, 降低因地震同相轴连续性差、反射弱所产生的干扰, 突出断裂与背景噪声的差异。

对比图 5a和图 5b可以看出, 增强断裂显示后地震体主要断裂特征得到有效改善, 地震同相轴明显平滑, 连续性更好。综合地震剖面统计分析, 滤波前背景噪声值为0.20~0.30, 断裂不相干值为0.33~0.45(图 5c); 滤波后背景噪声降低至0.07~0.20, 断裂不相干值为0.60~0.80, 减小了背景噪声干扰, 微、小断裂特征显著增强, 断裂预测明显改善, 连续性更好(图 5d)。

图 5(Fig. 5)

图 5 增强断裂构造滤波前、后4 900 m深度平面断裂识别效果 a 增强断裂构造滤波前地震剖面; b 增强断裂构造滤波后地震剖面; c 增强断裂构造滤波前断层似然属性; d 增强断裂构造滤波后断层似然属性

增强断裂构造滤波效果主要取决于Inline、Xline平滑参数及垂向采样点数。在保留断层滤波时, 根据研究目的层来选择参与计算的相干道数, Inline、Xline平滑参数越大, 同相轴平滑作用越强, 容易将有效信息滤掉, 宽度范围取7~11相同奇数道, 垂向采样点取11~21, 沿断面的三维倾斜中心窗进行滤波综合效果较好。

倾角、方位角可以全面直观地反映断裂的倾向与走向变化趋势, 对于研究区的构造特征具有一定的指导意义。在增强断裂构造滤波处理后, 计算地震体的倾角、方位角数据。在倾角、方位角数据体控制下, 对地震数据进一步进行构造导向滤波处理, 降低背景噪声, 平滑断面, 增强断面的分辨率, 改善断面的连续性, 突显断面特征。从构造导向滤波结果(图 6)可以看出, 断裂面更加聚焦、平滑, 断点更加干脆, 垂向连续性更好, 减小了主要断裂周围裂缝带对成像的影响。

图 6(Fig. 6)

图 6 构造导向滤波前(a)与构造导向滤波后(b)的地震剖面

在预处理地震数据体基础上, 计算地震数据体的倾角和方位角数据, 在倾角和方位角共同控制下, 计算断层似然属性及最大似然属性(图 7)。从断裂预测结果(图 7)可以看出, 断层似然属性基本刻画出了研究区北东向及北西向两组断裂分布特征, 但是断裂平面连续性较差, 断裂之间搭接关系受微、小断裂、背景噪声干扰严重; 最大似然属性表现出更多的断裂细节, 平面连续性较差, 主干断裂不突出, 断裂大小级别无法分辨, 不能有效指导断裂解释和优势构造裂缝发育区预测。

图 7(Fig. 7)

图 7 基于构造导向的断层似然属性(a)和最大似然属性(b)

利用最大似然属性Tracking数据体, 统计分析研究区主干断裂倾角为30°~60°, 走向为65°~105°; 微、小断裂倾角为30°~90°, 走向为90°~195°(图 8), 不同级别断裂之间倾角和走向具有明显的差异性, 微、小断裂整体倾角较大, 与主干断裂走向夹角成锐角。结合地震剖面及走向、倾角信息, 对断层似然属性及最大似然属性进行综合分析, 确定背景噪声范围为0~0.2, 为了突出主干断裂, 利用Clip Attribute对最大似然属性体进行处理, 直接消除背景噪声干扰, 但同时有效断裂信息也受到减弱。在此基础上, 采用边缘检测增强迭代处理技术, 针对主干断裂和微、小断裂采用不同参数的走向和倾角约束, 通过多次迭代处理, 达到分别刻画大断裂及微、小断裂的目的。在主干断裂倾角、走向参数共同约束下, 采用边缘检测增强技术进行3次迭代, 结果如图 9a所示, 与图 7a相比, 明显改善了主干断裂展布特征, 减小了微、小断裂对主干断裂的干扰, 突出了主干断裂信息。在微、小断裂倾角、走向参数共同控制下, 采用边缘检测增强技术进行2次迭代, 结果如图 9b所示, 与图 7b相比, 微、小断裂的分布更加清晰, 连续性更好, 减小了背景噪声的干扰, 直观地反映出北东向主干断裂之间发育协调性微小断裂, 能够有效指示构造裂缝发育优势区。高清似然属性(图 9)相对于常规的断层属性和最大似然属性(图 7)能够显著改善断裂识别效果, 减小噪声干扰, 可以有效刻画主干断裂和微、小断裂分布特征, 提高断裂识别精度。

图 8(Fig. 8)

图 8 主干断裂倾角及走向(a)及微、小断裂倾角及走向(b)

图 9(Fig. 9)

图 9 高清似然属性体主干断裂(a)及微、小断裂(b)预测结果

从地震剖面与不同技术方法预测的断裂融合对比结果(图 10)可以看出, 高清似然属性(图 10a)相对于基于构造导向的断裂似然属性(图 10b)能够清晰准确地反映出断裂发育位置, 断面成像更加聚焦; 与基于构造导向的最大似然属性(图 10c)相比, 垂向直线延伸、连续性好, 几乎无背景噪声干扰。倾角导向相干体(图 10d)能刻画出较大断裂, 不能识别微、小断裂, 垂向连续性较差, 背景噪声干扰严重, 不能清晰反映断面位置; 基于方差的蚂蚁体(图 10e)对断裂的预测效果整体较为杂乱, 背景噪声严重, 不能有效识别断裂。

图 10(Fig. 10)

图 10 地震剖面与不同技术方法预测的断裂的融合 a 高清似然属性; b 基于构造导向的断裂似然属性; c 基于构造导向的最大似然属性; d 倾角导向相干; e 基于方差的蚂蚁体

采用基于构造导向高清似然属性预测的白云凹陷研究区深层断裂见图 11, 相对于相干体、方差体属性(图 4a, 图 4b), 特征刻画清晰, 背景噪声小, 能够清晰刻画出北东和北西向展布的主干断裂及分支断裂, 断裂之间的搭接关系识别效果较好。微、小断裂的预测精度明显提高, 刻画出北东向主干断裂之间以及断裂搭接部位发育多条裂缝密集发育区, 这与背向叠置型转换带、断裂交汇位置为裂缝密集发育区的地质认识一致, 与已钻井B1和B2井区裂缝不发育相符合, 与地震剖面反射特征相同, 证实了该方法预测的正确性。常规的裂缝预测方法(图 4c, 图 4d)精度较低, 不能准确刻画出裂缝发育区, 不能有效地指导井位设计。在基于构造导向高清似然属性指导下, 完成研究区148条断裂精细解释(图 12), 落实4层15个断鼻-断块圈闭, 总面积近260 km²; 刻画了主干断裂间协调性微、小断裂, 突出了分支断裂及断裂带内部特征, 指导优选出3个深层构造裂缝发育优势区(图 11b), 为油气运聚研究和井位的优化设奠定了构造和储层基础。

图 11(Fig. 11)

图 11 采用基于构造导向高清似然属性预测的白云凹陷研究区深层(4 900 m)主干层断裂(a)及微、小断裂(b) (红色方框内为试验区, 红色圈内为构造裂缝发育优势区)

图 12(Fig. 12)

图 12 白云凹陷研究区深层目的层平面断裂组合前期成果(a)与最新成果(b)

1) 综合利用断裂增强构造滤波、构造导向滤波属性, 能够有效去除背景噪声, 保留断裂原始边界特征, 显著改善地震资料品质, 结合似然属性体和地震剖面确定随机噪声阈值, 采用Clip Attribute属性处理可以进一步去除噪声, 提升深层断裂和裂缝性储层的预测精度。

2) 基于最大似然属性提取的Clip Attribute属性体, 利用大断裂与微、小断裂的倾角和走向的差异性, 在不同参数控制下利用边缘检测增强迭代处理技术, 可以有效区分主干断裂及微、小断裂, 突出单一信息, 减小相互干扰, 提高深层断裂预测精度。解决了运用相干、似然属性技术预测深层断裂呈现出杂乱、连续性差无法准确区分主干断裂与微、小断裂的问题, 可以有效地指导断裂解释和裂缝发育区的优选。