Search Documents

By H. M. Strong

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

Jan 1, 1965

By Irwin Politziner F. M. Townsend, L. S. Reid

Recently published data indicate that the water vapor content of a gas, as determined by dew point measurement, is inaccurate when the gas has been dehydrated with diethylene glycol. Water vapor contents and dew point measurements of gases dehydrated with triethylene glycol have been obtained in this investigation in order to determine the magnitude of a similar error, if one exists. Experimental data show a very low concentration of triethylene glycol vapor in gases dehydrated at atmospheric ternperatures and pressures ranging from 500 to 2,500 psia and that the accuracy of dew point measurements is not impaired by the presence of triethylene glycol vapor. INTRODUCTION Nearly all natural gas transported by pipe line to northern and eastern markets must be dehydrated to a low residual moisture content in order to assure maximum transmission efficiency and continuous delivery under peak load conditions. Without dehydration, water vapor often condenses in pipe lines in quantities sufficient to restrict the flow3 and. under low atmospheric temperature conditions, gas hydrates may form and plug gathering. transmission and distribution facilities. Gas may be dehydrated to pipe line specifications by a number of methods. The most widely used methods are (1) adsorption of water vapor on a solid dessicant such as activated bauxite, alumina or silica gel,' (2) by absorption of water vapor by either diethylene or triethylene glycol-water solutions of high glycol concentration,11,13 and (3) by simultaneous expansion-refrigeration of very high pressure gas in which gas hydrates are purposely formed and then quickly decomposed.'"" Of these three fundamental methods, the adsorption and absorption processes provide the bulk of the dehydrated gas moving to market at the present time. The expansion-refrigeration process. a relatively new development, is gaining in favor because of its simplicity but is sharply limited in its application by available differentials between source and gathering system pressures. Where gas sale or purchase contracts contain a maximum water vapor content specification, the water vapor content is usually determined by measuring the dew point of the gas at system pressure. Knowing he dew point temperature and the system pressure, the water vapor content of the gas is determined from a graphical correlation of saturation, or dew point, temperature us the water vapor content of the gas which is customarily expressed as pounds of water vapor per million cubic feet of gas measured at 14.7 psia and 60°F. Dew points are usually measured with the U. S. Bureau of Mines Dew Point Tester. shown in Fig. 1, which is well suited for use at temperatures lower than 120°F and at pressurei; ranging from atmospheric to 3,000 psi. A number of correlations of dew point temperatures vs water vapor content are available,3,5,15 the most recent of which is that by McCarthy et al.4 Where natural gas is dehydrated by absorhing contact with concentrated diethylene glycol-water solutions. considerable difficulty is experienced in observing water dew points due to condensation of liquid hydrocarbon and glycol films on the mirrored surface of the dew point tester. Riesenfeld and Frazierl2 report that the dew point method for determining the water vapor content of a dietllylene glycol-treated natural gas is in error due to the glycol content of the condensate formed on the tester mirror, and that the water vapor content of the gas tested is actually lower than that indicated by the standard' dew point-water vapor content correlation. Extensive experience in testing gases dehydrated by concentrated triethylene glycol-water solutions lias not revesled difficulties of the nature cited above for diethylene glycol. Water dew points are clear and sharp and condensation of

Jan 1, 1951

By J. W. Rutter, K. T. Aust

THE object of the experiments to be described in this report was to determine, first, which grains, out of a large number introduced into a sample in which their growth could proceed, were able to grow most successfully; and, second, the factors which control their growth selection. In an earlier study of grain boundary migration in high-purity lead,' it was found that, under some conditions, the most successful grain in an experiment of this type bears a special orientation relationship to the grain into which it is growing. The present work represents an extension of the previous study, directed at the determination of the orientation relationships which appear under a wider variety of conditions. Experiments of this general type have been done using commercial aluminum2 with somewhat conflicting results. In the present experiments, in contrast to the above,' the sample composition was a principal controlled variable, resulting in a considerable clarification of the phenomena observed. EXPERIMENTAL TECHNIQUE The experimental technique used in this investigation was essentially that employed previously in studies of grain boundary migration kinetics.1'3 The material used was zone-refined lead with various small amounts of tin or silver added to study the effect of varying the sample composition. The specimens were in the form of single crystals, about 1/4 in. 1/4 in cross section by 3 to 4 in. in length, grown from the melt in a horizontal graphite boat. A crystal grown in this way generally contained a lineage or "striation" substructure of the type shown in Fig. 1 of Ref. 1. This substructure consists of an array of low-angle boundaries which partition the crystal into regions misoriented from each other by a few degrees. The substructure constitutes a source of driving energy, estimated to be about 4000 ergs per cc, for the migration of a large-angle grain boundary into the striated crystal. New grains were introduced by plastic deformation in compression of a localized reglon at one end of the sample. Re-crystallization occurs at room temperature in the deformed region, introducing a large number of new, striation-free grains into the specimen. These grains are able to grow into the melt-grown crystal on annealing at a suitable temperature. The stria-tions are removed by passage of the grain boundary, as the substructure energy is utilized to produce the boundary migration. Usually, only one or two of the recrystallized grains persist after growth has occurred for about 1 cm out of the deformed region. Laue back-reflection X-ray photographs were obtained to determine the orientation relationship between the melt-grown crystal and the recrystallized grain or grains which grew most successfully in each specimen. EXPERIMENTAL OBSERVATIONS Fig. 1 shows the results obtained for a series of specimens to which had been added tin in the range from 0 to 5 ppm by weight. In this figure is plotted, in the standard stereographic triangle, the single axis about which the most successful recrystallized grain or grains are related to the starting, melt-grown crystal by the smallest rotation. The range of angular rotations observed about such axes is indicated beside each stereographic triangle. The results shown in the stereographic triangle on the left, Fig. l(a), were obtained by localized deformation at one end of each specimen at room temperature, followed by annealing in a furnace which had been preheated to 300°C. The furnace atmosphere was argon. It will be seen that the orientation relationships observed between the melt-grown and recrystallized grains are quite random. On the right, Fig. l(b), are shown results for samples in which the de-

Jan 1, 1961

By David W. Cannon, Russell V. Mann

