Search Documents

By T. Y. Yan

SUMMARY Laboratory high pressure column tests have shown that the presence of 1-20 ppm of aluminum ion effectively prevents permeability loss during uranium leaching with leachates containing sodium carbonate. If added after permeability loss has occurred, aluminum ion can restore the permeability to nearly its original value. No deleterious effect was observed on uranium leaching performance and the technique should be quite compatible with all field operations. INTRODUCTION The recovery of uranium values from underground deposits by in situ leaching or solution mining has become economically viable and competitive with conventional open pit or underground mining/milling systems (Merrit, 1971). In situ leaching processes are particularly suitable for small, low-grade deposits located in deep formations and dispersed in many thin layers. Many such ore bodies occur along a broad band of the Gulf Coastal Plain (Eargle et. al., 1971). The advantages of the in situ leaching processes have been reviewed (Anderson and Ritchi, 1968). In the in situ leaching process, a lixiviant containing the leaching chemicals is injected into the subterranean deposit and solubilizes uranium as it traverses the ore body. The pregnant lixiviant or leachate is produced from the production well and is then treated to recover the uranium. The resulting barren solution is made up with the leaching chemical to form lixiviant for re-injection. Upon completion of the leaching operation, the formation is contaminated with leaching chemicals and other species made soluble in the leaching operation and has to be treated to reduce the concentration of these contaminants in the ground water to levels acceptable to the regulatory agencies (Witlington and Taylor, 1978). Restoration is accomplished by injecting a restoration fluid, which could be the fresh water or water containing chemicals, into the formation. As it traverses the leached formation, the restoration fluid picks up the contaminants and is then produced at the production well. This produced water is either disposed or purified for recycle. In both phases of operation, formation permeability or well injectivity is one of the most important parameters which determines the viability of the in situ leaching process. Low formation permeability limits production rates, leading to uneconomical operations. The formation is said to be sensitive if there is a sharp loss of permeability on contact with water and other fluids. Many uranium bearing formations, for example, the Catahoula formation of the Texas Coastal Plain, contain significant amounts of clay minerals which are water sensitive. Serious permeability losses can occur when the pH and chemical composition of the lixiviant is significantly different from that of the formation water. Jones has investigated the influence of chemical composition of water on clay blocking of permeability (Jones, 1964) and Mungan studied permeability reduction through changes in pH and salinity of the water (Mungan, 1965). Various mechanisms of permeability damage have been proposed and reviewed (Jones, 1964; Mungan, 1965; Gray and Rex, 1966; and Veley, 1969). When large amounts of swelling clays are present, a significant fraction of the flow channels in the formation can be reduced due to swelling. However, in most cases, swelling need not be the main cause of permeability losses. Particle dispersion and migration or clay sliming can be more important causes for formation damage. Clay particles entrained in the moving fluids are carried downstream until they lodge in pore constrictions. As a result, microscopic filter cakes are formed by these obstructions, plugging the pores, effectively restricting fluid flow and reducing the formation permeability. Moore found that as little as 1-4 percent clays present in a fine grained sandstone could completely plug the formation if they are contacted by incompatible injected fluids (Moore, 1960). It has been found that injection of NaHC03/Na2CO3 lixiviant into formations with significant clay content often leads to loss of formation permeability and well injectivity. To alleviate this problem a change of the lixiviant composition to KHC03/K2CO3 has been proposed. At present, however, many in situ leaching operations employ NH4HC03/(NH4)2C03 mixtures as a source of carbonates. This approach has been successfully used in South Texas by Mobil, Intercontinental Energy, Wyoming Minerals and U.S. Steel, etc. The use of ammonium carbonates solutions, however, contaminates the formation and requires a time-consuming restoration operation. The other approach to reduce the permeability loss is to pretreat the sensitive formation with chemicals which prevent clay dispersion and migration. Such chemicals include hydroxy-aluminum (Reed, 1972 and Coppel et. al., 1973), hydrolyzable zirconium salts (Peters and Stout, 1977), hydrolyzable metal ions in general (Veley, 1969) and polyelectrolyte polymers (Anonymous). Still another approach, is to minimize the "shock" caused by sudden injection by gradually changing the chemical composition of the injected fluids from that of the formation water. THE APPROACH Since permeability loss can be an important factor limiting the efficiency and economic viability of the in situ leaching process, a study was initiated on

Jan 1, 1982

By W. Philippoff

ALTHOUGH Gaudin1 and more recently Sutherland2 have calculated the probability of collision of a falling mineral particle with a rising bubble, there is no published information concerning the details of the mechanism of attachment of a collector-coated particle to a bubble. During the past year the writer has developed a theory for the mechanism of attachment, which has been substantiated experimentally.' Funds for the investigation and for some of the equipment used have been supplied by the Mines Experiment Station of the University of Minnesota. Motion picture studies of the phenomena involved in the collision between mineral particles and bubbles, such as those of Spedden and Hannan,3 show that the contact can be completed within 0.3 millisec. Formulas developed for rigid bodies have hitherto been used' for the calculation of the motion of bubble and particle, but it is obvious that a bubble cannot be regarded as a rigid body. On the contrary, Spedden and Hannan's pictures show a great degree of deformation during the collision. The time of attachment was calculated as the time the particle drifted past the bubble. Time of Collision The theory presented in this paper enables calculation of the time of collision; using the concept that the bubble, or more generally, a liquid-air interface, acts as an elastic body. The elasticity, defined as the restoring force on a mechanical deformation, is caused by, the surface tension and is the result of the principle of the minimum of free surface energy. It is well known that an elasticity together with a mass determines a frequency of vibration. The vibrations of jets and drops caused by the elasticity of the interface are known to comply exactly with the classical theory of capillarity.5 However, the vibrations of isolated bubbles, as distinct from foams, have not been investigated previously. The following equation, presented elsewhere,' has been deduced for these frequencies: [3fB = 9.20•'./V•Vn- (n-1) • (n+2) /8[1]] in which fB is the frequency of a harmonic of the bubble in cycles per second, V the volume of the bubble in cc, n a number determining the order of the harmonic, and n = 2 the basic vibration. The first (basic) harmonic describes a change of the spherical bubble to an ellipsoidal bubble. The higher harmonics are more complicated, for the circumference of the bubble is divided approximately into as many parts' as the order of the harmonic. As an example, Spedden and Hannan's published motion picture of, a vibrating bubble corresponds to the sixth harmonic. Eq 1 shows that only the first and third harmonics are simple multiples (1 and 3), all the others being irrational fractions of the basic frequency. This means that the shape of the vibration can change with time and is in general unsymmetric in respect to the time axis. Such conditions prevail when there is a distributed elasticity or mass, as in the case of vibrating membranes or rods. The constant 9.20 is valid for water at room temperature, but a general solution involving the physical constants of the liquid has not been found. The case of the floating particle is much easier to treat I than that of the bubble. It can be assumed that the elasticity is caused exclusively by the interface and that the mass is concentrated in the particle together with some adhering water. The following expression for the frequency of a system, of one degree of freedom can be applied: [1E/m[2] fP = 27] Here f, is the frequency of the particle vibration in cycles per second, E the elasticity in dynes per cm, and m the mass in grams. The classical theory of impact phenomena gives the time of collision during the striking of a spring (in this case the surface of the bubble) by a mass, as: [t~ = 2/f = 7r\/m/E[3]] It is now possible to develop an expression for the elasticity of a floating cylindrical particle. The force equilibrium of a cylinder floating end on at the air-liquid interface is given by the well-known equation (Poisson' 1831) [aP = 4 D2.pL•g•h +7rD•y sin a[4]] which accounts for the buoyancy and the action of the surface tension where P is the force acting on the particle in dynes (weight-buoyancy), D the diameter of the cylinder in cm, pL the density of the liquid in grams per cc, g the acceleration of gravity = 981 cm per sec2, h the depression of the cylinder below the surface of-the liquid in cm, y the surface tension in dynes per cm and a the supporting angle' or the one required to insure equilibrium, a being smaller than the contact angle ?. Although demonstrated by Poisson, it has not

