Natural Gas Technology - Performance Predictions for Gas Reservoirs Considering Two-Dimensional Unsteady-State Flow

Methods are presented for calculating the performance of multiwell gas reservoirs. These methods account for two-dimensional, unsteady-state flow of a non-ideal gas through a heterogeneous formation. The performance can be obtained for gas reservoirs with known physical properties and configuration subjected to specified development and production schedules. These methods may prove helpful in optimizing the exploitation of gas reservoirs. The numerical solrrtion procedures given for this non-linear problem are extensions of known techniques for the linear diffusion equation. Two of these procedures are based on a new difference equation given by Saul'ev. Limitations and advantages of these procedures are discussed. An example solution is presented and is used to compare the calculation methods. INTRODUCTION Li planning the operation of a multiwell gas reservoir, it is desirable to predict its performance for each alternate development and producing schedule under consideration. This paper presents three numerical procedures that can be used to calculate such gas reservoir performance using a digital computer. These methods account for two-dimensional transient flow of a non-ideal gas within the reservoir. Extension of these methods to three-dimensional flow is possible. The paper is divided into three parts. The first describes the problem for which a solution is wanted and three numerical procedures that can be used to obtain solutions. The second discusses the information necessary to set up a problem. This part will be of greatest interest to the field engineer. The third gives data and results for an example problem. PROBLEM AND SOLUTION METHODS In the following development, a gas reservoir is regarded as a sealed, porous, permeable, heterogeneous body in which a single, non-ideal gas phase is flowing. Flow is permitted in the x- and y-horizontal directions, but is considered to be negligible in the vertical or z-direction. The reservoir can be tapped by one or more wells. Methods to be described provide a basis for calculating the complete performance of the reservoir, including rates and cumulative production from the wells, and pressure distribution and decline in the reservoir. If the isothermal flow of a non-ideal gas obeying Darcy's Law occurs in two dimensions in a region R, the governing equation in this region is: