Search Documents

By C. D. Hall

This paper may be divided into two main parts: (1) analysis and data to show that Eq. 6 adequately correlates the authors' experimental data, and (2) the assertion that Eq. 20 "should permit scale-up with data taken from one pipe diameter for Newtonian, inelastic non-Newtonian, and viscoelastic fluids, such as dilute polymer fracturing fluids, which produce considerable pressure reductions." The authors might well recommend that this broad grouping of fluid types could be correlated by using Eq. 6 (requiring data from at least two pipe diameters). On the other hand, Eq. 20 is only shown to be adequate for conditions within the range of the present data set. Specifically, predictions by Eq. 20 for diameter scale-up to casing sizes would give larger deviations than implied in Table 2, and scale-up for some other fluids that give large drag reduction could be grossly in error. The description of a class of fluids "which produce a diameter family of straight parallel lines on a log (d?p/4L) vs log (8 v/d) plot. . ." defines a necessary condition for using Eq. 6; i.e., that s shall be constant for each fluid. The proposal that Eq. 20 is applicable to turbulent flow of such fluids is equivalent to saying that the spacing of these parallel lines, with A being constant, will be found to correspond to m = -1.20 for all such fluids. Table 1 shows the values of m derived from data on 16 non-Newtonian fluids. The key statement of this paper is that the diameter exponents "do not vary significantly from . . . -1.20." In statistical terms, this says that the m value for each fluid, which was obtained from all the measurements made on that fluid. has not been shown to be different from the arithmetic mean of the m values for the 16 fluids tested. This implies that the group mean, -1.20, is a more reliable value to use for an individual fluid than is the m value derived from the data on that fluid alone. No statistical test is described to support the key statement that the diameter exponents do not vary significantly from -1.20. Such a test would have to take into account that each individual m value is shown to correlate the pressure gradients for its fluid, with only about 3 percent standard deviation over the diameter range covered. That is, the m value is not the result of a single independent measurement. Thus, it would not be sufficient, or valid, to compare the deviation of one of the m values with the variance of the 16 values of m about their mean. An alternative interpretation of the key statement is that the diameter exponents do not vary appreciably from -1.20; i.e., that using -1.20 does not cause serious inaccuracy. Table 2 seems to be offered in support of such an assertion. However, such a statistical treatment (1) is unnecessary and (2) cannot lead to any valid scale-up generalization. The amount of error introduced by using a 1"wrong" value of diameter exponent — i.e., -1.20 instead of m — is determined by Eq. 6. The ratio of the pressures predicted by Eqs. 6 and 20 for scale-up from ds in small pipe to d, in casing is The pressure for 30-percent calcium carbonate slurry (m = -1.112) in 10-in. ID pipe would be predicted 30 percent higher by Eq. 6 than by Eq. 20, if based on tests in 0.5-in. ID pipe. The data set tested in Table 2 represents only a limited sample of pipe diameters (ratio of about 4:l) out of the population of diameters that would

Jan 1, 1970

By M. L. Hammond, G. M. Bowers

Silicon dioxide films have been prepared on silicon wafers by the oxidation of SiH, near 400°C. Results previously reported by Goldsmith and Kern7 have been qualitatively corrobovated using a system of different geometry, and significantly greater deposition rates have been achieved. The deposition rate is independent of thickness and nearly independent of the O2/SiH, ratio over a wide range; however, the rate is strongly temperature-dependent. The structure of the deposit depends upon the cleanliness of the substrate, and the electrical field strength of the deposited oxide depends upon the nature of the electrodes. MULTILEVEL interconnects in integrated circuits using aluminum metallization require an insulator that can be deposited at low temperature. Because the A1-Si eutectic is only 577°C, the maximum deposition temperature must be below about 500°C to prevent reaction of the aluminum with the silicon and with other materials in the integrated circuit. Many methods of depositing insulating films of SiO, and other dielectrics have been reported'-6 but the technique described by Goldsmith and Kern7 offers the simplest method of preparing SiO, films. The technique involves the reaction of SiH, and 0, to form SiO, films in the 300" to 500°C temperature range. Previous work using this gas system has been reported8 and the initial development work at Philco Microelectronics was done at the East Coast R&D Laboratories. 9 This paper corroborates some of the results of Goldsmith and Kern,7 shows that much higher deposition rates can be obtained than previously reported, and presents some observations on the structure of deposited SiO, films. EXPERIMENTAL CONDITIONS Fig. 1 is a schematic drawing of the apparatus used in this study. It differs from that of Goldsmith and Kern in that 1) the gases are mixed before they enter the reaction chamber, 2) there is no fritted glass diffuser, and 3) the wafers are stationary during deposition. The control thermocouple was in a well located 0.06-in. from the substrate, and three gas supplies were used: 1 pct SiH4 in N 2, O 2, and N 2. Water vapor must be excluded from the plumbing to prevent the formation of colloidal SiO, and the SiH4 must be diluted with an inert gas because SiH4-0, mixtures are extremely flammable above about 4 pct SiH The wafer temperature was determined by attaching a thermocouple to a test wafer and the temperature was maintained to +5 C. The conditions used in this investigation are listed in Table I. The temperature range was chosen for compatibility with the A1-Si system. The volume ratio of 0, to SiH4 is the actual volume of 0, delivered to that of SiH4, taking into account the initial dilution of the SiH4. This ratio, which was used by Goldsmith and Kern,7 is a useful indicator of the experimental conditions. The active/inactive ratio is the ratio of the volume of active gases, SiH4 and O2, delivered to the volume of dilution gas. The active/inactive ratio is an important parameter to specify for deposition rate determinations since it determines the amounts of reactants that arrive at the sample per unit time. Fairly high dilution ratios were necessary to achieve uniform deposits of SiO,. Most of the process characterizations were performed on the polished (111) face of 3 to 5 ohm-cm, n-type silicon wafers which were 0.009 in. thick; but the data approximately apply to deposition of SiO, on thin films of aluminum. It is difficult to obtain accurate deposition rate data over aluminum films because of coarsening the aluminum surface during exposure to 300" to 400°C. The deposited SiO, film fol-

Jan 1, 1969

By H. M. Manasevit, R. L. Nolder, L. A. Moudy

