Search Documents

By L. S. Reid, Irwin Politziner F. M. Townsend

Recently published data indicate that the water vapor content of a gas, as determined by dew point measurement, is inaccurate when the gas has been dehydrated with diethylene glycol. Water vapor contents and dew point measurements of gases dehydrated with triethylene glycol have been obtained in this investigation in order to determine the magnitude of a similar error, if one exists. Experimental data show a very low concentration of triethylene glycol vapor in gases dehydrated at atmospheric ternperatures and pressures ranging from 500 to 2,500 psia and that the accuracy of dew point measurements is not impaired by the presence of triethylene glycol vapor. INTRODUCTION Nearly all natural gas transported by pipe line to northern and eastern markets must be dehydrated to a low residual moisture content in order to assure maximum transmission efficiency and continuous delivery under peak load conditions. Without dehydration, water vapor often condenses in pipe lines in quantities sufficient to restrict the flow3 and. under low atmospheric temperature conditions, gas hydrates may form and plug gathering. transmission and distribution facilities. Gas may be dehydrated to pipe line specifications by a number of methods. The most widely used methods are (1) adsorption of water vapor on a solid dessicant such as activated bauxite, alumina or silica gel,' (2) by absorption of water vapor by either diethylene or triethylene glycol-water solutions of high glycol concentration,11,13 and (3) by simultaneous expansion-refrigeration of very high pressure gas in which gas hydrates are purposely formed and then quickly decomposed.'"" Of these three fundamental methods, the adsorption and absorption processes provide the bulk of the dehydrated gas moving to market at the present time. The expansion-refrigeration process. a relatively new development, is gaining in favor because of its simplicity but is sharply limited in its application by available differentials between source and gathering system pressures. Where gas sale or purchase contracts contain a maximum water vapor content specification, the water vapor content is usually determined by measuring the dew point of the gas at system pressure. Knowing he dew point temperature and the system pressure, the water vapor content of the gas is determined from a graphical correlation of saturation, or dew point, temperature us the water vapor content of the gas which is customarily expressed as pounds of water vapor per million cubic feet of gas measured at 14.7 psia and 60°F. Dew points are usually measured with the U. S. Bureau of Mines Dew Point Tester. shown in Fig. 1, which is well suited for use at temperatures lower than 120°F and at pressurei; ranging from atmospheric to 3,000 psi. A number of correlations of dew point temperatures vs water vapor content are available,3,5,15 the most recent of which is that by McCarthy et al.4 Where natural gas is dehydrated by absorhing contact with concentrated diethylene glycol-water solutions. considerable difficulty is experienced in observing water dew points due to condensation of liquid hydrocarbon and glycol films on the mirrored surface of the dew point tester. Riesenfeld and Frazierl2 report that the dew point method for determining the water vapor content of a dietllylene glycol-treated natural gas is in error due to the glycol content of the condensate formed on the tester mirror, and that the water vapor content of the gas tested is actually lower than that indicated by the standard' dew point-water vapor content correlation. Extensive experience in testing gases dehydrated by concentrated triethylene glycol-water solutions lias not revesled difficulties of the nature cited above for diethylene glycol. Water dew points are clear and sharp and condensation of

Jan 1, 1951

By R. H. Olds, B. H. Sage

The transient flow of gases in long pipe lines is a problem of industrial interest. The present discussion deals with the application of the conservation of momentum and material to the transient flow of gas in conduits which are long with respect to their diameter. Solutions of the resulting equations include consideration of variations in friction and specific volume of the fluid as functions of time and position along the conduit. The changes in state along the pipe line are established as a function of time. The methods proposed may be extended to other situations. The flow of compressible fluids under steady conditions' has been considered in detail. Primary emphasis was placed on the development of simple means of taking into account the non-uniform flow' realized in long conduits under such conditions. Binder2 discussed the essential differences between the flow of compressible and uoncompressible fluids. Bonilla" proposed direct methods of solution for equations describing isothermal flow in long pipes. Similar considerations limited to steady flow have been discussed by Ruth.' Palsgrove 5 sugrested a graphical means of solving specific problems relating to the flow of gases in conduits of uniform section. Some of the basic concepts of the flow of a compressible viscous fluid in straight tubes have been described by Kuo.6 Problems of this type require knowledge of the equation of state of the fluid. In many cases it liar been convenient to assume that the fluid follows the relationship ascribed to a perfect gas. The treatment of unsteady flow of compressible fluids is less complete. This state of affairs result* from the fact that the differential equations relating pressure and flow rate with time and position for unsteady flow are nonlinear in character.' Muskat and Wyckoff presented an excellent treatment of the situation in laminar Bow where inertia forces were not of controlling importance. Somewhat earlier Hurst9 discussed the uncteady flow of fluids in petroleum reservoirs. Helhering-ton. MarRoberts. and Huntington'" established experimentally the pressure distribution and accumulative flow rate of natural gas through an unconsolidated sand. Recently Joffe11,12 conhitlcretl the specific problem of the flow of natural gas in relatively long.,. pipe lines for unsteady conditions with emphasis upon the economic feusibility of the storage of gas in such a reservoir. These discussions were based upon approximations of the partial differential equations applying to the condititons of flow. The results were descriptive of the physical situation and indirated the effertiveness of large pipe lines as cyclic storage facilities for natural gas. The advent of analog computers" has increased materially the practicality of solving nonlinear partial differential expressions. Several iterative numerical and graphical methods14,15,16 have been developed to solve nonlinear equations. However, it appears that the pair of simultaneous nonlinear partial differential equations involved in problems of this nature are not directly susceptible to solution by presently available analog computer or formal iterative numerical techniques. Nevertheless. graphical analyses afford a means of obtaining data of engineering value without special equipment. MATHEMATICAL CONSIDERATIONS Since the length of the conduit is large in comparison to the cross section in the cases to be considered, it is possible to treat the characteristics of the flow as independent of position in the cross section of the conduit. This possibility follows since the energy- connected with the velocity profile is small compared to that required to overcome the friction associated with its movement throughout the length of the line. Average values of the pressure. velocity. kinetic energy. momentum, and specific volume may he assigned to the fluid at any particular cross section. From a mathematical viewpoint this may be expressed in the following way for velocity: µ = µ(?X)..........(1) The uncertainties introduced by such assumptions usually are negligible when the ratio of length to diameter is greater than the order of 1.000. On this basis a volume element ma!. be considered to occupy the entire cross section of the conduit. Such an element is presented in Fig. 1.

Jan 1, 1951

By R. H. Olds, B. H. Sage

