Search Documents

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the

By Spencer S. Shannon

The lanthanide elements from lanthanum (atomic number 57) to lutetium (71) plus yttrium (39) are called the rare-earth elements; scandium (21) is chemically similar to yttrium and the lanthanides, but differs in occurrence. The rare-earth metals were divided into three categories by the chemist, Berzelius, in the early 1800s; however, modern usage defines only two groups. The cerium group is named after the most abundant -element and comprises lanthanum (57), cerium (58), praseodymium (59), neodymium (60), promethium (61), samarium (62), europium (63), and gadolinium (64). Traces of natural promethium were discovered in apatite in Finland (Erametsa, 1965). The yttrium group, also named after its most abundant representative, comprises terbium (65), dysprosium (66), holmium (67), erbium (68), thulium (69), ytterbium (70), and lutetium (71). Yttrium most closely resembles holmium in chemical properties, and occurs most commonly with the heavy lanthanides. The rare-earth minerals may be divided into cerium-group, yttrium-group, or mixed-group classes. A number of rare-earth minerals contain thorium and uranium, but these are variable in content and are not essential to the composition of the minerals. In monazite, a mineral that is a principal source of rare earths, the worldwide average content of thorium oxide is 7.2 % ; however, varieties of monazite contain as little as 0.001% and as much as 31.5% thoria. Monazite is the chief ore of thorium and, until recently, thorium has been the main product; the rare-earth metals were in fact byproducts. However, because of the great demand for rare-earth metals by the petroleum refining and the iron and steel industries since the early 1970s, the roles have been reversed so that thorium is now the byproduct. Uses A 9 to 1 mixture of yttrium and thorium oxides is used for high-intensity incandescent lamps and for high-temperature-resistant windows and lenses. Yttrium oxysulfide and orthovanadate, doped with europium oxide, are used in color television tube phosphors. Yttrium, as well as cerium, praseodymium, and samarium, is alloyed with cobalt to form strong permanent magnets. It hardens magnesium alloys and retards the oxidation of iron-chromium alloys. Yttrium-iron garnets are used for transmitting low microwave frequencies. Yttrium-aluminum garnets are used as host crystals for lasers and as artificial gems. Yttria also forms strong, stable, high-temperature refractories. Gadolinium-gallium garnets are used in substrate wafers on which bubble-memory films are deposited. Rare-earth chlorides, particularly of lanthanum, neodymium, and praseodymium, are used as catalysts for petroleum cracking. Mischmetal is used for lighter flints and to improve the ductility and impact strength of iron and steel alloys. Mischmetal also reduces the fatigue resistance of magnesium and aluminum alloys. Rare-earth oxides are used for polishing glasses and lenses. Individual or specific combinations of high-purity oxides are used as additives in glass to maintain or change its color, dispersion, or refractive index. High-purity oxide combinations are also used to color tile and to opacify porcelain enamel. Rare-earth fluorides and oxides are used in carbon-arc electrodes to provide incandescent white light. Lanthanum or gadolinium oxysulfide phosphors, activated with terbium, are used as x-ray screen intensifiers. Hydrogen can be stored by sorption on lanthanum-nickel alloys (Anon., 1979). Thorium232 is bombarded with slow neutrons in breeder reactors to create fissionable thorium233 for nuclear fuel cells. The 330-MW, high-temperature gas-cooled reactor at Fort St.

Jan 1, 1983

By R. A. Morse, W. W. Owens, D. R. Parrish, T. M. Geffen

Laboratory studies of several factors affecting measurements of relative permeability were made using the three-section plastic-covered core technique. Results show that the core assembly, properly constructed, will perform as a single unit, and that the testing technique will, under suitable conditions of pressure gradient, gas expansion, and migration of partial water saturation, permit measurement of flow characteristics not affected by technique. Wettability equilibrium is readily established in cores exhibiting strong wetting preference to water or oil when initially saturated with water. Laboratory tests must be conducted so that saturation changes represent those that occur in the reservoir. Immediate implications of saturation history are (1:) that the possibility exists of increasing the displacement efficiency of solution gas drive reservoirs over the natural process, and (2) residual gas saturations following water flooding in gas or gas condensate reservoirs will be 15 to 50 per cent pore space rather than 1 to 11 per cent as generally believed. INTRODUCTION Solutions of petroleum reservoir problems pertaining to production performance require the use of true relative permeability characteristics. 'This relationship of fluid conductivity and saturation has been obtained by reservoir engineers in four ways, namely: 1. From past gross reservoir performance and the extrapolation of this data based on experience,1,2 2. By using published fluid flow relationships obtained in laboratory studies on general type porous materials,3, 4, 5 3. By attempting a mathematical derivation of flow behavior," using some experimentally obtained characteristics of reservoir rocks, and 4. By laboratory flow tests using representative rock samples of a reservoir. The first three methods listed above have shortcomings which make their use limited or questionable. Production characteristics of only certain processes are obtained from field data and these are not available at the beginning of a reservoir's producing life, at which time they are desirable. It is fortuitous if general fluid flow characteristics obtained experimentally have accurate application to specific field problems. Also, it is felt that at this time there is not sufficient knowledge of the flow behavior of oil, water, and gas in porous materials to enable applicable analytical description of this to be made based on other measured rock characteristics. Measurement of relative permeability in the laboratory offers the only direct method subject to adequate checking for determination of flow characteristics applicable to field problems. Primarily, this paper deals with laboratory experiments to establish the effects of several factors on the measurement of relative permeability and the practical significance of this knowledge. Laboratory flow tests are essentially model studies; the gross reservoir being represented by a small rock sample. It is necessary to simulate in experiments all the factors of importance which influence flow in the reservoir, as well a.; to eliminate effects peculiar to the model studied to assure measurement of valid flow behavior. Factors which affect the measurement and behavior of fluid flow in porous materials have been studied by other investigators.3,7,8,9,10,11 This analysis is not intended to evaluate fully all factors involved, but to establish first the validity of measurements by the technique used and then study flow behavior under various conditions. The studies reported are as follows: 1. Effects of pressure gradient, 2. Effects of gas expansion in tests conducted at near atmospheric pressure, 3. Effects of migration of a partial water saturation in oil-gas-water systems, 4. Effects If wettability characteristics of the system. and 5. Effects of the fluid saturation history. The first three are important only in laboratory measure-

Jan 1, 1951

By R. H. Hopkins, R. W. Kraft

