Search Documents

By R. V. Smith, E. J. Dewees, R. H. Williams

Before this paper was presented, the terms "f" and "v1/f in flow equations had been called "friction coefficient" and "friction factor," respectively. However, many authors have used the names for the terms interchangeably, with the result that considerable confusion has existed in the literature on flow of fluids. In an effort to establish a clear and concise terminology for these terms, the authors have designated f as "resistance coefficient" and dm as "transmission factor." Transmission factors for the flow strings of gas wells were determined by tests on wells in six fields of Oklahoma, Texas, and Louisiana. Data were collected on 14 gas wells while gas was flowing through 1 1/4, 2, 21/2, and 3-in. tubing and 51/2 and 7-in. casing. The transmission factors developed from these data were found to follow the resistance equation published by Nikuradse1 for turbulent flow in "rough" pipes with an average absolute roughness of 0.000,65 in. This average absolute roughness is believed to be representative generally of the flow strings of gas wells with the exception of those that have been exposed to corrosion or other conditions causing undue roughness. Also, use of an absolute roughness of 0.000,65 in. probably will result in transmission factors that are more representative of gas-well conditions than those in common use today, provided that the transmission factors are calculated by the methods outlined in this report. The authors believe that use of these transmission factors will permit calculation of more accurate subsurface pressures in gas wells than those in use today. Engineers use transmission factors in flow equations to compute subsurface pressures in flowing gas wells from measured wellhead pressures. These subsurface pressures are necessary for determining the productivity of gas wells from "back-pressure" tests. The methods of calculating pressures in flowing gas wells, as outlined by Rawlins and Schellhardt2 in Bureau of Mines Monograph 7, are based on the transmission factors in the well-known Weymouth formula" Natural gas engineers have realized that errors introduced by using transmission factors as given by the Weymouth formula are relatively unimportant in testing low-capacity gas wells; but they also know that such factors are important considerations in testing large-capacity wells. Accordingly, the research described in this report was undertaken to fulfill the need for more accurate means of calculating subsurface pressures in gas wells. This report is part of the Bureau's study of the flow of natural gas through pipe and commercial pipelines. Several of the conclusions regarding flow in the flow strings of gas wells were influenced by the results of the study of flow in horizontal pipelines. These instances are mentioned in the text of this report. METHODS OF COMPUTING TRANSMISSION FACTORS Methods and procedures for computing transmission factors for turbulent flow in circular pipes from test data taken on gas wells have been derived and explained in detail in a paper by Smith4. As these methods were used to calculate the transmission factors presented in this report, a summary is included. Starting with

Jan 1, 1955

By F. W. Jessen, R. W. Durie

The development in the laboratory of mall salt cavities illustrated that the rate at which salt dissolved was affected by the rate at which fresh water was injected into the cavity. Experiments carried out with small salt cores indicated, however, that the rate of salt removal with forced convection was not significantly different from that with natural convection, provided the flow was in the laminar range. The rate of salt removal was apparently controlled by the water salinity, which was indirectly determined by the rate at which fresh water was injected. With turbulent flow at the salt surface, the rate of salt removal could be increased several-fold with the same salinity in the cavity. A quantitative evaluation of the rate of salt removal under natural convection was attempted by using techniques developed for the analogous heat convection system. The values of the rate of salt removal obtained experimentally were considerably greater than the calculated values due to surface irregularities. To determine the effect of these irregularities, the rate of salt removal was found for a smooth surface by experimentally determining the rate of salt removal as an irregular surface developed and then extrapolating this trend back to initial time. Excellent agreement between the experimental and calculated values then were obtained. INTRODUCTION Underground cavities leached from massive salt formations have been used extensively for many years for storing LPG products. The development techniques for storing methane gases in the liquid state at a low temperature leads to speculation that underground salt cavities of controlled geometry may be quite important in this field. In any case, the increasing number of salt cavities used for storage indicates an impending need for a better knowledge of the washing process to efficiently develop maximum storage volumes within specified regions. It is the purpose of this paper to discuss some aspects of -the mechanism of the dissolution of salt in the formation of cavities of controlled geometry. Although no attempt is made to resolve the problems that may be associated with the storage of methane in the liquid state in salt formations, the fact that salt cavities can be formed at reasonable cost would appear to justify investigation in this direction. Further, underground cavities have the advantage of allowing storage at higher pressures (and hence higher temperatures) than surface installations. Although the thermal conductivity of salt is somewhat higher than frozen earths, 0.017 cal/(sec) (cm)(c) for salt at 32F vs 0.005 cal/(sec) (cm) () for ice at this temperature, this disadvantage is offset by the higher storage temperature in underground cavities. For example, at atmospheric pressure, the temperature of liquid methane would be —258F. In a storage cavity at 2,000 ft, the pressure would be approximately 1,000 psia, which would allow a storage temperature of -9OF. The formation temperature would be roughly 110F at this depth, hence the total temperature differential would be 200F. Under surface storage conditions the corresponding temperature difference would be some 333 F (75 F + 258 F). The results of this difference can be compared more readily by substituting the appropriate values of thermal conductivity, temperature difference, heat capacity and density into the solution of the diffusion equation for the unsteady linear conduction of heat. For the underground salt storage, the heat flux is 15cal/(sec) (sq cm), where t is the time in seconds from the initiation of the cooling. For the surface system, and using the thermal properties of ice at the average temperature, the corresponding figure is 21cal/(sec) (sq cm). Hence, the heat losses through ice under surface conditions would be some 30 per cent higher than the underground salt storage. This ratio would vary, considerablv. Of course, for other conditions, but it indicates that an advantage over surface facilities may be obtained by storage in underground salt cavities.

Jan 1, 1965

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By M. R. Tek, F. H. Poettmann, M. L. Grove

