Search Documents

By J. H. Brophy, M. H. Kamdar, J. Wulff

A constitutional diagram for the Ta-W-Re alloy system is presented. Rhenium dissolves in the complete range of solid solutions between tungsten and tantalum up to 48 wt pct in tantalum 'to about 30 wt pct in tungsten. Intermediate x and u phases only are found. X-ray diffraction and fluorescence, metallography and melting-point measurements were employed. THE constitution diagram of ternary alloys of tungsten, tantalum, and rhenium for temperatures above 1000°C is reported in this paper. Alloys in this system rich in tantalum and tungsten are of interest for high-temperature applications, and those rich in rhenium for electronic applications. No ternary phase diagram appears to be available in the literature. The three-component binary diagrams have been previously established1-3 and the present ternary data are consistent with these. EXPERIMENTAL TECHNIQUES Commercially pure tantalum (99.7 pct), rhenium (99.99 pet), and tungsten (99.9 or 99.95 pct) powders were purchased from Fansteel Metallurgical Corp., Chase Brass and Copper Co., and Fansteel or Westinghouse Electric Corp., respectively. These powders were pressed without binder and arc-melted under helium. The preparatory and analytical techniques employed were identical to those used for the Ta-Re system investigation carried out in this laboratory.324 These included X-ray diffraction and metallographic phase identification, X-ray fluorescent and wet chemical analyses, solidus determinations and heat treatments in a resistance-heated tantalum-filament vacuum furnace. For microstructural examinations standard mechanical polishing techniques were used through diamond laps. Polished speciments of all compositions were etched by swabbing with 4 parts HI?, 4 parts HNO3, and 1 part H20. Since the microstruc-tures of alloys containing one or two phases closely resembled their counterparts in the respective binary systems, typical examples are omitted in this presentation. In all the alloys examined, it was possible to observe the number of phases present. Particularly when there were small amounts of one phase involved, it was possible to identify it only by a combined conclusion based on X-ray diffraction, position of an alloy in a composition series, and the phase rule. A total of 100 different alloy compositions was arc-melted and studied in the as-cast condition metallographically and by X-ray diffraction. Specimens of each of these compositions were heat-treated at 2680, 2530, and 2020°C for 1, 2, and 4 hr, respectively. Those which were treated at 2020 were given a preliminary homogenization treatment at 2500°C for 1 hr. Specimens of 25 alloys were heat-treated at 1200o C for 168 hr. Solidus temperature measurements were made by the detection of incipient melting in eighteen alloy compositions. The metallographic and X-ray diffraction results are summarized in Table I. Graphically, this information appears in the isothermal plots of Figs. 2 and 3. Within the accuracy of the experimental methods, two such plots are sufficient to locate these data, since the results for 1200°C were similar to 2020 and those at 2530 resembled 2680

Jan 1, 1962

By P. Sanchez-Mejorada

GEOLOGISTS who are acquainted with the mines of northeastern Mexico realize that the most important lead-silver replacement deposits of the province are confined to a single formation. But since the formation is several hundred meters thick, the author has determined the relative position of the lead mantos in the belief that this information will prove useful in the search for new mines. This article describes the stratigraphic location of lead-silver replacement deposits known in several mines situated in a region comprising part of the states of Coahuila and Nuevo Leon. The area studied extends for 200 km* (125 miles) from east to west and 70 km (45 miles) north to south. Metallogenic Province: The most important metallic deposits of the region have the following common characteristics: 1) The orebodies occur as manto replacements of calcareous sediments. 2) The horizons that have contributed large tonnages of ore are in similar positions in the corresponding stratigraphic columns (as it will be shown later on). 3) Lead is the most important metal in the ores. Silver is commonly present in small amounts. Occurrence of zinc is dependent on the degree of oxidation of the primary sulfides. 4) A large proportion of the total production has been in the form of oxidized ores. 5) The deposits are located on anticlines. Rocks of Cretaceous age crop out over large areas of the states of Nuevo Leon and Coahuila. forming part of a belt that extends to the south along the eastern Sierra Madre. The sediments are cut by a small number of intrusive stocks ranging in characteristics from rhyolite to diorite. The most prominent structural feature of the region is a series of folds trending northwest-southeast. Faulting is not conspicuous in the region, although a major fault zone is found farther East.' The physiography closely reflects structural features, with mountains formed by anticlines and valleys conforming with intermediate synclines. Narrow valleys have also been developed along the axes of many anticlines, where solution of limy sedi- ments and erosion have proceeded along axial tension fractures. Stratigraphy: Rocks of the Upper Cretaceous are exposed in only a small percentage of the total area. Generally they have been eroded away from the mountains and are covered by recent talus and alluvium in the valleys. Humphrey' has published an excellent paper describing the Cretaceous formations of the Sierra de los Muertos, a short distance west of Monterrey. A brief summary of these formations is given below to convey an overall picture of the stratigraphic sequence that prevails just south of the area under consideration. Upper Cretaceous Indidura: Essentially thin shaly sediments. Cuesta del Cura: Limestone in thin wavy beds

Jan 1, 1959

By A. W. Schlechten, R. P. Abendroth

Jan 1, 1960

By Frederick Laist

AFTER a series of experiments covering a period of about three years, ending from the spring of 19iJ.2 to the spring of 1915, a 2,000-ton leaching plant for the treatment of the accumulated copper concentrator tailing was built and put into operation. During the experimental period, the first step was that of laboratory experiments or beaker leaches. The results on this small scale were so satisfactory that a small operating plant, capable of handling 5 tons of roasted tailing per day, was installed. Again, the results proved satisfactory and an 80­ton plant was built and operated continuously from August until February, 1913. In this plant, full-sized roasting and leaching units were

Jan 8, 1916

By L. E. Lynd

Ilmenite concentrates prepared from beach and dune sands of Florida, India, Australia, and other localities deviate to varying degrees from the theoretical composition of pure ilmenite, and generally contain more TiO2 than the formula FeTi03 would indicate, and some Fe2O3. This variation in composition is due to different amounts of altered ilmenite in the concentrates. (') The alteration has been shown to take place by oxidation and leaching of the iron in ilmenite, to yield residual products consisting of amaphous or very finely crystalline titanium dioxide and ferric oxide. Other studies(2) have stressed that ilmenite may alter first to an amorphous titanium-iron oxide, or to an amorphous mixture of titanium and iron oxides, with later alteration to crystalline titanium dioxide. There has been some speculation as to whether alteration of ilmenite is usually caused by weathering or by hydrothermal pr0cesses.(3,4,5 Three main methods of approach were used to study the mechanism and rate of ilmenite weathering. (1) Etch tests were used to learn the relationship between crystallographic orientation and the rate of chemical attack on ilmenite. (2) Altered ilmenites in New Jersey sands were compared with ilmenite in rocks of the area from which the sands were probably derived. (3) Experiments were carried out to determine at what rate, and by what agents, originally fresh ilmenite could be altered in the laboratory under simulated weathering conditions. ETCH FIGURES ON ILMENITE CRYSTAL PLANES A rather large ilmenite crystal from Sannidal, Norway, about 11/2 x 11/2 x 2 inches, bearing a few natural faces with imperfect surfaces, was prepared for etching by grinding and polishing crystal planes whose Bravais-Miller indices (referred to hexagonal axes) had been determined by X-ray methods to be 00.1, 10.0, and 10.2. Metallographic examination showed that the crystal was homogeneous, and uniform in orientation, except for the presence of twin lamellae which were oriented with their broadest surfaces parallel to the 00.1 planes of the large crystal. Since ilmenite is rather inert to standard etch reagents, strong solutions of hot sulfuric acid and cold hydrofluoric acid were used. Preliminary tests established that better etch figures were obtained using HF than with H2So4 solutions. The etch figures produced on 00.1 and 01.0 planes of ilmenite by etching 3 to 5 minutes in cold concentrated HF solution are shown in Figures 1, 2, and 3. The figures on 00.1 are "etch hillocks", rather than pits, and according to Honess(6), this phenomenon

Jan 1, 1961

By H. A. Domian, H. Biloni, G. F. Bolling

Jan 1, 1965

By R. A. Wattenbarger, H. J. Ramey

The trend in gas well testing has been to rely more on the early-time flow data of drawdown and buildup tests than on stabilized flow tests. The stabilized testing methods often are not adequate for complicated modern gas reservoirs because many of these reservoirs have extremely low permeabilities and the transient flow period of a well test can last for days or weeks. It has become common to augment the flow capacity of gas wells by hydraulically fracturing their producing formations. In deep reservoirs the induced fractures are generally vertical and tend to follow a single plane of weakness. The presence of a vertical fracture at the wellbore complicates the transient flow behavior of a low permeability gas well. The flow is further complicated when turbulence occurs near the wellbore. Russell and Truittl published transient drawdown solutions for vertically-fractured liquid wells. They developed methods of drawdown and buildup testing utilizing these solutions, which were based on numerical simulation. Clark2 applied the basic Russell-Truitt solutions to analyze fractured water injection wells by falloff tests. Field examples were given to substantiate the method of analysis. The Russell-Truitt solutions and well testing methods can be extended to gas well flow. Millheim and Cichowicz3 presented the drawdown equations for *Now with Scientific Software Corp., Denver. ideal gas flow, including the effects of formation damage and turbulence. Two actual field cases were presented, both of which exhibited extensive fractures and turbulent flow. For the pressure range of these well tests, the use of ideal gas equations was acceptable. These field cases were significant in being the first published data showing the occurrence of turbulence in fractured gas wells. We, also, observed such cases in practice. Our purpose here is to extend the theory of fractured gas well testing to the flow of real gases. To determine the effects of wellbore storage and turbulence on well test interpretation, we developed a finite-difference model to simulate well test conditions. Since gas wells usually are widely spaced and have high compressibility, the emphasis has been put on the early transient behavior before the effects of the outer boundary are noticeable at the wellbore. The Mathematical Model The geometry of our mathematical model is similar to that used by Prats.4 The well is centered in a circular uniform reservoir. A vertical fracture of infinite flow capacity penetrates the formation and passes through the wellbore. The wellbore itself is not important. Fig. 1 shows a sketch of the problem. Fig. 2 shows the idealized model. The equation for the flow of real gas in a porous medium with uniform porosity, permeability and thickness is5-7

Jan 1, 1970

By Thomas C. Wilder

Jan 1, 1965

By W. C. Maurer

A study of the mechanics of shear failure of rock under pressure has been made. The transition from brittle to ductile failure occurs when the friction along the fracture surfaces exceeds the shear strength of the rock. A machine has been developed for measuring both the shear strength of rock and the friction along rock fracture surfaces. In this shear machine, the fractures are forced along definite planes; as a result, the rock is stronger than in compression tests where the fractures choose the weakest paths through the specimens. Friction is composed of several forces: (1) the forces to slide irregulurities on surfaces over each other; (2) forces to shear some of these irregularities; (3) molecular attractions across the fractures; and (4) forces required to break bonds resulting from melting and welding at the tips of irregularities on the surfaces. Friction along rock fracture surfaces is not proportional to contact pressure as usually assumed; instead, it increases as the contact pressure raised to the 0.4 to 0.8 power. INTRODUCTION It is difficult to produce tensile stresses under the high compressive stresses present in the earth; therefore, in deep drilling and geological faulting, shear failure predominates. This shear failure mechanism is complicated because rock is a heterogeneous material containing pore spaces, microfractures, elastic discontinuities and other imperfections. A study of shear failure of rock under pressure has been made to obtain a better understanding of this complex mechanism. A machine has been developed and used to measure the shear strength of rock and friction along sheared surfaces. This machine has certain advantages over the compression test — shear properties can be measured under low normal stresses and shear stress and normal stress can be varied independently. The first part of this paper is a description of shear failure of rock under pressure. The role of shear strength and friction along fracture surfaces in the shear failure mechanism will be discussed. The second part consists of a description of the shear machine and a discussion of the shear strength and friction data obtained using this machine. It will be shown that shear strength depends upon the size of the stressed zone, and that friction does not increase linearly with contact pressure as usually assumed. Attempts are made to relate some of the new concepts to tectonic failure and to dril ling. SHEAR FAILURE IN ROCK TRIAXIAL TESTS The shear failure mechanism is illustrated by the triaxial test shown in Fig. 1. In this test, a cylindrical rock specimen is subjected to a confining pressure, PC, and an axial pressure, pa. These external pressures produce normal a and shear T stresses on inclined planes within the specimen equal to1 where 0 is the angle between the specimen axis and a normal to the inclined plane. Fig. 2 shows how triaxial specimens deform as the axial pressure is increased. The pressure is

Jan 1, 1966

By J. H. Swisher

An experimental study has been made of the possible mechanisms for foam stability in the system CaO-SiO2-Cr2O3, where Cr2O3is the foaming agent. The degree of lowering of surface tension by Cr2O3 was determined at 1600 "C for high -silica melts. Measurements were also made of foam stability under standardized conditions. The results of the two sets of measurements then were correlated by applying the Gibbs and Marangoni theories of film elasticity. It was demonstrated that the Marangoni elasticity effect is probably the largest single contributor to foam stability, although the high uiscosity of' silicate melts makes some contribution in controlling the rate of drainage of liquid from bubble lamellae. ThE tendency of metallurgical slags to foam has frequently been observed during steel-refining operations. In the open-hearth furnace,this tendency is considered undesirable because it slows down the rate of heat transfer from the burners to the molten-metal bath. On the other hand, it is sometimes desirable to have a "foamy" slag in the basic oxygen converter.' A more effective blanket over the steel bath is made by a "foamy" slag, and loss in yield of metal due to splashing is minimized. In addition, a flush slag can be removed more easily in this condition. Since slag-foaming problems in the field have been attacked largely on an empirical basis, it appeared desirable to obtain a better basic understanding of the mechanism of foam stability. The experimental work was directed almost entirely toward the behavior of Cr2O3 as the foaming agent in CaO-SiO,-Cr2O3 melts. In the first stage of the investigation, the degree of lowering of surface tension by Cr2O3 was determined at 1600°C. The second stage consisted of measuring foam stabilities under standardized conditions. The two sets of experimental results then were correlated by applying various theories of foam stability. For convenience in interpreting the surface-tension and foam-stability results reported here, a portion of the CaO-SiO2-Cr2O3 phase diagram2 is reproduced in Fig. 1. EXPERIMENTAL 1) Surface-Tension Measurements. Very little information is available in the literature on the lowering of surface tension of slags by foaming agents. Kozakevitch 3 mentioned that the addition of 1.5 pct Cr2O3 lowered the surface tension of FeO slightly. In another system, some data on the lowering of surface tension of foaming slags by P2O5 were reported by Cooper and Kitchener.4 The technique used for the surface-tension measurements was the maximum bubble-pressure method,' using a single tube immersed to various depths in the melt. This method has been used extensively in high-temperature systems, such as in liquid metals, slags, and glasses.6,7 It has the advantage that the contact angle between the liquid and the tube material need not be known. The equation used for calculating surface-tension values from maximum bubble-pressure measurements is the Schroedinger equatioq5 which can be written as follows:

Jan 1, 1964

By S. R. B. Cooke, H. S. Choi, I. Iwasaki

In a previous article1 the function of various fatty acids as collectors for iron ores was reported for the two alternate processes; (a) the flotation of iron-oxide minerals, and (b) the flotation of calcium-activated siliceous gangue. The selectivity of separation was shown to be markedly dependent not only on such common parameters as pH and the levels of collector and activator addition, but also on the structure and the degree of unsaturation of the acids used. The iron ores employed consisted primarily of hematite, goethite, and quartz and as the effectiveness of any one acid on the separation depends on the individual response of each constituent mineral toward the given collector, a comprehensive knowledge of the floatability of the individual minerals is desirable for the interpretation of the results obtained. The mechanism by which a collector ion becomes attached to a mineral-solution interface has been a subject of much investigation. From quantitative studies of the properties of the electrical double layer on mercury2 and on several solid surfaces3,4,5,6 , the adsorption of some flotation reagents at the mineral-solution interface is beginning to be understood. For example, it has been demonstrated from streaming potential measurements made on quartz5 and corundum6 that the interfacial electrical condition is strongly dependent on the pH. Thisoccurs presumably through the following electrolytic reaction at the interface, and in the presence of a flotation collector, controls the floatability of the solids. A concurrent investigation of the flotation behavior of goethite and of quartz carried out at the School of Mines and Metallurgy7 has shown a positive correlation with the electrokinetic behavior of these minerals. The purpose of the present investigation was to collect general information on the floatability of hematite, and to investigate the function of various 18-carbon aliphatic acids on the floatability of hematite, goethite, and of activated quartz. Electrophoretic mobility measurements, simplified flotation tests using a modified Hallimond tube, and contact angle and frothability measurements were employed. The results are of direct application in the interpretation of the batch flotation results on natural ores as

Jan 1, 1961

By F. A. Seward

A1though the expanding Eltran1 configuration has been used in induced polarization exploration for several years, there has been little application of this method to straight resistivity problems. A comparison of results of resistivity exploration by this technique versus those obtained from Wenner configuration surveys reveals a number of advantages in interpretative resolution as well as in field operation. The application of electrical resistivity to ground-water exploration is well known.= In the past a major portion of this work has been accomplished by use of one or more variations of a technique employing the Wenner3 electrode configuration. With the introduction of the frequency domain technique of induced polarization measurements, used primarily in the exploration for disseminated sulfide deposits, there arose the need for an electrode configuration by which potentials could be measured externally to the current electrode spread in order to minimize coupling between current sending and potential receiving circuits. It was also deemed desirable that the configuration would provide information laterally along profile as well as with depth. To reach these objectives the expanding dipole-dipole configuration (which probably first was used in the electrical transient method of exploration and from which the name 'Eltran' was derived) has been adapted to induced polarization exploration. The author has successfully applied the same techniques

Jan 1, 1963

Who's Who in Mineral Engineering 1977 - SME Membership Directory of the Society of Mining Engineers of AIME

Jan 7, 1977

By W. L. Saunders

WHEN work was begun in September, 1910, on the rock excavation for the foundation of the gigantic dam at Valhalla, N. Y., which is to convert Kensico lake into an important storage reservoir of the Catskill Aqueduct system, one of the first steps taken by the contractor, H. S. Kerbaugh, Inc., was the development of a quarry to supply the stone needed in the construction work.The importance of a bountiful supply of suitable rock in an accessible location is apparent when it is considered that nearly a million cubic yards of cyclopean masonry will be contained in the completed structure, not to mention dimension and facing stones, concrete blocks, and mass concrete. The dam structure will be over 1800 ft. long and about 170 ft. above the present river bed. However, for a distance of some 800 ft., it will be carried down to solid ledge rock, which will make the finished structure about 300 ft. in maximum height. It will be about 200 ft. thick at the base and directly beneath the coping it will be 28 ft. thick. This dam will raise the reservoir surface 110 ft. above old Kensico lake. The reservoir capacity will be 29,000,000,000 gal.-equal to New York's total needs for over two months. The contract price for this work is nearly $8,000,000 and is second only in amount to the cost, over $12,600,000, of the Ashokan Reservoir work. Fig. 1 is a view of the dam site.

Jan 2, 1914

By M. Prats

Most of the information available on the heat efficiency of hot fluid injection processes, both water and steam, has been obtained from calculated temperature distributions in the pay zone and adjacent formations. The usual approach has been to write the heat balance equations in terms of the temperatures, and then to introduce whatever simplifications are necessary to help obtain an analytical or a numerical solution. In either case, the assumption has often been made that the vertical thermal conductivity in the flooded interval is infinite, so that the temperatures within the flooded interval are then independent of vertical position. Because it was first made by Lauwerier,l and has been used extensively since, this assumption will be referred to as the Lauwerier assumption. Excellent reviews of the literature on the heat efficiency of hot fluid injection processes have been given by Spillette,2 Flock et al.,3 and Ramney.4 Of course, the most significant contributions take into account the effect of a finite vertical thermal conductivity on the vertical temperature profile. The most general analytic expression for the heat efficiency of a hot fluid injection process is that of Antimirov5 who considers injection of a heated incompressible fluid into a reservoir through an arbitrary number of wells. The rate of heat injection into the reservoir, as well as the injection temperature, is an arbitrary function of time. The geometry of the horizontal flow is arbitrary, although the reservoir is considered to be of uniform thickness and properties and to be of infinite areal extent. Heat transfer within the reservoir is by horizontal convection and conduc- tion, and by vertical conduction. In the formations adjacent to the flooded interval, heat transfer is by conduction in any direction. Thus, the available expression for the heat efficiency due to hot water injection has been developed for rather general conditions. This degree of applicability has not been obtained for other thermal recovery processes. By making the Lauwerier assumption we have been able to obtain expressions for the heat efficiency that are essentially independent of the thermal recovery process, be it steam, hot-water, or underground combustion. Clearly, then, the approach followed here is more restrictive than that of Antimirov5 in that it uses the Lauwerier assumption. At the same time less restrictive assumptions are made about the horizontal heat transfer mechanisms, either in the flooded interval or in the formations adjacent to it. The main difference between our approach and that of other investigators is that the heat balance in the pay zone is expressed for the pay zone as a whole rather than for a volume element within it. The Lauwerier assumption is introduced after the problem is developed along these lines as far as possible. The details of the development of the heat efficiency are given in the Appendix, while the main features of the model are discussed in the next section. Results and their implications are discussed later in a separate section. The Model Certain assumptions have been made to determine

Jan 1, 1970

By N. R. Morrow, C. C. Harris

The experimental points which describe capillary pressure curves are determined at apparent equilibria which are observed after bydrodynamic flow has ceased. For most systems, the time required to obtain equalization of pressure throughout the discontinuous part of a phase is prohibitive. To permit experimental points to be described as equilibria, a model of capillary behavior is proposed where mass transfer is restricted to bulk fluid flow. Model capillary pressure curves follow if the path described by such points is independent of the rate at which the saturation was changed to attain a capillary pressure point. Aemodified suction potential technique is used to study cyclic relationships between capillary pressure and moisture content for a porous mass. The time taken to complete an experiment was greatly reduced by using small samples. INTRODUCTION Capillary retention of liquid by porous materials has been investigated in the fields of hydrology,l soil science,2 oil reservoir engineering,3 chemical engineering," soil mechanics,5 textiles: paper making: and building materials In studies of the immiscible displacement of one fluid by another within a porous bed, drainage columnsl and suction potential techniques2 have been used to obtain relationships between pressure deficiency and saturation (Fig. 1). Except where there is no hysteresis of contact angle and the solid is of simple geometry, such as a tube of uniform cross section, there is hysteresis in the relationship between capillary pressure and saturation.2 The relationship which has received most attention is displacement of fluid from an initially saturated bed (Fig. 1, Curve R2), the final condition being an irreducible minimum fluid saturation Swr. Imbibition (Fig. 1, Curve A), further desaturation (Fig. 1, Curve R), and intermediate scanning curves have been studied to a lesser but increasing extent.2 ,3,9-12 This paper first considers the nature of the experimental points tracing the capillary pressure curves with respect to the modes and rates of mass transfer which are operative during the course of measurement. There are clear indications that the experimental points which describe these curves are obtained at apparent equilibria which are observed when viscous fluid flow has ceased; and any further changes in the fluid distribution are the result of much slower mass transfer processes, such as diffusion. Unless stated otherwise, this discussion applies to a stable packing of equal, smooth, hydrophilic spheres supported by a suction plate5 with water as the wetting phase and air as the nonwetting phase. The discussion can also be applied, in many

Jan 1, 1966

Who's Who in Mineral Engineering - The 1979 Directory of the Society of Mining Engineers of AIME

Jan 7, 1979

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

By R. G. Bentsen, D. A. T. Donohue

The cyclic steam injection process has become the most widely applied and most successful thermal recovery technique in use today.1-4 Normally, steam stimulation is repeated several times during the life of a project. The problem is to select the optimum number and size of steam treatments so that the profits realized over the entire life of the project are maximum. Obviously, operating costs rise directly with the frequency and amount of steam injection. Operating efficiency, on the other hand, can be adversely affected if the amount and frequency of stimulation are insufficient. Thus arises a multistage decision problem; that is, a decision must be made at each sub-period (day, for example), as to whether steam should be injected, and, if so, how much. It is desirable that these decisions be made in such a way as to secure maximum net profit before taxes. This can be achieved by dynamic programming. To use dynamic programming, we must be able to predict not only the production rate-time curve resulting from a given stimulation, but also how changes in stimulation policy will affect this curve. In this study, the effect of steam injection on well performance is estimated on the basis of a simple model that takes into account both steady-state and transient flow. Oil production rate is calculated for both types of flow, and the larger of the two values is used. This allows the calculation to change to radial steady-state flow at such time as this mechanism becomes more effective than the transient mechanism. We propose to show how dynamic programming can be used to optimize the steam soak process with respect to net profit. For this, a mathematical representation of the steam soak technique is necessary. A further objective, therefore, is the formulation of a mathematical model capable of realistically simulating the physical process. Oil Production Response Model The advantage of using a mathematical model is its generality and ease of manipulation. However, the model must be an accurate representation of the physical process, for a solution derived from a model can be no better than the model itself. The question now arises as to how well the model should be expected to simulate reality. In general, the simplest model should be tried first, for additional accuracy is apt to require additional cost and time. Also, the search for a high degree of accuracy may yield only greater complexity and difficulty without yielding significantly better results. With this in mind, we kept our model as simple as possible, yet generally representative of the actual physical process. The number of state variables is limited to three, and stochastic variations in the parameters are ignored. As a consequence, use of this model should be restricted to cases in which the underlying assumptions are not severely violated. In

Jan 1, 1970

By F. Seitz

STAFF: Editor, Gerhard Derge Carnegie Institute of Technology Schenley Park Pittsburgh 13, Pa. Editorial Assistant, M. A. Redmerski Managing Editor, James J. Burke THE METALLURGICAL SOCIETY: John Chipman, President; W. R. Hibbard, Jr., Post-President; C. C. Long, Vice-President; T. D. Jones, Treasurer; R. W. Shearman, Secretary. PUBLICATIONS COMMITTEE OF THE METALLURGICAL SOCIETY: R. Maddin, Chairman; Dennis Carney, J. C. Fulton, J. H. Hollomon, T. B. King, C. C. Long, W. 0. Philbrook, J. D. Sullivan, David Swan, and F.L. Vogel. Review and selection of manuscripts is under the supervisian of the following Publications Committees: Extractive Metallurgy D ivis ion: T. 0. King, Chairman; S. J. Dick inson, R. P. Hermsdorf, R. Schuhmann, Jr., and H. H. Kel logg. Institute of Metals Division: F. L. Vogel, Choirmon; R. L. Smith, Vice Chairman; R.W. Balluffi, J. H. Bachtold, A Servi, C.D.Graham, Jr., C. G. Goetzel, F. E. Joumot, Jr., G. E. Pellissier, Paul Gordon, Jock Washburn, Harold Margol in, A. U. Seybolt, I. R. Kromer, R. A. Meussner, R. W. Guard, M. Harmon, J. N. Hobstetter, Dean Starr, P. Shewmon, R. Gordon, J. 0. Betterton, and E. Morgan. Iron and Steel Division: W. 0. Philbrook, Chairman; F. C. Langenberg, C. W. Sherman, J. Wagstaff, and E. J. Whittenberger. Published bimonthly by the American lnstitute of Mining, Metallurgical, and Petroleum Engineers, Inc., 29 West 39th Street, New York 18, N. Y. Telephone: PEnnsylvania 6-9220. Subscription $20 per year for non-AIME members in United States and North, South, and Central America; $25, foreign; $5 for AlME members. Single copies, $4.00, single copies foreign, $5.00. The AlME is not responsible for any statement made or opinion expressed in its publications . . . Copyright 1959 by the Americon Institute of Mining, Metallurgical, and Petroleum Engineers, Inc. . . . Registered cable address, AlME New York ... Indexed in Engineering Index, Industrial Arts Index, and Chemical Abstracts . . .Second class postage paid at New York, N. Y. , ond Ann Arbor, Mich. Volume 215, Number 3

Jan 1, 1960