Search Documents

By R. C. Hassinger, D. V. von Rosenberg

Transverse dispersion has received considerably less treatment in the literature than has longitudinal dispersion. Different methods for determining transverse dispersion coefficients have been used in different investigations, and the results obtained have not been consistent enough to permit accurate generalizations as to the effect of various physical parameters on the magnitude of these coefficients. A numerical solution to the differential equation describing transverse dispersion in the absence of longitudinal dispersion was obtained to enable one to calculate the dispersion coefficient from experimental results. The more general dispersion equation including longitudinal dispersion also was solved numerically to give quantitative limits of a dimensionless group within which the assumption of negligible longitudinal dispersion is justified. Possible experimental procedures were examined, and one utilizing a cylindrical packed column was chosen for the determination of transverse dispersion coefficients. Values of these coefficients were determined for a system of two miscible organic fluids of equal density and viscosity, for two sizes of packing material over a wide range of flow rates in the laminar regime. The dispersion coefficient was found to decrease, for a constant value of the product of packing size and interstitial velocity, as the size of the packing material particles increased. INTRODUCTION Longitudinal dispersion has received extensive treatment in the literature, and consequently is better understood than its orthogonal counterpart, transverse dispersion. Many mathematical models of dispersion processes assume that transverse dispersion is rapid enough to damp out any radial** concentration gradients and therefore may be neglected. Laboratory and production results, however, indicate that this is a poor assumption.l,2 Various experimental procedures for determining transverse dispersion coefficients have been used in previous investigations, but the results have generally been expressed by similar correlations. The transverse dispersion coefficients obtained, however, have often varied considerably for given values of the correlation parameters. We feel that further experimental determinations of transverse dispersion coefficients will help alleviate some of the inconsistencies in these empirical correlations. One assumption implicit in all previous investigations is that of negligible longitudinal dispersion in the experimental system. An attempt to justify this assumption often is made using intuitive reasoning, but it is apparent that this reasoning must break down as the condition of zero flow rate is approached. A mathematical examination of the equations describing the system yields physical limits outside of which the assumption of negligible longitudinal dispersion is invalid. BACKGROUND In a porous medium, the "effective molecular diffusivity" De is less than the molecular diffusivity D measured in the absence of a porous medium, due to the tortuous path which a diffusing molecule must travel. Various authors3-5 have reported values of the ratio De/D in the range of 0.6 to 0.7. When there is fluid flow within the porous medium, mass transfer occurs by convective dispersion as well as by molecular diffusion. These are separate phenomena and can be treated as such on a microscopic scale. However, the mathematical complexity is such that only extremely simple geometries could be considered, and the results hardly would be applicable to the complex geometries existent in actual porous media. A more

Jan 1, 1969

Gentlemen of the American Institute of Mining Engineers.—As you well know an application is about to be made to Congress, by the American Society of Civil Engineers, for the appointment of a cornmission composed of persons skilled in the production and use of iron, steel, and other structural materials, whose duty it shall be to test the strength, elasticity, and other quatities of such materials ; not simply the substances used as such, but the members of large structures composed of such materials, of the forms and sizes in which they are actually used. As I understand, the papers that will be read here to-night will have for their object to show the necessity for such tests, and therefore the necessity for the commission asked for. As the representative of the Society of Civil Engineers, I have to thank you, gentlemen of the Mining Engineers, for the great zeal with which yon have taken up this subject, and seconded the efforts of the Civil Engineers; or, perhaps I ought rather to say, initiated the efforts to get a suitable commission appointed to obtain most vitally important information which it is admitted on all hands we do not now possess. In drawing up the memorial to Congress, in connection with some of the ablest and best informed of our Society, and in all the discossions on the subject, I have been struck with the emphasis with which the confession is made, that so little is known of the strength of the larger-sized members constantly used in building structures of the magnitude which the present requirements of the country now demand. I presume the remarks made to-night will bring this out very prominently. This is no disgrace, for with the means here-

Jan 1, 1882

By C. D. Stahl, S. M. Farouq Ali, W. E. Culham

One- and two-dimensional mathematical models have been developed that simulate transient, two-phase flow of hydrocarbon mixtures in porous media in a manner that accounts for interphase mass transfer. Numerical simulations of one-dimensional depletion-drive experiments using a two-component hydrucarbon fluid were used to establish the validity of the mathematical models. In addition, the experimental and numerical data were used to demonstrate that production rate had a relatively insignificant effect on the recovery of individual hydrocarbon components from the experimental system, and that attainment of equilibrium between phases is possible for a wide range of liquid and vapor velocities in reservoirs containing light hydrocarbon fluids. Results of some two-dimensional numerical simulations are also presented. INTRODUCTION This study was undertaken to develop a mathematical model that would simulate transient, two-phase flow of hydrocarbon mixtures in porous media under conditions that result in interphase mass transfer and to test the validity of the assumptions used to set up the model. In addition, the study was designed to determine if production rate influences the recovery of individual hydrocarbon components from reservoirs producing by depletion drive. Two-phase flow in porous media, with interchange of components between the two phases, is important in many petroleum recovery processes. Studiesl-5 conducted within the last 3 years have outlined methods of solving multiphase flow problems incorporating mass transfer. Some of these studies have also indicated the importance of accounting for mass transfer under various producing conditions. An earlier work6 first demonstrated the importance of combining relative permeability data with equilibrium ratios in compositional balance methods. The mathematical model presented in this paper is formulated so that a phase behavior package, as described in previous papers,5,17 is not required as an integral part of the routine employed to solve for the primary dependent variables. The finite difference formulation is designed so that all the primary variables can be solved for simultaneously. This is accomplished by utilizing one basic set of equations. These innovations, which are in contrast to other models2,3,5 but are similar in some respects to the approach used by Taylor,l render the total problem computationally simpler than any of the previously referenced formulations. The mathematical model was developed by combining Darcy's law with a continuity equation for each hydrocarbon component. The principal assumptions invoked in the formulation were (1) that capillary forces and diffusional effects are negligible, and (2) that thermodynamic equilibrium exists in the reservoir at all times. No assumption as to the type of vaporization process was made in formulating this model. Experimental data were required to complete this study. These were generated by conducting several depletion drive experiments. The experimental apparatus consisted of a sandstone core enclosed in a pressurized casing. The apparatus was designed in such a manner that the core could be charged with a liquid hydrocarbon mixture and depleted at different production rates. The experimental tests were designed to determine the effect of production rate on component recovery. In addition, direct comparison of experimental and mathematically predicted

Jan 1, 1970

By C. S. Matthews, M. Prats, R. I. Jewett, J. D. Baker

By mathematical analysis it was found that injectivity history of a uniform five-spot pattern can be calculated by rather simple formulas. These calculated injectivities were found to agree rather well with injec-tivities measured in a potentiometric analog. With this as a basis, a simple method was finally developed which allows prediction not only of rate of injection but also of rate and kind of production at the production well, for a uniform five-spot. Using this method formulas for the eflective injectivity and production behavior of a waterflooded reservoir having a wide range in permeability can be calculated by considering that the reservoir consists of several layers, each having a uniform but different permeability. Rather good agreement is shown by a comparison of the observed production history of a field in the Illinois basin and that calculated by the methods just described. INTRODUCTION The importance of being able to predict the injection rates for a water flood is well known, for upon injection rate depends the life of the flood, the size of pumping and treating facilities, and the rate of oil recovery. One method of determining rates of injection is through a pilot flood. However, it appears that even if a pilot flood is resorted to, a theoretical method for calculating injection rate will be valuable in extending field results to times after inter- ference and to locations of other five-spots. It is also important to be able to predict the production history of a flood. Such prediction must generally be done by theoretical means, since in general a small pilot flood will not furnish much quantitative information in this regard because of the distortion after oil-bank interference. Methods for predicting the behavior of five-spot water floods have been proposed by Yuster and Cal-houn' and by Hurst.? However, these methods do not consider the effect of the mobility of the different fluids in the reservoir. Other investigators3,4,5,6 have determined the effect of the water-to-oil mobility ratio on the production history by means of different experimental techniques. Most of the published work applies to homogeneous sand bodies, does not provide information for determining the injection requirements of a flood, and seldom considers the rate of build-up of the oil bank as it fills up the partially depleted reservoir.* The present work was undertaken because of the lack of predictive method which takes into account mobility ratio and the rate of buildup of oil bank resulting from the void space in a partially depleted reservoir. As water is injected into the reservoir, which contains oil, water, and gas in macroscopically homogeneous saturation, it forms several banks ahead of the injection well. In each bank or region there is assumed to be only one mobile phase—water in the water bank, oil in the oil bank and gas in the gas region. The saturation in each region and a plan view of the banks are shown in Figs. 1 and 4, respectively. FACTORS AFFECTING INJECTION AND PRODUCTION RATES Results in this report are determined as functions of the three parameters, F, M, and M, , and of a fourth parameter, r,/L** which relates the size of the well to the flood spacing. The parameter F is the displacement factor defined by Hurst in Ref. 2. It determines the rate of build-up of oil bank as it fills up the partially depleted reservoir. When gravity effects are neglected, the four parameters just given are sufficient to describe the waterflood behavior for horizontal and homogeneous reservoirs having incompressible liquids. Definition of InjectIvity The dinlensionless injectivity (or injection rate) used in this report is defined by

By Michael Tenenbaum

THE importance of semikilled steel as a high tonnage grade has long been recognized. The increasing severity of the applications for which semikilled steel is used makes it desirable to obtain further information regarding the features of open hearth practice and ingot structure that can affect steel performance in subsequent rolling and fabricating operations. Accordingly, an investigation is being carried out studying the structure of various types of semikilled steel ingots. This paper reviews some of the observations made and information gathered in this investigation. By definition, semikilled steel may be considered to include all nonrimming steels in which the natural shrinkage occurring during ingot solidification is offset to an appreciable degree by gas formation. Accordingly, this type of steel includes the entire range intermediate between rimmed and fully killed steel. The resultant ingot structures range from those of steels that have been capped to suppress a weak rimming action to those steels that have been killed in the furnace or ladle to the extent that no mold deoxidizer is required. A survey reported to the 1946 AIME Open Hearth Conference' reviewed the practices being used at 28 major plants for the deoxidation of semikilled steel. It was found that most plants divide the deoxidizing additions between the ladle and the mold many of the plants adding only minor amounts of deoxidizer to the bath or none at all. Silicon, aluminum and titanium were the only elements used for ladle and mold deoxidation. Mold deoxidizers were added as capping additions or as uniformly fed additions. Several practices were reported in which the steel was deoxidized in the ladle to the extent that little or no mold deoxidation was needed. No reference was made to grades capped by using special mold designs or special closures over the top of molds. Accordingly, specific examples of this type of practice will not be considered in this paper. It appeared from this survey of deoxidation practice that semikilled steel can be divided into three types, namely, capped semikilled, in which aluminum is added near the end of pour or after shut-off, intermediate semikilled steels, in which aluminum is fed uniformly throughout the pour, and semikilled steels requiring no mold deoxidation. In an attempt to cover the ingot structures obtained by the range of practices reported in the preceding survey, three series of experimental heats were made using increasing amounts of aluminum, silicon and titanium for ladle deoxidizers. The ladle deoxidizers were added as stick aluminum, 50 pct ferrosilicon and 20 pct ferrotitanium (4 pct carbon), respectively. Two methods of studying ingot structures were used. Ingots were removed from production, burned longitudinally and machined to allow observation of the central pipe and coarse blowholes. During the machining operation much of the finely porous subsurface structure was obliterated. Accordingly, observations were also made

Jan 1, 1947

By J. L. LeBlanc, R. L. Lewis

This study presents an analysis of annular backpressure variations associated with controlled gas kicks and their pronounced effect on casing .strings and exposed under lying formations. A mathematical model describing the volumetric behavior of an extraneous gas as it is transported from reservoir to .surface conditions under changing temperatures and pressures has been programmed in a Kingston FORTRAN II language for digital computer analysis. The gases under investigation typify Gulf Coast reservoir gases within a 0.6 to 0.7 .specific gravity range. The program output has been substantiated by actual field cases. of gas kicks encountered in Gulf Coart we1l.s. The development of empirical equations for calculating suitable gamy deviation factors for unique temperatures and pressures was incorporated in the program to provide realistic solution.. An output listing of annular backpressures and corresponding equivalent fluid densities resulting at a predetermined critical depth (casing setting depth) and total depth for selected .stages of circulation is provided in a chronological .sequence. Additional information including reservoir pressure and temperature, kill rnrid density, produced gas or surface volume of the expanded gas intro vion, drill pipe and annular volumes can he obtained from the model. This paper illustrates that a precise knowledge of the volumetric behavior of extraneous gases in annular flow and its effect on equivalent fluid densities at a critical depth is significant and should receive .serious consideration in controlling threatened blowouts and in the design of drilling programs. Surface pressures in excess of formation limitations are a threat to zones of lost circula/ion and are potentially injurious to productive intervals. A knowledge of annular backpressure and equivalent fluid density profiles for probable gas kicks aids in a technological accomplishment of drilling programs and provides a .sale tolerance in the event a threatened blowout is encountered. Introduction Drilling operations are frequently interrupted when the drill bit penetrates permeable gas sands with reservoir CtfuJ manuscript was received in Society of Petroleum Engineers ofice Am. 1 1967. Revised manuscript received JuIy 7. 1968. Paper (SPE 1860) kae presented at SPE 42nd Annual Fall Meeting held in Houston. Tex., Oct 1-4, 1967. @ Copyright 1968 American Institute of Mining, Metallurgical, and Petroleum Engineem, Inc. pressures greater than that exerted by the drilling fluid. The differential pressures resulting permit an extraneous influx of gas into the wellbore. A suspension in drilling progress is necessary to restore fluid equilibrium throughout the system. Whether formation gas kicks originate unintentionally or by design, the prospect of a threatened or actual blowout exists and a method assuring a safe and effective well control procedure must be observed. A significant contribution to well control technology was advanced by Records et a1.l in 1962. Using the concept of transmitting a constant equivalent formation pressure at the point of intrusion, Records et al. introduced a calculation technique providing the annular backpressures encountered in a well control environment as a func tion of the volumetric behavior of a 0.6 specific gravity natural gas. In essence, the procedure outlined an annular backpressure schedule in terms of fluid volume circulated at different stages of a well control operation. A number of other publications2-' proposing various techniques for controlling gas intrusions in a wellbore achieve pressure control essentially through maintenance of a constant bottom-hole pressure by surface choke adjustments. The subsequent pressure effects induced in the annulus unfortunately receive little emphasis. Due to the tedious and repetitive nature of annular backpressure computations, a theoretical solution by digital computer is introduced for predicting annular backpressure and equivalent fluid density profiles associated with controlled gas kicks. We point out the effects of volumetric behavior of extraneous gases in annular flow and related field phenomena on equivalent fluid densities at a critical depth. The investigation indicated that equivalent fluid densities at a critical depth are of significance and should receive consideration in the control of threatened blowouts and in the design of drilling programs. Theoretical Considerations The mechanism of vertical gas flow through an annulus is governed by the PVT properties of the fluid, the pressure distribution within the system, the fluid flow rates and the geometry of flow. Due to the numerous variables involved in this type of problem, certain assumptions were imposed in deriving the mathematical model and in establishing the solutions. Two gases, characterized by specific gravities of 0.6 and 0.7, were selected to typify Gulf Coast reservoir fluids. The gas intrustion entered the wellbore as an immiscible 'References given at end of paper. JOURNAL OF PETROLEUM TECHNOLOGY

