Search Documents

By W. E. Howlett, R. L. Slobod

Jan 1, 1965

By B. Habermann

Artificially consolidated sand models, representing one-quarter of a five-spot, have been developed and used to study factors aflecting misciblt. displacrmenr. Sweep efficiency at breakthrough, size of the mixing zone between two miscible liquids and per cent area contacted by drive after breakthrough were determined tor the high mobility ratios encountered in actual resvrvoirs. Quantitative relationships between the degree of viscous fingering and mobility ratio were obtained by measuring the length of the fluid interface. Scaled miscible-slug experiments, supported by field evidence reported in the literature, have shown that when the slug is followed by dry gas the process is less efficient than expected. In addition to low areal sweep efficiencies encountered for high mobility ratio displacements, the effectiveness of a miscible slug is greatly reduced. This is caused by an accelerated growth of the mixing zone between the driving and the displacing fluids. INTRODUCTION Considerable interest in miscible displacement as a secondary recovery process has been shown by the oil industry in recent years. In addition to extensive laboratory investigations, numerous field operations have been initiated so that 39 miscible displacement projects had been reported' in the United States by early 1959. By far the most popular approach in these projects is the use of an LPG slug followed by dry gas. Laboratory tests have indicated that with a small miscible slug, amounting to several per cent of the hydrocarbon pore volume, practically 100 per cent recovery of the oil in place can be achieved and that all of the slug material can be recovered by a subsequellt gas drive. However, these experiments were conducted on small-diameter long cores that can be considered as essentially uni-dimensional reservoirs that do not represent actual field conditions. It soon became apparent that the inherently adverse mobility ratios existing in the reservoir cannot be ignored in any miscible displacement field project. When a less viscous fiuid such as an LPG slug drives a more viscous reservoir oil, fingers develop that reduce the areal sweep efficiency. The slze of the miscible slug to be used is a subject of some controversy. It has been proposed that the size of siug needed be calculated by setting lengths of the mixing zone proportional to the square root of the distance traveled. This proposal is based on experimental evidence' and theoretical considerations"' in which fingering was not a factor. A prerequisite for successful application of the miscible displacement process to field operations is a better understanding of the following factors. L What is the effect of mobility ratio on the rate of deterioration of a miscible slug? 2. How large a slug should be used for given reservoir conditions? 3. What is the sweep efficiency? The purpose of this investigation was to provide some additional laboratory data that should help to answer these questions. More specifically, it was intended to study more quantitatively the effect of viscous fiinger-ing, mobility ratio and slug size on the efficiency of this relatively new secondary recovery process by means of experiments on thin, horizontal, porous medium models.

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By H. Groeneveld, C. A. Connally, P. J. Hoenmans, J. J. Justen, W. L. Mason

A miscible displacement pilot using a slug of LPG driven by separator gas was conducted in the Cardiurn reservoir of the Pembina field. The injection pattern was a 10-acre, inverted, isolated five-spot. Upon completion of the LPG-gar phase, an experiment was conducted using a slug of water followed by gas. Calculated performance of the pilot is compared with actual performance. Equations are developed to calculate the distribution of LPG into zones of varying permeability, to estimate the progress of the flood at different times in the various zones and to estimate gas rates after breakthrough. The analysis indicates that permeability stratification was a dominant factor in controlling oil recovery and that oil was completely displaced from the swept pore volume. The results of the pilot indicated that miscible flooding is a practical means of pressure maintenance in this reservoir. The total recovery from the pilot area was good in spite of the early breakthrough of LPG. The effects of stratification were reduced by injecting a slug of water into the partially swept reservoir. INTRODUCTION The Pembina field,' located in Alberta, is the largest oil field in Canada and one of the largest in the North American continent. The reservoir is a stratigraphic trap producing from the Cardium sand. Neither bottom water nor free gas has been found. The recovery of oil by the natural depletion mechanism has been estimated at 12.5 per cent. Pressure maintenance studies of various areas have indicated that the recovery can be increased 21/2 times by water flooding, and a large area of the field is presently under water flood. However, reservoir studies of the North Pembina area indicated that miscible flooding might be competitive with water flooding. A pilot test was conducted to evaluate the performance of a miscible flood. A 10-acre, inverted, isolated, five-spot pattern was selected for the pilot. The pattern area was large enough to minimize wellbore fracturing effects and contained sufficient oil to provide significant working numbers. The performance of each of the four producers could be evaluated individually and compared. In the event of breakthrough in one direction, the effect would be isolated from the other producers. The use of a single injector minimized the volume of LPG required, and, because of the high mobility of gas, one well was sufficient to inject the necessary daily volume to replace the high rate of production. With four producers, the test could be completed in time for results to be evaluated, additional engineering studies to be made and a unit to be formed before the reservoir pressure in the North Pembina area declined below the bubble point. The pilot was located in an area developed on staggered, 80-acre spacing. The injection well was drilled at a regular location, while the four producers were drilled 467-ft north, east, south and west of the injector. Each quadrant and its associated producer were identified according to their direction from the injector— that is, north, east, south or west. The eight surrounding producers on 80-acre spacing were shut in to isolate the pilot area and provide for reservoir pressure observation. The pilot wells were completed using permanent-type completion techniques. After coring, casing was run through the pay section and cemented. Inside 51/2-in. casing, 2 1/2-in. tubing was hung. The wells were perforated opposite the Upper Cardium sand and lightly fractured. Fracturing volumes, rates and pressures were low to minimize the extent of the fractures. The fracturing treatments average 1,000 lb of 20-40 mesh sand in 700 gal of a low fluid-loss sand-carrying agent. Feed rates and wellhead fracturing pressures averaged 5.5 bbl/min at 2,535 psig, respectively. After fracturing, the productivity index was measured in each of the five pilot wells. The average PI of the four producers was 0.41 BOPD/psig drawdown. The measured PI'S were approximately the same as PI'S calculated from core analysis data, indicating that the fracturing treatments were just sufficient to overcome