In a dispersion of particles in liquid a net charge will develop at the particle-liquid interface. This surface charge is usually due to the adsorption of charged material from solution. The existence of this surface charge gives rise to the formation of the electric double layer of counter-ions which surrond the charged particle. The particle surface charge and the electrostatic repulsion which exists between similarly charged particles is the primary stabilization mechanism for lyophobic colloids (1). The separation of charge which occurs at the particle-liquid interface gives rise to several dynamic phenomena associated with colloidal systems or with solid-liquid interfaces in general. These phenomena are known as electrokinetic phenomena and the four classic electrokinetic phenomena are; electrophoresis, electroosmosis, streaming potential, and sedimentation or Dorn potential. The actual driving force for electrokinetic phenomena is not the surface charge per se, but the charge at the interface between the liquid which is hydrodynamically bound to the particle surface and the bulk fluid. This interface is known as the slipping plane or the plane of shear and the potential at this interface is the zeta potential (2). The factors linking the electrokinetic phenomena is that they involve a relative motion between the liquid and the charged particle or solid surface and the driving force is the zeta potential of the solid. In addition to the four classic electrokinetic phenomena there are two additional electrokinetic effects in disperse systems; the electro-acoustic effects. When an alternating electric field is applied to a cooloidal dispersion the particles will move in the field due to their net zeta potential. If there is a density difference between the particles and the fluid this motion will result in the development of an acoustic wave. The effect was discovered at Matec and has been termed the Electrokinetic Sonic Amplitude or ESA (3). ESA is the pressure amplitude generated by the colloid per unit applied electric field strength and has SI units of pascals per volt per meter. When an alternating pressure field (acoustic wave) is applied to a colloidal dispersion the inverse of the ESA effects occurs. A density difference between the disperse phase and the continuous phase leads to a relative motion between the particles and the surrounding liquid. This means that there will be a periodic displacement between the charged particle and the oppositely charged counter-ions in the electric double year. This displacement results in the development of an alternating dipole moment at the frequency of the applied field. This effect is termed the Ultrasonic Vibration Potential or UVP and was first predicted for electrolyte solutions by Debye in 1933 (4). UVP is measured in units of volts per unit velocity amplitude of the applied acoustic field or volts per meter per second. In 1938, Rutgers (5) and Hermans (6) pointed out that the effect would also be present in colloidal dispersions. A detailed theory for UVP effects in colloids, also called CVP, was first presented by Enderby in 1951 (7). Extensive studies of the UVP in electrolytes have been carried out by Yeager et al (8). Recently, O'Brien (9) has developed a general theoretical treatment of electro- acoustic effects in colloids and has derived a reciprocal relation linking ESA and CVP effects. The most commonly studied electrokinetic phenomena is electrophoresis. Electrophoresis is the movement of charges particles in an applied electric field. The velocity of the particle divided by the applied electric field strength is the electrophoretic mobility of the particle. The zeta potential can be calculated from this mobility (2). The magnitude of both the ESA and CVP effects are directly proportional to the electrophoretic mobility of the particles. The mobility determined by the two electro-acoustic effects is the dynamic or AC mobility of the particles.

Jan 1, 1990

By L. S. Reid, Irwin Politziner F. M. Townsend

Recently published data indicate that the water vapor content of a gas, as determined by dew point measurement, is inaccurate when the gas has been dehydrated with diethylene glycol. Water vapor contents and dew point measurements of gases dehydrated with triethylene glycol have been obtained in this investigation in order to determine the magnitude of a similar error, if one exists. Experimental data show a very low concentration of triethylene glycol vapor in gases dehydrated at atmospheric ternperatures and pressures ranging from 500 to 2,500 psia and that the accuracy of dew point measurements is not impaired by the presence of triethylene glycol vapor. INTRODUCTION Nearly all natural gas transported by pipe line to northern and eastern markets must be dehydrated to a low residual moisture content in order to assure maximum transmission efficiency and continuous delivery under peak load conditions. Without dehydration, water vapor often condenses in pipe lines in quantities sufficient to restrict the flow3 and. under low atmospheric temperature conditions, gas hydrates may form and plug gathering. transmission and distribution facilities. Gas may be dehydrated to pipe line specifications by a number of methods. The most widely used methods are (1) adsorption of water vapor on a solid dessicant such as activated bauxite, alumina or silica gel,' (2) by absorption of water vapor by either diethylene or triethylene glycol-water solutions of high glycol concentration,11,13 and (3) by simultaneous expansion-refrigeration of very high pressure gas in which gas hydrates are purposely formed and then quickly decomposed.'"" Of these three fundamental methods, the adsorption and absorption processes provide the bulk of the dehydrated gas moving to market at the present time. The expansion-refrigeration process. a relatively new development, is gaining in favor because of its simplicity but is sharply limited in its application by available differentials between source and gathering system pressures. Where gas sale or purchase contracts contain a maximum water vapor content specification, the water vapor content is usually determined by measuring the dew point of the gas at system pressure. Knowing he dew point temperature and the system pressure, the water vapor content of the gas is determined from a graphical correlation of saturation, or dew point, temperature us the water vapor content of the gas which is customarily expressed as pounds of water vapor per million cubic feet of gas measured at 14.7 psia and 60°F. Dew points are usually measured with the U. S. Bureau of Mines Dew Point Tester. shown in Fig. 1, which is well suited for use at temperatures lower than 120°F and at pressurei; ranging from atmospheric to 3,000 psi. A number of correlations of dew point temperatures vs water vapor content are available,3,5,15 the most recent of which is that by McCarthy et al.4 Where natural gas is dehydrated by absorhing contact with concentrated diethylene glycol-water solutions. considerable difficulty is experienced in observing water dew points due to condensation of liquid hydrocarbon and glycol films on the mirrored surface of the dew point tester. Riesenfeld and Frazierl2 report that the dew point method for determining the water vapor content of a dietllylene glycol-treated natural gas is in error due to the glycol content of the condensate formed on the tester mirror, and that the water vapor content of the gas tested is actually lower than that indicated by the standard' dew point-water vapor content correlation. Extensive experience in testing gases dehydrated by concentrated triethylene glycol-water solutions lias not revesled difficulties of the nature cited above for diethylene glycol. Water dew points are clear and sharp and condensation of

Jan 1, 1951

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By R. F. Blakely, P. E. Potter

Any stratigraphic section or bedding sequence can be synthesized if there is a transition procedure from one lithology or bedding type to another, and if thickness distributions of the different lithologies are known. Stratigraphic sections of a fluvial sandstone body were synthesized with five bedding types: cross-bedding, massive beds, parting lineation, ripple mark and mudstone. The transition procedure from one bedding type to another used dependent, Markovian random processes which have a memory that extends one step backward in the depositional process. As observed in nature, median grain size and sand wave thickness (cross-bedding and ripple mark) decline upward in the synthesized sections as proportions of the different bedding types change. Grain size and permeability were also incorporated into the sections. By changing the transition procedures, bed thickness distributions, rate of upward decline or sand wave height and length, different types of sections can be synthesized, thus making it possible to model many different sedimentation problems. INTRODUCTION This paper describes a general method for synthesizing stratigraphic sections and bedding sequences of sedimentary, metamorphic or igneous origin. Synthetic generation is of interest for several reasons. Close correspondence between real and synthetic sections suggests that the factors used in the synthesizing model may indeed be the correct ones, thus giving the investigator a check on his assumptions. Rapid, inexpensive simulation of many stratigraphic sections permits one to synthesize a rock body (sandstone or carbonate reservoir) or, on a larger scale, the fill of a sedimentary basin. Harbaughl gives an example of mathematical simulation of a carbonate basin. He simulated the basin in the hope that improved prediction would follow better understanding of the depositional processes. From the petroleum engineer's viewpoint it seems reasonable to believe that the synthetic generation of rock properties and their distribution in a reservoir should be relevant in the study of reservoirs. Any stratigraphic section or bedding sequence can be generated provided there is a transition procedure from one lithology or bedding type to another, and provided the thickness distributions of the different units are known. The transition procedure involves random processes that are either independent or dependent. If the depositional process is independent, previous deposition will have no influence on present deposition. However, if it is dependent, past deposition will influence either present or future deposition. Such a dependent depositional process can be thought of as having a memory that extends backward in time through one or more pulses of deposition. A process with a memory can be described by a Markov process. Because the concept of memory or dependence appears to be in accord with our understanding of many depositional processes, Markov processes were used to synthesize the bedding sequences of this study (see Appendix). The above methods are perfectly general and are appropriate for any stratigraphic section or bedding sequence: bedding types in a beach deposit, an evolving carbonate bank or the changing lithologic fill of a thick geosyncline sequence. We chose to synthesize a vertical profile of a fluvial sandstone body because its characteristics were well documented, much was known about fluvial processes and fluvial-deltaic sandstone bodies constitute an important class of petroleum reservoirs. CHARACTERISTICS AND ORIGIN OF FLUVIAL CYCLE The fluvial cycle has been well documented in recent years by Bersier,2 Allen3-5and Visher.6,7 Deposits from fluvial cycles range from 10 to 150 ft or more in thickness and are characterized by a "fining upwards": coarse sandstones with occasional conglomerates grade upward into medium- to fine-grained sandstone, and hence into siltstone and mudstone. The dominant sedimentary structure is inclined bedding: thick cross-beds in the lower half

By D. L. Garthwaite, F. M. Stewart, F. K. Krebill

Pressure has been maintained in the Elk Basin Ten-deep reservoir since the initiation of inert gas injection in Sept., 1949. Oil is being produced under conditions favorable for gravity drainage, including high angle of dip, appreciable structural closure, and fairly good per~neability. Results being obtained are high per cent recovery of oil-in-place, sustained field productivity, reduced operating problems, and recovery of plant products. The high recovery has been calculated by alpplication of the gravity drainage theory described in an earlier AIME paper.' It is also confirmed from field performance by comparing nil recovery to date with gas cap spuce voided. GENERAL INFORMATION The Elk Basin field lies in Park County, Wyo., and Carbon County, Mont., approximately 50 miles east of Yellowstone Park. Located at the northern end of the Big Horn Basin. the field is situated on an elongated asymmetrical anticline (Fig. 1). The pressure maintenance project covered by this paper involves the Embar-Tensleep reservoir in this multi-pay field. Discovered in Nov., 1942, it has a proved productive area of over 6,300 acres. It consists of the 210-ft thick Tensleep sandstone of Pennsylvanian Age overlain by 40 ft of Embar dolomite of Permian Age. The two formations apparently are in communication with each other and are produced as a common source of supply. As 98 per cent of the reserves are contained in the Tensleep sandFtone, in this paper the Embar-Tensleep reservoir will he called simply the Tensleep. The reservoir is found at an average depth of 4,900 it, is approximately 7 miles long and 2 miles wide, and has a maximum oil productive closure of 2,330 ft. The strata dip an average of 21" on the west flank and 45" on the east. There are presently 132 Tensleep wells in the field drilled on 40-acre spacing; of these, 17 are shut in because of high gas-oil ratio (GOR) and eight are gas injection wells. Rock characteristics and related information for the pay section follow: Porosity, avg 10.7 per cent Permeability (air) from cores, avg 118 mds Permeability from PI'S 91 mds Connate water, avg 8 per cent Initial pressure (-400 ft) 2,234 psia Pressure April 1954 (-400 ft) 1.327 psia At high structural positions much of the better sandstone has connate water content of only 2 to 5 per cent as indicated by the coring of wells with oil in nearby Tensleep fields. In special imbibition tests, samples readily imbibed oil but repelled water. In tests on crushed samples, the sand grains rapidly settled in oil and carbon tetrachloride but formed a stable unsinking floc in water. This information indicates that the Tensleep sand is preferentially oil wet. The oil produced averages 29' API gravity and has a high sulfur content. Extensive sampling and testing have established the fact that characteristics of reservoir oil in this field vary greatly with structure, see Fig. 2. Additional detail on the variation of oil characteristics with structure at Elk Basin may be found in Reference 2. At original reservoir pressure the oil was undersatu-rated with gas, the bubble point being 1,250 psi at the crest but only 500 psi at the lowest elevation sampled. As noted on Fig. 2, other fluid properties—reservoir volume factor. solution ratio, gravity, viscosity, and hydrogen sulfide content—likewise vary greatly. These variations complicate mathematical analysis, for it is necessary to average these characteristics in order to enter them into reservoir engineering expressions. To further complicate analysis, sand-face pressure also varies greatly because of the increase in fluid head down structure. As an approximation, average properties under initial conditions may be read from Fig. 2 at the -450-ft elevation.

Jan 1, 1956

By Peter A. Allen

INTRODUCTION Lac Minerals Ltd., a significant North American producer of gold and other metals for over three decades, recently completed the financing for its third discovery in an eight year period. From 1975 to 1982, three important gold deposits in Quebec and Ontario: namely Bousquet, Doyon and Hemlo, were located and tested by company geologists. When current development plans are complete, Lac will have injected $400 million into these deposits. A discussion of these projects and their financing reveals some extraordinary challenges faced by top mining management. Responses to these challenges developed by Lac management may interest executives planning to proceed with new North American gold or base metal mines. MINE DE BOUSQUET In 1976, we explored a property containing a large low grade surface gold deposit. The geologial setting was exciting and we tested recent rethinking about the volcanigenic genesis of gold and base metal occurrences. These ideas helped pus to persist in drilling a nearby conductor and 600,000 ton, 0.22 ounce/ton orebody was outlined. With gold at $Cdn.150 per ounce the feasibility was approved and work started February 1978. Availability of a company custom mill 35 miles away reduced our immediate capital needs to $8.5 million. First of all, understanding the geological significance of results at Bousquet was critical. Mine and surface exploration geologists at Lac are continually and closely working together to better understand the nature of gold deposits in the Shield. This management emphasis on both pure and mine geology meant that a much larger ore reserve could confidently be predicted. (A proposal to quickly develop and ship the ore to a nearby company mine was discarded.) Instead of designing just for the initial plan of 500 tons per day we therefore included headframe, shaft and hoisting capacity for 1,700 tons and provided a compartment for future shaft deepening. These prebuilt additions were welcomed within two years. The type of collaboration whereby surface and mine geologists make study and field projects of each other's problems sets the foundation for this predictive ability. At Lac, the traditional responsibility for and supervision of the mine geologist by the mine manager is modified by this collaboration. This adds interest and scope to the geologist's work and lets him think beyond production requirements. Besides such management considerations, dedicated, professional geologists with integrity and lots of practical experience are needed to plan successful developments like Bousquet. Second, projecting future bullion revenue was critical and our experience selling production for future delivery gave us a way to ensure that loan repayment at term would not be jeopardized by falling metal prices. Output from a single mine is often reduced or curtailed for numerous reasons, so "selling forward" is done best and with least risk by companies like Lac with more than one or two producing locations for that metal. At first, a forward marketing program was considered unwise, speculative and unorthodox by analysts and other producers. Declining markets in the early 1980s prolonged gold project capital recoveries however, and critics were better able to appreciate the usefulness of such a tool. The notion that gold miners shouldn't sell forward is rooted in the belief that long-run gold prices will rise as modern social theory implementation in democracies inevitably erodes currency values. Selling against this trend goes against this industry belief or "ethic" and is at odds with the company shareholder who expects the gold price to rise not fall. Nonetheless, an efficient tool is available to the gold producer (not to most other metal producers) to quickly disengage from adverse price trends or to lock in spot prices instantaneously for two years or more of profits. Gold had recently retreated from U.S.$200 to $104 and now had reversed to $150. Lac had been forced to close several of its producers during the late 1960s because central bank-controlled

Jan 1, 1985

By Rajendra K. Devala, Sukumar Bandopadhyay

In order to be able to design an adequate mine ventilation system it is necessary to understand how pollutant masses disperse into turbulent ventilation air streams. The mixing of gases and dusts with ventilating air, and dispersion of these pollutants in mine airways, are important flow or transfer problems. For the solution of such problems, the finite difference technique is often used. In this paper, numerical approximation schemes for the solution of the mass transfer problem have been developed using several explicit finite difference schemes. Results from these schemes have been compared with the analytical solution for spatial and time increments using the concept of Root Mean Square error. The comparison of various schemes reveals that the explicit finite difference schemes analyzed here give more accurate results at lower velocities, though the results are dependent on the spatial and time steps considered. At higher velocities the schemes tends to oscillate about the analytical solution. INTRODUCTION In mining, most of the mine air pollutants, such as respirable dust, methane, blasting fumes and diesel exhaust, are produced at the mine face and in haulageways. Injection of any of these pollutants in the mine ventilating air causes disturbances of the equilibrium state of the pollutant concentration of the ventilating air. The mixing of these gases and dusts with the ventilating air are typical examples of mass transfer processes. Pollutants are diluted by the stream of ventilating air flow along airways and discharged into the outside atmosphere of the mine through a return shaft. The dispersions of these pollutants in mine airways are important flow or transfer problems. Knowledge about pollutant sources, emission and transport of these pollutants are essential for ventilation planning and design to ensure good engineering control over air quality, quantity, and environmental conditions. Modeling of mass transfer phenomena in mine ventilation is largely based on extensive studies of turbulent dispersion in wind tunnels and the lower atmosphere of the earth. The general approach in the formulation of models lies in the recognition that these mass transfer phenomena can be conceptualized as a flow problem that can be solved by the general application of Taylor's theory (1954) of turbulent flow and diffusion with modification appropriate static and dynamic parameters. If the face zone is initially contaminated by a pollutant over a long section (e.g. instantaneous exposure of large gas accumulation, rapid methane release etc.), the process of decontaminating the mine workings can be regarded as displacement and diffusion of the tailing edge of the cloud of contaminates in the form of a semi-infinite gas body, for which a one-dimensional mass transfer equation is represented by: [at + u • ax = Ex ex2 + f (x,t)(1.1)] In Eq. (1.1), the term [~~] is the rate of growth of the concentration in the differential element, while a [(5xx )] is the net gain of material due to convective transfer. The two terms balance the total loss of material due to turbulent dispersion that is represented by [Ex . 21] and f(x,t), the source term for the pollutant in the roadway. Several cases of mass transfer that occur in mining can be identified and modelled using the conservative volume element method. The governing equation, along with the associated boundary and initial conditions, source terms and their solutions are detailed elsewhere (Bandopadhyay and Ramani, 1988, Devala, 1989). Only one mass transfer situation is presented (Figure 1) here for the present analysis. Mixed Evolution of Pollutant (multiple pulses) Typical examples of these situations are emission of pollutants from multiple diesel vehicles in the face area, and gas emission from standing pillars as well as from broken coal seam during coal cutting with a continuous miner. The mathematical statement of the problem along with the initial and boundary conditions are given below, Governing Equation: [do+u- =E d2c(1.2) ataxx dx2] Initial conditions:

Jan 1, 1993

By Floyd C. Bossard

A method of hand-calculating the air-cooling requirements of a conceptualized underground mining operation is presented for the reader's orientation. Separate air heat load calculations were conducted for adiabatic compression, electromechanical equipment, wallrock, broken rock, groundwater, and blasting operations. The air heat sources were calculated for four mining levels under conditions representative of a typical planned mining operation at depth. The total heat gain on each level will approximate the level air-cooling requirements. The results of these hand calculations can be further modified by the use of mine ventilation computer software that refine the heat-source calculations, predict underground ambient air temperatures, and establish the air-cooling requirements of a mine. INTRODUCTION The principal mining method employed is a modified vertical crater retreat (VCR) blasthole operation. Typical scopes range from 20 to 40 feet wide, 75 to 100 feet long, and 100 feet in height. Six and one-half inch blastholes are drilled with an "in the hole" hammer drill. Eight-foot high horizontal rounds are blasted down. Mucking of ore from the undercut to orepass is done with LHD equipment operated from a remote control station. Backfill includes hydraulically placed tailings with cement, and waste rock when available. Stope access is from ramp sub-levels on 50-foot vertical intervals. Crosscuts are ramped down to the first stope cut 25 feet below the sub-level elevation. Then the crosscuts are raised by taking down the back when each stope cut is completed, until an elevation 25 feet above the sub-level is reached. The crosscuts are filled with tailings, and/or waste rock. See Figure 1. Typically, a two-pronged approach to defining the air-cooling requirements is conducted. First, the principal sources that make up the air heat load are individually hand-calculated. Second, projections of mine heat load are calculated by utilizing computer modeling techniques. This paper discusses the first method (hand calculations) of determining the individual components of heat flow into the mine. CALCULATED AIR HEAT LOADS Adiabatic Compression The plans for the mine include delivery of 300,000 cfm of air to ventilate the lower level operation. This is equivalent to approximately 200 cfm/ton of rock produced (300,000 cfm/1500 tpd of ore and waste). As air flows down a shaft, with no heat interchange between the shaft and air and no evaporation of moisture taking place, the air is heated in the same way as if it were compressed in a compressor. Dry air increases in temperature about 5.4°F per 1000 feet. One BTU is added to each pound of air for every 778 feet of decrease in elevation, or is subtracted for the same elevation increase. For dry air, the dry-bulb temperature change is 1/(0.24 x 778) = 0.00535°F/ft., or 1°F/187 ft. elevation. Auto-compression may be masked by the presence of other heating or cooling sources, such as shaft wallrock, groundwater, air and water lines, electrical facilities, etc. The major factors influencing the temperature of the air delivered underground by a shaft are (1) the night time cool air temperature's effect on the rock or lining of the shaft, (2) temperature gradient of ground rock related to depth, and (3) evaporation of moisture within the shaft which increases the latent temperature and decreases the sensible temperature. [Calculation For Adiabatic Heat of Compression] a. Assumptions: 1. Three hundred thousand cfm of fresh air at 3000 level has increased in temperature during the summertime to the point where it has little available cooling power. Air-cooling will be required on 3000 Level and below. 2. Elevation of 3000 Level is approximately +1000 ft. above sea level.

Jan 1, 1993

By R. A. Morse, C. R. Holmgren

The production of oil by water flooding can be substantially increased by the maintenance of free gas saturation in the reservoir during the flooding operation. This effect is accomplished by the alteration of oil relative permeability characteristics and the occupation by gas of pore space that would otherwise be filled with residual oil. The amount of reduction In residual oil can be calculated from appropriate water-oil relative permeability characteristics. This paper presents experimental data in support of the foregoing conclusions and an example of the calculations. The microscopic pore saturation concepts of the mechanism are discussed. A method of practical application to field floods is presented together with discussion of certain limitations. INTRODUCTION The presence of free gas has been reported by a number of investigators to significantly affect the oil recovery which can be obtained from sandstone flow systems by water flooding.1,2,3,4,5 The effect of gas, noted in every instance, has been to cause lower residual oil saturations than could be obtained by water flooding the same systems in the absence of free gas. The degree of improvement in recovery has been observed to vary widely, depending on the systems used and the conditions of the tests. The increased oil recovery obtained because of the presence of gas during a water flood has been variously attributed to changes in physical characteristics of the oil, selective plugging action of the gas, inclusion of oil mist in the free gas phase, and the additional sweeping or driving action of the free gas. All but the first of these suggestions imply changes in the displacement mechanism. The change in viscosity and inter-facial tension of the oil phase, within the pressure range used for all the experimental work, is certainly not sufficient to account for the differences in residual oil saturation noted dess there is a drastic change in the displacement process. One other effect which logically seems capable of causing differences in residual oil saturation of the magnitude noted in the experimental work is that of simple replacement. In a water-wet system containing oil, water, and gas, it is to be expected that the gas will exist inside the oil. This is the position of minimum free surface energy, since the gas-oil interfacial tension will be less than the gas-water interfacial tension. There is no apparent reason to expect that the existence of free gas within the oil phase should alter the saturation at which the non-wetting phase (now oil and gas) should become discontinuous and hence trapped so as to be unrecoverable by direct displacement by water. If this is the situation, then trapping of a certain percentage of gas saturation during water flood should result, at infinite water-oil ratio, in a like reduction of oil saturation below that attainable by flooding in the absence of free gas. It is visualized that the gas will exist as bubbles inside the discontinuous residual oil as illustrated in Fig. 1, with the size of the oil bubbles being substantially unchanged due to the presence of the gas. As a practical matter, it can be anticipated that the presence of a free gas saturation inside the oil phase will reduce the relative permeability to oil which will exist at any particular water saturation. This reduction will be caused by two factors — -the addition of the gas-oil interface, and the reduction of area available for oil flow in the pores containing gas. This reduction in oil permeability at any particular water saturation will result in water breakthrough at a lower water satu-

Jan 1, 1951

By P. T. Bail, J. E. Bobek

In reservoirs containing soft shales interbedded with oil-producing sands, it is possible to have water movement from the shales into the sands during reservoir depletion. This paper presents the results of scaled flow-model tests performed to investigate the efficiency of oil displacement by this vertical influx of water into oil-producing sands. In these model tests, the effects of rate of influx and rock wettability were considered. Oil recoveries by shale-wafer influx and by water flooding were cornpared after 0.4 hydrocarbon volume of water entered the producing sand. This comparison was made for water-wet and oil-wet systems. In water-wet sand, shale-water influx produces slightly less oil than conventional water flooding; in oil-wet sand, shale-water influx produces about one-half as much oil as conventional water flooding. In another type of comparison, it was found that shale-water influx followed by water flooding recovers less oil than water flooding alone. This is true for both water-wet and oil-wet sands. A comparison of tests run at different rates of shale-water influx shows that oil recovery is insensitive to rate in water-wet sand; in oil-wet sand, oil recovery decreases with decreasing rate of influx. INTRODUCTION Localized subsidence of the earth's surface has accompanied the depletion of some oil reservoirs, particularly in California,1, 2 along the Texas Gulf Coast,3 4 and in the Bolivar Coastal fields of Venezuela. Subsidence occurs when fluid withdrawals reduce the pressure in the reservoir. The reasons for subsidence accompanying a reduction in reservoir pressure are not known with certainty.10° However, there is considerable evidence that compaction of soft shale beds adjacent to producing oil sands can be a major cause of this subsi-dence.3, 6, 10-12 When the shales compact, a limited amount of water moves from the shales into the producing sands. A schematic illustration of this oil-producing mechanism is shown in Fig. 1. As the pressure in the sand declines due to oil production, a predominantly vertical pressure gradient develops at the sand-shale boundary. This causes water to move from the shale into the sand as indicated by the solid arrows. This paper reports the results of scaled flow-model tests to investigate the efficiency of oil displacement by the vertical influx of shale water into producing sands. The model tests were not performed with the aim of duplicating the recovery behavior of any specific reservoir. The primary purpose of the tests was to obtain a better understanding of oil displacement by shale-water influx. A knowledge of the oil- and water-production characteristics of even a simplified reservoir system should be of value to engineers estimating the performance of reservoirs where shale-water influx is occurring in conjunction with other production mechanisms. The following limitations and assumptions apply to the model work. 1. The sand of the prototype reservoir is homogeneous and horizontal. 2. The shale-water influx is uniformly distributed along the top and bottom of the sand member being studied. 3. No effects of edge-water drive or dissolved-gas drive that might accompany shale-water influx are considered. PROTOTYPE RESERVOIRS Scaled-model tests were performed to study oil displacement from three prototype reservoir sections. The reservoir sections modeled consist of single sand layers out of a sequence of -sand and -shale layers, as illustrated in Fig. 1. Each prototype represents a homogeneous, high-permeability, unconsolidated sand layer. In the field, the prototype is oriented in such a manner that one end is at the radius of drainage of a producing well and the other end is pointed toward the producing well (Fig. 2). The sand and fluid properties and the dimensions of the three prototypes are given in Table 1. The capillary pressure relationships for the three prototypes are presented in Figs. 3, 4 and 5. The oil-water per-

By Rex B. McLellan, W. A. Oates, William W. Dunn

Vapor transport experiments have been carried out in order to determine the saturation solubility of carbon in cobalt and nickel with respect to graphite over a large temperature range. Some of the uncertainty existing in the thermodynamics of these systems has been removed and accurate estimates of the heat and entropy of solution of carbon in nickel and cobalt have been deduced. The relatitle partial entropies of carbon are considered to be predominantly vibrational in origin. ALTHOUGH the thermodynamic properties of solid solutions of carbon in cobalt and in nickel have been investigated previously, the results of these investigations show a marked lack of agreement amongst themselves. In this investigation simple vapor transport experiments have been undertaken in order to measure the temperature variation of the saturation solubility of carbon in nickel and cobalt with the object of removing to some degree the uncertainties in the thermodynamics of these systems. In the case of the C-Co system, early measurements of the saturation solubility of graphite in Co show results at the same temperature differing by almost an order of magnitude.1-5 Recently Schenck, Frohberg, and Jaspert6 measured the equilibrium between CH4-H2 mixtures and C-Co solid solutions at five temperatures ranging from 850° to 1050°C. Their values of the activity of carbon in the solid solution, deduced from a knowledge of the chemical potential of carbon atoms in the gas, varies linearly with composition at each temperature and the authors concluded that the solutions were Henryan. The saturation solubilities deduced from Schenck, Frohberg, and Jaspert's measurements are in good agreement with the solubility determined at 1000°C by Rao and Nicholson7 and the value at 1000°C deduced from the measurements of smith8 on the equilibrium between CO/CO2 mixtures and C-Co solutions. The object of this investigation was to check the results of Schenck, Frohberg, and Jaspert and also to extend the range of solubility data to higher temperatures. In the case of the C-Ni system, there is fairly good agreement between the values given by schaefer9 for the saturation solubility of carbon C, in nickel, the value given at 1000°C by Rao and Nicholson,7 and that deduced from Smith's8 measurements of the equilibrium between C-Ni solutions and CO/CO2 mixtures at 1000°C. The values of C, deduced from Schenck, Frohberg, and Jaspert's6 measurements of the carbon activity in nickel are also in fairly good agreement. However, there is a large discrepancy between the above set of measurements and those of Lander, Kern, and Beach," who measured the carbon solubility of nickel in equilibrium with graphite by depositing an excess of graphite on the surface of Ni-foils from a gaseous mixture and subsequently diffusing in the carbon. These results are considerably higher than those mentioned above. The determinations6 of carbon activity a, in nickel from the CH4-H2 equilibrium showed that at all the temperatures investigated, the ac-values varied linearly with concentration up to about a, = 0.6 and then began to increase with concentration much more rapidly. On the basis of auxiliary experiments in which C-Ni solutions were equilibrated with graphite, Schenck, Frohberg, and Jaspert concluded that their CH4-H2 solid equilibrations at high values of a, were spurious and a, behaved linearly with composition up to the saturation limit (a, = 1). The short linear extrapolations of Schenck, Frohberg, and Jaspert's a, data yield values of C, in good agreement with those of Lander, Kern, and Beach." The solubility data for nickel discussed here is presented in Fig. 1 in the form of plots of ln[Cc/(l - 2cc)] vs 1/T. The discrepancies mentioned are clearly seen. In this investigation a vapor transport method is used to measure C, over as wide a temperature range

Jan 1, 1969

By Robert M. Lewis

Can wet-ground mica be produced from mica schist ores, and if so, would its properties be comparable with products now on the market? This was the problem which prompted the development of a flowsheet for the beneficiation of mica schist. The fine-ground mica industry relies on weathered pegmatites and alaskite ores as sources of raw material. The mica concentrates obtained from these ores are relatively clean and require moderate beneficiation to meet color and bulk-density specifications. Mica obtained from mica schist is iron-stained and decrepitated and requires additional beneficiation to produce a wet-ground mica product. Despite these adversities, there are several advantages to be gained by mining and processing mica schist ores. The mica content of these ores is considerably higher than that of the pegmatites and alaskite ores now being mined, often containing in excess of 40% mica. This higher-grade mica feed should promote a lower mining and milling cost. The mica schist ores with their higher mica content would mean less tailings to handle; however, there would have to be a means for disposing of acid and sulfate wastes. A flowsheet was developed whereby mica could be recovered from mica schist ores and beneficiated in a way that would result in a high-quality end product.' The essential phase of the process involves using reducing reagents to increase the magnetic susceptibility of iron contaminants which then can be removed with a magnetic separator. General Treatment Procedure The mica is floated by either of two methods, an amine float in an acid circuit' or the U.S. Bureau of Mines (USBM) amine-fatty acid float in an alkaline circuit' The flotation mica is ground in a pebble mill to liberate contaminated particles and expose iron-stained surfaces. The concentrate is then conditioned with a reducing reagent and passed through a Frantz Ferrofilter for the removal of magnetic material. The nonmagnetic product is acid-leached to further reduce the iron contamination. The leached material is then ground to product size and checked for specifications. Discussion of Flowsheet and Its Development Flotation: The mica is concentrated by flotation using either of two methods: an amine float in an acid circuit, or the USBM amine-fatty acid float in an alkaline circuit.' The writer does not wish to elaborate on the relative merits of the two flotation procedures; however, after investigating ores from various localities, it was felt that the amine float in an acid circuit was more consistent when using a standard procedure. Iron Reduction and Magnetic Separation: This phase of the flowsheet may be considered the key to the process. It was developed when it became apparent that hot acid leaching had just about reached its maximum effective¬ness for color improvement. S. C. Sun and William Hirsh in "Hydrochloric Acid Leaching of Iron From Pennsylvania Aluminous Clays,"` reported iron extraction for -20-mesh minerals of 82.20% for hematite and only 33.46% for magnetite. The test procedure employs a 19% hydrochloric acid solution, 104°C temperature, and 10% solids. It was felt that magnetic removal of some of the iron contaminants before leaching therefore would be advantageous. Since the iron minerals present as contaminants vary in magnetic susceptibility, it was decided to experiment with reducing reagents as a means of influencing their magnetic properties. Zinc hydrosulfite or sodium hydrosulfite were the reagents used throughout the tests. The procedure developed consists of conditioning the mica concentrate in a closed container at approximately 5% solids with a reducing reagent. Approximately 10.0 lb of zinc hydrosulfite per ton of mica concentrate is used. The ferric iron is reduced to the ferrous state (magnetic) and oxidation is retarded by isolation from the atmosphere in a closed container. After the desired conditioning time has elapsed, the material is passed through a Frantz Ferrofilter where the magnetic particles are retained and discarded as

Jan 1, 1972

By R. A. Morse, C. R. Holmgren

The production of oil by water flooding can be substantially increased by the maintenance of free gas saturation in the reservoir during the flooding operation. This effect is accomplished by the alteration of oil relative permeability characteristics and the occupation by gas of pore space that would otherwise be filled with residual oil. The amount of reduction In residual oil can be calculated from appropriate water-oil relative permeability characteristics. This paper presents experimental data in support of the foregoing conclusions and an example of the calculations. The microscopic pore saturation concepts of the mechanism are discussed. A method of practical application to field floods is presented together with discussion of certain limitations. INTRODUCTION The presence of free gas has been reported by a number of investigators to significantly affect the oil recovery which can be obtained from sandstone flow systems by water flooding.1,2,3,4,5 The effect of gas, noted in every instance, has been to cause lower residual oil saturations than could be obtained by water flooding the same systems in the absence of free gas. The degree of improvement in recovery has been observed to vary widely, depending on the systems used and the conditions of the tests. The increased oil recovery obtained because of the presence of gas during a water flood has been variously attributed to changes in physical characteristics of the oil, selective plugging action of the gas, inclusion of oil mist in the free gas phase, and the additional sweeping or driving action of the free gas. All but the first of these suggestions imply changes in the displacement mechanism. The change in viscosity and inter-facial tension of the oil phase, within the pressure range used for all the experimental work, is certainly not sufficient to account for the differences in residual oil saturation noted dess there is a drastic change in the displacement process. One other effect which logically seems capable of causing differences in residual oil saturation of the magnitude noted in the experimental work is that of simple replacement. In a water-wet system containing oil, water, and gas, it is to be expected that the gas will exist inside the oil. This is the position of minimum free surface energy, since the gas-oil interfacial tension will be less than the gas-water interfacial tension. There is no apparent reason to expect that the existence of free gas within the oil phase should alter the saturation at which the non-wetting phase (now oil and gas) should become discontinuous and hence trapped so as to be unrecoverable by direct displacement by water. If this is the situation, then trapping of a certain percentage of gas saturation during water flood should result, at infinite water-oil ratio, in a like reduction of oil saturation below that attainable by flooding in the absence of free gas. It is visualized that the gas will exist as bubbles inside the discontinuous residual oil as illustrated in Fig. 1, with the size of the oil bubbles being substantially unchanged due to the presence of the gas. As a practical matter, it can be anticipated that the presence of a free gas saturation inside the oil phase will reduce the relative permeability to oil which will exist at any particular water saturation. This reduction will be caused by two factors — -the addition of the gas-oil interface, and the reduction of area available for oil flow in the pores containing gas. This reduction in oil permeability at any particular water saturation will result in water breakthrough at a lower water satu-

Jan 1, 1951

By M. van Swaay, C. E. Birchenall

A. S. lBrling (Johnson, Matthey & Co. Ltd., Laboratories)— Because of its initial emphasis upon the production of membranes by powder metallurgy and by normal casting and rolling techniques, this paper is of considerable interest to those concerned with the design and application of industrial diffusion units. As the subsequent performances of the three grades of palladium were not compared in detail, it can perhaps be concluded that specific permeabilities were very similar. A short table in which the absolute diffusion rates obtained by the Authors with their three grades of palladium were compared with those of other investigators18-25 would, however, have been of particular value in view of the very variable results obtained with this membrane material. The Authors state early in the paper that the permeability of the thoriated samples decreased in much the same way as that of the nonthoriated samples. Later, they emphasize that sintered palladium without thoria was very susceptible to poisoning and did not respond as well as the other materials to the decontamination treatment. These results, which suggest that sintered membranes are less effective than those produced by normal methods, have considerable industrial implications, and amplified experimental details would be very welcome. The results of the Authors differ in one important respect from those of other investigators. Thus, no "bleeding" was required to maintain a constant rate of diffusion, although Davis, Darling, and De Rosset found this to be necessary. Davisz3 found that "bleeding" was necessary even when the hydrogen had already been purified by diffusion through a palladium tube, and it would be of value to know the purity of the gas stream emerging from the purification train shown on Fig. 1 of the Authors' paper. The Authors show that diffusion rates varied according to the square root of the pressure difference providing that the palladium membrane was not contaminated in any way. This finding is in agreement with the results of those workers2,18,19,23,24 who diffused their hydrogen through the membrane into a partly or completely evacuated system.Those workers who did not find this square-root law dependence did not have a partial vacuum on the downstream-side of the membrane.'' The difference between the two experimental con- ditions appears to be significant, as the absolute diffusion rates per sq cm of diffusion area are lower when the square-root law dependence is not observed. While the effect could be attributable to poisoning of the membrane it is conceivable that a positive pressure on the downstream side of the membrane might influence the mechanism of diffusion. In view of the capricious behavior of palladium, factors of this type, which are likely to influence its industrial significance, should be carefully considered. M. van Swaay and C. E. Birchenall (authors' 9-eply) — Our Eq. 191 may be compared with several others for hydrogen diffusion in palladium: D = 5.95 x exp (-5770/1)' D = 3.55 x 10"3 exp (-5550/)' D = 4.3 x 10"3 exp (-5620/1')' Other values are given in the literature, but in many cases surface poisoning may have affected the results. No analysis was performed on the hydrogen gas stream after purification. It should be noted that the presence of liquid nitrogen traps near the heated membranes may have set up convection currents and provided the equivalent of "bleeding." The question of whether a partial vacuum existed on the downstream side of the membrane is not clear. Does it imply that there is something unique about the permeability measured near one atmosphere pressure? We do not believe that there is. When square root pressure dependence is not observed, diffusion rates are not measured. Is "absolute diffusion rate per sq. cm. of diffusion area" equivalent to permeability? The terminology is new to us. In any case we doubt very much that pressures of the order of several atmospheres can affect the mechanism of diffusion.

Jan 1, 1962

By G. W. Doepel, W. P. Sibley

A three-dimensional, multilayer, mathematical model has been developed for predicting performance by mis-c.ible displacement. Areal and vertical coverage, total gas-oil ratio, propane-oil ratio and oil-recovery factor for the multilayer system are computed as functions of hydrocarbon pore volumes of solvent injected. The technique has been programed for use on the IBM-7090 digital computer and was used in calculating perfortnance for continuous miscible displacement of a five-spot, 50-layer .system. Results reported are based on the use of published laboratory data for a single, uniform five-spot layer. However, the technique is sufficiently generalized to calculate performance of any system of layers of any pattern configuration and mobility ratio, provided data are available to describe perfortilance of a single layer. The technique has also been adapted for predicting oil recovery resulting from various gas-driven slug sizes. The results presented for both continuous and slug-type injection show, that permeability variation has 0 significant effect on oil recovery but that solvent composition, within the range studied, has a lesser effect on oil recovery. INTRODUCTION Experience from many of the miscible-phase displacement projects currently in operation indicates that the injected solvent flows preferentially through the more permeable zones in the reservoir. This preferential solvent flow in the high-permeability streaks results in a reduction of vertical coverage (C,) and an attendant decrease in oil recovery. The multilayer technique presented in this paper uses areal coverage (C.,) data for a single layer from the literature and from laboratory measurements of displacement efficiency (herein referred to as microscopic coverage, C,,) to calculate areal coverage (C,) and vertical coverage (C,) for the multilayer system. Oil recovery for the multilayer system is obtained as the product of (Ca) (Cm) (Cv). It is assumed that C, in a single layer is 100 per cent; therefore, (C,) (Cm) is a measure of oil recovery from that layer. Data presented in the literature for a single layer generally deal only with areal coverage and/or microscopic coverage —but seldom specifically with oil recovery. Nevertheless, this technique is flexible in that either Ca and C,, of an individual layer or a direct measurement of oil recovery for the layer may be used. In developing this technique, the problem arose concerning a method of determining the relative amounts of solvent that will enter each layer during any increment of solvent injected. To solve the problem, a dimensionless equation was derived for expressing solvent injected into each layer as a fraction of total injection. A finite-difference method was applied in integrating the equation to determine oil recovery for the system of layers as a function of pore volumes injected. The technique subsequently was programed for solution on a digital computer. The method was applied to predict performance of a deep California reservoir under continuous miscible displacement and has been extended for calculating oil recovery as a function of slug size for gas-driven slug injection. The technique may be applied to any petroleum reservoir that can be adequately represented as a layered system, by considering permeability variation,1-3 permeability range and average permeability. Each layer is then assigned its appropriate permeability, porosity and hydrocarbon saturation. RESERVOIR VOLUME CONTACTED BY SOLVENT It has been recognized for some time that high microscopic coverage (C,) can be obtained by injecting a miscible agent, such as propane, into an oil reservoir. These high values of C,,, have been demonstrated in the laboratory and are normally between 80 and 100 per cent, depending on solvent and reservoir oil composition as well as on reservoir temperature and pressure. In miscible-phase displacement, because there is no interface between solvent and oil, effective permeability to solvent is the same as that to oil. (This assumes that oil saturation has not changed and that water is immobile.) Consequently, mobility ratio (M) simply reduces to that of viscosity ratio . Perrine4 and Habermann5 have shown that large

By H. Schwartzbart, J. R. Low

Annealed, .Poly crystalline, low carbon steel exhibits a phenomenon known as the "yield point." If such a steel is loaded in tension, the load increases steadily with elastic strain, drops suddenly, fluctuates about some constant load value for a time and then rises with further strain. The load at which the sudden drop occurs is called the "upper yield point"; the load at which deformation proceeds at constant load is known as the "lower yield point"; the amount of elongation at constant load is termed the "yield point elongation." Deformation which occurs during the yield point elongation is heterogeneous in nature, deformed metal existing next to undeformed metal. The drop in load from the upper yield point is concurrent with the formation of a band of deformed metal which usually appears at a stress concentration such as a fillet and at an angle of approximately 45" with the direction of pull. As deformation proceeds during the yield point elongation, the deformed region grows until the entire specimen has been strained an amount equal to the strain in the first band. A phenomenon very often associated with the yield point is that of strain aging. Aging in steels may be of two types: quench aging, or precipitation hardening, and strain-aging. Strain aging differs from precipitation hardening in that if the metal is merely strained instead of given a solution anneal and quenched, it then undergoes changes in properties similar to those observed in quench aging. The yield point effect and strain aging are most markedly exhibited by low carbon steel, but have been observed in other alloys. Numerous explanations for the mechanisms of both phenomena have been advanced, but the evidence for the acceptance of any of the hypotheses is incomplete. One of the explanations of the yield point phenomenon frequently proposed is that there exists a honeycomb of grain boundary material that supports the load above the yield strength of the ferrite and that when this network ruptures, the ferrite flows with no increase in load. This notion was first advanced by Dalbyl in 1913 and has been since put forward in more concrete form by many investigators. The testing of this hypothesis was, in part, the object of the present investigation. Similarly, the mechanism by which strain aging occurs remains to be explained. Presumably, it is a precipitation hardening phenomenon similar to quench aging, the precipitation being initiated by strain. However, it is not known whether the process is a solution and precipitation effect or a change in solubility due to strain. One obvious method of determining whether or not the yield point is a grain boundary phenomenon is to determine if single crystals exhibit a yield point. Low and Gensamer2 have shown that carbon and nitrogen are responsible for the yield point and aging in steel and that these effects can be eliminated by the removal of these elements through wet hydrogen purification. Single crystals that have been grown in wet hydrogen treated material must then be recarburized or re-nitrided in order to evaluate the effect of grain boundaries. In the present investigation, single crystals of iron which had been annealed in wet hydrogen to remove the carbon and nitrogen were cut into two specimens having the same orientation with respect to the tension axis. Of these, one specimen was pulled in tension in the wet hydrogen treated state and the other was carburized or nitrided and tested similarly. All the crystals were examined for a yield-point and strain aging. The grain boundary hypothesis for the explanation of the yield point in steel has been tested only indirectly previously. Arrowsmith,3 Andrew and Lee,4 Edwards and Pfeil,5 Fell,6 Ludwik and Scheu,7 Winlock and Leiter8 and others found a decrease in proportional limit and yield point effect with an in-

Jan 1, 1950

By J. B. Cohen

EARLY investigators1-4 reported that additions of beryllium to silver caused embrittlement. The in-termetallic phases which are now known to exist5,' were not known at the time. Solidification appeared to occur by a simple eutectic reation.4,7 (The primary beryllium-rich phase appeared globular, when the liquidus temperature was near the melting point of pure beryllium, and angular when the liquidus was much lower; in discussion of Sloman's work,4 it was suggested that this difference was due to an intermetallic compound which was not revealed by the polishing technique.) In the investigation reported here, it was found that if the compounds now known to exist were avoided, certain beryllium-silver alloys containing large amounts of the primary beryllium phase possessed sufficient ductility to be cold formed. Alloys containing 14 and 20 at. pct Ag were prepared from vacuum-cast beryllium containing 1.0 to 1.2 pct impurities and silver, 99.99+ pct pure. The elements were melted together by induction in beryllia or alumina crucibles under argon at less than 1 atm pressure. The cooling rate was adjusted by regulating the temperature of the melt and by reducing the power to zero without stopping the water flow in the induction coils, which were in contact with the crucible. Ingots were about 0.75 in. in diam with a grain size of 0.12 to 0.25 in. Microscopic investigation of transverse or longitudinal sections of an ingot cast from a temperature well above its melting point showed that, in general, the compounds did not form, presumably because of rapid cooling; the alloys consisted of beryllium-rich particles in a silver-rich matrix as shown in Fig. 1. A eutectic structure was not observed, probably .due to the low beryllium concentration of the eutectic composition and eutectic divorcement. Within a given grain, the small globules of the primary phase were all oriented in the same direction, as indicated by examination under polarized light; the orientation varied from grain to grain. This suggested that the globules were branches of the same dendrites. If cooling was too slow, a third constituent was observed in the microstructure; because it formed around the primary beryllium-rich globules it was Fig. 1—Ag-Be alloy, 18.6 at. pct Ag, as-cast. Bright field illumination. Etched with an aqueous solution of NH,OH and H2O2. X500. Reduced approximately 21 pct for reproduction. probably the 6 compound. A thin ring around an ingot often contained this constituent. When this compound was present in large quantities, or the silver segregated so that large regions occurred with the beryllium-rich globules in contact, sections from the ingot were brittle. When the compound was initially absent, it was possible to carry out reductions in thickness of 50 pct by cold rolling with only limited edge cracking; failure did not occur until reductions of 80 pct. A typical cross section after cold rolling is shown in Fig. 2. After simple bend tests, microcracks were found; these apparently were difficult to propagate in the two-phase structure, as the specimens continued to behave in a ductile manner after the cracks appeared. Similar cracks would probably be closed during rolling. It is apparent in Fig. 2 however, that the primary beryllium particles undergo extensive plastic deformation. Hardness values of the alloys indicated that some work-hardening occurred during cold working. Fig. 2—Ag-Be alloy, 18.6 at. pct Ag, reduced 83 pct in thickness by cold rolling. Bright field illumination. Etched with an aqueous solution of NH,OH and H2O2. X500. Reduced approximately 21 pct for reproduction.

Jan 1, 1961