X-ray and metallographic analysis reveal that the preferred cryslallograPhic relationships in direc-tionally solid~fied Pb-Sn eutectic specimens can be stated: interface growth direction, Experiments in which tin single-crystal seeds were used to nucleate the eutectic demonstrate that the normal lamellar structure observed at steady state can be forced to break down initially to a semide-generate structure when the crystallograPhic direction of the tin seed does not lie within the preferred eutectic interface plane. A normal structure eventually develops by preferred grain growth. X-ray topo-graphs and rocking curves substantiate that the perfection of eutectic samples increases as growth proceeds until steady-state conditions prevail. RECENT years have seen substantial progress made in the understanding of eutectic solidification through both intensive theoretical1'' and experimental investigations (see Ref. 7 for a review and extensive bibliography of both aspects). Several important observations concerning the nucleation and growth aspects of binary eutectic solidification have been brought to light. Briefly some of these include the nonrecipro-cal nucleation of one phase upon another,'j9 the tendency for a specific low-index direction of one or both phases of the eutectic to be aligned in the growth direction during unidirectional solidification,'0711 and the production of lamellar microstructures having a preferred interface plane between "11-13 phases. Questions quite naturally arise as to the mechanism by which preferred grain growth evolves during eutectic solidification and what relationships (if any) connect the initial nucleation of one phase upon another with the preferred growth directions and interface planes observed after a few centimeters of growth in many systems. This study was initiated in an attempt to answer these questions. An experiment was devised in which both the nucleation and growth steps in the crystallization of a binary eutectic, Pb-Sn, could be correlated with subsequently developed eutectic morphology and crystallography. This was accomplished by nucleating the eutectic upon a tin single crystal seed and then causing the eutectic to grow from the seed by directional freezing. The morphology of the structures produced in this manner was then monitored metallographically and the orientation relationships and substructural perfection determined by X-ray analysis. The information thus ob- tained was correlated to the growth history of the samples. EXPERIMENTAL DETAILS The eutectic of composition 38.1 wt pct Pb and 61.9 wt pct sn14 was prepared by accurately weighing lead and tin of 99.999+ (from ASARCO and Lytess Metals and Chemicals, respectively), melting in a clean Pyrex beaker under an argon blanket, and casting into the proper crucible for resolidification as described below. Tin seeds were grown by the method of ~halmers," saw-cut, oriented by the Laue method, and soldered to the pull rod of the Czochralski puller. All seeds were given a final electropolishing (to assure the removal of any deformed material from the cutting operation) and acetone degreasing immediately prior to the nucleation experiment in which they were to be used. Eutectic samples were unidirectionally solidified in a horizontal graphite boat1' or by pulling from the melt16 with tin (or graphite) seeds. A flowing argon atmosphere was employed in all experiments and growth rates and temperature gradients, constant for a given run, varied from 2 to 4 cm per hr and 8" to 15°C per cm over the course of all the experi-

Jan 1, 1969

By F. W. Warnaars

The Los Bronces copper deposit is located on the west side of the Andes Mountains in central Chile about 54 km (34 miles) northeast of Santiago. The deposit consists of a hydrothermal breccia complex that was emplaced in the western part of an earlier porphyry copper mineral system. This system exhibits propylitic, sericitic, silicic to potassic hydrothermal alteration with disseminated and stockwork copper-molybdenum sulfide mineralization over an area of about 12 km2 (5 sq. miles). The Rio Blanco Mine, about 800 m (0.5 miles) east of Los Bronces, is an underground block caving mine located in the central part of this system. The mineralization at Rio Blanco is characterized by one of the highest primary copper grades known for porphyry copper systems. Present mining at Los Bronces is from an open pit located in the northern part of a hydrothermal complex that is composed of at least seven different copper-bearing tourmaline breccia pipes. Together these pipe form one large contiguous body about 2 km (1.2 miles) long and 0.7 km (0.4 miles) wide at the surface. This body crops out at elevations of between 3.5 and 4.1 km (2 and 2.5 miles) above sea level. Each breccia pipe has a characteristic shape, matrix, type and degree of mineralization and alteration. Quartz, tourmaline, specularite, anhydrite, pyrite, chalcopyrite, sericite, chlorite, and rock flour are present in significant quantities in most breccia matrices. The clasts consist mainly of quartz monzonite, andesite, and minor amounts of quartz latite porphyry. The seven different breccia pipes are identified, from oldest to youngest, as Ghost, Central, Western, Infiernillo, Fine Gray, Anhydrite, and Donoso. The breccias coalesce, interfinger, or display gradational contacts. In places they form sharp contacts with each other and with the surrounding intrusive or volcanic host rocks. Early breccias are disrupted or cut by subsequent breccia events. Explosive emplacement was followed by collapse after the release of hydrothermal pressures. The Los Bronces breccia complex has sharp contacts with the surrounding intrusive stock and andesites. These contacts dip inward on the north, west, and south margins. The eastern contact is nearly vertical and/or dips slightly outward and suggests a westward tilting of the body of about 15° after emplacement. The breccia complex is not fully explored at depth, but is known to extend at least below the 3 km (2 mile) elevation in its central eastern part, about 1 km (0.6 miles) below the highest breccia outcrop. Primary copper sulfide minerals are most abundant in the Donoso and Infiernillo breccias where the matrix consists mostly of tourmaline, quartz, sulfides, specularite, and/or chlorite with minor amounts of rock flour, sericite, and clay. The volume of the matrix material is mostly smaller than in the other breccias. Secondary enrichment enhanced the primary copper grade in the southern two-thirds of the Los Bronces breccia complex and in much of the surrounding porphyry system. Primary grades of 0.3% to 0.5% copper have been enriched to approximately 0.9%. This enrichment may extend from surface outcrops to a depth of more than 500 m (1, 640 ft) in more permeable parts of the breccia complex. A leached and oxidized capping of variable thickness covers a large part of the enriched layer and exceeds 300 m (984 ft) at one location. The shape and depth of the enrichment blanket and the capping suggest that the enrichment process is still active.

Jan 1, 1983

By R. J. Stokes

This paper compares the room-temperature mechanical behavior of magnesium oxide crystals containing 'fresh' and 'grown in' sources. 'Fresh' dislocation sources introduced by cleavage or mechanical contact move and multiply to generate slip bands at 3000 psi. Macroscopically the crystals yield smoothly at 6000 psi and can be quite ductile. When 'fresh' sources are eliminated by chemical polishing the crystals have to rely on 'grown in' sources to initiate slip. They then deform elastically up to extremely high stresses (-30 to 50,000 psi) yielding with a sharp stress drop down to the level at which Ifresh' sources operate. Grown in' sources are considered to be associated with impurity wecipitate particles or inclusions. Annealing crystals at 2000°C dissolves these particles thereby eliminating sources for slip. The crystals then become consistently stronger (> 70,000 psi) approaching within 1/10 of the theoretical strength. Dislocations of the grown-in network do not move under stress to nucleate slip. It is proposed that they are locked electrostatically. DISLOCATIONS were originally conceived to resolve the large discrepency between the theoretical and real strength of solids. The idea that slip propagates by the movement of dislocations is now well confirmed both by electron transmission and etch pit studies. But this by itself is not sufficient to explain the low yield strength of crystals for there is still the need to understand where slip dislocations originally come from. Theoretical estimates place the stress to generate dislocations homogeneously in the lattice at approximately 1/30 of the shear modulus,' a value 104 to 105 times greater than the stress at which the very first manifestations of plastic deformation have been observed. Only an special cases can the theoretical strength be approached; in crystals of high perfection such as whiskers2 where dislocation sources do not preexist and in covalent crystals where the dislocations present at room temperature are almost completely immobilized by the bond character of the solid.3,4 For the most part it is necessary to assume that the yield strength of ionic and metallic solids corresponds to the motivation and multiplication of dislocation sources present beforehand. The question as to their origin appeared to be answered by the suggestion of Frank and Read5 that dislocations forming part of the grown-in structure could operate as a mill for the large scale production of dislocations. The Frank-Read mechanism is still widely accepted as the source of dislocations but a number of observations have been reported which are in serious disagreement with it. One of the most intriguing of these is that under normal testing conditions dislocations of the grown-in network appear to play no role in the plastic deformation of rock salt structure ionic solids. This was first reported in the etch pit studies on lithium fluoride by Gilman and Johnston6 and is true also for magnesium oxide. To explain the comparatively low yield stress of lithium fluoride it was proposed by a number of authors7,9 that small unresolved prismatic loops formed by vacancy condensation were acting either directly or indirectly as the dislocation sources. This proposal led Gilman10 to perform a series of detailed and careful experiments on lithium fluoride to try to identify the sources responsible for the initiation of slip. He concluded that prismatic loops played no general role and that other sources were more important. Dislocation sources in ionic crystals can be divided into two general categories. First, there

Jan 1, 1962

By D. R. Parrish, T. M. Geffen, R. A. Morse, W. W. Owens

Laboratory studies of several factors affecting measurements of relative permeability were made using the three-section plastic-covered core technique. Results show that the core assembly, properly constructed, will perform as a single unit, and that the testing technique will, under suitable conditions of pressure gradient, gas expansion, and migration of partial water saturation, permit measurement of flow characteristics not affected by technique. Wettability equilibrium is readily established in cores exhibiting strong wetting preference to water or oil when initially saturated with water. Laboratory tests must be conducted so that saturation changes represent those that occur in the reservoir. Immediate implications of saturation history are (1:) that the possibility exists of increasing the displacement efficiency of solution gas drive reservoirs over the natural process, and (2) residual gas saturations following water flooding in gas or gas condensate reservoirs will be 15 to 50 per cent pore space rather than 1 to 11 per cent as generally believed. INTRODUCTION Solutions of petroleum reservoir problems pertaining to production performance require the use of true relative permeability characteristics. 'This relationship of fluid conductivity and saturation has been obtained by reservoir engineers in four ways, namely: 1. From past gross reservoir performance and the extrapolation of this data based on experience,1,2 2. By using published fluid flow relationships obtained in laboratory studies on general type porous materials,3, 4, 5 3. By attempting a mathematical derivation of flow behavior," using some experimentally obtained characteristics of reservoir rocks, and 4. By laboratory flow tests using representative rock samples of a reservoir. The first three methods listed above have shortcomings which make their use limited or questionable. Production characteristics of only certain processes are obtained from field data and these are not available at the beginning of a reservoir's producing life, at which time they are desirable. It is fortuitous if general fluid flow characteristics obtained experimentally have accurate application to specific field problems. Also, it is felt that at this time there is not sufficient knowledge of the flow behavior of oil, water, and gas in porous materials to enable applicable analytical description of this to be made based on other measured rock characteristics. Measurement of relative permeability in the laboratory offers the only direct method subject to adequate checking for determination of flow characteristics applicable to field problems. Primarily, this paper deals with laboratory experiments to establish the effects of several factors on the measurement of relative permeability and the practical significance of this knowledge. Laboratory flow tests are essentially model studies; the gross reservoir being represented by a small rock sample. It is necessary to simulate in experiments all the factors of importance which influence flow in the reservoir, as well a.; to eliminate effects peculiar to the model studied to assure measurement of valid flow behavior. Factors which affect the measurement and behavior of fluid flow in porous materials have been studied by other investigators.3,7,8,9,10,11 This analysis is not intended to evaluate fully all factors involved, but to establish first the validity of measurements by the technique used and then study flow behavior under various conditions. The studies reported are as follows: 1. Effects of pressure gradient, 2. Effects of gas expansion in tests conducted at near atmospheric pressure, 3. Effects of migration of a partial water saturation in oil-gas-water systems, 4. Effects If wettability characteristics of the system. and 5. Effects of the fluid saturation history. The first three are important only in laboratory measure-

Jan 1, 1951

By N. J. Grant, C. Stein