Jan 1, 1952

By Francis Collins, E. R. Brownscombe

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

Jan 1, 1949

By L. L. Handy

Solvent floods using slugs of solvent have been found to show continuity in behavior from the vapor pressure of the solvent to the critical pressure for the two-component driving gas-solvent system. In the pressure region between the solvent vapor pressure and the critical pressure for the gas-solvent system, the gar and solvent are only partially miscible. Although complete miscibility cannot be obtained at these pressures, complete oil recovery is possible in principle. In two-phase solvent floods the solvent is propagated tllrough the reservoir, primarily, in the vapor phase. The carrier gas requirements constitute a significant factor in the economics of the process. A qualitative theory is proposed for estimating the amount of dry gas required to move the solvent through the reservoir. The theory shows that for two-phase solvent floods the total gas needed is a minimum at the vapor pressure of the solvent and at the critical pressure for the gas-solvent system, and is a maximum at some intermediate pressure. The predictions of the theory are supported by experimental studies using methane, butane and decane or methane, propane and decane in a natural sandstone core. INTRODUCTION Previously, solvent slug processes have been found effective for oil recovery in two pressure ranges. First, conventional miscible displacements are possible at pressures greater than the critical pressure for the gas-solvent system. Second, Jenks, et al,' have shown that, at pressures slightly in excess of the vapor pressure of the solvent, a solvent slug can be propagated through a reservoir by a gas essentially insoluble in the liquid solvent. The solvent bank displaces the oil ahead of it. Both of these processes, at least ideally, are capable of recovering all of the oil in the swept regions. Slug processes for which the gas and solvent are partially miscible have not been considered; that is, those systems for which the solvent and driving gas form two equilibrium phases in which the vapor phase contains a significant amount of solvent and the liquid phase an appreciable amount of the driving gas. Welge and Johnson' have shown that the gas needed to movc a solvent slug through the reservoir increases with increasing pressure above the vapor pressure of the solvent. It will be shown that solvent slug processes can, theoretically, recover all of the oil at any pressure greater than the vapor pressure of the solvent. But the amount of gas required to move the solvent through the reservoir depends very much on the pressure and temperature. In the present study a maximum in the gas requirements was both predicted theoretically and observed experimentally. This result has not been reported previously, and would not have been predicted from the Welge and Johnson model. The gas requirements are a minimum at the pressures corresponding to the vapor pressure of the solvent and again at the critical pressure for the gas-solvent system, and are a maximum at some intermediate pressure. AN APPROXIMATE THEORY FOR TWO-PHASE SOLVENT FLOODING The differences and similarities between conventional solvent floods and two-phase solvent floods are best understood by referring to concepts developed for miscible displacement in which miscibility is generated in the reservoir. In Fig. 1(A), a ternary diagram is shown for a hypothetical gas-solvent-oil system. To be rigorous the three components should each consist of a single molecular species. The pressure for Fig. 1(A) is greater than the critical pressure for the binary gas-solvent system at the specified temperature. Diagrams of this type are the ones most frequently referred to in discussions of enriched-gas drive and miscible displacement. A limiting tie line is shown tangent to the two-phase envelope and intersecting the gas-solvent line at Point A. To obtain generated miscibility with this type system, others have shown that, for an oil of Composition D, a mixture of gas and solvent must be injected which is richer in solvent than that composition indicated by A. An oil repeatedly contacted with a gas phase richer than A changes toward a composition which would be at equilibrium with the injected mixture, that is, a composition lying on a tie line which passes through the injected-gas composition. Since no such tie line exists, the oil is enriched to the point at which it becomes directly miscible with the injected mixture. At pressures lower than the critical pressure for the gas-solvent system, other types of phase diagrams are observed. The ones of interest in this paper are for pressures greater than the vapor pressure of the solvent, but less than the critical pressure of the gas-solvent system. Such a ternary diagram is shown on Fig. 1(B). In this case, two-phase behavior is observed not only for gas-oil mixtures, but also for certain compositions in the gas-solvent system. If a gas of Composition A (a dew-point vapor) is injected, once again the original oil is enriched by successive

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

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

Jan 1, 1965

By E. R. Brownscombe, Francis Collins

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

Jan 1, 1949

By Che-Yu Li, J. L. Gregg, W. F. Brehm

Oriented, oxygen-doped niobium bicrystals were tested in liquid lithium. The grain boundaries were attacked preferentially. The depth of the penetrated zone varies as (time)2. The penetration was aniso-tropic, had a high activation energy, and increased with the increased oxygen doping level. A possible model was proposed to account for the experimental observations. 1 HE grain boundary penetration of a metallic system by liquid metal has been studied by several investigators. Their results are summarized by Bishop.' Most of these works show that the penetration by liquid metal corresponds to the phenomenon of liquid metal wetting. In the case of a grain boundary, wetting will occur when twice the solid-liquid interfacial tension is smaller than the grain boundary tension resulting in the replacement of the grain boundary by two new solid-liquid interfaces. Other possibilities exist; for example, the atoms of the liquid metal may diffuse into the grain boundary region due to chemical potential gradient. The gradient can be produced by impurity segregation or simply be due to the increase in solubility in the grain boundary region. The penetrated grain boundary in these cases may remain solid at the test temperature. The Nb-Li system has been of considerable interest because of its possible technological applications. For fundamental interest it provides a possibility of studying the grain boundary penetration process which is not controlled by the wetting mechanism. The pure niobium is not attacked by the liquid lithium, but if niobium containing more than 300 to 500 ppm oxygen by weight is exposed to liquid lithium, corrosion occurs at the solid-liquid interface and preferentially at grain boundaries. Previous investigators2-' have proposed that this preferential corrosion at grain boundaries is caused by oxygen segregation there, with subsequent inward diffusion of lithium to form a Li-Nb-0 compound. These investigators also found that the corrosion could be retarded by adding 1 pct Zr to the niobium to precipitate the oxygen as ZrO2 upon proper heat treatment. However, there are no quantitative data on the kinetics of the grain boundary penetration process to test the validity of the proposed corrosion mechanism. In this work an investigation of this penetration process in oriented bicrystals was made as a function of the oxygen doping level in the bulk niobium and the grain boundary orientation. A possible model for the penetration process based on the experimental results was proposed. EXPERIMENTS Oriented niobium bicrystals were grown by arc-zone melting oriented single-crystal seeds.7 These bicrystals contained simple tilt boundary. The [001] directions in the two grains were tilted about a common [110]. The bicrystals were 31/2 in. long and 5 by 4 in. in cross section with the straight, symmetric, planar grain boundary longitudinally bisecting the crystal rod. The bicrystals were doped with oxygen by anodically depositing a layer of Nb2O on the surface in a 70 pct HNO solution at 100 v, using a stainless-steel cathode. The specimens were homogenized by annealing in evacuated quartz tubes at 127 5°C. Oxygen content of the niobium was measured from microhardness values, after DiStefano and Litmman.' Supplementary checks were made with vacuum-fusion analysis.7 Individual test specimens cut from the doped bi-crystal rods, about by by % in. in size, were tested inside double jacket sealed capsules. The inner jacket was niobium, the outer was stainless steel. The niobium inner jacket eliminated the problem of dissimilar-metal mass transfer.' The lithium (99.8 pct pure, obtained from Lithium Corp. of America) was handled only in a purified argon atmosphere in a Blickman stainless-steel glove box. After introduction of lithium, the capsules were sealed by welding. Further detailed experimental procedures are given in Ref. 7. The capsules were heat-treated in vertical Marshall resistance furnaces. Temperatures were controlled to When heating above 1100°C, it was necessary to seal the furnace work tube and flow argon through to prevent failure of the stainless-steel outer jacket of the capsule. Tests were made on 6" 2", 16" 2, and 33" i2" bicrystals at oxygen levels up to 2600 ppm by weight in the 6' and 16" crystals and with 1300 ppm oxygen in the 33' crystals. The oxygen levels were controlled to 100 ppm. Most of the quantitative data were obtained from 16" bicrystals between 800" and 1050°C. The capsules were quenched into water after the test and cut open with a water-cooled abrasive wheel. The capsules were then submerged in water, which dissolved the lithium and freed the specimen. Measurement of the depth of the penetrated zone in the grain boundary was done either on metallographically prepared surfaces or directly on the grain boundary plane after the specimen was fractured in tension in the grain boundary plane. The depth of penetration measured by both methods agreed well. Further details describing these techniques have been reported elsewhere.'p7

