Instability in Water Flooding Oil from Water -Wet Porous Media Containing Connate Water

1964 ◽  
Vol 4 (02) ◽  
pp. 133-148 ◽  
Author(s):  
H.H. Rachford

Rachford Jr., H.H., Member AIME, Humble Oil and Refining Co., Houston, Tex. Abstract This work presents a first-order analysis of the instability underlying viscous fingering in adverse viscosity-ratio water floods. It extends previous analyses of frontal instabilities, which were carried out with equations for parallel plate models, by including effects of the saturation transition zone observed behind the front in water floods in water wet systems. This zone tends to insulate incipient fingers from the high-mobility water; thus conditions for the onset of fingering differ from those in the parallel plate theory. Finite-difference solutions of the two-dimensional equations of displacement in porous media exhibited the predicted stability characteristics in six hypothetical field- and laboratory-scale floods in rectangular reservoirs. In contrast to results with parallel plate systems, this paper concludes that for water-wet reservoirs, laboratory models scaled by the usual criteria are also correctly scaled for frontal instability. Further, fingering in the systems studied can occur in any saturation range behind the front, and may occur at an intermediate saturation even though stability obtains both at the saturation corresponding to the Buckley-Leverett front and near residual oil saturation. Other points of contrast are that the likelihood of occurrence of fingering may not increase as flow rate or viscosity difference increases, but may be sensitive to changes in the relative permeability and capillary pressure functions. Introduction The recovery of oil by water flooding frequently involves displacing the oil by water of a lower viscosity. Displacement of a fluid by a less viscous one may lead to gross channeling or fingering like that observed in solvent floods, in which it severely lowers recovery efficiency. In addition to adverse effects on recovery, it has been suggested that the unstable movement causing fingering may interfere with interpretation of scaled model studies of proposed water floods, since the instability in the model might not be faithfully scaled to that in the reservoir prototype. In view of the serious implications of this possible breakdown of widely used scaled model techniques, it is the purpose of this paper to examine the question further. Instabilities in the solutions of systems of differential equations imply a loss of smooth dependence on initial and boundary conditions. Thus, the possibility exists that in using models whose scaling is based on the differential system there may arise size- and rate-dependent factors which are not properly scaled. This possibility was examined in detail by Chouke et al., who analyzed the instability of frontal advance in a related problem, water-oil displacement in parallel plate models, in which a moving interface separates two regions of constant, unequal mobilities. First-order perturbation theory predicts the existence of a critical wave length for the growth of perturbations: and a wave length of maximum instability of . The interpretation is that wave lengths in a perturbation which are longer than will grow. Thus, if the width of a two-dimensional channel is greater than, fingers will grow, and the spacing of the fingers which grow at the maximum rate will be approximately. It is important that the higher the velocity and/or the difference in flow resistance, the lower is, and thus the greater the number of fingers that can grow in a given model. In applying these conclusions to porous media, a pseudo-interfacial tension, was assumed for the invasion front. Since this would not necessarily be equal to the liquid-liquid interfacial tension, an unknown constant was substituted for in the foregoing expression for. JPT P. 133ˆ

1964 ◽  
Vol 4 (03) ◽  
pp. 231-239 ◽  
Author(s):  
A.S. Michaels ◽  
Arnold Stancell ◽  
M.C. Porter

MICHAELS, A.S., MASSACHUSETTS INSTITUTE OF TECHNOLOGY, CAMBRIDGE, MASS. MEMBER AIME STANCELL, ARNOLD, MASSACHUSETTS INSTITUTE OF TECHNOLOGY, CAMBRIDGE, MASS. PORTER, M.C., MASSACHUSETTS INSTITUTE OF TECHNOLOGY, CAMBRIDGE, MASS. Abstract 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. 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 through put), 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 aqueous bexylamine, is a result of transient changes in the oil wettability of the pore walls. If the am in e 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 am in e, a favorable "wettability 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 estimates 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, 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. 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, 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 Timmins 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. SPEJ P. 231^


