[發明專利]一種毛細管力及相對滲透率曲線的計算方法和裝置

有效

• 申請號: 202110843918.2
• 申請日: 2021-07-26
• 公開（公告）號: CN113569407B
• 公開（公告）日: 2022-08-26
• 發明（設計）人: 王秀坤;劉峻嶸;唐維宇;宮潤東;方立勤
• 申請（專利權）人: 中國石油大學(北京)
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F119/14
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 賈磊;劉飛

1.一種毛細管力及相對滲透率曲線的計算方法，其特征在于，包括：

確定孔隙和喉道的截面均為四角星形的巖心孔隙網絡模型；

對所述巖心孔隙網絡模型進行單相流模擬計算，確定巖心的絕對滲透率；

對所述巖心孔隙網絡模型依次進行多種驅替類型的驅替模擬計算，分別確定各驅替模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系；

將各驅替模擬計算對應的所述油相有效滲透率與巖心含水飽和度的變化關系、所述水相有效滲透率與巖心含水飽和度的變化關系分別與所述巖心的絕對滲透率相結合，確定各驅替模擬計算對應的油相的相對滲透率與巖心含水飽和度的變化關系、水相的相對滲透率與巖心含水飽和度的變化關系。

2.根據權利要求1所述的毛細管力及相對滲透率曲線的計算方法，其特征在于，所述對所述巖心孔隙網絡模型進行單相流模擬計算，確定巖心的絕對滲透率，包括：

對所述巖心孔隙網絡模型進行單相流模擬計算，確定所述巖心孔隙網絡模型中每一孔隙的傳導系數，以及相互連通的孔隙之間的傳導系數；

根據所述相互連通的孔隙之間的傳導系數，確定每一孔隙的壓力；

根據每一所述孔隙的壓力及其對應的傳導系數，確定所述巖心的絕對滲透率。

3.根據權利要求2所述的毛細管力及相對滲透率曲線的計算方法，其特征在于，所述對所述巖心孔隙網絡模型依次進行多種驅替類型的驅替模擬計算，分別確定各驅替模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系，包括以下兩個步驟：

第一步驟：

對所述巖心孔隙網絡模型進行油驅水的模擬計算，確定油驅水的模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系；

第二步驟：

對所述巖心孔隙網絡模型進行水濕巖心水驅油的模擬計算，確定水濕巖心水驅油的模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系；

或，對所述巖心孔隙網絡模型進行油濕巖心水驅油的模擬計算，確定油濕巖心水驅油的模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系；

或，對所述巖心孔隙網絡模型進行混合潤濕巖心水驅油的模擬計算，確定混合潤濕巖心水驅油的模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系。

4.根據權利要求3所述的毛細管力及相對滲透率曲線的計算方法，其特征在于，所述對所述巖心孔隙網絡模型進行油驅水的模擬計算，確定油驅水的模擬計算對應的毛細管力與巖心含水飽和度的變化關系、油相有效滲透率與巖心含水飽和度的變化關系、水相有效滲透率與巖心含水飽和度的變化關系，包括：

確定每一喉道的門限壓力值、初始狀態下的第一巖心含水飽和度、初始狀態下的第一毛細管力、終止狀態下的第一毛細管力以及第一驅替步長；

將初始狀態下的第一毛細管力確定為當前第一毛細管力；

循環執行對所有喉道和孔隙進行篩選，確定其中門限壓力值小于或等于當前第一毛細管力的喉道，以及與該喉道連通的所有孔隙；利用篩選出來的所述喉道和所述孔隙的含水飽和度更新初始狀態下的第一巖心含水飽和度，確定當前第一毛細管力對應的當前第一巖心含水飽和度；根據所述當前第一巖心含水飽和度，以及每一孔隙或喉道分別對應的第一油相傳導系數和第一水相傳導系數，確定當前第一巖心含水飽和度對應的第一油相相對滲透率和第一水相相對滲透率；對所述當前第一毛細管力進行更新，將所述當前第一毛細管力與所述第一驅替步長之和確定為更新后的當前第一毛細管力，逐一確定當前第一毛細管力對應的當前第一巖心含水飽和度，以及當前第一巖心含水飽和度對應的第一油相相對滲透率和第一水相相對滲透率；

直至當前第一毛細管力大于所述終止狀態下的第一毛細管力時結束循環。