Jan 1, 1969

By W. F. Chiao

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

Jan 1, 1969

By D. R. Parrish, T. M. Geffen, R. A. Morse, G. W. Haynes

Flow tests on small core plugs have indicated that a large amount of gas is trapped and not recovered by water flooding a gas sand. Instead of I to 15 per cent pore space, as is usually assumed, the residual gas saturation is 15 to 50 per cent pore space, and is thus of the same magnitude as residual oil after water flooding oil sands. A thorough investigation was made to ascertain that large amounts of residual gas actually remain in reservoirs after a water flood and that this condition is not merely a laboratory phenomenon. In field experiments, the amount of gas left in a watered-out gas sand was measured by use of a pressure core barrel and the residual gas saturation of two watered-out gas sands was determined by electric log evaluation. In the laboratory, an investigation was made of factors that could possibly cause the value of residual gas saturation as measured on small core plugs to differ from that in the reservoir, and the effect of these factors on the amount of residual gas saturation was studied. The factors studied include flooding rate, static pressure, temperature, sample size and saturation conditions before flooding All evidence established that a relatively high gas saturation is trapped in water flooded gas sands and that this residual gas saturation can be measured in the laboratory by tests on small core plugs. INTRODUCTION There. has been general agreement among engineers that very high recovery of gas could be obtained from natural reservoirs By water displacement. Gas recoveries of 80 to 95 per cent of the original gas in place have become the normal expectation in water drive fields. The assumption of high recovery has been based on: 1. low density and viscosity of gas compared with water; 2. the erroneous assumption that the flow relationship in a gas-liquid system where gas is the displaced phase will be the same as when it is the displacing phase. It has long been recognized that gas can flow at very low gas saturations (in the range of 1 to 15 per cent pore space) in systems where liquid is being displaced by gas. By assuming the reversibility of this process, the conclusion was reached that the residual gas saturation following water flooding of a gas reservoir would be the qame (1 to 15 per cent) as that at which gas first flowed continuously as a displacing phase. Recent laboratory relative permeability Studies have demonstrated that the flow characteristics are very different in gas-liquid systems, depending on whether gas is displacing or being displaced by, a liquid. Also, it has been shown that there is no difference between the flow characteristics of oil and water or gas and water in water wet porous rocks. The residual gas saturation that can be expected following water flooding of a gas reservoir then would be in the same range as the residual oil saturation normally expected after water flooding an oil reservoir! i.e., in the range of 15 to 50 per cent pore space, depending on the rock characteristics. Obviously, such a difference in residual gas saturation means very important differences in recoverable gas reserves from water drive reservoirs. For example, if the original gas saturation in a field were 70 per cent, and the residual gas following flooding were 35 per cent, only half of the gas in place could be recovered by complete water drive, compared to the previously expected 80 to 95 per cent. This is a situation in which complete pressure maintenance could result in very greatly reduced recovery, since straight pressure depletion recoveries from gas reservoirs can approach 80 to 90 per cent. Such a change in thinking must be based on more complete information than a series of small core tests. Hence the work reported herein was undertaken with the objective of determining whether or not residual gas saturations indicated from small core relative permeability tests at atmospheric pressure and room temperature are representative of the residual gas saturations which could be expected after water flood of natural reservoirs. A study was made both through laboratory and fields tests to determine any differences in residual saturation which might be occasioned by differences in pressure, temperature, rate of flooding, and original saturation conditions. Four separate types of laboratory experiments and two field mesurements were made in this investigation. The apparatus and testing methods of each will he discussed individually. LABORATORY EXPERIMENTS, APPARATUS AND PROCEDURE Relative Permeability Tests Steady state flow experiments Here conducted using the Penn State type apparatus. The equipment used has beer (described in a previous publication? "Irreducible" water satu-ration was established in the core' by imposing ,a capillary pressure of 45 psi before simultaneous flow of air and water

Jan 1, 1952

By D. E. Coates, G. R. Purdy, S. V. Subramanian

The morphological stability of the planar solid-liquid interface in dilute ternary alloys, undergoing steady-state unidirectional solidification, is analyzed in terms of both the constitutional supercooling principle and the perturbation methods recently developed by Mullins and Sekerka. First, various steady-state solutions for the two solute distributions ahead of a planar interface are examined. The nature of the solutions depends on the size and concentration dependence of the off-diagonal diffusion coefficients. W~thin the framework of the constitutional supercooling principle, a cumulative contribution to instability frorn the two solutes is found to exist in the absence of diffusional interaction. It is shown that the latter can produce a further enhancement of instability or can have a stabilizing influence, depending on the form of the liquidus surface and on the sign of the solute-solute interaction. A perturbation analysis, which ignores diffusional interaction, verifies the cumulative influence of lhe solute fields and demonstrates that the Mullins-Sekerka stability criterion for binary systems (with capillarity accounted for) can be readily extended for application to ternary systems. SOME time ago, Tiller et al.' calculated the solute concentration distribution ahead of the planar solid-liquid interface of binary alloys undergoing steady-state unidirectional solidification. An earlier qualitative proposal that the transition from planar to nonplanar growth morphologies is associated solely with the onset of constitutional supercooling in the liquid layer ahead of the moving interface2 was used in conjunction with this calculation to put the now well-known constitutional supercooling (C-S) stability criterion into quantitative terms. Mullins and Sekerka,3 in a recent and very elegant analysis, established a more complete criterion (hereafter referred to as the M-S criterion). Interfacial stability was investigated by determining the time derivative of the amplitude of a sinusoidal perturbation of infinitesimal amplitude which had been introduced into the originally planar shape of the moving interface. Of particular importance is the fact that capillarity was included in the boundary conditions of their calculation. The purpose of the present paper is to extend all of this earlier work on dilute binary systems for application to dilute ternary alloy solidification. The analysis is divided into three sections. In the first the two solute distributions ahead of a moving planar interface are considered. Mathematical solutions are de- termined for situations in which: a) diffusional interaction is negligible, 6) diffusional interaction must be considered but circumstances permit use of constant diffusion coefficients, and c) the concentration dependence of off-diagonal diffusion coefficients can be described by first-order dilute solution approximations. In the next section, a stability criterion analogous to the C-S criterion is developed and the influence of diffusional interaction on interface stability is analyzed. Finally, the perturbation formalism of Mullins and Sekerka, with capillarity included in the boundary conditions, is extended for analysis of ternary systems in which diffusional interaction is negligible. The study of interface stability in binary systems usually commences with the assumption that the equilibrium distribution coefficient and the slope of the liquidus line are constant at values corresponding to infinite dilution. Similar assumptions have not been introduced into the present treatment; that is, we do not assume planar solidus and liquidus surfaces joined by tie lines which yield constant distribution coefficients. The latter involves the assumption of no ther-modynamic interaction between solute species in both the solid and liquid. We consider a ternary phase diagram for which the solidus and liquidus surfaces are, in general, nonplanar and of course pass through the corresponding binary solidus and liquidus lines. These lines are not assumed to have constant slope. In the dilute regions we are concerned with, the following assumptions are made: i) The solidus and liquidus surfaces are of a form such that both the solidus and liquidus temperatures are monotonically varying functions of each solute concentration. ii) The tie lines are such that the equilibrium distribution coefficient of a given solute is greater than unity for every point on the solidus (or liquidus) surface or it is less than unity for every point. STEADY-STATE SOLUTE DISTRIBUTIONS IN THE LIQUID As will be demonstrated in the next section, a knowledge of the steady-state solute profiles is not a necessary prerequisite for the formulation of a ternary C-S stability criterion. However, in that details, such as the complete description of the equilibrium liquidus temperature profile, require an evaluation of the solute distributions, the overall treatment is enhanced if these distributions are determined. Consider a ternary system (solvent plus solutes 1 and 2) for which a planar solid-liquid interface is in unidirectional motion at constant velocity V. At this stage it is unnecessary to limit ourselves to dilute solutions. For a stationary frame of reference the generalized forms of Fick's equations are:

Jan 1, 1969

By G. Borelius

ABOUT the turn of the century, Gibbs' thermo-dynamic theory of heterogeneous equilibrium, on the one hand, and the experimental methods of thermal and microscopic analysis, on the other, gave to the physical metallurgist his first scientific tool, the equilibrium diagram. The classical equilibrium diagram of a binary alloy system shows the boundaries between ranges of homogeneous and heterogeneous equilibrium in their dependence of concentration and temperature. A homogeneous solid sohtion which on cooling passes such a boundary is assumed to precipitate, forming a mixture of two phases with different concentrations. The equilibriunl diagram and the equilibrium theory, however, give no information about the time scheme of the process or the intermediate states passed during precipitation. For this reason it satisfies neither the practical need of the metallurgist nor the curiosity of the physicist. As a matter of fact, in the heat treatment of alloys for technical use the objective very seldom is the equilibrium state. Thus good mechanical properties of construction material are connected, for the most part, with some intermediate state. As these intermediate states are thermodynamically unstable, there is, from a theoretical point of view, always to be expected a decay of the good properties with time; and it is a matter also of practical interest to know whether this natural life time of a material is of the order of, say, ten or thousands of years. Thus, for many reasons, there is a current demand to complete our knowledge of equilibrium through knowledge of the kinetics of the precipitation phenomena. From the point of view of the physicist, the most interesting question in this case is whether there are any general laws governing the kinetics. According to a generally accepted view, precipitation is ruled by two more or less independent phenomena, the formation of nuclei of a new phase and the growth of these nuclei. It is also commonly accepted that there is a tendency for the velocity of growth to increase with increasing temperature because of the increasing mobility of the atoms. There is also a tendency for the velocity of growth to decrease in the neighborhood of the two-phase boundaries. So far, however, very little is known quantitatively about this fundamental phenomenon in the case of solid metallic systems. In our laboratory attention has been directed especially toward the nucleation phenomena, and a series of measurements have been carried out with the guidance of a work- ing hypothesis (based on experiences from previous work on order-disorder transformations in alloys) about the influence on the nucleation of thermo-dynamic potential barriers. However, before discussing the experiments, the theoretical ideas will be considered. In a binary solid solution the arrangement of atoms on the lattice points approaches with increasing temperature a state of full randomness, as illustrated by the ball model of Fig. 1, that might represent a [111] plane of a face-centered alloy with 30 pct "black" and 70 pct "white" atoms. In reality the atoms are changing places continually with their neighbors so that the picture should rightly have been a moving one. On account of this thermal motion the concentration of black atoms within a certain group of, say, a hundred or a thousand lattice points fluctuates with time around the bulk concentration of 30 pct in a manner governed by statistical laws. With decreasing temperature two independent changes in this state grow more and more important. First, the mobility of the atoms decreases, and second, the forces between the atoms will have an increased influence on the fluctuations. In alloys with a tendency for precipitation, which are the concern of this lecture, the distribution function of concentration fluctuations will broaden, so that the relative probability of great local variations from the bulk concentration increases. Fig. 2 gives an example of such a fluctuation. When the alloy is supercooled below the solubility limit into the range of two-phase equilibrium, the fluctuations will now and then at some point give rise to a state that resembles the equilibrium state and thus will form a stable nucleus that is capable of growing by diffusion processes. In discussions with colleagues and in the literature, I have often encountered the idea that three or four atoms of the dissolved metal could form a nucleus of the new phase. A look at the ball model might be enough to indicate that this cannot be true. If it were true, there should be nothing but nuclei, whereas we know from experiment that nucleation must be a rather rare occurrence. In fact we have, as will be mentioned later, certain reasons to believe that the nuclei are formed by fluctuations containing some hundreds of atoms, which should be the order of the number of black balls in the fluctuating group of the figure, if it were extended into three dimensions. As a working hypothesis we have assumed that the fluctuations producing nuclei, though large and rare, still are ruled by the distribution laws of fluctuations of the supercooled solid solution in its initial state. Thus the probability of nucleation will be connected to the thermodynamic properties of the solid

Jan 1, 1952

By A. N. Holden

A DISLOCATION mechanism has been described by Cottrell' by which metals can yield locally, I. form Liiders bands, giving rise to a characteristic stress-strain curve with a sharp yield point and appreciable strain at constant or decreasing stress. It is undoubtedly the best mechanism that has been suggested to date." In its present development, however, the dislocation mechanism provides a more satisfying explanation for the sharp yield point than for the extensive localized flow occurring at the lower yield stress. The primary objective in this paper is to extend the dislocation mechanism to account for localized cataclysmic flow by a dislocation collision process and to give experimental evidence to support such a process. Only the yielding of iron containing carbon -will be discussed, although other metal-solute systems are known to behave similarly. Cottrell Mechanism In brief, Cottrell explains the yield point in the following way: The dislocations in iron which must propagate to produce slip usually lie at the center of local concentrations of carbon atoms, since segregation about these dislocatlons relieves some of the local stress resulting from them. A dislocation surrounded by a "cloud" of carbon atoms is thus anchored, and a higher stress is required to set it in motion than to move a free dislocation. Considering all available dislocatlons to be anchored in this fashion, the iron exhibits a yield point when the first dialocations break free and move through the lattice causing slip. This first breaking away of a dislocation enables other dislocations to break loose by "interaction" and the process becomes a cataclysm producing local deformation or Luders bands. The yield point in the stress-strain diagram for iron is absent in freshly deformed material, but returns gradually with time; the phenomenon is one aspect of what is called strain aging. The rate at which the yield point returns following straining depends on the temperature of aging. According to Cottrell the rate of return of the yield point in strained iron is limited by the rate of diffusion of carbon at the aging temperature, the mechanism is onr: of reforming the solute atmospheres around carbon-free dislocations that had stopped moving coincident with the removal of stress. If the specimen is retested immediately after straining and unloading, carbon will not have had time to diffuse to, and re-anchor, dislocations and the yield point will not occur. The carbon diffusion limitation for the rate of strain aging apparently applies if the criterion for strain aging is either the change in hardness" or the change in electrical resistance" of the strained speci- men with aging time. The possibility exists, however, that the yield point actually returns to strained iron at some rate other than that deduced from hardness or electrical resistance data. Therefore, as a preliminary experiment, the rate of yield point return in a rimmed sheet steel strained 6 pct in tension was measured at 27°, 77°, and 100°C. A plot of yield-point elongation for each of these temperatures against aging time appears in Fig. 1. The aging process is described by curves which rise to a plateau value of elongation that seems independent of temperature, but at a rate that depends on temperature. Very long times lead to a further rise in the yield-point elongation above the plateau value. However, if the later increase in yield-point elongation is ignored and the log of the time to reach half the plateau value of elongation is plotted against 1/T, a straight line results for which an activation energy of about 25 kcal pel- mol may be assigned. Within the accuracy of this sort of experiment this is approximately the activation energy for the diffusion of carbon in iron (20 kcal per mol), and the carbon diffusion limitation suggested for the yield-point return on strain aging is valid. The Cottrell mechanism thus explains in a qualitative manner the occurrence of a yield point in iron and its return with strain aging. It fails, however, to explain some of the other experimental observations that have been made of the yielding behavior of iron. For example, it is known that the yield point in iron becomes less pronounced with increasing grain size. Annealed single crystals of iron have very small yield-point elongations .if indeed they have any,' compared to a polycrystalline steel. If the only requirement for a yield point is that the dislocations in the lattice of the annealed. material be anchored by carbon atoms, the difference in the behavior of single crystals and polycrystals is not explained. That a dislocation mechanism may be entirely consistent with little or no yield point in an annealed single crystal will become apparent later when dislocation interaction is discussed. Strain aging produces a definite yield point even in single crystals. This accentuation of the yield-point phenomenon in single crystals after strain

Jan 1, 1953

By R. C. McFarlane, T. D. Mueller, J. E. Walstrom

In evaluating uncertainty, experiments are usually performed repeatedly and then conclusions are drawn from the distribution of results. With the advent of high-speed electronic computers, it is possible to perform experiments using mathematical models constructed to simulate complex experiments or operations. Statistical methods are then applied to the results of the simulated experiments. This procedure forms the busis of this paper. Demonstrated is the need for properly accounting for uncertainty in petroleum engineering problems. How uncertainty affects solutions is evaluated in three example illustrations. The method used to evaluate uncertainty in petroleum engineering studies is the Monte Carlo simulation procedure.'-" INTRODUCTION The solution to most technical problems may be derived from interrelationships among several quantities called variables or parameters. There may be only a few variables or several hundred. Interrelationships among parameters may be explicit or implicit, well established or only approximate. Some variables that fully or partially depend on the magnitude of others are called dependent variables. Input variables for most practical problems are not precisely known; there is usually an uncertainty in their value. The degree of uncertainty may vary from one variable to another. Variables that are known accurately are called determinates.' For instance, the gravity of crude obtained from a particular pool may be known precisely, and therefore is a determinate. The degree of precision with which a quantity can be determined increases as data describing the pool are accumulated during the development of the field and the producing life of the pool. The uncertainty of a parameter may result from difficulty in directly and accurately measuring the quantity. This is particularly true of the physical reservoir parameters which, at best, can only be sampled at various points, and which are subject to errors caused by presence of the borehole and borehole fluid or by changes that occur during the transfer of rock and its fluids to laboratory temperature and pressure conditions. Uncertainty may also result in attempting to predict future parameter values. This type of uncertainty is particularly evident in investment analyses involving future costs, prices, sales volumes and product demand. Uncertainty in the solution to investment problems is often called risk, and its study is called risk analysis.' Uncertainty also enters into biological and sociological analyses in which indeterminate factors are often important due to limited control of the experimental material. It is customary, in evaluating uncertainty, to perform repeated experiments and to draw conclusions from the distribution of the results of these experiments. With the advent of the high-speed electronic computer, it is possible to construct mathematical models which simulate complex experiments or operations and to perform the experiments repeatedly, utilizing the models. Statistical methods are then applied to the results of the simulated experiments This method forms the basis of the investigation reported here. PROBABILITY DISTRIBUTIONS FOR VARIABLES The uncertainty in the value of a variable may be indicated by a probabilistic description accomplished by expressing the quantity by a probability distribution. Many recognized probability distributions can be used to describe physical quantities. Recent studies used various types of distributions to describe core analysis data.',' However, for the examples in this paper, the uniform and triangular distributions are believed to reasonably approximate the data used (Fig. 1). The uniform distribution confines the variable between an upper and a lower limit. The variable may lie anywhere between the two limits. This distribution is used when no one range of values for a variable is more probable than any other, but information or intuitive reasoning indicates the variable will lie somewhere between the chosen limits. The triangular distribution is used for a variable when more data are available to indicate a central tendency of distribution. This allows postulating a "most likely" value to the distribution and upper and lower limits. In this case, as for the uniform distribution, the variable is not expected to assume a value less than the lower limit or greater than the upper limit. However, with improved quality of data it can be postulated that the variable will tend to assume a value close to the most likely value, and that there will be a decreasing probability for values away from the most likely value. The area under either of these probability distributions is equal to unity since it is assumed that there is a 100 percent probability that the variable will lie somewhere under the curve. An ordinate erected at any particular value of the variable divides the area under the curve into two parts: the area to the left of the ordinate represents the probability that the value of the variable will be equal to or less than the value of the variable at the position of the ordinate, and vice versa. The probability is zero that the variable will have any specific deterministic value. If two ordinates are drawn for any two values of the variable, the probability that the variables will have a value lying between these ordinates is equal to the area under the curve lying between the ordinates.

By R. E. Cech, T. D. Tiemann

It has been shown that hydrogen may be made to serve as a rapid and eflicient reducing agent for Cu, Ni, Co, and Fe sulfides if a scavenging agent for hydrogen sulfide is intimately mixed with the sulfide particles being reduced. Accelerated reduction kinetics are demonstrated for nickel sulfide. Copper, nickel, and cobalt sulfides, when treated at certain temperatures in a combined reducing agent-scavenging agent system, are converted to voluminous masses of fibrous metal product. Studies have been carried out to determine the conditions which lead, on the one hand, to irregular poly crystalline fibers and, on the other, to long single crystal filaments a few microns in diameter. A mechanism is proposed to account for the formation of single crystal filuments. The sulfide minerals of Cu, Ni, Co, and Fe are an important source of these metals yet there has been comparatively little scientific effort devoted towards understanding reduction mechanisms of these minerals. This may be, in part, due to the fact that the most convenient reducing agents for carrying out such studies, viz., hydrogen and carbon, do not react appreciably with sulfides. We have found that the reaction of hydrogen with metal sulfides can be markedly accelerated by placing a scavenging agent for hydrogen sulfide in close proximity to the metal sulfide. A brief series of experiments demonstrating relative reduction rates is reported in this paper to illustrate the effect. With the reduction process thus accelerated we have observed an unusual type of reduction behavior on some of the sulfides investigated. Under certain conditions the metallic product of the reduction reaction takes the form of filaments growing outward from the sulfide particles. The present paper deals largely with efforts to classify the various types of growth forms observed. This study has shown that filamentary growths from sulfides take a much greater variety of forms than has heretofore been reported by Ercker,1 Hardy,2 and Nabarro and Jackson3 in their reviews of metallic growths from copper and silver sulfides. THERMODYNAMIC CONSIDERATIONS The thermodynamics for hydrogen reduction of metal sulfides is quite unfavorable. For the sulfides considered here equilibrium constants typically range from 10-3 to 10-5. These low equilibrium constants impose severe kinetic limitations on reduction since hydrogen sulfide must be transported out of the system at concentrations of only a few hundred ppm. Unless extremely high gas flow rates are employed the atmosphere surrounding any sulfide particle will always be essentially in equilibrium with the sulfide. If, however, one places an efficient scavenging agent for hydrogen sulfide in close proximity to the metal sulfide particles the concentration of H2S near the metal sulfide will be held to a very low value. This would permit the reduction reaction to proceed with little or no inhibition from a buildup of reaction product gas. It is well known that calcium oxide is capable of removing hydrogen sulfide from a hydrogen gas stream of low dew point.4 If a sufficient quantity of calcium oxide is mixed with the metal sulfide particles the reaction: CaO+H2S=CaS+ H2O [l] will substitute moisture in place of hydrogen sulfide in the gas stream and this will not affect, in a direct manner, the reaction: MeS +H2=Me + H2S [2] A convenient method of considering the thermodynamics of the combined reducing agent-scavenging agent system is to consider the atmosphere when the partial pressure of hydrogen sulfide is the same over both the metal sulfide and the scavenging agent, i.e., pH2S (1) =pH2S (2). As a consequence: pH2O (1) pH2(2) =K1K2 The chemical driving force for reduction will depend inversely upon the moisture content of the gas and will be 0 when, in the system, pH2O = pH2.K1K2. Table I lists values of the equilibrium constants for reduction and H2S scavenging reactions for a number of sulfides at several temperatures. Data are taken from Rosenqvist4,5 and Kelly.6 The equilibrium constant products calculated from this data show that the limiting level of gaseous reaction product has been increased by a factor of 10' to l04 as a result of substituting a reducing agent-scavenging agent system for a simple reducing agent system. One possible side effect which must be considered is the possibility that the moisture evolved in the scavenging reaction might cause the atmosphere in the system to be sufficiently oxidizing to favor the formation of oxide rather than metal. This possibility was examined by comparing the equilibrium constant products listed in Table I with equilibrium constants for hydrogen reduction of the respective metal oxides. It was found that for copper, nickel, and cobalt the combined reduction-scavenging reactions could not develop a sufficiently high oxidizing potential in the

Jan 1, 1970

By O. B. Bennett

WHEN diamond drilling was introduced in the Rhokana mines in 1939 it was used principally for pillar removal and for completion of the upper portions of shrinkage stopes which were being affected by increasing pressure. This method of drilling long blastholes proved so successful that it was extended gradually to cover stoping, pillar recovery, and hanging cave work. BY 1949 virtually all the ~roduction of Mindola and Nkana was being obtained by this method. At the present time 87,500 ft are drilled each month by the 80 diamond drills in daily operation. Responsibility for control and issue of diamond drilling equipment and crowns, as well as tabulation of all performance figures, was taken over by a sPecially formed Roto drill department, which also investigated the problems encountered with this new method. To assist this department a fully equipped test chamber, Fig. 1, was established underground where performances of various types of machines and equipment could be studied under conditions as nearly uniform as possible. Since the establishment of this department, which was eventually taken over and incorporated into the study department, considerable experimental work has been done on every aspect of the subject. The problems can be classified broadly under four headings: improvement of drilling equipment, crown design, machines, and stoping layouts. One of the major problems with drilling equipment has been to eliminate vibration. Owing to flexing of rods in the hole, severe friction is set up on the back end of the 'Ore barrel and On any high spots in the rods, inducing harmonic vibration in the string of rods and causing the crown to chatter against the face. This not only causes premature crown failure but also reduces penetration speeds and increases wear on the machines and rods used. In the early days, when only holes of EX size were drilled, vibration was largely overcome by periodic greasing of rods and core barrel during each run, but with the change-over to the larger BX hole it became obvious that application of grease by hand was inefficient and time-consuming, and attempts were made to perfect a self-lubricating core barrel. A series of these core barrels was made up and tested and a number of the latest type were used under normal operating conditions, but although footages up to 120 ft were drilled without refilling the overall performance was inconsistent, and the idea was shelved in view of the success of the stabilizer rods referred to later in this paper. At the same time tests were made with barrels 5 ft and later 6 ft long instead of the normal 2 ft. Although a slight improvement was noticed, greasing was still necessary. It was found that rod vibration increased as the core barrel became worn, and in an early test chamber experiment crowns drilled with a worn core barrel averaged 95 ft with a diamond loss of 4.76 carats, whereas the same type of crowns with a new barrel averaged 228 ft with a diamond loss of 3.13 carats. until then all BX drilling had been done with B-sized rods, but during a test on a string of BX-sized rods it was noticed that vibration was negligible. Because of the larger surface area of metal bearing on the sides of the hole, however, the friction and resistance of rods of this size rendered them impracticable on any but the most powerful of the machines, The use of stabilizers spaced evenly along the rods was the next logical step, and for this B couplings, see Fig. 2, were set with three tungsten carbide inserts 1 in. long placed around the periphery equidistantly and at an angle of 45" with a right hand lead. These were placed immediately behind the core barrel and then at 12-ft intervals, as it was found that vibration still occurred when the stabilizers were more than 15 ft apart. The effect of these stabilizers was immediately noticeable; holes were drilled with a minimum of vibration, penetration speeds were improved, and as it was no longer necessary to grease the rods there was a marked decrease in the overall drilling time for each hole. While tests were being made with the stabilizer comeb periodic were taking place with a set of tapered threaded rods, and because there was marked improvement in efficiency it was decided to incorporate the stabilizers and tapered threading in all new rods ordered for Rhokana. The feature of these rods is that only four full turns are required to tighten the coupling as against nine for the present type of B rods. Also, as they are self-centering it is virtually impossible to crossthread them. Each rod has a male 5" tapered Acme thread, Fig, 3, on one end and a female at the other, so that separate couplings are unnecessary, and every fifth rod has an

Jan 1, 1955

By D. Turnbull, J. C. Fisher

ACCORDING to the "reaction-path" modell,2 of martensite nucleation, the shear angle of the embryonic martensite plate must be treated as a variable, and included in any calculation of nucleus critical size. Also, as can be deduced from this model, the interfacial free energy between austenite and martensite does not reach its final value until the shear is completed. It is zero for zero shear angle. However, in order to account for the kinetics of the martensite transformation, some sort of interfacial energy barrier appears to be necessary even with the reaction-path model, for otherwise the volume and the energy of formation of the critical size nucleus both collapse to zero.3 Cohen independently suggested that surface energy could be incorporated into the reaction-path model, with the overall free energy of a martensite embryo being a function of its volume and shear angle.' It is possible to estimate the energy associated with the formation of a critical-size martensite nucleus starting with the reaction-path model and including a surface free-energy barrier. As the dependence of interfacial free energy upon shear angle is unknown, a simple type of dependence will be assumed, with the belief that the true dependence would not lead to appreciably different results. Consider the work required to form a lenticular martensite plate with radius r, thickness t, and shear angle 8. There are three contributions; one being the interfacial free energy, one being the free energy change in the martensite plate, and one being the free energy increase in the surrounding austenite. The interfacial free energy u is assumed to depend upon the shear angle 0 according to the relationship s=s0(?/?0)n [1] where 8, is the equilibrium shear angle and n is an exponent that may lie in the range 0 n 2. The work required to form the interfaces of a martensite plate then is W. = 2pr² s0(?/?0)n [2] The free energy change per unit volume of martensite is composed of two parts, one the ordinary volume free energy ?f1. which is negative, and the other the elastic strain energy G?m²/2, where G is the shear modulus and 7, the shear strain relative to the martensite structure. This expression for the strain energy is valid only when the shear strain ym, is sufficiently small that the martensite is within its linear elastic range. There is no doubt that ym, lies beyond the linear elastic range for embryos that are considerably subcritical. However, for critical nuclei it will be shown that ym, is 1.5 pct or less, within the linear elastic range of martensite. For embryos of nearly critical size, then, the strain energy of the martensite is correctly given by G?m²/2. The shear strain in the martensite is ym, = 8, — 8, and the work required to form the strained martensite is Wm --= (pr²t/2) [?fv + G(?O - ?)²/2] [3] The free energy change in the austenite is entirely that due to elastic distortion. The elastic strain is not uniformly distributed in the austenite, being large near the martensite plate and small elsewhere. Approximately, however, the energy corresponds to a uniform shear strain ya= (?t/2)/r [4] throughout the volume 4pr³/3 surrounding the plate. The work required to strain the surrounding austenite then is Wa = (4pr³/3) (G?a²/2) = (G?²/6) prt² [51 For simplicity, the same shear modulus G is assumed for each structure. The total free energy for forming a plate then is W = W3 + Wm + Wa. = 2pr² s0 (?/p?0)n + (pr²t/2) [?fr+G(?0-?)²/2] + (G?²6) prt2 [6] This expression is correct for nuclei and for embryos of nearly critical size, where, as will be shown, the strain energy in the martensite is correctly given by the expression G (? — ?)². Having W as a function of r, t, and 8, as in Eq. 6, there is a saddle-point where W has a stationary value, W subsequently decreasing indefinitely as the nucleus volume increases along the reaction path. The stationary value of W is the energy of the critical nucleus. The critical nucleus has radius, thickness, and shear angle such that ?W/?r - awlat: = ?W/?p? = 0. Performing these differentiations and calculating the critical nucleus energy, W* = [8192p(G?/6)²;s/27 ?fv4] [7] where a= (?/?0)3n+1[l +G(8"-8)'/2af.]' [7a] and where 8 is to be determined from the equation (1 + 3n/4) + G8(6O - (9)/[Af. +G(6>o-6>)72] = 0 [8] For ?f, near —200 cal per mol or —10" ergs per cc, and 8, near 1/6, as for iron-base alloys, Eq. 8 gives ?0 - ? ~ - (4 + 3n) ?f1./4G0O [9] as the difference between the equilibrium shear angle and the actual shear angle for a critical nu-

