Search Documents

By D. M. Buck

greatly restricted and every effort is being made to do away with it where possible, and to substitute a lead-base babbitt or a babbitt with 50 per cent,. tin. The amount of solder having the composition 50 lead, 50 tin, used by the Pennsylvania Lines East, during 1917, was approximately 100,000 lb., but there is reason to believe that a large portion of this can be replaced by a 60-lead, 40-tin solder with satisfactory results, and instructions have been issued to this effect. Inregard to the quantity of new tin used, it is not possible to give the amount, except approximately and from calculations based on the 1917 consumption of bearing metals by the Pennsylvania Lines East only, it is estimated that about 770,000 Ib. of new tin were required in a total of about 21,380,000 lb. of all kinds of bearing metals turned out by the foundry or purchased in the market. No change has been made in the specifications for bearing metals for some years, as the metals used have been satisfactory. A large proportion of the bearing metals are made up from old material re-melted and brought to standard composition by some addition of new metal, and every effort is being made to utilize old material to the best advantage and reduce the amount of new metal of all kinds purchased. For a number of years no tin has been used in the lining metal of either passenger or freight car journal-bearings, except such small amounts as come in from re-melting old linings. No change has been made in phosphor bronze used for rod bushings, as we should expect some trouble from bushings pounding out of shape if a phosphor bronze were used which contained less tin or more lead than the present specifications call for. In this, as well as in the case of all other bearing metals, we expect to use our utmost endeavors to economize and to substitute for tin wherever possible. The Tin-plate Industry D. M. Buck,* Pittsburg, Pa.—During the first 5 months of 1918, approximately 11,000,000 Ib. per month of pig tin were consumed in the United States. Solder, bearing metals, bronzes, etc, used about 5,500,000 lb.; collapsible tubes a little more than 250,000 lb.; tin-foil about 500,000 lb.; and the tin- and terne-plate industry somewhat less than 5,000,000 Ib. In an effort to reduce this consumption and thus conserve our tin supplies, several methods of procedure suggest themselves: 1. Salvage. The most careful and systematic collection and re-use of tin and tin-bearing materials is economically important, in that we thus

Jan 1, 1919

By A. Vifian, I. lwasaki

Several samples of copper-nickel ore from the Duluth Gabbro were studied to relate their mineral-ogical characteristics with their amenability to concentration by flotation. The most common rocks in the samples were norites, gabbros, and troctolites. The texture of the rocks varied from coarse pegma-titic to very fine. The copper contents ranged from 0.2-0.7% and the nickel contents from 0.07-0.15%. The ore minerals were chalcopyrite, cubanite, pyrrhotite, and pentlandite. Most of the ore minerals occurred in interstices between laths of plagioclase, but some were included in the gangue grains. The samples responded readily to xanthate flotation. After a 65 to 100-mesh grind about 95% of the copper and 85% of the nickel were recovered in a bulk sulfide concentrate. Such recoveries remained constant even with grinds as fine as 270 mesh due to the extremely fine inclusions of sulfides in the silicate gangue. A finer grind than 270 mesh was necessary to free the copper-iron sulfides from the nickel-iron sulfides. A combination of time and British gum* was quite effective in producing concentrates analyzing 25-30% copper. The nickel concentrate averaged 1.5% nickel. Most of the nickel lost in the tailing was present as un-liberated sulfides and not as nickel in the olivine structure. The nickel content of the olivine was only about 0.05 or 0.08%. Western world nickel usage in 1967 is expected to be the second highest on record, amounting to 825 million lb. Despite a rise in mine production, however, demand is expected to exceed supply in 1968. ' Hence, exploration, development, and production expansion projects are being accelerated on a world scale. In Minnesota the copper-nickel bearing Duluth Gabbro is receiving much attention. During the past year 11 companies have been granted leases in the region, and the 1967 Minnesota legislature enacted a tax law designed to encourage the development of a copper-nickel industry. Although it is only during the past few years that mining companies have shown an active interest in the Minnesota copper-nickel deposit, copper-nickel minerals have been known to exist in northeastern Minnesota since 1899.3 At the 28th Mining Symposium, Dr. P. K. Sims, from the Geological Survey of Minnesota, gave a review of copper-nickel exploration in northeastern Minnesota,4 in which he emphasized that the potential of the Duluth Gabbro complex is large. Tests on samples from the Kawishiwi River area were reported by the United States Bureau of Mines in 1955.' In 1964, the Mines Experiment Station described tests on samples from the Gunflint Trail: and in 1966 the Minnesota Geological Survey published a geological map of the Gabbro Lake Quadrangle. The Mines Experiment Station published a study of the concentration of a Minnesota copper-nickel ore in 1967.7 The present paper is divided into three parts. The first describes the mineralogy of some rock samples from the Kawishiwi River area of northeastern Minnesota; the second summarizes the results of numerous flotation tests that were made on the samples, and the third explains the loss of nickel in the flotation tailing. MINERALOGY OF KAWlSHlWl ROCKS Recently five samples of mineralized gabbro from the Kawishiwi River area, a few miles east of Ely, were received at the Mines Experiment Station. One of the samples was a 175-lb composite, which was designated Ore 1556. The other four samples were 7-lb grab samples that had been obtained from four different exploration pits. They were designated Ore 1626, Lots 1 through 4. Basic Mineralogy: Many authors have reported that the rocks of the so-called gabbro of the Duluth Intrusion vary widely in composition.8-9 The samples on which our studies were based also contained a wide range of rocks. A gabbro, olivine-gabbro, anorthosite-gabbro, norite, and a troctolite were represented. The essential minerals in the samples were plagio-clase, pyroxene, and olivine, in variable proportions.

Jan 1, 1969

By H. A. Lesser, H. L. Stone, G. H. Bruce

