Solution of Two-Phase Flow Problems Using Implicit Difference Equations

Abstract

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. Introduction Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of particle differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow pressure and saturation. Thus, solution of the flow problems by difference equations involves solving problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems or problems of gas percolation. Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the net, level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. SPEJ P. 417

Enabling interactive and collaborative oil reservoir simulations on the Grid

“…A fully implicit approach for the solution of the immiscible flow equations was first introduced by Blair and Weinaug [21]. Aziz and Settari [22] performed a stability analysis of the linearized semi-implicit and fully implicit methods and concluded that they all satisfy material balance and are unconditionally stable.…”

Simulation of air storage aquifer by finite element model

“…Subscripts a and w refer to air and water, respectively. The relationship between air pressure, pa, water pressure, pw and capillary pressure, p c is given by : (2) can be expressed in terms of the dependent variables if it is assumed that (a) hysteresis effect is negligible, and (b) isothermal conditions prevail. Furthermore, it is assumed that deformation of the porous medium is very small so that the porosity may be kept constant.…”

Acquisition and analysis of transient data through unsteady-state core flooding experiments

“…Relative permeability values are computed directly from dynamic displacement data. However, these methods have two primary limits as follows: (a) analysis are limited to Buckley-Leverett model extended by Welge (1952) and not appropriate for low-flow rate experiments with significant effect of capillary pressure; (b) these methods require numerical or graphical differentiation of experimental data (Blair and weinaug 1969). The other implicit method was developed to overstep the limitations that are associated with the explicit approach.…”