Search Documents

By E. H. Crabtree, Neil S. Parker

The custom milling plant of the Eagle-Picher Mining and Smelting Co. is at Sahuarita, Arizona, approximately 20 miles south of Tucson. It is connected by a 2-mile railroad spur to the main line of the Southern Pacific Railroad between Tucson and Nogales. The capacity of the plant is 500 tons per day. The mill is conventional except perhaps in the rather wide variety of ores treated. These vary from straight copper ores, from which only one concentrate is produced, and lead-zinc or copper-zinc ores from which two concentrates are produced, to complex gold-silver-copper-lead-zinc ores, which require the simultaneous production of three or lour concentrates. Custom ores are received in lots varying in size from a few truckloads to several carloads. Similar ores, after crushing and sampling, may be commingled, or individual lots may be milled individually. The ores are first dumped either into one of two 200-ton railroad bins, or into one of the two 75-ton truck bins. From these bins the lot of ore is individually crushed in a 24 by 16-in. jaw crusher to 3-in. size, and then to ½-in. size in a 3-ft Symons shorthead crusher. After crushing to ½-in., the ore is conveyed to the automatic sampling plant. This consists of three Vezin samplers in series, each successively cutting out 10 pct, 10 pct and 5 pct of its feed. After No. I sampler, the material passes through a mixing barrel and is then crushed to -in. in a set of rolls before passing to No. 2 sampler. Between No. 2 and No. 3 samplers, the ore is again mixed in a mixing barrel. The final sample, consisting of one pound per ton of ore, is taken to the laboratory, where it is further crushed to 1/8-in. in a coffeemill, and then riffled down to the size required for grinding in a pulverizer. The capacity of the crushing and sampling plant is up to 100 tons per hour, depending upon the character of the ore. After sampling, the ore is stored in five ore-storage bins, consisting of four 200-ton bins and one 175-ton bin. Each of these bins is discharged by means of belt feeders and conveyor to the ball mills, so that ores may be composited for milling in any required ratio, or milled separately by themselves. Fine grinding is done in one 6 by 4-ft Allis-Chalmers ball mill and one 8-ft by 36-in. Hardinge mill. These two mills both discharge to one 54-in. Akins Simplex High weir classifier. The classifier sands can be returned to either or both mills, so that flexibility in grinding for different ores is easily obtained. Flotation equipment consists of four 66-in. Fagergren machines for copper-lead roughing and six similar machines for zinc roughing. Lead-copper bulk concentrate is cleaned and recleaned in eight No. 18-S Denver machines, and zinc concentrate is cleaned in three 66-in. Fager-grens. Copper-lead differential separation

Jan 1, 1949

By Arch J. Alexander

Coal is produced, in commercial quantities, in thirty-five of the fifty-five counties of West Virginia. Oil and gas are produced, commercially, in forty-two counties. So, you may readily see that coal, oil and gas arc closely related as to production areas. Many of the West Virginia coal mines are studded with oil and gas wells, as drilling has been done in thirty-two of the coal bearing counties. When the drilling of oil and gas wells was placed under the regulatory control of the Department of Mines on June I, 1929, it was estimated, by the State Geological Survey, that 64,000 wells had been drilled in West Virginia and that 34,000 wells were then being operated as oil and gas producers. No records of many of these remaining 30,000 abandoned wells exist and some of them remain a standing menace to safety as our coal extraction continues. Of a total of 13,100 new wells drilled since June 1929, approximately 5200 have been drilled through coal measures. Approximately 7000 wells, which were drilled through coal measures, have been abandoned in the same period of time. Many of these wells were drilled prior to the enactment of the oil and gas laws. Prior to the enactment of the Oil and Gas Law, which established the Oil and Gas Division of thc Department of Mines in 1929, there had been several mine ex-plosions attributable to cutting into un- charted wells and the improper plugging of other wells. The most serious explosion was that of Mines No. 47 and No. 49 of the Consolidation Coal Co. in Marion Co., when one of the so-called plugged wells leaked or seeped through some 80 ft of strata to the coal bed elevation, filling the mine with gas, the gas being ignited, killing a good many men and wrecking two coal mines. Legislation, to combat well hazards to coal mines, was slow in coming. In 1928, R. M. Lambie, then Chief of the Dept. of Mines, got sufficient support from coal, oil and gas interests to get the legislature to pass our present oil and gas laws. These statutes not only protected the coal industry but went a long way toward conserving the oil and gas resources of our state. Before drilling for oil or gas on any tract of land, known to be underlaid with one or more workable seams of coal, the well operator is required to submit a map of the area showing the location of the proposed well, and to send a copy to the Dept. of Mines and a copy to the coal operator. If no objections to the location of the proposed well are made by the coal operator within ten days, a drilling permit and a number for the well are then issued and work may proceed. If, however, the location of the proposed well is objected to by the coal operator, on the grounds that the location of the well will interfere with present mine operations, projected workings, or otherwise cause a dangerous condition, the operator may file a written objection with the Dept. of Mines.

Jan 1, 1949

By B. H. Caudle, R. E. Watson, I. H. Silberberg

The production history of an inverted nine-spot injection pattern was studied with a porous plate model similar to Habermann's. Colored fluids were used so the boundary movements could be photographed. By varying the production rates of the side and corner producing wells, the optimum pattern sweep efficiency can be obtained. Results of the laboratory tests indicate that interpolation of the pattern sweep and conductivity ratio data will permit application of the data to predicting fluid injection operations with mobility ratios between 0.1 and 10 and production rate ratios between I and the optimum values. INTRODUCTION The inverted nine-spot injection pattern in Fig. 1 has some unique advantages over other regular injection patterns (five spot, line drive, etc.) used in secondary recovery operations. It has three producing wells for each injection well, while the five-spot and line-drive patterns require an equal number of injectors and producers. Futher-more, of the eight producing wells surrounding each injection well, four are of the "side well" type which lie nearer the injection well than the remaining four "corner wells". This allows the production rates of these two types of wells to be adjusted independently of each other to obtain better pattern sweep efficiency. This latter advantage may well be the reason why the nine-spot patterns have not been studied as extensively as the patterns containing only one type of injection or producing well. The pattern sweep efficiency of the inverted nine spot is affected almost as much by operating with different producing rates of the comer and side wells as by differing mobility ratios. Possibly because of this added operating parameter, published studies of nine-spot systems have been limited in their coverage of the pertinent variables.'-' The model studies reported here covered the mobility-ratio range normally found in waterflood operations, as well as the variation in well production rates necessary to achieve optimum pattern sweep efficiency. Conductivity ratios — a measure of the change in total flow at constant pressure difference — were obtained during each model run. A porous-plate model with miscible fluids representing the in-place and injected fluids was used in this study. This model technique is the same as that used in both X-ray shadowgraph" and visual techniques4 described in the literature. Using miscible fluids eliminates residual oil in the model, and the permeability to the injected fluid is the same as to the in-place fluid. Thus the mobility ratio in the model is given by the ratio of the viscosity of the displaced fluid to the viscosity of the injected fluid. Because of the non-linear nature of the viscosity vs composition curve for two miscible liquids, the mobility change across the boundary between the two fluids is reasonably sharp. This technique thus represents the displacement of a region of one fluid mobility by a region of possibly different mobility in an isotropic system under steady-state flow conditions. Table 1 shows the dimensionless similarity groups used to scale the model conditions and results to

Jan 1, 1965

By D. B. Jugle, J. B. Darby

The equilibrium solubility limits of the transition elements Ti, V, Cr, Fe, and Co in liquid tin were determined in the temperature range from 827" to 1211°K. Since the equilibrium concentration of Ti, V, or Cr in liquid tin is moderately dilute in the temperature range investigated, the solubility data were employed to estimate the partial molar enthalpy and the partial molar excess entropy of solution for the three solutes in liquid tin. The reference state s for titanium or vanadium were selected relative to the solute in an intermetallic compound. The thermodynamic quantities for chromium were calculated relative to pure chromium. RECENT experiments1 have shown that the rate of solution of several transition elements in liquid aluminum can be qualitatively predicted from a knowledge of the temperature dependence of the solvent-rich liquidus. In order to extend the analysis to binary tin systems,' it was necessary to know the equilibrium solubility of several of the first-long-period transition elements in liquid tin. Since solubility data for the tin-rich portion of binary systems are meager, the present investigation was undertaken to determine the equilibrium solubility limit of Ti, V, Cr, Fe, and Co in liquid tin in the temperature range from 850" to 1200°K. EXPERIMENTAL PROCEDURE A useful technique3 to determine the equilibrium solubility of a solute in a liquid-metal solvent consists essentially of placing a pure solute in contact with a pure solvent and allowing sufficient time for equilibrium to be attained at the temperature of interest. Subsequently, the solvent is removed and the equilibrium solute concentration in the solvent is determined by an appropriate analytical procedure. The technique can be applied very successfully to relatively simple systems such as the copper-germanium system.3 However, in many binary systems of transition elements with tin, several intermediate phases are known to exist and precipitation of the phases from solution, in solubility experiments, can lead to an erroneously large solubility. The difficulty was avoided in the present work by equilibration of the most solvent-rich intermediate phase with the pure solvent over a range of temperatures. The compounds employed in this investigation, and the temperature range over which each compound is in equilibrium with liquid tin are listed in Table I. No stable intermediate phases are reported in the chromium-tin system4 and, therefore, pure chromium was utilized as the solute. The stoichiometric compounds were prepared by arc-melting on a water-cooled copper hearth in a helium-argon atmosphere. The elemental tin employed in.this investigation was 99.998 pct pure, and typical lot analyses of the transition elements were reported elsewhere.' A sample of an appropriate solute, in the form of a button and in sufficient quantity to exceed the equilibrium solubility, and the pure tin solvent were placed in a high-purity recrystallized-alumina crucible. The crucible and contents were placed in a horizontal furnace under a high vacuum (10-6 torr) and held at temperature for periods varying between 24 and 72 hr. A water-cooled quench tank, partially filled with diffusion-pump oil, was attached close to one end of the

Jan 1, 1970

By F. D. Wright

