Search Documents

By P. B. Baxendell

The performance of a water-drive reservoir having a gas cap depends primarily on the movement of the gar-oil and oil-water contacts. The movement of the contacts during production depends in turn on fluid withdrawals and how the reservoir pressure changes as fluids are produced from the reservoir; that is, on how effectively the aquifer maintains pressure by replacing withdrawals. inasmuch as pressure changes, fluid withdrawals, contact positions, and produced gas-oil and oil-water ratios are interdependent, the analysis and prediction of the performance of a reservoir produced in a given way must take into account this interdependence throughout depletion. This paper presents an analysis which, within the limitations of the assumptions made, yields an engineering approach to predicting future performance based on reservoir pressure and production history. The most significant assumption is the method of extrapolation of future gas-oil and water-oil ratios. The extrapolation procedure smoothly increases both ratios to preselected values as the remaining oil column undergoes a specified decrease in thickness. This preselection is made on the busis of previous field experience in depletion of similar reservoirs under similar conditions. For computing future performance a vdumebic balance L combined with the differential equation defining pressure distribution in the aquifer to obtain positions of water-oil and gas-oil contacts. From these positions are extrapolated produced water-oil and gas-oil ratios. Reservoir performance can be investigated when oil production rates are dependent upon various factors including the performance of the reservoir itself. Examples of practical application of the procedure are included. INTRODUCTION To predict performance of water-drive reservoirs with gas caps and thin oil columns, it is necessary to describe the motions of the fluids within the reservoir during the entire production period to depletion. These motions depend on the pressure changes and on the withdrawal of oil, gas and water. Production of oil from the oil zone primarily causes water to move in to take its place; production of gas from the cap tends to cause oil to migrate into its place. In addition, the volumes of oil and gas remaining in the reservoir depend on changes in the pressure, since a decrease in pressure causes fluid expansion, gas liberation and oil shrinkage. A method of relating the future pressure to total withdrawals is used to describe the motions of fluids within the reservoir under conditions arising in possible modes of production. The analysis is based on two relationships and will reduce the problem to one amenable to digital computation. The first of these concerns the dependence of the pressure distribution in the entire aquifer furnishing the water drive upon the total reservoir withdrawals. This dependence is dictated by the permeability distribution and the extent of the aquifer; at present, such information is most readily obtained from the performance history by means of the resistancecapacitance reservoir analyzer. The second relation involves withdrawals, pressure in the reservoir, and movement of oil, gas, and water within the reservoir. The analysis is subject to certain simplifying assumptions that are necessary to permit solution of the problem. A comparison of the methods of this paper with those presently practiced is pertinent. One method of analysis is to use the reservoir analyzer to predict reservoir behavior based on aquifer characteristics determined from production history. Inasmuch as the total withdrawal rate depends upon the gas-oil and water-oil ratios, which depend in turn upon, among other things, the positions of the gas-oil and water-oil contacts, these ratios must be assumed in advance. In this paper these ratios are, instead, related to the computed positions of gas-oil and water-oil contacts. Thus, our method can be applied to problems in which the oil production rate is limited by produced gas-oil ratio, or to problems in which it is desired to determine the variable gas injection rate that will maintain the gas-oil or water-oil contact stationary. These problems cannot be worked satisfactorily on the analyzer. Another method presently used is based on a paper by Hurst.' Since his procedure is dependent upon using the solution of the heat flow equation, which requires constant permeability within the aquifer, and our procedure recognizes variations of permeability in

By J. J. Arps

n approach to the study of the electrical properlies of drilling muds, their cakes and their filtrates was made by treating the drilling fluid as a porous medium. Lluring this study it was found that the weight pararneter was omitted from previous investigations on electrical properties. Patnode in 1949 made a similar approach to ,his investigation. His work showed that there was a variation of resistivity relationships with per cent Puid present in a given suspension. This present study ex-tends Patnode's investigation to one in terms of the solids fraction of the drilling fluids. Previous investigations were made on low weight muds where the effect of the rolids fraction was obscured. The following equationr were fitted to data taken from 94 field muds and 47 laboratory muds (see Tables I and 2). The parameter, C, is a complex function of the solidr content, or mud weight, of a drilling fluid. It was also found that for all drilling muds examined. A nomograph is presented combining Eqs. I and 2. INTRODUCTION It is generally recognized that the materials present in drilling muds affect the relationship between mud filtrate and mud resistivities. However, most electric log interpretations proceed with the assumption that the mud resistivity and temperature are the only important parameters needed to obtain a filtrate resistivity or a mud cake resistivity. Lamont1 suggests that only average relationships exist between mud cake, mud filtrate and mud resistivities. He indicates that there is limited accuracy when these intercorrelated factors are used since no other measure-able variable is apparent. A previous paper by Patnode indicates that drilling fluids may be considered to be a system of fluids in solid particles which may be correlated electrically in a manner similar to rock solids. 1 he equation for the resistivity factor. F = (per cent fluid) is satisfied for a given clay type in water, but the value of m varies as different materials are used in the water. The resistivity factor is not constant for a given amount of clay. As is suspected it varies as the solids fraction conducts electrical current. Since both weighting material and clay solids are present in commonly used drilling fluids, it is expected that both materials contribute to the variation of the m factor in Eq. 3. An extension of Patnode's work is presented in this paper to show some of the interrelations which exist between the resistivities of mud cake, mud filtrate and mud. Use of these data enable a more accurate calculation of mud filtrate and mud cake resistivities if only the mud resistivity and weight are measured. EXPERIMENTAL The measurements of mud weight, fluid and solids resistivity, alkalinity, and chloride content of mud samples were made in this research. Mud samples representing drilling fluids from 94 Louisiana and Texas wells (Table 1) and 47 laboratory-constructed systems (Table 2) were used in the experimentation. The measurements of the mud properties were made both in the field and in the laboratory. The tests on the properties of weight, alkalinity and chloride ion were standard tests as described by API specifications. The resistivity measuring devices were constructed in the laboratory for this study. For details on the resistivity equipment see Fig. 1. The fluid cell