The general problem of single-phase natural gas flow through porous media has been considered both by mathematical analysis and field experiments. Theoretical studies on the fluid dynamics of natural gas wells resulted in derivations of expressions which permit a ber-ter understanding of the flow mechanism for conditions under which Darcy's law is no longer valid. In applying these results to natural gas production and testing technology, special emphasis war placed on the study of low permeability, slowly stabilizing formations. These applications resulted in a method of predicting the performance coefficients of tight natural gas wells. The calculation of the performance coefficient requires the determination of the "test-index", pkh/µ, the drain- age radius, R, and the effective wellbore radius, rw. Methods for estimating these parameters from field data were also developed. Field data from 15 wells were collected and analyzed. The measured values of the performance coefficients were compared with the predicted values. This statistical evaluation gave a bias of 7.14 per cent and a standard deviation of 11.6 per cent. INTRODUCTION The performance of a natural gas well is determined by the physical properties of the reservoir rock, the extent and the geometry of the drainage area, the properties of the flowing fluid, and the conditions of pressure distribution within the drainage matrix. The relation between the quantity of gas delivered and the pressure drop within the reservoir is characteristic of the behavior of each well and is often referred to as the back-pressure behavior of the well. It is well established that for normal, single-phase gas wells the relationship between the gas delivery rates and the bottom-hole pressure take, in general, the form, Q = c(P3 8 -P2 1)n. ,.......(1) where C is the performance coefficient, and n is the exponent corresponding to the slope of the straight-line relationship between Q and (P2 s, — P2 l) plotted on logarithmic coordinates. The determination of a well's performance by the back-pressure behavior not only permits the production engineers to analyze many operating problems, but also provides necessary information useful and essential to predict the future development of the field.' Important engineering and production problems such as calculation of gas deliverability into a pipeline at a predetermined line pressure, design and analysis of gas-gathering lines, determinations of spacing and number of wells to be drilled during field development to meet future market or contract obligations, etc., all depend, at least in part, on the availability and use of reliable back-pressure curves.' In the past the behavior of gas wells was evaluated

Jan 1, 1958

By R. E. Schilson, F. H. Poettmann

The direct determination of the stabilized performance behavior of low capacity, slowly stabilizing gas wells is extremely time-consuming and wasteful of gas. From both field experience and theoretical considerations, a test procedure has been evolved by which the stabilized hack pressure behavior of such gas wells can he predicted without having to revert to long time flow tests. The method consists of using the isochrona1 test procedure to establish the slope of the back pressure curve, "n", and the short time variation of the performance coeficient, "C", with time. From this short time transient flow data and theoretical considerations, the value of C at large times can he established. By assuming the radius of drainage of a well to be half the distance between wells, one can calculate the stabilization time for various well spacing patterns. Once the stabilization time for a given .spacing has been determined, the value of C can be calculated and the stabilized back-pressure curve can he. establbrhed. The calculated perfortnance coefficient as a function of tinze was cotnpared to the experintentally measured values for a number of gas wel1.e. The deviation of the calc~*lated from the cxperimental res1tlt.s vary depending on the set of short time experimental points used to evabrate the parcrttzeters of the equation. The longer the tirne far the flow test data user1 in the calculations, the better was the agreement with the experimental results. The time necessary to obtain this data from well tests varies consitlerably, depending on the physical natrcre oj the rt3scrvoir under consideratiot~. INTRODUCTION For many years, the U. S. Bureau of Mines Monograph 7' has served as a guide for testing and evaluating the performance of gas wells by means of the back-pressure method. The back-pressure performance of a gas well is expressed by the following equation: Q~ CW-W.........(I) where the characteristics of the back-pressure equation are determined by C, the performance coefficient, and n, the exponent which corresponds to the slope of the straight line when Q and {P* - PN2) are plotted on logarithmic paper. Q is gas flow rate at standard conditions, and P, and P, are equalized and flowing bottom-hole pressures, respectively. Prior to the development of the back-pressure tcst, the "open flow" capacity method of testing a well was common. By this method, the wells were flowed wide open to the air and the flow rate measured. Such procedure was wasteful of gas and did not provide information on the deliverability of the gas to the pipe line. Monograph 7 Procedure The back-pressure method of testing wells was developed to overcome these shortcomings. Although much has been learned regarding the laws of the flow of gas through porous formations, the original development of the back-pressure relationship was based entirely on empirical methods. The back-pressure behavior provides the engineer with information essential in predicting the future development of a field. It permits him to calculate the deliverability of gas into a pipe line at predetermined line pressures, to design and analyze gas gathering lines, to determine the spacing and number of wells to he drilled during the development of a field to meet gas purchasers' requirements, and to solve many other technical and economic problems. As described in Monograph 7, the flow-after-flow method of back-pressure testing, when applied to fast stabilizing and usually high capacity wells, correctly characterized the behavior of the well. However, as the value of the gas at the wellhead increased, small capacity gas wells having slow rates of stabilization became economically operable. The flow-after-flow method of testing could not be used to describe the behavior of these slowly stabilizing wells. The procedure of Rawlins and Shellhardtl for establishing the back-pressure behavior of a gas well was based on the rcquirement that the data be obtained under stabilized flow conditions; that is, that C is constant and does not vary with time. C depends on the physical properties of the reservoir, the location, extent and geometry of the drainage radius, and the properties of the flowing fluid. In a highly permeable formation, only a very short period of time is required for the well to reach a stabilized condition, and, consequently, the requirements for the test procedure described in Monograph 7 are met. For a given well, n is also constant

By H. O. McLeod, J. M. Campbell

