Search Documents

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By J. C. Griffiths, D. W. Bennion

A mathematical model, which does not assume a priori that stratification exists, but was designed to test for the stratification was developed. The model segmented the reservoir horizontally into areas of common variance, then divided it vertically into strata (if strata were present). Next, trend surface techniques were used to determine the lateral extent and variation of each stratum. The model was tested on two reservoirs (sandstone and limestone) using porosity as the reservoir rock property. The sandstone reservoir contained approximately 60,000 samples from 2,000 cored wells, while the limestone reservoir had approximately 24,000 samples from 430 cored wells. Within the areas of common variance tested, the model was able to distinguish four separate lithological units or zones in the sandstone reservoir, four in the Marly section and seven in the Vuggy section of the limestone reservoir. INTRODUCTION To accurately predict the movement and production of fluids from a reservoir rock, it is necessary that the reservoir engineer have a knowledge of the rock properties and their distribution throughout the reservoir. For this study, a reservoir rock is defined as a solid containing interconnecting holes or voids which occur relatively frequently and are dispersed within the rock in a regular or random manner. The fraction of bulk volume that these interconnecting voids occupy is called porosity. The rock must also possess sufficient permeability to allow fluids to move through it and it may be composed of one or more strata or lithological units. A stratum or lithological unit as defined in this study is a body or volume of reservoir rock whose properties are so distributed that its lateral and vertical extent can be traced throughout or through a portion of the reservoir. Operationally, the existence of strata are question- able unless the variance within strata is less than among strata. Within the past few years computers have been widely used as a tool for calculating movement and production of fluids from reservoirs. Most prediction models require some type of a mathematical description of certain reservoir rock properties. In the past, most of these models have been relatively simple. Some have assumed vertical variation but no lateral variation, while others have assumed lateral variation but no vertical variation. In general, most of these models have not used functional relationships to predict reservoir rock properties as functions of position. The object of this study was to develop a model which would predict lateral and vertical variations of reservoir rock properties. Since samples are usually available from only a small portion of the total reservoir rock, it seemed logical that if measurements from these samples were to be used to infer the properties of the actual reservoir, the data should be treated statistically as a sample from the total population (reservoir). Therefore, the model developed in this study to predict reservoir rock properties was a stochastic model. The model was tested using porosity (a macroscopic reservoir rock property) measurements from two reservoirs — the first, a sandstone reservoir from which 60,000 samples were obtained; and the second, a limestone reservoir with 24,000 samples. Several investigators have proposed methods to determine the vertical and areal variation of a reservoir rock property. Stiles: Dykstra and Parsons 2 and Suder and Calhoun3 have all developed waterflooding prediction techniques which assume vertical variation but no lateral variation. Stiles developed a method to segment a reservoir arbitrarily using frequency distribution. Law4 suggested that porosity has a normal frequency distribution and that permeability has a log-normal frequency distribution. Elkins and Skov 5 proposed that when a reservoir is segmented vertically, the

Jan 1, 1967

By C. S. Matthews, G. L. Stegemeier

In one field in South Texas, approximately 72 per cent of the pressure build-up results show a characteris-i.rtic "hump" (i.e., the pressure builds up and then falls off) which makes interpretation by. standard methods impossible Correlation of size and time of hump with formation permeability. well productivitv index, and method of completion led to the tentativc, conclusion that the humps were caused by segregation of gas and oil in the wellbore after closing-in"'. This conclusion was confirmed by performance of simple laboratory bubble-rise experiments, by theoretical bubble-rise time calculations, and by a detailed calculation of PVT behavior in the well of a particular well Oil which accurate Surface and bottom-hole pressure measure-ments were untile. The hump behavior has since been found to occur in many other fields. The cause, however, is not the mine in all cases. In some of these the hump is trace-ahie to leaks in the tnbing which allow influx from the annulus after closing-in. In other cases the hump is traceable to leaks in the device separating pay horizons in dually completed wells. It is concluded that the recording of both surface and bottom-hole pressures is desirable in wells which show an anomalous build-up behavior. A number of field examples is discussed where use of both sets of measurements enables the cause for anomalous behavior to he found, (and a reasonable interpreation of bottom-hole pressure to be made. INTRODUCTION Theoretically the pressure build-up in an infinite reser-voir should be a linear function of In [(t + at)/nere t is the production time and At a the closed-in time. Some of the variations from this behavior arc well known, such as the curved portion immediately after shut-in which results from after-production and skin effect, and the flattened end portion which results from boundary effects in a limited reservoir. The effect of stratified producing zones and irregular geometrical drainage patterns may also contribute unusual characteristics to build-ups. However, the effect of still another phenomenon—that of movement of fluids within the wellbore—has been neglected in most build-up studies to date. Fluid movements within the wellbore after shut-in occur as a result of after-production, packer failures, leaks in the casing or tubing, or buoyancy of the gas phase when both gas and liquid are present in the well. Each of these movements can influence the pressure build-up, sometimes sufficiently to negate use of the data for computing permeability or stabilized build-up pressure. Examples of each of these phenomena are discussed below. FIELD OBSERVATIONS ON HUMP BUILD-UPS In approximately 75 per cent of the wells in a medium-sized field in South Texas the pressure build-up curves rise to 3 maximum and then decline to what appears to be a stabilized reservoir pressure as shown in Fig. I. Similar behavior has also been observed in wells in other fields in Texas and Louisian:,. Because of the possibility of such behavior being a mechanical failure of the pressure bomb. the bombs were run in tanden into one well which had been known to exhibit this behavior. When both bombs recorded identical curves with the characteristic hump, it was inferred that the unusual build-up is a result of well or forma-

By M. R. Tek, K. H. Coats, D. L. Katz

In natural gas storage operations, seasonal pressure fluctuations in the gas reservoir cause the water from the surrounding aquifer to flow into and out of the gas sand. The theory of unsteady-state liquid flow through porous media developed by Van Everdingen and Hurst' has been applied to predict the water movement into and out of the gar bubble for .several .postulated pressure cycles. Those cycles with as many pound-days above the original aquifer pressure as pound-days below, may cause the gas reservoir to slowly grow or shrink rather than hold to a constant volume. Applications to an actual field case study give the predicted gas reservoir monthly pressures and volumes as compared with the observed monthly pressures and volumes. INTRODUCTION A large number of natural gas storage reservoirs are bounded by or adjacent to large water-saturated formations called aquifers. The presence of these aquifers is usually evidenced by the production of water from wells delineating the gas-hearing sands. Some gas storage reservoirs in use today have been purposely developed on aquifers. Cyclic pressure variations in a gas storage reservoir cause water influx and efflux from the surrounding aquifer. This, in turn, results in a varying gas reservoir volume. The prediction of the effect of this aquifer fluid movement on the size and size variation of the gas reservoir can provide information valuable in the study of several reservoir engineering problems. The economics of gas storage operations are directly influenced by the influx of aquifer fluid into the reservoir since a shrinking storage reservoir requires increasing pressures for the storage of the same quantity of gas. Material balance and reserve or recovery calculations also obviously require knowledge of reservoir volume. Various other reservoir engineering calculations such as interpretation of well interference data, evaluation of physical characteristics of porous media, water coning, gas injection and pressure maintenance studies are typical problems where variations of volume due to edge or bottom-water encroachment becomes important. The solution of equations describing influx and efflux of water in the gas sand, application to gas storage reservoir behavior and interpretation of results are the main objectives of this paper. While this work was directed toward gas storage reservoirs, it applies equally well to repressuring of partially depleted oil reservoirs with an active water drive. MATHEMATICAL ANALYSIS As just mentioned, the cyclic pressure variations in an aquifer-surrounded gas storage reservoir cause water movement and varying gas reservoir volume. Prediction of this volume variation from the reservoir pressure cycles or the injection-production schedule requires adoption of a flow model, certain assumptions concerning reservoir and fluid properties, and formulation and solution of the governing differential equation. The flow model used in this study is a gas bubble located on a very large aquifer (Fig. la). Gas bubble radius and aquifer sand thickness are denoted by rr and h, respectively. Fig. lb shows the horizontal, radial two-dimensional porous reservoir attended by the unsteady, compressible motion of the aquifer liquid. Volume variation of the gas bubble is determined from aquifer water flow across the cylindrical boundary, r = r0

By R. V. Higgins

