Search Documents

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

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

Jan 1, 1970

By B. P. Nutter, M. Lebourg, J. A. McAlister

The Multi-Flow Evaluator (MFE) is a new system of tools providing an original approach in drill-stem testing. It improves control during the testing operation and gives a more accurate evaluation of the fluid recovered while providing additional pressure information for reservoir analysis. The tools are operated entirely by up and down motion of the drill stem. This up and down operation provides a positive means of control and offers easily observed surface indications of tool operating position. An unlimited number of shut-in and flow periods may be taken with this tool while in the hole. The key to the success of the MFE system of tools is the safety seal packer. Until this development, operation of the tools by up and down motion only had proven unreliable. The success ratio of the MFE now exceeds the ratio achieved by conventional tools. The equipment includes a 2,750-cc chamber in which a representative sample of the flowing formation fluid is trapped at the end of the last flow period and brought to the surface under pressure. The sample can be evaluated at the wellsite or transferred under pressure for laboratory analysis. The sampling feature of the tool allows a sample to be obtained from a reservoir which has suffered minimum influence from production. A representative sample can be obtained for laboratory or empirical analysis by employing a testing technique to minimize drawdown during flow periods. Interpretation methods to take advantage of this additional information are presented and supported by actual field examples. INTRODUCTlON The Multi-Flow Evaluator (MFE) represents the most recent technical advance in formation testing. Drill-stem testing techniques prior to the MFE have been limited in three areas—(1) operaiion, (2) recovery analysis, and (3) reservoir analysis. The tools in the MFE system allow maximum control during testing and provide test information superior to that obtainable with conventional tools. The system of tools and the operation procedures are described to illustrate the mechanical advantages over present tools. Methods of pressure build-up and fluid

Jan 1, 1966

By A. M. Khan, I. Fatt

A new type of steady-state thermal conductivity apparatus is presented and described. Results of measurements on Berea sandstone, Solehofen limestone, Bandera sandstone and Bakelite are presented. The measurements were made in a two-hour period at temperatures up to 80C and axial stress to 125 bars. Results show that the thermoelectric heat pump and tbermistor-controlled electric heater can be used in a steady-state apparatus for measuring thermal conductivity of mck. Steady-state is reached within two hours. With Urethane foam insulation, not more than 3 per cent of heat is lost by radial flaw. The apparatus can be used as a comparator or for absolute thermal conductivity determinations. INTRODUCTION Calculations of heat flow through the earth, both for scientific purposes and engineering design, require thermal-conductivity data on rock. However, the high variability of rock often makes handbook thermal-conductivity data of little value. For engineering purposes in particular, it is desirable to have measured data on samples of the rock mass under study. Laboratory methods for measuring thermal conductivity of rock use either time-dependent1-3 or steady - state temperatures .4,5 The time -dependent (transient) techniques give thermal diffusivity from which thermal conductivity can be calculated if the specific heat and density of the sample is known. Most transient techniques require lengthy and involved sample preparation. Also, the accuracy tends to be low. Steady-state techniques are more exact and are the most commonly used. The steady-state methods, however, also have their difficulties. A long time is often required to reach steady state, troubles have been encountered in establishing good thermal contact between the heating and cooling elements and the rock sample, and unaccounted-for heat losses are difficult to prevent. In this paper we present a description of a new type of steady-state thermal conductivity apparatus. With it, measurements have been made in two hours at temperatures up to 80C and axial stress to 125 bars. Sample preparation is simple and no thermal conductivity standard is needed.6 APPARATUS Our steady-state apparatus is an improvement . over the conduction pile used by Zierfuss and van der vliet4 and Somerton.5 A linear heat flux is established through a conduction pile consisting of the rock sample and a standard material of known thermal conductivity in series. The thermal conductivity of the sample is calculated from the ratio of the temperature drop across the sample

Jan 1, 1966

By K. Leutwyler, H. L. Bigelow

Jan 1, 1966

By R. G. Bentsen, D. A. T. Donohue

