[发明专利]一种最优方向边缘监测方法

权利要求书

1.一种最优方向边缘监测方法，其特征在于：所述方法利用结构导向平滑滤波器结合结构导向边缘检测滤波提取三维地震数据中的地质体边界信息。

2.根据权利要求1所述的最优方向边缘监测方法，其特征在于：所述方法包括以下步骤：

(1)求取三维地震数据中的每个样点的梯度grad=dfdxdfdydfdt,]]>利用奇异值分解计算所述三维地震数据的结构张量T(x)＝grad^Tgrad＝λ_uuu^T+λ_vvv^T+λ_www^T；

其中，f表示原始三维地震数据，x，y，t分别表示所述三维地震数据的三个方向，λ_u、λ_v和λ_w分别表示结构张量T的奇异值，u、v和w分别为特征向量；规定特征值降序排列，λ_u≥λ_v≥λ_w≥0，其中特征向量u代表梯度最大的方向，即垂直于地质体平面特征；特征向量w代表梯度最小的方向，平行于线性地质体；

特征向量v和w都平行于地质体平面；

(2)利用结构导向平滑滤波器计算结构导向边缘检测滤波的分子：

所述结构导向平滑滤波器为(2)式：

g(x)-α▿·D(x)▿g(x)=f(x)---(2)]]>

其中f(x)为输入的三维地震数据，g(x)为平滑输出结果；常数变量α为时间步长，α＝0时，f(x)＝g(x)，D(x)与T(x)具有相同特征向量，不同特征值；

所述结构导向边缘检测滤波为(14)式：

svw,u=<<f>vw2>u<<f2>vw>u---(14)]]>

(14)式中，s_vw，u为计算得到的边缘系数；

(3)利用结构导向平滑滤波器计算结构导向边缘检测滤波的分母；

(4)计算步骤(2)得到的分子与步骤(3)得到的分母的比值，该比值即为地震数据中的地质体边界信息。