Search Documents

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

By J. D. Rice, P. L. Wearden, W. L. Hensou

Solutions to the unsteady-state partial water-drive reservoir performance problem can be obtained through the use of analogue computers or high-speed electronic digital computers. The solutions that have previously been resolved for use on digital computers, however, demand a knowledge of the aquifer parameters. Generally, the analogue computers now in use do not require this knowledge. A solution is presented herein where the aquifer performance, expressed in terms of difference equations, is related to the reservoir performance as expressed by the modified Scbiltbuis material-balance equation. A numerical procedure for a medium-sized digital computer also is presented in which a solution to the set of equations defining the aquifer and reservoir performance is obtained and the aquifer parameters (permeability and sand thickness) are automatically optimized while simultaneously matching the known pressure behavior. Predicted pressure behavior can be calculated using rates from any assumed future production practice. The procedure provides an output format which presents the cumulative, incremental and average rate of natural water influx, the per cent gas-cap expansion, the calculated reservoir pressure, the measured reservoir pressure, and the difference between the measured and calculated pressures. Results of a test problem are presented in comparison with results obtained by the Bruce Analyzer, Ohio Oil Co. 's Pace General Purpose Analogue Computer, and Sun Oil Co.'s Single-Pool Electronic Reservoir Analyzer. These results indicate that unsteady-state reservoir perlormance for a single-pool system can be adequately simulated by a numerical method employing a digital computer and that the special-purpose analogue computer can be supplanted by this method. INTRODUCTION The numerical approach to the unsteady-state partial water-drive reservoir performance problem heretofore has been limited in that solutions have been resolved only for certain fixed parameters. Also, prior to the advent of high-speed electronic digital computers in the last few years, calculations were so tedious and time-consuming that they ordinarily were neglected, or some other means of obtaining the solutions were used. The analogue approach had been used for obtaining solutions to the transmission-of-heat problem and, consequently, means were developed for applying the electrical analogue to the unsteady-state reservoir performance problem. Bruce submitted a paper for AIME publication in Sept. 1942, in which he described an electronic device (now called the Bruce Analyzer) for analyzing the water-drive performance of any reservoir, based on either steady or unsteady-state flow characteristics. 1 Sun Oil Co. adopted the analogue approach in 1949 and now has two special-purpose high-speed electronic reservoir analyzers in use.2,3 Two main features of Sun's analyzers are (1) the high repetitive rate of reproducing the entire history of reservoir performance cyclically a number of times per second, and (2) the elimination of the necessity of knowing the properties of the associated aquifer. Since the advent of larger and faster digital computers, the possibility of incorporating this second feature into a numerical approach appeared feasible. Previous work has been done and material published covering the application of numerical methods to the behavior of the pressure with respect to time at any point in an associated aquifer? However, the aquifer parameters must be known from some source. The material presented in Ref. 4 has been used to a great degree in developing the numerical method described herein. DIFFERENTIAL EQUATION AND BOUNDARY CONDITIONS The problem of simulating the performance of a reservoir and the associated aquifer is to determine the characteristics of the aquifer such that the

By J. P. Heller, H. B. Bradley, A. S. Odeh

A potentiometric model technique is presented for determining the areal sweep efficiency of a five-spot well pattern, at and beyond breakthrough. A sharp interface between displaced and displacing fluids is assumed. Although the prototype svstem is referred to as a waterflood operation, obvious changes in the notation and computations will adapt the results to other displacement processes. The results of the study include the following. 1. Areal sweep efficiencies for a five-spot well pattern, at and beyond breakthrough, for mobility ratios (displacing to displaced fluid) of 4:1, 2:1, 1:1 and 1:4. 2. Extension of potentiometric analysis to the investigation of the areal sweep efficiency beyond breakthrough in a five-spot pattern by the application of conformal mapping and conductive-solid models. 3. The use of layers of conductive fabric in representing mobility ratio changes in potentiometric models. 4. The development of a probe mechanism for probing conductive solids. The results obtained by conductive-cloth models agree with earlier areal sweep efficiencies at breakthrough obtained by Aronofsky and Ramey 1 on the potentiometric analyzer using electrolytic-tank models. Results beyond breakthrough differ from those obtained by the X-Ray Shadowgraph technique. Data from this study show that, for mobility ratios greater than one, water cut rises rapidly as fluid is produced after breakthrough. However, for mobility ratios smaller than one, a large increase in area swept resulted with only a small increase in water cut. INTRODUCTION In calculating reservoir performance of waterflood-ing operations and other fluid-injection programs, it is necesssary to estimate the areal sweep efficiency * before and after injected-fluid breakthrough into production wells. Influence of mobility ratio on oil-production history, before and after breakthrough for a five-spot well pattern, has been studied by X-Ray Shadowgraph techniques 2-5 and by gelatin models.6 In none of these investigations was the transition zone controlled experimentally. Steady-state anaiog techniques assume that a vertical and discrete interface exists between the displacing and displaced phase. Thus, when using a potentiometric analysis, it is assumed that no transition zone exists. To determine the areal sweep efficiency beyond breakthrough for a five-spot well pattern using the potentiometric analyzer, four principal problems had to be resolved. These were the following. 1. The laborious method of representing changes in mobility ratio in electrolytic tanks by means of contoured wooden blocks. 2. Vapcrization of the electrolyte solution during the study causing concentration changes and, thereby, variable conductivity. 3. The scale-limitation problem 7 that is accentuated by abrupt eiectrolyte depth changes. 4. A wellbore geometry problem. EXPERIhlENTAL EQUIPMENT AND PROCEDURES CONDUCTIVE-SOLID MODEL The theory, equipment and procedural techniques of the electrolytic-tank potentiometric model have been treaced excellently by Lee8 and by Muskat.9 By the use of a conductive solid to represent the fluid conductivity of the porous media, the need for contouring wooden blocks to represent mobility ratio was eliminated. It was found that a carbon-black coated fabric, Uskon D-16,** could be employed as a fairly uniform conductive media. This material, when stacked in appropriate layers and pressed together, can be used to represent changes in mo-

By D. H. Henley, F. F. Craig, W. W. Owens

The oil recovery performance of systems producing entirely by bottom-water encroachment has been experimentally determined in a series of scaled laboratory-model tests. The effects of well spacing, fluid mobilities, rate of production, capillary and gravity forces, well penetration and well completion techniques on the oil recovery performance have been investigated. The laboratory tests were performed using two uniform, un-consolidated sand-pack models. The models have ratios of the interwell distance to the formation thickness of 12 and 2, respectively. Tests at constant total fluid production rate were performed simulating a range of uniform reservoir characteristics and operating conditions encountered in field operations. The performance was determined by material balance and by observation of the encroachment of dyed fluids into the models. The results of the model tests agreed with those obtained mathematically when the conditions previously considered in theoretical studies were simulated, that is, when the oil and water are of equal density and no capillary forces exist. The model study of bottom-water drive indicated that certain variables can aoect the oil recovery performance to a greater degree than can be predicted by present analytical methods. In one comparison, the oil recovery at a water-oil ratio of 20 (obtained at a wide well spacing) varied as much as threefold, depending upon the system's properties and the production rate. Lesser effect of mobility ratio and no eflect of capillary forces over the range studied were observed. The test 'results also showed that the deeper the well penetration into the oil column, the greater the total water production to a producing WOR of 20. However, the ultimate sweep efficiency, and so the oil recovery to this level of WOR, did not vary significantly with well penetration. Horizontal fractures at the top of the formation did not significantly change the sweep characteristics of the reservoir models when values of radius and fracture capacity encountered in actual reservoirs were used. Impermeable pancakes at the bottom of the well moderately increased the oil recovery efficiency both at water breakthrough and at high water-oil ratios. A method is outlined by which the oil recovery performance of other uniform bottom-water drive systems can be estimated from the information obtained in these model tests. INTRODUCTION When oil is produced from a well which partially penetrates an oil zone completely underlain by water, the water rises directly beneath the well in a symmetrical cone when the system is uniform. Two different flow mechanisms can cause the water cone to form—coning and bottom-water drive. In coning, the aquifer is relatively inactive and the cone is formed beneath the well by the pressure gradients associated with the oil flow to the well. The oil can be produced by a solution-gas drive, an edge-water drive or other driving forces in the interwell area. In a bottom-water drive, the driving force for oil production comes from an upward encroachment of the underlying active aquifer. Two papers have analyzed the theoretical performance characteristics of bottom-water drive reservoirs. In the initial mathematical investigation, Muskat' established the equations which determine the pressure distribution in this type of reservoir and solved these equations for certain conditions. Specifically, it was assumed that the water and oil had equal mobilities and equal densities, there were no capillary forces, the pressure throughout the oil zone remained above the bubble-point pressure, a constant pressure existed at the initial water-oil contact and the oil was completely displaced by the encroaching water. These assumptions were used in obtaining analytical solutions. In general, Muskat found that the sweep efficiency to initial water breakthrough to the well was larger for the thicker oil zones, the closer well spacings, the lower ratios of vertical to horizontal permeabilities, the smaller the penetration of the well into the oil zone and the smaller the bore size of the well. The production history after water breakthrough was expressed as a volumetric sweep efficiency at a given producing water-oil ratio. The results indicated that cumulative oil production at producing water-oil ratios of 10 is less affected by the well spacing than is the water-free production history. Muskat studied well spacing which today would be regarded as close. The maximum value of his dimensionless well spacing (ratio of interwell distance to formation thickness) was 4.3. This would require the development of a 50-ft-thick oil sand on less than 10-acre spacing if the vertical and horizontal permeabilities were equal, with

By F. Brons, W. C. Miller

Pressure information for use in material-balance calculations is obtained, where possible, from pressure build-up surveys in shut-in wells. Using proper extrapolation methods, static pressures are obtained which, by averaging, give the field static pressure. Often, particularly in fields with a large number of wells, only spot pressure readings are available. A graphical method is presented here which enables the determination of the corrections which must be applied to the spot readings to give an estinzate of the static pressure. Although the static pressures thus obtained are not so reliable as those from build-up surveys, they are more nearly correct than the spot readings themselves for use in material-balance calculations. INTRODUCTION The average static pressure in a bounded reservoir is determined from a number of individual-well pressure measurements. Where pressure build-up curves are available, extrapolation methods such as described in Ref. 3 may be used to determine the static pressure in the drainage area of each individual well. The possibility of using pressure build-up curves presupposes that each well being surveyed is closed in during a certain period and that, during this time, a pressure recorder is left suspended in the hole. For surveys in fields with a large number of wells, this is usually not feasible for economic and other reasons. Either a large number of pressure recorders would have to be used simultaneously or, using a limited number of recorders, the time interval needed to complete a field-wide survey would be too long. The data obtained for this incomplete coverage are insufficient to permit reliable averaging of the few calculated static pressures over the entire reservoir. In many such cases another type of pressure information is available, i.e., spot pressure readings. During this type of pressure survey all wells are closed-in simultaneously and, after a reasonable shut-in time (e.g., three days), pressure readings are taken by running the recorder in and out of each successive well. These surveys, which are less expensive, are often taken at regular intervals, say once or twice a year. In fact, such surveys are required by the state regulatory bodies in many portions of the united States. These spot readings are normally lower than the static (fully built-up) pressure due to the relatively short shut-in time. No direct means of extrapolating such a reading to the static pressure is available because no pressure build-up record is obtainable. Consequently, these data are inferior to pressure build-up information as a basis for estimating the stabilized reservoir pressure in the drainage area of a specific well. For the common case where abundant spot readings but only a few pressure build-up surveys are available, our studies suggest that it may be possible to design a program which will combine the two types of information and provide an estimate of the average reservoir pressure that is more reliable than that obtainable solely from either set of data. in time) gives a point marked p. The actual static pressure at the time the well was shut-in is marked by p.. The illustration shows that, even if a build-up curve is available and can be extrapolated, a correction

By J. D. Testerman

A statistical technique to identify and describe naturally occurring zones in a reservoir and to correlate these zones from well to well is described. The technique is particularly useful in describing a reservoir where cross flow between adjacent strata is important in determining reservoir behavior. Although it has been used primarily for permeability zonation, the technique is general and can be used to correlate any reservoir property or related data, such as the information contained in well logs. INTRODUCTION One of the first problems encountered by the reservoir engineer in predicting or interpreting fluid displacement behavior during secondary recovery processes is that of organizing and using the large amount of data available from core analysis. Permeabilities pose particular problems in organization because they usually vary by more than an order of magnitude between different strata. Due to the sheer volume, it is almost always necessary to group data and to use an average value to represent a number of measurements. Perhaps the most common method now used to group permeability data is the capacity-fraction technique, which ranks permeabilities in order of magnitude, regardless of the physical location of the permeabilities within the reservoir. The cumulative per cent capacity is plotted against cumulative per cent thickness. This plot is divided into an arbitrary number of zones, generally of equal thickness. Five zones (or averaged groups of data) usually are obtained, each of which is treated as homogeneous in subsequent calculations. The division so obtained has no physical meaning; strata in the same zone, calculation-wise, are usually not adjacent in the reservoir. Reservoir engineering techniques being developed will handle crossflow that occurs between adjacent communicating reservoir strata because of imbibition and gravity segregation. Since crossflow occurs between physically adjacent layers within the reservoir, a new zonation technique recognizing the actual location of strata within the reservoir is necessary. Similarly, the recognition of natural zones is important for predictions of oil recovery by processes involving diffusion. One such process is miscible displacement, where predictions of lateral diffusion within the reservoir must recognize the actual location of the invaded zones in relation to the rest of the formation. Natural zones must also be adequately recognized to account for heat transfer within the reservoir during thermal exploitation. Because of the complexity of the problem, statistics appear to offer the only practical hope of dividing a reservoir into physically-meaningful natural zones. This paper presents a statistical technique for identifying these natural zones and for ascertaining which ones are likely to be continuous between adjacent wells. The zones defined have minimum variation of permeability internally and a maximum variation between zones. The technique is general and can thus be applied to reservoir properties other than permeability. The method will guide the reservoir engineer in estimating which zones are likely to be continuous between wells. However, a statistical correlation based on permeabilities in two different wells is no guarantee that the zones so defined are, in fact, continuous. Rather, the assumption of continuity must be consistent with geological data concerning the depositional environment, as well as justified on the basis of engineering judgment in combination with statistics, just as judgment is required with conventional zonation methods. CALCULATION PROCEDURE The reservoir zonation technique is a two-step operation. The steps are individually described, and a sample calculation is presented in the Appendix. ZONATION OF INDIVIDUAL WELLS First, the set of permeability data at a single well is zone, into Zones. These zones are selected so that variation is minimized within the zones and maximized between the zones. The equations4,6 used to zone the data are where B = the variance between zones, , = the number of zones, i = the summation index for the number of zones, j = the Sumation index for the number of data within the zone, mi = the number of data in the ith zone, k,. = the mean of the permeability data in the ith zone, k . = the over-all mean of the data in the well, W = the pooled variance within zones, N = the total number of