Immiscible Displacement in the Interacting Capillary Bundle Model Part I. Development of Interacting Capillary Bundle Model

An interacting capillary bundle model is developed for analysing immiscible displacement processes in porous media. In this model, pressure equilibration among the capillaries is stipulated and capillary forces are included. This feature makes the model entirely different from the traditional tube bundle models in which fluids in different capillaries are independent of each other. In this work, displacements of a non-wetting phase by a wetting phase at different injection rates were analysed using the interacting capillary bundle model. The predicted evolutions of saturation profiles were consistent with both numerical simulation and experimental results for porous media reported in literature which cannot be re-produced with the non-interacting tube bundle models.     

Theoretical and Experimental Investigations of Miscible Displacement in Fractured Porous Media

In this paper a mathematical model for miscible displacement in fractured porous media is developed. The model takes into account mechanisms of mass transfer between fracture and matrix. The model is normalized by using the dimensionless parameters, which characterize the process, and the analytical solutions of the resulting system of equations are provided by utilizing the method of characteristics. For comparison the results of model with experimental results, laboratory displacement tests have been performed in fractured systems under miscible displacement. The porous media used were cylindrical Asmari cores from Iranian reservoirs containing an artificially vertical fracture. Normal heptane and kerosene were two miscible fluid used. There is very good agreement between experiments and model prediction.     

Review of Steady-State Two-Phase Flow in Porous Media: Independent Variables,Universal Energy Efficiency Map,Critical Flow Conditions,Effective Characterization of Flow and Pore Network

In many applications of two-phase flow in porous media, a wetting phase is used to displace through a network of pore conduits as much as possible of a non-wetting phase, residing in situ. The energy efficiency of this physical process may be assessed by the ratio of the flow rate of the non-wetting phase over the total mechanical power externally provided and irreversibly dissipated within the process. Fractional flow analysis, extensive simulations implementing the DeProF mechanistic model, as well as a recent retrospective examination of laboratory studies have revealed universal systematic trends of the energy efficiency in terms of the actual independent variables of the process, namely the capillary number, Ca, and the flow rate ratio, r. These trends can be cast into an energy efficiency map over the (Ca, r) domain of independent variables. The map is universal for all types of non-wetting/wetting phase porous medium systems. It demarcates the efficiency of steady-state two-phase flow processes in terms of pertinent system parameters. The map can be used as a tool for designing more efficient processes, as well as for the normative characterization of two-phase flows, as to the predominance of capillary or viscous effects. This concept is based on the existence of a unique locus of critical flow conditions, for which the energy efficiency takes locally maximum values. The locus shape depends on the physicochemical characteristics of the non-wetting phase/wetting phase/porous medium system, and it shows a significant mutation as the externally imposed flow conditions change the type of flow, from capillary- to viscosity-dominated. The locus can be approached by an S-type functional form in terms of the capillary number and the system properties (viscosity ratio, wettability, pore network geometry, etc.), suggesting that formative criteria can be derived for flow characterization in any system. A new, extended definition of the capillary number is also proposed that effectively takes into account the critical properties of all the system constituents. When loci of critical flow conditions pertaining to processes with different viscosity ratio in the same pore network, are expressed in terms of this true-to-mechanism capillary number, they collapse into a unique locus. In this context, a new methodology for the effective characterization of pore networks is proposed.     

Unstable Displacement of Non-aqueous Phase Liquids with Surfactant and Polymer

In this paper, we study two-phase multicomponent displacement of two immiscible fluids in both homogeneous and heterogeneous porous media. In many applications such as enhanced oil recovery, fluid mixing and spreading can be detrimental to the efficacy of the process. Here, we show that when an initially immobile phase is being displaced by a finite-size slug of solvents (surfactant and polymer), viscous fingering significantly enhances mixing and spreading of solvents. These effects are similar to those caused by medium heterogeneity and lead to poor displacement efficiency. We first quantify the displacement efficiency subject to different mobility ratios, Peclet numbers, and levels of medium heterogeneity. We observe a non-monotonic behavior in displacement efficiency as a function of mobility ratio, indicating that although stable frontal interface is desirable, miscible viscous fingering on the rear interface will eventually disintegrate the solvents slugs and reduce the displacement efficiency. Then, we show that miscible viscous fingering developing on the rear interface of the chemical slug could be greatly suppressed when viscosity contrast is gradually decreased using exponential or linear functions, leading to 10% increase in displacement efficiency while using the same amount of chemicals. To elucidate this low displacement efficiency, we study the evolution of mixing, spreading, and interfacial length and show that while higher viscosity ratios are quite effective in mobilizing the initially immobile phase in 1D displacements, they are in fact detrimental in 2D unstable displacements since they enhance mixing and spreading of solvents.

     

Pore-scale Simulation of Water Alternate Gas Injection

We use a three-dimensional mixed-wet random network model representing Berea sandstone to compute displacement paths and relative permeabilities for water alternating gas (WAG) flooding. First we reproduce cycles of water and gas injection observed in previously published experimental studies. We predict the measured oil, water and gas relative permeabilities accurately. We discuss the hysteresis trends in the water and gas relative permeabilities and compare the behavior of water-wet and oil-wet media. We interpret the results in terms of pore-scale displacements. In water-wet media the water relative permeability is lower during water injection in the presence of gas due to an increase in oil/water capillary pressure that causes a decrease in wetting layer conductance. The gas relative permeability is higher for displacement cycles after first gas injection at high gas saturation due to cooperative pore filling, but lower at low saturation due to trapping. In oil-wet media, the water relative permeability remains low until water-filled elements span the system at which point the relative permeability increases rapidly. The gas relative permeability is lower in the presence of water than oil because it is no longer the most non-wetting phase.     

Mixing and Mass Transfer in Multicontact Miscible Displacements

In this work, we investigate the accuracy of some physical models that are frequently used to describe and interpret dispersive mixing and mass transfer in compositional reservoir simulation. We have designed a quaternary analog fluid system (alcohol?Cwater?Chydrocarbon) that mimics the phase behavior of CO₂-hydrocarbon mixtures at high pressure and temperature. A porous medium was designed using PolyTetraFlouroEthylene (PTFE) materials to ensure that the analog oil acts as the wetting phase, and the properties of the porous medium were characterized in terms of porosity, permeability and dispersivity. Relative permeability and interfacial tension (IFT) measurements were also performed to delineate interactions between the fluid system and the porous medium. The effluent concentrations from two-component first-contact miscible (FCM) displacement experiments exhibit a tailing behavior that is attributed to imperfect sweep of the porous medium: A feature that is not captured by normal dispersion models. To represent this behavior in displacement calculations, we use dual-porosity (DP) models including mass transfer between flowing and stagnant porosities. Two 4-component two-phase displacement experiments were performed at near-miscible and multicontact miscible (MCM) conditions and the effluent concentrations were interpreted by numerical calculations. We demonstrate that the accuracy of our displacement calculations relative to the experimental observations is sensitive to the selected models for dispersive mixing, mass transfer between flowing and stagnant porosities, and IFT scaling of relative permeability functions. We also demonstrate that numerical calculations substantially agree with the experimental observations for some physical models with limited need for model parameter adjustment. The combined experimental and modeling effort presented in this work identifies and explores the impact of a set of physical mechanisms (dispersion and mass transfer) that must be upscaled adequately for field-scale displacement calculations in DP systems.     

Radial viscous fingering patterns in Hele-Shaw cells

Viscous fingering experiments were performed by injecting a liquid to radially displace a much more viscous liquid in a Hele-Shaw cell consisting of two parallel closely-spaced glass plates. Both smooth and etched plates were used to study the influence of plate roughness on the fingering mechanism. Effect of flow rate and interfacial tension was also demonstrated. The results show that symmetric dendritic finger patterns can form in the presence of anisotropy provided by an etched square network, for both miscible and immiscible fingers. Chaotic finger patterns can form both in a cell with smooth surfaces and in one having a network of randomly oriented channels etched on one plate. Due to interfacial tension, the immiscible finger patterns are less ramified than their miscible counterparts, are more sensitive to the flow rate and become compact as the flow rate decreases. Possible applications of two-phase displacement studies in Hele-Shaw cells are discussed, which include two-phase flow in porous media and acid fracturing of oil reservoirs.     

STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

IntroductionMisciblefloodingisadriveprocessbymixinginjectionfluid (solvent)andoil.Itsmainmechanismistodecreasetheresidualoilsaturationbyeliminatinginterfacialtensionbetweenphases.GasdrivehasanincreasingpercentageofEORprojectsinU .S .A .,Canadaandsomeothercountriesyearsbyyears.From 1960’s ,carbondioxideinjectionhasbeingstudiedinDaqingoilfieldofChina ,buttheprocessisslowforlackinggasresource .Inrecentyears,withtheinconsistentinreserve_productionequilibriumbecomingmoreseriousandthediscoveryo…     

Efficiency of Capillary Imbibition Dominated Displacement of Nonwetting Phase by Wetting Phase in Fractured Porous Media

The critical and optimum injection rates as well as the critical fracture capillary number for an efficient displacement process are determined based on the experimental and numerical modeling of the displacement of nonwetting phase (oil) by wetting phase (water) in fractured porous media. The efficiency of the process is defined in terms of the nonwetting phase displaced from the system per amount of wetting phase injected and per time. Also, the effects of injection rate on capillary imbibition transfer dominated two-phase flow in fractured porous media are clarified by visualizing the experiments. The results reveal that as the injection rate is increased, fracture pattern begins to become an effective parameter on the matrix saturation distribution. As the rate is lowered, however, the system begins to behave like a homogeneous system showing a frontal displacement regardless the fracture configuration.     

Numerical study of non-uniqueness of the factors influencing relative permeability in heterogeneous porous media by lattice Boltzmann method

The effects of capillary number Ca and viscosity ratio M on non-uniqueness of the relative permeability (RP) – wetting saturation (S_w) relationships of steady-state immiscible two-phase flow in the heterogeneous porous media are studied by a developed lattice Boltzmann method (LBM). In this work, it is suggested that the non-uniqueness of capillary number Ca and viscosity ratio M influencing RP–S_w curves is due to the presence of the heterogeneity of porous media. For the different Ca and M, the RP–S_w curves in the homogenous and heterogeneous porous media with the same porosity are discussed in detail. It is found that (1) the pore-size distributions in porous media have significant effect on the RPs of both wetting phase (WP) and non-wetting phase (NWP); (2) the heterogeneity of porous media has negative effect on the increment of the RP of NWP caused by the large capillary number. The lubricating effect is strongly dependent on pore-size distributions in porous media; and (3) the large viscosity ratio increases the RP of NWP, while it has little effect on the RP of WP. The effect on the RP of NWP may be magnified by the heterogeneity of porous media.     

Multiresolution Wavelet Scale Up of Unstable Miscible Displacements in Flow Through Heterogeneous Porous Media

We describe scale up of geological models of field-scale porous media using a new method based on the wavelet transformations. The porous media of interest contain broadly-distributed and correlated permeabilities. Wavelet transformation of the permeability field of such porous media coarsens the geological model from smallest to largest length scales, drastically reduces the number of equations to be solved, preserves the important information on the permeability field at all the relevant length scales, and yields numerical results for any fluid flow property that are as accurate as those that are obtained with the highly detailed geological model of the same porous media. To test this method, we carry out extensive computer simulations of unstable miscible displacement processes and the associated viscous fingering phenomenon in highly heterogeneous porous media, both with the fine-scale geological model and the coarsened model. Excellent agreement is found between the results of the two sets of simulations.     

STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

According to the research theory of improved black oil simulator, a practical mathematical model for C02 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance fast, streamline method is employed to solve this model as an alternative to traditional finite difference methods. Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method (BEM), an explicit total variation diminishing (TVD) method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less information need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an effective tool for early stage screening andmiscible oil field.reservoir dynamic management of the CO2 miscible oil field.     

Method of Negative Saturations for Modeling Two-phase Compositional Flow with Oversaturated Zones

We examine the two-phase flow through porous media of multicomponent partially miscible fluids. The composition of both the phases is variable in space and time and is assumed to be in local thermodynamic equilibrium. One of the basic problems in modeling such systems is related to the appearance of single-phase zones occupied by the fluid which is over- (or under-) saturated, i.e., it is significantly remote from the equilibrium two-phase region. In an oversaturated zone, the two-phase flow equations degenerate and can no longer be used, which provokes serious numerical problems. We propose to describe the two-phase and oversaturated single-phase zone by a uniform system of classic two-phase equations while extending the concept of the phase saturation so that it may be negative and higher than one. Physically this means that we consider the oversaturated single-phase states as the pseudo two-phase states which are characterized by a negative saturation of the imaginary phase. Such an extension of the concept of the phase saturation requires the development of some consistence conditions that ensure the equivalence between the pseudo two-phase equations and the true single-phase flow model in the oversaturated zones. This method allows using the existing numerical simulators of two-phase flow for modeling single-phase zones by adding a simple plug-in with no modification of the structure of the simulators. The method is illustrated by several examples of hydrogene-water flow in a waste radioactive storage and of CO₂ injection in an oil reservoir.     

Interfaces of Phase Transition and Disappearance and Method of Negative Saturation for Compositional Flow with Diffusion and Capillarity in Porous Media

The specific case of interfaces separating a single-phase fluid and a two-phase continuum appears in the theory of compositional flow through porous media. They are usually called the interfaces of phase transition (PT-interfaces) or the interfaces of phase disappearing (PD-interfaces). The principle of equivalence is proved which shows that a single-phase multi-component fluid may be replaced by an equivalent fictitious two-phase fluid having specific properties. The equivalent properties are such that the extended saturation of a fictitious phase is negative. This principle enables us to develop the uniform system of two-phase equations in the overall domain in terms of the extended saturation (the NegSat model), and to apply the direct numerical simulation. In the case with diffusion, the uniform NegSat model contains a new term proportional to the gradient of saturation in the relation for flow velocity. The canonical NegSat model represents a transport equation with discontinuous nonlinearities. The qualitative analysis of this model shows that the PT-interfaces represent the shocks of the extended saturation, or, in some cases, can transform into weak shocks. The diffusion and capillarity do not destroy necessarily the shocks, but change their velocity. The analytical technique is developed which allows capturing PT-shocks. The method is illustrated by several examples of miscible gas injection in oil reservoir. In two-dimensional case, the effects of multiple shock collisions in heterogeneous media are automatically modeled. In the case of immiscible fluids and a classic interface, the suggested method converges to the VOF-method.     

Analytical Solutions to Multiphase First-Contact Miscible Models with Viscous Fingering

In this paper, we analyze an empirical model of viscous fingering for three-component, two-phase, first-contact miscible flows. We present the complete range of analytical solutions to secondary and tertiary water-alternating-gas (WAG) floods. An important ingredient in the construction of analytical solutions is the presence of detached (nonlocal) branches of the Hugoniot locus, that is, curves in composition space that satisfy the Rankine–Hugoniot conditions but do not contain the reference state. We illustrate how, in water–solvent floods into a medium with mobile water and residual oil (immobile to water), the solvent front and the water Buckley–Leverett front may interact, resulting in a leading water/solvent shock that is stable to viscous fingering. The analytical solutions explain why in these miscible tertiary floods, oil and solvent often break through simultaneously. We discuss the implications of the new solutions in the design of miscible tertiary floods, such as the estimation of the optimum WAG ratio.     

Miscible displacement from fractured porous media

A model of miscible displacement of incompressible fluids from a fractured porous medium is proposed. The model describes the process of displacement of oil by solvents, the cycling process of displacement of aliphatic hydrocarbon gas by dry gas at low repressions on the formation, and other processes of single-phase multicomponent displacement from fractured porous media. Problems relating to the pumping of a neutral admixture and admixture slugs through a fractured porous reservoir are solved.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 100–110, November–December, 1989.The authors are grateful to K. S. Basniev, A. K. Kurbanov, V. I. Maron, and M. I. Shvidler for useful discussions.     

Simulation of Single-Phase Multicomponent Flow Problems in Gas Reservoirs by Eulerian–Lagrangian Techniques

Over the past two decades most discussions of the simulation of miscible displacement in porous media were related to incompressible flow problems; recently, however, attention has shifted to compressible problems. The first goal of this paper is the derivation of the governing equations (mathematical models) for a hierarchy of miscible isothermal displacements in porous media, starting from a very general single-phase, multicomponent, compressible flow problem; these models are then compared with previously proposed models. Next, we formulate an extension of the modified method of characteristics with adjusted advection to treat the transport and dispersion of the components of the miscible fluid; the fluid displacement must be coupled in a two-stage operator-splitting procedure with a pressure equation to define the Darcy velocity field required for transport and dispersion, with the outer stage incorporating an implicit solution of the nonlinear parabolic pressure equation and an inner stage for transport and diffussion in which the mass fraction equations are solved sequentially by first applying a globally conservative Eulerian–Lagrangian scheme to solve for transport, followed by a standard implicit procedure for including the diffusive effects. The third objective is a careful investigation of the underlying physics in compressible displacements in porous media through several high resolution numerical experiments. We consider real binary gas mixtures, with realistic thermodynamic correlations, in homogeneous and heterogeneous formations.     

考虑热流固耦合作用的多孔介质孔隙尺度两相流动模拟

为了研究深层油气资源在岩石多孔介质内的运移过程, 使用一种基于Darcy-Brinkman-Biot的流固耦合数值方法, 结合传热模型, 完成了Duhamel-Neumann热弹性应力的计算, 实现了在孔隙模拟多孔介质内的考虑热流固耦合作用的两相流动过程. 模型通过求解Navier-Stokes方程完成对孔隙空间内多相流体的计算, 通过求解Darcy方程完成流体在岩石固体颗粒内的计算, 二者通过以动能方式耦合的形式, 计算出岩石固体颗粒质点的位移, 从而实现了流固耦合计算. 在此基础上, 加入传热模型考虑温度场对两相渗流过程的影响. 温度场通过以产生热弹性应力的形式作用于岩石固体颗粒, 总体上实现热流固耦合过程. 基于数值模型, 模拟油水两相流体在二维多孔介质模型内受热流固耦合作用的流动过程. 研究结果表明: 热应力与流固耦合作用产生的应力方向相反, 使得总应力比单独考虑流固耦合作用下的应力小; 温度的增加使得模型孔隙度增加, 但当注入温差达到150 K后, 孔隙度不再有明显增加; 温度的增加使得水相的相对渗流能力增加, 等渗点左移.