Measurements of Physical Properties - Calculations of Unsteady-State Gas Flow Through Porous Media

The problem of unsteady-state gas flow through porous media leads to a second-order non-linear partial differential equation for which no analytical solution has been found. In this paper a stable numerical procedure is developed for solving the equation for production of gas at constant rate from linear and radial systems. An electronic digital computer is used to perform the numerical integration using an implicit form of an approximating difference equation. Solutions are presented in graphical form for various values of dimension-less parameters. The solutions are compared with the laboratory study of gas depletion in a linear system. INTRODUCTION Production of fluids from porous rock reservoirs is essentially a transient process. Transient gradients develop as soon as production begins, and further withdrawals continue to cause disturbances which propagate throughout the reservoir, each adding in some way to the prior ones. A correct mathematical analysis of this behavior is complicated by the fact that the transient or unsteady-state flow of compressible fluids must be described by difficult second-order partial differential equations. As a practical matter, three distinctly different cases arise: 1. Flow of single-phase liquid 2. Flow of gases 3. Multiphase flow The first of these has been found to give a linear second-order equation similar to the well-known heat flow equation. where 7 is fluid density, 8 is time, and a a constant for the system, provided d 7 = r cdp, where the compressibility, c, is constant over the range of pressure, p, considered. Solutions of Equation (1) for both linear and radial flow are available in several forms.2'3' On the other hand. the second case, which is the flow of gas, gives a non-linear second-order equation, the solution of which is not known. Although a number of approximate solutions have been proposed, each is limited in value by the associated simplifying assumptions. Inasmuch as the analysis of transient flow is limited to liquid systems, a solution of the second case is necessary if further progress is to be made in studying underground fluid movement. For this reason a solution of Equation (2) was undertaken by means of numerical integration of approximating difference equations. BASIC DIFFERENTIAL EQUATION Equation (2; is derived by combining the equation of continuity, the perfect gas law, and Darcy's law to give the basic differential equation Strictly speaking, the fluids present in a gas reservoir are not perfect gases. Furthermore, there are indications that over certain velocity ranges, Darcy's law is not applicable. It would not be practicable. however, to attempt to obtain numerical solutions which would take into account all the possible variations from these laws which could occur in actual systems. The