Search Documents

By Lee Hillard Meltzer, Ralph E. Davis

INTRODUCTION During the past few years emphasis has been placed upon methods of estimating the future expectancy of gas production from natural gas fields. Before technical methods were applied, the production expectancy over future years was based upon the knowledge of gas well behavior, learned through long experience and embedded in the "know-how" of men long in the gas producing business. It is doubtful that a technical study of future expectancy of a gas field or a group of fields was ever prepared for the preliminary planning of a natural gas pipe line system built prior to about five years ago. The decline in well production capacity was naturally recognized by all familiar with the business since its earliest beginnings more than 75 years ago. In 1953, the Bureau of Mines published Monograph Number 7, "Back-Pressure Data on Natural Gas Wells and Their Application to Production Practices," which gave to the industry the first technical analysis of the decline in production of individual gas wells. This method affords a means of estimating the future production in relation to decline in reservoir pressure. The demand for technical determination of expectancy of future gas productivity from fields or a group of fields led technical men to the application of the knowledge of well behavior to the problems. The decline in a well's ability to produce as pressures declined could be estimated by the use of the curve known as the "back-pressure potential curve" as developed by the Bureau of Mines. A field containing few, or even numerous, wells could be analyzed on the basis of the sum of potentials of all wells. In most studies of this nature, the problem is to estimate the rate of production that can be expected, not only from present wells but also, from wells that will in the future have to be drilled into the reservoir being studied. The "back-pressure potential" method requires that the following data be known or estimated: (1) Proved gas reserves. (2) Current shut-in pressures and rate at which shut-in pressures change with production. (3) Back pressure potential data on wells in the source of supply. (4) Ultimate number of wells which will supply gas, and their potential. (5) Limitations on productivity such as line pressures against which the wells will produce, friction drop in the producing string, and so forth. It is evident that the resulting estimate of gas available in each year for a future of say, 20 years, contains many uncertainties. While the method may have considerable merit for a field that is fully developed, it cannot be completely dependable in fields that are only partially developed. In such cases, some of the data upon which it is based can only be estimated or assumed. In the study of this problem during the past few years, a method has been developed which we believe has great merit, especially when applied to fields subject to substantial future drilling, and when applied to the study of fields which, on the average, appear to have characteristics similar, in general, to the average of the fields used in the development of the "yardstick" outlined herein. From an analysis of the production history of 49 reservoirs which are depleted, or nearly depleted, a curve has been constructed which shows the average performance of the reservoirs during the declining stages of production. When properly applied, this "average performance curve" can be used to determine the stage of depletion at which a reservoir or group of reservoirs will no longer be able to yield a given percentage of the original reserves. "AVAILABILITY" AND "AVAILABILITY STUDIES" The rate at which. a reservoir will yield its gas depends basically upon physical factors, such as the thickness and permeability of the sand, the effect of water drive, if any, and other conditions, and upon economic factors, such as the number of wells drilled. Within the ranges set by the physical conditions, a rate of delivery tends finally to become established. The rate (or range of rates) represents a balance between the interests of the operator, who desires the maximum return from his property and of the pipe line owner, who desires to maintain a firm supply for his market. This balance, which is influenced by the terms of the contract, determines the capacity which will be developed by the operator, and the time and rate at which the decline in production is permitted to occur. Thus the "availability" of gas

Jan 1, 1953

By J. S. Aronofsky, R. Jenkins

A simple means of predicting the flowing well pressure history in a natural gas reservoir has been developed. The differential equation for unsteady radial flow of gases through porous media was solved by numerical methods using electronic punch card machines. The results obtained from these calculations are presented in a generalized form which is believed to be most convenient for engineering use. These results are compared with those presented in a paper by Bruce, et al4. An effective radius of drainage has been defined which permits the equation of steady state gas flow to be used easily to predict well pressures in gas reservoir depletion. An equivalent effective radius of drainage is defined for liquid flow systems and a comparison is made of the gas solution to the liquid solutions. INTRODUCTION This paper presents a numerical method for describing the transient flow of gases radially through a porous medium from which the production rate is specified. The results obtained from these calculations are presented in a generalized form which is believed to be the most convenient for engineering use. The principal application of these calculations should be to flow in natural gas reservoirs, in the interpretation of pressure drawdown tests, in the estimation of reserves, and in calculation of injection pressures required in gas injection and cycling operations. The evaluation of natural gas wells is based upon back-pressure tests in which measurements of production rate and flowing well pressure are made for periods of one, two, or more days. The flow is assumed to have stabilized at this time, and the observed values of pressure and production rate are used with standardized formulae to determine the well's allowable production rate. Furthermore, the results obtained from such a test of short time duration are extrapolated in such a manner to allow a prediction of what the well pressure will be at some future time for a specified (but possibly untested) flow rate. The manner of conducting such back-pressure tests suggests several questions as follows: 1. 'What is the definition of stabilized flow and what factors determine the length of time required for the flow to become stabilized? 2. How may the recorded values of pressure and rate be utilized? 3. As the well is produced and the reservoir pressure falls, how will the flowing well pressure change? Can this flowing well pressure be predicted from the back-pressure test? 4. Are the results obtainable from the test influenced by a particular procedure used in obtaining the results? 5. Is it possible to estimate the gross reservoir permeability from back-pressure data? It is not within the scope of this paper to provide solutions Or even to discuss adequately all the above questions. However. in considering such matters, the

Jan 1, 1955

By K. H. Coats

This paper describes an approximate technique for handling the problem of percolation of evolved gas upwards through the oil column in computer simulation of natural depletion. This technique has been incorporated into a general model for simulating three-phase flow in one, two or three dimensions. The mathematical model for perfonning these calculations is described in detail in a separate article. 1 While vertical gas percolation occurs during natural depletion in a reservoir of any configuration, it is especially pronounced in the pinnacle reel or bioherm. The reel may have an areal extent less than one section with a thickness of up to 800 It. The upward percolation of evolved gas often limits the time step in calculations as little as one day. Using this small time step, computer expense for a single 30- to 40-year simulation on even a one-dimensional basis has exceeded 92,000. The method described here for handling the percolation allows time steps of 60 days or more (depending on reservoir size and production rate), resulting in a considerable reduction in computing expense. The method also allows calculations in which secondary gas caps build up in tight zones of the oil column below the main gas cap. The validity of the method is indicated in connection with an example pinnacle reef field. Calculated results using a 2-day time step (where the method in question is not needed and not invoked) are compared with results using the method and a 60-day time step. The comparison shows good agreement. Results from one- and two-dimensional simulations of the reef are presented along with the corresponding computing times on the CDC 6600 computer. A three-dimensional simulation was also performed and the required computing time is given. INTRODUCTION A general analysis for simulating three-dimensional, three-phase flow in reservoirs has been developed. 1 In applying this model to reservoirs undergoing natural depletion, a time-step restriction was encountered due to the flow of evolved gas upwards through the oil column and toward the gas cap. A remedy for this problem has been incorporated in the analysis and is described here. The time-step restriction is encountered only in calculations which include flow in the vertical or near vertical direction. The restriction occurs to some extent in reservoirs of any geometrical configuration, but it is especially pronounced in the pinnacle reef or bioherm where the ratio of thickness to areal extent is unusually large. This paper describes the problem, a method of handling it, incorporation of the method in the three-dimensional, three-phase model, and, finally, a test of the method's validity in an application to an example pinnacle reef reservoir. THE PROBLEM During the early stages of reservoir depletion, pressure falls below bubble point in progressively lower regions of the oil column. Continuing production results in evolution of dissoived gas throughout the oil column. This gas then percolates upwards toward the top of the reservoir. If the gas encounters a sufficiently low permeability zone in its travel upward, then it can accumulate and form a secondary gas cap. If no zones are encountered which are tight enough to hold gas against the gravity forces, then the gas travels on until it reaches the top of the sand or the main gas cap. Numerical simulation of multiphase flow in reservoirs is generally performed with time steps such that the flow of a fluid into or out of a block in one time step is a fraction (considerably less than one) of the total amount of that fluid present in the block. Use of a time step where the time step's

Jan 1, 1969

By G. Rowan, M. W. Clegg

Approximate analytical solutions for non-Darcy radial gas flow are derived for bounded and infinite reservoirs producing at either constant rate or constant pressure. These analytical solutions are compared with published results for non-Darcy flow obtained on digital and analogue computers, and the agreement is show to be very good. Some observations on the interpretation of gas well tests are made. INTRODUCTION The flow of gases in porous media is a problem that has been the subject of much study in recent years [l - 71 , and many methods have been proposed for solving the non-linear equations associated with it. The assumption that the flow satisfies Darcy's Law It has been observed, however, that the linear relationship between the flow rate and pressure gradient is only approximately valid even at low flow rates, and that as the flow rate increases the deviations from linearity also increase. It has been suggested by a number of authors that Darcy's Law should be replaced by a quadratic flow law of the form This form of equation was first suggested by Forchheimer [9] and, later, Katz and Cornell [10], and Irmay [lll, developed a similar equation. Houpeurt [12] derived this form of equation using the concept of an idealised pore system in which each channel consists of sequences of truncated cones giving rise to successive restrictive orifices along the channel. This type of representation leads to a quadratic flow law of type [4], for all fluids, but it is found that the quadratic term is only significant in the case of gas flow. The methods of Houpeurt for solving gas flow problems will be discussed further in another section of this paper. Solutions of the non-linear equation for Darcy gas flow 121 may be classified as either computer (digital and analogue), or approximate analytical ones. The former include the well-known solutions of Bruce et al. [13], and Aronofsky and Jenkins [14] , but the latter solutions, apart from the simple linearisation of equation (2) to yield a diffusion equation in p2 , are not so well-known. Amongst these approximate analytical solutions should be mentioned the work of Barenblatt [IS], who has applied the techniques of dimensional analysis to obtain similarity solution to a large class of non-linear equations; also, a number of Russian authors, including Barenblatt [16] , Kim [17] and others, who have used the method of moments for obtain-

Jan 1, 1965

By J. Husson, R. Iffly, M. Gondouin

A systematic variation of well deliverability, as reflected from isochronal back-pressure tests performed at regular intervals, has been observed in some gas condensate wells producing at high rates. The same effects have been obtained using a numerical model of gas and condensate flow which takes into account secondary gasoline deposited in the pore space as a result of pressure reduction, and non - Darcy flow of gas in the vicinity of the wells. Matching calculated values with previous test results has been possible, and future predictions have been obtained. An application of this method to the Hassi Er R'Mel gas-condensate field in Algeria is tentatively shown. INTRODUCTION Flow capacity of gas wells is generally derived from an analysis of back-pressure tests. The empirical equation q = C (?p2)n used by Rawlins and Schellhardtl can be derived rigorously assuming that steady-state radial flow of a dry gas of constant viscosity and compressibility is established during each flow period of the well tests. Furthermore, when Darcy's law applies in the entire flow region, the theory predicts that the exponent n is equal to 1. In low-permeability reservoirs, it was soon discovered that the time required to reach a stabilized flow often exceeded the duration of the flow periods normally available for testing wells. Consequently, transient gas flow had to be considered instead of the steady-state assumption previously used. This led to the isochronal testing procedure established by Cullender2 which has largely replaced conventional back-pressure testing. For dry gas fields, this method yields definite values of C and n equivalent to those of the empirical equation. These values should remain constant for each well as long as the permeability of the formation and the characteristics of the gas (viscosity and compressibility) do not change appreciably. This is the case when reservoir pressure remains -close to the original value and when the formation near the wellbore remains free of plugging. Under those conditions, stabilized flow potential curves of gas wells can be established from a single sequence of isochronal flow and shut-in periods. An analysis of a pressure build-up following a longer production period provides additional data on the transmissivity (k h/ u) of the reservoir, and eventually on the 'drainage radius rd of the well, which can be related to the value of C so that future performance of the well can be predicted using the concepts developed by A. Houpeurt.3 At high flow rates, Darcy's law no longer applies in the vicinity of the wellbore, and inertial effects in the high velocity gas flow introduce additional pressure drops. As a consequence, exponent n of the backpressure tests becomes smaller than 1, and a slight curvature of the log-log plot of ?p2 vs q can be predicted when going from very low to high rates of flow (Elenbaas and Katz).5 The effects of variations of viscosity and compressibility with pressure on the radial flow of dry gas in an infinite reservoir were taken into account by Jenkins and Aronofsky.7 Numerical solutions of the transient flow of an ideal gas in finite radial reservoirs were presented by Bruce, Peaceman, Rachford and Rice.8 In the case of gas condensate wells, however, the presence of gasoline in the pore space as soon as reservoir pressure is reduced below the dewpoint pressure further complicates the interpretation of flow tests so that the prediction of stabilized well performance becomes very difficult. Field observation shows that both C and n derived from isochronal tests vary in time, even when reservoir pressure has not changed appreciably. Such a variation cannot be attributed to any change of the gas characteristics, and must result from the effect of a gasoline saturation on the gas flow.

By R. G. Turner, M. G. Hubbard, A. E. Dukler

Gas phase hydrocarbons produced from underground reservoirs will, in many instances, have liquid phase material associated with them, the presence of which can affect the flowing characteristics of the well. Liquids can come from condensation of hydrocarbon gas (condensate) or from interstitial water in the reservoir matrix. In either case, the higher density liquid phase, being essentially discontinuous, must be transported to the surface by the gas. In the event the gas phase does not provide sufficient transport energy to lift the liquids out of the well, the liquids will accumulate in the wellbore. The accumulation of the liquid will impose an additional back pressure on the formation that can significantly affect the production capacity of the well. In low pressure wells the liquid may completely kill the well; and in the higher pressure wells there can occur a variable degree of slugging or churning of the liquids, which can affect calculations used in routine well tests. Specifically, the calculated bottom-hole pressures used in multirate backpressure tests will be erroneous if the well is not removing liquids on a continuous basis, and gas:liquid ratios observed during such a test may not be correct. Several authors1,3,8,14 have suggested methods to determine if the flow rate of a well is sufficient to remove liquid phase material. Vitter14 and Duggan' proposed that wellhead velocities observed in the field would be adequate for keeping wells unloaded. Jones8 and Dukler3 presented analytical treatments resulting in equations for calculating, from physical properties, the minimum necessary flow rate. An analysis of these studies indicates the existence of two proposed physical models for the removal of gas well liquids: (1) liquid film movement along the walls of the pipe and (2) liquid droplets entrained in the high velocity gas core. Although there probably is a continuous exchange of liquid between the gas core and the film, they will be treated separately for the purposes of this study. The development and comparison of these separate models with experimental data will permit the determination of which, if either, is the controlling mechanism for the removal of liquids from gas wells. The Continuous Film Model Liquid phase accumulation on the walls of a conduit during two-phase gas/liquid flow is inevitable due to the impingement of entrained liquid drops and the condensation of vapors. The movement of the liquid on the wall is therefore of interest in the analysis of liquid removal from gas wells. If the annular liquid film must be moved upward along the walls in order to keep a gas well from loading up, then the minimum gas flow rate necessary to accomplish this is of primary interest. The analysis technique used follows Dukler2 and Hewitt5 and involves describing the profile of the velocity of a liquid film moving upward on the inside of a tube. The minimum rate of gas flow required to move the film upward is then calculated.

Jan 1, 1970

By R. Al-Hussainy, H. J. Ramey

Previous gas well test analyses have been based mainly upon linearizations of ideal gas flow results, although a method for drawdown analysis based upon real gas flow results has been proposed. Linearizations based upon ideal gas flow require estimation of gas physical properties at Tome sort of average pressure, and implicitly involve the assumption that pressure gradients are small everywhere in the reservoir. A new real gas flow equation has been developed by means of a substitution which couples presTure, viscosity and gas law deviation factor. This sub-stitution has been called the real gas pseudo-pressure. Use of the real gas pseudo-pressure leads to simple equations describing real gas flow which do not contain pressure-dependent gas properties, and which do not require the assumption of small pressure gradients everywhere in the flow system. Equations required to determine flow capacity, well condition and static formation pressure from pressure drawdown and build-up tests with the real gas pseudo-pressure concept are presented. Also shown are applications of the real gas pseudo-pressure to back-pressure testing, the gas materials balance and rigorous determination of average gas properties for previous gas flow equations. Included are sample calculations for well test analyses. INTRODUCTION A recent study' of the flow of the transient flow of real gases in ideal radial systems showed that it is possible to consider gas physical property dependence upon pressure by means of a substitution called the real gas pseudo-pressure. The substitution has the advantage that it enables engineering solutions for steady and transient flow of real gases that are more accurate and general than those previously available. The solutions are particularly important for the case of gas flow in tight, high-pressure formations under large drawdowns. Even for more ideal flow conditions, the real gas pseudo-pressure concept leads to important rules for finding average gas physical properties pertinent for older well test analyses. Ref. I provides a detailed study of the theory behind the real gas flow correlations to be used in this paper. The purpose of this paper is to provide necessary engineering forms for use of the real gas flow results, and to illustrate applications. The following considers flow of real gases in an ideal radial flow system. It is assumed that: (1) formation thickness, porosity, water saturation, absolute permeability, temperature and gas composition are constant; (2) gas compressibility and density are functions of pressure as described by the gas law pv = znRT; (3) gas viscosrty is a function of pressure; (4) effective permeability to gas may be a function of pressure to account for liquid condensation; (5) condensate is immobile; and (6) the formation has no dip.

Jan 1, 1967

By P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By W. Hurst, R. E. Leeser, W. C. Goodson

Three aspects of gas deliverability are presented in this paper. The first treats with the gas deliverability or availability of a normal depletion-type dry gas field. Such encompasses not only the period of stabilized constant rate, but more so, the "tailings" when a fixed abandonment pressure is reached and the rate by necessity must decline. A comprehensive work plot is offered, developed from mathematics herein included, that removes the triai-and-errnr computations that attended such undertakings in the past. The second part treats with the discount factor of the open flow potential constant from what is observed initially in testing a gas well to what is evidenced when stabilization is reached. This prevails in tight formations, such as the Kansas Hugoton field which is offered as the example. The means of establishing this factor are pressure build-up curves which, as sustained by analytical deductions, reproduce this entire period of transient flow under conditions of a constant rate influx. Finally, what is offered in this paper is the deliverability performance of an exceedingly rich gas condensate field producing from a tight formation. The example shown is the Knox Bromide field in Oklahoma, producing from the Bromide formations. The results are ominous, showing early reduction in permeability to gas pow, due to the retrograde condensate forming in the pore space, with the attending early logging-up of these wells. The analytics of lowered permeability are incorporated in the gas deliverability formula along with the PVT data that gives the increased condensate liquid saturation as the gas flows to the well bore. This paper would not be complete without a critique oflered at the end. With the many gas wells now in production and those that have completed their life, there has been no factual information collected by any source as to what constitutes that permeability range where a gas well would be unimpaired in its gas deliverability by the presence of rich condensate content, and the lowered range where such would be harmful. This question confronts all producers. INTRODUCTION Various aspects of gas deliverability are presented in this paper that includes depletion-type reservoirs, deteriora- tion factor of the gas deliverability constant, and the performance of a rich gas condensate reservoir producing from a tight sand. With respect to the presentation of gas deliverability and its tailings for depletion-type gas reservoirs, one notes that this is essentially the information offered by every gas transmission company and producer appearing before the Federal Power Commission for Letters of Conveyance in the dedication of reserves. In the ordinary procedure, as many engage upon this study, trial-and-error calculations are included, particularly as apply to the tailings. For many years one of the writers has employed mathematical analyses to perform this step and avoid the complexities so associated. In the preparation of this paper these analyses have been amplified to include any slope n for the open flow potential relationship for which the tailings can be determined from Fig. 1. With reference to the deterioration or discount factor of the open flow potential constant as such occurs in the gas deliverability formula, this for the most part has been an unexplored subject. Although the issue first appeared in the Kansas Hugoton field, where such was surmised but only recently resolved, this situation of a deterioration of the gas deliverability constant can occur wherever dry gas production from a tight sand is encountered. The first concerted attacks upon this problem were the presentations of Hurst' and Goodson' before the Kansas Corporation Commission to show that transient fluid flow and unsteady-state flow formulas prevailed. This was amplified later before the Federal Power Commission3 to show that this deterioration factor could be identified from pressure build-up curves. This has been reported by McMahon.4 Its importance to the industry merits the review of these essential features in completing the program on the aspects of gas deliverability. Finally, as illustrated here, for a low permeability formation such as the Knox Bromide field where the gas is rich, representing some 165 bbl of condensate per MMcf of effluent gas, the gas deliverability can be of limited extent in the life of the field, leaving substantial amounts of condensate and gas unrecovered. In cases such as this, gas cycling is mandatory. This is particularly revealed by the fluid mechanics introduced here, employing factual field as well as laboratory data, to show this-restriction upon gas deliverability. PRESSURE DEPLETION What will now be offered is the study of gas deliver-

By A. J. Garnier, N. H. van Lingen

Rock downhole is known to be lesc. drillable than when brought to the surface. This must be ascribed mainly to the presence under downhole conditions of a pressure differential across already made chips, which hinders their being lifted. The pressure differential has partly a static and partly a dynamic origin. Balling-up of bits is another consequence of this pressure differerztial. Reduction in penetration rate owing to an increase in the strength of the rock is governed by the difference between the mud pressure and the pressure of the for/nation pore liquid. Rotationally symmetric geostatic stresses have 170 effect on drillability. Fracture of rock when drilling will be brille in most cases. The above is supported by laboratory drilling experiments with drag bill and roller bits an elevated mud, pore, and confining pressures On rocks differing in strength and permeability. INTRODUCTION In oil well drilling drillability of rock is found to decrease with increasing depth of the hole. Naturally deep rock will be more compacted and, therefore, harder to drill than shallow rock of the same type. However, apart from this the drillability of a sample of deep and compacted rock brought to the surface is generally much higher than in its original location downhole. In view of the economic implications of this reduction in drillability, it seems worthwhile to analyze its causes. The origin clearly has to be sought in the difference of environment. The only conceivable factors would seem to be the presence of mud under pressure, the pressure of the formation pore liquid, and the overburden of the rock. Down the borehole the rock is compressed triaxially by mud pressure and overburden. It is well known that the strength of rock is increased when confined by external pressure1. Various authors have, therefore, ascribed the difference in drillability mainly to the strengthening of rock by triaxial compression. Another factor, mentioned by Bobo and Hoch5, is that forces. including "pressure differential forces", tend to hold a dislodged particle in place. However, the conditions determining their magnitude are not clarified nor is their effect on drilling rate assessed. In the laboratory, drilling experiments on pressurized samples of rock yielded evidence that pressure differential forces holding the chips down are the major factor in reducing rate at depth. This paper describes the experiments and shows how the results enable both a qualitative and a quantitative interpretation of the factors determining the effect of chip hold-down on drilling rate. The implication with respect to balling-up and jet action is also discussed. The greatcr part of the experiments were performed in the pressure vessel shown in Fig. 1. The spacc between the sample of rock and the vesscl is divided by "0" rings into three separate chambers. The rock sample is confined laterally by pressurized oil in thc middle chamber. Penetration of oil into the sample is prevented by partly jacketing it with brass foil. Thc pressure of the drilling fluid in the hole can be adjusted via the upper chamber. Permeable rock specimens arc water saturated before being drilled. With a properly plastering mud as drilling fluid the pressure of the pore water can be adjusted independently via the lower chamber.

By J. G. Richardson, K. H. Coats

During the initial growth of a gas bubble in an aquifer storage reservoir the injected gas tends to override the water. The resulting low displacement efficiency and high rate of gas travel down-structure make it difficult to maintain water-free production. This paper illustrates the use of two-dimensional, two-phase calculations to simulate this gas-water displacement. Calculations were performed for a stratified aquifer and for several homogeneous sands differing in permeability, porosity and thickness. Results are discussed and compared with performance predictions obtained from the simpler Buckley-Leuerett and Dietz formulas. This comparison indicates the necessity for the two-dimensional calculations in realistically simulating the gravity override of the water by injected gas. INTRODUCTION In recent years natural gas has been stored near markets in aquifers where insufficient storage capacity is available in depleted fields. Operators of these aquifer storage reservoirs have encountered technical problems relating to the gas-water displacement accompanying initial growth of the gas bubble. Injected gas tends to override the water, with a resultant low displacement efficiency and high rate of gas travel down-structure toward spill points. Low displacement efficiency makes if difficult to sustain water-free gas production. Sometimes the fingering of gas down-structure may be so pronounced that the injection rate must be severely curtailed, which excessively lengthens the time required for bubble growth. This paper is concerned with estimating — for given aquifer characteristics and fluid properties — the displacement efficiency, rate of gas movement down-structure and rate of gravity drainage of water behind the gas front. Several published articles relate to simulative capability necessary to handle this problem. The Dietz formula1 and Buckley-Leverett2 method have some applicability to the problem; Woods and comer3 reported one-dimensional, radial calculations which accounted for the two-phase flow of gas and water. Most recent research directed toward understanding aquifer behavior in relation to gas storage has concentrated on single-phase flow in aquifer.4 Douglas, Peaceman and Rachford5 presented a method for calculating multidimensional two-phase flow in reservoirs; their method has been applied in work reported by Nielsen and Tek,6 Blair and peaceman7 and Goddin.8 Refs. 5 and 7 show comparisons between experimental data and calculations similar to those employed here. One purpose of this paper is to illustrate the use of two-dimensional, two-phase flow calculations in simulating the difficult problem of water displacement by gas in a vertical cross-section. The pronounced fluid density and viscosity differences combine with the disparity in horizontal and vertical dimensions of the cross-section to pose one of the more difficult problems in reservoir simulation. The calculations described account for capillary and gravity forces, relative permeability and reservoir heterogeneity. An example reservoir is described and calculated results are presented for a variety of injection rates and values of permeability, reservoir thickness and dip angle. A second purpose is to compare the two-dimensional calculations with results obtained from the Buckley-Leverett and Dietz formulas. For this reason most calculations were performed for a simplified example reservoir to which the latter techniques might reasonably be applied. However, the computer program employed applies equally well to cases involving arbitrary spatial variations of permeability, porosity, reservoir thickness and dip angle.

By J. M. Campbell, T. S. Buxton

The most widely used methods of predicting the volumetric properties of gas are based on the principle of corresponding states, which asserts that the compressibility factor is a universal function of reduced temperature and pressure. Previous studies have shown that the acentric factor, as proposed by Pitzer, l is an important addition to reduced pressure and reduced temperature as factors affecting the compressibility factor. Results of this study indicate that, if the pseudocritical temperature and pressure used to determine the reduced conditions are adequately predicted, characterization of natural gas-carbon dioxide mixtures with the acentric factor will allow reliable determination of the compressibility factor. Comparisons of predicted and experimental compressibility factors have shown that the pseudo-critical constant rules of Stewart, Burkhardt and Voo 2 are satisfactory for hydrocarbon mixtures. However, these rules fail to predict the pseduo-critical constants for hydrocarbon-carbon dioxide mixtures. Based on graphically determined pseudocritical temperatures for binary hydrocarbon-carbon dioxide mixtures, a correlation which gives the required correction to the Stewart, Burkhardt and Voo rules was prepared and a compressibility factor prediction technique was proposed. To test the proposed technique, compressibility factors for live mixtures of methane, carbon dioxide and either ethne or propane were experimentally determined at 100, 130 and 160F and pressures up to 7,026 psia. The predicted and experimental compressibility factors for these live mixtures had an average absolute deviation of 0.55 percent. INTRODUCTION Two-parameter generalized correlations based on the principle of corresponding states (PCS) are presently available for reliably predicting the compressibility factor of lean natural gas, i.e., natural gas with low concentrations of hydrocarbons heavier than methane. When these same correlations are used for mixtures of natural gas and carbon dioxide, large deviations between actual and predicted compressibility factors are observed. This study was undertaken to provide a means for accurately predicting the compressibility factor for such mixtures. In the past, two avenues of approach have been employed to extend the applicability of the PCS. The first has been to introduce additional parameters. The second has been to develop combination rules for predicting pseudocritical constants which do not suffer from the limitations of the commonly used molal average constants. In this study, both of these approaches have been considered and employed to arrive at a method for predicting the compressibility factor of mixtures of hydrocarbons and carbon dioxide. SELECTION OF A THIRD PARAMETER Failure of the original PCS to predict the compressibility factor for all pure gases, regardless of their mass, shape or polar moment, has led to the introduction of additional parameters into the PCS. Parameters have been introduced to correct for quantum deviations,3-5 non-spherical or globular shapes6-9 and for polar moments.8-12 In addition to these parameters which are based on some microscopic property and intended to correct for some specific cause of deviation, more general parameters based on some bulk property have also been introduced.1913-16 To use additional parameters to correct for each of the possible causes of deviation among the individual constituents of a gas mixture would require a rather complex form of the PCS. To maintain the PCS in as simple a form as possible, it is more desirable to have a single third parameter which is based on some bulk property influenced by several of the factors causing deviations. The third parameter selected for use in this study is a factor originally proposed by Pitzer.9

By P. Pollard

A method has been developed for evaluating acid treatments in fractured limestone fields by breaking down pressure drawdown into three component parts: (I) pressure differcntial across "skin" near the borehole face, (2) pressure differential due to flow resistance in the coarse communicating fissures and (3) pressure differential between the fine voids and the coarse fissures. It is apparent that in most successful acid treatments the first term, skin resistance, has been reduced or eliminated. Further, it is often possible to estimate the volume of coarse fissures associated with the second term, coarse fissure flow resistance. In cases where this volume is comparable with practical acid volumes it seems likely that this resistance also may be attacked with a suitably retarded acid. INTROD UCTION Acid treatment has been successfully applied as a general practice in the limestone wells in the Mara/ Maracaibo districts of Venezuela for some 10 years. At the same time the need has been felt for a more precise method of: (1) evaluating the effect of an acid job and (2) selecting wells which should benefit from an acid treatment. The following method, based on an analysis of build-up curves for wells producing from fractured limestone, provides a means of forecasting the effect acid will have on a well and in many cases the probable order of production rate increase which may be expected. In certain cases it is also possible to calculate the fissure volume through which the acid has to be displaced for best results. Various approximations are introduced in the method, but within the limitations of these the analysis gives an indication of how a well, whether previously acidized or not, may be expected to respond to acid treatment. Suitable build-up curves cannot be obtained on "tight" wells which do not give sustained production and the method is not applicable to such wells. Figs. 1 through 7 illustrate various types of build-up curves obtained from limestone wells in Western Venezuela. Curves are included showing (1) both skin and coarse fissure resistance, (2) skin resistance only, (3 and 4) reduction of skin resistance by acidization, (5) fissure resistance only, (6) reduction of coarse fissure resistance by acidization, and (7) acidization failure. This theory has been successfully applied for over three years during which 22 acid treatments have been carried out specifically to reduce skin resistance and three treatments to reduce coarse fissure resistance. The method has been found to be a useful aid in the selection of wells for acidization.

By David Cornell, Yusuf K. Sukkar

The fundamental differential equation for fluid flow has been rearranged, integrated numerically for natural gases, and the results presented in tabular and graphical form suitable for the direct calculation of vertical gas flow problems. The integration is rigorous except for the use of a constant average temperature. Both static and flowing bottom-hole pressures may be calculated by (his method without resort to trial and error procedures. INTRODUCTION Bottom-hole pressures can be determined either by direct measurement with a bottom-hole pressure gauge or by calculation. Since measurement with a bottom-hole pressure gauge is costly, calculation of bottom-hole pressures is to be preferred when possible. Calculation involves knowledge of the wellhead pressure, the properties of the gas, the depth of the well, the flow rate, the temperature of the gas, and the size of the flow line. These quantities must be combined in a suitable equation for vertical flow. Several equations have been developed for the purpose of calculating gas flow through both horizontal and vertical pipes1,2,4,5,6,7,8,9 These equations all have the basic differential equation for flow as their starting point. They differ only in the assumptions which are made in order to simplify the integration step. The bottom-hole pressures calculated using any of the better known equations are in good agreement with bottom-hole pressures calculated by a more rigorous but considerably more tedious stepwise integration method. All of these methods, however. require a trial and error solution for the calculation of flowing bottom-hole pressures. The purpose of this paper is to present a method which provides a direct solution of static and flowing bottom-hole pressures without use of a trial and error solution. At the same time the method involves fewer simplifying assumptions than any of the previous methods used, except lor the lengthy stepwise method. The derivation of any practical flow equation begins with the basic differential equation which is rigorous for almost all flow problems. The five terms in Equation 1 represent the potential energy, the potential energy of position, the kinetic energy, the shaft work done, and the energy losses due to motion of the fluid respectively. Equation 1 reduces to Bernoulli's theorem for the case of horizontal flow of an incompressible fluid, no work done, and no energy losses due to friction. In the present case, the term (v) (dv)/(g,.) can be shown to be quite negligible compared to the other terms, and is therefore neglected in the derivation. The term — dws includes work done by a pump or turbine in the system and is zero in this case. The frictional losses in the pipe are included in the term — d(lw) which is expressed in terms of the Moody friction factor as in Equation 2. (f) (v2) (dL) d(lW) = (f)(V2) (dl)/(2) (g) (d) The specific volume of a gas can be expressed in terms of the temperature, pressure, molecular weight, and compressibility factor. The velocity of the gas can be expressed in terms of the weight rate of flow, the specific volume, and the dimensions of the pipe.The reduced pressure Pr is to be used because z can be expressed as a function of Pr and Tr alone for most natural gases:

Jan 1, 1956

By M. W. Legatski, D. L. Katz

The best currently available description of the longitudinal mixing properties of a porous medium is an equation of the form which relates the effective longitudinal dispersion coefficient Dp to the molecular diffusion coefficient Do, the electrical resistivity factor F, the porosity <f> and a Peclet number. If the parameters dpa and m are determined for a porous medium of known porosity and electrical resistivity factor, then a dispersion coefficient may be estimated for a given flow rate and a given gas pair. A new method, featuring on-line gas analysis by thermal conductivity and on-line data reduction by analog computation, was developed and used to determine these mixing parameters for eight naturally occurring sandstones and two dolomite samples. The exponent m of the above equation was found to vary between 1.0 and 1.5. The characteristic length dpa in the above equation was found to vary between 0.25 and 1.9 cm, with an average value of 0.4 cm for sandstones. Measurements were made on two cores in which paraffin wax had been deposited by evaporation from a pentane solution. They indicated that the presence of an immobile phase such as connate water could increase the dispersion coefficients significantly. INTRODUCTION While the petroleum and chemical industries have studied the mixing of miscible liquids flowing in consolidated porous media and of miscible gases flowing in unconsolidated porous media, relatively little data have been presented to describe the mixing of gases flowing through consolidated porous media. Such data are of particular interest to the gas storage industry. For instance, the U. S. Bureau of Mines is storing large quantities of a rich helium-nitrogen gas in contact with a natural gas in a dolomite reservoir. Since the rich gas occupies only 15 percent of the total reservoir volume, it is essential that the extent of rich gas-natural gas mixing be predicted and understood as a function of rock properties, pressure and rate of movement. This investigation was concerned only with the determination of longitudinal dispersion coefficients. It is understood that a transverse dispersion coefficient, which characterizes mixing perpendicular to the direction of flow, may be an order of magnitude less than the coefficient characterizing mixing in the direction of bulk flow.5719 It should also be recognized that the use of any dispersion coefficient is in itself a simplification. It is necessary to assume that mixing in a porous medium may be characterized by the equation for flow in a single direction. A number of authors1 l8,19 have pursued the mixing problem, not in terms of the so-called "dispersion model" described by Eq. 1, but in terms of a "mixing cell model". This model supposes that a porous medium is constructed of a large number of small mixing chambers and that the concentration of the diffusing component within each mixing chamber is uniform. Fick&apos;s law (Eq. 1) assumes that there is no gross by-passing of one fluid by another, and that there are not stagnant pockets of gas in the system under consideration as discussed by Coats and Smith.8 These assumptions are not always valid for flow through porous media and it is important to recognize the limitations upon Eq. 1. The most serious limitation upon the assumption of Fick&apos;s law mixing in a porous medium is

By D. E. Marks, Arnold, C. W, R. J. Robinson, A. E. Hoffmann

The efficiency of operation of a fixed-bed adsorption unit is infEuenced both by the absolute adsorption capacity of the bed and by the rate of adsorption. This paper describer studies of adsorption rate which were conducted in an experimental unit designed such that conditions existing in the treatment of high-pressure natural-gas mixtures could be duplicated. Variables investigated included pressure, temperature, gas composition, adsorbent particle size, depth of packed bed and gas velocity. The adequacy of a simplified mathematical model for predicting the observed phenomena was tested. A correlation is preserited which relates adsorption rate to the process variables stlldied. This correlation is useful in combination with the matheinatical model. INTRODUCTION Of the techniques available for contacting adsorbent particles with fluid streams to be treated, fixed-bed adsorption columns offer definite advantages in simplicity and ease of operation. As a result, they are often used in preference to others for such petroleum industry applications as dehydration and purification of natural gas and hydrocarbon recovery. Fixed-bed adsorption units usually consist of two or more towers filled with a desired adsorbent and operated in a cyclic manner. While one is being used to process the main flow stream, the others are undergoing regeneration to remove the adsorbed phase. When the tower on stream becomes saturated with the preferentially adsorbed material, the roles of the towers are switched, and the freshly regenerated tower is placed on stream. Cacle duration is determined by the bed capacity under the process conditions and by the flow rate through the bed. The sharpness of separation which can be effected is a function of both the absolute capacity of the bed and the rate of adsorption in the bed. The effect of rate for a particular set of conditions is evidenced by the sharpness or diffuse-ness of the adsorption front as it advances through the bed. Since data needed for design of adsorption units to treat high-pressure natural-gas systems were not available, an experimental program was designed to investigate the effects of different variables upon adsorption rate in fixed beds. In the present paper, effects of gas composition, column length, temperature, pressure, adsorbent particle size and flow rate (actual linear flow rate of the gas) are shown, and utility of a simplified mathematical model for describing the process is discussed. As gas enters the top of a cool, clean bed of adsorbent, preferentially adsorbed materials are stripped from the main flow stream by the uppermost particle layers. As these layers become saturated with a particular component, new supplies of this component are carried further down the column until fresh adsorbent is encountered. An adsorption wave thus moves through the column as material is supplied to saturate succeeding elements of the bed. Adsorption from a Multicomponent gas stream occurs as a succession of such moving waves corresponding to the different components in the gas. The leading edge of an adsorption wave for a component of a natural-gas stream moving through a bed of a common commercial adsorbent such as silica gel would be sharp but for the influence of certain broadening fac tors. These factors include a nonuniform velocity profile in the bed, longitudinal dispersion or mixing in the main gas stream, and the time required for a molecule to migrate from the main gas stream and be adsorbed at a site within the body of an adsorbent particle. If packing is uniform and the ratio of column to particle diameter is greater than approximately 15:1, the first factor is relatively unimportant&apos; Longitudinal mixing is of importance only for the case of moderately high mass transfer with extremely slow flow rates.&apos; The sharpness of an adsorption front, therefore, is, primarily a function of the rate of adsorption or the time required to saturate a particle of zdsorbent. Two methods for defining adsorption rate are used in this work. The first is a normalized or relative rate which describes the rate of saturation of a differential element of the packed bed. This can be measured by observing the time required for the concentration of the preferentially adsorbed material in the effluent gas from the bed to rise from zero to a value equal to that in the inlet gas stream. The second definition describes the absolute rate of mass transfer from the gaseous to the adsorbed phase. This definition is used in a mathematical description of the adsorption process. If the concentration of a component in the gas strcam leaving an adsorption column is measured and plotted as a function of time, a curve such as that shown in Fig. I results. It is seen that for a period of time the effluent gas is devoid of the component under consideration. As the bed approaches saturation, a small percentage of this material will appear in the effluent gas. The concentration will then rise with time, or increasing cumulative gas flow, until it is equal to that in the inlet gas stream. If adsorp-

By K. L. Young

The general energy equation, including change in kinetic energy, was solved by numerical integration and used to evaluate simplifying assumptions and application practices over a wide range of conditions. When extreme conditions were encountered, sizable errors were caused by large integration intervals, application of Simpson&apos;s rule and neglecting change in kinetic energy. A maximum error of only 1.31 percent was caused by assuming temperature and compressibility constants at their average value. It was discovered that a discontinuity can develop in the integral for the injection case. This discontinuity indicates a point of zero pressure change and is an inflection point in the pressure traverse. INTRODUCTION When a pressure in a gas well is to be calculated, one of the first decisions is to select a method of calculation. In many instances, this selection becomes a problem because the literature, at best, provides an evaluation of any method for only a limited range of conditions. Once a method has been selected, a question often arises as to the size of the calculation interval which should be used. The question regarding calculation interval arises because an analytic solution is not obtainable and approximate solutions must be used. This paper presents an evaluation of major assumptions and application practices of probably the two most widely used methods for calculating steady-state single-phase gas well pressures. The two methods are Cullender and Smith1-" (numerical integration), and average temperature and com-pressibility. The Cullender-Smith method assumes that change in kinetic energy is negligible and is normally applied in two steps with a Simpson&apos;s rule correcton. The average temperature and compressibility method, in addition to neglecting kinetic energy change, asumes that temperature and compressibility are constant at their average values. This method is normally applied for wellhead shut-in pressures of less than 2,000 psi, and in one step. Computer programs were written to compute bottom-hole pressure with and without the assumptions, using various approaches. Values of input parameters investigated are shown in Table 1. Flow rate was limited to a maximum of 5,000 and 10,000 Mcf/D for tubing sizes of 1.610 and 1.995 in. ID, respectively. Flow rate was also limited to 10,000 Mcf/D for a tubing size of 2.441 in. ID when wellhead flowing pressure was 100 psia. These limitations were imposed on flow rate so as not to exceed sonic velocity. The z factor routine available necessitated limiting bottom-hole temperature to 240F and wellhead pressure to 3,000 psia. Pressures were compared on the basis of percent deviation from the trapezoidal integration of Eq. l or 2 at 100-ft intervals. A preliminary investigation indicated that a 1,000-ft interval solution would differ from a 50-ft interval solution by less than 0.25 percent; therefore, the 100-ft interval was chosen for a base. For the purpose of comparison, deviations less than 1 percent were considered insignificant. where 111.1 qeC2/d p = kinetic energy term. Eqs. 1 and 2 can be evaluated numerically at specific depths using the trapezoidal rule as shown by Cullender and Smith. If change in kinetic energy is neglected and temperature and compressibility are assumed constant at their average values, Eq. 1 can be integrated to give the average temperature and compressibility equation,

By D. R. Parrish, T. M. Geffen, R. A. Morse, G. W. Haynes

Flow tests on small core plugs have indicated that a large amount of gas is trapped and not recovered by water flooding a gas sand. Instead of I to 15 per cent pore space, as is usually assumed, the residual gas saturation is 15 to 50 per cent pore space, and is thus of the same magnitude as residual oil after water flooding oil sands. A thorough investigation was made to ascertain that large amounts of residual gas actually remain in reservoirs after a water flood and that this condition is not merely a laboratory phenomenon. In field experiments, the amount of gas left in a watered-out gas sand was measured by use of a pressure core barrel and the residual gas saturation of two watered-out gas sands was determined by electric log evaluation. In the laboratory, an investigation was made of factors that could possibly cause the value of residual gas saturation as measured on small core plugs to differ from that in the reservoir, and the effect of these factors on the amount of residual gas saturation was studied. The factors studied include flooding rate, static pressure, temperature, sample size and saturation conditions before flooding All evidence established that a relatively high gas saturation is trapped in water flooded gas sands and that this residual gas saturation can be measured in the laboratory by tests on small core plugs. INTRODUCTION There. has been general agreement among engineers that very high recovery of gas could be obtained from natural reservoirs By water displacement. Gas recoveries of 80 to 95 per cent of the original gas in place have become the normal expectation in water drive fields. The assumption of high recovery has been based on: 1. low density and viscosity of gas compared with water; 2. the erroneous assumption that the flow relationship in a gas-liquid system where gas is the displaced phase will be the same as when it is the displacing phase. It has long been recognized that gas can flow at very low gas saturations (in the range of 1 to 15 per cent pore space) in systems where liquid is being displaced by gas. By assuming the reversibility of this process, the conclusion was reached that the residual gas saturation following water flooding of a gas reservoir would be the qame (1 to 15 per cent) as that at which gas first flowed continuously as a displacing phase. Recent laboratory relative permeability Studies have demonstrated that the flow characteristics are very different in gas-liquid systems, depending on whether gas is displacing or being displaced by, a liquid. Also, it has been shown that there is no difference between the flow characteristics of oil and water or gas and water in water wet porous rocks. The residual gas saturation that can be expected following water flooding of a gas reservoir then would be in the same range as the residual oil saturation normally expected after water flooding an oil reservoir! i.e., in the range of 15 to 50 per cent pore space, depending on the rock characteristics. Obviously, such a difference in residual gas saturation means very important differences in recoverable gas reserves from water drive reservoirs. For example, if the original gas saturation in a field were 70 per cent, and the residual gas following flooding were 35 per cent, only half of the gas in place could be recovered by complete water drive, compared to the previously expected 80 to 95 per cent. This is a situation in which complete pressure maintenance could result in very greatly reduced recovery, since straight pressure depletion recoveries from gas reservoirs can approach 80 to 90 per cent. Such a change in thinking must be based on more complete information than a series of small core tests. Hence the work reported herein was undertaken with the objective of determining whether or not residual gas saturations indicated from small core relative permeability tests at atmospheric pressure and room temperature are representative of the residual gas saturations which could be expected after water flood of natural reservoirs. A study was made both through laboratory and fields tests to determine any differences in residual saturation which might be occasioned by differences in pressure, temperature, rate of flooding, and original saturation conditions. Four separate types of laboratory experiments and two field mesurements were made in this investigation. The apparatus and testing methods of each will he discussed individually. LABORATORY EXPERIMENTS, APPARATUS AND PROCEDURE Relative Permeability Tests Steady state flow experiments Here conducted using the Penn State type apparatus. The equipment used has beer (described in a previous publication? "Irreducible" water satu-ration was established in the core&apos; by imposing ,a capillary pressure of 45 psi before simultaneous flow of air and water

Jan 1, 1952

By J. S. Crump, C. R. Hocott, A. E. Hoffman

Planning of the efficient operation of a gas-condensate reservoir requires a knowledge not only of the gross phase behavior of the system but also of the equilibrium distribution of the various components between the gas and condensate phases. This equilibrium distribution can be calculated with appropriate equilibrium constants. In this paper are presented equilibrium constants determined experimentally for the oil and gas phases initially present in the same reservoir and for the gas and condensate phases of the gas cap material at a series of pressures below the original reservoir pressure. Also presented is a method for the correlation of the experimentally determined equilibrium constants. The utility of the correlation is demonstrated further by an example of its use in adjusting the equilibrium data to permit their application to another gas-condensate system of similar composition. INTRODUCTION Planning of the efficient operation of a gas-condensate reservoir requires a thorouugh knowledge not only of the gross phase behavior of the particular hydrocarbon system but also of the equilibrium distribution of the various components between the gas and condensate phases. At the initial conditions of reservoir temperature and pressure, the original hydrocarbon materials in a gas-condensate reservoir or in the gas-cap of an associated reservoir exist in a single, homogeneous vapor phase. However, some condensation of hydrocarbons to a liquid phase usually occurs in the reservoir as pressure declines incident to production. Because of this condensation. the produced gas changes composition continuously. The composition of the produced gas, as well as that of the condensed liquid, can be calculated from the composition of the original reservoir material through the application of appropriate equilibrium constants or "K"-values provided these values are known. Equilibrium constants have been used for many years in problems of surface separation of gas and oil and natural gasoline recovery. However, no satisfactory equilibrium data at reservoir conditions have been available for the higher boiling hydrocarbons which acquire abnormal volatility at high pressure and are therefore present in the gas phase in condensate reservoirs. Early work in this field indicated that the behavior of the higher boiling hydrocarbons largely determines the behavior of gas-condensate systems at high pressure.&apos; Consequently, a portable test unit was designed to permit determination of the distribution of the heavier hydrocarbons between the gas and liquid phases of a gas-condensate system. This paper describes the test equipment and presents the results of an investigation of the reservoir fluids from an associated oil and gas reservoir in the Frio formation in Southwest Texas. Also presented is a method for correlation of the experimentally determined equilibrium constants, together with an example of the use of the method in adjusting these data for application to another gas-condensate system of similar composition. FIELD TESTING AND SAMPLING In order to determine the compositions of the reservoir fluids, field samples were secured from two different wells, one completed in the oil zone, and one completed in the gas cap. The oil well was completed through perforations in 51/2-in. casing at a depth of 6,822-6,828 ft subsea. Pressure and temperature traverses made in the well indicated the pressure and temperature to be 3,831 psig and 201 °F, respectively, at a depth of 6.828 ft subsea. Subsurface samples were taken from this well after a four-hour shut-in period.

Jan 1, 1953

By F. R. Scauzillo

Equilibrium data which should be useful in the design and/or evaluation of glycol dehydration units were prepared from an analysis of various published data and the correlation of these data by the use of the thermo-dynamic equilibrium ratio. The equilibrium ratios of water are used to solve the glycol absorber problem; such solutions are necessary to define the number of trays and the glycol circulation rate needed to meet drying requirements. Activity coefficients were obtained which relate directly to the equilibrium ratios of water in the water-TEG-natural gas system. These activity coefficients have been used to calculate the equilibrium dew points for aqueous TEG concentrations of 60 to 99.9 weight per cent for the temperature range of 40° to 120° F. They also provide a means of calculating equilibrium ratios for water in the water-TEG-natural gas system; this applies to any desired TEG concentration and to the temperature range from 40° to 120°F. INTRODUCTION During the last several years, the drying of natural gas with aqueous triethylene-glycol (TEG) solutions has become very prominent. Most users of TEG as a drying agent have been satisfied with the performance of TEG solutions at the conditions used; however, there always has been some discussion of the drying ability of TEG solutions at conditions not commonly encountered such as temperatures below 50° or 60°F, but more particularly temperatures above 100°F, and pressures above 1,000 psia. Today, when higher wellhead temperatures such as 120°F are more commonly encountered many are skeptical of TEG&apos;s ability to dry sufficiently well to provide dew points around 32°F, which is generally the maximum tolerable when attempting to dry a gas to contract specifications of 7 Ib/MMscf. These views probably have evolved to some extent from the days when suppliers would not guarantee dew-point depressions in excess of 65° to 75 °OF. Also, the feelings about TEG drying may have arisen from the lack of information about the equilibrium relation of water in the drying operations. Porter and Reid&apos; have reported equilibrium data for a 95 per cent TEG solution, while Townsend2 reported equilibrium data for 95, 98 and 100 per cent TEG solutions. Townsend also presented a calculational method which obtains activity coefficients in the liquid phase and, subsequently, the equilibrium water content of a saturated gas over glycol from published3 atmospheric dew points. Wise, Puck and Failey&apos; presented activity data for the water-TEG system at atmospheric pressure. In the past, the equilibrium ratios published for water in a 95 per cent by weight aqueous TEG solution have been used indiscriminately by many for all concentrations of TEG encountered in gas drying operations. Since essentially all such gas operations require the use of more concentrated TEG solutions, a study was undertaken to correlate the existing equilibrium data for aqueous TEG and gas systems and to provide some means of calculating the equilibrium ratio of water in the natural gas-water-TEG system where TEG concentrations were other than 95 per cent. The data presented herein consider only the equilibrium drying ability of the TEG solutions and do not consider the effect of temperature and pressure on the tray eficiency of contactors. In other words, this paper is concerned primarily with the development of the equilibrium relationship between water in the dried natural gas and the water in the lean TEG entering the top tray of the absorber. BASIC EQUILIBRIUM RELATIONS A relationship which relates the K value of water in the natural gas-water-TEG system to the vapor pressure of water at the system temperature, the total pressure of the system and the activities of the liquid and gas phases has been evaluated. It is well known that, for any component of a mixture at equilibrium between a gas phase and a liquid phase,