Geothermal Energy最新文献

High-Temperature Aquifer Thermal Energy Storage (HT-ATES) can be used to reduce greenhouse gas emissions from heating. The thermal recovery efficiency is the main parameter indicating the performance of an HT-ATES system and it is influenced by multiple aquifer properties and storage characteristics. This study presents a method for estimating recovery efficiency through numerical modeling, data analysis, and curve fitting. This method shows the relation between the recovery efficiency and various storage conditions, such as aquifer properties and storage temperature. In addition, this research explores an analytical relationship between energetic efficiency and recovery efficiency and verifies that relationship with the generated data. The proposed method can be used for the purpose of initial screening to estimate the performance of an HT-ATES system and for efficiently using HT-ATES as a component in larger energy system models. This method uses the modified Rayleigh number in combination with aquifer thickness and injected volume and has a (R^2) of 85%. The analytical relation between energetic efficiency and recovery efficiency was shown to be accurate for all calculated energetic efficiency values above 60% and is less accurate with lower calculated energetic efficiency values.

The mechanism of heat transfer and accumulation in sedimentary basin geothermal anomaly belts is very complicated, resulting in a lack of in-depth understanding. The Bohai Bay Basin is the largest Mesozoic–Cenozoic rifted basin in China and taken as the research object in the present work, the relationship between geothermal anomaly characteristics and influencing factors, including the structural geometry of the bedrock bulge, the bedrock thermal conductivity, and the caprock thickness was quantitatively discussed. Based on a typical measured geothermal geologic profile in the central Bohai Bay Basin, the formation mechanism, controlling factors, and the guiding significance for geothermal resources exploration of geothermal anomaly belts were analyzed. The formation of a geothermal anomaly in the sedimentary basin was driven by the rapid transfer of uniform deep heat flow in bedrock with the high thermal conductivity and accumulation of that in caprock with the low thermal conductivity of shallow bedrock bulge belts. As a result, heat flow depletion in the deep bedrock and heat flow enrichment in the shallow caprock were formed, with the characteristics of a nearly symmetrical coupling variation pattern. The bedrock bulge amplitude and the bedrock thermal conductivity were the main factors controlling geothermal anomaly characteristics, and the sufficient caprock thickness was vital to the formation of the geothermal anomaly. Simulation results of theoretical models showed that the top surface temperature of bedrock geothermal reservoirs with a depth of 2 km can increase by up to 31.3 ℃, and positive anomaly amplitudes were up to 41.8%. Based on the difference of tectonic deformation patterns in the rifted basin and their determining role in a geothermal gradient of caprock, geothermal fields developing bedrock geothermal reservoirs in the Bohai Bay Basin were divided into seven types, including the high bulge type, the bulge type, the low bulge type, the sag type, the gentle slope type, the basin margin fault-uplift type, and the shallow burial type, with a decrease in turn of the caprock geothermal gradient from > 5.0 ℃/100 m to < 2.0 ℃/100 m. The high bulge-type and the bulge-type geothermal fields were favorable targets for the shallow-burial bedrock geothermal reservoir exploration. The low bulge-type and the sag-type geothermal fields were the main targets for exploring deep-burial bedrock geothermal reservoirs with high temperatures above 150℃.

This paper reports more than 500 newly measured values of the dynamic viscosity of pure and mixed sodium chloride (NaCl) and calcium chloride (CaCl₂) aqueous solutions of which about half are unique. The data were acquired with a commercial rolling ball viscometer at ambient pressure, temperatures between 293 and 353 K, concentrations up to 6.0 mol/kg (NaCl) and 5.3 mol/kg (CaCl₂) as well as five different mixing ratios for the ternary system. Compared to existing values found in the literature, the present data match within mostly 5% uncertainty. Potential sources of errors are thoroughly discussed. An empirical correlation was derived from the measured data, accurately reproducing the measured data with good precision. This permits the functional dependence of dynamic viscosity on concentration and temperature of pure and mixed NaCl and CaCl₂ solutions to be easily and reliably implemented in coupled thermal–hydraulic–mechanical–chemical (THMC) numerical models aiming to predict the evolution of geothermal reservoir dynamics. Moreover, numerical sensitivity analyses were conducted exemplarily for one deep geothermal (DG) as well as one aquifer thermal energy storage (ATES) system to constrain the effect of viscosity variations and/or uncertainty on the prediction of key reservoir and operational parameters. It is demonstrated that viscosity variations systematically affect the productivity and injectivity indices (both systems) and the required pumping power (ATES), though no effect was observed for the timing of thermal breakthrough (DG) and the temperature evolution at the warm well (ATES). Not least, the stated precision of the analytical data obtained in this study proves well sufficient for the quality of numerical reservoir simulations aiming at predictions in a geothermal context.

