Search Documents

By H. N. Hall

Reservoir and producing characteristics can govern the decision to use either a one-, two- or three-dimensional model for making predictions for gravity-drainage reservoirs. Examples of conditions requiring one-, two- and three-dimensional calculations are given. In 1961 the author presented a method for predicting one-dimensional gravity-drainage performance. This work has been exrendcd to obtain a three-dimensional model, which is described and for which a sample problem is presented. Introduction Many papers on gravity drainage or gas-cap drive'-3 have been presented for conditions where the reservoir is treated as a single areal segment so that, in effect, only one-dimen-sional flow is considered. Other authors'-Q ave given evidence showing that the one-dimensional concept would be unsatisfactory to predict reservoir performance for certain types of reservoirs. Experience has indicated that specific reservoir and producing conditions govern whether one-, two- or three-dimensional concepts must be used as a basis for reliable reservoir performance predictioils. Previously the author presented a paper outlining a method for predicting the performance of gravity-drainage reservoirs.' That paper treated the reservoir as a single areal segment so that it was applicable for only one-dimensional flow. The work has been extended to obtain a three-dimensional model that can be used to predict performance of gravity-drainage reservoirs. Although capable of being three-dimensional, it can be used equally well for one- and two-dimensional predictions. This paper discusses those reservoir and producing conditions governing the choice of model, and also describes a three-dimensional model suitable for predicting reservoir performance in the more complex situations. Discussion The Choice of Model To predict reservoir performance, it is necessary to represent the reservoir by a physical model. A mathematical model, based on the assumed physical model, then is used to predict the movement of oil and gas within the reservoir. The word "model" is used here to connote the combination of both a physical model and a mathematical model based on the physical model. Certain types of reservoir situations may be predicted suitably with a one-dimensional model, whereas two- and three-dimensional models may be required for other situations. This section illustrates reservoir conditions associated with the use of a specific type of model. Many types of reservoir structures are such that gravity drainage possibly could be an important factor in oil production. Production practices, as well as the nature of the structure. can influence the type of model that should be used. obviously, it is impossible to discuss all combinations of reservoir structure and production practices that might be encountered in gravity-drainage reservoirs. The following are typical examples of conditions where a choice must be made among a one-, two- or three-dimensional model to predict reservoir performance satisfactorily. Conditions Requiring a One-Dimensional Model Fig. la shows contours for an asymmetrical anticline in which gravity drainage might be expected to be important. Fig. Ib depicts a cross-section along Line A-A in Fig. la. Existence of a gas cap is a good clue that gravity drainage will be an important factor in producing oil from this type of reservoir. However, it is not essential since a gas cap will form as oil is produced from the reservoir. If the proper combination of reservoir permeability and withdrawal rates is present, the reservoir will behave like a

Jan 1, 1969

By H. W. Mossman

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

Jan 1, 1957

By W. S. Morris

THE East Texas oil field is the largest in the United States and perhaps the largest 'in the world ; likewise, it is one of the most interesting. The East Texas oil field is a water-drive field. Small amounts of salt water were produced from wells on the western side of the field from the bt-ginning of the development. By the end of 1933 (the field having been discovered in September 1930) water production had reached 9000 barrels per day. By the end of 1936, approximately 2000 wells were producing about 45,000 barrels of water per day, and the problem of salt water disposal attracted the serious attention of the operators. Salt water production in the field has steadily increased and today approximately 450,000 barrels of water per day is produced while the oil production averages about 375,000 barrels. The disposal of this salt water has therefore become one of the major operating problems. Various methods of salt water disposal have been tried in the East Texas field. Evaporation in open pits proved

Jan 1, 1944

By AIME AIME

THE 139th meeting of the Institute will be held Feb. 17-20, in the Engineering Societies Building, New York. The annual business meeting will be held on Feb. 18 at 10 a. m., and no technical sessions will be scheduled at that time, so that all members may participate in the election of officers, hear the presentation of the annual reports of officers and committees and discuss any matter which they may care to bring up at that time. The annual dinner will be held at the Conlmodore Hotel on the evening of Feb. 20, and the usual good time is promised. Arrangements for this and other social features of the meeting are, as usual, in hands of the New York Sec- tion, of which the chairman this year, Ralph M. Roosevelt, has a well-established reputation earned as president of the Zinc Institute, as host.

Jan 1, 1930

By Harold Margolin

The Ti-Al system certainly merits the investigative attention it has been receiving and this latest contribution by Blackburn is therefore to be welcomed. The titanium-rich end of the phase diagram shown here, as well as the diagram by crossley,17 Clark et a1.,18 and Tsujimoto and Adachi,29 reveal no other intermediate phase than Ti3A1, although a previous work by Ence and Margolin20 had reported an additional intermediate phase. Since then, Buerger Precession pictures of single crystals of 10.5, 11.5, and 12.5 wt pct A1 were made at New York university21 and the extinctions were found to be typical of a D6h a Ti structure. Some superlat-tice spots and/or streaks suggesting the presence of an Mg3Cd structure, i.e., Ti3Al, were also found. Since these specimens had been quenched from the y region20 and had revealed a Ti, it was concludedz1 that y does not exist and the structures interpreted as m + y could not be so considered. The absence of normal two-phase behavior in these structures had been indicated in the paper. Accordingly, the two-phase region y + 6 becomes a + Ti3Al and there is some agreement between the results of Blackburn and of Ence and Margolin20 as to both shape and existence above 900°C of an a + Ti3Al field. Crossley 17 proposed that Ti3Al is a line compound which decomposes congruently at 900°C. It is difficult to discount entirely the fact that clear two-phase structures, varying in amount with position in the two-phase field, could be detected above 900°C by Ence and Margolin.20 Precipitation during quenching could not have been so regular and the twin-phase syndrome17 does not apply. The heat treatment and quenching technique was essentially the same as was used by Crossley, whereas Blackburn, who used a rapid quenching procedure, also showed a two-phase a + Ti3A1 field above 900°C. If the two-phase field exists above 900° C, it is easy to picture a two-phase field which could account for the data of Crossley below 900°C without invoking a line compound existence of Ti3A1. All that would be necessary would be an 0 + Ti3Al/Ti3Al boundary which curved sharply to higher aluminum contents below 900° C. It is clear that the controversy must be solved by techniques which obtain data at elevated temperatures without quenching. It is conceivable that the difference in width of the a + Ti3N fields reported by Blackburn and by Ence and Margolin is due to the longer annealing times used by the latter investigators. The criterion for establishment of equilibrium in our work20 was the

Jan 1, 1969

By G. K. Manning, C. H. Lorig

TWO metallurgical tools have acquired wide use within the past several years as a means of studying the transformation characteristics of steel. One is a technique used first by Bain and Davenport for determining the transformation characteristics of a steel at constant temperature; the second is the end-quench hardenability test. Constant-temperature transformation curves like those of Bain and Davenport have done much to rationalize the metallography of steel; however, in commercial practices, constant-temperature transformation is seldom encountered, therefore only rarely can such curves be directly applied. The real reason for the determination of such curves lies in the fact that they provide a way of thinking about transformation that is both convenient and productive. The end-quench hardenability test is a more practical means of studying transformation, for, by correlating positions on the hardenability bar with shapes of different masses, an almost direct application of the results can be made. Unfortunately, the end-quench test, though practical, can do very little toward improving the metallurgists' way of thinking about transformation. These two metallurgical tools, dealing with precisely the same phenomenon, should be very closely associated in metallurgists' minds. Actually, they are frequently pushed into separate mental corners and considered as being more or less distinct and unrelated. The reason for this is that no strong connecting link between them has been developed. An attempt is sometimes made to relate these two tools by projecting the critical cooling rate on the isothermal transformation plot as illustrated in Fig. I. The critical cooling rate is considered as a constant rate-though in practical applications, cooling rates are never constant—and this constant cooling rate is given uniform significance from the Ae temperature down to the temperature of initial transformation. A little thought indicates the errors of such a procedure. Certainly the time interval spent just below point A of Fig. I cannot be as effective in promoting nucleation as is the time interval spent just above point B. Yet the procedure represented in Fig. I considers all the time spent between temperatures A and B to be just as effective in promoting nucleation as though all the time were spent at temperature B. Furthermore, change the shape of the isothermal curve between points B and C in any way desired, and the required critical cooling rate is not changed one iota. This cannot possibly be true. Such a way of thinking about critical cooling rates not only is incapable of giving any results that are quantitatively satisfactory but is likely to give very misleading qualitative results.

Jan 1, 1947

By G. K. Manning, C. H. Lorig

TWO metallurgical tools have acquired wide use within the past several years as a means of studying the transformation characteristics of steel. One is a technique used first by Bain and Davenport for determining the transformation characteristics of a steel at constant temperature; the second is the end-quench hardenability test. Constant-temperature transformation curves like those of Bain and Davenport have done much to rationalize the metallography of steel; however, in commercial practices, constant-temperature transformation is seldom encountered, therefore only rarely can such curves be directly applied. The real reason for the determination of such curves lies in the fact that they provide a way of thinking about transformation that is both convenient and productive. The end-quench hardenability test is a more practical means of studying transformation, for, by correlating positions on the hardenability bar with shapes of different masses, an almost direct application of the results can be made. Unfortunately, the end-quench test, though practical, can do very little toward improving the metallurgists' way of thinking about transformation. These two metallurgical tools, dealing with precisely the same phenomenon, should be very closely associated in metallurgists' minds. Actually, they are frequently pushed into separate mental corners and considered as being more or less distinct and unrelated. The reason for this is that no strong connecting link between them has been developed. An attempt is sometimes made to relate these two tools by projecting the critical cooling rate on the isothermal transformation plot as illustrated in Fig. I. The critical cooling rate is considered as a constant rate-though in practical applications, cooling rates are never constant—and this constant cooling rate is given uniform significance from the Ae temperature down to the temperature of initial transformation. A little thought indicates the errors of such a procedure. Certainly the time interval spent just below point A of Fig. I cannot be as effective in promoting nucleation as is the time interval spent just above point B. Yet the procedure represented in Fig. I considers all the time spent between temperatures A and B to be just as effective in promoting nucleation as though all the time were spent at temperature B. Furthermore, change the shape of the isothermal curve between points B and C in any way desired, and the required critical cooling rate is not changed one iota. This cannot possibly be true. Such a way of thinking about critical cooling rates not only is incapable of giving any results that are quantitatively satisfactory but is likely to give very misleading qualitative results.

Jan 1, 1947

By R. E. Goodman, M. A. Mahtab

The, state of stress around a vertical wellbore in rock following nonlinear stress-strain laws is examined by means of finite element analysis. The wellbore, is considered an axisymmetrier body with axisymmetric loading. The initial vertcal and borizontal stresses are "locked" in the rock elements around the wellbore and a new state of stress is generated by the displacements which occur around the borehole. A point-uise variation of the elastic moduli is made on the basis of the new stress state and the triaxial data. The initial stresses tire now reintroduced along with tile changed moduli and original boundary constraints. This procedure is repeatrd until convergent stresses are reached. The effect of nonlinearity on streases is examined for a 6,000-ft wellbore in a scbistose gnezss and Berea sandstone using results of laborratory triaxial compression tests. The results sbou that the effect is restrited to one well radius from the bottom periphery of the bole. Beyond a distance of one-quarter radius, the elfect of nonlinearity on stresscs IS almost always less than percent for the cases conszdered, The consideration of a static pressure inside the well does not magnify the effect of monlinearity on borehole stresscs. INTRODUCTION The terms " wellbore" and "borehole" here designate cylindrical openings in the ground with vertical axis and a circular cross-section. A knowledge of the stress redistribution that occurs on excavating a wellbore is important in understanding (1) the behavior of the lined or unlined hole, (2) hydraulic fracture response, and (3) the effect of stress redistribution on drillability; also it is important in predicting initial stresses in the virgin ground, and in analyzing the response of measuring instruments placed in the borehole. Our knowledge of the state of stress around a wellbore has been restricted to homogeneous, isotropic, elastic material and derives chiefly from the analysis by Miles and Toppingl and the photoelastic work of Galle and Wilhoit2 and Word and Wilhoit.3 In this investigation the state of stress is examined for a nonlinear elastic material by means of finite element analysis. Many rocks possess stress-strain curves that depart notably from straight lines in their initial or final portions. While the literature contains abundant stress-strain data from triaxial tests (axisymmetric loading) on cylindrical rock specimens, there is little information on rock deformability under nonaxisymmetric loading conditions such as occur at each point around the bottom of a wellbore. Although there is some knowledge of the effect of intermediate principal stress on rock strength, there is virtually nothing known about its effect on rock deformability; therefore, we have assumed here that the effect of intermediate principal stress can be ignored. A schistose gneiss4 and Berea sandstone5 were selected as representative rocks for this analysis. The traditional graphs of deviator stress (ol -o3) vs axial strain were reworked to give the tangent modulus as a function of the deviator stress for varying values of the minor principal stress. The result is a nesting family of skewed, bell-shaped curves for tile gneiss (Fig. 1A) and the sandstone (Fig. 2A). A similar replotting of the lateral strain data defines the variation of Poisson's ratio (v) with the deviator stress and confining pressure. These curves, shown in Fig. 1B for the gneiss and in Fig. 2B for the sandstone, are not so well ordered as the tangent modulus curves. However, all of these display an increase of v with deviator stress application, but the rate of increase diminishes with confinement. The ET and v curves for the two rock types are tabulated in Tables 1A and 1A for use in a digital computer so that material properties corresponding to a given state of stress can be assigned by interpolation. STATEMENT OF THE PROBLEM Given a state of stress in the ground prior to drilling a hole (initial stresses), we assume that

Jan 1, 1969

By D. C. Seidel, M. C. Fuerstenau, R. J. Roman

The existence of a surface charge on the alkali halide salts is shown, and a flotation mechanism involving surface charge and collector solubility is developed. Data suggest that KC1 is positively charged and that NaCl is negatively charged in saturated brines, and that the neutral molecule RNH,Cl is the collector in the amine-KC1 system, RSO, is the collector in the sulfonate-KC1 system, and RCOONa) is the collector in the fatty acid-NaCl system. A number of hypotheses have been advanced to explain the mechanism of soluble salt flotation; one of these involves exchange between collector ions and surface cations. With this mechanism Gaudin1 suggested that the differences in flotability of halite and sylvite with primary amines are due to the ionic size of the RNH;, Kt, and Na t ions. Fuerstenau and Fuerstenau expanded this hypothesis to include many other alkali halide salts. Several exceptions to this hypothesis are known. That is, this model does not explain why salts containing the same cation, or why salts of different degrees of hydration, respond differently to flotation. For example, why NaI is floated with an amine whereas NaI .2H,O is not. Further, no correlation between ionic size and flotation response can be made with sulfonate or car-boxylate as collectors. In addition, this hypothesis does not explain the effect of magnesium or lead activation, nor does it offer any explanation of the role of temperature in salt flotation. Rogers and Schulman3 realized some of the inadequacies of the ion exchange mechanism and proposed another mechanism based on the heat of solution of the salts. They concluded that weakly hydrated salts can be floated with alkyl amines and sulfonates, that strongly hydrated salts can be floated with fatty acids, and that excessively hydrated salts cannot be floated at all. Their work was not limited to the al-

Jan 1, 1969

By R. W. Heckel, R. L. DeGregorio

The DeHoff method of determining the size distribution of ellipsoid-shaped, second-phase particles has been applied to the spheroidization of cementite in a eutectoid steel. The surface area of Precipitates determined from the various size fractions in the distribution was correlated with surface-area measurements based on the number of intersected interfaces on a random lest line. The precipitate particles were found to be oblate ellipsoids with an axial ratio of about 0.90. The size distributions were found to be log-normal. A method is proposed whereby the shrinkage of small particles and growth of large ones can be determined from the experimental data. The experimental data are compared to previously proposed mathematical models describing diffusion-controlled kinetics and various types of interface-controlled kinetics. The experimental growth and shrinkage rates are considerably slower than those predicted by diffusion-controlled kinetics. The best fit is obtained for a model describing interface-controlled kinetics limited by the rate of formation of cementite at the growing interfaces where the interfacial reaction rate is proportional to the solute thermodynamic activity gradient across the surface. THE subject of spheroidization of second-phase particles has been considered previously by other investigators. Livingston1 has studied the precipitation of the cobalt-rich phase from copper alloys containing 0.7 to 3.2 wt pet Co. His results indicate that the average particle size increases as the cube root of the heat treatment time in accord with diffusion-controlled kinetics.'-' Komatsu and rant' in their studies of the growth of SiO2 in a dispersion-strengthened copper-silica alloy found that the initial growth of SiO2 proceeded by diffusion-controlled growth and later stages were limited by interface-controlled growth. The values of the activation energy obtained for the interface- controlled process led them to the conclusion that the process was limited by the dissociation of SiO2 at the shrinking interface. The transition from diffusion-controlled growth to interface-controlled growth was characterized by a particle size which varied with the heat-treatment temperature. It is also interesting to note that the particle-size distribution as found by Komatsu and Grant exhibits log-normal behavior when replotted on log-probability coordinates. Dromsky, Lenel, and Ansell9 have observed the growth of Al2O3 particles having a mean free path of about 1.5 to 13 µ. Their photomicrographs indicate that the Al2O3 particle sizes were larger than those observed by Komatsu and Grant. Dromsky, Lenel, and Ansell concluded that the Al2O3 coarsened by an interface-controlled growth mechanism limited by the solution of A12O3 at the shrinking interfaces. Bannyh, Modin, and Modin10 have studied the spheroidization of a eutectoid steel (0.83 wt pet C). Their spheroidization (tempering) treatments were carried out in the range from 210o to 700°C for times between 1.5 sec and 20 hr. Measurements of mean particle size as a function of time indicated a cube root of time dependence of size at 700°C, in agreement with previous analyses of diffusion-controlled kinetics.2°7 At lower temperatures, the time dependence was less than the one-third power. It is important to note that, although mathematical models of growth* considering particle-size distribution have been available, measurements of only mean particle size have been carried out. DeHoff11 has presented a quantitative metallography technique which is applicable to the determination of the size distribution of ellipsoids of constant shape. This method is applicable to both oblate and prolate ellipsoids of all ratios of minor to major axis and is based upon an extension of Saltykov's analysis12 of the distribution of spheres of varying size. DeHoff's analysis is based upon measurements of the minor axis of the particles on a random plane of polish in a unit area. This method provides a measure of the number of ellipsoids in various size ranges per unit area. As pointed out by DeHoff, such measurements may be used to obtain surface area and total precipitate volume data. Thus, the accuracy of the distribution analysis may be checked by comparison of the surface and volume

Jan 1, 1965

By L. D. Mullins, E. B. Elfrink, W. L. Wahl

In the past few years several articles and papers presenting results of solution gas-drive depletion calculations have appeared in the lit-erature. Such calculations are of interest to the oil industry, for investment decisions must often be made before much is known about a reservoir. At other times, an estimate of the possible benefits to be realized from alternate production methods is desirable, and theoretical depletion calculations can serve as a floor or reference level from which to work. In any case, an estimate of ultimate oil recovery based upon engineering data is commonly required. An engineer confronted with the problem of obtaining, for a specific reservoir system, an estimate of ultimate oil recovery by solution gas-drive depletion usually will be forced to perform the calculations himself. This is despite the quantity of data in the literature. Rarely will either experience or the literature provide results from a reservoir system similar in all important respects to the one under consideration, and calculated results are not so plentiful that satisfactory interpolation procedures can be devised. Performing the calculations, however, is a tedious, time-consuming task unless an electronic computer is available, and, in practice, time and manpower are not always available for this purpose. A quick, simple, consistent method was needed for reducing the uncer- tainty in estimated oil recovery from solution gas-drive reservoirs when only minimum information about the reservoir system is available. PROCEDURE Method of Calculation The usual requisite assumptions were made so that the material balance equation could be used to calculate data for the charts. The following assumptions were made: (1) the reservoir is homogeneous and isotropic; (2) oil recovery is due entirely to solution gas drive and neither a gas cap nor a water drive nor gravity drainage is present; (3) the initial reservoir pressure is the bubble-point pressure of the reservoir fluid; (4) initial total liquid saturation is 100 per cent of pore space; (5) interstitial water saturation remains at the initial value as the reservoir pressure declines from the bubble-point pressure to atmospheric pressure; (6) equilibrium gas saturation is 5 per cent of pore space; and (7) oil and gas saturations are uniformly distributed throughout the reservoir at all times. There are no saturation gradients due to a wellbore, nor is the geometry of the reservoir system considered. The material balance equation was written in the form of a differential equation' which was integrated to determine the change in oil saturation for an assigned pressure drop. Formal integration was not possible, so recourse was made to the Runge-Kutta method' of numerical integration. All computations were performed on IBM equipment. Numerical integration yielded the change in oil saturation within the reservoir as the pressure. declined from the bubble-point pressure to atmospheric pressure. The initial oil saturation minus the change in oil saturation yielded the oil saturation at atmospheric pressure. The oil originally in place was obtained by dividing the initial oil saturation by the initial formation volume factor (differential liberation) while the oil in place at atmospheric pressure was obtained by dividing the final oil saturation by the formation volume factor at atmospheric pressure. U1timate oil recovery, expressed as a percentage of the initial oil in place, was obtained by dividing the difference between the oil initially in place and the oil in place at atmospheric pressure by the oil initially in place and multiplying by 100. PVT Data Charts were based upon 135 solutions to the material balance equation. PVT properties of the reservoir fluids were variables in this equation and had to be known as functions of pressure. The required PVT data might have been obtained from either actual reservoir fluid systems or correlated data. However, correlated PVT data were developed and employed in the calculations for the following reasons: (1) it is doubtful that 135 sets of PVT data could have been obtained for the values of variables investigated in this study, and (2) results of recovery calculations from randomly obtained PVT data could not be correlated as well as rssults from selected PVT data. The PVT data used to develop the

By P. W. Rainier

THE Chivor emerald field is situated on the eastern slope of the Andes in the Department of Boyacá, at an elevation of about 8000 ft. above sea level. It overlooks the Llanos (plains) of the Orinoco and is watered by many streams flowing into the Rio Meta, the Orinoco's largest tributary. The Rio Meta is navigable from the Orinoco to a point within about 50 miles of the Chivor district and may eventually prove to be the reasonable outlet for the future development of this part of the country. The topography of the district is extremely rugged and mountainous (Fig. 1), the difference in elevation between two points on a mining claim 1 km. square being sometimes as much as 4000 ft. The ground is covered with dense forest growth and the available timber is suitable and sufficient for mining purposes, but it is becoming gradually depleted by Indian clearings as the country around the mines becomes settled. The many streams combined with the rough topography give good possi-bilities for hydroelectric power. The center of population nearest to Chivor is Bogota, the capital of the republic and a city of some 300,000 inhabitants; the distance is two and one-half days' travel. The first stage of the journey from Bogota, is accomplished by rail to Chocontá, some 4 hr. travel. From Chocontá the journey is made in the saddle, one night being spent in Guateque. The latter is a small town and the nearest post office to the mines, being some 10 hr. ride distant. The roads are mere bridle paths, in parts hazardous and almost impassable in the wet season. All supplies brought to the mines, other than those produced locally, are carried by pack mule; mule freight from Chocontá costs 3 centavos1 per pound. The climate of the Chivor district is healthful and bracing. The wet season begins in March and concludes in November, the remaining months being dry. The annual precipitation is about 150 in. The temperature throughout the year is nearly constant, varying from about 45° F. to about 65° F. shade temperature. The rains are cold and unpleasant, and for the greater part of each day in the wet season the country is shrouded in dense mist. A fire can with comfort be kept burning in a living room all the year round. There are a few slight

Jan 1, 1929

By Lewis T. Wright

IT is commonly believed by metallurgists that in copper-smelting, the copper in the slags, which is irreducible by continued smelting, is retained in the form of "prills" of matte. I have frequently held well-settled slag in a molten condition for a long period without being able thereby to reduce the copper-content. The slag acted as though it contained a minimum of dissolved or combined copper that could not be settled out by gravity. I have used reagents, but without satisfying myself in what form this copper existed. The same slags, by fine grinding and elutriation, could not be separated into portions containing more or less copper than the average content. By treating copper-slags in an electric furnace an impure button of copper. was produced, but the high temperature and the strong reducing-action of the furnace were influences that suggest an explanation of a result not obtainable with ordinary smelting-temperatures. If all the copper in the slag were in the form of copper-matte, and not existing as a dissolved compound with some of the elements of the slag, then the other metals in the matte, such as gold and silver, should occur in the slag in the same ratio to the copper as to the copper in the accompanying matte. This reasoning led me to study the ratio of metals in the products of smelting, and, apart from the issues discussed in this paper, the research has proved both interesting and illuminative. On the assumption of a similar ratio of metals in matte and slag, if the percentage of copper in the slag accompanying; a matte containing 50 oz. of silver and 1 oz. of gold per ton of copper was 0.3 per cent., or 96 oz. per ton of slag, there should be found 0.15 oz. of silver and 0.003 oz. of gold per ton of slag:; but, in my experience, there is less silver and much

Sep 1, 1909

By J. E. Walstrom

While intensive research continues to promote a more complete understanding of the potential and resistivity measurements that comprise the electric log, it is believed that consideration should also be given to translating these numerous and often widely separated findings into a coordinated and readable body of fundamental facts designed specifically for the petroleum engineer and geologist. Although provision is made through publication for a ready exchange of new theoretical concepts. it is also desirable to provide reviews and appraisals of the more established techniques and methods from the operating standpoint so that an economic and practical application may be realized concurrently with the theoretical progress. With these basic premises as a guide the author reviews the presnt state of electric log interpretation. The paper is directed not so much to the logging or research specialist as to the petroleum engineer and geologist to whom the electric log is only one of the many tools which he employs. Frequently, these persons do not have the time to follow in detail the many specialized contributions that appear and, as a consequence. are not in a position to place these contributions in proper relation to each other, or to the art as a whole. The paper reviews the basic steps in making quantitative determinations from the electric log of the amount of oil or gas present in subsurface formations and also discusses the degree of reliability of these determinations under various conditions. The paper also indicates the trend of future developments in electric logging systems and methods of interpretation. INTRODUCTION The electric log has been used about 20 years as a means for studying the formations penetrated by a well bore. The first half of this period is characterized by the development of suitable logging techniques and equipment. Although progress in this direction is continuing at a satisfactory rate, the last ten years are characterized more by an increasing interest in methods of electric log interpretation. During this period, a large number of fundamental papers have been published, expounding various logging techniques and particular phases of the interpretation problem. Many of these papers represent important contributions, and a few are classic. This paper is an effort to outline as concisely as possible and in simple terms the main course of progress in electric log interpretation. More specifically, it is the purpose of the paper to review the necessary elements and basic steps used in making quantitative determinations of water saturation from the electric log; and to point out the degree of reliability of these determinations under different conditions. It is strongly advised that the operating staffs of the drilling and exploration departments of oil companies cooperate wholeheartedly with both the electric logging service companies and research organizations in the testing and development of new logging systems and interpretation methods. One purpose of the paper is. however, to indicate the degree of caution which must be exercised in placing confidence in new techniques and interpretation methods that have not been thoroughly tested in the field. It is entirely possible to be cooperative in trying new methods and yet reluctant to believe in the results until the methods are firmly established. It is important to define the meaning of quantitative electric log interpretation. In the most general sense, an interpretation of the log has been made when the electrical characteristics of the formations, as portrayed on the log, have been translated into terms describing the formation geometry, rock type, or any other physical characteristics of the formations. The determination that the top of a sand is at a certain depth is an interpretation of the log. Structural determinations made by correlating electric logs from a given area are also interpretations of the logs. The term quantitative interpretation, however, will be used in this paper in the restricted sense to indicate the determination of the water saturation of a formation. This determination defines the fluid content of an oil and gas productive formation only if the porosity is known, and it assumes that the remainder of the pore space contains hydrocarbons. This assumption is believed to be true for most oil and gas productive formations. The quantitative electric log interpretation may he said to be a determination of the fluid content only to the extent which the water saturation, under the conditions given above. defines it. THE BASIC STEPS The fundamental steps in calculating water saturation from the electric log are: 1. Determination of the true resistivity of the formations from the apparent resistivities as recorded on the electric log.

Jan 1, 1952

By GEORGE D. RAMSAY

ABOUT three miles west of Fontana San Bernardino County, California, and fifty miles east of Los Angeles, the Kaiser Co., Inc., has built an integrated steel plant. By integrated, I mean that from its wholly owned iron ore and coal supply, the Kaiser Co. is producing in its own mills iron and steel products. Pig iron to meet the West Coast needs has been produced at Fontana since Jan. 1, 1943; steel ingots were first poured there May 14, 1943; and steel plate for ships was first rolled Aug. 15, just sixteen months after ground was broken. Briefly, the plant consists of: (I) Two batteries of by-product coke ovens, totalling ninety ovens, carbonizing about 1740 tons of coal per day at capacity; with by-product facilities for the production of tar, ammonia, and ammonium suIphate, light oil distillates, and tar distillates. (2) One blast furnace producing at

Jan 1, 1944

By C. E. Sims, C. A. Zapffe

MANY hundreds of publications have appeared during the past 78 years that treat the subject of hydrogen in iron and steel,105 but conclusions regarding the functions of hydrogen in causing some important defects in steel are still unsatisfactory to many; and other hydrogen-caused phenomena yet remain to be identified with this gas. The most widely discussed of these defects is hydrogen embrittlement, but the conception, perhaps the popular one, that embrittlement is due to hydride formation does not conform to the known facts.104 The true identity of hydrogen embrittlement seems instead to lie in a mosaic nature of metal crystals, which is a concept not yet accepted by most metallurgists. A study of hydrogen embrittlement therefore provides new viewpoints in regard to the crystalline substructure of steel. A second widely discussed hydrogen-caused defect is the "flake" or "shatter crack." Although it is generally agreed that the presence of hydrogen is a prerequisite for their appearance, it is also realized that internal stress plays a decisive part. The relationship of hydrogen to internal stress has not been explained very clearly; nor is it always clear what the relationship is between a "flake" and a "snowflake," or between a "fish eye" and a "shatter crack." More importantly, the relationships among hydrogen embrittlement, hydrogen-caused fissures such as "flakes" and the ultramicroscopic structure of steel have scarcely been recognized. The present investigation attempts to provide an explanation for the embrittlement of steel by hydrogen by showing cause for accepting the "block" concept of metal crystals, and then to show how the numerous hydrogen-caused defects are interrelated in the light of this concept. EVIDENCE OF SUBSTRUCTURE IN METAL CRYSTALS It appears that the concept of a crystal as an aggregate of smaller multi-atom units must be accepted in some form if the behavior of metals is to be understood. The abundant literature on the evidence of an ultramicroscopic "mosaic" structure in metals and other crystalline substances suggests that the physicist has provided the metallurgist with an invaluable new insight. The exact picturization that any one physicist provides may be open to question, but the fact that there is some such fundamental substructure seems undeniable. The present work makes no pretense of comprehensively treating the imperfection structure theories. Nor will conclusions be drawn here regarding the exact nature of the substructure of crystals. However, the acceptance of some imperfection theory of fundamental attributes seems so necessary to an understanding of the behavior of hydrogen in steel that the following brief

Jan 1, 1941

HE purpose of this chapter is to present a general outline of the basic open-hearth process for the benefit of students, practicing open-hearth operators, and metallurgists who wish to review the subject. The discussion is limited to a description of the sequence and purpose of the operations, illustrated by data and examples. Subsequent chapters deal with more detailed and fundamental studies of the subject. The function of the basic open-hearth furnace is to convert various types of ferrous material into finished steel of given composition and quality. In comparison with other steelmaking methods, the basic open-hearth process is one of the most versa- tile, both from the standpoint of the raw materials that can be used and the types of steel that can be produced. The process can be considered as a sequence involving melting, refining, and deoxidation. With some types of liquid-metal charges, the melting stage is not present, but refining and deoxidation in varying degrees are fundamental features of the process. GENERAL DISCUSSION OF BASIC OPEN-HEARTH PROCESS The first consideration in open-hearth practice is the type of charge on which the operation is based. This consists of three components-fluxing, oxidizing, and metallic materials. Fluxes. In the basic process the fluxing material is primarily lime, added in the form of either limestone (CaC03) or burnt lime (CaO). This is supplemented by magnesia (MgO) and lime from the furnace bottom and banks. The function of these fluxes is to form a slag with the silica (SiO2) , manganese oxide (MnO), iron oxides (FeO and Fe2O3), and phosphorus oxide (P2O5) that are produced by the oxidation of varying amounts of the silicon, manganese, iron, and phosphorus in the charge during melting and refining. The slag, composed of the basic fluxes and the neutral and acid oxides, floats on top of the steel bath when the charge is melted. The composition of this slag depends on the amount of each slag-forming element present in the charge and on the extent of oxidation of these elements during the course

Jan 1, 1951

By V. Leroy, A. G. Guy

Serious disagreements are often found between experimentally determined intrinsic diffusion coefficients and those calculated employing the usual form of the vacancy theory. In the new theory it is proposed that the total intrinsic flux, Ji, of component i, is the sum of a part, f, due to the usual random exclzanges of component i with the vacancies, and a second part, Ji, due to exchanges with the uacancies composing the net vacancy flux. The present treatment, while less powerful than that of Manning, has the advantage of easy uisualization and of facilitating the application of the vacancy-flux effect to complex systems. IT is becoming increasingly evident that there are serious deficiencies in the version of the vacancy theory of diffusion that has been widely used for the past 20 years. One type of evidence is the frequent lack of agreement between intrinsic diffusion coefficients and tracer diffusion coefficients, even taking account of the thermodynamic factor. A second kind of evidence is the observation of a Kirkendall shift larger than theoretically possible, that is, larger than can be accounted for without assigning a negative value to one of the two intrinsic diffusion coefficient.'- The thermodynamic factor could conceivably make both coefficients negative, but not just one. It is clear that a cause of these anomalies, apart from any inadequacy of the usual vacancy theory, might lie in an oversimplified treatment of the data. Adequate experimental techniques, including the use of moderate pressure during the diffusion anneal, are now available to insure that porosity, lateral expansion, and so forth, can be kept negligibly small in most cases. The effect of differences in atomic volume can be of major importance, and it is essential that one of the available methods4 be used to account for this factor. In the present treatment this is accomplished by the consistent use of moles per cubic centimeter as the unit of concentration. Of the various possible inadequacies of the vacancy theory, attention will be given here only to effects of the net vacancy flux. annin' has previously considered this question, beginning with an analysis of atomic jumping of tracer atoms. When he added the effect of a concentration gradient, new terms arose that could be associated with the flow of vacancies. The present treatment uses quite a different approach. The usual vacancy flux, J,, is introduced explicitly, and a simple analysis predicts major changes in the intrinsic diffusion coefficients from this cause. The usual assumptions are made that only a vacancy mechanism is operative, that the formation of voids can be neglected, and that changes in the partial molal volumes, vl and v2, are negligible. The significant diffusion coefficients for the present topic are Dl and D,, the intrinsic coefficients, which enter in the equations, where the flux Ji, moles per sq cm per sec, is that crossing the Kirkendall interface. The concentration, ci, is in units of moles per cu cm, and the concentration gradient, aci/ax, is evaluated at the Kirkendall interface. It will be recalled' that the calculation of Dl and D2 involves the measurement of areas on the diffusion curve with respect to the positions of the Kirkendall and Matano interfaces. In the case of the anomalies mentioned earlier, the Kirkendall shift is too large to be accounted for by the diiferetzce in fluxes (J2 -J1), given by Eqs. [I] and [2]. The logical inference is that the flux of the solvent atoms, J1, is actually in the same direction as the flux of the solute atoms, Jz. In terms of Eq. [I] this requires that Dl have a negative value. However, it would be somewhat misleading to state that the solvent atoms are diffusing up their own concentration gradient. The explanation that will be advanced here pictures competing processes producing the net flux of solvent atoms: 1) diffusion of the solvent atoms down their own gradient by random exchanges with vacancies; and 2) diffusion of solvent atoms in the opposite direction by exchanges with the net vacancy flux. ACTION OF THE NET VACANCY FLUX Theories of vacancy diffusion can be formulated with varying degrees of refinement, and the present theory has purposely been kept as simple as appeared adequate to explain the phenomenon in question. In particular the following aspects have been neglected: 1) the gradient of vacancy concentration in comparison to the gradient of the atomic species; 2) departure of vacancy concentration from the local equilibrium value; 3) variation of the jump frequency, LO, with the specific surroundings of the atom-vacancy pair being considered; 4) correlation effects. These and other refinements can be considered once the essential mechanism has been established. The essential idea of the present analysis is to calculate the total intrinsic flux, Ji, of component i as the sum, JlJ?±j{ [3] where J; is attributable to the usual random atomic jumping, and J{ is a contribution arising from the net vacancy flux, J,. The latter quantity, of course,

Jan 1, 1967

By J. S. Kirkaldy, D. H. Weichert, G. R. Purdy

Two-phase diffusion couples have been used to simulate the growth of proeutectoid ferrite in ternary Fe-C-Mn austenites. It has been shown, theoretically and expermentally, that the results fall into two classes: a high-super saturation class in which the major role of manganese additions is to influence the boundary conditions for carbon diffusion; and a low-supersaturation class in which the partition of manganese is required. In the latter case, a drastic inhibition of the reaction rate is to be expected. As an aid to the kinetic analysis, a portion of the iron-rich comer of the Fe-C-Mn constitution diagvam has been determined. The conclusions drawn from this work have a direct bearing on the kinetics of the grain boundary ferrite precipitation reaction and thus yield insight into the manner in which alloying elements can influence the rate of austenite-decomposition reactions. THE role of alloying elements in the austenite-decomposition reactions has provoked considerable interest, since it is apparent that an understanding of this phenomenon will yield insight into the harden-ability problem. In the present investigation, attention is focused upon the proeutectoid-ferrite transformation in the system Fe-C-Mn, and upon the manner in which manganese additions affect the growth rate of ferrite. When hypoeutectoid austenite is supersaturated, the ferrite first formed usually exhibits two distinct types of morphology. At low super saturations, grain boundary layers grow into the austenite grains with an approximately planar interface, and usually possess at least one incoherent phase boundary. At high supersaturations, Widmanstiitten plates appear, which bear an orientation relationship' to the parent austenite. In this case, there is a chance of partial lattice matching along the sides of the plates and an attendant semicoherent interface structure. Purdy and Kirkaldy2 have used a two-phase binary Fe-C diffusion couple to demonstrate that carbon volume diffusion controls the incoherent reaction, while Rouze and rube,' using thermionic-emission microscopy, have shown that the thickening of Widmanstatten ferrite plates is appreciably slower than expected for carbon-diffusion control. This suggests that the plates are bounded by partially coherent, low-mobility interfaces. The redistribution of alloying elements during transformation has been the object of studies by Aaronson4 and owmman,5 who have shown that no alloy partition occurs during ferrite precipitation in the systems Fe-C-Cr and Fe-C-Mo, respectively. More recently, Aaronson et al.6 have shown that some manganese partition occurs at low super-saturations in the system Fe-C-Mn, but that no partition occurs at higher super saturations. Pickle-simer et al.7 observed that no manganese partition occurred during the early stages of the pearlite transformation in a manganese steel, and concluded that volume diffusion of the alloying element is not the rate-deter mining factor. The diffusion coefficients pertinent to this investigation have been measured by Wells, Batz, and Mehl8 (carbon in austenite), R. P. smith9 (carbon in ferrite), Wells and Mehl 10 (manganese in austenite), and Kirkaldy and Purdy11 (diffusion of carbon on a manganese gradient in austenite). Kurdjumov 12 has shown that the diffusivity of substitutional elements in ternary austenites up-quenched from martensite may be greatly enhanced, due, apparently, to the presence of a persistent defect structure. The recent work of Krauss 19 showing high densities of tangled dislocations in reversed Fe-Ni martensite is in agreement with this conception. THEORETICAL Diffusion in ternary or higher-order systems may be described with the aid of Onsager's extenbion of Fick's first law,14 in which the flux of each of the n components is assumed to be a linear function of all concentration gradients:

Jan 1, 1964

By B. Ramaswami, U. F. Kocks, B. Chalmers

This paper describes the tempetrature, solute (gold), and orientalion dependence of the plastic rleforlirnlion of silver single crystals. The virgin flow stress, To, and its temperatutre dependence incrense with increase in solute concentration. el decreases with decrease in temperature below room teiirperature and with increase in solute coucentralion. The dependence of 011 on temperature arid solute sontenl is due to the dependence of Ts on the two variables. An increase in ts, due to either the temperature or the solute conlent, results in a decrease of e11 in crystals of single-slip orientation and an increase of 011 in crystals of' (100) find (111) orientation. is orientation-dependent and at all the temperatures and in all the materials studied. The Leveling off of the engineering stress-strain curve in crystals of (100) 07-ientotion at a stress tl, is associated with extensive cross slib. The temperature and strain -rate dependence of tl are similar to those of indicating that a rate process is the cause 01 this phenomenon. The plastic deformation of silrer crystals of (110) orientation is multiple to that of crystals of single-slip orientation. THE work hardening of fcc crystals has been the subject of intensive investigations in the past 15 years. The effect of orientation, deformation temperature, solutes, second phase or precipitates, strain rate, etc., has been studied in a number of fcc metals.'-5 In spite of the considerable amount of work reported on the dependence of work hardening on the orientation, temperature, and solute concentration, certain aspects of the plastic deformation of fcc crystals have not been clarified. For instance, with the exception of a few investigations,6-8 where crystals of multiple-slip orientations have been studied, almost all the studies have been on crystals of single-slip orientations. Hence, a major portion of this study has been devoted to the dependence of plastic deformation on the multiplicity of slip systems and the enhanced possibility of dislocation interactions resulting from such multiplicity. Further, the abnormally low rate of work hardening observed at large strains at certain temperatures in crystals of the (100) orienta-tion 6-8has not been explained satisfactorily. It is the aim of this study to clarify this aspect of the orientation and temperature dependence of the plastic deformation of fcc single crystals. The work-hardening characteristics of single crystals of aluminum of ( 110) orientation have been studied at room temperatur "'12 and above." They have not been studied at temperatures where the virgin flow stress is dependent on temperature or where, in Seeger's terminology is important. The work hardening of single crystals of (110) orientation is of interest for another reason. As shown in Table I, the dislocation reactions that are possible in the (110) orientation are similar in certain respects to those that are possible in the (100) orientation. Hence, the deformation of single crystals of silver of (110) orientation was studied in the temperature range 100° to 900°K at intervals of 200°K. The effect of solute concentration on the plastic deformation of fcc single crystals of single-slip orientation has been studied by several investiaators.> In all of these investigations—with the exception of Ref. 20—single crystals of identical orientation have not been used, with the result that the orientation dependence of plastic deformation is superimposed on the effect of the solute. In addition, the effect of the solute on the work-hardening characteristics of single crystals of multiple-slip orientations has not been studied by the previous

Jan 1, 1965