Search Documents

By T. J. Neil, C. B. Manula, O&apos

This paper describes a simulation model developed to analyze truck haulage problems as they are found in the open pit mining and construction industries. The model, using a computer, cycles trucks between their assigned loading and discharge points over a measured haulage route. Required parameters for model operation include the number and type of operating trucks and loading units, maximum truck acceleration and velocity, a realistic vehicular deceleration rate, and the ratio of ore to total material loaded at each origin. Equipment performance characteristics, system profile characteristics, and service time distributions are also included. Records are kept of all waiting times at the loading and discharge points and of any truck interference on the haulroads (moving queues), and a current journal of ore and waste production is maintained. In the pursuit of reduced truck haulage costs, management is confronted with three particular problems in truck and shovel mines. 1) What is the best assignment of trucks to shovels in the system? 2) If new equipment is needed, which new shovels and trucks will operate most effectively in the system? 3) Are any physical changes in the haulage layout warranted to increase productivity? For small open-pit mines the answers to these questions may be obvious. However, for larger operations complex problems are encountered with production and system configuration relationships usually being unknown. The mining simulator discussed here was developed to aid management in evaluating alternative pit-haulage schemes to reduce the risk involved in the selection and assignment of equipment in open-pit truck and shovel mining systems. There are two fundamental differences between this model and other truck haulage simulators: one is a difference of degree, the other a difference of kind. This model can simulate more complex systems than previous models by allowing transportation of material from multiple mine faces to multiple mine destinations and by being readily adaptable to a wide variety of mine layouts. The second difference that characterizes the model is the standard simulation of truck movements. This technique enables each vehicle to perform according to its mechanical capabilities and the physical profile of the haul road. Among the similarities between other truck haulage simulators and this model are the recognition of the probabilistic nature of certain service times in the system and the use of a digital computer to operate the model. The computer inherits the staggering computational and bookkeeping work loads and provides a means for compressing the time scale by allowing several days' operations to be simulated in a matter of minutes. DISTRIBUTION SYSTEM Basic aspects of the system analyzed here may be brought into focus by referring to Fig. 1, a schematic of the typical open pit truck-shovel operation programmed into the simulation model. A single stage materials handling scheme is employed to transport both ore and waste. Loaded trucks are directed to a crusher or stockpile if ore is being mined, while overburden is conveyed to a waste dump. Each shovel has a given number of assigned trucks and has two associated haulage routes, one for waste and one for ore. Up to 5 shovels and as many as 10 trucks per shovel (or a total of 50 trucks) can be operated in the model. Essentially the model, using a computer, cycles trucks between their assigned shovels and the proper discharge points over measured haulage routes. Required input parameters for model operation include the number of shifts to be simulated, operating time per shift, number and type of operating trucks and shovels, maximum truck acceleration and velocity, a realistic vehicular deceleration rate, and the ratio of ore to total material loaded at each shovel. Equipment performance characteristics, system profile characteristics, and service time distributions are

Jan 1, 1968

By J. H. Scaff, H. C. Theurer, E. E. Schumacher

The microwave region of the radio spectrum was effectively utilized in radar designs during the recent war and it has become of increasing interest in the field of communications. Work in this field has led to an important use for silicon—that of the point contact rectifiers†1 in the frequency converter of microwave (radar or radio) receivers—and has stimulated considerable interest in the electrical properties and preparation of silicon and its alloys. Silicon is an electronic semiconductor. Its conductivity at room temperature results principally from the presence of certain impurities. While for metals an increase in impurity content increases the resistivity, for semiconductors such as silicon the opposite occurs and, in general, the addition of impurities lowers the resistivity. Silicon materials may be classified into one of two groups depending upon the manner in which the impurities contribute to electrical properties. These have been termed for convenience p-type or n-type. P-type silicon develops a very large positive thermal emf against metals, has a Hall coefficient of positive sign and a low resistance direction in point contact rectification with the silicon positive with respect to the point. N-type silicon, on the other hand, develops a negative thermal emf against metals, has a Hall coefficient of negative sign and a low resistance direction in rectification with the silicon negative with respect to the point. Impurities which produce n-type silicon are called donors inasmuch as these elements contribute to electrical conductivity by donating electrons to an unfilled energy band in the silicon. On the other hand, elements which produce p-type silicon are known as acceptors as these impurities contribute to the electrical conductivity by accepting electrons from a filled energy band permitting what is known as conductivity by "holes" in which the sign of the carriers appears to be positive. A general treatment of the mechanism of conduction in semiconductors from the viewpoint of modern band theory has been given recently by Pearson,2 by Becker, Green and Pearson,3 and by Torrey and Whitmer.4 In this investigation boron and aluminum have been found to be acceptors, and phosphorus, arsenic, and antimony to be donors in silicon. Data on the effect of boron and phosphorus on the electrical properties, when present singly and in combination, have been acquired. These data are discussed in the present paper. Raw Materials Used and Methods for Adding Seeond Constituents Silicon from two sources was used in this work and these will be referred to as A and B respectively. Silicon A is a material of high purity supplied by the Electro Metallurgical Co. It is prepared by chemical purification of "commercial" silicon obtained from the arc furnace reduction of SiO2. It contains 99.8 pct silicon, minimum, with small amounts of calcium, iron, aluminum, boron, and phosphorus as the principal impurities. This material was extensively employed in these studies as well as in the commercial preparation of rectifier materials. Silicon B is a material of high purity from the du-Pont Co. It is prepared by a pyrolytic reduction of Sicl4 and is free of analytically detectable amounts of boron and phosphorus. Its use permitted study of the effect of boron and phosphorus individually on the properties. This material contains, however, spec-troscopic traees of a number of metals, some of which affect electrical and rectification properties. Typical analyses for the two grades of silicon are given in Table 1. To make controlled additions of boron to the charge it was necessary to employ a low boron content master alloy since the quantity of boron to be added was usually only a few thousandths of one percent. The alloy containing nominally one percent boron was prepared by melting chemically pure boron with silicon A using the melting techniques for 320 g silicon

Jan 1, 1950

By P. Chiotti, J. T. Mason

