射角范围限定后的各个方位角叠加数据的叠加次数应大体一 致，以保证剖面的信噪比及能量均衡。入射角的设定范围是-针对目的层其最大入射角原 则上不大于30°，最小入射角应大于5°。方位角的个数原则上应大于或等于3个，具体可 视实际情况进行相应的设计。
[0063] iii、利用测井资料及各个方位角偏移数据体对井上的裂缝型储层进行准确标定， 并进行目的层的选取，目的层选取原则上是针对裂缝型储层段的底部进行选取。并确定好 目的层的反射特征，如波峰、波谷或零相位等这些容易进行层位追踪的反射特征。
[0064] iv、利用选取好的目的层的反射特征及相同的解释参数对各个方位角偏移数据体 分别进行层位拾取解释，解释参数包括起止线、道号及设定网格、反射特征，并对目的层位 进行网格化计算内插、平滑滤波等处理。网格化计算内插、平滑滤波等处理为将解释网格 (通常解释网格大于1线Xl道）的层位的双程走时值根据网格化插值算法内插成1线Xl 道，并对网格化计算后的层位双程走时数据设定一定的滤波参数进行平滑滤波，使其消除 某些走时值突变点，更接近实际的目的层的解释。网格化插值算法可有最小平方法、加权平 均法、泛克里金法等，使用这些计算方法应根据实际情况而定。
[0065]V、对各个目的层双程走时数据进行归一化处理，将其数值归一化校正到n~m值 域。归一化处理亦即是利用加、减、乘、除等组合运算，使得双程走时值均有一定程度上的变 小并使各点上的最大双程走时值与最小双程走时值的倍数相对增大，以利于进行椭圆拟合 计算，每个方位角的走时值是对应椭圆计算坐标系中的方位角-走时值。坐标系中X、Y轴 既是代表归一化后走时值数据，也代表方位角值。方位角值范围的设定是指当Y轴正半轴 (方位角值为0° )按顺时针方向旋转一圈时，方位角值范围即为0° -360°。亦即该坐 标系中的每个数据点对应两个数据值-分别是归一化后走时值大小及方位角值，该点与原 点之间的长度值为归一化后的走时值，且与原点的连线与Y轴正半轴夹角的大小为方位角 值。再对该坐标上的各点进行相关的椭圆拟合计算。具体如图5所示。
[0066] 其中，设定网格上各点分别进行地震玫瑰图及离散系数计算，综合并完成裂缝参 数平面图，包括以下步骤：
[0067]i、根据实际情况设置网格参数，这个网格称为设定网格。网格参数包括网格间距 和网格数等，网格参数的大小可根据要预测的裂缝点的网格分布及精度需要情况确定，原 则上通常设定的网格间距越大，则预测的精度相对较低，会丢失一些特征信息；网格间距越 小，则预测的精度相对较高，所绘结果越详细。设定网格的大小应根据实际及裂缝预测需要 情况而定，设定的网格通常是规矩的。
[0068]ii、提取设定网格点上经归一化处理后的各个目的层的双程走时值（t^并对该 点进行离散系数计算。该点上的离散系数可作为该点综合裂缝发育强度参数。
[0069] 上面步骤（32)中的离散系数计算公式如下：
【主权项】
1. 一种预测储层裂缝发育参数的方法，其特征在于：所述方法包括： (1) 输入地震数据； (2) 利用不同方位角的地震数据提取目的层双程走时值，并进行归一化处理，得到经归 一化处理后的各个目的层的双程走时值； (3) 对设定网格上各点分别生成地震攻瑰图和综合裂缝发育强度平面图，将此两图叠 合后进行显示，获得裂缝参数平面图。
2. 根据权利要求1所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（1) 中输入的地震数据为经动校正后形成的CDP道集数据。
3. 根据权利要求2所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（2) 包括W下步骤： (21) 、根据地震资料及实际情况，对处理后的CDP道集数据进行分方位角、入射角范围 限定的叠加及偏移处理，得到各个分方位角偏移数据体； (22) 、根据地震资料的实际情况及目的层的反射特征，对目的层使用相同的解释参数 对各个分方位角偏移数据体进行拾取解释，得到各个方位角数据体的目的层位的双程走时 值，并对拾取得到的各个方位角数据体的目的层位的双程走时值进行网格化插值和平滑滤 波处理，得到各个方位角数据体的处理后的目的层双程走时数据； (23) 、对所述各个方位角数据体的处理后的目的层双程走时数据进行归一化处理得到 经归一化处理后的各个目的层的双程走时值。
4. 根据权利要求3所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（3) 包括W下步骤： (31) 、根据实际情况设置设定网格的网格参数；所述网格参数包括网格间距和网格 数； (32) 、提取设定网格点上所述经归一化处理后的各个目的层的双程走时值，并对该点 进行离散系数计算；将该点上的离散系数作为该点综合裂缝发育强度参数； (33) 、计算设定网格点上的裂缝走向及该走向的裂缝相对强度，得到地震攻瑰图； (34) 、利用设定网格点上的离散系数的值经插值计算后生成综合裂缝发育强度平面 图； (35) 、将所述综合裂缝发育强度平面图和所述地震攻瑰图叠合在一张平面图上进行显 示，从而获得裂缝参数平面图。
5. 根据权利要求4所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（33) 中是采用H点楠圆拟合来进行计算的，具体如下： 楠圆的短轴方向为裂缝走向，其长轴与短轴之比值作为该走向的裂缝相对强度；该些 参数构成该网格点上的地震攻瑰图，综合各个网格点上的地震攻瑰图得到平面上各个网格 点上的裂缝走向及裂缝相对强度的分布情况。
6. 根据权利要求5所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（34) 是该样实现的： 将所述设定网格点上的离散系数的值经网格化插值计算后，形成新的网格和插值点， 并对该些数据点进行等值线绘制或形成颜色充填显示，从而得到综合裂缝发育强度平面分 布图。
7.根据权利要求6所述的预测储层裂缝发育参数的方法，其特征在于：所述步骤（35) 是该样实现的： 在一张平面图上设置两个Z值，一个Z值对应综合裂缝发育强度平面图中综合裂缝发 育强度的值，调整该值的范围能够控制综合裂缝发育强度的显示效果；另一个Z值对应地 震攻瑰图中裂缝相对强度的值，调整该值的范围能够控制各个地震攻瑰图中裂缝相对强度 值的显示效果；其中，叠合平面图上的X、Y坐标与综合裂缝发育强度平面图和所述地震攻 瑰图两图中的X、Y坐标是一样的。
【专利摘要】本发明提供了一种预测储层裂缝发育参数的方法，属于石油地球物理勘探领域。本方法包括：(1)输入地震数据；(2)利用不同方位角的地震数据提取目的层双程走时值，并进行归一化处理，得到经归一化处理后的各个目的层的双程走时值；(3)对设定网格上各点分别生成地震玫瑰图和综合裂缝发育强度平面图，将此两图叠合后进行显示，获得裂缝参数平面图。利用本发明方法可以准确预测出裂缝型储层。
【IPC分类】G01V1-28
【公开号】CN104678434
【申请号】CN201310631520
【发明人】胡伟光, 谢红斌, 印峰, 何鑫, 杨鸿飞, 赵卓男, 肖伟
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司勘探南方分公司
【公开日】2015年6月3日
【申请日】2013年12月2日