A more complete determination of the chromium-rich portion of the Cr-Ni phase diagram was made, based in part on resistivity measurements. The existence of a eutectoid reaction at 1215°C,made,based in addition to the previously reported eutectic reaction at 1343°C, was confirmed. Metallographic studies, lineal analysis, and diffusion-couple evidence were utilized to support Metallographicand extend the findings of the resistivity measurements. Delineation of the a chromium solvus line shows that the maximum solubility of nickel in a chromium decreases rapidly from 25 wt pct Ni at the eutectoid temperature to 16 wt pct. Ni at 950°C. IN recent years alloys of high chromium and nickel content have become of more interest as better performance of metals at higher temperatures is sought. It has become increasingly important, therefore, that the Cr-Ni binary diagram, which is composed of two of the elements commonly found in many high temperature alloys, be known more accurately. In particular, the chromium-rich portion of this diagram has for a long time remained incompletely determined and has been reported in the literature with a great deal of conflicting evidence. In general, electrolytic chromium of nominal purity was used by most of the earlier workers, who, except for Jenkins et al., made no further attempt to purify the metal. In the present redetermination of the chromium-rich portion of the Cr-Ni phase diagram, every effort was made to use the purest component metals available and to retain this purity throughout the preparation of the alloys and their subsequent handling. Early work by Voss, Bain," and Phebus and Blake4 indicated that the Ni-Cr binary diagram was a simple eutectic system with rather extensive terminal solid solutions. Later, Nishigori and Hamasumi," Tasaki, Weaver and Jellinghaus,' Matsunaga, and Jette, Nordstrom, Queneau, and Foote attempted to determine the limits of the phase fields. Except for Jette and coworkers, the other investigators agreed that the temperature of the eutectic reaction occurred between 1320" and 1340°C and at a composition of 50 or 51 pct Cr. There was, however, considerable disagreement as to the melting point of the chromium-rich alloys and the limits of the solid-solution fields. As a result of the studies of Jenkins et al., the liquidus and solidus curves up to 80 and 90 pct Cr, respectively, were outlined. In addition, the solid-solution phase fields at higher temperatures had been investigated and showed only two phases. The boundary of the chromium-rich a phase was not determined precisely according to these authors, but indications were that it tended to confirm the earlier work of Nishigori and Hamasumi. A more recent investigation by Taylor and Floyd,'" using X-ray methods, essentially corroborated the diagrams proposed by Jenkins et al., and by Wise and Eash (see Fig. 1). Again, no indication of a third phase in the chromium-rich portion of the diagram appeared. The earlier observations, however, are not substantiated by the modifications to the Ni-Cr binary diagram proposed by Bloom and Grant." In the investigation carried out by Bloom and Grant, X-ray methods showed that an eutectoid reaction existed at approximately 1180°C, in addition to the eutectic reaction at 1343°C (see Fig. 1). Their determination of the melting point of pure chromium was 1930°C (later revised to 1903"C)." Experimental Procedure and Results The chromium used for this study was electrolytic chromium, made available by the United States

Jan 1, 1956

By J. S. Aronofsky, A. S. Lee

Results and experimental procedures are presented covering a preliminary laboratory investigation of the compaction of reservoir rocks and its effect on porosity and permeability. "Egective" compressibility was measured by subjecting core material to an overburden load and reducing pore pressure. Preliminary results showed relatively high effective compressibilities which varied from 1.8 X 10-' pore volume/pore volume /lb / in.~ for limestones to 7.7 X 10-6 pore volume/pore volume/lb/in.' and even higher for sandstone reservoirs conmining large amounts of fine clay material. It appeared from the initiol results that the effective compressibility of the material could not be correlated directly with porosity and that other factors, such as the amount of clay material present might have to be taken into account. The work to date has indicated that effective compressibility factors will have to be measured for the particular reservoirs heing studied. Oil permeability of core materid was measured under various effective overburden loads from 0 to 8,000 psi. It has been found than reduction of permeability in clean sands was relatively small and the results agreed with published data. For sandstones containing large amounts of clay, preliminary tests indicated very large reductions in permeability with increasing effective overburden pressures. It is planned to continue laboratory work on this phase because of its importance in low-permeability reservoirs. INTRODUCTION During the past few years increasing interest has been shown in the differences obtained between conventional laboratory measurements on cores at atmospheric pressure and those when reservoir conditions are simulated. Since these data are used as the basis for calculations and the prediction of reservoir behavior, their evaluation under reservoir conditions is most important. The effect of the compressibility of the rock upon the reservoir behavior, as fluid pressure declines, is an important subject which has been insufficiently studied. The total pressure on any plane through a reservoir is the result of the grain-to-grain rock pressure due to the overburden and the pressure of the interstitial fluid on the plane. As reservoir fluid pressure declines, the net load on the rock phase must increase, and the result is the compaction of the rock structure which produces a decrease in pore volume. Simultaneously, the rock grains expand into the pore space as the fluid pressure declines. The result of these two effects is defined as the effective compressibility of the rock and is expressed in units of pore volume/original pore volume/psi. Work was started by Imperial in 1954 on the effects of pressure on the porosity of reservoir rock. These studies established effective compressibility factors for some limestones and dolomites. Since that time effective compressibilities have been obtained for other Canadian and some U. S. reservoir rocks. In conjunction with the study of porosity, an investigation has also been initiated to determine the effects of pressure on permeability. One outcrop and several reservoir sandstone cores have been tested to date. ROCK COMPRESSIBILITY Apparatus anD Procedure A diagram of the apparatus used in these studies is shown in Fig. 1. The equipment was constructed so that a uniform pressure, equivalent to the reservoir overburden, could be applied by the oil pump and transmitted on all sides of the core through the Lucite sheath. The sample was completely saturated by circulating brine through the mounted core. Pressure in

By D. M. Boulin, E. S. Greiner

Single crystals of vanadium were bent at 293" and 77°K to a maximum fiber strain of about 1 pct. Analysis of slip line and etch markings of- crystals deformed at either temperature show the slip planes are {112}. The slip directions at 298° and 77°K, deterwined from the asterism axes of rotation of crystals compressed in a {110} direclion, were (111). Etch pits, associated with dislocations, were observed in the plastically deformed vanadium. There was no evidence of mechanical twinning. THIS paper presents results on single crystals of vanadium plastically deformed by bending and compressing at 298° and 77°K. Bcc metals have been observed to slip on {110}, (1121, and (123) planes while the slip direction has invariably been found to be ( 111) 1 Lindley and smallman2 have studied the deformation intension of polycrystalline vanadium between 298° and 20°K and reported that at and above 77°K the deformation is predominantly by slip, but at 20°K the deformation is associated with twinning. Later, Edington and smallman3 found that single crystals of vanadium compressed in a (111) direction at 77°K deformed by twinning on { 112}( 111) systems. Clough and pavlovic4 have reported that mechanical twins were observed in vanadium loaded by impact at 77°, 173°, and 195°K, respectively. PREPARATION OF SINGLE CRYSTALS The vanadium single crystals were made from commercial vanadium reported to contain 0.03 pct C, 0.03 pct 0, 0.02 pct N, 0.0X* pct Ca, 0.0X pct Si, proximately 0.8 cm diam, 4 cm long and were grown in (100) and (110) directions. A crystal for studies of etch pits was gound to have a face parallel to a {l00} plane by the Bond crystal-orienting technique.7 In order to avoid straining, the crystals were placed on the steel mount with a low-melting organic compound (glycol phthalate) before grinding. Crystals for the bending experiments were ground to have the four side faces parallel to ( 100) planes. The crystals were electropolished8 in 10 pct H2So4 using a current density of about 3 amp per sq cm to remove the strained surface material. The crystals were bent slowly by three-point loading in a stainless-steel jig designed by Sylwestrowicz.9 The end rollers of the jig were 3.8 cm apart. DEFORMATION AT 298°K The bent vanadium crystal, after removal from the jig, had a radius of curvature of 20.4 cm over the 2-cm length close to the center loading roller, where the maximum bending moment occurred. This radius of curvature corresponded to a maximum fiber strain of 0.9 pet. The side faces contained slip bands as shown in Fig. 1. These bands are straight where the strain is high, top and bottom of micrograph, but are curved on either side of the neutral axis in the center. This curvature of the slip bands in the low-stress area may be due to cross slip. The plastically deformed crystal was subsequently electropolished and then electroetched in 10 pct H2SO4 solution, using a current density of 0.25 amp per sq cm. The structure is shown in Fig. 2. The directions of the traces on two adjoining faces were plotted stereo graphically, Fig. 3.* Since there were two sets of

Jan 1, 1968

By H. S. Price, D. A. T. Donohue

The system of equations describing displacement of a hydrocarbon liquid by a hydrocarbon vapor in a porous medium where mass transfer takes place between the phases is solved numerically for a variety of gas injection processes. Even though the method of solution is quite general, only systems with three hydrocarbon components are considered. Computer simulations of displacement processes wherein mass transfer between phases is both considered and neglected are compared, and it is shown that neglecting mass transfer can give pessimistic displacement efficiencies. INTRODUCTION The role of the gas displacement process in the recovery of petroleum has been subjected to a series of detailed analyses; as a result, a number of predictive models have been published in the literature. However, because of major simplifying assumptions, most of these models do not completely represent the physical system. As a result, the effect of making the simplifying assumptions is unknown. Therefore, a complete representation of this process — one without major simplifying assumptions —should lead to a fuller understanding of the process, and perhaps to methods of improving it. The general method of developing a model for two-phase fluid flow in a porous medium is to solve simultaneously the continuity equation, the energy equations and the equation-of-state for each phase under the prescribed initial and boundary conditions. For an isothermal system, the energy equations reduce to the momentum equation. Darcy's law. However, since natural gas is the vapor state of the reservoir liquid, interphase mass transfer may take place with concomitant changes in both the intensive and extensive thermodynamic properties of each phase. It is this phenomenon that has often been omitted in previous mathematical models. An additional relation, then, which accounts for mass transfer between the phases, must be included with the other equations to specify a complete model. Completely formulating the equations to be solved is not a difficult task1 but obtaining their solution has been intractable up to now. Availability of large-memory, high-speed digital computers now makes an attack on this formidable problem possible. This paper presents a preliminary study of the problem. Since this investigation is intended to be exploratory, it is restricted to the linear, horizontal, isothermal, two-phase viscous flow of oil and gas in an oil reservoir. In the early development of predictive models of this process, the reservoir system was considered as a unit and various forms of the material balance equation were proposed.2-4 Pressure and saturation gradients were then added in the Buckley-Leverett model.5 The Buckley-Leverett formulation considered the fluids to be incompressible; thus, the mathematical model reduces to a steady-state system. In the 1950's, studies incorporating numerical techniques were being published.6-9 These mathematical models differed in the efficiency of finite difference techniques, the inclusion or exclusion of capillarity or the number of space dimensions considered. To solve these nonlinear, partial differential equations, each phase was considered to be homogeneous with time; therefore, mass transfer between phases was neglected. The effect of mass transfer on the gas displacement process was first reported by Attra. 10 He simulated the one-dimension flow system by a series of cells in each of which the fluids were equilibrated during a time step. In addition, the pressure throughout the system during each time step was predetermined and constant-phase velocities were calculated according to the Buckley-Leverett incompressible fluid flow model. Welge et al.11 developed a model for the displacement of oil by an enriched gas where composition is considered to be a dependent variable. Fluid

By Robert T. MacMillan

Of all the mineral substances utilized by man, salt or sodium chloride has one of the longest and most varied histories. Because all animal life is descended from marine organisms, sodium and chlorine became necessary to life, and as sodium chloride meets requirements most readily, common salt has been essential to human and animal diets from earliest times. As civilization progressed, man learned to evaporate sea water and brines of saline lakes and springs to obtain salt, and gradually ceased to depend on the flesh of animals and the occasional salt outcrops or "salt licks" they frequented. He also learned that salt had a multitude of uses, particularly as a preservative of meat and other foods. The importance of salt in the minds of early men may be readily seen from many allusions in the Bible and other ancient writings. Expressions such as "not worth his salt," and "the salt of the earth" are deeply embedded in the English language. Salt has symbolized friendship, fidelity, purity, and hospitality. It has been used as currency; had religious significance; caused wars; and has been a favorite commodity for taxation. The English word "salary" stems from the Latin word for salt, in which Roman soldiers received part of their pay. Although salt today retains its dietary importance, it has thousands of other uses. Of each 100 pounds produced in the United States less than three pounds is for table use. The other 97 pounds goes for agricultural, industrial, and chemical uses. Fortunately, reserves of salt are enormous and reasonably well distributed on the face of the earth. The sea is an inexhaustible source of salt which may be had at prices which vary with climate and other factors. In addition, large underground salt deposits and many salt lakes, springs and playas are scattered throughout the world. Most countries have access to a salt supply, although some are more favored than others. Properties Crystallized sodium chloride or common salt is known as halite. Its molecular weight is 58.454 and it contains, by weight, 39.34 pct sodium and 60.66 pct chlorine. Alternate positions in the cubic lattice are occupied, respectively, by negative chloride and positive sodium ions; the crystal form is isometric. Crystallization from water produces concentric shells, and when broken, a crystal shows conchoidal fracture. Salt has a hardness of 2.5; its specific gravity ranges from 2.1 to 2.6 for rock salt and is 2.165 for the pure crystal. The refractive index is 1.554; it melts at 800.8 ± 0.5°C and boils at 1413°C. Its heat of fusion is 123.59 cal/ gram. Salt is transparent or translucent and highly diathermous. It is colorless when pure and has a vitreous luster. Salt containing traces of iron oxides has varying shades of yellow, orange, brown, pink, or red. The peculiar shade of red sometimes found in that produced by solar evaporation of sea water has been attributed to tiny organisms. Blue mottlings sometimes seen in halite crystals have been attributed to colloidal sodium. Salt is soluble in water to the extent of 35.7 parts per 100 parts of water at 0°C and 39.8 parts per 100 of water at 100°C. The small difference in solubility with increasing temperature is important in the economics of producing salt by evaporating brine. Salt has a negative heat of solution which is 1.281 kg-cal/ mole at 18°C in 200 moles of water.

Jan 1, 1960

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge's method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey's theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey's method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge's6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

By L. S. Darken

A quadratic formalism is developed lor the representation of the excess free energy, and of the activity coefficients of each component of a ternary system in the vicinity of a single component selected as solvent. This formalism, in contrast to the e formalism, is thermodynamically consistent at .finite concentrations. Numerous specific systems are used as illustrations, with cross checks where possible; these systems are Fe-C-X systems in the vicinity of 1600°C and include practically all those for which reliable data are available. It is shown that this formalism is usually quite adequate to rebresent available data up to a solute concentration of 20 to 30 at. pet. It is hoped that the formalism here presented may prove useful for the treatment and tabulation of data on ternary and multi-component solutions. SINCE, as discussed in a prior paper,' we do not as yet have an adequate basic understanding of the nature of binary metallic solutions, we cannot hope to be in better position for ternary and multicomponent systems. However, the accumulation of thermodynamic data on ternary and multicomponent metallic solutions has led to an urgent need for a rational formalism for the presentation and use of such data. Notable attempts to represent the isothermal isobaric thermodynamics of ternary solutions were made by Benedict et 01.' and by Wohl.3 These attempts were directed toward developing an analytic expression adequate to encompass the entire ternary triangle, primarily for organic systems; they generally involve the power-series approach, which has been shown' to be inappropriate, in general, even for binary metallic systems. On the other hand, interest in the thermodynamics of ternary and multicomponent metallic solutions has tended to focus either on general thermodynamic relations4-10 over wide compositional ranges, or on cases where one or more components have been at rather low concentration. This latter approach has been followed by Wagner,= by chipman,11 and by Alcock and Richardson,12 and more recently in many publications by the same and other authors. It has already been shown' that the thermodynamic behavior of binary metallic systems is, in general, relatively simple only in the terminal regions. Hence, it would appear that any approach, at this time, to a simple formalism for ternary and multicomponent systems must, a fortiori, be limited to the vicinity of the major component or solvent. The present treatment is aimed to be useful for a given solvent with moderate solute concentrations ranging up to 20 or 25 at. pet and in exceptional cases to considerably higher concentrations. In order to achieve even this modest objective, it is necessary that any appropriate formalism have built into it 1) a thermodynamic consistency (at least within the desired degree of approximation), and 2) an approach at infinite dilution to Raoult's law for the solvent and to Henry's law for each solute. Let us first formulate the condition of thermodynamic consistency for a ternary system. Under isothermal isobaric conditions, we may write for the molal value G of any extensive function dG - C1 dN1 + G2 dN2 + G3 dN3 [l] where the G's are the corresponding partial molal quantities and the N's are the atom fractions. Taking component 1 as the solvent, this relation may be rewritten in terms of the independent variables N2 and N3

Jan 1, 1968

By C. S. Land

Relative permeability functions are developed for both two- and three-phase systems with the saturation changes in the imbibition direction. An empirical relation between residual nonwetting-phase saturation after water imbibition and initial nonwetting-phase saturations is found from published data. From this empirical relation, expressions are obtained for trapped and mobile nonwetting-phase saturations which are used in connection with established theory relating relative permeability to pore-size distribution. The resulting equations yield relative permeability as a function of saturation having characteristics believed to be representative of real systems. The relative permeability of water - wet rocks for both two- and three - phase systems, with the saturation change in the imbibition direction, may be obtained by this method after properly selecting two rock properties: the residual nonwetting-phase saturation after the complete imbibition cycle, and the capillary pressure curve. INTRODUCTION Relative permeability is a function of saturation history as well as of saturation. This fact was first pointed out for two-phase flow by Geffen et al. 1 and by Osaba et a1.2 Hysteresis in the relative permeability-saturation relation also has been reported for three-phase systems.3 Since saturations may change simultaneously in two directions in a three-phase system, four possible relationships arise between relative permeability and saturation for a water-wet system. The four saturation histories of this system were given by snel14 as 11, ID, DI and DD. I and D refer to the direction of saturation change (imbibition and drainage), with the first letter of the symbol indicating the direction of change of the water phase. As used in this paper, the second letter of the symbol refers to the direction of saturation change of the gas phase, i.e., D and I indicate an increase and decrease, respectively, in gas saturation. Only a few three-phase relative permeability curves have been published. Leverett and Lewis5 published three-phase curves for unconsolidated sand, and Snell4 reported results of several English authors for both drairrage and imbibition three-phase relative permeability of unconsolidated sands. Three-phase relative permeability curves for a consolidated sand were published by Caudle et a1.3 for increasing water and gas saturations (ID). Corey et a1.6 reported drainage (DD) three-phase relative permeability for consolidated sands. Recently, Donaldson and Dean7 and Sarem8 calculated three-phase relative permeability curves from displacement data on consolidated sands, also for saturation changes in the drainage direction. The only published three-phase relative permeability curves for consolidated sands with saturation changes in the imbibition direction (11) are those of Naar and Wygal.9 These curves are based on a theoretical study of the model of Wyllie and Gardner 10 as modified by Naar and Henderson.11 Interest in three-phase relative permeability has increased recently due to the introduction of new recovery methods and refinements in calculation procedures brought about by the use of large-scale digital computers. The scarcity of empirical relations for three-phase flow, and the experimental difficulty encountered in obtaining such data, have made the theoretical approach to this problem attractive. RELATIVE PERMEABILITY AS A FUNCTION OF PORE-SIZE DISTRIBUTION Purcell used pore sizes obtained from mercury-injection capillary pressure data to calculate the permeability of porous solids.12 Burdine extended the theory by developing a relative permeability-pore size distribution relation containing the correct tortuosity term.13 The works of Purcell and Burdine were combined by Corey14 into a form that has

