[发明专利]基于储层物性的油水相对渗透率表征处理方法

权利要求书

1.基于储层物性的油水相对渗透率表征处理方法，其特征在于，该基于储层物性的油水相对渗透率表征处理方法包括：

步骤1，选取能反映区块或油藏物性特征的多条实验室测得的油水相对渗透率曲线数据；

步骤2，建立油水相对渗透率曲线特征参数与岩心样品的空气渗透率K和孔隙度φ的函数关系式；

步骤3，对含水饱和度S_w、含油饱和度S_o、水相相对渗透率K_rw和油相相对渗透率K_ro进行标准化处理；

步骤4，建立标准化后的水相相对渗透率K_rwn与标准化后的含水饱和度S_wn以及标准化后的油相相对渗透率K_ron与标准化后的含油饱和度S_on之间的幂函数关系表达式，并建立幂指数A，B与空气渗透率K和孔隙度φ的函数关系式；

步骤5，利用不同的储层物性得到不同的油水相对渗透率曲线；

在步骤2中，建立的油水相对渗透率曲线特征参数与岩心样品的空气渗透率K和孔隙度φ的函数关系式为：

S_wi＝a₁ln K+b₁φ+c₁ (1)

S_ow＝a₂ln K+b₂φ+c₂ (2)

K_rwmax＝a₃ln K+b₃φ+c₃ (3)

ln K_wmax＝a₄ln K+c₄ (4)

ln K_omax＝a₅ln K+c₅ (5)

式中：

S_wi—束缚水饱和度，无因次；

S_ow—油水共渗饱和度区间，无因次；

K_rwmax—残余油饱和度下的水相相对渗透率，无因次；

K_wmax—残余油饱和度下的水相渗透率，mD；

K_omax—束缚水饱和度下的油相渗透率，mD；

K—岩心样品空气渗透率，mD；

φ—岩心样品孔隙度，％；

a_i、b_i、c_i(i＝1，2，……，5)—待定系数；

在此基础上，利用岩心样品数据进行线性回归，求得回归系数，建立以上回归关系式；

在步骤3中，对含水饱和度S_w、含油饱和度S_o、水相相对渗透率K_rw和油相相对渗透率K_ro进行标准化处理，其公式如下：

K_rwn＝K_rw/K_rwmax (9)

K_ron＝K_ro/K_romax (10)

式中：

S_w—含水饱和度，无因次；

S_wn—标准化后的含水饱和度，无因次；

S_wi—束缚水饱和度，无因次；

S_or—残余油饱和度，无因次；

S_on—标准化后的含油饱和度，无因次；

K_rw—水相相对渗透率，无因次；

K_rwmax—残余油饱和度下的水相相对渗透率，无因次；

K_ro—油相相对渗透率，无因次；

K_romax—束缚水饱和度下的油相相对渗透率，无因次；

K_rwn—标准化后的水相相对渗透率，无因次；

K_ron—标准化后的油相相对渗透率，无因次；

在步骤4中，首先建立标准化后的水相相对渗透率K_rwn与标准化后的含水饱和度S_wn以及标准化后的油相相对渗透率K_ron与标准化后的含油饱和度S_on之间的幂函数关系表达式，其表达式如下：

K_rwn＝S_wn^A (11)

K_ron＝S_on^B (12)

式中：

A、B—待定幂函数指数；

对式(11)和式(12)两边取自然对数可得：

ln K_rwn＝A ln S_wn (13)

ln K_ron＝B ln S_on (14)

对式(13)和式(14)进行线性回归，求得系数A、B；

通过实验数据分析发现幂函数指数A，B与空气渗透率K和孔隙度φ存在以下函数关系：

A＝a₆ln K+b₆φ+c₆ (15)

B＝a₇ln K+b₇φ+c₇ (16)

a_i、b_i、c_i(i＝6,7)—待定系数；

根据式(15)和式(16)对不同物性储层油水相对渗透率曲线数据的表征预测处理；

在此基础上，利用岩心样品数据进行线性回归，求得回归系数，建立以上回归关系式。