Search Documents

By W. C. Hauber

Techniques previously published to predict the production performance of a water flood when the mobility ratio is not unity have been primarily restricted to a five-spot well pattern. When other types of well patterns are considered, the prediction of performance has required either a model study or a complex numerical solution of a multiple-fluid system with a high-speed digital computer. Simple methods are developed in this report for predicting the areal sweep efficiency and injectivity history of a water flood for arbitrary well patterns and mobility ratios. The effect of an initial gas saturation can be incorporated into this method of prediction. With the layer superposition concepts published by Prats, et al.,1 these methods can be used to predict the production performance of a water flood in a reservoir having a wide range in permeability. A reasonably good agreement is shown by comparisons of the results calculated by this method with those obtained either by model studies or by complex numerical solution of multiple-fluid systems. Therefore, it has been concluded that a reasonably accurate method has been developed for predicting the production performance of a water flood employing an arbitrary type of well pattern for any given mobility ratio. A medium-speed digital computer can be used in applying this method. INTRODUCTION Production performance of multi-fluid five-spot water floods can be predicted by the techniques published by Prats, et al.1 This technique takes into account the effects of initial gas saturation, mobility ratio and vertical variations in horizontal permeability. A basic assumption in applying this technique is that unit displacement efficiency does not vary with time after the passage of the flood front. Prats, et al. demonstrated that this assumption is valid for normal waterflood applications with small mobility ratios. To extend Prats' technique to other types of well patterns has required either an experimental model study or a numerical solution of a multiple-fluid system to describe the performance of an individual homogeneous layer. Sheldon and Dougherty2 presented an exact numerical solution employing any arbitrary well pattern. This technique is complex to apply, requiring (1) numerical solution of the single-fluid system, (2) conformal mapping of the streamlines and iso-potential lines to form a new grid network, and (3) a numerical solution of the multi-fluid system on the new grid network. A method was developed in this paper to approximate areal sweep efficiency and injectivity as functions of time for an individual homogeneous layer for arbitrary well patterns and mobility ratios. Analytic expressions were derived to directly apply this method to five-spot and direct line drive well patterns. For other types of well patterns, this method is similar to the approach of Higgids and Leighton,6 in that it utilizes the stream tubes and iso-potential lines obtained by the solution of differential equations when the mobility ratio is unity. However, this method is easier to apply than the Higgins-Leighton method since it does not require knowledge of the variation between oil and water relative permeabilities and saturation. Using the relationship between areal sweep efficiency, injectivity and the volume of water injected into an individual homogeneous layer as obtained by the method presented in this paper, and Prats' technique for the superposition of individual layers, one can predict the production performance of any water flood. ASSUMPTIONS GENERAL Basic assumptions made in this report are that fluids are incompressible and the layers are horizontal, homogeneous and uniformly thick. It is assumed that as water is injected into the layer, three distinct regions are formed —a water bank, an oil bank, and an unflooded region. Each region is assumed to displace other regions in a pis-tonlike manner, with only one fluid flowing in each region, so that there is no change in unit displacement efficiency with time after the flood front has passed in the reservoir. Therefore, the residual oil saturation in the water bank and the oil saturations in the unflooded region are constant throughout the layer. The ratio (F) of the bulk volume of a layer occupied by the oil and water banks at any time before oil breakthrough, to the corresponding bulk volume of a layer occupied by the water bank, is given' by the equation BREAKTHROUGH AREAL SWEEP EFFICIENCY To predict the behavior of any type of waterflood pattern, one must know the breakthrough areal sweep efficiency Eb(Mc,o, F) for the pattern for any water:oil mo-

Jan 1, 1965

By T. D. Mueller, F. C. Miller, Jr. Earlougher R. C.

