Search Documents

By A. Satter, J. M. Campbell

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

By K. H. Kilgren

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

Jan 1, 1967

By J. W. Wolfe

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

Jan 1, 1967

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

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

Jan 1, 1952

By L. Cichowicz, K. K. Millheim

The constant-rate drawdown test performance for a low-permeability, verticany fractured gas well was investigated. A series of gar wells were tested by flowing each well at constant rate until the data could be analyzed using con-ventional radial flow theory. Each well was then shut in to build up. After a sufficient buildup was obtained, another flow Iest commenced but at a higher flow rate than Ihe first test. Again, the well was shut in when radial fiow was obtained. his procedure was repeated for three to four different flow rates. Two wells in the San Juan basin were tested using this procedure. Both wellere fractured after completion, cleaned up and then shut in until flow testing commenced. Test designs of both wells permitted investigation of the most realistic values of egective permeability, wellbore radius and turbulence factor. Also, being able to determine the eflective fracture flow area and vertical fraclure efficiency was inherent with this testing approach. It was observed that fractures in both wells influenced the Pressure behavior for approximately Is to 40 hours (depending on the flow rate) before radial flow was evident. After this time, drawdown data were analyzed using radial pow theory. When a low-permeability gas well has vertitally oriented, induced fractures, the early flow geometry ic. essentially linear. It will be shown how to determine when a flow test has been conducted long enough so lha' the most representative values of effective permeahiiity. wellbore radius and turbulence factw can be calculated. From the linear pressure data, valuable information about The fracture treatment, such as the effective flow area and vertical fracture efficiency, can be determined for vertically froctured wells. Introduction During tests on gas wells in the Sari Juan basin7 initial transient behavior lasted for many days because of the low permeability of some porous media. As a result, stabilized flow performance could not be obtained. If these wells received some type of stimulation treatment, early pressure behavior deviated from conventional theoretical radial Row. When conventional radial flow theory was used to analyze these low-permeability fractured gas wells, larger values of flow capacity and absolute open flow potentials (AOF) sometimes resulted. Wells were assigned open flow poten- tials that proved to be 3 to 10 times higher than the well would sustain over a longer period of production. In some cases where the wells had flowed for longer periods of time during a constant-rate drawdown test, it was noticed that the effective flow capacity appeared to be decreasing with time until a certain value was reached. The early nonradial pressure behavior can be explained If linear flow is assumed. Rusell and Truitt mathematic. ally investigated the vertically fractured well in a bounded area. They showed that early flow behavior was essentially linear and, for x,/x. approximately less than 0.10 radial flow was obtained after short periods of time. Then realistic values of effective permeability and skin could be determined. Scotta experimentaliy studied the vertically fractured well with a heat flow analog. He showed that earb flow was linear, Both studies indicate that, for small values of x,/x,, linear flow approaches radial flow if the well is tested long enough, To help prove this concept of early linear flow caused by induced vertical fractures, two low-permeability gas wells were tested. Both wells received large fracture treatments prior to testing. A vertical fracture was indicated from the analysis of fracture treatments, As anticipated, tests of both wells indicated early linear flow that was later followed by a period of radial flow, Data collected from each well were analyzed. From the well tests, plus other information on each well, the effective permeability, wellbore radius and turbulence factor were calculated. Effective fracture flow areas calculated from test analyses for each well proved to be approximately one-fourth the created area calculated from classic hydraulic fracturing theory: other fractured wells that were tested but not presented in this paper also indicated that the effective fracture flow area was one-fourth to one-third the created area predicted from hydraulic fracturing theory. The vertical fracturing efficiency was estimated from the calculated values of effective wellbore radius and fracture flow area. For the two wells tested, calculated fracture lengths x, were 112 and 105 ft, and the vertical fracturing efficiencies E, were 122 and 183 percent. Development of Flow Model Agnew' showed that most induced fractures below 1,500 ft are vertical. Anderson and Stah15 indicated that most of the fractures they studied were vertical. Thc model proposed for early flow in most vertically fractured gas wells is shown by Fig. 1. This model should approximate early flow behavior until radial flow is reached, at which time a radial model with an effective wellbore radius of 0.5 they will

Jan 1, 1969

By R. F. Hinds

A pressurizing system was designed and built to apply a radial pressure of 5.000 psi to rock samples. Samples of the Bradford, Weir and Kirkwood sandstones were subjected to radial pressures parallel and perpendicular to the bedding in the system and the changes in conductivity, porosity and permeability were deter-mined. Asymptotic decreases in conductivity, porosity and Klinkenberg permeability were observed over the 0- to 3,500-psi range. From these data, the increase of formation factor over the same range was calculated. The cementation exponent expressed in the Archie equation was found to increase with pressu,re. The rate of change varied from formation to formation. Only small differences were observed for samples taken parallel to the bedding as comnared to those taken perpendicular to it except for the Bradford sand. Above 3,500 psi the data are not good enough to tell if the properties continue to fall-off asymptotically or in a different manner. A possible explanation for the probable litnit of the asymptotic rate around 3,500 to 4,500 psi is that this represents the maximum pressure which the formations have unrdergone during their geological life as a result of burial. INTRODUCTION For many years measurements have been made in the laboratory of the properties of oil-bearing sandstones with little consideration being given to whether these properties had been changed when the samples were moved from their original position in the ground to the laboratory. Relatively little has been reported in the literature on the effect of pressure on the properties of sandstones. The changes in permeability for eight typical consolidated oil-bearing sandstone samples with overburden pressure (0- to 15,000-psi range) was presented by Fatt and Davis' in graphical form. In 1953 Fatt' also reported the relationship between porosity and overburden pressure for four typical consolidated oil-bearing sandstone samples (0- to 5,000-psi range). Fatt later described the changes of formation factor and porosity of 20 typical oil-bearing sandstone samples with overburden pressure changes (0- to 5,000-psi range) giving the values of the coefficients and exponents in the Archie' equation, and the more general Winsauer equation', where C, k and m are empirical constants, F is the formation factor, and + is the porosity. During the experiments of this investigation the core was contained in a neoprene rubber sleeve as shown in Fig. 1. It was held between two steel end members which were mounted in a steel cylinder. The unassembled pressure cylinder with its component parts is shown in Fig. 2. The core was insulated from the two steel end members by a Teflon sleeve with a small hole in the

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

By W. M. Rust, W. D. Mounce

THE almost universal acceptance of electrical logging by the petroleum industry calls for a critical examination of the physical bases of the common methods. This is particularly needed for the natural potential log. The conventional explanations are based on the phenomena of electrofiltration and electro-osmosis. However, experimental data show that they fail to give a basis that is quantitatively sound. Moreover, they lead to interpretations of details of electrical logs that are erroneous. It was found experimentally that when a plug of shale is placed in a chain of dissimilar electrolytes, a current flow results. Experiments on oil-well cores and fluids similar to those existing in wells show that potentials are set up which are in reasonable agreement with those encountered in electrical logging. Doubt- less other sources of potential are important but this explanation not only accounts for much of the potential but also clarifies some fairly familiar examples of misinterpretations.

Jan 1, 1944

By Gordon L. Draeseke

"Before getting into more specific aspects of my talk today I should put in one disclaimer - one week in Mainland China and one day with Premier Kosygin does not make me an expert on Sino-Soviet affairs nor on Canadian-Chinese relations. Perhaps though there is a message in the fact that one meeting took place in relative openness in China and the other was at least partly in a bullet-proof limousine in Canada. I thought it might be interesting and appropriate to explain something of the forest industry to you who are engaged in the mining industry. Our industry, like yours, is based on a natural resource and the basic activity that we engage in is the extraction and conversion of that resource into saleable materials. However, the fundamental difference between our industries is that our resource is on top of the ground - easily visible and easily measurable, whereas yours is under the ground - difficult to find and, having been found, expensive to measure. A second difference is that ours is renewable while yours is extractive, and you must keep finding more to keep up reserves."

Jan 1, 1971

By O. M. Wicken

Magnesite has the unique combination of properties that make it useful in a wide variety of metallurgical applications. The use of refractory magnesite has been increasing throughout the world, but the source of the "ore" differs. The United States has increased production of synthetic magnesites while the rest of the world is exploiting its resources, using natural magnesites. The history and methods of manufacturing synthetic magnesite are discussed as well as the practices for beneficiating natural magnesite. The properties, costs, and uses of the two types of magnesite are discussed and compared.

Jan 1, 1964