潜山储层的多尺度裂缝带预测方法、装置、介质及设备.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115963558A(43)申请公布日2023.04.14(21)申请号202211638994.0E21B49/00(2006.01)(22)申请日2022.12.20G01V1/32(2006.01)G06F18/25(2023.01)(71)申请人中海石油（中国）有限公司地址100010北京市东城区朝阳门北大街25号申请人中海石油（中国）有限公司北京研究中心(72)发明人杜昕范廷恩张显文樊鹏军马淑芳高云峰范洪军张晶玉乐靖(74)专利代理机构北京纪凯知识产权代理有限公司11245专利代理师谢斌(51)Int.Cl.G01V1/30(2006.01)权利要求书2页说明书7页附图5页(54)发明名称潜山储层的多尺度裂缝带预测方法、装置、介质及设备(57)摘要本发明涉及一种潜山储层的多尺度裂缝带预测方法、装置、介质及设备，利用高斯曲率属性预测宏观尺度的潜山风化缝优势富集区，利用f‑k变换技术刻画潜山内幕高角度断裂，预测大尺度的构造缝富集区，利用蚂蚁体属性预测潜山中尺度的断裂带富集区，基于小波变换技术提取曲率属性、f‑k变换属性、蚂蚁体属性这三种尺度裂缝带预测属性的低频与高频分量，设计低频端与高频端融合规则，实现三种尺度裂缝带预测属性的融合。本发明的方法能够综合表征潜山多尺度缝网发育特征，进而指导储量动用与井网设计。CN115963558ACN115963558A权利要求书1/2页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所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括多个低频分量的融合规则设计方法：基于拉普拉斯能量公式设计低频子带融合规则，低频子带融合规则可表达为：其中ML(i,j)为样点(i，j)处的低频子带系数融合结果，表示第k个属性的拉普拉斯能量总和，Lk(i,j)代表第k个属性在样点(i，j)处低频子带系数。6.根据权利要求5所述的一种潜山储层的多尺度裂缝带预测方法，其特征在于，还包括多个高频分量的融合规则设计方法：基于属性高频子带系数的均方根能量配置高频子带融合权重计算公式设计高频子带融合规则，高频子带融合规则可表达为：2CN115963558A权利要求书2/2页其中，Mh(i,j)，Mv(i,j)，Md(i,j)分别表示样点(i，j)处的水平方向、竖直方向、对角线方向的高频子带融合结果，表示第k个属性在样点(i，j)处沿n特定方向的高频子带融合权重，为第k个属性在样点(i，j)处沿n特定方向的高频子带系数，n代表水平h、竖直v、对角线d三个方向中的一种。7.一种分析装置，其特征在于，所述分析装置包括：第一处理单元，用于精细追踪解释潜山顶面时间层位，形成时间网格，提取时间网格的高斯曲率属性作为油田宏观尺度的潜山风化缝优势富集区预测结果；第二处理单元，用于利用f‑k变换技术刻画潜山内幕高角度断裂，提取f‑k变换属