Search Documents

By John Lamont, Walter Crafts

The principle formulated by Grossmannl for calculating hardenability of steel by multiplying the ideal diameter of "pure" iron-carbon alloys by factors for grain size and alloying elements has been confirmed on many steels and has proved useful in designing substitute alloy steels. Grossmann found that the typical behavior of the multiplying factors was to increase directly in proportion to the amount of the alloy present in a steel. However, the multiplying factor for the effect of silicon on hardenability was found by Grossmann to increase to a lesser degree as the amount of silicon was increased. In view of the increasing use of silicon in heat-treating steels, the effect of silicon on hardenability is becoming more important. Calculation of the hardenability of high-silicon steels made at the Union Carbide and Carbon Research Laboratories, Inc., indicated that the multiplying factor for silicon increased directly in proportion to the amount of silicon present. A more extended study of high-silicon steels covering fairly wide ranges of grain size, carbon, manganese, and aluminum has confirmed that up to 2 per cent silicon the multiplying factor is directly proportional to the silicon content. Manganese and aluminum factors have also been! determined. The results of the study offer strong confirmation of the general validity of Grossmann's multiplying principle. Testing Procedure Steels on which the study of the effect of silicon on hardenability was based were made in high-frequency induction furnaces using Armco iron and standard alloying materials. The ingots were forged to bars, and sections were normalized prior to machining Jominy hardenability test specimens. The Jominy test was carried out under standard conditions.2 The depth of hardness penetration for determination of ideal diameter at "half hardness" was taken at the hardness level indicated in Fig. 29 in Grossmann's paper by the solid line for 50 per cent martensite in plain carbon steel. The depth of penetration of hardness on the Jominy specimen was converted to ideal diameter by the relation given in Fig. 28 of Grossmann's paper. Chemical analysis was determined on a sample taken adjacent to the Jominy specimen. Aluminum was determined as "acid-soluble aluminum" rather than as total aluminum, as it is probably more representative of the aluminum that is effective as an alloying agent. In addition to the elements nominally present in the steels, analyses were made for residual phosphorus, sulphur, nickel, copper, molybdenum and chromium on one heat of those made with the same materials during the same period. Actual grain size was determined by microscopic examination of small samples quenched in water after receiving the same thermal cycle as the Jominy hardenability specimens. Chemical analyses, heat-treatment, grain size, and hardenability are shown in Table I representing three ranges of composition.

Jan 1, 1943

By C. V. Millikan

Temperature measurement in wells is an old practice and geothermal gradients have been of interest to geologists for many years.le2 Their application to the operation of oil wells is a more recent practice. It has long been recognized that temperature anomalies occur in drilling and producing wells but thermometers that could obtain a satisfactory record of the anomalies were not available. The normal geothermal gradient is con-sidered as being about one degree Fahren-heit for each 60 ft. of depth.4 This "nor-mal " gradient varies in different areas, but whatever it may be in an area there is little variation from the surface to as deep as the : drill penetrates. The normal gradient is appreciably changed by gas expanding ; from the reservoir into the borehole or by the movement of fluid through a borehole : while drilling, producing, or circulating. It is the resulting anomalies and abnormal gradients from which conditions in wells can be interpreted. Some of the applications of such interpretations are locating the position of oil or gas sands, permeable strata in a reservoir zone, gas-oil contact, source of water, casing leaks, and the top of cement behind casing. Thermometers A thermometer for use in wells must be sensitive and have a low thermal lag. The sensitivity should permit recording 0.2OF. or less for differential temperatures although the actual temperature need not be more accurate than 1°F. The thermal lag should be as low as rugged construction will permit in order to allow surveys to be made in a minimum of time. Of the several kinds of thermometers, these requirements are most easily obtained with the vapor-pressure, expansion, and electric-resistance types. Thermometers for use in wells may be self-contained—that is, the complete instrument is lowered into the well and a record made of the temperature—or the thermal element only may be lowered into the well and the temperature recorded at the surface. The vapor-pressure and expansion types are more adaptable for a self-contained instrument while the electric-resistance thermometer is more suitable for surface recording of temperatures. The self-contained instrument can be run against high pressure on a single wire line with only one operator required. The instrument and its hoisting equipment is easily portable and can be mounted on a small car. The disadvantages are a somewhat greater thermal lag and the necessity of completing a run before any temperature or temperature anomaly is known. A surface-recording thermometer has the advantage of being somewhat more sensitive and of having a lower thermal lag. The Operator can see what anomalies are encountered and can immediately change the running conditions of the instrument or the well conditions in order to obtain the most

Jan 1, 1941

By James E. Berryman

A method is presented for calculating the pressure depletion performance of a gas-condensate reservoir from the original reservoir fluid composition using the laboratory determined pressure-production history and retrograde liquid accumulation as guides. Equilibrium vaporization ratios for this calculation were modified from published data until the calculated dew point pressure and retrograde liquid accumulation agreed with the values observed in the laboratory. The calculated results are compared with the produced gas composition determined by low temperature analysis during constant volume laboratory experiments of pressure depletion behavior. This comparison shows that the calculation procedure is adequate for determination of gas phase composition. Data from field separation tests in low temperature separators are presented which verify the applicability of published equilibrium vaporization ratios for a wide range of separator conditions. INTRODUCTION A reliable prediction of the pressure depletion performance of a gas-condensate reservoir is necessary in determining reserves and evaluating field separation methods. The predicted performance is also used in planning future operations and studying the economics of projects for increasing liquid recovery. Because the flow of fluid in a volumetric gas-condensate reservoir is usually single phase flow, the pressure-depletion per-formake is primarily a function of the phase behavior of the reservoir fluid. The properties of a gas-condensate system of major importance are: (1) compressibility factor at reservoir temperature and pressure; (2) dew point pressure; (3) composition changes of the vapor phase with pressure depletion; (4) recovery of original in-place hydrocarbons at any pressure; and (5) retrograde condensate accumulation. These properties are most generally determined by a pressure depletion-type experiment on recombined separator gas and liquid samples in the laboratory. ldeally the pressure-depletion performance can be calculated from the composition of the original reservoir fluid if the reservoir temperature and initial pressure are known.1,2,3 It has been reported "... comparison of such computations with depletion-type ex-periments on recombined samples in the laboratory have yielded results that appeared to be entirely adequate.". The success of such calculations is determined by the equilibrium vaporization ratios chosen, especially those for the higher boiling point hydrocarbons. It is the purpose of this paper to present a method for computing the change in composition of the vapor phase during pressure depletion by production using laboratory determined properties of the system to guide the calculation and to illustrate the use of the method with data from the Cockfield "D" sand of Washington field, St. Landry Parish, La. Results of field tests and laboratory analysis are included. As a result of this reservoir fluid study, a cycling program has been evaluated, the fidd has been unitized, and the cycling started in July, 1956.

Jan 1, 1958

By John L. Lamont, Walter Crafts

THE principle formulated by Grossmann1 for calculating hardenability of steel by multiplying the ideal diameter of "pure" iron-carbon alloys by factors for grain size and alloying elements has been confirmed on many steels and has proved useful in designing substitute alloy steels. Grossmann found that the typical behavior of the multiplying factors was to increase directly in proportion to the amount of the alloy present in a steel. However, the multiplying factor for the effect of silicon on hardenability was found by Grossmann to increase to a lesser degree as the amount of silicon was increased. In view of the in¬creasing use of silicon in heat-treating steels, the effect of silicon on hardenability is becoming more important. Calculation of the hardenability of high-silicon steels made at the Union Carbide and Carbon Research Laboratories, Inc., indicated that the multiplying factor for silicon increased directly in proportion to the amount of silicon present. A more extended study of high-silicon steels covering fairly wide ranges of grain size, carbon, manganese, and aluminum has confirmed that up to 2 per cent silicon the multiplying factor is directly proportional to the silicon content. Manganese and aluminum factors have also been determined. The results of the study offer strong confirmation of the general validity of Grossmann's multiplying principle. TESTING PROCEDURE Steels on which the study of the effect of silicon on hardenability was based were made in high-frequency induction furnaces using Armco iron and standard alloying materials. The ingots were forged to bars, and sections were normalized prior to machining Jominy hardenability test specimens. The Jominy test was carried out under standard conditions.2 The depth of hardness penetration for determination of ideal diameter at "half hardness" was taken at the hardness level indicated in Fig. 20 in Grossmann's paper by the solid line for 50 per cent martensite in plain carbon steel. The depth of penetration of hardness on the Jominy specimen was converted to ideal diameter by the relation given in Fig. 28 of Grossmann's paper. Chemical analysis was determined on a sample taken adjacent to the Jominy specimen. Aluminum was determined as "acid-soluble aluminum" rather than as total aluminum, as it is probably more representative of the aluminum that is effective as an alloying agent. In addition to the elements nominally present in the steels, analyses were made for residual phosphorus, sulphur, nickel, copper, molybdenum and chromium on one heat of those made with the same materials during the same period. Actual grain size was determined by microscopic examination of small samples quenched in water after receiving the same thermal cycle as the Jominy hardenability specimens. Chemical analyses, heat-treatment, grain size, and hardenability are shown in Table I representing three ranges of composition.

Jan 1, 1943

By J. K. Stanley

It has been known for many years that the magnetic properties of a ferromagnetic material are very sensitive to internal strain. Any structure-sensitive property such as ferromagnetism, which is a function of the regularity of the atomic arrangement, and discontinuities such as grain boundaries must show deviations from an equilibrium state when the regularity of the lattice is disturbed. Experience has shown that such distortions have a marked effect on the magnetic properties of a ferromagnetic material. The purpose of this paper is to show the feasibility of using magnetic methods in detecting internal strains. More specifically, it will be shown how such properties as permeability, rema-nence, and coercive force change on cold-working of aluminum iron and how these magnetic properties change during the annealing below the recrystallization temperatures. This work was conducted for two reasons. One was the desire to see in what manner strains are relieved in aluminum iron at low temperatures—a subject of practical importance in soft magnetic materials—and the other was the feeling that a study of cold-working and strain relief might shed some light onour understandingof what takes place in the deformation of metals. Definition of Terms It is necessary to define the terminology used in this Paper because magnetic nomenclature is usually unfamiliar to the metallurgist, and certain other terms, such as recovery, have been used ambiguously by metallurgists and others. Magnetic terms are best explained by reference to the common hysteresis loop of ferromagnetic materials (Fig. I). In Fig. 1 the coordinates are magnetic induction, B, (ordinate) and magnetizing force, H. (abscissa). The c.g.s. unit for the magnetic induction is the gauss, the number of flux lines per square centimeter, and the unit for the magnetizing force is the oersted, the number of lines of force per centimeter in air or vacuum. If one starts with a completely demagnetized specimen and magnetizes it, a typical virgin magnetization curve beginning at the origin is obtained. This curve is also the locus of tips of a family of smaller hysteresis loops. From this curve it is to be noted that the induction is not proportional to the magnetizing force, and that it appears to approach a limiting induction called the saturation value, and does so for values of H of 1000 oersteds. The ratio of the magnetic induction (gauss) to the magnetizing force (oersteds) is called permeability, . Thus u= B/H and if the magnetizing force H (oersteds) is specified as ro or 100 oersteds, the permeability is written as or l00 This is the accepted notation, and will be used in this paper. The values of permeability for H = 10 or I00 are chosen as indices for 106

Jan 1, 1945

By J. E. Chappelear

We are concerned with the problem of calculating the forces produced by passing waves on a submerged object (in particular a drilling barge) located on the ocean floor. We employ the linearized theory of waves (the Airy1 theory), since solutions are generally more easily obtained in such a theory. It will be indicated how an approximate solution to the problem can be obtained for waves of finite height by use of a linear combination of Airy waves. We shall assume the submerged object has an infinite length, thus neglecting diffraction and making the problem two-dimensional. We find that there is not a considerable difference between our solution and the usual engineering solution. The force, as a function of time, is plotted in Figs. 4 and 5, together with the usual solution for easy comparison. So far, no experiments have been found to compare with the results. As will be shown, the calculations can be extended to waves of finite height (Stokes2-Struik3 waves) in an approximate fashion. This can be done by calculating the force produced by each Fourier component of the wave, separately, and omitting the coupling between the various harmonics produced by Bernoulli's theorem. The problem and its boundary conditions will be briefly stated. We seek a function, +, the velocity potential, which is harmonic, ?2?=0...........(1) and which satisfies the boundary conditions ??/?n =0 .......(2) (the normal derivative of + is zero) on the ocean floor and the sides of the barge, and s??/?y = ?..........(3) on the ocean surface, where s = g/w2........(4) The physical reasons for these choices and an interpretation of ? are given by Lamb.' As previously stated, we assume also that the problem is two-dimensional. We set up our coordinates such that x = 0 is the center of the barge and y = 0 is the ocean floor. The coordinates are shown in Fig. 1, together with the location of the barge on the bottom. Our solution may be decomposed into a function which is symmetric about x = 0,?e, and one which is antisymmetric, +,, ? = ? + ?o...........(5) Thus directly above the barge we have as boundary conditions that for x = 0, ?o = 0, ...........(6) and ??e/?x = 0........(7) Now, far from the barge, we have that (for x positive and assuming the waves are incident from the left) ?e ? Ae cosh ky cos k (x + a) , . . . (8) and ?e ? Ao cosh ky cos k (x + ß) . . . . (9) where + means ". .. approaches ... for x large." Let us put Ae, Ao, and k equal to one, and remember that, in accordance with the derivation of Eq. 3, a sinusoidal time factor is to be supplied. Then we may rewrite Eq. 8 for ?e, and we find that

Jan 1, 1958

By M. I. Sherman, J. D. H. Strickland

PRELIMINARY experiments in these laboratories showed that whereas pyrite1 produced only sul-fate the action of aqueous chlorine solutions on most other sulfide ores resulted in the formation of a mixture of sulfate and elemental sulfur, the production of the latter being favored by a very dilute solution of oxidant. The reactions taking place were obviously complex and only a detailed kinetic study could be expected to clarify the position and enable a prediction to be made of conditions leading to the rapid formation of metal salt and the economically desirable elemental form of sulfur. Galena was chosen for the first investigation as it is of widespread occurrence and is representative of a simple cubic lattice sulfide containing the monosulfide unit, S2-. Experimental Procedure The apparatus and general technique used in this work have already been described.' Analyses for chlorine and sulfate were carried out by the same methods as those used for pyrite. The lead in 1 ml aliquots of solution was determined polarographi-cally, after destroying chlorine with hydrazine hy-drochloride and making solutions 0.5 molar with hydrochloric acid. A large sample of Kansas galena, assaying at better than 98 pet PbS, was crushed and wet sieved to give fractions of —10 + 14, —20 + 28, —28 + 35, and —48 + 65 Tyler screen. The particles were almost perfect cubes and the apparent surface Per gram was calculated assuming each particle a cube Of side equal to the mean mesh opening. A volume shape factor, calculated as for pyrite, was exactly unity. The main apparatus was unsuitable for direct visual examination of the ore particles during the reaction. For this purpose a small cell from a Spekker Absorptiometer optical flats 1 cm apart) was l by a greased microscope slide after filling with oxidant solution. One or two ore particles were placed in the cell and the solution was stirred by a glass covered speck of iron wire activated by an external rotating magnet. The surface of an ore particle was examined under a binocular microscope at between 10 and 40 magnifications. Results Experiments were made under a wide variety of conditions, the form of presentation of the data being illustrated by Fig. 1, where the concentrations of chlorine, sulfate, and sulfur are plotted against time. The fictitious concentration of solid sulfur is used for convenience and is the molarity of sulfur which would be recorded if the elementary sulfur formed during the reaction were dissolved in the liquid iii the reaction vessel. As some lead sulfate is occluded with the sulfur on the ore particles the amounts of sulfur and sulfate formed during a reaction cannot be calculated

Jan 1, 1958

By B. F. Tillson

HAVING been asked to talk about recent devel-opments in mining equipment, I chose the sub-ject of small units of storage-battery locomo-tives which recently have been developed with the idea of motorizing transportation in the mines as far as possible, replacing the manual labor of hand-tramming. Metal mines are probably less advanced in this field than coal mines, where transportation of large ton-nages on any one level is quite common. Of course, metal mines, where large tonnages are handled on any one level, are as advanced as coal mines, and storage-battery locomotives for main haulages are used quite commonly in such metal mines, but the mines having levels at many intervals in the depth of the orebody, or a number of intermediate levels in which a consider-able amount of hand-tramming has to be done, and small metal mines in which the size of the orebody does not seem to warrant the initial capital expenditure for the driving of main haulage levels and installation of rock and ore chutes and raises for the delivery of the broken ore and rock to the cars have been at consider-able disadvantage in replacing manual labor by some mechanical means of transportation. Several years ago the field was canvassed among the manufacturers of electric-locomotive equipment and the conclusion was reached that the art had not been devel-oped to the advantage of metal miners who desire to replace small units of labor by locomotive equipment. The need for such small locomotives is all the more em-phasized by the development in the past few years of mechanical mucking-machines (loading equipment for filling the cars from muck piles on the level). Any me-chanical loader is a liability instead of an asset unless means are also provided for keeping the material away from the loading-machines. Therefore the transporta-tion system must be equal to the capacity of a loading-machine, otherwise the latter is idle such a large per-centage of its time as to prove an unprofitable invest-ment. So there is a need for portable, small haulage units to handle small mine cars in many metal mines, and, naturally, units are desired which can be so readily shifted from level to level as to facilitate their free use. Of course, in mines in which each level comes to the shaft it is not such a difficult problem, but many mines have a great many intermediate levels which must be approached through raises, timber slides or other small openings required to transfer ore and to get timber and other supplies to those levels.

Jan 5, 1927

By H. Becke, D. Stolnitz, D. Flatley, W. Kern

Zinc difhsions in gallium arsenide having surface concentrations as low as 5 x 10'' atoms per cu cm have been attained. A multiple-difhsion sequence is employed during which zinc enters the gallium arsenide from a silicon dioxide film in a solid-to-solid diffusion step. The technique is illustrated for the case of a zinc surface concentration of -1 x 1017 atoms per cu cm and a junction depth of 1.2 p. Radiotracer measurements of -0.1-p layers of gallium arsenide are used to delineate the resulting profile. A diffusion coefficient for zinc in pyrolytically deposited silicon dioxide is reported. ZINC is frequently used as an acceptor impurity in gallium arsenide. As examples, p-n junctions are formed using zinc-diffused layers in solar cells, lasers, and varactor diodes. For these purposes surface concentrations of 10'' to over loz0 zinc atoms per cu cm are desired and readily obtained from a vapor source. Diffusion of zinc into gallium arsenide to yield surface concentrations of about 1017 atoms per cu cm, however, is more difficult. For instance, oldstein' reports that large changes in the vapor density of zinc produce only small changes in surface concentration. Such low surface concentrations are necessary in forming the diffused base of a gallium arsenide double-diffused n-p-n transistor in order to obtain an adequate emitter efficiency. An alternative approach to zinc diffusion which can provide reduction in surface concentration of several orders of magnitude is the introduction of the diffusant from a solid source. This can be conveniently accomplished by using oxide films of the type generally employed for selective diffusion masking. Such oxides also serve to prevent the dissociation of the gallium arsenide surface at the high temperatures of diffusion. DIFFUSION PROCEDURE Fig. 1 illustrates the diffusion sequence which is employed for obtaining zinc surface concentrations in the range from 5 x lom to 1 X 10" atoms per cu cm. The description of technique will be given for the instance of -1 x 1017, which is typically utilized in double-diffused gallium arsenide transistor investigations at this laboratory. Tin- or germanium-doped gallium arsenide wafers drawn from crystals grown by horizontal Bridgman or Czochralski techniques are used, with bulk concentrations ranging from 1 to 3 X 10" carrier atoms per cu cm. Each wafer is chemically polished on the (iii) arsenic face using a mixture of five parts concentrated sul-furic acid, one part water, and one part 30 pct hydrogen peroxide at room temperature, followed by rinsing in deionized and quartz-distilled water. Pyrolytic silicon dioxide was chosen as the oxide coating because of the high quality of the films obtainable, the ease of their preparation, and the excellent control over film thickness and uniformity.

Jan 1, 1964

By W. C. Ellis, M. E. Fine

In the iron-cobalt system there are several property-composition relationships of theoretical importance. The alloys are ferromagnetic exhibiting a maximum saturation at approximately 33 at. pct cobalt;' the interatomic distances at room temperature show large positive deviations from Vegard's rule;2 and in the region of 50 at. pct cobalt there is an order-disorder transition2 which has a large effect on other properties. The thermal expansion data for the system have been reported only from room temperature to 100 C and for this temperature range a minimum in coefticient in the rcgion of 50 at. pct cobalt is well established.3,4,6 This paper reports the expansion properties in the iron-cobalt system for temperatures from 30 to 850°C, shows the correlation with other physical properties, and proposes an explanation for the tcrrelation. Experimental Procedure The thermal expansivities of iron-cobalt alloys in 10 pct steps in composition were determined from 30 to about 850°C. The major series of alloys was prepared by induction melting from Armco iron and cobalt rondelles with the addition of about 0.6 wt pct manganese. This series was supplemented by three compositions: high purity electrolytic iron, high purity electrolytic cobalt, and an alloy of 52.1 wt pct cobalt prepared from the high purity electrolytic iron and electrolytic cobalt by vacuum-melting. The compositions of both series of alloys are given in Table I. The specimens for measurement were prepared from the cast ingots by first hot-swaging followed in most instances by cold-swaging to 0.180-in. diam rod. Some alloys in the region of 5c pct cobalt, because of inherent cold-brittleness, were swaged hot to the final diameter. The expansion specimens were approximately 2 1/2 in. long with flat-ground and polished ends. Before measurement they were annealed in hydrogen for I hr at 900°C. In the apparatus used for the thermal expansivity measurements (Fig I), a modification of that designed by Hidnert and Sweeney,? the differential expausion between the specimen and vitreous silica was measured by a dial micrometer reading to 0.001 mm. An invar fixture supported the dial gauge, and this fixture, in turn, was connected to a split invar collar which supported the outer silica tube. A pyrex "T" was placed over the outcr silica tube and was sealed to the invar collar by Apiezon Q sealing compound. A movable hollow vitreous silica rod which rested on top of the test specimen slid through a flexible rubber diaphragm cemented to the top of the pyrex "T." The specimen rested On a vitreous silica pedestal fused into the outer tube and was kept centered by a silica sleeve. A nonoxidizing atmosphere was obtained by admitting hydrogen into the

Jan 1, 1949

By William Hurst

This is a presentation of the diffusivity theory for the calculation of the water drive on an oil reservoir in which the history of reservoir pressure with time are the essential parameters for the determination of the rate and cumulative water encroachment into a field. The main body of the paper is essentially the application of the principles underlined here to actual field data. It consists of the work plots of the radial-flow case, and includes a discussion of an illustrated problem for which the water-drive calculations are treated in conjunction with the volumetric-balance equation. The appendix is a mathematical treatment of the linear, radial, and spherical cases of water influx into a reservoir. Introduction The evaluation of oil reserves and the prediction of pressure behavior with different rates of production are related to the volumetric balance of the fluids entering and leaving the space occupied by the oil and gas originally present in the formation; that is to say, important information can be gleaned from a volumetric study of the oil, gas, and water produced from a formation as it affects the extent of change in volume of oil zone, the gas cap, and the amount of water influx from the adjoining water sands. It is not the purpose of this paper to enter into a discussion of the equation of volumetric balance, in which such pertinent information as subsurface sampling and the measurement of reservoir pressures play a predominant part, as this has been discussed in detail by others.1,2 However, there is an inherent difficulty that has limited the general application of this equation to a study of oil reservoirs; that is, the insufficient knowledge of the analytics for expressing the water drive on a field as a function of reservoir pressure. This to a large degree has been overcome and it is the purpose of this paper to present this theory in such form as to be applicable for reservoir calculations. Water Drive If the amount of oil originally in place and the gas space are known, it is a relatively simple matter to determine, by the volumetric equation, the cumulative water influx into a reservoir from the adjoining water sands. Often, however, the oil reserve as well as water influx must be 'determined, and in order to be able to make this calculation it is essential to assign a mathematical relationship for the varying water term that is compatible with conditions known to exist in the formation and gives an accuracy acceptable for engineering calculations. Various equations have been employed empirically by different investigators in making reservoir studies of water drive. One form used by Schilthuis2 expresses the rate of water influx into a reservoir at any time as proportional to the pressure

Jan 1, 1943

By J. D. Rice, H. H. Rachford, G. H. Bruce, D. W. Peaceman

The problem of unsteady-state gas flow through porous media leads to a second-order non-linear partial differential equation for which no analytical solution has been found. In this paper a stable numerical procedure is developed for solving the equation for production of gas at constant rate from linear and radial systems. An electronic digital computer is used to perform the numerical integration using an implicit form of an approximating difference equation. Solutions are presented in graphical form for various values of dimension-less parameters. The solutions are compared with the laboratory study of gas depletion in a linear system. INTRODUCTION Production of fluids from porous rock reservoirs is essentially a transient process. Transient gradients develop as soon as production begins, and further withdrawals continue to cause disturbances which propagate throughout the reservoir, each adding in some way to the prior ones. A correct mathematical analysis of this behavior is complicated by the fact that the transient or unsteady-state flow of compressible fluids must be described by difficult second-order partial differential equations. As a practical matter, three distinctly different cases arise: 1. Flow of single-phase liquid 2. Flow of gases 3. Multiphase flow The first of these has been found to give a linear second-order equation similar to the well-known heat flow equation. where 7 is fluid density, 8 is time, and a a constant for the system, provided d 7 = r cdp, where the compressibility, c, is constant over the range of pressure, p, considered. Solutions of Equation (1) for both linear and radial flow are available in several forms.2'3' On the other hand. the second case, which is the flow of gas, gives a non-linear second-order equation, the solution of which is not known. Although a number of approximate solutions have been proposed, each is limited in value by the associated simplifying assumptions. Inasmuch as the analysis of transient flow is limited to liquid systems, a solution of the second case is necessary if further progress is to be made in studying underground fluid movement. For this reason a solution of Equation (2) was undertaken by means of numerical integration of approximating difference equations. BASIC DIFFERENTIAL EQUATION Equation (2; is derived by combining the equation of continuity, the perfect gas law, and Darcy's law to give the basic differential equation Strictly speaking, the fluids present in a gas reservoir are not perfect gases. Furthermore, there are indications that over certain velocity ranges, Darcy's law is not applicable. It would not be practicable. however, to attempt to obtain numerical solutions which would take into account all the possible variations from these laws which could occur in actual systems. The

Jan 1, 1953

By P. A. Beck, A. Merlini

Quantitative texture determinations were made for rolled and for rolled and annealed strips of 3, 6, and 10 pct Zn brasses. The main components of the rolling texture gradually shift over a wide composition range from the orientation characteristic of copper to that of 70-30 brass. At least the main features of the corresponding gradual changes in the annealing texture may be accounted for qualitatively on the basis of the oriented growth theory of annealing textures. All the oriented nucleation-type theories are notably less successful in accounting for the observed textures. EARLIER investigations by von Goeler and Sachs,' Dahl and Pawlek,' and Brick, Martin, and Angier" suggested that the annealing texture of copper is more sensitive to the addition of solute elements than the rolling texture. It was reported3 that the cube texture changes over into a brass-type recrys-tallization texture on addition of about 1 pct Zn, while the deformation texture changes only at about 5 pct Zn. Such behavior would be difficult to reconcile with any theory of annealing texture formation that relates the annealing texture to the deformation texture. In particular, the oriented growth theory of the cube texture" is based on a certain orientation relationship between the cube texture and the main components of the deformation texture which could not be reconciled with the reported behavior of low brasses. It was, therefore, of interest to re-examine these textures by means of the quantitative methods now available. The deformation and annealing textures of rolled a brasses with approximately 3, 6, and 10 pct Zn have been determined in terms of (111) and (200) pole figures. Experimental Methods The chemical analyses of the Cu-Zn alloys of commercial purity used in the present work are given in Table I. Chill-cast bars of the tabulated compositions were rolled and annealed according to the following schedule: I—rolling from 1.2 to 1.02 in. thickness, 15 pct reduction area, annealing 15 min at 450°C; 2— rolling from 1.02 to 0.867 in. thickness, 15 pct reduction area, annealing 15 min at 450°C; 3—rolling from 0.867 to 0.600 in. thickness, 31 pct reduction area, annealing 2 hr at 540°C; 4—rolling from 0.600 to 0.400 in. thickness, 33 pct reduction area, annealing 2 hr at 480°C. In the final straight rolling of 96.25 pct reduction area, the strips were reversed end to end after each pass. The final annealing temperature was determined for each alloy by annealing specimens of the cold-rolled strips in salt bath for 5 min at varying temperatures. The results of the hardness measurements are given in Fig. 1, where the final annealing temperature used for each alloy in determining the annealing texture is indicated by a cross.

Jan 1, 1956

By William Hurst

This is a presentation of the diffusivity theory for the calculation of the water drive on an oil reservoir in which the history of reservoir pressure with time are the essential parameters for the determination of the rate and cumulative water encroachment into a field. The main body of the paper is essentially the application of the principles underlined here to actual field data. It consists of the work plots of the radial-flow case, and includes a discussion of an illustrated problem for which the water-drive calculations are treated in conjunction with the volumetric-balance equation. The appendix is a mathematical treatment of the linear, radial, and spherical cases of water influx into a reservoir. Introduction The evaluation of oil reserves and the prediction of pressure behavior with different rates of production are related to the volumetric balance of the fluids entering and leaving the space occupied by the oil and gas originally present in the formation; that is to say, important information can be gleaned from a volumetric study of the oil, gas, and water produced from a formation as it affects the extent of change in volume of oil zone, the gas cap, and the amount of water influx from the adjoining water sands. It is not the purpose of this paper to enter into a discussion of the equation of volumetric balance, in which such pertinent information as subsurface sampling and the measurement of reservoir pressures play a predominant part, as this has been discussed in detail by others.1,2 However, there is an inherent difficulty that has limited the general application of this equation to a study of oil reservoirs; that is, the insufficient knowledge of the analytics for expressing the water drive on a field as a function of reservoir pressure. This to a large degree has been overcome and it is the purpose of this paper to present this theory in such form as to be applicable for reservoir calculations. Water Drive If the amount of oil originally in place and the gas space are known, it is a relatively simple matter to determine, by the volumetric equation, the cumulative water influx into a reservoir from the adjoining water sands. Often, however, the oil reserve as well as water influx must be 'determined, and in order to be able to make this calculation it is essential to assign a mathematical relationship for the varying water term that is compatible with conditions known to exist in the formation and gives an accuracy acceptable for engineering calculations. Various equations have been employed empirically by different investigators in making reservoir studies of water drive. One form used by Schilthuis2 expresses the rate of water influx into a reservoir at any time as proportional to the pressure

Jan 1, 1943

By J. J. Arps, T. G. Roberts

Since the introduction of the relative permeability concept in the middle thirties1,2 various investigators have shown3,10,11,12,15 how the basic equations for the flow of oil and gas through porous media can he utilized to compute the recovery from depletion type or solution gas drive type reservoirs under certain conditions, when the necessary physical data pertaining to the reservoir rock and its reservoir fluids are known. In the study by Muskat and Taylor" in 1945, the effects of viscosity, gas solubility, shrinkage, and the permeability-saturation characteristics of the producing formation on the production histories and the recovery of gas drive reservoirs were analyzed. Each parameter was varied a limited number of times, while keeping the others constant. The work was not carried to the point where the data could be used directly to estimate the primary recovery in those cases where certain rock and fluid characteristics were known or could be assumed. Since the publication of this classic study, the work by Beal in 194616 and Standing in 1952'" has established general correlations of oil viscosity and shrinkage with gas solubility and oil gravity. In addition. a considerable number of relative permeability relationships for different types of reservoir rocks has been published2,5,13,17,19,20. Also, the burden of the numerical work required in the solutions of the differential equation has. during the last few years, been greatly reduced with the availability of modern electronic computing facilities. The purpose of this paper is therefore to present the results of a large number of such recovery computations covering the normal range of variation of the main parameters such as the type of rock, the API gravity of the oil, and the amount of gas in solution. As is to be expected, the ultimate recovery is found to increase with the oil gravity. except for the higher solution gas-oil ratios. Within the normal range of solubilities, the ultimate recovery generally appears to decrease with the amount of gas in solution for the higher oil gravities, and to show a slight increase for the lower oil gravities. Intermediate oil gravities seem to fall between the two trends. The type of reservoir rock, as identified by its relative permeability relationship, appears to have a very pronounced effect on the ultimate recovery. In general, sands and sandstones show higher recoveries than limestones and dolomites, although certain intergranular limestones seem to show a higher theoretical recovery than unconsolidated sand. The recovery from sands and sandstones generally decreases with increasing cementation and consolidation, while the recoveries from limestones and dolomites are highest for the intergranular type and lowest for the fracture type porosity.

Jan 1, 1956

By John Chipman, Raymond L. Orr

Recently reported high -temperatlnre thermal data have been incorporated in a reuieu' of the therinody-r~ar,zic propesties of iron. The recent data permit mot corlsisterlt and reliuble choices to be made for many of the pvoperties, resnlling in bettev eslablished tabrtlations of the thermodynamic functions. Values of Cp, h°t- - H°298. s°t - S2sa, and (G°T -H298/T are tabulated for -the solid and liquid phases of iron in both stable and melastable regions. Consistency of the tabulations has been maintained to the precision necessr~vg to yield the Gibbs energy change belween the bcc (a, 8) and fcc (y) phases of iron at temperatures pertinenl to alloy studies. THE importance of the metal iron has stimulated many investigations of its thermodynamic properties at elevated temperatures. Reviews and summaries of the data have appeared from time to time. Among these the 1951 paper of Darken and smith' deserves highest recognition for the quality of its criticism and for the precision and thermodynamic consistency of its tabulations. These and the more recent table of Kaufman, Clougherty, and weiss2 reproduce the Gibbs energy of the (a, 6)--Y transformation with sufficient sensitivity for alloy studies. Other recent compilations of equal validity3-' make no attempt to do this. Since this Gibbs energy change is not greater than 22 cal per g-atom over the temperature range 1100° to 1800°K, it is obviously desirable to tabulate values of (G°t - H298)/T to four places of decimals even though the absolute accuracy would not warrant such precision. The necessity for reconsideration of the problem is emphasized by the recent publication of several important researches on the enthalpy and heat capacity of iron. The older work is well-covered in the reviews cited and will be mentioned here only when actually employed in the tabulations. GAMMA IRON Transition temperatures are taken from the selec -tions of Elliott and Gleiser,5 based on the work of Boulanger7 and earlier workers. For the a-y (A3) transition at 1184°K we retain the value 215 cal per g-atom adopted by Darken and smith.' It has received confirmation in the recent work of Dench and Kubaschewski,' and of Braun8 using adiabatic calori-metry, and of Wallace, Sidles, and Danielson10 using a pulse heating technique. The enthalpy data of Olette and ferrier" and of Anderson and Hultgren12 at higher and lower temperatures may be extended to 1184°K to give the following values which are adjusted to the required value of ? 1184 Ha,1184 Ha 298 = 8030 cal per g-atom [11 hr,1184-ha298 = 8245 cal per g-atom [2]

Jan 1, 1968

By G. H. Alani, H. T. Kennedy

One of the major difficulties in predicting the performance of oil reservoirs from their early pressure history lies in the uncertainty of estimating the volume of the liquid hydrocarbons contained in them. As a first step in filling this need, an equation was developed to determine the Inolal liquid volume of pure hydrocarbons over a wide range of temperature and pressure. The second step consisted of adapting the equation to apply to mixtures, with the heary hydrocarbons expressed as c The equatious are similar in form to van der Waals equation. hut the constants a and b tire considered as furnctions of temperature. In addition to the gas constant R, there are four constants characteristic of each hydro- When comptired with experimental values jound in the literal the. the average absolute deviation in the calculated molal volumes is found to be a Maximum of 0.33 per -cem for any of the pure liquid hydrocarbons studied. This maximum deviation was that found when comparing the calculated and observed values over a temperature range of 86° to 482°F and a pressure range from the bubbl-point ro 10,000 psis. The equations expressing the correlation for. mixrurrc were developed from 647 experimental measure-rilents of volume on 47 bouom-hole samples covering ; temperature range of 72" to 250°F and a pressure znge from bubble-point ro 5,000 psig. The average absotute deviation was found to he 1.6 per cent with the maximum for any measurement of 4.9 per cent. INTRODUCTION Accuratc information of the pressure-volume-temper-iiture behavior of hydrocarbon liquids is of considerable importance in the field of both applied and theoretical science and, especially, in the solution of petroleum reservoir engineering problems. These PVT relationships can be expressed graphically, in tabular form or as equations of state. Equations of state for gases have been developed based on the ratio of pressure, temperature and volume to the pressure, temperature and volume at the critical points of the hydrocarbons. These equations have been successful in yielding values fairly accurate for reservoir engineering purposes. For liquids, however, relations of this type are more complex and have not yielded useful results. The composition of reservoir liquids is normally expressed in terms of the moles of hydrocarbons C,? C,. C,. elc, per mole of mixture. The composition must include the fraction of a mixture such as C,,, which cannot be separated with reasonable time and effort into specific hydrocarbons. Although there are several equations of state having variable degrees of accuracy applicable to pure liquid hydrocarbons, none includes a means for estimating the volume of mixtures as a function of temperature. pressure and composition expressed in this manner.

Since production curtailment for other than engineering reasons is gradually disappearing, and more and more wells are now producing at capacity and showing declining production rates, it was considered timely to present a brief review of the development of decline-curve analysis during the past three or four decades. Several of the commoner types of decline curves were discussed in detail and the mathematical relationships between production rate, time, cumulative production and decline percentage for each case were studied. The well-known loss-ratio method was found to be an extremely valuable tool for statistical analysis and extrapolation of various types of curves. A tentative classification of decline curves, based on their loss ratios, was suggested. Some new graphical methods were introduced to facilitate estimation of the future life and the future production of producing properties where curves are plotted on semilogarithmic paper. To facilitate graphical extrapolation of hyperbolic-type decline curves, a series of decline charts was proposed, which will make straight-line extrapolation of both rate-time and rate-cumulative curves possible. Introduction During the period of severe production curtailment, which is now behind us, production-decline curves lost most of their usefulness and popularity in prorated areas because the production rates of all wells, except those in the stripper class, were constant or almost constant. While production-decline curves were thus losing in importance for estimating reserves, an increasing reservoir consciousness and a better understanding of reservoir performance developed among petroleum engineers. This fact, together with intelligent interpretation and use of electric logs, core-analysis data, bottom-hole pressure behavior and physical characteristics of reservoir fluids, eliminated a considerable part of the guesswork in previous volumetric methods and put reserve estimates, based on this method, on a sound scientific basis' At the same time, a number of ingenious substitutes were developed for the regular production-decline curve, which made it possible to obtain an independent check on volumetric estimates in appraisal work, even though the production rates were constant. with the now steadily increasing demand for oil to supply the huge requirements of this global war, proration for reasons of than prevention of underground waste is gradually disappearing. More and more wells are, or will be, producing at capacity or at their optimum rates, as determined by sound engineering practice. With this trend, the character of producing wells seems to regain, more or less, its "individuality," and the old and familiar decline curve appears to have had a comeback as a valuable tool in the hands of the petroleum engineer. It may be timely, therefore, to retrace the development of decline-curve analysis in the past by presenting a brief chronological review of bulletins and papers published during the past three or four decades, which have

Jan 1, 1945

By T. H. Alden, H. W. Schadler

The strain rate dependence of the flow stress of the eutectoid Zn-Al alloy has been determined as a function of mechanical processing, microstructure, and temperature. The best superplastic properties result from the ultrafine, equiaxed structure produced by solution treatment and isothermal transformation at 0°C. Lamellar structures produced by isothermal transformation and coarse equiaxed structures produced by annealing also exhibit strain rate sensitivity in excess of 0.5, but at higher stresses and lower strain rates than the finest structure. Solution treatment and isothermal transformation eliminate the effect of prior processing history. Small composition variations in the binary alloy do not influence the mechanical properties. Below the eutectoid transformation temperature (.275C) increasing temperature lowers the flow stress at all strain rates. The data do not favor any of the present proposals for the mechanism of superplastic flow. ALLOYS of zinc and aluminum, near the eutectoid composition, have been investigated for: 1) the flow stress-strain rate and elongation behavior,1, 2) the transformation kinetics,3 and 3) structure,3"5 but the structure sensitivity of the mechanical properties is poorly understood. Using a near eutectoid alloy water-quenched from 380°C, Backofen et aL2 estab- high and increases with increasing strain rate. They found a correlation between strain rate sensitivity and tensile elongation as expected from the increased resistance to necking that a high strain rate sensitivity provides.296 The large elongation at high strain rates is attractive for metal forming. Since most other superplastic alloys7-l3 do not show high strain rate sensitivity at such a high strain rate (for the exception see Ref. 12) there is the possibility of a unique deformation mechanism for this alloy. The sensitivity of the mechanical properties of the eutectoid alloy to thermal history1 suggests a structure sensitivity that has not been demonstrated. Gar-wood and Hopkins3 identified a lamellar product for material transformed between 200" and 260°C. The size of the lamellae decreased with decreasing temperature. They also demonstrated that an unresolved transformation product formed below 200°C. A granular two-phase material was resolved for specimens quenched to below 150°C by Mitbauer and sauerwald4 who observed coarsening of the structure after annealing. The role of grain boundaries in superplastic deformation is suggested by the strong effect of grain refinement in encouraging this phenomenon at high strain rateS7-11, 13 and is emphasized in several proposals about rnechanism7, 8, 12 and in a theory14 of high-temperature deformation. Alden8 showed metallographic evidence that grain boundary sliding is occurring extensively at strain rates where m is large but less at high strain rates where normal plasticity is observed. Holt and Backofen10 have shown that sliding also oc-

Jan 1, 1969

By R. D. Pehke, J. F. Elliott

The solubility of nitrogen in liquid pure iron has been measured as a function of pressure and temperature. Sieverts' Law is obeyed at all pressures up to 1 atm and the temperature coefficient of solubility is 8 x 106 %/°C. The solubility of nitrogen in liquid iron alloys has been studied by the Sieverts' method and by the sampling method. The solubility is decreased by Al, C, Co Cu, Ni, 0. Si atzd Sn, and it is increased by Cr, Cb, Mn, Mo, Ta, W, and V. Interaction coefficients eN, for the effect of these elements on nitrogen dissolved in liquid iron, are reported. 1 HE effect of nitrogen upon the properties of steel may be desirable or not, depending upon the composition, processing treatment, and the use of the product. To understand the behavior of nitrogen in the various steelmaking processes, it is necessary to have reliable fundamental data pertaining to the solubilities, rates of solution, and chemical reactions of nitrogen in steel, particularly in the liquid state. In an effort to extend our knowledge in these areas, an investigation of the thermodynamics and kinetics of the solubility of nitrogen in liquid iron alloys was undertaken. The thermodynamic aspects of the study are presented in this paper. The solubility of nitrogen in liquid pure iron has received the attention of several investigators whose results are summarized in Table I. Reasonable agreement among the more recent researches1-l8 is found for the solubility of approximately 0.044 wt pct at 1600°C and 1-atm pressure of N,, for the temperature coefficient of solubility being small but positive, and for adherence to Sieverts' law for pressures up to 1 atm. The solubility of nitrogen in a number of binary liquid iron alloys has been studied: Fe-A1,3 Fe-As,16 Fe-C,3, 71 18, 19, 20 Fe-Cr ,5,7, 10, 11, 15 21 Fe-Co,'4, '6 Fe-Cu,'6 Fe-Mn,10, 11, 18 Fe-Mo,13, 16 Fe-Ni,10, 11, 13, 15, 16 Fe-O,16, 17, l9 Fe-p,7 Fe-S,l6 Fe-Sb,16 Fe-Se,lg Fe-Si,4) 97 '77 187 19>2° Fe-Sn,16 and Fe-v.51 '3 A few ternary systems have been studied, most of which were reviewed by ~angenberg" who applied the methods suggested by WagnerPz3 Chipman,24 Darken and G~rry, and Morris and Buehl26 for predicting solubilities in multicomponent alloys. Langenberg22 and Humbert and Elliott15 compared the solubilities computed from simple binary alloy systems with experimental measurements of nitrogen solubilities in liquid iron containing two or more alloying elements. The interaction coefficients for nitrogen in

Jan 1, 1961