Search Documents

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By R. F. Domagala, D. J. McPherson

Iodide zirconium was combined with calculated amounts of ZrO2 or master alloys and arc-melted. Annealing treatments were carried out at 21 temperature levels. Metallographic examination of the heat treated specimens permitted construction of the binary phase diagram from zirconium to ZrO2. Oxygen additions to zirconium raise the transformation temperature as well as the melting point. Features of the diagram include the peritectic formation of beta, the formation of alpha directly from the melt, an intermediate phase ZrOB with a' range of homogeneity, and a eutectic between alpha and Zr02. PHASE relationships in the Zr-0 system were determined in the range 0 to 66. 7 atomic pct 0 (ZrO,). Arc-melted alloys were annealed at temperature levels between 600" and 2000°C. Determination of the phase boundaries was accomplished by metallographic evaluation of specimens quenched from the various temperatures. Incipient melting techniques were used to determine solidus curves, and X-ray diffraction work was employed to study the lattice parameters of the a solid solution and the phase ZrO,. Experimental Procedures Westinghouse "Grade 1" iodide zirconium crystal bar (approximately 99.8 pct pure) was employed for the investigation. This material is substantially free of the impurity inclusions characteristic of most zirconium metal. Previous work has shown it to be the most satisfactory grade for phase diagram studies. The zirconium crystal bar, as-received, was coated with corrosion product from autoclave tests by which its grade designation is determined. The bars were sand blasted lightly, pickled for 1 min in a 20 pct HN0,-5 pct HF aqueous solution, rinsed in water and acetone, and dried. The bars were rolled to about 1/32 in. strip, cut into 10 in. lengths, and pickled again. The material was sheared to approximately Y4 in. squares, cleaned with acetone, and stored for use. Two pounds of specially prepared "chemically pure," and 100 grams of spectrographic grade ZrO, were obtained for the preparation of alloys. The "Specpure" ZrO, was obtained from Johnson, Mat-they Co.. Ltd.. and Titanium Alloy Manufacturing Div. of National Lead Co. supplied the other grade of oxide. The "chemically pure" material (henceforth referred to as C.P. ZrO,) was used for the preparation of all alloys employed in the determination of phase relationships. A limited number of alloys were prepared with the "Specpure" oxide to recheck certain phase boundaries. Analyses of the two grades of oxide are given in Table I. The C.P. ZrOI was received in powder form, not suitable for arc melting. Therefore, this material was compressed to wafers 1 in. in diameter and about Vn in. thick. These were segmented and stored for use. The "Specpure" material, in the form of small granules, was suitable for use in the as-received condition. A nonconsumable electrode arc-melting furnace was used to prepare alloys. A water-cooled copper crucible and tungsten-tipped electrode were employed. The material was placed in the crucible and rapidly melted under a protective atmosphere of helium gas by striking the arc on a tungsten stud and directing it upon the charge material. Details of operation, including drawings of this type furnace, have been published.'

Jan 1, 1955

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By Michael P. Robinson

It has been suggested that streaming potentials are not nomlally logged because the streaming potentials known to be generated across mud filter cakes are substantially cancelled by streaming potentials generated in contiguous shale beds. Laboratory measuments of streaming poterztials acrons hand specirnens of shale have been adduced as evidence to support the suggestion. The experiments show that streanzing potential depends on the resistivity of the equilibrating solution and on pressure drop. Meizsurenlents were made perpendicular to the bedding planes of the shales. The mechanism of streaming potential generation in shale which the suggested cancellation theory iniplies has been critically examined. Conclusions reached are tested by comparing self potentials logged before and immediately after changing the resistil~ity (but not the weight) of mud in a borehole. The resultant SP change is found to equal the theoretical electrochetnical potential change, i.e., cancellation of .strearning potential is independent of mud resistivity. Since streaming po/entials generateal across mud filter cakes depend critically upon mud reistivity, shale .strearlling potential in situ must he equally sensitive. Using plausible values for the pemleahility, porosiy arzd compressibility of shale, the van Ererdingerl-Hurst approach shows that the pressure drop between horehole pressure and formation pressure in a radially honzogeneous shale is unlikely to coincide with the zone of filtrate invasion. If, perchance. this coincidence occurs with one mud a chatzge 10 a second must result in non-coincidence, i.e., a .YP logged in the second mud should show a streaming potenlial. The dilernrna may he overcome by postrrlating: (I) perrneabilities of shale7 parallel to their beddine planes greatly exceed their cross-bedding permeabilities; (2) a resistivity-sensitive mud cake of low permeability forms either on the face of, or probnhly in the lanlinations of, many shales. These postulates are examined with reference to the pres.sure and SP histories of the McClosky and Aux Vases zones in the Illinois Basin and some srcpport for them found. It is concluded that a definitive solultion to the streaming potential problem does not exist. INTRODUCTION Recently Gondouin and Scala' published a contribution to the streaming potential problem which is notable particularly for the elegance of its experimental techniques. The conclusions drawn by Gondouin and Scala from their studies amplified those previously suggested by Schehck'. In essence they suggested that the streaming potentials known to develop across mud filter cakes" are not normally observed on SP logs because they are largely cancelled by streaming potentials generated in adjacent shales. This suggestion appears so reasonable that many workers in the logging field will be prone to accept it. We incline strongly to the view that the suggestion is basically correct. Nevertheless, it is important to realize that the cancellation hypothesis in the form put forward by Gondouin and Scala does not appear to meet all the known experimental facts. It is the principal intention of this note to call attention to what appears to be a flaw in the relevancy of the experimental evidence adduced by Gondouin and Scala and to suggest an alternative hypothesis to remove it. The consequences of the new hypothesis are compared with field data which, we believe, are unique. A second but not less important purpose of this note is to emphasize that even the new hypothesis does not fully explain all field observations. Hence we conclude that further study of the streaming potential problem is required before it can be said that a definitive solution to it has been attained. The bare experimental bones of laboratory data on streaming potentials appear to be as follows. 1. A streaming potential is developed across a mud filter cake. This streaming potential is virtually independent of cake thickness but depends on the mud resistivity. The relationship is of the form, 2. A streaming potential is developed across a shale. This potential is independent of shale thickness but depends on the resistivity of the permeating solution if the shale has been first equilibrated in this solution. The relationship is of the form, The relationship in (1) was determined first by Wyllie3, confirmed by Gondouin and Scala and qualita-

