Search Documents

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

By L. G. Jones

An approximate method of calculation is developed in this paper which allows duplication of radial unsteady-state gas flow computer results where Darcy's law applies, such as those reported by Aronofsky and Jenkins and Bruce, Peaceman, Rachford, and Rice. Moreover, the new calculation method can be used to obtain results for radial unsteady-state gas flow obeying the quadratic flow law proposed by Duwez and Green. Means are discussed for predicting well behavior at single or superimposed flow rates in finite or infinite reservoirs, determining reservoir rock properties from well-test data, reproducing and interpreting back-pressure test data, and determining the radial extent and reserves of gas reservoirs from well-test data. Example calculations are presented for gas flow following both Darcy's law and the quadratic flow law. INTRODUCTION Since the publishing of U. S. Bureau of Mines Monograph 7,1 most gas-well testing methods have been based on the equation q = y (pl2 - pw2)a, where q = production rate, pl and pw, are formation pressure and sandface pressure, respectively, and y and a are constants to be obtained from test data. These methods, used for predicting both deliverability and "open-flow" capacity of gas wells, have been useful and accurate in many cases but unsatisfactory in others. Even at best, however, they do not supply information about the formation or lead to an understanding of nonsteady-state gas flow in porous media. Many theoretically based studies2-4 of gas flow obeying Darcy's law have been made.' Since the partial differential equations which result from combining Darcy's flow law with the continuity equation are nonlinear, most of the published research consists of either numerical solutions or analytical solutions for linear approximating equations. Such solutions have been of limited value in field work due to their unhandy form and their failure to correlate most field data. There is evidence5,6 which indicates that Darcy's law is inadequate to describe gas flow at some flow rates of practical interest. A quadratic flow law: which reduces to Darcy's law at low rates, is more successful in accounting for experimentally observed behavior. This flow law has been applied successfully to a few hypothetical reservoir cases in work which has recently been published.10 However, these numerical solutions of the equations involved have been successful only on a special analog computer. Routine applications to field cases would be awkward and have not been attempted. The present paper describes an approximate method for computing nonsteady-state gas flow solutions which has been completely successful in predicting the results for both Darcy flow and quadratic flow obtained by elaborate numerical methods. The new calculation method allows determination of the observable variables in gas-well testing at constant rates. It is similar to the scheme of using a succession of steady states suggested by Muskat7 in that it makes use of steady-state and material-balance equations. It also is similar to the work of Aronofsky and Jenkins2 in that the new method includes Aronofsky and Jenkins' stabilized flow equation as a special case. It improves upon both of these calculation schemes in that it accurately describes all portions of reservoir history and suggests means of determining reservoir rock properties from well-test data. This paper deals only with production from a reservoir, in which case the rate is defined as being negative. The reservoir studied here is a homogeneous, disk-shaped porous body of uniform thickness, with all boundaries sealed except the inner radial boundary of the well.

By G. W. Nabor, M. Mortada

The effects of anisotropic or directional permeability on the areal sweep efficiency and the flow capacity are examined. The paper points out the importance of taking directional permeability into consideration in planning a flood. It analyzes the two-dimensional flow pattern associated with the skewed line drive for a unit mobility ratio. The direct and staggered line drives are treated as special cases of the skewed line drive. Analytical expressions are developed for the areal sweep efficiency at breakthrough and the flow capacity. They are related to the spacing between like wells, the distance between a row of injectors and the nearest row of producers, and the degree of skewness of the line drive. The latter quantity is defined such that it is equal to zero for the direct line drive and equals one-half for the staggered line drive. The areal sweep efficiency and the flow capacity depend also on the orientation of the flood pattern with respect to the principal axis of anisotropy. The paper provides a simple method for determining the areal sweep efficiency and the flow capacity for a formation in which the permeability in the bedding plane is anisotropic. INTRODUCTION Directional or anisotropic permeability is manifested by the ability of the formation to conduct fluids more readily along certain preferred directions. This situation occurs in many producing formationsl-4 and is usually attributed to depositional features in which the sand grains are oriented in a preferred direction. In some cases it results from the formation of a major and a minor fracture system. 4 Directional permeability should be taken into account in many phases of the production and exploration activities. Recognizing its existence in the formation of interest and planning accordingly can lead to increased recovery and substantial savings. For instance, the areal sweep efficiency in a water flood depends to a great extent on the orientation of the flood pattern with respect to the principal axis of permeability. Anisotropic5 permeability is specified by the directions of its three principal axes and the permeability along each axis. The principal axes of permeability are mutually perpendicular. This paper deals with the areal sweep efficiency at breakthrough and the flow capacity for formations with aniscltropic permeability. The flood pattern considered consists of alternate rows of injecting and producing wells. The rows of wells are parallel and form a developed, skewed line drive which is illustrated in Fig. 1. The staggered and direct line drives are treated as special cases of the skewed line drive. PRELLMINARY CONSIDERATIONS The formations considered in this paper are assumed to be of uniform thickness and porosity. The permeability is uniformly anisotropic, with each principal axis of permeability maintaining

By G. Rowan, M. W. Clegg

The basic equations for the flow of gases, compressible liquids and incompressible liquids are derived and the full implications of linearising then discussed. Approximate solutions of these equations are obtained by introducing the concept of a disturbed zone around the well, which expands outwards into the reservoir as fluid is produced. Many important and well-established results are deduced in terms of simple functions rather than the infinite series, or numerical solutions normally associated with these problems. The wide range of application of this approach to transient radial flow problems is illustrated with many examples including: gravity drainage of depletion-type reservoirs; multiple well systems; well interference. INTRODUCTION A large number of problems concerning the flow of fluids in oil reservoirs have been solved by both analytical and numerical methods but in almost all cases these solutions have some disadvantages — the analytical ones usually involve rather complex functions (infinite series or infinite integrals) which are difficult to handle, and the numerical ones tend to mask the physical principles underlying the problem. It would seem appropriate, therefore, to try to find approximate analytical solutions to these problems without introducing any further appreciable errors, so that the physical nature of the problem is retained and solutions of comparable accuracy are obtained. One class of problems will be considered in this paper, namely, transient radial flow problems, and it will be shown that approximate analytical solutions of the equations governing radial flow can be obtained, and that these solutions yield comparable results to those calculated numerically and those obtained from "exact" solutions. It will also be shown that the restrictions imposed upon the dependent variable (pressure) are just those which have to be assumed in deriving the usual diffusion-type equations. The method was originally suggested by Guseinov,1 who (a) postulated a disturbed zone in the reservoir, the radius of which increases with time, and (b) replaced the time derivatives in the basic differential equation by its mean value in the disturbed zone. In this paper it is proposed to review the basic theory leading to the equations governing the flow of homogeneous fluids in porous media and to consider the full implications of the approximation introduced in linearising them. The Guseinov-type approximation will then be applied to these equations and the solutions for the flow of compressible and incompressible fluids, and gases in bounded and infinite reservoirs obtained. As an example of the application of this type of approximation, solutions to such problems as production from stratified reservoirs, radial permeability discontinuities; multiple-well systems, and well interference will be given. These solutions agree with many other published results, and in some cases they may be extended to more complex problems without the computational difficulties experienced by other authors. THEORY In order to review the basic theory from a fairly general standpoint it is proposed to limit the idealising assumptions to the minimum necessary for analytical convenience. The assumptions to be made are the following: 1. That the flow is irrotational. 2. That the formation is of constant thickness. 3. Darcy's Law is valid. 4. The formation is saturated with a single homogeneous fluid.

By J. H. Henderson, J. Naar

The displacement of a wetting fluid from a porous medium by a non-wetting fluid (drainage) is now reasonably well understood. A complete explanation has yet to be found for the analogous case of a wetting fluid being spontaneously imbibed and the non-wetting phase displaced (imbibition). During the displacement of oil or gas by water in a water-wet sand, the porous medium ordinarily imbibes water. The amount of oil recovered, the cost of recovery and the production history seem then to be controlled mainly by pore geometry. The influence of pore geometry is reflected in drainage and imbibition capillary-pressure curves and relative permeability curves. Relative permeability curves for a particular consolidated sand show that at any given saturation the permeability to oil during imbibition is smaller than during drainage. ' Low imbibition permeabilities suggest that the non-wetting phase, oil or gas, is gradually trapped by the advancing water. This paper describes a mathematical image (model) of consolidated porous rock based on the concept of the trapping of the non-wetting phase during the imbibition process. The following items have been derived from the model. 1. A direct relation between the relative permeability characteristics during imbibition and those observed during drainage. 2. A theoretical limit for the fractional amount of oil or gas recoverable by imbibition. 3. An expression for the resistivity index which can be used in connection with the formula for wetting-phase relative permeability to check the consistency of the model. 4. The Iimits of flow performance for a given rock dictated by complete wetting by either oil or water. 5. The factors controlling oil recovery by imbibition in the presence of free gas. The complexity of a porous medium is such that drastic simplifications must be introduced to obtain a model amenable to mathematical treatment. Many parameters have been introduced by others in Iprogressing" from the parallel-capillary model to the randomly interconnected capillary models independently proposed by Wyllie and Gardner3 and Marshall4. To these a further complication must be added since an imbibition model must trap part of the non-wetting phase during imbibition of the wetting phase. Like so many of the previously introduced complications, this fluid-block was introduced to make the model performance fit the observed imbibition flow behavior. THE IMBIBITION MODEL ASSUMPTIONS The following assumptions are common to many two-phase models previously described. 1. The wetting characteristics of the medium are clearly defined. 2. When two fluids co-exist in the medium, the wetting phase will fill the smaller pores and the non-wetting fluid the bigger ones. 3. The pores have a circular cross section. 4. Apparent pore radii are related to the capillary pressure (drainage) by

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By H. N. Hall

Various factors must be considered in an engineering evaluation of gravity-drainage reservoirs. Among these are: (1) the effect of producing rate on total oil recovery; (2) the effect upon well productivity and ultimate recovery of the pressure level maintained during the producing life of the reservoir; (3) the economic advantage of full or partial pressure maintenance; and (4) estimate of the rate of gas production and injection and the possible purchase of gas under conditions of full pressure maintenance to ascertain compressor facilities needed. All of these factors can be evaluated only when a reliable method is employed for determining reservoir performance in gravity-drainage reservoirs. The purpose of this paper is to present a general method for calculating the performance of a gravity-drainage reservoir. This method is applicable for conditions of complete pressure maintenance, partial pressure maintenance and normal pressure depletion. Provisions are made to take into account variations throughout the reservoir of reservoir configuration, changes in permeability and fluid composition. Based on the method presented in this paper, an IBM 650 computer program has been developed. The past performance of an actual gravity-drainage reservoir producing under conditions of declining pressure and no gas injection was duplicated using this program. INTRODUCTION In tilted reservoirs the production of oil is influenced by drainage of oil from upstructure to downstructure locations. When this downstructure drainage of oil is sufficient to cause effective segregation of the gas and oil in a reservoir, the reservoir is usually classified as a segregation drive or gravity-drainage reservoir. (Discussion will be restricted to gravity-drainage reservoirs which have no encroachment of edge water.) The important feature in gravity-drainage reservoirs is the density difference between reservoir oil and gas. These phases tend to segregate in the reservoir with the result that in the gas cap the oil saturation is maintained at a higher level by drainage of oil from the gas-cap area. Oil can be produced from the oil zone at a low gas-oil ratio and reservoir energy is thereby conserved. The standard material balance in not adequate for predicting gravity-dramage reservoir performance because it does not take into account the difference in saturation above and below the gas-oil contact. Several authors'.' have presented methods for calculating the performance of gravity-drainage reservoirs in which reservoir pressure is maintained constant by gas injection into the gas cap. Using some simplifying assumptions, these methods can be employed with a desk calculator to give acceptable results. The problem of predicting the performance of gravity-drainage reservoirs under the conditions of declining reservoir pressure is many time more complex than that of constant pressure. fierefore, attempis to develop a method suitable for desk calculation have required excessively simplified assumptions. In the past several years, highspeed digital computers have become more widely available for reservoir engineering problems. These corn puters are well suited to problems such as the prediction of the performance of gravity-drainage reservoirs with pressure decline. Many of the simplifying assumptions necessary for hand computation can be eliminated so that a realistic approach to the gravity-drainage process can be made. CONCEPTUAL PICTURE OF OIL MOVEMENT IN GRAVITY-DRAINAGE RESERVOIRS Before attempting to develop an analytical treatment for conditions occuring in a gravity-drainage reservoir, a concept should be formed concerning the movement of fluids in the reservoir as oil is produced. A review of the literature'.' shows that it is customary to classify gravity-drainage operations into two categories—(1) with complete pressure maintenance, and (2) with declining pressure. The same line of reasoning will be followed in presenting the concept of the movement of fluids in the reservoir because it is easier to visualize the movement of fluids under conditions of complete pressure maintenance. After discussing complete pressure maintenance, an analogy will be made between that and the case of declining pressure. It should be kept in mind throughout that the final aim for the problem of solving gravity-drainage performance with digital computers will be to develop a general program for any kind of gravity-drainage reservoir. COMPLETE PRESSURE MAINTENANCE One feature which is generally common in gravity-drainage reservoirs is a gas cap located at the top of the structure. This is shown in Fig. 1(a). Fig. l(b) shows oil saturations that might occur through the reservoir. In the gas cap, oil saturation is lower than

By R. E. Cook

This work presents the development and application of equations of the form developed by Martin1 to describe gravity segregation performance during natural depletion. One-dimensional depletion analyses are applied to several hypothetical reservoirs. Considerable attention is given to the changes in saturation distribution in the dip direction resulting from depletion at various rates. Vertical gas saturation distribution within a sand having vertical permeability is studied to illustrate conditions which provide For low gus-oil ratio production from wells selectively completed away from sand tops, as observed from well performance. Results of performance calculations are shown graphically and illustrate the rate sensitivity of the gravity segregation mechanism and the influence of reservoir geometry and withdrawal distribution. It is shown that for sufficiently high depletion rates gas saturation distribution in the dip direction is characterized by the formation of two discontinuities, or saturation "fronts"; the updip front being the gas-oil contact of the secondary gas cap and the downdip front occurring at a position dependent upon the reservoir depletion rate, withdrawal distribution and reservoir geometry. The gas saturation in the region between fronts is shown to continually increase with time, while that downdip from the lower front remains very low and relatively stable. The effect of increased depletion rate is to reduce the region of stable gas saturation and thus enlarge the portion of the reservoir producing with continually increasing gas saturation. Special attention is given to the performance of reservoirs having sufficient vertical permeability to permit liberated gas to segregate vertically within the sand section before proceeding in the updip direction. It is shown that reservoirs of this type are far less sensitive to depletion rate due to the mechanical control of producing gas-oil ratios afforded by selective well completions away from sand tops. It is shown that gravity segregation performance is adversely af- fected by decreasing reservoir cross-section area in the updip direction. It is also shown that down-dip withdrawal concentrations can restrict gas segregation rates in the dip direction. INTRODUCTION The natural depletion performance of several large producing reservoirs in Western Venezuela has unmistakably shown the significant role that gravity forces play in governing the gas-oil ratio behavior of individual wells and the over - all pressure-production performance of the reservoirs. Analysis and prediction of the natural depletion performance of these reservoirs is of great importance in determining the economics of pressure maintenance operations. In 1957, a method for the analysis of pressure maintenance performance was presented by Martin1 which, in the original and more general form of the equations involved, provides a physically sound basis for the analysis of one-dimensional gravity segregation performance during natural depletion. The work reported herein is an application and extension of Martin's outstanding developments to a series of hypothetical reservoirs in an attempt to explain the general behavior of the gravity segregation mechanism and the influence of such parameters as depletion rate, reservoir geometry and withdrawal distribution. The study considers two types of flow: "distributed flow", where vertical permeability in the reservoir is zero and fluids flow only in the dip direction while uniformly distributed over the sand thickness; and "segregated flow", where vertical permeability is sufficient to permit gas to segregate against the sand top before proceeding updip. This latter type of flow has been recognized as an important mechanism active in the massive sand reservoirs of Western Venezuela. ANALYSIS OF ONE-DIMENSIONAL DISTRIBUTED FLOW SYSTEMS For an inclined reservoir having zero vertical permeability, the velocities of oil and gas flow at any point in the reservoir can be derived from Darcy's equations as

By H. Dykstra

Calculated wellbore pressures were obtained for parameters of radilcs ratio and permeability. In all cases bur two, after-production was allowed to occur for one day. The calculated pressure build-up data were compared with actual pressure build-up data from a condensate well. The paper discusses the comparison, gives reasons for, inability to match calculated (2nd actual data, and presents conc1usions derived from the comparison. INTRODUCTION Considerable information has been published on the decline in wellbore pressure for oil and gas fields resulting from production of fluids. Considerable information also has been published on pressure build-up for oil fields, but relatively little has been published for gas and condensate fields. Perrine' summarized the equations derived for pressure build-up in oil wells and showed how the different methods could be applied. On the other hand, very little theoretical information has been published on pressure build-up in gas or gas-condensate wells.2-4 Tracy,' by pointing out the similarity between the equation describing oil wellbore pressure decline and the equation describing gas wellbore pressure decline, showed how methods of pressure build-up analysis could be applied to a gas well. In analyzing pressure build-up data for gas or gas-condensate wells, it would be very desirable to compare actual pressure build-up data with calculated pressure build-up data. Such a comparison would lead to a better feeling for the quantitative picture of gas-well build-up curves and would help to establish a degree of confidence in a given build-up curve. This paper discusses the comparison and presents conclusions made from it. Calculated wellbore pressures are given for parameters of radius ratio and permeability. In all cases but two, after-production was allowed to occur for one day. The work was done in an attempt to evaluate the reserves of a gas-condensate field. As will be discussed later, it was not possible to make an evaluation. METHOD OF CALCULATION The method used for solving the depletion and pressure build-up behavior for a gas-condensate well was that developed by West, Garvin and Sheldon.5 or the IBM 704 computer program used, it was assumed that liquid condensing out of solution would have only a minor effect on the flow behavior. Therefore, the field was treated as having single-phase flow, with the gas saturation S, remaining constant at unity minus the connate-water saturation Sw. As will be discussed later, this may not have been a good assumption. It also was assumed that oil production could be considered as equivalent gas production by using a conversion factor based on an analysis of the liquid produced. An outer boundary condition of no-flow was assumed. BASIC DATA The data required for the study were reservoir properties, fluid properties and production data. Reservoir properties included the following: thickness h, 179 ft; porosity, 0.194; connate-water saturation, 0.43; gas saturation, 0.57; permeabilities k, 1.0, 0.5, 0.25 and 0.15 md; radius ratios re/rw, 800, 1,000 and 2,000; wellbore radius rw, 5 in.; initial reservoir pressure P, 6,529 psia; and reservoir temperature, 272°F. The selection of permeabilities was based on core-analysis data showing an average air permeability of about 2 md. With a relatively high connate-water saturation, the average effective gas permeability would be about one-half, or less, of the air permeability. The selection of radius ratios was based partly on the observed rapid decline in wellbore pressure during production and partly on the fact that the well was believed to be located in a relatively small fault block. Reservoir fluid properties were obtained from a reservoir fluid study on a recombined sample of gas and condensate. Gas formation-volume factors were calculated from pressure-volume measurements made at reservoir temperature. Gas viscosities were calculated by the method of Carr, et al,G from an analysis of the re-combined sample. These data are shown as a function of pressure in Table 1. The production schedule from the gas-condensate

By D. D. Smith, D. P. Squier, E. L. Dougherty

A system of differential equations describing the temperature behavior of fluid injected at constant surface temperature in a well is derived and .solved analytically. A formula for the fluid temperature at any time and depth is given, as well us a special formula valid for very large times. These formulas are used to calculate temperatures for several typical cases. The results indicate that, initially, the temperature of the water entering the formation is considerably lower than the injection temperature. This condition lasts for only a short period— less than three days for most cases of practical interest. Following this highly transient period, during which the temperature of the fluid entering the formation builds up to about 50 to 75 per cent of the injection temperature. the system enters a quasi-steady state in which the temperature changes are very slow. After severl years, the bottom-hole temperature will still be 50" to 100°F lower than the injection temperature, hilt the heat losses may he tolerable. INTRODUCTION Predicting the behavior of a hot-water flood requires that the temperature of the water entering the injection interval be estimated. This report describes the development and solution of a system of equations which describes the temperature behavior of the injected water in the wellbore with certain simplifying assumptions. The only previous means known to the authors for describing such a process is that of Moss and White.' Their results appear to be close to those obtained by our method in the practical cases which were compared; this agreement is largely due to the fact that in our method temperature soon approaches a quasi-steady state, as was assumed in their method throughout. However, our model covers all times, is continuous (whereas the Moss-White model depends on breaking the depth into discrete intervals) and. we feel. more closely describes the physical problem. FORMULATION OF THE PROBLEM PHYSICAL SYSTEM AND ASSUMPTIONS The injection procedure consists of pumping water at a fixed surface temperature T., down an infinitely long cylindrical well or tubing of inner radius Any material exterior to the water column such as mud, casing, or cement is regarded as part of the formation. The general behavior of the system may be described qualitatively as follows. When the hot water is first introduced into the system, the temperature difference between the formation and the water is large, resulting in a high rate of heat transfer. As a result, the temperature adjacent to the wellbore rises very quickly. Because the segment of the formation adjacent to the wellbore largely controls the heat transfer rate, the heat transfer rate will become relatively constant when this portion has reached a temperature close to that of the water opposite it. The temperature of the water and formation then increase very slowly with time. The length of the initial highly transient period and the temperature of the water at its conclusion will be functions of depth, injection rate, injection-string radius, surface injection temperature and the physical properties associated with the water-formation system. The following additional assumptions were made. 1. There is no heat transfer by radiation in the system. 2. There is no heat transfer by conduction in the vertical direction in either the injection stream or the formation. 3. The linear volumetric and mass flow rate of the water is constant throughout the injection stream. 4. No horizontal temperature gradient exists in the injection stream. 5. The product of density and heat capacity is constant for both the water and the formation, and the formation thermal conductivity is constant. 6. Initially, both the water in the wellbore and the reservoir are at a temperature given by the (constant) ambient surface temperature plus the product of depth and geothermal gradient (assumed constant). At large distances for the wellbore (r m), the formation will remain at this temperature. 7. The water temperature and the formation temperature at r — r,, are equal for all depths D. DERIVATION OF EQUATIONS The differential equation satisfied by the fluid temperature T,(D, t), which is obtained by writing a heat balance on a cylindrical differential of volume dV of the injection string between the depths D and D i dD, is

By C. F. Weinaug

In recent years, several methods have been proposed for correlation of physical properties and behavior of hydrocarbon systems. These correlations were found to be unsatis. factory for predicting differential gas liberation and shrinkage data now obtained experimentally. To achieve the desired accuracy, methods were developed to adjust the correlations with two experimental pieces of data. The required data are the saturation pressure at reservoir temperature and the specific volume or density at these conditions. The saturation pressure is used to adjust the vapor-liquid equilibrium-ratio correlations. The density measurement is used to determine the apparent density of liquid methane. With these adjustments, differential gas liberation and shrinkage are predicted within the experimental error. EQUILIBRIUM RATIOS Equilibrium ratios published by the National Gas Assn. of America' were used as a basis for the calculation. The convergence pressure for the sample was determined by using a hypothetical binary system, with methane as one component and the remainder of the system lumped together to form the other component. For the critical temperature of the heavy constituent, the weight average of the critical temperatures of the ethane and heavier components was used. The correlation of pseudo-critical temperature, density and molecular weight of Matthews, Roland and Katzl was used in this calculation for the hexane and heavier fractions. Points for the critical properties of methane and the hypothetical heavy constituent were placed on a chart of the critical loci of hydrocarbon systems. A critical loci for the hypothetical binary system was drawn, and the critical pressure of this system was determined at the reservoir temperature. This critical pressure became the convergence pressure for which equilibrium ratios were read from the NGAA4 charts. The NGAA equilibrium constants for C, through C, were used for the C,; through C, fractions since they were separated on the bases of normal paraffin boiling points. The correlation of Hoffman, Crump and Hocott5 as used to estimate the equilibrium constant for the C fraction. The NGAA eauilibrium-ratio charts are based primarily on paraffin systems and have to be corrected for the presence of naphthenes and aromatics. The correlations in the literature proved too unsatisfactory for determining these corrections. EXperimental data, the saturation pressure, is used with the correlation of Soloman2 to make this adjustment. The correction factors of Soloman were cross-plotted so that the factors for ethane and heavier components were made a function of the factor for methane. a2 = 0.7927 at + 0.1646 .........(1) where a, = correction coefficient for for methane NGAA equilibrium ratio, and a2t = correction coefficient for ethane and heavier NGAA equilibrium ratios. This equation and the equation for the bubble point are solved simultaneously for the correction coefficients. The bubble-point equation is anKnXn = 1.0 . . . . (2) The equilibrium constant for the saturation pressure at reservoir temperature is used. This procedure not only corrects for the relative volatility between components, but also assures that the results will be directly related to the experimental saturation pressure. METHANE DENSITY After determining a set of equilibrium ratios, the apparent density of condensed methane at 60°F and 14.7 psia was estimated. This was found to be necessary because the correlations of Standing3 ere found to be unsatisfactory. The law of additive volumes was employed to make this calculation.

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

Interest by petroleum engineers in the flow of three phases—oil, gas and water—in irregularly bounded porous media lies mostly in the performance calculation of water floods of reservoirs that have been partially depleted as the result of expansion of much of the originally dissolved gas. The authors present a method to forecast three-phase flow in complex geometry and explain the details by the use of a specific example of a five-spot water flood of a partially depleted stratified reservoir. In this example, the fluid and rock properties of a field given by Prats, et al,' were used. The computed results fit the field performance more accurately than Prats, et al, and Slider.' The time required for the high-speed digital computer to make the calculations, including the contributions from the different layered zones, is about one minute. INTRODUCTION The declining rate of discoveries of new oil fields in the United States makes the recovery of more oil from known reservoirs more attractive than previously. Water flooding has an excellent proved background for recovering additional oil economically. Accordingly, the main interest of this paper is in this type of recovery. In petroleum engineering studies, commercial interest in three-phase flow is mostly in the water flooding of reservoirs in which the oil has been partially produced by the expansion of dissolved gas. When water is pumped into these reservoirs, three-phase flow takes place. Although this paper is concerned with these conditions, the principles involved could be used for other conditions when and where they come to the fore. Many investigators have made contributions to non-empirical forecast methods using basic scientific engineering principles. Several of these used the oil in place at the start of the water flood and the oil remaining after a large quantity of water has passed through a core for key values in their calculations. Recently, Prats, et al,' and Slider have reduced assumptions by adding the third phase—gas—in their forecasting methods. Slider uses the mobilities in the immediate vicinity of inlet and outlet wells as an aid to simulate the resistance effect in the five-spat pattern of the flow of oil, gas and water. Prats, et al. minimize assumptions by using previously determined laboratory data for determining sweep efficiency in the five-spot pattern. In neither of these papers is the saturation profile continuously affected by permeability-saturation curves. Sheldon and Dougherty- recently described a method that employs continuously changing saturation profiles using permeability curves, has a minimum of assumptions and needs no prior sweep efficiency. The Higgins-Leighton method, described in this paper, has all of these features; however, many of the techniques are different from those of Sheldon and Dougherty. The Higgins-Leighton method, tested in May, 1961, is direct and easy to apply and requires very little computer time to calculate a forecast. The short computer time is especially helpful in the study of a reservoir containing many layers of different relative permeabilities. In the Higgins-Leighton method, the individual pressures do not have to be calculated, as the resistance to flow in each cell in the flow pattern is readily determined without the use of any iterative techniques. The saturation and permeability distributions are readily determined. These data and the shape factor, which is measured only once from the potentiometric model when mobility ratio is one, determine the resistance to flow in each cell. THEORY The authors showed in a previous paper4 that, as an aid to calculating performance, the reservoir can be divided into channels using the streamlines of a potentiometric model as a guide. See Fig. 1. This procedure also was used in this paper. The authors also showed that, by treating conduits as approximately one-dimensional and neglecting pressure gradients transverse to the main flow, the Buckley-Leverett equation may be expressed as The principles expressed by this equation are employed extensively in the three-phase flow, as they were in two-phase flow. In the three-phase flow, the channels (the size and shape of which are taken from a potentiometric model)