详情

原文传递一种混合边界条件渗吸作用湿相相对渗透率的计算方法

专利名称：	一种混合边界条件渗吸作用湿相相对渗透率的计算方法
摘要：	本发明涉及石油开发技术领域，具体公开了一种混合边界条件渗吸作用湿相相对渗透率的计算方法，主要包括以下步骤：S1确定一维混合边界渗吸时间与渗吸采出程度的函数关系、S2得到渗吸时间对渗吸采出程度的微分与渗吸采出程度的函数关系、S3得到渗吸时间对渗吸采出程度的微分与渗吸采出程度函数关系直线段的斜率和截距、S4得到E和D的表达式、S5得到混合边界条件渗吸作用湿相相对渗透率的表达式、S6计算混合边界条件渗吸作用的湿相相对渗透率值。本发明所述的一种混合边界条件渗吸作用湿相相对渗透率的计算方法，可以根据一维混合边界渗吸实验结果快速的求解渗吸作用湿相相对渗透率，无需经过多次假设试算，并且简化了混合边界条件渗吸作用湿相相对渗透率的计算步骤。
专利类型：	发明专利
国家地区组织代码：	湖北;42
申请人：	中国地质大学(武汉)
发明人：	孟庆帮;蔡建超;孙秋莹
专利状态：	有效
申请日期：	2019-06-26T00:00:00+0800
发布日期：	2019-11-01T00:00:00+0800
申请号：	CN201910562144.9
公开号：	CN110398450A
代理机构：	武汉知产时代知识产权代理有限公司
代理人：	龚春来
分类号：	G01N15/08(2006.01);G;G01;G01N;G01N15
申请人地址：	430000 湖北省武汉市洪山区鲁磨路388号
主权项：	1.一种混合边界条件渗吸作用湿相相对渗透率的计算方法，其特征在于，主要包括以下步骤： S1、基于一维正向渗吸前缘推进距离与渗吸时间的函数关系，确定一维混合边界渗吸时间与渗吸采出程度的函数关系，其表达式为：式中，R为渗吸采出程度，Swf为渗吸前缘平均含水饱和度； S2、基于S1中得到的一维混合边界渗吸时间与渗吸采出程度的函数关系,并对渗吸采出程度进行求导，得到渗吸时间对渗吸采出程度的微分与渗吸采出程度的函数关系，其表达式为： S3、将S2中得到的渗吸时间对渗吸采出程度的微分与渗吸采出程度的函数关系采用线性拟合方法，得到其函数直线段的斜率k和截距b，其表达式分别如下： S4、基于S3中得到的斜率k和截距b，即可分别得到E和D的表达式，其分别如下所示： D＝b·E·Swf/2 S5、基于S4中得到的E和D的表达式，以及S1中所述一维正向渗吸前缘推进距离与渗吸时间的函数关系，即可得到混合边界条件渗吸作用湿相相对渗透率，其表达式如下：式中，μw为湿相粘度，mPa·s；μnw为非湿相粘度，mPa·s；krw为湿相相对渗透率；krnw为非湿相相对渗透率，其值为1； S6、获取工区的生产数据，代入S5中得到的混合边界条件渗吸作用湿相相对渗透率的表达式中，以计算该工区混合边界条件渗吸作用的湿相相对渗透率值。 2.根据权利要求1所述的一种混合边界条件渗吸作用湿相相对渗透率的计算方法，其特征在于，S1中确定一维混合边界渗吸时间与渗吸采出程度的函数关系，还包括以下步骤： S101、确定所述一维正向渗吸前缘推进距离与渗吸时间的函数关系，其表达式如下：式中，xf为渗吸前缘推进距离，cm；L为岩心长度，cm；K为岩心渗透率，×10-3μm2；φ为岩心孔隙度；μw为湿相粘度，mPa·s；μnw为非湿相粘度，mPa·s；pc,f为渗吸前缘等效毛管压力，Pa；krw为湿相相对渗透率；krnw为非湿相相对渗透率，其值为1；Swf为渗吸前缘平均含水饱和度。 S102、确定所述渗吸采出程度的表达式，其表达式如下：式中，R为渗吸采出程度，Swi为初始含水饱和度； S103、假设S102中Swi为零，即可进一步得到所述渗吸采出程度的表达式，其表达式如下： S104、结合S101中所述一维正向渗吸前缘推进距离与渗吸时间的表达式(1)和S103中所述渗吸采出程度的表达式，即可得到所述一维混合边界渗吸采出程度与渗吸时间的函数关系，其表达式如下： 3.根据权利要求2所述的一种混合边界条件渗吸作用湿相相对渗透率的计算方法，其特征在于，S5中确认所述混合边界条件渗吸作用湿相相对渗透率的表达式，还包括以下步骤： S501、根据S4中得到的E和D的表达式，以及结合S101中所述D的表达式(2a)，即可得到所述混合边界条件渗吸作用湿相相对渗透率的表达式，其表达式如下： 4.根据权利要求1所述的一种混合边界条件渗吸作用湿相相对渗透率的计算方法，其特征在于，S6中所述生产数据包括：湿相粘度、非湿相粘度、岩心长度和渗吸前缘平均含水饱和度。
所属类别：	发明专利