[發明專利]適用于油藏數值模擬的三維矢量滲透率的獲取方法

有效

• 申請號: 201810356830.6
• 申請日: 2018-04-19
• 公開（公告）號: CN108804750B
• 公開（公告）日: 2021-11-12
• 發明（設計）人: 蘇海波;孫業恒;楊勇;曹緒龍;張世明;于金彪;汪勇;董亞娟;張波;宋勇;易紅霞;侯玉培
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司勘探開發研究院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06T17/00;G06F111/10
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 崔曉艷

1.適用于油藏數值模擬的三維矢量滲透率的獲取方法，其特征在于，該適用于油藏數值模擬的三維矢量滲透率的獲取方法包括：

步驟1，計算相鄰節點單元的滲流面積、滲流長度和確定驅替壓力梯度的方向；

步驟2，計算相鄰節點單元滲流面的法向量與水平X軸、Y軸和垂向Z軸之間的夾角余弦值；

步驟3，計算滲流驅替方向與相鄰節點單元滲流面法向量之間的夾角余弦值；

步驟4，根據節點單元的水平、垂向滲透率值和矢量滲透率模型，計算相鄰節點單元間在驅替壓力梯度方向上的矢量滲透率數值；

步驟1包括：

a，節點單元I和節點單元J是空間上相鄰的兩個節點單元，S_ABCD是節點單元I和節點單元J的共有滲流面積，計算節點單元I的體中心坐標I(xi,yi,zi)；

b，計算滲流面ABCD面中心O的空間坐標O(x_o,y_o,z_o)；

c，計算滲流距離L_IO，并根據滲流距離L_IO的計算方法計算滲流距離L_JO；

在步驟1中，ΔP_i為節點單元I的壓力，ΔP_j為節點單元J的壓力，當ΔP_i＞ΔP_j時，其驅替壓力梯度方向為從節點單元I到節點單元J，當ΔP_i＜ΔP_j時，其驅替壓力梯度方向為從節點單元J到節點單元I；

在步驟4中，獲取節點單元I在X方向、Y方向和Z方向的滲透率K_xi、K_yi、K_zi；獲取節點單元J在X方向、Y方向和Z方向的滲透率K_xj、K_yj、K_zj；計算節點單元I與相鄰節點單元J在驅替方向上的三維矢量滲透率；

當ΔP_i＞ΔP_j時，其驅替壓力梯度方向為從節點單元I到相鄰節點單元J，沿該驅替方向上的三維矢量滲透率K_pi，計算公式為

K_pi＝(K_xicos²α+K_yicos²β+K_zicos²ψ)cosθ (8)

當ΔP_i＜ΔP_j時，從節點單元J到相鄰節點單元I沿驅替壓力梯度方向的三維矢量滲透率K_pj，計算公式為

K_pj＝(K_xjcos²α+K_yjcos²β+K_zjcos²ψ)cosθ (9)

式中，K_xi、K_xj為分別為節點單元I和J水平方向X軸的滲透率，md；K_yi、K_yj分別為節點單元I和J在水平方向Y軸的滲透率，md；K_zj、K_zj分別為節點單元I和J在垂直方向Z軸的滲透率，md；α、β、ψ分別為節點單元I與相鄰節點單元J的共有滲流面的法向量與水平X軸、Y軸和垂向Z軸之間的夾角，弧度；θ為驅替壓力梯度方向與該滲流面法向量間的夾角。

2.根據權利要求1所述的適用于油藏數值模擬的三維矢量滲透率的獲取方法，其特征在于，在步驟a中，空間點A、B、C、D、E、F、H、G是構成節點單元I的空間坐標點，節點單元I的體中心坐標xi的計算公式為

依據節點單元I的體中心坐標xi的計算方法計算坐標點yi和zi的值。

3.根據權利要求1所述的適用于油藏數值模擬的三維矢量滲透率的獲取方法，其特征在于，在步驟b中，空間坐標O(x_o,y_o,z_o)的計算方法同節點單元I的體中心坐標xi的計算方法。