[發明專利]一種基于局部峰值頻率的陸相薄儲層解釋方法及系統

在審

• 申請號: 201910644703.0
• 申請日: 2019-07-17
• 公開（公告）號: CN110456415A
• 公開（公告）日: 2019-11-15
• 發明（設計）人: 許輝群;郭淑文;周淑慧;邢興;佘曉宇;周宗良
• 申請（專利權）人: 中國石油大港油田勘探開發研究院;長江大學
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/28
• 代理公司: 42231 武漢智嘉聯合知識產權代理事務所(普通合伙)
• 代理人: 江慧<國際申請>=<國際公布>=<進入國

本發明公開了一種基于局部峰值頻率的陸相薄儲層解釋方法及系統，屬于地震勘探技術領域，解決了對薄層進行精細刻畫較難的問題。一種基于局部峰值頻率的陸相薄儲層解釋方法，包括以下步驟：將地震數據進行Q補償，得到經Q補償后的地震數據；對所述經Q補償后的地震數據進行廣義S變換，得到時頻數據體；獲取所述時頻數據體的局部時頻峰值所有對應的時刻，由局部時頻峰值所有對應的時刻確定薄層位置。實現提高了薄層的解釋精度。

技術領域

本發明涉及地震勘探技術領域，尤其是涉及一種基于局部峰值頻率的陸相薄儲層解釋方法及系統。

背景技術

隨油田開發的時間的變化，以及數據量增多，這些地震解釋成果與地質認識不一致，以及地質認識與實際工程開發結果存在矛盾，很顯然，要提供可靠的地震解釋成果需要地震解釋人員結合開發階段的各類數據進行綜合分析；

而陸相湖盆烴類儲層以碎屑巖占絕對優勢，儲層非均質性強，盡管開發階段井網密度很高，但僅靠井資料對井間關系進行推斷仍非常不確定，而地震信息雖具備橫向可追蹤的特性，但常規地震資料的分辨率低，難以用于對薄儲層(薄層)進行精細刻畫。

發明內容

本發明的目的在于至少克服上述一種技術不足，提出一種基于局部峰值頻率的陸相薄儲層解釋方法及系統。

一方面，本發明提供一種基于局部峰值頻率的陸相薄儲層解釋方法，包括以下步驟：

采集地震數據，將地震數據進行Q補償，得到經Q補償后的地震數據；

對所述經Q補償后的地震數據進行廣義S變換，得到時頻數據體；

獲取所述時頻數據體的局部時頻峰值所有對應的時刻，由局部時頻峰值所有對應的時刻確定薄層位置。

進一步地，所述將地震數據進行Q補償，得到經Q補償后的地震數據，具體包括：

利用公式對地震數據進行Q補償，得到經Q補償后的地震數據S_Q(t)，其中，r(τ)為反射系數序列，a_Q(τ，f)＝e^{-πfτQ+iH(πfτ/Q)}為Q的傳遞函數，H為在τ時刻對頻率f的希爾伯特變換。

進一步地，對所述經Q補償后的地震數據進行廣義S變換，得到時頻數據體，具體包括：

利用公式對所述經Q補償后的地震數據進行廣義S變換，得到時頻數據體ST(τ,f)，其中，λ和q是控制廣義S變換時頻分辨率的兩個參數，λ>1，q>1。

進一步地，獲取所述時頻數據體的局部時頻峰值所有對應的時刻，具體包括：

求取時頻數據體ST(τ，f)中頻率的極大值，此時對應的時刻τ則為時頻數據體的局部時頻峰值所有對應的時刻。

另一方面，本發明還提供一種基于局部峰值頻率的陸相薄儲層解釋系統，包括地震數據采集模塊、補償模塊、時頻數據體獲取模塊、薄層位置確定模塊；

地震數據采集模塊用于采集地震數據；

補償模塊用于將地震數據進行Q補償，得到經Q補償后的地震數據；

時頻數據體獲取模塊用于對所述經Q補償后的地震數據進行廣義S變換，得到時頻數據體；

薄層位置確定模塊用于獲取所述時頻數據體的局部時頻峰值所有對應的時刻，由局部時頻峰值所有對應的時刻確定薄層位置。

進一步地，時頻數據體獲取模塊利用公式對地震數據進行Q補償，得到經Q補償后的地震數據S_Q(t)，其中，r(τ)為反射系數序列，a_Q(τ，f)＝e^{-πfτQ+iH(πfτ/Q)}為Q的傳遞函數，H為在τ時刻對頻率f的希爾伯特變換。