Search Documents

By D. K. Atwood

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

Jan 1, 1965

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

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

Jan 1, 1969

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

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

Jan 1, 1953

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

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

By A. G. Weber, K. H. Coats, M. H. Terhune, R. L. Nielsen

Two computer-oriented techniques for simulating the three-dimensional flow behavior of two fluid phases in petroleum reservoirs were developed. Under the first technique the flow equations are solved to model three-dimensional flow in a reservoir. The second technique was developed for modeling flow in three-dimensional media that have sufficiently high permeability in the vertical direction so that vertical flow is not seriously restricted. Since this latter technique is a modified two-dimensional meal analysis, suitably structured three-dimensional reservoirs can be simulated at considerably lower computational expenses than is required using the three-dimensional analysis. A quantitative criterion is provided for determining when vertical communication is good enough to permit use of the modified two-dimensional areal analysis. The equations solved by both techniques treat both fluids as compressible, and, for gas-oil applications, provide for the evolution of dissolved gas. Accounted for are the effects of relative permeability, capillary pressure and gravity in addition to reservoir geomety and rock heterogeneity. Calculations are compared with laboratory waterflood data to indicate the validity of the analyses. Other results were calculated with both techniques which show the equivalence of the two solutions for reservoirs satisfying the vertical communication criterion. INTRODUCTION Obtaining the maximum profits from oil and gas reservoirs during all stages of depletion is the fundamental charge to the reservoir engineering profession. In recent years much quantitative assistance in evaluating field development programs has been provided by computerized techniques for predicting reservoir flow behavior. Because of the spatially distributed and dynamic nature of producing operations, automatic optimization procedures, such as those now in use for planning refining operations, are not now available for planning reservoir development. However, present mathematical simulation techniques do furnish powerful means for making case studies to help in planning primary recovery operations and in selecting and timing supplemental recovery operations. A number of methods have been reported which simulate the flow of one, two or three fluid phases within porous media of one or two effective spatial dimensions>-4 However, applying computer analyses to actual reservoirs have been limited mostly to two-dimensional areal or cross-sectional flow studies for two immiscible reservoir fluids. To obtain a three-dimensional picture of reservoir performance using such two-dimensional techniques, it has been necessary to interpret the calculations by combining somehow the results from essentially independent areal and cross-sectional studies. To the authors' knowledge, the only other three-dimensional computational procedure, in addition to those presented here, was developed by Peaceman and Rachford8 to simulate the behavior of a laboratory waterflood. Two computational techniques which may be used to simulate three-dimensional flow of two fluid phases are described in this paper. The first method, called the "three-dimensional analysis", employs a fully three-dimensional mathematical model that treats simultaneously both the areal and cross-sectional aspects of reservoir flow. The second method, called the vertical equilibrium (VE) analysis", is applicable

By C. Perez-Rosales

A statistical method for determining simultaneously some of the basic geometrical characteristics of porous media such as porsity, specific surface* mean pore width, mean grain thickness and absolute permeability is presented. The proposed method is characterized by its simplicity and the minimum amount of laboratory equipment that it requires. The experimental procedure used to evaluate the various geometrical characteristics is based 'pan the 'surface analysis Of a a In view of this fact the app1icability Of the method is limited to homogeneous and isotropic materials. INTRODUCTION The behavior of a fluid flowing through a porous medium depends, among other things, upon the internal geometry of the medium; hence the importance of developing efficient methods to determine the geometrical properties of porous materials. This paper describes a statistical method to determine simultaneously some of the basic geometrical properties of porous materials such as porosity, specific surface, mean pore width, mean grain thickness and absolute permeability. The mathematical formulation has been developed that the data required to calculate the various geometrical characteristics of a sample can be easily measured by analyzing a section of the sample with an evenly spaced grid. THEORY A simple way of obtaining information about the internal geometry of a porous material is to throw a point at random over a cross-section of a sample. l-Through this procedure the porosity can be detemined.l,2 By considering that the material is homogeneous and isotropic, and that the point is dropped many times, the porosity 4 is given by where N is the total number of times the point is thrown and n is he number of times the point falls within pore areas. Another simple manner of analyzing the structural characteristics of a porous material is to superimpose an arbitrarily long line on a cross-section so that its length is evenly distributed over the surface area. Through this procedure the porosity can also be determined. If L is the total length of the line, and 1 is the sum of the lengths of the line segments within void spaces, the porosity is given by 4 = —.1.............. (2) The advantage of this type of analysis is that geometrical characteristics other than porosity can be obtained. Thus, if c represents the number of intersections between the line and the perimeter of pores, it can be shown3 that the specific surface, defined as he surface area of pores per unit bulk volume, is given by To give a proper description of the internal geometry of a porous medium, it is necessary to define quantities that characterize the notions of pore size and grain size. Unfortunately, because of the complexity of porous materials, it is difficult to give exact geometrical definitions of what is meant by the ConCePtS of "poresl' and "grains", especially when dealing with consolidated materials. Never theless, one often talks about the "size of pores" and the "size of grains" without defining accurately the terms. A possible solution to this problem follows. Assume that an arbitrarily long line is placed on a section of a sample so that the line is evenly distributed over the surface area (see Fig. 1). 1n a system like this, each line segment within a void space will be, by definition, the pore width in a given location and direction. Consequently, if the

