疏水缔合聚合物在孔隙介质中渗流阻力模拟分析.docVIP

疏水缔合聚合物在孔隙介质中渗流阻力模拟分析 　　摘要:对疏水缔合聚合物驱油机理及其在多孔介质流动中物化现象的进一步研究发现,该类聚合物渗流特征与普通高分子聚合物有较大不同。研究以收缩流道作为孔喉模型,采用渤海区块现场试验数据,使用数值计算方法,分析了缔合聚合物溶液流动阻力随流变参数和喉道直径的变化规律。并以此为理论基础,探讨其在多孔介质中的流动阻力特性,为疏水缔合聚合物驱的理论研究和现场分析提供依据。 　　关键词:疏水缔合聚合物;多孔介质;渗流阻力;数学模型;收缩流道 　　中图分类号:TE319文献标识码:A 　　 　　引 言 　　近年来,新型疏水缔合聚合物的研制成功,为高温高矿化度油藏的开采提供了技术保证。渤海油田某井进行的注聚矿场试验初步结果表明,疏水缔合聚合物驱将成为渤海稠油开采体系中的关键技术。对于聚合物溶液在孔隙介质中渗流特性的分析,主要针对线性聚合物溶液进行。大量室内研究表明,疏水缔合聚合物与传统高分子线性聚合物相比,无论是溶液结构还是渗流特性等方面都有明显不同,无法准确建立疏水缔合聚合物溶液在孔隙介质中渗流的数学模型,导致模拟和动态预测结果失真[1]。该研究结合物理模拟实验结果,应用数值模拟分析方法,建立疏水缔合聚合物溶液在孔隙介质中渗流阻力的数学模型。 　　1 模型建立 　　1.1 孔隙模型 　　油藏中孔隙类型一般较复杂,分为毛细管模型、网络模型、波纹管模型、扩张流道模型、扩张-收缩流道模型、带横向槽流道模型等[2]。该研究主要针对收缩流道孔喉模型。以驱替实验填砂管直管多孔介质模型作为物理实验模型。 　　1.2渗流数学模型 　　1.2.1 控制方程的建立 　　采用有限体积法,建立离散化数学模型。通过交错网格计算方法[3],分析不同控制体积的网格系统中压力p、x方向速度u、y方向速度v、应力T(包括剪切应力T??xy??、x方向的法向应力T??xx??、y方向法向应力T??yy??)之间的关系。 　　流动基本方程如下: 　　(1) 连续性方程(质量守恒): 　　??u??x+??v??y=0(1) 　　式中:u为x方向流速,m/s;v为y方向的流速m/s。 　　(2) 运动方程(动量守恒): 　　ρu??u??x+ρv??u??y=-??p??x+??T??xx????x+??T??xy????y 　　ρu??v??x+ρv??v??y=-??p??y+??T??xy????x+??T??yy????y(2) 　　式中:ρ为流体密度,?K/m??3;p为节点压力,Pa;T??xy??为剪切应力,Pa;T??xx??为x方向的法向应力,Pa;T??yy??为y方向法向应力,Pa。 　　因此,方程为: 　　T??xx??+τu??T??xx????x+v??T??xx????y-2??u??xT??xx??-2??u??yT??xy??=2η??u??x 　　T??xy??+τu??T??xy????x+v??T??xy????y-??u??yT??xy??-??v??xT??xx??=η??u??y+??v??x 　　T??yy??+τu??T??yy????x+v??T??yy????y-2??v??xT??xy??-2??v??yT??yy??=2η??v??y(3) 　　式中:τ为松弛时间,s;η为零剪切黏度,mPa#8226;s。 　　当τ=0时,上述控制方程会退化为描述牛顿流体(如常见的油田驱油剂水)在多孔介质中的渗流方程[4]。 　　 　　1.2.2 算法及实现 　　设x方向最后1个节点为(i,M1),倒数第2点为(i,M2),v为边界条件方程的基本变量,出口边界条件可近似表示为: 　　υ??i,M1??-υ??i,M2??δx=0(4) 　　得到: 　　υ??i,M1??=υ??*????i,M2??(5) 　　设υ??*????i,M2??为下一层次迭代计算的第1类边界条件。在此基础上建立方程组矩阵,设某个变量x的离散网格共有n个节点,则对n个节点建立的方程矩阵表达为: 　　a??11??,a??12??,…,a??1n?? 　　a??21??,a??22??,…,a??2n?? 　　…,…,…,… 　　…,…,…,… 　　a????n1??,a????n2??,…,a????nn??x??1x??2……x??n=s1s2……s??n(6) 　　式中:x??n(n=1,2,3,…,n)为变量x在离散网格中的分量。 　　再进行矩阵计算,按照SIMPLER算法(SIMPLE Revised)[5]修正和迭代,最终计算出各节点处p、u、v各值。 　　2 实验内容 　　通过文献调研和研究分析可知,影响流体在孔喉中流动阻力的