The geothermal exploration of Campanian volcanoes: Historical review and future development

Since Roman time, the heat produced by Neapolitan volcanoes was an appeal for people living in and outside the area, for the fruition of the famous thermal baths. This very large area, which spans from Campi Flegrei and Ischia calderas to Somma-Vesuvius volcano, is characterized by high temperature at shallow depth and intense heat flow, and is yet utilized for the bathing and spa treatment industry, while only in the middle of the 20th century a tentative of geothermal exploitation for energy production was performed. Pioneering researches of geothermal resource were carried out in Campanian region since 1930, until 1985, during which a large amount of geological data were collected. In this paper, we make for the first time a review of the history of geothermal explorations in the active Campanian volcanic area. By the analysis of a great amount of literature data and technical reports we reconstruct the chronology and the main information of the drillings performed since 1930 by the SAFEN Company and successively in the framework of the ENEL-AGIP Joint Venture for geothermal exploration. The available data are utilized to correlate the temperatures measured within the deeper wells with the possible sources of geothermal heat in the shallow crust, down to about 8-10 km of depth. Finally, we assess the geothermal potential of the hottest areas, Ischia Island and Campi Flegrei, which have shown the best data and favorable physical conditions for a reliable, and cost-effective, exploitation for thermal and electric purposes.     

Geochemical exploration of the three most significant geothermal areas of Lesbos Island, Greece

Lesbos Island has several thermal manifestations linked to extensional active faults that act as channels for the ascent of deep thermal fluids.The present work describes detailed geochemical exploration aimed at evaluating the potential of the Lesbos Island geothermal resource. Exploration was carried out on the three sites (Kalloni-Stipsi, Petra-Argenos and Polichnitos) that have the most favourable hydrogeological and structural setting on the island.Hydrogeochemical data reveal the presence between Kalloni and Stipsi of a shallow thermal aquifer with temperatures below boiling point, which coincides with carbon dioxide and temperature anomalies in the overlying groundwaters. All the thermal waters in the study area have fairly similar physico-chemical features; their geochemical temperature is in the range 115–125°C. While low-medium enthalpy geothermal resources are relatively abundant in the three selected areas of Lesbos Island, the presence of a high enthalpy geothermal system is still the subject of debate.     

Thermal and geochemical structure of the Uenotai geothermal system, Japan

The Uenotai geothermal area is located in southern Akita prefecture of northern Honshu Island. The Uenotai geothermal system is a liquid-dominated system with a central zone of aquifer boiling. The two-phase reservoir has evolved from liquid in the natural state due to exploitation. Gas composition of the vapor phase in the reservoir is nearly in equilibrium and correlates with the vapor fraction in the reservoir and with discharging steam quality. The marginal part of the Uenotai system has cooled with the drop in ground-water level. The chemical characteristics of the geothermal water indicate mixing of the immature high Cl source water with conductively heated or steam-heated shallow water or surface water, as well as boiling and steam gain.     

Modeling studies for production potential of Chingshui geothermal reservoir

The Chingshui geothermal power plant was decommissioned in 1993 due to a continued decline in production. Although some geothermal exploration and field investigation had been exercised, the production potential of the reservoir is still not well understood. In this paper, numerical modeling approaches for characterization of the geothermal reservoir, investigation of reservoir production performance, and evaluation of exploitation scheme design are presented. At first, a site-scale refined grid numerical model was developed for simulating the natural state of Chingshui geothermal reservoir. Through the model, the production potential of the geothermal reservoir was estimated and the availability of water resources was assessed. We further built production model to simulate the production history during 1981–1993. From the production model, we can conclude that the abnormal drop of the reservoir production capacity is mainly caused by carbonate scaling. Potential production schemes with different reinjection designs were evaluated through the model. Simulation results indicated that a sustainable hot water production capacity of Chingshui geothermal reservoir is about 200 t/h without reinjection, and 300 t/h or even higher with reinjection which is enough for a 3 MWe power plant. The simulation results indicate that reinjection provides an effective approach for maintaining reservoir pressure during hot water/steam production.     

Temperature evaluation of the Bugok geothermal system, South Korea

Using a variety of chemical geothermometers and statistical analysis, we estimate the temperature of a possible deeper geothermal reservoir at Bugok, Southern Korea. Shallow thermal aquifers (down to about 400 m depth) are under exploitation in this area; the temperatures (up to 78 °C) of the produced fluids are the highest found in South Korea. Based on hydrochemical data and occurrence, the groundwaters at Bugok can be classified under three groups: Na-SO₄ thermal groundwaters (CTGW) occurring in the central (about 0.24 km²) part of the area; Ca-HCO₃ cold groundwater (SCGW) found in shallow peripheral parts of the CTGW; the intermediate-type groundwater (STGW). The CTGW type is typical of the Bugok thermal waters; they have the highest discharge temperatures and contain very high concentrations of Na (75.1–101.0 mg/L), K (2.9–6.9 mg/L) and SiO₂ (62.0–84.5 mg/L) and are rich in sulfates.The major ion composition of the CTGW suggests that these waters are in partial equilibrium with rocks at depth. The application of various alkali-ion geothermometers yields temperature estimates in the 88–198 °C range for the thermal reservoir. Multiple-mineral equilibrium calculations indicate a similar but narrower temperature range (from about 100 to 155 °C). These estimates for CTGW are significantly higher than the measured discharge temperatures. Considering the heat losses occurring during the ascent of the waters, one can infer the presence of a deeper (around 1.8 km) thermal reservoir in the Bugok area that could be developed for district heating or other direct applications of geothermal heat.     

Groundwater changes in the Wairakei–Tauhara geothermal system

Overlying most of the Wairakei–Tauhara geothermal system is a sequence of shallow aquifers. Some of these are groundwaters heated by fluids escaping from the deep geothermal reservoir, others are cold; some are extensive and confined, others are perched, unconfined and of limited extent. Water levels of some of the shallow aquifers have declined at various times throughout the history of geothermal development. One of the mechanisms causing the observed changes is an increase in groundwater downflow into depressurizing steam zones. This occurs through natural fractures that once channelled the ascent of two-phase fluids. Another important mechanism is internal downflow within boreholes connecting permeable aquifers that were otherwise hydrologically separated. These cool downflows have had a significant impact on the deeper reservoir, and also affect the discharge characteristics of surface thermal features. However, in the Southern Tauhara sector, outflows of dilute chloride fluid from shallow thermal aquifers have not been affected by the deep pressure drawdown. This could be attributed to liquid or two-phase upflow continuing to recharge the chloride component in the south, or to significant quantities of residual chloride in a large geothermally heated aquifer still masking recharge depletion that has occurred during the last 50 years.     

Hot springs and geothermal gradients in northern Thailand

Chemical geothermometry of hot springs in northern Thailand indicates that many have reservoir temperatures in excess of 150°C and some in excess of 180°C. Measurements of temperatures in abandoned oil wells in Fang Basin indicate geothermal gradients of 70 – 130 mK/m. The high geothermal gradient may be the result of extensional tectonics in northern Thailand, caused indirectly by sea-floor spreading in the Andaman Sea. Relatively high reservoir temperatures and shallow reservoir depths suggest that hot spring areas in northern Thailand may be potential sources of geothermal energy.     

An updated numerical model of the Larderello–Travale geothermal system, Italy

Larderello–Travale is one of the few geothermal systems in the world that is characterized by a reservoir pressure much lower than hydrostatic. This is a consequence of its natural evolution from an initial liquid-dominated to the current steam-dominated system. Beneath a nearly impermeable cover, the geothermal reservoir consists of carbonate-anhydrite formations and, at greater depth, by metamorphic rocks. The shallow reservoir has temperatures in the range of 220–250 °C, and pressures of about 20 bar at a depth of 1000 m, while the deep metamorphic reservoir has temperatures of 300–350 °C, and pressures of about 70 bar at a depth of 3000 m. The 3D numerical code “TOUGH2” has been used to conduct a regional modeling study to investigate the production mechanism of superheated steam, the interactions between the geothermal field and the surrounding deep aquifers, and the field sustainability. All the available geoscientific data collected in about one century of exploration and exploitation have been used to provide the necessary input parameters for the model, which covers an area (4900 km²) about 10 times wider than the Larderello–Travale geothermal field (400 km²). The numerical model explains the origin of the steam extracted in about one century of exploitation and shows that, at the current level, the production is sustainable at least for the next 100 years.     

