Search Documents

By Howard Martens

HARDNESS behavior of commercial titanium alloys following various heat treating processes has been studied for some time. However, the hardness of such alloys following a definite measured cooling rate from the single phase region has not been reported. Therefore, a series of experiments was conducted on the three alloys: 150A, 130A, and 130B. The manufacturers gave the following compositions for these alloys: Ti 150A: 0.045 pct C, 0.069 pctN, 1.32 pct Fe, and 2.68 pct Cr; RC 130A: 0.13 pct C, and 7.9 pct Mn: RC 130B: 0.11 pct C, 3.5 pct Mn, and 3.2 pct Al. A method similar to that described by Greninger was used, and the cooling rate varied from 6" to 9200°F per sec. The cooling curves were recorded on a Speedomax or a recording oscillograph, and the hardness was measured on a Tukon hardness tester using a Knoop indentor with a 500 g load. The curves of Fig. 1 show the variation of the Knoop hardness with the cooling rate for these three alloys. The alloys 150A and 130B behave in a similar manner, as shown by the curves. In neither of these alloys was it; possible to retain the high temperature ß phases, even by the most rapid cooling rate. The slow cooling rates produced a typical ß structure in these two alloys. As the cooling was increased, these phalses became more finely dispersed and the hardness increased. This increase in hardness with cooling rate continued for the 150A alloy until a rate of approximately 1000°F per sec was reached. At this cooling rate, the structure showed the result of the dif-fusionless transformation of the ß solid solution to the supersaturated a solid solution which is referred as a'. At higher cooling rates, the structure showed no change and the hardness showed no great change, as indicated by the dotted portion of the curve in Fig. 1. The increase in hardness with cooling rate for the 130B alloy continued until a rate of approximately 300°F per sec was reached. At this rate the structure was produced and the maximum hardness was reached. Higher cooling rates caused no marked changes in hardness or structure. The alloy 130A behaved in a different manner, as shown by the curve of Fig. 1. The structure produced by the slow cooling rates was a typical a+ß structure, which became finely dispersed as the cooling rate increased. The hardness increased with cooling rate and reached a maximum at approximately 120°F per sec. At this maximum hardness, the structure showed the result of the decomposition of the ß phase into the ß + phase structure. As the cooling rate was increased, smaller amounts of the at phase were formed and the alloy became softer. For cooling rates of 3000°F per sec and above, the high temperature ß phase was retained and the hardness of the alloy did not change markedly. Acknowledgment This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. DA-04-495-0rd 18, sponsored by the Dept. of the Army, Ordnance Corps.

Jan 1, 1957

By M. T. Simnad, C. E. Birchenall, M. H. Davis

The transport mechanisms previously proposed in wustite and hematite have been confirmed by movement of inert markers during the growth of these phases. By similar measurements the mechanism for magnetite has been determined. The rates of growth of multilayer scales on iron, magnetite and hematite on wustite, and hematite on magnetite have been measured. All rates are found to be diffusion controlled under the conditions of the experiments. MANY studies have been made of the rate of oxidation of iron in both the low temperature, film-forming region and the high temperature scaling range. The interpretation of the results has been greatly complicated by the formation of multi-layered scales containing the various phases of iron oxides. In most cases no attempt has been made to determine the composition of the scale and its dependence on time of oxidation in a quantitative manner. However, several recent studies'" have taken this variable into account. It is generally agreed that the growth of thick scales on iron follows a parabolic rate law for the increase of total thickness. In the thin film region agreement has been claimed with parabolic equations of various forms and with logarithmic equations. The discussion here will be concerned only with the formation of scales of considerable thickness, greater than several microns, although some of the conclusions may be applicable to thin films. A brief survey will be given of the solid phases of the iron-oxygen system, the stability of wiistite in bulk and in thin films, the relative densities of iron in the metal and oxides, and the kinetics of forma- tion of oxide scales on iron. New experimental results will be introduced at several points during the discussion, and an attempt will be made to formulate a mechanistic picture consistent with these observations. The experimental techniques are described in an appendix. Iron-Oxygen System Fig. la shows the iron-oxygen phase diagram. Wüstite ("FeO"): The lowest of the oxides has a NaCl type cubic lattice formed by close-packing of the relatively large oxide ions with the smaller iron ions arranged in the interstices. The ionic radius for 0' is about 1.3 A, while Fe++ has a radius of about 0.8A and Fe"' about 0.7A. Its range of existence on the phase diagram is very wide for an ionic material of this type and does not seem to include the stoi-chiometric composition FeO. These deviations appear to be well established as vacancies on iron lattice sites by the comparison of lattice parameters and densities for a series of compositions by Jette and Foote.4 The limits of the field are disputed, but the work which seems most reliable" indicates vacancies in from 5 pct to more than 12 pct of the iron sites. In order for the crystal to be electrically neutral, there must be electron holes which correspond to the formation of trivalent iron ions equal to twice the number of vacancies. Magnetite ("Fe³O4"): Considerable disagreement also exists in the literature about the stable range of the Fe³O4 field. Some diagrams show the Fe³O4 composition to be lower in oxygen than the formula indicates, but Schenck6 gives a diagram from the

Jan 1, 1952

By R. C. Edwards, Will Mitchell, T. G. Kirkland

THREE years ago, when the Process Research Laboratory at Allis-Chalmers Manufacturing Co. sought a remedy for the increasing cost of disposing of great quantities of spent sands from foundries, R. C. Edwards developed a continuous inclined paddle scrubber which removed a substantial amount of the charred coatings from the grains with a surprisingly small amount of degradation of the sand itself.1 The satisfactory results obtained on foundry sand, using this piece of equipment in conjunction with classification, suggested that a similar procedure could probably be applied to the cleaning of other sands, such as those required by the glass and chemical industries. This paper describes the flowsheet and the results of tests performed on impure natural sands. The four materials tested were silica sands obtained from widely separated geographical areas. Three of the sands tested contained various amounts of iron oxide and the fourth an excessive amount of calcareous bonding material. The iron oxide was present predominantly in the form of grain coatings, which were scuffed by the paddle scrubber, and in all cases a product was manufactured which would meet one of the sets of specifications imposed by the glass industry' The Edwards paddle scrubber used in these tests consists of a rubber-lined steel cylinder 18 in. in diam and 48 in. long, mounted on a frame so that its slope can be varied 5 ° to 35º from horizontal. This variable slope permits controlled flow of material through the machine. Rubber paddles rotate inside the cylinder and are secured to a center shaft 2 in. square, which extends through the entire length and is turned at both ends so that it can be supported by sleeve bearings centrally fixed in the end plates of the shell. At the upper or feed end it extends through the bearing on the end plate into a sheave, to which the motor drive is connected by means of Texrope transmission. The four sets of paddles mounted on the shaft, two in each set, are made of rubber belting, 1/4 in. thick, 9 in. wide, and 16 in. long. The two paddles in each set are secured on opposite sides of the square shaft, but in the same longitudinal position along the shaft axis. In other words, the two paddles in the set operate 180° apart. Each successive set is mounted 90º from the preceding one. One end of each paddle is bolted to the square shaft and the other weighted by a steel bar, 1x1 ½ x9 in. long, which is attached about 1/ 2 in. from the outer tip of the paddle. By centrifugal force at proper shaft speed, these weights force the paddle surface to rub smoothly against the rubber lining of the cylinder. During operation, the sand forces its way between the paddle face and the surface of the cylinder, providing the necessary abrading action as the pulp flows by gravity down the length of the cylinder. The capacity of the machine is limited by the amount and the hardness of the coating material and by the quality of product desired. The degree of scrubbing is governed by the retention time in the scrubber, which in turn may be changed by varying the feed rate, water content of the feed, and

