Search Documents

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

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

Jan 1, 1966

By M. H. Cullender

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

Jan 1, 1956

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

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

Jan 1, 1955

By G. Thodos, D. E. Matschke

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

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

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

Jan 1, 1955

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

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

Jan 1, 1952

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

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

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

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

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

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

Jan 1, 1969

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

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

By R. V. Smith

The theory for unsteady radial flow of gas, as simplified by Aronofsky and Jenkins, has been reviewed and compared with actual well performance. This study indicated that the theory should be modified by the introduction of an empirical "rate of flow" function "Y(q)". The expansion of the theory on the flow of gas to include a Y(q) term bridges the gap between the theory of unsteady-state flow of gas and actual gas-well behavior. Apparently, the Y(q) term is a function only of the rate of flow for a given well. The completion factor or skin effect is associated with the Y(q) function in such a manner that at least two or more sets of drawdown or pressure build-up tests are needed to separate the Y(q) function and the completion factor. Since the Y(q) function used in this report represents an energy loss dependent on rate of flow near the well-bore of a gas well which is in addition to the loss required by Darcy's law, the Y(q) function is related to the exponent of the back-pressure curve for a gas well. Techniques and equations are presented which permit the estimation of stabilized and "short-time" deliver-abilities. INTRODUCTION The published theory of radial unsteady-state flow of gas through porous media may be divided into papers1-" which use Darcy's law as a premise and other papers4-7,16-10 which recognize that Darcy's law may not hold for flow into gas wells. Miller8 compared theory based on Darcy's law with actual gas-well behavior and concluded that flow into gas wells deviates from Darcy's law. Houpeurt4,5 in a thorough investigation of flow into gas wells concluded that the flow rate into a gas well is proportional to the difference between the squares of upstream and downstream pressures elevated to a power between 1.0 and 0.5. Houpeurt also concluded that the deviation from Darcy's law is caused by irreversible, kinetic energy exchange between the flowing gas and the porous media. Tek9 presented a complete evaluation of the flow problem in gas wells. Several investigators state that the inclusion of gas compressibility into flow equations based on Darcy's law does not begin to explain the flow problems suggested by gas-well behavior. For many years, natural-gas engineers have realized that gas wells by actual test exhibit behavior which cannot be explained entirely by Darcy's law. Many tests on gas wells have indicated that a value of 0.85 for the exponent in the back-pressure equation is a fair average for many wells. However, Darcy's law either for steady or unsteady-state conditions indicates that the exponent should be 1.0. The purpose of this paper has been to examine this deviation from Darcy's law and to present a reasonable mathematical means of evaluating the deviation. The method presented in this paper combines previous theory with an empirical function to provide a description of gas-well behavior. All equations in this report are for constant terminal-rate-of-flow conditions. While present methods of testing wells usually are not based on this concept, they are believed to approach the constant terminal-rate conditions closely in many tests. However, several gas wells were tested by holding the rates of flow constant in order to approach more closely the constant terminal rate of flow; the results confirm the concepts presented in this paper. REVIEW OF THEORY FOR UNSTEADY-STATE FLOW OF GAS The basic partial differential equation for the flow of natural gas has been given in technical literature; the equation for a cylindrical reservoir with the well in the center is

By David B. Burrows, Rlkl Kobayashi, Norman L. Carr