Jan 1, 1954

By M. J. Spendlove, H. W. St. Clair

The problem frequently arises, particularly in refining metals or smelting scrap metals, of separating metals in the metallie state. Many metals may be separated by taking advantage of their difference in vapor pressure. Such separations can be made at atmospheric pressure, but the separations are much more selective and can be carried out at considerably lower temperatures if the distillation is done at pressures of a few millimeters or less in an evacuated enclosure. Until recently, this has not been considered feasible as a metallurgical operation, but the recent improvemcnts that have been made in vacuum technology have broadened the applicability of vacuum processes and have prompted re-examination of low-pressurc distillation of metals as a practicable process. The distillation of zinc from lead is one separation that has already been reduced to practice.l This paper is the first of a series of studies being made on separation of nonferrous metals by distillation at low pressures. Although these experiments were confined to the separation of zinc from aluminum, the significance of the results is by no means confined to these two metals. The purpose has been to investigate a metallurgical technique rather than merely to devise a means of separating specific metals. The experimental work on distillation of zinc from zine-aluminum alloys at reduced pressure grew out of earlier work on distillation at atmospheric pressure.2 The earlier work indicated that it would not be practicable to decrease the zinc in the alloy much below 10 pct owing to the high temperature required. Therefore attention was turned to distillation ah low pressures, at which lower temperatures are required. After preliminary tests were made in a small, evacuated tube furnace, a larger furnace having a capacity of 100 to 150 Ib of metal per charge was constructed. Distillation tests were first made on pure zinc and then on aluminum-zinc alloys of various composition. Particular attention was given to the limit to which zinc could be reduced in the residual metal. Data were also taken on the rate of evaporation, and heat balances were made for both the crucible and the condenser. Distillation Furnace The vacuum-distillation unit is illustrated schematically in Fig 1. The major components are the induction furnace, the condenser, the vacuum system, and the power-conversion unit. Power is supplied to the induction furnace from a 50-kw 3000-cycle motor-driven alternator. The pressure in the furnace is reduced by a vacuum pump having a nominal pumping speed of 10 liters per sec. When in operation, the metal vapors travel upward from the furnace to the water-cooled condenser where they are collected in amounts of 50 to 100 lb. The condenser is removed with aid of an electric hoist. When the system is under vacuum, the condenser is made self-sealing by a rubber gasket between the smooth-faced, water-cooled flanges at the top of the furnace and the bottom of the condenser. The pressure of the atmosphere is more than sufficient to insure sealing. At the conclusion of an experiment, the residual metal can be removed from the furnace by removing the condenser and tilting the furnace with the electric hoist. The metal was cast into the molds carried on a mold truck. A photograph of the furnace and auxiliary equipment is shown in Fig 2. The details of the vacuum furnace are illustrated in Fig 3. The furnace proper is made vacuum-tight with rubber gaskets placed at each end of a fused quartz cylinder. A clamping plate at the bottom and a ring at the top are made to squeeze the rubber between the metal and the end of the quartz tube. A large graphite crucible placed inside the quartz cylinder is thermally insulated and physically supported by refractory insulating bricks. A thermocouple in a quartz protection tube is located at the bottom of the crucible: the leads pass through a rubber seal in the bottom plate. The supporting structure for the furnace is an angle iron frame with transite board sides. The condenser is made in the form of a water jacketed cylinder with an opening to the vacuum line at the top. The bottom has a projecting skirt inside the machined flange to provide additional cooling for the rubber gasket. Condenser sleeves are made in the form of two semicylindrical pieces of sheet metal that fit snugly inside the cooling jacket. The split sleeve facilitates removal of the condensate. Measurement of Temperatare and Pressure The metal temperature was measured by a platinum-platinilm rhodium thermocouple inserted in a well extending up into the bottom of the graphite crucible. During rapid evaporation there is a wide difference in temperature between the surface and the main body of metal in the crucible because of the large amount of heat that must be conducted to the surface to supply the heat of evaporation. The heat of

