Search Documents

By B. Ancker, E. R. Parker

The small-angle dislocation-boundary density of nickel and some of its alloys was investigated as a function of strength. It was found that the strength is a linear function of the density for pure nickel and for Ni-Ti alloys, but not for Ni-Co alloys. The results are interpreted in terms of dislocation theory. THE processes by which metals deform are still poorly understood, even though there has been a concentrated research effort in this field during the past fifteen years. Single crystals and individual grains in polycrystalline aggregates often develop a subgrain structure when deformed. Substructure formation occurs as a result of local plastic bending and its development is augmented by annealing at moderate temperatures.' Investigations made to date, largely qualitative in nature, have shown that the substructure is composed of a network of small-angle boundaries. The purpose of the present paper is to present some quantitative data on substructure (obtained by X-ray techniques) and to correlate the results with the strength of the metal in which the substructure existed. Nickel and some of its alloys were used for the experiments. The presence of certain kinds of foreign atoms greatly alters the formation and effect of the substructure. Polycrystalline tensile specimens were prepared from high purity nickel, two Ni-Ti alloys, and two Ni-Co alloys; the composition of these materials is given in Table I. Tensile specimens having a reduced section 1/4 in. diam by 3 in. long were machined from materials which were refined, cast, forged, and swaged in the authors' laboratories. The specimens were annealed for 1/2 hr at 1150°C, a treatment which produced grains having a diameter of approximately 0.5 mm. This rather large grain size was used to facilitate later examination of the substructure. Duplicate specimens of pure nickel and of each alloy were strained various controlled amounts in tension at room temperature and recovered for 1 hr at 800°C. (This strain induced prior to recovery will be referred to hereafter as the "prestrain.") One of each of the recovered specimens was pulled in tension at room temperature to obtain the stress-strain curves. The other bar was cross-sectioned in the gage length, carefully polished, and deeply etched to remove all cold-worked material. This specimen was then examined by the X-ray microscopic technique2 to determine the nature and magnitude of the substructure. The use of duplicate specimens permitted a direct correlation between tensile properties and substructure characteristics. X-ray microscopic studies were made with an apparatus constructed to have an average resolving power of about 4 microns. The monochromatic radiation reflected from grains which were oriented favorably was recorded on fine-grain photographic plates. These were examined at Xl00 which represents the true magnification, since there was essentially none in the X-ray system. This permitted detection of angular differences between subgrains of as little as 10 min. A very large number (over 100) of grains for each specimen was examined at this magnification and the number of substructure boundaries per grain was tabulated. An accurate relative measure was obtained under fixed experimental conditions, e.g., same grain size, resolving power, magnification, etc. The greatest source of error lay in counting the subboundaries, but this proved to be surprisingly small when repeated counts were made by different investigators. Furthermore, as the subsequent analysis indicates, whatever human error was introduced was not significant. Analysis of the Substructure Data When the number of substructure boundaries per grain was plotted against the frequency, skew (asymmetrical) distribution curves of the type shown in Fig. 1 were obtained. These are smoothed curves for the 1 atomic pct Ti in nickel alloy subjected to various prestrains and are representative of the curves obtained with all the metals. Since skew-distribution curves cannot be analyzed by elementary statistical methods, detailed investigations were undertaken in order to 1—test the consistency

Jan 1, 1955

By F. A. Crossley

The additions: o, Sn, Zr, Cb, Mo, and V were studied for their effects on the a/a + Ti3A1 boundary of Ti-A1 base ternary alloys. These additions were chosen because of their importance to commercial titanium alloy practice—oxygen as the major impurity (in terms of influence on mechanical properties) and the others as major alloying additions. The composition and temperature ranges of investigation varied according to the system as shown in Table I. The principal means of investigation was light-microscopy. The basis of the ternary studies was the recently determined titanium-rich end of the Ti-A1 system.' MATERIALS, METHODS, AND PROCEDURES Ingots were prepared in 15 g size for all alloys except those of the Ti-Al-Mo and Ti-A1-V systems which were prepared in 20 g size. These alloys were prepared from high-purity, electrolytically refined titanium obtained from the Bureau of Mines, Boulder City Metallurgical Laboratory, Boulder City, Col., and Chicago Development Corp., Riverdale, Md., for the 15 g and 20 g ingots, respectively. The alloying additions were obtained from commercial sources, and were at least of 99.9 pct purity, except the V-15 wt pct A1 master alloy for vanadium additions. This alloy was of 98.2 pct purity; the major impurities were 0.37 wt pct Fe and 0.38 wt pct Si. The ingots were prepared by multiple melting in a nonconsum-able-electrode arc furnace. The ingots were covered with Markal coating "CR-T22" (Markal Co., Chicago) to protect them from scaling. Markal coating is applied as a slurry of solids in a volatile carrier liquid. The coated article is air dried. Upon heating to the hot working temperature the coating fuses to form a glassy, protective film. The ingots were heated to 1000°C and rolled to 0.16 in. thickness, reheating after each pass. The glassy film was removed by sand blasting. Ten mils were removed from all surfaces by grinding and/or pickling (3 pct HF, 30 pct HN03, balance water). Grinding was always followed by pickling to remove traces of the grindmg medium. Hydrogen removal was accomplished by vacuum annealing at 1000°C for 2 hr to a final pressure of about 0.02 µ Hg. This was followed by annealing in evacuated Vycor bulbs (800°C, 2 hr, ice brine quenched (IBQ)) in order to maximize ductility. The materials were then lightly pickled and cold rolled to 0.12 in. thickness or cracking, whichever occurred first. After rolling, 1/2 by 1/4 in. specimens were cut for the equilibration anneals given in Table 11 and carried out as previously described.1 Final polishing of the specimens for metallographic examination was done by electrolytic means. Specimens were first examined without etching under bright field and polarized light illumination, for this is the only means of detecting the "two-phase" phenomenon.1,2 The term "two-phase syndrome" refers to Ti-A1 al-loys which appear on the bases of their chemical, phys-ical, and mechanical properties, to be two-phase, but are actually single phase and inhomogeneous in composition. Since processing and heat treatment of the alloys of the first group were completed before the two-phase syndrome was discovered, a number of the specimens were found to exhibit this phenomenon—80 pct of those of the Ti-Al-Zr system, 5 pct of those of the Ti-Al-Sn system, 30 pct of those of the Ti-Al-Zr system. The affected specimens were homogenized by heating above the ß transus temperature and quenching. The heat treatment and working schedules to accomplish this are given in Table 111. After rolling, specimens were pickled (3 pctHF, 30 pct HNO3 in water) to remove 1 to 2 mils from the

Jan 1, 1970

By Wolfgang O. Gluschke

