Search Documents

By E. M. Boatman, S. J. Pirson, R. L. Nettle

More than a decade ago some theoretically derived relationships were proposed that permitted the prediction of the relative ability of reservoir fluids (oil, gas, water) to flow simultaneously within the porous structure of rocks. While many successful predictions of fluid-flow ratios from wells have been made from electric logs, the theoretical deductions have not been subjected to laboratory verifications. For the relative permeability to liquid and gas, simultaneous flow of water and air was used, employing the familiar external-gas-drive technique. While in the process or water desaturation, the resistivities of the test cores (Woodbine sand) were measured. Their permeabilities varied from 10 to 280 md. The laboratory results indicated the need for changing the values of exponents in the theoretically derived formulas, whereas the form of the functions of resistivities and saturations remained as derived by theory. For the relative permeability to oil and water, the three-core dynamic technique was used on six cores. In this case it was found that the formulas were valid with but slight modification as derived by theory, to a sufficient degree of accuracy for engineering use and for well productivity prediction from electric logs. INTRODUCTION Conventional quantitative well log interpretation generally stops after reasonable values for formation porosities and water saturations have been established. Yet such in formation is insufficient to predict either the outcome or the desirability of well tests and completion, let alone the economic aspect of producing oil- and gas-bearing sections which might have been established by conventional log analysis. To give an engineering answer to such questions, it is necessary to predict the expected fluid flow ratios, the productivity index and the ultimate recovery from potentially productive wells. This problem requires in part for its solution a determination of the relative permeability characteristics of the potentially producing formations, both to oil (or water) and gas, and to oil and water.

Jan 1, 1965

By D. M. Kehn

A method has been devloped for chromatographic analysis Of the vapor and liquid phases Of a a system containing methane to components having 20 or more carbon atoms. The method uses a windowed equilibrium cell in which volumetric phase behavior of the system can be observed accurately and from which small samples of gas or liquid can be withdrawn for analysis. Analyses are made using two chromatographs, one for the lighter and one for the heavier components in a sample. Combination of the two analyses yields a detailed analysis of the gas or liquid sample. The complexity of the condensate heavy ends is evident from the chromatograms of these fractions, and the predominance of the paraffin hydrocarborn serves as a useful marker in interpreting the chromatograms. The K-values obtained in this analytical method are presented for a high-pressure condensate system and predict closely the observed volumetric behavior of the system. INTRODUCTION Quantitative analysis of hydrocarbons from natural gas reservoirs is necessary for several reasons—to calculate the amount of sales gas produced, to calculate the amount of natural gasoline produced, to plan a liquid recovery system, or to calculate the potential economic value of a reservoir produced under one or more of several different conditions. Analysis of natural gas fluids produced to the surface consists of identifying and computing the mol fraction of each component of the mixture. Although methane is the predominant component, varying amounts of ethane, propane, butane, pentanes and heavier components are also present. Materials containing up to 30 carbon atoms occur in amounts which decrease with increasing molecular weight. However, the quantities of components in the 20 to 30 carbon atom range are usually so small that their importance is negligible, and they are undetect-able in natural gas by ordinary analytical methods. All the components up to those having 20 carbon atoms may sigsficantly affect phase behavior, however. Commonly, only the methane-through-pentane fraction is analyzed quantitatively for each component, while components heavier than Pentane are lumped and repored as "hexane- plus". Expensive, tedious techniques are required for analysis of this fraction. Consequently the detailed analyses needed for prediction of reservoir behavior are usually undertaken only when major gas fields are being developed. The need for complete analyses of condensate systems is apparent when it is recalled that most gas fields are produced by pressure depletion. As the pressure declines, some of the heavier hydrocarbons are lost as liquids which are in the reservoir. In many instances the amount of liquid in equilibrium with the gas phase at high pressure constitutes only 1 or 2 mol per cent of the total system. Flash calculations generally must predict the actual amount of liquid with an accuracy of a few per cent in order to be useful. This retrograde condensation has been understood for years, but accurate correlation methods to permit quantitative prediction of phase behavior in the retrograde region are not presently available. The increasing importance of natural gas has made accurate prediction of phase behavior and composition of produced natural gas streams an economic necessity. The work reported here was undertaken to provide a rapid, economical method for obtaining the vapor-liquid equilibrium information needed to predict accurately the composition of the fluids produced from a gas reservoir throughout its life. TO develop this method, a pressure cell equipped with windows was designed and built for observing the volumes of liquid and gas present at reservoir pressures and temperatures. Use was made of established chromatographic methods for rapid and detailed analysis of both phases. This paper describes the equipment and techniques developed for obtaining vapor-liquid equilibrium data, presents the results of analyses of a condensate system, and indicates the usefulness of these data in predicting hydrocarbon phase behavior. DESCRIPTION OF EQUIPMENT USED The equipment used in obtaining the required information on phase behavior and the complete analysis of hydrocarbon mixtures through C, will be described first, followed by a discussion of the operation of the equipment. It will be helpful, however, to consider first a brief outline of the technique used. A sample of separator gas and liquid is charged to the windowed cell (see Figs. 1 and 2) where volumetric equilibrium phase behavior at reservoir pressures and temperatures can be determined. Then samples of the coexisting phases are withdrawn. The methane-through-pentane fraction is analyzed with a chromatograph equipped with a hot-wire detector, and the pentane-plus fraction is analyzed with a second chromatograph equipped

Jan 1, 1965

By V. S. Gupt, W. H. Somerton

