Search Documents

By C. J. Guntner, R. P. Reed

A commercial 18Cr-8Ni iron alloy (AISI 304L) was examined in tension at 300°, 76°, 20°, and 4°K. Continuous stress-strain recordings were made, X-ray analyses at periodic stress (strain) intervals were obtained, and the magnetic measurements were taken. From this data the percentage of martensitic products [bcc(a) and hcp (E)] was computed as a function of stress (strain). It was found thatup to 15 pct E phase forms at low temperntures. The amount of E formed increases to a maximum at about 5 pct strain, then decreases. This decrease indicates the additional transformation of E to a'. The total amount of E and a' was suppressed at constant stress (strain) at 4°K as compared to 76°K. It is proposed that the suppression of E and a' is associated with the decreased mobility of extended dislocations at very low temperatures. The yield strength decreased as the temperature was depressed below room temperature and then increased rapidly near 4°K. SOME ferrous alloys which are austenitic (fcc ?) at room temperature appear to be unique in that two martensitic products (hcp e and bcc a') may form on cooling to lower temperatures or on application of mechanical stress. The most common room-temperature austenitic ferrous alloys are 18Cr, 8Ni stainless steels. Most aspects of the spontaneous transformations have been previously described for these steels.' Several previous papers have described special aspects of the stress-induced transformations at low temperatures for the stainless steels, such as the existence of the hcp phase (c) after straining at 76oK,2-7 the morphology after straining using electron microscopy,7 and the decrease in E at higher strains at 76oK.4 However, for a complete representation, one must know the stress-strain characteristics and the dependence of both martensitic products on applied stress and temperature. It is the intention of this paper to provide that documentation. To accomplish this, continuous stress vs strain recordings were made at four temperatures: 300°, 76", 20°, and 4°K for annealed AISI 304L (a commercial 18Cr-8Ni alloy). At periodic stress intervals at each temperature the integrated X-ray line intensity of a selected peak for each phase (y, E, and a') was measured. In addition, photomicrographs of the strained surfaces were taken and magnetic measurements were made. The magnetic readings can be directly converted into percent a'.',e With these measurements the percentage of each phase may be plotted as a function of stress (or strain) and test temperature. It was found that up to 15 pct E phase forms upon stressing the AISI 304L alloy at low temperatures. The E percentage increases abruptly after the alloy yields, but then decreases gradually at higher stresses. The rapid increase in e at 76°K is associated with an "easy-glide" portion of the stress-strain curve. The total amount of a' + .G is suppressed below 76°K at a constant stress or strain. The yield strength decreases down to 76°K but increases rapidly below 20°K. EXPERIMENTAL PROCEDURE Tensile test specimens were cut parallel to the rolling direction from 0.1-in.-thick sheet. Continuous stress vs strain recordings were obtained at each test temperature (300°, 76o, 20°, and 4°K) using equipment and methods described elsewhere.' The specimens which were used in the X-ray analysis were stressed to successive increments of strain at each temperature, analyzed at room temperature, then restressed at the test temperature. This procedure was repeated until approximately ten X-ray analyses had been performed with approximately 1.0 pct strain increments. The specimens had a reduced section 1 in. long, 1.2 in. wide, and 0.1 in. thick. They were electro polished prior to testing and after each strain increment. Table I lists the chemical composition, grain size, and hardness for the alloy which was used. This is the same alloy for which extensive mechanical-property tests3 and morphological studies of the spontaneous transformations' have previously been made. For the low-temperature tests (76o and 4°K) below the Ms temperature the specimens were initially cooled to the test temperature, held for 1/2 hr, then warmed and X-rayed at room temperature. The results are listed in Table 11. From earlier work8 it was known that additional transformation on the second cycle would be considerably less (-0.1 pct

Jan 1, 1964

By M. Mills, M. W. Clegg

This Paper Presents an approximate analytical solution to the problem of the pressure distributions arising from the production of a compressible liquid in a partially penetrating well. The limits of validity of the approximation are discussed and it is shown that the basic approximation is the same as that used in previous numerical of the problem. The variation in pressure at the production well with time for different penetration ratios is investigated and is shown to agree with previous work. The solutiOns are used to assess the effects of partial penetration at large distances from the well and at the wellbore itself. In the latter case the effects of anisotropy are also studied for constant, but unequal, horizontal and vertical permeabilities. It is 'uggested that measurements of the Pressure drawdown below the producing interval in a well should give some indication of a thickness-vertical permeability product for the formation; this would chmacterize the effectiveness of the formation below the producing interval. INTRODUCTION In many oil and gas reservoirs the producing wells are completed as partially penetrating wells; that is, only a portion of the pay zone is perforated. This may be done for a variety of reasons, but the most common one is to prevent or delay the intrusion of unwanted fluids into the wellbore. This paper studies the effects of partial penetration on the pressure distribution resulting from the production of a compressible liquid at a constant rate. To simplify the mathematical analysis the case of a single well located in an infinite reservoir of constant thickness is considered, and it is assumed that this well is perforated over an interval adjacent to the upper (or lower) impermeable boundary of the formation. A well producing from an arbitrary interval within the bulk of the formation could be treated in a similar way, but this action would merely add to the complexities of the solution and contribute little to the appreciation of the over-all effects of partial penetration. In a recently published paper by Odeh,10 the more general problem of producing from an arbitrary interval within the pay zone was studied for the quasi-steady state flow of a compressible liquid, This is the case in which there is a a potential at the outer boundary and a constant flux across the same boundary. This contrasts with the problem studied here, which is concerned with the nonsteady-state behavior in an infinite reservoir. The exact solution of the partial penetration problem presents great analytical problems because the boundary conditions that solutions of the partial differential equations must satisfy are mixed; i.e., on one of the boundaries the pressure is specified on one portion and the flux on the other. This difficulty occurs at the wellbore, for the flux over the nonproductive section of the well is zero, and the potential over the perforated interval must be constant. In the case of constant rate production from the well, this uniform potential is time dependent and unknown, and the additional condition that where hl is the thickness of the perforated interval, must also b'e satisfied. This problem may be overcome in the case of constant rate production by making the assumption that the flux into the well is uniform over the entire perforated interval, so that on the wellbore the flux is specified over the total formation thickness. This approximation naturally leads to an error in the solution since the potential (pressure) will not

Jan 1, 1970

By S. Nakajima, H. Okazaki

Quantitatiue studies of the deformation texture in drawn tungsten wives were made by the X-vay dif-fractottletetr. Experimental results show that the diffraction Intensities are equal to tilose pvedicted from the (1 10). fiber lexlure but the angxla), spreads of. diffraction peaks in the pole distribution curres are different for different diffraction planes and directions. For this reason a modified (110) fiber lextuve model, in which a kind of anisotropy is assumed, is proposed to explain the results. According to this model the poles lying on a line directing front the (110) to the (110) poles in the (1 10) standard stereograpllic projection should show spreads which are different from those lyitlg on a line directing from the (001) to the (001) poles, which is confirmed by the experiments. The anisolvopy and the spveads of the pole positions are large at the outer part of the wires and decrease gradually lowards the inside of the wire. The possibilily of occurrence of such anisolropy in irrelals with fcc stvuctures is discltssed. THE deformation texture of drawn tungsten wires has been assumed by different investigators to be the simple ( 110) fiber texture.' Recently, however, Leber2,3 has shown that a swaged tungsten rod has a cylindrical texture. It changes gradually to the (110) fiber texture by drawing through dies. However, even after drawing to 0.25 mm in diam, the cylindrical texture can still be found in wires together with the (110) fiber texture. This was deduced from the pole figures obtained from the longitudinal section of these wires. Use was made also of quantitative measurements of the pole distribution curves. Leber stated that the angular spread of the pole distribution curves (henceforward called dispersions) are quite different for (400) 45 deg and (400) 90 deg: the former is always larger than the latter. This inequality is accompanied by deviations of the diffraction intensities from the theoretical values for the ( 110) fiber texture. Bhandary and cullity4 have reported similar results on iron wire and explained them by assuming a cylindrical texture. Both Leber3 and Bhandary4 used only the results of the (400) reflection for the determination of the dispersion. The pole figures found by Leber3 and by Rieck5 are largely different. The model given by Leber to explain the effects is in the authors' opinion in some respects unsatisfactory, especially if one looks at other than the (400) reflections. Intensities and dispersions of diffraction peaks are conclusive factors for the determination of the fine structure in wire textures. For this reason we studied them extensively to come to a model which is more suitable to fit the facts. In the following, after giving the experimental set-up, we report about measurements of X-ray diffraction on drawn tungsten wires. Different models to describe the experimental results will be discussed. EXPERIMENTAL GO-SiO2-A12O3 doped tungsten wires drawn to 0.18 mm in diam were used for the measurements. The wires were chemically etched to various diameters down to 0.03 mm. Measurements were carried out for the different wires in order to determine the dependence of the texture on the radius. The wires were cut to pieces of 10 mm length and fixed with paste closely against each other on a flat, polished glass plate. Parallelism of the wires with the surface of the glass plate should be adequate. For the diffraction studies three different X-ray sources were applied, respectively, giving the CuK,, FeK,, and FeKp emission. The measurements were carried out with a diffrac-tometer with a GM counter. The latter was fixed to a certain diffraction angle 20hkl and the diffraction intensity was recorded as a function of the angle of rotation of the specimen around the axis, lying in the specimen surface and perpendicular to the wire axis, as shown in Fig. 1. Measurements were also done with the detector at angles slightly deviating from the diffraction maxima The measured intensities in this case were taken to be equal to the background level. The deviations were chosen as small as possible but large enough to eliminate the influence of the diffraction maxima. The useful range of the rotation angle x of the specimen is generally limited by the wavelength of the X-rays. We have: where and cp is the angle between the wire axis and the normal of the diffraction plane. Intensity measurements were made to find the necessary corrections for counting loss of the GM counter and for distortion resulting from such effects as absorption of X-rays and from inclination of the reflection plane under study with respect to the surface of the specimen. The counting loss was esti-

Jan 1, 1968

By C. E. Taylor, L. V. Fegan

With on-stream electrical-resistance measurement as a basis, Bethlehem Steel developed a system for the completely automatic control of moisture in the feed to pelletizing. This system was later successfully adapted to the sintering process, with automatic control of moisture serving to maintain the air permeability required for effective burnthrough. Control of moisture content is an important factor in the agglomeration processes by which iron ores are converted into either pellets or sinter for use in blast furnace burdens. In both pelletizing and sintering the iron ore feed must be mixed with the proper amount of water so that capillary attraction and surface tension forces can be utilized to cause the ore particles to coalesce and form larger agglomerates. The success or failure of the entire pelletizing or sintering process depends on maintaining the proper moisture. Therefore, continuous automatic control of moisture is one of the conditions for the automation and optimum use of both pelletizing and sintering processes. In pelletizing, variation in moisture can cause ball size to change and is a common cause of cycling which results in a fluctuating feed to induration. Excessive moisture can cause lumps to form at the cutter bar and contributes to the formation of buildup on the green pellet screens. Green pellets with excess moisture deform and reduce bed permeability. In sintering, accurate control of moisture is also important because moisture must be within the limits that give the bed the required permeability for good burnthrough. In response to these problems, a moisture meter, based on the principle of electrical resistance measurement, was developed. This moisture meter has been a key factor in the development of systems for providing uniform feed of concentrates in most of Bethlehem's pelletizing and sintering plants. The automatic moisture control system for pelletizing was successfully plant-tested and resulted in the following: (1) Control of moisture in the green pellets to +0.2%. (2) Up to 50% reduction of bentonite consumption. (3) Control of unscreened green pellets so that 85% of the furnace product at Cornwall was within the size range of -5/8 +3/8 in. (4) Improved quality of product with a tumble test of 94 and a "Q" index of 90 as the result of feeding a consistently uniform feed to a shaft furnace. In sintering, where sinter-mix temperature change and variable feed could cause errors in moisture measurement, the system was designed to maintain the correct moisture under virtually all plant operating conditions. It was successfully plant-tested and is now in use in most of Bethlehem's sintering plants. The system achieves: (1) Control of moisture within the desired limits for optimum production. (2) Continuous measurement and control of moisture even with intermittent feeding of mix, so that attention by the operator is seldom required. The two automatic control systems provide smoother startup, minimize upsets and contribute to the consistent control of process variables in pelletizing and sintering. MOISTURE-CONTROL SYSTEM FOR PELLETIZING Fig. 1 shows the general arrangement of the system for automatic control of moisture in pelletizing.1 In this system a bin level controller, using a single probe which regulates filter speed, maintains a uniform low concentrate level in the bin.' The moisture meter, comprised of an analyzer and control components, continuously measures the moisture in the bed and adds water to adjust to the desired level. The meter is placed in an open-loop (feed-forward) control, i.e., water addition occurs downstream from the probes. Filter vacuum must be set so as to produce a filtercake moisture which is below the optimum for pelletizing. Bentonite is added after the

Jan 1, 1969

By N. A. D. Parlee

