Search Documents

By Benjamin C. Koskiniemi

Introduction When evaluating any ore body, one of the first questions concerns the ore reserves. In the case of an open pit mine, this is not possible to answer reliably until the ultimate (final) pit limits have been established. Techniques used in designing an ultimate pit are classed as (1) manual, (2) computer, and (3) combined manual-computer. This chapter will describe how manual techniques can be utilized in designing an ulti¬mate pit Certain economic and design criteria must be established before the actual design begins In order to begin designing an ultimate pit, it will be assumed that the engineer already has the following data available (1) vertical sections, (2) horizontal sections for each level, (3) stripping curve, (4) bench height, (5) bank slope angle between levels, (6) level berm width, (7) roadway width, (8) pit slope angles at ultimate pit limits (estimated average including roads and ramps and between roads and ramps), and (9) minimum width of pit bottom

Jan 1, 1979

By B. G. Matson, M. A. Whitfield, G. R. Dysart

This paper describes the development of u computer program to calculate changes that occur in the length of tubular goods due to temperature and pressure changes during stimulation operations. Due to the numerous variables involved and the uncertainty of all static and dynamic conditions that could exist, it becomes a staggering task for individuals charged with completions to perform the necessary mathematical calculations. The computer program permits advance calculations for several sets of conditions. INTRODUCTION In the Delaware basin of West Texas alone, 50 wells were contracted or drilled to 15,000 ft or deeper in 1965. Deep well activity is continuing in this and other areas on an expanding scale. Many of these deep wells require extensive stimulation for successful commercial production, and during these operations, pressures and temperatures are encountered that have a pronounced effect on the length of tubular goods. This length change during a large-volume, high-pressure stimulation treatment utilizing fluids considerably cooler than bottom-hole temperature can be of such a magnitude that permanent damage to casing and tubing will result unless mechanical design, pressures and fluid temperatures are evaluated and controlled. These pressure and temperature effects can be calculated. However, the process lends itself well to computer solutions because of the mathematical nature of the problem and the calculating hours involved in arriving at an answer. The engineering-hour demand becomes more severe as tapered strings are involved. On initial treatments on a given well, surface pressure and injection rate conditions are unknown, and offset well conditions have not proven to be a reliable method for making predictions. For these reasons, it has become rather standard procedure for operators to compensate for these uncertainties by placing unnecessary pressure and fluid temperature restrictions on stimulation design. On a number of occasions treating fluids have been preheated to as much as 160F as a means of compensating for thermal contmction resulting from pumping cool fluids. The maintenance of packer seals has been treated by Lubinski, Althouse and Logan',' and the problem of therma1 effects on pipe has been explored by Ramey." These works were expanded and the results made applicable to everyday oilfield terminology before submitting them to computer programming. The pressure and temperature effects on tubing movement previously mentioned occur simultaneously as fluid moves through the pipe. The pressure changes, for purposes of explanation, are categorized here as to the various effects these pressures have on a tubing string. These divisions are (1) the piston-like results of forces acting on horizontal surfaces exposed to pressure, (2) swelling or ballooning of the tubing along its entire length due to the forces of pressure acting against the tubing walls, (3) the elongation of tubing due to frictional drag and (4) corkscrewing of the pipe due to internal pressure. Thermal changes are also of great importance, as their results may be more significant than any of the pressure effects. Steel is an excellent conductor of heat and the earth is a relatively poor conductor. It has been calculated that pipe temperatures at depths of more than 20,000 ft approach within as little as 25" the temperfature of the surface fluid after pumping for 2 hours, or a drop in temperature in some treatments of more than 220F. The equations presented in this paper were developed for computer programming and simplicity of input information; therefore, numerical constants such as Young's modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson's ratio for steel (0.3) are included with unit conversion factors. The moment of inertia of tubing cross-sectional area with respect to its diameter was changed to a constant times (D' — d') where D is outer diameter and d is inner diameter. Units in the equations are length in feet, diameter in inches, density in pounds per gallon, pressure in psi, rate in barrels per minute and time in hours. PISTON-LIKE REACTIONS A change in tubing internal dimensions and the exposure of other horizontal surfaces to different pressures on the inside and outside of the tubing result in a reaction much like a piston under pressure. Such is the case when the internal diameter changes in a combination string of pipe, when seals of a slick joint assembly are subject to pressure and in the end effects of a tubing string. The change in tubing length due to the piston effects of a slick joint packer is affected by the various diameters involved, the tubing pressure Ap,, the casing pressure ,Ap,, length of pipe L, densities of fluid in the tubing before and during pump-

By R. P. Alger, M. P. Tixier, D. R. Tanguy

In the combination induction-electrical log used at present in the field, the induction logging tool is appropriate for the investigation of moderately invaded formations. A new induction sonde with a radius and investigation about twice as large has been developed recently for the case of deep invasion. It has very nearly the same vertical resolution as the present sonde so that thin beds are defined as accurately as before. The characteristics of the new tool are described, the corresponding interpretation charts are given and field examples are discussed. The design of the sonic logging tool has been modified to improve the calibration and the reliability. The fact that porosity can be accurately recorded by means of the sonic log has prompted new interpretation procedures for saturation estimation, wherein the data concerning the various permeable beds in a given well are correlated. One approach consists of plotting transit time vs true resistivity, with an appropriate scale. With this approach, saturations can be estimated conveniently even in cases where formation water resistivity is not well known. In another approach, a comparison is made of the values of the formation waters computed from the re- sistivity and sonic logs. Using the concept of continuity, this procedure makes possible a quick determination of zones of saturation in shaly sands and/or in case of appreciable variations of formation salinities with depth. It has been found that the comparison of porosity from the sonic log with the apparent porosity computed from a short-investigation resistivity log may reveal, in many cases, the presence of residual oil and thus detect potentially productive formations; this procedure is valuable when the true formation resistivity and the resistivity of the formation water are in doubt. INTRODUCTION During the past year, the efficiency of log interpretation has been vastly improved. The improvements have largely resulted from the introduction of a deep-investigation induction device and from the application of new interpretation techniques that utilize sonic vs resistivity readings. Since the new interpretation techniques depend, in part, upon good values of true formation resistivity, the new induction log will be discussed under Part I. The sonic interpretation techniques will be studied under Part 11. Early in 1959, the 6FF40 induction equipment was introduced in the field. This device was designed for a better approach to true formation resistivities in deeply invaded zones. The greatly improved radial investigation of the 6FF40 equipment has been achieved without sacrificing vertical resolution. The first combination induction-electrical log, the 5FF40, was introduced as a standard tool in 1956 for the logging of wells drilled with fresh muds. The tool has received wide industry acceptance in the United States. 'The 5FF40 induction log has a radial investigation sufficient to overcome average depths of mud filtrate invasion. At 5d invasion, for example, the 5FF40 induction log will read about 1.4 R. in a water sand where R., = 10R,. At 10d invasion, such induction log would read 2.45 R. in the same water sand. In either case, the effects of invasion would not be sufficiently great to cause a water sand to be mistaken for a shale-free oil- or gas-producing zone. Some formations, however, invade deeply — in excess of 10d. Such water zones could be mistaken for oil-or gas-saturated sands unless the porosity balance' can clearly make the distinction. It is for these deeply invaded formations that the 6FF40 was developed. CHARACTERISTICS OF THE 6FF40 Radial Investigation Characteristics To describe the comparative responses of the 5FF40

By A. R. Gregory, M. R. J. Wyllie

