Search Documents

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

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

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

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

Jan 1, 1966

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

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

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

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

Jan 1, 1966

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

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

By Homer N. Mead

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

Jan 1, 1957

By W. A. Burns

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

Jan 1, 1970

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

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

Jan 1, 1966

By A. M. Sarem, J. M. Campbell

Molecular refraction is introduced as a new and improved third parameter for prediction of the PVT behavior of hydrocarbon systems. This parameter, characterizing the complex as well as the pure hydrocarbon systems (1) is a direct measure of London dispersion forces which affect the PVT behavior; (2) can be measured directly, accurately and easily for complex as well as pure hydrocarbon systems; (3) varies significantly between the reference substances thus allowing accurate interpolation of PVT properties; and (4) has a fundamentally sound basis. Molal average molecular refraction was used with the pseudocritical mixing rule proposed by Leland and Mueller to predict the molal volume of liquid and vapor for 84 published data points on nC4-nCI0, Cl-nC5, Cl-nC4-nC10 systems in the region of coexisting vapor-liquid and at the critical point. The same technique was applied to predict the densities of liquid hydrocarbon mixtures in the two-phase and single-phase regions near the critical point for seven different systems each having different properties of the heptanes-plus fraction, but all simulating the fluids of either gas-condensate or volatile oil reservoirs. A modification of the Leland-Mueller mixing rule for the determination of pseudocritical properties is also presented which yields more accurate liquid density predictions than the original form when used with the proposed third parameter. INTRODUCTION The PVT behavior of multicomponent systems has received widespread attention in the literature in recent years. A large portion of the reported work has stemmed from attempts to overcome the inherent difficulties encountered when using the standard compressibility-factor chart which shows 2 as a function of reduced temperature T, and and reduced pressure p,. This chart assumes that z, is constant for all systems. The error in this assumption is not pronounced for pipeline-type gases high in methane and low in heavy ends, particularly when contained at moderate pressures. For such gases the error in the calculated z is usually less than 2 to 3 per cent. However, the reliability of this approach becomes somewhat nebulous for gas condensate and volatile oil systems at high pressure. One logical approach has been to recognize the variation of 2,. This idea was first suggested by Meissner and Seferianl and used by Lyderson, Greenkorn and Hougen2 to prepare extensive tables of thermodynamic properties. Other third parameters have been proposed by Stockmeyer, Kihara, Riedel, Lightfoot, Bloomer and Pitzer. These were tested by Satter and Campbel13 for many gas mixtures and it was shown that they were inter-related. All of these third parameters have a common fault — they cannot be directly evaluated for the heptanes-plus fraction. Consequently, one necessary requisite of an improved method for predicting the density of multicomponent hydrocarbon systems, was development of a third parameter that could be evaluated for that fraction. PREVIOUS MULTICOMPONENT DENSITY PREDICTION METHODS Sage and Lacey4 suggested a method of using partial molal volumes for computing the density of hydrocarbon gas mixture for pressures to 3,000 psi. In this technique, they approximate the complex multicomponent gaseous mixtures by a quarternary system of methane, ethane, propane and butane. Sage, Hicks and Lacey5 suggest a method of using partial molal volumes for computing the density of hydrocarbon liquids. The results agree within 3 per cent of the experimental values. However, the composition range is limited to about 10 per cent by weight of methane. Consequently, this correlation does not cover the low-molecular-weight liquid similar to natural gasoline and high-pressure gas condensates.

Jan 1, 1966

By C. S. Matthews, P. Hazebroek, H. Rainbow

It ha been suggested that lormation fractures created by well stimulation treatments will adversely affect sweep-out efficrency in injection operations. Fluid-flow model studies involving vertical fractures of various lengths and fluid systems of various mobility ratios have been carried out to study this subject. In addition, limited data have been obtained on one model containing a horizontal fracture. It was found that relatively long and highly conductive fractures (not generally obtained fracturing operations) were required to affect the sweep-out efficiency substantially. In a given case in the field an approximate distinction can be made between the presence of long conductive fractures and shorter or less conductive ones. This is done with pressure build-up analyses along with data on the relationship of fracture length and conductivity to well productivity. This type analysis shows that usually the fractures induced are either short or of limited conductivity and therefore do not damage sweep-out efficiency. INTRODUCTION Improved well productivity and in-jectivity can frequently be exploited in injection operations. Higher total throughput can yield improved economics. On occasion, the achievement of an increased productivity or injectivity in specific wells can bring about a more uniform sweep of the reservoir. Higher rates can be exploited particularly in water floods of "depleted" reservoirs where a rapid "fill-up" is desired and where low pressures contribute to low well productivity. The creation of fractures local to the wellbore is an excellent means for achieving these objectives. Even though fracturing has been employed in some floods with success,1"3 there still seems to be some reluctance to employ this tool for fear of undue damage to the flood pattern and ultimate recovery. We therefore need to examine the influence which fractures of varying length may have on flood performance and then determine the length of fractures which obtain in the field with conventional fracture treatments. A substantial influence of fractures on the recovery obtained at breakthrough of the injected fluid has been presented by Crawford and Collins for equal fluid mobilities for the line-drive pattern.4,5 In addition, the influence which a fracture has on the production performance after breakthrough and on the ultimate recovery warrants consideration in reaching a conclusion concerning the use of induced fractures in flooding operations. This report presents this type of data for the five-spot injection pattern for several fluid mobilities. In arriving at some conclusion concerning the fracture lengths obtained in field operations we must examine the performance characteristic most affected by the fracture. This is the change in the flow system as reflected in the change in productivity and pressure build-up behavior. A study of these changes is also presented in this report. RESERVOIR ANALOGS The X-ray shadowgraph technique, employing miscible displacement in porous models, has been used in the study of the influence of fractures on pattern sweep-out efficiency. The X-ray shadowgraph procedure is described in detail in an earlier report8. Fractures were represented by leaving the proper portion of the model surface exposed to either injection or production. This assumes the fracture resistance to be negligible compared to that of the formation. Actually the flow resistance in propped fractures obtained in the field is sometimes not negligible so that the results with this model indicate the maximum influence of the fracture. Two types of models were necessary to represent vertical fractures in a five-spot flood. These pattern elements are illustrated in Fig. 1. In studying the influence of the horizontal fracture, only one well spacing length to thickness ratio of 50 was used. The pool unit of a Carter Electric Analyzer was used in studying the influence of fractures on productivity and build-up behavior. The square drainage system of a well was represented by a network of 576 elements of equal volume and resistance. One-fourth of this drainage area was studied as a model unit using a network of 144 condensers and resistors. Vertical fractures were represented by a shunt directed from the well perpendicular to the drainage boundary. Horizontal fractures were represented by a circular shunt. Resistance of a shunt was varied in representing different conductivities for the fractures. FRACTURE DIRECTION AND LENGTH The direction at which a vertical fracture extends into the formation and its length and conductivity influence the sweep behavior. The two

By Michael P. Robinson

It has been suggested that streaming potentials are not nomlally logged because the streaming potentials known to be generated across mud filter cakes are substantially cancelled by streaming potentials generated in contiguous shale beds. Laboratory measuments of streaming poterztials acrons hand specirnens of shale have been adduced as evidence to support the suggestion. The experiments show that streanzing potential depends on the resistivity of the equilibrating solution and on pressure drop. Meizsurenlents were made perpendicular to the bedding planes of the shales. The mechanism of streaming potential generation in shale which the suggested cancellation theory iniplies has been critically examined. Conclusions reached are tested by comparing self potentials logged before and immediately after changing the resistil~ity (but not the weight) of mud in a borehole. The resultant SP change is found to equal the theoretical electrochetnical potential change, i.e., cancellation of .strearning potential is independent of mud resistivity. Since streaming po/entials generateal across mud filter cakes depend critically upon mud reistivity, shale .strearlling potential in situ must he equally sensitive. Using plausible values for the pemleahility, porosiy arzd compressibility of shale, the van Ererdingerl-Hurst approach shows that the pressure drop between horehole pressure and formation pressure in a radially honzogeneous shale is unlikely to coincide with the zone of filtrate invasion. If, perchance. this coincidence occurs with one mud a chatzge 10 a second must result in non-coincidence, i.e., a .YP logged in the second mud should show a streaming potenlial. The dilernrna may he overcome by postrrlating: (I) perrneabilities of shale7 parallel to their beddine planes greatly exceed their cross-bedding permeabilities; (2) a resistivity-sensitive mud cake of low permeability forms either on the face of, or probnhly in the lanlinations of, many shales. These postulates are examined with reference to the pres.sure and SP histories of the McClosky and Aux Vases zones in the Illinois Basin and some srcpport for them found. It is concluded that a definitive solultion to the streaming potential problem does not exist. INTRODUCTION Recently Gondouin and Scala' published a contribution to the streaming potential problem which is notable particularly for the elegance of its experimental techniques. The conclusions drawn by Gondouin and Scala from their studies amplified those previously suggested by Schehck'. In essence they suggested that the streaming potentials known to develop across mud filter cakes" are not normally observed on SP logs because they are largely cancelled by streaming potentials generated in adjacent shales. This suggestion appears so reasonable that many workers in the logging field will be prone to accept it. We incline strongly to the view that the suggestion is basically correct. Nevertheless, it is important to realize that the cancellation hypothesis in the form put forward by Gondouin and Scala does not appear to meet all the known experimental facts. It is the principal intention of this note to call attention to what appears to be a flaw in the relevancy of the experimental evidence adduced by Gondouin and Scala and to suggest an alternative hypothesis to remove it. The consequences of the new hypothesis are compared with field data which, we believe, are unique. A second but not less important purpose of this note is to emphasize that even the new hypothesis does not fully explain all field observations. Hence we conclude that further study of the streaming potential problem is required before it can be said that a definitive solution to it has been attained. The bare experimental bones of laboratory data on streaming potentials appear to be as follows. 1. A streaming potential is developed across a mud filter cake. This streaming potential is virtually independent of cake thickness but depends on the mud resistivity. The relationship is of the form, 2. A streaming potential is developed across a shale. This potential is independent of shale thickness but depends on the resistivity of the permeating solution if the shale has been first equilibrated in this solution. The relationship is of the form, The relationship in (1) was determined first by Wyllie3, confirmed by Gondouin and Scala and qualita-

By H. Kazemi

An ideal theoretical model of a naturally fractured reservoir with a uniform fracture distribution, motivated by an earlier model by Warren and Root, has been developed. This model consists of a finite circular reservoir with a centrally located well and two distinct porous regions, referred to as matrix and fracture, respectively. The matrix has high storage, but low flow capacity; the fracture has low storage, but high flow capacity. The flow in the entire reservoir is unsteady state. The results of this study are compared with the results of the earlier models, and it has been concluded that major conclusions of Warren and Root are quite substantial. Furthermore, an attempt has been made to study critically other analytical methods reported in the literature. In general, it may be concluded that the analysis of a naturally fractured reservoir from pressure transient data relies considerably on the degree and the type of heterogeneity of the system; the testing procedure and test facilities are sometimes as important. Nevertheless, under favorable conditions, one should be able to calculate in-situ characteristics of the matrix-fracture system, such as pore-volume ratio, over-all capacity of the formation, total storage capacity of the porous matrix, and some measure of matrix permeability. INTRODUCTION The analysis of flow and buildup tests for obtaining in-situ characteristics of oil and gas reservoirs has received considerable attention in the past decade. Most of the available techniques result in reliable conclusions in macroscopically homogeneous reservoirs or in the homogeneous reservoirs with only certain types of induced and/or inherent heterogeneity (such as wellbore damage, etc.). In general, the greater the degree of heterogeneity, the less the reliability of the information deduced from the pressure transient data. A commonly encountered heterogeneous system is a naturally fractured reservoir where two types of distinct porosities occur in the same formation. The region containing finer pores may have high storage and low flow capacities. This is called the matrix. The remaining region may have high flow capacity with low storage. The latter region is generally the set of interconnecting fractures and fissures of the rock, and for this reason it is called the fracture. Ordinarily, we wish to obtain the permeability and porosity of each region and perhaps the frequency of the fracture distribution in a reservoir. Such information is necessary for reservoir engineering. Other information, such as wellbore damage, will be useful in evaluating possible remedial work for such fields. Few authors have suggested theories to aid in calculating the in-situ characteristics of a naturally fractured reservoir similar to what we have described above. Pollardl suggested that a naturally fractured reservoir contained three distinct regions: a damaged or an improved region surrounding the wellbore, and the two remaining regions the same as described earlier. He suggested that the flow was taking place from the tight matrix into the highly conductive fractures, then into the we wllbore region and finally into the well column. He concluded that the average pressure buildup in each of these distinct regions can be expressed approximately in terms of an exponential decay function of time. He also hypothesized that the decay coefficients for each of these functions were significantly different from each other; consequently, each of these functions became dominant in turn, in the process of pressure buildup. Thus, by a proper graphic plot of the logarithm of wellbore pressure differences vs buildup time, each of these functions could be determined (see Fig. 1). Pollard suggested methods of determining the wellbore damage and fracture volume. Later Pirson and Pirson2 extended

