The geothermal system of Ischia Island (southern Italy): Critical review and sustainability analysis of geothermal resource for electricity generation

In this paper we analyze the main available data related to the geothermal system of Ischia Island, starting from the first geothermal exploration in 1939. Our aim is to define a conceptual model of the geothermal reservoir, according to geological, geochemical, geophysical and stratigraphic data. In recent times, the interest on geothermal exploitation for electricity generation in Italy is rapidly increasing and the Ischia Island is one of the main targets for future geothermal exploitation. Nowadays, one of the main economic resources of the island is the tourism, mainly driven by the famous thermal springs; so, it is crucial to study the possible interaction between geothermal exploitation and thermal spring activities. To this aim, we also analyze the possible disturbance on temperature and pressure in the shallow geothermal reservoir, due to the heat withdrawal for electric production related to small power plant size (1–5 MWe). Such analysis has been performed by using numerical simulations based on a well known thermofluid-dynamical code (TOUGH2^®). Obtained results show that such geothermal exploitation generates a perturbation of temperature and pressure field which, however, is confined in a small volume around the well. At shallow level (0–100 m) the exploitation does not produce any appreciable disturbance, and can be made compatible with thermal spring exploitation. Moreover, such results are crucial both for the evaluation of volcanological processes in the island and for the general assessment of geothermal resource sustainability.     

Ecological characteristics and management of geothermal systems of the Taupo Volcanic Zone,New Zealand

New Zealand has an array of geothermal systems with distinctive ecological features, with many occurring in the Taupo Volcanic Zone in the Central North Island. Associated with these geothermal features are characteristic geophysical and geochemical components, and distinctive terrestrial and aquatic ecosystems with many attributes that are common across a range of the biotic groups. Zonation amongst vegetation communities is closely related to soil temperature and these associations generally occur in a predictable sequence along the soil temperature gradient. Similarly, clear distinctions in aquatic flora and fauna occur longitudinally downstream from the source of thermal springs and vertically on geyser mounds.     

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.     

Rapid estimation of geothermal coverage by district-heating systems

An algorithm is developed for the quick estimation of the fraction of annual heat-load that is covered by geothermal energy in district-heating systems. The cases of direct use, heat exchange and heat-pump application are distinguished. The algorithm proved accurate enough for pre-feasibility study purposes and for preliminary considerations with regard to the management of geothermal energy. It is demonstrated that for coverage estimation, detailed heat-supply calculations can be substituted by the use of degree-days, which are quite often available for various base temperatures. Furthermore, the analytic expressions of the algorithm allow the quantitative evaluation of several parameters. Specifically, for low temperature geothermal fluids the coverage is an increasing second-degree function of the temperature of the fluids. Increasing network return temperatures has also a decreasing effect on the coverage. The required geothermal-flow is inversely proportional to the heat exchanger effectiveness, for the same coverage. Application of a heat pump results in a higher coverage that should otherwise require a higher temperature fluid. It is consequently defined as an equivalent geothermal-temperature increase – due to the heat pump – that is proved to be proportional to the compressor capacity. Last, the network supply-temperatures have a rather secondary effect upon the coverage, affecting the flow through the network and so the effectiveness of the heat exchanger.     

Numerical modeling of gland seal erosion in a geothermal turbine

Excessive erosion of the low-pressure rotor end gland seal of a 25 MWe geothermal turbine produced a partial loss of turbine vacuum that degraded cycle efficiency. This study used computational fluid dynamics (CFD) to identify the causes of erosion and the optimal steam seal system flow conditions for reducing the erosion problem. The predictions were based upon a numerical calculation using a commercial CFD code (Adapco Star-CD) to model the rotor end gland seal with a steam flow containing hard solid particles. The results confirmed that flow conditions play a major role in rotor gland seal erosion. By changing steam seal flow pressures to vary flow, it was confirmed that there is a threshold seal flow condition below which erosion does not occur, or is minimized. Optimizing the rotor end gland seal supply pressure and intercondenser pressure reduced the turbulent flow kinetic energy by 49%, with a corresponding decrease in the erosion rate of the rotor gland seal surface. The erosion rate is related directly to the particle velocity and turbulent flow kinetic energy. Recommendations are provided for adjusting the rotor end gland seal system to avoid erosion.     

Thermodynamic analysis of a hybrid geothermal heat pump system

A thermodynamic analysis of a hybrid geothermal heat pump system is carried out. Mass, energy, and exergy balances are applied to the system, which has a cooling tower as a heat rejection unit, and system performance is evaluated in terms of coefficient of performance and exergy efficiency. The heating coefficient of performance for the overall system is found to be 5.34, while the corresponding exergy efficiency is 63.4%. The effect of ambient temperature on the exergy destruction and exergy efficiency is investigated for the system components. The results indicate that the performance of hybrid geothermal heat pump systems is superior to air-source heat pumps.     

Piping network design of geothermal district heating systems: Case study for a university campus

Geothermal district heating system design consists of two parts: heating system and piping network design. District heating system design and a case study for a university campus is given in Yildirim et al. [1] in detail. In this study, piping network design optimisation is evaluated based on heat centre location depending upon the cost and common design parameters of piping networks which are pipe materials, target pressure loss (TPL) per unit length of pipes and installation type. Then a case study for the same campus is presented.     

Hyperspectral detection of geothermal system-related soil mineralogy anomalies in Dixie Valley, Nevada: a tool for exploration

Hyperspectral data analysis has been applied to the mapping of soil anomalies that may be related to present or past geothermal systems. Anomalous accumulations of certain soil minerals can indicate buried geologic structures and possible zones of elevated permeability. Hyperspectral data can be used to map these anomalies as part of geothermal exploration activities. The study area for this project was northern Dixie Valley, Nevada, which is host to a structurally-controlled deep-circulation hydrothermal convection system. Advanced visible and Infrared Imaging Spectrometer (AVIRIS) airborne hyperspectral imagery was used. Both supervised and unsupervised spectral unmixing methods were tested to separate minerals from other components in the image. Both methods produced useful spectral end-members leading to the detection of anomalous soil minerals that may be related to the geothermal system and buried geologic structures.     

Conjugate fluid flow and kinetics modeling for heat exchanger fouling simulation

Thermal treatment of fluid foods represents a major unit operation in the food industry, to ensure the product's safety and quality features. But during the thermal treatments of such sensible fluids in common plate heat exchangers, food constituents such as proteins can be thermally damaged and precipitated to form fouling that greatly affect the treatment efficiency and alter the product's desired features.Computational Fluid Dynamics simulations can then be successfully exploited, bringing forth temperature and velocity information that yield for deposit distributions when coupled to biochemical notations for thermal denaturation of fluid constituents.The present work exploits such modeling for a single-channel heat exchanger during pasteurization of milk. The model enforces a conjugate system of differential equations to a heat exchanger's corrugated plate to combine flow, heat transfer and local transport of β-lactoglobulin. A preliminary computation has been performed that could be applied to geometry optimization (different corrugation shape and orientation) and for a variety of biochemically evolutive products.     

Influence of different fracture morphology on heat mining performance of enhanced geothermal systems based on COMSOL

Formation of enhanced geothermal systems (EGS) is the necessary approach to obtain geothermal energy efficiently. In-situ stress, nature of reservoir physical properties and fracturing methods will affect the artificial fracture morphology after reservoir stimulation. A three-dimension thermal coupled seepage model of fractured media was established to simulate the influence of fracture morphology on heat mining performance of EGS, considering the pressure- and temperature-dependent physical properties of working medium. The results indicate that formation of complex fracture network is favorable for heat mining. Production mass flow in Case1 with complex fracture network enhances nearly 2.5 times comparing to the unenhanced model at exploitation beginning. The total net energy rate will up to 44 MW and be maintained above 10 MW for 5 years. The system impedance can be effectively reduced, however the sustainable heat mining duration decreased to 30 years. The increase in length and number of branch fractures is expected. While increasing the width of branch fractures deliberately has little effect on the exploitation of EGS. Finally, we investigate the adaptability of employing supercritical CO₂ in EGS with complex fracture network. Production mass flow will be enhanced 3–5 times compared with water, but the stability is poor, total net energy rate decrease from 90 MW to 3 MW over the 10-year operation period.     

Evaluation of heat exchange rate of GHE in geothermal heat pump systems

Total thermal resistance of ground heat exchanger (GHE) is comprised of that of the soil and inside the borehole. The thermal resistance of soil can be calculated using the linear source theory and cylindrical source theory, while that inside the borehole is more complicated due to the integrated resistance of fluid convection, and the conduction through pipe and grout. Present study evaluates heat exchange rate per depth of GHE by calculating the total thermal resistance, and compares different methods to analyze their similarities and differences for engineering applications. The effects of seven separate factors, running time, shank spacing, depth of borehole, velocity in the pipe, thermal conductivity of grout, inlet temperature and soil type, on the thermal resistance and heat exchange rate are analyzed. Experimental data from several real geothermal heat pump (GHP) applications in Shanghai are used to validate the present calculations. The observations from this study are to provide some guidelines for the design of GHE in GHP systems.     

Transient heat conduction in an infinite medium subjected to multiple cylindrical heat sources: An application to shallow geothermal systems

In this paper, we introduce analytical solutions for transient heat conduction in an infinite solid mass subjected to a varying single or multiple cylindrical heat sources. The solutions are formulated for two types of boundary conditions: a time-dependent Neumann boundary condition, and a time-dependent Dirichlet boundary condition. We solve the initial and boundary value problem for a single heat source using the modified Bessel function, for the spatial domain, and the fast Fourier transform, for the temporal domain. For multiple heat sources, we apply directly the superposition principle for the Neumann boundary condition, but for the Dirichlet boundary condition, we conduct an analytical coupling, which allows for the exact thermal interaction between all involved heat sources. The heat sources can exhibit different time-dependent signals, and can have any distribution in space. The solutions are verified against the analytical solution given by Carslaw and Jaeger for a constant Neumann boundary condition, and the finite element solution for both types of boundary conditions. Compared to these two solutions, the proposed solutions are exact at all radial distances, highly elegant, robust and easy to implement.     

Geochemical modeling of the near-surface hydrothermal system beneath the southern moat of Long Valley Caldera, California

Geochemical reaction path and mass balance modeling techniques were used to test the hypothesis that an eastwardly flowing plume of thermal water in the southern moat of the Long Valley caldera system reacts with hydrothermally altered intra-caldera tuffs and mixes with non-thermal groundwater. Our conceptual model is based on hypotheses in the literature and published geochemical and petrologic data. Mixing of thermal and non-thermal waters and reaction with wall rock were simulated using the reaction path code EQ3/6. Mass balance calculations were conducted to estimate the extent of water-rock interaction between the intra-caldera tuffs and fluids. A mixing ratio of 82% thermal and 18% non-thermal water reacting with altered tuff minerals closely matches Casa Diablo fluid compositions and minerals observed in petrographic studies. Results of this study show that the mineralogy and fluid chemistry observed in the shallow reservoir at Long Valley caldera are formed in an open system. Further, calcite precipitated in the system serves as a sink for high levels of CO₂ generated by the deeper magmatic system. Our study serves as an example that processes acting in a geothermal system can be effectively quantified using geochemical modeling and mass balance calculations.     

Study of calcite scaling and corrosion processes in geothermal systems

Prediction of scaling in geothermal systems is important in order to be able to take preventive action. Scaling and corrosion processes are correlated in geothermal systems so that the factors influencing both processes are discussed in this paper. Pilot-scale equipment to model dissolution and scaling, and investigate corrosion, has been constructed. The construction of the pilot-scale equipment and the design of the experiments were based on the results of earlier laboratory model experiments. In order to check the measurement results we used a computer programe that was developed earlier, and which is suitable for the calculation of the equilibrium solubilities. Based on the results of the equilibrium experiments, we developed a modelling method for scaling in the pilot-scale equipment, and studied corrosion during scaling on the surfaces of five different structural materials.     

Numerical computation of fluid flow and heat transfer in microchannels

Three-dimensional fluid flow and heat transfer phenomena inside heated microchannels is investigated. The steady, laminar flow and heat transfer equations are solved using a finite-volume method. The numerical procedure is validated by comparing the predicted local thermal resistances with available experimental data. The friction factor is also predicted in this study. It was found that the heat input lowers the frictional losses, particularly at lower Reynolds numbers. At lower Reynolds numbers the temperature of the water increases, leading to a decrease in the viscosity and hence smaller frictional losses.     

Sustainability analysis of the Ahuachapán geothermal field: management and modeling

The Ahuachapán geothermal field (AGF) is located in north western El Salvador. To date, 53 wells (20 producers and 8 injectors) have been drilled in the Ahuachapán geothermal field and the adjacent Chipilapa area. Over the past 33 years, 550 Mtonnes have been extracted from the reservoir, and the reservoir pressure has declined by more than 15 bars. By 1985, the large pressure drawdown due to over-exploitation of the resource reduced the power generation capacity to only 45 MW_e. Several activities were carried out in the period 1997–2005 as part of “stabilization” and “optimization” projects to increase the electric energy generation to 85 MW_e, with a total mass extraction of 850 kg/s.     

The geothermal characteristics of the ground and the potential of using ground coupled heat pumps in Cyprus

The ground can be used as an energy source, an energy sink, or for energy storage and for this reason ground characteristics should be available. Therefore the purpose of this paper is to present the recorded ground temperatures at eight representative sites of Cyprus, in relation to depth, time of year, geology and altitude and discuss the efficiency of Ground Coupled Heat Pumps. The ground temperature was recorded for a period of one year, from October 2009 to October 2010. According to the results obtained, in several locations in Cyprus the surface zone reaches a depth of 0.5 m. The shallow zone penetrates to 7-8 m and there after the deep zone follows in which the temperature remains constant throughout the year with a range between 18 and 23 °C. For the eight boreholes, additional geothermal data were also recorded like the type of ground and thermal conductivities of the various geologic layers. The data collected clearly indicate that there is a potential for the efficient use of Ground Coupled Heat Pumps (GCHPs) in Cyprus leading to significant savings in heating and cooling energy consumption.     

Use of solar assisted geothermal heat pump and small wind turbine systems for heating agricultural and residential buildings

The main objective of the present study is twofold: (i) to analyze thermal loads of the geothermally and passively heated solar greenhouses; and (ii) to investigate wind energy utilization in greenhouse heating which is modeled as a hybrid solar assisted geothermal heat pump and a small wind turbine system which is separately installed in the Solar Energy Institute of Ege University, Izmir, Turkey. The study shows 3.13% of the total yearly electricity energy consumption of the modeled system (3568 kWh) or 12.53% of the total yearly electricity energy consumptions of secondary water pumping, brine pumping, and fan coil (892 kWh) can be met by using small wind turbine system (SWTS) theoretically. According to this result, modeled passive solar pre heating technique and combined with geothermal heat pump system (GHPS) and SWTS can be economically preferable to the conventional space heating/cooling systems used in agricultural and residential building heating applications if these buildings are installed in a region, which has a good wind resource.     

Inverse modeling and forecasting for the exploitation of the Pauzhetsky geothermal field,Kamchatka, Russia

A three-dimensional numerical model of the Pauzhetsky geothermal field has been developed based on a conceptual hydrogeological model of the system. It extends over a 13.6-km² area and includes three layers: (1) a base layer with inflow; (2) a geothermal reservoir; and (3) an upper layer with discharge and recharge/infiltration areas. Using the computer program iTOUGH2 [Finsterle, S., 2004. Multiphase inverse modeling: review and iTOUGH2 applications. Vadose Zone J. 3, 747–762], the model is calibrated to a total of 13,675 calibration points, combining natural-state and 1960–2006 exploitation data. The principal model parameters identified and estimated by inverse modeling include the fracture permeability and fracture porosity of the geothermal reservoir, the initial natural upflow rate, the base-layer porosity, and the permeabilities of the infiltration zones. Heat and mass balances derived from the calibrated model helped identify the sources of the geothermal reserves in the field. With the addition of five make-up wells, simulation forecasts for the 2007–2032 period predict a sustainable average steam production of 29 kg/s, which is sufficient to maintain the generation of 6.8 MWe at the Pauzhetsky power plant.