The literature reveals that scant attention has been paid to the systematic experimental determination of the formation factors of unconsolidated porous media. No experiments appear to have been made on the effect of increasing cementation on the formation factor of an initially unconsolidated porous medium. Measurements have accordingly been made of formation factors as a function of porosity for aggregates of spheres in the porobity range 12-56 per cent and of cubes, cylinders, discs and triangular prisms in the porosity range 30-45 per cent. The results are examined in the light of the theoretical equations of Clerk Maxwell. Fricke and Slawinski. Aggregates of unconsolidated spheres and beach sands have been artificially cemented with silica and the formation factor-porosity relationships determined. A theory is outlined which seeks to explain the results obtained and which postulates that formation factor and porosity data for consolidated sandatones may be used to determine the original porosity of the unconsolidated sands from which the consolidated sandstones were derived. It follows also that for consolidated sandstones F' = C where F is formation factor, fractional porosity and C and k are constants. INTRODUCTION All schemes of electric log interpretation implicitly or explicitly require the use of formation resistivity factors or? as they are now generally called. formation factors. 'The originator both of the concept of formation factor and of the use of formation factor data in electric log interpretation was Archie.' More recently it has been suggested that formation factors are of value in describing fluid flow in porous media. A full discussion of this problem together with adequate reference.: has been given by Wyllie and Spangler.' In spite of the theoretical and practical importance of the subject, there appear to be remarkably few basic formation factor data in the petroleum literature. A wider search reveals that this situation is nut confined to the literature of the to industry. at various times over the last century, attempts have been made to compute the theoretical relationship between the formation factor and porosity of simple systems of unconsolidated particle.. However. the formulae derived are not all identical and the experimental evidence adduced in support of each i; not wholly convincitng. Most of the experimental work has been carried out on system. of unconsolidated spheres, although some attention has been given to systems of spheroids and cylinder. Other simple shapes have received neither theoretical nor experimental investigation. The writers feel that a better understanding of the formation factors of complex systems of practical importance requires a firm foundation of reliable experimental data applicable to simple systems. Accordingly, they have reviewed the relevant literature and have determined experimentally the relationship between the formation factor and porosity of unconsolidated aggregates of particles of regular geometrical shapes. Some experiments have also been carried out to investigate the effect of cementation on the formation factor-porosity relationship. PREVIOUS WORK Although the term formation factor was not used. the relationship between the resistivity of a solution containing dielectric solids and the resistivity of a similar solution not containing solids is an old and important problem. The formation factor of a system of dispersed spheres as a function of porosity was considered by Clerk Maxwell. Maxwell arrived theoretically at the following expression: where F is formation factor and fractional porosity. Lord Rayleigh' also examined the problem of the formation factor of systems of spheres and cylinders and concluded that for spheres Clerk Maxwell'.; expression gave a value of F whirl]. at low porosities particularly, was too small. However,

Jan 1, 1953

By F. Wald, A. J. Rosenberg

The general features of the constitution of the B-Pt system were determined using standard rnetal-lograph~c, thermoanalytic, and X-ray diffraction techniques. Three compound were found. Two of these, Pt3B and Pt,B, are formed by peritectic reactions at 523° and 890°C, respectively. The third, Pt3B,, is congruently melting with a flat maximum at 940°C but decomposes eutectoidally in to Pt,B ant1 boron nt - 600° to 650°C. THE low-temperature allomorph of boron (red, simple rhombohedra1 a boron) is of scientific and technological interest as an elemental semiconductor.' However, the studies of this material have been hampered by its reported instability above 1200"~ which precludes crystal growth from the melt (mp - 2200°C). Crystallization from platinum solutions has been suggested as an alternative crystal-growth technique, but has met with only limited success.' The technique depends upon the existence of a significant difference between the eutectic temperature and the transformation temperature of boron. In order to clarify the conditions for further crystal-growth experiments, we found it desirable to redetermine the main features of the B-Pt phase diagram since previous reports on the system1'5'6'7 are in marked disagreement. EXPERIMENTAL The experimental methods used were thermal analysis, metallography, X-ray analysis, and, to a lesser extent, measurements of microhardness. Most of the alloys were prepared from spectrograph-ically standardized boron obtained from Johnson-Matthey &Co., Ltd. (212 ppm impurities, exclusive of carbon and oxygen) and platinum powder obtained from F. Bishop & Co. (200 ppm impurities, mainly of other platinum group metals). Some alloys were also prepared with very high-purity, float-zone refined boron (99.9999 pct obtained from "Wacker Chemie" and extrahigh-purity platinum (99.999 pct) obtained from Johnson-Matthey & Co., Ltd. The reported results did not depend on the choices of these starting materials. Five-gram alloy specimens containing 10, 20, 25, 27.5, 30, 33.3, 34, 35, 37, 37.5, 38, 39, 40, 41, 42, 43, 45, 50, 55, 60, 70, and 80 at. pct B were made by melting the elements together in boron nitride crucibles using rf heating of a graphite susceptor, either in vacuum or under high-purity argon. All alloys were heated to at least 1800°C for -5 to 15 min. Most of the alloys did not wet the crucibles when the latter were outgassed by preheating under vacuum. In any event, no weight loss was detected after melting, and the nominal composition was assumed for all specimens. Thermal analysis on 2.5-g samples were carried out in boron-nitride crucibles under a vacuum of 5 x X torr. The apparatus was heated in a "Kan-thal A 1" wound furnace, which limited the maximum temperature to about 1100°C. The output of the indicator thermocouple was fed to a dc recorder with a 1-mv full-scale span and an adjustable zero. The apparatus was calibrated repeatedly, using the freezing points of high-purity aluminum, silver, and gold. The results justified the use of the NBS voltage vs temperature tables for Pt/Pt 10 pct Rh thermocouples. All thermal analyses were run at least twice and both the heating and cooling effects were recorded. Most of the alloys had a very strong tendency to supercool. However, the use of mechanical vibration permitted reproducibility within *5°C for all alloys, except in the region around 40 at. pct B. Only the cooling effects are plotted in Fig. 2, since they appear to be more reliable. For metallography, the alloys were cut with a diamond cutting wheel, cast in a polymethacrylate resin, ground and polished with diamond paste, and etched with dilute aqua regia, a common etch for platinum alloys. Both copper and molybdenum radiation were employed to obtain X-ray diffraction data using Debye-Scherrer cameras and a "Norelco" diffractometer Diffractometry with high scanning speeds (1 deg per min) using nickel filtered CuK, radiation was used to identify the main regions of the diagram. However, molybdenum radiation was used for the detection of boron, since the latter showed very strong absorption and fluorescence effects with CuK, radiation. RESULTS AND DISCUSSION Three intermediate compounds, corresponding to the compositions Pt3B, Pt2B, and Pt3B2, were found in the system. Fig. 1 reproduces their X-ray diffraction spectra, together with those of pure boron and pure platinum. As can be seen from the thermal-analysis data in Fig. 2, Pt3B and Pt2B are formed by

Jan 1, 1965

By Arnulf Muan, Klaus Schwerdtfeger