This paper presents the results of the data obtained in the first stage of a long-range study at high pressures of the system, vapor-hydrate-water rich liquid-hydrocarbon rich liquid. The data presented are for the three-phase systems in which no hydrocarbon liquid exists. Tests were performed on 10 gases at pressures from 1,000 to 10,000 psia. One of these was substantially pure methane, and the remainder were binary mixtures of methane with ethane, propane, iso-butane and normal butane. Several conclusions may be drawn from the data. 1. Contrary to previous extrapolations, the hydrocarbon mixtures tested form straight lines in the range of 6,000 to 10,000 psia which are parallel to the curves for pure methane, when the log of pressure is plotted vs hydrate formation temperature. 2. The hydrate formation temperature may be predicted accurately at pressures from 6,000 to 10,000 psia by using a modified form of the Clapeyron equation. The total hydrate curve may be predicted by using the vapor-solid equilibrium constants of Carson and Katz' to 4,000 psia and joining the two segments with a smooth continuous curve between 4,000 and 6,000 psia. 3. The use of gas specific gravity as a parameter in hydrate correlations is unsatisfactory at elevated pressures. 4. The hydrate crystal lattice is pressure sensitive at elevated pressures. INTRODUCTION Prior to 1950 many studies had been made of the hydrate forming conditions for typical natural gases to pressures of 4,000 psia.""'"'"" Most of these attempted to correlate the log of system pressure vs hydrate formation temperature, with gas specific gravity as a parameter. One of the more promising correlations was made by Katz, et al, which utilized vapor-solid equilibrium constants. The only published data above 4,000 psia are those of Kobayashi and Katz7 for pure methane to a pressure of 11,240 psia. In the intervening years, most published charts for the high-pressure range have represented nothing more than extrapolations of the low-pressure data, with the methane line serving as a general guide. The reliability of these charts has become increasingly doubtful (and critical) in our present technology as we handle more high-pressure systems. The portion of our high-pressure hydrate research program reported here was designed to: (1) investigate the reliability of existing charts; (2) obtain actual data on gas mixtures to 10,000 psia; and (3.) develop a simple hydrate correlation that was more reliable than those which simply used specific gravity as a parameter. Binary mixtures of methane and ethane, propane normal butane, or iso-butane were injected into a high-pressure visual cell containing an excess of distilled water. Hydrates were formed and then melted to observe the decomposition temperature of the hydrates at pressures from 1,000 to 10,000 psia. EQUIPMENT The equipment consisted of a Jerguson 10,000-lb high-pressure visual cell, a 10,000-1b high-pressure blind cell and a Ruska 25,000-1b pressure mercury pump. The visual cell was placed in a constant-temperature water bath controlled by a refrigeration unit and an electric filament heater. A Beckman GC-2 gas chromatograph was used in analyzing the gas mixtures after each run was completed. EXPERIMENTAL PROCEDURE After evacuating the gas system, the heavier hydrocarbon was injected into the high-pressure mixing cell to that pressure necessary to give the desired composition. This cell then was pressured to 1,100 to 1,200 psia by methane from a high-pressure cylinder. The mixing cell holding the gas contained a steel flapper plate and was shaken intermittently over a period of 15 minutes. After mixing, the valve to the high-pressure visual cell containing excess distilled water was opened, and the gas mixture was allowed to flow into the cell. The temperature in the water bath was lowered 10" to 15'F below the estimated hydrate decomposition point. As a first check, the temperature was increased at a rate of 1°F every six minutes to find the approximate point of decomposition. It was again lowered 1.5° to 5°F to form hydrates. The temperature was raised to within l° of the estimated decomposition point and then increased 0.2F every 10 to 15 minutes until the hydrates decomposed. This procedure was repeated at various pressures to obtain 7 to 13 points for each mixture between 1,000 and 10,000 psia. After completion of the hydrate decomposition tests, the gas mixture composition was analyzed with a calibrated gas chromatograph. These gas analyses have an estimated error of ± .1 per cent.

By Anthony Folger, Howard R. Lowe, J. F. Dougherty, E. G. Trostel, Joffre Meyer

The status of the natural gas industry in Alberta, Canada, is described with particular reference to the current extent of natural gas reserves and possibilities for additional development. Certain major fields are discussed as representative of the types of gas accumulations within various broad geographical divisions of the provinces. The estimation of the future performance of a gas reservoir is discussed as the basis for an estimate of the future availability of natural gas from presently known reserves. The future availability of the provincial reserves and the estimated future demand for gas within the province are described. INTRODUCTION This paper is a condensation of a detailed study of the natural gas reserves and the future availability of pipeline gas in the province of Alberta as prepared for Trans-Canada Pipe Lines Limited and presented in public hearings at Edmonton and Calgary before the Petroleum and Natural Gas Conservation Board of the Province of Alberta. The study began in October, 1950, and resulted in the preparation of five volumes of data at the expenditure of more than 20,000 man-hours. This paper is designed to afford a broad picture with detail limited to a type example of both natural ga-accumulation and future availability of gas supply. The province of Alberta has an area of 255.285 square miles which, for sake of comparison. is 8,359 square miles smaller than the state of Texas. Natural gas was first discovered in 1883 at the now abandoned Alderson well located some 20 miles southeast of the town of Princess. The first gas in commercial quantities was found at Medicine Hat in 1890. During the last 68 years. 187 localities have been discovered which are capable of gas production, and 73 per cent of these have been found since 1945. Sixty fields currently are producing gas, including 22 non-associated gas accumulations. Location of the oil and gas fields and prospects in the province of Alberta is shown in Fig. 1, together with a list of the localities of measurable gas, Fig. 1A. The estimated total provincial gas reserve as of Jan. 1, 1952, for the 104 fields estimated, is 11.7 trillion cu ft proved and probable, and for the 187 known areas capable of production is in excess of 16 trillion cu ft on a total proved, probable, and possible basis. Five fields (Pincher Creek, Leduc-Woodbend, Viking-Kinsella, Cessford and Medicine Hat) are estimated to have recoverable reserves of more than one trillion cu ft each. A sixth field, Turner Valley, was initially in this group, but its reserves have been depleted substantially. Twenty-four fields, of which the above mentioned six fields are a part, have estimated gas reserves (including the possible category) in excess of 100 billion cu ft. These 24 fields account for 79 per cent of the total provincial reserves. Comparative systemic analyses reflect that 56 per cant of Alberta's gas reserves are contained in sediments of Cretaceous age; 30 per cent in sediments of Mississippian age; 10 per cent in sediments of Devonian age; three per cent in sediments of Triassic age; and the remaining one per cent in Jurassic and Permo-Pennsylvanian sediments. It is expected that the relative importance of the Mississippian, Devonian, and Triassic sediments will be increased as exploration progresses. GENERAL GEOLOGY The western Canada basin comprises the sedimentary area between the Canadian Shield and the Rocky Mountain Front Range, the Arctic Ocean, and the International Boundary. The basin has an approximate area of 765,000 square miles and embraces parts of the provinces of Manitoba. Saskatchewan. Alberta, British Columbia and the Northwest Territories. Generalized columnar sections for the Province of Alberta are shown in Fig. 2. Substantially all of Alberta, except its northeast corner, has sedimentary cover. Sediments of all the geologic systems are represented. except those of Silurian age.

Jan 1, 1952

By H. J. Ramey

