Search Documents

By R. A. Morse, P. M. Bridges, L. E. Wilsey, Howard N. Hall, P. L. Terwilliger

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By L. E. Wilsey, R. A. Morse, P. L. Terwilliger, Howard N. Hall, P. M. Bridges

Theoretical and experimental investigations of a constant pressure gravity drainage system are reported. Experimental data are presented to show that recovery to gas breakthrough by gravity drainage is inversely proportional to rate. The gravity drainage reference rate, which is numerically equal to the so-called "maximum theoretical rate of gravity drainage" is shown to have no particular significance from a recovery standpoint. Before this rate can be used as a basis of comparison for recoveries, it is necessary that the relative permeability and capillary pressure characteristics and displacing fluid viscosities be identical for the systems compared. A method is presented by which accurate prediction of the performance of a gravity drainage system can be made. Close agreement between experimental and calculated drainage performance shows that steady state relative permeability and static capillary pressure data can be used to describe fluid displacement behavior. The very wide range of liquid recoveries before gas breakthrough which result from production rate variation alone demonstrates the importance of this factor in planning depletion of a gravity drainage reservoir. Calculated results are presented which show that little additional recovery can be expected from a high pressure gravity drainage system between the times of gas breakthrough and attainment of such high gas/liquid ratios as to make further pressure maintenance impractical. INTRODUCTION It bas been recognized for some time that gravity forces play an important part in the recovery of oil from some types of reservoirs. Field experience has shown that under certain conditions, gravity drainage can result in very high oil recov- eries. Qualitative reasoning has led most engineers to the general conclusions that: (1) Where gravity drainage is important, the reservoir pressure should be maintained by gas injection at the crest of the structure to prevent shrinkage of the oil in place and to keep a low viscosity so the oil can drain at the fastest possible rate. (2) Recovery by gravity drainage is rate sensitive. A survey of the literature indicates that while considerable work has been done on the effects of gravity in oil production problems, no satisfactory method of calculating the performance of gravity drainage reservoirs has been reported.1,2,3,4,5 In the absence of any proven method of calculating reservoir performance, the level at which pressure should be maintained and the rate of production for most efficient operation has been open to debate. The limits and general character of the recovery us rate relationships are apparent from the following: As rates of production approach zero, oil will drain by virtue of its own weight fast enough that the saturation distribution is always approximately the same as the original static capillary distribution. That is, at the gas-oil contact the liquid saturation goes from 100 per cent to a value of perhaps 10 to 30 per cent in a very short distance, so that as production continues. a relatively sharp gas-oil contact moves down-structure, leaving a very low residual oil saturation in the originally oil saturated zone. At very high rates of production, the oil that drains by virtue of its own weight is negligible and recovery will approach that of a horizontal gas drive, leaving a much higher residual oil saturation. At intermediate rates of production. recovery should be between these extreme*. Despite the lack of quantitative methods, tentative agreement had been reached on the question of optimum producing rate. It has been generally concluded that rather sharp decrease in recovery would occur at rates of production above the "maximum rate of gravity drainage" and hence this rate should not be exceeded. The "maximum rate of gravity drainage" is defined as the rate of production from a 100 per cent

Jan 1, 1951

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By Robert L. Hoss

The application of Muskat's differential equations for predicting the performance of a solution gas-drive reservoir to the Fullerton field indicates that a recovery may be expected of 14,030 bbl per acre under primary depletion, and 15,910 and I7,190 bbl per acre, respectively, under the two pressure-maintenance programs analyzed. In this analysis it is assumed that fluid properties are the same throughout the reservoir, that the core data used are typical of the Clear Fork limestone in the locality, and that injected gas is completely diffused through the pay zones. The economic study of the three recovery programs considered indicates that the one in which 60 pct of the produced gas is returned to the formation is the most profitable. Introduction The purpose of this paper is to consider the effect of pressure maintenance by gas injection on the ultimate recovery of oil from a solution gas-drive reservoir. This prediction of reservoir behavior is the result of calculations based on the equations presented by M. Muskat in the Journal of Applied Physics, March 1945. Reservoir The reservoir analyzed is the Fullerton field in Andrews County, West Texas. It is situated on a broad anticlinal structure with approximately 330 ft of closure, and extends approximately eight miles north and south and five miles east and west. The areal extent is taken as 16,642 acres. The producing horizon is the Lower Clear Fork limestone of Permian age, occurring at depths from 6650 to 7350 ft. It is divided into four zones. For the purposes of this study, the fluid characteristics in all zones are assumed to be the same. Therefore an equivalent zone is assumed of thickness equal to the sum of the thicknesses of the four zones, or an average of 216 ft. The reservoir contains 3,586,991 acre feet of producing limestone, with an initial volume of stock tank oil in place of 918,-403,100 bbl. The average connate water saturation is 24.0 pct. The permeability ranges from 6.7 md. to 12.6 md., and the porosity ranges from 7.7 pct to 11.4 pct, except for occasional cavernous streaks. Based on these values, the limestone may be considered intergranu1ar.l Production histories and material-balance studies have indicated that the driving mechanism is solution gas. Water levels have been established. Theory It is possible to express mathematically the change in residual oil saturation with pressure decline in a solution gas-drive reservoir as being equal to the change in the physical properties of the fluids present in the reservoir as this pressure decline occurs. Essentially these equations, in differential form as presented by M. Muskat,2 express quantitatively the instantaneous net effect

Jan 1, 1948

By Robert L. Hoss

The application of Muskat's differential equations for predicting the performance of a solution gas-drive reservoir to the Fullerton field indicates that a recovery may be expected of 14,030 bbl per acre under primary depletion, and 15,910 and I7,190 bbl per acre, respectively, under the two pressure-maintenance programs analyzed. In this analysis it is assumed that fluid properties are the same throughout the reservoir, that the core data used are typical of the Clear Fork limestone in the locality, and that injected gas is completely diffused through the pay zones. The economic study of the three recovery programs considered indicates that the one in which 60 pct of the produced gas is returned to the formation is the most profitable. Introduction The purpose of this paper is to consider the effect of pressure maintenance by gas injection on the ultimate recovery of oil from a solution gas-drive reservoir. This prediction of reservoir behavior is the result of calculations based on the equations presented by M. Muskat in the Journal of Applied Physics, March 1945. Reservoir The reservoir analyzed is the Fullerton field in Andrews County, West Texas. It is situated on a broad anticlinal structure with approximately 330 ft of closure, and extends approximately eight miles north and south and five miles east and west. The areal extent is taken as 16,642 acres. The producing horizon is the Lower Clear Fork limestone of Permian age, occurring at depths from 6650 to 7350 ft. It is divided into four zones. For the purposes of this study, the fluid characteristics in all zones are assumed to be the same. Therefore an equivalent zone is assumed of thickness equal to the sum of the thicknesses of the four zones, or an average of 216 ft. The reservoir contains 3,586,991 acre feet of producing limestone, with an initial volume of stock tank oil in place of 918,-403,100 bbl. The average connate water saturation is 24.0 pct. The permeability ranges from 6.7 md. to 12.6 md., and the porosity ranges from 7.7 pct to 11.4 pct, except for occasional cavernous streaks. Based on these values, the limestone may be considered intergranu1ar.l Production histories and material-balance studies have indicated that the driving mechanism is solution gas. Water levels have been established. Theory It is possible to express mathematically the change in residual oil saturation with pressure decline in a solution gas-drive reservoir as being equal to the change in the physical properties of the fluids present in the reservoir as this pressure decline occurs. Essentially these equations, in differential form as presented by M. Muskat,2 express quantitatively the instantaneous net effect

Jan 1, 1948

By R. I. Parsons, M. B. Standing, E. N. Lindblad

The recovery of the heavier components from a gas cap or retrograde pool is shown to be the greatest when the sand is cycled with a dry gas at a low pressure. This conclusion is in direct opposition to the belief that the most efficient production program is pressure maintenance and cycling at or near the dew point. The results are calculated from: (I) constant volume, variable composition pressure-volume-temperature tests on a mixture of trap gas and liquid from a producing well; (2) published equilibrium constant data and the measured composition of the dew-point material; and (3) the fact that in a sand section of homogeneous permeability, injected gas displaces reservoir gas nearly, quantitatively. The results are based on the simplifying assumption that variable permeability systems may be defined by the ratio of two statistical parameters, and that gas injected into an actual sand will behave as though the sand were composed of many layers, each of constant permeability. Introduction Pressure decline in gas-condensate type reservoirs is accompanied by the formation of a liquid phase throughout the reservoir. Over the past ten years the processing of the material from these types of reservoirs for the heavier hydrocarbon components and the return of the light fractions to the reservoir ("cycling") has become increasingly popular. It has often been stated that the purpose of such a program is to prevent the loss of the retrograde liquid phase formed in the reservoir. The purpose of this paper is to present the results of laboratory tests and computations concerned with several possible methods of producing a gas cap or condensate type of reservoir. The results show that the recovery of heavier hydrocarbons for this type of reservoir is not a maximum under conditions of cycling at or near the dew-point pressure. Instead, variations in permeability and the ability of the dry injected gas to revaporize condensate, point to cycling at a considerably reduced pressure as the optimum production-method. In addition the paper suggests a way of evaluating sands for their permeability variation. The calculations presented in the following sections of the paper are in terms of the production of the butanes and heavier fraction from the gas cap of a field which had an original pressure of approximately 3000 psi. Specifically, the calculations show: 1. To what extent the condensate formed in the reservoir after pressure decline can be recovered by the cycling of a dry absorber-plant gas through the sands. 2. Whether there would be greater recovery of hydrocarbons if cycling were instigated at the original reservoir pressure.

Jan 1, 1948

By M. B. Standing, R. I. Parsons, E. N. Lindblad

The recovery of the heavier components from a gas cap or retrograde pool is shown to be the greatest when the sand is cycled with a dry gas at a low pressure. This conclusion is in direct opposition to the belief that the most efficient production program is pressure maintenance and cycling at or near the dew point. The results are calculated from: (I) constant volume, variable composition pressure-volume-temperature tests on a mixture of trap gas and liquid from a producing well; (2) published equilibrium constant data and the measured composition of the dew-point material; and (3) the fact that in a sand section of homogeneous permeability, injected gas displaces reservoir gas nearly, quantitatively. The results are based on the simplifying assumption that variable permeability systems may be defined by the ratio of two statistical parameters, and that gas injected into an actual sand will behave as though the sand were composed of many layers, each of constant permeability. Introduction Pressure decline in gas-condensate type reservoirs is accompanied by the formation of a liquid phase throughout the reservoir. Over the past ten years the processing of the material from these types of reservoirs for the heavier hydrocarbon components and the return of the light fractions to the reservoir ("cycling") has become increasingly popular. It has often been stated that the purpose of such a program is to prevent the loss of the retrograde liquid phase formed in the reservoir. The purpose of this paper is to present the results of laboratory tests and computations concerned with several possible methods of producing a gas cap or condensate type of reservoir. The results show that the recovery of heavier hydrocarbons for this type of reservoir is not a maximum under conditions of cycling at or near the dew-point pressure. Instead, variations in permeability and the ability of the dry injected gas to revaporize condensate, point to cycling at a considerably reduced pressure as the optimum production-method. In addition the paper suggests a way of evaluating sands for their permeability variation. The calculations presented in the following sections of the paper are in terms of the production of the butanes and heavier fraction from the gas cap of a field which had an original pressure of approximately 3000 psi. Specifically, the calculations show: 1. To what extent the condensate formed in the reservoir after pressure decline can be recovered by the cycling of a dry absorber-plant gas through the sands. 2. Whether there would be greater recovery of hydrocarbons if cycling were instigated at the original reservoir pressure.

Jan 1, 1948

By D. T. MacRoberts

A simple disturbance in a gas reservoir travels with a finite velocity which is nearly independent of the amplitude of the disturbance. As a result very complex transients may be set up which seriously affect observation at a well. The time of transient is much longer than commonly supposed. INTRODUCTION A change in pressure at a gas well is not propagated through the reservoir instantaneously but with a definite fin- ite velocity. Knowledge of the period ! of this transient is of value because: 1. Observations of pressure and flow-rate during the transient may affect interpretation of back-pressure tests. 2. Interference effects between wells may be used to determine continuity of formation. 3. Additional information concerning the physical characteris- tics may possibly be deduced from the behavior of a well during the transient. MATHEMATICAL DEVELOPMENT It will be assumed that the well is initially closed in and that the pressure is constant throlughout the reservoir. As the initial disturbance passes outward from the well it is further assumed that a logarithmic pressure distribution is established between the most distant point reached and the well, so that: <l> = ((/)r-<t>w) I ------------ + (/) . (1) Manuscript received at office of the Branch September 15. 1948. Paper presented at Branch Fall Meeting, Dallas. Texas. Oct. 4-6, 1948. r. = distance to any point between r and rw r = radius of drainge iw= radius of well (/) = potential (P for incompressible fluids, P2 for perfect gases) Up to the time the effect reaches the point r the production from the well will be: W= (Pr-P,,,) . (2) i b W = volume produced t = sand thickness F = fractional porosity Pr = reservoir pressure Pm = average pressure between r and rw Pb = atmospheric Pressure All pressures will subsequently be stated in atmospheres and the constant Pb will be omitted. This will be produced at a rate: dw t k (Pr³-PW2) do d 0 - . (3) log — Pw = pressure at the well k = permeability 0 = time = viscosity

Jan 1, 1949

By D. T. MacRoberts

A simple disturbance in a gas reservoir travels with a finite velocity which is nearly independent of the amplitude of the disturbance. As a result very complex transients may be set up which seriously affect observation at a well. The time of transient is much longer than commonly supposed. INTRODUCTION A change in pressure at a gas well is not propagated through the reservoir instantaneously but with a definite fin- ite velocity. Knowledge of the period ! of this transient is of value because: 1. Observations of pressure and flow-rate during the transient may affect interpretation of back-pressure tests. 2. Interference effects between wells may be used to determine continuity of formation. 3. Additional information concerning the physical characteris- tics may possibly be deduced from the behavior of a well during the transient. MATHEMATICAL DEVELOPMENT It will be assumed that the well is initially closed in and that the pressure is constant throlughout the reservoir. As the initial disturbance passes outward from the well it is further assumed that a logarithmic pressure distribution is established between the most distant point reached and the well, so that: <l> = ((/)r-<t>w) I ------------ + (/) . (1) Manuscript received at office of the Branch September 15. 1948. Paper presented at Branch Fall Meeting, Dallas. Texas. Oct. 4-6, 1948. r. = distance to any point between r and rw r = radius of drainge iw= radius of well (/) = potential (P for incompressible fluids, P2 for perfect gases) Up to the time the effect reaches the point r the production from the well will be: W= (Pr-P,,,) . (2) i b W = volume produced t = sand thickness F = fractional porosity Pr = reservoir pressure Pm = average pressure between r and rw Pb = atmospheric Pressure All pressures will subsequently be stated in atmospheres and the constant Pb will be omitted. This will be produced at a rate: dw t k (Pr³-PW2) do d 0 - . (3) log — Pw = pressure at the well k = permeability 0 = time = viscosity

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.&apos; 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,&apos; 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&apos; material balance equation and Hurst&apos;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.&apos; 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,&apos; 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&apos; material balance equation and Hurst&apos;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;&apos; 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;&apos; 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 N. De Nevers

A calculation method has been developed for carbonated water flooding. This method takes into account the effects of oil viscosity reduction and oil swelling, due to carbon dioxide transferred to the reservoir oil from the carbonated flood water. It shows the effect of changes in carbonation pressure and carbonated water slug size. The method is based on a Buckley-Leverett-type linear flow model. The mathematical approach is similar to that developed by Welge, et al, for enriched gas drives. However, it contains a very important advance over the previous papers in its solution to the problem of injection of a "slug" of carbonated water rather than continuous injection of fluid of constant composition. The method also shows the chromatographic nature of the transport of Co 2 in a carbonated water flood. Sample calculations indicate that viscosity reduction is the most important source of incremental oil recovery. Swelling contributes to a lesser extent. Raising the carbonation pressure increases recovery, but this effect diminishes as pressure is increased. increasing the slug size also increases recovery, but again this effect diminishes as slug size is increased. INTRODUCTION In recent years there has been considerable interest in carbonated water flooding as a secondary recovery method. Several carbonated waterflood field tests are now in progress.1&apos;3 Yet there is no published method for predicting the results of a carbonated water flood. This paper presents a method for predicting the oil recovery for ; carbonated water flood. In a carbonated water flood. a carbonated water slug is injected into the reservoir. Then this slug is driven forward by injecting plain water behind it. The field operator must choose the carbonated water slug size and the CO2:water ratio in the slug. The calculation method presented here shows the effects of both slug size and CO2:water ratio on oil recovery. The calculation assumes that the addition of carbon dioxide to flood water increases the oil recovery by two mechanisms: 1. Carbon dioxide transfers from the flood water to the oil, thereby reducing the viscosity of the oil. At the lower viscosity, the oil flows more readily and is more easily swept out by flood water. 2. The carbon dioxide which transfers to the oil phase causes the oil phase to swell. Assume that the residual oil saturation (volume per cent of pore space) is the same whether the oil is carbonated or not. Then the same volume will contain less stock tank oil if it is carbonated, due to oil swelling. The plain water driving the carbonated water slug forward will strip the CO2 out of this residual oil., causing this oil to shrink. This swelling-followed-by-shrinking leaves behind a lower residual oil saturation than would have been left behind by a plain water flood. Since some oils swell to almost twice their original volume at reasonable CO2 pressures, this effect may be large. Two additional benefits sometimes claimed for carbonated water flooding are: 1. Adding CO2 to flood water greatly increases the water injectivity in some cases.1 (Injectivity = BPD of water injected/ft of sand x differential bottom-hole pressure, psi.) CO2 forms a weak acid in water; injecting carbonated water is, therefore, somewhat like acid treating. If a formation responds to acid treating, it will probably have a higher injectivity for carbonated water than for plain water. Although the increased injectivity increases the profitability of a carbonated water flood over a plain water flood, it does not change the production curve based on total fluid injected. Therefore, the method presented here does not show any effect of increased injectivity. 2. It is claimed in one patent4 that carbon dioxide will react with some constituents of the oil to form detergents. These detergents are expected to increase the oil recovery. There is little data available to support this idea. The calculation method presented here takes no credit for extra oil recovery due to the formation of any such detergents.

Jan 1, 1965

By C. R. McEwen

The displacement equations of Buckley and Lever-ett&apos; have been successfully applied to the prediction of oil recovery in frontal drives for a number of years. Commonly, the capillary pressure term is omitted from Leverett&apos;s fractional flow formula. The resulting satura-tion-vs-distance curves become multi-valued, however. Buckley and Leverett, recognizing that the neglect of capillarity was most likely the cause, proposed that this part of the solution he ignored. An ingenious method of propagating a discontinuity in place of the multi-valued portlon of the curve was devised by Cardwell using "non-capillary displacement theory". However, it would seem that a completc description of the displacement pracess would require the inclusion of capillary forces. Three publications on displacement theory have appeared in which capillary pressure was included. Ter-williger. et al&apos;, calculated saturation profiles which matched those which they observed when gas displaced water vertically downward. Based on their study, Jones-Parra and Calhoun5 roposed a method of computing saturation distributions sufficiently removed in time from the initial conditions so that all saturations in the flood front region were assumed to move with the same velocity. This latter assumption corresponds to the tangent construction on the flowing fraction-vs-satura-tion curve proposed by Welge." A numerical solution to the problem of water displacing oil from a water-wet porous medium was recently presented by Blair, ef at&apos;. The initial conditions included a uniform saturation slightly above interstitial water saturation. that is, the permeability to water ahead of the front was not zero. This appeared to result in a region of slowly rising water saturation ahead of the front, at least for the particular set of relative permeability and capillary pressure data used to illustrate the computation. The authors also postulated an "end-effect" at the outflow face of the system, whereby no water could pass out of the system until the water saturation had built up to a high value. The saturation profiles computed under these assumptions show water accumulating at the downstream end of the system while the flood front proper is still traveling toward the outlet. The purpose of this paper is to present the results of a numerical solution to the one-dimensional water-flooding equation in which no outlet end-effect is considered and where there is no flow of water ahead of the flood front. It should be noted that the contemporary paper of Fayers and Sheldon8 deals with the problem under discussion. For a thorough treatment including both capillarity and gravity the reader is referred to their paper.

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&apos;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&apos;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 B. H. Caudle, M. D. Witte

In predicting the performance of a pattern injection operation, the engineer needs to know both the amount of oil to be recovered and the rate at which the recovcry will take place. This paper- describes fluid flow model studies on the effect of mobility ratio on the rate of oil recovery in a five-spot. The results show the change in fluid conductivity (total flow rate/pressure drop) as the sweep-ont pattern increases for mobility ratios between 0.1 (and 10. There data, when cornbitted with a knowledge of reservoir permeabilities and sweep-orit pattern efficiencies, cata he used to predict the cLrmulntivr oil production as a function of time for homogeneous five-spot injection systems. INTRODUCTION The flow potential in a pattern flood element is a function of the mobility ratio and how far the injected flood front has gone.1,2 For a waterflood project having a mobility ratio of one, the conductance—which is the flow rate per unit pressure drop (q/Ap)—-will remain the same throughout the life of the flood. If the mobility ratio is favorable ( < I), the conductance of the project will drop continuously throughout the operation. If, on the other hand, the mobility ratio is unfavorahlc. the conductance ratio will increase continuously. Therefore, if the reservoir engineer wishes to predict either the rate at a specified pressure or the pressure at some rate for a secondary recovery project, he must know how the conductance will vary during the flood. Other investigators have studied this problem using potentiometric models&apos; or analytical solution&apos;. The

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&apos;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&apos;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