By M. M. Mebta, G. W. Dean, W. H. Somerton

With the advent of underground heating operations, interest has developed in the alteration of rock properties by high-temperature treatment. In the present work a number of sandstones were heated to temperatures in the range of 400°C to 800°C under both atmospheric and simulated reservoir pressures. Pertneabilities increased by at least 50 per cent and sonic velocities and breaking strerlgths decreased by an equivalent amount. Differential thermal expansion and other reactions of constituent min-era1 grains are the causes of these alterations. INTRODUCTION In the underground combustion of petroleum reservoirs, temperatures of the order of 600C are reported to have been reached in the combustion zone.' At this temperature rocks are subject to extensive thermal alteration. Temperatures of this magnitude and higher may also occur in subsurface formations when subjected to bottom-hole heating, thermal drilling operations, and underground nuclear explosions. Temperatures of this magnitude might also be generated by conventional rock drilling methods at points of bit-tooth contact. In earlier work, the permanent deformation of rocks resulting from heating was reported. Major structural damage of rocks occurs due to differential thermal expansion of mineral constituents. A number of mineral alterations including crystal inversions, loss of water of crystallization and dissociation, may also contribute to changes in physical structure and properties of rock. In the present work, samples of three typical sandstones were heated to several temperatures up to a maximum of 800C and then allowed to cool to room temperature. Heating was done under both atmospheric pressure and simulated reservoir pressure conditions. Physical properties of the samples were measured before and after heating and comparisons made. Measured properties included permeability, sonic velocity, breaking strength and fracture index. Changes in physical properties were compared to changes in mineralogical characteristics as determined by thin-section. X-ray diffraction and chemical tests. EQUIPMENT AND PROCEDURE Two outcrop sandstones (Bandera and Berea) and one sub-surface sandstone (St. Peter) were selected for the tests. These samples have a wide ranee in composition and physical properties as shown in Table 1. The first series of tests was made on 2-in. diameter by 5-in. long test specimens. Test specimens used in all later work were 3/4-in. diameter by 1 1/8-in. long, this being the specimen size required for heating at simulated reservoir pressures. After careful washing, the cores were oven dried at 100 ± 5C for a minimum of 24 hours before the tests were run. Test specimens were heated in an electric furnace at a constant rate of temperature increase of 3C per minute. When maximum temperature of the run was reached, the sample was allowed to soak for one hour. The furnace was then cooled to room temperature at the same rate of 3C per minute. The entire heating operation was designed for reproducibility without subjecting the test specimens to excessive thermal shock. For samples heated under simulated reservoir pressures, a pressure cell designed by Dean was used (Fig. 1).3 The core sample was inserted into a thin-walled (0.006 to 0.01-in.) copper cup which was then mounted in a high-pressure cell. Provisions were made for the application of internal pore pressure as well as confining pressure. Tests showed that the thin-walled copper cup closed tightly around the core and satisfactorily transmitted confining pressure to the core. The core was heated by placing the entire cell into the electric furnace. The heating program was the same as that used in the atmospheric pressure runs: tempera-ture rise of 3C per minute to maximum temperature of the test, soaking at maximum temperature for one hour, and cooling at a rate of approximately 3C per minute. The cell was designed to withstand 5,000 psi at 1,000C. However, since it was considered likely that repeated heating and cooling would in time weaken the steel, 2,000 psi at 850C was set as a working limit. In the present series of tests, the pore pressure was held constant at 750 psi and the confining pressure at 1.500 psi. The pressure source was a high-pressure nitrogen tank. The two pressures were controlled manually but are accurate well within ± 50 psi.

Jan 1, 1966

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By G. H. Alani, H. T. Kennedy

One of the major difficulties in predicting the performance of oil reservoirs from their early pressure history lies in the uncertainty of estimating the volume of the liquid hydrocarbons contained in them. As a first step in filling this need, an equation was developed to determine the Inolal liquid volume of pure hydrocarbons over a wide range of temperature and pressure. The second step consisted of adapting the equation to apply to mixtures, with the heary hydrocarbons expressed as c The equatious are similar in form to van der Waals equation. hut the constants a and b tire considered as furnctions of temperature. In addition to the gas constant R, there are four constants characteristic of each hydro- When comptired with experimental values jound in the literal the. the average absolute deviation in the calculated molal volumes is found to be a Maximum of 0.33 per -cem for any of the pure liquid hydrocarbons studied. This maximum deviation was that found when comparing the calculated and observed values over a temperature range of 86° to 482°F and a pressure range from the bubbl-point ro 10,000 psis. The equations expressing the correlation for. mixrurrc were developed from 647 experimental measure-rilents of volume on 47 bouom-hole samples covering ; temperature range of 72" to 250°F and a pressure znge from bubble-point ro 5,000 psig. The average absotute deviation was found to he 1.6 per cent with the maximum for any measurement of 4.9 per cent. INTRODUCTION Accuratc information of the pressure-volume-temper-iiture behavior of hydrocarbon liquids is of considerable importance in the field of both applied and theoretical science and, especially, in the solution of petroleum reservoir engineering problems. These PVT relationships can be expressed graphically, in tabular form or as equations of state. Equations of state for gases have been developed based on the ratio of pressure, temperature and volume to the pressure, temperature and volume at the critical points of the hydrocarbons. These equations have been successful in yielding values fairly accurate for reservoir engineering purposes. For liquids, however, relations of this type are more complex and have not yielded useful results. The composition of reservoir liquids is normally expressed in terms of the moles of hydrocarbons C,? C,. C,. elc, per mole of mixture. The composition must include the fraction of a mixture such as C,,, which cannot be separated with reasonable time and effort into specific hydrocarbons. Although there are several equations of state having variable degrees of accuracy applicable to pure liquid hydrocarbons, none includes a means for estimating the volume of mixtures as a function of temperature. pressure and composition expressed in this manner.

By S. J. Prison, C. R. Smith

The effect of fluid injection to control water coning in oil and gas wells was investigated. Analytical and model techniques were employed. The factors investigated were the position and length of the completion interval, the point of fluid injection, the viscosity of the injected fluid and the relative thickness of the oil and water sections. The resulting influence of these factors on the net producing water-oil ratio was determined. Several important conclusions can be drawn from the study. In general, it was found that the net producing water-oil ratio can be reduced by fluid injection. The magnitude of this reduction depended on the factors listed above. An important practical consideration is that the injection fluid may be either oil or water. If the injected fluid is less dense than the connate water of the reservoir, the fluid will not be lost. This fact is reassuring when valuable oil is being injected. Efforts to suppress water production were more successful when the injection fluid was more viscous than the reservoir oil, or when a zone of reduced permeability existed in the vicinity of the point of fluid injection. Under test conditions, little benefit was derived through the use of impermeable barriers or cement "pancakes". INTRODUCTION The occurrence of water coning has been known for at least 60 years. In thin oil or gas pay sections, the presence of an oil-water or gas-water contact hinders production and often causes early abandonment of the afflicted well if a completion is even attempted. Even when relatively thick pay sections are found, the encroachment of water when a water drive is present will eventually pose serious water coning problems. This water is often corrosive, expensive to separate from the oil or gas and is costly to dispose of. The theory of water coning has been discussed by a number of authors. 1,2,3,4 Briefly, water coning to the producing interval in a well is due to pressure gradients resulting from the production of fluid from the reservoir. These pressure gradients will cause a water cone to rise toward the bottom of the producing interval if a water-oil or water-gas contact exists. The tendency of the water to cone is offset or partially offset by gravity forces since the water has a higher specific gravity than the oil. A balance then exists between two forces, gravitational forces arising from the difference in specific gravities of the oil and water, and the pressure gradients causing the flow of fluids to the wellbore. If the pressure gradient exceeds the gravitational force, water coning to the wellbore occurs and water production results. Through the years considerable thought has been given to the water coning problem. More than 50 U. S. Patents have been granted to inventors on the subject. A relatively complete literature review of the water coning problem has been made. 5 A number of these patents hold considerable promise for the solution or the partial solution of the water and/or the coning problem. Very little has been written describing field tests of techniques for the suppression of water coning. A notable exception is the paperby West. 6 He reports success in reducing gas coning by a combination of gravel packing and oil injection above the oil-producing interval. He also describes a comparable method to prevent water coning, but provides no field examples. This study experimentally and analytically verifies the benefits of oil injection as a means of partially or completely suppressing the water cone. While the gas coning problem was not treated, it is anticipated that results comparable to those obtained in water suppression could be obtained with reduced oil injection since the viscosity contrast between oil and gas exceeds that between oil and water. For the purpose of this study, the conventional potential flow theory was applied to the water-coning problem. The experimental verification centered on both a radial and a linear model. The model study permitted the investigation of complex flow configurations and the use of fluids of differing densities and viscosities. No analytic expressions are available to permit a solution of the problem as stated (see Fig. 1).

By D. C. Lindley, D. H. Gaucher

The displacement of oil by water in a waterflood project is accomplished by the action of transient viscous, gravitational and capillary forces which drive fluid through interconnecting pore spaces toward production wells. The relative importance of each of these forces to the effectiveness of the displacement and the production history depends on properties of the reservoir itself and the manner in which it is operated. There has been considerable discussion'-3 concerning the importance of these factors and the extent to which they control the performance of individual waterflood projects. No data have been presented, however, to demonstrate the effect on waterflood behavior of variations in rock permeability, water-injection rates and mobility ratios in three-dimensional reservoir systems in which viscous, gravitational and capillary forces are allowed to assume the proper relative influence on fluid flow. Much information relevant to water flooding has been gained from the study of field case histories. But, the complexity of the systems involved, the difficulty of defining geometry and obtaining sufficient data, and the fact that the same reservoir is never flooded twice with the same initial conditions, all obscure the role of process variables such as rate. Further, mathematical calculations of waterflood behavior are restricted at this time by available mathematical techniques to highly simplified situations. Scaled models, however, do offer an approach at the present time to evaluation of the effects of the pertinent variables. Several investigators have predicted prototype reservoir behavior from observations of model waterflood performance under a variety of conditions. Two-dimensional aspects of waterflood pattern efficiency have been investigated by Muskat, Aronofsky, Dyes, et al, raig, et al,' and Rapoport, et al. he effects of gravity segregation, viscous fingering and permeability stratification have also been studied in two-dimensional models. Craig, et al," studied the importance of gravitational forces on frontal displacement in three-dimensional models with wells on a five-spot pattern. This work in- vestigates the tendency for water to segregate in homogeneous reservoirs because of gravity and to under-run the oil for a range of water-inject ion rates. The effects of permeability stratifications were investigated by these authors using miscible fluids to simulate water-channeling through permeable strata in systems of negligible capillary force. In this investigation, a three-dimensional model of a five-spot pattern flood was designed to represent typical prototype conditions, and the sand and liquid properties were selected so that capillary, gravitational and viscous pressure gradients were scaled. It was not, however, possible to achieve exact similitude regarding viscous fingering when the more viscous oil was used. The effects of rock stratification on waterflood behavior were studied by changing both the permeability ratio and the position of the more permeable sand stratum in a two-layered, horizontal model. For each type of stratification, constant-rate floods were performed and the rates were varied, in separate experiments, over a tenfold range. To observe the effect of mobility ratio on flood performance, two series of constant-injection-rate floods were performed on the most highly stratified reservoir at a water-oil viscosity ratio of one.

By R. V. Higgins, A. J. Leighton

A fast method is needed for calculating thoroughly the performance of two-phase flobr~ in reservoir rock with complex geotrietry. The authors present such a method and .slzow its accuracy by romparison with laboratory performrrnce of a five-spot pattern. In making the performance c.alculations, the octant of the five-spot pattern was divided into four channels. The computer tinle required for the culculation was one rllinrcte. This is rorighly about minute per channel. INTRODUCTION By using an electronic computer it is now possible to make reservoir performance calculations in which fewer assumptions and more significant variables are included than was possible with a desk calculator. These performance calculations include such variables as the properties of the reservoir rock and its contained fluids, properties of the rock-fluid system and even different well-spacing patterns. Usually, the more variables that are included the higher the computer cost—and time on the computer is expensive. As a result, a short procedure is needed so the many variables will be included and accurate performance calculations will result. If a reservoir is composed of alternate layers of permeable and impervioi~s rock, and if each conducting layer has a different relative permeability, the combined computer cost for calculating each layer by the present long computer programs (which includes different well-spacing patterns and other variables) would be expensive. The performances of complex spacing patterns in a field often are studied by a potentiometric model, but in these model studies the analogy is to the flow of only a single phase. This paper presents a procedure to calculate the performance of a five-spot pattern with two phases flowing using a potentiometric model as a guide for the dimensions of channels and other related data. The time required to calculate the complete performance is about one minute on the IBM 7090 or four minutes on the IBM 704. The procedure can be extended to calculate the performance of more complex well spacings or an entire field with little modification. The procedure has been tested on the performance of the laboratory water floods of a five-spot pattern reported by Douglas, et alThe range of viscosity ratios in the laboratory experiments is from 0.083 to 754, which is extensive. The accuracy of the performance calculations is excellent as shown by the plotted points in Fig. I. The recovery as influenced by the geometry of the five-spot, taking into account the deviation from true radial character, has been predicted by Aronofsky and Ramey." The authors used a potentiometric model to determine the sweep efficiency to breakthrough of the flooding agent. Dyes, Caudle and Erickson,' by means of a laboratory model, show the influence of mobility ratio on recovery after breakthrough of the injected material. Dyes, et ai, used the X-ray shadowgraph technique and miscible phases in their studies. Craig, Geffen and Morse from the analysis of the performance of their laboratory model correlate many variables, including average saturation of a flooding agent determined by permeability relationship. Craig, et al, used X-ray shadowgraphs to determine the sweep efficiencies in sandstone models. In this report the sweep efficiency is not needed, and the important effect of continuously changing saturation gradients before and after breakthrough

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst's horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter's edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the

By J. P. Heller, H. B. Bradley, A. S. Odeh

