Controls on the deep thermal field: implications from 3-D numerical simulations for the geothermal research site Groß Schönebeck

The deep thermal field in sedimentary basins can be affected by convection, conduction or both resulting from the structural inventory, physical properties of geological layers and physical processes taking place therein. For geothermal energy extraction, the controlling factors of the deep thermal field need to be understood to delineate favorable drill sites and exploitation compartments. We use geologically based 3-D finite element simulations to figure out the geologic controls on the thermal field of the geothermal research site Groß Schönebeck located in the E part of the North German Basin. Its target reservoir consists of Permian Rotliegend clastics that compose the lower part of a succession of Late Carboniferous to Cenozoic sediments, subdivided into several aquifers and aquicludes. The sedimentary succession includes a layer of mobilized Upper Permian Zechstein salt which plays a special role for the thermal field due to its high thermal conductivity. Furthermore, the salt is impermeable and due to its rheology decouples the fault systems in the suprasalt units from subsalt layers. Conductive and coupled fluid and heat transport simulations are carried out to assess the relative impact of different heat transfer mechanisms on the temperature distribution. The measured temperatures in 7 wells are used for model validation and show a better fit with models considering fluid and heat transport than with a purely conductive model. Our results suggest that advective and convective heat transport are important heat transfer processes in the suprasalt sediments. In contrast, thermal conduction mainly controls the subsalt layers. With a third simulation, we investigate the influence of a major permeable and of three impermeable faults dissecting the subsalt target reservoir and compare the results to the coupled model where no faults are integrated. The permeable fault may have a local, strong impact on the thermal, pressure and velocity fields whereas the impermeable faults only cause deviations of the pressure field.     

The influence of model structure on groundwater recharge rates in climate-change impact studies

Numerous modeling approaches are available to provide insight into the relationship between climate change and groundwater recharge. However, several aspects of how hydrological model choice and structure affect recharge predictions have not been fully explored, unlike the well-established variability of climate model chains—combination of global climate models (GCM) and regional climate models (RCM). Furthermore, the influence on predictions related to subsoil parameterization and the variability of observation data employed during calibration remain unclear. This paper compares and quantifies these different sources of uncertainty in a systematic way. The described numerical experiment is based on a heterogeneous two-dimensional reference model. Four simpler models were calibrated against the output of the reference model, and recharge predictions of both reference and simpler models were compared to evaluate the effect of model structure on climate-change impact studies. The results highlight that model simplification leads to different recharge rates under climate change, especially under extreme conditions, although the different models performed similarly under historical climate conditions. Extreme weather conditions lead to model bias in the predictions and therefore must be considered. Consequently, the chosen calibration strategy is important and, if possible, the calibration data set should include climatic extremes in order to minimise model bias introduced by the calibration. The results strongly suggest that ensembles of climate projections should be coupled with ensembles of hydrogeological models to produce credible predictions of future recharge and with the associated uncertainties.     

The stress regime in a Rotliegend reservoir of the Northeast German Basin

In-situ stresses have significant impact, either positive or negative, on the short and long term behaviour of fractured reservoirs. The knowledge of the stress conditions are therefore important for planning and utilization of man-made geothermal reservoirs. The geothermal field Groß Schönebeck (40 km north of Berlin/Germany) belongs to the key sites in the northeastern German Basin. We present a stress state determination for this Lower Permian (Rotliegend) reservoir by an integrated approach of 3D structural modelling, 3D fault mapping, stress ratio definition based on frictional constraints, and slip-tendency analysis. The results indicate stress ratios of the minimum horizontal stress S _hmin being equal or increasing 0.55 times the amount of the vertical stress S _V (S _hmin ≥ 0.55S _V ) and of the maximum horizontal stress S _Hmax ≤ 0.78–1.00S _V in stress regimes from normal to strike slip faulting. Thus, acting stresses in the 4,100-m deep reservoir are S _V  = 100 MPa, S _hmin = 55 MPa and S _Hmax = 78?100 MPa. Values from hydraulic fracturing support these results. Various fault sets of the reservoir are characterized in terms of their potential to conduct geothermal fluids based on their slip and dilatation tendency. This combined approach can be adopted to any other geothermal site investigation.     

Estimating surface runoff and groundwater recharge in an urban catchment using a water balance approach

Hydrogeology Journal - Land-use changes often have significant impact on the water cycle, including changing groundwater/surface-water interactions, modifying groundwater recharge zones, and...     

Modelling of fractured carbonate reservoirs: outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany

Fluid flow in low-permeable carbonate rocks depends on the density of fractures, their interconnectivity and on the formation of fault damage zones. The present-day stress field influences the aperture hence the transmissivity of fractures whereas paleostress fields are responsible for the formation of faults and fractures. In low-permeable reservoir rocks, fault zones belong to the major targets. Before drilling, an estimate for reservoir productivity of wells drilled into the damage zone of faults is therefore required. Due to limitations in available data, a characterization of such reservoirs usually relies on the use of numerical techniques. The requirements of these mathematical models encompass a full integration of the actual fault geometry, comprising the dimension of the fault damage zone and of the fault core, and the individual population with properties of fault zones in the hanging and foot wall and the host rock. The paper presents both the technical approach to develop such a model and the property definition of heterogeneous fault zones and host rock with respect to the current stress field. The case study describes a deep geothermal reservoir in the western central Molasse Basin in southern Bavaria, Germany. Results from numerical simulations indicate that the well productivity can be enhanced along compressional fault zones if the interconnectivity of fractures is lateral caused by crossing synthetic and antithetic fractures. The model allows a deeper understanding of production tests and reservoir properties of faulted rocks.     

An integrative geothermal resource assessment study for the siliciclastic Granite Wash Unit,northwestern Alberta (Canada)

This study explores the siliciclastic Granite Wash Unit in northwestern Alberta as a potential geothermal reservoir. The approach covers regional 3D structural geological modelling of a 90 km × 70 km area based on well log and legacy 2D seismic data. The fault strike was interpreted from lineaments, which were identified with the refined trend surface analysis method. The stress state of the Granite Wash reservoirs was determined by an integrated approach of 3D fault modelling, stress ratio definition based on frictional constraints, and slip tendency analysis. The results show that the best site for a geothermal application is located in the southwestern study area, where the highest temperatures (above 70 °C) coincide with the largest thickness (above 20 m) and zones of elevated porosity and permeability. The integrated stress analysis indicates an in situ stress regime from normal to strike-slip faulting maintaining a non-critically stressed reservoir or faults therein, assuming a friction coefficient of 0.7. The granite wash reservoirs could be used for heating of greenhouses, domestic warm water provision and district heating.