INTRODUCTION Virtually all countries have general legislation covering most aspects of mining and mineral processing, including investment and tax laws, safety and health regulation, and specific rules applying to exploration reporting procedures. Yet most developing countries negotiate individual agreements with foreign investors for major mineral projects, in particular when all or part of the expected output is exported. Such agreements are known under several different names, as for example concession agreement; technical assistance agreement or service contract; management agreement; production sharing contract; joint venture; or contract of work. The objectives and conditions governing such contracts are as diverse as these titles are. While a mining project is only one among a large mumber of similar undertakings in industrial countries, in most developing countries such projects usually have a dominating position within the country's economy. The general legislation is often inadequate and insufficiently detailed to cover all aspects of large projects with extensive foreign inputs and international markets. Therefore, the need for specific arrangements. Specific arrangements may also become necessary due to the complexity of markets for mineral products, the large number of participants - equity partners, suppliers, financial institutions, buyers, etc. - or economic and political conditions in the host country and fiscal regulations in the home country of the investors. Most agreements are concluded between governments, or government agencies, in developing countries and private industry headquartered in developed market economies - Canada and the United States of America in North America, Western European countries, Japan, South Africa and Australia - or, in some instances, with state-owned firms in these countries. There are fewer agreements between developing countries and the centrally planned economy countries - Eastern Europe and the USSR - which cannot be considered here because of absence of adequate information. Agreements can be concluded at any of the different stages from exploration through closing of a mine. While most contracts are negotiated at the time when a potentially viable ore body has been delineated and investigated sufficiently to anticipate a realistic time for commencement of production, other agreements start with "grass-roots" exploration, during operation of a mine, smelter or refinery or, exceptionally, after a temporary shutdown. While several major provisions are included in all types of agreements, others are covered in much more detail (or not at all) in some contracts reflecting the specific work programme of each project. While it is realized that subdivision and classification of complex contractual relations is necessarily arbitrary to a certain extent, such grouping is useful in practical terms, especially for the non-specialist. The following discussion, therefore, will start with a description of the major forms of contractual relation - concession, investment contract, joint venture, management agreement, service and production sharing contracts. Major issues covered in all, or almost all, agreements will then be dealt with. These include stages of development, corporate structure and control, fiscal arrangements, infrastructure development and environmental law, and dispute settlement and renegotiation. FORMS OF CONTRACT RELATIONS There has been marked change in the type of agreement entered into during the 1960s (and before) and the early 1970s and the period thereafter. While straight "concessions" - granting of exclusive rights to foreign investors over a certain deposit or areas with little restriction or well-defined obligations - were the rule during the first half of this century, more complex agreements have evolved since. Some of these cover every conceivable aspect, foreseen or only possible, and define rights and obligations of all partners to the last detail. Over the years, agreements have without doubt become more beneficial for the developing host countries, in part as a result of well-publicized nationalizations, as for example in Chile, and in part due to changed attitudes in all countries

Jan 1, 1985

By G. H. Enzian

THE effects of phosphorus and nitrogen on the properties and behavior of low carbon steels are important considerations to both the manufacturer and the user of such material. For one thing, these two elements are the principal points of distinction between the chemical compositions of Bessemer and open-hearth grades. Further, open-hearth steels themselves may vary considerably in nitrogen content depending on the raw materials and ..practice used. In addition, appreciable amounts of rephosphorized open-hearth steels are being used at the present time in some low alloy, high strength grades and in certain other cases. For these and other reasons, it was felt desirable to obtain more complete knowledge of the effects of simultaneous small variations in phosphorus and nitrogen content on such properties of low carbon steel as tensile strength, impact toughness, sensitivity to cold work, and strain aging. A long range investigational program designed to furnish such information, was carried out, and the present paper summarizes some of the data currently obtained. General Discussion The strengthening effect of relatively small amounts of phosphorus in steel has been recognized, and to some extent utilized, for many years. A number of investigators have pointed out, from time to time, that under proper conditions of control, phosphorus can be a useful alloying element in low carbon steels.1-3 Nevertheless, in spite of the volume of data on the subject, the use of phosphorus for this purpose has always been handicapped by widespread feeling that its presence makes the steel brittle, difficult to cold work, and unreliable at low temperatures. This belief on the part of some steel users is not entirely justified in many cases, and often stems from early experience with acid Bessemer steels which have appreciably higher phosphorus (and nitrogen) contents than basic open-hearth grades. Nitrogen in steel, and its effect on properties, has also received considerable attention in the literature.4 In 1906, Braune5 ascribed the brittleness of Bessemer rails to the nitrogen content, and proposed a schedule of nitrogen limits for steels for various uses. This was one of the first papers to recommend directly the use of nitrogen content as a quality control measure, but it did not meet with much favor or support; as a matter of fact, the Jernkon-toret was rather sharply criticized for publishing the paper. In the intervening years a number of papers have been published on the subject of nitrogen in steel. As a result, it is now generally recognized that nitrogen can have an important effect on the properties and behavior of steel and that the presence of this element should be considered in any systematic study of factors affecting steel quality. In 1932, Graham6 introduced the concept that the quality of various steels could conveniently be evaluated and compared in terms of "sensitivity", which was defined as "the characteristic reaction rate of a piece of steel toward brittleness or hardness." This concept formed the basis for a broad research program into factors causing variations in the quality of commercial steels. At an early stage of this work, it became apparent that the observed effect of either phosphorus or nitrogen on the sensitivity of steel depended, in no small measure, on how much of the other of these two elements was also present. Indeed, the combined effect of phosphorus and nitrogen in steel may often be more important than the specific effect of either. With a few exceptions, it is not possible to predict accurately the effect on steel properties of a change in the amount of either phosphorus or nitrogen in the steel without also specifying the approximate amount of the other element. It seems desirable, at this point, to discuss briefly the distinction between strain-sensitivity and strain aging as used in this paper. In the early stages of its development, the work-brittleness test7 (which will be described in more detail later) was used

Jan 1, 1951

By D. C. Hilty, W. Crafts

A qualitative pseudoternary solidification diagram for the Fe-S-O system modified by manganese is proposed and supported by experimental derivation of an isothermal section at 1475°C and substantially 1 pct Mn in the solid-metal phase. The solubility of oxygen in liquid iron containing 1 pct Mn and sulphur contents up to 0.9 pct has been measured at 1600°C and found to increase as the sulphur content exceeds 0.15 pct. AS has been pointed out previously,' all oxide and sulphide inclusions in steel, excluding mechanically entrained matter, result from modifications of the basic Fe-S-O system. Therefore, one method of attack on the problem of the mechanism of deoxida-tion and inclusion formation is to examine the influence of deoxidizers on the Fe-S-O liquidus diagram. This approach has been adopted in work on the subject being carried out at the Metals Research Laboratories of the Electro Metallurgical Co. and has resulted in the postulation of types of solidification diagrams for the Fe-S-O system modified with common deoxidizers. This method has been of considerable assistance in interpretation of the factors affecting the formation of inclusions.' For the most part, however, these diagrams have been highly speculative. They were constructed primarily by inference from the pertinent binary and ternary systems and from the appearance of inclusions produced by specific deoxidizers. Moreover, since the systems contain four or more components, simplifying assumptions and qualifications were necessary in order to permit their representation as modifications of the Fe-S-O diagram. Consequently, experimental confirmation of the type of construction employed for these diagrams is desirable, and the present paper is concerned with such evidence. The addition of manganese to the Fe-S-O system was selected for initial investigation. Manganese is an important component of practically all steel and is, therefore, of fundamental interest. Because of its great affinity for sulphur, manganese has a more profound effect on the sulphide side of the Fe-S-O diagram than any of the other common deoxidizing eletments. In addition, more quantitative information pertaining to the contiguous binary and ternary systems involving iron, manganese, sulphur, and oxygen is available than for the other deoxidizers. Fe(Mn)-S-O Solidification Diagram Many of the side diagrams of interest in connection with this system are known. The Fe-Mn binary has been well established by several investigators;% the Fe-O system has been clarified by Darken and Gurry;' and FeO-MnO diagrams have been published by Herty,5 Andrew, Maddocks, and Howat,6 and by Hay, Howat, and White.7 The Mn-O diagram does not appear to have been derived experimentally, but the portion of 'it that is of interest can be visualized easily from a consideration of the melting points of Mn and MnO and the fact that a broad miscibility gap must exist between them. The ternary Fe-Mn-O diagram has been discussed by Benedicks and Löfquist,8 Wentrup,9 and by Körber and Oelsen.10 Experimental study of the Fe-Mn-O system has been carried out by Körber and Oelsen and to a lesser degree by Chipman, Gero, and Winkler,11 Hilty and Crafts," and others. On the basis of published information, then, an outline of the Fe-Mn-O liquidus diagram can be constructed as illustrated by Fig. 1. In this, as in all succeeding diagrams, transformations such as the ?-? reaction in iron have been ignored for purposes of simplification. The arrows denote falling temperature. The distinguishing feature of this diagram is the broad region of liquid immiscibility extending

Jan 1, 1955

By Walter Crafts, D. C. Hilty

