Search Documents

By J. Wade Watkins, E. S. Mardock

The accurate evaluation of reservoir-performance characteristics in the secondary recovery of petroleum by water flooding requires use of a water tracer that may be injected into water-input wells and detected at oil-production wells to supplement data obtained from core analyses, wellhead tests, and subsurface measurements. Radioactive iodine has been used successfully as a water tracer in field tests to determine: (I) relative rates and patterns of flow of injected water between water-input and oil-production wells and (2) zones of excessive water entry into oil-production wells. Laboratory evaluations of potential water tracers, previous tracer studies, the value of using a radioactive tracer, general field procedures, and the use of surface and subsurface instruments for the detection of the emitted gamma radiation, are summarized. Data from the field tests are presented graphically and discussed in detail. It is concluded that the radioactive-tracer method, using radioactive iodine, may be used successfully to measure either the relative rates and patterns of flow or zones of excessive water entry into wells under conditions of comparatively rapid transit time between wells. INTRODUCTION Extensive use is made of data obtained from core analyses, wellhead tests, fluid characteristics, and subsurface measurements in evaluating the sweep efficiency of water injected into oil sands for the recovery of oil. Theoretical flow rates and patterns may be calculated from those data using radial-flow formulas, if it is assumed that the physical conditions in the productive formation are homogeneous. Unfortunately, homogeneous conditions rarely, if ever, exist in oil-productive formations. The use of a tracer that may be injected into an oil sand and detected quantitatively, or even qualitatively, at offsetting oil-production wells provides basic data that may be used in determining more accurately the subsurface rates and patterns of flow of injected water between wells than is possible by theoretical calculations based on assumed conditions. Consideration of the data obtained by using a water tracer assists in the application of remedial measures to water-input wells, such as plugging of channels, or selective plugging of highly permeable zones, thereby effecting a more uniform flood and a greater ultimate oil recovery. A water tracer should have the following characteristics: (1) low adsorption on solid reservoir material; (2) high solubility in water; (3) low or negligible solubility in crude petroleum; (4) a wide range of solubility of compounds that may be formed by chemical reactions with ions present in reservoir rocks or waters; (5) high detectability in low concentrations by portable apparatus; and (6) general availability at low cost. Additional desirable characteristics are that the tracer should be nonhazardous in nature under normal working conditions and quantitatively detectable in the well bore, as well as at the surface. Review of Previous Tracer Studies The need for a water tracer has been recognized for many years. As far back as 1906, Dole1 described the use of fluorescein in flow-rate studies made in France in 1901. In 1921, Ambrose" discussed in detail the use of organic and inorganic dyes, chlorides, nitrates, and other anions, lithium salts, and the Slichter electrical

Jan 1, 1955

By C. A. Bruch, J. F. Duffey

On the basis of metallographic and cooling-rate studies of 15 U-Zr alloys, a new version of the 6 phase region of the U-Zr phase diagram is proposed. IN several published diagrams of U-Zr alloy system:' there is good agreement that at elevated temperatures complete solid solubility exists between y uranium and ß zirconium. Similarly, there is reasonably good agreement concerning the uranium-rich and zirconium-rich ends of the diagram. However, there is no agreement concerning the nature of the diagram at intermediate compositions and at temperatures below about 650°C. The early work of Saller et al.8 hows an intermediate phase, 6. In a later diagram by Saller and Rough' shown in Fig. 1, the d phase field is drawn with dotted lines, indicating doubt that this phase exists. The diagram of Summers-Smith2 does not show the d phase. More recently Holden and Seymour' have shown that in the composition range 40 to 60 wt pet Zr, the high-temperature ? phase, which is body-cen-tered-cubic, transforms to an ordered structure, 8. They found that the order-disorder transformation occurs between 595° and 610°C, and presented a diagram based on the study of four different alloy compositions. Most recently Saller et a1.5 presented a new diagram which differs with respect to the composition limits of the 8 field from the others proposed by these authors. The purpose of the present investigation was to determine the boundaries of the d phase field. Experimental Procedures and Results Alloys—Buttons, weighing 80 g of 15 different alloy compositions were prepared by arc melting Mal-linckrodt uranium and iodide-process zirconium in a water-cooled crucible. A controlled atmosphere of argon and helium was maintained in the furnace. Prior to making each button, a sacrificial button of zirconium was melted to increase the purity of the furnace atmosphere. Each alloy was then melted six times to improve its homogeneity. Subsequently the alloy buttons were jacketed in copper and hot rolled to slabs at 775°C. From these slabs, specimens were machined for chemical, thermal, and metallographic analyses, Table I gives both the estimated chemical composition of the alloys, and the values which were obtained by analysis. The oxygen contents varied from 0.013 to 0.032 pet, nitrogen from 0.0014 to 0.014 pet, and carbon from 0.0013 to 0.011 pet. The purity of the alloys ranged from 99.947 to 99.975 pet. Thermal Analysis—The specimens for thermal analysis were blocks 11/4 x 1/4 x % in. The bead of a calibrated thermocouple was placed between two of these blocks and molybdenum wire was wrapped around the assembly. Each assembly was sealed in a quartz capsule having an argon atmosphere. The capsule was placed inside a refractory crucible which was heated in a furnace. Using a method described by Smith,' the emf of a differential thermocouple, with junctions inside and outside of the refractory crucible, was used to activate the

Jan 1, 1959

By D. Havlena, A. S. Odeh

The use of the straight-line method of solving the material balance equation is illustrated by means of six field cases. Also, the application of statistical criteria to arrive at the most probable answer is shown. The theory underlying the straight-line method of solution and the applicability of the statistical criteria was presented in a previous paper.' The field cases include saturated and undersaturated oil reservoirs with and without water drive. The aquifers discussed are: limited radial, infinite radial, very small aquifer and infinite linear. The field cases also include a gas reservoir producing under water drive. INTRODUCTION In a previous paper,' the authors presented the theory underlying the solution of the material balance as an equation of a straight line. The appropriate equations for various material balance cases as well as the methods of analysis and interpretation with comments and discussion were also included. To illustrate the various theoretical cases treated previously, five field cases are analyzed in this paper by employing the straight-line method of solving the material balance equation (MBE) and one example previously published is referred to. The use of statistical criteria to arrive at the most probable answer is also shown. All the field examples presented, except Case 2, are excerpts from complete reservoir studies. To illustrate the method, only sections specifically dealing with the material balance principles are included. Additional geologic information and basic data are reported to better acquire an understanding of the cases and thus to better follow the reasoning that suggested the successful application of the straight-line method of solving the MBE. The six cases are: (1) saturated reservoir, small gas cap, limited aquifer; (2) saturated reservoir, very small gas cap, infinite aquifer; (3) undersaturated-saturated reservoir, very small aquifer; (4) highly undersaturated reservoir, no water drive; (5) high undersaturated one-well reservoir, limited aquifer; and (6) gas reservoir, infinite linear aquifer. THE D4 SAND, GUICO FIELD, VENEZUELA The D, sand, which was discovered in 1943, is presently in a depleted state. Since its discovery it has produced under water drive, gas-cap-gas expansion, and solution gas drive. In Nov., 1947, water injection was initiated to arrest further pressure decline. When discovered, the D, sand was a saturated reservoir with a gas cap /oil zone volume ratio m estimated volu-metrically at 0.0731, an average permeability of 500 md, a porosity value of 25 per cent, and an oil viscosity at reservoir conditions of 0.3 cp. The volumetrically determined stock-tank oil initially in place was 23.1 million bbl. The volumetrically weighted physical data and production data available until Nov., 1953 are reported in Table 1. In Ref. 1, the effects on the straight-line plot of various values of ro/r for a constant ?to, or of various dimen-sionless times for a constant r8/r10 were theorized and were illustrated in Fig. 3A of that reference. In this field case, the previously theoretically predicted effects are established. Thus, the MBE calculations using Eq. 3c of Ref. 1 were performed for various r6/r10 and dimension-less time values. Eq. 3c of Ref. 1 is:

Jan 1, 1965

By B. S. Gottfried

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By Frank Woolley, John F. Elliott

The high-temperatzrre solulion calorimeter has been modified and an extension dynamic analysis of i/s transient behavior has provided an improved basis for interpretation of the experimental data. The partial molar heats of mixing, Him, of aluminum, copper, and silicon in the corresponding binary liquid-iron alloys haz,e been measured directly at 1600°C. The resulls in kcal per g-alom are: THE high-temperature1,2 solution calorimeter was planned and designed to obtain directly the partial molar heat of solution of solute elements in a variety of liquid metals that might serve as solvents. Results obtained in the temperature range of 1000° to 1200°C were described in the two earlier publications. However, a long-term goal in the program has been to obtain measurements at 1600°C of the heats of solution of elements in liquid iron. Hence, subsequent to the work last reported,2 and based on the experience of the first two studies, the calorimeter system has been modified and operated at 1600°C to obtain measurements on the Fe-Si, Fe-A1, and Fe-Cu systems. This paper reports the results of the work. THE MODIFIED SYSTEM It was clear from the earlier results that the most precise measurements can be obtained with the high-temperature solution calorimeter when the solution time is less than approximately 30 sec. This can be achieved only when the bath temperature, Tb, is above the melting (liquidus) temperature of the solute material being dropped. A second point of importance is that the refractories used earlier liberated considerable water so that, although the oxygen potential of the atmosphere in the system was satisfactory for work with nickel and copper, it was too high for work with iron, silicon, and so forth. Third, for work at 1550" to 1650°C, an improved temperature-sensing unit is needed. The modified calorimeter in its latest form is shown schematically in Fig. 1. The refractory brick and the Kanthal-Super heating elements which formed the furnace for heating the solute addition have been replaced by a smaller resistance-heated tube furnace which is insulated with bubble-grain alumina. A more elaborate gas purification system was installed to maintain a satisfactorily low oxygen potential in the atmosphere of the calorinieter chamber and of the solute addition furnace. As was the case in the earlier work, at the beginning of a run the calorimeter chamber is evacuated and then filled with purified argon so that a total pressure slightly in excess of 1 atm is maintained during the run. When either of the furnaces within the enclosure attains a red heat, 6 pct H2 is added to the atmosphere. Throughout the run, gas from the calorimeter is circulated through a de-oxidation system at the rate of approximately 5 standard liters per min; the changeover time for the gas within the calorimeter enclosure is approximately 1 hr. The essential stages of the deoxidation system are: a copper furnace at 600°C, a drying tube of anhydrite, and a titanium furnace at 800°C. The oxygen potential of the gas stream in the purification system can be monitored at several points by a stabilized zirconia oxygen cell. This cell is invaluable in detecting the development of leaks in the container vessel and the gas circulating system. The installation of the solute-addition furnace allowed either of two types of temperature-sensing units to be mounted on a water-cooled probe. The first was a small optical device using a silicon diode that was

Jan 1, 1968

By Frederic L. Kadey

Diatomite is a siliceous, sedimentary rock consisting principally of the fossilized skeletal remains of the diatom, a unicellular aquatic plant related to the algae. Thus, it has been formed by the induration of diatomaceous ooze, and consists mainly of diatomaceous silica, a form or variety of opal which is first formed in the cell walls of the living diatom. Diatomaceous silica is not generally regarded as a synonym or the equivalent for diatomite, although it has been so used at various times. Accurately, diatomaceous silica is the preferred name for the principal mineral component of which the rock, diatomite, is composed. The terms diatomaceous earth and kieselguhr are used as synonymous with diatomite. The designations tripoli, tripolite, infusorial earth, etc., were used at one time but are now obsolete. With the changing nomenclature, these terms that were at one time correct when proposed and used for generations would be considered incorrect if used today in the light of current knowledge. The designation diatomite is reserved for those accumulations of diatomaceous silica that are of sufficient quality, size, and minability to be considered of potential commercial value. Processed diatomite possesses an unusual particulate structure and chemical stability that lends itself to applications not filled by any other form of silica. Foremost among these applications is its use as a filter aid, which accounts for over half of its current consumption. Its unique diatom structure, low bulk density, high absorptive capacity, high surface area, and relatively low abrasion are attributes responsible for its utility as a functional filler and as an extender in paint, paper, rubber, and in plastics; and as an anti-caking agent; thermal insulating material; catalyst carrier; and chromatographic support; polish, abrasive, and pesticide extender to name a few representative applications. The United States is the principal producing country, although diatomite is found in numerous other locations. Geology Composition and Morphology Diatomaceous silica qualifies as a. mineral of organic origin in much the same way that aragonite and collophane do. The silica of the fossilized diatom skeleton closely resembles opal or hydrous silica in composition (SiO2 • nH2O) .4, 12, 203 The silica is of acute biological significance, not only for the cell wall component, but also for the basic life process.45, 51 Without silica, cell development ceases .31 In addition to bound water, varying between 3.5 and 8%, the siliceous skeleton may also contain, in solid solution, or as part of the Si02 complex, small amounts of associated inorganic components-alumina, principally-and lesser amounts of iron, 46, 49, 55 alkaline earths, alkali metals, and other minor constituents.59, 63, 67 Boron is reported to be an essential element for diatom growth .52, 58 Since diatomaceous silica is not pure hydrous silica but contains other intimately associated elements, there is good reason to consider it a distinct type or variety.4 Associated with the diatomaceous silica, and integrated as part of the diatomite, may be variable amounts of organic matter, soluble salts, and particles of rock-forming minerals that were syngenetically deposited or precipitated with the diatom frustules. Sand, clay, carbonate, and volcanic ash are typical common contaminants. Other contaminating minerals may be present, such as feldspar, mica, amphiboles, pyroxenes, rutile, zircon-the result of weathering, then transporting, and subsequent redeposition of surrounding land masses. Commercial diatomite may also contain fragments and particles of

Jan 1, 1983

By John H. Schemel

Large Zircaloy-2 hammer or press forged bars did not exhibit the uniform excellent corrosion resistance to steam normally expected of the alloy in wrought form. Weight gains of coupons cut from forged bars were 40 to 60 mg per sq dm compared to the 28 mg per sq dm obtained on hot rolled sheet and strip. The mechanism of this corrosion and its relation to microstructure is discussed. After these observations, a solution heat treatment followed by rapid cooling was tried on coupons from eight heats of forged bars. The basket-weave structure resulting from the quench did not show agglomeration of intermetallic compounds. These coupons showed uniformly good corrosion resistance with weight gains very close to the 28 mg per sq dm expected of Zircaloy-2 after the 14-day steam test. A large forging was solution heat-treated and tested for structure, mechanical properties, and corrosion resistance. Corrosion resistance was excellent and uniform throughout the section. The normally anisotropic mechanical properties were changed to completely iso-tropic. Strength levels were raised while there was a small loss in ductility at the service temperature. Tempering of the quenched structure below the ß transus did not improve the ductility at the test temperature of 600°F. ZIRCALOY-2 end caps for nuclear fuel elements Zircaloy-2 is an alloy of 1.5 pct Sn, 0.1 pct Fe, 0.1 pct Cr, and 0.05 pct Ni, balance hafnium-free zirconium patented by Westing-house Electric Co. used in pressurized water reactors are machined from large hammer-finished forged bars. The bars forged for this application range from 3 1/2 in. to 7 1/2 in. squares. One of the tests used to evaluate Zircaloy-2 for this service is a 14-day exposure to very high-purity steam at 75oF and 1500 psi pressure. Normally wrought products will exhibit a black lustrous film of zirconium oxide after the test and show a weight gain of approximately 28 mg per sq dm. Coupons representing the forged bars exhibited a wide variety of corrosion results ranging from acceptable black coupons to some covered with white crystals of zirconium oxide that had weight gains of nearly 100 mg per sq dm. One of the most common effects was a white corrosion product that looked like a stain to some observers and the outline of massive metal grains to others. Very careful specimen preparation prior to the corrosion test did not affect the result and metallographic examination did not reveal a structural feature of a size and shape that would correlate with the corrosion product. In addition to the relatively poor corrosion resistance, the forged product was not as ductile as desired. The problem of obtaining uniform and more acceptable properties in these heavy-forged sections seemed to be a function of the microstructure of the metal and, in particular, the distribution of the intermetallic compounds. None of the four alloying elements however are appreciably soluble in a zirconium and are present as compounds which usually appear scattered randomly throughout the structure. Moudryl showed that the distribution of these intermetallic compounds in Zircaloy could be related to a ostringero corrosion failure. Grozier et al. discussed another structural defect that causes a similar corrosion effect. In this case, an elongated gas-void was determined to be the cause. M. L. Picklesimer held in his discussion of the Grozier paper that at least a part of the observed =stringersn were caused by the distribution of the intermetallic compounds. He proposed a

Jan 1, 1962

By L. R. Duncan, O. M. Wicken

Magnesium, the eighth most abundant element in the earth's crust, is found widely distributed in a variety of minerals. Among the more commercially important ones are magnesite (MgCO3), brucite (Mg[OH]2), dolomite (CaMg[C03]2), and the salts magnesium sulfate and chloride often found in natural brines. Among other commercially important magnesium-containing minerals are olivine ([MgFe]2Si04), talc (H2Mg3[Si03]4), and serpentine (H4Mg3Si2O9). These are valuable because of desirable characteristics given them by the magnesium content and its placement in the particular crystal structure. These minerals are the raw materials for a host of products including magnesium metal and several grades of magnesia or magnesia-containing materials for refractories, fluxes, fillers, insulation, cements, decolorants, fertilizers, and chemicals. Considering only those minerals valued for their magnesia or magnesium content, over 85% of the material tonnages are processed to be used as refractories. Dolomite, which is used in enormous tonnages for its physical properties in construction, has an important use as a refractory raw material in its own right and, after calcination, as a precipitant for magnesium hydroxide, which is an intermediate for the production of magnesia or magnesium metal. Magnesia from the reaction of calcined dolomite (or limestone) with seawater or magnesium-containing brines has supplanted much of the production of magnesia from mined magnesite, particularly in the United States. The best known of the minerals directly and widely exploited for magnesia content is magnesite, one of the calcite group of rhombohedral carbonates which includes calcite (CaC03), siderite (FeCO3), and rhodochrosite (MnCO3) among others. The members of this group enter into a wide range of substitutional solid solutions when the positive ions have similar radii. The radii of magnesium and iron ions are within 6% of each other; hence, magnesite and siderite form a complete series of which breunnerite (ferroan magnesite) is a well-known end member. However, the radius of calcium ion is 36% larger than that of magnesium ion, and only limited substitution exists at each end of the MgCO3-CaCO3 series. Dolomite is not a member of the calcite group, but results when calcium and magnesium ions alternate in equal number in an ordered structure among carbonate ions. The result of these relationships is that calcite and dolomite, often found intermixed with magnesite, occur commonly as identifiable crystal entities which can often be separated to a varying degree from the magnesite by various beneficiation techniques. Magnesite, when pure, contains 47.8% MgO and 52.2% CO.. The pure mineral is sometimes, but rarely, found as transparent crystals resembling calcite, but, preponderantly, magnesite contains variable amounts of the carbonates, oxides, and silicates of iron, calcium, manganese, and aluminum. Magnesite may be either crystalline or amorphous (cryptocrystalline). The crystalline form has a hardness of 3.5 to 4.0. The color may range from white to black with shades of yellow, blue, red, or gray. The color is not a fundamentally significant indicator of purity, but in a given deposit, an experienced person can often roughly grade magnesite by assessing color and crystallinity. The crystalline form of magnesite occurs in

Jan 1, 1983

By R. F. Domagala, D. J. McPherson

Iodide zirconium was combined with calculated amounts of ZrO2 or master alloys and arc-melted. Annealing treatments were carried out at 21 temperature levels. Metallographic examination of the heat treated specimens permitted construction of the binary phase diagram from zirconium to ZrO2. Oxygen additions to zirconium raise the transformation temperature as well as the melting point. Features of the diagram include the peritectic formation of beta, the formation of alpha directly from the melt, an intermediate phase ZrOB with a' range of homogeneity, and a eutectic between alpha and Zr02. PHASE relationships in the Zr-0 system were determined in the range 0 to 66. 7 atomic pct 0 (ZrO,). Arc-melted alloys were annealed at temperature levels between 600" and 2000°C. Determination of the phase boundaries was accomplished by metallographic evaluation of specimens quenched from the various temperatures. Incipient melting techniques were used to determine solidus curves, and X-ray diffraction work was employed to study the lattice parameters of the a solid solution and the phase ZrO,. Experimental Procedures Westinghouse "Grade 1" iodide zirconium crystal bar (approximately 99.8 pct pure) was employed for the investigation. This material is substantially free of the impurity inclusions characteristic of most zirconium metal. Previous work has shown it to be the most satisfactory grade for phase diagram studies. The zirconium crystal bar, as-received, was coated with corrosion product from autoclave tests by which its grade designation is determined. The bars were sand blasted lightly, pickled for 1 min in a 20 pct HN0,-5 pct HF aqueous solution, rinsed in water and acetone, and dried. The bars were rolled to about 1/32 in. strip, cut into 10 in. lengths, and pickled again. The material was sheared to approximately Y4 in. squares, cleaned with acetone, and stored for use. Two pounds of specially prepared "chemically pure," and 100 grams of spectrographic grade ZrO, were obtained for the preparation of alloys. The "Specpure" ZrO, was obtained from Johnson, Mat-they Co.. Ltd.. and Titanium Alloy Manufacturing Div. of National Lead Co. supplied the other grade of oxide. The "chemically pure" material (henceforth referred to as C.P. ZrO,) was used for the preparation of all alloys employed in the determination of phase relationships. A limited number of alloys were prepared with the "Specpure" oxide to recheck certain phase boundaries. Analyses of the two grades of oxide are given in Table I. The C.P. ZrOI was received in powder form, not suitable for arc melting. Therefore, this material was compressed to wafers 1 in. in diameter and about Vn in. thick. These were segmented and stored for use. The "Specpure" material, in the form of small granules, was suitable for use in the as-received condition. A nonconsumable electrode arc-melting furnace was used to prepare alloys. A water-cooled copper crucible and tungsten-tipped electrode were employed. The material was placed in the crucible and rapidly melted under a protective atmosphere of helium gas by striking the arc on a tungsten stud and directing it upon the charge material. Details of operation, including drawings of this type furnace, have been published.'

Jan 1, 1955

By Robert B. Hoy

If the reader finds that the basic information in a commodity chapter is insufficient, he can consult the appropriate sources in this chapter to find more detailed or more up-to-date information. In these days, technology develops so rapidly that this 5th edition will be somewhat out of date even by the time it is distributed. A comparison of this edition with the 4th edition (1975) of Industrial Minerals and Rocks will clearly demonstate how rapidly the scene changes. One of the best sources of new information, of course, is attendance at technical meetings. For example, the papers presented at the national and sectional meetings of AIME usually discuss new developments in the respective divisions of the society. Frequently this information is available in Preprints. In addition some are printed in Mining Engineering and/or the Transactions of the Society of Mining Engineers of AIME. Furthermore, other professional and trade journals contain articles pertinent to the field. A journal introduced in 1967, Industrial Minerals, published by the Metal Bulletin, Ltd., London, is an excellent source of information. The Mining Journal, Ltd., also in London, publishes a weekly journal, a Mining Magazine (monthly), and a Mining Annual Review. The latter covers the field from three aspects: (1) by commodities, (2) through technical developments, and (3) geographically. The Industrial Minerals Directory 1977, 1st ed., Metal Bulletin Books, Ltd., London, provides addresses of organizations in the field. Other useful magazines and periodicals include Engineering and Mining Journal, Pit and Quarry, Rock Products, and World Mining, all published in the United States. In addition many new books related to industrial minerals have been published since the 4th edition of Industrial Minerals and Rocks. Most of these are included in the bibliography of each chapter in this volume. Another excellent information source is the proceedings of the Forum on Geology of Industrial Minerals. This group, which held its 18th annual meeting in 1982, is sponsored by a state geological survey or another mineral agency on a rotational basis. The proceedings have been published by the various sponsors and are available as indicated in Table 1. Two major sources of basic industrial mineral and rock data in the United States, and also a good source for many foreign countries, are the US Geological Survey and the US Bureau of Mines. The two agencies publish many reports and maps on a variety of geological and mining subjects. In addition, each agency has an indexed list of all their publications with a free service announcing new publications: New Publications of the Geological Survey, US Geological Survey, 329 National Center, Reston, VA 22092 and New Publications-Bureau of Mines, Publications Distribution Branch, US Bureau of Mines, 4800 Forbes Ave., Pittsburgh, PA 15213. The US Bureau of Mines (USBM) publishes annually Minerals Yearbooks which provide statistical data on mineral products including production, consumption, imports and exports, as well as information on new developments. Another very useful USBM series is called Mineral Facts and Problems. Editions were published in 1960, 1965, 1970, and 1975. Four major series of publications by the Bureau are Information Circulars, Reports of Investigations, Mineral Trade Notes, and Minerals and Materials. The last two are particularly good for news of foreign developments. The Bureau of Mines Employee Directory, available at all USBM offices, lists addresses and telephone numbers for all offices as well as same for all commodity and area specialists.

