[發明專利]少井區基于統計巖石物理建模的孔隙度預測方法及系統

有效

• 申請號: 201710876544.8
• 申請日: 2017-09-25
• 公開（公告）號: CN108037528B
• 公開（公告）日: 2019-08-30
• 發明（設計）人: 凡睿;李文成;蘇建龍;李宇平;繆志偉;楊鴻飛;陳靈君;孫均;趙卓男
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勘探分公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創大成知識產權代理有限公司 11614
• 代理人: 高爽

公開了一種少井區基于統計巖石物理建模的孔隙度預測方法及系統。該方法可以包括：基于地質資料與測井曲線，建立巖石物理模型；基于巖石物理模型，在預定范圍內將孔隙度增量設置為多個常數值，計算孔隙度加上孔隙度增量時速度、阻抗、彈性參數曲線；分別將孔隙度加相應的孔隙度增量定義為一類物相，建立速度、阻抗、彈性參數曲線交會圖，計算一類物相對應的概率密度函數；基于測井數據和疊前地震道集數據進行疊前地震反演，獲得速度、阻抗、彈性參數反演數據體，利用概率密度函數對反演數據體進行轉換，預測每一類物相所對應的孔隙度概率體。本發明在少井區有效開展孔隙度預測工作，提升少井區孔隙度預測精度，為勘探開發提供有力的技術支持。

技術領域

本發明涉及石油天然氣勘探開發領域，更具體地，涉及一種少井區基于統計巖石物理建模的孔隙度預測方法及系統。

背景技術

孔隙度是評價儲集性能的重要指標，因此對于儲層孔隙度的預測是地球物理勘探中面臨的重要科學難題。通?？紫抖鹊念A測方法主要是利用統計的方法開展，通過縱波速度、橫波速度、密度、波阻抗以及彈性參數測井曲線與孔隙度進行交會，優選出數據點分布較為集中，即兩者相關性較好的測井曲線，擬合兩者之間的關系式，通過對該參數進行反演以及利用該關系式進行轉換獲得孔隙度數據體，進而完成對高孔儲層預測的目的。但是該方法存在的問題是在縱波速度、橫波速度、密度、波阻抗以及彈性參數測井曲線影響的因素中孔隙度僅僅是一個方面，通過擬合公式的方法無法消除不同深度采樣點巖石組分差異、孔隙結構類型差異等因素的影響，因此該方法預測孔隙度的不確定性較高。其它孔隙度預測方法主要包括多元線性回歸方法和非線性預測(神經網絡、支持向量機等)方法，這一類的方法原理是建立井點位置處孔隙度曲線與井旁道地震屬性、反演數據體的線性和非線性關系，將該線性和非線性關系應用到整個地震、反演數據體預測得到孔隙度數據體，該類方法的缺陷在于對于鉆井樣本數有一定要求，只有在鉆井數滿足要求時才能夠起作用，因此此類方法不適用于少井區進行相關研究。因此，有必要開發一種少井區基于統計巖石物理建模的孔隙度預測方法及系統。

公開于本發明背景技術部分的信息僅僅旨在加深對本發明的一般背景技術的理解，而不應當被視為承認或以任何形式暗示該信息構成已為本領域技術人員所公知的現有技術。

發明內容

本發明提出了一種少井區基于統計巖石物理建模的孔隙度預測方法及系統，可以在少井區有效開展孔隙度預測工作，提升少井區孔隙度預測精度，為勘探開發提供更有力的技術支持。

根據本發明的一方面，提出了一種少井區基于統計巖石物理建模的孔隙度預測方法。所述方法可以包括：基于地質資料與測井曲線，建立巖石物理模型；基于所述巖石物理模型，在預定范圍內將孔隙度增量設置為多個常數值，并分別計算孔隙度加上孔隙度增量時相應的速度、阻抗、彈性參數曲線；分別將孔隙度加相應的孔隙度增量定義為一類物相，建立相應的速度、阻抗、彈性參數曲線交會圖，在交會圖上利用貝葉斯分類法計算所述一類物相對應的概率密度函數；基于測井數據和疊前地震道集數據進行疊前地震反演，獲得所述速度、阻抗、彈性參數反演數據體，利用所述概率密度函數對反演數據體進行轉換，預測每一類物相所對應的孔隙度概率體。

優選地，所述基于地質資料與測井曲線，建立所述巖石物理模型包括：利用羅伊斯-博伊特-希爾平均構建包含不同礦物混合物的巖石骨架模型；利用微分等效介質理論以及庫斯特-陶克瑞茲等效模量理論在所述巖石骨架模型中加入不同形狀孔隙，構建干巖石模型；通過蓋思曼流體替換在所述干巖石模型的孔隙中加入流體，獲取所述巖石物理模型。

優選地，所述孔隙包括鑄模孔、粒間孔、微裂縫和濕粘土柔性孔的至少其中之一。

優選地，所述測井曲線包括縱波速度曲線、橫波速度曲線、密度曲線、礦物組分測井解釋曲線、孔隙度測井解釋曲線及含水飽和度測井曲線。

優選地，還包括：基于所述縱波速度、橫波速度、密度測井曲線計算阻抗曲線與彈性參數曲線。

優選地，所述地震資料包括疊前地震道集與疊后數據。