The observation that silicon-on-sapphire crystallo-graphic relationships exist which are unrelated to one another implies that there are regions of sapphire orientations containing these relationships which are separated by boundaries. The presence of these boundaries is demonstrated, initially on rounded cylindrical surfaces md more definitively on flat sapphire substrates with orientations along primary crystallograPhic zones. By means of reflection electron and X-ray diffraction techniques the crystallographic relationships have been studied. Two new relationships have been found, viz., (111) Si (1014) sapphire and (001)Si (1012) sapphire. Also found were multicrystalline relationships. The approximated boundaries between the crystallographic regions for each relationship are shown by means of a stereo-graphic projection. These studies also clarify discrepancies in the literature concerning the silicon-sapphire relationships on the (1120) sapphire substrate. SINCE the first report of the deposition of single-crystal silicon on sapphire in 1963' this epitaxial system has been the object of study by many semiconductor laboratories. In particular, there has been strong interest in the crystallographic relationships between the two materials. In 19652,3 we demonstrated that at least four distinct modes of orientation relationships can occur, namely, (100)Si (0112) sapphire, (111)Si 11 (1124) sapphire, (111)Si (1120) sapphire, and (111)Si (0001) sapphire. Each relationship could be obtained when the substrate plane was at the particular orientation indicated. Other workers as well as our laboratory concurred in general with the initial results. However, two discrepancies were observed. First it was observed by Joyce et a1.4 that on a (1120) substrate a (110) silicon deposit could be obtained. This was actually consistent with the (111)Si (1124) sapphire mode. Second, two other laboratories,5,6 as well as our own,7 have found a deposition on the basal plane which is different than our initial findings. Even though a (111)Si deposit was obtained, the crystallographic relationships were not identical to those previously found. Subsequent work in our laboratories has confirmed that the form found in these later investigations of the basal plane is by far the more predominant, and the original finding was an anomaly. The relationships which have been previously estab1.ished are summarized in Table I. It is seen that at present a rather small number of possible substrate orientations has been investigated. Therefore, it seemed appropriate that a systematic investigation of silicon epitaxy on various sapphire substrate orientations would be of value. It was expected that the answers to three questions could be obtained from such an investigation. These are as follows: 1) Since we know of at least four distinct modes of orientation relationships, over what range of substrate orientation do these modes occur and what is the form of the boundary expected between the different regions in which each mode occurs? 2) Are there other modes which occur on substrate orientations not previously investigated? 3) What is the explanation for the discrepancy between our laboratory's findings and those results reported by Joyce et al.4 concerning the relationships occurring on (1120) substrates? APPROACHES AND EXPERIMENTAL METHODS Two different approaches were used. Initially, depositions were carried out on rounded cylindrical surfaces which lay perpendicular to the (0112) plane.* By this means the presence of boundaries and their approximate positions could be determined. Finally and more definitively, depositions on flat sapphire substrates with a variety of orientations, mostly along primary crystallographic zones and about the (1120) plane, were studied by means of reflection electron diffraction and X-ray diffraction (using back-reflection Laue patterns and full circle goniometer techniques) in order to determine crystallographic relationships. The sapphire orientations investigated by these means

Jan 1, 1969

By Brown, R. W., G. H. Neill, R. G. Loper

All facets of ball sealer behavior must be known and understood to design for their optimum use in well treatments. The down-hole factors including the inertial forces, drag forces and holding forces have been discussed. Equations have been developed and presented for these forces. The use of these equations in this process will contribute to the most efficient use of an economical and proven ball sealer process. The urocedure presented has been used to successfully analyze field results indicating problems involving the use of ball sealers. These problems have been isolated and corrected in future treatments by using the procedure presented to revise the perforating program, the injection rate, the design of ball sealer stages and/or the type and size of ball sealers. The procedure presented also indicates that dislodgement of ball sealers occurs primarily in ultra slim-hole completions. INTRODUCTION The perforation sealing process has been proven highly successful and economical since it was introduced to the oil and gas industry in early 1956.' Since it has proved so successful in multistage fracture treatments at a much lower cost than conventional packers, the ball sealer process has become a byword in well completions.' Not only has it changed the design of fracture treatments, it has changed the concepts of selectively perforating and has greatly aided the success of the single-point and limited-entry techniques." Yet with all the advantages, there are still facets of ball sealer behavior that are not widely recognized. These must be known and understood to better design and execute the optimum ball sealer treatment. The down-hole behavior of ball sealers is influenced by a number of factors including an "inertial force", a "drag force" and a "holding force". These are considered in two phases—the initial seating of the ball and the forces tending to unseat the ball after contact has been made (see Fig. 1). This paper presents a theoretical study of these factors and formulas for predicting their influence. The authors acknowledge the fact that the flow equations presented herein describe the behavior of Newtonian fluids and do not cover the so-called power-law fluids. The assumptions made, however, are reasonable and well within practical limitations of the conditions encountered during field application. For example, it is assumed that the perforations are round with a discharge coefficient of 0.82. While this may not be exactly true, it has been verified within practical field limitations.' The discharge coefficient would not vary significantly even if the perforation size or configuration varied pronouncedly. The approach of assuming uniform perforation size is not as strained as it might first appear. Normally, a perforator contains charges equally spaced with the same number in each vertical (or horizontal) plane. If the perforating gun is not centralized, the holes in one plane will be smaller than average, while the holes in opposing planes will tend to be larger than average. The net result is that the two effects tend to balance. FACTORS AFFECTING CONTACT WITH PERFORATlON Ball sealers are injected into well treating fluids for the purpose of contacting and sealing those perforations which are accepting the fluid flow. The efficiency of the sealers is primarily influenced by: (1) the velocity of the balls down the pipe, and (2) the fluid velocity through the perforations. To divert the sealer to the perforation, the inertial force of the ball must be overcome by the drag force created by the fluid velocity through the perforation. BALL VELOCITY The final, or stable, velocity of the ball is the sum of the actual fall velocity and the fluid velocity. This will be

By A. Gittens

In two recent papers stevens13,14 has compared the methods used by various investigators to calculate the contribution of grain boundary sliding to the total creep strain. In previous work Gittins15 has shown that the strain due to grain boundary sliding Egb may be estimated from the lateral, d, and vertical, v, displacements of a longitudinal marker line as it crosses a grain boundary. He obtained the formula, where nL is the number of grains per unit length measured on the specimens before creep, and \p and ?2 are, respectively, the angle between the grain boundary and the stress axis in the plane of the surface and normal to the surface (see Ref. 13, Fig. 1). The subscript L indicates that observations are made in a longitudinal direction parallel to the stress axis. It was generally found more convenient to measure the grain size of a specimen before creep; this was necessary to characterize the specimen. Stevens13,14 however prefers to define Egb in terms of NL, the number of grains per unit length after creep. Since nL - NL(l + €1) where cr is the total strain, Eq. 1 11 becomes The quantities d and ?, can be determined directly from the marker lines, while v and ?2 are found by sectioning normal to the surface. In the interior when symmetry considerations suggest that v = d and ?1 - ?2 the strain due to grain boundary sliding may be estimated from a suitable internal longitudinal marker line by the equation To determine cgb at the surface using either Eqs. |1| or [2] requires a large number of observations, particularly if many estimates of cgb are required, and also it is impossible to determine the variation of cgh with creep strain on the same specimen because the specimen must be sectioned to determine ?2. Gittins15 found that for a number of materials the sliding vectors d and v were related to each other so that provided one was measured the other could be estimated. Also during creep, grain boundaries tend to migrate to meet the surface at angles close to 90 deg,16 so that the term nL(v/tan?2)L was comparatively small when 4,- 90 deg. By determining the frequency distributions of ?2 on control specimens this term can be estimated without sectioning the creep specimen. In this way an empirical equation based on Eqs. |1| and |2| enables cgb to be determined with reasonable accuracy with fewer observations of the angles and sliding components. For B brass, the relationship between the sliding vectors indicated that ceb = 0.85nLVL . In Table I of Stevens paper this equation is stated to be incorrectly derived. However it should be pointed out that this is an empirical equation derived from Eq. |1|, based on a knowledge of the relationship between d and v and on the frequency distributions of ?1 and ?2 for a particular material under given experimental conditions. It was used only when these relationships were known, to reduce the number of observations required to estimate cb. Therefore the comment "derivation incorrect does not apply. Gittins did in fact make use of Eqs. |1| and 131. Author's Reply I must apologize to Dr. Gittins for the inadvertent misrepresentation of his work. Other slight errors in Table I have been pointed out to me. Langdon (No. 12 Table I), Gifkins (No. 7, Table I), and possibly other workers did not measure the number of boundaries per unit length after creep (A') but the corresponding quantity bejove creep (n). As Dr. Gittins points out, a completely analogous set of equations to those given in the paper can be derived by the same method: using the number of boundaries per unit length before instead of after creep. The basic equations corresponding to Eqs. [2] and 171 of the paper, from which all the others can be derived, are

