Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

Metamorphic and geochronogical study of the Triassic El Oro metamorphic complex, Ecuador: Implications for high-temperature metamorphism in a forearc zone

In the forearc of the Andean active margin in southwest Ecuador, the El Oro metamorphic complex exhibits a well exposed tilted forearc section partially migmatized. We used Raman spectroscopy on carbonaceous matter (RSCM) thermometry and pseudosections coupled with mineralogical and textural studies to constrain the pressure–temperature (P–T) evolution of the El Oro metamorphic complex during Triassic times. Our results show that anatexis of the continental crust occurred by white-mica and biotite dehydration melting along a 10 km thick crustal domain (from 4.5 to 8 kbar) with increasing temperature from 650 to 700 °C. In the biotite dehydration melting zone, temperature was buffered at 750–820 °C in a 5 km thick layer. The estimated average thermal gradient during peak metamorphism is of 30 °C/km within the migmatitic domain can be partitioned into two apparent gradients parts. The upper part from surface to 7 km depth records a 40–45 °C/km gradient. The lower part records a quasi-adiabatic geotherm with a 10 °C/km gradient consistent with an isothermal melting zone. Migmatites U–Th–Pb geochronology yielded zircon and monazite ages of 229.3 ± 2.1 Ma and 224.5 ± 2.3 Ma, respectively. This thermal event generated S-type magmatism (the Marcabeli granitoid) and was immediately followed by underplating of the high-pressure low-temperature (HP-LT) Arenillas–Panupalí unit at 225.8 ± 1.8 Ma. The association of high-temperature low-pressure (HT-LP) migmatites with HP-LT unit constitutes a new example of a paired metamorphic belt along the South American margin. We propose that in addition to crustal thinning, underplating of the Piedras gabbroic unit before 230 Ma provided the heat source necessary to foster crustal anatexis. Furthermore, its MORB signature shows that the asthenosphere was involved as the source of the heat anomaly. S-type felsic magmatism is widespread during this time and suggests that a large-scale thermal anomaly affected a large part of the South American margin during the late Triassic. We propose that crustal anatexis is related to an anomaly that arose during subduction of the Panthalassa ocean under the South American margin. Slab verticalization or slab break-off can be invoked as the origin of the upwelling of the asthenosphere.

Precambrian to Paleozoic zircon record in the Siviez-Mischabel basement (western Swiss Alps)

U–Pb zircon analyses from three meta-igneous and two metasedimentary rocks from the Siviez-Mischabel nappe in the western Swiss Alps are presented, and are used to derive an evolutionary history spanning from Paleoarchean crustal growth to Permian magmatism. The oldest components are preserved in zircons from metasedimentary albitic schists. The oldest zircon core in these schists is 3.4 Ga old. Detrital zircons reveal episodes of crustal growth in the Neoarchean (2.7–2.5 Ga), Paleoproterozoic (2.2–1.9 Ma) and Neoproterozoic (800–550 Ma, Pan-African event). The maximum age of deposition for the metasedimentary rocks is given by the youngest detrital zircons within both metasedimentary samples dated at ~490 Ma (Cambrian-Ordovician boundary). This is in the age range of two granitoid samples dated at 505 ± 4 and 482 ± 7 Ma, and indicates sedimentation and magmatism in an extensional setting preceding an Ordovician orogeny. The third felsic meta-igneous rock gives a Permian age of intrusion, and is part of a long-lasting Variscan to post-Variscan magmatic activity. The zircons record only minor disturbance of the U–Pb system during the Alpine orogeny.

The Basement of the Deccan Traps and Its Madagascar Connection: Constraints from Xenoliths

Paleogeographic reconstructions of India and Madagascar before their late Cretaceous rifting juxtapose the Antongil Block of Madagascar against the Deccan Traps of India, indicating that the Western Dharwar Craton extends below the Deccan lavas. Some recent studies have suggested that the South Maharashtra Shear Zone along the northern Konkan coast of India limits the northern extent of the Western Dharwar Craton, implying that the craton does not extend below the Deccan Traps, raising a question mark on paleogeographic reconstructions of India and Madagascar. The continuity of the Western Dharwar Craton north of the South Maharashtra Shear Zone below the Deccan Traps—or its lack thereof—is critical for validating tectonic models correlating Madagascar with India. In this study, zircons in tonalitic basement xenoliths hosted in Deccan Trap dykes were dated in situ, using the U-Pb isotope system. The data furnish U-Pb ages that define three populations at 2527 ± 6, 2456 ± 6, and 2379 ± 9 Ma. The 2527 ± 6 Ma ages correspond to the igneous crystallization of the tonalites, whereas the 2456 ± 6 and 2379 ± 9 Ma ages date metamorphic overprints. The results help to establish for the first time that the basement is a part of the Neoarchean granitoid suite of the Western Dharwar Craton, which extends northward up to at least Talvade in central and Kihim beach in the western Deccan. By implication, the South Maharashtra Shear Zone cannot be the northern limit of the Western Dharwar Craton. The granitoids are correlated with the Neoarchean felsic intrusions (2.57–2.49) of the Masaola suite in the Antongil Block of Madagascar, supporting the existence of a Neoarchean Greater Dharwar Craton comprising the Western Dharwar Craton and the Antongil-Masora Block.