Jan 1, 1950

By H. H. Kaveler, Z. Z. Hunter

Variation of the horizontal permeability (parallel to the bedding plane) in the vertical section of reservoir rocks has long been observed as a characteristic of a normally heterogeneous system which reservoir rock represent. The use of a recently developed water injection profile device offered opportunity to measure with a high degree of reliability the rate of inflow of water into Burbank sandstone in wells previously cored. Water injection profiles were not correlative with core permeability profiles in such wells. Apparently the vertical permeability substantially influences the flow between strata in a formation in a manner as to void the usual conclusions that have been drawn from consideration of the horizontal permeability measurements alone. The results obtained in comparing water injection profiles with horizontal permeability profiles suggest that many of the usual production operations based upon "selective" behavior or treatment of rock exposed in well bores need to he re-valuated and re-examined. INTRODUCTION Petroleum reservoir rock are heterogeneous systems. Heterogeneity exists in respect to lithologic character insofar as such rock are composed of distinguishable solid phases. Heterogeneity also exists in respect to certain properties, such as porosity and permeability, that vary due to variation of the physi-cal structure of the rock. Except in exceptional cases, both the horizontal permeability (measured parallel to the bedding planes) and the vertical permeability (measured perpendirularly to the bedding planes) exhibit significant variation in any common source of supply. The variation in horizontal permeability. as reflected by con. analyses. has drawn the greatest attention of petroleum technologists probably out of the general notion that the mass movement of fluids in a reservoir is predimonantly in the horizontal direction. Furthermore, in the usual case, the rock permeability measured in the horizontal direction is greater than that in the vertical. The variation of horizontal permeability of reservoir rock has been the basis for developing a number of operating practices and procedures intended to improve the petroleum production operation. Many such procedures are referred to as "selective" in the sense that the practice is intended to control the flow to a more. or less. permeable interval within the common source of Supply. It is often said that such practices are "tailored" to the permeability profile. The practices referred to involve, among others, the following: selective perforation of casing; selective shooting, acidizing and plugging: plugging back to intervals of low permeability; and, regulation of flow to prevent coning of water or gas, or irregular encroachment of water or gas. Certain expressed notions involving a concept of "by-passing," or "trappingl" that are held to be particularly harmful in causing the avoidable loss of recoverable petroleum have grown from observed variations in the horizontal permeability. Oftentimes estimates of the reserve of a common source of supply are tempered by conclusions relating variation in horizontal permeability to recover-ability of the oil-in-place. Certain conclusions attributed to the significance of the variation of the horizontal Permeabilitv often extend to the design and operation of pressure-maintenance projects involving both water flooding and gas-injection. Many advocate increasing the number of injection wells, advocate maintaining uniform and equidistant input-output well patterns, or advocate so-called "dispersed" gas-drive techniques rather than gas-cap injection because the permeability profile of cored wells is supposed to indicate that "by-passing" or "trapping" would otherwise exist. It is important, therefore, to have an opportunity to test whether the variation in the horizontal permeability found through core analyses of a typical reservoir rock is sufficient to establish the paths of fluid flow in a reservoir. It is particularly important to have an opportunity to determine whether flow at the sand face of a well conforms to the permeability profile as established by core analyses. In that manner, the merit of certain 50-called "selective" operating procedures and other notions may be evaluated. The purpose of this paper is to compare horizontal permeability profiles of wells in the Bartlesville (Bur-bank) sandstone with water injection profiles, for the purpose of showing that there is no correlation between the horizontal permeability of a core and the water intake characteristics of a typical sandstone. GENERAL CHARACTERISTICS OF BARTLESVILLE (BURBANK) SANDSTONE The Bartlesville sandstones of Northeastern Oklahoma are off-shore bar deposits.' Although other reservoirs had different processes associated with their deposition or with the formation of their porous, permeable structure, the l!artlesville sandstones on which these field Fields were made are, in every respect. typical petroleum reservoir rock. The permeability of the Bartlesville sandstones shows a typical variation in both the horizontal and vertical direction. Furthermore, the permeability profile logs of wells in any pool are not correlative, even as between wells as close as 660 ft and 330 ft apart.'. The same condition exists in such sand-tones as the Jones Sand at Shuler' and is the ordinary and usual characteristic of reservoir rock. THE FIELD DATA The data reported herein are those obtained from coring of nine wells on the center of ten-acre locations for the purpose of providing water-injection wells in the Bartlesville (Burbank)