ESSO PRODUCTION RESEARCH CO., HOUSTON, TEX. A mathematical model that represents the conduction heating of a rock formation of limited permeability is formulated. Heat is introduced by the injection of a hot condensing gas into horizontal fractures through the formation. The fractures are assumed to be equally spaced; the flow of fluid is linear. The model consists of three coupled partial differential equations and two functional relationships. Temperature histories for both fractures and formation are obtained by numerical procedures. Application of the model in predicting in situ heating of oil shale is illustrated with four sample cases. Influences of shale thermal diffusivity, fluid pressure in the fracture and fracture spacing upon heating rate are demonstrated. In all cases considered, several years were required to heat all the formation to 700F with 1,000F steam. Less than 32 per cent of the injected heat was utilized for heating the shale. To facilitate extensive calculations of this type, a set of seven dimension-less quantities are derived. With this dimensionless system, the effect of controllable parameters on economics can be explored. INTRODUCTION Vast deposits of hydrocarbons are contained in oil shale formations in the western United States and in other countries. These hydrocarbons, known as kerogen, are difficult to recover because they are immobile solids and because the rock formations in which they are found have very limited permeability. The kerogen can be converted through pyrolysis to recoverable liquids and gases if the formation temperature is elevated to the range 600 to 800F. To accomplish this heating, the formation might be artificially fractured and a hot fluid injected into the fractures. If the injection and production wells are vertical and the fractures horizontal, interwell communication will be established and the hot fluid will flow from an injection well to a production well. In this process, some of the heat in the fluid will be given up and conducted into the formation. Products of pyrolysis will flow from the heated rock matrix into the fractures and be produced, commingled with the injected fluid. The practicality of such a process is strongly influenced by the rate at which heat can be transferred from the hot fluids to the formation. An analysis of this heat conduction problem is the subject of this paper. Linear flow of a condensing gas through a set of horizontal, equally spaced fractures is considered. Temperature histories for both fractures and formation are obtained from numerical procedures. The analysis also applies to formations containing streaks of natural permeability instead of fractures. Heat flow problems similar to the one treated here have been studied by other workers. Thomas1 has solved a radial problem in which he considers the injection of a noncondensing gas into a single, horizontal fracture and approximates the temperature profile in the fracture by a step function. The solutions to other problems, which are less similar to the present problem than the one treated by Thomas, have been reviewed by Spillette.2 In addition to these mathematical studies, a field test of an in situ shale oil recovery process using hot natural gas as the injected fluid has been reported.3 The first part of this paper is concerned with the mathematical development of the problem. The differential equations, boundary conditions and functional relationships which describe the physical model are formulated and the equations are cast into finite difference form. The method of solution of the difference equations is discussed (Appendix B). The second part of this paper presents results computed for a typical oil shale formation for which steam at 1,000F is the injected fluid. Fracture spacing, horizontal and vertical thermal diffusivities

Jan 1, 1967

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967

By D. S. Rowley, F. C. Appl

Assuming that rock behavior, during cutting with a single diamond, may be approximated by that of a rigid, Coulomb, plastic material, a theory of single diamond cutting action has been developed. Using this theory, the stresses on the diamond cutting surface and the components of the cutting force have been determined. Theoretical results agree reasonably well with available experimental data. The theory of cutting of single diamonds may be used as the basis for subsequent study of the perfonnance of surface set diamond drill bits and other surface set diamond cutting tools. INTRODUCTION For more than twenty years, surface set diamond bits have been used successfully in the mining and petroleum industries to drill and core rock formations that range from medium hard to the hardest known formations. Most extremely hard rock formations cannot be economically drilled or cored by any other means. In view of the growing use and importance of surface set diamond drill bits, it is reasonable to study the way in which these tools cut, and the factors that influence wear of the diamonds. Such knowledge will contribute to better bit design and will improve operating procedures. To analyze the cutting performance of these tools it is appropriate to consider the way in which an individual diamond cuts. Once this is established, the over-all performance of the bits can be determined by considering the total effect of the individual diamonds. In recent years there has been a growing interest in rock mechanics and there now exist several publications dealing with the deformation and failure of rocks. Notable experimental work has been reported by Bredthauer,1 Handin and Hager,2, 3 and Gnirk and Cheatham.4 An extensive collection of data pertaining to rock properties at low confining pressures has been given by Wuerker.5 Cheatham and Gnirk have recently presented an excellent review of the state of knowledge regarding analytical and experimental rock deformation and failure.6 Generally rocks are extremely weak in uniaxial tension, and weak or only moderately strong in uniaxial compression. Both tensile failure and compressive failure are generally brittle in nature. This means that there is little or no deformation or "flow" of the material before shear or tensile fracture occurs. Thus, until fracture occurs, rock behaves very nearly like an elastic material in uniaxial stress. Significant changes in behavior of rock occur under conditions of triaxial compressive stress. There is generally a marked increase in compressive strength and in ductility. It has been demonstrated that the failure becomes plastic (in a macroscopic sense) at sufficiently high values of confining pressure. Sharp wedge chisel indentation studies conducted by Gnirk and Cheatham indicate that many rocks behave in a plastic manner at relatively low triaxial compressive stress states.4 If "failure" under triaxial compressive conditions is considered to be either complete rupture or macroscopic plastic flow, then the general failure characteristics of a rock can be reasonably well represented by the commonly known Mohr envelope. This certainly is not a perfect failure criterion (since the effect of the intermediate principal stress is neglected) but it has been found to be adequate for most engineering purposes. Mohr envelopes usually are constructed graphically using experimental triaxial test data. For purposes of analysis it is convenient to represent the envelope in equation form. Among the several approximate forms appearing in the literature, a parabolic form is suggested by Cheatham,7 and an exponential form is suggested by Paone and Tandanand.8 The exponential form has been adopted for this study, and based on published data, appropriate parameters have been established for 19 rock types. A typical plot of the data compared with the approximate envelope, and a summary of the rock properties and

Jan 1, 1969

By A. K. Snelgrove

This symposium is designed to marshal a state-of-the-art account of approaches and solutions to problems arising out of waterlogging and salinization in the world's largest irrigation system which was dismembered by legal beheading of the Eastern Rivers under the lndus Waters Treaty of 1960. In this introduction, the geological and historical backgrounds of surface waters and groundwaters in the Indus Basin are presented with a summary of remedial measures undertaken and proposed to combat waterlogging and salinization. At a time when water-resource planning is proceeding on a continental scale, a study devoted solely to the Indus River Basin may seem to require an apologia. None is offered, however, for three major reasons: (1) Although the Indus is a single river basin, it is one of great magnitude and importance. The Indus River, which comes from Himalayan sources to form part of the great Indo-Gangetic Plain drainage system, is vital to agriculture and hydroelectric power in the thickly populated northern subcontinent of Indo-Pakistan. (2) The Indus has become a plexus of extraordinary problems since the dismemberment of the world's largest irrigation system of which it was a part. The break-up of this system occurred with the legal beheading of the Eastern Waters under the Indus Waters Treaty of 1960. (3) Th River Basin has been the target of a Herculean effort to contain the twin encroachments of waterlogging and salinization that now threaten it with a shortage of fresh water. The British irrigation system of barrages and canal on the Indus was a venture admirable in conception. Unfortunately, it is now a legacy burdened with unwanted consequences due to interference with the hydrological regimen, particularly the rise of water table occasioned by leakage. To cope with this situation, there has been a possibly unprecedented concentration of domestic and international engineering and scientific talent on the Indus and its problems (about a score of organizations by recent count). To marshal a state-of-the-art account of their approaches and solutions is the purpose of this Symposium. The coordinated effort now being made on a variety of fronts under the general supervision of West Pakistan Water and Power Authority (WAPDA) is directed toward the modification and control of interactions and relationships between geological materials — soils, water and its salt balance, and air —in what has been termed a huge "underground lake" in an alluvial basin. The problem is an intriguing one in that a highly complicated case history of river-course changes visible in the present surface must be projected downwards through Recent and Pleistocene sediments in order to decipher the environmental controls of deposition in terms of vertical and lateral variations in lithology and patterns. Such studies must be extended to compaction, deformation, and fabrics of the alluvium. Geohydrology is often defined as groundwater hydrology. As employed herein, it includes the distribution, quality, and behavior of surface and subsurface waters, including their evolution, vagaries, and geological, climatological, and botanical milieux." GEOGRAPHY The Indus (Sanskrit, Sindhu = river) rises at nearly 17,000 ft above the sea in the glaciers on the northern slopes of Kailas Parbat, Tibet. It is about 1800 miles long and has a total drainage area estimated at 372,000 sq miles. The average discharge is 196,000 cfs at mouth. Leaving the Plateau of Tibet, the river flows at a level which is miles above the sea, and in places miles below Himalayan peaks, in a structurally controlled northwestward course. Joined in this course by six tributaries and swollen by the melt waters of giant glaciers, it swirls about the base of the famous

