Search Documents

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)

By R. W. Nelson

Methods for determining permeability presented by W. D. Kruger3 in a paper entitled "Determining Areal Permeability Distribution by Calculation" have considerable merit. They can be expected to contribute significantly to analysis of flow in heterogeneous systems. Some rather subtle yet important conditions must be met, however, to assure that the computed permeabilities are correct. The conditions which insure a unique permeability distribution may have been overlooked in the method as proposed and presented in the subject paper. The subtlety of the special requirements may permit good agreement between computed and measured pressures in the reservoir; but the computed permeability may still be in error. The general theory and discussion of the requirements for a unique determination of the permeability distribution are presented elsewhere in more detail.7'8 Accordingly, only major features will be presented here in the notation of Kruger. GENERAL DESCRIPTION OF REQUIREMENTS A brief mention of the route to be followed in showing the special requirements may be helpful. The equation to be solved is a quasilinear, first-order, partial differential equation in the unknown permeability. This equation can be solved through an extension of Lagrange's method by reduction to a system of subsidiary ordinary differential equations. Through consideration of the identity of one of the Lagrange subsidiary differential equations with the stream function, a special interrelationship can be shown. The interrelationship has special significance with respect to the boundary condition. If the boundary function satisfies part of the subsidiary differential equation, no unique solution for the permeability distribution exists. In the physical problem this requirement for uniqueness indicates that the permeability measurements used as a boundary condition can not be along a stream tube. The boundary condition must be along a line, or series of lines, which completely traverse the region and pass through every stream tube.

By E. W. Hough, G. L. Stegemeier

Empirical equations for surface tension of propane and normal butane as functions of reduced temperature are obtained from experimental data. Another correlation relating surface tension to enthalpy of vaporization is given for these two compounds. In addition, new paracbor numbers are calculated for the normal paraffin hydrocarbons. These numbers are utilized for the calculation of interfacial tension of two-component systems as functions of pressure and temperature, using a modified form of Weinaug-Katz equation. The experimental data for two binary systems are approximated by the correlation. From these results it is found that the interfacial tension in the high-pressure region remains extremely low at large pressure decrements below the critical pressure. Thus, it appears that condensate systems may have flow characteristics almost like single-phase conditions even though the pressure is within the two-phase region. Experimental data have shown that interfacial tension divided by density difference approaches zero as the critical pressure is approached. A calculation of wetting-phase saturations indicates that the saturation gradient at the two-phase contact becomes increasingly abrupt as the critical pressure is approached. DISCUSSION Prediction of the surface and interfacial tension of the light hydrocarbons and of two-component hydrocarbon mixtures at various temperatures and pressures may be made if other physical properties are known. Extensive experimental work on single-component and binary systemsl is the basis for the correlations outlined in this paper. Interfacial tension is defined as the specific surface-free energy between two phases of unlike fractional composition, while surface tension is defined as the specific surface-free energy between two phase:; of the same fractional composition. The usual definitions relating interfacial tension to a liquid-liquid interface and surface tension to a gas-liquid interface are not clearly defined when the critical region is included, and there is no sharp distinction between a gas and a liquid phase. Interfacial tension is probably the most important single force that makes one-half to one-third of the total oil actually in place in a reservoir rock unrecoverable by conventional gas-drive or waterflood methods. A rough estimate of this figure for the United States is 100 billion bbl. Interfacial tension presently is used by petroleum engineers in the estimation of saturation gradients at the gas-oil contact and at the oil-water contact. The data in this paper should prove useful for estimates of reserves involving gas-oil contacts. Relative permeability undoubtedly is influenced by interfacial tension, for sufficiently small values. These data should be useful in determining how small the values are In addition, these data should eventually add to our fundamental knowledge of surfaces. At the critical point, all surface excesses approach zero and the thickness becomes very large. SINGLE-COMPONENT SYSTEMS It has been observed1 that the following relationships are good approximations to the physical properties of propane and n-butane.

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By J. W. Lacey, F. H. Brinkman, J. E. Faris

This paper presents an analysis of the high-pressure, gas-driven LPG-slug process, based on fluid flow tests in areal models. Two types of tests were made. One series was made in low-pressure models which permitted observation of fluid movement. Three completely misci-ble analog fluids were used. A second series of tests was made in high-pressure models using methane, propane and a light refined oil saturated with methane at room temperature and 1,550 psig. Under the test conditions of room temperature and a pressure level of 1,550 psig, the phase diagram for the fluids used is similar to those for many of the field systems where the process is considered for use. A method for using these laboratory data to calculate field performance of the process is outlined. As a result of this work, it is concluded that small banks of LPG (5 per cent HV or less) are not effective in increasing oil recovery in horizontal reservoirs. Znstead, where small banks are used, the driving gas quickly penetrates the LPG bank because of fingering and channeling; and from this point on, the process behaves essentially as an immiscible gas-injection project. The validity of this conclusion was substantiated by: (1) laboratory studies of the effect of rate, model size and mobility ratio on miscible displacement in areal models; and (2) calculation of field recovery, which compared closely with actual field recovery. INTRODUCTION Field applications and pilot tests of the gas-driven, LPG-bank, oil-recovery process are on the increase. Most of these tests are employing small banks (2% to 5 per cent hydrocarbon volume) of LPG, which is miscible with both the driving gas and the oil in place, in an effort to attain an effective yet economical miscible displacement of oil by gas. The expectation of miscible displacement with small banks is based on the concepts that: (1) lengths of solvent-oil mixed zones, measured during miscible displacements in long slim cores, are representative of those that will occur in the field; and (2) areal sweep efficiencies5 measured in electrolytic model studies are applicable to miscible displacement in reservoirs. Our experimental evidence indicates that the mixed-zone lengths and sweep efficiencies mentioned are not applicable to miscible displacement in reservoirs. In this paper we present an evaluation of the gas-driven, LPG-bank oil-recovery process based on fluid flow experiments in areal models. These results are used to predict the performance of a field pilot test of this process, and the results are compared with the actual test results. THE EFFECTIVENESS OF LPG BANKS The effectiveness of LPG banks of various sizes in accomplishing miscible displacement of oil by gas was determined by displacement tests in a model representing one-quarter of a confined five-spot pattern. The model, 12 X 12 X 1/4 in. in size, was packed uniformly with glass beads. It was operated at a pressure of 1,550 psig and at room temperature; the fluids used were methane, propane and refined oil saturated with methane at 1,550 psig. The saturated oil had a viscosity of 1.2 cp at room temperature. The results of these tests are shown in Fig. 1, where recovery is plotted as a function of the volume of fluid injected for: (1) an immiscible-gas drive with an oil-to-gas viscosity ratio of 85; (2) three sizes of LPG banks; and (3) two completely miscible displacements with mobility ratios of 85:1 and 10:1. (The miscible displacement data plotted are the averages of several tests.) The results show that a 2½ per cent bank of LPG does not increase oil recovery over that obtained by immiscible-gas drive. The 7 and 17 per cent banks are, respectively, about 30 and 50 per cent effective (with reference to the M = 85 miscible-displacement curve) in increasing oil recovery at the point where 2½ hydrocarbon volumes (HV) of fluid have been injected. The percentage effectiveness of the banks at a given volume of fluid injected is defined as Recovery by Bank — Recovery by Gas Drive Recovery by Miscible Recovery by Flooding Gas Drive We conclude from these results that banks of LPG smaller than about 5 per cent HV in an individual stratum will cause from little to no increase in oil recovery. This finding was substantiated by work in low-pressure models, which permitted visual observation of fluid movements. In these tests, three completely miscible fluids having the same viscosity ratios as those used in

By P. Hazebroek, M. Prats, E. E. Allen

Well locations in water floods frequently do not coincide with the ideal well locations associated with the flood pattern which is used. A question of importance is. "How much oil is lost because of this?". This is the problem treated in this paper. It is assumed that each pattern in a five-spot flood has the proper number of production and injection wells but that the wells may be off-pattern. The distance by which a well is off-patern is given in terms of the pattern spacing, defined as the length of one side of a five-spot. The influence of off-pattern wells on the performance of a five-spot flood has been investigated mathematically for unit mobility ratio. All production and injection rates are assumed to remain equal. The loss in recovery resulting from pattern irregularitie: has been calculated at a cumulative water injection corresponding to a 98 per cent water-cut in a regular five-spot. Displacing a single production or injection well by about one-third the pattern distance reduces the total field recovery by 4.4/N per cent of the recovery obtainable with regular spacing, where N is the number of five-spots in the flood (assumed large). When off-pattern wells are separated by at least one regular five-spot and are displaced by less than one-third the pattern distance, the total loss in recovery in a large-scale flood is, at most, 1 per cent of the recovery obtainable from that field with a regular five-spot pattern. However, much higher losses may occur if the off-pattern wells are not separated by a regular five-spot, if the rates are imbalanced, or if the mobility of the injected water is higher than that of the crude. Although the loss in oil recovery due to pattern irregularity often may be low, the off-pattern case will always require injection of more water to reach the same final water-cut. The volume of water injected per unit volume of oil displaced has been calculated for each five-spot adjacent to a five-spot containing a single 08-pattern well. The cumulative water injection at the time of 98 per cent water-cut may be increased by as much as 37 per cent when a production well is displaced one-third the pattern distance toward an injection well. The existence of differences between actual subsurface locations and the idealized well patterns used in secondary recovery operations is well recognized. However, the methods commonly used for predicting waterflood behavior ignore such differences, even when subsurface locations are known. The reason for this is simple. The idealized well patterns are highly symmetrical and, as such, the production history of the pattern (and even of the field) can be determined from that of a few production wells—usually from only one. Consideration of the many irregular well spacings found in a specific field would entail a large amount of work, and then the results would be applicable only to that specific field. In an attempt to provide an answer to the more general problem, this paper considers a water flood having an idealized well pattern, except for one well which is not at its regular location. When several off-pattern wells are present, their contributions are assumed to be additive unless the off-pattern wells are close to each other. Only four positions of the off-center well (at x, y. of Fig. 1) are considered in this paper. These are shown