Search Documents

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 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 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 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 M. Shahsavar, J. Sattarvand, J. Gomez

Mining can be categorized as a hazardous activity, considering some factors such as environmental conditions with a considerable presence of humidity, suspended particles, or falling rocks have affected the severity and number of accidents compared with other economic sectors. The industry analyzes incident reports to narrow the rate of severe injuries and fatalities, conducting root cause analysis and identifying leading indicators. As the International Council on Mining and Metals noted, the vast trove of incident data is not analyzed as much as possible due to a lack of analytics expertise at mine sites. However, machine learning could solve the problem of analyzing all the incident data in a no-timeconsuming way, considering the abundant data, and without using expert personnel in data science. Thus, a Convolutional Neural Network and a Naïve Bayes model were introduced to perform classification in the MSHA database. The database from 2020 consists of 60 fields to describe safety incidents; these fields include mine I.D., accident date, subunit (mill, surface), material extracted, and other metadata.

Jun 25, 2023

By Parke A. Dickey

POTENTIAL differences between strata of shale and sandstone have been recognized for about 15 years, and they form the basis of the electrical logging of oil wells. Hitherto these potentials have been ascribed solely to electrochemical reactions caused by the disturbing effects of the water in the drill hole through which they are generally measured. Observations by the author in wells empty of water and in a mine shaft suggest that potential differences between sandstones and shales are natural and specific properties of the rocks, although they are modified to measurement by the disturbing effects of the water in the drill hole. The cause of these potentials probably is to be sought in the relative polar adsorptive capacities of quartz and clay particles.

Jan 1, 1944

By C Noble, J Webster-Brown

A survey of selected West Coast Streams was undertaken in Summer 2001, Summer 2002 and Winter 2002 to establish the degree of arsenic (As) contamination in waters and sediments downstream of old coal and gold mine workings on the West Coast. Of the catchments studied, Soldiers/Ford Creek and Rapid Stream showed the most significant elevated As concentrations, predominantly in dissolved form, and evidence of AMD or ARD contamination (high sulfate concentrations and/or acidity). The rate and processes of removal of dissolved arsenic from the water column (attenuation) have been studied in these catchments. The Rapid Stream and Soldiers/Ford Creek catchments are directly affected by adit drainage from the Sullivan and Blackball coal mines respectively. Both catchments have high As concentrations in the sediments, with up to 770 mg/kg in Soldiers Creek and 490 mg/kg in Rapid Stream. Dissolved As concentrations are significantly elevated in the neutral pH waters of Soldiers Creek; typically as high as 100 - 200 ¦g/kg, up to 500 m downstream of the main adit discharge (800 ¦g/kg As). Adit inflows to Rapid Stream were acidic and contained 20 - 30 ¦g/kg dissolved As, with relatively high concentrations of arsenic in the stream water in winter (up to 18 ¦g/kg), but not in summer. Analysis of concentration trends and As:SO4 ratios downstream of these obvious As sources confirmed that dissolved As had been actively removed from the water column, but varied seasonally with greater removal evident in summer. Adsorption of arsenic onto iron oxide in the stream sediments is the most likely mechanism of removal. The Snowy River catchment, on the other hand, has relatively high As concentrations in both the sediments and the water column, but there was little evidence of a distinct AMD or ARD input from the Blackwater gold mining area. Preliminary investigation suggests that dilution rather than attenuation is responsible for decreasing As concentrations downstream.

Jan 1, 2003

By M. Kahn

This paper traces the rise of South Africa's mining industry in the context of the different state forms that oversaw its development. It commences with the mining-led industrial revolution through four periods of conflict, each of which contributed to a deepening of the industrial base, with constant battles between domestic and foreign capital and labour, and the rise of both Afrikaner and African nationalism. The role of state-owned enterprises has been crucial in shaping this dynamic so that at its peak the public sector and leading mining houses constituted two massive centres of power. This was followed by a thirty-year period of neo-liberal experimentation. The current debate on the 'developmental state' is an attempt to reconstruct the state as a counterweight to capital. It is argued that while the country has moved toward inclusive political and economic institutions, these are both limited by deep 'extractive' characteristics. This transformation failure provides a fertile ground for populism that presents one of the biggest threats to inclusive growth. Cited examples point to the many ways that states have sought to deal with the bonanza that conversion of mineral non-resources into resources brings. Such greenfield resources on state land, or in areas where the state holds monopoly rights, do not present ownership issues. It is suggested that this is the route that South Africa should follow, as the alternative will be conflictual.

Jan 1, 2013

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