Search Documents

By Jun Hino, P. G. Shewmon, P. A. Beck

THE investigations of Crussard,' of Guinier and Tennevin,' and of Dunn and Daniels," indicate that the subgrains formed in a cold worked and annealed metal are capable of growing at each other's expense during annealing, if the temperature of annealing is sufficiently high and the time is long enough. The results of Dunn and Daniels are particularly convincing in showing that subgrain growth is essentially the result of the free surface energy associated with the subboundaries. In subgrain growth, as in ordinary grain growth, energy is gained as a result of the decrease in total subboundary surface area per unit volume. Recently, Wood and Scrutton4 found that the rate of subgrain growth upon annealing increased very considerably, if simultaneously a creep strain was applied to the specimen at a low strain rate. Working with 99.98 pct pure fine-grained aluminum, these investigators found that the continuous X-ray diffraction back-reflection rings of material strained at room temperature remained continuous after heating for ten days at 250 °C. However, when heated at the same temperature under a stress of 1000 psi even for only three days, the formerly continuous X-ray diffraction rings broke up into numerous dots, which were fairly clearly separated from each other. The continuous X-ray diffraction rings were interpreted as indications of a very small subgrain size, not resolved by the X-ray diffraction method used. The breaking up of the continuous X-ray diffraction ring during annealing under stress was taken as an indication of a great increase in subgrain size, so that the individual subgrains could then be resolved. The effect of simultaneous strain at a low strain rate in accelerating subgrain growth, discovered by Wood and Scrutton and designated by them as "cell growth," is of fundamental importance. The experiments described in this note were carried out in order to confirm Wood and Scrutton's results by direct metallographic observation. Also, information was sought as to the minimum creep strain necessary to produce this effect. A fine grained high purity aluminum strip was prepared by alternate 33 pct rolling and annealing treatments for 1 hr at 350°C. Specimens cut from this strip were subjected to a relatively fast creep strain of 7.2 pct in 3.5 min at 300°C under a constant load initially corresponding to 1185 psi. The subgrains set up were large enough (about 0.015 to 0.03 mm) to be clearly observed (Fig. 1) at X200 magnification with polarized light, after electrolytic polishing and anodic etching, producing a fine oxide film." The corresponding X-ray diffraction pattern is shown in Fig. la. After the fast creep strain treatment a portion of the specimen was subjected to a creep strain of 8.3 pct in 44 hr at 350°C under a constant load initially corresponding to 320 psi. Another portion, annealed under the same conditions, but not strained, served as control specimen. The subgrain size of a typical area of the specimen annealed under strain and of the one annealed without strain is shown in Figs. 2 and 3 (X-ray diffraction patterns, Figs. 2a and 3a). Annealing without strain produced clearly observable subgrain growth (subgrain size: about 0.03 to 0.05 mm). The effect of simultaneous strain was to increase greatly the rate of subgrain growth (resulting in a subgrain size of approximately 0.05 to 0.13 mm), in accordance with Wood and Scrutton. In another experiment, the effect of the amount of simultaneous strain was studied. The initial subgrain structure was set up by fast creep, as described previously. The specimen was then heated at 350 °C under a constant load initially corresponding to 320 psi, as above, for 2.2 hr (0.24 pct strain), and 8.8 hr (0.58 pct strain). Comparison with the corresponding unstrained control specimens by means of X-ray diffraction showed that 0.58 pct strain definitely had an effect, but the effect of 0.24 pct strain was doubtful.

Jan 1, 1953

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By Howard Martens, Pol Duwez

IN two papers published in 1949, alloys of chromium with the refractory metals tungsten, molybdenum, tantalum, and columbium were investigated in view of their possible use as high temperature resisting materials. For the Cr-Ta system, a partial phase diagram was presented and the only intermediate phase was identified at Ta2Cr3. A phase of the same composition was also observed in the Cb-Cr system. The X-ray diffraction data presented in these papers, however, were insufficient for crystal structure determination. It is shown in the present study that the only intermediate phase in both the Ta-Cr and the Cb-Cr systems corresponds to the ideal stoichiometric ratio TaCr2, or CbCr2. Both structures are cubic, MgCu, type. At high temperature, however, TaCr2 has a hexagonal MgZn, type structure, which can be retained at room temperature by fast cooling. The alloys were prepared by melting in a helium arc furnace on a water-cooled plate. The design of the furnace was essentially the same as that described in ref. 3. Some alloys were also obtained by sintering compacts made of the mixed powders pressed at 80,000 psi. The sintering was carried on for 4 hr at 1375°C. The tantalum and columbium powders were supplied by Fansteel Metallurgical Corp., North Chicago, 111. The tantalum powder was the reagent grade, with a particle size smaller than 400 mesh and a total impurity content less than 0.1 pct. The columbium powder was smaller than 325 mesh and contained approximately 0.1 pct C and traces of Fe, Ti, and Zr. The electrolytic chromium powder from Charles Hardy, Inc., New York, was smaller than 300 mesh and contained about 0.1 pct Na, 0.05 pct Ca, and traces of Cu, Al, Mg, Si, and Co. Powder diffraction patterns were obtained with a 14.32 cm camera, using copper Ka radiation filtered through nickel foil. The powder pattern of the TaCr2 alloy obtained by sintering at 1375'C was different from that obtained on the same alloy rapidly cooled from the melt. Contrary to this result, the powder pattern of CbCr2 was the same, whether the alloy was made by sintering at 1375°C or by melting, and was similar to that of the TaCr, sintered. It was also found that the structure of the TaCr2 specimen obtained by melting was retained after heating for 4 hr at 1590°C, but transformed into the structure found in the sintered specimen after heating for 4 hr at 1375°C. Hence, the structural change of TaCr2, appears to be a reversible polymorphic transformation. CbCr2 and ToCr2 Structure, Low Temperature Form By using large scale Hull-Davey charts, the powder pattern of CbCr, and of the low temperature form of TaCr2 were readily interpreted on the basis of a face-centered cubic lattice with a parameter of approximately 6.95 kX. The indices of the reflections together with the values of sin' 0 are given in Tables I and 11. From this list of observed reflections, it appears that the (200), (600), (024), (046), and (028) reflections are missing. The lack of (h00) reflections for h 4n indicates a four-fold screw axis. The missing (Okl) spectra for k + 1 An indicate the existence of a diamond glide d. The combination of these symmetry elements can be found in the O— Fd3m space group, which is therefore the most probable one. After having determined the approximate density of TaCr, by the immersion method, the number of molecules per unit cell was calculated and found to be nearly eight. This information, added to the fact that the most probable space group is O leads to the consideration of a structure of the MgCu2 type, in which the atoms have the following positions: 8 magnesium in a and 16 copper in d. On the basis of this structure, intensities were computed by means of the usual formula: 1 cos'20 I a sin2 cos where F is the structure factor; 8, the Bragg angle: and p, the multiplicity factor. As shown in Tables I

Jan 1, 1953

By P. R. Celmer, Leonard R. Weisberg

