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.

Long-term fluid circulation in extensional faults in the central Catalan Coastal Ranges: P-T constraints from neoformed chlorite and K-white mica

The neoformation of chlorite and K-white mica in fault rocks from two main faults of the central Catalan Coastal Ranges, the Vallès and the Hospital faults, has allowed us to constrain the P–T conditions during fault evolution using thermodynamic modeling. Crystallization of M1 and M2 muscovite and microcline occured as result of deuteric alteration during the exhumation of the pluton (290 °C > T > 370 °C) in the Permian. After that, three tectonic events have been distinguished. The first tectonic event, attributed to the Mesozoic rifting, is characterized by precipitation of M3 and M4 phengite together with chlorite and calcite C1 at temperatures between 190 and 310 °C. The second tectonic event attributed to the Paleogene compression has only been identified in the Hospital fault with precipitation of low-temperature calcite C2. The shortcut produced during inversion of the Vallès fault was probably the responsible for the lack of neoformed minerals within this fault. Finally, the third tectonic event, which is related to the Neogene extension, is characterized in the Vallès fault by a new generation of chlorite, associated with calcite C4 and laumontite, formed at temperatures between 125 and 190 °C in the absence of K-white mica. Differently, the Hospital fault is characterized by the precipitation of calcite C3 during the syn-rift stage at temperatures around 150 °C and by low-temperature fluids precipitating calcites C5, C6 and PC1 during the post-rift stage. During the two extensional events (Mesozoic and Neogene), faults acted as conduits for hot fluids producing anomalous high geothermal gradients (50 °C/km minimum).

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Complex geomorphologic assemblage of terrains in association with the banded terrain in Hellas basin, Mars

Hellas basin acts as a major sink for the southern highlands of Mars and is likely to have recorded several episodes of sedimentation and erosion. The north-western part of the basin displays a potentially unique Amazonian landscape domain in the deepest part of Hellas, called “banded terrain”, which is a deposit characterized by an alternation of narrow band shapes and inter-bands displaying a sinuous and relatively smooth surface texture suggesting a viscous flow origin. Here we use high-resolution (HiRISE and CTX) images to assess the geomorphological interaction of the banded terrain with the surrounding geomorphologic domains in the NW interior of Hellas to gain a better understanding of the geological evolution of the region as a whole. Our analysis reveals that the banded terrain is associated with six geomorphologic domains: a central plateau named Alpheus Colles, plain deposits (P1 and P2), reticulate (RT1 and RT2) and honeycomb terrains. Based on the analysis of the geomorphology of these domains and their cross-cutting relationships, we show that no widespread deposition post-dates the formation of the banded terrain, which implies that this domain is the youngest and latest deposit of the interior of Hellas. Therefore, the level of geologic activity in the NW Hellas during the Amazonian appears to have been relatively low and restricted to modification of the landscape through mechanical weathering, aeolian and periglacial processes. Thermophysical data and cross-cutting relationships support hypotheses of modification of the honeycomb terrain via vertical rise of diapirs such as ice diapirism, and the formation of the plain deposits through deposition and remobilization of an ice-rich mantle deposit. Finally, the observed gradual transition between honeycomb and banded terrain suggests that the banded terrain may have covered a larger area of the NW interior of Hellas in the past than previously thought. This has implications on the understanding of the evolution of the deepest part of Hellas.

UHP Metamorphism Documented in Ti-chondrodite- and Ti-clinohumite-bearing Serpentinized Ultramafic Rocks from Chinese Southwestern Tianshan

The southwestern Tianshan (China) metamorphic belt records high-pressure (HP) to ultrahigh-pressure (UHP) conditions corresponding to a cold oceanic subduction-zone setting. Serpentinites enclosing retrogressed eclogite and rodingite occur as lenses within metapelites in the UHP unit, which also hosts coesite-bearing eclogites. Based on the petrology and petrography of these serpentinites, five events are recognized: (1) formation of a wehrlite–harzburgite–dunite association in the mantle; (2) retrograde metamorphism and partial hydration during exhumation of the mantle rocks close to the seafloor; (3) oceanic metamorphism leading to the first serpentinization and rodingitization; (4) UHP metamorphism during subduction; (5) retrograde metamorphism during exhumation together with a second serpentinization. The peak metamorphic mineral assemblage of the serpentinized wehrlite comprises Ti-chondrodite + olivine + antigorite + chlorite + magnetite + brucite. A computed pseudosection for this serpentinized wehrlite shows that the Al content in antigorite is mostly sensititive to temperature but can also be used to constrain pressure. The average XAl = 0·204 ± 0·026 of antigorite (XAl = Al (a.p.f.u.)/8, where Al is in atoms per formula unit for a structural formula M48T34O85(OH)62, and M and T are octahedral and tetrahedral sites, respectively) included in Ti-chondrodite and average XAl = 0·203 ± 0·019 of antigorite in the matrix result in a well-constrained peak metamorphic temperature of 510–530°C. Peak pressures are less precisely constrained at 37 ± 7 kbar. The Tianshan serpentinites thus record UHP metamorphic conditions and represent the deepest subducted serpentinites discovered so far. The retrograde evolution occurs within the stability field of brucite + antigorite + olivine + chlorite and formation of Ti-clinohumite at the expense of Ti-chondrodite has been observed, suggesting isothermal decompression. The resulting P–T path is in excellent agreement with the metamorphic evolution of country rocks, indicating that the UHP unit in Tianshan was subducted and exhumed as a coherent block. To refine the metamorphic path of the ultramafic rocks, we have investigated the stability fields of Ti-chondrodite and Ti-clinohumite using piston-cylinder experiments. A total of 11 experiments were conducted at 25–55 kbar and 600–750°C in a F-free natural system. Combined with previous experiments and information from natural rocks we constructed a petrogenetic grid for the stability of Ti-chondrodite and Ti-clinohumite in F-free peridotite compositions. The formation of Ti-chondrodite in serpentinites requires a minimum pressure of about 26 kbar, whereas in Ti-rich systems it can form at considerably lower pressures. A key finding is that at UHP conditions, F-free Ti-chondrodite or Ti-clinohumite breaks down in the presence of orthopyroxene between 700 and 750°C, at temperatures that are significantly lower than those of the terminal breakdown reactions of these humite minerals. These breakdown reactions are an additional source of fluid during prograde subduction of serpentinites.