[發明專利]一種地質體邊界識別方法及裝置

在審

• 申請號: 201510445467.1
• 申請日: 2015-07-27
• 公開（公告）號: CN105093300A
• 公開（公告）日: 2015-11-25
• 發明（設計）人: 孫勤華;潘建國;張虎權;王宏斌;孫東;劉曉梅
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王濤

技術領域

本發明是關于石油地球物理勘探地震數據的解釋技術，具體地，是關于一種能夠識別小規模地質體邊界的地質體邊界識別方法及裝置。

背景技術

目前，相干技術被廣泛用于檢測河道、斷裂。雖然相干技術經歷了三代發展，但是都存在平均效應，難以凸顯窄河道及小斷距斷裂。同時，相干算法對噪音敏感，在河道、斷裂相互切割處，地震信號雜亂、噪音強，相干技術難以準確識別窄河道和小斷裂，只能識別出河道、斷裂的大致發育區域。

現有的廣義希爾伯特變換邊界檢測技術的抗噪能力強，能夠凸顯小規模地質體邊界。在實際運用流程中，采用窗函數從主測線及聯絡測線兩個方向開展廣義希爾伯特變換檢測不同地質體邊界。然而，實際地質體(河道、斷裂)的走向往往空間分布復雜，只對主測線及聯絡測線兩個方向進行邊界檢測不能全面地檢測地質體邊界。

發明內容

本發明實施例的主要目的在于提供一種地質體邊界識別方法及裝置，以更加愛全面地檢測地質體邊界。

為了實現上述目的，本發明實施例提供一種地質體邊界識別方法，所述的地質體邊界識別方法包括：設置一長度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中的最大值，并將所述最大值作為所述分析樣點的邊界檢測值。