Author(s):  
Sergey I. Zhavoronok

he extended plate theory of I.N. Vekua – A.A. Amosov type is constructed on the background of the dimensional reduction approach and the Lagrangian variational formalism of analytical dynamics. The proposed theory allows one to obtain the hierarchy of refined plate models of different orders and to satisfy the boundary conditions on plates’ faces exactly by introducing the corresponding constraint equations into the Lagrangian model of two-dimensional continuum. The normal wave dispersion in an elastic layer is considered, the convergence of the two-dimensional solutions to the exact one is studied for the locking phase frequencies, the dimensionless stress distributions across the thickness of a layer are shown.


2005 ◽  
Vol 03 (02) ◽  
pp. 195-221 ◽  
Author(s):  
KARIM TRABELSI

A new two-dimensional nonlinear membrane plate theory is derived via a formal asymptotic procedure for a family of hyperelastic nonlinear materials proposed by Ciarlet and Geymonat [11], whose stored energy function is polyconvex and becomes infinite, when the determinant of the deformation gradient tends to zero, and can be adjusted to arbitrary Lamé constants.


2017 ◽  
Vol 19 (3) ◽  
pp. 120-127
Author(s):  
S.Y. Kudryashov ◽  
L.A. Onuchak ◽  
D.A. Panarin

Generalization of chromatographic plate theory in case of two-dimensional non-sorbing tracer flow in a homogeneous porous media of uniform thickness with equal debit injection and producing wells are considered. It is shown that in this approach approximate calculation of tracer concentration at different times in producing well, and at a fixed time at various porous media points can be carried out.


2008 ◽  
Vol 8 (2) ◽  
pp. 101-115 ◽  
Author(s):  
B.N. CHETVERUSHKIN ◽  
N.G. CHURBANOVA ◽  
M.A. TRAPEZNIKOVA ◽  
A.A. SUKHINOV ◽  
A.A. MALINOVSKIJ

AbstractThis paper considers two-dimensional hierarchical locally-refined rect-angular meshes with dynamic adaptation to the solution. A parallel algorithm with dynamic load balancing has been developed. Adaptive meshes have been used for the problem of passive contaminant transport in an oil-bearing stratum at water flooding. A model neglecting the capillary and gravity forces (the Buckley — Leverett model) has been used.


1962 ◽  
Vol 2 (02) ◽  
pp. 165-176 ◽  
Author(s):  
H.D. Outmans

Abstract Present first-order theory for frontal stability and viscous fingering of immiscible liquids is improved by including the nonlinear terms in the equations describing conditions at the interface of the liquids. This leads to the revision of several conclusions which were based on first-order theory. They concern the growth rate and changing shape of a sinusoidal disturbance, particularly in relation to the wave number of this disturbance. Nonlinear aspects of effective inter facial tension are discussed and it is shown that this tension is not simply a positive proportionality constant in a linear relation between pressure difference and curvature at the interface. Scaling requirements are determined from the dimensionless groups which govern fingering. Gravity and interfacial tension invalidate a previously formulated conclusion that the shape of a finger after a given displacement is independent of the displacement velocity. Also, similarity of fingering (and, hence, of sweep efficiency in case of an unstable front) requires geometrical similarity of the initial disturbance in the model and the reservoir with a scale factor which is the same as the one for scaling down the dimensions of the reservoir. Introduction The critical velocity at which the transition zone in the vertical displacement of immiscible liquids becomes unstable was calculated by Hill. More recently, the same was done for nonvertical displacements and it was also shown that the stability is affected by an effective interfacial tension between the liquids. Although in the oil reservoir a stable front is always desirable from a recovery point of view, the necessary velocity may have to be so low that the corresponding production rate is not economical. In that case, a knowledge of the fingering unstable front is necessary for predicting recovery at breakthrough. Studies in this direction have not gone beyond the very moment at which fingering first occurs. Conclusions about stability and fingering in these references are all based on linear theory. In this theory the nonlinear terms in the equations describing conditions at the front are neglected. The usual justification for this lack of mathematical rigor is that the nonlinear terms can be made small relative to the linear terms and, supposedly, small causes produce only small effects. It has been recognized, however, that this is not necessarily true in the nonlinear problems of hydrodynamics. The hydrodynamic stability and fingering of the front between two liquids, accelerated in a direction normal to the front, for instance, is strongly influenced by the nonlinearity of the boundary conditions. Since these boundary conditions are not unlike those arising in the problem of stability and fingering of the front during a slow immiscible displacement in porous media, it was thought that there, too, the nonlinear terms should be taken into consideration. A summary of the results obtained in linear theory and a comparison with experimental data precede the nonlinear theory developed in this paper. This summary serves to introduce some quantities which will be used in the further development. It should also be pointed out that the method of higher-order approximation by which the nonlinear stability problem is solved has application in other reservoir studies. SUMMARY OF LINEAR THEORY A plane interface remains stable if its velocity is smaller than a critical velocity. ............................(1) This equation is derived for uniform flow in a thin layer inclined to the horizontal plane. The displacement takes place in a direction normal to the intersection of this layer and the horizontal plane. The undisturbed (line) interface is normal to the velocity. Fluid 1, the upper fluid (gas), is displacing Fluid 2 (oil). SPEJ P. 165^


1999 ◽  
Vol 2 (3) ◽  
pp. 251-262
Author(s):  
P. Gestoso ◽  
A. J. Muller ◽  
A. E. Saez

2021 ◽  
Vol 11 (8) ◽  
pp. 3421
Author(s):  
Cheng-Yu Ku ◽  
Li-Dan Hong ◽  
Chih-Yu Liu ◽  
Jing-En Xiao ◽  
Wei-Po Huang

In this study, we developed a novel boundary-type meshless approach for dealing with two-dimensional transient flows in heterogeneous layered porous media. The novelty of the proposed method is that we derived the Trefftz space–time basis function for the two-dimensional diffusion equation in layered porous media in the space–time domain. The continuity conditions at the interface of the subdomains were satisfied in terms of the domain decomposition method. Numerical solutions were approximated based on the superposition principle utilizing the space–time basis functions of the governing equation. Using the space–time collocation scheme, the numerical solutions of the problem were solved with boundary and initial data assigned on the space–time boundaries, which combined spatial and temporal discretizations in the space–time manifold. Accordingly, the transient flows through the heterogeneous layered porous media in the space–time domain could be solved without using a time-marching scheme. Numerical examples and a convergence analysis were carried out to validate the accuracy and the stability of the method. The results illustrate that an excellent agreement with the analytical solution was obtained. Additionally, the proposed method was relatively simple because we only needed to deal with the boundary data, even for the problems in the heterogeneous layered porous media. Finally, when compared with the conventional time-marching scheme, highly accurate solutions were obtained and the error accumulation from the time-marching scheme was avoided.


2021 ◽  
Vol 2021 (4) ◽  
Author(s):  
Clifford V. Johnson ◽  
Felipe Rosso

Abstract Recent work has shown that certain deformations of the scalar potential in Jackiw-Teitelboim gravity can be written as double-scaled matrix models. However, some of the deformations exhibit an apparent breakdown of unitarity in the form of a negative spectral density at disc order. We show here that the source of the problem is the presence of a multi-valued solution of the leading order matrix model string equation. While for a class of deformations we fix the problem by identifying a first order phase transition, for others we show that the theory is both perturbatively and non-perturbatively inconsistent. Aspects of the phase structure of the deformations are mapped out, using methods known to supply a non-perturbative definition of undeformed JT gravity. Some features are in qualitative agreement with a semi-classical analysis of the phase structure of two-dimensional black holes in these deformed theories.


Sign in / Sign up

Export Citation Format

Share Document