Search Documents

It has been established that molybdenum sesquisulphide, not molybdenite, is in equilibrium with molybdenum metal and sulphur vapor in the vicinity of 1100°C. The S2 pressure for this system has been obtained using the Knudsen orifice method and the standard free energy for the dissociation of S2. THIS work was undertaken to help in determining the thermodynamic properties of the molybdenum sulphide which is in equilibrium with metallic molybdenum and sulphur vapor in the temperature range 1025" to 1150°C. From the results of MontoroL and of Parravano and Malquori' molybdenite would be expected to be in equilibrium with molybdenum metal, while the results of Guichard indicate that a lower sulphide must be in equilibrium with the metal. Thus it was concluded that the problem warranted reconsideration, especially in view of the suspiciously low temperature independent term for the standard free-energy function for the reaction MoS² Mo + S,. The reaction was calculated and duly noted as being low by Richardson and Jeffes,' using the data of Parravano and Malquori.' Experimental Method The general method used in this investigation is that of the Knudsen orifice method for the determination of equilibrium pressure. The apparatus and technique are similar to those described for the determination of the vapor pressure of silver.% owever, a closed-end McDanel tube instead of fused silica was used to maintain the vacuum, a molybdenum effusion cell with a replaceable lid was employed instead of a silica cell, and a titanium getter maintained at run temperature was used to assure that no oxidation of the molybdenum or molybdenum sulphide was taking place. All runs were made at a residual pressure less than 10-5 mm Hg. The data are given in Table I. Runs 1 through 13 were obtained using molybdenite, analyzed 98.5 pct, from the Molybdenum Corp. of America. The remaining runs were done with a sample from the Climax Molybdenum Corp. which was their grade No. 3 and analyzed 59.87 pct Mo. The weight loss of the cell during heating up and cooling down was subtracted from the total weight loss of the cell during the run. Also, when the tem- perature was changed, the Knudsen cell was placed in the vacuum chamber and heated at temperature for several hours in order to allow the sulphide and the metal to adjust their compositions to the equilibrium values, to desorb or adsorb gases into or from the molybdenum cell itself, and to desorb any water from the molybdenum sulphide or molybdenum powder inside the cell. Water may have been adsorbed inside the cell because, after each run, the top of the cell was removed and the material in the cell was stirred to insure that a layer of molybdenum metal would not build up and obstruct the passage of sulphur vapor into the gas phase. It was established that the molybdenum cell alone decreased in weight an inappreciable amount when heated at the highest run temperature for long periods of time. There are two questions to be answered before the data in Table I can be used to calculate thermodynamic functions: 1—Which molybdenum sulphide is in equilibrium with molybdenum metal and sulphur vapor? 2—What sulphur species are present in the vapor phase? It can be assumed safely that no molybdenum is vaporized, since its vapor pressure is almost negligible at these temperatures. It will be shown subsequently that no sulphide of molybdenum is volatile. Question 1 was answered in the following way: The especially purified molybdenite, obtained from the suppliers mentioned previously, was charged

Jan 1, 1956

By C. L. Sollenberger, Will Mitchell, Ford F. Miskell

Heat treatment of ore prior to comminution reduces power requirements for grinding, reduces grinding media wear, and improves recovery of iron values from a typical Minnesota magnetic taconite. Test data demonstrating this, as well as an analysis of the economics of commercial application of the technique, are presented. THE taconites in general are hard, tough ores, difficult to grind. Liberation of iron mineral constituents usually is accomplished by grinding the ore through at least 100 mesh, and often it has been found necessary to grind substantially through 325 mesh to achieve satisfactory recovery and grade in the concentration process. Because of the fineness of grind required and the enormous tonnages of material contemplated for treatment in the future, costs resulting from grinding media wear and power consumption, together with capital investment required for comminution, approach astronomical figures. Economy in any one of these elements per ton of material ground could very well reflect a considerable saving to the ferrous industry in yearly costs. With this in mind, the Research Laboratories of Allis-Chalmers Manufacturing Co. have launched a program to discover means of effecting this economy. The initial phase of the work as described here deals with a heat treatment of crushed raw ore, followed by thermal shock in cooling, for the preparation of rod mill feed. Several investigators1-" ave noted an improvement in the grindabilities of ores treated in this manner. One investigator2 subjected low grade iron ores to heat treatment in an electric furnace. Basing his conclusions on screen analyses of crushed products, he observed that treated ore was more easily crushed than the untreated. He reported that observation of the treated and untreated ore through a microscope revealed cracks following the grain boundaries in treated ore, whereas no cracks were present in untreated ore. However, few if any have quantified the improvement in terms of total hp hours saved, reduction of wear of grinding media, or reduction of capital costs of grinding equipment involved. By means of Bond rod mill and ball mill grindabilitiesh nd by comparative wear tests, the conditions maintained during this investigation have been definitely evaluated. Heat treatment experiments have been made on a batch scale in a muffle furnace, followed by continuous scale experiments in a rotary kiln under various conditions of temperature, atmosphere, retention time, and quenching to determine the combination that would give the greatest improvement in the grindability of a taconite without affecting adversely the magnetic susceptibility of the magnetite. From these data, relative costs have been calculated for grinding both untreated and heat- treated ore on a basis of plant capacity of 120 tons per hour. Heat requirements for a rotary kiln, as well as the kiln capacity required for the treatment, have been estimated. For these tests 50 tons of ore were obtained from the Aurora mine of the Erie Mining Co., Hibbing, Minn. The ore was reported7 to have been selected from the magnetic portion of the lower chert horizons. An average chemical analysis of' this portion is shown in Table I and an approximate mineralogi-cal analysis is shown in Table 11. No attempt was made to differentiate between the various silicates because the compositions vary widely on the different horizons and depend on the degree of oxidation. Minnesotite, however, is the predominant silicate. The siderite grains were dispersed throughout the ore but were essentially associated with the rock-forming minerals. In a series of preliminary tests in which batches of —6 mesh ore were heat-treated in a muffle fur-

Jan 1, 1953

By D. W. Fuerstenau

