Search Documents

By Mohamed Mortada

A practical analytical method is developed in this paper which provides the practicing reservoir engineer with a handy method for analyzing oilfield interference problems. The procedure employed entails the principle of superposition of individual eflects. A set of charts has been developed and can be used to evaluate the individual effects. The charts represent non-steady state pressure distribution in an extensive aquifer due to a constant rate of production from an oilfield which can be approximated by circular geometry; they are constructed from solutions to the diffusivity equation considering a circular sink as distinguished from a point sink. The solutions are obtained partially from mathenlatical considerations and partially from analog and digital computations. By virtue of the method of images, the charts can be used also to treat the case of production near a fault. In addition, provisions are made for treating oil fields, which are produced by partial water-drive. Treating an actual field case is deferred to a subyequent writing; however, illustrative example is presented to demonstrate the method of application. The pressure distribution in linear and limited radial aquifers is currently under study. INTRODUCTION With the increasing difficulty of discovering new oil reserves, attention is necessarily focused on the efficient development and production of existing reservoirs. Efficient production depends on the availability of reliable means for predicting the pressure-production behavior of oil fields. The prediction is usually performed by analytical methods or computers. Analytical methods are often limited to systems which can be approximated by simple geometry. However, analytical methods have the distinct advantage of being readily performed without special equipment. Analog and digital computers lend themselves to a variety of problems, but are not always locally available and consequently their use is limited. One particular situation in which analytical techniques are needed involves studies of multiple, water-drive oil fields sharing a common aquifer. Such oil fields are in hydrodynamic communication. Therefore, production from any of the oil fields is accompanied by a pressure decline which is transmitted through the aquifer to the other oil fields and manifested as pressure interference. The rate of propagation of the pressure decline is such that the pressure may be significantly reduced many miles away from a producing pool. An example of pressure interference is found in the Woodbine formation in East Texas. There. a dozen or more pools, including several of major size, draw on the Woodbine formation for water drive. The pressure in the water-bearing portion of the formation has been significantly reduced over an area of several thousand square miles. This reduction of pressure has been such that when the Hawkins pool was discovered in 1940, the pressure in it had already been reduced 280 psi below the initial value in the aquifer. This resulted from prior production from other Woodbine sand pools, including, in particular, the East Texas field some 20 miles to the east of Hawkins.' The above example indicates that the interference effect can be significant and should be taken into con-

Jan 1, 1956

By A. M. Khan, I. Fatt

A new type of steady-state thermal conductivity apparatus is presented and described. Results of measurements on Berea sandstone, Solehofen limestone, Bandera sandstone and Bakelite are presented. The measurements were made in a two-hour period at temperatures up to 80C and axial stress to 125 bars. Results show that the thermoelectric heat pump and tbermistor-controlled electric heater can be used in a steady-state apparatus for measuring thermal conductivity of mck. Steady-state is reached within two hours. With Urethane foam insulation, not more than 3 per cent of heat is lost by radial flaw. The apparatus can be used as a comparator or for absolute thermal conductivity determinations. INTRODUCTION Calculations of heat flow through the earth, both for scientific purposes and engineering design, require thermal-conductivity data on rock. However, the high variability of rock often makes handbook thermal-conductivity data of little value. For engineering purposes in particular, it is desirable to have measured data on samples of the rock mass under study. Laboratory methods for measuring thermal conductivity of rock use either time-dependent1-3 or steady - state temperatures .4,5 The time -dependent (transient) techniques give thermal diffusivity from which thermal conductivity can be calculated if the specific heat and density of the sample is known. Most transient techniques require lengthy and involved sample preparation. Also, the accuracy tends to be low. Steady-state techniques are more exact and are the most commonly used. The steady-state methods, however, also have their difficulties. A long time is often required to reach steady state, troubles have been encountered in establishing good thermal contact between the heating and cooling elements and the rock sample, and unaccounted-for heat losses are difficult to prevent. In this paper we present a description of a new type of steady-state thermal conductivity apparatus. With it, measurements have been made in two hours at temperatures up to 80C and axial stress to 125 bars. Sample preparation is simple and no thermal conductivity standard is needed.6 APPARATUS Our steady-state apparatus is an improvement . over the conduction pile used by Zierfuss and van der vliet4 and Somerton.5 A linear heat flux is established through a conduction pile consisting of the rock sample and a standard material of known thermal conductivity in series. The thermal conductivity of the sample is calculated from the ratio of the temperature drop across the sample

Jan 1, 1966

By R. F. Krueger, L. C. Vogel