Jan 1, 1970

By C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By T. N. Dixon, J. E. Briggs

A study was made of numerical techniques for solving the large sets of simultaneous equations that arise in the mathematical mode ling of oil reservoir behavior. It was found that noniterative techniques, such as the Alternating Direction Implicit (ADI) method, as well as some other finite difference approximations, produce oscillatory or unsmooth results for large time steps. Estimates of time step sizes sufficient to avoid such behavior are given. A comparison was made of the Point Successive Over-Relaxation (PSOR), Two-Line Cyclic Chebyshev Semi-Iterative SOR (2LCC), and iterative ADI methods, with respect to speed of solution of a test problem. It was found that, when applicable, iterative ADI is fastest for problems involving many points, while 2LCC is preferable for smaller problems. INTRODUCTION With the advent of high speed, large memory, digital computers, there has been an increasing emphasis on the development of improved methods for simulating and predicting reservoir performance. Two-dimensional, three-phase reservoir models with various combinations of PVT effects, as well as gravity and capillary forces, are common throughout the industry. Such models are also available through consulting firms, to anyone desiring to use them. Three-dimensional models will probably be practical in only a few years. We conducted a study of some of the numerical methods used for solving the large sets of simultaneous equations that arise in such models. A typical set of equations for a reservoir model is shown below:

Jan 1, 1969

By A. S. Odeh

This paper analyzes the effect of limited entry to flow at the wellbore on the steady-state productivity of a well. Wells that have been opened to flow along a fraction of their productive interval are termed wells with limited entry. Previous work treated the cases of a partially penetrating well, a well producing from the central portion of the productive interval and a well in which several intervals equally spaced were open to flow. In this paper the open interval can be located anywhere within the productive interval. Thus, in a sense, it generalizes previous work. The finite cosine transform was used to arrive at a solution for steady-state flow of a slightly compressible fluid. The solution was programmed for a CDC 1604 computer. Numerical vaIues for rd = 660 ft, r, = 1/4 ft, and range of sand thickness of 20 to 200 ft are presented in graphical form. The effect of rd and r, values on the result is shown in a table. The correct calculation of skin and damage ratio in the presence of limited entry to flow is explained and illustrated by examples. Moreover, the paper shows how to calculate the net decrease in productivity due to the combined effect of limited entry and perforations. INTRODUCTION In some wells only a fraction of the productive interval is open to flow. Location of this fraction is usually dictated by formation characteristics and reservoir behavior. For instance, if a gas cap exists, the open interval is located away from the gas-oil contact to prevent any possible gas coning. Wells that intentionally have been opened to flow along a fraction of their productive formation are tened wells with limited entry. Obviously, unintentional completions of this type also exist. Limited entry to flow decreases well productivity. Magnitude of the loss depends on the fraction of the formation open to flow, on the thickness of the sand, on the location of the open interval and on the ratio of rd /r, where r , is well radius and rd is the drainage radius of the well. The use of pressure buildup data on producing wells to calculate the condition of the formation around the wellbore is an accepted practice. van Everdingenl and Hurst2 introduced the concept of she skin factor s considered to be due to a thin layer of different permeability immediately around the wellbore. These authors dealt with the case of a well of complete radial geometry, i.e., a well with open-hole completion that completely penetrates the formation. The presence of a low-permeability skin results in a loss of productivity, as does limited entry. Therefore, if pressure buildup data obtained on a well with limited entry are used to establish the presence or absence of skin (i.e., formation damage), and a correction is not made for this loss of productivity, the calculations would result in an erroneous skin value. They might indicate the presence of formation damage when in reality there is none, or they might indicate a value larger than the true value. This could lead to an incorrect basis for planning remedial measures. Muskat3 studied the problem of partially penetrating wells for the case of incompressible flow. He presented equations and figures which allow estimation of loss in productivity. Brons and Marting,4 using equations based on Nisle's work,5 studied the loss of productivity for three cases. The first was for a partially penetrating well; the second was for a well producing from only the central portion of a productive interval; and the third was for a well in which several intervals equally spaced were open to flow. Their work was for steady-state depletion-type reservoirs wherein the well radius of drainage is established and the fluid is considered to be slightly compressible. Considered in this paper is the problem of wells with limited entry in which the open intervals are located anywhere within the productive sand. The finite cosine transform is used to arrive at a

Jan 1, 1969

By M. P. Tixier, F. Segesman

Explanations bused on theory and laboratory tests tire offered for certain unusual features of SP logs which may befog their interpretation, namely the "sawtooth" SP and the reduction or increase of SP with time. The sawtooth SP is frequently seen in SP recordings of thick, very permeable watcr. sands of the Louisiana Gulf Coast. Using SP data from a resistor network analog, this anomaly is explained on the basis of a thin shale streak within the sand. under which the invading filtrate accumulates by gravity segregcitiorr. The electromotive force then appearing across the shale streak gives rise to the anomalous SP. Reduction of the SP on successive logging runs Ira often been observed in very permeable sands when as indicated by the resistivity curves, the invasion depth has been reduced to practically zero. Laboratory results verify that this phenomenon may be due, to the presence of air electromotive force, across the mod rake when the fluids on either side differ in activity. Also, there is, in some cases; an increase of SP deflection on successive logging runs, because the emf initially present across the mud cake disappears as the somewhat saline water just behind the mud cake is progressively flushed out by the invading filtrate.. Field example and laboratory data are presented. I N T R ODUCT1ON Some unusual features are frequently observed on certain SP logs which make interpretation difficult. The purpose of this paper is to offer explanations of several cases of interest. All of these arc related to the conditions of invasion. The cases to be discussed are as follows. 1. The sawtooth SP anomaly which is frequently seen within the thick, highly permeable saltwater sands of the Louisiuana Gulf Coast 1. Reduction of the SP deflection from the usual sand line in very permeable sections when the depth of invasion is practically nil. As will be explained, such a reduction is attributable to a membrane electromotive force across the mud cake. 3. An increase of the SP with time just after the iormation is drilled as a result of progressive flushing of the invaded zone. FUNDAMENTALS OF SP PHENOMENA An explanation of SP phenomena must begin with the following basic ideas. 1. The SP log derives its character from ohmic potential differences in the hole produced by currents flowing through the mud.' These currents are driven by electromotive forces which are primarily electrochemical in origin. Other possible types of SP-producing electromotive forces may be disregarded for the purposes of this paper. 2. The electrochemical emf's associated with permeable formations are now understood to consist of liquid-junction emf's and membrane emf's.',' As schematically shown by thc battery symbols in Fig. 1, the liquid-junction emf, E,,,, exists at the contact between waters of different chemical activities, a,,,, (mud filtrate) and tr,,. (formation water); the shale membrane emf, lS.sn. exists across the adjacent shales whcn they separate waters of different activities. The emf's just cited are additivc in producing SP currents (arrowed circles of Fig. I). Thus, neglecting for the moment any possible membrane emf across thc mud cake, the total electrochemical emf, E,, will be given by their sum, where K,., -! K,q,, = 12 + 59, or 71 at 75 °F. In the log interpretations the SP of a thick sand (i.e.. the static SP) is generally related to the activities of the filtrate and formation water through Eq. 1. by setting the static SP equal to - E,.

By Robert G. Nisle

The classic theory of pressure build-up in shut-in oil wells as developed by Horner and van Everdingen is based on two-dimensional radial symmetry in the well-reservoir system. Such symmetry does not exist in the case of a well which parfiially penetrates the producing formation. As a result of this lack of symmetry the use of the classic theory in such cases become questionable. In this paper the mathematical theory has been extended to include the case of partially penetrating wells. Numerical solutions illustrating the up-dication of the equations are presented. The effect of partial penetration on pressure build-up is shown by a comparison of synthetic pressure build-up curves derived from the numerical solution of the equation for partially penetrating wells for various degrees of penetration. It is shown that partial penetration is detectable from the characteristic shape of the pressure build-up curve and that formation productivity may be calculated from the pressure buildup data in a manner identical to that described by the classic theory. INTRO DUCTION The use of pressure build-up data on shut-in oil wells is a well-established technique for the measurement of reservoir productivity. The work of Homer' and van Everdingen' is well known. Since the publication of their pioneering work, others"' have elaborated and extended the technique. These authors have dealt exclusively with the case of wells that completely penetrate he producing formation. Two-dimensional radial symmetry is achieved, thereby, and the problem is greatly simplified. This ideal situation is seldom encountered in practice. But, in spite of this, the technique is often employed anyway. The question, therefore, arises: how much confidence can be placed in the results of such calculations when the data are taken on wells that only partially penetrate the producing formation? It is the purpose of this work to answer this question. Only the single-phase fluid case will be considered. van Everdigen' has presented a rigorous solution of the two-dimensional case, but because of mathematical complexity, this solution has not proved to be particularly useful. Horner' utilized Kelvin's point source solution in two dimensions and showed that it is a sufficiently good approximation for oil reservoir studies. The geometry of the well and the producing formation is illustrated in Fig. l. Geometrically the formation may be described as an infinite slab. The infinite slab is defined as the solid (the porous medium in this case) bounded by two parallel planes. Thus, the slab extends to very great distances in both the positive and negative directions of two of the coordinate axes, but is finite in

By J. J. Arps, T. G. Roberts

Since the introduction of the relative permeability concept in the middle thirties1,2 various investigators have shown3,10,11,12,15 how the basic equations for the flow of oil and gas through porous media can he utilized to compute the recovery from depletion type or solution gas drive type reservoirs under certain conditions, when the necessary physical data pertaining to the reservoir rock and its reservoir fluids are known. In the study by Muskat and Taylor" in 1945, the effects of viscosity, gas solubility, shrinkage, and the permeability-saturation characteristics of the producing formation on the production histories and the recovery of gas drive reservoirs were analyzed. Each parameter was varied a limited number of times, while keeping the others constant. The work was not carried to the point where the data could be used directly to estimate the primary recovery in those cases where certain rock and fluid characteristics were known or could be assumed. Since the publication of this classic study, the work by Beal in 194616 and Standing in 1952'" has established general correlations of oil viscosity and shrinkage with gas solubility and oil gravity. In addition. a considerable number of relative permeability relationships for different types of reservoir rocks has been published2,5,13,17,19,20. Also, the burden of the numerical work required in the solutions of the differential equation has. during the last few years, been greatly reduced with the availability of modern electronic computing facilities. The purpose of this paper is therefore to present the results of a large number of such recovery computations covering the normal range of variation of the main parameters such as the type of rock, the API gravity of the oil, and the amount of gas in solution. As is to be expected, the ultimate recovery is found to increase with the oil gravity. except for the higher solution gas-oil ratios. Within the normal range of solubilities, the ultimate recovery generally appears to decrease with the amount of gas in solution for the higher oil gravities, and to show a slight increase for the lower oil gravities. Intermediate oil gravities seem to fall between the two trends. The type of reservoir rock, as identified by its relative permeability relationship, appears to have a very pronounced effect on the ultimate recovery. In general, sands and sandstones show higher recoveries than limestones and dolomites, although certain intergranular limestones seem to show a higher theoretical recovery than unconsolidated sand. The recovery from sands and sandstones generally decreases with increasing cementation and consolidation, while the recoveries from limestones and dolomites are highest for the intergranular type and lowest for the fracture type porosity.

Jan 1, 1956

By D. Havlena, A. S. Odeh

The use of the straight-line method of solving the material balance equation is illustrated by means of six field cases. Also, the application of statistical criteria to arrive at the most probable answer is shown. The theory underlying the straight-line method of solution and the applicability of the statistical criteria was presented in a previous paper.' The field cases include saturated and undersaturated oil reservoirs with and without water drive. The aquifers discussed are: limited radial, infinite radial, very small aquifer and infinite linear. The field cases also include a gas reservoir producing under water drive. INTRODUCTION In a previous paper,' the authors presented the theory underlying the solution of the material balance as an equation of a straight line. The appropriate equations for various material balance cases as well as the methods of analysis and interpretation with comments and discussion were also included. To illustrate the various theoretical cases treated previously, five field cases are analyzed in this paper by employing the straight-line method of solving the material balance equation (MBE) and one example previously published is referred to. The use of statistical criteria to arrive at the most probable answer is also shown. All the field examples presented, except Case 2, are excerpts from complete reservoir studies. To illustrate the method, only sections specifically dealing with the material balance principles are included. Additional geologic information and basic data are reported to better acquire an understanding of the cases and thus to better follow the reasoning that suggested the successful application of the straight-line method of solving the MBE. The six cases are: (1) saturated reservoir, small gas cap, limited aquifer; (2) saturated reservoir, very small gas cap, infinite aquifer; (3) undersaturated-saturated reservoir, very small aquifer; (4) highly undersaturated reservoir, no water drive; (5) high undersaturated one-well reservoir, limited aquifer; and (6) gas reservoir, infinite linear aquifer. THE D4 SAND, GUICO FIELD, VENEZUELA The D, sand, which was discovered in 1943, is presently in a depleted state. Since its discovery it has produced under water drive, gas-cap-gas expansion, and solution gas drive. In Nov., 1947, water injection was initiated to arrest further pressure decline. When discovered, the D, sand was a saturated reservoir with a gas cap /oil zone volume ratio m estimated volu-metrically at 0.0731, an average permeability of 500 md, a porosity value of 25 per cent, and an oil viscosity at reservoir conditions of 0.3 cp. The volumetrically determined stock-tank oil initially in place was 23.1 million bbl. The volumetrically weighted physical data and production data available until Nov., 1953 are reported in Table 1. In Ref. 1, the effects on the straight-line plot of various values of ro/r for a constant ?to, or of various dimen-sionless times for a constant r8/r10 were theorized and were illustrated in Fig. 3A of that reference. In this field case, the previously theoretically predicted effects are established. Thus, the MBE calculations using Eq. 3c of Ref. 1 were performed for various r6/r10 and dimension-less time values. Eq. 3c of Ref. 1 is:

Jan 1, 1965

By J. E. Berry

AII evaluation is made of the acoustic velocity log for measurement of formation porosity. Plots of field-observer1 velocities vs core-measured porosities of sandstones and limestotnes with inter intergranular porosity show that the velocity log gives a useful measure of porotity, in agreement with published data. A new method of plotting electrical and velocity log data makes it easier to recognize hydrocarbon-water coritnct~ from log darn and to make serni-quantitative estinlates of hydrocarbon saturation. Two examples of this interpretation technique are given. This paper also discusses the various factors that affect acoustic velocity and shows how corrections can be made for some of these factors. INTRODUCTION Although the continuous velocity log was developed primarily as an aid to seismic interpretation, it is finding widespread use for measurement of porosity. The purpose of this paper is to evaluate this use of the log. Results are presented for the velocity-porosity relationship in many sandstones and one limestone. Practical applications of porosity values derived from the velocity log to improve interpretation of electrical surveys are presented. Various parameters which affect velocity in porous media are considered. Porosity, composition, cementa-tion, pressure difference (overburden pressure minus fluid pressure), fluid saturation, and wellability are discussed. Some of the ideas are opinions based on general considerations of well logging problems while others are supported by published data and our own observations. Major conclusions from these studies arc given in the body of the paper and detailed discussions of some of the various factors affecting porosity are given in the Appendix. POROSITY Sandstones The primary factor affecting the velocity of sound in porous media is porosity. Wyllie, et al,1 have reported a time average equation for the relationship between velocity and porosity.* This has found considerable acceptance in the industry. Other relationships between velocity and porosity have been proposed, none of which are entirely satisfactory. It is doubtful that any one relationship will completely describe all situations. We have tried an empirical relationship based on correlating log velocities with core-measured porosities. Measured porosities were averaged over a depth interval equal to the spacing used in recording the velocity logs. The depth range for the information used was from 2,000 to 12,000 ft. Velocity measurements were normalized for the effects of depth and overburden pressure. A detailed discussion of the basis for this normalization is given in the Appendix. The change is usually insignificant and when working with individual formations it can be neglected. All velocity measurements were normalized to a depth of 10,000 ft. The velocity-porosity data for sandstones in Fig. 1 came from 20 wells scattered throughout the Mid-Continent and Gulf Coast areas. The scatter of the data points is attributed to three causes: 1. There is no theoretical basis for expecting an exact velocity vs porosity relationship and it is not possible to correct precisely for the effects of such factors as impurities, grain shape, manner of cementation, etc. Some of these factors are discussed in the Appendix. 2. Porosities measured by core analysis are not necessarily representative of what the velocity log "sees". 3. Inaccuracies appear in the velocity log reading. This is particularly true of velocities from one-receiver logging tools.*" A linear-time curve (time-average equation) is plotted in Fig. 1 along with the best straight line through the points. Velocities used for the linear-time curves were 19,500 ft/sec for the matrix velocity (zero porosity) and 5,500 ft/sec for fluid velocity. For clean sandstones with velocities greater than 10,000 ft/sec, there is little difference between the velocity-porosity and the time-porosity curves. Both curves are well within the scatter of data points. The data points in Fig. 1 were restricted primarily to clean sandstones, using the SP and/or gamma-ray logs as clues to the lithology. Where a lithologic de-

By W. F. Kieschnick, and W. J. McGuire&apos, Eugene Harrison

This paper concerns the induction and extension of fractures into rock formations as involved in drilling, completing, and production stimulating operations on wells. Conclusions concerning formation breakdown are derived from (I) a review and extension of published analyses relating to mechanical theories of rock stress and the state of stress in the earth's crust and (2) a correlation of field data from fracturing operations. Conclusions concerning the mechanics of fracture extension, which indicate the relationship between fracture dimensions and rock properties, depth, and volume of injected fluid, are tentative and largely establish limits of relationships. These conclusions are derived from stress calculations, limited field data, and laboratory experimental studies. The experimental work involves the study of the stresses at the fracture boundaries and the geometry of pressurized fractures by means of photo-elastic modeling methods. Results of this investigation indicate that a large majority of pressure induced wellbore fractures are vertical, particularly in deeper wells; and variations in the pressures necessary to create and extend fractures can be explained largely on a basis of established rock properties. It is also shown that variations due to tectonic forces should usually be expected to be slight. Other results indicate that during that extension of fractures rather large fracture volumes are temporarily created by the parting of the formation. The purpose of this paper is to present the results of calculations and laboratory experiments concerning the mechanics of fracture induction and extension with a view to broadening existing knowledge relating to these phenomena. It is believed that continued progress in developing knowledge of this type is important to the further development of techniques for drilling and completing wells. When a wellbore is subjected to sufficiently high pressure, fracturing of the formation surrounding the well occurs. This phenomenon is a common cause of lost circulation during drilling operations, is a deliberately planned factor in most squeeze cementing operations, and is the controlling factor in many productivity stimulation operations. The exploitation of formation fracturing in recent years has been particularly impressive in sand-viscous fluid well stimulation treatments which were first introduced to the industry by the Stanolind research group1,2,3 several years ago. The borehole fracturing phenomenon has been analyzed mathematically by application of the principles defining the elastic and inelastic behavior of thick-walled cylinders of homogeneous, isotropic, and impermeable material.1,5,6 This basic argument is re-presented in this paper along with new considerations which strengthen it. In addition to these arguments relating to fracture induction some calculations and experimental results relating to the mechanics of fracture extension are presented. Since the nature of the pressure induced fracture at distances from the well is of great importance in considerations of well stimulation procedures there appears to be a need for an integrated approach which deals with the mechanics of fracture extension as well as with wellbore rupture.

Jan 1, 1955

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-