Search Documents

By Lyon F. Terry

IN contemplating the prospects of natural gas being transported from the fields where it is produced to such distant points as Philadelphia, New Jersey, New York City, and New England, let us review the advance that natural gas has made and the projects currently scheduled for construction. This may give us a perspective of the point in the long-range development of the natural gas business in which we now find ourselves, and some basis for judging the further expansion of the industry. Authentic statistics of the total marketed production of natural gas in the United States have been compiled by the Bureau of Mines since the year 1906. The growth of the business is indicated by the following data on marketed production for each ten-year interval since 1906:

Jan 1, 1947

By R. T. Elworthy

Canada ranks second in the list of world consumers of natural gas. In 1923 our production was nearly fifteen million thousand cubic feet valued at $5,800,000. Our output, however, appears very small compared with that of the United States, which in 1923 totalled over one thousand and eight million cubic feet, valued at $240,000,000. Natural gas is found in commercial quantities in three? provinces in Canada. New Brunswick: Natural gas in New Brunswick has been utilized commercially only since 1912,although the investigation of oil and gas possibilities of the province commenced as early as 1859. The annual production is a little more than one-twentieth of the total Canadian production and remains fairly constant from year to year. Ontario: In Ontario the earliest gas wells were drilled in Essex and Welland counties between the years 1885 and 1890. On account of the small demand in Canada the cities of Detroit and Buffalo were supplied with large quantities of gas from these wells, but this export was curtailed in 1898 and finally ceased in 1900. The production of natural gas in Ontario, mainly from fields in Lambton, Kent, Essex, Norfolk, Haldimand and Welland counties, increased rapidly from 1900 until the peak was reached in 1917, when over nineteen million thousand cubic feet valued at $3,641,587 was consumed.

Jan 1, 1925

By J. H. Henderson, H. J. Ledbetter, N. B. Gove, J. D. Griffith

INTRODUCTION The results of a series of laboratory flood tests using liquid Iso-butane to displace refined oils from test cores are pre- ented and interpreted on an empirical bask. The study Revealed the similarity of the miscible liquid displacement to that of the immiscible-fluid displacement mechanism. The efficiency of the iso-butane flood decreased markedly as the oil Viscosity increased. but the effect of injection rate on the effectiveness of the primary production stage was negligible over the range investigated. The presence of free-gas saturation prior. to iso-butane breakthrough increased the volume of iso- Butane required to recover a given percentage of the oil present. Injection of liquid i.e.-butane prior to water follicle sculpted in a marked improvement in oil recovery by the water flood. INTRODUCTION The production of petroleum is often enhanced by the injec-1i1)n of extraneous fluids into the reservoir. It has been suggeited that the injection or cycling of a wet gas or a liquefied petroleum gas would he instrumental in increasing the effic.ienr!; of oil recover! to a greater degree than the commonly injected fluids. dry pas and water. Although instances of injection of light liquid hydrocarbons are not unknown, a systrmatic recording of the result is unavailable and superfi- cially, the practice gives the impression of being economicaIIy unsound. For specific types of displacement, analytic: solutions have been advanced which aid in understanding the physical mechanisms involved. However, there is no generalized theoretical concept which permits adequate prediction of displacement efficiency from a knowledge of readily measured reservoir rock and fluid properties. The Buckley-Leverett approach to oil displacement by an immiscible phase is an admirable analysiy of that type of mechanism. Dykstra and Parsons have introduced a useful concept in their evaluation of the role of the mobility ratio and its merit in rationalizing displacement phenomena. Everett, Gooch and Calhoun3 in their report on the effect of viscosity ratio and miscibility have pointed out the lack of a comprehensive theoretical approach and the desirability of investigating the mechanism of miscible-miscible displacement. In lieu of a well-developed theory concerning miscible-111isci1)le displacement and in view of the scarcity of engineering information on actual LPG injection trials, the limited program of laboratory experimentation described below was undertaken. OBJECTIVES The original broad objectives of the program, to study the oil recovery process using LPG as the displacing phase, defined the direction and scope of the laboratory tests. The specific objective. of the experimentation were to evaluate

Jan 1, 1953

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 P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By A. J. Garnier, N. H. van Lingen

Rock downhole is known to be lesc. drillable than when brought to the surface. This must be ascribed mainly to the presence under downhole conditions of a pressure differential across already made chips, which hinders their being lifted. The pressure differential has partly a static and partly a dynamic origin. Balling-up of bits is another consequence of this pressure differerztial. Reduction in penetration rate owing to an increase in the strength of the rock is governed by the difference between the mud pressure and the pressure of the for/nation pore liquid. Rotationally symmetric geostatic stresses have 170 effect on drillability. Fracture of rock when drilling will be brille in most cases. The above is supported by laboratory drilling experiments with drag bill and roller bits an elevated mud, pore, and confining pressures On rocks differing in strength and permeability. INTRODUCTION In oil well drilling drillability of rock is found to decrease with increasing depth of the hole. Naturally deep rock will be more compacted and, therefore, harder to drill than shallow rock of the same type. However, apart from this the drillability of a sample of deep and compacted rock brought to the surface is generally much higher than in its original location downhole. In view of the economic implications of this reduction in drillability, it seems worthwhile to analyze its causes. The origin clearly has to be sought in the difference of environment. The only conceivable factors would seem to be the presence of mud under pressure, the pressure of the formation pore liquid, and the overburden of the rock. Down the borehole the rock is compressed triaxially by mud pressure and overburden. It is well known that the strength of rock is increased when confined by external pressure1. Various authors have, therefore, ascribed the difference in drillability mainly to the strengthening of rock by triaxial compression. Another factor, mentioned by Bobo and Hoch5, is that forces. including "pressure differential forces", tend to hold a dislodged particle in place. However, the conditions determining their magnitude are not clarified nor is their effect on drilling rate assessed. In the laboratory, drilling experiments on pressurized samples of rock yielded evidence that pressure differential forces holding the chips down are the major factor in reducing rate at depth. This paper describes the experiments and shows how the results enable both a qualitative and a quantitative interpretation of the factors determining the effect of chip hold-down on drilling rate. The implication with respect to balling-up and jet action is also discussed. The greatcr part of the experiments were performed in the pressure vessel shown in Fig. 1. The spacc between the sample of rock and the vesscl is divided by "0" rings into three separate chambers. The rock sample is confined laterally by pressurized oil in thc middle chamber. Penetration of oil into the sample is prevented by partly jacketing it with brass foil. Thc pressure of the drilling fluid in the hole can be adjusted via the upper chamber. Permeable rock specimens arc water saturated before being drilled. With a properly plastering mud as drilling fluid the pressure of the pore water can be adjusted independently via the lower chamber.

By J. G. Richardson, K. H. Coats

During the initial growth of a gas bubble in an aquifer storage reservoir the injected gas tends to override the water. The resulting low displacement efficiency and high rate of gas travel down-structure make it difficult to maintain water-free production. This paper illustrates the use of two-dimensional, two-phase calculations to simulate this gas-water displacement. Calculations were performed for a stratified aquifer and for several homogeneous sands differing in permeability, porosity and thickness. Results are discussed and compared with performance predictions obtained from the simpler Buckley-Leuerett and Dietz formulas. This comparison indicates the necessity for the two-dimensional calculations in realistically simulating the gravity override of the water by injected gas. INTRODUCTION In recent years natural gas has been stored near markets in aquifers where insufficient storage capacity is available in depleted fields. Operators of these aquifer storage reservoirs have encountered technical problems relating to the gas-water displacement accompanying initial growth of the gas bubble. Injected gas tends to override the water, with a resultant low displacement efficiency and high rate of gas travel down-structure toward spill points. Low displacement efficiency makes if difficult to sustain water-free gas production. Sometimes the fingering of gas down-structure may be so pronounced that the injection rate must be severely curtailed, which excessively lengthens the time required for bubble growth. This paper is concerned with estimating — for given aquifer characteristics and fluid properties — the displacement efficiency, rate of gas movement down-structure and rate of gravity drainage of water behind the gas front. Several published articles relate to simulative capability necessary to handle this problem. The Dietz formula1 and Buckley-Leverett2 method have some applicability to the problem; Woods and comer3 reported one-dimensional, radial calculations which accounted for the two-phase flow of gas and water. Most recent research directed toward understanding aquifer behavior in relation to gas storage has concentrated on single-phase flow in aquifer.4 Douglas, Peaceman and Rachford5 presented a method for calculating multidimensional two-phase flow in reservoirs; their method has been applied in work reported by Nielsen and Tek,6 Blair and peaceman7 and Goddin.8 Refs. 5 and 7 show comparisons between experimental data and calculations similar to those employed here. One purpose of this paper is to illustrate the use of two-dimensional, two-phase flow calculations in simulating the difficult problem of water displacement by gas in a vertical cross-section. The pronounced fluid density and viscosity differences combine with the disparity in horizontal and vertical dimensions of the cross-section to pose one of the more difficult problems in reservoir simulation. The calculations described account for capillary and gravity forces, relative permeability and reservoir heterogeneity. An example reservoir is described and calculated results are presented for a variety of injection rates and values of permeability, reservoir thickness and dip angle. A second purpose is to compare the two-dimensional calculations with results obtained from the Buckley-Leverett and Dietz formulas. For this reason most calculations were performed for a simplified example reservoir to which the latter techniques might reasonably be applied. However, the computer program employed applies equally well to cases involving arbitrary spatial variations of permeability, porosity, reservoir thickness and dip angle.

By P. Pollard

A method has been developed for evaluating acid treatments in fractured limestone fields by breaking down pressure drawdown into three component parts: (I) pressure differcntial across "skin" near the borehole face, (2) pressure differential due to flow resistance in the coarse communicating fissures and (3) pressure differential between the fine voids and the coarse fissures. It is apparent that in most successful acid treatments the first term, skin resistance, has been reduced or eliminated. Further, it is often possible to estimate the volume of coarse fissures associated with the second term, coarse fissure flow resistance. In cases where this volume is comparable with practical acid volumes it seems likely that this resistance also may be attacked with a suitably retarded acid. INTROD UCTION Acid treatment has been successfully applied as a general practice in the limestone wells in the Mara/ Maracaibo districts of Venezuela for some 10 years. At the same time the need has been felt for a more precise method of: (1) evaluating the effect of an acid job and (2) selecting wells which should benefit from an acid treatment. The following method, based on an analysis of build-up curves for wells producing from fractured limestone, provides a means of forecasting the effect acid will have on a well and in many cases the probable order of production rate increase which may be expected. In certain cases it is also possible to calculate the fissure volume through which the acid has to be displaced for best results. Various approximations are introduced in the method, but within the limitations of these the analysis gives an indication of how a well, whether previously acidized or not, may be expected to respond to acid treatment. Suitable build-up curves cannot be obtained on "tight" wells which do not give sustained production and the method is not applicable to such wells. Figs. 1 through 7 illustrate various types of build-up curves obtained from limestone wells in Western Venezuela. Curves are included showing (1) both skin and coarse fissure resistance, (2) skin resistance only, (3 and 4) reduction of skin resistance by acidization, (5) fissure resistance only, (6) reduction of coarse fissure resistance by acidization, and (7) acidization failure. This theory has been successfully applied for over three years during which 22 acid treatments have been carried out specifically to reduce skin resistance and three treatments to reduce coarse fissure resistance. The method has been found to be a useful aid in the selection of wells for acidization.

By K. L. Young

The general energy equation, including change in kinetic energy, was solved by numerical integration and used to evaluate simplifying assumptions and application practices over a wide range of conditions. When extreme conditions were encountered, sizable errors were caused by large integration intervals, application of Simpson's rule and neglecting change in kinetic energy. A maximum error of only 1.31 percent was caused by assuming temperature and compressibility constants at their average value. It was discovered that a discontinuity can develop in the integral for the injection case. This discontinuity indicates a point of zero pressure change and is an inflection point in the pressure traverse. INTRODUCTION When a pressure in a gas well is to be calculated, one of the first decisions is to select a method of calculation. In many instances, this selection becomes a problem because the literature, at best, provides an evaluation of any method for only a limited range of conditions. Once a method has been selected, a question often arises as to the size of the calculation interval which should be used. The question regarding calculation interval arises because an analytic solution is not obtainable and approximate solutions must be used. This paper presents an evaluation of major assumptions and application practices of probably the two most widely used methods for calculating steady-state single-phase gas well pressures. The two methods are Cullender and Smith1-" (numerical integration), and average temperature and com-pressibility. The Cullender-Smith method assumes that change in kinetic energy is negligible and is normally applied in two steps with a Simpson's rule correcton. The average temperature and compressibility method, in addition to neglecting kinetic energy change, asumes that temperature and compressibility are constant at their average values. This method is normally applied for wellhead shut-in pressures of less than 2,000 psi, and in one step. Computer programs were written to compute bottom-hole pressure with and without the assumptions, using various approaches. Values of input parameters investigated are shown in Table 1. Flow rate was limited to a maximum of 5,000 and 10,000 Mcf/D for tubing sizes of 1.610 and 1.995 in. ID, respectively. Flow rate was also limited to 10,000 Mcf/D for a tubing size of 2.441 in. ID when wellhead flowing pressure was 100 psia. These limitations were imposed on flow rate so as not to exceed sonic velocity. The z factor routine available necessitated limiting bottom-hole temperature to 240F and wellhead pressure to 3,000 psia. Pressures were compared on the basis of percent deviation from the trapezoidal integration of Eq. l or 2 at 100-ft intervals. A preliminary investigation indicated that a 1,000-ft interval solution would differ from a 50-ft interval solution by less than 0.25 percent; therefore, the 100-ft interval was chosen for a base. For the purpose of comparison, deviations less than 1 percent were considered insignificant. where 111.1 qeC2/d p = kinetic energy term. Eqs. 1 and 2 can be evaluated numerically at specific depths using the trapezoidal rule as shown by Cullender and Smith. If change in kinetic energy is neglected and temperature and compressibility are assumed constant at their average values, Eq. 1 can be integrated to give the average temperature and compressibility equation,

By John Riegle

To date. utility company underground storage of gas has generally been restricted to depleted dry gas fields. The Playa del Rey project is probably the first to successfully store gas in a partially depleted oil reservoir with the recovery of large volumes of pas at high rates as the main objective. Oil recovery has been a secondary consideration. Problems encountered unusual to those of storage in a dry gas zone were: (1) the removal and retardation of formation of emulsion. (2) the upstructure movement of fluid during withdrawal periods which formed fluid blocks, and (3) reservoir shrinkage resulting from encroachment of edgewater. The solution of these problems as outlined in the paper has resulted in increasing the withdrawal rate from the original design capacity of 4,000 Mcf per hour to the current 10.000 Mcf per hour. Additional increase is anticipated as the operations continue. Migration. reservoir performance. operational procedure and a historical record are included to make this a resume of the project from its war-time inception in 1942 to the first of 1952. INTRODUCTION The Playa del Rey underground gas storage project is probably the first project in which gas has been stored in a partially depleted oil reservoir, with the recovery of large volumes of gas at high rates as the main objective and the recovery of oil a secondary consideration. Problems unusual to those of storage in dry gas zones have made this a pioneer endeavor. The total storage capacity of approximately 1,500,000 Mcf of this reservoir is small in comparison with other underground gas storage projects. However. by study and experiment the deliverability into transmission lines. operating at pressures in excess of 150 lbs., has been nearly tripled in the past three years. to the current 10.000 Mcf per hour late, without the use of compression. It is anticipated this rate will gradually increase as additional fluid is removed from the formation. LOCATION AND HISTORY Playa del Key is located on Santa Monica Bay about two miles south of the beach commuity of Venice and 15 miles southwest of the renter of the City city Los Angeles. Fig. 1. The discovery well was completed at 6,194 ft in 1929. An upper zone at 3,905 ft was found on June 18. 1930. Townlot development was very rapid and by the end of 1930 there were 141 producing wells in the field. The Union Oil Co. completed a lower zone well south of the Townlot Field in May, 1931, which was believed to be an extension, until additional drilling of about 50 wells during 1934 and 1935 indicated this area to be a separate accumulation as shown on Fig. 2. The lower oil horizon, which is now the storage zone, had no original gas cap and is believed to be a continental detrital deposit of conglomerate that accumulated in valleys and depressions of a weathered schist high. which later subsided and was overlain by an impervious nodular shale. This type of formation has rapid variation of porosity. permeability and thickness due to lack of sorting action common to marine deposition. and the uneven surface of the original schist high. Production was limited on the upstructure sides by the pinching out of the conglomerate, and downstructure by edgewater. Unrestricted production, together with high permeabilities, led to large wastage of gas from this new area as facilities were not available for transport of gas to market. The decline was very rapid and by 1942 the portion of the new area south of Ballona Creek was nearly depleted. ACQUISITION OF PROPERTY The rapid increase in industrialization of Southern California during World war 11 forced upon the local gas com-

Jan 1, 1953

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By L. H. Wolfe, J. O. Thieme

The tensile and shear strengths of frozen soil and the compressive strengths of ice and frozen soil were measured These tests showed that the strength of ice and of frozen soil increased as the temperature was decreased. A method is presented for the measurement of thermal conductivity and specific heat of earth and ice. Using the parameters thermal conductivity, specific heat and density, the thermal diffusivity was calculated. Results of thermal property tests establish water as the major variable contributing to thermal property values. INTRODUCTION During the last few years, a considerable interest has developed in underground storage of cryogenic liquids such as liquid natural gas. Because of its economy and safety, underground storage is being used. Before underground caverns could be seriously considered, however, the strength and thermal characteristics of the soil had to be obtained. Methods are presented herein for measuring tensile, compressive and shear strengths, specific heats and thermal conductivities. The tests were performed on two soils at temperatures from ambient to -195C. Because it is difficult to find data in the literature that pertain to a particular soil, ice was used to evaluate some of the experimental procedures. Some of our ice data are compared to the published values. PHYSICAI, PROPERTIES SOIL SELECTION AND SAMPLE PREPARATIONS The physical and thermal tests were performed on recompacted samples of two typical soils that might be encountered in underground caverns: a gray, fire clay that had been pulverized and mixed with water into a pliable mud; and a brown, sandy silt that was dug from a dry river bank and subsequent1y mixed with tap water. The moisture content of the test specimens was 17 to 22 per cent (based on the wet weight of the soil). This amount of water nearly saturated the silt, but the clay was well below the saturation point. The soils are partially described by the sieve analyses (Fig. 1) which show the particle size distribution of the coarser than 44 micron portion of the soils. The curves also show that approximately 75 per cent of the clay and 20 per cent of the silt are finer than 44 microns. The particle size seems to affect some of the physical properties, and the finer than 44 micron portion is important in the thermal conductivity test. Physical specimens were cut with a bandsaw from frozen blocks of the soil. The specimens were then sanded and kept frozen until after they were tested. They were stored at the approximate test temperature for several hours. Then, to equilibrate the samples to the test temperature, they were stored at the temperature at least an hour before they were tested. TEST EQUIPMENT An Instron, which is a constant speed gear-driven instrument, applied force to the samples and drew stress-strain plots of the test. All of the samples were strained at 0.02 in./min. All of the physical tests were conducted in a temperature cabinet that controlled the temperature to 1C. TENSILE STRENGTH Briquette or "dog-bone" shaped tensile specimens were pulled with suitably shaped jaws. This shape of specimen, which is usually used for tensile tests of concrete, is satisfactory for testing frozen soils except that the samples tend to

Jan 1, 1965

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

By G. Thodos, R. B. Grieves

A method has been developed for the accurate calculation of the cricondentberm and cricondenbar temperatures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain any number of components including paraffins and isoparaffins, olefins, acetylenes, naphthenes and aromatics. This approach is based upon the mole fraction of the low-boiling component in the mixture and graphically presents the ratios of the cricondentherm and cricondenbar temperatures to the pseudocritical temperature as {unctions of a boiling-point parameter. The only requirements are the critical temperature, normal boiling point and the approximate vapor pressure behavior of each component. For mixtures of more than two constituents a stepwise calculation procedure is necessary where the pseudocritical temperature is based upon the critical temperature of the pure low-boiling component and upon the actual cricondentherm or cricondenbar temperature of the mixture of all the remaining higher-boiling components. From an analysis of 22 binary systems (123 compositionsl, the average deviation of calculated cricondentherm temperatures from reported values is 0.87 per cent, based on degrees Rankine; and for 10 multicomponent mixtures containing from three to six components, the average deviation is 0.98 per cent. From an analysis of 18 binary systems (104 compositions) the average deviation of calculated cricondenbar temperatures from reported values is 0. 79 per cent; and for 15 multicomponent mixtures, 1.47 per cent. Equations, derived from the graphical relationships, are presented which enable rapid calculations for both binary and multicomponent systems. INTRODUCTION With the increasing tendency towards the use of high pressure processes, the need for an accurate method of calculation of the properties of a multicomponent hydrocarbon mixture in the critical region is becoming more essential. Critical region phase relations are also significant in gas condensate reservoirs, for which a knowledge of the fluid phase boundaries and regions of retrograde condensation make it possible to evaluate optimum operating conditions. The possibility of liquefaction in the underground structure may be established, and information on the feasibility of the use of repressurizing may be determined. A typical phase diagram for a multicomponent hydrocarbon mixture is presented as Fig. 1. This figure also shows the regions of retrograde behavior encountered by following lines AB and AC. If accurate methods of calculating the temperatures and pressures at the critical, cricondentherm (maximum temperature) and cricondenbar (maximum pressure) points of multicomponent hydrocarbon mixtures of known composition were available, it would be possible to estimate closely the entire phase diagram of any mixture. Several methods have been presented in the literature for the estimation of the critical temperatures and pressures of multicomponent mixtures containing all types of hydrocarbons and including the fixed gases.4-8,21 Considering the cricondentherm and cricondenbar points, the work is more limited. Etter and Kay5 have developed rather complex equations for

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 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 A. L. Lee, M. H. Gonzalez, R. F. Bukacek

Experimental viscosity data for methane are presented for temperatures from 100 to 340F and pressures from 200 to 8,000 psia. A summary is given of the available data for methane, and a comparison is presented for that data common to the experimental range reported in this paper. Correlation of the data is presented, and predicted values are given for temperatures up to 460F and pressures up to 10,000 psia. INTRODUCTION The increasing ranges of temperature and pressure at which fluids are produced, transferred and processed in the petroleum and chemical industries stress the need for accurate information on physical properties, both for engineering calculation and for improving the methods used to estimate physical properties. A survey of the literature reveals that disagreements between published data on the viscosity of methane are common and that most investigations have been conducted over restricted temperature and pressure ranges. This situation made it desirable to undertake an investigation with a scope common to most of the literature data available. APPARATUS AND MATERIALS The apparatus used in this investigation was the capillary tube viscometer'described by Dolan,10 with modifications introduced in the general design and operation of the instrument. It has an effective pressure range from 14.7 to 10,000 psia and a temperature range from room temperature to 400F. The design of the viscometer is based on the establishment of a manometric head between two vessels containing the test fluid and a volume of mercury. The reservoirs are connected by a capillary tube through which the test fluid flows and a tube through which mercury flows. A pressure gradient is established by elevating one of the vessels above the other, the resulting flow of mercury displacing the test fluid through the capillary. The schematic diagram of the system (Fig. 1) shows the arrangement of the equipment auxiliary to the viscometer. The density cells assembly (E, Fig. 1) was not used in this investigation since reliable data are available on methane. 23 The assembly consists of a bank of eight 316 stainless steel pycnometers. Experimental values of the density of a natural gas19 and i-butane14 samples have been obtained for pressures up to 8,000 psia and temperatures up to 340F. The methane used was obtained from the Southern California Gas Co. Mass spectrometric analysis showed a composition of 99.8 percent CH4, 0.1 percent C2H6 and 0.1 percent C3H8. EXPERIMENTAL VALUES AND ACCURACY OF DATA Determinations were carried out at several temperatures. With the temperature in the cabinet controlled at a given level, a series of runs were made at various pressures. To establish pressures lower than that obtained from the methane cylinder, it was necessary to vent the gas slowly to the atmosphere until the desired pressure was reached. The consistency of the test run time for a particular condition can be used as a criterion for the accuracy of the data. A further check is obtained by comparison of the viscosity values obtained using different driving heads. To establish the instrument's reproducibility, a large number of runs were made at three completely different conditions of temperature and pressure. Data obtained were analyzed statistically to determine the most reliable value for the variance of the experimental data. It was found that the variances of the data taken at different experimental conditions were not statistically equivalent, but they were small. The largest variance was found at 8,000 psia and 100F, where the standard deviation was 0.385 micropoise. Therefore, 99.9 percent of the data can be expected to fall within 1.15 micropoises of the population mean, except for the intervention of systematic errors in operation.

By M. R. Tek, L. K. Thomas, D. L. Katz

Threshold displacement pressures are needed to determine how much overpressure can be used in storing natural gas. An experimental technique for determining threshold pressures by displacing water with gas from samples saturated with water is presented. Threshold pressures for eight low permeability samples were measured. Threshold pressure data obtained in this work, plus data on higher permeability samples reported in the literature, are correlated with porosity, permeability, surface tension and formation resistivity factor. Mercury injection pressures also were measured and correlated with air-water threshold pressures. Threshold pressures were found to be independent of time. An aerial photograph of gas emerging from the top of a core shortly alter its threshold pressure has been exceeded shows that the gas bubbles are uniformly distributed across the face of the core. Channeling will occur, however, when an increased gas phase permeability is reached, A porous medium can be re resealed alter its threshold pressure has been reached, provided it has not been desaturated below a fixed saturation. Its new threshold pressure will be lower than when the sample is 100 percent saturated with water. INTRODUCTION The degree of overpressure in excess of the discovery pressure that a gas storage reservoir can withstand is determined by the ability of its caprock to contain gas, providing the pressure does not exceed the structural limit of the reservoir. The retention of gas in a reservoir by a caprock saturated with water is a result of the capillary forces acting at the gas water interface. Without the presence of water in a caprock, gas would leak out of the reservoir at a rate determined by the permeability of the caprock to gas. The ability of a caprock to contain gas is expressed in terms of its threshold displacement pressure, which is defined as the minimum pressure needed to initiate the displacement of a wetting phase by a nonwetting phase from a porous medium 100 percent saturated with the wetting phase. The gas industry must be able to predict how much "overpressure" a gas reservoir can withstand before leakage occurs by gas displacing water from caprock overlying the reservoir. This is of economic importance since the storage capacity of a gas reservoir is proportional to its pressure. The discovery pressure for most reservoirs is found between the upper and lower limits1 shown in Fig. 1. Discovery pressures in excess of their hydraulic pressure gradient are attributed to the compaction of an isolated shale formation by the overburden.2.3 The lower limit of 0.433 psi/ft corresponds to the weight of the overburden, the gradient of the rock itself. Somewhere between

Jan 1, 1969

By B. E. Eakin, A. L. Lee, M. H. Gonzalez

Expeximental viscosity and density data of lour natural gases are presented for temperatures from WO to 340F and pressures from 100 to 8.000 psia. A correlation is also discussed and results reported. INTRODUCTION This investigation is one of several efforts by the authors to provide viscosity data for pure hydrocarbons and mixtures. Results of some pure hydrocarbons and their mix tures have been presented." Several predictive methods and correlations have also been reported.'," This paper presents the experimental viscosity and density data of four natural gases and confirms a correlation in a previous study. APPARATUS The viscometer was described previously. 2,6 A bank of stainless steel pycnometers was included to determine density in conjunction with the viscosity measurements (Fig. I). DATA AND MATERIALS The natural gases were furnished by the Atlantic Richfield Co. (Samples 1 and 21, the Continental Oil Co. (Sample 3) and the Pan American Petroleum Corp. (Sample 4). Table 1 shows the composition of these gases obtained by mass spectrometer analysis. The pure component viscosity data used in the correlation have been published."".'"', '" Experimental and calculated viscosity data for the natural gases are presented in Tables 2 through 5. CORRELATION Starling and Ellington"' have reported several semi-empirical expressions based to a certain extent on the theory of viscosity advanced by Born and Green.' The final expression presented by Starling and Ellingto'n is: i'. (micropoise) -- /t,, exp rX(7') p ..... (I) Eq. 1 was modified by Lee et 01." to represent mixture and pure component data simultaneously. This equation has the form: F = K exp [XpY] :.........(2) where K=- (7.77 i 0.0063M)Ti 122.4 + 12.9M +T .................(3) ' ' ' X= 2.57 + ---1914.5+ 0.0095M .... (4) Y = 1.11 -r 0.04 X ........(5) Over the pressure and temperature range studied in this investigation, this equation represents the data on methane, ethane, propane, n-butane and four methane-n-butane mixtures with a standard deviation* of ? 1.89 per cent. Density data by Sage and Lacey" were used to fit the equations. Eq. 2 gives reliable values of the viscosity of light hydrocarbons without prior knowledge of experimental viscosity values, but accurate density data must be available. Eq. 2 was used to predict viscosity values for the natural gases studied in this paper. The particular set of parameters contained in Eqs. 3 through 5 gave values that reproduced experimental data within 25 per cent. Needless to say, in using Eq. 2 better density values will give better viscosity results. However, those engineers who do not have accurate density values at hand or the time, or such facilities as a library, computer, etc. readily available, may feel uneasy in using Eq. 2. Therefore, the authors sought the easiest but by no means best density prediction method, which was reported by Kaf 30 years ago. A set of density values, which were calculated based on Kay's method, and a generalized compressibility factor chart were used to fit Eq. 2. The results are:

Jan 1, 1967

By G. Pizzi, G. M. Ciucci, G. L. Chierici

The use of the material balance equation to estimate the volume of hydrocarbons originally present in a reservoir, whose producing mechanism is partly due to water drive, has been discussed in the literature by several authors. There is no general agreement upon the possibility of obtaining reliable results by this method. Gas reservoirs in contact with an active aquifer are considered in this paper. Theoretical considerations based on the cybernetic principle of uncertainty (which states that the internal structure of a system cannot be uniquely determined from its observed external behavior) lead to the conclusion that the volume of gas originally present in a reservoir of this type cannot be uniquely determined from its past history. The range of values which encompasses the actual value of the reserves varies from case to case and must be determined by either numerical or analogical methods. Results obtained for six gas fields are reported. All these fields were produced with small fluctuations in their production rates, as is common practice for gas reservoirs; no gas storage fields were considered. Results obtained show that reserves values in a range of 1 to 2, associated with appropriate aquifers, allow the matching of the reservoir past history with mean-square deviations less than the experimental errors involved in pressure and produclion measurements. Similar results have been found in several other partial water drive gas reservoirs. From these results it is concluded that gas reserves cannot be uniquely determined from the past performance of partial water drive reservoirs, at least in cases where the reservoir has been submitted to small fluctuations in the production rates, and pressure data of normal accuracy are available. INTRODUCTION A number of authors have analyzed the problem of estimating the reserves originally present in a partial water drive reservoir from its past pressure-production performance. Literature which deals with this subject can be grouped, according to their conclusions, as follows. 1. A reliable value for the reserves can be obtained even if reservoir data (pressure and cumulative production) are affected by errors within the normal range.'-' 2. A reliable value for the reserves can be obtained only if reservoir data are very accurate'- or if past production performance has been subjected to abrupt variations in the production rate."' 3. No unique value for the reserves can be obtained from reservoir past production performance. This conclusion has been based upon theoretical considerations" and verified in several field cases. The purpose of this paper, which deals only with partial water drive gas reservoirs, is to test the above conclusions against actual field cases. Some theoretical considerations on this problem are also presented. THEORETICAL CONSIDERATIONS As the behavior of a gas reservoir communicating with an aquifer depends on both the aquifer and the reservoir characteristics, the physical system to be studied is the combined gas reservoir plus aquifer. The information which is available for studying the performance of such a system is the well production rates and bottom-hole pressures, all given as functions of time. The external behavior of the reservoir-aquifer system is therefore described by 2n input variables (Gp, Wp,) and n output variables (p,), n being the number of wells in the reservoir. In reservoir engineering it is common practice to consider the reservoir as a whole, disregarding the internal distribution of pressures and of producing wells. This practice is equivalent to substituting the above multi-variable system with a single-variable system, where the average reservoir pressure is the only output variable and the cumulative production G,(t) and W,(t) are the input variables. The internal structure of such a system, defined by initial gas reserves G, aquifer shape and dimensions, boundary conditions and petrophysical parameters distribution throughout the aquifer, is unknown. Therefore, from a cybernetic point of view the system is a "black-box It has been demonstrated" that for a blackbox the indetermination principle holds. Accordingly, the number of different internal structures (or set of parameters) which can account for the observed external behavior is infinite. As a consequence, the initial reserves cannot be uniquely determined from the reservoir past performance. When W,(t) is known, the determination of G from