By N. Mungan, J. L. Thompson

Laboratory displacement tests were performed to study oil recovery elficiency by gravity drainage in fractured systems under miscible conditions. The porous media used were cylindrical Berea and Blue Jacket sandstone cores containing a number of well-defined, artificially formed, vertical and subvertical fractures. Butane and Soltrol 130 were the two miscible fluids used. The purpose of the study was to examine the influences of the displacement rate, fracture density, fracture orientation, fracture permeability, matrix permeability, crossflow, core length and connate water on the oil recovery. It was found that displacement rate, matrix permeability and the subvertical fractures affected oil recovery most. The critical flow rate, based on the matrix permeability, was found to be a significant factor in the process. For displacement rates below the critical flow rate, the oil recovery efficiency appeared to be unaffected by the density of the subvertical fractures. At the high displacement rates, the fracture density becomes important, with the recovery being the most efficient in the core having the greatest number of the subvertical fractures. The magnitude of the fracture permeability, the fracture orientation, the core length and the connate water have little effect on the oil recovery efficiency. INTRODUCTION Gravity drainage under miscible conditions from relatively thick reservoirs can be a very efficient recovery process, especially at low flow rates where the gravity forces are dominant and, consequently, the adverse viscous fingers associated with the unfavorable viscosity ratio are minimized or eliminated. Such a process might involve injection of an LPG or some solvent bank at the crest and then driving the bank with dry gas. For some favorable combination of reservoir temperature, pressure and oil composition, there may be no need to inject any solvent bank, since enrichment of the dry gas by the light ends of the crude creates an in-situ solvent bank. Slobod and Howlett2 have studied the effects of gravity segregation in vertical unconsolidated porous media under miscible conditions in the laboratory. The main variables in their study were the viscosity ratio, the density differences and the rate of flow. In this study, the objective was to find the influence of fractures on gravity drainage under miscible conditions. A greater number of related variables were also studied. LABORATORY STUDY CORE DESCRIPTION Displacement tests were performed mostly in cylindrical Berea sandstone cores. In a few cases in which a much lower matrix permeability was desired, Blue Jacket sandstone cores were used. The vertical fractures were formed by cutting the cores with a thin circular blade or by parting them along a bedding plane. The vertical fracture plane always contained the axis of the cylindrical core. The subvertical fractures, cut with a saw, were inclined 45' from the horizontal. The line of intersection between the planes of the vertical and subvertical fractures was parallel to the circular faces of the cores in all cases except one, as depicted by VF3HF-2 in Fig. 1. Fracture geometry of the other cores is also given in Fig. 1. Tables 1 and 2 give the physical properties of the solid and fractured cores, respectively. The suffixes and prefixes shown for each core have been used throughout the paper so that the reader may discern the fracture geometry from the core numbers.

Jan 1, 1970

By G. T. Pruitt, T. L. Irby, P. B. Crawford, H. Ferrell

