Search Documents

By G. Nabi, W-K. Lu

Cylindrical specimens of natural dense hematite were reduced to magnetite at atmospheric pressure in H2-H2O mixtures of known composition over the temperature range 1084° to 1284°K. The rate of reduction was measured by the rate of movement of the interface between hematite and magnetite. The diffusion of gases through the gaseous boundary layer, the magnetite layer, and the interfacial chemical reaction were all considered in the interpretation of experimental data. The mass transfer coefficient through the boundary layer was calculated using accepted correlations. Values of the chemical reaction rate constant and the diffusivity of hydrogen in the magnetite phase were determined. THE present investigation is concerned with the reduction kinetics of natural hematite to magnetite by H2-H2O mixtures in the temperature range 1084" to 1284°K at atmospheric pressure. This reaction is the first step in the series of topochemical reactions in the process of reducing hematite to iron. Kinetic information of the simple steps such as hematite-magnetite transformation is necessary in order to have a better understanding of the complex processes of hematite reduction in iron-making. It also has direct industrial significance because magnetic roasting is one of the most important methods in benefication of lean ore.&apos; Although many technical papers have been published on the process of magnetic roasting and iron oxide reduction, very little information is available in the literature concerning the fundamental nature of hematite reduction to magnetite by reducing gases. Hansen et al.2 reduced the dense synthetic pellets of high-purity oxide in CO-CO2 mixtures and determined the reaction rate by weight-loss method. They were able to interpret most of their results by applying the interfacial area control theory developed by Mckewan.3 In contrast, Wilhelm and St. Pierre,4 who studied reduction of hematite to magnetite in H2-H2O mixtures by weight-loss method, stressed that the resistance of the porous magnetite layer to the diffusion of gases cannot be neglected in consideration of the overall reaction rate. In the present study the contributions of interfacial chemical reaction, diffusion of gases through the magnetite phase, and the gaseous boundary layer to the overall reaction rate will be considered. APPARATUS AND PROCEDURE Hematite Specimens Preparation. Natural hematite ore from Vermillon range of Northern Minnesota was selected for the present investigation because of its high purity and thermal stability. Chemical analysis of five samples gave the following average values: 67.52 pct total iron (96.62 pct Fe2O3, 0.28 pct FeO, 0.03 pct metallic iron), 2.53 pct SiO2, <0.07 pct MgO, 0.03 pct CaO, 0.05 pct combined mixture, 0.07 pct loss on ignition, and 0.34 pct other. Cylindrical specimens of 0.93 cm in diam and 2.7 cm in length were drilled from slabs of ore with a water-cooled diamond core drill. These specimens were heated to 1000°C and furnace-cooled. Specimens with silica pockets developed large cracks. The uncracked specimens were heated a second time, and their surfaces were carefully examined with a microscope. Those with hairline cracks or surface inhomoaenitv-- were rejected. Preparation of H2-H2O Mixtures. H2-H2O mixtures were prepared by the combustion of H2-O2, mixtures in a pyrex glass chamber in the presence of a catalyst. Alumina pellets coated with palladium, supplied by Englehard Industries, were used as the catalyst. Purified grades of hydrogen and oxygen were used which were repurified by usual techniques. Hydrogen before entering the combustion chamber was passed through an activated alumina H2O absorption bulb, with copper turning at the top. The cover of this bulb was not made pressure-tight so that any pressure development in the hydrogen line would cause the cover to blow off and also the copper turnings would act as a flame arrester in the case of a flashback from the combustion flame. Oxygen flow rates were measured with a bubble flow meter after purification with 1 pct accuracy. Hydrogen flow rates were measured by "precision wet test meter" and the amount of unburnt hydrogen was accurately measured by a bubble flow meter, after condensing water vapor in the gaseous stream. The Pyrex glass bulb contained concentric Vycor glass tubes as shown in Fig. 1. Oxygen was prevented from diffusing into the hydrogen line by threading platinum wire through pores at the combustion end of gas inlet tube. The glass bulb was heated with a Kanthal heating wire pasted in asbestos paper. The surface temperature of the bulb was measured with a thermocouple and adjusted to remain at approximately 350°C. The gaseous reaction chamber also served as a preheater for gases to avoid thermal segregation. The following sequence of operation was adopted. 1) Nitrogen was passed through the outer concentric tube to purge the catalyst bulb of oxygen. 2) Hydrogen was introduced through the inner tube until a steady flow was obtained. 3) Oxygen was then introduced into the nitrogen stream passing through the outer tube. 4) When combustion had commenced and a flame was visible over the platinum wire, the N2 was turned off.

