Search Documents

By G. A. Zeito

Detailed visua1 examination of outcrops was used to ob-tain data on the lateral extent of shale breaks. Thirty vertical exposures belonging to maritie, deltaic and channel depositiorral environrrrents were exatmind, surveyed and photographed. The dimensions of the outcrops ranged from 356- to 8,240-ft long and 25- to 265-ft thick. Shale breaks were found to extend laterally for significant distances. and in some sands terminates by joining other break v much more frequently than by disappearance. Consequently with regard to flaw, a gross sand consisted of both continuous and discontinuous subunits. The degree of continuity of shale breaks as well as the occurrence and spatial distribution of discontinuities were different for the three depositional environments. Statistical eva1uations were performed to determine the confidence level with which estimates derived from outcrops can be applied to reservoir sands. Results of these evaluations revealed that: (I) the lateral continuity of shale breaks in marine. sands is si~nificatit, and the estimates of lateral extent can he applied to reservoir sands with a high degree of confidence (80 to 99 per cent of the shale breaks continued more than 500 ft, with a confidence of 86 per cent); and (2) the tendency for adjacent shale breaks to converge upon each other over small distances in deltaic and channel sands is highly significant (62 to 70 per cent of the shale breaks converged in less than 250 ft, with a confidence of 50 per cent), hut the probable magnitude of the resulting sand discontinuities cannot yet he predicted with adequate confidence. INTRODUCTION Almost all of the efforts devoted to characterization of the variable nature of reservoir sands have been focussed on permeability variations. Among the widely used concepts that have emerged from these efforts are those of stratified permeabilities, random permeabilities, and communicating and noncommunicating layers of different permeabilities. This study is concerned with the presence of interbedded shales and silt laminations. These features are impermeable or only slightly permeable to flow. Therefore, knowledge of the extent to which they continue laterally and the manner in which they terminate within the bodies of gross sands is important for proper description of reservoir flow. Initial field observations made on outcrops revealed that shale breaks and the relatively thinner silt laminae have impressive lateral continuity. They appeared to divide sand sections into separate individual sand layers. Although most of the layers were continuous across the total lengths of the outcrops, some were discontinuous because the- bounding shale breaks converged. Furthermore, the discontinuous layers appeared more prevalent in channel and deltaic sands than in marine sands. Based on these initial findings, a detailed investigation was carried out to determine, quantitatively: (1) the degree of continuity of shale breaks in marine. deltaic and channel sands; and (2) the frequency and spatial distribution of discontinuities in the three environments. PROCEDURE The procedure used to obtain field data from outcrops included visual examination, surveying and photographing each outcrop. The photographs were examined carefully and important outcrop features were traced, measured and recorded. The selection of outcrops for this study was made on the basis that each outcrop should be exposed clearly to permit detailed visual examination of vertical lithology. and it should also be sufficiently long (over 200 ft) to provide useful data on the lateral continuity of lithology. Identification of the depositional environment for each outcrop was made on the basis of bedding characteristics, vertical sequence of lithology and the presence of indicative sedimentary features. Whenever possible, hand specimens of associated shales were collected to determine depositional origin. Almost one-half of the outcrops used in this study required environmental identification; the remainder had already been identified by previous investigators. Several photographs of each outcrop were usually required to cover the entire length of the outcrop. These photographs were taken from one station or several, depending on the terrain, size of the outcrop and distance to the outcrop. A Hasselblad camera, with a standard 80-mm lens and a 250-mm telephoto lens, was used. The telephoto lens permitted photographing outcrops as far as two miles away. Slow-speed films were used. either Panatomic-X or Plus-X. The final operation conducted in the field was that of surveying the outcrops. The distance of an outcrop from a point of observation was determined by a triangulation method using the plane table. The measured distance was then combined with the angle of view of the camera lens to establish a scale to be used on the photographs. Films were processed using standard processing techniques and 4.5X enlargements made. The enlargements of each outcrop were butted together to form a single panorama. Slides were also prepared on several outcrops; these were used whenever greater magnification (wall projection) was required to bring out maximum lithologic detail. The shale breaks and bedding planes in each outcrop were traced on transparent acetate film superimposed on

Jan 1, 1966

By A. Madrazo

The Standing-Katz method for predicting liquid densities of reservoir fluids has been tested using experimental data of 154 bottom-hole or recombined reservoir fluid samples. New pressure- and temperature-correction curves for the Standing-Katz type correlation are presented which improve the accuracy of prediction. INTRODUCTION The present study was undertaken to investigate the known methods of determining liquid densities and to select one of them which might yield a better correlation if additional data were employed. An appropriate correlation would incorporate the following characteristics: (1) ease of handling in the computations, (2) results within engineering accuracy and (3) data easily accessible. After careful investigation of the different methods, the Standing-Katz correlation1 appeared to have these characteristics. The Standing-Katz correlation for liquid densities of hydrocarbon systems is based upon limited data. The original work, which employs apparent densities for methane and ethane, was based on data of 15 saturated crude oils in equilibrium with natural gas. Therefore, a test of the correlation using additional experimental data would serve either to validate the Standing-Katz correlation or to form a basis for corrections if such were necessary. PROCEDURE AND GRAPHICAL SOLUTION The data of 154 bottom-hole or recombined reservoir fluid samples were employed to calculate the densities at 14.7 psia and 60°F by the method described in Table 1. Volume contributions for methane and ethane were assigned from the pseudo liquid density plot (Fig. I), but volume contributions for N2, CO2 and HzS were not considered. The densities computed at reference conditions of 14.7 psia and 60°F were then elevated to their respective temperatures and pressures by the use of the correction curves proposed by Standing.' The percentage differences between the experimental values and the calculated values were determined; these samples were grouped into temperature, pressure and density ranges to determine a possible trend. It appeared that, at temperatures above 160°F, the calculated density values were consistently larger than the experimental density values; however, no obvious trend was found for the density and pressure ranges. There is a relationship between the pressure-correction and temperature-correction curves. This relationship, embodied in the temperature-correction curve, can be expressed as follows. The experimental densities of 125 samples above their bubble points were known. From these isothermal data, it was possible to determine an experimental-density difference between two pressures; and (from Fig. 2*), a calculated density difference between the same two pressures was determined. The relative effect of temperature must be considered in the calculated-density values. The effect of temperature can be determined by entering Fig. 3* at the two densities in consideration and obtaining the change in density at constant temperature. This change, or delta density, must be added

By E. L. Dougherty

