Search Documents

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By F. G. Miller, T. D. Mueller, L. B. Davidson

A mathematical model of reservoir behavior during one cycle of the intermittent steam injection process is presented in this paper. A cycle of the process is considered to be composed of three stages: injection, condensation and production. The model concentrates on behavior during the production stage. The influence of several parameters on production improvement is described, using the results obtained with a numerical equivalent of the model. Reservoir system parameters considered are thickness, reservoir thermal conductivity, oil viscosity, solution gas-oil ratio and the change in oil effective permeability with temperature. Operating parameters investigated are heat injection rate, steam temperature and cumulative injection time. In the face of considerable uncertainty, the model appears to give a fair prediction of the behavior observed during an actual test. This makes the model a useful means of selecting projects. It is shown that the model may also be utilized to select an optimum set of operating conditions for a particular project. Results of this study show the extent to which production improvement depends on each of the parameters investigated and give an indication of the range of system conditions within which the process would be effective. INTRODUCTION In recent years the intermittent steam injection or steam soak process of oil recovery has become a popular method of improving recovery in reservoirs bearing high-viscosity hydrocarbons. Production improvement due to the process results from heating the reservoir near the wellbore. The process is applied on an individual well basis and is executed in a series of cycles, each consisting of three stages: injection, condensation (the soaking period) and production. During the injection stage, high-pressure steam is injected into the reservoir surrounding the well-bore. This study assumes that injection is carried out over the entire vertical section of the reservoir rock that can be exposed to fluid injection. The condition of steam entering the reservoir is considered to fall within the following limits: (1) steam pressures less than 1,000 psig, (2) steam temperatures between 300 and 700F and (3) steam quality less than 90 percent. The last restriction is due to equipment limitations. Normally, injection re maintained for less than 3 weeks. After termination of injection, the condensation stage begins. The well is shut in for several days, supposedly to make certain that the injected steam condenses. This condensation would insure the retention of injected heat. The production stage is initiated when the well is opened again. Reservoir response during this stage determines the effectiveness of the process (response is defined as the oil production rates, the water production rates and the wellbore temperatures observed during the entire production stage). The immediate response of the reservoir is sometimes quite good. Sets of case histories by Longl and by Owens and Suter2 give some typical figures for response. Duration of the production stage is equal to the time needed for the oil production rate to decline to some predetermined value. This usually requires 4 to 6 months. Termination of the production stage coincides with the start of the next cycle. The process is continued, cycle by cycle, until it becomes unprofitable to continue. One simple model2 and two more sophisticated models3,4 of the process have been presented. These are based on fundamentally different assumptions than those used in this investigation. Although the three models might be modified to fit the assumptions used in this study, it is felt that greater rigor would result by developing a

By C. Wachmann

The theory presented ill this paper makes possible examination of the phenomenological aspect of secondary oil recovery by alcohol flooding. Limitutions of the theory are contingent on three primary assumptions: (1) chemical equilibrium, (2) no fingering, and (5) no dispersion. The second assumption implies one-dimensional flow. In a geometric sense the theory extends the Buckley-Leverett construction to systems of two slightly miscible fluids and to systems of three components. This is made possible by considering single-phase fluids to be made up two phases, namely of an excess of one phase and a negative amount of the other phase. The procedure is discussed in the paper. The theory also permits direct determination of the relationship between the aqueous phase saturation and the concentration of water in the aqueous phase in a zone of variable composition. The velocity of planes of constant water concentration and of constant saturation can be computted likewise and the existence of saturation discontinuities inferred. Consequently the behavior of such a Zone, which may exist at some time between a displaced and a displacing fluid, May be determined. mathematical details are given in Appendices. Problems of experimental verification are pointed out on the basis of a supporting experiment. Implications of the theory and effects of violations of primary assumptions are discussed briefly. The major conclusion is that, if the primary assumptions are met, with pure alcohol and alcohol-rich .slugs piston-like displacement of oil and water is possible. INTRODUCTION The recovery of oil from porous media by alcohol displacement has been studied for several years. Extensive experimental investigations into alcohol flooding have been carried out in the laboratory although little has been done in the field. A large amount of practical information is available, the extent being indicated by the first three papers cited 1-3 and their respecrive bibliographies. This has permitted the Formulation of inferential theories. A more analytical approach in which theories are developed primarily by reasoning from known physical principles and relationships has, however, lagged behind. In fact, an analytically derived theory of alcohol displacement has not been available. The theory presented in this paper on oil and water displacement by alcohol from a porous medium is a result of the analytical approach. The theory is idealized, as certain assumptions had to be made to make the mathematics tractable and reduce complexity of argument. In particular, major assuptions concerning conditions of equilibrium, interfacial stability, longitudinal mixing, etc., led to the exclusion of fingering, for example, from all considerations. In many laboratory experiments reported in the literature fingering could not be avoided. There is, therefore, disagreement between theories based on experimental evidence available and the analytically derived theory presented here. This theory presented is based on known physical data and makes use of the Buckley-Leverett construction and its suggested extension. The mathematical derivation of the theory is given in Appendices. Three examples relevant to alcohol displacement are discussed to explain the use of this theory: (1) system of two immiscible fluids (the Buckley-Leverett case); (2) system of two slightly miscible fluids involving an extension of the Buckley-Leverett construction; and (3) system of three components, including injection of pure alcohol or alcohol-rich fluid as a special case, which involves a more general extension of the Buckley-Leverett construction. In these discussions strict assumptions were made to impose equilibrium conditions and exclude mixing and fingering. Effects of relaxing these assumptions and consequent extensions of the theory are indicated. The complexity of the alcohol displacement process makes experimental verification of the theory essential. Some relevant evidence is offered in this paper.

Jan 1, 1965

By B. Zondek, W. T. Cardwell, J. W. Sheldon

