Search Documents

By R. E. Goodman, M. A. Mahtab

The, state of stress around a vertical wellbore in rock following nonlinear stress-strain laws is examined by means of finite element analysis. The wellbore, is considered an axisymmetrier body with axisymmetric loading. The initial vertcal and borizontal stresses are "locked" in the rock elements around the wellbore and a new state of stress is generated by the displacements which occur around the borehole. A point-uise variation of the elastic moduli is made on the basis of the new stress state and the triaxial data. The initial stresses tire now reintroduced along with tile changed moduli and original boundary constraints. This procedure is repeatrd until convergent stresses are reached. The effect of nonlinearity on streases is examined for a 6,000-ft wellbore in a scbistose gnezss and Berea sandstone using results of laborratory triaxial compression tests. The results sbou that the effect is restrited to one well radius from the bottom periphery of the bole. Beyond a distance of one-quarter radius, the elfect of nonlinearity on stresscs IS almost always less than percent for the cases conszdered, The consideration of a static pressure inside the well does not magnify the effect of monlinearity on borehole stresscs. INTRODUCTION The terms " wellbore" and "borehole" here designate cylindrical openings in the ground with vertical axis and a circular cross-section. A knowledge of the stress redistribution that occurs on excavating a wellbore is important in understanding (1) the behavior of the lined or unlined hole, (2) hydraulic fracture response, and (3) the effect of stress redistribution on drillability; also it is important in predicting initial stresses in the virgin ground, and in analyzing the response of measuring instruments placed in the borehole. Our knowledge of the state of stress around a wellbore has been restricted to homogeneous, isotropic, elastic material and derives chiefly from the analysis by Miles and Toppingl and the photoelastic work of Galle and Wilhoit2 and Word and Wilhoit.3 In this investigation the state of stress is examined for a nonlinear elastic material by means of finite element analysis. Many rocks possess stress-strain curves that depart notably from straight lines in their initial or final portions. While the literature contains abundant stress-strain data from triaxial tests (axisymmetric loading) on cylindrical rock specimens, there is little information on rock deformability under nonaxisymmetric loading conditions such as occur at each point around the bottom of a wellbore. Although there is some knowledge of the effect of intermediate principal stress on rock strength, there is virtually nothing known about its effect on rock deformability; therefore, we have assumed here that the effect of intermediate principal stress can be ignored. A schistose gneiss4 and Berea sandstone5 were selected as representative rocks for this analysis. The traditional graphs of deviator stress (ol -o3) vs axial strain were reworked to give the tangent modulus as a function of the deviator stress for varying values of the minor principal stress. The result is a nesting family of skewed, bell-shaped curves for tile gneiss (Fig. 1A) and the sandstone (Fig. 2A). A similar replotting of the lateral strain data defines the variation of Poisson's ratio (v) with the deviator stress and confining pressure. These curves, shown in Fig. 1B for the gneiss and in Fig. 2B for the sandstone, are not so well ordered as the tangent modulus curves. However, all of these display an increase of v with deviator stress application, but the rate of increase diminishes with confinement. The ET and v curves for the two rock types are tabulated in Tables 1A and 1A for use in a digital computer so that material properties corresponding to a given state of stress can be assigned by interpolation. STATEMENT OF THE PROBLEM Given a state of stress in the ground prior to drilling a hole (initial stresses), we assume that

Jan 1, 1969

By J. C. Savins

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

By I. F. Jackson, H. L. Hartman

Hard-metal inserts have received considerable attention in the mining industry recently for a variety of excavating and penetrating purposes. Their use in percussive drilling and in coal cutting machines is well known, but attempts to utilize them in cutting materials such as dimension stone have been discouraging because of excessive tool wear and failure. Particular attention is given to those properties of hard-metal inserts that affect their resistance to wear. A consideration of the actual wear mechanism shows the importance of correct mounting of the inserts and of the selection of suitably powered driving motor so as to achieve minimum tool wear and maximum tool life. Comparisons are made between the performances of types of hard-metal inserts, relating their inherent properties with results obtained in laboratory tests. As part of a broad program to improve the quarrying and processing of industrial minerals, The Pennsylvania State University has been investigating the machinability of slate under laboratory and field conditions. It should be emphasized that the results of the research presented in this paper are also applicable to other natural structural materials, and particularly to the more workable ones such as limestone and sandstone. Improvements which accrue in cutting and sawing should prove of benefit to all the dimension stone industry. Further, it may be pointed out that the results reported here are applicable to rotary drilling (with drag bits) as well as to cutting and sawing of rock. INVESTIGATIONS IN RELATED FIELDS The investigation was initiated with an extensive survey of available literature related to metal cutting and to rotary drilling. The following similar features between metal cutting and slate cutting became evident: 1) Hard-metal inserts are widely used as a cutting tool in both cases. 2) The forces acting on the tool may be resolved into normal and tangential components. 3) Rate of penetration and tool wear are important criteria in each case. On the other hand, important differences also exist: 1) Most metals, exhibit plastic behavior, while slate (and most rock) is a brittle material under normal conditions. 2) In metal cutting the component of force acting parallel to the surface of the workpiece is considerably greater than the force holding the tool against the workpiece. 111 slate cutting, this situation is reversed. Similarities between drag-bit rotary drilling and slate cutting are as follows: 1) In both casts, a tool, usually of the hard-metal insert type, is used to penetrate a brittle material. 2) The drilling forces and related variables, such as rotation speed, are similar. 3) Chip formation in each case is essentially the same. 4) The shape of the tool is similar in rotary drilling and slate cutting. Important differences are also apparent: 1) The objective of drilling is to produce a directed hole, while in cutting it is to form a kerf. In core drilling, however, a kerf is also produced. 2) The horizonta1 component of force in rotary drilling decreases from the center to the outside of the tool, whereas in slate cutting the component is almost uniformly distributed across the width of the tool. 3) The power absorbed in the drill rod is greater in rotary drilling. 4) The clogging action of chips, if excessive power is applied to the tool and if flushing is insufficient, is more serious in rotary drilling. Because of the striking similarities between metal cutting and slate cutting, and between rotary drilling and slate cutting, it was considered that the findings of investigators in those fields would serve as a guide for experimental studies of slate cutting. The following is thus devoted to a discussion of those aspects of metal cutting and rotary drilling pertinent to the present investigations.

Jan 1, 1962

By A. Lubinski, J. E. Hansford

Rotating drill pipe passing through dog-legs suffers fatigue damage due to cyclic bending stresses. Curves of the cumulative fatigue damage incurred in either steel or aluminum drill pipe going through Such dog-legs are presented. These curves provide means for properly planning the rate of angle build-up in directionally drilled wells and allow discarding of drill pipe which has been fatigue-damaged in dog-legs. INTRODUCTION Drill pipe which is subjected to bending as it rotates, experiences a cyclically varying bending stress. If this bending stress is sufficiently large, the drill pipe is fatigue-damaged by each revolution, and ages as fatigue damage accumulates. Drill pipe can be bent intentionally, as in the case of directionally drilled holes, or unintentionally due to unde-sired changes in hole direction caused by variations in formation properties, weight on bit, etc. These changes in hole direction, whether intentional or unintentional, are called dog-legs. The theoretical limits of dog-legs which are permissible before the drill pipe begins to fatigue can be calculated.' However, in wells where the maximum permissible dogleg is exceeded, the concept of cumulative fatigue damage can be used to estimate damage suffered by the drill pipe. This concept has already been used to estimate the cumulative fatigue damage of drill pipe due to drilling vessel roll.' The purpose of this paper is to extend the usage of the cumulative fatigue damage concept to drill pipe in doglegs, and to draw conclusions which may be beneficial in preventing numerous frustrating and expensive drill pipe failures. FATIGUE OF STEEL DRILL PIPE IN IDEALIZED GRADUAL AND LONG DOG-LEGS It is taken for granted that, if necessary, the reader will become familiar with the definitions and concepts which follow by examining the reference cited in each case. For Ref. 1, it will suffce to read only through the section Abrupt Dog-Legs. 1. Concept of the two extreme kinds of idealized doglegs: gradual and long dog-legs and abrupt dog-legs,' 2. Concept of why the bending moment in drill pipe un- der tension is the largest in the vicinity of the tool joint,' 3. Concept of why dog-leg severity which drill pipe can withstand decreases with an increase in tension,' and 4. Definition of the rate of change of over-all hole angle (dog-leg severity)." Figs. 1 and 2 indicate the aging (fraction of life expended) of steel drill pipe passing through a gradual and long dog-leg. Both figures have a left-hand ordinate scale labled Length of Drill Pipe Below Dog-Legs. This scale can be used for cases in which the drill collar weight is applied to the bit and does not contribute to the tension in the drill pipe. It can be used for either 3 1/2-, 4 1/2- or 5-in. drill pipe. This is due to the fact that, if the drill collars do not contribute to the drill pipe tension, then the effect of pipe tension on bending due to a change in pipe size is nearly compensated for by a corresponding change in pipe rigidity. On the other hand, if drill pipe is rotated off bottom, then drill collars contribute to the tension in the drill pipe and this increase in tension is not compensated for. In this case, one of the three right-hand ordinates must be used (the one which corresponds to the correct drill pipe size). In the right-hand ordinates, tension should be understood as the weight in mud of the portion of the drilling string below the dog-leg, minus the weight on bit.

Jan 1, 1967

By G. W. Neighbor, M. S. Abrahams

The actice slip plane in InSb is found to be of the {111} type by using the method of two-tmce anulysis. Measurements of the rotation of the tensile axis with increasing plastic shear strain indicate that the tnacroscopic slip direction is a <110>. Recrystal-lization is first observed at a sheav strain of about 70 pct for a crystal deformed at a temperatzcre of 473° C and at a shear strain rate of 2.90 times 10- %ecc-I; the dislocation density corresponding to this value of the strain is 1 x 108 em-&apos;. At a shear strain of about 170 pct, the recrystallization process is complete. This appears to be the fimt time re-crystallization has heen observed in a semiconductor. V HE operation of the {Hi} glide plane in Ge was first noted by Gallagher.&apos; This finding was confirmed by Treuting who also extended the result to Si. Using X-ray techniques, Treuting3 demonstrated that the macroscopic slip direction in Ge was of the <110> type. The choice of the (111) <110> slip system in Ge is exected since in the diamond-type structure the {Ill7 plane is the plane of highest density and the <110> direction is the direction of closest packing.4 There is little reason to expect that the glide elements of InSb are different from those of Ge, since the zinc-blende structure characteristic of InSb is quite similar to the structure of Ge. With regard to InSb, Allen5 has indicated that the glide plane is the {111} type. There appears to be no quantitative determination of the glide direction, however. The present work was undertaken to establish both the glide plane and glide direction in InSb and also to see if recrystallization occurs at large shear strains. EXPERIMENTAL PROCEDURE Mechanical Deformation. The single crystals of InSb used in this investigation were deformed in uni-axial tension. The tensile axis of the specimens was [145] and the periphery was bounded by (111) and (321) planes as shown in Fig. 1. The metallographi-cally polished specimens were deformed in a dynamic atmosphere of He, using an Instron machine. The tensile tests were performed at temperatures in the range of 300°C to the melting point (525OC) and at shear strain-rates varying between 1.94 X 104 sec-&apos; to 1.94 x 103 sec-&apos;. Additional details concerning the preparation and deformation of the specimens have been presented previously.&apos; Analysis of the Glide Elements. The traces of the active slip plane in both (111) and (321) planes were analyzed according to the usual method of two-trace analysis." The macroscopic slip direction was determined by measuring the rotation of the tensile axis as a function of plastic strain,8 These measurements were made from Laue back-reflection X-ray photographs. The X-ray beam was perpendicular to the (111) plane. A given crystal was used for only one measurement of the rotation of the tensile axis after it had been strained to a predetermined value. Therefore, care was taken at all times to insure that the X-ray beam was parallel to the [Ill ] direction of all the samples. RESULTS AND DISCUSSION Glide Plane. After 3.4 pct plastic shear-strain, at a temperature of 369°C and a shear strain-rate of 1.94 x 104 sec-&apos;, the slip traces on the (111) and (3.21) planes appear as in Figs. 2(a) and 2(b), respectively. The results of an analysis of these glide markings demonstrate that the slip plane is (111). The operation of this plane is expected on the basis of the orientation of the specimen, since glide on this plane and in the [I011 direction will result in

Jan 1, 1962

By J. Geertsma

The relationship between porosity and the speed of propagation of acoustic waves in fluid-saturated porous rocks as measured by the Sonic log and by ultrasonic tecbniques is analyzed. Biot&apos;s continuum theory 1 is used to explain the difference in acoustic wave propagation between a dry and a liquid-saturated porous material. The porosity is a variable in this theory. However, the acoustic wave propagation in the dry rock depends too on porosity, and this dependence is not predicted by the theory. Frequently in dry sandstones, a nearly linear relationship between reciprocal acoustic wave velocity and porosity is observed in the low-porosity range. The physics behind this behavior is outlined. An empirical relationship of the form, 1/v = A + F, applies accordingly for many porous dry rocks, provided the porosity is the only variable. The presence of a liquid in the pores changes the value of B, and this change is found to be in agreement with the Biot theory. The time-average relation introduced some years ago2 results in an equation of the same type — 1/V = F/Vf + (1 - F)/Vr — but is not based on a sound physical picture. Still, this relation sometimes predicts approximately correct A and B values. Carbonate rocks with their complicated pore structures sometimes show a different relationship between wave velocity and porosity, unfavorable for log interpretation. Examples are presented. The simultaneous presence of calcite, dolomite and anhydrite, with their different grain densities and matrix compressibilities, complicates acoustic-log interpretation in carbonate rocks still further. Other complicating effects of acoustic-log interpretation are discussed. They are related to the influence of shale streaks and natural fractures on the average wave velocity observed by the logging tool and to the effect of adsorption phenomena on wave propagation in unstressed rocks particularly in sandstones. INTRODUCTION The velocity of propagation of sound waves in porous sediments as a function of depth is a quantity frequently measured. Velocity loggers of various d-sign may be used, for which velocity is defined as the. distance between a wave generator and a wave detector (one-receiver system) or the distance between two wave detectors (two-receiver system), divided by the shortest time required for a vibrational pulse to cross this distance. Velocity loggers ale used to assist in seismic prospecting, to differentiate between the various types of sedimentary rock layers, and also to determine rock porosity and pore fluid content. This paper is especially concerned with the relation between velocity and porosity and the effect of the pore fluid. Most velocity-log interpretation in terms of porosity is based on a formula advocated by Wyllie, et al,2 suggested earlier by Hughes and Jones,3 and known as the "time-average relation". This formula app1ies for a model consisting of a layered system of parallel alternating slices of a solid and a liquid, crossed by the wave path perpendicular to the interfacrbs. This must be an unsuitable model for the derivation of wave propagation properties of liquid-saturated porous media. It suggests that only rock matrix and fluid properties influence the wave velocity. The surprising fact, nevertheless, is that applications of this formula to clean sandstones under specific conditions (sufficient effective stress) results in porosity predictions which are close to reality. Other authors, such as Gassmann,4 White and Sengbush,5 Brandt,6 and Hicks and Berry,7 make use of more realistic models to arrive at a wave velocity formula for porous sediments. Their models consist of various liquid-saturated packings of frictionless spherical grains and are, therefore, also highly idealized. A number of experimental results obtained in the laboratory can be better explained with these model theories than by the time-average relation. A study of a general character, not depending a priori on a special model of the porous structure, was still lacking.

By J. W. Downey, R. A. Morris, M. V. Nevitt

An experimental survey has determined the occurrence of Ti2 Ni-type phases in binary and ternary alloys of Ti, Zr, or Hf with first, second, or third long-period transition metals and also in ternary systems in which 0 is the third component. The existence of eleven previously unreported isomorphs is cited and some approximate phase boundaries are given. Lattice parameters are reported for the new phases. The periodic table relationships involved in the occurrence of the phase are discussed. A previous paper&apos; reported the occurrence of a group of Ti2Ni-type phases (structure-type E93) containing Zr in combination with Rh, Ir, and Pt, and Hf in combination with Rh, Pd, Ir, and Pt, and cited evidence to show that all of the Zr phases and the Hf-Pd phase are not binary phases but occur only when 0 is present. More recently one of the present authors found that in the ternary systems of Ti and 0 with Cr, Mn, Fe, Co, and Ni the occurrence of the Ti2Ni-type phase could be correlated qualitatively with valency electron concentration, and that the role of 0 appeared to be that of an electron acceptor.&apos; This paper describes the results of a broader investigation of the occurrence of Ti2Ni-type phases in binary and ternary systems containing Ti, Zr, or Hf. EXPERIMENTAL PROCEDURE Alloys were prepared by arc melting. The purity of the metals Zr, Rh, Hf, Ir, Pd, and Pt was given in the previous paper.&apos; The additional materials used for this work were metals with a minimum purity of 99.9 pct and reagent grade ZrO2, TiO2, Cr2O3, and NiO. Chemical analysis of several alloys in each system indicated good agreement between intended compositions and cast-button compositions. For the heat treatment precautions were taken to avoid oxygen contamination. The specimens were wrapped in Mo foil and sealed in fused quartz tubes at a room temperature pressure of 1 x mm of Hg. These capsules were then annealed in a vacuum furnace at an external pressure of 5 X mm of Hg. Since the major constituent in most of the alloys studied was the brittle Ti2Ni-type phase, Debye-Scherrer patterns with sharp lines could be obtained with powder crushed from the annealed solid samples without being reannealed. Metallographic, X-ray diffraction, and density determination techniques previously described1, 2 were used. RESULTS AND DISCUSSION Zr-Rh-O, Zr-Ir-Oand Zr-Pt-O Systems-Figs. 1 through 3 show the approximate boundaries of the Ti2Ni-type phase fields at the temperatures indicated in the ternary systems Zr-Rh-O, Zr-Ir-O and Zr-Pt-0. In each isothermal section the Ti2Ni-type phase field is labeled E and the unidentified binary phase found in the previous work&apos; is labeled ?. The dotted lines, which are not phase boundaries, correspond to the composition (Zr2 B)1-xOx. In each system the Ti2Ni-type phase field is displaced to the high-Zr side of this composition line. However, the departures from Zr2B stoichiometry for the metal components are not as high as Zr3B, as was inferred from the previous work. Lattice parameters of the Ti2Ni-type phases in the alloys shown as triangular data points in Figs. 1 through 3 are given in Table IA. The probable compositions of the Ti2Ni-type phases have been estimated from the isothermal sections. Although there are a few exceptions that probably result from uncertainties in phase boundary location and from the effect of variation of oxygen content, the generalization can be

Jan 1, 1961

By W. C. Winegard, J. E. Gruzleski

The stages in the development of cells in the Sn-Cd eutectic have been studied by unidirectionally solidifying specimens under known conditions of growth rate, temperature gradient, and impurity concentration. By means of a quenching technique, the shape of the solid-liquid interface during solidification could be observed. Cells are shown to develop from defects or depressions which exist on the solid-liquid interface of even a pure eutectic. Th,e container wall, grain boundaries , and fault lines in the micro structure play an important role in cell development. A two-phase eutectic dendrite found at high impurity concentrations is described. L HE cellular or colony structure formed during eutectic solidification caused much confusion among early workers in their attempts to classify eutectic structure, and a major advance in our understanding of eutectic solidification was made when the colony structure was identified with a cellular interface caused by the presence of impurity elements.1&apos;2 To date, however, no detailed investigation of the stages in cell development in eutectics has appeared, and little work has been done on the cell morphology. The present paper presents the results of such an investigation. EXPERIMENTAL The Sn-Cd eutectic was used in this investigation with lead as the impurity element to promote cell formation. Materials of 99.999 pct purity were zone-refined to an approximate purity of 99.9999 pct. Alloys as close as possible to the eutectic composition, 67.75 wt pct Sn, were prepared using zone-refined tin and cadmium. The alloys were then zone-melted to the exact eutectic composition, as described by Yue and Clark,3 using a total of twenty to twenty-five passes on each eutectic bar. Only those portions of the bars which contained no primary phases were used in subsequent experiments. Master ternary alloys were prepared by melting together the required amounts of lead and the Sn-Cd eutectic in a sealed Pyrex tube under argon, and these alloys were then diluted with pure Sn-Cd eutectic as required. Unidirectional solidification was conducted vertically upward using the two-element resistance furnace shown schematically in Fig. 1. With this apparatus, it was possible to obtain a range of growth rates and temperature gradients so that the cell structure could be studied under different growth conditions. Growth rates in the range of 0.4 to 9.1 cm per hr and temperature gradients from 2.25° to 1475°C per cm were used. To ensure steady-state solidification, power to the furnace was stabilized, and the entire apparatus was enclosed in a large box in which the temperature was controlled to within ± 0.2°C. The apparatus was arranged so that the specimen could be quenched in situ during solidification by pouring a large quantity of water down the central vycor tube. Excellent delineation of the solid-liquid interface at the time of the quench was obtained because of the structural changes accompanying the rapid solidification of the remaining liquid. The quenching rate was always extremely rapid as no evidence of a transition in the lamellar spacing at the interface was ever observed. The alloys to be solidified were contained in degassed, high-purity graphite tubes, 1 cm OD and 0.6 cm ID. The specimens were about 18 cm in length. After the alloy had been cast into the tubes, three 0.8-mm holes were drilled in the tube to the center of the specimen. One hole was placed 1 cm from the bottom of the ingot, and the other holes were placed at 8 and 9 cm, respectively, from the base of the ingot. Thirty-four-gage chromel-alumel thermocouples were sealed in the holes. The lower thermocouple was used to position the solid-liquid interface at the start of each experiment,

Jan 1, 1969

1-Iron and Steel JOHN A MATHEWS, Chairman Ralph H SWEETSER Vice-chairman WILLIAM CAMPBELL, Vice-chairman MAX ROESLER, Secretary Iron Ore CHARLES B MURRAY, Chairman R C ALLEN W O HOTCHKISS MAX ROESLER ERNEST F BURCHARD THOMAS T READ DWIGHT E WOODBRIDGE M M DUNCAN WILLIAM KELLY CARL ZAPFFE C T FAIRBAIRN CHARLES F RAND

Jan 1, 1923

By John P. Hager, Igor A. Wilkomirsky

The thermodynamic properties of the Pb-Ag system have been determined between 775° and 950°C by the cell: Electrotransport measurements on silica-saturated PbO-SiO2 melts established that the conduction was entirely ionic and that lead is present only in the divalent state. The activities of lead and silver exhibit Positive deviations from ideal behavior. In dilute solutions the activities of lead and silver follow Henry&apos;s law. The values obtained for yhb and yig at 1200°K are 1.49 and 1.69, respectively. Determination of a portion of the liquidus curve agrees well with the established phase diagram. provides an accurate means for the thermodynamic study of liquid alloy systems. However, Wagner and wernerl have shown that cells of type I are limited to alloy systems for which the reaction at the alloy-electrolyte interface will be extremely small. The traditional electrolyte for cells of type I has been a low melting point mixture of group IA chlorides, to which a small amount of the compound ACIZa has been added. However, the choice of a chloride electrolyte may not always be practical, in that the Gibbs free energies of formation of the compounds ACIZ, and BCIZb may not be sufficiently different to minimize Reaction [I]. On the other hand, the Gibbs free energies of formation of the oxides of metals A and B may be sufficiently different to minimize Reaction [I]. In this case, the use of molten oxide electrolyte would be clearly advantageous. In the present study, the thermodynamic properties of the Pb-Ag system have been investigated using the cell: Previous galvanic cell studies on the Pb-Ag system2&apos;= are not reliable due to the displacement reaction between the alloys and the chloride electrolyte. Wagner and ~erner,&apos; using Laitinen and Liu&apos;s&apos; value at 450°C for the standard electrode potential of Reaction 111, have estimated the percent error in the activity of lead bPb) when using a chloride electrolyte. At 1000°C the error for xpb = 0.10 is 12 pct. Taking the Gibbs free energies of formation for PbClz and Ag cI,&apos; and the equations of Wagner and Werner,&apos; the error in apb at 1000°C and Xpb =0.10 was calculated to be 140 pct. However, using the Gibbs free energies of formation for Pb O&apos; and an estimated value for Ag 20, the error in apb when using a molten oxide electrolyte at 1000°C was estimated to be less than 1 pct for The virtual reaction of cell I1 for the passage of n Faradays of positive electricity from left to right in the cell is: where 5 is the Faraday equivalent (23,063 cal per v equivalent), E is the electromotive force of the cell in volts, and n is the number of Faradays of electricity required for the transfer of 1 g-atom of lead through the cell. The value of n should be equal to 2 since lead is generally divalent in oxide melts. However, the presence of either tetravalent lead ions or electronic conductivity through the molten oxide electrolyte will result in a value of n greater than 2. Bockris, Kitchener, and Davies&apos; have determined the amount of ionic conductance in a liquid PbO-SiO2 mixture (xPbO= 0.50) at 920°C by measurement of the volumes of CO and C& evolved at a carbon anode during electrolysis. From their data a value of 2.13 * 0.16 was calculated for n. Although their measurements indicate the possible presence of multivalent lead ions and/or electronic conductivity, their measurements were not at silica saturation and no definite conclusions can be drawn. However, as a result of the electrotransport measurements described in the next section the ;due of n can be taken as equal to 2.

