[发明专利]一种基于信息散度的地震数据属性提取方法

摘要

本发明公开了一种基于信息散度的地震数据属性提取方法，包括1)地震道的选择；2)地震道概率分布的计算；3)地震道Copula函数的计算；4)多道地震道信息散度的计算。本发明能实现地震属性中相干属性的提取，避免了在常规方法中的协方差矩阵的求解以及矩阵特征值的求解，降低以往相干属性提取方法在协方差矩阵和特征值求解时的计算时间。方便实现，能实现对地层地质的边缘、河道等不连续性检测和砂体的非均质性以及断层等地质特征的刻画，技术方案易于实现，可操作性强。

主权项

1.一种基于信息散度的地震数据属性提取方法，其特征在于，避免构建协方差矩阵，而是直接计算每个目标地震道数据与周围所划分的四个象限范围内各自两道及两道以上的地震道之间基于Copula函数的Renyi信息散度值，将所求值进行平均，并作为相干属性的值输出，在计算过程中避免构建协方差矩阵和矩阵特征值的求解，而是直接求取多地震道之间的信息散度，其具体包括如下步骤：1)地震道的选择选取以目标地震道为中心的3×3窗口的地震数据，以中心道为坐标原点进行划分为四个象限，将每个象限中的地震道包含中心道作为计算信息散度的数据源，以每个象限中的三角形的顶点为数据源；2)地震道概率分布的计算设地震道在窗内的时间采样点数为N，计算窗内每一道数据的概率分布P(x)：式中x表示地震道中的每一个数据，Number函数表示满足条件的个数；3)地震道Copula函数的计算C(u,v,w)表示联合概率密度函数，其中u，v，w分别表示地震道经过步骤1)得到的概率数据地震道，通过概率地震道将其值作为坐标，将三组进行对应，组合成三维坐标，将其对应的点投射到三维空间，其中每一维的坐标的范围为(0,1)，另外将其各自坐标分割为M，N，K块，统计每个块中所投射点的个数，计算出每个块中投射点个数与整体投射点个数的比值；4)多道地震道信息散度的计算根据Renyi信息散度公式(2)，可以得到公式(3)：其中R表示地震道之间的信息散度，α表示Renyi信息散度中的参数，f₁和f₂表示两个概率分布，M，N和K分别表示步骤3)中三维空间各自坐标上分块的个数，计算四个象限中各自的R(u,v,w)，并将其所有象限的值取平均后得到多道信息散度值作为目标道中目标点的相干属性值，通过上述步骤遍历整个地震数据体便可得到相干属性数据体。