The diffusion rate of hydrogen in liquid iron has been measured by a gas-liquid metal diffusion cell technique. The diffusion cell was formed by immersing an alumina tube containing hydrogen gas at 1 atm in a bath of stagnant liquid iron. The change in the composition of the melt in the cell was determined by measuring the rate of absorption of the gas in the cell. The appropriate solution to Fick's second law was used to examine the data and calculate diffusivi-ties. The absorption of hydrogen in stagnant pure liquid iron has been found to be diffusion-controlled. The results show that the chemical diffusion coefficient, D, of hydrogen in pure iron in the range of 1547" to 1726°C can be represented by the following Arrhenius relation: D(sq cnz per sec) = 3.2 x X exp(- 3300 i 1800/RT) where the uncertainty in the activation energy corresponds to the YO pct confidence level. Oxygen in the melt (above 0.015 pct 2) increased the apparent rate of absorption of hydrogen. The importance of diffusion data on liquid metals for predicting the rates of certain metallurgical processes has been recognized for a long time. Moreover, these data are much needed to test and develop theory for diffusion in liquid metals. Despite this practical and theoretical interest, however, relatively little reliable information about diffusion in liquid metals is available in the literature. This is particularly true for gas components such as hydrogen, oxygen, and nitrogen in liquid metals, where almost no data on chemical diffusion coefficients are to be found. This is probably due to a multitude of experimental difficulties particularly associated with high-temperature melts. In an effort to fill this gap in information, a research program was undertaken to study the diffusivities and rates of solution of gases in liquid metals. This paper presents the results of a study of the diffusion of hydrogen in liquid iron. EXPERIMENTAL METHOD Two methods for the study of the kinetics of dissolution of gases in liquid metals are being employed in this laboratory. Both involve the measurement of the volume of gas absorbed by the melt as a function of time and as such both avoid the uncertainties involved in chemical analyses of quenched samples for relatively small amounts of gas. In the first method, the gas dissolves in an inductively stirred melt and, in the absence of a slow surface reaction, the results are often interpreted in terms of mass transport across a liquid "boundary layer" between the homogeneous gas phase and well-stirred part of the melt. Other interpretations of the results of such experiments have also been described in the literature.1'5 In the second method a gas-liquid metal diffusion cell is used.' The gas dissolves in a cylindrical column of stagnant liquid metal and, in the absence of a slow surface reaction, the results are interpreted in terms of a non-steady-state diffusion solution to Fick's second law. The weakness of the first method is that while it gives information on the mechanism of absorption by stirred melts it yields an overall rate constant which even in the simplest cases depends on the nature and the thickness of the "mass transport layer" or "boundary layer". It yields no values of diffusion coefficients. The second method was used in this research because in many cases it is possible to determine the diffusion coefficient of the gas component in the liquid metal. In this research it has been utilized to measure diffusion coefficients of hydrogen in liquid iron. The apparatus used was essentially the same as that described by Mizikar, Grace, and par lee but certain modifications have been introduced to meet the elevated temperatures and special conditions of this research. Fig. 1 is a schematic drawing of the apparatus and Table I gives the identification of various parts in this figure. The diffusion cell, shown in detail in Fig. 2, was formed by immersing an impervious alumina tube (hereafter called absorption tube) in a bath of pure liquid iron contained in an alumina crucible. Two types of tubes were used, Morganite triangle RR and McDanel AP35. The crucible was contained in a vertical impervious alumina combustion tube (32 mm ID by 914 mm long) which was closed at both ends by water-cooled brass heads employing O-ring compression seals, Fig. 1. A protection tube enclosing a Pt, 5 pct Rh-Pt, 20 pct Rh thermocouple was introduced through the lower end of the combustion tube

Jan 1, 1968

By M. P. Tixier

In several log interpretation methods, water saturation is evaluated by use of the ratio of the readings of a short spacing resistivity device and a long spacing resistivity device plus information on the mud filtrate and connate water resistivities, which is often derivable from the SP. These methods are valid over a certain range of conditions, usually specified by moderate invasion and R.. greater than R,. Since these methods involve no explicit evaluation of the formation factor, F, the saturafiion so found may be used in the standard Archie equation to derive a computed formation factor which, as a check, may be compared with the formation factor known from other independent measurements. Discrepancies in the two values of formation factor generally indicate that the ratio method is not within its range of applicability. If the computed F is too low, the corresponding S, is too low; if the computed F is too great, the corresponding S, is too great. By means of this "porosity balance" check, and some knowledge of the probable conditions of invasion, the interpretation can often be improved. The poqosity balance check is discussed for the cases of the Induction-Electrical Log Interpretation Method, the Rocky Mountain Method, and the R.,/R, Method. For the first method, a discussion is also given for the case of shaly sands. INTRODUCTION Al though many logging tools are available to obtain a practical value of porosity, the only logging method reflecting the water saturation is electrical logging. Water saturation is commonly computed from electrical log data by means of the basic Archie formula, Where F is the formation factor, R, the resistivity of the formation water, and R, the true resistivity of the formation. One difficulty in this procedure is a possible error in the evaluation of R, from the logs— for example, in the case of very deep invasion. There is no possibility, inherent in the procedure, which would provide a check of the accuracy of R, and, hence, of the computed saturation. Other approaches have been proposed which are based essentially on the ratio of the readings made with two dzerent devices: one device, with a short radius of investigation, provides a value close to the resistivity of the invaded zone; the second device, with a long radius of investigation, gives a value close to the true resistivity of the bed. By these methods, the determination of saturation is made independently of any direct knowledge of the formation factor. Therefore, if the formation factor (or porosity) is known from other sources (core analysis, neutron log, MicroLog . . .), it becomes possible to take advantage of this additional information in order to check the value of saturation derived from these latter methods. Inasmuch as this proposed check procedure relies on a knowledge of the value of formation factor or porosity, it has been deemed appropriate to call this verification a porosity balance. The purpose of this paper is to discuss the principle and the application of the porosity balance, depending on the method of interpretation being used. Basis of the Method There are in existence three interpretation methods by means of which the water saturation can bc found without a knowledge of the formation factor (or porosity). As just stated, they are based on the fact that the ratio of the two resistivities (one measuring the resistivity close to the hole, the other measuring the resistivity far from the hole) is independent of the formation factor since it affects the two resistivity measurements in the same way. Chronologically, the resistivity ratio methods are: (1) the Rocky Mountain Method, which has been used since 1943 and was explained in a paper in 1949;' (2) the R../R, Method or Shaly Sand Method, presented in 1954, which permits the calculation of the water saturation in both clean and shaly sands2; and (3) the Induction-Electrical Log Method, presented in 1957 for the interpretation of the induction-electrical log combination." Each method requires certain definite conditions in order for the interpretation to be reliable. The porosity balance method will let us know whether these conditions are fulfilled. In clean formations we know that the fundamental Archie relation (Eq. 1) must always be verified. This relation can be rewritten:

By J. Larsen-Badse

The resistance to dry sliding abrasion us been measured for some hardened and tempered carbon steels. The wear resistance, defined as the reciprocal of the volume loss per unit sliding distance, is found to vary linearly with the hardness of the abraded surface. The data fit two lines, one, of high slope, for the pure metals and for the annealed steels and another, of smaller slope, for the hardened and tempered carbon steels. A relationship is developed showing the alwasion resistance to be proportional to the hardness of the material and to en, where n is the exponent of work hardening. The values of n and of the hardness FOR most materials a high resistance to abrasion is generally associated with a high value of indentation hardness. Brine11 even used the abrasion resistance as a measure of hardness in an attempt at expanding the hardness range covered by his indentation-hardness test.' When comparing different types of materials there is, however, only a very general relation between abrasion resistance and hardness. Sample brittleness seems to contribute to the deviation from a direct relationship. In general, increased brittle- are those relating to the heavily deformed abraded surface. However, for materials having n constant or varying m a regular manner over a wide range of strailts, the wear resistance will also be a linear function of the bulk hardness of the material. For tempered steels containing a dispersion of cementite in ferrite the increase in wear resistance due to the dispersion is found to be a linear function of the volume fraction of Fe3C and of the reciprocal of the mean free ferrite path. This indicates indirectly that the yielding of this material, when it is already heavily deformed, follows the Orowan criterion. ness leads to decreased abrasion resistance.' When comparing materials in the same group, e.g., the metals, the spread is less. Kruschov found that the abrasion resistance of pure metals and of annealed steels is directly proportional to the hardness of the metal. For hardened and tempered steels he obtained another linear relationship, one of smaller slope and intersecting the line for the pure metals at the location of the data for the annealed steels. He also found that the abrasion resistance of cold-worked metals is independent of the degree of cold working. This shows, as he concluded, that the process of abrasion work hardens the metal in the surface region to a very high degree.

Jan 1, 1967

By DeVerle P. Harris

Risk for most mineral investments is greater than that of average business ventures, because the physical and economic characteristics of the mineral deposit are never known with certainty. Inasmuch as the merits of a potential mineral investment cannot be fully appraised without consideration of the impact to the investor of the possible unfavorable events associated with the adoption of the project, methods of risk analysis are important in project evaluation and investment selection. The conceptual framework for risk analysis and the means by which some of the recommended practices might be implemented are examined. In addition to the normal business risk due to market forces, the physical and economic characteristics of the mineral deposit are never known with certainty, for only partial or indirect information is available, such as data from drill cores or surface and underground workings. Consequently, risk for most mineral investments is greater than that of average business ventures. Inasmuch as the merits of a potential mineral investment cannot be fully appraised without consideration of the impact to the investor of the possible unfavorable events associated with adoption of the project, methods of risk analysis are important in project evaluation and investment selection. This paper examines the conceptual framework for risk analysis and means by which might be implemented some of the recommended procedures. Utility Theory A discourse on risk and its impact upon investment selection quite often leads the reader into the concept of utility theory, for utility theory is a device constructed by economists to measure the satisfaction an individuai receives from a good; in the case of investment analysis, that good is the outcome of an investment. It briefly is reviewed here to describe the reaction of the individual to risk and as a stepping stone to other techniques more suited to risk analysis at the level of the firm. Utility theory proposes that an investor possess a utility map (curve), such as indicated in Fig. 1, which indicates the value in utiles of potential monetary rewards to an individual with some specific endowment of resources. Points on this map are obtained from the following relationship:

Jan 1, 1971

By S. K. Tarby, W. O. Philbrook

Limited experimental data and a critical review of the literature are given to indicate that the true equilibrium distribution of manganese between carbon-saturated iron and blast-furnace type slags has not yet been satisfactorily established as a function of the composition, temperature, and pressure of the system. THE distribution of manganese between carbon-saturated iron and blast-furnace type slags has been studied by a number of investigators.1-7 Generalization of the results of these works indicates that the weight percent distribution ratio, (%Mno)/[%Mn], decreases with increasing temperature and slag basicity, in accord with thermodynamic data and practical observation. However, where comparisons among these studies can be made for similar slag compositions and temperatures, the reported values of the manganese distribution ratio are discordant. The variations are attributed to failure of the experimenters to bring their systems completely to equilibrium with graphite and CO at 1 atm, as judged by departures from silicon equilibrium. Even in the most recent work,5-7 this problem was recognized but the authors made the questionable assumption that equilibrium could prevail for manganese distribution even though other reactions were still in progress. Manganese equilibrium values for three particular slag compositions were required incidental to a kinetic study8 and were determined experimentally. Special pains were taken to approach concurrent equilibrium with respect to silicon. Although the results were not completely definitive, they indicate that the true manganese equilibrium distribution ratios are lower than have heretofore been reported. The data are presented together with a brief commentary on the literature to suggest the need for further work in this area. EXPERIMENTAL METHOD Equilibrium runs were carried out in a stationary, induction-heated graphite crucible closed by a graphite lid and flushed continuously by a slow flow of carbon monoxide at 1 atm pressure. Temperature was controlled manually within ±3 deg C with the aid of immersed thermocouples in mullite tubes protected against slag attack by graphite sleeves. The desired temperature was established by a calibrated Pt vs Pt-10 pct Rh thermocouple and control was then maintained using a more stable W-Mo couple. Runs were made at 1500º and 1575ºC for durations of 7.5 to 10.5 hr. The charges were calculated to yield 400 g of metal and 100 g of slag. The primary slag compositions and designation numbers were as follows:

Jan 1, 1963

By Floyd E. Bliven

