A relationship between the fractal dimension and scaling groups of unstable miscible displacements

A detailed two-dimensional flow visualization study was performed to examine the dynamics of viscous fingering in miscible displacements. Detailed quantitative miscible displacement experiments using a microcomputer-based imaging workstation on a variety of oil recovery fluid systems were performed. The effect of two dimensionless scaling groups, namely gravity number and viscosity ratio, on the displacement behavior was investigated. Based on image analysis, the irregular fingering patterns of the flow visualization experiments were analyzed for fractal characteristics. Results indicate that the areal sweep efficiency of unstable miscible displacement follows a fractal scaling law with a fractal dimension and proportionality constant related to the gravity number and the viscosity ratio. The study shows that the fractal dimension decreases with decreasing gravity number and increasing viscosity ratio. This relationship was mapped by an artificial neural network model, which can be used to estimate the fractal dimension and the proportionality constant of miscible displacements as functions of the two scaling groups. These results have potential application in the mathematical modeling of unstable EOR displacements and in the scaling of laboratory displacements to field conditions. Received: 19 December 1999/Accepted: 26 January 2001     

Numerical investigation of viscous fingering in Hele-Shaw cell with spatially periodic variation of depth

Viscous fingering in a modified Hele-Shaw cell is numerically investigated. The cell allows periodic variation of depth in the lateral direction. The wavenumber n of the depth perturbation has great influence on fingering patterns. For n = 1, the fingering pattern due to the interface instability remains the same as that in the conventional Hele-Shaw cell, while the depth variation causes the steady finger to be a little narrower. For n = 2, four different fingering patterns are captured, similar to the available experimental observations in a modified Hele-Shaw cell containing a centered step-like occlusion. It is found that new fingering patterns appear as n further increases, among which, two patterns with spatial oscillation along both edges of the finger are particularly interesting. One is a symmetric oscillatory finger for n = 3, and the other is an asymmetric one for n = 4. The influence of capillary number on fingering patterns is studied for n = 3 and 4. We find that spatial oscillation of the finger nearly ceases at moderate capillary numbers and occurs again as the capillary number increases further. Meanwhile, the wide finger shifts to the narrow one. It is accompanied by a sudden decrease in the finger width which otherwise decreases continuously as the capillary number increases. The wavenumber and the amplitude of depth perturbation have little effect on the finger width.     

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.

     

Investigating the Temperature Dependency of Oil and Water Relative Permeabilities for Heavy Oil Systems

A look into the literature on the temperature dependency of oil and water relative permeabilities reveals contradictory reports. There are some publications reporting shifts in the water saturation range as well as variations in the relative permeability curves by temperature. On the other hand, some authors have blamed the experimental artifacts, viscous instabilities and fingering issues for these variations. We have performed core flooding experiments to further investigate this issue. Glass bead packs and sand packs were used as the porous media, and Athabasca bitumen with varying viscosities was displaced by hot water at differing temperatures. The unsteady-state method of relative permeability measurement was applied and the experimental data were history matched by a simulator that is tailor made to predict the relative permeabilities. The matches were obtained by varying the relative permeability correlation parameters. The results indicated that the initial water saturation has a direct relation with temperature, while residual oil saturation generally drops at higher temperatures. Although the water saturation range shifts, no direct and unique trend for either oil or water relative permeability is justified. The spread in relative permeabilities especially in the case of higher permeable cores suggests that viscous instabilities are present. As the same saturation shift happens by only changing the oil viscosity, the relative permeability variations with temperature can be attributed to oil to water viscosity ratio changes with temperature. Temperature dependency of relative permeabilities is more related to experimental artifacts, viscous fingering and viscosity changes than fundamental flow properties.     

Well-Posedness, Instabilities, and Bifurcation Results for the Flow in a Rotating Hele�CShaw Cell

We study the radial movement of an incompressible fluid located in a Hele–Shaw cell rotating at a constant angular velocity in the horizontal plane. Within an analytic framework, local existence and uniqueness of solutions is proved, and it is shown that the unique rotationally invariant equilibrium of the flow is unstable. There are, however, other time-independent solutions: using surface tension as a bifurcation parameter we establish the existence of global bifurcation branches consisting of stationary fingering patterns. The same results can be obtained by fixing the surface tension while varying the angular velocity. Finally, it is shown that the equilibria on a global bifurcation branch converge to a circle as the surface tension tends to infinity, provided they stay suitably bounded.     

Enhancing Immiscible Fluid Displacement in Porous Media by Capillary Pressure Discontinuities

Fluid displacement in porous media plays an important role in many industrial applications, including biological filtration, carbon capture and storage, enhanced oil recovery, and fluid transport in fuel cells. The displacement front is unstable, which evolves from smooth into ramified patterns, when the mobility (ratio of permeability to viscosity) of the displacing fluid is larger than that of the displaced one; this phenomenon is called viscous fingering. Viscous fingering increases the residual saturation of the displaced fluid, considerably impairing the efficacy of fluid displacement. It is of practical importance to develop suitable methods to improve fluid displacement. This paper presents an experimental study on applying the discontinuity of capillary pressure to improve immiscible fluid displacement in drainage for which the displacing fluid (air) wets the porous media less preferentially than does the displaced fluid (silicone oil). The concept involves using a heterogeneous packing system, where the upstream region features large pores and small capillary pressure, and the downstream region features small pores and large capillary pressure. The increase in capillary pressure prevents fingering from directly crossing the media interface, thus enhancing the displacement. The experimental apparatus was a linear cell comprising porous media between two parallel plates, and glass beads of 0.6 and 0.125 mm diameter were packed to compose the heterogeneous porous media. The time history of the finger flow was recorded using a video camera. Pressure drops over the model from the inlet to the outlet were measured to compare viscous pressure drops with capillary pressures. The results show that the fluid displacement was increased by the capillary discontinuities. The optimal displacement was determined through linear regression by adjusting the relative length of the large- and small-pore region. The results may assist in the understanding of fingering flow across the boundaries of different grain-sized bands for the gas and oil reservoir management, such as setting the relative location of the injection and production wells. The findings may also serve as a reference for industrial applications such as placing the grain bands in an adequate series to improve the displacement efficacy in biological filtration.     

粘性指进显示及其分维测量

考虑到液体驱油时粘性指进可能造成原油回采严重减少，利用图像采集系统对多孔介质中由粘性指进不稳定性导致的指进花样进行显示，测算其分维指数并研究聚合物浓度与分维的关系。     

Wetting front instabilities: a three-dimensional experimental investigation

Wetting front instability or fingering experiments were conducted in three-dimensional infiltration columns, featuring stratified fine-over-coarse texture porous media, to investigate the influences of various soil and wetting phase properties on finger diameter and propagation velocity. The system parameter varied in this study included permeability, system flux rate, media size and gradation uniformity, initial moisture content, viscosity, density, and surface tension. The influence of each parameter is discussed and compared, where applicable, to the finding of previous studies. Finger diameter and velocity data were acquired using a neutron radioscopy based, real-time imaging system. Through the use of the imaging system, a very accurate and reliable experimental data set was produced for three-dimensional fingering events.     

Numerical simulations of interfacial instabilities on a rotating miscible droplet in a time‐dependent gap Hele–Shaw cell with significant Coriolis effects

Interfacial instability of a rotating miscible droplet with significant Coriolis force in a Hele–Shaw cell is simulated numerically. The influences of the relevant control parameters are first discussed qualitatively by fingering patterns. More vigorous fingerings are found at higher rotational effects, a lower viscosity contrast and a weaker effective surface tension (Korteweg constant). For a time‐dependent gap Hele–Shaw cell, a higher cell lifting rate makes the rotating droplet bear an inward straining flow, which leads to fingering enhancement. On the contrary, a higher pressing rate provides more stable effects by additional squeezing outward flow. A quantitative analysis between the Coriolis effects and tilting angles of fingers is addressed. For arbitrary combinations of all relevant control parameters, the values of tilting angles follow a nearly linear relationship with the Coriolis effects. We estimate the correlation between the relevant control parameters (dimensionless Coriolis factor Re, viscosity parameter R, cell lifting rate a) and tilting angles (θ) of fingers that can be approximated as for significant Korteweg stresses. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Pore-Scale Investigation of Two-Phase Flows in Three-Dimensional Digital Models of Natural Sandstones

The results of numerical simulation of the processes of two-phase flow through a porous medium in three-dimensional digital models of the porous space of three natural sandstone samples are given. The calculations are carried out using the lattice Boltzmann equations and the digital field gradient model over a wide range of the capillary numbers and the viscosity ratios of injected and displaced fluids. The conditions of flow through a porous medium with capillary fingering, viscous fingering and with stable displacement front are revealed.     

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.     

Relative permeability and capillary pressure functions of porous media as related to the displacement growth pattern

Visualization experiments of the unsteady immiscible displacement of a fluid by another are performed on glass-etched pore networks of well-controlled morphology by varying the fluid system and flow conditions. The measured transient responses of the fluid saturation and pressure drop across the porous medium are introduced into numerical solvers of the macroscopic two-phase flow equations to estimate the non-wetting phase, k_rnw, and wetting phase, k_rw, relative permeability curves and capillary pressure, P_c, curve. The correlation of k_rnw, k_rw, and P_c with the displacement growth pattern is investigated. Except for the capillary number, wettability, and viscosity ratio, the immiscible displacement growth pattern in a porous medium may be governed by the shear-thinning rheology of the injected or displaced fluid, and the porous sample length as compared to the thickness of the frontal region. The imbibition k_rnw increases as the flow pattern changes from compact displacement to viscous fingering or from viscous to capillary fingering. The imbibition k_rw increases as the flow pattern changes from compact displacement or capillary fingering to viscous fingering. As the shear-thinning behaviour of the NWP strengthens and/or the contact angle decreases, then the flow pattern is gradually dominated by irregular interfacial configurations, and the imbibition k_rnw increases. The imbibition P_c is a decreasing function of the capillary number or increasing function of the injected phase viscosity in agreement with the linear thermodynamic theory.     

Miscible displacements in Hele–Shaw cells: Nonmonotonic viscosity profiles

The influence of nonmonotonic viscosity–concentration relationships on viscous fingering of neutrally buoyant, miscible fluids in a Hele–Shaw cell has been investigated. In a first step, quasisteady base states are obtained by means of nonlinear Stokes simulations. The properties of these base states are analyzed as a function of the Péclet number, the viscosity ratio, and the profile parameters. Subsequently, the stability of these base states is investigated by means of a linear stability analysis. Overall, the nonmonotonicity of the viscosity–concentration relationship is seen to have a much smaller influence on Hele–Shaw displacements than on corresponding Darcy flows. The reason for this difference lies in the nature of the respective base states. For Darcy flows, the base state is characterized by constant velocity and a diffusively decaying concentration (and hence viscosity) profile. This base viscosity profile is strongly affected by the nonmonotonicity. On the other hand, for Hele–Shaw displacements the quasisteady base states are convectively dominated and characterized by sharp fronts, so that their shape depends only weakly on the details of the viscosity–concentration relationship. Hence, for Hele–Shaw displacements both the eigenfunctions and the associated growth rates are quite similar for monotonic and nonmonotonic profiles, in contrast to the findings by [O. Manickam, G.M. Homsy, Stability of miscible displacements in porous media with nonmonotonic viscosity profiles, Phys. Fluids A 5 (1993) 1356–1367] for Darcy flows.     

Displacement of yield-stress fluids in a fracture

We consider a displacement of several yield-stress fluids in a Hele-Shaw cell. The topic is relevant to the development of a model for the flow of multiple phases inside a narrow fracture with application to hydraulically fracturing a hydrocarbon-bearing underground formation. Existing models for fracturing flows include only pure power-law models without yield stress, and the present work is aimed at filling this gap. The fluids are assumed to be immiscible and incompressible. We consider fluid advection in a plane channel in the presence of density gradients. Gravity is taken into account, so that there can be slumping and gravitational convection. We use the lubrication approximation so that governing equations are reduced to a 2D width-averaged system formed by the quasi-linear elliptic equation for pressure and transport equations for volume concentrations of fluids. The numerical solution is obtained using a finite-difference method. The pressure equation is solved using an iterative algorithm and the Multigrid method, while the transport equations are solved using a second-order TVD flux-limiting scheme with the superbee limiter. This numerical model is validated against three different sets of experiments: (i) gravitational slumping of fluids in a closed Hele-Shaw cell, (ii) viscous fingering of fluids with a high viscosity contrast due to the Saffman–Taylor (S–T) instability in a Hele-Shaw cell at microgravity conditions, (iii) displacement of Bingham fluids in a Hele-Shaw cell with the development of fingers due to the S–T instability. Good agreement is observed between simulations and laboratory data. The model is then used to investigate the joint effect of fingering and slumping. Numerical simulations show that the slumping rate of yield-stress fluid is significantly less pronounced than that of a Newtonian fluid with the same density and viscosity. If a low-viscosity Newtonian fluid is injected after a yield-stress one, the S–T instability at the interface leads to the development of fingers. As a result, fingers penetrating into a fluid with a finite yield stress locally decrease the pressure gradient and unyielded zones develop as a consequence.     

Zero-Dissipation Limit of Solutions with Shocks for Systems of Hyperbolic Conservation Laws

We consider the convergence of solutions of conservation laws with viscosity to solutions having shocks of hyperbolic conservation laws without viscosity as the viscosity tends to zero. Our analysis reveals a rich structure of nonlinear wave interactions due to the presence of shocks and initial layers. These interactions generate four different wave patterns: initial layers, shock layers, diffusion waves and coupling waves. We study the propagation and interactions of the four wave patterns by a detailed pointwise analysis. (Accepted February 19, 1998)     

Discrete particle modeling of granular Rayleigh–Taylor instability

The gravitational air–grain Rayleigh–Taylor (RT) flow instability in a Hele-Shaw cell was studied using a parallel three-dimensional discrete particle model (DPM). The onset of flow instability and the development of fingering flow structures were well captured by the model. Power spectra analysis of solid volume fraction field indicated the non-linear coarsening process of the fingering flow structures. The sensitivity of the flow patterns to the initial porosity, the Atwood number, and the ratio of particle size to the Hele-Shaw cell width was also demonstrated. The excellent agreement of DPM simulation results with the reported experimental observations proved the robustness and reliability of the numerical approach to model complex multiphase flows such as granular RT instability.     

Pore-Scale Simulations of Gas Displacing Liquid in a Homogeneous Pore Network Using the Lattice Boltzmann Method

