Residual trapping of carbon dioxide during geological storage-Insight gained through a pore-network modeling approach

To reduce emissions of the greenhouse gas CO2 to the atmosphere, sequestration in deep saline aquifers is a viable strategy. Residual trapping is a key containment process important to the success of CO2 storage operations. While residual trapping affects CO2 migration over large scales, it is inherently a pore-scale process. Pore-network models (PNMs), capturing such processes, are useful for our understanding of residual trapping, and for upscaling trapping parameters for larger scale models. A PNM for simulation of quasi-static two-phase flow; CO2 intrusion (drainage) followed by water flooding (imbibition) was developed. It accounts for pore-scale displacement mechanisms, and was used to investigate residual CO2 trapping. The sensitivity of the residual CO2 saturation to several parameters was studied, to validate a trapping behavior in agreement with earlier studies. Then the PNM was calibrated to core sample data and used to simulate drainage-imbibition scenarios with different turning point saturations. From these the initial-residual saturation curves of CO2 in Heletz sandstone were estimated, essential for future macroscopic-scale simulations. Further, the occurrence of different pore-scale mechanisms were quantified and the size distribution of the residual clusters was shown to exhibit a bimodal appearance. The findings improve the understanding of residual trapping in Heletz sandstone.