A century ago, Thomas Carlyle sketched the Emblem of the Wasting Candle in one of his papers. "Terar Dum Prosim" he wrote under the candle-" May I be wasted so that I be of use." The application of his reasoning is an excellent guide for a discussion of the scrap problems of today. The General Electric Co. has addressed itself to the need of proper and careful preparation of scrap metals for many years. A supervisor of salvage directs this work in each of the works. Generally he is an engineer and has had broad factory training. He gives general direction to the collection, preparation and disposition of scrap materials. He operates the extensive facilities which each works has provided for the centralization and control of scrap accumulations, and also initiates studies leading to the reduction of scrap accruing from machining and assembly operations. r We recognize that each industry, each factory, each locality faces a different set of conditions. Irrespective of size, in a plant employing 20 men or a large works with 15,00 employees-certain fundamental principles can be applied for organization and the collection of scrap. The manufacture of electrical apparatus is a highly complicated intricate business. In one area we make electronic tube and adjacent buildings house the manufacture of steam turbines. Both are instruments for power-different in size and in many cases serving different fields of application. Unlike the great industries that build automobiles, trucks, or engines where repetitive operations, once solved, set the stage, electrical manufacturers' plants are large jobbing shops. Therefore, we must carefully plan every step in handling the scrap that "in so wasting to produce the finished motor or locomotive-we can return this product to reprocessing and ultimate use." ESSENTIALS FOR SALVAGE EFFICIENCY The essential steps for an efficient salvage department are: I. Selection of Personnel.-An able foreman would be a man of broad experience and understanding-one who commands the respect of the entire foreman group as well as that of his own men. In the interest of clean scrap, it is necessary to acquaint the foremen in the separate manufacturing units with the need of careful segregation in the shop and the best methods of obtaining it. At least 40 per cent contact time should be devoted for this purpose. 2. The physical equipment of the salvage department should be devised and selected to handle complete, at one point, all scrap operations. Bulky materials as produced by some large shops can be economically placed in salable shape at the divisional point. Because salvage

Jan 1, 1943

By R. J. Blackwell, J. R. Rayne, W. M. Terry

This paper presents results of an experimental investigation of factors that control the efficiency with which oil is displaced from porous media by a miscible fluid. The study was made to elucidate the relevant processes both on microscopic level (within individual or between neighboring pore spaces) and on macroscopic level (within a large sand body). Mixing of miscible fluids on the microscopic level was studied in sand-packed tubes. It was found that molecular diffusion is the dominant dispersion mechanism for reservoir conditions of rate, length and pore sizes. Macroscopic channeling was studied for various mobility ratios in reservoir mode1s scaled to relate viscous, gravitational, and diffusional forces. The formation of channels was due to viscous firzgering, gravity segregation and variations in permeability. With adverse mobility ratios, it was found for reservoirs of realistic widths that diffusion will not be effective in preventing the formation and growth of fingers, even in homogeneous sands. At sufficiently low rates channeling was eliminated by gravity segregation in tilted reservoirs. The dependence of recovery on mobility ratio, length-to-width ratio, flow rate and angle of dip is presented. INTRODUCTION Oil recovery by solvent flooding is finding increasing application in the field. While the process promises high recoveries from the region swept by solvent, under adverse conditions only a small fraction of the reservoir volume may be swept at the time solvent breaks through to the producing well. Further, the high cost of the solvent encourages its use only as a bank whose size must be kept at a minimum. Thus, two important questions arise: (1) what fraction of the reservoir can be swept, practically, by solvent? and (2) what is the minimum size solvent bank that can be used to carry out the displacement? The answers to these questions require knowledge of both macroscopic channeling processes and microscopic mixing processes. The studies described here were carried out to gain this knowledge. Microscopic mechanisms which cause mixing will be discussed first, because an understanding of these mechanisms is necessary for proper interpretation of the experimental work on channeling described later. INVESTIGATION OF MICROSCOPIC DISPLACEMENT PHENOMENA FOR MISCIBLE FLUIDS The questions to be resolved concerning microscopic displacement behavior are: (1) how much mixing occurs between oil and solvent in the direction of flow? and (2) does the solvent completely flush the oil from the pore spaces in the region invaded by solvent? The classic work of Sir Geoffrey Taylor' and its extension by Aris2 have shown that it is possible from theory to describe the amount of mixing in single straight capillaries when solvent displaces a fluid of equal viscosity and equal density. Aris showed that the lengths of the mixing zones, 81, corresponding to the 10 per cent and 90 per cent concentration levels of the solvent can be calculated from the formula, 81 = 3.62 \/Kt, .........(1) where 8l is the length of the mixing zone, cm; t is the time, sec; and K is the effective dispersion coefficient, cm2/sec, given by K/D = 1 + (1/48) (au/D)2......(2) where a is the radius of the capillary, cm; D is the molecular diffusivity, cm2/sec; u is the average velocity

By A. A. Rauth

In this paper, the author derives a series of formulas from basic principles and illustrates the application of these formulas to practical grinding charge problems. The paper establishes the nearly proportional relationship of mill power to charge level and shows how this relationship can be used to advantage in the design, manufacture and operation of grinding mills. The procedures outlined offer a rapid means of anticipating the power requirements for a given charge or determining makeup charges for a given power reading. The theories which have been developed concerning the dynamics of grinding have been motivated largely by the desire to derive workable formulas for the accurate determination of grinding mill power. This accurate determination has always been handicapped by the large number of variables to be considered, some of which are mill diameter, length, speed, liner design, grinding charge and material density, volumetric mill loading, internal friction, and charge size and gradation. In efforts to determine mill power, attempts have been made to utilize torque concepts. The formulas derived consider the grinding mass as a rigid body and utilize factors that are based on empirical data so that the movement of the charge within the mill is compensated for. These formulas are of principal interest to the mill designer since he is charged with the responsibility of analyzing all of the variables and from them predicting a mill size and a motor size which are compatible both with each other and with the grinding job for which they were intended. The mill operator, on the other hand, has an entirely different problem. The mill dimensions are already established and so is the driving motor. The operator's main concern is to obtain maximum performance from the mill he already has. For him, the many variables which perplex the mill designer are fixed values. Mill diameter, length and speed are set and, to a large extent, so are the conditions inside the mill from which the torque formula advocates derive their empirical data. It is conceivable that if the major mill dimensions which affect power demand are held constant, and if the variables influencing empirical factors are held fairly constant, as they are in an operating mill, then the only major variable remaining which affects power is the level of the charge itself. Under this condition.

Jan 1, 1970

By W. A. Backofen, M. L. Ebner