A lattice Boltzmann high-density-ratio model, which uses diffuse interface theory to describe the interfacial dynamics and was proposed originally by Lee and Liu (J Comput Phys 229:8045–8063, 2010), is extended to simulate immiscible multiphase flows in porous media. A wetting boundary treatment is proposed for concave and convex corners. The capability and accuracy of this model is first validated by simulations of equilibrium contact angle, injection of a non-wetting gas into two parallel capillary tubes, and dynamic capillary intrusion. The model is then used to simulate gas displacement of liquid in a homogenous two-dimensional pore network consisting of uniformly spaced square obstructions. The influence of capillary number (Ca), viscosity ratio ( $M$ M ), surface wettability, and Bond number (Bo) is studied systematically. In the drainage displacement, we have identified three different regimes, namely stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number, viscosity ratio, and Bond number. Gas saturation generally increases with an increase in capillary number at breakthrough, whereas a slight decrease occurs when Ca is increased from $8.66\times 10^{-4}$ 8.66 × 10 - 4 to $4.33\times 10^{-3}$ 4.33 × 10 - 3 , which is associated with the viscous instability at high Ca. Increasing the viscosity ratio can enhance stability during displacement, leading to an increase in gas saturation. In the two-dimensional phase diagram, our results show that the viscous fingering regime occupies a zone markedly different from those obtained in previous numerical and experimental studies. When the surface wettability is taken into account, the residual liquid blob decreases in size with the affinity of the displacing gas to the solid surface. Increasing Bo can increase the gas saturation, and stable displacement is observed for $Bo>1$ B o > 1 because the applied gravity has a stabilizing influence on the drainage process.     

A boundary element method for the solution of finite mobility ratio immiscible displacement in a Hele‐Shaw cell

In this paper, the interaction between two immiscible fluids with a finite mobility ratio is investigated numerically within a Hele‐Shaw cell. Fingering instabilities initiated at the interface between a low‐viscosity fluid and a high‐viscosity fluid are analysed at varying capillary numbers and mobility ratios using a finite mobility ratio model. The present work is motivated by the possible development of interfacial instabilities that can occur in porous media during the process of CO₂ sequestration but does not pretend to analyse this complex problem. Instead, we present a detailed study of the analogous problem occurring in a Hele‐Shaw cell, giving indications of possible plume patterns that can develop during the CO₂ injection. The numerical scheme utilises a boundary element method in which the normal velocity at the interface of the two fluids is directly computed through the evaluation of a hypersingular integral. The boundary integral equation is solved using a Neumann convergent series with cubic B‐Spline boundary discretisation, exhibiting sixth‐order spatial convergence. The convergent series allows the long‐term nonlinear dynamics of growing viscous fingers to be explored accurately and efficiently. Simulations in low‐mobility ratio regimes reveal large differences in fingering patterns compared with those predicted by previous high‐mobility ratio models. Most significantly, classical finger shielding between competing fingers is inhibited. Secondary fingers can possess significant velocity, allowing greater interaction with primary fingers compared with high‐mobility ratio flows. Eventually, this interaction can lead to base thinning and the breaking of fingers into separate bubbles. Copyright © 2015 John Wiley & Sons, Ltd.     

Theoretical and experimental aspects of vacuum impregnation of porous media using transparent etched networks

Vacuum impregnation is a process method in which air and native solution are removed from the porous space of a given porous material and replaced by an external solution. Vacuum impregnation is divided into two steps: Firstly, the porous material is immersed in a liquid solution and exposed to subatmospheric pressure for a given time to ensure that air trapped in the porous materials will be removed; secondly, atmospheric pressure is re-established and the external solution penetrates the pore structure of the porous material. The objective of this study was to describe the hydrodynamic mechanisms involved in vacuum impregnation of porous materials as a function of capillary number and viscosity ratio. To achieve the objectives proposed in the present study, a transparent glass micromodel 7.7 cm × 7.4 cm was first constructed using the photolithographic technique. In addition, a stainless steel vacuum tank was built. The tank top was covered with a transparent reinforced glass plate. The whole system was connected to a vacuum pump, and a conventional video camera was adapted to record the experiments. Liquid saturation was determined through the image analysis process. Capillary number and viscosity ratio were determined for the drainage and imbibition processes. For the systems studied, we conclude that transport mechanisms ranged between stable displacement and capillary fingering during the vacuum step (drainage) while transport mechanisms ranged between continuous capillary and discontinuous capillary domains during the atmospheric step (imbibition). Earlier work indicated that our proposed process should be even more efficient for realistically large systems.