本征值分散，故Ec随之降低，反映子体的相似性变差。
[0062] (2)技术实现步骤
[0063] 技术实现流程如图1所示，步骤如下：
[0064] ①资料的预处理，主要针对前期的常规地震保幅处理工作
[0065] 资料的预处理主要针对前期的常规地震保幅处理工作，主要包括道编辑、带通滤 波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性补偿、叠前反權积及动校 正等。
[0066] ②宏面元抽取及不同方位均匀道集的形成
[0067] 该步骤主要为了考虑实际地震资料的叠前道集由于采集系统和成本，在叠前时差 提取时，总是存在着方位和偏移距不足的因素，因此，为了提高相邻道间信噪比，避免炮检 分布不匀带来的弊端，并保证在不同方位上有足够密度的不同炮检距道集分布、较一致的 叠加次数，通过扩大原CMP(共中必点）面元的手段建立CMP宏面元（一般为3X3或者 5X5)(对于一个地震道来说就是1X1，如果加上四周的两个道集，就是3X3,加上四周的 四个道集，就是5X5依次类推）。通过该方法可W对炮、检相邻的道进行部分叠加（可W选 择分方位或者分偏移距叠加，具体的方位大小、偏移距大小、W及方位和偏移距的间隔根据 实际地震资料有所不同）提高叠前资料信噪比，能够增强有效信号能量，并保证不同方位 叠加次数基本一致。实际宽方位采集的地震道集分布的偏移距-方位示意图如图2所示。
[0068] ⑨分方位全偏移距道集的叠加形成方位子体
[0069] 假如该地震工区均匀采集，一个CDP点（共反射点）的叠前地震道集的道数为N， 选取的宏面元尺度为3X3,则该CDP点的宏面元（就是通过步骤2获得的CMP宏面元）道 集道数为9XN道（即一个CDP点是N道，现在把周围的9个点全部算成该点的数据，送个 点就是9XN道，送是在假设采集均匀的情况下），该9XN道地震数据均对应一个偏移距和 一个方位角，因此可W得到一个9XN道的地震数据的偏移距-方位角分布图，如图3-1至 图3-6所示，在该图中，可W观察到最大满方位角覆盖的偏移距范围，因此通过选取若干个 方位角和全偏移距（方位子体如果为6个，就选择30间隔，方位子体为4个，就选择45度 间隔，偏移距选择一般为0米到满偏移距的最大值）进行叠加，进而得到数量不等（因为采 集观测系统的原因，方位角虽然是相等的，但是不同方位角满偏移距的最大值常常是不相 等的，从图2中可W看出，所W每个方位间隔从0米到满偏移距的地震道集的数量是不相等 的。）的分方位全偏移距的叠加（即把某个方位范围内的全偏移距叠前道集全部相加，再除 W道集的个数），就可W得到不同方位的地震叠加数据体，即地震方位子体数据（不同方位 的地震子体数据线号，道号均一样，地震道数，W及每个地震道的采样数均一样，只是由于 不同方位的叠加导致可能在某些方位的地震数据有所差异，图3-1至图3-6为6个方位地 震数据体，30°间隔）。生成的方位地震数据体如图3-1至图3-6所示。如果方位子体的个 数为Μ个，采集系统一般都是对称的，因袭考虑半圆圈180°即为全方位采集，此时具体的 方位间隔为180° /Μ。考虑到现实地震采集系统的限制，一般分为3个到6个不等不等。
[0070] ④地震方位子体分别的相干计算
[0071] 针对上一步骤形成的不同方位的地震数据子体数据，利用传统的Η维本征相干计 算其相干属性，形成若干个不同方位的地震数据相干体（不同方位的地震方位相干体的线 号，道号，地震道数W及采样均一样，只是由于不同方位的叠加导致可能在某些方位的相干 计算结果有所差异，附图中的方位子体为6个）。相干计算的原理为上述技术原理中利用本 征值相干的相干体计算方法（公式（1)到公式巧)），它的优点是：空间分辨率高，参与计算 的道集相互关系信息风度，具有较高的可信度和信噪比。而且不同级次的数据连续或者错 断，也表达了不同规模的多级断裂和裂缝的展布。生成的方位地震数据相干体如图4-1至 图4-6所示。
[0072] ⑤最优方位的地震数据相干体
[0073] 在若干个方位的地震数据相干体中，由于裂缝发育的延展性，可能在某些方位裂 缝的发育最为明显，但是在其他方位有可能无法得到展现，因此就要对全方位的地震数据 相干体进行寻优，可W采用的寻优算法有很多，比如最小二乘寻优，或者在各个方位子体 之间计算其楠圆属性然后找长轴的最大值寻优，对于相干算法而言，最简单也是最有效果 的一种寻优方法是在方位之间寻找最大值，具体如下：在若干个方位的地震数据相干体中 (比如有6个方位地震子体数据，该6个地震数据体的线号，道号均一样，地震道数，W及每 个地震道的采样数均一样，只是由于不同方位的叠加导致可能在某些方位的相干计算结果 有所差异），对每一道地震数据的每个点在各个方位之间寻找最大值（如果有六个方位，贝U 有六个点（因为方位的差异，每个方位计算出来的相干数据都是不相等），取送六个点中的 最大值即可），选择的最大值点地震道形成的地震方位最大相干数据体。相干数据体是裂 缝型储层中最有效也是最直接的属性计算方法，相干延伸或者变化较大的地方与地质中断 裂/裂缝发育的地方密切相关，因此计算得到的地震方位最大相干数据体就可W针对该地 震工区的裂缝发育情况进行预测。
[0074] 生成的方位地震数据相干体如图5所示。图5给出了利用本发明方法对某油田某 裂缝型储层发育区预测的结果图，从图中可W看出，其清晰地预测出了裂缝储层的整体发 育情况，经过与井中的FMI成像对比分析，吻合率较好。
[0075] 本发明针对实际油田勘探中面临的裂缝型储层的发育问题，采用的方法是针对实 际叠前地震资料开展特殊处理，然后利用叠前方位共中必点地震道集开展分方位全偏移距 进行叠加，叠加W后的不同方位的地震数据体再开展基于协方差矩阵的相干体计算程序， 再对计算得到的不同方位的地震相干体进行断裂识别，并进行寻优得到最大方位相干数据 体，进而得到最终的不同方位裂缝延展发育的预测成果图。
[0076] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种最优方位子体相干的裂缝预测方法，其特征在于：所述方法利用叠前地震道集 生成方位子体，然后分别计算不同方位之间的相干属性，最后通过寻优来获得最优方位的 地震数据相干体； 所述方法包括： (1) 叠前地震资料的预处理； (2) 宏面元道集抽取； (3) 分方位全偏移距道集的叠加形成方位子体； (4) 地震方位子体相干计算形成方位相干子体； (5) 计算最大方位相干子体，得到最优方位的地震数据相干体。2. 根据权利要求1所述的最优方位子体相干的裂缝预测方法，其特征在于：所述步骤 (1) 包括：道编辑、带通滤波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性 补偿、叠前反褶积及动校正。3. 根据权利要求2所述的最优方位子体相干的裂缝预测方法，其特征在于：所述步骤 (2) 是这样实现的： 将原CMP面元进行扩大形成CMP宏面元，具体如下：一个地震道是1X1,如果加上四周 的两个道集，就形成3X3的CMP宏面元，如果加上四周的四个道集，就形成5X5的CMP宏 面元，依此类推。4. 根据权利要求3所述的最优方位子体相干的裂缝预测方法，其特征在于：所述步骤 (3) 是这样实现的： 假如该地震工区均匀采集，一个CDP点的叠前地震道集的道数为N，选取的宏面元尺度 为3X3,则该⑶P点的通过步骤（2)获得的CMP宏面元的道集道数为9XN道，该9XN道地 震数据均对应一个偏移距和一个方位角，因此得到一个9XN道的地震数据的偏移距-方位 角分布图，在该图中，找到最大满方位角覆盖的偏移距范围，然后通过选取Μ个方位角和全 偏移距的地震叠前道集进行叠加，进而得到数量不等的地震叠前道集，就得到不同方位的 地震叠加数据体，即地震方位子体数据。5. 根据权利要求4所述的最优方位子体相干的裂缝预测方法，其特征在于：所述Μ的 取值为3个至6个，方位间隔为180°/Μ;所述偏移距取值为0米到满偏移距的最大值。6. 根据权利要求5所述的最优方位子体相干的裂缝预测方法，其特征在于：所述通过 选取Μ个方位角和全偏移距的地震叠前道集进行叠加，进而得到数量不等的地震叠前道集 是这样实现的： 把一个方位范围内的全偏移距地震叠前道集全部相加，再除以道集的个数就得到数量 不等的地震叠前道集。7. 根据权利要求6所述的最优方位子体相干的裂缝预测方法，其特征在于：所述步骤 (4) 是这样实现的： 利用三维本征相干计算方法计算不同方位的所述地震方位子体数据的相干属性，形成 若干个不同方位的地震数据相干体。8. 根据权利要求7所述的最优方位子体相干的裂缝预测方法，其特征在于：所述步骤 (5) 是这样实现的： 对步骤（4)得到的不同方位的地震数据相干体进行寻优，得到地震方位最大相干数据 体。9.根据权利要求8所述的最优方位子体相干的裂缝预测方法，其特征在于：所述寻优 的方法如下：在步骤（4)得到的不同方位的地震数据相干体中，对每一道地震数据的每个 点在各个方位之间寻找最大值，选择的最大值点地震道形成的地震方位最大相干数据体即 为最优方位的地震数据相干体。
【专利摘要】本发明提供了一种最优方位子体相干的裂缝预测方法，属于油田勘探领域。所述方法利用叠前地震道集生成方位子体，然后分别计算不同方位之间的相干属性，最后通过寻优来获得最优方位的地震数据相干体；所述方法包括：(1)叠前地震资料的预处理；(2)宏面元道集抽取；(3)分方位全偏移距道集的叠加形成方位子体；(4)地震方位子体相干计算形成方位相干子体；(5)计算最大方位相干子体，得到最优方位的地震数据相干体。
【IPC分类】G01V1/28
【公开号】CN105445787
【申请号】CN201410392355
【发明人】梁志强, 王世星, 张克非, 肖鹏飞
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年3月30日
【申请日】2014年8月11日