Thermal, X-ray, metallographic, and resistivity data were obtained to establish the phase diagram for the Yb-Zn system. Six intermetallic compounds and five eutectics have been identified. The compounds position at 695°C. No evidence for the compound YbZn,, was observed. Whether this compound is stabilized or inhibited by small amounts of impurities is uncertain. SEVERAL aspects of the Yb-Zn system have been reported in the literature. Much of this work deals with the crystallographic identification of some of the compounds. The YbZn compound has been reported by Iandelli and palenzonal and the YbZn, structure and thermal data for the composition range between YbZn and Yb3Znll have been reported by Michel et al. A compound of uncertain stoichiometry in the range 60 to 68 wt pct Zn with a cubic structure has been reported by Kuz'ma et aL4 The YbzZn17, YbZn,,, and YbZn13 compounds have been reported by both Kuz'ma et a1 and Iandelli and Palenzona.' YbnZn17 has also been identified by Laube. This investigation was undertaken to establish the phase diagram of this system. Ytterbium in many respects behaves as an alkaline earth rather than a rare earth metal6 and it is of interest to note the resemblance between the Yb-Zn system, the Ca-Zn,' and typical rare earth-zinc systems.g'10 MATERIALS AND EXPERIMENTAL PROCEDURES The materials used in this investigation were Ames Laboratory ytterbium and Bunker-Hill slab zinc, 99.99 pct pure. Analyses of the ytterbium showed 38 ppm Ta, 82 ppm 0, and 15 ppm N with Ca, Cu, C, Mg, Al, Si, Cr, Fe, Ni, Sc, Y, Er, Tm, and Lu all less than 50 ppm. The ytterbium was dehydrogenated by heating in tantalum under vacuum for 8 hr at 900°C. The metal was kept under argon except when being cut with shears and when being weighed. Tantalum was found to be a suitable container for alloys in the composition range between ytterbium and Yb2Zn17. With higher zinc concentrations the tantalum was penetrated by the alloy and failed after several days at temperatures of 600" to 700°C. Approximately 20 g of zinc and ytterbium metal were sealed under argon in *-in.-diam tantalum cruci- bles equipped with thermocouple wells. The crucibles were in turn sealed in stainless-steel jackets. The alloys were equilibrated in an oscillating furnace and given various heat treatments in a soaking furnace. The steel jackets were removed and differential thermal analysis data were obtained under an atmosphere of argon. The thermal analysis apparatus was calibrated against the melting points of pure zinc, aluminum, silver, and copper and found to give thermal arrests within 1°C of the accepted melting temperatures for these elements. The thermal analysis apparatus has been described in an earlier paper.g Alloys were further heat-treated as necessary and subjected to metallographic examination and X-ray diffraction analysis. Most alloys were etched with Nital (0.5 to 10 pct nitric acid in alcohol). EXPERIMENTAL RESULTS The results of the thermal analysis are summarized by the points on the phase diagram, Fig. 1. Six compounds and five eutectics were observed. The eutectic temperature and compositions in weight percent zinc are 491°C and 12 pct, 645°C and 32 pct, 687°C and 54 pct, 742°C and 67 pct, and 750°C and 77 pct. The compounds YbZn, YbZn2, Yb2Zng, Yb2Zn17, and YbZnll melt congruently at 650°, 751°, 75Z°, 754", and 755"C, respectively. The compound Yb3Znll undergoes a peritectic decomposition at 695°C. The transformation and melting temperatures for the ytterbium were observed to be 797" and 817°C. re-

Jan 1, 1969

By J. N. Ong

The oxidation of tantalum is assumed to occur by four simultaneous first-order chain reactions; solution of oxygen in the metal, nucleation and growth of a suboxide phase at the metal surface and two phase boundary processes which give rise to two different modifications of Ta2O5. A rate equation is then formulated which accounts for all principal details of oxidation behavior, unifies and correlatss the results of the principal investigations, and expresses the rate as a function of pressure, temperature, and time. For purposes of practical application, the expressions: may be used below and above 700T respectively to determine the regression rate of tantalum in cm per hr, where T is to be expressed in "K and P in atmospheres pressure of oxygm. The equations will reproduce within a factor of two all reported experimental rates in the temperature range 475" to 1400°C and pressure range 0.000026 to 40.8 atm 0. 1 HE oxidation kinetics of tantalum have been investigated thoroughly in the temperature range 50" to 1400 and pressure range 1.32 x l0-' to 40.8 atm pressure of 0,and as a result, the principal features of the oxidation of tantalum are known and have been described in considerable detail. The exposition of the characteristics of the behavior of tantalum in oxygen has been largely descriptive, with considerable attention being devoted to elucidating the structure, composition, and morphology of the resulting oxide products. The weight change vs time results have usually been reported graphically. With the exception of the linear rate below 700 C and the process of diffusion of interstitial oxygen into the metal,6'9 no effort has been made to characterize the reaction rates in functional form. It is the purpose of this paper to demonstrate that when the fundamental principles of chemical kinetics are rigorously followed, the results of the principal investigations will be correlated and unified and the oxidation rale may be expressed as an analytic function of pressure, temperature, and time. In addition to accounting for the temperature and pressure dependency of the oxidation rate, the following minor features will be accounted for as a result of the basic analysis: the initial period of nonline-arity in rate at low temperatures, which decreases with increasing temperature;3'5'7 the rate vs temperature cusps obtained at constant pressure"5 the concentration gradient of oxygen interstitially dissolved in the metal.5'7 The analysis will also be consistent with the overall chemistry of reaction and with all visual, microscopic, and X-ray observations. DEVELOPMENT OF RATE EQUATION As a fundamental basis for the treatment of metal-oxygen reaction kinetics, it will be regarded that the rate of reaction is determined by the concentration and not the activity or chemical potential. It is further assumed that all metal-oxygen reactions are first order complex chain reactions" whose rate is dependent on the concentration of an intermediate species which is a complex function of the oxygen concentration, or oxygen pressure, only. The basic rate equation is therefore: is the weight of oxygen consumed by the reaction, A is that portion of the area of the metal over which the reaction is occurring (cm2), t the time (sec), kf the rate constant for forward reaction, C the intermediate species concentration and f '([Oz]) and f (Po,) are complex functions of oxygen concentration and Oxygen pressure respectively. Reverse reactions may usually be neglected in metal-oxygen systems since a large negative free energy generally accompanies these reactions. Eq. [I] takes different forms when the rate of reaction is controlled by diffusion into or through a phase or at a phase boundary. For diffusion into the metal phase, the rate is given by1' Di is the diffusion coefficient of species diffusing into the metal phase (cm2/sec), and Q is the corre-

