A one-dimensional moving grid solution for the coupled non-linear equations governing multiphase flow in porous media. 2: Example simulations and sensitivity analysis

This paper presents numerical examples for the moving grid finite element algorithm derived in Part Ito solve the non-linear coupled set of PDEs governing immiscible multiphase flow in porous media in one dimension. Examples include single- and double-front simulations for two- and three-phase flow regimes and incorporating a mass sink. The modelling approach is shown to achieve significant savings in computation time and memory allocation when compared with fixed grid solutions of equivalent accuracy. This work includes sensitivity analyses for the parameters which are incorporated in the grid adaptation method, including the curvature weights, artificial viscosity and artificial repulsive force. It is found that the curvature weights are exponential functions of the negative ratio of the square root of the domain length to the number of discrete nodes. These weighting parameters are also shown to depend upon the shape of the front. On the basis of the examined simulations, it is recommended that artificial viscosity be neglected in the solution of the coupled non-linear set of PDEs governing multiphase flow in porous media. Similarly, use of a repulsive force is found to be unnecessary in simulations involving the migration of two liquid phases. For multiphase flows incorporating a gas phase it is recommended to use a non-zero value for the repulslive force to avoid development of an ill-conditioned nodal distribution matrix. An equation to evaluate the repulsive force under these circumstances is suggested.     

Ensemble phase averaged equations for multiphase flows in porous media. Part 1: The bundle-of-tubes model

A bundle-of-tubes construct is used as a model system to study ensemble averaged equations for multiphase flow in a porous material. Momentum equations for the fluid phases obtained from the method are similar to Darcy’s law, but with additional terms. We study properties of the additional terms, and the conditions under which the averaged equations can be approximated by the diffusion model or the extended Darcy’s law as often used in models for multiphase flows in porous media. Although the bundle-of-tubes model is perhaps the simplest model for a porous material, the ensemble averaged equation technique developed in this paper assumes the very same form in more general treatments described in Part 2 of the present work (Zhang, D.Z., 2009. Ensemble Phase Averaged Equations for Multiphase Flows in Porous Media, Part 2: A General Theory. Int. J. Multiphase Flow 35, 640–649). Any model equation system intended for the more general cases must be understood and tested first using simple models. The concept of ensemble phase averaging is dissected here in physical terms, without involved mathematics through its application to the idealized bundle-of-tubes model for multiphase flow in porous media.     

Ensemble phase averaged equations for multiphase flows in porous media. Part 2: A general theory

Most models for multiphase flows in a porous medium are based on a straightforward extension of Darcy’s law, in which each fluid phase is driven by its own pressure gradient. The pressure difference between the phases is thought to be an effect of surface tension and is called capillary pressure. Independent of Darcy’s law, for liquid imbibition processes in a porous material, diffusion models are sometime used. In this paper, an ensemble phase averaging technique for continuous multiphase flows is applied to derive averaged equations and to examine the validity of the commonly used models. Closure for the averaged equations is quite complicated for general multiphase flows in a porous material. For flows with a small ratio of the characteristic length of the phase interfaces to the macroscopic length, the closure relations can be simplified significantly by an approximation with a second order error in this length ratio. This approximation reveals the information of the length scale separation obscured during an averaging process and leads to an equation system similar to Darcy’s law, but with additional terms. Based on interactions on phase interfaces, relations among closure quantities are studied.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

剧变截面圆管内渗流的数值计算方法

对于剧变截面圆管的渗流问题写出不可压缩渗流的基本方程组,对直接求解原始变量(速度和压力)的数值计算方法作出改进。先由非主流方向的运动方程计算压力,后由主流方向的运动方程计算主流方向的速度分量,再由连续性方程计算非主流方向的速度分量。这样可以避免在一般的求解原始变量方法中由连续性方程计算压力时出现的困难和麻烦。根据本方法和剧变截面圆管的特点,采用半交错不等距非正交贴体混合网格系。本文详细写出差分方程和迭代计算公式,对剧变截面圆管内的渗流算例进行数值计算。本方法的优点是简单和实用,在工程上具有较大的应用价值。     

Flux calculation at the interface between two rock types for two-phase flow in porous media

When simulating two-phase flow in porous media, one has to consider the case where there is a discontinuity in the medium. There relative permeabilities and capillary pressure functions may change and we address the problem of calculating the convective part of the numerical flux at the interface between the two rock types. Several solutions are compared.