[發明專利]基于Zoeppritz方程的粘彈介質地震反演方法

有效

• 申請號: 202110411461.8
• 申請日: 2021-04-16
• 公開（公告）號: CN113050162B
• 公開（公告）日: 2022-04-15
• 發明（設計）人: 潘新朋;葛子建;張鵬飛;李磊;柳建新
• 申請（專利權）人: 中南大學
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/48;G01V1/50;G01V1/30;G06F17/11
• 代理公司: 濟南圣達知識產權代理有限公司 37221
• 代理人: 祖之強

1.一種基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：包括以下過程：

獲取地震數據和測井數據；

根據獲取的地震數據，得到角度依賴的地震子波；

根據獲取的測井數據，得到地質層位模型先驗數據；

根據獲取的先驗數據得到地質層位模型參數初始值；

根據得到的地質層位模型參數初始值、角度依賴的地震子波和預設粘彈介質Zeoppritz方程，得到地質層位模型參數，根據得到的模型參數得到反演結果；

根據角度依賴的地震子波得到含反Q因子的復值地震波速度，將含反Q因子的復值地震波速度引入到預設粘彈介質Zeoppritz方程中表征各向同性粘彈介質，得到復值的粘彈介質精確Zeoppritz方程：

其中，

其中，V_P和V_S分別表示實數縱波速度和實數橫波速度，ρ是密度，θ₁為入射角，R_PP和T_PP分別表示縱波反射和透射系數，R_PS和T_PS表示橫波反射和透射系數，1/Q_P,1/Q_S分別表示縱波反Q因子和橫波反Q因子。

2.如權利要求1所述的基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：

采用柯西分布函數作為先驗分布函數，根據獲取的測井數據，得到地質層位模型先驗數據。

3.如權利要求1所述的基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：

基于工區內所有測井數據統計各模型參數的均值，求取各參數的自相關系數和互相關系數，構建參數相關的協方差矩陣，形成符合該工區的模型參數先驗分布函數。

4.如權利要求1所述的基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：

對測井數據進行平滑處理并沿解釋層位外推，獲得地質層位模型參數初始值。

5.如權利要求4所述的基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：

利用三維空間插值方法得到地質層位模型參數初始值，包括以下過程：

利用散點插值的方法對各個層位的數據進行插值，完成地質層位建模；

根據地質層位進行模型參數橫向插值，將測井信息進行橫向插值，計算得到地下每個點上的模型參數值。

6.如權利要求1所述的基于Zoeppritz方程的粘彈介質地震反演方法，其特征在于：

基于貝葉斯原理、先驗數據和正演算子構建最大后驗概率下的反演目標函數，通過求目標函數關于模型參數的導數，得到地質層位模型參數。