A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage

CO2 injection and storage in deep saline aquifers involves many coupled processes, including multiphase flow, heat and mass transport, rock deformation and mineral precipitation and dissolution. Coupling is especially critical in carbonate aquifers, where minerals will tend to dissolve in response to the dissolution of CO2 into the brine. The resulting neutralization will drive further dissolution of both CO2 and calcite. This suggests that large cavities may be formed and that proper simulation may require full coupling of reactive transport and multiphase flow. We show that solving the latter may suffice whenever two requirements are met: (1) all reactions can be assumed to occur in equilibrium and (2) the chemical system can be calculated as a function of the state variables of the multiphase flow model (i.e., liquid and gas pressure, and temperature). We redefine the components of multiphase flow codes (traditionally, water and CO2), so that they are conservative for all reactions of the chemical system. This requires modifying the traditional constitutive relationships of the multiphase flow codes, but yields the concentrations of all species and all reaction rates by simply performing speciation and mass balance calculations at the end of each time step. We applied this method to the H2O-CO2-Na-Cl-CaCO3 system, so as to model CO2 injection into a carbonate aquifer containing brine. Results were very similar to those obtained with traditional formulations, which implies that full coupling of reactive transport and multi-phase flow is not really needed for this kind of systems, but the resulting simplifications may make it advisable even for cases where the above requirements are not met. Regarding the behavior of carbonate rocks, we find that porosity development near the injection well is small because of the low solubility of calcite. Moreover, dissolution concentrates at the front of the advancing CO2 plume because the brine below the plume tends to reach high CO2 concentrations quite rapidly. We conclude that carbonate dissolution needs not to be feared. © 2013 Elsevier Ltd.

Saved in:
Bibliographic Details
Main Authors: Saaltink, Maarten W., Vilarrasa, Víctor, De Gaspari, F., Silva Rojas, Orlando, Carrera, Jesús, Rötting, Tobias
Format: artículo biblioteca
Published: Elsevier 2013-12
Subjects:Multiphase flow, Density dependent flow, Reactive transport, Numerical models, CO2 sequestration,
Online Access:http://hdl.handle.net/10261/93001
Tags: Add Tag
No Tags, Be the first to tag this record!
id dig-idaea-es-10261-93001
record_format koha
spelling dig-idaea-es-10261-930012022-03-30T11:01:15Z A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage Saaltink, Maarten W. Vilarrasa, Víctor De Gaspari, F. Silva Rojas, Orlando Carrera, Jesús Rötting, Tobias Multiphase flow Density dependent flow Reactive transport Numerical models CO2 sequestration CO2 injection and storage in deep saline aquifers involves many coupled processes, including multiphase flow, heat and mass transport, rock deformation and mineral precipitation and dissolution. Coupling is especially critical in carbonate aquifers, where minerals will tend to dissolve in response to the dissolution of CO2 into the brine. The resulting neutralization will drive further dissolution of both CO2 and calcite. This suggests that large cavities may be formed and that proper simulation may require full coupling of reactive transport and multiphase flow. We show that solving the latter may suffice whenever two requirements are met: (1) all reactions can be assumed to occur in equilibrium and (2) the chemical system can be calculated as a function of the state variables of the multiphase flow model (i.e., liquid and gas pressure, and temperature). We redefine the components of multiphase flow codes (traditionally, water and CO2), so that they are conservative for all reactions of the chemical system. This requires modifying the traditional constitutive relationships of the multiphase flow codes, but yields the concentrations of all species and all reaction rates by simply performing speciation and mass balance calculations at the end of each time step. We applied this method to the H2O-CO2-Na-Cl-CaCO3 system, so as to model CO2 injection into a carbonate aquifer containing brine. Results were very similar to those obtained with traditional formulations, which implies that full coupling of reactive transport and multi-phase flow is not really needed for this kind of systems, but the resulting simplifications may make it advisable even for cases where the above requirements are not met. Regarding the behavior of carbonate rocks, we find that porosity development near the injection well is small because of the low solubility of calcite. Moreover, dissolution concentrates at the front of the advancing CO2 plume because the brine below the plume tends to reach high CO2 concentrations quite rapidly. We conclude that carbonate dissolution needs not to be feared. © 2013 Elsevier Ltd. This work was supported by the European Commission through the MUSTANG and the PANACEA projects (Seventh Framework Programme FP7/2007-2013 under Grant agreements nos. 227286 and 282900, respectively). We also want to acknowledge the financial support received from Instituto para la Diversificación y Ahorro de la Energía (IDAE, Spanish Government) and by the European Union through the “European Energy Programme for Recovery” and the Compostilla OXYCFB300 project. Peer Reviewed 2014-03-05T10:57:58Z 2014-03-05T10:57:58Z 2013-12 2014-03-05T10:57:58Z artículo http://purl.org/coar/resource_type/c_6501 doi: 10.1016/j.advwatres.2013.09.013 issn: 0309-1708 Advances in Water Resources 62(C): 431-441 (2013) http://hdl.handle.net/10261/93001 10.1016/j.advwatres.2013.09.013 #PLACEHOLDER_PARENT_METADATA_VALUE# #PLACEHOLDER_PARENT_METADATA_VALUE# info:eu-repo/grantAgreement/EC/FP7/227286 info:eu-repo/grantAgreement/EC/FP7/282900 none Elsevier
institution IDAEA ES
collection DSpace
country España
countrycode ES
component Bibliográfico
access En linea
databasecode dig-idaea-es
tag biblioteca
region Europa del Sur
libraryname Biblioteca del IDAEA España
topic Multiphase flow
Density dependent flow
Reactive transport
Numerical models
CO2 sequestration
Multiphase flow
Density dependent flow
Reactive transport
Numerical models
CO2 sequestration
spellingShingle Multiphase flow
Density dependent flow
Reactive transport
Numerical models
CO2 sequestration
Multiphase flow
Density dependent flow
Reactive transport
Numerical models
CO2 sequestration
Saaltink, Maarten W.
Vilarrasa, Víctor
De Gaspari, F.
Silva Rojas, Orlando
Carrera, Jesús
Rötting, Tobias
A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
description CO2 injection and storage in deep saline aquifers involves many coupled processes, including multiphase flow, heat and mass transport, rock deformation and mineral precipitation and dissolution. Coupling is especially critical in carbonate aquifers, where minerals will tend to dissolve in response to the dissolution of CO2 into the brine. The resulting neutralization will drive further dissolution of both CO2 and calcite. This suggests that large cavities may be formed and that proper simulation may require full coupling of reactive transport and multiphase flow. We show that solving the latter may suffice whenever two requirements are met: (1) all reactions can be assumed to occur in equilibrium and (2) the chemical system can be calculated as a function of the state variables of the multiphase flow model (i.e., liquid and gas pressure, and temperature). We redefine the components of multiphase flow codes (traditionally, water and CO2), so that they are conservative for all reactions of the chemical system. This requires modifying the traditional constitutive relationships of the multiphase flow codes, but yields the concentrations of all species and all reaction rates by simply performing speciation and mass balance calculations at the end of each time step. We applied this method to the H2O-CO2-Na-Cl-CaCO3 system, so as to model CO2 injection into a carbonate aquifer containing brine. Results were very similar to those obtained with traditional formulations, which implies that full coupling of reactive transport and multi-phase flow is not really needed for this kind of systems, but the resulting simplifications may make it advisable even for cases where the above requirements are not met. Regarding the behavior of carbonate rocks, we find that porosity development near the injection well is small because of the low solubility of calcite. Moreover, dissolution concentrates at the front of the advancing CO2 plume because the brine below the plume tends to reach high CO2 concentrations quite rapidly. We conclude that carbonate dissolution needs not to be feared. © 2013 Elsevier Ltd.
format artículo
topic_facet Multiphase flow
Density dependent flow
Reactive transport
Numerical models
CO2 sequestration
author Saaltink, Maarten W.
Vilarrasa, Víctor
De Gaspari, F.
Silva Rojas, Orlando
Carrera, Jesús
Rötting, Tobias
author_facet Saaltink, Maarten W.
Vilarrasa, Víctor
De Gaspari, F.
Silva Rojas, Orlando
Carrera, Jesús
Rötting, Tobias
author_sort Saaltink, Maarten W.
title A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
title_short A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
title_full A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
title_fullStr A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
title_full_unstemmed A method for incorporating equilibrium chemical reactions into multiphase flow models for CO2 storage
title_sort method for incorporating equilibrium chemical reactions into multiphase flow models for co2 storage
publisher Elsevier
publishDate 2013-12
url http://hdl.handle.net/10261/93001
work_keys_str_mv AT saaltinkmaartenw amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT vilarrasavictor amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT degasparif amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT silvarojasorlando amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT carrerajesus amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT rottingtobias amethodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT saaltinkmaartenw methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT vilarrasavictor methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT degasparif methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT silvarojasorlando methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT carrerajesus methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
AT rottingtobias methodforincorporatingequilibriumchemicalreactionsintomultiphaseflowmodelsforco2storage
_version_ 1777669254393364480