### Abstract

Language | English |
---|---|

Article number | 052103 |

Number of pages | 17 |

Journal | Physics of Fluids |

Volume | 27 |

Issue number | 5 |

Early online date | 27 May 2015 |

DOIs | |

State | Published - 31 May 2015 |

### Fingerprint

### Keywords

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

### Cite this

*Physics of Fluids*,

*27*(5), [052103]. DOI: 10.1063/1.4921611

}

*Physics of Fluids*, vol. 27, no. 5, 052103. DOI: 10.1063/1.4921611

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

Research output: Contribution to journal › Article

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 -