Jan 1, 1952

By K. G. Wikle, W. W. Beaver

The factors which control the rate of dissolution of pure gold in cyanide solution were studied both directly and through measurement of solution the current-potential curves for the anodic and cathodic portions of the reaction. The mechanism of dissolution is probably electrochemical the reaction in nature, and the rate is determined by the rate of diffusion of dissolved oxygen or cyanide to the gold surface, depending on their relative concentrations. The significance of the results and the effects of impurities are considered. ALTHOUGH the dissolution of gold in aerated cyanide solutions has been used as an industrial process for treatment of gold ores since the late nineteenth century, the factors which determine the rate of the reaction have never been identified unambiguously. Studies of the rate of dissolution by Maclaurin,1 White,2 Christy,3 Beyers,4 Thompson,6 and others are contradictory in their conclusions; some claiming that diffusion of the reactants to the gold. surface controls the rate, and others that the chemical reaction is inherently slow and related to high activation energy for the reaction. Christy3 and 'Thompson" both suggest that the reaction is electrochemical in nature and that the dissolution of gold proceeds at local anodic regions while the oxygen is reduced at cathodic regions on the gold surface. Although their studies are ingenious and do indicate an electrochemical reaction under the conditions of study, their experiments were of limited nature and failed to identify the rate-controlling process in the system. The importance from an industrial viewpoint of a knowledge of the mechanism and rate-controlling factors in gold dissolution can be illustrated as follows: If the rate is controlled by a slow chemical reaction rather than by diffusion of the reactants, then an increased temperature should have a marked accelerating effect; agitation of the slurry should have no effect on rate: and increased concentration of reactants should cause acceleration of the rate. If the rate is controlled by the diffusion of one or the other of the reactants to the gold surface, then increased agitation should increase the rate; increased temperature will increase the rate, but not as much as for the case of a slow chemical reaction; increased concentration of the reactant which is diffusion limited will increase the rate; and the concentration of other reactants should be without effect on the rate. It may be concluded that for design of a commercial process for gold leaching, the rate-controlling factors of the reaction should be understood so that an intelligent choice of the conditions of agitation, temperature, and reactant concentration may be made. The experiments described here lead to the unambiguous conclusion that in a system of pure gold and a pure aerated cyanide solution the rate of dissolution is controlled either by the rate of diffusion of dissolved oxygen or cyanide to the gold surface, depending on the relative concentrations of each. There is also ample, but not conclusive, evidence that the mechanism of the reaction is identical to that of electrochemical corrosion. The practical significance of these conclusions will be discussed later in the paper. Experimental The experimental method used in this work was to employ an electrolytic cell which performed the overall gold-dissolution reaction, and to study the anodic and cathodic reactions of this cell as to their nature and the rate-controlling factors. Simple experiments on the rate of dissolution and the potential of the dissolving specimen also were performed under conditions of agitation, temperature, and concentration identical to those used in the electrode studies. Analysis of the electrode studies by well established theories of electrochemical corrosion were made, and the results were found to bear a one-to-one relation with actual rate and potential measurements. Electrode Studies: The Anodic Reaction: The gold specimen used for all of the electrode studies and the rate determination consisted of a sheet of 99.99 + pct Au wrapped around a lucite rod and sealed at the edges with plastic cement, thus forming a cylinder of gold of known and constant area (8.0 sq cm). The lucite rod was threaded into a brass spindle which could be rotated at speeds of 100, 300, and 500 rpm. For the electrode studies electrical contact between the gold cylinder and the brass spindle was made by means of a gold strip covered with plastic. The anodic dissolution of gold was studied by immersing the electrode in a solution containing known concentrations of KCN and KAu(CN)2 but free of oxygen, and by passing an anodic current through the gold electrode. The pH of the solution was maintained between 10.5 to 11.0 in these and all other tests by addition of KOH. The pH was measured before and after each test by means of a glass-elec-