THE principle of fractional crystallization has been successfully used to prepare high-purity (99.999 pct) Ga. Hoffman and Scribnerl removed single crystals of gallium solidifying in a gallium melt, while Zimmerman2 grew single crystals from the melt by the Kyropolous technique. In contrast, attempts at purifying gallium by zone refining have been less successful. ichards, reported that despite the passage of 40 zones through a gallium ingot, there still remained 5 to 70 ppm each of Cu, Fe, Ca, Mg, Si, Al, and Ag. Previously, Detwiler and Fox4 detected only one impurity, Pb, segregating in zone-refined ingots. These surprising results prompted an investigation of the factors controlling impurity segregation in gallium. Possible reasons for this were insufficient diffusion of impurities in the melt; recontamination of the melt by its oxide film which is not affected by the passage of the zone; reaction of gallium with the boat; sudden freezing of gallium following supercooling, especially since gallium easily supercools to and trapping of impurities at grain boundaries. Impurity segregation tests were carried out by directionally freezing gallium using the Bridgman rnethd, modified in that the molten gallium is lowered out of a furnace into a slush of dry ice and trichloroethylene, thus minimizing supercooling. Since the gallium is contained vertically, the oxide film is in contact only with the tail end of the melt. It was found that Teflon makes an excellent crucible for gallium since it is quite pure, non-reactive, translucent, flexible, machinable, and is not wet by gallium. The gallium crystals could be grown at various speeds, and the melt could be vigorously stirred by a Teflon rod moving through the gallium in a vertical reciprocating fashion. Single crystals could be grown by placing a solid Teflon plug at the bottom of the melt drilled out in such a way to cause the solidifying gallium to follow a winding path. Thus, even if many crystals are originally nucleated, only one grain will predominate. The grain structure of the gallium crystals was revealed by an etchant composed of equal volumes of HC1, HNO3 and HF, diluted with water to half strength. Emission spectrographic analyses were carried out on samples removed from the front and tail ends of the resulting gallium ingots. Typical results of this study are summarized in Table I. The rate of freezing in all three cases was about 1 in. per hr. It can be seen that even though stirring of the melt does help, it is even more important to grow a single crystal of Ga in order to obtain good segregation of impurities. The effect of crys-tallinity on the segregation of impurities was previously observed6 in the directional freezing of germanium; however, in this case, the effect was much less pronounced. This dependence of impurity segregation on the crystalline perfection of Ga may be related to its thermal conductivity which is the most anisotropic of all metals.7 The anisotropic thermal conductivity can cause the solid-liquid interface to be nonuniform, thus leading to trapping of impurities during freeing, and therefore reduced segregation. In conclusion, it is indicated that zone refining of gallium would be more successful if seeding and similar precautions are taken to insure single crystal growth. The authors are indebted to Mr. H. H. Whitaker for the spectrographic analyses and to Drs. B. Abeles and F. D. Rosi for helpful advice and encouragement throughout the course of this work. This research was supported by the Electronics Research Directorate, Air Research and Development Command, under Contract No. AF33(616)-5029. REFERENCES 'J. 1. IToffmon and B. T. Scribrer: I. Research h'atl. Bur. Standards, 1935, "01. 15, p. 205. 'W. Zirnmerman: Science, 1954, vol. 119, p. 41. %J. L. Richards: Nature, 1956, vol. 117, p. 182. 'D. P. Detwiler and W. M. Fox: I. Metals, 1955, "01. 7, p. 205. 5P. W. Bridgman, Proc. Am. ilcod. Sci., 1925, vul. 60, pp. 305,385,423. "S. L1. Christian: private c ommuni cation. 'K. W. Powell: roy. Sac., 1951, vol. 209, p. 525. 'W. G. Pfann: Zone Re fining, p. 20. John Wiley and Sons, Inc., New York, L058.

Jan 1, 1962

By M. R. Tek, K. H. Coats, D. L. Katz

The nature and the limits of validity of Darcy's law US applied to the flow of natural gas through reservoirs has been considered in order to resolve some controversial aspects of the effect of turbulence on pressure drops. The equivalence between various concepts and viewpoints advanced in the past by several investigators to explain how and why a gas well does not necessarily perform according to Darcy's law is shown. Starting with generalized equations of flow of fluids through porous media, a partial differential equation has been derived which accurately represents the flow at all rates. This equation has been numerically solved using an IBM 704 digital computer. The results permit plots of unsteady radial pressure distribution curves from which specific isochronal backpressure curves may be constructed. These back-pressure curves show the effect of the ,8 factor on the slope of the back-pressure curve. The calculations further indicate that the drainage radius for a gas well in turbulent flow propagates at a rate dependent upon the rate of production at the wellbore. This is quite different from the case with liquid flow or natural-gas flow in laminar regime. Additionally, the effect of reservoir inhomogeneities and crossflow between layers of different permeability on the back-pressure performance of gas wells has been conridered. In light of the current numerical results the significance and linlitation of the rate of flow function Y proposed by Smith" has been discussed. INTRODUCTION The relationship between the pressure drop and flow rate in problems of fluid flow through porous media is known to be affected by the nature of flow through the porous matrix. It has been observed by many that, for a range of flow rates, the pressure drop remains proportional to the rate of flow. When some flow rate is reached, however, it is usually observed that the pressure drop gradually begins to increase more than proportionally to the flow rate. It is well known that this phenomenon was first observed by Osborne Reynolds in 1901 in experimenting with flow through pipes. In his classical experiments, Reynolds made visual observations on the condition of streamlines evidenced by injecting a dye into water flowing through glass tubes. In these experiments, the abrupt transition between steady, "streamline, laminar" flow and unsteady random turbulent flow was found to be a function of the dimensionless group (Dvp/u), now known as the Reynolds number. During these experiments, in addition to observations on the nature of flow regimes, the proportionality between flow rate and pressure drop in laminar flow was contrasted with the nonlinearity between these variables in turbulent flow. Fancher and Lewis' reported data on various consolidated and unconsolidated sands in 1933. Their conclusions were that ". . . the flow of fluids through these porous materials closely resembles that through pipes; that there is a condition of flow in porous systems which resembles viscous flow, another which corresponds to turbulent; that the change from one type to the other takes place at a definite and reproducible condition for each system". In 1947, Brownell and Katz published a method to predict the laminar and turbulent flow behavior from the particle size, bed porosity and the particle sphericity, employing the friction factor-Reynolds number charts for pipes. Several investigators have verified the work of Fancher, Lewis and Barnes and presented their data as friction factor-vs-Reynolds number plots.' The equation which would represent the pressure gradient over the whole range of velocity must have an added term over that represented by Darcy's law. Accordingly, the pressure gradient necessary to sustain flow at the velocity (v) through a porous medium may be represented by the following equation, suggested by Forscheimer." The nature and the range of validity of Darcy's law has been the subject of studies by many investigators over the past years. While everyone seemed to agree on the need for a quadratic correction term to Darcy's law to make it effective over the range of velocities, the concept of inception of turbulence and the use of the term "turbulent flow" remained controversial. Some fluid dynamicists define turbulence as a flow

By Byoung Whie Lee

In order to investigate the change in surface energy as a function of cold work ad environment, the contact angles between a liquid and the surface of elec-tropolislzed aluminum with various degrees of cold work were measured. Specinlens cut from aluminum rod (99.99 pct Al) were annealed at 600°C for 24 hr and then cold compressed parallel to the drawing di-rectzon of the rod in increments of approximately 10 pcl from 0 to 70 pct. The sections, taken transversely to the drawing direction, were mechanically and electro chemically polished. Contact angles between the sessile drop of distilled water, glycerin, and glycol and the specimen surface were measured optically. The differences in interfacial energy were calculated from the contact angle data. The effects of chemical washing and surface roughness of the specimen on the contact angle were also studied. The liquid contact angle was a linear function of the degree of cold working (D) ad changed markedly with differed liquids. THE plastic deformation of polycrystalline metals induces a preferred orientation in crystal structure. A preferred orientation changes the atomic density of a metal surface. This change in surface atomic density brings about the changes in surface interactions of the metals."' Also, the plastic deformation of metals raises the internal energy by increasing the stored energy.' Consequently, the total Helmholtz free energy of the metal increases and, in turn, the surface free energy of the metal. In order to study the effect of plastic deformation on solid surface phenomena, the liquid contact angle of a sessile drop made on cold compressed and elec-trochemically polished aluminum was measured. The preferred orientation in aluminum, due to uniaxial compression, is known to be [110] direction parallel to the compression axis. However, this description is only a first approximation of a rather complex orientation. In actuality, half the crystallites are more than 10 deg from this orientation, regardless of the amount of compression." For randomly oriented fcc crystals, the average surface atomic density can be calculated as follows: The number of atoms in a unit volume is 4/a3, where a is a lattice constant. Hence, the number of atoms per unit surface area is about (4/a3)2/3 = 2.52/a2 However, when the preferred orientation took place along [110] direction, the surface atomic density is as follows: The area of (110) plane of a unit cell, considering ridges and valleys, is BYOUNG WHlE LEE, Member AIME, is Research Metallurgist, Atomic Energy Research Institute of the Republic of Korea, Seoul, Korea. Manuscript submitted May 29, 1968. IMD 2 x J2- a x Ga = 6a2 Since this plane contains 4 atoms, the surface atomic density is 4i6a2 = 2.3/a2 Upon changing from random to preferred orientations, the surface atomic density decreases about 10 pct. The increase in the stored energy and the change in surface atomic density due to plastic deformation would certainly alter the surface tension of a metal. EXPERIMENT High purity aluminum rod of 6-mm diam (99.99 pct Al, obtained from Johnson, Matthey & Co., Ltd, London, England) was cut into sections approximately 1 cm long. The transverse cross section was mechanically polished. Then, the specimens were annealed in air at 600°C for 24 hr. Thus, eight annealed specimens were uniaxially compressed at room temperature from 0 to 70 pct of cold working with an increment of 10 pct. The extent of cold working was determined by the percentage of reduction in axial length to the original length. The specimens were further polished mechanically and, then, electrochemically with the Modern Battelle Type 1 electropolishing solution with a current density of 0.3 to 0.4 amp per sq cm and an electrode potential of 8 to 12 v at the solution temperature of 78" to 82°C for approxirr~ately 5 min. The Modern Whlte Diffusing Screen \t/ --------------- re Goniometer for Mercury Lamp ^-^ Sessile Drop ¦=• Condenser c=> Microscope l_*~_i bJ Mlcromster Filer with Protmcter Fig. 1—Schematic diagram of apparatus set-up for liquid contact angle measurement.

Jan 1, 1970

By G. L. Simard, D. J. Salley, J. Chupak

DITHIOPHOSPHATES and xanthates are the principal collectors for sulphide minerals, and consequently any knowledge of mineral-collector systems of this type is of value. In the present investigation an attempt was made to obtain information on the interaction of a typical dithiophosphate with galena. In carrying out this study, radioiso-topes, so much discussed in the past few years,', ' were extensively employed. By using radioactive dithiophosphate synthesized from radioactive phosphorus, a rapid and sensitive analytical procedure G. L. SIMARD is in the Research Division, Stamford Research Labs., American Cyanamid Co., Stamford, Conn.; J. CHUPAK, formerly in the Research Division, Stamford Research Labs., is now at Camp Detrick, Frederick, Md.; and D. J. SALLEY is in the Research Division, Stamford Research Labs. AlME New York Meeting, Feb. 1950. TP 2815 B. Discussion (2 copies) may be sent to Transactions AlME before April 30, 1950. Manuscript received Oct. 17, 1949. was at hand. This permitted determination of such important quantities as the rate of uptake of dithiophosphate by the mineral, the amount existing at equilibrium on the mineral surface and in the solution, and the desorption of the agent from the surface. Such measurements were possible even on single crystals of galena. In addition the exchange of dithiophosphate between the solution and the sorbed phase was examined, a matter which could be accomplished only by the use of isotopes. It is desired to point out at this time that throughout this paper the term "sorption" has been employed to designate the uptake of agent by mineral, without implication as to the nature of the process by which the uptake was accomplished. Experimental Materials: Galena: The galena was from the Tri-State district. For studies on ground mineral four preparations were used during the course of the investigation. These were prepared from selected large crystals by wet grinding in order to reduce surface oxidation. Alcohol was chosen as a convenient medium for this purpose. The ground mineral was then fractionated by sedimentation in alcohol, dried by evacuation, and stored in a nitrogen-filled desiccator. Data on the preparations are tabulated in table I. The size analysis indicated that the areas of the preparations were of comparable magnitude, even though the absolute values may be somewhat incorrect. Dithiophosphates: Nonradioactive dithiophbsphate (di-isopropyl or di-secondary butyl) was obtained by purification of a commercial product. An aqueous acid solution of the agent was extracted with petroleum ether, the ether layer dried, and the dithiophosphate precipitated as ammonium salt with anhydrous ammonia. Several repetitions of this process resulted in a nearly colorless, flaky product of good purity. (General formula (RO,)PSSNH,). Radioactive dithiophosphate was synthesized by heating radioactive elementary red phosphorus* • For early experiments (1943 to 1944), the elementary red phosphorus was obtained from the cyclotron group at the Crocker Radiation Laboratory of the University of California. Berkeley, through the courtesy of Dr. Joseph W. Hamilton. For more recent work (1946 to 1949). the elementary radioactive Phosphorus has been supplied by the Oak Ridge National Labarotories on allocation from the U. S. Atomic Energy Commission. with sulphur at 270" to 300°C to produce radioactive P,S,. The latter was then treated with the desired purified alcohol (isopropyl or secondary butyl) at 60" to 80°C to form the dithiophosphoric acid derivative. Purification was effected in the same manner as for the nonradioactive material. The effectiveness of the purification method was established by the isotopic dilution method." The same technique was used to show that only negligible decomposition of neutral or of carbonate solutions of the dithiophosphate took place over a period of a day; this was true whether or not galena was suspended in the solutions. Procedures: Radioactivity: The activity of a solution was obtained by counting with a small glass-jacketed, silvered Geiger counter, using a conventional scaling circuit.' Crystals and other solids were counted under a bell-shaped, mica-window counter tube. Sufficient counts were made so that the probable error of counting was of the order of +I-3 pct. The specific activity in terms of counts per minute

Jan 1, 1951

By Lyman H. Hart

Some of the significant facts that will affect the supply and demand for metals during the next few decades are given in this presentation. This is important because the only hope for intelligent guidance of the national destiny depends upon presenting a consistent and accurate picture of basic problems. The metals problems now fomenting are only slightly less important than those precipitating the energy crisis. Admittedly, it is difficult to generate serious thinking on this subject when the world is enjoying easy access to all metals at modest prices. This is possible because of the current world capacity to over-produce, and United States credit standing, at least up to now, has permitted buying freely on world markets. Although both these vital conditions are subject to change, there is little doubt that if the U.S. financial house is kept in order, we can muddle along for a while on a generally expanding, metal-based economy. But should our credit standing be permitted to collapse, or should there be a tightening of world supplies of metals, the nation would immediately find itself in trouble. Unfortunately, clear signals indicate that both of these conditions are threatening. The metal age started many centuries ago, and with the exception of the noble metals, barely a dent was made in world reserves until the present century. Then, as a result of the pioneering ideas of D. C. Jackling and Paul Gemmel, and the willingness of the McNeals, Penroses and Guggenheims to back ideas with money, a virtual explosion of metal production was propagated throughout the world. It is important to examine the causes of this expansion to determine whether it is a trend that may be sustained or whether it is anomalous and therefore attributable to circumstances not likely to be repeated or continued. This period of plenty is a once only phenomenon. Vast metal resources have been sitting around here for eons, while the world was unaware of or unconcerned with their existence. But when population began to grow exponentially, and ways were discovered to use metals to develop machines to do the work of men and animals, the exploitation of raw materials increased at a staggering rate. And now the realization is dawning that minerals are a wasting asset and are beginning to present supply problems that should begin to be of concern. Two main factors have caused this condition; One is the population growth, and the other is the expanding demand for irreplaceable raw materials as nations become more developed. A group of systems analysts, at Massachusetts Institute of Technology in a book Limits of Growth, has warned that at present rates of population growth and resource consumption together with environmental constraints, a point will be reached within 100 years when the world population will become totally incapable of supporting itself. The authors conclude that this end result could be avoided if a planned equilibrium could be established at a lower rate of population growth and resource consumption. Fortunately, there is evidence that the population problem is moderating. The growth rate in the U.S. has decreased from 2.3% per year in 1900 to 1.1% in 1970. But, even at the lower rate, population will double in 70 years. The world rate is greater, and if it continues, world population will rise from 3.5 billion to nearly 10 billion in 50 years. All of these projections are questionable because the dangers of growing population are now being recognized. A marked decrease in growth rate has been noted in the more developed countries-in fact, it was zero in the U.S. in 1972. The other factor affecting consumption of raw materials is the increase in per capita use of most items as individuals become more affluent. Many economists believe that developing countries such as China and India represent a vast future potential market; but there are others, with whom the author agrees, who believe over-population tends to stagnate development. Nevertheless, it is a fact that vigorously developing societies increase consumption at a rate far above their population growth rate.

Jan 1, 1974

By Laurance S. Reid, Joe A. Porter

