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.
Main Authors: | , , , , , |
---|---|
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 |