Jan 1, 1968

By T. M. Kegley, J. H. Frye, D. L. McElroy, M. L. Picklesimer, E. E. Stansbury

Measurements of rate of growth, thermodymmic quantities, and partitioning of Mn are reported for high-purity eutectoid Fe-C and Fe-C-Mn steels for the auistenite-pearlite reaction. Evaluztion of the cotnstants of a derived-rate equation with these data indicate that Mn additions cause the decrease observed by increasivg the activation energy and apparent activation entropy or the movenment of complexes at the interface. A diffusion-controlled mechanism is not required. Appreciable partitioning of Mn to cementite occurs entirely behind the austenite-pearlite interface for more than 30°C of subcooling. ThE growth of pearlite from austenite has been extensively studied over a period of years. but the rate-controlling step in the process is not known nor is the mechanisrn of the effect of alloying elements understood. A comprehensive review of the subject has been given by Mehl and Hagel. The present study has been concerned with further considerations of the effects of alloying elements on the rate of growth of pearlite. Analytical relatioriships for the rate of growth of pearlite have been considered by Zener.2 Turnbull,3 Machlin,4 Cahn.5 Frye,6 and Frye et a1.7 The derived expressions in each case involve the free energy of the transformation. the interlamellar spacing. and some rate term of the nature of the diffusion coefficient. mobility or frequency of occurrence of some rate-determining step. The treatments. of course. differ in detail. Further consideration is given in the present paper to the equation of Frye, Stansbury. and McElroy which was derived from modified absolute rate theory without in any way specifying the mechanism of the reaction. This equation is r=c ? T ? F exp -?E*/RT [1] where r = rate of growth in mm/sec c = a constant AT = equilibrium transformation temperature minus actual transformation temperature "C ?F = free energy per mol of austenite minus that of pearlite ?E* = activation energy R= gas constant ?F in Eq. [1] is computed from the following: ?F=?H?T/To - ?T' dT' ?T' ?CdT"/T" [2] To To where AH is the enthalpy per mol of austenite minus that of pearlite at the equilibrium temperature, To, and AC is the difference in specific heat between austenite and pearlite per mol. Subsequently, Darken, Fisher, and Galligan8 have reported that a layer of ferrite separates pearlite from austenite at the growing interface (i,e., the advancing face of Fe3C is in contact with ferrite, not austenite) and the ratio of the thickness of this layer to the spacing of the pearlite has not been found to vary with alloy content or transformation temperature. This observation does not require any alteration in the above equation: however, it may require that a different significance be attached to the factor, AT. The AT term appears in Eq. [I] as a consequence of a factor for the amount of growth per complex transformed at the interface.* It is based *This is not the sole reason for expecting such a relationship. See Ref. 7._______________________________________________________ on a l/2 ?T dependence of the interlamellar spacing. Zener2 has developed a theory which leads to Soa 1/?V where So = the pearlite interlamellar spacing. This relation is as consistent with the data on interlamellar spacing as other relationships, although, as pointed out by Mehl and Hagel,1 the data are not sufficiently precise to consider the relationship as confirmed. The work of Darken, Fisher, and Galligan8 requires that the above be written where f = the thickness of ferrite layer between austenite and pearlite. The present paper is concerned with iron-carbon-manganese eutectoid alloys and reports experimental measurements of the following kinds: 1) the effect of

Jan 1, 1961

By E. A. Hopkin

Laboratory tests conducted by the author. together with actual field experience in Canada. have indicated the magnitudc of some of the factors affecting ability of drilling mud to clean the hole. A correlation was observed between funnel viscosity and particle slip velocity. A relationship was also observed between the Bingham yield value of the mud and the particle .slip velocity. Laburatury tests indicated little change in slip velocity with funnel viscosities up to ahout 80 sec/qt or to yield values of aboul 15 lb/100 .sq ft. The lip velocity of the partic1e.s. studied vuried from 160 ft/min with water to neerly zero for high mud viscositie. lncreasing the mud density, creating laminar annular mud flow or rotating the drill pipe Mlay also improve the carrying capacity of mud. Field tests indicated that during fast upper-hole drilling, the rate of cutting removal must he sufficien to maintain the concentration of cuttings in the annulus at less than 5 percent by volume to prevernt balled-up drill collars and .stuck drill pipe. INTRODUCTION In the Mid-Continent area, hole cleaning during rotary drilling normally involves the use of low-density drilling muds with funnel viscosities varying from 26 to 80 sec/qt. However, frequently it is necessary to raise the mud viscosity in excess of 200 sec/qt to clean the hole. In addition, field experience has indicated that during very fast upper-hole drilling, annular mud velocities considerably in excess of the maximum slip velocity of the cuttings (normally considered to be 120 to 150 ft/min) are necessary to prevent balling the bit and drill collars. In view of these variations in mud properties and circulation rates necessary to clean the wellbore, a study of factors affecting hole cleaning and the relative magnitude of each was undertaken. The study was also prompted by the knowledge that a drilling mud with properties close to those of water would drill at rates several times that of a viscous mud. The carrying capacity of these near-water muds was not clearly understood, however, and many drillers were reluctant to use them. DISCUSSION EFFECT OF MUD PROPERTIES ON DRILLING RATE The effect of mud properties on drilling rate has been studied by a number of investigators. 1t is clear from their work that mud properties close to those of water will significantly improve drilling rate and bit life. If the water-type mud will also effectively clean the hole and cause no logging or formation impairment problems, the choice of an appropriate drilling fluid is clear—use water. Often, however, a water-type drilling fluid is not capable of effectively cleaning the hole. THEORY OF PARTICLE SETTLING The forces acting on a particle that is allowed to settle through a fluid quickly reach equilibrium, and the particle settles at some constant velocity referred to as the slip velocity. The slip velocity of any particle depends upon a number of factors such as density and viscosity of the fluid, the volume, specific gravity, shape and roughness of the particle, and upon the shape and area of the projected face of the particle at right angles to the direction of relative fluid solid movement. Brown' expressed this relationship for spheres in Newtonian fluids in the following equation. where v, = slip velocity D, — diameter of sphere p, = density of spheye p, = density of liquid /a = friction factor g = constant. The friction factor f, varies as the size, shape, surface roughness and slip velocity of the body, and as the density and viscosity of the fluid varies. WaddelV ound experimentally that where dynamic similarity exists. f, was a function of the particle Reynolds number (NRe) of the sphere (Fig. 1). Fig. 1 shows that f, for spheres varies from over 10 for low NKe numbers to 0.44 for N,,. over 600.

By S. R. B. Cooke, H. S. Choi, D. H. Harraway, I. Iwasaki

The mineralogy of natural iron ore samples and their slimes have been investigated and the nature of slime interference in iron ore flotation has been studied. The results have been correlated with the flotation behavior of a natural ore sample. Various conclusions have been reached as a result of these studies. In the flotation of iron oxide ores, desliming is generally considered to be an essential step to successful separation, irrespective of whether the iron oxide or siliceous gangue minerals are to be floated, although some workers have found empirically that certain collectors and reagent combinations are less sensitive to the presence of slime than others. The detrimental effect of slime is considered to be twofold. First, the presence of slime presupposes a high reagent consumption because of high specific surface. When clays are present and when cationic collectors are used abstraction by ionic exchange may augment this condition. Attempts to compensate by the addition of large amounts of reagent result in only partial correction. Second, slime coatings on granular particles tend to interfere with the bubble-mineral contact. The nature of slime coatings in flotation is quite complex and has been the subject of considerable study; an attempt was recently made by Fuerstenau et al.' to relate the slime-coating density with flotation recovery and with the properties of the electrical double layer at the mineral surfaces. Although major constituent minerals in the oxidized iron ores of the Mesabi and Cuyuna ranges are hematite and goethite (limonite) as the iron oxides, and quartz and chert as the siliceous gangue, very little is known of the mineralogy and of the chemical and physical characteristics of the slimes, and still less of the interaction between slimes, collector, and the mineral surfaces. For pragmatic reasons iron ore slimes are defined as — 20p material, based upon the calculated settling velocity of a quartz sphere of the stated diameter, and of all other equal-settling mineral particles. Active slime interference with flotation is usually attributed to particles 5p in size or finer, the extent of such interference appearing to depend not only on the mineralogy but also on the size distribution of each mineral in the slime fraction. The objective of this work was to acquire a more specific knowledge of the constitution of iron ore slimes and of their role in flotation. The mineralogi-cal composition of natural slimes was obtained by differential thernal analysis (hereinafter referred to as DTA), electron microscopy, and chemical and X-ray analyses. Some flotation tests were made on natural ore samples. The surface properties of several constituent minerals were investigated by electrophoretic measurements and by determination of the adsorption of some typical flotation collectors. Finally, simplified flotation tests were made on granular goethite and quartz in the presence of various artifically prepared slimes; the results were compared with laboratory flotation tests made on natural ores. MATERIAL AND METHODS Mineralogy of Some Iron Ore Slimes: Samples of two wash-ore tailing::, hereinafter identified as A* and B, *Sample A of the present investigation 18 the same am ore A in ref. 3. of a wash ore C, and of a paint rock D, all from the Mesabi range, were selected for investigation on the basis of prior information concerning their flotation behavior. The two wash-ore tailings, when deslimed, were known to respond well to flotation, whereas the wash ore was high in goethite, tended to be slimy, and responded poorly to flotation. The paint rock also exhibited sliming tendencies, and was suspected to contain relatively large quantities of clay minerals because of its high A1203 content. Gruner states that the paint rock of the Mesabi frequently contains kaolinite. To ascertain the relative sliming tendencies of each sample a standardized desliming method was employed. All samples were ground through 65 mesh using jaw and rolls crushers. A 500-g sample of each

Jan 1, 1962

By J. L. Lytton, T. E. Tietz

The recovery of flour stress of recrylstallized are-cast molyhdenum after 10 pet prestrain at 150°C was studied for recovery temperatures of 725o 801o and 901oC. The yield phenomenon observed after recovery was such that the yield drop diminished and disappeared as recovery proceeded. A fractional flow-stress recovery parameter fr, was developed and it was shown that, independent of recovery ternperatzwe, equal f, values of about 0.16 produced superposition of stress-strain curves except for the initial yield phenomenon. The activation energy for low degrees oj recovery was found to he 32,600 cal per mole, in close agreement with a reported valzie for difjzdsion of carbon in molybdenum. For values of fr greater than about 0.15 where recovery was thought to predominate over strain aging, the activation energy varied from 43,600 to 82,500 cal per mole and increased with increasing recovery and recovery temperature. A 10 pet pre-strain shifted the ductile-brittle transition temperature, as indicated by tensile elongation, from 10°C. to about -22"C. A 100-hr recovery anneal at 801oC was shown to shift the transition temperature to about —5oC, the major shift occurring during the first hour of recovery. ThE recovery of cold-worked metals is accompanied by changes in the arrangement of the dislocations generated during cold work. In recent years, the use of transmission electron microscopy has allowed direct observation of some of the dislocation configurations formed, such as "tangles", networks, cells, subboundaries, and loops. These observations have recently been reviewed by Keh and Weissmann.' Accordingly, some understanding has been gained concerning the nature of the dislocation rearrangements involved. The relationship between recovery and the decrease of mechanical strength, however, has only recently been explored. Previous efforts of the authors were directed toward investigating the decrease of tensile flow stress of aluminum and Al-1 pet Mg during recovery2 and the effects of concurrent elastic and plastic strain on the flow-stress recovery of aluminum.3 In addition, an investigation of the effect of four different solute additions on the tensile flow-stress recovery of aluminum has been made.4 Many important metals, however, have the bee crystal structure, e.g., iron, molybdenum, tungsten, vanadium, columbium, tantalum, and chromium. These metals might be expected to behave differently during recovery from those with other crystal structures for two reasons. First, the previous and present investigations of aluminum and its alloys have shown that tensile flow-stress recovery is strongly affected by foreign atoms in the lattice. In the bee metals, interstitial impurities are usually present and might be expected to have various effects upon both flow behavior and recovery. For example, heterogeneous yielding might occur during prestraining. The diffusion of these interstitials can also relock dislocations during recovery so that a yield point might be obtained upon restraining (strain aging). Second, some of the bee metals (e.g., iron, molybdenum, and tungsten) undergo the familiar structure-sensitive ductility transition and exhibit brittle fracture at suitably low temperatures. It is of considerable interest to determine the effects of recovery on the transition temperature of such metals. The purpose of the present investigation was to extend our knowledge of flow-stress recovery to the bee metals and to study the effects of recovery on the ductile-brittle transition temperature. Unalloyed recrystallized are-cast molybdenum was chosen as a test material, since it exhibits a well-defined ductility transition near room temperature. As in the previous studies, the degree of recovery was measured in terms of the decrease of tensile flow stress. Prestraining and restraining were carried out at 150°C to avoid the possibility of brittle fracture. EXPERIMENTAL PROCEDURE The are-cast molybdenum sheet used in this investigation was obtained from the Climax Molybdenum Co. of Michigan. It was 0.060 in. thick by 10

Jan 1, 1964

By Frederick C. Pearson, Albert Brown, Emil R. Ermert

EARLY in 1946 the Shen-Penn Production Co., a subsidiary of the Philadelphia and Reading Coal and Iron Co., was organized to operate the Shenandoah Stripping, one of the largest open-pit anthracite mines in the world. This mine is located in the Western Middle Field near Shenandoah, Schuylkill County, Pa., and involves the removal of approximately 60,400,000 cu yd of overburden to permit the recovery of 9,200,000 gross tons of rough cleaned coal from an area 8000 ft long and 2700 ft wide, with a maximum depth of 400 ft. Operations began in August 1946, and within the next few years 31,811,000 cu yd of overburden and 2,994,000 gross tons of rough cleaned coal had been removed. During the early days of open-pit mining operations, steam was the primary motive power. However, as the size of the equipment increased, so did the problems of fuel, water, labor, and maintenance. Because of its efficiency and convenience, electric power rapidly replaced steam. Later the modern diesel engine was developed, and today practically all stripping equipment is either electric, diesel-electric, or diesel-powered. In locations where electric power is not readily available or where operations are scattered over a large area, diesel power is favored. At Shenandoah operations are highly concentrated in a relatively small area, and this fact, coupled with the availability of ample electric energy, indicated the selection of electric-powered equipment. Since this was a new company, it was necessary to purchase every item, from the smallest tools to stripping shovels. Although the organization period was hampered by the postwar scarcity of materials, an excellent opportunity was afforded to set up a thoroughly modern plant. This paper will describe the power facilities, with particular emphasis on the more unusual features. The substation is a major item in a power distribution system, and a simple arrangement utilizing high quality equipment is an excellent investment. With this in mind, a 3000-kva unit type substation was selected and installed. This unit is made up of three sections, which are assembled simply by bolting together the various sections and connecting the power leads, see Fig. 1. The transformer proper, rated at 3000 kva, is oil-filled and nitrogen-sealed, having a 66,000-v Delta primary and 4160-v Wye secondary with a grounded neutral and 46-ohm ground resistor. The function of this resistor will be described later. Electrical energy is purchased from a public utility and is delivered to the main substation at 66,000 v over one or the other of two incoming lines. These lines terminate on a high voltage switching section on which high voltage fuses and Thyrite lightning arresters are installed. The outgoing feeder switching unit is made up of the following: 1-An auxiliary compartment consisting of potential and current transformers with necessary terminal blocks for connecting to power company metering equipment. 2-Another auxiliary compartment containing a control transformer, battery, and battery charger for the furnishing of power for tripping and reclosing of power circuit breakers. The assemblies in these auxiliary compartments are mounted on carriages that can be rolled in and out with ease. Potential transformers are readily accessible for inspection and fuse replacement. When withdrawn from the operating position they are completely de-energized and grounded. 3-Two feeder breaker compartments, each with a 1200-amp, 5000-v, solenoid-operated magneblast air circuit breaker, voltmeter, ammeter, overcurrent relays, and a ground relay. These power circuit breakers can be quickly and safely withdrawn or repaired. A built-in isolating mechanism is used for lowering the breaker to a transfer truck for repair or inspection or as a means of disconnecting the breaker from the source of power as an added precaution when working on power lines. Thus inspec-

Jan 1, 1952

By H. D. Kessler, R. J. Van Thyne

Phase diagrams of the titanium-rich portion of the ternary systems from 0 to 10 wt pct Al and 0 to 1 wt pct 0, N, and C were determined. Micrographic analysis of annealed high purity arc melted alloys was the principal method of investigation and was supplemented by X-ray diffraction. TITANIUM-ALUMINUM alloys exhibit excellent properties, particularly for elevated temperature use. For this reason, the Materials Laboratory, Wright Air Development Center, sponsored an investigation of the effect of the interstitially soluble contaminants, oxygen! nitrogen, and carbon, on the Ti-A1 system. Using high purity arc melted alloys and micrographic analysis of annealed samples as the chief method of investigation, titanium-rich partial phase diagrams were determined for the ternary systems. Experimental Procedure Materials: The titanium used in the preparation of the alloys was iodide crystal bar (99.9+ pct pure) produced by the New Jersey Zinc Co. and Foote Mineral Co. The aluminum was obtained from the Aluminum Co. of America in the form of sheet. The given analysis was: Si, 0.0006 pct; Fe, 0.0005; Cu, 0.0022; Mg, 0.0003; Ca, <0.0006; Na, <0.0005; and Al, 99.99. High purity titanium dioxide purchased from the National Lead Co. was used in the preparation of the oxygen-bearing alloys. The spectrographic analysis of this material was as follows: SiO,, 0.07 pct; Fe,O,,, 0.002; A1,0,, <0.001; Sb,O,, <0.002; SnO,,, <0.001; Mg, <::0.001; Cb, <0.01; Cu, 0.0004; Pb, 0.002; Mn, <0.00005; W, <0.01; V, <0.002; Cr, <0.002; Ni, <0.001; and Mo, <0.002. Special spectroscopic graphite rod purchased from the National Carbon Co. was used in the preparation of the Ti-A1-C alloys. Alloy Preparation: Alloy ingots weighing 10 grams were melted in a nonconsumable tungsten electrode arc melting furnace using water-cooled copper hearths and a helium atmosphere. The arc was struck on a tungsten stud in the copper block to minimize contamination from the electrode. Details of the techniques have been reported.&apos;, &apos; Each ingot was inverted and remelted a minimum of four times to insure homogeneity, without opening the furnace. Remelting small control ingots of iodide titanium under similar conditions resulted in no measurable hardness increase, indicating contamination-free melting conditions were used. The alloy charges and resultant ingots were weighed to the nearest milligram. Only small weight losses were obtained upon melting, indicating that the actual compositions were very close to the nominal compositions. Check analyses substantiated this. Oxygen and nitrogen were added during ingot preparation as master alloys. An oxygen master alloy containing 25 wt pct 0 was prepared by melting iodide titanium and pressed titanium dioxide (40 pct 0) powder. The Ti-25 pct 0 alloy is more easily handled during weighing and charging than the pressed TiO, powder. Chemical analysis indicated that titanium-nitride received from an outside source contained large amounts of impurities (2.7 pct Ca). Therefore, a master alloy containing approximately 12 pct N was prepared by melting nitrided sponge titanium. As the sponge melted with difficulty, the alloy was diluted by addition of titanium to lower the melting point. The chemical analysis of the final alloy was 6.69 pct N. Both master alloys are friable and were used as —8 +16 mesh lumps resulting from crushing the ingots. NIelting the Ti-A1-0 and Ti-A1-N compositions five times produced homogeneous ingots. The preparation of homogeneous carbon-bearing alloys caused considerable difficulty. Alloys containing less than 1 pct C were prepared using spectrographic graphite rod, 1/8 in. diameter by 1/8 in. long. However, the graphite did not dissolve even with long melting times and many remelts. The ?/s in. pieces of graphite were used rather than powder because the low density powder is troublesome in arc melting.

Jan 1, 1955

By Tomek Ulatowski

This presentation discusses the necessary ingredients for the creation of an acceptable credit structure, i.e., the structure that enables the borrower to attract the funding from different groups of lenders. As the amount of debt in the capitalization of new projects increases, it becomes quite critical that the mining venture sponsors obtain funding in a manner which maximizes the project economics. Attention normally is focused on the inherent operating cost characteristics of a venture such as the cost of mining, milling, transportation, smelting, and refining. The projected costs are derived from lengthy feasibility studies often taking millions of dollars and years to prepare. For a variety of reasons as the amount of equity directly provided by the sponsors declines, the balance must be made up from external credit sources. As a result, once the project becomes operational, debt service requirements assume crucial importance. For most new ventures, in the early years of production, the total amount of periodic debt amortization charges together with the projected interest expenses equals or exceeds the total operating cost. Utmost care must be taken, therefore, to reducing these costs which will then increase the chances for the implementation of the project. Project sponsors must carefully evaluate different funding sources in terms of their potential economic impact and, concurrently, evaluate alternative credit structures in order to be able to attract the selected funding vehicles. In this presentation I will discuss credit sources in the order of their relative appeal; i.e., from the standpoint of maximizing rates of return to the project owners. The tenor of the credit, that is the period from signing a loan agreement to the time of last repayment under the loan agreement, is undoubtedly the most important variable in evaluating the attractiveness of debt sources. Therefore, I will begin by discussing debt sources with the longest maturities and will end with a description of credits with relatively shorter tenors. The availability of different funding possibilities will naturally depend on the geographic location of the proposed venture, the creditworthiness of the borrower and the sponsors, the project economics, and finally the actual liquidity conditions in the capital markets at the time that the funds are committed. U.S. LONG-TERM DEBT MARKET - PUBLIC OFFERINGS AND PRIVATE PLACEMENTS The long-lived nature of most mining projects dictates that the majority of the mine financing be raised from long-term sources. In today&apos;s world capital markets U.S. debt issues probably carry the longest maturities available. The maturities can range from 5 to 40 years. Various types of long-term corporate debt instruments can be differentiated from one another by their terms, seniority as to claims against the assets of the issuer, legal form, and whether issued privately or through a public offering. The terms of these debt issues vary widely. For example, the interest rate which is typically fixed at the time of issuance, is determined by the credit quality of the issuer, the maturity, the average life, the seniority of the debt, the quality of control covenants or type of assets pledged as security, and other exogenous factors such as: economic policy, general level of interest rates and cur-rent level of financing activity. The entire principal amount of the debt may be due and payable at maturity (known as a "bullet" maturity) or, more typically, through periodic repayments in equal or staggered amounts. Almost all such debt issues allow for prepayment at the option of the issuer, but require varying amounts of repayment penalties. Disbursements are made over a relatively short period after signing of the debt documentation which may not coincide with the actual funding needs of the project. Financial covenants can be quite stringent, covering areas such as the issuer&apos;s financial leverage, working capital, payment of dividends and other distributions, merger and acquisition activities, and limitations on the disposition of property. Long-term debt issues may also be distinguished by the seniority of their claims against

Jan 1, 1985

By H. M. Strong

Nickel and carbon form a metastable nickel-nickel carbide eutectic system which was ohsen)able by freezing latent-heat arvests at pressures 210 khars. The eutectic freezing temperature Is pressure had a slope of -0.9°Cper kbar. At 50 khars, the nickel-nickel carbide eutectic was at 1297° * 5°C and 0.34 mole fraction Ni3C. The composition of the eutectic increased at an initial rate of 0.004 mol fraction per kbar. The carbon content of the a phase also increased with pressure. Seveval illustrations of nickel carbide microstruc-tures are given and data on the lattice expansion of nickel in terms of carbon content are included. IRON, cobalt, and nickel form a variety of carbides with compositions varying between MeC and MC. Of the three, the carbides of iron are most easily recovered. The nickel carbides are so unstable at temperatures above 900°C that their occurrence in Ni-C alloys is unusual. Under sufficient pressure, however, the temperature range for the stability of the ferrous metal carbides is likely to be increased, or their decomposition suppressed (e.g., see Hilliard&apos; on FesC), providing an opportunity to study the unstable carbides more closely. This was found to be true for nickel carbide which, under pressure, formed a eutectic with nickel. This paper describes the conditions of formation of nickel-nickel carbide eutectic compositions and their subsequent recovery. Despite the comparative instability of the nickel carbides, something is already known about them through various chemical reactions at temperatures not exceeding a few hundred degrees, and from their recovery by the very rapid quenching of molten Ni-C alloys from temperatures of about 1450°C or above. At least four carbides of nickel have been reported: A discussion of the state of knowledge about Ni3C is given by ansen.&apos; The studies of Morrogh and williams7 who obtained Widmanstatten nickel carbide structures and of sidorenko8 who first obtained nickel carbide eutectic structures are also of interest. Sidorenko used short heating times for molten nickel in a graphite crucible followed by very rapid quenching to obtain the eutectic structures. Long heating times produced only the nickel-graphite eutectic probably due to the presence of graphite crystals in the melt which nucleated graphite crystallization. According to the early work summarized by anssen,&apos; nickel carbide may form at temperatures far above the melting point and may be recovered by very rapid quenching. More recently, Giardini and Tydings, working on diamond synthesis, found evidence for "free nickel carbide&apos;&apos; in high-pressure reactions carried to temperatures far in excess of the limit for diamond growth, followed by rapid thermal quench. They stated that nickel carbide was not found in association with diamond growth. In the present work, observations were made on the melting and solidification temperatures of Ni-C alloys as a function of composition, temperature, and pressure. The nickel-nickel carbide eutectic was found at pressures above 10 kbar and its variation with pressure followed up to 60 kbar. Nickel-nickel carbide microstructures were recovered by quenching under pressure from temperatures of about 1400°C and higher. The formation of nickel carbide structures under ordinary diamond-growing conditions was usually observed but the role of carbide formation as an intermediate step in diamond growth needs more study. EXPERIMENTAL METHODS The melting and solidification temperatures of various nickel-graphite compositions were observed by latent-heat arrests in the sample holder illustrated in Fig. 1. The samples were composed of mixtures of carbonyl nickel powder and National Carbon Co.&apos;s SP-1 graphite packed to densities of about 6.4 g per cu cm. They were contained within

Jan 1, 1965

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 Hiroshi Kato, Shin’ichi Nagashima, Hiroshi Takechi

The preferred orientations detleloped during cold rolling, annealing, and hot rolling of low-carbon steel sheets have been investigated by means of an inverse pole figure technique. And X-ray line profile meas-urerrlents were made for the study of the factor which causes the textural changes between cold-rolled and recrystallized sheets. From ND (sheet normal divection) and RD (rolling direction) inverse Pole figures tlze axis detisity change with increasing cold reduction was measured and the Preferred orientations were analyzed. These results can be explained by crystal rotation due to slip during rolling deformation. The ND axis density change during annealing of cold-rollitig sheets is contrary to the change observed in the case of rolling. It was found that the orientation dependency of stored strain energy resulting from cold rolling apparently corresponds to the textural changes during recrystallization. The texture in hot-rolled sheets was found to be influenced by finishing tenlperatltre. In a sheet kol-rolled below the Ar, temperature a significant differernce in texture between the surface and central layers was obtained, the former having {hkl}(111} and {110 )(uvw) orientations and the latter having {pqr}( 110) orientation. A large number of investigations1 6 have been made on the textures of iron and low-carbon steel sheets, and, since it was shown7 13 that preferred orientations are closely related to deep drawability of steel sheets, studies14" 7 on the influence of production processes upon textures have been made from the standpoint of industrial interest. Generally, the textures in low-carbon steel sheets have several components of preferred orientations, but the overlapping of poles makes it difficult to analyze texture quantitatively from pole figures. On the other hand, an inverse pole figure directly expresses quantitative values of components of preferred orientations as the density of a certain sample axis with respect to the crystallite coordinates. Such a texture representation is desirable for a study of deep drawability. The purpose of the present work is to use both inverse and conventional pole figures, with emphasis on the former, to clar~fy the textural changes ~n low-carbon steel sheets by the production processes of hot rolling, cold rolling, and annealing. And the factor causing the textural changes during annealing is studied using an X-ray diffraction line profile technique. EXPERIMENTAL PROCEDURE Sheet specimens for study of cold rolling and annealing textures were taken from a hot-rolled strip of rimmed steel 2.7 mm thick. The chemical composition is shown in Table I. Although the study of cold-rolling textures requires a suitable material with random orientation, it is difficult to produce a perfect random one. But the used specimens, obtained from the hot-rolled strip with a finishing temperature higher than Ar3, were nearly random, such as shown in Fig. 15. These were cold-rolled to various reductions by a 4-high rolling mill with 150-mm-diam work rolls and 300-mm-diam backup rolls. The influence of hot-rolling conditions such as finishing and coiling temperatures on textures of hot-rolled strips was studied. The chemical compositions of these hot coils are given in Table I1 and the hot-rolling conditions are given in Table 111. inverse Pole ~i~ure -~easurement. Two methods for determining inverse pole figures with high precision have been proposed. The first, by Jetter et a1.,18 consists of obtaining mean densities of divided regions in a unit stereographic triangle by solving simultaneous linear equations to approximate the integral relation between conventional and inverse pole figures. The second, proposed by ~oe&apos;~-&apos;~ and Bunge,22&apos;23 consists of determining an expanded series of spherical surface harmonics which express axis densities from several sets of pole distribution data. Both methods require laborious numerical calculations which make them unsuitable for practical use. However, the axis densities of the locations corresponding to low Miller indices in a unit stereographic triangle can be obtained