Gas dehydration plays an important part in the production of natural gas. Effective dehydration prevents formation of gas hydrates and the accumulation of water in transmission Systems,2,6,7 insuring uninterrupted gas deliveries at maximum efficiency under the most adverse weather conditions. At the present time, most gas companies require a maximum water vapor content of seven lb per million standard cu ft of gas. so that virtually all gas tendered for sale must he dehydrated to meet this specification. For a number of years it has been common practice to produce gas and gather it at a common point for dehydration prior to discharge into the transmission system.1,3,11,16,17 However, higher transmission line pressures, long gathering lines and relatively low ground temperatures have made it necessary to dehydrate gas at. or near. individual. wells in order to gather gas from a number of newly developed fields without unusual difficulty. Where gas has been dehydrated at pressures ranging from 300 to 800 psi in the past. future trends indicate that these processes may be operated at pressures as high as 2,000 psi. Economics of gas dehydration are of great importance, partitularly where facilities must he provided to process relatively small quantities of gas, such as the production from an individual well. Although the adsorption of water vapor from gas on a granular sorbent material such as activated bauxite, activated alumina, or one of the alumina-silica gels is highly effective and produces virtually "bone dry" gas, the cost of a small unit of this type is substantially greater than that of an absorption process which, through proper selection of the absorbing liquid, will dehydrate the gas sufficiently to meet pipe line specifications. For this reason, a great deal of emphacis has been placed on the development of small, inexpensive dehydration units8,24 and the search for more effective absorb. ent liquids has been intensified. A wide variety of methods for dehydrating gas are known' and many of these have been used in industry. Earlier applications of the absorption process employed concentrated solutions of calcium and lithium chlorides as the absorbent. The severe corrosion problems inherent in handling these solutions and the relatively small dew point depressions obtained caused early abandonment in favor of, or conversion to. diethylene glycol when it was found that aqueous solutions of this organic liquid were more hygroscopic than the brines and were non-corrosive. Processes employing diethylene glycol-water solutions are widely used for gas dehydration at pressures ranging as high as 1,200 psi.13,14,15 At nominal pressures a dew point depression of 45° to 50°F may be be and the data of Russell et al." indicate that a minimum dew point is obtained from the effluent gas at a pressure of approximately 1,200 psi when the gas is in equilibrium contact with a 95 per cent by weight diethylene glycol solution. In a number of instances the dew point depression obtained with diethylene glycol-water solutions is not sufficient to produce a specification product without cooling the inlet gas. In a recent search for a better absorhent, triethylene glycol was used in a small commercial dehydration unit and subjected to rather exhaustive field tests.' The data obtained were encouraging and indicated that, at pressures ranging from 300 to 500 psi, triethylene glycol porduced a substantially greater dew point depression than diethylene glycol. These results led to an investigation of the system natural gas-water-triethylene glycol in an effort to obtain vapor-liquid equilibrium data, to determine pressure limitations, and to develop other data pertinent to the design of gas dehydration processes. A review of the literature has failed to reveal any data which permit reasonably accurate calculation of the vapor-liquid equilibrium conditions for a solution of water and triethylene glycol in contact with natural gas at high pressure. Since these constituents form a non-ideal system, the Poynting equation18,21 or the usual combination of Raoult's and Dalton's laws19,20,22 would not be valid. Correction of Raoult's and Dalton's laws by the use of activity coefficients2' is not feasible for available data are insufficient for the prediction of the actual increase in the ratio of the activity of one component in the vapor phase to its activity in the liquid. Therefore. experi-

Jan 1, 1950

By M. R. Tek

A new method for correlating the data on multiphase flow through vertical pipe is presented. The correlation is based on a "two-phase f factor" concept which was developed and successfully applied to horizontal multiphase flow by Bertuzzi, Tek and Poettmann. Field data previously published on several flowing and gas-lift wells have been used as the basis for the developed correlation. The application of the method to actual two-phase flow problems indicate that this method is capable of predicting the pressure distribution in vertical multiphase flow strings well within the accuracy range usually desired for common engineering and design calculations. A new working chart developed for calculation of two-phase pressure gradients and a graphical step-by-step procedure for the computation of pressure distribution are presented along with an example problem. INTRODUCTION Multiphase flow through vertical pipe is encountered in many engineering installations. In petroleum, chemical process, nuclear engineering and many other industries, problems associated with simultaneous flow of two or more phases through vertical pipe have been of interest for a long time. This interest has increased considerably during recent years due to applications to new processes in petroleum production and refining and to problems of steam generation and heat removal from nuclear reactors. One prominent example of vertical two-phase flow is provided by the gas-lift process where oil, water and gas flow simultaneously. If the pressure profile in a gas-lift well can be predicted within reasonable accuracy, it would be possible to get good estimates of the power required to lift the oil, the optimum depth, the pressure and the rate at which to inject gas. Furthermore, the effect of production rate and tubing sizes on these quantities can be evaluated before any design decision is made on the installation and operation of the flow string. The majority of experimental work available in the literature deals with two-component systems where individual phase flow rates in and out of the pipe remain constant. The general problem of prediction of pressure drops in multiphase flow systems is very complicated. The co-existence of numerous flow patterns of widely different geometry and mechanism, conditions of surface instability and the nature of force fields acting upon the system are among the major difficulties commonly encountered. The classical approach of fluid dynamics which would be based upon the formulation and solution of Navier-Stokes equations has been found by many investigators completely devoid of any hope of success-—not only because of inherent nonlinearities but also because of insurmountable analytical difficulties standirng in the way of setting up the boundary conditions. The presence and effect of interfacial forces on mu1tiphase flow systems further complicate the theoretical approach. For these reasons, many investigators choose to adopt semi-empirical if not purely empirical approaches in order to obtain solutions of engineering utility. A prominent practical solution has been given by Poettmann and Carpenter' in the form of an empirical correlation. In their paper, total flowing densities of fluids and solubility effects of gas in oil have been taken into account for the correlation of field data covering a wide range of operating conditions for oil wells. They treated the gas. oil and water as a single phase of combined propcrties and correlated the multiphase friction factor as a function of the product (pvd), i.e., density X velocity X diameter of the flow string. However, because the product (pvd) is dimensional while the frictior factor is dimensionless, the generality of their result's is somewhat restricted. It seems that the omission 01 viscosity effect may be one of the reasons for the scatter of data as shown in their correlation chart. The concepts of the "two-phase f factor" and the "two-phase Reynolds' number function" were recently developed and successfully applied to correlate horizontal multphase flow by Bertuzzi, Tek and Poettmann.' Recently, two new methods of correlation by Tek and Chan' have been presented on simultaneous flow of liquid and gas through vertical pipe. These two methods of correlation and the working equations in field units necessary for their application are included in this paper. The extension of the conccpt of two-phasc Reynolds' number function successfully developed for horizontal flow into vertical multiphase flow systems, development and evaluation of working charts permitting calculation of two-phase pressure gradients and

By J. P. Hirth, A. H. Clauer

