Search Documents

By F. G. Robb, D. N. Madge

This paper details some of the specifications and problems associated with implementing on-line control of a potash concentrator. The computer system described has the capability of three levels of control, including DDC with a single central processor. The system is also capable of on-line program compilation and modification through an extension of vendor utilities. The system is now in operation.

Jan 1, 1969

By J. L. Konchesky

The Bureau of Mines investigated the main reactions of coal gasification by steam. Coal (minus 70 mesh)-water slurry was fed into a laboratory-scale tube-coil gasifier to show the effects of residence time (0.6 to 1.1 sec), steam-coal ratio (2 to 11 lb/lb), and temperature (1,600° to 1,900° F), at about 1-1/2 atm. Devolatilization was very rapid: about 70 pet of the coal's volatile matter evolved within 0.6 sec at about 1,600° F; temperature increase released more volatile matter until about 95 pct was driven off at 1,900° F. Devolatilization was not increased by raising steam concentrations and was only slightly increased by extending residence time. Below 1,600° F, there was no evidence of the steam-carbon reaction, but at 1,900° F more than 80 pct of the coal's fixed carbon was converted to gas; the conversion was increased Significantly by increased residence time and steam concentration. The water-gas shift reaction occurred simultaneously with devolatilization, and the steam-carbon reaction, approaching equilibrium with increased residence time and temperature and decreased steam concentration.

Jan 1, 1969

By R. K. Datta, T. S. Davis

YVO4:EU, excited under 3650A radiation, shows a marked increase in total brightness with in~rease in temperature. At 25°C the response to 3650A radiation is about 48 pct of its maximum response which occurs at about 350°C, Fig. 1. In an attempt to explain the phenomenon, one of the present authors' (R.K.D.) investigated the diffuse reflectance spectra of the phosphor at different temperatures (up to -150°C) and was unable to observe any significant change to which the increased brightness could be reasonably ascribed. Since then diffuse reflectance spectra at elevated temperatures have been obtained in this laboratory using an improved device. A significant decrease in the diffuse reflectance of a phosphor in the ultraviolet region—especially in the longer wavelength region—at -300°C has been observed. The purpose of this note is to describe the device used for high-temperature reflectance measurements and to explain the increased brightness of YVO4:Eu phosphor at elevated temperatures under 3650A excitation. Experimental. YVO4 and YVO4:Eu samples were prepared as previously described.' The elevated-temperature reflectance spectra were obtained with a Cary 14 spectrophotometer equipped with a diffuse reflectance accessory and a ring collector. CaF2 was used as the standard. The cell used

Jan 1, 1969

By R. D. Nelson, S. D. Dahlgren