A design is presented for an electrical analog computer which can solve non-steady state heat transfer problems in an extensive radial forrnation containing a moving cylindrical source. The computer is used to simulate a .simplified thermal oil recovery process in which heat transfer from a moving, constant temperature source takes place radially by conduction only. Temperature distribution curves are .shown for several different assumed modes of travel for the heat source. The data are used to estimate the residual fuel requirements necessary to maintain a self-propagating isotherma1 front for the particular system being studied. Although the computer is designed to represent an unique system. conversion factors can he adjusted to show the effects of changing the assumed values of system constants, such as thermal conductivity of the formation and temperature of the source isotherm. Examples are given to illustrate the relative influence the assumed values of these variab1es will have upon the quantity of fuel required to maintain the source. The maximum effect of heat transfer by the flowing air stream is estimated mathematically, and full requirements for this system are shown to he markedly reduced. Inasmuch as the examples discussrd refer to a simplified process, the results are not directly applicable to a practical field operation. It is believed, however, that the data il1ustr-ate some general trends which are important in thermal recovery processes. Data from laboratory and field experimentation can he used to modify the computer to take into account the influence of heat transfer by various mechanisms, for example, by injected air and by fluids ahead of the hot zone. With .such refinements it would be possible to estimate more accurately the limiting values for the rate of travel of the high temperature front and for the required air- injection rates. Use of this computer should aid in evaluating the economic feasibility of oil recovery by various thermal processes. It is hoped this paper will stimulate further work by others to help accomplish this objective. As the costs of exploration for and development of new oil fields rise. increasing the recovery of oil from established fields becomes more and more important. In fields containing low gravity crude oil, ultimate recoveries by primary means are sometimes less than 10 per cent and often total only 10 to 20 per cent of the oil in place. Conventional secondary recovery methods are not always effective in increasing this total. Methods which have been proposed for increasing the recovery of viscous oils include thermal processes for increasing the mobility of the oil. One method is to provide a moving heat source generated within the formation to heat the oil and a flowing gaseous medium (for example, air) to sweep the oil to the producing wells. The possibility of such a process has intrigued engineers for many years.'. '." but little or no information has been available in the literature relative to the technical or economic feasibility of this method. Recently, however, interest has been revived by publication of laboratory and field test results indicating the feasibility of maintaining

Jan 1, 1956

By B. C. Craft, Murray F. Hawkins, John R. Rayne

A reservoir mechanism of sulfur recovery by the Frasch process is presented. Improving the economic, of recovery appears to be largely a well, rather than a reservoir problem. A most important factor is the limitation of the lateral extraction of .sulfur about the wells due to the almost vertical .flow of the injected hot water. Model studies are described which confirin the mechanism of sulfur recovery by presently used methods. Studies on a five-spot system indicate improved econo171ics regarding recovery, recovery rate, therinal efficiency, and well spacing. INTRODUCTION Sulfur mining by the Frasch process is one of the most important methods of sulfur production. Developed by Herman Frasch in 1894, it enjoyed a rapid growth and presently supplies 40 per cent of the world's sulfur and about 80 per cent of the U. S. sulfur." The Frasch process of mining the underground native sulfur deposits involves melting the sulfur in place by introducing super-heated water into the sulfur bearing formations, and producing the sulfur as a liquid. Many problems are encountered in mining these native sulfur deposits, and sulfur production depends chiefly upon the nature of the sulfur deposit and the underground thermal efficiency of the heat transfer system. THE RECOVERY MECHANISM Sulfur is typically found in salt dome cap rock as crystal aggregates occurring within the cavities, seams, and pores of limestone and gypsum formations.' Sulfur-bearing limestone is typically cavernous, vugular, and fractured. Porosity averages about 20 per cent and sulfur saturation averages approximately 25 per cent by weight.1,15 Vhe sulfur bearing formations are very heterogeneous in nature and physical characteristics vary widely within small distances. Total cap rock thickness may vary from 50 to more than 1,000 ft, and net productive thicknesses from a few feet to a few hundred feet. In most cases the upper portion of the limestone formation is barren of sulfur, varying from 5 to 200 ft in thickness.1,15,33 The immediately overlying sediments and the underlying dense anhydrite formation provide permeability barriers to the flow of water from the cap rock section.'" Sulfur wells are drilled with conventional oilfield drilling rigs. Fig. 1 illustrates the casing program of a completed sulfur well. The cap rock sections are drilled with core bits and continuous core logs are recorded. After coring, the cap rock section is underreamed to a diameter of 9 7/8 in. The bottom of the 8-in. casing is welded closed and rests on the bottom of the drill hole, which extends through the cap rock to a few feet above the bottom of the sulfur-bearing formation. The 8-in. pipe contains 5 ft of 3/4-in. perforations on 6-in. centers and 15 ft of 1-in. perforations on 1-ft centers. The 20-ft perforated section is set opposite the lower 20 ft of the sulfur-bearing formation. The 4-in. sulfur pipe is set inside the 8-in. water line on a seat ring at the top of the 5-ft section of 3/4-in. perforations. The mining water, heated to about 330°F under a pressure of 250 psi, is forced down the annulus between the 8-in. and 441-1. pipes and through the 1-in. perforations to the sulfur-bearing formation. Compressed air, at 400 psi pressure, is supplied through a 1-in. air line near the bottom of the well to provide air-lift of the liquid