Jan 1, 1969

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By K. I. Savage, S. C. Sun

This paper describes a process developed to recover coal, clays and pyrite from coal wastes. The process consists of fine grinding followed by coal and pyrite flotation which leaves the clays in the flotation pulp. A bituminous coal refuse containing 10% sulfur and 30% carbonaceous material was treated by this method to yield a coal product containing 4% sulfur, 10% ash; a pyrite product containing 45% sulfur (84% FeS2), 1% carbonaceous material; and a clay product containing 2% sulfur (3.5% FeS2). The coal yield was about 89%. The pyrite yield was about 77%. The process steps may be entirely flotation, or gravity separation (hydrocycloning) may be used to increase the pyrite : coal ratio in the flotation feed. Cost estimates for the process show a profit of $2.28 per ton of low pyrite grade refuse, but these do not include labor, maintenance, overhead and plant depreciation. The development of this process consisted of three parts: (1) exploratory tests, (2) op-timazation tests and (3) confirmatory tests. The objectionable qualities that sulfur imparts to coal have been commented upon from early times, and they have become more objectionable as the uses of coal have grown. In whatever form coal takes — raw, carbonized or gasified —the sulfur content remains objectionable, and therefore its compounds are removed as completely as possible. It has long been known that sulfur occurs in coal in different forms. In 1861, sulfur was said to exist in the state of sulfuric acid in combination with a base; in combination with iron as iron pyrites; as bisulfides of iron; and in combination with the organic elements of coal.' Pyritic sulfur is a term loosely used to cover the sulfur associated with iron in its various forms. The mineralogically recognized forms are pyrite (FeS2), pyrrhotite (Fe 1-x S) and marcasite (FeS2). The particle sizes of pyrites vary widely. Isolated grains of marcasite smaller than 15 microns have been found disseminated through coal.2 Hair-like "veins" of pyrite filling cracks in vitrain have been found. At the other extreme, lumps or nodules of pyrites large enough for removal by hand picking have been encountered. Organic sulfur, unlike pyritic sulfur, does not exist as discrete particles, but is instead intimately associated with the coal structure and thus it is impossible to remove it or reduce its concentration by physical or mechanical means.2 In the preparation of coal for its various markets, the pyrite minerals (pyrites) are separated from the raw coal feed. This separation process concentrates the pyrites in the tailings or other waste products. It would be desirable to recover these pyrites for three reasons: (1) Pyrites are potential sources of sulfur and iron. In 1967, for the fifth consecutive year, Free World consumption of sulfur exceeded production.3 Propelled by the shortage, the domestic price of sulfur has risen from $24 to $38 per long ton (bright). (2) The refuse, when placed in piles, becomes ignited. The pyrites (FeS2) bum, giving off SO,. Thus, the pyrites are a cause of air pollution. (3) Also, the pyrites undergo chemical reaction when exposed to air. The refuse is leached by waters which result in stream pollution due to the water-soluble iron and acidic reaction products. Coal refuse also contains coal minerals and clay minerals. Therefore, any process for recovering the pyrites must successfully separate them from the coal and clay minerals. In the study discussed here, ten different bituminous coal refuse samples were successfully upgraded in pyrite content. These samples represented a wide variety of coal waste materials from Pennsylvania and other states. The variabilities of the sulfur and coal contents are shown in Table I. The extremes are Sample J, a high sulfur-low coal, and Sample B, a low sulfur-high coal. Definitions of some of the symbols or terms used in this report are given below: Mesh size—Tyler standard mesh screen sieves with an opening based on the square root of two. Fe —Indicates iron, as determined by a stannous reduction-dichromatic oxidation method. S—Indicates sulfur, as determined by the ASTM "Eschka" method for sulfur in coal.

Jan 1, 1969

By R. L. Vaughn, H. R. Cooper

Effective procedures for automating X-ray assaying by use of a computer are described. These methods have advanced the reliability and accuracy of X-ray data, resulting in improved assay measurements for operations and process control. Assaying by X-ray fluorescence has been increasingly used in the mining industry, both for analyzing solid H.R. COOPER and R.L. VAUGHN, Members AIME, are Systems Manager-Process Control and Senior Systems Analyst, respectively, Scientific & Engineering Computer Center, Ken-necott Copper Corp., Salt Lake City, Utah. TP 68B236. Manuscript, Sept. 16, 1968. Discussion of this paper, submitted in duplicate prior to Dec. 15, 1969, will appear in SME Transactions, March 1970, and AIME Transactions, 1970, Vol. 247. samples for smelter, mill, and mine quality control and for continuous stream analysis of mill pulps. The important advantages of X-ray methods in laboratory assaying are reduced analysis cost, and for on-stream analysis, immediate assay information on concentrator streams. Of the numerous reports concerning X-ray analysis in mining operations, one is referenced which can be noted particularly in conjunction with this report.' Continuous X-ray analyzers have been installed at Kennecott and other mining companies for operations guidance and trend information for a number of years. However, complete operational acceptance of X-ray data has been slow, due* to the difficulty in maintaining

Jan 1, 1970

By F. D. Patchen

Changes in the properties of partland cement upon the addition of fine-ground silica are discussed. Data were collected from formulations cured for periods up to 60 days at temperatures varying from 180" to 350°F. Compressive strengths, permeabilities and thickening times are compared. It is demonstrated that addition of find-ground silica in proper proportion results in a hydrated cement with superior properties above 230°F; i.e., early compressive strengths are higher and increase with additional curing, permeabilities are significantly lnrver and thickening times are increased. X-ray analysis of hydrated cement samples have been used to determine more fully the complex chemistry of these systems. Dicalcium silicate alpha-hydrate was found in all neat cements cured above 260°F. A 10 weight per cent silica addition produced a maximum amount of this phase. Xonotlite and tobermorite have been identified as hydration products at 320°F in systems containing 35 or more weight per cent silica. Formation of these components appears to be primarily responsible for the gain in strength and reduction of permeability observed after the addition of silica. INTRODUCTION For several years it has been recognized that loss of compressive strength occurs when portland cement is cured at a temperature above 250°F. This phenomenon has been called "strength retrogression". More recently it has been recognized that increased permeability usually is associated with decreased compressive strength at these temperatures. At curing temperatures above 300°F, these factors reach serious proportions. A compressive strength loss of 50 per cent between curing times of 24 hours and 14 days has been reported,' and permeabilities as high as 10 md have been measured in neat cements cured at 320°F for seven days.' Earlier investigators,' although unsuccessful in eliminating strength retrogression, stated ". . . .it appears possible that by the proper selection of additives for cement, strength loss with age at high temperatures can be either eliminated or considerably lessened". One solution suggested was the use of a pozzolanic material, hydrated lime, and a retarder. Menze15 in 1934 discovered that portland cement containing finely divided silica possessed high compressive strengths when the curing temperature was 350°F. The use of such a composition is now an established practice in the manufacture of steam-cured cement blocks. It is reasonable to assume that the same beneficial effect can be used advantageously in high-temperature oilwell cementing operations. Ludwig first mentioned the use of sand for this purpose. Later, a portland cement containing moderate additions of silica flour and retarded with CMHEC, a cellulose derivative. was recommended for high-temperature application." These formulations have been successfully employed in the field."'" This paper presents the results of a research program to test the efficacy of finely ground silica in increasing the compressive strength and reducing the permeability of oilwell cement compositions containing portland cement. Included also are X-ray diffraction patterns of powdered samples which have been used to determine some of the important chemical reactions occurring in these systems. LABORATORY PROCEDURES AND MATERIALS All compressive-strength test specimens were prepared and cured from a slurry in accordance with API RP 10B, with two modifications. A final curing pressure of 5,000 psi was maintained for a one- or three-day period; but, for further curing, pressures were reduced to the vapor pressure of water at the curing temperature. Also, the water content was varied in samples containing silica to maintain the slurry viscosity within a range of 5 to 15 poise as measured on a consistometer. Test specimens at 180°T were cured in a water bath at atmospheric pressure. At higher temperatures the specimens were cured in a steam-jacketed high-pressure curing vessel at 5,000 psi for a one- or three-day period. cubes cured for longer periods were removed from the pressure-curing vessel after 24 hours, sealed in small bombs containing water and placed in an aging oven.

By L. Cichowicz, K. K. Millheim

The constant-rate drawdown test performance for a low-permeability, verticany fractured gas well was investigated. A series of gar wells were tested by flowing each well at constant rate until the data could be analyzed using con-ventional radial flow theory. Each well was then shut in to build up. After a sufficient buildup was obtained, another flow Iest commenced but at a higher flow rate than Ihe first test. Again, the well was shut in when radial fiow was obtained. his procedure was repeated for three to four different flow rates. Two wells in the San Juan basin were tested using this procedure. Both wellere fractured after completion, cleaned up and then shut in until flow testing commenced. Test designs of both wells permitted investigation of the most realistic values of egective permeability, wellbore radius and turbulence factor. Also, being able to determine the eflective fracture flow area and vertical fraclure efficiency was inherent with this testing approach. It was observed that fractures in both wells influenced the Pressure behavior for approximately Is to 40 hours (depending on the flow rate) before radial flow was evident. After this time, drawdown data were analyzed using radial pow theory. When a low-permeability gas well has vertitally oriented, induced fractures, the early flow geometry ic. essentially linear. It will be shown how to determine when a flow test has been conducted long enough so lha' the most representative values of effective permeahiiity. wellbore radius and turbulence factw can be calculated. From the linear pressure data, valuable information about The fracture treatment, such as the effective flow area and vertical fracture efficiency, can be determined for vertically froctured wells. Introduction During tests on gas wells in the Sari Juan basin7 initial transient behavior lasted for many days because of the low permeability of some porous media. As a result, stabilized flow performance could not be obtained. If these wells received some type of stimulation treatment, early pressure behavior deviated from conventional theoretical radial Row. When conventional radial flow theory was used to analyze these low-permeability fractured gas wells, larger values of flow capacity and absolute open flow potentials (AOF) sometimes resulted. Wells were assigned open flow poten- tials that proved to be 3 to 10 times higher than the well would sustain over a longer period of production. In some cases where the wells had flowed for longer periods of time during a constant-rate drawdown test, it was noticed that the effective flow capacity appeared to be decreasing with time until a certain value was reached. The early nonradial pressure behavior can be explained If linear flow is assumed. Rusell and Truitt mathematic. ally investigated the vertically fractured well in a bounded area. They showed that early flow behavior was essentially linear and, for x,/x. approximately less than 0.10 radial flow was obtained after short periods of time. Then realistic values of effective permeability and skin could be determined. Scotta experimentaliy studied the vertically fractured well with a heat flow analog. He showed that earb flow was linear, Both studies indicate that, for small values of x,/x,, linear flow approaches radial flow if the well is tested long enough, To help prove this concept of early linear flow caused by induced vertical fractures, two low-permeability gas wells were tested. Both wells received large fracture treatments prior to testing. A vertical fracture was indicated from the analysis of fracture treatments, As anticipated, tests of both wells indicated early linear flow that was later followed by a period of radial flow, Data collected from each well were analyzed. From the well tests, plus other information on each well, the effective permeability, wellbore radius and turbulence factor were calculated. Effective fracture flow areas calculated from test analyses for each well proved to be approximately one-fourth the created area calculated from classic hydraulic fracturing theory: other fractured wells that were tested but not presented in this paper also indicated that the effective fracture flow area was one-fourth to one-third the created area predicted from hydraulic fracturing theory. The vertical fracturing efficiency was estimated from the calculated values of effective wellbore radius and fracture flow area. For the two wells tested, calculated fracture lengths x, were 112 and 105 ft, and the vertical fracturing efficiencies E, were 122 and 183 percent. Development of Flow Model Agnew' showed that most induced fractures below 1,500 ft are vertical. Anderson and Stah15 indicated that most of the fractures they studied were vertical. Thc model proposed for early flow in most vertically fractured gas wells is shown by Fig. 1. This model should approximate early flow behavior until radial flow is reached, at which time a radial model with an effective wellbore radius of 0.5 they will

Jan 1, 1969

By R. A. Wolfe, H. W. Paxton

Self-diffilsion coefficients for cr51 and Fe55 in 12 pct Cr-Fe and 17 pct Cr-Fe for Fe55 in chromium, and for Cr51 in vanadium have been measured. The results are compared with other values for the Fe-Cr system, and with the various theories of diffusion in hcc metals. Some empirical correlations are discussed between Do and Q in hcc systems, or, expressed differently, the constancy of ?G*/T solidus for seveval bcc metals and alloys is noted. It appears very probable that a vacancy mechanism is operative in bcc metals, hut this cannot he stated with certainty. THE great bulk of work on diffusion in metals, both experimental and theoretical, was for many years concentrated on those with close-packed and, in particular, fcc lattices.1,2 There appears to be little doubt that the mechanism of diffusion in these solids is vacancy migration, leading to mass transfer and in substitutional solid solutions to a Kirken-dall effect.3,4 For bcc metals, the picture is much less clear. The Kirkendall effect certainly occurs in several alloys.5-10 However, attempts to understand the factors contributing to the pre-exponential in the usual expression for the diffusion coefficient D =D, exp {-Q/RT) by extension of ideas useful in close-packed lattices have not always been successful. Zener,11 Leclaire,12 and Pound, Paxton, and Bitlerl3 have suggested that various forms of ring diffusion may be important in some bcc metals. For close-packed metals, Do is usually about 1 sq cm per sec and Q - 35Tm kcal per mole (Tm = melting temperature in OK). The theory of Pound et al. suggests for ring diffusion that Do may be about 10-4 and Q, although difficult to calculate with any precision, would be significantly less than 35 T,. The experimental results on self and solute diffusion in ? uranium14,15 and ß zirconium,10 and for solutes in 0 titanium,17 and possibly for self-diffu- sion in chromium below about 0.75 T,," gave some credence to this theory. However, not all bcc materials display low values of DO and Q, and the exceptions were not predicted by any theory. Furthermore, it has recently become apparent that, in bcc materials, log D is not always linear with T-l if a sufficiently wide range of temperature is studied.16,18 This variation may be such that Q may increase18,19 or decrease20 with increasing temperature. The present work was undertaken in an attempt to provide further diffusion data on bcc metals, and to try to understand the factors which contribute to differences in behavior between the various elements. For part of this work, the Fe-Cr system was chosen since it is of considerable technological importance, and data on 12 pct Cr and 17 pct Cr alloys appeared well worthwhile to supplement that existing for the remainder of the stern.18,22 The diffusion of Fe55 in chromium was studied as an example of a more or less "normal" tracer element in a possibly abnormal host lattice. Finally, no data were available for vanadium, the neighbor of chromium in the periodic table, because of lack of a suitable isotope so cr55 was used as a tracer in a few preliminary experiments. For convenience, we shall refer to elements whose Do and Q are low compared to those predicted by Zener's theory as "anomalous". PROCEDURE This investigation determined self-diffusion rates by means of radioactive tracers and the integral-activity method first utilized by Gruzin.23 In this method a thin layer of radioisotope of the diffusing element is plated or coated onto a planar surface of the diffusion sample, which is then given an isothermal-diffusion annealing treatment. The determination of an activity-penetration curve involves measuring the residual activity of the specimen after each successive layer or section has been removed parallel to the original planar surface. The method used here is essentially the same as that used by Gondolf18 and Kunitake.21 Two radioactive tracers, cr51 and Fe55, were used in this investigation. Diffusion coefficients were determined for the diffusion of one or both of these tracers in four different materials, viz., Fe-12 wt pct Cr alloy, Fe-17 wt pct Cr alloy, chromium, and vanadium. The diffusion samples had nominal dimensions of 1.5 cm diameter and 0.5 cm thickness. The grain size was several millimeters for the Fe-Cr alloys and at least 1 mm for the chromium and vanadium samples. Accurately planar surfaces

Jan 1, 1964

By G. Derge, Ling Yang, G. M. Pound

The specific conductance and its temperature dependence were measured over the entire composition range of the molten Cu2S-CuCI system. At a typical temperature of 1200°C, 10 rnol pet of the ionically conducting CuCl reduced the specific conductance from about 77 ohm-lcm-l for pure Cu2S to about 32 ohm -1cm -1, and 50 mol pet CuCl reduced the conductance to that for pure CuCI—about 5 ohm 1cm1. The nature of electrical conduction in molten Cu2S, FeS, CuCI, and mixtures was studied by measuring the current efficiency of electrolysis at about 1100°C. The Cu2S, FeS, and mattes were found to conduct exclusively by electrons, but addition of 1 5 wt pet CUS to Cu2S produces a small amount of electrolysis. Addition of CuCl to Cu2S suppresses electronic conduction, and ionic conduction reaches almost 100 pet at a CuCl concentration of about 50 mol pet. These facts are interpreted in terms of electron energy level diagrams by analogy to the situation in solids. RESULTS of electrical conductivity studies on molten Cu-FeS mattes as a function of composition and temperature have been reported.' The specific conductances ranged from about 100 ohm-' cm-' for pure Cu2S to 1500 ohm-' cm-1 for pure FeS. This is in sharp contrast with the low specific conductance of molten ionic salts for which the transfer of electricity is by migration of ions in the field. In general, these ionically conducting molten salts, such as NaC1, KC1, CuC1, etc., have a specific conductance of the order of magnitude of 5 ohm-' cm-'. It was concluded on the basis of this evidence that molten FeS and Cu,S exhibit electronic conduction. Pure molten FeS has a small negative temperature coefficient of specific conductance, resembling metallic conduction, while pure molten Cu2S has a small positive temperature coefficient, resembling semi-conduction. The molten Cu2S-FeS mattes follow a roughly additive rule of mixtures, both with respect to specific conductance and temperature coefficient. Savelsberg2 has studied the electrolysis of molten Cu2S and Cu2S + FeS. He concluded that while molten Cu2S is an electronic conductor, there is some ionic conduction in molten Cu2S + FeS3 owing to the formation of the molecular compound 2Cu2S.FeS and its dissociation into Cu1 and FeS2-1 ions. The present work does not verify his results. Chipman, Inouye, and Tomlinson" have studied the specific conductance of molten FeO and report a high specific conductance, about 200 ohm-1 cm-1 of the same order of magnitude as that found for molten mattes, and a positive temperature coefficient. They interpret these results in terms of p-type semiconduction in the ionic liquid by analogy to the situation in solid FeO.1 imnad and Derne' detected appreciable ionization in molten FeO by means of electrolytic cell efficiency measurements. In order to verify the conclusion that electrical conduction in molten Cu2S and mattes is electronic, and to gain further insight into the structure of molten sulfides, the following investigations were carried out in the present work: 1) The specific conductance, s of the molten system Cu2S-CuC1 was measured as a function of temperature over the entire composition range. As discussed later, molten CuCl is an ionic substance. It was thought that if molten Cu2S were simply ionic in nature, addition of small amounts of CuCl might not have a catastrophic effect in lowering the high conductance of the Cu2S. On the other hand, if much electronic conduction occurs, addition of the ionic CuCl should have a large effect in destroying the electronic conduction. 2) The electrolytic cell efficiency of the following molten systems was measured at about 1100°C in specially designed cells: Cu3; Cu2S + FeS, 50:50 by wt; FeS; Cu2S + CuS, 15 wt pet; Cu2S + CuC1, 5.9 to 46.4 mol pet; and CuC1. This gives a direct measure of the fraction of current carried by ions in these melts. Further, the cell efficiency, extrapolated to zero ionic current, is given by cell efficiency = (s leasile + s elexstronic). [1] s lucile for molten CulS would be expected to be no greater than that for molten CuC1, whose s lonle is about 5 ohm-' cm-1, as will be seen in the following. u,.,,.,.,.......for molten Cu,S is of the order of 100 ohm-' cm-'.' Thus, a large increase in cell efficiency from 0 to values of 10 to 100 pet upon addition of CuCl to Cu2S would indicate destruction of the electronic conductance. Conductance Measurements Experimental Procedure—The apparatus and experimental method were the same as those described in detail in connection with the study of electrical conduction in molten Cu,S-FeS mattes.' A four terminal conductivity cell and an ac poten-

Jan 1, 1957

By O. D. Gonzalez, R. A. Oriani

The migration of hydrogen and of deuterium dissolved in a iron and in nickel induced by an applied electrical potential has been measured over a range of temperature. In all cases the intevstitial solute migrates to the cathode. The deduced charges of transport are positive, independent of or slightly dependent on temperature, and markedly larger for the beavier mass. These results are interpreted as showing predominating momenlum transfer between electron holes and the acticated complex of the elementary diffusion act. THE migration of components of a metallic system under the action of an applied electrical potential is called electromigration or electratransport. This phenomenon has been known since 1861 and is becoming of increasing interest to metal scientists because of the opportunity it offers to explore the connections between mass transport and electron transport. Until recently, reliable data have been very scarce because of the difficulties inherent in the measurements. In addition, theoretical understanding has been either nonexistent or rudimentary and at present this situation is only moderately improved. Whereas the prevalent notion had been (and is still widely held) that electromigration directly reflects the state of ionization of a component in its equilibrium state in the metal, it is now becoming clear that the major, if not the sole, factor is the momentum transfer between the charge carriers, themselves in motion because of the applied field, and the metallic components in their activated state. It becomes important therefore to obtain reliable data over a range of temperature on electrornigration in systems where only one component moves, that is, either in pure metals (vacancy motion) or in interstitial solid solutions. Our choice of systems was dictated by our measurements' of the Soret effect in solid solutions of hydrogen and deuterium in iron and nickel. Indeed, the main reason for making measurements on electromigration at all is that we suspected a mechanistic connection between thermal diffusion and electromigration. In our former paper' we proposed that what is common to both phenomena is the momentum transfer between the charge carriers in the metal and the atom as it jumps from one equilibrium position to another. In Soret diffusion, however, there is in addition the dissipation of the energy of the jumping atom via pho-nons. We have since found that similar ideas for this interconnection have been developed by Fiks.2 In the present paper we will discuss in more detail our understanding of the mechanism of electromigration itself, but we turn first to a presentation of the experiments and results. EXPERIMENTS AND RESULTS Let us consider a metallic system in which only one component, m , can diffuse at the temperature in question, and to which are applied simultaneously an electrical potential gradient, grad cp, and a chemical potential gradient, grad p,. One may write3 for the material flux Jm and the electronic charge flux Je by applying the condition, grad = 0, to Eq. [I] and evaluating grad(um/T) for the isothermal case and for activity coefficient independent of composition over the relevant composition range. ere, c is the concentration and D the diffusivity of the mobile component. The experimental technique that we have used to measure Z* for the isotopes of hydrogen individually dissolved in each of two transition metals has been developed from a method used by Wagner and Heller.4 (A somewhat similar method was also used by Herold.5) A wire or thin rod of the metal under study connects two volumes of gaseous hydrogen at identical pressures; hence grad um = 0. A direct current is made to pass through the specimen, serving to maintain the specimen at the desired temperature and also to produce the electromigration of the hydrogen dissolved in the metal in equilibrium with the gaseous hydrogen. Under these conditions, the flux of hydrogen from one volume to the other is given, from Eq. [I] for grad( um/T) =0by

Jan 1, 1968

By R. O. Williams

Measurements have been made of the internal energy of deformation in a Cu-A1 alloy and a Cu-Zn alloy as the deference between the work and the released heat. The method required the rapid compression of samples to reduce thermal interactions. The rate of energy storage is low initially but then becomes almost linear with strain. Both alloys stored about 60 cal per mole at a strain of 0.4 in the recrystallized condition. The fractional rate of storage increases at low strains reaching about 0.25 at a strain of 0.15 and then decreases slowly with increasing strains. As do pure metals, these alloys release appreciable energy immediately following deformation; in contrast with pure metals, however, the release continued for extended periods of time. This energy release is somewhat sensitive to the amount of deformation, the deformation temperature, and prior thermal & eatment. It is believed that this release is a composite of that observed in pure metals which has been attributed to dislocations and that due to the motion of vacancies which produces increased short-range order. IT is expected that changes brought about by deformation of a pure metal can be understood in terms of the three defects—dislocations, vacancies, and interstitials—along with the interactions possible between these defects. Frequently stacking faults are created and must also be considered. 1t is not immediately clear whether or not the stress fields which are created should be regarded as independent of dislocations. For alloys there is the important additional effect of atomic configuration, for example short-range order, which must also be considered and will, in general, interact with each of the above defects. While it is true that a complete understanding of pure metals is not yet possible, it would seem very valuable to establish the similarities and differences between pure metals and alloys. Such knowledge could contribute to the understanding of both. A study of the changes in internal energy resulting from deformation of two alloys is reported here. Although some information has been available for alloys almost as long as for pure metals, the volume of work has been substantially less. Sato' and Quinney and Taylor' were responsible for the earlier work, both using brass and annealing calorimetry. tizhnova' studied stored energy in Cu-Ni alloys by a method somewhat similar to the one reported here. Bever and his co-workers have made extensive studies of Au-Ag alloys4'5 and more recently CU-AU'U' using solution calorimetry. Clarebrough et a1.' have recently studied brass using annealing calorimetry. EXPERIMENTAL METHOD The present results were obtained by the rapid, adiabatic compression of a sample about 1/4 in, diam by 1/2 in. long between two tungsten carbide hammers. A vacuum provided thermal isolation and oil provided effective lubrication between the sample and hammers. The increase in internal energy is the difference between the mechanical energy and the heat liberated within the sample. The mechanical energy was determined from the velocity of the hammers (calculated from the height through which the hammers traveled) and was essentially constant for each deformation cycle. The heat was determined from the temperature increase of the sample which was measured by an embedded thermocouple. Specific heat data were calculated from Neumann's rule using tabulated data for pure metals.9 Since the fraction of the energy which is stored is large compared to that for pure metals, specific heat data of high accuracy are not required. Successive increments of deformation were carried out by repeating the procedure, the limit being determined mostly by the decreasing accuracy caused by the increasing corrections and the decreasing storage of energy. About five cycles were used to give a true strain of about 0.5 while the time interval between cycles was 1 to 2hr.

Jan 1, 1963

By R. E. Cook

This paper presents results of a study to determine to what extent errors in estimated free gas saturation affect the results of static pressure calculations from build-up curves in two-phase systems. Use is made of the method of pressure build-up curve analysis developed by Miller, Dyes and Hutchinson1 , with modifications for two-phase flow advanced by Perrine2. Cornhined fluid compressibility of gas-oil mixtures is included, with special attention given to the discontinuity that occurs in fluid compressibility at the saturation pressure and its effect on static pressure calculations from build-up curves. The sensitivity of static pressure calculations to errors in gas saturation is considered over wide ranges of formation and fluid properties and reservoir pressure conditions. It is concluded from this study that for analyses of two-phase pressure build-up curves (I) the use of under-saturated crude compressibility data, even when producing GOR's are at or near solution values, will result in very significant negative errors in calculated static pressrrre, and (2) except for large values of buildup curve .slope and/or very low crude gravities, rather large errors in calculated gas saturation can he tolerated without .significant effect on calculated static pressures. The conc1usions arc qualified as referring only to the effects of gas saturation and excludc possible errors in the nuumerous other factors that enter into the ancllysis of pressure build-up curves. INTRODUCTION As part of the analysis of basic data for several major Eocene reservoirs in the Bolivar Coastal field (Western Venezuela), it was necessary to analyze a large number of pressure build-up curves for oil permeability and static reservoir pressure. From performance data it became apparent that gravity segregation is very active in these reservoirs. As a result of this performance it was recognized that the free gas saturation within the oil leg of the reservoirs would be small but would be subject to certain errors in calculation due to the inaccuracy of relative permeability data in the low gas saturation region. Since all of the various methods of build-up curve analysis were originally derived for sin- gle-phase, constant compressibility systems, and since the presence of a small but rather indeterminate gas saturation would greatly affect the reservoir fluid compressibility, a study was undertaken to investigate the effect of gas saturation on the results of two-phase pressure build-up analyses. TWO-PHASE BUILD-UP CURVE ANALYSIS The method of analysis used in this investigation is basically that originally developed by Miller, Dyes and Hutchinson1, with modifications for two-phase flow advanced by Perrine2. Miller, et al, showed that the original equation for single-phase permeability away from the well bore, held for two-phase flow except that the calculated permeability is the effective permeability to oil away from the wellbore. This finding was acknowledged by Per-rine, who also presented the equation for effective gas permeability* as Since Eqs. 1 and 2 contain no terms influenced by fluid compressibility, no distinct problem should develop in calculating permeability from two-phase build-up curves. In regard to the calculation of static reservoir pressure from build-up curves in single-phase systems, the modified Miller, et al, method as presented by Perrine is briefly as follows. 1. From build-up curve data, PVT data for the reservoir fluid and the well spacing and sand thickness, calculate the climensionless shut-in time corresponding to the highest value of actual shut-in time (t, hrs) on the straight-line portion of a semi-log build-up curve. Dimenrsionless shut-in time is related to actual shut-in time by the relation, t- 0.000264 kt 2. From a graph of dimension less shut-in time vs dimensionless pressure difference (Fig. 1 1 ) find the value of dimensionless pressure difference corresponding to the value of 7 calculated under (1) above. Dimen-

By R. G. Agarwal, Ramey Jr. H. J., Al-Hussainy R.

Although it has long been realized that gas recovery from a water-drive gas reservoir may be poor because of high residual saturations under water drive, it appears that only limited infomlation on the subject has been available until recently. This study was performed to show the qiiantitative potential importance of water influx. Results indicate that gas recovery may be very low in some cases: perhaps as low as 45 per cent of the initial gas in place. Gas recovery under water drive appear to depend in an important was on: (I) the prodirction rate and manner of production; (2) the residual gas saturation; (3) aquifer propertie.); and (4) the volumetric displacement effciency of water invading the gas reservoir. The manner of estimating water-drive gas reservoir recovery can vary considerably. Examples are: the steady-state tnethorl. the Hurst modified steady-state method, and various unsteady-state methods such ac. those of van Ever-dingen-Hurst, Hurst, and Carter-Tracy. The Carter-Tracy water influx expression was used in this study. In certain cases, it appears that gas recovery can be increased significantly by controlling the production rate and manner of production. For this reason, the potential importance of water influx in particular gas reservoirs should he investigated early to permit adequtrtr planning lo optirtize the pay reserves. INTRODUCTION In recent years, the economic importance of natural gas production has become increasingly apparent. This has been evidenced by more intensive exploration efforts aimed at gas production, and exploitation of both deep, as well as low-permeability gas reservoirs. Technical developments such as deep-penetration fracturing have made development of such formations economically feasible. Unfortunately, water influx has forced abandonment of a number of gas reservoirs at extraordinarily high pressures. Although reservoir engineering methods for estimating water influx have long been available, it appears that application of these methods to the water-drive gas reservoir has been sporadic.'a Available methods for estimating water influx which can be applied to the water-drive gas reservoir problem include the steady-state method,1 the Hurst modified steady-state method and various unsteady-state methods such as those of van Everdingen-Hurst. Hurst, and Carter-Tracy. Interesting applications of these solu- tions to gas reservoir and the aquifer gas-storage problems have appeared recently.3,12,14 The experimental study of residual gas saturations under water drive by Geffen et al. in 1952 indicated that residual gas saturations could be extremely high. A value of 35 per cent of pore volume is often used in field practice when specific information is not available. The study of Geffen et al. showed that residual gas saturation might be much higher in some cases. Naar and Henderson concluded that the residual non-wetting phase saturation under imbibition should be about half of the initial non-wetting phase saturation. The Naar and Henderson result that residual gas saturation under water influx should be about half the original gas saturation is recommended as an estimate if laboratory measurements are not available. Thus, it is clear that a considerable portion of the initial gas in place might be trapped in a water-drive gas reservoir as residual gas at high pressure. A full water-drive would result in loss of residual gas trapped at initial reser.voir pressure. Consideration of transient aquifer behavior leads to the conclusion that high-rate production of water-drive gas reservoirs could result in improved gas recovery by reduction of the abandonment pressure. However, there appears to be little quantitative information on this possibility. One of the few advantages of water-drive gas production appears to be improved deliverability through water-drive support of the reservoir pressure. There may also be an advantage in higher condensate recovery caused by pressure maintenance for gas-condensate water-drive reservoirs. In view of the preceding, this study was made to assess the potential importance of water-drive in gas reservoir engineering. The Carter-Tracy approximate water-influx expression was used because this equation offers some advantages in hand-calculation which do not appear to have been generally recognized.' However, calculations were performed in the main with a high-speed digital computer to permit evaluation of the effect of water-drive under a large variety of conditions. CALCULATION METHOD Water-drive gas reservoir performance can be estimated in a manner completely analogous to oil reservoir calculations: a materials balance is written for the reservoir, and a water influx equation is written for the aquifer. Siniltaneous solution provides the cun~ulative water influx and the reservoir pressure. When reservoir performance data (gas produced and reservoir pressures) are available, it is usually possible to match performance data to determine the initial gas in place and the water influx parameters —

Jan 1, 1966

By J. Chipman, N. J. Grant, J. H. Walsh, T. B. King

WITH the development of adequate sampling and analysis techniques, much information has been obtained concerning the behavior of hydrogen in the steel bath during the course of steelmak-ing operations. From studies of this type it has become clear that the slag has an important role to play in determining the hydrogen content found in commercial products. A knowledge of the behavior of hydrogen in slags is necessary to round out our knowledge of hydrogen in steelmaking practice; with this information we are better able to assess practical methods of control. In this paper a method developed for analyzing slags for hydrogen content, studies of the solubility of hydrogen in synthetic slags, and studies of the variation of hydrogen content in slag and steel in industrial heats will be described. Previous Work—While very little has been done in the past to determine the hydrogen content of metallurgical slags, measurements have been made on such slag-like materials as geologic magmas and glasses. Goranson as made extensive studies of the solubility of water in natural and synthetic granites. Many workers have investigated the solubility of gases in synthetic and commercial glasses, but the recent work of Russell" is the most extensive. Russell found that the solubility of water in glass was proportional to the square root of the water vapor pressure in water vapor-oxygen mixtures. It was found that the solubility could either increase or decrease with temperature depending on the glass composition. Wentrup, Fucke, and Rief,' and Piper, Hagedorn, and Backes' have extracted hydrogen by holding a slag at 1300" and 800°C respectively in vacuum. Wentrup et al. reported hydrogen and water vapor as evolved together while Piper et al. passed the evolved gases from the hot extraction treatment over hot ferromanganese to reduce the evolved water vapor to hydrogen and, therefore, reported the analyses in terms of total dissolved hydrogen. Yavoiskii has reported hydrogen and nitrogen analyses on steelmaking slags. Dobrokhotov, Povol-skii, and Khan have studied the hydrogen content of basic open hearth steel under slags of varying viscosity. These authors reported that the hydrogen content of the steel bath increased as the fluidity of the slag was increased by the addition of bauxite. They recommended that a viscous slag be used and the steel tapped without deoxidation in the furnace to ensure a low final hydrogen content. Analysis of Slags for Hydrogen Development of Apparatus—The vacuum fusion analysis equipment described by Shields, Chipman, and Grant" was modified for the analysis of slags." This method depends upon the reduction of water in the slag to hydrogen, which is readily determined by the method used for hydrogen in steel. Preliminary experiments were made with a eu-tectic ferrosilicon melt as the reducing agent. Other experiments included the use of graphite and molybdenum crucibles. The ferrosilicon alloy bath was found to produce a thick coating on the furnace tube which could adsorb hydrogen. Graphite produced CO by the reduction of FeO, MnO, and SiO2 in such quantities that analysis for hydrogen was virtually impossible. Molybdenum, while giving a low blank value, could reduce the slag sufficiently to produce a volatile oxide which again was found to adsorb hydrogen. The arrangement finally adopted is shown in Fig. 1. The molybdenum sheet served as the induction-heated susceptor, and was separated from the escaping gases by alumina guide tubes. An outer jacket of alumina protected the quartz furnace tube, and alumina and quartz stools completed the assembly. The falling sample is shown wrapped in aluminum foil, 0.002 in. thickness, which acted as the reducing agent' This foil was careful'~ washed in acetone and an d and then dried in hot air before wrapping the weighed slag sample for charging into the apparatus. Investigation of the hydrogen content of the aluminum foil itself showed that a small quantity of hydrogen was released in vacuum, but had no relation to the weight of the aluminum charged. This quantity of hydrogen was

Jan 1, 1957

By Charles S. Barrett

THIS paper is a progress report covering metallographic applications of the electron microscope that have been made during the past year at Carnegie Institute of Technology. An account is presented of the techniques that have been most satisfactory, the difficulties encountered, the resolution obtained, and the peculiarities that have been noted in the microstructures. About 12m exposures have been made during the year on the instrument, which is an RCA type B2 transmission microscope,1 and a number of these have been selected to illustrate the results obtained. Prints are also included showing the common types of falsities-to borrow a term from the biological field, the "artifacts." In brief, the paper is a record of our attempt to determine the value of the electron microscope in the study of the structure of metals and alloys. n addition to comments on techniques and the quality of photomicrographs obtainable with them, the paper discusses the nature of the etch attack on a few metals, and the possible relation between the fine structure revealed by etching and the mosaic block structure of the grains. Since numerous papers and reviews are available on the subject of imperfection in crystals, the comments here may be very brief; in fact, are restricted almost entirely to points having to do with recent theories relating the strength of metals to imperfections and to spacings of slip lines. SPECIMEN PREPARATION AND PHOTOGRAPHS Mechanical polishing by hand was used throughout this work. Precautions that must be taken in microscopy at high magnification, which have been well explained by Vilella in his book on Metallographic Technique for Steel, are likewise important with the electron microscope. This applies particularly to the need for repeated polishing and etching (each successive polish becoming less severe and each etch lighter), since the final etch should in general be no heavier than is used for oil-immersion optical microscopy and in many instances should be much lighter (compare Figs. 6a and 7a, 6d and 7b). An extremely light etch, of course, is incapable of cutting through a flowed layer of appreciable thickness. Unquestionably it would be very advantageous in each instance to compare the structures obtained by electrolytic polishing with those obtained by mechanical polishing, but difficulties with the electrolytic method have prevented this in the following collection of prints. The choice of etchant is important, of course; for example, picral proved superior to nital for tempered martensite. With the exception of Figs. Ib and IC, the replicas used were of silica, deposited up0n polystyrene and floated from the polystyrene in ethyl bromide according to polystyrene in ethyl bromide according to the method of Heidenreich and Peck.2 The samples were mounted in Bakelite, polished and etched, then the polished Bakelite

Jan 1, 1943

HUDSON MAXIM, Brooklyn, N. Y.-Immediately after the creation of the Naval Consulting Board, Admiral Strauss, Chief of the Bureau of Ordnance of the Navy, told members of the Board that the most important problem for the scientists to solve was erosion in large guns. The Navy, for years, had been working on the problem with marked suceess, but there was much yet remaining to be accomplished. As Chairman of the Committee on Ordnance and Explosives of the Naval Consulting Board, I communicated with several experts in the manufacture of steel, principal among whom was Prof. Howe. He generously agreed to help in every way possible, and the work that he has done and the time that he has spent upon it can he appreciated only by reading his paper. Prof. Howe asked for some rings. from a, certain 14-in. gun which Admiral Strauss had mentioned in particular, and in which, he had informed us, at each shot the bore of the gun, was enlarged about 0.001 in.,' and that after a relatively small number of rounds it was considerably eroded. It has puzzled scientists for a long time to account for the peculiar erosive effect produced by the powder gases in a gun, arid: what Prof. Howe has done in his paper is worthy of Sherlock Holmes. He is a veritable Sherlock Holmes of science.

Jan 4, 1918