Search Documents

By W. T. Weller

In Part I, a study of pressure build-up curves calculated for conditions under which both oil and gas flow led to the conclusion that the presence of a dispersed free gas phase in an oil reservoir must be taken into consideration to estimate accurately average reservoir pressure and permeability from build-up curves. Familiar methods based on the assumption of no free gas can be extended to the two-phase case by using total compressibility and mobility in place of single-fluid compressibility and mobility. These methods give correct values for average pressure and permeability when gas saturation is small. Errors become larger as the gas saturation increases. However, for the use to which the results will be put, the methods are satisfactory for reservoir engineering purposes. An improved method of calculating the performance of depletion-type reservoirs is presented in Part 2. Because the mathematical relationships describing simultaneous flow of oil and gas are quite involved. simplifying assumptions are made to provide means of obtaining approximate solutions of reasonable accuracy. One such approximate method now in use is the constant-GOR solution. It involves the assumption that at any instant, the ratio of total gas flow rate (both free and disolved) to oil flow rate is the same at all points in the reservoir. The approach is not applicable unless the free gar saturation in the reservoir is everywhere greater than the critical gas saturation. This paper presents a modification which, by avoiding the constant-GOR assumption, makes the method applicable to all reservoir conditions, and so far appears to be more accurate than the constant-GOR solution and to he comparable in required compuctation time. PART I— BUILD-UP CURVE ANALYSIS INTRODUCTION Pressure build-up characteristics of shut-in wells have been used for many years by engineers to evaluate average reservoir pressure, effective permeability thickness of the pay section, effectiveness of well completion (skin effect) and reservoir size. A number of methods of analysis have appeared from time to time."" Without exception, these methods are based on the assumption that the reservoir contains but one fluid of constant small compressibility and constant mobility. It has been suggesteda" hat these single-fluid methods may be applied to data from reservoirs containing both oil and gas by substitution of some effective total properties of the multiphase system in place of the corresponding single-phase properties. The present investigation was undertaken primarily to evaluate this approach. METHOD A number of theoretical build-up curves were calculated for conditions of two-phase flow, under the assumption of certain reservoir and fluid properties, and were analyzed by single-fluid methods with appropriate total compressibility and total mobility values for the corresponding single-fluid properties. Results of the analyses were compared with the assumed conditions. The theoretical build-up curves were completed by procedures similar to those of West, Garvin and Sheldon." Since these calculations require a considerable amount of computer time, an attempt was made to derive an approximate calculation method. The attempt was unsuccessful for calculating build-up curves, but the effort did result in a new approximate method of calculating the performance of solution gas drive reservoirs, which appears to be an improvement over the constant-GOR method" used previously (see Part 2). The West, Garvin and Sheldon calculations involved the following assumptions: (1) the reservoir is circular and completely bounded, with a completely penetrating well at its center; (2) the porous medium is uniform and iso-tropic, with a constant water saturation at all points; (3) gravity effects can be neglected; (4) compressibility of rock and water can be neglected; (5) the composition and equilibrium are constant for oil and gas; (6) the same pressure exists in both the oil phase and the gas phase; and (7) no afterproduction occurs, i.e., the well is shut in at the sand face for build-up. These assumptions make it possible to describe two-phase flow of oil and gas by the partial differential equations:

Jan 1, 1967

By Fred C. Bond, Jen-Tung Wang

Comminution energy is principally energy of deformation before breakage, which appears as heat An empirical equation is presented which covers the entire comminution range. The new strain-energy theory considers comminution from the known principles of mechanics and the reduction ratio. Energy requirements according to the different theories the romnnred. THE present status of the theory of comminution is extremely unsatisfactory. The amount of various ores and rock products crushed and ground annually approaches the staggering sum of perhaps one half billion tons. Yet the basic concepts underlying an operation of such magnitude are unknown, and actual knowledge of what takes place in comminution is almost entirely empirical. The responsibility for this situation cannot be charged to apathy or lack of appreciation of the importance of the problem. Many attempts have been made to evolve a workable hypothesis, but none has been completely successful. General Considerations In order to break rock, it must be subjected to a stress which strains the rock beyond its critical breaking point. The stress imparts energy to the rock, most of which is released in the form of heat when the stress and resulting strain are removed. If the strain induced exceeds the critical strain, the energy is released by breaking; if the induced strain does not reach the critical point, the energy is released with removal of the stress as energy of resilience. Rock is a brittle material and it is assumed that it breaks approximately at its yield point, so that no permanent deformation results from strains below the critical. When rock is broken the total energy input is accounted for by the heat liberated, and by the surface energy of the new surface produced. A small and probably negligible amount is released as noise. Rock is commonly broken under compression and the applied stress is ordinarily compressive. However, for breakage to occur, it is only necessary that the induced strain exceed the critical value, and this strain may represent the resultant of compressive, shearing, and tensile forces. The stress-strain diagram of rocks under compression, and the modulus of elasticity as determined therefrom, should be fundamental considerations in the development of any theories of rock breakage. Crushing in jaw or gyratory crushers results primarily from a squeezing action, and reduction in a hammer mill is the result of impact. The action in rod and ball mills is a combination of impacting, squeezing, and wearing away by attrition, or rubbing. In ordinary crushing and grinding, the forces are applied at protruding points and are not distributed evenly throughout the rock. This is one of the features of comminution which theoretical considerations have not covered adequately. Another feature is the effect of impact velocities, which may greatly influence the total energy required for breakage. Impact velocities may vary from 1 or 2 fps in crushers to 10 or 20 fps for ball mills, and as high as 100 or more for hammer mills. However, according to the findings of geophysicists, the velocity of sound and of compression waves in stone is many times greater; the primary longitudinal compression waves travel at perhaps 15,000 fps, and the secondary transverse waves, which may reach the surface and cause cracks to form thereon, travel at perhaps 5000 fps. When a crack tip forms, the total surrounding stresses are concentrated in this tip, which rapidly extends throughout the rock particle.1 It would seem that the energy required to deform the rock beyond its critical strain, resulting in the formation of a crack tip, represents practically all of the energy required. Most of the primary crack tips presumably form on the surface. A simple method of designating feed and product size is necessary for the evaluation of crushing and grinding results. Taggart's suggestion2 that the 80 pet passing size be used has been found very practical. In this paper, the square opening screen size which 80 pct of the feed passes is designated as the feed size F, and the size which 80 pct of the product passes, or the product size, is designated as P. The reduction ratio at the 80 pet passing size, or F/P is designated as n. The feed and prod-

Jan 1, 1951

By Fred C. Bond, Jen-Tung Wang