A simple graphical method is presented for the determination of climb forces on dislocations exerted by uniaxial stresses. In conjunction with standard stereographic projections, the technique is applicable for any crystal system. Climb forces on glide dislocations in bcc crystals are determined as an application. In connection with a study1 of elevated temperature creep of molybdenum single crystals, a detailed analysis of the climb force on a dislocation was necessary. In general, such forces can be calculated for an arbitrary dislocation segment from the Peach-Koehler equation2 (or one of its variants3"6) where Fk are the components of the force per unit length of dislocation line, ?i of the unit sense vector tangent to the line, ojl of the stress tensor, bl of the dislocation Burgers vector, and is the Einstein permutation operator with components ?123 = ?231 = ?321 = —?132 = — ?321 = —?—213 = unity, and all other components zero. It is understood in the suffix notation of Eq. [I.] that repeated indices are to be summed over 1, 2, and 3. For the most general coordinate system, the matrix manipulations required to solve Eq. [I] are fairly complex. However, a simplification can be achieved by specifying a coordinate system based on the dislocation in question. The pertinent system involves cartesian coordinates xi, x;, xi, with attendant unit vectors e'i; ei being parallel to the Burgers vector b, and e'3; being parallel to the climb directions (b x ?), i.e., normal to the slip plane. In this coordnate system b = (bi, 0, 0) and ? = (?'1,?'2, 0). Hence, Eq. [I] yields for the climb force per unit length F3 = ?'2s'11b'1—?'1s'12b'1 [2] The first term on the right side of Eq. [2] is the climb force on the edge component of the dislocation, while the second is the (cross-slip) force normal to the glide plane acting on the screw component. Thus, the problem of determining the climb force reduces to that of resolving the stress components —il and s'12 in the xi coordinate system. Several authors7- l0 have shown, for the glide case, that such stress resolution can be performed most readily by graphical methods. Therefore, we present a similar method for the climb case. Of the suggested methods, we follow that of Hartley and Hirth,9 which can be performed directly with a standard stereographic projection for the crystal system in question. We treat the case of simple tension, or compression, this being the most common creep test in which climb is important. A set of cartesian coordinates (x1, x2, x3) is fixed with x3 parallel to the specimen axis. This set is connected to the x'i set by the factors defined in Table I. A perspective view of the coordinate systems relating them to the glide ellipse, is given in Fig. 1, while a stereographic projection of the coordinate systems is presented in Fig. 2. The transformation matrix connecting the xi and xi coordinates is x'i =aijXj [3] where an a11 a13 aij = a21 a22 a23 a31 a23 [4] s The corresponding transformation for the stresses is i [5] The only nonzero component of okl for the simple tension case is s33,. Hence, using Eqs. [4] and 151, we find that Eq. [2] becomes F3 = b'1s33 sin2 cos ?[cos ? sin a - sin ? cos a] [6] = m1b'1s033 sin a- rn2b's33 cos a where sin a = 5: and cos a = ?i, with a the angle measured from b to ?, i.e., a polar angle in the x1, x2 plane. The same equation with a sign change applies for compression. A graphical plot of the "Schmid" factor m1, relating to the climb force on the edge component, is presented in the stereographic projection of Fig. 3. In this projection the glide plane normal points toward the top of the page and the glide direction points normal to the plane of the page. To use the graph, one

Jan 1, 1970

By E. S. Machlin

An extension of the thermodynamic theory has been made for the case of irreversible growth processes occurring by the motion of an interface. The theory is applicable to such diverse phenomena as diffusion, growth in recrystallization, continuous grain growth, growth of carbides, etc., growth of eutectoid products, growth in solidification, recovery, "slipless" flow, etc. THE publication of a recent book1 has served to focus attention on a very powerful means of treating certain irreversible processes. This method which has been formally described as the thermodynamic theory of irreversible processes is applicable to processes which involve an approach to equilibrium and for which the deviations from equilibrium are small. The irreversible processes can be classified into three groups: chemical reactions, transport processes (diffusion, heat and electrical flows), and relaxation phenomena involving a degradation of internal energy to more stable states. The limits of applicability of the theory can be precisely defined1,2 for each type of irreversible process. For example, in a chemical reaction at constant temperature and pressure, the Gibbs free energy released per mol in the process should be less than the thermal energy, RT. Another type of irreversible process, which sometimes involves a combination of the first two of the above-mentioned groups, is treated in this paper. This process is one characterized by the motion of an interface separating two regions having different values of free energy and may be briefly described as a growth process. Prigogine2 and Herring3 have treated special cases of this type of process. Previous conscious applications of the theory in metallurgy have been limited to the field of diffusion, one of the transport processes. Darken,4 Bardeen,5 Prigogine2 and others have made significant contributions in this respect. Some other applications of the theory are described in this paper. The Theory The thermodynamic theory of irreversible processes is based on the work of Onsager,6 DeDonder,7 Prigogine2 and others. A resume of the theory has been given by Prigogine2 and DeGroot.1 A brief description of the theory follows, although for a complete understanding the reader is urged to read the references. It is first assumed that the change in entropy, even for a system removed from equilibrium, is given by dS = [dU + pdV - Si µi, dn, - Sk Pk dxk]T-1 Using this expression, the irreversible entropy production in the system is calculated by subtracting the contribution to the change in entropy of the system by transfer of heat, work, or matter from the surroundings. Thus, for example, the irreversible entropy production of a system at constant temperature and pressure is given by dtS = — — dG (DeDonder7) Now it has been shown by Prigogine,2 DeDonder,7 and DeGroot1 that the rate of irreversible entropy production diS/dt, when calculated in the manner suggested, can be written as a sum of products of conjugate forces XK and fluxes Jk, i.e., d,s/ = Sk Jk Xk [I-1] dt The assumption is next made that a linear relation exists between a given flux and the forces, obtained from Eq. I-1, i.e., Ji = Sk Lik Xk (i = 1, 2, . . . n) [I-2] Generally, there is a degree of freedom in choice of the fluxes and forces. However, for all the choices consistent with Eq. I-1 the Onsager8 relations, based on the principle of microscopic reversibility, namely, Lik = Lki, are valid. The coefficients Lik, depending upon the choice of the conjugate fluxes and forces, may be more or less dependent upon the parameters defining the state of the system. In general, the proper choice of fluxes and forces utilizing Eq. I-1, makes the coefficients Lik independent of time. Also, because many of the three classes of irreversible phenomena have already been treated,1,2 a proper choice can be made by analogy. In any case, the validity of Eq. 1-2, from whatever choice made using Eq. I-1, must be tested by experimentation. The theory, thus, comprises the following steps: 1—The irreversible entropy production is calculated to yield Eq. I-1. 2—The terms in Eq. I-1 are grouped

Jan 1, 1954

By J. Chupak, D. J. Salley, G. L. Simard

DITHIOPHOSPHATES and xanthates are the principal collectors for sulphide minerals, and consequently any knowledge of mineral-collector systems of this type is of value. In the present investigation an attempt was made to obtain information on the interaction of a typical dithiophosphate with galena. In carrying out this study, radioiso-topes, so much discussed in the past few years,', ' were extensively employed. By using radioactive dithiophosphate synthesized from radioactive phosphorus, a rapid and sensitive analytical procedure G. L. SIMARD is in the Research Division, Stamford Research Labs., American Cyanamid Co., Stamford, Conn.; J. CHUPAK, formerly in the Research Division, Stamford Research Labs., is now at Camp Detrick, Frederick, Md.; and D. J. SALLEY is in the Research Division, Stamford Research Labs. AlME New York Meeting, Feb. 1950. TP 2815 B. Discussion (2 copies) may be sent to Transactions AlME before April 30, 1950. Manuscript received Oct. 17, 1949. was at hand. This permitted determination of such important quantities as the rate of uptake of dithiophosphate by the mineral, the amount existing at equilibrium on the mineral surface and in the solution, and the desorption of the agent from the surface. Such measurements were possible even on single crystals of galena. In addition the exchange of dithiophosphate between the solution and the sorbed phase was examined, a matter which could be accomplished only by the use of isotopes. It is desired to point out at this time that throughout this paper the term "sorption" has been employed to designate the uptake of agent by mineral, without implication as to the nature of the process by which the uptake was accomplished. Experimental Materials: Galena: The galena was from the Tri-State district. For studies on ground mineral four preparations were used during the course of the investigation. These were prepared from selected large crystals by wet grinding in order to reduce surface oxidation. Alcohol was chosen as a convenient medium for this purpose. The ground mineral was then fractionated by sedimentation in alcohol, dried by evacuation, and stored in a nitrogen-filled desiccator. Data on the preparations are tabulated in table I. The size analysis indicated that the areas of the preparations were of comparable magnitude, even though the absolute values may be somewhat incorrect. Dithiophosphates: Nonradioactive dithiophbsphate (di-isopropyl or di-secondary butyl) was obtained by purification of a commercial product. An aqueous acid solution of the agent was extracted with petroleum ether, the ether layer dried, and the dithiophosphate precipitated as ammonium salt with anhydrous ammonia. Several repetitions of this process resulted in a nearly colorless, flaky product of good purity. (General formula (RO,)PSSNH,). Radioactive dithiophosphate was synthesized by heating radioactive elementary red phosphorus* • For early experiments (1943 to 1944), the elementary red phosphorus was obtained from the cyclotron group at the Crocker Radiation Laboratory of the University of California. Berkeley, through the courtesy of Dr. Joseph W. Hamilton. For more recent work (1946 to 1949). the elementary radioactive Phosphorus has been supplied by the Oak Ridge National Labarotories on allocation from the U. S. Atomic Energy Commission. with sulphur at 270" to 300°C to produce radioactive P,S,. The latter was then treated with the desired purified alcohol (isopropyl or secondary butyl) at 60" to 80°C to form the dithiophosphoric acid derivative. Purification was effected in the same manner as for the nonradioactive material. The effectiveness of the purification method was established by the isotopic dilution method." The same technique was used to show that only negligible decomposition of neutral or of carbonate solutions of the dithiophosphate took place over a period of a day; this was true whether or not galena was suspended in the solutions. Procedures: Radioactivity: The activity of a solution was obtained by counting with a small glass-jacketed, silvered Geiger counter, using a conventional scaling circuit.' Crystals and other solids were counted under a bell-shaped, mica-window counter tube. Sufficient counts were made so that the probable error of counting was of the order of +I-3 pct. The specific activity in terms of counts per minute

Jan 1, 1951

By E. J. Fasiska, H. Wagenblast

The dilalion of a ivon by interslilial carbon was measured by two independent techniques —dilatometric mesurements at 719 c and X-ray measurments of the urlil cell parameters a1 room temperature after quenching. The relative expansion per increase in cavbon content by both methods is (2.9 + 0.3) x 10-3 pev at. pcl C nrzd is temperature- independent within experimental error. This corresponds to (6.3 5 0.6) cu cm per g-atom for the partial gram-atomic volume of carbon in a iron-only slightly smaller than the atomic volume for iron in the same temperature vange. THE only previous quantitative study of the dilational effect of carbon dissolved in a iron was performed in 1934 by Burns1 which, at the time, generated some discussion of possible sources of experimental error.2 For a system of such widespread importance, we felt that a new investigation was merited. Both X-ray diffraction and dilation measurements were used to determine the expansion of the a iron lattice by dissolved carbon, avoiding as much as possible any previous experimental problems and deficiencies. The dilation method at solution temperature offers not only measurements which are free of residual strain but also, in conjunction with the room-temperature X-ray measurements, a method to detect any large temperature dependence of the partial gram-atomic volume of carbon. To insure that quenching strains did not affect the room-temperature X-ray measurements, wire specimens of constant carbon content but different diameter were examined for such an effect. SPECIMEN PREPARATION The material used for both experimental techniques was "Ferrovac E" iron received in the form of 19-mm-diam rod and stated as having the following impurity contents: C, Cr, Cu, Mn, P, and V, each in the range of 10 to 50 ppm; Co, O, Mo, S, and Si, 60 to 100 ppm; W, 200 ppm; Ni, 230 ppm; N, 4 ppm. The stock material was cold-swaged to 0.71 mm diam for the Debye-Scherrer X-ray camera specimens and portions were cold-rolled from 6.3 mm diam to 0.79- and 0.25-mm sheet for X-ray diffractometer and dilation measurements. The wires were annealed in wet hydrogen for 6 hr at 840°C and 15 hr at 720°C, and then quenched into 0°C water. A chemical analysis for carbon after this treatment gave 0.0046 + 0.0014 at. pct C. Three portions of these wires were subsequently held at 719°C in three different hydrogen + methane mixtures and then quenched, resulting in carbon concentrations of 0.0283, 0.0598, and 0.1067 at. pct C by chemical analysis. After carburizing, the wires were swaged to 0.48 mm and electroplated with silver to prevent carbon loss during subsequent heat treatment. The final heat treatment consisted of holding the wires at 72 1°C for 5 min in a He-2 pct H mixture followed by quenching into 0°C brine. The wires were held at room temperature for a few minutes to remove the silver plating using a phosphoric acid-hydrogen peroxide solution, and then stored in liquid nitrogen until the measurements were made. EXPERIMENTAL TECHNIQUE A) Dilation Measurement. The dilatometer consisted of a gas-atmosphere vertical tube furnace modified so that length changes of a ribbon-shaped specimen could be measured externally. This was done by installing a gas-tight mercury seal at the top of the furnace as shown schematically in Fig. 1. The specimen (21.0 cm long, 1.27 cm wide, and 0.25 mm thick) was suspended in the center of the furnace by 3-mm-diam quartz rods with the upper one passing through the cap of the mercury seal. Above the cap, the upper quartz rod was coupled to a lever exerting a load of about 25 g (-12 lb per sq in.) on the specimen and having a 10 times mechanical magnification. The vertical position of a marker at the other end of the beam was read with a traveling microscope with a precision of 0.01 mm. The temperature gradient of the furnace was meas-

