按纵横、斜交4个相交方向、按表1所示对应2阶、4阶、6阶、……，k阶滤 波算子连续计算，将获得不同尺度下利用不同范围相邻地震道计算中心点地震道的变阶次 方位属性变化率；它们分别表示：不同方位（4个正交方位）、不同阶次（k = 2,4,6,…）的
[0084] ④不同尺度与不同方位多阶最大属性变化率计算方法
[0085] 对上述形成的数据集，在每次每点计算所获得的成果中，选取4个方位、2~k阶 (如2~20阶）最大值作为被计算点属性变化率最终值，通过步骤1到步骤3在整个三维 数据体上进行计算，即可获得不同尺度方位多阶最大属性变化率的计算结果，同时记录取 得最大值对应的点阵、方位，即可获得该点最优势裂缝发育强度、发育尺度（规模）、发育主 体方向的三维空间数据体，分析不同油藏部位裂缝发育的主体规模。
[0086] ⑤不同尺度与不同方位多阶属性变化率逐次晋级计算方法
[0087] 在第k阶属性变化率的计算过程中，对4个方位计算结果取其平均，获得第k阶属 性变化率的平均值，表示一种裂缝尺度下裂缝的平均发育强度0。如果对所有不同尺度k 并具有一定发育强度（裂缝发育强度大于某个预设值E)的裂缝发育都感兴趣的话，通过大 于预设值不同尺度裂缝发育强度的累计求和，就可获得全尺度的裂缝预测结果。
[0090] ⑥属性变化率中稳定性处理
[0091] 计算属性可以选择地震处理获得的各种数据类型，有效的属性包括：地震振幅、频 率、相位，地震反演获得的速度、弹性参数与岩性、储层与其它变化量，在计算前为抑制噪音 和使异常具有一定的延续性，可附加各种处理（三瞬，道积分微分+滤波、平滑）以抑制数 据的局部跳跃和极大个别值。
[0092] 图3-1至图3-4为塔河油田某区2、4、12、120阶次方位最大振幅变化率属性图，4 幅不同阶次方位振幅变化率带一致准确地指示了该区溶蚀缝洞发育有利分布区，仔细观测 可以发现随着阶次的提高，大尺度缝洞附近的小尺度缝洞异常带愈加明显。图4为该区叠 后多阶（2~12阶）方位振幅变化率逐次晋级合成图，其中每阶4个方向平均振幅变化率 必须大于预设值1E8,才能被累计叠加进入合成图，通过2~12阶逐次晋级形成了连续块体 状异常，指示明显边界的异常区块，反映了溶蚀洞缝区，而呈小线状、枝状异常则为与方位 相关部分，一般为裂缝发育区，从图可以看出这些线状、枝状异常总体分布在溶蚀洞缝区近 邻，为与缝洞伴生的次一级裂缝的有效发育区。
[0093] 本发明针对裂缝型油气田勘探开发中有效裂缝发育带圈定问题，设计一种利用地 震属性变化率的逐次晋级来检测裂缝发育带的方法，采用的方法是针对实际叠后地震资料 开展具有地层倾斜校正的、方位发育方向的、不同尺度多阶算子的、逐次累计有效属性变化 率的一种逐次晋级属性变化率的计算，并最终将该步骤推广到整个工区，得到整个地震数 据的属性变化率逐次晋级的裂缝发育预测结果，可应用于油田勘探开发中裂缝型油气藏有 效延伸范围的圈定。
[0094] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种属性变化率逐次晋级的裂缝发育带检测方法，其特征在于：所述方法对叠后中 心点地震道展开变点阵的道集相关校正、对相邻点的地震道开展变尺度的高阶方位振幅变 化率计算、累计一定发育强度的不同尺度裂缝形成有效的裂缝发育区带，并最终推广到整 个工区，得到整个地震数据的最大方位高阶属性变化率属性、多阶逐次晋级的不同尺度裂 缝发育综合预测结果。2. 根据权利要求1所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述方法包括以下步骤： (1) 获取叠后三维地震数据体； (2) 第k阶计算所需道集方阵选择与坐标重置：对步骤（1)得到的叠后三维地震数据 体中的每一个地震道的每一个样点，取其为中心点，选择相邻纵横向（k+1) * (k+Ι)相邻道， 构建第k阶道集方阵，并将该中心点置为原点（0,0)，然后对其相邻各点的坐标进行重置； (3) 计算中心道与相邻道相关时差：计算中心道与相邻每一道的互相关，取其相关最 大值时的时移量进行道间时差校正，消除地层起伏的影响； (4) 利用第k阶滤波算子计算不同方位属性变化率； (5) 提取第k阶不同方位属性变化率最大值； (6) 提取第k阶不同方位属性变化率平均值； (7) 提取各阶不同方位属性最大值与平均值：在每个时间样点重复步骤（2)-?)，即可 获得k = 2,4,6, -·Ν各阶不同方位属性最大值与平均值；其中，N为用户选定的任意最大 整数； (8) 提取大于预设值的累计多阶不同方位属性平均值的和：在每个时间样点上对k = 2,4,6, 一N各阶按步骤（7)获得不同方位属性最大值与平均值的数据集合，再取最大值数 据集合中的最大值，并取平均值数据集合中大于某个给定预设值的各个平均值进行累计求 和，即逐次晋级，得到该道该时间样点上多阶不同方位属性平均值累计和； (9) 输出不同尺度裂缝预测成果：对三维数据中的每一道中的每一个时间样点，重复 步骤（2)-(8)完成多阶不同方位属性平均值逐次晋级提取，然后将得到的多阶不同方位属 性平均值累计和输出，即为不同尺度裂缝发育综合预测结果。3. 根据权利要求2所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述步骤（1)是这样实现的： 对按点距叠前采集的地震数据进行常规处理后得到叠后成像地震资料，形成纵向、横 向、深度三个方向的叠后三维地震数据体。4. 根据权利要求2所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述步骤（2)中的对其相邻各点的坐标进行重置是这样实现的： 按与所述原点纵横方向相对单位步长计量对相邻各点的坐标进行重置，k的起始值为 2 ;对坐标重置可提高计算寻址的效率。5. 根据权利要求4所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述步骤（4)是这样实现的： 将该中心点的地震道与步骤（3)得到的相关校正后的相邻点的各个地震道，利用一阶 导数的第k阶滤波因子系数，计算纵横、斜交4个不同方向上的第k阶属性变化率，得到该 中心点在该时间样点上四个方向上的地震道的第k高阶属性变化率。6. 根据权利要求5所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述步骤（5)是这样实现的： 取步骤⑷中得到的四个方向上的地震道的第k阶属性变化率中的最大值，即为中心 点地震道在该时间样点上的第k阶最大方位属性变化率。7. 根据权利要求6所述的属性变化率逐次晋级的裂缝发育带检测方法，其特征在于： 所述步骤（6)是这样实现的： 取步骤（4)中得到的四个方向上的地震道的第k阶属性变化率中的平均值，即为中心 点地震道在该时间样点上的第k阶方位属性变化率。
【专利摘要】本发明提供了一种属性变化率逐次晋级的裂缝发育带检测方法，属于油田勘探开发领域。本方法对叠后中心点地震道展开变点阵的道集相关校正、对相邻点的地震道开展变尺度的高阶方位振幅变化率计算、累计一定发育强度的不同尺度裂缝形成有效的裂缝发育区带，并最终推广到整个工区，得到整个地震数据的最大方位高阶属性变化率属性、多阶逐次晋级的不同尺度裂缝发育综合预测结果。
【IPC分类】G06F19/00, G01V1/30
【公开号】CN105572732
【申请号】CN201410524363
【发明人】王世星, 梁志强
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年5月11日
【申请日】2014年10月8日