图4为根据本发明实施例的多方向读取幅值过程的示意图；
[0022] 图5为根据本发明实施例的方向1的地震序列图；
[0023] 图6为根据本发明实施例的在两个数据体上沿志留系顶面做切片分析的分析结 果示意图；
[0024] 图7为根据本发明实施例的河道的地震剖面响应剖面示意图；
[0025] 图8为根据本发明实施例的地质体边界识别装置的结构示意图。
【具体实施方式】
[0026] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0027] 本发明实施例提供一种地质体边界识别方法及装置。以下结合附图对本发明进行 详细说明。
[0028] 本发明实施例提供一种地质体边界识别方法，如图1所示，该地质体边界识别方 法主要包括以下步骤：
[0029] 步骤SlOl :设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为N，获取2M+1 个采样点，并根据2M+1个采样点计算高斯窗函数G (η);
[0030] 步骤S102 :利用倾角扫描方法确定每个采样点所在的局部反射界面的最佳倾角 对；
[0031] 步骤S103 :分别将每一米样点作为分析样点，并以分析样点为中心，根据分析样 点的最佳倾角对从多个方向上分别读取2Μ+1个点的振幅值，生成多个地震信号序列S k (η)， k = 1，2, · · ·，i，η = 0, 1，· · ·，2Μ，i = 2Μ+1 ;
[0032] 步骤S104 :根据地震信号序列Sk (η)及高斯窗函数G (η)计算新地震序列Gsk (η);
[0033] 步骤S105 :计算新地震序列Gsk (η)的广义希尔伯特变换值GHk;
[0034] 步骤S106 :根据广义希尔伯特变换值GH1^成分析样点的边界检测值；
[0035] 步骤S107 :根据多个边界检测值生成一数据体，并对数据体进行沿层切片分析， 识别地质体边界。
[0036] 通过上述的步骤SlOl~步骤S107,通过窗函数获取采样点的最佳倾角对及多个 方向上的地震序列，并对该地震序列进行广义希尔伯特变换，从而获取边界检测值，识别地 质体边界。上述的识别方法具有较强的抑噪能力，能够更加准确、全面地识别各种地质环境 下的地质体边界。
[0037] 以下将结合上述各步骤，对本发明实施例的地质体边界识别方法进行详细说明。
[0038] 在上述的步骤SlOl中，设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为 N，在本发明实施例中，滑动窗口的长度的阶值N取1或2。通过该滑动窗口获取2M+1个采 Iv .Jc-N'fl：； 样点，并根据获取的2M+1个采样点计算高斯窗函数G(n)，其中，M为自 然数，η = 0, 1，. . .，2Μ，α为标准偏差的倒数，在本发明实施例中，α多3。
[0039] 上述的步骤S102,利用倾角扫描方法确定每个所述采样点所在的局部反射界面的 最佳倾角对。首先，根据Marfurt提出的相干计算方法，设置一个以采样点（X，y)为中心 点、长短轴分别为a、b的平面椭圆或者矩形的横向窗口，如图2所示，该平面窗口包含J道 地震道，分析点相干计算公式：
[0040] Λ-- i\. /.- 1
[0041] 其中，τ为采样点的时间（ms)或者深度（m) ;p，q分别为采样点所在的局部地质 界面沿X、y方向的视倾角，单位为ms/m(时间域）或m/m(深度域）；u(t, X，y)为横纵坐标 为（X，y)时间为t的地震数据点振幅值；uH (t, X，y)为地震道u (t, X，y)的希尔伯特变换；K 为长度包含w/Δ τ个样点的垂向窗口，该窗口大小为2w ms或者m，Δ τ为采样间隔，其单 位为ms或者m。
[0042] 在实际的层位解释工作中，设置线、道方向地层视倾角的极值P_min，P_max，Q_ min，Q_max，P_min〈p〈P_max，Q_min〈q〈Q_max。按一定的步长 Δ P、Δ Q 搜寻地质体反射界面 的视倾角P、Q，如图3Α及图3Β所示，按奈奎斯特（Nyquist)采样定律：Δ P < V(2af_)， Λ Q彡V(2bf_)，其中fmax为地震资料的最高频率。可以得到np*nq个倾角对（P，Q)及其 对应的相干值 c (t, P, Q)，其中 np= (P_max-P_min) / Δ P，n q= (Q_max-Q_min) / Δ Q。
[0043] 当存在c (t，PL，Qm)彡c (t，P，Q)时，即最大的相干值对应的（PL，Qm)为采样点所在 地震反射面对应的视倾角。由于按一定的步长（A P，AQ)离散采样取视倾角对（P，Q)，按步 长扫描过程中，即使c (t，Pb Qm)彡c (t，P，Q)，（Pb Qm)也并非是最佳的视倾角对。因此，在以 视倾角对（P^Qm)为中心点选取c*d个视倾角（PpQ n^P1= P上c*AP，Qn= QM±d*AP，及 其对应的相干值c (t，P1, QJ拟合成一个曲面G (P，Q)，G (P，Q)的最大值对应的倾角对（P，Q) 就是最佳的倾角对。
[0044] 上述的步骤S103 :分别将每一采样点作为分析样点，并以分析样点为中心，根据 分析样点的最佳倾角对从多个方向上分别读取2M+1个点的振幅值，生成多个地震信号序 列 Sk(η)。
[0045] 具体地，以采样点为中心，从九个方向分别读取2Μ+1个采样点的振幅值，形成9个 地震信号序列。并且将9个地震信号序列记为S k(n)，k= 1，2，...，9，η = 0，1，...，2Μ。其 中，9个方向具体可按如图4所示的方向获取：方向1、2为过采样点且与主测线和联络测线 平行；方向3为过采样点且与深度方向平行；方向4、5分别为过采样点的同一深度平面且 与主测线和联络测线呈45度夹角的方向；方向6、7为过采样点的且处于同一联络测线平面 上，与深度方向或者主测线方向呈45度夹角的两个方向；方向8、9为过采样点的且处于同 一主测线平面上，与深度方向或者联络测线方向呈45度夹角的两个方向。
[0046] 例如：在方向1上，窗口长度为9个采样点为例，如图5所示，中间点S1⑷为分析 样点，在方向1上，以分析样点为中心，在两侧分别选择4个采样点，并且从反方向第一个 点开始，读取整个窗口内的采样点的振幅，从而形成方向1上的地震信号序列S 1 (η)，η = 0，1，····，8〇
[0047] 需要说明的是，在本发明实施例中，是以从9个方向上分别读取采样点的振幅值 为例，但本发明并不限于此。
[0048] 在生成了上述多个方向上的地震信号序列Sk(n)及高斯窗函数G(n)后，通过上述 步骤S104,计算新的地震序列G sk (η)。具体地，该新的地震序列Gsk (n) = Sk (n) *G (n)，k = 1，2, · · ·，9, η = 0, 1，· · ·，2Μ。
[0049] 进一步地，执行步骤S105,计算新地震序列Gsk(n)的广义希尔伯特变换值GH k。 GHk为复数，对于9个复数分别计算其虚部的绝对值IGHk，k = 1，2, ...，9。