[发明专利]一种属性变化率逐次晋级的裂缝发育带检测方法

说明书

技术领域

本发明属于油田勘探开发领域，具体涉及一种属性变化率逐次晋级的裂缝 发育带检测方法。

背景技术

三维地震振幅属性的变化率作为一种非均匀储层描述的手段，已成为碳酸 盐岩地区缝洞型储层密集发育带及岩溶储层发育有利目标预测的重要技术手 段。但数学上连续的变化率计算由于被地质观测采样中不得不采用的离散抽样 的网格点所取代，不同的采样数据与采用不同的步长，使得导数的计算与横向 步长相关，以往采用的计算不能涵盖不同尺度分量的全体，而使很多有效区带 内圈定的振幅变化率区不能成带成片分布，不能准确描述不同尺度的计算结果， 而使裂缝有效发育带的圈定边界模糊、有效性不明确。

基于属性(或振幅)变化率的技术都是在地震叠后地震资料的基础上展开 的，但是传统的方法通常只是在相邻的几个方向、选用一套公式、在固定的网 格区间的地震道之间开展计算，这对于不同尺度、不同发育级别的变尺度裂缝 预测而言通常是不全面的、结果也是局限的，只代表大规模裂缝而不能获取从 小到大不同尺度的裂缝的全体。现有技术存在的缺陷是：1)要么只能在相邻的 一个方向的地震道之间开展计算，这对于具有一定方位的裂缝预测而言通常是 不准确的；2)要么只能选择一种高阶计算式(如4阶)，只能使用一种阵列(如 5x5点)、一列滤波算子计算裂缝发育，不能预测不同尺度的裂缝，尤其是小尺 度的裂缝，还要需要预知裂缝的发育大致尺度。

综上所述，属性变化率的计算都是在叠后地震资料上展开的，原先连续的 计算被有限的观测点限制而被离散化，传统的方法通常只能靠人工选取与观测 步长相关的一种固定的差分或高阶差分格式、在固定点阵的相邻地震道之间开 展计算，数学上这种处理技术对小尺度(高波数)的裂缝有压制作用，检测不 了不同尺度尤其是小尺度裂缝的相对变化，因此检测的成果体现不了裂缝的多 尺度性，在很多有效区带内只能断断续续检测结果，而不能检测成片成带的不 同裂缝的整体发育区域，不能精确描述平面空间裂缝发育的整体分布区域。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种属性变化率 逐次晋级的裂缝发育带检测方法，在继承传统的叠后地震资料的包括方位信息、 高阶变化率等信息的基础上，引入不同尺度、变高阶各种级次计算，并使这种 信息叠置起来，挖掘出从大到小不同级次裂缝发育信息，提供空间有效范围的 圈定，可以得到更为精确全面的裂缝发育有效带的预测结果。

本发明是通过以下技术方案实现的：

一种属性变化率逐次晋级的裂缝发育带检测方法，对叠后中心点地震道展 开变点阵的道集相关校正、对相邻点的地震道开展变尺度的高阶方位振幅变化 率计算、累计一定发育强度的不同尺度裂缝形成有效的裂缝发育区带，并最终 推广到整个工区，得到整个地震数据的最大方位高阶属性变化率属性、多阶逐 次晋级的不同尺度裂缝发育综合预测结果。

所述方法包括以下步骤：

(1)获取叠后三维地震数据体；

(2)第k阶计算所需道集方阵选择与坐标重置：对步骤(1)得到的叠后 三维地震数据体中的每一个地震道的每一个样点，取其为中心点，选择相邻纵 横向(k+1)*(k+1)相邻道，构建第k阶道集方阵，并将该中心点置为原点(0，0)， 然后对其相邻各点的坐标进行重置；

(3)计算中心道与相邻道相关时差：计算中心道与相邻每一道的互相关， 取其相关最大值时的时移量进行道间时差校正，消除地层起伏的影响；

(4)利用第k阶滤波算子计算不同方位属性变化率；

(5)提取第k阶不同方位属性变化率最大值；

(6)提取第k阶不同方位属性变化率平均值；

(7)提取各阶不同方位属性最大值与平均值：在每个时间样点重复步骤(2) -(6)，即可获得k＝2，4，6，…N各阶不同方位属性最大值与平均值；其中， N为用户选定的任意最大整数，可以选择任意高的阶数，但到达一定阶后，基本 就不变了，具体实施时，用N＝30效果较好；