Use of Two-Dimensional Methods for Calculating Well Coning Behavior
Authors: Welge Henry J. , Weber A.G. 
Published: December 2007

Abstract
A published calculation method for predicting incompressible, multidimensional fluid displacement has been adapted to the problems of water and gas coning in oil wells. Since depth and radial distance from the wellbore are the two key dimensions affecting the shape of a gas or water cone, coning calculations are well suited to the use of two-dimensional methods. The alternating direction implicit procedure (ADIP) was used for relaxation calculations of two-phase potentials in a two-dimensional grid. From the potentials and the capillary pressure relationship, saturation and pressure distributions were calculated which trace cone growth with time. Predictions have been made for well producing histories both before and after core breakthrough. To check validity of the method, two-dimensional calculations have matched the coning behavior and produced water-oil ratio history of a laboratory sand-packed model. They have also matched the coning behavior of several producing wells, for which the calculations were compared with produced water or gas cuts and logs showing water or gas cone movements.
Introduction
The application of two-phase, two-dimensional calculations using ADIP to various reservoir flow problems has been described in the literature. When adapted for computer solution, this method has proved to be a powerful tool for simulating well and reservoir behavior. This paper discusses the method as applied to well coning calculations. Single-well and coning calculations comprise an especially difficult class of two-dimensional problems which require special techniques for computer calculation and determination of reservoir characteristics. Refs. 4 through 8 describe previous approaches to the coning problem. Several examples of water and gas coning calculations, including studies on both laboratory models and producing wells, are presented here. The two-dimensional method accounts realistically for the most critical parameters affecting coning behavior, including production rate, formation stratification, horizontal and vertical permeabilities, depth of well penetration, gravity and capillary forces. The method considers the different densities and viscosities of the two phases and the relative permeability and capillary pressure characteristics of the rock and fluids. In addition to tracing cone growth in the vicinity of the wellbore, the method calculates the overall movement of the fluid interface throughout the well's drainage volume. Incompressible fluid flow is assumed to occur between the producing interval and the well' s limit of drainage. Calculations can be made for the producing history both before and after cone breakthroughA typical two-dimensional grid or array of blocks used to solve a coning problem contains about 400 blocks. The well's cylindrical drainage volume can be represented by about 20 radial subdivisions and the formation thickness by 20 vertical subdivisions. The grid spacing is normally smaller near the withdrawal interval to define the cone shape accurately. For this work we used an IBM 7074 digital computer having a core memory of 10,000 ten-digit words. A typical study, covering 5 to 10 years of well producing history, required from three to six hours of computing time.
MATHEMATICAL SIMULATION OF CONING BEHAVIOR
BASIC METHOD
In coning calculations, the reservoir volume drained by the producing well is represented by a two-dimensional system of blocks as shown in Fig. 1 for water coning studies. The horizontal dimensions of the blocks increase with radial distance from the well axis in geometric progression, i.e., the block size is small near the wellbore and large near the well's drainage radius (re).
SPEJ
P. 345ˆ