In water flooding peripheral, center-to-edge, line-drive or water-encroachment patterns the question has arisen, "when should a producing well be converted to a water-injection well?". It is realized that, if the producing well is shut in when water breakthrough occurs, some oil could be left in an unswept area. Yet, operation to high water-oil ratios may require processing much additional water. If a well is shut in, "where should the allowable be transferred?". A review of the literature indicates that very few data are available to the engineer to assist him in answering these questions. It is the purpose of this paper to present information which may be of value in selecting the flooding technique and the injection and producing wells and in scheduling injection and production rates. Certain phases of these problems can be studied in the laboratory. The Potentiometric model is well suited for studying many waterflooding problems. This model assumes that steady-state flow exists, that the mobility ratio is equal to one and that gravitational effects are neglected.'-' The three-phase sand model also has been recognized as an effective means for studying certain reservoir problems. Both models were used in making the study reported here. This paper discusses seven cases utilizing various waterflooding techniques. Data for the first five cases were obtained on the potentiometric model. Data for Cases VI and VII were obtained by constructing and flooding an unconsolidated sand pack (a three-fluid-phase system) and by determining the swept area and cumulative oil recovery for various types of operations. APPARATUS AND PROCEDURE The potentiometric model reservoir was 30 in. long X 10 in. wide by about 1/2 in. deep. Four wells were equally spaced in a line array along the lateral center line. (See Fig. 1.) Well No. 1 initially was an injection well for each of the patterns reported. The cases studied with the potentiometric model and sand model are described herein. The laboratory sand model reservoir was a 30- 10-in. rectangular Lucite vessel packed with sand to a depth of about 0.65 in. Four wells were equally spaced in a line array along the lateral center line. The reservoir material was an unconsolidated sand into which the reservoir fluids were mixed to yield fluid saturations of 55 per cent oil, 25 per cent water and 20 per cent gas. The water was added to the sand and the material mixed, after which the oil was added. Kerosene and water, the fluids used, gave a viscosity ratio of 1.64:1. In packing the reservoir, care was taken to assure constant fluid saturations and constant porosities. A new mixture was prepared for each run. During water flooding, a colored water was injected by a constant rate pump. The injection rate was 1.1 B/D/ft of sand thickness. For this series of tests, each of the producing wells was open to the atmosphere and, as in the potentiometric model study, Well No. 1 initially was an injection well for each of the patterns. KESULTS OF POTENTIOMETKIC MODEL STUDlES Two cases (Cases I and 11) involving equal pressure injections were studied, with the following results. Case Water injection was initiated at Well No. 1; the producing wells (Nos. 2, 3 and 4) were operated at constant and equal bottom-hole pressures. When water broke through at Well No. 2 and achieved the selected water-oil ratio, Well 2 was converted to a water injection well, resulting in two injection and two producing wells. The injection pressures at Wells 1 and 2 were maintained equal and constant; pressures at the producing wells were held equal, but they were lower than those at the injecting wells. A similar process was repeated when water broke through at Well No. 3. The flood was terminated when

By J. T. Michalak, H. W. Paxton

The recovery of the initial flow stress of poly-crystalline iron is characterized by a) a logarithmic time dependence; b) an increasing activation energy with increasing recovery; c) an increased ?,ate and fraction of recovery with decreased temperature of deformation; and d) an increased rate with increased deformation. Isochro)ml studies of the recovery of single crystals revealed more sluggish recovery. The original flou3 stress of single crystals could be regained by recovery anneals at 800"or 900°C. ThE great majority of the theoretical treatments and experimental studies of strain-hardening have been concerned with the close-packed lattices.' Very little effort has been expended on the bcc structure so that the theory of work-hardening of this lattice is essentially nonexistent. The limited amount of experimental data resulting from the direct studies on silicon-iron,2 low-carbon steel,3 and purified irons4"7 does not as yet permit a detailed analysis of the strain-hardening. Studies concerned with the effects of annealing following deformation have been few in number," have been somewhat limited in extent, and have had the disadvantage of using irons containing a substantial impurity content. The present research was undertaken to augment the limited knowledge we have of the characteristics of a bcc metal as regards strain-hardening and recovery from the strain-hardened state. The recovery of the flow stress has been employed since it was felt that the stress associated with the hardened state gives at least a fair qualitative description of the state of the imperfections causing the hardening. The study was concerned primarily with the recovery of flow stress following strain-hardening of polycrys-talline, zone-melted iron as a function of the amount of deformation, the temperature of deformation, and the time and temperature of the recovery anneal. A limited investigation of the recovery of the initial flow stress of single crystals of this same material as a function of the recovery temperature was also conducted in conjunction with the Lambot X-ray technique." MATERIAL The zone-melted iron used in this investigation was obtained through the courtesy of the American Iron and Steel Institute and Battelle Memorial Institute. The analysis, as given by the supplier indicated a maximum impurity content of 0.02 wt pct excluding those impurities listed in Table I. (The actual impurity content is very probably less than this since detection limits were considered as the amount of a particular impurity level.) Also listed in Table I are the interstitial impurity contents of the iron at various stages of experimental work indicating that these im-

Jan 1, 1962

By Carl E. Swartz

