利索能及
我要发布
收藏
专利号: 2023102791016
申请人: 广东石油化工学院
专利类型:发明专利
专利状态:已下证
更新日期:2025-12-30
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种基于断层模型的油气藏搜寻方法,其特征在于,所述方法包括以下步骤:

S100:利用所探测区域断层的三维地震数据计算多种敏感性瞬时谱属性体,并选出最优的敏感性瞬时谱属性体;

S200:基于所述最优的敏感性瞬时谱属性体进行相干计算获得瞬时谱相干属性体;

S300:对所述瞬时谱相干属性体进行频谱分析,按照频率的设定阈值将目的层段的瞬时谱相干属性体划分为低频和高频;

S400:从各个低频目的层段的瞬时谱相干属性体对应的地层切片中识别出低频阴影相关的区域作为待优化区域,对所述地层切片中的待优化区域进行优化获得优化后的地层切片图像。

2.根据权利要求1所述的一种基于断层模型的油气藏搜寻方法,其特征在于,在S100中,利用所探测区域断层的三维地震数据计算多种敏感性瞬时谱属性体,并选出最优的敏感性瞬时谱属性体的方法为:从三维地震数据中获取瞬时振幅谱属性体、瞬时频率谱属性体和瞬时相位谱属性体;根据三维地震数据进行地层切片,并分别计算探测区域井孔储层厚度和对应地层切片之间的相关性系数;选择出相关性系数最大的敏感性瞬时谱属性体作为最优的敏感性瞬时谱属性体。

3.根据权利要求1所述的一种基于断层模型的油气藏搜寻方法,其特征在于,在S200中,基于所述最优的敏感性瞬时谱属性体进行相干计算获得瞬时谱相干属性体的方法为:通过相干计算得到瞬时谱相干属性体的相干属性值后选取相干属性值小于预设的相干阈值内的瞬时谱相干属性体;

根据标准道信号与所有地震道信号的关系,相干属性值C的计算方式如下:

其中,Sj(i)表示三维地震数据中的第i个样本点的第j道地震信号,M为地震道信号总道数,j为地震道信号索引,N为计算窗的总样本点数,i为计算窗的样本点索引。

4.根据权利要求1所述的一种基于断层模型的油气藏搜寻方法,其特征在于,在S300中,按照频率的设定阈值将目的层段的瞬时谱相干属性体划分为低频和高频的方法为:按照频率的设定阈值将目的层段的瞬时谱相干属性体划分为低频和高频,具体为小于设定阈值为低频,大于设定阈值则为高频。

5.根据权利要求4所述的一种基于断层模型的油气藏搜寻方法,其特征在于,令设定阈值为10Hz、15Hz、20Hz、25Hz、30Hz、35Hz中的任一频率。

6.根据权利要求1所述的一种基于断层模型的油气藏搜寻方法,其特征在于,在S400中,从各个低频目的层段的瞬时谱相干属性体对应的地层切片中识别出低频阴影相关的区域作为待优化区域,对所述地层切片中的待优化区域进行优化获得优化后的地层切片图像,具体步骤如下:S401,把各个低频目的层段的瞬时谱相干属性体对应的地层切片转换成灰度图;

S402,对各个灰度图进行边缘检测以检测到的边缘线对各个灰度图依次分割出若干个图像区域;以G为灰度图中所有图像区域的集合,i1为集合G中图像区域的序号,i1∈[1,N1]且为整数,N1为图像区域总数,Gi1表示集合G中第i1个图像区域;设置i1的初始值为1;设置一个空的集合A作为待优化区域集合;

S403,筛选集合G中与低频阴影相关的图像区域作为待优化区域,具体为:

令图像区域Gi1所在的灰度图为当前灰度图;

遍历i1的取值范围,计算图像区域Gi1在除当前灰度图外的各灰度图上的相同位置中计算每个图像区域的像素灰度值总和作为Gi1位置对应的在对应图像区域的能量值Eng,按图像区域对应的频率从低到高进行排序构成Gi1位置对应的图像区域将对应的能量值构成列表EngList;把不同图像区域对应的频率按其频率值从低到高进行排序构成频率列表FreList;