Jan 1, 1952

By R. C. Gifkins

R. C. Gifkins (CSIRO)—In this paper evidence is put forward to support the idea of grain boundary shearing in aluminum at 4.2°K and the phenomenon is explained in terms of a low-temperature "equicohesive temperature". It is stated that these findings show that grain boundary shear cannot be ruled out in any situation by temperature arguments alone. We believe, from our own recent study of "apparent grain boundary sliding'' at temperatures below 0.45 Tm, that it can be ruled out on temperature arguments alone.45 Our experiments were made principally with magnesium and a Ms-A1 alloy ("Magnox A12"). Careful examination of offsets of marker lines or steps at grain boundaries produced by deformation at room temperature (-0.3 Tm) showed many of these to result from curvature over narrow zones along the grain boundaries. A significant residuum did, however, appear as a sharp offset or step. It was found that the mean value of this step rose to a small value during the first 0.5 to 1 pct strain and then remained at this level with further strain. Moreover, exactly the same values were obtained (at the same stress) at 75°C instead of 21°C. That is, in this range of temperature the "apparent sliding" was not temperature-sensitive. At 125°C and above there was sliding which increased with strain and which was temperature-dependent. When this temperature-dependent sliding was used to calculate the sliding to be expected at room temperature, this underestimated the "apparent sliding" by a factor of 104. The conclusion drawn is that the "apparent sliding" consists of shear localized in a zone of a width measured perhaps in microns, whereas it seems that true sliding occurs on the boundary surface.46 A model for such localized shear has been developed in general terms, following a similar model of Urie and Wain.47 Other strong evidence for low-temper- ature sliding has been based on the offsets of marker lines consisting of a substructure seen using the electron microscope to resolve these.48 We have repeated this work and found that the effects can be accounted for in terms of cracks in an oxide film which is produced when the substructure is etched up; these cracks form over zones of shear rather than at actual interfaces. However, in aluminum at room temperature (-0.3 Tm) there is no "apparent sliding" which is optically detectable. We suggest that, when the temperature is much lower (and possibly the rate of strain higher, too), as in the experiments of Chin et al., then aluminum behaves like magnesium at room temperature, and a narrow zone of shear develops to accommodate differences in strain across the grain boundary. Such a zone could well give rise to triple-point cracks. It is hoped to investigate these suggestions further. B. Y. Chin, W. F. Hosford, Jr., and W. A. Backofen (authors' reply)—It is interesting to have Dr. Gif-kin's comments, although it is somewhat difficult to reply in detail without the reference containing the work that he describes. Apparently though, he does find a low-temperature temperature-independent sliding that would seem to be related to the observations reported here. Perhaps some of the difficulty has semantic origins, growing out of the distinction between ('true sliding" and "apparent sliding".

Jan 1, 1965

By R. C. Troxell, S. C. Sun

EFFORTS to obtain definite quantitative data when employing the currently used bubble pick-up method1,2 as a pre-flotation investigation tool led to the adoption of the magnifying mirror method and the microscopic method modifications. Quantitative data obtained by these two methods were converted into bubble pick-up indices and bubble pick-up coefficients. Rough correlations were established between the calculated indices and the results of contact angle testing as well as actual flotation under the same conditions. The basic steps of the magnifying mirror method are as follows: 1—mixing 200 ml distilled water, sized and cleaned mineral particles of approximately 0.5 g, pH regulator, and flotation reagents in a glass beaker; 2—swirling the mineral particles to the center of the beaker to leave a clear circumference; 3—pressing an air bubble against the mineral particles for about 1 sec, and then moving the mineralized air bubble to the clear part of the beaker; and 4—dropping all the attached mineral particles of the air bubble by tapping the bubble holder, and counting the number of the dropped particles by looking at the magnified image in the mirror positioned beneath the beaker. This method is impractical for mineral particles which are transparent in water or finer than 200 mesh. The microscopic method is essentially similar to the above described magnifying mirror method, with the exception that the mineralized air bubble was drawn up into the bubble holder and subsequently discharged onto an externally placed watch glass, where the number of pick-up particles was counted with the aid of a microscope or hand lens. Drying was necessary when the discharged liquid interfered with the counting. The results of the bubble pick-up tests can be simply expressed in the form of the actual count. However, this can be misleading because the number of pick-up particles varies not only with the size of bubble but also with the size of particles. This short-coming is largely remedied by the use of the bubble pick-up index. as calculated from the following equation: PI = (pn/p) 100 [1] in which PI and pn are respectively the pick-up index and the number of pick-up particles under any test conditions, and p is the maximum number of pick-up particles of the same mineral under optimum flotation conditions, with the same size holder and particle. The bubble pick-up coefficient is a term coined to denote the difference between the bubble pick-up index of mineral particles in a liquid and that of the same mineral-liquid mixture plus a chemical reagent at the same pH value. For example, the difference between the pick-up index of a naturally floatable bituminous coal in a solution of water and kerosene and that of the same coal in water alone is a bubble pick-up coefficient. Compared with the contact angle method, the quantitative bubble pick-up method is inferior in accuracy but superior in sensitivity. Other advantages of the bubble pick-up method are that the apparatus can be easily assembled, the mineral sample is conventionally prepared, and the test procedure is relatively simple. References S. R. B. Cooke: The Flotation of Quartz Using Calcium Ion As Activator, Trans. AIME (1949) 184, p. 306. 'G. A. Last and M. A. Cook: Theory of Collector-Depressant Equilibrium at Mineral Surfaces. Paper presented at St. Louis Meeting, AIME (February, 1951).

Jan 1, 1954

LIBERIAN ore has been called the richest iron ore mined in the world, and thus far, the output from the Liberian Mining Co.'s Bomi Hills mine has lived up to its reputation. Iron content has been variously given as 68.47 pct to 68.825 pct. A recent United Nations survey gave it first place, surplanting the previously top-rated Swedish product. The first shipment of Liberian ore to the United States reached Baltimore in June, 1952 and went into blast furnaces at Republic Steel's Cleveland plant, where excellent results were reported. The physical make-up of the ore was hard, dense lumps, free of fines and with ideal size distribution. The Bomi hills mines were scheduled to be in full operation by the end of 1952. At the close of 1951, 200,000 tons of ore had been delivered to Baltimore. Liberian Mining Co. states that production is about 90,000 tons of ore per month at present. Between 1200 and 1500 men are employed at the mine site, with about 25 pct of them Americans. The Bomi Hill project was financed with the aid of the Export Import bank, which loaned Liberia Mining part of the sum needed to initiate and put the development into production. No definite expansion plans are in the making but extensive diamond drilling is being pursued. Original estimates of annual production for Bomi Hills ran in the neighborhood of 1.5 million tons, but the figure is no longer applicable. Liberia's potential as an iron producing country has been known for a long time, but the physical obstacles in the way of large scale production were such as to prohibit even contemplation of the project. Until recent years, the only route from the mine site to Monrovia was a foot trail. Once to the coast, over swamps and unbridged rivers, the ore was still useless because of the lack of loading facilities. Transportation Solved One problem was solved with the completion of the artificial harbor at Monrovia. Harbor construction was well under way when the mine project began in 1948. The island location of Monrovia also made the building of a bridge across the St. Paul River necessary. The 980 ft reinforced concrete structure, accommodating railroad and vehicular traffic, was opened in 1949. Meanwhile, preliminary work on railroad construction started. Grading was completed quickly. In April, 1951, the first train made the trip from Monrovia to the Bomi Hills. About half of the cost of the entire project went into railroad construction. More than 2200 men, often forced to resort to primitive methods, were engaged in clearing a path through the heavy forests between the mine and port. Swamps were drained and filled, and bridges, including the William V. S. Tubman, named for the Liberian President, were constructed. At the port, a Robbins Conveyor Belt Co., orehandling plant was installed. Ships can be loaded from the stockpile at the rate of 3000 tons per hr and vessels need spend only one day in port. Liberian Mining is operating under an 80 year lease, giving the company exclusive mining rights to all minerals, except gold, platinum, and diamonds, within a 40 mile radius of Bomi Hills. The contract provides a fixed royalty to the Liberian Government of five cents per ton of ore extracted, and a premium -based on the New York market price for Bessemer-grade pig iron. The entire project started because an American,

Jan 1, 1952

By W. J. Slatosky

As a means of effecting better control of endpoint temperatirres at the Jones & Laughlin basic oxygen furnace plant, a set of mathematical equations has been developed. The eqlutions are the product of a themlochemical anaysis of the process and aye designed to calculate the required scrap, lime, and hot metal additions in terms of a number of independent variables. Results of test heats have warranted adoption of this technique by the Prodrrction Department. BECAUSE of the autogeneous nature of the basic oxygen steel-making process, bath temperature can be controlled without an external fuel supply by charging the furnace with additions that are thermally balanced. The thermal requirements of the charge materials are such that, during the refining process, they throttle the heat generated by the metallurgical reactions in a manner designed to result in a speci-fied temperature at the completion of the heat. In the past, operating personnel at the basic oxygen furnace plant of Jones & Laughlin's Aliquippa Works relied on their experience and technical knowledge of the process to determine the quantities of charge additions needed to result in a finishing temperature in the range 2880"to 2920" F. (The charge consists primarily of 93 tons of scrap and hot metal plus an amount of lime sufficient to maintain a basicity ratio of 2.8 to 3.2). Estimates of these materials are based on a consideration of the effects on finishing temperature of 1) iron silicon content, having a variation of 0.8 to 1.8 pct; 2) iron temperature, ranging from 2250°to 2600°F; and 3)any excessive cooling of the furnace due to a production delay. The end temperature of the preceding heat also serves as a guide in that, if a heat was within the specified temperature range, the succeeding heat could be charged with materials of nearly the same proportions, provided the hot metal used in each of the two charges was of approximately the same temperature and composition. On the other hand, if a heat was outside the specified tapping range, or if the hot metal used in that heat was of different analysis and temperature from that of the iron to be charged, an adjustment in the proportion of additions is in order for the following heat. Due to the complex thermochemical behavior of the process and to the inexact and subjective nature of the described method of determining charge additions, consistently accurate temperature control was not to be expected. Therefore, those heats out- side the specified tapping range necessitated subsequent adjustments by either reblowing the cold heats for a suitable length of time so as to elevate the bath temperature to the desired level, or cooling hot heats with a proper amount of scrap. Because extra time is required to make these adjustments, production is delayed. In an attempt to devise a method for improving temperature control, an analysis of the thermochemistry of the process was undertaken. This, in turn, led to the development of a set of mathematical equations which enable the calculation of the quantities of scrap, lime, and hot metal needed to result in any specified tapping temperature range. The analysis was not intended to be a repetition of work done by others such as McMulkinl or ~hilbrook.' It was meant to be an extension of their work so that charge additions could be calculated not in terms of silicon alone but, rather, as a function of all independent variables. This paper presents the derivation of these relationships, their effectiveness in controlling bath temperatures, and a method of utilizing them on an operational basis. The Heat Balance—The first step undertaken in the analysis of the problem was the enumeration of the pertinent variables. A list is presented in Table I where it is noticed that these quantities have been separated into the following three categories: important variables, variables considered as constants, and variables to be neglected. The breakdown was an arbitrary one designed to facilitate the analysis; otherwise, the mathematical treatment would have been exceedingly cumbersome and complex. Fortunately, experience has shown that these simplifying assumptions do not seriously impair the accuracy of the calculations. These variables along with the limiting assumptions listed in Table n were then used to write a heat balance of the process by applying the equation of continuity, Rate of Rate of Rate of Increase = Income - Outgo PI ] of Heat of Heat of Heat.

Jan 1, 1962

By A. L. Mincher, W. F. Sheely

Mechanical properties of columbium have been studied over the temperature range of -196 to 1093oC. The decreased strengthening influence of cold-work at temperatures below ambient has been interpreted in terms of the Peierls-Nabarro effect. Maxima in the rate of strain hardening observed during tensile testing in the range 250-600°C. have been correlated with interstitial impurities to indicate the temperature ranges at which carbon, oxygen, and nitrogen, respectively, are responsible for strain aging. THE growing need for structural materials for use above the useful service temperatures of the iron-, nickel-, or cobalt-base alloys has caused the refractory metals to be considered as potential engineering materials. These metals, which include columbium, tantalum, molybdenum, and tungsten, are called refractory because the lowest melting point among them,that of columbium, is about 1000°C higher than the average melting temperatures of conventional high-temperature alloys. They are all body-centered cubic transition metals and, as such, their mechanical properties have basic characteristics which distinguish them from the face-centered cubic metals. For example, all show a much steeper rise in strength with decreasing temperature below room temperature than do the face-centered cubic metals, and their mechanical properties are strongly influenced by interstitially dissolved impurities. In order that these new metals may be used efficiently, it is necessary that their characteristics of behavior be fully known. In this paper, the mechanical properties of columbium will be examined over a wide range of temperatures. In particular, the influences of cold-work and individual species of interstitial impurity atoms on mechanical properties will be described, and basic mechanisms which may control the observed characteristics will be explored. EXPERIMENTAL The material used in this investigation was Union Carbide Metals Co. columbium roundels consolidated to four 4-in. diam ingots, three by consumable-electrode arc melting and one ingot by electron beam melting. Impurity contents of the ingots and methods of ingot conversion and treatment are summarized in Table I. The only metallic impurity occurring in any significant quantity was tantalum at about 0.1 pct. Iron, silicon, titanium, and zirconium were each less than 0.015 pct; boron was 1 ppm or less. This should have no appreciable influence on properties. The electron beam melted material, being the purest, will be used as the basis for comparison in the discussions to follow. Tensile tests were conducted from-196 to 1093oC, on both cold-worked and fully recrystallized arc-melted and electron-beam melted columbium using standard 1/4-in. diam, 1-in. long gage length test specimens. A strain-rate of 0.005 in. per in. per min was employed until the 0.2 pct yield strength was achieved and then the strain-rate was increased to 0.05 in. per in. per min for the balance of the test. Samples were protected in an inert atmosphere at tests above 300°C. The tensile properties obtained on the electron-beam melted columbium, E, in both the cold-swaged and recrystallized conditions are given in Fig. 1. The yield strength data of Dyson, et al.,' obtained on recrystallized electron beam melted columbium and the tensile strength data reported by Tottle2 on powder metallurgy columbium are included in Fig. 1. The material used by Tottle had been purified by vacuum sintering. There is excellent agreement between Dyson's data and those obtained in the present investigation. The tensile strengths obtained by Tottle were slightly greater than those obtained in this investigation on electron-beam melted columbium but varied with temperature in a similar manner. Tottle's data showed a maximum in tensile strength near 500°C, as did our data on electron-beam melted material, and also showed a small maximum at 300°C. The significance of these maxima will become evident later in the discussion. The tensile properties of cold-swaged and recrys-tallized arc melted columbium are plotted in Fig. 2. It was found that the properties of the recrystallized arc-melted columbium from all three heats showed very close agreement except at temperatures between about 500" and 800°C. A reason for this range of disagreement will be suggested in the discussion. The generally good agreement, however, attests to the ability of cold-working and subsequent recrystal-lization to erase the effects of the three different primary breakdown procedures and to produce nearly equivalent structures in the samples derived from the three different heats. wesse13 reported tensile data on columbium having interstitial impurity contents between those of the

Jan 1, 1962

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner' makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,' is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors'' derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel' terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones'" initiates steady-state

Jan 1, 1967

By Harry A. Wahl

Precent fracture design procedures are bared on the total fracture area created. A method to distinguish beI,,.ecn total area and [he propped or effective fracture area has not been available. This paper presents a solution to this problem, applicable to horizontal fractures. The difierences between effective fracture area and torn] area are demonstrared in example calculations. This work is hayed on experimentally determined transport efficiencier of solids in sand-liquid slurries. Newtonian and non-Ne~vtoninn systems are considered. INTRODUCTION Fifteen years after commercial introduction, hydraulic fracturing remains the most successful stimulation technique in the oil field. This success is primarily due to ability of induced fractures to penetrate and alter permeabilities deep within formations. Many fields producing today could not have been developed without the hydraulic fracturing process. Because of wide usage, fracture-treatment design has received a great deal of engineering and research effort. This work, resulting in improved equipment and materials, has increased the benefits from fracture treatments as well as the applicability of the process. A major contribution was the development of fluid-loss additives. Necessarily, the number of parameters to be considered in treatment design has steadily increased, resulting in more complicated design techniques. Almost all present design procedures are based on the precepts set forth by Howard and Fast. Relating the fluid volume lost into the formation, the volume required in extending the fracture, and the total slurry volume injected, they developed an expression for the total fracture area created in terms of pertinent treatment parameters. Fluid loss during treatment was expressed as a function of time for three flow mechanisms. Although modifications of fluid loss equations have been made, the total fracture area concept has remained unaltered. A vast amount of field data indicate that induced fractures must be propped and held open to be effective. A notable exception is the Mesa Verde formation in the San Juan basin. However, analysis of these treatments shows that improved well productivities are obtained when propping agents are incorporated in the treating fluid. Although propped fracture area has been recognized as an important design parameter, a method to distinguish between total area and effective fracture area has not been available. The necessary information on slurry-sand transport in fractures has been lacking. Interest in the propped region of induced fractures is not restricted to areal extent alone. The distribution of sand within fractures is important from the standpoint of fracture flow capacities. Flow capacity affects the increase in well productivity after stimulation. The work of Huitt and Darin4 hows that partial monolayers of sand have large flow capacities compared to thick* dense sand packs. It has been postulated that gelled fluids have the ability to transport sand within the fractures at the deired low concentrations. An early contribution in the area of sand placement in fractures was made by Kern et al.' They studied sand movement in a transparent vertical fracture model. It was observed that the sand tended to settle out in the bottom of the model before moving very far. When the fluid velocity exceeded a certain critical value, all of the sand injected began moving through the crack even though it settled to the bottom. This critical velocity was determined under several flow conditions. Some work on sand movement in horizontal fractures has been reported in Russian publications. Sand movement was studied by Izyumova and Shan'gin' using a transparent "pie-shaped" flow model to simulate a horizontal radial flow system. However, the data were limited, especially in a quantitative sense. Dorozhkin, Zheltov and Zheltow studied the behavior of sand-liquid slurries in a horizontal linear flow model. The quantitative data were restricted primarily to the thickness of sand deposits formed at the bottom of the fracture. An earlier paper provided basic data on the flow of sand in horizontal fractures. This study was designed to yield specific quantitative information on rate of advance of sand particles and pressure behavior under various flow conditions. A comprehensive photographic study was undertaken in a 10-ft windowed flow cell to provide the necessary qualitative and quantitative data. Since the number of potential variables far exceeded the capacity of the initial study, emphasis was placed on the effects of sand concentration, oil viscosity and oil flow rate. A detailed description of these experiments and the results are described in Ref. 9. However. the implications of this work on the fracture design calculations were not discussed. An analysis of these data as well as new data is provided in the following sections. EXPERIMENTAL RESULTS The primary objective of the experimental investigation was to provide information on the rate of advance of the solids in sand-liquid sturries. A 10%-ft long transparent

Jan 1, 1966

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.' ' Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.'" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By R. G. Agarwal, Ramey Jr. H. J., Al-Hussainy R.

Although it has long been realized that gas recovery from a water-drive gas reservoir may be poor because of high residual saturations under water drive, it appears that only limited infomlation on the subject has been available until recently. This study was performed to show the qiiantitative potential importance of water influx. Results indicate that gas recovery may be very low in some cases: perhaps as low as 45 per cent of the initial gas in place. Gas recovery under water drive appear to depend in an important was on: (I) the prodirction rate and manner of production; (2) the residual gas saturation; (3) aquifer propertie.); and (4) the volumetric displacement effciency of water invading the gas reservoir. The manner of estimating water-drive gas reservoir recovery can vary considerably. Examples are: the steady-state tnethorl. the Hurst modified steady-state method, and various unsteady-state methods such ac. those of van Ever-dingen-Hurst, Hurst, and Carter-Tracy. The Carter-Tracy water influx expression was used in this study. In certain cases, it appears that gas recovery can be increased significantly by controlling the production rate and manner of production. For this reason, the potential importance of water influx in particular gas reservoirs should he investigated early to permit adequtrtr planning lo optirtize the pay reserves. INTRODUCTION In recent years, the economic importance of natural gas production has become increasingly apparent. This has been evidenced by more intensive exploration efforts aimed at gas production, and exploitation of both deep, as well as low-permeability gas reservoirs. Technical developments such as deep-penetration fracturing have made development of such formations economically feasible. Unfortunately, water influx has forced abandonment of a number of gas reservoirs at extraordinarily high pressures. Although reservoir engineering methods for estimating water influx have long been available, it appears that application of these methods to the water-drive gas reservoir has been sporadic.'a Available methods for estimating water influx which can be applied to the water-drive gas reservoir problem include the steady-state method,1 the Hurst modified steady-state method and various unsteady-state methods such as those of van Everdingen-Hurst. Hurst, and Carter-Tracy. Interesting applications of these solu- tions to gas reservoir and the aquifer gas-storage problems have appeared recently.3,12,14 The experimental study of residual gas saturations under water drive by Geffen et al. in 1952 indicated that residual gas saturations could be extremely high. A value of 35 per cent of pore volume is often used in field practice when specific information is not available. The study of Geffen et al. showed that residual gas saturation might be much higher in some cases. Naar and Henderson concluded that the residual non-wetting phase saturation under imbibition should be about half of the initial non-wetting phase saturation. The Naar and Henderson result that residual gas saturation under water influx should be about half the original gas saturation is recommended as an estimate if laboratory measurements are not available. Thus, it is clear that a considerable portion of the initial gas in place might be trapped in a water-drive gas reservoir as residual gas at high pressure. A full water-drive would result in loss of residual gas trapped at initial reser.voir pressure. Consideration of transient aquifer behavior leads to the conclusion that high-rate production of water-drive gas reservoirs could result in improved gas recovery by reduction of the abandonment pressure. However, there appears to be little quantitative information on this possibility. One of the few advantages of water-drive gas production appears to be improved deliverability through water-drive support of the reservoir pressure. There may also be an advantage in higher condensate recovery caused by pressure maintenance for gas-condensate water-drive reservoirs. In view of the preceding, this study was made to assess the potential importance of water-drive in gas reservoir engineering. The Carter-Tracy approximate water-influx expression was used because this equation offers some advantages in hand-calculation which do not appear to have been generally recognized.' However, calculations were performed in the main with a high-speed digital computer to permit evaluation of the effect of water-drive under a large variety of conditions. CALCULATION METHOD Water-drive gas reservoir performance can be estimated in a manner completely analogous to oil reservoir calculations: a materials balance is written for the reservoir, and a water influx equation is written for the aquifer. Siniltaneous solution provides the cun~ulative water influx and the reservoir pressure. When reservoir performance data (gas produced and reservoir pressures) are available, it is usually possible to match performance data to determine the initial gas in place and the water influx parameters —

Jan 1, 1966

By R. E. Schilson, F. H. Poettmann

The direct determination of the stabilized performance behavior of low capacity, slowly stabilizing gas wells is extremely time-consuming and wasteful of gas. From both field experience and theoretical considerations, a test procedure has been evolved by which the stabilized hack pressure behavior of such gas wells can he predicted without having to revert to long time flow tests. The method consists of using the isochrona1 test procedure to establish the slope of the back pressure curve, "n", and the short time variation of the performance coeficient, "C", with time. From this short time transient flow data and theoretical considerations, the value of C at large times can he established. By assuming the radius of drainage of a well to be half the distance between wells, one can calculate the stabilization time for various well spacing patterns. Once the stabilization time for a given .spacing has been determined, the value of C can be calculated and the stabilized back-pressure curve can he. establbrhed. The calculated perfortnance coefficient as a function of tinze was cotnpared to the experintentally measured values for a number of gas wel1.e. The deviation of the calc~*lated from the cxperimental res1tlt.s vary depending on the set of short time experimental points used to evabrate the parcrttzeters of the equation. The longer the tirne far the flow test data user1 in the calculations, the better was the agreement with the experimental results. The time necessary to obtain this data from well tests varies consitlerably, depending on the physical natrcre oj the rt3scrvoir under consideratiot~. INTRODUCTION For many years, the U. S. Bureau of Mines Monograph 7' has served as a guide for testing and evaluating the performance of gas wells by means of the back-pressure method. The back-pressure performance of a gas well is expressed by the following equation: Q~ CW-W.........(I) where the characteristics of the back-pressure equation are determined by C, the performance coefficient, and n, the exponent which corresponds to the slope of the straight line when Q and {P* - PN2) are plotted on logarithmic paper. Q is gas flow rate at standard conditions, and P, and P, are equalized and flowing bottom-hole pressures, respectively. Prior to the development of the back-pressure tcst, the "open flow" capacity method of testing a well was common. By this method, the wells were flowed wide open to the air and the flow rate measured. Such procedure was wasteful of gas and did not provide information on the deliverability of the gas to the pipe line. Monograph 7 Procedure The back-pressure method of testing wells was developed to overcome these shortcomings. Although much has been learned regarding the laws of the flow of gas through porous formations, the original development of the back-pressure relationship was based entirely on empirical methods. The back-pressure behavior provides the engineer with information essential in predicting the future development of a field. It permits him to calculate the deliverability of gas into a pipe line at predetermined line pressures, to design and analyze gas gathering lines, to determine the spacing and number of wells to he drilled during the development of a field to meet gas purchasers' requirements, and to solve many other technical and economic problems. As described in Monograph 7, the flow-after-flow method of back-pressure testing, when applied to fast stabilizing and usually high capacity wells, correctly characterized the behavior of the well. However, as the value of the gas at the wellhead increased, small capacity gas wells having slow rates of stabilization became economically operable. The flow-after-flow method of testing could not be used to describe the behavior of these slowly stabilizing wells. The procedure of Rawlins and Shellhardtl for establishing the back-pressure behavior of a gas well was based on the rcquirement that the data be obtained under stabilized flow conditions; that is, that C is constant and does not vary with time. C depends on the physical properties of the reservoir, the location, extent and geometry of the drainage radius, and the properties of the flowing fluid. In a highly permeable formation, only a very short period of time is required for the well to reach a stabilized condition, and, consequently, the requirements for the test procedure described in Monograph 7 are met. For a given well, n is also constant

By E. Miller, J. M. Eldridge, K. L. Komarek

The liquidus curve of the Mg-Pb system was accurately redetermined. The compound Mg2Pb decomposes peritectically at 538.2° ± 0.3°C to liquid and to a compound p' which melts congruently at 35.0 at. pct Pb and 549.0° ± 0.3°C. The solidus curve of ß' was determined. X-ray diffraction studies indicate that 4' has an orthorhombic structure. Activity values of magnesium calculated from the phase diagram agree with those published in the literature. EXPERIMENTAL thermodynamic properties of binary metallic systems have to be consistent with values calculated from the phase diagram. In systems forming intermetallic compounds the shape of the liquidus curve near a compound is determined by the thermodynamic properties of the coexisting solid and liquid phases. Hauffe and Wagner' neglected the temperature dependence of the chemical potentials and obtained the potential differences of the components of the liquid alloys, relative to stoichiometric liquid. Their calculations were based on the liquidus curve and on the heat of fusion of the compound, and were only valid near the congruent melting point. Steiner, Miller, and Komarek2 developed equations which account for the temperature dependence and obtained the chemical potentials of liquid Mg-Sn alloys over the entire phase diagram from the liquidus and solidus curves and from enthalpy values with the pure components as the standard states. The Mg-Pb phase diagram has been studied by several investigators whose results have been compiled and critically evaluated by Hansen.3 Although the liquidus curve was poorly defined, the general features of the diagram, i.e., one congruent melting compound, Mg2Pb, of essentially stoichiometric composition, two eutectics, and limited terminal solid solubilities, seemed to be suitable for a similar thermodynamic analysis. A careful redeter-mination of the liquidus by thermal analysis revealed, however, the existence of another compound. The liquidus curve between the two eutectics was precisely delineated and the structure and solidus curve of the new compound were investigated. The revised phase diagram was thermodynamic ally analyzed to evaluate the activity of magnesium in the liquid alloys. EXPERIMENTAL PROCEDURE The magnesium metal (Dominion Magnesium Ltd., Toronto, Canada) had a purity of 99.99+ pct; lead (American Smelting and Refining Co.) contained 99.999 pct Pb. Most experiments were carried out in graphite crucibles. Several experiments were made in high-purity alumina (Triangle R.R., Mor-ganite, Inc.) and in Armco iron crucibles to test the inertness of the graphite crucibles. Chemical analysis of magnesium and detailed description of the procedure for thermal analysis have been given previously. For the determination of the solidus curve of the compounds, specimens of initial composition Mg2Pb were equilibrated in a closed isothermal system with magnesium vapor. The source of the magnesium vapor was an alloy which had a gross composition lying in the 0' + L field at the temperature of equilibration. As equilibrium was approached, the specimens lost magnesium to the two-phase reservoir thereby lowering the activity of magnesium in the specimens until activity and composition equaled that of the ß'/ß' + L boundary. Crucibles (1.9 cm ID by 2.2 cm OD by 4.1 cm high) and tightly fitting lids were machined from a molybdenum rod; small, shallow trays were fashioned from thin (0.005 in.) molybdenum sheet, and all the molybdenum components were degreased in hot carbon tetrachloride and then dried. The pieces were then degassed in vacuum at 950°C for about 6 hr. The two-phase alloy was placed at the bottom of the crucible and small specimens of the Mg2Pb compound, weighed on an analytical balance, were placed in two molybdenum trays above the two-phase alloy. The crucible was closed by forcing its lid on and then inserted in a titanium crucible. This crucible was evacuated, flushed twice with argon, and welded under argon. The specimens were equilibrated for about 1 week in a resistance furnace regulated by a Celectray controller, and the runs were terminated by water quenching. The specimens were again weighed and the equilibrium compositions were calculated on the basis that the weight losses were solely due to a loss of magnesium to the two-phase alloy. The structure of the B' phase was investigated by the Debye-Scherrer X-ray diffraction technique. Selected ingots from thermal-analysis experiments containing about 35 at. pct Pb were re-melted, slowly cooled, and crushed in an argon-filled glovebox until the entire ingot passed through a 50-mesh sieve. The powder was thoroughly

Jan 1, 1965

By G. L. Fisher, R. H. Geils

The effect of small amounts of columbium (<0.01 to 0.10 pct) on the ?-a transformation occurring during the continuous cooling of a low carbon Ni-Cu steel was investigated. Dilatometer specimens were aus-tenitized at 950" and 1068?C and cooled at 17? and 375 C° per min. Columbium caused a marked depression in the ?-a transformation temperature except when cooling at the slower rate from 950°C. The effec of columbium on the transformation temperature was greater the higher the austenitizing temperature and rate of cooling. A maximum depression of 92 C" was observed. Metallographic examination of specimens of <0.01 and 0.07pct Cb steels heated at 1200°C for 1 hr and cooled at various rates showed that columbium had a major effect on the ferrite morphology. The fer rite in the columbium -free steel remained equiaxed at cooling rates as high as 440 C? per min while the columbium-bearing steel exhibited mixed structures o equiaxed and bainitic ferrite at cooling rates as low as 130 C° per min. The ? grain boundaries in the columbium -free steel provided the ferrite nucleation sites in rapidly cooled specimens. There was a complete absence of nucleation at these sites in the colum bium-bearing steel. It is concluded that columbium depresses the transformation temperature by suppressing ferrite nucleation at the austenite grain bound-aries. In this respect the effects of columbium are analogous to those of boron in low C-Mo steels. It is well known that small columbium additions can substantially strengthen plain carbon steels. As little as 0.02 pct Cb can increase the yield strength of mild steels by 10,000 psi.1 A fine precipitate of CbC has been observed in columbium-bearing steels2 and is generally thought to be responsible for the strengthening. Little attention has been devoted to the effect of columbium on the ?-a transformation. Webster and woodhead3 have studied the effect of columbium on the isothermal proeutectoid ferrite reaction in mild steels. They found similar transformation behavior in steels both with and without columbium additions. However, as the austenitizing temperature increased, the incubation time for the start of the ferrite transformation became longer in the columbium-containing steel. Morrison1 found that the addition of 0.03 pct Cb to a C-Mn steel lowered the transformation temperature by 50 C° during cooling from 1200°C at a rate of 80 C" per min. The strengthening effect of columbium has recently been utilized in an age-hardenable, low-alloy steel containing copper and nickel.4 A small amount of columbium has a substantial effect on the as-rolled strength of this steel. By increasing the columbium level from <0.01 to 0.13 pct the as-rolled yield strength is increased by 15,000 psi. Columbium also significantly lowers the ?-to-a transformation temperature of this steel during continuous cooling from the austenitizing temperature. Because of the low carbon level in this steel (0.05 pct max), it is almost entirely ferritic. Thus, it offers the opportunity of studying the effect of small columbium additions on the proeutectoid ferrite reaction. Of particular interest in this study was the reason for the marked lowering of the transformation temperature by columbium during continuous cooling. EXPERIMENTAL PROCEDURE Materials. The compositions of the steels used in this investigation are shown in Table I. The steels were 30-lb air induction melts. They were forged to 4 by 8 by 1 in. plate at 1230°C, air cooled, and then reheated to 1230°C and cross-rolled in two passes to in. plates. Dilatometry. A Leitz Bollenrath dilatometer was used to record the transformation during continuous cooling from two different austenitizing temperatures. The dilatometer specimens were + in. in diam and 2 in. long. Oxidation and decarburization of the specimens was prevented by maintaining a small positive pressure of dry argon in the dilatometer furnace and by plating the specimens with 1 mil of Cu. For the lowest cooling rate, 17 C" per min, the temperature of the specimen was measured with a Pt-Pt 10 pct Rh thermocouple placed in a & in. diam well in the center of the specimen. During air cooling, 375 C° per min, this method of measuring the temperature interfered with the operation of the dilatometer. However, it was found that the temperature of the specimen could be measured accurately by placing a thermocouple in an identical specimen in a holder adjacent to the one being used to operate the dilatometer mechanism. The dilation-temperature curves were recorded on photographic film and then converted to volume percent ferrite-vs-temper-ature curves. The cooling rates obtained with the dilatometer are shown in Table 11. Cooling rates of

Jan 1, 1970

By B. A. Eaton, K. E. Brown, C. R. Knowles, D. E. Andrews, I. H. Silberberg

This paper presents the resitlts of an investigation of two-phase, gm-liquid flow in horizontal pipelines. Experimental data were taken in three field-size, horizontal pipelines, two of which were constructed for this purpose. The data were obtained using water, distillate and crude oil separately as the liquid phase, and natural gas as the second phase. Experimental pressure-length traverse, liquid holdup and flow-pattern data were collected for each set of flow rates. These data were used to develop three correlations that are presented herein. The liquid-holdup values correlated with various flow parameters without regard to the existing flow pattern. The same was true for the energy-loss factors. A new flow-pattern map is presented that appears to be quite reliable, but not required for the pressure-loss calculations. The liquid-holdup correlation and the energy-loss factor correlation are used in conjunction with a two-phase flow power balance, developed during this study, to predict the pressure losses that occur during gas-liquid flow in horizontal pipelines. A recommended calculational procedure is given, as well as a statistical analysis of the results. This procedure lends itself to computer application, since several small pressure decrements are needed to calculate a pressure-length traverse. The correlations are shown graphically, but may be curve fitted with existing curve-fitring computer programs. INTRODUCTION Due to the frequent occurrence of gas-liquid flow in pipelines and the desire to accurately calculate the pressure losses that occur in these lines, two-phase flow is of considerable interest to the petroleum, chemical and nuclear industries. In the petroleum industry, gas-liquid mixtures have been transported over relatively long distances in a common line due to the advent of centralized gathering and separation systems. Long two-phase flowlines are usually accompanied by large pressure drops which influence the design of the system. Gas-lift installations are designed on the basis of known tubing pressures at the wellheads. The horizontal flowline connecting the wellhead and the separator system must be correctly sized in order to minimize the horizontal flowline pressure losses and the wellhead tubing pressure. Practically all oilwell production design involves horizontal two-phase flow in pipeknes. All of the flow processes of oil and gas production must be studied simultaneously to insure good well design. Since the beginning of offshore oilfield development, long horizontal flowlines have been constructed. Because pressure losses greatly influence the performance of producing wells, a method is desired that can be used to predict such pressure losses and select optimum flowline size. Several types of gas-liquid flow exist, and many of these are discussed by Gouse The study of pressure gradients, fluid distributions and flow patterns that occur in horizontaI multiphase flow is made difficult by the great number of variables involved. The various flow regimes give rise to changing velocities of the fluid particles in all directions. These instabilities of the interface between the gas and liquid prohibit the determination of actual vector velocities of fluid particles in each phase. Also, it is practically impossible to arrive at correct sets of boundary conditions. Therefore, most investigators have concluded that a solution to the problem by the classical fluid dynamics approach, whereby the Navier-Stokes equationsM are formulated and solved, is far too complex. Other methods must be utilized to develop general correlations that will predict the behavior of gas-liquid horizontal flow systems. Multiphase flow studies have sought to develop a technique with which the pressure drop can be calculated. Pressure losses in two-phase, gas-liquid flow are quite different from those encountered in single-phase flow; in most cases an interface exists and the gas slips past the liquid. The interface may be smooth or have varying degrees of roughness, depending on the flow pattern. Therefore, a transfer of energy from the gaseous phase to the liquid phase may take place while energy is lost from the system through the wetting phase at the pipe wall. Such an energy transfer may be either in the form of heat exchange or of acceleration. Since each phase must flow through a smaller area than if it flowed alone, amazingly high pressure losses occur when compared to single-phase flow. Most investigators of horizontal two-phase flow phenomena have chosen to separate their experimental data into several groups of observed flow patterns or regimes.

By P. A. Witherspoon, R. W. Donovan, T. D. Mueller

The use of petroleum-barren aquifers for underground storage has become extremely important to the natural-gas industry. A critical problem in assessing the feasibility of a specific aquifer for such use is the permeability determination of the caprock over the proposed storage project. The approach used here is to conduct both static and dynamic field tests on the aquifer being analyzed. Valuable information on the possibility of communication between the storage aquifer and any other aquifers above can be obtained by measuring hydrostatic water levels and water analyses. Significant differences in such data give evidence of the lack of communication between the intended storage reservoir and other horizons. The dynamic approach requires that one well be pumped in the storage aquifer, and changes in fluid levels recorded in both the aquifer and its caprock. The interpretation of the data from such pumping tests involves the solution of nonsteady radial flow in an infinite aquifer and the influence on such flow of a leaky caprock. A finite-difference model has been used to investigate this problem, and the transient behavior has been solved numerically with a digital computer. It has been found that the pressure transients in the storage aquifer are not affected significantly by moderate caprock leakage. The pressure behavior of the caprock is a much better indicator of the degree of leakage, and generalized solutions for this behavior are included. Field data are presented to demonstrate both the static and dynamic approach. If is concluded that appropriate investigation of the groundwater hydrology in an aquifer-type gas-storage project can provide much valuable information for determining the effectiveness of the caprock to hold gas. INTRODUCTION Underground storage of natural gas in the United States has been developing at a rapid rate over the past few years. In 1955, the total gas-storage capacity was about 1.6 trillion cu ft; by 1961, this figure was almost 3.2 trillion cu ft, an increase of 100 per cent in six years.&apos; This trend un- doubtedly will continue because the economics favor the development of gas storage, as opposed to the construction of new pipelines, to meet the inherent cyclic demand for fuel in the metropolitan areas of this country.&apos; About 15 per cent of the current underground gas storage has been developed in petroleum-barren aquifers, i.e., geological domes or anticlines in which no commercial quantities of oil or gas had been produced prior to the storage operations. The necessity for using barren aquifers outside many metropolitan areas of this country has been due to the lack of depleted oil or gas fields that were near enough and large enough to meet the demands of such consuming areas. Pipeline companies have developed aquifer storage along their transmission lines to meet the fluctuating needs of their complex systems. Considerable thought has also been given to the problem of storing gas in a structureless aquifer, both in this country&apos; and in the Soviet Union outside the city of Leningrad.&apos;," Conditions such as these have led to the development of aquifer gas-storage projects in many parts of the U. S. Most of these developments have centered in the Mid-Continent area, and the greatest amount of activity has been concentrated in Illinois.6 Thus, the use of petroleum-barren aquifers for gas-storage purposes has become extremely important to the natural-gas industry. There are three basic problems in developing aquifer-type storage: (1) finding an adequate geologic structure, (2) finding a suitable storage reservoir within the structure and (3) determining the tightness of the caprock over the intended storage zone. The first two problems can be solved by applying conventional methods of exploration geology, but once these problems are solved, the question arises as to why no oil or gas is present in an otherwise favorable setting. Two situations are possible: (1) an adequate source bed was never present, or (2) a source bed was present but the petroleum seeped away because of a leaky caprock. Determining the tightness of the caprock is one of the most critical problems in assessing the feasibility of a specific aquifer for storage purposes. In attacking this problem, one usually takes cores of the caprock and subjects them to a rigorous investigation. Such core data are desirable, but they only detail the matrix properties and cannot be expected to reveal the gross characteristics of the caprock. Several gas-storage projects in the U. S. have had considerable leakage where

By J. R. Elenbaas, J. A. Vary, D. L. Katz

The development of gas fields, oil fields and aquifers for storing natural gas is treated from two main vieu.-points: (I) the volumetric storage capacity for gas in a given situation and (2) the prediction of the number of wells required for the delivery of gas. Other experiences in the design and operation of storagc fields are incluclerl INTRODUC TION Storage of natural gas in underground reservoirs near the terminus of long distance pipelines has been the prime factor in opening the space heating market to the natural gas industry. Storagc has permitted a major. increase in both the load and the load factor of pipe-lines; some are now operating at steady load throughout the year. Thus, underground storage has been responsible for the rapid increase in demand for natural gas in recent years. Three types of reservoirs have been used for gas storage: natural gas reservoirs, oil reservoirs, and waterbearing sands or aquifers. This paper presents the factors to be considered when developing gas storage reservoirs of these three categories. There are two prime considerations tor any storage reservoir: (1) the volume of gas which a given reservoir will store advantageously and (2) the number 01 wells needed to provide the required peak deliverability. These two problems will be considered for the three types of reservoirs just noted STORAGE IN PARTIALLY DEPLETED GAS FIELDS Early storage operations consisted of replenishing the natural gas in a depleted gas field situated adjacent to the market. Today, newly discovered fields near the market may be considered for storage, and this discussion applies equally to both types of reservoirs. For reservoirs originally containing gas or oil, the question of the impermeability of the cap rock nor-mally does not arise. However, such fields are likely to have many wells drilled either to or through the reservoir under consideration. Positive assurance must be obtained that such wells are or can be made mechanically tight. Corroded casings may need to be lined or permanently plugged. Abandoned wells should bc reopened and properly cemented. The volumetric capacity for gas storage depends upon space available in the porous rock as well as pressure and temperature of the gas in the reservoir. The production-pressure decline data on partially depleted gas reservoirs without water drive permit calculation of the reservoir space for gas. Isopachous maps of sand volume and porosity data for the reservoir rock provide an alternate method of calculating the pore volume for water-drive reservoirs. The pressure range selected for the storage cycle depends upon ()) the safe upper limit of pressure. 2) the flow capacity of wells and (3) compression requirements when injecting gas into the reservoir or delivering to market. Normally, gas and oil fields have pressures at discovery in the range of 0.43 to 0.52 psi/ft of depth. Pressures of around 1.0 to 1.2 psi/ft of depth appear to lift the overburden1-3 and invite uncontrolled movement of fluids in the porous rock. Some top pressure is normally selected for a storage reservoir ranging from below discovery pressure for deeper reservoirs to 0.65 psi/ft of depth for shallower reservoirs. Pressures to 0.66 psi/ft have been experienced without difficulty. The lower pressure limit is set by water intrusion accompanying low pressures, reduced flow capacity for wells at lower pressures and compression requirements. Depletion-type gas reservoirs often encounter water problems in the later stages of gas production. Such water intrusion may be due to movement from the surrounding aquifer. Accordingly, displacement of this water back into the aquifer by gas pressure and subsequent surges of water corresponding to the gas storage pressure cycle must be considered. Storage fields often produce in four months a volume of gas equal to its initial content. Rapid decreases in reservoir pressure occur, such as 20 psi/day. Accordingly, closed-in pressure observation wells which reflect the pressure in the bulk of the reservoir are required for following the operation of the reservoir. It has been found that a plot of observation wellhead pressures against gas content, Fig. 1, is very useful in observing operation of the field, checking the inventory and predicting future behavior. The plot is based on a given quantity of base or cushion gas in place. The injection and withdrawal curves may spread depending upon the homogeneity of the reservoir rock. permeability of the rock, well spacing and flow rates.

By W. E. Bonini, E. A. Hickok

IN the course of an investigation directed toward expanding ground-water facilities in Essex and Morris counties, New Jersey, the Board of Water Commissioners of the city of East Orange authorized a seismic-refraction survey&apos; for the purpose of de-lineating bedrock topography below unconsolidated overburden. Results of the survey were highly satisfactory and led to the preparation of a comparatively detailed bedrock contour map. Knowledge of the bedrock depth and configuration was an important aid in selection of sites for test drilling. The portion of the East Orange Water Reserve under consideration is in the flood plain of the Pas-saic River about 10 miles west of Newark, N. J. The flood plain is about 175 ft above mean sea level and is bordered by low hills rising to elevations of approximately 250 ft. The bedrock underlying the Water Reserve consists of sandstone and shale of the Triassic Brunswick formation and is covered everywhere by deposits of unconsolidated glacial outwash sand and gravel, lacustrine clay, and recent river silt as much as 150 ft thick. Yield of wells in the sandstone and shale averages 100 to 200 gpm. Since production wells constructed in the sand and gravel aquifer in the buried river valley shown on the contour map (Fig. 1) yield 300 to 1400 gpm, it was proposed to locate additional production wells in this buried valley, where the yields per well would be maximum. In 1939 and 1946 the East Orange Water Dept. had electrical-resistivity surveys made to determine depths to bedrock. From the resistivity data the exploration company prepared a bedrock contour map. A well field expansion program begun in 1955 utilized this information to locate sites for test wells along a predicted northward extension of the buried valley in which existing production wells are located. After several test wells (wells 201-205) had been drilled, it became apparent that the resistivity information was unreliable." For example, test well 201 recorded bedrock at a depth of 72 ft, whereas the resistivity depth determination was 130 ft. As a consequence, the test drilling program was temporarily suspended and a seismic survey was under- taken to determine the topography and extent of the buried valley known from well records to underlie the existing well field. In the first phase of this study, several seismic shot point locations were placed at sites where well logs had been obtained previously. This procedure is necessary in a new area to determine whether the seismic method is applicable and what degree of accuracy is to be expected. At the East Orange Water Reserve, depths obtained from the shot points near test wells 202, 203, and 204 were within 8 to 11 pct of the depths logged (Table I). With this assurance that accurate results could be obtained, additional seismic spreads were located on the Water Reserve. Using a portable refraction seismograph, in the fall of 1955 a crew of four men shot a total of 29 reversed seismic spreads in a period equivalent to six field days. Charges as heavy as 3 1b of 40 pct dynamite were necessary at a few places to overcome ground vibrations caused by traffic on nearby highways. At most other sites, a 1-1b charge was sufficient. Travel-time plots were made for all spreads, and depths and true velocities were calculated according to formulas for multiple sloping layers by Ewing, Woollard, and Vine.&apos; The plot of spread 7 (Fig. 2) is typical of the short spreads in which bedrock was shallow—about 50 ft in this case. Where there were not enough arrivals through the bedrock to define the high velocity bedrock line, the spreads were lengthened. This was done by placing shots on line several hundred feet away from each end of the line of geophones. It was then possible to construct complete reverse plots for both short and extended shot points (see spread 27, Fig. 3). Four individual depths were calculated from each extended spread. Three and in some cases four seismic layers were observed. The surficial layer had a velocity range of 900 to 1200 fps, the lowest velocity recorded. This seismic layer is above the water table and is interpreted as recent river silt. The bedrock had the highest velocities, which ranged from 10,600 to 16,400 fps. Intermediate velocities ranged from 4500 to 6800 fps. In every case the intermediate layer was within

Jan 1, 1959

By F. C. Bond

SIZE of grinding media is one of the principal factors affecting efficiency and capacity of tumbling-type grinding mills. It is best determined for any particular installation by lengthy plant tests with carefully kept records. However, a method of calculating the proper sizes, based on correct theoretical principles and tested by experience, can be very helpful, both for new installations and for guiding existing operations. As a general principle, the proper size of the make-up grinding balls added to an operating mill is the size that will just break the largest feed particles. If the balls are too large the number of breaking contacts will be reduced and grinding capacity will suffer. Moreover, the amount of extreme fines produced by each contact will be increased, and size distribution of the ground product may be adversely affected. If the balls added are too small, grinding efficiency is decreased by wasted contacts that are too weak to break the particles nipped; these largest particles are gradually worn down in the mill by the progressive breakage of corners and edges. Ball rationing is the regular addition of make-up balls of more than one size. The largest balls added are aimed at the largest and hardest particles. However, the contacts are governed entirely by chance, and the probability of inefficient contacts of large balls with small particles, and of small balls with large particles, is as great as the desired contact of large balls with large particles. Ball rationing should be considered an adjunct or secondary modification of the principle of selecting the make-up ball size to break the largest particle present. Empirical Equation In 19521,2 the author presented the following emerical equation for the make-up ball size: B - ball, rod, or pebble diameter in inches. F = size in microns 80 pct of new feed passes. Wi - work index at the feed size F. Cs - percentage of mill critical speed. S — specific gravity of material being ground. D == mill diameter in feet inside liners. K - 200 for balls, 300 for rods, 100 for silica pebbles. Eq. 1 was derived by selecting the factors that apparently should influence make-up ball size selection and by considering plant experience with each factor. Even though Eq. 1 is completely empirical, it has been generally successful in selecting the proper size of make-up balls for specific operations. But an equation based on theoretical considerations should be used with more confidence and have wider application. The theoretical influence of each of the governing factors listed under Eq. 1 was accordingly considered in detail, as described below, and a theoretical equation for make-up ball sizes was derived. Derivation of Theoretical Equation Ball Size and Feed Size: The basis of this analysis is that the largest ball in a mill should be just sufficient to break the largest feed particle into several pieces, excluding occasional pieces of tramp oversize. In this article the size F which 80 pct passes is considered the criterion of the effective maximum particle feed size. The smallest dimension of the largest particles present controls their breaking strength. This dimension is approximately equal to F. As a starting point for the analysis it is assumed that a 1-in. steel ball will effectively grind material with 80 pct passing 1 mm, or with F- 1000µ or about 16 mesh. The breaking force exerted by a ball varies with its weight, or as the cube of its diameter R. The force in pounds per square inch required to break a particle varies as its cross-sectional area, or as its diameter squared. It follows that when a 1-in. ball breaks a 1-mm particle, a 2-in. ball will break a 4-mm particle, and a 3-in. ball a 9-mm particle. This is in accordance with practical experience, as well as being theoretically correct. Confirmation of this reasoning is supplied by the Third Theory of Comminution," which states that the work necessary to break a particle of diameter F varies as F. Since work equals force times distance, and the distance of deformation before breaka4e varies as F it follolvs that the breaking force should vary as F½ These relationships are expressed in Table I, with a 1-in. ball representing one unit of force and breaking a 1-mm particle. This establishes theoretically the general rule used in Eq. 1 that the ball size should vary as the square root of the particle size to be broken. Ball Size and Work Index: The work input W required per ton" varies as the work index Wi, and the

Jan 1, 1959