Jan 1, 1969

By V. M. Dobrynin

Experimental data demonstrate that physical properties of porous rocks change under pressure. In this paper an assumption is made and proved that under pressure the changes of physical properties such as porosity, density, permeability, resistivity and velocity of elastic waves are controlled to a large extent by the pore compressibility of rocks. It is also shown that the pore compressibility of rocks can be determined, within the range of pressures from 0 to 20,000 psi, by knowing the maximum pore compressibility and the magnitude of the pressure. Mathematical equations were developed which permit one to define changes in physical properties of porous rocks under pressure. These equations were verified by experimental data obtained from the study of sandstones. INTRODUCTION In studying the behavior of porous rocks under pressure in the field of petroleum technology, the most interesting aspect is the observation of those properties which characterize the rocks as possible reservoirs — for example, porosity, permeability, resistivity, density and the velocity of elastic waves. The literature dealing with this problem mainly contains data concerning the study of only one or at most two of these parameters, but not of the group as a whole. An attempt is made in this paper to find general equations involving each of these parameters, which will permit the study of the behavior of rocks under pressure. All experimental data used here were obtained from the investigation of consolidated sandstones. EXPERIMENTAL In addition to the use of published experimental data, an experiment was carried out which studied the main physical properties of sandstones under pressure. Two homogeneous quartz sandstones were chosen for this purpose: (1) the Torpedo sandstone from Kansas, and (2) the Medina sandstone from Ohio. The porosity of the Torpedo sandstone was 20.2 per cent, and that of the Medina 8.7 per cent. Permeabilities were 45 md and less than 1 md, respectively. Each sandstone contained about 5 per cent clay minerals, consisting mostly of kaolinite and chlorite, which were distributed quite evenly throughout the samples. One cylindrical sample 2 in. in diameter and 5 in. in length was cut from each sandstone and then saturated in a vacuum with a 3N solution of NaC1. This high concentration was used in order to obtain true formation factors and to decrease the swelling of the clay minerals. The methods of mounting the samples and measuring the changes in porosity and resistivity were practically the same as those described by Fatt and Mann.1,2 Changes of resistivity under pressure were studied for sandstones with 100 per cent water saturation, and for sandstones with the irreducible water saturation. The irreducible saturation was obtained by enclosing the saturated rock samples in relatively fine silicate powder so as to remove the water by capillary action. This procedure is described by Orkin and Kuchinski.3 Changes of permeability with pressure were determined at room temperature using nitrogen as the flowing medium. In studying the effects of pressure, one series of measurements was made using an internal pore pressure Pi equal to the atmospheric pressure, while the overburden pressure Pe ranged from 0 to 20,000 psi. A second series of measurements was used over the same range of overburden pressure, but with an internal pore pressure of 1,800 psi. When the results were compared on the basis of net overburden pressure (Pe — 0.85 Pi ), there was practically no difference for these two sandstones. The origin of the factor 0.85 in the expression for net overburden pressure is given by Brandt,4 Fatt5 and Geertsma.6 The relatively large bulk volumes of the sample (about 15 cu in.) made it possible to observe a

By H. E. Cline, T. H. Alden

A study zvas made of rate-sensitive plastic deformation and superplasticity in alloys of lead and tin. These alloys (lead, PO-19 pct Sn, Sn-38 pct Pb, Sn-2 pct Pb, tin) were extruded, aged at room temperature, or melted and resolidified to produce a variety of ~nicrostrnctures. Flow stress us strain rate data were obtained frorn tension tests for strain rates between 10-6 and 1 min-1. The curve log o .us log £ was S-shaped jor fine-grained materials and approxinzately linear for coarse-grained nzaterials. At constant grain size, the behavior was relatively insensitive to conposition. High rate sensitivity (d log o/d log k) and superplasticity occurred at intermediate strain rates in fine-grained materials and are believed to be associated with a transition with strain rate in the mode of plastic flow. This transition was shifted to higher strain rates with increasing temperataire or decrensing grain size. The high rate rneclzanisrn is presunzed to be slip, while rnetallographic evidence sug-gests that the low rate lnechanism is pain boundary sliding. PEARSON&apos; discovered that extruded eutectic alloys of Pb-Sn and Bi-Sn could be elongated as much as 2000 pct in tension. The ability for stretching decreased with aging time or annealing to increase the grain size. At high strain rates, or if the material was in the cast condition, normal ductility was observed. Following deformation of the extruded Pb-Sn eutectic at low rates, Pearson observed an equiaxed granular structure on an initially flat surface, while rapid straining produced elongated grains. He attributed the large deformation in the extruded alloys to the persistence of extremely small grains to a temperature sufficiently high for inter granular flow. Avery and backofen2 measured the strain rate sensitivity of the flow stress, m = d In old In 6, in the Pb-Sn eutectic and established a correlation between m and the tensile elongation. In the extruded eutectic alloy the strain rate sensitivity was 0.1 to 0.2 at rapid rates and the elongation was small, while at low rates m increased to 0.4 to 0.6 and superplastic elongation was observed. They also noted a grain size effect. In the case of the cast Pb-Sn eutectic the strain rate sensitivity remained at 0.1 over the range of strain rates of their experiment (3 x 10"3 to 2 min-I). These authors suggested that diffusional (Nabarro-Herring) creep was the mechanism by which deformation occurred at low strain rates, while climb-controlled creep occurred at high rates. Two predictions of this model are that, at low strain rates, the strain rate sensitivity approaches 1, i.e., Newtonian viscosity, and that the strain rate at constant stress varies as D/L2, where D is the coefficient of self-diffusion and L the grain size. More recently, a mechanical and microstructural investigation of the dilute superplastic alloy, Sn-5 pct Bi, has been done.3 The results show that the high strain rate sensitivity is associated with a transition with strain rate of the mechanism of plastic flow. At low strain rates, grain boundary sliding and its associated accommodation processes were evident, while at high strain rates quasi-uniform slip was concluded to be the dominant deformation process. The following experiments were performed to extend our knowledge of the phenomena of superplasticity to a larger range of strain rates, compositions, temperatures, and microstructures in the Pb-Sn alloy system. The question of why a cast structure does not exhibit viscous behavior although it appears to have a fine microstructure2 was investigated.

Jan 1, 1968

By Wallace W. Wilson

This DaoerAA calls to attention for the first time, to the writer&apos;s knowledge, a purported recovery of oil by substantially horizontal gas cycling that is considerably in excess of 50 per cent of the oil originally in place. The average angle of dip along the path of travel of the gas is less than one degree of arc. The best recovery by gas cycling previously brought to the writer&apos;s attention1 occurred in the case of the Mile Six Field in Peru, where about 68 per cent of the original oil in place has so far been recovered. In this case, however. the recovery was greatly aided 1)y the high dip angle of about 17.5". An approximate estimate of recovery was made by the method developed by the writer&apos; for the case of the Cook Ranch Field, using an average relative permeability ratio function dictated by past experience with a number of good quality sandstones. A result of about 40 per cent recovery of the original oil in place wa* obtained after the 8.5 pore volumes of cumulative gas cycling undergone by the Cook Ranch Field. Tile writer desires to emphasize what appears to him as the sound engineering practiced at Cook Ranch, and also the care attending the taking of data and the high quality of the calculations made. It is impossible to take issue on factual grounds with any of the conclusions drawn in the paper. It is, however, desired to point out the difficulty of explaining the unexpectedly favorable oil recovery on ordinary theoretical grounds. It would be of interest to know whether relative permeability ratios have been measured on the Cook sand. The operator is inclined to attribute some of the excess benefit to settling of the oil across the bedding planes into the lower part of the sand. It is difficult for the writer to appreciate the importance of such an effect, inasmuch as the concomitant opening of excessive vacated pore space in the upper part of the sand would he expected to cause excessive bypassing of the gas. REFERENCE 1. Welge. H. J.: "A Simplified Method for Computing Oil Recovery by (Gas or Water Drive." Trans. AIME, (1952), 195, 91. AUTHOR&apos;S REPLY TO MR. WELGE Welge has poined out clearly the difficulty of applying material-balance procedures in calculating the Performance of low-pressure, dispersed gas injection operations. There can be no argument that even with the most favorable assumptions of K,/K,, and with cumulative gas cycling in excess of normal expectations, the calculated oil recovery will be considerably less than the 72.5 per cent recovery realized from this project to Jan. 1, 1951. There probably are at least two reasons for the unusually high recovery from this reservoir, the effects of which would be difficult or impossible to calculate by known methods: 1. The operator did not use the same injection wells for the entire operation, and exercised considerable control over producing GOR by changing the locations of injection wells whenever well performance indicated excessive bypassing. Fifty-six wells were used for injection, not over 29 of which were in use simultaneously. This procedure resulted in frequent shifting of the fluid movement patterns in the reservoir and undoubtedly increased the ultimate recovery. The writer knows of no satisfactory procedure for calculating the effect of this practice. 2. Vertical movement of oil across the bedding planes of the reservoir rock hy gravity drainage tended to maintain a relatively high oil saturation, and hence a favorable K,/K, ratio, in the lower sections of the reservoir. As pointed out by Welge, this condition has aggravated the tendency for high producing GOR, but at the same time it has resulted in increased oil recovery. As may be seen from the performance charts for the various battery areas, operating (GOR&apos;S have been quite high for the last several years of reported history.

Jan 1, 1952

By Matti J. Saarivirta

The relative machinabilities of the precipitation-hardenable Cu-0.9 pct Cr, Cu-0.25 pct Zr, and Cu-0.9 pct Cr-0.25 pct Zr alloys are improved from 18 to 80 pct by 0.75 to 1.25 pct Pb. The Lead also improves the tensile strength, hardness, ad resistance to softening of the Cu-Cr and Cu-Cr-Zr alloys. When the alloy containing 0.7 to 2.2 pct Cr, 0.10 to 0.25 pct Zr, 0.75 to 1.25 pct Pb, and balance copper is solution-annealed at 1000°C (1830°F) and then aged at 525°C (Y75°F), tensile strength of 65,000 psi, 16 pct elongation, 62Rb hardness, and 84 to 87pct IACS electrical conductivity are obtained. The copper base Cu-Cr, Cu-Zr, and Cu-Cr-Zr alloys have high electrical and thermal conductivities, high strength and resistance to softening at elevated temperatures, but their machinability is poor, i.e., close to that of the pure copper. In the manufacture of various components such as welding electrodes or electrical and thermal conductors it is ofter desirable that the alloy have good machinability for high production rates and long tool life. Thus, alloys that combine these properties may be used not only for existing applications, but for new ones where good machinability is a requirement. This investigation was undertaken to improve machinability without impairing other desirable characteristics of the Cu-Cr, Cu-Zr, and Cu-Cr-Zr alloys. These commercial alloys contain about 0.4 to 1.2 pct Cr, 0.10 to 0.25 pct Zr, and balance copper. Lead was selected as an element that should improve the machinability and possibly other characteristics of the alloy. WORK PROCEDURE This investigation began with a study of the effect of lead on the hardness, tensile strength, electrical conductivity, and machinability of the binary alloy containing 0.10 to 1.30 pct Cr and balance copper. This was followed by an investigation of the effect of lead on the Cu-0.1 to 0.4 pct Zr and Cu-Cr-Zr alloys. The alloys were made from initially oxygen-free copper and for comparison some castings were made from phosphorus-deoxidized Lake Copper. Commercially pure chromium, zirconium, and lead were used for alloying. In melting the chromium was added to molten copper at 1400°C (2550°F), the heat held for 15 min at temperature, and the temperature then lowered to 1350°C (2460°F) when the other elements were added. The alloys were held for 15 min at this temperature and cast into a cast-iron mold. One-kilogram, 3-in.-diam castings were made. An argon cover was maintained during melting and casting to prevent oxidation. The chemical compositions of the alloys are shown in Table I. The micro structure and properties of the alloys were studied in the cast and wrought conditions prior to and after various solution-annealing and heat-treating cycles. The relative machinability was determined as a function of the lead content and compared with the free-machining Cu-lpct Pb alloy. To determine the machinability, 13/16-in.-diam cast specimens were cut by a machine hack saw using a constant cutting speed of 50 fpm, a fixed feed force, and a new blade for each alloy specimen. Molybdenum highspeed steel saw blades, 0.032 in. thick with 14 teeth per in., were used in all tests. No cutting fluid was used. Five tests were made on each specimen and the average elapsed time per test was taken as a measure of the relative machinability for that alloy. The variation of cutting time for each alloy was less than 110 pct.

Jan 1, 1969

By D. Harker, W. E. Seymour

CERTAIN alloys of iron and nickel, when rolled and annealed, possess a preferred crystal orientation: (001) in the rolling plane and [loo] in the rolling direction, when recrystallized at 850" to 1050°C after approximately 98 pct cold reduction. Since the preferred orientation makes a direction of easiest magnetization coincide with the rolling direction, these alloys, especially a 50 pct iron in nickel alloy, have found wide application in the electrical industry&apos; for choke coils and special types of transformers. The rate of the recrystallization reaction was studied at temperatures ranging from 500° to 600°C. The heat of activation for the reaction was calculated from the observed rates, and crystal orientation determinations were made before and after re-crystallization. A vacuum melted iron-nickel alloy* analyzing 49.6 pct nickel, 0.018 pct carbon, and 0.30 pct manganese was used for the experiments. The alloy was cold-rolled 98 pct into strips 1/2 in. wide by 0.002 in. thick. Specimens 1 in. long were sealed under a vacuum (approximately 10-8 mm Hg) in 1/2 in. ID pyrex glass tubes. For heat treatment these tubes were fastened to Nicrohm wires and submerged in a molten salt bath controlled to ±2°C. (The maximum measurement error was within 3°C.) The time at temperature was varied logarithmically from one sample to another, and runs were made at 500°, 525", 550" and 575°C. A molten lead bath was used for a run at 600°C which was controlled to the same accuracy. To determine the extent of recrystallization after a particular time at a given temperature, a method was used suggested by the work of Decker, Asp, and Harker.3,4 This method employs an X-ray spectrometer whereby X rays are diffracted from crystallites whose diffracting atom planes are parallel to the rolling plane of the specimen. The intensities of the diffracted rays are measured by a Geiger counter. The counter can be rotated through a range of angle, and can be motor driven through this range at a speed of 2" per min. The measured intensities are plotted vs. angle by a potentiometer recorder. In a 50 pct iron-nickel alloy it was found that only (200) reflections were obtained from the rolling plane with recrystallized material, and that only (220) reflections were obtained from the rolling plane with unrecrystallized material. No other reflections were obtained under these conditions over a range of Bragg angle of 0 to 45" using either type of specimen. As a consequence, the intensity of the (200) reflections from the rolling plane was taken to indicate the amount of recrystallized ma-terial present in a specimen, and the corresponding intensity of the (220) reflections was taken to indicate the amount of unrecrystallized material present. To insure flatness and proper alignment of the reflecting specimen in the spectrometer, the sample was laid in a slot 1/2 in. wide x 0.001 in. deep which was machined in a 3/4x2x1/2 in. steel carrier. The specimen then was taped down with cellophane tape. This steel cradle was mounted vertically in the specimen holder of the spectrometer so that the bottom edge of the carrier rested against a horizontal shoulder. Intensity data used in plotting the percent of recrystallized material vs. time at temperature were obtained by directly counting the diffracted beam of nickel filtered copper Ka radiation used throughout the work. The X-ray spectrometer was employed also in obtaining pole figure data to be used for orientation determinations. At the beginning of this investigation a transmission pole figure holder was used." As a thickness of 0.002 in. of this alloy was found to be opaque to the copper Ka radiation, it was necessary to reduce the thickness of the pole figure specimen. A cold-worked sample was etched in a 50 pct HCl solution for 5 min and was found to transmit a

Jan 1, 1951

By Paul B. Crawford, Melvin B. Burton

A procedure is given which describes how the gelatin model may be-used to study the effect of mobility ratio on areas swept before and after breakthrough in fluid injection programs. Differences in mobility ratio are achieved through the proper selection of electrolytes for the invading and displaced fluids. The progress of the front is followed visually and photographed or traced onto plates (it desirrd positions before and after breakthrough. Example applications to standard flooding pattrrrls, vertically and horizontally fractured reservoirs and the in-situ combustion recovery process ore presentrtl. INTRODUCTION In water flooding, gas cycling and other fluid injection programs, the sweep efficiencies and swept areas after breakthrough are known to depend on the mobility ratio. Previous investigators have described the application of the X-ray technique, stepwise use of the poten-tiometric model and numerical procedures for studying these problems.1,3,4,6 Others have described the application of the electrolytic and gelatin models for studying sweep efficiencies at a mobility ratio near one.2,5,7,8,9 It is a purpose of this paper to describe a procedure for using the gelatin model to study swept areas at mobility ratios other than one. PROCEDURE The technique developed for gelatin model studies is composed of four phases: the construction of the model, the selection and preparation of electrolytes, the running procedure, and the recording of the swept areas. Each phase is discussed independently in order that the over-all procedure will not be confused. Construction of the Model Since some of the advantages of the gelatin model are its economy, and speed, it is desirable to have as small a model as is possihle without creating error through the inability to follow visually the progress of the colored front. For this work, it was found that a distance of 4 to 6 cm between salt bridges is convenient to use. Since the size is relatively small, it follows that the gelatin layer must be thin. For this study, 1/8 in. was used as the gelatin thickness in each case. The model was constructed on a flat 1/2-in. plastic (Lucite) base with boundaries constructed from 1/8-in. plastic strips. The well pattern to be studied was outlined with the 1/8-in. plastic strips on the base with narrow apertures (- - 2 mm) in the boundary to represent the wells. Behind each well opening, a rectangular area of some 8 to 12 sq cm was outlined by additional plastic strips. The plastic strips were fastened to the base by any suitable non-conducting cement so that the boundary system was liquid tight. The electrodes were constructed from sheet copper in a fashion that would allow easy removal. In this work, it was found that by simply forming a large U-shaped clip from a strip of the copper that would slip over the boundary was satisfactory and easily removed. The electrodes should be as wide as one of the sides of the rectangular area behind the well opening. If a fractured pattern is studied, fractures inside the well pattern may be represented by strips of copper of the desired shape mounted in some easily removable fashion. The reason that apertures (or salt bridges) were used for wells instead of metallic electrodes is that the plating of the electroactive ion at the cathode and the introduction of an electroactive ion at the anode has a detrimental effect on the swept area outline if these operations \rere permitted to be carried on within the well pattern itself. The anode may be used within the well pattern

Jan 1, 1957

By J. Fageant, W. C. Ellis

All major regions in a progressively solidified germanium ingot were related through successive orders of octahedral twinning. The occurrence of lineage structure and the generation and survival of orientations are discussed. WHEN metals solidify in a temperature gradient a preferred orientation, or casting texture, usually develops. For example, in face-centered and body-centered cubic metals, most of the solidifying crystals have a cube axis, <loo>, aligned approximately in the freezing direction.&apos; In the solidification of diamond cubic metals this is not the case, since there is a natural tendency to form twin orientations.&apos; . In this paper, twinning relationships in cast germanium are described; some conclusions with respect to nucleation and survival of new orientations are discussed. In the first experiment, an ingot of high purity germanium, % in. in diameter and 1 % in. in length, was formed by progressive solidification4 from the bottom. The ingot was cut into 21 slices, each circular and about 0.030 in. thick. The top surface of each slice was polished and etched to reveal the regions of different orientations. Orientations were determined by the Laue X-ray back-reflection method. The uncertainty in the determined orientation relationships was l" or 2" when separate slices and films were involved. Frequently, however, advantage was taken of illuminating across a boundary when the patterns of the two orientations were recorded on the same film. The uncertainty then was less than 1". Often the twin relation was evident from common zones or individual reflections. Twin relationships were recognized through comparison of the values of the nine angles between cube poles of the two orientations with those calculated for successive orders of octahedral twinning. The birth and survival of each region can be reconstructed by visual examination of the etched slices shown in Fig. 1. In slice 21, where solidification began, many regions are present; these regions were found to be twin-related. Solidification appears to have started from one nucleus formed probably at the wall of the crucible.* Subsequent ori- entations were developed by nucleation on existing solid surfaces, and only in the twin relationship. This is an example of a general process of nucleation on a surface of existing solid and may be descriptively termed oriented nuc1eation.t In germanium, and other diamond cubic substances, the orienting habit is that of twinning on an octahedral plane.&apos; " The prevalence of straight traces in the polished surfaces of succeeding slices (Fig. 1) suggests that the twinning observed in the first slice persisted throughout solidification. This was confirmed by X-ray determinations. The schematic diagram of Fig. 2 shows the twinning relationships of the major regions throughout the ingot. Since all observed relationships are describable as twinning, it is reasonable to conclude that no random nucleation occurred in the freezing of the ingot. An alternative possibility would be that nucleation of a random nature occurred, but only those nuclei in twin orientations survived in subsequent growth. This is unlikely for there would be as good a reason for the growth of randomly oriented nuclei as for the growth of the many different orientations present after several orders of twinning. The many twin components contain a wide distribution of orientations, and did grow. In fact, they are the only ones found in the solid. The conditions of solidification of this ingot were such as to discourage nucleation of crystals wholly in the liquid. The rate of solidification was 0.125 in. per min. This slow rate of heat removal afforded but small opportunity for supercooling within the liquid —a condition needed to provide enough decrease in free energy for such nucleation. When nucleation occurs on an existing solid surface, the interfacial energy between the existing and forming regions is important in determining the new stable orientation. The new orientation is expected to be one which minimizes this energy. In the diamond cubic structure, the interface for the

Jan 1, 1955

By J. Lohrenz, R. J. Francis, G. C. Clark

Compositional material balance methods have been extended to solution gas-gas cap-water drive reservoirs. Partial or complete pressure maintenance as well as reservoir repressuring by gas and/or water injection may be considered. Increrments of reservoir withdrawals and injections equivalent to a given time increment are used as the basic iterative step to allow the extension. An example field study shows how the reservoir volumes. extent of the water drive, and relative permeability ratio function may be estimated. Knowing these, the methods can be used to predict reservoir behavior. INTRODUCTION The compositional material balance is similar to conventional material balances in many aspects. The assumption of a pressure-homogenous reservoir with no transients is common to all conventional material balances1-1 although numerical methods are different. The compositional material balance also contains this assumption. Inclusion of phase equilibria directly in reservoir calculations is the key difference and improvement in the compositional material balance. A continuous account of the compositions of the reservoir fluids and effluents is determined. Calculations of fluid phases and properties can be based on these compositions. Excepting fluid phases and properties determined in this manner, the compositional material balance uses procedures exactly analogous to conventional methods. Brinkman and Weinaug5 and Jacoby and Berry6 have given compositional material balances for simple solution gas drive reservoirs; subsequently Jacoby and Berry&apos; extended their procedure to include pressure maintenance by gas injection at constant reservoir pressure. This paper extends compositional material balance techniques to a more general combination drive reservoir with various types of pressure maintenance operations. The procedure for this extended compositional material balance and an example of its field application are presented. PROCEDURE The numerical calculations of this compositional mate- rial balance are performed on a high-speed digital computer. The most comprehensive case the present compositional material balance program can handle is a solution gas-gas cap-water drive reservoir with (1) water production, (2) water influx, (3) separate or commingled production from the oil and gas zones, (4) any surface processing arrangement including LPG and gasoline production, and (5) gas injection into the gas and/or oil zones." This description of the procedure applies to this most comprehensive case. Appropriate deletions in the procedure can be made for more simple mechanisms. ASSUMPTIONS In common with conventional material balances, the compositional material balance makes the following assumptions: (1) the total reservoir is at a uniform pressure at a given time; (2) the volumetric ratio of gas to liquid withdrawal at reservoir conditions is directly proportional to their mobilities; and (3) instead of assuming a pressure-dependent total gas-oil equilibrium, the cornpbsitional material balance makes a sianificantly different and more rigorous assumption as follows. Thermodynamic equilibrium exists for each component in the liquid-gas mixtures in the reservoir and in surface processing. These three assumptions are fundamental in the compositional material balance procedure. Other assumptions involve arbitrary simplifications of the reservoir mechanisms. This compositional material balance uses the following simplifying assumptions: (1) no bulk mixing occurs between the gas and oil zone; (2) total mixing occurs between any injected gas and the contents of the injection zone; (3) no in situ hydrocarbon liquid is produced from the gas zone; (4) no reservoir gas is retained behind any net bottom-water influx; (5) the residual saturation of hydrocarbon liquid at reservoir conditions behind any net bottom-water influx is constant; and (6) fluids retained behind any net bottom-water influx exert no further influence on the reservoir. These six assumptions are not fundamental. Modifications of the present procedure can be made for different assumptions of the reservoir mechanisms. THE TIME INCREMENT Previous compositional material balances have used a pressure increment for simple solution gas drive reservoirs"" and an increment of injected gas for pressure maintenance.&apos; These increments are suitable for these simple cases. They are, however, inadequate for reservoirs with time dependent drive mechanisms including water drive. A time increment removes this inadequacy and is

By S. C. Miller, R. D. Carter, H. G. Riley

A merhod based on theoretical considerations is presented for estimating srabilized performance of gas wells from data obtained from short-term flow tests. Calculated results are presented for II gas wells evaluated with four-to five-hour flow tests. For seven wells on which comparison war available, satisfactory agreement between measured and calculated performance was obfained. Although the final results are favorable, diflerences in flow capacity k,h obtained from flow tests conducted at different rates were observed. This inconsistency did not seriously affect the calculated results in the test wells. INTRODUCTION In order to negotiate the sale of gas from newly drilled gas wells and to plan otherwise for the orderly exploitation of gas reserves, it is necessary to predict the long-term withdrawal rates physically possible for each well. Since low permeability gas wells do not stabilize except after a lengthy period of flow, the operator is unable, without wasting large amounts of gas, to measure directly the stabilized performance prior to contracting for sale of the gas and connecting the well to a pipeline. This problem is not new, but is one which has recently received considerable attention because of the rising importance of natural gas in the economy and new emphasis on the development of gas reserves. Thus, a method is needed for estimating stabilized performance from information obtained from short-flow tests in which stabilized flow is not attained. Such short tests would be suitable for wells not connected to a pipeline since the tests require only small amounts of time, labor and production of gas. Other requirements of such a method are short calculation time and use of a model of the gas well and its reservoir drainage area which is as consistent as possible with present theoretical, experimental and field knowledge. The method presented herein was developed in an effort to meet these requirements. It has been named the "two-Row method". In addition to the section on the evelopment and description of the method, a description of field application of the method and a discussion of well test results are in- cluded. Finally, there is a section on an additional consideration prompted by some of the test results. The Appendix contains certain details of the derivation of the two-flow method and a sample calculation for Test Well 1. MATHEMATICAL DEVELOPMENT AND DESCRIPTION OF THE METHOD The method is based upon the following assumptions: (1) The system can be described as a finite cylindrical gas reservoir produced by a single cylindrical well located at the center of the reservoir. Flow is strictly radial. (2) Flow obeys Darcy&apos;s law except very near the well where a quadratic flow law is assumed. (3) The system is also characterized by a laminar flow skin factor. The use of a quadratic or turbulent flow term is more fully justified in Refs. 1 and 2. When such a system is opened at the wellbore to constant rate flow from an initial condition of uniform pressure distribution, the radius or drainage quickly passes through the non-Darcy flow and "skin effect" region imrnediately surrounding the well. After the drainage radius or "transient" has passed beyond this region and before it reaches the external boundary, drawdown can be expressed by&apos;." After the drainage radius reaches the external boundary, stabilized flow exists and drawdown is given by&apos;," Eqs. 1 and 2 form the basis for the method. The validity of these equations will now be examined in the light of some well test data previously available. It has been found by actual measurements on gas wells that the relationship between drawdown and flow rate, that is the performance curve, is closely approximated by&apos; Eq. 3 is a straight line on log-log graph paper. Eq. 2 and, for any lixed value of time, Eq. 1 plot as a curve which is slightly concave upward on log-log graph

By C. H. Li, S. H. Bendel, J. A. Sartell

Small additions of Mg to a 62 Cu-38 Ni alloy reduce the oxidation yates drastically between 500° and 850°C. Metallographic examination of the oxide scale disclosed a continuous network, identified by X-ray diffraction as (Ni, Mg)O solid solutions, located along the former alloy grain boundaries. This network acts as a barrier to diffusion. It is thought to result from the segregation of mg to the alloy grain boundaries. Oxidation occurs preferentially in these regions because of the presence of this more active element. To substantiate these ideas, homogenized alloys ulere oxidized. No oxide network was produced and oxidation rates were the same as for the parent 62 Cu-38 Ni alloy. An interesting aspect of the oxidation of the 62 Cu-38 Ni alloy previously described by the present authors, is the relationship between oxidation kinetics and changes in the microstructure of the oxide.&apos; The oxide scale on these alloys consists of two layers: a thin layer of CuO outside a much thicker layer of Cu20 in which particles of NiO are dispersed. In addition, immediately beneath the oxide-metal interface there is an internally oxidized zone in which NiO particles are dispersed in a copper-rich alloy matrix. Experimental observations on the growth of the individual layers under isothermal conditions have led to the conclusion that the rate-determining step in the oxidation of this alloy is the diffusion of anions through the Cu20 layer. It has been found that this diffusion rate is increased by the recrystallization of the outer CuO layer. Willardson and Gorton have reported that the addition of magnesium to the 62 Cu-38 Ni alloy causes a decrease in the oxidation rate.2 The present investigation was undertaken in order to determine the mechanism by which magnesium additions improve the oxidation resistance of the 62 Cu-38 Ni alloy. EXPERIMENTAL PROCEDURE Fifty-gram ingots of Cu-Ni alloys were prepared by vacuum melting materials of 99.999+ pct purity. These ingots were then remelted in a He atmos- phere and magnesium added in the form of a previously prepared and analyzed Cu-Mg master alloy. This procedure was necessitated by the high reactivity of pure magnesium at high temperatures. The proportion of the alloying elements was adjusted so as to obtain alloys with the desired magnesium content while maintaining a Cu-Ni ratio of 62 :38 by weight. (All alloy compositions will be given in weight percent.) Since it was considered desirable to study the oxidation of only single-phase, solid-solution alloys, it was necessary first to determine the approximate limit of solid solubility of magnesium in the 62 Cu-38 Ni alloy. A series of melts with various magnesium contents was made, homogenized by annealing for 48 hr at 975OC, then cooled in the furnace. Metallographic samples were prepared, etched in a potassium dichromate solution and examined at XI500 for the presence of a second phase. The alloys were then chemically analyzed and it was found that the approximate limit of solid solubility of magnesium in the 62 Cu-38 Ni alloy was 1 pct Mg. The ingots containing less than 1 pct Mg (0.19, 0.29, 0.58, and 0.87 pct Mg) were cold-rolled down to thicknesses suitable for oxidation experiments and specimens for gravimetric studies were cut to dimensions of 3.0 by 1.5 by 0.05 cm, specimens for inert marker studies, 0.6 by 0.6 by 0.17 cm. These were heat-treated in a protective atmosphere at a temperature (975" C) well above the oxidation temperatures and furnace-cooled to insure structural stability of the alloy during oxidation. Immediately prior to oxidation, the specimens were chemically cleaned by immersion for 1 min in a solution of 2 parts nitric acid, 1 part acetic

Jan 1, 1961