The liquidus diagram for the iron field of the Fe-S-O system has been derived experimentally. The solubility of oxygen in molten Fe-S alloys has been measured at several temperatures and found first to decrease slightly and then increase rapidly with increasing sulphur content. IT is well-known that the types, distribution, and amount of inclusions in any given steel are closely allied with the oxygen and sulphur contents of the metal, so that the term "deoxidation" as conventionally applied in steelmaking refers not only to the control of oxygen and oxide inclusions but to collateral modification of sulphide inclusions as well. Inclusions commonly occur in steel as mixtures of oxides and sulphides, sometimes closely associated and sometimes apparently completely divorced. Their effects on the properties of the metal may vary widely depending on their types and distribution as well as on their quantity. The mechanism by which inclusions form, however, and the means for controlling it are exceedingly complex and not well understood. Nevertheless, as described by Benedicts and Lofquist' and by Wentrup' all oxide and sulphide inclusions, excluding mechanically entrained matter, must result from specific modifications of the basic equilibria of the Fe-S-O system. In an investigation of the mechanism of deoxidation and inclusion formation being carried out at the Union Carbide and Carbon Research Laboratories, it became evident that quantitative data for the fundamental Fe-S-O system were essential to a better understanding of the problem. Consequently, the liquidus surface of the portion of the Fe-S-O diagram of interest in steelmaking has been determined, and the solubility of oxygen in liquid iron containing sulphur has been measured. Fe-S-O Diagram The portion of the Fe-S-O diagram of specific interest in steelmaking is the Fe-FeS-FeO corner. Of the contiguous binary diagrams for this partial system, only those for Fe-FeO and Fe-FeS are known with reasonable certainty.3 Several attempts to derive the FeO-FeS diagram have been made, chief among them being the work of Giani as reported by Oberhoffer,4 and, more recently, that of Vogel and Fulling." Although both Giani and Vogel and Fulling published diagrams for the FeO-FeS system, their observations indicated that the oxide component approximated Fe3O4 rather than FeO, so that there is doubt whether the system is truly binary. Benedicts and Lofquist and also Wentrup published qualitative sketches of the ternary diagram based on the assumption that it was of the type described by Marsh for a ternary system in which two of the contiguous binary systems are simple eutectic systems and the third exhibits a solubility gap in the liquid. Apparently, however, the only previous effort to establish the diagram experimentally was that of Vogel and Fülling whose results are reproduced in Fig. 1. Although defining approximately the limits of the miscibility gap and the location of the ternary eutectic, point E in Fig. I, Vogel and Fülling's diagram must be considered only semiquantitative. It was derived from a limited number of observations mainly along two quasi-binary sections, and the directions of the tie-lines and the location of the critical point, K in Fig. 1, were estimated by computation. Experimental Procedure The experimental technique employed for the present investigation involved melting mixtures of ferrous sulphide and ferric oxide in contact with iron. Since it was considered that in this part of the system the vapor pressures of sulphur and SO, are so low that they would not affect the results significantly, the melting was carried out in a current of argon to sweep out products of any initial re-

Jan 1, 1953

By F. W. Jessen, N. Mungan

The plastic flow characteristics of clay water suspensions were first recognized by Binghaml in 1916 and further studied by Ambrose and Loomis' in 1931-1932. Many physical and chemical properties of clay suspensions are deeply influenced by particle size distribution as well as by the exchange cation. Separation of a number of fractions is particularly helpful in making a detailed study of a suspension. A limited amount of work has been conducted on the effects of particle size on specific properties of mud drilling fluids. A knowledge of particle size distribution is essential for a better understanding of fundamental properties of clay suspensions, particularly viscosity and gel strength. The purpose of this investigation was to study the effect of particle size on the viscosity and gel properties of clay suspensions and to determine any changes in mineral composition accompanying variation in size. EXPERIMENTAL PROCEDURE Numerous methods may be used to fractionate a poly-dispersed clay suspension into nearly mono-dispersed systems. Particles larger than 44 micron (325-mesh) may be obtained by wet screening. For sizes from 44 to 1 micron, the fractionation may be accomplished by sedimentation, elutriation plus microscopy,' or turbidimetry.6 When dealing with systems which contain a major portion of particles below l-mi-cron size, slit ultramicroscopy,7 centrifuge methods of Svedberg,5 eams,9 and McBain,l0 and X-ray methods" may be employed. The Sharples super centrifuge has been proved suitable for the economical separation of suspended particles of submicron range.12,13,14 Three-hundred and forty-five gal of 1 per cent by weight Aquagel suspension was mixed in a high-speed 5-gallon mud mixer and the suspension was allowed to stand in 5-gal glass bottles for 30 days to assure full hydration. All particles larger than 44 microns were removed by wet screening. Five fractions in the range of 44 to 1 micron were obtained by sedimentation. The Sharples super centrifuge method described by Hauser and Reed.13 Norton and Speil,14 Ha user and Lynn,15 and Hauser and Schachman,10 as employed to yield five fractions of submicron size. A super centrifuge operation chart prepared by Fancher, Oliphant, and Houssierre12 was used to determine operating conditions for the desired particle ranges. The clay suspension, after removal of all material above 1-micron size, was centrifuged. During rotation of the bowl, the particles were deposited on a thin plastic liner which fit closely inside the bowl. After passing 5 gal of suspension through the centrifuge, the liner was removed and the particles on the liner were scraped off, remixed and recentrifuged. This was repeated five times. Thus, Fractions H, I, K, and L were obtained. The finest fraction, M, remained in an approximately 0.4 per cent suspension. Water was evaporated under reduced pressure from this suspension to yield the finest particles. Results of the fractionation are included in Table 1. In order to achieve proper orientation and humidity conditions for X-ray diffraction analysis, the following procedure was employed. The dilute suspensions were placed on glass slides and allowed to dry at room temperature. This was repeated several times. Previous to measurement each sample first was placed in a water desiccator for half an hour and then allowed to stand for 15 minutes at room temperature. thus yielding an equilibrium condition of about 52 per cent relative humidity. The relative humidity of the room in which the X-ray diffraction was carried out was maintained at approximately 50 per cent so that there was no change in the spacings during the X-ray examination. Samples

By A. J. Griest, A. P. Young, P. D. Frost

A study was made of the effect of ß grain size on the tensile ductility of a-ß titanium alloys haet treated to strengths in the range 165,000 to 180,000 psi. It was concluded that the primary cause of $ embrittlement is the large grain size which is obtained so rapidly in many a-ß titanium alloys when heated into the ß field. It was also found that grain-boundary a and acicular intragranular a in the microstructure are not necessarily detrimental to the tensile ductility if the prior ß grain size is sufficientlv fine. BETA embrittlement occurs in a-ß alloys when the finish hot-working or stages of the subsequent heat treatment have been done at temperatures in the ß field. In order to avoid this embrittlement, for-gings must be finished at temperatures under the ß transus. For better workability, it would be desirable to forge at temperatures in the ß field. For this reason, research has been conducted to determine the cause of ß embrittlement. Although ß grain size has been considered as a possible contributing factor to ß embrittlement, only limited research has been done to evaluate its effect. The effect of addition elements on ß grain growth in unalloyed titanium has been studied by Liu.' Complete mechanical property data, however, were not obtained for correlation with grain size. Ogden and coworkers2 have investigated the effect of both a and ß grain size and grain shape in Ti-Cr-Mo alloys. Their data indicated that tensile ductility was relatively insensitive to ß grain size at annealed strength levels (140,000 psi ultimate). No study has apparently been made to determine the relationship between ß grain size and ß embrittlement in a-ß alloys heat treated to high strength (180,000 psi ultimate). The high-strength condition is that in which ß embrittlement is particularly severe. Beta grain size is not easily varied through the use of a range of annealing times and temperatures in the ß field, because of the rapidity with which ß grain growth occurs. Short-time heating techniques are needed to obtain the smallest ß grains. EXPERIMENTAL PROCEDURE In this research, the experimental procedure was as follows: 1) Tensile and microimpact specimen blanks were resistance heated to produce a range of ß grain sizes. 2) The tensile and microimpact blanks were then solution-treated in the a-ß field and aged to produce predetermined ultimate tensile strengths in the range 165,000 to 180,000 psi. Tensile and impact ductility were then correlated with a grain size for constant strength levels. Materials—The alloys Ti-6A1-4V, Ti-4A1-4Mn, and Ti-16V-2.5Al are representative of the heat-treatable a-ß types in which ß embrittlement can be severe. The compositions and ß transi of the heats studied are given in Table I. The Ti-16V-2.5Al heat was melted and fabricated to 1/4 in. round at Battelle; the other alloys were obtained as 1/4 in. round rolled from commercial heats. The micro-structures of the as-rolled alloys were fine a-ß mixtures. This structure is typical of titanium alloys hot-worked extensively in the a-ß field. Resistance-Heating Treatments to Establish ß Grain Size—For direit resistance heating, a 1/4-in.-diam s~ecimen was heated with a 3000-watt transformer rated for 1000 amp with adjustable secondary voltage. The apparatus included a spray quench fixture. Both heating and quenching periods could be automatically controlled. For temperature measurements, 28-gage Chromel-Alumel wires were tack-welded to the specimen at its midpoint. A high-speed recorder-controller was used to control the heating cycle.