For many years kettles used in the melting and refining of lead and other nonferrous metals and alloys have been made of cast iron. The logic of this probably lies in the fact that cast iron has been known as a material of construction for vessels of this type much longer than steel, and satisfactory fabrication from it has been more easily accomplished and better understood. As the capacity of lead refineries has increased, it has been only natural to increase the size of kettles. Today kettles of 150 to 250 tons are not uncommon. As the size of the kettles increased and better shapes were designed, the uniformity of material of construction became more important. Records kept over comparatively long periods of time have shown that the length of service of kettles has not been uniform. Many causes contribute to this and the importance of each has not yet been evaluated. Kettle design, kettle setting, method of firing and material of construction have been changed from time to time to obtain better service, but the variation in length of service is still great. A kettle may fail during the first few heats or may last through more than one thousand heats. The purpose of this article is not to present the conclusions of a finished piece of research on kettle manufacture and use, but to present a series of observations made over considerable time in the hope that others will publish data and that all may be correlated into a logical scheme which will result in better service and longer kettle life.

Jan 1, 1931

By Irving Hale

The superiority of electric power for mining purposes was recognized in a general way as soon as the electric motor be-

Jan 1, 1897

By W. A. Backofen, D. Lee

The deformation of textured sheet of titanium and several of its all-a (hcp) alloys was investigated with two kinds of experiments: tension testing along different directions in the plane of the sheet, and combined-stress loading along the principal axes of anisotropy (the rolling and transverse directions). Texture and plastic anisotropy were evaluated with pole figures and by measurement, as a function of orientation, of the tensile transverse-strain ratio, or R parameter. The results as a whole can be rationalized qualitatively with the present understanding of deformation mechanisms in titanium. Quantitative agreement between the orientation dependence of tensile yield strength and the Predictions of continuum ANISOTROPY in the properties of engineering materials is more the rule than the exception because anisotropy in structure is so common, and the one generally follows the other. Two broad categories of structural anisotropy are represented by mechanical fibering and crystallographic texture. The former is important in many aspects of fracture.' The latter has more significance for yielding and plastic deformation, theory, using the R measurements, was relatively poor. However, in tension-tension quadrants of the yield-locus diagrams, the difference between theory and experiment was hardly more than 10 pct. The maximum texture-hardening effect was found in a 4Al-1/4O2 alloy to be nearly 50 pct. In the compression-compression quadrant, {1012} deformation twinning prevented the locus from expanding much beyond the Position predicted for isotropic material. With increasing temperature the trend in R was downwards, approaching unity (isotropy) near the a-transus temperature; an accompanying development was large amounts of neck-free extension. which are of special concern in this paper. The tendency has been to view crystallographic texture more as a liability or curio than as an asset. The prominent exception is found in silicon-ferrite sheet for power-transformer cores, which is processed to develop a new optimum texture for high electrical efficiency. Yet it should be possible to find other examples, involving other kinds of properties. More specifically, there is good analytical reason to expect improved combined-stress yielding resistance in suitably textured material,2,3 and it has now been shown by experiment that texture hardening or softening is a potentially useful phenomenon in certain engineering devices and processes.'-I To date, work has included measurements of the transverse (width to thickness) strain ratios in the

Jan 1, 1967

By J. P. Matoney, E. J. Sandy, D. A. Dahlstorm

INTRODUCTION Before an understandable discussion concerning coal-dewatering processes can be initiated, it is necessary to define some terms of usage. The many different meanings and interpretations of terms preclude development of standard definitions so that the terms which are defined in this chapter should be regarded mainly as working terms intended to facilitate communication. A coal slurry ' might be defined as a watery suspension carrying particles of fine coal and other assorted material ranging upward in size from the finest clay particles. Some particles may be in the colloidal sue range and exhibit characteristics of suspended colloids. The dilution of the system generally is such that it is in a condition of independent particle subsidence. In some coal mining areas a slurry is regarded as a suspension of fine coal and other materials having a dilution that cannot be economically dewatered by filtration. A coal slime may be regarded as a slurry containing particles of such size range that 50 per cent or more (by weight) will pass a 200-mesh sieve. Slime suspensions generally are difficult to handle without special treatment in conventional sedimentation equipment. A coal sludge is defined as a slurry that has been partly dewatered by sedimentation usually to a thick or viscous consistency where further dewatering can be most readily facilitated by filtration. Culm and silt are terms characteristic of coal-preparation practice which should be defined. In anthracite terminology, culm is the generally coarse residue or waste accumulation of coal, bone and rock from preparation plants most of which are no longer functional. Culm has a wide range of sizes, some rather coarse, and except to differentiate it from silt, the term is of limited interest in this discussion. In anthracite terminology, silt is the accumulation of waste fine coal, bone and slate settled out from preparation plant water in excavated basins. Silt is generally made up of particles ranging in size from 3/32-in. round-opening screen to the very finest sues. Silt is also referred to at times as sludge, fines, slush, and mud and is the partly dewatered derivative of what has been previously defined as slurry. In the terminology of the bituminous coal producing areas, the word culm commonly corresponds to slurry or slime, depending upon the size distribution of the suspended solids. Editor's Note: Because of the numerous and varied contributions of the associate authors of this chapter, the text does not differentiate between sections authored by each individual as is the case in most chapters of this book.

