2015, Willeit, M., Ganopolski, A.

Permafrost influences a number of processes which are relevant for local and global climate. For example, it is well known that permafrost plays an important role in global carbon and methane cycles. Less is known about the interaction between permafrost and ice sheets. In this study a permafrost module is included in the Earth system model CLIMBER-2, and the coupled Northern Hemisphere (NH) permafrost–ice-sheet evolution over the last glacial cycle is explored. The model performs generally well at reproducing present-day permafrost extent and thickness. Modeled permafrost thickness is sensitive to the values of ground porosity, thermal conductivity and geothermal heat flux. Permafrost extent at the Last Glacial Maximum (LGM) agrees well with reconstructions and previous modeling estimates. Present-day permafrost thickness is far from equilibrium over deep permafrost regions. Over central Siberia and the Arctic Archipelago permafrost is presently up to 200–500 m thicker than it would be at equilibrium. In these areas, present-day permafrost depth strongly depends on the past climate history and simulations indicate that deep permafrost has a memory of surface temperature variations going back to at least 800 ka. Over the last glacial cycle permafrost has a relatively modest impact on simulated NH ice sheet volume except at LGM, when including permafrost increases ice volume by about 15 m sea level equivalent in our model. This is explained by a delayed melting of the ice base from below by the geothermal heat flux when the ice sheet sits on a porous sediment layer and permafrost has to be melted first. Permafrost affects ice sheet dynamics only when ice extends over areas covered by thick sediments, which is the case at LGM.

2019, Bauer, Johanna F., Krumbholz, Michael, Luijendijk, Elco, Tanner, David C.

Geothermal energy is an important and sustainable resource that has more potential than is currently utilized. Whether or not a deep geothermal resource can be exploited, mostly depends on, besides temperature, the utilizable reservoir volume over time, which in turn largely depends on petrophysical parameters. We show, using over 1000 (n=1027) 4-D finite-element models of a simple geothermal doublet, that the lifetime of a reservoir is a complex function of its geological parameters, their heterogeneity, and the background hydraulic gradient (BHG). In our models, we test the effects of porosity, permeability, and BHG in an isotropic medium. Furthermore, we simulate the effect of permeability contrast and anisotropy induced by layering, fractures, and a fault. We quantify the lifetime of the reservoir by measuring the time to thermal breakthrough, i.e. how many years pass before the temperature of the produced fluid falls below the 100 ∘C threshold. The results of our sensitivity study attest to the positive effect of high porosity; however, high permeability and BHG can combine to outperform the former. Particular configurations of all the parameters can cause either early thermal breakthrough or extreme longevity of the reservoir. For example, the presence of high-permeability fractures, e.g. in a fault damage zone, can provide initially high yields, but it channels fluid flow and therefore dramatically restricts the exploitable reservoir volume. We demonstrate that the magnitude and orientation of the BHG, provided permeability is sufficiently high, are the prime parameters that affect the lifetime of a reservoir. Our numerical experiments show also that BHGs (low and high) can be outperformed by comparatively small variations in permeability contrast (103) and fracture-induced permeability anisotropy (101) that thus strongly affect the performance of geothermal reservoirs.

2013, Schumacher, S., Schulz, R.

In Germany, many hydro-geothermal plants have been constructed in recent years, primarily in the region of Munich. As the host formation here mainly consists of carbonates, nearly all recently drilled wells have been acidized in order to improve the well yield. In this study, the effectiveness of these acid treatments is analyzed with respect to the amount of acid used and the number of acid treatments carried out per well. The results show that the first acid treatment has the largest effect, while subsequent acidizing improves the well only marginally. Data also indicate that continued acidizing can lead to degradation of the well. These findings may not only be important for geothermal installations in Germany but also for projects, for example, in Austria, France or China where geothermal energy is produced from carbonate formations as well.

2019, Mraz, Elena, Wolfgramm, Markus, Moeck, Inga, Thuro, Kurosch

The recent interest on environmentally friendly energy resources has increased the economic interest on the Upper Jurassic carbonate rocks in the North Alpine Foreland Basin, which serves as a hydrogeothermal reservoir. An economic reservoir use by geothermal fluid extraction and injection requires a decent understanding of porosity–permeability evolution of the deep laying Upper Jurassic strata at depths greater than 2000 m. The analysis of paleofluids caught in cements of the rock mass helps to determine the postdepositional reservoir evolution and fluid migration. Therefore, the high- and low-permeability areas of the Upper Jurassic in the North Alpine Foreland Basin referred to as Molasse Basin were analyzed by means of encountered postdepositional cements to determine the reservoir evolution. The cements were sampled at different hydrocarbon and geothermal wells, as well as at outcrops in the Franconian and Swabian Alb. To determine the composition and temperature of the paleofluids, fluid inclusions and cements of the Upper Jurassic carbonate rocks were analyzed by microthermometry and stable isotope measurements. Since drill cuttings are a rather available sample material compared to drill cores, a new microthermometry measurement method was achieved for the around 1 mm drill cuttings. Salinity and formation temperature of paleofluids in fluid inclusions and isotope data are consistent with previous studies and reveal a 5-stage evolution: the main cementation phases are composed of (I) the early diagenesis in limestones (200-400 m, 40-50°C), (II) early diagenetic dolomitization, and (III) burial dolomitization (1-2 km, II: 40-90°C; III: 70-100°C; 40 g/L NaCl equiv.), and (IV) late burial calcification (IIIa: 110-140°C, IIIb: 140-200°C) linked to tectonic features in the Molasse Basin. In the outcrop samples, a subsequent (V) cementation phase was determined controlled by karstification. In the southwest, an increase in salinity of the fluid inclusions in vein calcites, above the salinity of the Jurassic seawater, highlights the influence of basin fluids (diagenetic, evaporitic). In the other eastern wells, vein calcites have precipitated from a low saline fluid of around 10-20 g/L NaCl equiv. The low salinity and the isotope values support the theory of a continuous influence of descending meteoric fluids. Consequently, the Upper Jurassic seawater has been diluted by a meteoric fluid to a low saline fluid (<1 g/L), especially in areas with high permeability. Here, we show how a better understanding of cementation trajectory at depth can help to generate a better understanding of geothermal usability in deep carbonate reservoirs.