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

摘要

本发明公开了一种基于中心扩散加倾角属性的三维地震层位自动追踪方法，包括如下步骤：调整种子点位置；数据初始化准备；利用种子点的倾角属性确定映射点的位置，在映射点的时窗范围内寻找振幅极大值点，在每个振幅极大值点周围取一段波形数据，和种子点的波形数据进行相关系数计算，然后找出相关系数最大的映射点，并计算该映射点和种子点的振幅比；如果振幅比小于等于设定的门限值，则将该映射点确定为追踪点。本发明的积极效果是：利用中心扩散算法，有利于解决种子点提供不充分和种子点位置分布不均匀情况下，沿剖面追踪导致的层位大空洞；利用倾角属性在追踪平面的角度投影，更加精确了映射点的位置，提高了追踪的精确性效果。

主权项

一种基于中心扩散加倾角属性的三维地震层位自动追踪方法，其特征在于：包括如下步骤：步骤一、调整种子点位置：如果种子点位置不在局部振幅最大的位置则对种子点位置按如下方法进行调整：在种子点的时窗范围内寻找离种子点最近的极大值点，如果找到极大值点，则将种子点移动到极大值点的位置；如果未找到极大值点，则种子点位置不动；步骤二、数据初始化准备：首先，计算出整个三维地震数据体的梯度向量；其次，将各初始的种子点的位置放入队列queue中存放，并用一个数组array来标识每个种子点是否已经进行过相关追踪；步骤三、中心扩散追踪：1）计算种子点的倾角属性：利用步骤二计算出的梯度向量，计算出数据体中每个种子点的局部梯度方差矩阵，再求出该矩阵的特征向量，将特征向量投影到追踪平面上，得到特征向量在追踪平面的投影角度；2）确定映射点：根据第1）步确定的投影角度算出映射点的位置；3）确定追踪点：在第2）步确定的映射点的时窗范围内寻找振幅极大值点，在每个振幅极大值点周围取一段波形数据，和种子点的波形数据进行相关系数计算，然后找出相关系数最大的映射点，并计算该映射点和种子点的振幅比；判断振幅比是否大于设定的门限值，如果是，则放弃该映射点，并将该映射点所在的种子道标记为已追踪；如果否，则将该映射点确定为追踪点。