The cyclic steam injection process has become the most widely applied and most successful thermal recovery technique in use today.1-4 Normally, steam stimulation is repeated several times during the life of a project. The problem is to select the optimum number and size of steam treatments so that the profits realized over the entire life of the project are maximum. Obviously, operating costs rise directly with the frequency and amount of steam injection. Operating efficiency, on the other hand, can be adversely affected if the amount and frequency of stimulation are insufficient. Thus arises a multistage decision problem; that is, a decision must be made at each sub-period (day, for example), as to whether steam should be injected, and, if so, how much. It is desirable that these decisions be made in such a way as to secure maximum net profit before taxes. This can be achieved by dynamic programming. To use dynamic programming, we must be able to predict not only the production rate-time curve resulting from a given stimulation, but also how changes in stimulation policy will affect this curve. In this study, the effect of steam injection on well performance is estimated on the basis of a simple model that takes into account both steady-state and transient flow. Oil production rate is calculated for both types of flow, and the larger of the two values is used. This allows the calculation to change to radial steady-state flow at such time as this mechanism becomes more effective than the transient mechanism. We propose to show how dynamic programming can be used to optimize the steam soak process with respect to net profit. For this, a mathematical representation of the steam soak technique is necessary. A further objective, therefore, is the formulation of a mathematical model capable of realistically simulating the physical process. Oil Production Response Model The advantage of using a mathematical model is its generality and ease of manipulation. However, the model must be an accurate representation of the physical process, for a solution derived from a model can be no better than the model itself. The question now arises as to how well the model should be expected to simulate reality. In general, the simplest model should be tried first, for additional accuracy is apt to require additional cost and time. Also, the search for a high degree of accuracy may yield only greater complexity and difficulty without yielding significantly better results. With this in mind, we kept our model as simple as possible, yet generally representative of the actual physical process. The number of state variables is limited to three, and stochastic variations in the parameters are ignored. As a consequence, use of this model should be restricted to cases in which the underlying assumptions are not severely violated. In

Jan 1, 1970

By K. P. Fournier

This report describes work on the problem of predicting oil recovery from a reservoir into which water is injected at a temperature higher than the reservoir temperature, taking into account effects of vis cosity-ratio reduction, heat loss and thermal expansion. It includes the derivation of the equations involved, the finite difference equations used to solve the partial differential equation which models the system, and the results obtained using the IBM 1620 and 7090-1401 computers. Figures and tables show present results of this study of recovery as a function of reservoir thickness and injection rate. For a possible reservoir hot water flood in which 1,000 BWPD at 250F are injected, an additional 5 per cent recovery of oil in place in a swept 1,000-ft-radius reservoir is predicted after injection of one pore volume of water. INTRODUCTION The problem of predicting oil recovery from the injection of hot water has been discussed by several researchers. l-6,19 In no case has the problem of predicting heat losses been rigorously incorporated into the recovery and displacement calculation problem. Willman et al. describe an approximate method of such treatment. 1 The calculation of heat losses in a reservoir and the corresponding temperature distribution while injecting a hot fluid has been attempted by several authors.7'8 In this report a method is presented to numerically predict the oil displacement by hot water in a radial system, taking into account the heat losses to adjacent strata, changes in viscosity ratio with temperature and the thermal-expansion effect for both oil and water. DERIVATION OF BASIC EQUATIONS We start with the familiar Buckley-Leverettg equation for a radial system:* This is sometimes referred to as the Lagrangian form of the displacement equation. If the injection water is at some temperature Ti which is above reservoir temperature TR, then lw, the flowing fraction which is water, is a function of both water saturation and fluid temperature at any given time and position in the reservoir. The temperature dependence of lw results from the temperature dependence of the viscosity ratio µo/µm of which lw is a function. Thus:

Jan 1, 1966

By N. R. Morrow, C. C. Harris

The experimental points which describe capillary pressure curves are determined at apparent equilibria which are observed after bydrodynamic flow has ceased. For most systems, the time required to obtain equalization of pressure throughout the discontinuous part of a phase is prohibitive. To permit experimental points to be described as equilibria, a model of capillary behavior is proposed where mass transfer is restricted to bulk fluid flow. Model capillary pressure curves follow if the path described by such points is independent of the rate at which the saturation was changed to attain a capillary pressure point. Aemodified suction potential technique is used to study cyclic relationships between capillary pressure and moisture content for a porous mass. The time taken to complete an experiment was greatly reduced by using small samples. INTRODUCTION Capillary retention of liquid by porous materials has been investigated in the fields of hydrology,l soil science,2 oil reservoir engineering,3 chemical engineering," soil mechanics,5 textiles: paper making: and building materials In studies of the immiscible displacement of one fluid by another within a porous bed, drainage columnsl and suction potential techniques2 have been used to obtain relationships between pressure deficiency and saturation (Fig. 1). Except where there is no hysteresis of contact angle and the solid is of simple geometry, such as a tube of uniform cross section, there is hysteresis in the relationship between capillary pressure and saturation.2 The relationship which has received most attention is displacement of fluid from an initially saturated bed (Fig. 1, Curve R2), the final condition being an irreducible minimum fluid saturation Swr. Imbibition (Fig. 1, Curve A), further desaturation (Fig. 1, Curve R), and intermediate scanning curves have been studied to a lesser but increasing extent.2 ,3,9-12 This paper first considers the nature of the experimental points tracing the capillary pressure curves with respect to the modes and rates of mass transfer which are operative during the course of measurement. There are clear indications that the experimental points which describe these curves are obtained at apparent equilibria which are observed when viscous fluid flow has ceased; and any further changes in the fluid distribution are the result of much slower mass transfer processes, such as diffusion. Unless stated otherwise, this discussion applies to a stable packing of equal, smooth, hydrophilic spheres supported by a suction plate5 with water as the wetting phase and air as the nonwetting phase. The discussion can also be applied, in many