由于灰度图是不同频率的目的层段地层切片一一对应地转换得到的,因此灰度图数量与频率总个数相同;记j1为与Gi1位置对应的灰度图的每个图像区域的能量值在列表EngList中的序号或不同频率在列表FreList的序号,j1∈[1,D1]且为整数,D1为地层切片的频率总个数;EngList中的第j1个图像区域的能量值表示为EngListj1;FreList中的第j1个频率为与第j1个图像区域对应的频率表示为FreListj1;

记EngList中的第c个元素为图像区域Gi1的能量值,c为Gi1的能量值在EngList中的序号,图像区域Gi1的能量值为EngListj1=c,与其对应瞬时谱相干属性体的地层切片的频率为列表FreList的第c个元素FreListj1=c;

S404,根据频率从低到高按顺序逐一计算EngList中除了EngListc以外的其他能量值与EngListc之差的绝对值,记这些差的绝对值构成的列表为ChaList;

从ChaList中选出最大值ChaMax,之后,从ChaList中选出最小值ChaMin,记EngList中的第i2个能量值EngListi2和第j2个能量值EngListj2分别与EngListj1=c计算得到了ChaMin和ChaMax,则FreListi2和FreListj2分别为EngListi2和EngListj2的对应频率,其中i2为ChaMin对应的能量值在EngList中能量值的序号,j2为ChaMax对应的能量值在EngList中能量值的序号;

S405,通过ChaMin和ChaMax确定EngList中的频率变能区间,具体为:若i2<c或i2>c且j2‑i2>50,EngList的第i2到第j2个元素为频率变能区间,也就是与Gi1对应频率从FreListi2到FreListj2之间不同频率的图像区域的能量值;否则令i2的值为c,以EngList的第c到第j2个元素为频率变能区间,也就是与Gi1对应频率从FreListi2=c到FreListj2之间不同频率的图像区域的能量值;由EngList中频率变能区间之外的其他元素构成非频率变能区间;

计算频率变能区间/非频率变能区间的能量趋变系数η1/η2的方式为:

若η1>η2,则把图像区域Gi1放进待优化区域集合A中;

S406,当i1≤N1时,令i1的值增加1并转到步骤S403,否则输出待优化区域集合A,转到步骤S407;记k1为待优化区域集合A中图像区域的序号,则集合A中的第k1个图像区域为Ak1;

k1∈[1,B1]且为整数,B1待优化区域集合A中图像区域的总数;设置k1初始值为1;

S407,遍历待优化区域集合A中的每个图像区域Ak1,Ak1结合与Ak1位置相同的且在频率变能区间中的不同频率地层切片灰度图的图像区域以EngList中的能量值作为权重对包括Ak1在内的上述所有图像区域进行加权求和,以得到信息优化后的图像区域AOPk1,按以上步骤遍历完待优化区域集合A则完成了对筛选出的低频地层切片所对应灰度图的信息优化;

记T为包括Ak1在内的灰度图中的与Ak1位置相同的图像区域的列表;所以T中的第j1个元素可表示为Tj1,且T的第c个元素Tj1=c= Ak1;AOPk1的计算方式为:若i2≤c:

若i2>c:

其中, 表示能量值 对图像区域 中的像素值进行乘法运算;

S408,将信息优化后的优选的低频目的层段的瞬时谱相干属性体地层切片所对应的灰度图进行灰度变换以得到信息优化后的优选的低频目的层段的瞬时谱相干属性体地层切片,并输出该地层切片作为优化后的地层切片图像。

7.一种基于断层模型的油气藏搜寻系统,其特征在于,所述一种基于断层模型的油气藏搜寻系统包括:处理器、存储器及存储在所述存储器中并在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现权利要求1‑6中的任一项所述的一种基于断层模型的油气藏搜寻方法中的步骤。