[发明专利]一种油水两相相对渗透率的模拟方法

权利要求书

1.一种油水两相相对渗透率的模拟方法，其包括以下步骤：

获取岩样的孔喉数据，所述孔喉数据包括孔隙半径的频率分布数据和喉道半径的频率分布数据；

基于所述孔喉数据，构建孔隙网络；

在所述孔隙网络中建立油水两相流体流动模拟准则；

根据所述两相流体流动模拟准则，在所述孔隙网络中实施两相流体的流动模拟，以获取不同孔隙网络含水饱和度下的孔隙网络水相相对渗透率数据和孔隙网络油相相对渗透率数据；

其中，在所述孔隙网络中建立油水两相流体流动模拟准则包括以下过程：

获取所述孔隙网络在驱替过程中的驱替前缘，从与所述驱替前缘相连的所有喉道中选出净压力最大的喉道作为最优喉道；

将所述最优喉道，以及与所述最优喉道相连的所有孔隙的当前状态更新为被非润湿相充满的状态，同时，对所述孔隙网络的压力场的当前状态进行更新。

2.根据权利要求1所述的方法，其中，基于所述孔喉数据，构建孔隙网络包括以下过程：

基于所述孔隙半径的频率分布数据，对所述孔隙网络中的孔隙进行赋值处理；

基于所述喉道半径的频率分布数据，对所述孔隙网络中的喉道进行赋值处理。

3.根据权利要求2所述的方法，其中，基于所述孔隙半径的频率分布数据，对所述孔隙网络中的孔隙进行赋值处理包括以下过程：

过程1：基于所述孔隙半径的频率分布数据，构建孔隙半径的频率分布直方图；

过程2：以所述孔隙半径的频率分布直方图为基准，随机生成一个落在所述直方图中孔隙半径分布值内的数据，将当前随机生成的数据赋给相应的孔隙作为其半径值，并统计已赋值孔隙中各种半径值出现的次数，一旦当前生成的随机半径值在已赋值孔隙中出现的次数超过所述相应的孔隙的半径频数，则确定当前生成的随机半径值无效，重新生成一个新的随机半径值，直至不再出现当前生成的随机半径值在之前已赋值孔隙中出现的次数超过所述相应的孔隙的半径频数这一现象。

4.根据权利要求2所述的方法，其中，基于所述喉道半径的频率分布数据，对所述孔隙网络中的喉道进行赋值处理包括以下过程：

过程1：基于所述喉道半径的频率分布数据，构建喉道半径的频率分布直方图；

过程2：以所述喉道半径的频率分布直方图为基准，随机生成一个落在所述直方图中喉道半径分布值内的数据，将当前随机生成的数据赋给相应的喉道作为其半径值，并统计已赋值喉道中各种半径值出现的次数，一旦当前生成的随机半径值在之前已赋值喉道中出现的次数超过所述相应的喉道的半径频数，则确定当前生成的随机半径值无效，重新生成一个新的随机半径值，直至不再出现当前生成的随机半径值在之前已赋值喉道中出现的次数超过所述相应的喉道的半径频数这一现象。

5.根据权利要求1所述的方法，其中，对所述孔隙网络的压力场的当前状态进行更新包括以下过程：

将各孔隙的流量之和的方程联立成式1所示的矩阵方程组，并对所述矩阵方程组进行求解，即得到新的压力场

在式1所示的方程组中，且其中，N_i表示与孔隙i通过其周围喉道直接相连的孔隙集，r_ij表示孔隙i和孔隙j之间喉道的半径，L_ij表示孔隙i和孔隙j之间喉道的长度，μ_ij表示孔隙i和孔隙j之间喉道的粘度，P_n表示孔隙n的毛管力，b_n表示孔隙n的通过的流量。

6.根据权利要求1或5所述的方法，其中，获取所述孔隙网络在驱替过程中的驱替前缘包括以下过程：

获取所述孔隙网络在驱替过程中的所有被非润湿相充满的孔隙，根据这些孔隙即可获取所述孔隙网络在驱替过程中的驱替前缘。