Search Documents

By O. C. Holbrook, George G. Bernard

A new theoretical treatment has been obtained for the behavior of pattern waterflood injection wells when closed in. Two cases are treated: Case I where oil and water are assumed to have the same properties, and Case 2 where they arc different. In applying the method, one plots log (p — p,) vs closed-in time, where p is well-bore pressure at any tims and p, is static pressure. The value of p. is determined by trial and error as that value which makes the plot linear at large time. A value for the permeability-thickness product can be determined from the intercept of this linear part, and a value of the skin factor from the injection pressure at time of closing in. Application of the method to data from water floods at three fields seems to give reasonable results. For the case of unit mobility ratio, it is proved that this new method should give the same value for permeability-thickness product as the conventional pressure build-up method. In addition, the new method gives correct values for static pressure, whereas the conventional method does not, often indicating negative static pressures. The new method may be used in cases where the surface pressure persists after closing in as well as in cases where it does not. INTRODUCTION It is of considerable interest and importance to be able to determine the characteristics of the reservoir in an area surrounding a water injection well. Thus, if we can determine early in the life of an injection well that there is a considerable "skin effect", remedial measures can be started before a full-scale pattern flood begins. Similarly, if it can be shown that a gradual buildup of skin effect is occurring with time, measures to free the water of plugging material can be taken. Determination of static pressure in the water-injection well may show that the water is entering a thief zone and not the desired reservoir. Finally, determination of the permeability of the sand around the injection well will allow estimation of the future relation between injection pressure and rate. It should be possible to determine average reservoir permeability, skin effect and static pressure from pressure fall-off data. However, at the time we began work on this subject, it was thought that no adequate theory on which to base such determinations' was available. According to the conventional method which considers the reservoir to be filled with one fluid of small compressibility (see Van Everdingen, Joers2, and Nowak2), shut-in pressure is plotted vs log where is injection time, and At is closed-in time. The physical significance of injection time, may well be questioned in this case, since in a reservoir completely filled with a single fluid (as required by this theory) and with input and output rates equal, the pressure behavior after an initial transient is independent of t,. Attempts by our Tulsa area to use this theory led to negative values of static pressure in most cases. Because of these limitations of the method discussed above, it was decided to attempt to develop a new theory of pressure behavior in water injection wells, one which would apply when there is a gas saturation, as is so often the case in water floods. In the following treatment the assumptions and basic equations are given first, then the method of application of the equations. A complete example is given to clarify details of application. All difficult mathematics has been placed in the appendices so that the reader can follow the text without difficulty. However, if he wishes only to apply the results without knowing the basis for them, he can learn how to do this from reading only the sections entitled "Plotting of Experimental Results" and "Example." ASSUMPTIONS AND BASIC EQUATIONS Statement of Problem It will be assumed that a horizontal layer of constant thickness contains in its pore system a mixture of oil, gas and water. While water is being injected into this pore- system through a well at constant rate, an oil bank is built up, gas being expelled from the space taken by the oil as shown in Fig. 1. The saturations within each

By W. M. Robertson

The kinetics of grain boundary groove formation on copper surfaces immersed in liquid lead have been studied over the temperature range of 400° to 900°C. The groove widths were Proportional to the cube root of the annealing time, indicating that the diffusion of copper through the saturated liquid-lead solution is the mechanism by which grooves form. The rate constant for the process can be calculated by a theory due to Mullins. Using published data for the solubility of copper in liquid lead, for the interfacial energy of the solid copper-liquid lead interface, and for the diffusion coefficient of copper in liquid lead, the calculated rate constants almost exactly reproduce the measured values over the entire temperature range. The decay of scratches placed on copper surfaces has also been studied. The rate constant for scratch decay agrees with that for grain boundary groove growth, though scratch decay is sensitive to convection currents in the liquid. The interaction of solid and liquid metals has been studied extensively and has been found to be rather complex. Many of the studies have been complicated by a host of competing effects, including temperature gradients, concentration gradients, the formation of one or more intermetallic phases, natural or forced convection, and the presence in both the solid and liquid of considerable amounts of impurities and alloying elements. Recently theories have been developed of changes in surface profiles driven by capillary forces,' which allow the interaction of the solid and liquid to be analyzed quite directly. The present paper reports a study of the formation of grain boundary grooves on pure copper immersed in liquid lead, which allows the theory to be checked quantitatively. The Cu-Pb equilibrium diagram is quite simple: there are no intermetallic compounds, copper is just slightly soluble in liquid lead, and lead is almost completely insoluble in solid copper.2 The Cu-Pb system has been studied extensively,3-9 and all of the quantities—interfacial energy, diffusion coefficients, and equilibrium solubilities-necessary to calculate the rate of groove growth by a volume-diffusion mechanism are known. The kinetics of grain boundary grooving and scratch smoothing have been studied on the noble metals and the iron group metals in a reducing atmosphere or in vacuum.10 The main interest of these studies has been the determination of surface-diffusion coefficients at metal-gas interfaces using the theories developed by Mullins.1 The present paper uses similar methods for the study of diffusion processes at the interface between two condensed phases. There has been one previous study" of the kinetics of groove growth in a solid-liquid system, which, however, was troubled by the fact that neither interface energies nor diffusion coefficients were known in the Ni-S system studied. THEORY The theory of grain boundary groove growth has been developed by Mullins for the cases of volume diffusion," surface diffusion," and a combination of the two processes.14 He has also analyzed the kinetics of the flattening of a scratched surface by these processes.l5,l6 Several assumptions and approximations were used in the derivations of the surface profile shapes. The applicability of these assumptions and approximations to the present system will be considered in the discussion. At the intersection of a grain boundary with an initially flat interface the equilibrium between the interfacial energy and the grain boundary energy establishes an equilibrium groove angle. This induces curvatures in the interface. The chemical potential of material at a curved interface is higher than at a flat interface, so that material moves away from the region of the grain boundary. Mullins' calculations indicated that the groove profile would be similar to that shown schematically in Fig. 1. Small humps form above the level of the original flat surface. For a volume-diffusion mechanism the separation, w, of these humps at a time, t, is given by and where co is the equilibrium concentration of the solid in the liquid, y, the solid-liquid interfacial free energy, O the volume occupied by a solute atom in the liquid, D the diffusion coefficient of the solid in the liquid, and kT has its usual meaning. For a surface-diffusion mechanism the groove width is proportional to the fourth root of the time. For grooves forming by a combined surface and volume-diffusion mechanism the exponent of the

Jan 1, 1965

By G. W. Preckshot, J. W. Gorman

