Techniques for Fully Implicit Reservoir Simulation

Techniques For Fully Implicit Reservoir Simulation A.D.K. Au; A.D.K. Au Computer Modelling Group Search for other works by this author on: This Site Google Scholar G.A. Behie; G.A. Behie Computer Modelling Group Search for other works by this author on: This Site Google Scholar B. Rubin; B. Rubin Computer Modelling Group Search for other works by this author on: This Site Google Scholar P.K. Vinsome P.K. Vinsome Computer Modelling Group Search for other works by this author on: This Site Google Scholar Paper presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, September 1980. Paper Number: SPE-9302-MS https://doi.org/10.2118/9302-MS Published: September 21 1980 Cite View This Citation Add to Citation Manager Share Icon Share Twitter LinkedIn Get Permissions Search Site Citation Au, A.D.K., Behie, G.A., Rubin, B., and P.K. Vinsome. "Techniques For Fully Implicit Reservoir Simulation." Paper presented at the SPE Annual Technical Conference and Exhibition, Dallas, Texas, September 1980. doi: https://doi.org/10.2118/9302-MS Download citation file: Ris (Zotero) Reference Manager EasyBib Bookends Mendeley Papers EndNote RefWorks BibTex Search Dropdown Menu nav search search input Search input auto suggest search filter All ContentAll ProceedingsSociety of Petroleum Engineers (SPE)SPE Annual Technical Conference and Exhibition Search Advanced Search AbstractThis paper summarizes the methods that have been used to develop fully implicit and highly implicit reservoir simulators.The techniques developed emphasize simplicity and clarity. The governing equations, written in fully implicit form are treated as functions. These coupled nonlinear functions are solved by Newtonian iteration, the derivatives of the Jacobian matrix being evaluated numerically. A procedure is described for the efficient evaluation of the Jacobian.A technique is presented for avoiding variable substitution in variable bubble point problems. This is accomplished by means of a novel pseudo solution gas formulation.An efficient method for including two-point upstream into fully implicit simulators, without increasing the Jacobian bandwidth, is described. The scheme is based on a fully implicit upstream mobility plus an explicit correction term derived from the plus an explicit correction term derived from the two-point mobility formulation. In addition a new method called centralized upstream is given. The new method has less truncation error then two-point upstream, and is little more difficult to implement.A modified Crank-Nicholson method is described for reducing the time truncation error associated with large time steps in implicit models. The modified version does not exhibit the usual stability problems associated with the Crank-Nicholson method on nonlinear problems.A section on well models shows how a simple enhancement of conventional well models can be used to account for the position of a well within a grid It is also shown how to introduce a multiblock completion into an implicit simulator without troying the block-banded form of the Jacobian, or without increasing the bandwidth.An efficient iterative solution method is described for the sparse block-structured Jacobians. The method extends the SIP factorization to block form and uses minimizations and orthogonalizations to accelerate the convergence rate.Several applications and examples are provided. A water coning example is included to show the stability and convergence characteristics of the Newtonian iteration. A vertical cross-section oil resaturation problem with variable bubble point is used to illustrate the pseudo solution gas formulation. The modified Crank-Nicholson method is applied to a one-dimensional Buckley-Leverett type problem. problemIntroductionIt is well known that fully implicit methods are remarkably stable and can tolerate much larger time steps than those used in explicit formulations. However, fully implicit methods are not widely applied in large scale numerical problems such as oil reservoir simulation. The larger matrix bandwidth associated with fully implicit approach demands more computation effort per time step and computer core memory compared with explicit methods. Other considerations such as time truncation error associated with the large time step size and the difficulties in the implementation of higher order methods to reduce spatial truncation have also made the fully implicit method less attractive.A collection of techniques for resolving these difficulties in fully implicit simulation is contained in this paper.GENERAL CHARACTERISTICS OF NEWTONIAN ITERATIONThe Newton algorithm is a well known iterative process for the numerical solution of equations. process for the numerical solution of equations. Keywords: PVT measurement, bubble point, implicit simulator, spe 9302, reservoir simulation, Newtonian iteration, time step, iteration, Modeling & Simulation, two-point upstream Subjects: Fluid Characterization, Reservoir Simulation, Phase behavior and PVT measurements This content is only available via PDF. 1980. Society of Petroleum Engineers You can access this article if you purchase or spend a download.