Search Documents

By Willard P. Ward

Having been engaged during the past year in the manufacture of ferro-manganese in a blast furnace, I have imagined that some further information on this subject might be of interest to that large number of members of the Institute who are engaged in the iron and steel industries. The great question of the economical production of good steel, ingot metal, or homogeneous iron (we scarcely know now what to call it) from good materials, and at reasonable prices, has, thanks to the intelligence and energy of a few American engineers, been successfully solved in this country. Nearly all first-class railroads, and many that would not rank so high, have adopted the Bessemer rail. Martin boiler-plates are rapidly gaining ground, in competition with

By G. Mark Aylmore

Various processes have been developed using a combination of elevated temperature and chemical treatment processing to recover Li from silicate minerals. To facilitate further process development, a comprehensive understanding of the deportment of Li and associated minerals in ore bodies is essential to allow the industry to predict the response of ore reserves to metallurgical treatment options. This paper describes results from the integrated use of the John de Laeter Centre’s state of the art analytical and mass spectrometry techniques to characterise a selection of Li bearing ore bodies and determine their amenability to potential processing options for the extraction of Li. The mineralogy, mineral associations, and liberation characteristics of ore-bearing and gangue minerals were characterised using a combination of the TIMA and XRPD studies. The Li content and distribution within minerals were defined using LA-ICPMS and field emission scanning electron microscopy techniques (EBSD, ToF-SIMS) and atomic probe microscopy.

By Leiting Shen

The thermodynamics of tungsten ore decomposition in mineral acid and alkaline solutions were studied. The published thermodynamic data of tungsten minerals were collected and assessed. The Gibbs energies of CaWO4 (−1538.43 kJ/mol), FeWO4 (−1053.91 kJ/mol), MnWO4 (−1206.08 kJ/mol), H2WO4 (−1003.92 kJ/ mol), and Na2WO4 (−1455.58 kJ/mol, aq) at 25 °C were adopted in the calculation using HSC software. The results show that CaWO4 is decomposed more readily in Na2CO3 solution than in NaOH,while FeWO4 andMnWO4 aremore reactive inNaOH solutions. Froma thermodynamic point of view, tungsten ore decomposes easily in acid solutions despite most ΔGo T increasing slightly with temperature. Oxidizing Fe2+ to Fe3+ in acidic solutions facilitates decomposition of FeWO4, and the reaction of CaWO4 in H2SO4 solution occurs more easily than in other mineral acids, due to the formation of sparsely soluble CaSO4. The results fit well with the experimental data and industrial experience previously reported in the literature.

By Takahiro Furuzono

Saganoseki Smelter and Refinery of Pan Pacific Copper Co., Ltd. started treating copper anode slimes using a unique hydrometallurgical process in chloride media in 1997. The process comprises the following stages, decopperizing, wet chlorination, solvent extraction for gold, silver chloride reduction, and separate reduction stages that treat the gold extraction raffinate with sulfur dioxide gas. Precious metals and rare metals in the copper anode slimes are dissolved in an aqueous solution of hydrochloric acid by oxidative leaching, which is termed wet chlorination. Gold is efficiently extracted from the chlorination liquor by dibutyl carbitol in the solvent extraction stage. Gold is reduced by oxalic acid from goldloaded dibutyl carbitol and recovered as gold powder. The gold powder is cast directly into granules or ingots of more than 99.99% purity without electrorefining. Selenium is recovered as a selenium residue by reduction with sulfur dioxide gas and refined through vacuum distillation. The plant was designed to process 132 tonnes of copper anode slimes per month and has a current capacity of 206 tonnes per month due to several enhancements. Compared with the conventional process, the hydrometallurgical process described here provides some advantages, in particular, a higher product quality, shorter gold-retention time, and reduction in environmental load.

By William D. Monaghan, Michael A. Trevits

Over the period 2000 to 2003, roof falls have accounted for 4 to 14% of the fatalities in underground mining operations. The National Institute for Occupational Safety and Health (NIOSH) is conducting research to reduce the frequency, exposure, and risk of these events through an ongoing program of field and laboratory studies. One area of research involves the evaluation of polyurethane grouting technology that is commonly used to stabilize fractured mine roof strata. Concurrently, NIOSH is conducting work to evaluate the application of ground penetrating radar for advanced delineation of problematic mining areas. In this study, NIOSH partnered with Sub-Technical, Inc. who injected a polyurethane grout (also known as glue) into a roof area of NIOSH’s Safety Research Coal Mine. The mine roof area was scanned using ground penetrating radar technology before and after grout injection in an attempt to determine the extent of grout infiltration. A comparison of the pre-and post-grouting radar records showed a significant change at a mine roof depth of about four to five feet. The interpreted radar records were then compared with drill core information, borescope evaluation of roof-bolt holes in the study area, and underground observations. At this site, the interpretations of the radar records correlated with data obtained from the core holes, borescope evaluations and underground observations. This study showed that ground penetrating radar technology can be a useful tool for detecting changes in mine roof due to the injection of the grout.

