Search Documents

By R. I. Parsons, M. B. Standing, E. N. Lindblad

The recovery of the heavier components from a gas cap or retrograde pool is shown to be the greatest when the sand is cycled with a dry gas at a low pressure. This conclusion is in direct opposition to the belief that the most efficient production program is pressure maintenance and cycling at or near the dew point. The results are calculated from: (I) constant volume, variable composition pressure-volume-temperature tests on a mixture of trap gas and liquid from a producing well; (2) published equilibrium constant data and the measured composition of the dew-point material; and (3) the fact that in a sand section of homogeneous permeability, injected gas displaces reservoir gas nearly, quantitatively. The results are based on the simplifying assumption that variable permeability systems may be defined by the ratio of two statistical parameters, and that gas injected into an actual sand will behave as though the sand were composed of many layers, each of constant permeability. Introduction Pressure decline in gas-condensate type reservoirs is accompanied by the formation of a liquid phase throughout the reservoir. Over the past ten years the processing of the material from these types of reservoirs for the heavier hydrocarbon components and the return of the light fractions to the reservoir ("cycling") has become increasingly popular. It has often been stated that the purpose of such a program is to prevent the loss of the retrograde liquid phase formed in the reservoir. The purpose of this paper is to present the results of laboratory tests and computations concerned with several possible methods of producing a gas cap or condensate type of reservoir. The results show that the recovery of heavier hydrocarbons for this type of reservoir is not a maximum under conditions of cycling at or near the dew-point pressure. Instead, variations in permeability and the ability of the dry injected gas to revaporize condensate, point to cycling at a considerably reduced pressure as the optimum production-method. In addition the paper suggests a way of evaluating sands for their permeability variation. The calculations presented in the following sections of the paper are in terms of the production of the butanes and heavier fraction from the gas cap of a field which had an original pressure of approximately 3000 psi. Specifically, the calculations show: 1. To what extent the condensate formed in the reservoir after pressure decline can be recovered by the cycling of a dry absorber-plant gas through the sands. 2. Whether there would be greater recovery of hydrocarbons if cycling were instigated at the original reservoir pressure.

Jan 1, 1948

By M. B. Standing, R. I. Parsons, E. N. Lindblad

The recovery of the heavier components from a gas cap or retrograde pool is shown to be the greatest when the sand is cycled with a dry gas at a low pressure. This conclusion is in direct opposition to the belief that the most efficient production program is pressure maintenance and cycling at or near the dew point. The results are calculated from: (I) constant volume, variable composition pressure-volume-temperature tests on a mixture of trap gas and liquid from a producing well; (2) published equilibrium constant data and the measured composition of the dew-point material; and (3) the fact that in a sand section of homogeneous permeability, injected gas displaces reservoir gas nearly, quantitatively. The results are based on the simplifying assumption that variable permeability systems may be defined by the ratio of two statistical parameters, and that gas injected into an actual sand will behave as though the sand were composed of many layers, each of constant permeability. Introduction Pressure decline in gas-condensate type reservoirs is accompanied by the formation of a liquid phase throughout the reservoir. Over the past ten years the processing of the material from these types of reservoirs for the heavier hydrocarbon components and the return of the light fractions to the reservoir ("cycling") has become increasingly popular. It has often been stated that the purpose of such a program is to prevent the loss of the retrograde liquid phase formed in the reservoir. The purpose of this paper is to present the results of laboratory tests and computations concerned with several possible methods of producing a gas cap or condensate type of reservoir. The results show that the recovery of heavier hydrocarbons for this type of reservoir is not a maximum under conditions of cycling at or near the dew-point pressure. Instead, variations in permeability and the ability of the dry injected gas to revaporize condensate, point to cycling at a considerably reduced pressure as the optimum production-method. In addition the paper suggests a way of evaluating sands for their permeability variation. The calculations presented in the following sections of the paper are in terms of the production of the butanes and heavier fraction from the gas cap of a field which had an original pressure of approximately 3000 psi. Specifically, the calculations show: 1. To what extent the condensate formed in the reservoir after pressure decline can be recovered by the cycling of a dry absorber-plant gas through the sands. 2. Whether there would be greater recovery of hydrocarbons if cycling were instigated at the original reservoir pressure.

Jan 1, 1948

By C. S. Matthews, G. L. Stegemeier

