Vertically integrated dual-continuum models for CO2 injection in fractured geological formations

• Published in: Computational geosciences 2019-04, Vol.23 (2), p.273-284
• Main Authors: Tao, Yiheng, Guo, Bo, Bandilla, Karl W., Celia, Michael A.
• Format: Article
• Language: English
• Subjects: Aquifers; Brines; Carbon dioxide; Carbon sequestration; Computational efficiency; Computer Science; Computer simulation; Connecting; Continuum modeling; Dimensions; Earth and Environmental Science; Earth Sciences; Flow equations; Fractures; Geology; Geotechnical Engineering & Applied Earth Sciences; Hydrogeology; Integration; Mass transfer; Mathematical Modeling and Industrial Mathematics; Mathematical models; Migration; Modelling; Original Paper; Permeability; Porosity; Porous media; Rocks; Saline water; Segregation; Soil Science & Conservation; Three dimensional models; Transfer functions
• ISSN: 1420-0597; 1573-1499
• DOI: 10.1007/s10596-018-9805-x

Various modeling approaches, including fully three-dimensional (3D) models and vertical-equilibrium (VE) models, have been used to study the large-scale storage of carbon dioxide (CO 2 ) in deep saline aquifers. 3D models solve the governing flow equations in three spatial dimensions to simulate migration of CO 2 and brine in the geological formation. VE models assume rapid and complete buoyant segregation of the two fluid phases, resulting in vertical pressure equilibrium and allowing closed-form integration of the governing equations in the vertical dimension. This reduction in dimensionality makes VE models computationally much more efficient, but the associated assumptions restrict the applicability of VE models to geological formations with moderate to high permeability. In the present work, we extend the VE models to simulate CO 2 storage in fractured deep saline aquifers in the context of dual-continuum modeling, where fractures and rock matrix are treated as porous media continua with different permeability and porosity. The high permeability of fractures makes the VE model appropriate for the fracture domain, thereby leading to a VE dual-continuum model for the dual continua. The transfer of fluid mass between fractures and rock matrix is represented by a mass transfer function connecting the two continua, with a modified transfer function for the VE model based on vertical integration. Comparison of the new model with a 3D dual-continuum model shows that the new model provides comparable numerical results while being much more computationally efficient.