Search Documents

By M. Prats

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

Jan 1, 1970

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

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

Jan 1, 1969

By A. F. van Everdingen

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

Jan 1, 1953

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

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

By C. Chu

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

Jan 1, 1965

By I. Fatt, D. N. Saraf

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

By A. M. Sarem

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

Jan 1, 1967

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

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

Jan 1, 1970

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

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

Jan 1, 1965

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

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

Jan 1, 1969

By A. G. Spillette, R. L. Nielsen

The purpose of this paper is to further the understanding of reservoir response to hot-water injection by desuribing a two-dimensional, mathematical model of the process. Key assumptions are that no gas phase is present, and that the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. The resulting system of three partial differential equations is solved simultaneously through the use of a "leap-frog" application of standard alternating direction implicit methods for the solution of the mass-balance equations and the method of characteristics for solution of the energy-balance equation. The utility of the mathematical model is demonstrated by comparing numerical and analytic temperature distributions for hot-water bank injection and by comparing calculated with observed field behavior. Additional calculations show that hot waterfioding can recover significantly more oil than cold waterflooding, and that a hot-water bank recovers, with less energy input, nearly as much oil us continuous injection. Introduction Consistent field success with hot fluid injection processes requires good reservoir description and a thorough understanding of recovery mechanisms. The latter is fostered by a combination of well designed experimental studies and physically sound mathematical modeling. The purpose of this paper is to further the understanding of reservoir response to hot-water injection by describing a two-dimensional mathematical technique, indicating its validity, and demonstrating its utility for studying the effects of reservoir and operating parameters. Published experimental studies of hot fluid injection processes are few." Even if extensive data were available, two considerations discourage exclusive reliance on laboratory or field work. First, because of the severity of scaling requirements," laboratory results must be interpreted with care. Second, because of the one-shot nature of field experiments — injectivity/productivity tests and pilot floods—results seldom are available over the desired range of operating conditions. These factors emphasize the need for devoting attention to mathematical modeling to complement laboratory and field work. Through the use of mathematical models, scaling uncertainties can be bridged and response can be predicted for unique combinations of reservoir and operating conditions. Numerous mathematical models developed during recent years enable us to calculate temperature distributions, thermal efficiency (or conversely, fraction of heat lost), and oil recovery behavior for hot fluid injection. Ramey" and Spillette' recently have provided reviews of these methods. Most models have concentrated on predicting temperature distributions and thermal efficiency; few have been directed at predicting oil recovery.'."1° Although these techniques are useful for estimating effects of reservoir parameters and operating conditions on process performance, results are limited by the assumptions made and by the methods of coupling independently solved fluid-flow and energy-balance equations. A computer-based method is presented for predicting total reservoir response to hot-water injection that obviates most simplifying assumptions. The model simulates fluid flow and heat transfer in two dimensions within a vertical cross-section spanning the oil sand and adjacent unproductive strata. Known numerical procedures are used to solve the governing partial differential equations. The mathematical model handles the effects of reservoir heterogeneity, gravity, capillarity, relative permeability and temperature-dependent fluid properties. In addition, a wide range of operating conditions can be modeled, including hot-water followed by cold-water injection. Mathematical Description of Hot-Water Injection Key assumptions in the mathematical model for hot-water injection are that (1) no gas phase is present and (2) the injected fluid reaches thermal equilibrium instantaneously with the reservoir fluids and sand. Relative permeability and capillary pressure are assumed independent of temperature. Darcy's law for multiphase fluid flow relates the velocities of the oil and water phases to the phase pressure gradients by the following relations:

Jan 1, 1969

By H. J. Morel-Seytoux

