薄互层砂体厚度的预测方法及装置制造方法及图纸

本发明专利技术提供一种薄互层砂体厚度的预测方法及装置，方法包括：基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度。本发明专利技术能够有效提高薄互层砂体储层厚度的预测精度，进而提高勘探成功率，为后期勘探井位部署提供依据。

【技术实现步骤摘要】
薄互层砂体厚度的预测方法及装置
本专利技术涉及油气物探工程领域，具体涉及一种薄互层砂体厚度的预测方法及装置。
技术介绍
随着油田勘探开发进入到中后期的阶段，勘探难度逐渐增大。针对薄互层砂体的储层预测研究可以在一定程度上提高勘探成功率，减少勘探成本。薄互层砂体的储层厚度小于λ/4(λ为地震波长)时，薄互层砂体储层的不同反射界面对地震波的响应相互干涉，形成复合波，导致常规地震资料解释和反演对于薄互层砂体储层的厚度和范围预测精度较低。因此，亟需一种更高精度的薄互层砂体储层厚度的预测方式。
技术实现思路
针对现有技术中的问题，本专利技术提供一种薄互层砂体厚度的预测方法及装置，能够提高薄互层砂体储层厚度的预测精度，进而提高勘探成功率，为后期勘探井位部署提供依据。为解决上述技术问题，本专利技术提供以下技术方案：第一方面，本专利技术提供一种薄互层砂体厚度的预测方法，包括：基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度。进一步的，在所述基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量之前，还包括：获取地震数信号并对所述地震信号进行筛选处理得到筛选后地震信号；相对应的，所述基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量，包括：基于模态分量的预设数量对筛选后地震信号进行变分模态分解得到预设数量个模态分量。其中，所述对所述地震信号进行筛选处理得到筛选后地震信号，包括：剔除所述地震信号中频率大于预设频率的分量，得到筛选后地震信号。其中，所述根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据，包括：将所述瞬时频率域所述反射系数进行褶积处理，得到地震道数据。其中，所述基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度，包括：采用Teager能量算子计算所述地震道数据中各个地震道各自对应的薄互层的砂体厚度。第二方面，本专利技术提供一种薄互层砂体厚度的预测装置，包括：模态分解单元，用于基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；选取单元，用于选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；合成单元，用于根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；计算单元，用于基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度。进一步的，还包括：筛选单元，用于获取地震数信号并对所述地震信号进行筛选处理得到筛选后地震信号；相对应的，所述模态分解单元包括：模态分解子单元，用于基于模态分量的预设数量对筛选后地震信号进行变分模态分解得到预设数量个模态分量。其中，所述筛选单元，包括：筛选模块，用于剔除所述地震信号中频率大于预设频率的分量，得到筛选后地震信号。其中，所述合成单元，包括：褶积模块，用于将所述瞬时频率域所述反射系数进行褶积处理，得到地震道数据。其中，所述计算单元，包括：计算子单元，用于采用Teager能量算子计算所述地震道数据中各个地震道各自对应的薄互层的砂体厚度。第三方面，本专利技术提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的薄互层砂体厚度的预测方法的步骤。第四方面，本专利技术提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的薄互层砂体厚度的预测方法的步骤。由上述技术方案可知，本专利技术提供一种薄互层砂体厚度的预测方法及装置，通过基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度，能够提高薄互层砂体储层厚度的预测精度，进而提高勘探成功率，为后期勘探井位部署提供依据。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例中的薄互层砂体厚度的预测方法的第一种流程示意图。图2为本专利技术实施例中的薄互层砂体厚度的预测方法的第二种流程示意图。图3为本专利技术实施例中的薄互层砂体厚度的预测装置的结构示意图。图4为本专利技术实施例中的电子设备的结构示意图。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术提供一种薄互层砂体厚度的预测方法的实施例，参见图1，所述薄互层砂体厚度的预测方法具体包含有如下内容：S101：基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；可以理解的是，地震信号的信噪比较低，并且地震信号通常看作是由多种频率成分组成的复杂非平稳信号，需要对地震信号进行频谱分解，形成单频信号，通过频率能量与砂体厚度的关系，实现预测薄互层砂体厚度。在本步骤中，通过对地震信号进行变分模态分解(ViationalModeDecomposition，VMD)能够得到地震信号对应的各个模态分量。通过设置模态分量的预设数量，在进行变分模态分解能够得到预设数量个模态分量。其中，变分模态分解的表达式为：其中，μk(t)为地震信号的信号的第k个模态分量，{μk}为模态分量的集合，ωk为信号的第k个模态分量的中心频率，{ωk}为模态分量的中心频率的集合，代表对t求导，t为时间，δ(t)为狄拉克函数，π为圆周率，f(x)为地震信号，j为虚部。需要说明的是，预设数量的值越大，可以获得的模态分量个数越多，并且模态分量的中心频率越高。可以根据地震频带的最高值确定模态分量的预设数量。其中，在本实施例中，预设数量至少为两个最多为十个。S102：选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；在本步骤中，根据步骤S101中得到的预设数量个模态分量，确本文档来自技高网...

【技术保护点】
1.一种薄互层砂体厚度的预测方法，其特征在于，包括：/n基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；/n选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；/n根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；/n基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度。/n

【技术特征摘要】
1.一种薄互层砂体厚度的预测方法，其特征在于，包括：
基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量；
选取所述预设数量个模态分量中最大的瞬时频带能量并确定所述最大的瞬时频带能量对应的瞬时频率；
根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据；
基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度。


2.根据权利要求1所述的薄互层砂体厚度的预测方法，其特征在于，在所述基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量之前，还包括：
获取地震数信号并对所述地震信号进行筛选处理得到筛选后地震信号；
相对应的，所述基于模态分量的预设数量对地震信号进行变分模态分解得到预设数量个模态分量，包括：
基于模态分量的预设数量对筛选后地震信号进行变分模态分解得到预设数量个模态分量。


3.根据权利要求2所述的薄互层砂体厚度的预测方法，其特征在于，所述对所述地震信号进行筛选处理得到筛选后地震信号，包括：
剔除所述地震信号中频率大于预设频率的分量，得到筛选后地震信号。


4.根据权利要求1所述的薄互层砂体厚度的预测方法，其特征在于，所述根据所述瞬时频率和地震信号时间域的反射系数合成地震道数据，包括：
将所述瞬时频率域所述反射系数进行褶积处理，得到地震道数据。


5.根据权利要求1所述的薄互层砂体厚度的预测方法，其特征在于，所述基于所述地震道数据和Teager能量算子确定薄互层的砂体厚度，包括：
采用Teager能量算子计算所述地震道数据中各个地震道各自对应的薄互层的砂体厚度。


6.一种薄互层砂体厚度的预测装置，其特征在于，包括：
模态分解单元，用于基于模态分量的预设数量对地震信号...

【专利技术属性】
技术研发人员：吴一平，胡英杰，李龙，王晓辉，阎明，姜立，迟润龙，蓝阔，何浩瑄，周艳，李明，翁文胜，廖培淳，董颖鑫，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11