Jan 1, 1970

By H. W. Mead, C. E. Birchenall

The rates of growth of graphite nodules in cast irons are calculated for a model of a growing graphite sphere surrounded by a shell of austenite through which carbon and iron are diffusing. The carbon is supplied by a cementite-austenite mixture at the outer edge of the shell. The calculated rates from carbon diffusion agree fairly well with the measured values, but the rates based on iron self-diffusion are much too low. Because the iron must be removed for the process to continue, o plastic deformation mechanism is suggested. THERE has been an unfortunate tendency in the literature of growth and dissolution reactions to identify the slope of a plot of the logarithm of the growth or dissolution rate vs reciprocal of temperature with the activation energies determined by diffusion measurements in chemically similar systems. If the solubility limits at the growing or dissolving phase interfaces change appreciably with temperature, such procedures may lead to serious errors in specifying the mechanisms controlling the rates of these processes. A reaction involving the simultaneous growth of graphite from and dissolution of cementite in austenite will be discussed here to show that substantial oversimplifications have been made in earlier studies of this process. From a dilatometric study of the overall y-range graphitization process, Bunin and Danil'chenko' found an apparent activation energy of 90 kcal per g-atom. They concluded that, since this value was much nearer the approximately 68 kcal per g-atom observed for self-diffusion in pure iron than it was to 32 kcal per g-atom for diffusion of carbon in iron, the rate-controlling process must be the diffusion of iron away from the growth interface to provide space for graphite to grow. Bunin and Salli' bolstered this agreement with the observation that graphite forms preferentially at surfaces, cracks, and pores where free volume is present initially. Ziegler, Meinhart, and DeaconQ ad found an apparent activation energy for graphitization of 98 kcal per g-atom. This value, like that of Bunin and Danil'chenko, applies to nucleation and growth occurring simultaneously. It is not generally assumed that nucleation is entirely diffusion-controlled. Treatments such as prior cold-working and prequenching, which strongly affect the rate of graphitization, would be expected to have little effect on diffusion processes occurring at the graphitization temperature. In any case 90 kcal per g-atom is not particularly close to the 54 kcal or lower to be expected for diffusion of iron which is carbon-saturated.' The presence of silicon in iron is said to reduce the activation energy for diffusion,' so its presence in cast irons would be expected to have a similar effect. The growth rates of nearly spherical graphite nodules have been studied by Brown and Hawkes." The initial rate of growth before impingement of adjacent carbon-depleted areas was found to be parabolic, and rate constants were measured at several temperatures. Their apparent activation energy for the growth process only was roughly 45 kcal per mol. Burke and Owen' made growth-rate measurements on several alloys, recording appreciable induction periods, initial nucleation rates, and average radial growth rates based on an assumed linear rate of radial growth. They found that increased silicon content increases the ovel-all reaction rate and the initial rate of nucleation, and decreases the induction period. Assuming 1) spherical nodules. 2) constant volume rate of nucleation, 3) austenite in equilibrium with cementite at positions where the cementite is not depleted, and 4) a constant diffusion distance during the growth of a nodule, they calculated an apparent overall activation energy of 68 kcal per mol, which was related to an activation energy for nucleation of 79 kcal per mol and activation energy

Jan 1, 1957

By T. K. Perkins, O. C. Johnston

Because of the influence of dispersion on miscible-displacement processes, diffusion and dispersion phenomena in parous rocks are of current interest in the oil industry. This paper reviews and summarizes a great deal of pertinent information from the literature. Porous media (both unconsolidated packs and consolidated rocks) can be visualized as a network of flow chambers, having random size and flow conductivity, connected together by openings of smaller size. In such a porous medium, the apparent diffusion coefficient D is less than the molecular diffusion coefficient Do, as measured in the absence of a porous medium. For packs of unconsolidated granular material the ratio D/D, is about 0.6 to 0.7. For all porous rocks, both cemented and unconsolidated, the ratio of diffusion coefficients can also be represented as D/Do = 1/FF where F is the for- mation electrical resistivity' factor and 4 is the porosity. If fluids are flowing through the porous medium, dispersion may be greater than that due to diffusion alone. At moderate flow rates the porous medium will create a slightly asymmetrical mix zone (trailing edge stretched out), with the longitudinal dispersion coefficient approximately proportional to the first power of average fluid velocity (if composition is nearly equalized in pore spaces by diffusion). If the velocity in interstices is large enough, there will be insufficient time for diffusion to equalize concentration within pore spaces. In this region, longitudinal dispersion increases more rapidly than fluid velocity. At low velocities in interstices, transverse dispersion is characterized by a region in which transverse diffusion dominates. If the fluid velocity gets high enough, there will be a transition into a region where there is stream splitting with mass transfer but with insufficient residence time to completely damp-out concentration variations within pore spaces. There are several variables that must be con- trolled to get consistent longitudinal and transverse dispersion results, viz., (1) edge effect in packed tubes, (2) particle size distribution, (3) particle shape, (4) packing or permeability heterogeneities, (5) viscosity ratios, (6) gravity forces, (7) amount of turbulence, and (8) effect of an immobile phase. INTRODUCTION Diffusion and dispersion in porous rocks are of current interest to the oil industry. This interest arises because of the influence of dispersion on miscible-displacement processes. In a recovery process utilizing a zone of miscible fluid, there is the possibility of losing miscibility by dissipating the miscible fluid or by channeling or "fingering" through the miscible zone. Diffusion and dispersion are two of the mechanisms that may lead to mixing and dissipation of the slug. On the other hand, dispersion may tend to damp-out viscous fingers which may be channeling through the miscible slug. 58 Hence, dispersion may be detrimental or beneficial (if it prevents fingering through the miscible zone). Therefore, it is doubly important that we understand these processes. In this paper we review, summarize and interpret a great deal of information from the literature. In particular, we will briefly discuss molecular diffusion in miscible fluids. Then we will discuss what differences to expect for diffusion in a porous rock. If there is movement of the fluid through the rock, then there may be an additional mixing or "dispersion". Furthermore, the dispersion longitudinally (in the direction of gross fluid movement) and transverse to the direction of fluid movement will not be equal. We will discuss both types of dispersion as well as several variables which can affect dispersion (viscosity differences, density differences, turbulence, heterogeneity of media, etc.). This group of variables has sometimes led to difficulty when comparing literature data. DIFFUSION OF MISCIBLE FLUIDS If two miscible fluids are in contact, with an initially sharp interface, they will slowly diffuse into one another. As time passes, the sharp interface between the two fluids will become a diffuse

