[發明專利]疊前各向異性特征參數反演方法及計算機可讀存儲介質

有效

• 申請號: 201710507808.2
• 申請日: 2017-06-28
• 公開（公告）號: CN109143351B
• 公開（公告）日: 2020-07-24
• 發明（設計）人: 劉春園;陳天勝;陳雙全;季玉新;劉炯
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種疊前各向異性特征參數反演方法，包括：

基于疊前地震道集，獲得反射系數；

針對疊前地震數據進行疊前同時反演，分別獲得上層介質與下層介質的縱波速度、橫波速度與密度；

基于所述反射系數、所述上層介質與下層介質的縱波速度、橫波速度與密度，根據各向異性參數反演公式進行反演，獲得所述上層介質與下層介質的各向異性參數。

2.根據權利要求1所述的疊前各向異性特征參數反演方法，其中，基于疊前地震道集，獲得反射系數包括：

將所述疊前地震道集劃分為多個方位角道集；

基于速度模型將所述多個方位角道集轉換為多個入射角道集；

沿目標層針對每個入射角道集提取反射系數。

3.根據權利要求1所述的疊前各向異性特征參數反演方法，其中，所述針對疊前地震數據進行疊前同時反演，分別獲得上層介質與下層介質的縱波速度、橫波速度與密度包括：

對所述疊前地震數據進行疊前同時反演，獲得縱波速度體、橫波速度體與密度體；

將目標層劃分為所述上層介質與所述下層介質，在所述上層介質與所述下層介質中分別選定一個時窗，提取所述時窗所對應的縱波速度體、橫波速度體與密度體，并分別計算所述對應的縱波速度體、橫波速度體與密度體的均方根，作為所述上層介質與所述下層介質的縱波速度、橫波速度與密度。

4.根據權利要求1所述的疊前各向異性特征參數反演方法，其中，所述各向異性參數反演公式為：

其中，R_p表示反射振幅，i表示入射角，φ表示方位角，δ表示縱波變化系數，Δδ^(V)表示上下層介質的縱波變化系數之差，ε表示縱波各向異性參數，Δ ε^(V)表示上下層介質的縱波各向異性參數之差，γ表示橫波各向異性參數，Δγ表示上下層介質的橫波各向異性參數之差，表示縱波速度，Δα表示上下層介質的縱波速度之差，表示橫波速度，表示縱波垂直入射時的波阻抗，ΔZ表示縱波垂直入射時的波阻抗差，表示橫波的切向模量，ΔG表示橫波的切向模量差。

5.根據權利要求1所述的疊前各向異性特征參數反演方法，其中，還包括：基于所述上層介質與下層介質的各向異性參數，計算所述上層介質與下層介質的各向異性參數之差。

6.一種計算機可讀存儲介質，其上存儲有計算機程序，其中，所述程序被處理器執行時實現以下步驟：

基于疊前地震道集，獲得反射系數；

針對疊前地震數據進行疊前同時反演，分別獲得上層介質與下層介質的縱波速度、橫波速度與密度；

基于所述反射系數、所述上層介質與下層介質的縱波速度、橫波速度與密度，根據各向異性參數反演公式進行反演，獲得所述上層介質與下層介質的各向異性參數。

7.根據權利要求6所述的計算機可讀存儲介質，其中，基于疊前地震道集，獲得反射系數包括：

將所述疊前地震道集劃分為多個方位角道集；

基于速度模型將所述多個方位角道集轉換為多個入射角道集；

沿目標層針對每個入射角道集提取反射系數。

8.根據權利要求6所述的計算機可讀存儲介質，其中，所述針對疊前地震數據進行疊前同時反演，分別獲得上層介質與下層介質的縱波速度、橫波速度與密度包括：

對所述疊前地震數據進行疊前同時反演，獲得縱波速度體、橫波速度體與密度體；

將目標層劃分為所述上層介質與所述下層介質，在所述上層介質與所述下層介質中分別選定一個時窗，提取所述時窗所對應的縱波速度體、橫波速度體與密度體，并分別計算所述對應的縱波速度體、橫波速度體與密度體的均方根，作為所述上層介質與所述下層介質的縱波速度、橫波速度與密度。