Characterizing small-scale migration behavior of sequestered CO2 in a realistic geological fabric

For typical reservoir conditions, buoyancy and capillary forces grow dominant over viscous forces within a few hundred meters of the injection wells as the pressure gradient due to injection decreases, resulting in qualitatively different plume migration regimes. The migration regime depends on two factors: the capillary pressure of the leading edge of the plume and the range of threshold entry pressures within the rock at the leading edge of the plume. A capillary channel regime arises when these two factors have the same magnitude. Flow patterns within this regime vary from finger-like structures with minimal rock contact to back-filling structures with compact volumes of saturation distributed between fingers. Reservoir heterogeneity is one of the principal factors influencing CO2 migration pathway in the capillary channel regime. Here we characterize buoyancy-driven migration in a natural 2D geologic domain (1 m × 0.5 m peel from an alluvium) in which sedimentologic heterogeneity has been resolved at sub-millimeter (depositional) resolution. The relevant features of the heterogeneity are grain size distribution, which determines the mean and range of threshold pressures and correlation lengths of threshold pressures in horizontal and vertical directions. The relevant physics for this migration regime is invasion percolation, and simulations indicate that CO2 migrates through the peel in a few narrow pathways which cannot be captured by conventional coarse-grid simulations. The storage efficiency of the capillary channel regime would be low and consequently CO2 would also migrate greater distances than expected from models or simulations that neglect the capillary channel flow regime.