[發明專利]融合地震屬性顯示方法及裝置

在審

• 申請號: 201810730094.6
• 申請日: 2018-07-04
• 公開（公告）號: CN109061733A
• 公開（公告）日: 2018-12-21
• 發明（設計）人: 孟凡彬;郎玉泉;范雪輝;左衛華;吳建國;金學良;林建東;楊曙光;王剛;羅忠琴;劉鵬
• 申請（專利權）人: 淮北礦業（集團）有限責任公司;中國煤炭地質總局地球物理勘探研究院;中國石油集團科學技術研究院有限公司;新疆維吾爾自治區煤田地質局煤層氣研究開發中心
• 主分類號: G01V1/34
• 分類號: G01V1/34
• 代理公司: 北京超凡志成知識產權代理事務所(普通合伙) 11371
• 代理人: 蘇勝

1.一種融合地震屬性顯示方法，其特征在于，包括：

獲取關于目標地質體的頻率域數據體，所述頻率域數據體包括：所述目標地質體上多個目標層位參數點的頻率屬性數據，所述頻率屬性數據包括：頻率數據和位置數據；

對于每個所述目標層位參數點，將所述頻率數據分別代入預設的紅色頻帶計算公式、綠色頻帶計算公式和藍色頻帶計算公式，分別計算得到每個所述目標層位參數點的紅色頻帶、綠色頻帶和藍色頻帶；

根據每個所述目標層位參數點的所述紅色頻帶、所述綠色頻帶、所述藍色頻帶和預設的理論頻譜，確定每個所述目標層位參數點的融合頻譜；

顯示與每個所述目標層位參數點的所述位置數據對應的所述融合頻譜。

2.根據權利要求1所述的融合地震屬性顯示方法，其特征在于，所述根據每個所述目標層位參數點的所述紅色頻帶、所述綠色頻帶、所述藍色頻帶和預設的理論頻譜，確定每個所述目標層位參數點的融合頻譜，包括：

將每個所述目標層位參數點的所述紅色頻帶、所述綠色頻帶、所述藍色頻帶和所述理論頻譜代入預設的頻譜計算公式中，得到頻譜方程組；

根據所述頻譜方程組，確定頻帶權重向量；

對于每個所述目標層位參數點，將所述紅色頻帶、所述綠色頻帶、所述藍色頻帶和所述頻帶權重向量代入所述頻譜計算公式中，計算得到每個所述目標層位參數點的所述融合頻譜。

3.根據權利要求2所述的融合地震屬性顯示方法，其特征在于，所述根據所述頻譜方程組，確定頻帶權重向量，包括：

根據所述頻譜方程組，得到殘差平方和函數，所述殘差平方和函數的自變量為未知頻帶權重向量；

對所述殘差平方和函數關于所述未知頻帶權重向量進行微分且求最值，計算得到所述頻帶權重向量。

4.根據權利要求3所述的融合地震屬性顯示方法，其特征在于，所述獲取關于目標地質體的頻率域數據體，包括：

獲取關于所述目標地質體的地震數據體；

利用預設的構造導向濾波算法，對所述地震數據體進行濾波，得到濾波地震數據體；

利用預設的螞蟻追蹤算法對所述濾波地震數據體進行處理，得到螞蟻屬性數據體；

將所述螞蟻屬性數據體轉換為所述頻率域數據體。

5.根據權利要求4所述的融合地震屬性顯示方法，其特征在于，

所述紅色頻帶計算公式為：其中，k為常數，f為所述頻率數據，f_R為預設的紅色頻率窗的中心頻率，f_Bandwidth為預設的信號帶寬，b_R(f)為所述紅色頻帶；

所述綠色頻帶計算公式為：其中，f_G為預設的綠色頻率窗的中心頻率，b_G(f)為所述綠色頻帶。

6.根據權利要求5所述的融合地震屬性顯示方法，其特征在于，所述藍色頻帶計算公式為：其中f_B為預設的藍色頻率窗的中心頻率，b_B(f)為所述藍色頻帶。

7.根據權利要求6所述的融合地震屬性顯示方法，其特征在于，所述頻譜計算公式為：u(f)＝c_R*b_R(f)+c_G*b_G(f)+c_B*b_B(f)，將所述理論頻譜代入u(f)中。

8.根據權利要求7所述的融合地震屬性顯示方法，其特征在于，當所述目標層位參數點的個數為m時，其中m為大于1的整數，所述頻譜方程組為：其中，記記記C為所述未知頻帶權重向量，u(f_m)為第m個所述目標層位參數點的所述理論頻譜，b_R(f_m)為第m個所述目標層位參數點的所述紅色頻帶，b_G(f_m)為第m個所述目標層位參數點的所述綠色頻帶，b_B(f_m)為第m個所述目標層位參數點的所述藍色頻帶。