气水二相渗流数学模型.pptxVIP

1气、水二相裂隙（孔隙）流;式中;;上述方程（1）和方程（2）中的Pw、Pg、Sw、Sg分别满足下面两个附加方程 ;将方程（3）、（4）、（5）、（6）耦合起来，这样四个方程中只有4个未知数Pw、Pg、Sw、Sg，未知数个数和方程个数相等结合初边值条件，即可用数值法进行求解。 方程（3）和（4）中的气、水源汇项qg、qw，如给定气、水产量，则直接进入方程参与计算，如给定井底流动压力（或井中动液面位置），则在计算时段内用稳定流方法进行处理，在气井数值模拟研究中，通常给定井底流动压力（井中动液面位置）。 ;;式中; Pwfg 、Pwfw分别为气和水的井底流压，MPa；如给定井中的动液面位置，则需将其换算成井底流压:Pwfw =ρwgΔh+Patm , Δh为井中动液面至含水层中心点距离，Patm为大气压；Pwfg=Pc+ Pwfw。 Pg、Pw分别为对应等效供给半径re的气压和水压，MPa; 其余符号意义同前。 ;方程（7）和（8）中的等效供给半径re在有限差分计算中可用下式计算： ;2定解条件 ;3模型求解 ;（1）借助水气饱和度关系式,通过乘以适当系数，合并气方程和水方程，消去差分方程中的变量和,从而得到一个只含变量和的压力方程； （2）由毛管压力公式,可将表示为和的形式，代入上述方程，得到一个只含变量的压力差分方程； （3）将井点所在网格的水气产量项均作显式处理，由上一时间阶段的饱和度值计算出井点网格的水气产量； ;（4）方程左端达西系数用上一时间阶段值，同时毛管压力也用上一时间阶段值，即显式处理系数，由此形成一线性代数方程组，即可用迭代法或直接法进行求解； （5）在解出压力差分方程式之后，将解出的压力值作为已知条件代入气方程或水方程，用显式计算格式求出饱和度值。 ;3.2 全隐式解法 ;3.2.1基本原理 ;;当k=0时;当 时;对于非线性方程组:;按上述原理处理，就可以得到一组用牛顿迭代法求解的方程形式： ;3.2.2方程的建立 ;式中:;2.方程（3）右端项的展开 ;将式(14)代入式(13)得： ;将式(15)代入式(12)得 ;将方程(10)、(11)、(16)代入方程(3)得 ;再将流动项进行空间差分;;;;将式(18) 、(19)、 (20)、(21) 、(22)代入式(17)得下列方程： ;;;b.考虑方程（4）式，将代入方程(4)式得：;;;;3.2.3参数及边界条件的处理; 2)渗透率 由裂隙（孔隙）孔隙度的变化引起的岩土体渗透率的变化，可用下式表示： ;3) 流体密度 由于抽水降压所引起的气和水的密度变化，可用下式表示 ： ;4)流动系数 在水相、气相方程中，流动系数，中，渗透率是空间坐标的函数，相邻两块单元之间的渗透率取调和平均值，可用下式表示： ;， 的值取流动方向上的上游权衡,可用下式表示： ;;5)相对渗透率与毛管压力 相对渗透率与毛管压力是根据相对渗透率和饱和度的关系曲线，毛管压力与饱和度的曲线，由线性插值给出，可用下式表示 : ;b.边界条件的处理 1)内边界条件的处理 网格（i，j，k）上的气气井或水井，如已知产量或则可在该网格的方程中直接加上产量项；当已知井底压力时，通过等效供给半径，把用网格压力和井底流动压力表示，然后进入该网格的方程之中。 ;式中等效供给半径 用下式计算 ;2)外边界条件处理 : 本次研究差分网格采用块中心接点，外边界条件的处理是通过在边界网格处虚拟一排网格来实现，当边界为已知压力和饱和度时，采用线性插值法处理，可用下式表示： ;;式中:;当边界为封闭边界时，令边界网格和虚拟网格的流体势及饱的度值相等，可用???式表示： ;式中： ;3.2.4方程组的求解 ;经处理的线性代数方程组的系数矩阵是一个七对角系数矩阵，且每一对角线均有 和 两类待求未知数，因此，直接求解不但计算时间长，消耗计算机内存大，且计算复杂、难度大，有时计算根本无法进行。为此本次研究采用交替方向取隐式差分法，将每一个计算时间步分成相等的三个部分，先对x方向进行扫描，x方向取隐式，y、z方向取显式，将方程(21)和(23)处理成下列形式： ;;;;;;采用同样的方法，可求得k+2/3次的 和 的值。同理，再对z方向进行扫描，z方向取隐式，x，y方向取显式，将方程（42）和方程（43）处理成下列形式： ;;用同样的方法，可求得k+1次的 和 值，再由下列迭代过程计算出 和 的值，即： ;重复上述计算过程，直到 足够小，满足精度要求为止，此时的 的值即为所求时刻 的值，再由式（33）分别求得本