Equilibrium ratios C02/C0 of a gas phase coexisting with selected phase assemblages of the system Fe-Mn-0 have been determined in the temperature range 1000" to 1300°C. The oxygen pressure for the "hfnO" +hfn30, equilibrium and for the "(Fe,hTn)O" + (Fe,Mnh 0* equilibrium at high manganese contents has been determined by electromotive force measurements using stabilized zirconia as a solid electrolyte. The notstoichometry 01' "hTnO" and of "(Fe, iM1z)O" solid solutions has been determined by ther-mog-/avi?netry and by wet-chemical analysis. The data obtained are used to derive activity-composition relations in "(Fe,hfn)O" and (Fe,Mn),O4 solid solutions. WUSTITE "FeO" and manganosite "MnO" form a continuous series of solid solution at high temperatures,' and so do magnetite Fe304 and the high-temperature, cubic modification of Mn304 (Ref. 2) (high hausmannite, -1170). The oxides "FeO" and "MnO" are cation-deficient phases.495 The nonstoi-chiometry of "(Fe,Mn)O" solid solutions has been studied by Engell and ~ohl' at two selected C02/C0 ratios at 1250°C. The two oxide end members of the spinel solid solution, FesO4 and Mn,04, however, are known to be close to stoichiometric under the experimental conditions used in the present investigation.''' The oxygen pressures of "(Fe,Mn)07' solid solutions in equilibrium with iron have been determined by Schenck and coworkers,8 by Foster and welch," and by ~n~e1l.l' The two former groups equilibrated the condensed phases in C02-CO atmospheres of lmown compositions, whereas Engell" used a galvanic cell with stabilized zirconia as a solid electrolyte. The results of these investigators are not in good agreement. Activities of FeO in manganowiistite as calculated from the results of Foster and Welch show ideal behavior, those of Engell yield a pronounced positive deviation, and those of Schenck et 01. show a moderate positive deviation from ideality. In the present work oxygen pressures for the iron + manganowiistite and manganowustite + spinel equilibria and the nonstoichiometry of manganowiistites have been measured. The data were used to calculate activities in the manganowiistite and spinel solid solutions. EXPERIMENTAL METHODS The COz/CO ratios at which manganowustite and iron are in equilibrium were determined by thermo-gravimetric and quenching methods. Experimental details are described in a previous publication.'2 In the thermogravimetric technique, incipient reduction of manganowiistite pellets to metallic iron was observed as a break in the weight vs log COZ/CO curve. In the quenching technique, manganowiistite samples were partially reduced to metallic iron, or the metallic iron of manganowustite + metallic iron mixtures was partially oxidized to manganowustite, in atmospheres of constant C02/CO ratios. After quenching the composition of the oxide phase was determined by X-ray lattice parameter measurements and comparison with a standard curve obtained from oxide solid solutions of known compositions. The nonstoichiometry of "MnO" and "(Fe,Mn)07' solid solutions was determined by chemical analysis of samples equilibrated in C02-CO atmospheres and quenched to room temperature, as well as thermo-gravimetrically by reducing (Fe,Mn),04 or Mn304 to manganowiistite or manganosite. The equilibrium between manganowiistite and (Fe,Mn),04 was measured thermogravimetrically by reducing (Fe,Mn),04 solid solutions having composition in the range of %„ l(NFe +NM) from 0 to 0.63. No experiments could be performed with this technique at higher manganese contents, because the equilibrium C02/C0 ratios are too large for accurate control. An additional difficulty arises at the higher manganese contents due to the strong increase in oxygen content of the manganowustite phase with increasing log Py near the manganowiistite-spinel boundary. Consequently a sharp break in the weight loss vs log C02/CO curve cannot be observed at the phase boundary. At high manganese contents of the manganowiistite, e.g., (NMn/(NF~ + NMn) > 0.9, electromotive force measurements with stabilized zirconia as a solid electrolyte were made to determine the equilibrium oxygen partial pressure. Experimental details are described in a previous paper.* Mixtures of "(Fe,Mn)O" and (Fe,Mn),04 were pressed to pellets, and the oxygen pressure of the equilibrated samples was compared to that of Ni + NiO mixtures in the cell The composition of the manganowiistite in the equilibrated two-phase mixture was determined by lattice parameter measurements and comparison with known standards. The oxygen pressure for the Ni + NiO equilibrium was taken from available data.l3~l4 No reliable results were obtained with the electromotive force technique on iron-rich oxides. The electromotive force drifted strongly with time in this composition range. An additional difficulty arises from the partial de-

Jan 1, 1968

By M. E. Shank, J. Wulff

At the present time, the designer of dies for metal powder pressing is handicapped by relative ignorance of stress distribution and frictional effects at the interior surface of the die. Unckell was the first to develop a method for the study of wall friction. He used three Brinell balls on which the die rested during pressing. The total frictional wall force was determined by the size of impression these balls left on a soft metal plate. Since the method does not give radial pressures, or distribution of such pressures, coefficients of friction could not be determined. Although Unckel measured density distribution, he was not able to determine radial or shear stresses. Shaler2 has proposed theoretical expressions for the stress and density distribution within cylindrical compacts during pressing, in accordance with the experimental results of Kamm, Steinberg, and Wulff.3 By application of Siebel's method,4 Kamm et a13 plotted stress trajectories for two compacts. From the stress trajectories they calculated coefficients of friction from point to point along the die wall. As pointed out by Shaler in the discussion of Ref. 3, these values are based on progressive point-to-point calculations on finite size grid squares across the compact. In the region of the die wall such calculated values may therefore have considerable cumulative error. The purpose of the present paper is to develop an experimental method by which the nonhydrostatic pressures and shears acting on the interior wall of a cylindrical die can be measured. Such measurements can then he correlated with existing data to aid in the explanation of the pressing process. The method used is based on the elastic: properties of the thick-walled tube used as the die. The principle of super-position of force systems on an elastic body is assumed to hold. Electric strain gauges were mounted in adjacent positions on the exterior die wall in order to get an exact measurement of the variation of tangential strain over the length of the die during pressing. While in this paper, measurements in terms of only tangential strains are considered, it is well to note that similar calculations may be set up for axial strains. The latter are not preferred, since they tend to be smaller than the tangential strains and therefore permit less sensitive measurements. Discussion in this work is restricted to compacts pressed from both ends, since the elastic deformation of the die is then more amenable to analysis. Before choosing the electric strain gauge method, a more direct line of attack was considered and discarded. The discarded idea was the insertion of a pressure gauge through a hole in the die wall.* The gauge would have been in the form of a small piston. If pressure were exerted against such a gauge, it would move outward along a radius of the die. One disadvantage of the scheme is its inability to measure shears along the die wall. Another more serious disadvantage is the disturbance caused by the device itself. It would serve to change the forces it was designed to measure. No matter how small the movement of the gauge, when pressure is applied a discontinuity would exist in the wall surface at that point. Due to the stress concentration caused by the hole, abnormal deflections of the die wall would occur around the gauge. During pressing, powder would be forced into the resulting depression. The depression would then become larger with increasing compacting pressure. Powder, not being a fluid, is capable of supporting shear. The ease with which it would flow into the die wall depression to further move the piston is an indication, not of the radial pressure at that point, but of the state of shear retarding the movement. Thus the "pressure" gauge is really a criterion of flowability, and of the capability of the powder to support shear. For these reasons, it was decided that the electric strain method, herein employed, was more reliable, if more indirect. The gauges and lead wires, mounted on the external die wall do not in any way affect the behavior of the metal powder or the die during pressing. Theory of the Method THE EFFECT OF RADIAL PRESSURE ON THE DIE WALL Effect of a Single Small Band of Hydrostatic Pressure Consider a die which is a thick-walled cylinder of outer radius R. and inner radius Ri. If over a small finite length L there is a normal pressure P, a tangential strain distribution at the outer wall results. This is shown schematically in Fig 1. The exact shape of the curve may he predicted by an extension of the theory of a semi-infinite beam on an elastic foundation.6 This

Jan 1, 1950

By Seymour G. Epstein

