[发明专利]一种利用自适应分析时窗提取地震属性的方法、存储介质

说明书

本发明提出了一种利用自适应分析时窗提取地震属性的方法、存储介质，该方法包括获取目的层地震信号的平均瞬时频率；根据平均瞬时频率求取分析时窗的长度以及确定分析时窗的顶端和底面的原始层位面；利用趋势拟合目标函数对分析时窗的顶端和底面的原始层位面分别进行平滑处理，获得顶端和底面的优化层位面；保持分析时窗的长度不变，结合顶端和底面的优化层位面构建新的分析时窗；利用新的分析时窗作为顶底约束提取目的层的地震属性，以进行地震属性分析。该方法解决瞬时频率稳定性低、滤波插值损失精度的问题，提高了振幅类、层序类地震属性的分析精度。

技术领域

本发明属于石油天然气地震勘探与开发技术领域，尤其涉及一种利用自适应分析时窗提取地震属性的方法，还涉及一种计算机可读存储介质。

背景技术

由于地震属性中包含了储层物性变化、储层饱和流体成分、含油气储量等有关地震开发的重要信息，因此分析地震属性是储层预测工作中的一个非常重要的环节。

由于地震资料的频带范围是时变的，因此在对地震资料提取地震属性时，分析时窗的选择对于地震属性分析结果的影响很大。目前，常规的以计算点为中心的固定分析时窗法没有考虑地震资料的频率特征差异，导致在地震属性分析时会出现时窗跨越多个同相轴以及时窗长度不够且在同相轴零点附近会受到严重干扰的问题。即使时窗大小选择合适，也会受到相邻波组的影响。

近年来，有学者提出了求取地震自适应时窗的理念，并利用其开展相干体分析。该方法是利用不同的时频分析法，求得一个与瞬时频率相关的时窗长度，将其用于相干体分析中。但是，地震信号是不同子波的叠加之和，其瞬时频率无法精确求得，而且瞬时频率也容易受到噪声的影响而产生异常值，尽管通过简单的滤波和插值平滑可以在一定程度上去处这些异常值，但是又会降低自适应时窗的精度。虽然在相干体分析中允许使用精度降低的时窗，但是在振幅类、层序类地震属性的提取和分析过程中需要尽可能使用高精度的分析时窗。因此，如何提升瞬时频率的抗噪性能，同时减少时窗精度的损失就成为了获取合适的自适应时窗的关键问题。

发明内容

针对上述技术问题，本发明提出了一种利用自适应分析时窗提取地震属性的方法。所述方法包括：

S100，获取目的层地震信号的平均瞬时频率；

S200，根据平均瞬时频率求取分析时窗的长度以及确定分析时窗的顶端和底面的原始层位面；

S300，对分析时窗的顶端和底面的原始层位面分别进行平滑处理，获得顶端和底面的优化层位面；

S400，保持分析时窗的长度不变，结合顶端和底面的优化层位面构建新的分析时窗；

S500，利用新的分析时窗作为顶底约束提取目的层的地震属性，以进行地震属性分析。

根据本发明的实施例，上述步骤S100包括以下步骤：

S110，将目的层的地震信号转换为离散地震信号；

S120，利用小波变换对目的层的离散地震信号进行变换；

S130，基于下式求取变换后的目的层的离散地震信号的平均瞬时频率：

其中，FT(f，t)为变换后的目的层的离散地震信号，f为频率，t为时间，f_m为平均瞬时频率。

根据本发明的实施例，上述步骤S200中，根据下式确定分析时窗的长度：

其中，n为根据分析时窗的长度需求而选取的常数值，f_m为平均瞬时频率，T为周期，L为分析时窗的长度。

根据本发明的实施例，n的取值范围为0.1至2。