Jan 1, 1983

By R. W. Heckel

The growth and shrinkage rates of second-phase particles of various size distributions are analyzed in terms of five mathematical models. These models apply to various forms of diffusion-controlled growth (one model) and interjace-controlled growth (four models). The expression for diffusion-controlled growth assumes local equilibrium at the interface between the second Phase and the matrix and through the choice of the proper flux area perrnits the definition of minimum diffusion-controlled growth kinetics. The interface-controlled kinetics models consider the situations where the rate-controlling step occurs at either the interlace undergoing solution or the interface undergoing deposition. In addition, each of these interface-controlled kinetics models is analyzed for interface reaction rates which are independent of the thermodynamic activity gradient across the interface or proportional to the gradient. THE variation of the size distribution of second-phase particles in a solute-saturated matrix has been a stimulating topic for many years. Research has been carried out in this area for situations where the matrix was either liquid or solid. It is well-understood that the over-all process results in an increase in the average particle size, a decrease in the total number of particles per unit volume, and a net decrease in surface area. Greenwood1 has presented a mathematical model based upon the mass balance existing between growing and shrinking particles. His analysis is stated in terms of Fick's First Law using concentration gradients set up between particles of different size as a result of the effect of surface curvature on solubility. Greenwood's analysis indicates that the mean particle size should increase approximately as the cube root of time. Asimow2 has shown that several of the simplifying assumptions employed by Greenwood may be eliminated from this diffusion-controlled model with little increase in complexity. Lifshits and Slyozov3-5 have also considered the growth-shrinkage problem involving diffusion-controlled kinetics. Their analysis was extended to include the effects of interacting particles ("kinetic approximation"), particle shape, and elastic strains. These authors also show that the mean particle size varies as the cube root of the reaction time. wagner6 has formulated expressions for both diffusion-controlled kinetics and interface-controlled kinetics in the growth-shrinkage problem. Wagner's diffusion-controlled kinetics is based upon the same form of mass balance using Fick's First Law as used by others and is shown to be applicable to situations where kr >> D (k is the interface reaction constant, r is the particle radius, and D is the diffusion coefficient). This analysis also yields a mean particle size which varies as the cube root of time, in agreement with other workers. Wagner also developed interface-controlled kinetics to handle kr « D. This model was based upon the fact that the interface reaction rate is proportional to the interface area. Wagner assumed that the rate was also proportional to the thermodynamic activity difference, ?a, existing across the interface. Thus, the specific reaction rate at the surface of radius r was defined as k • 4nr2 . ?a. It is the purpose of the present paper to reconsider the processes of diffusion-controlled growth and interface-controlled growth in the variation of particle-size distributions as a function of time at constant second-phase volume. The analyses are formulated in a manner which permits data obtained

Jan 1, 1965

By Gareth Thomas

The analysis of Kikuchi pattersns of exct ovientalions from single cryslals and paired Kikuchi lines from single and overlapping crystals is shown to be useful and quanlitalve and is applied to Phase transfovnzcitions including ordering, spinodals, mavten-silic, and nuclealion and growth pyocesses. 112 pinciple, the analysis of exacl orientations enables the crystal system and the Bravais lattice of a crystal to he determined. The advantages of the davk-field imaging technigure for detecting vevy small precipitates are also described. ALTHOUGH Kikuchi electron-diffraction patterns were first observed nearly 40 years ago,' little systematic application seems to have been made of them, until fairly recently during electron microscopy, when their usefulness in contrast experimen and for determining exact orientations8'9 has been pointed out. With the availability of gonio-metric specimen tilting stages it has now become possible to make much wider use of diffraction patterns, particularly the Kikuchi pattern, which is the main subject of this paper. The treatment presented here is not exhaustive but it is hoped that it will stimulate more use of Kikuchi patterns during electron-microscopy investigations. The Kikuchi pattern is formed as a result of Bragg diffraction of the inelastically scattered electrons produced during the interaction of the beam and thick specimens, The important feature of these patterns is that they give an accurate representation of the symmetry of the crystal being investigated, so that it is possible to identify crystal systems and even the Bravais lattice. This means that new structures, e.g., formed in phase transformations, may be identified during normal electron microscopy, so that Kikuchi-pattern analyses considerably extend the uses of the electron microscope. Recent work has also shown that dark-field images are more informative than bright-field images, particularly, for observing small precipitates. The second part of this paper discusses some applications and advantages of dark-field imaging in studies of two-phase systems. 1) DIFFRACTION PATTERNS 1.1) Spot Patterns. Electron-diffraction spot patterns have their limitations because of the importance of the form factor on the intensities and shape of the reciprocal lattice points. Because of the extension of these points into relrods, reflections are possible over a large angular range (+5 deg) and the patterns from thin regions can be complicated because second-layer relrods intersect the reflecting sphere. Spot patterns can thus give only an approximate idea of the crystallography unless the foil is tilted into exact orientation. In this case the spot pattern is symmetrical, with equal numbers of spots on the positive and negative zone directions about the origin. Such cases are necessary for structural analyses and have been used recently to determine the crystal structure of the ordered Ta64C phase." Exact orientation means that the plane of the reciprocal lattice lies exactly normal to the incident beam as shown in Fig. 1(b). It should be noted that in exact orientations the angle of reflection is less than the exact Bragg angle 9 so that the reciprocal lattice points lie to the outside of the sphere, Fig. 1(b). This deviation is denoted by the parameter s and in the usual convention4 s is negative for exact orientations and, of course, zero at the exact reflecting position shown in Fig. l(a). In order to avoid secondary reflections from thickness relrods it is advantageous to work in thicker regions of foils. Diffraction from thick regions may also produce Kikuchi patterns and as dis-

Jan 1, 1965

By Langtry E. Lynd

