Search Documents

By J. M. Bugbee

THE engineer and the conservationist agree that effective water drive is the desirable reservoir production mechanism. Water drive may result either from the expansion of edge water, the reservoir water bordering an oil pool, or from the artesian flow of edge water from the outcrop. Thus, the extent of the reservoir has an important bearing upon both the rate and the ultimate quantity of water influx to an oil pool. Lens-type reservoirs are generally recognizable, but the complete structure of the more important domal types, having possible limitations from crestal and synclinal faulting, is not well understood. The delineation from cores, electrical logs and exploitation phenomena of the crestal portions of reservoirs is increasingly exact. This paper suggests that inference from these factual pictures, combined with the scattered well and geophysical data of the synclines and with consideration of the dynamics of structural formation, may be applied to solving the complete geologic structure. Since faulting on the domes and in the synclines may effectively seal and thus limit the size of a reservoir and bar additional water influx, consideration of geologic structure and its possible effects is essential to complete analysis and prediction of reservoir performance.

Jan 1, 1942

By I. Haskett

ABSTRACT The Redwater pool, in north-central Alberta, was discovered in October, 1948, and at the end of 1950 was about 85 per cent developed. Production is from a Devonian coral limestone reef and is found at an average depth of 3,150 feet. The pool has an estimated .productive area of 37,000 acres, and a maximum oil zone thickness of slightly over 200 feet. Wide variations exist in both type and values of porosity, with the average value determined for the reservoir being 6.16 per cent. The total stock tank oil in place is about 1.24 billion barrels. The crude has a saturation pressure of 485 psig, com-pared to the original value of 1,050 psig for the reservoir, and the principal recovery mechanism so far has been the natural water drive. Cumulative production to January 1st, 1951, has been about 15.5 million barrels. These withdrawals resulted in a reservoir pressure decline of slightly over 100 psi. Predictions of future behaviour indicate that the water drive will not be of sufficient magnitude to maintain the reservoir above saturation pressure past 1953, and that, with natural depletion, the pool will be produced largely by a solution gas drive. Pressure maintenance by water injection has been considered and the possibility of raising the pressure to the original value of 1,050 psig is also discussed. INTRODUCTION WITH the discovery of the Lewduc pool, early in 1947, the Devonian limestone became the major geological formation contributing to Western Canada's rapidly growing petroleum reserves. Redwater, discovered in October, 1948, is now the largest known Devonian pool in Canada, with a recoverable reserve estimated at approximately 500 million barrels. Until recent years, the knowledge of the nature and characteristics of this type of limestone reservoir has been comparatively limited.

Jan 1, 1951

By W. F. Kieschnick, C. C. Miller, E. R. Brownscombe