By H. A. Wriedt

Zinc-vapor pressures in equilibrium with bcc solutions of zinc and iron at 703.5", 757", 793", and 900°C were measured by an isopiestic method. Thermody-na?nic properties of zinc in these solutions were derived directly and those of iron with a Gibbs-Duhem integration. Zirzc and iron exhibit large, positive, and moderately temperature-sensitive deviations from ideality. Excess enthalpies and Gibbs energies are positive in the ranges studied. It has been establzshed in relating published phase equilibria to the Gibbs energies of forming solid solutions that positive departures from ideality also exist in the fcc solutions. Tlze stability, as defined by Darken, of the bcc solutions is small but positive at 703.5"C. Its decrease with decreasing temperature indicates that a miscibility gap may exist in the stable or metastable bcc solutions at lower temperatures. The changes of thertnodynanlic properties relatable to changes in magnetic state are minor compared to other changes due to alloying. The accuracy with which the relation between lattice parameter and composition of the bcc solutions is known has been improved. Zinc increases the lattice parameter of a iron by O.ZOZA per unit atom action zinc to at least NZn = 0.25. TO our knowledge, only Ciganl has previously investigated thermodynamic properties of the bcc solutions of the Fe-Zn system. He measured zinc-vapor pressures in the temperature range from about 320" to 360°C by an effusion technique. We have measured the vapor pressures of zinc in equilibrium with bcc solutions of zinc and iron at 703.5", 757", 793", and 900°C by an isopiestic technique. This investigation and the pursuant thermodynamic computations are part of a program for studying several properties of these solid solutions. In the appendix we report a remeasurement of the relation at room temperature between the lattice parameter and zinc content. Related work not reported herein embraces measuring the diffusivity2 of zinc in ferrite at 700" to 900°C and the Curie point temperature3 as functions of composition. EXPERIMENTAL METHOD AND APPARATUS Purified iron at a series of constant temperatures was saturated with zinc at a series of constant zinc-vapor pressures, the range of conditions being so limited that the equilibrated solid solutions were ferri-tic. The equilibration was conducted in long capsules, Fig. 1, that had been evacuated, gettered with titanium, and sealed. The purified iron, part of a vacuum-melted ingot of electrolytic-iron stock (Glidden Co. 1044-grade), had the impurity contents listed in Table I. It was used in the form of millings about 0.001 in. thick. A zinc of 99.999 pct purity (United Mineral and Chemical Corp.) and later one of 99.9999 pct purity (Cominco Products, Inc.) were used, with no observed difference in results. Sufficient zinc was placed in the capsule (separate from the iron and kept at a different temperature) to supply the vapor phase and the iron alloy, while leaving a pool of pure liquid zinc throughout the run. The imposed temperature of the residual zinc pool fixed the vapor pressure of zinc throughout the capsule; condensation was avoided at other locations by making this the lowest temperature in the capsule. Two common complications of nonisothermal systems were, we believe, effectively absent from ours: first, molecular species other than monatomic have no appreciable concentration in zinc vapor at 700" to 900~~: second, the development of zinc-pressure differentials in the capsules was suppressed by use of relatively large-diameter tubing and the absence of diluents in the zinc vapor other than desorbed impurities. The specimens and zinc pool were located at two isothermal zones (each ranging within 1°C over several inches) of a wire-wound, 37-in.-long vertical-tube furnace. The required temperature distribution was achieved by external adjustment of the currents in three electrically isolated windings. The temperatures of the isothermal zones were controlled independently by ther-

Jan 1, 1968

By H. Bartholomew, L. Bernstein

Experirmental bonds of Ag-In SLID to gold, copper, nickel, Kovar, Dumet, nickel-plated molybdenum, ALL-Pt vrzetallizing on ceramic, and nickel-plated Mo-Ti metallizing on Al2O3 ceramic have been made and exalnined. Metallographic taper sections were made of all bonds for Phase identification. In addition, failuve temperatltres were determined for bonds heated undev shear load. Sowe bonds were subjected to electvon-probe microanalysis. The ternary SLID system Ag-In-Sn was similarly evaluated. Plwtonticrogaphs and electron-probe data are discussed and related to over-all phase diagrams , diffusion effects, and wetting properties of the liquid. The rates of diffusion andor compound formation of gold, silver, and copper in indium and In-Sn liquids are found to descend in that order, as measured by the relatzve thickness of phase layers in the interfaces of bonds made at 200to 400°C. It is concluded that the SLID process enables one to produce bonds which are stable at high temperatures, after having been fabricated at low temperatures, and meet most requirements for application in interconnection, silicon bonding, or hermetic sealing of electronic devices. In a previous publication by one of the authors, the theory and principle of SLID bonding in simple binary systems was discussed. SLID bonding can be defined as a joining process using a preform with a laminated structure: a low-melting material clad on a high-melting core. During the bonding cycle a metallurgical transformation occurs such that the resultant structure can withstand, without reliquefying, temperatures higher than those at which the bond was made. The phases present in the joint interface region are high-melting compounds or, if equilibrium is reached, solid solution. Properties will depend on the phases present, ultimately approaching those of the solid substrate core of the SLID preform. The high-temperature strength of bonds made in these systems was correlated to phase identification in the interfacial bond regions from metallographic cross sections and electron-probe X-ray microanalysis data. When complex electronic components are assembled, multiphase systems result at the interfaces of bond regions. In these ternary and higher-order systems very little diffusion and phase-diagram data is available with which to predict the nature of the interfacia1 region that will result from conventional fabrication processes and subsequent field use of the component. Clarke and Rhines did early work in the field using the Al-Mg-Zn system as a vehicle. Guy and co-workers studied a portion of the Cu-Zn-Ni system and some mathematical4 and general5 aspects of multiphase diffusion. irkald' and coworkerse recently studied diffusion in the Cu-Zn-Sn system. No diffusion studies have been made on ternary systems directly applicable to this investigation but some phase diagrams are available: g-u-n, u-i-n, and CU-i-n. I) DESCRIPTION AND EXPERIMENTAL PROCEDURE The elements Au, Ag, Cu, Ni, In, and Sn are found in alloy systems commonly used in fabricating and assembling electronic components. Reference to the appropriate phase diagrams9 and the previous publication will indicate that the binary systems of any of the two low-melting elements, indium and tin, with any one of the other four higher-melting elements are capable of producing SLID bonds. It is not unreasonable to assume that any ternary or higher-order alloy system containing at least one of the elements gold, silver, copper, or nickel and one or both of the elements indium and tin could be used to produce bonds which are stable at high temperatures, after having been fabricated at low temperatures. Indeed, the liquidus contours and isothermal sections of those phase diagrams available7 do so indicate. The parameters affecting the bond and the equipment used in evaluating the bond were described in an earlier paper.' The current investigation used the