Jan 1, 1970

By P. B. Baxendell

The performance of a water-drive reservoir having a gas cap depends primarily on the movement of the gar-oil and oil-water contacts. The movement of the contacts during production depends in turn on fluid withdrawals and how the reservoir pressure changes as fluids are produced from the reservoir; that is, on how effectively the aquifer maintains pressure by replacing withdrawals. inasmuch as pressure changes, fluid withdrawals, contact positions, and produced gas-oil and oil-water ratios are interdependent, the analysis and prediction of the performance of a reservoir produced in a given way must take into account this interdependence throughout depletion. This paper presents an analysis which, within the limitations of the assumptions made, yields an engineering approach to predicting future performance based on reservoir pressure and production history. The most significant assumption is the method of extrapolation of future gas-oil and water-oil ratios. The extrapolation procedure smoothly increases both ratios to preselected values as the remaining oil column undergoes a specified decrease in thickness. This preselection is made on the busis of previous field experience in depletion of similar reservoirs under similar conditions. For computing future performance a vdumebic balance L combined with the differential equation defining pressure distribution in the aquifer to obtain positions of water-oil and gas-oil contacts. From these positions are extrapolated produced water-oil and gas-oil ratios. Reservoir performance can be investigated when oil production rates are dependent upon various factors including the performance of the reservoir itself. Examples of practical application of the procedure are included. INTRODUCTION To predict performance of water-drive reservoirs with gas caps and thin oil columns, it is necessary to describe the motions of the fluids within the reservoir during the entire production period to depletion. These motions depend on the pressure changes and on the withdrawal of oil, gas and water. Production of oil from the oil zone primarily causes water to move in to take its place; production of gas from the cap tends to cause oil to migrate into its place. In addition, the volumes of oil and gas remaining in the reservoir depend on changes in the pressure, since a decrease in pressure causes fluid expansion, gas liberation and oil shrinkage. A method of relating the future pressure to total withdrawals is used to describe the motions of fluids within the reservoir under conditions arising in possible modes of production. The analysis is based on two relationships and will reduce the problem to one amenable to digital computation. The first of these concerns the dependence of the pressure distribution in the entire aquifer furnishing the water drive upon the total reservoir withdrawals. This dependence is dictated by the permeability distribution and the extent of the aquifer; at present, such information is most readily obtained from the performance history by means of the resistancecapacitance reservoir analyzer. The second relation involves withdrawals, pressure in the reservoir, and movement of oil, gas, and water within the reservoir. The analysis is subject to certain simplifying assumptions that are necessary to permit solution of the problem. A comparison of the methods of this paper with those presently practiced is pertinent. One method of analysis is to use the reservoir analyzer to predict reservoir behavior based on aquifer characteristics determined from production history. Inasmuch as the total withdrawal rate depends upon the gas-oil and water-oil ratios, which depend in turn upon, among other things, the positions of the gas-oil and water-oil contacts, these ratios must be assumed in advance. In this paper these ratios are, instead, related to the computed positions of gas-oil and water-oil contacts. Thus, our method can be applied to problems in which the oil production rate is limited by produced gas-oil ratio, or to problems in which it is desired to determine the variable gas injection rate that will maintain the gas-oil or water-oil contact stationary. These problems cannot be worked satisfactorily on the analyzer. Another method presently used is based on a paper by Hurst.' Since his procedure is dependent upon using the solution of the heat flow equation, which requires constant permeability within the aquifer, and our procedure recognizes variations of permeability in

By J. J. Arps

n approach to the study of the electrical properlies of drilling muds, their cakes and their filtrates was made by treating the drilling fluid as a porous medium. Lluring this study it was found that the weight pararneter was omitted from previous investigations on electrical properties. Patnode in 1949 made a similar approach to ,his investigation. His work showed that there was a variation of resistivity relationships with per cent Puid present in a given suspension. This present study ex-tends Patnode's investigation to one in terms of the solids fraction of the drilling fluids. Previous investigations were made on low weight muds where the effect of the rolids fraction was obscured. The following equationr were fitted to data taken from 94 field muds and 47 laboratory muds (see Tables I and 2). The parameter, C, is a complex function of the solidr content, or mud weight, of a drilling fluid. It was also found that for all drilling muds examined. A nomograph is presented combining Eqs. I and 2. INTRODUCTION It is generally recognized that the materials present in drilling muds affect the relationship between mud filtrate and mud resistivities. However, most electric log interpretations proceed with the assumption that the mud resistivity and temperature are the only important parameters needed to obtain a filtrate resistivity or a mud cake resistivity. Lamont1 suggests that only average relationships exist between mud cake, mud filtrate and mud resistivities. He indicates that there is limited accuracy when these intercorrelated factors are used since no other measure-able variable is apparent. A previous paper by Patnode indicates that drilling fluids may be considered to be a system of fluids in solid particles which may be correlated electrically in a manner similar to rock solids. 1 he equation for the resistivity factor. F = (per cent fluid) is satisfied for a given clay type in water, but the value of m varies as different materials are used in the water. The resistivity factor is not constant for a given amount of clay. As is suspected it varies as the solids fraction conducts electrical current. Since both weighting material and clay solids are present in commonly used drilling fluids, it is expected that both materials contribute to the variation of the m factor in Eq. 3. An extension of Patnode's work is presented in this paper to show some of the interrelations which exist between the resistivities of mud cake, mud filtrate and mud. Use of these data enable a more accurate calculation of mud filtrate and mud cake resistivities if only the mud resistivity and weight are measured. EXPERIMENTAL The measurements of mud weight, fluid and solids resistivity, alkalinity, and chloride content of mud samples were made in this research. Mud samples representing drilling fluids from 94 Louisiana and Texas wells (Table 1) and 47 laboratory-constructed systems (Table 2) were used in the experimentation. The measurements of the mud properties were made both in the field and in the laboratory. The tests on the properties of weight, alkalinity and chloride ion were standard tests as described by API specifications. The resistivity measuring devices were constructed in the laboratory for this study. For details on the resistivity equipment see Fig. 1. The fluid cell

By A. Timur

Fluid flow properties of porous media have long been of interest in such varied disciplines as geology, geophysics, soil mechanics, and chemical, civil, and mechanical engineering. This interest has resulted in numerous models of porous media that have been proposed, tested, and found to be useful under, at best, only special circumstances. A critical review of most of these models is given by Scheidegger. 1 The application of each new tool to the study of porous media has helped to test some of the existing theories and to form new ideas based on the new parameters being measured. In the application of nuclear magnetic resonance (NMR) techniques to the study of the properties of fluids in porous media, the theoretical studies of Korringa et al.2 resulted in a model for the relaxation of spin polarization of protons in a hydrogenous fluid in a pore of a solid. Seevers3,4 verified this model, established a method for measuring the surface-to-volume-ratio distribution in a porous medium, and proposed a technique for determining specific permeability in sandstones. In the present investigation, the uses of NMR methods are extended and a simple model of porous media is developed from an analysis of pulsed NMR measurements. In this model, the pore spaces of a porous medium are divided into three groups on the basis of their surface-to-volume ratios. To evaluate predictions of the NMR model, laboratory measurements of spin-lattice relaxation time, porosity, permeability, and residual (irreducible) water saturation were conducted on more than 150 sandstone samples obtained from four oil fields in North America. Applications are discussed for estimating the volume of movable fluid, ?p, and the specific permeability, k, of sandstone samples. The former parameter, which can be considered as producible porosity, is expressed by ?p = [1 -- (Swr/100)]? , .... (1) where + is porosity in percent of bulk volume and Swr is the residual (irreducible) water saturation in percent of pore volume at a capillary pressure of 50 Psi. Fast, nondestructive laboratory methods are described for determining porosity, movable fluid, and permeability of sandstone samples from measurements with a pulsed NMR (spin-echo) apparatus. Procedures are given for adapting these methods to samples of irregular shape and small size. Theory General In the Korringa, Seevers, and Torrey2 (KST) model the observed decrease in the spin-lattice relaxation time, T1, of protons of a hydrogenous liquid contained in the pore spaces of a porous solid is attributed to an increase in the correlation time for the random motion of the water molecules and to the presence of paramagnetic centers at the liquid-solid interface. This KST model predicts the observed relaxation time of a hydrogenous liquid in a pore of a solid through

Jan 1, 1970

By John C. Martin

The discovery of a new gas reservoir demands that the planning of a sottnd well-spacing program be initiated early in the development stage. It is the purpose of this discussion to illustrate by actual field examples the application of basic well-spacing principles, previously developed for oil reservoirs, to the problem of well spacing in natural gas fields. These studies are presented for the field use of geologists and engineers who are concerned with the initial planning of the proper development of the newly discovered gas reservoir. INTRODUCTION The phenomenal growth of a vigorous natural gas industry emphasizes the increasing importance of natural gas as a source of energy, fuel, and raw materials to our nation's economy. Since 1945 marketed production of natural gas for the U. S. has increased 21/2 times to a record high of 10.6 trillion cu ft during 1957. As major participants in the gas industry, we share an added interest to develop and produce our natural gas reserves with constantly improved efficiency. The subject of well spacing is vitally important to the gas industry, for the well itself plays a significant role in the development of the natural gas reservoir and in control of the recovery process. Maximum utilization of wells is an integral part of sound conservation practices. The discovery of a new gas reservoir demands that careful choice of well location and well spacing be made and that the planning of a sound well spacing program be initiated early an the development stage. With the drilling of the first development well, efforts of the geologist and engineer must be directed toward acquisition of adequate technical evidence upon which a firm recommendation for a spacing program may be based. With this technical appraisal as a foundation, operators and state regulatory agencies jointly can go far in providing a framework for sound development of gas fields to achieve a program of conservation that avoids the unnecessary well. WELL-SPACING CONCEPTS Through laboratory and field investigations of the mechanism of the recovery of oil and gas, of fluid behavior, and of effective control of reservoir and well, a crystallization of ideas regarding reservoir behavior has emerged as a well-developed technology. Associated with a better understanding of the fundamental principles underlying reservoir and well behavior has been the growth of concepts concerning the role of wells and their spacing in the development and operation of an oil or gas reservoir. In addition to serving as outlets for the withdrawal of fluids from the reservoir, wells are recognized as having two other important functions: (1) providing access to the reservoir to obtain information concerning the characteristics of the reservoir and its fluids, and (2) serving as a means by which the natural or induced recovery mechanism may be effectively controlled. Beyond a minimum number of wells required to fulfill these two functions, additional wells will not increase recovery. With particular emphasis upon well spacing in oil reservoirs, many studies of the well spacing-recovery relationship have evolved the concept that the ultimate oil recovery is essentially independent of the well spacing.' These fundamental concepts are no different when regarding the role of wells and their spacing in the natural gas reservoir. They are equally applicable to the consideration of well spacing in gas reservoirs. For the gas reservoir, the problem of well spacing then revolves around the question of drainage and the degree or extent to which a well may drain gas from its surrounding reservoir environment. Theoretical and Experimental Work During the past 30 years, theoretical and experimental work carried on to study the physical principles involved in the flow of fluids through porous media has shed light upon the matter of drainage. Fundamental mathematical equations have been derived to describe the mechanism of flow of oil and gas through porous rocks. With the recent advent of high-speed digital computers, attempts have been made with mounting success to develop solutions, employing numerical techniques, to mathematical expressions that describe more rigorously the physical behavior and mechanism involved in the unsteady-state flow of compressible fluids, such as a gas, through porous rock. In 1953, Bruce, Peace-man, Rachford and Ricez published a stable numerical procedure for solving the equation for production of gas at constant rate. The results of these calculations are significant with respect to this matter of drainage, for they indicated (1) that depletion of the gas reservoir resulted in a drop in pressure at the extremity of the

By John P. Dolan, Charles A. Einarsen, GiIman A. Hill

his paper discusses how the following formation characteristics can be determined mathematically from drill-stem test pressure charts: true formation pressures, effective permeability of the entire section tested, well productivities, wellbore damage, and possible detection of barriers (faults, pinchouts, etc.). This paper also presents a practical method for im-mediate determination of eflective permeability and well-bore damage from successful double shut-in pressure tests. A list of recommendations for improving the reliability of drill-stem rest pressures is also presented. INTRODUCTION A drill-stem test is a temporary completion designed to sample the formation fluid and to establish the possibility of commercial production. Early pressure recording devices were used merely to verify proper operation of the testing tool. Until recently the accuracy of the pressure gauges has been insufficient for any reliable quantitative use of the recorded pressures. In view of the need for more reliable formation evaluation and as a result of the recent interest in exploration work employing the concept of hydrodynamic en-trapment,1,2,3,4 better pressure recording gauges are now in use. These devices can record pressure within 1 per cent above 1,000 psig and can detect differential pressures as low as 1/2 psig. In addition to formation pressure, several other reservoir characteristics can be determined from DST charts; namely, well productivity, formation permeability, wellbore damage, and the possible existence of bar- riers (faults, pinchouts, changes in permeability, etc.). This paper presents a practical method to interpret DST pressure charts for formation pressures and many other reservoir properties, a method that has been developed in analyzing approximately 4,000 DST charts during the last five years. The techniques used are a composite of published articles on drill-stem testing5,6, together with well-known pressure build-up analysis methods. THEORY It has been shown7,8 that the following equation may be used for analysis of pressure build-up curves: Po - Pw = 162.6 qu/kh log ( to + ?/?) .....(1) When this equation is applied to the curves obtained in drill-stem testing, the assumptions and boundary con-ditions are more nearly realized than in conventional flow and build-up tests on producing wells. Zak and Griffen9 have recently discussed in detail the use of this equation in analyzing DST charts. One of the problems with DST curves is the lack of reservoir data for precise analysis. Therefore, it is necessary to develop empirical rules and field methods for analyzing DST charts in quantity. For this reason, the empirical methods presented in this paper have been developed, and their derivation is found in Appendices A and B. METHOD USED FOR ACCURATE READING OF DRILL-STEM TEST CHARTS In order to apply the pressure build-up theory to DST charts, it is necessary to obtain a digitized expression for the pressure and time data recorded by the pressure gauge. These data may be either provided by the service company or interpolated by projecting a

Jan 1, 1958

By J. R. Jargon, H. K. van Poollen

Non-Newtonian fluids may be injected into a reservoir during secondary recovery operations. The non-Newtonian fluid used in this work is a power-law type of fluid that is, the viscosity of the fluid decreases as the flow rate or rate of shear increases. This paper presents equations for steady-state linear and radial flow of such fluids, transient behavior results from a finite difference model of a radial system, and transient behavior results from a field test. The equations that describe the flow of a non-Newtonian fluid are non-linear and are solved numerically. Finite difference solutions are presented as curves of dimensionless pressure drop at the wellbore us dimensionless time for a constant injection rate. Solutions were obtained for 5-percent, 10-percent and 100-percent PV of a non-Newtonian fluid for injection rates of 1, 10, 100 and 1000 cc/sec and for a 5 percent PV of non-Newtonian fluid located at T = Tw' 3, 10, 20, 50 and 100 ft for a flow rate of 1 cc/sec. The buildup curves do not exhibit a stmight-line portion as is the case for Newtonian flow through porous media. Correlations also are shown for the productivity index vs rate for the computer model study and the field tests. INTRODUCTION During various secondary recovery operations non-Newtonian fluids are injected. Such fluids, in general, include thickened water and gelled fluids. The term "non-Newtonian" implies that the viscosity is not only dependent upon temperature and pressure, but also on the rate of shear that is applied to the fluid. For example, water, which is a Newtonian fluid, will have essentially the same viscosity no matter what rate of shear is applied. A pseudo-plastic fluid exhibits a decreasing viscosity when higher rates of shear are applied; a dilated fluid has an increasing viscosity with increasing rates of shear. The objective of studies performed and described in this report is to obtain relationships and mathematical and empirical descriptions of the flow of non-Newtonian fluids through porous media. Simple mathematical relationships, computer studies that include the unsteady-state behavior of such fluids, and field tests are used. PREVIOUS LABORATORY STUDIES Laboratory studies have been performed by several investigators. The fluids used in one investigationl were surfactant-stabilized dispersions of water in hydrocarbons. Its Fig. 6 is reproduced as our Fig. 1. This effective-viscosity vs frontal-velocity diagram shows that this fluid is of the power-law type over a fairly large range of frontal velocities. Gogartyl developed an equation for the effective viscosity as a function of shear rate that reduces to Eq. 2 for frontal velocities greater than about 2 ft/D. The fluid characteristics used in the present study are similar to those reported by Gogarty. Christopher and Middleman2 performed experiments with power-law fluids flowing through pipes packed with glass beads. They used a modified Blake-Koteny

Jan 1, 1970

By Stephen G. Dardaganian

Pore volume compresibilitieu measured in the laborutory on core .samples are reported for typical reservoir .sarrrl.rtorres at reservoir pressure. These cornpressibilities ore different for each reseriloir sample and cannot be correlated to porosity. The compressihilities are a1so functions of presslire. In recent years several authors have pointed out the importance of including pore volume compressibility of the reservoir rock in certain reservoir engineering calculations. Hall' showed that in calculating the hydrocarbon volume of an undersaturated reservoir from the production per unit change in reservoir pressure, the neglect of rock compressibility could, in the extreme case of low porosity rock, lead to results in error by a factor of two. Hawkinsl' and Hobson and Mrosovsky showed that both rock and interstitial water compressibilities must be included to achieve satisfactory accuracy in this calculation. Although rock compressibility is an important reservoir property in these calculations, very few data are available. Hall' made compressibility measurements on a few samples in a small pressure range. No other direct pore volume compressibility measurements have been reported. Geertsma' discussed the elasticity of reservoir rock but did not present any new data. In this paper, rock compressibilities suitable for reservoir calculations are reported for a variety of sandstones typical of U. S. reservoirs. The measurements were made through a range of pressures which is also typical of reservoirs. EXPERIMENTAL PROCEDURE Plugs 1 in. in diameter and 2- to 3-in. long were diamond drilled, 14 from cores from seven oil-bearing or potentially oil-bearing sandstone, one from a quarried sandstone, and another from a shallow, oil-free, pure orthoquartzite. The plugs were extracted in toluene, vacuum dried, and then weighed. The plugs were then saturated with kerosene and reweighed. This gave pore volume. They were then again extracted in toluene. vacuum dried, and jacketed in 0.004-in. thick copper foil. A short length of -in high pressure tubing led from the jacket. The jacketed sample was then saturated with kerosene and placed in the hydraulic pressure cell as shown in Fig. 1. Movement of the mercury slug showed the change in pore volume of the sample under a given pressure change. Temperature was about 72" F for all measurements. Before taking any measurements, the external pressure on the sample was raised to 10,000 psig and kept there for several hours. This formed the thin copper jacket around the plug. 7-he pressure was then removed and the sample left at least overnight, during which time it recovered its original pore volume. An external pressure of 12,000 psig was then applied and held constant. The change in pore volume was noted as the internal pressure was raised in 1,000 psig increments to 10,000 psig. The internal pressure was then reduced to zero, again noting pore volume changes,

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970