Jan 1, 1962

By Gerald Sherman

MINING asbestos ore by caving at Thetford Mines, Quebec, has been described in the Transactions of the Canadian Institute of Mining and Metallurgy in papers presented by the staffs of The Asbestos Corp. Ltd., in 1934, and Johnson's Co. Ltd., in 1941. The production by caving prior to 1949, in round figures, has been 14,500,000 tons of ore by The Asbestos Corp. and 5,500,000 by Johnson's Co. They are the only caving operations in asbestos at the present time. The Jeffrey mine of the Canadian Johns-Manville Co. at Asbestos, 40 miles southwest of Thetford, is well along in preparations for caving, and The Bell Asbestos Mines Ltd. at Thetford Mines, will follow their examples later. GERALD SHERMAN, Member AIME, is Consulting Engineer, Johnson's Co. and Canadian Johns-Manville Co., New York, N. Y. AIME San Francisco Meeting, Feb. 1949, and New York Meeting, Feb. 1950. TP 2828 A. Discussion (2 copies) may be sent to Transactions AIME before May 31, 1950. Manuscript received Nov. 2, 1949. Nearly all of the asbestos produced by the western hemisphere is mined in those districts of Canada. Since it has been the subject of a number of technical papers, this contribution is offered as a review of progress. Some introductory information seems appropriate, for few members of the American Institute of Mining and Metallurgical Engineers are familiar with the mineral, or how asbestos deposits differ from others in caving characteristics. Asbestos is found in peridotite, a basic intrusive rock irregularly altered to serpentine and erratically fractured, well enough to cave in bulk, but leaving many blocks of solid ore of all sizes scattered through the deposit. Chrysotile is the variety of asbestos produced in Canada. Its fibers are the most valuable commercially, having a good proportion of long fibers, strong, silky, flexible, and suitable for spinning. The fibers occur in all lengths, from the microscopic to 3 in., or more, the shorter are more abundant. Their values depend on length combined with other qualities noted above. They occur as veinlets filling fissures in the peridotite, with crystals compactly arranged in parallel across the cracks. They are sometimes lightly frozen to the walls but the partings are distinct, often marked by particles of magnetite adhering to the fiber ends. Sometimes they break away from the rock in slabs an inch or more in thickness. In wide fissures the fibers are not always continuous, the "broken fibers" being separated by lines of parting with occasional scattered grains of magnetite, which is an injurious mineral to be elirninated when possible. In addition to the chrysotile ''Cross fiber," some "slip fiber" is found in fault fractures in limited quantities of less valuable character. With the exception of the Asbestos Corp, ore, in which the content is said to be more uniform, asbestos is eccentrically distributed through the peridotite. It is the general impression that the higher grades are more often found in hard rock. There are rich pockets and leaner areas which range down through various grades to minor barren sections. Wholesale mining is necessary to make sure of recovering all the high grade for there are no reliable indications of its existence. Ore selection only appears in constructing assay-grade limits for caving. The mill yield of the Thetford ore body in 1934 was 6 to 10 pct of the ore milled, made up of all classes, with occasional runs up to 15 pct or so. Quotations for Canadian asbestos* are as follows: Per Ton No. 1 Crude (% in. and longer) $896 @ $960 No. 2 Crude (Y8 to % in.) 350 @ 545 Spinning (less than % in.) 204 @ 378 Other classes, price diminishing with length. Shorts or "refuse" 24 @ 46 • From Eng. and Min. Jnl., Nov. 1948. Asbestos of even superior quality for electric insulation is found in Arizona limestone. It is white, containing very little combined iron, but occurs only in small quantities.

Jan 1, 1951

By E. W. Hough, B. B. Wood, M. J. Rzasa

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By J. Orkiszewski

A method i.s presented which can accurately predict, with a precision of about 10 percent, two-phase pressure drops in flowing and gas-lift production wells over a wide range of well conditions. The method is an extension of the work done by Grifith and Wallis and was found to be superior to five other published methods. The precision of the method was verified when its predicted values were compared against 148 measured pressure drops. The unique features of this method over most others are that liquid holdup is derived from observed physical phenomena, the pressure gradient is related to the geometrical distribution of the liquid and gas phase (flow regimes), and the method provides a good analogy of what happens inside the pipe. It takes less than a second to obtain a prediction on the IBM 7044 computer. INTRODUCTION The problem of accurately predicting pressure drops in flowing or gas-lift wells has given rise to many specialized solutions for limited conditions, but not to any generally accepted one for broad conditions. The reason for these many solutions is that the two-phase flow is complex and difficult to analyze even for the limited conditions studied. Under some conditions, the gas moves at a much higher velocity than the liquid. As a result, the down-hole flowing density of the gas-liquid mixture is greater than the corresponding density, corrected for down-hole temperature and pressure, that would be calculated from the produced gas-liquid ratio. Also, the liquid's velocity along the pipe wall can vary appreciably over a short distance and result in a variable friction loss. Under other conditions, the liquid is almost completely entrained in the gas and has very little effect on the wall friction loss. The difference in velocity and the geometry of the two phases strongly influence pressure drop. These factors provide the basis for categorizing two-phase flow. The generally accepted categories (flow regimes) of two-phase flow are bubble, slug, (slug-annular) transition and annular-mist They are ideally depicted in Fig. 1 and briefly described as follows. RUBBLE FLOW (FIG. 1A) The pipe is almost completely filled with the liquid and the free-gas phase is small. The gas is present as small bubbles, randomly distributed, whose diameters also vary randomly. The bubbles move at different velocities depending upon their respective diameters. The liquid moves up the pipe at a fairly uniform velocity and, except for its density, the gas phase has little effect on the pressure gradient. SLUG FLOW (FIG. 1R) In this regime, the gas phase is more pronounced. Although the liquid phase is still continuous, the gas bubbler coalesce and form stable bubbles of approximately the rame size and shape which are nearly the diameter of the pipe. They are separated by slugs of liquid. The bubble velocity is greater than that of the liquid and can be predicted in relation to the velocity of the liquid slug." There is a film of liquid around the gas bubble. The liquid velocity is not constant—whereas the liquid slug always moves upward (in the direction of bulk flow); the liquid in the film may move upward but possibly at a lower velocity, or it may move downward. These varying liquid velocities will result not only in varying wall friction losses, but also in a "liquid holdup" which will influence flowing density. At higher flow velocities, liquid can even be entrained in the gas bubbles. Both the gas and liquid phases have significant effects on the pressure gradient. TRANSITION FLOW (FIG. 1C) The change from a continuous liquid phase to a continuous gas phase occurs in this region. The liquid slug between the bubbles virtually disappears, and a significant amount of liquid becomes entrained in the gas phase. Although the effects of the liquid are significant, the gas phase is more predominant.