Comminution energy is principally energy of deformation before breakage, which appears as heat An empirical equation is presented which covers the entire comminution range. The new strain-energy theory considers comminution from the known principles of mechanics and the reduction ratio. Energy requirements according to the different theories the romnnred. THE present status of the theory of comminution is extremely unsatisfactory. The amount of various ores and rock products crushed and ground annually approaches the staggering sum of perhaps one half billion tons. Yet the basic concepts underlying an operation of such magnitude are unknown, and actual knowledge of what takes place in comminution is almost entirely empirical. The responsibility for this situation cannot be charged to apathy or lack of appreciation of the importance of the problem. Many attempts have been made to evolve a workable hypothesis, but none has been completely successful. General Considerations In order to break rock, it must be subjected to a stress which strains the rock beyond its critical breaking point. The stress imparts energy to the rock, most of which is released in the form of heat when the stress and resulting strain are removed. If the strain induced exceeds the critical strain, the energy is released by breaking; if the induced strain does not reach the critical point, the energy is released with removal of the stress as energy of resilience. Rock is a brittle material and it is assumed that it breaks approximately at its yield point, so that no permanent deformation results from strains below the critical. When rock is broken the total energy input is accounted for by the heat liberated, and by the surface energy of the new surface produced. A small and probably negligible amount is released as noise. Rock is commonly broken under compression and the applied stress is ordinarily compressive. However, for breakage to occur, it is only necessary that the induced strain exceed the critical value, and this strain may represent the resultant of compressive, shearing, and tensile forces. The stress-strain diagram of rocks under compression, and the modulus of elasticity as determined therefrom, should be fundamental considerations in the development of any theories of rock breakage. Crushing in jaw or gyratory crushers results primarily from a squeezing action, and reduction in a hammer mill is the result of impact. The action in rod and ball mills is a combination of impacting, squeezing, and wearing away by attrition, or rubbing. In ordinary crushing and grinding, the forces are applied at protruding points and are not distributed evenly throughout the rock. This is one of the features of comminution which theoretical considerations have not covered adequately. Another feature is the effect of impact velocities, which may greatly influence the total energy required for breakage. Impact velocities may vary from 1 or 2 fps in crushers to 10 or 20 fps for ball mills, and as high as 100 or more for hammer mills. However, according to the findings of geophysicists, the velocity of sound and of compression waves in stone is many times greater; the primary longitudinal compression waves travel at perhaps 15,000 fps, and the secondary transverse waves, which may reach the surface and cause cracks to form thereon, travel at perhaps 5000 fps. When a crack tip forms, the total surrounding stresses are concentrated in this tip, which rapidly extends throughout the rock particle.1 It would seem that the energy required to deform the rock beyond its critical strain, resulting in the formation of a crack tip, represents practically all of the energy required. Most of the primary crack tips presumably form on the surface. A simple method of designating feed and product size is necessary for the evaluation of crushing and grinding results. Taggart's suggestion2 that the 80 pet passing size be used has been found very practical. In this paper, the square opening screen size which 80 pct of the feed passes is designated as the feed size F, and the size which 80 pct of the product passes, or the product size, is designated as P. The reduction ratio at the 80 pet passing size, or F/P is designated as n. The feed and prod-

Jan 1, 1951

By L. W. Rowe, S. S. Cole

A laboratory bench scale unit is described whereby finely divided chlorinatable residues are held for a short period by a restraining bed of a coarse-grained ore of comparable composition to permit 'flash' chlorination of the high surface area powder. Results on chlorinating a finely divicled titania residue using rutile sand as the restraining bed showed that 92 pct utilization of chlorine and equivalent conversion of TiO2 in the fine residue to Ticl4 were attained with a dust loss of about 3 to 4 pct and an equivalent amount of ch1om'nation of the titanium values in the restraining bed. THE chlorination of coarse grain titaniferous ores; i.e., 20 to 150 mesh sizes in fluidized beds is practiced commercially.' One effective method found for the chlorination of a very fine material in a fluidized bed is the "flash" chlorination technique.' By charging the finely divided charge material into the bottom of an established bed of optimum sized rutile, ilmenite, titanium beneficiate, or even an inert material such as sand, chlorination occurs before the dust can pass through the bed. The major requirement is that the bed be less chlorinatable than the charge. In this process, the charge consists of a finely divided (between one and 50 µ) titania residue obtained from a titaniferous ore beneficiation process. The charge is blended with a 20 pct wt basis addition of 20 to 60 .mesh size low volatile coke (to minimize hydrogen chloride formation) and fed to the bottom of the bed with the chlorine gas stream. A ball check device is used to permit entry of the charge and prevent the bed from discharging when the gas flow is cut off. The coarse pieces of coke serve to reduce compaction in the charge mixture and permit sat]sfactory feeding without plugging or jamming. Feed and gas rates are adjusted to provide for complete reaction of the charge material, with 90 to 95 pct chlorine efficiency, and maintain a gas flow in the range of 0.4 to 0.6 ft per sec. The equipment consisted of a 3-in. diam mullite tube enclosed in an electric furnace. A bow-shaped 3/4 in. Saran tube fed the charge from an air lock type of feed hopper (to prevent the escape of the chlorine gas into feed bin) to the bottom of the bed. The charge rate was set by providing a small awiliary hopper below the feed hopper where a weighed portion of the charge could be discharged in a specified time interval into the main chlorine stream. The chlorine rate was measured by a flowmeter. The reaction tube was filled with a bed of Australian rutile ore, this restraining bed serving two purposes; i.e., it rapidly heated the fresh charge to the reaction temperature and retained the particles in the bed for a sufficient time for selective chlorination to proceed. Due to the fine particle size and open, porous structure of the charge, this reaction was almost instantaneous, and, with proper adjustment of feed and gas rates, very little fine material passed through the bed without reacting and only a very small fraction of the restraining bed was chlorinated. A typical reaction procedure was to fill the bench scale fluidized bed unit with a 1 kg charge of Australian rutile ore analyzing 95.5 pct TiO2, 0.6 pct FeO, a beach sand deposit 40 to 150 mesh size. When fluidized, the bed was approximately 2 ft in depth. The ore was fluidized with a 1:l mixture of chlorine and nitrogen at a rate of 0.40 ft per sec. This gas mixture was necessitated by the comparatively short bed depth; a larger unit, with a deeper bed, could be fluidized with chlorine alone without significant free chlorine loss. Feed rate was 8.0 g per min of beneficiated ilmenite analyzing 85.0 pct TiO2, 2.6 pct FeO plus 1.6 g per min of 20 to 60 mesh petroleum coke. The

Jan 1, 1962

By M. Garfinkle

LITERATURE data1'3 have suggested that the resolved shear stress necessary to initiate slip in a particular slip system in tungsten single crystals is strongly dependent on the tensile orientation of the crystals. This orientation dependence is contrary to the critical resolved shear stress criterion. Because the Schmid factors for slip on particular {110}( I l l), {112}(Ill), or {123}(111) slip systems in a (100) oriented tungsten crystal are virtually identical to the respective Schmid factors for a (110) oriented crystal, slip should initiate on the same system in both orientations, independent of the values of the critical resolved shear stresses. Thus, it would be expected that, for a critical resolved shear stress model to hold, the proportional-limit stress for these two orientations must be similar. The expectation of similar proportional -limit stresses has not been borne out by reported experimental observations. Rose et nl.' reported that the proportional-limit stress for crystals near the (110) orientation was 92,000 psi, approximately three times that of (100) oriented crystals (31,000 psi). Garlick and probst2 reported greater than a fourfold difference. Most recently, Beardmore and ~ull~ examined tungsten crystals with orientations along the edge of the standard stereographic triangle and showed a gradual increase in proportional-limit stress from 28,000 psi for the (100) orientation to 93,000 psi for the (110) orientation. Again, the difference was more than threefold. As has been emphasized by these investigators, the apparent gross variation of proportional-limit stress with orientation precludes any meaningful analysis of slip mechanics using a simple critical resolved shear stress relation. The purpose of the present communication is to report that, by use of refined strain-measuring techniques, the proportional-limit stresses for the (100) and (110) orientations are quite similar. Obviously, to determine accurately the value of the proportional-limit stress, which is by definition the stress at which nonlinear stress-strain behavior first occurs, the method of measuring strain must be capable of high accuracy. The investigators cited used as the measured elongation the crosshead displacement of the tensile machine. This method, however, introduces extraneous strain contributions from the load train. Because these extraneous contributions are considerable in the initial portions of loading (as has been well-documented in the literature) the proportional-limit stress measurement, dependent on the deviation of a line from linearity, can be significantly affected. During the course of a recent study, room-temperature strain was measured directly on several (100)-and (110)-oriented tensile specimens axially oriented within 5 deg of the desired direction. Each specimen was fitted with two foil strain gages mounted on opposite sides of the 0.08-in.-diam reduced section. The tensile specimens were loaded at a crosshead speed of 0.01 in. per min on an Instron Universal Testing Machine. Strain was monitored through an external bridge circuit in which the specimen strain gages comprised two arms. The bridge output was recorded against the applied load. A strain of 2 x 10"" could be readily detected by this method. Some specimens were subjected to conventional loading while the load for others was cycled between zero and continually increasing values until the first deviation from linearity was observed. The results of these methods were similar, and the stress at which nonlinear behavior occurred was considered the proportional-limit stress. The proportional-limit stresses measured in this study are approximately 5000 and 6500 psi for the ( 100) - and ( 110) -oriented crystals, respectively, a range of values almost a magnitude less than those reported previously. Representative curves are presented in Fig. 1. Because of the nature of the measurements made, the point of departure of a line from linearity is often debatable. While the load-strain curve for the (100) -oriented crystals made a relatively abrupt deviation from linearity, the curve for the (110) -oriented crystal was more gradual, so that the actual point of tangency may be at a lower load than reported. Thus, the difference in values for the proportional-

