Abstract A mathematical model is presented and solved for determination of the saturation distribution and pressure in a radial gas-storage reservoir. The model consists basically of two parts:the growing gas-bubble core, andthe surrounding aquifer. Since the total pressure at the injection well is a function of the two-phase flow in the gas bubble and the unsteady single-phase flow in the aquifer, the resistance to flow in both zones was taken into consideration. The assumptions involved for both the radial equivalent of the Buckley-Leverett two-phase-flow equation and the injection-pressure equation are as follows:the geometry is radial,the gas bubble is free to expand or contract,compression or expansion of the gas within the bubble may occur at the beginning of a time step,the fluids are immiscible,water is incompressible within the gas-storage region whereas it is compressible outside of this region,a stable gas-water interface exists andgas injection occurs at a constant rate or a series of constant rates. The mathematical model was solved numerically using an IBM 650 computer. A comparison is presented between the predicted results of the model, the results assuming steady-state flow and the actual initial injection-pressure history of an operating reservoir. Using the initial field pressure for a basis, the average deviation between the predicted pressure and the actual field pressure was less than 4.3 per cent. Introduction Underground storage of natural gas in abandoned oil fields, in abandoned coal mines, in caverns and in aquifers has had varying degrees of success. The lack of these first three facilities in the vicinity of most major gas marketing areas is leading more and more to storage in virgin aquifers. For the most part, this storage in aquifers has led to the study upon which this paper is based. Several excellent papers have been published discussing the problems involved in the underground storage of natural gas; however, to these writers' knowledge, none of them considered the two-phase flow of fluids in the gas bubble and the fact that the gas bubble will grow or shrink in size depending upon the injection and withdrawal history of the reservoir. Not only will both of these factors influence the required injection pressure, but also they will have a definite bearing upon the amount of water production occurring upon withdrawal from the storage area. To solve this problem on an IBM 650 computer, several simplifying assumptions were made. Only the drainage portion of the relative-permeability curve was considered, thereby neglecting any hysteresis effects introduced during an imbibition cycle. The quantity of gas withdrawn during the period of the field study was only a small percentage of the total gas in place; therefore, this assumption is justified. For the problem studied in this case, the combined effects of capillarity and gravity were assumed to be negligible. This assumption becomes less valid as the formation thickness or the average pore diameter increases markedly. The justification of these assumptions should be studied in each case. THEORY This mathematical model for the prediction of saturation distribution and pressure for a radial gas-storage aquifer is based upon the equations for radial two-phase fluid flow and for radial, unsteady-state, single-phase fluid flow. The two-phase flow is considered to take place in a "core" whose radius is equal to the maximum radius that the gas zone will attain, as indicated in Fig. 1. The injection-withdrawal history is approximated by a series of constant flow rates. Within the "core", gas is assumed to behave as a semi-compressible fluid; that is, the gas is assumed to have a constant density based on the average gas-zone pressure for the flow period. Between each constant-rate time increment, the gas density is allowed to change. This involves an iterative solution of the pressure equation. The liquid phase in this zone is considered to be incompressible, whereas it behaves as a compressible liquid outside of this region. JPT P. 1389^