The predominant tonnage use of titanium is as a white pigment, in the form of rutile or anatase, which are allotropic modifications of titanium dioxide. Because of its whiteness, high refractive index and resulting light scattering ability, titanium dioxide is unequaled as a white pigment for paints, paper, rubber, and a wide variety of other materials. Consumption of titanium dioxide pigments for these uses has been rising steadily for many years, United States consumption in 1956 amounting to about 486,000 tons, worth $250,000,000. Production was off a little in 1957 and 1958. An outstanding development in the titanium field in recent years is the expansion of titanium metal production. Virtually nonexistent before 1948, United States production of pure titanium metal sponge reached a high in 1957 of about 17,200 tons, worth about $77,000,000. United States titanium sponge metal production in 1958 was cut back to about 4,600 tons because of reductions in military requirements. Another major use of titanium, in the form of the mineral, rutile, is for welding rod coatings. United States consumption of rutile for this purpose in 1958 amounted to about 11,000 tons, worth $1,100,000. Titanium-Bearing Minerals Although several minerals contain titanium as a major constituent, only a few of these occur in sufficient quantities and concentrations to be of commercial importance. These are ilmenite (FeTiO3), rutile (TiO2), and alteration products of ilmenite which have been referred to as weathered ilmenite, leucoxene, or arizonite. The other polymorphs of titanium dioxide, brookite and anatase, occur in a few deposits, of which only one has been of even minor importance. Considerable concentrations of titanite or sphene (CaTiSiO5) and perovskite (CaTiO3) are known, but do not compete with ilmenite and rutile deposits in grade or size. Rutile is very widely distributed as an accessory mineral in most rocks, but large deposits of commercial importance are rare. A few rock deposits have been worked for rutile but only one, in Virginia, is in operation (1959). Development of one in Mexico was discontinued recently. At present, the main sources of rutile are the beach and dune sands of Australia and Florida. A major use of rutile is for the manufacture of titanium tetrachloride, which is the intermediate compound universally used in titanium metal manufacture. Another main use of rutile is for welding rod coatings. Although the whiteness and brightness of titanium dioxide pigment is due to the high refractive index of rutile or anatase, natural rutile is too impure to be used directly as a white pigment, is too expensive to use for manufacturing pure oxide pigment, and is too scarce to provide the required tonnages. Instead, the cheaper, more abundant, mixed oxide of iron and titanium, ilmenite, which is readily decomposed by sulfuric acid, is used for manufacture of rutile and anatase pigments. Ilmenite is mined from both rock and sand deposits. Rock deposits contain ilmenite in massive or heavily disseminated form along with magnetite or hematite and are almost invariably associated with anorthositic and gabbroic rocks. Ilmenite also occurs as an accessory mineral in a great many igneous and metamorphic rock types. Concentrations of

Jan 1, 1960

By J. F. Breedis

The morphology and crystallography of marten -site formed during quenclzing were examined by transmission electron microscopy in alloys whose compositions lie between Fe-19 wt pct Cr-11 wt pet Ni and Fe-33 wt pet Ni. Depending upon composition, four types of transformation products were observed in this portion of the Fe-Cr-Ni system. Three of these products are denoted as lath, plate, and surface martensite. The hcp structure found in association with martensite it? certain higlz-clzromium stairzless steels does not influence the transformation crystallogyaphy and cannot he Considered to act as a transitional structure in the transformation from austenite to martensite. DISTINCT changes in the characteristics of mar-tensitic transformations are produced by the partial replacement of nickel by chromium in austenitic iron alloys. Large, discrete martensite plates form in Fe-30 pct Ni, while aggregates of small crystals form in sheets parallel with (111)A austenite planes in stainless steels containing in the vicinity of 18 pct Cr and 8 pct Ni. These transformation structures are termed plate and lath martensite, respectively. In addition, an hcp structure (E) has been observed in certain high-chromium stainless steels in association with martensite. It has been suggested that the E structure acts as an intermediate step in the nucleation of martensitel-' and influences the morphology and crystallography of the final transformation product. However, this structure may be the consequence of the deformation of austenite arising from the shape distortion accompanying the transformation to bcc martensite.' An hcp stacking of atoms is produced in an fcc lattice by the introduction of an intrinsic stacking fault on alternate (111)A planes. At present, the arguments behind each explanation for the E structure are based on circumstantial evidence. This investigation seeks to determine the importance of the hcp structure to the martensitic transformation, and to examine the crystallography and morphology of the structures produced in bulk samples by quenching in alloys whose compositions lie between Fe-19 wt pct Cr-11 wt pct Ni and Fe-33 wt pct Ni. These aspects of the martensitic transformations occurring in these alloys have been studied by electron metallography and selected-area diffraction of foils examined in transmission in the electron microscope. EXPERIMENTAL METHODS The alloys have been prepared from thoroughly mixed powders of pure iron (99.6+ pct, 0.06 max pct C, 0.15 max pct N), chromium (99.9+ pct), and nickel (99.9+ pct). After sintering in an Ar-H2 atmosphere, the mixtures were melted by induction heating in recrystallized alumina boats under purified argon. The ingots were cold-rolled to 0.4-mm sheets (100 pct reduction in thickness) after the removal of at least 1.5 mm from their surfaces by grinding. Three-millimeter discs were punched from the cold-rolled sheet and annealed at 1150°C for approximately 24 hr in evacuated quartz capsules. The average grain size was 40 µ Thin foils for examination in transmission in the electron microscope were prepared from discs quenched to liquid-nitrogen temperature. A depression was first jet-machined into the center of a disc and a thin area (less than approximately 3000A in thickness) was produced at the base of the depression by electropolishing. Electron-microscope studies utilized the RCA EMU-3 and the Philips EM-100B microscopes operating at 100 kv. X-ray examination of the polycrystalline alloys was performed either in a diffractometer or in a

Jan 1, 1964

By F. D. Wright

Production of a coal preparation plant is governed by many restrictions, such as the tonnage of different products and blends that can be sold within a given period, capacities and output proportions of the cleaning and sizing units, blending proportions, quality specifications, and costs and prices for the various products. Determination of the tonnage of each product and blend that should be made in order to obtain the maximum profit can be difficult unless a systematic method such as linear programming is used. In this paper the basic method of linear programming is described briefly. The Old Ben No. 9 preparation plant is used as an example to illustrate in detail how eqwations can be written to form a linear programming model of a coal preparation plant. Three sample problems, each requiring 56 or more equations and 63 structural variables, were solved with an IBM 650 computer. Linear programming is one of several mathematical tools used for operations research. It has been applied to many fields and has been used by a number of industries either to maximize the profits of certain operations or to minimize costs. P. B. Nalle and L. W. weeks' have described the use of the method by the Riverside Cement Co. to minimize the cost of blending raw materials to make portland cement. Their problem is to obtain at minimum cost a mix with certain specifications from a number of possible materials which have various costs and various amounts of CaO, SiO2, Fe2O3, and other constituents. In a paper on the use of linear programming by the National Coal Board in England, K. B. Williams and K. B. Halley2 describe how the transportation method, a variation of linear programming, is used to minimize the cost of sending 37 grades of coal from 28 mines to seven central washing plants which produce coal for furnace coke and foundry coke. The various coals have different percentages of volatile matter, moisture, sulfur, ash, and phosphorous, so there is considerable choice in how they can be blended to meet the specifications of the two products. The purpose of this paper is to show how linear programming can be used to maximize the profits of a coal washing plant which produces individual final products as well as blends. The Old Ben Corp. furnished assumed sample data from their Old Ben Mine No. 9 preparation plant for this investigation. However, the data that have been used are entirely the author's responsibility. METHOD OF LINEAR PROGRAMMING Numerous articles and books have been written on the theory and applications of linear programming.3-5 However, since the method has not been widely applied by the mining industry, a brief, nontheo-retical discussion of its basic method seems to be in order. Linear programming (Table I) is used to determine the best possible solution to a number of interdependent activities. It is essentially a method for making systematic selections from a number of possible solutions to determine positively the optimum solution. There may be other equally good solutions but no better ones. Each activity must have linear coefficients and there must be a criterion to judge how good the solution actually is. Linear programming is different from the solution