By A. Timur

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

Jan 1, 1970

By John P. Dolan, Charles A. Einarsen, GiIman A. Hill

his paper discusses how the following formation characteristics can be determined mathematically from drill-stem test pressure charts: true formation pressures, effective permeability of the entire section tested, well productivities, wellbore damage, and possible detection of barriers (faults, pinchouts, etc.). This paper also presents a practical method for im-mediate determination of eflective permeability and well-bore damage from successful double shut-in pressure tests. A list of recommendations for improving the reliability of drill-stem rest pressures is also presented. INTRODUCTION A drill-stem test is a temporary completion designed to sample the formation fluid and to establish the possibility of commercial production. Early pressure recording devices were used merely to verify proper operation of the testing tool. Until recently the accuracy of the pressure gauges has been insufficient for any reliable quantitative use of the recorded pressures. In view of the need for more reliable formation evaluation and as a result of the recent interest in exploration work employing the concept of hydrodynamic en-trapment,1,2,3,4 better pressure recording gauges are now in use. These devices can record pressure within 1 per cent above 1,000 psig and can detect differential pressures as low as 1/2 psig. In addition to formation pressure, several other reservoir characteristics can be determined from DST charts; namely, well productivity, formation permeability, wellbore damage, and the possible existence of bar- riers (faults, pinchouts, changes in permeability, etc.). This paper presents a practical method to interpret DST pressure charts for formation pressures and many other reservoir properties, a method that has been developed in analyzing approximately 4,000 DST charts during the last five years. The techniques used are a composite of published articles on drill-stem testing5,6, together with well-known pressure build-up analysis methods. THEORY It has been shown7,8 that the following equation may be used for analysis of pressure build-up curves: Po - Pw = 162.6 qu/kh log ( to + ?/?) .....(1) When this equation is applied to the curves obtained in drill-stem testing, the assumptions and boundary con-ditions are more nearly realized than in conventional flow and build-up tests on producing wells. Zak and Griffen9 have recently discussed in detail the use of this equation in analyzing DST charts. One of the problems with DST curves is the lack of reservoir data for precise analysis. Therefore, it is necessary to develop empirical rules and field methods for analyzing DST charts in quantity. For this reason, the empirical methods presented in this paper have been developed, and their derivation is found in Appendices A and B. METHOD USED FOR ACCURATE READING OF DRILL-STEM TEST CHARTS In order to apply the pressure build-up theory to DST charts, it is necessary to obtain a digitized expression for the pressure and time data recorded by the pressure gauge. These data may be either provided by the service company or interpolated by projecting a

Jan 1, 1958

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

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

Jan 1, 1970

By R. G. Jones, W. L. Kinney, R. E. Schilson, R. S. Wilson, G. A. Clark, H. Suralo

This paper describes the field opeturluns at the Fry in situ combustion project in Crawford County, 111. Field operations may be divided into three phases—-preliminary engineering studies of 1960 and 1961, ignition, and, finally, operation of the test. The Fry combustion project was a pilot test of cocur-rent in situ combustion on a 3.3-acre inverted five-spot. In the summer of 1960, air injection tests were conducted using old wells in the area and a new air injection well drilled especially for the test. In the summer of 1961, four new oil wells were drilled on a 3.3-acre pattern for the pilot test. All old wells in the near vicinity of the pilot test were abandoned to preclude any possibility of gas escaping to adjacent formations. A compressor plant was also installed at the test site in 1961. Air injection tests and other engineering studies were made. At the conclusion of these tests, ignition was accomplished with a 40 kw electrical igniter. Operation of the Fry combustion test was little different from the usual oilfield operation with the exception of the high rate of air injection. No difficulty was found with the openhole completion. There were no unusual corrosion or emulsion problems until bottom-hole temperatures exceeded 150F. Four cores were drilled into the burned-out portion of the reservoir with no mechanical difficulties. INTRODUCTION A cocurrent or conventional in situ combustion test was conducted on the Fry unit of the Marathon Oil Co., in Crawford County, Ill. Besides testing combustion in the Fry type of reservoir, the test was intended to evaluate various operating procedures. This paper presents a discussion of the Fry combustion test field operations. The geology of the reservoir and the test performance are discussed in other papers.'.' HISTORY OF THE FRY UNIT The Fry unit consists of 11 leases totaling 337 acres in an area in which oil was originally discovered in 1906. Production from this reservoir, as from all Robinson sand reservoirs in the area, was very high during the first few years. As the dissolved-gas energy was dissipated, early attempts to retard the decline consisted of vacuum, air, and gas repressuring, and recently, most of the area has been subjected to water flooding. Following early air repressuring, the Fry waterflood unit was formed and water flooding initiated in July, 1952. In the area in which the Fry combustion test was conducted, a peripheral water flood of a 160-acre area was developed as shown in Fig. 1. Between 1956 and 1960, it became apparent that the water flood was unsuccessfd. In looking for a suitable site for an in situ combustion test, the Marathon Oil Co. selected the Fry unit as possessing those features most desirable for such a research project. Field operations associated with the Fry combustion test began in the spring of 1960.

