Elisa Heim

Ground-source heat pumps (GSHPs) coupled with borehole heat exchangers (BHEs) are energy-efficient technologies for heating and cooling buildings. However, these systems often fail to operate at their full potential due to discrepancies between the assumptions made during the design phase and the actual conditions during operation. To enhance overall GSHP performance, it is crucial to collect and analyze long-term monitoring data from operating BHE fields. To our knowledge, no long-term, high-resolution dataset of double U-tube BHEs is currently publicly available. Additionally, most studies typically monitor only the inlet and outlet of the entire ground heat exchanger rather than individual BHEs, hindering detailed performance analysis. With this data descriptor, we present a 6-year dataset from a BHE field comprising 40 BHEs, each with sensors for volume flow and inlet/outlet temperatures, recorded every 30 seconds. We believe this dataset will enhance understanding of individual BHE performance, provide validation for BHE models, and thus support better GSHP design and operation.

Model predictive control (MPC) is a promising approach for optimizing the performance of borehole heat exchangers (BHEs) in ground-source heat pump systems. The central element of MPC is the forward model that predicts the thermal dynamics in the ground. In this work, we validate the prediction accuracy of four BHE modeling approaches against real-world measurement data across various operational events and timescales. We simulate the fluid temperature leaving a BHE using a fully discretized 3-D numerical model, a resistance–capacitance model, a g-function model, and a hybrid model. The simulated temperatures are compared to measured temperatures using three validation metrics that quantify temperature offset, noise, and accuracy. The main reason for a mismatch between measured and modeled temperatures is a temperature offset of the simulated temperature. To remove this effect, the models were calibrated for their most sensitive parameter, the ground temperature, and their prediction accuracy over 4 years was evaluated. Thereby, model calibration seems to be a viable solution to account for an unknown load history. The results show that the resistance–capacitance model provides decent predictions in the short term and the g-function model in the long term. However, both models are strongly dependent on accurate calibration. The hybrid model provides the most accurate short and long-term predictions and is less dependent on calibration. Still, its integration into optimization syntax poses challenges compared to the other models. Although not yet applied in model predictive control, the hybrid model stands out as a promising choice for optimizing BHE field operations across various timescales.

Quantitative uncertainty analysis, 2-D and 3-D modeling of the subsurface, as well as their visualization form the basis for decision making in exploration, nuclear waste storage and seismic hazard assessment. Methods such as cross-section balancing are well established and yield plausible kinematic scenarios. However, they are based on geological data with errors and subject to human biases. Additionally, kinematic models do not provide a quantitative measure of the uncertainty of structures at depth. New 3-D modeling approaches have emerged that use computational interpolation, which are less dependent on human biases. Probabilistic extensions enable the quantification of uncertainties for the modeled structures. However, these approaches do not provide information on the time evolution of structures. Here, we compare classical cross-section balancing (2-D, kinematic modeling) with 3-D computational modeling to pave the way towards a solution that can bridge between these approaches. We show the strengths and weaknesses of both approaches, highlighting areas where probabilistic modeling can possibly add quantitative structural uncertainty information to improve section balancing. On the other hand, we show where probabilistic modeling still falls short of being able to cover the observed geometric complexities. We ultimately discuss how a workflow that iteratively combines results of the approaches can improve structural and kinematic constraints. As an example, we use the fold-and-thrust belt of the northern Alpine Foreland, the so-called Subalpine Molasse, focusing on the Hausham Syncline (Bavaria) and adjacent areas. We take advantage of the fact that here the stratigraphy as well as the tectonic history are well constrained. We show that shortening within the syncline progressively increases from west to east, independent from structural uncertainties. Two equally viable models can explain this. First, strain in the west is accommodated underneath the syncline in a triangle zone that progressively tapers out, or second, the strain difference is accommodated in more internal units. This highlights the importance of introducing uncertainty modeling also in kinematic restorations, as it enables identifying key regions, where different hypotheses can be tested.

Over 50% of the energy consumption in Europe is used for heating and cooling, of which the majority still comes from fossil fuels. For the heating and cooling supply of buildings, heat pumps using borehole heat exchangers (BHE) are an efficient and sustainable alternative. Design and performance of BHE fields strongly depend on the thermal properties of the subsurface. For the dimensioning of BHE fields, the subsurface properties are assessed using geothermal potential maps published by federal states or authorities.Such maps are created based on available geological and hydrogeological subsurface data. From these data, the stratigraphic sequence a BHE would intersect is inferred and thermal properties are assigned to the rock sequence. The mean thermal conductivity of the intersected sequence has a high uncertainty that is not quantified in most geothermal potential estimates. Moreover, geothermal potential maps are often not comparable because the used calculation methods differ in terms of complexity and calculated quantity. In Germany, for example, almost each of the federal states published its own geothermal potential map showing different quantities.For these reasons, we propose an approach to estimate the subsurface thermal conductivity together with its uncertainty. To enable applicability to datasets of different regions, e.g., all German federal states, the method takes only publicly available and standardized hydrogeological and geological data as input. Thermal conductivity estimation takes place on1) a large scale, to obtain a general estimate of thermal conductivity, and2) on a fine scale, to update the large-scale estimate in smaller target areas.For both estimations, we consider uncertainty by assigning not only the mean but also the minimum and maximum thermal conductivities of each rock type to the intersected sequence. Data handling and initiation of calculations is done via a geoportal. It connects the calculation method to geodata servers, instead of storing data locally, thus ensuring data up-to-dateness. This way, the approach does not only estimate the thermal conductivity, but also automatizes the cumbersome step of geological data gathering for all types of subsurface investigations.We demonstrate the approach using data from the federal state of Hamburg and compare our results to a thermal conductivity map published by the city of Hamburg.  The partly considerable differences between the two estimates illustrate the uncertainty of the geothermal potential, and thus the need to communicate it.