Jan 1, 1969

By V. O. Struebing, K. A. Johnson, F. H. Ellinger, C. C. Land

The system Pu-Ce-(Pu,Ce)Co2 has been investigated by thermal, micrographic, and X-ray diffraction methods. The principal characteristics of this system were found to be: 1) a "eutectic" valley extending across the phase diagram from the 6 Pu/Pu3 Co eutectic at 408°C to the 7 Ce/Ce24Co11 eutectic at 42Z°C, the maximum temperature in the valley being 443°C; 2) three ternary inter mediate phases, each having an appreciable mutual solid solubility for plutonium and cerium; 3) eleven four-phase equilibria comprising one eutectic, four peritectics, and six Class II reactions; 4) a pseudobinary section between the ternary phase A and the y Ce solid solution, and another such section between the 6 Pu solid solution and the 77 phase (Pu3 CO) containing cerium in solid solution. ALLOYS of plutonium that are completely molten at temperatures in the vicinity of 500°C are of interest for use in the liquid state as fast reactor fuels. Since exploratory work1,2 on Pu-Ce-Co alloys had indicated that a "eutectic" valley extended across the ternary system between the Pu-co[6 Pu/n (Pu3co)] and the Ce-CO[? Ce/? (ce24co11)] binary eutectics, thus allowing alloys of any plutonium concentration to be selected within the temperature limitation, the present investigation was begun for the purpose of accurately defining the "eutectic" valley. However, when the existence of three ternary intermediate phases and the consequent complexity of the liquid-solid equilibria were discovered, the investigation was expanded to cover the complete system, except for the corner between the isomorphous series PuCo2-CeCo2 and pure cobalt. 1) EXPERIMENTAL The purities of the alloy components are given in Table I. Alloy buttons, each weighing about 5 g, were prepared by arc melting the accurately weighed components (±0.1 mg) under a purified He-Ar atmosphere. Weight loss as a result of melting was usually less than 1 part in 1000; therefore, the buttons were not chemically analyzed. One half of each button was re-melted by induction heating in vacuo in a magnesia crucible and furnace-cooled, so as to produce a somewhat coarser microstructure than that of the rapidly cooled arc-melted button. About 20 g of each alloy to be examined by thermal analysis were remelted in vacuo in a tantalum crucible of the inverted thermocouple well type. The thermal, X-ray diffraction, heat treating, and micrographic methods used in this investigation were essentially the same as have been described in earlier reports.3,4 Two etchants for these alloys were found to be satisfactory: A) 8 (parts by volume) orthophosphoric acid, 5 ethanol, and 5 glycol; B) 7 tetraphosphoric acid, 36 water, and 57 2-ethoxyethanol. Etchant A was used at 5 to 6 v for 15 to 20 sec, and etchant B was used at 10 v for 10 sec. 2) THE LIMITING BINARY SYSTEMS 2.1) Plutonium-Cerium. The phase relationships of this system have been investigated by Ellinger, Land, and Cramer,5 and by Selle and Etter.6 The over-all agreement of their phase diagrams is good, but there are some discrepancies in the positions of some phase