Jan 1, 1968

By H. S. Thomsen, W. C. Duesterhoeft, R. E. Hartline

The attenuation and phase shift which formations produce in the electromagnetic field of an induction-type electrical well-logging instrument are great enough to substantially affect the response of the tool to formation conductivity under normal logging conditions. The application of the general equations for the propagation of an electromagnetic field in a conductive medium to the transmitter-receiver coil pair of an induction-logging-tool coil system gives an expression for the response which properly takes the propagation effects into account. A comparison of a calculation of this type with the response computed from the commonly used geometric factor concept, which does not include the propagation effect, leads to a factor Go representing the ratio of the responses computed by the two methods. The value of Go decreases not only with increasing formation conductivity, but also with increasing transmitter-to-detector coil spacing. For a single coil pair with a 40-in. spacing, the value of Go is 0.972 in a 20-millimho/m formation. The value is reduced to 0.915 and 0.812 for conductivities of 200 and 1,000 millimhos/m. As a result, the basic signal generated by the induction-logging coil system is not linearly related to formation conductivity as expected from the geometric factor concept. A uniform conductivity scale can be obtained only by adding a suitable nonlinear element to the recording system. The addition of auxiliary coils, having spacings less than the main-coil span, to achieve focusing results in an even greater departure of the basic signal from linearity with conductivity. The solution of the field equations near the interface between two formations of different conductivity gives the curve shape on an induction logging tool in crossing the interface. The addition of auxiliary coils to achieve focusing can add anomalous character to the curve shape in crossing the interface, which might be mistaken for lithological detail. Preliminary calculations which include the propagation effects in the determination of the conductivity correction for thin beds lead to correction factors which are substantially smaller than those obtained from geometric factor considerations. It is apparent that thin-bed corrections derived from geometric factor calculations are of doubtful validity. INTRODUCTION Although the electromagnetic induction type of electrical well log has been used for a number of years and has gained general acceptance in quantitative well-log interpretation, no technically complete investigation of the mechanism and of its operation has been published. The studies which have been presented are based upon the premise that neither the shape nor the intensity of the induction field of the logging tool is in any way affected by the electrical characteristics of the formation. From the inception of induction logging it has been recognized that under some formation conditions the performance of the instrument could be expected to be substantially different from that predicted by the simplified analysis. This anomalous performance has been called "skin effect" since it is the result of the same phenomenon which causes high-frequency alternating currents to flow only near the surface of metallic conductors. In the range of conductivities encountered in the sediments, the formation actually does not exhibit the sharply defined conducting "skin" shown by metals. The gradual modification of the distribution of the currents in the formation, produced by the induction field of the logging tool, is the direct result of phenomena associated with the propagation of the electromagnetic field through the formation. Therefore, it can be much more appropriately described as the "propagation effect". One purpose of this paper is to review the results of a solution of the equations for the performance of a single transmitter-receiver coil pair induction-logging system which takes into account the propagation effect, showing the manner in which the tool performance established by this analysis differs from that derived from the less comprehensive analyses currently in use. PRINCIPLE OF OPERATION OF THE INDUCTION LOGGING INSTRUMENT The high-frequency alternating current which is maintained at a constant value in the transmitter coil

