Fluid Distributions During Immiscible Displacements in Porous Media

1966 ◽  
Vol 6 (03) ◽  
pp. 261-266 ◽  
Author(s):  
L.L. Handy ◽  
P. Datta

Abstract For a wetting phase displacing a nonwetting phase from a porous medium the distribution of the residual fluid may depend on displacement conditions. Although this subject has been debated in the literature, only a few, experiments have been cited to support the various conclusions. Experimental results presented in this paper show that fluid distributions are dependent on imbibition procedures. Results agree qualitatively with predictions from the pore doublet model. if the rate of water imbibition is restricted, the nonwetting phase is trapped preferentially in the larger pores as expected. But if the rate of water imbibition is unrestricted, trapping occurs somewhat more in the smaller pores. These conclusions were deduced primarily from relative permeability measurements. Introduction Relative permeabilities are known to depend on the saturation history of the porous medium. For either continuously increasing or continuously decreasing wetting phase saturations, however, they have normally been assumed to be single-valued functions of saturation. Several investigators have compared relative permeabilities measured by different method and have reported acceptable agreement. Studies of simple models of the displacement process suggest that the distribution of residual fluids can be influenced by the displacement method. If this is so, relative permeabilities also should depend on the method of displacement. These predictions have not been supported by experimental data. The object of this paper is to investigate experimentally some of the predictions of these model studies. The particular question of importance is whether steady-state and unsteady-state methods should be expected to give the same values of relative permeability. Much of the reported data showing agreement between steady-state and unsteady-state methods have been obtained on unconsolidated sand. Johnson, Bossler and Neumann, however, compared steady-state and unsteady-state results for Weiler sandstone and found no significant difference. Levine made a detailed study of pressure and saturation distributions for a laboratory waterflood in a consolidated system but made no attempt to calculate performance from steady-state relative permeability data. Bail and Marsden in a somewhat similar set of experiments did attempt a comparison of observations with predictions but the results were inconclusive. Excluding gravity, two forces affect the distribution of wetting and nonwetting phases during an unsteady-state, immiscible displacement: viscous and capillary forces. If relative permeabilities are indeed the same when measured by steady- or unsteady-state methods, the fluid distributions at the same fluid saturations must be similar for both methods. Furthermore, the implication is that the relation between capillary and viscous forces is the same for both methods insofar as the effect on microscopic fluid distributions is concerned. Unsteady-state displacements are normally run at high rates to maximize the ratio of the viscous to capillary forces. The objective is to reduce the effect of capillary forces on macroscopic fluid distributions. For floods in which the phase which wets the porous medium displaces the nonwetting phase, capillary forces compete with viscous forces in determining fluid distributions. The residual nonwetting phase is discontinuous. Competitive aspects of viscous and capillary forces and the resulting effect on water-oil distribution in porous media have been illustrated in an elementary way using a pore doublet model. This model and predictions from it have been discussed at considerable length by Rose and Witherspoon, Rose and Cleary and by Moore and Slobod. Rose and Witherspoon show that so long as water invades the model at a rate equal to or greater than the free imbibition rate, the water will move through the larger capillary of a doublet first and trap oil in the smaller capillary. SPEJ P. 261ˆ

1979 ◽  
Vol 19 (01) ◽  
pp. 15-28 ◽  
Author(s):  
P.M. Sigmund ◽  
F.G. McCaffery

Abstract With typical heterogeneous carbonate coresamples, large uncertainties of unknown magnitudecan occur in the relative permeabilities derived using different methods. This situation can beimproved by analyzing the recovery and pressureresponse to two-phase laboratory displacement tests by a nonlinear least-squares procedure. Thesuggested technique fits the finite-differencesolution of the Buckley-Leverett two-phase flowequations(which include capillary pressure) to theobserved recovery and pressure data. The procedureis used to determine relative-permeability curves characterized by two parameters and their standarderrors for heterogeneous cores from two Albertacarbonate reservoirs. Introduction Several recent investigations have recognizedpossible problems when obtaining reliable two-phasedisplacement data from heterogeneous carbonate core samples. Huppler stated that waterfloodresults on cores with significant heterogeneitiescan be sensitive to flooding rate, core length, andwettability, and that these effects should beconsidered before applying the laboratory results atfield flooding rates. Brandner and Slotboomsuggested that realistic displacement results maynot be obtainable when vertically flooding aheterogeneous core with a nonwetting phase becauseof the fluid's inability to maintain a properdistribution when the sample length is less than the height of capillary rise. Ehrlich noted thatstandard relative-permeability measurement methodsusing core plugs cannot be applied when the media are heterogeneous. Archer and Wong reported that application of theconventional Johnson- Bossler - Neumann (JBN)methods for determining relative permeabilities froma waterflood test could give erroneous results forheterogeneous carbonate as well as for relativelyhomogeneous porous media having a mixed wettability (see Refs. 1, 6, and 7). The observedstepwise or humped shape of water relativepermeability curves mainly were attributed to theeffect of water breakthrough ahead of the main floodfront entering into the JBN calculation. Archer andWong suggested that such abnormally shapedrelative-permeability curves do not represent theproperties of the bulk of the core sample, and proposed the use of a reservoir simulator forinterpreting laboratory waterflood data. The work referred to above provides the majorbackground for this study involving the developmentof an improved unsteady-state test method tocharacterize the relative-permeability properties ofheterogeneous carbonate core samples. The methodcan be applied to all porous media, regardless ofthe size and distribution of the heterogeneities.However, the presence of large-scaleheterogeneities, especially in the form of vugs, fractures, and stratification, could cause the derivedrelative-permeability relations to be affected by viscosityratio and displacement rate. Remember also that extrapolation of any core test data to a field scaleis associated with many uncertainties, particularlyfor heterogeneous formations. The inclusion ofcapillary pressure effects permits the interpretationof displacement tests at reservoir rates. The proposed calculation procedure extends theapproach suggested by Archer and Wong in thatthe degree of fit between observed laboratory dataand simulator results is quantified. We suggest thatrelative-permeability curves for a variety of rocktypes can be expressed in terms of two adjustable parameters and their standard error estimates.To illustrate the method, the results of displacementtests performed on cores from Swan Hills Beaverhill Lake limestone oil reservoir and Rainbow F KegRiver dolomite oil reservoir are interpreted. SPEJ P. 15^


1969 ◽  
Vol 9 (02) ◽  
pp. 221-231 ◽  
Author(s):  
R. Ehrlich ◽  
F.E. Crane

Abstract A consolidated porous medium is mathematically modeled by networks of irregularly shaped interconnected pore channels. Mechanisms are described that form residual saturations during immiscible displacement both by entire pore channels being bypassed and by fluids being isolated by the movement of an interface within individual pore channels. This latter mechanism is shown to depend strongly on pore channel irregularity. Together, these mechanisms provide an explanation for the drainage-imbibition-hysteresis effect. The calculation of steady-state relative permeabilities, based on a pore-size distribution permeabilities, based on a pore-size distribution obtained from a Berea sandstone, is described. These relative permeability curves agree qualitatively with curves that are generally accepted to be typical for highly consolidated materials. In situations where interfacial effects predominate over viscous and gravitational effects, the following conclusions are reached.Relative permeability at a given saturation is everywhere independent of flow rate.Relative permeability is independent of viscosity ratio everywhere except at very low values of wetting phase relative permeability.Irreducible wetting-phase saturation following steady-state drainage decreases with increasing ratio of nonwetting- to wetting-phase viscosity.Irreducible wetting-phase saturation following unsteady-state drainage is lower than for steady-state drainage.Irreducible nonwetting-phase saturation following imbibition is independent of viscosity ratio, whether or not the imbibition is carried out under steady- or unsteady-state conditions. Experiments qualitatively verify the conclusions regarding unsteady-state residual wetting-phase saturation. Implications of these conclusions are discussed. Introduction Natural and artificial porous materials are generally composed of matrix substance brought together in a more or less random manner. This leads to the creation of a network of interconnected pore spaces of highly irregular shape. Since the pore spaces of highly irregular shape. Since the geometry of such a network is impossible to describe, we can never obtain a complete description of its flow behavior. We can, however, abstract those properties of the porous medium pertinent to the type of flow under consideration, and thus obtain an adequate description of that flow. Thus, the Kozeny-Carmen equation, by considering a porous medium as a bundle of noninterconnecting capillary tubes, provides an adequate description of single-phase provides an adequate description of single-phase flow. With the addition of a saturation-dependent tortuosity parameter in two-phase flow to account for flow path elongation, the Kozeny-Carmen equation has been used to predict relative permeabilities for the displacement of a wetting permeabilities for the displacement of a wetting liquid by a gas. It has long been recognized that relative permeability depends not only on saturation but permeability depends not only on saturation but also on saturation history as well. Naar and Henderson described a mathematical model in which differences between drainage and imbibition behavior are explained in terms of a bypassing mechanism by which oil is trapped during imbibition. Fatt proposed a model for a porous medium that consisted of regular networks of cylindrical tubes of randomly distributed radii. From this he calculated the drainage relative permeability curves. Moore and Slobod, Rose and Witherspoon, and Rose and Cleary each considered flow in a pore doublet (a parallel arrangement of a small and pore doublet (a parallel arrangement of a small and large diameter cylindrical capillary tube). They concluded that, because of the different rates of flow in each tube, trapping would occur in one of the tubes; the extent of which would depend upon viscosity ratio and capillary pressure. SPEJ p. 221


2020 ◽  
Vol 146 ◽  
pp. 03002
Author(s):  
Marios S. Valavanides ◽  
Matthieu Mascle ◽  
Souhail Youssef ◽  
Olga Vizika

The phenomenology of steady-state two-phase flow in porous media is recorded in SCAL relative permeability diagrams. Conventionally, relative permeabilities are considered to be functions of saturation. Yet, this has been put into challenge by theoretical, numerical and laboratory studies that have revealed a significant dependency on the flow rates. These studies suggest that relative permeability models should include the functional dependence on flow intensities. Just recently a general form of dependence has been inferred, based on extensive simulations with the DeProF model for steady-state two-phase flows in pore networks. The simulations revealed a systematic dependence of the relative permeabilities on the local flow rate intensities that can be described analytically by a universal scaling functional form of the actual independent variables of the process, namely, the capillary number, Ca, and the flow rate ratio, r. In this work, we present the preliminary results of a systematic laboratory study using a high throughput core-flood experimentation setup, whereby SCAL measurements have been taken on a sandstone core across different flow conditions -spanning 6 orders of magnitude on Ca and r. The scope is to provide a preliminary proof-of-concept, to assess the applicability of the model and validate its specificity. The proposed scaling opens new possibilities in improving SCAL protocols and other important applications, e.g. field scale simulators.


Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 266
Author(s):  
Péter German ◽  
Mauricio E. Tano ◽  
Carlo Fiorina ◽  
Jean C. Ragusa

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 158-177 ◽  
Author(s):  
Pål Østebø Andersen ◽  
Yangyang Qiao ◽  
Dag Chun Standnes ◽  
Steinar Evje

Summary This paper presents a numerical study of water displacing oil using combined cocurrent/countercurrent spontaneous imbibition (SI) of water displacing oil from a water-wet matrix block exposed to water on one side and oil on the other. Countercurrent flows can induce a stronger viscous coupling than during cocurrent flows, leading to deceleration of the phases. Even as water displaces oil cocurrently, the saturation gradient in the block induces countercurrent capillary diffusion. The extent of countercurrent flow may dominate the domain of the matrix block near the water-exposed surfaces while cocurrent imbibition may dominate the domain near the oil-exposed surfaces, implying that one unique effective relative permeability curve for each phase does not adequately represent the system. Because relative permeabilities are routinely measured cocurrently, it is an open question whether the imbibition rates in the reservoir (depending on a variety of flow regimes and parameters) will in fact be correctly predicted. We present a generalized model of two-phase flow dependent on momentum equations from mixture theory that can account dynamically for viscous coupling between the phases and the porous media because of fluid/rock interaction (friction) and fluid/fluid interaction (drag). These momentum equations effectively replace and generalize Darcy's law. The model is parameterized using experimental data from the literature. We consider a water-wet matrix block in one dimension that is exposed to oil on one side and water on the other side. This setup favors cocurrent SI. We also account for the fact that oil produced countercurrently into water must overcome the so-called capillary backpressure, which represents a resistance for oil to be produced as droplets. This parameter can thus influence the extent of countercurrent production and hence viscous coupling. This complex mixture of flow regimes implies that it is not straightforward to model the system by a single set of relative permeabilities, but rather relies on a generalized momentum-equation model that couples the two phases. In particular, directly applying cocurrently measured relative permeability curves gives significantly different predictions than the generalized model. It is seen that at high water/oil-mobility ratios, viscous coupling can lower the imbibition rate and shift the production from less countercurrent to more cocurrent compared with conventional modeling. Although the viscous-coupling effects are triggered by countercurrent flow, reducing or eliminating countercurrent production by means of the capillary backpressure does not eliminate the effects of viscous coupling that take place inside the core, which effectively lower the mobility of the system. It was further seen that viscous coupling can increase the remaining oil saturation in standard cocurrent-imbibition setups.


10.2118/93-25 ◽  
1993 ◽  
Author(s):  
D. Brant Bennion ◽  
Gurk Sarioglu ◽  
Mark Chan ◽  
Toshiyuki Hirata ◽  
Dave Courtnage ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document