【技术实现步骤摘要】
本专利技术涉及油藏注气开发，特别涉及一种考虑油藏co2气驱渗流时变的数值模拟方法。
技术介绍
1、目前国内部分油藏已进入高含水阶段，亟需探究一种提高采收率的手段。注气开发技术是目前最有前景的提高采收率技术，co2进入地层后可有效补充地层能量，并且对于部分储层致密的非常规油藏，co2能进入微小孔隙，提高原油采收率。因此，利用二氧化碳驱油技术，可以达到增油与减排的双重目的。
2、目前，国内外关于油藏注气数值模拟技术的研究多基于储层物性恒定不变的情况，然而，在油藏注co2开发过程中，由于注入气的长期冲刷，储层的宏观、微观结构、非均质性、流体组成及流体分布较开发初期均发生一定变化，尤其影响渗流的主要渗流参数发生了改变，这些参数对剩余油的分布存在影响，致使储层参数、油藏流场规律变化更加复杂。因此，亟需一种方法能实现油藏co2气驱数值模拟过程中的渗流时变表征，保证模拟结果与矿场实际开发情况保持一致。
技术实现思路
1、针对上述问题，本专利技术旨在提供一种考虑油藏co2气驱渗流时变的数值模拟方法。
2、本专利技术的技术方案如下：
3、一种考虑油藏co2气驱渗流时变的数值模拟方法，包括以下步骤：
4、s1：获取粘度与残余液相饱和度的函数关系一以及粘度与残余气相相对渗透率的函数关系二；
5、s2：建立油藏co2气驱数值模拟模型，并对油藏进行网格化以及赋予各个网格单元初始地层油粘度；
6、s3：进行数值模拟，获得下一个时间步各网格
7、s4：根据所述函数关系一和所述函数关系二，确定所述下一个时间步各网格单元的残余液相饱和度和残余气相相对渗透率；
8、s5：根据步骤s4的计算结果更新所述油藏co2气驱数值模拟模型中同时间步网格单元的残余液相饱和度和残余气相相对渗透率；
9、s6：重复步骤s3-s5，直至数值模拟时间结束，输出co2气驱数值模拟的累计产油量。
10、作为优选，步骤s1具体包括以下子步骤：
11、s11：获取目标油藏的地层流体和岩心，并测得所述地层流体的原始粘度；
12、s12：对所述岩心进行地层水饱和，获得饱和地层水的岩心；
13、s13：采用所述地层流体驱替所述饱和地层水的岩心，建立束缚水饱和度；
14、s14：采用氮气进行驱替，获得油气相对渗透率曲线；
15、s15：对所述岩心重新饱和地层水，恢复初始状态；
16、s16：将所述地层流体与不同倍数的co2进行混合，获得混合地层流体，并测得所述混合地层流体的粘度；
17、s17：采用所述混合地层流体重复步骤s13-s15，获得不同倍数co2条件下的油气相对渗透率曲线；
18、s18：根据步骤s14和步骤s17的油气相对渗透率曲线测试结果，获得不同原油粘度条件下残余液相饱和度和残余气相相对渗透率；
19、s19：对步骤s18的结果进行回归分析，获得所述粘度与残余液相饱和度的函数关系一以及粘度与残余气相相对渗透率的函数关系二。
20、作为优选，步骤s1中，所述函数关系一为：
21、（1）
22、式中：sle为残余液相饱和度，%；a和b均为粘度与残余液相饱和度的相关系数，无量纲；μx为注入x倍数co2后混合地层流体的粘度，mpa·s。
23、作为优选，步骤s1中，所述函数关系二为：
24、（2）
25、式中：krge为残余气相相对渗透率，无量纲；c和d均为粘度与残余气相相对渗透率的相关系数，无量纲；μx为注入x倍数co2后混合地层流体的粘度，mpa·s。
26、作为优选，步骤s2中，所述油藏co2气驱数值模拟模型包括数值模拟运动方程和数值模拟渗流方程，所述数值模拟运动方程包括：
27、油相数值模拟运动方程：
28、（3）
29、式中：vo为油相的渗流速度，cm3/s；k为绝对渗透率，d；kro(sg)为含气饱和度sg条件下的油相相对渗透率，无量纲；μo为油相粘度，mpa·s；po为油相压力，10-1mpa；x为坐标轴上的x方向；ρo为油相密度，g/cm3；g为重力加速度，10-2cm/s2；α为流动方向与平面的夹角，°；
30、气相数值模拟运动方程：
31、（4）
32、式中：vg为气相的渗流速度，cm3/s；krg(sg)为含气饱和度sg条件下的气相相对渗透率，无量纲；μg为气相粘度，mpa·s；pg为气相压力，10-1mpa；ρg为气相密度，g/cm3；
33、残余气相数值模拟运动方程：
34、（5）
35、式中：vge为残余气相的渗流速度，cm3/s；krge为残余气相相对渗透率，无量纲；
36、所述数值模拟渗流方程包括：
37、油相数值模拟渗流方程：
38、（6）
39、式中：为拉普拉斯算子；kro(so)为含油饱和度so条件下的油相相对渗透率，无量纲；μo(p)为油藏压力p条件下的油相粘度，mpa·s；bo(p)为油藏压力p条件下的油相体积系数，m3/m3；p为梯度压力，mpa；为地层孔隙度，%；so为含油饱和度，%；t为时间，d；
40、气相数值模拟渗流方程：
41、（7）
42、式中：krg(so)为含油饱和度so条件下的气相相对渗透率，无量纲；μg(p)为油藏压力p条件下的气相粘度，mpa·s；bg(p)为油藏压力p条件下的气相体积系数，m3/m3；swc为束缚水饱和度，%；rs(p)为油藏压力p条件下的气体的溶解气油比，m3/m3；
43、残余气相数值模拟渗流方程：
44、（8）
45、式中：sle为残余液相饱和度，%。
46、本专利技术的有益效果是：
47、本专利技术能够在co2气驱数值模拟过程中实现残余气相相对渗透率的时变过程，更加真实模拟co2注入地层后的渗流特征。采用粘度作为桥梁来关联相对渗透率，通过利用残余气相相对渗透率变化和残余液相饱和度变化表征气驱时变过程，预测采收率的结果更加可靠合理，能够很好的满足海上油藏利用co2气驱数值模拟技术预测气驱采收率的需要，为海上油田气驱提供技术支持。
本文档来自技高网...【技术保护点】
1.一种考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，步骤S1具体包括以下子步骤：
3.根据权利要求1或2所述的考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，步骤S1中，所述函数关系一为：
4.根据权利要求1或2所述的考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，步骤S1中，所述函数关系二为：
5.根据权利要求1所述的考虑油藏CO2气驱渗流时变的数值模拟方法，其特征在于，步骤S2中，所述油藏CO2气驱数值模拟模型包括数值模拟运动方程和数值模拟渗流方程，所述数值模拟运动方程包括：
【技术特征摘要】
1.一种考虑油藏co2气驱渗流时变的数值模拟方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的考虑油藏co2气驱渗流时变的数值模拟方法，其特征在于，步骤s1具体包括以下子步骤：
3.根据权利要求1或2所述的考虑油藏co2气驱渗流时变的数值模拟方法，其特征在于，步骤s1中，所述函数关系一为：<...
【专利技术属性】
技术研发人员：何佑伟，张皓川，汤勇，刘梦云，安岭，秦佳正，谢坤，芮振华，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。