Jan 1, 1967

By M. E. Wadsworth, A. S. Peck

Infrared spectra disclose that oleic acid will not adsorb on the surface of pure beryl unless the mineral is first activated with HF. The adsorption of oleic acid on HF activated beryl is attributed to the hydrogen bonding of oleic acid monomers to F-surface bridging sites. These sites are formed as a result of the reaction of HF with chemisorbed water or surface hydroxyl. Characteristic infrared spectra of physically adsorbed oleic acid are correlated with flotation at different pH at temperatures of 25°C and 75°C for colloidal and -60-mesh particles. The usual mode of occurrence of the beryllium mineral beryl, Be3A12Si6O18, is in granite pegmatites as pale green, white, or yellow crystals. A dark green gem variety, emerald, may be found in the wall rock of pegmatite veins or in mica schist. The most important beryl deposits in the United States are in the Black Hills District of South Dakota, and several mountain ranges of central and west-central Utah contain potential commercial sources of this mineral. Several papers have been published concerning the flotation of beryl with different conditioning reagents. Kennedy and O'Meara1 reported on the essential use of HF in the flotation of beryl from different ores. The use of calcium hypochlorite to separate HF activated beryl from feldspar after flotation with an amine collector was made by Runke.2 A report on the activation of beryl with different cations using sodium oleate collector was presented by Viswa-nathan and co-workers.3 Flotation of beryl with sulfonate collector and different activating salts was employed by Fuerstenau and Bhappu.4 The role of the surface charge or zeta potential in the flotation of beryl with different collectors and salts was also extensively investigated.5"7 Previous infrared investigations by the authors8 disclosed that the orthosilicate beryllium mineral, phenacite, reacts with oleic acid to form a chemisorbed oleate monolayer on the mineral surface. Thus, chemisorption is the effective mode of adsorption in the phenacite-oleate acid flotation system. In contrast, neither chemisorption or physical adsorption of oleic acid occur on the surface of pure beryl as evidenced in this study. Efforts were subsequently made to determine the mechanism of activation and collector adsorption for HF pretreated beryl. EXPERIMENTAL PROCEDURES Large, white, terminated crystals of beryl were obtained from Madagascar through Ward's Co., Monterey, Calif. The crystals were broken into -60-mesh fragments in a porcelain mortar. Part of this prepared material was ground to a powder with a Fisher Grinder using an agate mortar and pestle. The powder was transferred to a graduated cylinder and mixed with 500 ml of distilled water. After 5 min, the colloidal suspension was separated from the settled particles by decantation and used as a stock for infrared test determinations. Adsorption tests were performed by first transferring 20 ml of colloidal suspension to 50 ml glass centrifuge tubes. The solids were then centrifuged out and the liquid was discarded. Then, the solids were pretreated with one molar HC1 or one molar HF. After 15 min, the solids were again centrifuged out. The mineral solids were washed free of residual acid by repeated centrifuging and water washing operations. Then, 20 ml of distilled water was added to the solids, and the resultant suspension was used for test purposes. In the preparation of mineral for flotation tests, 1.5 g of -60-mesh beryl was deslimed by removing the colloidal fraction by three successive water washes and decantations. Then, the coarse particles were pretreated with either one molar HCl or one molar HF for 15 min. The mineral was next washed free of acid with distilled water, and 20 ml of distilled water was added to the beryl prior to collector conditioning. The conditioning of mineral with collector was performed in 50 ml centrifuge tubes to which either NaOH or HC1 was added for pH modification. One drop (6.0 mg) of U.S.P. reagent grade oleic acid or 2 ml of 0.01 molar sodium oleate was added to the pulp and stirred at about 1500 rpm for 3 min at a

Jan 1, 1968

By J. L. Taylor, Pol Duwez

THE present knowledge of the Ti-Al system is limited to the portion of the diagram extending from pure aluminum to the intermetallic compound TiAl3' A preliminary investigation of the titanium-rich Ti-A1 alloys revealed that the solubility of aluminum in titanium is quite extensive and that the transformation temperature of titanium is raised by alloying with aluminum. From an X-ray diffraction study of alloys containing up to 75 atomic pct Al (TiAl3), the phase boundaries at 750°C were located as follows: the a titanium solid solution extends from 0 to about 36 atomic pct Al, a two-phase region exists from 36 to 46 atomic pct Al, and an intermediate phase of varying concentration extends from 46 to 62 atomic pct Al. The purpose of the present paper is to describe the crystal structure of this new phase. The alloys were prepared by melting in a helium arc furnace on a water-cooled copper plate, using a furnace essentially the same as that described in ref. 2. The titanium metal was of the iodide type furnished by the New Jersey Zinc Co., and a piece of pure aluminum (99.99 pct) was obtained through the courtesy of the Aluminum Co. of America. After melting, the samples, weighing approximately 5 g, were sealed in evacuated fused silica tubes and homogenized for 4 hr at 1000°C. The specimens were then held for 10 days at 750°C and rapidly cooled to room temperature. Filings were then taken from each sample, sealed in evacuated vials, annealed at 750°C for 4 hr, and rapidly cooled to room temperature by quenching the vials in water. Powder diffraction patterns were obtained with a 14.32 cm diam camera, using K copper radiation filtered through a nickel foil. Structure Determination The X-ray diffraction patterns of the alloys containing 46, 50, 55. 60, and 62 atomic pct Al were identical, except for a slight shift in the positions of the reflections. The alloy containing 55 atomic pct Al was chosen far structure determination. The powder pattern contained 37 reflections and the a, doublets were well resolved in the back-reflection range. The relatively small number of reflections and the sharpness of the pattern were considered as an indication of a simple atomic arrangement in a relatively small unit cell of high symmetry, and an attempt was therefore made to solve the structure without single crystal work. A satisfactory fit was found on a large scale Hull-Davey tetragonal chart for an axial ratio of approximately 1.45. All the reflections were readily indexed and the tentative unit cell dimensions were computed to be a - 2.81 kX, c - 4.07 kX, and c/a = 1.448. The indices (hkl) and the values of d and sin2 A are given in Table I. The agreement between observed and computed sin2 is quite satisfactory. Once the unit cell was known, the details of the structure were readily found. First, it is obvious that the number of atoms per unit cell cannot exceed two, since more than two atoms would lead to a density greater than that of pure titanium. The next logical step is to assume that the two atoms occupy the corner and center of the cell. If this assumption is correct, the structure may be described as follows: space group D1, - P 4/mmm; one titanium atom in a: 000; one aluminum atom in d: 1/2 1/2 1/2: and no extinctions. Assuming this structure to be the correct one. intensities were computed by means of the usual equation: sin2 0 cos 0 where F is the structure factor; 8, the Bragg angle. and p, the multiplicity factor. The calculated values of intensities are compared in Table I with the visually estimated intensities. The agreement is quite satisfactory and the structure is therefore confirmed. The crystal structure of TiAl just described is the same as that of AuCu ordered. The AuCu ordered

Jan 1, 1953

By Michael Hoch, Walter C. Wyder

The Cr-Ni system was investigated between 900" and 1400°C and 45 to 100 wt pct Cr. The Cr-Fe system was studied between 1200" and 1385°C and 70 to 100 wt pct Cr. Powdered alloy compacts were produced from high-purity chromium, nickel, and iron poulders. The compacts were sintered in an argon atmosphere and then machined into small, 0.030-in. diam cylindrical specimens. These were then heated by induction in the high-temperature X-ray diJfractwn camem under a helium atmosphere and irradiated with copper Ka radiation. From the results of the high-temperature X-ray diffraction patterns it can be concluded that chromium exists only in the bcc form, and that no eutectoid reaction occurs in the Cr-Ni system in the temperature interval 900" to 1400°C. The X-ray diffraction patterns are in agreement with the Cr-Ni phase diagram of Bechtoldt and Vacher, and Williams. It is possible that chromium undergoes one or more allotropic transformations between room temperature and its melting point. Most of the evidence was gathered from the study of the chromium-rich ends of the Cr-Ni and Cr-Fe binary systems. All the investigations rely heavily on indirect methods: thermal analysis, electrical resistivity, and micro-graphic inspection of quenched samples. In the present work the chromium-rich side of the Cr-Ni and Cr-Fe binary systems is studied, using high-temperature X-ray diffraction techniques. Thus the coexistent structures can be observed directly at temperature. PREVIOUS INVESTIGATIONS ~ansenl shows two possible Cr-Ni phase diagrams 1) a simple eutectic [Liquid =a-Cr(bcc) + y-Ni(fcc)], 2) a eutectic [Liquid+ P-Cr(fcc) + y-Ni(fcc)] and euctectoid [p-Cr(fcc) ==a-Cr(bcc) + yNi(fcc)]. The latter is based on the work of Grant and co-workers.2 Stein and rant'" placed the eutectoid composition at 1215°C and 68 wt pct Cr; their diagram is shown in Fig. l. Similarly, Sully3 also discusses the two proposed diagrams. Price and Grant,4 working on an isothermal section of the Cr-Ni- Fe ternary diagram at 1300°C, made no mention of a ß-chromium phase (fcc). The ß-chromium structure was also reported by Misencik 5 in the Cr-Nb(Cb) system. Williams 6 studied the Cr-Ni system using precipitation techniques. His work included the chromium and nickel solvus lines below 1250°C. He reported no ß structure nor any eutectoid reaction. The solvus lines which he produced were very similar to those of Jette, et al 7 below 1000°C. He concluded that Grant and coworkers2 had misinterpreted their data, and that effects which they had reported were actually due to the rapid and extensive precipitation of the nickel-rich phase. Bechtoldt and vacher8 used sintered powdered alloy compacts in an attempt to find the ß phase. They used metallographic and room-temperature X-ray diffraction techniques and also failed to detect the fcc chromium phase. They attributed the low slope of the chromium solvus at approximately 1150°C as the probable cause of the abundant precipitation obtained. They also stated that the precipitation was the probable cause for the abrupt discontinuities in the thermal analysis and electrical resistivity data which led to the conclusion that a eutectoid reaction existed. Fig. 1 also shows their diagram. Grigor'ev and coworkers9 investigated the polymorphic changes of chromium in the Cr-Ni phase diagram. Using standard microstructural and thermal analysis techniques, they published a phase diagram which contained five polymorphic forms of pure chromium and four eutectoid reactions. They stated that the room-temperature bcc chromium transformed to a fcc structure at 930°C, to bcc at 1300°C, to hcp at 1650oC, and then finally to bcc at 1830°C. The eutectoid reactions were placed at 850°, 960o, 1140°, and 1220°C. Also to be noted was the fact that the uppermost phase, the bcc above 1830°C, was different from the upper phase of Grant and coworkers,2 the fcc above 1840°C. Their eutectoid reactions, temperatures and compositions fall

Jan 1, 1963

By D. A. Stevenson, R. A. Burmeister

Single-phase silver antimony telluride has been prepared by zone-melting techniques using initial compositions of A new phase appears upon prolonged annealing of this material, but the reaction does not appear to be a simple eutectoid decomposition. A complete analysis of the phase equilibria is complicated by the slow kinetics involved. The Hall coefficient, magneto -resistance, electrical resistivity, and Seebeck coefficient are all sensitive to the presence of second phases. The low Hall mobilities measured for single-phase material indicate that the usual band theory is inadequate to explain the observed transport properties in the system. Density anomalies of up to 2.5 pet between measured and theoretical density were observed but are not conclusive evidence for a defect structure. COMPOSITIONS in the Ag-Sb-Te system have been studied previously by several investigators.1"18 The interest in this system arises from potential thermoelectric applications of alloys on the Ag2Te-Sb2Tes vertical section. Although the composition corresponding to the formula AgSbTe, has received most attention, it has been found to consist of more than one phase.7'8 A thorough understanding of the properties of this heterogeneous material has been impeded by both lack of knowledge of the properties of the homogeneous phases comprising it and the problem of analysis of transport properties in an inhomogeneous system. The present work describes the preparation of the homogeneous ternary phase and the corresponding electrical properties of both homogeneous and heterogeneous material. MATERIAL PREPARATION Most specimens for this investigation were prepared by encapsulating the elements in evacuated quartz tubes after which they were melted and homogenized in the liquid state. The ambient temperature was then dropped to 500°C and the resulting ingot zone melted. The growth rate, width of zone, stoichiometry, and number of passes have an effect on the resultant microstructure. Grains several centimeters in length were easily produced by this method. Other solidification techniques were also used, including uniform slow cooling of the entire specimen and rapid freezing. The microstructures of specimens produced by these techniques frequently differed appreciably from similar compositions prepared by zoning. A variety of nonequilibrium microstructures characterized by long needlelike particles resulted from the rapid-freeze method. PHASE EQUILIBRIA The lack of information on phase equilibria is a major difficulty associated with a comprehensive study of the Ag-Sb-Te system. Considerable confusion has resulted from the use of the formula AgSbTe, to identify the cubic ternary intermediate phase even though it has been established that material of this stoichiometry normally contains AgzTe as a second phase.798 In this paper, silver antimony telluride will denote the cubic ternary intermediate phase comprising the major portion of AgSbTe,. The term "single phase" will denote material which consists of only the cubic phase (as evidenced by metallographic examination and X-ray diffraction) and the term heterogeneous will describe multiphase material containing AgzTe, SbzTe3, or other phases in addition to the cubic phase. The single-phase material may actually contain a variety of inhomogeneities—gradual changes in composition on a macroscopic scale, localized fluctuations in composition, clusters or other products of early stages of the precipitation process, and a variety of point and line defects—all of which will not be detected by the present techniques for determining homogeneity. Single-phase material has been prepared from compositions close to 59 mole pct SbzTe3 21 (this notation refers to the location on the Ag,Te-Sb2TeS vertical section). Zone widths of =2 cm and growth rates 51.2 cm per hr were used. The single-phase region at elevated temperatures extends off the vertical section to a composition which can be expressed approximately as Ag,,SbzgTes,.9 The latter two compositions as well as AgSbTe, are represented in the conventional Gibbs triangle in Fig. 1. It is not presently possible to ascribe exact values to the limits of the single-phase region on a given isotherm or vertical section due to the extremely

Jan 1, 1964

By J. C. Wilhoit, J. B. Cheatham

Although an oil well is a long cylindrical hole with an irregular bottom, it appears likely that the nature of the stress concentration at the bottom of the hole can be ascertained from an analysis of the stresses around a short cylindrical cavity with rounded comers and smooth bottom. Such a cavity is studied primarily because it leads more readily to a solution to the problem by the use of stress functions. In this paper the stress distribution around a short cylindrical cavity subjected to bit loading, overburden and drilling-fluid pressures is determined by means of an analytical solution which approximately satisfies the boundary conditions of the problem. From this solution the stresses at the corner of the hole are calculated to be about 35 per cent lower than comparable results obtained by photoelastic and relaxation analyses. This difference is apparently due to the large radius of curvature at the corner of the cavity in the present analysis. Since good agreement is obtained between the results of this analysis and the stresses calculated for a similar loading on a semi-infinite elastic solid, it is concluded that the bit action in the region near the center of the hole is not appreciably affected by the presence of the sides of the hole. INTRODUCTION Much has been written concerning drilling "under down-hole conditions" and pertaining to the stress distribution in the rock at the bottom of an oil well.1-5 For example, it is known that identical rocks can be drilled more rapidly at the surface than under subsurface conditions of pressure and stress.6 Information on the behavior of rocks under loading can be obtained from triaxial test data.7-9 From such tests it is found that rocks exhibit brittle failure when the confining pressure and pore pressure are equal, but the mode of failure may change to ductile as the difference between the confining pressure and the pore pressure is increased. Brittle failure implies that there is very little permanent deformation before fracture, whereas ductile failure indicates that permanent deformation takes place before fracture. Some rocks are ductile at differential pressures of 5,000 psi , but other rocks are brittle even at differential pressures of more than 50,000 psi. Cuttings embedded in mud at the bottom of the hole may act as a plastic mass which the bit teeth must penetrate in order to attack the virgin rock below. 10 During drilling, the mud pressure acts as the confing pressure, and in many cases the difference between the mud pressure and the formation pore pressure IS sufficiently low so that the rock most likely falls in a brittle manner. Very little penetratiori by the bit teeth may be required for brittle failure of the rock. For these brittle materials, the elastic stress distribution is of practical interest in determicing the possible effects of a given loading on failure. The penetration of bit teeth into ductile or plastic rock has been analyzed previously ll,12 and will not be considered further in the present work. During actual drilling, many teeth act at various points ovel an irregular hole bottom. In the present analysis a solution is obtained for an idealized problem of determining the elastic stress distribution caused by only one tooth acting alone at the center of a smooth hole bottom. It obviously is not proposed that any driller should put a special one-tooth bit on bottom, but it is hoped that extensions of the simpler problems can eventually lead to a better understanding of more complex actual drilling phenomena. To obtain a solution by the use of stress functions, a short cylindrical cavity with rounded corners and smooth bottom is studied. The solution to this problem should give an insight into the nature of the stress concentration at the bottom of an oil well which is in reality a long cylindrical hole with an irregularr bottom. After selection of the proper curvilinear co-ordinate system, the stress functions are expressed as series of Legendre polynomials. The coefficient of each term in the series is then determined to satisfy the boundary conditions at a finite number of discrete points on the boundary of the cavity with least square error. Results are compared with bottom-hole stresses obtained by photoelastic and relaxation analyses.

By T. Yamasaki, S. Usu

Adsorption of potassium dodecyl xanthate from aqueous solutions on artificial and natural zinc sul-fides was studied by means of infrared absorption spectroscopy. The adsorption species and their stabilities varied depending upon the conditions of the surfaces of the zinc sulfides. Adsorption products were analyzed after they had been removed from the surface with pyridine. Zinc dodecyl xanthate was confirmed for surfaces covered with oxidized products while dodecyl alcohol and carbon disulfide were detected on purified or weakly oxidized surfaces. It is of great significance to study the adsorbed species of collectors on the surface of sulfide minerals with a view of increasing our understanding of the mechanism of flotation. Compared with other methods, infrared spectroscopy promises to be of particular value in determining the nature of the species adsorbed at the mineral/solution interface since it gives direct identification of the adsorption products; and, furthermore, it has an advantage over the electron diffraction technique which has been employed in similar problems in that infrared spectroscopy utilizes radiation of very low energy so that the probability of changing the unstable surface films during examination should be correspondingly decreased. Infrared spectroscopic studies concerned with sulfide mineral flotation have already been reported by several investigators on zinc sulfide-hexanthi01' and lead sulfide-xanthate2,3 systems. The purpose of the present study was to investigate the adsorption of potassium dodecyl xanthate (KDX) at the zinc sulfide-aqueous solution interface by means of infrared spectroscopy. Rarely is the pure zinc sulfide-xanthate system encountered in practical flotation since zinc sulfide is usually activated with small amounts of certain metallic ions, such as copper ion. However, investigation of the pure zinc sulfide-xanthate system appeared promising for studying the fundamental mechanism of the adsorption of xanthate on the sulfide minerals. Although there are few cases where xanthates having an alkyl radical higher than hexyl are employed in practical operations, KDX was used in this study since it was reported that the adsorbability of xanthate on zinc sulfide was increased with an increasing number of carbon atoms of the alkyl radicals.4 EXPERIMENTAL Samples: Experiments were made on artificial zinc sulfide and natural sphalerite. The artificial zinc sulfide was obtained commercially from Dainihon Toryo Co. as a fluorescent pigment of the purest grade A. Considering the fact that there are often cases where the surface of sulfide minerals are apt to be oxidized so that the adsorption of xanthate at sulfide mineral/solution interfaces may be substantially affected, two different procedures were employed in preparing the samples of the artificial zinc sulfide. 1) The zinc sulfide was washed with hot concentrated ammonium chloride solution (300 g per liter) to remove surface oxidation products and followed by washing with deoxygenated hot distilled water by centrifugation until the constant electric conductivity of the solution was reached (Sample I) (specific conductivity of the solution = 8.6 x l0-5 ohm-' cm-'). 2) Sample I in a suspension of 5% by weight was oxidized with a 1.5% aqueous hydrogen peroxide solution for about one hour, and then was washed with distilled water as mentioned in the preparation of Sample I (Sample 11) (specific conductivity of the solution = 1.0 x 10"4 ohm-' cm-'). Both Samples I and II were dried under vacuum and stored in a vacuum desiccator. The specific surface areas of the samples determined by the BET method using carbon dioxide gas were 19 x l04 cm2 per g for Sample I and 35 x l04 cm2 per g for Sample 11. In addition, a sample which was derived from Sample I was also prepared. In 40 ml of 3 x 10-2 M zinc sulfate solution was placed 2 g of Sample I for 2 hours, washed with distilled water by successive centrifugation until no sulfate ion could be detected and dried under vacuum (Sample 111). The purpose of this treatment will be described later. Natural sphalerite was carefully selected by hand, ground in a mechanical agate mortar for about 40

Jan 1, 1965

By J. C. Sarace, S. M. Sze, C. R. Crowell

Thin films of tungsten 077 n-type germanium, silicon, and gallium arsenide were obtained by reacting tungsten hexafluoride with the semiconductor surface in an argom atmosplrere at temperatures between 325° and 400° C. Capacity-voltage, current-iloltage, and photoelectric measurements were used to investigate the characteristics of the tungsten -semiconductor diodes thus Produced. The junctions are shown to he very close to ideal Schottky barlp/ers with barrier heights measured with respect to the Fermi energy of 0.18, 0.65, and 0.78 1.1 jar W-Ge, W-Si, and W-GaAs, respectively. The electrical properties of the W-Si interface show no deterioration when heated to 1000°C in dry forming gas for 5 min. A theoretical value of the Richardson constant, A, appropriate to the semiconductor-hand structure has been used in evaluating the current-voltage characteristics. ThE W-Si surface-barrier diode was initially proposed for investigation because the eutectic temperature with silicon (1400°C) is much higher than that in the Si-Au system 1370°C).1 This would permit more flexibility in heat treatment and possibly provide greater reliability at elevated temperatures. The lower work function of tungsten (4.54 ev)2 compared with that of gold (4.78 ev)3 also suggested that a lower barrier height would be obtained with tungsten and hence a lower forward bias and lower minority carrier injection ratio for a given current density. The investigation was extended to include the characterization of W-Ge md W-GaAs surface-barrier diodes. The tungsten films have been produced by reacting WF6 with germanium, silicon, and GaAs surfaces in an argon atmosphere at temperatures from 300° to 500°C.4 This process is a very satisfactory alternative to the relatively difficult process of evaporating tungsten films in vacuo. To ensure an adequate electrical characterization of the tungsten-semiconductor interface, three types of barrier-height measurements have been performed. The mutually consistent results obtained lead to the conclusion that the tungsten-semiconductor junctions are indeed of the Schottky type. EXPERIMENTAL PROCEDURE The apparatus used for producing tungsten films is shown schematically in Fig. 1. It consists of an argon carrier gas line to which metered amounts of tungsten hexafluoride can be rapidly added. The mixture passes through a heated reaction tube containing the semiconductor slices and is exhausted to a hood. The argon is purified by passage through a 6-in. column of titanium turnings maintained at 800°C. The tungsten hexafluoride dispensing arrangement was designed by V. C. Garbarini and W. R. Bracht.4 A measured amount of liquid tungsten hexafluoride is injected into the argon stream and vaporizes. The mixture passes through a sodium fluoride absorption cell to remove traces of hydrogen fluoride, then into the nickel reaction tube. The tube is 12 in. long with an inside diameter of 1/2 in. The center section is heated by a furnace of the self-supporting wire-filament type. It was chosen for its rapid heatup and cooling. The wall thickness is 10 mils except for a 3-in. hot zone which is 30 mils thick to reduce thermal gradients along the length. The samples to be coated are placed on a sapphire plate and centered in this section. The tungsten is deposited with the following sequential steps: the loaded reaction tube is flushed with argon at a rate of 500 cu cm per min and heated to 370°C. Then 0.45 g of liquid tungsten hexafluoride is injected into the argon stream. The samples are held at this temperature for 2 min. The tube and samples are then cooled to room temperature and the samples removed. Tungsten films were grown on (111) faces of silicon and germanium polished with Linde A abrasive and lightly etched with HF-HNO3. The GaAs surfaces were (100) faces chemically polished with a H2SO4-H2O2 solution. The films have typical sheet resistances of 0.2, 8, and 15 52/0 when grown on germanium, silicon, and GaAs, respectively. After the tungsten deposition, ohmic contacts were alloyed on the back surfaces of the wafers. Ohmic contacts to the germanium and silicon were obtained by alloying Au-Sb at 370°C, ohmic con-

