Search Documents

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

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

Jan 1, 1970

By W. A. Klikoff, I. Fatt

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

By H. J. Ramey

Conventional calculation of total system isothermal compressibility for a system containing a free gas phase involves, among other things, evaluation of the change of oil and gas formation volume factors and the gas in solution with pressure. Preferably, this information should be obtained from laboratory measurements made with particular oils and gases. Often, experimental measurements are not available. In this case, it is necessary to obtain pressure-volume-temperature relationships from general correlations such as those of Standing' for California oils. In order to speed estimates of compressibility, generalized plots have been prepared of the change of both oil formation volume factors and gas in solution, with pressure from Standing's correlations. A generalized plot for estimating the change in the two-phase (oil and gas) formation volume factor with pressure is also presented. Usually, the effect of gas dissolved in reservoir water upon the total system compressibility is neglected for gas saturated systems, due to the low solubility of gas in water. Results of this study indicate that the increase in total system compressibility caused by solution of gas in water is often as large as the compressibility of wafer, and can be magnitudes larger for low pressure systems. Generalized results for estimating the change of gas in solution in water with pressure are presented in tabular and graphical form. INTRODUCTION During the past decade, pressure build-up and drawdown techniques have gained an important place in reservoir engineering. Build-up and drawdown analyses are only two special applications of the broad field of transient fluid-flow theory. All solutions of transient fluid-flow problems contain a parameter called the total system isothermal compressibility. This property of fluids and porous rock is a measure of the change in volume of the fluid content of porous rock with a change in pressure, and it may vary considerably with pressure. Evaluation of total system isothermal compressibility is not difficult, but it is tedious and time-consuming. Often compressibilities are estimated roughly, or transient flow methods are neglected completely. The benefits of using accurate system compressibility in properly-executed build-up or drawdown analyses are: 1. Better planning of pressure build-ups may be achieved to avoid unnecessary loss of revenue due to excessively long shut-in periods, or to shut-in periods too short to yield useable data. 2. Better and more reliable estimates of static formation pressures for reserves estimates and rate performance estimates. 3. Reliable information for evaluation of well completion effectiveness, and planning and interpretation of well stimulation efforts. The purpose of this paper is to clarify the nature of the total system isothermal compressibility, and to present useful methods for estimation of compressibility, particularly for systems containing a gas phase. DEVELOPMENT Numerous publications have presented solutions to transient single-phase flow of slightly compressible fluids, stressing pressure build-up applications. In transient flow, a compressibility* term arises to permit volume content of fluids in porous rock to change as pressure changes. The basic nature of the compressibility term is usually taken for granted. Problems arise in practical applications of transient fluid theory because most published works consider only one flowing fluid—in an ideal porous system containing only one fluid. In 1956, Perrine2 presented an intuitive extension of single-phase flow pressure build-up methods to multiphase flow conditions. Later, Martin- established conditions under which Perrine's multiphase build-up method had a theoretical foundation. Perrine has shown that improper use of single-phase build-up analysis in certain multiphase flow situations can lead to gross errors in estimated static formation pressure, permeability and well condition. It is likely that, much pressure build-up data for oil wells should be analyzed on the basis of multiphase flow. For both single-phase and multiphase build-up analysis, the isothermal compressibility term in dimensionless time groups often should be interpreted as the total system compressibility. All real reservoirs contain one or more compressible fluid phases. In addition, rock compressibility can contribute in an important way to the total system compressibility. The proper total system compressibility expression may contain terms for compressibility of oil, gas, water, reservoir rock and terms for the change of solubility of gas in liquid phases.

Jan 1, 1965

By J. D. Rice, S. C. Pitzer, C. E. Thomas

Interpretation of pressure build-up data obtained in the conventional manner has often been difficult because of the deviation from theoretical behavior. Major causes of this deviation have been attributed to damage and afterflow, and to fluid redistribution in the wellbore which, in extreme cases, can result in a pressure hump in the early portion of the build-up curve. Theoretical investigations show that bottom-hole pressure is definitely influenced by phase redistribution in the tubing column during surface shut-in tests, and that the magnitude of the effect is sensitive to the producing gas-oil ratio and stabilized rate of flow in the well. For field experiments a wire-line tubing packer, which can be run in the tubing against a stabilized flow rate, was developed for bottom-hole shut-in tests. By use of this bottom-hole .shut-in method, pressure humps previously observed in surface shut-in tests were completely eliminated and the effects of afterflow minimized. Analog and digital computer studies have been made to obtain theoretical curves for comparison with field curves, and remarkable agreement between the results of bottom-hole shut-in tests and theoretical curves has been obtained. INTRODUCTION Pressure build-up data from shut-in wells have been used by the petroleum industry to determine the permeability of the formation, to estimate wellbore damage, to evaluate the static reservoir pressure and to speculate reservoir volume. Calculation of these factors is based on methods of analysis developed for theoretical systems whose build-up curves have a charac- teristic shape when the wells are shut-in at the sand faCe.1,2,3,5,8,10 However, field curves obtained by conventional surface shut-in methods do not always exhibit this characteristic. Such factors as stratification, rock heterogeneities and irregular reservoir geometry can cause the character of a build-up curve to deviate from that predicted by theory for a simple system. In addition, all field build-up data are affected by the methods utilized in obtaining the measurements. Conventional surface shut-in techniques allow two effects to occur which contribute directly to the manner in which pressure builds up at the sand face. The first, afterflow, has been recognized for some time and methods are available for estimating the magnitude of its effect and for correcting its influence.' The second, phase segregation in the tubing column after shut-in, has been reported only recently and appears to exert a considerable influence on the character of field build-up curves.!' In extreme cases, phase segregation can produce "pressure humps" in the early portion of the build-up curve. Such curves have been considered anomalous and have defied analysis. It is the purpose of this paper to discuss the effect of phase segregation in the tubing string during build-up and to present a bottom-hole shut-in technique for obtaining field build-up data which minimizes the influence of afterflow and eliminates phase segregation effects. It will be shown that reliable reservoir information can be calculated from build-up measurements obtained using this method, even though conventional surface shut-in tests on the same wells yield anomalous data. DISCUSSION OF PHASE REDISTRIBUTION EFFECTS Some of the effects of phase redistribution on pressure build-up curves have been described previously by Matthews and Stegemeier" who have presented evidence that phase segregation is responsible for the pres-

By A. G. Weber, K. H. Coats, M. H. Terhune, R. L. Nielsen

Two computer-oriented techniques for simulating the three-dimensional flow behavior of two fluid phases in petroleum reservoirs were developed. Under the first technique the flow equations are solved to model three-dimensional flow in a reservoir. The second technique was developed for modeling flow in three-dimensional media that have sufficiently high permeability in the vertical direction so that vertical flow is not seriously restricted. Since this latter technique is a modified two-dimensional meal analysis, suitably structured three-dimensional reservoirs can be simulated at considerably lower computational expenses than is required using the three-dimensional analysis. A quantitative criterion is provided for determining when vertical communication is good enough to permit use of the modified two-dimensional areal analysis. The equations solved by both techniques treat both fluids as compressible, and, for gas-oil applications, provide for the evolution of dissolved gas. Accounted for are the effects of relative permeability, capillary pressure and gravity in addition to reservoir geomety and rock heterogeneity. Calculations are compared with laboratory waterflood data to indicate the validity of the analyses. Other results were calculated with both techniques which show the equivalence of the two solutions for reservoirs satisfying the vertical communication criterion. INTRODUCTION Obtaining the maximum profits from oil and gas reservoirs during all stages of depletion is the fundamental charge to the reservoir engineering profession. In recent years much quantitative assistance in evaluating field development programs has been provided by computerized techniques for predicting reservoir flow behavior. Because of the spatially distributed and dynamic nature of producing operations, automatic optimization procedures, such as those now in use for planning refining operations, are not now available for planning reservoir development. However, present mathematical simulation techniques do furnish powerful means for making case studies to help in planning primary recovery operations and in selecting and timing supplemental recovery operations. A number of methods have been reported which simulate the flow of one, two or three fluid phases within porous media of one or two effective spatial dimensions>-4 However, applying computer analyses to actual reservoirs have been limited mostly to two-dimensional areal or cross-sectional flow studies for two immiscible reservoir fluids. To obtain a three-dimensional picture of reservoir performance using such two-dimensional techniques, it has been necessary to interpret the calculations by combining somehow the results from essentially independent areal and cross-sectional studies. To the authors' knowledge, the only other three-dimensional computational procedure, in addition to those presented here, was developed by Peaceman and Rachford8 to simulate the behavior of a laboratory waterflood. Two computational techniques which may be used to simulate three-dimensional flow of two fluid phases are described in this paper. The first method, called the "three-dimensional analysis", employs a fully three-dimensional mathematical model that treats simultaneously both the areal and cross-sectional aspects of reservoir flow. The second method, called the vertical equilibrium (VE) analysis", is applicable

By C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By R. C. Earlougher

Jan 1, 1970

By B. S. Gottfried

