Search Documents

By W. J. West, T. D. Mueller, J. E. Warren

Two methods are used to obtain the relatiorlship between ratio of gas to oil permeability and gas saturation. One method is to calculate a curve from field measurements of reservoir behavior: the other is to measure the relationship on a core .sample in the laboratory. The curves resulting from these TWO methods arc often comparer1 and seldom found to agree. The most important reason for the nonagreement is often overlooked. This reason is that the performance k9/k0 curve is calculated with the assumption of uniform pressure and saturution throughout the reservoir. To show the differerence between performance type kg/k0 curves and laboratory type kg/k0 curves. two performance type kg/k0 curves for two different production rates were calculaled from the production perforrnance of a hypothetical reservoir for which a laboratory type curve was assumed. Production performance of the hypothetical reservoir bras calcu- lated by use of international Business Machine's 701 Computer. The perfortnance type kg/k0 curves differed from the laboratory curve for the hypothetical reservoir. The performance type crrrve for the lowest production rate most nearly coincided with the laboratory curve. INTRODUCTION Predictions of oil reservoir performance are often made with the aid of complex numerical calculations. An important quantity required in these calculations is the ratio of gas permeability to oil permeability (kg/k0 )at a particular gas saturation. Two methods arc used to obtain a relationship between this ratio and gas saturation. The curves resulting from these two methods are called performance kg/k0 curves and laboratory kg/k0 curves. The performance curve is obtained from field measurements of reservoir behavior. Particular kg/k0 values are calculated from the producing gas-oil ratio, the cumulative reservoir withdrawals, the original oil in place, the average reservoir pressure, and the PVT properties of the reservoir fluids. The particular kg/k0 values are then plotted against the average rcservoir gas saturations for the entire history period. In order to use this derived kg/k0 relationship for performance predictions, a curve through these points is extrapolated to higher gas saturations. The laboratory kg/k0 curve is obtained from measurements made on core samples. Saturation is maintained uniform throughout the core during these measurements. The kg/k0 curves obtained with these two independent methods are often compared and seldom are found to be the same. Various reasons have been offered for the lack of correlation. One is that insufficient core? were obtained to derive a curve representative of the reservoir rock. Another is that laboratory methods have failed to reproduce reservoir conditions. The most important reason for the noncorrelation has been often overlooked because of the difficulty in evaluating its effect. This reason is that the performance or field kg/k0 curve is calculated with the assumption of uniform pressure and saturation throughout the reservoir. In any

Jan 1, 1956

By C. R. Dodson, D. Goodwill, E. H. Mayer

The paper summarizes the historical background of the pressure-volume-temperature analyses of reservoir fluids, the errors involved in both the sampling and testing of reservoir fluids, the type of information required of a PVT determination, and the field conditions that limit the application of any one analysis. Particular emphasis is placed on the necessity for approximating as closely as possible the liberation sequence occurring in the producing formation, flow string, and surface separators. A combined differential and flash or "composite" liberation is suggested as the best means of approximating this liberation sequence. HISTORICAL BACKGROUND Petroleum reservoir engineering commonly is considered to be one of the newest fields of petroleum science, yet much work of a theoretical and intuitive nature was done many years before the modern techniques of reservoir analyses were developed. The period 1910-1924 saw considerable work in the field of reservoir behavior done by the U. S. Bureau of Mines. This work:.' although entirely theoretical, pointed out the importance of gas in the recovery of oil from the reservoir. The statement by J. 0. Lewis that 20 per cent or less of the oil originally in place in a pool was recovered from the ground under uncurtailed production conditions caused some operators to re-examine and modify their production methods. However, though these works showed methods of production that have subsequently resulted in greater oil recoveries, they in general went unnoticed, even though there was a fear that the nation's petroleum resources were being exhausted. The period 1924-1933 saw the industry take considerably more interest in reservoir behavior because of an unfounded fear of its inability to replace oil reserves, the possibility of government regulation, and the energy of one man, Henry L. Doherty. Doherty aroused heated discussions in the industry concerning conservation. To prove or disprove his theories the first experimental work on the reservoir behavior of petroleum was undertaken. The papers of Dow and Calkin,' Beecher and Parkhurst," and Mills and Heitheckere are classics although their experimental procedures were crude. These papers proved that the properties of petroleum in the reservoir are quite dissimilar to the properties measured at the well head. Gas dissolved in the oil phase was recognized as having considerable importance as had been predicted by the earlier theorists. Slightly later Miller and his co-workers assembled the first text correlating all the known data on the reservoir behavior of petroleum and used, these data to show the economic value of conservation. Quantitative application of pressure-volume-temperature analyses of reservoir fluids was given in the paper of Coleman, Wilde, and Moore"; this work showed that with sufficient laboratory and production data prediction of reservoir behavior was possible. The refined sampling and labora-

Jan 1, 1953

By I. Javandel, P. A. Witherspoon

The finite element method was originally developed in the aircraft industry to handle problems of stress distribution in complex airframe configurations. This paper describes how the method can be extended to problems of transient flow in porous media. In this approach, the continuum is replaced by a system of finite elements. By employing the variational principle, one can obtain time dependent solutions for the potential at every point in the system by minimizing a potential energy functional. The theory of the method is reviewed. To demonstrate its validity, nonsteady-state results obtained by the finite clement method are compared with those of typical boundary value problems for which rigorous analytical solutions are available. To demonstrate the power of this approach, solutions for the more complex problem of transient flow in layered systems with crossflow are also presented. The generality of this approach with respect to arbitrary boundary conditions and changes in rock properties provides a new method of handling problems of fluid flow in complex systems. INTRODUCTION Problems of transient flow in porous media often can be handled by the methods of analytical mathematics as long as the geometry or properties of the flow system do not become too complex. When the analytical approach becomes intractable, it is customary to resort to numerical methods, and a great variety of problems have been handled in this manner. One such method relies on the finite difference approach wherein the system is divided into a network of elements, and a finite difference equation for the flow into and out of each element is developed. The solution of the resulting set of equations usually requires a high speed computer. When heterogeneous systems of arbitrary geometry must be considered, however, this approach is sometimes difficult to apply and may require large amounts of computer time. The finite element method is a new approach that avoids these difficulties. It was developed originally in the aircraft industry to provide a refined solution for stress distributions in extremely complex airframe configurations. 27 Clough has recently reviewed the application of the finite element method in the field of structural mechanics.6 The technique has been applied successfully in the stress analysis of many complex structures.l, 27, 28 Recognition that this procedure can be interpreted in terms of variational procedures involving minimizing a potential energy functional7 leads naturally to its extension to other boundary value problems. In the field of heat flow, there recently have been introduced several approximate methods of solution that are based on variational principles.2-4, 17 By employing the variational principle in conjunction with the finite element idealization, a powerful solution technique is now available for determining the potential distribution within complex bodies of arbitrary geometry. In the finite element approximation of solids, the continuum is replaced by a system of elements. An approximate solution for the potential field within each element is assumed, and flux equilibrium equations are developed at a discrete number of points within the network of finite elements. For the case of steady-state heat flow, the technique is completely described by Zienkiewicz and Cheung.33 Since the flow of fluids in porous media is analogous to the flow of heat, Zienkiewicz et al. have employed the finite element method in obtaining steady-state solutions to heterogeneous and anisotropic seepage problems. 34 Taylor and Brown have used this method to investigate steady-state flow problems involving a free surface.25 The work of Gurtin has been instrumental in laying the groundwork for the application of finite element methods to linear initial-value problems.12 As a result, Wilson and Nickell have recently

Jan 1, 1969

By E. H. Timmerman, A. F. van Everdingen, J. J. McMahon

The prevent paper contains a method which combines the material balance equation' with the water influx equation' to obtain reliable values for the active oil originally in place and a quantitative evaluation of the cumulative water influx. The method is illustrated by an application to a reservoir without original gas cap. In the absence of an original gas cap. results may be obtained using only field production and pressure data, PVT analyses. and a minimum of subsurface information. CHARACTERISTICS OF RESERVOIR AND RESERVOIR FLUIDS Production in the field under review is obtained from the top of the Wilcox formation of Eocene age, at a depth of approximately 8,100 ft subsea. The structural map, Fig. 1, shows that the accumulation, half elliptical in shape with the long axis in east-west direction, is trapped to the north and on the upthrown side of a normal fault. A number of east-west faults with additional minor faulting in random directions are found in the general area, and seismic information and production performance indicate in this particular case the existence of a second fault a short distance downdip, with a strike more or less parallel to the fault controlling the accumulation. It is estimated that 1,830 acres were originally underlain by oil and that the maximum thickness of the original oil column, about 37 ft, was only slightly greater than the maximum thickness of the oil-bearing sand. Maximum net sand thickness was estimated at approximately 26 ft. From logging and subsurface information the gross sand volume has been placed at 37.400 acre-ft. which reduces to about 27,500 acre-ft after probable non-productive intervals have been deducted. The following table presents laboratory analyse; data for the Wilcox sand cores. Average Average Samples Considered Permeability Porosity All Samples 236 md ..19.9percent Samples with k>l md 256 md 20.5 per cent Samples with k>8 md 275 md 20.9 per cent Interstitial water was estimated at 15 per cent. The initial gas-oil ratio was approximately 900 cu ft/bbl. and the initial shrinkage factor was determined by the laboratory to

Jan 1, 1953

By W. J. Joslin

The frontal-advance equation can determine how the fluid withdrawal rate and subsurface operating pressure influence oil recovery from pressure-maintained reservoirs having characteristics favorable for efficient vertical segregation. For effective vertical segregation in formations having low dip angles, the requirements are: thick, massive sands, low viscosity oil and high vertical permeability with few discontinuous shale members. This approach has been successfully applied to the LL-370, which is a massive Eocene sandstone reservoir containing crude with an average gravity of 26" APT, a dip of only 3" and vertical permeability of ½ darcy. The average gas cap saturation was calculated for various rates of vertical gas-oil contact movement corresponding to a wide range of reservoir producing rates. This method also evaluated the effect of different pressure maintenance levels on oil recovery. Applicability of this procedure for calculating the present oil saturation in the LL-370 gas cap has been substantiated from field performance, where the oil saturation was computed by comparing oil recovery with the gas cap volume vacated. INTRODUCTION During the natural depletion history of the LL-370 reservoir, located in Lake Maracaibo, Venezuela, excellent vertical segregation resulted in the formation of a thin gas zone over a substantial portion of the oil-productive area. Subsequently, gas injection forced this thin gas-saturated zone slowly downward, reducing the oil saturation in the gas cap to a low value. The oil recovery has been raised considerably above earlier predictions where vertical segregation effects were under-estimated. In order to determine the maximum efficient producing rate and optimum operating pressure for this reservoir, a vertical gas frontal-advance analysis has been made through application of the fractional flow equation. LL-370 RESERVOIR CHARACTERISTICS The reservoir consists of a truncated monocline of massive. highly-permeable Eocene sandstone, called the B-b formation, covering an area of more than 13,500 acres. The original oil in place was 2.17 billion STB. Internal non-sealing faults subdivide the structure into four blocks called the B-6-x.6, 10, 11 and 13, which are effectively bounded by major faulting on the flanks, an unconformity updip to the west, and an immobile aquifer to the southeast. The producing formation occurs at an average depth of 5,250 ft subsea and has a closure of more than 1,500 ft. A typical producing section has a net oil sand thickness of 170 ft. Scattered throughout the section are a few randomly-positioned shale lenses which increase in thickness and areal extent downdip. The reservoir rock and fluid properties vary with depth. The permeability and produced oil gravity are 1600 md and 28" API at the top, compared to 300 md and 18° API at the water-oil contact. The productivity of updip wells is, therefore, much higher than the reservoir average. The reservoir oil was originally saturated at the crest, grading to a slight degree of undersaturation at the water-oil contact. There was no original gas cap; and material balance and the field performance indicate no water influx. Except for the northern B-6-x.6 block, the fluid and rock properties have already been published.' Average LL-370 reservoir properties are shown in Table 1. RESERVOIR PERFORMANCE UNDER NATURAL DEPLETION Fig. 1 shows production performance from 1939 to date. Early development was gradual, with 30 wells pro-

Jan 1, 1965

By E. L. Claridge, K. S. Lee

Areal sweep efficiency of oil displacement by enhanced-viscosity water exhibiting pseudoplastic behavior was measured in a Hele-Shaw model representing one-quarter of a five-spot pattern. The pseudoplasticity of polymer solutions and the velocity distribution in the five-spot pattern produced a condition under which the mobility ratio between the displacing and the displaced fluid could not be assigned a single value. Instead, the movement of the displacement front is governed by local mobility ratios which are also time dependent. The areal sweep at breakthrough with polymer solutions was poorer than the sweep obtained with Newtonian fluids of comparable viscosity. However, the areal sweep at 1 PV throughput was greatly improved as compared to flood water without polymer. It was also demonstrated that, even after the oil-cut had declined to a low value during a regular waterflood, switching to polymer flood efficiently swept out the oil remaining in the model. INTRODUCTION The behavior of fluid displacements in isotropic porous media for various patterns of injection and production wells has been extensively investigated.1-7 These investigations all concerned Newtonian fluids, i.e., the viscosity of each fluid was constant regardless of flow rate. The generally unfavorable influence on areal sweep efficiency of higher mobility of the displacing fluid as compared to the mobility of the displaced fluid has been established for both miscible and immiscible fluids. The principle was also established that a close correspondence exists between miscible and immiscible flood front behavior,2, 3 although oil recovery in a waterflood at unfavorable mobility ratio may be less than that observed in a miscible displacement at the same mobility ratio. This is true even when oil recovery is expressed on the basis of movable oil. The reason is that oil saturation only slowly achieves its final value behind the water-flood front in accordance with the Buckley-Leverett8 simultaneous flow relations. It is convenient to use miscible displacements for laboratory simulation of waterflood frontal advance since the interfacial tension forces which are negligible in proportion to viscous forces on a reservoir scale9 are thus made nonoperative in the laboratory model. For miscible displacements, the Hele-Shaw type of model adequately represents a porous medium so long as the appropriate scaling rules5-10 are observed in its design and operation. During simulation of waterflood front behavior in the laboratory by using miscible displacements, the behavior of connate water may ordinarily be disregarded since it is usually indistinguishable from flood water in this process. However, when the flood water is deliberately thickened to improve the mobility ratio between water and oil. the effect on the sweep efficiency due to generation of a connate water bank during the process must be considered. In a uniform porous medium, such a bank is generated11 and efficiently displaced by injection of thickened water.l2 The oil originally in-place at the start of the waterflood is then displaced by connate water followed by thickened water. If the flood water must be thickened to obtain a favorable mobility ratio, the mobility of the oil phase is appreciably less than that of the connate water. Hence, the oil phase is inefficiently displaced by the connate water bank, and a considerable proportion of the oil comes in contact with and is displaced by the thickened waterflood front' A recent study 12 of five-spot displacements using miscible Newtonian fluids in which connate water banks of various sizes were simulated shows that, for connate water banks as large as 40 percent of 1 PV, major improvements in sweep efficiency resulted from use of a thickened

Jan 1, 1969

By J. C. Trantham, J. W. Marx

From 1955 to 1958 the Phillips Petroleum Co. conducted a series of small scale counterflow combustion field tests in a tar sand about 60-ft deep and 6 to 12-ft thick near Bellamy, Mo. A total of seven different well patterns, including conventional five-spots and seven-spots plus a 15-well line drive pattern and a 10-well radial pattern, were employed. Electric heaters, gas burners and wellbore fuel packs were tested as counterflow ignition devices. Direct drive ignition followed by reversal to counterflow operation was not possible in this case because of formation plugging by the semi-solid tar during the direct drive phase. Sustained counterflow burning was achieved with wellbore fuel packs under conirolled conditions which permitted close correlation of fire front velocity, air velocity, air-oil ratios and yields for a wide variation in process parameters. The tar in place was about 10' APl with a reservoir viscosity greater than 500,000 cp. The composite oil prodlcct was 26' API with a viscosity in the range of 5 to 15 cp. Vertical sweep efficiency in the line drive pattern was essentially 100 per cent, with no evidence of gravity segregation. Original tar in place ranged between 800 and 1,000 bbllacre-ft, specific rock permeability averaged about 800 md before burning and effective air pemeability with residual oil and water in place (at the time of ignition) averaged 250 md. The volume of oil recovered during stabilized line drive counterflow burning represented about 67 per cent of the volume of tar originally in place in the burned out zones. A minimum was found when air-oil ratios were plotted against formation air velocities, and a limiting air velocity existed below which the counterflow fire front reversed and burned back over its own trajectory, using the residual coke left behind. INTRODUCTION Early in 1955 the Phillips Petroleum Co. began preparations to field test counterflow underground combustion as a method for producing oil from high viscosity hydrocarbon deposits, such as bituminous (tar) sands. At that time all published information on the use of underground burning io increase oil recovery was based on the direct drive process, in which the fire front moved through the formatim in the same direction as the injected air stream. Such direct burning is not applicable to many heavy oil sands or to most tar sands because the heat-thinned native hydrocarbon congeals in the cold rock ahead of the fire front forming a gas permeability block which smothers the combustion. Laboratory experiments first demonstrated that this permeability blocking would not occur if the fire front and the air stream moved in opposite directions. In this counterflow (reverse) burning process the heated heavy oil or tar is driven back through the combustion zone. Here thermal cracking produces relatively light hydrocarbon products which pass through the heated rock behind the fire front (largely in the vapor phase) and are incapable of causing a gas permeability block.' Exploratory wring revealed that consolidated bituminous sandstones with a combination of tar saturation, permeability, thickness and depth favorable for experimental field testing existed at several points in western Missouri. By the fall of 1955, a test site for determining technical feasibility of the counterflow process under field conditions had been selected near Bellamy, Vernon County, Mo. — about fifty miles north of Joplin. The pay zone was too thin and too shallow to constitute a commercial prospect, but these same factors permitted controlled experimental operation with a small compressor installation and with a large number of monitor wells. These were needed for adequate observation of the underground phenomena at reasonable cost. RESERVOIR CHARACTERISTICS The reservoir chosen for the experimental field tests was a 124 thick section of tar sand extending from 49 to 61 ft subsurface, with shale and siltstone laminations sealing the top and bottom. This section was part of a larger 30-ft thick tar sand deposit which extended both above and below the test zone. Core analyses showed that the 12-ft sand was effectively subdivided into two approximately equal sections, with the lower zone being consistently more permeable than the upper zone. The most significant laboratory core data are summarized in Table 1. The line drive test which is the principal subject of this report utilized only the lower 6-ft thick zone, with the upper zone regarded as part of the 55-ft overburden. WELL PATTERNS The over-all field test program involved a total of seven different well patterns, including conventional five-spots and seven-spots plus a 15-well line drive pattern and a 10-well radial pattern. In some of the tests the full 12-ft zone was employed; in other cases only the more perme-

Jan 1, 1967

By J. Lohrenz, C. R. Clark, B. G. Bray

Procedures to calculate the viscosities of in situ reservoir gases and liquids from their composition have been developed and evaluated. Given a composition expressed in methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide together with the molecular weight and specific gravity of the heptanes-plus fraction, the procedures are capable of calculating the viscosity of the gas or liquid at the desired temperature and pressure. The procedure for reservoir liquids was developed using the residual viscosity concept and the theory of corresponding states, and was evaluated by comparing experimental and calculated results for 260 different reservoir oils ranging from black to highly volatile. The average absolute deviation was 16 per cent. This is the first known procedure for calculating the viscosity of reservoir liquids from their compositions as normally available, i.e., including the heptanes-plus fraction. The procedure for reservoir gases uses a sequence of previously published correlations. Evaluation of the procedure was accomplished by comparison of 300 calculated and experimental viscosities for high-pressure gar mixtures in the literature. The average absolute deviation was 4 per cent. The calculations are useful for (I) determining viscosities in compositional material balance computations and (2) predicting the viscosity decrease which occurs when gases, LPG, or carbon dioxide dissolve in reservoir oils. INTRODUCTION Methods to predict viscosities of reservoir fluids from the normally available field-measured variables have been presented. Beal,1 Standing,2 and Chew and Connally3 orrelated oil viscosities with temperature, pressure, oil gravity and gas-oil ratio. Carr, Kobayashi, and Burrows4 and Katz et al.5 have presented correlations for reservoir gas viscosities as a function of temperature, pressure and gas gravity. Lie all intensive physical properties, viscosity is completely described by the following function: Eq. 1 simply states that viscosity is a function of pressure, temperature and composition. These previous correlations1- hay be viewed as modifications of Eq. 1, wherein one assumes more simple functions may be used. The assumptions are practical, because the composition is frequently not known. Further, the assumptions are sufficiently valid so that these correlations are frequently used for reservoir engineering computations. In compositional material balance"' computations, the compositions of the reservoir gases and oils are known. The calculation of the viscosities of these fluids using this composition information is required for a true and complete compositional material balance. For reservoir gases, Carr, Kobayashi and Burrows4 have presented a suitable compositional correlation. For reservoir oils, no correlation is available, and data from reservoir fluid analyses have been used7-9 for compositional material balance calculations.* From a theoretical point of view, this is entirely invalid. The reservoir fluid analysis, whether flash, differential, or other process, does not duplicate the compositions which occur during the actual reservoir depletion process, therefore the viscosities measured during reservoir fluid analysis are not those which occur in the reservoir. From a practical point of view, the "error" of using viscosities from reservoir fluid analysis is of varying and unknown significance. One can say qualitatively that the error is greatest where compositional effects are greatest, i.e., for volatile oil and gas condensate reservoirs and pressure maintenance operations. The first requirement to obtain a quantitative estimate of the significance of the error is to develop a reliable compositional correlation for the viscmities of reservoir oils. No such correlation has been available. Consistent with this requirement, the objective of this study was to develop a procedure to predict the viscosity of reservoir fluids from their compositions. Normally, the compositions of reservoir fluids are available expressed as mole fractions of hydrogen sulfide, nitrogen, carbon dioxide and the hydrocarbons methane through the heptane-plus fraction, with the average molecular weight and specific gravity of the latter. The final correlation was to use the composition in this form. While the more challenging objective of the study was the development of a correlation for the viscosities of reservoir oils, the viscosities of reservoir gases were also studied. The end result of the study was a procedure to calculate the viscosities of reservoir gases and liquids suitable

Jan 1, 1965

By C. S. Land

Relative permeability functions are developed for both two- and three-phase systems with the saturation changes in the imbibition direction. An empirical relation between residual nonwetting-phase saturation after water imbibition and initial nonwetting-phase saturations is found from published data. From this empirical relation, expressions are obtained for trapped and mobile nonwetting-phase saturations which are used in connection with established theory relating relative permeability to pore-size distribution. The resulting equations yield relative permeability as a function of saturation having characteristics believed to be representative of real systems. The relative permeability of water - wet rocks for both two- and three - phase systems, with the saturation change in the imbibition direction, may be obtained by this method after properly selecting two rock properties: the residual nonwetting-phase saturation after the complete imbibition cycle, and the capillary pressure curve. INTRODUCTION Relative permeability is a function of saturation history as well as of saturation. This fact was first pointed out for two-phase flow by Geffen et al. 1 and by Osaba et a1.2 Hysteresis in the relative permeability-saturation relation also has been reported for three-phase systems.3 Since saturations may change simultaneously in two directions in a three-phase system, four possible relationships arise between relative permeability and saturation for a water-wet system. The four saturation histories of this system were given by snel14 as 11, ID, DI and DD. I and D refer to the direction of saturation change (imbibition and drainage), with the first letter of the symbol indicating the direction of change of the water phase. As used in this paper, the second letter of the symbol refers to the direction of saturation change of the gas phase, i.e., D and I indicate an increase and decrease, respectively, in gas saturation. Only a few three-phase relative permeability curves have been published. Leverett and Lewis5 published three-phase curves for unconsolidated sand, and Snell4 reported results of several English authors for both drairrage and imbibition three-phase relative permeability of unconsolidated sands. Three-phase relative permeability curves for a consolidated sand were published by Caudle et a1.3 for increasing water and gas saturations (ID). Corey et a1.6 reported drainage (DD) three-phase relative permeability for consolidated sands. Recently, Donaldson and Dean7 and Sarem8 calculated three-phase relative permeability curves from displacement data on consolidated sands, also for saturation changes in the drainage direction. The only published three-phase relative permeability curves for consolidated sands with saturation changes in the imbibition direction (11) are those of Naar and Wygal.9 These curves are based on a theoretical study of the model of Wyllie and Gardner 10 as modified by Naar and Henderson.11 Interest in three-phase relative permeability has increased recently due to the introduction of new recovery methods and refinements in calculation procedures brought about by the use of large-scale digital computers. The scarcity of empirical relations for three-phase flow, and the experimental difficulty encountered in obtaining such data, have made the theoretical approach to this problem attractive. RELATIVE PERMEABILITY AS A FUNCTION OF PORE-SIZE DISTRIBUTION Purcell used pore sizes obtained from mercury-injection capillary pressure data to calculate the permeability of porous solids.12 Burdine extended the theory by developing a relative permeability-pore size distribution relation containing the correct tortuosity term.13 The works of Purcell and Burdine were combined by Corey14 into a form that has

Jan 1, 1969

By R. J. Wagner, Jim Douglas Jr., P. M. Blair

The calculation of the behavior of an oil reservoir during a water flood has long been an important problem to reservoir engineers. Buckley and Leverett derived the differential equation which describes the displacement of oil from a linear porous medium by an immiscible fluid, but this equation could not be solved by the methods of classic mathematics. Consequently, in order to integrate the equation over the length of the reservoir, they neglected the effects of capillary pressure. In the present paper, a numerical method has been developed for determining the behavior of a linear wafer flood with the inclusion of capillary pressure. The differential equation which was derived for the case of incompressible fluids is second order and non-linear. This differential equation was approximated by an implicit form of difference equation which is second order correct in both time and distance. An electronic digital computer was used to perform the numerical solution of the difference equation. INTRODUCTION The problem of calculating the flow and distribution of fluids in an oil reservoir subjected to a water flood has long challenged the reservoir engineer. The ability to solve this problem would provide a valuable tool for the design and study of field waterflooding programs. One of the first contributions in this field was made by Buckley and Leverett,'.' who developed a method of calculating waterflood performance in a linear reservoir. Their technique was limited by the practical necessity of excluding quantitative consideration of capillary pressure. It is the purpose of this paper to describe a method for calculating the behavior of a linear water flood with capillary pressure considered. This method, although limited to the linear case, should serve as a step toward the solution of the two- or three-dimensional waterflooding problem which would better describe actual reservoirs. The treatment of the problem has assumed that the reservoir is linear and homogeneous, that both oil and water -are incompressible, and that gravitational forces may be neglected, and that water is injected into one end of the reservoir and oil and water are produced from the other. Though these assumptions may not be directly applied to natural reservoirs, they are closely approached in many laboratory core tests. The physical problem has been represented mathematically by means of a differential equation and suitable boundary conditions. The inclusion of capillary pressure effects causes the differential equation to be second order and non-linear, and it is not amenable to solution in terms of known functions. Nevertheless, the problem may be solved by numerical methods. The authors have successfully calculated several cases from the start of water injection through to depletion of the hypothetical reservoir using realistic capillary pressure and relative permeability characteristics that were chosen so as to put the method to a test. DIFFERENTIAL EQUATION In two- or three-dimensional problems involving the two-phase flow of incompressible fluids through a porous medium, the equations describing the flow can be combined until only two dependent variables remain, the pressure in one phase, and the saturation of one phase. These variables appear in two simultaneous partial differential equations which are derived on the basis of the conservation of mass and Darcy's law as applied to two-phase flow. For a linear reservoir, this system of differential equations can be reduced to a single, nonlinear, second order, parabolic differential equation for one dependent variable, the saturation of one phase. For incompressible water and oil, it can be shown that the conservation of mass for the two phases requires that Where q, and q, are the water end oil flow rates per unit cross-section, s, and so their saturations, +, the porosity, and [ and T, dimensional space and time variables. Darcy's law, as extended to two-phase flow including the effect of capillary pressure,' is as follows:

By P. M. Blair

This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. ln the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. INTRODUCTION The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradients, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942.1 Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers.2"4 The influence of capillary forces in laboratory water floods has also been discussed by several authors.6 While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes7 pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured-matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been 'published on the experimental techniques involved?? Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. However, recent advances in the applica-

Jan 1, 1965

By N. W. Ratliff, P. J. Closmann

Production of oil by expansion from a cylindrical reservoir composed of two concentric regions of different properties has been determined as a function of time for a reservoir producing at constant terminal pressure. The parameters involved are permeability, porosity, compressibility and oil viscosity. Results agree with those of a generally accepted method for uniform reservoirs when all parameters are taken as uniform. Cumulative production at any given time is reduced below that of the uniform reservoir when a region of considerably lower permeability adjoins the well; cumulative production increases when a high-permeability region adjoins the well. Curves are presented illustrating quantitative effects of these variations. INTRODUCTION Most petroleum reservoirs can be exploited by release of pressure and consequent expansion of underground fluid. During part of the production history of certain reservoirs by this mechanism, fluid compressibility can be considered small and constant. Cumulative fluid produced by this means as a function of time was calculated for a uniform reservoir by van Everdingen and Hurst.l For certain cases where a production well has been damaged or where an acid treatment has been used, a zone of properties different from those of the reservoir may be created around the production well. In such a case, expansion of the reservoir fluid takes place in a system composed of two concentric cylindrica1 regions of different reservoir properties. Such composite systems were studied in connection with heat flow by Jaeger2 who presented a solution for temperature distribution in a radial system with an infinitely large outer radius. Similar heat flow problems have been studied by other authors Oil production from composite reservoirs was studied for constant production rate by Hurst,4 Loucks and Guerrero,5 and Carter.6 THEORY This work considers constant terminal pressure. Cumulative flow is calculated for a composite bounded system in which no flow takes place across the outer boundary (Fig. 1). Since the flow is considered to be purely radial, the problem involves only one space dimension and time. It is convenient to define the following dimen-sionless pressure drops (Fig. 1).

By Jim Douglas, H. H. Rachford, D. W. Peaceman

The problem of unsteady-state gas flow through porous media has been solved numerically only for the case of linear or radially symmetric reservoirs. A recently introduced numerical method for solving the unsteady-state heat flow equation in two dimensions is applied to the calculation of the depletion of a square region containing a perfect gas. Solutions are presented in graphical form for various values of dimensionless parameters. The solutions are compared with published .solutions for radial reservoirs. I NTR ODUCTION The problem of unsteady-state flow of gas through porous media gives rise to a second-order non-linear partial differential equation for which no analytical solution has been found. Numerical approximations to solutions of the gas flow problem have been obtained by the stepwise solution of an associated difference equation1,2. However, the methods so far developed have required that the reservoir be either linear or radially symmetric. This restriction in shape has been necessary so that only two independent variables be considered, namely, one distance variable and time. In order to deal with reservoirs having more realistic shapes, it is necessary to develop numerical procedures for the solution of the gas flow problem involving two distance variables. A numerical procedure, denoted as the alternating-direction implicit method, for the solution of the heat-flow problem in two dimensions has recently been introduced" By the use of this procedure, approximate solutions have been obtained for heat-flow problems in a square8, and in regions having various non-rectangular boundariesa. Because of the similarity between the equation for heat conduction in solids and the equation for gas flow in porous media, it is reasonable to expect that the alternating-direction implicit method should also be useful for solving the gas flow problem in various two-dimensional regions. Solutions have been calculated for the simplest two-dimensional region, a square reservoir with a single well in the center. While reservoirs of such simple geometry seldom, if ever, exist, the solution of this problem is of some practical importance because it is also the solution of another problem, that of determining the depletion history about an individual well in an infinite uniform reservoir containing wells spaced throughout on a square lattice and producing equally from each well. Finally, it is of interest to compare solutions for a square reservoir with those for a circular reservoir and to determine the effect of the shape of the boundary for that particular case. METHOD OF CALCULATION Basic Differential Equation By combining the equation of continuity, the perfect gas law,

Jan 1, 1956

By L. A. Schrider, J. A. Wasson

A method of predicting waterflood performance has been developed that combines certain facets of several previously published prediction techniques. The manner of combination has required the development and use of some new and some little-known relationships and has eliminated several of the weaker assumptions inherent in the original individual prediction methods. This approach, which was designed for computer solution, has removed the necessity of referring to plotted curves and permits the analytical prediction of water flood performance for a five-spot well pattern in either homogeneous or stratified reservoirs. The Predicted values are expressed in common oilfield units rather than in abstract or dimensionless terms. The calculation procedure has been programmed in FORTRAN IV and the entire program listing is available to potential users.* Introduction Prior to initiation of secondary recovery of oil by waterflooding, the possibilities of economic success or failure must be evaluated with the greatest practical accuracy. Although many prediction techniques are available, most, if not all, have serious limitations. In recent years, the Bureau of Mines has attempted to stimulate interest in the secondary recovery of oil from Appalachian pressure-depleted reservoirs. Waterflooding, because of its widespread acceptance and success, is the principal secondary-recovery method considered. During the life of the project, several performance prediction methods have been studied and applied. Because these methods, individually, have inherent assumptions and limi-tations, portions of some of them have been selected and combined in an effort to eliminate such shortcomings wherever possible. We believe that the resulting technique is as accurate as the necessary available data will permit and that it has several distinct advantages over previously published methods. Development of the Method The fractional flow equation developed by Buckley and Leverett' from Darcy's law is the fundamental mathematical statement describing immiscible fluid-fluid displacement in porous media. A generalized form of this relationship, whose derivation is readily available in the literature, is The formation dip, a, is usually very small in those petroleum reservoirs being considered for pattern waterflooding in the Appalachian area. Therefore. in many cases. a has an almist negligible value. In addition, Ap, the difference between the density of water and the density of oil, is relatively small. For these reasons the gravitational term may often be neglected. Capillary pressure is a function of saturation, which is in turn a function of distance from the upstream boundary of the system or the injection well. In this study it is assumed that the frontal zone of the waterflooding (or zone of capillarity) has negligible length and that only those saturations in the water-swept zone at breakthrough and later are of interest. As the distance in a waterflood pattern is comparatively large and the change in water saturation with distance is usually small, the change in capillary pressure with distance is negligible. The simplified form of the fractional flow equation is appropriate, therefore, in this application. Three variables must be considered in determining that fraction of the total fluid flow that consists of water: the viscosity ratio, saturation, and relative permeability ratio.

Jan 1, 1969

By N. H. Harrison, J. K. Rodgers, S. Regier

This paper presents comparisons of data obtained from a laboratory reservoir study and from a calculated behavior prediction with the actual production history of a condensate reservoir. A small non-commercial discovery was depleted under closely controlled conditions and the well fluids were sampled at frequent intervals. Data on the reservoir and production variables were accumulated on a fixed schedule. A laboratory reservoir study war made using the initial well fluid samples as charging stock. The production procedures and operating conditions were held constant throughout the study wherever possible and in general paralleled the field work. The well fluid compositions and the cumulative recoveries ar a function of the reservoir pressure were also calculated using conventional flash vaporization procedures and equilibrium constants. Comparisons based on the composition of the well fluid show good agreement, the laboratory study agreeing within experimental accuracy with the field work and the calculated data comparing equally well. The gas-oil ratios are also in good agreement, but with somewhat greater deviations at the higher pressures. In the overall picture, it is believed that a model st,~tiy can predict within experimental accuracy the production history of a condensate reservoir. Better equilibrium constants for the heavier hydrocarbons are needed in order to attain improved composition accuracy by calculation. INTRODUCTION In Aug., 1955, a gas condensate well was completed in San Juan County, Utah, that was initially thought capable of good commercial production. These conclusions were derived principally from core data and electric logs, which indicated good permeability, porosity and gas content. However, after the usual series of potential tests it was found that the reservoir pressure had declined some 22 per cent, and it was obvious that the zone tapped was but a small pocket or trap. It became apparent that, with a controlled depletion of a small reservoir, a unique opportunity was available to compare laboratory and calculated studies with an actual field depletion and to further the present knowledge of condensate reservoirs. FIELD WORK The Coalbed Canyon Well No. 1 was conventionally completed in the Paradox limestone formation to a total depth of 5,912 ft. The producing zone from 5,762 to 5,806 ft was perforated with four jet shots per ft. The wellhead and field equipment were also conventional, the major items consisting of a two-pass indirect fired line heater, a high- and a low-pressure separator with the necessary controls and accessories, gas meters, back-pressure regulators, flare stacks and condensate stock tanks. The initial testing of this well con-sisted of a series of flow potential and pressure build-up tests during which some 30 MMcf of gas was produced. The reservoir pressure declined from an estimated 2,300 to 1,782 psig during this period, from which it was concluded that the reservoir was very small. In order to approach steady-state conditions in the reservoir and so provide optimum conditions for making comparisons, the field depletion was programmed to approach, if possible, constant production conditions. Bi-hourly readings were taken of the tubing pressure, the pressure and temperature of the separators, oil and gas rates, and other pertinent operating data. The gas rate, as indicated by the orifice meter, was held constant by the adjustment of the choke in the line heater. The temperature of the first stage separator was held constant by adjustment of the line heater jacket temperature. Practical considerations of production made the maintenance oi a constant gas rate impossible. The test started with a gas rate of 4 MMcf/D and a separator pressure of 250 psig. This rate was maintained until the choke was fully opened. The gas rate then declined with the falling tubing pressure and production was continued until the rate was about 2

By J. Jones-Parra, R. S. Reytor

The porosities of fractured limestone reservoirs can be divided into two broad types in accordance with their effects on fluid distribution and fluid flow. In the coarse porosity, gravity segregation lakes place freely and the resistance to fluid flow is very small. In the fine porosity there is no segregation and a high resistance to flow, and it has relative permeability characteristics similar to tight sandstones. By analyzing the affect of the two porosities it is concluded that in some cases to recover the maximum amount of oil it is necessary to remove large quantities of gas from the reservoir by producing at high gas-oil ratios. In this manner the fine porosity is drained of its oil and the gas-oil contact drops slowly permitting higher production rates from the oil leg. Since this conclusion is contrary to widely accepted principles of conservation, a mathematical model (Fig. I) was constructed to duplicate the conditions desired. The behavior of the model indicates that under certain conditions it is possible to recover more oil by producing at high gas-oil ratios than by production at low gas-oil ratios, and that the rate of production is affected more by gas shut-offs than by the decrease in pressure (Figs. 2, 3, 4 and 5). INTRODUCTION It has been noted that certain limestone reservoirs behave in such a way as to indicate that there are two distinct types of porous spaces available for fluid storage and fluid flow. Regardless of the nature of these porous spaces the distinguishing characteristic is that in one type of porous space there is definite gravity segregation between the oil and the gas; while in the other, the gas evolved tends to remain distributed equally throughout the reservoir. For the sake of simplicity the porosity in which gravity segregation takes place at a fast enough rate to affect the behavior will be called "coarse po- rosity", while that in which this is not so will be termed "fine porosity". The behavior of a reservoir of this type would be different than that of a sandstone reservoir, especially if most of the oil is in the fine porosity. It has been observed that fluid flow to the wells takes place mostly through the coarse porosity. The large fractures which

By F. W. Jessen, G. F. Sears

Cavities in massive salt for the purpose of storage of liquid hydrocarbons have assumed a prominent position in recent years. This paper describes a program to facilitate leaching operations for the formation of specifically shaped storage cavities. Various forms and sizes of cavities may be possible through use of the techniques developed. INTRODUCTION The creation of large under ground storage facilities for natural gas and liquified petroleum gases has been practiced for many years. Use of cavities obtained through solution of salt for this purpose is a fairly new and novel approach, but has gained increasing importance due mainly to the economics of this type of mining operation opposed to hard rock mining or surface and pit storage.1 The idea of controlling the shape of any cavity dissolved from massive salt has not been a prime consideration of companies engaged in the formation of storage space. This has been due mainly to an insufficient knowledge of the mechanics of the leaching process and a dearth of published information dealing with both the desirability and ease of control possible for this type of operation. Two distinct advantages are readily ascertainable. From a stability standpoint (i.e., the ability of the completed cavity to withstand stress imposed by overburden pressure and lateral tectonic stresses) a controlled cavity may be generated which will yield the most favorable attitude to these external forces and remain in operable use for longer periods of time.2 A second advantage is that for a given volume, a sphere (which is one controlIed shape possible) represents the minimum surface area exposed. This may become more important in the future when dealing with refrigeration and product losses in underground cavities. The basic solution mining process, without regard for controlling the final shape, is quite simple in that the equipment and materials required to dissolve a cavity in a salt dome or layered salt section are a source of fresh water, a circulating pump, several strings of tubing and a means of disposal of the return brine. To add control measures involves the use of an inert blanket material above the position at which solution proceeds. Initially the annulus of the largest wash pipe is filled with this blanket down to the top of the proposed cavity. The bottom of the proposed cavity is determined by the depth of the original drilled hole or a blanking plug. The blanket material is added incrementally in stages and displaces the water in the enlarging cavity downward. Where the blanket material has displaced the water, no further solution takes place. The rate at which to add this controlling blanket material so as to be able to form specifically shaped cavities of any particular size is of considerable importance. THE PROBLEM It is desirable to have as much information as possible as to the behavior of the mining operation since visual observation obviously is impossible. It would be advantageous, for example, to have a step-by-step program showing how much salt is to be removed, at what rate and the time required. This type of program would facilitate procurement of surface equipment and provide for proper management of manpower requirements, as well as yield a host of intangible benefits. As mentioned earlier, only rule-of-thumb estimates were available until recently and made this type of pre-planned operation haphazardous at best. Recent work3-' has done much to clear up the ambiguities and inaccuracies attendant to solution mining of salt and has helped to place the operation on a more scientific basis. This study attempts to correlate much of the information thus far developed into a complete mining program; to extend the idea of controlled solution mining to include not only spheres, but solid conic sections of the ellipsoidal variety; and to refine the computation of the rate of removal of salt during the various stages of mining. Dommers3 and IIemson4 made extensive laboratory

Jan 1, 1967

By J. G. Giddings, R. Kobayashi

Residual viscosity, the ViSCOSitY at a given pressure and temperature minus the dilute gas viscosity at the same temperature, has been found to be independent of temperature for pure components and mixtures for simple molecules and light hydrocarbons. The residual viscosities for methane, ethane, propane and four mixtures were correlated as a function of reduced density and molecular weight. INTRODUCTION Earlier efforts in correlating the viscosity of light hydrocarbons and associated gases under pressure, primarily by some form of the corresponding states principle, have been reported elsewhere.7 The formal theory of corresponding states principle for viscosity from non-equilibrium statistical mechanical arguments for simple molecules has been developed by Helfand and Rice.13 The most serious obstacle to the development of a generalized corresponding states representation of viscosity data, even for simple molecules, has been the limited quantity of precise viscosity data of pure components and mixtures over a wide range of pressure and temperature conditions. Viscosity data on light hydrocarbons have recently been contributed by several laboratories.2,6,8,23,25,26 In the study of the viscosity of nitrogen, Michels and Gibson" observed that the viscosities of nitrogen as a function of density gave parallel isotherms. The same behavior has since been reported for other gases.".= The observed behavior indicates that the residual viscosity or the viscosity at a given pressure and temperature minus the dilute gas viscosity at the same temperature is a unique function of the density. Abas-Zade,1 Michels and Botzen,18 and Michels, Bot-Zen, Friedman and Sengers,19 observed similar behavior for thermal conductivity, the residual thermal conductivity being defined in a manner parallel to that of residual viscosity. Thodos and co-workers 4,15,24 applied the residual viscosity concept to the viscosity of pure monatomic and diatomic gases. Ellington and co-workers10,25 made the residual viscosity vs density plots for ethane and propane, but omitted the region below a residual viscosity of 15 micropoise and showed a systematic separation of the ethane isotherms in the low density region. Eakin correlated the residual viscosity of methane, ethane, propane and n-butane empirically as a function of density and molecular weight: Carmichael and Sage6 recently presented ethane residual viscosity vs density values in the gaseous liquid, and near critical regions as a single univariant function of density. Flynn, Hanks, Lemaire and ROSS" have contributed extremely accurate viscosity data to confirm that the residual viscosity concept is essentially correct, but may not be accurate within 0.1 per cent of the viscosity value for the gaseous state- The recent studies by Giddings12 indicate that the residual viscosity of methane, propane and each of four methane-propane mixtures are unique functions of density Over the temperature and the density ranges studied.

Jan 1, 1965

By W. D. Von Gonten, R. L. Whiting

Regression analysis was used to correlate the physical properties of 478 sandstone and 90 carbonate core samples. Porosity, permeability, electrical formation resistivity factor, capillary pressure and the sonic velocity of the shear and compressional waves were measured. Prediction equations for porosity, permeability and electrical formation resistivity factor were found which should be useful in understanding the relationships between the physical properties of porous media in formation evaluation. INTRODUCTION The physical properties of porous media are important to the petroleum engineer and geologist. To evaluate fully the potential and behavior of a subsurface formation such as a petroleum reservoir, certain physical properties of the porous medium must be known. The problem of accurately determining physical properties of subsurface formations has not been solved because many determinations must be made by indirect measurements and because of the difficulties caused by complex pore structure and presence of clays in most naturally occurring porous media. Due to the difficulty in measuring some of the physical properties of porous media, it would be advantageous to be able to predict a certain physical property of a rock from other physical properties of the rock which could be measured more easily and more accurately. Since most potential reservoir rocks are heterogeneous, relationships between the physical properties are very complex and thus far no satisfactory correlations based on theory or laboratory models have been developed. It appears that empirical relationships obtained by measuring the physical properties of a large number of samples of naturally occurring porous media and applying regression analysis to develop equations for one physical property in terms of other rock properties is the best approach. REVIEW LITERATURE The three physical properties used as dependent variables for correlating purposes were porosity, permeability and formation factor. Since formation factor is more difficult to determine, a brief review of the literature is provided on this property. The first work on determining formation factors was published by Archie1 in 1942. He defined this property of a porous medium as F - Ro/Rw......(1) where Ro is the resistivity of the porous medium when completely saturated with a brine of resistivity Rw. Archie found the best correlation between formation factor and porosity was the following equation, where is the porosity fraction and m is the constant characteristic of the rock. The value of m was 1.3 for unconsolidated sand packs, and ranged from 1.8 to 2.0 for consolidated sandstones. In 1950 Patnode and Wyllie2 and De Witte3 observed that the formation factor as determined by Eq. 1 was valid only when the porous medium contained no conductive solids such as clay or shale. When conductive solids are present the formation factor is also dependent on the resistivity of the saturating fluid. Therefore, in samples containing conductive solids the formation factor decreased as the resistivity of the saturating fluid increased. Because of this. the measured formation factor was called the apparent formation factor and was designated Fa. Patnode and Wyllie proposed the following equation,

By L. A. Wilson, P. J. Root

The relative costs of heating a reservoir by steam injection and by combustion have been examined. The comparison was based on a model similar to that proposed by Chu.' The cost of boiler feed water, the price of fuel, pressure and plant capacity were parameters in determin-ing the costs of air compression and steam generation. The analyses compare the cost of heating to the same radius by the two methods. Results suggest that the two primary factors for comparison are the price of fuel and the amount of crude burned during underground combustion. The cost of fuel has a greater effect on the cost of heat from steam than it does on its cost by combustion. As a result, analyses indicate that when the price of fuel is low, steam may be unequivocally cheaper than air. The influence of heat loss is such, however, that as the heated radius increases combustion becomes relatively more competitive depending upon the amount of crude burned. This implies that steam may be cheaper for small stimulation jobs (huff and puff) but combustion may be more economically attractive for heating large areas (flooding). INTRODUCTION Use of thermal methods of recovery is an accepted fact today. After an induction period of several years, processes are being widely used that involve reservoir heating to augment recovery. Of the several techniques, steam injection and forward combustion appear to be destined to dominate the field. Although the objectives of both are the same, the basic differences between generating heat in situ and injecting heat after surface generation influence the cost in different ways. This study compares the cost of heating a reservoir either by steam injection or by forward combustion. There has been no consideration of recovery. Presumably, recovery from the swept region would be high in either case. The sole consideration was the cost of heating to the same radial distance by either process. PROCEDURE THE MODEL The basis for comparison was a mathematical model similar to that used by Chu' for combustion. The model simulates a radial heat wave in two-dimensional cylindricaI coordinates. It includes heat generation, conduction and convection within the reservoir and conduction in the bounding formations. Thus, heat losses from the formation are considered. Three significant modifications were made. 1. Equal logarithmic increments rather than equal increments were used for the mesh spacing in the r direction. By this technique large distances were simulated with relatively few mesh spaces. 2. A backward difference approximation to the convection term was used to avoid troublesome oscillations which result from a central difference approximation when the convection term is large. 3. The radial increments of the combustion zone motion were not necessarily uniquely related to the mesh configuration. The cumbersome step function introduced by the heat of vaporization of steam was circumvented by assuming the enthalpy of the steam to be a linear function of temperature between reservoir temperature and steam temperature. This is equivalent to assuming an average heat capacity numerically equal to the difference between the enthalpy of saturated steam and the enthalpy of water at reservoir temperature divided by the difference between the two temperatures. Heat losses obtained by this model are in essential agreement with those obtained by the analytical solution of Rubenshtein.' A detailed description of the model is presented in the Appendix. Using the model, the times required to heat to particular radial distances were obtained as a function of injection rate and other physical parameters. For the steam case, injected fluid was assumed to be saturated steam at pressures of either 500, 1,000 or 1,500 psia. The corresponding temperatures are 467, 544 and 596F, respectively. Thickness ranged from 10 to 50 ft and injection rate ranged from 100,000 to 1 million Ib/D. Reservoir and overburden temperatures at the injection well were assumed to be that of saturated steam at the injection pressure. The effect of maintaining the overburden temperature at the well at a different temperature (initial reservoir temperature) was examined with no significant change in behavior. The influence of wellbore heat losses for the steam case was determined in the following manner. The rates of heat loss as a function of time were estimated using an approach similar to that suggested by Ramey." he data were based on injection through 2%-in. tubing in 7-in. casing. Integration of these data over the entire iniection period yielded the total heat loss. Total heat losses were then corrected to their equivalents in steam (this number resulted from dividing the total heat loss by the latent heat). This was considered additional steam required to accomplish the reservoir heating and the total cost was increased accordingly.

Jan 1, 1967