The conditions of temperature, strain rate, and total strain favoring deformation by grain boundary sliding, slip, or deformation with concurrent recrystallization were evaluated for alpha plutonium. Grain boundary sliding and slip were studied by examining structures formed on polished surfaces with 4.5 pct compressive deformation. The conditions required for recrystallization were determined from compressive deformation us time curves. At 110" C and at a strain rate of 10&apos;per min, deformation was almost exclusively by grain boundary sliding whereas slip was the predominant deformation mode at -10°C and lo-&apos; per min. Defarmation at intermediate tenzperatures and strain rates produced structures showing mixtures of grain boundary sliding and slip. Recrystallization occurred concurrently with deformation only after a critical strain was reached. About 14, 6, and 3 pct strain was required at 105", llO°, and lZO°C, respectively, before recrystallization started, irrespective of the strain rate in the range of lo-&apos;1 to 10&apos;4 per nzin. L HE temperature and strain-rate dependency of the modes of plastic deformation were evaluated for high-purity as-cast alpha plutonium. The techniques used to study slip and grain boundary sliding were similar to those employed previously for alpha plutonium by Bronisz and Gorum,&apos; and spriet.&apos; Recrystallization was investigated using the methods reported by Nelson.3 Bronisz and Gorum&apos; found that slip occurred on more than one slip system at room temperature, and suggested that grain boundary sliding also contributed to the deformation. spriet2 found that deformation was predominantly by slip at room temperature, and observed one to three orientations of slip traces in individual grains. In addition, he deformed polished samples at 10O° C, but concluded that the deformation character of alpha plutonium at 100° C and at room temperature were not essentially different. Several investigators reported that twinning was only occasionally observed.1&apos;2&apos;4 Nelson3 recently found in creep tests that high-purity alpha plutonium would recrystallize concurrently with compressive deformation at temperatures between 25" and 115°C. EXPERIMENTAL PROCEDURE Electrorefined plutonium having less than 300 ppm total impurities was received from Los Alamos Scientific Laboratory in the form of +-in. diarn cast rods. Major impurities were americium, <I00 ppm, and tungsten, (60 ppm. Less than 25 ppm each of other impurities were present. The rods were cut into right half-cylinders 0.35 in. long and 0.25 in. in diam, and the flat faces along the cylinder axes were metal-lographically polished. To avoid loss of the metallo-graphic polish by oxidation, the polished samples were deformed in a silicone fluid heat treating medium contained in a chamber evacuated to a pressure of lo-&apos; torr. Deadweight loading was used to compress the samples between 4.5 and 20 pct at strain rates between 10"4 and lo- &apos; min-&apos;. Strain rates were determined by dividing the total deformation by the elapsed time of the test, even though the strain rate was not constant during deformation. Curves of strain vs time for samples deformed with a constant compressive stress3 show that the strain rate changes by less than a factor of four during the testing of a given sample, which is small compared to the three orders of magnitude change in strain rate investigated. Two samples were loaded into an ordinary metallographic mounting mold and impacted by striking the mold piston with one blow of a 2-lb hammer. Deformation temperatures for the two impact samples were 25" and 90°C. The samples, following deformation, were washed in carbon tetrachloride and their polished and deformed surfaces were immediately inspected metallographi-cally. Subsequently, the deformation structures were replicated with acetate replicating tape for later electron microscope examination, and the samples were then repolished and etched to reveal the microstruc-tures of the deformed metal. Second-stage carbon replicas were prepared from decontaminated acetate replicas using the procedure developed by Miller, Bierlein, and astel.&apos; Additional samples, which were 0.35 in. long by 0.25 in. in diam, were used to obtain constant load deformation vs time curves. These samples were deformed 10 to 20 pct with loads of 15 and 20,000 psi at 120°C; 30, 45, and 55,000 psi at 110°C; and 20 and 30,000 psi at 105°C. Samples were polished and etched to determine whether or not recrystallization or twinning had occurred during deformation. RESULTS AND DISCUSSION The experiments showed that alpha plutonium deforms by grain boundary sliding, slip, and deformation with concurrent recrystallization. Deformation by twinning was insignificant. Metallographic examination of polished and deformed surfaces revealed grain boundary sliding and slip. The start of recrystallization was ascertained from constant load deformation vs time curves. Deformed surface structures of samples compressed 4.5 pct under four different conditions of temperature and strain rate, Fig. 1, show how the deformation mode changes from predominantly grain boundary sliding to predominantly slip with increasing strain rate and decreasing temperature. Deformation at 110°C and 10" 4 min-&apos;, Fig. l(a), was predominantly by grain boundary sliding. Only a few slip traces can be seen in Fig. l(a).

Jan 1, 1969

Jan 1, 1969

By John Ward Smith

The Bureau of Mines derived equations expressing the theoretical relationship between density of oil-shale rock and its organic content. Effects on theoretical shale density produced by variations in mineral composition and organic matter properties were evaluated. The theoretical relationship was extended to relate oil-shale density to oil yield in both oil weight and volume units by considering the fraction of organic matter converted to oil and the specific gravity of the oil produced. Corollary equations relating organic content by volume to organic content by weight and to oil yield were developed. The general relationships apply to an oil-shale deposit uniform in character. Application of the theoretical relationships to Green River Formation oil shales was demonstrated.

Jan 1, 1969

By Johnson M. W

Factual data about the sedimentation in the Newcastle Coal. Measures are presented. Data concerning the total coal measure sequence (isopachs, sand /shale ratio and total coal thickness variation), data about three coal seams (Borehole, Dudley and Victoria Tunnel) and data about interseam con- glomerate units are discussed. From this data conclusions are drawn about the coal measure environment. The western edge of the coal basin, the Lochinvar Dome, although a high feature was apparently stable during coal measure time. Despite its height it apparently supplied little sediment to the basin. The northern and north-eastern edge, which was unstable, suppl- ied most of the coarse detritus including the conglomerates to the basin. The conglomerates show patterns of fluvial origin but after a period of "normal" variations in pattern they become restricted to two parallel zones evidently tectonically controlled. Deposition of conglomerates is episodic and separated by periods of coal formation when deformation was apparently slight. Formation of each of the three seams considered followed a rather similar pattern - each shows development of thick coal around the basin edge and poor thin coal in the "axial" region of the present synclinal structure. Evidence of in-situ formation of some areas of the seams is discussed and periodic deformation of the basin during coal measure time is postulated.

Jan 1, 1969

By Robert M. Lewis

Chemically, mica is a complex silicate of sodium, potassium and aluminum. There are no important uses of mica where the chemical analysis is important, except that in some uses of ground mica, the content of iron oxide present as an impurity must be kept to a very low limit. Muscovite mica is frequently found containing black and red "stain" as a homogeneous part of the books or sheets. In most cases this stain or impurity is magnetite or hematite. Some mica contains impurities deposited by infiltration. This is commonly a red stain, less often brown or white. The best mica is free from all such impurities, and the greater their content, the less value does the mica have.(1) It is because of these undesirable impurities that a large segment of the fine ground mica industry relies on weathered pegmatites and alaskite ores as a source of raw material. The mica concentrates obtained from these ores are relatively clean and require moderate beneficiation to meet color and density specifications. One need only to observe the stained and decrepitated condition of the mica recovered from mica schist to appreciate the challenge that processing this material presents. Despite these adversities there are several advantages to be gained by mining and processing mica schist ores. The mica content of these ores are considerably higher than the pegmatites and alaskite ores now being mined, often containing in excess of forty percent mica. This-higher grade mica feed should promote a lower mining and milling cost with its resulting competative advantage. With more rigid government

Jan 1, 1969

By Charles R. Simcoe

The stress-strain curves in the microyield region for various Fe-C solid solutions shown by Davies and Ku are typical examples of parabolic curves defined by the expression a = kEn 111 These data, when replotted on log-log coordinates, as shown in Fig. 9, exhibit three essentially parallel straight lines. It is apparent that carbon increases the flow stresses, but does not affect the exponent n. In other words, the only factor in Eq. [1] affected by in in carbon content is the constant k. As Eq. [1] is written k is the stress at a strain of 1.0. This is unsatisfactory for examining data in the microyield region. Therefore, the equation can be altered to the form: to permit the use of the authors&apos; microyield stress at a strain of 10-6 as k. Extrapolating the curves in Fig. 9 to a strain of 10-6 gives stress values of 2000, 3200, and 5400 psi for 0.008, 0.014, and 0.02 pct C, respectively. These values do not agree exactly with the authors&apos; measured values, but this will not affect the following discussion. The effect of carbon content on the slopes of the stress-strain curves shown in Fig. 2 of the paper by Davies and Ku can be analyzed by differentiating Eq. [2] with respect to strain. Thus: where n is always less than 1.0. Since n is essen- tially constant as a function of carbon content for the data under discussion, the slopes of the stress-strain curves depend only on k and E. If the slopes are examined at various selected Plastic strains. i.e., 10-6, 10-5, and so forth, it follows from Eq. [3] that the dif-ferent slopes for the different carbon contents depend only on the different values of ke=10-6. Thus, the larger ke=10-6, the steeper the linear stress-strain curve will rise with increasing plastic strain. This

Jan 1, 1969

A pyrite flotation concentrate assaying 0&#183;065 oz/ton gold was investigated to determine the mode of occurrence of the gold. Ftom cyanidation and acid leaching experiments it was concluded that O&#183; 02-0&#183; 03 oz/ton goldwas in solid solution in pyrite and the balance of the gold was finely dispersed as grains throughout the pyrite.Autoradiography and electron-probe examination of possible areas in briquetted concentrate did not detect any grains of gold. Neutron irradiation and gamma spectrometry confirmed on 100 mg samples the fire assays obtained on 1 assay ton samples.Oxidation to hematite and reduction to magnetite at 600&#176;C made practically the whole of the gold extractable by cyanide. This technique could be of major significance for gold ores which contain a substantial amount of &quot;refractory&quot; gold.

Jan 1, 1969

