Geothermal Energy最新文献

The presence of non-condensable gases (NCGs) in a geothermal fluid disrupts the vacuum process in the condenser, reducing turbine efficiency and decreasing the total power output of the geothermal power plant (GPP). Therefore, to optimize the thermodynamic efficiency of a GPP, NCGs should be removed using a gas removal system. Since there is a substantial lack of design and simulation software for a GPP including NCG removal alternatives, in this study, we aimed to model and develop a software-based interface to simulate mass and energy balance involving an NCG fraction in a single-flash GPP as well as examine the thermodynamic performance of the gas removal system, which is the most important step in the planning and designing phase of a GPP. This software was validated using outputs of Kamojang GPP Units 2, 3, and 4 located at Kamojang geothermal field, Garut, West Java, Indonesia. Units 2 and 3 use two ejectors which are installed in series, and Unit 4 utilizes a hybrid system (HS) that is mostly a combination of vacuum pumps and a steam jet ejector (SJE). Our results showed that Kamojang GPP Units 2 and 3 generate 55.295 MW of power with an absolute error of 0.53%, whereas Unit 4 generates 60.218 MW of power with a 0.36% absolute error concerning the field data. These results correspond with the expected minimum error; therefore, our model’s parameters are considered valid and can be used for simulation. We found that using the simulation, the total steam saved by the HS at Kamojang GPP Units 2 and 3 was 534 kW. Furthermore, the net power production was reduced by 1.6% for the HS and 2.03% for the SJE with every 1% increase in the NCG fraction. The power requirement of the vacuum pumps remained less than the power generated by the motive steam which the ejector requires to dispose of the same amount of NCG, indicating that Kamojang GPP Units 2 and 3 will be more efficient if an HS is used.

In the foreseeable future, the geothermal exploitation from hot dry rocks (HDR) using a horizontal well will bear potential. Thus, in-depth studies should be conducted on the selection of injection-production scheme (IPS) and working fluid, design of reinjection parameters, optimization of wellbore structure and materials, and analysis of geological settings. This paper proposed a fully coupled model to study the above scientific questions. For Model A, the working fluid was injected into the annulus and then flowed out of the thermal insulation pipe (TIP). Its temperature passes through two stages of temperature rise and two stages of temperature decline. But for model B, the working fluid was injected into the TIP and then flowed out of the annulus. Its temperature undergoes five stages, four stages of temperature rise and one stage of temperature decline. The results show that the Model A is the best IPS owing to its high outlet temperature, stable thermal recovery, and low fluid injection volume. In Model A, when the working fluid was supercritical carbon dioxide and the liquid injection volume was 135.73 m³/d, the heat recovery ratio (HRR) was as high as 85.40%, which was 17.85% higher than that of the Model B whose working medium was water, and its liquid injection volume was only 25% of that. Meanwhile, over ten years of continuous production, the outlet temperature decreased by 7.5 °C and 18.38 °C in the latter. The optimal working fluid has a low volume heat capacity and thermal conductivity for any IPS. Sensitivity studies showed that for the area that met the HDR standard, the effect of reinjection temperature on the outlet temperature can be ignored. As for Model A, HRR drops sharply by 6.74–9.32% when TIP goes from completely adiabatic to nonzero thermal conductivity. Meanwhile, the horizontal segment length of the TIP is shorter when Model A obtains the optimal outlet temperature compared with Model B. In addition, the correlation between the outlet temperature and different formations of thermophysical properties was seriously affected by the IPS and exploitation period, which was summarized in detail.

Compared with conventional fossil fuel sources, geothermal energy has several advantages. The produced geothermal energy is safe for the environment and suitable for meeting heating power needs. Because the hot water used in the geothermal process can be recycled and used to generate more steam, this energy is sustainable. Furthermore, the climate change does not affect geothermal power installations. This study suggests a combined power generation cycle replicating using the EES software that combines a single flash cycle with a trans-critical carbon dioxide cycle. The findings demonstrate that, in comparison to the BASIC single flash cycle, the design characteristics of the proposed system are greatly improved. The proposed strategy is then improved using the Nelder–Mead simplex method and Genetic Algorithm. The target parameter is exergy efficiency, and the three assumed variable parameters are separator pressure, steam turbine outlet pressure, and carbon dioxide turbine inlet pressure. The system’s exergy efficiency was 32.46% in the default operating mode, rising to 39.21% with the Genetic Algorithm and 36.16% with the Nelder–Mead simplex method. In the final step, the exergy destruction of different system components is calculated and analyzed.

Graphical Abstract