Jan 1, 1955

By J. A. Barr, R. B. Burt, I. S. Tillotson

THE pumping of solids in water suspension is an important part of many metallurgical and mining operations. In most cases, it is still in the rule of thumb category for which no universal formula has been developed, and much research is needed. Because of the limited and incomplete data available, this article may be classed as an experience paper, which is presented with the hope that some contribution will be made toward the development of the so-called universal formula. This formula, if and when developed, may be evolved from several factors, many of which are not now available for general application. The designing engineer is interested in obtaining accurate forecasts on: 1—the minimum velocities needed to prevent choke-ups in the pipeline, which in turn dictates pipe sizes, 2—power required for pumping, 3—pump selection. The basic factors for a given problem will include: 1—weight per unit of time of solids to be handled, 2—specific gravity of solids, for calculation of volume, friction and power, 3—screen analysis of solids with the colloidal acting, i.e., the slime fraction, a very important factor, 4— shape of particle or some means of determining a friction constant, 5—effects of percentage of solids, 6—development of a viscosity factor to be used in the overall calculations, 7—calculation of the lower limits of pipeline velocities permissible, 8—calculation of total head, pump horsepower, and 9—setting up of pump specifications. In certain limited cases horsepower and total heads and minimum velocities may be computed and a suitable pump selected from basic data, but in many cases, as in mining of Florida pebble phosphate, experience rather than a hydraulic formula still should be used as a basis of selection. Pumping Florida Pebble Matrix Pumping at the Noralyn mine of International Minerals and Chemical Corp. will be used as an example. Other areas will vary as to the characteristics of the matrix, especially the slime content. A typical screen analysis of this matrix is: +14 mesh, pebble size,* 2.1 pct; —14 +35 mesh, 11.4 pct; -35 +I50 mesh, 60.5 pct; -150 mesh, 25.0; total, 100 pct; moisture in bank, 20.0 pct; weight per cu ft in bank, 120 lb. The —150 mesh fraction may increase to as much as 35 pct in adjacent areas. When thoroughly elutriated, the matrix has a relatively slow settling rate, which is an important factor in permitting lower pipeline velocities without choke-ups. Exact data is not available to evaluate settling rates. For a factor of 100 a suspension of clean building sand in water is suggested. When pumping long distances, a quick settling matrix allows the coarser solids to settle out along the bottom of the pipeline, causing drag, turbulence, and increased friction. With a slow settling matrix as at Noralyn, turbulence acts to keep the solids in suspension at a lower friction head, regardless of the pumping distance. When the pebble content of the matrix, i.e., the + 14 mesh fraction, is in excess of 10 pct of the total solids, trouble may be expected from settling out even in normal pumping distances. To prevent choke-ups and maintain tonnage, an additional pump must be added in the long runs, where one pump would otherwise be satisfactory. A typical pulp handled is: total volume, 7800 gpm; water, 4500; solids pumped per hr, 4200 lb; sp gr pulp, 1.4; percent solids in pulp, 46.; pipe size, 16-in. ID; pulp velocity, 12.85 fps; probable critical velocity, 10 fps, as below this minimum velocity choke-ups would be numerous. In calculating friction heads the Armco handbook is used where a roughness factor based on 15-year-old pipe is set up. Because the pipe used in pumping matrix is smooth and polished because of the scouring action of the phosphate and its silica content, the head losses in the Armco table for water are practically the same as in pumping the Noralyn matrix through smooth pipe, plus the fact that conditions vary widely over short periods, making accurate determinations difficult to obtain. New pumps and pump changes are being tested continuously and a wealth of data built up. This has resulted in a substantial improvement and lower relative costs in pumping matrix. The Florida phosphate industry is constantly seeking to offset higher wage and material costs with improved technique. Until a few years ago a 12-in. discharge pump was commonly used, with heads as low as 80 ft. Sizes have gradually increased and heads more than doubled. For example, the following pump was placed under test at the Noralyn mine: make, Georgia Iron Works; size, suction 16 in., discharge 14 in.; impeller, 39-in. diam; motor, 600 hp, slip ring; full load speed, 514 rpm. The results were increased head, higher capacity than the older design, with fewer pumps in the line from mine to washer.

Jan 1, 1953