By Gabriel S. Esterhuizen, Anthony T. Iannacchione

The roof rock in underground limestone mines in Northern Appalachia can be subject to high horizontal stresses in spite of the shallow depth of the workings. The high stresses can cause roof stability problems. The National Institute for Occupational Safety and Health staff have observed a distinctive asymmetrical mode of failure at the pillar-roof contact in two underground stone mines, which is different from the typical failure mode observed in response to excessive levels of horizontal stress. A dip of greater than 5° was identified as a possible cause of this mode of failure. Numerical model experiments showed that an increase in the dip of the workings can cause an increase in the stress at the up-dip corner of the roof beam. A case study is presented in which failure at the pillar-roof contact was observed where the dip of the workings was 7° in a high horizontal stress field. Numerical modeling showed that the relatively small dip of the workings could have induced stresses changes in the immediate roof that would explain the failures. The model results also showed that this mode of failure is less likely to occur if the limestone pillars contain weak bedding planes that provide stress relief. In addition, the results showed that for the conditions at the case study site, the stress in the roof beam was not sensitive to the thickness of the roof beam or to the excavation span. The high horizontal stresses at this site are an important contributing factor to the observed failures.

By I. Lowthian Bell

THERE has been probably no improvement introduced into the manufacture of iron which created more surprise in the minds of practical smelters and of scientific men than Neilson's discovery of the hot blast. In 1829, Messrs. Dunlop & Co. consumed at the Clyde Works, near Glasgow, nearly 79 tons of coal to make a ton of pig iron, of which about 20 cwts. was employed under the boilers of the blowing engines, leaving nearly 7 tons for the consumption of the furnace itself. In the year 1833, by heating the air to 612' F. (322' C.), they redoccd the 79 tons to 3 tons by the mere burning of 8 cwts. of coal in the apparatus for raising the blast to this temperature. The statenlent just given is strictly true in a commercial sense, but when we come to consider the question as one of physical science, it is necessary to eliminate some of the conditions which conduced to so extraordinary a result. Anterior to Neilson's time the fuel employed in smelting iron was coke, and it was supposed at that period, erroneously, however, that coal in its raw state, when burned with cold air, was totsally unfit for the purpose in question. This mineral, as it was obtained from the Lanarkshire coal-field, contained about 35 per cent. of volatile matter, but the process of charring it was performed in so unskilful a manner that instead of receiving 88 cwts. of coke from 135 cwts. of coal, 60 cwts. only was the product of the operation. Again, the reduced colisumption of carbon in the furnace itself was accompanied by a corresponding diminution of blast, and this was followed by an economy of half a ton in engine coal. After making a proper allowance, however, for all these collateral circumstances, we are within the mark when we admit that in actual coke, which is really the form the fuel has to assume before it is burnt in the blast-furnace, not less than 30 cwts. was saved by heating the blast to 615° F. (322° C.). This unexpected economy, be it remembered, was obtained by the combustion of the 8 cwts. of coal, somewhat wastefully applied in the hot-air apparatus. Numerous have been the opirlions advanced by different authorities in the scientific world to account for the apparent anomaly, bat none,

By T. Sterry Hunt

This remarkable mine, to which attention has lately been drawn, is situated not far from the New River, in Ashe County, North Carolina, on a spur of the Blue Ridge which lies between the main crest of this name and the so-called Iron Mountain, a part of the Unaka or Smoky Mountain range. It is on a true fissure-lode, which cuts the strata of gneiss and mica-schist of the region. These, which have the lithological characters of what I have called the Montalban or White Mountain series, have here a dip of about forty-five degrees to the southeast, while the lode is vertical, with a course north 60' east. Both the country-rock and the lode, as is usual in this region, are decomposed to considerable depths; and the latter exhibits a cap of hydrous peroxide of iron or gossan, blocks of which, left by the denudation of the softer inclosing rock, are scattered over the surface, and first called attention to the locality as a source of iron ore. In sinking upon this lode, it was found that the porous gossan at a certain depth became charged with carbonate and red oxide of copper, and lower down was replaced by sulphuretted ores of this metal of remarkable richness. This deposit was worked for copper, in an irregular manner, before the late civil war, but was afterwards abandoned, and only reopened during the present year, by Messrs. S. S. & J. E. Clayton, of Baltimore. The outcrop has been traced for 1900 feet, and the sulphuretted copper ores have been met with in five shafts in a distance of 661 feet, a drift having moreover been carried for the above distance through the solid ore. The breadth of the lode, so far as opened, varies from six to eight feet to fourteen feet in places, and that of the outcrop of gossan is, in some parts, twenty feet. After comparing this deposit with the copper mines of Ducktown, Polk County, Tennessee, and those of Carroll County, Virginia, it is evident that we have, at Ore Knob, to do with the same kind of black ores as were met with in these mines, and will be noticed further on, Beneath the gossan, at Ore Knob, is found, in some parts, an iron-black, friable, drusy, crystalline, sulphuretted ore, inclosing grains of quartz, garnet, and magnetite, besides a black non-magnetic mineral, not attacked by nitric acid. This ore contains about 36 per cent. of copper, and portions have the mineralogical charac-

By Stephen C. Tadolini, Steve Bessinger, James Pile, Christopher Mark

The San Juan Coal Mine, located near Farmington, New Mexico, supplies the San Juan Generating Station with more than 6 million tons of coal annually. To replace dwindling surface mine production, San Juan installed a longwall in 2002. The roof rock in much of the underground reserve consists of interbedded coals, carbonaceous shales, and mudstones. Roof support consists of fully-grouted roof bolts, both tensioned and non-tensioned, supplemented by modified cable trusses. Short-encapsulation pull tests conducted by San Juan’s geotechnical engineering department showed that some of the immediate roof layers provided variable bolt anchorage. To ensure that the roof bolts were obtaining the necessary anchorage, additional tests were conducted in a variety of locations. The results have been used to determine the optimum support design for different areas of the mine. Another series of tests evaluated the performance of several bolt and hole diameter combinations. The experience at San Juan Mine shows that short-encapsulation pull tests can provide the information necessary to improve ground control in variable roof conditions. It also illustrates how a pro-active geotechnical engineering program can be an essential component of a state of the art underground coal mine.

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

R. C. J. Goode (President): South Africa is known for its sunshine and mineral deposits-especially its gold and its diamonds-and as this year marks the 75th Anniversary of the establishment of the South African Institute of Mining and Metallurgy, it is appropriate that I say a few words about mining and metallurgy and the part this industry and our Institute has played in the life of this country. Africa, as the Dark Continent, with its fascinating tales of Prester John, Ophir, Queen of Sheba and its myths and legends has throughout the ages fired the public imagination. Pharaoh Necho, 600 years before the birth of Christ, sent an expedition from Egypt around the Cape to probe this strange land. The gold for King Solomon's throne was said to come from Monomotapa-the country now known as Rhodesia. Centuries later the Arabs and Portuguese searched for this mineral wealth. The Arabs found little gold but unfortunately established a lucrative trade in black slaves. Now let us take a brief look into the closing stages of the Stone Age when the indigenous population of the lower portion of this continent consisted of the Bushmen, a nomadic hunting people who knew no metals and kept no stock, and the Hottentots who were a pastoral folk with cattle and sheep. Then somewhere about the 11th century the Bantu started their long migration southwards and brought with them the knowledge of smelting of iron ore. One of their settlements was here in Johannesburg on the Melville Kopje at Emmarentia. These people were the first miners and metallurgists on the Witwatersrand and the beautifully preserved iron smelting furnace constructed nearly 1,000 years ago is a tribute to their technical ability. This ability to use iron to make better spears for hunting or to fashion hoes for cultivating the soil led to their establishment as a superior tribe. The earlier Bushmen were pushed back to the inhospitable swamps and deserts and the Hottentots fled to the south. The Bantu also knew the art of copper smelting. In the days of the Dutch East India Company, at the end of the seventeenth century, Bantu workers brought copper from the north-west Cape to Simon van der Stel, but it was another two hundred years before the white man turned this to advantage. Stories of the Bantu mining gold spurred the early European hunters to travel deeper into the heart of Africa in search of wealth, and we are told that Karel Kruger, in 1834 whilst leading an ivory hunting expedition, discovered gold on the Witwatersrand, and took samples to Cape Town; but when he returned with a larger expedition a few years later his party was attacked by Moselekatze's Matabele near Potchefstroom. Kruger was killed and only the legend survived. The real spur to mining in this country was the discovery of diamonds in Griqualand West in 1867. The field seemed fabulously rich and diggers flocked to the scene from the four corners of the earth. With them came Cecil Rhodes, Barnato, Beit and some of the finest mining and financial brains the country had yet seen. These giant entrepreneurs were not tempted by the small deposits of gold in the Barberton and Pilgrim's Rest areas when these were made a few years later but

By J. Durkin

The Bureau of Mines has conducted field studies in coal mines throughout the United States to determine the effectiveness of electromagnetic techniques in locating miners trapped underground following a mine accident. Data from these tests have been used to generate models of expected signal and noise distributions as found above these mines. These distributions have aided in placing the expected performance of a through-the-earth electromagnetic communications technique into a probabilistic framework. Results show that at a 10% false-alarm rate, the expected probability o detecting a miner's signal" from a depth of 1000 ft is 54% at 500 ft it is 95%. These depths exceed the actual depths of 90 and 50%, respectively, of United States coal mines. Sensitivity studies have shown that at depth of 1000 ft, the probability of detection will improve approximately 2% for each dB of increase in signal-to-noise (SNR) ratio.

By Alastair Grogan, Ron Elliott, Dale MacLean

Blasting for armour stone presents many unique challenges. Often, blasters get themselves into trouble when they assume that they can simply modify a production blast design to produce the required fragmentation size. They bid on the project using typical production blasting costs, and then find that they cannot produce the required fragmentation size and quantity with an acceptable waste factor. This paper outlines industry best practice guidelines for the production of riprap or armour stone. The paper examines the importance of initial geological mapping to determine the average in-situ block size, and the advantages to be gained from working with the geology to maximize production of the required block sizes. The paper also outlines the importance of proper explosives selection as well as blast designs that are unique to the production of riprap. The paper reviews challenges and experiences realized from multiple projects in differing rock types.

By Pete G. Kovalchik, Lynn A. Alcorn, Rudy J. Matetic, Adam K. Smith

Over-exposure to noise remains a widespread, serious health hazard in the U.S. mining industries despite 25 years of regulation. Most other categories of illnesses and injuries associated with mining have improved, with the exception of hearing loss. Mine Safety and Health Administration (MSHA) coal noise sample data collected from 2000 to 2002 show that 65% of the equipment whose operators exceeded 100% noise dosage is comprised of seven different types of machines. The continuous mining machine is first among all the equipment with 35% of the noise overexposures. The ational Institute for Occupational Safety and Health (NIOSH) is conducting research to reduce excessive exposure for operators of continuous mining machines and preventing additional cases of oise Induced Hearing Loss (NIHL) by developing low-cost retrofit noise controls for mining equipment. This paper describes a noise control for reducing the noise overexposures of continuous mining machine operators. Underground results show a 26.4% and 27.4% noise exposure reduction for the continuous mining machine operator. This research is providing the mining community with an additional noise control to be utilized on continuous mining machines, therefore reducing operator noise overexposure. Utilizing this newly developed noise control, along with previously proven controls for the continuous mining machine chain conveyor, will provide operators of these machines an opportunity to be within the MS HA-Permissible Exposure Limit (MSHA-PEL).

By Pingchao Ke

A crystalline polyferric sulfate (PFS) adsorbent was synthesized by oxidizing and precipitating ferrous ions under atmospheric conditions. The morphology, structure, composition and specific surface area (SSA) and adsorptive efficacy to As (III) of the adsorbent were characterized by SEM/EDS and TEM images, XRD patterns, FT-IR spectra, BET SSA analyses and adsorption experiments. The adsorbent presents as morphology of near-spherical aggregate and good crystallinity. A small amount of goethite co-precipitated with PFS in the case of the initial ferrous concentration of 1 mol/L and increased the SSA of the aggregates. The absorbent with co-precipitation of goethite showed good adsorption to As(III) and good filtering performance in the high As(III)-content solution of 10–100 mg/L under acidic, neutral and alkaline conditions (pH 2.09–9.01), and its maximal adsorption capacity and adsorption efficiency were 7911 μg/g and 96% separately. After Adsorption, the stability of the absorbents adsorbing As(III) were evaluated by TCLP tests, and the results indicated that the absorbents adsorbing As(III) were extremely stable, and the concentration of arsenic in all leaching solutions were less than 0.01 mg/L.

By W. R. Reed, J. A. Organiscak

Surface mining operations use large off-road haul trucks extensively to move material at mining properties. Past research, using the United States Environmental Protection Agency’s (EPA) emissions factors for unpaved haul roads, has shown that haul trucks generate the majority of dust emissions from surface mining sites, accounting for approximately 78%-97% of total dust emissions. Observations of dust emissions from haul trucks show that if the dust emissions are uncontrolled, they can be a safety hazard by impairing the operator’s visibility. This increases the probability for haul truck accidents (See photos above). However, the greatest long-term health hazard of dust generated from hauling operations is due to inhalation of the respirable dust [median diameter <4 micrometers (µm)] and tho­racic dust, which is equivalent to the EPA’s definition of PM10 [particulate matter with a median diameter <10 µm]. Exposure to respirable dust has long been considered a health hazard at surface mining operations, especially if silica dust is present. There are two legislative acts that regulate the air quality for mining operations: the Federal Coal Mine Health and Safety Act of 1969 and the Clean Air Act of 1970. The Federal Coal Mine Health and Safety Act of 1969 established the limits for dust in the work place for health and safety pur­poses. The Clean Air Act of 1970 regulates air emissions from facilities from an envi­ronmental perspective.

By C. F. Weinaug

In recent years, several methods have been proposed for correlation of physical properties and behavior of hydrocarbon systems. These correlations were found to be unsatis. factory for predicting differential gas liberation and shrinkage data now obtained experimentally. To achieve the desired accuracy, methods were developed to adjust the correlations with two experimental pieces of data. The required data are the saturation pressure at reservoir temperature and the specific volume or density at these conditions. The saturation pressure is used to adjust the vapor-liquid equilibrium-ratio correlations. The density measurement is used to determine the apparent density of liquid methane. With these adjustments, differential gas liberation and shrinkage are predicted within the experimental error. EQUILIBRIUM RATIOS Equilibrium ratios published by the National Gas Assn. of America&apos; were used as a basis for the calculation. The convergence pressure for the sample was determined by using a hypothetical binary system, with methane as one component and the remainder of the system lumped together to form the other component. For the critical temperature of the heavy constituent, the weight average of the critical temperatures of the ethane and heavier components was used. The correlation of pseudo-critical temperature, density and molecular weight of Matthews, Roland and Katzl was used in this calculation for the hexane and heavier fractions. Points for the critical properties of methane and the hypothetical heavy constituent were placed on a chart of the critical loci of hydrocarbon systems. A critical loci for the hypothetical binary system was drawn, and the critical pressure of this system was determined at the reservoir temperature. This critical pressure became the convergence pressure for which equilibrium ratios were read from the NGAA4 charts. The NGAA equilibrium constants for C, through C, were used for the C,; through C, fractions since they were separated on the bases of normal paraffin boiling points. The correlation of Hoffman, Crump and Hocott5 as used to estimate the equilibrium constant for the C fraction. The NGAA eauilibrium-ratio charts are based primarily on paraffin systems and have to be corrected for the presence of naphthenes and aromatics. The correlations in the literature proved too unsatisfactory for determining these corrections. EXperimental data, the saturation pressure, is used with the correlation of Soloman2 to make this adjustment. The correction factors of Soloman were cross-plotted so that the factors for ethane and heavier components were made a function of the factor for methane. a2 = 0.7927 at + 0.1646 .........(1) where a, = correction coefficient for for methane NGAA equilibrium ratio, and a2t = correction coefficient for ethane and heavier NGAA equilibrium ratios. This equation and the equation for the bubble point are solved simultaneously for the correction coefficients. The bubble-point equation is anKnXn = 1.0 . . . . (2) The equilibrium constant for the saturation pressure at reservoir temperature is used. This procedure not only corrects for the relative volatility between components, but also assures that the results will be directly related to the experimental saturation pressure. METHANE DENSITY After determining a set of equilibrium ratios, the apparent density of condensed methane at 60°F and 14.7 psia was estimated. This was found to be necessary because the correlations of Standing3 ere found to be unsatisfactory. The law of additive volumes was employed to make this calculation.

By H. Engelmann

The very extensive development of a brown coal formation, in the region of the Rocky Mountains, is well known to all of yon. The existence of these coals was known years ago, but they were of no practical importance until within a few years, when the mineral wealth of that country began to be developed and railroads were built. As their qualities became better known, their adaptability to domestic manufacturing and metallurgical purposes was more appreciated. They are still much undervalued, although they are mined at some points on a very extensive scale. A large portion of this coal formation is of cretaceous age, while it appears to extend into the tertiary epoch; but whatever difference there may be in the ages of these coals of different localities, they all have some common characteristics which distinguish them from our Eastern stone coals of the

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing