Search Documents

By J. H. Henderson, J. Naar

The displacement of a wetting fluid from a porous medium by a non-wetting fluid (drainage) is now reasonably well understood. A complete explanation has yet to be found for the analogous case of a wetting fluid being spontaneously imbibed and the non-wetting phase displaced (imbibition). During the displacement of oil or gas by water in a water-wet sand, the porous medium ordinarily imbibes water. The amount of oil recovered, the cost of recovery and the production history seem then to be controlled mainly by pore geometry. The influence of pore geometry is reflected in drainage and imbibition capillary-pressure curves and relative permeability curves. Relative permeability curves for a particular consolidated sand show that at any given saturation the permeability to oil during imbibition is smaller than during drainage. ' Low imbibition permeabilities suggest that the non-wetting phase, oil or gas, is gradually trapped by the advancing water. This paper describes a mathematical image (model) of consolidated porous rock based on the concept of the trapping of the non-wetting phase during the imbibition process. The following items have been derived from the model. 1. A direct relation between the relative permeability characteristics during imbibition and those observed during drainage. 2. A theoretical limit for the fractional amount of oil or gas recoverable by imbibition. 3. An expression for the resistivity index which can be used in connection with the formula for wetting-phase relative permeability to check the consistency of the model. 4. The Iimits of flow performance for a given rock dictated by complete wetting by either oil or water. 5. The factors controlling oil recovery by imbibition in the presence of free gas. The complexity of a porous medium is such that drastic simplifications must be introduced to obtain a model amenable to mathematical treatment. Many parameters have been introduced by others in Iprogressing" from the parallel-capillary model to the randomly interconnected capillary models independently proposed by Wyllie and Gardner3 and Marshall4. To these a further complication must be added since an imbibition model must trap part of the non-wetting phase during imbibition of the wetting phase. Like so many of the previously introduced complications, this fluid-block was introduced to make the model performance fit the observed imbibition flow behavior. THE IMBIBITION MODEL ASSUMPTIONS The following assumptions are common to many two-phase models previously described. 1. The wetting characteristics of the medium are clearly defined. 2. When two fluids co-exist in the medium, the wetting phase will fill the smaller pores and the non-wetting fluid the bigger ones. 3. The pores have a circular cross section. 4. Apparent pore radii are related to the capillary pressure (drainage) by

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By H. N. Hall

Various factors must be considered in an engineering evaluation of gravity-drainage reservoirs. Among these are: (1) the effect of producing rate on total oil recovery; (2) the effect upon well productivity and ultimate recovery of the pressure level maintained during the producing life of the reservoir; (3) the economic advantage of full or partial pressure maintenance; and (4) estimate of the rate of gas production and injection and the possible purchase of gas under conditions of full pressure maintenance to ascertain compressor facilities needed. All of these factors can be evaluated only when a reliable method is employed for determining reservoir performance in gravity-drainage reservoirs. The purpose of this paper is to present a general method for calculating the performance of a gravity-drainage reservoir. This method is applicable for conditions of complete pressure maintenance, partial pressure maintenance and normal pressure depletion. Provisions are made to take into account variations throughout the reservoir of reservoir configuration, changes in permeability and fluid composition. Based on the method presented in this paper, an IBM 650 computer program has been developed. The past performance of an actual gravity-drainage reservoir producing under conditions of declining pressure and no gas injection was duplicated using this program. INTRODUCTION In tilted reservoirs the production of oil is influenced by drainage of oil from upstructure to downstructure locations. When this downstructure drainage of oil is sufficient to cause effective segregation of the gas and oil in a reservoir, the reservoir is usually classified as a segregation drive or gravity-drainage reservoir. (Discussion will be restricted to gravity-drainage reservoirs which have no encroachment of edge water.) The important feature in gravity-drainage reservoirs is the density difference between reservoir oil and gas. These phases tend to segregate in the reservoir with the result that in the gas cap the oil saturation is maintained at a higher level by drainage of oil from the gas-cap area. Oil can be produced from the oil zone at a low gas-oil ratio and reservoir energy is thereby conserved. The standard material balance in not adequate for predicting gravity-dramage reservoir performance because it does not take into account the difference in saturation above and below the gas-oil contact. Several authors'.' have presented methods for calculating the performance of gravity-drainage reservoirs in which reservoir pressure is maintained constant by gas injection into the gas cap. Using some simplifying assumptions, these methods can be employed with a desk calculator to give acceptable results. The problem of predicting the performance of gravity-drainage reservoirs under the conditions of declining reservoir pressure is many time more complex than that of constant pressure. fierefore, attempis to develop a method suitable for desk calculation have required excessively simplified assumptions. In the past several years, highspeed digital computers have become more widely available for reservoir engineering problems. These corn puters are well suited to problems such as the prediction of the performance of gravity-drainage reservoirs with pressure decline. Many of the simplifying assumptions necessary for hand computation can be eliminated so that a realistic approach to the gravity-drainage process can be made. CONCEPTUAL PICTURE OF OIL MOVEMENT IN GRAVITY-DRAINAGE RESERVOIRS Before attempting to develop an analytical treatment for conditions occuring in a gravity-drainage reservoir, a concept should be formed concerning the movement of fluids in the reservoir as oil is produced. A review of the literature'.' shows that it is customary to classify gravity-drainage operations into two categories—(1) with complete pressure maintenance, and (2) with declining pressure. The same line of reasoning will be followed in presenting the concept of the movement of fluids in the reservoir because it is easier to visualize the movement of fluids under conditions of complete pressure maintenance. After discussing complete pressure maintenance, an analogy will be made between that and the case of declining pressure. It should be kept in mind throughout that the final aim for the problem of solving gravity-drainage performance with digital computers will be to develop a general program for any kind of gravity-drainage reservoir. COMPLETE PRESSURE MAINTENANCE One feature which is generally common in gravity-drainage reservoirs is a gas cap located at the top of the structure. This is shown in Fig. 1(a). Fig. l(b) shows oil saturations that might occur through the reservoir. In the gas cap, oil saturation is lower than

By R. E. Cook

This work presents the development and application of equations of the form developed by Martin1 to describe gravity segregation performance during natural depletion. One-dimensional depletion analyses are applied to several hypothetical reservoirs. Considerable attention is given to the changes in saturation distribution in the dip direction resulting from depletion at various rates. Vertical gas saturation distribution within a sand having vertical permeability is studied to illustrate conditions which provide For low gus-oil ratio production from wells selectively completed away from sand tops, as observed from well performance. Results of performance calculations are shown graphically and illustrate the rate sensitivity of the gravity segregation mechanism and the influence of reservoir geometry and withdrawal distribution. It is shown that for sufficiently high depletion rates gas saturation distribution in the dip direction is characterized by the formation of two discontinuities, or saturation "fronts"; the updip front being the gas-oil contact of the secondary gas cap and the downdip front occurring at a position dependent upon the reservoir depletion rate, withdrawal distribution and reservoir geometry. The gas saturation in the region between fronts is shown to continually increase with time, while that downdip from the lower front remains very low and relatively stable. The effect of increased depletion rate is to reduce the region of stable gas saturation and thus enlarge the portion of the reservoir producing with continually increasing gas saturation. Special attention is given to the performance of reservoirs having sufficient vertical permeability to permit liberated gas to segregate vertically within the sand section before proceeding in the updip direction. It is shown that reservoirs of this type are far less sensitive to depletion rate due to the mechanical control of producing gas-oil ratios afforded by selective well completions away from sand tops. It is shown that gravity segregation performance is adversely af- fected by decreasing reservoir cross-section area in the updip direction. It is also shown that down-dip withdrawal concentrations can restrict gas segregation rates in the dip direction. INTRODUCTION The natural depletion performance of several large producing reservoirs in Western Venezuela has unmistakably shown the significant role that gravity forces play in governing the gas-oil ratio behavior of individual wells and the over - all pressure-production performance of the reservoirs. Analysis and prediction of the natural depletion performance of these reservoirs is of great importance in determining the economics of pressure maintenance operations. In 1957, a method for the analysis of pressure maintenance performance was presented by Martin1 which, in the original and more general form of the equations involved, provides a physically sound basis for the analysis of one-dimensional gravity segregation performance during natural depletion. The work reported herein is an application and extension of Martin's outstanding developments to a series of hypothetical reservoirs in an attempt to explain the general behavior of the gravity segregation mechanism and the influence of such parameters as depletion rate, reservoir geometry and withdrawal distribution. The study considers two types of flow: "distributed flow", where vertical permeability in the reservoir is zero and fluids flow only in the dip direction while uniformly distributed over the sand thickness; and "segregated flow", where vertical permeability is sufficient to permit gas to segregate against the sand top before proceeding updip. This latter type of flow has been recognized as an important mechanism active in the massive sand reservoirs of Western Venezuela. ANALYSIS OF ONE-DIMENSIONAL DISTRIBUTED FLOW SYSTEMS For an inclined reservoir having zero vertical permeability, the velocities of oil and gas flow at any point in the reservoir can be derived from Darcy's equations as

By H. Dykstra

Calculated wellbore pressures were obtained for parameters of radilcs ratio and permeability. In all cases bur two, after-production was allowed to occur for one day. The calculated pressure build-up data were compared with actual pressure build-up data from a condensate well. The paper discusses the comparison, gives reasons for, inability to match calculated (2nd actual data, and presents conc1usions derived from the comparison. INTRODUCTION Considerable information has been published on the decline in wellbore pressure for oil and gas fields resulting from production of fluids. Considerable information also has been published on pressure build-up for oil fields, but relatively little has been published for gas and condensate fields. Perrine' summarized the equations derived for pressure build-up in oil wells and showed how the different methods could be applied. On the other hand, very little theoretical information has been published on pressure build-up in gas or gas-condensate wells.2-4 Tracy,' by pointing out the similarity between the equation describing oil wellbore pressure decline and the equation describing gas wellbore pressure decline, showed how methods of pressure build-up analysis could be applied to a gas well. In analyzing pressure build-up data for gas or gas-condensate wells, it would be very desirable to compare actual pressure build-up data with calculated pressure build-up data. Such a comparison would lead to a better feeling for the quantitative picture of gas-well build-up curves and would help to establish a degree of confidence in a given build-up curve. This paper discusses the comparison and presents conclusions made from it. Calculated wellbore pressures are given for parameters of radius ratio and permeability. In all cases but two, after-production was allowed to occur for one day. The work was done in an attempt to evaluate the reserves of a gas-condensate field. As will be discussed later, it was not possible to make an evaluation. METHOD OF CALCULATION The method used for solving the depletion and pressure build-up behavior for a gas-condensate well was that developed by West, Garvin and Sheldon.5 or the IBM 704 computer program used, it was assumed that liquid condensing out of solution would have only a minor effect on the flow behavior. Therefore, the field was treated as having single-phase flow, with the gas saturation S, remaining constant at unity minus the connate-water saturation Sw. As will be discussed later, this may not have been a good assumption. It also was assumed that oil production could be considered as equivalent gas production by using a conversion factor based on an analysis of the liquid produced. An outer boundary condition of no-flow was assumed. BASIC DATA The data required for the study were reservoir properties, fluid properties and production data. Reservoir properties included the following: thickness h, 179 ft; porosity, 0.194; connate-water saturation, 0.43; gas saturation, 0.57; permeabilities k, 1.0, 0.5, 0.25 and 0.15 md; radius ratios re/rw, 800, 1,000 and 2,000; wellbore radius rw, 5 in.; initial reservoir pressure P, 6,529 psia; and reservoir temperature, 272°F. The selection of permeabilities was based on core-analysis data showing an average air permeability of about 2 md. With a relatively high connate-water saturation, the average effective gas permeability would be about one-half, or less, of the air permeability. The selection of radius ratios was based partly on the observed rapid decline in wellbore pressure during production and partly on the fact that the well was believed to be located in a relatively small fault block. Reservoir fluid properties were obtained from a reservoir fluid study on a recombined sample of gas and condensate. Gas formation-volume factors were calculated from pressure-volume measurements made at reservoir temperature. Gas viscosities were calculated by the method of Carr, et al,G from an analysis of the re-combined sample. These data are shown as a function of pressure in Table 1. The production schedule from the gas-condensate

By D. D. Smith, D. P. Squier, E. L. Dougherty

A system of differential equations describing the temperature behavior of fluid injected at constant surface temperature in a well is derived and .solved analytically. A formula for the fluid temperature at any time and depth is given, as well us a special formula valid for very large times. These formulas are used to calculate temperatures for several typical cases. The results indicate that, initially, the temperature of the water entering the formation is considerably lower than the injection temperature. This condition lasts for only a short period— less than three days for most cases of practical interest. Following this highly transient period, during which the temperature of the fluid entering the formation builds up to about 50 to 75 per cent of the injection temperature. the system enters a quasi-steady state in which the temperature changes are very slow. After severl years, the bottom-hole temperature will still be 50" to 100°F lower than the injection temperature, hilt the heat losses may he tolerable. INTRODUCTION Predicting the behavior of a hot-water flood requires that the temperature of the water entering the injection interval be estimated. This report describes the development and solution of a system of equations which describes the temperature behavior of the injected water in the wellbore with certain simplifying assumptions. The only previous means known to the authors for describing such a process is that of Moss and White.' Their results appear to be close to those obtained by our method in the practical cases which were compared; this agreement is largely due to the fact that in our method temperature soon approaches a quasi-steady state, as was assumed in their method throughout. However, our model covers all times, is continuous (whereas the Moss-White model depends on breaking the depth into discrete intervals) and. we feel. more closely describes the physical problem. FORMULATION OF THE PROBLEM PHYSICAL SYSTEM AND ASSUMPTIONS The injection procedure consists of pumping water at a fixed surface temperature T., down an infinitely long cylindrical well or tubing of inner radius Any material exterior to the water column such as mud, casing, or cement is regarded as part of the formation. The general behavior of the system may be described qualitatively as follows. When the hot water is first introduced into the system, the temperature difference between the formation and the water is large, resulting in a high rate of heat transfer. As a result, the temperature adjacent to the wellbore rises very quickly. Because the segment of the formation adjacent to the wellbore largely controls the heat transfer rate, the heat transfer rate will become relatively constant when this portion has reached a temperature close to that of the water opposite it. The temperature of the water and formation then increase very slowly with time. The length of the initial highly transient period and the temperature of the water at its conclusion will be functions of depth, injection rate, injection-string radius, surface injection temperature and the physical properties associated with the water-formation system. The following additional assumptions were made. 1. There is no heat transfer by radiation in the system. 2. There is no heat transfer by conduction in the vertical direction in either the injection stream or the formation. 3. The linear volumetric and mass flow rate of the water is constant throughout the injection stream. 4. No horizontal temperature gradient exists in the injection stream. 5. The product of density and heat capacity is constant for both the water and the formation, and the formation thermal conductivity is constant. 6. Initially, both the water in the wellbore and the reservoir are at a temperature given by the (constant) ambient surface temperature plus the product of depth and geothermal gradient (assumed constant). At large distances for the wellbore (r m), the formation will remain at this temperature. 7. The water temperature and the formation temperature at r — r,, are equal for all depths D. DERIVATION OF EQUATIONS The differential equation satisfied by the fluid temperature T,(D, t), which is obtained by writing a heat balance on a cylindrical differential of volume dV of the injection string between the depths D and D i dD, is

By C. F. Weinaug

In recent years, several methods have been proposed for correlation of physical properties and behavior of hydrocarbon systems. These correlations were found to be unsatis. factory for predicting differential gas liberation and shrinkage data now obtained experimentally. To achieve the desired accuracy, methods were developed to adjust the correlations with two experimental pieces of data. The required data are the saturation pressure at reservoir temperature and the specific volume or density at these conditions. The saturation pressure is used to adjust the vapor-liquid equilibrium-ratio correlations. The density measurement is used to determine the apparent density of liquid methane. With these adjustments, differential gas liberation and shrinkage are predicted within the experimental error. EQUILIBRIUM RATIOS Equilibrium ratios published by the National Gas Assn. of America' were used as a basis for the calculation. The convergence pressure for the sample was determined by using a hypothetical binary system, with methane as one component and the remainder of the system lumped together to form the other component. For the critical temperature of the heavy constituent, the weight average of the critical temperatures of the ethane and heavier components was used. The correlation of pseudo-critical temperature, density and molecular weight of Matthews, Roland and Katzl was used in this calculation for the hexane and heavier fractions. Points for the critical properties of methane and the hypothetical heavy constituent were placed on a chart of the critical loci of hydrocarbon systems. A critical loci for the hypothetical binary system was drawn, and the critical pressure of this system was determined at the reservoir temperature. This critical pressure became the convergence pressure for which equilibrium ratios were read from the NGAA4 charts. The NGAA equilibrium constants for C, through C, were used for the C,; through C, fractions since they were separated on the bases of normal paraffin boiling points. The correlation of Hoffman, Crump and Hocott5 as used to estimate the equilibrium constant for the C fraction. The NGAA eauilibrium-ratio charts are based primarily on paraffin systems and have to be corrected for the presence of naphthenes and aromatics. The correlations in the literature proved too unsatisfactory for determining these corrections. EXperimental data, the saturation pressure, is used with the correlation of Soloman2 to make this adjustment. The correction factors of Soloman were cross-plotted so that the factors for ethane and heavier components were made a function of the factor for methane. a2 = 0.7927 at + 0.1646 .........(1) where a, = correction coefficient for for methane NGAA equilibrium ratio, and a2t = correction coefficient for ethane and heavier NGAA equilibrium ratios. This equation and the equation for the bubble point are solved simultaneously for the correction coefficients. The bubble-point equation is anKnXn = 1.0 . . . . (2) The equilibrium constant for the saturation pressure at reservoir temperature is used. This procedure not only corrects for the relative volatility between components, but also assures that the results will be directly related to the experimental saturation pressure. METHANE DENSITY After determining a set of equilibrium ratios, the apparent density of condensed methane at 60°F and 14.7 psia was estimated. This was found to be necessary because the correlations of Standing3 ere found to be unsatisfactory. The law of additive volumes was employed to make this calculation.

By R. V. Higgins, A. J. Leighton

Interest by petroleum engineers in the flow of three phases—oil, gas and water—in irregularly bounded porous media lies mostly in the performance calculation of water floods of reservoirs that have been partially depleted as the result of expansion of much of the originally dissolved gas. The authors present a method to forecast three-phase flow in complex geometry and explain the details by the use of a specific example of a five-spot water flood of a partially depleted stratified reservoir. In this example, the fluid and rock properties of a field given by Prats, et al,' were used. The computed results fit the field performance more accurately than Prats, et al, and Slider.' The time required for the high-speed digital computer to make the calculations, including the contributions from the different layered zones, is about one minute. INTRODUCTION The declining rate of discoveries of new oil fields in the United States makes the recovery of more oil from known reservoirs more attractive than previously. Water flooding has an excellent proved background for recovering additional oil economically. Accordingly, the main interest of this paper is in this type of recovery. In petroleum engineering studies, commercial interest in three-phase flow is mostly in the water flooding of reservoirs in which the oil has been partially produced by the expansion of dissolved gas. When water is pumped into these reservoirs, three-phase flow takes place. Although this paper is concerned with these conditions, the principles involved could be used for other conditions when and where they come to the fore. Many investigators have made contributions to non-empirical forecast methods using basic scientific engineering principles. Several of these used the oil in place at the start of the water flood and the oil remaining after a large quantity of water has passed through a core for key values in their calculations. Recently, Prats, et al,' and Slider have reduced assumptions by adding the third phase—gas—in their forecasting methods. Slider uses the mobilities in the immediate vicinity of inlet and outlet wells as an aid to simulate the resistance effect in the five-spat pattern of the flow of oil, gas and water. Prats, et al. minimize assumptions by using previously determined laboratory data for determining sweep efficiency in the five-spot pattern. In neither of these papers is the saturation profile continuously affected by permeability-saturation curves. Sheldon and Dougherty- recently described a method that employs continuously changing saturation profiles using permeability curves, has a minimum of assumptions and needs no prior sweep efficiency. The Higgins-Leighton method, described in this paper, has all of these features; however, many of the techniques are different from those of Sheldon and Dougherty. The Higgins-Leighton method, tested in May, 1961, is direct and easy to apply and requires very little computer time to calculate a forecast. The short computer time is especially helpful in the study of a reservoir containing many layers of different relative permeabilities. In the Higgins-Leighton method, the individual pressures do not have to be calculated, as the resistance to flow in each cell in the flow pattern is readily determined without the use of any iterative techniques. The saturation and permeability distributions are readily determined. These data and the shape factor, which is measured only once from the potentiometric model when mobility ratio is one, determine the resistance to flow in each cell. THEORY The authors showed in a previous paper4 that, as an aid to calculating performance, the reservoir can be divided into channels using the streamlines of a potentiometric model as a guide. See Fig. 1. This procedure also was used in this paper. The authors also showed that, by treating conduits as approximately one-dimensional and neglecting pressure gradients transverse to the main flow, the Buckley-Leverett equation may be expressed as The principles expressed by this equation are employed extensively in the three-phase flow, as they were in two-phase flow. In the three-phase flow, the channels (the size and shape of which are taken from a potentiometric model)

By R. W. Nelson

Methods for determining permeability presented by W. D. Kruger3 in a paper entitled "Determining Areal Permeability Distribution by Calculation" have considerable merit. They can be expected to contribute significantly to analysis of flow in heterogeneous systems. Some rather subtle yet important conditions must be met, however, to assure that the computed permeabilities are correct. The conditions which insure a unique permeability distribution may have been overlooked in the method as proposed and presented in the subject paper. The subtlety of the special requirements may permit good agreement between computed and measured pressures in the reservoir; but the computed permeability may still be in error. The general theory and discussion of the requirements for a unique determination of the permeability distribution are presented elsewhere in more detail.7'8 Accordingly, only major features will be presented here in the notation of Kruger. GENERAL DESCRIPTION OF REQUIREMENTS A brief mention of the route to be followed in showing the special requirements may be helpful. The equation to be solved is a quasilinear, first-order, partial differential equation in the unknown permeability. This equation can be solved through an extension of Lagrange's method by reduction to a system of subsidiary ordinary differential equations. Through consideration of the identity of one of the Lagrange subsidiary differential equations with the stream function, a special interrelationship can be shown. The interrelationship has special significance with respect to the boundary condition. If the boundary function satisfies part of the subsidiary differential equation, no unique solution for the permeability distribution exists. In the physical problem this requirement for uniqueness indicates that the permeability measurements used as a boundary condition can not be along a stream tube. The boundary condition must be along a line, or series of lines, which completely traverse the region and pass through every stream tube.

By E. W. Hough, G. L. Stegemeier

Empirical equations for surface tension of propane and normal butane as functions of reduced temperature are obtained from experimental data. Another correlation relating surface tension to enthalpy of vaporization is given for these two compounds. In addition, new paracbor numbers are calculated for the normal paraffin hydrocarbons. These numbers are utilized for the calculation of interfacial tension of two-component systems as functions of pressure and temperature, using a modified form of Weinaug-Katz equation. The experimental data for two binary systems are approximated by the correlation. From these results it is found that the interfacial tension in the high-pressure region remains extremely low at large pressure decrements below the critical pressure. Thus, it appears that condensate systems may have flow characteristics almost like single-phase conditions even though the pressure is within the two-phase region. Experimental data have shown that interfacial tension divided by density difference approaches zero as the critical pressure is approached. A calculation of wetting-phase saturations indicates that the saturation gradient at the two-phase contact becomes increasingly abrupt as the critical pressure is approached. DISCUSSION Prediction of the surface and interfacial tension of the light hydrocarbons and of two-component hydrocarbon mixtures at various temperatures and pressures may be made if other physical properties are known. Extensive experimental work on single-component and binary systemsl is the basis for the correlations outlined in this paper. Interfacial tension is defined as the specific surface-free energy between two phases of unlike fractional composition, while surface tension is defined as the specific surface-free energy between two phase:; of the same fractional composition. The usual definitions relating interfacial tension to a liquid-liquid interface and surface tension to a gas-liquid interface are not clearly defined when the critical region is included, and there is no sharp distinction between a gas and a liquid phase. Interfacial tension is probably the most important single force that makes one-half to one-third of the total oil actually in place in a reservoir rock unrecoverable by conventional gas-drive or waterflood methods. A rough estimate of this figure for the United States is 100 billion bbl. Interfacial tension presently is used by petroleum engineers in the estimation of saturation gradients at the gas-oil contact and at the oil-water contact. The data in this paper should prove useful for estimates of reserves involving gas-oil contacts. Relative permeability undoubtedly is influenced by interfacial tension, for sufficiently small values. These data should be useful in determining how small the values are In addition, these data should eventually add to our fundamental knowledge of surfaces. At the critical point, all surface excesses approach zero and the thickness becomes very large. SINGLE-COMPONENT SYSTEMS It has been observed1 that the following relationships are good approximations to the physical properties of propane and n-butane.

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By J. W. Lacey, F. H. Brinkman, J. E. Faris

This paper presents an analysis of the high-pressure, gas-driven LPG-slug process, based on fluid flow tests in areal models. Two types of tests were made. One series was made in low-pressure models which permitted observation of fluid movement. Three completely misci-ble analog fluids were used. A second series of tests was made in high-pressure models using methane, propane and a light refined oil saturated with methane at room temperature and 1,550 psig. Under the test conditions of room temperature and a pressure level of 1,550 psig, the phase diagram for the fluids used is similar to those for many of the field systems where the process is considered for use. A method for using these laboratory data to calculate field performance of the process is outlined. As a result of this work, it is concluded that small banks of LPG (5 per cent HV or less) are not effective in increasing oil recovery in horizontal reservoirs. Znstead, where small banks are used, the driving gas quickly penetrates the LPG bank because of fingering and channeling; and from this point on, the process behaves essentially as an immiscible gas-injection project. The validity of this conclusion was substantiated by: (1) laboratory studies of the effect of rate, model size and mobility ratio on miscible displacement in areal models; and (2) calculation of field recovery, which compared closely with actual field recovery. INTRODUCTION Field applications and pilot tests of the gas-driven, LPG-bank, oil-recovery process are on the increase. Most of these tests are employing small banks (2% to 5 per cent hydrocarbon volume) of LPG, which is miscible with both the driving gas and the oil in place, in an effort to attain an effective yet economical miscible displacement of oil by gas. The expectation of miscible displacement with small banks is based on the concepts that: (1) lengths of solvent-oil mixed zones, measured during miscible displacements in long slim cores, are representative of those that will occur in the field; and (2) areal sweep efficiencies5 measured in electrolytic model studies are applicable to miscible displacement in reservoirs. Our experimental evidence indicates that the mixed-zone lengths and sweep efficiencies mentioned are not applicable to miscible displacement in reservoirs. In this paper we present an evaluation of the gas-driven, LPG-bank oil-recovery process based on fluid flow experiments in areal models. These results are used to predict the performance of a field pilot test of this process, and the results are compared with the actual test results. THE EFFECTIVENESS OF LPG BANKS The effectiveness of LPG banks of various sizes in accomplishing miscible displacement of oil by gas was determined by displacement tests in a model representing one-quarter of a confined five-spot pattern. The model, 12 X 12 X 1/4 in. in size, was packed uniformly with glass beads. It was operated at a pressure of 1,550 psig and at room temperature; the fluids used were methane, propane and refined oil saturated with methane at 1,550 psig. The saturated oil had a viscosity of 1.2 cp at room temperature. The results of these tests are shown in Fig. 1, where recovery is plotted as a function of the volume of fluid injected for: (1) an immiscible-gas drive with an oil-to-gas viscosity ratio of 85; (2) three sizes of LPG banks; and (3) two completely miscible displacements with mobility ratios of 85:1 and 10:1. (The miscible displacement data plotted are the averages of several tests.) The results show that a 2½ per cent bank of LPG does not increase oil recovery over that obtained by immiscible-gas drive. The 7 and 17 per cent banks are, respectively, about 30 and 50 per cent effective (with reference to the M = 85 miscible-displacement curve) in increasing oil recovery at the point where 2½ hydrocarbon volumes (HV) of fluid have been injected. The percentage effectiveness of the banks at a given volume of fluid injected is defined as Recovery by Bank — Recovery by Gas Drive Recovery by Miscible Recovery by Flooding Gas Drive We conclude from these results that banks of LPG smaller than about 5 per cent HV in an individual stratum will cause from little to no increase in oil recovery. This finding was substantiated by work in low-pressure models, which permitted visual observation of fluid movements. In these tests, three completely miscible fluids having the same viscosity ratios as those used in

By P. Hazebroek, M. Prats, E. E. Allen

Well locations in water floods frequently do not coincide with the ideal well locations associated with the flood pattern which is used. A question of importance is. "How much oil is lost because of this?". This is the problem treated in this paper. It is assumed that each pattern in a five-spot flood has the proper number of production and injection wells but that the wells may be off-pattern. The distance by which a well is off-patern is given in terms of the pattern spacing, defined as the length of one side of a five-spot. The influence of off-pattern wells on the performance of a five-spot flood has been investigated mathematically for unit mobility ratio. All production and injection rates are assumed to remain equal. The loss in recovery resulting from pattern irregularitie: has been calculated at a cumulative water injection corresponding to a 98 per cent water-cut in a regular five-spot. Displacing a single production or injection well by about one-third the pattern distance reduces the total field recovery by 4.4/N per cent of the recovery obtainable with regular spacing, where N is the number of five-spots in the flood (assumed large). When off-pattern wells are separated by at least one regular five-spot and are displaced by less than one-third the pattern distance, the total loss in recovery in a large-scale flood is, at most, 1 per cent of the recovery obtainable from that field with a regular five-spot pattern. However, much higher losses may occur if the off-pattern wells are not separated by a regular five-spot, if the rates are imbalanced, or if the mobility of the injected water is higher than that of the crude. Although the loss in oil recovery due to pattern irregularity often may be low, the off-pattern case will always require injection of more water to reach the same final water-cut. The volume of water injected per unit volume of oil displaced has been calculated for each five-spot adjacent to a five-spot containing a single 08-pattern well. The cumulative water injection at the time of 98 per cent water-cut may be increased by as much as 37 per cent when a production well is displaced one-third the pattern distance toward an injection well. The existence of differences between actual subsurface locations and the idealized well patterns used in secondary recovery operations is well recognized. However, the methods commonly used for predicting waterflood behavior ignore such differences, even when subsurface locations are known. The reason for this is simple. The idealized well patterns are highly symmetrical and, as such, the production history of the pattern (and even of the field) can be determined from that of a few production wells—usually from only one. Consideration of the many irregular well spacings found in a specific field would entail a large amount of work, and then the results would be applicable only to that specific field. In an attempt to provide an answer to the more general problem, this paper considers a water flood having an idealized well pattern, except for one well which is not at its regular location. When several off-pattern wells are present, their contributions are assumed to be additive unless the off-pattern wells are close to each other. Only four positions of the off-center well (at x, y. of Fig. 1) are considered in this paper. These are shown

By J. R. Kyte, V. O. Naumann, C. C. Mattax

Imbibition and waterflooding tests at reservoir conditions show that conventional tests at room conditions sometimes lead to recovery estimates that are too low. Core samples which are of intermediate wettability at normal test conditions may be strongly water-wet at reservoir conditions. This was found to be true for samples from two of three reservoirs investigated. For the other reservoir, less difference was observed between the wettability of samples at normal test conditions and at reservoir conditions. Samples which became more water-wet at reservoir temperature and pressure also water flooded more efficiently at reservoir conditions. These results show that core analysis tests on intermediate wettability samples should be performed at reservoir conditions to obtain representative data. INTRODUCTION Throughout the oil industry, core analysis tests normally are performed using refined oils and at room temperature and pressure. This procedure has always been subject to question because the possibility is generally recognized that the reservoir environment may influence displacement behavior. Burkhardt, et al,' have performed waterflood tests at reservoir conditions on East Texas cores. They found little difference between results obtained at reservoir conditions and results of conventional flooding tests at room conditions. As they pointed out, however, this situation might not be general since they obtained data on only one rock and crude-oil system. This paper presents data which demonstrate that it is sometimes necessary to perform core analysis tests in a simulated reservoir environment to properly represent the reservoir wettability and displacement characteristics. The work is specifically concerned with reservoirs which were thought to be of intermediate wettability on the basis of conventional test results at room conditions, i.e., core samples imbibed only small amounts of either oil or water and exhibited no decided preferential wettability characteristics. Tests were performed on cores from three such reservoirs to determine their behavior in the reservoir environment. lmbibition and waterflood tests were used to demonstrate the effect of reservoir environment on reservoir rock wettability and displacement characteristics. Results of tests performed at reservoir conditions were compared with results from conventional tests. Imbibition tests were performed on samples from all three reservoirs. Waterflood tests were performed on samples from only one reservoir. Fluid properties for the three reservoirs studied are shown in Table 1. Tests at reservoir conditions were performed using bottom-hole crude oil at reservoir temperature and pressure. The aqueous phase used in these tests was a synthetic brine approximating the reservoir brine composition. Conventional tests were performed at room conditions with refined oils and synthetic reservoir brine. The oil-to-water viscosity ratios were adjusted to the reservoir values. The cores used in all the tests were taken from the reservoirs using either brine alone or a bentonite-car-boxymethylcellulose mud without additives as coring fluids. Previous tests show that these fluids have no significant effect on the wettability of reservoir rock samples.1 The cores were preserved in brine at the wellsite to prevent wettability changes during shipment and storage. Wellsite imbibition tests have been performed on fresh samples from several reservoirs (including Reservoir A) and compared with later laboratory imbibition tests on samples preserved in brine. The comparison showed that the wettability of the preserved cores was the same as that of fresh, surfaced cores at the wellsite.2 IMBIBITION TESTS Conventional imbibition tests were performed on preserved cores using the apparatus and procedures described by Bobek, et al.2 A sketch of the apparatus for imbibition tests at reservoir conditions is shown in Fig. 1. The preserved samples were first centrifuged to expel excess water and then placed in a stainless-steel container which served

By W. R. Strickler, P. Hazebroek, M. Prats

The pressure and production behavior of a homogeneous cylindrical reservoir producing a single fluid through a centrally located vertical fracture of limited lateral extent was determined by using mathematical methods to solve the appropriate differential equation. It is assumed that there is no pressure drop within the fracture — that is, that the fracture capacity is infinite. It was found that the production-rate decline of such a reservoir is constant (except for very early times) when the flowing bottom-hole pressure remains constant. The production-rate decline increases as the fracture length increases. Thus, the lateral extent of fractures can be determined from the production-rate declines before and after fracturing or from the decline rate after fracturing when the properties of the formation and fluids are known. The production behavior over most of the productive life of such a fractured reservoir can be represented by an equivalent radial-flow reservoir of equal volume. The effective well radius of this equivalent reservoir is equal to one-fourth the total fracture length (within 7 per cent); the outer radius of this equivalent reservoir is very nearly equal (within 3.5 per cent) to that of the drainage radius of the fractured well. The effective well radius of a reservoir producing at semisteady state is also very nearly equal to one-fourth the total fracture length. It thus appears that the behavior of vertically fractured reservoirs can be interpreted in terms of simple radial-flow reservoirs of large wellbore. INTRODUCTION An earlier report1 has considered the effect of a vertical fracture on a reservoir producing an incompressible fluid. That investigation of the fractured reservoir producing an incompressible fluid was started because of its simplicity. Thus, pertinent behavior of fractured reservoirs was obtained at an early date, while experience was being gained of value in the solution of more complicated fracture problems. One of these more complicated problems, and the one discussed in this report, considers the effect of a compressible fluid (instead of incompressible fluids) on the production behavior of a fractured reservoir. In the incompressible-fluid work mentioned, it was shown that the production rate after fracturing could be described exactly by an effective well radius equal to one-fourth the fracture length whenever the pressure drop in the fracture was negligible. Because of the simplification in interpretation, it is a matter of much interest to determine whether the production behavior of reservoirs producing a compressible liquid could be described in terms of an effective well radius which remains essentially constant over the producing life of the field. The details of the mathematical investigation are given in the Appendixes. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable circular cylindrical outer boundary of radius re. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. Fig. 1 indicates the general three-dimensional geometry of the fractured reservoir just described. When gravity effects are neglected, the flow behavior in the reservoir is independent of the vertical position in the oil sand. Thus, the flow behavior in the fractured reservoir is described by the two-dimensional flow behavior in a horizontal cross-section of the reservoir, such as the one shown in Fig. 2. Production is due to the expan-

By M. Prats

The effect of a sand-filled vertical fracture of limited radial extent and finite capacity (fracture capacity is the product of the permeability and width of the fracture) on the flow behavior of a cylindrical reservoir producing an incompressible fluid through a centrally located well has been investigated mathematically. The shape of the lines of equal pressure near the fracture is essentially independent of the size of the reservoir, provided that the field radius is of the order of the fracture length or larger. For reasonable values of production rates and of fluid, reservoir and fracture properties, the total pressure drop between the end of the fracture and the well is gelrerally negligible compared with the pressure drop in the reservoir. With regard to production response, the effect of vertical fractures can be represented by the production response of an equivalent or effective well radius. For a high-capacity fracture, the effective well radius is a quarter of the total fracture length, decreasing with the fracture capacity. When invasion effects are simulated by decreasing the width of the damaged zone with distance from the well, the effect of formation damage around a fracture on the production response is not so seriolls as indicated by the literature. This suggests that frac fluids with a conventional filter-loss response are better than high-spurt-loss frac fluids, provided the effective permeability of the damaged zone is the same. INTRODUCTION This paper considers the effect of the fracture capacity, as well as the formation damage which can result from fracture treatments, on the productive capacity of vertically fractured wells. Other publications, notably those of van Poollen,l,2 consider these same effects. In addition to providing more general results for vertical fractures than are available from the liter- ature, the present paper gives the equivalent well radius of fractures having different lengths and capacities and, also, includes pressure distributions in and around the fractures. The effect of a damaged zone around a fracture on the production response was not found to be so great as that reported by van Poollen. 2 This difference probably stems from the fact that we consider a damaged zone which is widest (but is still small) near the well and thins out toward the extremities of the fracture, whereas van Poollen considers a damaged zone having a uniform width for the entire fracture length. Simple, but adequate, equations which describe the effect of these variables on production response are presented (in Appendixes A and B). Thus, results can easily be extended to values of the variables not specifically considered here. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The liquid is incompressible. It is assumed that, after the well is fractured, the cylindrical outer "boundary" is at a uniform potential, provided that it is not too near the fracture. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below the producing layer (reservoir).3 Fig. 1 indicates the general three-dimensional geometry of the fractured reservoir described above. If gravity effects are neglected, the flaw behavior in the reservoir will be independent of the vertical position ill the oil sand. This means that the flow behavior in the fractured reservoir can be described by the two-dimensional flow behavior shown in the horizontal cross section of the reservoir, Fig. 2. PARAMETERS The effect of a fracture on the pressure distribution in the fracture itself and in an otherwise undamaged (no skin) reservoir indicates that the pressure distribution is a function of three parameters.

By M. D. Arnold, H. J. Gonzalez, P. B. Crawford

A method is presented for estimating the effective directional permeability ratio and the direction of maximum and minimum permeabilities in anisotropic oil reservoirs. The method is based on the principle that production from a well in an anisotropic reservoir results in elliptical isopo-tentials about the well, rather than circular. Bottom-hole pressure data from three observation wells surrounding a producing well are required to apply the method. The method involves fitting field pressure data to a set of general charts of isopotentials and making a few simple calculations until a solution is found. The method is based on a steady-state equation for homogeneorrs fluid pow. In addition to the method, a brief discussion of the theory underlying it is presented. INTRODUCTION The existence of a different permeability in one direction than another in oil reservoirs has been mentioned in several papers. Hutchinson' reported laboratory tests on 10 limestone cores and pointed out that one-half of them showed significant, preferential, directional permeability ratios, the average being about 16:1. Johnson and Hughesz reported a permeability trend in the Bradford field in the northeast-southwest direction with flow being 25 to 30 per cent greater in that direction. Barfield, Jordan and Moore -eported an effective permeability ratio of 144:1 in the Spraberry. Crawford and Landrum4 showed that sweep efficiencies could often vary by a factor of two to four, and sometimes considerably more, due to variations in flooding direction and patterns in anisotropic media. These findings indicate that the poss'bility of anisotropy may be worthy of consideration in the development of an oil field. In considering this, it should first be determined if anisotropy exists. If it does, the direction of the maximum and minimum permeabilities and the ratio of their magnitudes are quantities which can be of value in planning the most efficient well-spacing patterns. Past methods of determining these quantities have included analysis of oriented cores and analysis of flooding performance of pilot injection patterns. In recent work, Elkins and Skov5 resented an analysis of the pressure behavior in the Spraberry which accounted for anisotropic permeability. This work was based on the transient pres- sure distribution in a porous and permeable medium, with the solution expressed as an exponential integral function involving rock and fluid properties. The purpose of this study is to provide a method, based on steady-state equations, of estimating the direction and relative magnitude of permeabilities in an oil reservoir from field pressure data and well locations only. The method presented is based on work by Muskat6 which shows that Laplace's equation represents the steady-state pressure distribution for homogeneous fluid flow in homogeneous, anisotropic media if the co-ordinates of the system are shrunk or expanded by replacing x with it is desirable that data be obtained early in the history of a field because knowledge of an anisotropic condition would allow new wells to be spaced in such a manner that reservoir development and subsequent secondary recovery programs could be planned more efficiently. THEORETICAL CONSIDERATIONS A brief discussion of the theoretical basis on which the graphical solution was developed is presented in this section. Muskat's two-dimensional6 olution for the pressure distribution in an homogeneous, anisotropic medium with an homogeneous fluid flowing can be algebraically manipulated to show that the isobaric lines are perfect ellipses. The ratio of the major axis to the minor axis, a/b, is related to the permeability ratio, k,/k,, as follows. alb = dk,/k,--...........(1) It can also be shown that the pressure varies linearly with the logarithm of the radial distance from the producing well. However, the gradient along any ray is a function of the orientation of that ray, and a ..xiable is present when anisotropy exists which cancels out for a radial (isotropic) system. For a system such as that described, a dimensionless pressure-drop ratio was developed which is completely independent of the actual magnitude of the pressures. This was done by arranging Muskat's solution in such a way that aIl variables cancelled out except k,/k, and well positions. However, this solution depends on having a co-ordinate system with axes coinciding with the major and minor axes of the elliptical isobars. Thus, it was necessary to introduce a co-ordinate system rotation factor. The two unknown variables are then k,/k. and 0, and the two measured dimensionless pressure-drop ratios are related to the unknown variables as follows.

By R. A. Greenkorn

This project originated in a practical problem—we needed five tracers that could be used together to locate flow paths in a pilot flood. While tracers for subsurface liquids have been used since the turn of the century,1-15 none of those reported in the literature appeared to be either consistent or quantitative enough for our purposes. Most were used in field systems without controlled experiments to determine the accuracy and precision of analysis, and many were tracers that could not be used collectively. The ideal tracer, of course, would follow the fiuid of interest exactly, traveling at the same velocity as the fluid front. But the ideal is impractical to attain because adsorption-desorption effects cause the tracer to lag behind the front; these effects, plus diffusion-dispersion effects, cause the tracer front to spread more than the fluid front. Thus, our objective was not to locate a tracer that would be ideal for all circumstances but rather, to find one that would approximately follow the fluid, or one that under controlled conditions could be corrected to calculate the movement of the fluid front. We considered tracers satisfactory—(1) if they were easy to analyze; (2) if their breakthrough-elution curves were not too different from those for the chloride ion, a tracer believed to follow the fluid front closely; and (3) if we could calculate from the curves a material balance, at 1.25-PV (pore volume) injected, within 5 per cent of that calculated from chloride curves. Of a possible 35 materials, we selected 13 tracers that could be quickly and easily identified and whose analysis was claimed to be accurate within 5 per cent. Only one of these, tritiated water, was a radioactive tracer. Radioactive tracers are easy to detect even in small quantities, but they require special handling and special equipment. Also, those that can be used together are limited because special equipment is required to separate emissions from the various tracers. Two of the original 13 tracers were eliminated in static tests to determine how accurately and precisely they could be analyzed, and to check on gross adsorption. The remaining 11 were flowed through a 9-ft linear sandstone model, and breakthrough-elution curves were obtained. Finally, three tracers were field-tested as breakthrough tracers. These tests are described in the following sections. The 13 tracers considered in these experiments were EDTA (ethylene diamine tetra acetic acid), fluo-rescein, picric acid, salicylic acid and ammonium, boron (as borate), bromide, dichromate, iodide, nitrate and thiocyanate ions, plus chloride ion and tritiated water. All but the chloride ion and tritiated water were subjected to static sand tests to eliminate the tracers that could not be analyzed quickly and accurately and to eliminate those that showed excessive adsorption. All but two of the tracers, EDTA and salicylic acid, qualified for flow tests on this basis. The procedures used in the static tests (see Appendix A for analytical details) were as follows. Each tracer (in the amount shown in Table 1) was dissolved in 1,000 ml of water. Then, 800 ml of this solution was added to 500 gm of sand, and the remaining 200 ml was reserved as a control standard. The sand-tracer mixture was agitated once a day over a 10-week period, and during this period we analyzed five duplicate samples from this mixture, along with five duplicate samples from the control solution. From the control-sample analyses, we constructed a control chart on which the sand-tracer analysis results were plotted. The control chart, and its use to determine accuracy and precision of analytical results, as well as the amount of adsorption, is described in Appendix B. The results of the static tests are summarized in Table 1, which lists the concentration c (concentration at time t) divided by c, (initial concentration). If there was no adsorption, c/c, = 1. The acceptable deviation from this value varied from tracer to tracer, depending upon the limits of analytical precision established for the different tracers. The best results were obtained with boron, bromide and dichromate, with all values falling within limits of accuracy and precision. Ammonium, iodide, nitrate and picric acid also were satisfactory, although the ammonium results were erratic. Flourescein and thiocyanate were unsatisfactory ini-

By H. S. Price, J. E. Warren

Techniques for studying the performance characteristics of heterogeneous reservoirs have been developed. The effect of permeability variation on both the steady-state and the transient flow of a single fluid has been investigated. Limited comparisons with field and laboratory data have been made. The physical models studied consist of random three-dimensional arrays of homogeneous porous blocks. The permeabilities of the individual elements are assigned according to a specific distribution function; uniform anisotropy is introduced by varying the relative dimensions of the blocks. A particular model is perturbed simply by re-arranging its elements at random. The behavior of a physical model is determined by digitally solving its numerical analogue. Based on computational experiments, subject to the restrictions implied by the assumptions that were made, the following general conclusions have been drawn. 1. The most probable behavior of a heterogeneous system approaches that of a homogeneous system with a permeability equal to the geometric mean of the individual permeabilities. 2. The effects of flow geometry and anisotropy on the most probable value of the effective permeability of a heterogeneous medium are finite but not significant. 3. The permeability determined from a pressure build-up curve lor a heterogeneous reservoir gives a reasonable value for the effective permeability of the drainage area. 4. A qualitative measure of the degree of heterogeneity and its spatial configuration are obtained from a comparative study of core analysis and pressure build-up data. It has been indicated that these conclusions are predicated on the assumption that the core analysis and the pressure build-up data represent true reservoir characteristics. Common sources of error in these data have been discussed. One of the most significant problems of reservoir engineering is that of determining the nature and the disposition of the heterogeneities that inevitably occur in petroliferous formations. It is with this problem that this paper is concerned. Extensive research, both theoretical and experimental, has provided reasonable descriptions of the physical processes that are associated with primary and secondaryrecovery mechanisms. Similar progress in the fie1d of numerical analysis and the evolution of truly high-speed digital computers have made it feasible to calculate the behavior of these mechanisms under almost any environmental conditions. It appears to be logical, then, to conclude that the mechanistic performance of any reservoir can be predicted with some confidence. Unfortunate1y, a necessary condition for the validity of this conclusion is that the formation be adequately described in a macroscopic sense. This condition cannot be satisfied ill general; therefore, a frustrating situation currently exists — it is possible to compute the effect of reservoir heterogeneity on behavior, but it is not possible to specify the heterogeneity itself. In many cases, the predicted behavior of a project is so completely dominated by irregularities in the formation's physical properties that the gratuitous assumptior, of a particular form for the variation can reduce the solution of the problem to a mere tautological exercise. While all porous media are microscopically heterogeneous, only macroscopic fluctuations need be considered since the fundamental concepts of reservoir mechanics are based on macroscopic quantities. For the simplest case in which only one physical property varies from point to point, the degree of heterogeneity that is apparent depends on three distinct factors, i.e., the nature of the property's variation (frequency distribution), its disposition (spatial distribution) and the inherent stability of the mechanism being studied. Consequently, any attempt to devise an unambiguous scale of heterogeneity must take these factors into account. Since this investigation is intended to be exploratory, it will be restricted to single-phase flow because of the simplification implied; both steady-state and unsteady-state behavior will be considered. Furthermore, to keep the number of parameters

By J. E. Warren

The mathematical problem considered by the author1 can be given another physical interpretation which is of some practical significance. The alternative physical problem involves the approximate behavior of a single well which is producing at a constant rate through an axially symmetric, horizontal fracture of infinite flow capacity which is located at the center, top or bottom of a uniformly thick, horizontal, homogeneous, anisotropic reservoir of infinite lateral extent which contains a single, slightly compressible fluid. Using Coats' asymptotic results, the wellbore pressure in the fractured system is given by the following equation.