Experimental–numerical study of heat flow in deep low-enthalpy geothermal conditions

This paper presents an intensive experimental–numerical study of heat flow in a saturated porous domain. A temperature and a flow rate range compared to that existing in a typical deep low-enthalpy hydrothermal system is studied. Two main issues are examined: the effect of fluid density and viscosity on heat flow, and the significance and effect of thermal dispersion. Laboratory experiments on a saturated sand layer surrounded by two impermeable clay layers, subjected to different flow rates under cold and hot injection scenarios, and for both vertical and horizontal flow directions, are conducted. A temperature range between 20 °C and 60 °C is studied. The finite element method is utilized to analyze the experimental results. Backcalculations, comparing the numerical results to the experimental results, are conducted to quantify the magnitude of thermal dispersion. A constitutive model describing thermal dispersion in terms of fluid density, viscosity and pore geometry, taking into consideration different injection scenarios, is developed. This study demonstrates the importance of taking the variation of formation water density and viscosity with temperature into consideration in predicting the lifetime of deep low-enthalpy geothermal systems. It shows that if ignored, the lifetime of a system with hot injection will be overestimated, and that with cold injection, will be underestimated.     

Modelling of fluid flow and heat transfer to assess the geothermal potential of a flooded coal mine in Lorraine,France

The flooding of the Lorraine coal mines (France), representing a huge reservoir of about 154 × 10⁶ m³, began in June 2006. After attaining thermal equilibrium with the surrounding rocks, the water temperature in the deepest parts is expected to reach 55 °C, giving the opportunity for the extraction of low-enthalpy geothermal waters that may be suitable for district heating purposes. We present some numerical modelling results of the thermally driven convective flow in an open vertical shaft and in the entire mine reservoir. A dual permeability/porosity approach was used in the reservoir model, which includes open galleries and vertical shafts, coal panels backfilled with sand, and intact rock masses. Two scenarios of heat extraction with different flow regimes were investigated. A sensitivity analysis shows that the temperature decline in the production zone is highly dependent on the permeability of the surrounding porous rocks. Larger permeabilities result in higher water temperatures at the production shaft due to greater inflows of warm water from those rock masses.     

Geophysics in geothermal prospecting

Geophysical prospecting of high temperature geothermal reservoirs aims at identifying either fluid trapping structures or anomalies related to the properties of the hydrothermal fluid and rock to fluid interactions. Two types of reservoir environments can be characterized: (i) sedimentary reservoirs when a carbonate reservoir is generally capped by a dominantly argillaceous, hydraulically impervious and thermally insulating cover, and (ii) volcanic and volcano-sedimentary reservoirs associated with hydrothermally altered areas. Based on the aforementioned exploration goals and reservoir settings, a wide spectrum of geophysical methods can be applied whose selection is largely commanded by local geological conditions and expected reservoir morphology. Major geophysical techniques are reviewed and their potential, as to geothermal reservoir prospecting issues, discussed.     

The geothermal project Den Haag: 3D numerical models for temperature prediction and reservoir simulation

The proposed Den Haag Zuidwest district heating system of the city of The Hague consists of a deep doublet in a Jurassic sandstone layer that is designed for a production temperature of 75 °C and a reinjection temperature of 40 °C at a flow rate of 150 m³ h⁻¹. The prediction of reservoir temperature and production behavior is crucial for success of the proposed geothermal doublet. This work presents the results of a study of the important geothermal and geohydrological issues for the doublet design. In the first phase of the study, the influences of the three-dimensional (3D) structures of anticlines and synclines on the temperature field were examined. A comprehensive petrophysical investigation was performed to build a large scale 3D-model of the reservoir. Several bottomhole temperatures (BHTs), as well as petrophysical logs were used to calibrate the model using thermal conductivity measurements on 50 samples from boreholes in different lithological units in the study area. Profiles and cross sections extracted from the calculated temperature field were used to study the temperature in the surrounding areas of the planned doublet. In the second phase of the project, a detailed 3D numerical reservoir model was set up, with the aim of predicting the evolution of the producer and injector temperatures, and the extent of the cooled area around the injector. The temperature model from the first phase provided the boundary conditions for the reservoir model. Hydraulic parameters for the target horizons, such as porosity and permeability, were taken from data available from the nearby exploration wells. The simulation results are encouraging as no significant thermal breakthrough is predicted. For the originally planned location of the producer, the extracted water temperature is predicted to be around 79 °C, with an almost negligible cooling in the first 50 years of production. When the producer is located shallower parts of the reservoir, the yield water temperatures is lower, starting at ≈76 °C and decreasing to ≈74 °C after 50 years of operation. This comparatively larger decrease in temperature with time is caused by the structural feature of the reservoir, namely a higher dip causes the cooler water to easily move downward. In view of the poor reservoir data, the reservoir simulation model is constructed to allow iterative updates using data assimilation during planned drilling, testing, and production phases. Measurements during an 8 h pumping test carried out in late 2010 suggest that a flow rate of 150 m³ h⁻¹ is achievable. Fluid temperatures of 76.5 °C were measured, which is very close to the predicted value.     

Shallow submarine and subaerial, low-enthalpy hydrothermal manifestations in Punta Banda, Baja California, Mexico: Geophysical and geochemical characterization

The low-enthalpy geothermal system at Punta Banda (NW Baja California Peninsula, Mexico) has been studied because it might provide heat to future desalination plants in the city of Ensenada. Utilization of subaerial, intertidal and submarine hot springs is evaluated based on geochemical and geophysical data. The results of the geochemical studies show that the geothermal fluids have a major meteoric water component because seawater is not present at the subaerial springs and hot wells. The highest estimated reservoir temperature (140 °C) calculated using a silica geothermometer corresponds to the Agua Caliente intertidal manifestation, a promising area also identified by geophysics. Geothermometric calculations applied to the computed composition of the thermal end member yield a reservoir temperature of 137 °C. Cl/B ratios indicate that the thermal fluids discharged by the intertidal vents and subaerial springs are similar, but they differ from those of submarine vents. Geoelectrical models depict an anomalous conductive trend from the La Jolla well to the Agua Caliente manifestation, suggesting the presence of a fault that allows upflow of hot water from depth. Lastly, integration of geochemical and geophysical data identified the best site for future exploration drilling at Punta Banda.     

Geothermal resources in Algeria

The geothermal resources in Algeria are of low-enthalpy type. Most of these geothermal resources are located in the northeastern of the country. There are more than 240 thermal springs in Algeria. Three geothermal zones have been delineated according to some geological and thermal considerations: (1) The Tlemcenian dolomites in the northwestern part of Algeria, (2) carbonate formations in the northeastern part of Algeria and (3) the sandstone Albian reservoir in the Sahara (south of Algeria). The northeastern part of Algeria is geothermally very interesting. Two conceptual geothermal models are presented, concerning the northern and southern part of Algeria. Application of gas geothermometry to northeastern Algerian gases suggests that the reservoir temperature is around 198 °C. The quartz geothermometer when applied to thermal springs gave reservoir temperature estimates of about 120 °C. The thermal waters are currently used in balneology and in a few experimental direct uses (greenhouses and space heating). The total heat discharge from the main springs and existing wells is approximately 642 MW. The total installed capacity from producing wells and thermal springs is around 900 MW.     

A proper spent nuclear fuel management strategy could enhance the continuity of nuclear power in the Spanish energy mix

This article presents two economic analyses performed with the Mariño model, which was specially designed to analyse the costs of different spent nuclear fuel (SNF) management strategies in the real Spanish context. These analyses are: (a) a Monte Carlo study for those strategies and (b) the effects of a longer operational lifetime for the Spanish nuclear power plants (NPPs) on the costs of spent nuclear fuel (SNF) management. For the first analysis, a triangular distribution for the different unitary costs was assumed and the data and assumptions from numerous studies were used to obtain the values required for the distribution. The second analysis was performed for the current official shutdown dates for the NPPs, and the results were compared to other operational lifetime scenarios. The main assumption for these scenarios was a progressive shutdown of the reactors, in order to avoid numerous shutdowns in a few years. These scenarios were proposed for 40 to 60 years of mean operational lifetime of the reactors. The results show that, for all scenarios analysed, the additional electricity production due to longer operational lifetimes compensate the extra costs caused by the larger amount of SNF to be managed. Additionally, for the current SNF management strategy, a progressive shutdown at 40 years of mean operational lifetime has shown to entail lower costs than the official shutdown scenario. However, a strategy without a centralised interim storage facility would be the most economically favourable one for all the scenarios analysed.     

Production capacity estimation by reservoir numerical simulation of northwest (NW) Sabalan geothermal field, Iran

A three dimensional numerical model of the northwest (NW) Sabalan geothermal system was developed on the basis of the designed conceptual model from available field data. A numerical model of the reservoir was expressed with a grid system of a rectangular prism of 12 km × 8 km with 4.6 km height, giving a total area of 96 km². The model has 14 horizontal layers ranging in thickness between 100 m to 1000 m extending from a maximum of 3600 to −1000 m a.s.l. Fifteen rock types were used in the model to assign different horizontal permeabilities from 5.0 × 10⁻¹⁸ to 4.0 × 10⁻¹³ m² based on the conceptual model.Natural state modeling of the reservoir was performed, and the results indicated good agreements with measured temperature and pressure in wells. Numerical simulations were conducted for predicting reservoir performances by allocating production and reinjection wells at specified locations. Three different exploitation scenarios were examined for sustainability of reservoir for the next 30 years. Effects of reinjection location and required number of makeup wells to maintain the specified fluid production were evaluated. The results showed that reinjecting at Site B, immediate adjacent to production zone, is most effective for pressure maintenance of the system.On the base of existing data and assumptions the reservoir can sustain producing fluid equivalent to 50 MW_e of electricity for more than 30 years. The reservoir can produce the maximum amount of fluids equivalent to 90-100 MW_e for only 5 years, but the production capacity decreases to 50 MW_e after 20 years of operation because of pressure and enthalpy drop. The reservoir can sustain 50 MW_e over 100 years that can be defined as a sustainable production level of the field.     

Optimization of downhole pump setting depths in liquid-dominated geothermal systems: A case study on the Balcova-Narlidere field,Turkey

Downhole pumps are being used increasingly more often in low-enthalpy geothermal wells. The depth at which these pumps are installed depends on the physical and chemical characteristics of the geothermal fluid, the production flow rate, and the reservoir pressure and permeability. In this study we have investigated the factors affecting pump setting depths in low-temperature, liquid-dominated geothermal systems and defined the relationship between flow rate and pressure drawdown based on multi-rate test results. The methodology proposed was applied to the Balcova-Narlidere geothermal field, Turkey. It was found that the most important parameters for determining the capacity and setting depth of a downhole pump were flow performance, non-condensable gas concentration, and temperature. The implementation of the methodology is illustrated.     

The Podhale geothermal reservoir simulation for long-term sustainable production

The Podhale geothermal system, located in the southern, mountainous part of Poland, is the most valuable reservoir of geothermal waters discovered in the country to date and the one with the highest capacities in Central and Eastern Europe. Over 20 years of continuous operation has proved its stable operating parameters – a small drop in pressure and an unnoticeable temperature change. Production of over 500 m³/h of geothermal water with an 86 °C wellhead temperature is current practise, while drilling a new production well and reconstruction of an injection well allows for production that may significantly exceed 600 m³/h. To utilize these vast resources, a binary power cycle for electricity and heat production is considered by group of researchers. The results of numerical modelling of heat extraction from the Podhale reservoir are presented in the article as a preliminary step to the detailed analysis of combined heat and power production through a binary power cycle.     

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.     

The Hijiori Hot Dry Rock test site, Japan: Evaluation and optimization of heat extraction from a two-layered reservoir

The Hijiori hot dry rock (HDR) system, Japan, consists of a shallow and a deep reservoir, both in fractured granitic rocks. During a long-term circulation test (LTCT) lasting approximately 18 months and which tested different fluid production scenarios, large changes were observed in output fluid temperatures, pressures, and flow rates. A multi-reservoir, numerical model of the Hijiori HDR system was developed using the finite element heat and mass (FEHM) transfer code and applied to simulations of the LTCT. The model reproduced the pressure, temperature, and flow data observed during the test. Based on the modeling study, it was concluded that most of the produced fluid came from the shallow reservoir, that the permeability of the deep reservoir changed during the initial part of the LTCT, and that the redistribution of injected water between the two reservoirs had little impact on the relative amounts of deep and shallow fluid production. After validating the model on the LTCT, it was used to optimize injection rates in both reservoirs. The model was also used in simulations of reservoir performance where an additional heat transfer surface area has been created in the subsurface through hydraulic fracturing.     

A generalized non-isothermal tank model for liquid dominated geothermal reservoirs

We present a generalized non-isothermal tank model for predicting the pressure and temperature behaviors of liquid dominated geothermal reservoirs. A geothermal system is represented by a single or multiple tanks. These tanks can represent the reservoir, multiple reservoirs, aquifers or any other component of a geothermal system. The mass and energy balance equations for each tank are solved jointly. One of the main advantages of the model is that only a small number of tanks are used for modeling which avoids over parameterization of a geothermal system and results in faster run times when compared with fully discretized numerical simulators. Synthetic examples are used for studying the effects of heat conduction on reservoir performance, an analysis of the location of injection wells, recovery times of depleted geothermal fields and the benefits of using temperature data for a better characterization of the geothermal system.     

Lumped-parameter models for low-temperature geothermal fields and their application

The behavior of low-temperature geothermal reservoirs under exploitation is simulated using analytical lumped-parameter models. These models consider the effects of fluid production and reinjection, as well as natural recharge, on the pressures (or water levels) of low-temperature, liquid-dominated geothermal systems. The computed responses for constant production/injection flow rates are given in the form of analytical expressions. Variable flow rate cases are modeled, based on the Duhamel's principle. Reservoir parameters are obtained by applying a weighted nonlinear least-squares estimation technique in which measured field data are history matched to the corresponding model response. By using history-matched models, the future performance of the reservoir can be predicted for different production/injection scenarios in order to optimize the management of a given geothermal system.We demonstrate the applicability of the models by simulating measured data from the Laugarnes geothermal field in Iceland, and the Balcova–Narlidere field in Turkey.     

Evaluation of the reconnaissance results in geothermal exploration using GIS

Geographical Information System (GIS) is used to determine the spatial association between geophysical and geological evidence and production zones in a well-known geothermal field (Los Azufres, Mexico). Surface observations in Los Azufres were used to delineate areas characterized by high permeability and hot fluid transport from the reservoir: main faults, superficial fracture density, surface manifestations, contacts with the most recent rhyolite domes, and low values in the apparent resistivity surveys. Three knowledge-driven models were constructed based on a conceptual model of the field: a hydrothermal system in rough terrain with secondary permeability. Boolean, Index Overlay and Fuzzy scheme models were proposed and the results obtained show a good correlation with the location of the producing and non-producing wells that have been drilled in the field. The results obtained are useful for well siting (Boolean and Low-Risk Fuzzy models) or for planning further detailed exploration (Index Overlay and High-Risk Fuzzy models).     

Energetic and exergetic performance evaluations of a geothermal power plant based integrated system for hydrogen production

In this paper, we propose an integrated system aiming for hydrogen production with by-products using geothermal power as a renewable energy source. In analyzing the system, an extensive thermodynamic model of the proposed system is developed and presented accordingly. In addition, the energetic and exergetic efficiencies and exergy destruction rates for the whole system and its parts are defined. Due to the significance of some parameters, the impacts of varying working conditions are also investigated. The results of the energetic and exergetic analyses of the integrated system show that the energy and exergy efficiencies are 39.46% and 44.27%, respectively. Furthermore, the system performance increases with the increasing geothermal source temperature and reference temperature while it decreases with the increasing pinch point temperature and turbine inlet pressure.     

西藏羊易EGS开发储层温度场与开采寿命影响因素数值模拟研究

西藏羊易地区具有丰富的地热能,单井开发潜力接近10 MW,对其深部热储进行EGS开采,可缓解西部能源紧缺问题。本文建立二维理想EGS开发模型,探讨深层地热开采过程中开采流量、注采方式、注入温度等参数对热储温度场分布及开采寿命的影响。基于羊易温度信息设计了12个数值模型,对比研究发现,开采流量对EGS开采的影响较大,为保证开采50年内的商业利用价值,最大开采流量应控制在0.028 kg/s以下;考虑到钻井成本,注采方式的选择以高注高采和中注高采为最佳;注入温度对热储开采影响较小,可选择40℃ ~ 80℃之间任意温度的地热尾水进行回灌,实现地热资源梯级利用。     

System performance of a deep borehole heat exchanger

Deep borehole heat exchanger (BHE) systems, installed in abandoned boreholes, have been operative in Switzerland for several years now. The operational conditions of the 2302 m deep BHE plant at Weggis have been monitored continuously since 1994. In the first operational phase, lasting from October 1994 to May 1996, the plant was severely underused, as shown by the high production temperatures (40 °C). This behaviour was investigated by a numerical model accounting for the heat transport in the rock matrix and along the different tubing systems, with special emphasis on the heat transfer in a multi-layer insulated central pipe. Lacking detailed logging data or undisturbed temperature profiles, an axis-symmetrical model had to be used, assuming uniform rock parameters. Sensitivity studies highlighted the effect of varying flow rate or operation/recovery cycle lengths and helped to develop a strategy that allowed us to make an accurate calculation of the long-term Weggis production history. The initial model assumptions, based on this detailed treatment of the tubing system, could explain the operational data. By means of slight model variations that account only for the minor effects of metallic sleeves, the long-term production temperature history of the Weggis plant could be accurately fitted. These findings were confirmed by a detailed analysis of the May 1996 data. Due to the low degree of utilization, only numerical sensitivity analyses were able to highlight the potential of the deep BHE plant at Weggis. The results indicate that the low utilisation of 40 kW during the first operational phase could be increased to over 200 kW. The specific yield of deep systems is much higher than in conventional shallow BHE systems. Our simulation procedure proves that the heat transfer in a deep BHE system is well understood.     

Flexibility valuation in the Brazilian power system: A real options approach

This paper presents a valuation study of operational flexibility in the complex Brazilian Power System. Thermopower plants represent operational flexibility for the national system operator provided they can be dispatched in dry periods to supply part of the load when reservoir waters are too low. Deficit costs can be avoided as a result.

We have adopted a real options approach to calculate the fair value of a financial subsidy to be paid to thermal generators for their availability to the system. A financial subsidy is mandatory for their economic feasibility and, therefore, for increased thermopower capacity in the Brazilian Power System. This is why this policy is currently being studied by Brazil's federal government.

In order to illustrate our modeling we have run the model for the southeast subsystem. We found a flexibility value of US$4.52 billion, which represents US$497/MW per year. This means that a 100 MW thermopower plant should receive US$49,700 for each year of its economic lifetime as a fair premium incentive to investments.