[發明專利]巖石儲集層的電阻率計算方法、裝置及存儲介質

有效

• 申請號: 201910071330.2
• 申請日: 2019-01-25
• 公開（公告）號: CN111487685B
• 公開（公告）日: 2022-12-02
• 發明（設計）人: 何家歡;周克明;段勇;肖富森;楊雨;李農;謝邦華
• 申請（專利權）人: 中國石油天然氣股份有限公司;四川科力特油氣技術服務有限公司
• 主分類號: G01V3/38
• 分類號: G01V3/38
• 代理公司: 北京三高永信知識產權代理有限責任公司 11138
• 代理人: 賈敏

1.一種巖石儲集層的電阻率計算方法，其特征在于，所述方法包括：

建立目標儲集層樣品對應的網格圖，所述網格圖中包括多個網格；

根據所述目標儲集層樣品的孔隙度和所述網格圖中的網格總數，確定所述網格圖中的孔隙網格的數量；其中，所述孔隙網格是用于代表所述目標儲集層樣品中的孔隙的網格；

根據所述目標儲集層樣品的孔徑分布規律和所述孔隙網格的數量，確定所述目標儲集層樣品包含的k種孔徑各自對應的孔隙網格組中包含的網格數量，以及所述k種孔徑各自對應的孔隙網格組的數量；其中，所述孔徑分布規律用于表征所述k種孔徑的孔隙在所述目標儲集層樣品中的數量分布情況，所述孔隙網格組包括一個孔隙網格或者多個連接的孔隙網格，所述k為正整數；

根據所述k種孔徑各自對應的孔隙網格組中包含的網格數量，以及所述k種孔徑各自對應的孔隙網格組的數量，在所述網格圖中標注出所述k種孔徑各自對應的孔隙網格組；

根據標注后的所述網格圖中的各個節點對應的電壓值以及任意兩個相鄰的節點對應的電導率，計算所述目標儲集層樣品的電阻率；

所述根據所述k種孔徑各自對應的孔隙網格組中包含的網格數量，以及所述k種孔徑各自對應的孔隙網格組的數量，在所述網格圖中標注出所述k種孔徑各自對應的孔隙網格組，包括：

對于所述k種孔徑中的第i種孔徑，確定所述第i種孔徑在所述網格圖中的第一網格，當所述第一網格對應的值為第一預設值時，將所述第一網格對應的值設置為第二預設值，所述i為小于等于所述k的正整數；

把所述第一網格作為目標網格，判斷與所述目標網格共用邊長的相鄰網格中是否存在對應的值為所述第一預設值的第二網格，當所述目標網格的相鄰網格中存在所述第二網格時，將所述第二網格對應的值設置為所述第二預設值；

把所述第二網格作為所述目標網格，并再次從所述判斷與所述目標網格共用邊長的相鄰網格中是否存在對應的值為所述第一預設值的第二網格的步驟開始執行，直到在所述網格圖中存在m個連接的所述目標網格時，將所述m個連接的所述目標網格作為所述第i種孔徑對應的一個孔隙網格組；其中，所述m為所述第i種孔徑對應的孔隙網格組中包含的網格數量，所述m為正整數；

再次從所述確定所述第i種孔徑在所述網格圖中的第一網格的步驟開始執行，直到在所述網格圖中標注出所述第i種孔徑對應的n個孔隙網格組時，結束流程，所述n為正整數。

2.根據權利要求1所述的方法，其特征在于，所述根據所述目標儲集層樣品的孔隙度和所述網格圖中的網格總數，確定所述網格圖中的孔隙網格的數量，包括：

將所述目標儲集層樣品的孔隙度和所述網格圖中的網格總數相乘，得到所述網格圖中的孔隙網格的數量。

3.根據權利要求1所述的方法，其特征在于，所述孔徑分布規律包括所述k種孔徑各自對應的孔隙數量比值，其中，所述k種孔徑中的第i種孔徑對應的孔隙數量比值是指所述第i種孔徑在所述目標儲集層樣品中的孔隙數量與所述目標儲集層樣品的孔隙總數的比值，所述i為小于等于所述k的正整數；

所述根據所述目標儲集層樣品的孔徑分布規律和所述孔隙網格的數量，確定所述目標儲集層樣品包含的k種孔徑各自對應的孔隙網格組中包含的網格數量，以及所述k種孔徑各自對應的孔隙網格組的數量，包括：

對于所述第i種孔徑，根據所述第i種孔徑和基本孔徑的比值，得到所述第i種孔徑對應的孔隙網格組中包含的網格數量；其中，所述基本孔徑是指所述網格圖中的一個網格所代表的孔徑大小；

將所述第i種孔徑對應的孔隙數量比值和所述孔隙網格的數量相乘，得到所述第i種孔徑對應的孔隙網格組的數量。

4.根據權利要求1所述的方法，其特征在于，所述根據標注后的所述網格圖中的各個節點對應的電壓值以及任意兩個相鄰的節點對應的電導率，計算所述目標儲集層樣品的電阻率，包括：

根據高斯賽德爾迭代法和基爾霍夫電流定律，計算所述網格圖中的各個節點對應的電壓值；

根據所述網格圖中的各個節點對應的電壓值，計算所述網格圖中的各個網格對應的電勢差；

根據所述網格圖中的各個網格對應的電勢差以及任意兩個相鄰的節點對應的電導率，計算所述網格圖中的各個網格對應的電流值；

根據所述網格圖中的各個網格對應的電勢差以及所述網格圖中的各個網格對應的電流值，計算所述目標儲集層樣品的電阻率。