This paper presents a computer method to obtain the shape factors and equal cell volumes of the channels for any well spacing pattern from a potentiometric model. By using this program the authors have processed the data for the seven-spot, direct line-drive and the staggered line-drive patterns. The data for the five-spot pattern had been previously processed by a noncomputer method and are included for completeness. The shape factors and volumes for the channels are presented in tables for those who want to use them to process data using their own permeability relationships and viscosities of their reservoir oils. The authors have used the data and sets of representative permeability curves to process sample calculationr of waterflood performances. The comparison of the calculated results shows that the influence of well spacing is small. The permeabilities of the reservoir rock to oil and water had a greater influence on oil recovery for a given pore-volume throughput of water than the well spacing pattern. The more water-wet the reservoir rock, the better the possibility of permeabilities which are conducive to good recovery. The viscosity of the reservoir oil also influences the recovery more than the well spacing pattern. The reduction in the percentage recovery of oil with increase in viscosity of the reservoir oils is small when oil viscosities are in the range of 0.1 to 3. Above this range the reductions in recoveries are extensive. Sample comparisons of the time required for different patterns to recover the oil are presented. Results of an example calculation are given to show the effect of the permeability profile on recovery. INTRODUCTION The effect of well spacing pattern on the recovery of oil when flooding with either gas or water has been studied by many investigators. Muskat et al.1 presented an analysis using conductivity, sweep efficiency and unit mobility to the time of breakthrough. Dyes et al.2 used experimental techniques (X-ray shadowgraphs) and different mobility ratios. They presented quantitatively the relationship between mobility and sweep efficiency at and after breakthrough. Hauber3 presented a method to predict waterflood performance for arbitrary well spacing patterns and mobility ratios. Craig et al.,4 using techniques similar to Dyes et al. to determine sweep efficiency for a five-spot pattern, added the use of relative permeability curves at breakthrough and thereafter. Douglas et al.5 sed rela- tive permeabilities and continuously changing saturations throughout the entire five-spot flood pattern. In obtaining their solutions they used finite-difierence equations. Hig-gins and Leighton6,7 also used relative permeabilities and continuously changing saturations throughout the pattern before and after breakthrough. They employed techniques that process a flood-pattern calculation on the computer in about one minute. The methods of Douglas et al. and Higgins and Leighton both checked closely the laboratory results for a wide range of mobility ratios. This paper presents some sample performances calculated by the Higgins and Leighton method that show the effect on recovery of different permeabilities and viscosities using the seven-spot, the line-drive and the staggered line-drive, as well as the fivespot flood pattern. No previous paper has presented these data using different permeability curves and continuously changing saturations throughout the flood patterns. The paper also presents (1) the results and analyses of the flood-pattern prediction, (2) the computer techniques for determining the shape factors and volumes from the potentiometric models for the foregoing flood patterns, and (3) the shape factors and volumes of the channels of the flood pattern in the event reservoir engineers may like to process waterflood calculations using their own permeability curves and reservoir oils. DESCRIPTION OF METHOD VOLUMES AND SHAPE FACTORS The use of channels taken from a potentiometric model (see Fig. 1) to aid in calculating the performances of water floods of nonlinear patterns has been thoroughly explained in the literature.0,7 Therefore, very little theory, discussion, or proof regarding this phase will be repeated in this paper. The computer method presented in this paper to calculate the volumes and the shape factors of the channels of potentiometric models employs the trapezoidal rule for the volumes and the Pythagorean theorem (the hypotenuse equals the square root of the sum of the squares of the two sides of a right triangle) for the shape factors. In calculating the volume of a channel, the area of each cell in the channel is determined and then multiplied by the thickness to obtain the volume. In determining the areas of the cells, trapezoids are constructed whose vertical sides are spaced Ax apart, as shown in Fig. 2. The length of the sides is the difference between an ordinate cut off by the top and bottom of the cell — usually equipo-tentials. The coordinates of the points along an equipo-tential or streamline are obtained by Lagrange's8 equation of interpolation for which the constants are coordinate points at the intersections. Three intersections for con-

Jan 1, 1965

By R. A. Greenkorn, R. C. Smith

Hele-Shaw cells are used to model creeping flow through porous media (where Darcy&apos;s law is valid). The effects of inertia on flow about obstructions in a Hele-Shaw cell can be calculated by a perturbation method if one can determine a solution to Laplace&apos;s equation. Results of a computer solution for flow about circular, square and elliptical obstructions are presented. These results show that or a modified Reynolds number of less than I, the inertia terms are small; and for values of less than 3, the average streamline predicts the ideal flow. Therefore, the analogy might be used for studying flow in porous media up to a modified Reynolds number of at least 3. INTRODUCTION The nature of fluid flow in porous media is of interest in the fields of soil mechanics, ground water flow, petroleum production, filtration and flow in packed beds. Because it is very difficult to study the phenomenological behavior of flow in porous media, homologs and analogs are used to study flow characteristics. A Hele-Shaw model, made of two closely spaced plates - usually glass — is often used as an analogy to two-dimensional flow in porous media. Hele-Shaw2 showed experimentally that the streamline configuration for creeping flow around an obstacle located between two closely spaced parallel plates is the same as for two-dimensional, ideal flow about the same obstacle. Stokes7 verified these observations mathematically. The usual equation of motion for flow in porous media is Darcy&apos;s law. The form of the mathematical statement of Darcy&apos;s law is identical, within a multiplicative constant, to the expression for the average velocity over the plate gap in the plane of a Hele-Shaw model. These models may be used to describe flow in both homogeneous and heterogeneous porous media. In the mathematical proof of the Hele-Shaw analogy it is assumed that the convective terms in the Navier-Stokes equations are negligible and that the equations of motion degenerate to Laplace&apos;s equation, with pressure the dependent variable. Whenever a Hele-Shaw model is used as an analogy to flow in porous media, .the validity of this assumption is in question. Riegels5 showed that if convection is not neglected, the velocity distribution around a cylindrical obstruction in the flow field depends on a Reynolds number, the plate spacing, and a dimension characteristic of the obstacle. Riegels&apos; solution, a perturbation solution, uses the boundary condition that the flow rate into the obstacle averaged over the plate gap at any point on the obstacle is zero. The method requires that a solution to Poisson&apos;s equation for the perturbation pressure be found. Riegels evaluated this solution for the case of the cylindrical obstruction. His method may be simplified by eliminating the need for solving Poisson&apos;s equation for the perturbation pressure. Instead, an analytic expression for the perturbation pressure gradient is obtained (valid for arbitrary shapes) and used to eliminate pressure from the equations for the perturbation velocities. The results show, for symmetrical shapes, that if NRe < 1, the convective acceleration terms are small, and that the average velocities represent ideal flow up to at least NRe = 3, where: L is a characteristic dimension of the obstacle perpendicular to flow, b is the plate spacing, u is viscosity, v, is velocity of approach and p is density. THEORY* Incompressible fluid flow through a Hele-Shaw model and a porous medium are analogous because the equations for the average velocity components in two-dimensional flow are identical if:

Jan 1, 1970

By G. B. Thomas

Several methods of analyzing pressure build-up data in wells have been presented by various authors. This paper reviews the theory and method of D. R. Horner and presents example calculations performed on data obtained by testing several different types of wells. These calculations include, (1) graphical estimation of final static pressure, (2) determination of the productive capacity of the pay away from the well bore and. (3) the degree to which the formation adjacent to the well bore has been damaged by completion or other causes. The methods of testing and precautions which should be taken to assure the best data possible are discussed. Limitations and reliability of calculated results are also treated. INTRODUCTION Pressure testing of wells is generally limited to the determination of producing and static mean formation pressures. The so-called "static" pressure determination.. along with PVT, electric log and production data, enable the reservoir engineer to determine, within reasonable limits, the drive mechanism of the reservoir and in some cases. the amount of edge water encroachment. Producing pressure tests enable calculation of productivity indices and allow the engineer to plan the systematic production of a pool for optimum conservation of suhsurface energy. The radial flow formula advanced by Muskat&apos; has been based on the assumption of incompressible radial fluid flow. It has been known that reservoir fluids do not behave in an ideally incompressible manner. For example, incompressible flow theory indicates a simple logarithmic relationship between the difference of the instantaneous and static well pressures when plotted against time. &apos;The latter stages of this type of plot of pressure build-up data generally show a marked devia-, tion from the earlier straight line trend. which deviation may be shown to be due to the compressible flow of fluids toward the well bore. Shut-in times of 24, 48. 72. or at most 96 hours are currently in wide use for determining so-called "static" reservoir pressure. Due to the continuation of compressible flow of fluids into the well bore long after this arbitrarily taken shut-in time, the determination of static pressures has almost invariably resulted in lower than equilibrium values. Materials balance calculations made early in the life of a reservoir often result in a calculated reserve which later observations prove to be too low. Failure to obtain reliable "static" reservoir pressures within the prescribed 24 or 48-hour build-up period has undoubtedly been a major factor in obtaining these low estimates. Comparison of theoretical and actual productivity performance has indicated that formation damage and not being able to attain true static pressures have been partially responsible for the observed discrepancies. Research into the theory of compressible flow behavior has resulted in methods of applying these theories to the testing of wells. Application of the developed theories to pressure build-up performance indicates that in many cases the following can be estimated. 1. The static reservoir pressure. 2. The in-place effective formation permeability away from the well bore. 3. The degree or extent by which the formation has been damaged adjacent to the well bore, either through completion methods or subsequent damage due to fluid entry. REVIEW OF THEORY Horner2 has shown that the pressure build-up within a "point source" well is approximated by the following formula: Equation (1) is the approximate "point source" solution to the radial compressible flow equation advanced by Muskat. The solution assumes the following: 1. A point source well is producing at constant rate from the center of an infinite reservoir with a constant pressure at its external boundary. 2. The fluid flowing is present in one phase only. 3. The compressibility and absolute viscosity of the fluid remain essentially constant over the range of temperature and pressure variation encountered. 4. The well i.; shut-in at the sand face arid there is 110 after production into the well bore. 5. The formation permeability is homogeneous in the direction of flow. From Equation (1) it can be seen that if an ideal well were shut-in while producing from a reservoir under the conditions assumed, the pressure build-up would be a logarithmic func-

Jan 1, 1953

By G. W. Lester, T. J. Nowak

The paper presents a practical method of utilizing pressure fall-off data obtained when a water injection well is shut in for determination of: (I) the static reservoir pressure, (2) the potential water injective capacity (md-ft) of the formation, and (3) the formation damage (skin effect) surrounding the wellbore. In the development of the method, equations are presented which give The theoretical background for predicting the behavior of water injection wells. Certain physical factors which limit the application of the method are discussed. The procedure and the relevant precautions in obtaining the data are discussed. A sample analysis and computation from an actual pressure fall-off curve is illustrated. Fintrlly, general characteristics of pressure fall-off curve! obtained in several reservoirs are presented. INTRODUCTIO N Theoretically, in an injection well the flow of water into the formation will continue as long as the sand-face pressure is greater than the reservoir pressure. However, in practice, engineers have observed that upon stopping injection and opening the wellhead to atmospheric pressure, some wells backflow water for time intervals ranging from several seconds to several hours although the sand-face pressure, as calculated from the hydrostatic head, is considerably greater than the static reservoir pressure. Further, field observations indicated that upon shutting in a water input well, the wellhead pressure does not drop off immediately, but a pressure decline occurs lasting from several seconds to several days. These observations make it evident that when an input well is shut in, a residual pressure lingers in the vicinity of the wellbore which is greater than the average reservoir pressure. This residual pressure declines first rapidly, then more slowly, to the static reservoir pressure: therefore, it is called transient back-pressure or pressure fall-off. The first step in investigating the mechanism of pressure fall-off is to establish what general conditions are necessary for its occurrence, and to this end, Dickey and Andresen1 have presented a qualitative explanation describing the mechanism. They point out that the formation pressure gradient in the vicinity of the well does not disappear immediately when the input well is shut in, as if the system were filled with a totally incompressible fluid, but it is maintained by the expansion of the compressible fluids and by the elasticity of the formation. In such a compressible system the decrease in pressure is possible only to the extent that the water front continues to advance until the pressure is equalized. The principal aim of this paper is three-fold; namely. (1) to present a theoretical quantitative analysis of pressure fall-off data, (2) to present an application of the developed theory to field data, and (3) to discuss the general characteristics of actual pressure fall-off curves. Under (I), theoretical equations dealing with such factors as compressibility of the fluid-reservoir system, mobility of the fluids and certain physical restrictions that

Jan 1, 1956

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&apos;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:.&apos; 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&apos;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,&apos; 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&apos; with the water influx equation&apos; 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.&apos; 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&apos; 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&apos; APl with a reservoir viscosity greater than 500,000 cp. The composite oil prodlcct was 26&apos; 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.&apos; 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"&apos; 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,&apos;.&apos; 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&apos;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&apos;s law, as extended to two-phase flow including the effect of capillary pressure,&apos; 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 &apos;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).