Search Documents

By Verne E. Hooker, Robert M. Becker

The evaluation of in-situ rock stresses from measurements of the change in diameter of a borehole that is stress-relieved by overcoring has become a common practice. Generally the stress evaluation has been based on isotropic relations whereas rock in varying degrees is aniso tropic. A plane elasticity solution fora circular hole along one axis of an orthotropic medium has been adapted to evaluation of in situ stress. The necessary elastic parameters for this evaluation may conveniently be obtained from static tests on the overcore in conjunction with a consistent orthotropic thick-wall cylinder solution. The method and equations for the laboratory and field evaluations are summarized and illustrated. The determination of in-situ rock stresses from boreholc deformation measurements has gained widespread acceptance. Hastl measured the change in diameter of a borehole that was stress-relieved by drilling a larger concentric core; this overcoring method has produced a major portion of the field data for recent stress evaluations in rock. These evaluations are based on the elastic properties of the rock, which are conveniently determined by testing the overcore.6-8 The stress evaluation and the overcore testing procedures generally assume isotropic relations whereas rock is anisotropic in varying degrees. Recently, an anisotropic plane stress solution was adapted for evaluation of in-situ stress9 and a related plane strain solution was developed for the testing of stress-relief core. 10 The present paper unifies the stress evaluation and overcore testing procedures for a general class of anisotropic rock. Equations are presented for calculating the state of stress in a plane normal to the borehole, and a consistent thick-wall cylinder solution is given for determining the necessary elastic properties. In addition, a complete example of stress evaluation from laboratory and field data is presented, and the order of error in computing the state of stress from isotropic relations is illustrated. Several limitations of the results should be emphasized. All results assume small linearly elastic deformations of a homogeneous medium. Rock characterized by orthotropic, transversely isotropic, and isotropic elasticity is considered; however, the anisotropic results assume a special relation among certain elastic constants in addition to geometric orientation with respect to a plane of elastic symmetry. For some stress determinations in rock, these may be severe limitations. On the other hand, the anisotropic evaluation procedure involves little additional effort and is much more flexible than the corresponding isotropic evaluation procedure. The collection of field data is identical. The laboratory determination of elastic parameters is no more difficult and time-consuming for the anisotropic than for the isotropic model. Only the calculations require some additional effort for the anisotropic evaluation, and this additional effort is eliminated by use of a small computer. The flexibility of the anisotropic procedure is greater because more than two independent elastic constants are available in fitting rock data. For the special relationship assumed among orthotropic constants, five independent elastic constants are available for the determination of in-situ secondary principal stresses as compared with two independent elastic constants for the isotropic model. Much of the material presented is either directly from the literature or adapted from the literature through superposition or the usual transformations. Consequently, only enough detail for clarity of presentation is included. RELATIONS FROM ANISOTROPIC ELASTICITY Materials that have three mutually orthogonal planes of elastic symmetry are called orthotropic. By directing the axes of the coordinates perpendicular to these planes, Hooke's law for an orthotropic solid 11 is

Jan 1, 1970

By G. N. Goldberg, Robert A. Rapp

The steady-state kinetics and microstructures of simultaneous internal oxidation and external scale formation were investigated for the oxidation of Cb-Zr and Ch-Zr-Re alloys in pure oxygen at 1000°C. For each binary alloy, a constant rate of scale formation was observed at steady state; the internal oxidatton zone reached a steady-state, time-independent thickness. The dependence of the scaling rate on zirconium content exhibited a maximum at a low zirconium colttent; more concentrated alloys oxidized at a slower rate than pure columbium. This composition dependence of the scaling rate may be partially attributed to the effect of the internal oxidation process as an intenla1 sink for oxygen, with the formation of voluminous, and relatively impermeable ZrO2 precipitates. However, the morphologies of the ZrO2 internal oxide precipitates also affected the scaling rate. The internal oxide precipitates were distributed uniformly in the external scale. The steady-state thicknesses of the internal oxidation zones were not in agreement with those predicted theoretically from an idealized simple model. However, in this particular alloy systen1 the ideal model is not satisfied experitmentally. For ternary Cb- Zr - Re alloys with Nr10e = 0.02 or 0.05, the steady-state thicknesses of the internal oxidation zones were less than those for the corresponding binary Cb-Zr alloys. For ternary alloys with NRe(0) - 0.05, a more adherent and much less porous external scale was formed, and a reduction in the kinetics of scale formation was observed. A number of authors'-7 have reported linear oxidation kinetics for the reaction of pure columbium at 1000°C in air or in oxygen of about 1 atm pressure. A thin and tightly adherent oxide layer is found at the metal/oxide interface beneath a porous external layer of scale. Although the suboxides CbO and CbO2 represent thermodynamically stable phases at 100O°C, generally only Cb2O5* is identified as a reaction product in quenched specimens. Studies at 1000°C at lower oxygen pressures8,10 (Po2 - 10-4 Torr) show that both CbO and CbO2 do form and that CbO2 is a protective oxide. The growth of CbO2 results in parabolic kinetics until Cb2O5 is nucleated and grows. Apparently because of the large volume change associated with its formation, Cb2O5 is porous and does not serve as a diffusion barrier; linear kinetics prevail after the surface is covered with Cb2o5. On the basis of these observations, several authors5,9,10 have suggested that the linear oxidation of pure columbium at 1000°C in air or oxygen of 1 atm is limited by a Loriers-type mechanism,11 whereby the total rate or reaction is controlled by ionic diffusion through a thin layer of CbO2 which is maintained at a constant thickness throughout the linear oxidation. However, the concept of a Loriers model for a three-phase CbOlCbO2 Cb2O5 scale, with the only oxygen activity gradient across the NbO2 phase, is not self-consistent. Further, the suboxides CbO and CbO2 were not observed in the high-temperature X-ray diffraction study of columbium oxidation by Goldschmidt.8 From investigations in which the large PO2, dependence of the linear scaling rate was demonstrated, the importance of an oxygen adsorption or dissolution step has been suggested.4,5,7,10,12 Thus the mechanism for the linear oxidation of pure columbium remains quite controversial. Since further insight into the oxidation mechanism of pure columbium is not provided by this investigation, the authors wish to emphasize at the outset that their experimental results are not interpreted in terms of a particular rate controlling step for the oxidation of pure columbium. Previous investigations13-18 of the oxidation of Cb-Zr solid-solution alloys in air or oxygen at 1000°C suggest a remarkable dependence upon zirconium content. For an alloy of bulk zirconium mole fraction, N2r(0) equal to 0.05 or 0.10, the oxygen uptake is reported to be as high as four times greater than that for pure columbium;13,14 for alloys with NZr(0) 0.05 or 0.10 the oxygen uptake is reported to decrease with increasing NZr(0); to a minimum uptake of about one quarter that for pure columbium at NZr(0): = 0.50.13-15Since zirconium forms a more stable oxide (ZrO2) than the lowest columbium oxide (CbO), and since columbium exhibits a high solubility19,20 and diffusivity21,23 for oxygen, the internal oxidation of the zirconium component to ZrO2 is expected at a reaction front ahead of the advancing metal/scale interface. The external scale is then formed at the metal/scale interface by the inward migration of oxygen through the scale, probably through a series combination of molecular and ionic diffusion. Internal oxidation in conjunction with external scale formation has been investigated by Maak24 for the oxidation of dilute Cu-Be alloys in pure oxygen at 850°C. For cu-Be and many other binary alloys,24-27 very small, essentially uniaxial internal oxide particles are formed in the most dilute compositions; for somewhat more concentrated alloys. the internal oxide particles precipitate as platelets or needles. From pertinent solutions to the diffusion equation, Maak28