The wellbore acts as a storage tank during drawdown and build-up testing and causes the sand-face flow rate to approach the constant surface flow rate as a function of time. This effect is compounded if non-Darcy flow (turbulent flow) exists near a gas wellbore. Non-Darcy flow can be interpreted as a flow-rate dependent skin effect. A method for determining the non-Darcy flow constant using this concept and the usual skin effect equation is described. Field tests of this method have identified several cases where non-Darcy flow was severe enough that gas wells in a fractured region appeared to be moderately damaged. The combination of wellbore storage and non-Darcy flow can result in erroneous estimates of formation flow capacity for short-time gas well tests. Fortunately, the presence of the wellbore storage eflect permits a new analysis which can provide a reasonable estimate of formation flow capacity and the non-Darcy flow constant from a single short-time test. The basis of the Gladfelter, Tracy and Wilsey correction for wellbore storage in pressure build-up was investigated. Results led to extension of the method to drawdown testing. If non-Darcy flow is not important, the method can be used to correct short-time gas well drawdown or build-up data. A method for estimation of the duration of wellbore storage effects was developed. INTRODUCTION In 1953, van Everdingen and Hurst generalized results published in their previous paper3 concerning wellbore storage effects to include a "skin effect", or a region of altered permeability adjacent to the wellbore. Later, Gladfelter. Tracy and Wilsey4 presented a method for correcting observed oilwell pressure build-up data for wellbore storage in the presence of a skin effect. The method depended upon measuring the change in the fluid storage in the wellbore by measuring the rise in liquid level. To the author's knowledge, application of the Gladfelter, Tracy and Wilsey storage correction to gas-well build-up has not been discussed in the literature. It is, however, a rather obvious application. Gas storage in the wellbore is a conlpressibility effect and can be estimated easily from the measured wellbore pressure as a function of time. Several approaches to the wellbore storage problem have been suggested. As summarized by Matthews, it is possible to minimize annulus storage volume by using a packer, and to obtain a near sand-face shut-in by use of down-hole tubing plug devices. Matthews and Perrine have suggested criteiia for determining the time when storage effects become negligible. In 1962, Swift and Kiel' presented a method for determination of the effect of non-Darcy flow (often called turbulent flow) upon gas-well behavior. This paper provided a theoretical basis for peculiar gas-well behavior described previously by Smith. Recently, Carter, Miller and Riley observed disagreement among flow capacity k,,h data determined from gas-well drawdown tests conducted at different flow rates for short periods of time (less than six hours flowing time). In the original preprint of their paper, Carter et al. proposed that the discrepancy in flow capacity was possibly a result of wellbore storage effects. Results of an analytical study of unloading of the wellbore and non-Darcy flow were recorded by carter.14 In the final text of their paper, Carter et al.!' stated that they no longer believed wellbore storage was the reason for discrepancy in their kgh estimates. In view of the preceding, this study was performed to establish the importance of non-Darcy flow and well-bore storage for gas-well testing. In the course of the study. a reinspection of the previous work by van Everdingen' and Hurst' was made, and the basis for the Gladfelter, Tracy and Wilsey' wellbore storage correction was investigated and extended to flow testing. WELLBORE STORAGE THEORY As has been shown by Aronofsky and Jenkins,11-12 Matthews," and others, flow of gas can often be approximated by an equivalent liquid flow system. The following developnlent will use liquid flow nomenclature to simplify the presentation. Application to gas-well cases will be illustrated later. First, we will use the van Everdingen-HursP treatment of wellbore storage in transient flow to establish (1) the duration of wellbore storage effects, and (2) a method to correct flow data for wellbore storage. DURATION OF WELLHORE STORAGE EFFECTS When an oil well is opened to flow. the bottom-hole pressure drops and causes a resulting drop in the liquid level in the annulus. If V. represents the annular volume in cu ft/ft of depth, and p represents the average density of the fluid in the wellbore, the volume of fluid at reservoir conditions produced from the annulus per unit bottom-hole pressure drop is approximately: res bbl-- (V, cu ft/ft) (144 sqin./sq ft) psi -(5.615 cu ft/bbl)(pIb/cuft) ........(I)

Jan 1, 1966

By A. Satter, J. M. Campbell

Reported herein are the results of a careful and detailed study of the non-ideal behavior of pure gases and their mixtures. Included are: (1) new data on five ternary systems composed of methane, ethane and H2 S; (2) a simple compressibility factor correlation that is inherently superior to present correlations, particularly for gases containing H2S and CO2; and(3) a detailed study of combination rules and the effect of system composition on the choice thereof. This study makes use of the rather large mass of data already available in the literature. A complete re-examination of the data and ideas presented in the last 25 years was considered desirable as a prelude to our basic concern — the effect of diluents on gas behavior. A consideration of both the macroscopic and microscopic properties of gases provides a better insight which, in turn, gives a firmer basis for improved correlation techniques. Such a study has shown that expressing the compressibility factor Z as a function of acentric factor w, as well as reduced temperature and pressure, yields a correlation that is broader in scope. The study of various combination rules has shown that better results are obtained by "tailoring" the rule used to the system composition. To do so improves the basic reality of results by overcoming some of the anomalies often found when using Kay's rule alone. Tentative recommendations are made regarding the most reliable combination rule for use with a given class of gas. The data presented are useful for estimating the direction and magnitude of the expected deviation when using a given rule. Although more work is needed, particularly around the critical region and with CO2 mixtures, the advantage of the classification scheme proposed is apparent. INTRODUCTION When one attempts to write a PVT equation to fit the data for actual gases, greater precision is obtained by the use of a multiple number of empirical constants. This has lead to multiple-constant equations such as Benedict-Webb-Rubin, Beattie-Bridgman, Keyes, etc., which are capable of yielding very precise results for pure gases in a range for which data to get the constants are available. As a matter of practicality, though, the use of such equations for gas mixtures is limited. Because of the infinite number of gas analyses available, any attempt to compile the constants needed requires a prohibitive amount of experimental data. This could be overcome by the use of a combination rule, but there is no real advantage in doing so because the end result offers no practical impovement over the Z factor correlation. The most widely used method of predicting the volumetric properties of pure gases is based upon the "theorem of corresponding states". According to this theorem, "all pure substances have corresponding molal volume at corresponding temperature and pressure if the reference point of correspondence is the critical point". Generalized compressibility charts for gases were prepared first by Cope and associates1 in 1931 and later by Brown and co-workers2 in 1932. However, the most commonly used charts are those of Dodge,3 Nelson and Obert,4 Hougen and atsson: and Standing and Katz.6 The work of Katz and co-workers has provided us with basic data for the hydrocarbons most widely used today. Their original chart6 was compared with a relatively large amount of multi-component data for gases consisting almost entirely of normal paraffin hydrocarbons. A deviation of only + 1.2 per cent was obtained.39 In the 20 years following publication of this work it has been found that the behavior of most mixtures of paraffin hydrocarbons could be predicted by this correlation within at least 5 per cent. Where difficulty has been encountered it has largely involved one or more of the following circumstances: pressures above 4,000 psig, mixtures containing large amounts of heavy ends and/or aromatics, systems in the critical region and mixtures containing polar compounds and/or CO2. The abnormal error sometimes found with such gases, not too unexpected for this method, is

