Search Documents

By L. A. Rapoport

This paper is intended as an aid in the perfornzance and interpretation of experimental studies of multi-phase flow in porous tnedia. The mathenmatical formulation of incompressible, two-phase flow phenomena in a three-dimensional porous medium is presented. The resulting equations which account for the effects of gravitational and capillary forces are used to derive general scaling laws permitting accurate representation of water-oil displacement processes by means of experimental flow models. The significance of the scaling laws is examined in the light of the physical concepts of micro-behavior of fluids in porous media. Formulations of scaling con(litions obtained by application of dimensional analysis are reviewed. The principles of application and the limitations of flow model studies are discussed. INTRODUCTION Most of the problems facing a reservoir engineer involve three-dimensional, two-phase flow systems. influenced by the effects of gravitational and capillary forces. Present computational procedures generally do not permit the evaluation of such systems without introducing greatly simplifying assumptions. The actual reservoir configuration is usually replaced by one or a combination of several "equivalent" linear systems, and the effects of capillarity and gravity are considered to be negligible. While these assumptions may in many instances be justified, there are a number of important practical problems for which the results obtained on the basis of such simplifications remain open to question. An alternate approach to the problems that are not amenable to analytical treatment consists of seeking their solution by experimental laboratory investigations. The main difficulty of such investigations lies in the fact that multiphase flow phenomena subjected to gravity and capillary forces proceed differently at different rates and in systems of different dimensions. Therefore, unless great discrimination is exercised in selecting the conditions and interpreting the results of laboratory tests, such tests may be entirely misleading in regard to field applications. The purpose of this paper is to present a theoretical background for the performance of laboratory tests that are more nearly representative of field behavior. It will be shown that simulation of field behavior can be achieved by conducting laboratory tests under conditions selected in accordance with definite scaling laws. The formulation of conditions required for the scaling of two-phase flow experiments in porous media has been presented in the literature on several occasions.12,3,4 One of these formulations advanced by Leverett1, cannot be considered applicable to most practical cases inasmuch as it does not cover generally the regime of laminar flow. An adequate treatment, based on the consideration of dimensional analysis, was presented by Engelberts and Klinkenberg.2 Their treatment, however, is concerned primarily with the scaling of linear flow systems, and hence somewhat limited in scope. Accordingly, the formulation of more general scaling laws, applicable to three-dimensional systems, appeared to be warranted. It appeared, furthermore. desirable to establish this formulation on the basis of a mathematical treatment, rather than by application of general principles of dimensional analysis, to permit more explicit

Jan 1, 1956

By K. Leutwyler, H. L. Bigelow

Jan 1, 1966

By H. E. Ellis, O. L. Patterson, Granville Dutton

A large-scale water-injection program for the Holt Bryant reservoir of the Delhi field has been established by means of the new automatic multi-pool analyzer. The objectives of this program were to increase reservoir pressure and eliminate a cross-dip gradient of 200 psi. Pressure increases during the first 15 months of operation have allowed for increased oil production rates of 2,600 BID or 25 per cent. To determine the required water injection it was necessary to predict natural water influx ram. Elimination of the pressure gradients required the water-injection program to account for the fluid migration and pressure interaction between segments of the reservoir. An integral part of the ana1ysis included the determination of in-place volumes of oil and gas. The operation of the automatic multi-pool analyzer in determining the water-injection program is described. New features of the analyzer inc1ude an automatic maferial balance analog and a leart-squares computer for determination of unknown reservoir properties. Operational procedures are sufficiently simple to allow direct use by reservoir engineers on a routine basis. The rapidity and accuracy with which problems may be solved permit repeated ana1yses as a check upon reservoir performance. INTRODUCTION This paper describes the use of a new analog com- puter in determining the optimum injection program for a complex reservoir. The Holt-Bryant (Fig. 1) is the major reservoir in the Delhi field of Louisiana. Production is from the Tuscaloosa at subsea depths of 2,970 to 3,286 ft. The sand dips to the south and southwest from three to five degrees. The reservoir is approximately 12 miles long and varies in width from 1/2 to 2 miles. Various calculations have placed the initial volume of oil in-place in the wide range of 275 to 355 million STB. Two small gas caps were initially present in the reservoir. The Holt-Bryant reservoir was discovered in Dec.. 1944, and was fully developed with 188 producing wells by 1948. After development was completed, an east-west pressure gradient began to establish itself. This

Jan 1, 1957

By B. C. Craft, Murray F. Hawkins, John R. Rayne

A reservoir mechanism of sulfur recovery by the Frasch process is presented. Improving the economic, of recovery appears to be largely a well, rather than a reservoir problem. A most important factor is the limitation of the lateral extraction of .sulfur about the wells due to the almost vertical .flow of the injected hot water. Model studies are described which confirin the mechanism of sulfur recovery by presently used methods. Studies on a five-spot system indicate improved econo171ics regarding recovery, recovery rate, therinal efficiency, and well spacing. INTRODUCTION Sulfur mining by the Frasch process is one of the most important methods of sulfur production. Developed by Herman Frasch in 1894, it enjoyed a rapid growth and presently supplies 40 per cent of the world's sulfur and about 80 per cent of the U. S. sulfur." The Frasch process of mining the underground native sulfur deposits involves melting the sulfur in place by introducing super-heated water into the sulfur bearing formations, and producing the sulfur as a liquid. Many problems are encountered in mining these native sulfur deposits, and sulfur production depends chiefly upon the nature of the sulfur deposit and the underground thermal efficiency of the heat transfer system. THE RECOVERY MECHANISM Sulfur is typically found in salt dome cap rock as crystal aggregates occurring within the cavities, seams, and pores of limestone and gypsum formations.' Sulfur-bearing limestone is typically cavernous, vugular, and fractured. Porosity averages about 20 per cent and sulfur saturation averages approximately 25 per cent by weight.1,15 Vhe sulfur bearing formations are very heterogeneous in nature and physical characteristics vary widely within small distances. Total cap rock thickness may vary from 50 to more than 1,000 ft, and net productive thicknesses from a few feet to a few hundred feet. In most cases the upper portion of the limestone formation is barren of sulfur, varying from 5 to 200 ft in thickness.1,15,33 The immediately overlying sediments and the underlying dense anhydrite formation provide permeability barriers to the flow of water from the cap rock section.'" Sulfur wells are drilled with conventional oilfield drilling rigs. Fig. 1 illustrates the casing program of a completed sulfur well. The cap rock sections are drilled with core bits and continuous core logs are recorded. After coring, the cap rock section is underreamed to a diameter of 9 7/8 in. The bottom of the 8-in. casing is welded closed and rests on the bottom of the drill hole, which extends through the cap rock to a few feet above the bottom of the sulfur-bearing formation. The 8-in. pipe contains 5 ft of 3/4-in. perforations on 6-in. centers and 15 ft of 1-in. perforations on 1-ft centers. The 20-ft perforated section is set opposite the lower 20 ft of the sulfur-bearing formation. The 4-in. sulfur pipe is set inside the 8-in. water line on a seat ring at the top of the 5-ft section of 3/4-in. perforations. The mining water, heated to about 330°F under a pressure of 250 psi, is forced down the annulus between the 8-in. and 441-1. pipes and through the 1-in. perforations to the sulfur-bearing formation. Compressed air, at 400 psi pressure, is supplied through a 1-in. air line near the bottom of the well to provide air-lift of the liquid

Jan 1, 1957

By A. Lubinski, C. R. Stewart, K. A. Blenkarn

Storage porosity has been considered one of the important pore geometry characteristics of heterogeneous-porosity limestones. Storage pores are only containers for fluids, in contrast to flow channel pores which both contain fluids and provide continuity for fluid flow. The concept of another geometry characteristic, "poro-striction", is presented as a second pertinent variable in describing limestone pore space. In simple terms, poro-striction is a measure of the flow resistance between storage and flow channel pores. Alternating-flow core-testing theory provides flow relationships which can be used to divide the pore space of heterogeneous-porosity media into flow channel and storage pores and to measure the "porostriction" of the latter. Experimental application of this theory to naturally occurring heterogeneous limestones shows that "porostriction" and the ratio of storage to flow channel pores can be estimated. Porostriction and porosity ratio are microscopic characteristics which should influence oil-recovery efficiencies during certain types of displacement processes. A knowledge of pore geometry should be valuable in designing or selecting the most effective oil-recovery process for heterogeneous limestones containing storage porosity. INTRODUCTION The primary and secondary recovery of oil from limestone (carbonate) reservoirs has been and will continue to be a major source of the world's supply of petroleum. Oil recovery from this type of formation is strongly influenced by the size, shape and arrangement of the pore spaces which hold the formation fluids. Therefore, it is expected that better knowledge of pore geometry would lead to design of more effective recovery processes. The concept of storage and flow channel pores has been presented as an important characteristic of limestone porosity.1 In this concept, storage pores are only containers for fluids, whereas flow channel pores are those pores which both contain and transmit the fluids. The work reported herein involves the concept of an additional pore characteristic called "porostriction". This characteristic can be visualized as the flow resistance between the storage and flow channel pores. This resistance is an important pore geometry characteristic because it determines how easily gas and oil can replenish fluids removed from the flow channel pores. Fatt and co-workers at the U. of California have studied the effect of storage porosity and a factor similar to porostriction on transient-flow response of models of porous media during a pressure build-up or drawdown test.55 These models consisted of porous sandstone cores with artificially simulated storage pores. It is the purpose of this paper to describe tests which permit the storage porosity and porostriction of actual reservoir core samples to be estimated. The tests used for this purpose involved not transient tests but, rather, an alternating-flow technique in which pressure and flow are varied according to a sine-wave relationship. This technique may offer advantages over transient methods. THE POROSTRICTION CONCEPT Before discussing alternating-flow (A.F.) tests in core samples, it is necessary to define the porostriction parameter. For this purpose attention is focused upon an elemental volume v representative of a porous medium having both flow channel pores and restricted-access storage pores. For such an elemental volume, it is considered that the storage porosity is combined to form a single pore and that the access into the storage porosity may be represented by a single capillary tube connecting the single pore to the flow channels of the elemental volume. Let the radius of the capillary tube be designated as Rc and its length as L,. When fluid flows into or out of the pore, the following equation (Poiseuille's law) relates flow rate q to pressure drop p, and the physical characteristics of the tube and the fluid of viscosity jx.

By G. A. Croes, N. Schwarz

The dimensionless groups, to which the variables that govern the displacement of oil from reservoirs by liquids can be combined, are derived. Three types of displacement are considered, viz. cold-water drive, hot-water drive and solvent injection. The derivation of the dimerzsionless groups is carried out by means of the relevant basic equations (in-spectional analysis). The resulting sets of groups are afterwards completed by means of dimensional analysis. The form of the groups is given in such a way that they can be adapted to suit the various boundary conditions that are encountered in practice. The physical meaning of the groups is discussed. They have all been brought together on a chart, from which their mutual relation and their corresporzdence to related dimensionless groups in common use in other engineering sciences can be read off. The limitations of dimensionally scaled model experiments as a useful tool for studying liquid flow in porous media, as occurring in oil reservoirs, are discussed. The use of models to study fluid mechanics has an appeal for everyone endowed with natural curiosity. What active boy has not played with ship and airplane models, or crude models of dam and drainage systems Even in the most advanced technical engineering, such models play a fundamental and indispensable role. "And yet in few departments of the physical sciences is there a wider gap between theory and practice, between scientific knowledge and the state of art, than in the use of models to study hydrodynamic phenomena." Garrett Birkhoff1 INTRODUCTION Laboratory displacement experiments are extensively used to investigate, directly or indirectly, the production behavior of petroleum reservoirs. Such experiments are representative of the reservoirs as a whole, if they are carried out with models that are "properly scaled." The performance of oil reservoirs is governed by the values of a number of variables, which can be combined to dimensionless groups. Two general methods are available for the derivation of these groups. In the first of these, dimensional analysis, combination of the variables is done essentially by trial and error. Its only premise is knowledge of the complete set of relevant variables. In the second method, inspectional analysis, the dimensional homogeneity of the equations describing the behavior of the system to be studied is used. The resulting dimensionless groups can be divided into: (a) independent (variable) groups, such as dimensionless length, time, etc.; (b) dependent groups, giving the dimensionless form of the variables, such as recovery and pressure, that can, at least in principle, be measured during an experiment; (c) similarity groups, which are independent constant groups, such as ratio of length to height of the reservoir and ratio of

Jan 1, 1957

By O. C. Baptist

The Umiat oil field is in Naval Petroleum Reserve No. 4 between the Brooks Range and Arctic Ocean in far-northern Alaska. The Umiat anticline has been tested by 11 wells, six of which produced oil ; however, [lie productive capacity and recoverable reserves of the field are subject to considerable speculation because of unusual reservoir conditions and because several wells appear to have been .seriously damaged during drilling and completion. Oil is produced at depths of 275 to 1,100 ft; the depth to the bottom of the permanently frozen zone varies from about 800 to 1,100 ft, .so that most of the oil reserves are in the permafrost Reservoir pressures are estimated to range from 50 to 350 psi, increasing with depth, and the small amount of gas dissolved in the oil is the major source of energy for production. Laboratory tests were made on cores under simulated permafrost conditions to estimate oil recoverable by solution-gas expansion from low saturation pressures. The cores were also tested for clay content and susceptibility to productivity impaiment by swelling clays and increased water. content if exposed to fresh water. The results indicate that oil can be produced fronz reservoir rocks in the permafrost and that substantial amounts of oil can be produced from depletion-drive reservoirs by a pre.s.r~lrr drop of as little as 100 psi below the saturation pressure. Freezing of formation water reduces oil productivity much more than that due to increased oil viscosity: Failure of we1ls drilled with rtuter-base mud to produce is attributed to freezing of water in the urea immediately surrounding the wellbore. Swelling clays apparently contributed very little to the plugging of the wells. INTRODUCTION Naval Petroleum Reserve No. 4 lies between the Brooks Range and the Arctic Ocean in northern Alaska. The Umiat oil field is located in the southeastern part of the Reserve and is about 180 miles southeast of Point Barrow (the only permanent settlement in the Reserve and the primary supply point for drilling of the wells at Umiat). Eleven wells were drilled for the U. S. Department of the Navy, Office of Naval Petroleum and Oil Shale Reserves, between 1944 and 1953 to test the oil and gas possibilities of the Umiat anticline. Six of these wells produced oil in varying quantities and the best one pumped about 400 B/D.' Estimates of recoverable oil range from 30 to 100 million bbl. The main oil-producing zones are two marine sandstone beds in the Grandstand formation of Cretaceous age: these are referred to as the upper and lower sands. Good oil shows were found throughout the sand settions in the first three wells drilled on the structure, but the highest rate of oil production obtained on any 01 the many tests was about 24 BOPD. These first wells were drilled with conventional rotary methods using water-base mud; later wells were drilled either with cablc tools using brine or rotary tools using oil or oil-base mud. These experiments were successful as is shown by comparing the oil production from Well No. 2 with that from No. 5. These two wells are only 200 ft apart and are located at about the same elevation on the structure. Well No. 2. drilled with a rotary rig using water-base mud, was abandoned as a dry hole after all formation tests were negative. Well No. 5. drilled with cable tools and reamed with a rotary using oil, pumped 400 BOPD which was the maximum capacity of the pump and less than the capacity of the well. These field results indicated that the producing sands were extremely "water sensitive" and it was assumed that the cause of this sensitivity was the presence of swelling clays in the sands. Because of the very unusual reservoir conditions and the difficulties encountered in completing oil wells in the permafrost. the Navy asked the U. S. Bureau of Mines to make laboratory studies under simulated permafrost conditions to assist them in estimating the production potential of the field and the recoverable reserves. These tests were designed to determine the cause of the plugging of wells in the permafrost and to test oil recovery from frozen sand by solution-gas expansion with the oil gas-saturated at very low pressures. EXPERIMENTAL METHODS AND PROCEDURES Samples Analyzed Core samples were analyzed that represent the lower sand in Umiat Well No. 7, the upper sand in No. 3. and both the upper and lower sands in No. 9. These sands should be productive in all of the wells because of their location on the structure. Core samples from

By Walter J. Karplus

By means of a finite difference expansion a fluid flow field in cylindrical coordinates with axial symmetry, is simulated by a network of electrical resistors. A series of DC analog computing units, connected to the network bounrlary corresponding to the free surface effectively adjust the network parameterr to satisfy simultaneously Laplace's equation and the borrndary conditions. The system automatica1ly "re1axes" to the correct boundary shape in a fraction of a second. A detailed application of the method to the water coning problem in oil wells is presented. INTRODUCTION The water coning problem in petroleum production is a prime example of a class of problems known as free surface problems. A free surface arises in fluid flow regions when the shape of the boundary of a single phase flow region is affected in some manner by the pressure or potential gradients within the region. When there exists an interface of two fluids of different densities, the shape of this interface will be determined by the relative velocities of the two fluids and will therefore comprise a free surface. Other free surface problems of interest to petroleum engineers include the gas coning, water injection, and gravity drainage probIcms. In each case it is desired to predict the shape of the boundary of the oil zone as a function of the flow rate, the well penetration and the physical parameters and dimensions of the system. Such problems do not lend themselves to the classical analytical treatments of boundary value problems because in this case the boundary of the potential field and the potential distribution within the field are interdependent and must be calculated simultaneously. A further difficulty arises from geometric considerations. Since the fluid flow is of the radial type in most cases, such effective methods as conformal transformation and the hodograph are not applicable. It is the purpose of this paper to present a method whereby free surface problems of this type may be solved rapidly to the desired order of accuracy by means of a combination of a network of electrical resistors and conventional analog computer units. The water coning problem is treated in some detail. Then the procedure for extending the technique to other free surface problems is indicated. FORMULATION Most commercial oil reservoirs consist of an oil-saturated sand, bounded above by an impermeable strata and below by a water-saturated sand. Gravimetric separation results in a horizontal oil-water interface. When an oil well is drilled so that it partially penetrates the oil zone, and production is initiated at a constant rate, the oil-water interface will assume a new shape as a result of the pressure gradients in the oil zone. It should be emphasized that this is a steady phenomenon. After an initial transient period the water phase will be completely static and only the oil will be in motion. The new steady-state shape of the interface, as shown in Fig. 1, depends upon the depth of penetration, d, the thickness of the oil zone, I, the production rate, q, the hydraulic permeability zone of the sand, k, the viscosity of the oil, µ, and the densities of the oil and water, p.,. and pw.. The problem is to determine this relationship. Darcy's law + = - (P f g p z).....(1) it.

Jan 1, 1957

By Morris Muskat, J. M. McDowell

In both production and refining operations of the oil industry many processes are controlled by the gas-liquid phase relationships of the hydrocarbon mixtures of interest. The quantitative behavior of fractionating or distillation columns and stills depends on the changes in the equilibrium distributions of the hydrocarbon components between the gas and liquid phases as the pressure or temperature conditions may vary along the length of the column or as the gross operating parameters and compositions are varied. In oil producing operations the nature and amount of gas phase developed within an underground reservoir as it is being depleted and its pressure declines also involve the basic equilibrium gas-liquid phase interactions. The influence of surface conditions of temperature and pressure and various separation procedures on the nature and amount of stock tank oil and natural gas recovered from a given well stream is likewise determined by the same basic phase equilibrium characteristics of the hydrocarbon systems. Such equilibrium separation of the gas and liquid phases in a well stream is of importance both in establishing optimum separator conditions for obtaining maximum stock tank oil yields and in the general problem of crude stabilization. The phenomena involved in these problems and their quantitative aspects can, of course, be established in each individual instance by appropriate laboratory experimentation. It has been found, however, that by associating with the individual hydrocarbon components functional characteristics de- scribing their individual gas - liquid phase equilibrium behavior, the phase properties of the composite mixtures can be predicted by mathematical computation. These characteristics have been termed "equilibrium constants", or "equilibrium ratios"," and are defined by the general equation: where K, is the equilibrium ratio for the jth component, y, is the rnol fraction of that component in the gas phase, and xj is the corresponding mol fraction in the coexisting liquid phase. Considering the Kj as representing available and measurable properties of the individual components of a mixture, and as defined by Eq. (l), it is readily shown that in a coexisting gas-liquid system the rnol fraction concentrations of the individual components in the gas phase, y,, and those in the liquid phase, x,, are given by the expressions: where the nj are the corresponding rnol fraction concentrations in the composite system, and ng is the mol fraction of the whole which is in the gas phase. The basic unknown in Eq. (2) is n,: which gives the gross separation between the gas and liquid phases. This may be determined in principle by imposing the requirement that the sums of the mol fractions in both the gas and liquid phases are unity, i.e. by the equations: The solution of the equivalent Eqs. (3) lying between 0 and 1, when inserted into Eqs. (2), provides a complete description of the compositions of the coexisting gas and liquid phases. While the Eqs. (3) are essentially equivalent to polynomial equations, such transformations are of little value when dealing with multicomponent sys-

Jan 1, 1949

By J. M. McDowell, Morris Muskat

In both production and refining operations of the oil industry many processes are controlled by the gas-liquid phase relationships of the hydrocarbon mixtures of interest. The quantitative behavior of fractionating or distillation columns and stills depends on the changes in the equilibrium distributions of the hydrocarbon components between the gas and liquid phases as the pressure or temperature conditions may vary along the length of the column or as the gross operating parameters and compositions are varied. In oil producing operations the nature and amount of gas phase developed within an underground reservoir as it is being depleted and its pressure declines also involve the basic equilibrium gas-liquid phase interactions. The influence of surface conditions of temperature and pressure and various separation procedures on the nature and amount of stock tank oil and natural gas recovered from a given well stream is likewise determined by the same basic phase equilibrium characteristics of the hydrocarbon systems. Such equilibrium separation of the gas and liquid phases in a well stream is of importance both in establishing optimum separator conditions for obtaining maximum stock tank oil yields and in the general problem of crude stabilization. The phenomena involved in these problems and their quantitative aspects can, of course, be established in each individual instance by appropriate laboratory experimentation. It has been found, however, that by associating with the individual hydrocarbon components functional characteristics de- scribing their individual gas - liquid phase equilibrium behavior, the phase properties of the composite mixtures can be predicted by mathematical computation. These characteristics have been termed "equilibrium constants", or "equilibrium ratios"," and are defined by the general equation: where K, is the equilibrium ratio for the jth component, y, is the rnol fraction of that component in the gas phase, and xj is the corresponding mol fraction in the coexisting liquid phase. Considering the Kj as representing available and measurable properties of the individual components of a mixture, and as defined by Eq. (l), it is readily shown that in a coexisting gas-liquid system the rnol fraction concentrations of the individual components in the gas phase, y,, and those in the liquid phase, x,, are given by the expressions: where the nj are the corresponding rnol fraction concentrations in the composite system, and ng is the mol fraction of the whole which is in the gas phase. The basic unknown in Eq. (2) is n,: which gives the gross separation between the gas and liquid phases. This may be determined in principle by imposing the requirement that the sums of the mol fractions in both the gas and liquid phases are unity, i.e. by the equations: The solution of the equivalent Eqs. (3) lying between 0 and 1, when inserted into Eqs. (2), provides a complete description of the compositions of the coexisting gas and liquid phases. While the Eqs. (3) are essentially equivalent to polynomial equations, such transformations are of little value when dealing with multicomponent sys-

Jan 1, 1949

By Harold L. Overton, Leonard B. Lipson

The first paper in this series1 outlined practical methods for applying the theory of steady-state flow of an ideal Bingham plastic liquid through a circular pipe and axially through a stationary concentric annulus to the enginekring analysis of friction loss problems. In the present paper the properties of another useful rheological model, the pseudoplastic generalized Newtonian liquid, are considered. The terminology applied to the rheological models is derived from the following classification of rheological models. The term "generalized Newtonian" is a class designation applied to models which do not have yield points but which exhibit a dependence on shear rate. The term "pseudoplastic" applies to a sub-group of models within this class which exhibit a decrease in "viscosity" with increasing shear rate. The term "power model" refers to a mathematical model proposed for describing the behavior of the pseudoplastic liquid. This model is of interest since certain of the salt saturated, oil emulsion, and "low solids" drilling fluids, inverted emulsion- and oil-base drilling fluids, aqueous gels, gelled crudes and blocking agents employed in hydraulic fracturing operations are of this type. The present treatment considers the equations describing steady-state Poiseuille flow through a circular pipe and a stationary concentric annul us, and Couette flow between concentric rotating cylinders for a particular pseudoplastic model. As in the previous paper' it is demonstrated how these equations can be applied to engineering calculations of friction losses. Graphical procedures for recovering the various flow curves and cxample calculations are included. A method for simplifying turbulent-flow correlations is also proposed. PROPERTIES OF PSEUDOPLASTIC GENERALIZED NEWTONIAN LIQUID In contrast to the Bingham plastic the generalized Newtonian liquid is yield point-free. However, in the flow region the relation between stress and rate of shear is not linear as it is with ordinary Newtonian fluids. If the stress increases with shear rate at a less than linear rate the system is called pseudoplastic. While there is no yield point in this system, the Bingham advocates' coined the term pseudoplastic to distinguish this kind of behavior from Bingham plastic behavior. For the pseudoplastic a plot of flow rate vs friction loss or rpm vs torque, from pipe flow or rotational viscometer measurements, respectively, generally has three qualities: (1) it can be clearly extrapolated toward the origin; (2) there is very little slope at the origin. This accounts for the high apparent viscosities at low mixing and pumping rates; and (3) the flow curve possesses a continuous curvature which gradually diminishes as the shear increases. Rheological measurements on a pseudoplastic liquid can be easily misinterpreted as indicating Bingham plastic behavior. A typical case in point is illustrated in Fig. 1. These data were obtained from measurements on an oil-emulsion mud using a commercially available rotational-type viscometer (Model 35 Fann V-G meter). Applying the "two-point"* method3 gives a plastic viscosity of 17 cp and a yield point of 24 lb/100 ft2; at the same time weighting material settles very easily in this system. This illustrates that it is not possible to distinguish between the pseudoplastic and the Bingham plastic from only two data points. It is possible to determine qualitatively which model is more appropriate by a combination of dynamic and static measurements. Experience indicates that the simplest criterion for the Bingham plastic liquid is that the calculated dynamic yield point should be of the same magnitude as the measured gel strengths. For the oil-emulsion drilling fluid illustrated in Fig. 1, a value of 4 lb/100 ft2 was obtained for gel strengths measured

By P. van Meurs

This paper describes a technique which permits visual observation of oil displacement processes througlrout the interior of a porous structure as thick as 2 in. A model having glass walls is filled with finely powdered glass. Such a model becomes completely transparent if the powdered glass pack is saturated with an oil having the same refractive index as the glass. When water or gar, with a refractive index different from that of the oil (glass) is injected, the model becomes opaque in the space occupied by the water (gar). visual observation and photography of these displacement processes are thus rendered possible. The following processes have been studied: (I) displacement of oil by water (linear drive) in a homo-geneous formation and in a stratified formation, (2) water flooding (five-spot well pattern), (3) solution gas-drive process, and (4) oil-bank formation. INTRODUCTION Most of the scaled model experiments heretofore carried out in the Koninklijke/Shell Laboratorium, Amsterdam, have been done in steel containers, in which oil-saturated sand packs represented the oil-bearing formation. In such experiments, for instance those on the displacement of oil by water, the cumulative oil production was measured as a function of the cumulative water injection. In order to interpret the measured production data, a need has for some time been felt for a means of observing visually what happens in the course of such experiments. Some information of this kind was obtained earlier by observing the water-oil distribution in various cross sections of the sand after the conclusion of the experiment.' However, the technique used did not guarantee that the fluid distribution remained undisturbed. Moreover, the procedure was somewhat elaborate and the information obtained limited. Therefore a new technique has been developed, which is described below. DESCRIPTION OF EXPERIMENTAL TECHNIQUE AND APPARATUS Glass models simulating part of an oil reservoir are filled with ground Pyrex glass of a specific grain size, representing the porous structure of the prototype. Before an experiment is started, the porous formation has to be saturated completely with an oil having exactly the same refractive index as the Pyrex glass powder. All light-scattering surfaces of the glass grains become invisible and the formation is then transparent. If the oil is displaced by another fluid (e.g. by water in a water-drive experiment, or by gas liberated from the oil in a solution gas-drive experiment) the formation becomes opaque in the region where the oil is expelled from the rock.* The back of the apparatus is painted black and illumination is from the front; therefore, in the photographs of the experiments, the oil-bearing formation appears black, whereas the water- or gas-invaded areas are white.

Jan 1, 1958

By M. R. J. Wyllie, Michael A. Torcaso

With the continued deep drilling of today, increasing numbers of high pressure and high temperature gas-condensate reservoirs are being discovered. Correspondingly, the ranges of properties of gas-condensate reservoir fluids are also being increased. The literature1*'." contains several excellent reviews of the property ranges in presently known condensate reservoirs. When one examines these property ranges, it becomes apparent that the need for reliable data is great if the engineer is to evaluate properly or appraise a gas-condensate reservoir. It is the purpose of this paper to present phase-behavior data that will point out some of the important variables that influence recovery from gas-condensate reservoirs produced by pressure depletion. All these variables may be resolved into two major categories: the hydroc:arbon composition of the reservoir fluid and the physical properties of the reservoir. But the variables which we will consider are those which are usually determined from well tests shortly after discovery and are, consequently, readily available to the engineer. They are as follows: (1) the producing separator gas-stock-tank liquid ratio, (2) stock-tank liquid gravity, and (3) reservoir temperature. Laboratory depletion studies have been performed on several reservoir fluids to illustrate the effect of changes in the cronsidered variables upon ultimate recoveries. The deviation factors, Z, of a series of reservoir fluids were determined experimentally to illustrate their order of magnitude and their effect upon the calculation of in-place reserves existing at high pressurcs and temperatures. The procedures followed in conjunction with the performance of these studies and their associated calculations are illustrated in Fig. 1. For the purpose of this paper, the depletion recoveries are expressed only in terms of pcr cent of initial stock-tank liquid, with the abandonment pressure assumed to be 500 psig. VARIATION OF RECOVERY WITH INITIAL PRODUCING SEPARATOR GAS-STOCK-TANK LIQUID RATIO The producing separator gas-stock-tank liquid ratio has long been one of the principal criteria in assumptions and generalized correlations' used for the estimation of ultimate recovery from gas-condensate reservoirs. Fig. 2 presents a plot of the cumulative recovery of stock-tank liquid at abandonment as a function of the initial producing separator gas-stock-tank liquid ratio. The graphical illustration represents the results of 44 laboratory depletion studies performed on actual gas-condensate fluids. These data are presented not only to illustrate the general trend of ultimate recovery increase with initial producing ratio increase but to emphasize the considerable variation in possible recovery at any given producing ratio. Table 1 presents the properties of two actual reservoir fluids chosen to illustrate the approximate extremities of producing ratios for gas-condensate systems. The fluids designated as A and B exhibited initial producing ratios of 3,400 and 114,000 scf/bbl, respectively. Fig. 3 presents a graphical comparison of the cumulative stock-tank liquid recoveries, calculated from the well-stream compositions and produced volumes for each fluid during depletion. It is to be noted that Fluid A exhibited a recovery of only 16.5 per cent when produced to an assumed abandonment pressure of 500 psig while the corresponding recovery for Fluid B was calculated to be 62.2 per cent. This difference in recovery, which exemplifies the general trend of recovery increase with initial producing ratio increase, is readily explained by consid-

By J. S. Aronofsky, J. P. Heller

This paper presents a mathematical analysis of the fluid mixing which occurs during flow through porous media. The analysis is based on the well-known diffusion equation with mass transfer term. It is pointed out that the use of this equation is justified by a single general assumption which does not specify any particular mechanism for the mixing process. Formulas are given for tracer fluid concentrations for two different boundary and initial conditions. Calculated numerical values compare closely with some published results of experiments in which no viscosity or density difference existed between displaced and displacing fluids. INTRODUCTION A phenomenological theory of the mixing and diffusion process is presented for the flow in a porous medium of two miscible liquids of equal viscosity and density. This is based on the classical diffusion equation which has been used to explain such processes as Brownian motion and heat conduction. More recently, this same differential equation has been derived (in a manner analogous to the work in this paper) for the problem of heat transfer during fluid flow in porous media.' The purpose of this paper is to show that the available experimental evidence justifies the use of this simple diffusion equation. The recent data of Koch and Slobod' and of von Rosenberg' are used to test this diffusion model by comparing concentration curve shapes and calculating "effective diffusion coefficients". The same differential equation for miscible displacement, as well as some experimental confirmation, has been published in Japan by K. Yuhara,4 and recently in this country by Day.5 Yuhara arrives at his con-clusions through analogy with turbulent diffusion, de- scribing the microscopic velocities in terms of "eddy motion" — although unlike the ordinary eddies of hydraulics. Day proceeds directly from Scheidegger's6 concept of the dispersion of velocities. After a very clear description of Scheidegger's theory, he applies it to calculate the motion first of a drop and then of a thin layer of brine into a bed of sand. For the latter case, Day derives and solves the same differential equation as Eq. 1. In our presentation no detailed assumptions are made as to the mechanism of displacement. In fact, several microscopic processes may well be involved—the only general assumption made is that their cumulative microscopic effect may be described by the addition of a second-order term to the differentirl equation of continuity. DERIVATION OF EQUATIONS Let us consider the displacement in one spatial dimension of Liquid A from a sample of porous media by Liquid B. A and B are assumed completely miscible and of the same viscosity and density. Both fluids are also assumed to be incompressible and it is specified that the total flow is the same for any cross section. It is evident that the rate at which Liquid B would be carried across such a section would be given by the product of the fluid velocity, u, and the concentration, C, of B in the fluid mixture assuming the constancy of u and C over the cross section. Since the volume of Fluid B is conserved, this leads to the one-way wave equation - u ?C/?x = ? ?C/?l where x increases in the direction of u—the "Darcy velocity", Q/A (or velocity flux per unit area), and 9 is the fractional porosity. Consider the case where fluid B is displacing Fluid A and the concentration distribution is represented by a step funcLon (with a vertical discontinuity). The solution of Eq. 1 predicts no mixing

Jan 1, 1958

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

A generalized mathematical model is presented, which describes the thermal recovery of oil in linear systems with convective external heat loss. Three-phase fluid flow, conduction-convection heat transfer, chemical reaction between oxygen and oil and aqueous phase change are included, making the model applicable to a variety of thermal recovery processes. The complexity of the model requires that numerical solutions be obtained with a digital computer. The numerical approximations of the partial differential equations are presented and the algorithm is discussed in detail. The analysis is applied to the detailed simulation of a forward-combustion laboratory tube experiment. The results of the simulation exhibit with reasonable accuracy all of the principal thermal and bydrodynamic characteristics observed in the laboratory (i.e., the propagation of a combustion zone, a steam plateau, and oil and water banks). Numerical variations in the parameters appearing in the external heat-loss term, the aqueous phase-change term, and the Arrhenius combustion term are examined, and their effects on the predicted temperature profiles, oil and water production histories, and the air-oil ratio are presented. In spite of an oversimplification in the chemical reaction mechanism and the considerable computer time required to obtain solutions, the model represents a major advance in the ability to simulate the phenomena observed in thermal recovery experiments. INTRODUCTION Thermal oil recovery refers to a class of recovery processes where heat is supplied to a reservoir to provide the necessary expulsive energy. This thermal energy can be supplied externally as steam or hot water, or it can be generated in situ by forward or reverse combustion. In either case, however, thermal recovery processes are charac- terized by the simultaneous flow of two or three fluid phases in a variable - temperature field, accompanied by possible chemical reaction or phase-change effects.2-5, 7- 10, 14,16-18,20-22 Although a physical understanding of the thermal recovery processes is far from complete, it is possible to construct mathematical models which describe approximately all of the principal physical and chemical phenomena. However, attempts to solve such models, even with high-speed computers, involve formidable mathematical difficulties. Consequently, theoretical solutions have been obtained only for idealized cases in which important physical phenomena are neglected. For example, consider the process of forward in situ combustion. All such theories which have been developed consider only certain aspects of the process, such as heat transfer, 2-4,8,9,l6,l8 heat transfer with phase change,10 heat transfer with chemical reaction,7 or the hydrodynamics of three-phase flow.21 A general theory including all of the above phenomena has not been developed to date. This paper presents a unified theory of thermal recovery processes in linear systems. A mathematical model is developed which explicitly includes conduction-convection heat transfer with convective external heat loss, chemical reaction between air and oil, aqueous phase change, and the hydrodynamics of three-phase flow. A system of equations is developed which can be solved numerically on a high-speed digital computer, resulting in predicted temperature, pressure, and saturation histories in space and time. The model allows a more detailed simulation of thermal recovery tube experiments than had previously been possible. THEORETICAL DEVELOPMENT Consider the linear flow of gas, water and oil in a homogeneous porous medium. Assume that the oil will react with gaseous oxygen, and that mass is transferred between the water and gas phase by evaporation or condensation. In the absence of hydrocarbon phase change, intraphase diffusion,

Jan 1, 1966

By Kenneth F. Whinery, John M. Campbell

Basically, the Buckley-Leverett theory involves two systems which are similar in nature but are differentiated by time. These systems may be described by the fractional flow and frontal advance equations which essentially characterize the mechanics of oil movement while being expelled from the reservoir. The fractional flow equation originally developed by Leverettl may be expressed in the more usable form The development of this equation is based on Darcy's law describing fluid flow through porous media and applies to the flow at only one point (it is a point function). For simplification, the fractional flow equation is written in the above form because the capillary forces always increase the fractional flow of the displacing phase regardless of the direction of flow or the displacing phase. If, for simplicity, the effects due to gravity and capillary pressure differences are neglected, the fraction of the displacing fluid, f, at any point in the flowing stream is related to the properties of the system by 1 Thus, it is seen that in the absence of capillary and gravitational effects, fd for a given sand and fluids varies only with saturation and pressure. The magnitude of the viscosity ratio, —a has an effective range (range of about 30) in the system where gas is displacing oil and a much smaller range for water displacing oil. In order to make the fractional flow equation more versatile, it is necessary to connect the fractional flow at a given point and saturation with time. This problem was approached by Buckley and Leverett" who developed the frontal advance equation, Eq. 3 states that the rate of advance of a plane that has a fixed saturation, Sd, is proportional to the change in composition of the flow stream caused by a small change in the saturation of the displacing fluid. It is, essentially, a transformation of a material balance equation representing the net rate of accumulation of the displacing phase within a homogeneous sand block. This accumulation is proportional to the difference between the rate at which the displacing fluid enters the sand and that at which it leaves. Eq. 3 describes the velocity with which a plane of constant displacing phase saturation advances through a porous system. Buckley and Leverett," Babson,' Kern," Welge, and others have adequately discussed the basic mechanism and application of the fractional flow and frontal advance equations. APPLICATIONS OF THE FRONTAL ADVANCE THEORY TO PETROLEUM RECOVERY Two general applications of the Buckley-Leverett frontal advance theory involve the system in which the oil is being displaced by an expanding gas cap overlying the oil zone and that in which the oil is being displaced by water. Displacement by Gas in the Presence of an Immobile Water Saturation A system in which gas is the displacing phase may be thought of as having two forces effecting the displacement process. These forces are the gravitational force and that force exerted by the displacing gas. The gravitational effects control the displacing efficiency of the gas. The gravitational effect will be less at higher rates of flow, thereby reducing the effectiveness of the displacement of the oil by the gas. The more efficient displacements occur at flow rates which are less than the gravity free fall rate. Capillary forces can be neglected without materially changing the magnitude of the gas saturation. The Mile Six pool is used herein to illustrate the calculating procedures in evaluating gas drive-gravity drainage field perfomlance. These calculations represent the determination of the gas-oil contact when the distribution of the hydrocarbon pore volume is considered. Two methods

By R. Ehrlich, F. E. Crane

A consolidated porous medium is mathematically modeled by networks of irregularly shaped, interconnected pore channels. Mechanisms are described that form residual saturations during immiscible displacement both by entire pore channels being bypassed and by fluids being isolated by the movement of an interface within individual pore channels. This latter mechanism is shown to depend strongly on pore channel irregularity. Together, these mechanisms provide an explanation for the drainage-imbibition-hysteresis effect. The calculation of steady-state relative permeabilities, based on a pore-size distribution obtained from a Berea sandstone, is described. These relative permeability curves agree qualitatively with curves that are generally accepted to be typical for highly consolidated materials. In situations where interfacial effects predominate over viscous and gravitational effects, the following conclusions are reached. 1. Relative permeability at a given saturation is everywhere independent of flow rate. 2. Relative permeability is independent of viscosity ratio everywhere except at very low values of wetting-phase relative permeability. 3. Irreducible wetting-phase saturation following steady-state drainage decreases with increasing ratio of nonwetting- to wetting-phase viscosity. 4. lrreducible wetting-phase saturation following unsteady-state drainage is lower than fir steady-state drainage. 5. Irreducible nonwetting-phase saturation 101lowing imbibition is independent of viscosity ratio, whether or not the imbibition is carried out under steady- or unsteady-state conditions. Experiments qualithtiveIy verify the conclusions regarding unsteady-state residual wetting-phase saturation. Implications of these conclusions are discussed. INTRODUCTION Natural and artificial porous materials are generally composed of matrix substance brought together in a more or less random manner. This leads to the creation of a network of interconnected pore spaces of highly irregular shape. Since the geometry of such a network is impossible to describe, we can never obtain a complete description of its flow behavior. We can, however, abstract those properties of the porous medium pertinent to the type of flow under consideration, and thus obtain an adequate description of that flow. Thus, the Kozeny-Carmen equation, by considering a porous medium as a bundle of noninterconnecting capillary tubes, provides an adequate description of single-phase flow. With the addition of a saturation-dependent tortuosity parameter in two-phase flow to account for flow path elongation, the Kozeny-Carmen equation has been used to predict relative permeabilities for the displacement of a wetting liquid by a gas. 1-5 It has long been recognized that relative permeability depends not only on saturation but also on saturation history as we11.6,7 Naar and Henderson8 described a mathematical model in which differences between drainage and imbibition behavior are explained in terms of a bypassing mechanism by which oil is trapped during imbibition. Fatt9 proposed a model for a porous medium that consisted of regular networks of cylindrical tubes of randomly distributed radii. From this he calculated the drainage relative permeability curves. Moore and slobod,10 Rose and Witherspoon,11 and Rose and Cleary12 each considered flow in a pore doublet (a parallel arrangement of a small and large diameter cylindrical capillary tube). They concluded that, because of the different rates of flow in each tube, trapping would occur in one of the tubes; the extent of which would depend upon viscosity ratio and capillary pressure. They limited

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 A. Schild

The effect of heating a producing well on the rate of oil recovery has been analyzed in the simple case of a well producing oil by a radial drive and in the steady state. Differences of production rates with and without borehole heating have been calculated and plotted; these measure some of the benefits of the thermal recovery process. A striking feature is that the production increase tends rapidly to a constant value for large production rates. Some concrete examples, including a rough comparison with field data. are discussed. Heating requirements have been calculated and plotted on the assumption that the heaters used are 100 per cent efficient. Heating has been thought of as applying only to low volume wells producing heavy oils. However, the calculations suggest that the thermal process considered might be profitable even for some wells producing large quantities of light oil, especially in unconsolidated formations where fracturing is not possible. In the case of wells producing other fluids with the oil or of wells producing partly or wholly by gravity drainage, the theory developed does not apply in detail. It can be expected, however, that any general qualitative be-havior will remain valid. Under normal temperature conditions, high viscosity emulsions sometimes form around the wellbore, lowering the permeability in the vicinity of the well and impeding the production of oil. Borehole heating often prevents such emulsification. This important effect is not easily amenable to mathematical analysis and has not been included in the theory presented here. It should therefore be kept in mind that the benefits due to borehole heating may in fact exceed the values given by the following calculations. INTRODUCTION The recovery of viscous, low gravity oils can often be stimulated by applying heat to the production zone. The higher temperature lowers the viscosity of the oil, increasing its mobility, and therefore its rate of production. One thermal method of recovery 1,2,3,4. consists in operating a heater inside a producing well and, at the same time, producing oil from the well. In this paper a simple theoretical model is analyzed to study the effect of such heating on the recovery of oil. The reservoir mechanics of the process is the only aspect of the problem which is considered here. The details of the heating system are ignored, and it is simply assumed that the heater raises the temperature in the immediate vicinity of the producing well to a pre-assigned temperature T w. Although heating requirements are calculated, these refer only to the heat supplied to the formation and to the oil being produced; they do not include heat losses due to the inefficiency of the heater. A short outline of the theory follows: The heat supplied to the producing well flows out into the formation by means of thermal conduction. At the same time, the flow of oil carries heat by convection from the formation towards the producing well. After sufficient time has elapsed, a steady-state is reached where the temperature, the velocity of the oil, etc., no longer vary with time but are functions of position only. In the steady state the temperature decays like an inverse power from the high value Tw at the wellbore to the original reservoir temperature T. at some distance away from the well. The effect of the high temperature and low oil viscosity near the producing well is the same as if the oil were produced at reservoir temperature To but through a producing well of a larger radius. From this point of view, the energy and money being expended to provide heat

Jan 1, 1958