Search Documents

By J. J. Arps, T. G. Roberts

Since the introduction of the relative permeability concept in the middle thirties1,2 various investigators have shown3,10,11,12,15 how the basic equations for the flow of oil and gas through porous media can he utilized to compute the recovery from depletion type or solution gas drive type reservoirs under certain conditions, when the necessary physical data pertaining to the reservoir rock and its reservoir fluids are known. In the study by Muskat and Taylor" in 1945, the effects of viscosity, gas solubility, shrinkage, and the permeability-saturation characteristics of the producing formation on the production histories and the recovery of gas drive reservoirs were analyzed. Each parameter was varied a limited number of times, while keeping the others constant. The work was not carried to the point where the data could be used directly to estimate the primary recovery in those cases where certain rock and fluid characteristics were known or could be assumed. Since the publication of this classic study, the work by Beal in 194616 and Standing in 1952'" has established general correlations of oil viscosity and shrinkage with gas solubility and oil gravity. In addition. a considerable number of relative permeability relationships for different types of reservoir rocks has been published2,5,13,17,19,20. Also, the burden of the numerical work required in the solutions of the differential equation has. during the last few years, been greatly reduced with the availability of modern electronic computing facilities. The purpose of this paper is therefore to present the results of a large number of such recovery computations covering the normal range of variation of the main parameters such as the type of rock, the API gravity of the oil, and the amount of gas in solution. As is to be expected, the ultimate recovery is found to increase with the oil gravity. except for the higher solution gas-oil ratios. Within the normal range of solubilities, the ultimate recovery generally appears to decrease with the amount of gas in solution for the higher oil gravities, and to show a slight increase for the lower oil gravities. Intermediate oil gravities seem to fall between the two trends. The type of reservoir rock, as identified by its relative permeability relationship, appears to have a very pronounced effect on the ultimate recovery. In general, sands and sandstones show higher recoveries than limestones and dolomites, although certain intergranular limestones seem to show a higher theoretical recovery than unconsolidated sand. The recovery from sands and sandstones generally decreases with increasing cementation and consolidation, while the recoveries from limestones and dolomites are highest for the intergranular type and lowest for the fracture type porosity.

Jan 1, 1956

By W. R. Foster, J. M. McMillen, A. S. Odeh

The complete equations of average linear momentum balance for a single-phase fluid in an incompreisible, homogeneous, porous medium are derived. The derivation begins with Euler's equation of motion for a continuum and uses an integral transform recently developed by Slattery. 7 For steady flow of a compressible, Newtonian fluid, the usual equations of motion result. For transient flow, the space-time description of the pressure is determined in the lowest approximation by the telegrapher's equation. From the analysis a new phenomenological coefficient results which connects the viscous traction to the deriv ative of the linear momentum density. The magnitude of this coefficient determines the velocity of sound through the pore structure in this approximation to the pressure field. INTRODUCTION The modification of Darcy's law of momentum balance for steady, single-phase flow through porous media has been discussed for many years. The first modification was suggested by Forcheimerl who added terms of higher order in the velocity. These can be expected to appear because the underlying microscopic equations of momentum balance are themselves nonlinear in the point velocity field. The Reynolds tensor pvv, which represents the convective flux of momentum density, appears in the momentum balance equation. Only in rectilinear flow (parallel stream lines) does the divergence of this tensor vanish. Since the steady flow stream lines in most porous media are not parallel, nonlinear dependence of the pressure gradient on the velocity should naturally appear. This nonlinearity has nothing to do with turbulence in the ordinary sense of random fluctuations in the pressure and velocity fields. It arises simply because the stream lines converge and diverge, even for steady flows. Klinkenberg2 demonstrated that the permeability coefficient in Darcy's law depends on the absolute pressure or, alternatively, on the density field. However, because he neglected inertial terms of the Forcheimer type, his correction coefficient could not be represented by a constant but tended toward a constant as the velocity decreased. Forcheimer's and Klinkenberg's modifications can be combined in a rigorous way to account for both inertia and slip during steady flow. This will be shown in a future paper. The transient change of pressure in porous media has been described by the diffusion equation.3 This form results from eliminating the velocity and density fields from a combination of the equations of motion in the form of Darcy's law, the continuity equation and an equation of state. Fatt4 suggested that the cause of deviations from the prediction of the diffusion equation for pressure transients lies not in the choice of Darcy's law as the equation of motion but on the existence of dead-end pores which might invalidate the averaged equation of continuity. On the other hand, Oroveanu and Pascal5 noted chat the time derivative of the momentum density must be included in the equations of motion since this measures the local rate of change of momentum density. Their differential equation for pressure is the telegrapher's equation (neglecting gravity). However, their form of this equation predicts that the speed of pressure propagation through the pore structure is the same as that through the bulk fluid. M. K. Hubbert6 attempted a derivation of Darcy's law by volume averaging the Navier-Stokes equations. Since these equations represent momentum balance at a point within an open set of points containing the fluid itself, Hubbert's volume averaging cannot lead to tenns involving transfer of momentum between the fluid and the walls of the pores. Once these viscous tractions are lost by choosing a control volume containing only fluid, they cannot

By M. Prats

Most of the information available on the heat efficiency of hot fluid injection processes, both water and steam, has been obtained from calculated temperature distributions in the pay zone and adjacent formations. The usual approach has been to write the heat balance equations in terms of the temperatures, and then to introduce whatever simplifications are necessary to help obtain an analytical or a numerical solution. In either case, the assumption has often been made that the vertical thermal conductivity in the flooded interval is infinite, so that the temperatures within the flooded interval are then independent of vertical position. Because it was first made by Lauwerier,l and has been used extensively since, this assumption will be referred to as the Lauwerier assumption. Excellent reviews of the literature on the heat efficiency of hot fluid injection processes have been given by Spillette,2 Flock et al.,3 and Ramney.4 Of course, the most significant contributions take into account the effect of a finite vertical thermal conductivity on the vertical temperature profile. The most general analytic expression for the heat efficiency of a hot fluid injection process is that of Antimirov5 who considers injection of a heated incompressible fluid into a reservoir through an arbitrary number of wells. The rate of heat injection into the reservoir, as well as the injection temperature, is an arbitrary function of time. The geometry of the horizontal flow is arbitrary, although the reservoir is considered to be of uniform thickness and properties and to be of infinite areal extent. Heat transfer within the reservoir is by horizontal convection and conduc- tion, and by vertical conduction. In the formations adjacent to the flooded interval, heat transfer is by conduction in any direction. Thus, the available expression for the heat efficiency due to hot water injection has been developed for rather general conditions. This degree of applicability has not been obtained for other thermal recovery processes. By making the Lauwerier assumption we have been able to obtain expressions for the heat efficiency that are essentially independent of the thermal recovery process, be it steam, hot-water, or underground combustion. Clearly, then, the approach followed here is more restrictive than that of Antimirov5 in that it uses the Lauwerier assumption. At the same time less restrictive assumptions are made about the horizontal heat transfer mechanisms, either in the flooded interval or in the formations adjacent to it. The main difference between our approach and that of other investigators is that the heat balance in the pay zone is expressed for the pay zone as a whole rather than for a volume element within it. The Lauwerier assumption is introduced after the problem is developed along these lines as far as possible. The details of the development of the heat efficiency are given in the Appendix, while the main features of the model are discussed in the next section. Results and their implications are discussed later in a separate section. The Model Certain assumptions have been made to determine

Jan 1, 1970

By D. Havlena, A. S. Odeh

The use of the straight-line method of solving the material balance equation is illustrated by means of six field cases. Also, the application of statistical criteria to arrive at the most probable answer is shown. The theory underlying the straight-line method of solution and the applicability of the statistical criteria was presented in a previous paper.' The field cases include saturated and undersaturated oil reservoirs with and without water drive. The aquifers discussed are: limited radial, infinite radial, very small aquifer and infinite linear. The field cases also include a gas reservoir producing under water drive. INTRODUCTION In a previous paper,' the authors presented the theory underlying the solution of the material balance as an equation of a straight line. The appropriate equations for various material balance cases as well as the methods of analysis and interpretation with comments and discussion were also included. To illustrate the various theoretical cases treated previously, five field cases are analyzed in this paper by employing the straight-line method of solving the material balance equation (MBE) and one example previously published is referred to. The use of statistical criteria to arrive at the most probable answer is also shown. All the field examples presented, except Case 2, are excerpts from complete reservoir studies. To illustrate the method, only sections specifically dealing with the material balance principles are included. Additional geologic information and basic data are reported to better acquire an understanding of the cases and thus to better follow the reasoning that suggested the successful application of the straight-line method of solving the MBE. The six cases are: (1) saturated reservoir, small gas cap, limited aquifer; (2) saturated reservoir, very small gas cap, infinite aquifer; (3) undersaturated-saturated reservoir, very small aquifer; (4) highly undersaturated reservoir, no water drive; (5) high undersaturated one-well reservoir, limited aquifer; and (6) gas reservoir, infinite linear aquifer. THE D4 SAND, GUICO FIELD, VENEZUELA The D, sand, which was discovered in 1943, is presently in a depleted state. Since its discovery it has produced under water drive, gas-cap-gas expansion, and solution gas drive. In Nov., 1947, water injection was initiated to arrest further pressure decline. When discovered, the D, sand was a saturated reservoir with a gas cap /oil zone volume ratio m estimated volu-metrically at 0.0731, an average permeability of 500 md, a porosity value of 25 per cent, and an oil viscosity at reservoir conditions of 0.3 cp. The volumetrically determined stock-tank oil initially in place was 23.1 million bbl. The volumetrically weighted physical data and production data available until Nov., 1953 are reported in Table 1. In Ref. 1, the effects on the straight-line plot of various values of ro/r for a constant ?to, or of various dimen-sionless times for a constant r8/r10 were theorized and were illustrated in Fig. 3A of that reference. In this field case, the previously theoretically predicted effects are established. Thus, the MBE calculations using Eq. 3c of Ref. 1 were performed for various r6/r10 and dimension-less time values. Eq. 3c of Ref. 1 is:

Jan 1, 1965

By J. E. Berry

AII evaluation is made of the acoustic velocity log for measurement of formation porosity. Plots of field-observer1 velocities vs core-measured porosities of sandstones and limestotnes with inter intergranular porosity show that the velocity log gives a useful measure of porotity, in agreement with published data. A new method of plotting electrical and velocity log data makes it easier to recognize hydrocarbon-water coritnct~ from log darn and to make serni-quantitative estinlates of hydrocarbon saturation. Two examples of this interpretation technique are given. This paper also discusses the various factors that affect acoustic velocity and shows how corrections can be made for some of these factors. INTRODUCTION Although the continuous velocity log was developed primarily as an aid to seismic interpretation, it is finding widespread use for measurement of porosity. The purpose of this paper is to evaluate this use of the log. Results are presented for the velocity-porosity relationship in many sandstones and one limestone. Practical applications of porosity values derived from the velocity log to improve interpretation of electrical surveys are presented. Various parameters which affect velocity in porous media are considered. Porosity, composition, cementa-tion, pressure difference (overburden pressure minus fluid pressure), fluid saturation, and wellability are discussed. Some of the ideas are opinions based on general considerations of well logging problems while others are supported by published data and our own observations. Major conclusions from these studies arc given in the body of the paper and detailed discussions of some of the various factors affecting porosity are given in the Appendix. POROSITY Sandstones The primary factor affecting the velocity of sound in porous media is porosity. Wyllie, et al,1 have reported a time average equation for the relationship between velocity and porosity.* This has found considerable acceptance in the industry. Other relationships between velocity and porosity have been proposed, none of which are entirely satisfactory. It is doubtful that any one relationship will completely describe all situations. We have tried an empirical relationship based on correlating log velocities with core-measured porosities. Measured porosities were averaged over a depth interval equal to the spacing used in recording the velocity logs. The depth range for the information used was from 2,000 to 12,000 ft. Velocity measurements were normalized for the effects of depth and overburden pressure. A detailed discussion of the basis for this normalization is given in the Appendix. The change is usually insignificant and when working with individual formations it can be neglected. All velocity measurements were normalized to a depth of 10,000 ft. The velocity-porosity data for sandstones in Fig. 1 came from 20 wells scattered throughout the Mid-Continent and Gulf Coast areas. The scatter of the data points is attributed to three causes: 1. There is no theoretical basis for expecting an exact velocity vs porosity relationship and it is not possible to correct precisely for the effects of such factors as impurities, grain shape, manner of cementation, etc. Some of these factors are discussed in the Appendix. 2. Porosities measured by core analysis are not necessarily representative of what the velocity log "sees". 3. Inaccuracies appear in the velocity log reading. This is particularly true of velocities from one-receiver logging tools.*" A linear-time curve (time-average equation) is plotted in Fig. 1 along with the best straight line through the points. Velocities used for the linear-time curves were 19,500 ft/sec for the matrix velocity (zero porosity) and 5,500 ft/sec for fluid velocity. For clean sandstones with velocities greater than 10,000 ft/sec, there is little difference between the velocity-porosity and the time-porosity curves. Both curves are well within the scatter of data points. The data points in Fig. 1 were restricted primarily to clean sandstones, using the SP and/or gamma-ray logs as clues to the lithology. Where a lithologic de-