Jan 1, 1957

By M. R. J. Wyllie, Michael A. Torcaso

With the continued deep drilling of today, increasing numbers of high pressure and high temperature gas-condensate reservoirs are being discovered. Correspondingly, the ranges of properties of gas-condensate reservoir fluids are also being increased. The literature1*'." contains several excellent reviews of the property ranges in presently known condensate reservoirs. When one examines these property ranges, it becomes apparent that the need for reliable data is great if the engineer is to evaluate properly or appraise a gas-condensate reservoir. It is the purpose of this paper to present phase-behavior data that will point out some of the important variables that influence recovery from gas-condensate reservoirs produced by pressure depletion. All these variables may be resolved into two major categories: the hydroc:arbon composition of the reservoir fluid and the physical properties of the reservoir. But the variables which we will consider are those which are usually determined from well tests shortly after discovery and are, consequently, readily available to the engineer. They are as follows: (1) the producing separator gas-stock-tank liquid ratio, (2) stock-tank liquid gravity, and (3) reservoir temperature. Laboratory depletion studies have been performed on several reservoir fluids to illustrate the effect of changes in the cronsidered variables upon ultimate recoveries. The deviation factors, Z, of a series of reservoir fluids were determined experimentally to illustrate their order of magnitude and their effect upon the calculation of in-place reserves existing at high pressurcs and temperatures. The procedures followed in conjunction with the performance of these studies and their associated calculations are illustrated in Fig. 1. For the purpose of this paper, the depletion recoveries are expressed only in terms of pcr cent of initial stock-tank liquid, with the abandonment pressure assumed to be 500 psig. VARIATION OF RECOVERY WITH INITIAL PRODUCING SEPARATOR GAS-STOCK-TANK LIQUID RATIO The producing separator gas-stock-tank liquid ratio has long been one of the principal criteria in assumptions and generalized correlations' used for the estimation of ultimate recovery from gas-condensate reservoirs. Fig. 2 presents a plot of the cumulative recovery of stock-tank liquid at abandonment as a function of the initial producing separator gas-stock-tank liquid ratio. The graphical illustration represents the results of 44 laboratory depletion studies performed on actual gas-condensate fluids. These data are presented not only to illustrate the general trend of ultimate recovery increase with initial producing ratio increase but to emphasize the considerable variation in possible recovery at any given producing ratio. Table 1 presents the properties of two actual reservoir fluids chosen to illustrate the approximate extremities of producing ratios for gas-condensate systems. The fluids designated as A and B exhibited initial producing ratios of 3,400 and 114,000 scf/bbl, respectively. Fig. 3 presents a graphical comparison of the cumulative stock-tank liquid recoveries, calculated from the well-stream compositions and produced volumes for each fluid during depletion. It is to be noted that Fluid A exhibited a recovery of only 16.5 per cent when produced to an assumed abandonment pressure of 500 psig while the corresponding recovery for Fluid B was calculated to be 62.2 per cent. This difference in recovery, which exemplifies the general trend of recovery increase with initial producing ratio increase, is readily explained by consid-

By J. S. Aronofsky, J. P. Heller