SINCE the early work of Gough with Hanson and Wright,l-3 the study of fatigue has been characterized by experiments on single crystals only in recent times.9-10 Now, increasing attention is given to this aspect of fatigue research for the insight that it may provide into details of mechanism. The investigations have concentrated to a large extent on the development of deformation markings on fatigued crystals, and have shown the cracks to originate in slip bands possibly preceded or accompanied by slip-band extrusions. Experiments of special interest to the present work were conducted by paterson5 on copper crystals and involved both metallographic examination and measurement of change in flow stress. Crystals were cycled in alternating tens ion-compress ion with a constant plastic shear-strain amplitude of approximately 0.8 pct, and were particularly revealing for their demonstration of hardening with accumulated strain similar to that in unidirectional straining, through an easy-glide stage I followed by a stage II of rapid hardening; deformation was not continued beyond 40 cycles, however, so that the eventual course of the hardening curve could not be decided. For the conditions used by Paterson, surface slip markings were similar to those observed in unidirectional straining, but there were no X-ray asterisms and no deformation bands on the surface. In current thinking about the nature of fracture in fatigue, two views relative to mechanism are generally acknowledged, with the reservation that both could apply simultaneously in some measure. As one possibility, fracturing across a slip plane is regarded as a result of loss of cohesion from the creation of many point defects by dislocation movement under the cyclic loading.'' On the other hand, fracture has also been taken to follow as a consequence only of the geometry of slip at a free surface, consisting of offsets and crevices which eventually become fatigue cracks.12, l3 The work of McCammon and Rosenbergl4 showing fatigue in polycrystals at 4.2oK makes clear that any long-range diffusion of point defects is unnecessary, yet studies such as those of Forsyth and stubington15 show that accelerated diffusion may be a characteristic of deformation under fatigue loading. Interest in obtaining data of possible use for resolving such questions led to the experiments described below. Copper crystals were conveniently loaded in alternating four-point bending at constant deflection, a test condition shown to approximate a constant plastic strain amplitude for a wide range of axial orientations. The method of testing was not readily adapted to extensive study of temperature effects, but an investigation of the geometry of crack formation could simply be made by orienting,the slip direction at different angles to the surface. EXPERIMENTAL PROCEDURES Single crystals were grown by a modified Bridgman method in a stationary gradient furnace under an atmosphere of purified dry nitrogen. Purity of as-grown crystals was 99.999 pct as determined by spectrographic analysis plus vacuum fusion and gravimetric analyses for oxygen and sulfur, respectively. Crystals were of reasonable perfection as evidenced by a critical resolved shear stress at 10 deg from (110) of 60 g per sq mm and half-width of a (400) X-ray diffraction line from one crystal of

Jan 1, 1960

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 V. D. Perry

Mineralized breccia pipes, because of their widespread occurrence and close structural relations to some of the world's great ore bodies, are objects of unusual interest for mining engineers and geologists. The literature contains many references to them, but it is questionable whether their genetic significance and economic importance have been sufficiently emphasized. The purpose here is to stress these features, relating them to the field facts, for the particular benefit of younger generations of geologists who, confronted with and sometimes confused by the growing flood of geochemical, geophysical, and other specialized research approaches, may be reassured that mappable field relations remain a foremost guide to a better understanding of ore deposits. A mineralized breccia pipe is a pre-mineral, breccia structure which has controlled the circulation and deposition of subsequently introduced mineralization. It is composed of relatively rotated angular or rounded rock fragments, set in a mineralized matrix. A pipe in plan outline may be circular, oval or approach polygonal form, with a steep to vertical axis proportionately much greater than its horizontal dimensions. The pipe is a steeply plunging, chimneylike mass of brecciated rock cemented with later minerals. Rock breaks in a variety of ways and complete fragmentation often occurs without rotation of individual pieces. A finely broken rock mass may fit into a tight jig-saw pattern, each fragment having mutually concordant boundaries with its neighbors. The result is a stockwork of innumerable reticulating cracks that, once cemented by mineralization, forms a complicated intersecting network of individually insignificant but collectively important seams and veinlets. Stockwork fracturing among its many forms takes the shape of domes of subsidence, fracture pipes, and related peripheral zones around and over breccia columns, or a combination of any of these structures. The significance of the mineralized breccia pipe is that it represents the extreme or climactic expression of a structural type which has a variety of mutations including subsidence domes, fracture pipes, and other stockwork zones all with related ancestry and similar to dissimilar characteristics. These allied and associated structures hold an answer to the fundamental question of the origin of many important ore deposits. CANANEA-TYPE LOCALITY FOR BRECCIA PIPES The Cananea district is characterized by an unusual development of mineralized breccia pipes. It is an important copper producer located in Sonora, Mexico, a short distance southwest of Bisbee, Arizona, and at the southerly limit of the great porphyry copper belt of the southwestern U.S. Cananea's rocks consist of Paleozoic quartzite and limestone capped unconformably by a thick series of volcanics including andesitic flows, tuffs, and agglomerates. These rocks have been intruded by a deep-seated granite with related basic and acid differentiates including dikes and plugs of quartz mon-zonite porphyry. Mineralization coincides with a northwesterly trending belt of these intrusives which break upward into and through the sedimentary-volcanic rock sequence. The district has weakly defined tectonic alignments in a northwesterly direction with subordinate intersecting fracture elements, but lacks important faulting or fissuring to provide throughgoing avenues for the upward circulation of mineralizing fluids. Thus, as will be discussed in subsequent paragraphs, the alternate way in which late magmatic and hydro-thermal derivatives of the parent magma reached the near-surface zone was by excavating their own breccia pipe channelways. There are numerous stages of breccia pipe development, related both in time and space to magmatic activity. A compilation of similarities and differences in various pipes suggests that proximity or remoteness of a demonstrable or inferred magmatic source provides an orderly genetic basis for describing the following representative types. Cananea Duluth Type: There are no intrusive rocks within or close to the Cananea Duluth pipe; therefore, the existence of any deep-seated magma that may have been related to its formation must be inferred. The structure is an oval-shaped ring 1200 x 300 ft in plan dimensions, cutting steeply across low angle, bedded tuffs, and other volcanics; it has been developed by drill holes to a depth of 2000 ft below the surface. The ore follows the periphery of the pipe and is composed of intensely brecciated rock which is cemented by minor galena, sphalerite, chalcopyrite, quartz, carbonates, and adularia. There is a definite vertical zoning of sulfides with less galena, continuing sphalerite, and increasing chalcopyrite at deeper levels. Within the interior of the ore ring, the brecciation becomes progressively weaker and coarser, the whole indicating relatively gentle slumping with broken, thin tuff beds pre-

Jan 1, 1961

By J. C. Shyne, T. D. Gulden

The steady-state creep behavior, in compression, of indium containing a dispersion of atomized glass particles was studzed over a range of temperature, stress, and composition. The observed behavior was not consistent with a simple model of dislocation climb over dispersed particles. From room temperature to about 100°C the temperature dependence could be described by an activation energy about one and a half times that of self-dijfusion in indium while at higher temperatures a much higher ternperature dependence was observed. THE value of dispersion-hardened and precipitation-hardened alloys has long been recognized. The mechanism of hardening, however, is not completely understood, especially at elevated temperatures where dispersed particles increase creep resistance. Following the suggestion of Schoeck1 that the rate of creep in dispersion-hardened materials is controlled by dislocation climb over second-phase particles, Ansel and weertman2 derived steady-state creep equations for dispersion-hardened materials. Their derivation was based on dislocation climb over non-deformable spherical particles in a metallic matrix. For low stresses they predicted a linear relation between applied stress and steady-state creep rate. For stresses large enough to cause dislocations to pinch off around the precipitate particles, the theoretical creep rate was shown to be proportional to the fourth power of stress. For very high stresses an exponential stress dependence was predicted. In all cases the theoretical creep rate is linearly proportional to the diffusion coefficient. This provided the primary temperature-dependent factor. Subsequently, Ansel and Lenel3 reported some experimental data on SAP, a composite of aluminum and aluminum oxide, that was in good qualitative agreement with the above theories. Their experimental creep rates, however, were four orders of magnitude less than the theory predicts. This was attributed to a low density of active dislocation sources in the material; somehow the dispersed particles inactivated the dislocation sources instead of merely interferring with dislocation motion. Meyers and sherby4 have recently reported on the creep behavior of SAP. Their results were inconsistent with a simple dispersion-strengthening model and led them to conclude that a continuous network of oxide controls the high-temperature mechanical characteristics of SAP. Because of the rather complex mechanical behavior of SAP and the lack of agreement about the morphology of its oxide phase, it appeared desirable to study the creep of a dispersion-strengthened material of controlled morphology. For the present work a dispersion of spherical glass particles imbedded in a matrix of indium was chosen. The microstructure of this material was similar to the idealized uniform dispersion of hard spheres in a deforming matrix assumed by Ansel and Weertman in their theoretical treatment. SPECIMEN PREPARATION For experimental convenience, it was decided to use a low melting-point metal in this investigation. It was necessary that the metal not be oxidized readily up to its melting temperature and that it wet the second-phase particles. It was desired that the material chosen for the dispersed phase be available as a fine powder with closely controlled size and shape. The material chosen to satisfy these requirements consisted of a matrix of indium metal containing a dispersion of atomized particles of a soft soda-lime-silica glass ranging from 5 to 30 µ in diam. The glass particles had a highly regular spherical shape. The composite was made by stirring the glass powder into liquid indium maintained at a temperature just above the melting point of indium, 156°C. The creep specimens were prepared by a hot pressing operation. Five to ten g of the mixture were placed in a steel die and pressed at 125°C for 5 min under a pressure of 30,000 psi. The resulting specimens were cylinders about 3/8 in. in diam and ranged from 0.4 to 0.7 in. high. By following this procedure a uniform dispersion was obtained in the creep specimens. Typical microstructures of specimens containing twenty and forty volume percent glass powder are shown in Figs. 1 and 2. Unfortunately, grain size could not be determined by the metallo-graphic techniques used. CREEP TESTING PROCEDURE Creep tests were performed in compression under constant stress. The specimens were placed in

Jan 1, 1963

By Don L. Kendall, Martin G. Buehler

A simple physical model has been developed to fit carrier-removal data in silicon irradiated near room temperature with reactor spectrum neutrons. Commonly observed donor and acceptor defect energy levels are assumed to be introduced linearly with neutron fluence. The donor levels (in ev) are Ev + 0.16, Ev + 0.27, and Ev + 0.31 and the acceptor levels are Ec - 0.55, Ec - 0.40, and Ec - 0.1 7, where Ev and Ec are the valence and conduction band energies, respectively. The introduction rates of each level are adjusted to fit literature initial carrier-removal rate data. When normalized with respect to the Ev +0.27 level, the relative values of introduction rates are 5.3, 1.0, 3.1, 1.0, 2.0, and 20.0, respectively for the six levels indicated above. To fit p-f (hole concentration vs neutron fluence) and n-f (electron concentration us neutron fluence) data, the introduction rates are multiplied by a factor which preserves the relative values given above. This factor depends upon irradiation temperature, reactor energy spectrum, neutron fluence calibration, and oxygen content of silicon. An extensive study of the effect of neutrons on carrier-removal in silicon irradiated with reactor spectrum neutrons (E > 10 kev) has been given by Stein and Gereth1 (SG) and Curtis, Bass, and Germano' (CBG). They measured initial carrier-removal rates for both p- and n-type silicon over an impurity range typical of silicon devices. In this work, we attempt to fit a simple theory to this data to establish a usable relationship between hole and electron concentration, p and n, respectively, and neutron fluence f. The p-f and n-f relations are needed to assist in the design of neutron tolerant silicon devices and are needed to clarify presently used empirical resistivity-fluence relationships.3 Neutron damage in silicon produces a variety of defects ranging from simple point defects to defect clusters. For the purpose of this treatment, we assume that simple point defects dominate carrier-removal effects. In contrast to this view, stein4 has proposed that defect clusters are responsible for a significant portion of carrier-removal effects. In the following section, it is shown that the carrier-removal effect in n-type silicon with an electron concentration less than 1015 cm-3 can be explained adequately by assuming that the divacancy is the dominant defect and that its introduction rate is independent of the electron concentration. For electron concentrations greater than 1015 cm-= an additional acceptor defect center is needed, and for simplicity the A-center (vacancy-oxygen pair) has been chosen. Although the E-center (vacancy-phosphorus pair) can account for some of the results, the A-center was chosen because the E-center requires a more involved treatment which the presently available data do not justify. In p-type silicon three radiation-induced donor levels are assumed, namely the divacancy and two other centers of unspecified nature located at Ev + 0.16 ev and Ev to 0.31 ev. The donor divacancy at Ev + 0.27 ev is assumed to be introduced at the same rate in p-type as in n-type. However, this rate is too low to fit p-type initial carrier-removal data. The dominant centers in p-type silicon are assumed to be the Ev + 0.16 ev and Ev + 0.31 ev levels where the latter is not the divacancy. The introduction rates are chosen to fit initial carrier-removal rate data. Assuming that the introduction rates are independent of Fermi level, the ratio between them is fixed for subsequent p-f and n-f calculations. Using the same ratios, the initial carrier-removal rate data1,2 as well as p-f and n-f data1,5 can be fit provided the absolute value of the introduction rates are adjusted to account for irradiation temperature, reactor energy spectrum, neutron fluence calculation, and the oxygen content of silicon. THEORETICAL ANALYSIS This analysis is basically the same as that used by Hi116 to analyze electron damage in silicon except we express the degree to which an impurity level is ionized not in terms of the Fermi level, but in terms of carrier concentration. Landis and pearson7 have used the latter approach to analyze y-damage in silicon. Neutron-induced defects responsible for carrier-removal at room temperature are assumed to be simple point defects with no interaction between defects so that they may be represented by discrete energy levels. It is also assumed that no constituent of a defect complex is used up and defects stabilize shortly after irradiation. Defects are assumed to be introduced linearly with fluence according to the product Rtf where Rt is the defect introduction rate and f the neutron fluence. Taking into account the ionization of defects according to Fermi statistics, and considering charge neutrality where minority carriers are neglected, the n-f relation is where no is the preirradiation electron concentration. The parameter Nt is the electron concentration at which the ionized defect concentration is one-half the total defect concentration (Rtf) or where Et is the defect energy level. For silicon at 300°K, ni = 1.45 X 1010 cm-3 and Ei = Ev+ 0.542 ev which was determined using Ec — Ev = 1.11 ev and me* = 1.07 mo and mh* = 0.558m0. The spin degeneracy factor, which usually appears as a number multiplying the Nt/n term of Eq. [1], is taken as unity. In effect, this factor has been incorporated into the defect en-

