[发明专利]一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法、设备及可读存储介质

说明书

本发明公开一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法、设备及可读存储介质，其中，一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法包括计算待分析的地震数据的傅里叶变换频谱，确定待抽取的频率分量数量；计算待分析的所述地震数据的时频谱；根据所述频率分量数量和所述时频谱抽取所述地震数据的多频率属性；根据所述地震数据的大小确定融合多频率属性的参数；结合所述融合多频率属性的参数和t分布随机近邻嵌入算法融合所述多频率属性，以解释砂体厚度与空间展布。t分布随机近邻嵌入算法可以在缺少先验知识的情况下更好的拟合砂体厚度与多频率属性之间的复杂的非线性关系。

技术领域

本发明涉及油气勘探开发技术领域，特别涉及一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法、设备及可读存储介质。

背景技术

在油气开采技术领域中，砂体解释，尤其是砂体厚度解释和空间展布特性解释对于油气开采和井位部署至关重要。然而砂体性质与其地震响应之间的关系通常是复杂和非线性的，这增加了准确划定砂体边界、解释砂体厚度的难度。

在砂体解释方法中，地震属性分析技术被广泛应用。一些地震属性，例如振幅、频率、相位和纹理等属性，被应用于预测储层厚度和解释储层的地质结构。由于地震数据的非平稳性质和宽频带性质，仅仅使用单一地震属性很难准确圈定砂体分布并解释砂体厚度。为了解决使用单一地震属性分析的局限性，勘探地球物理学家开发和设计了许多基于多属性分析的技术和方法。红-绿-蓝(RGB)混合技术能够融合三种地震属性，进而可以用于解释不同厚度的储层。然而RGB混合技术是一个线性模型，并且只能混合三个地震属性，因此，RGB混合技术难以用于解决复杂的、非线性的砂体解释问题。

发明内容

基于此，有必要提供一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法、设备及可读存储介质，以解决现有技术中复杂的、非线性的砂体解释问题。

本发明提供的一种基于t分布随机近邻嵌入的多频率属性融合砂体解释方法，包括以下步骤：

计算待分析的地震数据的傅里叶变换频谱，确定待抽取的频率分量数量；

计算待分析的所述地震数据的时频谱；

根据所述频率分量数量和所述时频谱抽取所述地震数据的多频率属性；

根据所述地震数据的大小确定融合多频率属性的参数；

结合所述融合多频率属性的参数和t分布随机近邻嵌入算法融合所述多频率属性，以解释砂体厚度与空间展布。

进一步地，所述计算待分析的地震数据的傅里叶变换频谱，确定待抽取的频率分量数量的步骤包括：

读取三维的所述地震数据u(t，inline，xline)；

计算三维的所述地震数据的傅里叶频谱；

根据所述傅里叶频谱的带宽范围确定待抽取的频率分量数量M。

进一步地，所述计算待分析的所述地震数据的时频谱的步骤包括：

通过S变换计算三维的地震数据；

对于给定的三维的地震数据u(t，inline，xline)，其使用S变换转换得到的时频谱结果记为ST(τ，f)，S变换的定义为

其中，inline和xline代表三维地震数据的待处理的单道地震信号，t代表时间变量，f代表频率变量，τ揭示了高斯窗函数w(t)的位置信息。

进一步地，所述根据所述频率分量数量和所述时频谱抽取所述地震数据的多频率属性的步骤包括：