By R. D. Carter

Methods are presented for calculating the performance of multiwell gas reservoirs. These methods account for two-dimensional, unsteady-state flow of a non-ideal gas through a heterogeneous formation. The performance can be obtained for gas reservoirs with known physical properties and configuration subjected to specified development and production schedules. These methods may prove helpful in optimizing the exploitation of gas reservoirs. The numerical solrrtion procedures given for this non-linear problem are extensions of known techniques for the linear diffusion equation. Two of these procedures are based on a new difference equation given by Saul'ev. Limitations and advantages of these procedures are discussed. An example solution is presented and is used to compare the calculation methods. INTRODUCTION Li planning the operation of a multiwell gas reservoir, it is desirable to predict its performance for each alternate development and producing schedule under consideration. This paper presents three numerical procedures that can be used to calculate such gas reservoir performance using a digital computer. These methods account for two-dimensional transient flow of a non-ideal gas within the reservoir. Extension of these methods to three-dimensional flow is possible. The paper is divided into three parts. The first describes the problem for which a solution is wanted and three numerical procedures that can be used to obtain solutions. The second discusses the information necessary to set up a problem. This part will be of greatest interest to the field engineer. The third gives data and results for an example problem. PROBLEM AND SOLUTION METHODS In the following development, a gas reservoir is regarded as a sealed, porous, permeable, heterogeneous body in which a single, non-ideal gas phase is flowing. Flow is permitted in the x- and y-horizontal directions, but is considered to be negligible in the vertical or z-direction. The reservoir can be tapped by one or more wells. Methods to be described provide a basis for calculating the complete performance of the reservoir, including rates and cumulative production from the wells, and pressure distribution and decline in the reservoir. If the isothermal flow of a non-ideal gas obeying Darcy's Law occurs in two dimensions in a region R, the governing equation in this region is:

Jan 1, 1967

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

An experimental method and apparatus for the study of the low-temperature phase and volumetric behavior of volatile mixtures are described. The phase diagrams, including the critical points, and the volumetric behavior of a natural gas and a mixture of the natural gas with methane at temperatures from -40°F to -200°F and pressures to 1.200 11si are givrn. INTRODUCTION Knowledge of the phase and volumetric behavior of mixtures of the light hydrocarbons at low temperatures is becoming increasingly important as interest in low-temperature separation processes increases. Distillation at low-lemperatures has been used for some time to recover ethane and ethylene from natural gas and refinery gas for petrochemical uses. Removal of excess nitrogen from natural gas by distillation has been proposed as a means of increasing the capacity of gas transmission lines and of improving the marketability of the gas.1,2 Evaluation of the latter process has been hampered, however, by the scarcity of data from which the phase behavior of natural gases at low temperatures might be predicted. While existing correlations of vaporization equilibrium constants may be used to predict the behavior of such mixtures at lower pressures, high-pressure computations often require knowledge of the location of the critical point of the system involved. Although correlations have been developed for the prediction of the critical temperature and pressure of certain types of hydrocar-bon systems, they do not apply to very volatile systems such as nitrogen-bear-ing natural gases. This paper describes the initial stages of an investigation of the low-temperature behavior of a numher of natural gases and other gases of commercial interest. The experimental method and apparatus are described. and the phase and volumetric behavior of a natural gas and a mixture of the natural gas with methane are presented. PREVIOUS WORK The published low-temperature vapor-liquid equilibrium data for hydrocarbons and mixtures of hydrocarbons with nitrogen and helium are listed in Table I. Most of these data are for binary systems, and the data of Stutz-man and Brown" at 100 psi are the only data on a natural gas. The data of Bloomer and Parent3 on the methafie-nitrogen system, which have been phblished since this study was completed, are a valuable addition to the knowledge in this field. Prior to the publication of these data the only published critical point data on systems having critical temperatures below the ice point were those of Ruhemann.1 on the methane-

Jan 1, 1952

By K. H. Kilgren

At high formation pressures the distillate produced from a gas-condensate reservoir may be black in color. In this event the dense gas phase existing above the dew point is correspondingly dark. Volumetric phase data and an analysis of a reservoir fluid exhibiting these characteristics, together with a description of the visual equilibrium cell in which these observations were made, are presented in this paper. INTRODUCTION Previously the author, like many others in the oil and gas industry perhaps, tacitly assumed that the expressions black or dark oil system, crude oil system and bubble-point system were synonymous. Crude oil reservoir fluids are bubble-point systems and yield a black or dark stock-tank oil of relatively low API gravity. Conversely, a clear or amber colored trap distillate of high API gravity is assumed indicative of a dew point system or a gas-condensate reservoir fluid. This broad classification appears satisfactory for shallow reservoirs, but as the following study demonstrates, may be misleading when applied to deep reservoirs. Theoretically, there is no reason to exclude the possibility of producing a dark, low-gravity distillate from a gas-condensate reservoir. At sufficiently high values of pressure and temperature, heavy, dark-colored hydrocarbons may exist in the vapor state of a multi-component system. If enough dark-colored components are present in the reservoir vapor phase, the resulting condensate will be dark. The reservoir fluid investigated in the present study supports this contention. Stock-tank production was black in color and measured 29" APT gravity. From outward appearances, the liquid closely resembled a medium gravity crude oil. Experimental measurements proved the reservoir fluid was in reality a gas-condensate system. Volumetric phase data for this high-pressure system and a description of the visual cell in which the study was conducted successfully are presented in this paper. THEORY Phase behavior of a reservoir fluid can be predicted accurately with reference to a pressure-temperature phase diagram. If the reservoir temperature is lower than the critical temperature of the hydrocarbon fluid in place, bubble-point behavior will be observed. If the reservoir temperature lies between the critical and cricondentherm temperature, dew point behavior and retrograde condensation will occur. For reservoir temperatures above the cri- condentherm, only a single gas phase can exist in the reservoir regardless of pressure. Providing the composition of the reservoir fluid were known, it would be possible to predict the critical temperature and estimate the phase behavior from equilibrium relationships. However, the usual practice is to obtain a sample of reservoir fluid, subject it to varying pressures at the reservoir temperature and observe the phase behavior experimentally. The latter method was used to obtain the data reported here. WELL AND TRAPPING INFORMATION A summary of pertinent data for the well from which the reservoir fluid was sampled is presented in Table 1. This well is located offshore Louisiana. Except for the pressure which substantially exceeds hydrostatic pressure, the information does not appear unusual. Prior to the sampling program, the well was produced for 22 hours at a stock tank oil rate of 139 B/D. Average trapping conditions and gauging data for the six-hour test period that followed are summarized in Table 2. Samples of the first-stage trap gas and liquid were obtained during the latter portion of the test period. Ambient temperature remained 5 to 10F below trap temperature and presented no problem for sampling. Surface wind and moderate foaming of stock-tank oil presented some difficulty in obtaining accurate stock-tank gauges. SAMPLE ANALYSIS Compositions of the gas and liquid samples are shown in Table 3. The trap gas was analyzed by isothermal chromatography which revealed only a trace of heptane in the stream. The trap liquid was initially analyzed by low-temperature fractional distillation, yielding a bottom product of heptane and heavier components. Specific gravity of this fraction was measured and the mol weight was determined by freezing point depression. The hexane and lighter overhead gas collected during distillation was