Jan 1, 1961

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By C. C. Mattax, J. E. Bobek, M. O. Denekas

ABSTRACT investigations in recent years have shown that rock wettability can exert a profound influence on the displacement of oil by water from oil producing reservoirs. Core analyses frequently show oil recoveries from preferentially water-wet rock to be significantly greater than those from preferentially oil-wet rock. Thus, an accurate prediction of water-drive or waterflood oil recovery is dependent on the evaluation of reservoir rock wettability. Laboratory and field studies have been undertaken to examine and develop methods which may be used to measure rock wettahililiy, to investigate factors which may alter the wettability of reservoir rock, and to determine the wettability of specific reservoirs. It has been found that measurement of the rate and volume of spontaneous imbibition of the wet-ting phase by a rock is a reliable and reproducible test for semi-quantitarive determination of preferential rock wettability. Laboratory tests have shown that the wettability of reservoir rock may depend on both the crude oil composition and the rock type. Field and laboratory tests have indicated that coring flulids and core handling techniques can cause significant changes in the wettability of rock surfaces. However, a few fluids, brine in particular, appear not to affect core wettability and may be used when coring to determine reservoir wettability. Core handling and packing procedures have been developed which preserve core wettability during storage and laboratory testing. The wetfability characteristics of six reservoirs have been evaluated by wettabilify tests performed in the field on fresh cores cut with brine or special water-base muds. The tests indicated that all six reservoirs are water wet; however, the de-gree differed from one reservoir to another. INTRODUCTION Although surface properties of reservoir rock and interfacial characteristics of rock, crude oil, and water systems have been discussed for many years, the full significance of these properties to the production of oil is still in the process of evaluation. It is known, however, that one surface characteristic—preferential wettability of the reservoir rock—-is of paramount importance when the reservoir is being produced by water-flood or water-drive mechanisms. The importance of rock wettability has been pointed out by several authors'.',%* who have found that pref- erentially water-wet cores flood more efficiently than oil-wet cores; that is, more oil is recovered from water-wet cores in the early flooding stages than from oil-wet cores. Other authors'," indicate that waterflood oil recovery from cores of intermediate wettability may be greater than that from either strongly water-wet or strongly oil-wet cores. These recent investigations emphasize that the true wettability of a reservoir must be known before the performance of waterflood or water-drive reservoirs can be accurately evaluated. When reservoir wettability is to be determined by core analysis, precautions must be taken to maintain the original wettability which can be altered readily during coring5 nd testing.' In the coring operation the core may be partially or completely penetrated by the drilling fluid which, if it contains surface active materials, may drastically change the wettability of the core. Also, improper core handling during stcrage and testing may change its wettability because of evaporation of fluids or exposure to oxygen or surface active contaminants. One purpose of this paper is to discuss methods for obtaining, preserving, and testing cores without altering their original wettability. In addition, the results of wettability tests on fresh and preserved cores from six reservoirs are presented where special precautions were taken to obtain samples without altering their wettability during the coring operation.

By Mark J. Klein

An internal-friction profile induced by nitrogen in Cr-35 at. pet Re was studied as a function of nitrogen concentration and heat treatment. From these studies, the solubility of nitrogen in this alloy was deduced. The equilibrium solubility of nitrogen in chromium does not appear to be substantially changed by the addition of 35 at. pet Re to chromium. However the precipitation kinetics of nitrogen from solution are significantly decreased. IT has been found that rhenium has a beneficial effect on the low-temperature ductility of Group VI-A metals. As part of a study to investigate the mechanism of this increased ductility, it was of interest to compare the concentration of interstitials retained in solution in one of these ductile bee alloys with the concentration retained in solution in its base metal. In this regard, a study of the solubility of nitrogen in ductile Cr-35 at. pet Re was selected for this phase of the investigation because the solubility of nitrogen in the base metal, chromium, has been determined over a fairly large temperature span.1-3 In addition, nitrogen is known to have an embrittling effect on chromium,4 but does not appear to have a comparable embrittling effect on Cr-35 at. pet Re.5 The experimental method selected for this study was the measurement of internal friction. This technique allows interstitial-concentration measurements to be made in bee metals that are sensitive to small changes in the concentration of a particular interstitial, even when other interstitials are present. When subjected to a cyclic stress, interstitial atoms in solid solution give rise to an internal-friction peak whose origin is the stress-induced ordering of interstitials into preferred lat- tice sites. These internal-friction peaks signify a dissipation of elastic energy which is greatest when the period of oscillation is comparable with the relaxation time for the attainment of an equilibrium interstitial distribution. For dilute interstitial binary alloys, the height of the internal-friction peak is proportional to the concentration of interstitial atoms in solid solutions. snoek8 has shown that the internal friction, Q-1, may be expressed by the equation where A is the relaxation strength, ? is the angular frequency, and 7 is the relaxation time. The relaxation strength, A, is proportional to the concentration of interstitials taking part in the relaxation if no interaction occurs between solute atoms. Thus values of peak height, which can be determined as a function of temperature, give a measure of the relative concentration of interstitial atoms in solution at various temperatures. If equilibrium has been attained, the peak heights will be proportional to the interstitial solubilities. The constant that relates A to the interstitial concentration can be determined experimentally from peak-height measurements as a function of concentration for specimens of constant crystallographic orientation. When this has been done, Q-1 measurements can be used to determine solubilities according to the equation C - KQ-1 where C is the soluble interstitial concentration and K is a constant that can be experimentally determined. It is possible, however, that lattice imperfections may bind interstitials to certain lattice sites. Internal friction does not measure these interstitial~, but does measure interstitials which are free to move in the normal lattice sites when subjected to an applied oscillating stress. On alloying the pure metal with a substitutional

