Search Documents

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

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

Jan 1, 1969

By A. F. van Everdingen

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

Jan 1, 1953

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

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

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

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

Jan 1, 1965

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

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

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

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

By C. Chu

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

Jan 1, 1965

By I. Fatt, D. N. Saraf

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

By A. M. Sarem

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

Jan 1, 1967

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

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

Jan 1, 1970

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

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

Jan 1, 1965

By K. K. Clark

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

Jan 1, 1969

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

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

Jan 1, 1969

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

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

Jan 1, 1969

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&apos; 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.&apos;."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&apos;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 R. W. Parsons

wben a displacement front encounters an isolated permeability heterogeneity, it will be perturbed. The details of this perturbation will depend on the beterogeneity size and shape, the Permeability contrast, and the mobility ratio. This paper deals with the displacement front behavior for one, simple, mathematically tractable model: a circular permeability discontinuity and unit mobility ratio. The progress of an initially planar front is shown for four permeability contrasts: k2/k1 = 0, 5/2, 2, m. The front distortion is more pronounced than streamline and isopotential distortion and it approaches a constant shape downstream from the discontinuity. INTRODUCTION The impact of reservoir heterogeneities on oil displacing processes depends upon the variation of the heterogeneous property (frequency distribution), its disposition (spatial distribution), and the inherent stability of the displacement mechanism.l The interplay of these factors plus the impracticability, if not impossibility, of describing a natural rock system in any detail makes exact displacement front movements difficult to determine. It is useful, however, to be able to qualitatively visualize what heterogeneities will do to a displacement front. In this respect certain simple models that can be treated analytically are of value. The following analysis will show the displacement front behavior for one such model. MODEL DESCRIPTION The flow field is two-dimensional and infinite in extent. The fluid flow is single phase, steady state and unidirectional at x, y = 2 m, or before insertion of the discontinuity. A circular permea- bility discontinuity is placed at the origin. Porosity and thickness are the same within and exterior to the heterogeneity. This is exactly the same model that was used by Greenkorn et al. to check analytically Streamlines as determined by Heie-Shaw flow models and finite difference computer techniques.2 FRONT POSITION CALCULATIONS The results of displacement front calculations are presented in Figs. 1 through 4 for four different permeability contrasts: k2/kl = O, %, 2, m. The times corresponding to the front locations are the Same for the 0 and cases! and for the 1/2 and 2 cases. The determination of front positions at various times is based on known analytical expressions for velocity potentials ().2 and stream functions ().2 These expressions are in terms of the x and y

Jan 1, 1969

By G. M. Gay, R. L. Perrine

This paper describes a method of numerical computation for three-dimensional, unstable, miscible displacement behavior useful for heterogeneous systems, as well as for more ideal conditions. In the method, flow equations are first linearized by a perturbation approach. The basic flow process is separated and a solution for its behavior readily obtained. The remaining problem of deviations from the basic flow caused by non-ideal conditions is then subjected to numerical analysis. Results obtained from use of the method are also presented. Although conditions assumed in the test cal-culations were severe, results show the type of dispersing flow expected and appear quite satisfactory. The method has eliminated or reduced in importance problems of oscillating values near steep fronts, excessive computer smoothing, etc. A unique advantage of the method is that the source of variations from ideal behavior can be observed. ,The one serious drawback results from the algebraic complexity of the perturbation approach, and the need for second-order terms to be retained in calculations of interest. Fewer array points are available and more computer time is required than would be desired. However, these difficulties are also experienced with other approaches to the solution of three-dimensional dispIacement problems. INTRODUCTION Prediction of behavior of the miscible displacement process within a porous medium for any system of engineering importance is plagued by a number of difficulties. With typical fluid properties the displacing fluid has the lower viscosity, and there is a natural tendency toward flow instability. The problem of predicting instability is compounded by the fact that every real system is heterogeneous. Permeability will vary from point to point — not entirely systematically, and yet not in a random fashion leading to a readily defined average. Furthermore, permeability is unlikely to be isotropic. These permeability properties, which are characteristic of any real porous medium, accentuate the effects of flow instability. Another factor to be considered is that flow dispersion accompanies the displacement process in a porous medium. Mechanistically, dispersion is due to the fact that flow between any two points in the medium follows multiple tortuous paths, each characterized by slightly different flow properties. The fact that the coefficient characterizing dispersive properties of the real medium JS a tensor with variable elements complicates matters. Further difficulties arise because the parabolic influence, while not negligible, is small. Thus, there is a tendency toward steep and uneven displacement fronts along which dispersion smoothing must be represented accurately. Yet frequently used numerical analysis schemes may tend toward either instability, oscillation or excessive smoothing, none of which gives the desired accurate picture of flow behavior. Thus, while many experimental and analytical studies of the process have been made, predictions of actual performance are still subject to considerable uncertainty and possible improvement. This paper reports on part of a study which attempts to develop improved methods for solution of the flow equations describing the miscible displacement process. Of necessity, the calculations were performed on a large digital computer. A three-dimensional system is represented, and the coefficients defining dispersion and permeability can be varied in a manner representative of a real system. Because the differential equations are nonlinear, a first step is to linearize them using a perturbation approach. Objective of this approach has been to separate the basic flow which would occur in an ideal system and for which predictions are readily made. Then the numerical analysis, which is necessarily complex, is directed toward solution of the perturbation equations. These equations state how the basic flow is modified

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 F. W. Jessen, N. Mungan

The plastic flow characteristics of clay water suspensions were first recognized by Binghaml in 1916 and further studied by Ambrose and Loomis&apos; in 1931-1932. Many physical and chemical properties of clay suspensions are deeply influenced by particle size distribution as well as by the exchange cation. Separation of a number of fractions is particularly helpful in making a detailed study of a suspension. A limited amount of work has been conducted on the effects of particle size on specific properties of mud drilling fluids. A knowledge of particle size distribution is essential for a better understanding of fundamental properties of clay suspensions, particularly viscosity and gel strength. The purpose of this investigation was to study the effect of particle size on the viscosity and gel properties of clay suspensions and to determine any changes in mineral composition accompanying variation in size. EXPERIMENTAL PROCEDURE Numerous methods may be used to fractionate a poly-dispersed clay suspension into nearly mono-dispersed systems. Particles larger than 44 micron (325-mesh) may be obtained by wet screening. For sizes from 44 to 1 micron, the fractionation may be accomplished by sedimentation, elutriation plus microscopy,&apos; or turbidimetry.6 When dealing with systems which contain a major portion of particles below l-mi-cron size, slit ultramicroscopy,7 centrifuge methods of Svedberg,5 eams,9 and McBain,l0 and X-ray methods" may be employed. The Sharples super centrifuge has been proved suitable for the economical separation of suspended particles of submicron range.12,13,14 Three-hundred and forty-five gal of 1 per cent by weight Aquagel suspension was mixed in a high-speed 5-gallon mud mixer and the suspension was allowed to stand in 5-gal glass bottles for 30 days to assure full hydration. All particles larger than 44 microns were removed by wet screening. Five fractions in the range of 44 to 1 micron were obtained by sedimentation. The Sharples super centrifuge method described by Hauser and Reed.13 Norton and Speil,14 Ha user and Lynn,15 and Hauser and Schachman,10 as employed to yield five fractions of submicron size. A super centrifuge operation chart prepared by Fancher, Oliphant, and Houssierre12 was used to determine operating conditions for the desired particle ranges. The clay suspension, after removal of all material above 1-micron size, was centrifuged. During rotation of the bowl, the particles were deposited on a thin plastic liner which fit closely inside the bowl. After passing 5 gal of suspension through the centrifuge, the liner was removed and the particles on the liner were scraped off, remixed and recentrifuged. This was repeated five times. Thus, Fractions H, I, K, and L were obtained. The finest fraction, M, remained in an approximately 0.4 per cent suspension. Water was evaporated under reduced pressure from this suspension to yield the finest particles. Results of the fractionation are included in Table 1. In order to achieve proper orientation and humidity conditions for X-ray diffraction analysis, the following procedure was employed. The dilute suspensions were placed on glass slides and allowed to dry at room temperature. This was repeated several times. Previous to measurement each sample first was placed in a water desiccator for half an hour and then allowed to stand for 15 minutes at room temperature. thus yielding an equilibrium condition of about 52 per cent relative humidity. The relative humidity of the room in which the X-ray diffraction was carried out was maintained at approximately 50 per cent so that there was no change in the spacings during the X-ray examination. Samples