Jan 1, 1969

By Harry Littleton, Jr. Scarber, Wanliang Sun

"The Lost Foam Casting (LFC) process is relatively immature when compared to other casting processes. To make the process more robust, computer modeling is needed to reduce lead times and casting defects. In this study, metal filling of lost foam aluminum castings was visualized using state-of-art real time X-Ray technology. Metal filling of the lost foam casting was simulated using a commercial software package. The result of the computer simulation was compared with the results of the real time X-Ray visualization. Some of the computer simulation results agreed with the real time X-Ray visualization, but some did not. Based on other studies using real time X-Ray technology, new input variables, such as foam density variations, were developed for the computer modeling. With the new variables, modeling of the lost foam casting process becomes more realistic. IntroductionMold filling of the Lost Foam Casting (LFC) process is an important step toward making a quality part. In the LFC process, molten metal is poured against a refractory coated foam pattern surrounded by dry sand. At the pouring temperature of aluminum, the foam pattern is decomposed into 30% (wt) gaseous and 70% (wt) liquid foam pyrolysis products. Existence of the gaseous foam decomposition products significantly increases the back pressure on the metal stream compared with the conventional open cavity sand mold. Entrapment of the liquid pyrolysis products during the metal filling generates many casting defects, such as internal pores, fold and blisters. Existence of these casting defects sometimes seriously degrades the quality of the components manufactured by the LFC process. Understanding and controlling the metal filling of the LFC process will certainly make this process more robust."

Jan 1, 2004

By H. Huang, R. N. Gow

"The reactivity of an enargite mineral electrode was investigated spectroelectrochemically under conditions typically found in alkaline sulfide leach solutions, with pH 11-13 and Na2S additions of 50 and 250 gL-1. A bulk enargite electrode was cycled between 0 and -1,000 mV for each set of solution conditions. Raman spectroscopy was used to determine and compare surface formations against species in mass-balanced EH-pH diagrams generated using STABCAL. Previously, it was determined that arsenic could be selectively leached under alkaline conditions at an average solution potential of -300mV, across the pH range 11-13. Operating under such conditions allowed for the selective leaching of arsenic with elemental sulfur formation. From several of the voltammograms obtained during this current testing, the formation of elemental sulfur appeared at increasingly reducing conditions with increasing sulfide concentration. At high concentrations, copper started to dissolve, possibly as a bisulfide complex, reducing the selectivity of the leach, unlike previous studies. However, the presence of the higher sulfide concentration is also anticipated to prevent passivation by consuming any surface formation of elemental sulfur. This study provides insight into the operating parameters of alkaline sulfide leaching in regards to the selectivity of arsenic leaching, as well as guidance for control of reagents and solution potential and, of course, the need for water treatment afterwards.IntroductionEnvironmental limitations on the pyrometallurgical processing of arsenic-bearing concentrates have forced operations to consider other solutions. Attempts to hydrometallurgically process enargite (Cu3AsS4) fall into two categories: the co-dissolution of arsenic and copper under acidic conditions, and the selective leaching of arsenic under basic conditions. The authors believe that the selective leaching of arsenic is more amenable because it produces a solid copper phase, preferably chalcocite (Cu2S) or covellite (CuS) that could then be processed using conventional methods already in place in the copper and gold industries. Potential leaching methods include alkaline sulfide leaching (Achimovicovica et al., 1999; Anderson and Twidwell, 2008a,b; Tongamp et al., 2009 and 2010), hypochlorite leaching (Vinals et al., 2003; Curreli et al., 2005; Mihaljovic et al., 2007), ammoniacal leaching (Gajam and Raghavan, 1983), and acid leaching (Gow et al., 2014b) with or without fine grinding to mechanically activate the enargite surface and increase normally slow leach reaction kinetics (Achimovicova et al., 1999; Balaz et al., 2000; Balaz and Achimovicova, 2006). In this study, the reactivity of enargite under alkaline sulfide leach conditions was examined."