Numerical modeling of the Momotombo geothermal system, Nicaragua

The Momotombo geothermal field, Nicaragua, has been under commercial production since 1983; presently it has an installed capacity of 77 MWe. A three-dimensional, porous-media numerical model of the system was developed and calibrated using the code iTOUGH2, and was utilized to study the response of the geothermal reservoir under different fluid production and injection scenarios. Satisfactory agreements were obtained between measured and computed discharge enthalpies and flow rates for most of the shallow wells. The model also qualitatively reproduced the pressure drawdowns measured in selected wells.     

Simple modelling of the effects of exploitation on hot springs, Geyser Valley, Wairakei, New Zealand

Many hot springs in Geyser Valley declined and ceased flowing during well testing (1950–1958) or the early stages of development (1958–1964) of the Wairakei field, North Island, New Zealand. A simple model that considers mixing, geothermal reservoir pressure, and Darcy flow is applied to analyse the changes in flow rate and chemistry of spring SP18. The model suggests that prior to the testing the observed flow rate of 3 l/s was comprised of about 2.1 l/s of water from the deep reservoir and about 0.9 l/s from shallow groundwater. As the pressure in the reservoir decreased as a result of exploitation, the reservoir component declined but the groundwater component remained near constant until the spring had nearly ceased flowing. The model reasonably predicts the flow cessation date for four other springs (Group A) but poorly predicts that for seven other springs (Group B). These two groups of springs do not appear to be clustered at the ground surface. The pressure–elevation relationship suggests that the Group A springs feed from within the shallower Waiora Formation and Group B springs from within the deeper Wairakei Ignimbrite Formation. The observation that declines of chloride concentrations in Group A springs occurred at least three years before those in springs of Group B is consistent with the finding that Group A springs have a shallower feeder depth.     

System dynamics modelling of spring behaviour in the Orakeikorako geothermal field,New Zealand

System dynamics software STELLA is used to obtain mass and thermal balances of a spring in the Orakeikorako geothermal field, New Zealand, based on field measurements of water level, barometric pressure, rainfall and spring temperature. The model identifies the interactions of the principal influences on spring behaviour of rainfall, groundwater, geothermal steam and barometric pressure. The geothermal steam inflow estimated from the model, of about 0.022 kg/s, confirms the existence of a weak hydraulic connection with a deeper geothermal reservoir.     

Geochemical study of fluids on Lesbos island, Greece

Sixty-five water samples and seven associated gas samples have been collected on Lesbos island. The lithology and structural setting have resulted in two main types of hydrological circulation: a shallow circulation hosting low-salinity cold waters and a deeper one, hosting high-salinity hot waters that often emerge in thermal springs near the coast. The cold waters are characterized by Ca(Mg)-HCO₃(SO₄) composition, while the thermal waters generally have an Na-Cl composition. The chemical features of the former can be explained by their circulation in the ophiolite-bearing phyllitic basement and volcanic rocks. Waters circulating in the ultramafic layers of the basement are richer in Mg than the waters whose circulation is mainly within marble levels or volcanic rocks. The Na-Cl thermal waters are characterized by salinities ranging from 1910 to 35,700 mg/kg. As indicated by previous hydrogeochemical and isotopic studies, the Na-Cl composition of the thermal waters on Lesbos is the result of mixing between shallow meteoric waters and marine waters. While interacting with the minerals of the geothermal reservoir, the saline waters retain the Na/Cl sea water ratio but become enriched in Ca²⁺ and depleted in Mg²⁺ with respect to sea water.Processes of hydrothermal alteration at depth are activated by a gas phase enriched in CO₂, which reaches the geothermal reservoir by rising along the deep fractures of the basement. Thermodynamic calculations based on hydrothermal alteration processes occurring at the estimated temperatures of the geothermal reservoir (about 120 °C) indicate that the thermal waters of Lesbos are in equilibrium with talc and dolomite.     

San Jacinto-Tizate geothermal field,Nicaragua: Exploration and conceptual model

The results of the exploration of the San Jacinto-Tizate geothermal field during 1992–1995 included geological, hydrogeological and geophysical investigations (magnetotelluric, frequency soundings, subsurface temperature, and soil-gas surveys), and the drilling and testing of seven deep wells (728–2339 m). The geothermal field, located within a composite volcano-tectonic depression, can be divided into two main areas: San Jacinto and Tizate. The San Jacinto area shows evidence of a high-temperature (250–300°C) fossil geothermal system that at present has reservoir temperatures in the 180–190°C range. In the Tizate area there is an active geothermal system with temperatures of 250–285°C. An upflow zone with an excess pressure gradient exists in the central part of this area. Two hydraulically connected reservoirs exist: a shallow one at 550–1200 m depth, and a deeper one below 1600m. Two-phase conditions exist in the upper part of the shallow reservoir. Production tests demonstrate the commercial potential of both Tizate reservoirs.     

Reservoir engineering assessment of Dubti geothermal field,Northern Tendaho Rift,Ethiopia

Following on from surface exploration surveys performed during the 1970s and early 1980s, exploration drilling was carried out in the Tendaho Rift, in Central Afar (Ethiopia), from October 1993 to June 1995. Three deep and one shallow well were drilled in the central part of the Northern Tendaho Rift to verify the existence of a geothermal reservoir and its possible utilisation for electric power generation. The project was jointly financed by the Ethiopian Ministry of Mines and Energy and the Italian Ministry for Foreign Affairs. Project activities were performed by the Ethiopian Institute of Geological Surveys and Aquater SpA. The main reservoir engineering data discussed in this paper were collected during drilling and testing of the above four wells, three of which are located inside the Dubti Cotton Plantation, in which a promising hydrothermal area was identified by surface exploration surveys. Drilling confirmed the existence of a liquid-dominated shallow reservoir inside the Dubti Plantation, characterised by a boiling-point-for-depth temperature distribution down to about 500 m depth. The main permeable zones in the Sedimentary Sequence, which is made up of lacustrine deposits, are located in correspondence to basalt lava flow interlayerings, or at the contact between volcanic and sedimentary rocks. At depth, the basaltic lava flows that characterise the Afar Stratoid Series seem to have low permeability, with the exception of fractured zones associated with sub-vertical faults. Two different upflows of geothermal fluids have been inferred: one flow connected to the Dubti fault feeds the shallow reservoir crossed by wells TD-2 and TD-4, where a maximum temperature of 245 °C was recorded; the second flow seems to be connected with a fault located east of well TD-1, where the maximum recorded temperature was 270 °C. A schematic conceptual model of the Dubti hydrothermal area, as derived from reservoir engineering studies integrated with geological, geophysical and geochemical data, has been tested by numerical simulation, using the TOUGH2/EWASG code. Preliminary simulations, using a simple 3-D numerical model of the Dubti fault area, showed that measured temperature and pressure distribution, as well as evaluated non-condensable gas pressure at reservoir conditions, are compatible with the rise of geothermal fluid, at about 290 °C, along the sub-vertical Dubti fault from beneath the surface manifestations DB1, DB2 and DB3 located at the south-eastern end of the fault. According to the proven shallow field potential, development of this field could meet the predicted electricity requirements of Central Afar until the year 2015.     

油区地热能开发新工艺的研究

在油区地热资源开发过程中,应用井内换热和尾水回灌的对井供热系统可以提高地热资源开发效益、避免地热水资源浪费和实现可持续发展。文章使用了地热储二维分布模型和井筒传热模型．通过对地热储的数值模拟预测了尾水回灌过程中地热储的压力响应和冷却效应,通过对井筒的热力计算得到了井内换热器的最优设计方法。计算结果表明,地热储对尾水回灌的压力响应非常迅速,压力场能很快达到稳定;回灌冷水的影响区集中在回灌井的周围,冷区半径增长的速度越来越慢。井内换热器在最佳设置深度时,经济效益达到最大值;井内换热器设置深度不变,增大载热水的流量经济效益将增加,但是获得的热能温度下降。     

Forecast and evaluation of hot dry rock geothermal resource in China

Utilizing information from plate tectonics characteristics, volcanic activities, and geothermal anomaly, this paper identifies areas where hot dry rock (HDR) may exist as potential geothermal resource in China. Further investigations are also carried out in the paper based on results from regional tectonics, volcanic geology and lithology, as well as data from geothermal displays, geochemistry, geophysics, and shallow borehole temperature measurements. The study reveals several promising areas of HDR geothermal resource in China, including Tengchong of Yunnan province, Qiongbei of Hainan province, Changbaishan of Jilin province, Wudalianchi of Heilongjiang province, and the Southern Tibet area. A 3D static heat conduction model was developed to study the underground temperature gradient characteristics of the Rehai geothermal field in Tengchong and the Yangbajing geothermal field in Tibet. The model adopted is a geological block 10 km deep from the ground surface and 50 km wide in each of the horizontal directions (2500 km² area). The numerical simulation results in evaluations on the quantities of the HDR geothermal resource in Rehai and Yangbajing geothermal fields. The paper shows that there is abundant HDR geothermal resource with large exploitation value in China. If developed with a power capacity of 1×10⁸ kW, the Rehai and Yangbajing fields along would be able to generate electricity for 1560 years.     

Geochemistry in geothermal exploration

Techniques based on the variations in composition of water, gas and stable isotopes in the liquid and gas phases of the geothermal fluids have been applied for some time now in the major geothermal fields and are now also used regularly in geothermal exploration. There are numerous processes capable of modifying isotopic composition after infiltration of water from the surface, such as water-rock exchanges, formation of secondary minerals and exchange with the gaseous phase (CO₂ and H₂S). During ascent to the surface, the two main processes are steam separation and dilution and mixing with shallower waters. This paper also deals with the chemical characteristics of the waters, their classification and the water-rock interaction producing hydrothermal alteration. During exploration the chemical and isotopic geothermometers represent a unique method for investigating the deep system. The choice of geothermometer and interpretation of geothermometric data are two crucial steps in geothermal exploration. Finally, the paper discusses the geochemistry of gas mixtures, especially the origin of the gas species and the main chemical reactions that produce semi-empirical geothermometers and some recent non-empirical geothermometers based on models of a two-phase system in the reservoir. Gas-geothermometers can be developed to calculate the reservoir temperature for natural manifestations.     

Heat flow and deep temperatures in the Chaves geothermal system,Northern Portugal

The geological, geoelectrical, geochemical and temperature data related to the Chaves geothermal system have been integrated to obtain a better understanding of the Chaves basin. Geoelectrical surveys carried out in the basin reveal a low-resistivity zone (10 ohm m), associated with a shallow geothermal reservoir, in the central part of the graben, bounded by higher-resistivity rocks. The top of this zone varies between 400 and 200 m and its maximum thickness (1600 m) is located at the centre of the basin. Thermal models for the Chaves basin and for the region are presented using the structure obtained by geoelectrical methods and a mean heat flow value of 95 mW m^-2 derived from borehole measurements. The heat transfer takes place mainly by conduction, except near the faults, where convective flow is important. The medium is considered dishomogeneous and there is a great thermal conductivity contrast between the sediments in the basin and the surrounding rocks. The results obtained for the Chaves basin show that the mean temperature value in the shallow geothermal reservoir is 62°C. The maximum temperature value predicted to the bottom of this reservoir is 95°C. A regional forced convective-circulation model is suggested based on geomorphological, geochemical, isotopic data and to rmal models.     

Engineering evaluation of geothermal reservoirs

Reservoir engineering evaluation of geothermal systems attempts to provide answers on the extent of the reserves, their probable longevity and the deliverability and production rate of the reservoir. Economic decisions on the desirability of the exploitation of the particular reservoir hinge on these findings. This paper presents a set of calculational procedure and thinking sequences available to the geothermal reservoir engineer that would aid in an appropriate management decision. An exploitable geothermal system consists of a fluid as well as a heat reservoir. Since much of the heat is stored in the confining rock, reinjection strategies are outlines for the efficient mining of heat.