Jan 1, 1965

By D. W. Mitchell

INE operators have a choice of several classi- fications of mechanical loaders. Within each classification there are many types and makes available. Table I lists loaders on which manufacturing data as to operating characteristics are available. This paper discusses the conditions met in a mine as they affect these characteristics. It is assumed that management will provide satisfactory engineering, supervision, power, maintenance, loader service within a concentration of workings, and a balanced working cycle served by balanced production equipment. These are factors which affect the optimum operational efficiency of a loader. The conditions which determine loader choice are: Height of Vein—The main limitation to the use of a loading machine is the height to which the vein is mined; i.e. a loader cannot be higher than the place in which it is to work unless rock is taken. Since there are several low vein loaders being developed and successfully applied, this does not appear to be a requisite of an efficient mining operation. The maximum useable height of a loader should be equal to the working seam thickness less a working clearance for travel and operation. A—Working seam thickness = B—Artificial roof support thickness = C—Safe headroom below support = (generally about 6 in.) D—Height of roadway = E—Total of B + C f D =_______________________ A—E = maximum loader height F—Distance from top of roadway to bottom of roof support = G—-Height of loader's coal or ore line (see Table I) = H—Difference between F and G = Size of largest broken particle + a safety factor of 2 to 3 in. H is a loadability factor and is important because the clearance over the top of the conveyor chain to the roof or roof support must be greater than the largest broken particle or else the particle and/or the conveyor chain may be broken. From this standpoint it might be desirable to use a lower height machine even though it may have a smaller capacity, though this is warranted only when increased lump realization offsets lower man-day and machine pro- duction rates. Ivan Given has stated, "Loading should not be handicapped by securing lump only to grind it in the breaker." Width of Working Place—The width of a working place is the greatest width a place may be safely driven. The working place includes not only the face but the roads the loader must tram when mucking out more than one face. The width is limited by the system of placing props with sufficient room between them for the machine runner (about 2 ft). Maximum width of place in which the loader will operate at maximum efficiency is determined by the type of mounting. With track-mounted loaders, maximum width of place is determined by maximum angular swing of the loading head and proximity of track to the face. The use of double track leads to added expense and operational complexities. Several mines that have used double-track systems have found that a substantial saving is made by the use of machines that load out wider places. Trackless loaders may work in rooms of any maximum width. Minimum width of place in which a loader will operate is determined by width of the loader plus a safe movement area for the operator. Often a barodynamic study of the mine will show that either an increase in room width and/or a change in the system of propping may be made by using roof bolting, full or part-width room timbers, removable aluminum I-beams at the face, etc. Increasing the width of a working place will, in general, increase loader operating efficiency by providing more material per fall and by decreasing the preparation and moving time per ton. A change in the system of propping that would increase the minimum width might permit the use of a higher capacity loader with safer working conditions. Proper face tirnbering is necessary to insure safety from roof falls and because the speed of preparation and loading generally increases since the men work with greater assurance and there is less chance of kicking out posts. Maximum Room Length—With the exception of the scraper and the duckbill, the limit to room

Jan 1, 1952

By Gearge C. Wallick

This paper describes experimental procedures for the calibration of capillary tubes to be employed as comparison standards in gas flow-rate measurements and considers several types of flow which were observed in the calibration of seven capillaries of varying diameter. It is shown that under certain experimental conditions deviations from Poiseuille flow are observed which may not be attributed to turbulence, and that this flow behavior may be described empirically by introducing the classical kinetic energy correction into the Poiseuille equation. INTRODUCTION Techniques have been developed for the measurement of the volume rate of flow of nitrogen through glass capillary tubes as a function both of the differential pressure applied and of the absolute gas pressure. Capillary tubes thus calibrated are frequently used to determine the rate of gas flow through a porous sample in the measurement of slip-corrected permeability. The purpose of this note is to describe the calibration procedure and to discuss some observed deviations from the ideal flow behavior predicted by the Poiseuille equation. EXPERIMENTAL PROCEDURE The volume flow rate of nitrogen as a function of the differential pressure was determined experimentally for a series of seven glass capillary tubes of varying diameter. As is shown in Table 1, the capillary radii in the series ranged from 0.004 cm for K1 to 0.07 cm for K7. The capillary tubes were cut from capillary tubing purchased on special order from several manufacturers and the bore radii were determined by microscopic measurement. The gas volume flow rate was measured for each capillary tube at a series of differential pressures in the range from 0 to 3 cm of mercury and at average absolute pressures of 1.0, 2.5, 4.0, and 6.0 atm. Experimental Procedure for Small Capillaries A schematic diagram of the apparatus used to calibrate the three smallest capillaries, Kl, K2, and K3, is shown in Fig. 1. By means of a standard Ruska Instrument Corp. volumetric pump (equipped with a synchronous motor drive, two stages of gear reduction, and a quick-change lathe gear box), nitrogen was displaced by mercury from a metal cell immersed in a constant-temperature bath. The pump, the line to the cell, and a part of the cell were filled with mercury. As mercury was forced into the cell, it displaced nitrogen which passed through the capillary tube being calibrated and into a large container with a volume of approximately 30 liters. Since the volume of this container was very large compared to the total volume change of the gas resulting from the action of the displacement pump, the absolute pressure of the system was essentially constant throughout any one run. Prior to any determination this pressure could be set at any desired level by means of an auxiliary nitrogen cylinder. The differential pressure across the capillary was measured with a Meriam Red Oil manometer and recorded together with the corresponding pump displacement rate which was equal to the gas flow rate. By repeating this process for a series of different displacement rates at each of the absolute pressures of interest, the data necessary for plotting the calibration curves for capillaries K1, K2, and K3 were obtained. All calibration measurements were made in a constant temperature room to minimize errors due to temperature fluctuations. In each run. the flow was maintained until equilibrium was established. Experimental Procedures for Large Capillaries The flow rates through Capillaries K4, K5, K6, and K7 required to obtain the pressure drops desired were larger than could be obtained from available constant displacement pumps, and it was thus necessary to employ an alternative calibration procedure. The arrangement of apparatus used in calibrating these capillaries at an average pressure of one atmosphere is shown in Fig. 2. Each capillary was calibrated by flowing nitrogen through the capillary and recording the differential pressure as a function of the flow rate. 'The differential pressure was measured with a Meriam Red Oil manometer and the flow rate was determined by collecting the effluent gas over water at constant pressure and recording the increase in volume per unit time after an equilibrium had been established. Appropriate corrections were made for the

Jan 1, 1953

By J. E. Owen

