Search Documents

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

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

Jan 1, 1958

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

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

By 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 A. B. Cook, C. J. Walker, G. B. Spencer

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

Jan 1, 1970

By 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 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 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 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 L. J. E. Hofer

THIS BULLETIN analyzes and summarizes research at the Bureau of Mines dealing with the preparation, properties, and reactions of iron carbides in the period 1948-60, part of which was sponsored by Wright Air Development Laboratories. The objective was to develop better catalysts for the synthesis of methane and higher hydrocarbons from hydrogen and carbon monoxide. The investigations resulted in the first adequate characterization of epsilon iron carbide, which later proved to be an important intermediate in the tempering of martensite in steel. At temperatures near 170° C, the carbon in carbon monoxide reacts nearly quantitatively with finely divided alpha iron to form an iron carbide in which the metal atoms occupy the positions of the close packed hexagonal arrangement with lattice parameters of ao=2.754 A and co=4.349 A. The carbide has a composition within 2 atomic percent of that required for the formula Fe2C, a Curie point of about 370° C, a specific magnetization of 160.4 ergs per gauss-gram Fe at 50° C, and a ferromagnetic moment of 1.70 Bohr magnetons. Epsilon iron carbide can be converted to chi iron carbide by thermal treatment at 343° C for about 70 hours. The chi iron carbide so formed has a Curie point of 247° C, a specific magnetization of 142.5 ergs per gauss-gram Fe at 50° C, a ferromagnetic moment of 1.72 Bohr magnetons per iron atom, and a paramagnetic moment of 5.5 Bohr rnagnetons per iron atom. Direct synthesis of this carbide from carbon monoxide at 240° C indicates a composition within 7 atomic percent of that required for Fe2C. Chi iron carbide can be converted to cementite at 550° C in about 315 minutes. The theta iron carbide so prepared has a Curie point of 210° C, a specific magnetization of 147.3 ergs per gauss gram Fe at 50° C, a ferromagnetic moment of 1.72 Bohr magnetons per iron atom, and a paramagnetic moment of 3.9 Bohr magnetons per iron atom. Pure chi iron carbide reacts to form theta iron carbide by the elimination of free carbon. In the presence of an excess of free iron in intimate contact with the carbide, chi iron carbide reacts with the free iron to form theta iron carbide. The second reaction proceeds at an appreciable rate at 320° C and can be used to snythesize finely divided theta iron carbide with no free-carbon impurity at low temperatures. The ease of the reaction in the presence of free iron suggests that chi iron carbide is an intermediate in the tempering of epsilon iron carbide to theta iron carbide as it proceeds in steels. Experiments on carbon steel show that the magnetic moment of the iron atoms in martensite is about 2.2 Bohr magnetons per atom, a value approximately the same as that of the iron atoms in alpha iron. The Curie points of epsilon iron carbide, chi iron carbide, and theta iron carbide were found -in various specimens of tempered steel to be 10° to 20° C lower than normal values. This shift is probably the result of manganese in solid solution. Preparations of various compositions of epsilon iron carbonitrides show many effects of extensive isomorphism between epsilon iron carbide and epsilon iron nitride. The action of carbon monoxide on finely divided alpha iron at 400° C results in free carbon deposits containing filaments about 0.1 micron in diameter and several times as long. The filaments are characterized by nuclei of dense material located midway from each end of the filament. Although carbiding to near completion at 240° C over periods of about 300 hours produces a preparation containing no appreciable percentage of epsilon iron carbide, very brief periods of carbiding of the order of a minute result in the formation of detectable epsilon iron carbide even at temperatures as high as 325° C. The data indicate that epsilon iron carbide is an intermediate in the formation of chi iron carbide during its formation from carbon monoxide in the gas phase. Attempts to convert cementite to chi iron carbide by cold work increased the Curie point but did not change the X-ray powder diffraction pattern.

Jan 1, 1966

By Q. Liu, S. Liu, Y. Cui

To reduce the frequency of mining incidents and improve safety and production efficiency, it is important to develop technology for the autonomous positioning and navigation of mining equipment, without the assistance of satellite navigation signals. In this study we review the current technological environment for positioning and navigation underground and present several integrated navigation strategies. Existing underground positioning and navigation technology is described, and the performances of different methods based on this technology analysed and compared. The combination types and positioning effects of underground positioning and navigation strategies are discussed. Recent research directions for filtering methods and positioning estimation algorithms for underground positioning and navigation are reviewed. Finally, we summarize the future direction of underground positioning and navigation technology and related algorithms.

Jun 1, 2021

By Motohisa Kanda

The time and amplitude statistics necessary to adequately describe electromagnetic (EM) noise in mines are illustrated through computer software techniques. They are 1) Allan Variance Analyses (AVA), 2) Interpulse Spacing Distributions (ISD), 3) Pulse Duration Distributions (PDD), 4) Average Crossing Rates (ACR) , and 5) Amplitude Probability Distributions (APD) . These statistics are illustrated using data taken from a rather large store of raw analog data recorded in operational mines. The curves generated for the illustrations characterize the noise environment in the mines from which the corresponding data were taken, and should aid in the design of reliable communication systems for such mines.

By E. Lutz, M. Momayez, K. W. Brown Requist

