Search Documents

By C. F. Gates, H. J. Ramey

Recent literature shows that pronounced increases in oil recovery can result from the use of miscible systems in recovery operations. This literature also points out certain problems associated with maintaining miscibility, e.g., compositional changes resulting from alterations in pressure and temperature or from zone dilution. The purpose of the work described in this paper was to study miscible zones for possible use in waterflooding operations and examine the manner in which these zones break down. The fluid system selected for study was oil—terNary-butyl alcohol—water, which, because of its narrow range of miscibility, is ideally suited for investigation in short, linear systems. The laboratory-observed effects on oil recovery and zone stability of a number of variables are reported. Among the variables investigated were miscible zone size, viscosity ratio, length of travel, and interstitial water saturation. Results obtained with this system show that interstitial water adversely affects recoveries because of premature phase break. Salinity of the water further aggravates Ais situation. This effect, although pronounced for this system, would probably occur in any fluid system containing an alcohol. As expected, viscosity ratio has a decided influence on recovery efficiency and zone stability. The results also show that use of a suficieno zone size, even though the system has poor phase characteristics, yields higher oil recoveries than are obtained from straight water floods. Length-of-travel studies showed a square root relationship between zone growth and path length for favorable viscosity ratios. No such clear-cut dependence was observed under unfavorable viscosity conditions. INTRODUCTION As new oil reserves have become increasingly more difficult and expensive to find and develop, the oil industry has devoted more and more time and money to finding more efficient methods for exploiting known reserves. The increasing number of water floods is one manifestation of this attempt to improve oil recoveries from existing fields. Since considerable quantities of oil are by-passed by the waterflood process, it is only a partial answer to the problem. The search continues. Of late, several proposed processes have received the critical attention of investigators. Among them are gas drives, in situ thermal reactions and miscible displacements. The present paper is concerned with miscible displacements in conjunction with water floods. Recent literature1- reports that, under proper conditions, oil recoveries from LPG floods, propane sweeps, condensing gas drives, etc., approach 100 per cent of the oil in place in the region contacted by the flood. This literature also points out some of the problems associated with these processes, one of the main ones being the maintaining of miscibility. Compositional changes resulting from pressure and temperature variations or from zone dilution give rise to phase breaks and loss of miscibility, the conclusion being that, after loss of miscibility, much of the beneficial influence of the zone is lost.' The results of the subject work show that, if sufficient zone material exists ahead of a phase interface, considerable improvement in oil recovery can be obtained. In order to increase our understanding of the behavior of miscible systems as they degenerate into immiscible processes, a laboratory system having rapid deterioration characteristics was needed. The fluid system oil—ter tiary-butyl alcohol (TBA)—water possesses these properties and was the system investigated. Although these fluids are different from those generally employed in miscible studies, an understanding of their behavior can be applied to other miscible slug processes. However, under conditions where mass transfer maintains a miscible system through multiple contacts, e.g., enriched gas drive or high pressure gas drive, the observed behavior of alcohol zones would not be applicable. The variables of interest were zone size, viscosity ratio, length of travel, initial water saturation, rate and core characteristics. Of these, only the effects of the following on zone breakdown and oil recovery are reported in detail: (1) zone size, (2) path length, and (3) initial water saturation. The effects of viscosity ratio and core characteristics are discussed in conjunction with other variables. Although initial water saturation is not envisioned as important to the phase behavior of hydrocarbon systems, e.g., LPG, it is important in alcohol systems. The properties of the fluids and cores used are given in Table 1. As previously stated, this fluid system was selected because of the limited range of miscibility

By H. S. Price, J. C. Cavendish, R. S. Varga

This paper presents a new technique for solving some of the partial differential equations that are commonly used in simulating reservoir performance. The results of applying this technique to a simple problem show that one obtains accurate pressure values near wells, as well as accurate pressure pdients, which can be explicitly calculated. The method is completely rigorous in that convergence of the discrete numerical solution to the continuous solution for both pressure and pressure gradient is established. High-order, piecewise-polynomial approximations are used near the wells where pressure gradients are steep, while low-order, piecewise-polynomial approximations are used elsewhere to reduce greatly the calculation time. This combination is shown to give a uniformly good approximation to the solution. These approximations, obtained by using a Galerkin process with suitable Hermite subspaces, are shown to be theoretically and numerically superior to the usual approximations obtained from standard finite-difference techniques. Not only are much greater accuracies obtained, but computer times are also greatly reduced. The application of this technique to multiphase flow problems (e.g., single well coning problems) would have considerable practical interest, but such extensions of this technique with full mathematical rigor have not been made as yet. However, the numerical methods presented here are general, and in principle extend to multidimensional, multiphase flow. Moreover, the preliminary results given in this paper are sufficiently encouraging that we feel the effort in attempting these extensions is justified. INTRODUCTION The problem of obtaining accurate pressure distributions and pressure gradients around wells is of considerable importance in the numerical simulation of reservoir performance. The most common approach to solving this problem is to use finite difference techniques (see McCarty and Barfield9 or Peaceman and Rachfordl0). This approach, however, has many disadvantages, the major one being that many grid points are generally necessary for accurately describing the pressure distribution and the pressure gradient around wells. This need for a fine grid results in large computer times and often in prohibitively high costs. Besides investigating the method of finite differences, some authors, such as Welge and weber17 and Roper, Merchant and Duvall,13 have considered a combination of analytical and numerical techniques with some success. These approaches, however, are all nonrigorous and quite often cannot even be applied. In this paper, we present a numerical formulation of high-order accuracy, based on the Galerkin method, for solving this problem. We treat here only the partial differential equation that describes steady-state, single-phase flow. However, the methods presented are general and in principle extend to multidimensional, multiphase flow. Specifically we treat special cases of the problem in two dimensions described by: where G is a rectangle (with sides parallel to the coordinate axis) with boundary ?G, ?/?n denotes the outward normal, and y and B are non-negative constants such that y + B > 0.

Jan 1, 1970

By D. R. Shreve, L. W. Welch

An analysis for predicting the behavior of reservoirs exploited by a combined gas drive and gravity draintrge mechanism is presented. The method allows the prediction of gas-oil ratio, oil production rate, and cumlulative oil recovery for high relief pools in which pressure is maintained either by gas injection in the crest of the structure or by a combined gas and water injection program. The reservoir is represented by a model of communicating vertical prisms, with production or injection at the upper or lower boundaries of each prism. The theory of macroscopic gravitational segregation of two fluids within a reservoir is discussed in detail for all ranges of saturation. A description of an application of the analysis to a typical field problem is also presented. INTRODUCTION The importance of gravity drainage as a producing mechanism has become more evident during recent years. However, many of the reservoir engineering methods which have been used in the past for predicting the behavior of reservoirs producing by this mechanism have required such drastic, simplifying assumptions that the calculated performance in some cases does not resemble the actual pool behavior. The mechanism of fluid displacement in porous media has been described by Buckley and Leverett,1 and a simplified "one zone" method based on this development has been used by engineers to evaluate the performance of gravity drainage pools. Welge2 has modi- fied this approach so that the calculation time has been markedly decreased. In addition, many other authors3-9 have analyzed various phases of the gravity drainage producing mechanism during the past several years. In a high relief reservoir in which gravity may affect the recovery mechanism, there are two main types of exploitation programs which can be employed: 1. Reservoir pressure is maintained by injecting gas into the upper portions of the reservoir. In this case gas is not released from solution and, if dynamic pressure drops are small compared to the total pressure of the system, the gas and oil may be considered to be immiscible and to have constant fluid characteristics. 2. Reservoir pressure is not maintained and a "natural depletion" mechanism prevails. The term "natural depletion" can be used to describe the production history of a reservoir in which a gradual decline of average pressure and the associated release of solution gas continue throughout the producing life. In this mechanism gas is released throughout the reservoir and there is a continuous change in the characteristics of both the gas and the oil phases. It is the purpose of this paper to present the results of mathematical developments by which a method has been devised for evaluating the performance of reservoirs which are depleted by the first type of exploitation program discussed above. It should be mentioned that, although throughout the analysis gas and oil represent the two fluids existing in the reservoir, the method is equally applicable to gravity drainage in water-oil systems. STATEMENT OF PROBLEM In reservoirs of pronounced structural relief, where the closure of the oil section is great or the angle of

Jan 1, 1957

By R. Simon, D. J. Graue

This paper presents correlations for predicting the solu-lility, swelling and viscosity behavior of CO2-crude oil sys8i.m~. The correlations were developed from experimental data obtained by the authors. These data are also presented. The data were determined by measuring the properties of mixtures of CO, and nine different oils. Experiinental conditions covered a range of 100 to 250°F and pressures up to 2,300 psia. Properties predicted by the correlations have average deviations, expressed as per cent of experimental value, of 2 per cent for solubility, 0.5 per cent for swelling and 12 per cent for viscosity. INTRODUCTION Interest in CO2 injection as an oil recovery process has led to the development of performance prediction methods which can be applied to specific reservoirs.1 iS To use these performance prediction methods, it is necessary to know the solubility, swelling and viscosity properties of CO2-crude oil mixtures at reservoir conditions. Some information on these properties has appeared in the literature; however, this information did not cover the range of different oils and conditions needed to prepare generalized correlations for reservoir engineering purposes. Consequently, an experimental program was undertaken to collect the data needed. The data obtained and the correlations developed from the data are described in the following sections of this paper. SOLUBILITY OF CO2 IN CRUDE OILS CO2 solubility data in the literature come from six principal sources. The solubility prediction method of Welker and Dunlop3 is limited to 80F. The information in Ref. 4 is of two types: the first includes binary and ternary mixtures of CO, and light hydrocarbons (C1 to C6), and the second gives data for CO1 and heavy hydrocarbons for a temperature range of 40 to 90F. Ref. 5 contains a KCO2 chart for systems whose convergence pressure is 4,000 psia. The KCO2's are based mainly on CO2-natural gas mixtures. Poettmann's work covered CO2 solubility in one condensate and one crude oil6,7 Ja-coby and Rzasa measured CO? solubilities as a function of pressure and temperature for two natural gas-absorber oil mixtures and two natural gas-crude oil mixtures.'8 CO, concentration in these four systems was fixed at 5 mol per cent. The work reported in this paper extends CO, solubility data to a variety of different crude oil types in a temperature range from 110 to 250F and pressures up to 2,300 psia. The experimental procedure used by the authors to obtain the solubility data consisted of combining known amounts of pure CO, and crude oil in a visual cell at a fixed temperature and measuring the bubble point of the mixture. Measurements were made for a total of 40 different CO2-oil mixtures and the results are shown in Table 2. The mixtures included nine different oils (seven crude oils and two refined oils) whose properties are listed in Table 1. All nine oils had vapor pressures less than 1 atm at the experimental temperatures. Consequently, analysis of the bubble-point vapor showed a CO, concentration over 99 mol per cent. At no time during these experiments was a second, more dense, liquid phase observed. The solubility correlation which was developed from the data in Table 2 is presented in Figs. 1, 2 and 3. In these figures, solubility is expressed as xco2 the mol fraction of CO, in the CO2-oil mixture. Fig. 1 shows solubility as a function of CO2 fugacity6 and temperature. Fig. 2 shows the same solubility data expressed as a function of saturation pressure and temperature. The solubility shown in Figs. 1 and 2 is for an oil whose UOP characterization factor is 11.7. UOP characterization factors of crude oils can be determined from Ref. 10 if the viscosity and APT gravity of the oil are known. Fig. 3 gives the solubility correction factor for oils whose UOP characterization factors differ from 11.7. The solubility correlation in Figs. 1, 2 and 3 predicted

Jan 1, 1966

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