Jan 1, 1969

By David E. Nicholas

INTRODUCTION In this paper, a numerical process for selecting a mining method, with the emphasis on underground mass mining techniques, such as caving, induced caving, and stoping, is proposed. In the past, selection of a mining method for a new property was based primarily on operating experience at similar type deposits and on methods already in use in the district of the deposit. Then, the chosen method was modified during the early years of mining as ground conditions and ore character were bet- ter understood. Today, however, the large capital investment required to open a new mine or change an existing mining system make it imperative that the mining methods examined during the feasibility studies and the method actually selected have a high probability of attaining the projected production rates. Although experience and engineering judgment still provide major input into the selection of a mining method, subtle differences in the characteristics of each deposit, which may affect the method chosen or the mine design, can usually be perceived only through analysis of measured characteristics. The parameters that must be examined when choosing a mining method include: 1) geometry and grade distribution of the deposit; 2) rock mass strength for the ore zone, the hanging wall, and the footwall; 3) mining costs and capitalization requirements; 4) mining rate; 5) type and availability of labor; 6) environmental concerns; and 7) other site-specific considerations. This paper encompasses a detailed look at the first two parameters since they, plus mining costs, have the greatest impact on the selection of a mining method. The proposed method selection process is for a project where drilling has defined sufficient geologic reserves, but little or no underground development has been done. Since each deposit has its own characteristic geanetry/grade distribution, and rock mechanics properties, mining method selection should be at least a two-stage process. In Stage 1, the deposit is described in terms of geometry, grade distribution, and rock mechanics properties. Using these parameters, the mining methods can be ranked to determine which are most applicable; they can then be considered in general terms of mining and capitalization cost, mining rate, type and availability of personnel, environmental concerns, and other site-specific considerations. In Stage 2, the most likely mining methods are costed out, based on a general mine plan. Mining and capitalization costs are used to determine a cut-off grade fran which a minable reserve can be calculated; economic comparisons can then be made to determine the optimum mining method and economic feasibility. During the mine planning phase of Stage 2, rock mechanics information would be used to provide realistic estimates of underground opening size, amount of support, orientation of openings, and caving characteristics, and open pit slope angles. If ground control or operational problems should be encountered with the methods being considered, modifications could be made. Although planning on paper extends start-up time, it is cheaper to err on paper than to find the error after mining has begun. METHOD SELECTION - STAGE 1 The main purpose of Stage 1 is to select those mining methods which should be considered in greater detail. The simplest way to do this is by defining those characteristics required for each mining method and then determining whether the characteristics of the deposit are suitable. However, no one mining method is so restrictive that it can be used for only one set of characteristics, as indicated by the classification system proposed by Boshkov and Wright (1973). In the mining method selection proposed, geometry, grade distribution, and rock mechanics characteristics are ranked according to their acceptability for ten general mining methods.

Jan 1, 1981

By R. R. Harvey

Recent evidence has indicated that the films forming at crude oil-water interfaces are the result of naturally occurring su,rface-active components in the crude oil rather than oxidation products resulting from exposure of the oil to air. Since the presence of such films may strongly affect the displacement of oil from porous media, a knowledge of their composition and formatiorz is of fundamental importance. In this study, without the assumption of a specific mechanism, film formation is considered to involve the transfer of a weakly surface-active species in solution in the oil to the oil-water interface. The film gas law and the Gibbs adsorption equation may thus he applied to the system. The rigid films which form upon decrease of the oil-water interfacial area are considered a result of the very slow rate of re-solution of the film material into the oil. When the shrinkage of the interface is very rapid, the amount of film material leaving the interface is negligible, and the system may then be treated in the same manner as an insoluble film. The combination of these treatments leads to a new concept of the film factor in which it may he quantitatively related to the concentration of film-forming material at the interface. Suitable experiments can then yield the following information. 1. Original concentration of surface-active species in the oil. 2. The change of concentration with time. 3. Concentration of surface active species at the interface. 4. The change of this concentration with time. 5. Interfacial tension of the system at any time and its variation with time. 6. Equilibriunz concentrations for a given set of conditions. 7. Concentration requirements for solid film formation. 8. Molecular diameter of the adrorbed species. These calculations may be made solely from data taken from measurements on the original oil. While using the pendant drop method of interfacial tension measurement of crude oil-water systems, Bartell and Niederhauser&apos; remarked on the thin, transparent, pliable, solid films or membranes which were often formed at the interface between crude oil and water. As a result of their observations of these films, Bar-tell and Niederhauser concluded: 1. The films are formed by material from the oil adsorbed at the interface rather than by a reaction at the interface. 2. The adsorption of the film-forming material is a rather slow process. 3. There are good indications that the process is reversible; i.e., the film constituents may re-dissolve into the oil. 4. The film will spontaneously spread over any new surface formed. 5. The process of adsorption at the interface is accompanied by a decrease in the interfacial tension of the system. Although they had proved that film formation was not due to a reaction at the interface, Bartell and Niederhauser were unable to demonstrate conclusively whether these film-forming constituents were present in the original reservoir oil or whether they were formed as oxidation products after exposure of the oil to air. Recent evidence seems to point to their presence in the original oil. Although the true nature of the film-forming constituents is still uncertain, these conditions would be fulfilled by the assumption of the presence of a slightly soluble surface-active component in the original oil. This component would migrate to a freshly formed interface because it would thereby lower the free energy of the system. The component at the surface constitutes an expanded monolayer which would become compressed, upon retraction of drop, due to the slow rate of re-solution of the film components into the bulk of the oil. It is possible, by involved extraction processes, to separate or at least concentrate this film-forming fraction. However, using the assumptions made above, the technique of monolayer film study may be applied to reveal much of the nature of the adsorbing material without resorting to an actual separation of the film-forming constituents.

By H. C. Chang, N. J. Grant

Creep of very coarse grained, high purity aluminum was studied at 400° to 1100°F with an initial stress range of 50 to 1200 psi. The process of boundary sliding and migration was studied. The driving force of boundary migration under creep conditions is shown to be a combination of strain energy migrationand surface energy. A theory is presented regarding the strainenergyrole the grain boundaries play under creep conditions. From this thetheory intercrystalline failure of commercial alloys can be explained. Also from this theory an optimum grain size should exist for good high temperature properties of high purity materials. IN the investigation of the creep of very coarse grained, high purity aluminum previously reported,&apos; the cyclic nature of the process of grain boundary sliding and migration has been established by microscopic observation of the actual creep process as well as metallographically. It was also substantiated by measuring localized strains across the grain boundaries. In addition, it was shown how extensive subgrain and fold formation may result from boundary sliding and how the slid grain boundary may become sharply wavy as a result of subgrain formation. Since the publication of those results, additional information on the process of grain boundary sliding and migration has been obtained. It is the purpose of this paper to report on the direction and driving force of boundary sliding and migration, the effect of grain size, and the effect of temperature. Finally the effect of decreasing purity (alloys 2s and 3s) on the process of boundary sliding and migration is noted. From these results, an explanation is presented for the intercrystalline failure of alloys under creep conditions. The effects of grain size on the creep properties of impure and high purity materials are discussed in the same light. The experimental procedures and technique have been presented elsewhere.&apos; It is only necessary to state how the 2s and 3s aluminum specimens are prepared. Machined specimens were first annealed at 600°F for about 12 hr. In order to produce a large grain size the specimens were then stretched 3 to 8 pct. The specimens were marked with reference marks, electropolished, annealed at 900°F for 24 hr and at 1150°F for 12 hr, and repolished. The grain size of 2s aluminum produced by this treatment is comparable to that of high purity aluminum in the direction of its width of the specimen, i.e., there are 2 to 4 grains across the width of the specimen. However, the grains of 2s aluminum are 3 to 10 times longer in the direction of the specimen axis. On the other hand, the grain size of coarse grained 3s aluminum ranges from 10 to 20 grains across the width of the specimen due to the treatment just outlined and therefore is much smaller than that of both high purity and 2s aluminum. As in 2s aluminum, the grains are about 3 to 10 times longer in the direction of the specimen axis than in the direction of the width of the specimen. The arrangements of the grains in 2s and 3s aluminum were generally irregular; therefore it was difficult to follow the course of the grain boundaries. Results Method of Presentation: In order to show the arrangement of the grains, especially the direction and positions of the grain boundaries, traces of the grain boundaries on the four surfaces of the specimens were unfolded and mapped out after macroetching. Figs. 1 and 2 show the structural development drawings of specimens P-3 and P-8 respectively. The two flat surfaces (of which the front one faces the microscope) and the round-edged surfaces will be called the "front surface" and "back surface" and