Jan 1, 1968

By John T. Norton, Robert J. Anderson

The purpose of this paper is to report evidence, regarding metal structures, that is contradictory to the amorphous-metal hypothesis of Beilby, and particularly evidence that is opposed to the proposition that a vitreous amorphous layer is produced by polishing a piece of metal. According to the theory of Beilby, which has been developed by Rosenhain and others, polishing a piece of metal greatly affects the condition of the metal at or near the surface polished, making that surface entirely different in character from the metal beneath. A vitreous amorphous surface film is said to be formed. The polishing process, on the basis of Beilby's explanation, is well known and does not require discussion here, but it is advisable to seek the precise meaning of the term "amorphous." According to Beilby, "amorphous means non-crystalline in the most intimate sense of the word; the molecules arc not marshalled in crystalline order and orientation." According to Rosenhain, the amorphous condition is that assumed by a metal when the crystalline arrangement is completely broken up so that the atoms exist in a state of disorderly arrangement similar to that which exists in the liquid state. Thus, from the point of view of the phase rule, the "amorphous phase" is regarded as being identical with the liquid phase. Amorphous metal is often likened to an undercooled liquid. A crystalline substance is one in which the atoms are arranged in space according to some definite repeating geometrical pattern. The existence of amorphous metal has not been definitely proved. The evidence offered here that amorphous metal probably does not exist is of a direct nature and is afforded by x-ray diffraction effects given by a number of metals and alloys when polished, cold-worked, and annealed. Considerable of the experimental work was carried out with substantially pure aluminum and the aluminum alloy duralumin, but many of the observations made on these materials were checked by similar observations on copper, iron, tin, zinc, lead, mild steels, and

Jan 1, 1925

By D. T. Peters, C. R. Cupp

Aging of commercial 18 pct Ni maraging steel and of related binary, ternary, and quaternary alloys has been studied by electrical resistivity, X-ray diffrac-tion, and hardness measurements. The aging process is in three stages: 1) recovery of martensite; 2) two precipitation reactions; 3) formation of nickel- rich austenite. The main objective of this work was to discover the nature of the Co-Mu interaction. Addition of cobalt to the Fe-18 pct Ni-5pct Mo alloy results in acceleration of the pvecipitation reaction by increasing the super-sat~lration 0f molybdenum. The cobalt-containing al- loys show no incubation period for the precipitation reaction, but the alloy without cobalt shows an initial slow reaction. These results are interpreted in terms of precipitation on dislocations. At low aging temperatures, the presence of cobalt leads to the appearance of a second type of precipitate which is a particularly effective hardener. Cobalt was found to retard the formation of austenite. Overaging is a combination of classical overaging and the formation of austenite in the structure. he- activation energy for the precipitation reaction determined from resistivity ,measurements is 42.8 kcal per mole and the growth of the precipitate is found to obey at law at high tenzperatures and a tU3 law at low temperatures. HE mechanism of strengthening of the 18 pct Ni maraging steels has been studied extensively since 1962.'-? The strengthening is due to precipitation hardening of the martensitic matrix. The main precipitate is an inter metallic compound involving molybdenum, but there is disagreement about its identity.

Jan 1, 1967

By C. L. Bulow

In this discussion of stress-corrosion cracking of widely used copper-base alloys, no attempt has been made to prepare an extensive review or bibliography since this has been quite ably covered by other investigators. An excellent bibliography was published by H. P. croft(1)2 in connection with his studies of stress-corrosion cracking in mercurous nitrate solutions and his analyses of the magnitude and distribution of stresses in brass parts. We are reporting on a few of the results which have been obtained to date in connection with an ammonia vapor test and an alternate immersion test. These are accelerated corrosion tests for determining the susceptibility of materials to failure from season or stress-corrosion cracking due to the presence of moist ammonia and oxygen so prevalent in nature. Our experiences with the ammonia test seem to indicate that it is a more searching test, than the mercurous nitrate test.

Jan 1, 1945

By T. E. Leontis

