Search Documents

By R. W. Parsons

The over-all apparent single-phase permeability of fracture-rock systems was studied using two different two-dimensional models. In a strict sense the results are applicable only to these models, yet may give some insight into real reseruoir behavior. The regular fracture-matrix model can have any number of regular fracture sets. All fractures within each set have the same orientation, width and spacing. The gross behavior is amlytically equivalent to that of an anisotropic permeable medium. In the heterogeneous fracture systems, individual fracture conductivities vary within a two-dinensioml network pattern. A computer program calculates the overall permeability after the conductivities are randomly placed in the network. The calculated permeability depends upon the distribution of these conductivities and the particular fracture pattern into which they are placed. Although detailed information on actual reseruoir fracture systems is not available, this study shows how some of the possible variables may affect obserued gross reservoir characteristics. INTRODUCTION Characteristics of reservoir rocks are determined by direct (core analysis) or indirect (logging and production tests) examination. Both methods, singly and together, are widely used in reservoir evaluations. For fractured reservoirs a problem exists in relating characteristics of cores to the in situ reservoir properties. Knowing fracture characteristics can be of value in determining production mechanisms. Because of the limitations of observation, the intimate details of a fractured rock system will never be known. One approach in analyzing fractures is to study the behavior of some conceptually simple models. The two models analyzed contain features which may be realistic, yet are mathematically tractable. Flow is considered in two dimensions only, with the plane of interest at right angles with the fractures. This concept represents the horizontal flow in a reservoir system with the fractures oriented vertically. Only single-phase, laminar flow is considered. The regular fracture-rock system permeability is calculable from a simple analytical expression. The heterogeneous system requires a computer program for calculating overall permeability.

Jan 1, 1967

By C. F. Weinaug, R. W. Farley, J. F. Wolfe

A rapid, accurate method for predicting the dew points of gas condensate systems and their subsequent normal and retrograde phase behavior with pressure decline has been developed. The method predicts dewpoint pressure, volume percent liquid, condensate content of produced gas, and the liquid and gas compositions at pressures below the dewpoint pressure. Use of a modified form of the BWR equation of state, which was established by regression calculations on experimental data from several gas condensates, provides depletion behavior predictions that are significantly improved over predictions using K-values based on convergence pressures and Standing's correlation for density. INTRODUCTION Prediction of complex hydrocarbon-mixture thermodynamics is fundamental to reservoir and gas process engineering. In particular, accurate descriptions of phase behavior are needed to plan production operations for gas-condensate reservoirs. Although analysis of reservoir behavior during producing and cycling operations can be made for simple reservoir models, extensions of these studies to more complex reservoir systems are still difficult, both with regard to numerical techniques and the phase behavior of the fluids. This report presents a mathematical method for handling the complex phase relationships encountered in retrograde condensates; it represents an advance in reservoir behavior technology in that it will aid the development of closed-form mathematical representations of total production systems. A typical phase diagram for a complex hydrocarbon mixture is shown in Fig. 1. When the reservoir temperature T, is higher than the critical temperature Tc, and lower than the critical-condensation temperature Tct' the mixture is a gas condensate. Decreasing the pressure on a condensate while maintaining nearly constant temperature will result in two phases forming at the dew point. Further decrease in pressure will cause additional liquids to be formed. This phenomenon is called retrograde condensation. As pressure continues to decrease, normal phase behavior regions are encountered (liquid vaporizes as the pressure decreases) until the low-pressure dew point is encountered. When a gas-condensate reservoir is produced, gas is withdrawn and the volume of fluids in the reservoir remains nearly constant as the pressure declines. Heavier components tend to accumulate in the liquid phase as a result of retrograde condensation, while lighter components tend to concentrate in the vapor phase. Thus, as the gas is produced the composition of the gas and the condensate

Jan 1, 1970

By H. N. Hall

Reservoir and producing characteristics can govern the decision to use either a one-, two- or three-dimensional model for making predictions for gravity-drainage reservoirs. Examples of conditions requiring one-, two- and three-dimensional calculations are given. In 1961 the author presented a method for predicting one-dimensional gravity-drainage performance. This work has been exrendcd to obtain a three-dimensional model, which is described and for which a sample problem is presented. Introduction Many papers on gravity drainage or gas-cap drive'-3 have been presented for conditions where the reservoir is treated as a single areal segment so that, in effect, only one-dimen-sional flow is considered. Other authors'-Q ave given evidence showing that the one-dimensional concept would be unsatisfactory to predict reservoir performance for certain types of reservoirs. Experience has indicated that specific reservoir and producing conditions govern whether one-, two- or three-dimensional concepts must be used as a basis for reliable reservoir performance predictioils. Previously the author presented a paper outlining a method for predicting the performance of gravity-drainage reservoirs.' That paper treated the reservoir as a single areal segment so that it was applicable for only one-dimensional flow. The work has been extended to obtain a three-dimensional model that can be used to predict performance of gravity-drainage reservoirs. Although capable of being three-dimensional, it can be used equally well for one- and two-dimensional predictions. This paper discusses those reservoir and producing conditions governing the choice of model, and also describes a three-dimensional model suitable for predicting reservoir performance in the more complex situations. Discussion The Choice of Model To predict reservoir performance, it is necessary to represent the reservoir by a physical model. A mathematical model, based on the assumed physical model, then is used to predict the movement of oil and gas within the reservoir. The word "model" is used here to connote the combination of both a physical model and a mathematical model based on the physical model. Certain types of reservoir situations may be predicted suitably with a one-dimensional model, whereas two- and three-dimensional models may be required for other situations. This section illustrates reservoir conditions associated with the use of a specific type of model. Many types of reservoir structures are such that gravity drainage possibly could be an important factor in oil production. Production practices, as well as the nature of the structure. can influence the type of model that should be used. obviously, it is impossible to discuss all combinations of reservoir structure and production practices that might be encountered in gravity-drainage reservoirs. The following are typical examples of conditions where a choice must be made among a one-, two- or three-dimensional model to predict reservoir performance satisfactorily. Conditions Requiring a One-Dimensional Model Fig. la shows contours for an asymmetrical anticline in which gravity drainage might be expected to be important. Fig. Ib depicts a cross-section along Line A-A in Fig. la. Existence of a gas cap is a good clue that gravity drainage will be an important factor in producing oil from this type of reservoir. However, it is not essential since a gas cap will form as oil is produced from the reservoir. If the proper combination of reservoir permeability and withdrawal rates is present, the reservoir will behave like a

Jan 1, 1969

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

By G. Thodos, R. J. Smialek

A method has been developed that permits the correlations of Mehra and Thodos8 and Dastur and thodos2 for binary systems to be used for the prediction of the vapor-liquid equilibrium constants of the light and heavy constituents, Kland Kh, of a ternary hydrocarbon mixture. Furthermore, the unique relationship observed in this study, i.e., a plot of log Ke/Ki vs log Kp/K h is linear and passes through the origin, permits the calculation of the K-value for the intermediate component. This approach has been used to establish K-values for the constituents of the three ternary hydrocarbon systems used to develop this method and has also been applied to two additional ternary hydrocarbon systems to yield an average deviation of 7.7 per cent resulting from a total of 294 points for the five ternary systems considered. INTRODUCTION The rigorous calculation of vapor-liquid equilibrium constants for systems containing two or more components is a difficult problem and becomes extremely complex at elevated pressures and for conditions approaching the critical point. Although thermodynamic relationships can be developed to account for the vapor-liquid equilibrium behavior of simple systems, their use is limited because the necessary PVT data for these systems are not generally available. To establish a method for calculating K constants directly, without any recourse to the PVT behavior of a system, Mehra and Thodos8 applied the principle of corresponding states to the vapor-liquid equilibrium constants of binary hydrocarbon systems. They developed a series of charts which made possible the prediction of K-values for all binary hydrocarbon systems not containing methane. This work has been extended by Dastur and Thodos2 to include binary hydrocarbon systems containing methane. The charts developed in these two studies cover the high-pressure region and are applicable to conditions in the vicinity of the critical point. The ability of this approach to predict vapor-liquid equilibrium constants of binary hydrocarbon systems suggests that this method may be extended to the prediction of the K-values of the components of ternary hydrocarbon systems. In this investigation, the charts developed for both the methane and non - methane binary systems,2,8 together with mixture rules developed in this study, have been used to predict the K-values of the light and heavy components of a ternary hydrocarbon system. These calculated values are then used to establish the K-value of the intermediate component through the use of a novel relationship developed in the present study. DEVELOPMENT OF METHOD FOR PREDICTING K-VALUES OF TERNARY SYSTEMS The application of the corresponding state charts developed for binary hydrocarbon systems is restricted to those systems containing two components. The extension of these correlations to three-component systems requires that the ternary system be reduced to an equivalent binary, which can be treated in a manner consistent with the corresponding state development. To establish such an approach, the vapor-liquid equilibrium constant of the light component K is first calculated by assuming that the ternary mixture consists of the light component and a pseudo-heavy component of the intermediate and heavy constituents of the original ternary. To establish the vapor-liquid equilibrium constant of the heavy component Kh, the light and intermediate constituents of the ternary are combined to form a pseudo-light component of another equivalent binary system. The calculated values Kl and Kb are then utilized to calculate the vapor-liquid equilibrium constant of the intermediate component EQUILIBRIUM CONSTANT FOR LIGHT COMPONENT K To calculate the K value of the light (most

Jan 1, 1965

By J. J. Cosgrove, J. E. Warren

A general model which approximates the effect of cross- flow has been developed to give a practical method for predicting the waterflood behavior of a stratified reservoir. The model is based on a modification of Dietz's theory and allows for variations in both the permeability and hydrocarbon pore volume. In the particular cases considered, it is assumed that the permeability can be characterized by a log-normal distribution and the hydrocarbon pore volume by a normal distribution. A simple graphical method which enables the practicing engineer to predict the behavior of a stratified system is presented The results obtained by the proposed method are compared with those obtained by the Dykstra-Parsons method. As a result of this study the following conclusions have been drawn: (1) The effect of cross-flow in a stratified system can be appreciable, particularly at very favorable or very unfavorable mobility ratios; (2) Under normal conditions, the effect of variations in the hydrocarbon pore volume can be neglected; and (3) The failure to use all of the available permeability data can lead to large errors in the prediction of the behavior of a stratified reservoir. INTRODUCTION Various methodshavebeenproposedto characterize and to predict the waterflood behavior of a stratified reservoir. Most of the methods assume that the reservoir is composed of discrete homogeneous continuous layers. With this model, the degree of stratification can be measured by several parameters based on core-analysis data. Among these are the Lorenz coefficient1 and the variation of a log-normal permeability distribution? The behavior of stratified systems is usually predicted by the Stiles, Dykstra-Parsons method or some modification of these. In both of these methods the reservoir is divided into discrete homogeneous layers with no cross-flow between layers. In the Stiles method the mobility ratio is assumed to be equal to unity while in the Dykstra-Parsons method it is allowed to take on any value. Other predictive methods have been proposed; method allows for cross-flow between the beds and a method by Schmalz and Rahme6 correlate.s the recovery directly with the Lorenz coefficient. In this new approach, which is essentially a continuous analog of Hiatt's method, the effects of both mobility ratio and cross-flow between the beds have been included. It is assumed that the permeability can be represented by a log-normal distribution and the hydrocarbon pore volume or porosity by a normal distribution. These types of distributions have been observed by several authors, 1,2,4 and most field data seem to confirm their observations; e.g., if these distributions are assumed and there is a 1:l correspondence between porosity and permeability samples, the commonly encountered linear relationship between porosity and the log of permeability is obtained. Since, in many field cases the permeability and poro sity data are truncated or cut-off at predetermined upper and/or lower values, the effect of discarding part of the data was investigated with the proposed method. A better technique for truncating the data is suggested. THEORY In the derivation of the method to be described in this paper the following assumptions are made: 1. Capillary forces are negligible per se; their effects are only manifest in the relative permeability curves. 2. The fluids are immiscible, incompressible and homogeneous. 3. The reservoir is horizontal and uniformly thick; it is initially liquid saturated i.e., connate water and oil. 4. The reservoir is composed of discrete layers, each having its own permeability, porosity and connate water saturation; each layer is homogeneous,

Jan 1, 1965

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas'-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller's results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By R. C. Earlougher, F. G. Miller, T. D. Mueller, H. J. Ramey

There are many studies of flow in radial systems that can be used to interpret unsteady rerervoir flow problems. Although solutions for systems of infinite extent can be used to generate solutions fur finite ow systems by superposition, application is tedious. In this paper a step is made toward simplifying calcu1ations of such solutions for finite flow systems. Superposition is used to produce a tabulu-lion of the dimension1e.s.s pressure drop function at several locations within a bounded square that has a well at its center. The square system provides a useful building block that may be used to generate flow behavior for any rectangular shape whose sides are in integral ratios. Values of the tabulated dimensionless pressure drop function are simply added to obtain the dimensionless pressure drop function for the desired rectangulm system. The rectangular system may contain any number of wells producing at any rates. Furthermore, the outer boundaries of the rectangular system may be closed (no-flow) or they may he at constant pressure. Mixed conditions also may be conyidered. Tables of the dimensionless pressure drop function for the square system are presented and various applications of the technique are illustrated. Introduction In 1945, van Everdingen and Hurst' published solutions for the problem of water influx into a cylindrical reservoir. Since this problem is mathematically identical with the depletion of a cylindrical reservoir with a well at the origin, the van Everdingen-Hurst solution may be used to study the depletion problem. In their analysis, they assumed that the fluid had a small, constant compressibility such that flow was governed by the diffusivity equation For a constant production rate q starting at time zero, van Everdingen and Hurst showed that the unsteady pressurc distribution for both finite and infinite systems could be expressed in terms of a dimensionless pressure 'Sabulations of the dimensionless pressure drop for a unit value of r, were provided by van Everdingen and Hurst,' and later by Chatas.' Others also presented values in graphical or tabular form." If the radius of the well becomes vanishingly small, r,+ O, the line source solution may be used for Eq. 2 when infinite systems are considered. Eqs. 5 and 6 are excellent approximations for Eq. 2 under certain conditions: In 1954, Matthews, Brons and HazebroekV emonstrated that solutions such as Eq. 5 can be superposed to generate the behavior of bounded geometric shapes; i.e., the behavior of a bounded single-well system can be calculated by adding together the pressure disturbances caused by the appropriate array of an infinite number of wells producing from an infinite system. These wells are referred to as image wells. Matthews, Brons and Hazebroek considered systems containing a single well producing at a constant rate. This superposition can be represented analytically as

Jan 1, 1969

By H. K. van Poollen, H. C. Bixel

A treatment is given of the transient pressure behavior of a well located at the center of a circular region surrounded by a radial discontinuity. On either side of the discontinuity, the values of permeability, viscosity, compressibility and porosity are uniform but may be different from those on the other side of the discontinuity. The results are obtained by solving a pair of finite difference equations. The numerical solution to these equations is obtained using a digital computer. The results are compared with previously published, approximated analytical solutions to the same problem. Solutions, presented graphically, show pressure decline for constant rate fluid production. The range of variables studied include dimensionless time from '0.001 to 100 and storage capacity ratio from 0.001 to 1,000. The well behaves as if it were in an infinite reservoir for dimensionless times less than 0.25. Reservoir properties near the well can be estimated in the usual manner. An overlay technique is used to match an experimental curve with one of the theoretical curves. It is possible to estimate the distance to a discontinuity by substituting the actual time t and the corresponding dimensionless time tD at which a match occurs into the equation a = (0.000264 tk1/tD , µ1 cI)I/2 where kI/ I µI cI is the diffusivity near the well, and may be estimated from data taken at early time. Several buildup curves are computed. These curves show that for early shut-in times, correct values for transmissibility are obtained from conventional analysis. However, erroneous values of static reservoir Pressure are obtained unless data at large shut-in times are taken. INTRODUCTION A mathematical treatment of the transient pressure response of a well located at the center of a region bounded by a circular discontinuity is given. Within a region (Fig. 20) the properties of both the rock and the fluid are considered to be constant on either side of the discontinuity. While these properties are considered constant, they may be different on opposite sides of the discontinuity. A discontinuity of this type could be a fluid-fluid contact or a sudden change in rock characteristics such as thickness, porosity or permeability. Analytical solutions to this problem are available. 1-3 However, they are so involved that they are of little practical use. An approximate solutionl is available but the range of times over which it is valid has not been specified. A numerical solution to a pair of finite difference equations is used to obtain the solution given in this paper. PREVIOUS WORK One of the first solutions to the problem was published by William Hurst.1 Hurst considered unsteady-state flow of fluids through two sands in series with different mobilities in each sand. He used Laplace transforms to obtain a solution for a single well located at the center of the circle enclosing the first of the two sands. In this case, the solution for the pressure change in Region I (Fig. 20) is2

By D. M. Beeson, G. D. Ortloff

Water-propelled banks of carbon dioxide recovered both high- and low-viscosity crude oil substantially in excess of that recovered by water flood in linear flow model experiments. The increase in oil recovery is explained by oil swelling and viscosity reduction. Additional oil was recovered by solution gas drive at the end of the carbon dioxide floods when the model systems were depres-sured. This study was an investigation of the use of a water-driven carbon dioxide bank to recover crude oil. Techniques involving the injection of carbon dioxide have aroused con-siderable interest in recent years among methods considered for increasing oil recoveries. A number of investigators have studied the process in the laboratory1-' and in the field.7,5 Laboratory studies referenced indicate that carbon dioxide in conjunction with water flooding recovers more oil than does water flooding alone. These re- ports diverge widely, however, in stating magnitude of improvement. We believed that a water-driven carbon dioxide bank would increase recovery of crude oils by two mechanisms: viscosity reduction and oil swelling. Due to reduction in the oil viscosity caused by solution of carbon dioxide, displacement efficiency should be improved. In addition, the residual oil should be highly swollen by the carbon dioxide held in solution. Thus, the final residual oil saturation, in terms of stock-tank oil, would be lower than in the absence of carbon dioxide. To test this hypothesis, experi-

By A. B. Dyes

In the operation of a water-driven reservoir, a free gas saturation can he established by maintaining production rates fast enough to cause the reservoir pressure to decline below the bubble point. The benefit of such a procedure on the displacement efficiency of the oil by water is illustrated for two types of sandstone samples from one reservoir. Those rock samples showing the poorest recovery to water drive in the absence of a free gas saturation give the most improvement in the presence of a free gas saturation. Employing one type of reservoir fluid (low .shrinkage) the benefit of evolved gas on oil displacement is calculated at several reservoir pressures below the bubble point. This demonstrates the presence of air optimum pressure for water displace,rnent which is .several hundred pounds per .square inch below the bubble point. Displacement of oil at the optimum pressures results in the removal of 7 to I? per cent more oil than would be obtained by displacement at the bubble point pressure. In some instance. this magnitude of increase may be worthy of consideration in establishing the maximum efficient rate of production from many water-driven fields. INTRODUCTION For many years the belief that the use of gas could be combined with that of water to result in a better recovery operation than the use of either of these phases alone has been prevalent in the oil industry. This is evidenced by the several field trials which have employed the injection of gas and water in combination.1.2.v In most cases the small amount of any additional recovery, in combination with the limited knowledge available on the older reservoirs, has left uncertain the influence of gas on recovery in these earlier attempts to combine the use of gas with water. Research studies in recent years have demonstrated through laboratory. floods the magnitude of the benefit of the presence of gas on the displacement of oil from sandstones by water.5,6,7,8,9,10,11 In addition, studies of the distributions of the phases which are present during the displacement process have developed concepts concerning the mechanism by which this benefit is obtained."." It now remains for the engineer to incorporate these findings into the general technology of reservoir operations. It is the purpose of this paper to illustrate that the production of a water-driven reservoir at pressures below the bubble point is one such operation whereby improved oil displacement efficiency can be obtained as a result of gas evolution. The establishment of a gas saturation by evolution of dissolved gas was one of the techniques used by several investigators in the laboratory studies. The possibility of obtaining the benefit of gas evolution in water-driven reservoirs by producing at pressures less than the bubble point has 3150 been pointed out previously." The desirability of this operation has undoubtedly

Jan 1, 1955

By H. J. Ramey

Conventional calculation of total system isothermal compressibility for a system containing a free gas phase involves, among other things, evaluation of the change of oil and gas formation volume factors and the gas in solution with pressure. Preferably, this information should be obtained from laboratory measurements made with particular oils and gases. Often, experimental measurements are not available. In this case, it is necessary to obtain pressure-volume-temperature relationships from general correlations such as those of Standing' for California oils. In order to speed estimates of compressibility, generalized plots have been prepared of the change of both oil formation volume factors and gas in solution, with pressure from Standing's correlations. A generalized plot for estimating the change in the two-phase (oil and gas) formation volume factor with pressure is also presented. Usually, the effect of gas dissolved in reservoir water upon the total system compressibility is neglected for gas saturated systems, due to the low solubility of gas in water. Results of this study indicate that the increase in total system compressibility caused by solution of gas in water is often as large as the compressibility of wafer, and can be magnitudes larger for low pressure systems. Generalized results for estimating the change of gas in solution in water with pressure are presented in tabular and graphical form. INTRODUCTION During the past decade, pressure build-up and drawdown techniques have gained an important place in reservoir engineering. Build-up and drawdown analyses are only two special applications of the broad field of transient fluid-flow theory. All solutions of transient fluid-flow problems contain a parameter called the total system isothermal compressibility. This property of fluids and porous rock is a measure of the change in volume of the fluid content of porous rock with a change in pressure, and it may vary considerably with pressure. Evaluation of total system isothermal compressibility is not difficult, but it is tedious and time-consuming. Often compressibilities are estimated roughly, or transient flow methods are neglected completely. The benefits of using accurate system compressibility in properly-executed build-up or drawdown analyses are: 1. Better planning of pressure build-ups may be achieved to avoid unnecessary loss of revenue due to excessively long shut-in periods, or to shut-in periods too short to yield useable data. 2. Better and more reliable estimates of static formation pressures for reserves estimates and rate performance estimates. 3. Reliable information for evaluation of well completion effectiveness, and planning and interpretation of well stimulation efforts. The purpose of this paper is to clarify the nature of the total system isothermal compressibility, and to present useful methods for estimation of compressibility, particularly for systems containing a gas phase. DEVELOPMENT Numerous publications have presented solutions to transient single-phase flow of slightly compressible fluids, stressing pressure build-up applications. In transient flow, a compressibility* term arises to permit volume content of fluids in porous rock to change as pressure changes. The basic nature of the compressibility term is usually taken for granted. Problems arise in practical applications of transient fluid theory because most published works consider only one flowing fluid—in an ideal porous system containing only one fluid. In 1956, Perrine2 presented an intuitive extension of single-phase flow pressure build-up methods to multiphase flow conditions. Later, Martin- established conditions under which Perrine's multiphase build-up method had a theoretical foundation. Perrine has shown that improper use of single-phase build-up analysis in certain multiphase flow situations can lead to gross errors in estimated static formation pressure, permeability and well condition. It is likely that, much pressure build-up data for oil wells should be analyzed on the basis of multiphase flow. For both single-phase and multiphase build-up analysis, the isothermal compressibility term in dimensionless time groups often should be interpreted as the total system compressibility. All real reservoirs contain one or more compressible fluid phases. In addition, rock compressibility can contribute in an important way to the total system compressibility. The proper total system compressibility expression may contain terms for compressibility of oil, gas, water, reservoir rock and terms for the change of solubility of gas in liquid phases.

Jan 1, 1965

By L. H. Reiss, J. Groult, L. Montadert

Many fields, where the reservoir is composed of sandy layers, show great complexity because of the lack of continuity which results from a particular type of seditnentation. This complexity may be a factor of prime importance in studying reservoir drainage, especially in the displacement of oil by water or gas. This paper describes work done to study this heterogeneity and its influence in the exploitation of the fields of the lower Devonian sandy shale in the Fort Polignac Basin in Algeria. Two lines of research were followed. 1. On the border of the hasin, a detailed sedimentological study established rules for the distribution of the sandy zones, which furnished hypotheses necessary for devising a geologic model of the reservoir. A certain number of sandy zones were recognized and characterized by their sedimentation. 2. In the interior of the basin, production logs deterrnined the productivity profile. Study of these point data led to a description of the distribution of oil and gas phases in the reservoir so that the continuity of different zones between the wells could be deduced. These techniques were applied to the study of the lower Devonian sandy shale reservoir of Zarzaitine. The sedimentological study guiding the correlation work, with the aid of the conventional static data (logs and cores), led to the establishment of a useful geologic model of the reservoir. Production logs, thanks to the contribution of dynamic elements, permitted corrobora-tion of the model. In particular, the logs confirmed separation, as proposed by the sedimenfological study, of two comparable sandy groups by revealing the continuity of a thin impermeable layer, invisible on conventional logs. INTRODUCTION The heterogeneity of sandstone reservoirs considerably influences the distribution of fluids, as well as their displacement in the reservoir. Rational development of a field and maximum recovery of hydrocarbons are closely tied to an exact and precise knowledge of the intern architecture. We will try to show how it is possible I refine the knowledge of the reservoir by seeking qui different and varied sources of information. Productic geology, the group of techniques which leads to idei concerning reservoir architecture and its geometric descris tion, furnishes the framework of this undertaking. Of these techniques. those used so far are qui varied. Descriptive as well as quantitative, the technique partake as much of classical geology as of reservo science. They include (1) geologic studies of the enti basin, especially a stratigraphic and sedimentologic study of the reservoir member at its outcrop, detail1 observations of cores taken from wells and, finall oriented cores; (2) interpretation of geophysical dat (3) interpretation of electric logs, nuclear and sonic lo and dipmeters: (4) petrophysical study of the reserve on samples taken from wells; (5) reservoir dynamics (stu of the production history and, in particular, interpretati of production logs); and (6) analysis of the fluids a thermodynamic studies. This incomplete list could also include (1) study the migration of the fluids, (2) geochemistry (chemi characterization of the oils), (3) regional hydrodynamii and (4) injection and detection of tracers. If most of these techniques are classic, attention i nevertheless be called to applying two of the techniql in conjunction: the study of the reservoir member its outcrop and the study of dynamic distribution fluid phases in the reservoir through use of product logs, used in the study of the sandy shale reservoirs the lower Devonian of the Fort Polignac Basin, Alge in particular, their application to the study of the 2 zaitine field where they have led to a new and satisfact model of the reservoir will be described. THE ZARZAITlNE FIELD Situated in the interior of the Fort Polignac 8 about 150 kilometers (km) north of the lower Devol outcrop, the Zarzaitine Devonian oil field was discovc On Ian. 28, 19'' (Figs. 1 and l8). It is in the fern a monocline truncated by two perpendicular fault: great throw, with a downthrown section to the north\: Erosion at the top of the structure has caused one of it to disappear. It is a reservoir with a gas cap an m coefficient on the order of 1:5, open to the aql On the northeast, with the oil zone covering an Of 106 sq km and with initial oil reserves of 222 mi

Jan 1, 1967

By W. T. Weller

In Part I, a study of pressure build-up curves calculated for conditions under which both oil and gas flow led to the conclusion that the presence of a dispersed free gas phase in an oil reservoir must be taken into consideration to estimate accurately average reservoir pressure and permeability from build-up curves. Familiar methods based on the assumption of no free gas can be extended to the two-phase case by using total compressibility and mobility in place of single-fluid compressibility and mobility. These methods give correct values for average pressure and permeability when gas saturation is small. Errors become larger as the gas saturation increases. However, for the use to which the results will be put, the methods are satisfactory for reservoir engineering purposes. An improved method of calculating the performance of depletion-type reservoirs is presented in Part 2. Because the mathematical relationships describing simultaneous flow of oil and gas are quite involved. simplifying assumptions are made to provide means of obtaining approximate solutions of reasonable accuracy. One such approximate method now in use is the constant-GOR solution. It involves the assumption that at any instant, the ratio of total gas flow rate (both free and disolved) to oil flow rate is the same at all points in the reservoir. The approach is not applicable unless the free gar saturation in the reservoir is everywhere greater than the critical gas saturation. This paper presents a modification which, by avoiding the constant-GOR assumption, makes the method applicable to all reservoir conditions, and so far appears to be more accurate than the constant-GOR solution and to he comparable in required compuctation time. PART I— BUILD-UP CURVE ANALYSIS INTRODUCTION Pressure build-up characteristics of shut-in wells have been used for many years by engineers to evaluate average reservoir pressure, effective permeability thickness of the pay section, effectiveness of well completion (skin effect) and reservoir size. A number of methods of analysis have appeared from time to time."" Without exception, these methods are based on the assumption that the reservoir contains but one fluid of constant small compressibility and constant mobility. It has been suggesteda" hat these single-fluid methods may be applied to data from reservoirs containing both oil and gas by substitution of some effective total properties of the multiphase system in place of the corresponding single-phase properties. The present investigation was undertaken primarily to evaluate this approach. METHOD A number of theoretical build-up curves were calculated for conditions of two-phase flow, under the assumption of certain reservoir and fluid properties, and were analyzed by single-fluid methods with appropriate total compressibility and total mobility values for the corresponding single-fluid properties. Results of the analyses were compared with the assumed conditions. The theoretical build-up curves were completed by procedures similar to those of West, Garvin and Sheldon." Since these calculations require a considerable amount of computer time, an attempt was made to derive an approximate calculation method. The attempt was unsuccessful for calculating build-up curves, but the effort did result in a new approximate method of calculating the performance of solution gas drive reservoirs, which appears to be an improvement over the constant-GOR method" used previously (see Part 2). The West, Garvin and Sheldon calculations involved the following assumptions: (1) the reservoir is circular and completely bounded, with a completely penetrating well at its center; (2) the porous medium is uniform and iso-tropic, with a constant water saturation at all points; (3) gravity effects can be neglected; (4) compressibility of rock and water can be neglected; (5) the composition and equilibrium are constant for oil and gas; (6) the same pressure exists in both the oil phase and the gas phase; and (7) no afterproduction occurs, i.e., the well is shut in at the sand face for build-up. These assumptions make it possible to describe two-phase flow of oil and gas by the partial differential equations:

Jan 1, 1967

By D. K. Atwood

Experiments resulted in a satisfactory laboratory method for restoring permeability to clay-containing cores damaged by fresh water. Clay contents of a number of field cores were measured, and permeabilities of plugs from these same cores were then deliberately reduced with fresh water. This damage is attributed to swollen and dispersed clays occupying the pore space. After damaging, a number of experiments were performed to meaJure the amount of damage and to establish some means by which permeability could be restored. The experiments included flooding the damaged cores with water-miscible fluids such as salt water, acetone, isopropyl alcohol and ethanol. Permeability was not successfully restored in these experiments. However, part of the damage was repaired by flooding with oil; when water was removed by distillation in the presence of immiscible fluids such as air or toluene, permeability was completely restored. This evidence suggested that swollen and dispersed clays could be collapsed to their original volume by strong interfacial and capillary forces. It was further postulated that the required forces could be generated by flooding the damaged cores with a solvent partially miscible with water. The flooding experiments were repeated using n-hex-an01 as the partially miscible solvent. Permeability was restored to five of six damaged cores and substantially increased in the sixth. A large fraction of the restored permeability was retained even after water saturation was raised to its original value with 12 per cent salt water. INTRODUCTION Sharp reductions in permeability often occur when relatively fresh water contacts clay-containing formations during drilling and workover operations. These permeability losses are caused by removing inorganic ions from the environment surrounding the clay, and consequent swelling and/or dispersion of clay minerals into the available pore space.' This phenomenon is generally termed clay damage, fresh-water damage, or simply formation damage; it causes large losses in current revenue by preventing oil wells from making their allowable production. Attempts to repair the damage and restore permeability by flowing salt water solutions or brines through clay-damaged cores containing montmorillonite have been unsuccessful.' This irreversibility is thought to result from formation of brush-heap, or edge-to-face, structures when the dispersed clay is flocculated. The brush-heap structures occupy much more space than the close packed domains present before damage.' One solution of the problem is to destroy the clay-water brush-heap and thus restore permeability. Because no satisfactory method existed for restoring permeability to clay-containing formations damaged by fresh water, the work described in this paper was under taken. The laboratory experiments generally consisted of deliberately damaging fresh cores containing clay and then attempting to repair this damage by various means. Results indicate that generating strong interfacial forces within the pore space of damaged cores collapses the clay brush-heap and restores permeability. These forces are most conveniently generated by flowing partially water-miscible solvents, such as n-hexanol, through a core. THEORY OF THE DAMAGE PROCESS The most common clay mineral groups known to cause permeability damage to formations are the mont-morillonites, kaolins, chlorites and illites. These clays are constructed of particles which can adsorb water on their surfaces and edges and, in the case of montmorillonite, between layers of the basic particle itself. This adsorption increases as water salinity decreases. At low salinities the particles disperse into the aqueous phase. When the clays present in the formation are kaolin, chlorite and illite, dispersion accounts completely for permeability damage to porous media. However, unlike the other clays, montmorillonite particles can imbibe water and adsorb ions between layers of sub-particles, or platelets. These platelets have net negative charges on their faces and are held together by exchangeable (or removable) cations such as Na and Ca decrease in ion concentration (salinity) in the fluid surrounding a particle causes migration of water into the clay layers and disperses the basic particle, while diffusion removes the original exchangeable ions from between the platelets. Once these ions are removed, the facing negative platelets repel each other, causing the montmorillonite to swell until, for all practical purposes, the individual platelets are dispersed. For this reason, fresh-water* damage is much more severe in sands containing montmorillonite than it is in sands containing other clays. Many investigators have shown that even trace amounts of montmorillonite can be responsible for marked reduction in the permeability of reservoir sands in the presence of fresh water." ." Monaghan and others have shown that fresh-water damage in montmorillonite-containing cores cannot be

Jan 1, 1965

By W. L. Martin, J. N. Dew, H. B. Steves, M. L. Powers

This paper presents and discusses the results obtained during a pilot test in the Loco field in southern Okla homa. The test was conducted in a 2%-acre pattern that was part of a 20-acre conventional waterflood pilot area. The conventional flood was well past its economic limit when the hot waterflood was initiated to obtain technical information and operating experience. Temperature data from nine observation wells showed that about 75 percent of the pattern was affected by heat, and that heat losses were severe in the pilot pattern. About 60 percent of the injected heat was lost to overburden and underburden zones during hot-water injection. Wellbore heat losses were held at tolerable levels at the shallow depth of the test by providing a low-pressure air annulus between the injection tubing and well casing. Hot water provided water injectivity increases of 200 to 400 percent. The hot water channeled across the lower portion of the oil sand in three directions through zones of relatively high water saturation. There was no conclusive evidence that natural fractures or pressure partings aflected the flow of fluids in the pilot pattern; however response undoubtedly was aflected by localized pressure gradients and by injection-production rate ratios. The results showed that hot waterflooding increased oil recovery in a reservoir containing 600 cp oil. The total tertiary oil production from the pilot pattern area was 3,896 bbl or about 156 bbl/acre-ft swept by heat. The corre.sponding WOR was about 34:I. Description of the Hot Waterflood Process To hot waterflood an oil reservoir, water is injected that has been heated to a temperature substantially higher than the original reservoir temperature, but lower than the vaporization temperature of water at the prevailing pressures. In the reservoir the hot water flows continuously into cooler sand and rapidly loses heat to the sand until it has been cooled to the original reservoir temperature. Thus, a heated zone and a region or "bank" of cooled water begins to accumulate immediately after hot water injection is started. This bank of cooled water continues to grow ahead of the heated zone, which also grows, but at a slower rate. This occurs because heat transfer is almost instantaneous, and the ratio of heat capacities of water to rock is such that two or three unit PV of hot water must be injected to heat a given unit bulk volume of the reservoir. The primary displacement mechanism is the same for both hot and conventional cold waterfloods; i.e., "piston" displacement occurs at the original reservoir temperature. The incremental benefits of hot waterflooding usually occur long after the breakthrough of cold water at producing wells, and the increased oil recovery necessarily is accomplished with high WOR's (water-oil ratios). Heat decreases the viscosity and density of oil and water. These effects result in more rapid and increased recovery of secondary or tertiary oil. If the cost of the required heat is low enough, the ultimate oil recovery of a hot waterflood should be increased substantially over what would be expected at the economic limit of a conventional cold waterflood. The economic benefits of any hot-fluid injection project depend primarily upon the cost of the heat required to produce more oil at an increased rate. This cost depends in part upon the amount of heat lost to surrounding formations. Heat loss depends upon reservoir thickness, water injection rate and temperature, the depth of the formation, well spacing, and the characteristics of the reservoir rock and surrounding formations. In general, percentage heat losses decrease as injection rate and reservoir thickness increase. Although it is an old idea, hot-water injection has not received widespread field application in the oil industry as a drive process. Much of the original field work with hot water was done in Pennsylvania fields where water permeabilities and injection rates are low. In these cases, hot-water injection was used primarily as a means of increasing in-jectivities — not as a recovery process. Ramey' recently has published an excellent review of the development of hot fluid injection processes. Recent publications"' and our own experiences now indicate that hot water or steam injection also can effectively increase the oil recovery from reservoirs containing viscous crudes. However, the economics of these methods as displacement processes have not been established. Several theoretical predictive techniques have been published (for a review see Ramey'), but simplifying assumptions make

Jan 1, 1969

By W. G. Bearden, P. P. Scott, G. C. Howard

This paper concerns the rupture or breakdown of rock formations as related to drilling, completing, and stimulating production of wells, and comprises data compiled from a study of literature and records of treatment of oil and gas wells. and from tests conducted in bores drilled into rock cores and outcrops of rock. Results of the investigation indicate that the internal pressure to rupture cylinders of rock and to breakdown rock formations surrounding a bore in the earth is dependent upon the extent of intrusion of fluids. the position of bedding planes. the ratio of internal to external diameter, the tensile strength of rock, and magnitude of confining pressure, and is independent of the size of bore, degree of fluid saturation. and temperature of rock within practical limits. It is concluded that the mathematical relationship of pressure in bores and stresses in the surrounding rock must not be limited by the simplifying assumptions of homogeneity. isotropy, and impermeability; that the incidence of lost circulation of drilling fluids to induced fractures may be reduced by preventing intrusion of fluids into the small intrinsic fractures along weak bedding planes; and that the magnitude of the breakdown pressure of wells to be treated may be lowered by removal of mud cake. INTRODUCTION The purpose of this paper is to present and discuss the results of tests which may serve to broaden the understanding of the phenomena of rupture or breakdown of rock and thus contribute to the improvement of the techniques for drilling and completing wells. including such operations as preventing lost circulation, stimulating production, and placing cement. Since the early recognition of the possibility of rupturing rock adjacent to a well by fluid pressure, the distribution and magnitude of stresses around a well and the internal pressure to cause failure have been expressed mathematically by application of the principles defining the elastic and inelastic behavior of thick-walled cylinders of a homogeneous. isotropic, impermeable material. In either the elastic or plastic state, if the conditions of homogeneity. isotropy, and impermeability were the normal characteristics of rock. there would be no reason to question the validity of the above principles when applied to rock. However. since most rock formations are characteristically permeable to some degree, contain bedding planes and small intrinsic fractures. and are heterogeneous, rendering invalid the assumptions of homogeneity, isotropy, and impermeability, the application of the thick-walled cylinder theories to the behavior of rock formations penetrated by a bore appears illogical. Also the failure to consider the effect of the properties of the drilling fluids, which certainly intrude in various degrees into most rock formations. is further cause for questioning the application of these principles. Even in the case of cylinders of hardened chrome-nickel steel, which would he considered relatively homogeneous isotropic. and impermeable, as compared to rock, it has been observed that the internal pressure capable of being withstood by the cylinders varied with the sizes of molecules comprising the rupturing fluid used. Because of these apparent weaknesses in the theories of the behavior of thick-walled cylinders when applied to bores in rock, and the importance of knowing the true behavior of rock under the influence of fluid pressure when planning drilling and completing procedures for wells, tests were undertaken to determine by observation the actual rupturing pressure of cylinders of rock and the effects of such variables as environmental conditions and characteristics of rock formations, fluid properties, and bore dimensions. It was intended that by these tests the validity of the thick-walled cylinder theories when applied to rock would be determined and, if found invalid, mathematical expressions would be derived for predicting stress and rupturing pressure of rock formations under various conditions. Since the derivation of mathematical expressions has been only partially completed, all of the results of the latter phase are not included in this paper. PROCEDURE The investigation was initialed with a study of the theories pertaining to the rupturing of thick-walled cylinders of homogeneous. isotropic, impermeable material and a study of the

Jan 1, 1953

By H. H. Rachford, J. Douglas, D. W. Peaceman

A numerical solution of equations describing two-phase flow in porous media shows promise in providing a technique for predicting the displacermet from satlds of oil by water or gas. The description includes the influence of relative permeability, fluid viscosities and densities, graviry, and capillary pressure, and, tholrgh tested ortly for two-dimensional cases, should be eqlmlly applicable in three-dimensional geometry. Two memorical methods ore presented: the first method is quite general in its applicability; the second method can be used only with certain types of boundary conditions but requires less computing. Comparisons of computed results with data from Inhoratory models are presented. These data were taken on a water flood of a strarified model and on water floods of a five-spot model for favorable and unfavorable mobility ratios. On the stratified model, excellent agreement with recovery at breakthrough was obtained; agreement with recovery after breakthrough was poor. In the five-spot model, good agreement was obtainecl with recovery at and after breakthrough. The purpose of this paper is twofold: first, to present a reservoir engineering method requiring knowledge only of rock geometry and the normally measured rock and fluid properties for calculating the multi-dimensional flow of water displacing oil from porous, water-wet rock containing connate water; and, second, to investigate the validity of the method by comparing results of calculations with previously observed displacements in laboratory models. The availability of such a verified technique of reservoir analysis would afford major advantages. First, it would demonstrate that the fundamental macroscopic concepts of two-phase fluid mechanics, i.e., relative permeability and capillary pressure, yield an adequate description of the physical process. Also, since most reservoir rocks appear to be water-wet and contain connate water along with oil, it would provide a method of technical and economic value for calculating the course of oil displacement by water directly from measured reservoir and fluid properties. It is important to emphasize that the calculation is subject only to the limitations of detail and precision of reservoir information, but not to limitations introduced by simplifying assumptions. Further, it would providc almost the only means of examining the influence of factors such as the size and extent of reservoir inhomogeneities and uncertainties in the basic reservoir data. Knowledge of this influence is essential in stating quantitatively the detail required to define a reservoir and in establishing the relation between the uncertainty in the definition of reservoir properties and the reliability of predicted performance. Moreover, in the process of solving a particular problem much detailed information would become available about the displacement process. For example, at each point in the reservoir for all stages of the displace ment. the calculations would yield not only the water and oil saturations but also the direction and magnitude of fluid velocities and the local fluid mobilities. Such dctail is potentially of great value in giving insight into the mechanics of particular displacements. The development of such a method has long been the goal of research in the application of numerical an-alysis to petroleum reservoir enginecrinz. The most complete treatment of the displacement process published to date is that of Douglas. Blair, and Wagner,' which, however, was limited to flow in a single dimension. The method presented in the present work is hascd on the numerical solution of a finite difference analogue of the multi-dimensional differential system describing the displacement process. Although current work has considered only displacement of oil by water from water-wet sands, the differential system for other immiscible displacements, such as gas displacing oil with which it is in phase equilibrium, is quite similar. It should therefore be expected that the technique described will be applicable to displacement by gas as well. The description of the simultaneous flow of fluids through porous media in terms of relative perrncability and capillary pressure has been adequately discussed in the literature and standard textbooks. Sec, for exampje, Muskat. Chapter VII.'

By J. D. Rice, S. C. Pitzer, C. E. Thomas

Interpretation of pressure build-up data obtained in the conventional manner has often been difficult because of the deviation from theoretical behavior. Major causes of this deviation have been attributed to damage and afterflow, and to fluid redistribution in the wellbore which, in extreme cases, can result in a pressure hump in the early portion of the build-up curve. Theoretical investigations show that bottom-hole pressure is definitely influenced by phase redistribution in the tubing column during surface shut-in tests, and that the magnitude of the effect is sensitive to the producing gas-oil ratio and stabilized rate of flow in the well. For field experiments a wire-line tubing packer, which can be run in the tubing against a stabilized flow rate, was developed for bottom-hole shut-in tests. By use of this bottom-hole .shut-in method, pressure humps previously observed in surface shut-in tests were completely eliminated and the effects of afterflow minimized. Analog and digital computer studies have been made to obtain theoretical curves for comparison with field curves, and remarkable agreement between the results of bottom-hole shut-in tests and theoretical curves has been obtained. INTRODUCTION Pressure build-up data from shut-in wells have been used by the petroleum industry to determine the permeability of the formation, to estimate wellbore damage, to evaluate the static reservoir pressure and to speculate reservoir volume. Calculation of these factors is based on methods of analysis developed for theoretical systems whose build-up curves have a charac- teristic shape when the wells are shut-in at the sand faCe.1,2,3,5,8,10 However, field curves obtained by conventional surface shut-in methods do not always exhibit this characteristic. Such factors as stratification, rock heterogeneities and irregular reservoir geometry can cause the character of a build-up curve to deviate from that predicted by theory for a simple system. In addition, all field build-up data are affected by the methods utilized in obtaining the measurements. Conventional surface shut-in techniques allow two effects to occur which contribute directly to the manner in which pressure builds up at the sand face. The first, afterflow, has been recognized for some time and methods are available for estimating the magnitude of its effect and for correcting its influence.' The second, phase segregation in the tubing column after shut-in, has been reported only recently and appears to exert a considerable influence on the character of field build-up curves.!' In extreme cases, phase segregation can produce "pressure humps" in the early portion of the build-up curve. Such curves have been considered anomalous and have defied analysis. It is the purpose of this paper to discuss the effect of phase segregation in the tubing string during build-up and to present a bottom-hole shut-in technique for obtaining field build-up data which minimizes the influence of afterflow and eliminates phase segregation effects. It will be shown that reliable reservoir information can be calculated from build-up measurements obtained using this method, even though conventional surface shut-in tests on the same wells yield anomalous data. DISCUSSION OF PHASE REDISTRIBUTION EFFECTS Some of the effects of phase redistribution on pressure build-up curves have been described previously by Matthews and Stegemeier" who have presented evidence that phase segregation is responsible for the pres-