Multiphase Reactive Transport in a Heterogeneous Flow Field: Channel Formation in Hydrocarbon-Bearing Carbonate Rock
Abstract
Carbon capture, utilization and storage (CCUS) is a key approach for reducing anthropogenic CO2 emissions. Currently the vast majority of CO2 stored is injected into depleted hydrocarbon reservoirs. Reactive transport in hydrocarbon-bearing carbonate reservoirs is controlled not only by the balance between reactant delivery and surface reaction, but also by the heterogeneous pore-scale flow field created by pore structure and fluid distribution. Here, we used time-resolved micro-CT imaging, pore-network analysis, and direct numerical simulation to investigate channelized dissolution during injection of CO2-saturated brine into oil-bearing Ketton limestone at 0.5 mL/min. The system remained in a high Pe, low Da regime throughout the 180 min injection. Despite the advection-dominated flow regime, dissolution became strongly localized. Preferential flow pathways were already present originally because the heterogeneous pore structure and remaining-oil occupancy restricted brine flow to a subset of the connected pore space. Dissolution progressively amplified these pathways. The associated reduction in hydraulic resistance further focused flow, producing a continuous tortuous channel by 180 min. Meanwhile, the effective reaction rate was only 1.6*10^5 mol/m2s,, approximately one order of magnitude lower than the batch reaction rate, showing that rapid advective transport did not translate directly into rapid overall dissolution under multiphase conditions. These results demonstrate that channel formation arose from the coupling between strong reactant delivery and pore-scale flow-field heterogeneity, rather than from the bulk Pe Da regime alone. Accounting for flow-field heterogeneity is therefore important for predicting reactive transport and pore-structure evolution during CO2 injection into hydrocarbon-bearing carbonate reservoirs.
Related
- Solute dispersion in magnetically influenced multiphase flow through a porous tube: axial transport and microrotational effects
- Probing Structure and Ionic Transport in Molten Lithium Carbonate
- Phase-field analysis of fracture in heterogeneous wellbore systems: effects of casing eccentricity and cement-formation interface strength
- Optimal Scaling of Unitary Design Formation in U(1)-Symmetric Random Circuits: A Bottleneck Slower than Charge Transport
- Biomimetic apatite coatings—Carbonate substitution and preferred growth orientation
BibTeX
@misc{ma2026multiphase,
title = {{Multiphase Reactive Transport in a Heterogeneous Flow Field: Channel Formation in Hydrocarbon-Bearing Carbonate Rock}},
author = {Qianqian Ma and Zhuangzhuang Ma and Rukuan Chai and Yanghua Wang and Martin J. Blunt and Branko Bijeljic},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.01328},
eprint = {2609.01328},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.01328},
}
RIS
TY - UNPB
TI - Multiphase Reactive Transport in a Heterogeneous Flow Field: Channel Formation in Hydrocarbon-Bearing Carbonate Rock
AU - Qianqian Ma
AU - Zhuangzhuang Ma
AU - Rukuan Chai
AU - Yanghua Wang
AU - Martin J. Blunt
AU - Branko Bijeljic
PY - 2026
DA - 2026/09/01
PB - arXiv
DO - 10.48550/arxiv.2609.01328
UR - https://arxiv.org/abs/2609.01328
AB - Carbon capture, utilization and storage (CCUS) is a key approach for reducing anthropogenic CO2 emissions. Currently the vast majority of CO2 stored is injected into depleted hydrocarbon reservoirs. Reactive transport in hydrocarbon-bearing carbonate reservoirs is controlled not only by the balance between reactant delivery and surface reaction, but also by the heterogeneous pore-scale flow field created by pore structure and fluid distribution. Here, we used time-resolved micro-CT imaging, pore-network analysis, and direct numerical simulation to investigate channelized dissolution during injection of CO2-saturated brine into oil-bearing Ketton limestone at 0.5 mL/min. The system remained in a high Pe, low Da regime throughout the 180 min injection. Despite the advection-dominated flow regime, dissolution became strongly localized. Preferential flow pathways were already present originally because the heterogeneous pore structure and remaining-oil occupancy restricted brine flow to a subset of the connected pore space. Dissolution progressively amplified these pathways. The associated reduction in hydraulic resistance further focused flow, producing a continuous tortuous channel by 180 min. Meanwhile, the effective reaction rate was only 1.6*10^5 mol/m2s,, approximately one order of magnitude lower than the batch reaction rate, showing that rapid advective transport did not translate directly into rapid overall dissolution under multiphase conditions. These results demonstrate that channel formation arose from the coupling between strong reactant delivery and pore-scale flow-field heterogeneity, rather than from the bulk Pe Da regime alone. Accounting for flow-field heterogeneity is therefore important for predicting reactive transport and pore-structure evolution during CO2 injection into hydrocarbon-bearing carbonate reservoirs.
ER -