Jan 1, 1960

By H. Majima, E. Peters

The oxidation rates of pyrite, pyrrhotite, chalcopy-rite, chalcocite, covellite, bortzite, galem, sphalerite, and stibnite have beet2 carefully compared at 120 oC, using aqueous phosphate solutions buffered at pH 2.7, 7.1, and 11.2 as zvell as in 1.0 M NaOH and in am~onia-containing solutions. Under conditions of these studies, variations in rates of more than a factor of two per unit of mineral surface area were not observed for different minerals in phosphate buffered solutions, except that pyrrhotite showed unusually rapid oxidation rates In acid solutions, and stibnite showed very rapid oxidation rates in basic solutions. In caustic solutions, oxidation vates varied considerably, the most rapid being for pyrite and the slowest for sphalerite. In the presence of ammonia, only copper minerals showed more vapid oxidation rates in neutral solutions. Oxidation rates were measured on the basis of oxygen-consunzption rates per unit surface area of exposed mineral, using equivalent surface areas, mineral mesh sizes, pulp densities, agitation, and oxygen pressure. OXIDIZING pressure leaching of sulfide minerals has been a commercial method of treating concentrates of nickel, cobalt, and copper for over a decade, and recently the treatment of lead and zinc concentrates,' and even the recovery of sulfur from pyrrhotite by such processes, have been shown to be technically feasible.' Furthermore, the leaching of uranium ores by solutions acidified through the oxidation of pyrite and other sulfides present in the ore is commercially practicable.3 5 The mechanisms by which sulfide minerals oxidize in aqueous environments are not yet understood. However, certain generalities can be made for systems containing no unusual reagents, in which oxygen is the oxidizing agent. In acid solutions the oxidation of many sulfide minerals leads to the formation of elemental sulfur.8'7 The attack of some sulfides, notably galena (PbS), sphalerite (ZnS), and pyrrhotite (FeS), at temperatures below about 110°C leads to an almost quantitative conversion of the sulfide mineral to elemental sulfur according to the equation: On the other hand pyrite (FeS2) and copper minerals do not yield elemental sulfur quantitatively. Not more than 50 pct of the sulfur in pyrite is converted to the elemental form in acid oxidation, and this percentage is lowered by higher temperatures and possibly higher oxygen pressures.8'9 Copper minerals sometimes yield elemental sulfurlo'" but, since sulfur can reduce cupric ions with the precipitation of CuS or Cu2S, it is not likely that elemental sulfur is obtained in good yield. In acid solutions sulfur that is not converted to the elemental form always appears as sulfate salts or sulfuric acid. Sulfide minerals that yield elemental sulfur nearly quantitatively on oxidation in acid solutions tend to generate H2S when the conditions are not oxidizing. However, although it has been suggested that H2S is an essential intermediate in the mechanism of oxidation," this is doubtful because it is oxidized very slowly by oxygen. In basic solutions elemental sulfur is never formed, as predicted by thermodynamical considerations, summarized in the form of a potential-pH diagram in Fig. 1. However in neutral and basic solutions, unstable forms of sulfur not shown in the diagram appear during oxidation, notably thiosulfate (SOT-) and polythi-onates (S,O;-, x = 2 to 6).l3,l4 On prolonged oxidation only hexavalent sulfur as sulfate (SOi-) or sulfamate (NHzSO;) is found to remain, the latter appearing only in the presence of ammonia. The oxidation kinetics for thiosulfate by oxygen are the subject of a separate paper,'5 and were studied on the presumption that thiosulfate may, in fact, be a primary oxidation product for sulfide minerals, at least in neutral and basic solutions. A survey of the literature indicated that comparative oxidation rates for different sulfide minerals were needed to establish the role of the mineral structure and its thermodynamic properties on the oxidation rates and mechanisms.

Jan 1, 1967

By M. R. J. Wyllie, Walter D. Rose

