[發明專利]一種離散裂縫中基質與裂縫間的傳質計算方法

有效

• 申請號: 201810287748.2
• 申請日: 2018-04-03
• 公開（公告）號: CN108595782B
• 公開（公告）日: 2020-09-04
• 發明（設計）人: 曹仁義;方思東;程林松;賈品;黃世軍;薛永超
• 申請（專利權）人: 中國石油大學(北京)
• 主分類號: G06F30/20
• 分類號: G06F30/20
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 王濤

1.一種離散裂縫中基質與裂縫間的傳質計算方法，其特征在于，包括：

獲取滲流區域的包含基質和裂縫的油藏信息，并生成用于表示所述基質的背景網格和用于表示所述裂縫的嵌入式單元，并確定背景網格和嵌入式單元中基質與裂縫之間的幾何信息，包括：將油藏剖分為若干用于表示所述基質的背景網格，并為所述背景網格的各節點和中心點進行編號，確定每個背景網格的局部節點編號對應關系；在背景網格中生成用于表示裂縫的嵌入式單元，并判斷每條裂縫與所述背景網格是否相交，獲得所述裂縫與所述背景網格相交的節點；在所述背景網格中判斷是否存在裂縫相交，若存在裂縫相交，則以裂縫相交的交點劃分裂縫單元；

根據所述幾何信息，采用Peaceman推導數值模擬中的井指數方法，確定背景網格下的基質和裂縫竄流的傳導系數關系，包括：

根據所述背景網格的各節點、中心點、所述裂縫與所述背景網格相交的節點以及所述裂縫單元，采用Peaceman推導數值模擬中的井指數方法，確定背景網格下的基質和裂縫竄流的傳導系數關系；所述傳導系數關系為裂縫與基質竄流量和裂縫與基質壓差的關系：

其中，k為地層滲透率；h為儲層厚度；μ_o為油相粘度；B_o為油相地層體積系數；p_ed為網格邊界的壓力；p_f為裂縫單元壓力；q_omf為基質到裂縫中油相的竄流量；n是網格邊界單位法向向量，nf是背景網格內裂縫單元數；G為2D拉普拉斯方程基本解；[G_fm]、[G_ff]、[Gf]和為與G相關的矩陣；

根據預先設置的流量和壓力連續性條件，將基質和裂縫的有限體積流動解和基于所述傳導系數關系的基質和裂縫竄流的邊界元解進行耦合，并進行流動模擬，生成流動模擬數據，包括：

將基質和裂縫的有限體積流動解和基于所述傳導系數關系的基質和裂縫竄流的邊界元解進行耦合，獲得油組分竄流量表達式：

獲得水組分竄流量表達式：

獲得氣組分竄流量表達式：

其中，q_wmf為基質到裂縫中水相的竄流量；q_gmf為基質到裂縫中氣相的竄流量；k_ro為油相相對滲透率；k_rw為水相相對滲透率；k_rg為氣相相對滲透率；μ_w為水相粘度；B_w為水相地層體積系數；μ_g為氣相粘度；B_g為氣相地層體積系數；p_o為基質中油相壓力；p_w為基質中水相壓力；p_g為基質中氣相壓力；p_of為裂縫中油相壓力；p_wf為裂縫中水相壓力；p_gf為裂縫中氣相壓力；R_s為溶解氣油比；[G_conv]為與G相關的矩陣；

根據預先設置的流量和壓力連續性條件，對所述油組分竄流量表達式、水組分竄流量表達式和氣組分竄流量表達式進行求解，并進行流動模擬，生成流動模擬數據；

根據所述流動模擬數據，實時獲取并顯示壓裂面縫和復雜裂縫網格條件下地層內可壓縮流體的不穩定滲流的滲流場分布數據。