Search Documents

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 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. 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' 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,' 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

By A. G. Weber, H. J. Welge

A published calculation method for predicting incompressible, multidimensional fluid displacement has been adapted to the problems of water and gas coning in oil wells. Since depth and radial distance from the wellbore are the two key dimensions affecting the shape of a gas or water cone, coning calculations are well suited to the use of two-dimensional methods. The alternating direction implicit procedure (ADIP) was used for relaxation calculations of two-phase potentials in a two-dimensional grid. From the potentials and the capillary pressure relationship, saturation and pressure distributions were calculated which trace cone growth with time. Predictions have been made for well producing histories both before and after core breakthrough. To check validity of the method, two-dimensional calculations have matched the coning behavior and produced water-oil ratio history of a laboratory sand-packed model. They have also matched the coning behavior of several producing wells, for which the calculations were compared with produced water or gas cuts and logs showing water or gas cone movements. INTRODUCTION The application of two-phase, two-dimensional calculations using ADIP to various reservoir flow problems has been described in the literature.1-3 When adapted for computer solution, this method has proved to be a powerful tool for simulating well and reservoir behavior. This paper discusses the method as applied to well coning calculations. Single-well and coning calculations comprise an especially difficult class of two-dimensional problems which require special techniques for computer calculation and determination of reservoir characteristics. Refs. 4 through 8 describe previous approaches to the coning problem. Several examples of water and gas coning calculations, including studies on both laboratory models and producing wells, are presented here. The two-dimensional method accounts realistically for the most critical parameters affecting coning behavior, including production rate, formation stratification, horizontal and vertical permeabilities, depth of well penetration, gravity and capillary forces. The method considers the different densities and viscosities of the two phases and the relative permeability and capillary pressure characteristics of the rock and fluids. In addition to tracing cone growth in the vicinity of the wellbore, the method calculates the overall movement of the fluid interface throughout the well's drainage volume. Incompressible fluid flow is assumed to occur between the producing interval and the well's limit of drainage. Calculations can be made for the producing history both before and after cone breakthrough. A typical two-dimensional grid or array of blocks used to solve a coning problem contains about 400 blocks. The well's cylindrical drainage volume can be represented by about 20 radial subdivisions and the formation thickness by 20 vertical subdivisions. The grid spacing is normally smaller near the withdrawal interval to define the cone shape accurately. For this work we used an IBM 7074 digital computer having a core memory of 10,000 ten-digit words. A typical study, covering 5 to 10 years of well producing history, required from three to six hours of computing time. MATHEMATICAL SIMULATION OF CONING BEHAVIOR BASIC METHOD In coning calculations, the reservoir volume drained by the producing well is represented by a two-dimensional system of blocks as shown in Fig. 1 for water coning studies. The horizontal dimensions of the blocks increase with radial distance from the well axis in geometric progression, i.e., the block size is small near the wellbore and large near the well's drainage radius (re). Vertically, the blocks are bounded by horizontal planes located at different depths through the

Jan 1, 1965

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 J. M. McDowell, E. C. Doty, F. Morgan

Improvements in the X-ray method of measuring liquid saturation and saturation distribution are presented. Two identical direct current amplifiers have been added to measure continuously the intensities of the X-ray tube output and the beam which has passed through the core under test. The ratio of these intensities is automatically recorded on a modified Brown Electronik Recorder. Arrangements have been made to drive the core automatically and in synchronism with the recorder chart so as to give a continuous curve showing the variation of liquid saturation along the length of the core. By determining the ratio of the two X-ray beam intensities simultaneously the measurements are made essentially independent of small fluctuation in voltage and current supplied to the tube. This improved technique has been used for investigating the fluid distribution along the length of cores and core assemblies when one fluid is displaced by another. Typical results of these studies are shown. INTRODUCTION In a previous publication' a method was described for measuring the oil saturation in cores which was based on the absorption of X-rays by the fluid content of the core. While that technique is capable of yielding accurate results, the original apparatus was not especially convenient for determining the fluid distribution along the length of a core. This was primarily due to the slow response of the amplifier and galvanometer and the necessity of making a point by point survey of the material under investigation. It was, therefore, not only possible for sections of the core having a variable or unusual saturation distribution to escape observation entirely, but changes in output of the X-ray generator between readings on the unknown and the monitor introduced errors in the measurements that could be avoided only by further stabilization of the source or by taking a large number of readings. This paper describes a method and apparatus for recording continuously and automatically the intensity of the X-ray beam that is transmitted by a core. The variation of saturation with time at a given section of the core may thus be found, or by scanning the core by driving it in synchronism with the recorder the saturation distribution may readily be determined. As indicated in the previous paper, the X-ray method possesses an advantage over the weighing technique in relative permeability experimentation in that the core need not be removed from the holder for the purpose of making a saturation measurement. Furthermore, the actual distribution of fluid and the effects upon that distribution of inflow and outflow heads may be investigated. And errors due to gas bubbles or loss of liquid that are frequently sources of error in volumetric methods of saturation determination are avoided. PRINCIPLE OF THE METHOD The energy in an X-ray beam that has passed through an absorbing material is expressed by the well known formula: E = E0, e-cd (1) where E, is the energy in the incident beam, µ the mass absorption coefficient, c the density or concentration of absorbing material, and d is the length of path of the beam in the absorber. The coefficient µ varies with the wavelength of the radiation as well as with the absorbing material. The output of an X-ray tube varies in an approximately linear manner with the tube current, but to the second or higher power of the anode voltage. Moreover, the quality or wavelength distribution of the radiation varies with the voltage.

Jan 1, 1950

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 C. K. Eilerts, J. D. Ham

Laboratory research has been conducted to evaluate the characteristic effects of condensate saturation on the mobility of gas in typical reservoir rocks. A pump with two pistons provided phases in equilibrium at controlled flow rates and constant liquid-vapor volume ratios. Core samples were tested in a holder that permitted sealing steady-state phases in the pore spaces under flow test pressures for weighing with sufficient precision to give an accurate measure of the attained saturations. Parameters of the study were pressure, apparent velocity, flowing liquid-vapor volume ratio, fluid composition, core material and length-to-diameter ratio. Tests were conducted at 515 and 1,515 psia to determine the effect of pressure on the mobility-saturation relationship. Influence of the apparent velocities 0.30, 0.15 and 0.07 cm/sec and of five liquid-vapor volume ratios from 0.0001 to 0.01 was determined at 515 psia The principal fluid system was nitrogen and and separator liquid from a Gulf Coast condensate, but the effect of a condensate of different viscosity was also determined. The significance of a possible concentration of liquid near the outlet end of the test cores was investigated. Core materials included a consolidated sandstone and a low-permeability limestone. Relative mobility and liquid-vapor volume ratio relationships are concluded to be dependent on pressure, saturation and, to a lesser extent, on velocity. For each porous medium and fluid there is a minimum saturation essential to two-phase flow, and high velocities of flow have only limited effect on saturations in this range. INTRODUCTION It is now recognized that knowledge of the availability of a reservoir gas — the period that gas may be recovered from a reservoir at a given rate — is as important as the estimate of gas in place. When reservoir conditions are such that pressure decline owing to production will result in retrograde condensation, liquid saturation can increase in the structure around the production well in a way that may restrict flow capacity, and hence future availability. l Based on a numerical solution of transient radial flow of single-phase gas-condensate fluids, 2 a program for computing the accumulation of condensate in a reservoir and the effect of this condensate on deliverability has been developed. During this development a need evolved for data that described the influence of pressure, velocity of flow and condensed liquid on mobility as recovery progressed. This paper presents the results of experimental investigation to obtain insight into the mobility relationship at reservoir pressures of gas-condensate fluids. A core holder for flow tests at reservoir pressure was developed to weigh the core and its fluid contents under pressure to determine saturation. A specially constructed pump was used to drive a two-phase fluid through the core until a steady state prevailed. The core was used as a differential device, and properties measured were intensive rather than extensive. APPARATUS Fig. 1 shows a schematic representation of the major pieces of equipment used in the mobility determination. The weighable core holder is connected to a two-piston, constant-displacement pump. Attainment of steady-state flow is indicated by stability of the flow differential read from the high-pressure manometer. WEIGHABLE CORE HOLDER An important part of the data needed was the mobility-saturation relationship. This required that mobilities and saturations be measured at pressures, flow velocities and liquid-vapor volume ratios characteristic of condensate reservoirs.

By R. J. Sandrea, S. M. Farouq Ali

The results of an experimental and theoretical study of the effects of rectilinear impermeable barriers and highly permeable channels on the sweep efficiency and conductivity of a five-spot network are presented. The study was carried out using Hele-Shaw, and conductive sheet analogs for both. normal and inverted flood patterns. A functional relationship was developed which provides quantitative prediction of the above parameters; i.e., sweep efficiency and conductivity, for the nonhomogeneous system. To evaluate this, it is necessary only to measure the interference modulus of the obstacle using a simple expression. The results indicate that impermeable barriers, with an interference modulus greater than zero, always cause a decrease in the conductivity and the sweep efficiency of the pattern. On the other hand, permeability channels do not have a significant effect on these parameters, unless they are located along the streamlines. The study also discusses a method of formulating the difference equations pertaining to the two-dimensional Laplacian in a system containing rectilinear barriers. In addition, a technique is described for tracing the progress of a free surface, using marked particles, in a stream containing discontinuities SuCh as those occurring at the extremities of barriers. The above computational scheme was used in a digital computer program to simulate results of the Hele-Shaw analog. The computed values of sweep efficiency mere in good agreement with the experimental values. INTRODUCTION At the present time waterflooding remains the most widely used technique for the secondary recovery of oil. Among the numerous factors that affect the performance of a waterflood, the areal sweep efficiency constitutes one of the most critical parameters. It is a function only of the geometric distribution of the wells for a homogeneous medium and unit mobility ratio. Geological considerations, however, such as the presence of sealing faults or solution channels in the rock matrix, would generally exert an adverse effect on the areal sweep efficiency of a waterflood network. Surprisingly little information is presently available on the quantitative implications of such permeability interferences. It was the object of this investigation to consider the presence of rectilinear, impermeable and highly permeable interferences in an otherwise homogeneous five-spot pattern, and to analyze their effects on the resulting sweep efficiency and pattern conductivity. The experimental phase of the study was conducted using two different models: the Hele-Shaw flow cell and the field plotter. For experimental convenience, the geometric boundaries of the five-spot were assumed to constitute streamlines even in the case where the interference was asymmetrically placed with respect to a repeated five-spot pattern. A second object of the study was to develop a computational scheme for numerically calculating the pressure distribution in a five-spot flow pattern containing rectilinear, impermeable barriers, and to use this distribution to compute the sweep efficiency. EXPERIMENTAL PROCEDURE HELE-SHAW ANALOG The Hele-Shaw analog was used both as a visual model and for determining the sweep efficiency for systems involving impermeable barriers. The physical limitations of this model as a means of experimentally studying the motion of a fluid in porous medium have been adequately discussed elsewhere.2,7 Likewise, it can be easily verified that the analogy remains valid when a portion of the region between the plates contains an obstacle of the same thickness as the distance between the plates. The model used in this study was composed of

By J. N. Breston, R. V. Hughes

Conclusions drawn by previous research workers with reSPect to the relation between Pressure gradients and/or velocity and oil recovery obtained by laboratory water flood tests have been in disagreement probably due to variable procedures and unnatural conditions and materials. The Bradford Laboratory of the Pennsylvania Grade Crude Oil Association as part of its secondary recovery research program has conducted nineteen water floods on two long cores of widely differing characteristics in an attempt to clarify this relationship and make it an aid in predicting flooding pressures in the field. Unlike previous research procedures the present experiments were conducted with the aim of duplicating field conditions as closely as possible by using long unextracted consolidated cores, a live crude, and natural brines for both flooding and, connate water content. Also, the pressure gradients and flooding velocities were representative of field conditions where similar sands were being flooded. Eleven floods on one core and eight floods on the other core showed increased recoveries and lower residual oil saturation with increased flood pressure gradients and flood velocities. A marked decrease in recovery was obtained from both cores at very low flood velocities. This pressure versus recovery relationship is shown to hold up to the point of water breakthrough and also up to the 100 and 1 produced water to oil ratio point. INTRODUCTION The possibility of water flooding oil sands was suggested by Carl1 of the Pennsylvania Geological Survey in 1880. It is not known when the practice was tried intentionally for the first time, but its beneficial effects were noted in the annual production rate of the Bradford field as early as 1907. The practice was illegal in Pennsylvania until 1921. Early water floods in the Bradford field usually consisted of shooting or splitting the casing secretly to permit subsurface waters to enter the producing sands under hydrostatic head. As it was noted that the benefits of water flooding seemed to be proportional to the quantity of water dumped into the well many also began to utilize surface sources after the practice became legal. It was probably during the middle 20's before many producers realized that the pressurehead of the water upon the producing sand determined the rate and quantity of water that would enter the sand. Hence, rate and quantity of production appeared to be a direct function of input pressures. By 1927 a few producers had ventured the installation of pressure pumps in order to increase water-input rates and production through the combination of hydrostatic and hydraulic pressures. The adoption of pressure flooding and the "five-spot" drilling pattern in the Bradford field were essentially simultaneous. Water-input pressures in 1930 seldom exceeded 300 p.s.i. at the well head or 1100 p.s.i. at the sand face. Since that time, water-flood producers in the Bradford-Allegany fields have gone to higher and higher pressure until today 600 p.s.i. at the well heads is called a low pressure flood. Many high pressure floods now operate at 1300-1400 p.s.i. at the well head. The limiting and advisable pressure at the sand face has been pronounced as that pressure just under what is required to lift the overburden or to cause formation parting: According to this rule any water flood operation utilizing well-head pressures nearly equal in pounds per sq. in. to 1.1 times the average depth in feet to the top of the producing sand would be considered as a high-pressure flood. Despite the higher pressure trends in the Bradford-Allegany field operations and the results of early laboratory water flooding research, the desirability and benefits of high input pressures are still questioned by many operators, particularly in midwestern water flood operations. The present paper recounts a series of 19 laboratory water floods using two long, consolidated cores of widely differing characteristics saturated with live crude and flooded with oil field brines in an effort to simulate field conditions under • various pressure gradients and flooding velocities. For both cores, higher recoveries and lower residual oil saturations were obtained at higher pressure gradients and flood velocities. This relationship is shown to hold up to the water breakthrough point and also up to the 100 to 1 produced water to oil ratio point.

Jan 1, 1949

By R. V. Hughes, J. N. Breston

Conclusions drawn by previous research workers with reSPect to the relation between Pressure gradients and/or velocity and oil recovery obtained by laboratory water flood tests have been in disagreement probably due to variable procedures and unnatural conditions and materials. The Bradford Laboratory of the Pennsylvania Grade Crude Oil Association as part of its secondary recovery research program has conducted nineteen water floods on two long cores of widely differing characteristics in an attempt to clarify this relationship and make it an aid in predicting flooding pressures in the field. Unlike previous research procedures the present experiments were conducted with the aim of duplicating field conditions as closely as possible by using long unextracted consolidated cores, a live crude, and natural brines for both flooding and, connate water content. Also, the pressure gradients and flooding velocities were representative of field conditions where similar sands were being flooded. Eleven floods on one core and eight floods on the other core showed increased recoveries and lower residual oil saturation with increased flood pressure gradients and flood velocities. A marked decrease in recovery was obtained from both cores at very low flood velocities. This pressure versus recovery relationship is shown to hold up to the point of water breakthrough and also up to the 100 and 1 produced water to oil ratio point. INTRODUCTION The possibility of water flooding oil sands was suggested by Carl1 of the Pennsylvania Geological Survey in 1880. It is not known when the practice was tried intentionally for the first time, but its beneficial effects were noted in the annual production rate of the Bradford field as early as 1907. The practice was illegal in Pennsylvania until 1921. Early water floods in the Bradford field usually consisted of shooting or splitting the casing secretly to permit subsurface waters to enter the producing sands under hydrostatic head. As it was noted that the benefits of water flooding seemed to be proportional to the quantity of water dumped into the well many also began to utilize surface sources after the practice became legal. It was probably during the middle 20's before many producers realized that the pressurehead of the water upon the producing sand determined the rate and quantity of water that would enter the sand. Hence, rate and quantity of production appeared to be a direct function of input pressures. By 1927 a few producers had ventured the installation of pressure pumps in order to increase water-input rates and production through the combination of hydrostatic and hydraulic pressures. The adoption of pressure flooding and the "five-spot" drilling pattern in the Bradford field were essentially simultaneous. Water-input pressures in 1930 seldom exceeded 300 p.s.i. at the well head or 1100 p.s.i. at the sand face. Since that time, water-flood producers in the Bradford-Allegany fields have gone to higher and higher pressure until today 600 p.s.i. at the well heads is called a low pressure flood. Many high pressure floods now operate at 1300-1400 p.s.i. at the well head. The limiting and advisable pressure at the sand face has been pronounced as that pressure just under what is required to lift the overburden or to cause formation parting: According to this rule any water flood operation utilizing well-head pressures nearly equal in pounds per sq. in. to 1.1 times the average depth in feet to the top of the producing sand would be considered as a high-pressure flood. Despite the higher pressure trends in the Bradford-Allegany field operations and the results of early laboratory water flooding research, the desirability and benefits of high input pressures are still questioned by many operators, particularly in midwestern water flood operations. The present paper recounts a series of 19 laboratory water floods using two long, consolidated cores of widely differing characteristics saturated with live crude and flooded with oil field brines in an effort to simulate field conditions under • various pressure gradients and flooding velocities. For both cores, higher recoveries and lower residual oil saturations were obtained at higher pressure gradients and flood velocities. This relationship is shown to hold up to the water breakthrough point and also up to the 100 to 1 produced water to oil ratio point.

Jan 1, 1949

By A. C. England, William F. Haney, Paul B. Crawford, Bobby L. Landrum

The creation of unsymmetrical vertical or horizontal fractures in a producing well will result in a unique distortion of the normal shut-in pressures existing at the four surrounding wells. An electrical model study has been made to show how the shut-in pressures at the wells surrounding the fractured producing well may he used to estimate the size and orientation of vertical or horizontal fractures. The pressure differences at the four surrounding wells are recorded during production from the fractured well. The resulting pressure dato are matched from a table of several hundred possible fracture combinations. The table refers to the figure showing the effective size and orientation of the commercially created fracture. INTRODUCTION The creation of commercial fractures about oil and gas wells has led to considerable speculation as to the size and shape of the fractures. It is important that the size and shape be known in certain instances, especially when pressure maintenance programs are to be initiated. It was shown that if knowledge of a vertical fracture could be found that such information may be utilized in planning a waterflooding or gas cycling program and that in some cases the fractures might be utilized very effectively.' Other work showed that it may be desirable to create "pancake" fractures in certain reservoirs but that they should exist in a desirable range.' If the fractures are too large, they should be plugged prior to initiating a pressure maintenance program. If an irregular shaped "pancake" fracture is created, it is desirable to know the orientation of this fracture because similar but unreported studies have shown that "pancake" fractures may or may not substantially affect the sweep efficiencies and performance of the flooding or cycling program. Consequently, it may be desirable in some instances to estimate the effective size and orientation of the fracture. ANALOG An electrical model was used in making this study. The electrical model consisted of a circular reservoir approximately 8 in. in diameter, 1 in. thick and was filled with a dilute salt solution. A copper band formed the periphery of the reservoir. A vertical fracture was simulated by a copper strip soldered to the wire well located in the center of the reservoir. Vertical copper strips of various lengths and orientation were soldered to the wire well to represent vertical fractures in the reservoir. In studies reported here on vertical fractures it was arbitrarily assumed that two fractures originated at the well and extended out into the reservoir. An electrical potential was imposed between the wire well and the periphery of the tank. Plots of the iso-potentials are shown in Figs. 1 through 9. This model provides the steady-state isopotentials existing in a uniform, homogeneous circular reservoir having a vertically fractured well at its center. The permeability of the fracture is much greater than the permeability of the matrix proper. The effects of gravity are neglected. In preparing Tables 1 and 2, four symmetrically located unfractured wells were assumed to exist around

Jan 1, 1956

By Alton B. Cook, G. B. Spencer, F. P. Bobrowski

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By F. P. Bobrowski, G. B. Spencer, Alton B. Cook

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By John C. Calhoun, Raymond G. Loper

Theoretical calculations for the decline of pressure and the variation of instantaneous producing gas-oil ratio with increased cumulative production have been made for reservoir systems under various producing drawdown values and for various spacings. These variables have been included by a modification of previously Published techniques. This consists essentially of utilizing a radial system flowing gas and oil in a succession of steady state calculations. The results of the calculations do not show any variation in the predicted behavior of the internal gas drive reservoir when either the drawdown or the spacing are changed. The prediction of performance on an internal gas drive reservoir has been presented by various authors. ',"' These predictions are in the form of a pressure decline and an instantaneous producing gas-oil ratio as a function of cumulative oil withdrawal. In making these predictions a basic assumption has been that the saturation of the reservoir is uniform throughout at any time, or that the pressure differential imposed on the reservoir system is zero. Because of these assumptions the validity of such predictions has been questioned. Their application to actual systems, where saturations are not uniform and where producing differentials exist, has been a hazardous procedure. The authors cited have used different methods of approach to the derivation of equations for making their predictions and for arriving at a solution of both equations, but basically the same relationships are utilized in each case. In all instances the instantaneous producing gas-oil ratio is calculated from the equation: In utilizing this equation an average reservoir pressure and an average reservoir saturation must be known, the former to define ß, v, s, and uo/ug and the latter to define kg/ko. Predictions can be made, therefore, only for the conditions of a uniform reservoir saturation and a uniform reservoir pressure. This restriction of requiring a uniform oil saturation can be resolved by applying two-phase steady state flow for the radial system to the reservoir's history. Radial steady state flow is defined when the ratio of phases flowing,

Jan 1, 1949