Production of a coal preparation plant is governed by many restrictions, such as the tonnage of different products and blends that can be sold within a given period, capacities and output proportions of the cleaning and sizing units, blending proportions, quality specifications, and costs and prices for the various products. Determination of the tonnage of each product and blend that should be made in order to obtain the maximum profit can be difficult unless a systematic method such as linear programming is used. In this paper the basic method of linear programming is described briefly. The Old Ben No. 9 preparation plant is used as an example to illustrate in detail how eqwations can be written to form a linear programming model of a coal preparation plant. Three sample problems, each requiring 56 or more equations and 63 structural variables, were solved with an IBM 650 computer. Linear programming is one of several mathematical tools used for operations research. It has been applied to many fields and has been used by a number of industries either to maximize the profits of certain operations or to minimize costs. P. B. Nalle and L. W. weeks' have described the use of the method by the Riverside Cement Co. to minimize the cost of blending raw materials to make portland cement. Their problem is to obtain at minimum cost a mix with certain specifications from a number of possible materials which have various costs and various amounts of CaO, SiO2, Fe2O3, and other constituents. In a paper on the use of linear programming by the National Coal Board in England, K. B. Williams and K. B. Halley2 describe how the transportation method, a variation of linear programming, is used to minimize the cost of sending 37 grades of coal from 28 mines to seven central washing plants which produce coal for furnace coke and foundry coke. The various coals have different percentages of volatile matter, moisture, sulfur, ash, and phosphorous, so there is considerable choice in how they can be blended to meet the specifications of the two products. The purpose of this paper is to show how linear programming can be used to maximize the profits of a coal washing plant which produces individual final products as well as blends. The Old Ben Corp. furnished assumed sample data from their Old Ben Mine No. 9 preparation plant for this investigation. However, the data that have been used are entirely the author's responsibility. METHOD OF LINEAR PROGRAMMING Numerous articles and books have been written on the theory and applications of linear programming.3-5 However, since the method has not been widely applied by the mining industry, a brief, nontheo-retical discussion of its basic method seems to be in order. Linear programming (Table I) is used to determine the best possible solution to a number of interdependent activities. It is essentially a method for making systematic selections from a number of possible solutions to determine positively the optimum solution. There may be other equally good solutions but no better ones. Each activity must have linear coefficients and there must be a criterion to judge how good the solution actually is. Linear programming is different from the solution

Jan 1, 1961

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 M. Semchyshen, A. P. Coldren, W. G. Scholz

A survey of the Fe-Mo-B system was made to determine the extent to which boron might affect the microstructure and strength properties of iron-rich Fe-Mo alloys. Seventeen vacuum-induc tion melted ingots were prepared with approximate mo-lybednum contents 01 2, 4, 9, and 19 wt pct and approximate boron contents of 5, 14,50, and 630 ppm by weight. Phases extracted from annealed specimens and identified by X-ray diffraction analysis included E Fe3MO2, Mo2FeB2, and an oxide hazling a diffraction pattern practically identical to that of MnFez04. Three types of alloy strengthening were observed: a) Fe-2 pct Mo-B alloys were hardenable by a bainitic-type transformtion of austenite to ferrile; b) Fe-4 pct Mo-B alloys were hardened by the solid-solution mechanism; and C) Fe-9 pct Mo-B and Fe-19 pct Mo-B alloys were hardened by precipitation of the E Fe3Moz phase. Boron was observed to have a very strong hardening effect in the Fe-2 pct Mo alloys and a mild strengthening effect on the 1200°F properties of the Fe-4 pct Mo alloys. In the Fe-9 pct Mo alloys there was no consistent effect of boron on strength at either room temperature or 1200° F. The Fe-19 pct Mo alloys were so brittle that meaningful tensile or creep-rupture data could not he obtained. No dispersion hardening from borides was recognized in any of the alloys. THE high-temperature properties of iron-rich Fe-Mo alloys were studied by Reiter and Hibbard.' They found that the strength of the alloys increased markedly as the molybdenum content was raised from 0 to 15.9 pct by weight. The authors explained the various degrees of strengthening in terms of a) solid-solution strengthening, b) strengthening from a martensitic-type transformation, and C) precipitation strengthening. The objective of the present investigation was to conduct a survey of the iron-rich portion of the Fe-Mo-B system to determine whether the beneficial effects of molybdenum on the strength of iron can be enhanced by the presence of boron. It was thought that perhaps fine dispersions of molybdenum- rich borides would be found which could improve the high-temperature strength of ferritic alloys beyond the strengths attainable by molybdenum additions alone. Also, it was considered worthwhile to study the effects of boron in alloys which are essentially free of carbon. Most, if not all, of our present knowledge about boron in high-temperature alloys pertains to alloys containing carbon and carbides. EXPERIMENTAL PROCEDURES Alloy Preparation. Seventeen 8-lb ingots were prepared from electrolytic Plastiron, unalloyed molybdenum chips, and boron powder by a vacuum-induction melting procedure employing hydrogen as the main deoxidizer. The alloys were melted in alumina crucibles using 35-lb melts, with each melt being split four ways to achieve varying boron contents within four Fe-Mo base compositions. The results of chemical analyses conducted on chips machined from the chilled ends of the ingots are presented in Table I. The boron analyses were performed by a spectrographic method utilizing primary standards prepared from boric acid. They indicated that in ingots to which no boron was

Jan 1, 1964

By Frank Woolley, John F. Elliott