The viscosity of hydrocarbon mixtures, whether in the gas or liquid phase, is a function of pressure, temperature, and phase composition. This paper presents methods for the prediction of the viscosity of the gas or less dense fluid phase over the practical range of pressure, temperature, and phase compositions encountered in surface and subsurface petroleum production operations. The correlation necessary to predict the effect of pressure on viscosities is presented in Part I. Serious discrepancies in high pressure gas viscosity data in the literature are discussed. The application of the correlation to predict absolute viscosities is discussed in Part 11. Auxiliary correlations are presented to enable calculations of viscosities from a knowledge of the pressure, temperature, and gravity of the gas phase. INTRODUCTION A knowledge of the viscosity of hydrocarbon fluids is needed to study the dynamical or flow behavior of these mixtures through pipes, porous media, or more generally wherever transport of momentum occurs in fluid motion. Since flow is predominantly in the laminar region in petroleum reservoirs, the influence of fluid viscosity on this flow is especially important. As early as 1894, Onnes1 and Onnes and de Haas' noted that the viscosities of homologs under corresponding states could be correlated. The theorem of corresponding states has been further developed and applied to the viscosity of pure, nonpolar gases under pressure by Comings, Mayland, and Egly.3 It was demonstrated that the viscosity ratio could be expressed rather closely as a function of reduced pressures and temperatures, 1.e.: µ/µ1 = f( TR2 PR) ........(1) temperature, absolute units where: TR =------------------------------ ------ critical temperature. absolute units _pressure, absolute units critical pressure, absolute units µ = viscosity of gas at reduced temperature TR and reduced pressure PR µ1 =. viscosity of gas at atmospheric pressure and at temperature TR Serious discrepancies in the viscosity of pure hydrocarbon gases at high pressures have been called to Our attention by Comings, and and Egly. They made a careful analysis of the following methods commonly used to measure gas viscosities: 1. Oscillating disc viscometer4,5 2. Rolling ball viscometer6,7 3. Capillary tube viscometer12 On the basis of their analysis of the problem and their experimental data on the viscosities of methane, ethylene, carbon dioxide, and propane, they preferentially selected data obtained by the capillary tube viscometer to develop their viscosity correlation. Since then a need has existed to verify the experimental technique of Comings, Mayland, and Egly for the determination of viscosities at high pressure; to extend the data to mixtures of hydrocarbons; and to extend the range of their correlation to include higher pressures. Using the correlation procedure of Comings, May-land, and Egly, the prediction of gas viscosities is resolved into the correlation of the effect of pressure, temperature, and composition on the viscosity of hydrocarbons and their mixtures and the prediction of viscosities of mixtures of hydrocarbons and other natural gas components at one atmosphere pressure.

Jan 1, 1955

By B. E. Eakin, A. L. Lee, M. H. Gonzalez

Expeximental viscosity and density data of lour natural gases are presented for temperatures from WO to 340F and pressures from 100 to 8.000 psia. A correlation is also discussed and results reported. INTRODUCTION This investigation is one of several efforts by the authors to provide viscosity data for pure hydrocarbons and mixtures. Results of some pure hydrocarbons and their mix tures have been presented." Several predictive methods and correlations have also been reported.'," This paper presents the experimental viscosity and density data of four natural gases and confirms a correlation in a previous study. APPARATUS The viscometer was described previously. 2,6 A bank of stainless steel pycnometers was included to determine density in conjunction with the viscosity measurements (Fig. I). DATA AND MATERIALS The natural gases were furnished by the Atlantic Richfield Co. (Samples 1 and 21, the Continental Oil Co. (Sample 3) and the Pan American Petroleum Corp. (Sample 4). Table 1 shows the composition of these gases obtained by mass spectrometer analysis. The pure component viscosity data used in the correlation have been published."".'"', '" Experimental and calculated viscosity data for the natural gases are presented in Tables 2 through 5. CORRELATION Starling and Ellington"' have reported several semi-empirical expressions based to a certain extent on the theory of viscosity advanced by Born and Green.' The final expression presented by Starling and Ellingto'n is: i'. (micropoise) -- /t,, exp rX(7') p ..... (I) Eq. 1 was modified by Lee et 01." to represent mixture and pure component data simultaneously. This equation has the form: F = K exp [XpY] :.........(2) where K=- (7.77 i 0.0063M)Ti 122.4 + 12.9M +T .................(3) ' ' ' X= 2.57 + ---1914.5+ 0.0095M .... (4) Y = 1.11 -r 0.04 X ........(5) Over the pressure and temperature range studied in this investigation, this equation represents the data on methane, ethane, propane, n-butane and four methane-n-butane mixtures with a standard deviation* of ? 1.89 per cent. Density data by Sage and Lacey" were used to fit the equations. Eq. 2 gives reliable values of the viscosity of light hydrocarbons without prior knowledge of experimental viscosity values, but accurate density data must be available. Eq. 2 was used to predict viscosity values for the natural gases studied in this paper. The particular set of parameters contained in Eqs. 3 through 5 gave values that reproduced experimental data within 25 per cent. Needless to say, in using Eq. 2 better density values will give better viscosity results. However, those engineers who do not have accurate density values at hand or the time, or such facilities as a library, computer, etc. readily available, may feel uneasy in using Eq. 2. Therefore, the authors sought the easiest but by no means best density prediction method, which was reported by Kaf 30 years ago. A set of density values, which were calculated based on Kay's method, and a generalized compressibility factor chart were used to fit Eq. 2. The results are:

Jan 1, 1967

By D. B. Robinson, J. Lukacs

A knowledge of the equilibrium water content of hydrocarbon systems under pressure is important to the natural gas industry. The information available on the solubility of water in hydrocarbon, hydrogen sulfide, and carbon dioxide systems is reviewed in this paper and the influence of the more important variables such as temperature, pressure and molecular structure on solubility in liquids and gases is discussed. A suitable cbromatographic technique has been developed for determining low concentrations of water. Tailing of the water peaks has been elim inated by adding water to the carrier gas stream The method is applicable for both gas or liquid samples and is effective in the presence of hydrogen sulfide. The experimental study of water solubility in methane-hydrogen sulfide systems at a temperature of 16oF has shown that the presence of hydrogen sulfide causes only a modest increase in water content at pressures up to 1,400 psia. Theoretical considerations and data on pure hydrogen sulfide and carbon dioxide suggest that the effect of both these compounds will be greater at higher pressures and in the liquid phase. INTRODUCTION Before transporting or processing natural gases and gas condensates, it is usually necessary to dry them using suitable dehydration equipment. The design and operation of this equipment requires a knowledge of the amount of water present in the fluid at the reservoir and operating conditions. This is influenced by temperature, pressure and composition, particular1y when certain non-hydrocarbon components are present. Field experience indicates that hydrogen sulfide and carbon dioxide, for example, alter the usual water solubility relationships appreciably. However, an extensive search of the literature does not reveal any quantitative data on such systems. For sweet natural gases, generalized empirical correlations such as the one proposed by Katz, et al,6 can be used to predict water solubility with confidence at most temperatures and pressures of interest. However, existing theoretical relationships do not permit a calculation of the deviation from these curves when polar substances like hydrogen sulfide are present in the system. Thus one must resort to an experimental approach to obtain the necessary information. The fact that laboratory experimental methods frequently involve the use of mercury which reacts with hydrogen sulfide in the presence of water, and that hydrogen sulfide interferes with many chemical reactions specific for water has contributed to the difficulty of studying water solubility in systems containing hydrogen sulfide. In this investigation the water content of a limited member of methane-hydrogen sulfide mixtures was determined using a special technique with gas chromatography. REVIEW OF PUBLISHED DATA Experimentally determined water solubility data have been reported for methane, 14 ethane, 15 propane, 9 n-butane, 16 1-butene, 11 hydrogen sulfide 18 and carbon dioxide. 19 These studies report the effect of pressure, temperature and molecular structure on water solubility in single component gases and liquids. Fig. 1 compares experimentally determined water solubilities in several hydrocarbons, carbon dioxide and hydrogen sulfide at l00F as taken from this

By G. Pizzi, G. M. Ciucci, G. L. Chierici

The use of the material balance equation to estimate the volume of hydrocarbons originally present in a reservoir, whose producing mechanism is partly due to water drive, has been discussed in the literature by several authors. There is no general agreement upon the possibility of obtaining reliable results by this method. Gas reservoirs in contact with an active aquifer are considered in this paper. Theoretical considerations based on the cybernetic principle of uncertainty (which states that the internal structure of a system cannot be uniquely determined from its observed external behavior) lead to the conclusion that the volume of gas originally present in a reservoir of this type cannot be uniquely determined from its past history. The range of values which encompasses the actual value of the reserves varies from case to case and must be determined by either numerical or analogical methods. Results obtained for six gas fields are reported. All these fields were produced with small fluctuations in their production rates, as is common practice for gas reservoirs; no gas storage fields were considered. Results obtained show that reserves values in a range of 1 to 2, associated with appropriate aquifers, allow the matching of the reservoir past history with mean-square deviations less than the experimental errors involved in pressure and produclion measurements. Similar results have been found in several other partial water drive gas reservoirs. From these results it is concluded that gas reserves cannot be uniquely determined from the past performance of partial water drive reservoirs, at least in cases where the reservoir has been submitted to small fluctuations in the production rates, and pressure data of normal accuracy are available. INTRODUCTION A number of authors have analyzed the problem of estimating the reserves originally present in a partial water drive reservoir from its past pressure-production performance. Literature which deals with this subject can be grouped, according to their conclusions, as follows. 1. A reliable value for the reserves can be obtained even if reservoir data (pressure and cumulative production) are affected by errors within the normal range.'-' 2. A reliable value for the reserves can be obtained only if reservoir data are very accurate'- or if past production performance has been subjected to abrupt variations in the production rate."' 3. No unique value for the reserves can be obtained from reservoir past production performance. This conclusion has been based upon theoretical considerations" and verified in several field cases. The purpose of this paper, which deals only with partial water drive gas reservoirs, is to test the above conclusions against actual field cases. Some theoretical considerations on this problem are also presented. THEORETICAL CONSIDERATIONS As the behavior of a gas reservoir communicating with an aquifer depends on both the aquifer and the reservoir characteristics, the physical system to be studied is the combined gas reservoir plus aquifer. The information which is available for studying the performance of such a system is the well production rates and bottom-hole pressures, all given as functions of time. The external behavior of the reservoir-aquifer system is therefore described by 2n input variables (Gp, Wp,) and n output variables (p,), n being the number of wells in the reservoir. In reservoir engineering it is common practice to consider the reservoir as a whole, disregarding the internal distribution of pressures and of producing wells. This practice is equivalent to substituting the above multi-variable system with a single-variable system, where the average reservoir pressure is the only output variable and the cumulative production G,(t) and W,(t) are the input variables. The internal structure of such a system, defined by initial gas reserves G, aquifer shape and dimensions, boundary conditions and petrophysical parameters distribution throughout the aquifer, is unknown. Therefore, from a cybernetic point of view the system is a "black-box It has been demonstrated" that for a blackbox the indetermination principle holds. Accordingly, the number of different internal structures (or set of parameters) which can account for the observed external behavior is infinite. As a consequence, the initial reserves cannot be uniquely determined from the reservoir past performance. When W,(t) is known, the determination of G from

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

