[发明专利]基于中心扩散加倾角属性的三维地震层位自动追踪方法

权利要求书

1.一种基于中心扩散加倾角属性的三维地震层位自动追踪方法，其特征在于：包括如下步骤：

步骤一、调整种子点位置：

如果种子点位置不在局部振幅最大的位置则对种子点位置按如下方法进行调整：在种子点的时窗范围内寻找离种子点最近的极大值点，如果找到极大值点，则将种子点移动到极大值点的位置；如果未找到极大值点，则种子点位置不动；

步骤二、数据初始化准备：

首先，计算出整个三维地震数据体的梯度向量；其次，将各初始的种子点的位置放入队列queue中存放，并用一个数组array来标识每个种子点是否已经进行过相关追踪；

步骤三、中心扩散追踪：

1）计算种子点的倾角属性：

利用步骤二计算出的梯度向量，计算出数据体中每个种子点的局部梯度方差矩阵，再求出该矩阵的特征向量，将特征向量投影到追踪平面上，得到特征向量在追踪平面的投影角度；

2）确定映射点：根据第1）步确定的投影角度算出映射点的位置；

3）确定追踪点：

在第2）步确定的映射点的时窗范围内寻找振幅极大值点，在每个振幅极大值点周围取一段波形数据，和种子点的波形数据进行相关系数计算，然后找出相关系数最大的映射点，并计算该映射点和种子点的振幅比；判断振幅比是否大于设定的门限值，如果是，则放弃该映射点，并将该映射点所在的种子道标记为已追踪；如果否，则将该映射点确定为追踪点。

2.根据权利要求1所述的基于中心扩散加倾角属性的三维地震层位自动追踪方法，其特征在于：所述用数组array来标识每个种子点是否已经进行过相关追踪的方法是：初始化为0，表示没有进行过相关追踪；如果进行过相关追踪，则标识为1，表示进行过相关追踪。

3.根据权利要求1所述的基于中心扩散加倾角属性的三维地震层位自动追踪方法，其特征在于：根据第1）步确定的投影角度计算映射点位置的方法为：

首先按如下公式计算映射点和种子点相距的位置点数d：

d=h*tan(α)

其中，h表示种子道和相邻道点的位置距离，α为特征向量在追踪平面的投影角度；则，映射点的位置为a₀±d，其中，a₀为种子点的位置。

4.根据权利要求1所述的基于中心扩散加倾角属性的三维地震层位自动追踪方法，其特征在于：所述相关系数的计算公式为：

r=E([X-E(X)][Y-E(Y)])D(X)*D(Y)]]>

其中，X和Y分别是每个映射点的振幅值，E(X)和E(Y)分别是每段波形的平均振幅值，D(X)和D(Y)分别表示X和Y向量的方差。