Jan 1, 1967

By L. H. Wolfe, J. O. Thieme

The tensile and shear strengths of frozen soil and the compressive strengths of ice and frozen soil were measured These tests showed that the strength of ice and of frozen soil increased as the temperature was decreased. A method is presented for the measurement of thermal conductivity and specific heat of earth and ice. Using the parameters thermal conductivity, specific heat and density, the thermal diffusivity was calculated. Results of thermal property tests establish water as the major variable contributing to thermal property values. INTRODUCTION During the last few years, a considerable interest has developed in underground storage of cryogenic liquids such as liquid natural gas. Because of its economy and safety, underground storage is being used. Before underground caverns could be seriously considered, however, the strength and thermal characteristics of the soil had to be obtained. Methods are presented herein for measuring tensile, compressive and shear strengths, specific heats and thermal conductivities. The tests were performed on two soils at temperatures from ambient to -195C. Because it is difficult to find data in the literature that pertain to a particular soil, ice was used to evaluate some of the experimental procedures. Some of our ice data are compared to the published values. PHYSICAI, PROPERTIES SOIL SELECTION AND SAMPLE PREPARATIONS The physical and thermal tests were performed on recompacted samples of two typical soils that might be encountered in underground caverns: a gray, fire clay that had been pulverized and mixed with water into a pliable mud; and a brown, sandy silt that was dug from a dry river bank and subsequent1y mixed with tap water. The moisture content of the test specimens was 17 to 22 per cent (based on the wet weight of the soil). This amount of water nearly saturated the silt, but the clay was well below the saturation point. The soils are partially described by the sieve analyses (Fig. 1) which show the particle size distribution of the coarser than 44 micron portion of the soils. The curves also show that approximately 75 per cent of the clay and 20 per cent of the silt are finer than 44 microns. The particle size seems to affect some of the physical properties, and the finer than 44 micron portion is important in the thermal conductivity test. Physical specimens were cut with a bandsaw from frozen blocks of the soil. The specimens were then sanded and kept frozen until after they were tested. They were stored at the approximate test temperature for several hours. Then, to equilibrate the samples to the test temperature, they were stored at the temperature at least an hour before they were tested. TEST EQUIPMENT An Instron, which is a constant speed gear-driven instrument, applied force to the samples and drew stress-strain plots of the test. All of the samples were strained at 0.02 in./min. All of the physical tests were conducted in a temperature cabinet that controlled the temperature to 1C. TENSILE STRENGTH Briquette or "dog-bone" shaped tensile specimens were pulled with suitably shaped jaws. This shape of specimen, which is usually used for tensile tests of concrete, is satisfactory for testing frozen soils except that the samples tend to

Jan 1, 1965

By R. V. Smith, M. H. Cullender