J. Chipman (Massachusetts Institute of Technology, Cambridge, Mass.)—The fact that the experimental work has been applied to copper rather than iron and that the paper is presented to the Iron and Steel Division, I regard as rather significant. It shows the unity of metallurgy and the fallacy of trying to cut it up by metals. This result for the solubility of sulphur in molten copper correlates with Professor Schuhmann's finding that the published data on the other side of the copper-sulphur miscibility gap are also in error. I should like to ask the author to say a little bit more about the sulphur capacity of the slag. T. Rosenqvist (author's reply)-—I hope that Dr. Chipman will find the derivation of the expression for sulphur capacity more clearly explained in the printed version of the paper than in the oral discussion at the meeting. I feel that this quantity, which actually is the ratio of two activities, can be measured more easily than the individual activities. Even if the ratio CaO/CaS is chosen as the standard state, the expression can be used for any slag, even for slags completely free of lime, and it represents a way to put the desulphurizing power of all slag constituents into one bag. Some doubt has been expressed as to whether oxygen ions really exist in calcium oxide and in molten slags. From a thermodynamic view point that question is of minor importance. The term oxygen ion activity, or any activity for that matter, is defined rigorously by the equation: activity = exp u/RT, where u is the change in free energy connected with the transfer of one mol of ions from the standard state into the slag. Whatever happens to the ion in the slag is of no concern to the thermodynamicist. Regardless of whether the ion is "free" in the slag or not, or whether it is present in a very small amount, its activity can always be expressed, and for a thermo- dynamic calculation that is all we need. However, ionic activities will only be of some real value if they are simple functions of the slag composition, or can be measured easily. Concerning the real nature of the oxygen in the slag, my feeling is that the oxygen atom has a rather multiplex nature depending on how strongly it is tied by covalent forces or polarized by the other atoms or ions present. The oxides of iron, cobalt, and nickel differ from calcium oxide and blast furnace slags as to the amount of free electrons that can give rise to electronic conductivity. In slags we know that the conductivity is mostly ionic. The fact that reversible emf's can be obtained with oxygen electrodes in certain salt melts, indicates a significant amount of oxygen ions in these melts. But extended work, e.g. polaragraphic studies and measurements of transference number; are needed to obtain quantitative information about the real structure of the slags. D. E. Babcock (Republic Steel Corp., Youngstown, Ohio)—-With reference to the ion, it might be well to remember Dr. Moses Gomberg. All of his life he had no use for the ionization theory and he contributed greatly to the field of chemistry on the assumption there was no such thing as ions. I do not think we have to worry about whether the oxygen is ionic or not. I think one thing specifically should be brought to your attention and this I think is one of the important contributions of Dr. Rosenqvist. He pointed out what we know as oxygen potentials or what is described as oxygen potentials. I have used this concept for a long period of time and I want to state that if this concept is properly applied, it vitiates much of what we have in the literature, or makes our usual ideas regarding oxidation seem primitive. That one thing is more valuable than almost all the rest of the discussion as a fundamental basis on which to build a reasonable in-

Jan 1, 1952

By A. U. Seybolt

Experiments on both suljidation and hot corrosion have been carried out using ternary Ni-Cr-X alloys and commercial nickel-base superalloys. It has been shown that lhere are certain micro structural similarities between sulfidation using H,S/H2 mixtures and hot corvosion using Na,S04 in "crucible" experiments or hot corrosion tests using a flame tunnel apparatus. It is proposed that two key processes in the hot cowosion process are: 1) formation of chromium sulfides which deplete the alloy matrix, and 2) oxidation which can be accelerated because of chromium depletion. Other processes play an important role in the hot corrosion process, such as accelerated oxidation attack due to the presence of a N& S04 or N&SO,- NaC2 film which inhibits the formation of a protective external oxide scale. However, this latter phenomenon was not studied in this work. THE term hot corrosion has come to mean a combination of sulfur and oxygen attack on nickel-base or cobalt-base superalloys used as gas turbine blades or partitions. This attack can be initiated by any combustion process in air involving fuel containing sulfur if the sulfur level is sufficiently high. Levels of 5 pet SO2l for example, in the combusting gas have been observed to produce hot corrosion attack on nickel-base alloys. The important immediate corroding medium in aircraft gas turbines, however, has been established to be Na2S04, formed from reaction between NaCl from sea air and sulfur in the fuel.' The Na2S04 concentrates the sulfur sufficiently so that superalloy parts heated in the range 1400" to 1900°F can become sulfided by reducing the Na2S04 film which forms on the parts, and this reaction allows sulfur to gain access to the alloy. The sulfur thus released diffuses into the superalloy and forms primarily chromium sulfides of various compositions, but frequently Cr2S, has been identified. Figs. 1 and 2 show the morphology typical of hot corrosion of a cobalt-base superalloy X-45 (25 Cr, 10 Ni, 7.5 W, 0.25 C) and a nickel-base alloy, Udimet 500 (19 Cr, 19 Co, 4 Mo, 3 Ti, 2.9 Al, 0.07 C). Both alloy samples were subjected to the same hot corrosion conditions: 616 hr at 1600°F in a combusted 2.9 pct S distillate oil containing 125 ppm Na as NaC1. It will be noted at once that the nickel-base alloy is much more severely attacked. The depth of corrosion is much greater and the number and size of the sulfides present are much greater. Both alloys show a duplex oxide scale, which is more noticeable in the X-45 alloy. The external scale is probably a spinel in both cases, with CrQ3 lying next to the alloy. Many investigators1 -"ve been concerned with hot corrosion phenomena, but there has been a lack of unanimity of opinion as to the details of the hot corrosion mechanism. More than one definition of hot corrosion is possible. Some investigators con- sider as a measure of hot corrosion only the weight gain due to an accelerated oxidation. However another point of view is that the total depth of affected alloy may be used as a measure. While it is recognized that surface oxidation attack is greatly accelerated by the presence of salts like Na2S04 or Na2S04-NaC1 mixtures, sulfur penetration into the alloy also undoubtedly produces very important effects. The region of sul-fide penetration which is so prominent in Fig. 2 is one of the outstanding characteristics of hot corrosion. For highly stressed parts, such as first-stage buckets, this sulfur-penetrated region can be more significant in mechanical property loss than the reduction of cross section produced by the formation of the massive external scale. In addition, evidence will be presented that this sulfided region can be less resistant to subsequent oxidation attack. Since it appears that the sulfided region below the surface oxide scale plays an important role in the overall corrosion process,

Jan 1, 1969

By C. E. Kemp, A. B. Dyes

Results of a field research project on the thermal recovery of oil by movement of a combustion front are presented. This field test was conducted in the South Belridge field, This The war paitern was a Southsing1e 2.5-acre five-spot pattern in the 700-ft deep Tulare sand. Results of this project indicate the technical feasibility of producing high viscosity, low-grallity crude oils by thermal means. It was demonstrated that heavy oils could be moved rapidly over considerable distances by thermal drive, and that high percentage recoveries of oil could be effected. The bulk of the oil recovery came after arrival of the burning front at a production well. However, rapid corrosion and high-temperature failure of conventional pipe and equipment both below and above surface occurred. Stainless steels gave more satisfectory performance under these conditions than did conventional steels. Oxygen was seldom observed in gas from wells ahead of the burning front, and it appears that the front was continuous throughout the life of the experiment. INTRODUCTION During late 1953 and early 1954, industry-wide attention was focused on the thermal recovery of oil by movement of a combustion front through an oil zone. Two separate Oklahoma field tests were described in publications by Kuhn and Koch, and Grant and Szasz. To obtain basic information needed for making engineering and economic appraisals of commercial recovery of heavy oil by this thermal recovery method, a new field research projecta3,4 was started May 31, 1955. Twelve oil companies contributed to this project with General Petroleum Corp. as operator. Since combustion recovery experiments* with heavy oils in other areas had been conducted in small patterns,"' an important purpose of this test was to determine whether the process could be operated in a pattern of near-commercial well spacing. The south Belridge field, located 30 miles north of Taft, Calif., was chosen because its oil sands are similar to many heavy oil sands which typically have low primary recovery and leave a large amount of oil in place. such fields are not particularly suited for common secondary recovery methods. DESCRIPTION OF TEST SITE The test interval was a Tulare sand at a depth of about 700 ft. Nine test wells were drilled, cored and completed from Oct., 1955, through Jan., 1956. The initial pattern shown on the map, Fig. 1, consisted of an injection well (II)', four producing wells (1P-4P), and four observation wells (IT-4T). An old well (81) was also recompleted for use as an observation well. The area enclosed by the original four production wells was 2.5 acres; the distance from the injection well to each production well was 233 ft. During the combustion period, Well 5P was drilled to replace Well 2P, and Well 6P was drilled to replace Well 3P. The new wells were 25 ft beyond the old production wells along a diagonal from the injection well, and resulted in an increase in pattern area to 2.75 acres. Isopachs are shown in Fig. 1. The average thickness of the oil sand in the test pattern was 30 ft. Fig. 2 presents two cross sections through the test pattern. Oil sand is represented by white and claystone by gray. The dashed lines indicate layers of conglomerate grading to claystone. The sand thickness was fairly uniform throughout the test pattern except in the 1P area. The sand body dips gently from the northwest to southeast, about 3 degrees from horizontal. The sand is unconsoli-dated, and grades from very fine material to pebbles as large as 2 in. in diameter. Table 1 presents the average initial properties of the formation. The in-place porosity of 37 per cent was estimated from the porosity of coke-consolidated cores taken after combustion. Original core analyses data were adjusted to 37 per cent, giving pore volume saturations of 60 per cent oil, 37 per cent water, and 3 per

