[發明專利]一種致密砂巖儲層孔隙度的預測方法

在審

• 申請號: 202210077982.9
• 申請日: 2022-01-24
• 公開（公告）號: CN114488301A
• 公開（公告）日: 2022-05-13
• 發明（設計）人: 肖曦;韓剛;何峰;葉云飛;王志紅;牛聰;黃饒;王健花;張益明;王迪
• 申請（專利權）人: 中海石油（中國）有限公司;中海石油（中國）有限公司北京研究中心
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京紀凱知識產權代理有限公司 11245
• 代理人: 謝斌

1.一種致密砂巖儲層孔隙度的預測方法，其特征在于，包括：

確定砂巖的縱波阻抗三維數據體；

確定劃分薄儲層和厚儲層的縱波阻抗門檻值b；

根據低于所述劃分薄儲層和厚儲層的縱波阻抗門檻值b的厚儲層所對應的縱波阻抗三維數據體，采用分斜率線性擬合的方式，獲得厚儲層的孔隙度的預測結果；

根據高于所述劃分薄儲層和厚儲層的縱波阻抗門檻值b的薄儲層所對應的縱波阻抗三維數據體，采用貝葉斯判別分析和協模擬的方式，獲得薄儲層的孔隙度的預測結果；

將所述厚儲層的孔隙度和所述薄儲層的孔隙度的預測結果合并成一個數據體，得到全儲層的孔隙度的預測結果。

2.根據權利要求1所述的致密砂巖儲層孔隙度的預測方法，其特征在于，所述確定砂巖的縱波阻抗三維數據體具體包括以下步驟：

根據目標區的測井曲線繪制縱橫波速度比與縱波阻抗的交會圖，以巖相曲線為色標，確定砂巖對應的縱橫波速度比的門檻值a，低于所述砂巖對應的縱橫波速度比的門檻值a為砂巖，高于所述砂巖對應的縱橫波速度比的門檻值a為泥巖；

通過疊前約束稀疏脈沖反演技術計算縱波阻抗三維數據體和縱橫波速度比三維數據體；

將所述縱波阻抗三維數據體和所述縱橫波速度比三維數據體進行交會，并根據交會結果確定縱橫波速度比低于所述砂巖對應的縱橫波速度比的門檻值a時所對應的縱波阻抗的三維空間分布范圍，并將其輸出的數據確定為砂巖的縱波阻抗三維數據體。

3.根據權利要求1所述的致密砂巖儲層孔隙度的預測方法，其特征在于，所述確定劃分薄儲層和厚儲層的縱波阻抗門檻值b具體包括以下步驟：

根據目標區的測井數據對縱波阻抗進行低通濾波處理，保留反演頻帶信息，將濾波后的縱波阻抗曲線與孔隙度進行交會，并以巖相曲線為色標，得到孔隙度測井曲線與濾波后的縱波阻抗測井曲線的交會圖，并從中獲得I類儲層、II類儲層、干層的縱波阻抗分布范圍；再根據孔隙度進行劃分，得到I類儲層、II類儲層與孔隙度之間的擬合關系；

參考孔隙度測井曲線與濾波后的縱波阻抗測井曲線的交會圖，并結合深度域顯示的原始縱波阻抗曲線、濾波后的縱波阻抗曲線和巖相曲線，確定一個劃分薄儲層和厚儲層的縱波阻抗門檻值b。

4.根據權利要求3所述的致密砂巖儲層孔隙度的預測方法，其特征在于，在所述獲得厚儲層的孔隙度的預測結果的步驟中，還應參考所述孔隙度測井曲線與濾波后的縱波阻抗測井曲線的交會圖中I類儲層、II類儲層與孔隙度之間的擬合關系。

5.根據權利要求2所述的致密砂巖儲層孔隙度的預測方法，其特征在于，在所述獲得薄儲層的孔隙度的預測結果的步驟中，還應根據所述縱橫波速度比三維數據體。

6.一種致密砂巖儲層孔隙度的預測裝置，其特征在于，包括

第一處理單元，用于確定砂巖的縱波阻抗三維數據體；

第二處理單元，用于確定劃分薄儲層和厚儲層的縱波阻抗門檻值b；

第三處理單元，用于根據低于所述劃分薄儲層和厚儲層的縱波阻抗門檻值b的厚儲層所對應的縱波阻抗三維數據體，采用分斜率線性擬合的方式，獲得厚儲層的孔隙度的預測結果；

第四處理單元，用于根據高于所述劃分薄儲層和厚儲層的縱波阻抗門檻值b的薄儲層所對應的縱波阻抗三維數據體，采用貝葉斯判別分析和協模擬的方式，獲得薄儲層的孔隙度的預測結果；

第五處理單元，用于將所述厚儲層的孔隙度和所述薄儲層的孔隙度的預測結果合并成一個數據體，得到全儲層的孔隙度的預測結果。

7.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，所述計算機程序被處理器執行時實現權利要求1-5中任一項所述的方法的步驟。

8.一種計算機設備，包括存儲器、處理器以及存儲在存儲器上并可在處理器上運行的計算機程序，其特征在于，所述處理器執行所述計算機程序時實現權利要求1-5任一所述的方法的步驟。