[发明专利]一种地质体边界识别方法及装置

权利要求书

1.一种地质体边界识别方法，其特征在于，所述的地质体边界识别方法包括：

设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为N，获取2M+1个采样点，并根据所述2M+1个采样点计算高斯窗函数G(n)，M为自然数；

利用倾角扫描方法确定每个所述采样点所在的局部反射界面的最佳倾角对；

分别将每一所述采样点作为分析样点，并以所述分析样点为中心，根据所述分析样点的最佳倾角对从多个方向上分别读取2M+1个点的振幅值，生成多个地震信号序列S_k(n)，k＝1,2,...,i，n＝0,1,...,2M，i＝2M+1；

根据所述地震信号序列S_k(n)及高斯窗函数G(n)计算新地震序列G_Sk(n)；

计算所述新地震序列G_Sk(n)的广义希尔伯特变换值GH_k；

根据所述广义希尔伯特变换值GH_k生成所述分析样点的边界检测值；

根据多个所述的边界检测值生成一数据体，并对所述数据体进行沿层切片分析，识别地质体边界。

2.根据权利要求1所述的地质体边界识别方法，其特征在于，利用倾角扫描方法确定每个所述采样点所在的局部反射界面的最佳倾角对，包括：

设置一个以采样点(x，y)为中心点、长短轴分别为a、b的横向窗口；

通过所述横向窗口获取所述采样点的地震参数；

将所述地震参数代入采样点相干计算公式，获取所述采样点的相干值；

分别设置线方向、道方向地层视倾角的极值P_min，P_max，Q_min，Q_max；

按一定的步长ΔP、ΔQ确定地质体反射界面的视倾角P、Q；

根据奈奎斯特采样定律获取n_p*n_q个倾角对(P,Q)及对应的相干值c(t,P,Q)，其中，n_p＝(P_max-P_min)/ΔP，n_q＝(Q_max-Q_min)/ΔQ；

以视倾角对(P_L,Q_M)为中心点选取c*d个视倾角(P_l，Q_m)及其对应的相干值c(t,P_l,Q_m)拟合成一个曲面G(P,Q)，其中，P_l＝P_L±c*ΔP，Q_m＝Q_M±d*ΔP；

所述曲面G(P,Q)的最大值对应的倾角对(P,Q)为所述采样点的最佳倾角对。

3.根据权利要求1所述的地质体边界识别方法，其特征在于，所述的高斯窗函数为：

其中，α为标准偏差的倒数。

4.根据权利要求1所述的地质体边界识别方法，其特征在于，所述的新地震序列G_Sk(n)＝S_k(n)*G(n),k＝1,2,...,i。

5.根据权利要求1所述的地质体边界识别方法，其特征在于，根据所述广义希尔伯特变换值GH_k生成所述分析样点的边界检测值，包括：

计算所述广义希尔伯特变换值GH_k的虚部的绝对值IGH_k；

判断多个所述绝对值IGH_k中的最大值，并将所述最大值作为所述分析样点的边界检测值。