[發明專利]基于反射率法的多波AVO儲層彈性參數反演方法及系統

有效

• 申請號: 201510024810.5
• 申請日: 2015-01-19
• 公開（公告）號: CN104614763B
• 公開（公告）日: 2017-06-06
• 發明（設計）人: 李景葉;劉洪星;陳小宏;侯波;王守東;馬繼濤;劉國昌
• 申請（專利權）人: 中國石油大學(北京);中國石油天然氣集團公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三友知識產權代理有限公司11127
• 代理人: 湯在彥

1.一種基于反射率法的多波AVO儲層彈性參數反演方法，其特征是，所述的方法包括：

采集工區內的地震疊前道集、測井數據以及實際井旁角度域地震道集；

基于所述的地震疊前道集、測井數據以及實際井旁角度域地震道集確定地震子波以及振幅縮放因子；

采集工區內測井數據的統計模型參數；

根據所述的統計模型參數確定符合該工區的模型參數先驗分布函數；

采集地震構造解釋資料；

根據所述的地震構造解釋資料以及測井數據建立深度域的初始彈性參數模型；

根據所述深度域的初始彈性參數模型確定PP波與PS波反演殘差，該步驟包括：將所述深度域的初始彈性參數模型作為輸入；在頻率域通過傳遞矩陣算法確定每一個角度的PP波與PS波總反射率矩陣，所述PP波與PS波總反射率矩陣包含反射系數與傳播時間、層間多次波與透射損失；對所述PP波與PS波總反射率矩陣進行傅里葉反變換；將傅里葉反變換后的PP波與PS波總反射率矩陣與地震子波進行褶積，得到角度域的PP波與PS波角道集；采集實際角度域地震道集；將所述角度域的PP波與PS波角道集與實際角度域地震道集作差，得到PP波與PS波反演殘差；

構建最大后驗概率意義下的反演目標函數；

根據所述的反演目標函數以及PP波與PS波反演殘差確定最優彈性參數模型。

2.根據權利要求1所述的方法，其特征是，基于所述的地震疊前道集、測井數據以及實際井旁角度域地震道集確定地震子波以及振幅縮放因子包括：

根據所述的地震疊前道集以及測井數據、采用統計方法提取依賴于入射角度的地震子波；

將所述的測井數據作為輸入模型；

采用反射法正演模擬角度域的PP波與PS波道集；

將所述的PP波與PS波道集與所述的實際井旁角度域地震道集進行對比，得到振幅縮放因子；

將所述的振幅縮放因子應用于所述的地震子波。

3.根據權利要求1或2所述的方法，其特征是，根據所述的統計模型參數確定符合該工區的模型參數先驗分布函數包括：

所述的統計模型參數滿足三變量高斯分布；

分析所述的測井數據，得到所述統計模型參數的均值；

確定所述統計模型參數的自相關系數以及互相關系數；

根據所述的自相關系數以及互相關系數構建三參數相關的協方差矩陣；

根據所述的三參數相關的協方差矩陣形成符合該工區的模型參數先驗分布函數。

4.根據權利要求3所述的方法，其特征是，根據所述的地震構造解釋資料以及測井數據建立深度域的初始彈性參數模型包括：

根據所述的地震構造解釋資料、基于沉積模式建立地質模型；

將所述的測井資料按照構造模式進行插值和外推，得到每條測線的時間域的初始彈性參數模型，所述的初始彈性參數模型包括縱波速度模型、橫波速度模型以及密度模型；

通過時深轉換關系將所述時間域的初始彈性參數模型轉換為深度域的初始彈性參數模型。

5.根據權利要求4所述的方法，其特征是，構建最大后驗概率意義下的反演目標函數包括：

利用所述模型參數先驗分布函數對反演結果進行約束；

根據Bayesian原理、綜合反演似然函數以及先驗分布函數得到后驗概率分布函數；

根據所述的后驗概率分布函數構建反演目標函數。