Jan 1, 1966

By F. F. Skiba, P. J. Root

A mathematical model has been used to investigate crossflow effects in a stratified reservoir during an idealized displacement process. A process in which one incompressible fluid displaces completely another incompressible fluid of the same density and viscosity was considered. The reservoir was assumed to be composed of homogeneous layers which are discrete but interconnected. Approximate five-spot flow geometry was simulated by arranging two pie-shaped cylindrical (two-dimensional) wedges in series. Some of the common production methods rely, for their effectiveness, on the fact that the formation consists of isolated strata. This is particularly true of selective plugging and single-zone production-injection methods. The effectiveness of these techniques for several cases in which adjacent strata are in communication has been evaluated. The usual production-injection schemes for fluid injection processes involve using the same completion interval in both the production and the injection wells. The efficiency of this, as well as alternative productioz-injection procedures has been examined for some cases involving communicating strata. INTRODUCTION The performance of many fluid injection projects is marred by early breakthrough and by production of a high percentage of the displacing phase after breakthrough. This channeling effect, if it occurs in reservoirs which have large permeability variations, is often explained on the basis of a model composed of discrete, isolated strata. Corrective or preventive measures based on this same model are also common. The behavior of stratified models has been the subject of considerable study. An excellent review of the early literature was given by Seba.1 Studies of the depletion performance have since been made for continuous, unrestricted communication between strata2-4 and for production from one permeable bed in continuous but restricted communication with an adjacent permeable bed.5,6 Model studies of displacement processes in stratified systems with crossflow also have been reprted.7-12 Finally, methods which account for crossflow have been developed for computing displacement behavior in stratified systems.13,14 Most methods for controlling or preventing channeling involve some change in the interval open to production and/or injection. Although many variations of the selective completion problem have been studied in regard to the depletion of a reservoir, little attention has been given to the behavior of a displacement process under similar conditions. Sebal used a resistance network model to study the performance of a process in which one incompressible fluid displaces an identical fluid from a linear, stratified reservoir. He investigated the possibility of improving the vertical coverage by limiting the completion interval of the production well. This investigation deals with techniques in which selective completion of the production and/or the injection wells is used in an attempt to control channeling. At one extreme, selective plugging and single-zone isolation (assuming perfect execution) should be successful if the reservoir is actually made up of independent strata. At the other extreme, they can hardly be expected to be successful if the channeling is caused by a mechanism in which stratification effects are negligible even though there are extreme variations in the formation permeability and/or the injectivity profile of the wells is not uniform. By contrast, it is difficult to anticipate the effect that the amount of communication between the beds of a stratified system will have on the efficiency of selective completion techniques. The objective of this study is the quantitative determination of the effect of communication on the behavior of some particular stratified systems which will at least lead to qualitative conclusions.

Jan 1, 1966

By W. A. Burns

Vertical flow is an important mechanism in many petroleum reservoirs. Yet no adequate method has heretofore been proposed for determining the in-situ vertical permeability in a reasonable length of time. Core measurements of vertical permeability do not necessarily represent effective in-situ values. Cores may be representative of only localized pockets of high or low permeability, and total reliance on cores might cause one to overlook the presence of vertical fractures altogether. Conventional single-well tests, interference tests, and pulse tests commonly are used to determine horizontal transmissibility, khh/u, but generally are not designed to measure vertical trans-missibility or other groups containing vertical permeability, ku. A new well test described here provides a measure of in-situ vertical permeability as well as the previously measurable horizontal permeability. This "vertical test" determines the vertical diffusivity, kv/gcu, and horizontal transmissibility, khd1/u (where dl is the length of injection or production perforations). The test is basically a short-term interference or pulse test normally lasting less than a day. Conventional down-hole equipment is used between two isolated sets of perforations in the same well, as shown in Fig. 1. A "layer-cake" description of a reservoir can be obtained by running two or more vertical tests in the same well, each over different vertical intervals. Each test can be analyzed for the effective vertical and horizontal properties of the interval tested. In addition, a vertical well test can determine the effective- ness of a shale barrier or tight zone by testing across it. Field use to date of vertical testing has confirmed its utility in reservoir description. Mathematical Basis Interpretation of vertical tests for ku and kh is based on comparison of measured pressure response with computer-generated response. Pressure behavior during a vertical test in a homogeneous, anisotropic, infinite-acting reservoir (infinite in radial extent and bounded in both directions vertically) is given by the diffusion equation for slightly compressible fluids with negligible gravity effects.l A schematic diagram of the problem is presented in Fig. 2, and an analytic solution, based on the image