By H. Biloni, B. Chalmers

Metallographic studies of the segregation of solute in and near the surfaces of chill-cast alloys have shown that, under sufficiently severe conditions, nucleation on the cold mold surface is followed by the formation of a disc of solid of the same composition as the liquid. This disc is surrounded by two concentric rings in which the concentration first decreases and then increases. The discs and rings are surrounded by a region in which arms of low solute concentration separate regions of high concentration. These arms start at the periphery of the outer ring and are initially radial; they make a gradual curved transition from the radial direction to a crystallographic direction. A normal dendritic substrcture is established beyond the position at which the arms have become crystallographically aligned. When the cooling is less severe, the disc becomes less irregular in shape and has a concentration minimum at the center. The term "Predendritic" is proposed for these transitions. A qualitative explanation is outlined for the phenomena observed at the surface and in the equiaxed region of the ingots. It is well-established1-' that the normal mode of solidification for alloys under most conditions is dendritic. A characteristic of dendritic solidification is the growth of rods or plates in specific crys-tallographic directions. The remaining liquid solidifies later, usually with a higher concentration of solute than the original liquid. Very little is known, however, about the transition from a crystal nucleus until the dendritic pattern is established. It is clear that there must be a transition because the nucleus is much smaller than the characteristic spacings found in the "steady-state" structure. And it also seems clear that the initial conditions must be important since the undercooling existing during the nucleation of the grains has a remarkable influence upon the grain size4 and the segregation pattern. The present paper reports a study of the micro-segregation in regions where crystals have nucleated, and the conclusions that can be drawn regarding the early stages of growth, in connection with the details of the solidification process. Aluminum of 99.993 purity, A1-Cu alloys containing 0.5 to 20 wt pct Cu, and Pb-Sn alloys with 5 wt pct Sn were cast and examined metallographically at the as-cast surface, close to the surface, and in the interior after casting. The alloys were cast from controlled temperatures against surfaces at different temperatures ranging from room temperature to approximately the liquidus temperature of the alloy. The surfaces were, in various experiments, graphite, stainless steel, copper, and aluminum. The methods used for metallographic examination were as follows. 1) The Nomarski phase contrast interferometer6 was used to obtain information on the topography of the as-cast surface. 2) The distribution of any solute present in A1 99.993 pct and copper in the A1-Cu alloys was studied by the anodic film technique originated by Lacombe and Mouflard7 and previously applied in the study of the solidification of aluminum by Biloni.8 The electropolished surface of the specimen is anodized under controlled conditions so that the thickness of the oxide film and therefore the interference colors under white light are sensitive to the local composition. This technique was applied to the as-cast surfaces and to metallographically prepared surfaces from the interior. 3) If in Al-Cu alloys, the surface is anodically oxidized for several minutes, instead of seconds as in the color-film technique, the copper segregation is outlined by a different level with respect to the matrix. This method is less sensitive but complementary to that for the color films. 4) The distribution of copper revealed and mapped by the metallographic techniques was confirmed by electron probe. 5) In samples with high concentration of copper, evidence for the copper distribution could be seen optically after electropolishing.

Jan 1, 1965

By T. G. R. Rawlins, L. E. Brosselard

Direct platinum carbon replicas have been used to study substrates prior to growth and after initial nu-cleation of the layer. Replicas have been directly stripped and correlations have been made with other replica techniques. Mechanically polished substrates are very flat with occasional areas of "dirt", which is apparently silicon. When such substrates are etched in hydrochloric acid gas, scratches and hexagonal defects (thought to be polishing compound) are revealed, together with areas of apparent redeposition. Chemically polished substrates are much less flat although much more homogeneous in appearance. Heat treat)nents have no effect. Nuclei deposited on these surfaces, using a normal silicon tetrachloride system, have been examined. Nuclei on mechanically polished subsrates occur in random areas only, whereas on chemically polished surfaces they are more uniformly distributed. THE interest in the study of thin films and now their rapidly expanding use in many fields has increased over a period of years. Much of this current interest, particularly in the field of single-crystal or epitaxial thin films, has derived its stimulus not only from possibilities for application but also from fundamental electron-microscope studies which have indicated the manner of growth of such layers from the initial nuclei.1-5 Previous work involving study of growth processes may be divided into three categories. First, vacuum deposition, primarily of metals, on a variety of substrates has been studied by electron microscopy and electron-diffraction techniques. These have shown the three-dimensional aspect of nucleation and the remarkable liquidlike behavior of the growing nuclei which coalesce to form a complete layer. More recently this work has been expanded in two ways: by growth inside the electron microscope this has obvious advantages for kinetic studies; and by growth in ultrahigh vacuum—much of this work has been done on semiconductors.8-11 It is of interest to note that metals, specifically gold deposited on rock salt, will form only polycrystalline layers in ultrahigh vacuum unless the substrate surface is contaminated.12 In general the mode of growth for vacuum-deposited semiconductor films is the same as that for metal films. The second method of growth that has been studied is that of electrodeposition.13,14 The little that has been published suggests that the mode of growth is similar to the vacuum case, except that irregularities may perturb the electric field and disturb the growth process. The third and last primary method of growing epitaxial layers is by use of chemical methods. These techniques have generated much interest and much has been published on them. Considerable electron-microscope work has been done on the defects in the layers; however comparatively little work has been undertaken on the initial nucleation and early stages of growth. Confining ourselves now exclusively to the case of silicon, Charig et a1.15 have considered the growth of silicon on silicon substrates by hydrogen reduction of trichlorosilane, and studied the films by means of replica and transmission techniques. Subsequently they have studied the same system for growth on quartz,16 and with the addition of pyrolysis of silane on quartz and alumina.17 Also Manasevit el al. have studied deposits using silane and silicon tetrachloride on sapphire.18-14 Finally, in this volume studies of deposition on beryllia and sapphire are reported. Although Revesz and Evans20 have examined complete layers, grown by means of the widely used silicon tetrachloride reaction, no detailed examination similar to that published by Charig and Joyce21 for the trichlorosilane system has been given. Therefore this investigation has examined substrate surfaces prior to growth and then nucleation on these surfaces for the silicon tetrachloride system. EXPERIMENTAL In this work the surfaces have been studied by means of standard replica techniques. Because of the very real problems associated with artifacts introduced,22 different replicas and different cleaning

Jan 1, 1967

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'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 M. E. Nicholson

The effect of boron on the austenitic transformation rate of iron is smaller than on low carbon steels. The influence of austenitizing temperature on B-Fe is the reverse of its influence on steels. THERE have been many studies made of the influence of boron on the hardenability of steel, but apparently none has been made of its influence on the y-a transformation of pure iron. There are several facts suggesting that boron should have a pronounced effect on the transformation rate of iron. First, Simcoe et a1.l indicated that boron reduces the nucleation rate of ferrite and bainite. Second, they showed, along with others,', " that the lower the carbon content of steel the greater the influence of boron on hardenability. If this trend extends to pure iron, the influence on rate of ferrite formation should be considerably greater for iron than for medium carbon steels. According to Rahrer and Armstrong," the multiplying factor for boron in pure iron should be about 2.40, compared with 1.50 in a 0.35 pet C steel. According to an extrapolation made by Grange and Mitchell,' the effect should be even greater, and the multiplying factor should be about 3.5 in iron. Many investigations of the boron influence'.7, 2, 5 indicated that boron is effective only when in solid solution. On the other hand, Chandler and Bredig" suggested that boron increases hardenability by combining with nitrogen, thereby preventing formation of AlN, which catalyzes the decomposition of austenite. If their hypothesis is valid, boron should have no effect on an alloy free from nitrogen. Digges et al." showed that so far as low nitrogen high purity Fe-C alloys are concerned this is not the case. No experiments, however, have been made on low nitrogen iron. To determine whether or not boron has an effect on low nitrogen iron and, if an effect exists, to determine whether it is as large as predicted from the investigations of boron treated steels, the following work was performed. Because of the rapidity with which iron transforms from austenite to ferrite at even a moderate degree of supercooling, it was decided to study the influence of boron on the transformation rate of iron by determining the effect of quenching rate on the temperature of the beginning of transformation. The technique, essentially that used by Greninger' and Duwez,H consisted of heating small specimens, 0.1 in. sq and 0.010 in. thick, in an atmosphere of purified helium with a resistance heater coil, and then quenching the sample in a jet of helium gas. At the time the jet was applied, the heating coil was removed from around the sample. The sample was suspended in the furnace on 30 gage chrome1 and alumel wires, one welded to each side of the specimen, which also served as a thermocouple for measuring the specimen temperature. Specimens were prepared from vacuum-melted iron obtained from the National Research Corp. which was rolled to a thickness of 0.060 in. and then decarburized in wet hydrogen for 51/2 hr at 700°C. This treatment also removed any boron present in the iron. The resulting iron contained 0.001 pet C. 0.0002 pet N, and 0.005 pet 0. One half of the piece was then boronized by surrounding it with 325 mesh boron powder and heating it for about 2 hr at 900°C in H,. The boronized iron was then homogenized at 1100°C for 24 hr. The iron boride case formed by boronizing was removed by grinding 0.010 in. from each side of the plate, and the piece was finally rolled to a thickness of 0.010 in. The resulting alloy contained

Jan 1, 1957

By R. B. Adair, J. S. Browning

The laboratory batch and continuous flotation pilot plant tests demonstrated the technical feasibility of recovering high grade mica concentrates from weathered mica pegmatite ores of Alabama and Georgia. The research indicated that combinations of anionic and cationic collectors may be used effectively for flotation of fine size mica from weathered pegmatite ores. In continuous tests, concentrates containing 98.5% mica were obtained from the Georgia pegmatite ore; the Alabama pegmatite ore concentrates contained 98.4% mica. The recoveries were 91 and 89% respectively. INTRODUCTION The principal uses of fine ground mica are as a filler in wallboard joint cement, as a filler and surface coating for roofing, as an ingredient in paints, and in oil well drilling mud. The mineral has other uses in the manufacture of rubber, wallpaper, plastics, welding rods, electric insulation, house insulation, and textiles, and as an annealing agent in metal treatment. In recent years, more than 99% of the domestic mica produced has been scrap and flake mica (mica which does not meet specifications for sheet mica and is used for producing fine ground mica). There has been a continued increase in mica production for several years, the 1963 production of scrap and flake mica totaling 117,251 tons.1 Ground mica is obtained primarily by crushing and milling pegmatites and schists. To a lesser extent, mica is produced as a byproduct of kaolin washing and feldspar and spodumene flotation operations. The processes used in recovering mica by crushing and milling pegmatites are generally simple, consisting of various combinations of alternate roll crushers and trommel screens that separate the mica and gangue at screen sizes coarser than 6-mesh. As these processes are designed to recover only coarse mica, high losses in the plant rejects are common. Large tonnages of tailings from the crushing and screening plants have been accumulated in a number of areas. Methods for treating such products were developed by the Bureau of Mines in 1941,2 but have not been generally applied. More recently, mica flotation research has been completed and published by the Bureau of Mines.3 These methods required complete removal of 150- to 200-mesh materials from the flotation feed with consequent fine mica losses. Later, the Bureau of. Mines investigated methods for-recovering fine size mica from pegmatite ores after desliming sufficiently to remove clay materials, but not so drastically as to remove the fine sands. This report summarizes the results of these studies. The process developed was effective on pegmatite ores from two locations and should be applicable to the commercial treatment of other mica-bearing pegmatite ores and fine rejects that have been accumulated at various mica-milling operations. DESCRIPTION OF ORES The ores used in the investigation were obtained from the Dixie Mines, Inc., Heflin, Alabama, and the Ruberoid Corporation, Hartwell, Georgia. The sample from Alabama contained muscovite and quartz, with a high percentage of clay, and minor amounts of biotite, kaolin, limonite and tourmaline. The mica in the ore was essentially all minus 4-mesh in size and was free of attached mineral grains. The Georgia sample contained muscovite and quartz, with minor amounts of biotite, kaolin, and limonite. The mica in the ore, which was essentially all minus 4-mesh, was liberated. Petrographic analyses of the two samples are given in Table I. THE ANIONIC-CATIONIC MICA FLOTATION METHOD Previous investigators2,3 have reported that complete desliming of mica ores at 150- to 200-mesh was required prior to flotation with cationic collectors to obtain satisfactory selective separation of the mica from the other mineral components. Numerous tests were made at the Tuscaloosa Metallurgy Research Center to determine if some reagent combination could be used to selectively float finer size mica without complete desliming. The investigation led to development of a process using a simple reagent com-