Jan 1, 1967

By J. E. Fox Jr., J. L. Lummus, J. P. Gallus

Fresh water or brine drilling fluids may be kept free of suspended drilled solids by the addition of a water soluble acrylamide-carboxylic acid copolymer at the flowline. Addition of from .01 to 0.2 lblbbl of the polymer solution to drilling fluids containing less than 5 per cent clay solids by weight causes the solids to flocculate and settle rapidly to the bottom of the pits. Solutions of the polymer injected at the pump suction reduce the loss of water to permeable formations without impairing the permeability to oil. The use of drill-itlg fluids clarified with the polymer has resulted in increased drilling rates, extended bit life, and lower drilling fluid costs. INTRODUCTION Field experience and the results of a number of investigations have established that the accumulation of solids in drilling fluid reduces drilling rate.'-' Drilling with water makes possible higher drilling rates than can be achieved with mud but in many areas drilling with water may not be feasible because of the possibility of encountering high formation pressures or by-dratable formations. To combat these conditions, the drilling fluid must possess certain physical properties such as viscosity, density, fluid loss, etc., which require the presence of solids. The current trend in drilling fluid treatment is to maintain as low a percentage of solids as is possible without harmfully affecting the mud properties needed for drilling. Mechanical separating devices such as wet cyclones and decanting-type centrifuges are proving to be an economical means of controlling the percentage of solids and thereby maintaining low solids muds. Also: special muds have been developed for the purpose of retarding the hydration of drilled solids so that they may be more easily removed by shale shakers or settling pits. In some areas, such as West Texas, the characteristics of some formations will permit drilling with water. However, because of the dispersion of drilled solids, it is usually not possible to maintain clear water even with the aid of mechanical separators, or by utilizing large settling pits. For this reason the commonly used term "water drilling" is a misnomer since the water usually contains some solids. Because of the improved drilling rates and other benefits which might be effected, a laboratory and field study was made of the use of flocculating agents for maintaining clear water for drilling. PROCEDURE Laboratory—Solids Settling Two techniques were employed to determine the effectiveness of flocculating agents in settling solids from drilling mud. First, the rate of fall of solids was measured in graduated beakers containing mud treated with various flocculating agents. Natural gum, oil-in-water emulsions prepared with ethoxylated nonyl phenols and diesel oil. and a powdered synthetic acrylamide-carboxylic acid copolymer* (hereinafter referred to as "polymer") were tested as flocculating agents in both fresh water and salt saturated muds containing 4 per cent solids. The polymer was not added in the powder form directly to the muds but was first dissolved in water in a concentration of 3.5 lb/bbl. The polymer solution was then added directly to the mud. The concentration of the flocculating agents, including the polymer, was varied. High rates of fall of the solids indicated high flocculating efficiency. A second technique consisted of measuring the time in minutes required to obtain clear water at the pump suction after addition of the various flocculating agents at the flowline. These tests were made using a model mud circulating system having a total volume of approximately 10 gal. The natural fresh water and brine muds used contained 4 per cent solids. A sketch of the circulating system appears in Fig. 1. The circulating rate

By F. F. Skiba, J. E. Warren

An idealized miscible displacement process in a three-dimensional, heterogeneous porous medium has been studied via experimental computation based on a "Monte Carlo"' model. The macroscopic dispersion which results solely from variations in the permeability of the porous system is related to the scale of the heterogeneity as well as the distribution function for the permeabilities. The interpretation of laboratory data is often subjective and may bt? misleading since ambiguities are not always recognized. Furthermore, an experiment performed on a conventional oil-field core does not yield a valid measure of Macroscopic dispersion for reservoir engineering purposes since the scale of heterogeneity that is significant in the laboratory is not significant in the reservoir. The analysis of field-performance data is discussed and it is shown that a unique characterization of the reservoir can not be obtained from this information alone. A possible method based on a minimum of data, for estimating the spatial distribution of permeabilities in a reservoir is postulated INTRODUCTION Many studies concerned with diffusion and dispersion in porous media have been undertaken because of the potential importance of miscible displacement processes for increasing the recovery of oil. The results from a large number of these studies have been published; in a recent paper Perkins and Johnstonl presented a comprehensive survey of the pertinent literature and excellent appraisal of the "state of the art". In this investigation, dispersion resulting from macroscopic variations in the properties of the porous medium is of primary interest. Dispersion, in this sense, is caused by fluctuations in the velocities of the individual fluid elements as they move through the porous system. This phenomenon is the result of inhomogeneities in the permeability and/or porosity, but it can be amplified by differences in the characteristics of the fluids involved. The objectives of this study are the following: (1) to determine the qualitative manner in which macroscopic rock properties affect the observed dispersion coefficient; (2) to evaluate the possible influence of macroscopic dispersion on laboratory experiments; and (3) to appraise the significance of macroscopic dispersion with regard to reservoir performance. Investigations of this type are essential since many miscible displacement processes, while seemingly understood from a mechanistic point of view, do not behave in a predictable manner under field conditions. Although a plausible explanation is that the process is dominated by the environment, little serious effort has been devoted either to describing the reservoir itself or to determining the probable effect of the reservoir properties on the process. This study is intended to be exploratory in nature; however, it is hoped that it will provide additional insight into the physical problem and that it will suggest areas suitable for further investigation. EXPERIMENTAL COMPUTATION MODEL To limit the observed dispersion to that which results solely from variations in the properties of the porous medium, an idealized process in which one incompressible fluid displaces another identical fluid is simulated. This simulation is achieved by means of a modified Monte Carlo method — by tracking a number of mathematical particles (tracer) through the system according to the steady-state, single-phase velocity distribution. Then the system is characterized by the distribution of the residence times of the particles. Since any number of particles can be injected simultaneously (a true delta-function distribution) and only pressure conditions are imposed on the boundaries, the model is completely determinate. Furthermore molecular diffusion, microscopic dispersion, adsorption and/or reaction are automatically precluded from the model. The heterogeneous medium is represented by a rectangular, three-dimensional array of uniform, homogeneous porous blocks which are randomly arranged; it is assumed that any two opposing faces of the system are held at fixed potentials