Renewable energy sources are key to achieve the transition toward clean energy system. Among them, the geothermal energy has a production whose effectiveness requires sufficient understanding of the temperature distribution and fluid circulation at depth, as well as of the lithological and petrophysical properties of the crust. The focus of this paper is twofold: first, we summarize the main advances in the development of new methodologies and numerical codes to characterize the properties of the thermal lithosphere in terms of its, temperature, density and composition; second, based on the compilation of available thermal modelling results, we present the depth of the thermal Lithosphere–Asthenosphere Boundary (LAB) of the Iberian Peninsula and the temperature distribution at crustal depths of 5, 10, and 20 km, in addition to at Moho level. At 5 km depth, the temperature is above 110 °C with local anomalies (> 130 °C) located in the Iberian Massif and Cenozoic volcanic provinces. A similar pattern is observed at 10 and 20 km depth, where temperatures are above 190 °C and 350 °C, respectively. At 20 km depth, anomalies above > 500 °C, delineate the SE and NE Cenozoic volcanic provinces. At Moho depths, temperature ranges from 450 to 800 °C with hot regions mainly located along the Iberian Massif and the SE and NE volcanic provinces. The compiled results do not show any lithospheric anomaly that could give rise to high temperatures at shallow depths, but they do show an acceptable exploitation potential at intermediate depths. With regard to the direct use of district and greenhouse heating and for industrial processes, the potential is great throughout the Peninsula, the main challenges being the availability of groundwater and drilling costs.

One key aspect in the energy transition is to use the deep geothermal energy stored in sedimentary basins as well as in igneous and metamorphic basement rocks. To estimate the variability of deep geothermal potentials across different geological domains as encountered in the Federal State of Hesse (Germany), it is necessary to understand the driving processes of fluid flow and heat transport affecting subsurface temperature variations. In this study, we quantify the stored energy in a set of geological units in the subsurface of Hesse with the method of “heat in place” (HIP, sensu Muffler and Cataldi in Geothermics 7:53–89, 1978)—HIP is one proxy for the geothermal potential of these units controlled by their temperature configuration as derived from a series of coupled 3D thermo-hydraulic numerical models. We show how conductive, advective and convective heat transport mechanisms influence the thermal field and thereby the HIP calculations. The heterogeneous geology of the subsurface of Hesse ranges from locally outcropping Paleozoic basement rocks to up to 3.8 km thick Cenozoic, porous sedimentary deposits in the tectonically active northern Upper Rhine Graben. The HIP was quantified for five sedimentary layers (Cenozoic, Muschelkalk, Buntsandstein, Zechstein, Rotliegend) as well as for the underlying basement. We present a set of maps allowing to identify geothermally prospective subregions of Hesse based on the laterally varying thermal energy stored within the units. HIP is predicted to be highest in the area of the northern Upper Rhine Graben in the Cenozoic unit with up to 700 GJ (text {m}^{-2}) and in the Rotliegend with up to 617 GJ (text {m}^{-2}). The calculations account for the variable thicknesses and temperatures of the layers, density and heat capacity of the solid and fluid parts of the rocks as well as porosity.

The Groß Schönebeck site in the North German Basin serves as research platform to study the geothermal potential of deeply buried Permian reservoir rocks and the technical feasibility of heat extraction. The structural setting of the site was investigated in more detail by a newly acquired 3D-seismic survey to improve the former conceptual model that was based on several old 2D seismic lines. The new data allow a revision of the geological interpretation, enabling the setup of a new reservoir model and providing base information for a possible further site development of Permo-Carboniferous targets. The 3D seismic allows for the first time a consistent geological interpretation and model parameterization of the well-studied geothermal site. Main reflector horizons and the corresponding stratigraphic units were mapped and the structural pattern of the subsurface presented in the 8 km × 8 km × 4 km large seismic volume. Attribute analysis revealed some fracture and fault patterns in the upper Zechstein and post-Permian units, while formerly hypothesized large offset faults are not present in the Rotliegend reservoir. However, a well-established graben-like structure at the top of the Zechstein succession is most likely related to broken anhydritic brittle intra-salt layers of some meter of thickness. Most reflectors above the salt show a rather undisturbed pattern. The main reservoir sandstone of the Dethlingen Formation (Rotliegend) was mapped and characterized. The base of the underlying Permo-Carboniferous volcanic rock sequence and hence its thickness could not be depicted reliably from the geophysical data. Based on the seismic data and the available reconnaissance drilling, logging, and laboratory data of the Groß Schönebeck research site, the thickness and distribution of the sedimentary Rotliegend (with emphasis of the sandy reservoir section) and of the volcanic rock sequence was modelled and stochastically parameterized with petrophysical properties guided by seismic facies pattern correlation, providing a more realistic reservoir description. Properties include total and effective porosity, permeability, bulk density, thermal conductivity, thermal diffusivity, and specific heat capacity. The data and interpretation constitute the basis for a better understanding of the thermo and hydromechanical processes at the site and for future measures. Further site development could include a deepening of one well to provide evidence on the volcanic rock sequence and consider deviated wells into favourable zones and the design of a fracture-dominated utilization approach.

