Search Documents

By L. K. Nemeth, H. T. Kennedy

The investigation presented in this paper was performd to develop a relationship between the dewpoint pressure of a hydrocarbon reservoir fluid and its composition, temperature and characteristics of the heptanes-plus fraction. This is part of an over-all program to develop correlations for the essential properties of reservoir fluids with the composition, temperature and pressure. In all, 579 dewpoint pressures, determined experimentally on 480 different condensate systems along with their compositions and temperatures, were obtained from the literature and from previously unpublished sources. Compositions were expressed in terms of mol fractions of methane through heptanes-plus, nitrogen, carbon dioxide and hydrogen sulfide, and the molecular weight and specific gravity of the heptanes-plus fraction. Distribution of the samples assured a good representation of most of the major oil-producing areas in the world. A multiple regression analysis program written for the IBhl 7094 digital computer was used to compare numerous equations and to calculate their coefficients. The final equation containing 11 coefficients resulted in an average absolute deviation of 7.4 percent for all the samples. Segregation of data on the basis of one or more variables did not materially improve the results. The correlation developed is useful in predicting the dewpoint pressure of condensate reservoir fluids from variables normally measured in the laboratory. INTRODUCTION In any discussion of the dewpoint of hydrocarbon mixtures, it should be noted that there are two kinds of dewpoints which must be distinguished. The first type, or normal dewpoint, usually occurs at low pressures as dry gas is compressed to the point where liquid first forms. The second type, sometimes called the retrograde dewpoint, occurs when a gas mixture containing heavy hydrocarbons in solution is depressured until a liquid is formed. The first type, as applied to condensate mixtures, would normally be below atmospheric pressure and would be of no interest in reservoir performance. The second type is of great significance in interpreting reservoir behavior and is the property referred to throughout this paper. Laboratory measurement of dewpoints is laborious, expensive and subject to many errors, and it is desirable to find a reasonably accurate method to predict dewpoint pressures from the composition and other readily available properties normally measured on the fluids. A possible method of calculating dewpoint pressures would be to employ a set of equilibrium ratios (K-values). 1 However, the K-values themselves depend on composition as well as temperature and pressure; presently their use adds nothing to the accuracy of predicting dewpoints. If future correlations provide a significantly better method of estimating K-values for heavy hydrocarbons, they may be of greater use. REVIEW OF LITERATURE In 1944, Olds, Sage and Lacey 2 studied the behavior of fluids from the Paloma field. Properties of seven mixtures were studied at temperatures ranging from 100 to 250 F and at pressures up to 5,000 psia. From the information obtained, the influence of composition and temperature upon the dewpoint was determined. It was found that removal of the intermediate molecular weight components from the mixture resulted in considerable increases in the dewpoint pressures. Results indicated that the effect of changes in temperature was relatively minor when compared with the effect of modifying the mixture by removing the intermediate components. Thus, removal of propane and butanes resulted in an increase in dewpoint pressure more than 10 times greater than a change in temperature from 100 to 250 F. The authors stated that althouah

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. Husson, R. Iffly, M. Gondouin

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

By 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 David Cornell, Yusuf K. Sukkar

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

Jan 1, 1956

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

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

By K. L. Young

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

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

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

Jan 1, 1952

By F. R. Scauzillo

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

By P. A. Witherspoon, S. P. Neuman