Jan 1, 1965

THE Woman's Auxiliary now has its own publica-tion, a monthly news letter that is mailed direct to all its members. It is a lively and readable sheet. Any member of the Auxiliary failing to receive it should communicate with Mrs. J. E. Johnson, Jr., 49 Church Lane, Scarsdale, N. Y., and get on the mailing list. The detailed reports of the meetings will hereafter ap-pear in the News Letter rather than in MINING AND METALLURGY. The regular Board meeting held on March 3, with Mrs. Barron presiding, was well attended and the re-ports of the committees, and from local sections, showed a keen interest in the work. of the organization. Steps are being taken to endow a "Mary Young Westervelt Scholarship," libraries are being built up in mining towns, and interesting reunions .are being held by the members in the various local sections. The Auxiliary is filling an important place and doing it well. Women of the families of Institute members not already affili-ated with it would find it pleasant to get in touch.

Jan 4, 1927

By T. E. Seidel, A. U. MacRae

The dependence of the electrical and crystalline properties of silicon containing ion implanted boron atoms have been studied as a function of the incident dose, substrate temperature, and annealing tempera-ture. Hall effect studies show that the compensation due to defects is negligible after annealing at and above 500°C for a dose of 1013 ions per sq cm. Steps occur in the isochronal annealing curve for the sheet resistance at 300°C and 800°C and the mechanisms responsible for these steps are discussed. The range of the ions at 150 kev is 3500 * 5001. Electron diffraction studies on samples implanted with 150 ker ions show the existence of damaged (polycrystalline) silicon closer to the surface than the range of the ions. Etching and subsequent annealing of etched and unetched specimens suggest that the poly crystalline silicon is closer to the surface than the bulk of the distribution of ions. THE purpose of this study is to investigate some electrical and crystalline properties of ion implanted boron in silicon. Sheet resistance, Hall measurement, depth distribution by chemical staining, thermal probe and electron diffraction have been used to study the properties of samples implanted in a random direction with ions having energy up to 150 kev. The boron-silicon system has been the subject of considerable study'-7 due to its possible application in the fabrication of electronic devices. The implanted ions produce considerable crystalline damage to the substrate crystal before they lose most of their primary energy and come to rest. Fortunately, most, if not all, of this damage can be removed by subsequent annealing or by holding the substrate at an elevated temperature during the implantation process. Sheet resistivity measurements of the p-type layers produced by boron implantation into room-temperature silicon, reveal that the annealing occurs in a two-step process: the first occurring at -300°C and the second at -800°C.' Electron-microscope and low-angle electron diffraction studies of this damaged material indicate that the bombarding ions produce an amorphous silicon region about the trajectory of each ion.' At high doses, these damaged regions overlap and an amorphous surface layer results. Furthermore, the 800°C annealing step, appears to be associated with the elimination of this gross damage. The depth distribution of the as-implanted and also the annealed samples have been studied by an angle lap and staining technique as well as by Hall measurements.3-7 We find the experimentally determined ranges of the boron are less than those calculated by the Lindhard-Scharff-Schiott theorys (hereafter referred to as the LSS theory). Our results are in essential agreement with pre- vious results. In addition, we have made Hall measurements on samples bombarded with a low dose (2 x 1013 ions per sq cm) to minimize impurity banding effects, and find that deep ionized levels influence the carrier concentration and room-temperature hole mobility until the sample is annealed at temperatures above 500°C. Above this temperature, the data corresponds to that obtained when silicon is doped with boron by conventional techniques. We have also made a depth distribution determination of the crystalline damage due to the implanted boron by using glancing angle electron diffraction techniques. It was found that the gross damage depth distribution of heavily implanted samples is closer to the surface than the depth distribution of the electrically active boron. By removing this gross damage, without affecting the distribution of the implanted boron, we have determined that the 800°C annealing stage is not exclusively due to the improvement in the crystalline perfection of the amorphous material but also appears to be associated with the annealing of some point defects. These results and their relation to previous results are discussed in the following sections. IMPLANTATIONS Boron ions were implanted into silicon samples using an accelerator capable of operating from 10 to 150

Jan 1, 1970

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By W. B. Nelligam, J. Tittman

Refinements in radiation logging techniques during recent years have involved increasing usage of scintillation ditectors. These detectors produce voltage pulses whose heights are related to the energies of the gamma rays which initiate them. Analysis of the gamma-ray spectrum, as indicated by the pulse heights, yields information about the chemical elements composing the formations surveyed. Refined scintillation counter techniques can furnish chemical information concerning earth formations in situ, from a study of the gamma-ray spectra emitted by the formation either naturally or as a result of neutron bombardment."' Accompanying the rising interest in gamma-ray scintillation spectroscopy, there has been increased activity in the development of accelerator-type neutron sources (in contrast to encapsulated chemical-mixture sources). Such neutron generators are attractive for several reasons: (1) they greatly reduce radiation hazards to personnel; (2) there is a great reduction in contamination danger if they are lost in the hole; (3) they can produce larger neutron intensities than can conveniently available encapsulated sources; and (4) they are capable of being pulsed, thus permitting new techniques in logging. In the past, both accelerator and encapsulated neutron sources have been used by others in conjunction with scintillation-detector pulse-height analysis. The results have not been too encouraging, due to the interference among different gamma-ray spectral lines and to the fact that the gamma-ray peaks were not too clearly distinguishable above the large and ill-defined background "noise".'.' This paper is a status report on laboratory studies of a technique using a borehole accelerator as a neutron source, which gives an improved scintillation spectrum, thus permitting more accurate chemical analyses of the formations penetrated. The Schlumberger-accelerator neutron source' is presented; the origins of inelastic and of thermal-neutron, capture gamma rays are dis- cussed, and results are given for some laboratory measurements performed in borehole geometry. THE NEUTRON GENERATOR-TUBE ACCELERATOR In the attempt to develop an accelerator neutron source which would have the desirable properties outlined in the "Introduction" and which would operate satisfactorily under borehole logging conditions, a small-diameter, cylindrical, neutron generator tube has been developed. This tube utilizes the principle of accelerating deuterons (heavy hydrogen nuclei) by a high voltage so that they bombard a tritium (heavy, heavy hydrogen) target. The resulting reactions produce large numbers of neutrons of 14 Mev energy. Furthermore, the tube is permanently sealed and, thus, the use of pumping techniques in the sonde is avoided. The tube consists of a pressure control, an ion source in which the deuterons are stripped of their electrons, an accelerating gap down which the deuterons are sped by the high voltage, a secondary electron suppressor and a tritium-loaded target. Obviously, the tube requires auxiliary circuitry in a sonde for the control of the ion source, tube pressure, high voltage and other operating variables. Neutron yields, both continuous and pulsed, have been produced under simulated field conditions in the range between 1 and 10 times those conventionally used in neutron logging applications. For the experiments discussed in this paper, neutron-pulse repetition rates in the range between 500 and 5,000 pulses/sec are adequate. Furthermore, by making suitable adjustments in operating conditions, one can vary the pulse width from a minimum of several microseconds up to dc. As will be seen later, the present application to well logging does not require that one have available pulses of neutrons which are appreciably shorter than the time it takes for fast neutrons to slow down to thermal energy in water. As a consequence, much of our interest has been directed towards the operating characteristics of the tube with pulses longer than about 10 microseconds in duration. Our experience to date shows that tubes can be operated in this manner with reasonably constant, average neutron outputs.

