Search Documents

By A. F. Bertuzzi, M. R. Tek, F. H. Poettmann

A method is presented for predicting pressure drop for two-phase fluid flow in horizontal pipes. A set of 267 experimental measurements randomly sampled from approximately 1,000 measurements from various literature sources was used. Pressure gradients calculated by the procedure developed and compared with the experimental values showed a bias of +0.82 per cent and a standard deviation of 20.8 per cent. The advantatges of this method over other available methods of predicting two-phase flow pressure drop are (a) its comparative simplicity of application, (b) its relative independence of flow patterns, (c) its accuracy, which on the basis of a statistical evaluation predicts pressure drop closer than other available methods, and, (d) its ability to satisfactorily correlate laboratory data from various sources while the correlations from these sources do not appear to agree with one another. The method has been reduced to a simplified graphical procedure, suitable for field use. A two-phase f factor is defined and correlated with parameters involving the flowing gas-liquid mass ratio, a Reynolds number for the gas phase, and a Reynolds number for the liquid phase. The choice of parameters allows the correlation to reduce to the usual f factor plot for the limiting conditions of (111 gas or all liquid. INTRODUCTION The mechanics and characteristics of two-phase flow systems have been of interest throughout the industry for some time. In numerous engineering installations such as pipe lines, chemical reactors, and heat exchangers, two-phase flow conditions are of every day occurrence. In oil production operations it has been desirable, in some cases, to consider transporting gas and oil together in a common pipe from oil field to process plant. The trend toward centrally located stock tank batteries in oil fields has resulted in longer gathering pipelines in which more than one fluid phase is flowing. The increase in the producing capacity of oil wells due to new production tech- niques has created the need for review and re-design of many surface gathering lines for properly handling the increased production. Optimum pipe size for the situations described above has become an important factor. The problem which is of interest here involves the ability to predict the relationship between pressure drop, fluid properties, fluid rates, pipe diameter, and pipe length. Although the literature1,5,10,13,18,19 contains various articles dealing with specific cases of the problem, Lock-hart and Martinelli17 proposed the most general solution. Their method has been later modified by Baker:1 Alves2 demonstrated that different flow patterns are possible for a flow mechanism defined by Martinelli such as gas turbulent, liquid turbulent. Baker explained deviations experienced with Martinelli's correlations on the basis that different flow patterns could exist for the same flow mechanism. Using Martinelli's method of correlatioil, Baker correlated data for each flow pattern for the turbulent-turbulent flow mechanism. Bergelin and Gazley6,7 presented a correlation for predicting gas phase

Jan 1, 1957

By R. D. Carter

Nurnerical .solutions are presented to some problems of unsteady-state radial flow of gas in which Darcy's law holds. These solutions are intended to aid in explaining the observed behavior of gas wells during drawdown tests. The variation of viscosity and Z-factor with pressure is accounted for in the solutions. Results are presented for three types of problems: (I) unfractured wells producing at a constant rate, (2) fractured wells producing at a constant rate and (3) unfractured wells producing through a critical flow prover. In the critical-pow-prover solutions, a correction of the calculated well-bore pressure to account for non-Darcy flow is made. A good correlation is shown between the solutions to the problem of Type 1 and the corresponding diffusivity-equation solution. Using this correlation, the diffusivity solution may be used to obtain gas-well solutions. The solutions to the problems of Type 2 illustrate the difference in drawdown-curve behavior to be expected between fractured and unfractured gas wells. Methods for estimatting fracture radius and flow capacity are presented and verified from solutions of Type 2. The solutions to the problems of Type 3 illustrate the differences to be expected between critical-flow-prover drawdown tests and constant-rate tests. An equation is presented for correcting critical-flow-prover data to constant-rate conditions. A description of the methods used in obtaining the solulions appears in the Appendix. INTRODUCTION It is necessary for many engineering purposes, such as predicting future deliverability performance, to be able to relate bottom-hole flowing pressure, formation pressure and flow rate for conditions of stabilized flow in gas wells. This information is customarily given as a plot of on logarithmic co-ordinates. It is desirable to be able to predict stabilized performance curves for gas wells based on data obtained from drawdown tests of short duration. Tests of long duration may be wasteful of gas, and low-permeability wells do not stabilize in a short-duration test. Although it is now commonly recognized that non-Darcy flow plays a significant part in gas-well behavior, it is believed that an understanding of transient behavior of systems obeying Darcy's law is still important in the interpretation of short-term test data, and the prediction of stabilized performance curves therefrom. In this connection, an equation for correcting Darcy flow It seems likely that solutions of the type presented herein, combined with the non-Darcy flow correction of the type of Eq. 1 and the usual laminar-flow skin factor, could provide a model which would adequately explain transient gas-well behavior for times after the non-Darcy flow region (which should be confined to a small area surrounding the well) has been established. The solutions presented in this paper are for three types of problems: (1) unfractured wells producing at a constant rate, (2) fractured wells producing at a constant rate and (3) unfractured wells producing through a critical flow prover. In the problems of Types I and 2, no corrections for non-Darcy flow were made in the results presented, since it was reasoned that the constant-rate boundary condition would facilitate such a correction if it was desired later. For the problems of Type 3, however, a condition of varying flow rate required that a non-Darcy flow correction of the type of Eq. 1 be included. This was done. Based on the numerical results presented, a number of approximations and a correlation are presented which may be useful in analyzing gas-well test data. It is assumed that the system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial.

By R. C. Craze

This paper discusses the characteristics of the velocity logs now available to the petroleum industry, and some of their advantages and limitations. The velocity log was designed as an aid to the geophysicist to determine velocity layering in the earth and to provide fast, accurate and economical vertical travel time data. After a number of logs had been made, studies were carried out on their use ia correlating strata in adjacent wells. It was found that adjacent strata usually show sufficient velocity contrast to make it possible to correlate these from well to well, even in areas where electric or radioactive logs fail.' In addition, velocity logs have been useful in the determination of porosity and thickness of reservoir rock.' These uses for velocity logs require an understanding of the characteristics of the log in response to the changes in velocity and in borehole conditions encountered normally in logging. GENERAL DISCUSSION OF INSTRUMENTATION The measurement made in velocity logging is more accurately described as a measurement of time rather than velocity. Ideally the log is obtained by measuring the length of time required for compressional sound waves to travel through a fixed distance in the formation. In practice, the time measurement may be made accurately but the distance through which the waves have traveled may vary slightly. Changes in borehole diameter are responsible largely for the difficulty in maintaining a fixed distance, especially if the diameter variations are so abrupt and frequent that the borehole instrument cannot follow contours of the hole. Several service companies have velocity logging equipment available or under development which have different transmitters, receivers, and measuring systems. While these differences are important, there are actually only two systems in use at the present time (Fig. I)."' The first system initiates a pulse at a trans-mitter and detects it at a single receiver. The transmitter and receiver are separated by a fixed distance, called the spacing. The time measured is that between the initiation of the pulse and the first arrival of acoustical energy at the receiver. The second system uses two receivers and a transmitter. The receivers are separated by a fixed distance and the time between arrivals of energy at the two receivers is measured. Here, the spacing is the distance between receivers. The two systems have much in common for it is possible to record both logs simultaneously. Usually the transmitter is an electro-mechanical transducer which produces discrete

By M. P. Tixier, K. S. Kunz

A method for the interpretation of temperature curves recorded in gas producing wells is described. One essen-tial feature of the method is a simple graphical construction which conveniently reflects the amount of gas flowing from a given producing formation. This method, applied in conjunction with the Induction Log and radioactivity logs, is valuable to determine the points of entry of gas and to estimate the over-all thickness of the producing zones and their respective contributions to the total production. The method in particular is helpful in evaluating the effect of fracturing operations. An outstanding result of the application of the method is to bring to the light the relative importance of sand-shale interfaces in the production of gas. The paper is illustrated with field examp1es. INTRODUCTION The electrical logging tools involving conventional systems are not applied in wells which do not contain water base mud (empty holes, or wells filled with oil base mud), because of the absence of electrical connection between the electrodes and the surrounding formations. It is, therefore, common practice in this case to use other methods which do not require the presence of a conductive fluid in the borehole: resistivity log with scratcher electrodes, Induction Log, radioactivity logs, temperature log. . . . Such conditions are encountered in some regions, as the San Juan Basin, N. Mex., where a technique of drilling gas wells, using gas as a drilling fluid, has been introduced in the last two years. By the middle of 1953, an appropriate method for the logging of these wells was devised on the basis of field tests and of theoretical studies. It is the purpose of the present paper to give a description of this method. In these wells the logging operations are performed during the process of gas production, so that the temperatures in the boreholes are strongly affected by the cooling due to the expansion of the gas. On the other hand, a special thermometer, with a high sensivity and a short time constant, has been recently made available, which has enabled much more detailed and accurate recordings of the borehole temperatures than with the conventional instruments. As a result of these two circumstances, the temperature log constitutes the essential element of the method, in making possible the detection of the gas producing zones, the definition of their boundaries, and the estimation of the respective contributions of each separate zone to the total production. The method also involves the radioactivity and Induction Logs, as auxiliary documents. These logs in the present instance are helpful in the determination of the reservoirs. However, the information they give in this respect is limited. This comes from the fact that, due to the absence of the MicroLog and/or MicroLaterolog and of the SP log, the estimation of formation factors and connate water salinities is somewhat conjectural. Furthermore, little help is obtained from the Neutron log for the evaluation of porosity and of formation factor because the variations of gas saturation affect the Neutron log in a way similar to the variations of porosity. It should be said also that even an exact knowledge of the gas saturation of the reservoir beds would not

Jan 1, 1956

By L. Cichowicz, K. K. Millheim

The constant-rate drawdown test performance for a low-permeability, verticany fractured gas well was investigated. A series of gar wells were tested by flowing each well at constant rate until the data could be analyzed using con-ventional radial flow theory. Each well was then shut in to build up. After a sufficient buildup was obtained, another flow Iest commenced but at a higher flow rate than Ihe first test. Again, the well was shut in when radial fiow was obtained. his procedure was repeated for three to four different flow rates. Two wells in the San Juan basin were tested using this procedure. Both wellere fractured after completion, cleaned up and then shut in until flow testing commenced. Test designs of both wells permitted investigation of the most realistic values of egective permeability, wellbore radius and turbulence factor. Also, being able to determine the eflective fracture flow area and vertical fraclure efficiency was inherent with this testing approach. It was observed that fractures in both wells influenced the Pressure behavior for approximately Is to 40 hours (depending on the flow rate) before radial flow was evident. After this time, drawdown data were analyzed using radial pow theory. When a low-permeability gas well has vertitally oriented, induced fractures, the early flow geometry ic. essentially linear. It will be shown how to determine when a flow test has been conducted long enough so lha' the most representative values of effective permeahiiity. wellbore radius and turbulence factw can be calculated. From the linear pressure data, valuable information about The fracture treatment, such as the effective flow area and vertical fracture efficiency, can be determined for vertically froctured wells. Introduction During tests on gas wells in the Sari Juan basin7 initial transient behavior lasted for many days because of the low permeability of some porous media. As a result, stabilized flow performance could not be obtained. If these wells received some type of stimulation treatment, early pressure behavior deviated from conventional theoretical radial Row. When conventional radial flow theory was used to analyze these low-permeability fractured gas wells, larger values of flow capacity and absolute open flow potentials (AOF) sometimes resulted. Wells were assigned open flow poten- tials that proved to be 3 to 10 times higher than the well would sustain over a longer period of production. In some cases where the wells had flowed for longer periods of time during a constant-rate drawdown test, it was noticed that the effective flow capacity appeared to be decreasing with time until a certain value was reached. The early nonradial pressure behavior can be explained If linear flow is assumed. Rusell and Truitt mathematic. ally investigated the vertically fractured well in a bounded area. They showed that early flow behavior was essentially linear and, for x,/x. approximately less than 0.10 radial flow was obtained after short periods of time. Then realistic values of effective permeability and skin could be determined. Scotta experimentaliy studied the vertically fractured well with a heat flow analog. He showed that earb flow was linear, Both studies indicate that, for small values of x,/x,, linear flow approaches radial flow if the well is tested long enough, To help prove this concept of early linear flow caused by induced vertical fractures, two low-permeability gas wells were tested. Both wells received large fracture treatments prior to testing. A vertical fracture was indicated from the analysis of fracture treatments, As anticipated, tests of both wells indicated early linear flow that was later followed by a period of radial flow, Data collected from each well were analyzed. From the well tests, plus other information on each well, the effective permeability, wellbore radius and turbulence factor were calculated. Effective fracture flow areas calculated from test analyses for each well proved to be approximately one-fourth the created area calculated from classic hydraulic fracturing theory: other fractured wells that were tested but not presented in this paper also indicated that the effective fracture flow area was one-fourth to one-third the created area predicted from hydraulic fracturing theory. The vertical fracturing efficiency was estimated from the calculated values of effective wellbore radius and fracture flow area. For the two wells tested, calculated fracture lengths x, were 112 and 105 ft, and the vertical fracturing efficiencies E, were 122 and 183 percent. Development of Flow Model Agnew' showed that most induced fractures below 1,500 ft are vertical. Anderson and Stah15 indicated that most of the fractures they studied were vertical. Thc model proposed for early flow in most vertically fractured gas wells is shown by Fig. 1. This model should approximate early flow behavior until radial flow is reached, at which time a radial model with an effective wellbore radius of 0.5 they will

Jan 1, 1969

By G. Thodos, R. B. Grieves

A method has been developed for the accurate calculation of the cricondentberm and cricondenbar temperatures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain any number of components including paraffins and isoparaffins, olefins, acetylenes, naphthenes and aromatics. This approach is based upon the mole fraction of the low-boiling component in the mixture and graphically presents the ratios of the cricondentherm and cricondenbar temperatures to the pseudocritical temperature as {unctions of a boiling-point parameter. The only requirements are the critical temperature, normal boiling point and the approximate vapor pressure behavior of each component. For mixtures of more than two constituents a stepwise calculation procedure is necessary where the pseudocritical temperature is based upon the critical temperature of the pure low-boiling component and upon the actual cricondentherm or cricondenbar temperature of the mixture of all the remaining higher-boiling components. From an analysis of 22 binary systems (123 compositionsl, the average deviation of calculated cricondentherm temperatures from reported values is 0.87 per cent, based on degrees Rankine; and for 10 multicomponent mixtures containing from three to six components, the average deviation is 0.98 per cent. From an analysis of 18 binary systems (104 compositions) the average deviation of calculated cricondenbar temperatures from reported values is 0. 79 per cent; and for 15 multicomponent mixtures, 1.47 per cent. Equations, derived from the graphical relationships, are presented which enable rapid calculations for both binary and multicomponent systems. INTRODUCTION With the increasing tendency towards the use of high pressure processes, the need for an accurate method of calculation of the properties of a multicomponent hydrocarbon mixture in the critical region is becoming more essential. Critical region phase relations are also significant in gas condensate reservoirs, for which a knowledge of the fluid phase boundaries and regions of retrograde condensation make it possible to evaluate optimum operating conditions. The possibility of liquefaction in the underground structure may be established, and information on the feasibility of the use of repressurizing may be determined. A typical phase diagram for a multicomponent hydrocarbon mixture is presented as Fig. 1. This figure also shows the regions of retrograde behavior encountered by following lines AB and AC. If accurate methods of calculating the temperatures and pressures at the critical, cricondentherm (maximum temperature) and cricondenbar (maximum pressure) points of multicomponent hydrocarbon mixtures of known composition were available, it would be possible to estimate closely the entire phase diagram of any mixture. Several methods have been presented in the literature for the estimation of the critical temperatures and pressures of multicomponent mixtures containing all types of hydrocarbons and including the fixed gases.4-8,21 Considering the cricondentherm and cricondenbar points, the work is more limited. Etter and Kay5 have developed rather complex equations for

By A. M. Rowe

A considerable quantity of experimental hydrocarbon K-factor data has been correlated as a function of component identity, temperature, pressure and convergence pressure. To utilize these correlations effectively, convergence pressure must be determined accurately, particularly for volatile mixtures near their critical states. This paper presents phase diagrams that illustrate physically the meaning of convergence pressure. A new method, referred to hereafter as the "critical composition method", will be outlined for calculating convergence pressure. An example calculation has been included to illustrate how to use this new technique. INTRODUCTION The principle of hydrocarbon phase composition calculations as applied to such diverse problems as optimizing separator performance or predicting fluid compositions at various stages of reservoir depletion is the same. Usually, temperature, pressure and the numbers of moles of the various components of the fluid contained in a given volume are known. Questions answered by the phase calculation are what fraction of the total mass of fluid exists in each of the equilibrium phases, and what are the mole fractions of the various components in the two phases? To answer these questions the Natural Gasoline Supplymen's Assoc. (NGSMA) has correlated a considerable quantity of K-factor data as a function of temperature, pressure, component identity and convergence pressure.l To use these correlations to obtain the best answers possible, one must be able to calculate the convergence pressure. This is particularly true for over-all fluid compositions in the neighborhood of the critical state. STATEMENT OF THEORY Use of convergence pressure as a correlating parameter is based on a postulate similar to the law of corresponding states used in correlating PVT data of hydrocarbon gases. This postulate, which proposes convergence pressure as a correlating parameter, has been stated as follows: "The equilibrium vaporization constant for one component in a complex system is the same as the equilibrium constant at the same temperature and pressure for the same number or kind of components, providing only that the convergence pressures of the two systems are exactly the same at the same temperature and that the components are of the same homologous series".2 This law, as with all laws of physics, cannot be proven theoretically. It can only be justified by experimental data supporting its premises. Arguments have been made that this law violates Gibb's phase rule. For example, consider a four-component system. By Gibb's phase rule, which is thermodynamically rigorous, f = c — P + 2 =4, for a four-component, two-phase system. Thus, four independent intensive variables must be specified to establish completely all the intensive variables of the equilibrium phases. On the other hand, according to convergence pressure theory, only three variables need be specified for any mixture containing four different components. These variables are temperature, pressure and convergence pressure. Thus, the two laws appear to be in conflict. However, the convergence pressure postulate is more restrictive than Gibb's phase rule. It applies only to mixtures of the same homologous series. Hence, these two concepts are not in disagreement. PHASE DIAGRAMS DESCRIBING CONVERGENCE PRESSURE This paper presents phase diagrams grading from the simple two-component system to the more complex four-component system to illustrate convergence pressure. A less abstract definition of convergence pressure will be stated that is

By R. E. Holmes, T. W. Leland

The solution to the unsteady-state mass- and heat-transfer equations describing the adsorption of a dilute component from a gas stream flowing through a packed bed is readily applicable to the design of hydrocarbon recovery units. The solution to these equations depends on the assumption of a linear equilibrium equation to relate concentrations on the surface to concentrations in the gas. This assumption is shown to be justifiable for components which are very dilute in the gas stream. This approach permits the prediction of the effects of bed shape and size, cycle time, and gas mass velocity on the recovery of each component. It also allows the time and gas flow rate needed for regeneration to be determined. The equations involve a constant for each component and another constant characterized by the type of adsorbent. These constants must be evaluated from operating data. Values of these constants determined for silica gel are presented. This paper points out the type of experitrrental data needed for a thorough design of hydrocarbon adsorption units. INTRODUCTION The growing importance of adsorption processes to recover hydrocarbons from natural-gas streams has created an urgent need for better methods to design and predict the performance of these processes. Since the original use of solid adsorbents in natural-gas processing was for dehydration, most of the current design methods for hydrocarbon recovery re extensions of those developed earlier for dehydration. These methods are based on an adsorptivity or bed loading factor defined as the amount of a particular component adsorbed per unit weight or per unit volume of the packed bed. The recovery of a given component is then expressed as an empirical function of its adsorptivity on a given adsorbent. In spite of the simplicity and initial utility of these methods, they are not completely satisfactory in that they do not predict the effect of the dimensions of the bed in relation to flow rate and cycle time, and do not predict the distribution of various components along the length of the bed. It is possible, however, to make a considerable improvement in the design of these processes by applying the basic equations of mass and heat transfer applicable to gases containing small concentrations of adsorbing components flowing through granular beds. These equations have been developed and solved by Anzelius,' Furnas,% nd Hougen and Marshall." where G = mass velocity of the inert constituents (assumed to be methane), lb/hr-sq ft, y = weight of an adsorbing component in the gas per unit weight of inert, lb of a given constituent per lb of methane, x = bed length, ft, p, = packed bed density, lb/cu ft, P, = gas density, Ib/cu ft, f = void fraction in the bed, t = time, hours, and w = weight of a given component adsorbed per unit weight of adsorbate free solid, Ib/lb pure adsorbent. The term on the left side of Eq. 1 represents the flow rate per unit bed volume of a component entering less the flow rate per unit bed volume of that component leaving the differential bed length. The terms on the right side represent, respectively, the accumulation rates per unit bed volume on the solid and in the gas phase within the differential bed section. The last term is negligible for this application since pgfv is quite small in comparison with p, and G. The partial differential is determined by the mass-transfer rate between gas and solid in the differential section.

By R. Al-Hussainy, P. B. Crawford, H. J. Ramey

The effect of variations of pressure-dependent viscosity and gas law deviation factor on the flow of real gases through porous media has been considered. A rigorous gas flow equation was developed which is a second order, non-linear partial differential equation with variable coeficients. This equation was reduced by a change of variable to a form similar to the diffusivity equation, but with potential-dependent diffusivity. The change of variable can be used as a new pseudo-pressure for gas flow which replaces pressure or pressure-squared as currently applied to gas flow. Substitution of the real gas pseudo-pressure has a nrtmber of important consequences. First, second degree pressure gradient terms which have commonly been neglected under the assumption that the pressure gradient is small everywhere in the flow system, are rigorously handled. Omission of second degree terms leads to verious errors in estimated pressure distributions for tight formations. Second, flow equations in terms of the real gas pseudo-pressure do not contain viscosity or gas law deviation factors, and thus avoid the need for selection of an average pressure to evaluate physical properties. Third, the real gas pseudo-pressure can be determined numerically in term of pseudo-reduced pressures and temperatures from existing physical property correlations to provide generally useful information. The real gas pseudo-pressure was determined by numerical integration and is presented in both tabular and graphical form in this paper. Finally, production of real gas can be correlated in terms of the real gas pseudo-pressure and shown to be similar to liquid flow as described by diffusivity equation solutions. Applications of the real gas pseudo-pressnre to radial flow systems under transient, steady-state or approximate pseudo-steady-state injection or production have been considered. Superposition of the linearized real gas flow solutions to generate variable rate performance was investigated and found satisfactory. This provides justification for pressure build-up testing. It is believed that the concept of the real gas pseudo-pressure will lead to improved interpretation of results of current gas well testing procedures, both steady and unsteady-state in nature, and improved forecasting of gas production. INTRODUCTION In recent years a considerable effort has been directed to the theory of isothermal flow of gases through porous media. The present state of knowledge is far from being fully developed. The difficulty lies in the non-linearity of partial differential equations which describe both real and ideal gas flow. Solutions which are available are approximate analytical solutions, graphical solutions, analogue solutions and numerical solutions. The earliest attempt to solve this problem involved the method of successions of steady states proposed by Muskat.' Approximate analytical solutions' were obtained by linearizing the flow equation for ideal gas to yield a diffusivity-type equation. Such solutions, though widely used and easy to apply to engineering problems, are of limited value bemuse of idealized assumptions and restrictions imposed upon the flow equation. The validity of linearized equations and the conditions under which their solutions apply have not been fully discussed in the literature. Approximate solutions are those of Heatherington et al.. MacRobertsl and Janicek and Katz.' A graphical solution of the linearized equation was given by Cornell and Katz. Also, by using the mean value of the time derivative in the flow equation, Rowan and Clegg' gave several simple approximate solutions. All the solutions were obtained assuming small pressure gradients and constant gas properties. Variation of gas properties with pressure has been neglected because of analytic difficulties. even in approximate analytic solutions. Green and Wilts8 used an electrical network for simulating one-dimensional flow of an ideal gas. Numerical methods using finite difference equations and digital computing techniques have been used extensively for solving both ideal and real gas equations. Aronofsky and Jenkins"I " and Bruce et al.11 gave numerical solutions for linear and radial gas flow. Douglas et al." gave a solution for a square drainage area. Aronofsky 13 included the effect of slippage on ideal gas flow. The most important contribution to the theory of flow of ideal gases through porous media was the conclusion reached by Aronofsky and Jenkins" that solutions for the liquid flow case'" could be used to generate approximate solutions for constant rate production of ideal gases. An equation describing the flow of real gases has been solved for special cases by a number of investigators using numerical methods. Aronofsky and Ferris 10 onsidered linear flow, while Aronofsky and Porter 17 considered radial gas flow. Gas properties were permitted to vary as linear functions of pressure. Recently, CarteP 18 proposed an empirical correlation by which gas well behavior can be estimated from solutions of the diffusivity equation using instantaneour values of pressure-dependent gas

Jan 1, 1967

By R. F. Hinds

A pressurizing system was designed and built to apply a radial pressure of 5.000 psi to rock samples. Samples of the Bradford, Weir and Kirkwood sandstones were subjected to radial pressures parallel and perpendicular to the bedding in the system and the changes in conductivity, porosity and permeability were deter-mined. Asymptotic decreases in conductivity, porosity and Klinkenberg permeability were observed over the 0- to 3,500-psi range. From these data, the increase of formation factor over the same range was calculated. The cementation exponent expressed in the Archie equation was found to increase with pressu,re. The rate of change varied from formation to formation. Only small differences were observed for samples taken parallel to the bedding as comnared to those taken perpendicular to it except for the Bradford sand. Above 3,500 psi the data are not good enough to tell if the properties continue to fall-off asymptotically or in a different manner. A possible explanation for the probable litnit of the asymptotic rate around 3,500 to 4,500 psi is that this represents the maximum pressure which the formations have unrdergone during their geological life as a result of burial. INTRODUCTION For many years measurements have been made in the laboratory of the properties of oil-bearing sandstones with little consideration being given to whether these properties had been changed when the samples were moved from their original position in the ground to the laboratory. Relatively little has been reported in the literature on the effect of pressure on the properties of sandstones. The changes in permeability for eight typical consolidated oil-bearing sandstone samples with overburden pressure (0- to 15,000-psi range) was presented by Fatt and Davis' in graphical form. In 1953 Fatt' also reported the relationship between porosity and overburden pressure for four typical consolidated oil-bearing sandstone samples (0- to 5,000-psi range). Fatt later described the changes of formation factor and porosity of 20 typical oil-bearing sandstone samples with overburden pressure changes (0- to 5,000-psi range) giving the values of the coefficients and exponents in the Archie' equation, and the more general Winsauer equation', where C, k and m are empirical constants, F is the formation factor, and + is the porosity. During the experiments of this investigation the core was contained in a neoprene rubber sleeve as shown in Fig. 1. It was held between two steel end members which were mounted in a steel cylinder. The unassembled pressure cylinder with its component parts is shown in Fig. 2. The core was insulated from the two steel end members by a Teflon sleeve with a small hole in the

By R. G. Agarwal, Ramey Jr. H. J., Al-Hussainy R.

Although it has long been realized that gas recovery from a water-drive gas reservoir may be poor because of high residual saturations under water drive, it appears that only limited infomlation on the subject has been available until recently. This study was performed to show the qiiantitative potential importance of water influx. Results indicate that gas recovery may be very low in some cases: perhaps as low as 45 per cent of the initial gas in place. Gas recovery under water drive appear to depend in an important was on: (I) the prodirction rate and manner of production; (2) the residual gas saturation; (3) aquifer propertie.); and (4) the volumetric displacement effciency of water invading the gas reservoir. The manner of estimating water-drive gas reservoir recovery can vary considerably. Examples are: the steady-state tnethorl. the Hurst modified steady-state method, and various unsteady-state methods such ac. those of van Ever-dingen-Hurst, Hurst, and Carter-Tracy. The Carter-Tracy water influx expression was used in this study. In certain cases, it appears that gas recovery can be increased significantly by controlling the production rate and manner of production. For this reason, the potential importance of water influx in particular gas reservoirs should he investigated early to permit adequtrtr planning lo optirtize the pay reserves. INTRODUCTION In recent years, the economic importance of natural gas production has become increasingly apparent. This has been evidenced by more intensive exploration efforts aimed at gas production, and exploitation of both deep, as well as low-permeability gas reservoirs. Technical developments such as deep-penetration fracturing have made development of such formations economically feasible. Unfortunately, water influx has forced abandonment of a number of gas reservoirs at extraordinarily high pressures. Although reservoir engineering methods for estimating water influx have long been available, it appears that application of these methods to the water-drive gas reservoir has been sporadic.'a Available methods for estimating water influx which can be applied to the water-drive gas reservoir problem include the steady-state method,1 the Hurst modified steady-state method and various unsteady-state methods such as those of van Everdingen-Hurst. Hurst, and Carter-Tracy. Interesting applications of these solu- tions to gas reservoir and the aquifer gas-storage problems have appeared recently.3,12,14 The experimental study of residual gas saturations under water drive by Geffen et al. in 1952 indicated that residual gas saturations could be extremely high. A value of 35 per cent of pore volume is often used in field practice when specific information is not available. The study of Geffen et al. showed that residual gas saturation might be much higher in some cases. Naar and Henderson concluded that the residual non-wetting phase saturation under imbibition should be about half of the initial non-wetting phase saturation. The Naar and Henderson result that residual gas saturation under water influx should be about half the original gas saturation is recommended as an estimate if laboratory measurements are not available. Thus, it is clear that a considerable portion of the initial gas in place might be trapped in a water-drive gas reservoir as residual gas at high pressure. A full water-drive would result in loss of residual gas trapped at initial reser.voir pressure. Consideration of transient aquifer behavior leads to the conclusion that high-rate production of water-drive gas reservoirs could result in improved gas recovery by reduction of the abandonment pressure. However, there appears to be little quantitative information on this possibility. One of the few advantages of water-drive gas production appears to be improved deliverability through water-drive support of the reservoir pressure. There may also be an advantage in higher condensate recovery caused by pressure maintenance for gas-condensate water-drive reservoirs. In view of the preceding, this study was made to assess the potential importance of water-drive in gas reservoir engineering. The Carter-Tracy approximate water-influx expression was used because this equation offers some advantages in hand-calculation which do not appear to have been generally recognized.' However, calculations were performed in the main with a high-speed digital computer to permit evaluation of the effect of water-drive under a large variety of conditions. CALCULATION METHOD Water-drive gas reservoir performance can be estimated in a manner completely analogous to oil reservoir calculations: a materials balance is written for the reservoir, and a water influx equation is written for the aquifer. Siniltaneous solution provides the cun~ulative water influx and the reservoir pressure. When reservoir performance data (gas produced and reservoir pressures) are available, it is usually possible to match performance data to determine the initial gas in place and the water influx parameters —

Jan 1, 1966

By M. H. Cullender

The performance characteristics 01 gas wells producing from formations which fail to stabilize within a relatively short period of time are obscured by the interrelation of the coefficient (C) and the slope "n" of the conventional back-pressure curve q = C(P2 f- p2s)n An empirical method is presented whereby the characteristic or true slope (11) of the back-pressure curve may be determined for a particular gas well. For those gas wells which do not stabilize within a relatively short period of time, the coefficient (C) is considered as a variable with respect to time and as a constant only with respect to a specific time. Thus, the back-presrure performance of such a well is presented as a series of parallel curve5, each curve representing the performance of the well at the end of a given time interval. Through use of the method presented, a simple presture gradient is developed and maintained within the drainage area around a producing gas well during the test period in order that the variation of the coefficient (C) with respect to time does not obscure the true value of the lope. INTRODUCTION The back-pressure method of testing gas wells as set forth by Schellhardt and Rawlins1 is dependent upon the requirement that a series of flow rates and corresponding pressure data be obtained under stabilized conditions. These data arc then plotted on logarithmic coordinates of flow rate (Q) vs difference in squared pressures (P2f — P2s) in order to determine the constants (C) and (n) for Equation 1. With respect to this method of presenting the performance characteristics of a gas well, it should be pointed out that although Equation 1 is an empirical relationship, the form of the equation was justified to a considerable extent by an independent development by Muskat and Botset.2 This development was theoretical in nature; however, certain assumptions with respect to the type of flow existing (viscous vs non-viscous) were madc which may or may not be correct. As the use of the method presented by Schellhardt and Rawlins spread through the industry, it became evident that the method of testing was applicable for those wells which approached stabilized prbducing conditions within a relatively short period of time. Performance characteristics could not be determined by this method, however, for wells which approached stabilized protlucing conditions slowly over a considerable period of time. This characteristic of "slow stabilization" has generally been associated with gas wells producing frorn reservoirs of low permeability. As the demand for natural gas continued to grow through the years, a great number of wells were completed in reservoirs of low permeability which exhibited "slow stabiliz:~tion" characteristics. Consequently, the problem of determining the true performancc characteristics of these wells became more apparent. It is with this problem in mind that the isochronal performance method of determining the flow characteristics of gas wells is presented.

Jan 1, 1956

By F. Kurata, A. F. Bertuzzi, P. C. Davis, T. L. Gore

The phase and volumetric properties of 10 very volatile mixtures are presented for temperatures from — 200°F to above the critical points. These mixtures consisted of natural gases and of mixtures of natural gases with methane and nitrogen. Existing correlations of compressibility factors have been extended to apply to these mixtures. Existing correlations of critical pressures and temperatures do not apply to these mixtures. Methods have been developed for estimating critical temperatures of these mixtures to ± 4°F and critical pressures to 2 5 per cent. These methods apply to other volatile and very volatile mixtures. The results of this work are particularly applicable to the separation of nitrogen from hydrocarbons by low-temperature distillation. INTRODUCTION Knowledge of the low-temperature phase and volumetric properties of natural gases and of other mixtures of the light hydrocarbons is becoming increasingly important as interest in low-temperature separation processes increases. Low-temperature distillation has been used to separate ethane and ethylene from natural gas and from refinery gases for petrochemical uses. Removal of nitrogen from natural gas by low-temperature distillation has been proposed 7,20 as a means of increasing the capacity of existing gas transmission lines and of improving the marketability of the gas. Further application of low-temperature separation processes has been retarded by lack of data on the phase and volumetric behavior of mixtures at low-temperatures and at high pressures. This paper presents the results of determination of the phase and volumetric properties of 10 very volatile mixtures at temperatures from - 200 °F to above the critical point. These mixtures consisted of natural gases with methane and nitrogen. Generalized methods are presented for estimating the critical temperatures, critical pressures, and the densities of such mixtures. EXPERIMENTAL APPARATUS AND PROCEDURES The apparatus and the experimental procedures have been described in an earlier paper13. The method consists of displacing the mixture from a reservoir which is maintained at constant temperature and pressure into an agitated glass equilibrium cell which is maintained at constant temperature. The liquid level within the glass equilibrium cell was visually observed, and the pressure was measured to within & 3 psi with a calibrated Bourdon gage. COMPOSITIONS OF MIXTURES STUDIED The compositions of the mixtures studied are summarized in Table 1. The measured properties included in Table 1 will be discussed later in this paper. The sources of these gases are described below: Gas "A-1'' was obtained from the Phillips Petroleum Co. Its composition was reported in connection with earlier work." The properties of this gas were reported in an earlier paper.13 Gases "A-2, A-3, A-4" were prepared by mixing nitrogen with gas "A-1." Gas "AB-1" was prepared by mixing methane with gas "A-1." The properties of this gas were reported in an earlier paper.'" Gases "AB-2" and "AB-3" were prepared by mixing nitrogen with gas "AB-1." Gas "B" was obtained by W. W. Bodle of J, F. Pritchard & Co. from the primary separator of a well in western Kansas. Gas "C" was obtained from a Cities Service Co. pipeline at Lawrence, Kans. Gas "D" was prepared by mixing gas "AB" with gas "C."

Jan 1, 1955

By G. Thodos, D. E. Matschke

Cansiderable time and effort frequently are expended to establish, with a degree of confidence, the PVT behavior of pure substances. In particular, a great deal of experimental information contributed by a number of investigators is available in the literature. The data of each investigator, although significant in themselves, cannot be used with a high degree of confidence until they are compared with the continuum of information presented by others. Current interest in thermodynamic and transport properties, for not only the gaseous state but the con-densed state as well, requires that PVT infarmation be presented for both regions. The density dependence on temperature and pressure, when presented as a continuum, is of considerable value in itself, especially since the available equations of state fail to properly account for the PVT behavior in the condensed state. For this reason and to incorporate PVT data recently presented in the literature, this investigation was carried out. In their studies concerned with the thermodynamic properties of methane, Matthews and Hurd in 1946 utilized the data of Kvalnes and Gaddy and Olds, Reamer, Sage and Lacey to establish the PVT behavim of methane in the superheated region and the saturated liquid densities of Keyes, Taylor and Smithl in conjunction with the Clapeyron equation to develop the densities for the saturated envelope. Current activity sponsored by API Project No. 44 under the direction of F. D. Rossini has led to a tabular presentation of compressibility factors for gaseous methane in the superheated region.la These values were obtained by using smoothed residuals, developed from experimental data and the Ben edict-Webb-Rubin equation,' to establish internally consistent compressibility factors for the real gas. This approach offers probably the most precise presentation of PVT information for gaseous methane to date. The PVT behavior of methane also can be represented accurately by equations of state such as the Benedict-Webb-Rubin equation or the Martin-Hoa equation.' These equations, although valuable for analytically representing the PVT behavior, lack precision in the compressed-liquid region. The PVT data available in the literature for methane were utilized to obtain reduced densities that are presented in Figs. 1 and 2 as functions of reduced temperature with parameters of reduced pressure. The critical constants utilized for methane are those reported by Kobe and Lynn." These critical values, T, = 191.1°K, PC = 45.8 atm and p, = 0.162 gm/cc, produce a critical compressibility factor, z, = 0.289. The PVT information, as obtained from the literature, in many instances was not presented in an immediately

By E. B. Hunt, T. M. Geffen, C. R. Stewart, V. J. Berry

Laboratory data .show that the gas-oil ratio performance of non-uniform porosity limestones produced by solution gas drive is sensitive to producing rate and to fluid properties. Nan-uniform porosity limestones are those for which laboratcry solution and external gas drive tests yield considerably different relative permeability ratio characterislics. The oil recovery performance by solution gas drive depends directly on the number of gas bubbles formed. Laboratory rates of pressure decline, which are 100 to 10,000 times greater than normal field rates. cause the formation of an unusually large number of gas bubbles. This results in abnormally high oil recovery efficiencies. Since it is impractical to reproduce the number of bubbles formed under field conditions, laboratory solution gas drive data on non-uniform porosity limestones are therefore not directly applicable to field operations. However, certain laboratory data can be used to make a conservative estimate of field performance. The concepts presented in this paper indicate the possibility that increased field oil recoveries may he obtained from non-uniform porosity limestones by rap- idly reducing reservoir pressure for a short interval of time. It has yet to be established that significant improvements in oil recovery from such reservoirs can he realized by varying the pressure decline rate within limits possible in the field. However, the possibility that recovery may be increased in this manner warrants further study. INTRODUCTION Limestone reservoir rocks can be divided into two general classes according to the nature of their pore space. One class has a comparatively uniform pore system composed mainly of voids between grains of the rock (intergranular type porosity). The other class, in addition to having intergranular porosity, has a secondary pore system composed of combinations of solution cavities, fractures, etc., and is considered to have non-uniform porosity. It has been shown in a previous publication' that the laboratory measured gas-oil flow behavior of uniform type porosity limestones is essentially the same for a test simulating a field solution gas drive as for one simulating an external gas drive. For non-uniform porosity type limestones a significant difference in gas-oil flow behavior was found. More efficient behavior was

Jan 1, 1955

By C. R. Sandberg, N. Stein, L. D. Wiener, E. B. Elfrink

this paper presents the results of a study of the slip velocity of gases rising through liquids in vertical tubes, inclined tubes, and vertical annuli. The data were obtained in gas-liquid systems which included combinations of air, propane and natural gas (over 97 per cent methane) with water, lubricating oils, and crude oils. The 214 data points obtained in this study along with 11 data points reported in the literature are incorporated in an empirical correlation which relates the mean slip velocity of gases flowing through liquids with the parameters gas rate, tube size, ratio of liquid viscosity to liquid density, gas density, liquid density, and the angle of the tube from the vertical. The average numerical deviation of the measured slip velocity data from values obtained from the correlation is 9.2 per cent. The average algebraic deviation between measured and Predicted data is 0.31 per cent, an indication that the correlation is a satisfactory representation of the data. The correlation presented in this paper will be useful primarily in the design of subsurface gas-oil separation equipment for increasing the efficiency of oil-well pumping installations. but may perhaps be extended to other situations of gases rising through liquid columns. INTRODUCTION The production of oil by pumping is often complicated by the presence of free and dissolved gas in the oil at bottom-hole conditions. This gas can be drawn into the pump barrel and result in "vapor locking" the pump, thus reducing the efficiency of the pumping operation. Large amounts of free gas may be excluded from the pump by the use of a gas-anchor, which is designed to allow the separation of free gas from the oil before it is drawn into the pump. The design of a suitable gas-anchor depends on a knowledge of the difference between the velocity of the free gas and that of the oil as the oil travels down toward the pump intake. This rate of gas rise through the oil is termed "slip velocity." AS a first step toward the full understanding of this roblem, the mean slip velocities of gases flowing through static

Jan 1, 1952

By A. L. Lee, M. H. Gonzalez, R. F. Bukacek

Experimental viscosity data for methane are presented for temperatures from 100 to 340F and pressures from 200 to 8,000 psia. A summary is given of the available data for methane, and a comparison is presented for that data common to the experimental range reported in this paper. Correlation of the data is presented, and predicted values are given for temperatures up to 460F and pressures up to 10,000 psia. INTRODUCTION The increasing ranges of temperature and pressure at which fluids are produced, transferred and processed in the petroleum and chemical industries stress the need for accurate information on physical properties, both for engineering calculation and for improving the methods used to estimate physical properties. A survey of the literature reveals that disagreements between published data on the viscosity of methane are common and that most investigations have been conducted over restricted temperature and pressure ranges. This situation made it desirable to undertake an investigation with a scope common to most of the literature data available. APPARATUS AND MATERIALS The apparatus used in this investigation was the capillary tube viscometer'described by Dolan,10 with modifications introduced in the general design and operation of the instrument. It has an effective pressure range from 14.7 to 10,000 psia and a temperature range from room temperature to 400F. The design of the viscometer is based on the establishment of a manometric head between two vessels containing the test fluid and a volume of mercury. The reservoirs are connected by a capillary tube through which the test fluid flows and a tube through which mercury flows. A pressure gradient is established by elevating one of the vessels above the other, the resulting flow of mercury displacing the test fluid through the capillary. The schematic diagram of the system (Fig. 1) shows the arrangement of the equipment auxiliary to the viscometer. The density cells assembly (E, Fig. 1) was not used in this investigation since reliable data are available on methane. 23 The assembly consists of a bank of eight 316 stainless steel pycnometers. Experimental values of the density of a natural gas19 and i-butane14 samples have been obtained for pressures up to 8,000 psia and temperatures up to 340F. The methane used was obtained from the Southern California Gas Co. Mass spectrometric analysis showed a composition of 99.8 percent CH4, 0.1 percent C2H6 and 0.1 percent C3H8. EXPERIMENTAL VALUES AND ACCURACY OF DATA Determinations were carried out at several temperatures. With the temperature in the cabinet controlled at a given level, a series of runs were made at various pressures. To establish pressures lower than that obtained from the methane cylinder, it was necessary to vent the gas slowly to the atmosphere until the desired pressure was reached. The consistency of the test run time for a particular condition can be used as a criterion for the accuracy of the data. A further check is obtained by comparison of the viscosity values obtained using different driving heads. To establish the instrument's reproducibility, a large number of runs were made at three completely different conditions of temperature and pressure. Data obtained were analyzed statistically to determine the most reliable value for the variance of the experimental data. It was found that the variances of the data taken at different experimental conditions were not statistically equivalent, but they were small. The largest variance was found at 8,000 psia and 100F, where the standard deviation was 0.385 micropoise. Therefore, 99.9 percent of the data can be expected to fall within 1.15 micropoises of the population mean, except for the intervention of systematic errors in operation.

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.' ' Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.'" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By M. R. Tek, L. K. Thomas, D. L. Katz

Threshold displacement pressures are needed to determine how much overpressure can be used in storing natural gas. An experimental technique for determining threshold pressures by displacing water with gas from samples saturated with water is presented. Threshold pressures for eight low permeability samples were measured. Threshold pressure data obtained in this work, plus data on higher permeability samples reported in the literature, are correlated with porosity, permeability, surface tension and formation resistivity factor. Mercury injection pressures also were measured and correlated with air-water threshold pressures. Threshold pressures were found to be independent of time. An aerial photograph of gas emerging from the top of a core shortly alter its threshold pressure has been exceeded shows that the gas bubbles are uniformly distributed across the face of the core. Channeling will occur, however, when an increased gas phase permeability is reached, A porous medium can be re resealed alter its threshold pressure has been reached, provided it has not been desaturated below a fixed saturation. Its new threshold pressure will be lower than when the sample is 100 percent saturated with water. INTRODUCTION The degree of overpressure in excess of the discovery pressure that a gas storage reservoir can withstand is determined by the ability of its caprock to contain gas, providing the pressure does not exceed the structural limit of the reservoir. The retention of gas in a reservoir by a caprock saturated with water is a result of the capillary forces acting at the gas water interface. Without the presence of water in a caprock, gas would leak out of the reservoir at a rate determined by the permeability of the caprock to gas. The ability of a caprock to contain gas is expressed in terms of its threshold displacement pressure, which is defined as the minimum pressure needed to initiate the displacement of a wetting phase by a nonwetting phase from a porous medium 100 percent saturated with the wetting phase. The gas industry must be able to predict how much "overpressure" a gas reservoir can withstand before leakage occurs by gas displacing water from caprock overlying the reservoir. This is of economic importance since the storage capacity of a gas reservoir is proportional to its pressure. The discovery pressure for most reservoirs is found between the upper and lower limits1 shown in Fig. 1. Discovery pressures in excess of their hydraulic pressure gradient are attributed to the compaction of an isolated shale formation by the overburden.2.3 The lower limit of 0.433 psi/ft corresponds to the weight of the overburden, the gradient of the rock itself. Somewhere between

Jan 1, 1969

By F. G. Miller, R. C. McFarlane, T. D. Mueller

During the process of gas storage in pressure-depleted oil reservoirs, it has been observed that in some instances additional liquid oil is recovered and that the composition of the storage gas is materially altered. A mathematical study was made of the dynamic behavior of such a depleted oil reservoir undergoing gas injection. The important variable considered in this study, not included in previously published work, was that of compositional effects on the phase behavior of two-phase flow. Pressure, saturation and component composition profiles were developed for a linear, horizontal and homogeneous porous medium containing oil and gas but undergoing dry gas injection. Special new techniques were developed to overcome the problems of numerical smoothing which arise in the solution of the equations representing such systems. The method of solution includes the development of partial differential equations describing the behavior of the system, representing these equations by finite difference approximations, making certain simplifying assumptions and, finally, applying methods of numerical analysis with the aid of a high-speed digital computer. In an example calculation, results using the mathematical model are compared with field observations made on a gas storage project in Clay County, Tex. This field project involved a depleted oil reservoir used 'for gas storage and gas cycling purposes. As a result of these processes, the reservoir yielded substantial amounts of secondary oil, both in the form of stock tank oil and as vaporized products in the produced gas. The methods derived in this study may be applied to a variety of oil reservoir problems which are dependent on compositional effects. In recent years the number of oil reservoirs being used for gas storage purposes has increased greatly, and there has been at least one published account of additional oil recovery resulting from gas cycling a depleted oil reservoir after repressur-ing with dry gas for storage purposes. Additional oil recovery from oil reservoirs resulting from gas storage operations could become an important secondary recovery process. This is especially true since the use of natural gas in large metropolitan areas continues to increase and more gas storage volume near these areas is needed. These facts provided the motivation for the work reported here. This paper reports on a study of the inter-relations of composition, saturation and pressure changes which occur when hydrocarbon gas is injected into an oil reservoir system. From an understanding of the process, prediction methods may be developed for use in forecasting the secondary recovery products from gas storage operations in oil reservoirs and, consequently, .the economics of such projects can be developed. The process of storing hydrocarbon gas in an oil reservoir and the cycling process both involve the transient multi-phase flow of oil reservoir fluids with changing fluid compositions. There is a scarcity of work published on the subject. One of the first laboratory models developed to study a problem similar to the one posed here was published in 1948 by Standing, Lindblad and Parsons. l They developed a model to study gas condensate recovery by gas cycling, including the effect of heterogeneous reservoir permeability. Their results showed the compositional effects of injecting a gas of known composition into a two-phase condensate reservoir of known composition. A major finding of their study indicated that maximum oil recoveries from gas cycling projects result when low cycling pressures are used. Previous to their work it was generally considered