Search Documents

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

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

Jan 1, 1970

By C. G. Dunn

James, Jones, and Leak (JJL) conclude that growth-rate and orientation-distribution data obtained in a steep thermal gradient should be used with caution to account for isothermally produced recrystallization textures. I agree. However, it is not clear whether they conclude a) that the orientation dependence of rate of growth of secondaries in (commercial) silicon iron depends on both temperature and a temperature gradient or b) that, in applying rate data obtained in a steep thermal gradient, the isothermal annealing temperature must be close to the temperature for which the rate data apply. Also, it is not clear what tem-

Jan 1, 1967

By . M. Muskat

A THEORETICAL consideration of the mechanics of various water-flooding systems with the viewpoint of comparing their relative merits involves the study of two fundamental properties of such systems. On the one hand, it is important to determine the efficiency of the flood as measured by the percentage of the total area that is swept out by the encroaching water before water actually reaches the output wells. This is really a measure of the regional tendency of the water to "bypass" the oil. On the other hand, the various systems must be compared from the viewpoint of their resistances which determine the rate of fluid production obtainable for given pressures maintained at the water input wells. The first of these questions has been studied and very definite answers have been obtained by means of simple electrolytic models.1 These models showed that, in an ideal system, the percentages of the areas flooded out before water reaches the output wells are about 57 per cent for the direct line-drive flood (the input and output wells being equally spaced and placed in alternate parallel lines), 74 per cent for the regular five-spot network, and 79 per cent for the seven-spot network. With respect to the efficiency of the usual flooding arrangements, they are, therefore, given quite accurately by the above numbers. With regard to the resistances of the flooding networks, electrical conduction models have also been used successfully2 and are well adapted to the problem of determining experimentally the pressure distribution and streamlines, or over-all resistances of specific flooding networks. However, such models do not permit readily a study of the effect of changing dimensions or configuration within a network, since a new model must be made up for each set of dimensions considered.

Jan 1, 1933

By A. W. Ruff

The seismic wave produced by an explosive is very important in blasting. A true understanding of the wave is only important when considering possible structural damage to buildings located near the blast, but it is also of primary importance when considering rock fragmentation. Because of the importance of the seismic wave to fragmentation, a finda-mental study of the seismic wave and its effect was made at the Cananea Consolidated Copper Co.'s mine at Cananea, Sonora, Mexico. Strain gages were used to measure the strain in the rock when nearby explosive charges were detonated. The variance in the circuit voltage caused by changes in the resistance of the strain gage was measured with an oscilloscope and recorded with a Polaroid camera. CONFINED DRILLHOLE CHARGES To record the seismic wave from a confined drillhole explosive charge, a gage was cemented to the bottom of a 4-in diam drillhole. The horizontal gage hole and nine horizontal charge holes were drilled perpendicular to the center line of the 2-39 drift and 4 ft above the floor. The gage hole was 29 in. deep, and the charge holes were drilled so that each bottomed 1 in. deeper than the gage hole. The gage was installed with its axis horizontal and parallel to the 2-39 drift. The charge holes had the following horizontal distances from the gage: hole 3, 3 ft; hole 4, 4 ft; hole 5, 5 ft; hole 9, 9 ft; hole 10, 10 ft; hole 14, 14 ft; and hole 15, 15 ft. Figs. 1 and 2 show the relationships of the gage and the charge holes. Charges of 80 pct special gelatin, which weighed 50, 100, and 150 g were exploded in the charge holes and the seismic waves were recorded. All charges were stemmed with 200 g of mud and detonated with an electric blasting cap. Fig. 3 shows the recorded seismic waves. The reflected components of the waves were too weak to be objectional on these records. Table I shows the data taken from the wave forms. The approximate average velocity of a wave with a 9-ft travel distance was 12,900 fps. Since different charge weights were used, it was necessary to scale the data. The scaled maximum compress ion-strain data and the scaled maximum tens ion-strain data are shown in Figs. 4 and 5, respectively. Figs. 6 and 7 show the strain-distance data from Figs. 4 and 5 converted to a 1-lb charge. The results showed that the wave consisted of an initial compression pulse followed by a tension pulse. The amplitude of the tension pulse approached the amplitude of the compression pulse as the distance between the gage and the shot increased. The decay rate of the seismic wave is very important in the evaluation of the blasting characteristics of a rock. When high-amplitude stresses are produced, as in rock fragmentation, the rate of decay is much more rapid than when small amplitude stresses are produced. In underground mine blasting, the burdens are usually not greater than 6 ft and the decay is very rapid. The rapid seismic-wave decay rate in the fragmentation range means that in practice, if the minimum weight of powder is used to break a specified burden, only a very

Jan 1, 1961

By A. J. Terrones

THE classical papers on contact metasomatic deposits by Lindgren on the Clifton-Morenci district,' by Barrell on Marysville, Mont.,' and by Gold-schmitt on the Oslo district, Norway,' laid the foundation for understanding of the importance of separating the metamorphic process from the subsequent metasomatic stage. Knopf' made a big step forward by recognizing that sulfide ore minerals are invariably younger than silicates and by further stressing the importance of structural control" in ore finding. Some years later Harrison Schmitt" fought the idea that contact metasomatic deposits are erratic, irregular, bunchy, economically marginal, and of small size, concluding that these descriptions have little factual basis. He emphasized the importance of detail mapping and careful geological study as fundamental requirements for a better understanding of these deposits. This article presents a comparative description of the structural features at the Magistral, Antamina, Morococha, Vilca, Ferrobamba, Katanga (Quibio), Tintaya, and Huarca contact deposits (Fig. 1). Further, it will describe how these structural features controlled the emplacement of the intrusives, their shape, distribution of metamorphic and metasomatic halos, and sulfide deposition in post-skarn structural channelways. Finally, the modifications introduced by a superimposed hydrothermal stage recognized at some localities will be discussed. Primary exploration and careful sampling of recognized structural targets were instrumental in outlining the copper ore potential at the Antamina and Tintaya districts, thus proving the benefits and reward of careful geological examinations in this type of deposit. The Magistral District: This district lies in the Province of Pallasca, Department of Ancash, in northern Peru, at 8o 16' S and 77" 47' W. High to medium grade copper orebodies in garnetite, bordering a quartz monzonite stock, are exposed at the bottom of a glaciated valley in the western slope of the Andes at 4250 m (14,000 ft)* elevation.

Jan 1, 1959

By R. F. Hinds

A pressurizing system was designed and built to apply a radial pressure of 5.000 psi to rock samples. Samples of the Bradford, Weir and Kirkwood sandstones were subjected to radial pressures parallel and perpendicular to the bedding in the system and the changes in conductivity, porosity and permeability were deter-mined. Asymptotic decreases in conductivity, porosity and Klinkenberg permeability were observed over the 0- to 3,500-psi range. From these data, the increase of formation factor over the same range was calculated. The cementation exponent expressed in the Archie equation was found to increase with pressu,re. The rate of change varied from formation to formation. Only small differences were observed for samples taken parallel to the bedding as comnared to those taken perpendicular to it except for the Bradford sand. Above 3,500 psi the data are not good enough to tell if the properties continue to fall-off asymptotically or in a different manner. A possible explanation for the probable litnit of the asymptotic rate around 3,500 to 4,500 psi is that this represents the maximum pressure which the formations have unrdergone during their geological life as a result of burial. INTRODUCTION For many years measurements have been made in the laboratory of the properties of oil-bearing sandstones with little consideration being given to whether these properties had been changed when the samples were moved from their original position in the ground to the laboratory. Relatively little has been reported in the literature on the effect of pressure on the properties of sandstones. The changes in permeability for eight typical consolidated oil-bearing sandstone samples with overburden pressure (0- to 15,000-psi range) was presented by Fatt and Davis' in graphical form. In 1953 Fatt' also reported the relationship between porosity and overburden pressure for four typical consolidated oil-bearing sandstone samples (0- to 5,000-psi range). Fatt later described the changes of formation factor and porosity of 20 typical oil-bearing sandstone samples with overburden pressure changes (0- to 5,000-psi range) giving the values of the coefficients and exponents in the Archie' equation, and the more general Winsauer equation', where C, k and m are empirical constants, F is the formation factor, and + is the porosity. During the experiments of this investigation the core was contained in a neoprene rubber sleeve as shown in Fig. 1. It was held between two steel end members which were mounted in a steel cylinder. The unassembled pressure cylinder with its component parts is shown in Fig. 2. The core was insulated from the two steel end members by a Teflon sleeve with a small hole in the

By J. C. Osmond

One of the least known geologic regions in the U.S. is the area now called the Basin and Range Province. It is paradoxical that so little geologic information has been compiled for a province that has yielded billions of dollars worth of metals and is currently producing a small amount of oil. Regional stratigraphic studies of parts of this large territory have been published only recently, and regional structural interpretations of even parts of the province are rare in the literature. This, of course, does great injustice to the earliest workers: Gilbert,1,2 Russell,9 Le Conte,4 Spurr,5 Louderback,8 and others, who recognized the area as distinctively different from the rest of North America. Their simplified picture of the structure of the ranges and the province as a whole is the groundwork for present understanding and has never been superseded. Since their early works, however, geologic data from the province has consisted largely of detailed reports of widely scattered mining districts. Each district has been described as though it were an island in an unknown sea. More recent syntheses of the geologic history have been prepared by Lindgren,7 Ransome,8 Nolan,9 Bill-ingsley and Locke,30 Longwell,11 and Eardley.12 The tentative nature of the present interpretation, which treats the structural evolution of the Utah and Nevada area, must be considered in light of the scarcity of geologic information from the province as a whole. The methods of interpreting tectonic history from stratigraphy are based largely on the concepts of Marshal Kay.23 Discussions with Hoover Mackin also have aided in organizing the Tertiary history. TECTONIC AND STRATIGRAPHIC SETTING Two approaches to the subject of the tectonic history are possible; they correspond to looking through opposite ends of a pair of binoculars. One involves an expansion and extension of detailed work and the other an interpretation based on regional geological relationships. It is believed that the second approach, used here, is not in disagreement with detailed observations. The area of the Basin and Range Province is shown in Fig. 1. Within this area middle Tertiary to Recent high-angle faulting has controlled development of mountain ranges and valleys. Only the part of the area north of Las Vegas, Nev., and south of the Snake River in Idaho will be considered in this article. Within this part of the province most of the ranges trend north-south. Approximately two thirds of the area is covered by late Tertiary to Recent valley fill and lake deposits.

Jan 1, 1961

By David Cole

THE appearance of the Wilfley table in 1897 marked an epoch in the art of concentration of ores. The table has merited and received an almost unprecedented measure of public approval, lasting through its whole patent life. It has been very little improved in itself, or improved upon by competitors. The new machines bidding for popular favor have been of the Wilfley family, but without ability to do markedly better work than the parent m:-chine, or to (to more work except by the use of superimposed decks, the latter having obvious disadvantages. FIG. 1.-The- Butchart RIFFLE System. During the past three years there has been developed at Morenci a new type or arrangement of riffles applicable to the Wilfley type of concentrating table which corrects many of the objectionable features or limitations of the older system and obviates most of the difficulties encountered with it, particularly when handling ores having a low ratio of concentration. The new type is known as the Butchart1 riffle system, and its general arrangement is shown in Figs. 1 and 2.

Jan 2, 1915

CONTENTS PAGE ANDERSON, ARVID E., and CAMERON, FRANK K.-Recovery of Copper by Leaching, Ohio Copper Co. of Utah. Discussed by Thomas B. Brighton, H. C. Goodrich, Arvid E. Anderson, Thos. P. Billings, Frank K. Cameron 2 BEESON, J. J.-Mining Districts and Their Relation to Structural Geology. Discussed by Charles Keyes, Chas. A. Porter 4 DOBBEL, CHAS. A.-Deep-hole Prospecting at the Chief Consolidated Mines. Discussed by G. W. Crane 10 CRANE, G. W.-Notes on the Geology of East Tintic. Discussed by G. W. Crane, H. E. Havernor, N. O. Lawton 11 SACKETT, B. L., BARDWELL, CARLOS, JACOBSON, SIMON, and JENSEN, N. H. Lead Smelting in Utah. Discussed by D. W. Jessup, B. L. Sackett, E. H. Laws 13 BARDWELL, E. S., and MILLER, Roy H.-Pulverized Coal as Fuel for Copper Refining Furnaces. Discussed by L. D. Anderson, E. S. Bardwell, C. P. Crawford 15 WIGTON, G. H.-The Chief Consolidated Volatilization Process and Mill. Discussed by Thomas Varley, S. G. Olmstead, V. O. Lawton, O. C. Ralston.... 16 JESSUP, D. W.-Evolution in the Preparation of Ores for Lead Blast Furnaces. Discussed by W. Spencer Reid, D. W. Jessup, Carlos Bardwell 19 ALLEN, H. P., and MADGE, WM. C.-Chloridizing Mill of the Standard Reduction Co. Discussed by Thomas Varley, H. R. Hanley, H. P. Allen, George H. Dern 20 TUCKER, E. L., and HEAD, R. E.-Effect of Cyanogen Compounds on Floatability of Pure Sulfide Minerals. Discussed by Blarney Stevens 24 Mix, B. E., and BARKER, L. M.-Method of Unloading Ores and Coarse-crushing Practice at Magna Plant of Utah Copper Co. Discussed by W. C. Madge,B. E. Mix 25 KEIGLEY, C. T.-The New Byproduct Coke Plant of the Columbia Steel Corporation. Discussed by W. M. Grant, C. T. Keigley, A. C. Feilder 25 JACKSON, P. W.-Blast Furnace Plant of the Columbia Steel Corporation. Discussed by W. R. Phibbs, H. M. Henton 26 OLDRIGHT, G. L., KEYES, H. E., and WARTMAN, F. S.-Production of Ferric Sulfate and Sulfuric Acid from Roaster Gas. Discussed by G. D. Van Arsdale William E. Greenawalt, G. L. Oldright 27 MAZANY, M. S.-Braden Copper Company Caletones Smelter. Discussed by A. B. Young, G. L. Oldright 36 HARRINGTON, D.-Safety Methods in Utah Coal Mines. Discussed by Edward K. Judd, Otto Herres, Benedict Shubart, A. C. Watts, D. Harrington 36

Jan 12, 1925

By T. Kunitake, H. W. Paxton

The self-diffusion coefficient of chromium in various alloys in the iron-chromium system has been measured. A variation in Dofrom 10-4 for pure chromium to a maximum of 102 near 60 pct Cr appears with a range of activation energy from about 50 to 80 kcal per mol. The general pattern of behavior observed is not consistent with present theories of diffusion by a vacancy mechanism. Suggestions are made for further clarifying experiments. RECENT work on the self-diffusion of chromium,' ?-uranium2-4 and ß-zirconium,23, 26 has indicated that these elements have very low Do (10-3 to 10-4) compared to that conventionally expected for vacancy diffusion. Pound et al.5 have suggested that this may be due to ring diffusion being operative, and have proposed a model which gives Do of the correct order of magnitude. The essential difference between this model and Zener's earlier model6 is that whereas Zener assumes tight coupling between atoms in a ring to establish a vibrational frequency, Pound's model allows the atoms to vibrate essen- tially each in their own normal modes. The probability of each atom moving in the correct direction to cause a ring of n atoms to rotate by 2p/n must then be calculated as a "synchronal entropy" which may be -20 e.u. for a 4 atom ring. In iron, the latest determinations of self-diffusion yield Do in the range 18 to 522.7-9 The most accurate is probably that of Borg and Birchenall8 which is 118 ± 3. This value is not inconsistent with that expected for a vacancy mechanism. Since iron and chromium are very closely similar in size, are both transition metals, and their solution is thermodynamically ideal at temperatures of interest in diffusion, it appeared that a determination of self-diffusion coefficients, Do and Q as a function of composition would be of substantial interest to diffusion theory. The following measurements were made with this purpose in mind.

Jan 1, 1961

By L. F. Elkins

In 1946 Littlefield, Gray and Godbold published a thorough geologic description of the West Edmond Hunton Lime reservoir, located in Central Oklahoma, and discussion of its early performance.l They stressed implications of the nature of the reservoir rock on performance to be expected with production both by natural means and by pressure maintenance with gas injection. The essence of their analysis was that presence of an extensive interconnected system of fractures and solution channels constituting about 10 percent of the total reservoir void space in otherwise very low permeability reservoir rock would result in severe channeling of naturally encroaching water or injected gas through the fractures with little or no benefit to ultimate recovery of oil. Considerations were in progress to unitize the reservoir for pressure maintenance by gas injection. Some other engineers involved in the project had considerably different opinions at the time regarding continuity of fractures and regarding effects of such fractures on performance. The field was unitized Oct. 1, 1947. Subsequently, pilot gas injection tests were conducted in four areas of the field; large scale waterflooding was conducted in two parts of the field in the Bois d'Arc or Upper Hunton; and pilot waterflooding was conducted in the Chimney Hill or Lower Hunton. Full scale water-flooding is now in progress in the Chimney Hill section where it is oil productive. The primary purpose of this paper is to divine some of the secrets of this reservoir through analysis of all aspects of its performance, with the hope that it will be of help to engineers and operators in developing and operating other low permeability fractured reservoirs. In part, this considers variance of actual performance from some widely used conventional reservoir engineering assumptions. For brevity, other aspects of unit operation and its benefits to the owners are not discussed. Oil and Gas in Place For unit participation, net pay in the Bois #Arc defined by the gross SP anomaly on electric logs totalled 2,100,000 acre-ft in the 30,000-acre reservoir. Development of West Edmond predated modem logging techniques that now would be used to determine "net" or "productive" pay. In Fig. 1 this log pick is compared with detailed core analysis and core descripfion by Littlefield et a1.l for one well identified as A on the index map in Fig. 2. Not all of the gross Bois #Arc section had oil staining, but all cores with measureable intergranular permeability were stained. (The committee selection of the Bois d'Arc SP interval also includes the Frisco section as identified from cores by Littlefield.) The Bois d'Arc is oil productive throughout the entire 600-ft monoclinal accumulation from the eastern upstructure erosional pinchout to the west-side oil-water contact. (See Littlefield et al.l and McGee and Jenkins for geologic detail.)

Jan 1, 1970

By G. Thodos, D. E. Matschke

Cansiderable time and effort frequently are expended to establish, with a degree of confidence, the PVT behavior of pure substances. In particular, a great deal of experimental information contributed by a number of investigators is available in the literature. The data of each investigator, although significant in themselves, cannot be used with a high degree of confidence until they are compared with the continuum of information presented by others. Current interest in thermodynamic and transport properties, for not only the gaseous state but the con-densed state as well, requires that PVT infarmation be presented for both regions. The density dependence on temperature and pressure, when presented as a continuum, is of considerable value in itself, especially since the available equations of state fail to properly account for the PVT behavior in the condensed state. For this reason and to incorporate PVT data recently presented in the literature, this investigation was carried out. In their studies concerned with the thermodynamic properties of methane, Matthews and Hurd in 1946 utilized the data of Kvalnes and Gaddy and Olds, Reamer, Sage and Lacey to establish the PVT behavim of methane in the superheated region and the saturated liquid densities of Keyes, Taylor and Smithl in conjunction with the Clapeyron equation to develop the densities for the saturated envelope. Current activity sponsored by API Project No. 44 under the direction of F. D. Rossini has led to a tabular presentation of compressibility factors for gaseous methane in the superheated region.la These values were obtained by using smoothed residuals, developed from experimental data and the Ben edict-Webb-Rubin equation,' to establish internally consistent compressibility factors for the real gas. This approach offers probably the most precise presentation of PVT information for gaseous methane to date. The PVT behavior of methane also can be represented accurately by equations of state such as the Benedict-Webb-Rubin equation or the Martin-Hoa equation.' These equations, although valuable for analytically representing the PVT behavior, lack precision in the compressed-liquid region. The PVT data available in the literature for methane were utilized to obtain reduced densities that are presented in Figs. 1 and 2 as functions of reduced temperature with parameters of reduced pressure. The critical constants utilized for methane are those reported by Kobe and Lynn." These critical values, T, = 191.1°K, PC = 45.8 atm and p, = 0.162 gm/cc, produce a critical compressibility factor, z, = 0.289. The PVT information, as obtained from the literature, in many instances was not presented in an immediately

By Albert Zima

ACCORDING to recently published statistics, it is predicted that as much oil must be produced during the next 16 years as has been produced during the past 75, in order to satisfy the high rate of consumption of petroleum products. During 1934, over 17,000 oil wells were completed in the United States, and it is expected that over 20,000 additional wells will be drilled during the current year. Were it not for the remarkable progress made during recent years in the improvement of drilling technique and in the development of equipment used in the production phase of the oil industry, it is doubtful whether sufficient oil could be produced to meet present requirements. A few years ago, wells seldom were deeper than 3500 to 4000 ft., owing to the limitations of the drilling methods and equipment employed. Today, however, through the application of scientific drilling methods and the adaptation of modern equipment and improved materials, many wells range from 6000 to 8000 ft. in depth, and a few recent completions exceed 10,000 feet.

Jan 1, 1935

By Leo F. Reinartz

It was with humility that I accepted the invitation of the Board of Directors of our Institute to deliver the Twentieth Howe Memorial Lecture. Many previous lecturers could speak from personal experience of the great work done by Professor Howe. It was not my privilege to be one of his students but, indirectly, I came under the influence of his work many years ago. My mentor in Metallurgy at Carnegie Institute of Technology, Prof. Fred Crab tree, was an ardent admirer of Professor Howe and his work. Professor Crabtree taught me early to appreciate Professor Howe's love of truth and keen desire to find a scientific reason for every metallurgical phenomenon. In this early search, Professor Howe encountered many unexplainable facts. Nevertheless, in the classic text, "The Metallurgy of Steel," published in 1891, he indicated that he believed Research could light the way. Under the heading of "The Treachery of Steel" he said: "In the early use of Bessemer and Open Hearth Steel, many then unexplained and hence mysterious failures occurred. With our present knowledge, an easy explanation of most of them would probably have been seen. But even today certain unexplained and apparently inexpli- cable failures occur in steel; inexplicable in spite of full and intelligent investigation." Then Professor Howe continued to talk about "the trustworthiness of steel" in the following manner: First, the steel-maker has learned by experience. He knows today far better the effects of a high percentage of carbon or of phosphorus; of cracks, pipes and blow-holes; of segregation and imperfect mixing; of over-heating; of finishing too hot or too cold. Knowing, he guards against them more effectually, and keeps at home much steel that he would formerly have sent into the market. Those who would not, or could not learn and do, have been driven out of the business; the conditions necessary for producing good sound steel by the acid Bessemer and the acid Open Hearth processes have been mastered. By and by, the basic process came along with new conditions, new liabilities to unsoundness, a great hue and cry about mysterious fractures followed, and Lloyd's register provisionally forbade the use of basic steel. Much basic steel was irregular and brittle throughout; too much carbon, too much phosphorus, too cold teemings—imperfect mixing. Still there was little doubt that, with further experience, these difficulties of the basic process would be mastered as those of the acid process had been. This shows that Professor Howe could not take anything for granted. Believing that careful research would find the answer to all mysterious phenomena, he was several generations ahead of the practical acceptance and adoption of this principle.

Jan 1, 1943

By Leo F. Reinartz

It was with humility that I accepted the invitation of the Board of Directors of our Institute to deliver the Twentieth Howe Memorial Lecture. Many previous lecturers could speak from personal experience of the great work done by Professor Howe. It was not my privilege to be one of his students but, indirectly, I came under the influence of his work many years ago. My mentor in Metallurgy at Carnegie Institute of Technology, Prof. Fred Crab tree, was an ardent admirer of Professor Howe and his work. Professor Crabtree taught me early to appreciate Professor Howe's love of truth and keen desire to find a scientific reason for every metallurgical phenomenon. In this early search, Professor Howe encountered many unexplainable facts. Nevertheless, in the classic text, "The Metallurgy of Steel," published in 1891, he indicated that he believed Research could light the way. Under the heading of "The Treachery of Steel" he said: "In the early use of Bessemer and Open Hearth Steel, many then unexplained and hence mysterious failures occurred. With our present knowledge, an easy explanation of most of them would probably have been seen. But even today certain unexplained and apparently inexpli- cable failures occur in steel; inexplicable in spite of full and intelligent investigation." Then Professor Howe continued to talk about "the trustworthiness of steel" in the following manner: First, the steel-maker has learned by experience. He knows today far better the effects of a high percentage of carbon or of phosphorus; of cracks, pipes and blow-holes; of segregation and imperfect mixing; of over-heating; of finishing too hot or too cold. Knowing, he guards against them more effectually, and keeps at home much steel that he would formerly have sent into the market. Those who would not, or could not learn and do, have been driven out of the business; the conditions necessary for producing good sound steel by the acid Bessemer and the acid Open Hearth processes have been mastered. By and by, the basic process came along with new conditions, new liabilities to unsoundness, a great hue and cry about mysterious fractures followed, and Lloyd's register provisionally forbade the use of basic steel. Much basic steel was irregular and brittle throughout; too much carbon, too much phosphorus, too cold teemings—imperfect mixing. Still there was little doubt that, with further experience, these difficulties of the basic process would be mastered as those of the acid process had been. This shows that Professor Howe could not take anything for granted. Believing that careful research would find the answer to all mysterious phenomena, he was several generations ahead of the practical acceptance and adoption of this principle.

Jan 1, 1943

Jan 1, 1969

By J. H. Swisher

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

Jan 1, 1964

By D. F. Malott

The Michigan Dept. of State Highways makes extensive use of geophysics for subsurface surveying which would be applicable for uses in other fields. Examples of resistivity surveys are given which include a proposed large-volume excavation with suggestions for land use after excavation, subsurface evaluation for engineering design, and a subsurface survey for land appraisal. An example of a combined refraction seismic and resistivity survey is shown where the resistivity is used for outlining the unconsolidated soils, while the seismic method is used to outline bedrock and to evaluate it for rippability. A refraction seismic survey at the bottom of bridge caissons for evaluation of bedrock uniformity and a rehtive aid in assessing rock strength is reported. These past several years have been a period of unprecendented construction and commercial expansion. Vast private manufacturing complexes and government road, power, and office projects are being built or are in planning stages. Nearly everywhere, municipalities are expanding outward and upward. Many towns are developing into cities and many cities are displaying new skylines, and the building boom shows little sign of slowing down.

Jan 1, 1970

By H. K. VAN POOLLEEN, W. Hurst, H. C. Bixel

In an earlier publication* I showed the development of the instantaneous point source solution for a well producing at a constant rate at the center of a system of two radial, adjoining sands of different permeabilities. I also showed how a wellbore radius could be incorporated if this were the development required. For the purpose of that paper and the point I wished to make, this solution served. It yielded a simplifying solution of sands in series identified by Eqs. VI-26 and VI-27. (These equations are reproduced in part by Bixel and van Poollen as their Eq. 3, Page 302.) I am well aware that the simplified solutions are approximations. The facts that these Ei-functions satisfy the diffusivity equation and meet the boundary conditions both at the well and at the interface between sands, and that the earliest time considerations yield zero pressure drop, are evidence that the solutions are good approximations. This was the purpose for these formulas, which I interpret that where the permeabilities are of finite order of magnitudes in the respective media, the relations give satisfactory answers. However, from the authors' point of view there are reservations, which leads me to consider their entire generalization as illustrated in their Figs. 1 through 14. What I understand from their storage capacity, relationship Fs, which is the ratio of the product of porosity and compressibility of Sand II to Sand I, that is employed here by the authors as a correlative means of identifying one figure from another, there is nothing in the exact mathematics that supports such a term as a means of correlation. This is taken out of context, unsustained that this applies. Here the mathematics are known, where the true correlation is the diffusivity term of one sand to another. This entails permeability that refers back again to the term M in their paper for the transmis-sibility of fluid from one sand to the other that should have also been a part of this storage capacity. In Figs. 8 through 14, which purported to show the simplified equation following this correlation, there are apparent discontinuities for M equal to or less than 0.1000 that differ from the trends for greater values of M. In this connection there is no reference radius given for this cumulative pressure drop in the paper that would permit me to make these calculations. However, what is clear for M equal to 0.1000 or less we are treated to almost a closed barrier at the interface; a condition entered by the writers without qualification. This across-the-board treatment of these simplified formulas, whatever the motive of the authors, enters into a range of permeabilities that even I as the contributor would have reservations that such could apply. Under comparable conditions I would be looking at a closed reservoir or some facsimile for such. The extent to which the approximation or the exact equations pertain is not evident in the generality offered by Bixel and van Poollen in their paper. Only when they or other writers show the comparison of the pressure distributions within the formations can one be able to judge the rational limits to subscribe to the simplified formulas.

Jan 1, 1969

By K. E. Brown, J. P. Brill, A. R. Hagedorn, T. C. Doerr

A recent correlation for multiphase vertical flow by Hagedorn and Brown1,7 has been used to examine the effect of numerous variables such as pipe size, gas-liquid ratio, liquid flow rate, liquid viscosity, liquid surface tension, oil API gravity, water cut, etc., on flowing pressure gradients. Most previous correlations were derived by using a form of the general energy equation which assumed negligible changes in kinetic energy. Under certain flow conditions this can lead to serious errors in predicting pressure traverses. Examples of traverses calculated both with and without the kinetic energy term are included here. The Chew and Connully correlation for calculating changes in viscosity with changing pressure and temperature has been extended so that the effect of a changing viscosity with depth on pressure gradients can also be predicted. Examples illustrating the effect of tubing size show the importance of this correlation in determining pipe sizes for initial well completions. Limitations on production rates in small pipe sizes from excessive friction losses can be predicted with greater accuracy. The combined use of a vertical-flow correlation with a horizontal-flow correlation provides a method for determining the maximum flow rate possible from both flowing and gas-lift wells. Examples are included to illustrate application of the combined use of vertical and horizontal-flow correlations in optimizing production rates and minimizing excessive costs from over injection of gas in gas-lift wells. INTRODUCTION The problems of multiphase vertical and horizontal flow are immediately encountered in any producing well. Fig. 1' shows three stages of production: flow in porous medium, vertical flow and horizontal flow. All three stages affect the production rate from a well. If we are confronted with a particular well and consider only the variables involved in the vertical and horizontal stages of Fig. 1, we should be capable of calculating the flowing bottom-hole pressure necessary to produce a particular flow rate. A change in any of the variables results in a new solution to the problem. Recent approaches in developing horizontal and vertical-flow correlations have improved to the extent that the solution can now be obtained with more confidence than previously. VERTICAL AND HORIZONTAL FLOW CORRELATIONS Since Poettmann and Carpenter's' original work in multiphase vertical flow, there have been many additional vertical-flow correlations developed which have greatly improved the accuracy of pressure-loss calculations. These

Jan 1, 1967