By H. S. Price, J. E. Warren

Techniques for studying the performance characteristics of heterogeneous reservoirs have been developed. The effect of permeability variation on both the steady-state and the transient flow of a single fluid has been investigated. Limited comparisons with field and laboratory data have been made. The physical models studied consist of random three-dimensional arrays of homogeneous porous blocks. The permeabilities of the individual elements are assigned according to a specific distribution function; uniform anisotropy is introduced by varying the relative dimensions of the blocks. A particular model is perturbed simply by re-arranging its elements at random. The behavior of a physical model is determined by digitally solving its numerical analogue. Based on computational experiments, subject to the restrictions implied by the assumptions that were made, the following general conclusions have been drawn. 1. The most probable behavior of a heterogeneous system approaches that of a homogeneous system with a permeability equal to the geometric mean of the individual permeabilities. 2. The effects of flow geometry and anisotropy on the most probable value of the effective permeability of a heterogeneous medium are finite but not significant. 3. The permeability determined from a pressure build-up curve lor a heterogeneous reservoir gives a reasonable value for the effective permeability of the drainage area. 4. A qualitative measure of the degree of heterogeneity and its spatial configuration are obtained from a comparative study of core analysis and pressure build-up data. It has been indicated that these conclusions are predicated on the assumption that the core analysis and the pressure build-up data represent true reservoir characteristics. Common sources of error in these data have been discussed. One of the most significant problems of reservoir engineering is that of determining the nature and the disposition of the heterogeneities that inevitably occur in petroliferous formations. It is with this problem that this paper is concerned. Extensive research, both theoretical and experimental, has provided reasonable descriptions of the physical processes that are associated with primary and secondaryrecovery mechanisms. Similar progress in the fie1d of numerical analysis and the evolution of truly high-speed digital computers have made it feasible to calculate the behavior of these mechanisms under almost any environmental conditions. It appears to be logical, then, to conclude that the mechanistic performance of any reservoir can be predicted with some confidence. Unfortunate1y, a necessary condition for the validity of this conclusion is that the formation be adequately described in a macroscopic sense. This condition cannot be satisfied ill general; therefore, a frustrating situation currently exists — it is possible to compute the effect of reservoir heterogeneity on behavior, but it is not possible to specify the heterogeneity itself. In many cases, the predicted behavior of a project is so completely dominated by irregularities in the formation's physical properties that the gratuitous assumptior, of a particular form for the variation can reduce the solution of the problem to a mere tautological exercise. While all porous media are microscopically heterogeneous, only macroscopic fluctuations need be considered since the fundamental concepts of reservoir mechanics are based on macroscopic quantities. For the simplest case in which only one physical property varies from point to point, the degree of heterogeneity that is apparent depends on three distinct factors, i.e., the nature of the property's variation (frequency distribution), its disposition (spatial distribution) and the inherent stability of the mechanism being studied. Consequently, any attempt to devise an unambiguous scale of heterogeneity must take these factors into account. Since this investigation is intended to be exploratory, it will be restricted to single-phase flow because of the simplification implied; both steady-state and unsteady-state behavior will be considered. Furthermore, to keep the number of parameters

By R. Wachtell

IN order to study better the phenomena at work in various phases of diffusion of the Fe/Si system when compounded and alloyed by powder metallurgy methods, several attacks have been planned. Electrical, hardness and other measurements have been discussed e1sewhere.l,2 The metallographic phenomena will be considered here. Concerning the metallography of the ferrosilicon alloys, some initial discussion is in order. As this laboratory discovered to its chagrin, there are many traps that lie in wait for even the experienced metal-lographer when he deals with alloys of high silicon content. Chief of these is probably the anomalous behavior of these materials under conventional etching techniques. We have before us, for instance, Corson's paper" of 1928, wherein is mentioned and illustrated (probably for the first time) the peculiar "barley shell" structure sometimes found in these alloys. Again in 19412 he same author reports the structure in 5-14 pct silicon irons. Houghton and Becker mention its presence' but make only tentative effort to explain it. In 1943, Hurst and Riley6 discussed the "barley shell" and declared it to be a false structure, resulting from the peculiar and characteristic film-forming propensities of the alloy. At about the same time, Wrazej7 proved correct the contentions of Hurst and Riley, establishing pretty well the mechanisms of formation and identifying the constituents involved. Hurst and Riley6 a1so describe in their paper (1943) a "cracked film" structure, also a pseudo-morph, the precise nature of which they do not suggest. Careful examination of such a film reveals that it not only cracks but begins to curl away from the surface of the metal and peel, much in the manner of flaking paint. The careful microscopist will be able to focus on the topmost edge of the flake, and, by optical sectioning, follow it down to the metal surface. We have been able, by oblique illumination of such specimens, to observe the raised segment of the film, the shadow that it casts, and the mating segment on the metal into which it fits. As late as 1946, Hurst and Riley8 found occasion to correct misinterpretation of this "cracked film" structure in a work by others on a 12 pct Si/Fe alloy. The "barley shell" and "cracked film" pseudo-morphs are not the only ones encountered in studies of these materials. A third pseudomorph, which may, however, have some structural significance is a fine striation of the surface. Closely spaced and parallel striae extend entirely across single grains. Each grain shows a different direction for the striae, and different closeness of the spacing. The regularity of this structure, and the fact that the striae do not cross from one grain to another, suggest that the structure is a film cracking controlled in some way by the crystal plane orientation of the surface being etched. In general it is wise for the metallographer dealing with these alloys to cultivate a feeling of uncertainty not only of the more complex questions of interpretation, but also of the basic question of whether that which he sees is indeed a true representation of the structure. We have taken, as a first test of the validity of a structure, its reproducibility "in situ" after repolishing and re-etching of the samples. Thus we knew the "barley shell" to be false (fig. 1) even before encountering the definitive works on the subject. The Metallography of Silicon Diffusion in Laminate Bars: As a first step in the study of the metallography of the diffusion process, we have pressed a series of laminate bars. These bars, or "sandwiches" as we have termed them, were pressed with covering layers of iron powder, and a "filler" of the Fe/Si master alloy under study. Bar dimensions were approximately %, x 3 x 1/4 in. The total