By Michael P. Robinson

It has been suggested that streaming potentials are not nomlally logged because the streaming potentials known to be generated across mud filter cakes are substantially cancelled by streaming potentials generated in contiguous shale beds. Laboratory measuments of streaming poterztials acrons hand specirnens of shale have been adduced as evidence to support the suggestion. The experiments show that streanzing potential depends on the resistivity of the equilibrating solution and on pressure drop. Meizsurenlents were made perpendicular to the bedding planes of the shales. The mechanism of streaming potential generation in shale which the suggested cancellation theory iniplies has been critically examined. Conclusions reached are tested by comparing self potentials logged before and immediately after changing the resistil~ity (but not the weight) of mud in a borehole. The resultant SP change is found to equal the theoretical electrochetnical potential change, i.e., cancellation of .strearning potential is independent of mud resistivity. Since streaming po/entials generateal across mud filter cakes depend critically upon mud reistivity, shale .strearlling potential in situ must he equally sensitive. Using plausible values for the pemleahility, porosiy arzd compressibility of shale, the van Ererdingerl-Hurst approach shows that the pressure drop between horehole pressure and formation pressure in a radially honzogeneous shale is unlikely to coincide with the zone of filtrate invasion. If, perchance. this coincidence occurs with one mud a chatzge 10 a second must result in non-coincidence, i.e., a .YP logged in the second mud should show a streaming potenlial. The dilernrna may he overcome by postrrlating: (I) perrneabilities of shale7 parallel to their beddine planes greatly exceed their cross-bedding permeabilities; (2) a resistivity-sensitive mud cake of low permeability forms either on the face of, or probnhly in the lanlinations of, many shales. These postulates are examined with reference to the pres.sure and SP histories of the McClosky and Aux Vases zones in the Illinois Basin and some srcpport for them found. It is concluded that a definitive solultion to the streaming potential problem does not exist. INTRODUCTION Recently Gondouin and Scala' published a contribution to the streaming potential problem which is notable particularly for the elegance of its experimental techniques. The conclusions drawn by Gondouin and Scala from their studies amplified those previously suggested by Schehck'. In essence they suggested that the streaming potentials known to develop across mud filter cakes" are not normally observed on SP logs because they are largely cancelled by streaming potentials generated in adjacent shales. This suggestion appears so reasonable that many workers in the logging field will be prone to accept it. We incline strongly to the view that the suggestion is basically correct. Nevertheless, it is important to realize that the cancellation hypothesis in the form put forward by Gondouin and Scala does not appear to meet all the known experimental facts. It is the principal intention of this note to call attention to what appears to be a flaw in the relevancy of the experimental evidence adduced by Gondouin and Scala and to suggest an alternative hypothesis to remove it. The consequences of the new hypothesis are compared with field data which, we believe, are unique. A second but not less important purpose of this note is to emphasize that even the new hypothesis does not fully explain all field observations. Hence we conclude that further study of the streaming potential problem is required before it can be said that a definitive solution to it has been attained. The bare experimental bones of laboratory data on streaming potentials appear to be as follows. 1. A streaming potential is developed across a mud filter cake. This streaming potential is virtually independent of cake thickness but depends on the mud resistivity. The relationship is of the form, 2. A streaming potential is developed across a shale. This potential is independent of shale thickness but depends on the resistivity of the permeating solution if the shale has been first equilibrated in this solution. The relationship is of the form, The relationship in (1) was determined first by Wyllie3, confirmed by Gondouin and Scala and qualita-

By I. B. Cadoff, J. C. Bilello, R. Rosenberg

Diffiusion of magnesium into the surface of LiF crystals to controlled depths and subsequent heat treatments provided a wide range of surface zone harahesses and structure, The bend strength of the LiF crystals was increased by as much as an or-dev of magnitude. Ductility was achieved when dislocation generation occurred in the diffusion zone or when dislocations penetrated to the surface from the intevior. A critical surface hardness of 130 to 140 kg per sq mm was found helozu which generation could take place in the diffusion zone and ahoue which the zone was impenetrable, This hardness was obtainable by several methods, among them the aging of quenched MgF2 -LiF solutions to produce MgF, precipitation. Maximum hardness was ohtained in quenched specimens with no visihle evidence of MgF,. Diffusion-zone formation followed a parabolic rate law and an activation energy of 20.9 kcal per mole was obtained for the process. RECENTLY, the properties of ionic crystals as related to surface condition have been receiving much attention, specifically the transitions between ductile and brittle behavior. Originally Joffe 1 showed that NaCl crystals could be made ductile by immersion in water and related this to the elimination of surface microcracks. Aerts and DeKeyser 2 and Gorur 3' have subsequently shown that ionic crystals are inherently ductile and are embrittled through contact with air. Machlin and Murray4 hypothesized that a layer of NaCIO3 produced by contact of ozone with NaCl induced embrittlement by acting as a barrier to outward dislocation flow. westwood,' Rosenberg and Cadoff,9 and Bilello and cadoff' have reported surface strengthening of LiF crystals by coating with a magnesium compound and then heat treating for adherence. The major effect of the coat was to inhibit dislocation-slip lines from reaching the specimen surface. westwoods showed microcrack formation and fracture to be caused by slip-band interactions at the surface. The material presented in this paper is an extension of the work reported earlier by Bilello, Rosenberg, and cadoff'6,7 and illustrates the wide range of surface properties and bulk behavior obtainable by use of heat-treated magnesium-diffused surface regions in LiF crystals. EXPERIMENTAL PROCEDURE The LiF single crystals were obtained from the Harshaw Chemical Co. Some batch to batch variation was observed; therefore all specimens for a given test series were cleaved from the same crystal. The typical dimension used was 1 by 0.1 by 0.40 in. Surface damage resulting from cleavage was removed by chemically polishing in a 2 pct NH4OH solution. Each group of specimens was given a vacuum anneal at 700°C for 4 hr to provide a base standard for measuring comparative effects of various surface treatments. To produce the reacted surface zone, the annealed specimens were immersed in a boiling suspension of MgF, in doubly distilled HzO, agitated slbwly for 30 sec, removed, and dried at room temperature. Uniform coatings of MgF, were deposited with a thickness of approximately 5 mils. It should be noted that this technique can be modified for use with crystals which are soluble in water by using boiling absolute alcohol as the dissolving medium. This was found effective for the coating of NaCl with MgF2. The diffused surface zone was obtained by annealing the coated samples at elevated temperatures in a vacuum of lob4 mm Hg. Penetration depth was controlled by varying the annealing time from 1/2 to 28 hr. After heat treatment, the samples were tested for bend strength and hardness. Load was applied by four-point bending in a hard-beam testing jig. Four-point rather than three-point bending was used to provide a wide area of constant stress and to minimize the effect of localized inhomogenities in the specimen. The deflection rate was 8 x min-' and the distance between knife edges was 1/4 in. Load-time curves were obtained from a chart recorder coupled to the machine and converted to resolved shear stress on the shear plane vs deflection, as plotted in the figures. The unstressed portions of the sample outside of the two outer knife edges were used for the microhardness studies. Microhardness measurements were made with a Bergsman tester attached to a Reichert metallograph. All hardness impres-

Jan 1, 1964