In one field in South Texas, approximately 72 per cent of the pressure build-up results show a characteris-i.rtic "hump" (i.e., the pressure builds up and then falls off) which makes interpretation by. standard methods impossible Correlation of size and time of hump with formation permeability. well productivitv index, and method of completion led to the tentativc, conclusion that the humps were caused by segregation of gas and oil in the wellbore after closing-in"'. This conclusion was confirmed by performance of simple laboratory bubble-rise experiments, by theoretical bubble-rise time calculations, and by a detailed calculation of PVT behavior in the well of a particular well Oil which accurate Surface and bottom-hole pressure measure-ments were untile. The hump behavior has since been found to occur in many other fields. The cause, however, is not the mine in all cases. In some of these the hump is trace-ahie to leaks in the tnbing which allow influx from the annulus after closing-in. In other cases the hump is traceable to leaks in the device separating pay horizons in dually completed wells. It is concluded that the recording of both surface and bottom-hole pressures is desirable in wells which show an anomalous build-up behavior. A number of field examples is discussed where use of both sets of measurements enables the cause for anomalous behavior to he found, (and a reasonable interpreation of bottom-hole pressure to be made. INTRODUCTION Theoretically the pressure build-up in an infinite reser-voir should be a linear function of In [(t + at)/nere t is the production time and At a the closed-in time. Some of the variations from this behavior arc well known, such as the curved portion immediately after shut-in which results from after-production and skin effect, and the flattened end portion which results from boundary effects in a limited reservoir. The effect of stratified producing zones and irregular geometrical drainage patterns may also contribute unusual characteristics to build-ups. However, the effect of still another phenomenon—that of movement of fluids within the wellbore—has been neglected in most build-up studies to date. Fluid movements within the wellbore after shut-in occur as a result of after-production, packer failures, leaks in the casing or tubing, or buoyancy of the gas phase when both gas and liquid are present in the well. Each of these movements can influence the pressure build-up, sometimes sufficiently to negate use of the data for computing permeability or stabilized build-up pressure. Examples of each of these phenomena are discussed below. FIELD OBSERVATIONS ON HUMP BUILD-UPS In approximately 75 per cent of the wells in a medium-sized field in South Texas the pressure build-up curves rise to 3 maximum and then decline to what appears to be a stabilized reservoir pressure as shown in Fig. I. Similar behavior has also been observed in wells in other fields in Texas and Louisian:,. Because of the possibility of such behavior being a mechanical failure of the pressure bomb. the bombs were run in tanden into one well which had been known to exhibit this behavior. When both bombs recorded identical curves with the characteristic hump, it was inferred that the unusual build-up is a result of well or forma-

By G. W. Lester, T. J. Nowak

The paper presents a practical method of utilizing pressure fall-off data obtained when a water injection well is shut in for determination of: (I) the static reservoir pressure, (2) the potential water injective capacity (md-ft) of the formation, and (3) the formation damage (skin effect) surrounding the wellbore. In the development of the method, equations are presented which give The theoretical background for predicting the behavior of water injection wells. Certain physical factors which limit the application of the method are discussed. The procedure and the relevant precautions in obtaining the data are discussed. A sample analysis and computation from an actual pressure fall-off curve is illustrated. Fintrlly, general characteristics of pressure fall-off curve! obtained in several reservoirs are presented. INTRODUCTIO N Theoretically, in an injection well the flow of water into the formation will continue as long as the sand-face pressure is greater than the reservoir pressure. However, in practice, engineers have observed that upon stopping injection and opening the wellhead to atmospheric pressure, some wells backflow water for time intervals ranging from several seconds to several hours although the sand-face pressure, as calculated from the hydrostatic head, is considerably greater than the static reservoir pressure. Further, field observations indicated that upon shutting in a water input well, the wellhead pressure does not drop off immediately, but a pressure decline occurs lasting from several seconds to several days. These observations make it evident that when an input well is shut in, a residual pressure lingers in the vicinity of the wellbore which is greater than the average reservoir pressure. This residual pressure declines first rapidly, then more slowly, to the static reservoir pressure: therefore, it is called transient back-pressure or pressure fall-off. The first step in investigating the mechanism of pressure fall-off is to establish what general conditions are necessary for its occurrence, and to this end, Dickey and Andresen1 have presented a qualitative explanation describing the mechanism. They point out that the formation pressure gradient in the vicinity of the well does not disappear immediately when the input well is shut in, as if the system were filled with a totally incompressible fluid, but it is maintained by the expansion of the compressible fluids and by the elasticity of the formation. In such a compressible system the decrease in pressure is possible only to the extent that the water front continues to advance until the pressure is equalized. The principal aim of this paper is three-fold; namely. (1) to present a theoretical quantitative analysis of pressure fall-off data, (2) to present an application of the developed theory to field data, and (3) to discuss the general characteristics of actual pressure fall-off curves. Under (I), theoretical equations dealing with such factors as compressibility of the fluid-reservoir system, mobility of the fluids and certain physical restrictions that

Jan 1, 1956

By F. W. Jessen, G. F. Sears

Cavities in massive salt for the purpose of storage of liquid hydrocarbons have assumed a prominent position in recent years. This paper describes a program to facilitate leaching operations for the formation of specifically shaped storage cavities. Various forms and sizes of cavities may be possible through use of the techniques developed. INTRODUCTION The creation of large under ground storage facilities for natural gas and liquified petroleum gases has been practiced for many years. Use of cavities obtained through solution of salt for this purpose is a fairly new and novel approach, but has gained increasing importance due mainly to the economics of this type of mining operation opposed to hard rock mining or surface and pit storage.1 The idea of controlling the shape of any cavity dissolved from massive salt has not been a prime consideration of companies engaged in the formation of storage space. This has been due mainly to an insufficient knowledge of the mechanics of the leaching process and a dearth of published information dealing with both the desirability and ease of control possible for this type of operation. Two distinct advantages are readily ascertainable. From a stability standpoint (i.e., the ability of the completed cavity to withstand stress imposed by overburden pressure and lateral tectonic stresses) a controlled cavity may be generated which will yield the most favorable attitude to these external forces and remain in operable use for longer periods of time.2 A second advantage is that for a given volume, a sphere (which is one controlIed shape possible) represents the minimum surface area exposed. This may become more important in the future when dealing with refrigeration and product losses in underground cavities. The basic solution mining process, without regard for controlling the final shape, is quite simple in that the equipment and materials required to dissolve a cavity in a salt dome or layered salt section are a source of fresh water, a circulating pump, several strings of tubing and a means of disposal of the return brine. To add control measures involves the use of an inert blanket material above the position at which solution proceeds. Initially the annulus of the largest wash pipe is filled with this blanket down to the top of the proposed cavity. The bottom of the proposed cavity is determined by the depth of the original drilled hole or a blanking plug. The blanket material is added incrementally in stages and displaces the water in the enlarging cavity downward. Where the blanket material has displaced the water, no further solution takes place. The rate at which to add this controlling blanket material so as to be able to form specifically shaped cavities of any particular size is of considerable importance. THE PROBLEM It is desirable to have as much information as possible as to the behavior of the mining operation since visual observation obviously is impossible. It would be advantageous, for example, to have a step-by-step program showing how much salt is to be removed, at what rate and the time required. This type of program would facilitate procurement of surface equipment and provide for proper management of manpower requirements, as well as yield a host of intangible benefits. As mentioned earlier, only rule-of-thumb estimates were available until recently and made this type of pre-planned operation haphazardous at best. Recent work3-' has done much to clear up the ambiguities and inaccuracies attendant to solution mining of salt and has helped to place the operation on a more scientific basis. This study attempts to correlate much of the information thus far developed into a complete mining program; to extend the idea of controlled solution mining to include not only spheres, but solid conic sections of the ellipsoidal variety; and to refine the computation of the rate of removal of salt during the various stages of mining. Dommers3 and IIemson4 made extensive laboratory

Jan 1, 1967

By B. K. Larkin, H. R. Bailey

One of the most complicated and potentially one of the most promising secondary recovery methods is that of underground combustion. A number of field tests1,2,3 have been performed recently, apparently indicating that the underground combustion process is technically feasible. Relatively little has been published concerning laboratory or theoretical investigations of the underground combustion process. A few results of laboratory tube runs have been published4,5 and these data are valuable in understanding the process. The tube run data also emphasize the complexity of the thermal recovery process. Under various conditions, problems in heat transfer, multiphase fluid flow and combustion rates arise. For example, at low temperatures the combustion process proceeds so slowly that the front cannot be propagated. It is a problem in heat transfer to predict what injection rate of air is necessary to sustain combustion in a particular reservoir. Certain fluid flow aspects of underground combustion have been studied by Wilson, et a13. In the development of some of the flow equations, assumptions of adiabatic combustion were made. Vogel and Krueger1 have devised an analog model of the heat conduction process for a constant temperature cylindrical source moving radially outward in a homogeneous conducting medium. Jenkins and Ramey7, in a discussion of the Vogel and Krueger paper, present an analytical solution of a heat conduction problem related to problems considered by Vogel and Krueger. In the present paper an idealized model of the heat conduction problem is considered, which is more general than the model of Vogel and Krueger or Jenkins and Ramey. In particular, initial well heating, vertical heat losses and arbitrary frontal velocities are included. Equations for the temperature are obtained in terms of integrals. These integrals have been evaluated numerically for the case of frontal movement corresponding to a constant air injection rate and the case of constant frontal velocity. These results are presented graphically. From these graphs one may estimate conditions under which nearly adiabatic temperatures are obtained. The application of these results requires some knowledge of reservoir thermal properties, combustion rates and residual fuel content. Assuming no vertical losses it is possible to obtain explicit evaluations of these integrals in the linear con- stant velocity case and in the radial case corresponding to a constant air injection rate. PRESENTATION OF RESULTS The geometry of the radial system is illustrated in Fig. 1 where the reservoir is represented as a slab of thickness h. By making an energy balance for an infinitesimal volume, one obtains the following equation. T= x T+F(z.r).........(1) where F (z,r) =HWv1/pc (r-r1), if z h/2 = 0, if z >h/2 . In Eq. 1, V- represents the Laplacian operator and 8 is the Dirac delta function. Eq. 1 is the well-known diff-sion equation with a source term and may be applied to a linear system. For example, if the source is moving in the direction of the x axis, then the argument of the delta function is x - 1,. In deriving Eq. 1 the following approximations were made concerning the actual combustion process: (1) all heat transfer is by conduction in a homogeneous medium, (2) the formations bounding the reservoir have the same thermal properties as the reservoir, (3) the combustion fuel is consumed instantaneously as the front arrives, (4) there is a constant fuel concentration

