Research Areas

Multiphase flow in porous media
A wine stain spreading on tablecloth or oil percolating into wet sand are examples of a fluid displacing another in a porous material. The way by which fluids move and the associated enchanting invasion patterns are also scientifically and technologically important, controlling the efficiency of processes such as oil production, removal or storage of contaminant in the subsurface, or fluid mixing in microfluidics. Despite this importance and the great scientific and technological interest, modelling of how fluids flow in porous media remains elusive. A major challenge is that many of the processes occur at the scale of individual pores (sizes of the order of millimetres), whereas the phenomena of interest are often at the scale of meters and even kilometres.

Nonequilibrium flows
Multiphase and reactive fluid flow in porous media is often unstable, and highly heterogeneous: inevitable microstructural heterogeneity leads to the emergence of preferential pathways, where most of the flow is focused in a small portion of the medium. Furthermore, strong hysteresis and rate-dependence are frequently observed. This complex, out-of-equilibrium behavior is further corroborated by interactions between the fluids and the solid matrix, e.g. fracturing and dissolution. We combine experiments, simulations and theory to understand the pore-scale physics underlying nonequilibrium flows in porous media, and its implications on engineering problems, mainly in the field of energy and the environment.



Wetting-Dewetting Hysteresis

Immiscible flows in porous media involve a complex sequence of pore-scale events, from the smooth, reversible displacement of interfaces to abrupt interfacial reconfigurations and rapid pore invasion cascades (called Haines jumps). These jumps present a fundamental mechanism for irreversibility, or hysteresis, which persists at larger scales in the form hysteretic pressure-saturation relationship. Hysteresis has crucial implications on many processes such as soil wetting and drying, enhanced oil recovery and carbon geosequestration. Furthermore, pressure-saturation relation is a key component of continuum approaches to multiphase flow. Despite its importance, existing models of retention hysteresis rely on phenomenological or semi-empirical formulations.

We present here a physically-grounded (with no adjustable parameters) model of two-phase flow in disordered media, which we validate against experiments and analytical solutions. Not only that our model accounts naturally for large scale hysteresis, it also reproduces the memory properties of partial imbibition-drainage cycles and estimates of the energy dissipated in Haines jumps. This allows us to establish a direct link between pore-scale capillary instabilities and collective front displacements (avalanches), identify the mechanisms of energy dissipation in imbibition and drainage, and explain the physical origin of hysteresis and memory of partial trajectories in macroscale retention curves.

Wettability effects on invasion stability

Fluid displacement in porous media is often unstable, in the form of thin, convoluted fingers that penetrate only a small portion of the pore space. Instability associated with is unfavourable viscosity ratio--when the invading fluid is less viscous than the one displaced, e.g. gas displacing liquid--has been thoroughly investigated since the seminal Saffman-Taylor experiment (1958). Yet, the fundamental effect of wettability (the relative affinity of the fluids to the solid), that is the remarkable difference between imbibition and drainage, exposed in the classic experiments of Stokes and coworkers (1986), remain largely unexplained. Stokes et al. showed that if the invading, less viscous fluid is non-wetting (drainage), the displacement occurs through preferential flow paths leaving much of the original fluid behind, whereas if the invading fluid is more wetting (imbibition), the front becomes more compact with a more efficient displacement. We explain these classical yet intriguing observations by computer simulations, showing that wettability changes the microscopic mechanisms of pore invasion, which in turn changes the global pattern, an effect which gets weaker as the flow rate increases.

High Ca, strong drainage High Ca, weak imbibition Low Ca, strong drainage Low Ca, weak imbibition
Viscous fingering: destabilization of the entire interface, where high defending fluid pressure in the “gulfs” between fingers allows only the finger tips to advance (screening) Capillary fingering driven by disorder. Quasi-static (IP-like) advancement by invasion into the largest accesible pore along the interface Compact displacement resulting from dominance of non-local, cooperative pore-filling events ("overlaps")

For further information, see:
  • Phys. Rev. Lett. 2015 for a presentation of our novel pore-scale model, explaining the transition between the different displacement regimes.
  • Nature Sci. Rep. 2016 showing how wettability effects interacts with those of rate and pore size heterogeneity.
  • PNAS 2019 for a comprehensive comparison of pore-scale models for multiphase flow in porous media.

Effects of microstructural heterogeneity

The microstructure of a porous medium, characterized by both the distribution and the spatial arrangement of pores of various sizes, has a substantial impact on fluid displacement patterns. We combine pore-scale simulations with microfluidic experiments to quantify the impact of pore-scale heterogeneity on fluid displacement. We explain how quenched disorder (pore sizes which are spatially uncorrleated) and its interplay with rate affects drainage patterns in dedformable porous media such as bead packs (Phys. Rev. E 2010), as well as the cross-coupling with wettability (Nature Sci. Rep. 2016).

An important feature of the microstructure of many types of porous media such as soils and rocks is the existence of spatial correlations in pore sizes, such as when pores of similar size are clustered together, creating distinct regions with different hydraulic properties. We study how spatially-correlated heterogeneity affects both forced displacement in partially-wettable medium (Adv. Water Resour. 2019). We also investigated how spatial correlations in pore sizes affects isothrmal drying, where fluid displacement is driven by evaporation, numerically (Water Resour. Res. 2017) and experimentally (Nature Sci. Rep. 2017). We have used our data together with an invasion-percolation model to explain the link between microstructure and the temporal stochastic statistics (avalanches) of the pressure signal (Phys. Rev. Fluids 2018).

Fluid displacement in deformable granular media

Opening of fracture-like conduits in wet particulate materials is a phenomenon we encounter on a daily basis, for instance in drying paints or soils. It also plays a crucial role in many processes in which gas flows in the subsurface, including methane venting from lakes and marine sediments and enhanced oil recovery. We demonstrate experimentally, and explain by pore-scale simulations and scaling analysis, the emergence of fracturing as a mode of gas invasion in granular media (Phys. Rev. E 2010; Phys. Rev. Lett. 2012; see also Synopsis). We denote this regime 'capillary fracturing' because it is distinctly different from other types of fracturing. In contrast with hydraulic fracturing, which is the result of single-phase flow at high injection rates, capillary fracturing relies on interfacial phenomena (the presence of a fluid-fluid interface) and can occur at vanishingly small injection rates. Capillary fracturing allows fast exchange of mass and heat and is therefore critical to the water, carbon and energy budgets in the biosphere. One dramatic natural example is gas seepage and the creation of pockmarks in submarine sediments. The video shows our experimental reconstruction of pockmark formation in the laboratory, providing crucial insights into the unerlying mechanisms.

Reactive transport under stress

We developed a novel pore-scale model that captures coupling of dissolution and weakening-induced compaction during injection of a reactive fluid into porous media. We use numerical simulations to demonstrate how chemo-mechanical deformation of stressed rocks, namely the coupling of pore opening by dissolution with weakening-induced compaction, inhibits permeability evolution. We show that the underlying mechanism is stress concentration at the undissolved (hence stiffer) downstream region, resulting in a bottleneck effect. At high Damköhler number (Da), stress acts to reduce transport heterogeneity, promoting wormhole competition. Read more in Earth Planet. Sci. Lett. 2018.