MOLTEN metals offer an excellent medium for the study of molecular diffusion and interphase transfer. In the absence of intermetallic compound formation, solutions of molten metals are solutions of elements. Therefore, complicating factors such as molecular structure, dissociation, association, polarization, and so on, are generally absent. Further, metals remain in the liquid state over a wide range of temperature. Thus, the effect of temperature on transport processes can be studied in more detail. Although some data are available for the diffusion of metals in the molten state at high temperatures,"""".'!' the most useful and accurate data are those for the diffusion of various elements in mercury at room temperature. The effective radius of diffusion for the various elements in mercury may be calculated from the Stokes-Einstein relation.'-These calculated radii are compared with the atomic and ionic radii calculated from X-ray data on crystals." The effective radii of diffusion of the alkali and alkali-earth metals correspond to those of the neutral atom; those of cadmium, silver, lead, zinc, thallium, bismuth, and tin correspond to ionic radii; and those of copper and gold indicate larger radii due to association with mercury or more likely, faulty diffusion data. Independent evidence"' supports the notion that alkali metals are uncharged in amalgams. The concept of metals moving free of their conductance electrons in the liquid state has been discussed by Eyring, et al,'" in connection with studies on the viscosity of pure molten metals. Considering the approximate nature of the tabulated atomic and ionic radii and the variance in the diffusion measurements, the application of the Stokes-Einstein equation to the diffusion of metals in mercury gives gratifying results. It is noted that a number of these metals have article radii eaual to or less than those of the solvent mercury.".'*,"' Soba1,'" sing the inaccurate Roberts-Austin20 data, indicated that the effective radius of diffusion might well be a function of temperature. Additionally in studies of mass transfer between phases, many conflicting opinions have been expressed concerning the assumption of equilibrium at the interface." ,","," The results of the present work on the diffusion of solid copper in molten lead throws light on the tempera- ture variation of the diffusion radius as well as the extent of the interfacial resistance. The theoretical background and experimental procedure are presented here in detail. Results on other systems will be presented in subsequent papers. Nomenclature The following nomenclature is used in the paper: A—cross-sectional area of diffusion column, sq cm C, C,—concentration of solution and saturated solution respectively, g per cc D,D'—-diffusivity, sq cm per sec e—Naperian base h—ratio, k/D k—mass-transfer coefficient, cm per sec L—diffusion column length, cm n—series index q—flux, g per sq cm-sec r—radius, cm or A R—-gas constant t—time, sec T—temperature, K V—volume, cc w—weight of copper-lead alloy on copper cylinder WT—-total copper transferred to lead diffusion, column, g Z—distance from top of diffusion column, cm 2'—distance from copper-lead interface, cm (2' = L-2) u,,—roots of diffusion equation Mathematical Apparatus The model shown in Fig. 1 applies to the diffusion of solid copper into molten lead. Here Fick's first law may be written where all terms are defined in the nomenclature. It is assumed that the volume of the lead column is constant. This is justified since experimental densities compared favorably with calculated densities based on perfect interstitial mixing. Thus there is no counter diffusion of lead. It is clear that 2 in Equation [I] should be measured relative to the interface which moves relative to the diffusion cell but not to the liquid phase. Hence, the frame of reference is the plane of the copper-lead interface. The diffusivity calculated at constant volume, D", is that defined by Hartly and Crank" and is designated herein as D. No Interfacial Resistance—The two situations of the presence and absence of interfacial resistance will be treated. The simplest is the latter which re-

Jan 1, 1959

By E. A. von Hahn, T. R. lngraham

The rates of cementation of pd11 on electropolislzed copper cylinders were studied in aqueous perchloric acid solutions at pdII concentrations from 0.02 to 0.1 mM. The cylinders were rotated at various high speeds to reduce the diffusion layer to a minimum thickness. In 0.1 M HClO4 solutions, the cemented palladium is dense, adherent, and shiny. The rate data indicate that there are two stages in the deposition. In the initial stage, the rate of cementation is first order with respect to the pd11 concentration. The second stage is much slower. The first stage is consistent with rate control by the diffusion of pd11 ions to the copper surface, and/or chemical reaction at the surface, and the second stage is consistent with rate control by the diffusion of copper ions from the cop-per surface through the deposit of cemented palladium, out to the main body of solution. In 0.001 M HClO4 solution, only the first stage is evident and the rate is more rapid than at higher acidities. This rate enhancement is attributed to PdOH+ ions that predominate at low acidity and aye more reactive than unhydrolyzed pdII. The deposit is porous and loose. All of the cemented deposits are Pd-Cu alloys rather than pure palladium. Activation energies are 9.5 kcal per mole in 0.1 M HClO4 and 7.4 kcal per mole in 0.001 M HClO4. CEMENTATION or displacement reactions occur between aqueous solutions of metal salts and immersed metals according to the general equation: where M1 is electrochemically the more noble metal. Although these reactions find considerable application in metals recovery processes (e.g., the cementation of copper192 or gold3), in electrorefining processes (e.g., solution purification before electrolysis4,5), and in plating and metal finishing processes,6 few studies have been made of the kinetics of such reactions.7-9 The rates of cementation reactions will depend on one or more of the following factors: a) chemical reactions at the metal-solution interface; b) ionic transport to or away from the interface;'-' and c) the character and adherence of the cemented deposit on the substrate metal,6a,10 because the deposit will inhibit the rate of transport of dissolving ions (M2m+) into the solution. In this investigation, the kinetics of the early stages of cementation were studied to obtain an understanding of the reaction mechanism under conditions in which the deposit was sufficiently thin that its inhibiting effect could be disregarded. To achieve this, very dilute solutions in Mn+ were used. Cementation rates are often fast in the initial stages and transport-contr01led.9 To find conditions under which transport control would shift to chemical control, very high stirring velocities were used to minimize the diffusion-layer thickness. Of the many possible cementation reaction systems, the palladous perchlorate-copper system was chosen for this investigation because it was believed to involve simple ion-for-ion exchange. In addition, there are no interfering side reactions, such as the reduction of hydrogen ions, and anion effects are usually absent in perchlorate solutions. EXPERIMENTAL Materials. Reagent-grade chemicals and redistilled water were used in all experiments. Palladous perchlorate stock solutions (0.02 M) were prepared by dissolving 99.5 pct Pd sponge (Johnson Matthey and Mallory Ltd.) in concentrated nitric acid, adding concentrated perchloric acid, evaporating twice slowly (to prevent PdO precipitation) to a small volume, and diluting in a volumetric flask. To prevent gradual hydrolysis of PdII, the stock solutions were made 0.4 M in HC104. Tests for pdIV (Ref. 11) and C1- were negative. Copper strips were cut accurately from 0.025-in. electrolytic tough pitch sheet, Sample 1 (Mines Branch stock or Canadian Copper Refineries Ltd., Montreal), or from oxygen-free, high-conductivity sheet, Sample 2 (American Metal Climax O.F.H.C. brand, sold by Utility Brass and Copper Corp., Brooklyn). The strips were annealed for 1 hr at 335°C in a stream of purified nitrogen (Union Carbide, Linde Division) which had

Jan 1, 1967

By Peter Soo

Pvestrained Ferrovac E iron has been neutron-irradiated at approximately 90°C to an integrated flux of 1020 nut (E > 0.82 mev]. The irradiation was found to produce an incveased temperature dependence of the flow stress in addition to a greatly increased athemal stress. Measurements of the flow stress and stress relaxation, from which the activation volume and activation energy for slip were deduced, show that neutron irradiation changes the rate -controlling slip process to one based on dislocation interactions with tetragonal distortions which are Produced around submicroscopic interstitial loops in the lattice. The study indicates that without prestraining prior to irradiation the chances of detecting a change in the rate -controlling slip process are greatly reduced because in the initial stages of slip a substantial fraction of the radiation defects are swept out of the slip plane by gliding dislocations. Thus, activation parameters which are subsequently measured are representative of a greatly reduced defect density and would not differ appreciably from those for unirradi-ated material. The large increase in the athermal component of the flow stress is probably connected with the presence of depleted zones in the lattice which are introduced by irradiation. ALTHOUGH fast neutron-irradiation has not been observed to markedly alter the activation parameters for slip in bcc metals,' small but significant changes do occur. Most experimenters agree that irradiation predominantly increases the athermal component of the yield stress.'-= In addition to this, Laidler and smidt7 have shown that in iron irradiated to 5 X 10" nvt and molybdenum irradiated to 10" nvt, changes occur in the activation volumes for slip. A similar conclusion has been reached by Milasin and Malkin8 for irradiated iron. Work by Ohr et a1.5 shows that for Ferrovac E iron, irradiated to 1.2 X 1016 nvt, small increases in the activation energy for slip also occur. So far these changes in the activation parameters have not been explained on a firm theoretical basis. One important factor which would minimize the chances of detecting a change in the slip mechanism upon irradiation is the presence of "channeling" which has been observed in molybdenum,9 niobium,10 and iron.11These channels are formed by gliding dislocations which sweep irradiation defects out of the active slip planes and thereby create zones in which continued dislocation motion is encouraged. The activation parameters for the dislocations gliding in the defect-free channels would, therefore, be similar to those for unirradiated iron and a change in the rate-controlling slip process would be difficult to detect. In the present work, an attempt has been made to reduce the effect of uneven deformation on the measured activation parameters for slip in neutron-irradiated Ferrovac E iron polycrystals, so that a more realistic assessment of the effects of neutron-irradiation could be made. Primarily, the experiments involve the irradiation of specimens which had been prestrained to 9 pct elongation at room temperature prior to insertion into the reactor. It was hoped that the introduction of a large number of evenly distributed dislocations would substantially decrease any channeling effect which might otherwise occur. MATERIAL AND EXPERIMENTAL PROCEDURE The starting material was vacuum-melted Ferrovac E iron, an analysis of which is given in Table I. The standard tensile specimen had a gage length of 1.125 in., a cross-sectional diameter of 0.120 in., and a re-crystallized grain size of 1.2 x 10-3 in. All tensile tests were conducted on a floor model "Instron" tensile machine at a strain rate of 3 x 10-4 per sec. The irradiation of the prestrained specimens was performed in the Brookhaven High Flux Beam Reactor to an integrated flux of 1020 nvt (E > 0.82 mev) at a temperature of about 90°C. All specimens were excap-sulated in high-purity aluminum sheaths which were lightly swaged around the samples to ensure good thermal contact. Subsequent measurements on the irradiated specimens showed that within experimental accuracy the swaging had not deformed them. EXPERIMENTAL RESULTS Fig. 1 shows the flow stresses for a series of unirradiated control samples. In order to produce a comparable dislocation substructure throughout the test sm range, all specimens were prestrained

Jan 1, 1970

By Donald B. Evans, Robert D. Pehlke

The solubility of nitrogen in liquid Fe-B alloys has been measured up to the solubility limit for the formation of boron nitride. The activity coefficient of nitrogen increases with increasing boron content in the range 0 to 7 wt pct B. From experimen -tal data, values have been calculated for the B-N interaction parameter e3 at temperatures in the range 1550" to 1750°C. A value of 0.038 has been estimated for the boron self-interaction parameter eg at 1550°C. The standard free energy of decomposition of boron nitride into the elements dissolved in liquid iron has been determined to be: ?F° = 45,900 -21.25T in the range from 1550° to 1750°C. The nitride is assumed to be of composition BN. BORON nitride has an unusual combination of properties which make it appear attractive in a wide range of engineering applications. Some of its more important and most recent applications are in the nuclear area, particularly in connection with the liquid-metal cooled reactor concept now receiving considerable emphasis. Boron nitride has a high degree of stability at elevated temperatures. It also has excellent ma-chinability and the ease with which its crystals deform suggests applications as a lubricant. These properties stem from a hexagonal layer-type structure similar to the structure of graphite. One of its primary uses to date has been for seals in liquid-metal pumping systems. It is also used in nuclear reactors as an insulating layer to separate two solid metals which are not themselves compatible under the conditions of temperature and atmosphere in which they are used. Its inertness to liquid metals has also suggested use as a mold-release agent in casting processes. In addition to its excellent machinability and reported inertness to liquid metals such as iron, silicon, aluminum, copper, and zinc, boron nitride has high thermal conductivity and excellent thermal shock resistance. This combination of properties would make it appear ideal as a refractory crucible material for refining of high-purity liquid metals, for example high-quality steels. However, since it is known that concentrations of boron as low as 50 ppm can have a marked effect on the physical properties of certain steels,' in particular on the creep and stress-rupture properties, an investigation was undertaken to define accurately the chemical equilibrium among boron, nitrogen, and liquid iron in the range of steelmaking temperatures. EXPERIMENTAL PROCEDURE Two experimental approaches to this problem were employed: a Sieverts' method and a quenching method. In the first method, the Sieverts' technique was used to measure the equilibrium nitrogen solubility in liquid Fe-B alloys of 0 to 7 pct B as a function of nitrogen gas pressure over the melt. The solubility limit of the boron nitride phase formed was determined by the point of departure of the nitrogen absorption from Sieverts' Law. This technique has been applied to liquid Fe-Ti alloys by Rao and parlee,' to liquid Fe-A1 alloys by Evans and Pehlke,3 and to solid Fe-V alloys by Fountain and Chipman.4 In the second method a melt of liquid iron was held in a crucible of boron nitride under a known partial pressure of nitrogen gas. After thermodynamic equilibrium was attained, the melt was quenched in a stream of helium and then analyzed by wet-chemical methods for boron and nitrogen. The Sieverts' apparatus used in the first method was essentially of the same design as the one described by Pehlke and Elliott.5 The charge materials were vacuum-melted high-purity iron (Ferro-vac E) supplied by the Crucible Steel Co. and -325 mesh boron powder supplied by Cooper Metallurgical Associates of Cleveland, Ohio. The boron contained less than 0.02 wt pct O, according to supplier's analysis. Recrystallized alumina crucibles were used to contain the melt. Examination of solidified melts showed these crucibles to be satisfactory with no evidence of any reaction or physical penetration of the crucible wall by the melt. The melt temperature was measured by a disappearing filament-type optical pyrometer sighted vertically downward on the melt surface through a 1/4-in.-diam sight hole in the crucible lid. The pyrometer was calibrated against the melting point of pure iron in the same apparatus, taking the emissivity of

Jan 1, 1964

By S. R. Titley

The recent emphasis placed upon the use of molybdenum as a geochemical indicator has stimulated considerable inquiry into the behavior of molybdenum in the zone of oxidation. This paper represents a summary of recent and continuing experimental work and calculations in an investigation of a portion of the system Mo-Fe-S-H20 and treats the stability of various iron and molybdenum compounds with respect to oxidation potential and pH. Using experimentally determined free energy values for ferrimolybdite, diagrams of the stability fields of this mineral have been calculated with respect to ilsemannite (?), molybdenite, pyrite and limonite. This theoretical approach is the one utilized by Garrels3'4 in studies of the relationships of various minerals in the zone of oxidation. Use of this approach is particularly suited to the study of geochemistry in the zone of supergene processes because of the low temperatures, low pressures and other relatively easily determined and duplicated factors that exist in the near surface environment. Such a study of molybdenum, however, is complicated by several factors. One of these is the fact that the strong amphoteric behavior of molybdenum makes possible a wide variety of substances and complexes whose exact chemical natures are unknown. The second factor is the dearth of thermodynamic data for many of the simplest of the possible compounds and ions important to the study. For the purpose of this study, the problem of complexes has been approached through considering only the simpler possibilities realizing that the complexes and systems so treated may be gross oversimplifications. All the complexes considered, however, have been suggested in the chemical 2,7,9,12 or geologi-cal literature. EXPERIMENTAL WORK Studies of the behavior of molybdenum ions and ferrimolybdite have been carried out in aqueous sys- tems, and reactions have been followed by simultaneous measurement and recording of pH and oxidation potential. A major problem involved in this study was the long periods of time necessary for equilibrium to be established in some of the systems considered; in a few instances, these periods exceeded four weeks, that is, four weeks during which measured conditions changed continuously. Another problem, in some instances, was the indeterminate nature of the reactions that have taken place. Because of the tendency of molybdenum to complex, reactions which take place under some of the more extreme imposed conditions of oxidation and acidity or alkalinity are unknown. However, useful data for the present work were obtained from studies of the molybdate ion and ferrimolybdite between pH 2 and 8 and under only moderately oxidizing to nonoxidizing conditions. Instruments used in the experiments were a Hewlitt-Packard Model 412A VTVM, a modified Beckman Zeromatic pH Meter and two Varian model 11 recorders. pH was measured between glass and calomel electrodes and Eh was measured between calomel and platinum electrodes. Conditions of pH were varied with 1N H2SO4 and 1N NH4OH and 1N NaOH. Eh was varied with gaseous nitrogen, oxygen or aqueous H202. Ambient temperature was essentially uniform at 24ºC and was increased in the reaction vessel by the positioning of a magnetic stirrer whose temperature was sufficient to bring the temperature to 2S°C in the reaction. The temperatures of reaction did not vary more than 1/2ºfrom 25°C. Two reactions were investigated, the ferrimolybdite-limonite-water reaction and the reduction of molybdate ion. Studies of both reactions yielded useful data that were used in preparation of the Eh-pH diagrams. The traces of the reactions are shown on Fig. 1. Ferrimolybdite-Limonite- Water: Studies of this reaction involved analysis of the solution of ferrimolybdite and simultaneous formation of ferric hydroxide. The data are summarized on Fig. 1. In this reaction, equilibrium was attained at various pH's when coexistence of the two solid phases (ferrimolybdite and ferric hydroxide) was attained. At the proper pH, controlled by the amount of ferrimolybdite and thus the dissolved content of molybdate species, both phases were observed and reaction times to complete the transition were on the order of weeks. The reac-

Jan 1, 1963

By R. Schuhmann

Available thermodynamic data applicable to copper smelting systems are collected and tabulated, and the important gaps are pointed out. A few examples are given of estimations which can be made from the available data. An experimental research program is proposed to supply the thermodynamic data that appear most essential to better quantitative understanding of the chemistry of copper smelting. The proposed program is designed also to shed specific light on the practical problems of slag losses and magnetite behavior. OPPER smelting, from flotation concentrates to V-/ blister copper, is conspicuous among the large scale chemical processes which are conducted with only relatively incomplete knowledge of the physical chemistry involved. For example, no common metallurgy text explains adequately why copper enters the matte phase while iron enters only to the extent that sulphur is left over after satisfying the copper. Explaining this important phenomenon as a manifestation of the greater affinity of sulphur for copper than for iron is unsound because the affinities of copper and iron for sulphur are about the same at smelting temperatures. As will be shown, other affinities are really decisive in the relatively clear-cut separation of copper in the matte. In contrast, thermodynamic studies have contributed much to our understanding of copper refining, zinc oxide reduction, magnesium production, steelmaking, etc. For each of these processes, the important reactions are clearly recognized, and fair to good quantitative values of the free-energy changes and equilibrium constants are available. Such data have proved to be of constantly increasing practical value in process development and improvement. Work in other fields has furnished much thermodynamic data applicable to copper smelting systems. In fact, several promising starts have been made in applying such data to specific smelting problems. Kelleyl made an appraisal of the possibilities of recovering elementary sulphur from low grade matte. His calculations and compilations of thermodynamic data have represented the starting point for much of the work in this field, including the present survey. Huang and Hayward2 nd Aksoy3 used thermodynamic methods in the study of copper losses in reverberatory slags. Peretti4 used thermodynamic data in explaining the chemistry of converting. The recent publications of Darken and Gurry"." and of Darken,' dealing with the iron-oxygen and iron-silicon-oxygen systems, respectively, present equilibrium data that are applicable to copper smelting systems. Also additional data were reported recently on the affinity of sulphur for copper, manganese, and iron8 and on the sulphur pressures of iron-sulphur melts." The survey presented in this paper was made as the basis for planning an experimental program on the thermodynamics of copper smelting. Few researches on the chemistry of copper smelting have been reported in recent years, so that a reappraisal and coordination of old and new data are essential if further work is to take the directions of maximum value. The experimental program is now in progress, and the plan of attack is outlined at the end of this paper. Progressive Oxidation and Desulphurization of Copper-bearing Liquid Phases: The chemical activities of sulphur and oxygen are two of the most important thermodynamic yardsticks to be applied to copper smelting processes. Virtually the entire smelting and refining sequence involves a series of systems characterized by decreasing sulphur activity and increasing oxygen activity. In this section, therefore, an attempt is made to define and explain these activities in terms of equilibrium partial pressures of SO2, s2, and O,. Also, estimates of these quantities are presented, the estimates being based largely on calculations presented in a later section of this paper. The sequence of steps from raw flotation concentrate to fully oxidized copper ready for poling involves progressive and controlled oxidation. Iron is oxidized and enters the slag. Sulphur is oxidized and leaves in the gas. Table I summarizes several important features of this oxidation and desulphurization sequence, starting at the beginning of the matte blow in the converter. The top line gives in order the principal smelting and refining stages up to fully oxidized copper, plus the additional step, not used commercially, of oxidizing all the way to Cu,O. In the second line are shown the principal copper-bearing liquid phases which characterize the process. Through most of the sequence the copper

Jan 1, 1951

By D. J. Rose, W. M. Armstrong, A. C. D. Chaklader

The wetiability of single crystals of MgO by specimens of vacuum-cast iron was studied using the sessile drop technique in vacuo at 1550ºC. Formation of FeO at the liquid-vapor interface caused the contact angle (6) to decrease from 117 to 65 deg during the first minute. After cooling, all specimens possessed a peculiar annular interfacial deposit of Fe,O,. Within the annulus the interface showed no sign of chemical attack. Chemical reaction occurred where the iron was alloyed with Ti, V, Cr, Nb(Cb), Ta, cmd Zr. Vnriation of 6 with alloy concentration was studied. Although vanadium md chromium improved the wettability of MgO by iron, the effect of Zr, Ti, Nb, and Ta was indeterminable because the 8 derived from sessile drop considerations was that of the metal against a restrictive peripheral volume of liquid oxide wetting the substrate. INTEREST in the nature of metal-ceramic interactions has been stimulated by progressive development of metal-ceramic combinations. One of the more valuable methods used for bond investigation is the sessile drop technique. Recent attempts to improve wettability through solute additions to the drop have revealed that the solute may a) react with the oxide creating new compounds at the solid-liquid interface, b) adsorb at the interface in a monolayer formation, or c) distribute throughout the drop uniformly causing no wettability variation. This work, part m in a series1,2 of investigations of interface reactions between metals and ceramics, embraces a study of the Fe-MgO system. MgO possesses a large negative free energy of formation (-173 kcal per g mole 0, at 1550ºC) compared to FeO at this temperature (-73 kcal per g mole 02). Hence the electronegativity difference between the drop and substrate allows selection from an extensive group of metal solutes with intermediate electronegativity differences that would, in the absence of chemical reactions, be expected to adsorb pre- ferentially at the solid-liquid interface. However, the high oxygen pressure of MgO in vacuo creates an oxidation environment which influences the solute behavior. Recent studies of the Fe-MgO systemJ74 have been restricted by chemical reactions that obscured observations. In view of the current interest in liquid-phase sintering of metal-ceramic combinations under partial oxidizing conditions,5 a more comprehensive study of the Fe-MgO system seemed beneficial. EXPERIMENTAL PROCEDURE 1) Specimen Preparation. Optical-grade MgO single crystals and Ferrovac vacuum-cast iron were used in all experiments. A spectrographic analysis of the iron indicated 0.008 C, 0.05 V, 0.005 Mo, 0.01 Ti, 0.002 S, 0.01 Cr, 0.001 Mg, 0.01 Si, 0.003 A1 (wt pct). The MgO crystals were cleaved along (100) planes into plates approximately 20 by 20 by 1 mm. Following annealing in vactto at 1100°C for 3 hr, surface irregularities were removed by a chemical etch in phosphoric acid at 100°C. The iron rod was machined into small cylinders approximately 0.250 in. in diam and length and these were carefully cleaned in organic and acidic solutions to remove surface impurities. All specimens were weighed to the nearest 0.1 mg. 2) Apparatus. The apparatus described in detail by previous authors1,2 was modified for this work. A Polaroid Land Camera was incorporated into the optical system so that drop measurements at ten-second intervals after melting were possible with ultra-high-speed self-develop ing film. 3) Procedure. The iron cylinders were placed upright on the MgO plates and positioned in the molybdenum susceptor. After furnace assembly the system was pumped down to a vacuum of 1 µ and flushed with H2 at 800°C. The system was again evacuated to 5 x 10-5 mm of mercury and the temperature raised to 1550°C within 2 min. The molten drops were photographed at appropriate time intervals. TWO min after melting the iron vapor pressure caused gas discharge, or corona, within the tube. The specimens were then slowly cooled to room temperature, sectioned and examined metallographi-cally. X-ray identification of interfacial reaction products was attempted by the powder diffraction technique. Alloying elements (Ti, V, Cr, Zr, Ta, Nb) were added to the iron in various concentrations up

Jan 1, 1963

By A. D. Rovig, D. W. McGlashan, Donald M. Podobnik

Crystal-chemical and stntctural properties of sulfide minerals are considered. The information gained is to be used to interpret (I ) freshly broken mineral surfaces, (2) modifications of the mineral surfaces, and (3) reactions at the mineral surfaces. From these basic disciplines, concepts with regard to the changes that surfaces undergo, reactions that might take place, the geometry of the interface, the state of different atoms and ions in the interface, and other physical and chemical properties of an interface must be developed, weighed and applied. The authors first deal with these conceptual considerations from which hypotheses are set forth to describe the environ men tal-interfacial relationships for several sulfide minerals. Qualitative and quantitative explanation of par-ticulate solid separations are unknown entities because of the lack of adequate models and mathematical relationships to explain the activity occurring between the solid and liquid phases. It is a transitional region in which it is difficult to ascertain the mechanisms of adsorption of ions, molecules, etc., on mineral surfaces, as well as other secondary reactions which occur in interfacial regions. Thus, in this paper the authors deal with conceptual considerations from which hypotheses are set forth to describe the environmental-interfacial relationships for sulfide minerals. GENERAL CONSIDERATIONS Interfacial Reactions:Interfacial reactions are the cruxes of flotation schemes as well as other processes such as thickening, filtration and hydrometallurgy. However, as noted by Klassen and Mokrousov: 1 "The problems concerning the surfaces of natural minerals, the laws governing simultaneous adsorption from aqueous solution on these surfaces of a whole series of reagents, the laws of the surface reactions and the properties of water layers separating the minerals, are all known to a first approximation only." An investigator has the privilege of postulating methods of solid-liquid interfacial reactions. REACTIONS WITH WATER - Water consists of hydronium (H3 O 4) and hydroxyl (OH-) ions in the ionized state. That this is so forms the basis for the postulated reaction of the hydrated hydrogen ion with net-negative mineral surfaces as depicted in Fig. 1. In this case water is bonded to mineral surfaces through hydrogen-bridge-type bonds - possibly hydrogen bonds. Although bonding of the van der Waal type may be responsible for this reaction, it is most likely that stronger bonds are involved. Once firmly bonded to the surface, the water layer is known to be quite tenacious. It is further speculated that a shear plane exists at some distance (see Fig. 1) away from the mineral surface. Exact location of this plane is not known, but in all probability it will be positioned across a weak bond of the hydrogen-bond type where the surface-attached water is coordinated to the original hydrated hydrogen ion. REACTION WITH METAL IONS - Klassen and Mokrousov ' state: "The presence of an ion in water leads to an immediate formation around that ion of a highly condensed atmosphere of water dipoles and thus to hydration of the ion." The fact is known that multivalent cations are strongly hydrated, usually by six or eight molecules of water, and that anions are not so strongly hydrated. Conceding the fact that thermodynamics, concentration, pH, and physical considerations are of utmost importance in truly explaining a mechanism of metal ions reaction with a hydrated mineral surface, it does seem logical that the mechanism illustrated in Fig. 2 is feasible. In this reaction, the metal ion (M) is coordinated in one dimension - to the mineral-surface hydration layer. Note also that explanation of this reaction requires that the shear plane move to a less stable bond configuration; that is, the shear plane has moved to a position between the metal ion and the other coordinated water molecules which are more free to dissociate. Reactions as depicted in Fig. 2 should cause the formation of an apparent mineral surface which is net-positive. Immediately, it must be noted that measurements of apparent-mineral surfaces serve - within limits - to indicate a degree of ion-surface reaction capability. The major limiting factor he re is one related

Jan 1, 1970

By R. E. Zimmerman, Michael Sokaski, E. R. Palowitch, M. R. Geer, H. F. Yancey, W. Deurbrouck, S. C. Sun

PART 1: DENSE MEDIUM SEPARATION by M. SOKASKI, M. R. GEER and H. F. YANCEY INTRODUCTION In the early days of coarse-coal treatment by the dense-medium process in Europe, loess was one of the materials used for medium solids. Newly developed cyclones were used to reclaim and thicken :he loess suspension. On occasion, when the cyclone thickener at the Maurits mine in the Netherlands plugged, the overflow was found to be filled with clean coal free of impurity. From this observation stemmed the concept of the cyclone as a cleaning device. Development took place during World War 11 in the laboratory of the Dutch State Mines under the direction of M. G. Driessen. A 15-ton-per-hour pilot plant was in operation in the laboratory by 1945, when, in September, Driessen made the first public disclosure of the process in a paper presented to the Institute of Fuel in London. In April 1945,2 operation of this pilot plant was observed by one of the authors, who was so impressed with the potential of the cyclone process that, on his return to the United States, a research program on the cyclone was started by the Bureau of Mines immediately. This work led to an AIME paper in 1946 in which the first details of cyclone operation on American coals was made available. Commercial adoption of the cyclone proceeded in Europe, initially in the Netherlands, Germany, and France. In the United States, where coal is less costly and most of the coals are so easy to clean they respond well to simpler cleaning systems, the first cyclone plant was not built until 1961.3 As of the spring of 1967, a total of 31 commercial cy- clone plants had been installed or were under construction in the United States; this does not include eight Dynawhirlpool plants. The capacity of the cyclone plants ranges from 50 to 1300 tons per hour and totals 6600. THEORY Before discussing the theory of a dense medium cyclone, a brief description may be helpful to some readers who are not familiar with its basic operation. In a typical dense-medium cyclone, illustrated in figure 10-1, the mixture of medium and raw coal enters tangentially near the top of the cylindrical section, thus forming a strong vertical flow. The refuse moves along the wall of the cyclone and is discharged through the under- flow orifice. The washed coal moves toward the longitudinal axis of the cyclone and passes through the vortex finder to the overflow chamber.

Jan 1, 1968

By Francis Collins, E. R. Brownscombe

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By G. W. Ostroot, S. Shryock

Cementing deep, high-temperature oil wells where static temperatures range from 350 to 400F has become routine in the part decade. In the United States there were 271 wells drilled deeper than 15,000 ft during 1963. Many of these wells had static temperatures higher than 400F. Bottom-hole static temperatures near 700F are now realities in the geother-mal (steam producing) wells of California's Salton Sea area. The detailed planning initiated prior to drilling the wells is discussed together with the methods, materials and equipment used in solving the cementing problems which are encountered. Data are also presented that lead to development of cementing compositions that provide adequate thickening time, do not retrogress in strength, and maintain low permeability under these extreme temperature conditions. Field data include the cementing programs used on eight relatively trouble-free geothermal steam wells in the Salton Sea area. INTRODUCTION Not too many years ago cementing oil wells with temperatures in the range of 300F caused considerable anxiety. In some areas of the United States it is now fairly common to cement wells having bottom-hole static temperatures in excess of 400F. We are now confronted with the problem of cementing wells with temperatures ranging from 500 to 700F. Temperatures in this order of magnitude are often found in geothermal steam wells. From where does this extreme heat emanate? There are many theories as to the source of this steam flow. The most widely held views are: (1) heat- ing of ground water fairly close to the surface by an intrusive mass of hot rock; (2) steam generation from a reservoir of metamorphic rock, normally found below 25,000 ft and not at the shallower depths of the Salton Sea reservoir; and (3) high-temperature gases (water vapor) escaping and migrating from molten or semi-molten rock (magma) at a considerable depth. Of these. No. 3 seems to be the most generally acceptable explanation. Heat from springs and fumaroles has been used for years as a means of heating and cooking; however, significant progress in harnessing the vast power of underground steam reservoirs has been relatively slow. The first large-scale attempt to use the heat generated by steam from wells was made in Italy around the beginning of the 20th century. In excess of 250,000 kw of electrical power is now being produced from holes around Larde-rello, Italy. Another very active drilling program was initiated in the volcanic area of New Zealand in 1949.' Natural steam for power projects in the United States began in the early 1920's. An early commercial steam field is located at the Geysers, approximately 75 miles north of San Francisco, an area discovered in 1847 and used for many years as a health resort. Steam originates from 15 wells that have been drilled since 1957. The present output from this project is 25,000 kw. Success of the Geysers operation has been responsible for several companies taking a careful look at the feasibility of producing steam for power generation in the Salton Sea area of Southern California's Imperial Valley. Geothermal steam activity in this latter area began in 1961 when O'Neill, Ashmun and Hilliard completed Sportsman No. 1, at that time the hottest wellbore in the world.' Since its References given at end of paver. completion seven additional wells have been successfully completed in this area. Many problems encountered in drilling steam wells had to be overcome to make the ventures successful. Formation temperatures encountered in the Salton Sea seemed to be a straight-line function (a gradient of 13F per 100 ft of depth).' This imposed severe conditions on all aspects of drilling and completion. This varied, to some extent, from gradients in the older geothermal areas. Not to be overlooked is the effect of these temperatures on casing creep or elongation by thermal expansion (Table I), because standard API flanged wellhead equipment makes no provision for this kind of performance. Floating equipment was redesigned, and changes in types of downhole equipment were made in an effort to eliminate anticipated problems. In the later completed wells, standard oil-well cementing equipment has been used. During the early development of geothermal steam wells there were some problems resulting from blowouts. However, these were eliminated in the deeper Salton Sea wells and no problems were encountered with the drilling mud. A sodium surfactant mud was used on the Sportsman No. 1 to drill from 2,690 to total depth. Nevertheless, it was necessary that a cooling system be added and the mud cooled before circulating it back into the well. The difficulty in evaluating the size of the steam area and its potential in terms of pounds of steam and years of productivity still has not been resolved. Economic complexities have also entered into the steam play in the Salton Sea. The wells at the Geysers were drilled at a cost of $15,000 to $20,-000, whereas the Salton Sea wells will cost more than $150,000. This cost differential has to some extent been accounted for because of the heavily

Jan 1, 1965

By R. E. Smallman, I. Milne

The deformation behavior of single crystals of Nb-Mo alloys has been investigated with particular reference to the influence of composition, orientation, and temperature. Strong solid-solution hardening was observed reaching a maximum at the equiatomic cotrlposition and can be attributed to the difference in atomic size between niobium and molybdenutrz. Changes in the form of stress-strain curve, as shown by a high work-hardening rate and restricted elongation to fracture, were observed at a composition of Nb-85 pct Mo and are attributed to the presence of MozC DreciDitate. Conjugate slip was only extensive in dilute alloy samples; at the 50/50 composition deformation rnainly occurred by primary slip, and the onset of conjugate slip gave rise to failure by cleavage on (100). The variation of yield stress of Nb-50 pet Mo with orientation was consistent with slip on (011)(111) slip systems. The temperature deperndence of the yield stress between -196" and 250°C was similar to that of pure bcc metals, but at a much higher stress level; no evidence for twinning %as found. IN recent years the deformation behavior of various pure metals in groups VA and VIA has received considerable attention, but surprisingly little work has been carried out on binary alloys made by mixing metals from the two groups. Such an investigation would be of interest since single crystals of metals of group VA have been shown to deform characteristically with a multistage deformation curve1"3 while a parabolic type of deformation curve has been reported for most of the group VIA metals.4'5 It has been suggested by Law ley and Gaigher~ that the difficulty encountered in obtaining multistage deformation curves for molybdenum in group VIA was possibly because of the presence of a microprecipitate of MozC which they observed even at carbon contents as low as 11 ppm. Recently a multistage deformation curve has been reported for molybdenum ," although the stages are not so definitive as those for group VA metals. The binary alloys of the particular refractory metals which have been investigated in single-crystal form include Ta-w,' Ta- Mo,' and Nb- Na." While a large amount of hardening was observed for alloys of the Ta-W and Ta-Mo systems, associated with room-temperature brittleness for alloys approaching the equiatomic composition, Ta-Nb remained ductile over the complete composition range with little or no solution hardening. Other systems have been investigated by hardness measurements on polycrystalline material and a discussion of the hardening of these alloys has been presented by ~udman." The purpose of the present investigation was to examine the deformation behavior of Nb-Mo alloys in detail, with particular reference to alloy composition and single-crystal orientation. In this way it was hoped to shed some light upon the restricted ductility of these alloy specimens. 1) EXPERIMENTAL PROCEDURE The starting materials were obtained in the form of beam-melted niobium rod and sintered molybdenum rod of suitable dimensions. Since niobium and molybdenum form a complete solid-solution series at all temperatures, alloy single crystals were produced by melting the two constituents together in an electron bombardment furnace (EBM). To produce specimens free from segregation a molten zone was passed over the length of each rod six times in alternate directions at a speed of 10 in. per hr. Typical specimens were analyzed for interstitial impurities by gas analysis and for metallic impurities by spectrographic analysis. The results of this analysis are shown in Table I. Many of the tensile specimens were also analyzed (after testing) by scanning the gage length in an electron beam microanalyzer, from which it was found possible to predict the approximate composition of a specimen from the original proportions of each element in the EBM. The tensile specimens were made with a gage length of 0.5 in. and diameter of 0.075 in., using a Servomet Spark machine. By careful machining on the finest range for the final i hr of this technique, surface cracks could be reduced to the level where they were easily removed by electropolishing in a solution of nitric and hydrofluoric acids. The specimens were strained at a rate of 10 4 sec-' using friction grips designed to prevent accidental straining and maintain a good alignment before straining. The orientations of the individual specimens tested are shown in Fig. 1 and the corresponding compositions listed in Table I1 together with collated experimental data. 2)RESULTS a) General Deformation Behavior. The effect of composition on the room-temperature deformation curves of similarly oriented specimens is shown in Fig. 2. The yield stresses of the pure constituents, while not the lowest reported to date, were at least comparable with existing data. Although the solution hardening was large for alloys at either end of the phase diagram, and comparable with the Ta-W solution-hardening data of Ferris et a1.,8 the low work-hardening rate characteristic of niobium was sustained until a composition of Nb-85 pct MO had been reached. Associated with the peak yield stress ob-

Jan 1, 1969

By N. S. Stoloff, R. G. Davies

A study has been made of the efject oj different dislocation-precipitate interactions upon the temperature dependence of the flow stress of aged Ni-14 at. pct A1 alloy. It is observed that when the dislocations bow between widely spaced (-20004 coherent Ni3Al particles the flow stress decreases with increasing temperature in the normal way. However, when the dislocations cut closely spaced (-5004 particles the flow stress is independent of temperature from -100 to 600°C, due to a balance between softening of the matrix and an increase in strength of the particles with increasing temperature. The retention of strength at high tempera-tures of commercial nickel-base alloys, which are strengthened by the precipitation of a phase based upon Ni3Al, is thought to be due to the unusual strength properties of Ni3Al. The flow stress of Ni3Al increases continuous1y from -196"C to a maximum at -600"C. It is concluded from a series of thermal-mechanical tests that the sevenfold increase in flow stress over this temperature interval is due to a lattice effect and is not diffusion-controlled. The flow stress of precipitation- or dispersion-hardened materials depends on the resistance to dislocation motion within the matrix and the extra energy required for dislocations to bow between or to cut particles. If the dislocations bow between the particles or if the strength of the cut particles is constant with temperature, then the flow stress of the precipitation-hardened alloy must decrease with increasing temperature due at least to the decrease in elastic modulus of the material. There will be softening also from thermally activated cross-slip or climb, offering an additional degree of freedom for dislocations to avoid particles. For example, in the case of nickel containing a dispersion of thoria,' which most probably deforms by dislocations bowing between particles, the flow stress decreases by about 50 pct between 25" and 650°C. In A1-Cu alloys2 aged to produce the 8" precipitate, dislocations cut the particles, and the flow stress decreases by about 20 pct between -269" and 25°C. However, many commercial high-temperature nickel-base alloys, for example Inconel-X and Udimet-700, exhibit little or no decrease in flow stress with increasing temperature up to about 700°C. A characteristic feature of these alloys is that they are strengthened by the precipitation of a phase based upon Ni3A1. Guard and westbrook4 and flinn' have shown that Ni3Al (and alloys in which a third element such as molybdenum or iron is substituted for part of the aluminum) is unusual in that the hardness and flow stress increase with temperature to a maximum at about 600°C. For the flow stress of a precipitation-hardened alloy to be independent of temperature we propose that the particles must be cut by dislocations moving through the matrix and that the strength of the particle must increase with increasing temperature. Theories of precipitation hardening do not take into account the flow stress of the dispersed particles that are cut during deformation; the only dissipative process usually considered7 is the creation of interface within the particle and between the precipitate and matrix. The purpose of the present investigation has been to study in detail the temperature dependence of the flow stress of a nickel-base alloy strengthened by the precipitation of Ni3Al in two structural conditions such that when deformation occurs it does so by dislocations a) bowing between the particles and b) cutting the particles, respectively. A simple binary Ni-14 at. pct A1 alloy was chosen because considerable information is already available for this system concerning phase equilibria and precipitation reactions and rates.' Dislocation-precipitate interactions in the binary alloy should be similar to those in the more complex commercial alloys. In addition, the mechanical and physical properties of NisAl were studied in detail in the hope of elucidating the mechanism by which the strength increases with increasing temperature up to 600°C. EXPERIMENTAL PROCEDURE For the study of the effect of precipitation of Ni3A1 upon the temperature dependence of the flow stress, an alloy containing 14 at. pct A1 was utilized; a Ni-8 at. pct A1 solid-solution alloy was employed as a comparison material. Vacuum-cast ingots were hot-rolled at 1000°C and cylindrical compression samples, 0.20 in. diam by 0.40 in. high, were prepared from the 1/4-in.-diam rod. Specimens were recrystallized and solution-treated at 1000°C for 1/2 hr and then water-quenched. A preliminary study revealed that, when the Ni-14 at. pct A1 alloy was aged for 1 hr at 700°C, significant precipitation hardening was obtained, and that the structure was free from grain boundary discontinuous precipitation; an overaged condition was produced by annealing the aged specimens at 850°C for 1 hr. To circumvent the difficulties involved in the hot rolling and swaging of Ni3A1, compression samples,

Jan 1, 1965

By E. R. Brownscombe, Francis Collins

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By D. G. Westlake

Polycrystalline tensile specimens of various Zr-0-H alloys have been tested at 298°, 178°, and 77°K. Solute oxygen and hydride precipitates in quenched alloys made individual contributions to the yield strength at 0.2 pct strain which combined to produce a resultant strength increment, a,., Ductility changes which were ohserved can he interpreted in terms of the various oxygen and hydrogen concentrations, testing tem -peratures, and dispositions of the hydride. ADDITIONS of oxygen in solid solution were known to increase the yield and tensile strengths of polycrystalline zirconium as early as 1951.' More recently, the critical resolved shear stress (CRSS) for prism slip in zirconium single crystals was also shown to be affected by the solute oxygen impurity.' This latter work also demonstrated that large increments of strength could be contributed by the finely dispersed zirconium hydride precipitates that are present in quenched Zr-H alloys.3 It was concluded that the combined strengthening due to alloying could be expressed by where to is the increase in the CRSS due to solute oxygen alone and TH is the increase due to finely dispersed hydride precipitates. Eq. [I] is analogous to one used to express the combined strengthening effects of work hardening and neutron radiation damage.4 Eq. [1] was verified only indirectly and for only small amounts of the impurities—up to 0.14 at. pct 0 and 0.63 at. pct H. The present investigation was undertaken to obtain a more direct verification of the validity of the form of Eq. [1] for this system and also to determine the combined effects of oxygen and finely dispersed hydride precipitates on the tensile strength and ductility of polycrystalline zirconium. EXPERIMENTAL PROCEDURE Tensile specimens were machined from the same rolled billet of Kroll zirconium used in the earlier study.' These measured 38 by 4.7 by 0.5 mm and had 10-mm gage lengths which were 2.8 by 0.5 mm. Each specimen was ß-annealed in vacuo at 1173°K for 15.5 hr and a-annealed at 1073°K for 4 hr to D. G. WESTLAKE, Member AIME, is Associate Metal l ur-gist, Metallurgy Division, Argonne National Laboratory, Argonne, III. Manuscript submitted July 17, 1964. IMD______________ give an equiaxed structure with grain diameters averaging 0.06 mm. Oxygen was added by allowing the metal to react with a known quantity of oxygen during the 0 anneal and known quantities of hydrogen were added during the a anneal. Each alloy was encapsulated in Pyrex under vacuum, annealed at 873°K for 4 hr, quenched into ice water, and polished by immersion in a solution of 46.75 vol pct H2O, 46.75 vol pct concentrated HNO3, and 6.5 vol pct HF (49 pct) at 298°K. Special heat treatments given to a few specimens are described in the results below. Tensile tests were done on an Instron machine and were begun within 20 min after the quench, except where specified otherwise. Tests at 298°K were in air, at 178°K in acetone, and at 77°K in liquid nitrogen. All tests were at a strain rate of 8x sec-1. RESULTS AND DISCUSSION Yield Stress at 298°K. The compositions of alloys and the corresponding yield stresses (0.2 pct strain) are given in Table I. A plot of the yield stresses of the oxygen alloys, A, B, C, and D, indicates that varies linearly with CO1/2, where Co is the oxygen concentration, Fig. 1. This is in accord with Fleischer's6 theory for solution strengthening if the oxygen atoms do not cluster, or the cluster size remains constant with increasing oxygen concentration. In Fig. 1, it appears that if one could prepare some oxygen-free zirconium its yield stress would be very low. Therefore, we shall assume that for the oxygen alloys is equivalent to O0, the strength increment contributed by the presence of oxygen. The relationship between0.2and Co is expressed by 0.2 = 31.3 CO1/2, when the yield stress is in kg per sq mm and the concentration is in at. pct. Each of the hydrogen alloys, Al, A2, A3, and A4, contained 0.081 at. pct 0 as an impurity. In Fig. 1, it appears that this small amount of oxygen makes a significant contribution to the strength which cannot be ignored when we evaluate the contribution of the finely dispersed hydride. Let us assume the validity of the following equation: a0.2 = (a2o+a2R)1/2 [2] which is analogous to Eq. [I] for single crystals, and calculate values of UH for the hydrogen alloys by using the experimental values of 0.2 and o (0.081 at. pct) = 8.9 kg per sq mm. For 0.36 at. pct H, oH = 6.47; for 0.72 at. pct H, OH = 11.30; for 2.16 at. pct H, OH = 19.4; and for 3.60 at. pct H,

Jan 1, 1965

By W. F. Chiao

Ultrasonic attenuation, a, measurements in the frequency range of 5 to 55 mc per sec have been studied to determine their quantitative relationship with the following three variables of mixed microstructures in steels: 1) the volume percent, XF, of polygonal fer-rite in mixed structures of martensite and polygonal ferrite in Fe-Mo-B alloys: 2) volume percent, XA, of retained austenite plus martensite aggregates in high-carbon steel; and 3) substructural differences between 100 pct bainitic ferrite structures formed at various temperatures. The quantitative relationship obtained in the first two conditions by plotting a us the known structural parameters can be expressed, respectively, as: where al, a 2 and C1, Cz are constants. In the third condition the nature of the attenuation depends on the state of dislocations generated at the transformation temperatures and also on the alloy composition. From these measured results, the mechanism of ultrasonic attenuation caused by these mixed microstructures can also be studied. MUCH interest has recently been shown in the application of ultrasonic attenuation and wave velocity measurements to the study of the microstructural characteristics of steels. The general aims of most of the investigations in this field can be grouped into two categories: one is to study the mechanisms of ultrasonic losses caused by the characteristic phases in the microstructure of steel,''' and the other is to develop nondestructive test methods and applications for quality control.~' 4 Apparently no work has been done on the evaluation of ultrasonic attenuation meas -urements as a means of quantitative determination of a given phase in the microstructure of a steel. It is well-established that the decomposition of austenite results in four main microstructural constituents—polygonal ferrite, pearlite, bainite, and martensite—and that each phase has different mechanical properties. Thus, when a steel consists of mixed microstructures, the mechanical properties can often be related to a quantitative measure of the volume percent of each phase present. This study relates ultrasonic attenuation measurements to: 1) the volume percent of polygonal ferrite in mixtures of martensite and polygonal ferrite in Fe-Mo-B alloys; 2) the substructural differences between 100 pct bainitic ferrite structures formed at various temperatures; and 3) the vol- ume percent of austenite in austenite plus martensite aggregates in a high-carbon steel. The choice of the specimen materials was based on the laboratory stocks which were suitable to produce the required mixed microstructures for this study. EXPERIMENTAL PROCEDURES Materials and Heat Treatment. Polygonal Ferrite Plus Martensite Structures. This mixture of phases was produced in a vacuum-melted Fe-Mo-B alloy. The alloy was hammer-forged at 1900" ~ to a -f-in.-sq bar. By isothermally heat treating the alloy at 1300° F for various times and then water quenching, variations in the amount of polygonal (or proeutectoid) ferrite can be controlled in a microstructure in which the balance of the material is martensite. In the present work, four different times of isothermal transformation were adopted; after heat treatment, the four specimens were machined for ultrasonic measurements. The compositions, heat treatments, and dimensions of the four specimens are listed in Table I. 100 pct Bainite Structures Formed at Different Temperatures. It has been well-established by Irvine et al.= that the presence of molybdenum and boron in ferrous alloys can retard the formation of polygonal proeutectoid ferrite and expose the bainitic transformation bay, so that a more acicular or bainitic ferrite can be obtained over a wide range of cooling rates. Their investigation6 also showed that the mechanical properties of fully bainitic steels are usually closely dependent on the substructural characteristics of the steels. For studying the substructural characteristics in completely bainitic structures, six Fe-Ni-Mo alloys, of which five were free from carbon addition and one with 0.055 pct C addition, were selected so that a wide range of hardness values for 100 pct bainitic ferrite structures could be produced by normalizing at 1900" F followed by air cooling. The different bainitic transformation temperatures were recorded during air cooling. All of the alloys were vacuum-melted and then forged at 1900" F to square bars. Data on the six specimens of these structure series are summarized in Table 11. Austenite Plus Martensite Structures. The high-carbon steel used to study austenite plus martensite structures was vacuum-melted and then forged into Q-in.-sq bar. The series of mixed structures of austenite plus martensite was produced by quenching the specimens from the austenitizing temperature to room temperature and then refrigerating them at various temperatures within the range of martensite transformation to produce different amounts of retained austenite. Data on the four specimens of this series are listed in Table 111. Quantitative Analysis of the Microstructures. The microstructures containing martensite plus polygonal ferrite were analyzed by the point-counting technique.

Jan 1, 1969

By Thomas P. Forbath

THE sudden realization that known sulphur reserves amenable to mining by the Frasch hot water process are nearing exhaustion focused attention on widely scattered surface deposits throughout the world. These deposits are not necessarily of lower sulphur content than ores located underneath Louisiana or Texas salt domes which usually average about 30 pct sulphur disseminated in limestone matrix. Their near surface occurrence, however, renders exploitation by the Frasch process impossible. As is well known, the Frasch process depends on the presence of 500 to 1000 ft of overburden and cap rock above the sulphur deposits to permit melting underground sulphur in place by diffusing hot water under pressures of 200 to 600 psig in the formation and raising the molten sulphur to surface by air lift. This process renders possible the production of pure sulphur which is 99.5 pct pure without any subsequent treatment. Surface deposits contain sulphur in the same range of concentrations as the salt dome deposits, i.e., from 10 to 50 pct sulphur, associated with various gangue materials such as silica, limestone, and gypsum. The pirincipal distinction, then, does not lie in the percentage of sulphur contained in the ore, but in the geological nature of the deposit. A recent study' of the world sulphur supply situation estimated 1950 sulphur production in the free world countries at 5.6 million long tons, of which the United States produced 5.2 million tons, or 93 pct of the total. While America's domestic needs alone would have been covered by national production, about 1.4 million tons were exported during the same year. Despite all the steps which are being taken to restrict use of elemental sulphur and to force the fullest possible development of alternate sulphur sources here and abroad, the deficit in elemental sulphur production will probably increase with time. As a result of intensive prospecting for oil throughout the Gulf Coast area discovery of significant new salt domes is held unlikely. With the growing scarcity of sulphur and what appears to be an inevitable rise in price, recovery from deposits not amenable to Frasch-process mining assumes greater economic importance. Untapped Reserves The most important deposits in this category are located in Sicily, where elemental sulphur occurs in Miocene limestone and gypsum formation. Sulphur content of these ores ranges from 12 to 50 pct with an estimated average of 26 pct. Although quantitative estimate of these reserves is not available it is held that they exceed 50 million tons of sulphur. Similar deposits occur also on the mainland which contribute about one-third of Italy's total current annual production of 230,000 tons, but these are known to be nearing exhaustion. Significant surface deposits of volcanic origin are located in South America, Japan and western United States, silica being characteristic gangue con-stituent. The largest of these deposits are in South America. More than 100 extend over a zone 3000 miles long, paralleling the west coast of South America. 'Total sulphur content of these deposits has been estimated to be as high as 100 million tons. The main islands of Japan also possess at least 40 known volcanic sulphur deposits with probable reserves of 25 to 50 million tons.' Prospected reserves in western United States might amount to 2 million long tons, principal deposits being located in the northwestern part of Wyoming, southern Utah, and eastern California. Volcanic deposits of lesser importance are found around the Mediterranean, in Turkey and Greece, and in Africa, Egypt, Abyssinia, and Somaliland. Beneficiation Methods Different methods of beneficiation have been used in these various locations. In Italy the Calcarone kiln and Gill regenerative furnaces are used exclusively. Both utilize heat liberated by burning part of the sulphur in the ore to liquify or vaporize the remaining sulphur, which is recovered by solidification or condensation. The Calcarone kiln is of conical shape, generally 35 ft in diam at base and 18 ft high. A kiln of 25,000 cu ft capacity burns for about two months and yields about 200 tons of sulphur. The Gill furnace consists of a series of chambers with domed roofs. Sulphur is burned and melted in one chamber at a time and the hot combustion gases are used to preheat the ore charge in the subsequent cell. These furnaces operate on a cycle of 4 to 8 days. The recovery yield of both systems is about 65 pct. Sulphur losses amount to 25 pct through the combustion to sulphur dioxide; about 10 pct is retained in discarded calcines. Ores containing less than 20 pct are not considered suitable as furnace feed. These methods are not only wasteful because of the low recovery obtained, but represent a serious atmospheric pollution problem. Sulphur produced ranges from 96 to 99 pct purity and thus does not match Texas or Louisiana sulphur. Owing to the present shortage, sulphur in the Middle East sells

Jan 1, 1954

By H. B. Probst

Mechanical twins were produced in zone-refined tungsten single crystals by explosive working at room temperature. These twins are parallel to (112) planes and have irregular boundaries rather than the classical plane twin boundaries. These boundaries aye grooved surfaces in which the grooves themselves are parallel to a <111> direction and the sides of the grooves appear to be par-allel to (110) planes. TWINS were produced in tungsten single crystals by explosive working at room temperature. These twins differ in character from any previously reported for tungsten; however, they are similar to those found in molybdenum after compression at -196°C.1 Deformation twins "resembling Neumann bands in ingot iron" have been observed in tungsten by Bech-told and Shewmon.2 This observation was made with sintered polycrystalline tungsten pulled in tension to fracture at 100°C and using a strain rate of 2.8 x 10-4 sec-1. More recently Schadler3 found deformation twins in zone-refined tungsten single crystals pulled in tension at -196"&apos; and -253°C. These tests were conducted using a strain rate of 3.3 x l0-4 sec-1, and the twin bands were found to be parallel to a (112) plane. Deformation twins in tungsten&apos;s sister metal, molybdenum, were observed by Cahn.4 These twins were produced by compressing small (0.7 mm) vapor-deproducedposited molybdenum single crystals at -183°C. The compression was performed &apos;by impact." By the use of precession X-ray techniques, Cahn was able to identify the twin plane as {112} and the shear direction as <1ll>. Mueller and Parker1 produced deformation twins in polycrystalline electron-beam-melted molybdenum by compression at -196°C. Their "loading rate" was 5000 psi per min which, judging from their stress-strain curve, corresponds to a strain rate of approximately 0.3 x 10-4 sec-1. These twin bands were found to be parallel to (1 12) planes; however, they differed in appearance from previously observed twins. In place of straight and parallel twin boundaries they were found to be irregular, jagged, and sawtoothed. The sides of the saw teeth were identified as (110) planes and irrational planes of a (111) zone. The twins observed in the present work in tungsten single crystals are similar in appearance to those of Mueller and Parker in polycrystalline molybdenum. The starting material used in this investigation was 3/16-in. diam commercial tungsten rod produced by powder-metallurgy techniques. This material was converted to a single crystal by the electron-bombardment floating-zone technique.= The process was carried out in a vacuum of 10-5 mm of Hg using a traversing speed of 4 mm per min. Segments (=2 in. long and 3/16 in. in diam) of two crystals (A and B) produced in this manner were studied. Crystal A received one zoning pass, while crystal B received two passes. The two crystals were explosively worked at Bat-telle Memorial Institute in the following manner. A 1/2-in.-thick layer of plastic was applied to the crystals to serve as a buffer in an attempt to prevent cracking. The composite, crystal and buffer, was then wrapped with 1/8-in.-thick DuPont sheet explosive EL506A2 and detonated in water at room temperature. Metallographic samples of the worked crystals were prepared, and back-reflection Laue X-ray patterns were obtained using unfiltered molybdenum radiation. RESULTS AND DISCUSSION Blasting the crystals as described above failed to prevent cracking. The crystals fractured into several fragments about 3/16 to 1/2 in. long; however, the fragments were of sufficient size to be useful for the subsequent study. The diamond pyramid hardness of the crystals after blasting was in the range 430 to 450 as compared with 340 for the as-melted material, which shows a definite hardening resulting from plastic deformation. These hardness values were obtained using a 1000-g load and taking readings only in sound portions of the crystals free of cracks. The crystals exhibited profuse twinning as shown in Fig. 1. No such structure is present in the as-melted condition. Most of these twins have jagged twin boundaries and are similar in appearance to those found in molybdenum by Mueller and Parker. The twins in both crystals were found to be parallel to {112} planes. This identification was made by using the conventional two-trace method. Subsequent efforts to describe these twins more fully were carried out on crystal A. If the longitudinal axis of crystal A is placed in the (001)-(011)-(Il l) basic triangle of the standard cubic stereographic projection, as in Fig. 2, then the two sets of twins shown in Fig. 1 are parallel to the (112) and (121) planes. Fig. 3 shows a schematic representation of a twin with jagged boundaries. This type of twin with a <111>

Jan 1, 1962