Elsevier

Geothermics

Volume 89, January 2021, 101947
Geothermics

Geothermal fluid circulation in a caldera setting: The Torre Alfina medium enthalpy system (Italy)

https://doi.org/10.1016/j.geothermics.2020.101947Get rights and content

Highlights

  • New conceptual model of Torre Alfina geothermal field based on THOUGH2 simulations.

  • Reservoir geometry and structural discontinuities control heat and fluid flow patterns.

  • Modelling results discount the presence of local magmatic source and meteoric recharge.

  • We propose Torre Alfina is a “blind” geothermal system recharged by high enthalpy fluids from Bolsena caldera.

Abstract

The Torre Alfina medium enthalpy geothermal field is located about 10 km north of the Bolsena caldera (Italy). The reservoir is a buried structural high consisting of fractured Meso-Cenozoic carbonate sequences and sealed by clayey flysch successions and Pliocene marine clays. We performed TOUGH2 numerical simulations, testing different model designs based on contrasting conceptual models.

Results indicate that deep circulation is forced by the geometry of the reservoir and by the applied T and P gradients. We interpret the Torre Alfina field as a "blind" system, mostly recharged by lateral advection of heat and fluids from the Bolsena caldera deep high-enthalpy system, through the permeable caldera faults.

Introduction

With an estimated electrical potential of 10–100 times the current production, geothermal energy currently contributes to less than 1% of the world’s primary energy supply (World Energy Resources, Geothermal/World Energy Council, 2016), and has the potential to grow up to 8.3 % of the global electricity production and serve 17 % of the world’s population by 2050 (Bertani, 2016).

Although the heat supply is abundant, the distribution of the geothermal resources is not homogeneous and is localized mainly in active geodynamic regions such as volcanic or intrusive systems and active tectonic zones (Henley and Ellis, 1983; Moeck, 2014).

The majority of the known geothermal systems was localized thanks to the presence of surface manifestations. However, many resources show little or no surface expression and the identification of a new geothermal reservoir can be challenging (Hanson et al., 2014). The amount of “blind” geothermal systems - whose presence gives no signature at the surface (Forson et al., 2014) – remains uncertain. Considering the potential of these undiscovered resources and the high costs of drilling programs, the search for blind systems in diverse geologic environments has been carried out with cost-effective methods including geochemical, geophysical and spatial analysis/structural geology surveys (Bibby et al., 1994; Blackwell and Baag, 1973; Deon et al., 2015; Faulds et al., 2015, 2016; Fowler et al., 2017; Garg et al., 2010; Hanson et al., 2014; Kratt et al., 2009).

Geothermal systems in volcanic regimes are usually characterized by the highest enthalpy and occur in different geological settings such as plate boundaries, continental rifts and hot spot environments (Saemundsson et al., 2009).

Among the volcanic structures, calderas host some of the most consistent geothermal reservoirs. Due to the slow cooling of the magma bodies, the deep heat source can be available for up to 2 million years after the last volcanic activity. Its interaction with groundwater can generate hydrothermal systems developing in permeable rock volumes (e.g. Giordano et al., 2014; Wohletz and Heiken, 1992).

However, despite the high geothermal potential, investigation and exploitation in the interior of these structures can encounter environmental and logistical complications, which can represent an obstacle to the geothermal development itself. Volcanic regions, in fact, are often of historic/environmental significance and/or protected tourist destinations (Newsome and Dowling, 2018). Furthermore, these areas are often sites of intense urbanization with extremely high population density: Campi Flegrei, Italy, hosts more than 550.000 inhabitants within the area of the caldera, and the Kagoshima Bay in Japan has a population of more than 850.000 people. In such densely urbanized settings, any industrial development requires a careful hazard assessment which, in case of seismic areas (McNutt and Roman, 2015), should also include the effects of possible induced seismicity (Carlino et al., 2016; Gaucher et al., 2015; Lermo et al., 2008).

Besides the environmental and geologic issues, there are also social challenges, such as the skepticism of local communities and administrations (Manzella et al., 2018), and technical problems related to the high chemical aggressiveness of hydrothermal fluids on the materials of the plant (Nogara and Zarrouk, 2018).The presence of lakes (e.g. Colli Albani and Monti Vulsini, Italy; Mt Pinatubo, Philippines; Crater Lake, Oregon; Askja, Iceland) and seas submerging the caldera structures in part (e.g. Campi Flegrei, Italy; Santorini, Greece; Aira, Japan) or entirely (Kulo Lasi, Tonga; Kaikata, Japan) (Branney and Acocella, 2015; Wohletz and Heiken, 1992) can also create problems.

Considering these complications, the study of geothermal fluid migrations away from the main volcanic structures can be a successful approach for the identification of new strategic resources. As already mentioned, calderas are areas of intense hydrothermal circulation which is not only influenced by the different permeability of the volcanic deposits (e.g. Heap et al., 2014), but also by the complex geometry of structural discontinuities such as collapse ring faults, fractures, caldera-crossing faults and regional faults which may act as preferential channels for fluid migration (Garden et al., 2017; Giordano et al., 2013, 2014; Urbani et al., 2020; Vignaroli et al., 2013) and promote lateral advection of hot fluids away from the central part of the volcano.

Studies performed in non-caldera settings showed that faults can act both as high permeability conducts or as hydraulic barriers, depending on mineralization or the comminution of fault core and fracturation of damage zone (e.g. Faulkner and Rutter, 2001; Giordano et al., 2013; Rowland and Sibson, 2004; Vignaroli et al., 2013). However, only a few studies deal with the permeability of caldera collapse faults (Carlino, 2018; Garden et al., 2017; Jasim et al., 2015), because exposed outcrops of active structures are rare (Garden et al., 2017; Geyer and Martí, 2014).

The connection between caldera structural discontinuities and the anisotropic permeability and consequent fluids channelization in the internal margin of the volcano is now confirmed by drilling data (Carlino, 2018), numerical modelling (Jasim et al., 2015) and field structural surveys (Garden et al., 2017). The role of these structures in influencing hydrothermal fluids circuits towards peripheral areas outside the caldera margin remains unexplored.

The main aim of this study is to investigate the hydrothermal circulation where a caldera collapse structure intersects a regional reservoir, to define the role of caldera faults in the deep hydrothermal circuits and to assess their contribution in promoting lateral advection to peripheral areas.

To this purpose, numerical modelling was applied to the Torre Alfina geothermal system, for which a rich database is available from previous geologic studies and whose characteristics make it a suitable case study for modelling the paths of fluid migrations.

Results contribute to the comprehension of the mechanisms of fluid circulation and heat supply from the caldera to peripheral geothermal systems, with relevant implications for geothermal exploration planning and for targeting new distal resources.

Section snippets

Geological framework

The Torre Alfina geothermal area is located in Central Italy, along the Tyrrhenian margin of the Apennine Mts (Fig. 1a), and at the northernmost tip of the Quaternary alkali potassic volcanoes of the Vulsini Volcanic District. The evolution of the Apennine chain was characterized, at first, by Miocene-Pliocene compressional tectonic phases, with the eastward migration of thrust fronts and foredeep basins in a classical piggyback sequence; then, post-orogenic extension followed and is still

Previous conceptual model of Torre Alfina geothermal system

The geothermal reservoir of Torre Alfina was discovered after some preliminary regional studies aimed at identifying potential geothermal resources. Since 1973, the area has been investigated by geological, geophysical, geochemical and hydrogeological surveys (e.g. Barberi et al., 1994; Buonasorte et al., 1988, 1991; Chiarabba et al., 1995; Chiodini et al., 1995; Costantini et al., 1984; Doveri et al., 2010; Duchi et al., 1992).

The drilling of deep boreholes (Fig. 1) allowed the definition of

Methodology

We tested different conceptual models of the Torre Alfina geothermal system by applying the TOUGH2 geothermal simulator (Pruess et al., 1999). This numerical simulation program, based on an integral finite difference method, describes the coupled flow of heat and fluid through heterogeneous porous media. The fluid propagation is described according to a multi-phase version of Darcy’s law that accounts for water phase changes with latent heat effects. Phase interference is taken into account

Discussion

The first outcome of the simulations is the influence of the geometry of the permeable domain in hydrothermal fluid circulations: in SIM1 Case A: Closed system (Fig. 3), despite the absence of external perturbations, convective cells of variable extensions form and the asymmetry of the domain produces asymmetric circuits. This highlights the role of the structural setting of the geothermal reservoir as well as the contrasting permeabilities of the different rock units; the geometry of these

Conclusions

This study allowed us to understand some aspects of Torre Alfina geothermal system with implications also for processes that can take place in geothermal reservoirs.

The main considerations about geothermal system dynamics are:

  • -

    the presence of structural highs of reservoir rocks can promote the development of thermal anomalies, even in areas with no heat sources, because they act as a “trap” of deeper and warmer groundwater in shallower areas of the crust;

  • -

    the domain geometry and the distribution

Author statement

G.G. designed the research; I.T. and M.T. performed the numerical modelling; all the authors discussed the data and wrote the paper. All together revised and approved the final form of the manuscript.

Declaration of Competing Interest

The authors report no declarations of interest.

References (74)

  • S. Conticelli

    The effect of crustal contamination on ultrapotassic magmas with lamproitic affinity: mineralogical, geochemical and isotope data from the Torre Alfina lavas and xenoliths, Central Italy

    Chem. Geol.

    (1998)
  • S. Conticelli et al.

    Petrology of final-stage Latera lavas (Vulsini Mts.): mineralogical, geochemical and Sr-isotopic data and their bearing on the genesis of some potassic magmas in central Italy

    J. Volcanol. Geotherm. Res.

    (1991)
  • S. Conticelli et al.

    Shift from lamproite-like to leucititic rocks: Sr–Nd–Pb isotope data from the Monte Cimino volcanic complex vs. the Vico stratovolcano, Central Italy

    Chem. Geol.

    (2013)
  • M. Doveri et al.

    Revision, calibration, and application of the volume method to evaluate the geothermal potential of some recent volcanic areas of Latium, Italy

    Geothermics

    (2010)
  • V. Duchi et al.

    Chemical composition of thermal springs, cold springs, streams, and gas vents in the Mt. Amiata geothermal region (Tuscany, Italy)

    J. Volcanol. Geotherm. Res.

    (1987)
  • V. Duchi et al.

    Chemical relationship between discharging fluids in the Siena-Radicofani graben and the deep fluids produced by the geothermal fields of Mt Amiata, Torre Alfina and Latera (Central Italy)

    Geothermics

    (1992)
  • A.P. Fowler et al.

    Identification of blind geothermal resources in surprise Valley, CA, using publicly available groundwater well water quality data

    Appl. Geochem.

    (2017)
  • E. Gaucher et al.

    Induced seismicity in geothermal reservoirs: a review of forecasting approaches

    Renew. Sustain. Energy Rev.

    (2015)
  • G. Giordano et al.

    Structural control on geothermal circulation in the Cerro Tuzgle–Tocomar geothermal volcanic area (Puna plateau, Argentina)

    J. Volcanol. Geotherm. Res.

    (2013)
  • G. Giordano et al.

    Incorporating surface indicators of reservoir permeability into reservoir volume calculations: application to the Colli Albani caldera and the Central Italy Geothermal Province

    Earth Sci. Rev.

    (2014)
  • M.C. Hanson et al.

    Identifying blind geothermal systems with soil CO2 surveys

    Appl. Geochem.

    (2014)
  • R.W. Henley et al.

    Geothermal systems ancient and modern: a geochemical review

    Earth Sci. Rev.

    (1983)
  • A. Jasim et al.

    Impact of channelized flow on temperature distribution and fluid flow in restless calderas: insight from Campi Flegrei caldera, Italy

    J. Volcanol. Geotherm. Res.

    (2015)
  • A. Manzella et al.

    Environmental and social aspects of geothermal energy in Italy

    Geothermics

    (2018)
  • F. Marra et al.

    40Ar/39Ar single crystal dates from pyroclastic deposits provide a detailed record of the 590–240 ka eruptive period at the Vulsini Volcanic District (central Italy)

    J. Volcanol. Geotherm. Res.

    (2020)
  • S.R. McNutt et al.

    Volcanic seismicity

    The Encyclopedia of Volcanoes

    (2015)
  • I.S. Moeck

    Catalog of geothermal play types based on geologic controls

    Renew. Sustain. Energy Rev.

    (2014)
  • D. Newsome et al.

    Geoheritage and Geotourism

    Geoheritage

    (2018)
  • J. Nogara et al.

    Corrosion in geothermal environment part 2: metals and alloys

    Renew. Sustain.Energy Rev.

    (2018)
  • G. Poli et al.

    Geochemical characteristics of the South Tuscany (Italy) volcanic province: constraints on lava petrogenesis

    Chem. Geol.

    (1984)
  • G. Vignaroli et al.

    Structural compartmentalisation of a geothermal system, the Torre Alfina field (central Italy)

    Tectonophysics

    (2013)
  • V. Acocella et al.

    Caldera structure, amount of collapse, and erupted volumes: the case of Bolsena caldera, Italy

    Geol. Soc. Am. Bull.

    (2012)
  • F. Barberi et al.

    Plio-pleistocene geological evolution of the geothermal area of Tuscany and Latium

    Memorie Descrittive della Carta Geologica d’Italia

    (1994)
  • D.D. Blackwell et al.

    Heat flow in a “blind” geothermal area near Marysville, Montana

    Geophysics

    (1973)
  • C. Boni et al.

    Carta idrogeologica, scala 1: 500.000. Schema idrogeologico dell’Italia centrale

    Mem. Soc. Geol. It

    (1986)
  • G. Buonasorte et al.

    Ricerca ed esplorazione nell’area geotermica di Torre Alfina (Lazio-Umbria)

    Bollettino della Società Geologica Italiana

    (1988)
  • G. Buonasorte et al.

    The Alfina 15 well: deep geological data from northern Latium (Torre Alfina geothermal area)

    Bollettino della Società geologica italiana

    (1991)
  • Cited by (3)

    View full text