Evaluating the permeability of a caprock overlying a potential gas storage reservoir is a very critical problem. Pumping water from the reservoir can be used as an evaluation tool in analyzing this problem. Fluid level changes that occur in the aquifer us well as in the caprock can be measured with appropriately placed wells. If the leakage of water from the caprock into the aquifer is considerable, the effects will be apparent in the aquifer. If the leakage is slight, however, it will not be possible to detect it with certainty from observations in the aquifer alone. Fluid level measurements in the caprock must be relied upon. and improved methods of analyzing such effects have been developed which are based on a theoretical analysis of fluid flow through a caprock of infinite thickness. An example applying these methods to field data is discussed. INTRODUCTION One of the most critical problems in evaluating an aquifer gas storage project is determining the tightness of the caprock overlying the formation to be used as the storage reservoir. A formation that has previously held oil or gas obviously has a suitable caprock, but an aquifer that contains only water gives no such assurance. A number of aquifer projects in the United States have been troubled by gas leaking out of the intended storage zone, and the ensuing difficulties have led to the development of new evaluation methods. One of these new methods is pump testing wherein water is removed from the aquifer at some controlled rate prior to injection of gas. This fluid withdrawal causes a pressure drop to move out through the aquifer for considerable distances in a matter of days or weeks. Depending on the properties of the caprock, a pressure transient can also pass upward (as well as downward) through the caprock layers adjacent to the aquifer. Thus, if the operator has placed observation wells at appropriate distances from the pumping well, the rapidity with which the pressure transients reach different points in the system can be used to investigate the fluid transport properties of both the aquifer and its caprock. The usefulness of pump testing has been recognized by groundwater hydrologists for many years as a means of determining the potential yield and properties of aquifers used in water supply They have introduced the term "leaky aquifer" for a system in which an aquifer is overlain (or underlain) by semipermeable caprock layers. The ease with which water leaks into the aquifer during pumping can, of course, be very beneficial in bringing additional water to the pumped well. Hydrologists have therefore devoted considerable attention to this prob-lem. From the gas storage standpoint, however, the tighter the caprock layers that overlie the intended storage reservoir, the better are the conditions for minimizing or eliminating any vertical migration of gas. Thus, after a suitable geologic structure has been found, the emphasis in aquifer storage projects is in determining that the caprock is tight. Attention has recently been focused on the use of pump testing as one approach to solving this problem.23,24 This paper presents a further development on evaluating the permeability of a slightly leaky caprock when the caprock is of infinite thickness. From the practical standpoint, this means that the caprock layers are thick enough that pressure transients do not reach the outer boundaries of the system during the pumping test. In a subsequent paper, an analysis of the case where the caprock is of finite thickness will be presented. PREVIOUS WORK ON LEAKY AQUIFERS Jacob" developed a partial differential equation describing the flow of water in an aquifer of permeability k that is overlain by a leaky caprock of permeability k'. Fig. 1 shows a schematic cross-section of the system under consideration. One of his principle assumptions was that if k > > k', the direction of flow is essentially vertical in the caprock and horizontal in the aquifer. Neuman19 confirmed the validity of Jacob's assumption using a mathematical model. Another assumption was that a permeable source layer overlies the caprock (Fig. 1) and is able to maintain a constant hydraulic head at the upper boundary of the caprock. By neglecting the effects of compressibility within the caprock, Jacob1* developed a solution for a bounded circular aquifer. Later, Hantush and Jacobx5 used the same assumptions to solve the case of an infinite radial aquifer that is pumped at a constant rate. Their solution may be expressed in di-mensionless parameters by

By Donald L. Katz, Leslie J. Noaker

Experimental data are presented for hydrate formation conditions for mixtures of hydrogen sulfide and methane. Vapor-solid equilibrium constants for hydrogen sulfide are also presented. Equilibrium constants, defined here as the ratio of the mole fraction of hydrogen sulfide in the vapor phase (dry basis) to the mole fraction of hydrogen sulfide in the solid phase (dry basis), were calculated using the vapor-solid equilibrium constants for methane presented in an earlier paper. The conditions at which a natural gas will form gas hydrates with water may be computed from its composition, when hydrogen sulfide is present in addition to carbon dioxide and the paraffin hydrocarbons. INTRODUCTION Many gases are known to form solid hydrates in the presence of liquid water at temperatures above 32°F. This hydrate formation is of importance to the natural gas industry inasmuch as most of the constituents found in natural gases will form hydrates. The attention of the industry was called to hydrates by Hammer-schmidtl0 in 1934. A number of papers have been written giving the temperatures and pressures at which these hydrates will form3,4,5,6,7,8,9,10,13,14,15,18,19 and a bibliography on the subject has been compiled13. Deaton and Frost4,5,6,7,8 report work for various pure gases and for some natural gases. Katz and co-workers3,4,13,14,18,19 have reported data for various gases at high pressure. They developed a method of predicting conditions for hydrate formation from equilibrium constants3,19 and from gab gravity." The effects of anti-freeze agents on the formation of hydrogen sulfide hydrate were reported by Bond and Nelson.' The hydrogen sulfide-water system, including hydrate, has been the subject of recent investigation by Sage and co-workers." Mixtures of hydrates behave as solid solutions and an equilibrium exists between the vapor phase and the solid phase that appears to be comparable to the equilibrium between vapor and liquid phases.3,10 By means of the vapor-solid equilibrium constants which have been published for methane, ethane, propane, iso-butane, and carbon dioxidel6 it is possible to determine the temperature and pressure at which hydrates will form for a specific gas mixture when water is present. Hydrogen sulfide which is frequently found in natural gases has a very strong tendency to form hydrate and must enter the solid solution when present in a natural gas forming a hydrate. To determine the effect of hydrogen sulfide when computing the conditions for gas hydrate formation, it is necessary to have the equilibrium constant giving the ratio of the concentration of hydrogen sulfide in the gas phase to the concentration in the solid phase. This paper presents data on hydrate formation for mixtures of hydrogen sulfide and methane and presents the equilibrium constants which were calculated from these data. EXPERIMENTAL PROCEDURE The apparatus used to determine the conditions of hydrate formation for the methane-hydrogen sulfide system was similar to that employed by previous investigators.3,18,19 A charge gas cylinder was used to supply the hydrogen sulfide-methane mixtures, called the charge gas, and to supply the pressure up to 1.000 psia. This