Data are presented in this paper to show that addition of zirconium to sand-cast Mg-Th alloys effects a marked decrease in the grain size of these alloys which is accompanied by a significant increase in the mechanical properties over the entire range of thorium content investigated. The beneficial effect of zirconium on the strength properties is maintained at elevated temperatures up to 700°F. In addition the alloys exhibit exceptionally high creep resistance at temperatures up to 600°F. Zirconium does not greatly improve the strength or creep characteristics of extruded Mg-'Th alloys. Cerium is not a desirable addition to Mg-Th-Zr alloys. IN a previous paper,' it was shown that addition of thorium to magnesium imparts exceptionally high resistance to creep at elevated temperatures both in the sand-cast and in the extruded states. The large grain size of the sand-cast alloys, however, renders them of questionable value as engineering materials. Sauerwald² showed that the grain size can be reduced effectively by the addition of zirconium. Results of a detailed study of the properties of Mg-Th-Zr alloys over a wide range of thorium content are presented in this paper. In accordance with the preliminary findings of Sauerwald,² it is shown here that Mg-Th-Zr alloys possess the unusual combination of fine grain size and a high level of strength and creep resistance at elevated temperatures. Some information on the effect of cerium on the properties of Mg-Th-Zr alloys is presented also. Effect of Th on Solubility of Zr in Molten Mg Sand-cast test bars and 3 in. diam extrusion billets of Mg-Th-Zr and Mg-Th-Ce-Zr alloys* were pre- pared according to the procedures described in the previous paper.' The zirconium was added in the form of sponge metal as supplied by the Bureau of Mines.' The maximum amount of zirconium that usually can be introduced into magnesium at the customary alloying temperature of 1350° F is about 0.7 pct. The analyses in Table I clearly show that thorium in- creases the solubility of zirconium in liquid magnesium. The tabulated zirconium contents are designated as "acid soluble;" this serves to distinguish these values from the small amount of zirconium that is not dissolved when the alloy is subjected to the action of a dilute acid solution. The "acid-insoluble" zirconium is believed to be ineffective as an alloying ingredient. The acid-insoluble zirconium contents ranged between 0.03 and 0.16 pct in the present alloys, the tendency was toward higher values in alloys with higher thorium content. In contrast to the effect of thorium, cerium decreases the solubility of zirconium in magnesium. This is shown by the three Mg-Th-Ce-Zr alloyst in- cluded in Table I and stands in agreement with the observations of Nelson and Stricter".' on Mg-misch-metal-Zr alloys. Effect of Zr on Grain Size of Mg-Th Alloys The potent effect of zirconium as a grain refiner has been demonstrated for many types of magnesium alloys.3-8 Of greatest importance to date are Mg-Zn-Zr alloys5,6,8 and Mg-mischmetal-Zr alloys.3,4,7 Zirconium has also been found to exert a marked grain-refining action on Mg-Th alloys. All the Mg-Th-Zr alloys cast in the form of ½ in. diam test bars have an exceedingly fine grain size (Table I). Reference to the previous paper' will show that binary Mg-Th alloys containing up to 6 pct Th are almost completely columnar even in the M in. diam reduced section of a test bar, and that 10 to 50 pct Th alloys have an equiaxed grain structure of 0.03 to 0.08 in. The data in Table I also show that addition of cerium to Mg-Th-Zr alloys has no effect on the grain size. In addition to increasing the strength properties, the fine grain size of Mg-Th-Zr alloys

Jan 1, 1953

By A. E. Worthington, R. F. Meldau

An analog computer study was made of the reduction of the magnitude of the SP by formation geometry and resistivity effects. Practical field .situations were studied. The analog data were converted to the form of departure charts. These charts may he used to determine the reduction in the SP due to the above eflects pith good precision Quantitative interpretation of the SP log in terms of formation water resistivity is presently based on the premise that the observed SP is equal to the electrochemical potential developed across a shale between the mud filtrate and the formation water. This procedure inherently assumes that, first, no electro kinetic (streaming) potential is present, and, second, that the SP is not reduced by resistivity or geometry effects in the formation. When, however, such a reduction of the observed SP does occur, one may write, SP = b,E, b, E, (1) where SP is the observed self-potential opposite the center of the bed under consideration. In this case, E, is the electrochemical potential, E, is the electrokinetic potential, b, is the electrochemical reduction factor, and h, is the electrokinetic reduction factor. The first assumption above is equivalent to assuming E, = O in Eq. 1 and the second assumption is equivalent to h, = b, = 1. Examination of electric logs shows that for many formations of interest the observed SP is reduced below its maximum value. The object of this work is to aid quantitative SP interpretation in at least some of these situations by obtaining values for b, and b, in Eq. 1 as functions of realistic combinations of bed thickness, depth of invasion, and formation resistivity. The only previously published work in this field known to the authors is by H. G. Doll' in 1948. He evaluated b, under conditions permitting exact analytical solution. His work clearly established the importance and qualitative nature of formation geometry and resistivity effects upon the SP, but the necessity of obtaining analytical solutions restricted his study to situations of limited practical interest. In order to extend such studies to practical combinations of formation parameters which do not lend themselves to feasible analytical solutions, the results reported in this paper were obtained by use of a special analog computer. THE ANALOG MODEL The analog model used for this study was the Guyed Analog Computer which originally was designed and built for resistivity log departure curve studies. It consists of a two-dimensional mesh of resistors, with the vertical mesh corresponding linearly to vertical distance along the borehole, and the horizontal mesh corresponding to radial distance from the borehole axis. By proper choice of horizontal resistor values, a three-dimensional, axially symmetric, borehole-formation problem can be modeled with this two-dimensional array. These and other details of the model are contained in a paper by H. Guyod'. The geometry used for the SP studies with the analog is depicted in Fig. 1. This configuration has a horizontal plane of symmetry, and, hence, the model, and the following discussion, need only to be concerned with one quadrant such as shown in Fig. 2. A further simplification arises by considering that both the natural potential sources and those used in the model have fixed internal impedance. Therefore, the current flows produced by these sources can be separated and studied independently ( i.e., the superposition theorem can be applied). SP = SP, i SP, = h,E, + b,E, ,(7)