Jan 1, 1951

By R. W. Parsons, P. R. Chaney

Oil recovery by the imbibition mechanism can be important in fractured carbonate reservoirs with a bottom water drive. Laboratory experiments were performed on water-wet carbonate rocks to model this process. The results are interpreted in light of the applicable scaling laws. The imbibition production behavior of a pillar of rock subjected to a slowly rising water table can be synthesized from data on total immersion experiments on small rock samples. If the rate of water table rise is slow enough and if the physical properties are such that the imbibition zone is of a small vertical extent, then only the final oil recovery on the small samples need be known. Applying laboratory imbibition data to a real reservoir should be done with caution because of possible wettability problems and the unknown behavior of two-phase flow in fractures. INTRODUCTION The spontaneous taking up of a preferential wetting fluid into a porous medium with the simultaneous expulsion of the contained fluid is termed imbibition. One of the many physical imbibition processes is the production of oil by the taking up of water into a water-wet reservoir rock. In certain types of heterogeneous reservoirs, e.g., fractured, it is realized that imbibition alone may be the dominant mechanism. Theoretical studies of imbibition oil recovery have followed several lines of attack. The behavior of certain mathematically tractable and conceptually simple models have been analyzed. Birks1 considered a bundle of vertical capillaries model, while Perotti et al.2 and Kelemen3 used a model of a series of vertically oriented fractures. If the oil production behavior from a block of reservoir rock is known, then the cumulative behavior of a vertical stack of these blocks subjected to a rising water table can be determined by a method outlined by Aronofsky et al. 4 Knowing the capillary pressure and relative permeability functions should uniquely set the imbibition behavior for a given sample. Solutions of the partial differential equations describing the process, however, present some formidable problems. Digital computer techniques have been applied to displacement processes in which the capillary pressure has been included.5"8 More recently, similar techniques have been used by Blair 9 to calculate countercurrent imbibition behavior. The fact that imbibition behavior is described by the usual fluid flow equations (Darcy's Law and continuity equation) implies that the derived scaling laws should apply equally to this situation. Use of model experiments to study imbibition has been employed by Graham and Richardson 10 and Mattax and Kyte,11 who actually tested the scaling law applicability. This study is concerned with scaled model imbibition experiments on water-wet carbonate rocks. Of particular interest is the slowly rising water table experiment. This simulates a highly fractured reservoir with the water table rising uniformly through the vertical fracture system. Experiments were performed on different sized samples and with different rates of rise of the water-oil interface surrounding the rock pillar. Total immersion imbibition experiments on smaller rock samples were performed to test the basic linear scaling laws and to see if these simpler tests could be used to synthesize the rising water table experimental results. THEORY IMBIBITION SCALING LAWS Scaling laws for the flow of two incompressible immiscible fluids in a porous medium have been presented by Rapoport.12 These are derived directly from Darcy's Law for the individual phases and the continuity equation. Results from a laboratory model experiment will duplicate those from some reservoir prototype if a constant proportion is maintained among the three forces acting on the fluids — the capillary pressure gradient, the gravitational gradient and the flowing pressure gradient. Operation of the model (i.e., fluid injection or withdrawal) must be conducted in accordance with specific equations. If the

Jan 1, 1967

By A. S. Michaels, M. C. Porter, A. Stancell

Previous laboratory studies have demonstrated that the injection of small quantities of reverse wetting agents during water displacement can increase oil recovery from unconsolidated porous media.9 In the present investigation, an attempt has been made to determine more fully the effects of reverse wetting treatments and to clarify the mechanism by which increased oil recovery is effected. Water-oil displacements were performed in beds of 140-200 mesh silica sand. Hexylamine slugs (injected after 0.25 pore volume of water throughput), when adequate in size and concentration, were effective in promoting additional oil recovery. Their effectiveness increased with the quantity of amine injected. However, slugs of sufficient size and concentration to stimulate oil production at water flow rates of 34 ft/day did not do so at 4 ft/day. Visual studies in a glass grid micromodel have shown that the stimulation of oil production, via ayueous hexylamine, is a result of transient changes in the oil wettability of the pore walls. if the amine slug is of sufficient size and concentration to induce significant changes in the adhesion-tension, large continuous oil masses will be formed. if the superficial water velocity is high enough to result in rapid desorption of the amine, a favorable "wettobility gradient" may be established across the masses; under such conditions, high oil mobility is observed, and increased oil recovery results. INTRODUCTION It is generally agreed that the efficiency of oil displacement by water in porous media is limited in part by capillary forces which cause the retention of isolated masses of oil — resulting in the so-called "irreducible minimum oil saturation". Recent estimates10 indicate that there are about 220 billion bbl of petroleum in United States reservoirs which are not economically recoverable with present techniques (such as water flooding). This amounts to almost five times the known recoverable reserves. It has been recognized for some time that a suitable alteration in the water-oil interfacial tension and/or the contact angle, as measured between the water-oil interface and the solid surface, should result in better displacement efficiency. Surface active agents can be used as interfacial tension depressants to accomplish this objective,1,2,7 but unfortunately, the additional oil recovery is seldom commensurate with the treatment cost. In contrast to interfacial tension depressants, the effect of contact angle alterations on water-oil displacements has received relatively little attention in the literature. It is known that the wettability affects the displacement process.3,4,11 Displacements in water-wet systems generally result in lower residual oil saturations than those in oil-wet systems. The effect of "'transient" wettability alterations concurrent with the displacement process have been investigated by Wagner, Leach and coworkers,4,11 wherein it has been demonstrated that the establishment of water-wet conditions during water flooding of oil-wet, oil-saturated porous media is accompanied by significant increase in oil displacement efficiency. Michaels and Timminss studied the effects of transient contact angle alterations resulting from chromatographic transport of reverse wetting agents through unconsolidated sand. It was demonstrated that chromatographic transport of short-chain (C4 through C8) primary aliphatic amines can improve oil recovery and that the recovery increases with the quantity of amine injected (i.e., with either the amine concentration or the volume of the slug injected). Circumstantial evidence indicated that the increased displacement efficiency resulted primarily from transient changes in wettability of the porous medium. In the Present investigation. additional information has been obtained on the effects of reverse wetting treatments and the mechanism by which increased oil recovery is accomplished. Hexylamine slugs of varying size and concentration were chromato-

Jan 1, 1965