Rigorous equations for calculating subsurface pres.tures in flowing and static gas wells, and pressures along horizontal pipelines are presented in this paper. These general equations, based on the mechanical encrgy balance, contain no assumptions regarding temperature and can be used with either straight-line or curved temperature gradients. The friction factors recommended in this report are based on an absolute roughness of 0.0006 in. Flow is always considered to he turbulent. Although these general flow equations are solved by numerical means, the methods are as convenient as many of those used today to calculate pressures in gas wells and pipelines. These numerical methods are illustrated for flowing and static columns of gas in wells and for flow in pipelines. The authors believe that the use of these methods will result in more accurately calculated pressures in gas wells and pipelines and that [he methods are more flexible than any in use at the present time. Also, the methods are easily adaptable to automatic computers. INTRODUCTION The problem of calculating subsurface pressures in gas wells and along flowing pipelines has been studied by many investigators with the results that two widely-used methods have been presented in technical literature. One approximates conditions in a gas well or pipeline by assuming that temperature and compressi- bility are constants tor the entire gas column. The other method assumes that temperature is constant at some average value but permits compressibility to vary with pressure at the constant temperature. While these approximations may be justified for relatively shallow and low-pressure wells, they become more and more unrealistic as the calculations are applied to the deeper high-pressure wells. In wells or pipelines where there is an appreciable change in temperature between the inlet and outlet ends of the flow string, it is equally as important to consider the change in compressibility with temperature as it is to consider the change with pressure. Consequently, the methods presented in this paper were developed to provide a workable procedure for calculating pressures in gas wells and pipelines which makes no assumptions regarding either temperature or compressibility. The friction factors recommended for use with the equations were the result of study of published literature and extensive experimental work with gas wells. The F,. values in Appendix IV which contain the friction factor are equivalent to an absolute roughness of 0.0006 in. as calculated by the Nikuradse equation. The absolute roughness of 0.0006 in. was previously published by Cullender and Binckley2 for gas wells; it is in agreement with 0.00065 in. of Smith, Dewees, and Williams1 for gas wells and 0.00065 in. of Baker and Swerdloff4 for pipelines; it is in reasonable agreement with 0.0007 in. reported by the Bureau of Mines and the American Gas Assn.5 for large transmission lines. However, any type of friction factors may be used in the equations by a change in the F, values in Appendix [V].

Jan 1, 1957

By J. W. Wolfe

The Benedict-Webb-Rmbin equation of state was used in digital computer programs to make rapid determinations of natural gas equilibrium phase compositions. Mixture components were the nine hydrocarbons, methane through heptane. Comparative tests were made with data from field separators, plant exchangers and literature sources. For most tests made, the B-W-R equations predrcted component K-vallles, mixture enthalpies and phase aensities within tile accuracy of available correlations for rune-component mixtures: significant deviations were noted when mixtures contained components heavier than heptune. INTRODUCTION Rapid, accurate methods for calculating the properties of hydrocarbon mixtures are needed for predicting the behavior of natural gas processes. For example, because of the common use of two-phase separation methods in gas proccssing, calculations must be made for the density and composition of equilibrium vapor and liquid phases. In addition, predictions must often be made for the temperature and pressure changes associated with the flow of single- and two-phase natural gas fluids through heat exchangers, chokes and turbines. The complex nature of non-ideal mixtures makes rigorous calculations for these natural gas properties difficult and time-consuming. Natural gas properties have conventionally been determined with generalized data correlations. These correlations, such as the K-values presented in the NGAA Equilibrium Ratio Data Book, are usually based on some average property of the mixture and have proven to give satisfactory estimates for the properties of most oil and gas systems. Digital computers, however, have not only increased the utility of these correlations, but also have made practical the use of equations of state for calculating the properties of natural gases. Potentially, if basic pure component data are available, all the desired thermodynamic properties of mixtures can be derived from a single equation of state describing the mixture. Currently. the best equation of state for describing the behavior of pure, light hydrocarbons over single-and two-phase regions, both below and above critical pressure, is that of Benedict. Webb and Rubin.' The equation contains eight constants which are determined by empirically fitting the equation to pure-component PVT data. Benedict and his co-workers presented the clevelopment of their equation as early as 1940.2 The B-W-R equation is a modification of the Beattie-Bridgeman equation of state and gives better PVT predictions than the Beattie equation at densities above the critical density. Benedict et al. made the equation applicable to mixtures by relating the mixture constants to composition and the constants for pure compounds. These authors, in a serics of articles in 1951,' showed the accuracy of the B-W-R equation when applied to light hydrocarbons and their simple mixtures, and demonstrated the equation's utility for making vapor-liquid equilibria calculations. Since its inception, the B-W-R equation has been the subject of many papers. Numerous additions and modifications have been proposed,'-',"'-':' the purpose of most being to improve the equation's ability to predict PVT properties of specific mixture systems. In 1959, Opfell, Pings and Sage summarized much of the work on the B-W-R equation in their book on equations of state." Since natural gases are primarily mixtures of light hydrocarbons, the B-W-R equation is a natural choice for an equation of state to compute both phase equilibria and thermodynamic properties of natural gases. Consequent-

Jan 1, 1967

By David Cornell, D. L. Katz

Procedures for computing turbulent flow of gas in steady state near the well bore and a graphical method for predicting unsteady state laminar flow at distances from the well have been combined to compute pressure gradients in gas reservoirs. Methods are discussed for predicting single and multiple transients at constant flow rate in an infinite reservoir, predicting constant How rate. in a finite reservoir, reproducing and interpreting back pressure test data, and prediction of the behavior of a closed-in gas well. An example of the graphical method is given for a single transient and the results are compared to a published analytical solution. A typical calculation of the pressure gradient in a reservoir and the production of gas is made .starting with data from a hack pressure test. INTRODUCTION The calculation of pressure gradients throughouse a ga-reservoir at any point in its production history is a complex problem involving turbulent flow near the well bore and unsteady state flow of the gas from the reservoir. The problem is complicated further by the variety of boundary conditions that may be imposed upon the flow. Such boundary conditions include finite or infinite reservoirs. constant or variable rates of production, constant or varying bottom hole pressures, complex initial pressure distributions throughout the reservoir, and pressure maintenance through cycling. In addition to these problems there are the practical difficulties associated with natural gas reservoir analysis These might include: (1) variation. in the permeability and porosity throughout the formation. (2) variations in the thickness of the formation, (3) communication of the producing formation with other producing zones. (4) radial variations in permeability due to influx of drilling fluid. acidizing. or the buildup of liquid in the formation around the well, (5) partial penetration of the producing zone. (6) water flood, and many others. These problems will be considered eventually. In the mean- time, it is necessary to present methods of handling the case of radial flow through a homogeneous, regular, producing stratum to a single well uncomplicated by other factors. The idealized case of steady state, radial flow from a natural gas reservoir was studied by Elenbaas and Katz5 and by Mac-Koberts.' An analysis of the unsteady flow of gases through porous media has been made by Aronofsky and Jenkins' for the one dimensional laminar case with the boundary conditions of constant downstream pressure and uniform initial pressure. Solutions for the partial differential equation for unsteady state radial flow of liquids through porous media have been given by Van Everdingen and Hurst" for the constant bottom hole pressure and constant production rate cases for laminar flow in a reservoir initially at a constant pressure. The analogous heat transfer equation has been treated by Perry and Berggren9 for the case of quenching a cylindrical hole in an infinite medium to a constant temperature. The back pressure curve, which involves the variation of the bottom hole pressure as it is determined by the behavior of the reservoir as a whole. has been studied by Rawlins and Schellhardt,10 Binck-lev,' Baumel and Breitung,' and others. No complete. adequate treatment has been given for the unsteady state, radial flow of gases into a cylindrical hole from a porous medium with or without the presence of turbulent flow and for any set of boundary conditions. If one relies only on analytical solutions of the basic equations, a new solution must be obtained for each new boundary condition. Furthermore, the complexity of the mathematical expressions for boundary conditions other than the very simple cases limits the use of analytical solutions. Graphical methods exist, however. that are general in nature, accurate, and readily employed. The analysis of steady and unsteady, laminar and turbulent, radial flow for various boundary conditions with application to specific natural gas well problems by means of graphical procedures forms the content of this paper. STEADY STATE RADIAL FLOW EQUATION FOR LAMINAR AND TURBULENT FLOW A steady state radial flow equation for laminar and turbulent flow through unconsolidated sands has been given by Muskat as discussion of the work by Elenbaas and Katz.5 This equation is based on the properties of unconsolidated

Jan 1, 1953

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970

By L. G. Jones

The objective of this paper is to present a better method for determining reservoir reserves of gas from extended drawdown test data. Means are proposed for conducting the test with or without the requirement of prior build-up to a static reservoir pressure. The test is based upon solutions to unsteady-state gas flow developed in an earlier paper and is believed to be superior to previous methods for determining reserves from flow-test data. The limitatiom of such tests are discussed, with emphasis on the subject of test duration. The test can supply useful information if the mechanics and limitations of it are understood. This type of information can be particularly valuable in deciding whether or not to drill offsets to producing wells. INTRODUCTION In the life of any gas or oil field, the engineer in charge must decide if certain wells need stimulation and whether or not to drill more wells in the field. Well production tests can aid in making these decisions. For instance, the value of stimulation can be estimated by determining formation permeability and the amount of well damage from drawdown or build-up tests. These tests are freely used in the oil industry. A less well-known type of analysis uses data from extended drawdown tests for determining reservoir reserves. Into this category fits the "reservoir limit test" of Jones,'-3 and the "reservoir reserve test" presented in this paper. The chief value of these tests is their ability to "see" reservoir boundaries and thus prevent costly dry holes. Methods for determining reserves of both gas and oil reservoirs were proposed in the work of Jones. His proposals appear sound for single-phase liquid flow. The generality of his work on gas flow is open to question, however, since the analysis is based upon solutions to liquid-flow equations. The test for gas reserves presented here is derived from unsteady-state gas-flow solutions presented in a previous paper.' In deriving this reservoir reserve test, it has been assumed that all of the reservoir energy is provided by expansion of the reservoir fluid and that the reservoir is radially symmetrical. No mention has been made of the effects of deviations from Darcy's law such as quadratic (turbulent) flow and well damage. Since the terms representing these effects are function of rate only, since the rate is constant and since the data are plotted in a differential form, it seems reasonable that these terms will not affect the test results in any way. No attempt has been made to estimate the effects of water drive, faults, or other departures from an idealized and isolated reservoir. For a qualitative discussion of some of these effects, the reader is referred to a paper by Jones3

By C. E. Turner, J. R. Elenbaas, R. D. Grimm, J. A. Vary, D. L. Katz

A technique is presented by which well samples and core plugs of dolomite formations are classified by microscopic examination into seven different porosity grades. Quantitative values of porosity and permeability are determined for each grade by a statistical correlation of the core plug test data with the porosity grading system. These quantitative values are applied directly to the grades exhibited in the well samples for the purpose of estimating the reservoir void space for wells that were not cored. The procedure is described for estimating the gas reserves per unit area lor the South Hugoton gas field. but a reserve estimate for the field is not given. INTRODUCTION The miscroscopic examination of well sample; and the graphic recording of their lithologic qualities and other distinguishing characteristics of various geologic formations drilled is both a science and an art of long standing and wide application. Usually the primary objective of a geologist who "sits on the well" and examines the samples are: to identify the formation being drilled, determine the total depth, casing point. and completion interval. In most cases the porosity is described. if done at all, in general terms. such as: trace, scattered, fine, poor, fair, medium. good, excellent, or in some other relative terms. In fields where various geologists have examined samples and recorded observations on many wells considerable variations in lithologic terms and porosity descriptions occur unless there is primary effort to establish uniformity of logging observations and standards of recording observable porosity. When an estimate of the pore volume of a reservoir is made a geologic concept of the processes that control the magnitudes of porosity and permeability is developed by microscopic examination of well samples. The characeristics and appearances are then mentally related to rather general quantitative units of porosity based on physical core data from the same reservoir or on such data or experience in other reservoirs that have similar qualities. The reliability of such estimates depends largely on the variations of the lithology of the formations, the geometric properties of its void system. the extent of comparisons of sample appearances with porosity data, as well as the uniform recording of all relevant characteristics. This statement is particularly significant for dolo-mitized limestone formations of substantial thicknesses and heterogeneity such as the Permian Dolomites of the Hugoton gas field. Jn this field, as well as in most of the Permian Dolomite fields, the producing formations are of relatively great thicknesses in which the porosity and permeability of the reservoir varies substantially in all directions, depending on the crystalline structure. degree and kind of impurities, kind of fossils anti cementation thereof, degree of dissolution. and fracturing. The variations of the lithologic texture of the dolomites and post deposition alterations have resulted in porosities and permeabilities of such magnitudes that only a part of the gross thickness can be counted as "pay." At the time of this study insufficient gas production had been experienced to apply the pressure decline production method in the South Hugoton Field and the electric logs are not definitive enough. The problem of estimating gas reserves in the south part of the Hugoton Field is primarily one of determining the pay thickness and porosity from well samples and core data. The area studied embraced all that part of the field lying south of an east-west line through Guvmon, Okla., anti containing approximately 1.000,000 acres. This paper describes a technique of correlation of physical core data with well samples so that quantitative values of pay thickness, porosity. Permeability, and connate water may be assigned to well samples that are representative of a given interval, and thereby permitting the estimation of gas reserves lor a given unit area. The procedure was developed by a uniform microscopic qualitative porosity grading of the dolo. mite core plugs.. and relating these grades to the respective physical core data on a statistical basis The well sample-were also graded in a Similar manner in order that the quantitative values established lor the core plugs could be applied to the well Sample for wells that were not cored. GRADING OF DOLOMITE A group of experienced geologists was given the assignment of examining the samples on all wells in South Hugoton in order that they could log their observations in a uniform and standardized manner and grade the observed porosity so that it could be related quantitatively to the core data. The group initiated the study 011 chips from cores which bad been tested for porosity and permeability. This study continued until all of the geologists developed a common knowledge of lithologic terms and of the characteristic appearances of the samples and their relations to measured porosity. The characteristic appearances of the dolomite samples under twelve-power magnification as related to their qualitative porosities afforded a classification of the dolomite into :even grades of porosity, ranging from dolomite of no-visible porosity under twelve-power magnification to dolomite of excellent porosity. The assigned grade for a specific 10-ft interval is a weighted average of all visible grades of porosity exhibited by the cuttings representing that interval. The porosity characteristics were recorded by a color graph adjacent to the lithology column in conjunction with a numerical system for further definition of relative porosity as shown in Fig. 1. The three vertical lines to the right of the lithology column each represent 33 1/3 per cent. which lines were used to record the percentage of the samples, for any particular interval. that showed porosity under the microscope. The colors were used to denote actual pore size. i.e., orange. blue and I-ed for pore diameter of one-fourth or less. one-fourth to one-Ilalf. and greater than one-half millimeter, respectively. The area colored 1)). one or more colors represents the percentages of the samples exhibiting pores of the respective size or sizes. The numerals from one to six inclusive shown on the log in

Jan 1, 1952

By E. G. Woods, A. G. Comer

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