A great range of contacting equipment is used in solvent extraction. However, a tower filled with a regular packing offers many advantages in the direct extraction from slurries. A study was made of the flow of a dense slurry through a packed tower in which the continuous phase was kerosene. The grid type packing, which was preferentially wetted by the slurry, had slats inclined at an angle of 30¦, to the horizontal. On the sur- face of the slats there were uniformly spaced longitudinal indentat- ions. These were of a height slightly less than the film thickness and divided the slurry discharging from the bottom of the slats into a number of uniform streams. When slurry entered the tower it was distributed by means of a troughs and weirs type of distributor. Because of the relatively high slurry velocities that occur at the entry region, it was necess- ary to have the distributor and the first two grids above the solvent level. The interface at the bottom of the tower was located at the mid-point of the bottom grid, so that the slurry entered the interface as a continuous film across the full width of the slats. Using the above system, stable flow with low entrainment was achieved at slurry flow rates of 3.4 mis/cm/sec. This would be equivalent to a flow rate of 177 gallons per minute through a tower of cross-section 4 ft 4 in. by 4 ft 10 in., or 1200 tons per day of solids of specific gravity 2.6 when the pulp density was 65 per cent. Previous tests indicated that a total tower height of about 30 ft would be required for a recovery of better than 99 per cent.

Jan 1, 1969

By R. E. Reed-Hill, W. H. Hartt

High-purity magnesium single crystals, oriented with basal plane parallel to stress axis, were deformed in tension at room temperature so as to form second-order (1011)- (1012) twins. Investigation by optical and replica electron microscopy revealed that these twins can support very large internal strains approaching 1000 pct and that fracture follows as a direct consequence of these strains. At low strains the basal plane within (10i1)- (10i2) twins rotates with increasing deformation until it becomes approximately parallel to the twin boundaries. This rotation is not consistent with deformation by basal slip in the double twin, but may be rationalized in terms of nonbasal slip. Observations also indicate that basal slip in the primary (1011) twin, prior to second-order twinning, may contribute to this rotation. Deformation markings that are difficult to interpret in terms of slip on previously observed systems were noted in the twin bands. A possible means of rationalizing these markings is to relate them to grain boundary shear. The fracture associated with these twins is thought to initiate by formation of voids at the twin boundaries. These grow into microcracks, with ultimate fracture occurring by tearing of the interconnecting regions. THE significance of deformation and fracture along second-order (10i1)-(10i2) twins in polycrystalline magnesium1-&apos; and many magnesium alloys4 was discussed in a previous paper.5 It was noted that this double twinning mode effectively controls the room-temperature deformation and fracture of magnesium single crystals, oriented with the basal plane parallel to the tensile axis. Reed-Hill and Robertson6 observed that, although these crystals may fracture at a macroscopic strain of less than 1 pct, shear strains up to 1000 pct could occur within (1011)- (10i2) twin bands. They concluded that the fracture was "ductile" and resulted from concentrating the deformation into a very small volume. The present research was undertaken to obtain a better understanding of the deformation and fracture associated with these twins. EXPERIMENTAL TECHNIQUES The experimental procedure involved straining high-purity magnesium single crystals in tension at room temperature. Specimens were oriented so that the applied stress direction was within 2 deg of the [1010] crystal axis. The specimen cross section was rectangular with faces closely parallel to (1210) and (0002). The specimen preparation and testing procedures have been described in detail elsewhere.5 Second-order (1011)- (1072) twins nucleated during straining were studied using optical and replica electron microscopy techniques. For the latter, cellulose acetate-carbon double replicas were employed. Observations were made on the (1210) crystal surface, which is the plane of shear for twins of this type. In addition to the single crystals, some investigations were carried out on longitudinal polycrystalline specimens. These were obtained from high-purity magnesium plate with a texture in which the basal planes of the grains tended to be nearly parallel to the specimen axis. All. polycrystalline specimens were annealed prior to testing in order to produce a coarse grain structure, permitting Laue back-reflection X-ray photographs to be obtained from grains of interest. EXPERIMENTAL RESULTS a) Deformation Within ( 101 l )-( 1072) Twins. Fig. 1 is an optical micrograph demonstrating the extent of the plastic deformation that can occur in (1011)-(10i2) twin bands. Notice at the lower left corner of the photograph the large displacement of the dark band (a polishing step) where it crosses the twins. The shear strain in the twin band at this position was computed to be 700 pct, using the ratio of step width to twin thickness. It was previously shown5 that (1011)- (1012) twin bands are often composed of small, separate twins aligned close to the macroscopic habit plane. Fig. 2 is an electron micrograph from near the upper end of the twin band in Fig. 1. This illustrates that the twin habit often does not coincide exactly with the macroscopic band habit, but may be inclined at a slightly smaller angle to the matrix basal plane. Figs. 3 through 6 are electron micrographs of the twin band in Fig. 2 illustrating successive steps in the development of deformation inside the band. A twin band such

Jan 1, 1969

By C. A. Stout, S. W. Nicksic

The materials that stubilize emulsions of some Southern California crude oils were isolated from the produced crude. These substances were separated into acidic and nonionic fractions by ion exchange chromatography. Each fraction was further broken into groups of homogeneous molecular weight composition by gel permeation chromatography. Both ionic character and molecular weight of the individual fractions are reflected in the properties of the original crude oil emulsion. Specifically, the more acidic and the higher molecular weight fractions appear to be the most effective stabilizers. INTRODUCTION The treatment of oilfield emulsions is a continuing problem for the oil industry. In almost all cases, these are water-in-oil emulsions. Much more is known about the formation and stabilization of oil-in-water (O/W) emulsions than of the water-in-oil (W/O) type. To understand W/O emulsions, more information is needed on the materials responsible for their formation and stabilization. This is especially true in the case of oilfield emulsions where the stabilizing materials are complex mixtures of large asphaltene-like molecules. To study the chemical and physical properties of these emulsion stabilizers, it is necessary first to separate them from crude oil. This separation was achieved in the case of two Southern California crude oils. The materials that were isolated appear at first to be little different from an asphaltene mixture. However, if a mixture of distilled water and mineral oil containing about 1 percent of the natural emulsion stabilizers is passed through a colloid mill, the W/O emulsion formed may be stored in the laboratory for many months without any sign of breaking. This could not be done using undifferentiated asphaltenes from the same wells, indicating that there is a real difference between these materials and ordinary asphaltenes. The colloid mill used to generate the emulsion was required because tile natural emulsion stabilizers had very little effect on oil-water interfacial tension. The one characteristic — complexity of composition — that emulsion stabilizers share with asphaltenes makes virtually impossible the exact structure determination of the components of the mixture. For this reason, it was decided to investigate first the influence of molecular weight on the ability of these materials to stabilize emulsions. The emulsion stabilizers were first separated into acidic and neutral fractions according to whether or not the material was retained on a weakly basic ion exchange column. The extraction procedure was not designed to remove basic components from tile oils, and no attempt to isolate such a fraction was made. Following this separation, selected fractions were analyzed by gel permeation chromatography. Gel permeation separates according to molecular size, largely independently of chemical or polar characteristics. A heterogeneous mixture is separated into a number of fractions having much more homogeneous molecular weight distribution. A number average molecular weight of each fraction is therefore close to the true molecular weight of the principal component of that fraction. The fraction weights, expressed as weight percent of the sample, and the number average molecular weight of a representative selection of the fractions, will be referred to in this paper as the molecular weight profile of the original sample. DESCRIPTION OF GEL PERMEATION ANALYSIS The first reported gel permeation analyses were for separating partially hydrolyzed protein molecules in an aqueous system using expanded starch granules. Most of the theoretical work on gel permeation has been done by biochemists.1 This theoretical treatment is still incomplete, but the

Jan 1, 1969

By Barry J. Hansen

In the past six years, the Hanna Mining Company has increased its pellet producing capabilities greatly by the construction and expansion of six pelletizing plants in North America. At the present time Hanna managed pelletizers have a nominal capacity of 18 million long tons per year. This capability makes Hanna one of the world&apos;s largest pellet producers. To create and maintain such an operation has required the skill and dedication of many people in all aspects of the company&apos;s activities, in the various fields of exploration, operations and research. The Research and Development Group consists of the staffs of two laboratories situated on the Mesabi Range in northern Minnesota. These laboratories are well equipped to carry out work in the fields of beneficiation and extractive metallurgy especially iron ore treatment. Facilities exist for the batch and semi- continuous investigation of the unit operations common to iron ore beneficiation and pelletizing. The R and D Division is divided into groups according to specialized areas of interest. Presented in this paper today are some aspects of the activities of the Pelletizing Group of the Research & Development Organization.

Jan 1, 1969

By E. H. Boswell, E. M. Cushing

The Mississippi embayment east of the Mississippi River is part of a vast geologic and hydrologic province. The abundance of water available is largely responsible for the agricultural history of the area and is perhaps the key factor in the rapidly accelerating trend toward industrialization. The multiple artesian aquifer systems that underlie the region and the alluvial aquifer along the Mississippi River are virtually undeveloped sources of fresh water in that they are yielding only a fraction of the water available. The potential yield of the 12 aquifer systems is estimated to be about 25,000 mgd (million gpd) of which only about 600 mgd is presently being withdrawn. The amount of water in storage in these systems is about 1600 trillion gal. Ground water having a dissolved-solids content of less than 500 ppm is generally available at depths of less than 1000 ft, and water having a dissolved-solids content of less than 1000 ppm is available in some places to depths of more than 3000 ft. Iron is the most common troublesome chemical constituent in the ground water. Water in varying amounts is also available from streams within the region. The amount of water which originates within the region and which leaves it as streamflow during a year averages about 50 million acre-ft (about 45,000 mgd). An additional 375 million acre-ft (about 340,000 mgd), which originated outside the region, leaves the region as streamflow during an average year. Some streams have considerable flow during periods of no precipitation and the water during these periods generally has a low dissolved-solids content. Large quantities of water from streams and waterbearing sands (aquifers) are available in the Mississippi embayment east of the Mississippi River. Most of the region is underlain by aquifer systems that will yield large quantities of water to wells, making ground water the most readily available source of fresh water. About 90% of the water used in the region is obtained from wells. Streams have been used primarily to dispose of, and dilute, wastes that have been dumped into them. The Mississippi, Ohio and Tennessee rivers, and the lower parts of the Tombigee and Yazoo rivers also serve as a means of transportation. This nominally water-rich region receives from 48 to 56 in. of precipitation in average years and is crossed by several large rivers. The largest, the Mississippi River, is the western boundary of the region discussed in this paper. The economy of the region is basically agricultural, but industry has become increasingly important since World War 11. The diversification and mechanization of agriculture have resulted in a surplus of manpower, most of which is being absorbed by the expansion of industry. Further industrial expansion probably will occur because of the available manpower, natural resources and economical transportation. Increasing demands on available water supplies will be made and the future economic growth of the region will be largely dependent upon wise utilization and management of the water resources.

Jan 1, 1969

By U. Paakkinen, P. Eerola

In connection with theoretical studies of the dynamics and control of the flotation process, flotation cells and circuits at the Pyha- computer control of selective Cu - Zn - Fe sulfide flotation in the latter half of 1969. In this paper, the practical conclusions of these studies, conducted through simulation techniques, are summarized. The application of these results for control purposes implies a knowledge of the dependence of the flotation rate coefficients and the grade of concentrate on several factors, including aeration rate, froth depth, and reagent consumption. The relationships found in laboratory and plant scale tests are described. Off-line optimization studies, based on the relation- ships, are also described. The proposed steps in executing the computer control schedule are presented. In the first stage, a stabilizing control will be carried out using feed-back loops. A feed-forward control system which utilizes simple process models will also be tested. After these elementary stages, a proposed mathematical model of the circuit will be checked out and implemented with computer- run, plant-scale process tests. Off-line simulation will be continued in advance development of the control system. As a final step in the program optimal control of the process will be carried out through this model.

Jan 1, 1969

By Robert W. Ageton

The Bureau of Mines is currently conducting a comprehensive study to evaluate all aspects of recovering gold from large, low-grade deposits in the Golden Sunlight area, Whitehall mining district, Jefferson County, Mont. Prior to World War II this area produced significant quantities of gold. The study includes geologic mapping, diamond drilling, assaying, computer analysis of assay data, metallurgical research and ore testing, mine and mill planning, and an economic evaluation of all factors bearing on the establishment of a producing gold mine. Preliminary evaluation of the deposit indicates about 31 million tons of submarginal gold resources. American Exploration and Mining Co. (AMEX) has assisted the Bureau in this project by providing drill site preparation, control assays, and geological aid for the Bureau&apos;s exploration work and by furnishing comprehensive assay information of the property and copies of base and geologic maps. This report presents the geology and mining history of the area, together with a summary of Bureau work in progress.

Jan 1, 1969

By R. M. Geisenheyner

We are in the midst of a transportation revolution. From the Wright Brothers to the SST. From the Toonerville Trolly to the 120-mile-an-hour Northeast Corridor train. From Christopher Columbus to Neil Armstrong. There is another side to this revolution and it is going on every day at our mines and factories. We call it bulk distribution. Perhaps ordinary sounding, but every bit as exciting as those other transportation innovations. Basically, we have gone from bag to bulk, and in the process reached that Last Frontier in manufacturing cost savings. By revitalization of our distribution system, we have made it possible for our mines and factories to realize fully the vast economies of spage-age manufacturing techniques and accomplishments. The pacesetters in this distribution revolution have been you gentlemen, the producers and shippers of dry mineral products. In this area bulk distribution grew up. And here most of the current innovations are being introduced. Management of your industry was early to realize the economic advantages of bulk distribution. Distribution costs of your products comprises a sizeable, portion of the selling price and management is sensitive of any developments which may even slightly reduce the cost of transportation, marketing or the exploration of new frontiers.

Jan 1, 1969

By C. S. Land

Relative permeability functions are developed for both two- and three-phase systems with the saturation changes in the imbibition direction. An empirical relation between residual nonwetting-phase saturation after water imbibition and initial nonwetting-phase saturations is found from published data. From this empirical relation, expressions are obtained for trapped and mobile nonwetting-phase saturations which are used in connection with established theory relating relative permeability to pore-size distribution. The resulting equations yield relative permeability as a function of saturation having characteristics believed to be representative of real systems. The relative permeability of water - wet rocks for both two- and three - phase systems, with the saturation change in the imbibition direction, may be obtained by this method after properly selecting two rock properties: the residual nonwetting-phase saturation after the complete imbibition cycle, and the capillary pressure curve. INTRODUCTION Relative permeability is a function of saturation history as well as of saturation. This fact was first pointed out for two-phase flow by Geffen et al. 1 and by Osaba et a1.2 Hysteresis in the relative permeability-saturation relation also has been reported for three-phase systems.3 Since saturations may change simultaneously in two directions in a three-phase system, four possible relationships arise between relative permeability and saturation for a water-wet system. The four saturation histories of this system were given by snel14 as 11, ID, DI and DD. I and D refer to the direction of saturation change (imbibition and drainage), with the first letter of the symbol indicating the direction of change of the water phase. As used in this paper, the second letter of the symbol refers to the direction of saturation change of the gas phase, i.e., D and I indicate an increase and decrease, respectively, in gas saturation. Only a few three-phase relative permeability curves have been published. Leverett and Lewis5 published three-phase curves for unconsolidated sand, and Snell4 reported results of several English authors for both drairrage and imbibition three-phase relative permeability of unconsolidated sands. Three-phase relative permeability curves for a consolidated sand were published by Caudle et a1.3 for increasing water and gas saturations (ID). Corey et a1.6 reported drainage (DD) three-phase relative permeability for consolidated sands. Recently, Donaldson and Dean7 and Sarem8 calculated three-phase relative permeability curves from displacement data on consolidated sands, also for saturation changes in the drainage direction. The only published three-phase relative permeability curves for consolidated sands with saturation changes in the imbibition direction (11) are those of Naar and Wygal.9 These curves are based on a theoretical study of the model of Wyllie and Gardner 10 as modified by Naar and Henderson.11 Interest in three-phase relative permeability has increased recently due to the introduction of new recovery methods and refinements in calculation procedures brought about by the use of large-scale digital computers. The scarcity of empirical relations for three-phase flow, and the experimental difficulty encountered in obtaining such data, have made the theoretical approach to this problem attractive. RELATIVE PERMEABILITY AS A FUNCTION OF PORE-SIZE DISTRIBUTION Purcell used pore sizes obtained from mercury-injection capillary pressure data to calculate the permeability of porous solids.12 Burdine extended the theory by developing a relative permeability-pore size distribution relation containing the correct tortuosity term.13 The works of Purcell and Burdine were combined by Corey14 into a form that has

Jan 1, 1969

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969