Recently R. J. Charles1 showed that comminution of brittle or semi-brittle materials in batch operations is described more appropriately by a variable energy relationship than by the specific relationships proposed by Rittinger,2 Kick,3 and Bond.' That is, the value of n in the empirical differential operation which applies to crushing and grinding (Eq. 1) is variable. Independently, Holmes' developed an identical generalized equaton for size reduction. The experimental data of Charles and Holmes show that the value of n is determined by the material and the comminuting machine and is constant for a given combination of these two conditions. The general objective of the present article is to extend the Charles-Holmes concepts to continuous ball milling. The specific objective is to analyze the operation of a vibratory ball mill dry-grinding chromite in a continuous open circuit and to develop an equation for the comminution of chromite in this mill. Basic Principles: In a batch mill the length of time during which the material is ground is a direct measure of the energy expended on the material, except possibly for exceedingly fine grinding. However, in a continuous operation, the fraction of the total expended energy which is useful for grinding will vary according to feed rate and amount of material in the mill. If the amount of material within a continuously operating vibratory mill is low enough so that only part of the expended energy is used for grinding, the extent of size reduction at constant power input must depend upon the length of time the material is retained within the mill, namely, upon retention time. The following analysis pertains to a continuous mill operating under these conditions. The empirical differential equation that applies to the comminution of brittle materials is:" dE = -Cdx/x" [1] where dE is an infinitesimal energy expenditure, C a constant, dx an infinitesimal size change, x the particle size, and n a constant. For continuous grinding in an open circuit at constant power input, P, Eq. 1 must be written as (dtr)P = -C dx/x" [la] where dt, is an infinitesimal time during which a particle of size x is in the mill, and C' is a constant. A ground product never consists of particles of single size but always comprises a distribution of sizes. In comminution, the expenditure of energy causes a lateral shift in the size distribution.', ' " The general objective of any engineering study of comminution is to correlate the magnitude of this size distribution shift with the energy expended, the ultimate aim being to cause the greatest shift with the least amount of expended energy. Although a number of criteria can be used to measure this shift, Schuhmann's size modulus has been found convenient.' Mathematically, the Schuhmann size distribution can be expressed as follows:" y = l00(x/k) [2] or in differential form: dy = [F- xa-' dx ] K, where y is percentage of material finer than size x, a is a constant that determines weight distribution among particles of various sizes, and k the size modulus that denotes the theoretically maximum-sized particle in any size distribution. By considering the energy required to reduce an element of weight of material, dy, from size x, to size x, the element of time during which the material must remain in the mill at at power will be (dt..).= J'-Cax!(dy) [3] Summing up all the elements of weight from zero to 100 for a product that conforms to the Schuhmann size distribution yields C C f Cdx,\( 100«x°-'dx\

Jan 1, 1960

By J. W. Evans

ONE of the most important and frequent calculations that the extractive metallurgist is called upon to make is that of the standard free energy change of a reaction (?F°). For many reactions of metallurgical interest, the lengthy arithmetical calculations have been eliminated by the use of free energy-temperature diagrams which have been constructed by Ellingham,1 Richardson and Jeffes,2 Kellogg,3 and Osborn4 (e.g., for oxides, sulphides, halides) from critical surveys of existing data. The speed and simplicity of making calculations with these graphs make them of great practical value and, this should be stressed, in the majority of cases no greater accuracy is obtained by the individual calculations from specific heat data, entropies, and heats of reaction. With the increasing industrial application of vacuum methods to the winning of metals, notably at present magnesium and the alkaline earth metals, it is necessary to calculate low metal vapor pressures above reaction mixtures. The existing free energy diagrams usually give values where the solid and liquid metals are at unit activity. To obtain the free energy change associated with the reaction Metal (solid or liquid) = Metal (gas) it is customary to calculate this from either vapor pressure equations or from the appropriate free energy functions. This paper presents graphically the free energies of vaporization of a number of metals of economic importance which can be used in conjunction with other free energy diagrams. Certain metals have been excluded on the following grounds: titanium, tungsten, molybdenum, and platinum have very high heats of vaporization and boiling points and the uncertainty of these values arising from difficulties of high temperature measurements makes free energy values of doubtful accuracy; the data for cobalt, vanadium, and zirconium are uncertain; finally, arsenic, antimony, and bismuth could not be conveniently included. because of the presence of mixtures of polyatomic and mona-tomic vapors in the temperature range considered. Free Energy Diagram In Fig. 1, the standard free energies of vaporization per mol (?F°T) have been plotted for the temperature range 0" to 2000°C for the reaction M (solid or liquid) = M (gas) The data have been obtained from the free energy functions — (F-H298)/T for the solid, liquid, and gaseous elements and the heats of vaporization at 298°K, i.e., ?H298 (vapor). The values of these functions have been taken without change from the recent critical survey by Brewer;' the bulk of this data is from Kelley7 with adjustments in the light of more recent work. Values of the function — (F-H298) /T at T = 298°, 500°, 1000°, 1500°, and 2000°K were employed. Where possible, a further independent value was obtained by using the boiling point, i.e., where ?F°T = 0. All calculations are based on the assumption that the metal vapors are monatomic and that the boiling points refer to those temperatures at which the partial pressure of the monatomic vapor above the metal is 760 mm. For the alkali metals, e.g., sodium, potassium, and lithium, diatomic gas molecules are present as minor species; details of these are discussed by Brewer." On the basis of the above values and allowing for the change in slope at the melting point, straight line graphs were constructed and extrapolated to 2000°C. The deviation from a straight line did not exceed ±500 cal, which in all cases appears less than the probable error in the ?F°T values.

Jan 1, 1954

By S. Shapiro, G. Krauss

RECENT investigations1 3 of the products of marten-sitic transformation in Fe-Ni and Fe-Ni-C alloys have made use of light microscopy to describe mar-tensitic fine structure. The application of conventional metallographic techniques to this problem was possible because of the similarity of distributions of parallel striations developed by etching polished surfaces to distributions of transformation twins which are positively identified in thin films by transmission electron microscopy. The surface striations frequently are not well-defined, and could also be due to, for example, rows of etch pits which are not resolved by the light microscope. This note describes the etching effects revealed by examination of surface replicas of mar-tensite in the electron microscope. The application of a replica technique both extends the observations of fine structure made with the light microscope and allows the results of conventional metallography to be used with more confidence. Samples of an Fe-33 wt pct Ni alloy were used for this study because of the detailed knowledge of the fine structure of martensite which has been developed'-4 for alloys of about this composition. Samples transformed to martensite at -120°C, and transformed samples heated to 340°C to effect a small amount of reversal of martensite to austenite, were examined. Surfaces were first prepared for replication by elec-tropolishing in a chromium oxide-acetic acid electrolyte at a current density of 0.1 amp per sq cm and a temperature of 15°C or a perchloric acid-acetic acid electrolyte at a current density of approximately 2 amps per sq cm and a temperature of 15°C. The polished surfaces were then etched for 4 to 7 sec by immersion in a solution of equal parts HC1, HNO,, and H20. Application of the etchant for this period of time produced discrete attack of certain features without loss of resolution due to overlapping of areas most susceptible to etching. The etched surfaces were replicated by 1) preshadowing with Pt-Pd alloy at an angle of 20 to 30 deg, 2) depositing a thin carbon film, and 3) stripping the film and shadowing material by immersion in a 5 pct nital solution. Two types of etching effects were observed by examination of the surface replicas in the electron microscope: l) continuous striations centered at the bottom of long grooves or the tops of pyramidal elevations on the sample surface, and 2) high densities of fine etch pits. Examples of both types of features are shown in Figs. 1 to 4. The striations are in each case continuous and have the same characteristic configurations as do transformation twins in martensite. Fig. 1 shows that the density of grooves decreases toward the edge of a martensite plate and Fig. 2 shows a region where the striations are of varying length. Similar examples of transformation twin arrangements determined by transmission electron microscopy of thinned foils of martensite have been shown in the literature,'~ and, based on these similarities, the striations and their associated grooves or mounds are considered to be evidence for the presence of transformation twins. The transf~rmation twins in martensite are on the order of l00 A thick, much too fine for resolution by the light microscope. It appears that orientation-dependent etching of the matrix martensite about the twins results in grooves or surface elevations which