By J. O. Brittain, E. P. Lautenschlager, F. Felberbauer

This investigation was undertaken to evaluate the effect of a dispersed phase in a iron upon the strain hardening and the stress dependency of dislocation velocity as inferred from the strain-rate sensitivitv. Experiments were made in a-iron specimens with controlled amounts of carbon and/or A12O3. They were tested at various temperatures between 195' and 373 'K under a variety of conditions oj-concentration, strain rate, and quench aging. The resulting yield strength, flow strength, and strain hardening increased with addition of Al2O3. The parameter m, measured from changes in strain rate, was xsed to describe the stress dependency of the dislocation velocity. m was found to increase with additions of A12O3, increasing carbon in solution, and with increasing temperature. When applied to a recent model for yielding based upon dislocation multiplication and velocity characteristics, the values of m alone did not successfully predict yielding for the materials of this investigalion, but had to be adjusted with a consideration of the number of unlocked dislocations nucleated heterogeneously at discontinuities or inclusions. The change of yield strength with Al2O3 Interpar-ticle spacing appeared to obey a theory of Orowan. EVER since the concept of dislocations has been introduced to explain plastic deformation in crystalline materials, the origin of dislocations and the source of the number of dislocations required to account for plastic deformation have been persistent questions. Ruling out the possibility of homogeneous nucleation,' Gilman2 drew attention to the various cases of heterogeneous nucleation as the source of dislocations, whereas Frank and Read3 introduced the Frank-Read source. However, Gil-man found that the Frank-Read sources did not play a dominant role in the deformation of LiF crystals, but that small precipitates are often, if not always, associated with the nucleation of dislocations. In LiF crystals, however, most of the dislocations arose through the multiplication of moving dislocations in a mechanism described by Koehler4 and Orowan,5 Thus, relatively few precipitates could account for a large number of dislocations. A similar observation of dislocation nucleation at inclusions of precipitates was made by Stein and LOW' on Si-Fe crystals. Their crystals contained

Jan 1, 1964

By K. H. Coats

This paper describes an approximate technique for handling the problem of percolation of evolved gas upwards through the oil column in computer simulation of natural depletion. This technique has been incorporated into a general model for simulating three-phase flow in one, two or three dimensions. The mathematical model for perfonning these calculations is described in detail in a separate article. 1 While vertical gas percolation occurs during natural depletion in a reservoir of any configuration, it is especially pronounced in the pinnacle reel or bioherm. The reel may have an areal extent less than one section with a thickness of up to 800 It. The upward percolation of evolved gas often limits the time step in calculations as little as one day. Using this small time step, computer expense for a single 30- to 40-year simulation on even a one-dimensional basis has exceeded 92,000. The method described here for handling the percolation allows time steps of 60 days or more (depending on reservoir size and production rate), resulting in a considerable reduction in computing expense. The method also allows calculations in which secondary gas caps build up in tight zones of the oil column below the main gas cap. The validity of the method is indicated in connection with an example pinnacle reef field. Calculated results using a 2-day time step (where the method in question is not needed and not invoked) are compared with results using the method and a 60-day time step. The comparison shows good agreement. Results from one- and two-dimensional simulations of the reef are presented along with the corresponding computing times on the CDC 6600 computer. A three-dimensional simulation was also performed and the required computing time is given. INTRODUCTION A general analysis for simulating three-dimensional, three-phase flow in reservoirs has been developed. 1 In applying this model to reservoirs undergoing natural depletion, a time-step restriction was encountered due to the flow of evolved gas upwards through the oil column and toward the gas cap. A remedy for this problem has been incorporated in the analysis and is described here. The time-step restriction is encountered only in calculations which include flow in the vertical or near vertical direction. The restriction occurs to some extent in reservoirs of any geometrical configuration, but it is especially pronounced in the pinnacle reef or bioherm where the ratio of thickness to areal extent is unusually large. This paper describes the problem, a method of handling it, incorporation of the method in the three-dimensional, three-phase model, and, finally, a test of the method's validity in an application to an example pinnacle reef reservoir. THE PROBLEM During the early stages of reservoir depletion, pressure falls below bubble point in progressively lower regions of the oil column. Continuing production results in evolution of dissoived gas throughout the oil column. This gas then percolates upwards toward the top of the reservoir. If the gas encounters a sufficiently low permeability zone in its travel upward, then it can accumulate and form a secondary gas cap. If no zones are encountered which are tight enough to hold gas against the gravity forces, then the gas travels on until it reaches the top of the sand or the main gas cap. Numerical simulation of multiphase flow in reservoirs is generally performed with time steps such that the flow of a fluid into or out of a block in one time step is a fraction (considerably less than one) of the total amount of that fluid present in the block. Use of a time step where the time step's