Jan 1, 1955

By R. J. Burgess, A. R. Ramey, A. R. Adams

During interpretation of pressure buildup tests on gas wells in a tight dolomite gas reservoir, peculiar behavior was noticed. Two straight lines were apparent. Effective permeability to gas taken from either straight line was about the same, and the Miller-Dyes-Hutchinson dimensionless time check for the straight line was proper for both straight lines. Geological data indicated the likelihood of scattered trending fractures in the reservoirs. Since the first straight Iine yielded permeability values close to the geometric mean permeability from core analyses, it was postulated that the reservoir model was that of an acidized well completed in the tight dolomite, but that widely scattered hairline fractures caused the mean permeability of the reservoir distant from the well to be higher than the matrix permeability. Because all other studies of fractured reservoirs to the authors' knowledge assumed that the fracture matrix was dense enough to communicate directly with the well, no interpretative methods were available. The Hurst line-source solution for a radial change in permeability for interference between oil reservoirs was adapted to pressure buildup testing. The result indicated that the first straight line should yield the proper matrix permeability and wellbore skin effect. The second straight line may be extrapolated to obtain static pressure. The time of bend between the straight lines was used to estimate distance to a fracture. Application to field test data is shown. It is believed that the methods developed and the case history presented will add to present tools available for pressure buildup interpretation. Introduction Since the pioneer studies by Miller, Dyes, and Hutchin-son1 and Horner' in 1950 and 1951, well test analysis has become recognized as one of the most powerful tools available to both production and reservoir engineers. Well test analysis serves as a logical basis for well stimulation and completion analysis, and for long-term reservoir engineering. Since the early 19501s, much effort has been placed on the development of well-test analytical methods. Reservoir and well conditions of increasing complexity have been considered systematically to provide the analyst with a catalog of causes and effects. Matthews and Russella state that some 200 papers dealing with this subject have been published in the last 35 years. Developments in well test analysis appear to have originated in one of two ways. Either a physically realistic field condition was anticipated and analytical solutions for the condition achieved, or anomalous field test behavior was recognized and interpretative methods sought for the anomaly. In recent years, it has appeared that the latter has inspired an increasing number of studies. The analyst today finds an increasing number of known cause and effect studies available for well test analysis, the classic of which is that of finding the specific flow problem that generated the answer — the well behavior. Although it may be impossible to achieve this goal uniquely, the analyst often is able to select a useful interpretation that combines all known performance and geologic data — or to show that various logical alternatives would not significantly affect the interpretation. During a recent reservoir study, we observed gas well test behavior that did not appear to fit behavior described previously. Although it cannot be said that we have found a unique interpretation, we shall present in this paper the peculiar behavior observed, and describe the reservoir and interpretative methods developed. Reservoir Description The subject gas reservoir is a 9-mile-long, narrow dolomite reservoir lying within a limestone of Ordovician age. (See Fig. 1.) The dolomitized rock in the field consists of dark brown to buff, dense to coarsely crystalline, vugular dolomite containing numerous hairline fractures, many of which may have been closed in the reservoir and parted when cores were brought to the surface. Larger fractures are also apparent in core, but usually are filled and sealed with euhedral dolomite crystals. Portions of the north flank of the reservoir are known to be cut by a sealing fault downthrown to the north. Gas wells located near the fault have higher open flow potentials than those more distant from the fault. This is believed to be a result of higher permeability near the fault due to more extensive and open fractures. Detailed coring and core analysis have been performed on several of the wells in this reservoir. Fig. 2 presents permeability variation' plots for both horizontal and vertical

Jan 1, 1969

By R. A. Wattenbarger, H. J. Ramey