Jan 1, 1967

By Haydn H. Murray, Sam H. Patterson

The term clay is somewhat ambiguous unless specifically defined, because it is used in three ways: (1) as a diverse group of fine-grained minerals, (2) as a rock term, and (3) as a particle-size term. Actually, most persons using the term clay realize that it has several meanings, and in most instances they define it. As a rock term, clay is difficult to define because of the wide variety of materials that compose it; therefore, the definition must be general. Clay is a natural earthy, fine-grained material composed largely of a group of crystalline minerals known as the clay minerals. These minerals are hydrous silicates composed mainly of silica, alumina, and water. Several of these minerals also contain appreciable quantities of iron, alkalies, and alkaline earths. Many definitions state that a clay is plastic when wet. Most clay materials do have this property, but some clays are not plastic; for example, halloysite and flint clay. As a particle-size term, clay is used for the category that includes the smallest particles. The maximum-size particles in the clay-size grade are defined differently on various grade scales. Soil investigators and mineralogists generally use 2 micrometers as the maximum size, whereas the widely used scale by Wentworth (1922) defines clay as material finer than approximately 4 micrometers. Some authorities find it convenient to use the term clay for any fine-grained, natural, earthy, argillaceous material (Grim, 1968). When used this way, the term includes clay, shale, or argillite, and some argillaceous soils. Even though no standard definition of the term clay is accepted by geologists, agronomists, engineers, and others, the term is generally understood by those who use it. Clay is an abundant natural raw material, and it has an amazing variety of uses and properties which will be discussed in this chapter. Clays rank among the leading industrial minerals in both tonnages produced and total value. Their importance is indicated in Table 1, which lists the 1978, 1979, and 1980 production and values according to the US Bureau of Mines (Ampian, 1980). Statistical data summarizing world production of kaolin, bentonite, and fuller's earth are compiled by the US Bureau of Mines (Ampian and Polk, 1981). The total world bentonite production in 1980 was 6,625,204 st of which the United States produced 4,184,619 st. The total fuller's earth produced in the world in 1980 was 1,952,600 st of which 1,534,000 st was produced in the United States. The total world kaolin production in 1980 was estimated at 22,971,000 st. Approximately 13 million tons of this total was produced in two countries, the United States and the United Kingdom. World production data on ball clay, fire clay, and miscellaneous clays are not available. As industrial minerals, clays are a complex group that consists of several mineral commodities, each having somewhat different mineralogy, geologic occurrence, technology, and uses. In this chapter these commodities are organized as follows: (1) bentonite and fuller's earth; (2) kaolin, ball clay, halloysite, and refractory clays; and (3) miscellaneous clay and shale. Bentonite and fuller's earth are grouped together because they are closely interrelated. Though bentonite is a term based on mineral composition and fuller's earth is a term based on use, the two are virtually inseparable, because much, but by no means all, clay sold as

Jan 1, 1983

By J. P. Zannaras