A three-dimensional hydraulic-thermal model has been developed for the island of Pantelleria in order to investigate the driving heat transfer mechanisms and the role of rocks and fluids thermal and hydraulic properties in controlling the temperature distribution. The model simulates transient conditions accounting for changes with time of the background steady-state thermal field over a period of 50,000 years, which coincides with a period of major morphological/structural stability on the island. The comparison between observed and computed temperatures was used to validate the model. Model outcomes confirm temperature up to 490 °C at 2 km depth in the central part of the Cinque Denti caldera. The proposed thermal modeling agrees with the trend of the thermal anomaly on the island, clearly evidencing how the center of the thermal anomaly is located within the young caldera, while the north-eastern part of the island is characterized by much lower temperatures than its central-southern part. It is shown that a coupled thermal–hydraulic 3D model can provide a clear perspective on surface and subsurface conditions that drive heat transfer, thus helping to understand the characteristics of geothermal systems in the area.

Exergoeconomics is a vital complementation of thermodynamic performance analysis. In this study, a comprehensive exergoeconomic analysis of a binary geothermal power plant in southwestern Anatolia is conducted to determine improvement potentials in the plant configuration. By utilization of cost allocation rules of three different exergoeconomic methods (Moran, Specific Exergy Costing (SPECO), Exergy Cost Theory) plant is analyzed in terms of exergetic costs and possible optimization areas. Levelized cost of electricity (LCOE) estimated by 3 different methods vary within a 3.6% range, between 7.81 c$/kWh and 8.1 c$/kWh. It has also been determined that 51.5% of LCOE is constituted by waste/residual costs. Components especially including a thermal phase change or energy conversion, whose exergoeconomic factors below 0.5 warrant investment and optimization for performance improvement despite their higher individual exergetic efficiencies reported in previous studies. This phenomenon highlights the importance of considering exergetic efficiency and exergoeconomic factors together as plant design parameters. By using advanced materials or by optimizing the temperature gradient between the geothermal brine and the working fluid, heat transfer efficiency can be enhanced in heat exchanger devices. Turbines generally have mechanical losses which can be enhanced by optimizing blade design, reducing friction, and enhancing the thermodynamic cycle (i.e., use re-heat stages or improve steam conditions).

The development of geothermal energy through enhanced geothermal systems (EGS) often encounters challenges such as fluid short-circuiting, water loss, and insufficient connectivity. This study presents a time-dependent seepage and heat exchange model for the formation–wellbore–fluid system during the heat extraction process. Taking the Fenton Hill HDR project as a case study, this paper investigates the influence of formation characteristics, wellbore design, and injected fluid properties on heat transfer efficiency. Furthermore, a multi-well EGS utilizing multiple fracturing horizontal wells (MFHW) is proposed, and its production temperature is compared with two types of double-well EGS. The findings reveal that within the horizontal segment of the double-well EGS, an optimal output of 3.4 MW can be achieved at an injection rate of 30 kg/s. Additionally, the extraction temperature shows a positive correlation with factors such as heat production and electrical power generation. In the MFHW project, optimizing heat production potential can be accomplished by increasing the number of perforation fractures, enhancing artificial fracture spacing, improving the perforation angle, extending the horizontal segment, reducing well diameter, and employing a longer vertical heat insulation pipe with lower thermal conductivity. Finally, a comparative analysis of various development models indicates that two-injection-one-production multi-well EGS model exhibits superior performance, with its heat production being twice as efficient as that of one-injection-one-production double-well EGS model.