In analyzing pressure buildup tests for field wells producing both oil and gas, the common practice is to use a modification of single-phase flow theory. Validity of such an approximation has been demonstrated for single-well solution-gas-drive systems. This paper indicates that such approximations are also valid for two-well solution-gas'-drive systems. From this it is inferred that the technique can be used for multiple-well systems. A computer was used to simulate the behavior of a two-well solution-gas-drive reservoir to test the validity of the above type of analysis. Simulation results indicate that pressure buildup tests in such a system can be analyzed within engineering accuracy for formation permeability and pressure. A rule of thumb is given for estimating the length of time a well must be shut in for the pressure in a nearby producing well to increase significantly. To observe such an increase, the shut-in well would have to be left shut in much longer than normal. INTRODUCTION An important technique for obtaining data concerning a producing petroleum reservoir is the pressure buildup test. Such a test, when properly conducted and analyzed, provides information on the average reservoir pressure and the permeability in the major drainage area of a well. The theory upon which the analysis of a pressure buildup test is based assumes that the behavior of the fluid in the reservoir is adequately described by the diffusivity equation.1 Use of the diffusivity equation implies that a single fluid of small and constant compressibility is flowing. To analyze a pressure buildup test in situations which involve more than one fluid phase, the single-fluid analysis is extended. This paper verifies the extension of pressure buildup analyses to two-phase, two-well systems. Miller, Dyes and Hutchinson1 have presented a technique for analyzing pressure buildup tests in circular bounded reservoirs. Their method requires that the reservoir be producing at pseudo-steady state (constant pressure gradients) prior to shut-in. An alternate reservoir model, presented by Horner,2 assumes that the reservoir is infinite. The Homer model does not assume a Pseudo-steady state prior to shut-in, but the assumption that the reservoir is infinite implies that the average pressure can be determined only if the total production prior to shut-in is small compared to the total fluid originally present. Limitations imposed by the pseudo-steady state and infinite reservoir assumptions are avoided by the technique proposed by Matthews, Brons and Hazebroek.3 Their method can be used to calculate both the Permeability and the average pressure in bounded single-well systems which are not necessarily at steady state. They also suggested determining the average pressure of a multi-well reservoir producing from pseudo-steady state by volumetrically averaging the static pressures of the individual wells. Matthews and Lefkovits4 used numerical simulation techniques to verify that this method does produce adequate results for single-phase, multiple-well systems. Perrine5 proposed modifications to the single-phase theory so that pressure buildup tests in multiple-phase systems could be analyzed. He suggested replacing the single-phase mobility (k/µ) by the sum of the mobilities of the individual phases, and replacing the fluid compressibility by an average compressibility weighted by the saturations of the separate phases. By using numerical techniques to solve the two-phase flow problem for a single-well system, he concluded that this approximation gives results within engineering accuracy. More recently, Weller6 has shown that, for a gas-oil system with a single well at the center of a circular reservoir, pressure buildup tests can be analyzed adequately by using the modification proposed by Perrine. Feller's results also indicate that, as the gas saturation increases, this analysis becomes less accurate. To date, no evidence is available to show that the Perrine approximations can be combined with the Matthews-Brons-Hazebroek technique to calculate the average pressure of a multi-well,

By R. C. Earlougher, F. G. Miller, T. D. Mueller, H. J. Ramey

There are many studies of flow in radial systems that can be used to interpret unsteady rerervoir flow problems. Although solutions for systems of infinite extent can be used to generate solutions fur finite ow systems by superposition, application is tedious. In this paper a step is made toward simplifying calcu1ations of such solutions for finite flow systems. Superposition is used to produce a tabulu-lion of the dimension1e.s.s pressure drop function at several locations within a bounded square that has a well at its center. The square system provides a useful building block that may be used to generate flow behavior for any rectangular shape whose sides are in integral ratios. Values of the tabulated dimensionless pressure drop function are simply added to obtain the dimensionless pressure drop function for the desired rectangulm system. The rectangular system may contain any number of wells producing at any rates. Furthermore, the outer boundaries of the rectangular system may be closed (no-flow) or they may he at constant pressure. Mixed conditions also may be conyidered. Tables of the dimensionless pressure drop function for the square system are presented and various applications of the technique are illustrated. Introduction In 1945, van Everdingen and Hurst' published solutions for the problem of water influx into a cylindrical reservoir. Since this problem is mathematically identical with the depletion of a cylindrical reservoir with a well at the origin, the van Everdingen-Hurst solution may be used to study the depletion problem. In their analysis, they assumed that the fluid had a small, constant compressibility such that flow was governed by the diffusivity equation For a constant production rate q starting at time zero, van Everdingen and Hurst showed that the unsteady pressurc distribution for both finite and infinite systems could be expressed in terms of a dimensionless pressure 'Sabulations of the dimensionless pressure drop for a unit value of r, were provided by van Everdingen and Hurst,' and later by Chatas.' Others also presented values in graphical or tabular form." If the radius of the well becomes vanishingly small, r,+ O, the line source solution may be used for Eq. 2 when infinite systems are considered. Eqs. 5 and 6 are excellent approximations for Eq. 2 under certain conditions: In 1954, Matthews, Brons and HazebroekV emonstrated that solutions such as Eq. 5 can be superposed to generate the behavior of bounded geometric shapes; i.e., the behavior of a bounded single-well system can be calculated by adding together the pressure disturbances caused by the appropriate array of an infinite number of wells producing from an infinite system. These wells are referred to as image wells. Matthews, Brons and Hazebroek considered systems containing a single well producing at a constant rate. This superposition can be represented analytically as

Jan 1, 1969

By H. K. van Poollen, H. C. Bixel

A treatment is given of the transient pressure behavior of a well located at the center of a circular region surrounded by a radial discontinuity. On either side of the discontinuity, the values of permeability, viscosity, compressibility and porosity are uniform but may be different from those on the other side of the discontinuity. The results are obtained by solving a pair of finite difference equations. The numerical solution to these equations is obtained using a digital computer. The results are compared with previously published, approximated analytical solutions to the same problem. Solutions, presented graphically, show pressure decline for constant rate fluid production. The range of variables studied include dimensionless time from '0.001 to 100 and storage capacity ratio from 0.001 to 1,000. The well behaves as if it were in an infinite reservoir for dimensionless times less than 0.25. Reservoir properties near the well can be estimated in the usual manner. An overlay technique is used to match an experimental curve with one of the theoretical curves. It is possible to estimate the distance to a discontinuity by substituting the actual time t and the corresponding dimensionless time tD at which a match occurs into the equation a = (0.000264 tk1/tD , µ1 cI)I/2 where kI/ I µI cI is the diffusivity near the well, and may be estimated from data taken at early time. Several buildup curves are computed. These curves show that for early shut-in times, correct values for transmissibility are obtained from conventional analysis. However, erroneous values of static reservoir Pressure are obtained unless data at large shut-in times are taken. INTRODUCTION A mathematical treatment of the transient pressure response of a well located at the center of a region bounded by a circular discontinuity is given. Within a region (Fig. 20) the properties of both the rock and the fluid are considered to be constant on either side of the discontinuity. While these properties are considered constant, they may be different on opposite sides of the discontinuity. A discontinuity of this type could be a fluid-fluid contact or a sudden change in rock characteristics such as thickness, porosity or permeability. Analytical solutions to this problem are available. 1-3 However, they are so involved that they are of little practical use. An approximate solutionl is available but the range of times over which it is valid has not been specified. A numerical solution to a pair of finite difference equations is used to obtain the solution given in this paper. PREVIOUS WORK One of the first solutions to the problem was published by William Hurst.1 Hurst considered unsteady-state flow of fluids through two sands in series with different mobilities in each sand. He used Laplace transforms to obtain a solution for a single well located at the center of the circle enclosing the first of the two sands. In this case, the solution for the pressure change in Region I (Fig. 20) is2

By D. M. Beeson, G. D. Ortloff

Water-propelled banks of carbon dioxide recovered both high- and low-viscosity crude oil substantially in excess of that recovered by water flood in linear flow model experiments. The increase in oil recovery is explained by oil swelling and viscosity reduction. Additional oil was recovered by solution gas drive at the end of the carbon dioxide floods when the model systems were depres-sured. This study was an investigation of the use of a water-driven carbon dioxide bank to recover crude oil. Techniques involving the injection of carbon dioxide have aroused con-siderable interest in recent years among methods considered for increasing oil recoveries. A number of investigators have studied the process in the laboratory1-' and in the field.7,5 Laboratory studies referenced indicate that carbon dioxide in conjunction with water flooding recovers more oil than does water flooding alone. These re- ports diverge widely, however, in stating magnitude of improvement. We believed that a water-driven carbon dioxide bank would increase recovery of crude oils by two mechanisms: viscosity reduction and oil swelling. Due to reduction in the oil viscosity caused by solution of carbon dioxide, displacement efficiency should be improved. In addition, the residual oil should be highly swollen by the carbon dioxide held in solution. Thus, the final residual oil saturation, in terms of stock-tank oil, would be lower than in the absence of carbon dioxide. To test this hypothesis, experi-

By A. B. Dyes

In the operation of a water-driven reservoir, a free gas saturation can he established by maintaining production rates fast enough to cause the reservoir pressure to decline below the bubble point. The benefit of such a procedure on the displacement efficiency of the oil by water is illustrated for two types of sandstone samples from one reservoir. Those rock samples showing the poorest recovery to water drive in the absence of a free gas saturation give the most improvement in the presence of a free gas saturation. Employing one type of reservoir fluid (low .shrinkage) the benefit of evolved gas on oil displacement is calculated at several reservoir pressures below the bubble point. This demonstrates the presence of air optimum pressure for water displace,rnent which is .several hundred pounds per .square inch below the bubble point. Displacement of oil at the optimum pressures results in the removal of 7 to I? per cent more oil than would be obtained by displacement at the bubble point pressure. In some instance. this magnitude of increase may be worthy of consideration in establishing the maximum efficient rate of production from many water-driven fields. INTRODUCTION For many years the belief that the use of gas could be combined with that of water to result in a better recovery operation than the use of either of these phases alone has been prevalent in the oil industry. This is evidenced by the several field trials which have employed the injection of gas and water in combination.1.2.v In most cases the small amount of any additional recovery, in combination with the limited knowledge available on the older reservoirs, has left uncertain the influence of gas on recovery in these earlier attempts to combine the use of gas with water. Research studies in recent years have demonstrated through laboratory. floods the magnitude of the benefit of the presence of gas on the displacement of oil from sandstones by water.5,6,7,8,9,10,11 In addition, studies of the distributions of the phases which are present during the displacement process have developed concepts concerning the mechanism by which this benefit is obtained."." It now remains for the engineer to incorporate these findings into the general technology of reservoir operations. It is the purpose of this paper to illustrate that the production of a water-driven reservoir at pressures below the bubble point is one such operation whereby improved oil displacement efficiency can be obtained as a result of gas evolution. The establishment of a gas saturation by evolution of dissolved gas was one of the techniques used by several investigators in the laboratory studies. The possibility of obtaining the benefit of gas evolution in water-driven reservoirs by producing at pressures less than the bubble point has 3150 been pointed out previously." The desirability of this operation has undoubtedly

Jan 1, 1955

By H. J. Ramey

Conventional calculation of total system isothermal compressibility for a system containing a free gas phase involves, among other things, evaluation of the change of oil and gas formation volume factors and the gas in solution with pressure. Preferably, this information should be obtained from laboratory measurements made with particular oils and gases. Often, experimental measurements are not available. In this case, it is necessary to obtain pressure-volume-temperature relationships from general correlations such as those of Standing' for California oils. In order to speed estimates of compressibility, generalized plots have been prepared of the change of both oil formation volume factors and gas in solution, with pressure from Standing's correlations. A generalized plot for estimating the change in the two-phase (oil and gas) formation volume factor with pressure is also presented. Usually, the effect of gas dissolved in reservoir water upon the total system compressibility is neglected for gas saturated systems, due to the low solubility of gas in water. Results of this study indicate that the increase in total system compressibility caused by solution of gas in water is often as large as the compressibility of wafer, and can be magnitudes larger for low pressure systems. Generalized results for estimating the change of gas in solution in water with pressure are presented in tabular and graphical form. INTRODUCTION During the past decade, pressure build-up and drawdown techniques have gained an important place in reservoir engineering. Build-up and drawdown analyses are only two special applications of the broad field of transient fluid-flow theory. All solutions of transient fluid-flow problems contain a parameter called the total system isothermal compressibility. This property of fluids and porous rock is a measure of the change in volume of the fluid content of porous rock with a change in pressure, and it may vary considerably with pressure. Evaluation of total system isothermal compressibility is not difficult, but it is tedious and time-consuming. Often compressibilities are estimated roughly, or transient flow methods are neglected completely. The benefits of using accurate system compressibility in properly-executed build-up or drawdown analyses are: 1. Better planning of pressure build-ups may be achieved to avoid unnecessary loss of revenue due to excessively long shut-in periods, or to shut-in periods too short to yield useable data. 2. Better and more reliable estimates of static formation pressures for reserves estimates and rate performance estimates. 3. Reliable information for evaluation of well completion effectiveness, and planning and interpretation of well stimulation efforts. The purpose of this paper is to clarify the nature of the total system isothermal compressibility, and to present useful methods for estimation of compressibility, particularly for systems containing a gas phase. DEVELOPMENT Numerous publications have presented solutions to transient single-phase flow of slightly compressible fluids, stressing pressure build-up applications. In transient flow, a compressibility* term arises to permit volume content of fluids in porous rock to change as pressure changes. The basic nature of the compressibility term is usually taken for granted. Problems arise in practical applications of transient fluid theory because most published works consider only one flowing fluid—in an ideal porous system containing only one fluid. In 1956, Perrine2 presented an intuitive extension of single-phase flow pressure build-up methods to multiphase flow conditions. Later, Martin- established conditions under which Perrine's multiphase build-up method had a theoretical foundation. Perrine has shown that improper use of single-phase build-up analysis in certain multiphase flow situations can lead to gross errors in estimated static formation pressure, permeability and well condition. It is likely that, much pressure build-up data for oil wells should be analyzed on the basis of multiphase flow. For both single-phase and multiphase build-up analysis, the isothermal compressibility term in dimensionless time groups often should be interpreted as the total system compressibility. All real reservoirs contain one or more compressible fluid phases. In addition, rock compressibility can contribute in an important way to the total system compressibility. The proper total system compressibility expression may contain terms for compressibility of oil, gas, water, reservoir rock and terms for the change of solubility of gas in liquid phases.

Jan 1, 1965

By L. H. Reiss, J. Groult, L. Montadert

Many fields, where the reservoir is composed of sandy layers, show great complexity because of the lack of continuity which results from a particular type of seditnentation. This complexity may be a factor of prime importance in studying reservoir drainage, especially in the displacement of oil by water or gas. This paper describes work done to study this heterogeneity and its influence in the exploitation of the fields of the lower Devonian sandy shale in the Fort Polignac Basin in Algeria. Two lines of research were followed. 1. On the border of the hasin, a detailed sedimentological study established rules for the distribution of the sandy zones, which furnished hypotheses necessary for devising a geologic model of the reservoir. A certain number of sandy zones were recognized and characterized by their sedimentation. 2. In the interior of the basin, production logs deterrnined the productivity profile. Study of these point data led to a description of the distribution of oil and gas phases in the reservoir so that the continuity of different zones between the wells could be deduced. These techniques were applied to the study of the lower Devonian sandy shale reservoir of Zarzaitine. The sedimentological study guiding the correlation work, with the aid of the conventional static data (logs and cores), led to the establishment of a useful geologic model of the reservoir. Production logs, thanks to the contribution of dynamic elements, permitted corrobora-tion of the model. In particular, the logs confirmed separation, as proposed by the sedimenfological study, of two comparable sandy groups by revealing the continuity of a thin impermeable layer, invisible on conventional logs. INTRODUCTION The heterogeneity of sandstone reservoirs considerably influences the distribution of fluids, as well as their displacement in the reservoir. Rational development of a field and maximum recovery of hydrocarbons are closely tied to an exact and precise knowledge of the intern architecture. We will try to show how it is possible I refine the knowledge of the reservoir by seeking qui different and varied sources of information. Productic geology, the group of techniques which leads to idei concerning reservoir architecture and its geometric descris tion, furnishes the framework of this undertaking. Of these techniques. those used so far are qui varied. Descriptive as well as quantitative, the technique partake as much of classical geology as of reservo science. They include (1) geologic studies of the enti basin, especially a stratigraphic and sedimentologic study of the reservoir member at its outcrop, detail1 observations of cores taken from wells and, finall oriented cores; (2) interpretation of geophysical dat (3) interpretation of electric logs, nuclear and sonic lo and dipmeters: (4) petrophysical study of the reserve on samples taken from wells; (5) reservoir dynamics (stu of the production history and, in particular, interpretati of production logs); and (6) analysis of the fluids a thermodynamic studies. This incomplete list could also include (1) study the migration of the fluids, (2) geochemistry (chemi characterization of the oils), (3) regional hydrodynamii and (4) injection and detection of tracers. If most of these techniques are classic, attention i nevertheless be called to applying two of the techniql in conjunction: the study of the reservoir member its outcrop and the study of dynamic distribution fluid phases in the reservoir through use of product logs, used in the study of the sandy shale reservoirs the lower Devonian of the Fort Polignac Basin, Alge in particular, their application to the study of the 2 zaitine field where they have led to a new and satisfact model of the reservoir will be described. THE ZARZAITlNE FIELD Situated in the interior of the Fort Polignac 8 about 150 kilometers (km) north of the lower Devol outcrop, the Zarzaitine Devonian oil field was discovc On Ian. 28, 19'' (Figs. 1 and l8). It is in the fern a monocline truncated by two perpendicular fault: great throw, with a downthrown section to the north\: Erosion at the top of the structure has caused one of it to disappear. It is a reservoir with a gas cap an m coefficient on the order of 1:5, open to the aql On the northeast, with the oil zone covering an Of 106 sq km and with initial oil reserves of 222 mi

Jan 1, 1967

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