Jan 1, 1967

By B. D. Smith, K. H. Coats

Experiments in which calcium chloride displaced sodium chloride from four cores showed the extent of asymmetry in the resulting effluent concentration profiles. These results provided a check on how validly the mixing process is modeled by a differential (i.e., not finite-stage) capacitance mathematical model. The effluent concentration profile from two consolidated cores exhibited considerable asymmetry, while two unconsolidated cores yielded nearly symmetrical profiles. All runs resulted in breakthrough of the 50 per cent concentration significantly before one pore volume was injected. In addition, velocity appreciably affected the effluent concentration profile from a Torpedo sandstone core. The differential capacitance model matched the data significantly better than a simple diffusion model. The capacitance model allows determination of the amount of dead-end pore space in a porous matrix and the effect of velocity on the rate of diffusion into this space. An experimental program yielding insight into the physical validity of the capacitance effect is described. INTRODUCTION Axial dispersion — the mixing accompanying the flow of miscible fluids through porous media — has been the subject of many relativeiy recent studiesl-13 and a comprehensive review of the topic has been given by Perkins and Johnston.14 This dispersion is of practical interest in studies of the miscible displacement process, fixed-bed chemical reactors, and the adsorption of solutes from a flowing stream onto the surface of a porous medium. In the latter case, the effect of dispersion must be considered when adsorption parameters are determined from the nature of concentration profiles. In general, early studies of dispersion assumed applicability of a simple diffusion equation and were concerned with correlation of the experimentallv determined "effective" diffusion coefficient with system properties over a large range of the latter. Recent have been concerned with the deviations between the asymmetrical effluent concentration profiles observed and the symmetrical ones predicted by the diffusion model. In the present study, effluent concentration profiles were obtained from consolidated and unconsolidated cores. These profiles were compared with those predicted by a differential (i.e., not finite-stage) capacitance model. Solutions to the simple diffusion model, for three sets of boundary conditions, were compared with one another and with the experimental profiles. SUMMARY OF PREVIOUS WORK The reader is referred to Perkins and johnston14 for an extensive review of studies of dispersion in porous media. Many investigators have employed the simple diffusion model characterized by Eq. 1 below: The dispersion coefficient D for unconsolidated systems is correlated by u is interstitial velocity and dp is particle diameter. Since heterogeneity of the sand pack affects the mixing, this equation is also expressed as where a is proportional to the degree of heterogeneity and is about 3.5 for random packs of unconsolidated sand.l4 Eq. 3 holds for 2 For a Do homogeneous (regular) type of packing, should be 1 or less.6 Aris and Amundson6 and Carberry and Bretton consider o to be the number of particle lengths per mixing cell in the finite-stage model. Data from consolidated cores indicate dp to be about 0.36 cm for outcrop rocks, Torpedo sandstone having a reported value of 0.17 cm.15 For v /Do less than about 2, the effect of molecular diffusion must be considered and the total dispersion coeffi-

Jan 1, 1965

By F. F. Craig

A series of both welter and gas puttern Hoods was made in the laboratory to study the oil recovery performances of .such operations. These tests were contlrrctt~tl on consolidated sandstone models, using oil. water, and gas. The model floods were scaled to reprotlrrce field performance under gas arid water five-spot injection. X-ray shadowgraphs permitted observation of the gross fluid movement within the models. A method was developed for applying the mobility ratio concept to water flooding and dispersed gas drives in a five-spot well pattern. The area1 sweep efficiency. tit breakthrough for dispersed gas drives is much higher than previorsly expected, Iying in the range of 50 to 100 per cent. A method is presented for predictitig the water-oil ratio performance of five-spot pattern water flood., in uniform sands. This method is verified experimentally for the cotidition of 120 free gas initially present and for va1ues of gas saturation normally encountered in fields following depletion operations. Prodrrction perfort?zut~cc for pattern gas injection is also predictable by this method. INTRODUCTION In field pattern hater flooding or dispersed gas injection operations, the displacing fluid (water or gas) is injected and oil is produced through wells which are, area-wise. small openings in a large container or rescrwir of oil. As a result, not all of the area between the injection and producing wells is necessarily contacted by the injected fluid by the time it first reaches the producing wells. The question arises in predicting oil recovery. as to what traction of the pattern area involved is contacted by either the injected water or gas due to the relative position of the injection and producing wells. This is the areal sweep efficiency. It is also desirable to know how well arrangement affects the oil production performance. A solution to this problem has been attempted by many means. Mathematical analysis¹,²,³, and electrolytic models', have been used to establish a breakthrough areal sweep efficiency of 73 per cent for a five-spot pattern in which the ability of the oil to flow ahead of the invading fluid is equal to the ability of the invading fluid to flow in the invaded zone. All of these studies simulate only one of many field conditions. Improved studies of areal sweep efficiency at breakthrough were made using the fluid mapper5, poten-tiometric models" and X-ray shadowgraph techniques. In these it was possible to study the effect of variations in mobility of the injected and displaced fluids, and so to cover a wide range of field conditions. The limitation of these studies is that, with these models, it was not possible to simulate the saturation gradient behind the flood front when oil is displaced by gas or water. In the most recent published technical paperS on the subject, a method was presented for predicting the oil recovery performance after breakthrough in a five-spot pattern water flood. The method involved correlations which were obtained experimentally using miscible fluids. One limitation of this work was that a saturation gradient, usually found in immiscible fluid displacement, was absent in these studies. Also the laboratory studies did not cover the situation of a water flood following oil recovery by depletion drive. as is the situation in most field water floods. The tests reported in this paper were made using oil. water, and gas, thus eliminating the simplifications of former studies. Also reported are laboratory studies on

Jan 1, 1956

By Robert B. Hall, Charles W. Bauer

