Search Documents

By G. E. Perry, J. F. Bruskotter, D. G. Russell, J. H. Goodrich

The depletion performance of gas wells has been investigated by mathematical simulation techniques. The gas well model which was studied consisted of a single well located in the center of a bounded, cylindrical, homogeneous reservoir. Dependency of gas compressibility and viscosity on pressure was considered in studies of well performance at both constant mass flow rate and constant flowing pressure conditions. To carry out the investigation, the nonlinear, second-order, partial differential equation which describes Darcy flow of a nonideal gas through porous media was solved numerically. Some of the previous investigations of gas well performance have been of limited general use, because assumptions were introduced to simplify either the gas properties or the basic differential equation. Other studies have been rigorous in these respects but have presented a very limited set of calculated results. The present study was attempted to present a rigorous theoretical model and sufficient numerical results to permit meaningful conclusions to be drawn. It was found that all terms must be retained in the partial differential equation to make accurate predictions. The neglect of higher-order terms, e.g., terms of the order of the "gradient squared", leads to serious material balance errors at large times and to conservative estimates of gas well performance. The higher the gas flow rate and/or the lower the permeability-thickness product of the formation, the more pronounced are these deviations. For example, in a well draining 640 acres in a 25-md-ft formation (8,120 MMcf gas in place) at a constant rate of 993 Mcf/D, the rigorous solution predicts a bottom-hole pressure decline from 4,000 to 1,000 psia in 8.7 years. If higher-order terms are neglected in the differential equation, this decline in pressure is predicted to occur in 5.3 years. With the results of the numerical solution of the differential equation as a basis, simple, easy-to-use approximations for predicting gas well performance for Darcy flow conditions have been developed. These simple approximations are based on the familiar equations for flow of a single, slightly compressible fluid. The approximate methods possess a high degree of accuracy and enable the prediction of long-term gas well performance to be made quickly and accurately without the use of a digital computer. Both transient and stable flow period approximations were developed. INTRODUCTION In recent years income from the sale of natural gas and associated products has represented an ever-incre as ing fraction of the industry's total revenue from operations. To meet the surge in demand for natural gas, the industry has depended heavily upon established reserves and has actively pursued development of new reserves. The search has progressively led to reservoirs which yesterday were too tight and/or deep to yield the desired return on invested capital. More than ever before, evaluation accuracy is now required to forecast the criteria upon which engineering recommendations and management decisions are based. Considerable effort has been expended by both research and operations personnel on the development and application of methods for analyzing and predicting the performance of gas wells. Fundamentally, the problem is the familiar one of extracting data during the drilling, testing and early production life of a well and applying these data within an accurate simulation model to predict long-term behavior. During the past 30 or more years a voluminous literature dealing specifically with gas field problems has been generated. A recent book' lists a comprehensive bibliography of published material through 1959. Over 1,200 references are cited. Since then 39 additional articles on natural gas technology have been published in Transactions volumes of the Society of Petroleum Engineers of AIME. Most existing theory for predicting gas well performance requires that one or more idealizations (e.g., steady-state flow, ideal gas of constant viscosity, small and constant compressibility and constant-viscosity fluid) be applied. Although existing theory may apply directly or be adapted by various artifices to describe specific gas well and reservoir behavior, no widely applicable method is available, and existing methods appear to be subject to appreciable error unless better limits of applicability are defined. The objectives of this paper are (1) to present numerical solutions to the partial differential equation describing gas flow under conditions of general interest in gas well performance prediction work, (2) to present solutions which possess a high order of accuracy for both transient and stabilized flow periods of a well producing at constant rate or constant pressure, and (3) to develop and present simplified relationships which can be used as high-order approximations to the exact numerical results for fast and accurate predictions of gas well performance at the operating level. Combined, these objectives are designed gen-

Jan 1, 1967

By G. S. Cole, H. Biloni, G. F. Bolling

An experimental investigation of sovlze low .melting point alloys sJtows that a substvucture of isolated depressions can be the first manvestation of constitutional supercooling on solid-liquid interjaces veuealed by decanting. Electron-tni cvop vobe and wletallo gvaplic esanzinations, in tlze bulk belzind the interjace, oj the segregation associated with these isolated areas substantiate tlzei&apos;v depressed nature, since a solute of ko < 1 is enriched, and a solute of ko > 1 depleted. In contrast, the pox structuve, a set of projections often veported in the literature, leaves no trace oj. segvegation. These obserl;atims, accovlrpanied by a brief review of recent literature, point to inconsistencies between experirrental obsevvation and the idea that the fornzation of a projection is a causal step in the development of a cellular substructure. An argument is presented to show instead how it is plausible for substantial depvessiom to form in the pvesence of constitutional supercooling at dislocations threading the solid-liquid interjace. THE development of constitutional supercooling during growth from the melt leads to the formation of the cellular solidification substructure. This well-founded association between structure and instability has been basic in understanding cellular substructure and micro segregation; however, the initial formation of structure seems unclear. Rutter and Chalmers,&apos; in definitive experiments and theory, noted that in the presence of constitutional a planar interface might break down: "resulting in the formation of a small projection on an initially plane or uniformly curved interface." That is, the breakdown from a planar to a cellular interface was implied to be initiated via a projection into the unstable liquid. Later, Walton et (11. found that a structure of isolated projections, termed "pox", appeared at solid-liquid interfaces decanted under growth conditions near the onset of constitutional supercooling; the pox were taken as the indication of the instability promoted by the supercooling. Tiller and Rutter4 in their extensive work studied the shape transitions at decanted interfaces which were generally observed to proceed as— pox, "irregular cells", elongated cells, regular (hexagonal) cells, and so forth. The pox varied in size from lo-&apos; to 1CT4 cm, and tended to disappear as cells increased in number and regularity, but as noted,4 the first real array of cells did not seem to be a development from the pox. In fact these authors implied a lack of connection because they stated that the pox are denser on "irregular cells", and as cell boundaries increase in number (i.e., the cells become smaller) there is less need for the pox which do dis- appear. Thereafter, most authors dealing with either experiment or theory have accepted the reality of pox and have used them as a criterion for the onset of constitutional supercooling. In contrast, Spittle, Hunt, and smiths have now suggested that pox are irrelevant artifacts comprised of such things as entrapped oxide. This proposal invokes the observations of weinberg6 and chadwick7 each of whom have shown that the act of decanting leaves a residual liquid on a decanted interface; the remnant solid layer of the order 10 p may thus contain particles that might have been transported from the external surfaces, or elsewhere, during decanting. With the incentive of this suggestion,= some further experiments and a reexamination of the literature have been conducted, in order to question the validity of pox as evidence of an instability and to examine the initial development of the cellular substructure. 1) EXPERIMENTS Single crystals of zone-refined tin (-99.9999 pct) were grown from the melt in a controlled fashion with various, small concentration additions of lead and antimony, for which ko < 1 and > 1, respectively. The crystals were decanted at conditions near the onset of constitutional supercooling and were thus appropriate for observation of slight perturbations. It was possible to observe two types of small departure from smooth or "planar" interfaces in both cases of lead or antimony additions. Some were projections and others, if in regular array of any type, were depressions. The crystals were etched with suitable reagents progressively dissolving the decanted interface surface; projections left no record, but depressions were continuously associated with spotlike areas contrasting with the rest of the interface. Traverses were made with the beam of an electron microprobe across the regions of contrast; with lead addition the persistent spots were lead-rich, and with antimony addition the persistent spots were antimony-poor. This is consistent only with a dominant role for depressions, because if the projections had left spots but were incorrectly catalogued, a reversed observation should have been made; that is, the Pb(ko < 1) should have been depleted and the Sb(ko > 1) enriched. In the work of Cole and inegard, and elewhere, regular arrays of structure associated with the initial stage of instability have been shown, in photographs and represented as pox or projections. We believe this to be erroneous, by inference, since whenever a regular array was observed, in the present examination, it consisted of depressions, regardless of the nature of the solute, ko 1. Fig. 1 is reproduced8 as an ideal example of the possible optical illusion involved; the observer can satisfy himself from the distribution of illuminated areas that the markings are depressions. Fig. 2 from the present investigation is an interference photograph of an interface similar to that in Fig.

Jan 1, 1967

By J. Gurland

The wmk is concerned with the estimate of the degee of continuity of a particulate phase dispersed in a matrix. The first section is a verieu: of the parameters nzeasurable by quantitative rnetallogvaphy, namely the contact area and, under certain conditions, the number of contacts between particles of one phase. In the second section, the frequency of particle contacts is related to the degee of continuity of a randomly embedded phase orz the basis of the lattice-percolation model of Hamrnersley and others. THE extent of continuity of the constituent phases may greatly influence the structure-sensitive properties of phase mixtures, particularly if the phase properties differ significantly. An example from the literature1 is shown in Fig. 1. The marked decrease of the electrical resistivity of sintered silver-alumina bodies near 25 pct by volume of silver presumably occurs because the more conductive phase changes from a discontinuous to a continuous structure. Other example of electrical- and thermal-conductivity transitions are cited by kingery In view of these observations, a knowledge and understanding of the connectivity of the microstructure may be of considerable practical importance. This study is concerned with one aspect of the genera1 problem, namely the estimate of the degree of continuity of a particulate phase dispersed in an opaque matrix. The first section reviews the parameters measurable on a polished cross section by auantitative metallography, such as the contact area and under certain conditions, the number of contacts between particles of one phase. In the second section, an attempt is made to relate the frequency of particle contacts to the degree of continuity of a randomly embedded phase. For the purpose of this discussion, the aggregate has two constituents: the particulate phase and the matrix. The latter is considered to be homogeneous. The term "particle" is used with the meaning of "grain", i.e., a minute unit of matter which is usually a single crystal of a readily distinguishable substance within the confines of a simply connected boundary. While many of the concepts which are used here to describe the structural variables characterizing the microstructure of phase mixtures are qualitative and ill-defined, they are presented with the intent of stimulating further work of a more general and rigorous nature. I) PARTICLE CONTACT, CONTIGUITY, AND MEAN FREE PATH 1) Measurement of Contiguity. The fraction of the total internal surface of a phase a shared with particles of that same phase in a two-phase mixture of and $ is indicated by the contiguity ratio:3 whereis the area of the interface between o particles, and Sap is the area of the interface between a particles and matrix P, per unit volume of mixture. The contiguity ratio ranges from 0 to 1 as the distribution of the a phase changes from a completely dispersed to a completely agglomerated structure. The interfacial areas may be obtained from simple intercept measurements on a random plane, after Smith and Guttman where and N are the numbers of intercepts between a random line of unit length on a plane of polish and the interparticle and interphase interfaces, respectively. The equations are valid for any distribution of particle sizes and shapes.

Jan 1, 1967

By P. H. Thornton

MaNY intermediate phases are solid solutions which have a solute-solubility range greater than 3 at. pct. These phases, which are generally very strong and often have simple crystal structures, are usually extremely brittle so that specimens of simple shape cannot be produced by the usual mechanical-deformation techniques. Furthermore, single crystals of reasonably uniform composition are difficult to grow by the Bridgman technique because the phases generally form by a peritectic reaction, which causes pronounced segregation. In this note, a technique is described for the preparation of thin-rod specimens of brittle intermediate phases, which may be used for physical-property studies such as resistivity changes&apos; and for the conversion of such rods into single crystals. The method has been used for the preparation of rod single crystals of the phases and Ingots, of about 100 g weight and 11 mm diam, were made by melting together the appropriate amounts of high-purity materials (usually 99.999 pct) in sealed

Jan 1, 1967

By M. S. Burton, B. W. Christ, G. V. Smith

Evidence of the effects which high-angle grain boundaries and veining subgrain boundaries have on the yield phenomenon of zone-refined iron is reported. Room-temperature tensile tests on 50-mil wive at several gram sizes produced upon recrystallization at diffevent temperatures below 910C (a-y phase-transfornzatzon temperature) show that upper and lower yield stress decrease with increasing grain size. The upper yield point effectively disappears at grain sizes of 81 and 102 u in samples of this diameter. However, samples of 129 p grain size. produced during recrys-tallization at 950°C, exhibit a pronounced upper yield point and an increased lower yield stress. Possible causes of the reversal in the grain-size dependence of yielding are discussed. Veining substructure, detected with picval-HC1 etchant, but not with 2 pct nital, is formed during the a-y transformation, and is believed to account for the observed effecls. The strengthening associated with high-angle boundaries in zone-refined tron is discussed in terms of: 1) the stress dependence of slzp multiplicity, and 2) the stress dependence of dislocation emission from grain boundaries. WIDE use of transmission electron microscopy and dislocation etch-pitting techniques has led to a growing awareness of the presence of dislocation substructures in annealed pure metals. Whereas much speculation&apos;-4 has developed about the influence which subgrain structures have on the mechanical behavior of metals, particularly the yield strength, only limited experimental evidence correlating subgrain structures with lower yield strength has been presented. 2,3,5-7 It is the purpose of this note to present experimental evidence of the effects which grain boundaries and subgrain boundaries have on both the upper and lower yield strength of zone-refined iron. Effects due to high-angle boundaries formed upon recrystallization below the a-y phase transformation are compared with effects due to veining subgrain boundaries formed upon "recrystallization" above the a-y phase-transformation temperature. These effects are further compared with the reported influence 2,3,5-7 which low-angle boundaries have on the lower yield strength of bcc transition metals. EXPERIMENTAL PROCEDURE The starting material for these experiments was a 1-in.-diam zone-refined ingot prepared at Battelle Memorial Institute. Initial purification of electrolytic iron by vacuum-induction melting was improved by a floating zone-melting treatment in ultrapure hydrogen.&apos; Chemical analysis of this zone-refined iron appears in Table I. Tensile samples of 50-mil wire were prepared by rolling and drawing portions of the ingot. An intermediate recrystallization anneal at 625°C after rolling to a reduction in area of 85 pct was carried out in an evacuated Vycor capsule (10-=5 mm Hg). Further reduction in area of 91 pct was followed by recrystallization in evacuated Vycor capsules at temperatures between 600" and 950°C, see Table 111. All samples

Jan 1, 1967

By S. E. Szasz, E. J. Moore, B. F. Whitney

An accurate value of formation water resistivity R, is essential in calculating formation porosity and fluid saturation from electrical well logs. In the cases where R, has not been measured directly, it must be obtained from other data, e.g., the SP curve. This paper deals with another approach: how to calculate R, from the chemical analysis of the formation water. INTRODUCTION It is known that the resistivity of aqueous solutions of pure salts depends on their concentration and on the temperature; the concentrations are given in MPL (mg of solute per liter of solution), or sometimes in ppm (mg of solute per kg of solution): MPL = ppm X specific gravity. Values for different pure salts are available in the literature, but not for solutions of mixtures which are of practical interest. The major component of the dissolved material in almost all formation waters being sodium chloride, it is customary to express the resistivity of formation waters in terms of equivalent sodium chloride concentration, i.e., the concentration of a solution of pure NaCl which has the same resistivity at a given temperature as that of the formation water under consideration. Thus, the problem of calculating R, from the chemical anaylsis can also be stated as how to convert the other constituents of the solute into equivalent NaCl concentration. Salts dissolved in water are at least partly dissociated into ions, and do not conserve their identity. If known amounts of several salts are dissolved in water, the solution does not necessarily contain the same salts in the original proportion, but perhaps some other combination of the ions, along with free ions in solution. This is why the chemical analysis of formation waters is often given in terms of ions, as if all dissolved salts were completely dissociated. Our problem then boils down to how to convert the concentrations of the various ions to equivalent concentrations of Na&apos; and C1-. Dunlap and Hawthorne&apos; have proposed to convert the concentration of all other ions to equivalent Na&apos; and C1-concentrations by means of constant multipliers; e.g., 0.95 for Ca"; 2.0 for Mg"; 0.27 for HCO 3-; 0.5 for SO, -, etc. Their factors were based on measurements made at 68F on 26 formation water samples from the Texas Gulf Coast, ranging in concentration from 1,500 to 75,000 ppm. The Dunlap method is widely used in electric log interpretation, and is often extrapolated beyond its original concentration range. A comparison of R, values calculated by this method and values actually measured on formation water samples has shown large discrepancies, especially at higher concentrations. Therefore, two new methods were developed at Sinclair Oil Corp.&apos;s Tulsa Research Center to calculate equivalent sodium chloride concentration from the chemical analysis of formation water samples. FUNDAMENTAL CONSIDERATIONS The resistivity of a solution, or its reciprocal the conductivity, at a given temperature is determined by the charge, concentration and mobility of the ions actually present. Monovalent ions such as Na&apos; or C1- always carry the same charge. Compounds of polyvalent ions, however. may show incomplete dissociation, e.g., NaSO; + Na&apos; instead of SO,-- + 2Na&apos;. This happens especially in more concentrated solutions. Only very dilute solutions are completely dissociated, as assumed in the chemical analysis report. At higher concentration, the degree of dissociation depends not only on the nature and concentration of the particular salt under consideration but also on the nature and concentrations of the other solutes. Mobility of the ions depends on the viscosity of the solution. It also depends on the degree of hydration of the ions, which in turn is a function of the nature and the charge of the ions and also of the amount of free water available per ion, i.e., the total ionic concentration. The net effect is that the conductivity increases slower than proportional to the concentration, even if a solution contains only one salt such as NaC1, and is different for different salts (Fig. 1). Conductivity can even decline with a further increase in concentration, e.g., if additional salt is little dissociated but ties up some of the free water and/or causes an increase in viscosity. In solutions containing more than one salt, the contribution of one salt to the total conductivity depends not only on the fractional concentration of this same salt, but also on the concentration of all other solutes. A perfect method would give the conductivity or resistivity of a solution as a function of the concentrations of all solutes present. This is so complicated as to be impractical, and a simpler method must be found which is of acceptable accuracy. The Dunlap method, on the other hand. is too simple because it askmes that at any concentration the contrih-

Jan 1, 1967

By B. C. Allen

Grain-boundary gvoocing was studied irz chronziu?n. molybdenum, tungsten, and the solid-solution alloys, Cr-35Re. Mo-33Re, and W-25Re at 0.6 to 0.9 of the absolute liquidus temperature under an inevt atmosphere , foveign nzetal zlapor, and liquid metal. The controlling mechanism of matter transport was deduced frolr groove size and shape. kinetics of growth, and cotnparisoz with theoretical relatiotzships derived by Mullins. Grooves in chowzizon formed prittzarily by volume diffusion tkrongk kelium or argon at 1 atm pressure and by e1)apovation-condensatiotz in a static vacuum. Groozles in the other nzetals and alloys forrned priniarily by sirface diffusion, the coefficients of which were calculated. Carbon was found to decrease the acti7:ation energy for surface diffusion on molybdenun. Sitrface diffusizities weve sinzilav for iWo-0.0ZC and Mo-33Re. The presence of about 1 -?nm pressure of silver or tin vapor did not affect the neckanis)n of grooz;e formation. In liquid titz or silz:er, pain bolindary gvooves formed by zolume diffusion in the liquid. BASED on theories of ullins," the mechanisms of grain-boundary grooving have been determined for low-melting metals such as copper in an inert atmospheres or in liquid lead.4 The purpose of this investigation was to study grain-boundary grooving in the refractory Group VI-A metals, chromium, molybdenum, and tungsten, and the bcc solid-solution alloys, r-35e, 0-33e,= and -25e,&apos; at 0.6 to 0.9 of included to see whether Mullins&apos; theories were applicable to them as well. EXPERIMENTAL WORK The materials chromium, molybdenum, tungsten, Cr-35Re, Mo-33Re, and W-25Re were in the form of 0.050-cm-thick sheet assaying 99.98 pct or better. The nominal carbon level was generally 30 ppm. The oxygen + nitrogen + hydrogen content was -50 ppm. Mo-33Re and W-25Re were prepared by powder-metallurgy methods and supplied by Chase Brass and Copper Co. The remaining materials were arc-cast. Chromium and Cr-35Re were prepared at Battelle by arc melting, extruding, swaging, and rolling. Mo-0.003C was supplied by Climax Molybdenum Co., and Mo-0.02C and tungsten by Universal Cyclops Corp. As a source of liquid metals, tin and silver granules assaying 99.98 pct were used. Spectrographic analyses indicated the presence of -100 ppm each of iron and silicon. Vacuum melting reduced the combined gas content to c5 ppm (2.60,, 0.9N2, 0.4H2 in tin and l.602, 0.7N2, and c0.04H2 in silver). Specimens were prepared for grain-boundary grooving. Coupons 0.5 by 1.5 cm with a small hole in one end, or substrates 1.5 by 1.5 cm, were cut from sheet. The specimens were annealed in containers of the same material to provide a vapor-solid equilibrium as well as to protect them from gaseous impurities. Outgassing was done at 2 x 10"5 mm for 15 to 30 min at 1150°C for chromium and Cr-35Re, 1750°C for molybdenum and Mo-33Re, and 2600°C for tungsten and W-25Re. Then recrystallization to a stable grain size of about 0.5 mm was accomplished by annealing at least 1 hr at 1600°C for chromium, 1800°C for Cr-35Re, 2350°C for molybdenum, and 2300°C for Mo-33Re in a rhenium-element resistance furnace under 1 atm of static argon or helium. The atmosphere was gettered by tantalum radiation shielding. Tungsten and

Jan 1, 1967

By K. Mukherjee

A quasi-thermodynamic calculation is presented for estirt~ating vacancy Parameters in fcc metals. For fcc metals the equilibrium vacancy concentration at the melting point can be approximated by a constant mole fraction. The Gibbs free energy of formation of the vacancies can then be expressed as a simple function of the absolute melting temperature alone. From this description the formation energy Ef, the vibrational entropy Sv, and the formation volume of a vacancy can be calculated in terms of known thermodynamic parameters. The calculated values are in good agree-rnent with available experimental results. IT appears that many of the physical parameters associated with lattice defects can be empirically related to thermodynamic melting parameters.&apos;-7 Besides the implicit correlations between melting phenomenon and point defects, the melting temperature itself at atmospheric pressure can be used as a reference temperature for the purpose of comparison of vacancy concentrations in a given crystal class. The concentration at the melting point then corresponds to the maximum vacancy super saturation in the crystalline state at a fixed pressure. For pure fcc metals such as Au, Cu, Al, Ag, Pb, and Pt the mole fraction of vacancies at the melting point8 is of the order of 10*"3 to 10"1. The equation expressing the equilibrium vacancy concentration at the melting point contains three parameters: the formation energy, the vibrational entropy, and the absolute melting temperature, see Eq. [2]. Since these three parameters cover a wide range of values8 in the above-mentioned metals, it can be assumed that in other fcc metals, for which no reliable data are available, the mole fraction of vacancies at the melting point most likely lies in the range 10~3 to lO"4. This is then the basis of the calculations described in this article. Within this range the mixing entropy can be linearized and the vibrational entropy, the formation energy, and the vacancy formation volume can be calculated in terms of known thermodynamic constants. 1} LINEARIZATION OF MIXING ENTROPY The Helmholtz free energy of a crystal containing n vacancies at the melting point can be writteng as where Ef, T,, k, and S, are the formation energy of a vacancy, the absolute melting temperature, Boltz-mann&apos;s constant, and the vacancy formation entropy, respectively. After applying Stirling&apos;s approximation in Eq. [l] and minimizing with respect to n where (n/N), is the equilibrium vacancy concentration in mole fraction at the melting temperature at 1 atm pressure. If instead of setting (6~/6n) = 0 immediately after Stirling&apos;s approximation in Eq. [I], the configurational partition function is simplified, then some useful and significant results can be obtained by utilizing the experimental inforrnation After applying Stirling&apos;s approximation and rearranging is plotted as a function of (n/N) in Fig. 1. The points, not all of which are shown, are fitted by a straight line substituting from Eq. [5] into Eq. [I] and minimizing with respect to n combining Eqs. [6] and [2] It must be emphasized that Eq. [7] is only a crude approximation based on an oversimplified model of

Jan 1, 1967

By P. L. Rittenhouse, M. L. Picklesimer

A rapid and seniquantitative method of determining prefered orientation on large numbers of. Zircaloy-2 specimens was desired. knoop microhardness measurerrzetzls were irvestigated as a solldtion to this pro6lem. The variation of Knoop microhardness measurerlerzts on selected planes as a function of &apos;irzdenle axis relutilje lo cryslallograplzic or fabrication divectals were used to corzstrrcl a polar coordinate hardness contour map. With use of an empirical relationslip between the single-crystal hardnesses and those of the polycrystalline material conentional pole figures could be constructed which compare favorably also obtained. To determine preferred oriention qualitatively from hardness data requires a minimum of twelve measurements per plane on three, preferattention to grain size, specimen prepartion, and in-rreinetzt, and analjlsis is of the order of- 45 to 60 rrin. THE design of structures from Zr alloys requires consideration of preferred orientation and the resulting anisotropy of mechanical properties. Rapid and semi-quantitative methods of evaluating anisotropy and of determining preferred orientation are needed for quality control and for examining the large number of test specimens required in development programs. The variation of Knoop microhardness with indenter orientation has been studied in single crystals of several hcp metals including titanium.&apos; beryllium,&apos; magneium, and zinc.4 The maximum hardness observed on any crystallographic plane other than the basal plane invariably occurred when the long axis of the Knoop indenter was either parallel or perpendicular to the projection of the [ OOOI.] direction on the plane of examination. When the major slip mode was (0001)(i2i0j the maximum hardness occurred at the parallel positio! but when the operating slip system was {10i0) (i210) the hardness was a maximum perpendicular to the [0001] projection. A minimum hardness was always observed at a rotation of 90 deg from the maximum. It follows that the orientation of the projection of the [0001] can be determined on any non-basal plane by making hardness measurements at a number of indenter orientations. The determination of the orientation of the [0001] projection on two non-parallel, preferably orthogonal, surfaces of the specimen will allow location of the [0001] direction in the specimen. It seems possible that such measurements could be used to examine, at least qualitatively, the preferred orientation existing in polycrystalline hcp materials. An investigation of the microhardness anisotropy in Zircaloy-2 was undertaken to ascertain whether these measurements could be used for this task. EXPERIMENTAL PROCEDURE Single crystals of Zircaloy-2 were grown by an a-8-a annealing sequence using electron-beam heating.= The orientation of the crystals was determined by a conventional back-reflection Laue technique. The crystals were then mounted on a goniometer head and the desired crystallographic faces were milled and chemically polished. Polycrystalline Zircaloy-2 specimens were prepared from fabricated sheets or plates. Inverse pole figures for these materials were obtained using the X-ray diffraction techniaue described by Jetter, McHargue, and illiams.&apos; A Wolpert-Greis Micro-Reflex hardness-testing machine was used to make the Knoop microhardness measurements. The single-crystal and polycrystalline specimens were loaded to 0.5 and 2.0 kg, respectively. Determination of the scatter of hardness numbers as a function of applied load for several indenter orientations on several specimen surfaces showed that the loads selected were the lightest consistent with minimum scatter. Heavier loads did not appreciably decrease the scatter and produced hardness impressions too large to be conveniently measured with the equipment used. Seven crystallographic planes of the single crystals were examined, while six planes of examination were used in studying Zircaloy-2 polycrystals. The specimens were rotated 10 to 15 deg after each measurement and from four to eight impressions were made at each angle of rotation. RESULTS The two angles which were used to relate crystallographic directions and planes in the hcp cell to the planes of examination of the single crystals are shown in Fig. 1. 8 is the angle between the c axis. [0001]. and the normal, N, to the plane of examination. a is the angle between the long diagonal of the Knoop indenter and the projection of the c axis on the plane of examination. For exampIe, if P = 90 deg and a = 0 deg, the plane of examination is of the family

Jan 1, 1967

By J. H. Brophy, A. I. Prill

The authors have previously investigated the solid-state activated sintering of tungsten powder with various Ni-Cu additions.&apos; Log shrinkage vs log time plots showed changes in slope from 4 to x at intermediate times during the sintering process. This was attributed to a change in the rate-controlling mass-transport process from solution of tungsten into the tungsten-alloy interface to diffusion of tungsten in this interface. It was postulated that the melting of the activating layer should not affect this transport process. However, the superposition of rearrangement due to the presence of a liquid phase complicates the data analysis. It is the purpose of this note to apply the analysis of Prill, Hayden, and rohy&apos; to account for "rearrangement shrinkage&apos;? in liquid-phase sintering of the W-Ni-Cu system. This results in "solution-precipitation stage" shrinkage plots which were comparable to those obtained in solid-phase activated sintering of the same system. If one assumes total shrinkage to be the sum of rearrangement and solution precipitation, the effects can be isolated easily.2 This results in the general equation:

Jan 1, 1967

By W. V. Long, D. H. Reynolds, R. B. Bhappu

Equilibrium distribution coefficients were determined for some monovalent and polyvalent ions for several absorbing and exchanging media as functions of pH and initial ion concentration. Among ions studied were NH4 (I), K(I), Ni(II), Ce(III), Th(IV), and Mo(VI). Absorbents used were: ion exchange resins; a tertiary, and a quaternary amine in an inert solvent; and some activated carbons. The ions studied were in aqueous HCI, HNO3, H2SO4, or HC104 solutions. The classical mass action expression proved of little value in describing the variations observed. The data indicate that the determining factors in absorption and ion exchange may be the ionic species in the aqueous solution, which in turn depend upon pH, dilution, and perhaps other conditions. Three sorption processes have become prominent in recovering values from leach solutions. These are: 1) ion exchange, 2) solvent extraction (liquid ion exchange), and 3) physical adsorption by materials of high specific surface. In the first two of these processes the exchange reaction is considered chemical in nature, and is discussed in terms of the mass action principle. The third is considered a surface phenomenon and is discussed in terms of surface areas and various forms (Gouy, Stern, Grahame, etc.) of the Helmholtz double-layer theory. Generally speaking, we have found in chemical analysis and hydrometallurgical studies that there are close similarities between the three sorption processes, and also differences which are difficult to reconcile with any of the usual theoretical treatments. That others have experienced the same troubles is shown by the fact that usable processes, in chemical analysis as well as in hydrometallurgy, generally are discovered by empirical procedures. EXPERIMENTAL CONSlDERATlONS The gaps between theory and practice have allowed many conflicting statements to appear in the litera- ture. To learn, in particular, the behavior of some metallic ions on a cation exchange resin at a constant pH, equilibration studies were made over a wide range of concentrations.&apos; Results are shown in Fig. 1 for four metals in hydrochloric acid solution on a strongly acidic cation exchange resin, AG 50W-X8 (Bio-Rad). The distribution coefficient, defined as D = amount of metal extracted per g of resin/amount of metal unextracted per ml of solution is seen to vary in a way to support all manner of reports about its behavior with changes in solution concentration. The metals used were chosen from those least likely to form stable chloro-complexes. The changes of slope in the various isotherms may reflect formation of unsuspected chloro-complexes of low stability (ion-pairing for example) whose concen-

Jan 1, 1967

Jan 1, 1967

By Raynor Linzey, Karl M. Busen

The formation of thin films of silicon oxide by anodic oxidation of silicon and the subsequent removal of these films by an etch is a process often used for the evaluation of concentration distributions Profiles) in silicon layers by the differential sheet conductance method. The accuracy of the resulting profile is very strongly influenced by the uniformity of the thickness of the reacted silicon. Normally, it would be expected that for a constant number of coulombs passed the thickness would be the same from oxidation to oxidation. Investigations show that in certain electrolytes, for a given number of coulombs passed through an n-type silicon sample, the thickness of the reacted silicon increased with increasing resistivity. Even for the same resistivity the thickness varied sometimes by a factor of 1.5. When an electrolyte was used which consisted of 10 pct water by volume in ethylene glycol with 4.0 g KNO per 1000 ml of solution, anodiza-tion at 5 ma per sq cm led to satisfying results. Short-time anodizations gave oxide layers of a higher apparent density than those experienced from thermally gown silicon oxides. THE functioning and the electrical characteristics of semiconductor devices are based upon the incorporation of "impurities" into a single-crystalline body of suitable material and on the concentration distribution (profile) of this impurity. The incorporation can be achieved by well-known processes as, for example, by diffusion, epitaxial growth, alloying, or ion implantation. Often the profile resulting from these processes is not known. A powerful tool to learn about a concentration distribution is given by the method of differential sheet conductance which employs successive four-probe measurements on a layer subjected to stepwise removal of thin sublayers. Differential sheet conductance vs position (or penetration depth) of a sublayer is plotted and a smooth curve is drawn through the data. From this curve the profile is then calculated. When the semiconductor material is silicon, the sublayers are removed suitably by anodic oxidation of the silicon and subsequent dissolving of the formed silicon oxide by an etch. The accuracy of the resulting profile is very strongly influenced by the uniformity of the sublayer thickness. Normally, it would be expected that, for a constant number of coulombs passed, the thickness would be constant from oxidation to oxidation. Investigations showed that for sublayers several hundred angstroms thick the reproducibility can be rather poor. Therefore, efforts were made to obtain a reliable technique for uniform removal. The present paper describes such a technique and certain phenomena which were encountered during the investigations. APPARATUS AND TECHNIQUE The first report on the investigation of concentration distributions, where for differential sheet conductance measurements thin silicon sublayers were removed from a diffused layer by anodic oxidation, was given by Tannenbaum. The author reports the removal of sublayers which for the most part were 400 thick. More advanced device designs now require much narrower layers. When it was tried in these laboratories to determine profiles within such layers, difficulties were encountered with respect to sublayers which by necessity had to be thinner than the ones reported by Tannenbaum. The apparatus used for the investigations is sketched in Fig. 1. Two cylindrical containers connected by a wide tube are filled with an electrolyte. The left container receives the cathode whereas the right container is closed at the bottom end by a sample support consisting of a Teflon base and a copper pedestal. The silicon sample (1 by 1 cm) is mounted to the pedestal embedded flush in the Teflon base using silver paint (Degussa) for electrical contact. Pyseal is applied to the edges of the sample to protect the copper pedestal from the electrolyte. The base is mounted to the container using a water-tight silicon rubber gasket. Fig. 2 gives a view of the sample support. The copper pedestal is connected to the positive side of a power supply operating at a constant current output (Kepco Model ABC 425M). The electrolyte which has been reported by Duffek et al. consisted of either 2 or 10 pct water by volume in ethylene glycol and 4.0 g KNO3 in 1000-ml solution. The electrolyte is best prepared

Jan 1, 1967

By R. N. Upadhyay, M. Maleki, I. Fatt

Conventional laboratory core analysis tests on samples of two limestone reservoir rocks indicate that about 20 per cent of PV is in dead-end pores. These tests (electric logging formation factor, mercury injection capillary pressure and miscible displacement) were carried out on 1/4-in. diameter test plugs. Test results show a clear dilference between these samples and sandstone or homogeneous limestone reservoir rock. Although the amount of dead-end pore space can be only roughly estimated, the presence of such pore space seems clearly indicated. Pressure transient studies also show presence of deadend PV. Although they do not give quantitative results, pressure transient data yield a reasonable estimate of the size of the neck connecting dead-end pores to the main flow channels. INTRODUCTION Equations conventionally used to describe reservoir flow behavior contain the implicit assumption that all connected pore spaces contributed to both porosity and permeability. Several authors3-7 have pointed out the changes in pressure transient behavior and in electric log interpretation that may result if this assumption is incorrect and, instead, dead-end or cul-de-sac pores are present. There is a need for laboratory tests that can detect presence of dead-end pores in core samples. With such information on hand the petroleum engineer can make more rational use of the mathematical tools now available for analysis of reservoir flow behavior. This paper describes laboratory studies designed to detect and, if possible, give a quantitative measure of dead-end PV in laboratory-size core plugs. Three reservoir rocks were used, two of which were limestones suspected of having deadend pore spaces and a well-known sandstone, used as a comparison standard, in which there is believed to be little or no dead-end pore space. All the studies were designed to measure the natural dead-end PV; i.e., the pore space which is dead-ended because of rock structure. During multiphase flow in a rock without dead-end pores, some parts of one of the phases can become surrounded by the other, thereby giving (for certain flow behavior) an effective dead-end pV8,9. Such behavior will not be described here.

Jan 1, 1967

By A. T. English

X-ray dif action techniques were used to investigate ordering and domain coalescence for the B2 (CsC1) superlattice in this alloy. Essentially complete disorder is obtainable by drastic quenching. isothertmal ordering occurs very rapidly at temperatures between 565° and 680°C, only a few seconds being required for development of substantial long-range order. A very fine antiphase dormain structure, found in samples annealed for short times, undergoes coalescence with further annealing. Domain "size" is found proportional to (t)1/2 . This Yate law can be derived by adapting the formalism of Hillert&apos;s analysis 01- metallurgical grain growth. For this Purpose the material is descvibed in terms of a "Swiss-cheese " structure, consisting of two multiply connected, mutually intertwined domain volumes. The observed activation enevgy for coalescence, 70 kcal per g-atom, agrees with an estit~lated value ,for volume diffusion in disordered rrzaterial. It is suggesled that the experimental value vejlects the substantially disordered environment of mobile vacancies located in the moving antiphase domain boundary. The alloy Fe-Co-2V is bcc below about 900&apos;C and orders into a B2 (CsC1) superlattice below about 730c It has important applications in magnetic components, including many requiring extensive cold rolling.* example, Marcinkowski and Brown6 observed a re1atively small as-quenched domain size of about 250Å in B2 Fe3A1. Since there appear to be no published data on domain coalescence in B2 alloys, this process was studied in detail. MATERIALS AND EXPERIMENTAL TECHNIQUES A commercial high-purity alloy, nominally 49 Fe-49Co-2V, was obtained from the Arnold Engineering Co. in the form of cold-rolled sheet 0.004 in. thick. Samples were annealed for 1 hr at 800°C in dry H2 to produce disorder and quenched into iced brine. They were then annealed in a nitrate-nitrite salt bath for various times at temperatures between 525° and 625°C. Estimates of the LRO parameter S and the APD size D were made using established X-ray diffraction techniques. Superlattice peaks are best observed using CoKa radiation, for which the value of (fFe-fCo) is maximized. Peak-to-background intensity ratio was improved by filtering CoKß from the incident beam to eliminate sample fluorescence, and by using a pulse-height analyzer. The LRO parameter S was computed from the ratio of integrated intensities I 100/I 200- Owing presumably to systematic errors associated with the presence of crystallographic texture and X-ray extinction, only rough agreement was obtained between experimental and theoretical (absolute) intensity ratios. For this reason certain experimental refinements (e.g., use of balanced filters) were considered superfluous and were omitted. Emphasis was placed on maximizing intensity in order to reduce the time taken to accurately determine the contour of the very small (loo), peak. This was done by point counting 105 counts at appropriate intervals. Background intensity averaged about 1200 cps.

Jan 1, 1967

By M. J. Mullikin, J. B. Conway

for use in describing combined first- and second-stage creep data. In the above ?o represents the instantaneous strain on loading, ?t the limit of transient creep, ?S the steady-state creep rate, t the time, and ? a constant leading to the exhaustion rate of transient creep. It was noted that this equation was similar in form to that proposed originally2 by McVetty. This type of equation was also suggested by Andrade and Aboav3 in studies of the shear strain of cadmium. Calculation of the constants ?o, ?t, ?, and ?s, in Eq. [ l] is somewhat involved due to the slight complexity of this relationship. In the paper by Garofalo a computer approach based on a least-squares analysis was recommended. Such an approach involves a trial and error solution but when a computer facility is available the difficulties in this calculation are of little consequence. It does mean however that the evaluation of this type of equation would become quite tedious if computer facilities were not available. Some recent studies have identified a fairly straightforward method for calculating the constants in the Garofalo equation. While not leading to a least-squares evaluation the approach is accurate and is recom-

Jan 1, 1967

By M. Sagheer

This paper gives a report of the laboratory work investigating the fundamental flotation characteristics of chromite and its commonly associated mineral, serpentine. Floatability curves, showing floatability of chromite at different concentrations of sodium oleate and pH, were obtained by pick-up method. One of the difficulties encountered in flotation of chromite in presence of serpentine has been attributed to dissolved metal ions. The influence of metal ions, dissolved in pulp filtrate, on floatability of pure chromite and serpentine has been studied and it has been shown that Mg, Ca and Fe ions depress chromite. A considerable amount of work applying bubble pick-up method was carried out to counteract the depressing action of metal ions. Various methods and reagents were applied in this connection. Validity of each method and reagent has been discussed. Probably due to the abundance of high-grade chrome ores, very little has been published concerning the flotation of chromite. As improvements in mineral beneficiation are realized, more attention is being directed to low-grade ores previously considered unusable. Currently, gravity concentration is the principal process applied to upgrade chrome ores, but this cannot be used for finely disseminated ore because of the substantial losses in the fines,1,2,3 To avoid these losses, flotation can be applied to material which is too fine for satisfactory gravity concent ration, and it is probable that a combination

Jan 1, 1967

By John F. Elliott, David L. Schroeder

SURFACE depletion can be a significant source of error in the use of the effusion method for determining the thermodynamic properties of alloys, particularly in the case of solid alloys.&apos; However, as yet this problem has not been treated quantitatively by experimentalists in the field.1-3 The method necessitates diffusive transport of material out of the source alloy to obtain a relevant measurement. The general character of the system is shown in Fig. 1. The diffusive flux (nD) at the alloy surface is equal to the effusive flux (nE): where D is the diffusion coefficient of the diffusing substance, C is its concentration, X is the distance from the surface into the source alloy, p is the density of the alloy, and AS is the actual available surface of the alloy. The concentration of the diffusing substance at the surface is CS, and Here fE is the Clausing factor, AE is the area (sq cm) of the orifice, y is the activity coefficient of the substance in the alloy, PO is the vapor pressure (Torr) of the pure substance, M is the molecular weight of the effusing species, and T is the absolute temperature. The differential equation for diffusive transport in the source alloy is Fick&apos;s second law. Its solution is dependent on the geometry of the source alloy. Two boundary conditions must be established to obtain the proper solution, the first being Eq. [1], or its equivalent. The solutions for several important forms of the source alloy are considered below.

Jan 1, 1967

By P. Marier, T. R. Ingraham

When cylindrical powder compacts of nickelous oxide (NiO) are heated in an equilibrated stream of SO3, SO,, and 02, in the temperature range 500°to 800°c, tlickelous suljate (NiSO4) is formed as an adherent layer on the oxide. The reaction is diffusion-controlled, with an activation energy of 13 ± 3 kcal per mole. When cylindrical powder compacts of NiSO4 are heated in a stream of dry nitrogen in the temperature range 1000 0to 11500C, NiO is formed at a well-defined interface between the sulfate and oxide. The interjhce migrates into the sulfate at a constant rate at a constant temperature and gas flow. The reaction rate increases with increasing gas flow and decreases with increasing partial pressure of SO3 in the gas stream. Resolution of the partial-pressure data with a Lang-itzliir Adsorption Isothern suggests that about 80 pct of the NiO interface is covered with SO3. The heat of adsovption is about 68 kcal per mole. The activation energy of the decomposition process is 61 ± 3 kcal per mole. MOST of the work reported on the formation and decomposition of nickel sulfate is related to thermodynamics.1-3 No quantitative work has been reported on the decomposition process. The only kinetic work reported on the formation is that by Fletcher and shelef.4-6 whose main purpose was to study the effects of alkali metal sulfate additions on increasing the rate of nickel sulfate formation. They found that, in the absence of an additive, only small amounts of sulfate were formed and the kinetic data were difficult to interpret. However, they suggested that the sulfa-tion of nickel oxide is a diffusion-controlled process for which the activation energy is approximately 17 kcal per mole.6 In this paper, both the formation and decomposition reactions will be examined; and based on the effects of temperature, flow rate, interfacial area, and partial-pressure of sulfur trioxide, mechanisms and activation energies will be suggested for both processes. EXPERIMENTAL Material. The nickel source material used in all experiments was Baker and Adamson Reagent Grade NiSO4. 6H2O, for which the following analysis was supplied by the manufacturer: C1, 0.001 pct; NO3, 0.005 pct; Co, 0.02 pct; Cu, 0.005 pct; Pb, 0.002 pct; Fe, 0.001 pct: insol, 0.005 pct. The anhydrous sulfate was prepared by heating the hydrate at 420°C in a muffle furnace for 5 hr. Nickel oxide was prepared by decomposing the sulfate at 950°C over a period of 10 hr. Method. The technique used for measuring the rates of formation and decomposition was developed earlier in these laboratories.7-8 The method is based on the observation that the kinetic parameters of a reaction can be established in a flow system, with specimens prepared by compacting powder to a uniform pellet shape. The results are identical with those established with solid polycrystalline specimens or single

Jan 1, 1967

By H. D. Colson, Robert A. Rapp

When a binary alloy is oxidized to form simultaneously an internal oxidation zone and an external scale which grows at a constant rate, the oxidation should reach a condition of steady shale, in which the internal oxidation zone is maintained at a time-independent thickness. Assuming an ideal model for the steady-state oxidation, the diffusion equation is solved for the steudy-state concentration profiles for oxygen and alloying- element in the metal. A theoretical prediction of&apos; the steady-state thickness of the internal oxidation zone is made in terms of the experimental parameters. FOR the oxidation of binary alloys, Rhines, Johnson, and Anderson,1 as well as Mask: have derived equations describing the kinetics for the simultaneous formation of an internal oxidation zone and an external scale which was assumed to mow according to a parabolic rate law. Experimental determinations of internal oxidation zone thicknesses beneath a Cu2O scale for Cu-Be and Cu-A1 alloys3,4 have shown some deviation from Maak&apos;s theoretical expressions, because the ideal model assumed for the calculation is not met experimentally. Some pure metals, such as columbium and tantalum oxidize at elevated temperatures and high oxygen pressures to form nonprotec-tive external scales which thicken at a constant rate. If such a base metal is alloyed with a solute element which may form a more stable oxide than the solvent oxide, then internal oxidation of solute may occur beneath the growing external scale. A review of the internal oxidation of alloys has been published recently.5 The purpose of this paper is to analyze theoretically for a binary alloy the kinetics for the simultaneous formation of an internal oxidation zone and an external scale which grows at a constant rate. Consider an alloy A-B in which A is the more-noble solvent and B is the less-noble solute. The solute oxide is assumed to precipitate in internal oxidation as stoichiometric particles of BO, which are assumed to have a negligible solubility product in the A matrix. It is further assumed that the diffusion of oxygen through the zone of internal oxidation is not impeded by previously precipitated BO, particles. The specimen is assumed to be sufficiently thick so that the solute content in the center of the specimen does not change. For consistency: the following notation of Maak2 has been used:

Jan 1, 1967