Jan 1, 1970

By J. E. 214-000-000-010 Warren, M. R. J. Wyllie, T. Meidav, L. Scharon, R. Uhley, A. J. deWitte

By considering the stoichiometry of the underground combustion process, an equrztion was derived relating the point velocity of the combustion front as a function of the air flux, fuel content, efficiency' of oxygen uitilization, hydrogen:carbon ratio, and ratio of carbon dioxide In carbon monoxide produced. In addition, an equation for cnlcrllating the amount of injected air necessary 10 produce the oil from I Cu. ft of reservoir space was derived. This equation has a direct bearing on evaluating compression costs for the process. Values calculated from these equations are in good agreement with experimentnlly obscrvcd values obtained from laborator) tube runs reported in the literature. The equations are also u\eful for atlcrlyzinf: experinlental data for inconsistencie.r.. I NTRODUCTION The injection of energy in the form of heat into an oil-bearing formation as a means of increasing the recovery of oil has been the subject of much discussion and experimentation as far back as 1920. The ideas have ranged from the injection of steam to the direct use of electricity as a means of getting thermal energy into a formation. However, most laboratory and field studies have been confined to the process whereby a portion of the crude in the reservoir is burned to generate the necessary heat to increase the mobility of the oil by reducing its viscosity. A recent paper by W. L. Martin, et al.', described the results of a number of laboratory experiments in which the combustion process was studied in a sandpack. Studies such as this indicate that with the temperatures involved, the oil is partially coked in this process and it is this coke which provides the fuel for combustion. The combustion of this residual coke can be analyzed theoretically by applying simple stoich-iometric principles to the combustion process. It is the purpose of this paper to compare the results reported in these experiments with those obtained from theoretical considerations. The variables studied by Martin, et al, among others, include combustion front velocity, air flux, oxygen utilization, fuel concentration, ratio of carbon dioxide and carbon monoxide in combustion gases, and the hydrogen:carbon ratio of the residual fuel. All these variables come into consideration when the stoichiometry of the combustion process is studied. The experimental conditions and results are tabulated in Tables 3 and 4 of Ref. 1. For purposes of comparison with calculated results, it was necessary to change the units of a number of the variables. UNDERGROUND COMBUSTION STOICHIOMETRY The combustion or oxidation of the carbonaceous residue or "coke" on the sand can be described in terms of the equation: where n is atomic ratio of hydrogen to carbon, m is ratio of moles of CO2 to CO produced, and CH. is taken as the unit molecular weight, realizing that actually the coke has a large range of high molecular weights. Assuming that all the fuel laid down in the sand is burned, the ratio of moles of oxygen consumed to moles of coke can be written as follows. Moles 0MolesMolescoke In terms of moles O2 per pound of coke, where 12 + n is the molecular weight in pounds per mole. If the 0, consumed is measured in standard cubic feet (60°F and 1 atm) and the coke in pounds per cubic foot of formation, 2, then:

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By Mark J. Klein

An internal-friction profile induced by nitrogen in Cr-35 at. pet Re was studied as a function of nitrogen concentration and heat treatment. From these studies, the solubility of nitrogen in this alloy was deduced. The equilibrium solubility of nitrogen in chromium does not appear to be substantially changed by the addition of 35 at. pet Re to chromium. However the precipitation kinetics of nitrogen from solution are significantly decreased. IT has been found that rhenium has a beneficial effect on the low-temperature ductility of Group VI-A metals. As part of a study to investigate the mechanism of this increased ductility, it was of interest to compare the concentration of interstitials retained in solution in one of these ductile bee alloys with the concentration retained in solution in its base metal. In this regard, a study of the solubility of nitrogen in ductile Cr-35 at. pet Re was selected for this phase of the investigation because the solubility of nitrogen in the base metal, chromium, has been determined over a fairly large temperature span.1-3 In addition, nitrogen is known to have an embrittling effect on chromium,4 but does not appear to have a comparable embrittling effect on Cr-35 at. pet Re.5 The experimental method selected for this study was the measurement of internal friction. This technique allows interstitial-concentration measurements to be made in bee metals that are sensitive to small changes in the concentration of a particular interstitial, even when other interstitials are present. When subjected to a cyclic stress, interstitial atoms in solid solution give rise to an internal-friction peak whose origin is the stress-induced ordering of interstitials into preferred lat- tice sites. These internal-friction peaks signify a dissipation of elastic energy which is greatest when the period of oscillation is comparable with the relaxation time for the attainment of an equilibrium interstitial distribution. For dilute interstitial binary alloys, the height of the internal-friction peak is proportional to the concentration of interstitial atoms in solid solutions. snoek8 has shown that the internal friction, Q-1, may be expressed by the equation where A is the relaxation strength, ? is the angular frequency, and 7 is the relaxation time. The relaxation strength, A, is proportional to the concentration of interstitials taking part in the relaxation if no interaction occurs between solute atoms. Thus values of peak height, which can be determined as a function of temperature, give a measure of the relative concentration of interstitial atoms in solution at various temperatures. If equilibrium has been attained, the peak heights will be proportional to the interstitial solubilities. The constant that relates A to the interstitial concentration can be determined experimentally from peak-height measurements as a function of concentration for specimens of constant crystallographic orientation. When this has been done, Q-1 measurements can be used to determine solubilities according to the equation C - KQ-1 where C is the soluble interstitial concentration and K is a constant that can be experimentally determined. It is possible, however, that lattice imperfections may bind interstitials to certain lattice sites. Internal friction does not measure these interstitial~, but does measure interstitials which are free to move in the normal lattice sites when subjected to an applied oscillating stress. On alloying the pure metal with a substitutional

Jan 1, 1965