Using a modified capillary-reservoir technique, electromigration of dilute solutions of silver in liquid bismuth has been measured with relatively good precision. The effects of experiment duration, temperature. current density, and silver concentration have been investigated. Chemical diffusion coefficients of silver in liquid bismuth were also measured at temperatures ranging from 300°to 600°C. The effective valence of silver in bismuth was calculated as a function of temperature from a comparison of electric mobilities and diffusion coefficients. In an alloy containing 1 wt pct Ag, the electric mobility of silver increased from 1.65 x X sq cm per v-sec at 300°C to 2.38 x 10-4at 600°C; the diffusion coefficient in-creased from 3.0 to 17.5 x 10-5 sq cm per sec over the same temperature range. Calculated values of effective valence decreased from 0.28 at 300°C to 0.10 at 600°c. ALTHOUGH there has been considerable progress in the theoretical treatment of electron-transport properties in liquid metals, the mechanisms of mass-transport phenomena in liquid metals are still not clearly understood. This is especially true for the case of electro migration, alternatively termed elec-trotransport or electrodiffusion, which is mass transport induced by an applied electric field. A thorough knowledge of the mechanisms governing electromigration will considerably aid in understanding the electronic state of an impurity atom in a liquid metal. Electromigration has been long observed and studied in liquid metals and alloys. However, no satisfactory theory to explain the observed effects has been developed. The disparity in the present theoretical treatments, summarized by Verhoeven,1 has been attributed to the difficulty in describing the momentum transfer between electrons and ions and to the lack of precise measurements in liquid metals, especially at elevated temperatures. Prior to the last decade most of the electromigration studies were qualitative, in that only the direction of migration of alloy constituents was determined with any degree of reliability. These results have been summarized by schwarz,2 who was able to measure the electric mobility of several elements in mercury. More recently electromigration in liquid metals has been measured by a variety of techniques, with a varying degree of success. Mangelsdorf3,4 made precise measurements on dilute amalgams at room temperature, using a technique relating changes in composition to changes in alloy resistivity, but was unable to make measurements at elevated temperatures. Using a modification of this technique, Verhoeven and Hucke5 measured electromigration and resistivities in liquid Bi-Sn alloys. However, the resistivity technique is not sensitive to small composition changes and is impractical for very dilute alloys. Verhoeven and Hucke6 also measured electromigration of several solutes in liquid bismuth by adapting the capillary-reservoir method for determining diffusion coefficients in liquids. The reported precision of these measurements was poor, however, primarily due to the large uncertainties in the chemical analyses to determine concentration changes. Belashchenko7 has summarized the recent Russian investigations of electromigration, many of which, however, are of questionable accuracy. The present study was initiated to develop an experimental technique which would be both versatile and sufficiently sensitive to yield precise measurements of electromigration of solutes in dilute metal alloys at elevated temperatures. Using this technique, a comprehensive study of silver in liquid bismuth was made to determine the effects of experimental variables on the apparent electric mobilities. Chemical diffusion coefficients of silver in liquid bismuth were also measured over the same temperature range for comparison with the electromigration data. EXPERIMENTAL TECHNIQUE Electromigration Measurements. Electromigration of silver in liquid bismuth was measured by the capillary-reservoir technique, similar to that used by Verhoeven and Hucke,6 but modified to yield more

Jan 1, 1967

By J. L. Taylor, P. Duwez

PARTIAL phase diagrams of titanium with iron, cobalt, and nickel have been established by previous investigators.1-3 These diagrams seem to be reliable, at least for concentrations of titanium ranging from 0 to about 50 at. pct. More recently, Long and his collaborators have presented a tentative titanium-nickel diagram covering the titanium-rich side.' These various studies indicate that in each system three intermediate phases exist and correspond to the compositions Tia, TiX, and TiX2 (or TiX3). The purpose of this paper is to review the available information on the crystal structure of these phases and to give the results of additional X-ray diffraction studies. Previous Work A summary of the work of Wallbaum and his collaborators"-' is presented in table I. The first intermediate phase on the titanium-rich side of the three systems, with iron, cobalt, and nickel, is of the Ti2X type and has a face-centered cubic structure. Although it is stated that the unit cell contains 96 atoms: no lattice parameters are given. The second phase is centered around the compositions TiX for the three systems. The crystal structure of this phase is mentioned in a short note" as being body-centered cubic (CsCl type), but no parameters are given. The analogy between the three systems, which is perfect as far as the two phases of the type Ti2X and TiX are concerned, disappears at the titanium-poor end of the diagrams. In the system titanium- iron, TiFe2 has an hexagonal (MgZn2 type) struc-ture with the parameters given in table I. The existence of TiFe3 has also been reported8 (table I) but was questioned by other investigators2 and is also refuted by the present study. In the system titanium-cobalt, a phase corre-sponding to the composition TiCo2 has been investi-gated by Wallbaum and Witte6 who claim that the structure is cubic (MgCu2 type) for alloys slightly deficient in cobalt and hexagonal (MgNi2 type) for alloys slightly deficient in titanium. The parameters of the two phases are given in table I. In the system titanium-nickel, the nickel-rich phase corresponds to the composition TiNi3 and has been described as hexagonal7 with the parameters given in table I. The present investigation was initiated for the purpose of filling in the gaps in the knowledge con-cerning the structure of the various phases sum-marized in table I. Since the titanium used by Wallbaum and his collaborators was only 95 pct pure,' it was found desirable to redetermine the lattice parameters of the alloys prepared with a

Jan 1, 1951

By K. H. Hiller

A rotational viscometer has been designed which perrnits the measurement of the rheological properties of drilling muds and other non-Newtonian fluids under conditions equivalent to those in a deep borehole (350F, 10,000 psi). The important mechanical features of this instrument are described, and its design criteria are discussed. The flow equations for the novel configuration of the viscometer are derived and the calibration procedures are described. The data and their interpretation, resulting from measurement of the flow properties and static gel strengths of homoionic montmorillonite suspensiom at high temperatures and pressures, are presented. Data are also presented for the flow behavier of typical drilling fluids at high temperatures and pressures. The pressure losses in the drill pipe and the annulus depend critically upon the flow parameters of the drilling fluid. This work demonstrates the need to measure these parameters under bottom-hole conditions in order to obtain a reliable estimate of the pressure losses in the mud system. INTRODUCTION The rheological properties of drilling fluids are affected by temperature and pressure, but the extent of these effects on the dynamic flow properties is not well known. Measurements of changes of the flow properties of clay-water drilling muds with temperature have been reported by Srini-Vasan and Gatlin.1 The temperatures reported did not exceed 200F, a limitation imposed by the apparatus used by these authors. The rheological properties of clay suspensions were measured at temperatures up to 100C by Gurdzhinian.' Neither the nature of the exchange ions in the clay suspensions nor the degree of purity were defined in his work, nor were the measurements extended to currently used drilling fluids. The lack of systematic measurements of dynamic flow properties at high temperatures and pressures seems the more surprising since during the last decade the importance of the control of the hydraulic properties of drilling fluids has come to be widely recognized. Very good mathematical treatments of the friction losses in drill pipe and annulus have been developed.3 4 These treatments are based on the assumption that drilling fluids behave as Bingham plastic fluids. Quite often this assumption is justified, while in other cases a power law equation pro- duces better fit than the Bingham model does. For convenience in applying viscometer data to pressure-drop calculations, the Bingham plastic flow equation is preferable and, therefore, has been applied to the data reported in this paper, although other equations may fit these data more accurately. In a Bingham plastic fluid the relationship between the shearing stress 7 and the rate of shear D is given by the following equation: where is the plastic viscosity and 4 the yield point. If 4 = 0, the equation for simple Newtonian flow, 7 = pD, is obtained. Two empirical constants are required for the description of laminar flow of a Bingham plastic fluid, and calculations of the flow behavior at high temperatures and pressures cannot be better than is permitted by the accuracy with which these constants are known. For this reason a high-pressure, high-temperature rhe-ometer has been designed to measure the plastic viscosity the yield point +, and the static gel strength S, at pressures up to 10,000 psi and temperatures up to 350F. The important features of its design will be described. The results of measurements on homoionic clay slurries will be discussed insofar as they are relevant to an understanding of the general flow behavior of clay-water drilling fluids. The results of measurements on some typical drilling fluids will be presented also, and their practical implications will be briefly discussed. DESCRIPTION OF EQUIPMENT MECHANICAL FEATURES A viscometer designed to measure the plastic viscosity, yield point and gel strength of non-Newtonian fluids must permit the measurement of the shearing stress t at any given rate of shear D. This is possible only if t and D are approximately uniform throughout the entire sheared sample. A Couette apparatus is the most convenient method of realizing this condition, as has been pointed out by Grodde." The "high-pressure, high-temperature rheometer" described in this paper is basically a rotational Couette viscometer that is immersed in a cell in which pressure and temperature can be controlled over the range of interest. Fig. 1 shows schematically the important features of the pressure cell and associated equipment. The heart of the instrument is the rotating cup. It is shown more clearly in Fie. 2. which revresents the lower one-third of the pressure cell (below the input drive shaft shown in Fig. 1), and it is shown in detail in Fig. 3. For measurements of dynamic flow properties, the rotating cup is driven by a 1/2-hp electric motor, which operates through a Vickers

By E. J. Pasierb, H. W. Paxton

Mavkers ill a diffusion couple in the iron-chromium system move towards the high chromium side indicating Dcr > DFe at the marker composition. By a combituztion of measurements of D& and D a: a function of composition, this is also shown independently. DFe is calculated as a function of composition. and shown to be in reasonable agreement with measured values by other workers. MEASUREMENTS of marker movements in diffusion couples have been used from time to time to demonstrate that a vacancy or interstitial mechanism may be operative. Since there is a possibility that ring diffusion1 may be operative in chromium, the present experiments are of interest in showing some non-ring mechanism in iron-chromium alloys around the equiatomic composition. EXPERIMENTAL PROCEDURE The movement of markers originally at the interface between 1-in. sq sections of pure Cr (6 mm thick) and outer slabs of 30 Cr-70 Fe alloy was followed as a function of time at temperatures from 1150" to 1425"C, in a Kanthal-wound furnace. Temperature control was +5°C over extended periods. The analyses of these alloys are given in Table I. Selection of suitable markers presented a considerable experimental problem; the eventual solution was to draw quartz rod down to 0.01 cm and seat this in V-grooves 0.01 cm deep on each side of the chromium prior to welding. A preliminary experiment in which the groove was machined in the alloy gave identical results. Welding was accomplished by passing a heavy current through specimens under moderate pressure in a helium atmosphere. This method produced good welds with a planar interface as shown by subsequent metallo-graphic investigation. Any markers which did not adequately fill their groove and showed subsequent contamination from the atmosphere during heat treatment were not considered reliable. The composite sample was ground and metal-lographically polished on the four sides perpendicular to the interface. Fiducial marks (deep milled slots) were made outside the diffusion zone to check on total volume change. All necessary measurements were made on a metallograph reading to +0.0005 cm. The samples were wrapped in tantalum foil and given diffusion anneals for times to 200 hr in carefully purified dry hydrogen in a sealed McDanel porcelain tube. Chromium-zirconium chips were introduced near the specimen as an additional getter. The specimens were cooled relatively slowly after the diffusion anneal to preserve the McDanel tube but no significant error in effective time at temperature was thereby introduced. At least 1 mm was ground off the surface before measurements were taken. It was observed that the grain size following short diffusion anneals was several mm and was appreciably larger in the 30 Cr alloy than in the pure Cr. No significant grain boundary penetration was present as indicated by the etch limit. The results are presented in the next section. Measurements of D were made by standard sectioning techniques with chemical analysis. Machining difficulties limited the

Jan 1, 1961

By B. A. Eaton, K. E. Brown, C. R. Knowles, D. E. Andrews, I. H. Silberberg

This paper presents the resitlts of an investigation of two-phase, gm-liquid flow in horizontal pipelines. Experimental data were taken in three field-size, horizontal pipelines, two of which were constructed for this purpose. The data were obtained using water, distillate and crude oil separately as the liquid phase, and natural gas as the second phase. Experimental pressure-length traverse, liquid holdup and flow-pattern data were collected for each set of flow rates. These data were used to develop three correlations that are presented herein. The liquid-holdup values correlated with various flow parameters without regard to the existing flow pattern. The same was true for the energy-loss factors. A new flow-pattern map is presented that appears to be quite reliable, but not required for the pressure-loss calculations. The liquid-holdup correlation and the energy-loss factor correlation are used in conjunction with a two-phase flow power balance, developed during this study, to predict the pressure losses that occur during gas-liquid flow in horizontal pipelines. A recommended calculational procedure is given, as well as a statistical analysis of the results. This procedure lends itself to computer application, since several small pressure decrements are needed to calculate a pressure-length traverse. The correlations are shown graphically, but may be curve fitted with existing curve-fitring computer programs. INTRODUCTION Due to the frequent occurrence of gas-liquid flow in pipelines and the desire to accurately calculate the pressure losses that occur in these lines, two-phase flow is of considerable interest to the petroleum, chemical and nuclear industries. In the petroleum industry, gas-liquid mixtures have been transported over relatively long distances in a common line due to the advent of centralized gathering and separation systems. Long two-phase flowlines are usually accompanied by large pressure drops which influence the design of the system. Gas-lift installations are designed on the basis of known tubing pressures at the wellheads. The horizontal flowline connecting the wellhead and the separator system must be correctly sized in order to minimize the horizontal flowline pressure losses and the wellhead tubing pressure. Practically all oilwell production design involves horizontal two-phase flow in pipeknes. All of the flow processes of oil and gas production must be studied simultaneously to insure good well design. Since the beginning of offshore oilfield development, long horizontal flowlines have been constructed. Because pressure losses greatly influence the performance of producing wells, a method is desired that can be used to predict such pressure losses and select optimum flowline size. Several types of gas-liquid flow exist, and many of these are discussed by Gouse The study of pressure gradients, fluid distributions and flow patterns that occur in horizontaI multiphase flow is made difficult by the great number of variables involved. The various flow regimes give rise to changing velocities of the fluid particles in all directions. These instabilities of the interface between the gas and liquid prohibit the determination of actual vector velocities of fluid particles in each phase. Also, it is practically impossible to arrive at correct sets of boundary conditions. Therefore, most investigators have concluded that a solution to the problem by the classical fluid dynamics approach, whereby the Navier-Stokes equationsM are formulated and solved, is far too complex. Other methods must be utilized to develop general correlations that will predict the behavior of gas-liquid horizontal flow systems. Multiphase flow studies have sought to develop a technique with which the pressure drop can be calculated. Pressure losses in two-phase, gas-liquid flow are quite different from those encountered in single-phase flow; in most cases an interface exists and the gas slips past the liquid. The interface may be smooth or have varying degrees of roughness, depending on the flow pattern. Therefore, a transfer of energy from the gaseous phase to the liquid phase may take place while energy is lost from the system through the wetting phase at the pipe wall. Such an energy transfer may be either in the form of heat exchange or of acceleration. Since each phase must flow through a smaller area than if it flowed alone, amazingly high pressure losses occur when compared to single-phase flow. Most investigators of horizontal two-phase flow phenomena have chosen to separate their experimental data into several groups of observed flow patterns or regimes.

By A. Satter, J. M. Campbell

Reported herein are the results of a careful and detailed study of the non-ideal behavior of pure gases and their mixtures. Included are: (1) new data on five ternary systems composed of methane, ethane and H2 S; (2) a simple compressibility factor correlation that is inherently superior to present correlations, particularly for gases containing H2S and CO2; and(3) a detailed study of combination rules and the effect of system composition on the choice thereof. This study makes use of the rather large mass of data already available in the literature. A complete re-examination of the data and ideas presented in the last 25 years was considered desirable as a prelude to our basic concern — the effect of diluents on gas behavior. A consideration of both the macroscopic and microscopic properties of gases provides a better insight which, in turn, gives a firmer basis for improved correlation techniques. Such a study has shown that expressing the compressibility factor Z as a function of acentric factor w, as well as reduced temperature and pressure, yields a correlation that is broader in scope. The study of various combination rules has shown that better results are obtained by "tailoring" the rule used to the system composition. To do so improves the basic reality of results by overcoming some of the anomalies often found when using Kay's rule alone. Tentative recommendations are made regarding the most reliable combination rule for use with a given class of gas. The data presented are useful for estimating the direction and magnitude of the expected deviation when using a given rule. Although more work is needed, particularly around the critical region and with CO2 mixtures, the advantage of the classification scheme proposed is apparent. INTRODUCTION When one attempts to write a PVT equation to fit the data for actual gases, greater precision is obtained by the use of a multiple number of empirical constants. This has lead to multiple-constant equations such as Benedict-Webb-Rubin, Beattie-Bridgman, Keyes, etc., which are capable of yielding very precise results for pure gases in a range for which data to get the constants are available. As a matter of practicality, though, the use of such equations for gas mixtures is limited. Because of the infinite number of gas analyses available, any attempt to compile the constants needed requires a prohibitive amount of experimental data. This could be overcome by the use of a combination rule, but there is no real advantage in doing so because the end result offers no practical impovement over the Z factor correlation. The most widely used method of predicting the volumetric properties of pure gases is based upon the "theorem of corresponding states". According to this theorem, "all pure substances have corresponding molal volume at corresponding temperature and pressure if the reference point of correspondence is the critical point". Generalized compressibility charts for gases were prepared first by Cope and associates1 in 1931 and later by Brown and co-workers2 in 1932. However, the most commonly used charts are those of Dodge,3 Nelson and Obert,4 Hougen and atsson: and Standing and Katz.6 The work of Katz and co-workers has provided us with basic data for the hydrocarbons most widely used today. Their original chart6 was compared with a relatively large amount of multi-component data for gases consisting almost entirely of normal paraffin hydrocarbons. A deviation of only + 1.2 per cent was obtained.39 In the 20 years following publication of this work it has been found that the behavior of most mixtures of paraffin hydrocarbons could be predicted by this correlation within at least 5 per cent. Where difficulty has been encountered it has largely involved one or more of the following circumstances: pressures above 4,000 psig, mixtures containing large amounts of heavy ends and/or aromatics, systems in the critical region and mixtures containing polar compounds and/or CO2. The abnormal error sometimes found with such gases, not too unexpected for this method, is

By George J. Young

(San Francisco meeting, October, 1911.) THE nature of slimes handled in the treatment of gold- and silver-ores has been discussed in technical literature to a considerable extent. The subject of slime-filtration from the practical worker's stand-point has also received much comment, and scattered through the literature of the subject are descriptions of many slime-filtration installations. Articles of this nature serve a valuable purpose and assist materially in the design of new and the improvement of old plants. The subject of the physics of slime-filtration has been touched upon to only a slight extent. The underlying principles are worthy of more intensive study and experimentation than they have received, and the main purpose of this paper is to present the results of such study and experimental work as will serve to make clear in part at least many of the principles which control the filtration of slime. Nature of Slime. Much has bean written concerning an accurate definition of the term ? slime," but no comprehensive definition seems to be generally accepted. The reason for this is clear. A slime consists of at least three different substances, each, when separated, possessing distinctly different physical and to a certain extent chemical properties. These substances are extremely fine sand, a colloidal material which may be and generally is in a coagulated condition, and a colloidal material which is in a non-coagulated condition. Suspending a slime in a relatively large volume of water by shaking and allowing sedimentation to take place, results in the fine sand settling out with comparative rapidity, followed by the coagulated material, which settles much more slowly and finally a certain portion remains

Nov 1, 1911

By E. H. Snow, B. E. Deal

Instability effects in semicanductor devices have long been attributed to the motion of charges on or within oxide layers on the surface. These effects are of critical importance in metal-insulator-semiconductor MIS) field-effect devices. For this reason, the capacitance-voltage or the conductance-voltage characteristic of these devices can be used as a sensitive detector to study charge transport and polarization effects in the oxide or insulator. A review is given of the results of such studies of thermally grown SiO2 as well as of other insulating films of importance in silicon technology. Three types of effects are distinguished. The first of these is the drift of mobile cations within the dielectric, examples being thermally grown SiO, which is often contaminated with sodium ions, and a variety of other glasses in which the mobile ions are a part of the glass composition. The second effect is a dipole-type polarization which occurs in phosphosilicate glass films obtained by reacting P2O5 with thermally grown SiO2. The third effect involves the transfer of charges between the dielectric and the silicon electrode. This occurs in silicon nitride and other deposited dielectrics. It is concluded that MIS studies have provided a powerful technique joy the study of charge transport and polarization effects in insulating films. The knowledge gained from these studies has led to an understanding of surface effects on conventional transistors and diodes as well as making possible stable MIS transistors. THE metal-insulator-semiconductor field-effect transistor is conceptually the oldest type of active semiconductor device.' The earliest attempts at making this device were frustrated because of high surface state densities at the interface between the semiconductor and the gate insulation.' However, by using a silicon substrate with thermally produced silicon dioxide as the gate insulation, this problem was solved and metal-silicon dioxide-silicon devices with good characteristics were made as early as 1960. 3 Yet it was still over 5 years before these devices became a commercial reality. This delay was largly due, not to surface states, but to stability problems associated with polarization effects within the insulating layer which caused the threshold voltage of the device to drift under temperature and bias treatments. The solution to these problems has not only made possible stable MIS devices, but it has added immensely to our understanding of failure mechanisms in conventional bipolar transistors and has added to the reliability of ali types of planar devices. In this review, we shall first describe the effects of various types of polarization phenomena on MIS device characteristics. Then, since thermally grown SiO, is by far the most important insulator used in these devices, we shall review historically the type of instabilities which have been observed in thermal oxides, the attempts at understanding and eliminating them, and the present status of the problem. We shall then turn our attention to the various deposited insulators which have been used, including lead glasses, phosphosilicate glass, vapor-deposited silicon oxide, and silicon nitride. Interestingly enough, many of these materials show polarization effects which are quite different from those generally observed in thermally grown SiO2. THE EFFECTS OF POLARIZATION PHENOMENA ON MIS CHARACTERISTICS The simplest MIS device and the one which has been most frequently used in the study of polarization effects is the MIS capacitor. Two modifications of this structure with single- and double-layer dielectrics are illustrated in Figs. l(a) and (b), respectively. The capacitance of this structure as a function of voltage applied to the metal gate electrode is plotted in Fig. 2 for the case of an n-type silicon substrate. When the silicon surface is accumulated (positive bias) the measured capacitance is just that of the insulating layer C. When the surface is inverted (negative bias), the capacitance is that of the insulator and a silicon depletion layer in series CoCs/(Co + Cs). Indicated on the horizontal axis of Fig. 2 is the voltage VT at which the silicon surface becomes strongly inverted. This voltage corresponds to the threshold or turn-on

Jan 1, 1969

By H. Doi, Y. Fujiwara, K. Miyake

In order to throw light on the mechanism of plastic deformation of WC-Co alloys, compressive tests of WC-(7 to 43) vol pct Co alloys have been carried out at room temperature, and stress-micro strain relation has been investigated in detail. The analysis of the factors affecting the yield stresses reveals that the yield stresses can be predicted by modified Oro-wan's theory if one properly estimates the planar in-terfiarticle spacings. Conzpressive straining of some of the alloys by 0.066 to 0.17pct increases the decrements by a factor of as much as 3.4 to 14, whereas the corresponding increase in the electrical resistivities is less than 10 pct. The analysis of the decrement data in terms of -Gramto and Lücke theory shows that the marked increase is attributed to increased dislocation darnping itt the binder (cobalt) phase. By cornbilling the decrement data and the conzjwession duta, one obtains the relation between flow stress in shear (?t) and increase in dislocation density (p): At = const . v6 . This is interHeted to mean that the mechanism of strain hardening of CirC-Co alloys is essentially sarne as the one for dispersion strengthened alloys. The possible effect of bridge formations between the carbide particles has also been examined. OWING to the combination of hardness, strength, and other physical and chemical properties, WC-Co alloys have opened the way for unique fields of applications, the recent ones being, for instance, anvils for super-high-pressure generation apparatuses. In such applications, the alloys are frequently subjected to very high compressive stresses: these stresses may cause the alloys to deform plastically and eventually to fail. However, much remains obscure regarding the nature of the plasticity of the alloys. Evidently, the alloys owe their high strength to the hard carbide particles which frequently occupy as much as 80 to 90 pct in volume fraction, whereas the ductility required for practical applications is provided by the small amount of the binder phase between the carbide particles. When the volume fraction of the carbide phase is not very large, deformation behavior of the alloys may be described by some of the current dispersion strengthening theories. However, greatly increasing the carbide phase is thought to lead to some carbide skeleton structure or bridge formations owing to the increased chances for direct contacts between the carbide particles;1,2 this may appreciably affect the plasticity of the alloys. Regarding the effect of formation of the carbide skeleton structure, it is interesting to note the work by Ivensen et al.3 on compression tests of the alloys containing somewhat large carbide particles; they observe extensive generation of slip bands in the carbide particles after application of some preliminary compressive stresses. They interpret the results in terms of plastic deformatiot: of the carbide particles which are supposed to have formed a skeleton structure; the binder phase plays only a passive role, at least in the early stages of the deformation. That carbide crystals exhibit microplasticity at room temperature is apparent from the work of Takahashi and Freise4 and French and Thomas5 on indentation of WC single crystals. On the other hand, Dawihl and coworkers6-10 maintain that even when volume fraction of the carbide phase is very large (for instance, more than 90 pet), a very thin binder layer generally exists between the carbide particles. They interpret the results of the extensive mechanical tests in terms of the plasticity of such a layer. Gurland and Bardzil11 point out that decrease in ductility of the alloys with increase in the carbide phase is caused by the effect of plastic constraint exerted by the dispersed carbide particles. Drucker12 further develops this concept from a continuum-mechanics approach on an assumption that a continuous thin binder layer separates the carbide particles. A common feature of the studies reported so far on the plasticity of the alloys is that the information deduced is invariably qualitative in nature. Thus, very few systematic experiments for obtaining reliable and sufficiently detailed stress-strain curves of the alloys varying widely in the microstructural features have been carried out. In particular, it may be of special interest to investigate in detail the early stages of the plastic deformation of the alloys in order to shed light on the strengthening mechanism. However, such work appears to be extremely rare. Doi et al.13 recently reported a first brief account of the results of some quantitative analysis of the plasticity of the alloys in terms of dislocation theory. Their experiment was rather limited in the composition range covered (volume fraction occupied by the carbide phase: 79 to 83 pct), and thus they could not necessarily elucidate the controlling mechanism of plastic deformation of the alloys of a more general composition range. Consequently, in the present investigation, deformation behavior and some other physical properties of the alloys were investigated and discussed in more detail over a much wider composition range. SPECIMEN PREPARATION WC-Co alloys used in this experiment were prepared in cylindrical or rectangular form by sintering in vacuo compressed mixtures of tungsten carbide and cobalt

Jan 1, 1970

By K. D. Jacob, W. L. Hill

Fluorine recovery in the United States has been restricted chiefly to manufacture of ordinary superphosphate and wet-process phosphoric acid. However, there is an expanding use of fluorine by industry in the form of hydrofluoric acid. Improved methods for recovering fluorine from phosphate rock would yield an adequate supply of fluosilicates and a substantial contribution to the hydrofluoric acid requirement. THE bulk of natural phosphates is comprised of calcium phosphates, which are usually apatites;' calcium aluminum phosphates such as pseudowavel-1ite; and aluminum phosphates, which occur in extensive deposits, notably on Grand Connetable Island." Nearly all deposits of current commercial value are apatite, often with admixed material of other classes. In the trade coarsely crystalline material, mainly fluorapatite of igneous origin, moves as apatite, whereas the fine-grained material of sedimentary origin is marketed as phosphate rock. The latter classification includes the phosphate shales in western North America. Although several genetic varieties of phosphate rock are recognized, most known reserves are of marine origin. Calcium phosphate is concentrated from Florida land pebble and Florida hard rock deposits and from Tennessee brown rock deposits. Phosphate Rock as a Fluorine Carrier: More than two decades ago Reynolds and Jacob surveyed the fluorine content of commercial phosphate rock from various parts of the world. Later these samples were re-analyzed by a new and more reliable method of determining flourine. The revised results with interim additions of new samples, published a few years ago in a brief summary," are shown in Figs. 1-3. Through the points representing the respective varieties of rock medial straight lines are drawn. In some instances the lines are well defined by the points; in others, where scattering renders the position of the best line uncertain, the greater weight is given points near the end of the covered range. Wherever distinct varieties are not recognized among regional deposits the custom is to identify the sample with the geographical location of the deposit. Although all rock varieties cannot be classified with certainty, the plotted results exhibit three distinct types of rock in which the F to P2O5 ratio diminishes, remains sensibly constant, and increases, respectively, as the grade of rock increases. Existence of two types of Florida land pebble, Fig. 1, is supported by additional results recently supplied by a producer." Generally speaking, phosphate rock from continental deposits carries considerably more fluorine, and apatite and insular phosphates less, in proportion to the phosphorus content than is required by the formula for fluorapatite. Exceptions to this are apatite from Canada and certain samples of Florida waste-pond phosphate. The manner in which excess fluorine is held and the nature of the carbonate in carbonate apatite are discussed in recent articles. High-grade calcium phosphates with F to P2O5 ratios less than 0.05 have been found in Curacao, Christmas Island," Mexic,13 and elsewhere. The calcium aluminum phosphates and aluminum phosphates aften carry notable amounts of fluorine. The F to P2O, ratio for one sample of pseudowavel-lite concentrate was 0.05." Aluminum phosphate from Grand Connetable Island is very low in fluorine, whereas some aluminum phosphate minerals" carry as much fluorine as highly fluorinated phosphate rock. Fluorine in Domestic Phosphate: The curves shown in the figures permit interpolation of the F to P2O5 ratio corresponding to the midpoint of selected ranges in grade of any one of the several varieties of phosphate. Multiplication of the interpolated ratio by the mid-range grade (P2O5) yields a figure for the percentage of fluorine in this grade. Then, with the use of the mid-range grade, the amount of fluorine held in phosphate reserves can be estimated, see Table I. Accordingly, on Jan. 1, 1950, the fluorine complement of known domestic reserves exceeded 420,000,000 long tons of the element. A previous figure," deduced in 1940 by a simpler procedure

Jan 1, 1955

By James K. Stanley

PREFERRED orientation, particularly in irons and low-carbon steel, is a phenomenon that is both of considerable importance and theoretical interest. At times it is a liability and at other times an asset. In deep-drawing operations, preferred orientation is extremely undesirable because it necessitates trimming operations and increases scrap loss. Furthermore, oriented metals are likely to crack more easily than randomly oriented material because the forming of symmetrical and unsymmetrical shapes out of preferentially oriented metals and alloys increases the probability of rupture in regions of high stress, owing to differential ductility.1 In the electrical industry, on the other hand, preferentially oriented silicon irons are in demand for transformer application because of the improved magnetic permeabilities in the rolling direction.2 This phenomenon of preferred orientation in soft magnetic materials has been responsible for the tremendous improvement in transformer design and construction in recent years. Directional properties in metals have received a great deal of attention in the metallurgical literature,3-9 in which much of the information reported about orientation has been obtained by conventional methods, mechanical properties and X-rays. A wealth of data has been obtained and often at a great expenditure of time. The recent introduction of the magnetic torque measurements on ferromagnetic materials provides a very simple and rapid method of evaluating orientation. The judicious use of the method can supply important in- formation, but it should be borne in mind that such studies must be supplemented by X-rays, optical methods,10 and etch and pressure figures" to get a rational under- standing of the orientation phenomena. No claims are made that the magnetic torque of evaluating orientation will reveal the behavior of the metal in deep drawing. For this purpose the accepted tests for deep-drawing quality must be used. For a good drawing 'stock, which must have low impurities, desirable grain size, ductility, strength, hardness, good surface, random orientation, and good recrystallization behavior, it cannot be expected that any simple test will be satisfactory for appraising all the factors. The only factor on which the magnetic method can be used is orientation, which is, of course, of considerable importance. When directionality can cause a 10 to 15 per cent change in tensile strength and a IOO to 300 per cent change in tear length,12 it goes without saying that a reliable as well as rapid method for detecting these changes would be desirable. There are certain inherent advantages in testing orientation by the torque method: The method is rapid with the proper type of equipment; no elaborate preparation of specimens is required-it is only necessary punch or machine a disk of the sheet material and insert it into the testing device.

Jan 1, 1943

By L. R. Raymond, J. L. Hudson

This paper proposes a comparatively short-term (8 to 10 hours) well test for detecting and characterizing well-bore damage and for measuring mean formation permeability. The proposed test is made by injecting fluid at constant pressure, recording injection rate as a function of injection time. After one to four hours of injection, the well is shut in and fall-off of bottom-hole pressure is obtained as a function of shut-in time. Formation permeability is estimated by an iterative technique. First, a value of formation permeability is assumed. Then, a plot of the recorded injection rate as a function of dimensionless time is made, using the assumed pertneability value. From the slope of the injection-rate curve. a new value of formation permeability is calculated. If the new value agrees with the original assumed value, the assumption was the correct formation permeability. If the values do not agree, the process is repeated using the new permeability value in the calculation. Convergence is rapid, and a reliable permeability value results. Pressure fall-off data are used to check the result. Graphs of pressure and injection rate us functions of time given in the paper show that changes in permeability of the formation in the neighborhood of the wellbore are disclosed by this technique. Thus, the short-term test can he used to detect formation damage. Also, a rough measure of the radial extent of damage can be inferred, which is helpful in designing stimulation treatments. The mathematical model used for this work was a single-zone, horizontal reservoir with a damaged zone in which permeability decreased continuously as radial distance to the wellbore decreased. This model is more realistic than the usual two-zone, discontinuous permeability model used in published works; calculations indicate the realistic model is valid. Vertical variations in horizontal permeability were studied with this model, and results indicate that the permeability measured by the short-term test is the mean horizontal permeability for the vertical interval tested. The proposed short-term test thus should be useful in detecting and characterizing formation damage and in measuring formation permeability needed in calculating reservoir transmissibility. INTRODUCTION To plan the most efficient production or injection schedule for a well and to design or evaluate the optimal stimu- lation treatment, it is necessary to know the properties of the reservoir adjacent to the well, particularly the reservoir transmissibility and characteristics of a damaged zone, if one exists. Several techniques for determining reservoir transmissi-bility from well tests have been presented in the literature. 1,2,3,4 All these techniques rely on conducting constant-rate well tests that often are difficult to execute. A constant-pressure well test is generally easier to carry Out. and this paper contains the first available method for the analysis of constant-pressure well tests. Determination of wellbore damage from transient well tests has been the subject of several papers."" From these studies it is apparent that information necessary for determination of the characteristics of a damaged zone is available shortly after the transient well test is initiated. Consequently, it may not be necessary to carry out an extensive well test (for example, a pressure build-up test) if the primary purpose of the test is to detect the existence of wellbore damage. All previous studies of well testing to determine wellbore damage have been based on the two-zone perrneability model. In this model the damaged zone has a permeability k,, extending to a radius r,,, and the formation permeability k obtains from r, to the drainage radius r,.. Consequently, there is a discontinuity in permeability at r = r,,. This discontinuity can be eliminated by assuming a continuous variation in permeability through the damaged zone. The effect of this assumption on transient well tests is discussed in following sections of this paper. In addition, all formations have within them vertical permeability variations associated with lithology changes throughout the zone of interest. This paper also considers the effect of these variations on transient well tests. ANALYSIS OF CONSTANT-PRESSURE WELL TESTS The mathematical analysis associated with the injection of fluid at constant wellbore pressure into a single-zone, horizontal reservoir completely filled with a fluid of small and constant compressibility and constant viscosity is given in Appendix A. In this analysis it is assumed that the well is located at the center of an undamaged, circular drainage area. From this analysis, the formation permeability can be obtained as follows. 1, Estimate a value for the formation permeability k. 2. Prepare a plot of injection rate q vs

Jan 1, 1967

By A. M. Sarem, J. M. Campbell

Molecular refraction is introduced as a new and improved third parameter for prediction of the PVT behavior of hydrocarbon systems. This parameter, characterizing the complex as well as the pure hydrocarbon systems (1) is a direct measure of London dispersion forces which affect the PVT behavior; (2) can be measured directly, accurately and easily for complex as well as pure hydrocarbon systems; (3) varies significantly between the reference substances thus allowing accurate interpolation of PVT properties; and (4) has a fundamentally sound basis. Molal average molecular refraction was used with the pseudocritical mixing rule proposed by Leland and Mueller to predict the molal volume of liquid and vapor for 84 published data points on nC4-nCI0, Cl-nC5, Cl-nC4-nC10 systems in the region of coexisting vapor-liquid and at the critical point. The same technique was applied to predict the densities of liquid hydrocarbon mixtures in the two-phase and single-phase regions near the critical point for seven different systems each having different properties of the heptanes-plus fraction, but all simulating the fluids of either gas-condensate or volatile oil reservoirs. A modification of the Leland-Mueller mixing rule for the determination of pseudocritical properties is also presented which yields more accurate liquid density predictions than the original form when used with the proposed third parameter. INTRODUCTION The PVT behavior of multicomponent systems has received widespread attention in the literature in recent years. A large portion of the reported work has stemmed from attempts to overcome the inherent difficulties encountered when using the standard compressibility-factor chart which shows 2 as a function of reduced temperature T, and and reduced pressure p,. This chart assumes that z, is constant for all systems. The error in this assumption is not pronounced for pipeline-type gases high in methane and low in heavy ends, particularly when contained at moderate pressures. For such gases the error in the calculated z is usually less than 2 to 3 per cent. However, the reliability of this approach becomes somewhat nebulous for gas condensate and volatile oil systems at high pressure. One logical approach has been to recognize the variation of 2,. This idea was first suggested by Meissner and Seferianl and used by Lyderson, Greenkorn and Hougen2 to prepare extensive tables of thermodynamic properties. Other third parameters have been proposed by Stockmeyer, Kihara, Riedel, Lightfoot, Bloomer and Pitzer. These were tested by Satter and Campbel13 for many gas mixtures and it was shown that they were inter-related. All of these third parameters have a common fault — they cannot be directly evaluated for the heptanes-plus fraction. Consequently, one necessary requisite of an improved method for predicting the density of multicomponent hydrocarbon systems, was development of a third parameter that could be evaluated for that fraction. PREVIOUS MULTICOMPONENT DENSITY PREDICTION METHODS Sage and Lacey4 suggested a method of using partial molal volumes for computing the density of hydrocarbon gas mixture for pressures to 3,000 psi. In this technique, they approximate the complex multicomponent gaseous mixtures by a quarternary system of methane, ethane, propane and butane. Sage, Hicks and Lacey5 suggest a method of using partial molal volumes for computing the density of hydrocarbon liquids. The results agree within 3 per cent of the experimental values. However, the composition range is limited to about 10 per cent by weight of methane. Consequently, this correlation does not cover the low-molecular-weight liquid similar to natural gasoline and high-pressure gas condensates.

Jan 1, 1966