The trend in gas well testing has been to rely more on the early-time flow data of drawdown and buildup tests than on stabilized flow tests. The stabilized testing methods often are not adequate for complicated modern gas reservoirs because many of these reservoirs have extremely low permeabilities and the transient flow period of a well test can last for days or weeks. It has become common to augment the flow capacity of gas wells by hydraulically fracturing their producing formations. In deep reservoirs the induced fractures are generally vertical and tend to follow a single plane of weakness. The presence of a vertical fracture at the wellbore complicates the transient flow behavior of a low permeability gas well. The flow is further complicated when turbulence occurs near the wellbore. Russell and Truittl published transient drawdown solutions for vertically-fractured liquid wells. They developed methods of drawdown and buildup testing utilizing these solutions, which were based on numerical simulation. Clark2 applied the basic Russell-Truitt solutions to analyze fractured water injection wells by falloff tests. Field examples were given to substantiate the method of analysis. The Russell-Truitt solutions and well testing methods can be extended to gas well flow. Millheim and Cichowicz3 presented the drawdown equations for *Now with Scientific Software Corp., Denver. ideal gas flow, including the effects of formation damage and turbulence. Two actual field cases were presented, both of which exhibited extensive fractures and turbulent flow. For the pressure range of these well tests, the use of ideal gas equations was acceptable. These field cases were significant in being the first published data showing the occurrence of turbulence in fractured gas wells. We, also, observed such cases in practice. Our purpose here is to extend the theory of fractured gas well testing to the flow of real gases. To determine the effects of wellbore storage and turbulence on well test interpretation, we developed a finite-difference model to simulate well test conditions. Since gas wells usually are widely spaced and have high compressibility, the emphasis has been put on the early transient behavior before the effects of the outer boundary are noticeable at the wellbore. The Mathematical Model The geometry of our mathematical model is similar to that used by Prats.4 The well is centered in a circular uniform reservoir. A vertical fracture of infinite flow capacity penetrates the formation and passes through the wellbore. The wellbore itself is not important. Fig. 1 shows a sketch of the problem. Fig. 2 shows the idealized model. The equation for the flow of real gas in a porous medium with uniform porosity, permeability and thickness is5-7

Jan 1, 1970

By Jorge Gibson, P. Chris Pistorius, Megha Jampani

Increased utilization of natural gas in blast furnace ironmaking can decrease both the cost and the carbon intensity of ironmaking, given current US natural gas prices. In this paper, three ways to utilize natural gas are compared: tuyère injection, prereduction of iron ore, and shaft injection. The basis for comparison includes coke replacement ratios, carbon intensity and furnace productivity. These were calculated using relevant mass and energy balances, a blast furnace productivity correlation based on the bosh gas flow rate, and measured and modeled prereduction kinetics. Of the natural gas utilization methods, prereduction has the highest effective coke replacement ratio (and hence the largest advantage in raw material cost), but is likely to have the highest capital requirement.

Mar 1, 2017

By K. Adham

For use in many existing and emerging applications, the purification of natural graphite is required to achieve +99.9% carbon content with minimum metallic impurities. Currently, the hydrofluoric acid process is used, which carries certain environmental and workplace health and safety impacts. While thermal purification of graphite at temperatures over 2500 °C is a known alternative, such high-temperature furnaces are expensive to build and operate. Using chlorine at lower temperatures to purify graphite is also known and patented, but is currently limited to treating solid synthetic graphite shapes in small-scale batch furnaces. Chlorine treatment of natural flake graphite resources can be used commercially, if certain drawbacks are addressed through process improvements that can be helped with the use of fluidized bed reactor technology.

By S. S. Mauney

Natural and synthetic colors are used in almost every facet of daily life as a means of communication and expression. Natural ochres, umbers, hematites and magnetites were once used in their crude or semiprocessed state for their pigment value in construction products and low-end coatings applications. In recent years, though, they have yielded to synthetically produced iron oxides for more sophisticated products. Today, through innovative mining processing and marketing techniques, natural iron oxide pigments are once again gaining strong inroads into products and markets requiring precise chemical and physical characteristics. With the focus on environmental consciousness and resource conservation, value added natural iron oxides provide an economical and responsible alternative to color needs.

Jan 1, 1995

By S. S. Mauney

Natural and synthetic colors are used in almost every facet of our daily lives as a means of communication and expression. Natural ochres, umbers, hematites, and magnetites were once used in their crude or semi-processed state for their pigment value in construction products and low-end coatings applications, while yielding, in recent years, to synthetically produced iron oxides for more sophisticated products. Today, through innovative mining, processing, and marketing techniques, natural iron oxide pigments are once again gaining strong inroads into products and markets requiring precise chemical and physical characteristics. With our focus on environmental consciousness and resource conservation, value added natural iron oxides provide an economical and responsible alternative to our color needs.

Jan 1, 1993