Jan 1, 1966

By R. G. Jones, W. L. Kinney, R. E. Schilson, R. S. Wilson, G. A. Clark, H. Suralo

This paper discusses the results of the Fry conventional or cocurrent in situ combustion test, which was conduct-ed in a 3.3-acre inverted five-spot. The depth of the formation was between 880 and 936 ft; the oil had a specific gravity of 28.6° APl and a viscosity of about 40 cp at the reservoir temperature of 65°F. Preceding the combustion test, air injection tests were conducted which, in conjunction with geological studies, were used to evaluate the characteristics of the reservoir. Combustion was initiated on Oct. 13, 1961, with ignition being accomplished by a 40 kw electrical heater. The test phase of the project ended on Oct. 1, 1963. During the test, the average air injection rate was 1,520.-000 scf/D. Throughout the test, production of all fluids— gas, oil, and water—was monitored. Cumulative oil production credited to the project was 100,586 bbl. The cumulative air-oil ratio was 11,500 scf/bbl oil, and the oxygen uti1ization efficiency was 87 per cent. INTRODUCTION The Marathon Oil Co. conducted a successful in situ combustion test, beginning Aug. 22, 1961, at the Fry unit, Crawford County, 111. The purpose of the project was to test the feasibility of cocurrent in situ combustion as a means of oil recovery in the Fry type reservoir. Interest in in situ combulstion as an oil recovery tool has been stimulated mostly by the existence of large reserves of heavy viscous crudes with low expected recovery, usually less than 10 per cent. These are the so-called unrecoverable reserves, and most combustion tests to date have been conducted in this type of reservoir.1-4 In contrast, the Fry combustion test was conducted in a reservoir with a relatively high gravity oil having a relatively low viscosity. This paper discusses the performance of the test. The geology of the reservoir and the field operations are discussed in separate papers.5 TEST SITE The Fry combustion test was carried out in a 3.3-acre inverted five-spot portion of a lenticular body of Robinson sandstone. Net sand was 50 ft thick, porosity averaged 19.7 per cent, oil saturation was 68 per cent of pore volume, and water saturation was 20 per cent of pore volume. The oil in place was estimated at 1,040 bbl/acre-foot, or 171,600 bbl within the 3.3-acre pattern. The water in place was 326 bbl/acre-foot, or 53,800 bbl. The oil has a specific gravity of 28.7° API and a viscosity of 40 cp at the reservoir temperature of 65°F. AIR INJECTION PERFORMANCE Air injection took place in two phases, the phases separated by ignition of the reservoir. In the pre-ignition phase, air injection tests were conducted in the summers of 1960 and 1961. These indicated that the Fry reservoir was confined, and a high return rate of injected air could be expected. This proved to be an outstanding characteristic of the project, as cumulative gas production was 95.3 per cent of the air injected. The difference between the cumulative air injected and the total gas produced can be largely accounted for by the quantity of air stored in the reservoir. An estimated 16 X 10 cf of air remained in the burned-out portion of the reservoir and an indeterminate amount of gas was stored in the unburned but pressurized reservoir. Hence, a complete material balance of gases would account for nearly 100 per cent of the air injected. At times during the test, the daily gas production rate was 98 per cent of the daily air injection rate. Fig. 1 presents the air injection history of the Fry com-

Jan 1, 1966

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

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

Jan 1, 1970

By James E. Berryman

A method is presented for calculating the pressure depletion performance of a gas-condensate reservoir from the original reservoir fluid composition using the laboratory determined pressure-production history and retrograde liquid accumulation as guides. Equilibrium vaporization ratios for this calculation were modified from published data until the calculated dew point pressure and retrograde liquid accumulation agreed with the values observed in the laboratory. The calculated results are compared with the produced gas composition determined by low temperature analysis during constant volume laboratory experiments of pressure depletion behavior. This comparison shows that the calculation procedure is adequate for determination of gas phase composition. Data from field separation tests in low temperature separators are presented which verify the applicability of published equilibrium vaporization ratios for a wide range of separator conditions. INTRODUCTION A reliable prediction of the pressure depletion performance of a gas-condensate reservoir is necessary in determining reserves and evaluating field separation methods. The predicted performance is also used in planning future operations and studying the economics of projects for increasing liquid recovery. Because the flow of fluid in a volumetric gas-condensate reservoir is usually single phase flow, the pressure-depletion per-formake is primarily a function of the phase behavior of the reservoir fluid. The properties of a gas-condensate system of major importance are: (1) compressibility factor at reservoir temperature and pressure; (2) dew point pressure; (3) composition changes of the vapor phase with pressure depletion; (4) recovery of original in-place hydrocarbons at any pressure; and (5) retrograde condensate accumulation. These properties are most generally determined by a pressure depletion-type experiment on recombined separator gas and liquid samples in the laboratory. ldeally the pressure-depletion performance can be calculated from the composition of the original reservoir fluid if the reservoir temperature and initial pressure are known.1,2,3 It has been reported "... comparison of such computations with depletion-type ex-periments on recombined samples in the laboratory have yielded results that appeared to be entirely adequate.". The success of such calculations is determined by the equilibrium vaporization ratios chosen, especially those for the higher boiling point hydrocarbons. It is the purpose of this paper to present a method for computing the change in composition of the vapor phase during pressure depletion by production using laboratory determined properties of the system to guide the calculation and to illustrate the use of the method with data from the Cockfield "D" sand of Washington field, St. Landry Parish, La. Results of field tests and laboratory analysis are included. As a result of this reservoir fluid study, a cycling program has been evaluated, the fidd has been unitized, and the cycling started in July, 1956.