A model of a porous body is presented in which the pore space consists of a system of voids and interconnecting tubes. Relationships between porosity and resistivity formation factor are determined partly by calculation, partly by experiment. Con triction effects characteristic of the model are shown to be sufficient to account for high formation factors. It is shown that constriction may be combined with moderate amounts of tortuosity to give model pore systems exhibiting to a first approximation porosity and resistivity properties similiar to those of natural porous bodies. INTRODUCTION The relationship between the electric resistivity of a fluid-filled porous body and the geometry of its pore space is so complex that the calculation of the resistivity of a natural porous rock is a practical impossibility. Both the resistivity of a body and its porosity are measurable quantities, however, and previous successes at relating them have been reached by an empirical approach. Efforts at obtaining theoretically derived formulae relating them have generally been unsatisfactory. One of the reasons for this may lie in the pore geometry that has been assumed. THE TORTUOSITY CONCEPT A Parameter called the formation factor is useful in dis-cussing the resistivity of a fluid-filled porous body. This parameter is the ratio of the resistivity of a fluid saturated porous body to the resistivity of the saturating fluid. Formation factors are often available from measurements on cores or from electric logs, and many attempts have been made to correlate formation factors and porosities of geological formations. Whenever a successful correlation is found, the engineer working with electrical logs has a useful tool for the determination of porsities of pay section?. One of the more successful formulae applicable to these correlations is the familiar equation empirically obtained by Archie.' which F is the formation factor. $ is the porosity, and rn is an exponent called the cementation factor. When the for- mula applies, the cementation factor usually is found to be between 1.3 and 2.2. The values for formation factors experimentally obtained are often higher than simple pore geometry would lead one to expect. In an effort to account for such high values certain formulae have been derived based on a so-called "tortuosity concept." In deriving these formulae a synthetic porous body is usually assumed in which the solid material is an electrical non-conductor. and in which the pore system consists of three sets of fluid-filled tubes of uniform diameter connecting opposite faces of the body which, for convenience, is considered to be cubical in shape. The three sets of tubes account for the whole of the effective porosity of the body, and usually, it is specified that they do not interconnect. By considering that the pore tubes are not straight but tortuous, their resistance to the flow of electric currents can be made as high as needed to explain high formation factors. Such an explanation has some basis in fact, but it appears that the tortuosity concept is often incorrectly applied when other factors are largely responsible for observed high resistivities. Recently, Wyllie and Spangler have recognized that tortuosity as calculated by conventional formulae has little if any physical significance.' RESISTIVITY AND THE CONSTRICTION CONCEPT Any explanation of high formation factors which depends solely on tortuosity of uniform pore paths necessarily ignores the effect that variations in the cross-sectional area of the conducting paths have on the resistivity of a body. Although, as previously pointed out, the calculation of such paths for an actual body is impossible, it will he shown that a synthetic pore network can be devised which will yield to analysis, and lead to results in agreement with the experimental data represented by Equation (1). The porous body to be considered is assumed to be homogeneous and isotropic or, for present purposes, identical in its characteristics in the three directions parallel to its coordinate axes. It will he assumed to be built of identical unit cubes, each of which contains a single pore network connecting all faces of the unit cube. A unit of such a pore network is shown

Jan 1, 1952

By T. H. Alden

BaCKOFEN, Avery, and Turner have shown"2 that the large tensile elongation in superplastic metals is correlated with a high strain-rate sensitivity of the flow stress. At present, the reported superplastic materials are complex eutectic or eutectoid systems with a large solubility range in at least one of the terminal phases. This note describes the behavior of a solid-solution alloy (at 22°C)5,= Sn-1 wt pct Bi, which has been processed in order to produce a fine grain size.2,3,6 The chill-cast, 1/2-in.-diam ingot was homogenized for 1 week at 130°C, extruded at O°C to 0.083 in. diam, quenched within 15 sec, and stored at -196°C. Samples were warmed to room temperature 16 hr prior to mechanical testing or metallography. During tension testing at room temperature, it was found that the flow stress reached a steady-state value which depended on the imposed elongation rate and was independent of strain to an error of about ±5 pct. In Fig. 1 is plotted log a vs log e where a is the measured steady-state flow stress and i the imposed strain rate. These quantities were calculated using "instantaneous" area and length to an accuracy in the sample diameter of ±5 x 10-4 in. The log a-log 1 curve is sigmoidal. The strain-rate sensitivity, which may be given by d log old log 6,' has a maximum value of about 0.48 between i = 5 x 10-4 and 5 x 10-3 per min. The value is well above frequently measured values for metals at creep temperatures (0.1 to 0.3) and comparable to the rate sensitivity for Pb-Sn solder: which has been shown to be superplastic.2,3 Note, however, that below or above the strain rate of 10-3 per min the rate sensitivity of Sn-1 pct Bi decreases and at high rates is more nearly like that of conventional metals. The grain size of this material is relatively uniform, Fig. 2, and may be specified by an average linear intercept i.' Grain growth occurs at room temperature but after 16 hr, L = 5. Tensile elongation was measured on finer-grain-size specimens, which had been at room temperature 15 min prior to testing. After 15 min, the samples are fully recrystallized. Three strain rates, 5 x 10-3, 10-2, and 2.5 x 10-2 per min, were employed for which the measured elongations varied from 425 to 500 pct. These results would seem to justify the term superplastic, applied to this alloy, although the ductility is inferior to that of the Sn-Bi eutectic.3 The Sn-1 pct Bi alloy, which is a single-phase material with equiaxed, fine grains, is an excellent one on which to make an estimate of the Nabarro-Herring9 creep rate. According to the model of Avery and Backofen2 a substantial part of the deformation must be accomplished by N-H creep in the strain rate-temperature regime in which d log a/d log 1 = 0.5. The maximum N-H creep rate is given by9 L2kT where v is the atomic volume and D the self-diffusion coefficient. For this test v = 3.O x 10-23 cu cm, D = 1.4 e-23,300/RT = 2.5 x 10-l7 sq cm per sec,7 L = 5x 10+4 cm, T = 295°K, and a = 2000 lb per sq in. = 1.38 x 10' dyne per sq cm. The calculated creep rate is E= 4 x 10-11 per sec which may be compared with the measured creep rate at this stress, Fig. 1, of 6.7 x 10-5 per sec. The N-H creep rate is seen to be negligible. Earlier consideration2 of this question was in error because incorrect diffusion coefficients were used. The assistance of R. R. Russell in metallography is acknowledged gratefully.

Jan 1, 1967

By Louis P. Gignac

