Search Documents

By H. T. Kennedy, S. G. Dardaganian, J. L. Clanton

This paper presents equilibrium ratios measured on oil and gas samples taken from the Kelly-Snyder field at the reservoir temperature of 128°F. From these data smoothed curves were constructed for the reservoir fluid up to the bubble-point, 1,874 psia. A material balance calculation, based on the use of equilibrium ratios instead of the usual solubility and shrinkage meayr~rements is made on the field and compared to the con\.entional method. DETERMINATION OF EQUILIBRIUM RATIOS The equilibrium ratios, or K values were obtained by recombining separator samples of oil and gas taken when the reservoir was above its bubble point, and thus represented the reservoir fluid before any change had occurred due to the evolution of gas. The well sampled had a gas-oil ratio of 617 cu ft of high pressure sepalator gas per barrel of separator oil. and a bottom-hole temperature of 128°F. Samples were brought to equilibrium by agitating for 24 hours in a thermostat, maintained at the proper temperature within 0.2ºF, mercury being added when necessary to maintain the pressure constant. Samples of oil and gas were then taken by constant pressure displacement, and analyzed by low temperature frac-tionation. Equilibrium ratios were determined at 128°F and 500, 800, 1,000 and 1,400 psig. At the latter pressure, there were indications that the gas samples were contaminated with small quantities of oil, and only data taken at the first three pressures were included in Fig. I. This figure has been extended to higher and lower pressures by the use of the Fluor charts for 3,000 psi convergence pressure, on which the experimental points fit well. MATERIAL BALANCE EQUATIONS The possible errors involved in conventional material balance equations have been analyzed by Muskat1. Where the errors occur in field measurements, or in obtaining effective averages from variables such as gas-oil ratio and pressure, they will evidently be involved in any attempt to estimate the content of the reservoir. On the other hand, substantial errors may be introduced in conventional material balance calculations by: (1) the neglect of the amount of oil dissolved in gas at high pressures and temperatures in the reservoir, but which condenses under separator conditions; (2) the assumption that the separation of oil and gas in the reservoir may be duplicated by a differential liberation in the

By H. N. Mead

In a reservoir when gas comes out of solution and rises, additional pressure is created because of the change in position of this gar in the bounded volume. If this pressure effect is not taken into account as a pseudo-influx ill mate rial balance calculations on reservoirs in which there is evolving gas, an error is introduced that is directly proportional to the length of oil column and amount of gas evolved and inversely proportional to the reservoir pressure. In reservoirs with small oil column this pseudo-influx is of little importance. The larger the oil column the more important the effect becomes. The most important application, however, will be for reservoirs with declining pressure in which increasing water entry has been indicated by material balance calculations. It could be that this indication is the result of the cited pseudo-influx and not water entry at all. This paper attempts to explain this phenomenon in relation to material balance calculations and prevents an approximate method of determining tile pseudo-influx. INTRODUCTION Only recently has the pressure-volume relationship of gas rising in a liquid appeared in engineering literature, although its effect has been discussed for several years—especially among engineers engaged in drilling operations. The writer published a paper in 1957 in which this phenomenon was examined in relation to blowouts in drilling wells.' Also, in 1957 Stegemeier and Matthews' authored a paper in which the relationship was presented with respect to pressure build-up tests when shutting in producing wells. In this paper the effect will be called "pressure inversion" which is descriptive of the phenomenon. Fundamentally, the effect just mentioned is the result of changing the position of a quantity of energy contained within the gas itself (as a result of its compression) within the bounded volume in which the gas exists. If the gas were allowed to expand "beyond all bounds" there could be no pressure inversion effect. This is because all of the energy contained within the gas as a result of compression would have been re- PRESSURE INVERSION A phase that has not been discussed in the literature is the effect of pressure inversion upon material balance calculations in reservoirs in which the reservoir pressure drops with reservoir withdrawal and is below the bubble point. This action releases solution gas which, after the critical gas saturation has been reached, will tend to rise and collect at the top of the reservoir. The pseudo-influx caused by this effect could be great enough in reservoirs of large oil column and high permeability to render invalid any analyzer studies made on them — if pressure inversion has not been taken into account. Most suspect would be those reservoirs with large oil column in which increasing water entry has been predicted but the bottom wells have not yet produced water. Example 1 Simply stated, this phenomenon exists in oil reservoirs because of the difference in density and expansibility between gas and oil. For example, if a bubble of gas exists at the bottom of a 1,000-ft column of oil having a gradient of 0.30 psi/ft with the pressure at the top of the column zero, the bubble of gas will be compressed to a pressure of 300 psi. If temperature and the combined volume of the oil column and the gas bubble remain constant, as the gas bubble rises in the oil column it will have the same volume and consequently the same pressure at any point. Under these conditions when thc bubble reaches the top of the oil column, the pressure will be approximately 300 psi at the surface and approximately 600 psi at the bottom of the oil column. If it is assumed that the oil is incompressible the pressures will be exactly 300 and 600 psi, respectively.

By J. G. Eenink, R. A. Cunningham

One phase has been completed of a laboratory invesrigation of formations with relatively high permeability under conditions of overburden, formation and mud coltrmn pressures. The following statements are based on these tests within the limitations described. Drilling rate decreased when mud column pressure was greater than formation pressure. The decrease was primarily due to a layer of cuttings and mud particles held to the hole bottom by the difference in pressure. Much of the force of the bit was wasted in this layer and was unavailable for penetrating virgin formation. Adequate jet velocities helped clean the filter cake and chips away and resulted in increased drilling rate — the higher the jet velocity the faster the drilling rate. Drilling rate increased slightly when formation pres sure was greater than mud column pressure. The i11(reuse resulted from cleaning the hole bottom as formarion fluid flowed into the borehole. Overburden pressure had practically no effect on drilling rate. INTRODUCTION Some formations are difficult to drill in a wellbore hole but are easily drilled under atmospheric conditions in the laboratory. The reasons for this apparent difference in drillability are not fully understood. An increase in formation compressive strength causes a reduction in drilling rate. It has been shown by Griggs1, Handin2 and others3-1 that formation strength increases under conditions of high hydrostatic stress. Based on these results, a decrease in drilling rate would be expected with increased drilling fluid pressure in formations of low permeability having low formation fluid pressure. Such a decrease has been shown by Murray. et a1 and Eckel" This effect holds even when using nitrogcn under high pressure as the drilling fluid (unrc-ported tests from the Hughes laboratory). Poor bottom-hole cleaning may account for decreased drilling rates. Even though it may not be easily explained, evidence of poor cleaning can easily bc seen. One evidence is the increased drilling rate that accompanies decreased balling obtained with addition of oil in some muds used to drill shales.4,7,3 Many variables can affect drilling rates. These may include mud composition and properties, formation properties, bit types and operating conditions, formation fluid, mud column and overburden pressures. These must be studied a few at a time under carefully controlled conditions which simulate actual conditions as closely as possible. The purpose of this work was to investigate effect of overburden, formation and mud column pressures on drilling rates in permeable formations. These tests were suggested by conditions shown in Fig. I. Represented is a borehole drilled into a formation with relatively high permeability. In the vicinity of the hole bottom there exists a mud column pressure, P a formation fluid pressure, P and an overburden pressure, The mud has low water-loss properties and sufficient weight to more than overbalance the formation pressure. Laboratory investigation included a range with formation pressure greater than mud column pressure. NOMENCLATURE Pm mud column pressure = pressure of drilling fluid or mud at hole bottom Pf, formation pressure = pressure of fluid in intersticies of formation

By J. R. Elenbaas, J. A. Vary, D. L. Katz

The development of gas fields, oil fields and aquifers for storing natural gas is treated from two main vieu.-points: (I) the volumetric storage capacity for gas in a given situation and (2) the prediction of the number of wells required for the delivery of gas. Other experiences in the design and operation of storagc fields are incluclerl INTRODUC TION Storage of natural gas in underground reservoirs near the terminus of long distance pipelines has been the prime factor in opening the space heating market to the natural gas industry. Storagc has permitted a major. increase in both the load and the load factor of pipe-lines; some are now operating at steady load throughout the year. Thus, underground storage has been responsible for the rapid increase in demand for natural gas in recent years. Three types of reservoirs have been used for gas storage: natural gas reservoirs, oil reservoirs, and waterbearing sands or aquifers. This paper presents the factors to be considered when developing gas storage reservoirs of these three categories. There are two prime considerations tor any storage reservoir: (1) the volume of gas which a given reservoir will store advantageously and (2) the number 01 wells needed to provide the required peak deliverability. These two problems will be considered for the three types of reservoirs just noted STORAGE IN PARTIALLY DEPLETED GAS FIELDS Early storage operations consisted of replenishing the natural gas in a depleted gas field situated adjacent to the market. Today, newly discovered fields near the market may be considered for storage, and this discussion applies equally to both types of reservoirs. For reservoirs originally containing gas or oil, the question of the impermeability of the cap rock nor-mally does not arise. However, such fields are likely to have many wells drilled either to or through the reservoir under consideration. Positive assurance must be obtained that such wells are or can be made mechanically tight. Corroded casings may need to be lined or permanently plugged. Abandoned wells should bc reopened and properly cemented. The volumetric capacity for gas storage depends upon space available in the porous rock as well as pressure and temperature of the gas in the reservoir. The production-pressure decline data on partially depleted gas reservoirs without water drive permit calculation of the reservoir space for gas. Isopachous maps of sand volume and porosity data for the reservoir rock provide an alternate method of calculating the pore volume for water-drive reservoirs. The pressure range selected for the storage cycle depends upon ()) the safe upper limit of pressure. 2) the flow capacity of wells and (3) compression requirements when injecting gas into the reservoir or delivering to market. Normally, gas and oil fields have pressures at discovery in the range of 0.43 to 0.52 psi/ft of depth. Pressures of around 1.0 to 1.2 psi/ft of depth appear to lift the overburden1-3 and invite uncontrolled movement of fluids in the porous rock. Some top pressure is normally selected for a storage reservoir ranging from below discovery pressure for deeper reservoirs to 0.65 psi/ft of depth for shallower reservoirs. Pressures to 0.66 psi/ft have been experienced without difficulty. The lower pressure limit is set by water intrusion accompanying low pressures, reduced flow capacity for wells at lower pressures and compression requirements. Depletion-type gas reservoirs often encounter water problems in the later stages of gas production. Such water intrusion may be due to movement from the surrounding aquifer. Accordingly, displacement of this water back into the aquifer by gas pressure and subsequent surges of water corresponding to the gas storage pressure cycle must be considered. Storage fields often produce in four months a volume of gas equal to its initial content. Rapid decreases in reservoir pressure occur, such as 20 psi/day. Accordingly, closed-in pressure observation wells which reflect the pressure in the bulk of the reservoir are required for following the operation of the reservoir. It has been found that a plot of observation wellhead pressures against gas content, Fig. 1, is very useful in observing operation of the field, checking the inventory and predicting future behavior. The plot is based on a given quantity of base or cushion gas in place. The injection and withdrawal curves may spread depending upon the homogeneity of the reservoir rock. permeability of the rock, well spacing and flow rates.

By L. D. Mullins, E. B. Elfrink, W. L. Wahl

In the past few years several articles and papers presenting results of solution gas-drive depletion calculations have appeared in the lit-erature. Such calculations are of interest to the oil industry, for investment decisions must often be made before much is known about a reservoir. At other times, an estimate of the possible benefits to be realized from alternate production methods is desirable, and theoretical depletion calculations can serve as a floor or reference level from which to work. In any case, an estimate of ultimate oil recovery based upon engineering data is commonly required. An engineer confronted with the problem of obtaining, for a specific reservoir system, an estimate of ultimate oil recovery by solution gas-drive depletion usually will be forced to perform the calculations himself. This is despite the quantity of data in the literature. Rarely will either experience or the literature provide results from a reservoir system similar in all important respects to the one under consideration, and calculated results are not so plentiful that satisfactory interpolation procedures can be devised. Performing the calculations, however, is a tedious, time-consuming task unless an electronic computer is available, and, in practice, time and manpower are not always available for this purpose. A quick, simple, consistent method was needed for reducing the uncer- tainty in estimated oil recovery from solution gas-drive reservoirs when only minimum information about the reservoir system is available. PROCEDURE Method of Calculation The usual requisite assumptions were made so that the material balance equation could be used to calculate data for the charts. The following assumptions were made: (1) the reservoir is homogeneous and isotropic; (2) oil recovery is due entirely to solution gas drive and neither a gas cap nor a water drive nor gravity drainage is present; (3) the initial reservoir pressure is the bubble-point pressure of the reservoir fluid; (4) initial total liquid saturation is 100 per cent of pore space; (5) interstitial water saturation remains at the initial value as the reservoir pressure declines from the bubble-point pressure to atmospheric pressure; (6) equilibrium gas saturation is 5 per cent of pore space; and (7) oil and gas saturations are uniformly distributed throughout the reservoir at all times. There are no saturation gradients due to a wellbore, nor is the geometry of the reservoir system considered. The material balance equation was written in the form of a differential equation' which was integrated to determine the change in oil saturation for an assigned pressure drop. Formal integration was not possible, so recourse was made to the Runge-Kutta method' of numerical integration. All computations were performed on IBM equipment. Numerical integration yielded the change in oil saturation within the reservoir as the pressure. declined from the bubble-point pressure to atmospheric pressure. The initial oil saturation minus the change in oil saturation yielded the oil saturation at atmospheric pressure. The oil originally in place was obtained by dividing the initial oil saturation by the initial formation volume factor (differential liberation) while the oil in place at atmospheric pressure was obtained by dividing the final oil saturation by the formation volume factor at atmospheric pressure. U1timate oil recovery, expressed as a percentage of the initial oil in place, was obtained by dividing the difference between the oil initially in place and the oil in place at atmospheric pressure by the oil initially in place and multiplying by 100. PVT Data Charts were based upon 135 solutions to the material balance equation. PVT properties of the reservoir fluids were variables in this equation and had to be known as functions of pressure. The required PVT data might have been obtained from either actual reservoir fluid systems or correlated data. However, correlated PVT data were developed and employed in the calculations for the following reasons: (1) it is doubtful that 135 sets of PVT data could have been obtained for the values of variables investigated in this study, and (2) results of recovery calculations from randomly obtained PVT data could not be correlated as well as rssults from selected PVT data. The PVT data used to develop the

By R. C. McFarlane, T. D. Mueller, J. E. Walstrom

In evaluating uncertainty, experiments are usually performed repeatedly and then conclusions are drawn from the distribution of results. With the advent of high-speed electronic computers, it is possible to perform experiments using mathematical models constructed to simulate complex experiments or operations. Statistical methods are then applied to the results of the simulated experiments. This procedure forms the busis of this paper. Demonstrated is the need for properly accounting for uncertainty in petroleum engineering problems. How uncertainty affects solutions is evaluated in three example illustrations. The method used to evaluate uncertainty in petroleum engineering studies is the Monte Carlo simulation procedure.'-" INTRODUCTION The solution to most technical problems may be derived from interrelationships among several quantities called variables or parameters. There may be only a few variables or several hundred. Interrelationships among parameters may be explicit or implicit, well established or only approximate. Some variables that fully or partially depend on the magnitude of others are called dependent variables. Input variables for most practical problems are not precisely known; there is usually an uncertainty in their value. The degree of uncertainty may vary from one variable to another. Variables that are known accurately are called determinates.' For instance, the gravity of crude obtained from a particular pool may be known precisely, and therefore is a determinate. The degree of precision with which a quantity can be determined increases as data describing the pool are accumulated during the development of the field and the producing life of the pool. The uncertainty of a parameter may result from difficulty in directly and accurately measuring the quantity. This is particularly true of the physical reservoir parameters which, at best, can only be sampled at various points, and which are subject to errors caused by presence of the borehole and borehole fluid or by changes that occur during the transfer of rock and its fluids to laboratory temperature and pressure conditions. Uncertainty may also result in attempting to predict future parameter values. This type of uncertainty is particularly evident in investment analyses involving future costs, prices, sales volumes and product demand. Uncertainty in the solution to investment problems is often called risk, and its study is called risk analysis.' Uncertainty also enters into biological and sociological analyses in which indeterminate factors are often important due to limited control of the experimental material. It is customary, in evaluating uncertainty, to perform repeated experiments and to draw conclusions from the distribution of the results of these experiments. With the advent of the high-speed electronic computer, it is possible to construct mathematical models which simulate complex experiments or operations and to perform the experiments repeatedly, utilizing the models. Statistical methods are then applied to the results of the simulated experiments This method forms the basis of the investigation reported here. PROBABILITY DISTRIBUTIONS FOR VARIABLES The uncertainty in the value of a variable may be indicated by a probabilistic description accomplished by expressing the quantity by a probability distribution. Many recognized probability distributions can be used to describe physical quantities. Recent studies used various types of distributions to describe core analysis data.',' However, for the examples in this paper, the uniform and triangular distributions are believed to reasonably approximate the data used (Fig. 1). The uniform distribution confines the variable between an upper and a lower limit. The variable may lie anywhere between the two limits. This distribution is used when no one range of values for a variable is more probable than any other, but information or intuitive reasoning indicates the variable will lie somewhere between the chosen limits. The triangular distribution is used for a variable when more data are available to indicate a central tendency of distribution. This allows postulating a "most likely" value to the distribution and upper and lower limits. In this case, as for the uniform distribution, the variable is not expected to assume a value less than the lower limit or greater than the upper limit. However, with improved quality of data it can be postulated that the variable will tend to assume a value close to the most likely value, and that there will be a decreasing probability for values away from the most likely value. The area under either of these probability distributions is equal to unity since it is assumed that there is a 100 percent probability that the variable will lie somewhere under the curve. An ordinate erected at any particular value of the variable divides the area under the curve into two parts: the area to the left of the ordinate represents the probability that the value of the variable will be equal to or less than the value of the variable at the position of the ordinate, and vice versa. The probability is zero that the variable will have any specific deterministic value. If two ordinates are drawn for any two values of the variable, the probability that the variables will have a value lying between these ordinates is equal to the area under the curve lying between the ordinates.

By B. L. Landrum, J. Simmons, J. M. Pinson, P. B. Crawford

Patentiometric model data have been obtained to estimate the effect of vertical fractures on the areas swept after breakthrough in water flooding and miscible displacement programs such as gas cycling where the mobility is near one, The data are presented for the case of the fire-spot pattern in which the cemer well is fractured various lengths and orientations, the data indicate that for 10-acre spacing, fractures extetidirrg over 1300 ft in either directior1 from the fractured well may re.srrlt in reductions in sweep efficiencics from 72 to approximately 34 per cent. However. the area swept after break through may be quite largr and only 10 or 12 per cent 1ess than would be obtained if the reservoir were trot fractured. For the specific case when the volume of fluid injected is equivalent to 100 per cent of the pattern vol-unie, the swent area may vary from 80 to 88 per cent, depending on the lenght of the fracture. The former value is that which occurs when the break through or sweep efficiency was orrly 34 per cent and the latter figrrre of 88 per cent is that which is obtained if the reservoir were unfrac-ttm'd. It is pointed out that although the sweep efficiency may he very low in vertically fractured five-spot patterrz.s, the area swept at low water-oil ratios may be only 5 to 10 per cent less than those achieved if the reservoir were unfractured. INTRODUCTION Since the initiation of commercial reservoir fracturing techniques it has been desirable to determine the effect of fractures on the areas swept after breakthrough. Most water flooding or gas cycling projects are continued for substantial periods after the brcakthrough of the injected fluid. Although the sweep efficiency serves as one criterion for rating various flooding patterns. the area swept after breakthrough for various water-oil ratios or percentage wet gas, if cycling. is of perhaps more importance than the sweep efficiency alone. Sweep efficiency data on the vertically fractured five-spot have been presented3. Previous work on the line-drive pattern has shown the effect of vertical fractures on the area swept after breakthrough for the case in which the distance between injection and producing wells divided by the distance between adjacent input wells was equivalent to 1.5 (see lief. 2). The data indicated that for the line-drive pattern it may be desirable to flood or cycle substantially perpendicular to the fractures in order to achieve the greatest recovery for the smallest volume of fluid injected. For this study the center well of a five-spot is assumed as the fractured well. All fractures were assumed to originate at this well and extend into the reservoir for various distances and orientations. All the fractures are straight and are of large permeability compared to the matrix proper. These data are presented to aid the engineer in estimating fractured five-spot pattern performance. ANALOGY The potentiometric model was used in making this study. The model used was 20 20 in. by approximately 1-in. deep. For certain portions of the study one corner of this model was considered to be an injection well and the opposite corner a production well. To simulate vertical fractures a copper sheet was soldered to the wire well and made to conform to the desired length and orientation. In other studies the same model was used except that the four corners of the model might be considered as the corner wells of a five-spot pattern and a fifth well was placed in the center of the model. The well placed in the center of the model was fractured. The total fracture length is L and the well spacing. d. The complimentary fracture angles will be obvious from Figs. 3 and 4. The data obtained on the potentio-metric model assumes the pay to be uniform and homogeneous, the mobility ratio is one, steady-state conditions exist and gravity effects arc neglected. The permeability of the fractures is very great compared to that of the matrix proper. The po-tentiometric model has been used widely both in water flooding and gas cycling projects, and may be used for miscible displacement; how-ever. it is believed that the poten-tiometric model data are more properly applicable to gas cycling than water flooding because the model as-

By R. F. Jones, N. E. Truitt, M. Prats

In most published discussions and theories of in situ combustion, the combustion fronts are assumed to be vertical. However, evidence from field tests leaves no doubt that combustion fronts often advance more rapidly along the top than near the bottom of a formation as a result of difference in density between injected air and formation liquids. The approximation proposed in this paper to determine the movement of the resultant tilted combustion surfaces states that the vertical rate of movement of combustion surfaces is proportional to the horizontal oxygen flux. Where self-ignition is possible, the proposed method demands that a secondary combustion surface exist around production wells which produce some oxygen. These secondary combustion surfaces may be formed long before the primary combustion surface can advance from injection to production wells. Heat liberated near production wells at these secondary combustion surfaces can contribute to an early increase in production rate. Results indicate that significant oil recoveries cannot be obtained from the usual flood patterns (five-spots, seven-spots, etc.) without producing large volumes of unused oxygen. Ideally, to increase oxygen- consumption efficiency, well patterns should allow oil production from a first line of production wells and gas production from more distant lines of producers. However, it may be desirable to produce some gas at all wells to support (and benefit from) active secondary combustion surfaces. Results indicate that the well spacing through which combustion can be advanced is larger than that predicted by other methods. A large number of production wells may still be desirable to take quick advantage of gravity drainage. From a comparison with results at South Belridge field, California, it appears that this method adequately describes oxygen concentration and temperature histories and combustion-front shapes. However, this method does not accurately locate the most advanced point of the combustion surface. There is some field evidence to substantiate the actual presence of secondary combustion surfaces at South Belridge. Use of the proposed method appears warranted at this time when lay-over of the combustion surface can be anticipated. INTRODUCTION The assumption of vertical combustion fronts has been embodied in all previous publicationsl-6 which use the movement of combustion fronts away from injection wells to determine the temperature and fluid distributions in the reservoir. The only paper concerned with a mathematical model of the combustion process in which a nonvertical combustion front is used was written by Gottfried.11 Actually, nonvertical combustion fronts have been observed in most in situ combustion field tests for which adequate data are available.7, 8, 6 In practice, the typical vertical extent of the burned zone decreases with distance from the injection well, and this burned zone is at or near the top of the sand body. In some cases, such as at South Belridge field near Taft, Calif., the combustion surface is almost horizontal over a very large area.8 Thus, for some years an obvious and serious gap has existed between theory (vertical fronts) and practice (tilted fronts). This is indicated in Fig. 1. Tilted combustion fronts such as observed at South Belridge8 sometimes result from the natural tendency of injected gases to rise to the top of an oil sand. An initially present gas cap, even if thin, would make apparent the effect of gravity segregation from the first, particularly since formations considered for in situ combustion field operations are generally selected for their

Jan 1, 1969

By R. E. Cook

This paper presents results of a study to determine to what extent errors in estimated free gas saturation affect the results of static pressure calculations from build-up curves in two-phase systems. Use is made of the method of pressure build-up curve analysis developed by Miller, Dyes and Hutchinson1 , with modifications for two-phase flow advanced by Perrine2. Cornhined fluid compressibility of gas-oil mixtures is included, with special attention given to the discontinuity that occurs in fluid compressibility at the saturation pressure and its effect on static pressure calculations from build-up curves. The sensitivity of static pressure calculations to errors in gas saturation is considered over wide ranges of formation and fluid properties and reservoir pressure conditions. It is concluded from this study that for analyses of two-phase pressure build-up curves (I) the use of under-saturated crude compressibility data, even when producing GOR's are at or near solution values, will result in very significant negative errors in calculated static pressrrre, and (2) except for large values of buildup curve .slope and/or very low crude gravities, rather large errors in calculated gas saturation can he tolerated without .significant effect on calculated static pressures. The conc1usions arc qualified as referring only to the effects of gas saturation and excludc possible errors in the nuumerous other factors that enter into the ancllysis of pressure build-up curves. INTRODUCTION As part of the analysis of basic data for several major Eocene reservoirs in the Bolivar Coastal field (Western Venezuela), it was necessary to analyze a large number of pressure build-up curves for oil permeability and static reservoir pressure. From performance data it became apparent that gravity segregation is very active in these reservoirs. As a result of this performance it was recognized that the free gas saturation within the oil leg of the reservoirs would be small but would be subject to certain errors in calculation due to the inaccuracy of relative permeability data in the low gas saturation region. Since all of the various methods of build-up curve analysis were originally derived for sin- gle-phase, constant compressibility systems, and since the presence of a small but rather indeterminate gas saturation would greatly affect the reservoir fluid compressibility, a study was undertaken to investigate the effect of gas saturation on the results of two-phase pressure build-up analyses. TWO-PHASE BUILD-UP CURVE ANALYSIS The method of analysis used in this investigation is basically that originally developed by Miller, Dyes and Hutchinson1, with modifications for two-phase flow advanced by Perrine2. Miller, et al, showed that the original equation for single-phase permeability away from the well bore, held for two-phase flow except that the calculated permeability is the effective permeability to oil away from the wellbore. This finding was acknowledged by Per-rine, who also presented the equation for effective gas permeability* as Since Eqs. 1 and 2 contain no terms influenced by fluid compressibility, no distinct problem should develop in calculating permeability from two-phase build-up curves. In regard to the calculation of static reservoir pressure from build-up curves in single-phase systems, the modified Miller, et al, method as presented by Perrine is briefly as follows. 1. From build-up curve data, PVT data for the reservoir fluid and the well spacing and sand thickness, calculate the climensionless shut-in time corresponding to the highest value of actual shut-in time (t, hrs) on the straight-line portion of a semi-log build-up curve. Dimenrsionless shut-in time is related to actual shut-in time by the relation, t- 0.000264 kt 2. From a graph of dimension less shut-in time vs dimensionless pressure difference (Fig. 1 1 ) find the value of dimensionless pressure difference corresponding to the value of 7 calculated under (1) above. Dimen-

By E. L. Anders

The Mile Six pool is located on the La Brea-Parinas Cullcession of International Petroleum Co., Ltd., in northwestern Peru on the west coast of South America. The reservoir pressure in this pool has been maintained within 200 psi of its initial value throughout its history, and gravity drainage has played an important role in the production behavior. It has now produced 95 per cent of its estimated ultimate recovery. It is estimated that this interesting oil pool will ultimately produce 67 per cent of the initial oil in place and that the resulting residual oil saturation may be as low as 19 per cent of the pore volume (29 per cent of the hydrocarbon pore volume). An evaluation of reservoir rock and fluid characteristics and ultimate oil recovery is presented. INTRODUCTION This study of Mile Six pool was made to evaluate its performance according to latest available information. The production performance of this pool has been discussed in various articles in the past. and the reported behavior has been used as an example for application of computation procedures for gravity drainage depletion' and as an illustration of field behavior under gravity drainage or expanding gas cap drive. There have been wide variations in reported values of initial oil in place, reservoir oil volume factor. connate-water sauration, volume of effective sand, and ultimate recovery because of the paucity of reliable basic data. These various factors have been determined as accurately as practicable with the latest available information, and this evaluation is presented herein. The production history of Mile Six is an excellent example of gravity drainage depleion with effective pressure maintenance by gar injection. GENERAL Mile Six pool was discovered by cable-tool drilling in November. 1927. when well 1996 was completed in the Parinas sand. After slow development with cable tools and sporadic production. the pool was opened to continuous pro(iuction in November. 1933. and develpment was completed with rotary rigs. Pressure maintenance was started in December, 1933. by returning gas to upstructure wells. Most of the development was 'completed by 1937, but some additional wells were drilled in the period 1939-1947. and several old wells were deepened. A total of 46 oil and gas wells and 4 dry holes were drilled on approximately 7-acre spacing. Of the producers. 21 are now flowing. 2 are pumping. 44 are gas input wells. 3 are abandoned. I is a gas well shut in. and 15 are shut in because of non-commercial production or high gas-oil ratio. The locations of all wells are shown on the map of Fig. I. Total oil production on Dec. 31. 1952, was 30,867,373 bbl: cumulative gas production was 22,023,777 Mcf; and 26,410,946 Mcf of gas had been returned to the reservoir. These figures do not include oil and gas lost ill a blowout in January. 1940. GEOLOGICAL DESCRIPTION Mile Six pool is located on the northern end of a structural spur projecting from the La Brea-Negritos uplift.' The spur is probably a reflection of a basement structure. It plunges gently to the north, i broken into a complex series of fault blocks. and contain.; the Verdun Alto. Section Sixteen. and Mile Six pools. The Parinas handstone (lower Eocene!. which is the producing formation in Mile Six. occurs at an average depth of 2,200 ft in the pool and dips north and east at from 15' to 20". The pool covers an area of approximately 350 acres. Mile Six is downfaulted about 600 ft from Section Sixteen pool to the soutb. and a major fault forms its western boundarv. The north and east boundaries are formed by the intersection of the sand top with the water-oil contact which occurs at approximately 2.440 ft subsea. An original gas-oil contact probably existed at about 1.875 ft subsea. Fig. 1 presents the latest structural interpretation of the pool. and Fig. 2 is an isopach map showing thickness of the total Parinas formation above the original water-oil contact. The heavy lints of Fig. 1 are contours on the sand top, and the fine lines are contours on the fault planes. This type of straight-line structural map was developed 1)) International's geologists to reflect structural conditions where the bedding planes dip and have no curvature. 'The La Brea-Parinas Concession is highly faulted by normal fault.. The beds are flat wherever exposed. The Parinas formation is approximately 635 ft thick. and it is etimated that 62 percent II the formation is effective sand. The original oil zone was about 565 ft thick. Fig. 211 presents an electric log showing typical Parinas sand development ill Mile Six pool. The Parinas band in Mile Six i. a well-sorted. medium- to-coal-se-grained. cross-bedded sand with minor lenses of shale and small lenses and pockets of pebble conglomerate. The sand grains are subangular to rounded and consist chiefly of quartz with feldspars. biotite hornblende, and augite as accessory minerals. Because of faulting of the Parinas- formation to the east and north of the pool. there is probably l possibility of a significant, natural water drive in Mile six. The faults within the pool. as indicated in Figs. 1 and 2. are of smaller disI,laceInent and seem to act only , partial barriers to fluid movement within the reservoir. RESERVOIR CHARACTERISTICS Core analysis data are available from five wells. The data were obtained from three wells (Nos. 3401. 3586 and 3719) at the time of their completion and from well 1996 when the original liner was sidetracked and the well was deepened ill 1946. Data from well 2779 were obtained in 1943 from old cores taken when the well was deepened in 1934. From these core analyses. the average porosity was estimated to be 22.6 per cent. and the average permeability to dry ail. was estimated to be 780 and Measured productility indices varied from 3.1 to 71.4 B/D per psi differential. Specific productivity indices varied approximately from 0.1 to 0.3 B./D per psi per ft of sand.

Jan 1, 1953

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By J. C. Martin

Equations for three-phase, three-dimensional, compressible flow (including capillarity) are reduced to two-dimensional relations by a partial integration. This reduction allows three-dimensional flow problems to be treated with mathematics for only two spatial dimensions. The results can be used to formulate flow equations for two-dimensional reservoir simulators in which the effects of capillarity and fluid segregation in the third dimension are represented. Such reservoir simulators would retain many of the aduantages of two-dimensional simulators while simulating three-dimensional effects. The principal restriction of the method is that the thickness of the reservoir should be small, compared to the distance across the reservoir. INTRODUCTION In recent years, computers have been used to calculate the performances of many reservoirs. Most of the detailed calculations, however, are based on finite difference solutions of the flow equations, and present day computers are seldom able to handle a sufficient number of cells to produce entirely satisfactory solutions, even for reservoirs represented by two-dimensional arrays of cells. The situation becomes much worse when one wishes to approximate the reservoir by a three-dimensional array. A great saving in computation or a more detailed solution can be obtained for many reservoirs by using the partial integration of the equations of flow presented in this paper. The integration reduces the three-dimensional equations to two-dimensional relations; and, for studies of two-dimensional flows in vertical cross-sections, the equations are reduced to one-dimensional relations. Most reservoir performance calculations currently are based on one- or two-dimensional flow re1ations.l-3 In some cases flow in the third dimension is approximated by assuming a particular type of vertical saturation distribution, such as gravity segregation.l.2. The relations developed in this paper approach those for segregated flow as the capillary pressures approach zero, and they approach those for uniformly distributed saturations as the capillary pressures are increased. For this analysis, the ratio of the reservoir's thickness to the maximum distance across it must be small. The capillary pressures between the oil and water should also be small compared to the maximum pressure difference in the reservoir. This requirement is met by most reservoirs. It is assumed that the capillary-pressure curves are well defined, whether or not hysteresis effects are included. Also, the reservoir must have sufficient vertical permeability to allow the fluids to segregate. The results presented here provide a firm theoretical foundation to Coats' et aL5 assumption of vertical equilibrium and extend the relations to three-phase flow. Coats' assumption of vertical equilibrium, which he verified by calculations and experiment, is developed here mathematically from basic flow equations. DISCUSSION SATURATION AND PRESSURE DISTRIBUTIONS Appendix A presents a mathematical analysis of fluid flow in reservoirs where the ratio of thickness to maximum distance across the reservoir is small. The results indicate: (1) that the fluids along any line perpendicular to such a reservoir's upper surface are in static capillary equilibrium (vertical equilibrium); (2) that, to a first approximation, the fluid pressures and properties are functions of only areal position in the reservoir and time; and (3) that hydrostatic pressure gradients exist along any line perpendicular to the reservoir's upper surface. The results might be expected after studying several physical considerations. First, no flow is allowed normal to the upper and lower reservoir boundaries, which are relatively close together.

Jan 1, 1969

By P. B. Crawford, W. L. Huskey

A study has been made to determine the effects of random vertical fractures existing in a reservoir matrix on the effective matrix properties and producing characteristics of a well. The study indicated that short vertical fractures may distort the isopotentials by 1 percent in the vicinity of the fracture. A correlation was developed which showed the effect of fracture density on the producing capacity of a well. The correlation indicated that fracture shape had little effect on producing capacity, but that total fracture length was closely correlated with the producing capacity. Effective values of reservoir permeability were found to correlate closely with fracture density when measured parallel to the direction of flow. INTRODUCTION In reviewing the petroleum literature, it is found that in most analytical solutions of hydrocarbon reservoir behavior the assumption is made that the rock is uniform throughout although core analysis and well logging surveys indicate that substantially all reservoirs are heterogeneous and some are fractured. Only in recent years have investigators studied the effects of reservoir hetrogeneity. Russell1 discussed the effect of fractures in reservoirs. The literature indicates that fractures greatly affect the exploitation of man reservoirs. This is described by Elkins and Skov2 and Little-field et al. 3 Several types of reservoir heterogeneity have been studied by unsteady-state methods. In 1959, Landrum et al. 4 showed that transient phenomena in hydrocarbon reservoirs could be studied by thermal models. pickering5 studied the flow of heat in stratified linear reservoirs consisting of plates of different metals separated by sheets of rubber which were suddenly exposed to a low-temperature source at one end. Cotman6 studied both laterally and vertically heterogeneous thermal model reservoirs with special reference to cross flow. Another type of reservoir heterogeneity involves the presence of highly permeable, circular lenses around or near the wellbore which was studied by Miesch7 using a thermal and a steady-state model. Many studies have been made of the effects of induced hydraulic fractures on well performance. 8-10 Work regarding vertical fractures that exist initially in the reservoir was reported by Givens, et al.11 Since there appears to be little published work on the importance and influence on reservoir performance of fractures that exist out in the rock matrix, this study was conducted to determine the effects of vertical fractures on pressure distribution, producing capacity and the effective permeability of the reservoir rock. The analogy between the flow of electrical current through a conducting body and the flow of fluid in a porous medium provides a convenient method of studying various aspects of petroleum reservoirs. For electrical radial steady-state flow,12 For steady-state fluid flow, With the above analogy, it is only necessary to establish the ratio of the values of the electrical units in terms of the corresponding flow units to make the analogy quantitative. The ratios are shown below.

By C. F. Weinaug, R. W. Farley, J. F. Wolfe

A rapid, accurate method for predicting the dew points of gas condensate systems and their subsequent normal and retrograde phase behavior with pressure decline has been developed. The method predicts dewpoint pressure, volume percent liquid, condensate content of produced gas, and the liquid and gas compositions at pressures below the dewpoint pressure. Use of a modified form of the BWR equation of state, which was established by regression calculations on experimental data from several gas condensates, provides depletion behavior predictions that are significantly improved over predictions using K-values based on convergence pressures and Standing's correlation for density. INTRODUCTION Prediction of complex hydrocarbon-mixture thermodynamics is fundamental to reservoir and gas process engineering. In particular, accurate descriptions of phase behavior are needed to plan production operations for gas-condensate reservoirs. Although analysis of reservoir behavior during producing and cycling operations can be made for simple reservoir models, extensions of these studies to more complex reservoir systems are still difficult, both with regard to numerical techniques and the phase behavior of the fluids. This report presents a mathematical method for handling the complex phase relationships encountered in retrograde condensates; it represents an advance in reservoir behavior technology in that it will aid the development of closed-form mathematical representations of total production systems. A typical phase diagram for a complex hydrocarbon mixture is shown in Fig. 1. When the reservoir temperature T, is higher than the critical temperature Tc, and lower than the critical-condensation temperature Tct' the mixture is a gas condensate. Decreasing the pressure on a condensate while maintaining nearly constant temperature will result in two phases forming at the dew point. Further decrease in pressure will cause additional liquids to be formed. This phenomenon is called retrograde condensation. As pressure continues to decrease, normal phase behavior regions are encountered (liquid vaporizes as the pressure decreases) until the low-pressure dew point is encountered. When a gas-condensate reservoir is produced, gas is withdrawn and the volume of fluids in the reservoir remains nearly constant as the pressure declines. Heavier components tend to accumulate in the liquid phase as a result of retrograde condensation, while lighter components tend to concentrate in the vapor phase. Thus, as the gas is produced the composition of the gas and the condensate

Jan 1, 1970

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas'-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller's results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By R. C. Earlougher, F. G. Miller, T. D. Mueller, H. J. Ramey

There are many studies of flow in radial systems that can be used to interpret unsteady rerervoir flow problems. Although solutions for systems of infinite extent can be used to generate solutions fur finite ow systems by superposition, application is tedious. In this paper a step is made toward simplifying calcu1ations of such solutions for finite flow systems. Superposition is used to produce a tabulu-lion of the dimension1e.s.s pressure drop function at several locations within a bounded square that has a well at its center. The square system provides a useful building block that may be used to generate flow behavior for any rectangular shape whose sides are in integral ratios. Values of the tabulated dimensionless pressure drop function are simply added to obtain the dimensionless pressure drop function for the desired rectangulm system. The rectangular system may contain any number of wells producing at any rates. Furthermore, the outer boundaries of the rectangular system may be closed (no-flow) or they may he at constant pressure. Mixed conditions also may be conyidered. Tables of the dimensionless pressure drop function for the square system are presented and various applications of the technique are illustrated. Introduction In 1945, van Everdingen and Hurst' published solutions for the problem of water influx into a cylindrical reservoir. Since this problem is mathematically identical with the depletion of a cylindrical reservoir with a well at the origin, the van Everdingen-Hurst solution may be used to study the depletion problem. In their analysis, they assumed that the fluid had a small, constant compressibility such that flow was governed by the diffusivity equation For a constant production rate q starting at time zero, van Everdingen and Hurst showed that the unsteady pressurc distribution for both finite and infinite systems could be expressed in terms of a dimensionless pressure 'Sabulations of the dimensionless pressure drop for a unit value of r, were provided by van Everdingen and Hurst,' and later by Chatas.' Others also presented values in graphical or tabular form." If the radius of the well becomes vanishingly small, r,+ O, the line source solution may be used for Eq. 2 when infinite systems are considered. Eqs. 5 and 6 are excellent approximations for Eq. 2 under certain conditions: In 1954, Matthews, Brons and HazebroekV emonstrated that solutions such as Eq. 5 can be superposed to generate the behavior of bounded geometric shapes; i.e., the behavior of a bounded single-well system can be calculated by adding together the pressure disturbances caused by the appropriate array of an infinite number of wells producing from an infinite system. These wells are referred to as image wells. Matthews, Brons and Hazebroek considered systems containing a single well producing at a constant rate. This superposition can be represented analytically as

Jan 1, 1969

By D. M. Beeson, G. D. Ortloff

Water-propelled banks of carbon dioxide recovered both high- and low-viscosity crude oil substantially in excess of that recovered by water flood in linear flow model experiments. The increase in oil recovery is explained by oil swelling and viscosity reduction. Additional oil was recovered by solution gas drive at the end of the carbon dioxide floods when the model systems were depres-sured. This study was an investigation of the use of a water-driven carbon dioxide bank to recover crude oil. Techniques involving the injection of carbon dioxide have aroused con-siderable interest in recent years among methods considered for increasing oil recoveries. A number of investigators have studied the process in the laboratory1-' and in the field.7,5 Laboratory studies referenced indicate that carbon dioxide in conjunction with water flooding recovers more oil than does water flooding alone. These re- ports diverge widely, however, in stating magnitude of improvement. We believed that a water-driven carbon dioxide bank would increase recovery of crude oils by two mechanisms: viscosity reduction and oil swelling. Due to reduction in the oil viscosity caused by solution of carbon dioxide, displacement efficiency should be improved. In addition, the residual oil should be highly swollen by the carbon dioxide held in solution. Thus, the final residual oil saturation, in terms of stock-tank oil, would be lower than in the absence of carbon dioxide. To test this hypothesis, experi-

By W. T. Weller

In Part I, a study of pressure build-up curves calculated for conditions under which both oil and gas flow led to the conclusion that the presence of a dispersed free gas phase in an oil reservoir must be taken into consideration to estimate accurately average reservoir pressure and permeability from build-up curves. Familiar methods based on the assumption of no free gas can be extended to the two-phase case by using total compressibility and mobility in place of single-fluid compressibility and mobility. These methods give correct values for average pressure and permeability when gas saturation is small. Errors become larger as the gas saturation increases. However, for the use to which the results will be put, the methods are satisfactory for reservoir engineering purposes. An improved method of calculating the performance of depletion-type reservoirs is presented in Part 2. Because the mathematical relationships describing simultaneous flow of oil and gas are quite involved. simplifying assumptions are made to provide means of obtaining approximate solutions of reasonable accuracy. One such approximate method now in use is the constant-GOR solution. It involves the assumption that at any instant, the ratio of total gas flow rate (both free and disolved) to oil flow rate is the same at all points in the reservoir. The approach is not applicable unless the free gar saturation in the reservoir is everywhere greater than the critical gas saturation. This paper presents a modification which, by avoiding the constant-GOR assumption, makes the method applicable to all reservoir conditions, and so far appears to be more accurate than the constant-GOR solution and to he comparable in required compuctation time. PART I— BUILD-UP CURVE ANALYSIS INTRODUCTION Pressure build-up characteristics of shut-in wells have been used for many years by engineers to evaluate average reservoir pressure, effective permeability thickness of the pay section, effectiveness of well completion (skin effect) and reservoir size. A number of methods of analysis have appeared from time to time."" Without exception, these methods are based on the assumption that the reservoir contains but one fluid of constant small compressibility and constant mobility. It has been suggesteda" hat these single-fluid methods may be applied to data from reservoirs containing both oil and gas by substitution of some effective total properties of the multiphase system in place of the corresponding single-phase properties. The present investigation was undertaken primarily to evaluate this approach. METHOD A number of theoretical build-up curves were calculated for conditions of two-phase flow, under the assumption of certain reservoir and fluid properties, and were analyzed by single-fluid methods with appropriate total compressibility and total mobility values for the corresponding single-fluid properties. Results of the analyses were compared with the assumed conditions. The theoretical build-up curves were completed by procedures similar to those of West, Garvin and Sheldon." Since these calculations require a considerable amount of computer time, an attempt was made to derive an approximate calculation method. The attempt was unsuccessful for calculating build-up curves, but the effort did result in a new approximate method of calculating the performance of solution gas drive reservoirs, which appears to be an improvement over the constant-GOR method" used previously (see Part 2). The West, Garvin and Sheldon calculations involved the following assumptions: (1) the reservoir is circular and completely bounded, with a completely penetrating well at its center; (2) the porous medium is uniform and iso-tropic, with a constant water saturation at all points; (3) gravity effects can be neglected; (4) compressibility of rock and water can be neglected; (5) the composition and equilibrium are constant for oil and gas; (6) the same pressure exists in both the oil phase and the gas phase; and (7) no afterproduction occurs, i.e., the well is shut in at the sand face for build-up. These assumptions make it possible to describe two-phase flow of oil and gas by the partial differential equations:

Jan 1, 1967