Jan 1, 1958

By William Hurst

Muskat's depletion performance equation is here derived considering the expansion behavior of the reservoir hydrocarbon system and a simple fractional-flow equation. This nietkod of derivation leads logically to the two extensions that follow. The first of these is concerned with gravity segregation in a depletion-drive reservoir. The second is concerned with including an empirically determined tern? for water influx in the performance equation. The more general equation for gravity segregation when there is a primary gas cap and empirically-determined water influx is stated for completeness. These equations have been found useful in reservoir performance calculations in Eastern Venezuela. A discussion on the methods of solving these equations follows, and considers firstly the effect of taking finite intervals in the numerical integration, and sec-ondly, methods of incorporating the time functions involved in segregation in with the expansion behavior. The paper concludes with a brief general discussion on further extensions to the depletion performance equation. INTRODUCTION The two fundamental sources of energy by which oil is produced from a reservoir result from pressure depletion inside the boundaries of the reservoir and fluid encroachment across the boundaries of the reservoir. The wells in either case form low pressure outlets through which oil and gas may be produced by the expansive force of the reservoir fluids and associated encroaching fluids. When the reservoir pressure is higher than the bubble-point pressure of the oil, so that there is no free gas in the reservoir, these expansive forces are the only ones available for the production of oil. However, when the reservoir pressure is less than the bubble-point pressure of the oil, free gas is vaporized as the pressure falls. With both oil and free-gas phases present, the additional forces of gravity and capillarity may operate on the gas-oil system, as they have previously operated on the oil-interstitial water system. Gravity tends to segregate the free gas from the oil due to their density difference. Capillarity opposes and eventually balances gravity as the more extreme free gas and oil saturations are reached, preventing the independent move- ment of free gas until it is above a certain saturation, and the independent movement of oil when it is below a certain saturation. The type of depletion performance equation chosen for predicting the future performance of a reservoir depends on the amount of past history available. When the reservoir is somewhere past the halfway mark in depletion, some form of decline curve is often used. With less past history, material balance equations which incorporate empirical factors based on the past performance are often used. When, however, the amount of past history is small, the Muskat depletion performance equation will usually be used. The distinguishing feature of this type of equation is that empirical factors based on the over-all or macroscopic reservoir behavior are almost or entirely absent. Each parameter affecting the reservoir performance is ascribed an independent set of values based on measurements made on laboratory samples; that is, incorporating microscopic empirical factors. In establishing Muskat-type depletion performance equations, it is necessary to consider the reservoir as consisting of a number of associated blocks, in each of which the saturations and pressures may be considered uniform, and in each of which all substances have uniform pressure-volume characteristics. Thus, a primary gas cap can usually be considered as one block and an aquifer as another. Gravity segregation may be negligible for practical purposes when the rock and oil properties are adverse and/or the dip or thickness of the reservoir is too small. In this case the whole oil leg may be considered as one block, except in very large reservoirs. In very large reservoirs the fluid and rock properties may vary enough, particularly in the dip direction, for it to be necessary to divide the oil leg into a number of blocks, in each of which the relevant quantities may be considered uniform. When gravity segregation of the oil and free gas is not negligible, it is necessary to consider the space occupied by the initial oil leg as divided into two blocks, a secondary gas cap and an effective oil leg, in each of which saturations may be considered to be uniform. The total volume of these two blocks is thus constant, but the secondary gas cap grows continuously at the expense of the oil leg. Muskat1 derived depletion performance equations for the basic case of an oil reservoir with closed boundaries and without segregation, and for the case of an oil

By M. Sheffield

This paper presents a technique lor predicting the flow of oil, gas and water through a petroleum reservoir. Gravitational, viscous arid capillary lorces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. INTRODUCTION Since the introduction of high speed computing equipment, one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance, 193-5,7,8 The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. DISCUSSION The method of predicting three-phase compressible fluid flow. in a reservoir can be summarized briefly by the following steps. 1. The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. 2. Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations.1,4,8 The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary forces are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and differenceequations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method of solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. The section

Jan 1, 1970

By V. J. Kniazeff, S. A. Naville

The problem of unsteady-state condensate-gas flow through porous media leads to a set of second-order non-linear partial differential equations. Such a set of equations is numerically solved in the case of radial two-phase flow around a well, taking into consideration both the thermodynamical properties of the fluid and the mechanical properties of the reservoir. The fluid properties, reflecting the PVT relationship of the gaseous and liquid phases, are expressed by using the partial specific masses of the two main separator products in these phases. The flow properties of the reservoir rock are expressed by the generalized Darcy's law for the liquid phase and by a quadratic relationship between the rate of flow and the pressure gradient for the gaseous phase. The numerical solution of the equations for pressure and saturation US radius and time is worked out through programs written for a computer. The evolution of bottom-hole pressures, well productivities or- deliverabilities and effluent compositions with the depletion of the reservoir is easily derived. The application to the Saharian gas-condensate field Hassi R'Mel led to a better understanding of the drainage mechanism. A zone of fairly high liquid saturation develops around the wells, reducing the effective permeability, and represents a loss of condensible products in addition to the PVT-like retrograde condensation. lnside this zone, near the well, the deviation from Darcy's law for the flow of the gaseous phase governs the well deliver-ability. A back-pressure test has been computed and correlates with the field results. INTRODUCTION Two-phase flow of volatile hydrocarbons, like condensate gas or light crude oil, may be treated as the flow of a binary mixture by an arbitrary division of the chemical components into two groups. This is translated into two equations of mass continuity, which constitute a set of relation- ships for the pressure and the saturations vs the space coordinates and the time. The equations contain the laws governing the composition and the motion of the phases. The problem so defined is solved with the assumption that the compositions of the phases at any pressure are respectively the same as those observed in a PVT measuring cell under differential liberation. In a first series of computations, it was assumed that the flow obeys the generalized Darcy's law. A satisfactory representation of the retrograde condensation around the well was thus obtained. In addition, the trend toward decreasing effective permeabilities was obtained, and the computed composition of the effluent checked the laboratory values. However, it has not been possible within this basic assumption to represent the non-linear relationship between the production rate and the bottom-hole pressure drawdown as observed for gas wells in the field. Following the advice of A. Houpeurt it was decided to consider the relative permeability to gas as a function of the velocity of the gas Phase.1 The necessary physical determinations were made by E. COstaséque using a method devised by A. Houpeurt and R. Iffly. As numerical processing of the equations progressed, several difficulties were encountered which were overcome through collaboration with the computer manufacturer. This mathematical model of two-phase flow in porous media had been primarily intended for and extensively applied to the case of the Hassi R'Mel gas-condensate field, operated by SN Repal for SEHR, a joint subsidiary of SN Repal and CFP (A). The programs have also found their applications to forecast the behavior of several fields in the Sahara area containing light and volatile hydrocarbons. BASIC EQUATIONS We will consider a zone in the porous reservoir where the flow properties and the in situ composition of the fluids can be assumed to be uniform. A part of the equations of transient flow can be written with the specific gravity of the fluids being taken equal to the sum of the contributions due to the

Jan 1, 1966

By L. J. Health

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

Jan 1, 1966

By N. Mungan, J. L. Thompson

Laboratory displacement tests were performed to study oil recovery elficiency by gravity drainage in fractured systems under miscible conditions. The porous media used were cylindrical Berea and Blue Jacket sandstone cores containing a number of well-defined, artificially formed, vertical and subvertical fractures. Butane and Soltrol 130 were the two miscible fluids used. The purpose of the study was to examine the influences of the displacement rate, fracture density, fracture orientation, fracture permeability, matrix permeability, crossflow, core length and connate water on the oil recovery. It was found that displacement rate, matrix permeability and the subvertical fractures affected oil recovery most. The critical flow rate, based on the matrix permeability, was found to be a significant factor in the process. For displacement rates below the critical flow rate, the oil recovery efficiency appeared to be unaffected by the density of the subvertical fractures. At the high displacement rates, the fracture density becomes important, with the recovery being the most efficient in the core having the greatest number of the subvertical fractures. The magnitude of the fracture permeability, the fracture orientation, the core length and the connate water have little effect on the oil recovery efficiency. INTRODUCTION Gravity drainage under miscible conditions from relatively thick reservoirs can be a very efficient recovery process, especially at low flow rates where the gravity forces are dominant and, consequently, the adverse viscous fingers associated with the unfavorable viscosity ratio are minimized or eliminated. Such a process might involve injection of an LPG or some solvent bank at the crest and then driving the bank with dry gas. For some favorable combination of reservoir temperature, pressure and oil composition, there may be no need to inject any solvent bank, since enrichment of the dry gas by the light ends of the crude creates an in-situ solvent bank. Slobod and Howlett2 have studied the effects of gravity segregation in vertical unconsolidated porous media under miscible conditions in the laboratory. The main variables in their study were the viscosity ratio, the density differences and the rate of flow. In this study, the objective was to find the influence of fractures on gravity drainage under miscible conditions. A greater number of related variables were also studied. LABORATORY STUDY CORE DESCRIPTION Displacement tests were performed mostly in cylindrical Berea sandstone cores. In a few cases in which a much lower matrix permeability was desired, Blue Jacket sandstone cores were used. The vertical fractures were formed by cutting the cores with a thin circular blade or by parting them along a bedding plane. The vertical fracture plane always contained the axis of the cylindrical core. The subvertical fractures, cut with a saw, were inclined 45' from the horizontal. The line of intersection between the planes of the vertical and subvertical fractures was parallel to the circular faces of the cores in all cases except one, as depicted by VF3HF-2 in Fig. 1. Fracture geometry of the other cores is also given in Fig. 1. Tables 1 and 2 give the physical properties of the solid and fractured cores, respectively. The suffixes and prefixes shown for each core have been used throughout the paper so that the reader may discern the fracture geometry from the core numbers.

