Search Documents

By T. A. Pollard, C. R. Sandberg, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler', and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl' have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By C. R. Sandberg, T. A. Pollard, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler', and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl' have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By Francis Collins, E. R. Brownscombe

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By E. R. Brownscombe, Francis Collins

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By A. D. K. Laird, John A. Putnam

This paper describes the application of x-ray theory to design procedures in connection with fluid saturation determinations during fluid flow experiments with porous media. A reliable and rapid method for calibrating the x-ray apparatuy is described. Extension of the method to fluid saturation determinations in three-fluid systems is described. INTRODUCTION In rerearch on oil production problems a method is required which will give quickly the quantity of each component of a fluid flow system present at any cross-section of a porous medium. The sample of porous medium under investigation is usually referred to as a core. The ratio of the volume of one component to the total fluid volume is defined as the saturation of the porous medium by that component. This ratio is generally given as per cent saturation. Some means of measuring saturation which have received consideration include: electrical conductivity of the fluids;1,2 emissions from radioactive tracers dissolved in the fluids; the radioactivity of silver caused by reflection of neutrons from hydrogen atoms in the fluids;' the attenuation of a microwave beam. the diminution and phase shift of ultrasonic wave trains.4,5 and the reduction in intensity of x-ray beams in passing through the fluids. X-rays have already been used with some success. Since every material has a different power to absorb x-rays, the reduction in intensity of an x-ray beam as it passes through a core depends on the fluids present. The strength of the emergent beam can be found by converting its energy into a measurable form such as heat or ionic current. or by its effect on a photographic plate or fluorescent screen. The beam strengths could be interpreted as quantities of known fluids in the core if, previously, these beam strengths had been identified with a known combination of the same fluids. With some fluid cornbinations it might be desirable to dissolve powerful x-ray absorbing materials in one or more of the fluids, to increase the differences in the beam strengths for various fluid saturations. Boyer, Morgan and Muskat6 have described a method of measuring two component fluid saturation. One component was air or water; the other. minerat seal oil in which was dissolved 25 per cent by weight of iodobenzene to increase its absorbing power. The x-ray source was a tungsten target tube operated at 43 kv potential. The beam emerging from the core was measured as ionic current flowing across an air-filled ionization chamber by means of an amplifying circuit and galvanometer. Another portion of the beam from the x-ray tube was passed through a metal plate and measured in another ionization chamber. This portion, called the monitor beam, was used as an indication of the performance of the x-ray tube. The galvanometer readings were calibrated against air-oil core saturations, gravimetrically determined. The method was apparently established by experimental means. In the present investigation the available theory of x-radia-tion was surveyed with a view to extending the usefulness of the method and to developing design procedures for its application to measurement of fluid saturation in porous media. Application of the theory permits prediction of relative meter readings to be expected for any combination of porous matrix, various saturating fluids and auxiliary filtering media. It is thus possible to calibrate the equipment in terms of fluid saturation by an indirect but rapid technique. The results of calculations based on x-ray theory indicate. and results of the saturation calibration technique confirm. that a valid measurement of the saturation of the core can be made for any two components and in some cases for three components. THEORY The strength of an x-ray beam, after it has passed through a distance. 1, of matter of density, p, and mass absorption coefficient, µ at a given wavelength, A, may be expressed by the absorption formula I = I0 e ...........(1) where I, represents the intensity of the incident x-ray beam and I is the intensity of the emergent beam. The expression e is called the transmission factor of the material. The variation of I,, with wavelength depends upon the materials through which the x-ray beam has previously passed and upon the spectral distribution of energy at the source of the x-radiation. A group of curves. called spectra. which show the variation of intensity with wavelength and x-ray tube voltage are given in Fig. 1. These curves represent the general radiation from a tungsten target tube. When the tube voltage is greater than 69.3 kv, the characteristic radiation of the tungsten is emitted and is superposed on the general radiation. At a given voltage the minimum wavelength A,,,,, at which energy can be emitted by an x-ray tube is given by the formula 12,340 xml. = ——..........(2) volts where A,.,,.. is in Angstrom units. The wavelength at which the spectra have maximum intensity a1so decreases with increasing x-ray tube voltaue. The area under each curve represents to an arbitrarv scale the total energy emerging from the x-ray tube for that voltage. The variation of µ with wavelength has been determined for many substances and may be found in such references as those by Compton and Allison7 and by Hodgman.8 The phenomenon of absorption is composed chiefly of the capture of photons by the atoms of the absorbing material with associated displacement of electrons, and of the scattering, or the deflection, of the photons by the atoms. Curves of these mass absorption coefficients show jump discontinuities. or absorption edges. at wavelengths which are short enough for the photons,

Jan 1, 1951

By John A. Putnam, A. D. K. Laird

This paper describes the application of x-ray theory to design procedures in connection with fluid saturation determinations during fluid flow experiments with porous media. A reliable and rapid method for calibrating the x-ray apparatuy is described. Extension of the method to fluid saturation determinations in three-fluid systems is described. INTRODUCTION In rerearch on oil production problems a method is required which will give quickly the quantity of each component of a fluid flow system present at any cross-section of a porous medium. The sample of porous medium under investigation is usually referred to as a core. The ratio of the volume of one component to the total fluid volume is defined as the saturation of the porous medium by that component. This ratio is generally given as per cent saturation. Some means of measuring saturation which have received consideration include: electrical conductivity of the fluids;1,2 emissions from radioactive tracers dissolved in the fluids; the radioactivity of silver caused by reflection of neutrons from hydrogen atoms in the fluids;' the attenuation of a microwave beam. the diminution and phase shift of ultrasonic wave trains.4,5 and the reduction in intensity of x-ray beams in passing through the fluids. X-rays have already been used with some success. Since every material has a different power to absorb x-rays, the reduction in intensity of an x-ray beam as it passes through a core depends on the fluids present. The strength of the emergent beam can be found by converting its energy into a measurable form such as heat or ionic current. or by its effect on a photographic plate or fluorescent screen. The beam strengths could be interpreted as quantities of known fluids in the core if, previously, these beam strengths had been identified with a known combination of the same fluids. With some fluid cornbinations it might be desirable to dissolve powerful x-ray absorbing materials in one or more of the fluids, to increase the differences in the beam strengths for various fluid saturations. Boyer, Morgan and Muskat6 have described a method of measuring two component fluid saturation. One component was air or water; the other. minerat seal oil in which was dissolved 25 per cent by weight of iodobenzene to increase its absorbing power. The x-ray source was a tungsten target tube operated at 43 kv potential. The beam emerging from the core was measured as ionic current flowing across an air-filled ionization chamber by means of an amplifying circuit and galvanometer. Another portion of the beam from the x-ray tube was passed through a metal plate and measured in another ionization chamber. This portion, called the monitor beam, was used as an indication of the performance of the x-ray tube. The galvanometer readings were calibrated against air-oil core saturations, gravimetrically determined. The method was apparently established by experimental means. In the present investigation the available theory of x-radia-tion was surveyed with a view to extending the usefulness of the method and to developing design procedures for its application to measurement of fluid saturation in porous media. Application of the theory permits prediction of relative meter readings to be expected for any combination of porous matrix, various saturating fluids and auxiliary filtering media. It is thus possible to calibrate the equipment in terms of fluid saturation by an indirect but rapid technique. The results of calculations based on x-ray theory indicate. and results of the saturation calibration technique confirm. that a valid measurement of the saturation of the core can be made for any two components and in some cases for three components. THEORY The strength of an x-ray beam, after it has passed through a distance. 1, of matter of density, p, and mass absorption coefficient, µ at a given wavelength, A, may be expressed by the absorption formula I = I0 e ...........(1) where I, represents the intensity of the incident x-ray beam and I is the intensity of the emergent beam. The expression e is called the transmission factor of the material. The variation of I,, with wavelength depends upon the materials through which the x-ray beam has previously passed and upon the spectral distribution of energy at the source of the x-radiation. A group of curves. called spectra. which show the variation of intensity with wavelength and x-ray tube voltage are given in Fig. 1. These curves represent the general radiation from a tungsten target tube. When the tube voltage is greater than 69.3 kv, the characteristic radiation of the tungsten is emitted and is superposed on the general radiation. At a given voltage the minimum wavelength A,,,,, at which energy can be emitted by an x-ray tube is given by the formula 12,340 xml. = ——..........(2) volts where A,.,,.. is in Angstrom units. The wavelength at which the spectra have maximum intensity a1so decreases with increasing x-ray tube voltaue. The area under each curve represents to an arbitrarv scale the total energy emerging from the x-ray tube for that voltage. The variation of µ with wavelength has been determined for many substances and may be found in such references as those by Compton and Allison7 and by Hodgman.8 The phenomenon of absorption is composed chiefly of the capture of photons by the atoms of the absorbing material with associated displacement of electrons, and of the scattering, or the deflection, of the photons by the atoms. Curves of these mass absorption coefficients show jump discontinuities. or absorption edges. at wavelengths which are short enough for the photons,

Jan 1, 1951

By H. T. Kennedy, J. E. Little

An empirical equation for the prediction of the viscosity of several pure paraffin hydrocarbons and nitrogen is presented. It involves temperature, pressure and six constants of the material, and it applies reliably to both liquids and gases. The equation is similar in form to van der Waal's equation of state. For the paraffin hydrocarbons methane through r-hexane and nitrogen, an average absolute deviation of 1.9 percent was obtained on 1,006 data points described in the literature by 14 authors. When this equation is extended to complex, liquid hydrocarbon mixtures, a correlation was obtained with an average absolute deviation of 9.9 percent. INTRODUCTION Equations describing the flow of gas and liquid through porous media contain the viscosity coefficient of the fluid. If other pertinent variables remain constant, the volume rate of flow is inversely proportional to this coefficient. In dealing with condensate fluids and volatile oils, however, the compositional effects resulting from changing pressure materially affect the viscosity. The effect of compositional changes also may be significant in certain secondary recovery or pressure maintenance processes, notably miscible displacement or gas injection. Early attempts to describe the performance of reservoirs utilized a volumetric material balance method wherein gas and liquid in the reservoir were identified as produced gas and liquid at the surface. This method of analysis proved adequate for reservoirs at moderate temperature and pressure that contained gas with moderately low amounts of condensable materials. The volumetric material balance procedures for "black oil" reservoirs leave much to be desired when applied to condensate and volatile oil reservoirs because phase behavior and compositional changes are relatively more important in these cases. The alternative is a compositional material balance, which in turn, requires a correlation of properties of the reservoir fluid with composition. This paper supplies this correlation in regard to viscosity, for reservoir crude oils. REVIEW OF LITERATURE The literature contains many empirical equations describing the effects of composition, temperature and pressure on the viscosities of pure liquids and binary liquid mixtures. However, the applicability of a majority of these equations is limited to very low pressures and to a small number of systems. Most of them, when applied to complex hydrocarbon systems, are of little value. The lack of utility of the majority of equations results from the fact that they were developed to show the separate effect of temperature, pressure or composition on viscosity, but not to predict the viscosity as a function of all three variables. And with the few exceptions noted below, they were developed to apply to much simpler systems than hydrocarbon mixtures. Eakin1 developed an equation for the viscosity of pure light hydrocarbons applicable in the gas, liquid and dense fluid regions, and in its most general form is written as µ - µo =A(ebP - e-cpk) The pure component viscosity at atmospheric pressure and temperature T, µo, was calculated using the Sutherland equation: µo = BT3/2/T + S The constants in Eakin's equation permitted the term e-CPk to become quite small at high densities (above 0.215 gm/cc), and the general equation would simplify to

Jan 1, 1969

By A. S. Emanuel

A correlation has been developed for calculating the phase behavior of, gas-condensate systems at reservoir conditions. The correlation is based on the principle of corresponding states and has been coded for an IBM 7094. Experimental K-values were determined for several gas-condensate systems at reservoir conditions to evaluate various semi-empirical parameters of the correlation. The approximate range of application of the correlation is 150 to 300F and 1,500 to 6,000 psi. INTRODUCTION The rapid development of digital computers during the past several years has made feasible the calculation of hydrocarbon phase behavior by methods based on rigorous thermodynamic principles. Good correlationsl have been developed for low to moderate pressures, but these techniques have not yet been extended successfully to reservoir fluids at high pressures. Consequently, the determination of phase behavior of oil and gas systems at reservoir conditions is still based almost entirely on generalized data correlations or on experimental studies of the fluid in question. While these methods have been used successfully many times, they do have inherent limitations that restrict their applicability. Generalized correlations, such as the NGShlA K-charts, are limited to the range of pressure, temperature and components for which the data were determined. The accuracy of these correlations is often questionable because the effect of total system composition is not well defined. Experimental studies offer a reliable method for determining phase behavior, but usually the studies are costly and time consuming. Recently, Leland and co-workers presented a new approach to calculating phase behavior from the principles of corresponding states.6-8 Corresponding states methods determine the thermodynamics properties of a given system by comparison with a reference substance whose properties are known. The accuracy of this approach depends on close chemical and structural similarity between the reference substance and the system in question and between components within the system itself. For high accuracy, it is usually necessary to correct for chemical and structural dissimilarities. In principle, however, the corresponding states method should be no less accurate at high pressures than at low pressures, provided reference substance properties are known. This paper describes an empirical modification of the basic correlation proposed by Leland, et al.6 for the specific purpose of calculating the phase behavior of gas-condensate fluids at reservoir conditions. The modified correlation, which has been programmed for an IBM 7094, may be used for either approximate or precise determination of fluid behavior depending on the amount of analytical and, experimental data available for the system. BASIC THEORY The basic theory of the corresponding states phase equilibria correlation was first published by Leland, Chappelear and Gamson.6 Subsequently, Leland, Chappelear and Leach published methods for improving the accuracy of the original theory.7-9 The aim of the correlation is to calculate the K-value of each component of a given system as a function of pressure, temperature, and over-all composition, where Once the K-values are known, the phase behavior may be determined directly by an appropriate flash calculation. The basic equation for calculating component K-values was taken from the work of Joffe. 4 For any component i of a mixture, the K-value is given by

Jan 1, 1969

By A. L. Barnes, A. M. Rowe

A heat transfer study was made of hot gas injection into oil shale through wells interconnected by vertical fractures. This analysis involved the simultaneous numerical solution of a nonlinear, second-order partial differential equation that describes two-dimensional conduction heat transfer in oil shale and a nonlinear first-order partial differential equation that describes convection heat transfer in the fractures. Three nonlinear, temperature-dependent coefficients were used in this work; they are thermal conductivity, thermal capacity and retorting endothermic heat losses of oil shale. Vertical fractures were considered to be of finite height. Although vertical conduction heat transfer was not considered, an estimate of the error resulting from this limitation was made. The effects of injected gas temperature, injection rate, system geometry, cyclic injection and time upon retorting efficiency were investigated. Results from this study show that the rate of retorting oil shale is a direct function of both injection temperature and rate, and the theoretical producing air-oil ratio:(AOR) is an inverse function of temperature. Retorting rates are constant until "breakthrough" of the 700 F isotherm at the producing well, assuming constant injection parameters. Retorting rates for bounded systems are higher than the analogous unbounded systems and likewise AOR's are less. The use of an alternating injection-soak routine with high injection rates is less efficient than continuous injection at lower rates. These results indicate that injection temperatures on the order of 2000F or greater may give theoretical AOR's in the economic range. INTRODUCTION Over half of the known oil shale reserves are located in the U.S., and most of them lie in the Piceance Creek basin of Western Colorado. The Colorado oil shale outcrops on the edges of the Piceance Creek Basin. At the outcrops the shale beds are relatively thin, from 25 to 50 ft thick. In the center of the basin the oil shale is as great as 2,000 ft thick and is covered with 1,000 ft of overburden. It has been estimated that there are over 1,000 billion bbl of oil in shales having an oil content over 15 gal/ton in this basin. Oil shale does not contain free oil but an organic matter called kerogen. Kerogen yields petroleum hydrocarbons by destructive distillation. It must be heated to approximately 700F, at which temperature it decomposes into shale oil, gases and coke. The U.S. Bureau of Mines and, more recently, oil companies have conducted considerable research on surface retorting methods to economically recover oil from this shale. Another approach to exploit the oil shale deposits, in particular that portion having 1,000 ft of overburden, is to retort the oil shale in place and produce the liquid and gaseous hydrocarbons through wells drilled into the shale. Some research has been done on this approach.l,2 There are several variations to the in situ retorting approach. These variations fall into one of two groups, depending upon the geometry of the system: (1) retorting in a highly fractured or broken up matrix; (2) retorting from single fractures between production and injection wells. The latter is the group studied. Several investigators,3-6 using various assumptions, have studied flow of heat through horizontal systems. The objective of this work was to make a heat transfer study of in situ retorting oil shale by hot gas injection through wells interconnected by single vertical fractures of finite height. The oil shale thermal conductivity, thermal capacity and retorting endothermic heat losses were considered to be functions of temperature. A knowledge of the

Jan 1, 1969

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By F. G. Miller, T. D. Mueller, L. B. Davidson

A mathematical model of reservoir behavior during one cycle of the intermittent steam injection process is presented in this paper. A cycle of the process is considered to be composed of three stages: injection, condensation and production. The model concentrates on behavior during the production stage. The influence of several parameters on production improvement is described, using the results obtained with a numerical equivalent of the model. Reservoir system parameters considered are thickness, reservoir thermal conductivity, oil viscosity, solution gas-oil ratio and the change in oil effective permeability with temperature. Operating parameters investigated are heat injection rate, steam temperature and cumulative injection time. In the face of considerable uncertainty, the model appears to give a fair prediction of the behavior observed during an actual test. This makes the model a useful means of selecting projects. It is shown that the model may also be utilized to select an optimum set of operating conditions for a particular project. Results of this study show the extent to which production improvement depends on each of the parameters investigated and give an indication of the range of system conditions within which the process would be effective. INTRODUCTION In recent years the intermittent steam injection or steam soak process of oil recovery has become a popular method of improving recovery in reservoirs bearing high-viscosity hydrocarbons. Production improvement due to the process results from heating the reservoir near the wellbore. The process is applied on an individual well basis and is executed in a series of cycles, each consisting of three stages: injection, condensation (the soaking period) and production. During the injection stage, high-pressure steam is injected into the reservoir surrounding the well-bore. This study assumes that injection is carried out over the entire vertical section of the reservoir rock that can be exposed to fluid injection. The condition of steam entering the reservoir is considered to fall within the following limits: (1) steam pressures less than 1,000 psig, (2) steam temperatures between 300 and 700F and (3) steam quality less than 90 percent. The last restriction is due to equipment limitations. Normally, injection re maintained for less than 3 weeks. After termination of injection, the condensation stage begins. The well is shut in for several days, supposedly to make certain that the injected steam condenses. This condensation would insure the retention of injected heat. The production stage is initiated when the well is opened again. Reservoir response during this stage determines the effectiveness of the process (response is defined as the oil production rates, the water production rates and the wellbore temperatures observed during the entire production stage). The immediate response of the reservoir is sometimes quite good. Sets of case histories by Longl and by Owens and Suter2 give some typical figures for response. Duration of the production stage is equal to the time needed for the oil production rate to decline to some predetermined value. This usually requires 4 to 6 months. Termination of the production stage coincides with the start of the next cycle. The process is continued, cycle by cycle, until it becomes unprofitable to continue. One simple model2 and two more sophisticated models3,4 of the process have been presented. These are based on fundamentally different assumptions than those used in this investigation. Although the three models might be modified to fit the assumptions used in this study, it is felt that greater rigor would result by developing a

By B. H. Caudle, L. G. Davies, I. H. Silberberg

This paper presents a method of predicting the recovery and performance of a five-spot steam injection project, in which a realistic approach to pattern sweepout efficiencies is made. Published methods for radial systems were modified for the five-spot pattern by approximating the stream lines with straight lines radiating from the injection well and then converging to the producing well. In each radial segment, the position of the steam front and the temperature profile ahead of the steam front were determined by heat balance equations, which included an estimation of heat losses to surrounding formations. The location of the saturations behind the cold water front was determined from a Buckley-Leverett solution to the material balance equation. Results from this program show steamflood recovery in a five-spot pattern to be considerably less than that predicted for true linear or radial flow systems. For a specific reservoir containing 900-cp oil, a steamflood in a purely radial flow system was predicted to recover more than 75 percent of the original oil in place when 2 PV of water had been injected as steam. A five-spot steamflood with otherwise identical properties was predicted to recover 10 percent of the original oil in place when 0.15 PV of water had been injected as steam and to recover almost no oil thereafter. A cold water five-spot flood in this system was predicted to recover approximately W percent of the oil in place with I PV of water injected. For a five-spot pattern in an example reservoir with 10-cp oil, steam injection similarly showed lower ultimate recovery than water injection but no improvement in recovery rule. Introduction The thermal recovery method considered in this study is steam injection in a five-spot pattern. Pattern steam injection has been given little attention in the past, possibly because of some rather obvious disadvantages. Unless steam temperature is maintained throughout the entire swept area, the process will revert to simple waterflood with all heat being lost prior to reaching the oil bank. Also, the oil viscosity reduction takes place near the steam front and not around the producing wellbore, so that low producing rates must still be endured. This is the reason for the success of the steam stimulation, or cyclic injection, in which the fluids are produced from the same well used for injection. On the other hand, steamflooding can offer some advantages. Oil displacement is by four methods: (1) mechanical displacement by the condensed water, (2) viscosity reduction of the crude oil, (3) swelling of the crude oil, and (4) distillation of the crude oil in the steam zone. Although dependent upon the crude, laboratory experiments have shown displacements up to 80 percent by steamflooding. In addition, steam is a good heat transport medium, since it is cheap and has a high heat content. Previous Investigations Previous publications have reported methods for predicting recovery in a steamflood for linear and radial systems. Since no pattern sweep efficiency is taken into consideration, the recovery even from a radial system must be greater than from more realistic geometries such as the very common five-spot patterns. Probably the broadest coverage is given by Willman who reported experimental results and offered a method of predicting recovery for a radial flow system. They concluded that both hot water and steam injection recover more oil than an ordinary waterflood, and that steam injection could yield recoveries "as much as 100 percent greater than by water flood." They also concluded that both the heat requirements for a reservoir and the residual oil remaining in the reservoir after steam injection were independent of the amount of oil originally in place, that short exploitation times were desirable, and that a high percentage of net sand in the reservoir with a high initial oil saturation was desirable. The method assumes that the flood occurs in three concentric cylindrical zones: (1) an inner steam zone, (2) a central hot water zone, and (3) an outer cold waterflood zone. Displacement in the steam zone is based on laboratory-determined residual oil saturations while the hot water and cold water zones use the conventional Buckley-Leverett equations. Although the results shown by Willman et al. are for a radial system flowing out from the well to an assumed external circular boundary, their equations did

Jan 1, 1969

By R. H. Jacoby, V. J. Berry, R. C. Koeller

The experimental phase behavior of several field gas-condensate systems, one field volatile oil system, and a series of synthetic systems having gas-oil ratios from 2,000 to 20,000 scf/bbl stock tank oil was measured. The data were used to calculate the depletion performance of the field systems at their respective reservoir temperatures and of the synthetic systems at various temperatures. The results of the performance calculations were used to prepare correlations of total stock tank oil and separator gas in-place, and the amount of each recovered by primary depletion. Data needed to use the correlations are initial gas-oil ratio, initial tank oil gravity, reservoir temperature and reservoir pressure. The correlations presented should be useful for estiinating recovery from reservoirs producing volatile oils or rich gas condensates having gas-oil ratios from about 2,000 to 30,000 scf/bbl of stock tank oil. INTRODUCTION Prediction of depletion performance for gas-condensate reservoirs usually requires extensive laboratory work or lengthy phase-behavior calculations. For many small gas-condensate bearing zones the expense of such work is often not justified. Furthermore, where large holdings are involved, it is useful to have estimates of future oil and gas production in advance of the results of a laboratory analysis, even if such work is ultimately planned. For these reasons, it would be of value to have a method for rapidly estimating the depletion performance of gas-condensate reservoirs from the type of field data usually available. Few data on the reservoir depletion performance are available for the purpose of developing correlations. The increasing number of gas-condensate discoveries and the advent of high-speed computers have stimulated the development of methods for predicting performance. Such calculated performance data have been shown to be reliable and are becoming accepted as a basis for reserve estimates and for engineering studies. It therefore seems satisfactory to use results of performance predictions themselves for the purpose of developing the desired correlations. Such an approach is described in this paper. Several rich natural gas-condensate systems (gas-oil ratios of from about 3,600 to 60,000 scf/bbl) and one natural volatile oil system were studied (2,363 scf/bbl). In addition, to provide a more systematic basis for correlating purposes, a series of related, synthetic reservoir fluid mixtures was also studied (GOR's from about 2,000 to 25,000 scf/bbl). The purpose of this paper is to show how the recovery of oil and gas by primary depletion may be estimated from surface gas-oil ratio, reservoir temperature and pressure, and tank oil gravity. procedure Experiments Companion samples of gas and oil were collected from the high pressure separators in nine different fields at the time of their discovery. These samples were each recombined in accordance with the respective producing GOR to reproduce the reservoir fluids. Eight of these fluids, designated F-2 through F-9, were rich gas-condensate systems at the respective reser-

By M. R. Tek, F. F. Craig, J. O. Wilkes, C. S. Goddin

The waterflood performance of a water-wet, stratified system with crossflow is computed by a finite difference procedure. The effects of five dimensionless parameters on tile oil displacement efficiency, water saturation con-tour.7 and crossflow rates are evaluated in the absence of gravity forces. Crossflow due to viscous and capillary forces is shown to exert a significant effect on oil recovery in a field-scale model of a two-layered. water-wet sandstone reservoir. The crossflow is at a maximum in the vicinity of the front advancing in the more permeable layer. Under favorable mobility ratio conditions, the comparted oil recovery with crossflow always is interrnerliate between that predicted for a uniform reservoir and that for a layered reservoir with no crossflow. INTRODUCTION The important erects of reservoir heterogeneity on waterflood performance are commanding increased attention in the technical literature. Much of this attention is centered on two categories of layered reservoirs: those in which layers are non-communicating and those in which crossflow of fluids occurs between the layers. In the first category, the reservoir is assumed to consist of discrete layers, each uniform within itself and differing from the others only in such properties as thickness, porosity and absolute permeability. The performance within each layer is calculated by one-dimensional flow theory, and the performance of the total reservoir is obtained by summing individual layer performances. Capillary and gravity effects usually are not considered. Representative publications dealing with thi5 type of reservoir are those of Stiles,' Dykstra and Parsons,' Hiatt,3 Warren and Cosgrove' and Higgins and Leighton." Prediction of performance for reservoirs in the second category is considerably more difficult since viscous, capillary and gravitational forces all play important roles in causing crossflow between layers. A number of authors have investigated the simpler problem of two-dimensional displacement flow in a stratified system with a mobility ratio of unity and negligible capillary and gravity effects.'; Others have considered two-dimensional, non-steady-state flow of a single, slightly compressible fluid in a stratified reservoir. A limited number of laboratory oil displacement tests in layered models with crossflow have been reported. Miscible floods (with resultant zero capillary forces) in layered five-spot models were conducted by Dyes and Braun," who studied the effect of mobility ratio with zero gravity forces, and by Craig et al. 12 who studied the effect of gravity forces at constant mobility ratio. Waterfloods in layered five-spot models (with cross-tlow due to capillary, viscous and gravity forces) were conducted by Gaucher and Lindley,"' who showed the effect of gravity forces in causing underrunning of the injected water and by Carpenter, Bail and Bobek, 14 who demonstrated the reliability of Rapoport's" dimension-less parameters for scaling layered systems. Waterfloods in rectangular layered models were conducted by Richardson and Perkins."' who investigated the effect of velocity at constant mobility ratio and with zero gravity forces, and by Hutchinson," who studied the effects of varying mobility, layer permeability and layer thickness ratios. The differential equations which rigorously describe waterflooding in a heterogenous porous medium are non-linear and do not facilitate analytical solution. By using finite difference approximations it is possible to obtain a solution to any desired degree of accuracy. Such a solution, using an alternating direction implicit procedure (ADJP), is described by Douglas, Peaceman and Rachford.18 In the present study, a computer program using ADJP explores systematically the effects of important parameters on waterflood performance of a two-dimensional, two-layered, field-scale model of a water-wet sandstone system. Particular attention is given to evaluation of the water saturation contours and crossflow rates at the interface between layers to gain improved understanding of the crossflow mechanism. PROCEDURE BASIC: FLOW EQUATIONS The basic flow equations for two-dimensional, two-phase, immiscible, incon~pressible flow in a porous medium are:

Jan 1, 1967

By C. Perez-Rosales

A simplified method for evaluating specific surface of porous media is described. The theory takes as its point of departure a previous statistical method. Due to its simplicity, it is considered that the proposed method has certain advantages over methods currently used. By making some simplifications, the procedure for estimating specific surface is reduced to superposing, on an enlarged reproduction of a section of the sample to be analyzed, an evenly spaced grid and counting the number of intersections of the grid lines with the perimeter of pores. An indirect test of the theory is presented and the advantages of the method are discussed. INTRODUCTION The specific surface of a porous medium, defined as the surface area of the pores per unit bulk volume, is an important parameter with regard to flow of fluids through porous materials. In particular, it has long been recognized that specific surface is closely related to permeability. Since the internal surface of any porous material is extremely complex, the measurement of specific surface can only be accomplished by indirect means. It is considered that the best method' in current use for determining specific surface is the statistical method of Chalkley, Cornfield and Parka who originally applied it for estimating volume-surface ratios of morphologic components, such as cells or nuclei. This method will be outlined later. A method which has been used extensively in determining specific surface is based on the adsorption of a vapor by the internal surface of porous media. However, this method does not yield the surface pertinent to fluid flow, since it includes the tiny molecular interstices which do not influence the permeability of the porous material.' This method has been employed by Brooks and Purcell measure internal surface area of sedimentary rocks. The Kozeny equation, which relates the specific surface to the permeability of a porous medium, has found wide use in determination of specific surface. However, the KOzeny equation is not strictly correct. Large discrepancies have been found between theory and experiments.' Consequently. the specific surface estimated by this method is subject to serious uncertainties. Finally, there are other methods for specific surface measurements which are based upon a variety of physical and chemical properties: the heat conduction by a gas, ionic adsorption phenomena and direct chemical reactions of the surface with other substances. This paper describes a simplified statistical method for determining specific surface of porous materials. It is considered that, due to its simplicity, it has certain advantages over the methods presently used. For its validity, it is required that the porous structure of the material under study be of a spatially random nature. This is not a serious limitation since the naturally occurring porous media have this characteristic. THEORY The method proposed in this paper takes as its point of departure the interesting result obtained by Chalkley et al The method shows that if a line of length r is dropped repeatedly at random over a section of a porous medium, and counts are made of the number of times the two end points fall within pore areas (denoted by h for "hits") and the number of times the line intersects the perimeter of pores (denoted by c for "cuts"), then in a very large number of throws where S is the internal surface of the sample (surface of pores) and V, the volume of the pores (Fig. 1). Since the specific surface S, and the porosity 4 are given, respectively, by S. = Sl V, and 4 = V,l V,, where V, is the total volume of the sample, it follows from Eq. 1

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By C. S. Matthews, G. L. Stegemeier

In one field in South Texas, approximately 72 per cent of the pressure build-up results show a characteris-i.rtic "hump" (i.e., the pressure builds up and then falls off) which makes interpretation by. standard methods impossible Correlation of size and time of hump with formation permeability. well productivitv index, and method of completion led to the tentativc, conclusion that the humps were caused by segregation of gas and oil in the wellbore after closing-in"'. This conclusion was confirmed by performance of simple laboratory bubble-rise experiments, by theoretical bubble-rise time calculations, and by a detailed calculation of PVT behavior in the well of a particular well Oil which accurate Surface and bottom-hole pressure measure-ments were untile. The hump behavior has since been found to occur in many other fields. The cause, however, is not the mine in all cases. In some of these the hump is trace-ahie to leaks in the tnbing which allow influx from the annulus after closing-in. In other cases the hump is traceable to leaks in the device separating pay horizons in dually completed wells. It is concluded that the recording of both surface and bottom-hole pressures is desirable in wells which show an anomalous build-up behavior. A number of field examples is discussed where use of both sets of measurements enables the cause for anomalous behavior to he found, (and a reasonable interpreation of bottom-hole pressure to be made. INTRODUCTION Theoretically the pressure build-up in an infinite reser-voir should be a linear function of In [(t + at)/nere t is the production time and At a the closed-in time. Some of the variations from this behavior arc well known, such as the curved portion immediately after shut-in which results from after-production and skin effect, and the flattened end portion which results from boundary effects in a limited reservoir. The effect of stratified producing zones and irregular geometrical drainage patterns may also contribute unusual characteristics to build-ups. However, the effect of still another phenomenon—that of movement of fluids within the wellbore—has been neglected in most build-up studies to date. Fluid movements within the wellbore after shut-in occur as a result of after-production, packer failures, leaks in the casing or tubing, or buoyancy of the gas phase when both gas and liquid are present in the well. Each of these movements can influence the pressure build-up, sometimes sufficiently to negate use of the data for computing permeability or stabilized build-up pressure. Examples of each of these phenomena are discussed below. FIELD OBSERVATIONS ON HUMP BUILD-UPS In approximately 75 per cent of the wells in a medium-sized field in South Texas the pressure build-up curves rise to 3 maximum and then decline to what appears to be a stabilized reservoir pressure as shown in Fig. I. Similar behavior has also been observed in wells in other fields in Texas and Louisian:,. Because of the possibility of such behavior being a mechanical failure of the pressure bomb. the bombs were run in tanden into one well which had been known to exhibit this behavior. When both bombs recorded identical curves with the characteristic hump, it was inferred that the unusual build-up is a result of well or forma-

By G. B. Thomas

Several methods of analyzing pressure build-up data in wells have been presented by various authors. This paper reviews the theory and method of D. R. Horner and presents example calculations performed on data obtained by testing several different types of wells. These calculations include, (1) graphical estimation of final static pressure, (2) determination of the productive capacity of the pay away from the well bore and. (3) the degree to which the formation adjacent to the well bore has been damaged by completion or other causes. The methods of testing and precautions which should be taken to assure the best data possible are discussed. Limitations and reliability of calculated results are also treated. INTRODUCTION Pressure testing of wells is generally limited to the determination of producing and static mean formation pressures. The so-called "static" pressure determination.. along with PVT, electric log and production data, enable the reservoir engineer to determine, within reasonable limits, the drive mechanism of the reservoir and in some cases. the amount of edge water encroachment. Producing pressure tests enable calculation of productivity indices and allow the engineer to plan the systematic production of a pool for optimum conservation of suhsurface energy. The radial flow formula advanced by Muskat' has been based on the assumption of incompressible radial fluid flow. It has been known that reservoir fluids do not behave in an ideally incompressible manner. For example, incompressible flow theory indicates a simple logarithmic relationship between the difference of the instantaneous and static well pressures when plotted against time. 'The latter stages of this type of plot of pressure build-up data generally show a marked devia-, tion from the earlier straight line trend. which deviation may be shown to be due to the compressible flow of fluids toward the well bore. Shut-in times of 24, 48. 72. or at most 96 hours are currently in wide use for determining so-called "static" reservoir pressure. Due to the continuation of compressible flow of fluids into the well bore long after this arbitrarily taken shut-in time, the determination of static pressures has almost invariably resulted in lower than equilibrium values. Materials balance calculations made early in the life of a reservoir often result in a calculated reserve which later observations prove to be too low. Failure to obtain reliable "static" reservoir pressures within the prescribed 24 or 48-hour build-up period has undoubtedly been a major factor in obtaining these low estimates. Comparison of theoretical and actual productivity performance has indicated that formation damage and not being able to attain true static pressures have been partially responsible for the observed discrepancies. Research into the theory of compressible flow behavior has resulted in methods of applying these theories to the testing of wells. Application of the developed theories to pressure build-up performance indicates that in many cases the following can be estimated. 1. The static reservoir pressure. 2. The in-place effective formation permeability away from the well bore. 3. The degree or extent by which the formation has been damaged adjacent to the well bore, either through completion methods or subsequent damage due to fluid entry. REVIEW OF THEORY Horner2 has shown that the pressure build-up within a "point source" well is approximated by the following formula: Equation (1) is the approximate "point source" solution to the radial compressible flow equation advanced by Muskat. The solution assumes the following: 1. A point source well is producing at constant rate from the center of an infinite reservoir with a constant pressure at its external boundary. 2. The fluid flowing is present in one phase only. 3. The compressibility and absolute viscosity of the fluid remain essentially constant over the range of temperature and pressure variation encountered. 4. The well i.; shut-in at the sand face arid there is 110 after production into the well bore. 5. The formation permeability is homogeneous in the direction of flow. From Equation (1) it can be seen that if an ideal well were shut-in while producing from a reservoir under the conditions assumed, the pressure build-up would be a logarithmic func-

Jan 1, 1953

By C. R. Dodson, D. Goodwill, E. H. Mayer

The paper summarizes the historical background of the pressure-volume-temperature analyses of reservoir fluids, the errors involved in both the sampling and testing of reservoir fluids, the type of information required of a PVT determination, and the field conditions that limit the application of any one analysis. Particular emphasis is placed on the necessity for approximating as closely as possible the liberation sequence occurring in the producing formation, flow string, and surface separators. A combined differential and flash or "composite" liberation is suggested as the best means of approximating this liberation sequence. HISTORICAL BACKGROUND Petroleum reservoir engineering commonly is considered to be one of the newest fields of petroleum science, yet much work of a theoretical and intuitive nature was done many years before the modern techniques of reservoir analyses were developed. The period 1910-1924 saw considerable work in the field of reservoir behavior done by the U. S. Bureau of Mines. This work:.' although entirely theoretical, pointed out the importance of gas in the recovery of oil from the reservoir. The statement by J. 0. Lewis that 20 per cent or less of the oil originally in place in a pool was recovered from the ground under uncurtailed production conditions caused some operators to re-examine and modify their production methods. However, though these works showed methods of production that have subsequently resulted in greater oil recoveries, they in general went unnoticed, even though there was a fear that the nation's petroleum resources were being exhausted. The period 1924-1933 saw the industry take considerably more interest in reservoir behavior because of an unfounded fear of its inability to replace oil reserves, the possibility of government regulation, and the energy of one man, Henry L. Doherty. Doherty aroused heated discussions in the industry concerning conservation. To prove or disprove his theories the first experimental work on the reservoir behavior of petroleum was undertaken. The papers of Dow and Calkin,' Beecher and Parkhurst," and Mills and Heitheckere are classics although their experimental procedures were crude. These papers proved that the properties of petroleum in the reservoir are quite dissimilar to the properties measured at the well head. Gas dissolved in the oil phase was recognized as having considerable importance as had been predicted by the earlier theorists. Slightly later Miller and his co-workers assembled the first text correlating all the known data on the reservoir behavior of petroleum and used, these data to show the economic value of conservation. Quantitative application of pressure-volume-temperature analyses of reservoir fluids was given in the paper of Coleman, Wilde, and Moore"; this work showed that with sufficient laboratory and production data prediction of reservoir behavior was possible. The refined sampling and labora-

Jan 1, 1953

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969