Jan 1, 1969

By Å. Sterten, R. Tunold, J. Brun, K. Dalatun

The subjects of this study are two corrosion phenomena familiar to operators of aluminum plants employing Soderberg anodes of the vertical type, namely the sul-fide scale formation observed at the steel contact spikes and the corrosion wear observed at the lower edge of the gas collecting skirts. The causes of these corrosion phenomena are discussed on the basis of observations from practice and measurements on commercial cells as well as on laboratory experiments. It is concluded that the anode spike corrosion is mainly due to the attack by sulfur compounds contained in the anodically evolved gases and that this attack is greatly promoted by cracks and pores in the anode carbon embedding the lower spike ends. The corrosion of the gas collecting skirt is mainly attributed to anodic dissolution of the iron during occasional contacts with the melt. ANODE SPIKE CORROSION The visible result of this corrosion is the formation of a brittle loosely adherent scale on the lower ends of the mild steel spikes maintaining the electric contact to the carbon anode, Fig. 1. The thickness of the scale usually varies from 0.5 to 3 mm, but in some cases scale of about 10 mm thickness has been observed by the writers. During the spike pulling, usually repeated every 3 weeks, part of the scale often loosens and remains in the anode carbon, resulting in contamination of the aluminum with iron. This corrosion has been attributed to reaction of the iron with gaseous sulfur compounds originating from the anode paste materials.'-4 This paste usually consists of 70 pct precalcined petroleum coke with 1.2 to 1.4 pct S and 30 pct coal tar pitch containing 0.5 to 0.7 pct S. During operation the gasification of the anode materials takes place in two ways: First, gases are evolved by the gradual baking of the anode paste added on top of the anode. Due to the sealing effected by the upper semifluid zone of the paste, these gases, containing sulfur compounds from the pitch, are forced downward through pores and cracks in the rigid anode carbon. At the temperatures of the zone of the lower ends of the spikes, 700" to 900°C, the original components of the baking gas must undergo decomposition, similar to that taking place with the production of coal gas, resulting in a gas mixture consisting mainly of H2 and CH4 with small percentages of H2S as the main sulfur compound.5 Second, the anodic reaction results in a mixture of C02 and CO, the total quantity of which is about 15 times that of the baking gas. These anode gases must carry the entire sulfur content of the carbon gasified by the anodic reaction. Henry and Holliday6 have shown that COS is the main sulfurous component of the anode gases from Soderberg cells. Soderberg anodes usually exhibit a network of vertical cracks, frequently traversing from one spike to the next one and to the outside of the block and widening downward to the working face of the anode. In addition the anode carbon, particularly that below the spike ends, is rather porous, as indicated in Fig. 1. It therefore seems probable that the anode gases also may partly penetrate to the spike ends. As the anode is partly immersed in the melt, the static pressure of this must forward this penetration. I) SPIKE SCALE COMPOSITION A great number of freshly extracted spikes from industrial plants were inspected and samples of the scale collected for examination. Before analyzing the samples a thin iron oxide layer, formed by air oxidation of the outside of the red hot scale immediately after spike pulling, was removed. On the average the Fe/S ratio found by these analyses corresponded to I the composition FeSl. oo7. X-ray powder patterns showed no oxide and otherwise agreed with those found by Haraldsen7 for iron sulfide with 50.00 at. Pct

Jan 1, 1970

By Charles S. Barrett, Frederick N. Rhines, Robert F. Mehl

The solid solutions which aluminum forms with copper, and with magnesium and silicon, are not extensive, and accordingly could not be expected to form Widmanstatten figures profusely nor with great ease. A study of the Widmanstatten figures in these systems seemed desirable, however, because of the interest which these alloys bear towards age-hardening theories resulting from the necessary relationship between the age-hardening phenomena and the atomic processes leading to the formation of Widmanstatten figures.' Alloys of Aluminum with Copper The most accurate determination of the solid solubi1it.y of copper in aluminum seems to be that by Dix and Richardsoq2 whose results are represented in Fig. 1. According to this work, aluminum dissolves 5.6 weight per cent of copper at 548", and less than one per cent at room temperature. Judging from the constitutional diagram, this solid solu-lion on cooling tends to precipitate the compound CuA12, and this conclusion has been confirmed experimentally by Schmid and Wassermann.

Jan 1, 1932