Jan 1, 1967

By F. W. Jessen, H. Kazemi

A washing technique has been developed to form a spherical cavity in massive salt. The technique is, basically, a process of controlling the fluid motion in the cavity, the concentration distribution of brine, the rate of dissolution at the walls of the cavity, and the particular settings of the wash pipe, or pipes (inlet and outlet pipe) employed. For a particular arrangement of tubing and casing, fresh water enters from the casing annulus at a position above the bottom of the tubing, and the dissolved material leaves the cavity through the tubing near the bottom of the projected hole. The initial configuration of the cavity is cylindrical, representing the drilled hole situated at the axis of the projected cavity. As washing progresses the original hole is enlarged radially (any direction perpendicular to the axis of the hole). An inert fluid (generally a hydrocarbon, natural gas, or air) which serves as a blanket is injected through the annular space between the original hole and the casing. The controlled downward motion of this blanket maintains the upper edge of the brine region that is expanding on the surface of the projected sphere. A large spherical glass container and a radial model were used to investigate the flow patterns and the concentration distribution of brine occurring during the progress of solution of the salt. Aluminum flakes and dye were injected and flow patterns were observed and recorded on motion picture film. Based on the results of these model studies and other work carried out in this laboratory, the controlled washing technique was formulated. INTRODUCTION Underground cavities leached in salt beds and salt domes have been used extensively for storing hydrocarbons. Little attention has been focused on the shape of the cavities as long as large volumes were obtained and the cavities were structurally stable. The general process involves circulating fresh water and removing brine. Two distinctly different circulating systems can be employed, each of which yields characteristically different shapes. By far the most generally accepted method of washing storage caverns in salt domes is the direct circulation system. Water is introduced through the innermost pipe of a concentric pipe system, and brine is withdrawn through the annulus. In the reverse-circulation method, water enters through the annulus and is discharged through the smallesc tubing string. The size of the operation, and hence the size of tubing and casing used, depends on such factors as cavity volume, length of hole, availability of water, disposal of brine and power requirements. Underground caverns are from 100 to 1,000 ft or more in length, and when completed vary from 40 to 100 ft in diameter, providing storage of 100,000 bbl or more. The volumes of water required per volume of salt removed vary from 6 to 7 for the reverse circulation method to 10 to 11 for the direct circulation method. Up to 80,000 BWPD are often circulated. In the direct circulation method, because water enters near the bottom of the projected cavity, solution of salt is greatest at this point, particularly in view of the fact that when the washing is begun turbulent flow undoubtedly exists in the immediate vicinity of the discharge of the inlet water. As the diameter of the hole increases, flow rates also are increased, so perhaps the condition of fluid turbulence is maintained for some time. During the later stages of leaching, when capacity of equipment may be the limiting factor in maintaining flow rates, there may still exist some turbulence adjacent to the point of influx; however, for the most part laminar flow conditions prevail. There appears, therefore, rather good evidence to indicate that the state of fluid motion is responsible for the typically shaped cavity which results from such washing — namely a cylindrical form, generally somewhat larger at the bottom than at the top. Depending on the number of so-called stages which might be employed, the walls or surfaces of the salt diminish rather uniformly from bottom to

Jan 1, 1965

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 P. Datta, L. L. Handy

For a wetting phase displacing a nonwetting phase from a porous medium the distribution of the residual fluid may depend on displacement conditions. Although this subject has been debated in the literature, only a few experiments have been cited to support the various conclusions. Experimental results presented in this paper show that fluid distributions are dependent on imbibition pmcedures. Results agree qualitatively with predictions from the pore doublet model. If the rate of water imbibition is restricted, the nonwetting phase is trapped preferentially in .the larger pores as expected. But if the rate of water imbibition is unrestricted, trapping occurs somewhat more in the smaller pores. These conclusions were deduced primarily from relative permeability measurements. INTRODUCTION Relative permeabilities are known to depend on the saturation history of the porous medium. For either continuously increasing or continuously decreasing wetting phase saturations, however, they have normally been assumed to be single-valued functions of saturation. Severat investigators have compared relative permeabilities measured by different methods and have reported acceptable agreement.1-5 Studies of simple models of the displacement process suggest that the distribution of residual fluids can be influenced by the displacement method. If this is so, relative permeabilities also should depend on the method of displacement. These predictions have not been supported by experimental data. The object of this paper is to investigate experimentally some of the predictions of these model studies. The particular question of importance is whether steady-state and unsteady-state methods should be expected to give the same values of relative permeability. Much of the reported data showing agreement between steady-state and unsteady-state methods have been obtained on un consolidated sand. Johnson, Bossler and Naumann, however, compared steady-state and unsteady-state results for Weiler sandstone and found no significant differences.5 Levine made a detailed study of pressure and saturation distributions for a laboratory waterflood in a consolidated system but made no attempt to calculate performance from steady-state relative permeability data Bail and Marsden in a somewhat similar set of experiments did attempt a comparison of observations with predictions but the results were inconclusive.7 Excluding gravity, two forces affect the distribution of wetting and nonwetting phases during an unsteady-state, immiscible displacement: viscous and capillary forces. If relative permeabilities are indeed the same when measured by steady- or unsteady-state methods, the fluid distributions at the same fluid saturations must be similar for both methods. Furthermore, the implication is that the relation between capillary and viscous forces is the same for both methods insofar as the effect on microscopic fluid distributions is concerned. Unsteady-state displacements are normally run at high rates to maximize the ratio of the viscous to capillary forces. The objective is to reduce the effect of capillary forces on macroscopic fluid distributions. For floods in which the phase which wets the porous medium displaces the nonwetting phase, capillary forces compete with viscous forces in determining fluid distributions.8 The residual nonwetting phase is discontinuous. Competitive aspects of viscous and capillary forces and the resulting effect on water-oil distribution in porous media have been illustrated in an elementary way using a pore doublet model. This model and predictions from it have been discussed at considerable length by Rose and Witherspoon, Rose and Cleary and by Moore and Slobod.9-13 Rose and Witherspoon show that so long as water invades the model at a rate equal to or greater than the free imbibition rate, the water will move through the larger capillary of a doublet first and trap oil in the smaller capillary. On the other hand, if a negative pressure is applied to the water

