Search Documents

By H. J. Ramey

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

Jan 1, 1965

By W. T. Weller

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

Jan 1, 1967

By S. T. Yuster

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

Jan 1, 1953

By R. G. Grieves, G. Thodos

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

Jan 1, 1965

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By K. E. Gray, M. S. Seth

In Part I of this work,1 equations of elasticity were formulated for transversely isotropic, axisymmetric, homogeneous, porous media exhibiting pore fluid pressure. Equations of elasticity and the thermal analogy method were used to determine transient horizontal, tangential, and vertical stresses and radial displacement in a semi-infinite cylindrical region when either a constant rate of pressure or a constant rate of flow is maintained at the wellbore. In this paper, the approach presented earlier is extended to finite reservoirs for the cases of (1) steady-state flow, (2) constant pressures at the well bore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Results of this work show that radial and tangential stress gradients are high near the wellbore but diminish rapidt; away from the well; the vertical stress gradient behaves in the same way but is less severe. Radial stresses are compressive or neutral, whereas tangential and vertical stresses may be tensile, neutral or com-pressive, depending upon the boundary conditions, the physical Properties of the system and the radial distance involved (vertical stresses are always compressive in an unbounded system1. For constant boundary pressures, both radial and tangential stresses increase with time whereas they both decrease for a closed outer boundary and constant pressure at the wellbore. The vertical stress decreases with time for both systems. For steady-state systems, radial displacement may be positive or negative, depending upon the dimensions of the system, the pressure differential and the porosity. Radial displacement may be positive or negative for a closed outer boundary but is positive for constant pressures at both boundaries. INTRODUCTION The importance, utility and complexity of a realistic appraisal of the stress state at and local to a wellbore were indicated in Part I. In this paper the analytical approach presented earlier is extended to finite, cylindrical reservoir geometry for the cases of (1) steady-state flow, (2) constant pressures at wellbore and outer boundary and (3) constant pressure at the wellbore and no flow at the outer boundary. Other than the outer boundary of the reservoir being finite, the physical system and assumptions pertinent thereto are the same as before. The reader may wish to review the mathematical development through Eq. 49 of Part I before proceeding here. SOLUTION FOR GENERAL BOUNDARY CONDITIONS THERMAL STRESSES AND DISPLACEMENTS Since normal stresses do not exist at the free boundaries of a system subjected to thermal stresses, the boundary conditions are

Jan 1, 1969

By B. H. Caudle, W. J. Bernard

Factors which influence the success or failure of a waterflood can seldom be determined in the laboratory. For this reason pilot waterfloods are initiated in a repreventative portion of the oil reservoir in question. For a pilot flood to predict quantitatively the recovery to be expected in a field-wide waterflood operation, the pilot area must behave as though it were confined (surrounded by similar areas). In this study, laboratory fluid-flow models were used to determine the simplest pilot pattern, for particular conditions of mobility ratio and initial gas saturation, that would behave as though it were confined. Pilot patterns studied ranged in complexity from a single inverted five-spot to a grouping of nine regular five-spots. Only the innermost producing well in each pattern was studied. Model results showed that the optimum number of wells in the pattern depends upon the oil-water mobility ratio and the expected oil-bank size. Unfavorable mobility ratios will, in general, require more wells in the pilot pattern than will favorable mobility ratios. Pilot patterns in reservoirs which contain a dispersed, flowable, free gas saturation will require fewer wells than for the under-saturated case. The single inverted five-spot pattern was found to be unsatisfactory for predicting behavior of fully developed waterfloods. In particular, it is possible that, in reservoirs which contain a flowable, dispersed gas phase, the oil bank will never be observed at the producers due to the large amounts of free gas which continue to be produced with the oil. INTRODUCTION One method which has been used to predict the performance of a waterflood is the pilot flood. The pilot waterflood is a flood which involves only a small cluster of the reservoir wells and is located in a small, representative portion of the reservoir. The object is that oil produced from the pilot can, in some way, be related to the oil recovery to be expected from a field-wide expansion of the waterflood. However, these field pilot waterfloods have often been unreliable in the prediction of oil recovery in a fully developed waterflood. This unreliability has also been demonstrated in several laboratory studies of pilot floods. Some of the investigators have shown that there are situations in which the pilot flood oil production is far too optimistic with respect to the oil recovery in the fully de- veloped flood. Others4-G have shown that the pilot results can also be pessimistic, especially if the pilot waterflood is initiated in an oil reservoir which has been depleted by primary recovery and is at very low pressure. The major reason for this unreliability of pilot water-floods is the migration of fluids into or out of the pilot area. By the well-known method of images, if straight lines can be drawn to represent vertical planes of symmetry in a porous medium which contains pressure sources and sinks (injectors and producers), then these lines are invariant streamlines, or lines across which there is no potential gradient, and therefore no flow. In an actual reservoir, these lines of symmetry can never be established exactly because of reservoir inhomogeneities and irregular reservoir boundaries. However, if the reservoir is relatively large and contains wells in repetitive patterns, these lines of symmetry are commonly assumed to exist for the pattern units sufficiently far removed from the reservoir boundary. Lines of symmetry for the five-spot injection pattern are shown in Fig. 1. Each five-spot unit in this figure can be considered confined with respect to flow across its boundary. In pilot floods this is not the case. The lines of symmetry for the pilot patterns investigated in this study are shown in Fig. 2. It is obvious that the fluid within these pilots is not confined and is therefore able to migrate into or out of the pilot area. Intuitively, one can see that, if more wells are added to the pilot, the innermost unit tends to behave more and more like the confined pattern. However, there is a practical limit to the number of wells which should be placed in the pilot. This limit is usually determined by economic factors. It was the purpose of this study to use laboratory fluid-flow models to determine which of the previously mentioned pilot patterns will force the innermost producing well to behave as it would in a fully developed waterflood. Since fluid migration is influenced by initial saturation conditions and the mobility ratio, these factors were included in the study. The ultimate objective of this study was to develop data which would allow the operator to choose a pilot pattern and operating conditions that will yield a production history which can be applied directly as an estimate of the performance of each production well of the fully developed waterflood. BASIS FOR THE STUDY The basic problem of field pilot floods is the migration of the reservoir fluids into or out of the pilot area. This problem has been the subject of previously reported model studies on pilot floods. These studies have been concerned mainly with the development of arbitrary "correction factors" to be applied to the simple, unconfined pilot systems such as the single five-spot. The correction factors were intended to adjust the production history of the un-

By R. J. Sandrea, S. M. Farouq Ali

The results of an experimental and theoretical study of the effects of rectilinear impermeable barriers and highly permeable channels on the sweep efficiency and conductivity of a five-spot network are presented. The study was carried out using Hele-Shaw, and conductive sheet analogs for both. normal and inverted flood patterns. A functional relationship was developed which provides quantitative prediction of the above parameters; i.e., sweep efficiency and conductivity, for the nonhomogeneous system. To evaluate this, it is necessary only to measure the interference modulus of the obstacle using a simple expression. The results indicate that impermeable barriers, with an interference modulus greater than zero, always cause a decrease in the conductivity and the sweep efficiency of the pattern. On the other hand, permeability channels do not have a significant effect on these parameters, unless they are located along the streamlines. The study also discusses a method of formulating the difference equations pertaining to the two-dimensional Laplacian in a system containing rectilinear barriers. In addition, a technique is described for tracing the progress of a free surface, using marked particles, in a stream containing discontinuities SuCh as those occurring at the extremities of barriers. The above computational scheme was used in a digital computer program to simulate results of the Hele-Shaw analog. The computed values of sweep efficiency mere in good agreement with the experimental values. INTRODUCTION At the present time waterflooding remains the most widely used technique for the secondary recovery of oil. Among the numerous factors that affect the performance of a waterflood, the areal sweep efficiency constitutes one of the most critical parameters. It is a function only of the geometric distribution of the wells for a homogeneous medium and unit mobility ratio. Geological considerations, however, such as the presence of sealing faults or solution channels in the rock matrix, would generally exert an adverse effect on the areal sweep efficiency of a waterflood network. Surprisingly little information is presently available on the quantitative implications of such permeability interferences. It was the object of this investigation to consider the presence of rectilinear, impermeable and highly permeable interferences in an otherwise homogeneous five-spot pattern, and to analyze their effects on the resulting sweep efficiency and pattern conductivity. The experimental phase of the study was conducted using two different models: the Hele-Shaw flow cell and the field plotter. For experimental convenience, the geometric boundaries of the five-spot were assumed to constitute streamlines even in the case where the interference was asymmetrically placed with respect to a repeated five-spot pattern. A second object of the study was to develop a computational scheme for numerically calculating the pressure distribution in a five-spot flow pattern containing rectilinear, impermeable barriers, and to use this distribution to compute the sweep efficiency. EXPERIMENTAL PROCEDURE HELE-SHAW ANALOG The Hele-Shaw analog was used both as a visual model and for determining the sweep efficiency for systems involving impermeable barriers. The physical limitations of this model as a means of experimentally studying the motion of a fluid in porous medium have been adequately discussed elsewhere.2,7 Likewise, it can be easily verified that the analogy remains valid when a portion of the region between the plates contains an obstacle of the same thickness as the distance between the plates. The model used in this study was composed of

By Alton B. Cook, G. B. Spencer, F. P. Bobrowski

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By F. P. Bobrowski, G. B. Spencer, Alton B. Cook

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

Jan 1, 1951

By Frank G. Miller

This report presents developments of fundamental equations for describing the flow and thermodynamic behavior of two-phase single-component fluids moving under steady conditions through porous media. Many of the theoretical considerations upon which these equations are premised have received little or no attention in oil-reservoir fluid-flow research. The significance of the underlying flow theory in oil-producing operations is indicated. In particular, the theoretical analysis pertains to the steady, adiabatic, macroscopically linear, two-phase flow of a single-component fluid through a horizontal column of porous medium. It is considered that the test fluid enters the upstream end of the column while entirely in the liquid state, moves downstream an appreciable distance, begins to vaporize, and then moves through the remainder of the column as a gas-liquid mixture. The problem posed is to find the total weight rate of flow and the pressure distribution along the column for a given inlet pressure and temperature, a given exit pres5ure or temperature and given characteristics of the test fluid and porous medium. In developing the theory, gas-liquid interfacial phenomena are treated. phase equilibrium is assumed and previous theoretical work of other investigators of the problem is modified. Laboratory experiments performed with specially designed apparatus. in which propane is used as the test fluid, substantiate the theory. The apparatus. materials and experimental procedure are described. Comparative experimental and theoretical results are presented and discussed. It is believed that the research findings contributed in this * paper should not only lead to a better understanding of oil-reservoir behavior, but also should be suggective in regard to future research in this field of study. INTRODUCTION In recent years much time and effort has been consumed in both theoretical and experimental studies of the static and . dvnamic behavior of oil-reservoir fluids in porous rocks. Although lack of sufficient basic oil-field data, principally concerning the properties and characteristics of reservoir rocks and fluids, largely precludes quantitative application of research results to oil-field problems, qualitative application has become common practice. In effect. oil-reservoir engineering research is serving as a firm foundation for oil-field development and production practices leading to increased economic recoveries of petroleum. This province of research. however, still poses many perplexing problems. The thermodynamic behavior of two-phase fluids moving through porous media constitutes one facet of reservoir-fluid-flow research that has not received the attention it deserves. This report embodies a theoretical discussion of this subject and a description of a series of related laboratory experiments. The significance of the problem to oil field operations is indicated but in articular the report centers around a theory and method for analyzing the steady. macroscopically linear, two-phase flow of a fluid (a single molecular species) through a horizontal column of porous medium. For simplicity in showing how the thermodynamic behavior of two-phase fluids moving through porous media affects oil-reservoir performance problems, attention is focused temporarily on a particular well producing petroleum from an idealized water-free solution-gas drive reservoir, the reservoir rock being a horizontal, thin, fairly homogeneous sandstone of large areal extent confined between two impermeable strata. The flowing hydrocarbon fluid is considered to exist entirely as a Iiquid at points in the reservoir remote from the well; however. the decline in fluid pressure in the direction of the well causes vaporization of the hydrocarbon to begin at a radial distance r from the well. Upstream from r the fluid moves entirely as a liquid and downstream from r it moves either entirely as a gas or as a gas-liquid mixture depending on the properties of the hydrocarbon and on the thermodynamic process it follows during flow. The distance r would be variable under transient flow conditions. but for purposes of analysis the flow is considered to l~e steady at the particular instant of observation during the flowing life of the well of interest. If the flow were isothermal and the hydrocarbon a pure substance, the fluid would be entirely gaseous downstream from r. Thus, this isothermal flow process for a pure substance would require that the heat of vaporization be supplied at r. over zero length of porous medium, at the precise rate necessary to maintain the constant temperature. This means that the solid matrix of the porous medium (reservoir rock) and the surroundings (impermeable strata confining the reservoir rock) would have to serve as infinite heat sources. Heat-transfer requirements would be somewhat less severe for the isothermal flow of a multicorn-ponent hydrocarbon as bubble and dew points at the same temperature correspond to different pressures. In this instance isothermal conditions would be sustained without complete vaporization of the fluid over zero length of porous medium. Nevertheless. as the flow is in the direction of decreasing

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

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 C. Özgen Karacan

Degasification and methane control in underground coal mining is an important area augmenting conventional ventilation in order to prevent possible fires and explosions due to excessive methane emissions. Coal bed reservoir engineering has been applied for many years to coalbed methane production. In recent years, modeling the flow and draining methane, not only from coal bed itself but also from the overlying strata, have been recognized as important areas of expertise for supplementing underground ventilation. This paper demonstrates the applications of coal bed reservoir engineering and modeling techniques for optimizing and controlling methane emissions in longwall and continuous miner sections. Practical applications and some considerations are presented for reservoir modeling, borehole performances, and face emission predictions for various coal bed and borehole configurations. A section on field reservoir studies and well bore log analyses for methane control are also presented as a complementary way of controlling methane and optimizing ventilation.

Jan 1, 2009

By A. H. Heim

A method and apparatus for measuring gas porosities of rocks are described. The apparatus can be assembled from commercially available components. In principle, measurements are made by volume substitution at constant pressure. The maximum error is not more than 0.3 porosity per cent. Typical results are given. INTRODUCTION Determining the porosity of rock samples is one of the most important and yet most varied types of measurement in core analysis. Among the many techniques devised are the so-called "gas porosity" methods. An old and well known example is the Washburn-Bunting method.' The U. S. Bureau of Mines2-' described and later improved the apparatus for a now widely used method generally known as the "Boyle's law" method. In the present form of the Washburn-Bunting method,' the volume of air in the pores of a rock sample at atmospheric pressure is extracted and then collected in a graduated burette at atmospheric pressure. The volume of air is read directly as the pore volume of the sample. The absolute error in reading the collected volume of gas is independent of the total volume; thus, the relative error is larger when the volume is small, as it is for rocks of low porosity. In addition, the sample after measurement contains mercury, which limits its use for other analyses. The Bureau of Mines (or Boyle's law) method measures directly the solids volume of a sample from which the pore volume and porosity are derived, using a separate measurement of the bulk volume. Gas at a few atmospheres pressure is introduced into a sample chamber of known volume containing the rock sample. The pressure is accurately measured. Following, the gas is expanded into a burette at 1 atm, and the gas volume is read directly. From the initial pressure p, and the final pressure p2 and volume v,, the initial gas volume v1 is calculated using Boyle's law; that is, p1v1 = p2v2. Volume v, minus the volume of the empty sample chamber is the solids volume of the sample. The accuracy of the method is limited, unless corrections are made, by deviations of the gas from the "ideal" gas-law behavior assumed in the simple form of Boyle's law. The purpose of the present paper is to describe a method for measuring the gas porosity of a rock which avoids many of these difficulties. Gas volumes are measured directly with the same accuracy as the bulk volumes. Pressures of at least an order of magnitude larger than those of previous methods are employed to insure rapid penetration of the gas into the sample. While special equipment may be built to apply the method, the porosimeter may be constructed as well from commercially available components. For simplicity, the apparatus described will be referred to as the "Constant-Pressure gas porosimeter". THE CONSTANT-PRESSURE METHOD Fig. 1 shows schematically the arrangement of components comprising the present Constant-Pressure porosimeter. Briefly, the method is one of volume substitution and may be considered a null measurement. Omitting (for the present) some of the operational details, the method of measurement consists of the following three steps. 1. After evacuation, the volume of the measuring system (a ballast chamber, a manifold, two gauges and their connections) up to the sample chamber is filled with gas to a high pressure (- 1,000 psi). A sample of the gas at this pressure is trapped in one side of a sensitive differential pressure gauge to serve as the reference pressure for subsequent steps. 2. The evacuated sample chamber containing the rock sample is opened to the measuring system. As the gas expands into the chamber, the resulting decrease in pressure unbalances the differential pressure gauge. 3. The pressure is restored by means of a mercury volumetric pump. The volume of mercury injected exactly equals the free or void volume of the sample chamber (volume of empty chamber minus the solids volume of the rock within). From the injected volume and the known empty chamber volume, the solids volume is obtained and the porosity calculated. The pressure and the volume occupied by the gas are the same before and after opening the sample chamber. Expansion and compression of the gas are incidental operations and do not enter into the calculation of porosity. By the pressure balancing or nulling, the free volume of the sample chamber is merely substituted by an equal and measured volume of mercury. Since the measurements are at constant pressure, there are no compressibility corrections necessary for the sample chamber.

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 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 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 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