Hydraulic stimulation of pre-existing fractures and faults plays a significant role in improving hydraulic conductivity of the fracture network around injection and production wells in deep geothermal reservoirs. In present work, a three-dimensional distinct element method (3DEC, Itasca) is used to simulate the year 2000 hydraulic stimulation of GPK2 well of Soultz-sous-Forêts geothermal reservoir, where several major hydraulic stimulations have been performed and are well documented. The field scale numerical model of the reservoir (about 6000 × 4500 × 4500 m³) includes an explicit description of the main fault (FZ4770), was developed to constrain the large-scale hydromechanical properties of the fault, in particular, its behavior in terms of non-linear elastic response related to fault aperture changes. The first phase of the stimulation is modelled as a constant flow rate of 30 ls⁻¹ of water injection into the center of a deformable fault at the depth of approximately 4.7 km. We observed that the fluid pressure front migration from the injection point along the fracture follows, under the in-situ stress condition and the moderate injection pressure, a pseudo-diffusion behavior as power-law function of time with a 0.5 exponent (√t) when the injection flow rate is constant. It is demonstrated that the dynamic evolution of aperture opening due to fluid injection into the fracture is responsible for the pressure propagation behavior, owing to a hydraulic aperture change rather than a fluid pressure diffusion process. This numerically observed propagation process is compatible with a high fault effective diffusivity of 13 m²/s as that observed in the field. In case of a linear increase of the injection flow rate, the pseudo-diffusion process disappears leading to a time dependent power-law migration of the pressure front with exponent of 0.75. The pressure propagation is shown to be strongly influenced by the injection scheme.

Deep geothermal closed-loops have recently gained attention because of their advantages over classical geothermal applications (e.g., less dependence on the geology, no risk of induced seismicity) and technological advantages (e.g., in the drilling process, use of alternative to water fluids). This paper deals with the repurposing of an existing well in Mol, Belgium, by numerically evaluating the closed-loop concept. Two numerical tools are used to predict the evolution of the temperature and the produced energy over a period of 20 years considering the vertical coaxial well and the complete geological morphology. Full-scale simulations are initially carried out to estimate the maximum capacity of the well and to highlight the need to control the output of the well by adjusting the inlet conditions. Simulations are then performed either to deliver a constant power or to cover irregular thermal energy demands of two buildings by applying in both cases three process control operations. Through controlling the inlet temperature, the injected flow rate or successively both, the production of excess energy, resulting from the overdesign of the existing wellbore for the specific application, is limited. The simulations showed that continuous adjustments to the injection temperature and/or flow rate are needed to restrict the rapid drop in outlet temperature and consequent thermal depletion of the rocks, caused by the highly transient nature of the diffusive heat transfer from the rocks to the wellbore, as well as to supply a specific heat demand, constant or irregular, over the long term. In fact, the combination of both controls could be the ideal strategy for supplying the demand at the highest COP.

Geothermal installations often face operational challenges related to scaling which can lead to loss in production, downtime, and an increase in operational costs. To accurately assess and minimize the risks associated with scaling, it is crucial to understand the interplay between geothermal brine composition, operating conditions, and pipe materials. The accuracy of scaling predictive models can be impacted by uncertainties in the brine composition, stemming from sub-optimal sampling of geothermal fluid, inhibitor addition, or measurement imprecision. These uncertainties can be further increased for fluid at extreme conditions especially high salinity and temperature. This paper describes a comprehensive method to determine operational control strategies to minimize the scaling considering brine composition uncertainties. The proposed modelling framework to demonstrate the optimization under uncertainty workflow consists of a multiphase flow solver coupled with a geochemistry model and an uncertainty quantification workflow to locally estimate the probability of precipitation potential, including its impact on the hydraulic efficiency of the geothermal plant by increasing the roughness and/or decreasing the diameter of the casings and pipelines. For plant operation optimization, a robust control problem is formulated with scenarios which are generated based on uncertainties in brine composition using an exhaustive search method. The modelling and optimization workflow was demonstrated in a geothermal case study dealing with barite and celestite scaling in a heat exchanger. The results showed the additional insights in the potential impact of brine composition uncertainties (aleatoric uncertainties) in scaling potential and precipitation location. Comparing the outcome of optimization problem for the deterministic and fluid composition uncertainties, a change of up to 2.5% in the temperature control settings was observed to achieve the optimal coefficient of performance.