The influence of patten, geometry on assisted oil recovery for a particular dispIacement mechanism is the object of investigation in this paper. The displacement is assumed to be of unit mobility ratio and piston-like. Fluids are assumed incompressible and gravity and capillary effects are neglected. With these assumptions it is possible to calculate by analytical methods the quantities of interest to the reservoir engineer for a great variety of patterns. Specifically, this paper presents (1) very brieily, the methods and mathematical derivations required to obtain the results of engineering concern, and (2) typical results in the dorm of graphs or formulae that can be used readily without prior study of the methods. Results of this work provide checks for solutions obtained from programmed numerical techniques. They also reveal the effect of pattern geometry and, even though the assumptions of piston- like displacement and of unit mobility ratio are restrictive, they can nevertheless be used for rather crude but quick, cheap estimates. ,These estimates can be refined to account for non-unit mobility ratio and two-phase flow by correlating analytical results in the case M = 1 and the numerical results for non-piston, non-unit mobility ratio displacements. In an earlier paper1 it was also shown that from the knowledge of closed form solutions for unit mobility ratio, quantities called "scale factors" could be readily calculated, increasing considerably the flexibility of the numerical techniques. Many new closed farm solutions are given in this paper. INTRODUCTION BACKGROUND Pattern geometry is a major factor in making water-fIood recovery predictions. For this reason many numerical schemes have been devised to predict oil recovery in either regular patterns or arbitrary configurations. The numerical solutions, based on the method of finite difference approximation, are subject to errors often difficult to evaluate. An estimate of the error is possible by comparison with exact solutions. To provide a variety of checks on numerical solutions, a thorough study of the unit mobility ratio displacement process was undertaken. To calculate quantities of interest to the reservoir engineer (oil recovery, sweep efficiency, etc.), it is necessary to first know the pressure distribution in the pattern. Then analytical procedures are used to calculate, in order of increasing difficulty: injectivity, breakthrough areal sweep efficiency, normalized oil recovery and water-oil ratio as a function of normalized PV injected.

Jan 1, 1967

By A. T. Chatas

A description of the geometrical characteristics of spherical reservoir systems, a discussion of unsteady-state flow of such systems and examples of engineering applications are presented as backgmund material. The fundamental differential equation, a description of average spherical permeability and the introduction of the Laplace transformation serve as theoretical foundations. Engineering concepts are investigated to indicate particular solutions of interest, which are analytically obtained with the aid of the Laplace transform. These are numerically evaluated by computer, and presented in tabular form. INTRODUCTION A tractable mathematical analysis of unsteady fluid flow through porous media generally requires incorporation of a geometrical symmetry. The simplest forms include the linear, cylindrical (radial) and spherical. Most analytical endeavors have concentrated on cylindrical symmetry because it occurs more often in petroleum reservoirs. Nevertheless, some reservoir systems do exist that are better approximated by spherical geometry. Review of technical literature revealed but a single reference to unsteady spherical flow in petroleum reservoirs.l The motive and purpose of the present work was to remove this gap in technical information, and to provide the practicing engineer with some useful analytical tools. The mathematical details associated with the particular solutions of interest involved use of the Laplace transformation. Hurst and van Everdingen previously demonstrated the efficacy of this operational technique, and in many respects the present treatment was patterned after their earlier work.2 PRELIMINARY CONSIDERATIONS GEOMETRICAL CHARACTERISTICS Geometrically, a spherical reservoir system is defined at any instant of time by two concentric hemispheres whose physical properties of interest vary only with the radial distance. Every physical property is thus restricted to be a space function of only one variable: the distance along a radius vector emanating from the center. Such a system is composed of an outer region and an inner region, separated by a defined internal boundary. The inner region simply extends inward from this boundary, whereas the outer region extends outward from it to an external boundary. The position of the internal boundary is presumed fixed, so that the size of the inner region remains constant. On the other hand, the position of the external boundary at any given instant of time is determined by the distance into the system that a sensible pressure reaction has occurred, Thus, the external boundary may change position with time. It initially emerges from the inner region and advances outward to its ultimate position. When this ultimate position coincides with a geometric limit, the reservoir system is said to be limited. When it coincides with points subject to pressure gradients furthest removed from the internal boundary, yet short of a geometric limit, the system is said to be unlimited. In this investigation two different boundary conditions are imposed at the ultimate boundaries of limited systems. The first requires that no fluid flow occur across this boundary; the second that the pressure remain fixed at this boundary.3-5 UNSTEADY-STATE FLOW In a strict sense virtually all flow phenomena associated with a reservoir system are unsteady-state. The transient behavior of these phenomena requires accounting, however, only when time must be introduced as an explicit variable. Otherwise, steady-state mechanics may be used. Analytically, steady-state conditions prevail in a reservoir system only over that portion of its history when this relation is satisfied:

Jan 1, 1967

By G. E. Warner

This paper presents a review and analysis of a highly stratified Burbank sand waterflooding project in Osage County, Okla. Permeability values in this reservoir range from less than 0.1 md to nearly 3 darcys, and 90 percent of the permeability capacity is contained in approximately 26 percent of the reservoir volume. To date, the project has recovered 103 bbl/flood-acre-ft since waterflood initiation, which represents 50 percent of the total preflood cumulative and 12 percent of the original stock-tank oil in place (OSTOIP). The outstanding feature of this project has been the economical recovery of a sizable quantity of oil at water-oil ratios higher than those normally considered prohibitive. Nearly 70 percent of the recovery to date has been achieved at water-cut values greater than 95 percent. An ultimate waterflood recovery of 138 bbl/flood-acre-ft, or 67 percent of the preflood recovery and 16 percent of the OSTOIP, is predicted. The ultimate primary and secondary recovery will be 302 bbl/acre-ft and approximutely 40 percent of the OSTOIP. The actual performance is compared with the predicted performance by both the Stiles and the Dykstra-Parsons methods, and applicabiliry of these methods to this particular reservoir is evaluated. The performance of this project indicafes highly stratified reservoirs can be successfully and economically waterflooded when either reservoir or economic conditions preclude the use of mechanical and chemical means of mdbility control. Introduction The concept of the stratified reservoir has been used by reservoir engineers for many years to describe the vertical variations in permeability within a reservoir. Numerous techniques have been advanced to predict the expected waterflood performance of the stratified reservoir, and in most cases they prove reasonably accurate. The advisability of undertaking a prospective waterflooding project is determined, based on these predictions. This paper presents a review of a waterflood project in a highly stratified reservoir that normally would be considered a marginal prospect, at best, when evaluated with any prediction technique properly adjusted for reservoir conditions. This kind of project can be successful if the predictions are used to design a waterflooding system capable of economically recovering an oil cut of 1 to 2 percent. Experience gained in the early stages of the project's development was profitably applied in the later development of the entire reservoir and could prove beneficial to other projects with similar reservoir stratification or fracture characteristics. The East Burbank pool waterflooding operations have been developed in three separate stages, each of which has been a separate project: Mid-Burbank Unit, Flood A, and Stanley Waterflood. When referred to collectively, these projects will be called the Composite Project. The Skelly Oil Co.'s Barber lease — the only portion of the pool not included in the three projects — is operated as a separate facility with a cooperative injection-well agreement. It contains less than 80 productive acres and would not significantly affect the over-all performance picture.

Jan 1, 1969

By C. Cronquist

Since April, 1962, Shell Oil Co. has operated a peripheral line-drive waterflood of the 10.750-ft lower Tuscaloosa (Cretaceous) Denkman sand in the Little Creek field. Located in southwestern Mississippi in the lower Tuscaloosa trend, the field was discovered by Shell in Jan., 1958. Among the largest in the trend, the accumulation is controlled stratigraphically and occurs in an extensive series of alluvial-point bar deposits across a nose that dips gently to the south. The crude oil is 39" API gravity and was highly undersaturated at initial conditions. Early production behavior indicated a depletion drive with slight water influx; primary recovery wm estimated to be 25 million bbl, or 24.5 percent of stock-tank oil originally in place. A.s a result of a favorable mobility ratio and remarkably uniform rock properties, the volumetric sweep eficiency in the waterflood has exceeded 90 percent. Despite the high connate water saturation (56 percent), the displacement efficiency by water injection has been quite efficient, contrary to traditional concepts. Based on observed field flood-out performance, ultimate recovery is calculated to be 46 million bbl. As anticipated, the production rate began declining during 1964 due to a continu ing reducrion in the number of producing wells as flood fronts advanced across the field. Cumulative recovery as of Jan., 1968, was 44.8 million bbl; waterflood operations ~hould be complete by 1970. Introduction Most waterfloods described in reservoir engineering literature involve pattern operations, the most common of which is the five-spot pattern. Only occasiona'lly are linear drives described, even though there probably are increasing numbers being put into operation. Increasing use of linear floods may reflect the accelerating trend towards early application of fluid-injection programs when linear drives are more efficient, as was the case at Little Creek. However. the Little Creek operation is unusual in several respects. It is one of a few successful waterfloods conducted in sands with connate water saturations approaching 60 percent, and it is among the deepest successful water-floods being operated anywhere. Development The Little Creek field was discovered in Jan., 1958, in the lower Tuscaloosa trend in southwestern Mississippi (Fig. 1). The discovery well, Shell-Lemann No. 1, was drilled on a seismic closure about 2 miles southeast of lower Tuscaloosa production at Sweetwater. The well was completed in the Denkman sand at 10,752 ft, and on initial potential it flowed 588 BOPD of 39" API crude on a %+-in. choke; psi was 730 and GOR was 442:l. Development proceeded rapidly on 40-acre spacing, and by the end of 1958 the field was producing approximately 9,100 B/D from 56 wells. Almost all this production came from the northern half of the field. The southern part was discovered in Nov., 1958, at a location almost 2 miles from previously established production. Subsequent drilling established an apparent conneotion about half a mile wide between the twd areas. Development of the field was practically complete by 1961 with the drilling of 155 producers and 35 dry holes. Some early wells were completed open hole but moslt are completed with casing cemented through the productive interval and perforated for production. Geology The lower Tuscaloosa structure in this area is a north-south nose with maximum dips of about 2" on the flanks. There is no discernible faulting at the producing horizon. Productive limits of the field are shown in Fig. 2, which is contoured on a marker just above the main productive zone. A common water level at 10,415 ft subsea was logged at both the northwest and southeast ends of the field. There is a gross oil column of nearly 100 ft (Fig. 2). Since there is less than 40 ft of closure, it is apparent that bhe accumubtion is controlled stratigraphically, being limited by a sand pinch-out across the creslt of the nose. The total produc-

Jan 1, 1969

By B. B. Wood, M. J. Rzasa, E. W. Hough

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By E. W. Hough, B. B. Wood, M. J. Rzasa

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By B. H. Caudle, W. J. Bernard

Factors which influence the success or failure of a waterflood can seldom be determined in the laboratory. For this reason pilot waterfloods are initiated in a repreventative portion of the oil reservoir in question. For a pilot flood to predict quantitatively the recovery to be expected in a field-wide waterflood operation, the pilot area must behave as though it were confined (surrounded by similar areas). In this study, laboratory fluid-flow models were used to determine the simplest pilot pattern, for particular conditions of mobility ratio and initial gas saturation, that would behave as though it were confined. Pilot patterns studied ranged in complexity from a single inverted five-spot to a grouping of nine regular five-spots. Only the innermost producing well in each pattern was studied. Model results showed that the optimum number of wells in the pattern depends upon the oil-water mobility ratio and the expected oil-bank size. Unfavorable mobility ratios will, in general, require more wells in the pilot pattern than will favorable mobility ratios. Pilot patterns in reservoirs which contain a dispersed, flowable, free gas saturation will require fewer wells than for the under-saturated case. The single inverted five-spot pattern was found to be unsatisfactory for predicting behavior of fully developed waterfloods. In particular, it is possible that, in reservoirs which contain a flowable, dispersed gas phase, the oil bank will never be observed at the producers due to the large amounts of free gas which continue to be produced with the oil. INTRODUCTION One method which has been used to predict the performance of a waterflood is the pilot flood. The pilot waterflood is a flood which involves only a small cluster of the reservoir wells and is located in a small, representative portion of the reservoir. The object is that oil produced from the pilot can, in some way, be related to the oil recovery to be expected from a field-wide expansion of the waterflood. However, these field pilot waterfloods have often been unreliable in the prediction of oil recovery in a fully developed waterflood. This unreliability has also been demonstrated in several laboratory studies of pilot floods. Some of the investigators have shown that there are situations in which the pilot flood oil production is far too optimistic with respect to the oil recovery in the fully de- veloped flood. Others4-G have shown that the pilot results can also be pessimistic, especially if the pilot waterflood is initiated in an oil reservoir which has been depleted by primary recovery and is at very low pressure. The major reason for this unreliability of pilot water-floods is the migration of fluids into or out of the pilot area. By the well-known method of images, if straight lines can be drawn to represent vertical planes of symmetry in a porous medium which contains pressure sources and sinks (injectors and producers), then these lines are invariant streamlines, or lines across which there is no potential gradient, and therefore no flow. In an actual reservoir, these lines of symmetry can never be established exactly because of reservoir inhomogeneities and irregular reservoir boundaries. However, if the reservoir is relatively large and contains wells in repetitive patterns, these lines of symmetry are commonly assumed to exist for the pattern units sufficiently far removed from the reservoir boundary. Lines of symmetry for the five-spot injection pattern are shown in Fig. 1. Each five-spot unit in this figure can be considered confined with respect to flow across its boundary. In pilot floods this is not the case. The lines of symmetry for the pilot patterns investigated in this study are shown in Fig. 2. It is obvious that the fluid within these pilots is not confined and is therefore able to migrate into or out of the pilot area. Intuitively, one can see that, if more wells are added to the pilot, the innermost unit tends to behave more and more like the confined pattern. However, there is a practical limit to the number of wells which should be placed in the pilot. This limit is usually determined by economic factors. It was the purpose of this study to use laboratory fluid-flow models to determine which of the previously mentioned pilot patterns will force the innermost producing well to behave as it would in a fully developed waterflood. Since fluid migration is influenced by initial saturation conditions and the mobility ratio, these factors were included in the study. The ultimate objective of this study was to develop data which would allow the operator to choose a pilot pattern and operating conditions that will yield a production history which can be applied directly as an estimate of the performance of each production well of the fully developed waterflood. BASIS FOR THE STUDY The basic problem of field pilot floods is the migration of the reservoir fluids into or out of the pilot area. This problem has been the subject of previously reported model studies on pilot floods. These studies have been concerned mainly with the development of arbitrary "correction factors" to be applied to the simple, unconfined pilot systems such as the single five-spot. The correction factors were intended to adjust the production history of the un-

By John S. Bell, J. M. Shepherd

pressure difference of 280 psi measured initially on opposite sides of the major fault in the Hawkins Field led to pressure determinations in the Woodbine sand throughout the East Texas basin. Considering the large area, there is unusually good correlation between pressure losses and points of major withdrawals, these losses indicating that the basin pressure has been materially reduced by withdrawals, primarily from the East Texas Field and secondarily from the old fault-line fields of Limestone and Navarro counties, with the pressure loss varying inversely with the log of the distance from these two main areas of basin withdrawals to the measurement points. Apparently faults located throughout the basin serve as local blocks to continuous pressure gradients. INTRODUCTION The pressure anomalies found in the Hawkins Field, Wood County, Texas,' shortly after discovery have led to an area-wide investigation of the pressure distribution within the Woodbine sand of the East Texas basin. These pressures were obtained by drill-stem tests made for that sole purpose while drilling wildcat wells to deeper horizons. The Woodbine sand of Gulf Cretaceous Age is or has been productive of oil and gas in 42 different fields, with the 28 oil fields currently active having produced 118,000,000 bbl of oil during 1949. With one exception, producing energy for all Woodbine reservoirs has been from water drive provided predominantly by expansion of the tremendous volume of water contained in the Woodbine sand. This producing mechanism was recognized at an early date and was described mathematically by

Jan 1, 1951

By John S. Bell, J. M. Shepherd

pressure difference of 280 psi measured initially on opposite sides of the major fault in the Hawkins Field led to pressure determinations in the Woodbine sand throughout the East Texas basin. Considering the large area, there is unusually good correlation between pressure losses and points of major withdrawals, these losses indicating that the basin pressure has been materially reduced by withdrawals, primarily from the East Texas Field and secondarily from the old fault-line fields of Limestone and Navarro counties, with the pressure loss varying inversely with the log of the distance from these two main areas of basin withdrawals to the measurement points. Apparently faults located throughout the basin serve as local blocks to continuous pressure gradients. INTRODUCTION The pressure anomalies found in the Hawkins Field, Wood County, Texas,' shortly after discovery have led to an area-wide investigation of the pressure distribution within the Woodbine sand of the East Texas basin. These pressures were obtained by drill-stem tests made for that sole purpose while drilling wildcat wells to deeper horizons. The Woodbine sand of Gulf Cretaceous Age is or has been productive of oil and gas in 42 different fields, with the 28 oil fields currently active having produced 118,000,000 bbl of oil during 1949. With one exception, producing energy for all Woodbine reservoirs has been from water drive provided predominantly by expansion of the tremendous volume of water contained in the Woodbine sand. This producing mechanism was recognized at an early date and was described mathematically by

Jan 1, 1951