Multi-scale Investigations of Geological Carbon Sequestration in Deep Saline Aquifers

Geological carbon dioxide (CO2) sequestration (GCS) in deep saline aquifers is viewed as a viable solution to dealing with the impact of anthropogenic CO2 emissions on global warming. The trapping mechanisms that control GCS include capillary trapping, structural trapping, dissolution trapping, and mineral trapping. Wettability and density-driven convection play an important role in GCS, because wettability significantly affects the efficiency of capillary trapping, and density-driven convection greatly decreases the time scale of dissolution trapping. This work focuses on the role of wettability on multiphase flow in porous media, density-driven convection in porous media, and their implications for GCS in deep saline aquifers. Wettability is a critical control over multiphase fluid flow in porous media. However, our understanding on the wettability heterogeneity of a natural rock and its effect on multiphase fluid flow in a natural rock is limited. This work innovatively models the heterogeneous wettability of a rock as a correlated random field. The realistic wetting condition of a natural rock can be reconstructed with in-situ measurements of wettability on the internal surfaces of the rock. A Bentheimer sandstone was used to demonstrate the workflow to model and reconstruct a wettability field. Relative permeability, capillary pressure-water saturation relation are important continuum-scale properties controlling multiphase flow in porous media. This work employed lattice Boltzmann method to simulate the displacement process. We found that pore-scale surface wettability heterogeneity caused noticeable local scCO2 and water redistributions under less water-wet conditions at the pore scale. At the continuum scale, the capillary pressure-water saturation curve under the heterogeneous wetting condition was overall similar to that under the homogeneous wetting condition. This suggested that the impact of local wettability heterogeneity on the capillary pressure-water saturation curve was averaged out at the entire-sample scale. The only difference was that heterogeneous wettability led to a negative entry pressure at the primary drainage stage under the intermediate-wet condition. The impact of pore-scale wettability heterogeneity was more noticeable on the relative permeability curves. Particularly, the variation of the scCO2 relative permeability curve in the heterogeneous wettability scenario was more significant than that in the homogenous wettability scenario. Results showed that higher wettability heterogeneity (i.e., higher standard deviation and higher correlation length) increased the variations in the CO2/brine relative permeability curves. Dissolution of CO2 into brine is a primary mechanism to ensure the long-term security of GCS. CO2 dissolved in brine increases the CO2-brine solution density and thus can cause downward convection. Onset of density-driven instability and onset of convective dissolution are two critical events in the transition process from a diffusion-dominated regime to a convection-dominated regime. In the laboratory, we developed an empirical correlation between light intensity and in-situ solute concentration. Based on the novel and well-controlled experimental methods, we measured the critical Rayleigh-Darcy number and critical times for the onset of density-driven instability and convective dissolution. To further investigate the impact of permeability heterogeneity on density-driven convection, a three-dimensional (3D) fluidics method was proposed to advance the investigation on density-driven convection in porous media. Heterogeneous porous media with desired spatial correlations were efficiently built with 3D-printed elementary porous blocks. In the experiments, methanol-ethylene-glycol (MEG), was used as surrogate fluid to CO2. The heterogeneous porous media were placed in a transparent tank allowing visual observations. Results showed that permeability structure controlled the migration of MEG-rich water. Permeability heterogeneity caused noticeable uncertainty in dissolution rates and uncertainty in dissolution rates increases with correlation length. To sum up, this work comprehensively employed novel experimental methods and large-scale direct simulations to investigate the sequestration of CO2 in saline aquifers at a pore scale and a continuum scale. The findings advanced our understanding on the role of wettability heterogeneity and permeability heterogeneity on GCS in deep saline aquifers.