Search Documents

By P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By J. G. Richardson, K. H. Coats

During the initial growth of a gas bubble in an aquifer storage reservoir the injected gas tends to override the water. The resulting low displacement efficiency and high rate of gas travel down-structure make it difficult to maintain water-free production. This paper illustrates the use of two-dimensional, two-phase calculations to simulate this gas-water displacement. Calculations were performed for a stratified aquifer and for several homogeneous sands differing in permeability, porosity and thickness. Results are discussed and compared with performance predictions obtained from the simpler Buckley-Leuerett and Dietz formulas. This comparison indicates the necessity for the two-dimensional calculations in realistically simulating the gravity override of the water by injected gas. INTRODUCTION In recent years natural gas has been stored near markets in aquifers where insufficient storage capacity is available in depleted fields. Operators of these aquifer storage reservoirs have encountered technical problems relating to the gas-water displacement accompanying initial growth of the gas bubble. Injected gas tends to override the water, with a resultant low displacement efficiency and high rate of gas travel down-structure toward spill points. Low displacement efficiency makes if difficult to sustain water-free gas production. Sometimes the fingering of gas down-structure may be so pronounced that the injection rate must be severely curtailed, which excessively lengthens the time required for bubble growth. This paper is concerned with estimating — for given aquifer characteristics and fluid properties — the displacement efficiency, rate of gas movement down-structure and rate of gravity drainage of water behind the gas front. Several published articles relate to simulative capability necessary to handle this problem. The Dietz formula1 and Buckley-Leverett2 method have some applicability to the problem; Woods and comer3 reported one-dimensional, radial calculations which accounted for the two-phase flow of gas and water. Most recent research directed toward understanding aquifer behavior in relation to gas storage has concentrated on single-phase flow in aquifer.4 Douglas, Peaceman and Rachford5 presented a method for calculating multidimensional two-phase flow in reservoirs; their method has been applied in work reported by Nielsen and Tek,6 Blair and peaceman7 and Goddin.8 Refs. 5 and 7 show comparisons between experimental data and calculations similar to those employed here. One purpose of this paper is to illustrate the use of two-dimensional, two-phase flow calculations in simulating the difficult problem of water displacement by gas in a vertical cross-section. The pronounced fluid density and viscosity differences combine with the disparity in horizontal and vertical dimensions of the cross-section to pose one of the more difficult problems in reservoir simulation. The calculations described account for capillary and gravity forces, relative permeability and reservoir heterogeneity. An example reservoir is described and calculated results are presented for a variety of injection rates and values of permeability, reservoir thickness and dip angle. A second purpose is to compare the two-dimensional calculations with results obtained from the Buckley-Leverett and Dietz formulas. For this reason most calculations were performed for a simplified example reservoir to which the latter techniques might reasonably be applied. However, the computer program employed applies equally well to cases involving arbitrary spatial variations of permeability, porosity, reservoir thickness and dip angle.

By P. Pollard

A method has been developed for evaluating acid treatments in fractured limestone fields by breaking down pressure drawdown into three component parts: (I) pressure differcntial across "skin" near the borehole face, (2) pressure differential due to flow resistance in the coarse communicating fissures and (3) pressure differential between the fine voids and the coarse fissures. It is apparent that in most successful acid treatments the first term, skin resistance, has been reduced or eliminated. Further, it is often possible to estimate the volume of coarse fissures associated with the second term, coarse fissure flow resistance. In cases where this volume is comparable with practical acid volumes it seems likely that this resistance also may be attacked with a suitably retarded acid. INTROD UCTION Acid treatment has been successfully applied as a general practice in the limestone wells in the Mara/ Maracaibo districts of Venezuela for some 10 years. At the same time the need has been felt for a more precise method of: (1) evaluating the effect of an acid job and (2) selecting wells which should benefit from an acid treatment. The following method, based on an analysis of build-up curves for wells producing from fractured limestone, provides a means of forecasting the effect acid will have on a well and in many cases the probable order of production rate increase which may be expected. In certain cases it is also possible to calculate the fissure volume through which the acid has to be displaced for best results. Various approximations are introduced in the method, but within the limitations of these the analysis gives an indication of how a well, whether previously acidized or not, may be expected to respond to acid treatment. Suitable build-up curves cannot be obtained on "tight" wells which do not give sustained production and the method is not applicable to such wells. Figs. 1 through 7 illustrate various types of build-up curves obtained from limestone wells in Western Venezuela. Curves are included showing (1) both skin and coarse fissure resistance, (2) skin resistance only, (3 and 4) reduction of skin resistance by acidization, (5) fissure resistance only, (6) reduction of coarse fissure resistance by acidization, and (7) acidization failure. This theory has been successfully applied for over three years during which 22 acid treatments have been carried out specifically to reduce skin resistance and three treatments to reduce coarse fissure resistance. The method has been found to be a useful aid in the selection of wells for acidization.

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

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By 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 G. Thodos, R. B. Grieves