A systematic study has been made of the application of the real gas pseudo-pressure m(p) to short-time gas well testing. The m(p) function can be used in real gas flow problems to account for the variation of viscosity and z-factor with pressure. A mathematical model was solved numerically to generate solutions of various real gas flow problems. The effects of turbulence, formation damage and wellbore storage were included in the model. The analysis of simulated well tests showed that the interpretation methods normally used for liquid flow are generally accurate when the m(p) is used. For practical rates without turbulence, the solutions are not rate-sensifive, and the flow capacity.. kh and skin effect can be determined accurately from either a buildup or a drawdown test. However, the kh calculated from a drawdown test can be significantly low when turbulence is present. A case was simulated in which this error wa.s 36 percent. Turbulence does not affect the determination of kh from buildup tests. The proper determination of kh from p and pa buildup and drawdown plots is developed by analyzing their relationship to the m(p) method. A simple equation is given that can be used for long-range gas well performance forecasting. This expression is compared with a method pre-sended by Russell et al." Introduction The formation flow capacity and wellbore damage condition can be determined for liquid-producing tests by means of buildup and drawdown tests. These tests make use of short-time pressure transient data rather than stabilized flow tests that are often used for gas well testing. Tracy' presented a method of gas well testing that was based on ideal gas equations and utilized p2 buildup plots. This method was shown to be good for low pressure wells. MatthewsZ suggested plotting p and using an average slope and average gas properties. This method was more successful on high pressure wells. Al-Hussainy, Ramey and Crawfordh howed that variation in gas properties could be simplified by using the real gas pseu do-pressure m(p). Some cases showed that the use of the m(p) function Original manuscript received in Society of Petroleum Engineers office Aug. 8. 1967. Revised manuscript received June 11, 1968. Paper (SPE 1835) was presented at SPE 42nd Annual Fall Meeting held in Houston Tex., Oct. 1-4. 1967, and at SPE 4th Annual Eastern Regional ~eetin$ held in Pittsburgh. Pa., Nov. 2-3. 1967. 0 Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. 'References given at end of paper. provided an accurate method of interpreting buildup and drawdown tests. Because the work of Al-Hussainy et al. included only a few actual cases, there was a need to explore the use of the m(p) method for a greater variety of flow conditions. The occurrence of turbulent flow around the wellbore often is an important factor in gas well testing. Swift and Kiel,' and Carter et a1.' treated turbulence in gas well testing, but did not consider the variation of viscosity and z-factor with pressure. This paper presents the results of an investigation of the application of the m(p) method to buildup and drawdown testing. Vhe investigation included the effects of turbulence, formation damage and wellbore storage. Gas Flow Equations To formulate the mathematical model, many of the assumptions usually used in well testing theory are applied. The system has radial geometry with a closed outer boundary and is composed of a horizontal porous formation that has uniform and isotropic rock properties and uniform thickness. Allowance is made, however, for a radial region of reduced permeability near the wellbore. This region represents formation damage. The geometry of the system is shown in Fig. 1. Darcy's law does not always apply to gas flow. A more general expression is needed for non-Darcy, or "turbulent" REGION OF DAMAGED PERMEABILITY SEALED BOUNDARIES ClRCULAR WELL BORE Fig. 1—Radial flow model. AUGUST, BT7

Jan 1, 1969

By Jens P. Nielsen, W. H. Justice

This paper describes an improved acidizing technique which has been applied in acidizing gas wells in the La Gloria Field. Wells acidized in this manner exhibited a much greater increase in deliverahility than wells acidized in the conventional manner. The acid is injected into the well in small slugs separated by a small volume of high pressure gas. The acid and gas are displaced from the tuhing by high pressure gas to reduce the hydrostatic head on the well when the acid is recovered from the well. The well is opened for production as soon as the acid and gas are displaced. This procedure permits a very rapid recovery of the acid which results in large increases in the deliverability of the gas wells. It has been found advantageous in the acidization of gas injection wells to inject the high pressure gas directly after the acid without backflowing the acid out of the well. This practice has made it possible to inject the gas in the injection wells of this field with a much lower pressure differential. This reduces the horsepower required to inject the gas and also decrease:, the number of injection wells required per reservoir. INTRODUCTION The results obtained by acidizing gas wells ill the La Gloria Field have varied considerably in the past. For this reason. a planned study of the acidizing procedure was undertaken in order to obtain more consistent results. The Frio Sands of the La Gloria Field are rather permeable and it was not necessary to acidize many of the wells when they were originally completed. At the completion of an eight-year cycling program in several of the gas reservoirs, it became necessary to recomplete many wells in other reservoirs. It was found that the opan Flow potential of the wells recompleted in different sands was not as large as wells that were originally completted in that reservoir. A study of the drilling records of several of these wells indicated that the sands had originally been subjected to considerable hydrostatic drilling fluid overload while the well was being drilled to the deeper sands of the field. It is indicated that during this eight-year period, the drilling mud and water that were lost into the formation when the well was being drilled had formed a relatively impermeable Eone around the well bore. This impermeable zone around the well bore made it very difficult to obtain high deliverability wells when the original wells are recompleted in other zones. Muskat' has shown that in uniformly permeable and porous formations a zone of partial plugging one ft in radius around the well would drastically decrease the productivity of the well. The analysis for this problem shows that the production capacity of a well is very sensitive to the value of the permeability of the zone immediately surrounding the well bore. This is to be expected in view of the highly localized character of the pressure drop in a radial flow system about the well center. The reduction of the deliverability is not proportional to the thickness of the mud-affected zone as the major portion of the reduction in deliverahility IS caused hy the

Jan 1, 1952

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