The use of electrical well logs for the quantitative determination of such reservoir parameters as connate water saturation, formation permeability and connate water salinity has recently been attracting the attention of a number of workers1,2,3,4,5. While the theory of the determination of connate water salinity from the self potential S.P. log has received rather detailed treatment, relatively little attention has been paid to the theoretical aspects involved in the quantitative interpretation of resistivity data. It is clear that if electrical well logs can be used for the quantitative evaluation of physical characteristics of reservoir rock, they will provide a valuable tool to supplement cheaply information obtained by more laborious core analysis. In certain cases it is conceivable that the coring program could be considerably curtailed if the electric log could be relied upon to give reasonably accurate quantitative information. It is our object in this paper to examine the theoretical basis of quantitative log interpretation as expressed in such well-established logging concepts as formation factor and cementation factor. It is also our object to investigate the physical aspects of the rela- tionship which is presumed to exist between resistivity index and brine saturation in reservoir rock. In particular, we will endeavor to draw attention to the fact that it is possible to express these logging concepts in terms of capillary pressure-saturation relationships, permeability and tortuosity, parameters which we will consider in this paper to be fundamentally indicative of rock texture. The probability of being able to obtain from log recordings alone the data theoretically essential to permit quantitative log interpretatiou will be examined, and consideration will also be given to the problem of formulating simple semi-empirical relationships for use in the field. THEORETICAL CONSIDERATIONS The concept of formation resistivity factor, or as it is now commonly called, formation factor, appears to have been introduced by G. E. Archie",". Formation resistivity factor as defined by Archie is the resistivity of a rock 100 per cent brine-saturated divided by the resistivity of the brine. This relationship had previously been used by physical chemists and the concept is, for example, implicit in an early treatment by Fricke7 of the conductivity of aqueous slurries. Archie was able to establish the fact that although the relationship between formation factor and porosity did not change too greatly for different rock textures, the relationship between formation factor and permeability was markedly dependent on rock texture. The relationship found by Archie between formation factor and porosity was expanded by H. Guyod8. Archie had suggested that the formation factor, F, could be expressed in terms of the porosity, F, as follows: F =F-m .... (11 where, F = effective porosity fraction of the sand, and, m = an exponent greater than about 1.3. Guyod introduced the term cementation factor for the exponent, m, since it may be qualitatively shown that formation factor for any given porosity tends to increase as the sand becomes more cemented. It would appear that the concept of m as an exponent representing the cementation in a sand is strictly only of qualitative value, since the quantitative estimation of m by means other than the determination of F and F has not been possible. Guyod suggests that m can be estimated from a visual examination of a sample of rock, but this method appears to be both arbitrary and unreliable. It can be said, however, that in intergranular reservoirs of practical interest, the range of m is from about 1.4 to 3.0. A possible explanation for the fact that m is found to lie between these limits will be offered below. The determination of formation factor is an essential preliminary to the

Jan 1, 1950

By R. J. Wagner, Jim Douglas Jr., P. M. Blair

The calculation of the behavior of an oil reservoir during a water flood has long been an important problem to reservoir engineers. Buckley and Leverett derived the differential equation which describes the displacement of oil from a linear porous medium by an immiscible fluid, but this equation could not be solved by the methods of classic mathematics. Consequently, in order to integrate the equation over the length of the reservoir, they neglected the effects of capillary pressure. In the present paper, a numerical method has been developed for determining the behavior of a linear wafer flood with the inclusion of capillary pressure. The differential equation which was derived for the case of incompressible fluids is second order and non-linear. This differential equation was approximated by an implicit form of difference equation which is second order correct in both time and distance. An electronic digital computer was used to perform the numerical solution of the difference equation. INTRODUCTION The problem of calculating the flow and distribution of fluids in an oil reservoir subjected to a water flood has long challenged the reservoir engineer. The ability to solve this problem would provide a valuable tool for the design and study of field waterflooding programs. One of the first contributions in this field was made by Buckley and Leverett,'.' who developed a method of calculating waterflood performance in a linear reservoir. Their technique was limited by the practical necessity of excluding quantitative consideration of capillary pressure. It is the purpose of this paper to describe a method for calculating the behavior of a linear water flood with capillary pressure considered. This method, although limited to the linear case, should serve as a step toward the solution of the two- or three-dimensional waterflooding problem which would better describe actual reservoirs. The treatment of the problem has assumed that the reservoir is linear and homogeneous, that both oil and water -are incompressible, and that gravitational forces may be neglected, and that water is injected into one end of the reservoir and oil and water are produced from the other. Though these assumptions may not be directly applied to natural reservoirs, they are closely approached in many laboratory core tests. The physical problem has been represented mathematically by means of a differential equation and suitable boundary conditions. The inclusion of capillary pressure effects causes the differential equation to be second order and non-linear, and it is not amenable to solution in terms of known functions. Nevertheless, the problem may be solved by numerical methods. The authors have successfully calculated several cases from the start of water injection through to depletion of the hypothetical reservoir using realistic capillary pressure and relative permeability characteristics that were chosen so as to put the method to a test. DIFFERENTIAL EQUATION In two- or three-dimensional problems involving the two-phase flow of incompressible fluids through a porous medium, the equations describing the flow can be combined until only two dependent variables remain, the pressure in one phase, and the saturation of one phase. These variables appear in two simultaneous partial differential equations which are derived on the basis of the conservation of mass and Darcy's law as applied to two-phase flow. For a linear reservoir, this system of differential equations can be reduced to a single, nonlinear, second order, parabolic differential equation for one dependent variable, the saturation of one phase. For incompressible water and oil, it can be shown that the conservation of mass for the two phases requires that Where q, and q, are the water end oil flow rates per unit cross-section, s, and so their saturations, +, the porosity, and [ and T, dimensional space and time variables. Darcy's law, as extended to two-phase flow including the effect of capillary pressure,' is as follows:

By P. H. Turnock, R. D. Pehlke

The effects of mutual interactions between elements on the solubility of nitrogen in liquid iron alloys containing chrormium, columbium, molybdenurn, nickel, or silicon have been determined. The equilibrium nitrogen solubilities in the liquid iron alloys were measured by the Sieverts method. The second-order effects caused by the presence of chrormium, columbium , or silicon in the melt were found to be significant. The solubility of nitrogen in liquid iron alloys containing several alloying elements has also been measured as a function of melt temperature and nitrogen pressure. The heat of solution of nitrogen in microcomponent iron alloys has been found to be a function of the logarithnz of the activity coefficient of nitrogen, irrvespectilje of the composition of the melt. The heat of solution of nitrogen in pure liquid iron was determined to be 1265 ± 450 cal per g-atom of N over the temperatxue range 1600° to 1800°. Sieverts ' law was obeyed for all melt compositions studied in the pressure interval 0.4 to 1.0 atm. The solubility of nitrogen in liquid pure iron and in many dilute binary liquid-iron alloys has been the subject of many investigations. Pehlke and Elliott' made a comprehensive investigation on this subject in which they summarized and compared previous researches. However, the practical application of these solubility studies requires prediction of the nitrogen solubility in the multicomponent alloys encountered in operating practice. The interaction parameter as defined by wagner2 has been used to predict the behavior of a solute in a complex alloy, and a modification of this approach has met with reasonable success in several instances.'1"5 The necessity for extending the nitrogen solubility data on binary alloys is caused by the lack of direct measurements of the solubility of nitrogen in molten alloy steels. A direct comparison between measured solubilities and those predicted for multicomponent alloys based on data in binary systems has been made in only a few instances, and most of these were reviewed by Langenberg,6 who developed a graphical technique for predicting nitrogen solubilities in multicomponent iron alloys. The solubility of nitrogen has also been studied over the entire Fe-Cr-Ni ternary system by Humbert and Elliott.? Both studies compared the solubilities computed from binary iron alloy data with those measured in ternary alloys in the liquid state at 1600°C. Most of the nitrogen solubility measurements referred to above were for iron alloys at 1600°C; thus, the ability to predict nitrogen solubilities in molten iron alloys has been restricted primarily to alloys at 1600°C. Recently, Nelsson8 devised an approximate method for calculating the change in nitrogen solubility in molten iron alloys as a function of temperature, but an assumption of regular solution behavior limits the applicability of the method. Chipman and Corrigan9,10 have developed an empirical correlation. based on the activity coefficient of nitrogen in the iron alloy at 1600°C, which facilitates the estimation of nitrogen solubilities at higher temperatures. This method also has limitations, due to the necessity of calculating the activity coefficient for nitrogen in a multicomponent iron alloy at 1600°C from the nitrogen solubility data for binary iron alloys. The shortcoming of such a straightforward approach to the problem of predicting nitrogen solubilities is the disregard of mutual interactions between elements in a complex iron melt. This investigation was conducted in an attempt to provide a more reliable method of predicting nitrogen solubilities in multicomponent iron alloys. The effects of the mutual interactions between chromium, nickel, silicon, molybdenum, or columbium in iron were studied, using the Sieverts' technique1,7 to measure the effects of alloy composition, temperature, and pressure on the solubility of nitrogen in ternary liquid iron alloys.

Jan 1, 1967

By L. A. Rapoport, W. J. Leas, C. W. Carpenter

A program of scaled flow model experiments has been undertaken to study the performance of five-spot water floods. The modeling procedures are discussed and the construction and operation of the flow models are described. The tests were specifically designed to investigate the effects of injection rate and of oil-to-water viscosity ratio on oil recovery. The areal sweep of five-spot water floods was also studied, and direct comparisons between five-spot and linear flooding behavior have been established to permir better utilizatfon of core analysis results for purposes of field waterflood evaluations. INTRODUCTION For a long time reservoir engineering procedures have been confined to the consideration of idealized, one-dimensional or linear systems.11* However, direct application of the linear waterflood mechanism to the evaluation of actual field performance is justifiable only for the purposes of rough approximation. To obtain a realistic analysis of field behavior, it is necessary to account for reservoir geometry and well distribution, and to take into consideration the configuration of the flow of fluids between injection and producing wells. In the past it has been attempted to account for the complex nature of the fluid movement in the reservoir by associating each well pattern with a certain "sweep-efficiency constant", e.g., a sweep efficiency of 72 per cent for the conventional five-spot pattern.' More recently, methods of waterflood evaluation have been improved by the introduction of variable sweep-efficiency factors, defining the areal coverage as a function of cumulative water injection and viscosity or "mobility" ' ratio.""' These variable sweep-efficiency factors, derived from miscible fluid displacement experiments, require the assumption of piston-like displacement of oil by water. This implies that only one fluid, either oil or water, is moving at any one point of the reservoir ht any moment. Such an assumption, however, is incompatible with the bulk of the experimental evidence and also theoretical considerations which indicate that usually in the area invaded by water, both oil and water are flowing in varying proportions, in accordance with the relative permeability characteristics of the porous medium. Another recent approach toward the evaluation of five-spot flooding behavior involves correlations of the areal sweep efficiency (observed in immiscible fluid experiments) with a "mobility ratio," based on the consideration of the "average" relative permeability to the displacing fluid. Combination of these correlations with a modification of the Buckley-Leverett frontal-drive equation leads to a computational procedure for predicting the recovery performance of areal water floods. While of considerable interest, this approach has no complete theoretical justification and, therefore, does not entirely clarify the features of five-spot water floods. Rigorous, theoretical treatment of waterflood behavior in two- or three-dimensional systems as encountered in the reservoir, including proper consideration of relative permeability and capillary-pressure functions, is as yet beyond the scope of analytical or computational methods. The, behavior of such water floods can, however, be studied by means of scaled flow model ekperiments, which is the approach followed in these investigations. The flow model experiments described in this paper simulate a uniform, horizontal formation of constant thickness with fully penetrating wells. The system under consideration is entirely liquid saturated and not influenced by compressibility effects, which is representative of the situation achieved after gas fill-up under practical conditions. Gravity segregation effects are assumed to be negligible, which implies a relatively thin reservoir formation and/or low vertical permeability. THEORETICAL CONSIDERATIONS For a flow system of uniform thickness with no appreciable effects of gravity segregation, the analysis of water-oil displacement can be confined to a two-dimensional treatment based on the consideration of a "slice" of unit thickness. Referring to previously established derivations," the displacement of oil by water in such a system can be described by the following, dimen-sionless equations:

By J. C. Savins

The practical analysis of the hydrodynamics of the wellbore has long been a subject of interest to engineers. This paper presents a simplified solution to the problem of computing the pressure drop for the flow of drilling mud in the annulus of the wellbore. This solution is, however, an exact and rigorous solution under the assumptions which have been imposed. These assumptions are that the drilling fluid is a Bingham plastic fluid* and that the annulus is formed by two concentric, stationary, cylindrical pipes. It is further assumed that the fluid is incompressible and that its motion is isothermal and in a steady state. This problem under the same assumptions has been attacked by previous authors. Beck, Nuss and Dunn' proposed that the equation for the flow of a Bingham plastic fluid in a cylindrical pipe could be applied to an annulus if the pipe radius were replaced in the equation by the hydraulic radius. This equation, known as the Buckingham-Reiner equation' (see Appendix 1), was also used in an approximate form. Van Olphen pointed out that even for a simple or Newtonian- fluid the pipe equation (Poiseuille's law) could not be converted to the Lamb equation escriptive of flow in an annulus (see Appendix 1) by using the hydraulic radius. Van Olphen further attempted to give a solution for the annular flow of a Bingham plastic fluid by introducing approximations similar to those which have been used in the case of the Buckingham-Reiner equation. Other attempts to provide approximate or exact solutions have been made by Grodde' and by Mori and Ototakeq. The present authors some years ago in unpublished work derived the correct expressions relating the pressure drop and flow rate for this problem. It was found that the solution consisted of two simultaneous equations, one of which contained a logarithmic term. Thus, obtaining numerical results for any particular case of interest involves very tedious trial-and-error computations. Very recently Laird presented the correct derivation of the two equations which are given in full detail in Appendix 1. In order to reduce the amount of calculation time which would be involved in providing a complete tabular or graphical solution to the problem, a high-speed electronic digital computer has been utilized. For this purpose the two simultaneous equations were transformed into more compact expressions by introducing reduced variables. These expressions are given in the following theoretical section. A similar procedure in this problem has been developed by Fredrickson and Bird1" Their tabular results, however, are very incomplete in the range of practical interest for problems of wellbore hydrodynamics. We have furthermore been able to express our graphical results in terms of convenient and familiar dimensionless groups. THEORETICAL DEVELOPMENT Use of Reduced Variables In terms of reduced variables the two simultaneous equations just discussed take the following form, The reduced variables q, x, a and z are defined in terms of the various measured quantities, where Q is volumetric flow rate, AP/L is pressure gradient, Dl is OD of inner pipe, D, is ID of outer pipe, is plastic viscosity, and is yield point. Thus, we have a dimensionless volume flux, a dimensionless reciprocal pressure gradient, and the ratio of the pipe diameters Before introducing the fourth reduced variable, z, it is of interest to consider the physical significance of the parameter x. As may 'be seen from the velocity profile of Fig. 1 the Bingham plastic fluid has the interesting property that a portion of the stream flows at a uniform velocity without shearing action. This section of the stream is situated approximately in the center of the conduit and is known as the "plug flow" region. Its existence is due to the fact that the shearing stresses within the region do not exceed the yield point, which is one of the two flow properties characterizing the fluid. The parameter x then turns out to be the ratio of the

By Robert Maes, Robert de Strycker

The Fe-Ni-As phase diagram has been established by the study of about a hundred alloys, by microscopic observation, and by thermal analysis, with arsenic contents up to 50 pct. The iron and nickel arsenides present extensive solid-solution fields, owing to the substitution of nickel by iron and vice versa; the extent of the solubility field of each compound has been determined with an accuracy of ±1 pct. In the investigated range of compositions, the solidification reactions were established, and the temperatures of the invariant reaclions detevmined with a preciston of iZ°C in the most favorable and ±5°C in the least favorable cases. The isothermal lines of the liquidus surface have also been drazum, with an accuracy estimated at i5°C. Reactions in the solid state, which take place for the formation or the decomposition of certain phases , were investigated in detail. ThE constitution diagram of the Fe-Ni-As system is fundamental for the understanding of the properties of nickel speiss; these by-products of the extraction of certain nonferrous metals indeed often contain the three elements iron, nickel, and arsenic as main constituents. The Fe-Ni-As system has already been the object of earlier investigations. In 1932, Guertler and Savels-berg' have presented some elements of the ternary phase diagram, for arsenic contents up to 55 pct (all percentages in this paper are given in weight percent), including the vertical section FezAs-Ni&s2. This investigation is, however, incomplete and certain anomalies suggested the necessity to verify the results published by these authors: the section Fe2As-Ni5AsZ, for example, is presented as quasi-binary, with a field of complete miscibility in the solid state, even though these compounds do not have the same structure. Recently, ~useck' established an isothermal section at 800°C in the Fe-Ni-As system by X-ray diffraction and microscopic examination of water-quenched alloys. These techniques are sometimes inaccurate for the determination of the fields where alloys are liquid at the investigated temperature, and they may lead to erroneous conclusions when a compound exists at the investigated temperature but is not stable at room temperature and cannot be maintained by quenching. LIMITING BINARY PHASE DIAGRAMS The Fe-As constitution diagram has been the object of several investigations which have been reviewed by Hansen and Anderko.3 For the Ni-As phase diagram, a recent study has been effected by Yund.4 In this system, Heyding and calvert5 have determined the existence of a compound of unidentified structure at arsenic contents slightly higher than those corresponding to Ni5As2 and at temperatures lower than about 200°C; by analogy with iron and cobalt arsenides, this compound could correspond to the formula Ni2As, as suggested by Kulle-rud; although this is not definitely established. In the region of arsenic contents from 35 to 55 pct, Fried-rich7 detected anomalies in the solidification reactions; an interpretation of these anomalies was given by Hansen ,8 which assumed that the solidification reactions in practice are not in equilibrium, but are meta-stable. The main features of the Fe-Ni constitution diagram are the existence of a complete miscibility field, at least at high temperatures, for the fcc phase (which will be designated by My in the remainder of this work) and of a limited solubility field for the bcc phase (designated by Ma). EXPERIMENTAL METHODS The fundamental technique used in this investigation was microscopic observation, which allowed the determination of the reactions occurring during solidification of the alloys, and possible reactions in the solid state.

Jan 1, 1968

By O. D. Gonzalez, R. A. Oriani

The migration of hydrogen and of deuterium dissolved in a iron and in nickel induced by an applied electrical potential has been measured over a range of temperature. In all cases the intevstitial solute migrates to the cathode. The deduced charges of transport are positive, independent of or slightly dependent on temperature, and markedly larger for the beavier mass. These results are interpreted as showing predominating momenlum transfer between electron holes and the acticated complex of the elementary diffusion act. THE migration of components of a metallic system under the action of an applied electrical potential is called electromigration or electratransport. This phenomenon has been known since 1861 and is becoming of increasing interest to metal scientists because of the opportunity it offers to explore the connections between mass transport and electron transport. Until recently, reliable data have been very scarce because of the difficulties inherent in the measurements. In addition, theoretical understanding has been either nonexistent or rudimentary and at present this situation is only moderately improved. Whereas the prevalent notion had been (and is still widely held) that electromigration directly reflects the state of ionization of a component in its equilibrium state in the metal, it is now becoming clear that the major, if not the sole, factor is the momentum transfer between the charge carriers, themselves in motion because of the applied field, and the metallic components in their activated state. It becomes important therefore to obtain reliable data over a range of temperature on electrornigration in systems where only one component moves, that is, either in pure metals (vacancy motion) or in interstitial solid solutions. Our choice of systems was dictated by our measurements' of the Soret effect in solid solutions of hydrogen and deuterium in iron and nickel. Indeed, the main reason for making measurements on electromigration at all is that we suspected a mechanistic connection between thermal diffusion and electromigration. In our former paper' we proposed that what is common to both phenomena is the momentum transfer between the charge carriers in the metal and the atom as it jumps from one equilibrium position to another. In Soret diffusion, however, there is in addition the dissipation of the energy of the jumping atom via pho-nons. We have since found that similar ideas for this interconnection have been developed by Fiks.2 In the present paper we will discuss in more detail our understanding of the mechanism of electromigration itself, but we turn first to a presentation of the experiments and results. EXPERIMENTS AND RESULTS Let us consider a metallic system in which only one component, m , can diffuse at the temperature in question, and to which are applied simultaneously an electrical potential gradient, grad cp, and a chemical potential gradient, grad p,. One may write3 for the material flux Jm and the electronic charge flux Je by applying the condition, grad = 0, to Eq. [I] and evaluating grad(um/T) for the isothermal case and for activity coefficient independent of composition over the relevant composition range. ere, c is the concentration and D the diffusivity of the mobile component. The experimental technique that we have used to measure Z* for the isotopes of hydrogen individually dissolved in each of two transition metals has been developed from a method used by Wagner and Heller.4 (A somewhat similar method was also used by Herold.5) A wire or thin rod of the metal under study connects two volumes of gaseous hydrogen at identical pressures; hence grad um = 0. A direct current is made to pass through the specimen, serving to maintain the specimen at the desired temperature and also to produce the electromigration of the hydrogen dissolved in the metal in equilibrium with the gaseous hydrogen. Under these conditions, the flux of hydrogen from one volume to the other is given, from Eq. [I] for grad( um/T) =0by

Jan 1, 1968

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By A. Raman, P. P. Singh

The constitution of the La-Rh system was studied by powder X-ray diffraction, metallopaphic, and differential thermal analysis techniques and an equilibrium diagram is presented. Eleven intermediate phases occur in the system and the crystal structural data for nine of them were determined. La3Rh crystallizes in an orthorhombic structure of undetermined type, whose unit cell is obtained by doubling the 'a; and 'c,,' edges of an FesC type unit cell. The other intermediate phases of the system are LarRh-3( undetermined structures also occur in the system. LaRh, undergoes a polymorphic phase transformation at 1240°C. LaRh3 and La2Rh7 also exhibit polymorphisnz. The phases Laah and LazRh7 melt congruently. The latter undergoes a eutectoid transformation into LaRh, and Rh at 1205°C. Laah3 is formed by a peritectoid reaction between Laah and La,Rh,,. The other Phases result from peritectic reactions between the liquid and the adjacent rhodium-rich phases. The intermediate Phases of the La-Rh system are compared with those of the La-Co and La-Ni systems. DURING the course of a detailed investigation to study the occurrence of CrB, FeB, A1B2, and related structures in the rare earth alloys it was found that much information is lacking for the rare earth noble metal systems. Although the structures of several rare earth alloys containing the noble metals at the AB and AB2 stoichiometries have been reported, the occurrence of related structures at other stoichiometries has not been studied. We have initiated a project to study the crystal structural features of selected rare earth-rhodium alloys and to map the equilibrium diagrams of representative systems with conventional methods. The results of our investigation in the La-Rh system are presented in this paper. Two phases were known in the La-Rh system. LaRh has the CrB-type structure.' LaRhz is a MgCu2-type Laves phase.z EXPERIMENTAL PROCEDURE Alloys weighing less than 1 g were prepared from commercially pure lanthanum (99.9 pct +), supplied by Lunex Company, Pleasant Valley, Iowa, and rhodium (99.92 pct +), supplied by Engelhardt Industries, Newark, N.J., in a conventional arc melting furnace under argon atmosphere. The buttons were turned upside down and remelted three times to insure homogeneity in the samples. Since negligible loss of material was encountered during melting, a chemical analysis of the alloy buttons was not undertaken. Powder specimens for X-ray diffraction studies in the as cast state were then prepared. The buttons were wrapped in thin molybdenum foils and homogenized by heating in vacuum at suitable high temperatures for more than 1 week. They were then broken into three or four pieces for annealing experiments. The pieces were wrapped in molybdenum foils and annealed at various temperatures in evacuated quartz capsules. The annealing was carried out for 2 hr at or above 1200°C, 1 day at temperatures close to llOO°C, 2 days at 1000°C, and for 1 week at temperatures below 1000°C. After annealing the alloy pieces were again broken and powder specimens for X-ray diffraction were prepared. The powders of the lanthanum rich alloys with more than 80 at. pct La were prepared by filing. The filings were sealed in molybdenum tubings and stress-relieved at 600°C in vacuum. It was not deemed necessary to stress-relieve the powders of the other alloys, since the alloys were very brittle and were ground easily. POWDER X-RAY DIFFRACTION X-ray diffraction photographs of powders (-325 mesh size) of the alloys in the as cast and annealed states were prepared in a Guinier-de Wolff focussing camera with copper K, X radiations. These patterns were studied to identify the stoichiometries and the crystal structures of the intermediate phases. The lattice parameters of the phases were calculated after minimizing the differences between the observed sin2 6 values, calculated from the diffraction angles 8, and the sin2 8 values, calculated using approximate lattice constants obtained from a few lines. These differences were minimized manually to less than 0.0005. The latLice constants are judged to be accurate to *0.005A for values less thp about 10A and to k0.01~ for values greater than 10A. The relative intensities of the lines were calculated using a computer program written by Jeitschko and Parthk.~ No attempt was made to refine the atomic positional parameters in the phases. METALLOGRAPHY The phase equilibria in the investigated alloys in the as cast and annealed states were also studied by metallographic examination. The polished specimen surfaces were etched with 10 pct picric acid in alcohol (alloys up to 25 pct Rh), concentrated picric acid (from 25 to 37.5 pct Rh), 2 pct nital (40 to 50 pct Rh), 10 pct nital (from 50 to 66.7 pct Rh) and with concentrated 48 pct HF for the other rhodium-rich alloys. Selected microstruture~ were then photographed using a Po-laroid Land camera. THERMAL ANALYSIS Differential thermal analysis of the alloys was carried out in DTA-668 Stone differential thermal ana-

Jan 1, 1970

By J. H. Swisher

The rate of internal sulfidation of austenitic Fe-Ti alloys in H2S-H2 gas mixtures is controlled primarily by sulfur diffusion, with counterdiffusion of titanium playing a minor role. At temperatures below 1100°C, enhanced diffusion along grain boundaries becomes important. The rate of internal sulfidation at 1300°C is approximately equal to the rate computed from the sulfur diffusion coefficient. The diffusion coefficient of titanium in y iron has been determined from electron microprobe traces in the base alloy near the subscale interface. The solubility of sulfur in Fe-Ti alloys has been measured in the temperature range from 1150° to 1300°C. The equilibrium sulfur content is found to increase with titanium content, due to the large effect of titanium on the activity coefficient of sulfur. The Ti-S interaction becomes stronger as the temperature decreases. TITANIUM as an alloying element in stainless steels is an effective scavenger for interstitial impurities, carbon in particular. Titanium is known to form stable sulfides; however extensive thermodynamic data on the Ti-S system are not available. Schindlerova and Buzek1 have shown that the Ti-S interaction in liquid iron is moderately strong. There have been no previous studies of the Ti-S interaction in solid iron. Internal sulfidation of Fe-Mn alloys was the subject of a recent investigation by Herrnstein.2 He found the rate of internal sulfidation to be an order of magnitude greater than predicted from available solubility and diffusivity data. A satisfactory explanation for the discrepancy could not be given. In the present study, the solubility of sulfur in austenitic Fe-Ti alloys was measured using a standard gas equilibration technique. Fe-Ti alloy specimens were also internally sulfidized. The rate of internal sulfidation was measured as a function of temperature and alloy composition. Supplementary electron micro-probe measurements were made to provide additional information on the nature of the internal sulfidation process. EXPERIMENTAL The starting materials were alloys containing 0.12, 0.24, 0.38, and 0.54 wt pct Ti. The alloys were made in an induction furnace by adding titanium to electrolytic iron that previously had been vacuum-carbon-deoxidized. The major impurity in the alloys as determined by chemical analysis was carbon. The carbon content of the alloys averaged about 100 ppm; metallic impurities were presented in concentrations of 50 ppm or less. Specimens were made in the form of flat plates, 0.03 by 2 by 4 cm for the equilibrium measurements and 0.5 by 1.5 by 3 cm for the rate measurements. The experiments were performed in a vertical resistance furnace wound with molybdenum wire and containing a recrystallized alumina reaction tube. In the gas train, flow rates of the reacting gases were measured using capillary flow meters. The source of H2S was a mixture of approximately 2 pct H2S in H2, which was obtained in cylinders from the Matheson Co. A chemical analysis was provided with each cylinder. The H2-H2S mixture was diluted with additional hydrogen to obtain the desired ratio of H2S to H2, and the resulting mixture was diluted with 30 pct Ar to minimize thermal segregation of H2S in the furnace. Argon was purified by passage over copper chips at 350°C and subsequently over anhydrone. Hydrogen was purified by passage over platinized asbestos at 450°C and then over anhydrone. The H2-H2S mixture was purified by passage over platinized asbestos and then over Pas. The samples used in the solubility measurements were analyzed for sulfur by combustion and iodometric titration. The subscale thickness in the internally sulfidized samples was measured on a polished cross section, using a microscope with a micrometer stage. Electron microprobe traces for titanium in solution were made on several samples that had been internally sulfidized. A Cambridge microanalyzer was used, and the known titanium content at the center of the specimen was used as a calibration standard. The procedure for the microprobe measurements will be described further when the results are presented. RESULTS AND DISCUSSION Equilibrium Data. Fig. 1 shows the sulfur concentration as a function of gas composition for three alloys equilibrated at 1300°C. The dashed line is based on data published by Turkdogan, Ignatowicz, and pearson3 for pure iron. The breaks in the curves are the saturation points for the alloys. The fact that the initial slope decreases with increasing titanium content indicates that titanium interacts strongly with sulfur in solution. To obtain information on the composition of the precipitating sulfide phase, the measurements described in Fig. 1 were extended to higher sulfur partial pressures. These results are shown in Fig. 2. (The initial portions of the curves are reproduced from Fig. 1.) The highest PH2s /pH2 ratio used is believed to be below the ratio required for the formation of a liquid sulfide phase. Time series experiments were used to study the approach to equilibrium in the samples. It was found that equilibrium with the gas phase was reached in less than 4 hr at 1300°C.

Jan 1, 1969

By Morris Muskat

A theory is presented for calculating the performance history of complete water-drive systems producing from idealized stratified formations. The general equations are applied to systems where the permeability stratification is either of the exponential or linear type. Calculations were carried through for different degrees of permeability stratification, but with special emphasis on the effect of the mobility ratio between the produced oil and the invading water on the resultant performance. These results are also expressed graphically as curves for the initial water breakthrough recovery, for the different degrees of stratification, as a function of the mobility ratio, and of the composition of the produced fluid stream as a function of the cumulative oil recovery. For several typical cases the latter has also been plotted as a function of the cumulative oil and water throughflow. The general result is that when the mobility of the oil is lower than that of the invading water the channelling tendency resulting from the permeability stratification becomes aggravated as the higher permeability zones become flooded out. Situations of this type would obtain when producing low gravity or highly viscous oils. Conversely, if the mobility of the oil is high compared to that of the invading water, the flooding of the high permeability zones will lead to a retarding and choking effect, and the gross bypassing phenomena will be partially suppressed. These conditions would correspond to those of flooding high gravity or low viscosity oils. A discussion is given of the various basic assumptions made in the analysis, including that of ignoring the stripping phase of the production history as implied by relative permeability concepts. INTRODUCTION The physical ultimate recoveries from oil reservoirs are basically determined and limited by the physical oil displacement processes associated with the reservoir producing mechanism. In practice, however, the economic ultimate recoveries are further limited by the mobility of the reservoir fluids and the uniformity and continuity of the producing formation. In fact, it is the differential depletion between the component parts of the composite reservoir which ultimately determines the total recovery at the time of field abandonment. While this observation applies to both the solution gas drive and gravity drainage mechanisms, in which use is made only of the energy contained within the original oil-bearing reservoir, it is of even more paramount importance under operations wherein the energy associated with extraneous fluids provides the ultimate oil expulsion mechanism. Whether the invading fluid is the water from an edgewater drive, water injected for pressure maintenance, gas injected for pressure maintenance, or gas returned to the formation in a cycling program, it is often the continuity and uniformity of the producing section which will control the economic efficiency of the operations. The importance of the problem of reservoir non-uniformity does not, of course, lessen its complexity or the difficulties of its solution. In fact, these are inherently such that the concept of a "general" solution is virtually meaningless. About all that can be reasonably hoped for is the analysis of specific and well-defined types of non-uniformity which may give some degree of approximation to actual reservoir conditions. Since variations in the nature of the reservoir which depend only on the position along the streamlines will not lead directly to major differential depletion development within the reservoir, the types of non-uniformity considered thus far have involved stratification assumptions. That is, the producing section has been replaced by a multi-layer "sandwich," each uniform areally, and differing from the others only in its basic physical constants as to thickness, porosity and permeability. The fluid motion in the composite system is thus approximated by a parallel superposition of the independent fluid movements in the individual strata. For the specific application to cycling operations the theory of the effect of permeability stratification has been developed for both discontinuous' and continuous types of permeability stratification. Among the latter, treatments have been given of systems in which permeability distributions are governed by exponential, linear' or probability3 functions. In all these studies complete dynamical equivalence was assumed between the injected dry gas and the displaced wet gas. The overall effective permeability of each stratum was therefore considered as constant and independent of the degree of invasion of the injected fluid. In the case of the displacement of oil by water, the assumption of dynamical equivalence between the water and oil will be strictly valid only by accident. Even if the oil viscosity should be the same as that of the water, the effective permeability to the oil in the presence of the connate water will in general be quite different from that of the water behind the water-oil interface flowing past the trapped residual oil. As a result the effective permeability for the stratum as a whole and rate of water invasion will change with time as the intrusion continues. The differential fluid motion in the individual strata will thus also vary with time. Qualitatively, it is easy to predict the resultant effects. If the permeability to viscosity ratio of the invading fluid exceeds that of the fluid displaced, the stratification and bypassing effect of the perme-

Jan 1, 1950