The high-temperatzrre solulion calorimeter has been modified and an extension dynamic analysis of i/s transient behavior has provided an improved basis for interpretation of the experimental data. The partial molar heats of mixing, Him, of aluminum, copper, and silicon in the corresponding binary liquid-iron alloys haz,e been measured directly at 1600°C. The resulls in kcal per g-alom are: THE high-temperature1,2 solution calorimeter was planned and designed to obtain directly the partial molar heat of solution of solute elements in a variety of liquid metals that might serve as solvents. Results obtained in the temperature range of 1000° to 1200°C were described in the two earlier publications. However, a long-term goal in the program has been to obtain measurements at 1600°C of the heats of solution of elements in liquid iron. Hence, subsequent to the work last reported,2 and based on the experience of the first two studies, the calorimeter system has been modified and operated at 1600°C to obtain measurements on the Fe-Si, Fe-A1, and Fe-Cu systems. This paper reports the results of the work. THE MODIFIED SYSTEM It was clear from the earlier results that the most precise measurements can be obtained with the high-temperature solution calorimeter when the solution time is less than approximately 30 sec. This can be achieved only when the bath temperature, Tb, is above the melting (liquidus) temperature of the solute material being dropped. A second point of importance is that the refractories used earlier liberated considerable water so that, although the oxygen potential of the atmosphere in the system was satisfactory for work with nickel and copper, it was too high for work with iron, silicon, and so forth. Third, for work at 1550" to 1650°C, an improved temperature-sensing unit is needed. The modified calorimeter in its latest form is shown schematically in Fig. 1. The refractory brick and the Kanthal-Super heating elements which formed the furnace for heating the solute addition have been replaced by a smaller resistance-heated tube furnace which is insulated with bubble-grain alumina. A more elaborate gas purification system was installed to maintain a satisfactorily low oxygen potential in the atmosphere of the calorinieter chamber and of the solute addition furnace. As was the case in the earlier work, at the beginning of a run the calorimeter chamber is evacuated and then filled with purified argon so that a total pressure slightly in excess of 1 atm is maintained during the run. When either of the furnaces within the enclosure attains a red heat, 6 pct H2 is added to the atmosphere. Throughout the run, gas from the calorimeter is circulated through a de-oxidation system at the rate of approximately 5 standard liters per min; the changeover time for the gas within the calorimeter enclosure is approximately 1 hr. The essential stages of the deoxidation system are: a copper furnace at 600°C, a drying tube of anhydrite, and a titanium furnace at 800°C. The oxygen potential of the gas stream in the purification system can be monitored at several points by a stabilized zirconia oxygen cell. This cell is invaluable in detecting the development of leaks in the container vessel and the gas circulating system. The installation of the solute-addition furnace allowed either of two types of temperature-sensing units to be mounted on a water-cooled probe. The first was a small optical device using a silicon diode that was

Jan 1, 1968

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 B. S. Gottfried

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By Arnulf Muan, Carlos Landolt

Equilibria involving Cu-Pt alloys, copper oxides, and atmospheres of known ox2gen pressures have been determined at 1000' and 1200 C. The data obtained have been used to calculate activity-composition relations in the alloys. A moderate negative deviation from ideality is found. At high copper concentrations, the results of the present study are in good agreement with those of previous investigators, whereas at low copper concentrations the activities obtained in this investigation are slightly higher than those reported previously. The present study is part of a research program dealing with the determination of phase relations and thermodynamic properties of copper-containing systems at high temperatures. The information available for such systems is very scanty, particularly for copper oxides and their combinations with other oxide components. Platinum is commonly used as a container in studies of oxide equilibria. To evaluate quantitatively the extent of reaction between the oxide and the container, a knowledge of activity-composition relations in Cu-Pt alloys is necessary. Such data are not available for Cu-Pt alloys in the temperature range of particular interest in the present research program, and particularly for alloys low in copper. Studies of various aspects of the systems Cu-Pt and Cu-0 have been reported in previous literature. The system Cu-Pt is known1 to form a continuous solid-solution series above 812.C, whereas ordered phases are present in certain composition ranges at lower temperatures. Thermodynamic data for the system have been reviewed by Hultgren et al.,' who presented free energies of mixing at 65Oo C, mainly derived from electromotive force cell measurements by Weibke and Matthes.3 Enthalpies of mixing at 640°C have been obtained by Oriani and Murphy4 with liquid tin calorime-try. Their results differ considerably from those derived from the temperature dependence of the data of Weibke and Matthes. Thermodynamic data for the disordered alloys have been obtained by Assayag,5 and more recently by Schmahl and Minzl,6 by Myles and Darby,7 and by McCormack, Myers, and Saxer.8 The two latter sets of data7'' became available after the present study was completed. In the system Cu-0+ two oxides are stable at high temperatures, CuO and CuzO+ each having only minor deviations from stoichiometry. Equilibria in the liquidus temperature region of the system Cu2O-CuO have been studied by Roberts and systemsmyth,9'10CuzO-CuO by Vogel and ocher,I1 by Osterwald,12 and recently the the present authors.13 The system is of the eutectic type with a eutectic temperature of 1083.C The composition of the eutectic liquid corresponds to Xcu0 = °.44, and the oxygen pressure of the coexisting gas phase is Po2 = 0.55 atm.13 Equilibrium oxygen pressures for coexistence of solid CuO and solid Cu2O have been determined by Schmahl and Muller,l4 by Gadalla, Ford, and white,l5 by Bidwell,16 and recently by the present authors.13 There is reasonably good agreement among the various sets of data. Free energies of formation of Cu2O(s) from the elements have been derived by Tretyakov,17 by Bidwell,16 by Kiukkola and wagner,18 by Pal'guev and Neuimin," and by Steele and Alcock,20 There is particularly good agreement among the latter three1'-20 sets of data. THEORETICAL BASIS Consider a subsolidus isothermal section of the system Cu-Pt-0, as shown in Fig. 1. Two phase as-

Jan 1, 1970

G. D. Kneip, Jr., and J. 0. Betterton, Jr. (Union Carbide & Carbon Corp., Oak Ridge, Tenn.)—The authors have contributed to the theory of zone melting by considering the effects of the solidification of the final zone on the distribution curves for the finite length bar. However, they do not consider in sufficient detail the restrictions imposed by the particular metals involved. Eq. 4 has the analytical solution which is equivalent to Pfann's' equation for normal solidification of the final zone. For distribution coefficients less than one, it is apparent from this form of the equation that, as the solidification approaches completion, the concentration attains extremely large values which are inconsistent with the density of physical materials. The necessary restrictions then are that the solute concentration cannot exceed the density of the alloy, or more frequently, the solubility limit in the solid phase. Secondly, the distribution coefficient cannot be constant from 0 to 100 pet solute unless the liquidus and solidus coincide over this whole region. It would thus be more correct to impose limits on the solubility such as would be indicated by a typical eutectic diagram. Furthermore, in this case the assumption that the distribution coefficient remains constant is more likely to be realized. Eq. 4 should then be replaced by the following C,,(x) =k/L-x [ ?11 3 C11 (x)Ddx - ?11 4 Cn(x)dx= k[Cn (L-l)/k] [ L-x/l]x-1 L-l = x = L-l [kCenterlie/Cn (L-l)]1/k-1 The typical shape for the distribution curve in the final zone is indicated schematically in Fig. 12 for the simple eutectic case. The solidification proceeds according to Pfann's' expression until the concentration reaches the maximum solubility in the solid phasc. At this point, the concentration changes abruptly to the eutectic composition and remains at this value for the duration of the solidification. The width of the flat region, and the back reflection of this effect into the concentration curve on subsequent passes, depends upon the ratio of the eutectic concentration to the original concentration, upon the distribution coefficient, and upon the number of passes. Similar effects on the concentration curves in the first zone length would be expected in many systems for which the distribution coefficient is greater than one, and where the maximum solubility is not too much larger than the original concentration. The discussers agree with the authors that the assumption that the liquid is uniform in concentration should be considered with caution. Kneip and Better-ton' have shown, however, that in floating zone refining of zirconium, using induction heating, the distribution of iron agrees with theory which assumes a uniform liquid composition. Hence, in this case, the experimentally realized distribution coefficient agrees with the present phase diagram within the experimental uncertainty. L. Burris, Jr., C. H. Stockman, and I. G. Dillon (authors' reply)—The authors are happy to have received the comments of G. D. Kneip, Jr., and J. 0. Betterton, Jr., which provide additional insight into the zone refining process. It is true that the equations developed are inapplicable if the solute solubility in the solid phase is exceeded. However, this restriction on the use of the equations was clearly stated in the paper. A closer approach to

Jan 1, 1957

By J. H. Reynolds

The possibility of using a semiconductor, metal, semiconductor structure as a light detector is discussed. A brief theoretical argument is presented which predicts that this structure should have pho-toresponse characteristics similur to those of a metal-semiconductor diode, but with the magnitude of the current increased by a transistor gain factor. Measurements made at low frequencies on a particulur S-M-S device agree u~th this prediction aizd shoul that significant gain factors are attainable. The current gain is substantially higher than that predicted for the metal base transistor. It is proposed that the operation of these units is more closely approxinzated by the metal-gate model which predicts higher current gains. A qualitative discussion of its high-frequency performance indicates that detection of light modulated at microwave frequencies by a S-M-S deuice is possible. ADVANCES in the technology of chemical vapor deposition'-4 have made it possible to form a S-M-S structure consisting of a thin metallic layer sandwiched between layers of single-crystal semiconductor. It has been proposed5-7 that, if the metal layer is made sufficiently thin, this structure could be used in an amplifier having several advantages over presently available transistors including higher-frequency performance. Calculations8 predict that the maximum gain-bandwidth product is nearly an order of magnitude higher than that of a bipolar transistor having the same geometry. Maximum oscillating frequencies in the microwave range are predicted. Electrical measurements made on S-M-S devices have shown that, although the theoretical frequency characteristics are not presently attainable, current and power gain is obtained. Another possible application of the S-M-S structure is its use in a light detector. This structure is a three-layer relative of the metal-semiconductor diode. Light detection using M-S diodes is the subject of considerable present interest. They are known to have a short response time and the ability to detect light having a photon energy below the semiconductor bandgap. In this paper, the mechanism by which S-M-S photocurrent is generated is discussed. An expression is derived which shows that the response characteristics are similar to those of a M-S diode, but that the magnitude of the current is increased due to the transistorlike behavior of the device. The results of low-frequency response measurements on a particular S-M-S structure are presented. These agree with the theoretical predictions and show that significant gain factors are attainable. A qualitative discussion of the high-frequency characteristics is presented which predicts that detection of light at microwave frequencies by a S-M-S device is possible. THEORY In this discussion, the photocurrent, iT, measured across terminals attached to the semiconductor emitter and collector layers is related to the current pho-toemitted from the metal base layer. The latter will be referred to as the primary photocurrent. The behavior of the primary photocurrent is determined by the characteristics of the two metal-semiconductor diodes. The diode photoemission is illustrated in Fig. 1, which shows the emitter and collector potential barriers. The metal thickness has been greatly exaggerated. As will be discussed later, the actual barrier shape can be much more complex, but this diagram illustrates the essential features of the photoemission process. An absorbed photon having an energy greater than the barrier height can emit an electron into the emitter or collector. Since this process is discussed

Jan 1, 1968

By F. Forscher

Pronounced porosity, decreasing with distance from the fracture surface, is found in the necked region of tensile specimens tested at room temperature or liquid nitrogen temperature. A hydrogen solution treatment followed by a quench prevents pore formation in tests at —196°C, but has little effect on tests made at room temperature. Reduction of area values at —196°C are strikingly improved by the same treatment, but show no improvement in room temperature tests. The experimental results support the hypothesis that strain-induced porosity results with hydride initially present or, in material containing less than 50 ppm H, with hydride precipitated by strain aging. ZIRCONIUM'S advent as an engineering material in nuclear reactor construction has brought forth many investigations concerned with the properties of this metal and its alloys. Some attention has been given to the ductility of the metal.1, 2 Particularly, the effect of hydrogen on the impact strength was noted," indicating a pronounced dependence of ductility on the hydrogen content and its distribution in the metal. Micro graphic:' and X-ray diffraction studies' have shown that a precipitated hydride phase causes embrittlement in zirconium,' which implies a distinctly different mechanism from that causing hydrogen embrittlement in steels."-' In the latter case, it is assumed that atomic hydrogen diffuses to rifts, where it forms molecular hydrogen under such internal pressures that it causes the rifts to enlarge, leading to fracture. Thermo-dynamic studies of the Zr-H system8, 9 also indicate that the internal pressure of molecular hydrogen would be too low to cause any damage to the zirconium at room temperature, and the formation of a zirconium hydride (d phase) is favored." It is noteworthy that this phase has a face-centered-cubic structure and nearly 18 pet larger molar volume than close-packed-hexagonal a-zlrconlum' Other features of the ductility of zirconium are the absence of a typical cup-cone failure upon tensile testing, a pronounced strain-induced porosity in the necked region of a tensile bar, Fig 1, and the associated spongy appearance of the fracture surface A tremendous improvement in ductility at liquid nitrogen temperature can be observed following a heat treatment that retains the commonly present quan- tities (10 to 50 ppm) of hydrogen in solution. Such tests show a typical cup-cone failure and the absence of strain-induced porosity, Fig. 2. The present paper reports results of an investigation into the causes of strain-induced porosity in zirconium. Room temperature tensile fractures of mild steel, aluminum, and titanium have also been examined for strain-induced porosity. In each case porosity could be detected, but the amount was orders of magnitude less than was found in zirconium. Of the three metals, titanium showed the strongest effect. Material, Specimens, and Preparation The greater part of this investigation was done on Westinghouse crystal bar zirconium, are-melted, forged, and hot-rolled. Chemical analyses of the hot-rolled plates are given in Table I, along with those of the Foote crystal bar zirconium used in this work. The specimens were machined from the hot-rolled plate (1550°F) of 1/2 in. thickness and subsequently heat-treated according to one of the treatments shown in Table IS. All specimens were machined round, the longitudinal axis in the rolling direction, with a gage

Jan 1, 1957

By Roland de Wit

Conclusions reached in a paper by weertmanl are amplified in a mathematical and graphical way. It is shown that in a stressed crystal a straight dis-location may be in a position of unstable equilibrium with respect to helix development, and how a small perturbation may start the dislocation off on its helical path. Helix development occurs by diffision; bulk diffusion is necessary for development from a screw, but core diffision is sufficient for development from an edge dislocation. For a mixed dislocation helix development can take place by core diffusion if the helix also climbs; a diffusion reversal takes place when the slope of the helix equals the slope of the Burgers vector. The work done on a crystal by the external forces is graphically shown to be consistent with the direction of motion of the dislocation in helix development. Finally a detailed mechanism for tangle formation from helices is presented. In the course of helix development a glide situation may be reached where a segment of each helix loop lies in a slip plane in which it can expand by glide. Thus is it visualized how a helix can deteriorate into a tangle. A paper by weertmanl proposes dislocation-tangle formation by a helical dislocation mechanism. The present paper amplifies Weertman's arguments in a mathematical and graphical way and presents some further ramifications of the proposed model. The importance of dislocation tangles is discussed in Weertman's paper and the references cited therein. The mechanism of their formation is not yet understood. Weertman proposes that in cold-worked crystals conditions exist suitable for the conversion of straight dislocations into helices and subsequently into tangles. Helix formation occurs by a diffusion process and Weertman proposes that at low temperatures only core diffusion is important. The present paper necessarily repeats some of Weertman's work. Since the mathematical development and notation used are slightly different from his, this was thought necessary for clarity. GEOMETRY OF THE HELIX weertman2 has shown that the equilibrium form of any dislocation line is a helix if it is assumed that the energy per unit length (line tension) of the dis- location is constant. For the special case that no external forces act on the dislocation, the helix reduces to a straight line. We shall assume in this paper that the dislocation has a helical form even when it is not in equilibrium. The parametric equation of a helix is, see Fig. l: where a is the radius of the cylinder tangent to the helix, and 270 the pitch of the helix. The pitch is related to v, the number of turns per unit length (along the axis) by the relation: When a = 0 Eq. [I] gives a straight line in the z direction. The growth of a straight dislocation into a helix is now described mathematically by letting a vary from zero to its final value. The tangent to the helix is given by where is the path length along the helix. This gives the direction of t and since its magnitude is unity we must have So that dcp/ds = (a2 +p2)-"'. The curvature of the helix is given by where n = -(cos + sin) is a unit vector normal to the helix. Note that the vectors k and n always point toward the axis of the helix, Fig. 1. FORCE ON THE HELIX The external force which makes the straight dis-

Jan 1, 1963

By M. A. Torcaso, P. Raimondi

To study mass transport in systems simulating oil recovery processes, different porous media were saturated with a mobile (carrier phase) and a stationary phase. Slugs of carrier phase containing a small amount of solute were displaced with pure carrier phase. By analogy to the chromato-graphic processes, the velocity of the solute can be predicted from a knowledge of the partition coefficient and the saturation provided that equilibrium between the two phases exists. Equilibrium was found to exist for different porous media, solutes and rates. The conditions were varied over the range normally encountered in the laboratory and in the field. The longitudinal dispersion of a solute undergoing interphase mass transfer was also investigated. INTRODUCTION The production of hydrocarbons by gas cycling, enriched gas drive and C02 or alcohol displacement involves, among other factors, relative motion between two phases and compounds, hereafter called solute, which are soluble in both phases. The solute is carried forward by the faster flowing phase at a lower velocity than the average velocity of that phase. Retardation of the solute is caused by chromatographic absorption and desorption in the slower flowing phase and by the degree of departure from equilibrium. At equilibrium the concentration of solute in the two phases can be related by the equation* CSw = K Cso ,............(1) where C, and C are the concentration of solute in the aqueous and oleic phases respectively and K is the equilibrium ratio, or partition coefficient. Displacement theories must contain an explicit assumption with regard to equilibrium, i.e., whether the compositions can be related by Eq. 1. The existance of equilibrium depends, in general, on the relative velocity between the phases. Unfortunately, other factors such as gravity segregation and viscous fingering, also depend on velocity. For this reason, whenever effects of rate on displacement were observed,2-4 it was practically impossible to discern what caused them — lack of equilibrium or the factors mentioned above. Equilibrium between phases has been the subject of extensive studies in fields such as extraction or chroma to graphy. It has received only small attention in flow through the type of porous media encountered in oil production.5 For this reason a method was developed which makes it possible to study the movement of a solute as it is affected by rate, type of porous media, partition coefficient and carrier phase, but in the absence of segregation or fingering. The information obtained enables one to determine when the assumption of equilibrium can be made. Briefly, the method consists of (1) saturating the core with a mobile and an immobile phase, (2) injecting a slug made up of the same fluid as the mobile phase and a small concentration of mutually soluble solute, (3) measuring the lag and the peak height of the slug at arrival and (4) correlating these variables with fluid properties such as partition coefficient and mixing constants of the medium. PROPOSED MECHANISM The principles of chromatography are combined with the equation of longitudinal mixing to predict the velocity of a solute slug compared to the bulk velocity and the peak height of a slug.6-9 The equation so obtained is valid under equilibrium conditions only. Therefore, a comparison between experimental and predicted results will give a measure of departure from equilibrium. This work was done with either the oleic or the aqueous phase being immobile. For simplicity, the following development is based on the case where the oleic phase is immobile. However, the treatment is the same in either case.

Jan 1, 1966

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge's method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey's theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey's method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge's6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

By A. O. P. Casbarian

An analytical method is developed to determine the variation of the ultimate lateral resistance of a plate in a cohesionless soil with depth. This analysis is based on a modification of Rankine's classical earth pressure theory and the theory of plasticity as applied to soils. The ultimate resistance is defined as the product of the effective stress at the midpoint of the plate, the area of the plate and a dimensionless variable termed the ultimate resistance factor. This variable has been plotted us the depth ratio: i.e., the ratio of of the depth of embedment us the height of the plate. The resistance of a plate may then be calculated using the values af the ultimate resistance factor from the chart provided or the equation may be programmed for use in an analysis of anchor systems in cohesionless soils. It is emphasized that the analysis is semitheoretical. The theory has been compared with experimental results reported in the literature and results indicate general agreement. Actual field tests are necessary to further verify this theory. INTRODUCTION With the exploration for oil offshore in waters all over the world, it is of importance to determine the behavior of the various soils in relation to their ultimate resistance to deformation. Examples of such problems are the holding capacity of anchors and the lateral resistance of piles. Very little information is available in the published literature on the design and performance of anchors, in either cohesive or cohesionless soils. The analysis developed in this report is the first step in obtaining a solution for the determination of anchor holding capacity in a cohesionless material. DISCUSSION The theory used in the analysis is based on the ultimate strength of the soil and is the maximum resistance developed by the plate against further movement. In such a state the elastic deformations are disregarded in comparison with the plastic deformations. Hence, the plate can be considered as completely rigid. The theory of plasticity determines the three unknown stresses at any point by means of two equilibrium conditions for a small earth element in combination with the failure condition. However, the exact solutions can only be carried out in a few simple cases such as, for example, when the rupture or failure lines are straight (Rankine theory) or with spiral and straight rupture lines (Prandtl theory).l Kiitter 2 derived a single equation expressing the variation of the stress in any given rupture line. To utilize this equation, it is necessary to know the stress in the rupture line at a certain point. Unfortunately, this is difficult to obtain unless the rupture line intersects the free surface at a certain angle or when the earth is cohesion-less and unloaded. A method to overcome this is to consider only the boundary conditions at both ends of a rupture line without investigating the equilibrium of the earth above the rupture line. This assumes that the rupture lines meet the surface at statically correct angles, so that boundary stresses may be determined. As Kiitter's equation furnished a relation between these stresses, the unit earth pressure may be calculated at the point where the rupture line meets the wall. Another way in which Kotter's equation may be applied is in investigating the equilibrium of a soil mass above a rupture line. This method assumes that the failure or rupture line is known and that the boundary stresses in the rupture line at the ground surface can be determined. In this case it is possible to determine the earth pressure from the equations of equilibrium.

Jan 1, 1967

By J. C. Rengel

Why and how the nuclear industry entered the electric power generation business is discussed in terms that nuclear energy was an undoubtedly additional energy resource and that it had promise of becoming a virtually inexhaustible generating power if breeder reactors could be devised. The growth of nuclear power, its total generation capacity, are related to the need of increasing energy requirements. The use of atomic energy for power generation will grow if it can compete economically with other fuels. The rate of growth, the time scale and direction of technical development is dependent on ore supplies. It is predicted that by 1980, 50% of all new power installations will be nuclear stations. INTRODUCTION Ten years ago, the nuclear industry was virtually unborn, the Atomic Energy Act of 1954 had just been signed, and peaceful applications of atomic power was little more than a phrase. Serious discussion of the practical application of nuclear energy to electric utility installations in the West was highly unlikely. The questions: Why and how did the nuclear industry enter the electric power generation business at all?; What is the predicted growth of nuclear power?; What technological changes do we predict?; How will the above changes affect the utilities? will in part be answered by this paper. WHY NUCLEAR POWER? The only justification for the West as for any other section of the country for atomic power is economics. Whatever energy source will provide the lowest generation cost will prevail; the only electric utilities which have shown an interest in nuclear power are those to whom it means lower generation costs. It is true that the Government sponsored the development of this technology long before there was any evidence of better economics; but a Government must take a very long-range view and it is the only type of organization that can afford to look more than twenty or thirty years into the future. The Government looked at projections such as shown in Fig. 1, which shows various rates of energy consumption, and matched these against estimates of fossil fuel resources. This showed that we will exhaust our low-cost fossil fuel supplies within 75 to 100 years and our total known fuel supplies within 150 to 200 years. The decision to invest in nuclear energy rested, therefore, on two arguments: a) It was undoubtedly an additional energy resource; and b) It has promise of becoming a virtually inex-haustable resource if breeder reactors could be devised if, while generating power, the reactors would produce more fissile materials than they consumed. This decision produced an unexpected bonus: nuclear power became competitive with fossil power much earlier than had been expected and, as shown in Fig. 2, forced the coal industry to take a hard look at its own position. I believe that the drop of coal prices shown here is the direct result of the entry of nuclear power into the energy supply business. At this point, it ought to be explained why and under what conditions nuclear power can be cheaper than fossil power. Fossil fuel cost to the utility is composed of two factors: a) Mining; and b) Transportation. As a result, the cost of fossil fuel to the utilities varies from 126, 106 Btu to 306, 106 Btu and higher. Nuclear energy is highly concentrated and the transportation factor becomes insignificant. Today the cost of nuclear fuel to the utility ranges between 156

Jan 1, 1967

By W. T. Cardwell

Physically absurd, triple-valued saturations appear in the straight-forward solution of the Buckley-Leverell equations for the displacement of oil by water or gas. From an engineering viewpoint, the triple value causen no difficulty. It is well known how it may be compensnted in order to obtain physically meaningful, numerical results. From a scientific viewpoint, the question still arises: What did the triple value mean? This paper explains how and why the triple value arose in non-capillary Buckley-Leverett theory. The discussion should serve as a background for the understanding and use of the modern method of characteristics in di.~pkzc(>ment theory. INTRODUCTION Displacement theory was introduced to petroleum technologists in 1 941 by Buckley and Leverett.' Their first equation, which they wrote down without derivation, was equivalent to the following: ?s/?t= - qr/F ?f/?r (1) where S = saturation of displacing fluid (a fraction) r = time q,. = total volumetric velocity of both the displacing phase and the phase being displaced (usually, but not necessarily, assumed to be a constant) F = porosity (a fraction) f - fraction of the total volumetric velocity that is the volumetric velocity of the displacing phase (assumed to he a function of S only x = distance. Buckley and Leverett called Eq. I "a material balance equation". Actually it is derived from both a continuity equation, or material balance equation, and Darcy's law. Buckley and Leverett jumped from their Eq. 1 immediately to their Eq. 2, which was the equivalent of the following: (?x/?t)x = qr/F df/ds = qr/F f'(s) (2) The jump from Eq. 1 to Eq. 2 involves apparently simple mathematics, but in that apparently simple mathematics lies a subtle point that is part of the key to the meaning of the triple value to be discussed here. Eq. 2 may be readily interpreted as saying that a point of constant saturation (on a saturation-vs-dis-tance curve) moves with a constant velocity that is proportional to the total volumetric rate, inversely proportional to the porosity, and is otherwise a function of the saturation itself. So that if one knows the derivativc with respect to saturation of the fractional flow function, f(S), for each saturation, one knows the velocity of each point of the moving saturation-vs-distance curve. The function f(S) can he determined experimentally and its derivative f'(S) can be calculated. It turns out that the experimentally derived f'(S) is not a monotonic function of S, but instead it has a definite maximum and declines away at both high and low saturations. This in turn means, in accordance with Eq. 2, that on a moving saturation-vs-distance curve, ccrtain intermediate saturations will travel faster than the saturations either higher or lower. The result is indicated in Fig. 1. In Fig. 1, the abscissas represent distances along a column of porous medium, which column is assumed to be uniform in cross-section perpendicular to the ab-scissal direction. The ordinates represent fractions of pore space occupied by the displacing phase, which may be either water or gas. It is assumed that the saturation is uniform over all planes perpendicular to the abscissal direction. The saturation is a function only of time and one space variable. Time variation of the saturation is represented by change in shape of a curve such as curve AGJ, which represents the initial saturation-vs-