一种地质体边界识别方法及装置的制造方法
【技术领域】
[0001] 本发明是关于石油地球物理勘探地震数据的解释技术，具体地，是关于一种能够 识别小规模地质体边界的地质体边界识别方法及装置。
【背景技术】
[0002] 目前，相干技术被广泛用于检测河道、断裂。虽然相干技术经历了三代发展，但是 都存在平均效应，难以凸显窄河道及小断距断裂。同时，相干算法对噪音敏感，在河道、断裂 相互切割处，地震信号杂乱、噪音强，相干技术难以准确识别窄河道和小断裂，只能识别出 河道、断裂的大致发育区域。
[0003] 现有的广义希尔伯特变换边界检测技术的抗噪能力强，能够凸显小规模地质体边 界。在实际运用流程中，采用窗函数从主测线及联络测线两个方向开展广义希尔伯特变换 检测不同地质体边界。然而，实际地质体（河道、断裂）的走向往往空间分布复杂，只对主 测线及联络测线两个方向进行边界检测不能全面地检测地质体边界。
【发明内容】
[0004] 本发明实施例的主要目的在于提供一种地质体边界识别方法及装置，以更加爱全 面地检测地质体边界。
[0005] 为了实现上述目的，本发明实施例提供一种地质体边界识别方法，所述的地质体 边界识别方法包括：设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为N，获取2M+1 个采样点，并根据所述2M+1个采样点计算高斯窗函数G(n)，M为自然数；利用倾角扫描方法 确定每个所述采样点所在的局部反射界面的最佳倾角对；分别将每一所述采样点作为分析 样点，并以所述分析样点为中心，根据所述分析样点的视倾角从多个方向上分别读取2M+1 个点的振幅值，生成多个地震信号序列S k (n)，k = 1，2,. . .，i，η = 0, 1，. . .，2M，i = 2M+1 ; 根据所述地震信号序列Sk (η)及高斯窗函数G (η)计算新地震序列Gsk (η);计算所述新地震 序列Gsk (η)的广义希尔伯特变换值GHk;根据所述广义希尔伯特变换值GH 1<生成所述分析样 点的边界检测值；
[0006] 根据多个所述的边界检测值生成一数据体，并对所述数据体进行沿层切片分析， 识别地质体边界。
[0007] 在一实施例中，上述的利用倾角扫描方法确定每个所述采样点所在的局部反射界 面的最佳倾角对，包括：设置一个以采样点（x，y)为中心点、长短轴分别为a、b的横向窗口； 通过所述横向窗口获取所述采样点的地震参数；将所述地震参数代入采样点相干计算公 式，获取所述采样点的相干值；分别设置线方向、道方向地层视倾角的极值P_min，P_max， Q_min，Q_max ;按一定的步长Δ Ρ、Δ Q确定地质体反射界面的视倾角P、Q ;根据奈奎斯特采 样定律获取np*nq个倾角对（P，Q)及对应的相干值c(t，P，Q)，其中，n p= (P_max-P_min)/ Δ P，nq= (Q_max-Q_min) / Δ Q ;以视倾角对（P D Qm)为中心点选取c*d个视倾角（P1, QJ及 其对应的相干值c (t，P1, QJ拟合成一个曲面G (P，Q)，其中，P1 = Pj c* Δ P，Qn= QM± d* Δ P ; 所述曲面G(P，Q)的最大值对应的倾角对（P，Q)为所述采样点的最佳倾角对。
[,/2 -r
[0008] 在一实施例中，上述的尚斯窗函数为：，其中，α为标准偏差的 倒数。
[0009] 在一实施例中，上述的新地震序列Gsk (n) = Sk (n) *G (n)，k = 1，2, · · ·，i。
[0010] 在一实施例中，上述的根据所述广义希尔伯特变换值6&生成所述分析样点的边 界检测值，包括：计算所述广义希尔伯特变换值64的虚部的绝对值IGH 〇判断多个所述绝 对值IGHk中的最大值，并将所述最大值作为所述分析样点的边界检测值。
[0011] 本发明实施例还提供一种地质体边界识别装置，所述的地质体边界识别装置包 括：高斯窗函数计算单元，用于设置一长度L为2M+1的滑动窗口，并设置长度L的阶值N， 获取2M+1个采样点，并根据所述2M+1个采样点计算高斯窗函数G (n)，M为自然数；最佳倾 角对确定单元，用于利用倾角扫描装置确定每个所述采样点所在的局部反射界面的最佳倾 角对；地震信号序列生成单元，用于分别将每一所述采样点作为分析样点，并以所述分析样 点为中心，根据所述分析样点的最佳倾角对从多个方向上分别读取2Μ+1个点的振幅值，生 成多个地震信号序列S k (n)，k = 1，2,. . .，i，η = 0, 1，. . .，2Μ，i = 2Μ+1 ;新地震序列计算 单元，用于根据所述地震信号序列Sk (η)及高斯窗函数G (η)计算新地震序列Gsk (η);广义 希尔伯特变换值计算单元，用于计算所述新地震序列Gsk (η)的广义希尔伯特变换值6&;边 界检测值生成单元，用于根据所述广义希尔伯特变换值6&生成所述分析样点的边界检测 值；地质体边界识别单元，用于根据多个所述的边界检测值生成一数据体，并对所述数据体 进行沿层切片分析，识别地质体边界。
[0012] 在一实施例中，上述的最佳倾角对确定单元具体用于：设置一个以采样点（X，y) 为中心点、长短轴分别为a、b的横向窗口；通过所述横向窗口获取所述采样点的地震参数； 将所述地震参数代入采样点相干计算公式，获取所述采样点的相干值；分别设置线方向、道 方向地层视倾角的极值P_min，P_max，Q_min，Q_max ;按一定的步长Δ Ρ、Δ Q确定地质体反 射界面的视倾角P、Q ;根据Nyquist采样定律获取np*nq个倾角对（P, Q)及对应的相干值 c(t,P,Q)，其中，np= (P_max-P_min)/AP，nq= (Q_max-Q_min)/AQ ;以视倾角对（PL,QM)为 中心点选取c*d个视倾角（P1, QJ及其对应的相干值c (LP1, QJ拟合成一个曲面G (P，Q)， 其中，P1= P上c*AP，Qn= QM±d*AP ;所述曲面G(P，Q)的最大值对应的倾角对（P，Q)为 所述采样点的最佳倾角对。
[0013] 在一实施例中，上述的高斯窗函数为：,其中，α为标准偏差的 倒数。
[0014] 在一实施例中，上述的新地震序列Gsk (n) = Sk (n) *G (n)，k = 1，2,. . .，i。
[0015] 在一实施例中，上述的边界检测值生成单元具体用于：计算所述广义希尔伯特变 换值64的虚部的绝对值IGH k;判断多个所述绝对值IGH k中的最大值，并将所述最大值作为 所述分析样点的边界检测值。
[0016] 本发明的有益效果在于，具有较强的抑噪能力，能够更加准确、全面地识别各种地 质环境下的地质体边界。
【附图说明】
[0017] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述 中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些 实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些 附图获得其他的附图。
[0018] 图1为根据本发明实施例的地质体边界识别方法的流程图；
[0019] 图2为根据本发明实施例的相干计算方法中的窗口示意图；
[0020] 图3A及图3B为根据本发明实施例的倾角扫描示意图；
[0021]