Search Documents

By T. S. Hutchinson, C. E. Kemp

The bottom water drive analysis presented by Muskat has been extended to include fields with wider well spacings. Curves are presented from which volumetric sweepout, water-oil ratio, and rate relations may be calculated. The effects of non-unity ratios on sweepout eficiency and productivity decline have been studied by an approximate method. Also, the effect of gravity on volumetric sweepout has been approximated for a particurlar geometry. INTRODUCTION The purpose of this paper is to report work done in applying the analysis of bottom water drive reservoirs given by Muskat1 to a larger class of reservoirs, to extend this analysis to include the effects of gravity and non-unity mobility ratios in an approximate manner, and to investigate the rate of decline of productive capacity. The bottom water drive type reservoir is one where the oil is underlain by a water zone of sufficient thickness that water movement is essentially vertical. This system is defined in more detail in the appendix and a one well element is shown diagrammatically in Fig. 1. Muskat obtained solutions to the equations describing the volumetric displacement efficiency of idealized bottom water drive systems by assuming that gravity effects due to differences in density betwecn oil and water were negligible and that the mobilities of oil and water were the same. If either of these assumptions is not made, the complexity of the system makes an analytical solution well nigh impossible at the present stage of mathematical development. Coles and Stade' developed an approximation to consider these factors. The authors have refined approximations, starting with the same basic assumption, which should give a somewhat quantitative picture of bottom water drive performance when mobilities are not equal and the effects of gravity are not negligible. The results of these approximate calculations can then be used to obtain the decline in well productivity as the field is depleted.

Jan 1, 1957

By H. J. Morel-Seytoux

Methods of predicting the influence of pattern geometry and mobility ratio on water flooding recovery predictions are discussed. Two methods of calculation are used separately or concurrently. The analytical method yields exact solutions in a convenient form for a unit mobility ratio piston-like displacement. A few typical pressure distributions, sweep efficiencies and oil recoveries are presented for various patterns. For non-unit mobility ratio, one may resort to a numerical method, such as that of Sheldon and Douherty. 1,2 Because the domains of applicability of the analytical and numerical techniques overlap, the exact solutions provide estimates of the errors in the numerical procedures. The advantages of the analytical and numerical methods can be combined. To develop a numerical technique as independent of geometry as possible, the physical space is transformed into a standard rectangle. The entire effect of geometry is rendered through one term, the "scale-factor", derived from mapping relations. The scale factor can be calcu-lated from the exact unit-mobility ratio solution for the particular pattern of interest. By this means recovery performances for arbitrary mobility ratio con be obtained for many patterns. A sample of tesults obtained in this manner is presented. INTRODUCTION Pattern geometry and mobility ratio are two major factors in making a waterflood recovery prediction. Because assisted recovery has become increasingly important to the oil industry, pattern configuration and mobility ratio also assume a greater significance in the assessment of the economic value of recovery projects. The influence of pattern geometry and mobility ratio in shaping a recovery curve and on the other quantities of interest to the reservoir engineer is the main subject of this paper. Much effort has already been spent on estimating quantitatively the influence of either pattern or mobility ratio or both on oil recovery. The literature reports many investigations of this nature. 3-9 However, many results or methods of recovery prediction presented in the literature cannot be considered fully satisfactory. Even for unit mobility ratio and piston-like displacement, where analytical solutions are available, the literature shows discrepancies. For non-unit mobility ratio, the divergence in the results is extreme. For infinite mobility ratio in a repeated five-spot, depending on the investigator, the sweep efficiency ranges from 0 per cent to 60 per cent. With respect to the influence of pattern on recovery, only the repeated five-spot has received much attention. Other confined patterns and pilot configurations have received very little attention. Two calculation methods are presented in this paper, either separately or concurrently: the analytical method of potential theory and the numerical method of finite-difference approximation. The analytical method is more restricted in scope than the finite-difference method, but it has the definite advantage of providing exact solutions within its range of applicability. If a unit-mobility ratio piston-like displacement is assumed, the analytical approach is possible. A few typical results are reported in this paper; the detailed description of the general method and of a great variety of results will be the subject of other articles. Fornon-unity mobility ratio, we must resort to a numerical scheme. The numerical technique is that which was described by Sheldon and Dougherty.l,2 It is not limited to piston-like displacement. However, mainly single interface results will be presented here. Because the respective domains of applicability of the analytical and the numerical method overlap, useful comparisons of exact and numerical solutions can be made for a variety of patterns. The advantages of the analytical and numerical approaches can be combined. The reason for the success of this analytical-numerical approach can be summarized in the following two points: 1. The numerical solution for arbitrary mobility ratio can be programmed most efficiently when the physical space in which the displacement actually takes place is transformed into a standard shape; and 2. This can be done with remarkable simplicity whenever an analytical expression for the pressure

Jan 1, 1966

By C. S. Matthews, H. C. Lefkovits

Drilling progress is often delayed by sticking of the drill string. The development of preventive and remedial methods has been hampered by incomplete understanding of the sticking mechanism. A recent lahorntory investigation hns indicated that one type of sticking may be attributed to the difference in pressure between the borehole and formation. This paper shows, by means of soil mechanics, that the primary cause for differential pressure sticking is cessation of pipe movement, whereas diflerential pressre and stanrtding time determine the severity of the sticking. The analysis stresses the importance of using low-weight muds with low solids content and low water loss to alleviate diflerential pressure sticking and describes why packed hole drilling, long strings of drill collars, and a large deviation from the vertical are conducive to sticking. Finally, preventrve and remedial methods ore evaluated, and a theory is presented on the release of stuck pipe by spotting oil. INTRODUCTION Since drilling with long strings of oversize drill collars has become standard practice in many areas, the incidence and severity of the stuck pipe problem has increased. It has been noticed that in the majority of these cases the sticking could not possibly be attributed to key seating or caving of shales. It appeared that, due to the differential pressure between the mud column and the formation fluid, the collars were pressed into the wall and so became "wall stuck". Points to note about differential pressure sticking are: (1) sticking is restricted to the drill collars, (2) the collars become stuck opposite a permeable formation, (3) the sticking occurs after an interruption of pipe movement, (4) circulation, if interrupted, can be restarted after the sticking is noticed, and (5) no large amounts of cuttings are circulated out after restarting circulation. Helmick and Longleyl investigated pipe sticking by differential pressure in the laboratory and found an empirical relationship between the differential pressure, the sticking time and the required pull-out force. In this paper an explanation of the mechanism is given based on Terzaghi's theory of clay consolidation. A qualitative description is given in the following paragraphs while the derivation of fonnulas is given in Appendices. This paper is a first attempt to explain pressure differential sticking and many points will require additional theoretical and practical investigation before the problem can be fully understood. PRESSURE DIFFERENTIAL STICKING AS A CONSOLIDATION PROBLEM In any borehole, where the mud pressure is higher than that exerted by the formation fluids, a mud cake is formed opposite the permeable sections of the hole and a continuous flow of filtrate takes place from the mud, through the cake and into the formation. This radial flow pattern requires a certain distribution of the hydraulic and the effective (grain-to-grain) stresses inside the mud cake. Any quantitative or qualitative change in the external pressure conditions will produce a change in the flow pattern and, consequently, also in the internal stress distribution inside the cake. In view of the low permeability and the high compressibility of a clay mud cake, the adjustment of the internal stress distribution is slow and is accompanied by a change in volume. Time dependent stresses are thus created which gradually diminish as the new state of equilibrium between internal and external pressures is approached. Some 30 years ago, Terzaghi developed his "Theory of Consolidation" to account for the time-dependent stresses and settling of clay formations under the influence of external loads. He derived a differential equation by which the time-dependent hydraulic stress and the consolidation can be computed for any point inside the layer during the consolidation process. His theory is based on the assumption that the change in stress is solely due to a change in water content and it may only be applied to one-dimensional consolidation phenomena. Other investiga-tors5,10 have expanded his theory to include processes of more than one dimension. The difference between the external pressures on the mud cake before and after sticking is a qualitative one (isolation of part of the cake by the static contact with the drill collars after pipe movement has been stopped)', and the time-dependent stresses thus created may be investigated by means of Terzaghi's theory. By this analysis the changes in the nature of the contact surface between the drill collars and the mud cake during the sticking can be explained; and the friction force between the two may be computed as a function of the sticking time, the borehole dimensions and the mud cake characteristics.

By R. D. Taylor, Jim Douglas Jr., H. H. Rachford Jr., P. M. Dyke

A major obstacle to the use of wetting agents in .secondary recovery by water flooding is the adsorption of the agents on the sand. As a result of adsorption, the surfactant always lags behind the floodwater front. Consideration of the chromatographic theory of adsorption indicates that the detergents will not lag as much if used in very high concentrations. An investigation was made of the possibility of using high concentrations economically by flowing slugs of wetting agents followed by normal flood water. The experiments consisted of adsorption studies on Alundum powder and Berea sandstone. Flow rests on a 12-in. Alundum core and 22-in. Berea core were used to determine rate of detergent movement. The results of the flow experiments indicate that the relative rate of surfactant advance is, indeed, sensitive to the concentration of the agent. A 10 per cent slug moved with a rate that war 78 to 95 per cent as fast as the rate of advance of the flood water. By contrast, one with 25 ppm (the number of parts of commercial detergent in a million parts of water on a weight basis) concentration moved less than one-fourth as fast as the flood water, and calculations indicate that in very long porous systems the rate of movement of the lower concentrations will be a small fraction of the rate of advance of the flood front. The results. of the adsorption studies were utilized to calculate the rate of advance of the detergent when only the initial concentration was known. The calculated rates showed substantial agreement with the experimental flow tests in the high concentration ranges. The adsorption results were also used to estimate the cost of the materials for a slug-type surfactant flood in the field. In addition to the faster rates of movement, the concentrated detergent slugs removed much more oil than the dilute solutions. However, the effectiveness of the slug process depends on many variables. The quantity of oil removed is increased markedly by increasing the flooding rate. The efficiency is also influenced by the type of crude, type of reservoir rock and initial water saturation. Therefore, a careful analysis of each reservoir system is required before the economic value of the process can be determined. INTRODUCTION It is well known that the displacement of oil by invading water during water flooding is far from complete. It is generally agreed that the unrecovered oil is retained in the porous medium by the capillary forces which may be relatively large compared to the forces generated by the flowing water. Therefore, it was logical that some early workers should turn to surface-active materials to reduce the capillary forces to facilitate the release of oil. As early as 1927,' a patent was granted for the use of surface-active materials in water flooding. In 1932, when soap solutions were passed through Bradford and Venango sands, it was reported that the results were inconclusive, erratic and that "further investigation is needed to determine exactly the function of the solution and to obtain a clearer insight into the phenomena involved."' Some of the modern scientific reports conclude with a similar statement,' showing that the lack of agreement on the mechanism of oil removal by wetting agents is still very widespread even though several comprehensive studies have been reported.'." Although there is a lack of agreement as to the general effectiveness of the detergents for water flooding, most investigators do agree that all of the common detergents are strongly adsorbed onto the solid surfaces of the reservoir. In the early calculations it appeared that all additives would be lost before reaching much of the formation area which contained the additional oil to be removed. Experiments indicated that if the usual small waterflood concentrations of wetting agents were used, the rate of advance of detergent through the formation would be only a small fraction of the rate of advance of the flood front. Indeed, some investigators4 felt that the use of wetting agents would never be economically feasible because of their adsorption. For example, DunningG estimated that the wetting agent in concentrations of 25 ppm, would advance only 0.05 times as fast as the flood front. Ojeda, et al,' found that a surfactant in a concentration of 10 ppm moved less than 0.01 times as fast as the flood front. It is significant, however, that both investigators found that increased concentrations of wetting agents moved faster, relative to the flood front, than solutions at the lower concentrations. Ojeda showed that an extrapolation of his data indicated a relative rate of 0.5 at 1 per cent concentration, while Dunning6 estimated a relative rate of 0.46 for a 1 per cent concentration. It was obvious that these concentrations could not be used for continuous injection because the cost of the injected detergent would far exceed the value of additional oil produced. Traditionally, detergents are used in very low concentrations for they show good

By W. W. Yankie, A. T. Chatas

This paper is written as a discussion of the paper, "Buckling of Tubing in Pumping Wells, Its Effects and Means for Controlling It" by Arthur Lubinski and K. A. Blenkarn, which war published in the March, 1957, issue of Journal oP Petroleum Technology. An analogy is drawn between the buckling of tubing in pumping wells and a fiube, fitted with frictionless pistom, linked by an inelastic rod, and which is subject to an internal pressure. A simple derivation shows that the critical internal pressure exerts a hydrostatic force on the pisfon which is close to the Euler load for compression buckling of the same tube. While the results agree with those of Ref. I, the basis of calculations has been changed and clarified. The presence of internal upsets at the pistons is discussed. INTRODUCTION In the Appendix, Lubinski and Blenkarn state that a tube with ends closed by frietionless pistons connected by an inelastic rod (Fig. 1) will behave as a column loaded with a force equal to the tension force in the rod. They note that, while the fluid may not resist a shear force, the moment in every cross section is piston force multiplied by the elastic line deflection (c.f., the loaded column). Since the equations of the elastic lines are similar in the two cases, the solution is carried over from the loaded column. On the contrary it is suggested that while the method of solution remains unchanged from the column case, the actual stress distribution is different, a fact which is taken into account in this note. A simple method of approximating the buckling pressure is given, and their reasoning in our mathematical language is continued. GENERAL EQUATIONS The equation of an elastic line, initially straight in the x direction (Fig. 2), is given by In the case of a loaded column, this may be written: For the pressurized tube it was shown, by analogy with the loaded column, that tion has a solution that, for P, = EI/P, any deflection is stable (a criterion for buckling). An attempt will be made to describe more precisely, the forces acting on the tube, and to analyze the buckling condi,tions. DISCUSSION If the tube assumed a deflected position, the volume between the pistons would increase so that in deflecting, the energy of the fluid would decrease. When the increase in elastic energy of the tube (due to bending) is equal to the decrease of fluid pressure energy, buckling conditions exist. Analyzing the forces on the tube (Fig. 3), two transverse forces, T,, T,, exerted by the pistons, are balanced by a resultant of the pressure elements which are linearly dependent on the local curvature of the elastic line (Fig. 4). Unlike most buckling cases, the forces depend on the shape of the elastic line, a condition that makes a rigorous mathematical treatment difficult. Several approximations may be used to clarify the problem.

By H. Hooykaas

A few years ago Marsal published a method of predicting oil recovery as a function of time for an edge-water-drive reservoir with several rows of wells.' The method is based solely on oil- and water-production data obtained during a short period for the first row of wells and the gross productions (past and future) of all wells. This paper presents a streamlined version of the original publication and experimental evidence showing the reliability of this version. For this purpose scaled experirnents were performed in a large model (750 X 500 X 23 cm, 24.4 ft X 16.2 ft X 8.9 in.) in which edge-water-drive reservoirs with three rows of production wells were simulated. From the results of the experiments it is concluded that the method gives a fairly accurate prediction of the future production behavior of the reservoir and enables us to select the most profitable production plan. Results obtained with the method were always slightly conservative. INTRODUCTION The calculation method developed by Marsal aims at forecasting future reservoir behavior.' It can be applied to floods in which a number of down-dip wells are rapidly going to water. Essentially it is a method for predicting the future behavior of the up-structure wells on the basis of the observed behavior of the down-structure wells. As such, the method is especially useful for predicting behavior of edge-water drives toward several rows of wells. The requirement to be met for justifying the approach is identity of production mechanism in all regions of the reservoir. The attraction of the method then lies in the fact that it applies even when we have no complete description of the prevailing production mechanism. The original publication was in the German language and for this reason, perhaps, did not attract widespread attention. In our opinion it deserves to do so because it does offer in a simple way results which are in fair agreement with observed data. The aim of the present paper is to bring the method to the attention of a wider audience. The method in its most simplified form is reviewed and experimental data in support of it are presented. PRINCIPLES To show Marsal's principle. we will first write the continuity equation for very simple conditions, linear flow (in the positive x-direction) in an incompressible porous body [cross-section A(x), porosity +(x)], saturated with incompressible fluids (oil saturation S). When the oil in this formation is displaced by water, injected at x = 0 (inlet face), we get where q(t) is the total gross displacement rate and f the fractional oil flow. The equation can be solved if we make the assumptions: (1 ) that the fractional oil flow (f) is a univalent function of saturation (S); and (2) that the function is the same throughout the porous medium. Then: f = f(S) Solving this equation with the method of characteristics, we obtain

Jan 1, 1966

By E. J. Bonet, P. B. Crawford

It is fairly common practice to reinject water into the aquifer near the oil-water interface in water-drive reservoirs. There have been many studies of aquifer behavior without reinjection,l-5 but apparently no analytical studies have been made to indicate performance when injection wells are present in the aquifer. The diffusivity equation may be used to describe the pressure behavior in an aquifer. When sources are present it can be written as A general solution of Eq. 1 may be obtained as the sum of two solutions. One is a general solution of the homogeneous equation, obtained by dropping the source term in Eq. 1 and still satisfying the required boundary conditions. Another is a general solution of Eq. 1 satisfying homogeneous boundary conditions.1° The first solution, which corresponds to the aquifer without reinjection, has been given quantitative and qualitative treatment in the literature.l-5,8,9 The latter solution will be treated here. Solutions with variable injection rates are obtained by superimposing solutions with constant rates, so we have considered only the latter. Five cases have been treated here: 1. A linear aquifer of finite width and extent. Water is injected into the aquifer at point xo, yo at constant rate. The pressure along the line x = 0 is maintained at zero pressure. (See Fig. 1.) Equations showing the pressure distribution, rates, velocity and cumulative influx have been developed. Certain numerical results are shown. 2. A linear aquifer of finite width and extent. Water is injected into the aquifer at point xo, yo at constant rate. The rate of production along the oil-water water contact is maintained constant. Equations showing the pressure distribution inside the aquifer are presented. 3. A linear aquifer of finite width and extent. Constant pressures are maintained at opposite ends of the aquifer. Water is injected at constant rates at point xo, yo. Equations showing the pressure distribution are shown. 4. A linear aquifer of finite width and infinite extent. The pressure is maintained constant at the oil-water contact. Water is injected at a constant rate at point xo, YO. 5. A linear aquifer of finite width and finite extent. The pressure at the oil-water contact is maintained constant. Water is injected at constant rate per unit length along a line x,. The Linear Aquifer The linear closed aquifer, performing at constant terminal pressure was selected as the first case to be analyzed. The development of the equations is based on uniform injection per unit of thickness, and gravity effects are neglected. In practical application this condition is probably most closely achieved when the

Jan 1, 1970

By J. S. Lasater

Resu1ts of experinmental measurernents of heat capacities and thermal conductivities of some typical porous rocks are presented. Measured heat capacities agree closely with va1ues calculated front known chenlical . compositions of the rocks. On the basis of this agreement, heat capacitir.7 of fluid-satttrated rocks were calc~ililted. Therrnal conductivities of the rock samples were measured under varioits conditions of fluid saturation. From these data therrnal dif-fvsivities were calculated. A significant variation of thermal diflusivity with trnlperature is indicated. Knowledge of the thermal characteristics of fluid-bearing porous rocks has become increasingly important with the development of thermal oil-recovery processes. Reliable thermal data for petroleum reservoir rocks are not available and approximate values are generally assumed for reservoir calculations. 'The effects of temperature, pressure and fluid saturation on thermal properties of porous rocks have not been carefully studied. The purpose of the present work was to measure and report the thermal characteristics of some typical porous rocks. A group of eight sedimentary rock samples, including limestone, shale, sandstone and silt-stone, was selected for the study. Properties reported include heat capacities, thermal conductivities and thermal diffusivities. The effect of fluid saturation on these properties was also considered. Determinations of thermal properties require precise measurements which are time consuming and difficult to reduce to a routine basis. Methods have not been developed for measurement of thermal proper-ties under conditions of pressure, temperature, and fluid saturation, which may be encountered in thermal recovery processes. As one approach to these problems, some preliminary work on methods of estimating thermal properties from other more readily measurable characteristics of the rock-fluid system is reported. The samples selected for the studies were oilfield cores from within, or closely associated with oil-productive zones. Descriptions of individual samples are given in Table 1. Clay minerals present in the samples were determined by differential thermal analysis (DTA). Quartz content was estimated from DTA cooling curves. For all but the thermal conductivity tests, a representative 200-gm sample of each core was disaggregated, thoroughly extracted in diox-ane and oven dried at 100°C. The cleaned and dried sample was care- fully split into appropriate portions for the several tests. Chemical analyses of the samples were obtained to permit calculations of heat capacities by methods which are outlined later. Resuits of chenlical analyses are given in Table 2. Heat contents of the test samples were measured by a Runsen-type calorimeter. The heat content of the unknown sample is measured relative to the known heat content of platinum. The apparatus is referred to as a relative error type because of the comparative measurements, and has the advantage of canceling systematic errors. Experimental results are reproducible within 0.5 per cent. A temperature base of 298.16 K is taken for heat content measurements. Measurements were made at five temperatures (260, 440, 620, 800 and 980° F). It was necessary to preheal the sample above the maximum experimental temperature to eliminate volatile constituents. Assuming the weight loss on pre-heat, ing was mostly water, a correction to the observed heat content was made to account for the water originally present in the sample. This latter procedure was followed in the determination of calculated heat contents as will be described later. The apparatus and the method of measure-

By T. D. Mueller, H. Dykstra

This paper describes a method for obtaining compositions of gas and oil phases in equilibrium with each other at a given reservoir temperature and pressure as lean or enriched gas is injected into a reservoir oil not in equilibrium with the gas. The method is based on calculation of equilibrium constants, K, from the equation log Kp = A + B F, where p is pressure and F is component characterization factor. A and B are constants that vary with a composition parameter. The paper explains how A and B could be determined with a minimum amount of data. The method consists of calculating equilibrium constants as a function of liquid composition and using the K's to calculate changes in compositiort of injection gas as it successively contacts original resetvoir oil at the gas-oil front and the change in composition of reservoir oil as it is successively contacted by the injection gas in the vicinity of an injection well. An iterative process is required. Two examples of the method are given which show how oil and gas compositions and properties are calculated with the aid of a computer. INTRODUCTION The calculation of recovery from lean-gas or enriched-gas injection into a saturated or under-saturated reservoir oil requires a knowledge of how phase properties change in a reservoir between injection well and producing well. To determine the phase compositions or properties experimentally would be difficult and time consuming, especially if one wanted to look at the effect of various injection gas compositions. For some fields an engineer would want to make preliminary calculations on the economics of gas injection before he makes the decision to spend a corrsiderable sum of money for laboratory work to obtain the necessary fluid properties. A method for calculatine compositions and properties based on a minimum of laboratory data is described. The minimum amount of data required is injection gas and reservoir oil composition and injection gas solubility as a function of pressure for an undersaturated oil, or the compositions of mixtures of reservoir fluid with two or more enriched gases, which mixtures have a bubble-point pressure that is the same as the expected average displacement pressure. On the other hand, if the economics appear to be attractive, more laboratory data can be obtained in the form of equilibrium constants for three or more mixtures of the reservoir fluid and injection gases. The method is based on the calculation of K's from the equation, log Kp = A + B F, where K is equilibrium constant, p is pressure and F is component characterization factor. This method of plotting K data by means of the above equation is described by Hoffman, Crump and Hocott.1 A and B are constants that vary with a composition parameter. It is shown how experimental data are used to determine how A and B vary with composition. BASIS OF METHOD EQUILIBRIUM CONSTANTS Equilibrium constant (K) data are required for this method. The K of each component is obtained from an equation that gives the product of K and pressure as a function of component characterization factor. The characterization factor F = b(l/Tb — l/T), where b is a constant characteristic of a given component, Tb is the boiling point in OR at atmospheric pressure of the given component, and T is the temperature of the system in OR. This method of plotting K data is described by Hoffman, Crump and Hocott.l The straight-line equation of the log Kp - F plot is log hp = log Kop +BFi, where Kop is at F =0, i designates the component, and B is the slope of the line. The factor F can be

Jan 1, 1966

By J. E. Kirby, L. B. Schnitz, H. E. Stamm III

The production history of a gas-cap-drive reservoir was reproduced by calculations, and predictions were made for operations under primary depletion, pressure maintenance by gas injection, and pressure maintenance by water injection. Although the general equations used in this analysis have been available to the industry for a number of years, the procedure has not had extensive use because of the unusually large time requirements for developing numerical results. The use of digital computing equipment for processing the calculation work generally reduced the time usually required for an extensive analysis of this type. Calculation techniques were developed that were equally applicable for solving the gar-cap-drive problem, whether digital computing equipment or desk calculators are used. These techniques are outlined and the procedures are illustrated by the field problem. INTRODUCTION The analysis of a gas-cap-drive res-ervoir requires the simultaneous solution of the volumetric balance1 and displacement2 equations. Equations for the analysis of combination drive reservoirs were presented by Wooddy and Moscrip in a previous paper." These equations permit the evaluation of possible operating processes in order that the optimum method of pro- ducing the reservoir may be determined. Although these equations represent the basic technique known for the detailed analysis of gas-capdrive reservoirs, the calculation processes are quite complex and time consuming; hence, the procedure has not had extensive use. Recently, these equations were satisfactorily applied with the aid of a digital computer to investigate the completion history of a gas-cap-drive reservoir under the following operating methods: primary depletion, depletion under full pressure maintenance by gas injection, and depletion under full pressure maintenance by water injection. The use of digital computing equipment substantially reduced the time necessary for making the computations and permitted a more concise investigation of the reservoir's depletion history than was possible when the desk calculator was used. This more rigorous investigation also permitted a detailed analysis of the calculation procedures necessary in the solution of the gas-cap-drive problem. This paper outlines the procedure of analyzing a gas-cap-drive reservoir and illustrates the process with an actual field problem. CHARACTERISTICS OF THE RESERVOIR ANALYZED Production in the field under consideration is obtained from a Frio formation of Oligocene age at a depth of approximately 5,740 ft sub-sea. The accumulation, which is bounded on the flanks by strand lines, is semirectangular in shape with gentle narrowing toward the base of the structure as shown by the sand limits and completion map in Fig. 1. The reservoir development indicates an oil column of 319 ft that is defined by gas-oil and oil-water contacts. Basic data relative to the characteristics of the formation and the formation fluids are shown in Table 1. The formation exhibits considerable shaliaess and thinning near the oil-water contact indicating that little or no benefit will be received from natural water influx. This is substantiated by the fact that wells completed near the original oil-water contact have not produced water. A sizeable gas cap is present; however, considerable faulting is in evidence, and it is believed that only a limited portion of the gas cap will aid in oil displacement and pressure maintenance. Although the reservoir was in an initial stage of depletion, it was believed that as soon as the equilibrium gas saturation in the oil zone was exceeded and a number of upstructure wells were invaded by the expanding gas cap, the resulting high gas-oil ratio production would rapidly dissipate the natural drive and lead to limited oil recoveries and inefficient reservoir operation. This analysis evaluated the effectiveness of the primary displacement mechanism and

Jan 1, 1958

By N. R. Morrow, C. C. Harris

The experimental points which describe capillary pressure curves are determined at apparent equilibria which are observed after bydrodynamic flow has ceased. For most systems, the time required to obtain equalization of pressure throughout the discontinuous part of a phase is prohibitive. To permit experimental points to be described as equilibria, a model of capillary behavior is proposed where mass transfer is restricted to bulk fluid flow. Model capillary pressure curves follow if the path described by such points is independent of the rate at which the saturation was changed to attain a capillary pressure point. Aemodified suction potential technique is used to study cyclic relationships between capillary pressure and moisture content for a porous mass. The time taken to complete an experiment was greatly reduced by using small samples. INTRODUCTION Capillary retention of liquid by porous materials has been investigated in the fields of hydrology,l soil science,2 oil reservoir engineering,3 chemical engineering," soil mechanics,5 textiles: paper making: and building materials In studies of the immiscible displacement of one fluid by another within a porous bed, drainage columnsl and suction potential techniques2 have been used to obtain relationships between pressure deficiency and saturation (Fig. 1). Except where there is no hysteresis of contact angle and the solid is of simple geometry, such as a tube of uniform cross section, there is hysteresis in the relationship between capillary pressure and saturation.2 The relationship which has received most attention is displacement of fluid from an initially saturated bed (Fig. 1, Curve R2), the final condition being an irreducible minimum fluid saturation Swr. Imbibition (Fig. 1, Curve A), further desaturation (Fig. 1, Curve R), and intermediate scanning curves have been studied to a lesser but increasing extent.2 ,3,9-12 This paper first considers the nature of the experimental points tracing the capillary pressure curves with respect to the modes and rates of mass transfer which are operative during the course of measurement. There are clear indications that the experimental points which describe these curves are obtained at apparent equilibria which are observed when viscous fluid flow has ceased; and any further changes in the fluid distribution are the result of much slower mass transfer processes, such as diffusion. Unless stated otherwise, this discussion applies to a stable packing of equal, smooth, hydrophilic spheres supported by a suction plate5 with water as the wetting phase and air as the nonwetting phase. The discussion can also be applied, in many

Jan 1, 1966

By H. H. Rachford, D. S. Howard

In the withdrawal of fluids from tilted aquifers it is of value to be able to predict pressure patterns during the course of the pressure decline. As an example of this, in the displacement of fluids from oil reservoirs being produced by water-drive, prediction of the reservoir pressure is an important factor in fixing production rate. The prediction of pressure patterns in the depletion of an inclined aquifer has conventionally been made by solving the well known flow equations for the pressure Patterns in a horizontal aquifer, and correcting these patterm for the pressure difference due to the hydrostatic head associated with the vertical distance between a point in the horizontal aquifer and the corresponding point in the tilted aquifer. Inasmuch as the equation describing flow in tilted strata is different from that for horizontal flow. there was question as to the validity of the simple corrective procedure current1y in use. To examine the magnitude of error made using this procedure, an inclined aquifer was chosen for which the error was expected to he large, and the equation for flow in this tilted stratum was solved numerically. The results of the conventional procedure of correcting the solution of the equation for horizontal flow was compared with the solution of the correct equation by studying the difference between predicted pressures by the two methods. Three forms of the hydrostatic head correction were studied. The results for all three forms of the head correction are similar. The error in the approximate solution is quite dependent on the position of the horizontal aquifer with respect to the inclined one. The error is greatest at the producing end ad increases with amount of fluid withdrawn. The magnitude of the errors is not large if a good choice is made for the position of the horizontal datum. The most satisfactory position for the datum in all cases is near or slightly about the mid-point of the inclined stratum. Use of the best approximate method studied with the best datum tested yielded predicted pressures in error no greater than 7 psi. INTRODUCTION The movement of water by expansion from an aquifer to displace fluids from a contiguous oil reservoir was first recognized over 20 years ago.9 Since that time much work has been done to relate this mechanism to reservoir behavior mathematically. The result of this effort on the part of many workers has led to an understanding of underground fluid movement which has enabled us to follow and predict the behavior of the water-drive reservoir. One reason that good progress has been made in the application of mathematics to the study of pressure patterns and water movements in aquifers has been that the differential equation which describes this process is of a relatively simple form. The equation for the isothermal, horizontal, viscous flow of water through homogeneous sandstone is identical with that for the diffusion of heat through homogeneous solids, and this equation and means of obtaining solutions for it with many types of boundary conditions have been studied extensively. That this work has been highly successful is clearly emphasized by the large number of solutions developed for problems of interest both by methods of formal mathemat-

Jan 1, 1957

By Albert S. Trube

Increasing emphasis is being placed on the necessity for obtaining reasonably accurate estimates of the physical properties of reservoir fluids well in advance of more accurate laboratory data. One such factor is the isothermal coefficient of expansion of an undersaturated hydrocarbon liquid which may be contained in a particular reservoir. This coefficient, or liquid "compressibility", has often been assumed to be relatively insensitive and nearly constant. Although this assumption may be nearly correct in the case of high specific gravity liquids, it cannot apply in the case of medium to low specific gravity liquids. Any treatment of the nature of liquid compressibilities must give consideration to the variable nature of the isothermal expansion coefficient and the fact that it can be both pressure sensitive and temperature sensitive. The isothermal coefficient of expansion, or compressibility, of a substance is defined as -c = 1/V (?V/?p)T......(1) The value of Eq. 1 may be approximated for a finite pressure change by using the well-known relation, c = V1 - V1/Ve (p2-p1), ........(2) in which p1 > p1 and V2 > V2 and V1 + V2/2. The instantaneous compressibility may also be estimated from the slope of a curve on which log V is plotted vs pa. In this instance, it has been shown that c = In V1/V2/p2-p1 ...........(3) Eq. 3 may be used for variable compressibilities by plot-ting V vs p on semi-log paper. In this case, the instan-taneous or point compressibility is c = 2.303 m.......... . (4) in which m is the slope of the curve at a point in log cycles per psi. ISOTHERMAL COMPRESSIBILITY OF UNDERSATURATED LIQUID In the case of gases, the compressibility has been defined by Muskat' and others as Co = 1/p - ?z/z?p ..........(5) This was later generalized through application of the theorem of corresponding states' to co = 1/pr - (?z/z?pr) T8 ........(6) in which c, is defined as the pseudo-reduced compres-sibility and is numerically equal to (c.p,) and the sign is arbitrarily changed for mathematical convenience. Eq. 6 is general and applies to gases and liquids. In the case of liquids, c. = c1,pe, but if Eq. 6 is used, values of z and ?z/?p, must be determined for the liquid in question. The difficulty of obtaining accurate values of ?z/?p,, particularly in the region where T, is approaching 1.00, was pointed out in the previous paper'. This method of approach is unreliable for extension into computations of liquid compressibilities. On the other hand, c, can also be defined as follows: - Cr = 1/Vr (?Vr/?pr) Tr ......(7) and from which, Cr = 1/pr (?pr/?pr) Tr .....(8) In accordance with the theorem of corresponding states, Vr = V/Vr, p, = p/pr and pr = p/pr. Ref. 2 pro-vides an excellent source of reduced density, reduced pressure and reduced temperature data for a wide variety of fluids. For the purpose of this investigation, the set of fluids2 having ze = 0.27 was chosen because it in-cludes most of the heavier hydrocarbons usually found in those complex mixtures loosely defined as "reservoir oils". Data on undersaturated liquids in Ref. 2 were rearranged to conform with Eq. 8. To simplify the computation of cr, the isothermal p, data were plotted vs the log of pr, in which case for any straight line or tangent cr = 1/M pr pr .....(9)

Jan 1, 1958

By F. F. Skiba, P. J. Root

A mathematical model has been used to investigate crossflow effects in a stratified reservoir during an idealized displacement process. A process in which one incompressible fluid displaces completely another incompressible fluid of the same density and viscosity was considered. The reservoir was assumed to be composed of homogeneous layers which are discrete but interconnected. Approximate five-spot flow geometry was simulated by arranging two pie-shaped cylindrical (two-dimensional) wedges in series. Some of the common production methods rely, for their effectiveness, on the fact that the formation consists of isolated strata. This is particularly true of selective plugging and single-zone production-injection methods. The effectiveness of these techniques for several cases in which adjacent strata are in communication has been evaluated. The usual production-injection schemes for fluid injection processes involve using the same completion interval in both the production and the injection wells. The efficiency of this, as well as alternative productioz-injection procedures has been examined for some cases involving communicating strata. INTRODUCTION The performance of many fluid injection projects is marred by early breakthrough and by production of a high percentage of the displacing phase after breakthrough. This channeling effect, if it occurs in reservoirs which have large permeability variations, is often explained on the basis of a model composed of discrete, isolated strata. Corrective or preventive measures based on this same model are also common. The behavior of stratified models has been the subject of considerable study. An excellent review of the early literature was given by Seba.1 Studies of the depletion performance have since been made for continuous, unrestricted communication between strata2-4 and for production from one permeable bed in continuous but restricted communication with an adjacent permeable bed.5,6 Model studies of displacement processes in stratified systems with crossflow also have been reprted.7-12 Finally, methods which account for crossflow have been developed for computing displacement behavior in stratified systems.13,14 Most methods for controlling or preventing channeling involve some change in the interval open to production and/or injection. Although many variations of the selective completion problem have been studied in regard to the depletion of a reservoir, little attention has been given to the behavior of a displacement process under similar conditions. Sebal used a resistance network model to study the performance of a process in which one incompressible fluid displaces an identical fluid from a linear, stratified reservoir. He investigated the possibility of improving the vertical coverage by limiting the completion interval of the production well. This investigation deals with techniques in which selective completion of the production and/or the injection wells is used in an attempt to control channeling. At one extreme, selective plugging and single-zone isolation (assuming perfect execution) should be successful if the reservoir is actually made up of independent strata. At the other extreme, they can hardly be expected to be successful if the channeling is caused by a mechanism in which stratification effects are negligible even though there are extreme variations in the formation permeability and/or the injectivity profile of the wells is not uniform. By contrast, it is difficult to anticipate the effect that the amount of communication between the beds of a stratified system will have on the efficiency of selective completion techniques. The objective of this study is the quantitative determination of the effect of communication on the behavior of some particular stratified systems which will at least lead to qualitative conclusions.

Jan 1, 1966

By J. J. Arps

This paper reviews the methods currently in use for estimating primary oil reserves and discusses the principles on which these methods are based. Particular emphasis is placed on how these methods change with the type of information available during the life cycle of an oil property. This paper contains various novel estimating methods and shortcuts heretofore unpublished. INTRODUCTION Estimating oil reserves is one of the most important phases of the work of a petroleum engineer since the solutions to the problems he deals with usually depend on a comparison of the estimated cost in terms of dollars, with the anticipated result in terms of barrels of oil. His recommendations to management regarding the best course of action are therefore normally based on the most favorable balance between these two. Specific engineering problems which require such a knowledge of recoverable oil reserves and a projection of future rates are: [a.] the exploitation and development of an oil reservoir; [b.] the construction of gasoline plants, pipelines and refineries; [c.] the division of ownership in unitized projects; [d.] the price to be paid in case of a sale or purchase of an oil property, and the magnitude of the loan which it will support; [e.] the proper depreciation rate for the investment in oil properties; and [f.] evaluation of the results of an exploration program. This discussion will be confined to the various methods and tools which are currently in use for estimating oil reserves to be obtained during the primary phase of an oil-producing reservoir and for a projection of the future production rates. Reserves which may be obtained by secondary recovery methods or fluid injection programs and gas and gas condensate reserves will not be discussed in this paper. Unfortunately, reliable oil reserve figures are most urgently needed during the early stages when only a minimum amount of information is available. Management's interest in the oil recovery from a property—aside from its use for accounting purposes—usually declines when the property approaches its economic limit, just at the time when the reliability of the estimates is at its best. To illustrate this general idea a chart is presented (Fig. 1) showing the three periods in the life of an imaginary oil property. Time is shown on the horizontal axis, while the cumulative production and estimated ultimate recovery are plotted vertically. No particular units are used and this schematical chart is not to scale. During the first period, before any wells are drilled on a property, any estimates will of necessity be of a very general nature, based on experience from similar pools or wells in the same area, and usually expressed in barrels per acre. This will

Jan 1, 1957

By R. D. West

Miscible displacenzent recovers all oil in the area contacted by the injected .fluid, whereas water or immiscible gas drives usually leave substantial amounts of oil as residual. However, the Door mobility ratios associated with a gas-driven miscible displacement cause the sweep pattern efficiency to be much lower than that obtained with water flooding. One way in which the sweep eficiency in a miscible displacement process can be increased is by decreasing the mobility behind the flooding front. This can be achieved by injecting water along with the gas which drives the miscible slug. This water reduces the relative permeability to gas in this area and thus lowers the total mobility. The main operating conditions for the simultaneous injection -vrocess are that a zone of gas exists between the miscible slug and the leading edge of the water and that a su,@cient amount of gas be injected with the water to form the pas volume which is being left in the water zone. Laboratory model studies have shown that the ultimate sweep pattern efficiency can be as high as 90 per cent for a five-spot flooding system. If gar alone is used as the driving medium an ultimate sweep-out efficiency of about 60 per cent would be obtained in the same system. I INTRODUCTION The miscible displacement processes are a step towards total oil recovery. Conventional gas or water drives usually leave 25 to 5.0 per cent of the oil as residual in the swept portion of the reservoir. This residual can be eliminated if the oil is driven by a fluid with which it is miscible. At some reservoir conditions natural gas will become miscible with the oil. This is the "high pressure gas process".' More often, the oil does not contain enough light hydrocarbons to cause the gas to become miscible with the oil at reasonable pressures. In these cases a small band of fluid which is miscible both with the oil and gas must be kept between them2. Less than 2 per cent of the reservoir volume of the slug material is needed to keep the displacement miscible. Both processes work in the same manner, recovering all of the oil in the portion of the reservoir contacted by the injected fluids. The only difference is the manner in which the miscibility between the oil and the injected gas is obtained. Previous publications have contained detailed descriptions of these processes.1,2,3,4 However, total displacement of the oil in the swept region does not guarantee an efficient recovery process. The amount of oil to be recovered is also determined by the fraction of the reservoir contacted by the flood. This fraction is largely determined by the mobilities of the fluids. (The fluid mobility is the permeability of the rock to that fluid divided by the fluid's viscosity, k/p). This. dependence of the fraction swept on the mobility ratio has been shown in previous studies. Fig. 1 shows the ultimate fraction swept in a five-spot system as a function of the mobility ratio. The small drawings show the location of the areas left unswept for two different mobility ratios. The ultimate fraction of the reservoir swept is here considered to be attained when the producing stream contains less than 5 per cent oil at reservoir conditions. THE GAS-DRIVEN MISCIBLE DISPLACEMENT Since there is no oil left in the swept region after miscible displacement the mobility in this region is very high. It is often 50 times the mobility in the unswept regions. This means that the fraction of the reservoir contacted by the injected fluid will be less for a gas-driven miscible displacement than for a conventional water or gas drive. For a five-spot injection system, water would contact the entire reservoir volume, and the low pressure gas would contact about 90 per cent of this volume, while a gas-driven miscible displacement would only contact about 65 per cent of the reservoir. This poor sweep efficiency often offsets the benefits obtained through miscible displacement. Fig. 2 shows what the recovery curves for the three processes might look like for a five-spot system. The curves show the fraction of the in-place oil recovered as a function of reser-

By C. F. Gates, H. J. Ramey

Recent literature shows that pronounced increases in oil recovery can result from the use of miscible systems in recovery operations. This literature also points out certain problems associated with maintaining miscibility, e.g., compositional changes resulting from alterations in pressure and temperature or from zone dilution. The purpose of the work described in this paper was to study miscible zones for possible use in waterflooding operations and examine the manner in which these zones break down. The fluid system selected for study was oil—terNary-butyl alcohol—water, which, because of its narrow range of miscibility, is ideally suited for investigation in short, linear systems. The laboratory-observed effects on oil recovery and zone stability of a number of variables are reported. Among the variables investigated were miscible zone size, viscosity ratio, length of travel, and interstitial water saturation. Results obtained with this system show that interstitial water adversely affects recoveries because of premature phase break. Salinity of the water further aggravates Ais situation. This effect, although pronounced for this system, would probably occur in any fluid system containing an alcohol. As expected, viscosity ratio has a decided influence on recovery efficiency and zone stability. The results also show that use of a suficieno zone size, even though the system has poor phase characteristics, yields higher oil recoveries than are obtained from straight water floods. Length-of-travel studies showed a square root relationship between zone growth and path length for favorable viscosity ratios. No such clear-cut dependence was observed under unfavorable viscosity conditions. INTRODUCTION As new oil reserves have become increasingly more difficult and expensive to find and develop, the oil industry has devoted more and more time and money to finding more efficient methods for exploiting known reserves. The increasing number of water floods is one manifestation of this attempt to improve oil recoveries from existing fields. Since considerable quantities of oil are by-passed by the waterflood process, it is only a partial answer to the problem. The search continues. Of late, several proposed processes have received the critical attention of investigators. Among them are gas drives, in situ thermal reactions and miscible displacements. The present paper is concerned with miscible displacements in conjunction with water floods. Recent literature1- reports that, under proper conditions, oil recoveries from LPG floods, propane sweeps, condensing gas drives, etc., approach 100 per cent of the oil in place in the region contacted by the flood. This literature also points out some of the problems associated with these processes, one of the main ones being the maintaining of miscibility. Compositional changes resulting from pressure and temperature variations or from zone dilution give rise to phase breaks and loss of miscibility, the conclusion being that, after loss of miscibility, much of the beneficial influence of the zone is lost.' The results of the subject work show that, if sufficient zone material exists ahead of a phase interface, considerable improvement in oil recovery can be obtained. In order to increase our understanding of the behavior of miscible systems as they degenerate into immiscible processes, a laboratory system having rapid deterioration characteristics was needed. The fluid system oil—ter tiary-butyl alcohol (TBA)—water possesses these properties and was the system investigated. Although these fluids are different from those generally employed in miscible studies, an understanding of their behavior can be applied to other miscible slug processes. However, under conditions where mass transfer maintains a miscible system through multiple contacts, e.g., enriched gas drive or high pressure gas drive, the observed behavior of alcohol zones would not be applicable. The variables of interest were zone size, viscosity ratio, length of travel, initial water saturation, rate and core characteristics. Of these, only the effects of the following on zone breakdown and oil recovery are reported in detail: (1) zone size, (2) path length, and (3) initial water saturation. The effects of viscosity ratio and core characteristics are discussed in conjunction with other variables. Although initial water saturation is not envisioned as important to the phase behavior of hydrocarbon systems, e.g., LPG, it is important in alcohol systems. The properties of the fluids and cores used are given in Table 1. As previously stated, this fluid system was selected because of the limited range of miscibility

By H. S. Price, J. C. Cavendish, R. S. Varga

This paper presents a new technique for solving some of the partial differential equations that are commonly used in simulating reservoir performance. The results of applying this technique to a simple problem show that one obtains accurate pressure values near wells, as well as accurate pressure pdients, which can be explicitly calculated. The method is completely rigorous in that convergence of the discrete numerical solution to the continuous solution for both pressure and pressure gradient is established. High-order, piecewise-polynomial approximations are used near the wells where pressure gradients are steep, while low-order, piecewise-polynomial approximations are used elsewhere to reduce greatly the calculation time. This combination is shown to give a uniformly good approximation to the solution. These approximations, obtained by using a Galerkin process with suitable Hermite subspaces, are shown to be theoretically and numerically superior to the usual approximations obtained from standard finite-difference techniques. Not only are much greater accuracies obtained, but computer times are also greatly reduced. The application of this technique to multiphase flow problems (e.g., single well coning problems) would have considerable practical interest, but such extensions of this technique with full mathematical rigor have not been made as yet. However, the numerical methods presented here are general, and in principle extend to multidimensional, multiphase flow. Moreover, the preliminary results given in this paper are sufficiently encouraging that we feel the effort in attempting these extensions is justified. INTRODUCTION The problem of obtaining accurate pressure distributions and pressure gradients around wells is of considerable importance in the numerical simulation of reservoir performance. The most common approach to solving this problem is to use finite difference techniques (see McCarty and Barfield9 or Peaceman and Rachfordl0). This approach, however, has many disadvantages, the major one being that many grid points are generally necessary for accurately describing the pressure distribution and the pressure gradient around wells. This need for a fine grid results in large computer times and often in prohibitively high costs. Besides investigating the method of finite differences, some authors, such as Welge and weber17 and Roper, Merchant and Duvall,13 have considered a combination of analytical and numerical techniques with some success. These approaches, however, are all nonrigorous and quite often cannot even be applied. In this paper, we present a numerical formulation of high-order accuracy, based on the Galerkin method, for solving this problem. We treat here only the partial differential equation that describes steady-state, single-phase flow. However, the methods presented are general and in principle extend to multidimensional, multiphase flow. Specifically we treat special cases of the problem in two dimensions described by: where G is a rectangle (with sides parallel to the coordinate axis) with boundary ?G, ?/?n denotes the outward normal, and y and B are non-negative constants such that y + B > 0.

Jan 1, 1970

By D. R. Shreve, L. W. Welch

An analysis for predicting the behavior of reservoirs exploited by a combined gas drive and gravity draintrge mechanism is presented. The method allows the prediction of gas-oil ratio, oil production rate, and cumlulative oil recovery for high relief pools in which pressure is maintained either by gas injection in the crest of the structure or by a combined gas and water injection program. The reservoir is represented by a model of communicating vertical prisms, with production or injection at the upper or lower boundaries of each prism. The theory of macroscopic gravitational segregation of two fluids within a reservoir is discussed in detail for all ranges of saturation. A description of an application of the analysis to a typical field problem is also presented. INTRODUCTION The importance of gravity drainage as a producing mechanism has become more evident during recent years. However, many of the reservoir engineering methods which have been used in the past for predicting the behavior of reservoirs producing by this mechanism have required such drastic, simplifying assumptions that the calculated performance in some cases does not resemble the actual pool behavior. The mechanism of fluid displacement in porous media has been described by Buckley and Leverett,1 and a simplified "one zone" method based on this development has been used by engineers to evaluate the performance of gravity drainage pools. Welge2 has modi- fied this approach so that the calculation time has been markedly decreased. In addition, many other authors3-9 have analyzed various phases of the gravity drainage producing mechanism during the past several years. In a high relief reservoir in which gravity may affect the recovery mechanism, there are two main types of exploitation programs which can be employed: 1. Reservoir pressure is maintained by injecting gas into the upper portions of the reservoir. In this case gas is not released from solution and, if dynamic pressure drops are small compared to the total pressure of the system, the gas and oil may be considered to be immiscible and to have constant fluid characteristics. 2. Reservoir pressure is not maintained and a "natural depletion" mechanism prevails. The term "natural depletion" can be used to describe the production history of a reservoir in which a gradual decline of average pressure and the associated release of solution gas continue throughout the producing life. In this mechanism gas is released throughout the reservoir and there is a continuous change in the characteristics of both the gas and the oil phases. It is the purpose of this paper to present the results of mathematical developments by which a method has been devised for evaluating the performance of reservoirs which are depleted by the first type of exploitation program discussed above. It should be mentioned that, although throughout the analysis gas and oil represent the two fluids existing in the reservoir, the method is equally applicable to gravity drainage in water-oil systems. STATEMENT OF PROBLEM In reservoirs of pronounced structural relief, where the closure of the oil section is great or the angle of

Jan 1, 1957

By R. Simon, D. J. Graue

This paper presents correlations for predicting the solu-lility, swelling and viscosity behavior of CO2-crude oil sys8i.m~. The correlations were developed from experimental data obtained by the authors. These data are also presented. The data were determined by measuring the properties of mixtures of CO, and nine different oils. Experiinental conditions covered a range of 100 to 250°F and pressures up to 2,300 psia. Properties predicted by the correlations have average deviations, expressed as per cent of experimental value, of 2 per cent for solubility, 0.5 per cent for swelling and 12 per cent for viscosity. INTRODUCTION Interest in CO2 injection as an oil recovery process has led to the development of performance prediction methods which can be applied to specific reservoirs.1 iS To use these performance prediction methods, it is necessary to know the solubility, swelling and viscosity properties of CO2-crude oil mixtures at reservoir conditions. Some information on these properties has appeared in the literature; however, this information did not cover the range of different oils and conditions needed to prepare generalized correlations for reservoir engineering purposes. Consequently, an experimental program was undertaken to collect the data needed. The data obtained and the correlations developed from the data are described in the following sections of this paper. SOLUBILITY OF CO2 IN CRUDE OILS CO2 solubility data in the literature come from six principal sources. The solubility prediction method of Welker and Dunlop3 is limited to 80F. The information in Ref. 4 is of two types: the first includes binary and ternary mixtures of CO, and light hydrocarbons (C1 to C6), and the second gives data for CO1 and heavy hydrocarbons for a temperature range of 40 to 90F. Ref. 5 contains a KCO2 chart for systems whose convergence pressure is 4,000 psia. The KCO2's are based mainly on CO2-natural gas mixtures. Poettmann's work covered CO2 solubility in one condensate and one crude oil6,7 Ja-coby and Rzasa measured CO? solubilities as a function of pressure and temperature for two natural gas-absorber oil mixtures and two natural gas-crude oil mixtures.'8 CO, concentration in these four systems was fixed at 5 mol per cent. The work reported in this paper extends CO, solubility data to a variety of different crude oil types in a temperature range from 110 to 250F and pressures up to 2,300 psia. The experimental procedure used by the authors to obtain the solubility data consisted of combining known amounts of pure CO, and crude oil in a visual cell at a fixed temperature and measuring the bubble point of the mixture. Measurements were made for a total of 40 different CO2-oil mixtures and the results are shown in Table 2. The mixtures included nine different oils (seven crude oils and two refined oils) whose properties are listed in Table 1. All nine oils had vapor pressures less than 1 atm at the experimental temperatures. Consequently, analysis of the bubble-point vapor showed a CO, concentration over 99 mol per cent. At no time during these experiments was a second, more dense, liquid phase observed. The solubility correlation which was developed from the data in Table 2 is presented in Figs. 1, 2 and 3. In these figures, solubility is expressed as xco2 the mol fraction of CO, in the CO2-oil mixture. Fig. 1 shows solubility as a function of CO2 fugacity6 and temperature. Fig. 2 shows the same solubility data expressed as a function of saturation pressure and temperature. The solubility shown in Figs. 1 and 2 is for an oil whose UOP characterization factor is 11.7. UOP characterization factors of crude oils can be determined from Ref. 10 if the viscosity and APT gravity of the oil are known. Fig. 3 gives the solubility correction factor for oils whose UOP characterization factors differ from 11.7. The solubility correlation in Figs. 1, 2 and 3 predicted

Jan 1, 1966