A potentiometric model technique is presented for determining the areal sweep efficiency of a five-spot well pattern, at and beyond breakthrough. A sharp interface between displaced and displacing fluids is assumed. Although the prototype svstem is referred to as a waterflood operation, obvious changes in the notation and computations will adapt the results to other displacement processes. The results of the study include the following. 1. Areal sweep efficiencies for a five-spot well pattern, at and beyond breakthrough, for mobility ratios (displacing to displaced fluid) of 4:1, 2:1, 1:1 and 1:4. 2. Extension of potentiometric analysis to the investigation of the areal sweep efficiency beyond breakthrough in a five-spot pattern by the application of conformal mapping and conductive-solid models. 3. The use of layers of conductive fabric in representing mobility ratio changes in potentiometric models. 4. The development of a probe mechanism for probing conductive solids. The results obtained by conductive-cloth models agree with earlier areal sweep efficiencies at breakthrough obtained by Aronofsky and Ramey 1 on the potentiometric analyzer using electrolytic-tank models. Results beyond breakthrough differ from those obtained by the X-Ray Shadowgraph technique. Data from this study show that, for mobility ratios greater than one, water cut rises rapidly as fluid is produced after breakthrough. However, for mobility ratios smaller than one, a large increase in area swept resulted with only a small increase in water cut. INTRODUCTION In calculating reservoir performance of waterflood-ing operations and other fluid-injection programs, it is necesssary to estimate the areal sweep efficiency * before and after injected-fluid breakthrough into production wells. Influence of mobility ratio on oil-production history, before and after breakthrough for a five-spot well pattern, has been studied by X-Ray Shadowgraph techniques 2-5 and by gelatin models.6 In none of these investigations was the transition zone controlled experimentally. Steady-state anaiog techniques assume that a vertical and discrete interface exists between the displacing and displaced phase. Thus, when using a potentiometric analysis, it is assumed that no transition zone exists. To determine the areal sweep efficiency beyond breakthrough for a five-spot well pattern using the potentiometric analyzer, four principal problems had to be resolved. These were the following. 1. The laborious method of representing changes in mobility ratio in electrolytic tanks by means of contoured wooden blocks. 2. Vapcrization of the electrolyte solution during the study causing concentration changes and, thereby, variable conductivity. 3. The scale-limitation problem 7 that is accentuated by abrupt eiectrolyte depth changes. 4. A wellbore geometry problem. EXPERIhlENTAL EQUIPMENT AND PROCEDURES CONDUCTIVE-SOLID MODEL The theory, equipment and procedural techniques of the electrolytic-tank potentiometric model have been treaced excellently by Lee8 and by Muskat.9 By the use of a conductive solid to represent the fluid conductivity of the porous media, the need for contouring wooden blocks to represent mobility ratio was eliminated. It was found that a carbon-black coated fabric, Uskon D-16,** could be employed as a fairly uniform conductive media. This material, when stacked in appropriate layers and pressed together, can be used to represent changes in mo-

By D. H. Henley, F. F. Craig, W. W. Owens

The oil recovery performance of systems producing entirely by bottom-water encroachment has been experimentally determined in a series of scaled laboratory-model tests. The effects of well spacing, fluid mobilities, rate of production, capillary and gravity forces, well penetration and well completion techniques on the oil recovery performance have been investigated. The laboratory tests were performed using two uniform, un-consolidated sand-pack models. The models have ratios of the interwell distance to the formation thickness of 12 and 2, respectively. Tests at constant total fluid production rate were performed simulating a range of uniform reservoir characteristics and operating conditions encountered in field operations. The performance was determined by material balance and by observation of the encroachment of dyed fluids into the models. The results of the model tests agreed with those obtained mathematically when the conditions previously considered in theoretical studies were simulated, that is, when the oil and water are of equal density and no capillary forces exist. The model study of bottom-water drive indicated that certain variables can aoect the oil recovery performance to a greater degree than can be predicted by present analytical methods. In one comparison, the oil recovery at a water-oil ratio of 20 (obtained at a wide well spacing) varied as much as threefold, depending upon the system's properties and the production rate. Lesser effect of mobility ratio and no eflect of capillary forces over the range studied were observed. The test 'results also showed that the deeper the well penetration into the oil column, the greater the total water production to a producing WOR of 20. However, the ultimate sweep efficiency, and so the oil recovery to this level of WOR, did not vary significantly with well penetration. Horizontal fractures at the top of the formation did not significantly change the sweep characteristics of the reservoir models when values of radius and fracture capacity encountered in actual reservoirs were used. Impermeable pancakes at the bottom of the well moderately increased the oil recovery efficiency both at water breakthrough and at high water-oil ratios. A method is outlined by which the oil recovery performance of other uniform bottom-water drive systems can be estimated from the information obtained in these model tests. INTRODUCTION When oil is produced from a well which partially penetrates an oil zone completely underlain by water, the water rises directly beneath the well in a symmetrical cone when the system is uniform. Two different flow mechanisms can cause the water cone to form—coning and bottom-water drive. In coning, the aquifer is relatively inactive and the cone is formed beneath the well by the pressure gradients associated with the oil flow to the well. The oil can be produced by a solution-gas drive, an edge-water drive or other driving forces in the interwell area. In a bottom-water drive, the driving force for oil production comes from an upward encroachment of the underlying active aquifer. Two papers have analyzed the theoretical performance characteristics of bottom-water drive reservoirs. In the initial mathematical investigation, Muskat' established the equations which determine the pressure distribution in this type of reservoir and solved these equations for certain conditions. Specifically, it was assumed that the water and oil had equal mobilities and equal densities, there were no capillary forces, the pressure throughout the oil zone remained above the bubble-point pressure, a constant pressure existed at the initial water-oil contact and the oil was completely displaced by the encroaching water. These assumptions were used in obtaining analytical solutions. In general, Muskat found that the sweep efficiency to initial water breakthrough to the well was larger for the thicker oil zones, the closer well spacings, the lower ratios of vertical to horizontal permeabilities, the smaller the penetration of the well into the oil zone and the smaller the bore size of the well. The production history after water breakthrough was expressed as a volumetric sweep efficiency at a given producing water-oil ratio. The results indicated that cumulative oil production at producing water-oil ratios of 10 is less affected by the well spacing than is the water-free production history. Muskat studied well spacing which today would be regarded as close. The maximum value of his dimensionless well spacing (ratio of interwell distance to formation thickness) was 4.3. This would require the development of a 50-ft-thick oil sand on less than 10-acre spacing if the vertical and horizontal permeabilities were equal, with

By F. Brons, W. C. Miller

Pressure information for use in material-balance calculations is obtained, where possible, from pressure build-up surveys in shut-in wells. Using proper extrapolation methods, static pressures are obtained which, by averaging, give the field static pressure. Often, particularly in fields with a large number of wells, only spot pressure readings are available. A graphical method is presented here which enables the determination of the corrections which must be applied to the spot readings to give an estinzate of the static pressure. Although the static pressures thus obtained are not so reliable as those from build-up surveys, they are more nearly correct than the spot readings themselves for use in material-balance calculations. INTRODUCTION The average static pressure in a bounded reservoir is determined from a number of individual-well pressure measurements. Where pressure build-up curves are available, extrapolation methods such as described in Ref. 3 may be used to determine the static pressure in the drainage area of each individual well. The possibility of using pressure build-up curves presupposes that each well being surveyed is closed in during a certain period and that, during this time, a pressure recorder is left suspended in the hole. For surveys in fields with a large number of wells, this is usually not feasible for economic and other reasons. Either a large number of pressure recorders would have to be used simultaneously or, using a limited number of recorders, the time interval needed to complete a field-wide survey would be too long. The data obtained for this incomplete coverage are insufficient to permit reliable averaging of the few calculated static pressures over the entire reservoir. In many such cases another type of pressure information is available, i.e., spot pressure readings. During this type of pressure survey all wells are closed-in simultaneously and, after a reasonable shut-in time (e.g., three days), pressure readings are taken by running the recorder in and out of each successive well. These surveys, which are less expensive, are often taken at regular intervals, say once or twice a year. In fact, such surveys are required by the state regulatory bodies in many portions of the united States. These spot readings are normally lower than the static (fully built-up) pressure due to the relatively short shut-in time. No direct means of extrapolating such a reading to the static pressure is available because no pressure build-up record is obtainable. Consequently, these data are inferior to pressure build-up information as a basis for estimating the stabilized reservoir pressure in the drainage area of a specific well. For the common case where abundant spot readings but only a few pressure build-up surveys are available, our studies suggest that it may be possible to design a program which will combine the two types of information and provide an estimate of the average reservoir pressure that is more reliable than that obtainable solely from either set of data. in time) gives a point marked p. The actual static pressure at the time the well was shut-in is marked by p.. The illustration shows that, even if a build-up curve is available and can be extrapolated, a correction

By J. D. Testerman

A statistical technique to identify and describe naturally occurring zones in a reservoir and to correlate these zones from well to well is described. The technique is particularly useful in describing a reservoir where cross flow between adjacent strata is important in determining reservoir behavior. Although it has been used primarily for permeability zonation, the technique is general and can be used to correlate any reservoir property or related data, such as the information contained in well logs. INTRODUCTION One of the first problems encountered by the reservoir engineer in predicting or interpreting fluid displacement behavior during secondary recovery processes is that of organizing and using the large amount of data available from core analysis. Permeabilities pose particular problems in organization because they usually vary by more than an order of magnitude between different strata. Due to the sheer volume, it is almost always necessary to group data and to use an average value to represent a number of measurements. Perhaps the most common method now used to group permeability data is the capacity-fraction technique, which ranks permeabilities in order of magnitude, regardless of the physical location of the permeabilities within the reservoir. The cumulative per cent capacity is plotted against cumulative per cent thickness. This plot is divided into an arbitrary number of zones, generally of equal thickness. Five zones (or averaged groups of data) usually are obtained, each of which is treated as homogeneous in subsequent calculations. The division so obtained has no physical meaning; strata in the same zone, calculation-wise, are usually not adjacent in the reservoir. Reservoir engineering techniques being developed will handle crossflow that occurs between adjacent communicating reservoir strata because of imbibition and gravity segregation. Since crossflow occurs between physically adjacent layers within the reservoir, a new zonation technique recognizing the actual location of strata within the reservoir is necessary. Similarly, the recognition of natural zones is important for predictions of oil recovery by processes involving diffusion. One such process is miscible displacement, where predictions of lateral diffusion within the reservoir must recognize the actual location of the invaded zones in relation to the rest of the formation. Natural zones must also be adequately recognized to account for heat transfer within the reservoir during thermal exploitation. Because of the complexity of the problem, statistics appear to offer the only practical hope of dividing a reservoir into physically-meaningful natural zones. This paper presents a statistical technique for identifying these natural zones and for ascertaining which ones are likely to be continuous between adjacent wells. The zones defined have minimum variation of permeability internally and a maximum variation between zones. The technique is general and can thus be applied to reservoir properties other than permeability. The method will guide the reservoir engineer in estimating which zones are likely to be continuous between wells. However, a statistical correlation based on permeabilities in two different wells is no guarantee that the zones so defined are, in fact, continuous. Rather, the assumption of continuity must be consistent with geological data concerning the depositional environment, as well as justified on the basis of engineering judgment in combination with statistics, just as judgment is required with conventional zonation methods. CALCULATION PROCEDURE The reservoir zonation technique is a two-step operation. The steps are individually described, and a sample calculation is presented in the Appendix. ZONATION OF INDIVIDUAL WELLS First, the set of permeability data at a single well is zone, into Zones. These zones are selected so that variation is minimized within the zones and maximized between the zones. The equations4,6 used to zone the data are where B = the variance between zones, , = the number of zones, i = the summation index for the number of zones, j = the Sumation index for the number of data within the zone, mi = the number of data in the ith zone, k,. = the mean of the permeability data in the ith zone, k . = the over-all mean of the data in the well, W = the pooled variance within zones, N = the total number of

By N. T. Cotman, C. H. Pickering, P. B. Crawford

A beat-conduction model has been developed to study the flow of fluids in a stratified oil reservoir which is being subjected to unsteady-state depletion. To simulate stratification, plates of different metals were joined together so that cross flow could occur between strata. A description of the model, plate construction and necessary instruments is included. The data indicate that stratified oil reservoirs experiencing cross flow may be treated as a single, homogeneous, producing sand having properties intermediate between those of the layers making up the system. A technique for averaging core-analysis data to arrive at the proper average value of rock permeability and porosity is presented INTRODUCTION Although oil has been produced from reservoirs for 100 years and reservoir injection processes have been undertaken for 40 years, the industry is still inadequately informed on the problem of reservoir inhomogeneities. Stratification is perhaps the best-known type of reservoir heterogeneity. Urenl discussed the possible significance of this nonuniformity in his 1927 paper on the theoretical aspects of water flooding. Since this time the problem of vertical permeability variations has been the subject of much discussion and numerous papers. Levorsen2 simply describes the origin of layered beds in his Geology of Petroleum. Some beds are field-wide and may be traced from well to well, while others apparently pinch-out and become discontinuous. In addition, the physical characteristics of a layer may vary from place to place. For these reasons, a stratified system is not easily defined. In calculations involving oil reservoirs in which permeabi1ities vary, methods presently used contain many simplifying assumptions. The reservoir frequently is considered as consisting of permeable layers, each of uniform vertical and lateral permeability and each behaving as if separated from adjacent layers by thin impermeable layers. Other assumptions frequently made in reservoir calculations are steady-state flow conditions, unit mobility ratio, linear geometry, average or homogeneous porosity, and negligible capillary-pressure and gravity effects. Muskat3 has. developed solutions to the stratification problem for mobility ratios other than unity for both linear and exponential permeability distributions. He assumes lateral uniformity and continuity of all productive strata and does not consider porosity. This study will take into account transient effects, cross flow and varying porosity, as well as the variations in permeability, by means of a heat analogy. Muskat5 has shown the similarity of the equations which describe the flow of fluids in a porous medium to the flow of heat through a solid. Landrum, et al,6 made an analogy between the flow of fluid in an unsymmetrical reservoir and the conduction of heat in a similarly shaped metal plate. This analogy is summarized briefly as follows. If fluid is produced from an oil reservoir above the bubble point, the resulting pressure distribution will be identical to the temperature which would result in a metal or heat-conducting plate if heat were similarly removed. This means that a study of heat flow in stratified thermal models will behave like fluid flow in stratified oil reservoirs. For mathematical correspondencies, set- Ref. 6. The fluid mobility k/p is proportional to the thermal conductivity, and the porosity-times-compressibility product is proportional to the density times specific heat of the thermal plate. EXPERIMENTAL EQUIPMENT AND PROCEDURE The stratified reservoir model used in this study is illustrated schematically in Fig. 1. A Leeds and Northrup Speedomax Type G temperature recorder

By C. D. Stahl, E. E. Templeton, R. F. Nielsen

It has been customary, in predicting saturation changes, to use the Leverett fractional flow formula", obtained by eliminating the unknown pressure gradient from the generalized Darcy equations for the separate phases. The formula presents difficulties in the case of counterflow, since the "fractional" flow may be negative, greater than unity, or, in the case of a closed system, infinite. Recently, it has been shown by several authors that the corresponding equations (with capillary pressure and gravity terms) for actual flow of the phase may be used just as well. These equations are in agreement with Pirson's1 statement that, if the two mobilities differ considerably from each other in a closed system, the flow is largely governed by the lower value. The present study was undertaken because of an apparent lack of experimental data on gravity counter-/low with which to test the theory. A 4-ft sand-packed tube in a vertical position was employed. Electrodes for determining saturations by resistivity were spaced along the tube, one phase being always an aqueous salt solution. Air, heptane, naphtha, or Bradford crude oil was used for the other phase. A reasonably uniform initial saturation was set up by pumping the phases through the system, alter which the tube was shut in and saturation profiles obtained at definite intervals. Cumulative flows over certain horizontal levels were obtained by integration of the distributions; hence, differentiation of the cumulative flows with respect to time gave instantaneous flow rates. To compare experimental and theoretical flow values, capillary pressures were assumed given by the final saturation-distribution curve. The upper part corresponds to the "drainage" region and the lower part to the "imbibition" region, where trapping of the nonwetting phase occurred. While calculations indicated that the capillary pressure saturation function and, probably, the relative permeability saturation functions changed during the segregation, the relation of the measured rates to saturation distributions are in general accord with the frontal-advance equation. It appears that the Darcy equations, as modified for the separate phases, are generally valid for counterflow due to density differences. The usual method of predicting saturation changes, which involves a continuity equation and the elimination of the unknown pressure gradient from the flow equations, should therefore be applicable. However, the need for advance knowledge of drainage and and imbibition "capillary pressures" and relative permeabilities during various stages presents difficulties. INTRODUCTION The present study was undertaken because of a seeming lack of experimental data relating to vertical counterflow of fluids of different densities in porous media. In particular, it was desired to determine whether data obtained from these laboratory tests were in accordance with certain mathematical treatments of counterflow which have been proposed. The gravity "correction" has been incorporated into the flow equations (and, hence, into displacement theory) nearly as long as both have been used. field4 and laboratory5 data have generally borne out the validity of the theory as applied, for instance, to downward displacement by gas, with all fluids moving downward. However, the modifications for counterflow have only recently been pointed out. 2-3 It has been customary to use fractional flow rates instead of actual flow rates in displacement calcuIations. In the case of counterflow, this results in negative values, values greater than unity and, when rates are equal and opposite, in infinite values. As pointed out by Sheldon, et a1,3 and by Fayers and sheldon,2 actual flow rates may be used just as well. The fact that these may be of opposite signs for the two fluids does not present any difficulty.

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By E. F. Johnson, F. H. Brinkman, H. J. Welge, A. L. Hicks

A method is presented for predicting the character of gas or water displacement in a radial system, which can be either horizontal or inclined. The latter case would comprise cone-shaped or dome-shaped symmetry. The method conriders the one dimension of radial distance, taken parallel to the bedding planes. Approximate allowance can be made for segregation of the two fluids perpendicular to the bedding planes. A complete saturation history can be obtained in about a day's time with the use of a desk calculator. Results of an example calculation agree well with the performance of simulated gas drive in a large model representing a section of a cone. This pie-shaped, sand-packed model was 8-ft long and 4-in. thick. The method of solution is useful in handling other partial differential equations of first order and first degree. INTRODUCTION A new analytical method has been developed for predicting the performance of gas or water drive in a reservoir having radial or cone-shaped symmetry (see Fig. 1). In addition, this new development offers a general method of handling partial differential equations of first order and first degree. This work extends those methods developed previously'- "or predicting the displacement of oil by gas or water in linear systems. The method of calculating oil recovery in a radial system involves two major parts. The first part is a calculation of the distribution of saturations at various times of depletion. The second part is the calculation of oil recovery and producing gas-oil or water-oil ratios through volumetric balance on the saturation distributions obtained in the first part. This second part differs only slightly in principle from similar calculations made with linear systems, for instance. The basic assumptions used in the calculations are as follows. 1. One-dimensional radial flow. An approach is explained later for approximating a second vertical dimension. (With- drawals of oil, and of the driving fluid after breakthrough, are assumed uniform around the periphery of the dome.) 2. Constant-volume fluid flow. Thus, in a gas drive there is no change in gas solubility, or in oil or gas volume due to viscous and gravity pressure drops. 3. Capillary forces are negligible. 4. Constant rate of injection. A program has been written for carrying out the calculations needed on a digital computer. However, a prediction can also be arrived at using a desk calculator; this latter method of work requires about a day's time for each prediction made. CALCULATION OF SATURATION DISTRIBUTION For definiteness, we describe the development of our method for the case of gas drive from the crest of the structure. The reader will note, however, that the method is just as applicable to an advancing radial water drive. The equation describing gas flow over an infinitesimal