[發明專利]一種頁巖油儲層中夾層孔隙有效三維表征的方法及裝置

在審

• 申請號: 202011314324.4
• 申請日: 2020-11-20
• 公開（公告）號: CN114518309A
• 公開（公告）日: 2022-05-20
• 發明（設計）人: 徐二社;張文濤;錢門輝;李志明;陶國亮;申寶劍;劉偉新
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G01N15/08
• 分類號: G01N15/08
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，包括：

獲取甄選后的巖石樣品的實測孔隙度；

將所述甄選后的巖石樣品制備成與CT設備對應的巖石樣品；

獲取與所述CT設備對應的巖石樣品的灰度圖像；

對所述灰度圖像進行濾波降噪處理；

設置圖像處理的灰度閾值，對所述濾波降噪處理后的灰度圖像進行分割和二值化處理，初步識別所述巖石樣品的孔隙和骨架；

基于初步識別的孔隙和骨架，獲得所述巖石樣品的計算孔隙度；

基于所述計算孔隙度與所述實測孔隙度值的差值，以所述實測孔隙度值為標定，調整灰度閾值，獲取最終識別的孔隙和骨架；

基于所述最終識別的孔隙和骨架，獲得所述巖石樣品的孔隙度和孔徑分布特征。

2.根據權利要求1所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，所述基于所述計算孔隙度與所述實測孔隙度值的差值，以所述實測孔隙度值為標定，調整灰度閾值，獲取最終識別的孔隙和骨架包括：

計算所述計算孔隙度和實測孔隙度的差值，基于所述差值，以所述實測孔隙度值為標定，調整灰度閾值，獲得新的計算孔隙度；重復此過程直至所述差值小于預設閾值時；

在最后一次調整的灰度閾值下，獲取最終識別的孔隙和骨架。

3.根據權利要求2所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，所述基于所述最終識別的孔隙和骨架，獲得所述巖石樣品的孔隙度和孔徑分布特征包括：

基于所述最終識別的孔隙和骨架，建立所述巖石樣品的三維數字模型；

基于所述三維數字模型，獲得所述巖石樣品的孔隙度和孔徑分布特征。

4.根據權利要求1所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，將通過薄片鑒定或掃描電鏡選取的均一性巖石樣品作為甄選后的巖石樣品。

5.根據權利要求1所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，通過氦氣法孔隙度測試儀獲取甄選后的巖石樣品的實測孔隙度。

6.根據權利要求1所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，所述將所述甄選后的巖石樣品制備成與CT設備對應的巖石樣品包括：

根據所述CT設備的測量范圍，將所述甄選后的巖石樣品制備成與所述CT設備對應的尺寸。

7.根據權利要求1所述的頁巖油儲層中夾層孔隙有效三維表征的方法，其特征在于，所述CT設備為微米CT設備。

8.一種頁巖油儲層中夾層孔隙有效三維表征的裝置，其特征在于，包括：

氦氣法孔隙度測試儀，獲取甄選后的巖石樣品的氦氣法實測孔隙度；

CT設備，獲取與所述CT設備對應的巖石樣品的灰度圖像；

處理器，所述處理器與所述CT設備連接，所述處理器執行以下步驟：

對所述灰度圖像進行濾波降噪處理；

設置圖像處理的灰度閾值，對所述濾波降噪處理后的灰度圖像進行分割和二值化處理，初步識別所述巖石樣品的孔隙和骨架；

基于初步識別的孔隙和骨架，獲得所述巖石樣品的計算孔隙度；

基于所述計算孔隙度與所述實測孔隙度值的差值，以所述實測孔隙度值為標定，調整灰度閾值，獲取最終識別的孔隙和骨架；

基于所述最終識別的孔隙和骨架，獲得所述巖石樣品的孔隙度和孔徑分布特征。

9.根據權利要求8所述的頁巖油儲層中夾層孔隙有效三維表征的裝置，其特征在于，所述基于所述計算孔隙度與所述實測孔隙度值的差值，以所述實測孔隙度值為標定，調整灰度閾值，獲取最終識別的孔隙和骨架包括：

計算所述計算孔隙度和實測孔隙度的差值，基于所述差值，以所述實測孔隙度值為標定，調整灰度閾值，獲得新的計算孔隙度；

計算所述新的計算孔隙度與所述實測孔隙度值的差值，直至所述新的計算孔隙度與所述實測孔隙度值的差值小于預設閾值時，在最后一次調整的灰度閾值下，獲取最終識別的孔隙和骨架。

10.根據權利要求8所述的頁巖油儲層中夾層孔隙有效三維表征的裝置，其特征在于，將通過薄片鑒定或掃描電鏡選取的均一性巖石樣品作為甄選后的巖石樣品。