The successful exploitation of geothermal reservoirs relies upon the understanding of fluid circulation in the subsurface. However, large-scale fluid flow modelling often assumes that the permeability of the layers of rock within the model are isotropic. We present here a laboratory study in which we assessed the permeability anisotropy of seven Buntsandstein sandstone cores taken from the geothermal reservoir at Soultz-sous-Forêts (France) in the Upper Rhine Graben. The porosity and permeability of our samples, cored parallel and perpendicular to bedding, ranged from 5.2 to 16.3% and from 2.48 × 10⁻¹⁸ to 7.66 × 10⁻¹⁴ m², respectively. Our data show that permeability anisotropy can be up to four orders of magnitude in sandstones from the Buntsandstein, and that permeability anisotropy increases as a function of increasing porosity. Quantitative microstructural analysis combined with permeability modelling shows that the permeability anisotropy is the result of fine-grained and low-permeability laminations that are parallel or sub-parallel to bedding. We suggest, based on our data, that permeability anisotropy should be considered in future fluid flow modelling at geothermal sites within the Upper Rhine Graben.

Tight carbonate rocks are important hydrocarbon and potential geothermal reservoirs, for example, in CO₂-Enhanced Geothermal Systems. We report a study of outcrop samples of tectonically undeformed tight carbonates from the upper Jurassic “Malm ß” formation in Southern Germany near the town of Simmelsdorf (38 km NE of Nuremberg) to understand bulk petrophysical properties in relation to microstructure and to compare models for permeability prediction in these samples. We applied Archimedes isopropanol immersion, Helium pycnometry, mercury injection, gamma density core logging, and gas permeability measurements, combined with microstructural investigations and liquid metal injection (LMI-BIB-SEM). In addition, ultrasonic velocity was measured to allow geomechanical comparison of stratigraphically equivalent rocks in the South German Molasse Basin (SGMB). Results show only small variations, showing that the formation is rather homogeneous with bulk porosities below 5% and argon permeabilities around 1.4E−17 m². The presence of stylolites in some of the samples has neither a significant effect on porosity nor permeability. Pores are of submicron size with pore throats around 10 nm and connected as shown by Mercury injection and Liquid Metal injection. Samples have high dynamic Young’s Modulus of 73 ± 5 GPa as expected for lithified and diagenetically overmature limestones. Moreover, no trends in properties were observable toward the faults at meter scale, suggesting that faulting was post-diagenetic and that the matrix permeabilities were too low for intensive post-diagenetic fluid–rock interaction. Petrophysical properties are very close to those measured in the SGMB, illustrating the widespread homogeneity of these rocks and justifying the quarry as a reasonable reservoir analog. Permeability prediction models, such as the percolation theory-based Katz-Thompson Model, Poiseuille-based models, like the Winland, the Dastidar, the capillary tube, and the Kozeny-Carman Models, as well as several empirical models, namely, the Bohnsack, the Saki, and the GPPT Models, were applied. It is shown that the capillary tube Model and the Saki Model are best suited for permeability predictions from BIB-SEM and mercury injection capillary pressure results, respectively, providing a method to estimate permeability in the subsurface from drill cuttings. Matrix permeability is primarily controlled by the pore (throat) diameters rather than by the effective porosity.

Indonesia has high geothermal potential comprising 40% of the world’s potential geothermal energy, volcanic and non-volcanic systems. Volcanic systems have witnessed more exploration activities for geothermal resources compared to non-volcanic systems. A high potential non-volcanic system in Indonesia is located in the northern part of Konawe, Southeast Sulawesi. Previous research had identified surface temperature anomaly (high temperature) and some surface manifestations for this area, specifically in the northeast part of Wawolesea. However, the source of surface manifestations and permeable zones as an implication of a good reservoir are still unknown. Therefore, this research aims to investigate the permeable zones and geothermal potential in the non-volcanic geothermal system of north Wawolesea by applying lineaments analysis and the fault fracture density (FFD) method. A total of 1694 major and minor lineaments were manually delineated using ArcGIS based on Digital Elevation Model Nasional (DEMNAS). FFD map and rose diagrams displayed the orientation of all lineaments and structures with the major lineaments trending NNE–SSW, whereas the minor lineaments showed irregular distribution and orientation. Field measurements also show the same azimuth orientation for the mapped fractures. Five zones were characterized by high FFD values (2.81–4.54 km/km²). One of the extensively fractured zones (Zone C) is located between Meluhu and Lembo, covering an area of around 19.39 km². This area is interpreted to be highly permeable and suggestive of a recharge area that contributes to surface manifestation in the Wawolesea. Therefore, the area between Meluhu and Lembo in the northern part of Konawe shows high geothermal potential due to its planar morphology and high FFD values. This study allows an improved understanding of how fracture geometry, distribution and density control the permeability in geothermal reservoirs.