Jan 1, 1967

By F. H. Bunge, W. H. Jr. Dailey

This paper covers an investigation of the magnetic roasting characteristics of several iron ores when treated in the form of small pebbles (1/8" to 1/2" in size), and in the form of balled fines. The roasting treatment was carried out with apparatus capable of very close control of time, temperature, face velocity, and per cent combustible in the reaction gas. The effect of rapid initial heating was not studied although it is believed that this has a major effect on the results when processing solid pebbles. The test equipment is a fixed bed type, with such a wide range of gas velocities, temperature, combustible content and etc., that the results obtained can be made to approximate those from various commercial roasting processes. The data presented in this report covers reaction rates, effect of processing variables on silica in the concentrate, flow resistance of beds, and to a limited extent, the effect of particle size. The work is far from complete and raises many more problems than it gives answers for. We believe, however, that much of the information will be useful to those engaged in investigation of magnetic roasting processes. We measured the time required for various degrees of completion of the reaction of hematite to magnetite, for finely pound material in balled form and for carefully sized crushed ore. Measurements were made for a variety of Minnesota and Labrador ores. This report covers an ore requiring fine ginding (-325 M) for liberation (Lower Kevin) and an ore requiring coarser grinding (-100 M) for liberation (Upper Olson). These ores represent approximately the range of characteristics for a group in the west end of the Mesabi Range. Our data shows that the reaction rate is very fast for the finely ground material and that tests are readily reproducible. The reaction rate of the larger size material, in the range of 1/8" to 1/2", is much slower and the results more erratic. We found that for the ores tested, the percentage of silica in the Davis tube magnetic concentrate was a function of iron recovery. Normally the concentrate had a low and relatively uniform percentage of silica up to a critical value of recovery. Above this critical recovery value the percentage of silica increased very rapidly. For example, for one particular ore roasted in the form of 3/8" diameter balls made from -28 mesh fines, the silica remained at about 7% to 8% until the recovery exceeded about 93%. Increasing the recovery beyond 93% resulted in over 25% silica with the additional recovery. The degree of conversion of hematite to magnetite (% completion of roast) required to obtain 93% recovery for balled material, varied inversely with the gas flow rate. Percentage silica in the concentrate resulting from roasting 3/16" and 1/2" solid material was relatively erratic. However, for varying degrees of iron recovery the lowest silica percentages corresponded to those obtainable when roasting fines. It is thought the probable reason for the erratic values was the variation in shattering or cracking of the pebbles during initial heating and roasting. There is some indication that very fast initial heating will greatly increase the permeability of the individual particles of ore to gas flow so that the reaction rate for coarse sizes will approach that obtainable with balled fine material.

Jan 1, 1961

By T. C. Boberg, R. B. Lantz

This paper presents a method for calculating the producing rate of a well as a function of time following steam stimulation. The calculations have proved valuable in both selecting wells for stimulation and in determining optimum treatment sizes. The heat transfer model accounts for cooling of the oil sand by both vertical and radial conduction. Heat losses for any number of productive sands separated by unproductive rock are calculated for the injection. shut-in and production phases of the cycle. The oil rate increase caused by viscosity reduction due to heating is calculated by steady-state radial flow equations. The response of successive cycles of steam injection can also be calculated with this method. Excellent agreement is shown between calculated and actual field results. Also included are the results of several reservoir and process variable studies. The method is best suited for wells producing from a multiplicity of thin sands where the bulk of the stimulated production comes from the unheated reservoir. The flow equations used neglect gravity drainage and saturation changes within the heated region. INTRODUCTION This paper presents a calculation method which can be used to predict the field performance of the cyclic steam stimulation process. The calculation method enables the engineer to select reservoirs that have favorable characteristics for steam stimulation and permits him to determine how much steam must be injected to achieve favorable stimulation. While the calculation represents a considerable simplification of physical reality and the results are subject to numerous assumptions which must be made about the reservoir, it has been found that realistic calculations can be made of individual well performance following steam injection. The duration of the stimulation effect will depend primarily on the rate at which the heated oil sand cools which, in turn, is determined by the rate at which energy is removed from the formation with the produced fluids and conducted from the heated oil sand to unproductive rock. A complete mathematical solution to this problem is a formidable task, and finite difference techniques would undoubtedly have to be used. The calculation method pre- sented here utilizes analytic solutions of simple related heat transfer and fluid flow problems. The method is sufficiently simplified that it can be used as a hand calculation, although the calculations are somewhat lengthy and laborious. For that reason, the analysis was programmed for an IBM 7044 digital computer. Well responses observed at the Quiriquire field in eastern Venezuela&apos; have been matched using this program after making suitable approximations for reservoir and wellbore conditions. One of the most valuable uses of this calculation method is to assess the effect of reservoir and proc-cess variables on the stimulation response. This paper contains results of several studies made of key reservoir and process parameters. Among the most important of these is the influence of prior wellbore permeability damage. If a well is severely damaged prior to stimulation, a higher stimulation response will be observed than if it is undamaged. If a portion of this damage is removed, a permanent rate improvement will occur. THEORY I)ES(:KJJ&apos;TION OF CALCULATION METHOD The process of cyclic steam stimulation is essentially one of reducing oil viscosity around the wellbore by heating for a limited distance out into the formation through the injection of steam. Suitable modifications of the calculation technique presented here can be made so that stimulation of wells by hot gas injection or in situ combustion can also be calculated. A schematic drawing of the heat transfer and fluid flow considerations included in the calculation method is shown in Fig. 1. In brief, the calculation assumes that the oil sand is uniformly and radially invaded by injected steam. For wells producing from several sands, each sand is assumed to be invaded to the same distance radially. In calculating the radius heated rn energy losses from the wellbore and conduction to impermeable rock adjacent to the producing sands are taken into account. After steam injection is stopped, heat conduction continues and oil sands with r < ra cool as previously unheated shale and oil sand at r > r, begin to warm. The effect of warming of oil sand out beyond r, has little effect on the oil production rate compared to the effect of cooling of the oil sand nearer the wellbore than ra. Thus, in computing the oil production rate, an idealized step function temperature distribution in the reservoir is assumed where the original temperature exists for r > rn and where an average elevated temperature exists for r < rn. The average temperature in the oil sand for the region r < rn is computed as a function of

Jan 1, 1967

By H. D. Outmans

The paper presents a theoretical approach to the drilling problem based on rock mechanics and drilling fluid hydraulics at the bottom of the hole. The volume of the fractured rock around the vevtically penetrating bit tooth is expressed as a function of the tooth pressure, included tooth angle, drilling fluid pressure and formation characteristics on the assumption that fracturing is preceded by plastic flow. This assumption is verified by comparing the stress-strain state in the rock around the tooth to that in equivalent tri-axial compression tests described in the literature. An expression similar to the above is derived for the tooth penetration in the cuttings which have been retrained at the bottom. Combination of these two expressions leads to a relation between depth of tooth penetration in the virgin rock, the amount of cuttings retained and the other variables. By relating the cutting volume to the drilling rate, the hydraulic bit horsepower and the drilling fluid pressure, it is possible to arrive at formulas which express the drilling rate for viscous or turbulent flow at the cutting surface, for instant clearance, for retained cuttings, or for false tooth foundering as functions of such variables as the drilling fluid pressure, hydraulic horsepower at the bit, rotary speed and bit weight. Plotting the results as drilling curves and comparing them to published field and laboratory data justify the conclusion that the drilling model presented in this paper approximates the actual mechanism. The drilling curves serve to explain some characteristic features of the rotary drilling operation, and the equations may be used for numerical evaluation of the effect of changes in the magnitude of some of the variables on the drilling rate. INTRODUCTION Published laboratory data have established that consistent and relatively simple relations exist between the drilling rate and any one of the important variables. For instance, the drilling rate usually increases at an increasing rate with the bit weight. It increases at a decreasing rate with the rotary speed, decreases at a decreasing rate with increased drilling fluid pressure and plots as an S-shaped curve against the hydraulic horsepower at the bit. This relative simplicity seems difficult to understand considering that the drilling mechanism is a continuous process of generation and removal of innumerable individual cuttings, under different conditions for each. On the other hand, the drilling rate (although determined by this complicated process) only expresses the average rate at which cuttings are produced and cleared away. Ignorance of the drilling mechanism as far as individual cuttings are concerned, therefore, does not exclude the possibility of arriving at a concept of the drilling rate provided the mechanism is properly evaluated for the average cutting. It then may be concluded from the first sentence that this concept also should be rather simple in its final form, although not necessarily in its derivation. Two aspects of the drilling mechanism are difficult to evaluate — (1) the effect a penetrating drilling fluid has on the stresses around a bit tooth entering the formation and (2) the volume of rock fractured by a penetrating tooth which moves both vertically and horizontally. Thus, the scope of this paper has been limited to an evaluation of the drilling mechanism in formations which have a low permeability and are drilled with hard-formation bits. In the following sections, the drilling mechanism is analyzed in three steps. First, the effect of tooth pressure on the volume of fractured rock around the tooth is evaluated, assuming initially that the bottom of the hole is free of cuttings and, then, determining the effect of a layer of these cuttings. In the second section, the thickness of this layer is related to the drilling rate and the chip clearance time, and the latter to the hydraulic horsepower at the bit under conditions of laminar and viscous flow at the cutting surface. In the third section. the previous results are combined with the rotary speed into drilling equations which are then plotted and compared to experimental drilling curves. This section also contains some numerical examples of drilling rate computations. Some mathematical details are given in the appendixes. Finally, it should be mentioned that, although some drilling equations are already available, these equations (with the exception of one&apos; which attempts to relate drilling rate to drilling fluid pressure) are all mathematical expressions for observed laboratory data; none considers the effect of delayed drill-cutting removal, an essential feature of the drilling mechanism.

By Ulf S. Landergren

Diffusion coefficients and marker movements have been determined in brass using welded couples. Three different concentration ranges were employed at 750°C, while a fourth concentration range was measured at 500°, 600°, 700°, and 800°C. Marker displacements toward and porosity development in the high zinc side of the couple were observed in all cases. The results were interpreted as favoring a vacancy diffusion mechanism. IN the large amount of diffusion literature which has accumulated during the past half century, very few experimental data have been presented on body-centered-cubic metallic phases. The impressive combination of experimental and theoretical studies which has given the vacancy mechanism a fairly secure position as the preferred mechanism for most face-centered-cubic metals has no counterpart for body-centered-cubic phases, although the theoretical picture seems to be in a better state than the experimental work. The attempts to select diffusion mechanisms by the theoretical calculation of the activation energy for the various conceivable unit diffusion processes and the energy of formation of the type of lattice defect required have been based on models of two kinds. The first kind utilizes potential functions fitted for certain properties of pure crystalline copper; the second kind is fitted for sodium. The latter is an example of an open metal in which the metallic ion is small compared with interatomic distance,&apos; while the former is characteristic of the larger class of metals in which the ions have diameters about as large as the interatomic distance. brass should resemble copper in this fundamental property. It should be evident that results on the mechanism of diffusion obtained for one of these classes need not necessarily apply to the other. Huntington and Seitz, using a model of the copper type, showed that the activation energies to be expected for interstitial migration or a two-atom direct interchange were much higher than that for vacancy diffusion. Zener pointed out that a four-atom ring rotation would involve an activation energy considerably lower than the two unlikely mechanisms although still higher than the vacancy mechanism in the face-centered-cubic structure. Furthermore, he argued, the four-atom rotation would be relatively easier in a body-centered-cubic lattice and might be competitive with the vacancy mechanism. Paneth showed that a body-centered-cubic crystal-like sodium could form interstitial defects of special type, called a crowdion, by crowding an extra atom into < 111 > directions where the defects would affect a line of about eight atoms. The defect was constrained to move along the defining line. Unlike the face-centered-cubic case, therefore, the theoretical calculations for the body-centered-cubic mechanism do not provide a very clear-cut answer, and the absence of adequate experimental studies is the more serious. This investigation was undertaken for the purpose of obtaining useful information on diffusion in body-centered-cubic phases. brass was chosen for several reasons: 1—the phase field is sufficiently wide at high temperatures to permit the use of diffusion couples with adequate concentration ranges for satisfactory chemical analysis, 2-—thermodynamic activity data are available for the calculation of mobilities from diffusion coefficients, 3—the phase is as close to the theoretical model of the copper -type on which the calculations have been done as any body-centered-cubic metallic phase is likely to be. Measurements of the movements (if found) of inert markers during the diffusion process would indicate that a vacancy or Or mechanism played a major role in diffusion in brass. If no marker movement should be found, a ring mechanism would be indicated but not assured. Experimental Procedure The alloys used in these experiments, listed in Table I with analyses, were supplied by the Naval Research Laboratory. They were from ingots of 3x3 in. cross-section, heated to 760°C, reduced under a

Jan 1, 1957

By W. A. Riggs

In a single year the total transportation cost equals nearly 30 pct of the value of mineral commodities, the largest single cost from the deposit to consumer. The magnitude of this economic factor calls for more complete understanding of cost and operational problems between producer and carrier than now exists. Details are given of the costs of providing transportation as well as freight rates for selling transportation. In 1958, transportation of mineral commodities in the United States required 600 billion ton-miles to move 1 3/4 billion tons of products of mines comprising solid fuels, crude oil, industrial minerals and rocks, and ores. Possibly 17 pct was handled by more than one mode of transport, since combinations of several modes of transport may, at times, achieve better service/cost results than any single mode. Total transportation cost of nearly $5 billion equaled roughly 30 pct of the value of the mineral commodities, the largest single cost from the mineral deposit to the consumer. The tremendous magnitude of this transportation undertaking demands much more complete understanding than now exists between mineral producers and carriers about the economic and operational problems of each other. A substantial portion of the daily work of the writer is spent in trying to act as a combination interpreter and catalyst between potential shippers and carriers, hoping to direct both to points of mutual understanding and advantage. For the purpose of this paper, it is necessary to ask mining people to accept some terminology and to try to appreciate the mental approach of the freight rate making people of the carriers. This paper will consider costs of providing transportation as well as freight rates for selling transportation—although these are fields of limited certainty. There are no accurate statistical breakdowns for all modes of transport and types of carriers indicating commodities, volume of traffic, distances, costs, and rates. Much information is available for railroads, and some for the regulated portions of air, motor, pipeline, and water transportation. Unregulated carriers, handling about half of the transportation of products of mines, are not required to report, and very little is known of their activities. Figures used herein are national in scope, unless otherwise indicated, and have generally been derived from Interstate Commerce Commission reports and files. Tonnage figures have been reduced to net tons of 2000 lb. The very generalized comparison of average mileage freight rate levels for various modes of transport provided by Fig. 1 affords advantageous perspective for the entire field of transportation costs. Dashed guide lines show cents per ton-mile, and dotted guide lines show percentage of first class railroad freight rate. While costs of other than railroad transportation more nearly follow constant cents per ton-mile functions, costs of railroad transportation show characteristic taper, or decrease in cents per-ton mile with increasing length of haul. This fundamental difference in cost is due largely to the fact that railroads and pipelines provide their own fully taxed roadway facilities while highway, waterway, and air carriers use government-provided tax-free roadway at no cost other than fuel, license, and certain excise taxes. Thus railroad and pipeline transportation might be said to carry heavy threshold costs in addition to those transportation costs which are directly variable with ton miles of transportation produced. The original basis of railroad freight rate structures lies in the classification of all articles of commerce into railroad class ralings ranging from 7 to 40 pct of first class (100 pct) rate level, dependent on value, hazard, loading characteristics, service requirements, and whether carload or less than carload. Classification tables and mileage-based (but arbitrarily grouped), rate scales are published as tariffs subject to approval by regulatory bodies. Used together, these determine class rates in cents per 100 lb between points in the U. S. Where local commercial or operating conditions render it advantageous for carriers to depart from

Jan 1, 1961

By H. W. Mossman

Three aspects of magnetite smelting are discussed. The first is the working out of equilibrium conditions for eliminating sulfur. The second is the influence of magnetite solubility on the difficulty of tapping the reverb matte. The third is an approximation of the equilibrium conditions in the reverb gases which govern whether magnetite is mode or reduced in the reverb slag by these gases and by any iron sulfide in the slag. MAGNETITE has had a varied history in the Hurley Smelter since its start in 1939. Magnetite determinations on the smelter products are made regularly only on the monthly composite samples. Variations on the monthly averages are shown in Table I. Magnetite which drops from the slag and matte in the reverb has some slight bottom buildup which comes and goes, but no substantial accumulation from this source has been found at the end of a normal nine months&apos; furnace campaign. However, there has been some low grade magnetite bearing material mixed with considerable A1,0,, which has slid down from the bottom of the sloping flue between the reverberatory furnace and the waste heat boilers. This accretion has required drilling and blasting near the skimming end of the furnace. The magnetite has interfered with tapping at times. When the smelter was first started, tapping trouble from magnetite was extremely severe. Increasing the reverberatory furnace temperature by putting in an air preheater and a Dutch oven has helped greatly, although there still is occasional tapping trouble. When the present series of physical chemistry articles on copper smelting started coming out in 1950, they were read with interest, but no immediate application was seen for them. Results of some laboratory work in 1952 aroused a much stronger interest in this physical chemistry. A series of melts was made on some converter slags, which had magnetite in very large grain sizes, with the object of reducing the grain sizes in the slag, as it was known that it was easier to handle in the reverb in that condition. Anything done in the tests greatly reduced the grain sizes—even in the controls, where nothing was done except melt the slag and cast it. There was more magnetite in the slags after the tests than before, and with wide variations. There were no obvious reasons lor much of what happened in these tests. Much of the base material published in English in this field was made available for study. Recalculations were made on many of the type problems, and part of the data was reduced to local temperatures and compositions. Explanations were found for what happened in the 1952 series of tests on converter slags, and the same principles turned out to be a description of much of what magnetite does in the reverb. This article is to present the results of that study, from the viewpoint of applying the technical material in definite numerical form to the operating conditions in both the converters and the reverberatory furnaces at the Hurley smelter. Table I. Magnetite Variations on Monthly Averages, 1939 to 1955 Pet Magnetite Lowest Highest Average Converter slag 13.6 43.3 25.4 Roverb slaa 2.7 20.9 8.7 Matte 28 15.9 98 In general it was found that magnetite is made or reduced in both the converters and the rever-beratory furnace, depending on variations of temperature, matte composition, and reverb gas composition occurring in ordinary plant operation. Within reasonable limits, the field conditions for formation or reduction can be predicted, and probably can be set up and maintained. Converter conditions affecting magnetite formation can be put into numerical values better than for the reverb from purely technical calculations. The converter can be operated so as to keep the magnetite in the slag down to between 12 and 14 pct and still give satisfactory life for the converter brick. This depends upon having converter flux available which will make a slag with a good separation without raising the temperature too high. In the Hurley reverb and others with similar conditions, it is likely that a compromise of conditions will give a reasonably good control of combustion and still keep the magnetite from building up on the bottom. This discussion consists of three main parts. The first is the working out of the equilibrium conditions in the converter for determining in which direction the reaction 3 Fe,,O, (s) + FeS (1) F? 10 FeO (1) + SO, will go under actual converter operating conditions. The second deals with the influence of the solubility of magnetite in the slag and matte in the reverb on the difficulty of tapping matte. The third is an approximation of the equilibrium conditions in the

Jan 1, 1957

By R. I. Jaffee, C. M. Craighead, G. A. Lenning

The effects of increasing hydrogen content, introducing a notch, and changing the strain rate on properties of titanium and one of its alloys were investigated over a range of testing temperatures from —196° to 200°C. Both high purity and commercial purity A-55 titanium were used as representative a materials, while a commercial Ti-8 pct Mn alloy was used for an a-P alloy. It was found possible to analyze the data, using the ductile-to-brittle transition temperature concept. Increasing hydrogen, the presence of a notch, and increasing the testing speed raised the transition temperature for the a materials. The presence of hydrogen and notches raised the transition temperature of the a-B alloy also. However, increasing the testing speed generally decreased the transition temperature of the a-p alloy. IN previous papers, the authors have shown that hydrogen is most detrimental to the notch-bend impact toughness of a-titanium&apos; and to the tensile ductility of a-B titanium alloys.&apos; These general effects have been confirmed by Kessler and co-workers," and by Kotfila and Erbin." Kessler showed that the slower strain rate in notch-rupture tests was a more sensitive measure of hydrogen embrittle-ment of a-B alloys than the tensile test. In the present paper the authors cover, with tensile and notch-bend tests, the effect of testing speed, notches, and testing temperature on hydrogen em-brittlement. High purity and commercial purity titanium, examples of u materials, and a commercial a-,f3 Ti-8 pct Mn alloy were investigated. The purpose of the investigation was to correlate the factors producing embrittlement. Materials The high purity titanium for this investigation was obtained by arc-melting crystal bar from Foote Mineral Co. into 1 Ib ingots. As-cast hardnesses for these ingots ranged from 60 to 90 Brinell. The arc-melted ingots were forged to % in. diam rods at 1470°F. These rods were swaged to & in. diam at room temperature prior to machining of the specimens for hydrogenation and testing. The commercial A-55 titanium and the Ti-8 pct Mn alloy were obtained from Rem-Cru Titanium Inc., as % in. diam rod. The A-55 titanium was hot swaged to 1/4 in. diam at 1300°F prior to machining of the specimens for hydrogenation and testing. The a-p C-11OM alloy was hydrogenated prior to hot swaging to 1/4 in. roc1 at 1300°F. Hydrogenation prior to fabrication permitted testing of the a-P alloy in the equilibrated and stabilized condition. The equilibration and stabilization treatment used was an air anneal for 1 hr at 1300°F, followed by furnace cooling to 1100°F and then air cooling. Hydrogenation was done in a Sieverts apparatus described previously.&apos; The hydrogenation treatment consisted of heating for 1/4 to 1/2 hr at 820 °C in purified hydrogen, slow cooling to 680"C, holding 1 hr for homogenization, and rapid cooling to room temperature. Intended hydrogen levels for the a-tita-nium materials were 0.05, 0.5, 1.0, 1.5, and 2.0 atomic pct—about 15 to 430 ppm. For the Ti-8 pct Mn alloy, the intended hydrogen levels were 0.05, 0.5, 1.0, 2.5, and 5.0 atomic pct—about 15 to 1100 ppm. The material for the 0.5 atomic pct level of hydrogen was vacuum annealed prior to hydrogenation to increase the accuracy of the addition. Vacuum annealing to remove hydrogen consisted of heating for 6 hr at 820°C, and resulted in a hydrogen content of about 10 to 15 ppm. When hydrogen contents were above 0.5 atomic pct—100 ppm—the hydrogen was added without prior vacuum annealing.

Jan 1, 1957

By H. J. Morel-Seytoux

Methods of predicting the influence of pattern geometry and mobility ratio on water flooding recovery predictions are discussed. Two methods of calculation are used separately or concurrently. The analytical method yields exact solutions in a convenient form for a unit mobility ratio piston-like displacement. A few typical pressure distributions, sweep efficiencies and oil recoveries are presented for various patterns. For non-unit mobility ratio, one may resort to a numerical method, such as that of Sheldon and Douherty. 1,2 Because the domains of applicability of the analytical and numerical techniques overlap, the exact solutions provide estimates of the errors in the numerical procedures. The advantages of the analytical and numerical methods can be combined. To develop a numerical technique as independent of geometry as possible, the physical space is transformed into a standard rectangle. The entire effect of geometry is rendered through one term, the "scale-factor", derived from mapping relations. The scale factor can be calcu-lated from the exact unit-mobility ratio solution for the particular pattern of interest. By this means recovery performances for arbitrary mobility ratio con be obtained for many patterns. A sample of tesults obtained in this manner is presented. INTRODUCTION Pattern geometry and mobility ratio are two major factors in making a waterflood recovery prediction. Because assisted recovery has become increasingly important to the oil industry, pattern configuration and mobility ratio also assume a greater significance in the assessment of the economic value of recovery projects. The influence of pattern geometry and mobility ratio in shaping a recovery curve and on the other quantities of interest to the reservoir engineer is the main subject of this paper. Much effort has already been spent on estimating quantitatively the influence of either pattern or mobility ratio or both on oil recovery. The literature reports many investigations of this nature. 3-9 However, many results or methods of recovery prediction presented in the literature cannot be considered fully satisfactory. Even for unit mobility ratio and piston-like displacement, where analytical solutions are available, the literature shows discrepancies. For non-unit mobility ratio, the divergence in the results is extreme. For infinite mobility ratio in a repeated five-spot, depending on the investigator, the sweep efficiency ranges from 0 per cent to 60 per cent. With respect to the influence of pattern on recovery, only the repeated five-spot has received much attention. Other confined patterns and pilot configurations have received very little attention. Two calculation methods are presented in this paper, either separately or concurrently: the analytical method of potential theory and the numerical method of finite-difference approximation. The analytical method is more restricted in scope than the finite-difference method, but it has the definite advantage of providing exact solutions within its range of applicability. If a unit-mobility ratio piston-like displacement is assumed, the analytical approach is possible. A few typical results are reported in this paper; the detailed description of the general method and of a great variety of results will be the subject of other articles. Fornon-unity mobility ratio, we must resort to a numerical scheme. The numerical technique is that which was described by Sheldon and Dougherty.l,2 It is not limited to piston-like displacement. However, mainly single interface results will be presented here. Because the respective domains of applicability of the analytical and the numerical method overlap, useful comparisons of exact and numerical solutions can be made for a variety of patterns. The advantages of the analytical and numerical approaches can be combined. The reason for the success of this analytical-numerical approach can be summarized in the following two points: 1. The numerical solution for arbitrary mobility ratio can be programmed most efficiently when the physical space in which the displacement actually takes place is transformed into a standard shape; and 2. This can be done with remarkable simplicity whenever an analytical expression for the pressure

Jan 1, 1966

By Waynw P. Myers

Electricity, as normally thought of by a layman&apos;s definition, is a man- made force that has no color, no odor, is not visible, cannot be heard, yet man can control it and make it perform his work for him. A person dare not cheat on it, for it then will master him. It has laws, rules, and regulations that must be strictly adhered to if it is to perform in the safest, most reliable, and economical manner. In the early days of coal mining the power source was man himself. Later man employed steam and compressed air to reduce his muscle power requirements needed in the extraction of coal. By the turn of the 20th century, electricity had become the most common source of power for the coal industry. Steam and compressed air are still found in use in isolated cases today for specific applications; however, they are rapidly being replaced by electricity. Electrical laws and regulations are quite extensive and somewhat complex for the proper generation and use of electricity. Man-made laws, rules, and regulations for the care and handling of electricity, designed for man&apos;s safety, are rigidly enforced by Mine Safety and Health Ad- ministration (MSHA), as well as the various state departments of mines. BACKGROUND Until recent years most underground coal mines used direct current exclusively for all movable and for some stationary equipment. Exceptions were for permanent pumps, tipple, fan, and outside buildings where alternating current was normally employed. However, over 60 years ago there were mines with dual systems. Alternating current was used for cutting machines, drills, and pumps ; direct current for locomotives. There were no other machines as these were hand-loading mines. Later, during early mechanization in the 1920s, alternating current was also applied to loading machines at some mines. Use of dual electric systems then declined until fairly recently when

Jan 1, 1981

By H. B. Pulsifer

One of the most common basic open-hearth furnace products is a simple carbon steel with a carbon range from 0.05 to 0.15 per cent. The material is widely used for sheets, tubes, bars, wire and the innumerable special objects of secondary fabrication. The properties of the material vary over wide ranges, depending on the more exact composition and details of manufacture or treatment that affect its structure. The present paper is a review of those properties and structures that have been found by studying this steel as used in the cap screw industry. At first glance this might be thought a narrow application of the steel, but, since the industry uses the material in the cold-worked, annealed and heat-treated conditions, the title of the paper may be excused for its comprehensiveness. Yet, in studying the wide ranges of properties and structures that preclude a more specific title, there can be no suggestion of completeness or finality. Even a summary of a few years&apos; work is no more than a preliminary survey of a field that will be fertile in scientific data and industrial products for many years to come. An outline of the structures will be given first, then a list of the previous contributions that are at hand and finally some details of the work and related problems, especially concerning this low-carbon steel in the surface-hardened condition. Wire Used in Manufacture of Cap Screws The wire that is used for making cold-upset and roll-threaded cap screws has been cold-drawn from hot-rolled rods that were pickled and limed. The common sizes used in the industry range from 0.185 to 0.750 in. dia., the exact diameters depending on the plant practice and having the usual tolerances of 0.002 in. for the smaller sizes and 0.003 in. for the larger sizes. The carbon range, by pit analysis, is usually between 0.07 and 0.12 per cent. Manganese should be within the usual 0.30 to 0.50 per cent. range. Sulfur rarely exceeds 0.045 per cent. and phosphorus seldom exceeds 0.030 per cent. Silicon is usually below 0.03 per cent. but in the silicon-killed steels may be in the customary 0.10 to 0.25 per cent. range. Rimmed steel is not considered most suitable for cap screws that are cold-headed or for parts that are to be surface-hardened.

Jan 1, 1931

By W. A. Young

Polynomial functions of temperature were obtained for the specific heats, thermal diffusivities, and thermal conductivities of zirconium hydrides containing 4 at. pct U. Three hydrides (H/Zr atom ratios of 1.58, 1.65, and 1.70) were studied over the range a" to 900°C and a fourth (H/Zr = 1.81) was studied over the range 0° to 760°C. The specific heats were determined from enthalpy measurements which were obtained using a unique drop calorimeter specifically designed for use with materials in which high temperature phase transitions and/or high dissociation pressures occur. Thermal diffusivities were measured by the flash method using a pulsed laser. The thermal conductiuities were obtained as the product of specific heat, thermal diffusivity, and density. The specific heats agree, within 10 pct, with values derived using a theoretical model in which the hydrogen and zirconium atoms are treated as Einstein and Debye oscillators, respectively. RELIABLE values of the thermophysical properties of the fuel are required to predict the operating temperatures and temperature response of SNAP nuclear reactors. Among the most important of these properties are the thermal conductivity, specific heat, and thermal diffusivity. A considerable number of investigations1-4 have been made of these properties for the Zr-H and Zr-H-U systems.* However, little of the drides, however, this direct method cannot yield meaningful results, since the hydrogen will redistribute under the influence of the thermal gradient, thus forming a concentration gradient; hence, one has a spectrum of compositions, rather than a homogenous alloy. Although the "average" composition of the material may be identical to the initial uniform concentration, the directly measured value of conductivity will be dependent on the thickness of the specimen, due to the highly sensitive dependence of transport properties on hydrogen content. This dependence is strikingly illustrated by the work of Bickel,5 who found that the electrical conduction of zirconium hydrides ranges from primarily hole conduction to primarily electronic conduction, depending upon the hydrogen content. Fortunately, the direct measurement of thermal conductivity is unnecessary, since it can be expressed as the product of the specific heat, thermal diffusivity, and density, all of which can be directly measured with considerable accuracy. EXPERIMENTAL Specimen Preparation. The combined fuel-moderator material used in SNAP reactors is a hydrided zirconium-uranium alloy containing -10 wt pct U. The alloy used in this work was representative of that used in nuclear reactors except that normal uranium was substituted for the enriched uranium required for reactor usage. It was produced by a triple-arc-melt and double-extrusion process. All specimens were prepared from a single cylindrical extrusion which contained 10.30 pet U, 89,35 pct Zr, and 0.35 pct impurities, The specimens for each composition were hydrided simultaneously with ultrapure hydrogen (10 ppm total impurities) using standard fuel production techniques which routinely yield homogeneous, crack-free fuel with negligible increases in the impurity levels. The hydrogen content of each specimen was determined from its weight gain and the density was measured by liquid displacement, Chemical analyses yielded hydrogen concentrations which agreed with the weight gain data within ±0.02 in H/Zr atom ratio) the concentrations of all other elements agreed almost exactly with the initial values after adjustment for the added hydrogen. The specimens used for the determination of specific heat were centerless ground to 2.00 cm diam after hydriding. A thin slice was carefully removed From each end for metallographic examination. In every case, this examination revealed a uniform structure as evidenced by the appearance and distribution of the two phases present in the fuel at the hydrogen concentrations used. TWO specimens (H/Zr = 1.600 and 1.632) appeared to be entirely 6 phase with equi-axed grains; the specimen with H/Zr = 1.756 showed

Jan 1, 1970

By A. H. Geisler, J. H. Keeler

Preferred orientations in rolled and annealed titanium sheets were determined by the Geiger counter spectrometer X-ray diffraction technique. Five annealing textures dependent upon the temperature range of annealing were found, and in order of increasing annealing temperature pendent upon the temperature range are: 1—a deformation like texture, 2—a rotated inorder a-recrystallization temperature texture, 3-a retained u-recrysraIlization texture, on annealing at lower temperatures of the ß-region, 4—a transformation texture based on recrystallized a and predicted by the Burgers&apos; relationship, and 5—a ,ß-cube texture. These results are examined in terms of current theories of recrystallization textures. UMEROUS investigators have described the tex- ture obtained by cold rolling the hexagonal metals, titanium, zirconium, and beryllium, which have c/a ratios less than that of ideal packing, 1.633. The basal planes are rotated out of the rolling plane, about the rolling direction, so that the basal poles are tilted toward the transverse direction as shown schematically in Fig. la. In all instances but one,&apos; it was also reported that the [1010] direction was parallel to the rolling direction (see Fig. lb). Hot rolling has been reported as causing a similar tilt of the basal poles in the transverse direction (see Fig. la) and causing the [1010] direction also to be parallel to the rolling direction as shown schematically in Fig. lb. Annealing after deformation does not appreciably change the tilt of the basal poles in the transverse direction." Beryllium2-7 continues to have the [1010] direction in the rolling direction after annealing, and similar observations for titanium and zirconium&apos; . have been reported for annealing at fairly low temperatures, again as in Fig. lb. At higher annealing temperatures, however, the recrystallized grains of titanium" and zirconium have an orientation such that the [1120] direction is approximately in the rolling direction, although the basal poles are still inclined in the transverse direction. Figs. la and lc show the resulting orientations schematically. This change in orientation has been described as a nominally ±30° rotation of the hexagonal crystallites about the basal poles of the cold rolled texture and is apparent from the results which are summarized in Table I for investigations with the X-ray diffraction technique employing film. The angles y, , and ß are indicated in Fig. 2 which represents the stereographic projection of (1070) poles for the mean orientation of a pole figure. Texture determinations for titanium using the Geiger counter spectrometer have provided similar results except that in some instances additional components of the texture were proposed, as shown by the summary of data in the upper half of Table 11. On the other hand, the spectrometer technique, when applied to zirconium,* has revealed a splitting Recently completed studies of the textures of annealed zirconium", show zirconium to possess textures very similar to those reported here for titanium. Therefore, much of this discussion will include zirconium by virtue of its close similarity to titanium in pref erred orientations. of the intense areas of the pole figure for samples annealed at 600°C. This splitting could be described by a 7" rotation of the tilt axis about the normal to the rolling plane. Such a splitting for the annealed texture relative to the cold rolled texture was not observed in other determinations for either zirconium or titanium using the less sensitive film X-ray methoe and makes the relationship between the two types of texture more complex than the simple rotation about the (0001) pole based on film work. The more precise investigations on zirconium permit the descriptions in the lower part of Table 11, which show that the texture depends quantitatively on the temperature of annealing. When zirconium is annealed at temperatures up to 400°C, the texture is similar to the cold rolled texture, while annealing in the range 500" to 900°C produces a texture which is only approximately described as [11%] in the rolling direction. More precisely described results for zirconium show that the two types of splitting ( 1—about an axis in the rolling plane through an angle given in the second column in Table II and 2—about the normal to the rolling plane through an angle given in the third column of Table 11) depend on annealing temperature. The [1120 is the rolling direction only when the annealing temperature is in the vicinity of 900°C

Jan 1, 1957

By V.K. Zworykin

Throughout its development the science of electronics, like so many other branches of science and industry, has been indebted to the metallurgist. Metallurgy has provided the electronic engineer with conductors to furnish passage for electrons outside of his evacuated tubes: it has developed tungsten and tantalum wires for emitting electrons and various refractory metals and alloys suitable for electrodes, which receive them and guide them within the tube; finally it has supplied ferromagnetic materials both for shielding electronic apparatus from magnetic fields and for forming highly concentrated fields for deflecting electron . Today I hope to show you that electronics has not been altogether ungrateful for this assistance. More specifically, I wish to indicate the role that the electron microscope, a relatively recent product of electronic engineering, is playing even now in the study of metals. Every metallographer is aware that his standard metallographic microscope will not reveal structural detail involving separations smaller than a certain fraction —at the very best a tenth—of a micron. This limitation, he knows, is a consequence of the fact that the wave length of the light employed for forming the image has this order of magnitude. It is, thus, a property of all light microscopes. TO overcome it, a different imaging medium must be employed, The existence of such a medium became evident when it was demonstrated1 that electrons were affected by magnetic and electric fields with axial symmetry just as light is affected by glass lenses; these fields, or "electron lenses," curve the paths of electrons reaching them from an object located on the axis of symmetry in such a manner that, at some distance from the lens, a true image of the object is formed. Furthermore, although a wave length may . be ascribed to electron beams just as to light, this wave length is much shorter; for electrons accelerated through a difference of potential of 50 kilovolts, the wave length is only about one hundred-thousandth of that of visible light. Actually, the elcctron microscopes of today resolve detail one hundred times as fine as light microscopes; if the maximum useful magnification of the latter is taken to be a thousand, that of electron microscopes is one hundred thousand. The Electron Microscope Since electrons are readily scattered and absorbed by matter, specimens to be observed in the electron microscope must be prepared to be less than one micron in thickness; furthermore, the interior of the microscope must be carefully evacuated to prevent collisions between the clectrons and the molecules of the air. Finally, as already mentioned, the lenses are formed by axially symmetric magnetic or electric fields, and not by any material medium such as glass. Apart from these basic differences, the electron microscope closely

Jan 1, 1943

By V. K. Zworykin

Throughout its development the science of electronics, like so many other branches of science and industry, has been indebted to the metallurgist. Metallurgy has provided the electronic engineer with conductors to furnish passage for electrons outside of his evacuated tubes: it has developed tungsten and tantalum wires for emitting electrons and various refractory metals and alloys suitable for electrodes, which receive them and guide them within the tube; finally it has supplied ferromagnetic materials both for shielding electronic apparatus from magnetic fields and for forming highly concentrated fields for deflecting electron . Today I hope to show you that electronics has not been altogether ungrateful for this assistance. More specifically, I wish to indicate the role that the electron microscope, a relatively recent product of electronic engineering, is playing even now in the study of metals. Every metallographer is aware that his standard metallographic microscope will not reveal structural detail involving separations smaller than a certain fraction —at the very best a tenth—of a micron. This limitation, he knows, is a consequence of the fact that the wave length of the light employed for forming the image has this order of magnitude. It is, thus, a property of all light microscopes. TO overcome it, a different imaging medium must be employed, The existence of such a medium became evident when it was demonstrated1 that electrons were affected by magnetic and electric fields with axial symmetry just as light is affected by glass lenses; these fields, or "electron lenses," curve the paths of electrons reaching them from an object located on the axis of symmetry in such a manner that, at some distance from the lens, a true image of the object is formed. Furthermore, although a wave length may . be ascribed to electron beams just as to light, this wave length is much shorter; for electrons accelerated through a difference of potential of 50 kilovolts, the wave length is only about one hundred-thousandth of that of visible light. Actually, the elcctron microscopes of today resolve detail one hundred times as fine as light microscopes; if the maximum useful magnification of the latter is taken to be a thousand, that of electron microscopes is one hundred thousand. The Electron Microscope Since electrons are readily scattered and absorbed by matter, specimens to be observed in the electron microscope must be prepared to be less than one micron in thickness; furthermore, the interior of the microscope must be carefully evacuated to prevent collisions between the clectrons and the molecules of the air. Finally, as already mentioned, the lenses are formed by axially symmetric magnetic or electric fields, and not by any material medium such as glass. Apart from these basic differences, the electron microscope closely

Jan 1, 1943

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&apos;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 B. K. Larkin, H. R. Bailey

A model Of heat flow in an underground combustion process is studied. This model includes convection effects and thus is more general than previous studies which considered conduction as the only mechanistn for heat transfer. Both linear (tube run) and radial (field application) geometries are considered. The effects of ignition heaters, vertical losses and finite source width are considered for the linear case. The results are in the form of equations and are presented in graphical form for a number of cases. Convection effects increase frontal temperatures about 25 per cent over those computed for conductive transfer for typical field operating conditions. This increase in temperature is a result of heat being transported toward the front by the injection gas. Even greater temperature increases are realized as the per cent oxygen in the injection gas decreases. It is well known that compression costs are of considerable importance in estimating the economic feasibility of underground combustion. By assuming an ignition temperature for the combustion fuel, predictions of limiting conditions on fuel density and injection rate necessary to sustain the combustion zone are made. For typical field conditions, at least 0.75 lb/cu ft of fuel are needed with air as the injection gas. If the injection gas is 10 per cent oxygen and 90 per cent nitrogen, this figure is 0.69 lb/cu ft. INTRODUCTION The possibility of increasing oil recovery by underground combustion has been considered for many years. Recent field tests indicate that the underground combustion process is technically feasible. The economic feasibility depends to a large extent on the amount of air which must be injected to sustain combustion. Prediction of the success of employing underground combustion in a particular reservoir must be based on existing field tests, laboratory tube runs"&apos; and solutions of mathematical or analog models. Vogel and Krueger5 devised an electrical analog of the heat transfer problem in an underground combus-ion process. They considered the problem of heat conduction from a cylindrical source with increasing radius assuming no vertical losses. Bailey and Larkin and Ramey solved the corresponding problem including vertical losses by considering a mathematical model of the conduction process. Ref. 9, 10 and 11 also present information related to the subject study. The present paper generalizes these results to include both conduction and convection mechanisms for heat flow. A model of the conduction-convection process is described. The partial differential equations governing this model are written and solved for a number of cases of interest. The formulas are evaluated and the results are presented graphically. A CONDUCTION-CONVECTION MODEL Description of the Model A number of models are employed in this paper, each of which approximate a particular physical situation depending on the geometry, i.e., linear or radial flow, the type of vertical heat losses, the nature of the source, etc. The geometry for both the linear and radial case is shown in Fig. 1. A particular linear model is described in this section, and variations from this model are described in later sections of the paper. Consider a homogeneous, semi-infinite porous prism with lateral sides insulated as shown in Fig. 1. The left end, x = 0, is heated for a time t, with a plane heater supplying a constant heat flux and during this time; gas (nitrogen) is passed through the pores of the prism at a constant mass flux. At time t, the heater is turned off and air is injected in place of the nitrogen. This causes combustion to take place, and a uniform source of width 1 is formed which moves at a constant velocity v, to the right. The frontal velocity v, is related to the gas velocity v,, through stoichiometric considerations

By J. A. Herbst, D. W. Fuerstenau

This paper examines the zero order production phenomenon in the context of the size discretized batch grinding model. A restrictive interrelationship between the selection and breakage parameters of the model, which is mathematically sufficient to ensure zero order behavior, is delineated. Based on this interrelationship, a scheme is developed for predicting the values of the selection and breakage parameters for all size fractions from a minimum of experimental data. To test this scheme, dolomite was ground in a laboratory batch ball mill. Successful simulation of the comminution behavior of this system was achieved using the batch grinding model with the parameter values obtained from the proposed scheme. Historically, a preponderance of comminution research has centered on attempts to relate breakage energy to the performance of comminution machines. Concomitantly, grinding data were interpreted almost exclusively in terms of inherently empirical energy-size reduction relationships1-3 or "laws of comminution"4- 6 which were based on highly oversimplified descriptions of the fracture process. In some instances, these relationships provide a crude basis for the correlation of experimental data, but, invariably, this approach is inadequate for meaningful process simulation. The control and optimum design of comminution circuits require a mathematical model capable of depicting the size reduction behavior of every size fraction for grinding conditions of technological importance. Energy-size relationships do not provide this detailed information. For example, an energy-size analysis for batch ball-milled dolomite has recently been completed in the authors&apos; laboratories.7 This analysis provided a satisfactory correlation between the hypothetical size modulus, which characterizes only the finer portion of the size distribution, and the input energy. However, this relation did not constitute an adequate description of the complete system. In general, one of two distinct viewpoints toward simulation has been adopted in recent comminution research. The first approach focuses on the fracture of single mineral specimens, with the essential aim of representing the over-all process in terms of the breakage characteristics of individual particles and the characteristics of the stress field which the particles experience within the particular size reduction device. As illustrated by Harris&apos;s review,&apos; the actual incorporation of single-specimen fracture information into a description of the behavior of a multiparticle comminution system has seldom been attempted. The only model to accomplish this incorporation was the one developed by Schönert.8 The Schonert model treats single-passage grinding machines in terms of a distribution of effective loads acting on a particle, a distribution of required particle breakage energies, and the breakage product distribution. However, the use of single-particle fracture behavior to derive models of the size reduction process is at present limited to machines in which the particle residence time is approximately equal to the time required to apply stress. As a result, Schönert&apos;s analysis cannot be used to advantage to simulate the complex environment which prevails within a tumbling mill. For multiple-passage systems a second approach to the description of comminution processes has been more fruitful. This approach entails the formulation of a mathematical model which is phenomenological in nature in that it lumps together the entire spectrum of stress-application events which prevail in a system under a given set of operating conditions. The appropriately defined average of these individual events is then considered to characterize the over-all breakage properties of the device. Thus, to analyze the performance of a tumbling mill, the manner in which the particles of a particular size (or size fraction) are stressed need not be distinguished. Instead, a single parameter is assumed to represent the resistance of that size to fracture, given the average grinding environment which exists in the mill. The isolation of such a parameter and a related set of quantities, which constitute the breakage product size distribution for the average event in this size fraction, allows the formulation of physically meaningful descriptive equations capable of yielding precise and detailed information for simulation.

Jan 1, 1969

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

The rate of FeO and MnO reduction from silicate and aluminate slags by carbon-saturated iron is dependent on both slag composition and temperature. Owing to variable stirring rules during- the course of reaction, the reduction processes occur in two stages in stationary graphite crucibles. In the first stage the slirring caused by CO evolution results in forced convection conditions in the slag. As the boiling action subsides during reduction, the flow conditions within the system become those defined by natural convection. Analysis of the data by an unsteady-state penetration model indicates that the rate of reduction is controlled by the rate of cocur-rent flow of&apos; cations and anions from the bulk slag to the interface during both stages of reduction. The hearth reactions of the blast furnace have been extensively investigated from the viewpoint of equilibrium by deductions from controlled laboratory experiments. However, the literature contains relatively few studies that have been directed toward the area of slag-metal reaction kinetics, and the majority of these have been concerned with the rate and mechanism of sulfur transfer between iron and slag. More recently, attention has turned to the kinetic factors in the reduction of oxides from liquid slags. One of the first investigations in this field was conducted by Dancy1 on the reduction of pure liquid FeO and pure liquid Fe3O4 by carbon-saturated iron. The integrated form of the rate equation indicated that the reduction of FeO was of the first order up to 80 pct reduction. Over the initial 30 pct of reduction, the magnetite reaction was also interpreted as a first-order process. The rates of reaction were extremely rapid as indicated by a 30 pct reduction of the liquid magnetite in a time interval of about 1 to 2 sec. Phil-brook and irkbbride&apos; studied the reduction of FeO from a lime-alumina-silica slag by carbon-saturated iron and solid graphite in stationary crucible assemblies. With the use of the differential form of the rate equation, the rate of the reduction reaction was found to be proportional to the second power of the concentration of the reactant FeO for both the slag-metal and slag-graphite reactions. These authors presented a number of comments in an effort to explain the difference between the molecularity of the above reaction and the observed second-order relation. One such argument, recently extended by Wagner,3 was that the rate-limiting step may well be one of transport control. A complete discussion of the proposed mass-transport control mechanism of this reduction reaction will be presented in a later section of this paper. Kinetic studies of the reduction of other oxide species, namely chromous oxide, titania, and silica, have also been reported. McCoy and philbrook4 used rotating crucible assemblies to investigate chromium reduction from Ca0-SiO2-A12O3 slags. First-order kinetic law was obeyed for the slag-metal reduction reaction as determined by the integration method of data analysis. Due to scatter of the data, any dependence of the rate constant on temperature or slag composition was obscured. Concentration-time data for the reduction of titania from blast-furnace type slags under reducing conditions were obtained by Delve, Meyer, and Lander.&apos; The data were too few, however, for a formal kinetic interpretation. Kinetic data for the very slow reaction of silica reduction have been observed by McCoy and Philbrook6 and Fulton and chipman.7 The former experimenters found the reaction-rate constant for silica reduction by carbon-saturated iron to be 20 to 60 times less than the rate constant for chromium reduction. In their work on the reduction of SiO2 in mechanically stirred systems, Fulton and Chipman also found small values for the specific reaction rate and derived a high value for the energy of activation for the reaction. In both of these investigations, the reduction process was assumed to follow first-order behavior. schuhmann8 proposed that the rate-limiting step in this reaction is the diffusion of oxygen from the interface through a boundary-layer film to the bulk metal phase, and Rawling and Elliott9 have reported experimental confirmation of this hypothesis for temperatures below 1600°C. Turkdogan et a1.10 have also concluded that the rate of reduction of silica is a slow process and controlled by the diffusion of oxygen in the metal, but only in the absence of carbon monoxide bubbles at the slag-metal interface. In the presence of bubbles, achieved either by injecting carbon monoxide at the slag-metal interface or by blowing it through the metal and slag, these investigators found a rather rapid reduction of silica which appeared to be controlled by an interfacial reaction involving the de-sorption of silicate ions from the slag-metal interface to the metal phase as silicon and oxygen atoms. In view of the unresolved nature of the kinetics of FeO reduction and the lack of kinetic data for the manganese reaction,

Jan 1, 1968