Jan 1, 2015

By R. L. Robinson

An analytical framework is formulated to assist mining companies in assessing the relative merits of mineral exploration and acquisition strategies. The anticipated costs, risks, and returns of each strategy are conceptualised using decision trees. Economic criteria - expected net present value, expected return per unit of investment, and the probability of success - for each strategy are derived(from the appropriate decision tree network. Application of the analytical framework is illustrated using the example of a medium-size company in the Canadian base metal sector.

Jan 1, 1987

By Gregory D. Savage, Bruce M. Everett

Introduction The growing role of Colombia as a major supplier of internationally traded steam coal has been widely discussed in the industry. Although Colombian coal has been produced and exported in small quantities for a number of years, the country&apos;s position as a major coal exporter began in earnest in 1985 with the first commercial exports from the Cerrejon North Block Project. The mine is located on the Guajira peninsula in northeast Colombia, about 145 km (90 miles) from the coast (Fig. 1). In addition to the mine itself, the project involves construction of a dedicated railroad and port. This forms a modern, efficient, and integrated system. The project is embodied in a 50-50 association contract between International Colombia Resources (Intercor), a wholly-owned subsidiary of Exxon, and Carbocol, the Colombian state coal company. Intercor and Carbocol share investment and operating expenses equally. Intercor has the additional responsibility of managing operations on behalf of the partners. Each partner is entitled to 50% of the coal produced, and each markets its share separately. The project is a multi-seam, truck and shovel operation, based on large reserves (2 to 3 Gt or 2.2 to 3.3 billion st) of high quality bituminous coal. Exports began in the first quarter of 1985, at modest levels through what is called Early Coal Facilities. The permanent coal handling system began operation in February 1986. Production will increase rapidly in each of the next few years, reaching the design capacity of 15 Mt/a (16.5 million stpy) in 1989. Based on available resources, however, the resource will support expansion, and continued growth in production into the 1990s is feasible. Importance of transportation The success of a new coal export venture, such as Cerrejon, depends heavily on efficient transportation. In the international coal industry, overland and ocean transport can in some cases account for 40% to 50% of the final value of the product. Many coals, such as western US subbituminous coals, cannot reach export markets at all because of their low heat content and great distance from export markets. By the same token, low-grade resources are often commercially attractive for local markets even in the face of competition with much higher quality coals from distant sources. Figure 2 shows the growth in steam coal imports into Europe, the largest trading market. Through

Jan 6, 1986

By C B. McCampbell, R W. Jr Bickes

Sandia National Laboratories has developed a semiconductor bridge (SCB) igniter for the ignition of a variety of explosive materials. When subjected to a low-energy current pulse, the bridge bursts into a plasma causing rapid ignition of the explosive material pressed against the bridge. Despite the low energies required for ignition, SCB components have proven to be explosively safe, meeting both no-fire and electrostatic discharge (ESD) safety requirements.

Jan 1, 1991

By G. T. Brow

An experimental method is described for testing diffusion of gases in crude oils at pressures and temperatures up to 3,000 psia and 200° F. A sampling technique was developed to withdraw by constant-pressure mercury displacement about 60 µl of oil from the diffusion cell through any of several 1/2-in vertically spaced capillary tubes and to trap about 4 µl of the oil in a specially designed microsampler. Oil trapped in the microsampler can be completely inserted in the sample inlet of a chromatograph. From analysis in a three-stage chromatograph (GLC), gas concentration in a sample is determined from the ratio of gas peak heights to area of peaks of reference components (hexanes or heptanes) in the oil. A diffusion coefficient is calculated from measured concentration profiles with known equations for diffusion in finite or semi-infinite systems, depending on whether the gas has or has not reached the end of the oil column. Results are given for a test of methane diffusing at 2,560 psia and 160° F in crude oil from the Weber Sandstone reservoir in the Rangely field, Colorado. A constant diffusion coefficient of 3.34 X 10-5 cm2 per sec was obtained which gives a good fit between calculated and experimental profiles determined at six times after start of diffusion. The maximum standard error of gas concentration determinations experienced with this method was ±0.88 ml per ml of gas-free oil.

Jan 1, 1970

By Shaolong Chen

Feeding dry sulfide concentrate into smelting furnaces enhances process performance and reduces the energy consumption and exhaust gas flow. Direct heated rotary dryers and flash dryers have being replaced by rotary steam dryers featured by very low wear rate of the pipe heating elements and long maintenance interval. Impact of the drying processes on the environment receives more and more attention as the emission legislation tightens. In this regard, environmental factors such as SO2, CO2, dust and gas flow among rotary dryers, flash dryers and steam dryers were compared. Selection of a rotary steam dryer for copper and nickel concentrate drying would be an ideal solution to reduce the pollution of exhaust gas, as drying is achieved by indirect contact to the steam at relatively low drying temperature and as the off-gas contains almost only water vapor and air. The emissions can be significantly reduced and only small and simple off-gas handling units are needed.

Jan 1, 2006

By L. B. Phelps

This paper focuses on the use of an integrated mine planing software package entitled "REMINE" which was constructed to provide a project evaluation tool for analysis of remining projects. Coal remining constitutes the affecting of sites which have had a past history of mining and have a pre-existing pollutional discharge. Such sites require a Best Professional Judgement (BPJ) analysis. The program was written to provide such an analysis with flexibility in the mine plans which can be analyzed. The process required for a BPJ analysis will be detailed and mine plan simulation performed to demonstrate the analysis.

Jan 1, 1990

IRON AND STEEL COMMITTEE ALBERT SAUVEUR, Chairman. A. A. STEVENSON, Vice-Chairman. HERBERT M. BOYLSTON, Secretary, Abbot Bldg., Cambridge, Mass. John Birkinbine, William Kelly, J. S. Unger, William H. Blauvelt, Charles Kirchhoff, Felix A. Vogel, James Gayley, Richard Moldenke, Leonard Waldo, Henry D. Hibbard, Joseph W. Richards, William R. Walker, Henry M. Howe, C. F. W. Rys, William R. Webster, Robert W. Hunt, E. Gybbon Spilsbury, Frederick W. Wood. J. E. Johnson, Jr.,

Jan 1, 1914

By Chase E. Sutton

The porosity of a sand is usually the measure of its capacity for holding fluid, although porosity is sometimes selective; that is, a small percentage of the actual pore space (voids), particularly in fine-grained sandstones, may be isolated by the cementing material which prevents the entrance of fluids under normal conditions. Obviously, however, porosity determinations are essential to an engineering study of underground conditions of an oil field, and a simple and reliable method of obtaining a reasonably accurate figure for average porosity is desirable. This report describes a method suggested by Ruasell,3 which was used by the writer in the preparation of a report on the Texhoma-Gose pool, Archer County, Tex.4 The conditions in this pool were particularly favorable for obtaining an average porosity figure, in that some 30 to 35 core samples of the pay sand were available for analysis.

Jan 1, 1928

By S. Benedict Levin

In the spring of 1944, during the emergency In the supply of strategic mica, the Bureau core-drilled six aggregating 913 feet, at the Blister mica Trim. Cheshire County, N.H. These holes indicated a downward extension or the mica shoot which was being mined at that time and developed information essential for economic mining of the deposit.

Jan 1, 1949

By G. H. K. Pearse

Long range transport of pollutants, particularly of sulphur dioxide and the associated acidic precipitation, has received much attention in the press and in international discussion during the last few years. Studies revealing increasing acidity of lakes and the dying out of fish populations over broad regions of the eastern half of North America have created alarm. In this paper, I will not attempt a description of the acid rain phenomenon. For this r recommend the balanced treatment of the subject by Gene Likens et al (1979) in Scientific Americana For my purposes I accept that the problem is serious and direct my attention to exploring a solution. Ultimately, the key lies in the answer to the question; "What do we do with the enormous tonnages of pollution abatement byproducts, chiefly sulphur compounds, that would arise from pollution control?"&apos; Little thought has been given to this most limiting factor to environmental clean-up. This is symptomatic of a basic flaw in the [philosphy] that under-ins environmental policy. Although this paper is to propose a solution to the abatement byproduct problem, the basic flaw - the firmly entrenched, universally adopted, &apos;polluter pays [philosphy] as it is currently [understood] - needs some comment.

Jan 1, 1981

In addition to capital and operating costs, the accurate estimation of ore reserves is a key component to the successful outcome of any feasibility study. In recent years, there have been many project failures caused by the improper estimation of ore reserves. The accurate estimation of an ore reserve relies heavily on an understanding of the geology and the underlying assumptions unique to the estimation techniques being used. These assumptions are reviewed, and the shortcomings of the estimation techniques and their possible impact on a given project&apos;s economics are discussed. The latest developments in geostatistics in overcoming existing reserve-estimation problems are presented.

Jan 1, 1998

By E. C. Lane

During the past 15 years the Bureau of Mines has made and published several hundred analyses of typical petroleums and has developed a system of interpreting the analyses and of classifying crude oils. In describing crude petroleum frequent reference is made to the "base" of the oil, a usage the petroleum industry has found convenient. As part of its general study of the characteristics of crude oils, therefore, the Bureau of Mines has developed a system of classifying them according to "base." This system as first published in 19274 was limited to four classes and in its original form is no longer considered adequate. Since the system was first published the Bureau has analyzed a large number of crude oil-, including representative samples from nearly every important producing field in the world, and has found that the system of classification could be revised and amplified to advantage. The revised system is presented herein. Bulletin 2916 states that many oil producers and refiners use the word "base" much as it is used in pharmaceutical work. For example, the "base" or "basis" of an ointment or salve is the fatty material (lard, petroleum, or lanolin) used as a vehicle or carrier for the active medicinal principles. Apparently many of the pioneers in the petroleum industry considered crude oil similar in nature to these medicinal products in that it consists of a semisolid, nonvolatile base and 1 volatile portion can be distilled off to make gasoline, kerosene, and other products. Many crude oils when chilled slowly deposit a waxy material (paraffin); this fact may have been one reason for the use of the term "paraffin base" in referring to such oils. Almost all of the petroleum first discovered in the United States was of this nature; not until the oil industry was extended to California and along the Gulf coast was petroleum found that contained no wax but apparently did contain asphaltic material Accordingly, these oils were known as "asphalt-base" oils. The use of this term, however, seems unfortunate, as better-grade oils of this type contain little if any asphalt; the descriptive term "naphthene base" now seems to be a more desirable term than "asphalt base." The term "intermediate base" is also used commonly in referring to oils intermediate in oharacteristios between those of paraffin and naphthene base.

Jan 1, 1935

By Thang V. Nguyen, Lisheng Shao, Gary E. Taylor

"The Community Memorial Hospital of Ventura, California, proposed to build a new state-of-the-art Cancer Treatment Center with radiation shielding walls and roof. The loose silty soils at the site had the potential to liquefy during an earthquake, which would result in significant settlement under the heavy footing loads. Treatment options included various deep foundation alternatives, as well as several ground improvement methods. The choice of wet soil mixing minimized the disturbance to the adjacent hospital during construction and potential settlement. Wet soil mixing allowed the structure to be supported on a raft foundation, meeting the stringent building total and differential settlement criteria. The authors also performed unit cell 3-D finite element analyses to evaluate the behavior of the wet soil mixing system in seismic loading conditions. This paper highlights the ground improvement design, construction, QA/QC, and the research related to evaluating the effectiveness of soil mixing for liquefaction mitigation.INTRODUCTIONThe Ventura Community Memorial Hospital proposed the construction of a new Cancer Treatment Center to replace the existing structure, which was built in 1901. The hospital complex is located in a wellestablished residential neighborhood. The soil conditions at the site of the proposed structure required ground modification to support the very heavy radiation treatment section of the building and limit settlement. Originally, 90-foot-long precast piles were considered for foundation support. Based on supplemental Cone Penetration Tests (CPTs) performed, the deepest CPT did not encounter competent pile bearing strata at a depth exceeding 30.5 m (100 ft). This could result in failure of the piles under long-term conditions, especially if hydro-consolidation occurs. The more realistic pile penetration depth would have been 36.6 to 45.7 m (120 to 150 ft). In addition, the excessive noise in association with pile driving would have been unacceptable to the community, patients, and staff of the hospital, so an alternative ground improvement solution was required.The specialty geotechnical contractor recommended wet soil mixing, an economical solution of creating a soil-cement product known as soilcrete that would significantly reduce the construction related noise and provide a more complete deep foundation for the structure.SOIL CONDITIONS AND SETTLEMENT CONCERNSThe subsurface soils were interbedded (Fig. 1) consisting of loose to medium dense silty sands, and soft to medium stiff silty clays to a depth of 30.5 m (100 ft). Groundwater was encountered at depths of 9.1 m (30 ft) and historical high groundwater was at 8.5 m (28 ft) below the ground surface."

Jan 1, 2015

By Bulters J, Sheppard I

Two blasthole density logging programmes have been conducted on the Mt Whaleback iron ore mine to evaluate the use of this type of logging as an aid to blast pattern design. The first programme also involved the setting up of a Density Test facility with which to calibrate the density probe. The programmes indicated that about 22-23 blastholes can be consistently logged per day at a cost of $1.30 per metre. Even though this is only half the rate required to log the average number of blastholes drilled daily on Mt Whaleback, testwork indicated that no significant loss of accuracy results if only half the holes in a pattern are logged. The major results of the programmes were: 1. The density range for black pyritic Shale is 2.2-2,9 g/cc. 2. The density range for soft to medium hard BIF is 2.1-3.0 g/cc. 3. The density range for hard BIF is 3.0-3.6 g/cc. 4. The front row blastholes show consistently lower densities than the main body holes due to blast damage. 5. A low density zone extending to a maximum depth of 6 m occurs at the top of most blastholes and results from blast damage from the bench above. These density results have been used to implement savings in drilling and blasting via pattern expansion and charge optimisation of blastholes. Eight ore

Jan 1, 1986

By Edward Thomas

Note. - Securing mine roof through consolidating thin or otherwise incompetent strata to increase overall thickness of the immediate roof beam is a new method of roof support and control which has aroused great interest amongst mining engineers and which is now being experimented with and thoroughly investigated by the United States Bureau of Mines. A paper on Roof Bolting, illustrated by photographs and motion-picture films, will be presented by Edward Thomas, of the U. S. Bureau of Mines, at the forthcoming Annual Meeting of the Mining Society of Nova Scotia. In the following pages, Mr. Thomas describes the methods used in roof bolting as illustrated by the photographs that accompanied his paper. The photographs are supplied through the courtesy of the United States Bureau of Mines, whose willing cooperation with miners and mining engineers in Canada is very well known and even more greatly appreciated. This method of publishing Mr. Thomas&apos;s paper is adopted because few, if any, Canadian mining engineers have actually seen installations of roof bolts and it is felt that the best possible presentation of such an interesting subject through the medium of the Bulletin will be attained by publication of the photograps over explanatory captions. - Ed.

Jan 1, 1950

There are references to the early and continuous production of coal since 1851, but few definite figures. Some data were obtained from Bancroft, but nearly all the tonnages in Table 80, prior to 1869, are estimated. [ ]

Jan 1, 1942

By E. G. Woods, A. G. Comer

A mathematical model is presented and solved for determinution of the saturation distribution and pressure in a radial gas-storage reservoir. The model consists basically of two parts: (I) the growing gas-bubble core, and (2) the surrounding aquifer. Since the total pressure at the injection well is a function of the two-phase flow in the gas bubble and the unsteady single-phase flow in the aquifer, the resistance to flow in both zones was taken into consideration. The assumptions involved for both the radial equivalent of the Buckley-Leverett two-phase-flow equation and the injection-pressure equation are as follows: (I) the geometry is radial, (2) the gas bubble is free to expand or contract, (3) compression or expansion of the gas within the bubble may occur at the beginning of a time step, (4) the fluids are immiscible, (5) water is incompressible within the gas-storage region whereas it is compressible outside of this region, (6) a stable gas-water interface exists and (7) gas injection occurs at a constant rate or a series of constant rates. The mathematical model was solved numerically using an IBM 650 computer. A comparison is presented between the predicted results of the model, the results assuming steady-state flow and the actual initial injection-pressure history of an operating reservoir. Using the initial field pressure for a basis, the average deviation between the predicted pressure and the actual field pressure was less than 4.3 per cent. INTRODUCTION Underground storage of natural gas in abandoned oil fields, in abandoned coal mines, in caverns and in aquifers has had varying degrees of success. The lack of these first three facilities in the vicinity of most major gas marketing areas is leading more and more to storage in virgin aquifers. For the most part, this storage in aquifers has led to the study upon which this paper is based. Several excellent papers have been published discussing the problems involved in the underground storage of natural gas; however, to these writers&apos; knowledge, none of them considered the two-phase flow of fluids in the gas bubble and the fact that the gas bubble will grow or shrink in size depending upon the injection and withdrawal history of the reservoir. Not only will both of these factors influence the required injection pressure, but also they will have a definite bearing upon the amount of water production occurring upon withdrawal from the storage area. To solve this problem on an IBM 650 computer, several simplifying assumptions were made. Only the drainage portion of the relative-permeability curve was considered, thereby neglecting any hysteresis effects introduced during an imbibition cycle. The quantity of gas withdrawn during the period of the field study was only a small percentage of the total gas in place; therefore, this assumption is justified. For the problem studied in this case, the combined effects of capillarity and gravity were assumed to be negligible. This assumption becomes less valid as the formation thickness or the average pore diameter increases markedly. The justification of these assumptions should be studied in each case. THEORY This mathematical model for the prediction of saturation distribution and pressure for a radial gas-storage aquifer is based upon the equations for radial two-phase fluid flow&apos;.&apos; and for radial, unsteady-state, single-phase fluid flow.&apos; The two-phase flow is considered to take place in a "core" whose radius is equal to the maximum radius that the gas zone will attain, as indicated in Fig. 1. The injection-withdrawal history is approximated by a series of constant flow rates. Within the "core", gas is assumed to behave as a semi-compressible fluid; that is, the gas is assumed to have a constant density based on the average gas-zone pressure for the flow period. Between each constant-rate time increment, the gas density is allowed to change. This involves an iterative solution of the pressure equation. The liquid phase in this zone is considered to be incompressible, whereas it behaves as a compressible liquid outside of this region. The equations of incompressible flow in the "core" may be written for the gas phase