Jan 1, 1965

By C. D. Stahl, S. M. Farouq Ali

This investigation attempts to describe and simulate the alcohol displacement process by means of a cell model, as employed in chemical engineering practice. The proposed model is more simple than previously proposed models. and utilizes parameters chosen on a theoretical basis. The model successfully reproduced the formation of the stabilized bank and the breakthrough of alcohol, the latter depending on one of the model parameters, which may be correlated with the length of the porous medium. Moreover, the effects of the phase behavior of the liquid system involved, as observed in experimental studies, were reproduced. Several variations of the basic model were devised and tested on a digital computer. These included the cases in which: (1) the actual value of fractional flow was used in cell-to-cell computations; (2) the number of cells was varied within the same run; and (3) incomplete rather than complete phase equilibrium was assumed within each cell. The proposed cell model clarifies the basic mechanism of the process. Detailed concentration profiles obtained for each cell, for instance, showed the mechanism of bank formation in relation to the phase behavior characteristics. The results obtained indicated a varying degree of phase equilibrium concommitant with changes in the velocities of the phases in an actual alcohol displacement. This condition was approximated by changing the number of cells during the simulation. Interesting information was obtained on the influence of path length on the efficiency of alcohol displacement, which has been the subject of some controversy. Certain limitations preclude the use of the proposed model as a substitute for experimental studies. The results obtained were, nevertheless, of value in interpreting the experimental results. INTRODUCTION During recent years considerable effort has been directed toward an understanding of alcohol displacement, the process whereby oil and water are recovered from a porous medium by the continuous injection of a solvent. The complex nature of the physical process involved has so far defied a complete mathematical treatment. Other methods of approach, amounting to an overall material balance, have been proposed, yielding useful information on certain aspects of the process.1-3 Taber et al, in particular, defined the displacement mechanism in terms of the phase behavior of the alcohol-oil-brine system invoIved.2 Wachmann reported a mathematical treatment of alcohol displacement subject to certain simplifying assumptions? Donohue proposed the use of a "cell model" for simulating alcohol disPlacement.5 The nature of the assumptions involved limited the utility of the model. The present work attempts to examine the variables involved in the simulation of alcohol displacement and discusses several possible versions of the basic cell model. Under certain conditions the model results are similar to the experimental results. In particular, the spontaneous formation of of the stabilized bank and the effects of the system phase behavior were successfully reproduced. PREVIOUS WORK ON CELL MODELS Cell models and the theoretical plate concept are often used in solving chemical engineering problems in which an explicit mathematical solution may be difficult or impossible to obtain. Examples of such applications occur in distillation, gas-liquid chroma tography,6 reactor technology, absorption, etc. In petroleum engineering, such a model was used by Attra7 to describe non-equilibrium gas drive, and by Higgins and Leigh ton8 to calculate sweep efficiency in water flooding. Aris and Amundsen pointed out the equivalence between the diffusion model and perfectly mixed cells connected in series.9 Deans proposed a three-parameter cell model to simulate two-component

Jan 1, 1966

By G. R. Speich, Morris Cohen

The growth rate of bainite has been determined by hot-stage metallography in five hypereutectoid high-purity iron-carbon, iron-carbon-chromium, and iron-carbon-nickel alloys. The studies have been confined essentially to lower bainite. The growth rate is independent of time at constant temperature, but both the edgewise and sidewise growth rates increase with increasing temperature. Increasing the carbon content, or adding chromium or nickel, decreases the growth rate without greatly affecting the temperature dependence. The results lead to separate diffusion-controlled models for the edgewise and sidewise growth rates. ALTHOUGH bainite has been subject of many detailed investigations,1,2 the growth mechanism of this intermediate transformation product of steel remains unclear. One difficulty has been that practically all the kinetic measurements have been restricted to studies of the overall reaction rate. Only a limited amount of information3, 4-6 is available on the growth process itself. The purpose of the present work was to investigate the effect of time, temperature and composition on the growth rate of bainite using hot-stage metallography. This method makes it possible to observe continuously the isothermal growth of individual bainitic units because of the surface upheavals accompanying their formation.3 This information is needed to test the existing growth models and to obtain a better understanding of the factors controlling the bainitic transformation. EXPERIMENTAL PROCEDURE Alloys—The compositions of the alloys used in the present work are shown in Table I. The five alloys were vacuum-melted and cast, and were prepared from high-purity materials. The three iron-carbon alloys were available from the work of Roberts et al.7 in the form of 3/8-in. diam rods.* To facilitate machining, the rods were ing. The rods were then tempered for 3 hr at 650°C to allow machining. The iron-carbon alloys in Table I were selected for investigating the effect of carbon on the bainitic growth rate. Due to the rapidity of the transformations involved, it was not feasible to measure the growth rates in iron-carbon alloys having carbon contents much less than 0.9 pct C. The iron-carbon-chromium and iron-carbon-nickel alloys were selected to study the influence of chromium and nickel on the growth rate at a carbon level of about 1.1 pct C, since this carbon content was intermediate among the iron-carbon alloys. The concentrations of chromium and nickel were chosen to match isothermal trans formation diagrams in the literature.8'9 Hot-Stage—Since the iron-carbon alloys had such low hardenabilities, it was necessary to quench the hot-stage specimens drastically from the austentiz-ing temperature to the isothermal transformation temperature to prevent the formation of pearlite or proeutectoid constituents. A rapid-quenching metal-lographic hot-stage was built for this purpose. A maximum temperature of about 1100°C and quenching rates up to 500°C per sec could be attained with this equipment. The details of the apparatus are described in a recent publication. After being polished and mounted on the hot-stage, the specimen (0.020 in. thick) was protected from oxidation and decarburization during austenitizing by a vacuum of 10-4 to 10-5 mm. The quench to the reaction temperature was accomplished with a helium blast, and the isothermal run at this temperature was conducted in a helium atmosphere. A fine-wire thermocouple spot-welded to the sample and a high-speed

Jan 1, 1961