Jan 1, 1970

By A. Timur

Fluid flow properties of porous media have long been of interest in such varied disciplines as geology, geophysics, soil mechanics, and chemical, civil, and mechanical engineering. This interest has resulted in numerous models of porous media that have been proposed, tested, and found to be useful under, at best, only special circumstances. A critical review of most of these models is given by Scheidegger. 1 The application of each new tool to the study of porous media has helped to test some of the existing theories and to form new ideas based on the new parameters being measured. In the application of nuclear magnetic resonance (NMR) techniques to the study of the properties of fluids in porous media, the theoretical studies of Korringa et al.2 resulted in a model for the relaxation of spin polarization of protons in a hydrogenous fluid in a pore of a solid. Seevers3,4 verified this model, established a method for measuring the surface-to-volume-ratio distribution in a porous medium, and proposed a technique for determining specific permeability in sandstones. In the present investigation, the uses of NMR methods are extended and a simple model of porous media is developed from an analysis of pulsed NMR measurements. In this model, the pore spaces of a porous medium are divided into three groups on the basis of their surface-to-volume ratios. To evaluate predictions of the NMR model, laboratory measurements of spin-lattice relaxation time, porosity, permeability, and residual (irreducible) water saturation were conducted on more than 150 sandstone samples obtained from four oil fields in North America. Applications are discussed for estimating the volume of movable fluid, ?p, and the specific permeability, k, of sandstone samples. The former parameter, which can be considered as producible porosity, is expressed by ?p = [1 -- (Swr/100)]? , .... (1) where + is porosity in percent of bulk volume and Swr is the residual (irreducible) water saturation in percent of pore volume at a capillary pressure of 50 Psi. Fast, nondestructive laboratory methods are described for determining porosity, movable fluid, and permeability of sandstone samples from measurements with a pulsed NMR (spin-echo) apparatus. Procedures are given for adapting these methods to samples of irregular shape and small size. Theory General In the Korringa, Seevers, and Torrey2 (KST) model the observed decrease in the spin-lattice relaxation time, T1, of protons of a hydrogenous liquid contained in the pore spaces of a porous solid is attributed to an increase in the correlation time for the random motion of the water molecules and to the presence of paramagnetic centers at the liquid-solid interface. This KST model predicts the observed relaxation time of a hydrogenous liquid in a pore of a solid through

Jan 1, 1970

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

By H. Kazemi, G. W. Thomas, M. S. Seth

The double-porosity model of Warren and Root7 for examining pressure drawdown and buildup phenomena in naturally fractured reservoirs has been extended to intetpret interference test results. Both analytical and numerical solutions are presented. The practicality of some of the assumptions in the Warren-Root model was also investigated. The standard procedure for interpreting interference tests in naturally fractured reservoirs beretofore employed an equivalent homogeneous reservoir model. A comparison of our results with the equivalent homogeneous model indicates the latter is inadequate for early time responses. For late time responses, however, the simpler equivalent homogeneous model yields favorable results. Finally, the effects of well spacing and the magnitude of the interporosity flow parameter are briefly cited. INTRODUCTION Interference tests can be used to provide reservoir information not obtainable from conventional pressure drawdown and buildup tests. The technique involves measuring the response in one or more wells arising from a pressure perturbation introduced into another well. It has been. successfully applied to single-layer, homogeneous, unfractured, isotropic and anisotropic formations to arrive at values of well storage capacity and a measure of directional permeability.l The techniques of interference testing usually can be extended to heterogeneous systems, provided the permeability contrast between any existing rock layers is not excessive.9 In such cases, the practical interference effects in an observation well of an interconnecting heterogeneous reservoir can be predicted using an equivalent homogeneous model. The calculated properties are then average values. Russell and Prats2 have demonstrated the utility of this approach for a single-phase, two-layer system. For naturally fractured reservoirs, some different approaches usually are necessary. The major departure from homogeneous or nearly homogeneous reservoir treatment is to characterize the fractured system as an entity having "double-porosity"; namely, we suppose the existence of both primary and secondary porosity, where the former is contained within a systematic array of identical rectangular matrix blocks and the latter is contained within an orthogonal system of continuous, uniform fractures (see Fig. 1). It is sufficient to assume that the permeability of the rock matrix is so small that movement of reservoir fluids to the producing wells is entirely through the fracture network. This assumption has been validated by Barenblatt3 and .employed by others?-7 Warren and Root7 used the double-porosity model to develop an interpretation of pressure drawdown and pressure buildup tests in naturally fractured

Jan 1, 1970

By L. J. Health

Synthetic rock with predictable porosity and permeability has been Prepared from mixtures of sand, cement and water. Three series of mixes were investigated primarily for the relation between porosity and permeability for certain grain sizes and proportions. Synthetic rock prepared of 65 per cent large grains, 27 per cent small grains and 8 per cent Portland cement, gave measurable results ranging in porosity from 22.5 to 40 per cent and in permeability from 0.1 darcies to 6 darcies. This variation in porosity and permeability was caused by varying the amount of blending water. Drainage-cycle relative permeability characteristics of the synthetic rock were similar to those of natural reservoir rock. INTRODUCTION The fundamental behavior characteristics of fluids flowing through porous media have been described in the literature. Practical application of these flow characteristics to field conditions is too complicated except where assumptions are overly simplified. The use of dimensionally scaled models to simulate oil reservoirs has been described in the literature. 1-3 These and other papers have presented the theoretical and experimental justification for model design. Others have presented elements of model construction and their operation. In most investigations the porous media have consisted of either u ncon soli da ted sand, glass beads, broken glass or plastic-impregnated granular substances—materials in which the flow behavior is not identical to that in natural reservoir rock. The relative permeability curves for unconsoli-dated sands differ from those for consolidated sandstone. The effect of saturation history on relative permeability measurements is discussed by Geffen, et a1.4 Wyga15 has shown quite conclusively that a process of artificial cementation can be used to render un consolidated packs into synthetic sandstones having properties similar to those of natural rock. Many theoretical and experimental studies have been made in attempts to determine the structure and properties of unconsolidated sand, the most notable being by Naar and Wygal.6 Others have theorized and experimented with the fundamental characteristics of reservoir rocks. This study was conducted to determine if some general relationship could be established between the size of sand grains and the porosity and permeability in consolidated binary packs. This paper presents the results obtained by changing some of the factors which affect the porosity and permeability of synthetically prepared sandstone. In addition, drainage relative permeability curves are presented. EXPERIMENTAL PROCEDURE Mixtures of Portland cement with water and aggregate generally are designed to have certain characteristics, but essentially all are planned to be impervious to water or other liquids. Synthetic sandstone simulating oil reservoir rock, however, must be designed to have a given permeability (sometimes several darcies), a porosity which is primarily the effective porosity but quantitatively similar to natural rock, and other characteristics comparable to reservoir rock, such as wettability, pore geometry, tortuosity, etc. Unconsolidated ternary mixtures of spheres gave both a theoretically computed and an experimentally observed minimum porosity of about 25 per cent.7 By using a particle-distribution system, one-size particle packs had reproducible porosities in the reproducible range of 35 to 37 per cent.6 For model reservoir studies of the prototype system, a synthetic rock having a porosity of 25 per cent or less and a permeability of 2 darcies was required. The rock had to be uniform and competent enough to handle. Synthetic sandstone cores were prepared utilizing the technique developed by Wygal.5 Some slight variations in the procedure were incorporated. The sand was sieved through U. S. Standard sieves. A 18-20 sand denotes that the sand passed the Number

Jan 1, 1966

By P. van Meurs, C. van der Poel

From observations of the linear displacement of oil by water from a porous medium as visualized in transparent models, new insight into the mechanism of vircous fingering, as occurring in the case of unfavorable oil-water viscosity ratios, has been obtained. A mathematical description of the observed phenomena leads to formulas. expressing both oil production and pressure drop across the formation as a func-tion of cumulative water injection with the oil-water viscosity ratio as a parameter. Theoretical results are compared with those obtained from scaled laboratory experiments in sand-filled rubes in which water was injected at the one end and production was obtained from the other end. In the second part of the paper the same principles are used to calculate a more realistic production scheme in which production of oil and water is obtained from wells regularly distributed over the field. Results show satisfactory agreement with production data of an actual oil field. Application of the theory is extremely simple as only a limited number of readily available parameters have to be known and results are presented in an analytical form. INTRODUCTION When oil is displaced from a porous medium by water having a viscosity lower than the oil. the oil-water interface is essentially unstable and has the tendency to break up into what are called "fingers" or "streamers". This Viscous fingering has been noted before by Engelberts and Klinkenberg' and has recently been demonstrated in transparent models by van Meurs-. It is highly questionable whether it is at all possible to perform relative permeability measurements In the laboratory reproducing the same distribution of water and oil in the cores as Present in practice. Therefore, we based our descriptive theory on an idealized model which was "tailored" to the fingering mechanism as observed in transparent models under scaled conditions. This approach was the more attractive, as a highly simplified model had already yielded an adequate description of the process. Though we are aware our theoretical approach can be extended in several ways, it is given here in its simplest form to explain the fundamental principles and to avoid complicated formulas. SIMPLEST CASE OF A LINEAR DRIVE Observations and Theoretical Model In Fig. 1 a series of photographs shows how water penetrates into a formation saturated with an oil which is considerably more viscous than the water. The technique for obtaining these photographs, using transparent models, has already been described?. The fingers can be clearly observed and are seen to be mainly parallel to the flow direction. Note that the contours of these fingers are not smooth but show protrusions of water enclosing oil pockets.

By E. F. Johnson, V. O. Naumann, D. P. Bossler

A method is presented for calculating individual gas and oil or water and oil relative permeabilities from data obtained during a gas drive or a waterflood experiment performed on a linear porous body. The method has been tested and found both rapid and reliable for normal-sized core samples. INTRODUCTION Individual oil and gas or oil and water relative permeabilities are required for a number of reservoir engineering applications. Chief among these is the evaluation of oil displacement under conditions where gravitational effects are significant, such as a water drive or crestal gas injection in a steeply dipping oil reservoir. Numerous proposed methods of obtaining relative permeability data on reservoir core samples have been too tedious and time consuming for practical use, or have yielded questionable and sometimes inconsistent results. A method has been developed by which the individual relative permeability curves can be calculated from data collected during a displacement test. The method is based on sound theoretical considerations. Using this method, with a properly designed experimental procedure, relative permeability curves can be obtained using core samples of normal size (i.e., 2 to 3 in. in length and 1 to 2 in. in diameter) within a few days after receipt of the core. In a recent publication D. A. Efros1 describes an approach to the calculation of individual relative permeabilities that is based on the same theoretical considerations. We believe the approach described in the present paper is more adaptable to practical application than the method implied by Efros. In addition, comparisons with independently determined relative permeabilities are furnished to substantiate the reliability of the new method. DERIVATION Previously the theory of Buckley and Leverett as extended by Welge2 as been used to calculate the ratio of relative permeabilities. In the derivation which follows, this theory is further extended to permit calculation of the individual relative permeabilities. The theory assumes two conditions which must be achieved before the method is applicable. They are that the flow velocity be high enough to achieve what has been termed stabilized displacement,3 and that the flow velocity is constant at all cross sections of the linear porous body. In stabilized displacement the flowing pressure gradient is high compared with the capillary pressure difference between the flowing phases. The high pressure drop insures that the portion of the core in which capillary effects predominate will be compressed to a negligibly small fraction of the total pore space. The assumption of constant flow velocity at all

By R. E. Gilchrist, R. F. Nielsen

When a gas is displaced by mother in a porous medium, and there is a relative immobile liquid present. there is a transition zone in which the gas composition varies from essentially that of the original gas to that of the injected gas. The change in extent of this zone with time may be governed by a mixing process in the direction of flow and by the rate of establishment of gas-liquid equilibrium, if the gases are appreciably soluble in the liquid. The relative importance of these two factors in deterri~ining the growth of the transition zone is discussed, with appropriate equations and experimental data. INTRODUCTION The various solution methods of oil recovery involve a transfer of the components of a mixture, whether within a single phase or between a liquid and gas. The similarity of a natural porous medium to systems used in separation and chromatographic processes makes it seem reasonable that certain concepts used by chemical engineers in "mass transfer" might be applicable to petroleum reservoirs. One approach, the HETP (height equivalent to a theoretical plate), involves a distance or length of time over which equilibrium might be con. sidered as established. Another approach involves direct use of a rate of mixing or rate of establishment of equilibrium. The latter approach has been preferred in relation to production problems. DISCUSSION AND PROCEDURE When a gas or liquid in a porous medium is displaced by one of different composition but miscible with it there will be set up a transition zone in which the composition varies from cssentially that of the original fluid to essentially that of the injected fluid. If the fluid being displaced is a gas, and there is a relatively immobile liquid present, the transition zone may be

By D. R. Parrish, F. F. Craig

In some respects, forward combustion (alias: fie-flooding, underground combustion, in-situ combustion, dry combustion) is an efficient oil recovery method. For example, its displacement efficiency (as much as 95 percent in reservoirs with high porosity and high oil saturation) is good. In its use of heat, however, forward combustion is grossly inefficient. Most of the heat generated during forward combustion is wasted. About one-fifth of the heat is picked up by the incoming air and re-used (regenerated); the remainder is left behind the combustion zone and ultimately is lost. Compressed air is expensive, and heat generated from compressed air via underground combustion is too valuable to be thrown away. Improvement is needed. Since the early 1950's, improvements in forward combustion by addition of water injection have been discussed, particularly in the patent literature, with increasing frequency.'-'' These concepts are based on the fact that relatively inexpensive water can be used to pick up much of the heat that is ordinarily wasted during forward combustion, and transport it forward to a more useful place, thereby reducing the amount of air required to move a heat front through the reservoir. Several approaches have been proposed. Merriam and Squires1 taught that a combustible (fuel) gas should be injected along with the air and water. This, however, would give at least a tendency toward reverse combustion instead of forward combustion, and would reduce the velocity of the heat bank, defeating the primary purpose of the water injection. Pelzer2 described primarily a batchwise method of forward combustion and water flooding. He advocated burning part way through a reservoir with forward combustion, then switching from air to water injection alone to push the heat toward the producing wells. Though possible, this method has serious disadvantages, the major one being that, by the time water injection is started, much of the heat has been irretrievably lost to the overlying and underlying formations. We have described what is considered a practical Combination Of Forward Combustion And Water-flooding (COFCAW).3,4 Air and water are injected simultaneously (or alternately) after a small heat bank has been formed by forward combustion. Some of the injected fluid fills up the region behind the heat bank; the rest enters the heat bank. Upon encountering heated rock, the water is converted into steam, which then flows ahead of the combustion zone and displaces much of the oil in place. Only enough air is injected to provide heat for vaporization of the water and to compensate for reservoir heat losses. Since the size of the heat bank is minimized, heat losses also are minimized. Because of more efficient oil displacement ahead of the combustion zone, fuel consumption is reduced. Too, unlike a forward combustion zone that cannot move until all fuel at a point has been consumed, COFCAW leaves some

Jan 1, 1970

By J. R. Kyte, L. A. Rapoport

The material balance equation for partiai or complete 1:trter-drivc reser1,oir.s been re-arranged to include a pressure irrferferencr ternr. This pressure interfernce term was ohtained from the theory established hy Morinda. A .successful progranl deissed to solve the resrclritlg c,qrration in a few hours on a nleelitr .sizeti, digital compuier. No prior knowdelge of the properties of the eqrrifer is required in the solution. 7' ri~rl calculations areatle using various assumed con!binuiions of value of the constant expretsing the qquifer properties correspoinds to the best fit ohtained between observed and calculated field history. Details of the method of calculaion are illus trated by an example. INTKCDUCTICN A method 01' calculating the material balance equation for a partial water-drive reservoir completed in an aquifer has been described in the literature." Pressure interference i'rc,m other pools iocated in the same nquiucr is not taken into account in this solution. When present, the pres- sure interference effect must be included in order to make a complete performance analysis. A correction for this effect has been developed using the theory cstahlished by Mor-tada.' Manual calculation of the material balance for water-drive reservoirs is a formidable task. Such manual ca1culation when interference eflects are includcd would be impossible to completc in any reasonable time period. However, a medium-sized stored program, digltal conlputer can he used satisfactorily. The method is described in the following sections. THECKETICAL CONSIDEKATIOhS The material balance equation for a reservoir producing by partial or total water drive, as proposed by Van Everdingen, This relationship is linear for Z/a from which thc slopc, Z, and the intercept, N, can be obtained by the least-squares method from the data points. It is not necessary lo know the aquifer properties to uktain oplimum values tor 2nd IV. Various assumptions are made tor At,,: where 11 is production time, seconds (oilfield between the equal successive time intervals into which the history is divided. t,, h It, times the number of equal time periods from i.e.. the total dimen-sionless time interval from It is noted that the dinlensionless time interval, At,,, involves all the relevant aquifer properties. Accordingly, any assunlption of is an assumption of the combined aquifer properties which pertain to the solution of the equation. Deviations of the data points from a straight line arc calculated for each assumption of and the value of A, correspond ing to the minimunl deviation is chosen as the optimum solution. This method is described in Ref. 1 and 2. Thc material balance solution is usually made with an initial assumption that the reservoir being studied is located in an infinite aquifer from which there are no other withdrawals. Deviation of the points by amounts greater than can be attributed to inaccuracies in field data may be the result of: (1) boundary effects of a limited aquifer, (2) pressure interference from othcr pools completed in the .sarne aquifer, and (3) a combination of these effects. Differentiation of these effects can be made only on the basis of geologic data. Pressure interference from other fields will result in a pressure drop at the original oil-water boundary having an interference component, The measured instantaneous pressure drops at the oil-water

By L. B. Davidson

The petroleum literature contains many reports on the relative permeability properties of porous media. However, only recently have studies of relative permeability at elevated temperatures been made.1 -5 Each of these studies presents relative permeability data on water-oil systems in linear sandpacks or sandstone cores under isothermal conditions. The results of these investigations are dissimilar. Edmondsonl found that in Berea sandstone a temperature-dependent effect was present, while Shilobod,2 using the same material, found no temperature dependence after normalizing for initial water saturation. Workers at Texas A&M U.3-5 found that the ratio kw/ko was temperature dependent in an unconsolidated system. The differences in results of these studies are probably a consequence of differences in the material employed, in the experimental procedures followed, and in the properties of the porous systems used. This paper describes an experimental investigation of permeability ratio temperature dependence, conducted to clarify some of the results obtained by previous workers. Each experiment described in this paper involved the isothermal displacement of Chevron No. 15 white oil by either nitrogen, steam, or distilled water and displacement of water by nitrogen. All permeability ratios were calculated using the Welge procedure. Values of the water-oil permeability ratio were obtained over the temperature range 75OF to 540°F. Results indicate the ratio is temperature dependent at low water saturations, probably due to interfacial effects. However, as the saturation is increased, temperature dependence decreases until, above some minimum saturation, the ratio becomes insensitive to temperature. Previously published work indicates that the dependence reappears at very high water saturations, a region of the ratio curve not completely examined in the present study. It appears that water-oil permeability ratios are temperature dependent at the extremes of the ratio curve where changes in interfacial tension with temperature induce variations in residual saturations. Between the extremes, there is a region in which the permeability ratio is insensitive to temperature. The size and location of this region is probably a function of the characteristics of the porous system. Gas-oil permeability ratios calculated from nitrogen-white oil displacement experiments indicate a definite temperature dependence over the range 75OF to 500°F and over the entire gas saturation range studied. Evidence discussed in this report indicates that this dependence is a consequence of molecular slippage in the gas phase. This is an extension of the Klinkenberg effect6 to two-phase systems at high temperatures (the existence of slippage in two-phase systems at room temperature has been demonstrated7). Evidence for a similar source of temperature dependence in superheated steam-white oil experiments

Jan 1, 1970

By W. W. Krech, T. E. Perkins

As fractures are propagated through rocks, energy is absorbed near the extending crack tip. Apparent surface energies for several rocks have been measured by cleavage under dynamic con-ditions. At nominal crack velocities from 0.5 to 500 in (min, measurements showed that fractures propagated in discrete jumps. Calculated surface energies and moduli were relatively insensitive to nominal rate of cleavage. in another set of experiments, rocks were cleaved under high confining stresses. The rocks were submerged in low leak-off fluids which formed a filter cake on the freshly cleaved surfaces (similar to the hydraulic lracturing process). Apparent surface energies were increased substantially as the surrounding fluid pressure was increased. Moduli in bending increased significantly upon application of the first 1,000 psi but were insensitive to stress level at greater pressures. INTRODUCTION For almost 20 years, hydraulic fracturing processes have been utilized effectively to stimulate oil and gas wells. During this period, some process improvements have resulted from studies of fracture orientation, mechanics of fracturing, areas generated, conductivities of cracks, etc. Yet many questions remain concerning the conditions and pressures needed during fracture propagation. In this paper we will report additional studies of the mechanics of fracture extension. It was shown previously3 that large rock samples could be cleaved under controlled conditions so as to measure the apparent surface energy (that amount of energy absorbed per unit area of new surface created). In this paper we consider the effects of two additional factors on surface energies, viz.: (1) effect of cleavage rate and (2) effect of confining stress level. THEORY Cleavage experiments were conducted on rock samples similar to that illustrated in Fig. 1. Blocks of rock several inches wide, 2- to 3-in. thick, and up to 3 ft in length were grooved longitudinally with shallow guide slots. A crack was initiated and allowed to extend along the web as the top of the rock specimen was pulled (or pushed) apart. Auxiliary equipment permits the measurement of (1) force applied at the top, (2) separation at the top and (3) crack length. (Further experimental details will be given in the next section.) The rock beams created by the crack are considered ton be cantilever beams. The deflection (or separation of the rock beams) at any point is calculated 4,5 by the beam Eq. 1.

Jan 1, 1967