Search Documents

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

By B. H. Caudle, A. B. Dyes

In a recent publication' it was shown that wells with a free surface in a homogeneous gravity-drainage reservoir have a hyperbolic decline with index n '. This paper reports efforts to extend this theory to practical field cases. As a first step in making the extension, a large number of lease production curves from gravity-drainage fields in the stripper stage were fitted with the general hyperbolic decline equation, q/q, = 1/(1 + y) y and n in each case being determined by statistical methods. It was found that this equation is reliable for decline curve extrapolation for periods up to at least 25 years. The resrilts obtained were disturbing in the sense that values for the index, n, were often smaller than one, which fact leads to the prediction that a well would produce all infinite amount of oil in infinite time. A theoretical investigation based on the results of the previous paper provides a suitable explanation for these low values of n. It suggests that these are obtained if production occurs from two or more layers of diflerent permeabilitv and thicktress or two lnyers having different skin effects. It is predicted that as the more permeable layers are exhausted, the index of decline, n, will rise from a value smaller than 2, to n = 2. A practical method is given for fitting the general hyperbolic decline Eq. 6 to production decline data. When this method is used, the index of decline observed in gravity-drainage reservoirs in the stripper stage when wing only data over early years 01 decline usually agrees quite closely with that determined using data over 25 years or more. This indicates that the hyperbolic decline equation is reliable for predicting behavior of such reservoirs. A method is also suggested for taking into account any changes in the value of the index INTRODUCTION As a result of both experimental work with models and theoretical deductions, it was recently concluded' that wells having a free surface in a homogeneous reservoir producing by pure gravity drainage should have a hyperbolic decline with n = 2. The present report details efforts to extend this work on ideal homogeneous reservoirs to actual field cases. It should be remembered that both the previous work and the present extension are based on the assumption that the reservoir is at very low pressure, and hence, nearly all the potential available for driving fluid to the well is the gravity head. Attempts to apply conclusions of this study to reservoirs which have considerable pressures supplementing the gravity drive may lead to difficulty. WELLS WITH NO FREE SURFACE In a dipping gravity-drainage reservoir such as that shown in Fig. 1, some of the down-dip wells have no free surface. The previous study' showed that the up-dip wells, which have a free surface, should follow Eq. 1. The following approximate analysis, of a type which was justified experimentally for up-dip wells in the work just mentioned, leads to results for the down-dip wells.

By Potts N. L., Eilerts C. K., Summer E. F.

The second-order, nonlinear, partial- dillerential equation representing the transient radial flow of gas - condensate fluids in reservoirs has been integmted by using finite -dillerence equations and electronic computers. Effect was given to pressure -dependent permeability, viscosity, and compressibility and to distance-dependent permeability. The influence of a second-degree velocity term in the Darcy equation was investigated. Implicit methods were used and practical, convergent solutions were obtained with material balance to less than 6 x 10-4 for recovery of one- hall the reserve at constant flow rate. Integration results provide the productive period of a reservoir for a given constant rate and the fraction of the fluid initially in place that can be recovered in that period. The properties of a lean and a rich fluid are represented in a set of integrations designed to demonstrate the effect of different constant-recovery rates and significant variables emphasized one at a time. INTRODUCTION A method is needed for computing the availability or reserve of a gas in a formation, making use of all technical and engineering information that is pertinent. As a step in this direction, a program for computing transient linear flow1 was developed that utilizes principles of earlier finite-difference computing for a similar purpose2-4 and gives effect to significant pressure- and distance-dependent properties of reservoirs and their contents. The present paper pertains to transient, radial flow of gas-condensate fluids. In the partial-differential equation of flow, effect is given to the variables viscosity and compressibility,5 and also to permeability. Included in the study is the quadratic form of the Darcy equation of flow that has been the subject of field tests6 and that has been applied to a gas with constant properties in transient flow computing. 7 Inclusion in the differential equation of variable coefficients to represent properties greatly complicates the higher derivatives of the equation. Because it was impractical to make a required improvement in certain of the finite-difference derivatives by Taylor-series expansion, five-term derivatives were used in the implicit computing. Related methods were developed that can improve general facility in manipulating finite-difference equations. BASIC EQUATIONS PARTIAL DIFFERENTIAL EQUATION The basic partial differential equation for transient radial flow in the direction of decreasing radius is l_ a(vpr) _? .......(1) r gr at where r = radius, L v = apparent velocity aIong radius, L/t p = density, M/L3 t = time, t Ø = porosity, dimensionless. In this equation, only the porosity is regarded constant. The density is where p = pressure, m/Lt2 M = molecular weight of fluid, m T = constant temperature of fluid, T z(p) = pressure-dependent compressibility factor, dimensionless R = gas Constant, mL2/t2t. The physical boundary of the reservoir is at re, and the producing boundary is defined to be at tq with rq = re /200. By definition,

Jan 1, 1966

By M. H. Gaskell, D. M. Kehn, G. T. Pyndus

A correlation of the bubble point pressure for black oil systems is developed using the standard physical-chemical equations of soiutions. The correlation is based on 158 experimentally measured bubble point pressures of 137 independent systems and is expressed in terins of the usually measured field pnrameters— flash separation gm-oil ratio, rank oil gravity, total gas gravity, and reservoir temperature. The data were obtained on systems produced in Canada, Western and Mid-Continent01 United States, and South America. The average error (algebraic) in the representation is 3.8 per cent, and the maximum error encountered is 14.7 per cent. INTRODUCTION In the absence of experimentally measured properties of reservoir fluids, it is often necessary for the field engineer to make estimates regarding the fluid properties based on the usually measured producing parameters. To aid in these estimations, various correlations have appeared in the literature in recent years. Among the pertinent properties of interest is the bubble point pressure. A correlation for this parameter has been reported by Standing'. However, this correlation was based essentially on California produced crudes and this limitation was pointed out with its presentation. The correlation presented in this paper utilized data on crude oil systems from Canada, Wes- tern and Mid-Continental United States, and South America. CORRELATION DEVELOPMENT The basic assumption used in this development is the same as employed by Standing', There is a wide variety of ways to combine these parameters; however, in this instance the combination was made on the basis of Henry's law2. Accordingly, the relationship proposed is Although Eq. 2 defines an individual system, it is of limited value since H' is a function of gas-phase composition and the system temperature. It was observed that for the systems where the bubble point was measured at several temperatures that the ratio of the bubble point pressures and the ratio of the corresponding absolute temperatures (R) were practically identical. Thus, for correlation purposes the bubble point pressure may be taken as a direct function of the absolute temperature. This relation is valid only for those systems that are not near the critical point. Accordingly, this correlation will be inadequate for systems in the region of the critical point. The solubility of the various hydrocarbons found in the gas phase increases with the molecular weight. Thus, the saturation pressure should be inversely related to the gas gravity. Applying these principles to Eq. 2 and rearranging terms gives: The variables 02 the left side of Eq. 3 were designated as the "hubble point pressure factor". The number of mols of tank oil per barrel is a function of the "mclecular weight" of the tank oil. Although the tank oil is a complex mixture, it was assumed for the purposes of this correlation that a unique molecular weight could be assigned to a given oil. This was designated as the "effective molecular weight", and was related to the oil gravity, This empirical relationship was developed simultaneously with the correlation by assuming values of M,. and working to obtain a smooth curve for both the correlation and the effective molecular weight. The relationship between the oil gravity and the effective molecular weight used in this correlation is shown in Fig. 1. The effective molecular weight is somewhat higher than the molecular weight of the C fraction. The difference between these values is largest for the low-gravity systems. It is noted that this effective molecular

By B. H. Caudle, L. H. Silberberg

Reservoir depletion by natural water drive is typified by the movement of water from an aquifer into the adjacent oil-bearing formation. Prior studies of this type 01 water movement have generally neglected the resistance to flow in the aquifer and its effect on the movement of water into the oil bearing zone. A method for designing and operating scaled models of such reservoir systems is presented. Experimental data on a model of an edge-water-drive reservoir are shown and discussed. INTRODUCTION Most petroleum reservoirs derive at least a part of their productive capacity from water influx. This water may be injected from the surface or it may come from an aquifer adjacent to the oil zone. In either case, the reservoir engineer must be able to estimate the advance of the water as a function of either the elapsed time or the fluid produced. This paper describes a type of fluid flow model which includes the effect of viscous and gravitational forces in the reservoir and the surrounding aquifer. Natural water influx can be divided into three general (but widely overlapping) categories according to the direction of flow in the aquifer. These are shown schematically in Fig. 1. Fig. l(a) illustrates an edge-water-drive mechanism. In this case, water advances updip along the stratum but the movement of water is essentially horizontal and very little of the oil is actually underlain by water. The bottom-water drive (Fig. l(b)) is characterized by a thick aquifer underlying the oil zone. The water movement is generally vertical in the aquifer. The third category is illustrated in l(c). Referred to as the "thin oil column", this type of oil accumulation consists of oil over water in a relatively thin stratum. Fluid movement is horizontal in both the oil and water zones, except close to the producing well. One characteristic which all three of these types have in common is that part of the water flow takes place in the water invaded region of the oil reservoir while the remainder of the water flow occurs in the aquifer. The resistance to the flow of water usually will not be the same in these two regions. This occurs because the porous rock normally contains only water in the aquifer, while the water-invaded region always contains microscopic globules of bypassed oil (residual oil) which interfere with the flow of water through the rock. The production history of any type of a water-drive system is a function of two phenomena — where the water goes, and how it displaces the oil in the area invaded by the water. It is usually impractical to study these two phenomena at the same time. In most cases one type of analysis is used to predict the gross movement of water in the reservoir, while a different type of analysis is employed to determine the amount of oil to be recovered from the gross volume contacted by the water. The portion of the oil reservoir which is invaded by water is mainly a function of the resistances to fluid flow in the several parts of the flow system. This gross water influx in the three types of natural water drive typically results in water cusping to the production wells either from the side or from below (water coning). In 1747, Muskat described a greatly simplified mathematical model for predicting the area invaded by water in a bottom-water-drive reservoir.1 This model used the Laplace equation and suitable boundary conditions to describe the isopotentials and streamlines in the flow system. The use of this model assumes that (1) the reservoir rock is homogeneous in nature, (2) the oil and water mobilities are equal, (3) there is no oil flowing in the water invaded region, and (4) external forces — such as gravity — do not affect the flow. This method has been used for both natural water drives and water injection projects. In this, as in most models, the method of images is used to reduce the size of the model necessary to describe the reservoir. The assumption of equal viscosities and permea-

Jan 1, 1966

By L. H. Simons, H. H. Spain

I NT RODUCTION Theory indicates that linear water-floods should exhibit scaling and stabilization properties in both oil-wet and water-wet porous media'. Experimental verification of these proper-ties in oil-wet media was obtained some time ago', and more recently, Perkins', employing a high permeability unconsoli dated sand and a low oil-to-water viscosity ratio, has obtained flooding data in a water-wet medium which follow the pattern of scaling and stabilization. Also, Root and Calhoun3 have found similar trends in the displacement of gas by liquids from unconsolidated sands. In general, however, the experimen-tal evidence appears still somewhat conflicting since scaling and stabilization of water floods in water-wet media has not been observed in most of the studies to date'.".'. The purpose of this paper is to provide a more comprehensive picture of waterflood behavior in water-wet media and to corroborate the validity of the scaling and stabiliza-tion concepts. The need for differentiating between the intrinsic nature of water-oil displacements inside a porous medium and perturbating secondary end effects is discussed, pertinent experimental procedures are outlined and the results of flooding tests conducted at low and at high oil-to-water viscosity ratios on con- solidated, water-wet Alundum cores are analyzed. END EFFECTS IN WATER-WET MEDIA Outlet End Effects The peculiar feature of the outlet end effect in water-wet cores is that it not only results in an excessive wetting phase saturation at the outflow face,'.' but that it also retards the moment of water breakthrough: i.e., causes an undue delay in the appearance of water in the produced effluent. Conditions existing in a water-wet core when water first arrives at the outflow face are schematically shown in Fig. 1. At this flooding stage the saturation distribution, Fig. la, has not yet been distorted by the outlet end effect. A water-invaded pore at the outlet core face, opening into the oil-filled void space between the end of the core and the end plate of the flooding cell, is illustrated in Fig. lb. The curvature of the water-oil interface is concave toward the outlet, and the pressure on the water side of the interface is lower than on the oil side by an amount equal to a certain capillary pressure (Fig. lc). Before the water (of this particular pore) can be produced from the core, it is apparent that the pressure in the water phase must exceed that in the oil outside the core. Accordingly. when water first reaches the outlet face, it will not leave the core but will accumulate inside, and only oil will be produced. Eventually, as the water injection is continued. sufficient pressure is built up in the water to reverse the curvature of the interfaces in the outlet pores, and only then will water be produced from the core. This stage, illustrated in Fig. 2. corresponds to the conventional "water breakthrough" observed in flooding experiments. For purposes of discussion, it is convenient to distinguish between the flooding stage at which water first arrives at the outlet face of a core (Fig. 1) and the stage at which water is first produced from the core ( Fig. 2). The former will be referred to as "water arrival" and the latter water breakthrough. In oil-wet media

By Homer N. Mead

This report develops equations for decline curve analysis based upon the premise that the rate of change of the reciprocal of decline for succeeding time intervals is constant when the reservoir is produced under a fixed set of conditions. A method is shown for predicting maximum production rate against cumulative production for any reservoir. A method is also presented for predicting future production rate by an analysis of past production performance after decline has been esrahlished. INTRODUCTION No significant contribution to the analysis of decline curves by the loss-ratio method has been made since J. J. Arps paper1. The concept of loss-ratio that was developed by Arps and his contemporaries has been redefined in another way and the concept of instantaneous loss-ratio at time zero has been added. In the mathematical developmerit by Arps, production rate was related as a continuous function with time and his equatlons were developed on that basis. In this paper, production is considered to be a series of segments for equal time intervals and equations have been developed based upon these finite differences. It is realized that decline curve analysis is not the answer to all predictions of reservoir behavior. However, as production must decline from an initial maximum rate to zero in any reservoir, if such decline can be expressed as an infinite series, this series should accurately predict production. Decline curve analysis should be considered a valuable tool that may be used in conjunction with predictions of future recoveries by other methods. Various uses of the decline curve method will be discussed in this paper. Equations predicting production rate at any time and cumulative production by exponential and hyperbolic decline will be developed. Characteristic values of b related to various types of drive will be discussed. After that a hypothetical reservoir will be studied and ultimate recovery by natural depletion and pressure maintenance compared, then actual field examples will be shown and discussed. The rate-cumulative curve introduced by H. N. Marshs in 1928 is an important one for use in estimating future production rate and ultimate recovery. This graphical method is the one that is used in conjunction with the equations to be developed for predicting production behavior. EXPONENTIAL DECLINE In exponential decline each succeeding production rate per unit of time is a constant percentage of the production rate just before it. Since (1 — r) is the common ratio, where r is the constant decline rate, this condition may be expressed as a geometrical progression; thus Pn = P (1 - r)n-1 . ......(1) Where r is less than unity, r Since P1 (1 - r)n = Pn+1 P1 - Pn+1 C =P1-Pn+1/r........(2) LOSS-RATIO In this paper the calculation for loss-ratio has been changed from the method presented by Johnson and Bollens2. The first loss-ratio to be calculated is a,; it is the initial production per unit of time after decline has set in divided by the difference between the initial and the second production rate. This loss-ratio is for the first time interval. After b (which is the constant change

Jan 1, 1957

By W. A. Burns

Vertical flow is an important mechanism in many petroleum reservoirs. Yet no adequate method has heretofore been proposed for determining the in-situ vertical permeability in a reasonable length of time. Core measurements of vertical permeability do not necessarily represent effective in-situ values. Cores may be representative of only localized pockets of high or low permeability, and total reliance on cores might cause one to overlook the presence of vertical fractures altogether. Conventional single-well tests, interference tests, and pulse tests commonly are used to determine horizontal transmissibility, khh/u, but generally are not designed to measure vertical trans-missibility or other groups containing vertical permeability, ku. A new well test described here provides a measure of in-situ vertical permeability as well as the previously measurable horizontal permeability. This "vertical test" determines the vertical diffusivity, kv/gcu, and horizontal transmissibility, khd1/u (where dl is the length of injection or production perforations). The test is basically a short-term interference or pulse test normally lasting less than a day. Conventional down-hole equipment is used between two isolated sets of perforations in the same well, as shown in Fig. 1. A "layer-cake" description of a reservoir can be obtained by running two or more vertical tests in the same well, each over different vertical intervals. Each test can be analyzed for the effective vertical and horizontal properties of the interval tested. In addition, a vertical well test can determine the effective- ness of a shale barrier or tight zone by testing across it. Field use to date of vertical testing has confirmed its utility in reservoir description. Mathematical Basis Interpretation of vertical tests for ku and kh is based on comparison of measured pressure response with computer-generated response. Pressure behavior during a vertical test in a homogeneous, anisotropic, infinite-acting reservoir (infinite in radial extent and bounded in both directions vertically) is given by the diffusion equation for slightly compressible fluids with negligible gravity effects.l A schematic diagram of the problem is presented in Fig. 2, and an analytic solution, based on the image

Jan 1, 1970

By A. M. Rowe, I. H. Silberberg

A computer program was written to predict the phase behavior generated by the enriched-gas-drive process. This program is based, in part, on a new concept of convergence pressme, which is then used to select vapor - liquid equilibrium ratios (K-factors) for performing a series of flash calculations. The results of these calculations are the equilibrium vapor and liquid phase compositions which define the phase envelopes. The program was used to predict phase envelopes for 11 different hydrocarbon systems on which published experimental data were available. The predicted and experimental results compare favorably. INTRODUCTION The reservoir engineer is frequently faced with the problem of predicting what will happen if gas is injected into a reservoir. One aspect of this, general problem is predicting the phase changes that will occur when a non-equilibrium gas displaces a reservoir fluid. In particular, a "dry" gas, upon displacing a volatile oil, will pick up intermediate components from the oil. On the other hand, a "wet" gas, containing a high concentration of intermediates, will lose some of these components to a relatively low-gravity, non-eouilibrium1- crude. It is this latter Drocess. occurring in the enriched-gas displacement, which is treated in this paper. In the past, these phase changes have been determined experimentally and the results incorporated into various modifications of the Buckley-Leverett analysis.112 Such experimental work is time consuming, and the results are sensitive to numerous experimental errors. With the wide availability of high-speed digital computing equipment and numerous correlations pertaining to the vapor-liquid equilibria of hydrocarbon systems, it is now practical to calculate such phase behavior. This paper describes a computer program for performing these calculations. THE ENRICHED GAS DISPLACEMENT PROCESS Experimental results have shown that oil recovery can be significantly increased by enriching the displacing gas with intermediate hydrocarbon c0m~onents.3 The essential features of the phase behavior generated by this enriched-gas-drive process are commonly illustrated with ternary diagrams such as Fig. 1.4 In this figure, Gas D, which contains a high concentration of intermediate hydrocarbons with respect to the undersaturated Crude A, is injected into the reservoir. When D contacts A, gas goes into solution until the oil becomes saturated (Point. B). Further contacting of Gas D and saturated Oil B results in a Mixture C which separates into Vapor Y(c) and Liquid X(c). Liquid X(c) is contacted by additional Gas O, resulting in Mixture E which separates into Vapor Y(e) and Liquid X(e). Repeated contacts of the liquid by the injected gas will eventually result in Liquid X(4 of maximum enrichment existing in equilibrium with Gas Y(d). The equilibrium tie-line X(4 Y(4, when extended, passes through the Point D representing the enriched injection gas. For systems of more than three components, the predicted equilibrium states are dependent upon not only reservoir temperature and pressure, but also the compositions of the crude oil and injected gas. If the gas is sufficiently enriched, a miscible displacement is generated. Line If is tangent to the phas,e envelope at the critical point (Point Z) and represents the limiting slope of the tie-lines as the critical state is approached. Point 1 therefore represents the minimum enrichment of injection gas required to generate a miscible displacement. Point G represents the minimum enrichment required for initial miscibility of the injection gas with Crude A. Accra has presented a method to be used for prediction of oil recovery by the enriched gas displacement process.l To develop the phase behavior data needed, he designed the experimental procedure described in the following quotation from his paper: The original liquid was contacted by a volume of displacing gas and allowed to come to equi-

Jan 1, 1966