Search Documents

By W. T. Weller

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

Jan 1, 1967

By D. K. Atwood

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

Jan 1, 1965

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

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

Jan 1, 1969

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

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

Jan 1, 1953

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

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

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

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

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 M. W. Clegg

The application of the Laplace transformation to problems in the flow of compressible fluids in porous media has provided a large number of exact solutions. For plane radial flow, however, these solutions are either complex integrals or infinite series and are of little value to the field engineer. In the case of production at constant well pressure, the available approximate solutions me valid for large times only. In this paper it is shown that an approximate inversion formula for the Laplace transform, developed for the solution of viscoelastic problems, is applicable to radial flow problems and provides simple analytical solutions to constant terminal pressure problems. The method may be used to obtain approximate solutions to many problems, including media with radial permeability discontinuities, multi-layer formations and pressure buildup in wells after shut-in. The results are compared with the few available computer solutions as well as the large time solutions, and it is shown that this approximate method greatly extends the time interval over which a simple analytical solution is acceptable. INTRODUCTION The study of transient problems in the flow of fluids through porous media has benefited greatly from the application of transform methods. The use of the Laplace transformation for solving parabolic equations has been widely discussed in the field of heat conduction and diffusion as well as in the petroleum literature. Removal of the time variable with the Laplace transformation generally reduces the problem to a boundary value problem which may be solved by standard techniques. A much more formidable problem then faces the engineer, however, for frequently the transform does not possess a simple inverse. The result is that the general inversion integral must be used and this leads to either an infinite integral or an infinite series, both of which are difficult to handle from a computational standpoint. Asymptotic approximations for the inverse have been known for some time and these yield approximate inverse functions that are valid for very large or very small times — but frequently the times of interest lie somewhere between these two extremes. Therefore, some acceptable approximation valid over a larger interval of time is desirable. During the past few years a number of methods for achieving this have been developed and some of these are discussed briefly in this paper. The relative merits of the various methods are not evaluated here, but some general conclusions reached by other authors are given. One of these methods has been applied to problems associated with the radial flow of compressible liquids to producing wells. In the case of production at constant well pressure, the method leads to simple analytical solutions for a number of standard problems; e.g., homogeneous formation, permeability discontinuities, pressure buildup. These solutions greatly extend the range of validity of the asymptotic ones (valid for large times only) and should be of value in studying the behavior of wells producing under constant pressure conditions. METHOD FOR THE APPROXIMATE INVERSION OF LAPLACE TRANSFORMS Much of the work that has been done on the approximate inversion of Laplace transforms has originated in viscoelastic problems. Apart from the asymptotic approximations,1 one of the first approximate methods was due to Alfrey.2 If (s) is the Laplace transform of f(t) defined by

By R. C. Earlougher

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 B. S. Gottfried

This paper is concerned with possible transport mechanisms which occur during segregated burning (i.e., burning in an oil reservoir in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas). ,Two analytical models are presented. In the first, burning is assumed to occur at the gas-oil interface, and the dominant mechanism is taken to be the rate at which oxygen diffuses to the interface. The second model assumes that combustion occurs within the clean zone, with the dominant mechanism being the mechanism of fuel deposition from the oil-bearing formation to the clean zone. ,Comparison of the results of the analyses with limited experimental data indicates that the latter mechanism is much more plausible. Both models are idealized representations of the actual process and are somewhat conjectural, however, owing to the lack of more definitive laboratory or field data. INTRODUCTION In recent years there has been considerable theoretical and experimental research directed toward understanding the recovery of crude oil by underground combustion. Although two different recovery schemes have been proposed, known as forward and reverse combustion, it is only forward combustion that shows promise of becoming a practical recovery technique.17 Most of the analyses of the forward combustion process have considered it to be pistonlike; i.e., oil is displaced by a vertical combustion zone which, with its ensuing gases, pushes die oil forward as the flame front is propagated through the reservoir.l-8,ll,l2,14-16 However, information gained in field pilot studies indicates that, at least to a certain extent, forward combustion proceeds in a segregated manner in oil reservoirs.10 An ideal segregated system is one in which the oil-bearing formation is overlain by a "clean" porous zone containing only gas. The location of the combustion zone in such a system is not immediately evident, although combustion is believed to occur either at the gas-oil interface or within the clean zone. In either case, the combustion zone may be horizontal and it moves in parallel rather than in series with the produced oil. The purpose of this paper is to explore the likely mechanisms by which segregated burning can occur, and to determine the sensitivity of the resulting models to changes in various physical parameters. Two governing mechanisms are considered: the case of a horizontal combustion zone at the gas-oil interface with the rate of combustion governed by the rate at which oxygen diffuses to the interface, and the case of combustion within the clean zone with the combustion process controlled by the transport of fuel from the oil-bearing formation to the clean zone. Both models are somewhat hypothetical and involve considerable simplification of the actual process. However, the results allow some preliminary understanding of segregated burning which is useful in the absence of more definitive laboratory or field pilot data.

