潜山储层的多尺度裂缝带预测方法、装置、介质及设备制造方法及图纸

本发明专利技术涉及一种潜山储层的多尺度裂缝带预测方法、装置、介质及设备，利用高斯曲率属性预测宏观尺度的潜山风化缝优势富集区，利用f

【技术实现步骤摘要】
潜山储层的多尺度裂缝带预测方法、装置、介质及设备


[0001]本专利技术涉及渤海海域潜山油气田开发
，具体涉及一种潜山储层的多尺度裂缝带预测方法、装置、介质及设备。

技术介绍

[0002]裂缝发育程度是控制潜山储层品质的最主要因素，潜山内幕发育小至微裂缝、大至断层的不同尺度裂缝系统，各尺度裂缝相互耦合，形成潜山缝网系统，影响潜山储层品质与连通性。裂缝预测结果对潜山类型油气田储量动用与井位部署具有重要影响，然而现有裂缝预测方法主要描述单一尺度的裂缝带，无法描述潜山多尺度裂缝带整体平面展布形态，造成利用单一尺度裂缝带预测结果指导储量动用与井位部署存在较大不确定性。

技术实现思路

[0003]本专利技术的目的是提供一种潜山储层的多尺度裂缝带预测方法、装置、介质及设备，以解决现有技术中利用单一尺度裂缝带预测结果指导储量动用与井位部署存在较大不确定性的问题。
[0004]为实现上述目的，本专利技术采取以下技术方案：
[0005]本专利技术提供一种潜山储层的多尺度裂缝带预测方法，包括：
[0006]精细追踪解释潜山顶面时间层位，形成时间网格，提取时间网格的高斯曲率属性作为油田宏观尺度的潜山风化缝优势富集区预测结果；
[0007]利用f
‑
k变换技术刻画潜山内幕高角度断裂，提取f
‑
k变换属性的沿层振幅切片，作为油田大尺度的构造缝富集区预测结果；
[0008]以油田区域上断裂发育方向与角度作为约束，利用蚂蚁体属性从原始地震数据中追踪解释微小断裂，提取蚂蚁体属性的沿层振幅切片，作为油田潜山中尺度断裂带预测结果；
[0009]利用小波变换技术分别将提取的高斯曲率属性、f
‑
k变换属性、蚂蚁体属性分解为低频分量与高频分量，通过融合规则将分解的多个低频分量和多个高频分量分别融合，并通过小波反变换得到三种属性的融合结果，作为油田潜山多尺度裂缝带综合预测结果。
[0010]进一步地，还包括对高斯曲率属性中小于阈值常数的数值处理方法：采用阈值处理公式对高斯曲率属性中小于阈值常数的数值进行阈值处理，阈值处理公式为：
[0011]其中，A为阈值后的高斯曲率属性，C为阈值前的高斯曲率属性，C0为设计的阈值常数，且阈值常数通过实际地震剖面标定获取。
[0012]进一步地，还包括利用f
‑
k变换技术刻画潜山内幕高角度断裂的方法：利用f
‑
k变换技术将原始地震数据变换至f
‑
k域，在f
‑
k域中对非高角度反射信号进行压制，通过反f
‑
k变换技术将压制的非高角度反射信号的结果反变换至时
‑
空域，得到高角度断裂刻画数据
体。
[0013]进一步地，还包括非高角度反射信号的压制方法：非高角度反射信号在f
‑
k域表现为与高角度反射信号在分布区域上存在差异，利用信号压制公式在f
‑
k域压制高角度反射信号分布区域之外的信号对非高角度反射信号进行压制，信号压制公式为：
[0014]其中，B为在f
‑
k域压制非高角度反射信号之后的结果，C为原始地震数据变换至f
‑
k域的结果，I为高角度信号变换至f
‑
k域后分布的区域。
[0015]进一步地，还包括多个低频分量的融合规则设计方法：基于拉普拉斯能量公式设计低频子带融合规则，低频子带融合规则可表达为：
[0016][0017]其中M
L
(i,j)为样点(i，j)处的低频子带系数融合结果，表示第k个属性的拉普拉斯能量总和，L
k
(i,j)代表第k个属性在样点(i，j)处低频子带系数。
[0018]进一步地，还包括多个高频分量的融合规则设计方法：基于属性高频子带系数的均方根能量配置高频子带融合权重计算公式设计高频子带融合规则，高频子带融合规则可表达为：
[0019][0020][0021][0022]其中，M
h
(i,j)，M
v
(i,j)，M
d
(i,j)分别表示样点(i，j)处的水平方向、竖直方向、对角线方向的高频子带融合结果，表示第k个属性在样点(i，j)处沿n特定方向的高频子带融合权重，为第k个属性在样点(i，j)处沿n特定方向的高频子带系数，n代表水平h、竖直v、对角线d三个方向中的一种。
[0023]基于上述的潜山储层的多尺度裂缝带预测方法，本专利技术还提供该方法的分析装置，包括：
[0024]第一处理单元，用于精细追踪解释潜山顶面时间层位，形成时间网格，提取时间网格的高斯曲率属性作为油田宏观尺度的潜山风化缝优势富集区预测结果；
[0025]第二处理单元，用于利用f
‑
k变换技术刻画潜山内幕高角度断裂，提取f
‑
k变换属性的沿层振幅切片，作为油田大尺度的构造缝富集区预测结果；
[0026]第三处理单元，用于以油田区域上断裂发育方向与角度作为约束，利用蚂蚁体属性从原始地震数据中追踪解释微小断裂，提取蚂蚁体属性的沿层振幅切片，作为油田潜山中尺度断裂带预测结果；
[0027]第四处理单元，用于利用小波变换技术分别将提取的高斯曲率属性、f
‑
k变换属性、蚂蚁体属性分解为低频分量与高频分量，通过融合规则将分解的多个低频分量和多个高频分量分别融合，并通过小波反变换得到三种属性的融合结果，作为油田潜山多尺度裂缝带综合预测结果。
[0028]基于上述的潜山储层的多尺度裂缝带预测方法，本专利技术还提供一种计算机可读存
储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的潜山储层的多尺度裂缝带预测方法的步骤。
[0029]基于上述的潜山储层的多尺度裂缝带预测方法，本专利技术还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述的潜山储层的多尺度裂缝带预测方法的步骤。
[0030]本专利技术由于采取以上技术方案，其具备以下有益效果：
[0031]本专利技术利用高斯曲率属性预测宏观尺度的潜山风化缝优势富集区，利用f
‑
k变换技术刻画潜山内幕高角度断裂，预测大尺度的构造缝富集区，利用蚂蚁体属性预测潜山中尺度的断裂带富集区，基于小波变换技术提取曲率属性、f
‑
k变换属性、蚂蚁体属性这三种尺度裂缝带预测属性的低频与高频分量，设计低频端与高频端融合规则，实现三种尺度裂缝带预测属性的融合，从而能够综合表征潜山多尺度缝网发育特征，进而指导储量动用与井网设计。
附图说明
[0032]通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本专利技术的限制。在整个附图中，用相同的附本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种潜山储层的多尺度裂缝带预测方法，其特征在于，所述多尺度裂缝带预测方法包括：精细追踪解释潜山顶面时间层位，形成时间网格，提取时间网格的高斯曲率属性；利用f
‑
k变换技术刻画潜山内幕高角度断裂，提取f
‑
k变换属性；以油田区域上断裂发育方向与角度作为约束，利用蚂蚁体属性从原始地震数据中追踪解释微小断裂，提取蚂蚁体属性；利用小波变换技术分别将提取的高斯曲率属性、f
‑
k变换属性、蚂蚁体属性分解为低频分量与高频分量，通过融合规则将分解的多个低频分量和多个高频分量分别融合，并通过小波反变换得到三种属性的融合结果，作为油田潜山多尺度裂缝带综合预测结果。2.根据权利要求1所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括对高斯曲率属性中小于阈值常数的数值处理方法：采用阈值处理公式对高斯曲率属性中小于阈值常数的数值进行阈值处理，阈值处理公式为：其中，A为阈值后的高斯曲率属性，C为阈值前的高斯曲率属性，C0为设计的阈值常数，且阈值常数通过实际地震剖面标定获取。3.根据权利要求1所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括利用f
‑
k变换技术刻画潜山内幕高角度断裂的方法：利用f
‑
k变换技术将原始地震数据变换至f
‑
k域，在f
‑
k域中对非高角度反射信号进行压制，通过反f
‑
k变换技术将压制的非高角度反射信号的结果反变换至时
‑
空域，得到高角度断裂刻画数据体。4.根据权利要求3所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括非高角度反射信号的压制方法：非高角度反射信号在f
‑
k域表现为与高角度反射信号在分布区域上存在差异，利用信号压制公式在f
‑
k域压制高角度反射信号分布区域之外的信号对非高角度反射信号进行压制，信号压制公式为：其中，B为在f
‑
k域压制非高角度反射信号之后的结果，C为原始地震数据变换至f
‑
k域的结果，I为高角度信号变换至f
‑
k域后分布的区域。5.根据权利要求1所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括多个低频分量的融合规则设计方法：基于拉普拉斯能量公式设计低频子带融合规则，低频子带融合...

【专利技术属性】
技术研发人员：杜昕，范廷恩，张显文，樊鹏军，马淑芳，高云峰，范洪军，张晶玉，乐靖，
申请(专利权)人：中海石油中国有限公司北京研究中心，
类型：发明
国别省市：