Jan 1, 1967

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. Sandoz, D. A. Meyn

The fracture surfaces of specimens of titanium a1loys which exhibited susceptibility to subcritical crack growth in a wide variety of environments, including aqueous solutions, alcohols, hydrocarbon gases, carbm tetrachloride, and dry air, were examined. The dominant and characteristic fracture mode was cleavage. which became mixed with an increasing porportion of ductile fracture modes (dimples. and so forth) as the applied KI was increased. The cleavage plane in all the alloys was oriented at 15 deg from (0001) of the a phase. No indices could be assigned to the 15 deg cleavage Plane because of uncertainty as to which zone it belonged to. SLOW crack growth and delayed fracture take place in certain titanium alloys at stress intensity (KI) lower than the critical values for complete fracture (KIc) found in short-time tests for plane strain fracture toughness.' Such crack growth may take place in very dry air on sustained loading if the alloy contains a sufficient amount of hydrogen and if the stress intensity factor is higher than a threshold value designated KIH.2 Similarly, a number of environments such as water, alcohols, and hydrocarbons may cause slow crack growth provided the stress intensity factor exceeds the threshold value KISCC. In this paper it is intended to examine some fracto-graphic and crystallographic details of this subcritical crack propagation process as related to alloy composition, environment, and stress intensity. It has been previously reported that subcritical crack growth takes place in titanium alloys in aqueous environments by cleavage.2"6 The observations to be described here indicate that in the alloys investigated the mode of subcritical crack propagation is more dependent on the stress intensity during cracking than on the environment. The fracture mode is illustrated with photographs of the fracture surfaces at a wide range of magnifications. The crystallographic character of the cleavage mode is also presented and the results discussed. EXPERIMENTAL APPROACH The studies were made on four titanium alloys— Ti-8A1-1Mo-lV, Ti-5A1-2.5Sn, Ti-7Al-1Mo-lV, and Ti-0.350. Of these four alloys the first two were mill annealed and fine grained. The latter two were noncommercial and coarse grained, having been prepared for other experimental purposes. The crystallographic studies required a large-grain size. Consequently the fine-grained alloys were an- nealed in a vacuum at 1950°F and furnace cooled to produce samples for these studies. This treatment did not alter the microstructure significantly except for grain size. All the aluminum-containing alloys showed primary a with a network of ß-phase particles within the a grains; Figs. 1 and 2 show this microstructure for the specific case of fine-grained Ti-8A1-1Mo-1V. The cracks in the figures will be discussed later. Fractographic studies were made on all the alloys, but the electron fractographs displayed herein are from the Ti-8A1-1Mo-1V alloy in the fine-grained condition. Emphasis was placed on this commercial alloy because considerable information was available on the stress-corrosion cracking characteristics in a variety of environments. Sufficient supplementary studies of the alloy after the anneal to produce large grains were conducted to establish that the stress-corrosion cracking susceptibility and fractographic characteristics were not changed significantly. Thus it could be inferred that the crystallographic data obtained exclusively on large-grained specimens would apply also to fine-grained specimens. Fracture surfaces for study were obtained using precracked cantilever beam specimens loaded as described by Brown. The precracked section of a specimen was immersed in the environment of interest. Values of initial stress intensity factor Kli were chosen slightly higher than KISCC in order to obtain a fairly extensive region of slow crack growth. Fracture surfaces representative of a variety of stress intensity levels could then be examined because in the cantilever beam test KI increases as the crack propagates.

Jan 1, 1970

By K. H. Ribe

Water often enters an oil reservoir during completion or workover operations on a well and forms a partial "water block" to oil production. A mathematical study of radial two-phase flow, neglecting capillary effects, has been employed to study the formation of such a water block and subsequent re-moval by production from the well. The effects of reduced oil permeability about the well on the well productivity were studied. The fluid saturation distributions about the well during formation and removal of the water block have also been computed. Several relative permeability relations and viscosity ratios were employed. If water has invaded the formation, its influence through relative permeability effects alone can cause the following. 1. Oil productivity will be depressed for extended periods after production is resumed and will build up only gradually as the water is removed. 2. Oil injected for treatment of water blocking will delay rather than promote restoration of full well productivity by enlarging the region invaded by water. Thus, unless the specific action of chemicals contained in the oil is needed, oil injection appears undesirable. INTRODUCTION During oilwell workover operations, water may enter the oil-bearing formation from the wellbore. When production is resumed, oil must flow through the region invaded by this water. The presence of this region can cause both well productivity and oil production rate to be low and oil to be produced with high water-oil ratio for some time after production is initiated. This situation is sometimes described as a water block. The introduction of water into the formation may result in other actions which also lead to reduction in well productivity and which are also usually included in the connotation of the broad term, water block. Often considered, for example, are the possibilities of clay swelling by contact with fresh water and the formation of emulsions with the formation oil. If it is suspected that such specific actions have taken place, remedial treatments are undertaken which usually involve the injection of chemicals in oil. Since the introduction of water, even in the absence of specific interactions with the formation or oil, will cause a temporary water block (which might be misinterpreted as evidence of a more severe situation), it is of importance to evaluate the magnitude and duration of this blocking which results purely from the reduction in relative permeability to oil in the vicinity of the wellbore. It is also of interest to evaluate the effect of oil injection on the productivity of a well blocked by water in this manner. Inasmuch as this unfavorable condition may persist for some time, it may lead to premature condemnation of a workover or premature abandonment of a potentially productive pay zone. A quantitative evaluation of the influence of water entry on the oil productivity through changes in relative permeability was made by solving a radial form of the Buckley-Leverett equation. The distribution of water saturation around the wellbore during the entry of water was calculated and was followed by a similar calculation of the saturation distribution during the period of resuming production. At any stage in the removal of the invading water, knowledge of the distribution of its saturation permitted calculating the attendant loss in oil productivity. The influence of the shape of the relative permeability relationships was also evaluated by carrying out the calculations for two hypothetical cases. Further, the effect of the oil-water viscosity ratio was examined by repeating the calculations, for several ratios of unity and greater, with the same relative permeabilities. Fi-nally, results are presented to show how the length of time a well must be swabbed to resume production depends on the length of time it has been subjected to invasion by water. STATEMENT OF THEORY Differential Equations Water is assumed to enter a producing formation which is initially at the connate-water saturation and contains no gas. The water and oil are treated as in-