Jan 1, 1954

By R. J. Wagner, Jim Douglas Jr., P. M. Blair

The calculation of the behavior of an oil reservoir during a water flood has long been an important problem to reservoir engineers. Buckley and Leverett derived the differential equation which describes the displacement of oil from a linear porous medium by an immiscible fluid, but this equation could not be solved by the methods of classic mathematics. Consequently, in order to integrate the equation over the length of the reservoir, they neglected the effects of capillary pressure. In the present paper, a numerical method has been developed for determining the behavior of a linear wafer flood with the inclusion of capillary pressure. The differential equation which was derived for the case of incompressible fluids is second order and non-linear. This differential equation was approximated by an implicit form of difference equation which is second order correct in both time and distance. An electronic digital computer was used to perform the numerical solution of the difference equation. INTRODUCTION The problem of calculating the flow and distribution of fluids in an oil reservoir subjected to a water flood has long challenged the reservoir engineer. The ability to solve this problem would provide a valuable tool for the design and study of field waterflooding programs. One of the first contributions in this field was made by Buckley and Leverett,&apos;.&apos; who developed a method of calculating waterflood performance in a linear reservoir. Their technique was limited by the practical necessity of excluding quantitative consideration of capillary pressure. It is the purpose of this paper to describe a method for calculating the behavior of a linear water flood with capillary pressure considered. This method, although limited to the linear case, should serve as a step toward the solution of the two- or three-dimensional waterflooding problem which would better describe actual reservoirs. The treatment of the problem has assumed that the reservoir is linear and homogeneous, that both oil and water -are incompressible, and that gravitational forces may be neglected, and that water is injected into one end of the reservoir and oil and water are produced from the other. Though these assumptions may not be directly applied to natural reservoirs, they are closely approached in many laboratory core tests. The physical problem has been represented mathematically by means of a differential equation and suitable boundary conditions. The inclusion of capillary pressure effects causes the differential equation to be second order and non-linear, and it is not amenable to solution in terms of known functions. Nevertheless, the problem may be solved by numerical methods. The authors have successfully calculated several cases from the start of water injection through to depletion of the hypothetical reservoir using realistic capillary pressure and relative permeability characteristics that were chosen so as to put the method to a test. DIFFERENTIAL EQUATION In two- or three-dimensional problems involving the two-phase flow of incompressible fluids through a porous medium, the equations describing the flow can be combined until only two dependent variables remain, the pressure in one phase, and the saturation of one phase. These variables appear in two simultaneous partial differential equations which are derived on the basis of the conservation of mass and Darcy&apos;s law as applied to two-phase flow. For a linear reservoir, this system of differential equations can be reduced to a single, nonlinear, second order, parabolic differential equation for one dependent variable, the saturation of one phase. For incompressible water and oil, it can be shown that the conservation of mass for the two phases requires that Where q, and q, are the water end oil flow rates per unit cross-section, s, and so their saturations, +, the porosity, and [ and T, dimensional space and time variables. Darcy&apos;s law, as extended to two-phase flow including the effect of capillary pressure,&apos; is as follows:

By V. J. Kniazeff, S. A. Naville

The problem of unsteady-state condensate-gas flow through porous media leads to a set of second-order non-linear partial differential equations. Such a set of equations is numerically solved in the case of radial two-phase flow around a well, taking into consideration both the thermodynamical properties of the fluid and the mechanical properties of the reservoir. The fluid properties, reflecting the PVT relationship of the gaseous and liquid phases, are expressed by using the partial specific masses of the two main separator products in these phases. The flow properties of the reservoir rock are expressed by the generalized Darcy&apos;s law for the liquid phase and by a quadratic relationship between the rate of flow and the pressure gradient for the gaseous phase. The numerical solution of the equations for pressure and saturation US radius and time is worked out through programs written for a computer. The evolution of bottom-hole pressures, well productivities or- deliverabilities and effluent compositions with the depletion of the reservoir is easily derived. The application to the Saharian gas-condensate field Hassi R&apos;Mel led to a better understanding of the drainage mechanism. A zone of fairly high liquid saturation develops around the wells, reducing the effective permeability, and represents a loss of condensible products in addition to the PVT-like retrograde condensation. lnside this zone, near the well, the deviation from Darcy&apos;s law for the flow of the gaseous phase governs the well deliver-ability. A back-pressure test has been computed and correlates with the field results. INTRODUCTION Two-phase flow of volatile hydrocarbons, like condensate gas or light crude oil, may be treated as the flow of a binary mixture by an arbitrary division of the chemical components into two groups. This is translated into two equations of mass continuity, which constitute a set of relation- ships for the pressure and the saturations vs the space coordinates and the time. The equations contain the laws governing the composition and the motion of the phases. The problem so defined is solved with the assumption that the compositions of the phases at any pressure are respectively the same as those observed in a PVT measuring cell under differential liberation. In a first series of computations, it was assumed that the flow obeys the generalized Darcy&apos;s law. A satisfactory representation of the retrograde condensation around the well was thus obtained. In addition, the trend toward decreasing effective permeabilities was obtained, and the computed composition of the effluent checked the laboratory values. However, it has not been possible within this basic assumption to represent the non-linear relationship between the production rate and the bottom-hole pressure drawdown as observed for gas wells in the field. Following the advice of A. Houpeurt it was decided to consider the relative permeability to gas as a function of the velocity of the gas Phase.1 The necessary physical determinations were made by E. COstaséque using a method devised by A. Houpeurt and R. Iffly. As numerical processing of the equations progressed, several difficulties were encountered which were overcome through collaboration with the computer manufacturer. This mathematical model of two-phase flow in porous media had been primarily intended for and extensively applied to the case of the Hassi R&apos;Mel gas-condensate field, operated by SN Repal for SEHR, a joint subsidiary of SN Repal and CFP (A). The programs have also found their applications to forecast the behavior of several fields in the Sahara area containing light and volatile hydrocarbons. BASIC EQUATIONS We will consider a zone in the porous reservoir where the flow properties and the in situ composition of the fluids can be assumed to be uniform. A part of the equations of transient flow can be written with the specific gravity of the fluids being taken equal to the sum of the contributions due to the

Jan 1, 1966