As real-time air quality monitoring becomes more prevalent in US underground mines, it is important to provide the highest data reliability with the fewest possible sensors. Real-time sensors remain costly, and these costs are not exclusively financial; the time required to install, calibrate, and maintain real-time sensors poses a large barrier to widespread implementation. Current atmospheric monitoring systems typically rely on displaying point-specific values. This requires operators to infer real-time airborne contamination distributions. Monitoring and control software utilizing mine ventilation network (MVN) solvers has been implemented in limited cases because of their ability to simulate ventilation systems quickly, but these solvers use a one-dimensional representation of the mine, limiting spatial resolution of estimated distributions. Computational fluid dynamics (CFD) has likewise been considered as a means to improve spatial resolution, but processing times prevent its use as a basis for monitoring and control. For the real-time monitoring of airborne contamination distributions, we propose a spatial interpolation method that can estimate the distribution of airborne contaminants in near-real time. This method provides a middle ground between fast processing times and increased spatial resolution. With the use of a pathfinding algorithm and optimization through absolute percentage error minimization, this method outperforms spatial interpolation with a Euclidean distance. By providing contamination distribution information to operators, this method and its derivatives stand to outperform current atmospheric monitoring systems.

Apr 16, 2024

By John Heilig, Patrick Andrieux

"Near-field monitoring of blast induced seismic vibrations is a practical and useful way of assessing blast field behaviour for diagnosis and optimization purposes. It is also necessary in order to provide the information from which empirical vibration prediction equations, such as the one proposed by Holmberg (11 for example, can be derived. (Strictly speaking, the “near field” extends from the blasthole itself to about one charge length away from it; however, throughout this article, “near field”will.most often refer to the outer limit of this zone, from roughly 0.75 to 1.50 charge length away.) In order to obtain reliable data, proper field monitoring techniques are of paramount importance. Once the data are properly recorded, analysis can be done either in the time domain (describing vibration amplitude vs. time) or in the frequency domain (usually describing vibration amplitude vs. frequency)."

Jan 1, 1995

By Rajat Gangrade, Haotian Zheng, Jarrett Smith, Daniele Nebbia, Mike Mooney

The application of artificial intelligence (AI) to improve tunnel boring machine (TBM) operations requires a significant amount of information to train robust performance prediction models, e.g., for tunneling-induced settlement, productivity, tool wear, etc. This paper presents the AI framework applied to the North-East Boundary Tunnel (NEBT), the largest component of the Clean Rivers Project in Washington, DC. The paper focuses on required elements for real-time AI use to estimate tunneling-induced ground and building deformation and cutter tool wear. Keys to the framework include defining the necessary quantities of data required as well as the process of windowed training, validation and model testing. To date, the models developed for ground deformation and tool wear show satisfactory performance when compared to observed field observations during tunneling.

Jun 13, 2021

By J. M. Mutmansky

The usual method employed to determine the direction of flow in a connecting airway in advance is to conduct a survey between the two existing nodes being connected by means of the new airway. It is then assumed that airflow will take place from the node with higher pressure to the node with lower pressure. The purpose of this paper is to examine this direction-of-flow problem --using properties of the system characteristic curve for a two-terminal network --and specify some conditions that ensure that flow will take from the higher-pressure node to the lower-pressure node. It is concluded that the only requirement is that all characteristic curves of fans in the network should be strictly decreasing. Five numerical examples illustrate some important properties of the direction-of-flow problem.

Jan 1, 1999

By KeBler H-W, Seeliger A

Vibration analysis for diagnosis of machines is a powerful instrument for condition monitoring. Especially for diagnosis of antifriction bearings various equipment and techniques have been developed and are available. The envelope analysis is the most reliable method for bearing diagnosis. Different bearing failures lead to different patterns of the envelope spectrum. Checking the amplitudes of characteristic bearing frequencies is mostly not sufficient for a reliable diagnosis. Therefore human experts compare the pattern of the spectrum with typical patterns of bearing defects. In this paper a diagnosis system is presented which imitates this action of human experts. Using a neural net the system not only distinguishes typical bearing defects but also detects unbalance, beating of disengaged machine parts and other kinds of machine failures.

Jan 1, 1995

By D. H. Rubisov

The neural networks modelling approach is one of the most elegant black-box type methods to model industrial processes. In this paper, this approach is applied to model the performance of the CEZ zinc electrowinning plant. The model that was built was found to exhibit an excellent predictive ability for voltage and current density. The model was also used to investigate the sensitivity of the process to various key variables such as temperature, current density, zinc and sulphuric acid concentration, etc. Optimal regimes for zinc production were identified from the power consumption standpoint.

Jan 1, 1996

By C. Cho

One problem with the application of neural networks to the multivariable control of mineral and extractive processing is deciding whether and how best to use them. The objective was to compare neural network control to more conventional strategies and find any advantages in terms of set point tracking, rise time, settling time, disturbance rejection, and other criteria. The procedure was to develop neural network controllers using both historic plant data and simulation models. We tried various neural network control patterns including both inverse and direct neural network plant models. These were compared to state space controllers that were, by nature, linear. For grinding and leaching circuits, a nonlinear neural network-based model predictive control strategy was superior to a state space-based linear quadratic Gaussian controller. We also learned the importance of incorporating state space into neural networks by making them recurrent--feeding certain output state variables into input nodes in the neural network. We concluded that neural network controllers can have better disturbance rejection, set point tracking, rise time, settling time, lower set point overshoot, and can be more foolproof and easy to implement in complex, multivariable plants.

Jan 1, 1996