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

说明书

技术领域

本发明是关于石油地球物理勘探地震数据的解释技术，具体地，是关于一种能够识别小规模地质体边界的地质体边界识别方法及装置。

背景技术

目前，相干技术被广泛用于检测河道、断裂。虽然相干技术经历了三代发展，但是都存在平均效应，难以凸显窄河道及小断距断裂。同时，相干算法对噪音敏感，在河道、断裂相互切割处，地震信号杂乱、噪音强，相干技术难以准确识别窄河道和小断裂，只能识别出河道、断裂的大致发育区域。

现有的广义希尔伯特变换边界检测技术的抗噪能力强，能够凸显小规模地质体边界。在实际运用流程中，采用窗函数从主测线及联络测线两个方向开展广义希尔伯特变换检测不同地质体边界。然而，实际地质体(河道、断裂)的走向往往空间分布复杂，只对主测线及联络测线两个方向进行边界检测不能全面地检测地质体边界。

发明内容

本发明实施例的主要目的在于提供一种地质体边界识别方法及装置，以更加爱全面地检测地质体边界。

为了实现上述目的，本发明实施例提供一种地质体边界识别方法，所述的地质体边界识别方法包括：设置一长度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生成所述分析样点的边界检测值；

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

在一实施例中，上述的利用倾角扫描方法确定每个所述采样点所在的局部反射界面的最佳倾角对，包括：设置一个以采样点(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)为所述采样点的最佳倾角对。

在一实施例中，上述的高斯窗函数为：其中，α为标准偏差的倒数。

在一实施例中，上述的新地震序列G_Sk(n)＝S_k(n)*G(n),k＝1,2,...,i。

在一实施例中，上述的根据所述广义希尔伯特变换值GH_k生成所述分析样点的边界检测值，包括：计算所述广义希尔伯特变换值GH_k的虚部的绝对值IGH_k；判断多个所述绝对值IGH_k中的最大值，并将所述最大值作为所述分析样点的边界检测值。