INTRODUCTION The selection of trackless equipment for underground mining can be a complex engineering problem due to the wide range of equipment sizes and operating modes. Computer simulation is particularly useful in estimating the performance of different systems for specific rock- handling problems. A hybrid simulator, incorporating some of the features of deterministic and stochastic simulation, was developed in order to handle not only queuing delays at loading and dumping points, but a1 so traffic interferences. Any number of vehicles can operate in any one of four basic modes (LHD, LAF, FEL, HFC) in parallel or in series. If the units use common roadways, loading and dumping points, certain operating delays will occur and be registered by the simulator and thus give a better evaluation of the marginal productivity of each additional unit. Based on a typical layout of drawpoints and ore passes for the sub-level stoping method, productivity and operating costs of different rock-handling systems will be examined. 1 . COMPUTER SIMULATION Numerous applications of computer simulations are reported in the 1iterature for various mining problems. Depending on the complexity of the system to be studied, simulation models were conceived with different degrees of sophistication. Three different types of simulators are generally recognized: stochastic, deterministic, and hybrid. Stochastic or Monte Carl o simulation randomly generates items, transactions, or events from some population defined by a frequency distribution and produces some expected future situations. Because this type of simulation is governed by the input of probability distributions, it requires a detailed knowledge of the system to be simulated; it implies expensive and time-consuming studies and reports to gene- rate this input information. A major short- coming of stochastic techniques is in new equipment evaluation, where the lack of data is unavoidable, and in new system design where the conditions are outside the range of the known historical behavior of the equipment. However, probabilistic simulation is almost essential for the study of cyclic queues and traffic problems. Deterministic simulation studies a system by generating performances on the basis of the mechanical capabilities of the vehicles and the physical limitations of the mining scheme. It is based on the engineering principle that the engine converts its energy into a rimpull at the wheels, which is in turn opposed by the rolling resistance and the grade of the ground; the machine is accelerated or decelerated until the tractive and resistive forces are in equilibrium, at which point it moves at constant speed. The information required by this technique, such as rimpull charts and equipment weight, is readily available from equipment suppliers. However, equipment performance at the mining site is also dictated by human and environmental conditions and changes with time and usage. For this reason, deterministic simulation generally overestimates the system capabilities; these must then be adjusted by efficiency factors based on observation and experience. Recently, hybrid simulation models, using both stochastic and deterministic techniques, have been built with some of the events generated stochasticly and others being deterministic. This compromising approach originated at the Pennsylvania State University. O'Nei1 (1966) designed a simulator for truck-and- shovel operations that allows for transportation from multiple faces to mu1tiple destinations. Each truck performs according to its mechanical capabilities while its loading and dumping time and its load fluctuate according to specific probability distributions. A major advance in the simulation of mining systems is due to Sanford's model (1969) of underground coal mining operations. The originality of the model is in the use of an Executive System Control which sets up the initial system, de- fines the operation conditions , and instructs continuously four sub-assemblies representing shuttle cars, trains, continuous miners , and conveyors, which in turn generate a feed-back of their movement to the Executive Control System for further instructions. The model has enough flexibility to simulate simultaneous and sequential jobs that characterize any dynamic system. Sanford's early work evolved slowly to what is known today as the Under- ground Materials Handling Simulator for coal mining (Manula, 1974).

Jan 1, 1981

By A. L. Bailey, B. A. Landry

All who are even generally aware of the tremendous rate of increase in coal washing operations must realize the growing importance of the float-and-sink test. I believe it is conservative to estimate that, in the past decade, the dollar volume of float-and-sink testing has increased tenfold. It is a simple matter of economy, then, to examine the factors that determine the cost of adequate float-and-sink testing. When the Coal Preparation Section of the Bureau of Mines entered upon a greatly expanded program of such work in connection with the synthetic liquid fuels investigations, it seemed advisable to examine these factors experimentally. The principal consideration that differentiates float-and-sink test sampling from general purpose sampling, is that the original particle size must be preserved. Therefore, the total cost of the test will be directly affected by any standard that might be proposed to limit sample bulk reduction at any given particle size. For this reason, the relationship of sample size to variability of results was the first factor to be studied experimentally. Of course, the matter is rendered complex by the circumstance that the float-and-sink test, not a simple analytical measurement but a process test, comprehends a number of more-or-less independent items of fundamental data; and as shown in the report, the' variability of the samples differs with respect to these different items. This condition and the wide variety of situations in which float-and-sink test data are used, in combination with other factors, to study complex process operations, indicate the difficulty of setting up fixed standards for float-and-sink sampling and testing. At this stage at least, it is the intent rather to obtain experimental data on the principal Factors involved so that the reader may arrive more intelligently at a procedure adapted to his problem. The authors of this paper have presented experimental data showing the relationship between size of sample and particle size for different variability tolerances with respect to percentage of sink. In further studies, data are being collected to appraise also the variability of the samples with respect to float-ash content and size consist. The scope of this work will be broadened to cover particle sizes up to 4 in., and a third series of tests has yielded similar data for a much cleaner type of raw coal. Thus, the further studies will make available a fairly comprehensive meas- ure of variability with respect to size consist, percentage of sink, percentage of middlings, and percentage of ash in the float, for three coals ranging from 3.87 to 17.68 pct in refuse content (heavy sink material) and from 4.18 to 17.52 pct in middlings. Introduction At present there are no published standards for float-and-sink test sampling. During the rapid expansion of float-and-sink test work, varying procedures have been based on adaptations of the ASTM standards for sampling coal for analyses. For this reason, a special study of gross sample reduction has been undertaken to determine the limits for this step in the operation where no reduction in particle size is to be made before testing. Float-and-sink tests are made whenever a thorough study of coal characteristics is desired. The tests may be made on samples from coal-cleaning units such as jigs or tables, or coal samples may be tested which are taken from a loading boom, railroad car, or the coal seam itself. The resultant gross sample may be large and pose a problem of sample reduction. The question is, then, how much can the sample be reduced and still fall within preassigned limits of accuracy of the original gross sample of coal? Coal from channel samples may be crushed to liberate impurities and then separated into various gravity fractions from which washability curves are drawn; from these curves, it is possible to determine the cleaning characteristics of the coal. However, coal samples from coal-cleaning units cannot be

Jan 1, 1950

By P. Majani, F. P. Kokesh, M. Grosmangin

Determination of the quality of cementation of casing in oil wells in the past has involved inflow and circulation tests to insure that the producing zones were adequately sealed off from the adjacent zones. Existing logging methods, such as temperature and radioactivity surveys, may detect the presence of cement behind the casing. However, the qualities of the cement (i.e., its hardness and particularly its bond to the casing) are not indicated. The new logging method described in this paper operates on the principle that the attenuation of a sonic pulse transmitted by a casing is greatly increased when that casing is bonded to an outer annulus of hard material (such as set cement) which has an appreciably smaller sonic-wave velocity than that of the casing. The down-hole tool contains a source of recurrent sound pulses which are detected by a receiver spaced a few feet from the source The amplitude of the detected casing-borne pulse is measured, and the resulting signal is transmitted to the surface where it is recorded vs depth. Since amplitude is a function of attenuation, the log is readily interpreted. Laboratory studies have shown straightforward relationship between attenuation and such variables as source-detector spacing and per cent of circumference bonded. It is shown that cement not set or not bonded to the casing has compara- tively little attenuating effect. Field examples show not only the cement top, but also the variation in cementation quality below the top. Further, the increase of bonding with time and after squeeze cementation is depicted. The detection of poor cement jobs is confirmed by production tests and by formation-test results. It is anticipated that the method will have wide application in evaluating cementation quality prior to formation testing in completions and recompletions. The analysis it affords may aid in further improving cementation techniques. INTRODUCTION The main purpose of oilwell cementing is to isolate a production zone from other undesirable zones. To investigate whether this purpose has been accomplished, several logging methods have been used such as temperature logs, radioactive tracer logs, etc. While these logs all respond to the presence of cement behind the casing, they do not indicate the degree of bonding of the cement to the casing. Early in the application of sonic logging, it was noticed that considerable attenuation of sound signals takes place in cemented pipe and is often made evident on the standard Sonic log by cycle skipping.' The development of a circuit capable of continuously recording the amplitude of the casing-borne sound signal has made possible an extensive series of laboratory and field tests, which gave the following results. The amplitude of a sound signal after it has traveled in a firmly cemented pipe is only a small fraction of that recorded by the same device in free pipe. This provides a wide spectrum of energy levels; for given local conditions, empirical values of the amplitude can be correlated with the quality of cementation. Interpretations made in this manner generally have been confirmed by production tests, circulation tests and squeeze cementation. The purposes of this paper are to give a general description of the new logging method and to present some laboratory and field results. CYCLE SKIPPING Early attempts to study the quality of the cement behind the casing were performed with a standard Sonic log, which measures the transit time At, and were based on the well known phenomenon of cycle skipping.' Cycle skipping normally is interpreted as being a manifestation of weak signals at the receivers. The log of Fig. 1 was run with the recording instrumentation adjusted to enhance cycle skipping. A mirror-image presentation was used for better visual interpretation. The transit time of sound in steel is about 58 microseconds/ft (corresponding to a velocity of 17,000 ft/ sec), and the portions of the log where this value is recorded are interpreted as zones with no cement bond. Where cycle skipping produces a higher value of At, weak signals (or a high rate of attenuation) are indicated, and a good cement bond may be present. The occurrence of cycle skipping, however, depends too much on instrument adjustment to give a uni-

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.