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

有效

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

本公開提供了一種基于Zoeppritz方程的粘彈介質地震反演方法，獲取地震數據和測井數據；根據獲取的地震數據，得到角度依賴的地震子波；根據獲取的測井數據，得到地質層位模型先驗數據；根據獲取的先驗數據得到地質層位模型參數初始值；根據得到的地質層位模型參數初始值、角度依賴的地震子波和預設粘彈介質Zeoppritz反射系數方程，得到地質層位模型參數，根據得到的模型參數得到反演結果；本公開基于地震疊前粘彈介質AVO精確Zeoppritz方程進行反演，反演計算精度高，反演參數目的性強，避免了分步反演或間接計算帶來的累積誤差。

技術領域

本公開涉及地震反演技術領域，特別涉及一種基于Zoeppritz方程的粘彈介質地震反演方法。

背景技術

本部分的陳述僅僅是提供了與本公開相關的背景技術，并不必然構成現有技術。

地震反演是獲得地下介質內部圖像、對儲層進行精細描述的有效方法，也是高分辨率地震勘探的最終表現形式，地震數據反演很大程度上提高了儲層表征的價值。隨著地震研究的重點由勘探逐漸向開發轉移，以及由常規油氣向非常規油氣尤其是頁巖氣的轉移，通過地震反演等手段來揭示地下油氣藏的精細分布特征，對油氣藏儲層進行精雕細刻的技術也受到越來越多的關注。在巖石物理理論的指導下，與儲層流體異常特征敏感的彈性參數可以表征為縱橫波衰減參數。從目前的研究來看，地震疊前反演研究主要涉及基于近似和精確Zoeppritz方程的AVO反演方法。

發明人發現，基于近似Zoeppritz方程的線性反演方法能穩定地解決病態反演問題或局部極值問題，但其在推導過程中的各種假設條件造成反射系數在中、大角度范圍內存在較大誤差，未能充分利用豐富的疊前信息，在反演精度上不能滿足實際油藏儲層精細表征要求，尤其對于實際三維大偏移距地震數據。針對各向同性粘彈介質，傳統的基于Zoeppritz方程近似反射系數是一個包含彈性參數和衰減因子的表達式，雖然通過等價改寫能重新表示為五個模型參數的線性表達式，但是五個未知參數造成的病態反演問題限制了方法的應用。并且，常規反演獲得的彈性參數難以直接對儲層流體導致的衰減特征進行刻畫，分步反演或間接計算會引入累積誤差。

發明內容

為了解決現有技術的不足，本公開提供了一種基于Zoeppritz方程的粘彈介質地震反演方法，基于地震疊前粘彈介質AVO精確Zeoppritz方程進行反演，反演計算精度高，反演參數目的性強，避免了分步反演或間接計算帶來的累積誤差。

為了實現上述目的，本公開采用如下技術方案：

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

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

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

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

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

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

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

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

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

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