Search Documents

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

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 C. Scala, F. Bernstein

General expressions are derived relating the streaming potential to the electrochemical potential of a permeable junction separating two electrolytic solutions. By the methods of irreversible thermodynamics, these potentials are related to the flux of mobile species arising from differences in pressure and chemical composition. The flux is conveniently expressed as a transference number, the equivalents of a mobile species transferred through the junction per faraday of electricity. At constant pressure, the flux equation gives the electrochemical potential across the permeable junction. The streaming potential across the junction is obtained by imposing the condition of constant chemical composition. The familiar expression for the borehole SP is shown to be a particular solution of the general flux equation. The derived exprersions were verified in experiments with seven shale samples. Both the theory and experimental data provide additional proof of the presence and magniturle of the streaming potential in shales. INTRODUCTION In the interpretation of electric logs, there has been much concern over the correction to be made to the recorded SP for an electrokinetic or streaming potential originating from the difference in pressure between the mud column and the formation. It was long recognized by M. R. J. Wyllie that such an electrokinetic potential source mas in the mud cake opposite permeable sands.' More recently, M. Gondouin and C. Scala2 showed an additional but opposing electrokinetic source of potential to be in shale beds. Their work was supported by the very extensive and excellent work of Hill and Anderson.3 n the present contribution, the general relationships between the electrochemical potential and the corresponding streaming or electrokinetic potential are examined and clarified. A notable advance towards understanding the relation between the various sources of electrical potential in a borehole has been the application of the now common methods of irreversible thermodynamics to describe the passage of ions and water through charged electrolyte junctions such as in shales.' The advantages of these methods are twofold. Unlike classical thermodynamics, they can be applied to non-equilibrium sys- tems, i.e., systems in which flow takes place, provided the flow processes have reached a steady state. They also permit conventional thermodynamic properties to be related to the steady-state dynamic properties. These methods arc used here specifically to derive expressions relating the potential behavior of shales in terms of ion and water transport. Experimental verification of the expression for the shale streaming potential is given which in turn offers theoretical support to the existence of a streaming potential in shales.' Importantly, the relative contribution of transport parameters to the SP is clarified. The theoretical treatment is essentially that of Lori-mer, Boterenbrood and Hermans3 ut has been extended to derive the general expression for the SP.

By J. M. McDowell

A study has been made to deter~r~irre how the behavior of a water drivee reservoir changes as a function of the permeability of the formation and as a function of the size of the aquifer. The effect of pressure maintenance programs on the rote of nat11rn1 water influx is also studied, (IS well as the influence on the influx of the position of water injection we1ls with reference to the oil-water contact. The results are determined by the use of an electrical reservoir analyzer and are plotted as families of curves of pressure, natural water influx, and water injection vs time or cumulalive production. Two size aquifers are studided. One is limited with a 36-inile radius, and the, other approximates an infinite reservoir. The field remains the same in all respects for the completed. study except that the permeability is changed to agree with that assumed for the aquifer. 1 N TR ODUCTI ON this study was undertaken in -dcr to determine how the pres- sure decline and water influx of a water drive field would change as a function of the rate of fluid withdrawals, as a function of the permeability of the formation, and as a function of the size of the aquifer. The effect of pressure maintenance programs on the rate of natural water influx was also studied as was the position of water injection wells with reference to the oil-water contact. The information has been determined by the use of an electrical reservoir analyzer1, sometimes referred to as a Carter type analyzer. This instrument is an analog type computer which can be set up to represent a water-drive oil field. It can reproduce the past performance of the field in terms of analogous electrical quantities when the proper circuit analog of resistors and condensers has been determined, also, it can determine the future behavior for any method of production which may be selected, or permit the comparison among several different methods. Although most applications of the reservoir analyzer"" have been to thc study of specific field performance, it is intended to use the analyzer here to help understand how various reservoir characteristics effect the final performance of the field and to sec under what conditions water injection will be of benefit and to what extent. Within the oil industry the performance of water-drive reservoirs was first described analytically by LV. Hurst' and M. Muskat2 for certain types of boundary conditions. The solutions of the equations developed involved products of exponential functions and Bessel functions. Later publications by Hurst4 and van Ever-clingen" and the hooks of Muskat6 gave further information on reservoirs with other boundary conditions. They also presented rather extensive curves and tables to help evaluate the analytical expressions involved. A rccent series of publications by Chatas. bascd on the work of Hurst and van Everdingen, gives quite complete charts and tables of data to evaluate the pressure decline or water influx in infinite or limited radial reservoirs. and in linear reservoirs. The curves in all of these publications are given in dimensionless form

Jan 1, 1956

By K. Leutwyler, H. L. Bigelow

Jan 1, 1966

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 R. G. Bentsen, D. A. T. Donohue

The cyclic steam injection process has become the most widely applied and most successful thermal recovery technique in use today.1-4 Normally, steam stimulation is repeated several times during the life of a project. The problem is to select the optimum number and size of steam treatments so that the profits realized over the entire life of the project are maximum. Obviously, operating costs rise directly with the frequency and amount of steam injection. Operating efficiency, on the other hand, can be adversely affected if the amount and frequency of stimulation are insufficient. Thus arises a multistage decision problem; that is, a decision must be made at each sub-period (day, for example), as to whether steam should be injected, and, if so, how much. It is desirable that these decisions be made in such a way as to secure maximum net profit before taxes. This can be achieved by dynamic programming. To use dynamic programming, we must be able to predict not only the production rate-time curve resulting from a given stimulation, but also how changes in stimulation policy will affect this curve. In this study, the effect of steam injection on well performance is estimated on the basis of a simple model that takes into account both steady-state and transient flow. Oil production rate is calculated for both types of flow, and the larger of the two values is used. This allows the calculation to change to radial steady-state flow at such time as this mechanism becomes more effective than the transient mechanism. We propose to show how dynamic programming can be used to optimize the steam soak process with respect to net profit. For this, a mathematical representation of the steam soak technique is necessary. A further objective, therefore, is the formulation of a mathematical model capable of realistically simulating the physical process. Oil Production Response Model The advantage of using a mathematical model is its generality and ease of manipulation. However, the model must be an accurate representation of the physical process, for a solution derived from a model can be no better than the model itself. The question now arises as to how well the model should be expected to simulate reality. In general, the simplest model should be tried first, for additional accuracy is apt to require additional cost and time. Also, the search for a high degree of accuracy may yield only greater complexity and difficulty without yielding significantly better results. With this in mind, we kept our model as simple as possible, yet generally representative of the actual physical process. The number of state variables is limited to three, and stochastic variations in the parameters are ignored. As a consequence, use of this model should be restricted to cases in which the underlying assumptions are not severely violated. In

Jan 1, 1970

By Kenneth F. Whinery, John M. Campbell

Basically, the Buckley-Leverett theory involves two systems which are similar in nature but are differentiated by time. These systems may be described by the fractional flow and frontal advance equations which essentially characterize the mechanics of oil movement while being expelled from the reservoir. The fractional flow equation originally developed by Leverettl may be expressed in the more usable form The development of this equation is based on Darcy's law describing fluid flow through porous media and applies to the flow at only one point (it is a point function). For simplification, the fractional flow equation is written in the above form because the capillary forces always increase the fractional flow of the displacing phase regardless of the direction of flow or the displacing phase. If, for simplicity, the effects due to gravity and capillary pressure differences are neglected, the fraction of the displacing fluid, f, at any point in the flowing stream is related to the properties of the system by 1 Thus, it is seen that in the absence of capillary and gravitational effects, fd for a given sand and fluids varies only with saturation and pressure. The magnitude of the viscosity ratio, —a has an effective range (range of about 30) in the system where gas is displacing oil and a much smaller range for water displacing oil. In order to make the fractional flow equation more versatile, it is necessary to connect the fractional flow at a given point and saturation with time. This problem was approached by Buckley and Leverett" who developed the frontal advance equation, Eq. 3 states that the rate of advance of a plane that has a fixed saturation, Sd, is proportional to the change in composition of the flow stream caused by a small change in the saturation of the displacing fluid. It is, essentially, a transformation of a material balance equation representing the net rate of accumulation of the displacing phase within a homogeneous sand block. This accumulation is proportional to the difference between the rate at which the displacing fluid enters the sand and that at which it leaves. Eq. 3 describes the velocity with which a plane of constant displacing phase saturation advances through a porous system. Buckley and Leverett," Babson,' Kern," Welge, and others have adequately discussed the basic mechanism and application of the fractional flow and frontal advance equations. APPLICATIONS OF THE FRONTAL ADVANCE THEORY TO PETROLEUM RECOVERY Two general applications of the Buckley-Leverett frontal advance theory involve the system in which the oil is being displaced by an expanding gas cap overlying the oil zone and that in which the oil is being displaced by water. Displacement by Gas in the Presence of an Immobile Water Saturation A system in which gas is the displacing phase may be thought of as having two forces effecting the displacement process. These forces are the gravitational force and that force exerted by the displacing gas. The gravitational effects control the displacing efficiency of the gas. The gravitational effect will be less at higher rates of flow, thereby reducing the effectiveness of the displacement of the oil by the gas. The more efficient displacements occur at flow rates which are less than the gravity free fall rate. Capillary forces can be neglected without materially changing the magnitude of the gas saturation. The Mile Six pool is used herein to illustrate the calculating procedures in evaluating gas drive-gravity drainage field perfomlance. These calculations represent the determination of the gas-oil contact when the distribution of the hydrocarbon pore volume is considered. Two methods

By R. Ehrlich, F. E. Crane

A consolidated porous medium is mathematically modeled by networks of irregularly shaped, interconnected pore channels. Mechanisms are described that form residual saturations during immiscible displacement both by entire pore channels being bypassed and by fluids being isolated by the movement of an interface within individual pore channels. This latter mechanism is shown to depend strongly on pore channel irregularity. Together, these mechanisms provide an explanation for the drainage-imbibition-hysteresis effect. The calculation of steady-state relative permeabilities, based on a pore-size distribution obtained from a Berea sandstone, is described. These relative permeability curves agree qualitatively with curves that are generally accepted to be typical for highly consolidated materials. In situations where interfacial effects predominate over viscous and gravitational effects, the following conclusions are reached. 1. Relative permeability at a given saturation is everywhere independent of flow rate. 2. Relative permeability is independent of viscosity ratio everywhere except at very low values of wetting-phase relative permeability. 3. Irreducible wetting-phase saturation following steady-state drainage decreases with increasing ratio of nonwetting- to wetting-phase viscosity. 4. lrreducible wetting-phase saturation following unsteady-state drainage is lower than fir steady-state drainage. 5. Irreducible nonwetting-phase saturation 101lowing imbibition is independent of viscosity ratio, whether or not the imbibition is carried out under steady- or unsteady-state conditions. Experiments qualithtiveIy verify the conclusions regarding unsteady-state residual wetting-phase saturation. Implications of these conclusions are discussed. INTRODUCTION Natural and artificial porous materials are generally composed of matrix substance brought together in a more or less random manner. This leads to the creation of a network of interconnected pore spaces of highly irregular shape. Since the geometry of such a network is impossible to describe, we can never obtain a complete description of its flow behavior. We can, however, abstract those properties of the porous medium pertinent to the type of flow under consideration, and thus obtain an adequate description of that flow. Thus, the Kozeny-Carmen equation, by considering a porous medium as a bundle of noninterconnecting capillary tubes, provides an adequate description of single-phase flow. With the addition of a saturation-dependent tortuosity parameter in two-phase flow to account for flow path elongation, the Kozeny-Carmen equation has been used to predict relative permeabilities for the displacement of a wetting liquid by a gas. 1-5 It has long been recognized that relative permeability depends not only on saturation but also on saturation history as we11.6,7 Naar and Henderson8 described a mathematical model in which differences between drainage and imbibition behavior are explained in terms of a bypassing mechanism by which oil is trapped during imbibition. Fatt9 proposed a model for a porous medium that consisted of regular networks of cylindrical tubes of randomly distributed radii. From this he calculated the drainage relative permeability curves. Moore and slobod,10 Rose and Witherspoon,11 and Rose and Cleary12 each considered flow in a pore doublet (a parallel arrangement of a small and large diameter cylindrical capillary tube). They concluded that, because of the different rates of flow in each tube, trapping would occur in one of the tubes; the extent of which would depend upon viscosity ratio and capillary pressure. They limited

Jan 1, 1970

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 T. S. Hutchinson, C. E. Kemp

The bottom water drive analysis presented by Muskat has been extended to include fields with wider well spacings. Curves are presented from which volumetric sweepout, water-oil ratio, and rate relations may be calculated. The effects of non-unity ratios on sweepout eficiency and productivity decline have been studied by an approximate method. Also, the effect of gravity on volumetric sweepout has been approximated for a particurlar geometry. INTRODUCTION The purpose of this paper is to report work done in applying the analysis of bottom water drive reservoirs given by Muskat1 to a larger class of reservoirs, to extend this analysis to include the effects of gravity and non-unity mobility ratios in an approximate manner, and to investigate the rate of decline of productive capacity. The bottom water drive type reservoir is one where the oil is underlain by a water zone of sufficient thickness that water movement is essentially vertical. This system is defined in more detail in the appendix and a one well element is shown diagrammatically in Fig. 1. Muskat obtained solutions to the equations describing the volumetric displacement efficiency of idealized bottom water drive systems by assuming that gravity effects due to differences in density betwecn oil and water were negligible and that the mobilities of oil and water were the same. If either of these assumptions is not made, the complexity of the system makes an analytical solution well nigh impossible at the present stage of mathematical development. Coles and Stade' developed an approximation to consider these factors. The authors have refined approximations, starting with the same basic assumption, which should give a somewhat quantitative picture of bottom water drive performance when mobilities are not equal and the effects of gravity are not negligible. The results of these approximate calculations can then be used to obtain the decline in well productivity as the field is depleted.

Jan 1, 1957

By E. J. Bonet, P. B. Crawford

It is fairly common practice to reinject water into the aquifer near the oil-water interface in water-drive reservoirs. There have been many studies of aquifer behavior without reinjection,l-5 but apparently no analytical studies have been made to indicate performance when injection wells are present in the aquifer. The diffusivity equation may be used to describe the pressure behavior in an aquifer. When sources are present it can be written as A general solution of Eq. 1 may be obtained as the sum of two solutions. One is a general solution of the homogeneous equation, obtained by dropping the source term in Eq. 1 and still satisfying the required boundary conditions. Another is a general solution of Eq. 1 satisfying homogeneous boundary conditions.1° The first solution, which corresponds to the aquifer without reinjection, has been given quantitative and qualitative treatment in the literature.l-5,8,9 The latter solution will be treated here. Solutions with variable injection rates are obtained by superimposing solutions with constant rates, so we have considered only the latter. Five cases have been treated here: 1. A linear aquifer of finite width and extent. Water is injected into the aquifer at point xo, yo at constant rate. The pressure along the line x = 0 is maintained at zero pressure. (See Fig. 1.) Equations showing the pressure distribution, rates, velocity and cumulative influx have been developed. Certain numerical results are shown. 2. A linear aquifer of finite width and extent. Water is injected into the aquifer at point xo, yo at constant rate. The rate of production along the oil-water water contact is maintained constant. Equations showing the pressure distribution inside the aquifer are presented. 3. A linear aquifer of finite width and extent. Constant pressures are maintained at opposite ends of the aquifer. Water is injected at constant rates at point xo, yo. Equations showing the pressure distribution are shown. 4. A linear aquifer of finite width and infinite extent. The pressure is maintained constant at the oil-water contact. Water is injected at a constant rate at point xo, YO. 5. A linear aquifer of finite width and finite extent. The pressure at the oil-water contact is maintained constant. Water is injected at constant rate per unit length along a line x,. The Linear Aquifer The linear closed aquifer, performing at constant terminal pressure was selected as the first case to be analyzed. The development of the equations is based on uniform injection per unit of thickness, and gravity effects are neglected. In practical application this condition is probably most closely achieved when the

Jan 1, 1970

By J. E. Kirby, L. B. Schnitz, H. E. Stamm III

The production history of a gas-cap-drive reservoir was reproduced by calculations, and predictions were made for operations under primary depletion, pressure maintenance by gas injection, and pressure maintenance by water injection. Although the general equations used in this analysis have been available to the industry for a number of years, the procedure has not had extensive use because of the unusually large time requirements for developing numerical results. The use of digital computing equipment for processing the calculation work generally reduced the time usually required for an extensive analysis of this type. Calculation techniques were developed that were equally applicable for solving the gar-cap-drive problem, whether digital computing equipment or desk calculators are used. These techniques are outlined and the procedures are illustrated by the field problem. INTRODUCTION The analysis of a gas-cap-drive res-ervoir requires the simultaneous solution of the volumetric balance1 and displacement2 equations. Equations for the analysis of combination drive reservoirs were presented by Wooddy and Moscrip in a previous paper." These equations permit the evaluation of possible operating processes in order that the optimum method of pro- ducing the reservoir may be determined. Although these equations represent the basic technique known for the detailed analysis of gas-capdrive reservoirs, the calculation processes are quite complex and time consuming; hence, the procedure has not had extensive use. Recently, these equations were satisfactorily applied with the aid of a digital computer to investigate the completion history of a gas-cap-drive reservoir under the following operating methods: primary depletion, depletion under full pressure maintenance by gas injection, and depletion under full pressure maintenance by water injection. The use of digital computing equipment substantially reduced the time necessary for making the computations and permitted a more concise investigation of the reservoir's depletion history than was possible when the desk calculator was used. This more rigorous investigation also permitted a detailed analysis of the calculation procedures necessary in the solution of the gas-cap-drive problem. This paper outlines the procedure of analyzing a gas-cap-drive reservoir and illustrates the process with an actual field problem. CHARACTERISTICS OF THE RESERVOIR ANALYZED Production in the field under consideration is obtained from a Frio formation of Oligocene age at a depth of approximately 5,740 ft sub-sea. The accumulation, which is bounded on the flanks by strand lines, is semirectangular in shape with gentle narrowing toward the base of the structure as shown by the sand limits and completion map in Fig. 1. The reservoir development indicates an oil column of 319 ft that is defined by gas-oil and oil-water contacts. Basic data relative to the characteristics of the formation and the formation fluids are shown in Table 1. The formation exhibits considerable shaliaess and thinning near the oil-water contact indicating that little or no benefit will be received from natural water influx. This is substantiated by the fact that wells completed near the original oil-water contact have not produced water. A sizeable gas cap is present; however, considerable faulting is in evidence, and it is believed that only a limited portion of the gas cap will aid in oil displacement and pressure maintenance. Although the reservoir was in an initial stage of depletion, it was believed that as soon as the equilibrium gas saturation in the oil zone was exceeded and a number of upstructure wells were invaded by the expanding gas cap, the resulting high gas-oil ratio production would rapidly dissipate the natural drive and lead to limited oil recoveries and inefficient reservoir operation. This analysis evaluated the effectiveness of the primary displacement mechanism and

Jan 1, 1958

By N. R. Morrow, C. C. Harris

The experimental points which describe capillary pressure curves are determined at apparent equilibria which are observed after bydrodynamic flow has ceased. For most systems, the time required to obtain equalization of pressure throughout the discontinuous part of a phase is prohibitive. To permit experimental points to be described as equilibria, a model of capillary behavior is proposed where mass transfer is restricted to bulk fluid flow. Model capillary pressure curves follow if the path described by such points is independent of the rate at which the saturation was changed to attain a capillary pressure point. Aemodified suction potential technique is used to study cyclic relationships between capillary pressure and moisture content for a porous mass. The time taken to complete an experiment was greatly reduced by using small samples. INTRODUCTION Capillary retention of liquid by porous materials has been investigated in the fields of hydrology,l soil science,2 oil reservoir engineering,3 chemical engineering," soil mechanics,5 textiles: paper making: and building materials In studies of the immiscible displacement of one fluid by another within a porous bed, drainage columnsl and suction potential techniques2 have been used to obtain relationships between pressure deficiency and saturation (Fig. 1). Except where there is no hysteresis of contact angle and the solid is of simple geometry, such as a tube of uniform cross section, there is hysteresis in the relationship between capillary pressure and saturation.2 The relationship which has received most attention is displacement of fluid from an initially saturated bed (Fig. 1, Curve R2), the final condition being an irreducible minimum fluid saturation Swr. Imbibition (Fig. 1, Curve A), further desaturation (Fig. 1, Curve R), and intermediate scanning curves have been studied to a lesser but increasing extent.2 ,3,9-12 This paper first considers the nature of the experimental points tracing the capillary pressure curves with respect to the modes and rates of mass transfer which are operative during the course of measurement. There are clear indications that the experimental points which describe these curves are obtained at apparent equilibria which are observed when viscous fluid flow has ceased; and any further changes in the fluid distribution are the result of much slower mass transfer processes, such as diffusion. Unless stated otherwise, this discussion applies to a stable packing of equal, smooth, hydrophilic spheres supported by a suction plate5 with water as the wetting phase and air as the nonwetting phase. The discussion can also be applied, in many

Jan 1, 1966

By F. F. Skiba, P. J. Root

A mathematical model has been used to investigate crossflow effects in a stratified reservoir during an idealized displacement process. A process in which one incompressible fluid displaces completely another incompressible fluid of the same density and viscosity was considered. The reservoir was assumed to be composed of homogeneous layers which are discrete but interconnected. Approximate five-spot flow geometry was simulated by arranging two pie-shaped cylindrical (two-dimensional) wedges in series. Some of the common production methods rely, for their effectiveness, on the fact that the formation consists of isolated strata. This is particularly true of selective plugging and single-zone production-injection methods. The effectiveness of these techniques for several cases in which adjacent strata are in communication has been evaluated. The usual production-injection schemes for fluid injection processes involve using the same completion interval in both the production and the injection wells. The efficiency of this, as well as alternative productioz-injection procedures has been examined for some cases involving communicating strata. INTRODUCTION The performance of many fluid injection projects is marred by early breakthrough and by production of a high percentage of the displacing phase after breakthrough. This channeling effect, if it occurs in reservoirs which have large permeability variations, is often explained on the basis of a model composed of discrete, isolated strata. Corrective or preventive measures based on this same model are also common. The behavior of stratified models has been the subject of considerable study. An excellent review of the early literature was given by Seba.1 Studies of the depletion performance have since been made for continuous, unrestricted communication between strata2-4 and for production from one permeable bed in continuous but restricted communication with an adjacent permeable bed.5,6 Model studies of displacement processes in stratified systems with crossflow also have been reprted.7-12 Finally, methods which account for crossflow have been developed for computing displacement behavior in stratified systems.13,14 Most methods for controlling or preventing channeling involve some change in the interval open to production and/or injection. Although many variations of the selective completion problem have been studied in regard to the depletion of a reservoir, little attention has been given to the behavior of a displacement process under similar conditions. Sebal used a resistance network model to study the performance of a process in which one incompressible fluid displaces an identical fluid from a linear, stratified reservoir. He investigated the possibility of improving the vertical coverage by limiting the completion interval of the production well. This investigation deals with techniques in which selective completion of the production and/or the injection wells is used in an attempt to control channeling. At one extreme, selective plugging and single-zone isolation (assuming perfect execution) should be successful if the reservoir is actually made up of independent strata. At the other extreme, they can hardly be expected to be successful if the channeling is caused by a mechanism in which stratification effects are negligible even though there are extreme variations in the formation permeability and/or the injectivity profile of the wells is not uniform. By contrast, it is difficult to anticipate the effect that the amount of communication between the beds of a stratified system will have on the efficiency of selective completion techniques. The objective of this study is the quantitative determination of the effect of communication on the behavior of some particular stratified systems which will at least lead to qualitative conclusions.

Jan 1, 1966

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

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

By W. A. Burns

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

Jan 1, 1970

By A. Timur

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

Jan 1, 1970

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

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

Jan 1, 1970

By H. Kazemi, G. W. Thomas, M. S. Seth

The double-porosity model of Warren and Root7 for examining pressure drawdown and buildup phenomena in naturally fractured reservoirs has been extended to intetpret interference test results. Both analytical and numerical solutions are presented. The practicality of some of the assumptions in the Warren-Root model was also investigated. The standard procedure for interpreting interference tests in naturally fractured reservoirs beretofore employed an equivalent homogeneous reservoir model. A comparison of our results with the equivalent homogeneous model indicates the latter is inadequate for early time responses. For late time responses, however, the simpler equivalent homogeneous model yields favorable results. Finally, the effects of well spacing and the magnitude of the interporosity flow parameter are briefly cited. INTRODUCTION Interference tests can be used to provide reservoir information not obtainable from conventional pressure drawdown and buildup tests. The technique involves measuring the response in one or more wells arising from a pressure perturbation introduced into another well. It has been. successfully applied to single-layer, homogeneous, unfractured, isotropic and anisotropic formations to arrive at values of well storage capacity and a measure of directional permeability.l The techniques of interference testing usually can be extended to heterogeneous systems, provided the permeability contrast between any existing rock layers is not excessive.9 In such cases, the practical interference effects in an observation well of an interconnecting heterogeneous reservoir can be predicted using an equivalent homogeneous model. The calculated properties are then average values. Russell and Prats2 have demonstrated the utility of this approach for a single-phase, two-layer system. For naturally fractured reservoirs, some different approaches usually are necessary. The major departure from homogeneous or nearly homogeneous reservoir treatment is to characterize the fractured system as an entity having "double-porosity"; namely, we suppose the existence of both primary and secondary porosity, where the former is contained within a systematic array of identical rectangular matrix blocks and the latter is contained within an orthogonal system of continuous, uniform fractures (see Fig. 1). It is sufficient to assume that the permeability of the rock matrix is so small that movement of reservoir fluids to the producing wells is entirely through the fracture network. This assumption has been validated by Barenblatt3 and .employed by others?-7 Warren and Root7 used the double-porosity model to develop an interpretation of pressure drawdown and pressure buildup tests in naturally fractured

Jan 1, 1970