An advanced reactive transport simulation scheme for hydrothermal systems modelling. (March 2019)

• Record Type: Journal Article
• Authors: Yapparova, Alina; Miron, George D.; Kulik, Dmitrii A.; Kosakowski, Georg; Driesner, Thomas
• Is Part Of: Geothermics. Volume 78(2019)
• Page Start: 138
• Page End: 153
• Publication Date: 2019-03
• Subjects: Reactive transport modelling -- Gibbs energy minimization -- Multiphase flow -- boiling -- Partitioning of gases -- Hydrothermal systems; Hydrogeology -- Periodicals; Geothermal resources -- Periodicals; Énergie géothermique -- Périodiques; GEOTHERMAL ENGINEERING; GEOTHERMAL ENERGY; GEOTHERMAL EXPLORATION; Geothermal resources; Hydrogeology; Periodicals; Electronic journals; 621.44
• DOI: 10.1016/j.geothermics.2018.12.003
• Languages: English
• ISSNs: 0375-6505

Abstract

Highlights: Using non-ideal solid solution models instead of pure mineral phases in reactive transport simulations has a significant impact on the pH of the aqueous solution and the overall reaction progress. CSMP++GEM coupled code can be used with confidence in modelling natural geothermal systems up to 350 °C and can perform simulations at hydrothermal conditions (e.g., 600 °C, 1 kbar), which is difficult or impossible with other codes. CSMP++GEM coupled code can model two-phase vapour + liquid flow with partitioning of chemical species between the two phases. Preferential partitioning of CO2 to vapour has a major effect on the pH of a boiling solution. A highly acidic zone may develop at the border between the vapour-dominated and the two-phase zones. Abstract: Advanced reactive transport simulations of hydrothermal systems require the modelling of diverse physico-chemical phenomena, many of which have only partially been addressed in previous studies. Utilizing newly implemented capabilities of the CSMP++GEM code, we show the relevance of three phenomena: (1) At the example of alkali feldspars, we investigate how employing non-ideal solid solution models instead of pure mineral phases has a significant impact on the resulting chemical evolution of the system. For example, pH when considering solid solution differs by about 0.5 compared to a simulation with pure feldspar phases. Differences in aqueous K/Na ratios may be relevant for interpretation of the K/Na geothermometer and may, in addition, be biased by discrepancies in published thermodynamic databases as observed in this study. (2) Reactive transport at supercritical conditions (within the limits imposed by the HKF thermodynamic model) is illustrated by how a magmatic-hydrothermal fluid, equilibrated with granitic rock at 600 °C and 200 MPa, is reacted along a diffusion path with a similar rock at 400 °C and 50 MPa. This simple conceptual representation of hydrothermal alteration at the Butte porphyry copper deposit in Montana results in a reactive transport model that closely mimics the observed alteration zonation. (3) Emergent behaviour is observed for a simulation of two-phase vapour + liquid flow with partitioning of chemical species between the two phases. In this example, hot vapour front condenses into and displaces liquid, with a two-phase boiling/condensation zone developing in between. At the boundaries between steam and boiling zone, and between boiling zone and liquid zone, the interplay of partitioning of volatile species and different transport properties in the single- and two-phase regions leads to distinct local and transient effects: volatile species such as CO2, H2 S and HCl may accumulate or deplete at strongly varying rates, affecting pH, redox state and activity of sulfide. The correct simulation of these effects will be relevant in understanding hydrothermal processes such as ore formation, mineral alteration, and the evolution of high enthalpy geothermal areas.