Search Documents

By R. P. Harvey, J. E. Walstrom, A. A. Brown

In view of the large number of wells directionally drilled from offshore platforms and urban drillsites, greater interest is being focused on the directional survey and on the uncertainty inherent in calculating the position of the bottom of the hole, or other benchmarks in the wellbore. Of the two methods presented here for evaluating this uncertainty, the first is analytical and uses the principles of probability and statistics. The second method uses Monte Carlo simulation to generate a pattern of the scatter of possible bottom-hole locations. Both methods are based on a model for which the usual assumption is made that the curvilinear axis of the wellbore is adequately represented by a large number, M, of linear segments joined end to end from the top to the bottom of the hole and closely following the course of the hole. The errors due to this approximation are assumed to be negligible and are not studied here. One form of the model is to consider that for a given hole there exists a perfect set of data describing it, of the above form. Each linear segment has an angle of inclination with the vertical 7, in radians, an azimuthal angle from true north xi in radians, and a length Si, in feet. The subscript variable, i, extends from 1 to M. Each of the M linear segments representing the axis of the wellbore has three components in the X, Y, and Z directions. These components are calculated from the equations: The bottom-hole coordinates as determined from the reference set of data for the wellbore are given by the following for i = 1 to M: X=SX, Y = Sy, Z=SZ; . (2) A directional survey instrument, lowered into the wellbore to the M segments, or stations, measures the values of the 7; and the 3,. In general, because of errors from various sources, the resulting angular measurements will not be precisely equal to the 7" and xi. Random variables l, and Ai are therefore considered with associated probability density functions fi(li) and gi(Ai), respectively, which represent the distributions of potential survey measurements for this hole. Since it is assumed that the values Si are accurately measured they are treated deterministically. The random variables corresponding to the surveyed bottom-hole coordinates are given, for i = 1 to M, by:

Jan 1, 1970

By Samuel Hoyt

SOME of the more important properties of finished materials are strength, ductility, toughness, resistance to alternating and repeated stresses, etc: Of these, the property that appears to have received the least consideration, at least in this country, is toughness, which is clue, possibly, to our regarding toughness as a qualitative property or, at any rate, as one that, parallels ductility. This arises out of our greater experience with the simple carbon steels in the rolled or forged (non-heat-treated) condition which, qualitatively, are tough if ductile and are not tough if brittle. In this way we have erroneously become, accustomed to judging the toughness of a material from its behavior in the tensile or similar test. The point of view presented in the present paper is that toughness, like hardness or tensile strength, should be regarded as an independent property and of sufficient importance to require, in so far as that may be possible, quantitative determinations. If such be true, it becomes at once necessary to devise experimental means for measuring, or valuating toughness, and the notched-bar impact test (the Charpy test) is advanced as the most logical test that has so far been developed for this purpose. It is further advanced, with particular emphasis, that we have two kinds of toughness to deal with and, accordingly, they will be dealt with independently. CLASSIFICATION OF MATERIALS ACCORDING TO TOUGHNESS Toughness has been defined somewhat as follows: Tough materials are those that offer considerable resistance to permanent deformation but which, once such resistance has been overcome, may' be deformed plastically, but only by the expenditure of considerable energy. In other words; tough materials may be deformed plastically but they absorb a considerable amount of work in the process. This kind of toughness may be called "static" toughness when the rate of loading is reasonably slow or "dynamic" toughness when the rate of loading is

Jan 2, 1919

By K. C. Mills, E. T. Turkdogan

A number of investigators1 6 have noted the presence of a liquid miscibility gap in the Fe-Sn binary system. However, the first attempt to measure the

Jan 1, 1964

By Sir Robert Hadfield

SINCE the writer has been intimately connected with the development of manganese steel for many years, some remarks upon the early work with regard to the rolling and forging of this material might he of interest, and particularly as to the question of rails of special quality and of superior wearing qualities. In a paper on Iron Alloys, with Special Deference to Manganese Steel, read before this Institute at the International Engineering Congress in Chicago in 1893,1 the writer mentioned that, about 20 years ago, his firm in Sheffield produced railroad material, such as tires and axles, of forged and rolled manganese steel, thus showing that this product could be manipulated and worked up into the various required forms. The tires illustrated in the paper referred to above formed a particularly difficult class of work to produce in a special steel, but no special obstacles were met with Although there was no difficulty at that time in making and producing manganese steel in various forms, forged, rolled or pressed; nevertheless, on account of its higher expense, sufficient encouragement

Jan 2, 1914

By Pervez Patel, Martin G. Frohberg, Demetrios Papamantellos

Model experiments were carried out on the mass transfer of iodine between n-hexane and a 4 2 X 10-3 M aqueous sodium hydroxide solution. The influence of mechanical stirring by means of small gas bubbles (ca. 0.01 ml in volume) on mass transfer was examined in independent experiments. The results were described with the help of the penetration theory and in the case of bubble stirring with the help of a mathematical model based on this theory. S LAG-METAL reactions, especially those in steel-making processes, often occur under conditions where gas is evolved in the form of bubbles that break through the slag-metal interface. It is known that stirring by means of bubbles increases the mass transfer coefficient, but there has been very little quantitative work done in this field. Because of the difficulties that attach to experimental work at high temperatures this study was confined to model experiments at room temperature. EXPERIMENTAL APPARATUS The transport of iodine from an n-hexane layer into a weak sodium hydroxide solution was studied a) with stirring by means of bubbles and b) with mechanical stirring. The apparatus for the experiments with bubbles consisted of the cell, two Dreschel bottles filled with hexane in order to saturate the nitrogen flow, one flow meter, and a nitrogen gas cylinder fitted with a fine-regulating valve. Fig. 1 shows the cell consisting of a glass cylinder 33 cm in height with an overall diameter of 7.5 cm embedded in a copper flange. The flange is connected to a Plexiglas plate and to the gas inlet. The Plexiglas plate has holes of 0.2 mm diam bored into it which are spaced 1 cm from each other. For the experiments with mechanical stirring the bottom of the cell was sealed off by a Plexiglas plate. The upper layer was stirred by means of a glass stir-rer and the bottom layer with the help of a magnetic stirrer in a counter rotatory direction. Both stirrer speeds could be varied continuously. EXPERIMENTAL PROCEDURE Iodine dissolves in n-hexane producing a stable violet color thus enabling the iodine concentration in

Jan 1, 1970

By H. L. Miaw, A. M. Gaudin, D. W. Fuerstenau

A quantitative method for evaluating density of slime coatings has been developed and applied to formation of iron oxide slime coatings on quartz and on corundum. Slime coating density is related to flotation recovery and to properties of the electrical double layer at mineral surfaces. IN spite of considerable study, the nature of slime coatings in flotation is still not completely understood. However, phenomena that control floc-culation and dispersion of colloidal systems are now interpreted in terms of electrical double layers.",' Colloids will flocculate even though the particles carry the same charge, but when two dispersed colloids of opposite charge are mixed. flocculation ensues even more rapidly. Analogous phenomena seem responsible for slime coating on minerals in flotation systems. To shed new light on the issue it was considered desirable to devise a tool that would measure slime coatings. The usual method has been to examine photographs of presumably representative surfaces, after. an intervening washing in running water. But such an examination is non-quantitative, and the changes in chemical environment represented by the washing may have removed the slime coatings that existed during flotation. To provide the possibility for a correlation with electrical double layer phenomena, the experiments were made with quartz, corundum, and iron oxide, for which elec-trokinetic data are available.'-" Since the properties of the electrical double layers of these oxides depend on the pH of the solution, pH will be the most important variable. Furthermore, because the zero point of charge for silica, iron oxide, and alumina occurs at pH 3.7, 8, and 9.5, respectively, the relative charge on the slime and mineral can be changed over wide ranges. Because of its practical importance in quartz flotation, iron oxide powder was chosen as the slime. Also, its uniform composition, contrasting color, and ease of analysis make experimentation simple. In this investigation, correlation was sought between the effect on flotation recovery of slime addition and slime coating density. Materials and Methods Details are available in a thesis by Miaw." The quartz was selected from pure Brazilian crystalline quartz, crushed, sized (65/150 mesh) acid-cleaned, washed, and stored in conductivity water. The corundum (synthetic alumina of high purity prepared by Linde Co., Chicago) was crushed, sized (65/100 mesh) deslimed, boiled in nitric acid, washed, and stored in conductivity water. The iron oxide slime, reagent-grade ferric oxide (—400 mesh) from Allied Chemical & Dye Corp., was heated to 160°C, cooled, and stored in a vacuum desiccator. In the electron microscope the slime particles appeared to be about 0.3µ in size, but their large specific surface (87,900 cm per g, by gas-adsorption method, using krypton) suggests that the actual particles are smaller than 0.3r. Had they been perfect nonporous spheres, they would have had a diameter of 0.13µ.

Jan 1, 1959

By D. R. Parrish, F. F. Craig

In some respects, forward combustion (alias: fie-flooding, underground combustion, in-situ combustion, dry combustion) is an efficient oil recovery method. For example, its displacement efficiency (as much as 95 percent in reservoirs with high porosity and high oil saturation) is good. In its use of heat, however, forward combustion is grossly inefficient. Most of the heat generated during forward combustion is wasted. About one-fifth of the heat is picked up by the incoming air and re-used (regenerated); the remainder is left behind the combustion zone and ultimately is lost. Compressed air is expensive, and heat generated from compressed air via underground combustion is too valuable to be thrown away. Improvement is needed. Since the early 1950's, improvements in forward combustion by addition of water injection have been discussed, particularly in the patent literature, with increasing frequency.'-'' These concepts are based on the fact that relatively inexpensive water can be used to pick up much of the heat that is ordinarily wasted during forward combustion, and transport it forward to a more useful place, thereby reducing the amount of air required to move a heat front through the reservoir. Several approaches have been proposed. Merriam and Squires1 taught that a combustible (fuel) gas should be injected along with the air and water. This, however, would give at least a tendency toward reverse combustion instead of forward combustion, and would reduce the velocity of the heat bank, defeating the primary purpose of the water injection. Pelzer2 described primarily a batchwise method of forward combustion and water flooding. He advocated burning part way through a reservoir with forward combustion, then switching from air to water injection alone to push the heat toward the producing wells. Though possible, this method has serious disadvantages, the major one being that, by the time water injection is started, much of the heat has been irretrievably lost to the overlying and underlying formations. We have described what is considered a practical Combination Of Forward Combustion And Water-flooding (COFCAW).3,4 Air and water are injected simultaneously (or alternately) after a small heat bank has been formed by forward combustion. Some of the injected fluid fills up the region behind the heat bank; the rest enters the heat bank. Upon encountering heated rock, the water is converted into steam, which then flows ahead of the combustion zone and displaces much of the oil in place. Only enough air is injected to provide heat for vaporization of the water and to compensate for reservoir heat losses. Since the size of the heat bank is minimized, heat losses also are minimized. Because of more efficient oil displacement ahead of the combustion zone, fuel consumption is reduced. Too, unlike a forward combustion zone that cannot move until all fuel at a point has been consumed, COFCAW leaves some

Jan 1, 1970

By Charles H. Jacoby

The solution mining of evaporites has a history of approximately 5,000 years. The insitu leaching of metalliferous ore will have many of the same characteristic problems which have confronted industry for years in the production of evaporite minerals. Basically we will be converting the less soluble forms of valuable metals into their more soluble salts and producing these in much the same manner as, for example, sodium chloride. There are a number of problems which we are faced with that are more complex and complicated and, therefore, will require a higher degree of sophistication than previously employed. One of our most important undertakings will be the clearcut definition of the problems involved. While a great deal of expertise can be drawn from the solution mining of evaporites and the oil industry, problems peculiar to the hydrometallurgical field will need to be solved. One of these is the introduction of permeability into an impermeable rock so as to expose the sought after values of metals which by volume are infinitesimal when compared with sodium chloride, potash or other type water soluble mineralization. After establishing this permeability which will be miniscule when compared to a solution cavity, there will be the maintenance of this permeability. In general, salt deposits are of a somewhat uniform chemical composition and thus their chemistry was relatively simple. Each metalliferous ore body is somewhat unique in its chemistry and will need to be treated accordingly.

Jan 1, 1974

By L. F. Pease, J. H. Brophy

A phase diagram for the Ta-Zr system is presented. The system is of the minimum-melting point type with the 0-zirconium phase decomposing monotectoidally at 785°C and 95.5 at. pct Zr. The minimum solidus temperature is 1875°C and the temperature of the top of the monotectoid loop is 1775°C. Techniques employed included X-ray diffraction and fluorescmce, electron microprobe analyses, diffusion couples, metallography, and melting-point measurements. Errors due to interstitial solutes and in melting observations are particularly common in this system, and their occurrence is discussed. SEVERAL versions of the Ta-Zr phase diagram have been reported in the literature.lm7 Williams et al.' report a diagram of the minimum-melting point type with a monotectoid reaction at 800°C and 12 wt pct Ta. The minimum in the solidus occurs at 1820°C and 25 wt pct Ta, while the top of the monotectoid gap is located at 1780°C and 80 wt pct Ta. The experimental techniques employed were metallography, electrical resistance variation with temperature, dilatometry, X-ray diffraction, and direct observation of melting by the technique of Pirani and ~lterthum.' Experience in the Nb-Zr phase diagram has shown that small changes in the interstitial content can cause large shifts in the positions of the phase boundaries. Small amounts of inter -stitials can cause similar changes in the Ta-Zr system and may be responsible for the wide variations in the phase diagrams previously reported.'-' The authors15 have found that alloys containing 0.75 wt pct 0 and 1.0 wt pct Ni with tantalum contents from 19 to 85 at. pct often contain a lamellar or cellular precipitate after annealing above 1400°C. This precipitate is very similar to that reported by pollack,12 Ralls," and Donlevy17 for the Ni-Zr system. In the present investigation, a wider variety of experimental techniques and improved alloy purity were applied to the resolution of errors in diagrams for this system as they now appear in the literature. New results have been obtained in the areas of alloy homogeneity and purity. Errors in the common techniques of direct melting observation in binary alloys have been demonstrated. In the commercially important high-tantalum region, the more accurate metallographic method of solidus determination has been employed. EXPERIMENTAL PROCEDURE Nominal analyses of the materials employed are given in Table I. The zirconium was low hafnium reactor grade iodide crystal bar, and the tantalum was in the form of high-purity sheet. Powders were impossible to use because of gas content. Alloys were melted on a water-cooled copper hearth, with a thoriated tungsten-tipped nonconsumable electrode in 2/3 atm of Zr "gettered" argon. Alloys were

Jan 1, 1963

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 W. L. Saunders

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

Jan 2, 1914

By R. W. Parsons

The over-all apparent single-phase permeability of fracture-rock systems was studied using two different two-dimensional models. In a strict sense the results are applicable only to these models, yet may give some insight into real reseruoir behavior. The regular fracture-matrix model can have any number of regular fracture sets. All fractures within each set have the same orientation, width and spacing. The gross behavior is amlytically equivalent to that of an anisotropic permeable medium. In the heterogeneous fracture systems, individual fracture conductivities vary within a two-dinensioml network pattern. A computer program calculates the overall permeability after the conductivities are randomly placed in the network. The calculated permeability depends upon the distribution of these conductivities and the particular fracture pattern into which they are placed. Although detailed information on actual reseruoir fracture systems is not available, this study shows how some of the possible variables may affect obserued gross reservoir characteristics. INTRODUCTION Characteristics of reservoir rocks are determined by direct (core analysis) or indirect (logging and production tests) examination. Both methods, singly and together, are widely used in reservoir evaluations. For fractured reservoirs a problem exists in relating characteristics of cores to the in situ reservoir properties. Knowing fracture characteristics can be of value in determining production mechanisms. Because of the limitations of observation, the intimate details of a fractured rock system will never be known. One approach in analyzing fractures is to study the behavior of some conceptually simple models. The two models analyzed contain features which may be realistic, yet are mathematically tractable. Flow is considered in two dimensions only, with the plane of interest at right angles with the fractures. This concept represents the horizontal flow in a reservoir system with the fractures oriented vertically. Only single-phase, laminar flow is considered. The regular fracture-rock system permeability is calculable from a simple analytical expression. The heterogeneous system requires a computer program for calculating overall permeability.

Jan 1, 1967

By K. Leutwyler

The key to realistic casing stress analysis in thermal recovery installations is accurate knowledge of the temperatures involved. Much information leading to prediction of heat losses from tubing string to wellbore has been published. This .same information permits calculation of casing temperature under certain conditions. This paper reviews the basis for most of the published work, the line source solution to the diflusivity equation and examines some of its application criteria in the light of unsteady-state conditions. An alternate finite differences method is presented for correlation, and experimental data pertaining to the heat transfer mechanism in the annulus is reported. In the light of the information presented in the paper, it appears that the assumption of a line source for the tubing-wellbore system is valid even for reasonably short periods of heat flow, .such as "steam soak" operations. The empirically modified equation presented earlier matches computer derived results of the finite differences approach lor periods as short as 50 hours. Previously predicted cas ing temperatures and the associated stress levels are thus verified based on these correlative computations and the reported experimental information. Fall out from the computer output resulting in some studies of the influence of insulating cement properties on casing temperature, and experimentally studied effects of nitrogen placed in the annulus at various pressures., are described. INTRODUCTION Generally speaking, predictions of average casing temperatures have been based on an idealized model involving a centralized tubing string at uniform constant temperature radiating energy towards the casing under steady-state conditions.'. ' Heat is transferred away from the casing to the formation by unsteady-state conduction across thermal barriers consisting of cement and the mud sheath. Heat balances written on the casing permit prediction of instantaneous casing temperatures. During early periods of injection, casing temperature varies appreciably as a function of time; however, for long injection times this variation decays as heat flow from the wellbore to the formation stabilizes, and it becomes virtually negligible as conditions surrounding the wellbore become quasi-steady-state (Fig. 1). Carslaw and Jaeger -escribe a solution to the diffusivity equation for a line source (Fig. 2). A constant supply of heat is introduced into the medium from the line placed into the origin of the x-y coordinate system and parallel to the z-axis. Heat spreads outward in the infinite cylindrical solid and the temperature increase on a line going through point (x, y) parallel to the z-axis at a distance r can be predicted by: The following criterion must be considered: the Carslaw and Jaeger equation is obtained on the basis of a large

Jan 1, 1967

By E. S. Wheeler

The rapid increase in the use of monel metal and malleable nickel in the form of sheet, rod and tube has resulted in a similar increase in the demand for these metals in the form of castings. These castings range in size from small valves and pipe fittings weighing a few ounces to large centrifugal pumps, gate valves and marine propellers, some of the last weighing as much as 24,000 lb. A group of such castings is shown in Fig. 1. A large nickel agitator is shown in Fig. 2. Monel metal castings are widely used in accessories for steam power plants and in a variety of chemical plant equipment where the corrosion resistance, high strength and toughness characteristic of this metal render it particularly valuable. Nickel castings find extensive use in the construction of corrosion-resistant apparatus, and are often used in conjunction with wrought nickel sheet and tube in the fabrication of equipment for handling a variety of chemical and food products, particularly dairy products. In response to this demand for castings, numerous foundries have undertaken the production of castings of these metals and have met with varying degrees of success. The successful production of castings of monel metal and nickel requires careful attention to molding practice, proper melting equipment

Jan 1, 1930

Jan 1, 1970

Executive JOHN M. LOVEJOY, Chairman ERLE V. DAVELER G. B. WATERHOUSE WISER JUDSON WILLIAM WRAITH J. V. W. REYNDERS, Consultant Finance HENRY KRUM, Chairman PAUL D. MERICA R. M. ROOSEVELT H. G. MOULTON, Consultant Admissions. RALPH M. ROOSEVELT, Chairman GEORGE D. BARRON, Vice-chairman JOHN L. CHRISTI PHILIP KRAFT W. A. SINSHEIMER

Jan 1, 1936

By S. H. Ward

Since the advent of the first airborne electromagnetic system, it has been evident that such systems were inherently limited to shallow depths of exploration of the orderof 100 to 200 feet. Hence in 1948 the author and his associates commenced a study of the feasibility of applying natural audio frequency magnetic fields to investigation of the electrical properties of the subsurface in order to provide a system capable of overcoming this limitation. This research culminated, in the summer of 1958, in the development of a new airborne prospecting device — AFMAG. This device has detected surface massive sulphide bodies while flown at heights as great as 3000 feet. Hence, by inference, it is deduced that it could find similar sulphide bodies if they were buried 2500 feet below surface and the aircraft flown at 500 feet. Thus a considerable improvement over conventional airborne electro-magnetic systems is obtainable with AFMAG. Airborne exploration in mountainous terrain involving deep valley fill now becomes possible. Examples are presented of the results of actual airborne AFMAG test surveys over known deposits of massive sulphide mineralization. Interpretation of the data is discussed with particular reference to the knowledge of size, dip, depth, depth extent, and conductivity of causative bodies as deduced from the flight records. The possibility of mapping major geologic structures, for both mining and oil exploration purposes, is discussed with the aid of examples of airborne AFMAG data. INTRODUCTION Electrical prospecting with natural audio frequency electromagnetic waves as a source has, within the last three years, become a practical reality. Both airborne and ground survey instruments under the trademark "AFMAG" (for audio frequency magnetics) have been used successfully in the search for conductive mineral deposits. It is the object of this paper to discuss the principles of airborne AFMAG and to illustrate its application by means of examples obtained over known geological situations. WHY AFMAG? In view of the large number of good airborne electromagnetic systems available for purchase or contract surveying, one may well ask, "Why another new system?". Primarily the answer lies in the fact that AFMAG does not suffer from the very limited depth of exploration available with, and inherent in, any conventional airborne electromagnetic system. AFMAG is not conventional insofar as an artificial transmitter is not required. Only a receiving and detecting device is employed. The electromagnetic fields measured are those from millions of electrical discharges in the atmosphere (e.g;..- lightning strokes). Their random distribution in time and space make these fields difficult to employ and it is only with recent advances in the electronics of geophysical instrumentation that this has become possible. A second reason for AFMAG and one of its

Jan 1, 1961

By J. C. Martin

Stmplified equations are developed for the flow of fluids in gas drive reservoirs in which the effects of gravity can be neglected. The results show that the pressure distribution is governed by a nonlinear heat flow type of equation, and the saturation distributions are related to the pressure by the same equations QS those developed by Muskat for the average pressure and saturations in gas drive reservoirs. Unrler certain conditions the equation for the pressure can be approximated by the linear form of the heat flow equatlon. This equation tk analogous to the equation for the pressure in single-phase compressible flow The analysis reveals that the equation for the pressure in n~ultiphase flows can be obtained directly from the equation for the pressure in slngle-phase flows by simply replacing the single-phase compressibility by the total compressibility and the single-phase mobility by the sum of the mobilities of all the fluids present. Perrine's method of pressure buildup analysis for multiphme flows is based on these same substitutions. Thus, the development presented in this paper constitutes a theoretical justification for Perrine's method of pressure buildup analysis. One of the most difficult problems associated with the production of oil is the determination of the flow of the fluids within the reservoirs. One purpose of this paper is to present the development of simplified equations for the flow in gas drive reservoirs. The other purpose is to present some of the needed theoretical foundation for the methods currently being used in the analysis of multiphase pressure buildup data. The analysis presented in this paper makes use of a partial linearization of the equations governing multiphase flows. It is felt that the process of linearizing or partially linearizing the nonlinear equations of multiphase flows can lead to useful results since this process has proved extremely valuable in the study of other fields of science and engineering where nonlinear equations are encountered. Pressure buildup data of shut-in wells have been used for many years to estimate the static reservoir pressure. These data have also been used in studying the aver- age permeability both away from and close to the well-bores. The majority of the buildup data are obtained from wells producing two or three fluids. As yet the theoretical foundations for the methods of buildup analysis used in connection with rnultiphase buildup data have not been established. A number of methods of analysis of pressure buildup data have been developed." Most of these methods apply only to single-phase flows, and almost all of the theoretical analysis is confined to single-phase flows. Perrine' developed a method of analysis for rnultiphase flows from the method for single-phase flows of Miller, Dyes, and Hutchinson.' He did this by simply replacing the single-phase compressibility by the multiphase compressibility and replacing the single-phase mobility by the sum of the mobilities of the fluids in the rnultiphase flows. Obviously, such a procedure lacks a certain amount of theoretical justification, and such practices can lead to erroneous results. Fortunately, Perrine's results have a theoretical justification, and this justification is presented in the following discussion. It is assumed that gravitational effects and the cornpressibility of the rock can be neglected. Under these assumptions the equations governing the simultaneous flow of gas, oil and water in porous media areo* The preceding equations can be expanded into the following expressions.

By R. L. Whiting, H. J. Ramey

The basic study from which this paper was prepared was started as the result of the growing need throughout the world for increasing quantities of energy in all forms. Quite obviously, natural forms of energy that are readily available at low development cost are those in greatest demand. The underdeveloped countries — and particularly those having little or no petroleum resources — are the countries in which the most interest is being shown in the newer energy sources. One of the least expensive energy sources is natural geothermal steam.'-* Although this form of energy has been recognized for centuries, it has been only during the past 20 years that serious efforts have been made to harness it. Natural geothermal steam energy in Italy, New Zealand, Mexico, Japan, and California is now being produced through wells to drive turbines and generate electricity. Further, active exploration for natural geothermal steam is being conducted in Hawaii, Fiji, Taiwan, Chile, Russia, Greece, and Katanga. It is surprising to find that most geothermal steam exploration is in the "steam seep" stage. That is, surface studies are made and exploratory wells are drilled in the general area of steam seeps. However, in the larger geothermal steam areas, there have been efforts to apply the most modern geological and reservoir engineering principles in order to define the reservoir parameters, particularly those relating to estimates of reserves and future productivity. This paper is concerned with the development of appropriate equations and techniques to facilitate these estimates. Production of natural steam or hot water presents problems different from those experienced in the production of oil and gas. For example, steam or hot water systems may be essentially single-component systems, while hydrocarbon systems are most frequently multicomponent fluid systems. Heat effects are much larger for water than for hydrocarbon systems; and the natural steam production may or may not be isothermal, while production of petroleum reservoirs is considered normally to be isothermal. Petroleum reservoir engineering principles may be applied to natural steam or hot water reservoirs if the inherent differences in the systems are considered. The basic considerations involved in geothermal steam reservoir engineering are: thermodynamics, physical and thermal properties of water, materials and energy balances, fluid influx, and performance matching and predicting. Thermodynamics Fig. 1 is a pressure-temperature diagram for the liquid-vapor region for pure water, showing the critical point and five other points representing possible initial conditions for a geothermal steam reservoir. To generate a physical appreciation for the mechanics of natural steam production, it is helpful to consider

Jan 1, 1970

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970