By W. F. Kieschnick, and W. J. McGuire&apos, Eugene Harrison

This paper concerns the induction and extension of fractures into rock formations as involved in drilling, completing, and production stimulating operations on wells. Conclusions concerning formation breakdown are derived from (I) a review and extension of published analyses relating to mechanical theories of rock stress and the state of stress in the earth's crust and (2) a correlation of field data from fracturing operations. Conclusions concerning the mechanics of fracture extension, which indicate the relationship between fracture dimensions and rock properties, depth, and volume of injected fluid, are tentative and largely establish limits of relationships. These conclusions are derived from stress calculations, limited field data, and laboratory experimental studies. The experimental work involves the study of the stresses at the fracture boundaries and the geometry of pressurized fractures by means of photo-elastic modeling methods. Results of this investigation indicate that a large majority of pressure induced wellbore fractures are vertical, particularly in deeper wells; and variations in the pressures necessary to create and extend fractures can be explained largely on a basis of established rock properties. It is also shown that variations due to tectonic forces should usually be expected to be slight. Other results indicate that during that extension of fractures rather large fracture volumes are temporarily created by the parting of the formation. The purpose of this paper is to present the results of calculations and laboratory experiments concerning the mechanics of fracture induction and extension with a view to broadening existing knowledge relating to these phenomena. It is believed that continued progress in developing knowledge of this type is important to the further development of techniques for drilling and completing wells. When a wellbore is subjected to sufficiently high pressure, fracturing of the formation surrounding the well occurs. This phenomenon is a common cause of lost circulation during drilling operations, is a deliberately planned factor in most squeeze cementing operations, and is the controlling factor in many productivity stimulation operations. The exploitation of formation fracturing in recent years has been particularly impressive in sand-viscous fluid well stimulation treatments which were first introduced to the industry by the Stanolind research group1,2,3 several years ago. The borehole fracturing phenomenon has been analyzed mathematically by application of the principles defining the elastic and inelastic behavior of thick-walled cylinders of homogeneous, isotropic, and impermeable material.1,5,6 This basic argument is re-presented in this paper along with new considerations which strengthen it. In addition to these arguments relating to fracture induction some calculations and experimental results relating to the mechanics of fracture extension are presented. Since the nature of the pressure induced fracture at distances from the well is of great importance in considerations of well stimulation procedures there appears to be a need for an integrated approach which deals with the mechanics of fracture extension as well as with wellbore rupture.

Jan 1, 1955

By H. T. Kennedy, S. M. Avasthi

An equation developed for gaseous hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 0.73 percent when applied to 264 natural gas and condensate systems including 2,043 PVT points. Another equation developed for liquid hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 1.12 percent when applied to 346 crude oil systems including 1,759 PVT points. Both equations require composition of the mixture to be expressed as mole fraction of methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide, together with the characteristics of the heptanes-plus fraction in addition to the temperature and pressure. The equations cover wide ranges of the variables involved, and their accuracy is considerably better than that of other available methods. The equations were differentiated to allow calculation of the coefficients of isothermal compressibility and isobaric thermal expansion. (In this paper the coefficient of isothermal com-pressibility and the coefficient of isobaric thermal expansion will be expressed as compressibility and thermal expansion coefficient, respectively.) Equations to calculate these quantities are presented. INTRODUCTION Calculations of reservoir performance for petroleum reservoirs require accurate knowledge of the volumetric behavior of hydrocarbon mixtures, both liquid and gaseous. Compressibilities are required in transient fluid flow problems, and thermal expansion coefficients are important in thermal methods of production. An accurate laboratory investigation of the PVT behavior of each reservoir fluid encountered would be costly and time consuming. For this reason various correlations for predicting fluid properties have been developed and recorded in recent literature. Correlations have been presented in the form of graphs, tables and equations. Since an increasing number of studies are being conducted with the aid of electronic computers, recent efforts have been directed toward development of correla-tions suitable for computer programming. Application of computers permits the use of more complex correlations which otherwise are not feasible. Moreover, methods for predicting reservoir performance, particularly those based on the compositional material balance, depend upon the capability of accurately expressing the molal volumes and other fluid properties as functions of pressure, temperature and composition. The coefficient of isothermal compressibility c is defined by

Jan 1, 1969

By A. F. van Everdingen

