skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation

Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Hydrology
Additional Journal Information:
Journal Volume: 522; Journal Issue: C; Related Information: CHORUS Timestamp: 2016-09-04 14:09:26; Journal ID: ISSN 0022-1694
Country of Publication:

Citation Formats

Aghaei, Arash, and Piri, Mohammad. Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation. Netherlands: N. p., 2015. Web. doi:10.1016/j.jhydrol.2015.01.004.
Aghaei, Arash, & Piri, Mohammad. Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation. Netherlands. doi:10.1016/j.jhydrol.2015.01.004.
Aghaei, Arash, and Piri, Mohammad. 2015. "Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation". Netherlands. doi:10.1016/j.jhydrol.2015.01.004.
title = {Direct pore-to-core up-scaling of displacement processes: Dynamic pore network modeling and experimentation},
author = {Aghaei, Arash and Piri, Mohammad},
abstractNote = {},
doi = {10.1016/j.jhydrol.2015.01.004},
journal = {Journal of Hydrology},
number = C,
volume = 522,
place = {Netherlands},
year = 2015,
month = 3

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jhydrol.2015.01.004

Citation Metrics:
Cited by: 11works
Citation information provided by
Web of Science

Save / Share:
  • Permeability contrasts exist in multilayer geological formations under consideration for carbon sequestration. To improve our understanding of heterogeneous pore-scale displacements, liquid CO2 (LCO2) - water displacement was evaluated in a pore network micromodel with two distinct permeability zones. Due to the low viscosity ratio (logM = -1.1), unstable displacement occurred at all injection rates over two orders of magnitude. LCO2 displaced water only in the high permeability zone at low injection rates with the mechanism shifting from capillary fingering to viscous fingering with increasing flow rate. At high injection rates, LCO2 displaced water in the low permeability zone with capillarymore » fingering as the dominant mechanism. LCO2 saturation (SLCO2) as a function of injection rate was quantified using fluorescent microscopy. In all experiments, more than 50% of LCO2 resided in the active flowpaths, and this fraction increased as displacement transitioned from capillary to viscous fingering. A continuum-scale two-phase flow model with independently determined fluid and hydraulic parameters was used to predict SLCO2 in the dual-permeability field. Agreement with the micromodel experiments was obtained for low injection rates. However, the numerical model does not account for the unstable viscous fingering processes observed experimentally at higher rates and hence overestimated SLCO2.« less
  • Carbon sequestration in saline aquifers involves displacing resident brine from the pore space by supercritical CO 2 (scCO 2 ). The displacement process is considered unstable due to the unfavorable viscosity ratio (logM < 0). The unstable mechanisms that affect scCO 2 - water displacement under reservoir conditions (i.e., 41 °C, 9 MPa) were investigated in a homogeneous micromodel. A wide range of injection rates (logCa = -7.61~-4.73) was studied in two sets of experiments: discontinuous-rate injection, where the micromodel was first cleaned and saturated with water before each injection rate was imposed, and continuous-rate injection, where the rate wasmore » increased after quasi-steady conditions were reached for a certain rate. For the discontinuous-rate experiments, capillary fingering and viscous fingering are the dominant mechanisms for low (logCa <= -6.61) and high injection rates (logCa >= -5.21), respectively. Crossover from capillary to viscous fingering was observed for logCa = -5.91~-5.21, resulting in a large decrease in scCO 2 saturation. The discontinuous-rate experimental results confirmed the decrease in nonwetting fluid saturation during crossover from capillary to viscous fingering predicted by numerical simulations by Lenormand et al. (1988).1 Capillary fingering was the only mechanism that dominates all injection rates in the continuous-rate experiment, and resulted in monotonic increase in scCO 2 saturation.« less
  • Injection of anthropogenic carbon dioxide (CO{sub 2}) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO{sub 2} 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,more » 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 (S{sub nw}) is found to increase nearly linearly with logCa for 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 S{sub nw}. In either pore networks, the specific interfacial length is linearly proportional to S{sub nw} 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 S{sub nw} for either favorable (M ≥ 1) or unfavorable (M < 1) displacement, and the slope is slightly higher for the unfavorable displacement.« less
  • This work discusses the controversy, found in the published literature, over the microscopic displacement mechanism and the displacement efficiency in the pore doublet model. To this end, published notions are critically reviewed and the weaknesses of the assumptions used in deriving various criteria for the motion of the fluid-fluid interfaces are identified. Theoretical considerations of the motion of fluid-fluid interfaces of different curvatures, interacting in pore doublet arrangements, led to the realization of new mechanisms of entrapment of the displaced phase. These new mechanisms, along with the pore structure parameters and the displacement conditions prevailing at the pore level thatmore » jointly account for entrapment, were established with the help of appropriate displacement experiments recorded on photographs depicting the advance of fluid-fluid interfaces in transparent pore doublet models.« less
  • A mathematical model developed for SiO{sub 2} deposition in porous Vycor glass using SiCl{sub 4} hydrolysis describes reaction, diffusion and evolution of the pore structure due to accumulation of the solid product. The deposition reaction is described by transient heterogeneous kinetics in terms of the concentrations of silanol and chloride groups in the product layer as well as those of the gaseous reactants. For typical deposition conditions the pseudo steady-state approximation for surface species could lead to erroneous predictions. Pore structure evolution is modeled by incorporation results of percolation theory. For this purpose the porous glass is represented by amore » Bethe lattice with coordination number 3 and alternatively by a decorated Bethe lattice in which each bond is replaced by a composite bond consisting of two bonds in series. The second network can capture the effect of pore radius variation along a single pore. For the decorated lattice, pore connectivity interruption at a higher void fraction leads to thinner deposits and shorter deposition times for pore plugging compared to the corresponding ones for the simple lattice.« less