The fundamental equations that are used to describe two-phase fluid flow in porous media are Darcy's law for each phase and an equation of continuity for each component. The special case of one-dimensional, incompressible, two-phase flow has received much attention in the petroleum engineering literature.1-10 The basic paper on the subject is that of Buckley and Leverett.2 Buckley and Leverett showed that one could eliminate the pressure and obtain a single partial differential equation for the saturation. They solved this equation for the case when gravity and capillary forces are negligible. A frequently encountered property of the solution of the basic partial differential equation is that, as time progresses, the saturation becomes a multiple-valued function of the distance coordinate, x. Buckley and Leverett interpreted the formation of multiple values as an indication that the saturation-distance curve had be-come discontinuous. They developed a method of determining the position of the discontinuity from a material balance. We solve the Buckley-Leverett partial differential equation using the method of characteristics and the concept of shocks. Our treatment is analogous to treatments used in the theory of supersonic compressible flow." Although the numerical results obtained are the same as those obtained by previous methods, we believe that this approach is more logical and has broader applicability. It can, for example, be generalized to treat the case when the fluids are compressible. It is not clear how one would go about doing this using the method of Buckley and Leverett. To illustrate the power of the method of characteristics with shocks, we solve a problem in gravity segregation which has not been treated previously in the literature. A special case of the general method developed here was used by Welge, to solve a particular fluid displacement problem. We consider two-phase incompressible flow along the vertical direction x in a porous medium. We assume that there is no flow transverse to .x. We measure positive in the upward direction. The fundamental equations are obtained from Darcy's law and the mass conservation law. Since the two phases are incompressible, mass conservation is equivalent to volume conservation.

By C. S. Matthews, P. Hazebroek, F. Brons

A method has been developed for calculating the average pressure in a bounded reservoir. The reservoir is first divided into the individual drainage volumes of each well, by using the criterion that at steady state each individual drainage volume is proportional to a well's production rate. The average pressure in each drainage volume is then calculated by a method developed in the report. By volumetrically averaging these individual drainage volume pressures, the average pressure in the entire reservoir is obtained. To calculate the average pressure in each drainage volume, a correction is applied to the ordinary extrapolated pressure, i.e., the pressure obtained by extrapolating to infinite time the linear portion of the graph of closed-in pressure versus log[?t/(t + ?t)], where At is the closed-in time and t the production time. The correction, which is a function of the production time, is presented in graphical form for different shapes of the drainage area (horizontal cross section of the drainage volume). INTRODUCTION It is important to be able to find the volumetric average pressure in a reservoir so that the size of the reservoir may be determined from material balance calculations. It is also desirable to be able to find the approximate distribution of pressure within a reservoir for detection of fluid movement. 'The purpose of this paper is to present a method for calculating both the average reservoir pressure and the approximate distribution of pressure within a bounded reservoir — that is, a reservoir with no water drive. In reservoirs where the pressure builds up rapidly after wells are shut in, the determination of average pressure generally poses little problem, for one often need only average the final buildup pressures. It is when pressure buildup is slow that difficulties arise. For practical and economical reasons, the time allowable for closing in wells is limited. If at the maximum allowable closed-in time the pressure has not reached a con- stant value (and this is more often the case than is generally realized), calculation of average pressure presents difficulties. One method which has been used in the past for obtaining an "average" pressure involves the extrapolation to infinite closed-in time of curves of pressure as a function of closed-in time. Contour maps of these extrapolated pressures are used in arriving at an average pressure. In some cases wells are shut in for arbitrary Lengths of time, and contours of these arbitrary shut-in pressures are used for calculation. All these methods are somewhat arbitrary and suffer from lack of fundamental theoretical basis. For example, no precise meaning can be given either to the extrapolated or to the arbitrary shut-in pressure. Furthermore, there seems to be no reason why either set of pressures should be contourable. It was because of the vagueness of these old methods and their lack of a sound theoretical basis that, the present method for determining average reservoir pressure was developed. SYMBOLS A = drainage area of a well, cm' A, = drainage area of well j, cm2 A , = total drainage area of all wells in a reservoir, cm2 c = coefficient of compressibility of fluid at reservoir temperature, atm-1 f = hydrocarbin-filled porosity, fraction h = net thickness of formation, cm k = permeability of formation, darcy p = pressure at a well, atm p0 = initial reservoir pressure, atm p* = pressure obtained by straight-line extrapolation of linear portion of plot of p versus log [At/(t+ At)] p = volumetric average pressure inside drainage volume of a well q = volume flow rate at a well, cc/sec a: prevailing reservoir conditions q, = total volumetric flow rate from a reservoir, cc/sec at prevailing reservoir conditions rw = radius of well bore, cm T = kt/fµcA t = corrected time of production of a well, seconds, defined by

Jan 1, 1955

By R. H. Jacoby, V. J. Berry

When dry gas is injected into a reservoir containing a volatile crude oil, a significant amount of the reservoir liquid phase may become vaporized. The resultant rich gas phase, when subsequently produced, con-tributes to tank oil production. This contribution assumes greater importance, the more volatile the oil in-olved. Oil recovery may be sub-stantially greater than that predicted by conventional frontal-drive methods, which do not consider the vaporization equilibrium between the reservoir oil phase and the injected gas. A calculation method has been (developed to account for vaporiza-. tion of the reservoir liquid phase during gas-injectiorz operations, and for the tank oil production which results from this factor. Recovery performance calculations are presented for a reservoir containing a highly volatile oil. Tank oil recovery is calculated to be about twice that predicted by the use of the conven-rional frontal-drive equations. In contrastto usual pressure maintenance performance results, in which the gas-oil ratio rises at an increasing rate after gas breakthrough, the pre-dicated gas-oil ratio rises rapidly to about 12,000 scf/bhl and then rises much less rapidly. During gas inje-tion, most of the reservoir liquid phase contacted is evaporated by the dry injection gar. The gas-oil ratio during this period is dependent upon reservoir pressure. The higher the operating pressure, the lower the gas-oil ratio. The predicted behavior i., in accordance with laboratory PVT tests made to .simulate the vaporization behavior. In addition to recovery performance predictions, results of the calculation procedure include complete wellstrearn composition data of value in the design of gasoline plant facilities often used in con-rrction with gas-injection operations. INTRODUCTION In the cycling of gas-condensate reservoirs, dry gas is injected to maintain reservoir pressure during wet gas production and to thereby eliminate or reduce ultimate loss of liquids due to retrograde condensation within the reservoir. Gas injected into crude oil reservoirs has a dual function. It displaces oil to the producing wells and at the same time serves to partially or fully maintain reservoir pressure. Oil shrinkage which would occur upon pressure reduction is thereby minimized or eliminated. Accepted calculation methods are available'= for predicting recovery performance of either gas-condensate reservoirs or crude oil reservoirs which are being subjected to gas injection. In gas-condensate reservoirs, any retrograde liquid formed does not flow, and it is necessary to account only for the vaporization equilibrium between this liquid and the injected gas. Conversely, in normal crude oil reservoirs, both the oil and gas phases flow, but it has not been considered necessary to account for any vaporization of the reservoir liquid which might occur upon contact with the dry injection gas. Recently, high shrinkage reservoir fluids known as "volatile oils" have been found in increasing amounts. These oils are characterized by tank oil gravities above 4.5" API, solution gas-oil ratios above 1,000 scf/bbl, and reservoir volume factors above two. Special techniques have been devised for predicting depletion performance of reservoirs containing such oils.74.5.1; One of the characteristics of reservoirs producing volatile oils is that the reservoir gas phase carries a significant amount of oil which is recoverable as stock-tank liquid. This unusual vaporization be-havior implies that an appreciable amount of reservoir liquid would be vaporized upon contact of the oil with dry injection gas. In a gas-injection operation, tank oil recovery would he obtained not only through frontal displacement of the reservoir liquid by the injected gas, but also through production of the rich gas phase. This means that not only are improved methods needed to predict recovery of such oils from reservoirs undergoing gas injection, but it would also be expected that high oil recoveries might be obtained by such operations. When relatively dry gas is injected in to a volatile oil reservoir, phase equilibrium between the injected gas and the reservoir oil will tend to bc established. Initially the most volatile components, such as propane, the butanes and the pentanes, will account for most of the material transferred from the oil phase to the gas phase. As the partially stripped oil phase is contacted with additional dry injection gas, the heavier intermediate components, such as the hexanes, the heptanes and the octanes, will gradually transfer to the gas phase in increasing amounts. This is because the supply of lighter components in the oil phase dwindles due to the stripping action of the injected gas. This stripping action will usually continue to be effective down

By B. H. Caudle, L. G. Davies, I. H. Silberberg

This paper presents a method of predicting the recovery and performance of a five-spot steam injection project, in which a realistic approach to pattern sweepout efficiencies is made. Published methods for radial systems were modified for the five-spot pattern by approximating the stream lines with straight lines radiating from the injection well and then converging to the producing well. In each radial segment, the position of the steam front and the temperature profile ahead of the steam front were determined by heat balance equations, which included an estimation of heat losses to surrounding formations. The location of the saturations behind the cold water front was determined from a Buckley-Leverett solution to the material balance equation. Results from this program show steamflood recovery in a five-spot pattern to be considerably less than that predicted for true linear or radial flow systems. For a specific reservoir containing 900-cp oil, a steamflood in a purely radial flow system was predicted to recover more than 75 percent of the original oil in place when 2 PV of water had been injected as steam. A five-spot steamflood with otherwise identical properties was predicted to recover 10 percent of the original oil in place when 0.15 PV of water had been injected as steam and to recover almost no oil thereafter. A cold water five-spot flood in this system was predicted to recover approximately W percent of the oil in place with I PV of water injected. For a five-spot pattern in an example reservoir with 10-cp oil, steam injection similarly showed lower ultimate recovery than water injection but no improvement in recovery rule. Introduction The thermal recovery method considered in this study is steam injection in a five-spot pattern. Pattern steam injection has been given little attention in the past, possibly because of some rather obvious disadvantages. Unless steam temperature is maintained throughout the entire swept area, the process will revert to simple waterflood with all heat being lost prior to reaching the oil bank. Also, the oil viscosity reduction takes place near the steam front and not around the producing wellbore, so that low producing rates must still be endured. This is the reason for the success of the steam stimulation, or cyclic injection, in which the fluids are produced from the same well used for injection. On the other hand, steamflooding can offer some advantages. Oil displacement is by four methods: (1) mechanical displacement by the condensed water, (2) viscosity reduction of the crude oil, (3) swelling of the crude oil, and (4) distillation of the crude oil in the steam zone. Although dependent upon the crude, laboratory experiments have shown displacements up to 80 percent by steamflooding. In addition, steam is a good heat transport medium, since it is cheap and has a high heat content. Previous Investigations Previous publications have reported methods for predicting recovery in a steamflood for linear and radial systems. Since no pattern sweep efficiency is taken into consideration, the recovery even from a radial system must be greater than from more realistic geometries such as the very common five-spot patterns. Probably the broadest coverage is given by Willman who reported experimental results and offered a method of predicting recovery for a radial flow system. They concluded that both hot water and steam injection recover more oil than an ordinary waterflood, and that steam injection could yield recoveries "as much as 100 percent greater than by water flood." They also concluded that both the heat requirements for a reservoir and the residual oil remaining in the reservoir after steam injection were independent of the amount of oil originally in place, that short exploitation times were desirable, and that a high percentage of net sand in the reservoir with a high initial oil saturation was desirable. The method assumes that the flood occurs in three concentric cylindrical zones: (1) an inner steam zone, (2) a central hot water zone, and (3) an outer cold waterflood zone. Displacement in the steam zone is based on laboratory-determined residual oil saturations while the hot water and cold water zones use the conventional Buckley-Leverett equations. Although the results shown by Willman et al. are for a radial system flowing out from the well to an assumed external circular boundary, their equations did

Jan 1, 1969

By J. S. Aronofsky

The first step in a quantitative analysis ot the mechanism of oil displacement by water in a fractured reservoir is usually conceded to be the solution of the differential equation describing the saturation distribution of two immiscible fluids flowing in a porous medium, where the capillary pressure is taken into account. In such a system the production mechanism may consist of displacement of oil both by the flow of water due to natural or artificially imposed pressure gradients and by imbibition, which implies a flow of water not due to external pressure gradients. Owing to the presence of the two oil displacement mechanisms, the mathematical model given by the differential equation intended to describe the system may not properly represent the behavior of the physical system. In fact, in the reservoir the rate of water advance may be very slow, and in the case of a fractured reservoir with a great number of large fractures, the pressure difference determining the flow of water through the matrix may be much less than 1 lb/'psi over lengths of a few feet. In such a case, imbi-bition (the exchange between oil in the matrix and water in the fractures resulting from capillary forces) may become, with time. a significant element of the production mechanism. It occurred Lo the authors, however, that without going into a physical analysis of the process of production, it might be possible by means of simple abstract reasoning to throw some light on the variation of recovery with time under conditions occurring in a highly fractured oil reservoir with rising water table. The object of this paper is to present both the reasoning and its application to a reservoir of the highly fractured type. Specifically, the analysis given here was undertaken to try to explain the increase of recovery (as defined later) with time as observed in this reservoir, without having to assume unlikely variations in the reservoir parameters with depth. This attempt has been successful as will become clear upon comparison of the computed recoveries with the actual field data. ABSTRACT MODEL Let us consider a small volume. of porous matrix saturated with oil at time, t = 0. Let the process of oil displacement by water start at time, t = t,,. At some time, t, the process will have terminated. Then a volume of oil equal to or smaller than the original oil contained in the matrix will have been produced. The first basic assumption that describes the model and guides the forthcoming reasoning is that the oil production from the small volume. dv, is a continuous monotonic function of time and that it converges to a finite limit. Such an assumption is not inconsistent with the results of laboratory waterflood tests as well as with results of imbibition tests where this is, in fact, observed. The second basic assumption is that none of the properties which determine the rate of convergence change sufficiently during the process to affect this rate or the limit. Let it be assumed that the form of the function of time relative to production from the matrix volume. dv, is given by where V,(t) is the volume of oil produced up to that time t. R is the limit toward which the recovery converges, A is a constant giving the rate of convergence, and V,(t) is the volume of oil originally in place in the volume, dv. It should be noted that recovery at time I will be understood here to he It follows that r. tends to R as t tends to infinity. CONSTRUCTION OF RESERVOIR FROM ABSTRACT MODEL Let the reservoir consist of a series of identical blocks of porous matrix stacked vertically and separated by fractures. Let each of these blocks satisfy the conditions of our abstract model. These conditions are: (1) that recovery is a continuous mono-tonic function of time converging to finite limit, and (2) none of the properties that determine the rate of convergence change sufficiently to affect the rate or the limit. Let water be rising in fractures so that oil production from any part of the block starts when the water comes in conLact with it. TOTAL RECOVERY FROM THE RESERVOIR COMPARED TO RECOVERY FROM A SINGLE MODEL ELEMENT As stipulated in Eq. 1, in the case of the abstract model,

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

A numerical method for computing the dynamical behavior of fluid- fluid interfaces is described. Results of studies made to assess the accuracy and economy of the method on a computer are reported. It is concluded that it is practical to obtain digital computer solutions for fluid- fluid interface problems with two space dimensions. PROBLEM DEFINITION We consider two-dimensional flow in a rectangle as shown in Fig. 1, where a fluid-fluid interface divides the rectangle into two regions. The pressure is specified along the lower boundary (the injection boundary) and upper boundary (the withdrawal boundary). There is to be no flow of fluid across the side boundaries, which are walls of symmetry for the flow, and px = 0 along these boundaries. The fluids are incompressible, the absolute permeability and porosity is constant, and each fluid has its own constant mobility. Thus, the pressure satisfies Laplace's equation in each region. Across the interface, pressure is continuous and the normal component of fluid flux is 1. These conditions are shown in Fig. 1. The velocity of a point on the interface in the direction of the normal from Fluid 1 into Fluid 2 is U,/ Q, where Un is the normal component of volumetric flux at the point on the interface, i.e., For brevity, we call the problem defined above "Problem A". Problem A is a simple example of a dynamical fluid-fluid interface problem. However, the same numerical methods which we use to solve Problem A can be extended to the solution of a large number of fluid-fluid interface problems which are mathematical models for various oil recovery processes,l hence our interest in Problem A. We define the interface by prescribing the x, y coordinates of a finite set of points. These points are called "interface mesh points" and are illustrated by triangles or asterisks in the figures. Fig. 2 is a procedural (flow) chart displaying the main steps in the solution. We comment briefly on some of these steps in this section and then give further details in subsequent sections. For mobility ratio greater than one it is well known that Problem A is unstable.*2 Small amplitude,

Jan 1, 1965

By R. H. Jacoby, V. J. Berry, R. C. Koeller

The experimental phase behavior of several field gas-condensate systems, one field volatile oil system, and a series of synthetic systems having gas-oil ratios from 2,000 to 20,000 scf/bbl stock tank oil was measured. The data were used to calculate the depletion performance of the field systems at their respective reservoir temperatures and of the synthetic systems at various temperatures. The results of the performance calculations were used to prepare correlations of total stock tank oil and separator gas in-place, and the amount of each recovered by primary depletion. Data needed to use the correlations are initial gas-oil ratio, initial tank oil gravity, reservoir temperature and reservoir pressure. The correlations presented should be useful for estiinating recovery from reservoirs producing volatile oils or rich gas condensates having gas-oil ratios from about 2,000 to 30,000 scf/bbl of stock tank oil. INTRODUCTION Prediction of depletion performance for gas-condensate reservoirs usually requires extensive laboratory work or lengthy phase-behavior calculations. For many small gas-condensate bearing zones the expense of such work is often not justified. Furthermore, where large holdings are involved, it is useful to have estimates of future oil and gas production in advance of the results of a laboratory analysis, even if such work is ultimately planned. For these reasons, it would be of value to have a method for rapidly estimating the depletion performance of gas-condensate reservoirs from the type of field data usually available. Few data on the reservoir depletion performance are available for the purpose of developing correlations. The increasing number of gas-condensate discoveries and the advent of high-speed computers have stimulated the development of methods for predicting performance. Such calculated performance data have been shown to be reliable and are becoming accepted as a basis for reserve estimates and for engineering studies. It therefore seems satisfactory to use results of performance predictions themselves for the purpose of developing the desired correlations. Such an approach is described in this paper. Several rich natural gas-condensate systems (gas-oil ratios of from about 3,600 to 60,000 scf/bbl) and one natural volatile oil system were studied (2,363 scf/bbl). In addition, to provide a more systematic basis for correlating purposes, a series of related, synthetic reservoir fluid mixtures was also studied (GOR's from about 2,000 to 25,000 scf/bbl). The purpose of this paper is to show how the recovery of oil and gas by primary depletion may be estimated from surface gas-oil ratio, reservoir temperature and pressure, and tank oil gravity. procedure Experiments Companion samples of gas and oil were collected from the high pressure separators in nine different fields at the time of their discovery. These samples were each recombined in accordance with the respective producing GOR to reproduce the reservoir fluids. Eight of these fluids, designated F-2 through F-9, were rich gas-condensate systems at the respective reser-

By M. R. Tek, F. F. Craig, J. O. Wilkes, C. S. Goddin

The waterflood performance of a water-wet, stratified system with crossflow is computed by a finite difference procedure. The effects of five dimensionless parameters on tile oil displacement efficiency, water saturation con-tour.7 and crossflow rates are evaluated in the absence of gravity forces. Crossflow due to viscous and capillary forces is shown to exert a significant effect on oil recovery in a field-scale model of a two-layered. water-wet sandstone reservoir. The crossflow is at a maximum in the vicinity of the front advancing in the more permeable layer. Under favorable mobility ratio conditions, the comparted oil recovery with crossflow always is interrnerliate between that predicted for a uniform reservoir and that for a layered reservoir with no crossflow. INTRODUCTION The important erects of reservoir heterogeneity on waterflood performance are commanding increased attention in the technical literature. Much of this attention is centered on two categories of layered reservoirs: those in which layers are non-communicating and those in which crossflow of fluids occurs between the layers. In the first category, the reservoir is assumed to consist of discrete layers, each uniform within itself and differing from the others only in such properties as thickness, porosity and absolute permeability. The performance within each layer is calculated by one-dimensional flow theory, and the performance of the total reservoir is obtained by summing individual layer performances. Capillary and gravity effects usually are not considered. Representative publications dealing with thi5 type of reservoir are those of Stiles,' Dykstra and Parsons,' Hiatt,3 Warren and Cosgrove' and Higgins and Leighton." Prediction of performance for reservoirs in the second category is considerably more difficult since viscous, capillary and gravitational forces all play important roles in causing crossflow between layers. A number of authors have investigated the simpler problem of two-dimensional displacement flow in a stratified system with a mobility ratio of unity and negligible capillary and gravity effects.'; Others have considered two-dimensional, non-steady-state flow of a single, slightly compressible fluid in a stratified reservoir. A limited number of laboratory oil displacement tests in layered models with crossflow have been reported. Miscible floods (with resultant zero capillary forces) in layered five-spot models were conducted by Dyes and Braun," who studied the effect of mobility ratio with zero gravity forces, and by Craig et al. 12 who studied the effect of gravity forces at constant mobility ratio. Waterfloods in layered five-spot models (with cross-tlow due to capillary, viscous and gravity forces) were conducted by Gaucher and Lindley,"' who showed the effect of gravity forces in causing underrunning of the injected water and by Carpenter, Bail and Bobek, 14 who demonstrated the reliability of Rapoport's" dimension-less parameters for scaling layered systems. Waterfloods in rectangular layered models were conducted by Richardson and Perkins."' who investigated the effect of velocity at constant mobility ratio and with zero gravity forces, and by Hutchinson," who studied the effects of varying mobility, layer permeability and layer thickness ratios. The differential equations which rigorously describe waterflooding in a heterogenous porous medium are non-linear and do not facilitate analytical solution. By using finite difference approximations it is possible to obtain a solution to any desired degree of accuracy. Such a solution, using an alternating direction implicit procedure (ADJP), is described by Douglas, Peaceman and Rachford.18 In the present study, a computer program using ADJP explores systematically the effects of important parameters on waterflood performance of a two-dimensional, two-layered, field-scale model of a water-wet sandstone system. Particular attention is given to evaluation of the water saturation contours and crossflow rates at the interface between layers to gain improved understanding of the crossflow mechanism. PROCEDURE BASIC: FLOW EQUATIONS The basic flow equations for two-dimensional, two-phase, immiscible, incon~pressible flow in a porous medium are:

Jan 1, 1967

By C. Perez-Rosales

A simplified method for evaluating specific surface of porous media is described. The theory takes as its point of departure a previous statistical method. Due to its simplicity, it is considered that the proposed method has certain advantages over methods currently used. By making some simplifications, the procedure for estimating specific surface is reduced to superposing, on an enlarged reproduction of a section of the sample to be analyzed, an evenly spaced grid and counting the number of intersections of the grid lines with the perimeter of pores. An indirect test of the theory is presented and the advantages of the method are discussed. INTRODUCTION The specific surface of a porous medium, defined as the surface area of the pores per unit bulk volume, is an important parameter with regard to flow of fluids through porous materials. In particular, it has long been recognized that specific surface is closely related to permeability. Since the internal surface of any porous material is extremely complex, the measurement of specific surface can only be accomplished by indirect means. It is considered that the best method' in current use for determining specific surface is the statistical method of Chalkley, Cornfield and Parka who originally applied it for estimating volume-surface ratios of morphologic components, such as cells or nuclei. This method will be outlined later. A method which has been used extensively in determining specific surface is based on the adsorption of a vapor by the internal surface of porous media. However, this method does not yield the surface pertinent to fluid flow, since it includes the tiny molecular interstices which do not influence the permeability of the porous material.' This method has been employed by Brooks and Purcell measure internal surface area of sedimentary rocks. The Kozeny equation, which relates the specific surface to the permeability of a porous medium, has found wide use in determination of specific surface. However, the KOzeny equation is not strictly correct. Large discrepancies have been found between theory and experiments.' Consequently. the specific surface estimated by this method is subject to serious uncertainties. Finally, there are other methods for specific surface measurements which are based upon a variety of physical and chemical properties: the heat conduction by a gas, ionic adsorption phenomena and direct chemical reactions of the surface with other substances. This paper describes a simplified statistical method for determining specific surface of porous materials. It is considered that, due to its simplicity, it has certain advantages over the methods presently used. For its validity, it is required that the porous structure of the material under study be of a spatially random nature. This is not a serious limitation since the naturally occurring porous media have this characteristic. THEORY The method proposed in this paper takes as its point of departure the interesting result obtained by Chalkley et al The method shows that if a line of length r is dropped repeatedly at random over a section of a porous medium, and counts are made of the number of times the two end points fall within pore areas (denoted by h for "hits") and the number of times the line intersects the perimeter of pores (denoted by c for "cuts"), then in a very large number of throws where S is the internal surface of the sample (surface of pores) and V, the volume of the pores (Fig. 1). Since the specific surface S, and the porosity 4 are given, respectively, by S. = Sl V, and 4 = V,l V,, where V, is the total volume of the sample, it follows from Eq. 1

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By J. C. Griffiths, D. W. Bennion

A mathematical model, which does not assume a priori that stratification exists, but was designed to test for the stratification was developed. The model segmented the reservoir horizontally into areas of common variance, then divided it vertically into strata (if strata were present). Next, trend surface techniques were used to determine the lateral extent and variation of each stratum. The model was tested on two reservoirs (sandstone and limestone) using porosity as the reservoir rock property. The sandstone reservoir contained approximately 60,000 samples from 2,000 cored wells, while the limestone reservoir had approximately 24,000 samples from 430 cored wells. Within the areas of common variance tested, the model was able to distinguish four separate lithological units or zones in the sandstone reservoir, four in the Marly section and seven in the Vuggy section of the limestone reservoir. INTRODUCTION To accurately predict the movement and production of fluids from a reservoir rock, it is necessary that the reservoir engineer have a knowledge of the rock properties and their distribution throughout the reservoir. For this study, a reservoir rock is defined as a solid containing interconnecting holes or voids which occur relatively frequently and are dispersed within the rock in a regular or random manner. The fraction of bulk volume that these interconnecting voids occupy is called porosity. The rock must also possess sufficient permeability to allow fluids to move through it and it may be composed of one or more strata or lithological units. A stratum or lithological unit as defined in this study is a body or volume of reservoir rock whose properties are so distributed that its lateral and vertical extent can be traced throughout or through a portion of the reservoir. Operationally, the existence of strata are question- able unless the variance within strata is less than among strata. Within the past few years computers have been widely used as a tool for calculating movement and production of fluids from reservoirs. Most prediction models require some type of a mathematical description of certain reservoir rock properties. In the past, most of these models have been relatively simple. Some have assumed vertical variation but no lateral variation, while others have assumed lateral variation but no vertical variation. In general, most of these models have not used functional relationships to predict reservoir rock properties as functions of position. The object of this study was to develop a model which would predict lateral and vertical variations of reservoir rock properties. Since samples are usually available from only a small portion of the total reservoir rock, it seemed logical that if measurements from these samples were to be used to infer the properties of the actual reservoir, the data should be treated statistically as a sample from the total population (reservoir). Therefore, the model developed in this study to predict reservoir rock properties was a stochastic model. The model was tested using porosity (a macroscopic reservoir rock property) measurements from two reservoirs — the first, a sandstone reservoir from which 60,000 samples were obtained; and the second, a limestone reservoir with 24,000 samples. Several investigators have proposed methods to determine the vertical and areal variation of a reservoir rock property. Stiles: Dykstra and Parsons 2 and Suder and Calhoun3 have all developed waterflooding prediction techniques which assume vertical variation but no lateral variation. Stiles developed a method to segment a reservoir arbitrarily using frequency distribution. Law4 suggested that porosity has a normal frequency distribution and that permeability has a log-normal frequency distribution. Elkins and Skov 5 proposed that when a reservoir is segmented vertically, the

Jan 1, 1967

By C. S. Matthews, G. L. Stegemeier

In one field in South Texas, approximately 72 per cent of the pressure build-up results show a characteris-i.rtic "hump" (i.e., the pressure builds up and then falls off) which makes interpretation by. standard methods impossible Correlation of size and time of hump with formation permeability. well productivitv index, and method of completion led to the tentativc, conclusion that the humps were caused by segregation of gas and oil in the wellbore after closing-in"'. This conclusion was confirmed by performance of simple laboratory bubble-rise experiments, by theoretical bubble-rise time calculations, and by a detailed calculation of PVT behavior in the well of a particular well Oil which accurate Surface and bottom-hole pressure measure-ments were untile. The hump behavior has since been found to occur in many other fields. The cause, however, is not the mine in all cases. In some of these the hump is trace-ahie to leaks in the tnbing which allow influx from the annulus after closing-in. In other cases the hump is traceable to leaks in the device separating pay horizons in dually completed wells. It is concluded that the recording of both surface and bottom-hole pressures is desirable in wells which show an anomalous build-up behavior. A number of field examples is discussed where use of both sets of measurements enables the cause for anomalous behavior to he found, (and a reasonable interpreation of bottom-hole pressure to be made. INTRODUCTION Theoretically the pressure build-up in an infinite reser-voir should be a linear function of In [(t + at)/nere t is the production time and At a the closed-in time. Some of the variations from this behavior arc well known, such as the curved portion immediately after shut-in which results from after-production and skin effect, and the flattened end portion which results from boundary effects in a limited reservoir. The effect of stratified producing zones and irregular geometrical drainage patterns may also contribute unusual characteristics to build-ups. However, the effect of still another phenomenon—that of movement of fluids within the wellbore—has been neglected in most build-up studies to date. Fluid movements within the wellbore after shut-in occur as a result of after-production, packer failures, leaks in the casing or tubing, or buoyancy of the gas phase when both gas and liquid are present in the well. Each of these movements can influence the pressure build-up, sometimes sufficiently to negate use of the data for computing permeability or stabilized build-up pressure. Examples of each of these phenomena are discussed below. FIELD OBSERVATIONS ON HUMP BUILD-UPS In approximately 75 per cent of the wells in a medium-sized field in South Texas the pressure build-up curves rise to 3 maximum and then decline to what appears to be a stabilized reservoir pressure as shown in Fig. I. Similar behavior has also been observed in wells in other fields in Texas and Louisian:,. Because of the possibility of such behavior being a mechanical failure of the pressure bomb. the bombs were run in tanden into one well which had been known to exhibit this behavior. When both bombs recorded identical curves with the characteristic hump, it was inferred that the unusual build-up is a result of well or forma-

By M. R. Tek, K. H. Coats, D. L. Katz

In natural gas storage operations, seasonal pressure fluctuations in the gas reservoir cause the water from the surrounding aquifer to flow into and out of the gas sand. The theory of unsteady-state liquid flow through porous media developed by Van Everdingen and Hurst' has been applied to predict the water movement into and out of the gar bubble for .several .postulated pressure cycles. Those cycles with as many pound-days above the original aquifer pressure as pound-days below, may cause the gas reservoir to slowly grow or shrink rather than hold to a constant volume. Applications to an actual field case study give the predicted gas reservoir monthly pressures and volumes as compared with the observed monthly pressures and volumes. INTRODUCTION A large number of natural gas storage reservoirs are bounded by or adjacent to large water-saturated formations called aquifers. The presence of these aquifers is usually evidenced by the production of water from wells delineating the gas-hearing sands. Some gas storage reservoirs in use today have been purposely developed on aquifers. Cyclic pressure variations in a gas storage reservoir cause water influx and efflux from the surrounding aquifer. This, in turn, results in a varying gas reservoir volume. The prediction of the effect of this aquifer fluid movement on the size and size variation of the gas reservoir can provide information valuable in the study of several reservoir engineering problems. The economics of gas storage operations are directly influenced by the influx of aquifer fluid into the reservoir since a shrinking storage reservoir requires increasing pressures for the storage of the same quantity of gas. Material balance and reserve or recovery calculations also obviously require knowledge of reservoir volume. Various other reservoir engineering calculations such as interpretation of well interference data, evaluation of physical characteristics of porous media, water coning, gas injection and pressure maintenance studies are typical problems where variations of volume due to edge or bottom-water encroachment becomes important. The solution of equations describing influx and efflux of water in the gas sand, application to gas storage reservoir behavior and interpretation of results are the main objectives of this paper. While this work was directed toward gas storage reservoirs, it applies equally well to repressuring of partially depleted oil reservoirs with an active water drive. MATHEMATICAL ANALYSIS As just mentioned, the cyclic pressure variations in an aquifer-surrounded gas storage reservoir cause water movement and varying gas reservoir volume. Prediction of this volume variation from the reservoir pressure cycles or the injection-production schedule requires adoption of a flow model, certain assumptions concerning reservoir and fluid properties, and formulation and solution of the governing differential equation. The flow model used in this study is a gas bubble located on a very large aquifer (Fig. la). Gas bubble radius and aquifer sand thickness are denoted by rr and h, respectively. Fig. lb shows the horizontal, radial two-dimensional porous reservoir attended by the unsteady, compressible motion of the aquifer liquid. Volume variation of the gas bubble is determined from aquifer water flow across the cylindrical boundary, r = r0

By R. V. Higgins

This paper presents a computer method to obtain the shape factors and equal cell volumes of the channels for any well spacing pattern from a potentiometric model. By using this program the authors have processed the data for the seven-spot, direct line-drive and the staggered line-drive patterns. The data for the five-spot pattern had been previously processed by a noncomputer method and are included for completeness. The shape factors and volumes for the channels are presented in tables for those who want to use them to process data using their own permeability relationships and viscosities of their reservoir oils. The authors have used the data and sets of representative permeability curves to process sample calculationr of waterflood performances. The comparison of the calculated results shows that the influence of well spacing is small. The permeabilities of the reservoir rock to oil and water had a greater influence on oil recovery for a given pore-volume throughput of water than the well spacing pattern. The more water-wet the reservoir rock, the better the possibility of permeabilities which are conducive to good recovery. The viscosity of the reservoir oil also influences the recovery more than the well spacing pattern. The reduction in the percentage recovery of oil with increase in viscosity of the reservoir oils is small when oil viscosities are in the range of 0.1 to 3. Above this range the reductions in recoveries are extensive. Sample comparisons of the time required for different patterns to recover the oil are presented. Results of an example calculation are given to show the effect of the permeability profile on recovery. INTRODUCTION The effect of well spacing pattern on the recovery of oil when flooding with either gas or water has been studied by many investigators. Muskat et al.1 presented an analysis using conductivity, sweep efficiency and unit mobility to the time of breakthrough. Dyes et al.2 used experimental techniques (X-ray shadowgraphs) and different mobility ratios. They presented quantitatively the relationship between mobility and sweep efficiency at and after breakthrough. Hauber3 presented a method to predict waterflood performance for arbitrary well spacing patterns and mobility ratios. Craig et al.,4 using techniques similar to Dyes et al. to determine sweep efficiency for a five-spot pattern, added the use of relative permeability curves at breakthrough and thereafter. Douglas et al.5 sed rela- tive permeabilities and continuously changing saturations throughout the entire five-spot flood pattern. In obtaining their solutions they used finite-difierence equations. Hig-gins and Leighton6,7 also used relative permeabilities and continuously changing saturations throughout the pattern before and after breakthrough. They employed techniques that process a flood-pattern calculation on the computer in about one minute. The methods of Douglas et al. and Higgins and Leighton both checked closely the laboratory results for a wide range of mobility ratios. This paper presents some sample performances calculated by the Higgins and Leighton method that show the effect on recovery of different permeabilities and viscosities using the seven-spot, the line-drive and the staggered line-drive, as well as the fivespot flood pattern. No previous paper has presented these data using different permeability curves and continuously changing saturations throughout the flood patterns. The paper also presents (1) the results and analyses of the flood-pattern prediction, (2) the computer techniques for determining the shape factors and volumes from the potentiometric models for the foregoing flood patterns, and (3) the shape factors and volumes of the channels of the flood pattern in the event reservoir engineers may like to process waterflood calculations using their own permeability curves and reservoir oils. DESCRIPTION OF METHOD VOLUMES AND SHAPE FACTORS The use of channels taken from a potentiometric model (see Fig. 1) to aid in calculating the performances of water floods of nonlinear patterns has been thoroughly explained in the literature.0,7 Therefore, very little theory, discussion, or proof regarding this phase will be repeated in this paper. The computer method presented in this paper to calculate the volumes and the shape factors of the channels of potentiometric models employs the trapezoidal rule for the volumes and the Pythagorean theorem (the hypotenuse equals the square root of the sum of the squares of the two sides of a right triangle) for the shape factors. In calculating the volume of a channel, the area of each cell in the channel is determined and then multiplied by the thickness to obtain the volume. In determining the areas of the cells, trapezoids are constructed whose vertical sides are spaced Ax apart, as shown in Fig. 2. The length of the sides is the difference between an ordinate cut off by the top and bottom of the cell — usually equipo-tentials. The coordinates of the points along an equipo-tential or streamline are obtained by Lagrange's8 equation of interpolation for which the constants are coordinate points at the intersections. Three intersections for con-

Jan 1, 1965

By R. A. Greenkorn, R. C. Smith

Hele-Shaw cells are used to model creeping flow through porous media (where Darcy&apos;s law is valid). The effects of inertia on flow about obstructions in a Hele-Shaw cell can be calculated by a perturbation method if one can determine a solution to Laplace&apos;s equation. Results of a computer solution for flow about circular, square and elliptical obstructions are presented. These results show that or a modified Reynolds number of less than I, the inertia terms are small; and for values of less than 3, the average streamline predicts the ideal flow. Therefore, the analogy might be used for studying flow in porous media up to a modified Reynolds number of at least 3. INTRODUCTION The nature of fluid flow in porous media is of interest in the fields of soil mechanics, ground water flow, petroleum production, filtration and flow in packed beds. Because it is very difficult to study the phenomenological behavior of flow in porous media, homologs and analogs are used to study flow characteristics. A Hele-Shaw model, made of two closely spaced plates - usually glass — is often used as an analogy to two-dimensional flow in porous media. Hele-Shaw2 showed experimentally that the streamline configuration for creeping flow around an obstacle located between two closely spaced parallel plates is the same as for two-dimensional, ideal flow about the same obstacle. Stokes7 verified these observations mathematically. The usual equation of motion for flow in porous media is Darcy&apos;s law. The form of the mathematical statement of Darcy&apos;s law is identical, within a multiplicative constant, to the expression for the average velocity over the plate gap in the plane of a Hele-Shaw model. These models may be used to describe flow in both homogeneous and heterogeneous porous media. In the mathematical proof of the Hele-Shaw analogy it is assumed that the convective terms in the Navier-Stokes equations are negligible and that the equations of motion degenerate to Laplace&apos;s equation, with pressure the dependent variable. Whenever a Hele-Shaw model is used as an analogy to flow in porous media, .the validity of this assumption is in question. Riegels5 showed that if convection is not neglected, the velocity distribution around a cylindrical obstruction in the flow field depends on a Reynolds number, the plate spacing, and a dimension characteristic of the obstacle. Riegels&apos; solution, a perturbation solution, uses the boundary condition that the flow rate into the obstacle averaged over the plate gap at any point on the obstacle is zero. The method requires that a solution to Poisson&apos;s equation for the perturbation pressure be found. Riegels evaluated this solution for the case of the cylindrical obstruction. His method may be simplified by eliminating the need for solving Poisson&apos;s equation for the perturbation pressure. Instead, an analytic expression for the perturbation pressure gradient is obtained (valid for arbitrary shapes) and used to eliminate pressure from the equations for the perturbation velocities. The results show, for symmetrical shapes, that if NRe < 1, the convective acceleration terms are small, and that the average velocities represent ideal flow up to at least NRe = 3, where: L is a characteristic dimension of the obstacle perpendicular to flow, b is the plate spacing, u is viscosity, v, is velocity of approach and p is density. THEORY* Incompressible fluid flow through a Hele-Shaw model and a porous medium are analogous because the equations for the average velocity components in two-dimensional flow are identical if:

Jan 1, 1970

By G. B. Thomas

Several methods of analyzing pressure build-up data in wells have been presented by various authors. This paper reviews the theory and method of D. R. Horner and presents example calculations performed on data obtained by testing several different types of wells. These calculations include, (1) graphical estimation of final static pressure, (2) determination of the productive capacity of the pay away from the well bore and. (3) the degree to which the formation adjacent to the well bore has been damaged by completion or other causes. The methods of testing and precautions which should be taken to assure the best data possible are discussed. Limitations and reliability of calculated results are also treated. INTRODUCTION Pressure testing of wells is generally limited to the determination of producing and static mean formation pressures. The so-called "static" pressure determination.. along with PVT, electric log and production data, enable the reservoir engineer to determine, within reasonable limits, the drive mechanism of the reservoir and in some cases. the amount of edge water encroachment. Producing pressure tests enable calculation of productivity indices and allow the engineer to plan the systematic production of a pool for optimum conservation of suhsurface energy. The radial flow formula advanced by Muskat&apos; has been based on the assumption of incompressible radial fluid flow. It has been known that reservoir fluids do not behave in an ideally incompressible manner. For example, incompressible flow theory indicates a simple logarithmic relationship between the difference of the instantaneous and static well pressures when plotted against time. &apos;The latter stages of this type of plot of pressure build-up data generally show a marked devia-, tion from the earlier straight line trend. which deviation may be shown to be due to the compressible flow of fluids toward the well bore. Shut-in times of 24, 48. 72. or at most 96 hours are currently in wide use for determining so-called "static" reservoir pressure. Due to the continuation of compressible flow of fluids into the well bore long after this arbitrarily taken shut-in time, the determination of static pressures has almost invariably resulted in lower than equilibrium values. Materials balance calculations made early in the life of a reservoir often result in a calculated reserve which later observations prove to be too low. Failure to obtain reliable "static" reservoir pressures within the prescribed 24 or 48-hour build-up period has undoubtedly been a major factor in obtaining these low estimates. Comparison of theoretical and actual productivity performance has indicated that formation damage and not being able to attain true static pressures have been partially responsible for the observed discrepancies. Research into the theory of compressible flow behavior has resulted in methods of applying these theories to the testing of wells. Application of the developed theories to pressure build-up performance indicates that in many cases the following can be estimated. 1. The static reservoir pressure. 2. The in-place effective formation permeability away from the well bore. 3. The degree or extent by which the formation has been damaged adjacent to the well bore, either through completion methods or subsequent damage due to fluid entry. REVIEW OF THEORY Horner2 has shown that the pressure build-up within a "point source" well is approximated by the following formula: Equation (1) is the approximate "point source" solution to the radial compressible flow equation advanced by Muskat. The solution assumes the following: 1. A point source well is producing at constant rate from the center of an infinite reservoir with a constant pressure at its external boundary. 2. The fluid flowing is present in one phase only. 3. The compressibility and absolute viscosity of the fluid remain essentially constant over the range of temperature and pressure variation encountered. 4. The well i.; shut-in at the sand face arid there is 110 after production into the well bore. 5. The formation permeability is homogeneous in the direction of flow. From Equation (1) it can be seen that if an ideal well were shut-in while producing from a reservoir under the conditions assumed, the pressure build-up would be a logarithmic func-

Jan 1, 1953

By G. W. Lester, T. J. Nowak

The paper presents a practical method of utilizing pressure fall-off data obtained when a water injection well is shut in for determination of: (I) the static reservoir pressure, (2) the potential water injective capacity (md-ft) of the formation, and (3) the formation damage (skin effect) surrounding the wellbore. In the development of the method, equations are presented which give The theoretical background for predicting the behavior of water injection wells. Certain physical factors which limit the application of the method are discussed. The procedure and the relevant precautions in obtaining the data are discussed. A sample analysis and computation from an actual pressure fall-off curve is illustrated. Fintrlly, general characteristics of pressure fall-off curve! obtained in several reservoirs are presented. INTRODUCTIO N Theoretically, in an injection well the flow of water into the formation will continue as long as the sand-face pressure is greater than the reservoir pressure. However, in practice, engineers have observed that upon stopping injection and opening the wellhead to atmospheric pressure, some wells backflow water for time intervals ranging from several seconds to several hours although the sand-face pressure, as calculated from the hydrostatic head, is considerably greater than the static reservoir pressure. Further, field observations indicated that upon shutting in a water input well, the wellhead pressure does not drop off immediately, but a pressure decline occurs lasting from several seconds to several days. These observations make it evident that when an input well is shut in, a residual pressure lingers in the vicinity of the wellbore which is greater than the average reservoir pressure. This residual pressure declines first rapidly, then more slowly, to the static reservoir pressure: therefore, it is called transient back-pressure or pressure fall-off. The first step in investigating the mechanism of pressure fall-off is to establish what general conditions are necessary for its occurrence, and to this end, Dickey and Andresen1 have presented a qualitative explanation describing the mechanism. They point out that the formation pressure gradient in the vicinity of the well does not disappear immediately when the input well is shut in, as if the system were filled with a totally incompressible fluid, but it is maintained by the expansion of the compressible fluids and by the elasticity of the formation. In such a compressible system the decrease in pressure is possible only to the extent that the water front continues to advance until the pressure is equalized. The principal aim of this paper is three-fold; namely. (1) to present a theoretical quantitative analysis of pressure fall-off data, (2) to present an application of the developed theory to field data, and (3) to discuss the general characteristics of actual pressure fall-off curves. Under (I), theoretical equations dealing with such factors as compressibility of the fluid-reservoir system, mobility of the fluids and certain physical restrictions that

Jan 1, 1956