A phenomenological theory for a one-dimensional unstable miscible displacement similar in type to the Buckley-Leverett model but including the effects of mixing is proposed An equation giving the fractional flow of pure solvent through an oil phase containing dissolved solvent is derived including the effects of gravity and nonbomogeneities. Numer ical integration of the resulting pair of nonlinear hyperbolic partial differential equations gives volumes of dissolved and undissolved solvent as functions of space and time. Parameters in the model which characterize mixing due to dispersion were determined from experimental data; the parmeters correlated with viscosity ratio. The results indicate that the rate of dispersive mixing is proportional to volumetric flow rate Incorporated into our equations were Koval's heterogeneity factor H. which characterizes mixing due to channeling. This allowed satisfactory predictions of displacement behavior observed in horizontal floods in nonhomo-geneous cores; but in most cases the values of H which we used were considerably larger than those used by Koval. INTRODUCTION It has been shownl,2,3 that for favorable viscosity ratios, the diffusion equation with convection satisfactorily describes the behavior of a miscible displacement in a porous media. Numerous attempts to develop a satisfactory mathematical description for the case of unfavorable viscosity ratio have been reported, but they have met with only partial success. These attempts, which are reviewed by oval! have taken two approaches: (1) development of a flow model akin to the Buckley-Leverett approach for immiscible displacement neglecting the effects of mixing5,6 and (2) simultaneous application of Darcy's Law for flow and the diffusion equation with a convection term for mass transfer.7,8,9 We set out to construct a mathematical model of Type 1 including, though, the effects of mixing. Based upon a set of hypotheses, we derived a pair of nonlinear hyperbolic partial differential equations which describe in a one-dimensional system the combined effects of flow and mixing. To work with these equations it was necessary to develop a fractional flow formula. The combined system was integrated numerically using the method of characteristics. In our equations the mixing process is characterized by four parameters. Three of these, labeled ß, P1 and P2, account for dispersive type mixing. The fourth is the heterogeneity factor H, proposed by Koval to account for channeling due to nonhomogeneities in the porous media. Values of the parameters which would cause agreement between theory and experiment were determined by trial-and-error for miscible floods conducted in horizontal cores of both homogeneous and heterogeneous materials. Calculations were also performed for vertical floods in homogeneous cores. The purpose of this paper is (1) to present the details of the mathematical analysis, (2) to present the results of the calculations, (3) to consider what light the results shed on the mixing process in an unstable miscible displacement, and (4) to provide a firmer foundation for the correlation technique developed by Koval for predicting the behavior of unstable miscible floods. STATEMENT OF PROBLEM The problem is to construct a mathematical model which describes in one dimension the observed behavior in an unstable miscible displacement. The approach is phenomenological in that the equations are based on assumptions which violate certain physical precepts known to apply to the displacement phenomenon. However, the assumptions do allow us to account mathematically for the more essential phenomena. The work is of value if we can quantitatively predict experimental results which heretofore could not be predicted, even though the synthetic nature of the model belies complete explanation of the observations. We assume that the system is comprised of two contiguous flowing phases. One phase, which we call free solvent, has the properties of pure solvent. We call the fraction of the pore volume occupied by this phase the solvent saturation, designated s. The oil phase consists of oil containing dissolved

By G. R. Countryman, G. H. Thomas, I. Fatt

Miscible displacement in both thewetting and non-wetting phase has been studied in two-phase systems. Experimental data show that dispersion is a function of saturation and cannot be predicted from dispersion theory developed for single-phase miscible displacement. Some of the effect of saturation on dispersion is a result of the development of dead-end pores or a dendritic pore structure in multiphase systems. The presence of dead-end pores or dendritic structure is confirmed by pressure transient studies. INTRODUCTION Displacement by flooding with a miscible liquid is a possible means for recovering the estimated two-thirds of the oil that remains behind after primary production. The immense economic importance of a process that can recover such large quantities of oil has led to extensive laboratory studies of miscible displacement. As usual in production research, most laboratory studies of miscible displacement have not attemped to reproduce all of the conditions existing in a petroleum reservoir. In the early stages of investigation the need has been for information on the broad, basic principles of the phenomena involved in miscible displacement. The effects of overburden pressure, reservoir temperature, and wettability have been considered of secondary importance. However, one property of a petroleum reservoir which is expected to be of major importance and yet has been omitted from many laboratory studies is the presence of interstitial water. Two possible effects of interstitial water on the displacement mechanism in the hydrocarbon phase immediately come to mind. First, from the general1y accepted theory that capillarity governs the distribution of oil and water in porous rock one would expect that for water-wet rock the water will be in the small pores and oil in the large pores. A miscible displacement of oil carried out in the presence of water is operating in a pore size distribution shown in Fig. lb, whereas if the same test had been performed with only one phase present the pore size distribution is as shown in Fig. la. Although there is not yet a theory of miscible displacement which explains in detail the effect of pore size distribution, one would expect the differences between Figs. la and lb to influence the displacement efficiency. A second factor which may make a multiphase system different from a single-phase system is the presence in the multiphase system of dead-end pores or dendritic structure. Experiments of various kinds on reservoir rock have led to the belief that all pores in the network structure of a porous rock take part in conducting fluid during single phase fluid flow.1,2 There are then no dead-end pores and no fingers or dendritic structures containing stagnant fluid. In a multiphase system, however, the second phase may

By L. L. Handy

Solvent floods using slugs of solvent have been found to show continuity in behavior from the vapor pressure of the solvent to the critical pressure for the two-component driving gas-solvent system. In the pressure region between the solvent vapor pressure and the critical pressure for the gas-solvent system, the gar and solvent are only partially miscible. Although complete miscibility cannot be obtained at these pressures, complete oil recovery is possible in principle. In two-phase solvent floods the solvent is propagated tllrough the reservoir, primarily, in the vapor phase. The carrier gas requirements constitute a significant factor in the economics of the process. A qualitative theory is proposed for estimating the amount of dry gas required to move the solvent through the reservoir. The theory shows that for two-phase solvent floods the total gas needed is a minimum at the vapor pressure of the solvent and at the critical pressure for the gas-solvent system, and is a maximum at some intermediate pressure. The predictions of the theory are supported by experimental studies using methane, butane and decane or methane, propane and decane in a natural sandstone core. INTRODUCTION Previously, solvent slug processes have been found effective for oil recovery in two pressure ranges. First, conventional miscible displacements are possible at pressures greater than the critical pressure for the gas-solvent system. Second, Jenks, et al,' have shown that, at pressures slightly in excess of the vapor pressure of the solvent, a solvent slug can be propagated through a reservoir by a gas essentially insoluble in the liquid solvent. The solvent bank displaces the oil ahead of it. Both of these processes, at least ideally, are capable of recovering all of the oil in the swept regions. Slug processes for which the gas and solvent are partially miscible have not been considered; that is, those systems for which the solvent and driving gas form two equilibrium phases in which the vapor phase contains a significant amount of solvent and the liquid phase an appreciable amount of the driving gas. Welge and Johnson' have shown that the gas needed to movc a solvent slug through the reservoir increases with increasing pressure above the vapor pressure of the solvent. It will be shown that solvent slug processes can, theoretically, recover all of the oil at any pressure greater than the vapor pressure of the solvent. But the amount of gas required to move the solvent through the reservoir depends very much on the pressure and temperature. In the present study a maximum in the gas requirements was both predicted theoretically and observed experimentally. This result has not been reported previously, and would not have been predicted from the Welge and Johnson model. The gas requirements are a minimum at the pressures corresponding to the vapor pressure of the solvent and again at the critical pressure for the gas-solvent system, and are a maximum at some intermediate pressure. AN APPROXIMATE THEORY FOR TWO-PHASE SOLVENT FLOODING The differences and similarities between conventional solvent floods and two-phase solvent floods are best understood by referring to concepts developed for miscible displacement in which miscibility is generated in the reservoir. In Fig. 1(A), a ternary diagram is shown for a hypothetical gas-solvent-oil system. To be rigorous the three components should each consist of a single molecular species. The pressure for Fig. 1(A) is greater than the critical pressure for the binary gas-solvent system at the specified temperature. Diagrams of this type are the ones most frequently referred to in discussions of enriched-gas drive and miscible displacement. A limiting tie line is shown tangent to the two-phase envelope and intersecting the gas-solvent line at Point A. To obtain generated miscibility with this type system, others have shown that, for an oil of Composition D, a mixture of gas and solvent must be injected which is richer in solvent than that composition indicated by A. An oil repeatedly contacted with a gas phase richer than A changes toward a composition which would be at equilibrium with the injected mixture, that is, a composition lying on a tie line which passes through the injected-gas composition. Since no such tie line exists, the oil is enriched to the point at which it becomes directly miscible with the injected mixture. At pressures lower than the critical pressure for the gas-solvent system, other types of phase diagrams are observed. The ones of interest in this paper are for pressures greater than the vapor pressure of the solvent, but less than the critical pressure of the gas-solvent system. Such a ternary diagram is shown on Fig. 1(B). In this case, two-phase behavior is observed not only for gas-oil mixtures, but also for certain compositions in the gas-solvent system. If a gas of Composition A (a dew-point vapor) is injected, once again the original oil is enriched by successive

By A. K. Csazar, L. W. Holm

This paper describes our investigation of a post-water-flood, oil recovery process which consists of injecting a slug of propane followed by water. Also described are the results obtained by applying a modification of the process in which gas was injected ahead of the water. Under the conditions of the latter experiments, misci-bility was not achieved between the propane and gas. Preliminary experiments or) uniform, watered-out sandstone cores showed that an oil bank could be formed and produced by applying this recovery process. However, since reservoirs are not uniform in structure, the process also was applied to porous media containing irregular porosity and to stratified sand systems. As a supplenzerrt to the experinlental work, a mathernatical procedure was developed for calculating the performance of the recovery process in a bounded, layered, porous system with crossflow between layers. As a specific example, the method was applied to predict the perforrnance of the recovery process in a 6-ft long, two-layer, stratified, unconsolidated sand model for comparison with experinlental data. The calculations were programed for the ZBM 704 computer. The equations and calcula-tional procedure presented can be extended to systems containing any number of randomly distributed permeability variations or any number of parallel layers. INTRODUCTION The problem of recovering the oil that remains in a reservoir which has been waterflooded is receiving considerable attention now as an increasing number of water floods reach an economic limit. A large number of the waterflood projects are in shallow reservoirs which are at pressures below 1,000 psi. It has been demonstrated in the laboratory that post-waterflood oil can be recover-ered by miscible displacement, but the LPG-gas, miscible flood and the enriched gas drive cannot be applied effectively at pressures below 1,000 psi. Only a few reports have appeared in the literature2-4 on low pressure, partially miscible recovery methods. However, it is possible to use LPG in a partially miscible displacement process in a reservoir where pressures of 200 to 1,000 psi can be achieved. Under these Pressures and at normal reservoir temperatures, propane is miscible with the oil; but, of course, gas or water used to drive the propane slug would not be miscible with the propane. Because of the lack of complete miscibility, it has generally been concluded that excessive amounts of propane would be required to recover oil and that such a recovery method would not be economical; however, we have found that under conditions present in certain reservoirs, an imrniscible recovery process can be applied effectively. The oil saturation in reservoirs at the economic limit of waterflood projects is usually in the range of 20 to 35 per cent of the pore space." A certain portion of this oil is left trapped by water in various size pores of the rock, but a good part of this so-called "residual" oil can be present in the less permeable lenses or layers of the reservoir rock which were by-passed to some degree by the water. The oil in these permeability traps can be produced only if favorable pressure gradients are formed in the reservoirs between adjacent zones of high and low permeabilities. A low viscosity liquid, miscible with the oil in place, which is driven by water through a stratified or heterogeneous porous system can aid in the development of these favorable pressure gradients. The oil that is released thereby from the permeability traps can be recovered by the subsequent water flood. Studies were made to determine how much oil could be recovered from homogeneous and stratified cores and models, which had been water flooded, by injecting a slug of propane and driving it with water. The effect of injecting a slug of gas ahead of the water was also determined. Most of the work described herein was done with the propane-water combination; unless otherwise specified, no gas was injected. The principal objectives of the investigation were to determine (1) if an oil bank could be formed and (2) what ratio of oil recovered to propane injected would be obtained. A further objective was to develop a method for calculating fluid-flow performance in stratified systems which would account for fluid transfer between zones in hydrodynamic communication but of different permeabilities. THEORETICAL ANALYSIS In a theoretical study of the recovery process, analytical expressions were derived to calculate the pressure distribution, the fluid flux in longitudinal (parallel to layers) and transversal (across the layers) directions, and the fluid distribution at any point in the system. The equations were developed for a two-layer porous system in which it was assumed that the fluids in the system were incompressible and that capillary and gravity effects were

By M. Mortada

Nonsteady-state flow of slightly compressible liquids in porous media of non-uniform properties has been the subject of a number of recent studies. Most of these studies considered one-dimensional flow in systems with cylindrical geometry and with the properties of the porous medium varying with radial distance only. A. Houpeurt' studied the problem in its general form by examining the solutions of the differential equation, In Eq. 1, the permeability k and the porosity + are arbitrary functions of radial distance. Houpeurt concluded that, except for special forms of the functions k and the problem is not well suited for mathematical treatment, but lends itself to high-speed digital computation. Carslaw and Jaeger' discussed the solution of Eq. 1 for the case in which the permeability varies as a power of the radial distance. P. Albert3 studied analytically the pseudo steady-state solutions in systems of finite radial extent and non-uniform properties. He also obtained some numerical solutions for nonsteady-state flow using a high-speed computer. William Hurst4 examined the interference pressure drop due to a point sink located in an infinite system consisting of two permeability regions in series. In this case, the permeability changed at a radial distance r from the point sink. INTERFERENCE BETWEEN OIL FIELDS IN A NON-UNIFORM AQUIFER The interference pressure drop for oil fields located in an extensive aquifer of uniform properties is discussed in a number of places in the literature. This work considers the interference pressure drop for oil fields located in a non-uniform extensive aquifer comprising two regions of different properties as shown in Fig. 1. Region I of the aquifer extends from which is the radius of Oil Field A, to ar,,, where a>1. It has permeability k,, porosity +, and effective compressibility c,. Oil Field B is located in Region II which extends from ar, to infinity and has different properties from those of Region I, but the same formation thickness h. Region II has permeability k2, porosity and effective compressibility c2. The choice of this type of aquifer should not imply that it is of frequent occurrence. It is chosen primarily to illustrate the effects of non-uniform aquifer properties on oilfield interference. It is of interest to point out that a non-uniform radial aquifer in which both the permeability and the porosity are inversely proportional to the radial distance can be treated in a manner similar to a linear aquifer of uniform properties. The situation where only the permeability varies as the inverse of the radial distance is also amenable to analytical treatment. Both aquifers behave quite differently from a radial aquifer with uniform properties so far as the water influx and the interference pressure drop are concerned. METHOD OF SOLUTION The expression for the interference pressure drop in Oil Field B due to a constant rate of water influx in Oil Field A is developed in the Appendix. This is accomplished by the simultaneous solution of two partial differential equations describing the pressure in Regions I and 11. The two differential equations are solved subject to the following conditions: (1) at r,, the pressure gradient is constant; (2) at ar the pressure and

By G. W. Tracy, R. D. Carter

Numerical calculations were made to determine the behavior of reservoirs with high-pressure drawdown and wide well spacing where the initial productivity is low and the wells are completed by hydraulic fracturing. The two-phase flow equations were solved for the flow into a single well. This well was assumed to be producing from a reservoir with hydraulically created horizontal fractures (four different systems with fractures were studied). For comparison purposes, additional two-phase flow calculations were made assuming a reservoir with uniform rock properties. The two-phase flow results were also compared with the conventional calculation methods, which do not include the effect of saturation gradients resulting from a simultaneous flow of oil and gas which are normal to this type reservoir. It was found that the conventional methods predict (1) a high and too optimistic value of ultimate recovery, (2) a high producing rate and a high reservoir pressure at a given oil recovery and (3) a low trend of gas-oil ratio with oil recovery. Included in the two-phase flow calculations were provisions to control the oil production rate by an allowable rate and, also, by a gas-oil ratio penalty rule. For the systems with hydraulic fractures, the producing rate was controlled by the gas-oil ratio penalty rule for most of the life. This is in contrast to the system with uniform rock properties which went "on decline" almost immediately. An unexpected characteristic of the systems which included fractures was the early rise in producing gas-oil ratio from 730 cu ft/bbl to approximately 1,200 cu ft/bbl, followed by a "leveling off" before the normally expected gas-oil ratio rise began. Additional features which are a result of hydraulic fracturing are (I) greater ultimate recovery, (2) higher average producing rates and (3) a lower average reservoir pressure at a given oil recovery. INTRODUCTION Some oil fields discovered during the past few years are producing from certain volumetric ally controlled reservoirs (often referred to as solution or internal gas- drive reservoirs) which are characterized by high-pressure drawdown at the wells. Since the available pressure drawdown at a well is limited by the static reservoir pressure and the producing rate is controlled by the available drawdown, wells completed in this type of reservoir usually produce at a rate less than the allowable from the time of completion. Because of this, this type of reservoir is referred to as a low productivity reservoir. Economic considerations require the use of wide well spacing and well stimulation by hydraulic fracturing to make commercial wells in this type of reservoir. Performance predictions for volumetrically controlled reservoirs have been made using a combination of two standard equations. 1. The "Schilthuis" or "Muskat" type material balance equation is used to relate the average reservoir pressure and the cumulative oil recovery. 2. The results from the material balance equation and the productivity factor as described by Pirson' are used to relate the cumulative recovery with producing rate and time. The material balance equation assumed uniform pressure and liquid saturation conditions throughout a reservoir. The steady-state radial flow formula allows for a pressure gradient toward a well but assumes uniform liquid saturation. These calculation methods are adequate for application to reservoirs wherein the drawdown at the well to realize satisfactory producing rates is small compared to the total pressure. In low productivity, volumetrically controlled reservoirs, the pressure drawdown at the well is large cornpared to the total pressure. Although a precise number cannot be given for the magnitude of a large pressure drawdown, values in excess of 1,000 psi would definitely be included. For practical considerations, this usually occurs when the formation flow capacity is less than about 100 md-ft. However, this limit of formation flow capacity will vary with the well producing rate. The low pressure in the neighborhood of the well which results from a high drawdown causes evolution of large volumes of gas. This causes the gas saturation to be higher near the well than at a greater distance—-hence, a non-uniform gas saturation. Also, the relationship between the relative permeability to oil (K/K) and gas saturation is nonlinear but decreases approximately in an exponential way with increases in gas saturation. Because of this, the following chain reaction is established.

By D. D. Dunlop, J. R. Welker

The growing interest in the use of CO, in crude oil recovery increases the need for data on the effect of CO, on hydrocarbon physical properties. Data are presented on the solubility of CO, in various dead oils, the swelling changes in CO2-oil solutions and the effect of CO, on dead oil viscosity. This latter property shows the most pronounced effect, with viscosity reductions up to 98 per cent of the uncarbonated viscosities. An empirical method of estimating the viscosity of carbonated oils is presented. The apparatus and procedures used are described in sufficient detail to allow others to make similar studies. INTRODUCTION The effect of dissolved carbon dioxide on the swelling and viscosity reduction of specific hydrocarbon oils has been observed and recorded by a number of investigators.'- me object of this paper is to offer a means of predicting these effects for crude oils free from natural gas, using the dead state viscosity and gravity of the crude oils. The CO, solubility and swelling of numerous crude oils were determined in a visual cell at various pressure levels. The viscosity of the oils carbonated to various pressure levels was then determined by measuring the pressure drop across a capillary tube. From these data, the physical properties were correlated empirically. The resulting correlations allow the prediction of CO, solubility, swelling and viscosity reduction if the dead state gravity and viscosity of the oils are known. SOLUBILITY AND SWELLING MEASUREMENT EQUIPMENT AND PROCEDURE A high pressure visual cell was installed in a constant temperature cabinet. A test gauge was attached at the top of the cell for pressure measurement, and a line was run through the cabinet wall to a wet test meter which was used for volumetric measurement of the gas. The first step in making a test run was to put the oil in the cell up to a level about half to two-thirds of the total volume. This required about 50 to 65 ml of oil. carbon dioxide was then bubbled up through the oil for a time during which the pressure of CO2 in the cell was kept above 800 psia. Saturation of the oil with CO2 at this pressure and ambient temperature was confirmed by slowly bleeding CO2 through a valve to the atmosphere. If the oil was completely saturated with CO2, bubbles of gas would form in the oil at the first small decrease in pressure. If the oil was under-saturated, no bubbles formed until the pressure was decreased to the saturation pressure existing in the oil. If this saturation pressure was lower than that desired, more CO2 was bubbled through the oil until the desired level was reached. After saturation at ambient temperature was completed, the cabinet temperature was adjusted to the desired level and the cell was allowed to reach temperature equilibrium. After temperature equilibrium was reached, the pressure was again decreased slightly, and the oil again checked for full CO2 saturation at the cell pressure. The pressure now had changed because of the difference in solubility of the CO, in the oil at higher temperatures and the expansion of CO2 as the temperature increased. The outlet tube from the cell was then connected to the wet test meter and the CO2 was allowed to flow slowly out of the cell and through the wet test meter at ambient temperature and pressure. The water in the wet test meter had previously been saturated with CO2 at ambient temperature and pressure by allowing CO2 to flow continuously through it lor a period of several hours. The gas flow was stopped at several pressures during the run and the cell was allowed to come to equilibrium; this made possible the measurement of solubility and swelling data at the intermediate pressures. The volume of the oil in the cell was recorded at each of the equilibrium pressures in order to obtain swelling data. DATA AND RESULTS The solubility of CO2 in the oil was calculated by the relationship V — V, where R. = solubility of CO, in crude oil, cubic feet of CO, measured at 60F and 1.0 atm/ bbl of dead state oil at the temperature under which solubility was measured, V, = volume of gas released from the cell between the saturation pressure and zero pressure, corrected to 60F and 1.0 atm, cu ft, V, = volume of CO: contained in the gas space above the oil, corrected to 60F and 1.0 atm, cu ft, and V, = volume of the dead oil in the cell in bbl at the temperature of the run. The volumetric data of Sage and Lacey' were used to calculate V., from the volume of CO2 at high pressures. The swelling factor was calculated as where V, is the volume of the C0,-

By K. E. Brown, G. H. Fancher

An 8,000-ft experimental field well was utilized to conduct flowing pressure gradient tests under conditions of continuous, multiphase flow through 2 3/8-in. OD tubing. The well was equipped with 10 gas-lift valves and 10 Maihak electronic pressure recorders, as well as instruments to accurately measure the surface pressure, temperature, volume of injected gas and fluid production. These tests were conducted for flow rates ranging from 75 to 936 B/D at various gas-liquid ratios from 105 to 9,433 scf/bbl. An expanding-orifice gas-lift valve allowed each flow rate to be produced with a range of controlled gas-liquid ratios. From these data an accurate pressure traverse has been constructed for various flow rates and for various gas-liquid ratios. A comparison of these tests to Poettmann and Carpenter's correlation indicates that deviations occur for certain ranges of flow rates and gas-liquid ratios. Numerous curves are presented illustrating the comparison of this correlation with the field data. Poettmann and Carpenter's correlation deviates some for low flow rates and, in particular, for gas-liquid ratios in excess of 3,000 scf/bbl. These deviations me believed to be mainly due to the friction-factor correlation. However, Poettmann and Carpenter's correlation gives excellent agreement in those ranges of higher density. This was as expected and predicted by Poettmann. He pointed out that their method was not intended to be extended to those ranges of low densities whereby an extreme reversal in curvature occurs.l As a result of these experimental tests, correlations using Poettmann and Carpenter's method were established between the friction factors and mass flow rates which are applicable for all gas-liquid ratios and flow rates. Definite changing flow patterns do not allow any one correlation to be accurate for all ranges of flow. INTRODUCTION The ability to analytically predict the pressure at any point in a flow string is essential in determining optimum production string dimensions and in the design of gas-lift installations. This information is also invaluable in predicting bottom-hole pressures in flowing wells. Although this problem is not new to industry, it has by no means been solved completely for all types of flow conditions. Versluys,2 Uren, et al, 3 Gosline,4 May,5 and Moore, et al,6,7 were all early investigators of multiphase flow through vertical conduits. However, all of these investigations and proposed methods were very limited as to their range of application. Likewise, many are extremely complicated and therefore not very useful in the field. Only in the last decade have any significant methods been proposed which are generally applicable. The most widely accepted procedure in industry at the present time is a semi-empirical method developed from an energy balance, proposed by Poettmann and Carpenter8 in 1952. Their correlation is based on actual pressure measurements from field wells. Accurate predictions from this correlation are limited to high flow rates and low gas-liquid ratios. Although this method will be discussed in detail later, it should be pointed out that two important parameters, namely the gas-liquid ratio and the viscosity, were omitted in their correlation. The viscosity was justifiably omitted since their data was in the highly turbulent flow region for both phases, and most wells fall in this category. The gas-liquid ratio was incorporated to some extent in the gas-density term. In 1954, Gilbert9 presented numerous pressure gradient curves obtained from field data for various flow rates and gas-liquid ratios for the determination of optimum flow strings. However, no method is presented for predicting pressure gradients except by comparison to these curves. Baxendell, et al,10 proposed a method based on Poettmann and Carpenter's procedure

By F. Selig, A. S. Odeh

A second-order approximation to the exact solution of the diffusivity equation corresponding to the pressure build-up of a well producing at a variable rate is derived. This approximation is applicable when the well's shut-in time is larger than the total time elapsed since the well was first produced. The resulting equations are compact in form and easy to use. Thus, the need for Horner's' theoretically precise but rather laborious solution to the above problem is eliminated. In addition, these equations apply where the use of Horner's widely known approximate method is questionable. From a practical point of view, the reported method is best suited for analysis of drill-stem tests and short production tests conducted on new wells. INTRODUCTION The utility of drill-stem and short production tests in reservoir studies has long been recognized by the reservoir engineer. If interpreted correctly they could lead to a wealth of information upon which may depend the success or failure of reservoirs' analyses. Initial reservoir pressure and the average flow capacity are two quantities that are normally sought from a drill-stem and/or a short production test analysis. Pressures are the most valuable and useful data in reservoir engineering. Directly or indirectly, they enter into all phases of reservoir engineering calculations. Therefore, their accurate determination is of utmost importance. The flow capacity kh of the reservoir is indicative of its commercial capability. In addition, it can indicate the presence of a darnaged zone around the wellbore and, thus, the necessity for remedial measures. Of the several methods used to analyze drill-stem and short production tests, Horner's' method is by and large the most common. It applies to an infinite reservoir and or a limited reservoir where the effect of production has not been felt by the boundary. Horner's method makes use of the so-called "point-source" solution of the diffusivity equation. The point-source solution is approximated by a logarithmic function and the superposition theorem is utilized to give the familiar pressure build-up equation where is the shut-in time, q is in reservoir barrels per day and the rest of the symbols conform with AIME nomenclature. Eq. 1 was derived for a well which produced at a constant rate q from time zero to time t and was then shut in. In actuality, such a constant rate of production does not normally obtain. Therefore, a correction must be applied to Eq. 1 to account for the varying rates of production. Horner suggested two methods. The first, which results in a theoretically accurate solution, is rather lengthy and laborious and, thus, it is not suited for routine analysis. The second which has been termed a "good working approximation" is the one used by the majority of the reservoir engineers. In the second method, Eq. 1 is modified by simply introducing a corrected time t, and writing where q is the last established production rate prior to shut-in, and t, is obtained by dividing the total cumulative production by the last established rate. Horner's original paper does not give any indication that this method of correction is based on any theoretical justification. In addition, there is a question as to what constitutes the last established rate. In case of a drill-stem test some engineers use the average rate obtained by dividing the total fluid produced by the total flow time, while others calculate the average rate by dividing the total fluid produced by the last flow-period time. Obviously, different results obtain for the different flow rates used. Because of this, a simple method to the varying-rate case was developed which is theoretically sound and which defines clearly the flow rate and its associated time to be used in the calculations. The final equation arrived at is where q* and t* are a modified rate and time, respectively, and can be easily calculated. In addition, it is shown theoretically that Horner's approximate method, if used for a variable-rate case, gives the correct pressure but would not be expected to give the correct flow capacity. MATHEMATICAL ANALYSIS The general equation governing the flow of slightly compressible fluid in porous media may be written as The elementary solution to Eq. 4, representing an instantaneous withdrawal of Q units volume of fluid at the origin at t = 0, is known as the instantaneous sink

By L. G. Jones, A. S. Odeh

A theoretical developtment is presented which provides a straightforward method of handling the drawdown analysis for both oil and gas wells flowing at variable rates. In the past our inability to analyze variable-rate drawdowns has been a major obstacle in using them for formation evaluation. This method should permit wider use of these tests in the future. INTRODUCTION Pressure build-up analyses are widely used to obtain reservoir data such as the effective flow capacity of the formation. On the other hand. although pressure drawdowns have been run for a long time, they have not been used extensively to evaluate reservoirs because they are difficult to analyze. One of the difficulties in analyzing pressure drawdowns has been that the theory developed until now has required the test to be run at a constant rate. For both oil and gas wells this is a difficult requirement to satisfy. especially during the early period of a well's life when drawdown tests are usually run. In spite of this limitation, drawdown tests are often run with build-ups as a check. In this paper a theory, which handles the drawdown analysis for variable rates, is presented for both oil and gas wells. Thus, one of the biggest deterrents to running drawdowns is removed, and they should become a much more useful engineering tool in the future. THEORETICAL CONSIDERATIONS 011, RESERVOIRS The general equation used for describing the unsteady-state radial flow of slightly compressible fluids in homogeneous porous media can be written in cylindrical coordinates as where r is the radial distance in cm, t is the time in seconds, p is the pressure at r and t in atmospheres, + is the fractional porosity, p is the viscosity in cp, c is the compressibility in vol/vol/atm, and k is the permeability in darcies. There are a number of solutions of Eq. 1 for various boundary conditions. One of these is the so-called "point-source" solution, which approximates the case of an infinite oil reservoir with a well located at r = 0 and produced at a constant rate. The point-source solution of Eq. 1 is

Jan 1, 1966

By T. F. Mueller, P. A. Witherspoon

For the case of an infinite radial system operating at constant terminal rate, the reservoir engineer often uses the "point source" solution of the diffiusivity equation to study pressure interference effects. At early times and at short distances from the inner boundary these solutions are invalid. The amount of error is often not precisely defined because of mathematical difficulties. Work presented here shows that, if dimensionless time is defined appropriately, previous solutions of the pressure equation can be displayed .as a family of curves on one chart. These solutions include the point source solution (referred to in the field of hydrology as the Theis solution) and other solutions obtained with digital computer methods. With these curves, an exact evaluation of the pressure drop within a reservoir or an aquifer can be made by the engineer. Examples of field problem solutions are presented. In most reservoirs the error involved when the Theis solution is employed is often negligible; whereas, in the calculation of interference effects in an aquifer, a substantial error can occur through such an approach. INTRODUCTION Flow equations are used in petroleum engineering to study the behavior of individual wells and reservoirs. In the case of wells, the pressure response at the wellbore face is the major point of interest; whereas, in the case of reservoirs, the pressure response at the interface of the aquifer boundary is sought. To aid in such studies, the flow equations have been solved in terms of the behavior at these two inner boundaries. Only limited work has been published in regard to the pressure conditions away from these points, i.e., within the reservoir or aquifer. Theis and Mortada are among the few who have reported on this problem. The Theis approach employs the exponential integral and is valid for pressure conditions that occur some distance away from the flow disturbance. It is derived from the concept of a point source, as opposed to a flow across some finite area. The Mortada results, on the other hand, are valid at all points within the reservoir or aquifer. They are presented in terms of dimension less ratios of the radius where the pressure is desired to the radius where the flow rate is measured. Their main use, in the past, has been in aquifer studies. The published results are presented in the form of graphs that are limited to a maximum radius ratio of 64. These graphical results are cumbersome to interpolate at non-integral radius ratios, so that one may be forced to utilize a rather involved analytical expression presented by Mortada. BASIC EQUATIONS The solutions of Mortada and Theis are both based on the diffusivity equation as applied to the case of an infinite radial system subject to a constant terminal rate. The equation is obtained by combining the material balance equation with Darcy's flow equation. The assumptions implicit in the use of this equation are as follows: (1) a single fluid is present that occupies the entire pore volume; (2) the reservoir is horizontal, homogeneous, uniform in thickness, and of infinite radial extent; (3) compressibility and viscosity of the fluid remain constant at all pressures; and (4) fluid density obeys the equation Using the diffusivity equation in situations where the above conditions do not hold will result in errors. These errors (not discussed here) but only the errors which arise in the solution of the equation itself. The diffusivity equation for the homogeneous reservoir conditions cited above can be written in cylindrical coordinates, as To obtain a dimensionless equation, so a single solution can be used for applications of different porosity, permeability and fluid properties, the following transformations are usually made: After these transformations are made, Eq. 3 .can be written in dimensionless form, as MORTADA SOLUTION One solution of Eq. 6 has been given by Mortada.

Jan 1, 1966

By A. N. Crownover, H. T. Kennedy, E. P. Miesch, C. H. Bowman

The paper presents correlations of (l) molar volume of gaseous hydrocarbon mixtures with pressure, temperature, composition and properties of the C,-plus fraction; (2) shrinkage of oils during flash and differential liberation of gas, including the calculation of formation volume factor under various conditions; and (3) bubble-point pressure with temperature, composition and characteristics of C, plus. The data on which the correlations are based comprise 1,615 measurements on 900 hydrocarbon systems, including numerous systems containing nitrogen, hydrogen sul-fide and carbon dioxide. In each correlation, the number of data points covered and the accuracy is substantially greater than in previously available work. Thus, the equation yielding molar volumes of gases has an average deviation of 2.04 per cent, applied to mixtures having temperatures up to 313F and pressures up to 9,800 psis, compared to 2.37 per cent for the Benedict-Webb-Rubin equation applied to the same data, and 4.53 per cent for the method based on the law of corresponding states. he equations resented are all explicit in the dependent variable, and require no iteration on the digital cow puter. INTRODUCTION The ease and accuracy of determining the composition of hydrocarbon mixtures, compared to the difficulty of measuring their properties under reservoir conditions, makes it desirable to utilize composition as the key to physical behavior to the greatest possible extent. As a result, there are available correlations between composition, or easily measured characteristics dependent on composition, and practically every important engineering property of reservoir fluids. The task confronting us is one of finding more exact relationships between important variables rather than extending correlations to new properties. This paper describes new correlations of molar volumes of gases, formation volume factors, and bubble-point pressures with composition, temperature and where possible, pressure. Each correlation is obtained by employing a sufficiently large amount of data so the calculated properties are probably as least as accurate as the measurement on which they are based. MOLAR VOLUME OF GASES Although many equations of state have been proposed for pure gases, only a few methods are applicable to hydrocarbons at conditions comparable to those in petroleum reservoirs. Still fewer arc useful in describing the behavior of mixtures, with which the petroleum engineer is largely concerned. The correlation presented here involves procedures similar to those of Alani and Kennedy The van der Waals equation is modified to make a and b functions of temperature instead of being constants for each material.This is a cubic equation, which may have either one or three real roots. The lowest root corresponds to liquid volume, while the highest applies to gas. When the procedure was applied to 703 pressure-temperature points of 164 gases of known composition and volume, the average deviation was 12.08 per cent, and the standard deviation 8.15. The above calculations were made using the constants derived for liquids by Alani and Kennedy and the mixture equation developed by them. A closer approximation to measured molar volumes is obtained by employing different sets of constants for different areas on the pressure-temperature chart, and by changing the relationship between a,,, and b,,, for mixtures and the a, and b, for individual hydrocarbons. The various sets are shown in Table 1, and the areas for which they are recommended are plotted in Figs. 1 through 7. In these figures, the area designated by zero is in the critical region for pure materials, and the values obtained for them may be unreliable. For normal mixtures, however, the area may be added to any of the bounding areas without substantial error. The new relationships employed are

Jan 1, 1966

By A. L. Barnes

Residual oil in watered-out reservoirs is a tremendous reserve which has been unrecoverable by established production methods. A study of the new recovery methods indicated that the forward combustion urocess might- recover oil from such reservoirs; however, no thermal recovery operating experience in a watered-out system was available. The Delaware-Childers pilot thermal test was undertaken to test the feasibility of thermal recovery in this watered-out reservoir. The pilot test consisted of a 2.22-acre, inverted five-spot in the 600-ft deep Bartlesville sand. The reservoir in the pilot area had a porosity of 20.6 per cent, a permeability of 118 md and an average sand thickness of 45 ft. The reservoir contains a 33' API, 6-cp oil. Combustion was started Nov. 22, 1960. The initial air injection capacity was 750 Mscf/D, but it was eventually increased to 2,000 Mscf/D. The test was surrounded by an active water flood; therefore, water production was initially high, but decreased as the heat wave moved toward the producing wells. Oil-bank arrival at an individual well was indicated by a drop in GOR and WOR, and an increase in oil production. Combustion-front arrival was evident at three of the pilot producers, and they were plugged. Cumulative oil production from pilot area wells war over 12,000 bbl. Operational difficulties were negligible and only conventional equipment was necessary. The combustion efficiency of this test averaged over 80 per cent. Results from coring showed that the leading edge of the combustion front tended to be wedge-shaped but a nearly complete sweep of the reservoir was eventually obtained. An isopach map based on evidence from 10 core holes and the existing wells showed that 126 acre-ft had been swept by the heat wave. Using this swept volume, an air requirement of 15.7 MMscf/acre-ft was calculated. It was calculated that 275 ,STB/acre-ft was consumed bv the heat wave. INTRODUCTION There is a large amount of oil remaining in reservoirs that have been water flooded. A study of ways to recover this oil showed that the forward combustion process might be applicable. Results from a number of forward combustion tests have been reported in the literature, but none of these tests were conducted in a watered-out system. The Delaware-Childers pilot thermal test was initiated in 1960 to define the operating characteristics of underground combustion in this watered-out Bartlesville sand reservoir, The purpose of this paper is to present the pilot test results in detail. FORWARD COMBUSTION PROCESS The forward combustion process consists of initiating combustion in the formation surrounding an injection well and driving this heat wave through the formation toward offset producing wells. As the combustion front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase, and recon-densed in the cooler part of the formation. These distilled liquids, water of combustion and gaseous combustion products, form a bank or three-phase region ahead of the burning front. This bank pushes mobile reservoir fluids toward the production wells. The rate of movement of the combustion front is controlled by the rate at which the nondistill able residue which serves as process fuel can be completely burned off the sand. The production performance from a heat wave conducted in a watered-out oil reservoir will differ from one conducted in a dissolved gas-depleted reservoir because of the difference in fluid saturations. The primary depleted reservoir contains connate water saturation, relatively high oil saturation and some gas saturation. The watered-out reservoir contains a highly mobile water saturation, residual oil saturation and little gas saturation. As the combustion front moves forward in either type of depleted reservoir, it drives the distillable oil, the water in place and the water of combustion ahead, and burns the nondistillable portion of the oil. The heat-wave process in a watered-out reservoir differs from the process in a

Jan 1, 1966

By B. K. Larkin

Solutions to the diffusion equation for a line source located anywhere in a region bounded by a circular discontinuity are presented. Such solutions follow automatically from the Green&apos;s function presented by Jaeger.4 Petroleum industry applications of this problem would include transient compressible fluid flow in porous media bounded by a discontinuity. Since several recent papers have treated the subject in less general form, it would appear that many workers in this area of the petroleum industry are not aware of Jaeger&apos;s work. STATEMENT OF THE PROBLEM In several trearments1,2,3 of the transient flow of compressible fluids in porous media it was convenient to consider a region to be bounded by a circular discontinuity. In such a region (Fig. I), the properties in Region 1, r < a, are constant but different from the properties in Region 2, r > a. Hurst, l using the Laplace transform, obtained solutions for a line source located in the center of such a porous medium. The purpose of this work is to point out that a much earlier paper by Jaeger4 presents the solution for a line source located anywhere inside the discontinuity. The assumption is made that a line source may be used to approximate a well. In the following material it is convenient to use the terms source and well as being equivalent; however, the reader should remember the distinction in applications. SOLUTIONS Jaeger presented the Green&apos;s function for a line source, or well, in the composite region shown in Fig. 1. In the most general form, the solution for a system of j wells in Region 1, each producing at a time varying rate, qi (t), is

By L. A. Wilson, N. H. Harrison, R. R. Clay, R. L. Reed. D. W. Reed

Experiments have been performed in a linear near-adiabatic system for the purpose of extending data on reverse and forward combustion from atmospheric pressure to 1,000 psig. Results obtained from reverse combustion appear to conform qualitatively with the existing description of the process when reaction kinetics are suitably modified to account for the increase in pressure, i.e., increasing tbe pressure decreases the peak temperature and increases the combustion zone velocity. Forward combustion appears to be a fuel dominated process wherein peak temperature and combustion zone velocity are not very sensitive to changes in pressure. The moderate effects of pressure that do exist at low flux virtually disappear at high flux providing all oxygen is consumed. With this provision. increasing the pressure decreases the frontal velocity and increases the peak temperature. Results me shown graphicaIIy which demonstrate the effects of presswe on peak temperature, rate of advance, oil recovery, air/oil ratio, carbon oxides produced and temperature distributions. INTRODUCTION Numerous field tests of oil recovery using the technique of underground combustion are now in progress. Operating pressures used are always substantially in excess of atmospheric pressure. Nevertheless, the literature contains no laboratory data pertaining to the effects of pressure except for those of Martin, Alexander and Dew.1 Unfortunately, these data appear to reflect heat losses which may have obscured some of the effects of increased pressure. The underground combustion processes are exceedingly complex, and general concepts relating to them have necessarily been described in relatively simple terms. Particularly with regard to forward combustion, the interaction between multiphase fluid flow, heat transport and various rate mechanisms is so involved that to date a rigorous treatment is not available. In this paper a linear near-adiabatic system is described which was found capable of operating at pressures up to 1,000 psig and temperatures of at least 1,100 F. Using this equipment, experiments on both forward and reverse combustion were accomplished at elevated pressures and the data presented as a function of air flux or temperature. It will be assumed that the reader is familiar with the papers of Wilson, Wygal, Reed, Gergins and Henderson on forward combustion and that of Reed, Reed and Tracht3 on reverse combustion, so that it is not necessary to repeat the details of the two processes. EQUIPMENT AND PROCEDURE It was decided to study the effects of pressure over the interval 0-1,000 psig using a near-adiabatic system. Since satisfactory operation of such a system requires very thin tubing walls, it proved necessary to insert the entire combustion tube and heater assembly into a high pressure jacket in such a way that the internal tube pressure could be counterbalanced by pressurizing the annulus. This meant that the numerous heater and thermocouple wires had to pass through gas-tight seals that required considerable effort to maintain effective. An unexpected difficulty arose in connection with the nature of the annulus pressuring gas. It turned out that the composition of this gas critically affected the functioning of the adiabatic control system. Operation of the equipment at elevated pressures was not particularly difficult. The problems that did arise were a direct consequence of the large numbers of connections and seals involved. EQUIPMENT The equipment consisted of the combustion tube which fitted within the pressure jacket. This jacket was in turn mounted on a wheeled carriage which facilitated movement and provided for rotation of the

By F. E. Crane, H. A. Kendall, G. H. F. Gardner

Two effects of a difference in density may be observed in miscible displacement. As the difference is increased from zero there is first a change in mode from multiple fingers to a single finger. A further increase results in a more rapid growth of the single finger. Miscible displacement with multiple fingers has been discussed by Blackwell, et al, whereas miscible displacement with a single finger has been discussed by Craig, et al. In the present paper experiments are described which bridge the gap between these extreme cases. INTRODUCTION Craig, et al,1 investigated the effect of gravity on frontal displacements for circumstances in which the gravity forces were significant. Blackwell, et a1,2 dealt with the other extreme in which gravity forces were negligible. The considerable difference between their results is illustrated in Fig. 1 where the breakthrough recoveries are plotted against the mobility ratio. The experiments described in the present paper were conducted to examine how the transition from the one extreme case to the other occurs when the gravity forces are gradually reduced. The transition is well illustrated by the production history after breakthrough. Fig. 2 shows how the number of pore volumes which must be injected at a constant rate to displace 98 per cent of the original fluid depends on the gravity force. For each mobility ratio there is a transition range, and a very significant change in the number of pore volumes occurs when this range is traversed. This range of the gravity force is physically linked with the formation of viscous fingers. A single finger is predominant for larger gravity forces; a multiplicity of fingers is formed for smaller gravity forces. Improved recovery and the occurrence of multiple fingers is concomitant with a large increase in the amount by which the fluids are mixed within the porous medium. The transition range depends on all the dimension-less groups which characterize the problem; but this dependence has not been exhaustively studied as yet and only fragmentary results are discussed in this paper.

By P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

By A. W. Adamson

A model for transport processes in liquid mixtures is discussed which supposes that the elementary act involves a position exchange between two species and that the exchange is so confined by the solvent cage as to occur nearly isosterically. The rate-determining step, thus, is likened to a bi-molecular reaction and is so treated, using absolute rate theory. The cage model has been applied to diffusion, thermal diffusion, sedimentation and viscosity, but only the first and last of these phenomena are emphasized in the present paper. The model leads to semi-empirical relationships between the absolute value for a digusion coefficient and the activation energy for diffusion, between mutual and self-diffusion coefficients and for the variation of the viscosity of a binary mixture with composition. These are discussed in relation to experimental data for various systems, including hydrocarbon mixtures. It is shown that the proposed viscosity equation and seven other commonly used ones all may be regarded as special cases of a single general relationship; a brief critical analysis is made of the basis of selection of one or the other for data fitting or interpolation. INTRODUCTION AND GENERAL THEORY The present paper covers a brief discussion of a cage model for transport processes in liquid mixtures and how this model may be useful in treating the diffusional behavior and the viscosity of such systems. Since diffusion requires the more detailed treatment, it will be taken up first, and the model then applied to viscosity. There are two types of diffusion coefficients that may be measured experimentally, apart from thermal diffusion quantities. The first is the mutual or binary diffusion coefficient, D which may be defined in terms of Fick&apos;s first law. This states that the permeation, or flux P, is proportional to the concentration gradient. In the usual experiment, P is measured relative to a frame of reference fixed with respect to the medium (e.g., the diaphragm in a diffusion cell); as a consequence, the same value of D is obtained regardless of whether P and C refer to Component 1 or to Component 2; i.e., there is only one independent mutual diffusion coefficient for a binary system. In addition to D there will be various self-diffusion coefficients. defined in terms of the gradient in labelled species i and its permeation in an otherwise uniform medium. The thermodynamic approach to mutual diffusion supposes that the actual driving force is the gradient of the chemical potential, i.e., that In the case of a dilute solution of solute, Eqs. 1 and 3 lead to the Einstein equation, If the solution is ideal and the friction coefficient is taken to be then the familiar Stokes- Einstein equation results. Mutual and self-diffusion coefficients can not be related on general thermodynamic grounds; it is necessary to invoke some additional assumptions, i.e., a model; several such have been proposed. Hartley and Crank&apos; supposed the existence of separate, intrinsic diffusion coefficients (Dl and D2) for each component, essentially corresponding to the two self-diffusion coefficients. The two flows can not be independent, however, but must be coupled through the usual restriction that there be no net volume flow. For an ideal solution. one then obtains&apos; Glasstone, et al&apos; treated diffusion in terms of absolute rate theory, but their approach otherwise resembled the previously mentioned one in that each species was considered to move with respect to the general medium in a manner determined by its individual jump distance and specific rate constant. For other than dilute solutions, a coupling of flows leading to an equation such as Eq. 6 would again be present. However, as required by Eq. 6, one does expect that the self-diffusion coefficient for the solute and the mutual-diffusion coefficient for the system become identical at infinite dilution. Lamm4 recognized that there should be three distinctive interactions in a two-component system-1-1, 1-2 and 2-2 — and, therefore, proposed three rather than two fundamental friction coefficients. Mutual diffusion resulted from 1-2 interactions only, and self-diffusion resulted from 1-2 plus either 1-1 or 2-2 interactions. Again, a collective coupling between all motions was imposed to meet the condition of no net volume flow. Laity&apos; has shown how to convert the Onsager equations to a form very similar to Lamm&apos;s. Cage Model For Diffusion Work in this laboratory on diffusion in aqueous sucrose solutions made it apparent that three, rather than two, interactions were indeed needed," but considera-