Search Documents

By H. Y. Jennings

An extensive field and laboratory experimental program was carried out to compare the waterflood behavior of carefully preserved soft and unconsolidated cores with measurements on the same cores after extraction. Results obtained from using idealized consolidated and unconsolidated porous media in which wettability could be carefully controlled were contrasted with the preserved core data. The controlled tests on idealized porous media investigated the effect of wettability, flood rate, core length, core permeability and consolidation on the displacement of high viscosity oils. It was concluded from these studies that waterfloods are more favorable when carried out with crude oil in preserved soft and unconsolidated cores than with the same cores after they have been extracted and re-saturated. Waterfloods are usually more favorable when carried out with crude oil in extracted .soft and unconsolidated cores than with refined oil of the .same viscosity in the same cores. The less favorable behavior of extracted soft and unconsolidated cores compared to preserved cores is due to alteration of the core by extraction. Preserved cores saturated with native water and oil should be used for laboratory displacement experiments because they more accurately reflect true reservoir behavior INTRODUCTION The demand for low gravity crude oil created by refinery modernization has focused attention on increasing the production of this viscous crude oil. Billions of barrels are in place in fields that are depleted, or nearly depleted, by primary production mechanisms. Since low gravity reservoirs are relatively recent, geologically, the solid matrix material is usually soft and unconsolidated sand. Such formations are also characterized by a high clay content. Evaluation of sophisticated oil recovery processes with the associated high capital investments has increased the demand for special core analysis tests on material from these soft and unconsolidated sand reservoirs. Data in this paper have resulted from an extensive field and laboratory experimental program. The initial objective vras to provide soft and unconsolidated sand cores for laboratory measurements with as little alteration as possible from their reservoir condition. A secondary objective was to compare the waterflood behavior of these carefully preserved cores with measurements on the same cores after they had been extracted and resaturated. When the comparison showed markedly different behavior the final objective was to attempt to explain the difference by making measurements on idealized porous media free of clay in which initial wettability could be carefully controlled. EXPERIMENTS MATERIALS Core Material The preserved cores used in this study were obtained with as little alteration as possible from their reservoir conditions. Special techniques were developed to satisfy this objective. The cores were cut with native crude whenever the crude had the necessary properties to satisfy the minimum drilling fluid requirements. A pure hydrocarbon chromatographic tracer was added to the crude to provide a simple, safe and inexpensive method to distinguish between oil filtrate and formation oil. Details of the tracer technique used in this study and the procedures used to handle and process the soft, unconsolidated cores have been published.' The core samples were then preserved and packaged at the well site with a dip-applied strippable plastic, or by using a rubber sleeve from a rubber-sleeve core barrel. The final step was to insure that the carefully taken and preserved samples were not altered by processing in the laboratory. Some soft and poorly consolidated sands were carefully shaped in the lab and encased in a protective mounting without disturbing their three-dimensional integrity. Many samples were so poorly consolidated that they literally flowed from the package when the seal was broken. A technique was developed to obtain cylindrical plugs from these samples by using liquid nitrogen as the drilling fluid in a conventional core-drill, drill-press assembly. The idealized porous media consisted of sintered rods of Alundum RA 139, outcrop sandstone identified as Al-hambra sandstone and packs of No. 130 Nevada sand. This sand is 95 per cent 140-200 mesh. (Physical properties of these porous media with those representative of preserved cores are in Table 1.) The natural cores and most of the idealized porous media were 1 1/2-in. in diameter and 3-in. long. The samples used to study the effect of core length were 1 1/2-in. in diameter and 12-in. long. Liquids Water and oil produced from the formations being studied were protected from the atmosphere and collected in carefully cleaned glass or plastic-lined containers. The production was free of chemical additives such as emulsion breakers and corrosion inhibitors. The crude oil was

Jan 1, 1967

By D. L. Archer, W. W. Owens

Conventional waterflooding often is uneconomic in highly fractured reservoirs because of the gross bypassing of the reservoir oil by injected water. Imbibition and pressure pulse flooding have been used in at least one fractured reservoir in an attempt to achieve better oil recovery performance. This paper presents the results of laboratory flow tests conducted on large cores to evaluate the possible applicability of these methods (particularly pressure pulse flooding) to different types of reservoir systems. Test data were obtained on both water-wet and oil-wet systems, and on systems having two widely different levels of compressibility and flow capacity. Results indicate that pressure pulse recoveries from fractured reservoirs will likely not exceed 5 to 10 per cent pore space with maximum response achieved during the first pulse cycle. Improved recovery by this method is possible from both oil-wet and water-wet reservoirs. Comparable saturation distributions during imbibition and pressure pulse production suggest that an initial pressure pulse cycle to speed production response would not interfere with subsequent imbibition flooding in water-wet reservoirs. INTRODUCTION There is a steady increase in the number of fluid injection projects being initiated each year to improve oil and gas recovery over that obtained by primary production mechanisms. Experience gained in different geographical areas and with different recovery methods is teaching that reservoir anatomy is one of the more important factors controlling the success or failure of such projects. Fractured formations appear to be the rule rather than the exception, especially in carbonate reservoirs. Conventional methods of fluid injection, whether it is waterflooding, gas injection or miscible flooding, have limited applicability to highly fractured reservoirs because of the severe bypassing of reservoir fluids. The result is early breakthrough of the injected fluid and rapidly increasing ratios of injected to in-place fluid with an undesirable effect on return on investment. Thus, innovations to conventional methods must be developed if recovery from fractured reservoirs is to be optimized. One new method proposed for application to highly fractured reservoirs is pressure pulse waterflooding. Although pressure pulsing with gas has been suggested' and rested in at least one reservoir; the Sohio Oil Co. is generally credited with the development and testing of pressure pulsing in conjunction with waterflooding. Publications" of Sohio's waterflood operations in the Spra-berry Trend indicate that the idea for pressure pulse flooding evolved from a critical analysis of the disappointing early performance of their flood initiated in a portion of the trend in April, 1961. This flood was planned to take advtantage of imbibiltion and high pressure gradients across the reservoir matrix blocks (unfractured Mocks of the reservoir rock which are visualized to be generally surrounded 'by the fracture system) evaluated in pilot tests by the Atlantic Richfield Co.5 and Humble Oil & Refining Co.' However, the early recovery and pressure performance of this flood indicalted that high pressure gradients induced between the fracture system and the centers of the matrix blocks were forcing water into the periphery of the blocks and temporarily interfering with the counter-current flow of oil due to imbibition. Fill-up in the blocks was occurring with re-solution of the free gas phase as pressures climbed above the original reservoir bubble point. It appeared that cessation of water injection to permit the capillary forces to become dominant and that expansion of the rock and its contained fluids during pressure reduction might aid in expulsion of oil from the rock matrix into the fractures. Subsequent production performance. aflter cessation of injection, proved this hypotheslis correct. Thus, pressure pulse waterflooding was born. Pressure pulse flooding appears to have several advan'tages over imbibition type flooding in highly fractured reservoirs. First, all wells may be used for water injection which should hasten fill-up and achieve a more rapid increase in reservoir pressure. Second, because of increased injection pressures and rates, flooding gradients will force water into the reservoir matrix more rapidly than would be achieved by imbibition alone. Third, during the pressure depletion cycle of the flood, the compressibility of the system, which has been increased by resolution of the free gas phase, provides energy for the displacement of oil at a rate greater than would correspond to countercurrent flow during imbibition. Fourth, during depletion all wells may be put on production, thus contributing to higher oil withdrawal rates on a reservoir-wide basis. One possible disadvantage of pressure pulse flooding as compared to imbibition floodling, however, is that the outer periphery of each matrix block is flooded to a residual or near residual oil saturation during the Water lnjection stage. This zone of reduced permeability to oil m'ay offer greater restriction to the production of

Jan 1, 1967

By J. Lohrenz, R. J. Francis, G. C. Clark

Compositional material balance methods have been extended to solution gas-gas cap-water drive reservoirs. Partial or complete pressure maintenance as well as reservoir repressuring by gas and/or water injection may be considered. Increrments of reservoir withdrawals and injections equivalent to a given time increment are used as the basic iterative step to allow the extension. An example field study shows how the reservoir volumes. extent of the water drive, and relative permeability ratio function may be estimated. Knowing these, the methods can be used to predict reservoir behavior. INTRODUCTION The compositional material balance is similar to conventional material balances in many aspects. The assumption of a pressure-homogenous reservoir with no transients is common to all conventional material balances1-1 although numerical methods are different. The compositional material balance also contains this assumption. Inclusion of phase equilibria directly in reservoir calculations is the key difference and improvement in the compositional material balance. A continuous account of the compositions of the reservoir fluids and effluents is determined. Calculations of fluid phases and properties can be based on these compositions. Excepting fluid phases and properties determined in this manner, the compositional material balance uses procedures exactly analogous to conventional methods. Brinkman and Weinaug5 and Jacoby and Berry6 have given compositional material balances for simple solution gas drive reservoirs; subsequently Jacoby and Berry' extended their procedure to include pressure maintenance by gas injection at constant reservoir pressure. This paper extends compositional material balance techniques to a more general combination drive reservoir with various types of pressure maintenance operations. The procedure for this extended compositional material balance and an example of its field application are presented. PROCEDURE The numerical calculations of this compositional mate- rial balance are performed on a high-speed digital computer. The most comprehensive case the present compositional material balance program can handle is a solution gas-gas cap-water drive reservoir with (1) water production, (2) water influx, (3) separate or commingled production from the oil and gas zones, (4) any surface processing arrangement including LPG and gasoline production, and (5) gas injection into the gas and/or oil zones." This description of the procedure applies to this most comprehensive case. Appropriate deletions in the procedure can be made for more simple mechanisms. ASSUMPTIONS In common with conventional material balances, the compositional material balance makes the following assumptions: (1) the total reservoir is at a uniform pressure at a given time; (2) the volumetric ratio of gas to liquid withdrawal at reservoir conditions is directly proportional to their mobilities; and (3) instead of assuming a pressure-dependent total gas-oil equilibrium, the cornpbsitional material balance makes a sianificantly different and more rigorous assumption as follows. Thermodynamic equilibrium exists for each component in the liquid-gas mixtures in the reservoir and in surface processing. These three assumptions are fundamental in the compositional material balance procedure. Other assumptions involve arbitrary simplifications of the reservoir mechanisms. This compositional material balance uses the following simplifying assumptions: (1) no bulk mixing occurs between the gas and oil zone; (2) total mixing occurs between any injected gas and the contents of the injection zone; (3) no in situ hydrocarbon liquid is produced from the gas zone; (4) no reservoir gas is retained behind any net bottom-water influx; (5) the residual saturation of hydrocarbon liquid at reservoir conditions behind any net bottom-water influx is constant; and (6) fluids retained behind any net bottom-water influx exert no further influence on the reservoir. These six assumptions are not fundamental. Modifications of the present procedure can be made for different assumptions of the reservoir mechanisms. THE TIME INCREMENT Previous compositional material balances have used a pressure increment for simple solution gas drive reservoirs"" and an increment of injected gas for pressure maintenance.' These increments are suitable for these simple cases. They are, however, inadequate for reservoirs with time dependent drive mechanisms including water drive. A time increment removes this inadequacy and is

By R. Simon, L. Yarborough

This paper presents a correlation for predicting the critical pressures of gas-solvent-reservoir oil systems of the type encountered in gas and solvent injection processes. Experimental data from 14 systems were used to develop the correlation which expresses critical pressure as a function of composition, C7+ mol weight, and C7+ characterization factor. Critical pressures predicted by the correlation for the 14 gas-solvent-reservoir oil systems have an average deviation of 2.8 per cent of the experimental values. The correlation war used to study a wide variety of other hydrocarbon mixtures reported in the literature whose critical conditions are in the range of interest in oil recovery work (100 to 280F and 2,000 to 6,000 psia). Critical pressures predicted by the correlation had an average deviation, expressed as per cent of experimental value, of 5 per cent for 20 miscellaneous mixtures from the literature, and 6.4 per cent for I7 methane-butane-decane mixtures reported by Sage and Lucey. INTRODUCTION Predicting critical points of hydrocarbon mixtures is an important part of phase behavior calculations. If the critical point can be predicted for a specific mixture, the separation between the bubble-point and dew-point regions will be defined, and physical properties of the mixture can be calculated by using the principle of corresponding states. Numerous correlations for predicting critical points of hydrocarbon mixtures have been published. Six of these correlationsl-8 were for light hydrocarbons, and two correlations7,8 ere for high boiling-point petroleum distillates. Organick9 And Eilerts" presented correlations which combined the published data from both light and heavy mixtures, and these correlations became the only available methods for predicting critical pressures for wide boiling range mixtures. In the course of the authors' work on secondary recovery by gas and solvent injection, experimental data were obtained for 14 mixtures of dry gas, ethane-propane-butane solvents, and various reservoir oils. The critical pressures and critical temperatures of these gas-solvent-reservoir oil systems were in the range of 1,900 to 4,800 psia and 180 to 265F. The critical pressures predicted for the 14 mixtures by the correlations of Organick and Eilerts had an average deviation from the measured values of about 13 per cent. The authors felt that it was necessary to develop a better correlation for engineering calculations and the resultant correlation, which is presented on Figs. 1 and 2, predicts the critical pressures of the 14 mixtures with an average deviation of 2.8 per cent. The validity of the correlation was tested by applying the correlation to other hydrocarbon mixtures reported in the literature. This subject is discussed later in the paper. The authors feel that there is also a need for an improved critical temperature correlation for gas-solvent-reservoir oil systems. However, developing a critical temperature correlation is a significant study in itself and is not included in this paper. EXPERIMENTAL DATA The compositions and critical points obtained by the authors for 14 gas-solvent-reservoir oil mixtures are shown on Table 1. The methane concentrations in these mixtures varied from 20 to 63 mol per cent. The intermediate hydrocarbon group, consisting of ethane, propane, and butane, varied from 20 to 64 mol per cent, and the concentration of the C,, fraction from 10 to 20 mol per cent. The physical properties of the C,, fraction covered the range from light, paraffinic, nonasphaltic (Mixtures 5, 6, 7, 8, 9 and 10) to heavy, aromatic, and asphaltic (Mixtures 13 and 14).

By H. A. Deans

A three-parameter mathematical model for one-dimensional flow in porous media is developed. The objective of the model is to predict accurately the longitudinal dispersion associated with the flow of either gases or liquids over a wide range of Reynolds number. A qualitative analysis of the model parameters is given. Published experimental pulse-response curves are compared with results predicted by the model. Several general types of problems are given for which the model can be used. INTRODUCTION The problem of multicomponent single-phase flow through porous media is encountered in the study of petroleum reservoirs, gas chromatographic and ion-exchange columns, industrial fixed-bed con-tacters and elsewhere. These particular examples span the wide range of flow conditions possible; Reynolds numbers of less than l0-4 are not unusual in oil-production problems, while values in excess of 103 are common in large fixed-bed operations. The study of flow-dependent transport phenomena is complicated both by the changes in the character of the flow over this range, and by the irregularity of the flow boundaries inherent to porous media. Dispersion is one of the important phenomena known to depend fundamentally on flow conditions as well as on fluid and medium properties. As used in this paper, the term "dispersion" refers to the observed mixing of fluid elements of different composition which occurs in flow systems. The actual mechanism may be one or more of a number listed below. Only dispersion in the direction of the mean flow (referred to as axial or longitudinal dispersion, or mixing) is considered here, although lateral dispersion arises as part of certain coupled mechanisms. MECHANISMS FOR LONGITUDINAL DISPERSION A number of distinct mechanisms are known3,5,7,8 to contribute to the phenomenon of longitudinal dispersion. The more important of these are as follows. 1. Molecular diffusion in the flow direction. 2. Turbulent (cell) mixing. 3. Lateral transport processes coupled with velocity and/or residence time distributions, including: (a) "Taylor" diffusion caused by the interaction of velocity profiles in individual voids with lateral molecular diffusion; (b) separation and remixing or interdiffusion of streams having different velocities around particles; and (c) the couplipg of gross velocity profiles, caused by viscous instability or inhomogeneous porosity, with lateral dispersion. 4. Finite mass-transfer rate between a porous matrix and the flowing phase, and finite diffusion rate inside elements of the porous matrix. It is a well documented fact that these mechanisms, whether acting individually or in combinations, produce essentially similar integral effects. The observed dispersion can be described approximately by solutions of the diffusion equation with a properly chosen dispersion coefficient. Necessary conditions for equivalence between dispersion and diffusion are discussed by Klinkenberg and Sjenitzer.7 Apparently these conditions are not satisfied exactly in many real systems. "Best fit" solutions of the diffusion equation often depart appreciably from measured breakthrough curves for the usual pulse- or step-forced experimental systems. In particular, pulse response curves usually show some degree of asymmetry and "tailing" which cannot be reproduced by the one-parameter diffusion model. The three-parameter model proposed herein attempts to produce these effects in a realistic manner. Asymmetry and tailing are predicted for certain sets of parameter values, while for other sets the new model agrees with the diffusion model. VELOCITY DEPENDENCE OF DISPERSION Molecular diffusion apparently controls at sufficiently low Reynolds number in both gas and liquid systems.*l3 At Re > 20, turbulent cell mixing is dominant for gases flowing at low pressure; the data of McHenry and Wilhelm8 show that the mixing cell is approximately 1 particle diameter long in random sphere packs. on the other hand, the data

By E. J. Koval

Practical miscible displacement processes will be characterized by lingering of the solvent into the oil. The lingering process is brought on by viscosity differences, and can be accentuated by channeling and longitudinal dispersion. The effects of these factors on the efficiency of unstable completely miscible displacements are combined in what is called the K-factor method. This method, analogous to the Buckley-Leverett method, predicts recovery and solvent cut as a function of pore volumes of solvent injected. Experimental data are included and show excellent agreement with theory for a wide variety of sandstone cores and viscosity ratios. INTRODUCTION Theoretical considerations,1 laboratory experiments,2 and pilot tests3 lead to the conclusion that miscible displacements in the field will be unstable In an unstable miscible displacement, the solvent fingers through the oil. This fingering leads to early breakthrough of the solvent and an extended period during which both oil and solvent are produced. For such a system, there appear to be four principal factors which bring about or accentuate the effects of instability: longitudinal dispersion (including geometrical effects), channeling, viscosity differences and gravity differences. Other factors, such as diffusion and flooding rate, can also influence the effects of instability, but at the flooding rates considered in this report, circa 15 ft/D, they are unimportant. Longitudinal dispersion can be thought of as a spreading of the solvent front caused by the presence of microscopic inhomogeneities. Channeling of the solvent occurs when a porous medium has macroscopic inhomogeneities; i.e., gross permeability variations. Viscosity differences lead to fingering of the less viscous solvent. This difference in viscosity accelerates the growth of fingers along paths previously developed because of permeability variations. Gravity differences lead to overriding by the usually less dense solvent. Al- though gravity effects are generally small at a flooding rate of 15 ft/D, they would, nevertheless, unnecessarily complicate the interpretation of flooding experiments. For this reason, all the experiments reported herein were done with matched density fluids. Fingering and the resultant poor areal sweep were recognized early as the dominant influences on the efficiency and the economics of miscible displacement processes. Much research effort has been spent on ways to minimize fingering and increase areal sweep such as the use of graded viscosity slugs4 or water slugs. Some researchers attempted to work out ways to prevent fingering completely; i.e., to achieve a stable displacement through gravity or latitudinal dispersion stabilization. l Others did not attempt to control fingering but obtained an economic process by merely recycling the solvent and sweeping pattern by pattern.6 During this period, all aspects of fingering came under close scrutiny. Some researchers reported on how fingering looks2 and how it is affected by viscosity ratio,7 geometry,7 and slug size.2,* Peaceman and Rachfordg suggested a mathematical approach to the prediction of unstable miscible displacements in relatively homogeneous sand packs, but their work cannot be extended conveniently to heterogeneous systems. Hence for a heterogeneous system, no method is presently available for predicting solvent cut and recovery as functions of pore volumes of solvent injected. The purpose of this investigation was to attempt to fill in the gap in our knowledge concerning the prediction of performance of unstable miscible dis placements. Necessarily, the system selected for study was a relatively simple one. The restrictions placed on the system were: (1) The system was linear; (2) The solvent was miscible in all proportions with the oil in place; (3) The solvent was continuously injected into the porous medium; (4) Gravitational effects were eliminated by matching densities and (5) All the flood rates were high and constant at 15 ft/D to avoid any small rate effect and to minimize any diffusion effects. To simplify and to indicate that both longitudinal dispersion and channeling arise from permeability variations, the effects which they cause or influence have been termed heterogeneity effects. The use of

By R. W. Olson, A. L. Benham

Viscous fingering was studied as it occurred in an open Hele-Shaw model (1 ft x 4 ft x 1/16 in.); it was also studied in the same model packed with 80-mesh glass beads during miscible displacements under unfavorable viscosity ratio conditions. It was determined that the lengths of 20 elements of the front were distributed normally around the average position of the combined elements under all conditions. Represented by the normal distribution, the length of the viscous fingers grew linearly with the distance traveled from the point of finger formation, increased with displacement velocity, and increased with increasing mobility ratio. Results obtained during the first few inches of displacement were of little or no use in predicting finger growth and/or finger length throughout the 4-ft model since the point of finger initiation cannot be predicted and would often occur several inches from the point of injection with either positive or apparent negative coordinates. The fingering occurring during a miscible slug displacement was much greater than would be predicted based upon the actual mobility ratios between in-place fluid and slug and between slug and following fluid, using the results described above for miscible displacement in the absence of a slug. Many of the experiments in the packed model showed that the rate of growth of the viscous fingers was diminishing toward the end of the displacement in the 4-ft long model, indicating that microscopic mixing, such as diffusion or dispersion, was decreasing the viscous fingering effect. INTRODUCTION Viscous fingering is a manifestation of a finger-shaped interface between displaced and displacing fluids occurring during typical miscible displacement projects for oil recovery. Its cause may be traced to the instability of a viscous fluid being displaced by a more mobile fluid. Viscous fingering takes on important significance in the miscible slug process 1.2 where it may be a dominant factor in determining minimum slug size. Experimental1-3 and theoretical4-7 studies of viscous fingering have been made by other investi- gators. However, an exhaustive study of the variables affecting viscous fingering had not been made. The present study was undertaken in an effort to determine the effects of some of the more obvious variables—such as mobility ratio, displacement velocity, distance displaced, and packing—upon viscous finger length and growth in a small laboratory model. It does not necessarily follow that conclusions reached by studying the results of this model study may be applied to the field. Future studies would have to evaluate the effect of model size on extension of these results. EQUIPMENT AND EXPERIMENTAL PROCEDURES A flow diagram for the equipment used in this study is shown in Fig. 1. Basically, it consisted of (1) a constant-rate pump for the injection of displacing fluid into (2) a model made up of two flat transparent plates spaced a small distance apart and (3) provisions for the production of fluid into a calibrated graduate. The displacing fluid contained a dye which allowed visual observation of finger spacing and length. The pump performed at rates from 2.5 cc/hr to 480 cc/hr, while pumping a light mineral oil into a bladder in a closed glass container holding the displacing fluid. The model was made up of two plates of 1-in. Plexiglass with dimensions 1-ft wide by 4-ft long. These were spaced 1/16 in. apart so that the internal model dimensions were 1-ft wide x 4-ft long x 1/16

By K. E. Starling

Phase equilibria data were used to develop an equation-of-state correlation for complex hydrocarbon mixtures, thereby circumventing difficulties associated with use of pressure-volume-temperature data. Equilibrium phase data for two condensate reservoir fluids were used to determine equation-of-state parameters for hydrocarbons as heavy as C22H46 in a modified form of the Benedict-Webb-Rubin equation of state. Comparative tests of K-values from the resultant correlation were made with data for condensate reservoir, separator and gas plant absorber mixtures. Generally, for temperatures above OF and computed liquid densities below 0.55 Ib-mole/cu ft, the modified BIVR equation predicted K-values in close agreement with the experimental data. INTRODUCTION Accurate predictions of thermodynamic behavior for complex hydrocarbon mixtures are necessary for many calculations in the petroleum industry. Because of the relatively high cost of extensive experimental data, many correlations for prediction of phase behavior have been developed. Some of these, such as the correlation of K-values presented in the NGAA Equilibrium Ratio Data Book, are totally empirical.l Others, such as the Benedict-Webb-Rubin (BWR) equation-of-state method, are semitheoretical.2-4 Unfortunately, published correlation methods often do not accurately predict the phase behavior of complex hydrocarbon mixtures, principally because of inadequate representation of the effects of components heavier than decane. This paper presents a new approach to this problem in which basic equation-of-state relations including heavy hydrocarbon effects are applied. Equilibrium phase data for two condensate reservoir fluids containing hydrocarbons as heavy as C22H46 are used to correlate component K-values. The BWR equation was chosen as the prototype equation of state for this study because of its proven capability for accurately predicting phase behavior and thermodynamic properties of light - hydrocarbon mixtures? Research was directed toward development of BWR parameters for the heavy hydrocarbons and modifications of the mathematical form of the BWR equation for application to complex hydrocarbon mixtures. The new equation-of-state approach presented differs considerably from previous methods in that phase equilibria rather than PVT data are used for determination of equation-of-state parameters. It is an explicit approach since the parameters are determined directly from mixture data. As such, it does not encounter problems inherent in the implicit method used by Benedict, Webb and Rubin 2-4 and numerous other investigators—one in which mixture parameters were postulated to be functions of the pure component parameters. The pure component BWR parameters, in turn, were determined from experimental PVT data. This method has been limited to mixtures containing components lighter than decane because of lack of vapor phase PVT data. Studies have been reported in which BWR parameters have been determined explicitly from PVT data for binary mixtures.6-8 However, since small concentrations of heavier components have only a minor effect on PVT behavior, it is doubtful that these explicit methods would yield useful results for condensate systems. Ellington and Eakin9 have shown that the accuracy of K-values predicted by an equation of state developed from mixture PVT data probably would be more than an order of magnitude lower than the accuracy of the PVT data. On the other hand, an equation of state utilizing parameters developed from phase equilibria data should predict K-values with accuracy comparable to the accuracy of the experimental phase compositions. This work applies this explicit approach with the objective of improving hydrocarbon mixture phase behavior predictions. Results are evaluated by comparing predicted and experimental K-values at reservoir, separator and gas plant

Jan 1, 1967

By B. A. Eaton, R. H. Jacoby

This paper presents the depletion-performance correlations developed using data from PVT studies and well-completion tests of 27 rich gas-condensate and volatile-oil reservoir fluids. The PVT behavior for each fluid war determined experimentally in the laboratory, and the depletion performance of each fluid was then calculated. Well-completion test data consisting of reservoir pressure and temperature, initial total gas-oil ratio, and stock-tank oil gravity were used with the depletion-performance data to develop general performance correlations. Three correlations are presented. The first equation shows the empirical relation of gas in place per barrel of hydrocarbon pore space with reservoir pressure and temperature, and initial total gas-oil ratio. The second equation is the correlation of stock-tank oil in place per barrel of hydrocarbon pore space. The third equation relates stock-tank oil production (to 500-psia depletion pressure) with pressure, temperature and oil gravity. All three correlations are also presented in convenient graphical form. Gas production to 500-psia depletion pressure is approximately 92.6 per cent of the gas in place for each fluid, so no correlation was developed for gas production. INTRODUCTION The rate of discovery of shallow, low-pressure and low-temperature oil reservoirs in the United States has been declining in recent years. For this reason, the search for oil has gone to greater depths where high pressures and temperatures are encountered. A high-pressure and high-temperature reservoir is the environment in which many rich gas condensates and highly volatile oils are found. Because of this trend, the performance prediction of such reservoirs has become an important task of the reservoir engineer. The need and usefulness of a performance correlation based on field data that are available from the usual well-completion tests have previously been established. These earlier attempts have shown that the phase behavior of rich gas-condensate systems and volatile-oil systems can be correlated with fluid properties, pressure and temperature. However, very little information on gas-condensate depletion performance5 was available at the time of these earlier attempts. The availability of more recent data has made possible the development of an improved correlation. The data, as listed in Tables I and 2, are for 25 gas-condensate and two volatile-oil systems that were studied. The problem was to correlate the data for in-place fluids and the production data with the initial total gas-oil ratio, initial reservoir pressure and temperature, and stock-tank oil gravity. These correlations were obtained by multiple regression analysis, whereby the data were fitted to an equation by the method of least squared error. BASIS OF THE CORRELATIONS 1. The data used in the development of these correlations were obtained from depletion-performance calculations6 and well-test or production measurements. 2. All of the independent variables used in the correlations are available from the usual well-test completion data. 3. The reservoir volume considered is 1 bbl of hydrocarbon pore space (HCPS). 4. The gas-oil ratio used is the total standard cubic feet of gas (at 14.7 psia, 60F) coming off the entire separator train per barrel of stock-tank oil. This minimizes the effect of separator conditions which might vary from system to system. 5. The correlated pressures are the actual reservoir pressures which might or might not be the saturation pressure of the reservoir fluid.

Jan 1, 1966

By H. O. Jahns

This paper describes the application of regression analysis for obtaining a two-dimensional areal description of heterogeneous reservoirs from short-term pressure-time data such as that obtained in interference tests. The method replaces the time-consuming trial-and-error procedure commonly used to match field data on an electric analyzer or digital computer with a systematic search which is programmed for a computer. The computer program adjusts the properties of a reservoir model automatically until a least-squares fit is obtained between observed and calculated pressure data. The reservoir is simulated by a single-phase, compressible, two-dimensional model. It is divided into a number of homogeneous blocks whose transmissibility (kh/F) and storage (Fch) values are varied to obtain the least-squares fit. The reliability of these values is determined from their standard deviations and correlation coefficients. Although the method is rigorously applicable to single-phase flow only, multiphase flow can be handled provided saturation changes are small during the test. Possibly the method can also be used to obtain a reservoir description from pressure-production history, but this application is outside the scope of this work. The paper includes, in addition to a description of the numerical procedure, a discussion of some of the problems associated with the method. Rules are given to help in selecting the number of homogeneous blocks and deciding upon their arrangement. The uniqueness of a reservoir description is considered. Finally, the use of the method is illustrated by the interpretation of field data from two interference tests. INTRODUCTION Pressure data from short-term transient tests, such as single-well and interference tests, are widely used to obtain reservoir properties. These tests are usually analyzed by assuming a simple reservoir model; very often, a homogeneous one is used. As a result, analysis of the transient data from each well frequently gives different values for reservoir properties. The problem then arises to combine all these differing results into a more detailed picture of the reservoir. One technique is to simulate the reservoir with a digital computer or with an electrical analyzer and to adjust the reservoir parameters by trial and error until the simulated pressure data are in reasonable agreement with the observed pressure data for all wells. Although this method has been used for both transient tests and pressure-history data, it is time-consuming and subjective. A second technique uses regression analysis to replace the trial-and-error procedure with a systematic search that can be programmed for a digital computer. Use of regression analysis in reservoir description was proposed recently by Jacquard and Jain.1 They divided the reservoir into a number of homogeneous blocks whose properties are varied until a least-squares fit is obtained between observed and calculated pressures. However, they did not consider their technique to be operational, mainly because of ".. .(I) the lack of experience in using the method .. . notably for the improvement of convergence; and (2) limitations imposed by the insufficiency of available computers". 1 While the analysis presented in this paper applies the same general principle used by Jacquard and Jain, the specific method is significantly different. Some differences are (1) the regression problem is solved in a different way which requires less computer time in most cases; (2) a stepwise solution, in which the detail in the reservoir description is increased from step to step, is used to improve convergence; and (3) the reliability of the estimated reservoir .properties.. as measured by their standard deviation and correlation coefficient, is estimated'

Jan 1, 1967

By T. K. Perkins, O. C. Johnston

Because of the influence of dispersion on miscible-displacement processes, diffusion and dispersion phenomena in parous rocks are of current interest in the oil industry. This paper reviews and summarizes a great deal of pertinent information from the literature. Porous media (both unconsolidated packs and consolidated rocks) can be visualized as a network of flow chambers, having random size and flow conductivity, connected together by openings of smaller size. In such a porous medium, the apparent diffusion coefficient D is less than the molecular diffusion coefficient Do, as measured in the absence of a porous medium. For packs of unconsolidated granular material the ratio D/D, is about 0.6 to 0.7. For all porous rocks, both cemented and unconsolidated, the ratio of diffusion coefficients can also be represented as D/Do = 1/FF where F is the for- mation electrical resistivity' factor and 4 is the porosity. If fluids are flowing through the porous medium, dispersion may be greater than that due to diffusion alone. At moderate flow rates the porous medium will create a slightly asymmetrical mix zone (trailing edge stretched out), with the longitudinal dispersion coefficient approximately proportional to the first power of average fluid velocity (if composition is nearly equalized in pore spaces by diffusion). If the velocity in interstices is large enough, there will be insufficient time for diffusion to equalize concentration within pore spaces. In this region, longitudinal dispersion increases more rapidly than fluid velocity. At low velocities in interstices, transverse dispersion is characterized by a region in which transverse diffusion dominates. If the fluid velocity gets high enough, there will be a transition into a region where there is stream splitting with mass transfer but with insufficient residence time to completely damp-out concentration variations within pore spaces. There are several variables that must be con- trolled to get consistent longitudinal and transverse dispersion results, viz., (1) edge effect in packed tubes, (2) particle size distribution, (3) particle shape, (4) packing or permeability heterogeneities, (5) viscosity ratios, (6) gravity forces, (7) amount of turbulence, and (8) effect of an immobile phase. INTRODUCTION Diffusion and dispersion in porous rocks are of current interest to the oil industry. This interest arises because of the influence of dispersion on miscible-displacement processes. In a recovery process utilizing a zone of miscible fluid, there is the possibility of losing miscibility by dissipating the miscible fluid or by channeling or "fingering" through the miscible zone. Diffusion and dispersion are two of the mechanisms that may lead to mixing and dissipation of the slug. On the other hand, dispersion may tend to damp-out viscous fingers which may be channeling through the miscible slug. 58 Hence, dispersion may be detrimental or beneficial (if it prevents fingering through the miscible zone). Therefore, it is doubly important that we understand these processes. In this paper we review, summarize and interpret a great deal of information from the literature. In particular, we will briefly discuss molecular diffusion in miscible fluids. Then we will discuss what differences to expect for diffusion in a porous rock. If there is movement of the fluid through the rock, then there may be an additional mixing or "dispersion". Furthermore, the dispersion longitudinally (in the direction of gross fluid movement) and transverse to the direction of fluid movement will not be equal. We will discuss both types of dispersion as well as several variables which can affect dispersion (viscosity differences, density differences, turbulence, heterogeneity of media, etc.). This group of variables has sometimes led to difficulty when comparing literature data. DIFFUSION OF MISCIBLE FLUIDS If two miscible fluids are in contact, with an initially sharp interface, they will slowly diffuse into one another. As time passes, the sharp interface between the two fluids will become a diffuse

By G. W. Thomas

A mathematical tnodel of forward combustion in an oil reservoir is treated in this paper. The model describes a radial system having a vertical section of essentially infinite thickness, all of which is permeable to gas flow. Combustion, however, is presumed initiated over a limited thickness of the total vertical section. In the interval supporting cotnbustion, the mechanisms of radial conduction, convection and heat generation are taken into account. Above and below the burning interval, heat transport in the radial direction is by cottduction and convection. Vertical heat losses from the ignited interval are accounted for by conduction alone. A general solution is presented for the temperature distribution caused by radial movement of the combustion front. The results show that no feedback of heat occurs into the ignited interval when convection and conduction are acting in the bounding media. Peak temperatures are also 5 to 10 per cent higher than in the case where heat transport in the bounding media is by conduction alone. We arbitrarily define vertical coverage to be that fraction of the total ignited interval which is at 600F above atnbient, or greater, at any given time. The radial distance at which the vertical coverage becomes zero is the propagation range of the combustion front. It was found that an increase in vertical coverage results when the oxygen concentration, fuel concentration or gas-injection rate is increased. Moreover, the combustion front can be propagated 10 to 15 per cent further than in the case where only conduction is acting above and below the ignited interval. INTRODUCTION In the theoretical treatment of forward combustion in a radial system, one of the problems encountered is the determination of the transient temperature distributions caused by an expanding cylindrical heat source. Bailey and Larkin' and Ramey' simultaneously presented analytical solutions to the problem assuming heat transport by conduction alone. In a subsequent publication, Bailey and Larkin3 included the effects of both conduction and convection while treating linear and radial models. In this latter work, however, vertical heat losses were largely neglected. Selig and Couch' dealt with a radial model in which both conduction and convection were acting. Only a limiting case involving vertical heat losses was considered, however. Namely, temperatures on the boundary of the bed of interest were set equal to zero. Solutions thus obtained were representative of a system having a maximum vertical heat flux. Chu5 recently treated a more general case in which a permeable bed was considered bounded by impermeable media. Conduction and convection took place within the bed, and only conduction outside of the bed. The effects of vertical heat losses were included in his study. Solutions were obtained by numerical techniques. This paper is an extension of the theoretical work of other authors pertaining to forward combustion in a radial system. In particular, a mathematical model of the process is treated in which heat generation occurs over a small vertical interval of a larger permeable section. In the interval supporting heat generation, and above and below this interval, the mechanisms of radial conduction and convection are also presumed acting. Heat losses from the ignited interval are accounted for by vertical conduction. An analytical solution for the temperature distribution caused by radial movement of the burning front is presented. The effects of certain process variables are indicated and comparisons with Chu's results are made. THEORY To render the mechanism of forward combustion tractable to mathematical treatment, we idealize the problem to the extent of assuming continuous reservoir media possessing homogeneous and isotropic properties. The following additional assumptions are implicit in this analysis. 1. The thermal parameters, i.e., heat capacities, thermal conductivities and thermal diffusivities are invariant with temperature and pressure. Moreover, the bounding media possess the same thermal properties as the bed of interest. 2. The temperatures of the porous media and its contained fluids at any point and at any time are equal. 3. The reaction rate between the oxidant gas and the fuel is infinite. This assumption implies that the incoming oxygen concentration instantaneously goes to zero within an infinitesimal distance, i.e., the width of the combustion zone is negligible. 4. The rate of gas injection is constant and corresponds to the average rate throughout the lifetime of the project. 5. The fuel concentration is constant throughout the volume of rock swept out by the burning zone. 6. There is complete burnoff of fuel. This assumption demands that the rate of propagation of the burning front equals the rate of fuel burnoff. In a radial system, with a

By R. F. Nielsen, D. E. Menzie

The object of this study was to determine if crude oil could be produced successfully by a process of crude oil vaporization using carbon dioxide repressuring. This process appears to have application to highly fractured formations where the major oil content of the reservoir is contained in the non-fractured porosity with little associated permeability. Crude oil was introduced into the windowed cell and carbon dioxide was charged to the cell at the desired pressure. A vapor space was formed above the oil, and the crude oil-carbon dioxide mixture was allowed to come to equilibrium. The vapor phase was removed and the vaporized oil collected as condensate. Samples of all produced and unproduced fluids were analyzed. Tests were also performed to evaluate the amount of vaporized oil that can he produced by rocking from a high to a lower pressure. The carbon dioxide repressuring process was applied to a sand-filled cell to investigate the performance in a porous medium. A test was performed to evaluate how the condensate recovery changes as the size of the gas cap in contact with the oil changes. INTRODUCTION This study has been directed toward a relatively new process of vaporization of crude oil designed to increase ultimate production of hydrocarbons through the application of carbon dioxide to an oil reservoir. Suggested advantages of carbon dioxide repressuring of a petroleum reservoir are: (1) reduction in viscosity of liquid hydrocarbons due to the solubility of carbon dioxide in crude oil, (2) swelling of the reservoir oil into a larger liquid-oil volume with a resulting increase in production and decrease in residual oil saturation due to an increase in the relative permeability to oil, (3) displacement of more stock-tank oil from the reservoir since the residual liquid is a swelled crude oil, and (4) gasification of some of the hydrocarbons into a carbon dioxide-hydrocarbon vapor mixture. Balanced against these advantages are several detrimental factors which must be evaluated; i.e., high compression costs and corrosion of well equipment and flow lines. Some of the more outstanding contributions to the study of carbon dioxide injection have been reviewed in order to furnish a basis for a continuation of research pertaining to this method. The literature reviewed1-8 has been limited to that dealing with carbon dioxide repressuring processes or with carbon dioxide-crude oil-natural gas phase behavior. Articles relating to carbonated water injection and literature published on the use of low pressure carbon dioxide gas injection in water flooding have not been included in this study. In 1941 Pirson5 suggested the high pressure injection of carbon dioxide into a partially depleted reservoir for the purpose of causing the reservoir oil to vaporize and thus produce the oil as a vapor along with the carbon dioxide gas. By reducing the pressure on this produced mixture of hydrocarbons and carbon dioxide at the surface, it was proposed to separate the hydrocarbons from the carrier gas. He theorized that essentially all the oil in a reservoir could be produced by simply injecting enough carbon dioxide to vaporize the residual oil. This present investigation deals with the vaporization of a crude oil by carbon dioxide, the molecular weight and gravity of the vaporized oil product and the characteristics of the residual oil after several repressuring cycles with carbon dioxide. An attempt is made to evaluate the merits of a vaporization process for the crude oil rather than a flow process where the oil recovery is determined by relative permeability considerations. Such a vaporization of crude oil by carbon dioxide repressuring appears to have possible use in a highly fractured formation where the major oil content of the reservoir is contained in the non-fractured porosity with little permeability. The carbon dioxide flows into the fractures, contacts the crude oil in the matrix and vaporizes part of the crude oil; this vaporized oil is produced and recovered and the carbon dioxide is reinjected again. The specific problem of this study is to attempt to answer this question; Can crude oil be produced successfully (technically, but without economic considerations) from a petroleum reservoir by a process of vaporization of the crude oil by carbon dioxide repressuring? DEFINITION OF TERMS AS APPLIED IN THIS STUDY Carbon Dioxide Contact: One cycle in which carbon dioxide was injected and bled off. Condensate: The hydrocarbon liquid which was condensed out of the mixture of hydrocarbon-carbon dioxide vapor upon reduction of the pressure of the vapor. Hydrocarbons Produced (HCP): All the hydrocarbon!, which were vaporized by the carbon dioxide repressuring process and were removed from the cell during any specific cycle or carbon dioxide contact. Hydrocarbons Unproduced (HCU): All the hydrocarbons which were not vaporized by the carbon dioxide

By R. L. Perrine

This paper presents results obtained by consider ation of unstable miscible displacement within a porous medium as a quasi-turbulent displacement process, that is, a process characterized by fluctuations around a simpler, mean flow behavior. The instability usually results from the miscible displacing fluid having a viscosity lower than that of the displaced fluid. Such instability is evidenced in suitable experiments by the presence of viscous fingering. A unified theoretical relationship applicable to both stable and unstable displacement is derived. Parameters in the theory are obtained by empirical means. Although the theory describes an idealized conceptual model, agreement with experiment is quite reasonable up to moderately adverse viscosity ratios, such as 17:1. INTRODUCTION Many experimental studies of miscible displacement in porous media have been conducted with prediction of reservoir displacement efficiency as the ultimate objective. Most such studies have utilized lower displacing than displaced fluid viscosity, scaled to potential reservoir fluid pairs. Theoretical approaches have been largely limited to unit viscosity ratios, however, in spite of the necessity for an understanding of the mechanism of the unstable, adverse viscosity ratio displacement process. An obvious reason is the difficulty of describing in mathematical form the viscous fingering characteristic of these conditions. Observation of experiments conducted with dyed fluids in transparent systems suggests an analogy between unstable miscible displacement and turbulent flow in a pipe. In both cases there are fluctuations around a simpler, mean flow behavior. An important difference is that flow behavior of interest in the porous medium is entirely transient, contained within a transition zone between displacing and displaced fluids. Transient behavior complicates description in that coefficients in the equations become variables rather than constants. In the study reported here, the analogy with turbulent flow has been used in treating unstable miscible displacement as a quasi-turbulent displacement process. The purpose has been to derive, even if restricted to an idealized conceptual model, a unified theoretical relationship applicable to both stable and unstable displacement. A relationship meeting these specifications up to moderately adverse viscosity ratios, such as 17:1, is presented. One fluctuation parameter in the theory and dispersion coefficients are obtained by empirical means. The idealized theory is compared with experimental results. UNSTABLE MISCIBLE DISPLACEMENT AS QUASI-TURBULENT DARCY FLOW The miscible displacement of one fluid by another within a porous medium is considered in this study. Flow conditions are such that Darcy's law is applicable. It is further assumed that, by the stability criterion proposed by perrine,1 initial flow conditions lie well within the regime of instability. Thus substantial viscous fingering is assured. We wish to show how this flow regime can be represented as quasi-turbulent. That is, the Reynolds number for the established flow conditions is below that for inertial or turbulent flow within a porous medium, and lies within the regime for which Darcy's law is valid. Yet flow behavior can be described as the combination of some relatively simple average result, and characteristic fluctuations that are superimposed on the simpler behavior. Stated another way, flow behavior includes the movement of layers of fluids with different velocities past or over one another. Such descriptions are characteristic of turbulent pipe flow. The fluctuations in the present case are a direct consequence of the local viscous fingering which accompanies the unstable displacement process. Should displacement become stable, fluctuations would die out. It is of particular importance to note that fluctuations such as these can interact in a way that contributes to material transport by the basic flow. The basic transport equations required for engineering calculations must be modified to reflect this fact.

By W. B. Braden

This paper presents a suitable method for predicting gas-free oil viscosities at temperatures up to 500F knowing only the API gravity of the oil at 60F and the viscosity of the oil measured at any relatively low temperature. The API pravity and the one viscosity value are used as parameters to determine the slope of a straight line on the ASTM slanaord viscosity-temperature chart. Then, knowing the slope of the line and one point on the line, the vrscosities at higher temperatures can be determined. The line slope correlations were developed at I00 and 210F since viscosity data are frequently measured at these temperatures. A procedure is given for predicting line slopes from measurements at other tetnperatures. A nomogram is furnished for solving the relationship. The correlation has been evaluated at temperatures up to 5OOF for oils varyzng in gravity from 10 to 33 " API. The correiution is applicable only to Newtonian fluids. Comparison at 500F of true viscosities and those predicted from values at 100F shows an average deviation of 3.0 per cent (maximum deviation of 6.0 per cent). Predictions from the values at 21 0F for the same oils how an average deviation of 1.5 per cent (maximum deviation of 3.4 per cent). INTRODUCTION Correlations have been developed by Beal' and by Chew and Connally' for predicting viscosities of gas-saturated oils at reservoir conditions. Each of these correlations requires a knowledge of the solution gas-oil ratio and the viscosity of the gas-free oil at the reservoir temperature. For temperatures below 350F, measurements of the gas-free oil viscosities can be made easily using commercially available equipment. In thermal recovery processes, however, reservoir temperatures well in excess of 350F are encountered. Viscosity measurements at such conditions are more difficult and time consuming and require modification of existing equipment or the construction of new equipment. Measurements are further complicated by the difficulty of handling highly viscous oils associated with thermal recovery processes. Therefore, it is desirable to have a correlation which allows accurate prediction of viscosities at high temperatures. A commonly used technique for predicting viscosities at high temperatures is to measure the viscosities at two lower temperatures, plot the values on ASTM standard viscosity-temperature charts and extrapolate to the temperatures desired. If either of the values is slightly in error, the extrapolated value can be significantly in error. To justify an extrapolation, three points are actually necessary. This procedure can consume much time, particularly with heavy oils. Considering the cost of viscosity measurements, it would be desirable to eliminate the need for direct measurements by having correlations which would allow viscosity predictions from other physical or chemical properties. Beal1 investigated the possibility of correlating viscosity with oil gravity at temperatures from 100 to 220F. While showing that a general relationship exists, he also found significant deviations. It is possible that correlations will be developed based on oil composition as more information becomes available. While not eliminating the need for viscosity rneasurements, the method presented herein requires that only one viscosity measurement be made. The API gravity must also be known. The theory is based on the fact that the viscosity of paraffins (high gravity) changes less with temperature than does the viscosity of naph-thenes or aromatics (low gravity). The gravity. therefore, is used as a parameter to determine the slope of a straight line on the ASTM standard viscosity-temperature charts. The correlation is applicable only to Newtonian oils, and deviations due to thermal decomposition and nonhomo-geneity cannot be predicted. Oils containing additives have not been evaluated. PROCEDURE Fifteen oils were used in developing the correlation; eight were crudes and seven were processed oils. Oil gravities varied from 9.9" API (naphthene base) to 32.7' API (paraffin base). The temperature range studied was 81 to 516F. Each oil used had a minimum of three viscosity measurements and each plotted essentially as a straight line on the ASTM charts. In all, 91 viscosity measurements were used in the correlation. Saybolt, rolling ball and capillary tube viscometers were used for the measurements. Viscosity data for Samples 1, 2, 4, 7, 10, 11 and 14 were obtained in Texaco, Inc. laboratories. The data for Samples 3, 5, 6, 8, 9, 12 and 15 were from Fortsch and Wilson,3 and data for Sample 13 were from Dean and Lane.' All data points used in the correlation are plotted in Fig. 1. It is seen that some of the viscosity values deviated slightly from the straight-line plots at the higher temperatures. Properties of the oils after exposure to the

Jan 1, 1967

By P. M. Blair, D. W. Peaceman

The shape and position of the gas-oil transition zone during downdip displacement of oil by gas has been calculated using flow equations which include the effects of gravity, relative permeability, capillary pressure and compressibility of the fluids. The calculations treat the problem in two space dimensions, and results are compared with data from a laboratory model tilted at 30° and 60° from the horizontal on displacements near and above the maximum rate at which gravity segregation prevents channeling of the gas along the top of the stratum. The good agreement between calculated and experimental results demonstrates the validity of the tecbnique as well as that of the flow equations. INTRODUCTION Knowledge of the fluid distribution and movement in an oil reservoir is important in producing operations and estimation of reserves. The history of the oil industry has included steady progress in improving the accuracy of calculations which provide the required knowledge. The earliest method of calculating reservoir performance consisted of material-balance equations based on the assumption that all properties were uniform throughout a reservoir.l For many reservoirs such a simple formulation is still the most useful. However, when large pressure and saturation gradients exist in a reservoir, the assumption of uniform values throughout may lead to significant error. To reduce these errors, Buckley and Leverett2 introduced a displacement eqation which considers pressure and saturation gradients. Methods available at that time permitted solutions to the Buckley-Leverett equation in one space dimension; these solutions have been very useful in solving many problems related to the production of oil. However, the one-dimensional methods are not adequate for systems in which saturations vary in directions other than the direction of flow. An example of such a system is the case of gas displacing oil down a dipping stratum in which the gas-oil contact becomes significantly tilted. Of course, the Buckley-Leverett displacement method cannot predict the tilt of the gas-oil contact. Recent improvements3,4 of the one-dimensional Buckley-Leverett method achieve some success in predicting the tilt of the gas-oil contact at sufficiently low flow rates. However, at rates high enough that the viscous pressye gradient nearly equals or exceeds the gravity gradient, even these improved one-dimen,sional methods incorrectly predict the shape and velocity of the contact. Further progress in estimating such fluid movements in a reservoir appears to require consideration of the problem in more than one space dimension. The recent two-dimensional method of Douglas, Peaceman and Rachford5 appears adaptable to calculate changes with time of the saturation distribution in a vertical cross-section of a reservoir. The movement of saturation contours should represent the moving fluid contacts and include the effects of crossflow due to gravity, as well as variations in the rock and fluid properties. The nonlinear nature of the equations used in the method has prevented proof of the validity of the solutions. Douglas, Peaceman and Rachford made some comparisons with experiment but did not include cases in which gravity was important nor cases involving displacement by the nonwetting phase.For these reasons, a test of the two-dimensional method for a case in which these factors are included would be very desirable. The test selected was a comparison of calculated results with those from a carefully controlled laboratory experiment on a model with measured physical properties. The model selected was one in which gas displaced oil down a tilted, rectangular sand pack. The model can be thought of as representing a vertical cross-section taken parallel to the dip of a reservoir. The displacement thus simulates gas displacing oil downdip that might result from gas-cap expansion or gas injection. This report outlines the calculation method used

By D. M. James, J. C. Martin

The results are presented of a study of the application of analytical methods to the solution of two-phase flow into single wells. Approximate analytical expressions for the pressure distribution in two-phase flow are found ;;; a number of conditions. The results obtained from the analytical solutions me found to be in good agreement with results obtained by finite difference techniques using a high speed digital computer. Mathematical solutions for four sets of boundary conditions are presented. All of these solutions are composed of a short-term transient plus a steady or quasi-steady state. The rates of decay of the short-lived transients are analyzed. It is found that the durations of the short-term transients may be characterized by a parameter defined as the time constant which can be determined from simple relations. It is shown also that if the outer radius is much greater than the radius of the well. the short term transients decay at rates which are proportional to the square of the exterior radius, and the rates of decay are only slightly dependent upon the radius ratio. Numerical solutions based on finite difference techniques are presented for a number of conditions. The numerical solutions are in good agreement with the predictions based on the theoretical analysis for small and moderate drawdowns. Examples involving large drawdowns indicate that the nonlinearities in the equations of flow do not appreciably alter the longevity of the short-term transients. In all cases the time required for the short-tern transients to disappear is predicted satisfactorily. INTRODUCTION The mechanism by which oil and gas flow into a single well is of vital interest to the petroleum industry. The fundamental equations of two-phase flow which describe this mechanism are nonlinear partial differential equations. Numerical solutions of these equations describing pressure transients have been obtained with the aid of electronic computers. Although solutions obtained in this manner take into account a large number of effects, the reduction of this information to useful generalities is difficult. One method of obtaining generalities is the use of linearized approximations of the nonlinear equations. Since it is possible to obtain explicit solutions of the linearized equation, general properties of the role of pressure in the flow mechanism may be ascertained. Results obtained from this approach are limited to some extent by the linearizing assumptions. The severity of these limitations may be evaluated by comparing solutions of the linear equation with numerical solutions of the more exact nonlinear equations of two-phase flow. In the past considerable amount of work has been devoted to studying pressure build-up using the single-phase flow theory. Unfortunate1y, most prep sure build-up tests involve multiphase flow. A small amount of work has been done studying pressure build-up where the flow is two-phase. The encouraging results of these studies suggest that useful results may be found from additional studies of not only pressure build-up, but also the rapid transients associated with placing a well on production. This paper presents the basic theory of the pressure transients associated with placing a well on production and with closing it in. The paper is concerned chiefly with two-phase compressible flow; however, the results also apply to single-phase flow. The results are based on analytical solutions of the flow equations, and they are verified by numerical solutions using finite difference techniques. Much of the previous work on compressible flow into wells has been confined to single-phase flow. Some work has been done on compressible two-phase steady state flow, and solutions of the equations of flow have been found by finite difference techniques using high-speed computers. Muskat1 presents some rather general solutions to the equations of single-phase compressible flow into wells. Much work has been done on pressure build-up in wells (see, for example, Refs. 2-6). Almost all of the work on pressure build-up concerns single-phase flow with the exception of Ref. 5 and part of Ref. 2. Some work has been done on pressure fall-off in injection wells.718 Muskat9 presents the solution of the equations for radial steady state two-phase compressible flow. Handyl0 extends Muskat's solution to obtain

By L. T. Stanley

A method is described whereby the behavior of the front is approximated when gas is injected into a thick reservoir having reasonably homogeneous properties. The method is applied to the Abqaiq field of Saudi Arabia. The producing member is the Arab-D, a high-permeability calcarenite which has no continuous betlding planes. Field performance indicates the injected gas to be partially miscible with the reservoir oil because gas volume factors experienced are lower than those predicted by the gas law. Eqlrilibril~m calculations confirnz ihis and also show that the volume factor may vary appreciably, depending on the relative quantities of oil and gas in contact in the reservoir; consequently, volutnetric balances are made for the purpose of determining the behavior of the gas volume factor during injection history. The effect of gravity segregation on position and shape of the gas front is approximated by dividing the producing section into three zones vertically and tracing the movement of the average interface in each zone, with counterflow of oil and gas taking place betweell zones behind the front. This results in a gas front that advances more rapidly along the top of [he formation than along the base, producing an "umbrella" effect in cross section that becomes more pronounced as the front progresses. The calculations cover the gas-injection period of Abqaiq field history through 1958, and the position of the gas front is plotted at one-year intervals. Finally, some comparisotls between the calcrlated fronial position and the gas-oil contact as Measured by neutron logs are made. INTRODUCTION Fig. 1 is a structure contour map of the Abqaiq field. For analytical purpcses, the field has been divided into two areas designated A and B, as shown on the map. Gas injection into Area A was begun in early 1954 into two wells located near the crest of the structure on the north-south axis of the anticline. During the period covered by this paper, the average injection rate has been 122 MMscf/D with a maximum sustained rate of 160 MMscf/D. The reservoir pressure has been maintained at 2,450 psig since start of gas injection. The reservoir originally was undersaturated but, prior to pressure maintenance, the pressure was drawn below the bubble point in most of the Area A section, creating an in-place free-gas saturation of 8 per cent at the crest of the structure. The average critical gas saturation is 13 per cent established by laboratcry core tests. That the entire reservoir remained below critical gas saturation prior to gas injection is borne out by the fact that none of the crestal wells have exhibited high gas-oil ratios during production history. The Area A reservoir is a carbonate section about 200-ft thick having average porosity and permeability of the order of 22 per cent and 500 md, respectively. The productive horizon is the Arab-D member of Jurassic age and is encountered at an average depth of

By W. E. Showalter

This paper discusses some of the results of combstion-drive tests which were made in a test cell using a sand bed 10 in. in diameter x 10-ft long. The test method is illustrated and described. The relationship between the API gravity of the in situ oil and the amount of air required for combustion drive is discussed in detail. Other things constant, the air requirement for combustion drive increases as the APl gravity of the in situ oil decreases. If the test results apply to actual reservoirs, the lowest-priced oils may cost the most to recover by this method. Information is shown which indicates that the effect of pressure on the amount of hydrocarbon burned is not large. A method of predicting air requirements from the API gravity of the in situ oil is presented. INTRODUCTION Combustion drive is the term used to identify the process of interstitial or in situ burning as an oil recovery method. Part of the in situ oil is burned to generate the energy needed to produce the remainder of the oil. Combustion drive as an oil recovery mechanism remains an economic uncertainty in spite of all the work that has been done by the industry in both laboratory and field?-6 This paper will show some of the results of tests which were made in a test cell for the purpose of studying the nature of the combustion-drive process. It will present data which indicate that the API gravity of the in situ oil is a significant indicator of the amount of air required to drive a burning front through oil sand. Air requirement varies inversely as the API gravity of the in situ oil. EXPERIMENTAL The tests were performed in a cell which utilized a cylindrical sand section 10 in. in diameter x 10-ft long. The thin-walled metal pipe which held the sand was wound with twenty external electrical resistance heaters which, by means of an automatic controller, maintained adjacent sections of the wall of the pipe at temperatures equal to the temperatures of the contained sand. Each heater covered 6 linear in. of the pipe. By this means lateral heat loss from the sand section was minimized, thereby causing the sand section to simulate more closely a horizontal increment of a combustion-drive reservoir. Fig. 1 shows a schematic diagram of the test assembly. Thermocouples to measure the temperature in the sand were located every 6 in. along the length of the sand section. The pipe containing the sand was enclosed in a cell designed for an operating pressure of 500 psig. The inlet air pressure was controlled at the inlet, and the gas flow rate was controlled and measured at the outlet of the cell. The oil sand used for the tests was prepared by mixing first water and then oil with the noncon-solidated sand using a closed mixer similar to a cement mixer. Table 1 shows a screen analysis of the sand. Ninety percent of the sand was 100 mesh or finer. This sand was a mixture of 80 per cent No. 120 Nevada White Sand and 20 per cent Tennessee Hi-Fusion Moulding Sand No. 3. The Nevada sand was a clean silica sand. The Tennessee sand con-

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967