[发明专利]适用于油藏数值模拟的三维矢量渗透率的获取方法

权利要求书

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的计算方法。