An Improved Method for Simulating Reservoir Pressures Through the Incorporation of Analytical Well Functions
Authors: Hales Hugh B. 
Published: December 2007

Abstract
The simulation of reservoir performance is particularly difficult, compared with the simulation of similar processes, because of the extreme behavior of reservoir pressures and saturations. In particular, the reservoir pressures are very non-linear and exhibit near-singularities at the wells as a result of injection and production. Finite difference methods, which are used almost universally for reservoir simulation, are troublesome with such highly non-linear solutions. Accuracy is reduced and the solutions become more time consuming.
This paper describes an attempt to eliminate these nonlinearities from the finite difference solutions by incorporating analytical solutions in the pressure solution. The reservoir pressures are the sum of analytical functions and traditional finite difference solutions. The analytical functions are based on well test equations which describe the pressures around the wells in an infinite, homogeneous system. Although heterogeneities and reservoir geometries make these analytical functions grossly inaccurate for predicting actual reservoir pressures, they do exhibit nonlinearities much like the actual solutions. As a result, the finite difference components become much more linear, improving both the accuracy and the speed of the solution. Several other advantages of this new technology are also demonstrated:
- Empirical well equations, which relate cell pressures to well pressures, are unnecessary.
- The technique is applicable to Cartesian grids as well as other grid types.
- Wells need not be centered in the cells for best results.
- Accurate simulations of early-time, pressure transients make well test simulations practical.

Introduction
Despite the similarity of reservoir simulation with other engineering problems such as laminar fluid flow, conductive heat transfer, and convective/diffusive mass transfer, reservoir simulation is uniquely difficult. These difficulties must be due, at least in part, to the large non-linearities inherent in the solutions. For example, the reservoir pressures develop near-singularities at the wells as illustrated in Figure 1. Unlike the other engineering problems, reservoir simulation results depend on flow through very small areas, the well bores. The wells act much as line sources resulting in the very sharp pressure gradients near them.
The accuracy of finite difference solutions deteriorates when the solutions are highly non-linear, as can be easily demonstrated through the Taylor's series derivation of the finite difference approximations to the partial differential terms. The solution speed is also affected. This work investigates the elimination of the near-singularities in the finite difference equations, through the combination of analytical functions. As illustrated in Figure 1, an analytical function which accurately represents the pressure gradients around the well bores is added to the finite difference solution. The analytical solution may not accurately represent the actual reservoir pressure, particularly at large distances from the wells, but it does reduce the non-linear behavior of the finite difference component of the solution. Aftab has tried a similar approach to eliminate both the pressure singularities and the saturation discontinuities.