一種地質體邊界識別方法及裝置的制造方法
【技術領域】
[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]