Search Documents

By A. M. Khan, I. Fatt

A new type of steady-state thermal conductivity apparatus is presented and described. Results of measurements on Berea sandstone, Solehofen limestone, Bandera sandstone and Bakelite are presented. The measurements were made in a two-hour period at temperatures up to 80C and axial stress to 125 bars. Results show that the thermoelectric heat pump and tbermistor-controlled electric heater can be used in a steady-state apparatus for measuring thermal conductivity of mck. Steady-state is reached within two hours. With Urethane foam insulation, not more than 3 per cent of heat is lost by radial flaw. The apparatus can be used as a comparator or for absolute thermal conductivity determinations. INTRODUCTION Calculations of heat flow through the earth, both for scientific purposes and engineering design, require thermal-conductivity data on rock. However, the high variability of rock often makes handbook thermal-conductivity data of little value. For engineering purposes in particular, it is desirable to have measured data on samples of the rock mass under study. Laboratory methods for measuring thermal conductivity of rock use either time-dependent1-3 or steady - state temperatures .4,5 The time -dependent (transient) techniques give thermal diffusivity from which thermal conductivity can be calculated if the specific heat and density of the sample is known. Most transient techniques require lengthy and involved sample preparation. Also, the accuracy tends to be low. Steady-state techniques are more exact and are the most commonly used. The steady-state methods, however, also have their difficulties. A long time is often required to reach steady state, troubles have been encountered in establishing good thermal contact between the heating and cooling elements and the rock sample, and unaccounted-for heat losses are difficult to prevent. In this paper we present a description of a new type of steady-state thermal conductivity apparatus. With it, measurements have been made in two hours at temperatures up to 80C and axial stress to 125 bars. Sample preparation is simple and no thermal conductivity standard is needed.6 APPARATUS Our steady-state apparatus is an improvement . over the conduction pile used by Zierfuss and van der vliet4 and Somerton.5 A linear heat flux is established through a conduction pile consisting of the rock sample and a standard material of known thermal conductivity in series. The thermal conductivity of the sample is calculated from the ratio of the temperature drop across the sample

Jan 1, 1966

By A. E. Anderson, H. J. Hilll

Laboratory data for 155 field muds and 77 .shales have been used to develop correlations to estimate the net streaming potential component of the SP. Analysis of these data and comparison with field tests show that streaming potential corrections must be made to obtain an accurate estimate of the electrochemical potential for use in the determinntion of formation water resi-iiviy. The self-potential curve is related to formation water resistivity, and considerable effort within the oil industry has been directed toward clarifying this relationship. The Schlumbergers and Leonardon1,2 recognized that two types of potentials, electrochemical and electro-kinetic, resulted in the net potential variation in the borehole. Most of the later papers have dealt with the electrochemical component and its controlling parameters. The investigations of Wyllie,2" Salisch' and Moore5 provide the most extensive sources of data on streaming potentials across mud cakes. These authors observed streaming potentials of sufficient magnitude to be of critical importance in SP interpretation, if it is assumed that streaming potentials develop only opposite permeable beds. Streaming potentials across shales have been reported by Schenck and by Gondouin, et al.7 Since SP magnitude is the difference in potential opposite shale and sand beds, these authors suggested that the net streaming potential component may be negligible. The present investigation was initiated to evaluate empirically some of the parameters controlling the magnitude of the net borehole streaming potential and to determine whether useful relations between mud properties and net streaming potential could be established. Field muds were used in this investigation. These muds were characterized by the amount and type of chemical treatment utilized for five days before Sampling. Filter loss, mud resistivity, mud density and filtrate resistivity were determined for each mud with conventional methods and apparatus. Mud filtrate pH was estimated with Hydrion pH paper for the appropriate range. Streaming PotentIals Across MuD Cakes Commercial calomel-saturated KC1 electrodes are subject to potential changes attributed to contamination of the very small fiber-tip bridge by particles of clay from the drilling mud.3 To minimize this difficulty, a large diameter bridge (Fig. 1) was constructed with saturated KC1-agar gel and a filter disc of microporous porcelain. Experiments showed that streaming potentials across this bridge are less than 0.2 mv at pressure dif-

Born, March 28, 1853 Died July 23, 1916 On July 25, 1916, the Board of Directors adopted the following minute and directed that it be entered upon the records of the Board & published in the Bulletin of the Institute, and thereof be sent to the family of Mr. Kirchhoff: Charles Kirchhoff was one of the oldest, most widely known and most highly esteerned members of the American Institute of Mining Engineers. On his return in 1875 from a European course of study, to begin a highly successful professional career, he joined the Institute, then just commencing to acquire both national and international recognition. He became at once, and-for more than 40 years continued to be, its earnest , and active supporter, serving it as Manager, Councilor or Director for terms aggregating 9 years,-as Vice-President for 2 years, and twice as President for a year. During the critical period of the erection of the United Engineering Society's building, and the transfer of the Institute to permanent head-quarters therein, Mr. Kirchhoff was its efficient representative in the United Engineering Society, and on the Building Committee. In 1910, amid the congratulations of colleagues, employees, and innumerable friends, Mr. Kirchhoff retired from active business as Vice- President and General Manager of a great technical publishing company. Yet he sacrificed of his well-earned leisure yet one more year, to serve the Institute as President through 1911.

Jan 9, 1916

Jan 1, 1969

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

By E. J. Ripling, M. W. Toaz

Tensile tests were conducted on pure magnesium and on four commercial alloys over a variety of temperatures and strain rates. The high positive slope of the ductility vs testing temperature curves that is found over a short range of testing temperatures for these materaturerials was shown to result from the fact that the microstructure was not stable at these testing temperatures. Pure magnesium did not display a transition temperature, although its sub-room temperature ductility was low. This poor ductility resulted from a grain boundary weakness. The aluminum-bearing commercial magnesium alloys had a higher ductility than the pure magnesium. ness.Thealuminum-bearingcommercialThese alloys did exhibit a transition temperature and grain boundaries that were stronger than the grain material. Both pure magnesium and the commercial alloys began initiating cracks at strains equal to one half of their fracture strain. MAGNESIUM base alloys' enormous potential as structural materials lies in their high strength to weight ratios. At present, the high strength magnesium base alloys such as AZ-80, which will be discussed later, have about as high a tensile strength to specific gravity ratio as the highest strength wrought aluminum alloy, 75s-T6, and, of course, a much higher strength to weight ratio than structural steel. Further, since the magnesium base alloys develop their high strength to weight ratios with moderate strengths and extremely light weights, the ratio of load carrying capacity to the weight of many structures is relatively higher for the magnesium base alloys than their strength to weight ratios would indicate. This results because judicious increases in section size, such as increasing the depth of bending members or adding reinforcing ribs, can frequently increase the weight of a structure linearly while the strength and stiffness increase as the second or third power of the added weight. Strength and rigidity must be designed into lightweight structures by those techniques which establish multiaxial stress states. Consequently, strong, lightweight structures frequently also possess a greater propensity toward brittle failures. In order to insure against these catastrophic failures, an extra premium must be placed on ductility and toughness in light alloys. Current knowledge of the ductility and toughness behaviors—notched behaviors—of magnesium base alloys is limited.* Most of the magnesium base alloys * The state of present knowledge regarding the effect of testing temperature on mechanical properties has recently been summar-ized by Pritchard.' crystallize in the hexagonal system, so that they are expected to exhibit a tough and ductile to brittle transition temperature. Actually, a transition temperature behavior has been reported for commercially pure magnesium in the literature, but the characteristics of this transition temperatureh uggest that it may not be a transition in the usual meaning, but simply the onset of some thermal softening, a performance much like that which face-centered-cubic lead might exhibit slightly below room temperature. If transition temperatures do exist, such important questions as the manner in which they are affected by internal variables such as chemical composition and grain size, or the effect of external variables such as rate of straining or notch severity have never been investigated. The investigation reported upon herein was undertaken in order to define some of the mechanical behaviors of the magnesium base alloys better, to gain a better understanding of the influence of some of the above listed variables, and particularly to understand their effect on ductility and toughness better. Material and Procedure Material-—Pure magnesium and four commercial alloys were used in this investigation. The alloys were selected so that chemically they made up a series with approximately the same minor alloying

Jan 1, 1957

By Ogden Malin

IN conducting the manufacture of extruded brass rods and brass forgings it has been noticed that there is considerable variation in the physical properties, particularly the machinability of different rods and forgings, even though the composition, working temperatures and method of manufacture remain the same. It is the purpose of this paper to attempt to establish the cause and extent of this variation, and also to determine whether variations in composition within certain defined limits influence hardness and machinability to any marked degree. The writer believes that this variation probably is due to the variation in grain size of the metal, therefore the microstructure of several compositions under various conditions of working were examined in an effort to establish this as a fact.

Jan 1, 1932

By L. R. Raymond

With the trend toward deeper and consequently hotter holes, measurements of drilling mud properties at atmospheric temperatures are becoming increasingly inadequate.' Both the prediction and control of down-hole mud properties depend in part upon our knowledge of temperatures in the wellbore. Consequently, a better understanding of factors that affect temperatures during circulation and trips could improve drilling operations. It was toward the goal of obtaining this understanding that this study was directed. Other investigators2 have studied the response of fluid temperature in the borehole during trips. However, all of these studies have used an assumed formation temperature profile at the conclusion of circulation and have provided no means to calculate this profile directly. In this study, the formation profile during circulation was calculated. Other investiga-tors3-5 have also studied bottom-hole temperature during circulation. However, they have developed neither general techniques for calculating the entire temperature profile in the system nor generalized methods for predicting bottom-hole fluid temperatures during circulation. Such techniques and methods have been developed and are discussed briefly in the following pages. Temperature Behavior During Circulation Circulation of fluid during the drilling operation is represented schematically in Fig. 1. The process of circulation has three distinct phases: (1) fluid enters the drill pipe at the surface and passes down the drill pipe; (2) fluid exits the drill pipe through the bit and enters the annulus at the bottom; and (3) fluid passes up the annulus and exits the annulus at the surface. To simulate the thermal behavior of the fluid in the system, each of the phases of circulation must be described mathematically. In Phase 1, the fluid enters the drill pipe at a specified temperature, TDo. As the fluid passes down the pipe, its temperature is determined by the rate of heat convection down the drill pipe, the rate of heat exchange between the drill pipe and the annulus, and time. Phase 2 of the circulating process merely requires that the fluid temperature at the exit of the drill pipe be the same as the fluid temperature at the entrance of the annulus; i.e., TD(L,t) = TA(L,l). Thus in Phase 3, the fluid enters the annulus at TD(L, t). As the fluid flows up the annulus, its temperature is determined by the rate of heat convection up the annulus, the rate of heat exchange between the annulus and the drill pipe, the rate of heat exchange between the formation adjacent to the annulus and the fluid in the annulus, and time. These rates of heat exchange and the time dependency of mud temperature are described by well known heat-flow equation.6 Consequently, as shown in the Appendix, the temperature profiles in the drill pipe, annulus, and formation can be obtained by solving Eqs. 1 through

Jan 1, 1970

By J. Jones-Parra, R. S. Reytor

The porosities of fractured limestone reservoirs can be divided into two broad types in accordance with their effects on fluid distribution and fluid flow. In the coarse porosity, gravity segregation lakes place freely and the resistance to fluid flow is very small. In the fine porosity there is no segregation and a high resistance to flow, and it has relative permeability characteristics similar to tight sandstones. By analyzing the affect of the two porosities it is concluded that in some cases to recover the maximum amount of oil it is necessary to remove large quantities of gas from the reservoir by producing at high gas-oil ratios. In this manner the fine porosity is drained of its oil and the gas-oil contact drops slowly permitting higher production rates from the oil leg. Since this conclusion is contrary to widely accepted principles of conservation, a mathematical model (Fig. I) was constructed to duplicate the conditions desired. The behavior of the model indicates that under certain conditions it is possible to recover more oil by producing at high gas-oil ratios than by production at low gas-oil ratios, and that the rate of production is affected more by gas shut-offs than by the decrease in pressure (Figs. 2, 3, 4 and 5). INTRODUCTION It has been noted that certain limestone reservoirs behave in such a way as to indicate that there are two distinct types of porous spaces available for fluid storage and fluid flow. Regardless of the nature of these porous spaces the distinguishing characteristic is that in one type of porous space there is definite gravity segregation between the oil and the gas; while in the other, the gas evolved tends to remain distributed equally throughout the reservoir. For the sake of simplicity the porosity in which gravity segregation takes place at a fast enough rate to affect the behavior will be called "coarse po- rosity", while that in which this is not so will be termed "fine porosity". The behavior of a reservoir of this type would be different than that of a sandstone reservoir, especially if most of the oil is in the fine porosity. It has been observed that fluid flow to the wells takes place mostly through the coarse porosity. The large fractures which

By W. D. Johnston, T. P. Thayer

THE minerals that collectively are known as chromite form an isomorphous series of the general formula (Mg,Fe) 0. (Cr,Al,Fe) 203. So wide is the range in chemical composition in this group that chrome ore consisting entirely of the mineral chromite may range from less than 30 pct to more than 60 pct Cr203. The most comprehensive published data on the chemical com- position of chromite are by Stevens,62 who made complete analyses of 52 samples of chromite concentrate from many localities in North, Cen- [ ] tral, and South America. These samples were mechanically purified to remove associated silicates and might be regarded as the cleanest chromite that could be obtained by any mechanical milling process. Fig 1 shows the predominant zone of isomorphism, within the triangular prism of composition defined by the mineralogical end members of the series, wherein his analyzed samples fell.

Jan 1, 1949

By NONE N

Jan 1, 1963

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 A. M. Gaudin, J. H. Brown

In a survey of the effects of heat treatment on the comminution of various rocks, it was found that unusual size distributions often characterize the broken products, and that this is related to a great increase in intergranular fracture1. It is the purpose of the present paper to contribute to an understanding of this effect. In a thesis by one of use, critical consideration was given to many mechanisms which could conceivably be responsible for the increased intergran-ular fracture. It was shown that neither gross chemical changes nor intergranular corrosion are required. Among non-chemical mechanisms the most likely appeared to be (1) Uneven expansion due to thermal gradients, giving spalling (2) In a manner similar to that prevailing in metals, the appearance at high temperature of intergranular weakness. (3) Unequal thermal expansion of adjacent grains. Experiments were designed and conducted to permit a choice among the many alternative mechanisms and in particular among the three listed here. EXPERIMENTAL MATERIALS Tests were made on the materials listed in Table 1. Of each material, only lumps similar in composition and grain size were used. The three marbles listed in Table 1 were analyzed to assure that their chemical compositions did not differ greatly (Table 2). On the basis of the chemical analyses, of differential thermal analyses, and of X-ray powder diffraction patterns, it appeared that the marbles were very similar except for their grain sizes.

Jan 1, 1961

By Richard Miller

RECENT studies of creep have made it apparent that plastic deforma-tion may occur in metals under stresses less than the elastic limit as deter-mined from short-time tests. In summarizing conclusions drawn from creep tests by various investigators, Tapsell1 said: "It is evident that if a material has a true elastic range, the elastic limit must represent the maximum stress for no creep." The presence of a. definite elastic limit, represented in single crystals by a critical or threshold shearing stress, is claimed by Schmid2, Elam3 and others who have investigated the plastic deformation of single crystals. However, this "law of critical shearing stress" was deter-mined2 from short-time tests, and the measurement of the deformation was made with apparatus primarily designed to observe large amounts of plastic deformation, not to detect elastic changes. Experimental data recently obtained with highly sensitive instru-ments by A. W. Hanson4 on zinc and by Podaschowsky5 on rock salt indicate that the elastic range in single crystals is considerably narrower than formerly was supposed. Moreover, a theoretical possibility of permanent deformation, under the influence of even the smallest shear stress in single crystals above a certain temperature (Td), has been developed by G. I. Taylor6. This evidence seemed to warrant an examination of the elastic limit in large single crystals by means of creep tests. The first metal chosen for examination was zinc.

Jan 1, 1936

By M. Cross

INTRODUCTION The results of the Japanese dissective investigations have [l-4] transformed our understanding of how the blast furnace operates. These studies revealed an internal structure of the furnace burden with five main zones, as illustrated in Figure 1. In the granular zone the charged material retains its structure whereas in the cohesive zone the iron ore particles are at advanced stages of reduction and form a semi-molten mass with a relatively high flow resistance. In the active coke zone molten iron and slag flow through a matrix of loosely packed coke. The raceway is a relatively small void space generated by displacement of coke due to the hot blast momentum, whilst the hearth is packed with coke, slag and molten iron. [ ] It is now well appreciated that furnace control is largely achieved and maintained by using the burden charge distribution to determine the gas glow distribution within the furnace and, hence, the shape and extent of the cohesive zone. In fact, Wakayama et a1 [5] of Nippon Steel Corporation state that optimum furnace operation is synonomous with the formation of a desirable cohesive zone shape, which in turn is almost solely determined by the burden charge distribution and its mechanical properties [6]. Much of the development work on operating furnaces since the dissective investigations has been oriented towards the identification of external furnace measures which reliably predict the shape of the cohesive zone. Currently the most popular measures are top gas temperature, velocity and CO/CO, ratio profiles [5]. More recently pressure distributions as evaluated from a series of wall measurements down the height of the furnace are also being utilised [7] to identify the position and shape of the cohesive zone. A great deal of emphasis has been placed upon achieving a consistent quality in the burden charge materials, establishing a suitable set of operating conditions and only making changes when forced to, ie., the operation for some reason moves outside a desired operating range or the production requirements change. Operational control schemes based upon this philosophy (eg., Kawasaki's GO-STOP system [8]) have been remarkably successful at achieving both high production rates and relatively low fuel costs. Mathematical models have long been used as tools to assist in the analysis of the blast furnace process. Although, the models have been useful in post-operative analysis and for identifying trends, etc., they have not proved reliable for detailed design criteria or as the basis for effective control schemes. Such mathematical models have moved from simple statistical correlation, through overall and stagewise heat-and-mass balance calcu-

Jan 1, 1984

By C. S. Samis, D. P. Seraphim

In the presence of oxygen, galena is oxidized in an aqueous medium containing ammonium acetate in accordance with the following reaction: PbS + 1/2 0, + 2 NH~Ac -» PbAc, + So + 2 NH: + H2O. This reaction was studied in an autoclave under oxygen pressure by polarographic measurement of the concentration of lead ion in solution. Oxidation rate is parabolic. The reaction products are lead acetote and elemental sulfur, the latter forming a film on the galena crystal. Reaction kinetics appear to be determined by the diffusion of lead from the galena phase to the solution interface through the film of sulfur. RECENT applications of pressure leaching on a commercial scale have emphasized the metallurgical importance of the oxidation of sulfide minerals in aqueous media by oxygen under pressure. Research experience has shown that sulfates, thio-nates, and metal oxides are the products of humid oxidation of sulfide minerals in basic solutions. The metal oxides may be solubilized in the presence of a complexing ion.' Certain oxide minerals are also amenable to treatment by this type of oxidation.' Several kinetic studies have been made of reactions of this type. Stenhouse' found that pyrite (FeS?) oxidized in a sodium hydroxide solution to produce a soluble sulfate and an insoluble layer of iron oxide and that the kinetics of the reaction were controlled by oxygen diffusion through the oxide layer. Andersen' studied the oxidation of galena (PbS) in a sodium hydroxide solution under oxygen pressure. He found the lead and sulfur to be simultaneously oxidized to produce plum bite ions and sulfate ions, and that the reaction rate was proportional to the surface area of the galena. The kinetics appeared to be determined by the reaction of oxygen and water with the galena surface. In each of the above reactions the metal and sulfur components of the mineral were oxidized simultaneously and at equal rates. In the present study the alkaline medium has been replaced by an approximately neutral solution of ammonium acetate in which lead sulfate is soluble. Under the conditions no sulfate was produced and no sulfur could be detected in the solution. The metal component of the mineral was solubilized and the liberated sulfur remained as a film of elemental sulfur on the surface of the galena crystal. A kinetic study of this reaction was undertaken to establish the kinetic mechanism and to determine the rate controlling step. Specific reaction rates were determined by oxidizing crystals of galena of measured surface areas in an autoclave maintained at the desired temperature and oxygen pressure. The rate of reaction was followed by measuring the concentration of lead in the ammonium acetate solution with a cathode ray polarograph, employing stationary platinum electrodes incorporated in the autoclave. Experimental Materials and Equipment Crystals of galena obtained from Violamac Mines, Retallack, B. C., were selected for oxidation because

Jan 1, 1957

By D. F. Hewett

CONTENTS PAGE Introduction 1 Geologic features of eastern Oregon 2 Geology of the Sumpter quadrangle 2 Rocks older than the ore deposits 3 Rocks younger than the ore deposits 8 Lodes of the Sumpter quadrangle 8 Areal distribution 8 Geologic distribution 9 Major features of the veins 9 Mineralogy and texture of the veins 13 Ore deposits zonally disposed around the Bald Mountain batholith . . . . 13 Summary of mineralogy and texture 32 Zonal relations 33 Ore deposits zonally disposed around Greenhorn Mountain batholith . . . 36 Zonal relations 37 Relation of textural features to zonal arrangement of veins 37 Summary of zonal relations 38 Comparison with similar districts 39 Outlook for the district 40 INTRODUCTION IN an earlier report on the Sumpter quadrangle1 it was briefly stated that the quartz lodes of that region might be separated into several groups or belts on the basis of the kind and abundance of the minerals of which they are composed, and that these belts were disposed as areal zones around the principal intrusive mass of the region. When the opportunity to revisit the region was offered in 1929, as a result of a cooperative agreement between the State Mining Board of Oregon and the United States Geological Survey, it seemed advisable to pay particular attention to this feature of the lodes, in the hope that if the idea of zonal relations was confirmed it might serve as a guide in exploitation.

Jan 1, 1931

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

By R. G. Turner, M. G. Hubbard, A. E. Dukler

Gas phase hydrocarbons produced from underground reservoirs will, in many instances, have liquid phase material associated with them, the presence of which can affect the flowing characteristics of the well. Liquids can come from condensation of hydrocarbon gas (condensate) or from interstitial water in the reservoir matrix. In either case, the higher density liquid phase, being essentially discontinuous, must be transported to the surface by the gas. In the event the gas phase does not provide sufficient transport energy to lift the liquids out of the well, the liquids will accumulate in the wellbore. The accumulation of the liquid will impose an additional back pressure on the formation that can significantly affect the production capacity of the well. In low pressure wells the liquid may completely kill the well; and in the higher pressure wells there can occur a variable degree of slugging or churning of the liquids, which can affect calculations used in routine well tests. Specifically, the calculated bottom-hole pressures used in multirate backpressure tests will be erroneous if the well is not removing liquids on a continuous basis, and gas:liquid ratios observed during such a test may not be correct. Several authors1,3,8,14 have suggested methods to determine if the flow rate of a well is sufficient to remove liquid phase material. Vitter14 and Duggan' proposed that wellhead velocities observed in the field would be adequate for keeping wells unloaded. Jones8 and Dukler3 presented analytical treatments resulting in equations for calculating, from physical properties, the minimum necessary flow rate. An analysis of these studies indicates the existence of two proposed physical models for the removal of gas well liquids: (1) liquid film movement along the walls of the pipe and (2) liquid droplets entrained in the high velocity gas core. Although there probably is a continuous exchange of liquid between the gas core and the film, they will be treated separately for the purposes of this study. The development and comparison of these separate models with experimental data will permit the determination of which, if either, is the controlling mechanism for the removal of liquids from gas wells. The Continuous Film Model Liquid phase accumulation on the walls of a conduit during two-phase gas/liquid flow is inevitable due to the impingement of entrained liquid drops and the condensation of vapors. The movement of the liquid on the wall is therefore of interest in the analysis of liquid removal from gas wells. If the annular liquid film must be moved upward along the walls in order to keep a gas well from loading up, then the minimum gas flow rate necessary to accomplish this is of primary interest. The analysis technique used follows Dukler2 and Hewitt5 and involves describing the profile of the velocity of a liquid film moving upward on the inside of a tube. The minimum rate of gas flow required to move the film upward is then calculated.

Jan 1, 1970

By S. G. Gibbs, A. B. Neely

Jan 1, 1967