Jan 1, 1970

By J. C. Melrose, W. R. Foster, J. G. Savins, E. R. Parish

Many differences can be imagined between gas-oil flow in which the gas is supplied at the face of the core and gas-oil flow in which the flowing gas was originally dissolved in the oil. If capillary pressure characteristics and flow requirements control gas saturation distribution, the gas would be expected to be located at preferred sites within the porous medium as determined by pore sizes. On the other hand, during solution gas drive the gas first appears as bubbles through a nuclea-tion process. Nothing in self-nucleation theory specifies at which sites the first bubbles should be formed. In all probability they will be randomly distributed throughout the porous medium. Furthermore, it is not at all certain that even at low rates of production the gas will redistribute itself after nucleation to the channels normally occupied by gas in simple gas flow. Stewart, et at, have shown that at least for some limestone samples, oil recoveries could not be predicted for all rates of production using any one set of relative gas and oil permeabilities. An important factor in controlling recoveries during solution gas drive was the rate of bubble formation, higher rates giving higher recoveries. Stewart, et al, attributed the increase in recovery to a better distribution of the gas phase in heterogeneous limestone samples than is obtained by simple external gas drive. Differences in recovery from these causes were not reported for sandstone cores. In the experiments to be reported here, oil recovery, pressure and producing GOR history were measured during solution gas drive for a 5-ft sandstone core. The results were compared with predictions from the Muskat method for computing solution gas-drive behavior using external gas-drive relative permeability. The effects of changing the rate of production and oil viscosity were studied. At high laboratory rates of average pressure decline, two observations were made which would not have been predicted by Muskat's depletion theory: (1) oil recovery increased with increasing rate of production for a given viscosity oil, and (2) oil recovery increased with increasing oil viscosity for a given high rate of production. Both of these observations are explained as consequences of diffusion control of gas saturations superimposed on the normal gas-oil flow requirements, Fur- thermore, discontinuous gas phase flow appears to be significant during solution gas drive. The laboratory tests were performed at rates of average pressure decline many times greater than the maximum possible rate of average pressure decline in an actual oil field. It is, therefore, not possible to draw any direct conclusions regarding the effect of rate on recovery for the solution gas-drive mechanism under actual field conditions. However, at the lower laboratory rates, recoveries were nearly independent of rate and could be predicted by the Muskat method, using external gas-drive relative permeability data. These results suggest that at normal oilfield rates the effect of rate on recovery for the solution gas-drive mechanism is negligible. EXPERIMENTAL PROCEDURES The core material for the pressure depletion studies was Bandera sandstone from an outcrop in the Mid-Continent. This sandstone was selected because of its low permeability (about 10 md), which would permit the development of substantial pressure gradients in the corn at moderate flow rates. The core was 5-ft long and 2-in. in diameter. Its properties are listed in Table 1. Relative gas-to-oil permeability ratios were measured by an external gas-drive method2. The results are shown for a short 2-in. core and for the 5-ft core in Fig. 8. Oils used in the pressure depletion experiments were kerosene and a highly refined white oil (standard white oil No. 3) with gas-free viscosities of 1.8 and 25 cp, respectively. The gas was a naturally occurring methane from Gough field, Inglewood, Calif. The oil viscosities, gas solubilities and formation volume factors are plotted as functions of pressure at 75°F in Figs. 1 and 2. Methane viscosities and compressibilities were obtained from the literature"'. A core mounting was required which could withstand up to 2,500 psi internal pressure. This was obtained by first encasing the core completely in a plastic resin (Scotch Cast, manufactured by Minnesota Mining & Manufacturing Co.) The plastic covered core was then inserted into a steel pipe equipped with screw caps so that the plastic coating could be pressured from

By P. van Meurs, C. van der Poel

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

By D. R. Parrish, W. E. Lamoreaux, W. W. Owens

Several methods are now being used by the industry for determining the gas-oil flow characteristics of reservoir rock samples. Most of the laboratory experirnerltal rlzethods can be classified either as steady-state, unsteady-state, or Hassler-type tests. Although the validity of the results obtained by any specific method is not questiorzed, in general, the complexity of these methods leaves much to be desired. This investigation dealt primarily with the development of a simplified unsteady-state (gas drive) laboratory test method for measuring kg/ko re1ationships of small core samples. A procedure for selecting uniform and representative sandstone rock samples was developetl since it was observed that the reliability of the test resu1ts was dependent upon having a sample which was uniform in porosity and permeability. The importance of factors known to influence the results of this type test, such as oil viscosity, gas expansion, and pressure differential were eva/uated. The method of calculating the test results was found to give valid data over a wide range of gas expansion. The limits of pressure differential and oil viscosity which would insure reliable results were defined. In cornparisorz with other presently used testing methods, the gar drive method requires considerably less apparatus, is less time consuming, and is more readily adaptable to a routine type core test. INTRODUCTION The detcrmination of the gas-oil flow characteristics of reservoir rocks has been given concentrated study by the oil industry for the past decade. Both analytical and experimental methods have been considered. The analytical methods have generally relied on a simplified characterization of porous materials in order that the fundamental laws of fluid flow and capillarity could be utilized. More direct experimental methods have been proposed by numerous production research laboratories with the primary goal being a simple and inexpensive, yet reliable method for measuring the gas-oil flow characteristics of reservoir rock samples. The experimental methods which have been reported in the literature generally can be grouped into one of the following classifications: 1. Steady state1,2,3,4 in which a fixed ratio of flowing fluids is forced through the test sample until equilibrium conditions of saturation and pressure gradient are established. 2. Hassler or stationary liquid method5,6,7 in which a uniform liquid saturation is attained, usually through capillary desaturation, then gas or oil is flowed through the sample at a pressure drop insufficient to disturb the pre-established saturation. :. Non-steady state method1,8,9,10 In which only the gas phase is injected into the sample, causing displacement of the oil phase and resulting in a system of continually changing average saturation and saturation gradient. The steady-state and Hassler methods have been most popular because, no doubt, of the simplicity of the actual relative permeability calculations. The problems involved in these methods were limited to apparatus design and measurement techniques. In the steady-state method, for instance, a constant flowing fluid ratio is established which, at equilibrium, results in an essentially uniform saturation over the entire length of the test sample. Therefore, the basic flow characteristics could be obtained by correlating the average saturation of the test sample with the kg/ko determined at various flowing gas-oil ratios.

Jan 1, 1957

By J. G. Richardson, J. W. Graham

The purpose of this paper is to present the results of theoretical and experimental studies of water imbibition. The imbibition processes are involved in recovery of oil from stratified and fractured-matrix formations in natural water drives and water flooding. An understanding of the role of inhibition in implementing the recovery of oil from such formations is deemed essential to proper control of these reservoirs to achieve maximum recovery. The theoretical studies involved development of the differential equations which describe the spontaneous imbibition of water by an oil-saturated rock. The dependence of the rate of water intake by the rock on the permeability, interfacial tension, contact angles, fluid viscosities and fluid saturatiorls is discussed. A few experiments were performed using core samples to determine the effects of core length and presence of a free gas suturation. The role of water imbibition in recovery of oil from a fractured-matrix reservoir by water flooding was investigated by use of a laboratory model. This model was scaled to represent one element of a frac-tured-matrix formation. Water floods were made at various rates with several fracture widths. Interpretations were made of the behavior expected in a system containing many matrix blocks. The presence of a free gas sntu.ration was found to reduce the rate of water imbibition. In the reservoir prototype of the fractured-matrix model, water imbibition rather than direct displacement by water was the dominant mechanism in the recovery of oil at low rates. INTRODUCTION Imbibition may be defined as the spontaneous taking up of a liquid by a porous solid. The spontaneous process of imbibition occurs when the fuid-filled solid is immersed or brought in contact with another fluid which preferentially wets the solid. In the process of wetting and flowing into the solid, the imbibing fluid displaces the non-wetting resident fluid. Common examples of this phenomenon are dry bricks soaking up water and expelling air, a blotter soaking up ink and expelling air and reservoir rock soaking up water and expelling oil. As increasingly better lithological descriptions have been made of the characteristics of petroleum-bearing formations, it has become obvious that imbibition phenomena which were once considered laboratory curiosities are of practical importance. For instance, in reservoirs composed of water-wet sand strata of different permeability in intimate contact, the tendency of water to channel through the more permeable stratum is offset by the tendency for water to imbibe into the tight sand and expel oil into the coarse sand. Also, in fractured-matrix formations the tendency of water to channel through the fractures is offset by water-wet matrix blocks. As some imbibition of the water into the of the largest fields in the world are fractured-matrix reservoirs, it has become increasingly important to understand all the factors involved in the imbibition process. Examples of fractured-matrix reservoirs are the Spraberry field in West Texas which produces from a fractured sandstone', the giant Kirkuk field in Iran', the Dukhan field in Qatar, Persian Gulf2, and the Masjid-I-Sula-main and the Haft-Kel fields in Southwestern Iran, which produce from fissured limestone3. Research into recovery of oil from fractured-matrix formations was stimulated by the rapid decline of oil productivity of wells in the Spraberry formation. One result of this research was the water imbibition process developed by the Atlantic Refining Co.4 Another idea was that much of the Spraberry oil could be recovered by conventional water-flooding procedures5. Subsequently, pilot floods were conducted in this field to test the feasibility of these ideas. It was felt that an understanding of the role played by imbibition processes in displacement of oil from a fractured-matrix reservoir could not be obtained from field data alone because of the many complicating factors and uncertainties involved. Therefore, theoretical and laboratory studies were undertaken to provide this understanding. Study of the equations which describe the linear, countercurrent imbibition process provided an insight into the role of various factors in the process, such as the permeability of rock and inter-facial tension. In addition to the theoretical studies, imbibition experiments were conducted with core samples to determine the effect on the rate of imbibition of such variables as core length and free gas saturation. The principal experimental studies were conducted by water flooding a scaled model of an clement of a frac-tu red-matrix reservoir to evaluate

By C. S. Brooks

Several of the commonly used methods for deter-mination of effective grain densities of minerals have been evaluated for accuracy and precision, and their applicability to various types of petroleum reservoir minerals tested. Consideration was given to the effect of consolidation, swelling clay, and wettability on these grain-density determinations. Procedures with helium displacement, a modified low-pressure Boyle's law porosimeter, and liquid displacement (water, n-decane, and n-octanol) were tested. The helium-displacement procedure was used primarily as an independent reference for which the effects of fluid-solid interaction, e.g., adsorption and wettability, were negligible. The procedures with the Boyle's law porosimeter and with liquid displacement are the methods in common use. Minerals investigated included calcite, West Texas limestone core samples, California sandstone core samples, and Wyoming bentonite. The number of samples tested, while not large, was adequate to permit application of statistical methods. Grain-density measurements on calcite (Iceland spar) provided a basis of comparison for the accuracy and precision of the procedures under investigation for a material with a high degree of chemical and physical homogeneity and for which fluid-solid interactions were at a minimum. Liquid-displacement procedures (water or n-decane) provided the greatest reliability and precision for measurement of the grain densities of the limestone tested. Gas-displacement procedures (helium or Boyle's law porosimeter) provided the greatest reliability and experimental convenience, if not the greatest precision, for the measurement of the grain densities of the swell- ing clay, Wyoming bentonite, and the shah sandstones which -contained appreciable amounts of swelling, in-terstitial clays. These experiments were designed to evaluate the applicability of the grain-density procedures in common we to some of the less tractable minerals—tight limestones, shaly sandstones, and swelling clays. The results should provide a useful basis for a selection of the method most suitable for application to the various petroleum reservoir minerals encountered and especially for use in the estimation of formation porosity. INTRODUCTION These experiments were undertaken to evaluate the accuracy and precision of several commonly used methods for determination of effective grain densities of minerals, and to test their applicability to various types of petroleum reservoir minerals. Consideration was also given to the effect of consolidation, clay, and wettability on these grain-density determinations. Procedures used included measurements with helium displacement, with a modified low-pressure Boyle's law porosimeter, and with liquid displacement (water, n-decane, and n-octanol). The minerals investigated included calcite, West Texas limestone core samples, California sandstone core samples, and Wyoming bentonite. The accurate determination of the porosity of petroleum reservoir formations, in particular those of low porosity, has received much attention. Both bulk-volume and grain-volume measurements necessary in the evaluation of porosity have been examined by many investigators with the objectives of improving both the accuracy and the precision of these determinations. In this investigation statistical criteria for selection of the most reliable and suitable procedure for use with a given type of mineral are demonstrated. The grain-volume measurement procedures appear to be the most susceptible to improvement.

Jan 1, 1958

By P. E. Baker

An experimental study of heat flow in steam flooding was carried out with steam displacing water in a plane-radial fluid-flow model. Temperature distributions in the model reservoir, overburden, and substratum were measured at frequent intervals during flooding by means of a fixed array of thermocouples. Heat losses to overburden and substratum and the division of heat between the steam and hot-water zones in the reservoir were calculated from the temperature data by performing numerical integrations over the appropriate volumes. Temperature profiles and heat distribution results are presented. Some results were (1) the fraction of injected heat that is lost to overburden and substratum, when expressed as a function of time, did not depend on injection rate; (2) a significant portion of the injected heat was contained in the hot-water zone ahead of and underlying the steam front in the reservoir. The division of heat between the steam and hot-water zones did depend on injection rate; and (3) considerable gravity override was observed. The experimental results are compared with those of several published theoretical studies of heat transfer in hot-fluid injection; approximate agreement on heat losses is found with some of the most idealized theories. A valid comparison cannot be made in regard to steam us hot-water zone; basic interpretation of the hot-water zone, in terms of a vertical steam front, is clouded by gravity override in the experiments. INTRODUCTION The experimental study reported here primarily concerns two problems of heat transfer in steam flooding that have economic and engineering importance. These are (1) the fraction of injected heat lost to the overburden and substratum, and (2) the division of the remaining heat between a zone containing steam and a zone of hot water at lower-than-boiling temperature. The volume of the steamed region, which has a direct effect on oil recovery, can be calculated readily from these results. Several valuable theoretical studies have been published on this problem, but no experimental tests of the theories have previously appeared. Tests were needed because the effects of the necessary simplifications in the theories could not be determined otherwise. The basic, untested theoretical problem is the effect of "decoupling" heat and mass transfer, and the resulting inability to include the effects of heat transfer and condensation on fluid flow. The theoretical problem solved so far has been a heat transfer problem in which fluid is present as a heat carrier and flows in a stipulated manner, either linear or plane-radial. The simplifications that have sometimes been made in regard to thermal conductivity by, for example, Marx and Langenheim,l Lauwerier2 and Thomas3 apparently do not have an important effect on heat loss calculations, although the results of the different methods disagree to a certain extent, as shown by Spillette4 and by Ramey.5 A recent paper by the author6 compared the results of experiments in hot-water flooding with several theoretical results and showed a rather poor agreement. It will be seen that the experimental heat loss results for steam flooding, however, show acceptable agreement. These comparisons are discussed in a later section. EXPERIMENTAL APPARATUS AND PROCEDURE The system studied was one in which a sand pack, initially water saturated and at a uniform temperature, was flooded with steam at a constant rate and temperature. The temperature distribution in the "reservoir" and adjacent "strata" was measured with an array of thermocouples, and heat losses were determined from these measurements. The sand pack contained no oil because the only concern of the experiments was heat flow, and it was important to have the flow pattern as simple

Jan 1, 1970