By J. W. Pugh, J. D. Nisbet

THIS study of the ternary has been made as one phase of a metallurgical investigation which began nearly four years ago in the General Electric Company's Research Laboratory in Schenectady, N. Y. The objective of this program is the discovery of new metallurgical information which will lead to the development of better high-temperature materials. Combinations of the four pure base elements—iron, chromium, nickel, and cobalt—are being studied at the present time. It is essential in an investigation such as this to know as much as possible about the constitutional diagrams involved. The study of the iron-chromium-nickel system has been made to this end. Preliminary Explanation of Experimental Procedure: In all, fifty-five alloys at steps of 10 at. pct were made for the investigation of this ternary. Several ternary alloys in the chromium rich corner had to be omitted because their extreme brittleness made testing impractical. All alloys were vacuum melted with hydrogen reduction and centrifugally cast. The apparatus and technique of this process has been described in detail by Nisbet.1,2 The purity of the alloys prepared in this way is considered to be quite good. Impurities are listed as follows: Pct Carbon.............................0.02 Oxygen................................0.02 Nitrogen.........................0.005 Magnesium..................0.03-0.05 Sulphur...........................trace Hydrogen....................trace Phosphorous..................trace Silicon........................ trace After casting, all samples were given a homo-genization treatment which consisted of holding them for 15 hr at 1150°C (2100°F) and water quenching. Testing was begun with the samples in this condition. The tests employed in the study of this system were (1) dilatometer, (2) hardness versus aging temperature, (3) tensile strength and elongation versus temperature, (4) microstructure analysis, and (5) electrical resistance versus temperature. Dilatometer tests were made at the rather rapid heating rate of 1093°C (2000°F) in one hour on a Bristol-Rockwell type instrument. Hardness data were taken on specimens cooled from 204", 427", 649°, 760°, 871°, 982°, and 1093°C under conditions which are thought sufficient to bring the specimens to equilibrium in all but the cases of the very sluggish transformations. Electrical resistance data for several specimens were taken for both heating and cooling conditions in a vacuum furnace especially designed for this work by D. W. Bainbridge, formerly of this laboratory. Provisions were made for heating and cooling standard specimens at a constant rate, while autographic records of resistance and temperature were made simultaneously. Micrographs of all the alloys were made in the quenched condition. Some question may exist as to how such physical values as hardness, tensile strength, and elongation were interpreted to indicate a change in phase. Original data were recorded on physical property versus temperature graphs, each of which was made from the data of a single alloy (fig. la). From these, another series of graphs were plotted with the physical property as a function of composition for constant temperatures (fig. 1b). Sharp deviations in the slope of these composition versus hardness curves often indicate a change in phase. The example of fig. 1 will serve to describe this experimental technique. The plotting of a single point from graph "a" to graph "b," and finally to the phase diagram "c," is illustrated. The location of the point on graph "b" is indicated as point (1). If all the alloys of the A-B binary system were solid solutions at temperature X, the hardness in this curve would be expected to rise at a fairly even rate as A is diluted with B to a maximum at some intermediate value between A and B, as

Jan 1, 1951

By T. D. Mueller

Many investigators have used the response of the "dimensionless aquifer" to a unit pressure drop or a unit fluid-withdrawal volume to calculate the performance of an aquifer in supplying water influx to an oil reservoir. In the past, these response functions have been calculated with the aid of the Laplace transform. With the advent of ultra highspeed digital computers, it becomes practical to solve for the response functions with finite-difference techniques. The computer method also permits extension of the dimensionless-aquifer concept to include the nonhomogeneous aquifer wherein the permeability and other properties vary as a function of the space co-ordinates. This paper gives results of calculating the response functions for a series of nonhomogeneous aquifers. Response functions are presented for both linear and radial aquifers whose thickness, permeability-viscosity ratio and porosity-compressibility vary. These functions are new and should prove useful to the petroleum engineer in analyzing the behavior of nonbomogeneous aquifers. Results are presented in the form of charts that can be easily used by the field engineer. INTRODUCTION Aquifers which surround many oil and gas reservoirs have the ability to supply water influx to such reservoirs as oil and gas are withdrawn. This water influx, called natural water drive, provides one of the most effective driving mechanisms for the production of oil and gas. In producing a reservoir, therefore, it behooves one to make the maximum use of natural water drive. To achieve the maximum use, the reservoir engineer must be able to predict the performance of an aquifer under a variety of production schemes that may be proposed for the reservoir. Unfortunately, the physical properties which dictate aquifer behavior often are known only within limits. Seldom do wells penetrate the porous strata of the aquifer. Even when they do, quantitative information regarding porosity, permeability and water compressibility is seldom available. It is known, however, that the water efflux from most aquifer systems is governed by a single, relatively simple, linear, partial differential equation. Also, the general physical location of the aquifer boundaries often are known. A technique originall proposed by Hurst1 and van Everdingen and Hurst2 has been found useful in analyzing reservoirs in this situation. The idea was later expanded by van Everdingen, Timmerman and McMahan3 to include the mathematical technique of least-squares fitting. This latter method will be referred to as the VTM method. The basic assumptions of the VTM method include the following. 1. The location of the physical boundaries of the aquifer are known. 2. The flow conditions at these physical bound-daries are also known. 3. The aquifer is homogeneous; e.g., thickness, permeability, water compressibility and porosity are constant throughout. In the VTM method, a material balance is made on the fluids entering and leaving the reservoir. In the balance equations, the water-influx term is represented as the product of the water influx from an arbitrarily-selected, dimensionless aquifer system times an unknown conversion number. This balance can be formed as many times as there are data points in the history of the reservoir. Each time, the conversion number can be evaluated. If the reservoir engineer has picked the correct dimensionless aquifer to represent the real aquifer, the conversion number remains constant for all balances that have been made over the history period. If such a situation occurs, the reservoir engineer can then predict the performance of the reservoir for any type of production scheme by using the function associated with that particular dimensionless system and the derived conversion number. These functions will be referred to as the response function of the aquifer. If the aquifer is nonhomogeneous (e.g., if the porosity, permeability, thickness, porosity, or