Modeling Diffusion and Buoyancy-Driven Convection with Application to Geological CO2 Storage
ProgramEarth Sciences and Engineering
KAUST DepartmentPhysical Sciences and Engineering (PSE) Division
MetadataShow full item record
AbstractABSTRACT Modeling Diffusion and Buoyancy-Driven Convection with Application to Geological CO2 Storage Rebecca Allen Geological CO2 storage is an engineering feat that has been undertaken around the world for more than two decades, thus accurate modeling of flow and transport behavior is of practical importance. Diffusive and convective transport are relevant processes for buoyancy-driven convection of CO2 into underlying fluid, a scenario that has received the attention of numerous modeling studies. While most studies focus on Darcy-scale modeling of this scenario, relatively little work exists at the pore-scale. In this work, properties evaluated at the pore-scale are used to investigate the transport behavior modeled at the Darcy-scale. We compute permeability and two different forms of tortuosity, namely hydraulic and diffusive. By generating various pore ge- ometries, we find hydraulic and diffusive tortuosity can be quantitatively different in the same pore geometry by up to a factor of ten. As such, we emphasize that these tortuosities should not be used interchangeably. We find pore geometries that are characterized by anisotropic permeability can also exhibit anisotropic diffusive tortuosity. This finding has important implications for buoyancy-driven convection modeling; when representing the geological formation with an anisotropic permeabil- ity, it is more realistic to also account for an anisotropic diffusivity. By implementing a non-dimensional model that includes both a vertically and horizontally orientated 5 Rayleigh number, we interpret our findings according to the combined effect of the anisotropy from permeability and diffusive tortuosity. In particular, we observe the Rayleigh ratio may either dampen or enhance the diffusing front, and our simulation data is used to express the time of convective onset as a function of the Rayleigh ratio. Also, we implement a lattice Boltzmann model for thermal convective flows, which we treat as an analog for CO2 storage modeling. Our model contains the multiple- relaxation-time scheme and moment-based boundary conditions to avoid the numer- ical slip error that is associated with standard bounce-back. The model’s accuracy and robustness is demonstrated by an excellent agreement between our results and benchmark data for thermal flows ranging from Ra = 103 to 108. Our thermal model captures analogous flow behavior to that of CO2 through fluid-filled porous media, including the transition from diffusive transport to initiation and development of convective fingering.