Rock may undergo great changes in physical properties when heated to high temperatures and then cooled, The temperature and intensity of reactions causing rock alterntiorl.s can he controlled by introducing certain chemicals during heat treament. Three typical outcrop sandsone samples were saturuted with common salt solutions, then heated to several maximum temperatures. After cooling, it was found that per-rrreuhilities had increased much more for salt-saturated samples than for samples not saturated with salt but heated to the same temperatures. This was only true, however. for samples heated above the melting points of the particrrlar salts. Potassium chloritle was particularly effective with Bandera sandstone. Samples saturated with potassium chlorirlr. solution and heated to 900C showed an 11-fold increa.se in perrrieahility. Samples without potasrirdrn chloritle but heated to the .same temperature. showed only a 2.7-jold incresre in permecahility. In application, it seems possible that injecting chemicalc into the formations from a wellbore followed by applying intensive borehole heating might promote reactions which would greatly improve permeability of the formalions. INTRODUCTION Great changes in the physical properties of rocks which have been heated to temperatures in the range of 600 to 800C have been reported earlier. Permeability increases of 50 per cent, and equivalent decreases in sonic velocity and breaking strength have been observed. Although there might be some suppression of certain reactions responsiblc for these changes when rock samples are heated under simulated reservoir pressure conditions. recent work has shown that these are more than offset by the increased importance of other reactions at high pressures. The reactions considered responsible for alteration of rock properties by heating include differential thermal expansion. dehydration, phase changes and dissociation of mineral constituents. Ceramists have long known that temperatures at which reactions occur can be raised or lowered by the presence of certain impurities in the sys. tcnl. Thus the possibility exists that by saturating rocks with appropriate solutions. desired reactions might be accelerated and unwanted reactions might be suppressec! when the rock is heated. Purpose of the present work was to investigate the effects of several common salts on the thermal alteration of sandstones. Types of reactions which might occur are reviewed and the findings of a number of thernlochemical alteration tests are presented. THEKMOCHEMICAL ALTERATION OF SANDSTONE MINERALS The most abundant constituent of most sandstones is quartz. Quartz can exist in at least eight different forms, but in the temperature range of immediate interest (to 900C) only two forms are of importance—alpha quartz below 573C and beta quartz above this temperature. Although the alpha-beta inversion of quartz is very important in thermal alteration of rocks. it is not particularly important here because the inversion temperature cannot be changed significantly by the presence of impurities or by the application of pressure. The transition temperature of quartz to tridynlite is 867C. but the transition is very sluggish. Tridymite is considered to be a metastable transition phase between two stable phases: quartz?cristobalite. However. the transition temperature for tridymite to cristobalite is 1,470C.4 The quartz-tridymite-cristobalite transition can be accelerated substantially by the presence of fluxing agents. Finely divided calcium, magnesium and titanium oxides accelerate the conversion while A12O may retard it. A Mixture of NaCl and CsCl reduces the temperature for cristobalite formation by as much as 420°C. Feldspars constitute the second most important mineral in sandstones. The melting points of feldspars range between 1,000 and 1,700°C depending upon the variety. No information available indicated any important role of fluxing agents in the thermal alteration of this group of minerals. Carbonate minerals are subject to dissociation at temperatures within the range of the present investigation. The dissociation of magnesite can start at a temperature of 373°C, but the reaction is sluggish and might not occur untll a temperature of 500°C or ,higher is reached. Dolomite dissociates in two stages at 500°C and 890°C, whereas calcite dissociation temperature is about 885C. Because CO2 is released in the dissociation of carbonates, ail such reactions are somewhat dependent upon CO2 partial pressure. In the absence of CO: in the surrounding atmosphere, the dissociation of calcite starts at 500°C.However. when 1 atm of CO surrounds the sample, the dissociation does not start until 900C. The other carbonates are apparently much less sensitive to a change in partial pressure of CO2. The aragonite calcite transformation can occur an!. where within the temperature range of 357 to 488°C. depending upon the presence of impurities such as barium. strontium. lead and perhaps zinc. The differential thermal analysis (DTA) curves of both magnesitc and dolomite vary with the presence of 'impuritics. The presence of iron in particular seems to affect the magnesitc DTA curve

Jan 1, 1966

By N. Mungan

Laboratory water floods were performed in oil-wet and waterwet alundum and Torpedo cores, displacing a refined oil with n-hexylamine or Triton X-100 solution. Also, some floods were performed in which a sucrose solution was displaced with n-butyl alcohol. The Purpose ot the tests was to see if the oil recotlery could be increased and to examine what role the inter-facial tension reduction and wettability change play in the recovery mechanism. it was found that the chromatographic transport of amine was influenced by core wettability. In oil-wet cores, the rate of advance of the amine band could be predicted from equilibrium chromato-graphy, while in water-wet cores the amine band moved faster than predicted, indicating a non-equilibrium process. By reducing the interfacial tension to 1.1 dyne/ cm, oil recovery was increased. More oil was recovered from Dri-filmed cores than from water-wet cores. Reversing the wettability of the porous media from oil-wet condition also resulted in some additional oil recovery. Neutral wettability floods did not increase oil recovery. INTRODUCTION Interfacial forces in petroleum reservoirs are responsible for retention of large quantities of residual oil. Increasing exploration costs have created an incentive to attempt recovery of this residual oil by new and novel processes. One novel process involves changing the interfacial forces by introducing a chemical into the reservoir during water flooding. Recently, several investigators have carried out laboratory displacement tests using amines. lt3.4 The amines adsorbed onto initially water-wet core surfaces, changing them to oil wet. Additional oil was recovered when the cores were flooded in a manner reversing the. wettability from oil-wet conditions. In all these studies, the wettability reversal was considered responsible for increased oil recoveries although it was noted that amines also reduced oil-water interfacial tension. In the present study, displacement experiments were performed in water-wet and Dri-filmed cores. The purpose was to separate the effects of lowering the interfacial tension from changing the wettability. LABORATORY STUDY All experimental data given in this paper have been obtained at a temperature of 70F + 1 and prevailing atmospheric pressure. FLUIDS AND MATERIALS The oil used in these experiments* was washed with H2S04 and passed through two 100-200 mesh silica-gel beds to remove unsaturates and surface active impurities. Normal hexylamine** of practical purity was used, and the alundum cores*** were composed of nearly pure A1203. Fused quartz (SiOi) and Lucalox (A12 )t plates were used in the contact angle measurements. DISPLACEMENT TESTS Displacement tests were performed in alundum and Torpedo cores. Before displacement runs, cores were acidized with 0.1 normal HCI, flushed with distilled water and fired at 1,200F to make them water wet. This procedure also reduced water sensitivity of the Torpedo cores. Oil-wet cores were prepared from water-wet cores by flowing several pore volumes of 1 per cent by weight Dri-film SC-87 solution in heptane through the cores, flushing the cores with nitrogen, evacuating, and finally heating the cores at 300F for 4 or 5 hours. Between floods, cores were cleaned by conventional procedures and treated again as above.

Jan 1, 1965

By F. M. Perkins, R. H. Jamison

Publications concerning scaling laws for laboratory flow models of oil reservoirs indicate that the relative permeability and capillary pressure relations must be the same functions of saturation in the model and its prototype. In this paper, the relative permeabilities and saturations are redefined in a way which permits one to have diflerent relative permeability and capillary pressure relations in the model and prototype. The development of the new scaling criteria is demonstrated by deriving scaling laws for one simple flow problem. INTRODUCTION The successful application of information gathered through bench-scale model studies to the prediction of multi-phase fluid flow behavior in petroleum reservoirs requires careful consideration of the scaling criteria. The requirements which must be met for a small-scale model to duplicate the behavior of a much larger reservoir prototype have been discussed by Leverett, et al,' Rapoport,2 and Geertsma, et al.3 According to these authors, the relative permeabilities to water, oil and gas and the dimensionless capillary pressures must be the must have the same values in the model and prototype. Furthermore, they point out that the viscosity ratios must have the same values in the model and prototype. In addition to these requirements, other scaling criteria are given which, when considered, usually dictate a much higher permeability in the model than exists in the reservoir. Since the more permeable model sands may have connate-water and residual oil saturations quite different from those in the prototype, the scaling requirements pertaining to the dependence of relative permeabilities and dimension l ess capillary pressures on saturation usually are not met. Also, if the model fluids are selected to give equal viscosity ratios in the model and prototype, one may find greatly different mobility ratios* in the two cases and, therefore, may expect grossly different flow patterns. For this reason, Craig, et al,' correlated the results of their model studies with mo- bility ratio rather than with viscosity ratio. As yet, however, no satisfactory method has been advanced which permits use of anything other than the scaling requirements outlined in the previously mentioned papers. The purpose of this paper is to present a method of developing a set of scaling criteria which permits different relationships between saturation and relative permeability or capillary pressure in the model and prototype. Rather than attempt to report scaling criteria for a variety of practical situations, this paper demonstrates by use of one simple example how the modified scaling criteria are derived.

