Search Documents

By A. S. Michaels, M. C. Porter, A. Stancell

Previous laboratory studies have demonstrated that the injection of small quantities of reverse wetting agents during water displacement can increase oil recovery from unconsolidated porous media.9 In the present investigation, an attempt has been made to determine more fully the effects of reverse wetting treatments and to clarify the mechanism by which increased oil recovery is effected. Water-oil displacements were performed in beds of 140-200 mesh silica sand. Hexylamine slugs (injected after 0.25 pore volume of water throughput), when adequate in size and concentration, were effective in promoting additional oil recovery. Their effectiveness increased with the quantity of amine injected. However, slugs of sufficient size and concentration to stimulate oil production at water flow rates of 34 ft/day did not do so at 4 ft/day. Visual studies in a glass grid micromodel have shown that the stimulation of oil production, via ayueous hexylamine, is a result of transient changes in the oil wettability of the pore walls. if the amine slug is of sufficient size and concentration to induce significant changes in the adhesion-tension, large continuous oil masses will be formed. if the superficial water velocity is high enough to result in rapid desorption of the amine, a favorable "wettobility gradient" may be established across the masses; under such conditions, high oil mobility is observed, and increased oil recovery results. INTRODUCTION It is generally agreed that the efficiency of oil displacement by water in porous media is limited in part by capillary forces which cause the retention of isolated masses of oil — resulting in the so-called "irreducible minimum oil saturation". Recent estimates10 indicate that there are about 220 billion bbl of petroleum in United States reservoirs which are not economically recoverable with present techniques (such as water flooding). This amounts to almost five times the known recoverable reserves. It has been recognized for some time that a suitable alteration in the water-oil interfacial tension and/or the contact angle, as measured between the water-oil interface and the solid surface, should result in better displacement efficiency. Surface active agents can be used as interfacial tension depressants to accomplish this objective,1,2,7 but unfortunately, the additional oil recovery is seldom commensurate with the treatment cost. In contrast to interfacial tension depressants, the effect of contact angle alterations on water-oil displacements has received relatively little attention in the literature. It is known that the wettability affects the displacement process.3,4,11 Displacements in water-wet systems generally result in lower residual oil saturations than those in oil-wet systems. The effect of "'transient" wettability alterations concurrent with the displacement process have been investigated by Wagner, Leach and coworkers,4,11 wherein it has been demonstrated that the establishment of water-wet conditions during water flooding of oil-wet, oil-saturated porous media is accompanied by significant increase in oil displacement efficiency. Michaels and Timminss studied the effects of transient contact angle alterations resulting from chromatographic transport of reverse wetting agents through unconsolidated sand. It was demonstrated that chromatographic transport of short-chain (C4 through C8) primary aliphatic amines can improve oil recovery and that the recovery increases with the quantity of amine injected (i.e., with either the amine concentration or the volume of the slug injected). Circumstantial evidence indicated that the increased displacement efficiency resulted primarily from transient changes in wettability of the porous medium. In the Present investigation. additional information has been obtained on the effects of reverse wetting treatments and the mechanism by which increased oil recovery is accomplished. Hexylamine slugs of varying size and concentration were chromato-

Jan 1, 1965

By G. C. Bernard, L. W. Bernard, L. W. Holm, W. L. Jacobs

The effect of loam on the permeability of porous media to water was studied as a {unction of foaming agent concentration, specific permeability, pressure gradient, length of a porous medium and its oil saturation. At a given fluid saturation in a porous medium, the permeability to water was found to be the same whether loam was present or not. Foam decreases the permeability to water by developing a higher trapped gas saturation than that obtained by water flooding without loam present, increasing the concentration of foaming agent increased the trapped gas saturation and thereby decreased the permeability to water. The presence of oil reduced the capability of most foaming agents to decrease the permeability of a porous medium to water. A few surfactants were found to be effective loaming agents even in the presence of oil. These results are similar to those reported in a previous paper on the effect of loam on the permeability of porous media to gas. The effect of loam was found to persist in long porous media at moderately high reservoir temperatures and during the passage of many pore volumes of surfactant-free water. INTRODUCTION This paper describes part of a study on a novel approach in the use of surfactants for oil recovery; the use of foam rather than water to displace oil. Previously it was found that foam can displace oil which normally is not displaced by water.4 The foam is formed by successively injecting a suitable surfactant solution and gas into a porous medium. Foam appears to have at least two uses in the field: (1) it shows promise as a superior oil recovery agent, and (2) it shows promise as a selective permeability reducing agent. Foam may be very useful in water floods, or in other oil recovery processes, where highly permeable streaks or unfavorable mobility ratios are a problem. A previous paper reported the effect of foam on the permeability of porous media to gas.1 In the present study the effect of foam on the permeability of porous media to water is reported. The specific objectives of the study were to determine: (1) the effect of foam on the permeability to water in porous media of various specific permeabilities, (2) the effect of foam on the permeability to water in the presence of oil, (3) the effect of foam and crude oil on the trapped gas saturation, (4) the effect of foam on permeability to water at trapped gas saturation, (5) the effect of pressure gradient on the permeability to water under foaming conditions, (6) the persistence of foam during the passage of surfactant-free water through the porous medium, and (7) the effect of various foaming agents, length of the porous medium and temperature on the permeability reduction caused by foam. EXPERIMENTAL PROCEDURES EQUIPMENT AND MATERIALS The experimental apparatus consisted of consolidated and unco n soli dated porous media, wet test meters and constant delivery pumps. The porous media consisted of consolidated sandstone cores (6 to 36 in. long), and unconsolidated sand packs (3 to 30 ft long). The consolidated cores had permeabilities of 32 and 1,000 md and porosities of about 20 per cent. The sand packs had permeabilities of 3,500 to 211,000 md and porosities of about 40 per cent. (Throughout this report a term such as 100 md sand" is used. This term means that the porous medium had a dry, nitrogen permeability of 100 md.) Fluids used in the experiments were distilled water, 1 per cent NaCl solution, aqueous solutions of foaming agents, nitrogen gas, air and crude oil.

Jan 1, 1966

By E. P. Miesch, P. B. Crawford

A study was made of the effect of permeable and impermeable lenses in a reservoir on the production capacity of a well. Both steady-state and unsteady-state data were obtained. An electrical resistivity model was used to obtain the steady-state data and thermal models. were constructed to obtain the unsteady-state data. The productivity of a well is affected very greatly only when the lenses are close to the well. The effect of circular lenses on the productivity ratio can be correlated with the distance from the center of the lens to the center of the well divided by the radius of the lens. When this dimensionless distance is equal to six or greater, the effect of the lenses on production capacity will be negligible. The pseudo steady-state productivity of a heterogeneous reservoir can be predicted using steady-state data. INTRODUCTION Many analytical solutions of reservoir behavior assume that reservoir rock is uniform and homogeneous. Although this assumption is used, all of the data from core analyses and well logging indicate that the reservoirs are heterogeneous. Very little work has been done on the performance of heterogeneous reservoirs. The work of Landrum,l et al. showed that transient phenomena in oil reservoirs could be studied with thermal models. Pickering and Cotman2 used thermal models to study flow in stratified reservoirs and investigated the effect of inhomo-geneities in oil reservoirs on transient flow performance. Loucks3 made a mathematical study of the pressure build-up in a system composed of two concentric regions of different permeability. Root, Silberberg and Pirson4 studied the effect of the growth of the flooded region on water influx predictions using a thermal model consisting of three concentric cylindrical regions of different thermal properties which simulated the aquifer, the flooded region and the unflooded portion of the original hydrocarbon region. Tomme, et a1.5 made a mathematical study of vertical fractures. The object of this investigation was to study the effect of highly permeable and impermeable lenses in the vicinity of the wellbore on the pressure depletion history of the well. Steady-state data were obtained for both conductive and nonconductive lenses that completely penetrated the formation. The lenses were symmetrically located at various distances from the wellbore. The unsteady-state data were obtained on seven thermal models. EXPERIMENTAL EQUIPMENT AND PROCEDURE STEADY-STATE DATA The steady-state data were obtained from an electrical resistivity model 30 in. in diameter and approximately 11/2 in. deep. The outside of the model was lined with a 30-in. diameter copper strip, which served as the outer boundary of the reservoir. The bottom was covered with a sheet of plexiglass so that it would be nonconductive. The model was filled with a slightly saline solution. The well size was varied from an 0.064-in. diameter copper wire to a 10-in. diameter copper cylinder. Readings were taken with an impedance bridge using AC current to prevent polarization at the contacts. Copper and wax lenses were used to represent infinitely conductive and nonconductive lenses, respectively. The resistance was first measured for each well diameter with no lenses in the reservoir. Then the conductive and nonconductive lenses were spaced symmetrically at various distances from the well and the resistance read from each lens location. The diameters of the conductive lenses were 3, 1.022 and 0.624 in., and those of the nonconductive lenses were 3, 2.25 and 1.563 in. The ratio of each resistance reading to the resistance reading with no lenses in the reservoir

Jan 1, 1965

By D. M. Waldorf

Steam permeability measurements have been made in the laboratory on several samples of natural reservoir materials. The steam temperatures and pressures were selected to simulate conditions which might exist in a reservoir during the injection of steam. For each sample tested, the experimental permeability to superheated steam was comparable to that measured with air and no evidence of plugging was detected. Some samples were exposed to water at various temperatures and plugging was found to occur in materials which contained significant quantities of monmorillonite clay. Temperature had little effect on the degree of plug-ning between 75 and 325 F. The measured pemeabilities tended to increase slightly with temperature, but the changes were small compared with the initial loss of per~neability on wetting. Sequential pemzeability measurements were made on two samples using air, water, steam, water and air, in that order. Both samples were water-sensitive and plugged extensively after the initial injection of water. Upon exposure to superheated steatm the samples dehydrated and their permenbilities to superheated steam were comparable to those initially measured with air. The remaining measuretnetzts with water and air confirmed that the water plugging was reversible and that the samples were not seriorrsly damaged during the tests. INTRODUCTION The swelling of water-sensitive clays during water floods has long been recognized as a potential source of reservoir damage. The recent extensive application of steam injection and stimulation has compounded this problem since both hot water and steam (as well as fresh water at reservoir temperatures) are, at sume time, in contact with the producing zone adjacent to the bore of a steam injection well. The purpose of this paper is to present data which compare the sensitivity of some natural sedimentary rock samples to water at various temperatures, and to super-heated steam. Some properties of montmorillonite clay are briefly reviewed, and comparisons are drawn between empirical data and the predicted behavior of the montmorillonite known to be present in the samples. PROPERTIES OF MONTMORILLONIT E CLAY Water initially adsorbs on dry N a -montmorillonite clay in discrete layers in the interlaminar space between clal platelets. The platelet spacing, which is 9.6 A (angstroms) for a dehydrated clay, has been observed to expand in discrete steps to 12.4, 15.5, 18.4 and 21.4 A spacings, indicating the formation of four discrete layers of regularly oriented water molecules.' The first two layers are easily formed by hydrating a dry sample to equilibrium in an atmosphere with carefully controlled humidity. The formation of the higher layers is more difficult. The usual X-ray diffraction patterns of the more highly hydrated samples indicate a gradual increase in the average spacing betwcen 15.5 and 19.2 A, followed by a discontinuous expansion to 31 A when the weight ratio of water to dry clay is between 0.5 and 1.2.' Platelet expansion above 31 A proceeds monotonically as the moisture is increased and no regular arrangement of the platelets ib observed. Water-sensitivity in sedimentary rocks is usually associated with Na-montmorillonite clay when it is in the noncrystal-line state. Mering3 found that the average lattice spacing of sodium montmorillonite hydrated at 68 F and 70 per cent relative humidity was 15.5 A, and that the spacing, at 92 per cent humidity was 16.5 A. The water adsorbed at the higher humidity has the same free energy as liquid water at 65.6 F. Kolaian and Low' used a tensiometer to measure the thermodynamic properties of water in diffuse suspensions of montmorillonite clays relative to pure water. They observed that water in suspensions as dilute as 6 per cent clay became partially oriented when left undisturbed. The bonding associated with this orientation was not extensive because the free energy difference between the water in suspension and pure water was only a few millicalories per mole. They also found that the measured free energy difference decreased rapidly with temperature and became negligible above 100 F. This evidence indicates that montmorillonites contained in sedimentary rocks would dehydrate to a crystalline structure when exposed to superheated steam, and that the rock permeability measured with steam would be equivalent to that measured with air. The effect of elevated temperatures on the swelline of montmorillonite clays in aqueous suspensions has not been investigated. The Gouy-Chapman diffuse-ion-layer theory has been used to predict the swelling pressure of clay suspensions in dilute salt solutions at room temperature with reasonable success. theory also correctly predicts the direction of the thermal response of Na-mont-morillonite swelling pressures in dilute salt suspensions, 9 Over the temperature range of 33 to 68 F, an increase in

Jan 1, 1966

By R. A. Thomas, R. L. Solbod

The importance of transverse diffusion on the finger development in a miscible-phase displacement at an adverse mobility ratio of tbree was studied in a porous plate 1/4-in. thick, 3-in. wide and 18-in. long. Fast displacement rates (29 ft/D) and slow rates (1.6 ft/D) were used to determine the effect of residence time on the geometry of the fingers. The shape of the fingers was observed directly by use of the X-ray technique. At fast rates numerous narrow fingers were observed, but at slow rates a single somewhat bulging finger was produced. The amount of material moved transversely by diffusion across the plate was sufficient to modify the finger geometry in the slow-rate run because of the long residence time. These results are in contradiction to some of the postulates in the literature. The composition of the effluent stream, however, was not affected by the flow rate. This result is not inconsistent with the observed change in the shape of the finger in a short model, but it seems likely that a short model does not offer adequate and proper scaling of the reservoir. The model used was probably a valid one for studying the effect of transverse diffusion on the finger geometry, but a longer model would be needed for proper scaling of the effect of the change in the finger shape on the efficiency of displacement as measured by the composition of the effluent stream. INTRODUCTION Fingering can be defined as the uneven advance of the injected phase as it moves into a porous medium displacing the resident phase from the pore spaces of the rock. The use of this term is usually restticted to the situation in which the displacing phase is less viscous or more mobile than the fluid being displaced. Under these conditions, not only are fingers formed, but the length and width of the fingers grow with distance traveled in the porous medium. This subject has become one of great interest to the oil industry because of the present trend toward the use of various forms of miscible-phase displacement to increase oil recovery. Since in nearly all of the known modifications of the miscible-phase displacements an unfavorable mobility ratio exists (the displacing phase has a lower viscosity than that of the crude oil), the conditions are proper for fingering to develop. An appreciable amount of fingering appears to be a severe handicap to these processes for it increases the volume of agent required for the process to be a success, and such an increase puts a severe strain on the economics of the proposed processes. In some cases, such as for a mobility ratio of 200 unfavorable, it has already been demonstrated that the proposed process would not be economic if the fingering in the field were to be of the same magnitude as that observed in the laboratory. A number of aspects of fingering have been studied and reported in the literature. While the phenomenon of fingering cannot be regarded as a completely understood subject, considerable information exists on the effect of the path length and the mobility ratio on the growth of fingers. Less-complete data are available on the effect of the diameter of the flow path on the character and amount of the fingering, and even less agreement in results exists on the effect of rate of flow on the nature of fingering. This paper deals with one aspect of this latter subject. OBJECTIVE The objective of this study was narrowed down to one rather specific feature of the behavior of fingers in miscible-phase displacement in porous media. The variable studied was the effect of rate of flow on the nature and the development of fingers. It should be made clear at this point that, while rate was the apparent variable, the real variable was residence time; that is, at low rates the fluids are present at a given spot in the porous medium for a longer time interval than at fast rates. The purpose of the study, therefore, was to determine the changes which occur in the fingering and the possible benefits which might accrue from a longer residence time during that period when fingers are

By H. O. Jahns, R. A. Greenkorn, R. E. Haring, L. K. Shallenbarger

This paper is a study of the effects of heterogeneity on flow in an analog of porous media, the Hele-Shaw model. l-3 A set of experiments in heterogeneous Hele-Shaw models showed streamlines through and around heterogeneities of various sizes, shapes and levels. (A level, we define as the ratio between the transmissibility of the heterogeneity and that of the rest of the model.) The heterogeneities were restrictions or expansions of the flow stream analogous to variations in the transmissibility of porous media. The experimental data agreed well with numerical results and with an analytical solution, which we derived for a circular heterogeneity in an infinite field. This study considers the flow-stream distortion due to the shape, size and level of heterogeneities. Size and level are much more important than shape provided the heterogeneity is not long and narrow. Our analytical solution shows that a circular heterogeneity in a large field can be replaced by an equivalent circle of either zero or infinite permeability. The radius of the. equivalent circle is a simple function of size and level of the actual circle. INTRODUCTION With the availability of high-speed computers and numerical procedures to predict reservoir behavior, we are faced with an important question. How much information about the reservoir do we need to justify the cost of the computer in any given case? To answer this question, we have to know how various reservoir parameters affect flow behavior. Reservoir heterogeneity is one of these parameters. In this study, we used a simple, two-dimensional model of porous media, the hele-shaw model, to investigate the effect of heterogeneity on flow behavior. We restricted ourselves to linear, single-phase, steady-state flow in a rectangular field with a single heterogeneity at its center. The analogy between flow in Hele-Shaw models and in porous media is easily verified. Let us first consider a homogeneous Hele-Shaw model with constant plate separation h. (The Hele-Shaw model is constructed by placing two plates, usually glass, very close together and allowing liquid to flow between them. Streamlines are made visible by introducing colored fluid into the space between the plates at a number of points across the model.) We assume a cartesian coordinate system with its origin in the middle between the plates and the z axis directed perpendicular to the glass plates (Fig. 1). The fluid flow is always in a direction parallel to the glass plates and varies from a maximum value to zero in the very small distance from the middle (z = 0) to either plate (z = h/2). For slow motion of an incompressible fluid, neglecting inertia and body forces, we have the viscous flow equation: where p is the fluid pressure, p the viscosity, and v the velocity vector with components u, v and w in the x, y and z directions, respectively. In our rase. and the derivatives of 3 with respect to x and y are small as compared with the derivative in the z direction. Therefore, approximately, with Vp and v being two-dimensional vectors in the x-y plane. Two integrations with respect to z give

Jan 1, 1965

By D. J. Pye

Certain high molecular weight synthetic polymers in very dilute solutions decrease water mobility in porous media 5 to 20 times more than would be expected from the solution viscosity. This indicates that increasing water viscosity to reduce adverse mobility ratios in water floods is economically feasible in many situations. Laboratory results are consistent with accepted theory with respect to improved areal sweep, permeability distribution and displacement efficiencies as would be predicted from the solution mobility measured in the core. The laboratory results have been confirmed by field pilot floods which recovered 80 per cent and 100 per cent more oil than comparative water-flood pilots in fields containing 16 cp and 130 cp oil respectively. Economic analyses based on pilot data indicate the method to be profitable. INTRODUCTION The desirability of improving the water-oil mobility ratio in waterflood-ing operations has long been appreciated, especially for high viscosity crudes. Methods of reducing the oil viscosity by application of heat, gas resaturation and miscible drives have met with some success. Conversely, it has been generally conceded that increasing the water viscosity would also be effective, but the amount of thickening agent required to effect an appreciable increase in viscosity is discouragingly high. The use of such materials as glycerin, sugar, or glycols is completely out of the question economically, so that considerable attention has been given to the use of the much more efficient synthetic water-soluble polymers as viscosity improvers; but even these materials are not economically attractive. If increased water viscosity could be economically realized, marked improvements in areal sweep efficiency, as discussed by Caudle and Wittel and others, would be realized in field operations. Similarly, the improved mobility ratio would bring about increased recoveries by correcting permeability distribution problems as discussed by Stiles2 and Dykstra and Parsons,' and improvements in displacement efficiency would be obtained as recognized by Buckley and Leverett4 and Welge.5 LABORATORY STUDIES THE RESISTANCE EFFECT In the course of research to develop a more efficient water-soluble polymer viscosity builder, we discovered that a very few types of water-soluble polymers exhibited a most unusual and very interesting property not previously observed. Among the earliest materials showing such activity were polymers containing acry-lamide. Fig. 1 shows a typical plot of the solution viscosity vs concentration for one of these polymer solutions as determined in an Ostwald viscometer. As usual, such a plot is an approximate straight line on semi-logarithmic paper. Now, an inspection of Darcy's law, shows that a viscosity value can also be determined in a formation sample if the pressure, flow rate and permeability are known. If the core permeability is previously determined with wa- ter or brine, then the viscosity of another aqueous fluid may be determined from the equation. The viscosities so determined for glycerin (glycol or polyvinyl alcohol solutions, for example) agree with the viscometer determinations. However, in the case of the water-soluble polymer solutions in which we are interested, the viscosities measured in the formation sample depart very markedly from the viscometer values. A typical relationship for a polymer solution is also shown in Fig. 1. This unusual departure from the expected response is referred to as a resistance property of the polymer and is quantified as "resistance factor". Since the resistance factor R can best be determined in a core under study, the term has been defined on the basis of the ratio of the brine mobility to the polymer solution mobility at residual oil saturation so that:

Jan 1, 1965

By F. O. Jones

The capabilities of small proportions of divalent cations, such as calcium or magnesium, for controlling clay blocking are reported. Potentially sensitive formations can be exposed to fresh water if at least one-tenth of the salts dissolved in both the native water and invading fresh water are calcium and magnesium salts. Often, even less suffices. The behavior evidently depends upon the cation exchange properties of the clays, which favor adsorption of calcium and magnesium over sodium. Clays having sufficient adsorbed calcium and magnesium resist dispersion by water and consequent blocking of permeability. An allied phenomenon, the dependency of clay blocking on salinity contrast, is also reported. Abrupt change from highly saline to fresh water can cause blocking which may not occur if salinity is lowered slowly. A process related to osmosis is thought responsible for this behavior. Applications of findings to drilling, fracturing and water flooding are discussed. INTRODUCTION Many formations, when exposed to fresh water, can lose permeability through clay effects.',' For example, drilling-mud filtrate can cause oil permeability decreases, which persist long enough to interfere with drill-stem testing and well completion.3 Ultimate productivity may suffer when clay damage is severe, especially if drawdown pressures are low.' The injection pressures and the time required for water flooding can increase if clay blocking occurs. Earlier work indicates that clay blocking is caused by obstruction of flow channels by clays or other mineral fines dispersed by fresh water.'!" It is well known that increasing salt concentrations in water tends to prevent clay blocking. Also it is well established, although perhaps not as widely known, that the nature of the dissolved solids is also important. If, for instance, formation clays are exposed to calcium chloride solution and then exposed to weaker solutions of calcium chloride or distilled water, considerably less permeability damage results than if the clays had been exposed instead to sodium chloride solutions. The reasons for this behavior are: (1) clays are cation-exchange, or base-exchange, materials similar in this respect to zeolites or exchange resins, and (2) clays in the calcium form do not disperse easily in fresh water, whereas clays in the sodium form do. The base-exchange form of a clay is easily altered by flowing a solution through it. For example, a sodium clay can be changed to a calcium clay simply by passing a calcium chloride solution through the clay bed. There is general agreement that the reason calcium and magnesium clays do not disperse easily in fresh water is because the calcium ions do not ionize easily from the clay surfaces.' Sodium clays, however, are believed to allow the sodium ions to ionize from the clay particle surfaces. In fresh water, this allows the particle to accumulate a net negative electrical charge great enough for the particles to repel one another and, therefore, disperse more easily. The foregoing suggests that, even in the presence of a large excess of sodium, a relatively small ratio of divalent cation, such as calcium or magnesium, might serve to prevent clay dispersion. The divalent ion probably adsorbs to a disproportionately large extent, because of being adsorbed more tightly, with the result that the proportion of divalent cation on the exchange positions of the clay might then be great enough for the clays to resist dispersion. Consequently, an investigation was conducted to find how mixtures of salts in solution influence clay blocking. APPARATUS AND PROCEDURE Water permeability tests were conducted using conventional techniques. The equipment is illustrated schematically in Fig. 1. Flow rates of solutions driven through small core samples at known temperatures and pressures were measured. Regulated compressed air, usually at 10 to 15 psig, was used to drive the solutions. Core samples 0.75 in. in diameter and 1 in. long were used for the most part. Cores were saturated by evacuating and then admitting a de-aerated solution. Tight core samples were first saturated with carbon dioxide gas and then evacuated and saturated with the test solution. After saturation, 140 psig pressure was imposed for several hours. Cores were tested to insure absence of gas saturation by a Boyle's

Jan 1, 1965

By B. B. Sandiford

It has been known for many years that the efficiency of a water flood can be improved by lowering the water-oil mobility ratio in the system. Such a change leads to better sweep efficiency and also to more efficient oil displacement in the swept zone. Data from our laboratory water-flood tests of both small cores and long sand packs are presented which show that water mobility can be reduced and oil recovery increased by the addition of certain polymer solutions to flood water. The reduction in mobility, in many cases, is greater than would be expected from conventional viscosity measurements. These solutions, however, do not cause significant reductions in oil mobility. The over-all effect of these mobility changes is increased waterflood oil recovery. Encouraged by results of our laboratory work, we expanded our study to include pilot field tests of floods with such solutions. One such test, made in the West Cat Canyon field, Santa Barbara County, Calif., is described in detail in this paper. Three other field tests are also discussed. INTRODUCTION Oil production from most reservoirs following primary depletion and/or water flooding is often less than 50 per cent of the original oil in place. Heavy oil reservoirs seldom yield over 15 per cent of their original oil. With new reservoirs becoming harder to find, the improvement of oil recovery efficiency is one of our very important problems. We describe here some of our attempts to increase the efficiency of oil displacement by adding a water-soluble polymer, partially hydrolyzed polyacrylamide, to flood water. This technique will be termed "polymer solution flooding". The concept of using high-viscosity water to increase the efficiency of water flooding is not new. In 1944 Detling (Shell Development Co.) obtained a patent covering the use of several additives for viscous water flooding.' His objective was to improve water-oil mobility ratios by increasing the viscosity of the flood water. Other patents2-27 have been granted covering specific water-SO~LIble polymers or specific conditions of viscous water flooding. Barnes'" described his laboratory model study of the injection of a viscous water slug into a reservoir which had been partially invaded by bottom water. He concluded that, for this type of reservoir, "the cost of viscous water should not exceed a few cents per barrel for viscous water slug injection to be economically feasible" Our studies have led us to a somewhat different conclusion in a number of cases where hydrolyzed polyacrylamide solutions have been injected into reservoir models or actual reservoirs. Possible reasons for this difference are discussed in this paper. Our studies have shown that polymer solutions may lead to an increase in oil recovery over that from an ordinary water flood by (1) improving sweep efficiency, (2) improving microscopic displacement efficiency, or (3) a combination of these mechanisms. In the work of Barnes, only the benefit of improved sweep efficiency was considered. Also, our work has shown that there are marked differences in the effectiveness of different water-soluble polymers as flood water additives. Partially hydrolyzed polyacrylamide is better than many other water-soluble polymers we have tested because, even in very low concentrations, it can lead to increased oil recovery. This is an important advantage when either a dilute polymer solution is injected continuously or a relatively concentrated slug is injected followed by water. In the latter case! portions of the slug become diluted and function in the formation as very dilute solutions. As dilution takes place the effective slug size will increase which, in turn, will reduce the cost per barrel of the effective flooding medium. The reason that partially hydrolyzed polyacrylamide solutions are more efficient at low concentrations than certain other polymer solutions of equivalent viscosity (when measured in conventional viscometers) is not fully understood. We do know that the shapes and sizes of macro-molecules dissolved or suspended in liquids influence the flow properties of their solutions or suspensions. Solutions of partially hydrolyzed polyacrylamide cause greater reductions in water mobility than would be expected from conventional viscosity measurements. LABORATORY STUDIES Laboratory water floods were run in linear and radial systems with different water-soluble polymers and under varying conditions of flow. including reservoir conditions of temperature, pressure and fluid composition. Some of these runs are considered in this section. OIL DISPLACEMENT IN LINEAR MODELS In this group of runs the sweep efficiency approached 100 per cent because the linear sand packs used were as nearly uniform as possible. Results reflect primarily the micro-scopic displacement efficiencies. The laboratory models were unconsolidat-ed sand packs the lengths of which varied from about 4 in, to 40 ft. Further information on the flow models used is listed in Table 1. Using conventional procedures, water floods were run on sand packs containing either refined or crude oil at restored state. Frequently, a repeat run, similar to the first, was made as a check. Then a third run was made on the restored-state model using either polymer solution or a combination of polymer solution (slug) and water. Details of these runs are also reported in Table 1. In Figs. 1 and 2 results of two different groups of runs are shown in which the oil used was a 62-cp refined oil and the flooding medium

Jan 1, 1965

By R. L. Dalton, L. A. Rapoport, C. W. Carpenter

Procedures have been developed to study the performance characteristics of unconfined pilot water floods using scaled laboratory models. The effects of operating conditions on pilot performance for variozrs well patterns have been investigated. The limitations of the laboratory model studies are outlined. Relationships between the oil recoveries of pilots and large-scale floods have been established. The applications of model studies to the prediction of large-scale waterflood performance and pilot flood design are discussed. INTRODUCTION Pilot water floods are undertaken to evaluate operational procedures and to gain advance information about the performance of an extensive water flood. In general, however, the performance of a pilot flood is not directly indicative of what can be expected from a large-scale pattern water flood. This difference stems from the fact that the over-all flow configurations are not the same for the two systems. A water flood normally consists of an array of identical well patterns. In such an array, the perimeters of the well patterns are axes of symmetry and act as impermeable boundaries. Thus, an extensive, pattern water flood can be visualized as a repetition of "confined" floods. In contrast, a pilot water flood involving only one or a few well patterns is "unconfined". Here the well patterns are not balanced by other flood units; hence, the perimeter of the pilot area does not act as an effective boundary. Accordingly, only a portion of the fluids injected (and also of those originally contained within the pilot area) is captured by the producers of the pilot. The rest escapes into the surrounding reservoir. This usually causes the amounts of oil and water produced in the pilot to be different from those recovered by a producer in an extensive water flood. Moreover, the recovery observed in the pilot can be expected to be greatly influenced by operating conditions, i.e., by the bottom-hole pressures at the injection and producing wells. The above considerations show that full utilization of pilot test results and reliable prediction of large- scale water floods require a knowledge of the relationship between the recovery performance of confined and unconfined well patterns. To date, only limited basic information is available on the behavior of pilot water floods. Craig, et al,' have investigated the effects of initial gas saturation, oil-to-water mobility ratio, and well damage on the oil recovered from a single, "normal" five-spot pilot pattern (Fig. 1). Rosenbaum and Matthews2 studied the effect of initial gas saturation and mobility ratio on the ratio of production rate to injection rate for various groupings of five-spot patterns. Paulse113 investigated the effect of mobility ratio and distance from the pilot to the edge of the reservoir on the areal sweep in a single "inverted" five-spot (Fig. 1). However, none of these studies considered the effect of operating conditions on pilot flood behavior. Moreover, hardly any of the work to date dealt explicitly with the oil-recovery performance of the pilot floods. The purpose of this paper is to introduce the concept of interrelation between pilot flood behavior and operating conditions, to illustrate this interrelation by means of experimental data and to compare the recovery performances of dif-

By T. K. Perkins, R. N. Hoffman, O. C. Johnston

A simplified theory of viscous fingering in miscible systems has been developed. It predicts the correct functional relationship between pertinent variables and permits the calculation of the order of magnitude of fingering behavior for simple systems, such as linear or radial. The theory is based on the following four observations: (1) cross-flow takes place only near the ends of fingers; (2) the relative finger width is about 0.5; (3) fingers can be suppressed by transverse dispersion, the suppression being quantitatively described by a critical ratio of dispersion times; and (4) the widths of extending fingers are close to the minimum size finger that can grow at any stage of displacement. Fingering is studied in two-dimensional linear and diverging radial systems, both theoretically and experimentally. For linear systems, the length of the fingered region is proportional to mean displacement; the finger width is proportional to square root of mean displacement; and there is a small initial region void of fingers because of suppression by longitudinal dispersion. For the radial system, two limiting conditions are recognized. If the mean displacement is small compared with the wellbore radius, the length of the fingered region is proportional to the mean displacement. The width is proportional to the square root of mean displacement. If the mean displacement is large compared with wellbore radius, length of the fingered region is proportional to mean displacement, but the number of fingers approaches a constant value. Also, in the radial case there is a small initial region void of fingers because of longitudinal dispersion. INTRODUCTION The behavior of viscous fingers in miscible systems has been of interest to the oil industry for many years. Previous studies have clearly shown the existence and nature of fingers in small models.1,2,3,4 Engineering techniques for extrapolating to reservoir situations have been proposed. 5,6 Still, because of the lack of a real understanding of the mechanics of fingering, there remains uncertainty and disagreement as to the best way to scale or calculate fingering behavior in the reservoir. In this paper we discuss a study of the mechanics of fingering in miscible systems (i.e., why do fingers behave as they do?). A simplified theory is developed which we believe will (1) predict the correct functional relationship between pertinent variables, (2) permit us to calculate the order of magnitude of fingering behavior for simple systems such as linear and radial, and (3) give further insight into the problem of scaling or otherwise extending the results to more complicated reservoir conditions. The paper includes the following sections: (1) brief summary of the mixing behavior of miscible fluids in linear and radial systems; (2) four fundamental observations of fluid flow (under fingering conditions); (3) based on these four observations, a simplified theory of fingering in linear and radial systems is developed; and (4) the theoretical equations are compared with experimental fingering data measured in laboratory models. A REVIEW OF DIFFUSION AND DISPERSION As will be shown later, the behavior of viscous fingers in miscible systems is controlled in large degree by the mixing between the two fluids. A quantitative description of fingers will first require a quantitative description of mixing. Fortunately, much work has been done to clarify this subject; a fairly comprehensive review has been presented in a previous paper. 7 For present purposes let us (1) briefly summarize the quantitative description of diffusion and dispersion (both longitudinal and transverse) within a differential element; (2) present a simplification of the "width of mixed zone"; and (3) describe the effect of geometry on width of mixed zone. If one fluid is displaced from a porous medium by another miscible fluid (in the absence of fingers)

Jan 1, 1966

By F. W. Jessen, H. Kazemi

A washing technique has been developed to form a spherical cavity in massive salt. The technique is, basically, a process of controlling the fluid motion in the cavity, the concentration distribution of brine, the rate of dissolution at the walls of the cavity, and the particular settings of the wash pipe, or pipes (inlet and outlet pipe) employed. For a particular arrangement of tubing and casing, fresh water enters from the casing annulus at a position above the bottom of the tubing, and the dissolved material leaves the cavity through the tubing near the bottom of the projected hole. The initial configuration of the cavity is cylindrical, representing the drilled hole situated at the axis of the projected cavity. As washing progresses the original hole is enlarged radially (any direction perpendicular to the axis of the hole). An inert fluid (generally a hydrocarbon, natural gas, or air) which serves as a blanket is injected through the annular space between the original hole and the casing. The controlled downward motion of this blanket maintains the upper edge of the brine region that is expanding on the surface of the projected sphere. A large spherical glass container and a radial model were used to investigate the flow patterns and the concentration distribution of brine occurring during the progress of solution of the salt. Aluminum flakes and dye were injected and flow patterns were observed and recorded on motion picture film. Based on the results of these model studies and other work carried out in this laboratory, the controlled washing technique was formulated. INTRODUCTION Underground cavities leached in salt beds and salt domes have been used extensively for storing hydrocarbons. Little attention has been focused on the shape of the cavities as long as large volumes were obtained and the cavities were structurally stable. The general process involves circulating fresh water and removing brine. Two distinctly different circulating systems can be employed, each of which yields characteristically different shapes. By far the most generally accepted method of washing storage caverns in salt domes is the direct circulation system. Water is introduced through the innermost pipe of a concentric pipe system, and brine is withdrawn through the annulus. In the reverse-circulation method, water enters through the annulus and is discharged through the smallesc tubing string. The size of the operation, and hence the size of tubing and casing used, depends on such factors as cavity volume, length of hole, availability of water, disposal of brine and power requirements. Underground caverns are from 100 to 1,000 ft or more in length, and when completed vary from 40 to 100 ft in diameter, providing storage of 100,000 bbl or more. The volumes of water required per volume of salt removed vary from 6 to 7 for the reverse circulation method to 10 to 11 for the direct circulation method. Up to 80,000 BWPD are often circulated. In the direct circulation method, because water enters near the bottom of the projected cavity, solution of salt is greatest at this point, particularly in view of the fact that when the washing is begun turbulent flow undoubtedly exists in the immediate vicinity of the discharge of the inlet water. As the diameter of the hole increases, flow rates also are increased, so perhaps the condition of fluid turbulence is maintained for some time. During the later stages of leaching, when capacity of equipment may be the limiting factor in maintaining flow rates, there may still exist some turbulence adjacent to the point of influx; however, for the most part laminar flow conditions prevail. There appears, therefore, rather good evidence to indicate that the state of fluid motion is responsible for the typically shaped cavity which results from such washing — namely a cylindrical form, generally somewhat larger at the bottom than at the top. Depending on the number of so-called stages which might be employed, the walls or surfaces of the salt diminish rather uniformly from bottom to

Jan 1, 1965

By C. R. Johnson, R. A. Greenkorn, R. E. Haring

Three confined five-spot miscible displacements at unity, favorable, and unfavorable mobility ratios were conducted in a shallow, water-saturated sandstone of Pennsylvan-ian age near Chandler, Okla. These studies, plus associated laboratory experiments, were designed to measure miscible displacement performance in a controlled natural system, using known scaling criteria to develop an approach to modeling the heterogeneous field system. We have concluded from these studies that: (I) displacement efficiency in the field is a pronounced function of mobility ratio, indicating that miscible fingering observed in simple laboratory models occurs in the field: (2) field displacements can be quantitatively predicted by scaled laboratory models if the degree and location of field permeability variations are preserved in the models: and (3) arbitrary simplifications of heterogeneity will not necessarily predict observed displacement efficiency, and the simpler the model, the more optimistic the prediction. INTRODUCTION Many processes to achieve miscible displacement of reservoir oil by injected fluids have been conceived and field tested by the oil industry. Among the better known are high-pressure gas, enriched gas, and LPG banks. The simplest form of miscible displacement—one fluid miscibly displacing another fluid of different viscosity but the same density—has been studied extensively in homogeneous laboratory models. Observations of unstable fingering have been made which explain the significant decrease in displacement efficiency as the mobility ratio (ratio of displacing fluid mobility to displaced fluid mobility) increases. General industry experience with field tests of miscible displacement projects, mostly LPG banks, has been premature solvent breakthrough and lower than predicted production rate increases. These results have been attributed to either unstable fingering, unusual or unexpected permeability stratification, or both. Miscible displacement data in a controlled natural system have never been reported, however. Also, it has not been shown that properly designed and constructed laboratory models quantitatively predict field-scale behavior. The purpose of the combined field and laboratory ex- periments reported in this paper is twofold. The first was to measure miscible displacement performance at different mobility ratios in natural rock approaching field size under precise, controlled conditions. The second purpose was to utilize known scaling criteria plus several approaches to heterogeneity to model the field. Comparison of model and actual field results should then determine whether or not the laboratory phenomena (manifested by miscible displacement efficiency) are exhibited in large, natural rock systems. We carried out our program by first locating a shallow, water-saturated reservoir whose rock properties were representative of oil-bearing reservoirs. Detailed reservoir description by core analysis and interference testing showed the field site to be heterogeneous. A sequence of controlled, aqueous-phase miscible displacements was conducted at unity, favorable and unfavorable mobility ratios. A central, confined pattern was used to obtain the displacement data. A laboratory program using sand-packed models was conducted to determine the modeling criteria necessary to simulate field behavior of miscible displacement in a heterogeneous system. SCALING THEORY The detailed derivations and descriptions of the scaling laws that apply to laboratory models of reservoirs are adequately described elsewhere, so the following discussion will be restricted to facets of importance in this study. For a displacement in which one liquid miscibly displaces another, the following dimensionless groups are required to have the same numerical value in the model as in the field: The model also must be geometrically similar to the field, be spatially oriented the same as the field (same dip angle), and have the same initial and boundary conditions as the field (same initial fluid saturations and same injection-production well arrangement). When these conditions are satisfied, the theory predicts that at dny dimensionless time (pore volumes of produced fluids) the dimen-sionless flow potential and dimensionless fluid concentrations will be identical at all dimensionless spatial locations within the model and the field. If this prediction is correct, then the local dimensionless velocities must be identical, thus the instantaneous fraction of displaced fluid produced and the cumulative recovery expressed as fraction of original fluids in place must be identical at all dimensionless times. The theory outlined above has been obtained by either dimensional analysis or inspectional analysis of differential

Jan 1, 1966

By J. F. Wilson

An experimental investigation has been made of gas-driven slug displacements in a system of high gas saturation to evaluate the process for use in a California reservoir. Fluid compositions, temperature, pressure and core permeability duplicated reservoir conditions as closely as possible. The temperature was above the critical of the slug materials. The slug composition required for 100 per cent oil recovery in the linear flow systems was found to agree with that determined from equilibrium phase studies. A qualitative theory for the slug displacement of two-phase systems is proposed and found to be in general agreement with the observed flow behavior. The theory predicts that the length of the reservoir fluid-slug transition Zone will increase with increased initial gas saturation and decrease with increased oil swelling accompanying the composition change from reservoir oil to the critical composition. The precipitation of small amounts of an asphaltic phase during the formation of the transition zone had no adverse effects on the flow behavior of otherwise miscible systems. INTRODUCTION Despite the large number of recent publications relating to oil recovery by miscible displacement, little information is available concerning the behavior of crude oil-natural gas systems at conditions likely to be encountered in field application. The laboratory data reported here were obtained during an investigation of the possibility of using a gas-driven LPG slug in a partially depleted California reservoir. Seven gas-driven slug displacements employing slugs of different compositions and sizes were made in 26- to 27-ft long sandstone cores. The equilibrium phase behavior of the systems was also determined. The features of primary interest in this study are the following. 1. Data were obtained on the displacement of reservoir fluids at reservoir conditions. 2. A high gas saturation existed in the long cores at the start of the displacements. 3. The temperature of the studies was above the critical of the slug materials. 4. A comparison between the flow behavior and equilibrium phase behavior, as illustrated by a pseudo-ternary phase diagram, is afforded for a system of considerable complexity. THEORY Phase Behavior The use of the pseudo-ternary phase diagram to illustrate the phase behavior in miscible displacements has been discussed by Clark, Schultz and Shearin. Sie-vert, Dew and Conley and Arnold, Stone and Luffel' have presented both phase equilibrium and flow data for true ternary systems. For the true ternary, the application of the theory is straight-forward; for a reservoir system represented as a ternary, the theory undoubtedly gives a qualitatively correct picture. One purpose of this investigation was to see if the pseudo-ternary phase diagram could be used to predict with reasonable accuracy the conditions necessary for miscible displacements with actual reservoir systems. MiscIBLE Displacement of Two-Phase System A large number of papers "I3 have treated the miscible displacement of a single-phase fluid from both an experimental and a theoretical viewpoint. However, the normal field application of the process will be in a partially depleted reservoir. Often the state of the system will be such that the mobility of the gas phase is greater

By G. T. Pruitt, T. L. Irby, P. B. Crawford, H. Ferrell

In water flooding peripheral, center-to-edge, line-drive or water-encroachment patterns the question has arisen, "when should a producing well be converted to a water-injection well?". It is realized that, if the producing well is shut in when water breakthrough occurs, some oil could be left in an unswept area. Yet, operation to high water-oil ratios may require processing much additional water. If a well is shut in, "where should the allowable be transferred?". A review of the literature indicates that very few data are available to the engineer to assist him in answering these questions. It is the purpose of this paper to present information which may be of value in selecting the flooding technique and the injection and producing wells and in scheduling injection and production rates. Certain phases of these problems can be studied in the laboratory. The Potentiometric model is well suited for studying many waterflooding problems. This model assumes that steady-state flow exists, that the mobility ratio is equal to one and that gravitational effects are neglected.'-' The three-phase sand model also has been recognized as an effective means for studying certain reservoir problems. Both models were used in making the study reported here. This paper discusses seven cases utilizing various waterflooding techniques. Data for the first five cases were obtained on the potentiometric model. Data for Cases VI and VII were obtained by constructing and flooding an unconsolidated sand pack (a three-fluid-phase system) and by determining the swept area and cumulative oil recovery for various types of operations. APPARATUS AND PROCEDURE The potentiometric model reservoir was 30 in. long X 10 in. wide by about 1/2 in. deep. Four wells were equally spaced in a line array along the lateral center line. (See Fig. 1.) Well No. 1 initially was an injection well for each of the patterns reported. The cases studied with the potentiometric model and sand model are described herein. The laboratory sand model reservoir was a 30- 10-in. rectangular Lucite vessel packed with sand to a depth of about 0.65 in. Four wells were equally spaced in a line array along the lateral center line. The reservoir material was an unconsolidated sand into which the reservoir fluids were mixed to yield fluid saturations of 55 per cent oil, 25 per cent water and 20 per cent gas. The water was added to the sand and the material mixed, after which the oil was added. Kerosene and water, the fluids used, gave a viscosity ratio of 1.64:1. In packing the reservoir, care was taken to assure constant fluid saturations and constant porosities. A new mixture was prepared for each run. During water flooding, a colored water was injected by a constant rate pump. The injection rate was 1.1 B/D/ft of sand thickness. For this series of tests, each of the producing wells was open to the atmosphere and, as in the potentiometric model study, Well No. 1 initially was an injection well for each of the patterns. KESULTS OF POTENTIOMETKIC MODEL STUDlES Two cases (Cases I and 11) involving equal pressure injections were studied, with the following results. Case Water injection was initiated at Well No. 1; the producing wells (Nos. 2, 3 and 4) were operated at constant and equal bottom-hole pressures. When water broke through at Well No. 2 and achieved the selected water-oil ratio, Well 2 was converted to a water injection well, resulting in two injection and two producing wells. The injection pressures at Wells 1 and 2 were maintained equal and constant; pressures at the producing wells were held equal, but they were lower than those at the injecting wells. A similar process was repeated when water broke through at Well No. 3. The flood was terminated when

By B. H. Caudle, R. E. Watson, I. H. Silberberg

The production history of an inverted nine-spot injection pattern was studied with a porous plate model similar to Habermann's. Colored fluids were used so the boundary movements could be photographed. By varying the production rates of the side and corner producing wells, the optimum pattern sweep efficiency can be obtained. Results of the laboratory tests indicate that interpolation of the pattern sweep and conductivity ratio data will permit application of the data to predicting fluid injection operations with mobility ratios between 0.1 and 10 and production rate ratios between I and the optimum values. INTRODUCTION The inverted nine-spot injection pattern in Fig. 1 has some unique advantages over other regular injection patterns (five spot, line drive, etc.) used in secondary recovery operations. It has three producing wells for each injection well, while the five-spot and line-drive patterns require an equal number of injectors and producers. Futher-more, of the eight producing wells surrounding each injection well, four are of the "side well" type which lie nearer the injection well than the remaining four "corner wells". This allows the production rates of these two types of wells to be adjusted independently of each other to obtain better pattern sweep efficiency. This latter advantage may well be the reason why the nine-spot patterns have not been studied as extensively as the patterns containing only one type of injection or producing well. The pattern sweep efficiency of the inverted nine spot is affected almost as much by operating with different producing rates of the comer and side wells as by differing mobility ratios. Possibly because of this added operating parameter, published studies of nine-spot systems have been limited in their coverage of the pertinent variables.'-' The model studies reported here covered the mobility-ratio range normally found in waterflood operations, as well as the variation in well production rates necessary to achieve optimum pattern sweep efficiency. Conductivity ratios — a measure of the change in total flow at constant pressure difference — were obtained during each model run. A porous-plate model with miscible fluids representing the in-place and injected fluids was used in this study. This model technique is the same as that used in both X-ray shadowgraph" and visual techniques4 described in the literature. Using miscible fluids eliminates residual oil in the model, and the permeability to the injected fluid is the same as to the in-place fluid. Thus the mobility ratio in the model is given by the ratio of the viscosity of the displaced fluid to the viscosity of the injected fluid. Because of the non-linear nature of the viscosity vs composition curve for two miscible liquids, the mobility change across the boundary between the two fluids is reasonably sharp. This technique thus represents the displacement of a region of one fluid mobility by a region of possibly different mobility in an isotropic system under steady-state flow conditions. Table 1 shows the dimensionless similarity groups used to scale the model conditions and results to

Jan 1, 1965

By A. L. Jr. Pozzi, A. O. Garder, D. W. Peaceman

A new numerical method is proposed for the solution of multidimensional miscible displacement problems. Besides the usual stationary grid associated with numerical procedures, the method uses moving points for which positions and concentrations are computed each time step. Based on the method of characteristics for treating combined transport and dispersion, the method applies equally well to any number of dimensions and, in contrast with other numerical procedures, properly takes into account the physical dispersion, no matter how small. Extensive tests show that accurate solutions of one-dimensional problems can be obtained over a wide range of values of the dispersion coefficient, including zero. They also show that the moving points do not need to be uniformly spaced, and that increasing the number of moving points beyond 2/grid interval does not significantly improve the accuracy of the solution. Two-dimmsional calculations using the new method were carried out to simulate laboratory model displacements in which gravity caused overriding of solvent. Good agreement between observed and computed profiles and measured and computed concentrations of produced fluid was obtained for a viscosity ratio of 1.89; at higher viscosity ratios, fair agreement was obtained Computed results were found to be significantly affected by vertical dispersion, and negligibly affected by dispersion in the direction of flow, in agreement with conclusions based on experimental observations. INTRODUCTION Although miscible displacement processes are potentially of great economic significance, limited progress has been made in developing mathematical methods for predicting the outcome of a solvent flood. Because of the complexity of the partial differential equations which describe multidimensional miscible displacement, numerical methods are a natural approach to their solution. The task of finding suitable numerical approximations to one of the partial differential equations, which involves both dispersion and convective transport, has been particularly difficult. Dispersion in the absence of transport is described by the heat flow equation, a second-order equation of parabolic type. This equation has been very successfully treated by numerical methods. Transport in the absence of dispersion is described by a first-order equation of hyperbolic type. This equation has been treated numerically with some success in one dimension, both by Lagrangian and Eulerian techniques, 1.2 but extension to two or more dimensions has not been satisfactorily accomplished. Usually one of two things happens: either the numerical solution develops oscillations or it becomes smeared by artificial dispersion resulting from the numerical process. When transport and dispersion are considered simultaneously, this numerical dispersion may dominate the low physical dispersivicy that generally characterizes miscible displacement. Stone and Brian2 have developed a new difference equation for treating combined transport and dispersion in one dimension. By suitable choice of certain parameters, they succeeded in markedly reducing both oscillation and numerical dispersion. However, a satisfactory extension to higher dimensions has not yet been found. Peaceman and Rachford3 presented a centered-in-time and centered-in-distance difference equation combined with a 1,transfer of overshoot" procedure which was demonstrated to work well in one dimension. The technique was applied to the two-dimensional miscible displacement problem, but subsequent calculations with zero dispersivity showed little change in results compared with the calculations they presented. This indicates that their method involved a numerical dispersion of the same order as the physical dispersion. The present paper describes an alternate method for solving the combined transport-dispersion

Jan 1, 1965

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 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