By R. E. Cook

This paper presents results of a study to determine to what extent errors in estimated free gas saturation affect the results of static pressure calculations from build-up curves in two-phase systems. Use is made of the method of pressure build-up curve analysis developed by Miller, Dyes and Hutchinson1 , with modifications for two-phase flow advanced by Perrine2. Cornhined fluid compressibility of gas-oil mixtures is included, with special attention given to the discontinuity that occurs in fluid compressibility at the saturation pressure and its effect on static pressure calculations from build-up curves. The sensitivity of static pressure calculations to errors in gas saturation is considered over wide ranges of formation and fluid properties and reservoir pressure conditions. It is concluded from this study that for analyses of two-phase pressure build-up curves (I) the use of under-saturated crude compressibility data, even when producing GOR's are at or near solution values, will result in very significant negative errors in calculated static pressrrre, and (2) except for large values of buildup curve .slope and/or very low crude gravities, rather large errors in calculated gas saturation can he tolerated without .significant effect on calculated static pressures. The conc1usions arc qualified as referring only to the effects of gas saturation and excludc possible errors in the nuumerous other factors that enter into the ancllysis of pressure build-up curves. INTRODUCTION As part of the analysis of basic data for several major Eocene reservoirs in the Bolivar Coastal field (Western Venezuela), it was necessary to analyze a large number of pressure build-up curves for oil permeability and static reservoir pressure. From performance data it became apparent that gravity segregation is very active in these reservoirs. As a result of this performance it was recognized that the free gas saturation within the oil leg of the reservoirs would be small but would be subject to certain errors in calculation due to the inaccuracy of relative permeability data in the low gas saturation region. Since all of the various methods of build-up curve analysis were originally derived for sin- gle-phase, constant compressibility systems, and since the presence of a small but rather indeterminate gas saturation would greatly affect the reservoir fluid compressibility, a study was undertaken to investigate the effect of gas saturation on the results of two-phase pressure build-up analyses. TWO-PHASE BUILD-UP CURVE ANALYSIS The method of analysis used in this investigation is basically that originally developed by Miller, Dyes and Hutchinson1, with modifications for two-phase flow advanced by Perrine2. Miller, et al, showed that the original equation for single-phase permeability away from the well bore, held for two-phase flow except that the calculated permeability is the effective permeability to oil away from the wellbore. This finding was acknowledged by Per-rine, who also presented the equation for effective gas permeability* as Since Eqs. 1 and 2 contain no terms influenced by fluid compressibility, no distinct problem should develop in calculating permeability from two-phase build-up curves. In regard to the calculation of static reservoir pressure from build-up curves in single-phase systems, the modified Miller, et al, method as presented by Perrine is briefly as follows. 1. From build-up curve data, PVT data for the reservoir fluid and the well spacing and sand thickness, calculate the climensionless shut-in time corresponding to the highest value of actual shut-in time (t, hrs) on the straight-line portion of a semi-log build-up curve. Dimenrsionless shut-in time is related to actual shut-in time by the relation, t- 0.000264 kt 2. From a graph of dimension less shut-in time vs dimensionless pressure difference (Fig. 1 1 ) find the value of dimensionless pressure difference corresponding to the value of 7 calculated under (1) above. Dimen-

By W. N. Hiatt

The mathematical theory of the flow of two-phase, incompressible fluid through porous media is clarified, and the development of a general fluid displacement equation for two-phase, incompressible vertical flow through porous media is outlined. The relationship between the three components of the equation — the "fluid drive effect", the "gravity effect", and the "capillary effect'' — is discussed, as well as implications of the analysis in connection with the study of certain types of petroleum reservoirs. INTRODUCTION The purpose of this paper is to clarify and summarize the mathematical basis of the theory of two-phase incompressible fluid flow through porous media, particularly as it pertains to vertical flow, and to present a qualitative* exposition of the controlling factors involved in circumstances where this type of flow may exist. Although the generally accepted theory of heterogeneous fluid flow through porous media was developed a number of years ago by Muskat,1 and Buckley and Leverett2, 3 using the basic premise of Darcy, (and since then expounded upon by many others) there is a need for an over-all summary and clarification of this theory, especially regarding certain types of applications in connection with so-called "gravity-drainage-type reservoirs". The source and derivation is outlined for what is referred to here as the "General Fluid Displacement Equation for Incompressible Vertical Two-Phase Flow through Porous Media''. This equation would apply under the condition where the pressure drop along the path of flow is negligible in comparison with the pressure level on the system, or where compressibility and solubility effects are negligible. These conditions exist ordinarily in gravity-drainage-type reservoirs. Specifically, it will be shown that this general equation may be developed directly from the basic equations, which, according to present theory, govern the general case of viscous heterogeneous flow through porous media. This general equation has appeared in the literature previously, but in most cases its development from basic theory is not indicated clearly, or the equation is presented in terms of groups of constants and variables involving reduced or dimensionless quantities. This presentation is expedient in cases of experimental or computer applications; but it may make it difficult for the practicing petroleum engineer to identify all the factors involved or to understand the relationships between them. Therefore, in the following development, only ordinary mathematical parlance will be used, and no terms involving groups of other equations or reduced or dimensionless quantities will be used. DERIVATION The following equations for two-phase oil and gas flow result from analogy with Muskat's general equations,** if solubility and compressibility effects are neglected, and if all volumes are considered in terms of reservoir conditions, and the water phase, if present, is considered to be immobile (or effectively a part of the porous medium.) In this case, the sum of the oil and gas saturations may be taken as unity, and the saturation of each phase is in terms of the total hydrocarbon saturation. Note that in Eqs. 1 and 2, and in the derivation following, z is taken as positive in the vertically downward direction,

Jan 1, 1969

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By W. A. Klikoff, I. Fatt

The concept of fractional wet wattability is examined. Fractional water wettability of a reservoir rock is defined as the fraction of the internal surface urea that is in contact with water. Capillary pressure and relalive permeability of unconsolidated sand are shown to be functions of fractional wettability. INTRODUCTION The petroleum industry has long recognized that wettability of reservoir rock has an important effect on multiphase flow of oil, water and gas through reservoirs. As early as 1928 the American Petroleum Institute sponsored a study of wettability as part of API Project 27 at the U. of Michigan.' Despite 30 years of research, there is still little exact knowledge of the wettability of reservoir rocks. There are two parts to the wettability problem. After agreeing to a uniform nomenclature in regard to wettability,' the first question to be answered is, "What is the in situ wettability of a given reservoir rock?" If this can he answered the next question is, "What part does wettability play in determining the characteristics of multiphase fluid flow through the rock?" This paper represents an oblique attack on the problem of wettability. No attempt is made here to answer the basic question of wettability in situ. instead the consequences of the concept of fractional wettability are examined. Multiphase flow in sandpacks is shown to he highly influenced by fractional wettability. Jennings' has given 3 definition of wettability and the other terms used in discussing wettability. These terms must be applied to the physical situation existing In reservoir rock. A survey of the pertinent literature from 1928 to 1956 indicates that the concept of a contact angle was applied to reservoir rock in the same way it would be applied to a flat, homogeneous surface. Attempts were made to state quantitatively the wettability of a reservoir rock in terms of a contact angle which was presumably constant at all points on the very rough and heterogeneous interior surface of a porous rock. Calhoun, et al,3-6 prepared synthetic consolidated and un-consolidated porous media in which they claimed there was a known uniform contact angle. They then showed the effects on the capillary pressure and relative permeability characteristics of varying this angle. The API Project 47 at the U. of Texas' and others' have made extensive studies of an indirect approach to the contact angle through the use of heat of wetting data. Even if successful. however, this approach also states the wettability of porous rock in terms of a contact angle which is uniform over the entirc surface. If the angle varies from one part to another on the internal surface, there is no way of determining from the measurements the area distribution of contact angles. In 1956 Brown and Fatt5 suggested that the concept of a contact angle, as applied to reservoir rock, be abandoned. This suggestion was made because it is known that the internal surface of most reservoir rocks is composed of many different minerals, cach with a different surface chemistry and a different capacity to adsorh surface active materials from reservoir fluids. Furthermore, the operation of a contact angle in determining the form of a fluid-fluid interface is difficult to picture in the very complex geometry of a pore. Brown and Fatt proposed that the wettability of reservoir rock be stated in terms of the fractional internal surface area that is in contact with water or oil. All surfaces on which there is water are called water-wet; surfaces on which there is oil are called oil-wet. The fractional water wettability is then stated as a number which represents the fraction of the internal surface that is in contact with water. A symmetrical statement can be made for the fractional oil wettability. The concept of a fractional wettability as previously stated has in its favor the recognition of the heterogeneous mineral composition of most reservoir rocks. Another point in its favor is that fractional wettability can be measured quantitatively with relative ease. Hol-brook and Bernard10 se a simple dye adsorption test to obtain fractional wettability of reservoir rocks. Amott11 uses a combination of imbibition and displacement to arrive at a wettability index of reservoir rocks which seems related to fractional wettability in the range 0.25 to 0.75 fractional water wettability. Jennings" has shown the changes in relative permeability that take place when a porous material is changed from unity to zero fractional water wettability. He also shows that reservoir rocks in their natural state, but at room temperature and atmospheric pressure, have relative permeability characteristics which would indi-

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

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

Jan 1, 1969

By C. Perez-Rosales

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

Jan 1, 1970

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

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

Jan 1, 1970

By J. E. Berry

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

By E. C. Barfield, D. G. McCarty

Two-dimensional analyses offer considerable promise in providing the basic information required to effect more precise control of petroleum reservoir performance. This paper describes a method for conducting some engineering analyses of this type using a high-qreed digital computer. The general approach is to (I) develop a representative numerical model of the reservoir and (2) employ a suitable numerical technique to solve the basic equation of flow and to develop the required engineering information for the particular case represented by the model. This technique has been used to calculate pattern performance in connection with several field projects involving water injection into oil reservoirs. This type of analysis involves sequentially (I) calculating the pres-' sure distribution in the reservoir and (2) tracking the progress of the interface between the displaced and displacing fluids in a step-wise manner to provide a depletion history of the operation. The particular analysis presented in this discussion is subject to the restriction that the mobility of the displaced and displacing fluids be equal and assumes that the fluids are incompressible and that gravity and capillary pressure do not affect the shape of the flood pattern. Complete two-dimensional .flexibility is maintained with regard to definition of reservoir rock and fluid characteristics, placement of physical restrains and boundary conditions, investigation of flow characteristics in the reservoir, and movement of the displaced and displacing fluid interface. The results of these studies indicate that the highspeed digital computer is well suited for conducting reservoir performance studies in two-space dimensions and that dependable numerical techniques are available for making such analyses. INTRODUCTION One of the problems facing petroleum engineers is that of providing more precise engineering control over reservoir operations at all stages during the depletion cycle. The most direct approach to this problem is to simulate the reservoir with some type of physical or mathematical model and use the information developed from the behavior of the prototype to predict the performance of the actual reservoir. This principle is, of course, not new, and a variety of simulation techniques have been developed to provide information on the behavior of oil and gas reservoirs. One straightforward refinement of present reservoir analysis techniques is the progress from one-dimensional reservoir models to those where the reservoir can be represented in two- and three-space dimensions. Multidimensional reservoir analyses are desired in order to show not only what is happening in the reservoir but in addition to show where in the reservoir the various phenomena of interest are occurring. Thus far, two-dimensional analyses have been confined predominantly to potentiometric and electrolytic model studies. However, when large numbers of wells are involved and where large-scale models are required, such problems become very tedious and instrumentation problems become acute." It is possible through use of numerical methods and general purpose digital computers to avoid many of the instrumentation difficulties and to relax some of the restrictions of previous techniques. The purpose of this paper is to discuss some recent experiences with the use of high-speed digital computers for predicting pattern performance in fluid injection operations and to describe the numerical procedures which have been used for making such analyses. THE NUMERICAL MODEL The problem of representing the reservoir with an appropriate numerical model is composed of at least two important considerations. The first requires that an adequate description of the reservoir mechanism be formulated in the form of the basic differential equation describing the process. The second consideration is definition of the physical characteristics of the reservoir. The latter requires that representative basic data be obtained, examined, and interpreted to provide sufficient detailed definition of the pertinent physical properties throughout the system. The Differential Equation One of the basic assumptions made for the purpose of performing the two-dimensional studies being dis-

By R. W. Parsons

wben a displacement front encounters an isolated permeability heterogeneity, it will be perturbed. The details of this perturbation will depend on the beterogeneity size and shape, the Permeability contrast, and the mobility ratio. This paper deals with the displacement front behavior for one, simple, mathematically tractable model: a circular permeability discontinuity and unit mobility ratio. The progress of an initially planar front is shown for four permeability contrasts: k2/k1 = 0, 5/2, 2, m. The front distortion is more pronounced than streamline and isopotential distortion and it approaches a constant shape downstream from the discontinuity. INTRODUCTION The impact of reservoir heterogeneities on oil displacing processes depends upon the variation of the heterogeneous property (frequency distribution), its disposition (spatial distribution), and the inherent stability of the displacement mechanism.l The interplay of these factors plus the impracticability, if not impossibility, of describing a natural rock system in any detail makes exact displacement front movements difficult to determine. It is useful, however, to be able to qualitatively visualize what heterogeneities will do to a displacement front. In this respect certain simple models that can be treated analytically are of value. The following analysis will show the displacement front behavior for one such model. MODEL DESCRIPTION The flow field is two-dimensional and infinite in extent. The fluid flow is single phase, steady state and unidirectional at x, y = 2 m, or before insertion of the discontinuity. A circular permea- bility discontinuity is placed at the origin. Porosity and thickness are the same within and exterior to the heterogeneity. This is exactly the same model that was used by Greenkorn et al. to check analytically Streamlines as determined by Heie-Shaw flow models and finite difference computer techniques.2 FRONT POSITION CALCULATIONS The results of displacement front calculations are presented in Figs. 1 through 4 for four different permeability contrasts: k2/kl = O, %, 2, m. The times corresponding to the front locations are the Same for the 0 and cases! and for the 1/2 and 2 cases. The determination of front positions at various times is based on known analytical expressions for velocity potentials ().2 and stream functions ().2 These expressions are in terms of the x and y

Jan 1, 1969

By J. N. Breston, R. V. Hughes

Conclusions drawn by previous research workers with reSPect to the relation between Pressure gradients and/or velocity and oil recovery obtained by laboratory water flood tests have been in disagreement probably due to variable procedures and unnatural conditions and materials. The Bradford Laboratory of the Pennsylvania Grade Crude Oil Association as part of its secondary recovery research program has conducted nineteen water floods on two long cores of widely differing characteristics in an attempt to clarify this relationship and make it an aid in predicting flooding pressures in the field. Unlike previous research procedures the present experiments were conducted with the aim of duplicating field conditions as closely as possible by using long unextracted consolidated cores, a live crude, and natural brines for both flooding and, connate water content. Also, the pressure gradients and flooding velocities were representative of field conditions where similar sands were being flooded. Eleven floods on one core and eight floods on the other core showed increased recoveries and lower residual oil saturation with increased flood pressure gradients and flood velocities. A marked decrease in recovery was obtained from both cores at very low flood velocities. This pressure versus recovery relationship is shown to hold up to the point of water breakthrough and also up to the 100 and 1 produced water to oil ratio point. INTRODUCTION The possibility of water flooding oil sands was suggested by Carl1 of the Pennsylvania Geological Survey in 1880. It is not known when the practice was tried intentionally for the first time, but its beneficial effects were noted in the annual production rate of the Bradford field as early as 1907. The practice was illegal in Pennsylvania until 1921. Early water floods in the Bradford field usually consisted of shooting or splitting the casing secretly to permit subsurface waters to enter the producing sands under hydrostatic head. As it was noted that the benefits of water flooding seemed to be proportional to the quantity of water dumped into the well many also began to utilize surface sources after the practice became legal. It was probably during the middle 20's before many producers realized that the pressurehead of the water upon the producing sand determined the rate and quantity of water that would enter the sand. Hence, rate and quantity of production appeared to be a direct function of input pressures. By 1927 a few producers had ventured the installation of pressure pumps in order to increase water-input rates and production through the combination of hydrostatic and hydraulic pressures. The adoption of pressure flooding and the "five-spot" drilling pattern in the Bradford field were essentially simultaneous. Water-input pressures in 1930 seldom exceeded 300 p.s.i. at the well head or 1100 p.s.i. at the sand face. Since that time, water-flood producers in the Bradford-Allegany fields have gone to higher and higher pressure until today 600 p.s.i. at the well heads is called a low pressure flood. Many high pressure floods now operate at 1300-1400 p.s.i. at the well head. The limiting and advisable pressure at the sand face has been pronounced as that pressure just under what is required to lift the overburden or to cause formation parting: According to this rule any water flood operation utilizing well-head pressures nearly equal in pounds per sq. in. to 1.1 times the average depth in feet to the top of the producing sand would be considered as a high-pressure flood. Despite the higher pressure trends in the Bradford-Allegany field operations and the results of early laboratory water flooding research, the desirability and benefits of high input pressures are still questioned by many operators, particularly in midwestern water flood operations. The present paper recounts a series of 19 laboratory water floods using two long, consolidated cores of widely differing characteristics saturated with live crude and flooded with oil field brines in an effort to simulate field conditions under • various pressure gradients and flooding velocities. For both cores, higher recoveries and lower residual oil saturations were obtained at higher pressure gradients and flood velocities. This relationship is shown to hold up to the water breakthrough point and also up to the 100 to 1 produced water to oil ratio point.

Jan 1, 1949

By R. V. Hughes, J. N. Breston

Conclusions drawn by previous research workers with reSPect to the relation between Pressure gradients and/or velocity and oil recovery obtained by laboratory water flood tests have been in disagreement probably due to variable procedures and unnatural conditions and materials. The Bradford Laboratory of the Pennsylvania Grade Crude Oil Association as part of its secondary recovery research program has conducted nineteen water floods on two long cores of widely differing characteristics in an attempt to clarify this relationship and make it an aid in predicting flooding pressures in the field. Unlike previous research procedures the present experiments were conducted with the aim of duplicating field conditions as closely as possible by using long unextracted consolidated cores, a live crude, and natural brines for both flooding and, connate water content. Also, the pressure gradients and flooding velocities were representative of field conditions where similar sands were being flooded. Eleven floods on one core and eight floods on the other core showed increased recoveries and lower residual oil saturation with increased flood pressure gradients and flood velocities. A marked decrease in recovery was obtained from both cores at very low flood velocities. This pressure versus recovery relationship is shown to hold up to the point of water breakthrough and also up to the 100 and 1 produced water to oil ratio point. INTRODUCTION The possibility of water flooding oil sands was suggested by Carl1 of the Pennsylvania Geological Survey in 1880. It is not known when the practice was tried intentionally for the first time, but its beneficial effects were noted in the annual production rate of the Bradford field as early as 1907. The practice was illegal in Pennsylvania until 1921. Early water floods in the Bradford field usually consisted of shooting or splitting the casing secretly to permit subsurface waters to enter the producing sands under hydrostatic head. As it was noted that the benefits of water flooding seemed to be proportional to the quantity of water dumped into the well many also began to utilize surface sources after the practice became legal. It was probably during the middle 20's before many producers realized that the pressurehead of the water upon the producing sand determined the rate and quantity of water that would enter the sand. Hence, rate and quantity of production appeared to be a direct function of input pressures. By 1927 a few producers had ventured the installation of pressure pumps in order to increase water-input rates and production through the combination of hydrostatic and hydraulic pressures. The adoption of pressure flooding and the "five-spot" drilling pattern in the Bradford field were essentially simultaneous. Water-input pressures in 1930 seldom exceeded 300 p.s.i. at the well head or 1100 p.s.i. at the sand face. Since that time, water-flood producers in the Bradford-Allegany fields have gone to higher and higher pressure until today 600 p.s.i. at the well heads is called a low pressure flood. Many high pressure floods now operate at 1300-1400 p.s.i. at the well head. The limiting and advisable pressure at the sand face has been pronounced as that pressure just under what is required to lift the overburden or to cause formation parting: According to this rule any water flood operation utilizing well-head pressures nearly equal in pounds per sq. in. to 1.1 times the average depth in feet to the top of the producing sand would be considered as a high-pressure flood. Despite the higher pressure trends in the Bradford-Allegany field operations and the results of early laboratory water flooding research, the desirability and benefits of high input pressures are still questioned by many operators, particularly in midwestern water flood operations. The present paper recounts a series of 19 laboratory water floods using two long, consolidated cores of widely differing characteristics saturated with live crude and flooded with oil field brines in an effort to simulate field conditions under • various pressure gradients and flooding velocities. For both cores, higher recoveries and lower residual oil saturations were obtained at higher pressure gradients and flood velocities. This relationship is shown to hold up to the water breakthrough point and also up to the 100 to 1 produced water to oil ratio point.

Jan 1, 1949

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951

By Fred H. Poettmann

The vaporization characteristics of carbon dioxide in a League City natural gas - Billings crude oil system were studied at three temperatures, 38°. 120°, and 202°F and for pressures ranging from 600 to 8,500 psi. Variation of carbon dioxide concentration up to 12 mole per cent in the composite showed no effect on the equilibrium vaporization ratios (K values) of the hydrocarbon constituents or on the K value of carbon dioxide itself. It was shown that carbon dioxide is more soluble in crudes than in distillates which is contrary to the behavior of methane. A working chart of carbon dioxide K values is presented. INTRODUCTION The study of the equilibrium vaporization ratios of mixtures of paraffin hydrocarbons has been rather thorough.2,6,7,8,9 In the past few years considerable attention has been paid to the vaporization characteristics of the so-called noncondensable gases such as nitrogen, carbon dioxide, and hydrogen sulfide in mixtures of hydrocarbons. since they usually occur to some extent in most crude oils and natural gases.1,3,4,5 Knowledge of this behavior is useful to both the production and refining phases of the petroleum industry. This paper reports the equilibrium vaporization ratios (K's) of carbon dioxide in a mixture of League City natural gas and Billings crude oil, and compares them to those obtained in a natural gas-distillate system. The equilibrium vaporization ratios for the hydrocarbon components in this system had previously been studied by Roland.' In addition to the determination of the K values for carbon dioxide, the K values for methane and ethane were also determined in order to observe what effect, if any, the presence of carbon dioxide had on these K values. The concentration of carbon dioxide was also varied in order to observe the effect of this variable on the carbon dioxide K values. EXPERIMENTAL PROCEDURE The apparatus used in this study cotlsisted of a stainless steel equilibrium cell of about 2 liters capacity. The cell was mounted on trunions permitting rocking in a thermostatically controlled oil bath. Two high pressure valves fitted with steel tubing were mounted on the top of the cell. one was used for sampling the equilibrium gas phase and the other for sampling the equilibrium liquid phase by means of an induction tube within the cell. Stainless steel tubing from the bottom of the cell led to a mercury reservoir and manifold which was connected to a free-piston type pressure gauge manufac- lured by the American Instrunlent Ctr. and to a volumetric. putrip. The temperature of the oil bath was measured by means of a ralibrated mercury-in-glass thermometer. The recorded temperatures are believed to be accurate to ±0.5 °F. The pressures are correct to 22 psi. The crude oil used in this study was stock tank oil obtained from the Wilcox formation in the Billings Field, Noble County. Okla. The natural gas was obtained from the League City Field. Galveston County, Tex. The oil was treated with anhydrous calcium sulfate in order to remove the last traces of water. To insure a supply of constant composition gas at room temperature the cylinders of League City gas were cooled to about 30°F, inverted, and the condensed liquid was allowed to drain from the cylinders. The analysis of the gas and crude are tabulated in Table I. The carbon dioxide came from Pure Carbonic, Inc., and was .stated to have a purity of 99.5 per cent or better. The procedure used to obtain samples of the equilibrium liquid and vapor was similar to that employed by others making use of the rocking type equilibrium cell.6,7,8 The equilibrium cell was evacuated and calculated quantities of carbon dioxide, natural gas, and crude oil were charged to the cell to the desired pressure. In charging the equilibrium cell an attempt was made to maintain the ratio of the natural gas to crude oil as close as possible to that employed by Roland. After the cell was charged, samples of

Jan 1, 1951