Vertically integrated dual-continuum models for CO2 injection in fractured geological formations
AffiliationUniv Arizona, Dept Hydrol & Atmospher Sci
KeywordsGeologic CO2 storage
Vertically integrated models
MetadataShow full item record
CitationTao, Y., Guo, B., Bandilla, K.W. et al. Comput Geosci (2019) 23: 273. https://doi.org/10.1007/s10596-018-9805-x
Rights© Springer Nature Switzerland AG 2019.
Collection InformationThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at firstname.lastname@example.org.
AbstractVarious 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 (CO2) in deep saline aquifers. 3D models solve the governing flow equations in three spatial dimensions to simulate migration of CO2 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 CO2 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.
Note12 month embargo; published online: 12 January 2019
VersionFinal accepted manuscript
SponsorsCarbon Mitigation Initiative at Princeton University; U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) [DE-FE0023323]; DOE/NETL