The transient flow of gases in long pipe lines is a problem of industrial interest. The present discussion deals with the application of the conservation of momentum and material to the transient flow of gas in conduits which are long with respect to their diameter. Solutions of the resulting equations include consideration of variations in friction and specific volume of the fluid as functions of time and position along the conduit. The changes in state along the pipe line are established as a function of time. The methods proposed may be extended to other situations. The flow of compressible fluids under steady conditions' has been considered in detail. Primary emphasis was placed on the development of simple means of taking into account the non-uniform flow' realized in long conduits under such conditions. Binder2 discussed the essential differences between the flow of compressible and uoncompressible fluids. Bonilla" proposed direct methods of solution for equations describing isothermal flow in long pipes. Similar considerations limited to steady flow have been discussed by Ruth.' Palsgrove 5 sugrested a graphical means of solving specific problems relating to the flow of gases in conduits of uniform section. Some of the basic concepts of the flow of a compressible viscous fluid in straight tubes have been described by Kuo.6 Problems of this type require knowledge of the equation of state of the fluid. In many cases it liar been convenient to assume that the fluid follows the relationship ascribed to a perfect gas. The treatment of unsteady flow of compressible fluids is less complete. This state of affairs result* from the fact that the differential equations relating pressure and flow rate with time and position for unsteady flow are nonlinear in character.' Muskat and Wyckoff presented an excellent treatment of the situation in laminar Bow where inertia forces were not of controlling importance. Somewhat earlier Hurst9 discussed the uncteady flow of fluids in petroleum reservoirs. Helhering-ton. MarRoberts. and Huntington'" established experimentally the pressure distribution and accumulative flow rate of natural gas through an unconsolidated sand. Recently Joffe11,12 conhitlcretl the specific problem of the flow of natural gas in relatively long.,. pipe lines for unsteady conditions with emphasis upon the economic feusibility of the storage of gas in such a reservoir. These discussions were based upon approximations of the partial differential equations applying to the condititons of flow. The results were descriptive of the physical situation and indirated the effertiveness of large pipe lines as cyclic storage facilities for natural gas. The advent of analog computers" has increased materially the practicality of solving nonlinear partial differential expressions. Several iterative numerical and graphical methods14,15,16 have been developed to solve nonlinear equations. However, it appears that the pair of simultaneous nonlinear partial differential equations involved in problems of this nature are not directly susceptible to solution by presently available analog computer or formal iterative numerical techniques. Nevertheless. graphical analyses afford a means of obtaining data of engineering value without special equipment. MATHEMATICAL CONSIDERATIONS Since the length of the conduit is large in comparison to the cross section in the cases to be considered, it is possible to treat the characteristics of the flow as independent of position in the cross section of the conduit. This possibility follows since the energy- connected with the velocity profile is small compared to that required to overcome the friction associated with its movement throughout the length of the line. Average values of the pressure. velocity. kinetic energy. momentum, and specific volume may he assigned to the fluid at any particular cross section. From a mathematical viewpoint this may be expressed in the following way for velocity: µ = µ(?X)..........(1) The uncertainties introduced by such assumptions usually are negligible when the ratio of length to diameter is greater than the order of 1.000. On this basis a volume element ma!. be considered to occupy the entire cross section of the conduit. Such an element is presented in Fig. 1.

Jan 1, 1951

By Edmund C. Bovey

"AS I begin my remarks this morning, please let me say how honoured and gratified I was to receive your invitation to address you here during your deliberations in historic Quebec City. The high esteem in which your Institute is held is reflected in the respect and deference accorded to it by industry across Canada. I feel privileged to be here with you today.I would like to discuss with you some applications of natural gas in the mining industry and suggest to you, in closing, some possible avenues of future development.Natural gas first became avail-able for industrial purposes in Northern Ontario in late 1958. Acceptance of natural gas, as a source of energy by the industry of that area, has been both widespread and gratifying. The mining industry of Northern Ontario has in-creased its annual use of natural gas by approximately 35 per cent each year since 1958. We believe that this trend will continue, al-though we will probably have to accept a more modest 10 per cent annual gain as the market approaches saturation. I would like to explain how we were able to do this and where all this gas is used. Our first sales were largely a matter of price-and the savings we offered were substantial."

Jan 1, 1966

By Francis J. Quinn

The Natural Gas Act, of 1938, which put interstate pipelines under the jurisdiction of the Federal Power Commission, did not place a protective competitive shield around such companies - nor does any provision of the act guarantee them a profit. Natural-gas pipelines are competitive with each other and also with the coal and oil industries. Gas distributing companies must compete not only with the coal and oil industries, but also with the electric industry. Natural gas has captured markets from oil, which had taken them from coal. Oil men contend that gas prices have: been, held to an artifically low level by regulation. The electric industry is moving into the space-heating market, which pleases the coal people because the coal industry's largest customer since losing the railroads to dieselization has been the electric companies. The coal industry is attempting to monopolize the industrial and electric generating market through government intervention, and this is being done under the guise of concern for the national welfare through advocation of a "National Fuels Policy purporting to, further the cause of fuel conservation. Such a policy is contrary to American principles and would result in the individual's loss of freedom of choice and ultimately, in less reserves of preferred fuels because of diminished exploration efforts.

Jan 1, 1961

By Louis V. Bender

THIS paper gives a short review of the installation for and the use of gas, as a. fuel, at the Anaconda Reduction Works of the Anaconda Copper Mining Co. Before putting in gas tile fuels used were pulverized coal, lump coal, coke distillate and oil. These fuels are now almost entirely eliminated. For nearly fifty years bituminous coals from Wyoming and Montana fields had been used. These coals were of a good grade, analyzing 11,300 B.t.u. and having an ash content of 12 per cent. When hydro-electric power replaced steam-driven equipment, the coal requirements were reduced; and, on account of improvements in metallurgy and/ or in furnace design, the quantity was lessened. Later, oil and distillate were introduced and used in some of the operations. All of these things tended toward the lessening of the amount of coal used, and had it not been for the increase in the size of the plant, the tonnage used would have been small when the change was made to gas.

Jan 1, 1932

By E. R. MARBLE

THE introduction of a new fuel, such as natural gas, necessitates careful study where it has not been used previously. At the El Paso smelter natural gas required the installation of apparatus with which no one at the plant was familiar, either as to the type of equipment needed or as to the burning of the gas in it. Observation of existing installations of a similar character and assistance from gas engineers served as a basis to determine what was advisable to use. The gas line entering the plant goes to a metering station equipped with a duplicate system of pressure regulators, a temperature recorder and the gas company and smelter orifice meters in tandem. There are three distribution lines that can be inter-connected at a later date, if necessary, to assist in equalizing the plant load. All the larger lines are of heavy

Jan 1, 1930

By Lyon F. Terry

IN contemplating the prospects of natural gas being transported from the fields where it is produced to such distant points as Philadelphia, New Jersey, New York City, and New England, let us review the advance that natural gas has made and the projects currently scheduled for construction. This may give us a perspective of the point in the long-range development of the natural gas business in which we now find ourselves, and some basis for judging the further expansion of the industry. Authentic statistics of the total marketed production of natural gas in the United States have been compiled by the Bureau of Mines since the year 1906. The growth of the business is indicated by the following data on marketed production for each ten-year interval since 1906:

Jan 1, 1947

By R. T. Elworthy

Canada ranks second in the list of world consumers of natural gas. In 1923 our production was nearly fifteen million thousand cubic feet valued at $5,800,000. Our output, however, appears very small compared with that of the United States, which in 1923 totalled over one thousand and eight million cubic feet, valued at $240,000,000. Natural gas is found in commercial quantities in three? provinces in Canada. New Brunswick: Natural gas in New Brunswick has been utilized commercially only since 1912,although the investigation of oil and gas possibilities of the province commenced as early as 1859. The annual production is a little more than one-twentieth of the total Canadian production and remains fairly constant from year to year. Ontario: In Ontario the earliest gas wells were drilled in Essex and Welland counties between the years 1885 and 1890. On account of the small demand in Canada the cities of Detroit and Buffalo were supplied with large quantities of gas from these wells, but this export was curtailed in 1898 and finally ceased in 1900. The production of natural gas in Ontario, mainly from fields in Lambton, Kent, Essex, Norfolk, Haldimand and Welland counties, increased rapidly from 1900 until the peak was reached in 1917, when over nineteen million thousand cubic feet valued at $3,641,587 was consumed.

Jan 1, 1925

By E. D. Glass, R. C. Hessing

The volume of gas in storage reservoirs may be compuled from estimates of hydrocarbon pore volume and gas density. However, both are dificult 10 estimate accurrately. Further, no adequate method has been presented for estimating reservoir performance during operation for Volumetric gas-storage reservoirs. Normally, gas-storage reservoirs exhibit small pressure gradients in the area conr(zilliitg the storage wells, even under conditions of maximum injection or withdrawal. A zone of low permeability usally surrounds, and is in communication with, tile permeable zone which causes the storage reservoir to exhibit a definite pressure lag, as evidenced by lack of pressure stabilization after extended periods of shut-in. An analysis is presented for determining pore volume and gas in place, and for predicting future reservoir behavior of such reservoirs. The mathematical development is presented and a solution is shown for an operating storage reservoir. INTRODUCTION Depleted gas reservoirs arc often used for gas storage. The reservoir usually contains a permeable area in which numerous storage wells arc located. The permeable zone permits high deliverability and injectivity with relatively small pressure loss. However, pressure equilibrium may not be attained even after long periods of shut-in because of the fact that the permeable zone is in communication with an adjacent tight section. Most of these gas-storage reservoirs are sand lenses, such as ancient beaches or offshore bars. The areal geometry of these reservoirs is normally such that flow in the tight zone will be approximately linear, and the solutions presented here should apply. DISCUSSION For purposes of this paper, a gas-storage reservoir is considered to consist of a volume of permeable rock in communication with tight rock, all of which is surrounded by an impermeable barrier. Pressure gradients across the permeable zone are small due to low flow resistance and well density. Hence, the permeable zone is considered to be at a uniform pressure which is only a function of time as gas is injected or withdrawn. The tight zone is considered to be linear with uniform permeability and porosity, and the pressure in the tight zone is a function of both position and time. l'wo solutions for flow in this system are used to define the reservoir and determine its properties. The first solution is for an infinite linear tight zone in communication with the permeable zone. The expression for constant gas-withdrawal rate is Eqs. I and 2 ore superposed to account for varying rates of production or injection and are used to fit gas-storage field reservoir history. The solution for the infinite case, Eq. 1, is first applied to field data for early times following a shut-in period. Constants BDa and c are ohtained by a least-squares fit of pressure and rate data. After the pressure wave from gas injection or withdrawal reaches the external boundary of the tight zone, the infinite solution will no longer apply and pressure behavior should be more nearly represented by Eq. 2 for a gas reservoir in communication with a finite linear tight zone. The values of the constants found from curve-fitting Eq. 1 to the early time data may then be used in Eq. 2 that only the coefficient- L containing the unknown

By Frank Marotta, Louis B. Lesem, Frank Greytok, John J. McKetta

Although one of the primary variables in the calculation of the flowing bottom-hole pressure in gas wells from surface measurements is the temperature at any point and its distribution in the flow-string, only few experimental data are available in the literature and little attention has been given to analysis of the problem. Virtually all of the recently published methods of calculating flowing bottom-hole pressures depend on the assumption that either the temperature is constant at some average value or that the variation is linear with depth. The purposes of this work are to analyze the problem theoretically and to verify the analysis by comparison with experimental data so that practical problems in the analysis of the behavior of gas wells can be solved with greater accuracy, reliability, and ease. MATHEMATICAL OUTLINE OF PROBLEM Assuming that: 1. The mass velocity and chemical composition of the gas stream are constant and in normal gas well operations, the change in linear velocity in the entire flow-string is trivial, 2. The product of the density and heat capacity of the gas is constant, 3. No horizontal temperature gradient exists in the gas stream, 4. Net flow of heat by conduction within the formation and in the gas stream in the vertical direction is trivial in magnitude and can be neglected, 5. The regional vertical geothermal gradient is constant, and 6. The temperature of the gas entering the borehole is constant, and equal to that of the reservoir, a set of two simultaneous, linear partial differential equations with appropriate boundary conditions was derived to describe the temperature distribution in the gas stream and the surrounding formation. These equations were solved by operational techniques for the distribution of temperature in the gas stream. The resulting integrals were evaluated numerically on an IBM 604 at the Machine Accounting Div. of the Railroad Commission of Texas. DERIVATION OF FUNCTIONS FOR THE DISTRIBUTION OF TEMPERATURE The physical system is a circular hole of radius a in an infinite medium with thermal conductivity K, density p8, and heat capacity c,. Under shut-in conditions the well is in thermal equilibrium with the surroundings and the increment in temperature is given by - ?T3/L x. Gas of density pg and heat capacity c.. flows upwardly through this hole in the direction of increasing x with a linear velocity W. Writing a heat balance on an element of gas ?x in thickness, convective heat transfer + conductive heat transfer + energy required to lift a unit mass of fluid = change in energy content. pa2W po Co To - pa2Wpo Co To x+xx - 2apK?x ?T8/?r r=a - Wpa2 po?x/778 = CoPopa2?x?To/?t... Dividing by (copopa2?x) and taking the limit as ?x tends to zero, -W (?To/?x + 1/778Co) + 2 K/a Po Co ?T3/?r r=u = ?To/?t

Jan 1, 1958

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 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. 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.