Jan 1, 1967

By S. T. Yuster

An analysis has been made of the factors responsible for tilted fluid contacts in petroleum reservoirs. The factors are static and dynamic with the former being controlled by those variables responsible for the capillary rise of fluids, and the latter include rate of formation tilt and hydraulic circulation. The explanation that rate of formation tilt may be responsible for tilted contacts appears untenable. The role of some of these factors in the migration, fractionation. accumulation and possibly the prospecting for oil is indicated. INTRODUCTION While it is generally considered that the segregation of gas, oil, and water in a petroleum reservoir results in approximately horizontal planes separating these phases. there are factor; which modify this simplified picture to a very great extent. Some of these factors are static. and the departure of the fluid interfaces from horizontal may represent a static equlibriun. Some of these factors may be dynamic and the interface positions may be non-equilibrium and represent a steady .state condition. Of course combinations of these two situations may also exist. The reasons most commonly given for this behavior are hydraulic circulation of fluids, change in character of the formation, and a rate of tilt of the formations more rapid than the rate at which the fluids may re-establish an equilibrium. Russell's Principles of Petroleum Geology' contains an extensive discussion of tilted fluid interfaces and their importance in the migration and accumulation of petroleum. Russell cites examples of tilted water tables at Lance Creek, Wyo.; Rock River, Wyo.: Cut Bank. Mont.; Cushing, Okla.: and Graham; Okla. Also given was the Wasson Field of West Texas which has saucer shaped oil-water and gas-oil interfaces with the saucer being deeper in the case of the oil-water contact. The Bradford Field in Pennsylvania has a tilted water table. San Ardo. Coalinga None and some other California field.; have tilled contacts. A recent paper by King Hubbert discuses the effect of hydraulic circulation of the water underlying the oil in a given stratun and it- implications in the accumulation of and drilling for the oil. The foregoing would indicate that tilted fluid interfaces are not uncommon and that it would be of importance to underbtand the forces which are responsible. Further. it would be of interest to examine certain idealized models of situations which may cause tilting with the thought in mind of deriving relatinnships between the important variables. These relationships may only hold exactly for the model* assumed but the results may apply semi-quantitatively to actual formations and at least indicate trends. The discussion will be divided into the effect of static factors and the effect of dynamic factors. STATIC FACTORS It has already been mentioned that one reason given for tilted contacts is a variation in formation properties. Specifically this is a change in the pore size of the formation capilaries. If the internal surface of the reservoir is preferetially wet by water, capillarity will cause the water 10 rise to a greater extent above the free water table in a small capillary than in a large one. Consequently the water-oil contact will be higher in the lower permeability sands than in the higher permeability sand.;. The free water table may be defined as the position the oil-water interface would take if capillarity were absent. Since the phenomenon is the capillary rise of water in an oil environment, a derivation of this shold give a relationship between the variables. Fig. 1 is an idealized model of the situation with a uniform bore circular capillary ubstitttting for the tortuous non.circular formation pore. The upward force.; on the column of water will be equated to the forces downward to give the equilibrium situation. forces upward = force downward where r is the radius of the capillary in ceritimetera. a is the interfacial ten4on oil-water in dynes per centimeter, 0 is the

Jan 1, 1953

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 P. J. Closmann

By P. J. Closmann

Steam zone growth as a function of time has been calculated for the case of constant rate steam injection into a preheated reservoir. To simplify the calculation a linear temperature profile has been assumed in the cap and base rock at the start of steam injection. The results indicate that at early times augmentation of steam zone growth due to preheating should be greatest. At longer times the steam zone development becomes close to that calculated with no preheating. INTRODUCTION With the increasing application of thermal recovery processes to recover viscous oil, cyclic steam injection has become important in many large-scale projects. In this process, the cycle of steam injection followed by oil production is repeated a number of times. At the beginning of the second and later cycles, steam is injected into a reservoir that has already been heated but that has lost part of its heat both in produced fluids and by conductive heat loss away from the injection zone. In such a case, the temperature level of the injection zone and the temperature distribution of the surrounding rock will affect the growth of the steam zone developed during the subsequent steam injection. Knowledge of the size of steam zone developed is important in determining the amount of oil displaced and the extent of heating in the reservoir. It is also useful to be able to compute the size of the steam zone for cases where thief zones take most of the injected steam. This paper presents a fairly straightforward method of estimating the steam zone developed in a preheated formation, based on certain simplifying assumptions. Some cases of steam injection into a reservoir at its original temperature have already been considered elsewhere.1-3 Because it is difficult or impossible to obtain an accurate representation of the temperature distribution in the reservoir some time after initial heating has taken place, in this work, a linear temperature profile in the cap and base rock is assumed (Figs. 1A and 1B). The linear profile has the advantage of yielding heat flow equations that

Jan 1, 1969

By R. G. Grieves, G. Thodos

A method is presented for the accurate calculation of the cricondentherm and cricondenbar pressures of multicomponent hydrocarbon mixtures of known composition. The mixtures may contain six and quite possibly any number of components including paraffins, isoparaffins, olefins, acetylenes, naphthenes and aromatics. The approach is similur to that used for calculating critical pressures and cricondentherm and cricondenbar temperatures. The critical pressure, the normal boiling point and the approximate vapor pressure behavior of each component are all that art required. A stepwise culculation procedure is necessary for mixtures containing more than two components. From an analysis of 22 binary systems and 118 mixtures, the average deviation of calculated cricondentherm pressures from reported values is 2.4 per cent. For nine multicomponent mixtures the average deviation is 1 per cent. Considering 19 binary systems and 108 mixtures, the average deviation of calculated cricondenbar pressures from reported values is 1.7 per cent. For 15 multi-component mixtures, the average deviation is 2.2 per cent. INTRODUCTION A knowledge of the phase behavior in the critical region of multicomponent hydrocarbon mixtures is of value both in industrial processing operations and for the optimum operation of gas condensate reservoirs. Accurate methods of calculating the critical temperatures and critical pressures4,6,7,8,9,23 and cricondentherm and cricondenbar temperatures6p10 of multicomponent hydrocarbon mixtures are available in the literature. If the cricondentherm and cricondenbar pressures could be calculated with equal accuracy, the entire phase diagram of a multicomponent hydrocarbon mixture could be well approximated. The work of Etter and Kay6 is limited to systems containing the normal paraffins and has not been tested on systems containing a heavier component than heptane. In addition, the development of their multicomponent equations is based upon a limited number of mixtures containing from three to six components. Silverman and Thodos have considered systems containing both paraffinic and non-paraffinic hydrocarbons, but their correlation is limited to binary systems and is highly inaccurate for methane systems. Eilerts has done extensive work on the cricondenbar pressure; he has produced an excellent correlation for binary systems. However, his procedure for multicomponent mixtures is chiefly useful for highly complex mixtures requiring a knowledge of the vapor-liquid equilibrium behavior of the mixtures; he has not considered the cricondentherm pressure. The objective of this study was the development of an accurate and rapid method for the calculation of the cricondentherm and cricondenbar pressures of multicomponent mixtures containing all types of hydrocarbons and having a wide volatility range. The approach that was adopted is similar to that used by Grieves and Thodos for cricondentherm and cricondenbar temperaturesl0 and for critical pressures.9 However, methane systems had to be considered separately and a modified stepwise calculation procedure was utilized for the cricondertherm pressure. The correlations were developed in a manner similar to those for critical pressures and cricondentherm and cricondenbar temperatures. Based upon binary data reported in the literature it was observed that the ratios of cricondentherm and cricondenbar pressure to the pseudocritical pressure (molar average), pt/ppc and pp/ppc, respectively, in two-component systems depended upon the mole fraction of the low-boiling component and upon the diversity in properties of the two components. A dimensionless boiling-point parameter Tb,/Tb was chosen to represent the diversity in properties of the components. For a binary system, Ti,'is the molar average of the normal boiling points of the two components involved. For a multicomponent system it is either the molar average of the normal boiling points of all of the components involved or the molar average of the normal boiling point of the pure low-boiling component and of the atmospheric boiling point of the low-boiling-free mixture of

Jan 1, 1965

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