The pressure drop in a well per unit rate of flow is conrolled by the resistance of the formation, the viscosity of the fluid. and the additional resistance concentrated around the well bore resulting from the drilling and completion technique employed and, perhaps, from the production practices used. The pressure drop caused by this additional resistance is defined in this paper as the skin effect. denoted by the symbol S. This skin effect considerably detracts from a well's capacity to produce. Methods are given to determine quantitatively (a) the value of S, (b) the final build-up pressure, and (c) the product of average permeability times the thickness of the producing formation. INTRODUCTION Equations which relate the pressure in a well producing from a homogeneous formation with pressures existing at various distances around the well are generally used within the industry. The relation ii quite simple when the fluid flowing is assumed to be incompressible. It becomes somewhat more complicated when the flowing fluid is considered compressible so that the duration of the flow can he considered. In each case the major portion of the pressure drop occurs close to the well bore. However analyses of pressure build-up curves indicate that the pressure drop in the vicinity of the well bore is greater than that computed from these equations using the known, physical characteristics of the formation and the fluids. In order to explain there excessive drops it is necessary to assume that permeability of the formation at and near the well bore is substantially reduced as a result of drilling. completion and, perhaps. production practice. This possibility has been recognized in the literature. A method to compute the pressure drop due to a reduction of the permeability of the formation near the well bore. which is designated as the skin effect. S, is given in the following paragraphs. To start, equations normally used to describe flow in the vicinity of a well are given without considering this effect. These equations then are modified to include the effect of a skin on the pressure behavior. Finally a method is given to estimate the effect of the skin on the pressure and production behavior of a well. PRESSURE EQUATIONS Incompressible Fluid Flow If p is defined as the flowing pressure in a well of radius the pressure at distance r from the well has been shown to be:" The total pressure drop between the drainage boundary, and the well bore is given by These equations are valid only if the flow towards the well occurs in a horizontal homogeneous medium and the fluids are incompressible. The assumptions imply that all fluid taken from the well enters the system at r a condition rarely encountered in practice. Compressible Fluid Flow, Steady State A more realistic equation is obtained if it is assumed that the compressibility, c, of the flowing fluids is small and has a constant value over the pressure range encountered. After the well has been producing for some time so that its rate has become constant and steady state is reached, the pressures throughout the drainage area are falling by the same amount per unit of time, and the pressure differences between a point in the drainage area and the well are constant. When these conditions are met. the rate of production, q, from a well is equal to where dp/dt is the pressure drop per unit time. The fluid flowing at a distance from the center of the well is equal to From the last equation and from Darcy's law it can be shown that The equation holds for a depletion-type reservoir of radius drained by a well located in its center, provided the compressibility of the fluid per unit pressure drop is small and constant, and no fluid moves across the boundary Compressible Fluid Flow — Nonsteady State Table 111 of reference (5) shows the relationship between the pressure at the well bore and the reduced time, The pressure-drop function, p represents the drop below the original reservoir pressure, p caused by unit rate of production for several values of R, the ratio of drainage boundary radius to well radius, r In most reservoirs the values of approach infinity. and under these conditions the values of p shown in Table I of reference (5) can be used where p then signifies the difference between the pressure in the well and the prevailing reservoir pressure per unit rate of flow. The total pressure drop below prevailing reservoir pressure amounts to where the factor converts the cumulative pressure drop per unit rate of production to cumulative pressure drop for actual rate. q. For values of T > 100 the P function may be written (equation VI-15 of reference 5) as Using the time conversion the difference in pressure between reservoir and well becomes If values for the physical constants of the formation and the fluids are inserted, it is found that T exceeds 100 after a few seconds of production (or closed-in time), so that the approximation becomes valid almost at once. A simple relation between the pressure in the well and in the reservoir can also be derived by considering the well as a point source"" '" instead of a unit circle source, that is, by using Lord Kelvin's solution instead of the unit circle source

Jan 1, 1953

By S. S. Marsden, S. H. Raza

An experimental study of the flow of line-textured. aqueous foams through Pyrex tubes is described. The foams range in quality F (ratio of gas volume to total volume) from 0.70 to 0.96 and behave like pseudoplastic fluids. At lower flow rates they exhibit laminar flow ad have apparent viscosities which increase with quality, and which cover a range of 15 cp to 255 poise for tubes of 0.25- to 1.50-mm radius rj. At higher flow rates a plug-like type of flow is developed, the extent of which increases with both T and ri. When the same foams flow through either open or packed Pyrex tubes, remarkably high streaming potentials ?E are often generated. These can easily reach 50v if nonionic foaming agents are used, but are at least an order of magnitude less for ionic foaming agents. A linear relationship between ?E and the pressure differential ?p is observed; this usually extrapolates to positive values of ?p at ?E of zero. The slope of the line increases with both F and ri. An equation was derived to describe the streaming potential of non- Newtonian fluids in circular tubes and was used to correlate experimental results. The calculated < potentials are of the right order of magnitude. INTRODUCTION Foams are both unusual and inniguing in their physical properties, and have been the subject of many scientific studies. However, present knowledge of foams is still fragmentary, specific and often contradictory. Apparent viscosity of foam is the physical property of greatest interest to both rheo1ogists and engineers. sibreel reported that the apparent viscosity decreased with increasing shear rate in a manner similar to some non-Newtonian fluids. Penny and Blackman2 reported that fire-fighting foams had both a limiting shear stress and a tensile yield stress. There is little doubt that some foams at least behave like non-Newtonian fluids, and have apparent viscosities considerably higher than those of either constituent phase. The high apparent viscosity of foam with its concomittant effect on mobility ratio and sweep efficiency no doubt prompted several attempts by research groups to use foam as a displacing agent in porous media. Based on recent experience, most of these groups probably succeeded in completely blocking fluid flow in the porous media and then abandoned their efforts. Two groups apparently found the successful combination of experimental parameters at about the same time.3,7 Others have recently added to our knowledge of foam flow in porous media and its use as a displacing agent.8-13 An experimental problem encountered by Fried3 was a transient blockage of foam flow in porous media when distilled water was used to prepare the foam-producing solution. Fried surmised that this was due to an electrokinetic effect and he eliminated it by using electrolytes in preparing foaming solutions. He also measured the streaming potential of a number of foams in capillary tubes which he found to be appreciably higher than those obtained when the constituent liquid flowed under comparable conditions. This paper presents results of a more comprehensive study of the streaming potential generated by aqueous foam flowing in both open and packed Pyrex tubes. It also adds to knowledge of the rheology of these foams as deduced from their flow behavior in open tubes. APPARATUS AND PROCEDURE A diagram of the apparatus used is shown in Fig. 1. Details of its construction, testing and use are described elsewhere.14 Careful selection of materials, extreme cleanliness and rather elaborate electrical insulation and shielding were necessary to obtain reproducible results (±5 percent). Both streaming potential and streaming current were measured with an electrometer. The design of the foam generator developed for this work is novel (Fig. 2). It will produce uniform fine-textured foams at pressures sufficiently

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By E. C. Barfield, D. G. McCarty

Two-dimensional analyses offer considerable promise in providing the basic information required to effect more precise control of petroleum reservoir performance. This paper describes a method for conducting some engineering analyses of this type using a high-qreed digital computer. The general approach is to (I) develop a representative numerical model of the reservoir and (2) employ a suitable numerical technique to solve the basic equation of flow and to develop the required engineering information for the particular case represented by the model. This technique has been used to calculate pattern performance in connection with several field projects involving water injection into oil reservoirs. This type of analysis involves sequentially (I) calculating the pres-&apos; sure distribution in the reservoir and (2) tracking the progress of the interface between the displaced and displacing fluids in a step-wise manner to provide a depletion history of the operation. The particular analysis presented in this discussion is subject to the restriction that the mobility of the displaced and displacing fluids be equal and assumes that the fluids are incompressible and that gravity and capillary pressure do not affect the shape of the flood pattern. Complete two-dimensional .flexibility is maintained with regard to definition of reservoir rock and fluid characteristics, placement of physical restrains and boundary conditions, investigation of flow characteristics in the reservoir, and movement of the displaced and displacing fluid interface. The results of these studies indicate that the highspeed digital computer is well suited for conducting reservoir performance studies in two-space dimensions and that dependable numerical techniques are available for making such analyses. INTRODUCTION One of the problems facing petroleum engineers is that of providing more precise engineering control over reservoir operations at all stages during the depletion cycle. The most direct approach to this problem is to simulate the reservoir with some type of physical or mathematical model and use the information developed from the behavior of the prototype to predict the performance of the actual reservoir. This principle is, of course, not new, and a variety of simulation techniques have been developed to provide information on the behavior of oil and gas reservoirs. One straightforward refinement of present reservoir analysis techniques is the progress from one-dimensional reservoir models to those where the reservoir can be represented in two- and three-space dimensions. Multidimensional reservoir analyses are desired in order to show not only what is happening in the reservoir but in addition to show where in the reservoir the various phenomena of interest are occurring. Thus far, two-dimensional analyses have been confined predominantly to potentiometric and electrolytic model studies. However, when large numbers of wells are involved and where large-scale models are required, such problems become very tedious and instrumentation problems become acute." It is possible through use of numerical methods and general purpose digital computers to avoid many of the instrumentation difficulties and to relax some of the restrictions of previous techniques. The purpose of this paper is to discuss some recent experiences with the use of high-speed digital computers for predicting pattern performance in fluid injection operations and to describe the numerical procedures which have been used for making such analyses. THE NUMERICAL MODEL The problem of representing the reservoir with an appropriate numerical model is composed of at least two important considerations. The first requires that an adequate description of the reservoir mechanism be formulated in the form of the basic differential equation describing the process. The second consideration is definition of the physical characteristics of the reservoir. The latter requires that representative basic data be obtained, examined, and interpreted to provide sufficient detailed definition of the pertinent physical properties throughout the system. The Differential Equation One of the basic assumptions made for the purpose of performing the two-dimensional studies being dis-

By F. J. Fayers, R. L. Perrine

For many years the problem of increasing ultimate recovery of oil from a reservoir has been a subject of interest to the oil industry. At present, a standard secondary recovery technique is to flood the reservoir with water once natural forces of the system are depleted. The permeability characteristics of the rock, together with interfacial forces which exist between the two fluids and the porous material. allow only part of the oil in the rock formation to be recovered by the water-flooding process. Laboratory experiments have shown that when detergents are added to the flooding water, the flow properties of oil in a porous rock arc improved. This improvement results in greater oil recovery. Even though higher recoveries are obtained using detergents, the extra cost associated with the method may make it economically unsatisfactory. To establish economic feasibility requires experimental laboratory measurement of fluid and rock properties together with a calculational procedure to predict field results. Pilot tests in the field may also be required. This paper presents a mathematical procedure for predicting detergent flood behavior. A number of reasonable simplifying assumptions are made to obtain a tractable system of equations. The resulting method may be regarded as a generalization of Buckley-Lev-crett&apos; theory to solve simultaneously for water saturation and detergent concentration profiles. Consequently this type of calculation can be used in much the same manner as Buckley-Leverett theory is used to evaluate a water flood. The mathematical method employed is the so-called method of characteristics for solving hyperbolic partial differential equations. This application of the method is instructive because several other problems in reservoir analysis may be analyzed in a similar manner. FLUID FLOW EQUATIONS The flow behavior of oil and water without detergent through porous reservoir rock is usually expressed by Darcy&apos;s law and a continuity equation for each phase. If gravitational forces are neglected, vrg = kr (s)µr v p rg.........(1) vrg = kr (s)µr v p rg.........(1) ? pa va = F ?/?t (sd) where v is the volumetric flux, k represents permeability, S is saturation, and P the pressure. Fluid viscosity is denoted by 11, p is the density, 4 is porosity, and I the time. The subscripts w and o refer to water and oil. The capillary forces in the system lead to a difference in pressure between phases. This may be expressed by the relation: Pm -P0 =Pc(S).........(5) The capillary pressure? P,, is determined experimentally as a function of saturation. DETERGENT BALANCE EQUATION Addition of detergent to the water phase affects the flow properties of the system. Capillary pressure and relative permeability become functions of detergent concentration, c, as well as fluid saturation. That is: P--P (S,c) and k,,,.- - -- k,,,. At low oil saturation values, the presence of detergent reduces capillary pressure and reduces the relative permeability ratio. Detergents are surface-active and adsorb on rock surfaces. The quantity adsorbed, a, will depend on the concentration c of detergent solution and may also depend on the saturation; thus, a = a(S,c). In this analysis we have assumed that the detergents used are water soluble but not oil soluble, and their movement is therefore restricted to the water phase. The dispersion of detergent in water is neglected by comparison with transfer by flow within the system. Detergent movement is determined by the water flux and adsorption: ? cv m = -F ?/?t (cs + a) We will assume that the quantity of detergent in solution is too small to significantly change water density or viscosity. When the permeability and capillary pres-

By C. Chu

A theoretical investigation has been made of the forward combustion process using a one-dimensional linear mathematical model, taking into consideration the effect of the vaporization-condensation which occurs on the leading edge of the heat wave. This work involves the solution of five coupled partial differential equations. Besides the vaporization-condensation phenomenon, these equations account for conduction, convection, combustion, heat loss, diffusion and bulk fluid flow. For the one-dimensional linear model studied, the vaporization-condensation phenomenon does not induce appreciable change in the temperature at the combustion front; and its primary effect is to create a steam plateau and to increase the length of the heated zone ahead of the combustion front. This effect becomes more pronounced at lower pressures, higher porosities or reduced gas saturations. The peak temperature and the temperature profile on the leading edge of the heat wave stabilize after a certain period. The length of the steam bank remains practically constant, although the length of the water bank increases as the heat wave advances. INTRODUCTION The existence of the vaporization-condensation phenomenon in the heat-wave process and the important role played by the phenomenon have been recognized by several investigators. Kuhn and Kochl stated that steam plateaus were frequently observed on the temperature records of their experiments. The steam plateaus were attributed primarily to the vaporization and subsequent condensation of the interstitial water existing in the oil sand. Szasz2 suggested that both lighter hydrocarbons and water are vaporized on the leading edge of the heat wave, carried forward in the gas stream, and then condensed to create banks of oil and water. Martin et al suggested that the vaporization-condensation phenomenon is one of the main mechanisms of the heat-wave process, along with thermal expansion and viscosity reduction. Wilson et al. reported the existence of a steam plateau several inches in length in their small-scale tube-run experiments. However, this important phenomenon has never been taken into consideration in the numerous theoretical analyses by various authors. 5 The purpose of this work was to study the thermal aspects of a linear heat wave, taking into consideration the vaporization-condensation on the leading edge of the wave, to determine the effect of this phenomenon on the temperature profile of the reservoir, and to investigate how this effect varies when other process variables are changed. THEORY We consider a reservoir of porous medium of cross-sectional area A, extending from x = 0 to x = L. This reservoir contains, aside from the solid matrix itself, a gas phase and a "combined liquid phase" which is a combination of two immiscible liquid phases — namely, an oil phase and a water phase. The oil present in the reservoir is assumed to consist of three fractions, a noncondensable gas, a nondistillable residuum, and a vaporizable oil fraction which may be present in both the gas phase and liquid phase. Before the heat-wave process begins, preheating has taken place and has imparted an initial temperature distribution T,(x) to the reservoir. At the start of the process, a stream of oxygen-bearing gas is introduced through the face at x = 0. This gas supports the combustion of the residual fuel and supplies the heat throughout the process. In order to develop a tractable system of equations, the model properties are postulated as follows: (1) the system is isotropic and homogeneous; (2) there is a uniform fuel content throughout the reservoir; (3) combustion takes place in a moving vertical planar region of infinitesimal thickness; (4) the temperatures of the fluids and the solid matrix are identical; (5) radiation is negligible; (6) the gas and liquid saturations are uniform throughout the reservoir and remain unchanged

Jan 1, 1965

By I. Fatt, D. N. Saraf

A method is described for measuring two- and three-phase relative permeabilities in sandstones or sand packs using a nuclear magnetic resonance (NMR) technique to determine fluid saturations. Two- and three-phase relative permeabilities have been determined on Boise sandstone using the NMR technique of saturation measurement. Three-phase relative permeability to water was found to depend only on the water saturation, whereas three-phase permeability to oil depended on both the water and oil saturations. Relative permeability to gas in three-phase flow was found to depend only on the total liquid saturation. INTRODUCTION Three-phase relative permeabilities are extremely useful in calculating field performance for reservoirs being produced by simultaneous water and gas drives. Three-phase relative permeability data are also needed for analyzing solution gas-drive reservoirs which are partially depleted and are being produced by water drive. Some thermal recovery processes involve three-phase flow which require three-phase relative permeability data for predicting reservoir behavior. Unfortunately, three-phase relative permeability measurements have rarely been made. Also, because of the scarcity of three-phase data, it has not been possible to date to relate other measured rock characteristics to the relative permeabilities with a great certainty. Leverett and Lewis,l Reid2 and Snell3 have reported three-phase relative permeability data on unconsolidated sands. Leverett and Lewis used ring electrodes spaced along the length of the sand sample to measure the resistivity of the sample which was assumed to be monotonically related to brine saturation. Gas saturation was determined from pressure-volume measurements. Oil saturation was obtained by material balance on the cell containing the sand sample. This method is involved and time consuming. Another difficulty arises from the fact that the resistivity of the sand is a function not only of saturation of brine but also of the distribution and saturation history of the brine in the pore spaces. Reid2 used a gamma ray absorption technique for measuring liquid saturation. This method has the disadvantage that total liquid saturation rather than oil or brine saturation is all that can be measured and still another method is required to determine the saturations of individual components. snel13 used a neutron bombardment method which also required a separate determination of the individual component saturations. Caudle et al.4 measured three-phase relative permeability on consolidated sandstones using vacuum distillation for determining fluid saturations. Distillation after each reading makes this technique very lengthy and time consuming. Corey et a1. 5 and Naar and Wygal 6 measured three-phase relative permeability on sandstones by the capillary-pressure method. Semipermeable diaphragm assemblies were used at each end of the core specimen to keep the water phase in the core. Gravimetric methods were used to determine fluid saturations. Sarem7 recently reported an unsteady-state method for measuring three-phase relative permeability on sandstones. This method is an extension of Welge&apos;s method8 for measuring two-phase relative permeability. Although Sarem&apos;s method is simple and comparatively fast, the assumptions involved may oversimplify the problem. Sarem&apos;s assumption, that in all rocks relative permeability to each fluid will depend only on the saturation of that fluid, seems to be rather unrealistic. Neglecting capillary effects at the end of the core is also a weak assumption. Donaldson and dean9 measured three-phase relative permeability using a method similar to Sarem&apos;s.

By A. M. Sarem

For the performance prediction of multiphase oil recovery processes such as steam stimulation, there is an acute need for three-phase relative permeability data. No fast and simple experimental technique, such as the unsteady state method proposed by Welge for two-phase flow, is available for the three-phase flow. In this paper, an unsteady - state method is presented for obtaining three-phase relative permeability data; this method is as fast and easy as Welge&apos;s method for two - phase flow. Analytical expressions are derived by extension of the Buckley-Leverett theory to three-phase flow to express the saturation at the outflow face for all three phases in terms of the known parameters. It is assumed that the fractional flow and relative permeability of each phase are a function of the saturation of that phase. Other simplifying assumptions made include the neglect of capillary and gravity effects. The effect of saturation history upon relative permeability is acknowledged and attainment of similar saturation history in laboratory and field is recommended. The required experimental work and computations are simple to perform. The test core is presaturated with oil and water, then subjected to gas drive. During the test, required data are the rates of oil, water, and gas production, together with pressure drop and temperature. The ordinary gas-oil unsteady-state relative permeability apparatus can be readily modified to measure the required data. The proposed technique was applied to samples of a Berea and a reservoir core. The effect of saturation history upon relative permeability was studied on one Berea core. It was found that increase in initial water saturation has a similar effect upon three-pbase relative permeability as it does in two-phase flow. INTRODUCTION In the light of increasing demand for three-phase relative permeability data for predicting the performance of thermal and other multiphase-flow recovery processes, a simple and accurate method of experimental determination of such data is extremely desirable. Leverett and Lewis 1 described the simultaneous flow method of obtaining three-phase relative permeability data. However, Caudle et al.2 reported that this method is very time consuming and cumbersome. Corey3 proposed calculating the three-phase relative permeability from measured krg data. Corey&apos;s theory is based on simplified capillary pressure curves,4 assuming a straight line relationship between I/Pc2 and saturation. Also, Corey&apos;s method assumes a preferentially water-wet system. The simplest and quickest method of obtaining three-phase relative permeability data is the unsteady-state method where, for instance, oil and water are displaced by gas. However, in such a test the correlation of average saturation with relative permeability does not give a valid relationship because the rates of oil, water and gas flow in the sample change continuously from the upstream to downstream end. This difficulty in calculating valid relationships was solved by Welge for two-phase flow by deriving an expression for determining the outflow face saturation from Buckley and Leverett frontal advance equations.5,6 In this paper, relations are established to determine the outflow face saturation and relative permeability to all phases in a three-phase flow displacement experiment. PROPOSED METHOD The fundamentals established by Buckley and Leverett5 for two-phase flow were extended to three - phase flow and used as a basis for the derivation of saturation equations. This approach is comparable to Welge&apos;s6 use of Buckley and Leverett theory in arriving at expressions to determine the outflow face saturation of the displacing fluid in a two-phase flow system.

Jan 1, 1967

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.&apos; presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

By N. E. Truitt, D. G. Russell

The transient pressure behavior of a well which produces a single compressible fluid through a singte-plane wrticat fracture has been investigated mathematically. The fracture is assumed to possess infinite flow capacity, to be of limited mdial extent, and&apos; to penetrate the producing formation completely in the vertical direction. Previous studies of vertically fractured wells have been concerned primarily with production rate performance or semisteady-state pressure behavior. This study was undertaken to ascertain the influence of vertical fractures on transient pressure tests such as pressure build-ups and flow tests. In a vertically fractured system, flow in the region nearest the fracture is practically linear, whereas farther away from the fracture essentially radial flow prevails. Thus, transient pressure analyses based on radial flow theory are somewhat inaccurate. As fracture penetration increases radially, kh values calcutated from pressure build-up and flow test curves become increasingly larger than true values. Failure to consider the effect of fracture penetration also introduces inaccuracies into the catculation of fracture length from the apparent skin factor and into the determination of average reservoir pressure. If the total length of the fracture is 20 per cent, or greater, of the drainage radius of the well, corrections must be made to pressure build-up and flow test results. Methods for correcting such results are discussed in this paper. For wells with prefracturing pressure build-up or flow test data, it is possible to estimate fracture length by comparison with postfracturing build-up or flow test results. In new wells or wells without prefracturing build-up or flow test data, fracture length must be estimated to correct the values obtained from analysis of pressure tests after fracturing. Fracturing efficiency calculations should be made whenever possible to provide an estimate of fracture length. Tables of the dimensionless pressure drop as a function of time and fracture penetration are included in this paper. Using these values should permit analysis of other types of transient pressure behavior in vertically fractured wells. INTRODUCTION Hydraulic fracturing has been used quite successfully for over a decade as a completion and stimulation technique in oil and gas wells completed in low-permeability reservoirs During this period a considerable amount of theory has evolved on the performance of hydraulically fractured reservoirs and on more efficient means of artificial fracturing. Although theory has been developed, no rigorous investigation has been made of pressure build-up and flow test behavior in such wells. Prats et al.1 first discussed the performance of vertically fractured reservoirs for the case of a compressible fluid. Their work was primarily concerned with production performance at constant flowing pressure. These authors also considered large-time (semisteady-state) constant production rate behavior for vertically fractured wells; however, transient pressure behavior at constant rate was not investigated. McGuire and Sikora10 and Dyes, Kemp, and Caudle2 employed an electrical analog to investigate the influence of artificial vertical fractures on well productivity and pressure build-up. They found that fractures which extend beyond 15 per cent of the drainage radius away from the well alter the position and slope of the straight-line portion of the build-up curve. They concluded that these effects must be considered both in the determination of the effective permeability of the formation and in any calculations of final build-up pressure. Although these authors did not undertake an exhaustive study of the influence of vertical fractures on pressure build-up performance, their limited results were quite interesting from the standpoint of the effects they demonstrated. In a more recent paper, Scott- reported the results of an investigation of the effect of vertical fractures on pressure behavior, which was conducted with a heat flow model. Scott&apos;s results appear to be consistent with those reported in Refs. 1 and 2. However, the effects of different fracture lengths on performance were not investigated. Pressure build-ups and transient flow tests are among the most diagnostic tools available to the reservoir engineer or production engineer. Since a very high percentage of present-day well completions incorporate the hydraulic fracturing technique, a definite need exists for information on the effect of fractures on transient pressure performance. For these reasons we have undertaken a rigorous study of pressure build-up and flow test behavior in vertically fractured reservoirs. The objectives of this study were (1) to obtain synthetic pressure build-up and flow test

Jan 1, 1965

By K. K. Clark

Excessive injection pressures in water injection wells may create deeply penetrating fractures, or may open up existing reservoir fractures. If these fractures are oriented toward offset producing wells, oil recovery at injected water breakthrough and ultimate recovery will be impaired. A method of calculating fracture lengths from pressure fall-off test data is presented. The method is based on a linear flow model that simulates conditions present during the early-time period after shutting in an injection well. Fracture lengths can be calculated directly if formation permeability is known. A graphical technique is presented that provides fracture lengths based on permeabilities calculated from normal fall-off tests, where those permeabilities are adjusted for flow geometry. Test data from four injection wells and results obtained from applying the method to these wells are discussed. Introduction Rates of water injection into nonstimulated wells in low permeability reservoirs frequently fall below economically desirable levels. Therefore, some form of stimulation such as hydraulic fracturing often is performed on these wells. Usually these fractures penetrate a short distance from the well — short as compared with interwell spacing — and therefore should not create problems of water channeling to offset producing wells. In the past, some operators have been tempted to increase injection pressures continually to maintain a specified injection rate. But pressures cannot be increased indefinitely. Field data indicate that continued injection above fracture opening pressure can cause excessively long fractures to approach between well distances. If these long fractures intersect or closely approach offset producing wells, premature water breakthrough can result. It is difficult to seal effectively these fractures by various workover procedures after premature water breakthrough has occurred. This paper describes testing procedures that are potentially capable of detecting premature water breakthrough before it occurs. Knowledge of fracture lengths and fracture opening pressures can assist field engineers in selecting optimum operating conditions for water injection wells. Fracture lengths should be calculated from pressure fall-off tests on injection wells, and maximum permissible iniection pressures should be determined from step-rate injectivity tests. Equipped with this information, the field engineer is in a better position to recommend for a particular injection well the optimum operating conditions that will lessen the chances of channeling to offset producers. Historically, pressure transient tests have been conducted primarily on producing wells and analysis procedures have been based on radial flow concepts. Several authors have suggested that these same procedures be applied to tests conducted on water injection wells. This application is justified for wells that are not connected to extensive fracture systems. However, if the formation is severely fractured, the interpretation methods must be modified. Several authors have presented the pressure response of various systems that are intended to simulate reservoir conditions. A paper by Russell and Truitt* on vertically fractured svstems is probablv the most useful from the standpoint bf the fieid engineer. Russell and Truitt presented the pressure-vs-time behavior of a well connected to vertical fractures with various fixed lengths. Their permeability adjustment technique is used in this paper to refine methods of calculating fracture length and interwell permeability. Theory and Definitions Three flow-system geometries are considered in analyzing injection well test data. These various geometries are applicable during specific time intervals of transient pressure tests. Schematic drawings of the flow models, and the time intervals over which they are applicable, are shown in Fig. 1 in the chronological order that they are applied to the pressure response of water injection wells in fractured formations. The very early response is simulated by an infinite linear system as shown in Fig. 1A. This representation permits the calculation of matrix surface area exposed to the fracture. The model shown in Fig. 1B is identical with the one studied by Russell and Truitt. They evaluated the transient pressure response of this vertically fractured system by finite difference techniques. The results were presented as tabulations of dimensionless pressure drop functions for a range of fracture lengths. A semilog plot of these data shows that the apparent matrix permeability can be related to the true matrix permeability by a simple exponential function. Their adjustment factor is used in this paper to improve the accuracy of matrix-permeability and fracture-length calculations.

Jan 1, 1969

By K. E. Gray, M. S. Seth

Equations of elasticity for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure were formulated. Using an analogy between thermal and porous body stresses, it was shown that the solution for a transversely isotropic porous body may be obtained by incorporating body forces and the stresses due to a boundary load into the corresponding solution for the thermal stress problem. Equations of elasticity ad the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant pressure or a constant rate of flow is maintained at the wellbore. The vertical and tangential displacements are zero from the conditions of the problem. A numerical analysis was made of the solutions obtained by using a digital computer to determine the relative influence of each system variable. Considering rock as a porous body with internal fluid pressure generally gives results significantly different than considering the rock to be nonporous; the directional character of rocks leads to significant differences as compared to results based upon the common assumption of isotropy. Stress gradients are high near the wellbore but die out away from the well. Radial stresses are compressive or neutral, whereas tangential stresses are tensile, neutral or compressive, depending upon the boundary conditions and physical properties of the system Vertical stresses are compressive for an unbounded system. For constant wellbore injection rate, the vertical stress is proportional to the rate of fluid injection and decreases with time, whereas the radial and tangential stresses increase with time. At a given location, the radial displacement generally is very dependent upon time. INTRODUCTION A realistic appraisal of the state of stress in subsurface rock formations would be of considerable interest and use to the petroleum industry. For example, knowing the state of stress in proximity to a wellbore would be of fundamental importance in designing a fracturing operation or, more important, of clearly understanding the conditions necessary to produce rock failure of desired dimensions and geometry. Understanding conditions necessary for rock failure at the wellbore would also be of utility in a preventive sense. For example, borehole stability is an important consideration for many rock formations, and knowledge of the stress state at and near the wellbore under conditions of substantial pressure gradient due to fluid flow would be of great value. When fluid flows through a porous body which is initially at some uniform stress level, the following forces generate stresses at any point in the body. (1) Forces due to nonuniform pressure distribution. With increasing pressure the elements of a body are compressed. Such compression cannot proceed freely in a continuum when the pressure is not uniform throughout, and thus, stresses due to flow of fluid are set up. (2) Pore fluid pressure. This gives rise to normal stresses whose value at any point is the product of areal porosity and the fluid pressure. Although the fluid exerts uniform pressure, the stresses it creates in an anisotropic body may not be the same in all directions since the areal porosity in an anisotropic porous body is a direction-dependent quantity. This consideration leads to the concept of directional porosity. (3) Body Forces. The flow of the moving fluid imparts a force to the rock matrix, the magnitude of which depends upon the pressure

Jan 1, 1969

By K. E. Gray, M. S. Seth

In Part I of this work,1 equations of elasticity were formulated for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure. Equations of elasticity and the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant rate of pressure or a constant rate of flow is maintained at the wellbore. In this paper, the approach presented earlier is extended to finite reservoirs for the cases of (1) steady-state flow, (2) constant pressures at the well bore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Results of this work show that radial and tangential stress gradients are high near the wellbore but diminish rapidt; away from the well; the vertical stress gradient behaves in the same way but is less severe. Radial stresses are compressive or neutral, whereas tangential and vertical stresses may be tensile, neutral or com-pressive, depending upon the boundary conditions, the physical Properties of the system and the radial distance involved (vertical stresses are always compressive in an unbounded system1. For constant boundary pressures, both radial and tangential stresses increase with time whereas they both decrease for a closed outer boundary and constant pressure at the wellbore. The vertical stress decreases with time for both systems. For steady-state systems, radial displacement may be positive or negative, depending upon the dimensions of the system, the pressure differential and the porosity. Radial displacement may be positive or negative for a closed outer boundary but is positive for constant pressures at both boundaries. INTRODUCTION The importance, utility and complexity of a realistic appraisal of the stress state at and local to a wellbore were indicated in Part I. In this paper the analytical approach presented earlier is extended to finite, cylindrical reservoir geometry for the cases of (1) steady-state flow, (2) constant pressures at wellbore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Other than the outer boundary of the reservoir being finite, the physical system and assumptions pertinent thereto are the same as before. The reader may wish to review the mathematical development through Eq. 49 of Part I before proceeding here. SOLUTION FOR GENERAL BOUNDARY CONDITIONS THERMAL STRESSES AND DISPLACEMENTS Since normal stresses do not exist at the free boundaries of a system subjected to thermal stresses, the boundary conditions are

Jan 1, 1969