黏弹性聚合物溶液流动特性的数值模拟研究.docxVIP

黏弹性聚合物溶液流动特性的数值模拟研究 黏弹性流体主要驱油效果 在石油工业中，存在许多混合弹性流的问题，尤其是原油和聚合物溶液通过油藏孔道的流动，这一直是混合弹性流特性研究的重点。黏弹性流体可以改变油藏孔隙中残余油的流场和受力状况,大量岩心实验和微观实验证实,黏弹性聚合物(如水解聚丙烯酞胺等)可以较纯黏聚合物(如黄原胶等)驱替出更多的残余油,但目前对此现象尚无一致理论解释:文献认为,黏弹性聚合物溶液在孔喉处引起的高压降和高流速加大了其对油块、油滴或油膜的驱替和冲洗作用,有助于提高驱油效率;文献通过进行水、甘油、聚丙烯酞胺溶液驱替各类残余油对比实验,发现黏弹性流体主要以剥离、拉拽、携带等方式提高微观驱替效率;文献以T型流道模拟油藏孔喉和盲端,发现黏弹性聚合物在该类变截面孔隙中能够形成较强的黏弹涡,黏弹涡的形成和发展可使部分残余油成为可驱动原油,从而提高聚合物波及效率。 探索黏弹性流体在油藏孔隙中的微观流动特性有助于深化聚合物溶液微观驱油机理研究。为了进一步从理论上说明黏弹性流体在孔隙中的流动特征,利用黏弹性流体的上随体Maxwell本构方程描述油藏条件下以第一法向应力差为主要特征的聚合物溶液的流变性。但现有计算方法只能得到Re数和We数较小(Re≤8,We≤0.8)时的结果,因此,利用流体有限元分析软件Polyflow对收缩流道内的黏弹性聚合物溶液的流动过程进行了等温数值模拟分析,得到收缩流道中Re数和We数较大时的流动结果,确定了流体黏弹性对流场和应力场的影响规律。 1 规则流道模型 黏弹性聚合物溶液真实流道类型较多,但多数情况下可以简化为收缩的规则流道模型,如图1所示。通过Gambit2.2.30建立收缩流道的几何模型(图2),可生成Polyflow前处理器文件。 2 流体流动特性 在拟定流场时,考虑收缩段工艺条件及聚合物特性,假设流体为黏弹性不可压缩流体,等温、定常流动,不计重力,并采用上随体Maxwell本构方程来描述第一法向应力差为主要特征的黏弹性聚合物溶液流变性,则流动基本方程如下: 连续方程为: ?v=0 (1) 运动方程为: ρV·?v=-?p+?T(2) 本构方程为: Tij+λδδtTij=η0Aij(3)Τij+λδδtΤij=η0Aij(3) 式中:ρ为流体密度,kg/m3;v为流体流速,m/s;V为流体体积,m3;p为流体受到压力,Pa;T为偏应力张量,Pa;λ为松弛时间,s;t为流体流动时间,s;η0为临界剪切粘度,mPa·s;A为变形速度张量,Pa;i,j为空间向量基。 应力张量上随体导数则为: δδtTij=??tTij+vm?vm?xm??vm?xiTmj??vm?xjTim(4)δδtΤij=??tΤij+vm?vm?xm-?vm?xiΤmj-?vm?xjΤim(4) 式中:vm、vi、vj均为v的空间物理分量,m/s;Tij、Tmj、Tim均为T的空间物理分量。 引入无量纲参数韦森博格数(We),定义其为第一法向应力N1与2倍黏性剪切应力τ的比值,即:We=N1/(2τ) (5) 且: N1=2ληγ2(6) τ=ηγ(7) 式中:We为韦森博格数;N1为第一法向应力,N;τ为黏性剪切应力,N;η为剪切粘度,mPa·s;γ为剪切速率,s-1。 式(5)反应了弹性和黏性对黏弹性流体流动状态的影响:当We很大时,流体流动特征主要由第一法向应力差决定,即流体弹性起主要作用;当We很小时,流动特征则主要由流体黏性力决定,即流体黏性起主要作用。 将式(6)、(7)代入式(5),可得We的计算公式: We=λU/L(8) 式中:U为流体特征速度,m/s;L为流体特征长度。 而对于黏弹性聚合物在次流道内的流动,其特征速度可以取作流体在流道中的平均速度,且雷诺数为: Re=ρUL/η(9) 3 处流动的poaseiache作为 采用有限元法,假设流体在流道壁面处速度为零,流道壁壁面无滑移条件,且入口处流动是充分发展的Poiseuille流动。利用Polyflow软件模拟黏弹性聚合物溶液通过收缩流道,通过Polyflow自带的后处理模块Flpost进行图像处理,可以得到黏弹性对流道中聚合物流速和应力影响的可视化结果。 3.1 不同雷诺数时流体流态分析 通过数值模拟,可以得到相同Re数条件下,z坐标为零时,牛顿流体(We=0)和黏弹性聚合物(以We=4为例)在收缩断面处速度等值线图(图3)以及牛顿流体(We=0)和黏弹性聚合物(以We=4为例)在流道中线处的速度图(图4)。 由图3可以看出,当雷诺数较大(Re=13.8)时,黏弹性对流体流态具有明显影响,靠近中心线处黏弹性对流体流速的影响较大,而靠近壁面处黏弹性对流体流速的影响较小。由图4还可看出,当流道半径收缩时,牛顿流体流速增大到一定值