This paper is concerned with possible transport mechanisms which occur during segregated burning (i.e., burning in an oil reservoir in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas). ,Two analytical models are presented. In the first, burning is assumed to occur at the gas-oil interface, and the dominant mechanism is taken to be the rate at which oxygen diffuses to the interface. The second model assumes that combustion occurs within the clean zone, with the dominant mechanism being the mechanism of fuel deposition from the oil-bearing formation to the clean zone. ,Comparison of the results of the analyses with limited experimental data indicates that the latter mechanism is much more plausible. Both models are idealized representations of the actual process and are somewhat conjectural, however, owing to the lack of more definitive laboratory or field data. INTRODUCTION In recent years there has been considerable theoretical and experimental research directed toward understanding the recovery of crude oil by underground combustion. Although two different recovery schemes have been proposed, known as forward and reverse combustion, it is only forward combustion that shows promise of becoming a practical recovery technique.17 Most of the analyses of the forward combustion process have considered it to be pistonlike; i.e., oil is displaced by a vertical combustion zone which, with its ensuing gases, pushes die oil forward as the flame front is propagated through the reservoir.l-8,ll,l2,14-16 However, information gained in field pilot studies indicates that, at least to a certain extent, forward combustion proceeds in a segregated manner in oil reservoirs.10 An ideal segregated system is one in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas. The location of the combustion zone in such a system is not immediately evident, although combustion is believed to occur either at the gas-oil interface or within the clean zone. In either case, the combustion zone may be horizontal and it moves in parallel rather than in series with the produced oil. The purpose of this paper is to explore the likely mechanisms by which segregated burning can occur, and to determine the sensitivity of the resulting models to changes in various physical parameters. Two governing mechanisms are considered: the case of a horizontal combustion zone at the gas-oil interface with the rate of combustion governed by the rate at which oxygen diffuses to the interface, and the case of combustion within the clean zone with the combustion process controlled by the transport of fuel from the oil-bearing formation to the clean zone. Both models are somewhat hypothetical and involve considerable simplification of the actual process. However, the results allow some preliminary understanding of segregated burning which is useful in the absence of more definitive laboratory or field pilot data.

Jan 1, 1967

By P. J. Closmann

Steam zone growth as a function of time has been calculated for the case of constant rate steam injection into a preheated reservoir. To simplify the calculation a linear temperature profile has been assumed in the cap and base rock at the start of steam injection. The results indicate that at early times augmentation of steam zone growth due to preheating should be greatest. At longer times the steam zone development becomes close to that calculated with no preheating. INTRODUCTION With the increasing application of thermal recovery processes to recover viscous oil, cyclic steam injection has become important in many large-scale projects. In this process, the cycle of steam injection followed by oil production is repeated a number of times. At the beginning of the second and later cycles, steam is injected into a reservoir that has already been heated but that has lost part of its heat both in produced fluids and by conductive heat loss away from the injection zone. In such a case, the temperature level of the injection zone and the temperature distribution of the surrounding rock will affect the growth of the steam zone developed during the subsequent steam injection. Knowledge of the size of steam zone developed is important in determining the amount of oil displaced and the extent of heating in the reservoir. It is also useful to be able to compute the size of the steam zone for cases where thief zones take most of the injected steam. This paper presents a fairly straightforward method of estimating the steam zone developed in a preheated formation, based on certain simplifying assumptions. Some cases of steam injection into a reservoir at its original temperature have already been considered elsewhere.1-3 Because it is difficult or impossible to obtain an accurate representation of the temperature distribution in the reservoir some time after initial heating has taken place, in this work, a linear temperature profile in the cap and base rock is assumed (Figs. 1A and 1B). The linear profile has the advantage of yielding heat flow equations that

Jan 1, 1969

By Robert G. Nisle

The classic theory of pressure build-up in shut-in oil wells as developed by Horner and van Everdingen is based on two-dimensional radial symmetry in the well-reservoir system. Such symmetry does not exist in the case of a well which parfiially penetrates the producing formation. As a result of this lack of symmetry the use of the classic theory in such cases become questionable. In this paper the mathematical theory has been extended to include the case of partially penetrating wells. Numerical solutions illustrating the up-dication of the equations are presented. The effect of partial penetration on pressure build-up is shown by a comparison of synthetic pressure build-up curves derived from the numerical solution of the equation for partially penetrating wells for various degrees of penetration. It is shown that partial penetration is detectable from the characteristic shape of the pressure build-up curve and that formation productivity may be calculated from the pressure buildup data in a manner identical to that described by the classic theory. INTRO DUCTION The use of pressure build-up data on shut-in oil wells is a well-established technique for the measurement of reservoir productivity. The work of Homer' and van Everdingen' is well known. Since the publication of their pioneering work, others"' have elaborated and extended the technique. These authors have dealt exclusively with the case of wells that completely penetrate he producing formation. Two-dimensional radial symmetry is achieved, thereby, and the problem is greatly simplified. This ideal situation is seldom encountered in practice. But, in spite of this, the technique is often employed anyway. The question, therefore, arises: how much confidence can be placed in the results of such calculations when the data are taken on wells that only partially penetrate the producing formation? It is the purpose of this work to answer this question. Only the single-phase fluid case will be considered. van Everdigen' has presented a rigorous solution of the two-dimensional case, but because of mathematical complexity, this solution has not proved to be particularly useful. Horner' utilized Kelvin's point source solution in two dimensions and showed that it is a sufficiently good approximation for oil reservoir studies. The geometry of the well and the producing formation is illustrated in Fig. l. Geometrically the formation may be described as an infinite slab. The infinite slab is defined as the solid (the porous medium in this case) bounded by two parallel planes. Thus, the slab extends to very great distances in both the positive and negative directions of two of the coordinate axes, but is finite in

By M. Prats

Most of the information available on the heat efficiency of hot fluid injection processes, both water and steam, has been obtained from calculated temperature distributions in the pay zone and adjacent formations. The usual approach has been to write the heat balance equations in terms of the temperatures, and then to introduce whatever simplifications are necessary to help obtain an analytical or a numerical solution. In either case, the assumption has often been made that the vertical thermal conductivity in the flooded interval is infinite, so that the temperatures within the flooded interval are then independent of vertical position. Because it was first made by Lauwerier,l and has been used extensively since, this assumption will be referred to as the Lauwerier assumption. Excellent reviews of the literature on the heat efficiency of hot fluid injection processes have been given by Spillette,2 Flock et al.,3 and Ramney.4 Of course, the most significant contributions take into account the effect of a finite vertical thermal conductivity on the vertical temperature profile. The most general analytic expression for the heat efficiency of a hot fluid injection process is that of Antimirov5 who considers injection of a heated incompressible fluid into a reservoir through an arbitrary number of wells. The rate of heat injection into the reservoir, as well as the injection temperature, is an arbitrary function of time. The geometry of the horizontal flow is arbitrary, although the reservoir is considered to be of uniform thickness and properties and to be of infinite areal extent. Heat transfer within the reservoir is by horizontal convection and conduc- tion, and by vertical conduction. In the formations adjacent to the flooded interval, heat transfer is by conduction in any direction. Thus, the available expression for the heat efficiency due to hot water injection has been developed for rather general conditions. This degree of applicability has not been obtained for other thermal recovery processes. By making the Lauwerier assumption we have been able to obtain expressions for the heat efficiency that are essentially independent of the thermal recovery process, be it steam, hot-water, or underground combustion. Clearly, then, the approach followed here is more restrictive than that of Antimirov5 in that it uses the Lauwerier assumption. At the same time less restrictive assumptions are made about the horizontal heat transfer mechanisms, either in the flooded interval or in the formations adjacent to it. The main difference between our approach and that of other investigators is that the heat balance in the pay zone is expressed for the pay zone as a whole rather than for a volume element within it. The Lauwerier assumption is introduced after the problem is developed along these lines as far as possible. The details of the development of the heat efficiency are given in the Appendix, while the main features of the model are discussed in the next section. Results and their implications are discussed later in a separate section. The Model Certain assumptions have been made to determine

Jan 1, 1970

By H. T. Kennedy, S. M. Avasthi

An equation developed for gaseous hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 0.73 percent when applied to 264 natural gas and condensate systems including 2,043 PVT points. Another equation developed for liquid hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 1.12 percent when applied to 346 crude oil systems including 1,759 PVT points. Both equations require composition of the mixture to be expressed as mole fraction of methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide, together with the characteristics of the heptanes-plus fraction in addition to the temperature and pressure. The equations cover wide ranges of the variables involved, and their accuracy is considerably better than that of other available methods. The equations were differentiated to allow calculation of the coefficients of isothermal compressibility and isobaric thermal expansion. (In this paper the coefficient of isothermal com-pressibility and the coefficient of isobaric thermal expansion will be expressed as compressibility and thermal expansion coefficient, respectively.) Equations to calculate these quantities are presented. INTRODUCTION Calculations of reservoir performance for petroleum reservoirs require accurate knowledge of the volumetric behavior of hydrocarbon mixtures, both liquid and gaseous. Compressibilities are required in transient fluid flow problems, and thermal expansion coefficients are important in thermal methods of production. An accurate laboratory investigation of the PVT behavior of each reservoir fluid encountered would be costly and time consuming. For this reason various correlations for predicting fluid properties have been developed and recorded in recent literature. Correlations have been presented in the form of graphs, tables and equations. Since an increasing number of studies are being conducted with the aid of electronic computers, recent efforts have been directed toward development of correla-tions suitable for computer programming. Application of computers permits the use of more complex correlations which otherwise are not feasible. Moreover, methods for predicting reservoir performance, particularly those based on the compositional material balance, depend upon the capability of accurately expressing the molal volumes and other fluid properties as functions of pressure, temperature and composition. The coefficient of isothermal compressibility c is defined by

Jan 1, 1969

By F. J. Fayers, R. L. Perrine

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett' theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy's law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

By K. E. Gray, M. S. Seth

In Part I of this work,1 equations of elasticity were formulated for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure. Equations of elasticity and the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant rate of pressure or a constant rate of flow is maintained at the wellbore. In this paper, the approach presented earlier is extended to finite reservoirs for the cases of (1) steady-state flow, (2) constant pressures at the well bore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Results of this work show that radial and tangential stress gradients are high near the wellbore but diminish rapidt; away from the well; the vertical stress gradient behaves in the same way but is less severe. Radial stresses are compressive or neutral, whereas tangential and vertical stresses may be tensile, neutral or com-pressive, depending upon the boundary conditions, the physical Properties of the system and the radial distance involved (vertical stresses are always compressive in an unbounded system1. For constant boundary pressures, both radial and tangential stresses increase with time whereas they both decrease for a closed outer boundary and constant pressure at the wellbore. The vertical stress decreases with time for both systems. For steady-state systems, radial displacement may be positive or negative, depending upon the dimensions of the system, the pressure differential and the porosity. Radial displacement may be positive or negative for a closed outer boundary but is positive for constant pressures at both boundaries. INTRODUCTION The importance, utility and complexity of a realistic appraisal of the stress state at and local to a wellbore were indicated in Part I. In this paper the analytical approach presented earlier is extended to finite, cylindrical reservoir geometry for the cases of (1) steady-state flow, (2) constant pressures at wellbore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Other than the outer boundary of the reservoir being finite, the physical system and assumptions pertinent thereto are the same as before. The reader may wish to review the mathematical development through Eq. 49 of Part I before proceeding here. SOLUTION FOR GENERAL BOUNDARY CONDITIONS THERMAL STRESSES AND DISPLACEMENTS Since normal stresses do not exist at the free boundaries of a system subjected to thermal stresses, the boundary conditions are

Jan 1, 1969

By A. G. Spillette, R. L. Nielsen

The purpose of this paper is to further the understanding of reservoir response to hot-water injection by desuribing a two-dimensional, mathematical model of the process. Key assumptions are that no gas phase is present, and that the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. The resulting system of three partial differential equations is solved simultaneously through the use of a "leap-frog" application of standard alternating direction implicit methods for the solution of the mass-balance equations and the method of characteristics for solution of the energy-balance equation. The utility of the mathematical model is demonstrated by comparing numerical and analytic temperature distributions for hot-water bank injection and by comparing calculated with observed field behavior. Additional calculations show that hot waterfioding can recover significantly more oil than cold waterflooding, and that a hot-water bank recovers, with less energy input, nearly as much oil us continuous injection. Introduction Consistent field success with hot fluid injection processes requires good reservoir description and a thorough understanding of recovery mechanisms. The latter is fostered by a combination of well designed experimental studies and physically sound mathematical modeling. The purpose of this paper is to further the understanding of reservoir response to hot-water injection by describing a two-dimensional mathematical technique, indicating its validity, and demonstrating its utility for studying the effects of reservoir and operating parameters. Published experimental studies of hot fluid injection processes are few." Even if extensive data were available, two considerations discourage exclusive reliance on laboratory or field work. First, because of the severity of scaling requirements," laboratory results must be interpreted with care. Second, because of the one-shot nature of field experiments — injectivity/productivity tests and pilot floods—results seldom are available over the desired range of operating conditions. These factors emphasize the need for devoting attention to mathematical modeling to complement laboratory and field work. Through the use of mathematical models, scaling uncertainties can be bridged and response can be predicted for unique combinations of reservoir and operating conditions. Numerous mathematical models developed during recent years enable us to calculate temperature distributions, thermal efficiency (or conversely, fraction of heat lost), and oil recovery behavior for hot fluid injection. Ramey" and Spillette' recently have provided reviews of these methods. Most models have concentrated on predicting temperature distributions and thermal efficiency; few have been directed at predicting oil recovery.'."1° Although these techniques are useful for estimating effects of reservoir parameters and operating conditions on process performance, results are limited by the assumptions made and by the methods of coupling independently solved fluid-flow and energy-balance equations. A computer-based method is presented for predicting total reservoir response to hot-water injection that obviates most simplifying assumptions. The model simulates fluid flow and heat transfer in two dimensions within a vertical cross-section spanning the oil sand and adjacent unproductive strata. Known numerical procedures are used to solve the governing partial differential equations. The mathematical model handles the effects of reservoir heterogeneity, gravity, capillarity, relative permeability and temperature-dependent fluid properties. In addition, a wide range of operating conditions can be modeled, including hot-water followed by cold-water injection. Mathematical Description of Hot-Water Injection Key assumptions in the mathematical model for hot-water injection are that (1) no gas phase is present and (2) the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. Relative permeability and capillary pressure are assumed independent of temperature. Darcy's law for multiphase fluid flow relates the velocities of the oil and water phases to the phase pressure gradients by the following relations:

Jan 1, 1969

By G. E. Warner

This paper presents a review and analysis of a highly stratified Burbank sand waterflooding project in Osage County, Okla. Permeability values in this reservoir range from less than 0.1 md to nearly 3 darcys, and 90 percent of the permeability capacity is contained in approximately 26 percent of the reservoir volume. To date, the project has recovered 103 bbl/flood-acre-ft since waterflood initiation, which represents 50 percent of the total preflood cumulative and 12 percent of the original stock-tank oil in place (OSTOIP). The outstanding feature of this project has been the economical recovery of a sizable quantity of oil at water-oil ratios higher than those normally considered prohibitive. Nearly 70 percent of the recovery to date has been achieved at water-cut values greater than 95 percent. An ultimate waterflood recovery of 138 bbl/flood-acre-ft, or 67 percent of the preflood recovery and 16 percent of the OSTOIP, is predicted. The ultimate primary and secondary recovery will be 302 bbl/acre-ft and approximutely 40 percent of the OSTOIP. The actual performance is compared with the predicted performance by both the Stiles and the Dykstra-Parsons methods, and applicabiliry of these methods to this particular reservoir is evaluated. The performance of this project indicafes highly stratified reservoirs can be successfully and economically waterflooded when either reservoir or economic conditions preclude the use of mechanical and chemical means of mdbility control. Introduction The concept of the stratified reservoir has been used by reservoir engineers for many years to describe the vertical variations in permeability within a reservoir. Numerous techniques have been advanced to predict the expected waterflood performance of the stratified reservoir, and in most cases they prove reasonably accurate. The advisability of undertaking a prospective waterflooding project is determined, based on these predictions. This paper presents a review of a waterflood project in a highly stratified reservoir that normally would be considered a marginal prospect, at best, when evaluated with any prediction technique properly adjusted for reservoir conditions. This kind of project can be successful if the predictions are used to design a waterflooding system capable of economically recovering an oil cut of 1 to 2 percent. Experience gained in the early stages of the project's development was profitably applied in the later development of the entire reservoir and could prove beneficial to other projects with similar reservoir stratification or fracture characteristics. The East Burbank pool waterflooding operations have been developed in three separate stages, each of which has been a separate project: Mid-Burbank Unit, Flood A, and Stanley Waterflood. When referred to collectively, these projects will be called the Composite Project. The Skelly Oil Co.'s Barber lease — the only portion of the pool not included in the three projects — is operated as a separate facility with a cooperative injection-well agreement. It contains less than 80 productive acres and would not significantly affect the over-all performance picture.

Jan 1, 1969

By B. H. Caudle, W. J. Bernard

Factors which influence the success or failure of a waterflood can seldom be determined in the laboratory. For this reason pilot waterfloods are initiated in a repreventative portion of the oil reservoir in question. For a pilot flood to predict quantitatively the recovery to be expected in a field-wide waterflood operation, the pilot area must behave as though it were confined (surrounded by similar areas). In this study, laboratory fluid-flow models were used to determine the simplest pilot pattern, for particular conditions of mobility ratio and initial gas saturation, that would behave as though it were confined. Pilot patterns studied ranged in complexity from a single inverted five-spot to a grouping of nine regular five-spots. Only the innermost producing well in each pattern was studied. Model results showed that the optimum number of wells in the pattern depends upon the oil-water mobility ratio and the expected oil-bank size. Unfavorable mobility ratios will, in general, require more wells in the pilot pattern than will favorable mobility ratios. Pilot patterns in reservoirs which contain a dispersed, flowable, free gas saturation will require fewer wells than for the under-saturated case. The single inverted five-spot pattern was found to be unsatisfactory for predicting behavior of fully developed waterfloods. In particular, it is possible that, in reservoirs which contain a flowable, dispersed gas phase, the oil bank will never be observed at the producers due to the large amounts of free gas which continue to be produced with the oil. INTRODUCTION One method which has been used to predict the performance of a waterflood is the pilot flood. The pilot waterflood is a flood which involves only a small cluster of the reservoir wells and is located in a small, representative portion of the reservoir. The object is that oil produced from the pilot can, in some way, be related to the oil recovery to be expected from a field-wide expansion of the waterflood. However, these field pilot waterfloods have often been unreliable in the prediction of oil recovery in a fully developed waterflood. This unreliability has also been demonstrated in several laboratory studies of pilot floods. Some of the investigators have shown that there are situations in which the pilot flood oil production is far too optimistic with respect to the oil recovery in the fully de- veloped flood. Others4-G have shown that the pilot results can also be pessimistic, especially if the pilot waterflood is initiated in an oil reservoir which has been depleted by primary recovery and is at very low pressure. The major reason for this unreliability of pilot water-floods is the migration of fluids into or out of the pilot area. By the well-known method of images, if straight lines can be drawn to represent vertical planes of symmetry in a porous medium which contains pressure sources and sinks (injectors and producers), then these lines are invariant streamlines, or lines across which there is no potential gradient, and therefore no flow. In an actual reservoir, these lines of symmetry can never be established exactly because of reservoir inhomogeneities and irregular reservoir boundaries. However, if the reservoir is relatively large and contains wells in repetitive patterns, these lines of symmetry are commonly assumed to exist for the pattern units sufficiently far removed from the reservoir boundary. Lines of symmetry for the five-spot injection pattern are shown in Fig. 1. Each five-spot unit in this figure can be considered confined with respect to flow across its boundary. In pilot floods this is not the case. The lines of symmetry for the pilot patterns investigated in this study are shown in Fig. 2. It is obvious that the fluid within these pilots is not confined and is therefore able to migrate into or out of the pilot area. Intuitively, one can see that, if more wells are added to the pilot, the innermost unit tends to behave more and more like the confined pattern. However, there is a practical limit to the number of wells which should be placed in the pilot. This limit is usually determined by economic factors. It was the purpose of this study to use laboratory fluid-flow models to determine which of the previously mentioned pilot patterns will force the innermost producing well to behave as it would in a fully developed waterflood. Since fluid migration is influenced by initial saturation conditions and the mobility ratio, these factors were included in the study. The ultimate objective of this study was to develop data which would allow the operator to choose a pilot pattern and operating conditions that will yield a production history which can be applied directly as an estimate of the performance of each production well of the fully developed waterflood. BASIS FOR THE STUDY The basic problem of field pilot floods is the migration of the reservoir fluids into or out of the pilot area. This problem has been the subject of previously reported model studies on pilot floods. These studies have been concerned mainly with the development of arbitrary "correction factors" to be applied to the simple, unconfined pilot systems such as the single five-spot. The correction factors were intended to adjust the production history of the un-

By W. E. Stiles

A method is presented for predicting the performance of water flooding operations in depleted, or nearly depleted, petroleum reservoirs. The method makes use of permeability variations and the vertical distribution of productive capacity. From these two parameters can be calculated the produced water cut versus the oil recovery. Derivations of the mathematical analogy is shown and sample calculations and curves of prediction are presented. Comparison is made of the predicted and actual performance of a typical 5-spot in an Illinois water flood. INTRODUCTION The use of water as a flooding medium in both depleted and "flush" oil reservoirs is gaining greater recognition and acceptance. Many of the shallower fields, depleted by primary production, have been and are being subjected to water injection in order to obtain some part of the large volume of oil remaining after primary production. Some of the earlier water flood installation proved highly discouraging and the value of water flooding was often questioned. Many of these earlier floods were haphazardly selected and developed as little was known of the physical characteristics and contents of the producing formations. The prior evaluation of the flood performance was impossible. During the past decade the development of the required reservoir engineering tools-—core analysis, reservoir fluid analysis, electric logs, fluid flow formulae, etc.—has allowed the engineer to construct and apply the methods which are presently being used to evaluate the economic and mechanical susceptibility of a reservoir to flooding. This discussion will present a method for taking into account the effect of permeability variations in predicting the performance of water floods in depleted reservoirs. PERMEABILITY AND CAPACITY DISTRIBUTION It is generally agreed by most investigators that in a single phase system fluid will flow in a porous and permeable medium in proportion to the permeability of the medium. Producing formations are usually highly irregular in permeability, both vertically and horizontally. However, zones of higher or lower permeability are often found to exhibit lateral continuity. Thus, while structurally comparable stringers in adjacent wells may differ several fold in permeability values, they usually bear resemblance as being part of a general continuous higher or lower permeability section. It is generally agreed that where such stratification of permeability exists, injected water sweeps first the zones of higher permeability, and it is in these zones that "break-through" first occurs in the producing well. It is a basic assurnption of the presently described method that penetration of a water front follows the individual permeability variations as if such variations were continuous from input to producing well. This is admittedly not rigorously true, but can be justified as making possible a simplifying mathematical approach to an otherwise extremely complicated three dimensional flow problem. As a basis for study of the lateral flow of fluids in formations of irregular permeability, the irregularities may be conveniently represented by a permeability distribution curve and a capacity distribution curve. In obtaining these curves, the permeability values, regardless of their structural position in the formation, are rearranged in order of decreasing permeability. If these permeability values so arranged are plotted against the cumulative thickness, a permeability distribu- tion curve is obtained. This curve may then be likened to a "smoothed" permeability profile of the formation. In making comparison between different distribution curves it is convenient to state the permeabilities in terms of the ratio of the actual permeability values to the average permeability of the formation. These ratios termed "dimensionless permeabilities", are used in this paper rather than the permeabilities in terms of millidarcys. The capacity distribution curve is a lot of the cumulative capacity (starting with the highest permeabilities) versus the cumulative thickness. The capacity and thickness are given as fractions of the total capacity and thickness. Mathematically, the capacity distribution is the intergration of the permeability distribution curve. In practice it is convenient to first obtain the capacity distribution curve and derive from it a smoothed dimen-sionaless permeability curve. The method of obtaining the capacity distribution curve is illustrated in the successive column of Table 1, in which capacity and thickness are derived as fractions of their respective totals. If only a small number of permeability values are available, it is generally desirable to smooth the resultant curve. This has been done to give the capacity distribution curve shown in Figure 1. The differentiation of the capacity distribution curve to obtain the permeability distribution curve is shown in Table 2. Here, the capacity values are read from the smoothed curve at intervals of cumulative thickness, and the increments of capacity are divided by the increments of thickness to obtain the dimensionless permeability, K. Due to this stepwise procedure of calculation these premeability values must be plotted at the midpoints of the successive increments of thickness. The curve so obtained from these data is shown in Figure 1. The total area under the K' curve is equal to unity.

Jan 1, 1949

By W. E. Stiles

A method is presented for predicting the performance of water flooding operations in depleted, or nearly depleted, petroleum reservoirs. The method makes use of permeability variations and the vertical distribution of productive capacity. From these two parameters can be calculated the produced water cut versus the oil recovery. Derivations of the mathematical analogy is shown and sample calculations and curves of prediction are presented. Comparison is made of the predicted and actual performance of a typical 5-spot in an Illinois water flood. INTRODUCTION The use of water as a flooding medium in both depleted and "flush" oil reservoirs is gaining greater recognition and acceptance. Many of the shallower fields, depleted by primary production, have been and are being subjected to water injection in order to obtain some part of the large volume of oil remaining after primary production. Some of the earlier water flood installation proved highly discouraging and the value of water flooding was often questioned. Many of these earlier floods were haphazardly selected and developed as little was known of the physical characteristics and contents of the producing formations. The prior evaluation of the flood performance was impossible. During the past decade the development of the required reservoir engineering tools-—core analysis, reservoir fluid analysis, electric logs, fluid flow formulae, etc.—has allowed the engineer to construct and apply the methods which are presently being used to evaluate the economic and mechanical susceptibility of a reservoir to flooding. This discussion will present a method for taking into account the effect of permeability variations in predicting the performance of water floods in depleted reservoirs. PERMEABILITY AND CAPACITY DISTRIBUTION It is generally agreed by most investigators that in a single phase system fluid will flow in a porous and permeable medium in proportion to the permeability of the medium. Producing formations are usually highly irregular in permeability, both vertically and horizontally. However, zones of higher or lower permeability are often found to exhibit lateral continuity. Thus, while structurally comparable stringers in adjacent wells may differ several fold in permeability values, they usually bear resemblance as being part of a general continuous higher or lower permeability section. It is generally agreed that where such stratification of permeability exists, injected water sweeps first the zones of higher permeability, and it is in these zones that "break-through" first occurs in the producing well. It is a basic assurnption of the presently described method that penetration of a water front follows the individual permeability variations as if such variations were continuous from input to producing well. This is admittedly not rigorously true, but can be justified as making possible a simplifying mathematical approach to an otherwise extremely complicated three dimensional flow problem. As a basis for study of the lateral flow of fluids in formations of irregular permeability, the irregularities may be conveniently represented by a permeability distribution curve and a capacity distribution curve. In obtaining these curves, the permeability values, regardless of their structural position in the formation, are rearranged in order of decreasing permeability. If these permeability values so arranged are plotted against the cumulative thickness, a permeability distribu- tion curve is obtained. This curve may then be likened to a "smoothed" permeability profile of the formation. In making comparison between different distribution curves it is convenient to state the permeabilities in terms of the ratio of the actual permeability values to the average permeability of the formation. These ratios termed "dimensionless permeabilities", are used in this paper rather than the permeabilities in terms of millidarcys. The capacity distribution curve is a lot of the cumulative capacity (starting with the highest permeabilities) versus the cumulative thickness. The capacity and thickness are given as fractions of the total capacity and thickness. Mathematically, the capacity distribution is the intergration of the permeability distribution curve. In practice it is convenient to first obtain the capacity distribution curve and derive from it a smoothed dimen-sionaless permeability curve. The method of obtaining the capacity distribution curve is illustrated in the successive column of Table 1, in which capacity and thickness are derived as fractions of their respective totals. If only a small number of permeability values are available, it is generally desirable to smooth the resultant curve. This has been done to give the capacity distribution curve shown in Figure 1. The differentiation of the capacity distribution curve to obtain the permeability distribution curve is shown in Table 2. Here, the capacity values are read from the smoothed curve at intervals of cumulative thickness, and the increments of capacity are divided by the increments of thickness to obtain the dimensionless permeability, K. Due to this stepwise procedure of calculation these premeability values must be plotted at the midpoints of the successive increments of thickness. The curve so obtained from these data is shown in Figure 1. The total area under the K' curve is equal to unity.

Jan 1, 1949