If we assume that testing of the model started at time t1, that time t2 was the instant at which the elastic limit of the material was passed at the points of the maximum stress, and that at time t3 the failure of the model was completed as shown in the typical photograph (Fig. I), then actually the photograph shows the final results of what took place between time t2 and t3. By the principle of similitude all conclusions drawn from the behavior of the model were applicable and valid for actual mine roofs up to the time t2 after time t2 motions and friction came into play during the destruction of the model. The events between t2 and t3 and the photograph itself are not representative of the events expected in an actual mine roof unless corrections imposed by the principle of similitude are made. It can therefore be stated that all conclusions and observations of the author from the photographs are applicable only to his small-scale experiments and not to actual mine roofs. Without investigation of the true stresses or the combined maximum stress induced in the beds it has been assumed by the author that failure occurred solely due to bending stress caused by loading of the beds. (The author states that D was taken to be the theoretical strain given by elementary beam theory.) Fig. 1 shows that the beds were clamped together. This clamping was necessary to prevent slipping of the beds due to the horizontal shear. However, this clamping caused an uncontrolled and undertermined com-pressive stress S at the end of the beds, and this com-pressive stress caused a true tensile stress ?S (? PoiS-son ratio). This true tensile stress combined with the tensile stress produced at the end of the beds due to loading (fixed beam uniformly loaded), and the horizontal shear which is maximum at the end appears to be the most probable maximum stress produced in the beds. It is therefore evident that the author in performing this experiment has introduced stresses in the model not existing or dissimilar to actual mine roofs, and therefore his conclusions may be applicable to his small-scale experiments but inapplicable to actual mine roofs. Louis A. Panek (author's reply)—The discusser's difficulty is related to mechanics of materials rather than similitude. The effect of time did not enter into the tests in any way, being simply excluded from consideration in this study, because all quantitative results entering into the design equations are based on measured (between time t1 and time t2) bending strains that are well below the breaking strain (i.e., in the elastic range) for the model material. Time does not become a significant factor until the rock stress reaches 80 to 90 pct of the breaking stress." The objective of this investigation was not to predict that bolted roof will fail after time t, but was instead to develop a safe roof, that is, roof in which the bending strain is much less than (say less than 50 pct of) the breaking strain. The models tested to failure had no direct bearing on the design equations. Although for some rocks time may influence the magnitude of stress or strain at which fracture occurs, breaking stress or breaking strain values were not employed for any purpose in this investigation. Tests to failure were made only to obtain additional evidence regarding the behavior of the bolted laminae' as evinced by the locations of the cracks. As stated in the last paragraph of the paper, if one wishes to predict the time or stress conditions for which the roof will fail, additional information is required. The reason for clamping the laminae is to make them behave liked clamped beams, which closely approximate the clamped-plate condition of the actual mine roof beds (Ref. 5, p. 1; Ref. 6). The writer has previously demonstrated that the bending stresses induced by centrifugal loading in a clamped model beam are in agreement with those predicted by theory;' hence clamping has no significant effect on the bending ,trains measured during a test. Moreover, actual ,in, roof beds are likewise subject to a clamping effect due to weight of superincumbent strata carried by the pillars or ribs on each side of the opening. It is emphasized that this study was restricted to evaluating only the reinforcing effect produced by the bolts. Procedures employed in the experiments and in the data analysis were such as to deliberately exclude the effects of other factors. The objections offered are therefore invalid,

Jan 1, 1957

By E. M. Porbansky

DURING the investigation of the electrical transport properties of copper, it became necessary to prepare large single crystals of the highest obtainable purity. In an effort to meet these demands, single crystals of copper have been grown by the conventional pulling technique—as has been used for the growth of germanium and silicon crystals.' Low-temperature resistance measurements made on these crystals show that, as far as their electrical properties are concerned, they are generally of significantly higher purity than the original high-purity material. The use of these pure single crystals with very high resistance ratios has made possible the acquisition of detailed information regarding the electron energy band structure of copper2-' and has stimulated widespread effort on Fermi surface studies of a number of other pure metals. It is the purpose of this note to describe our method of preparing very pure copper crystals by the Czochralski technique. Precautions were taken to prevent contamination of the melt from the crystal growing apparatus. A new fused silica growing chamber was used to prevent possible contamination from previous groqths of other materials such as germanium, silicon, and so forth. A new high-purity graphite crucible was used to contain the melt. This crucible was baked out in a hydrogen atmosphere at -1200°C for an hour, prior to its use in crystal growth. Commercial tank helium, containing uncontrolled traces of oxygen, was used as the protective atmosphere. A trace of oxygen in the atmosphere appears to be necessary for obtaining high-purity copper single crystals. A 3/8-in-diam polycrystalline copper rod of the same purity as the melt was used as a seed. The copper rod was allowed to come in contact with the melt while rotating at 57 rpm. When an equilibrium was observed between the melt and the seed (that is, the seed neither grew nor melted), the seed was pulled away from the melt at a rate of 0.5 mils per sec. As the seed was raised, the melt temperature was slowly increased, so that the grown material diminished in diameter with increasing length. When this portion of the grown crystal was -1 in. long and the diameter reduced to less than 1/8 in., the melt was slowly cooled and the crystal was allowed to increase to - 1-1/4 in. diam as it was grown. By reducing the diameter of the crystal in this manner, the number of crystals at the liquid-solid interface was decreased until only one crystal remained. Fig. 1 shows a typical pulled copper single crystal. The purity of the starting material and the crystals was determined by the resistance ratio method: where the ratio is taken as R273ok/R4.2ok. The starting material, obtained from American Smelting and Refining Co., was the purest copper available. Most of the pulled copper crystals had much higher resistance ratios than the starting material. The highest ratio obtained to data is 8000. Table I is an example of the data obtained from some of the copper crystals. Note that Crystal No. 126 had a lower resistance ratio than its starting material and this might be due to carbon in the melt. The melt of this crystal was heated 250" to 300°C above the melting point of copper. At this temperature it was observed that copper dissolved appreciable amounts of carbon. The possible presence of carbon at the interface between the liquid and the crystal will result in reducing conditions and negate the slight oxidizing condition required for high purity as discussed below. The possible explanations of the improvement in the copper purity compared to the starting material are: improvement in crystal perfection, segregation, and oxidation of impurities. Of these, the latter seems to be most probable. A study of the etch pits in the pulled crystals showed them to have between 107 and 108 pits per sq cm. The etch procedure used was developed by Love11 and Wernick.10 The resistivity of the purest copper crystal grown was 2 x 10-10 ohm-cm at 4.2oK; from the work of H. G. vanBuren,11 the resistivity due to the dislocations would be approximately 10-l3 ohm-cm, which indicates that. the dislocations in the copper crystals would contribute relatively little to the resistivity of the crystals at this purity level. Segregation does not seem likely as the reason for purification of the material, since the resistivity of the first-to-freeze and the last-to-freeze portions are approximately the same, as was observed on Crystal No. 124. On most of the crystals that were examined, the entire melt was grown into a single crystal. If the

Jan 1, 1964

By E. J. Roberts

Prior to 1916, thickening was an art, and any accurate decision as to what size of machine to install to handle a given tonnage of a specific ore must have been one of those intuitive conclusions, based on both intimate and extensive acquaintance with thick-ners and ore pulps. Then in 1916 "knowledge of acquaintance," became "knowledge about" with the publication of the Coe and Clevenger paper.' The unit operation of thickening had graduated to the status of an engineering science. The fundamental similitude relationships for the two major phases of the operation were defined so clearly that batch tests on models as small as liter cylinders could serve to specify protypes as large as 325 ft in diameter. It is quite apparent from reading the literature that Coe and Clevenger's contribution is not generally appreciated. In so far as the basic engineering relationships are concerned, the only real advance which has occurred in the 30 odd years which have elapsed since the Coe and Clevenger paper is the recognition of the effect of the rakes on the thickening process. Bull and Darby2 noted this in 1926, and the extensive use of the "gluten type" thickener, in which the effect is magni-fied, bears witness to its importance. Comings3 further verified this effect of the rakes. As a matter of fact, a number of papers show an apparent regression from the Coe paper in that the area determinations are made on the basis of a single test from One concentration of solids. Coe and Clevenger amply demonstrated that this is unsafe, since the controlling zone may be one other than that of the feed dilution. Comings3 neatly demonstrated this without apparently realizing it. Of course there have been significant advances in the application of the operation to industry. Open tray thickeners were introduced to save area; balanced tray thickeners, washing thickeners, and multifeed clarifiers were developed with all of their special hydraulic and mechanical problems. Combinations of all kinds have been introduced, such as combination agitators and thickeners, combination flocculators and clarifiers, combination thickeners and filters. With the establishment of the operation on a firm engineering foundation, installation was facilitated and expansion proceeded. There are still problems, of course, functional as well as mechanical. Sometimes the moisture in the underflow obtained in practice is not as low as is expected on the basis of the test data. Sometimes the underflow is so "thick " that its discharge and subsequent handling requires special attention. Island formation plagues some operators. The use of the thickener as a surge basin and blending tank in the cement industry poses unusual problems. Design of rakes and the drive mechanism must be continually im-proved. Corrosion problems must he overcome. Power requirements for raking the settled solids occasionally is the controlling factor as it was in the case of the all American Canal desilting installation. Other similitude relationships and design problems come into the picture when we enter the field of clarification or nonline settlement. We have an energy dissipation problem in introducing the feed and any models must satisfy the Froude model relationships. Autoflocculation requires detention which involves the same similitude laws that we encounter in the compression zone. Approach to an Exact Science The next step beyond having control of the similitude relationships is to understand the why of these relationships right back up the line to first principles. The ultimate might be that, if given the mineralogical composition of the solids and their size distribution together with an analysis of the suspending liquid, we could calculate the entire thickening behavior of the system. Then we could say we had reduced the operation to an exact science. True it might be more trouble getting this basic analytical data than to make our empirical determinations for area and volume, and we would need an ENIAC to calculate the results, but that does not detract from the desirability of such understanding. Considerable work has been done by the chemical engineers with this end in view. Comings,3 Egolf,4 Work,5 Kam-mermeyer,6 Steinour,7 and others have studied the problem. The writer has no final answer to the thickening story but would like to propose a picture of the mechanics of the two phases of thickening which has been found useful in understanding the subject and which leads to some convenient relationship in treating the compression step and arriving at the compression depth.

Jan 1, 1950