By J. M. Dumore

If in a vertical, downward miscible displacement, the transition zone between the displacing and displaced fluids is neglected, a criterion for stable displacement can be obtained by considering a small hypothetical protrusion of one of the fluids into the other. This criterion leads to the definition of the well-known critical rate, uc = kg p/p. The consideration is further extended by taking into account the transition zone that develops as a result of diffusion and mixing. A generalization of the previous criterion leads to the definition of another characteristic rate, the stable rate, which in actual miscible drives will be less than the critical rate. In such drives, the entire transition zone is stable at rates less than the stable rate. At rates between the stable and critical rates, the displacement is only partly stable, i.e. part of the transition zone adjacent to the displaced fluid is in an unstable position. From that part of the transition zone viscous lingers will develop. At rates greater than the critical rate the entire displacement is unstable and viscous lingers will develop more strongly. Results of laboratory experiments are in agreement with the expected behavior based on the theoretically deduced stability of the displacement. INTRODUCTION The simplest form of miscible drive in an oil-bearing formation is the injection of a fluid that is completely miscible* with the oil under reservoir conditions. In general, such a fluid, a solvent for example, is less dense and less viscous than the oil present in the formation. If it is injected into a horizontal homogeneous layer, gravitational forces will lead to the formation of a gravity tongue of solvent in the upper part of the layer and the adverse solvent-oil viscosity ratio will cause viscous fingers to develop. If, however, the solvent is injected up-structure into ,a dipping layer, gravity has a favorable effect, because it tends to keep the less dense solvent up-structure. Tongue formation and viscous fingering are consequently reduced and it is even possible that they will be suppressed completely. Viscous fingering and gravity tonguing are the consequences of the instability of the displacement. A stable displacement cannot result in growing viscous fingers and/or growing gravity tongues. Since large amounts of oil can be bypassed if there is viscous fingering and/or gravity tonguing, the stability of a miscible drive is very important with respect to the recovery efficiency of the drive. The stability is of particular importance in miscible-slug drives, as it determines how quickly the miscible slug between the displaced and displacing fluids will be distorted and broken up, after which the drive is no longer completely miscible. Stability is thus a most important factor in determining the success of a miscible drive, and it is considered that the aspects of stability considered in this paper will make a useful contribution to existing theories. Consideration is given only to vertical downward displacements, in which no gravity tongues can develop and which are therefore simpler than downward displacements in sloping layers. THEORY CRITERION FOR STABLE DISPLACEMENT Let us consider a vertical column of a homogeneous permeable medium saturated with oil. A solvent less dense and less viscous than the oil is injected at the top of the column at a constant rate; in the column the velocities of oil and solvent are assumed to be uniform. Diffusion and mixing of oil and solvent will at first be neglected, i.e. we assume that solvent and oil are divided by a sharp interface. An additional assumption is that the interface is horizontal. Continuity of pressures exists at the interface, because no capillary forces are present between the two miscible fluids. Consequently, at the interface (Z = I)):** Ps = Po - P-where Ps is the pressure at the solvent side and

Jan 1, 1965

By R. L. Slobod, C. Gatlin

This .study defines the basic mechanism of the mis-cible displacerrzent of oil and writer from porous Medici by various water-driven alcohol .slugs. Three distinct alcohol slug processes were .studied. Considerable data concerning the quanity of alcohol required for oil recovery were also obtained. All data were obtained in a 1-in. diameter, 100-ft long, unconsolidated core. The porosity of this system was 35 per cent. and the permeability was approxitrlately 4 darcies. Total core pore volume wax 5,716 cc. All displacements were conducted at a constant injection rate of 5 to 6 cc//min, which corresponds to a frontal advance of 5 ro 6 ft/hr. The first portion of this paper is concerned with 111 use of one alcohol—isopropyl—a\. the slug material. Isopropyl alcohol (IPA) is completely miscible with both oil rind water; however, miscibility of the threc-component systems, oil-water-IPA, requires a relutively. high concentration of IPA. Hence, the displocement is not of the misoible type unless the IPA concentration is maintained above some critical value. A slug of IPA equal to only 13.5 per cent of the pore volume was jound to be sufficient to obtain complete recover?, of residua1 naphtha. In later studies two distinct process variations were developed. The firsr of these utilized methyl alcohol (MA) and IPA as slug materials. It was shown that methyl alcohol may he substituted for IPA at the front and rear of the slug with no loss of oil recovery. A .slug of 4 per cent MA-4 per cent IPA-4 per cent MA was sufficient for complere oil recovery. Because MA is considerably cheaper than IPA, this represents an important step toward economic application. A .second process variation used normal butyl alcohol (nBA) and MA as the composite slug, the nBA segnzent being injected first. This technique requires the smallest total slug size (approximately 10 per cent) of all processes studied. The high cost of nBA, however, precludes commercial application. It is possible that this basic process, subject to changes of alcohol type, may lead to a commrercicil process. INTRODUCTION Within the last 10 years, considerable study has been devoted to the general mechanics of miscible phase dis- placements in porous media. These studies have. in general, dealt with two-component displacernent. All reported field trials to date have been modified gas drives of one type or another. These utilize either lean gas at high pressure, rich gas at lower pressure or an LPG slug. The irreducible water content of the reservoir remains immobile during these processes. Although these techniques hold great promise. certain drawbacks do exist. Poor areal sweep efficiency, inherent in any displacement having a highly unfavorable mobility ratio, is one disadvantage. Furthermore, the miscibility of all processes which use gas is dependent on pressure. These pressure limitations may prohibit application to shallow reservoirs. There are also many areas where large quantities of LPG and natural gas are not readily available. These factors emphasize the need for improved techniques. The alcohol slug process is also a miscible-phase displacement process. It differs from the previously mentioned two-component techniques, however, in that both the reservoir oil and water are displaced by the slug. This behavior is a consequence of the miscibility of certain oil-water-alcohol systems. In the simplest case, a relatively small volume (slug) of an alcohol (such as isopropyl) is injected into the system. Water is then used to drive the slug through the medium. Thus. three components—-alcohol. oil and connate water-—exist at the front side of the slug. Miscibility is obtained at some alcohol concentration, this being dependent on the solubility of the particular system. Miscibility is maintained within the slug until the alcohol concentration falls below that value required to maintain miscibility. When miscibility is lost, the process reverts to a water flood. Certain inherent benefits of this technique are apparent. High pressures are not necessary to obtain miscibility in these liquid systems. Furthermore, water is a more desirable driving agent than gas because of the improved mobility ratio and the attendant improvement in areal sweep efficiency. The main disadvantage is, of course, the relatively high cost of alcohols. The use of an alcohol slug to recover reservoir oil is not a new idea." ' However, this process has received only limited study, presumably due to the lack of industrial interest caused by the seemingly prohibitive cost of alcohols. The main purpose of this study was to investigate the mechanism of the displacement of oil and water by various alcohol slugs and to develop process modifications aimed toward possible commercial application.

By W. F. Rogers, J. A. Stewart, P. J. Kalish, E. O. Bennett

Detailed descriptions are given of materials, apparatus and the experimental procedure used to study the effect of bacteria on sandstone permeability. The factors affecting permeability during injection of bacterial suspensions which have been investigated are: (I) concentration of bacteria; (2) core permeability and median pore size; (3) species of bacteria, mode of aggregation and relative size; (4) injection rate or pressure differential; (5) mean pressure; and (6) depth of penetration of bacteria. The investigation demonstrated that bacteria cannot reduce core pertneability to zero and that their effect on permeability is subject to definite limitations. Remedial or permeability restoration studies also were made. Acidization in combination with reverse flow war found to be an effective method for restoring permeability In cores partially plugged with bacteria. The relationship between the bacterial quantities in the laboratory tests and those found in field operations is discussed. The linear laboratory flow data have been translated into terms of field radial systems; these data indicate the most practical methods of maintaining injection rates in the presence of bacteria are to increase injection pressures or hydraulically fracture the formation. INTRODUCTION It has long been known that various bacteria flourish in the above-ground components of oil field water-injection systems. The presence of these bacteria has led to the suspicion that they might enter into the causes for reduction of injection well permeability, but only recently has an attempt been made to describe the effects quantitatively.1 Biocides are used in many cases without a firm knowledge as to whether the bacteria are killed, or if they are, whether the dead cells affect formation permeability. This project was undertaken as part of a study to determine the effect of bacteria, residual oil and precipitated solids such as iron sulfide or calcium carbonate on the permeability of sandstones and to learn whether any of these materials adversely affect brine injectivity in secondary recovery or disposal operations. This report deals with that portion relating to the effect of bacteria. MATERIALS CORES Three different permeability ranges of Berea sandstone cores were used in the study. The range of initial, absolute brine permeabilities kj were: high — 278 to 400 md: medium—130 to 162 md; and low—17.7 to 48.3 md. The permeability data for the cores used in the tests are included in Table 1. All cores were cut to a nominal 1 in. diameter. Following air permeability determinations, they were molded with an epoxy resin in 1.5 in. ID aluminum sleeves. The sleeves had drilled and tapped holes for intermediate pressure connections along the length of the cores. After molding, the high- and medium-permeability cores were trimmed to 4 in. in length and the low-permeability cores to 2 in. These lengths resulted in pore volumes of approximately 10.9, 10.4 and 4.4 ml, respectively, for the high-, medium- and low-permeability cores. Pressure-transmitting channels were drilled through the plastic into the core proper, using the sleeve holes as drill guides. Width and depth of channel penetration into the core proper were approximately 0.03 in. and 0.05 in., respectively. Three pressure taps were used with all cores; for the high- and medium-permeability cores the taps were spaced 1 in. apart; for the low-permeability cores the taps were spaced 0.4, 0.8 and 1.2 in. from the inlet end. Pore size distributions were calculated from restored-state capillary pressure curves for representative cores in each of the three permeability ranges. In each case the pore size range with the greatest percentage of the pore space included the very small pores up to 0.5 micron in radius. The median pore radii were in the 5.5 to 6 micron range for the high-permeability cores; in the 4.5 to 5 micron range for the medium-permeability cores; and in the 3.5 to 4 micron range for the low-permeability cores.

Jan 1, 1965

By R. Wiesenthal

A method is developed for improving the low recovery efficiency which results when viscous oils are flooded by water. Viscous oil has been diluted with a lighter liquid miscible in it in any ratio which produces a mixture possessing a reduced viscosity and which is then displaced by water. Light gasoline was used in these tests to reduce the viscosity of the oil in place. Oils possessing viscosities in the range of 1.28 to 324 cp at the test temperature of 30C were displaced from a cylindrical unconsolidated sand core, with a length of 2 m (6.6 ft) and a diameter of 10 cm (3.9 in.) at a flood of 6.27 cc/sq cm/hr. In two test series up to I PV of light gasoline was injected into the core before water flooding. In the first series the recovery of a low-viscosity oil (1.28 cp) increased as the volume of light gasoline was increased. However, only the oil which could not be recovered by water flooding alone could be replaced by the light gasoline. In the other series, the recovery of a more viscous oil (13.5 cp) was not increased even if large quantities of light gasoline were injected. An oil bank built up ahead of the water front so that the flood water met with equivalent saturation conditions, as would be the case in flooding solely with water. However, the recovery of more viscous oils could be increased substantially when the core was first flooded with water, followed by a slug of light gasoline, and then flooded with water again. Furthermore, an oil bank ,built up ahead of the second water front. The formation of this oil bank is particularly favorable for the displacement process; the greater part of the light gasoline injected was pushed out of the core ahead of the oil and recovered. Using a light oil of low viscosity, it was possible to replace only the residual oil that could not be recovered by water flooding by the light gasoline. As the basic study had shown that highly viscous oils could be recovered by flooding first with water, then with light gasoline, followed by a second injection of water, the study was extended by tests made with reservoir fluids and cores taken from the low-gravity oil reservoir of the Wietze field, Germany. A great increase in recovery over flooding with water alone was obtained in these tests. The results may show a way to recover high-viscosity oils by a comparatively simple flooding procedure. INTRODUCTION In fields where oil recovery has not benefited from natural water drive, increased recovery may be obtained by the injection of water into the reservoir, especially in the case of fields which produce high-gravity, low-viscosity oil. There is a rather close relation between the viscosity of the reservoir fluid and the recovery that can be obtained by water flooding. As the viscosity increases the oil recovery decreases proportionately, so that an oil viscosity about 30 times that of water is generally considered to be the upper limit for an economically successful water flood.' This limitation has been recognized theoretically by Buckley and Leverett,' and has been verified by numerous laboratory flooding tests, such as by Croes and Schwarz,V n the viscosity ratio range of oil to water from 1 to 500. The reason for the inferior recovery of viscous oils by water flooding can best be understood by considering the mechanism of flood tests. In such tests water is introduced into one end of a sand-packed tube which is saturated with oil of the viscosity to be tested, and oil is produced from, the other end. Oil, with a viscosity lower than that of water, is pushed ahead of the advancing water with considerable uniformity, and irregularities in the displacement process tend to be eliminated as the oil saturation decreases and the water saturation increases. However, oil recovery is not complete, for the water forms interfaces with the oil and traps residual oil droplets in certain of the pore spaces. These residual oil droplets are distributed rather uniformly through the sand, and they are not susceptible to recovery regardless of the amount of water passed through the test cylinder. In the case of an oil possessing a viscosity greater than that of water, the injected water tends to move irregularly through the sand in the test cylinder. At the place in the sand where initial water movement takes place, resistance to the movement of water is lowered and the water finger so formed will grow perceptibly; it will extend to the outlet end of the test cylinder so that a breakthrough of water takes place long before much of the oil-saturated volume of the sand has come in contact with water. Once the breakthrough of water takes place, the water-oil ratio is gradually approaching the economic limit, even though much of the oil in the test cylinder remains unaffected by water. Disconnected drops of oil are trapped in the portion of the sand that has been invaded by water, just as when an oil of low viscosity is flooded. The poor recovery of the more viscous oils, therefore, is due to the irregularity of water movement through the test cylinder, and the rate of water movement through restricted sections of the oil-saturated sand increases with the more viscous oils. There are three ways whereby the recovery of viscous oils by water flooding can be improved:' (1) adding materials which will increase the vis-

Jan 1, 1965

By W. E. Howlett, R. L. Slobod

Jan 1, 1965

By B. Habermann

Artificially consolidated sand models, representing one-quarter of a five-spot, have been developed and used to study factors aflecting misciblt. displacrmenr. Sweep efficiency at breakthrough, size of the mixing zone between two miscible liquids and per cent area contacted by drive after breakthrough were determined tor the high mobility ratios encountered in actual resvrvoirs. Quantitative relationships between the degree of viscous fingering and mobility ratio were obtained by measuring the length of the fluid interface. Scaled miscible-slug experiments, supported by field evidence reported in the literature, have shown that when the slug is followed by dry gas the process is less efficient than expected. In addition to low areal sweep efficiencies encountered for high mobility ratio displacements, the effectiveness of a miscible slug is greatly reduced. This is caused by an accelerated growth of the mixing zone between the driving and the displacing fluids. INTRODUCTION Considerable interest in miscible displacement as a secondary recovery process has been shown by the oil industry in recent years. In addition to extensive laboratory investigations, numerous field operations have been initiated so that 39 miscible displacement projects had been reported' in the United States by early 1959. By far the most popular approach in these projects is the use of an LPG slug followed by dry gas. Laboratory tests have indicated that with a small miscible slug, amounting to several per cent of the hydrocarbon pore volume, practically 100 per cent recovery of the oil in place can be achieved and that all of the slug material can be recovered by a subsequellt gas drive. However, these experiments were conducted on small-diameter long cores that can be considered as essentially uni-dimensional reservoirs that do not represent actual field conditions. It soon became apparent that the inherently adverse mobility ratios existing in the reservoir cannot be ignored in any miscible displacement field project. When a less viscous fiuid such as an LPG slug drives a more viscous reservoir oil, fingers develop that reduce the areal sweep efficiency. The slze of the miscible slug to be used is a subject of some controversy. It has been proposed that the size of siug needed be calculated by setting lengths of the mixing zone proportional to the square root of the distance traveled. This proposal is based on experimental evidence' and theoretical considerations"' in which fingering was not a factor. A prerequisite for successful application of the miscible displacement process to field operations is a better understanding of the following factors. L What is the effect of mobility ratio on the rate of deterioration of a miscible slug? 2. How large a slug should be used for given reservoir conditions? 3. What is the sweep efficiency? The purpose of this investigation was to provide some additional laboratory data that should help to answer these questions. More specifically, it was intended to study more quantitatively the effect of viscous fiinger-ing, mobility ratio and slug size on the efficiency of this relatively new secondary recovery process by means of experiments on thin, horizontal, porous medium models.

By F. W. Jessen, R. W. Durie

The dissolution of salt in the development of salt cavities is controlled by free-convection boundary layer flow along the salt surfaces. It is the purpose of this paper to expand upon results published previously, and primarily to illustrate the influence of surface features on the boundary layer and hence the dissolution process in the development of underground salt cavities. General dissolution equations are presented that apply to a surface inclined to the vertical. A salt surface inclined so as to be overlain by water is dissolved by boundary layer flow, while an inverted surface is dissolved by a combination of boundary layer and cellular flow. Detailed investigation of the boundary layer was limited to laminar free-convection flow. Truly laminar free-convection boundary layer flow may be encountered only in laboratory solution cavities, but the boundary layer will control dissolution even where cavity dimensions are such that free-convection turbulence will occur. The factors found to significantly affect the laminar boundary layer will affect dissolution with a turbulent boundary layer. INTRODUCTION Underground salt cavities have been used extensively for the storage of LPG products, and the use of these cavities is expanding each year. The increasing requirements for such cavities suggest the need for a better understanding of the mechanisms involved, for the purpose of the control of cavity geometry and the most effective use of the injected water. A previous paper outlined a quantitative approach to the determination of salt dissolution rates in the laminar free-convection dissolution of salt. It is the purpose of this paper to expand upon the results previously published and to point out significant factors in the application of these results to the dissolution of field cavities. The dissolution of salt in laboratory scale cavities is controlled by the free-convection boundary layer, which is the region of flow along the salt surface caused by the difference in buoyancy forces resulting from the salt dissolution. The transfer of salt into solution at the salt face causes an increase in salinity, and hence a greater density, in the water at the salt surface. The region of higher density tends to fall through the lighter, less saline water and a flowing boundary layer is the result. The water remains stationary at the contact with the salt surface,l so the resulting velocity profile is such that the influence of gravity is opposed by a viscous shear force at the salt surface. For salt to be transferred into solution from a solid surface to a boundary layer of fluid, it must diffuse through a stationary film of water at the salt contact. According to Fick's first law, the rate of mass transfer by diffusion is the product of the diffusion coefficient and the concentration gradient, with diffusion transfer in the direction of decreasing salinity. Once the salt has passed into the boundary layer, it will be transferred by a combination of molecular diffusion and convection. The convection transfer will in turn affect the rate of dissolution by influencing the gradient of the salt concentration at the surface. In the previous investigation by the authors to determine the basic mechanism of the salt dissolution process, an equation was developed that agreed very well with the experimental rates of free-convection salt dissolution from a smooth vertical surface. The technique was taken from the work of Eckert2 in the analogous heat transfer system and will be referred to as the integral method. The rate of salt transfer from a vertical section of height H was found to be:

Jan 1, 1965

By H. Groeneveld, C. A. Connally, P. J. Hoenmans, J. J. Justen, W. L. Mason

A miscible displacement pilot using a slug of LPG driven by separator gas was conducted in the Cardiurn reservoir of the Pembina field. The injection pattern was a 10-acre, inverted, isolated five-spot. Upon completion of the LPG-gar phase, an experiment was conducted using a slug of water followed by gas. Calculated performance of the pilot is compared with actual performance. Equations are developed to calculate the distribution of LPG into zones of varying permeability, to estimate the progress of the flood at different times in the various zones and to estimate gas rates after breakthrough. The analysis indicates that permeability stratification was a dominant factor in controlling oil recovery and that oil was completely displaced from the swept pore volume. The results of the pilot indicated that miscible flooding is a practical means of pressure maintenance in this reservoir. The total recovery from the pilot area was good in spite of the early breakthrough of LPG. The effects of stratification were reduced by injecting a slug of water into the partially swept reservoir. INTRODUCTION The Pembina field,' located in Alberta, is the largest oil field in Canada and one of the largest in the North American continent. The reservoir is a stratigraphic trap producing from the Cardium sand. Neither bottom water nor free gas has been found. The recovery of oil by the natural depletion mechanism has been estimated at 12.5 per cent. Pressure maintenance studies of various areas have indicated that the recovery can be increased 21/2 times by water flooding, and a large area of the field is presently under water flood. However, reservoir studies of the North Pembina area indicated that miscible flooding might be competitive with water flooding. A pilot test was conducted to evaluate the performance of a miscible flood. A 10-acre, inverted, isolated, five-spot pattern was selected for the pilot. The pattern area was large enough to minimize wellbore fracturing effects and contained sufficient oil to provide significant working numbers. The performance of each of the four producers could be evaluated individually and compared. In the event of breakthrough in one direction, the effect would be isolated from the other producers. The use of a single injector minimized the volume of LPG required, and, because of the high mobility of gas, one well was sufficient to inject the necessary daily volume to replace the high rate of production. With four producers, the test could be completed in time for results to be evaluated, additional engineering studies to be made and a unit to be formed before the reservoir pressure in the North Pembina area declined below the bubble point. The pilot was located in an area developed on staggered, 80-acre spacing. The injection well was drilled at a regular location, while the four producers were drilled 467-ft north, east, south and west of the injector. Each quadrant and its associated producer were identified according to their direction from the injector— that is, north, east, south or west. The eight surrounding producers on 80-acre spacing were shut in to isolate the pilot area and provide for reservoir pressure observation. The pilot wells were completed using permanent-type completion techniques. After coring, casing was run through the pay section and cemented. Inside 51/2-in. casing, 2 1/2-in. tubing was hung. The wells were perforated opposite the Upper Cardium sand and lightly fractured. Fracturing volumes, rates and pressures were low to minimize the extent of the fractures. The fracturing treatments average 1,000 lb of 20-40 mesh sand in 700 gal of a low fluid-loss sand-carrying agent. Feed rates and wellhead fracturing pressures averaged 5.5 bbl/min at 2,535 psig, respectively. After fracturing, the productivity index was measured in each of the five pilot wells. The average PI of the four producers was 0.41 BOPD/psig drawdown. The measured PI'S were approximately the same as PI'S calculated from core analysis data, indicating that the fracturing treatments were just sufficient to overcome

By M. M. Mebta, G. W. Dean, W. H. Somerton

With the advent of underground heating operations, interest has developed in the alteration of rock properties by high-temperature treatment. In the present work a number of sandstones were heated to temperatures in the range of 400°C to 800°C under both atmospheric and simulated reservoir pressures. Pertneabilities increased by at least 50 per cent and sonic velocities and breaking strerlgths decreased by an equivalent amount. Differential thermal expansion and other reactions of constituent min-era1 grains are the causes of these alterations. INTRODUCTION In the underground combustion of petroleum reservoirs, temperatures of the order of 600C are reported to have been reached in the combustion zone.' At this temperature rocks are subject to extensive thermal alteration. Temperatures of this magnitude and higher may also occur in subsurface formations when subjected to bottom-hole heating, thermal drilling operations, and underground nuclear explosions. Temperatures of this magnitude might also be generated by conventional rock drilling methods at points of bit-tooth contact. In earlier work, the permanent deformation of rocks resulting from heating was reported. Major structural damage of rocks occurs due to differential thermal expansion of mineral constituents. A number of mineral alterations including crystal inversions, loss of water of crystallization and dissociation, may also contribute to changes in physical structure and properties of rock. In the present work, samples of three typical sandstones were heated to several temperatures up to a maximum of 800C and then allowed to cool to room temperature. Heating was done under both atmospheric pressure and simulated reservoir pressure conditions. Physical properties of the samples were measured before and after heating and comparisons made. Measured properties included permeability, sonic velocity, breaking strength and fracture index. Changes in physical properties were compared to changes in mineralogical characteristics as determined by thin-section. X-ray diffraction and chemical tests. EQUIPMENT AND PROCEDURE Two outcrop sandstones (Bandera and Berea) and one sub-surface sandstone (St. Peter) were selected for the tests. These samples have a wide ranee in composition and physical properties as shown in Table 1. The first series of tests was made on 2-in. diameter by 5-in. long test specimens. Test specimens used in all later work were 3/4-in. diameter by 1 1/8-in. long, this being the specimen size required for heating at simulated reservoir pressures. After careful washing, the cores were oven dried at 100 ± 5C for a minimum of 24 hours before the tests were run. Test specimens were heated in an electric furnace at a constant rate of temperature increase of 3C per minute. When maximum temperature of the run was reached, the sample was allowed to soak for one hour. The furnace was then cooled to room temperature at the same rate of 3C per minute. The entire heating operation was designed for reproducibility without subjecting the test specimens to excessive thermal shock. For samples heated under simulated reservoir pressures, a pressure cell designed by Dean was used (Fig. 1).3 The core sample was inserted into a thin-walled (0.006 to 0.01-in.) copper cup which was then mounted in a high-pressure cell. Provisions were made for the application of internal pore pressure as well as confining pressure. Tests showed that the thin-walled copper cup closed tightly around the core and satisfactorily transmitted confining pressure to the core. The core was heated by placing the entire cell into the electric furnace. The heating program was the same as that used in the atmospheric pressure runs: tempera-ture rise of 3C per minute to maximum temperature of the test, soaking at maximum temperature for one hour, and cooling at a rate of approximately 3C per minute. The cell was designed to withstand 5,000 psi at 1,000C. However, since it was considered likely that repeated heating and cooling would in time weaken the steel, 2,000 psi at 850C was set as a working limit. In the present series of tests, the pore pressure was held constant at 750 psi and the confining pressure at 1.500 psi. The pressure source was a high-pressure nitrogen tank. The two pressures were controlled manually but are accurate well within ± 50 psi.

Jan 1, 1966

By F. I. Stalkup

Vapor-liquid phase equilibrium experiments have been conducted in a static equilibrium cell on mixtures of a light, 45 API stock- tank gravity reservoir fluid and a rich hydrocarbon gas containing approximately 55 mole per cent of intermediate hydrocarbons. Both apressure-vs-composition study of the gas and a simulated reservoir fluid, and a multiple-batch contact simulation of the condensing-gas-drive oil recovery process were performed. In the latter experiments equilibrium gas and liquid compositions were analyzed. Also, conventional, "condensing-gas-drive", long-tube displacement experiments of the reservoir fluid and gases of various richness were performed. The results of these experiments could not be satisfactorily interpreted by the conventional pseudo-ternary-diagram representation of multicom-ponent phase behavior. The results seem to be explained better by considering a bubble-point surface and a dew - point surface joined in a plait-point locus. Portions of the plait-point locus cannot be "seen" directly by the rich hydrocarbon gas because of curvature of the dew-point surface. In such a system, continuous injection of the rich gas over stationary reservoir fluid might form a zone of contiguously miscible compositions from pure rich gas to pure reservoir fluid by: (1) saturating the reservoir fluid with injected gas to the bubble-point surface; (2) creating by mass transfer with fresh injected gas a path of contiguously miscible compositions along the bubble - point surface to the plait-point locus; and (3) creating by mass transfer with additional injected gas a path of gas compositions along the dew-point surface up to the point where direct miscibility results between dew-point fluid and the injected rich gas. INTRODUCTION The use of the pseudo-ternary-phase diagram to illustrate miscible displacement phase behavior has been discussed by several authors.5,7 Such a representation of phase behavior is not rigorous, but the ternary diagram nevertheless gives a qualitative picture of what actually occurs in a miscible displacement process. Fig. 1 is a typical illustration of miscible displacement phase behavior by a ternary diagram. The multicomponent hydrocarbon system is divided into three pseudo-components: a light fraction containing methane and nitrogen, an intermediate fraction containing ethane through hexanes plus carbon dioxide, and a heavy fraction containing heptane and heavier components. A two-phase region is bounded by a dew-point curve and a bubble-point curve, which are joined at the critical point. The concept deduced from such a representation for miscible displacement by a condensing-gas-drive process is as follows: a rich gas G, which lies to the right of the limiting tie line through the critical point C, is injected into the reservoir and contacts reservoir fluid L, saturating the reservoir fluid to give bubble-point fluid L1 and equilibrium dew-point gas GI. Continued injection of rich gas changes the composition of the saturated liquid L1 through a series of liquid compositions lying along the bubble-point curve, until the critical composition C is reached, at which point direct miscibility with the rich gas is achieved. Some equilibrium gas with compositions lying along the dew-point curve from G1 to C is also formed in this process. Thus, continuous injection of rich gas forms a

Jan 1, 1966

By G. H. Alani, H. T. Kennedy

One of the major difficulties in predicting the performance of oil reservoirs from their early pressure history lies in the uncertainty of estimating the volume of the liquid hydrocarbons contained in them. As a first step in filling this need, an equation was developed to determine the Inolal liquid volume of pure hydrocarbons over a wide range of temperature and pressure. The second step consisted of adapting the equation to apply to mixtures, with the heary hydrocarbons expressed as c The equatious are similar in form to van der Waals equation. hut the constants a and b tire considered as furnctions of temperature. In addition to the gas constant R, there are four constants characteristic of each hydro- When comptired with experimental values jound in the literal the. the average absolute deviation in the calculated molal volumes is found to be a Maximum of 0.33 per -cem for any of the pure liquid hydrocarbons studied. This maximum deviation was that found when comparing the calculated and observed values over a temperature range of 86° to 482°F and a pressure range from the bubbl-point ro 10,000 psis. The equations expressing the correlation for. mixrurrc were developed from 647 experimental measure-rilents of volume on 47 bouom-hole samples covering ; temperature range of 72" to 250°F and a pressure znge from bubble-point ro 5,000 psig. The average absotute deviation was found to he 1.6 per cent with the maximum for any measurement of 4.9 per cent. INTRODUCTION Accuratc information of the pressure-volume-temper-iiture behavior of hydrocarbon liquids is of considerable importance in the field of both applied and theoretical science and, especially, in the solution of petroleum reservoir engineering problems. These PVT relationships can be expressed graphically, in tabular form or as equations of state. Equations of state for gases have been developed based on the ratio of pressure, temperature and volume to the pressure, temperature and volume at the critical points of the hydrocarbons. These equations have been successful in yielding values fairly accurate for reservoir engineering purposes. For liquids, however, relations of this type are more complex and have not yielded useful results. The composition of reservoir liquids is normally expressed in terms of the moles of hydrocarbons C,? C,. C,. elc, per mole of mixture. The composition must include the fraction of a mixture such as C,,, which cannot be separated with reasonable time and effort into specific hydrocarbons. Although there are several equations of state having variable degrees of accuracy applicable to pure liquid hydrocarbons, none includes a means for estimating the volume of mixtures as a function of temperature. pressure and composition expressed in this manner.

By S. J. Prison, C. R. Smith

The effect of fluid injection to control water coning in oil and gas wells was investigated. Analytical and model techniques were employed. The factors investigated were the position and length of the completion interval, the point of fluid injection, the viscosity of the injected fluid and the relative thickness of the oil and water sections. The resulting influence of these factors on the net producing water-oil ratio was determined. Several important conclusions can be drawn from the study. In general, it was found that the net producing water-oil ratio can be reduced by fluid injection. The magnitude of this reduction depended on the factors listed above. An important practical consideration is that the injection fluid may be either oil or water. If the injected fluid is less dense than the connate water of the reservoir, the fluid will not be lost. This fact is reassuring when valuable oil is being injected. Efforts to suppress water production were more successful when the injection fluid was more viscous than the reservoir oil, or when a zone of reduced permeability existed in the vicinity of the point of fluid injection. Under test conditions, little benefit was derived through the use of impermeable barriers or cement "pancakes". INTRODUCTION The occurrence of water coning has been known for at least 60 years. In thin oil or gas pay sections, the presence of an oil-water or gas-water contact hinders production and often causes early abandonment of the afflicted well if a completion is even attempted. Even when relatively thick pay sections are found, the encroachment of water when a water drive is present will eventually pose serious water coning problems. This water is often corrosive, expensive to separate from the oil or gas and is costly to dispose of. The theory of water coning has been discussed by a number of authors. 1,2,3,4 Briefly, water coning to the producing interval in a well is due to pressure gradients resulting from the production of fluid from the reservoir. These pressure gradients will cause a water cone to rise toward the bottom of the producing interval if a water-oil or water-gas contact exists. The tendency of the water to cone is offset or partially offset by gravity forces since the water has a higher specific gravity than the oil. A balance then exists between two forces, gravitational forces arising from the difference in specific gravities of the oil and water, and the pressure gradients causing the flow of fluids to the wellbore. If the pressure gradient exceeds the gravitational force, water coning to the wellbore occurs and water production results. Through the years considerable thought has been given to the water coning problem. More than 50 U. S. Patents have been granted to inventors on the subject. A relatively complete literature review of the water coning problem has been made. 5 A number of these patents hold considerable promise for the solution or the partial solution of the water and/or the coning problem. Very little has been written describing field tests of techniques for the suppression of water coning. A notable exception is the paperby West. 6 He reports success in reducing gas coning by a combination of gravel packing and oil injection above the oil-producing interval. He also describes a comparable method to prevent water coning, but provides no field examples. This study experimentally and analytically verifies the benefits of oil injection as a means of partially or completely suppressing the water cone. While the gas coning problem was not treated, it is anticipated that results comparable to those obtained in water suppression could be obtained with reduced oil injection since the viscosity contrast between oil and gas exceeds that between oil and water. For the purpose of this study, the conventional potential flow theory was applied to the water-coning problem. The experimental verification centered on both a radial and a linear model. The model study permitted the investigation of complex flow configurations and the use of fluids of differing densities and viscosities. No analytic expressions are available to permit a solution of the problem as stated (see Fig. 1).

By D. C. Lindley, D. H. Gaucher

The displacement of oil by water in a waterflood project is accomplished by the action of transient viscous, gravitational and capillary forces which drive fluid through interconnecting pore spaces toward production wells. The relative importance of each of these forces to the effectiveness of the displacement and the production history depends on properties of the reservoir itself and the manner in which it is operated. There has been considerable discussion'-3 concerning the importance of these factors and the extent to which they control the performance of individual waterflood projects. No data have been presented, however, to demonstrate the effect on waterflood behavior of variations in rock permeability, water-injection rates and mobility ratios in three-dimensional reservoir systems in which viscous, gravitational and capillary forces are allowed to assume the proper relative influence on fluid flow. Much information relevant to water flooding has been gained from the study of field case histories. But, the complexity of the systems involved, the difficulty of defining geometry and obtaining sufficient data, and the fact that the same reservoir is never flooded twice with the same initial conditions, all obscure the role of process variables such as rate. Further, mathematical calculations of waterflood behavior are restricted at this time by available mathematical techniques to highly simplified situations. Scaled models, however, do offer an approach at the present time to evaluation of the effects of the pertinent variables. Several investigators have predicted prototype reservoir behavior from observations of model waterflood performance under a variety of conditions. Two-dimensional aspects of waterflood pattern efficiency have been investigated by Muskat, Aronofsky, Dyes, et al, raig, et al,' and Rapoport, et al. he effects of gravity segregation, viscous fingering and permeability stratification have also been studied in two-dimensional models. Craig, et al," studied the importance of gravitational forces on frontal displacement in three-dimensional models with wells on a five-spot pattern. This work in- vestigates the tendency for water to segregate in homogeneous reservoirs because of gravity and to under-run the oil for a range of water-inject ion rates. The effects of permeability stratifications were investigated by these authors using miscible fluids to simulate water-channeling through permeable strata in systems of negligible capillary force. In this investigation, a three-dimensional model of a five-spot pattern flood was designed to represent typical prototype conditions, and the sand and liquid properties were selected so that capillary, gravitational and viscous pressure gradients were scaled. It was not, however, possible to achieve exact similitude regarding viscous fingering when the more viscous oil was used. The effects of rock stratification on waterflood behavior were studied by changing both the permeability ratio and the position of the more permeable sand stratum in a two-layered, horizontal model. For each type of stratification, constant-rate floods were performed and the rates were varied, in separate experiments, over a tenfold range. To observe the effect of mobility ratio on flood performance, two series of constant-injection-rate floods were performed on the most highly stratified reservoir at a water-oil viscosity ratio of one.

By R. V. Higgins, A. J. Leighton

A fast method is needed for calculating thoroughly the performance of two-phase flobr~ in reservoir rock with complex geotrietry. The authors present such a method and .slzow its accuracy by romparison with laboratory performrrnce of a five-spot pattern. In making the performance c.alculations, the octant of the five-spot pattern was divided into four channels. The computer tinle required for the culculation was one rllinrcte. This is rorighly about minute per channel. INTRODUCTION By using an electronic computer it is now possible to make reservoir performance calculations in which fewer assumptions and more significant variables are included than was possible with a desk calculator. These performance calculations include such variables as the properties of the reservoir rock and its contained fluids, properties of the rock-fluid system and even different well-spacing patterns. Usually, the more variables that are included the higher the computer cost—and time on the computer is expensive. As a result, a short procedure is needed so the many variables will be included and accurate performance calculations will result. If a reservoir is composed of alternate layers of permeable and impervioi~s rock, and if each conducting layer has a different relative permeability, the combined computer cost for calculating each layer by the present long computer programs (which includes different well-spacing patterns and other variables) would be expensive. The performances of complex spacing patterns in a field often are studied by a potentiometric model, but in these model studies the analogy is to the flow of only a single phase. This paper presents a procedure to calculate the performance of a five-spot pattern with two phases flowing using a potentiometric model as a guide for the dimensions of channels and other related data. The time required to calculate the complete performance is about one minute on the IBM 7090 or four minutes on the IBM 704. The procedure can be extended to calculate the performance of more complex well spacings or an entire field with little modification. The procedure has been tested on the performance of the laboratory water floods of a five-spot pattern reported by Douglas, et alThe range of viscosity ratios in the laboratory experiments is from 0.083 to 754, which is extensive. The accuracy of the performance calculations is excellent as shown by the plotted points in Fig. I. The recovery as influenced by the geometry of the five-spot, taking into account the deviation from true radial character, has been predicted by Aronofsky and Ramey." The authors used a potentiometric model to determine the sweep efficiency to breakthrough of the flooding agent. Dyes, Caudle and Erickson,' by means of a laboratory model, show the influence of mobility ratio on recovery after breakthrough of the injected material. Dyes, et ai, used the X-ray shadowgraph technique and miscible phases in their studies. Craig, Geffen and Morse from the analysis of the performance of their laboratory model correlate many variables, including average saturation of a flooding agent determined by permeability relationship. Craig, et al, used X-ray shadowgraphs to determine the sweep efficiencies in sandstone models. In this report the sweep efficiency is not needed, and the important effect of continuously changing saturation gradients before and after breakthrough

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst's horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter's edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer