In geothermal reservoir systems, changes in pore pressure due to production (depletion), injection or temperature changes result in a displacement of the effective stresses acting on the rock matrix of the aquifer. To compensate for these intrinsic stress changes, the rock matrix is subjected to poroelastic deformation through changes in rock and pore volume. This in turn may induce changes in the effective pore network and thus in the hydraulic properties of the aquifer. Therefore, for the conception of precise reservoir models and for long-term simulations, stress sensitivity of porosity and permeability is required for parametrization. Stress sensitivity was measured in hydrostatic compression tests on 14 samples of rock cores stemming from two boreholes of the Upper Jurassic Malm aquifer of the Bavarian Molasse Basin. To account for the heterogeneity of this carbonate sequence, typical rock and facies types representing the productive zones within the thermal reservoir were used. Prior to hydrostatic investigations, the hydraulic (effective porosity, permeability) and geomechanical (rock strength, dynamic, and static moduli) parameters as well as the microstructure (pore and pore throat size) of each rock sample were studied for thorough sample characterization. Subsequently, the samples were tested in a triaxial test setup with effective stresses of up to 28 MPa (hydrostatic) to simulate in-situ stress conditions for depths up to 2000 m. It was shown that stress sensitivity of the porosity was comparably low, resulting in a relative reduction of 0.7–2.1% at maximum effective stress. In contrast, relative permeability losses were observed in the range of 17.3–56.7% compared to the initial permeability at low effective stresses. Stress sensitivity coefficients for porosity and permeability were derived for characterization of each sample and the different rock types. For the stress sensitivity of porosity, a negative correlation with rock strength and a positive correlation with initial porosity was observed. The stress sensitivity of permeability is probably controlled by more complex processes than that of porosity, where the latter is mainly controlled by the compressibility of the pore space. It may depend more on the compaction of precedented flow paths and the geometry of pores and pore throats controlling the connectivity within the rock matrix. In general, limestone samples showed a higher stress sensitivity than dolomitic limestone or dolostones, because dolomitization of the rock matrix may lead to an increasing stiffness of the rock. Furthermore, the stress sensitivity is related to the history of burial diagenesis, during which changes in the pore network (dissolution, precipitation, and replacement of minerals and cements) as well as compaction and microcrack formation may occur. This study, in addition to improving the quality of input parameters for hydraulic–mechanical modeling, shows that hydraulic properties in flow zones largely characterized by less stiff, porous limestones can deteriorate significantly with increasing effective stress.

In the field of the efficiency of very shallow geothermal energy systems, there is still a significant need for research activity. To ensure the proper exploitation of this energy resource, the decisive geophysical parameters of soil must be well-known. Within this study, thermal conductivity, as a fundamental property for evaluating the geothermal potential of very shallow geothermal systems, was analyzed and measured with a TK04 device. A dataset, consisting of various geophysical parameters (thermal conductivity, bulk density, water content, and porosity) determined for a large range of different textural soil classes, was collated. In a new approach, the geophysical properties were visualized covering the complete grain size range. The comparison between the measured and calculated thermal conductivity values enabled an investigation with respect to the validity of the different Kersten equations. In the course of this comparison, the influence of effective bulk density was taken into account. In conclusion, both Kersten formulas should be used as recommended and regular bulk density corresponded better to the reference dataset representing the outcomes of the TK04 laboratory measurement. Another objective was to visualize the relation of thermal conductivities within their corresponding textural classes and the validity of Kersten formulas for various bulk densities, depths, and soils. As a result, the accessibility to information for expedient recommendations about the feasibility of very shallow geothermal systems will be improved. Easy, accessible know-how of the fundamentals is important for a growing renewable energy sector where very shallow geothermal installations can also cover heating and cooling demands.

Abstract. The Franconian Alb of SE Germany is characterized by large-scale exposures of Jurassic shallow marine limestones and dolostones, which are frequently considered to be outcrop analogues for deep geothermal reservoir rocks in the North Alpine Foreland Basin farther south. However, the burial history of the Franconian Alb Jurassic strata is not well known as they were affected by emersion, leading to extensive erosion and karstification with only remnants of the original Cretaceous and Cenozoic cover rocks preserved. To estimate the original thicknesses of the post-Jurassic overburden we investigated the petrophysical properties and the thermal history of Lower and Middle Jurassic mudstones to constrain their burial history in the Franconian Alb area. We measured mudstone porosities, densities, and maturities of organic material and collected interval velocities from seismic refraction and logging data in shallow mudstone-rich strata. Mudstone porosities and P-wave velocities vertical to bedding were then related to a normal compaction trend that was calibrated on stratigraphic equivalent units in the North Alpine Foreland Basin. Our results suggest maximum burial depths of 900–1700 m, 300–1100 m of which is attributed to Cretaceous and younger sedimentary rocks overlying the Franconian Alb Jurassic units. Compared to previous considerations this implies a more widespread distribution and increased thicknesses of up to ∼900 m for Cretaceous and up to ∼200 m for Cenozoic units in SE Germany. Maximum overburden is critical to understand mechanical and diagenetical compaction of the dolostones and limestones of the Upper Jurassic of the Franconian Alb. The results of this study therefore help to better correlate the deep geothermal reservoir properties of the Upper Jurassic from outcrop to reservoir conditions below the North Alpine Foreland Basin. Here, the Upper Jurassic geothermal reservoir can be found at depths of up to 5000 m.

Abstract. The Franconian Alb of SE Germany is characterized by large-scale exposures of Jurassic shallow marine limestones and dolostones which are frequently considered as outcrop analogues for deep geothermal reservoir rocks in the North Alpine Foreland Basin farther south. However, the burial history of the Franconian Alb Jurassic strata is not well known as they were affected by emersion, leading to extensive erosion and karstification with only remnants of the original Cretaceous and Cenozoic cover rocks preserved. To estimate the original thicknesses of the post-Jurassic overburden we investigated the petrophysical properties and the thermal history of Lower and Middle Jurassic mudstones to constrain their burial history in the Franconian Alb area. We measured mudstone porosities, densities, and maturities of organic material and collected interval velocities from seismic refraction and logging data in shallow mudstone-rich strata. Mudstone porosities and P-wave velocities vertical to bedding were then related to a normal compaction trend that was calibrated on stratigraphic equivalent units in the North Alpine Foreland Basin. Our results suggest maximum burial depths of 900–1700 m of which 300–1100 m are attributed to Cretaceous and younger sedimentary rocks overlying the Franconian Alb Jurassic units. Compared to previous considerations this implies a more widespread distribution and increased thicknesses of up to ~900 m for Cretaceous and up to ~200 m for Cenozoic units in SE Germany. Maximum overburden is critical to understand mechanical and diagentical compaction of the dolostones and limestones of the Upper Jurassic of the Franconian Alb. The results of this study therefore help to better correlate the deep geothermal reservoir properties of the Upper Jurassic from outcrop to reservoir conditions below the North Alpine Foreland Basin. Here, the Upper Jurassic geothermal reservoir can be found at depths of up to 5000 m.

Tight carbonate rocks are important hydrocarbon and potential geothermal reservoirs, for example, in CO2-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 m2. 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.