Search Documents

By C. Dirksen

This study investigates whether oxygen consumption during segregated forward combustion may be affected by natural convection. Linear theory indicates that thermal instability occurs in a horizontal porous medium when the modified Rayleigh number N 'Ra exceeds a critical value of about 40. Most experimental results, including those described here, indicate this value to be about an order of magnitude smaller. N ', is derived from "differential similarity" and the proper time scale factor is also obtained. The critical value of N'ra for a horizontal, water-saturated porous medium at atmospheric pressure subjected to a uertical temperature gradient was found to be about 3.0. The ambiguity in the value of N'R, arising when fluid properties cannot be assumed constant is indicated. Natural convection was observed for oblique systems below N 'Ra , while simultaneous horizontal flow did not affect N Racrit. The expected range of values for the various parameters pertaining to segregated forward combustion is indicated and the corresponding values of N 'Ra calculated. Only under very favorable conditions, such as high permeability, high pressure and low top temperature, can the critical value of N 'R, be exceeded. Thickness of the convection layer should not become too large, since the time required to obtain significant mixing is proportional to the square of the thickness. For an oblique combustion frant the conditions for significant mixing may be less strenuous. LNTRODUCTION Laboratory1,2 and field314 experiments have shown that a markedly upgraded oil can be produced by in situ forward combustion. Under certain conditions injected air may move through a high permeability, gas-saturated channel over the oil-saturated layer, while combustion takes place, exclusively, at the horizontal interface (Fig. 1). Such a channel may be a natural gas cap, or it may result from a forward combustion process in which injected air and combustion gases fingered through the top of the formation, causing a premature breakthrough at the production well. In this particular case of segregated forward combustion, the air stream moves parallel to the combustion front. Molecular diffusion and transverse mechanical dispersion will be insufficient to cause all the oxygen in the air stream to reach the combustion front. When the oxygen in the air stream is not completely consumed in the combustion process, it will reach the production well and create a serious fire hazard. In addition, the reduced efficiency may make the process economically unprofitable while also serious corrosion problems may be experienced in the production we11.3 This undesirable situation should, therefore, be prevented. However, in a situation as depicted in Fig. 1, an unstable temperature gradient exists in the gas phase. The combustion is taking place at the lower boundary at a nearly constant temperature that may be anywhere from 600 to 1,200F, depending on the conditions. The upper boundary of the gas phase is overlying, impermeable rock; its temperature might be around 100F, but may increase to 300F or higher due to conduction and convection as the combustion process continues. This represents a steep temperature gradient and it is conceivable that under those conditions free or natural convection will take place in the gas

Jan 1, 1967

By I. H. Silberberg, R. L. Reed, O. K. Kimbler

interfacial films have frequently been observed at interfaces between certain crude oils and water. Several investigators have postulated that the presence of these films should influence the efficiency of oil recovery in water drive or waterflood operations. They may also influence the stability of emulsions which are sometimes a problem in petroleum production, and may be a factor in the fonation of paraffin deposits in oil well tubing and flow lines. This paper presents a technique with which a modified Langmuir film balance may be used to study the compressibility and collapse pressure of these natural interfacial films. Experimental data are presented for several nude oil-water systems. Data developed are used to infer the phase state of the film as a function of such variables as rate of reduction of interfacial area, ionic composition of the subtrate and pH of the subtrate. A film of known physical characteristics is shown to have a sinnificant effect on oil recovery from an unconsolidated sand pack. Possible applications of these results to petroleum production are discussed. INTRODUCTION The use of water to displace petroleum from reservoir rocks is of major importance both as a primary and a secondary recovery process. As water invades the rock, oil is completely displaced from some pores and left as a discontinuous phase in other pores. The manner in which water moves from pore to pore is strongly influenced by capillary forces. In view of the complexity of reservoir fluid systems, there can be little doubt that complicated interactions take place at both the liquid-solid and oil-water interfaces. One of the more interesting, and least understood, of the phenomena which take place at the oil-water interface is the formation of interfacial films. These films are believed to result from the adsorption of high molecular weight polar molecules at the interface.l.2 Presence of such molecules may cause a striking alteration in interfacial tension. When the oil-water interfacial area of certain crudes is rapidly reduced, a thin region (film) about the interface assumes the appearance of a solid membrane, and striations, wrinkles and gross distortions may occur. If such a film is solid, it should greatly alter the interfacial tension normally assumed to exist between the oil and water phases. If the membrane is continuous, a solid phase would separate the oil and water. Interfacial films between crude oil and water were observed in 1949 by Bartell and Niederhausers who commented upon the apparent rigidity of the films and their possible importance in the petroleum industry. Morrell and Egloff4 had earlier attributed the extreme stability of emulsions of sea water in fuel oil to very stable asphaltic films. Numerous investigators have observed rigid films in the course of crude oil-water interfacial tension determinations by the pendent drop method. Several investigators5,6,2 have separated inter-facially active materials from crudes and attempted to characterize them chemically. Reisberg and Doscher,2 using Ventura crude, showed the interfacial tension against water (as measured by the pendent drop method) to be affected by aging, contraction and expansion of the interface, and the pH of the water. These investigators attributed the adhesion of oil to a water-wetted surface and the distortions of flow paths in glass capillaries to the presence of rigid films. Dodd7 has studied the interfacial viscosity of adsorbed films and found them to be non-Newtonian in behavior. Craighead and Harvey8 reported a series of displacements in tubes packed with 60 mesh glass beads. They interpreted the results as indicating an effect of stearic acid films on waterflood recovery and imply that natural films may produce similar results.

Jan 1, 1967

By D. E. Baldwin

A method has been devetoped for predicting the produced concentration profile for a miscible slug in a five-spot pattern. The technique consists of dividing the five-spot into radial etements and applying an approximate radial solution of the dispersion equation to each element. To test the method, laboratory experiments were conducted in five-spot models. The models were 2-ft X 2-ft X 2-in. Berea sandstone slabs simulating one quarter of a five-spot. An aqueous system was used with tritiated water as the tracer fluid. Effiuent concentration profiles were obtained for slug sizes ranging from 0.7 to 10.2 per cent of the pore volume. Agreement between predicted and experimental profiles was excellent. A prediction was also made for a field tracer test. For this case the predictive technique was modified to account for stratification by use of a layer cake model. INTRODUCTION Dispersion in porous media is of growing interest to the petroleum industry because of the increasing importance of secondary recovery operations. One of the controlling factors in the recovery of oil by miscible displacement or by the use of waterflood additives is the degree of mixing between the fluids of interest. Consequently, one of the prerequisites of any recovery prediction is an adequate method of predicting dispersion. Also, any quantitative interpretation of waterflood tracer tests is dependent upon the ability to determine the dispersion of the tracer as it flows through the reservoir. During the past few years there have been reported in the literature the results of a large number of investigations, both theoretical and experimental, concerning dispersion in porous media. (An excellent review of these results has been presented by Perkins and Johnston." Only a few of these investigations, however, deal with non-lin-ear flow systems. Raimondi et al.9 have presented an approximate solution of the dispersion equation for radial, diverging flow. Lau et al. 7 and Bentsen 1 have presented laboratory data for radial systems and have shown that the dispersion can be adequately predicted by Raimondi's solution. Bentsen has also presented a small amount of data for a radial, converging system. Brigham and Smith' have applied Raimondi's solution, with approximating assumptions, to a five-spot pattern to predict the behavior of a field tracer test. The afore-mentioned investigations have been concerned almost entirely with radial diverging flow and have neglected the effect of converging flow. For this reason the present investigation was undertaken to develop a method of predicting dispersion in a five-spot pattern including the effects of both diverging and converging flow. The method which was developed consisted of dividing the five-spot into diverging-converging radial elements and applying an approximate solution of the disperson equation to each element. DEVELOPMENT OF PREDICTION METHOD The development of the prediction method consisted of two phases: construction of a radial model of a five-spot pattern and derivation of equations to describe the dispersion in the radial model. FIVE-SPOT MODEL A radial model which simulated the flow characteristics of a five-spot pattern was constructed by dividing the five-spot into diverging-converging radial elements (Fig. 1). Because of symmetry the model represented only an octant of a complete five-spot. Each element had an included angle of l° (making a total of 45 elements) and a radius r, such that the sum of the individual element areas was equal to the area of the five-spot octant. The normal five-spot fractional flow behavior was determined by averaging the experimental five-spot data of Caudle and Witte,' Dyes, Caudle and Erickson 5 and Fay and Prats. 4 The resulting curve is shown as a solid line in Fig. 2. To match this curve with the radial model it was

Jan 1, 1967

By J. F. Evers, S. Sadiq, F. W. Preston, G. W. Swift

Proper scaling of unsteady-state flow in laboratory gas reservoir models can be achieved if the permeability/fractional porosity ratio is less than 0.5 md. Porous media meeting this specification can be obtained from mixtures of sand and Portland cement provided (1) the water/cement ratio is in the range of 0.35 to 0.40, (2) the sand/cement ratio is in the range of 1/1 to 2/1, (3) the sand grain size is in the range of 48-100 to 18-28 mesh, (4) the mortar curing time is 7 days or more and (5) the cured mortar is dried under vacuum to reduce the temperature level employed during the drying process. The areal homogeneity of the mortars tested is reasonably good, although significant trends in vertical permeability were noted. The Klinkenberg and non-Darcy coefficients of mortar are such that proper scaling of these effects in the laboratory model could be expected. INTRODUCTION A long term deliverability prediction for a gas well is generally based on data determined from flow tests that are of short duration relative to the producing life of the well. These data, together with an estimate of gas reserves in place, form the basis for establishing the parameters of the mathematical model used to predict performance at future times. Results obtained from a mathematical model represent only an interpretation of the physical process being modeled. Experimental verification of the production process from discovery to abandonment is needed where production occurs under carefully controlled conditions. Ideally, one might consider such an experiment to be conducted in the field. Even if the cost of such an experiment could be justified, the time required would make it impractical to consider. However, if the time scaling problem could be handled properly, an experimental model of the gas reservoir production process could provide the necessary data to study the mechanisms of the production process in a reasonable time period. An experimental model has been designed and constructed using the data of this study as a basis. The model reservoir system consists of one producing well located at the center of a right-circular cylinder of porous material. Ten to 15 observation wells are located at strategic offset positions from the producing well. These wells, and the producing well, are fitted with pressure transducers so that the pressure-time history can be recorded at various radial and azimuthal positions in the drainage area during the entire production life of the system. Time is scaled in the model so that the time required to produce gas from discovery pressure to abandonment pressure is about 2 hours. This compares with a 20-year producing life in the model's field counterpart and gives a time scale factor of approximately 1:100,000. DESIGN CRITERIA To establish design criteria for the experimental model, the reservoir is considered to be a finite right-circular cylinder of radius rr, and height h with a producing well of radius rw located at the center and completed through the producing formation. This reservoir consists of a porous medium which has a uniform permeability k and gas porosity and the gas contained in the reservoir has density p, viscosity µ and molecular weight M. Initially, the pressure at all points within the reservoir is uniform at po. and the reservoir temperature is T. Since the reservoir production process has been demonstrated to occur in two distinct transient modes for either Darcy or non-Darcy flow,l it is

By W. B. Gogarty

With flow of non-Newtonian fluids in porous media, effective viscosities are needed for use in the Darcy equation. These viscosities depend on the rock parameters and flowing conditions. In this investigation, rheological and flow results were correlated to develop an expression for calculating effective viscosities. Surfactant-stabilized dispersions of water in hydrocarbon were used with consolidated sandstone cores of different permeabilities. The average shear rate in a core was related to the permeability and porosity of the core and the frontal velocity through the core. By using the correlation obtained for determining effective viscosity, experimental and theoretical values were compared for each fluid system by determining average and maximum errors. For the fluid system showing the largest average error, the correlation fit 38 experimental points with an average and maximum percentage error of 4.3 and 14.8, respectively. For the fluid system showing the least error, average and maximum percentage errors of 1.8 and 4.3, respectively, were calculated using 29 experimental points. INTRODUCTION Use of non-Newtonian fluids in the petroleum industry is not new. Fluids of this type have been used for many years as fracturing agents and drilling muds. Recently, attention has focused on the use of high molecular weight polymer solutions for secondary recovery. Results have been reported on studies relating to the flow behavior of polymer solutions in porous media.1-4 Each study gives the viscosity in the Darcy equation as a function of the rheological properties of the fluid, the characteristics of the porous medium and the pressure gradient. The functional forms of the reported expressions differ somewhat in each study. The first three studies used unconsolidated porous media while the fourth used a single consolidated core. Because of the limited information available on the flow of non-Newtonian fluids in porous media, further investigation was undertaken in this area. We selected non-Newtonian fluids that did not contain high molecular weight polymers in order to minimize core plugging. The over-all objective was to gain a better understanding of the flow behavior of non-Newtonian fluids in consolidated porous media. Specifically, our goal was to determine the relationship between the average shear rate in a core and the frontal velocity. Because of this goal, it was necessary to determine the effect of permeability on the shear rate-frontal velocity relationship. We wanted to replace the viscosity in the Darcy equation with an expression relating to the average shear rate in the core and the rheological properties of the non-Newtonian fluid. THEORY For a time-independent, non-Newtonian fluid, the apparent viscosity is a function of the shear rate. This relationship is given by µ = F(v')...........(1) With Eq. 1, a single value of shear rate determines the apparent viscosity at that shear rate. With flow in porous media, a spectrum of shear rates occurs at any average throughput. An average shear rate v' can be used to represent the spectrum of shear rates that exists within the core. By using v' and the rheological data given by Eq. 1, an effective viscosity is determined for the flowing fluid. If v' is the correct average, µe becomes the appropriate viscosity for use in the Darcy equation. Under these conditions, flow rate-pressure drop relationships can be determined for non-Newtonian fluids. HOW the correct average shear rate might be related to the frontal velocity through the core and the characteristics of the porous media will now be considered. Later, the experimental data are presented that indicate the validity of these

By N. Mungan, F. W. Smith, J. L. Thompson

Adsorption of polymers and transport, rheology and oil recovery efficiency of their solutions were studied in the laboratory to evaluate the use of polymers in waterflood-ing. While a tenfold mobility reduction was obtained with polymer concentrations as low as 0.05 per cent by weight, the mobility reduction depended on the type of polymer, molecular weight, salinity and pH of water, crude oil and capillary properties of the porous media. Choice of a suitable polymer and a workable concentration will have to be tailored for each application. Little reduction in the residual oil saturation can be expected from polymer flooding. Improvement in the volumetric sweep efficiency is possible hut the extent of the improvement can best be evaluated by properly designed field testing. Some aspects of the field use of polymer floods are discussed. INTRODUCTION Waterflooding is a simple, inexpensive secondary recovery method and is being used widely. Innumerable laboratory studies have been made to unravel the fundamentals of the displacement of oil by water and to find the ways of most efficient oil recovery. These studies and a great many field case histories have revealed that the prime cause of poor oil recovery is the inefficient and incomplete sweep of reservoir volume by the injected water. Sweep efficiency is affected by many factors of which the mobility ratio is an important one. Mobility ratio M is defined here as the ratio of water to oil mobilities: M = (k»/y,r)/(k,JJJJ........(1) In Eq. I, the permeabilities are the effective permeabilities and depend on fluid saturations and, hence, change during the different depletion stages in a flood. A wide practice is to use the effective water permeability at residual oil saturation and the effective oil permeability at interstitial water saturation in Eq. 1. If the mobility ratio is greater than one, the mobility ratio is unfavorable and water, being more mobile than oil, would finger through the oil zone resulting in poor oil recovery efficiency. If the mobility ratio is favorable (one or less) the displacement of oil by water occurs more or less in a pistonlike fashion. In some waterfloods. the mobility ratio is unfavorable and any additives by which the mobility of water can be decreased would favor more efficient oil recovery. The thing to bear in mind, however, is whether or not the improvement in oil recovery is sufficient to more than pay out the cost of the additives needed. For example, materials like sugars, alcohols and glycerine reduce water mobility by thickening the water, but the cost of material requirement precludes any field application. For an additive to be useful in water-flooding, it must bring about a large reduction in water mobility at low concentrations; it must be adsorbed only negligibly; and it must not completely plug up the formation. Some synthetic organic polymers have shown promise of meeting these requirements and have been used in the field.'-W owever, no in-depth studies of the rheological, adsorption and oil displacement characteristics of polymer solutions have been reported. The present work is a study of these properties. EXPERIMENTAL In this work, concentrations are given on a weight per volume basis; 0.5 per cent concentration means 0.5 gm of polymer is dissolved in enough water (or NaCl solution) to make 100 ml. A bactericide, usually 0.1 to 0.2 per cent by volume of 38 per cent formaldehyde solution, was used in the polymer solutions. The NaCl solution was 30,000 ppm. Some properties of the polymers studied are given in Table I. Physical properties of all cores used are in Table 2. Flow behavior of polymer solutions was studied by three consecutive flow tests in cores. First, water (or brine) was injected at constant rate of about 1 ft/D to obtain the water mobility. Then, filtered polymer solution (prepared in water or brine) was flowed through the core. Since the rate was constant, increase in the pressure drop across the core reflects decrease in the mobility. Finally, the core was flooded with water (or brine) to study recovery of mobility. The Alundum cores which were used in

Jan 1, 1967

By G. W. Thomas

Jan 1, 1967

By D. L. Kelley

High-viscosity water injection has been proposed for use in reservoirs containing high-viscosity crude oils. Previous publications have largely ignored the possible effects of the connate water on the proposed process. This paper describes experimental work which indicates that the connate water will be forced ahead of the injected water to form a bank of low-viscosity water. This decreases the oil recovery which would be expected if such a bank were not formed. These effects are shown for a range of fluid mobilities and connate-water saturations for a five-spot injection system. In general, oil recoveries using viscous water are significantly greater than for untreated water even though they are less than would be expected if no connate water bank were formed. INTRODUCTION The effect of mobility ratio on the oil recovery of wa-terfloods has been known for many years. Muskat first pointed out that the fluid mobilities (k/µ) in the oil and water regions would affect the performance of the water-flood, and he estimated the general effect of these variables.' Since this early work, studies of the effect of mobility ratio on secondary recovery have been reported where mathematical,' potentiometric3 and scaled flow models' were used. These studies have shown that a reduction in the mobility ratio between the oil and the displacing fluid would cause additional oil recovery when water-flooding reservoirs containing viscous crude oils. Studies reported by Pye- nd Sandiford 8 have indicated that chemicals to increase injection water viscosity are now available and can be used to reduce the over-all mobility ratio of a waterflood. Where mobility ratios are controlled by the injection of viscous fluids, the connate water of the reservoir can play an important part in the displacement of the reservoir oil. The purpose of this study was to determine the effect of the connate-water saturation in waterfloods where viscous waters are used for injection. DISPLACEMENT OF THE CONNATE WATER Russell, Morgan and Muskat7 were the first to recognize the mobility of connate waters in waterflooding. They conducted waterfloods on oil-saturated cores containing 20 and 35 per cent irreducible water saturations, and found that from 80 to 90 per cent of the "irreducible" water was produced after only one pore volume of water was injected. However, their experiments were conducted at rates of flow significantly higher than those ordinarily occurring in waterfloods. Also, the cores were only from 4.0 to 8.5 cm long. Brown 4 studied a 100-cm linear sand pack which had been prepared to contain connate water and oil. He used 140- and 1.8-cp oils with injection water of essentially the same viscosity as the connate water. He found that all of the connate water was displaced by the injection water in both cases. However, the injection volumes required for complete displacement of the connate water were considerably higher in the case of the more viscous oil. To verify the results of the foregoing experiment, a 10-ft-long linear model was constructed by packing 250-300 mesh sand in a 1/2-in. diameter nylon tube. The model was evacuated, saturated with a brine of 1-cp viscosity, and flooded with a 41-cp mineral oil to the irreducible water saturation of 10.9 per cent. The model was then waterflooded by the injection of a water solution which had an apparent viscosity of 42.6 cp. The solution consisted of 0.5 per cent methylcellulose in distilled water. The viscosities of the oil and connate water were measured with an Ostwald viscosimeter. The viscosity of the polymer solution was calculated by Darcy's law using pressures measured during actual flow conditions. The ratio of the mobility in the oil region to the mobility in the inject ion-water region was approximately 0.32. The mobility ratio of the oil region to the connate-water bank was approximately 14. The mobility ratio between the connate-water bank and the injection water region was 0.024. Approximately 84.5 per cent of the recoverable oil was produced before water breakthrough. Immediately following breakthrough, oil and connate water were produced at an increasing water-oil ratio until the viscous injection water broke through. At viscous-water breakthrough, 96 per cent of the original connate water had been produced. After breakthrough of the viscous water, there was no additional production of connate water or oil. The near-

Jan 1, 1967

By S. S. Marsden, S. A. Khan

Externally generated foam was injected continuously into short porous media. Both flow rate and pressure drop were measured. Liquid saturation was determined by electrical conductivity. Foam yuality I, expressed as the ratio of gas volume to total volume, was varied from 0.70 to 0.96. As measured with a modified fann VG Meter, apparent viscosity of this foam µa decreases with increasing shear rate but usually falls within the range of 50 to 500 cp. At a given shear rate, µa increases almost linearly with I-. When measured with a Bendix Ultraviscoson, kinematic pa is independent of V but absolute pa increases with T from about 3 to 8 cp. The effective permeability -apparent viscosity ratio ke/µa decreases almost linearly with V for porous media of high permeability, but the rate of decrease becomes less for tighter ones. The relative permeability-apparent viscosity ratio kr/µa us V data does not fall on a single line. The kr/µa ratio increases with liquid saturation in the porous medium and with surfactant concentration. Estimates of pa for foam in porous media vary from 30 to 100 cp. INTRODUCTION Although research on the development of a foam-drive, oil recovery process has been going on for almost a decade, most of the significant publications have appeared within the last several years. This illustrates well the rate at which interest in this process is accelerating. Bond and Holbrook 1 were the first to describe the use of foam to improve oil recovery in their patent of 1958. They proposed that an aqueous foaming agent slug be injected into the formation and that this be followed by gas to produce a foam in situ. Fried 2 studied the injection of foam into porous media which had already been subjected to conventional gas or water drives and found that gas could be used to drive a foam bank which would, in turn, displace additional oil in the form of an oil bank. He attributed the increased oil recovery to the high effective viscosity of foam flowing in porous media. His microscopic observations showed the importance of foam generation and regeneration within the porous medium. By injecting both air and aqueous surfactant solution, Bernard3 generated foams within the porous medium in which oil displacement was being studied. In a separate empirical test, he also measured the dynamic foaming characteristics of the same surfactants in water and/or oil. With some exceptions and for the seyen surfactants studied, there seems to be a qualitative relationship between the efficiency of liquid rlisplacrment and the dynamic foaming testtsed. This relationship was not consistent enough to eliminate the necessity of actual foam flood tests in porous media for surfactant selection. In a study basic to gas storage in aquifers, Bennett4 described the displacement of brine by foam in consolidated porous media. Among other things, he stated that the ability of a surfactant solution to foam is more important than the stability of its foam. The presence of a foam bank between the displacing air and the displaced brine improved both breakthrough and ultimate recovery. In a continuation of this work Kolb5 attributed the great reduction in surfactant solution production rate as displacement by air progressed to a decrease in relative permeability to gas. These several effects reported by both Bennett and Kolb can all be attributed to the high apparent viscosity of foam which was obviously flowing in the porous media. Bernard and Holm6 observed that gas permeability in the presence of foam was much less than in the absence of foam when both were measured at the same gas saturations. This reduction in permeability was proportionately greater for more permeable sands and sandstones than for the less permeable. Viscosities of the surfactant solution and air were used for their calculations rather than apparent foam viscosities. The foam stability increased with decreasing permeability, which was probably due to the natural regenerative properties of porous media, and the permeability reduction increased

Jan 1, 1967

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