Alunite, KA13(SO4)2(OH)6, was used from the 15th century until early in the present century as a source of potash alum and aluminum sulfate. Since early in this century its use for this purpose has declined to virtually nil with acid treatment of bauxite or clay replacing alunite as a source of alum. During the First World War alunite served as an emergency source of potassium sulfate fertilizer in the United States and Australia. More recently, alunite has been investigated as a nonbauxite ore of aluminum, with potassium sulfate and sulfuric acid as recoverable byproducts. The Soviet Union established a commercial-scale plant in the mid-1960s producing cell-grade alumina from alunitic ore and recovering byproduct potassium sulfate and sulfuric acid. Interest in nonbauxite ores, including alunite, is expected to continue in countries that presently are dependent on imports of bauxite or alumina to support their aluminum industries. The US has the largest aluminum producing and fabricating industry in the world but imports more than 90% of its aluminum-bearing raw materials. This chapter reviews the status of alunite as a potential source of aluminum with byproduct fertilizer, and the resources available to support this industry, whenever this becomes economically feasible. Historical Background Alunite has received little attention until recently, yet for centuries it held a higher place than many of the industrial minerals considered essential today. The deposits of alunite at Tolfa, northwest of Rome, were mined almost continuously from 1462 until about the middle of the present century (DeLaunay, 1907; Lombardi, 1977). The potash alum made from the alunite, called Roman alum, was a valued commodity in international trade, being used in the making of paper, in tanning leather, and as a mordant in textile dyeing. The industry at Tolfa is said to have been controlled by the Vatican and was a source of revenue for the Papal State until fairly recent times (Lombardi et al., 1977). Since early in the present century, the preferred method for making alum has been by acid treatment of bauxite, clay, or other aluminous material, so that alunite has fallen into near-total disuse as a source of alum chemicals. Appreciable quantities of mixed alunite-kaolin rock still are quarried in the Tolfa district for use by the Italian cement industry (Lombardi and Mattias, 1979), but the role of this material in the making of Italian cement is not clear. A large deposit of alunite-bearing rock in Chekiang province, mainland China, is said to have yielded a "great quantity" of alum over a period of several centuries (Yih, 1931). During the 19th century, alunite deposits in France, Spain, Australia, and other countries were exploited as raw material for making potash alum and aluminum sulfate. Interest in alunite was revived during the First World War, not for alum, but as an emergency source of potassium sulfate fertilizer. Large veins of nearly pure alunite had just been discovered on Alunite Ridge, 11 km southwest of Marysvale in Piute County, Utah (Butler and Gale, 1912), and these were quickly brought into production for making potash fertilizer to replace supplies cut off from the famous deposits at Stassfurt, Germany. An estimated 225 kt of Marysvale alunite was mined and processed at that time (Callaghan, 1973). Similarly, Australian alunite deposits were mined as a source of potash fertilizer during both world wars (Hall, 1978).

Jan 1, 1983

By Donald E. Grimes

In conjunction with Inland Steel&apos;s development of lead-bearing steels possessing improved machin-ability because of tellurium and/or selenium additions, it was decided to determine liquidus and solidus temperatures in the Pb-PbTe-PbSe sub-ternary system. With the exception of alloys in the PbTe-PbSe system, it is believed that no previous study of this subternary has been performed. Thermal analysis was employed in a general survey of the subternary. In addition to PbTe, PbSe, and pure lead, ten ternary compositions were selected at intervals of ten at. pct. Sixty-gram alloys were prepared by sealing an appropriate mixture of the elements in quartz crucibles with a central thermocouple-well under a vacuum of less than 10 p. All three elements were 99.999+ pct pure with only trace impurities detected by spectroscopic examination. After the alloys had been heated at least 20°C above the liquidus in a Kanthal furnace, cooling and heating curves were obtained. Two or more cooling curves were employed both for the determination of each alloy liquidus and for the solidus of those al- loys containing greater than 50 at. pct Pb. Cooling rates were 0.5° to 3.5°C per min during the 5 min immediately prior to the start or end of solidification. For those ternary alloys containing 50 at. pct Pb, the solidus was determined from heating curves. Temperatures were measured with a Pt/Pt-10 pct Rh thermocouple which was calibrated in the same apparatus against the freezing points of zinc (A. D. MacKay—99.999 pct), antimony (A. D. MacKay-99.999 pct), and copper (ASARCO—99.999+ pct) both before and after the experimentation. The only interruptions in the continuous strip-chart recording of the temperature occurred for brief periods at significant points of interest when direct readings were made with a Rubicon Type B potentiometer sensitive to * 1 µv. The average liquidus and solidus temperatures obtained in this investigation are listed in Table I together with probable errors based on the accuracy of the thermocouple calibration and the reproduci-bility of the data. No chemical analyses of the alloys were performed because the synthetic or as-charged compositions were felt to be reliable. The melting points obtained for PbTe (924.7° ± 0.7°C) and PbSe (1084.0° ± 0.6°C) compare favorably with recent findings of other investigators.&apos;-&apos; As indicated in Table I, both the liquidus and solidus temperatures increase as selenium replaces tellurium for a constant atomic percent lead. As a result, the ternary liquidus and solidus temperatures for a given percent lead lie between the corresponding binary temperatures published in the literature. It should be noted that those ternary alloys containing 50 at. pct Pb exhibit a relatively narrow liquid-plus-solid region (<80°C) when compared with alloys containing greater than 50 at. pct Pb. Alloys with 60, 70, or 80 at. pct Pb not only have a much wider liquid-plus-solid region (>500°C) but also exhibit considerably lower and approximately constant solidus temperatures which thus form a virtually isothermal solidus "plane" slightly below the melting point of pure lead (327.3°C). For these alloys

Jan 1, 1965

By Will Mitchell, Ford F. Miskell, C. L. Sollenberger

Heat treatment of ore prior to comminution reduces power requirements for grinding, reduces grinding media wear, and improves recovery of iron values from a typical Minnesota magnetic taconite. Test data demonstrating this, as well as an analysis of the economics of commercial application of the technique, are presented. THE taconites in general are hard, tough ores, difficult to grind. Liberation of iron mineral constituents usually is accomplished by grinding the ore through at least 100 mesh, and often it has been found necessary to grind substantially through 325 mesh to achieve satisfactory recovery and grade in the concentration process. Because of the fineness of grind required and the enormous tonnages of material contemplated for treatment in the future, costs resulting from grinding media wear and power consumption, together with capital investment required for comminution, approach astronomical figures. Economy in any one of these elements per ton of material ground could very well reflect a considerable saving to the ferrous industry in yearly costs. With this in mind, the Research Laboratories of Allis-Chalmers Manufacturing Co. have launched a program to discover means of effecting this economy. The initial phase of the work as described here deals with a heat treatment of crushed raw ore, followed by thermal shock in cooling, for the preparation of rod mill feed. Several investigators1-" ave noted an improvement in the grindabilities of ores treated in this manner. One investigator2 subjected low grade iron ores to heat treatment in an electric furnace. Basing his conclusions on screen analyses of crushed products, he observed that treated ore was more easily crushed than the untreated. He reported that observation of the treated and untreated ore through a microscope revealed cracks following the grain boundaries in treated ore, whereas no cracks were present in untreated ore. However, few if any have quantified the improvement in terms of total hp hours saved, reduction of wear of grinding media, or reduction of capital costs of grinding equipment involved. By means of Bond rod mill and ball mill grindabilitiesh nd by comparative wear tests, the conditions maintained during this investigation have been definitely evaluated. Heat treatment experiments have been made on a batch scale in a muffle furnace, followed by continuous scale experiments in a rotary kiln under various conditions of temperature, atmosphere, retention time, and quenching to determine the combination that would give the greatest improvement in the grindability of a taconite without affecting adversely the magnetic susceptibility of the magnetite. From these data, relative costs have been calculated for grinding both untreated and heat- treated ore on a basis of plant capacity of 120 tons per hour. Heat requirements for a rotary kiln, as well as the kiln capacity required for the treatment, have been estimated. For these tests 50 tons of ore were obtained from the Aurora mine of the Erie Mining Co., Hibbing, Minn. The ore was reported7 to have been selected from the magnetic portion of the lower chert horizons. An average chemical analysis of&apos; this portion is shown in Table I and an approximate mineralogi-cal analysis is shown in Table 11. No attempt was made to differentiate between the various silicates because the compositions vary widely on the different horizons and depend on the degree of oxidation. Minnesotite, however, is the predominant silicate. The siderite grains were dispersed throughout the ore but were essentially associated with the rock-forming minerals. In a series of preliminary tests in which batches of —6 mesh ore were heat-treated in a muffle fur-

Jan 1, 1953

By P. Chiotti, E. R. Stevens

The electromotive force between pure magnesium and Mg-Zn alloys in a fused KC1-LiCl-MgC12 cell was measured over the temperature range 360° to 730°C and for alloy compositions of 0.052 to 0.635 atom fraction magnesium. The calculated activity coefficients of magnesium in the liquid alloys may he represented by the relation These data and the phase diagram were employed in calculating the standard free energy of formation for the compounds MgZn, MgZn2, and MgZn2, In the temperature range from 25° to 590°C the relation ZINC-rich magnesium alloys have been shown to be very effective in reducing various solute chlorides or oxides from fused salts such as KC1-LiC1 eutectic. Such reactions are of interest in the developinent of pyrometallurgical methods for reprocessing nuclear reactor fuels.&apos; As the conceniration of the magnesium is increased beyond 5 to 10 wt pct the extent to which solutes, such as yttrium, lanthanum, cerium, and other rare earths, can be reduced or removed from the salt decreases. This effect is due to the resulting change in the activity of the solutes in the Zn-Mg alloy as well as the change in the activity of both the zinc and mag- nesium. The purpose of the present investigation was to determine the thermodynamic properties of the Zn-Mg system which, in addition to yielding information on the alloying behavior of these two elements, would permit a more quantitative understanding of the above reactions. Since this investigation was initiated several publications have appeared in the literature giving results on the thermodynamic properties for the liquid alloys of this system. 2-4 Kozuka, Moriyama, and Kushima3 measured the vapor pressure of zinc over four different alloy compositions by a carrier-gas technique. Eremenko and Tukashenko2 and Terpilowski4 measured the electromotive force between pure magnesium and various alloys in a fused salt cell. Earlier work has been summarized by Hultgren, Orr, Anderson, and Kelley.5 In the present investigation electromotive-force measurements were made on a cell of the type The temperature range covered for the lower-melting alloys was roughly 375° to 730°C. These data were employed in calculating the thermodynamic properties for the liquid alloys and in conjunction with the phase diagram,&apos; reproduced in Fig. 1, were used in calculating the standard free energy of formation of the compounds MgZn, MgZn2, and Mg2 Zn11. The results obtained are compared with available data from the literature. MATERIALS AND EXPERIMENTAL PROCEDURES Materials employed were: Bunker Hill slab zinc, 99.99 pct pure; Ames Laboratory magnesium containing 13 ppm N, 195 ppm C, and trace amounts of calcium, copper, nickel, and silicon; reagent grade KC1, LiC1; and high-purity MgC12 obtained from Dow Chemical Co. The electrolyte consisting of KC1-LiCl eutectic (44.4 wt pct LiC1) to which was added 15 wt pct MgCl2, was dried in a quartz tube by heating to 340° C under vacuum for a period of 4 to 5 hr, and then passing HC1 through the salt as the temperature was slowly raised to above the melting point of the salt. After about 3 hr the HC1 was shut off and the system flushed with argon and then vacuum applied for 1 hr. The salt was finally siphoned through a sintered glass filter and cast into a Pyrex container which was sealed and

Jan 1, 1965

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

On the basis of metallographic analysis, incipient melting data, thermal analysis work, and X-ray diffraction, phase relationships in the 0 to 50 atomic pct Cr region were carefully resolved. Phase relationships in the 50 to 100 atomic pct region were outlined. A single intermediate phase, ZrCr,, was established in the system. A eutectic and a eutectoid were found in the zirconium-rich region of the system, while a eutectic was located in the chromium-rich region. Limited solubility of chromium in both a and a zirconium was observed. THIS is another in a series of papers on zirconium-base binary phase diagrams based on a program sponsored by the Atomic Energy Commission. The Zr-Sn,&apos; Zr-Si,&apos; and Zr-Mo and Zr-W" systems have already been published. Previous work on the Zr-Cr system includes a partial diagram by Hayes, Roberson, and Davies&apos; and a complete diagram by McQuillan." The phase relationships reported by these authors were determined with magnesium-reduced zirconium, whereas those reported herein are based on iodide metal. Some X-ray work on the intermediate phase in this system was also reported by Hayes et al., as well as by Wallbaum.&apos; Comparisons with these data are made in the appropriate sections. Thermal analysis, metallographic and incipient melting techniques were employed to accurately resolve phase relationships in the 0 to 50 atomic pct Cr region, while the remainder of the diagram is outlined on a basis of cast structures and thermal analysis. A description and analyses of the metals used in the preparation of alloys for this study are included in Table I. A "low-hafnium" zirconium crystal bar (Westinghouse "Grade 3") was used. This material was produced by the decomposition of a volatile iodide onto a hot zirconium filament. The as-received crystal bar was coated with corrosion product from a standard autoclave test by which its grade designation is determined. A sand-blasting and subsequent HF-HNO, pickling treatment was employed to remove this contamination. The bars were cold-rolled to approximately 1/32-in. sheet, pickled again, sheared to 1/4-in. squares, washed with acetone, and stored for use. Ingot chromium metal from National Research Corp. and electrolytic chromium flakes from Johnson, Matthey and Co., Ltd., were available as alloy additions. The ingot chromium was received in granular form. The large granules were sized by crushing, treated for iron removal, rinsed with acetone, and stored. The flakes were broken to approximately 1/8 in. with mortar and pestle, cleaned with acetone, and stored for charging. Experimental Procedure A nonconsumable electrode arc furnace identical to that described by Hansen, Kessler, and McPher-son,&apos; was used in this work. A 400 amp dc welding generator was the source of power. To insure homogeneity, ingots were melted in a water-cooled spun copper crucible, a total of four times without opening the furnace. No difficulties were encountered in the preparation of these binary alloys. Analyzed chromium contents were found to agree very well with the intended content. Melts of pure zirconium were interspersed in the course of preparing the alloys and their hardness was measured to keep a constant check on the melting technique. Vycor and quartz bulbs were used to contain the samples for isothermal annealing treatments. Vycor bulbs were employed for temperatures up to 110O°C, while quartz was used above this temperature. For treatments up to 950°C the bulbs were evacuated before sealing, but for higher temperatures a partial pressure of argon was admitted to the

Jan 1, 1954

By H. H. Rachford

This work presents a first-order analysis of the instability underlying viscous fingering in adverse viscosity-ratio water floods. It extends previous analyses of frontal instabilities, which were carried out with equations for parallel plate models, by including effects of the saturation transition zone observed behind the front in water floods in water wet systems. This zone tends to insulate incipient fingers from the high-mobility watet;. thus conditions for the onset of fingering differ from those in the parallel plate theory. Finite-difference solutions of the two-dimensional equations of displacement in porous media exhibited the predicted stability characteristics in six hypothetical field- and laboratory-scale floods in rectangular reservoirs. In contrast to results with parallel plate systems, this paper concludes that for waterwet reservoirs, laboratory models scaled by the usual criteria are also correctly scaled for frontal instability. Further, fingering in the systems studied can occur in any saturation range behind the front, and may occur at an intermediate saturation even though stability obtains both at the saturation corresponding to the Buckley-Leverett front and near residual oil saturation. Other points of contrast are that the likelihood of occurrence of fingering may not increase as flow rate or viscosity difference increases, but my be sensitive to changes in he relative permealdlity and capillary pressure functions. INTRODUCTION The recovery of oil by water flooding frequently involves displacing the oil by water of a lower viscosity. Displacement of a fluid by a less viscous one may lead to gross channeling or fingering like that observed in solvent floods, in which it severely lowers recovery efficiency.1 In addition to adverse effects on recovery, it has been suggested that the unstable move- ment causing fingering may interfere with interpretation of scaled model studies of proposed water floods, since the instability in the model might not be faithfully scaled to that in the reservoir prototype. In view of the serious implications of this possible breakdown of widely used scaled model techniques4-6 it is the purpose of this paper to examine the question further. Instabilities in the solutions of systems of differential equations imply a loss of smooth dependence on initial and boundary conditions. Thus, the possibility exists that in using models whose scaling is based on the differential system there may arise size- and rate-dependent factors which are not properly scaled. This possibility was examined in detail by Chouke et al., who analyzed the instability of frontal advance in a related problem, water-oil displacement in parallel plate models,2 in which a moving interface separates two regions of constant, unequal mobilities. First-order perturbation theory predicts the existence of a critical wave length for the growth of perturbations: and a wave length of maximum instability of A, = The interpretation is that wave lengths in a perturbation which are longer than A, will grow. Thus, if the width of a two-dimensional channel is greater than A,, fingers will grow, and the spacing of the fingers which grow at the maximum rate will be approximately It is important that the higher the velocity and/or the difference in flow resistance, the lower is A,, and thus the greater the number of fingers that can grow in a given model. In applying these conclusions to porous media, a pseudo-interfacial tension, , was assumed for the invasion front. Since this would not necessarily be equal to the liquid-liquid interfacial tension, an unknown constant was substituted for the foregoing expression for A, This was determined

Jan 1, 1965

By W. J. Helfrich, C. L. Sollenberger

Many present direct reduction processes utilize iron ore concentrates for the production of sponge iron and the sponge iron is usually preferred as an agglomerate. Pelletizing a high grade iron oxide concentrate prior to reduction is a simple approach to the high capacity production of a uniform and agglomerated sponge iron product. The present paper describes the effect of different atmospheres and physical characteristics on the reduction of pelletized iron oxide concentrates. The kinetics of iron oxide reduction have been I studied by numerous investigators1-4 and a rate equation for the hydrogen reduction of hematite and magnetite has been proposed. Mc Kewan&apos; obtained kinetic data for the reduction of hematite spheres sintered at 13&apos;70°C (2498&apos;F&apos;). The kinetics of magnetite reduction was studied by Quets. He compacted the oxide at a pressure of 1 ton per sq cm and sintered the specimens to densities of 95 to 96 pct of theoretical. These studies demonstrated that the reduction of iron oxide is a topochemical process and that the rate of reduction is proportional to the receding surface area of remaining iron oxide. Thus, the measurement of the kinetics of iron oxide reduction depends upon the development of a measurable and well defined interface between metal and reacting oxide. It is evident that the reduction rate will be constant only when the rate of formation of the reaction product layer is constant. It was recognized by Joseph&apos; and Edstrom&apos; that the physical character of the oxide will determine the configuration of the interface between oxide and metal. The reduction of oxides which exhibit few pores or no interconnected pores should proceed by topochemical processes. However, if a massive or dense crystalline ore cracks during the reduction of the oxide, thus allowing the reduction to Proceed away from a well defined oxide/metal interface, then the observed reduction rate will vary or change during reduction. From a process standpoint, it would be difficult to apply absolute reduction rates determined for crystalline or nonporous iron oxide to the direct reduction of a porous pellet of iron ore concentrate. Most present direct reduction processes are attempts to utilize iron ore concentraates in the production of sponge iron. The sponge iron is usually preferred in an agglomerated form for subsequent melting practices. Pelletizing a high grade iron oxide concentrate prior to reduction is a simple approach to the high capacity production of a uniform and agglomerated sponge iron product. One process for preparing such a product has already been de-scribed. This paper describes the effect of different atmospheres and physical characteristics on the reduction of pelletized iron oxide concentrates. ORE CONCENTRATES STUDIED A magnetite concentrate from Reserve Mining Co. and a hematite concentrate from U.S. Steel Corp. were used in preparing pellets for these experiments. A head analysis of these concentrates is given in Table I. To compare the reduction rates of pelletized concentrates, it was necessary that each pellet contain equal amounts of oxygen which could be removed by reduction, regardless of the oxidation state, i.e., Fez% or Fe304. Also, because of the difficulty in preparing spherical pellets of uniform size, weight, density, and porosity for experimental purposes, pellets of a cylindrical shape were formed in a press.

Jan 1, 1961

By H. H. Rachford, D. W. Peaceman

A system of partial differential equations describing miscible displacement of fluids in porous media is derived The system takes into account the influence of gravity, the spatial distribution of permeability, diffusion, and fluid viscosities and densities. A numerical procedure for approximating solutions to the differential systems has been tested for a horizontal two-dimensional geometry. In end-to-end displacements of oil with less-viscous solvents, the numerical solutions exhibited fingering qualitatively similar to that observed in laboratory models. Small random spatial variations in permeability about the mean value are sufficient to initiate fingering. Quantitative comparisons of computed results with laboratory data show good agreement. INTRODUCTION Miscible flooding of oil by solvent is receiving increasing consideration for field use.1-4 Unfortunately, such floods potentially present severe problems in loss of recovery through the by-passing of oil by expensive solvent. Consequently, their economic evaluation requires sound techniques for predicting recovery. The purposes of this work are to present a finite-difference method for calculating the multidimensional displacement of oil by solvent and to investigate the validity of the method by comparing results of calculations with data from displacements in laboratory models.5 The formulation of the method to simulate the model experiments treats a case of limited scope in the description of solvent flooding in the reservoir. The model experiments were carried out with fluids that were assumed to form an incompressible, ideal, two-component system with constant diffusivity. Establishing the validity of the method even in its present form provides a major step toward the goal of quantitatively evaluating individual solvent flooding projects. First, it demonstrates the fea- sibility of calculating the course of displacements which are dominated by an inherent macroscopic instability, i.e., viscous fingering. In addition, the method provides almost the only practical means of examining the effect of the size and extent of reservoir inhomogeneities on the development and propagation of the fingers, taking into account the important influences of diffusion and gravitational segregation. It might be inferred that such studies cannot yield results of practical value because in displacements dominated by macroscopic instability the pattern of finger development and the resulting performance should depend critically on small variations in the geometric distribution of reservoir inhomogeneities. The present method offers the capability inherent in computational techniques of predicting performance reproducibly with any arbitrary distribution of reservoir properties and, as such, provides a means of evaluating the sensitivity of behavior to uncertainties in reservoir definition and, thus, assessing the reliability of prediction of performance. THE PHYSICAL PROBLEM DATA USED FOR QUANTITATIVE TESTS OF CALCULATIONS In this work, a calculation for treating the miscible-displacement process is tested by comparing calculated results with experimental data. The experiments chosen for comparison have been described in detail,5 and the quantitative features are summarized in a later section. Briefly, oil was flooded by solvent of equal density from a thin rectangular channel in Lucite packed with uniform Ottawa sand. Under these conditions, two-dimensional geometry is considered adequate to represent the process. Inasmuch as the present work is concerned with testing of the method by comparison with specific two-dimensional data, the system of equations presented in the succeeding section will be oriented to defining the physical system specifically. THE DIFFERENTIAL EQUATIONS Suppose x, y to be a Cartesian co-ordinate system and define h(x, y) to be the height of a point above a horizontal reference plane. It will be assumed

By F. D. Rosi, C. A. Dube, B. H. Alexander

WHEN a deformed polpcrystalline metal is heated to a sufficiently high temperature, a recrystallized structure develops which consists of small, essentially stress-free grains. This transformation is referred to as primary recrystallization. Under certain conditions, on prolonged annealing of the specimens, a type of grain coarsening occurs in which a small number of grains commence to grow at various points by devouring the neighboring finegrained primary structure, which, in itself, exhibits a marked stability toward normal grain growth. This grain coarsening process has been referred to by Burgers&apos; as "secondary recrystallization" or "exaggerated grain growth." Although much work has been done on recrystallization and growth processes, present knowledge of the factors which determine the occurrence and growth of secondary grains is inadequate for the advancement of a theory which accounts for the observed characteristics of this recrystallization phenomenon. In this respect, it would be important to make measurements of the velocity of growth of these large secondary grains. The early work of Feitnecht2 on aluminum showed that the size of secondary grains increased with temperature and time, suggesting that such measurements would be possible. More recently, Burgers1 found that the growth velocity of these grains in aluminum was constant over a large time interval, as is generally found for the growth velocity of primary grains in a uniformly deformed test piece.3-5 It is reasonable to expect that the determination of the activation energy of the secondary growth process might lead to a better understanding of the underlying atomic mechanism. This investigation, therefore, is one of a series primarily designed to measure the growth velocity of secondary grains at various temperatures for silver, copper, and aluminum. It is hoped that these studies along with those of orientation relationships, will establish a sound basis for theorizing as to the origin and nature of the secondary transformation. Experimental Procedure The metal used in the present investigation was high-purity silver (999.9 fine), which was supplied in the form of cold-rolled plates 0.250 in. thick. These plates were annealed at 500°C for 1 hr and then straight-rolled to a thickness of 0.04 in. in three reductions of approximately 45 pct, followed in each case by an anneal at 450°C for 1/2 hr. The 0.04 in. sheets were given a final reduction of approximately 50 pct, which resulted in a thickness of 0.02 in., thereby insuring two-dimensional crystal growth. Individual specimens of this material about 1x3 cm were used for all heat treatments, unless otherwise noted. The annealing was carried out under an atmosphere of argon in an Inconel tube furnace whose temperature was controlled to within -+2°C. Metallographic examination of the secondary transformation was accomplished by mechanically polishing the specimens on one side to approximately one half their thickness through 2/0 metallographic emery paper, and then electrolytically polishing in a solution of potassium ferrocyanide and sodium cyanide. With a current density of approximately 20 amp per sq dm, the time required to obtain a satisfactory surface varied from 3 to 5 min. Standard potassium bichromate was used to etch the electrolytically polished surfaces. A heavy etch with this solution produced a sharp contrast between the secondary grains and the fine-grained primary structure, darkening the latter while leaving the large secondary grains bright. All area measurements of the secondary crystals were made with a planimeter on enlarged tracings of the grains.

Jan 1, 1953

By B. H. Caudle, A. R. Khan

The oil production performance of thin-oil-column type reservoirs producing entirely by natural water drive are discussed in this paper. The dimensionless production histories were obtained in a series of scaled porous-medium laboratory-model tests under varying effects of rate of production, mobility ratio across moving boundary, ratio of oil- and water-sand thickness and mobility ratio across fixed boundary including the effects of different residual oil saturations in the water-invaded regions. This scaled-model study of thin oil columns indicated that greater recovery of oil results for lower values of the following: rate of production, mobility ratio across the moving boundary, water-sand thickness relative to the oil-sand thickness, and mobility ratio across the fixed boundary or residual oil saturation in the water-invaded region of oil sand The results of the study are given as dimensionless production history plots which may be used to predict the performance of the field they represent. INTRODUCTION Oil fields producing under natural water drive can be classified into two major groups: (1) an "edge-water drive" reservoir group in which the water is encroaching from the reservoir periphery, and (2) the group of reservoirs in which the oil-sand is entirely underlain by an aquifer. Reservoirs in the second group are further subclassified as "bottom-water drive" reservoirs and "thin-oil-column" reservoirs. The latter reservoirs are characterized by their relatively equal oil- and water-sand thicknesses and by having both the oil and water zone relatively thin. When oil is produced from a water-drive reservoir by a partially penetrating well, water breaks through sooner or later before all of the displace-able oil is recovered. Many analyses of this general situation have been made with varying degrees of realism in their approaches, both theoretical and experimental5.6 in nature. Pertaining to the scaled laboratory model studies of natural water drives using porous media, detailed investigations have been made by Henley, et al.7 on bottom-water drive and by Caudle8 on edge-water drive reservoirs. No laboratory model studies, however, have been made to investigate the performance characteristics of thin-oil-column reservoirs. This paper presents the results of a study on scaled models of thin-oil-column reservoirs, operated under varying conditions corresponding to the actual field practices. SCALING PARAMETERS The object of scaling is to produce a one-to-one relationship between physical and dynamic parameters of the model -and prototype. The most comprehensive treatment of scaling techniques for oil reservoirs is given by Geertsma, et a1.,9 which covers two methods for arriving at dimensionless groups for scaling. Table 1 lists the scaling parameters suitable for creating kinematic and dynamic similarities between the model and the prototype in the variables considered for the thin-oil-column reservoirs. In order to achieve similarity it is essential that the dimensionless groups listed in Table 1 have equal values in the model and the protorype simulated. Since it has been established that the role of interfacial forces is neglible in gross fluid movement in a reservoir no attempt was made to scale capillary forces- Although there are reports describing Scaling criteria for these forces,10,ll consideration of the necessary but restrictive conditions for

Jan 1, 1970

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 R. H. Doremus, E. F. Koch

The first carbide phase that precipitates at 120°C from a iron containing about 0.02 wt pct, C was studied with the electron microscope. In both strained and unstrained material the carbide particles appeared to form on the dislocations in the iron, in agreement with the results of the kinetic study in part II. The structure of the carbide was tentatively identified as an expanded a-iron lattice (body-centered tetragonal) with about 0.7 wt pct C. This phase appeared to be metastable and apparently transformed into cementite or epsilon carbide when the iron was over-aged. WhEN a iron containing less than about 0.025 wt pct C is quenched from 700°C to room temperature the carbon is retained in solid solution and precipitates out slowly with time. The morphology and structure of the precipitating phase are still in doubt, two different carbide phases, epsilon carbide and cementite, have been identified, but there is considerable disagreement about the time and temperature at which they appear.1-4 Epsilon carbide seems to form at lower temperatures than cementite, although one group of investigators4 found only the latter phase over a wide range of aging temperatures. Usually these two phases have been identified after aging times much longer than needed to drain the carbon from supersaturated solution, at least for aging temperatures below 200oC. The precipitated carbides have been observed with the electron microscope,&apos;-5 but again most results have come from samples aged much longer than necessary to reduce the carbon concentration in solution to a low value. In this study the carbide phases were observed in earlier stages of precipitation with the electron microscope, in order to learn more about the shape, size, and distribution of the carbide particles. Samples strained small amounts before aging, as well as unstrained samples, were studied, and the structure of the carbide phase was tentatively identified. EXPERIMENTAL "Ferrovac E" vacuum-melted iron, from which the wires for kinetic studies were also made (Refs. 6 and 7 and the following paper) was used for this work. An analysis of this material is given in Table I; it agrees well with the manufacturer&apos;s values. The number of inclusions in the iron was very small, as can be seen in the optical micrograph in the following paper. Pieces of iron about 1 mm thick and 5 to 10 mm in diam were used for all the studies except for those in which the sample was strained before aging; then the specimen was a rod 3 to 6 mm in diam and 20 cm long. The samples were held at about 700°C in wet hydrogen for at least an hour to remove residual carbon and nitrogen, after which they were heated to 950" (in the y region) and cooled to room temperature in the furnace. This treatment gave an average grain diam of more than 1 mm, so that the effect of grain boundaries on the precipitation process was negligible. The samples were carburized at 710°C for about 2 hr in hydrogen that had passed over toluene held at -30°C; then they were rapidly quenched into water. From chemical analysis this treatment gave a carbon concentration of from 0.012 to 0.016 wt pct, which is less than the solubility of carbon at this temperature. After quenching, the samples were aged at 120 "C in an oil bath for vari-

Jan 1, 1961

By M. P. Tixier, M. Gondouin, G. L. Simard

In the quantitative interpretation of the SP logs, the electrochemical component is generally taken equal to — K log Rm /Rw where K has the theoretical value corresponding to solutions of pure sodium chloride. This method may be misleading when relatively large quantities of salts other than NaCl are present in form-rrtion waters—as is generally the case for low salinities —or in gyp-base muds. In such cases recourse is made in the field to changing the K value or to adding a correction term to the equation on the basis of local experience. An investigation has been made in the laboratory of the influence of HCO3, SO4=, Ca++, and Mg++ on the amplitude of the SP deflection, which included in par-ricular determinations of the activity coefficients of Ca++ and Mg++. The theory, the experimental techniques and a tentative method for applying the results are described. So far, applications to actual field cases where chemical analyses of waters were available have provided an excellent confirmation of the proposed method. The case of the SP in very salty brines is also considered and it is shown that the use of activity data for sodium gives satisfactory results. HISTORICAL In early times of electrical logging, the SP curve was used exclusively as a tool for the location of permeable beds and the definition of their boundaries. Later, with the introduction of methods of quantitative analysis, attention was called to the possibility of deriving from the SP log some information on the formation water resistivity, which constitutes one essential element for the computation of water saturation from log data. The first laboratory experiments&apos; in 1932 had shown that, as a result of electrochemical phenomena arising at the contact of the niud and the formation, the amplitude of the deflection of the SP curve was, at least partly, a function of the salt concentration of interstitial water and accordingly of its resistivity. It was indicated that the amplitude of the electrochemical component of the SP should be equal to K log R1/R2, R1 and R2 being the resistivities of the mud and of the formation water respectively. However, the exact value of the K coefficient was not determined at that time because the surface clay used for a laboratory model was far from an ideal shale membrane. Consequently, the way was not yet open towards quantitative interpretation. In the early forties, Tixier reported2, on the basis of field log observations, the following equation for the case of NaCl solutions, SP = -71 log f1c1/f2c2, c1 and c2 being the concentrations and f1, f2 the mean activity coefficients. The local variations observed in the empirical values of K, chiefly in the Rocky Mountains, led him to advance the opinion that, besides the relative amounts, the nature of the salts present in the solutions had a bearing on the magnitude of the SP. At about the same time Mounce and Rust- ad published laboratory results which emphasized the part played by the shale in the generation of the SP. On the basis of laboratory tests and theoretical considerations, Wyllie in 1948&apos; published the following equation of the electrochemical SP for the case of NaCl solutions, SP= - 71 log (a1)m/(a2)m

Jan 1, 1958