A method has been developed for the accurate calculation of the cricondentberm and cricondenbar temperatures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain any number of components including paraffins and isoparaffins, olefins, acetylenes, naphthenes and aromatics. This approach is based upon the mole fraction of the low-boiling component in the mixture and graphically presents the ratios of the cricondentherm and cricondenbar temperatures to the pseudocritical temperature as {unctions of a boiling-point parameter. The only requirements are the critical temperature, normal boiling point and the approximate vapor pressure behavior of each component. For mixtures of more than two constituents a stepwise calculation procedure is necessary where the pseudocritical temperature is based upon the critical temperature of the pure low-boiling component and upon the actual cricondentherm or cricondenbar temperature of the mixture of all the remaining higher-boiling components. From an analysis of 22 binary systems (123 compositionsl, the average deviation of calculated cricondentherm temperatures from reported values is 0.87 per cent, based on degrees Rankine; and for 10 multicomponent mixtures containing from three to six components, the average deviation is 0.98 per cent. From an analysis of 18 binary systems (104 compositions) the average deviation of calculated cricondenbar temperatures from reported values is 0. 79 per cent; and for 15 multicomponent mixtures, 1.47 per cent. Equations, derived from the graphical relationships, are presented which enable rapid calculations for both binary and multicomponent systems. INTRODUCTION With the increasing tendency towards the use of high pressure processes, the need for an accurate method of calculation of the properties of a multicomponent hydrocarbon mixture in the critical region is becoming more essential. Critical region phase relations are also significant in gas condensate reservoirs, for which a knowledge of the fluid phase boundaries and regions of retrograde condensation make it possible to evaluate optimum operating conditions. The possibility of liquefaction in the underground structure may be established, and information on the feasibility of the use of repressurizing may be determined. A typical phase diagram for a multicomponent hydrocarbon mixture is presented as Fig. 1. This figure also shows the regions of retrograde behavior encountered by following lines AB and AC. If accurate methods of calculating the temperatures and pressures at the critical, cricondentherm (maximum temperature) and cricondenbar (maximum pressure) points of multicomponent hydrocarbon mixtures of known composition were available, it would be possible to estimate closely the entire phase diagram of any mixture. Several methods have been presented in the literature for the estimation of the critical temperatures and pressures of multicomponent mixtures containing all types of hydrocarbons and including the fixed gases.4-8,21 Considering the cricondentherm and cricondenbar points, the work is more limited. Etter and Kay5 have developed rather complex equations for

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 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 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 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 Dennis P. Bryan

The southwestern United States has an abundance of natural lightweight aggregate (NLA) deposits of volcanic origin that are being used in the construction industry within the region. New Mexico, Arizona, Nevada and California utilize NLA in portland cement concrete and concrete masonry units for weight reduction and related economic reasons. The geological environment of the western United States is amenable to the occurrence of natural lightweight rocks because it contains Tertiary and Quaternary volcanic products that include low-density glass and vesicular material. There are two major classical types of natural lightweight materials of igneous origin; pumice and cinder. These account for the vast majority of volcanic rocks that are presently utilized as NLA. Pumice is an ejecta product generally found as part of tuffaceous units associated with nearby siliceous volcanic centers. Cinder, also an ejecta product, is associated with mafic volcanics and is most often found as a cinder cone. Both are highly vesicular. A third type of natural lightweight material, though also associated with siliceous volcanics, is less vesicular than pumice. It is best described as being a slightly pumiceous rhyolite or rhyolite-tuff breccia that occurs as a siliceous dome or tuft unit that is most likely to be glassy. Uses of these natural lightweight rocks center on the construction industry where they are utilized in lightweight portland cement concrete and lightweight concrete masonry units. The major reason for usage in portland cement concrete is cost saving effected by weight reduction in high-rise construction. The major usage in concrete masonry units, or "block", is for weight reduction that enables a mason to lay more lock.

Jan 1, 1987

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

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

Jan 1, 1944

By 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 A Black, J H. Youngson

At least 85 000 tonnes of mine tailings were discharged into Deepdell Creek, east Otago between 1890s and 1940s. The tailings contained about 1 wt % arsenic, so about 850 tonnes of arsenic in the form of arsenopyrite was discharged. Annual floods have removed nearly all tailings within c. 60 years. Small amounts of arsenic-bearing sediments (<250 ppm) remain in Deepdell Creek. About 7 kg/year dissolved arsenic was discharged in groundwater from mineralised Hyde-Macraes Shear Zone rocks in the vicinity of the historic mine sites. Adsorption of this arsenic on to iron oxyhydroxides in the stream sediments may have contributed to the elevated As content of those sediments. Even less arsenic remains in the main Shag valley below the Deepdell Creek confluence (15 - 25 ppm As). Rare relict arsenopyrite grains are trapped in bedrock crevices 25 km from source. No tailings deposition is detectable in theShag River estuary, and the tailings have presumably been discharged into the sea. Natural rehabilitation of the Shag catchment has been rapid and effective.

Jan 1, 2004

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 Winfield G. Wright, Philip L. Verplanck, Philip L. Hageman, Douglas B. Yager, Dana J. Bove, M. Alisa Mast

Geological, hydrological, and geochemical information synthesized in a Geographical Information System (GIS) for water and rock surrounding South and Middle Forks of Mineral Creek, northwest of Silverton Colorado, was analyzed to distinguish between the natural and the mining-related sources of metals to dace waters in the watershed An important natural source of metals to surface water emanates from a porphyry molybdenum deposit south and upslope from the Middle Fork of Mineral Creek. Interaction of surface and ground water with fractured, altered rocks and permeable, Quatemay-age surficial deposits produces downstream water quality that does not meet current State of Colorado water quality standards for metals and pH. GIS site characterization and three-dimensional modeling distinguish surficial deposits as volumetrically substantial, significant sources of naturally occurring metals that strongly influence water quality and complicate ongoing mine-site mitigation efforts. Our findings indicate that regional hydrologic and geologic information needs to be evaluated before realistic water quality guidelines are established and mitigation recommendations are confidently implemented.

Jan 1, 2000