This paper presents a mathematical analysis of the fluid mixing which occurs during flow through porous media. The analysis is based on the well-known diffusion equation with mass transfer term. It is pointed out that the use of this equation is justified by a single general assumption which does not specify any particular mechanism for the mixing process. Formulas are given for tracer fluid concentrations for two different boundary and initial conditions. Calculated numerical values compare closely with some published results of experiments in which no viscosity or density difference existed between displaced and displacing fluids. INTRODUCTION A phenomenological theory of the mixing and diffusion process is presented for the flow in a porous medium of two miscible liquids of equal viscosity and density. This is based on the classical diffusion equation which has been used to explain such processes as Brownian motion and heat conduction. More recently, this same differential equation has been derived (in a manner analogous to the work in this paper) for the problem of heat transfer during fluid flow in porous media.' The purpose of this paper is to show that the available experimental evidence justifies the use of this simple diffusion equation. The recent data of Koch and Slobod' and of von Rosenberg' are used to test this diffusion model by comparing concentration curve shapes and calculating "effective diffusion coefficients". The same differential equation for miscible displacement, as well as some experimental confirmation, has been published in Japan by K. Yuhara,4 and recently in this country by Day.5 Yuhara arrives at his con-clusions through analogy with turbulent diffusion, de- scribing the microscopic velocities in terms of "eddy motion" — although unlike the ordinary eddies of hydraulics. Day proceeds directly from Scheidegger's6 concept of the dispersion of velocities. After a very clear description of Scheidegger's theory, he applies it to calculate the motion first of a drop and then of a thin layer of brine into a bed of sand. For the latter case, Day derives and solves the same differential equation as Eq. 1. In our presentation no detailed assumptions are made as to the mechanism of displacement. In fact, several microscopic processes may well be involved—the only general assumption made is that their cumulative microscopic effect may be described by the addition of a second-order term to the differentirl equation of continuity. DERIVATION OF EQUATIONS Let us consider the displacement in one spatial dimension of Liquid A from a sample of porous media by Liquid B. A and B are assumed completely miscible and of the same viscosity and density. Both fluids are also assumed to be incompressible and it is specified that the total flow is the same for any cross section. It is evident that the rate at which Liquid B would be carried across such a section would be given by the product of the fluid velocity, u, and the concentration, C, of B in the fluid mixture assuming the constancy of u and C over the cross section. Since the volume of Fluid B is conserved, this leads to the one-way wave equation - u ?C/?x = ? ?C/?l where x increases in the direction of u—the "Darcy velocity", Q/A (or velocity flux per unit area), and 9 is the fractional porosity. Consider the case where fluid B is displacing Fluid A and the concentration distribution is represented by a step funcLon (with a vertical discontinuity). The solution of Eq. 1 predicts no mixing

Jan 1, 1958

By B. S. Gottfried

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By K. P. Fournier

This report describes work on the problem of predicting oil recovery from a reservoir into which water is injected at a temperature higher than the reservoir temperature, taking into account effects of vis cosity-ratio reduction, heat loss and thermal expansion. It includes the derivation of the equations involved, the finite difference equations used to solve the partial differential equation which models the system, and the results obtained using the IBM 1620 and 7090-1401 computers. Figures and tables show present results of this study of recovery as a function of reservoir thickness and injection rate. For a possible reservoir hot water flood in which 1,000 BWPD at 250F are injected, an additional 5 per cent recovery of oil in place in a swept 1,000-ft-radius reservoir is predicted after injection of one pore volume of water. INTRODUCTION The problem of predicting oil recovery from the injection of hot water has been discussed by several researchers. l-6,19 In no case has the problem of predicting heat losses been rigorously incorporated into the recovery and displacement calculation problem. Willman et al. describe an approximate method of such treatment. 1 The calculation of heat losses in a reservoir and the corresponding temperature distribution while injecting a hot fluid has been attempted by several authors.7'8 In this report a method is presented to numerically predict the oil displacement by hot water in a radial system, taking into account the heat losses to adjacent strata, changes in viscosity ratio with temperature and the thermal-expansion effect for both oil and water. DERIVATION OF BASIC EQUATIONS We start with the familiar Buckley-Leverettg equation for a radial system:* This is sometimes referred to as the Lagrangian form of the displacement equation. If the injection water is at some temperature Ti which is above reservoir temperature TR, then lw, the flowing fraction which is water, is a function of both water saturation and fluid temperature at any given time and position in the reservoir. The temperature dependence of lw results from the temperature dependence of the viscosity ratio µo/µm of which lw is a function. Thus:

Jan 1, 1966

By A. Schild

The effect of heating a producing well on the rate of oil recovery has been analyzed in the simple case of a well producing oil by a radial drive and in the steady state. Differences of production rates with and without borehole heating have been calculated and plotted; these measure some of the benefits of the thermal recovery process. A striking feature is that the production increase tends rapidly to a constant value for large production rates. Some concrete examples, including a rough comparison with field data. are discussed. Heating requirements have been calculated and plotted on the assumption that the heaters used are 100 per cent efficient. Heating has been thought of as applying only to low volume wells producing heavy oils. However, the calculations suggest that the thermal process considered might be profitable even for some wells producing large quantities of light oil, especially in unconsolidated formations where fracturing is not possible. In the case of wells producing other fluids with the oil or of wells producing partly or wholly by gravity drainage, the theory developed does not apply in detail. It can be expected, however, that any general qualitative be-havior will remain valid. Under normal temperature conditions, high viscosity emulsions sometimes form around the wellbore, lowering the permeability in the vicinity of the well and impeding the production of oil. Borehole heating often prevents such emulsification. This important effect is not easily amenable to mathematical analysis and has not been included in the theory presented here. It should therefore be kept in mind that the benefits due to borehole heating may in fact exceed the values given by the following calculations. INTRODUCTION The recovery of viscous, low gravity oils can often be stimulated by applying heat to the production zone. The higher temperature lowers the viscosity of the oil, increasing its mobility, and therefore its rate of production. One thermal method of recovery 1,2,3,4. consists in operating a heater inside a producing well and, at the same time, producing oil from the well. In this paper a simple theoretical model is analyzed to study the effect of such heating on the recovery of oil. The reservoir mechanics of the process is the only aspect of the problem which is considered here. The details of the heating system are ignored, and it is simply assumed that the heater raises the temperature in the immediate vicinity of the producing well to a pre-assigned temperature T w. Although heating requirements are calculated, these refer only to the heat supplied to the formation and to the oil being produced; they do not include heat losses due to the inefficiency of the heater. A short outline of the theory follows: The heat supplied to the producing well flows out into the formation by means of thermal conduction. At the same time, the flow of oil carries heat by convection from the formation towards the producing well. After sufficient time has elapsed, a steady-state is reached where the temperature, the velocity of the oil, etc., no longer vary with time but are functions of position only. In the steady state the temperature decays like an inverse power from the high value Tw at the wellbore to the original reservoir temperature T. at some distance away from the well. The effect of the high temperature and low oil viscosity near the producing well is the same as if the oil were produced at reservoir temperature To but through a producing well of a larger radius. From this point of view, the energy and money being expended to provide heat

Jan 1, 1958

By H. J. Morel-Seytoux

Methods of predicting the influence of pattern geometry and mobility ratio on water flooding recovery predictions are discussed. Two methods of calculation are used separately or concurrently. The analytical method yields exact solutions in a convenient form for a unit mobility ratio piston-like displacement. A few typical pressure distributions, sweep efficiencies and oil recoveries are presented for various patterns. For non-unit mobility ratio, one may resort to a numerical method, such as that of Sheldon and Douherty. 1,2 Because the domains of applicability of the analytical and numerical techniques overlap, the exact solutions provide estimates of the errors in the numerical procedures. The advantages of the analytical and numerical methods can be combined. To develop a numerical technique as independent of geometry as possible, the physical space is transformed into a standard rectangle. The entire effect of geometry is rendered through one term, the "scale-factor", derived from mapping relations. The scale factor can be calcu-lated from the exact unit-mobility ratio solution for the particular pattern of interest. By this means recovery performances for arbitrary mobility ratio con be obtained for many patterns. A sample of tesults obtained in this manner is presented. INTRODUCTION Pattern geometry and mobility ratio are two major factors in making a waterflood recovery prediction. Because assisted recovery has become increasingly important to the oil industry, pattern configuration and mobility ratio also assume a greater significance in the assessment of the economic value of recovery projects. The influence of pattern geometry and mobility ratio in shaping a recovery curve and on the other quantities of interest to the reservoir engineer is the main subject of this paper. Much effort has already been spent on estimating quantitatively the influence of either pattern or mobility ratio or both on oil recovery. The literature reports many investigations of this nature. 3-9 However, many results or methods of recovery prediction presented in the literature cannot be considered fully satisfactory. Even for unit mobility ratio and piston-like displacement, where analytical solutions are available, the literature shows discrepancies. For non-unit mobility ratio, the divergence in the results is extreme. For infinite mobility ratio in a repeated five-spot, depending on the investigator, the sweep efficiency ranges from 0 per cent to 60 per cent. With respect to the influence of pattern on recovery, only the repeated five-spot has received much attention. Other confined patterns and pilot configurations have received very little attention. Two calculation methods are presented in this paper, either separately or concurrently: the analytical method of potential theory and the numerical method of finite-difference approximation. The analytical method is more restricted in scope than the finite-difference method, but it has the definite advantage of providing exact solutions within its range of applicability. If a unit-mobility ratio piston-like displacement is assumed, the analytical approach is possible. A few typical results are reported in this paper; the detailed description of the general method and of a great variety of results will be the subject of other articles. Fornon-unity mobility ratio, we must resort to a numerical scheme. The numerical technique is that which was described by Sheldon and Dougherty.l,2 It is not limited to piston-like displacement. However, mainly single interface results will be presented here. Because the respective domains of applicability of the analytical and the numerical method overlap, useful comparisons of exact and numerical solutions can be made for a variety of patterns. The advantages of the analytical and numerical approaches can be combined. The reason for the success of this analytical-numerical approach can be summarized in the following two points: 1. The numerical solution for arbitrary mobility ratio can be programmed most efficiently when the physical space in which the displacement actually takes place is transformed into a standard shape; and 2. This can be done with remarkable simplicity whenever an analytical expression for the pressure

Jan 1, 1966

By C. S. Matthews, H. C. Lefkovits

Drilling progress is often delayed by sticking of the drill string. The development of preventive and remedial methods has been hampered by incomplete understanding of the sticking mechanism. A recent lahorntory investigation hns indicated that one type of sticking may be attributed to the difference in pressure between the borehole and formation. This paper shows, by means of soil mechanics, that the primary cause for differential pressure sticking is cessation of pipe movement, whereas diflerential pressre and stanrtding time determine the severity of the sticking. The analysis stresses the importance of using low-weight muds with low solids content and low water loss to alleviate diflerential pressure sticking and describes why packed hole drilling, long strings of drill collars, and a large deviation from the vertical are conducive to sticking. Finally, preventrve and remedial methods ore evaluated, and a theory is presented on the release of stuck pipe by spotting oil. INTRODUCTION Since drilling with long strings of oversize drill collars has become standard practice in many areas, the incidence and severity of the stuck pipe problem has increased. It has been noticed that in the majority of these cases the sticking could not possibly be attributed to key seating or caving of shales. It appeared that, due to the differential pressure between the mud column and the formation fluid, the collars were pressed into the wall and so became "wall stuck". Points to note about differential pressure sticking are: (1) sticking is restricted to the drill collars, (2) the collars become stuck opposite a permeable formation, (3) the sticking occurs after an interruption of pipe movement, (4) circulation, if interrupted, can be restarted after the sticking is noticed, and (5) no large amounts of cuttings are circulated out after restarting circulation. Helmick and Longleyl investigated pipe sticking by differential pressure in the laboratory and found an empirical relationship between the differential pressure, the sticking time and the required pull-out force. In this paper an explanation of the mechanism is given based on Terzaghi's theory of clay consolidation. A qualitative description is given in the following paragraphs while the derivation of fonnulas is given in Appendices. This paper is a first attempt to explain pressure differential sticking and many points will require additional theoretical and practical investigation before the problem can be fully understood. PRESSURE DIFFERENTIAL STICKING AS A CONSOLIDATION PROBLEM In any borehole, where the mud pressure is higher than that exerted by the formation fluids, a mud cake is formed opposite the permeable sections of the hole and a continuous flow of filtrate takes place from the mud, through the cake and into the formation. This radial flow pattern requires a certain distribution of the hydraulic and the effective (grain-to-grain) stresses inside the mud cake. Any quantitative or qualitative change in the external pressure conditions will produce a change in the flow pattern and, consequently, also in the internal stress distribution inside the cake. In view of the low permeability and the high compressibility of a clay mud cake, the adjustment of the internal stress distribution is slow and is accompanied by a change in volume. Time dependent stresses are thus created which gradually diminish as the new state of equilibrium between internal and external pressures is approached. Some 30 years ago, Terzaghi developed his "Theory of Consolidation" to account for the time-dependent stresses and settling of clay formations under the influence of external loads. He derived a differential equation by which the time-dependent hydraulic stress and the consolidation can be computed for any point inside the layer during the consolidation process. His theory is based on the assumption that the change in stress is solely due to a change in water content and it may only be applied to one-dimensional consolidation phenomena. Other investiga-tors5,10 have expanded his theory to include processes of more than one dimension. The difference between the external pressures on the mud cake before and after sticking is a qualitative one (isolation of part of the cake by the static contact with the drill collars after pipe movement has been stopped)', and the time-dependent stresses thus created may be investigated by means of Terzaghi's theory. By this analysis the changes in the nature of the contact surface between the drill collars and the mud cake during the sticking can be explained; and the friction force between the two may be computed as a function of the sticking time, the borehole dimensions and the mud cake characteristics.

By H. Hooykaas

A few years ago Marsal published a method of predicting oil recovery as a function of time for an edge-water-drive reservoir with several rows of wells.' The method is based solely on oil- and water-production data obtained during a short period for the first row of wells and the gross productions (past and future) of all wells. This paper presents a streamlined version of the original publication and experimental evidence showing the reliability of this version. For this purpose scaled experirnents were performed in a large model (750 X 500 X 23 cm, 24.4 ft X 16.2 ft X 8.9 in.) in which edge-water-drive reservoirs with three rows of production wells were simulated. From the results of the experiments it is concluded that the method gives a fairly accurate prediction of the future production behavior of the reservoir and enables us to select the most profitable production plan. Results obtained with the method were always slightly conservative. INTRODUCTION The calculation method developed by Marsal aims at forecasting future reservoir behavior.' It can be applied to floods in which a number of down-dip wells are rapidly going to water. Essentially it is a method for predicting the future behavior of the up-structure wells on the basis of the observed behavior of the down-structure wells. As such, the method is especially useful for predicting behavior of edge-water drives toward several rows of wells. The requirement to be met for justifying the approach is identity of production mechanism in all regions of the reservoir. The attraction of the method then lies in the fact that it applies even when we have no complete description of the prevailing production mechanism. The original publication was in the German language and for this reason, perhaps, did not attract widespread attention. In our opinion it deserves to do so because it does offer in a simple way results which are in fair agreement with observed data. The aim of the present paper is to bring the method to the attention of a wider audience. The method in its most simplified form is reviewed and experimental data in support of it are presented. PRINCIPLES To show Marsal's principle. we will first write the continuity equation for very simple conditions, linear flow (in the positive x-direction) in an incompressible porous body [cross-section A(x), porosity +(x)], saturated with incompressible fluids (oil saturation S). When the oil in this formation is displaced by water, injected at x = 0 (inlet face), we get where q(t) is the total gross displacement rate and f the fractional oil flow. The equation can be solved if we make the assumptions: (1 ) that the fractional oil flow (f) is a univalent function of saturation (S); and (2) that the function is the same throughout the porous medium. Then: f = f(S) Solving this equation with the method of characteristics, we obtain

Jan 1, 1966

By J. S. Lasater

Resu1ts of experinmental measurernents of heat capacities and thermal conductivities of some typical porous rocks are presented. Measured heat capacities agree closely with va1ues calculated front known chenlical . compositions of the rocks. On the basis of this agreement, heat capacitir.7 of fluid-satttrated rocks were calc~ililted. Therrnal conductivities of the rock samples were measured under varioits conditions of fluid saturation. From these data therrnal dif-fvsivities were calculated. A significant variation of thermal diflusivity with trnlperature is indicated. Knowledge of the thermal characteristics of fluid-bearing porous rocks has become increasingly important with the development of thermal oil-recovery processes. Reliable thermal data for petroleum reservoir rocks are not available and approximate values are generally assumed for reservoir calculations. 'The effects of temperature, pressure and fluid saturation on thermal properties of porous rocks have not been carefully studied. The purpose of the present work was to measure and report the thermal characteristics of some typical porous rocks. A group of eight sedimentary rock samples, including limestone, shale, sandstone and silt-stone, was selected for the study. Properties reported include heat capacities, thermal conductivities and thermal diffusivities. The effect of fluid saturation on these properties was also considered. Determinations of thermal properties require precise measurements which are time consuming and difficult to reduce to a routine basis. Methods have not been developed for measurement of thermal proper-ties under conditions of pressure, temperature, and fluid saturation, which may be encountered in thermal recovery processes. As one approach to these problems, some preliminary work on methods of estimating thermal properties from other more readily measurable characteristics of the rock-fluid system is reported. The samples selected for the studies were oilfield cores from within, or closely associated with oil-productive zones. Descriptions of individual samples are given in Table 1. Clay minerals present in the samples were determined by differential thermal analysis (DTA). Quartz content was estimated from DTA cooling curves. For all but the thermal conductivity tests, a representative 200-gm sample of each core was disaggregated, thoroughly extracted in diox-ane and oven dried at 100°C. The cleaned and dried sample was care- fully split into appropriate portions for the several tests. Chemical analyses of the samples were obtained to permit calculations of heat capacities by methods which are outlined later. Resuits of chenlical analyses are given in Table 2. Heat contents of the test samples were measured by a Runsen-type calorimeter. The heat content of the unknown sample is measured relative to the known heat content of platinum. The apparatus is referred to as a relative error type because of the comparative measurements, and has the advantage of canceling systematic errors. Experimental results are reproducible within 0.5 per cent. A temperature base of 298.16 K is taken for heat content measurements. Measurements were made at five temperatures (260, 440, 620, 800 and 980° F). It was necessary to preheal the sample above the maximum experimental temperature to eliminate volatile constituents. Assuming the weight loss on pre-heat, ing was mostly water, a correction to the observed heat content was made to account for the water originally present in the sample. This latter procedure was followed in the determination of calculated heat contents as will be described later. The apparatus and the method of measure-

By H. H. Rachford, D. S. Howard

In the withdrawal of fluids from tilted aquifers it is of value to be able to predict pressure patterns during the course of the pressure decline. As an example of this, in the displacement of fluids from oil reservoirs being produced by water-drive, prediction of the reservoir pressure is an important factor in fixing production rate. The prediction of pressure patterns in the depletion of an inclined aquifer has conventionally been made by solving the well known flow equations for the pressure Patterns in a horizontal aquifer, and correcting these patterm for the pressure difference due to the hydrostatic head associated with the vertical distance between a point in the horizontal aquifer and the corresponding point in the tilted aquifer. Inasmuch as the equation describing flow in tilted strata is different from that for horizontal flow. there was question as to the validity of the simple corrective procedure current1y in use. To examine the magnitude of error made using this procedure, an inclined aquifer was chosen for which the error was expected to he large, and the equation for flow in this tilted stratum was solved numerically. The results of the conventional procedure of correcting the solution of the equation for horizontal flow was compared with the solution of the correct equation by studying the difference between predicted pressures by the two methods. Three forms of the hydrostatic head correction were studied. The results for all three forms of the head correction are similar. The error in the approximate solution is quite dependent on the position of the horizontal aquifer with respect to the inclined one. The error is greatest at the producing end ad increases with amount of fluid withdrawn. The magnitude of the errors is not large if a good choice is made for the position of the horizontal datum. The most satisfactory position for the datum in all cases is near or slightly about the mid-point of the inclined stratum. Use of the best approximate method studied with the best datum tested yielded predicted pressures in error no greater than 7 psi. INTRODUCTION The movement of water by expansion from an aquifer to displace fluids from a contiguous oil reservoir was first recognized over 20 years ago.9 Since that time much work has been done to relate this mechanism to reservoir behavior mathematically. The result of this effort on the part of many workers has led to an understanding of underground fluid movement which has enabled us to follow and predict the behavior of the water-drive reservoir. One reason that good progress has been made in the application of mathematics to the study of pressure patterns and water movements in aquifers has been that the differential equation which describes this process is of a relatively simple form. The equation for the isothermal, horizontal, viscous flow of water through homogeneous sandstone is identical with that for the diffusion of heat through homogeneous solids, and this equation and means of obtaining solutions for it with many types of boundary conditions have been studied extensively. That this work has been highly successful is clearly emphasized by the large number of solutions developed for problems of interest both by methods of formal mathemat-

Jan 1, 1957

By J. J. Arps

This paper reviews the methods currently in use for estimating primary oil reserves and discusses the principles on which these methods are based. Particular emphasis is placed on how these methods change with the type of information available during the life cycle of an oil property. This paper contains various novel estimating methods and shortcuts heretofore unpublished. INTRODUCTION Estimating oil reserves is one of the most important phases of the work of a petroleum engineer since the solutions to the problems he deals with usually depend on a comparison of the estimated cost in terms of dollars, with the anticipated result in terms of barrels of oil. His recommendations to management regarding the best course of action are therefore normally based on the most favorable balance between these two. Specific engineering problems which require such a knowledge of recoverable oil reserves and a projection of future rates are: [a.] the exploitation and development of an oil reservoir; [b.] the construction of gasoline plants, pipelines and refineries; [c.] the division of ownership in unitized projects; [d.] the price to be paid in case of a sale or purchase of an oil property, and the magnitude of the loan which it will support; [e.] the proper depreciation rate for the investment in oil properties; and [f.] evaluation of the results of an exploration program. This discussion will be confined to the various methods and tools which are currently in use for estimating oil reserves to be obtained during the primary phase of an oil-producing reservoir and for a projection of the future production rates. Reserves which may be obtained by secondary recovery methods or fluid injection programs and gas and gas condensate reserves will not be discussed in this paper. Unfortunately, reliable oil reserve figures are most urgently needed during the early stages when only a minimum amount of information is available. Management's interest in the oil recovery from a property—aside from its use for accounting purposes—usually declines when the property approaches its economic limit, just at the time when the reliability of the estimates is at its best. To illustrate this general idea a chart is presented (Fig. 1) showing the three periods in the life of an imaginary oil property. Time is shown on the horizontal axis, while the cumulative production and estimated ultimate recovery are plotted vertically. No particular units are used and this schematical chart is not to scale. During the first period, before any wells are drilled on a property, any estimates will of necessity be of a very general nature, based on experience from similar pools or wells in the same area, and usually expressed in barrels per acre. This will

Jan 1, 1957

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 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 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 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 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