Lattice Boltzmann simulation of immiscible fluid displacement in porous media: homogeneous versus heterogeneous pore network

Haihu Liu, Yonghao Zhang, Albert J. Valocchi

Research output: Contribution to journalArticle

  • 28 Citations

Abstract

Injection of anthropogenic carbon dioxide (CO2) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO2 that can be captured and its long term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely stable displacement, capillary fingering and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (Snw) is found to increase nearly linearly with log Cafor each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in Snw. In either pore network, the specific interfacial length is linearly proportional to Snw during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on Snw for either favorable (M≥1) or unfavorable (M<1) displacement, and the slope is slightly higher for the unfavorable displacement.
LanguageEnglish
Article number052103
Number of pages17
JournalPhysics of Fluids
Volume27
Issue number5
Early online date27 May 2015
DOIs
StatePublished - 31 May 2015

Fingerprint

Immiscible Fluids
Lattice Boltzmann
Porous Media
Porous materials
porosity
Fluids
fluids
Wetting
Simulation
simulation
Fluid
Viscosity
wetting
Gas emissions
Greenhouse gases
Microfluidics
Carbon Dioxide
Drainage
Linearly
Carbon dioxide

Keywords

  • pore-scale simulations
  • lattice Boltzmann model
  • porous media
  • multiphase flow
  • fingering
  • heterogeneity

Cite this

@article{c9bf520d21584194afe48df7dcae6fa1,
title = "Lattice Boltzmann simulation of immiscible fluid displacement in porous media: homogeneous versus heterogeneous pore network",
abstract = "Injection of anthropogenic carbon dioxide (CO2) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO2 that can be captured and its long term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely stable displacement, capillary fingering and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (Snw) is found to increase nearly linearly with log Cafor each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in Snw. In either pore network, the specific interfacial length is linearly proportional to Snw during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on Snw for either favorable (M≥1) or unfavorable (M<1) displacement, and the slope is slightly higher for the unfavorable displacement.",
keywords = "pore-scale simulations, lattice Boltzmann model, porous media, multiphase flow, fingering, heterogeneity",
author = "Haihu Liu and Yonghao Zhang and Valocchi, {Albert J.}",
year = "2015",
month = "5",
day = "31",
doi = "10.1063/1.4921611",
language = "English",
volume = "27",
journal = "Physics of Fluids",
issn = "1070-6631",
number = "5",

}

Lattice Boltzmann simulation of immiscible fluid displacement in porous media : homogeneous versus heterogeneous pore network. / Liu, Haihu; Zhang, Yonghao; Valocchi, Albert J.

In: Physics of Fluids, Vol. 27, No. 5, 052103, 31.05.2015.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Lattice Boltzmann simulation of immiscible fluid displacement in porous media

T2 - Physics of Fluids

AU - Liu,Haihu

AU - Zhang,Yonghao

AU - Valocchi,Albert J.

PY - 2015/5/31

Y1 - 2015/5/31

N2 - Injection of anthropogenic carbon dioxide (CO2) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO2 that can be captured and its long term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely stable displacement, capillary fingering and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (Snw) is found to increase nearly linearly with log Cafor each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in Snw. In either pore network, the specific interfacial length is linearly proportional to Snw during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on Snw for either favorable (M≥1) or unfavorable (M<1) displacement, and the slope is slightly higher for the unfavorable displacement.

AB - Injection of anthropogenic carbon dioxide (CO2) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO2 that can be captured and its long term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely stable displacement, capillary fingering and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (Snw) is found to increase nearly linearly with log Cafor each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in Snw. In either pore network, the specific interfacial length is linearly proportional to Snw during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on Snw for either favorable (M≥1) or unfavorable (M<1) displacement, and the slope is slightly higher for the unfavorable displacement.

KW - pore-scale simulations

KW - lattice Boltzmann model

KW - porous media

KW - multiphase flow

KW - fingering

KW - heterogeneity

UR - http://scitation.aip.org/content/aip/journal/pof2

U2 - 10.1063/1.4921611

DO - 10.1063/1.4921611

M3 - Article

VL - 27

JO - Physics of Fluids

JF - Physics of Fluids

SN - 1070-6631

IS - 5

M1 - 052103

ER -