High-resolution simulation and characterization of density-driven flow in CO2 storage in saline aquifers

Simulations are routinely used to study the process of carbon dioxide (CO2) sequestration in saline aquifers. In this paper, we describe the modeling and simulation of the dissolution–diffusion–convection process based on a total velocity splitting formulation for a variable-density incompressible single-phase model. A second-order accurate sequential algorithm, implemented within a block-structured adaptive mesh refinement (AMR) framework, is used to perform high-resolution studies of the process. We study both the short-term and long-term behaviors of the process. It is found that the onset time of convection follows closely the prediction of linear stability analysis. In addition, the CO2 flux at the top boundary, which gives the rate at which CO2 gas dissolves into a negatively buoyant aqueous phase, will reach a stabilized state at the space and time scales we are interested in. This flux is found to be proportional to permeability, and independent of porosity and effective diffusivity, indicative of a convection-dominated flow. A 3D simulation further shows that the added degrees of freedom shorten the onset time and increase the magnitude of the stabilized CO2 flux by about 25%. Finally, our results are found to be comparable to results obtained from TOUGH2-MP.