The possibility has been mentioned that large pressure gradients in a solution gas driven field caused by high production rates might lead to a reduction in the ultimate recovery obtainable compared to that which would be obtained by a very slow rate of production. In the present study the reservoir conditions accompanying a high rate of production and the corresponding ultimate recovery were compared with those obtained if the reservoir were produced at some marginal rate throughout its entire life. The method of calculation involved the application of fluid flow - material balance analysis to a series of successive steady state conditions in the reservoir. In the case studied very little difference in recovery was obtained at the abandonment pressure. An analysis of the basic factors involved indicates that the sane results would hold for considerable variation in properties of reservoir fluids, permeability of the formation or well spacing. Even if the reservoir has a water drive it appears that no harm will be done by- open flow production if the rate is cut back before any appreciable free gas is produced. However, any condition which leads to disproportionate withdrawal rates and which causes large pressure differences over large areas might result in substantial loss in ultimate recovery. INTRODUCTION When the pressure declines in a reservoir, gas is evolved from solution causing a shrinkage in the reservoir oil so that a stock tank barrel of oil with its residual dissolved gas occupies a smaller volume than originally. For a given reservoir oil saturation, therefore. stock tank oil content will be higher at the lower pressure. At the same time the evolved gas which remains in the reservoir reduces the reservoir oil saturation, and this would be expected perhaps to cause a reduction in stock tank oil content. However, around the well the shrinkage occurs so rapidly that stock tank oil may actually accumulate in this region in spite of the simultaneous reduction in reservoir oil saturation. If this accumulation should continue, spreading into a larger zone, the loss might be appreciable and because larger pressure gradients accompany higher rates the accumulation might possibly be large if wells are produced at too high a rate. It is the purpose of this paper to determine the effect of rate of flow on ultimate recovery by solution gas drive. The problem under consideration involves the unsteady state flow of oil in a closed radial system from which oil is produced by the evolution of gas as a finite drawdown is estalllished and the reservoir pressure declines. A somewhat similar unsteady state production problem dealing, however, with a linear system has been considered by Muskat and Meres, but these authors do not discuss the effects which rate may have on recovery. Other methods have been reported for estimaling recovery as a function of pressure decline when producing at an infinitesimal rate, and Muskat and Taylor' have shown the effect of fluid and reservoir characteristics on recovery under these conditions. Assuming steady state to exist, Babson has given a procedure for computing the decline in well productivity as the field is depleted, and Moyer has combined these methods to estimate abandonment conditions and the resulting ultimate recovery. These latter approaches to the problem either do not consider the effect of rate or assume steady state can be established. In detail this paper presents the behavior of the reservoir producing under unsteady state conditions and considers the effect of production rate during this unsteady state flow on the recovery attained. The analysis also shows that there will exist at abandonment a type of steady state, and on this basis a method has been clevelopetl for rapid calculation of abandonment conditions.

Jan 1, 1949

By C. C. Miller, W. F. Kieschnick, E. R. Brownscombe

The possibility has been mentioned that large pressure gradients in a solution gas driven field caused by high production rates might lead to a reduction in the ultimate recovery obtainable compared to that which would be obtained by a very slow rate of production. In the present study the reservoir conditions accompanying a high rate of production and the corresponding ultimate recovery were compared with those obtained if the reservoir were produced at some marginal rate throughout its entire life. The method of calculation involved the application of fluid flow - material balance analysis to a series of successive steady state conditions in the reservoir. In the case studied very little difference in recovery was obtained at the abandonment pressure. An analysis of the basic factors involved indicates that the sane results would hold for considerable variation in properties of reservoir fluids, permeability of the formation or well spacing. Even if the reservoir has a water drive it appears that no harm will be done by- open flow production if the rate is cut back before any appreciable free gas is produced. However, any condition which leads to disproportionate withdrawal rates and which causes large pressure differences over large areas might result in substantial loss in ultimate recovery. INTRODUCTION When the pressure declines in a reservoir, gas is evolved from solution causing a shrinkage in the reservoir oil so that a stock tank barrel of oil with its residual dissolved gas occupies a smaller volume than originally. For a given reservoir oil saturation, therefore. stock tank oil content will be higher at the lower pressure. At the same time the evolved gas which remains in the reservoir reduces the reservoir oil saturation, and this would be expected perhaps to cause a reduction in stock tank oil content. However, around the well the shrinkage occurs so rapidly that stock tank oil may actually accumulate in this region in spite of the simultaneous reduction in reservoir oil saturation. If this accumulation should continue, spreading into a larger zone, the loss might be appreciable and because larger pressure gradients accompany higher rates the accumulation might possibly be large if wells are produced at too high a rate. It is the purpose of this paper to determine the effect of rate of flow on ultimate recovery by solution gas drive. The problem under consideration involves the unsteady state flow of oil in a closed radial system from which oil is produced by the evolution of gas as a finite drawdown is estalllished and the reservoir pressure declines. A somewhat similar unsteady state production problem dealing, however, with a linear system has been considered by Muskat and Meres, but these authors do not discuss the effects which rate may have on recovery. Other methods have been reported for estimaling recovery as a function of pressure decline when producing at an infinitesimal rate, and Muskat and Taylor' have shown the effect of fluid and reservoir characteristics on recovery under these conditions. Assuming steady state to exist, Babson has given a procedure for computing the decline in well productivity as the field is depleted, and Moyer has combined these methods to estimate abandonment conditions and the resulting ultimate recovery. These latter approaches to the problem either do not consider the effect of rate or assume steady state can be established. In detail this paper presents the behavior of the reservoir producing under unsteady state conditions and considers the effect of production rate during this unsteady state flow on the recovery attained. The analysis also shows that there will exist at abandonment a type of steady state, and on this basis a method has been clevelopetl for rapid calculation of abandonment conditions.

Jan 1, 1949

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 H. E. Stamm, H. H. Spain, R. C. Rumble

This paper presents a reservoir analyzer study of the performance of the Woodbine formation in the East Texas basin. The study was made possible by the compilation of available information on the configuration. thickness and the pressure drawdown of the formation. The investigation was made of the pressure variations in the basin incident to production from the several Woodbine reservoirs. In addition, the apparent compressibility of Woodbine water was evaluated so that the potential water yield of the formation could be determined. The distribution of permeability of the Woodbine sand and the interference between producing areas were also investigated. In conjunction with this basic study of the Woodbine sand, an acceptable match of the production-pressure relationship in the East Texas Field was established on the analyzer. INTRODUCTION When production is taken from a reservoir contiguous to an aquifer, the resultant pressure gradient causes a water influx into the reservoir. For some years it has been the practice to predict the performance of water drive reservoirs with mathematical equations that relate the water influx into reservoirs to their pressure behavior. A rigorous solution of these equations is quite complex: hence, the conventional method requires the simplifying assumptions that the formation have constant thickness, permeability and poiosity, and that it have a known and regular shape. It is necessary to use performance records of the reservoir in question to determine the equation constants that make the results of the mathematical calculations duplicate the reservoir history. The values of the constants determined in this manner sometimes differ so greatly from theoretical or measured values that the constants appear unreasonable. The equation constants comprise such factors as the thickness, porosity and permeability of the formation, and the viscosity and compressibility of water. Usually the thickness, porosity and viscosity are known within a reasonable degree of certainty. whereas the formation permeability and the water compressibility are known with considerably less certainty. Hence. permeability and compressibility values are arbitrarily selected that will make the mathematical equations duplicate the reservoir behavior. In many reservoir studies the value of the compressibility of water thus selected has exceeded the reported value by a factor of from severalfold to manyfold. Whether this excessive value results from the inadequacy of the mathematical approach to take into account the shape irregularities and the heterogeneities of the aquifer, or is due to some type of formation compaction and decrease of porosity has not been established. Because of the significance of the potential water yield of formations with reference not only to oil reservoirs but to water resources as well, it is very important that a representative value for the effective water compressibility be determined. One of the best known examples of a field for which it has been possible in the past to predict the future reservoir pressure is the East Texas Field. Until about 1943 it was possible by means of mathematical equations to duplicate past pressure behavior and to predict with a high degree of accuracy the reservoir pressure of the field. Since 1943, however, the calculated pressures have deviated progressively from the observed pressures to such a degree that in recent years, it has

Jan 1, 1951

By H. H. Spain, R. C. Rumble, H. E. Stamm

This paper presents a reservoir analyzer study of the performance of the Woodbine formation in the East Texas basin. The study was made possible by the compilation of available information on the configuration. thickness and the pressure drawdown of the formation. The investigation was made of the pressure variations in the basin incident to production from the several Woodbine reservoirs. In addition, the apparent compressibility of Woodbine water was evaluated so that the potential water yield of the formation could be determined. The distribution of permeability of the Woodbine sand and the interference between producing areas were also investigated. In conjunction with this basic study of the Woodbine sand, an acceptable match of the production-pressure relationship in the East Texas Field was established on the analyzer. INTRODUCTION When production is taken from a reservoir contiguous to an aquifer, the resultant pressure gradient causes a water influx into the reservoir. For some years it has been the practice to predict the performance of water drive reservoirs with mathematical equations that relate the water influx into reservoirs to their pressure behavior. A rigorous solution of these equations is quite complex: hence, the conventional method requires the simplifying assumptions that the formation have constant thickness, permeability and poiosity, and that it have a known and regular shape. It is necessary to use performance records of the reservoir in question to determine the equation constants that make the results of the mathematical calculations duplicate the reservoir history. The values of the constants determined in this manner sometimes differ so greatly from theoretical or measured values that the constants appear unreasonable. The equation constants comprise such factors as the thickness, porosity and permeability of the formation, and the viscosity and compressibility of water. Usually the thickness, porosity and viscosity are known within a reasonable degree of certainty. whereas the formation permeability and the water compressibility are known with considerably less certainty. Hence. permeability and compressibility values are arbitrarily selected that will make the mathematical equations duplicate the reservoir behavior. In many reservoir studies the value of the compressibility of water thus selected has exceeded the reported value by a factor of from severalfold to manyfold. Whether this excessive value results from the inadequacy of the mathematical approach to take into account the shape irregularities and the heterogeneities of the aquifer, or is due to some type of formation compaction and decrease of porosity has not been established. Because of the significance of the potential water yield of formations with reference not only to oil reservoirs but to water resources as well, it is very important that a representative value for the effective water compressibility be determined. One of the best known examples of a field for which it has been possible in the past to predict the future reservoir pressure is the East Texas Field. Until about 1943 it was possible by means of mathematical equations to duplicate past pressure behavior and to predict with a high degree of accuracy the reservoir pressure of the field. Since 1943, however, the calculated pressures have deviated progressively from the observed pressures to such a degree that in recent years, it has

Jan 1, 1951

By Gray L. L., Max Littlefield, Godbold A. C.

The West Edmond pool of Central Oklahoma, a limestone reservoir, has an area in excess of 29,000 acres and as of Sept. 15, 1946, had produced 53 million barrels of oil from 731 wells at an average depth of 6900 ft. Water has encroached into the reservoir along the west side of the pool and although the area of water invasion is large the net volume of water influx is small, because encroachment has been primarily into a fracture system. The most significant part of this study concerns the character of the Hunton reservoir. Geological study of cores of the producing section was made to supplement core-analysis data. It was determined that approximately 90 pct of the pore volume was contained in intergranular or sand-like porosity and 10 pct in fractures (intermediate porosity). The porous limestone is divided into blocks by the fractures, although in the strongly dolomitic parts of the producing section, fractures are less well developed. The low-permeability intergranular porosity, largely subsidiary to the high-permeability fractures, will produce oil into the fractures by evolution of solution gas. The fractures in effect serve as drainage channels for production of oil from the inter-granular porosity. Because of such a large difference in permeabilities of the two com-ponents, a high degree of by-passing will occur, and, accordingly, the economics involved in undertaking full-scale high pressure gas injection operations have been seriously questioned. Introduction The West Edmond Hunton pool, which covers portions of Canadian, Logan, Oklahoma, and Kingfisher Counties in Central Oklahoma, is the largest reservoir producing oil from a limestone in the state. The total area of the pool comprises about 29,240 acres and development is considered to be essentially complete. In all, 731 producing wells have been completed in the Hunton reservoir at an average depth of approximately 6900 ft, and have yielded a cumulative oil production in excess of 53 million barrels as of Sept. 15, 1946. In addition to wells completed in the Hunton limestone, 21 wells have been completed in the Bartlesville sand at a depth of approximately 6500 ft, which have produced a total of 730,000 bbl of oil. One small producer has been completed in the Cleveland sand at a depth of 5700 ft. However, a study of only the Hunton limestone reservoir will be presented in this paper. The most significant phase of this study concerns the special type of geological core examination that was made in an effort to ascertain the type and degree of porosity in the rock in accordance with lithologic types. Results of this geological study showed the west Edmond pool to be of such a complex nature as to appreciably affect the behavior normally associated with homogeneous reservoirs. Recognition of this condition in the

Jan 1, 1948

By Max Littlefield, Godbold A. C., Gray L. L.

The West Edmond pool of Central Oklahoma, a limestone reservoir, has an area in excess of 29,000 acres and as of Sept. 15, 1946, had produced 53 million barrels of oil from 731 wells at an average depth of 6900 ft. Water has encroached into the reservoir along the west side of the pool and although the area of water invasion is large the net volume of water influx is small, because encroachment has been primarily into a fracture system. The most significant part of this study concerns the character of the Hunton reservoir. Geological study of cores of the producing section was made to supplement core-analysis data. It was determined that approximately 90 pct of the pore volume was contained in intergranular or sand-like porosity and 10 pct in fractures (intermediate porosity). The porous limestone is divided into blocks by the fractures, although in the strongly dolomitic parts of the producing section, fractures are less well developed. The low-permeability intergranular porosity, largely subsidiary to the high-permeability fractures, will produce oil into the fractures by evolution of solution gas. The fractures in effect serve as drainage channels for production of oil from the inter-granular porosity. Because of such a large difference in permeabilities of the two com-ponents, a high degree of by-passing will occur, and, accordingly, the economics involved in undertaking full-scale high pressure gas injection operations have been seriously questioned. Introduction The West Edmond Hunton pool, which covers portions of Canadian, Logan, Oklahoma, and Kingfisher Counties in Central Oklahoma, is the largest reservoir producing oil from a limestone in the state. The total area of the pool comprises about 29,240 acres and development is considered to be essentially complete. In all, 731 producing wells have been completed in the Hunton reservoir at an average depth of approximately 6900 ft, and have yielded a cumulative oil production in excess of 53 million barrels as of Sept. 15, 1946. In addition to wells completed in the Hunton limestone, 21 wells have been completed in the Bartlesville sand at a depth of approximately 6500 ft, which have produced a total of 730,000 bbl of oil. One small producer has been completed in the Cleveland sand at a depth of 5700 ft. However, a study of only the Hunton limestone reservoir will be presented in this paper. The most significant phase of this study concerns the special type of geological core examination that was made in an effort to ascertain the type and degree of porosity in the rock in accordance with lithologic types. Results of this geological study showed the west Edmond pool to be of such a complex nature as to appreciably affect the behavior normally associated with homogeneous reservoirs. Recognition of this condition in the

Jan 1, 1948

By W. R. Davis

This paper presents a sunimary of the results of a study of the Dornick Hills-Springer sand reservoir of the Velma Pool. Stephens County, Oklahoma. The data are of interest because this type of reservoir is unusual in the mid-continent area. The reservoir is characterized by very steep dips and a relatively thick producing section. Its up-dip productive limit is defined by truncation of the productive sands, and the down-dip limit is defined by deterioration of sand quality rather than by the presence of an oil-water contact. The crude oil is of 29" API gravity, is relatively viscous and was saturated with dissolved gas at the original reservoir pressure. Consideratiqn of the physical characteristics of the reservoir indicates that natural water encroachment will be negligible, and that the producing mechanism will be a dissolved gas drive augmented by gravity drainage. Performance to date tends to confirm this. INTRODUCTION It is the purpose of this paper to present some of the features of the Dor-nick Hills - Springer sand reservoir which are unusual in mid-continent operations. The study was undertaken in an attempt to define the producing mechanism or mechanisms that will control reservoir performance, as an aid in determining well completion and production practices. GEOLOGY The reservoir under consideration is a series of sand bodies of Pennsylva-nian age found in the Lower Dornick Hills and/or Springer formations. For the purpose of this paper, the sand series is referred to as the Dornick Hills-Springer sand. This sand series lies on the northeast flank of the Velma anticline, and dips very sharply to the northeast in TWP 1S-5W, the average dip being about 45º. However, as shown in Fig. 1. the trend of the pool swings eastward in TWP 1S-4W with a coincident de- crease in dip angle. In this portion of the pool the average dip is about 25". Faults in Sections 30-1S-4W and 32-1S-4W are shown in Fig. 1. These faults show only minor displacements for the horizon contoured in Fig. 1, but they do show a pronounced displacement of the immediately overlying beds. Similar faults are thought to exist in TWP 1S-5W; however, since the detailed geological studies necessary to properly define these faults were not available to the author, no attempt to show them in Fig. 1 has been made. The up-dip productive limit throughout the pool is determined by the truncation of the Lower Dornick Hills-Springer sand. The down-dip productive limit does not appear to be defined by edge water, but by a gradual reduction in permeability of the sands. No dry hole has produced a significant volume of water from this sand section. There was a small gas cap in TWP 1S-5W. Most of the wells in TWP 1S-SW encountered an oil sand 70 ft to 150 ft thick, a shale 20 ft to 80 ft thick, and another oil sand 50 ft to 130 ft thick. A few scattered wells in this area encountered other and bodies below these

Jan 1, 1951

By W. R. Davis

This paper presents a sunimary of the results of a study of the Dornick Hills-Springer sand reservoir of the Velma Pool. Stephens County, Oklahoma. The data are of interest because this type of reservoir is unusual in the mid-continent area. The reservoir is characterized by very steep dips and a relatively thick producing section. Its up-dip productive limit is defined by truncation of the productive sands, and the down-dip limit is defined by deterioration of sand quality rather than by the presence of an oil-water contact. The crude oil is of 29" API gravity, is relatively viscous and was saturated with dissolved gas at the original reservoir pressure. Consideratiqn of the physical characteristics of the reservoir indicates that natural water encroachment will be negligible, and that the producing mechanism will be a dissolved gas drive augmented by gravity drainage. Performance to date tends to confirm this. INTRODUCTION It is the purpose of this paper to present some of the features of the Dor-nick Hills - Springer sand reservoir which are unusual in mid-continent operations. The study was undertaken in an attempt to define the producing mechanism or mechanisms that will control reservoir performance, as an aid in determining well completion and production practices. GEOLOGY The reservoir under consideration is a series of sand bodies of Pennsylva-nian age found in the Lower Dornick Hills and/or Springer formations. For the purpose of this paper, the sand series is referred to as the Dornick Hills-Springer sand. This sand series lies on the northeast flank of the Velma anticline, and dips very sharply to the northeast in TWP 1S-5W, the average dip being about 45º. However, as shown in Fig. 1. the trend of the pool swings eastward in TWP 1S-4W with a coincident de- crease in dip angle. In this portion of the pool the average dip is about 25". Faults in Sections 30-1S-4W and 32-1S-4W are shown in Fig. 1. These faults show only minor displacements for the horizon contoured in Fig. 1, but they do show a pronounced displacement of the immediately overlying beds. Similar faults are thought to exist in TWP 1S-5W; however, since the detailed geological studies necessary to properly define these faults were not available to the author, no attempt to show them in Fig. 1 has been made. The up-dip productive limit throughout the pool is determined by the truncation of the Lower Dornick Hills-Springer sand. The down-dip productive limit does not appear to be defined by edge water, but by a gradual reduction in permeability of the sands. No dry hole has produced a significant volume of water from this sand section. There was a small gas cap in TWP 1S-5W. Most of the wells in TWP 1S-SW encountered an oil sand 70 ft to 150 ft thick, a shale 20 ft to 80 ft thick, and another oil sand 50 ft to 130 ft thick. A few scattered wells in this area encountered other and bodies below these

Jan 1, 1951

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