(6.1)--66-1多组分渗流理论 - 副本.ppt

引言注气等开发过程中涉及的多相多分组系统是一个很复杂的物理化学系统，因此无论在对系统本身的物理化学性质的研究还是对于流动规律的研究，包括对物理化学过程的描述和流动规律的描述，都遇到极为困难的问题。即使有可能建立起基本微分方程，其求解也是相当困难的，往往不得不借助计算机求出其某些具体条件下的数值解。在本部分中我们仅就多相多组分渗流的基本规律、数学模型建立以及其求解方法作一一叙述。一、油气相平衡方程及闪蒸计算方法二、常用的状态方程三、求解PT相图中泡点压力方程的建立及求解方法四、求解PT相图中露点压力方程的建立及求解方法五、多相多组分流体物性参数计算方法六，多相多组分渗流数学模型及求解方法本章内容提要一、油气相平衡方程及闪蒸计算方法求解问题：已知一摩尔的流体混合物，其组成为zj（j=1、2、3、…、Nc），Nc为流体混合物的组分数，设体系温度为T，压力为P，达到流体相平衡，试建立描述这一问题的数学模型，并简述其求解方法。求解参数：流体相平衡时油相摩尔量L和气相摩尔量V，油相组成为xj以及气相组成为yj（j=1、2、3、…、Nc）,由此可见，求解的未知量为2Nc+2个。1、总的物质平衡方程：2、组分i的物质平衡方程：3、组分的约束方程：4、组分i的相态平衡方程：—i组分液相逸度—i组分气相逸度(1)(2)(3)(4)（一）数学模型—i组分液相逸度系数—i组分气相逸度系数1、L-x迭代法：这种方法适用于当V0.5，即混合物中气相占优这种情况。此时求解的未知数则为xi和L。2、V-y迭代法：这种方法适用于当L0.5，即混合物中液相占优这种情况。此时求解的未知数则为yi和V。思考：L-x迭代法、V-y迭代法—为什么？？？（二）迭代求解方法牛顿迭代法是解非线性代数方程组方法之一，它是同类方法中收敛速度最快的方法（平方级收敛）。设有n个函数f1,f2,…fn,n个未知量x1,x2,…xn求解的非线性代数方程组为：求解方法是：(1)给初始近似值：(2)计算:（二）迭代求解方法—牛顿迭代方法(3)计算雅可比矩阵：其中：(5)计算：(4)求解线性代数方程组：(6)用代替重复以上过程直至收敛时为止。将方程(1)带入(2)：(5)对方程(5)求导：当给定混合物组成时，zi=常数，即zi与xi和L无关，因此有：（三）迭代求解过程方程(4)和方程(3)还可以写成以下余量形式：(7)由(6)和(7)可写出Nc+1个方程组成的方程组。利用Newton-Raphson方法求解。(6)Newton-Raphson方法求解要点是形成Jacobi矩阵元素：下面推导矩阵Jacobi元素aij计算公式：另外，则由逸度系数公式求得。求解方程为：其中：由此Jacobi矩阵形式为：求解：首先对行向量消元，将矩阵转化为三角形矩阵，然后回代求解，得到xi和L的增量，再以新的迭代步的值代入，重复这一过程直到收敛时为止。收敛原则：(ε一般取为10-6)在求得xi和L以后，可以求得: