On Capillary Slowdown of Viscous Fingering in Immiscible Displacement in Porous Media

Hydrodynamic instability in immiscible porous media flows in the presence of capillarity is investigated here. The analysis and arguments presented here show that the slowdown of instabilities due to capillarity is usually very rapid which makes the flow almost, but not entirely, stable. The profiles of the stable and unstable waves in the far-field are characterized using a novel but very simple approach.     

Immiscible Displacement in Cross-Bedded Heterogeneous Porous Media

Direct insight into the mechanisms of flow and displacements within small-scale (cm) systems having permeability heterogeneities that are not parallel to the flow direction (cross-bedding and fault zones) have been carried out. In our experiments, we have used visual models with unconsolidated glass bead packs having carefully controlled permeability contrasts to observe the processes with coloured fluids and streamlines. The displacements were followed visually and by video recording for later analysis. The experiments show the significance that heterogeneities have on residual saturations and recovery, as well as the displacement patterns themselves. During a waterflood, high permeability regions can be by-passed due to capillary pressure differences, giving rise to high residual oil saturations in these regions. This study demonstrates the importance of incorporating reservoir heterogeneity into core displacement analysis, but of course the nature of the heterogeneity has to be known. In general, the effects created by the heterogeneities and their unknown boundaries hamper interpretation of flood experiments in heterogeneous real sandstone cores. Our experiments, therefore, offer clear visual information to provide a firmer understanding of the displacement processes during immiscible displacement, to present benchmark data for input to numerical simulators, and to validate the simulator through a comparison with our experimental results for these difficult flow problems.     

A New Model for Immiscible Displacement in Porous Media

Immiscible displacement is regarded as the superposition of forward flows of both water and oil, due to injection of water into the medium, and of additional forward flow of water coupled with reverse flow of oil, caused by the existence of capillary pressure gradients. The model has been evaluated numerically for the prediction of the evolution of saturation profiles in waterfloods covering a wide range of water injection rates. In agreement with experimentation, saturation profiles ranging from a completely flat shape to piston-shape, depending on the injection rate, have been obtained. Also in agreement with experimentation, numerical evaluation of the model for the case of a closed system with an initial step-function saturation profile has predicted a gradual spreading of the piston front into S-shaped profiles with an increasing variance. The final profile corresponds to uniform saturation everywhere in the medium.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

A Dynamic Model of Capillary Hysteresis in Immiscible Fluid Displacement

The capillary hysteresis in a dynamic and quasi-static two-fluid flow in a porous medium is discussed. Thermodynamic background is presented. It is shown that physically acceptable constitutive relations satisfying the thermodynamic conditions can be constructed in terms of Preisach hysteresis operators.     

Lattice Boltzmann Simulation of Immiscible Displacement in Porous Media: Viscous Fingering in a Shear-Thinning Fluid

Transport in Porous Media - In this work, we investigate immiscible displacement in porous media with the displaced fluid being shear-thinning. We focus on the influence the heterogeneous viscosity...     

Large Scale Mixing for Immiscible Displacement in Heterogeneous Porous Media

We derive a large scale mixing parameter for a displacement process of one fluid by another immiscible one in a two-dimensional heterogeneous porous medium. The mixing of the displacing fluid saturation due to the heterogeneities of the permeabilities is captured by a dispersive flux term in the large scale homogeneous flow equation. By making use of the stochastic approach we develop a definition of the dispersion coefficient and apply a Eulerian perturbation theory to determine explicit results to second order in the fluctuations of the total velocity. We apply this method to a uniform flow configuration as well as to a radial one. The dispersion coefficient is found to depend on the mean total velocity and can therefore be time varying. The results are compared to numerical multi-realization calculations. We found that the use of single phase flow stochastics cannot capture all phenomena observed in the numerical simulations.     

Effects of Precursor Wetting Films in Immiscible Displacement Through Porous Media

A computer-aided simulator of immiscible displacement in strongly water-wet consolidated porous media that takes into account the effects of the wetting films is developed. The porous medium is modeled as a three-dimensional network of randomly sized unit cells of the constricted-tube type. Precursor wetting films are assumed to advance through the microroughness of the pore walls. Two types of pore wall microroughness are considered. In the first type of microroughness, the film advances quickly, driven by capillary pressure. In the second type, the meniscus moves relatively slowly, driven by local bulk pressure differences. In the latter case, the wetting film often forms a collar that squeezes the thread of oil causing oil disconnection. Each pore is assumed to have either one of the aforementioned microroughness types, or both. The type of microroughness in each pore is assigned randomly. The simulator is used to predict the residual oil saturation as a function of the pertinent parameters (capillary number, viscosity ratio, fraction of pores with each type of wall microroughness). These results are compared with those obtained in the absence of wetting films. It is found that wetting films cause substantial increase of the residual oil saturation. Furthermore, the action of the wetting films causes an increase of the mean volume of the residual oil ganglia.     

LBM Simulation of Viscous Fingering Phenomenon in Immiscible Displacement of Two Fluids in Porous Media

This article presents the lattice Boltzmann simulation of viscous fingering phenomenon in immiscible displacement of two fluids in porous media. Such phenomenon generally takes place when a less viscous fluid is used to displace a more viscous fluid, and it can be found in many industrial fields. Dimensionless quantities, such as capillary number, Bond number and viscosity ratio between displaced fluid and displacing fluid are introduced to illustrate the effects of capillary force, viscous force, and gravity on the fluid behaviour. The surface wettability, which has an impact on the finger pattern, is also considered in the simulation. The numerical procedure is validated against the experiment about viscous fingering in a Hele-Shaw cell. The displacement efficiency is investigated using the parameter, areal sweep efficiency. The present simulation shows an additional evidence to demonstrate that the lattice Boltzmann method is a useful method for simulating some multiphase flow problems in porous media.     

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.     

Modeling and Numerical Simulations of Immiscible Compressible Two-Phase Flow in Porous Media by the Concept of Global Pressure

A new formulation is presented for the modeling of immiscible compressible two-phase flow in porous media taking into account gravity, capillary effects, and heterogeneity. The formulation is intended for the numerical simulation of multidimensional flows and is fully equivalent to the original equations, contrary to the one introduced in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986). The main feature of this formulation is the introduction of a global pressure. The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic (the global pressure equation) and a nonlinear diffusion–convection one (the saturation equation) which can be efficiently solved numerically. A finite volume method is used to solve the global pressure equation and the saturation equation for the water and gas phase in the context of gas migration through engineered and geological barriers for a deep repository for radioactive waste. Numerical results for the one-dimensional problem are presented. The accuracy of the fully equivalent fractional flow model is demonstrated through comparison with the simplified model already developed in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986).     

Scaling Immiscible Displacements in Porous Media with Horizontal Wells

Despite the increase in horizontal well applications, scaling fluid displacement in porous medium with horizontal wells is yet to be fully investigated. Determining the conditions under which horizontal wells may lead to better oil recovery is of great importance to the petroleum industry. In this paper, a numerical sensitivity study was performed for several well configurations. The study is performed in order to reveal the functional relationships between the scaling groups governing the displacement and the performance of immiscible displacements in homogeneous reservoirs produced by horizontal wells. These relationships can be used as a quick prediction tool for the fractional oil recovery for any combinations of the scaling groups, thus eliminating the need for the expensive fine-mesh simulations. In addition, they provide the condition under which a horizontal well configuration may yield better recovery performance. These results have potential applications in modeling immiscible displacements and in the scaling of laboratory displacements to field conditions.     

Immiscible Displacement of a Wetting Fluid by a Non-wetting One at High Capillary Number in a Micro-model Containing a Single Fracture

Most reservoirs in Iran are heterogeneous fractured carbonate reservoirs. Heterogeneity causes an earlier breakthrough and an unstable front which leads to a lower recovery. A series of experiments were conducted whereby the distilled water displaced n-Decane in strongly oil-wet glass micro-models containing a single fracture. Experimental data from image analysis of immiscible displacement processes are used to modify the Buckley?CLeverett and fractional flow equations by a heterogeneity factor. It is shown that the heterogeneity factor in the modified equations can be expressed as a function of fracture length and orientation.     

A Monte Carlo Algorithm for Immiscible Two-Phase Flow in Porous Media

We present a Markov Chain Monte Carlo algorithm based on the Metropolis algorithm for simulation of the flow of two immiscible fluids in a porous medium under macroscopic steady-state conditions using a dynamical pore network model that tracks the motion of the fluid interfaces. The Monte Carlo algorithm is based on the configuration probability, where a configuration is defined by the positions of all fluid interfaces. We show that the configuration probability is proportional to the inverse of the flow rate. Using a two-dimensional network, advancing the interfaces using time integration, the computational time scales as the linear system size to the fourth power, whereas the Monte Carlo computational time scales as the linear size to the second power. We discuss the strengths and the weaknesses of the algorithm.     

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.     

Modelling the Effect of Surface Tension on Immiscible Displacement in a Saturated Porous Medium

Comments are made about the model employed by Chen and Vafai for forced convection in a porous medium channel, with a surface tension effect at the moving interface between two fluids, when one fluid is displaced by the other. A simple situation is analysed, and the circumstances under which surface tension effects are important in this case are clarified.     

Relations Between Seepage Velocities in Immiscible,Incompressible Two-Phase Flow in Porous Media

Based on thermodynamic considerations, we derive a set of equations relating the seepage velocities of the fluid components in immiscible and incompressible two-phase flow in porous media. They necessitate the introduction of a new velocity function, the co-moving velocity. This velocity function is a characteristic of the porous medium. Together with a constitutive relation between the velocities and the driving forces, such as the pressure gradient, these equations form a closed set. We solve four versions of the capillary tube model analytically using this theory. We test the theory numerically on a network model.     

An Ising-Based Simulator for Capillary Action in Porous Media

Multiphase flows in porous media are encountered in several contexts—e.g., hydrocarbon recovery operations, battery electrodes, microfluidic devices, etc. Capillary-dominated flows are interesting due to the complex interplay of interfacial properties and pore geometries. Conventional hydrodynamic flow solvers are computationally inefficient in the capillary-dominated regime, particularly in complex pore structures. The algorithm developed here specifically targets this regime to reduce simulation times. We minimise the fluid–fluid and fluid–solid interaction energies through an approach inspired by the ferromagnetic Ising model. We validate the algorithm on (1) model pore geometries with analytical solutions for capillary action, and (2) rocks with available mercury porosimetry data. We validate its predictions for model geometries and sandstones using (1) curvatures calculated from theories developed by Mayer–Stowe–Princen, Ma and Morrow, and Mason and Morrow; (2) predictions from GeoDict, a commercial software package, which also includes a state-of-the-art drainage simulator; (3) mercury porosimetry data. Drainage capillary pressure curves predicted for Bentheimer and Fontainebleau rocks reasonably match porosimetry data.     

Capillary Heterogeneity Trapping and Crossflow in Layered Porous Media

We examine the effect of capillary and viscous forces on the displacement of one fluid by a second, immiscible fluid across and along parallel layers of contrasting porosity, and relative permeability, as well as previously explored contrasts in absolute permeability and capillary pressure. We consider displacements with wetting, intermediate-wetting and non-wetting injected phases. Flow is characterized using six independent dimensionless numbers and a dimensionless storage efficiency, which is numerically equivalent to the recovery efficiency. Results are directly applicable to geologic carbon storage and hydrocarbon production. We predict how the capillary–viscous force balance influences storage efficiency as a function of a small number of key dimensionless parameters, and provide a framework to support mechanistic interpretations of complex field or experimental data, and numerical model predictions, through the use of simple dimensionless models. When flow is directed across layers, we find that capillary heterogeneity traps the non-wetting phase, regardless of whether it is the injected or displaced phase. However, minimal trapping occurs when the injected phase is intermediate-wetting or when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. A dimensionless capillary-to-viscous number defined using the layer thickness rather than the more commonly used system length is most relevant to predict capillary heterogeneity trapping. When flow is directed along layers, we show that, regardless of wettability, increasing capillary crossflow reduces the distance between the leading edges of the injected phase in each layer and increases storage efficiency. This may be counter-intuitive when the injected phase is non-wetting. Crossflow has a significant impact on storage efficiency only when high-permeability layers contain a smaller moveable volume of fluid than low-permeability layers. In that case, capillary heterogeneity traps the wetting phase, regardless of whether it is the injected or displaced phase.