ALBERT

Description: The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-East Atlantic, the Mediterranean, and connected Seas (NEAM). In this online data product, the hazard results are provided by hazard curves calculated at 2,343 Points of Interest (POI), distributed in the North-East Atlantic (1,076 POIs), the Mediterranean Sea (1,130 POIs), and the Black Sea (137 POIs) at an average spacing of ~20 km. For each POI, hazard curves are given for the mean, 2nd, 16th, 50th, 84th, and 98th percentiles. Maps derived from hazard curves are Probability maps for Maximum Inundation Heights (MIH) of 1, 2, 5, 10, 20 meters; Hazard maps for Average Return Periods (ARP) of 500, 1,000, 2,500, 5,000, 10,000 years. For each map, precalculated displays are provided for the mean, the 16th percentile, and the 84th percentile. All data are also made accessible through an interactive web mapper and through Open Geospatial Consortium standard protocols. The model was prepared in the framework of the European Project TSUMAPS-NEAM (http://www.tsumaps-neam.eu/) funded by the mechanism of the European Civil Protection and Humanitarian Aid Operations (grant no. ECHO/SUB/2015/718568/PREV26).

Description: PyGMT is a library for processing geospatial and geophysical data and making publication quality maps and figures. It provides a Pythonic interface for the Generic Mapping Tools (GMT), a command-line program widely used in the Earth Sciences.

Description: Skarn deposits are found on every continents and were formed at different times from Precambrian to Tertiary. Typically, the formation of a skarn is induced by a granitic intrusion in carbonates-rich sedimentary rocks. During contact metamorphism, fluids derived from the granite interact with the sedimentary host rocks, which results in the formation of calc-silicate minerals at the expense of carbonates. Those newly formed minerals generally develop in a metamorphic zoned aureole with garnet in the proximal and pyroxene in the distal zone. Ore elements contained in magmatic fluids are precipitated due to the change in fluid composition. The temperature decrease of the entire system, due to the cooling of magmatic fluids and the entering of meteoric water, allows retrogression of some prograde minerals. The Hämmerlein skarn deposit has a multi-stage history with a skarn formation during regional metamorphism and a retrogression of primary skarn minerals during the granitic intrusion. Tin was mobilized during both events. The 340 Ma old tin-bearing skarn minerals show that tin was present in sediments before the granite intrusion, and that the first Sn enrichment occurred during the skarn formation by regional metamorphism fluids. In a second step at ca. 320 Ma, tin-bearing fluids were produced with the intrusion of the Eibenstock granite. Tin, which has been added by the granite and remobilized from skarn calc-silicates, precipitated as cassiterite. Compared to clay or marl, the skarn is enriched in Sn, W, In, Zn, and Cu. These metals have been supplied during both regional metamorphism and granite emplacement. In addition, the several isotopic and chemical data of skarn samples show that the granite selectively added elements such as Sn, and that there was no visible granitic contribution to the sedimentary signature of the skarn The example of Hämmerlein shows that it is possible to form a tin-rich skarn without associated granite when tin has already been transported from tin-bearing sediments during regional metamorphism by aqueous metamorphic fluids. These skarns are economically not interesting if tin is only contained in the skarn minerals. Later alteration of the skarn (the heat and fluid source is not necessarily a granite), however, can lead to the formation of secondary cassiterite (SnO2), with which the skarn can become economically highly interesting.

Description: Skarn-Lagerstätten befinden sich auf allen Kontinenten und wurden zu unterschiedlichen Zeiten vom Präkambrium bis zum Tertiär gebildet. Typischerweise wird die Bildung eines Skarns durch die Intrusion eines Granits in karbonatreiche Sedimentgesteine induziert. Während der Kontaktmetamorphose reagieren die Fluide aus dem Granit mit dem sedimentären Wirtgestein, was zur Bildung von Kalksilikaten auf Kosten von Karbonaten führt. Diese neu gebildeten Minerale entwickeln sich im Allgemeinen in einer metamorph zonierten Aureole mit Granat im proximalen und Pyroxen im distalen Bereich. Erzelemente die in magmatischen Fluiden enthalten sind werden aufgrund der veränderten Fluidzusammensetzung ausgefällt. Die Temperaturabsenkung des gesamten Systems, hervorgerufen durch die Abkühlung von magmatischen Fluiden sowie durch das Eindringen meteorischen Wassers, führen zu teilweisen oder vollständigen Umwandlung prograder Minerale. Die Skarn-Lagerstätte Hämmerlein hat eine mehrstufige Geschichte mit Skarnsbildung während der regionalen Metamorphose und Retrogression der primären Skarn-Minerale während der Intrusion von Graniten. Zinn wurde während beiden Ereignissen mobilisiert. Die 340 Ma alten zinnhaltigen Skarnminerale zeigen, dass Zinn in Sedimenten bereits vor dem Graniteintrag vorhanden war, und dass die erste Sn-Anreicherung während der Bildung des Skarns durch Fluide der Regionalmetamorphose stattfand. In einem zweiten Schritt um 320 Ma wurden Zinn-haltige Fluide durch die Intrusion des Eibenstockgranits freigesetzt. Diese Fluide überprägten den Skarn. Das freisetzen und das neu zugefügte Zinn ist in Kassiterit gebunden und führten dem System zusätzliches Zinn zu, wobei Zinn aus den Skarn-Kalksilikaten remobilisiert wurde. Im Vergleich zu Tonstein oder Mergel sind die Skarn mit Sn, W, In, Zn, und Cu angereichet. Diese Metalle sind während der Regionalmetamorphose und der Granitplatznahme zu unterschiedlichen Teilen zugeführt worden. Darüber hinaus zeigen die verschiedenen isotopen und chemischen Daten der Skarn-Proben, dass der Granit selektiv einige Elemente wie Sn hinzugefügt, und dass es keinen sichtbar granitischen Beitrag zur sedimentären Signatur des Skarns gab. Das Beispiel Hämmerlein zeigt, dass es möglich ist einen zinnreichen Skarn ohne zugehörigen Granit zu bilden, wenn Zinn von zinnhaltigen Sedimenten während einer Regionalmetamorphose mit wässrigen metamorphen Fluiden transportiert worden ist. Diese Skarne sind wirtschaftlich uninterssant wenn das Zinn nur in den Skarn-Mineralen enthalten ist. Spätere Umwandlung des Skarns (die Quelle der Wärme und Fluiden ist nicht unbedingt ein Granit) kann jedoch zur Bildung von sekundärem Kassiterite (SnO2) führen, womit der Skarn plötzlich wirtschaftlich hoch interessant sein kann.

Description: Earthquake swarms are characterized by large numbers of events occurring in a short period of time within a confined source volume and without significant mainshock aftershock pattern as opposed to tectonic sequences. Intraplate swarms in the absence of active volcanism usually occur in continental rifts as for example in the Eger Rift zone in North West Bohemia, Czech Republic. A common hypothesis links event triggering to pressurized fluids. However, the exact causal chain is often poorly understood since the underlying geotectonic processes are slow compared to tectonic sequences. The high event rate during active periods challenges standard seismological routines as these are often designed for single events and therefore costly in terms of human resources when working with phase picks or computationally costly when exploiting full waveforms. This methodological thesis develops new approaches to analyze earthquake swarm seismicity as well as the underlying seismogenic volume. It focuses on the region of North West (NW) Bohemia, a well studied, well monitored earthquake swarm region. In this work I develop and test an innovative approach to detect and locate earthquakes using deep convolutional neural networks. This technology offers great potential as it allows to efficiently process large amounts of data which becomes increasingly important given that seismological data storage grows at increasing pace. The proposed deep neural network trained on NW Bohemian earthquake swarm records is able to locate 1000 events in less than 1 second using full waveforms while approaching precision of double difference relocated catalogs. A further technological novelty is that the trained filters of the deep neural network’s first layer can be repurposed to function as a pattern matching event detector without additional training on noise datasets. For further methodological development and benchmarking, I present a new toolbox to generate realistic earthquake cluster catalogs as well as synthetic full waveforms of those clusters in an automated fashion. The input is parameterized using constraints on source volume geometry, nucleation and frequency-magnitude relations. It harnesses recorded noise to produce highly realistic synthetic data for benchmarking and development. This tool is used to study and assess detection performance in terms of magnitude of completeness Mc of a full waveform detector applied to synthetic data of a hydrofracturing experiment at the Wysin site, Poland. Finally, I present and demonstrate a novel approach to overcome the masking effects of wave propagation between earthquake and stations and to determine source volume attenuation directly in the source volume where clustered earthquakes occur. The new event couple spectral ratio approach exploits high frequency spectral slopes of two events sharing the greater part of their rays. Synthetic tests based on the toolbox mentioned before show that this method is able to infer seismic wave attenuation within the source volume at high spatial resolution. Furthermore, it is independent from the distance towards a station as well as the complexity of the attenuation and velocity structure outside of the source volume of swarms. The application to recordings of the NW Bohemian earthquake swarm shows increased P phase attenuation within the source volume (Qp 〈 100) based on results at a station located close to the village Luby (LBC). The recordings of a station located in epicentral proximity, close to Nový Kostel (NKC), show a relatively high complexity indicating that waves arriving at that station experience more scattering than signals recorded at other stations. The high level of complexity destabilizes the inversion. Therefore, the Q estimate at NKC is not reliable and an independent proof of the high attenuation finding given the geometrical and frequency constraints is still to be done. However, a high attenuation in the source volume of NW Bohemian swarms has been postulated before in relation to an expected, highly damaged zone bearing CO 2 at high pressure. The methods developed in the course of this thesis yield the potential to improve our understanding regarding the role of fluids and gases in intraplate event clustering.

Description: Erdbebenschwärme zeichnen sich durch eine große Anzahl an Ereignissen in einem relativ kleinen Zeitraum und Volumen aus. Im Gegensatz zu tektonischen Sequenzen ist in der Regel keine signifikantes Muster von Vor- und Nachbeben erkennbar. In Abwesenheit aktiven Vulkanismusses, kommen Erdbebenschwärme innerhalb kontinentaler Platten häufg an kontinentalen Verwerfungen vor, wie Beispielsweise im Bereich des Egergrabens im nordböhmischen Becken (Tschechien). Eine übliche Hypothese verbindet den Erdbebenentstehungsprozess mit Hochdruckfluiden. Der exakte kausale Zusammenhang ist jedoch häufig enigmatisch, da die zugrundeliegenden geotektonischen Prozesse im Vergleich zu tektonischen Sequenzen relativ langsam sind. Die hohe Erdbebenrate während aktiver Phasen stellt hohe Anforderungen an etablierte seismologische Routinen da diese häufg für Einzelereignisse konzipiert sind. So können sie einen hohen Aufwand bei manueller Selektion seismischer Phasen (picking) bedeuten oder rechenerisch aufwändig sein wenn volle Wellenformen verarbeitet werden sollen. Im Rahmen dieser methodologischen Thesis werden neue Ansätze zur Analyse seismischer Schwärme, sowie des zugrundeliegenden seismogenen Volumens entwickelt. Der Fokus liegt hierbei auf der gut untersuchten und überwachten nordböhmischen Schwarmregion. Ich entwickle und teste in dieser Arbeit einen innovativen Ansatz zur Detektion und Lokalisation von Erdbeben basierend auf einem tiefen konvolvierenden neuronalen Netzwerk. Diese Technologie bietet großes Potential da sie es erlaubt große Datenmengen effizient zu verarbeiten was durch die zunehmenden Datenmengen seismologischer Datenzentren immer weiter an Bedeutung gewinnt. Das entwickelte tiefe neuronale Netzwerk, trainiert auf Aufnahmen nordböhmischer Erdbebenschwärme, ist in der Lage 1000 Eregnisse in weniger als 1 Sekunde bei Verwendung voller Wellenformen zu lokalisieren und erreicht eine Präzision die vergleichbar ist mit der Genauigkeit eines Katalogs, der mittels Doppelte Differenzen Methode relokalisiert wurde. Eine weitere technologische Neuheit ist, dass die trainierten Filter der ersten Schicht des tiefen neuronalen Netzwerkes als Mustererkennungsfilter umfunktioniert werden und damit als Ereignisdetektor dienen können, ohne, dass zuvor explizit auf Rauschdaten trainiert werden muss. Für die weitere technologische Entwicklung stelle ich ein neues, automatisiertes Werkzeug für die synthetisierung realistischer Erdbebenschwarmkataloge, sowie hierauf basierender synthetischer voller Wollenform vor. Die Eingabeparameter werden durch die Geometrie des Quellvolumens, der Nukleationscharakteristik und Magnitude-Häufigkeitsverteilung definiert. Weiter können Rauschsignale realer Daten verwendet werden um äußerst realistische synthetische Daten zu generieren. Dieses Werkzeug wird verwendet um die Vollständigkeitmagnitude eines Detektors für volle Wellenformen anhand synthetischer Daten zu evaluieren. Die synthetisierten Daten sind Motiviert durch ein Hydrofrackingexperiment in Wysin (Polen). Des Weiteren stelle ich einen neuen Ansatz vor, der die Effekte der Wellenausbreitung zwischen Erdbeben und Stationen ausblendet und die Bestimmung der Dämpfung unmittelbar im Quellvolumen von Schwarmerdbeben erlaubt. Diese neue Methode benutzt das hochfrequente spektrale Verhältnis von Ereignispaaren mit gemeinsamen Strahlenwegen. Synthetische Tests zeigen, dass die Methode in der Lage ist die Dämpfung innerhalb des Quellvolumens mit hoher räumlicher Genauigkeit zu bestimmen. Weiter ist sie im Einzelnen unabhängig von der Entfernung zwischen Ereignis und Station als auch von der Komplexität der Dämpfungs und Geschwindigkeitsstruktur außerhalb des Quellvolumens. Die Anwendung auf Daten des nordböhmischen Erdbebenschwarms zeigt eine erhöhte P Phasen Dämpfung im Quellvolumen (Qp 〈 100) basierend auf Daten einer Station in der Nähe des Dorfes Luby (LBC). Die Wellenformen einer Station in unmittelbarer epizentraler Nähe, bei Novy Kostel (NKC), weisen eine relativ hohe Komplexität auf, was darauf hindeutet, dass seismische Wellen, die diese Station erreichen relativ stark gestreut werden im Vergleich zu anderen Stationen. Das hohe Maß an Komplexität destabilisiert die Methode und führt zu ungenauen Schätzungen an der Station NKC. Daher bedarf es einer weiteren unabhängigen Validierung der hohen Dämpfung bei gegebenen geometrischen und spektralen Voraussetzungen. Nichtsdestoweniger wurde bereits eine hohe Dämpfung im Quellvolumen der nordböhmischen Schwärme postuliert und erwartet, insbesondere im Zusammenhang mit einer Zone hoher Brüchigkeit die CO2 bei hohen Drücken beinhaltet. Die Methoden die im Rahmen dieser Thesis entwickelt werden haben das Potential unser Verständnis bezüglich der Rolle von Fluiden und Gasen bei Erdbebenschärmen innerhalb kontinentaler Platten zu verbessern.

Description: Partial melting is a first order process for the chemical differentiation of the crust (Vielzeuf et al., 1990). Redistribution of chemical elements during melt generation crucially influences the composition of the lower and upper crust and provides a mechanism to concentrate and transport chemical elements that may also be of economic interest. Understanding of the diverse processes and their controlling factors is therefore not only of scientific interest but also of high economic importance to cover the demand for rare metals. The redistribution of major and trace elements during partial melting represents a central step for the understanding how granite-bound mineralization develops (Hedenquist and Lowenstern, 1994). The partial melt generation and mobilization of ore elements (e.g. Sn, W, Nb, Ta) into the melt depends on the composition of the sedimentary source and melting conditions. Distinct source rocks have different compositions reflecting their deposition and alteration histories. This specific chemical “memory” results in different mineral assemblages and melting reactions for different protolith compositions during prograde metamorphism (Brown and Fyfe, 1970; Thompson, 1982; Vielzeuf and Holloway, 1988). These factors do not only exert an important influence on the distribution of chemical elements during melt generation, they also influence the volume of melt that is produced, extraction of the melt from its source, and its ascent through the crust (Le Breton and Thompson, 1988). On a larger scale, protolith distribution and chemical alteration (weathering), prograde metamorphism with partial melting, melt extraction, and granite emplacement are ultimately depending on a (plate-)tectonic control (Romer and Kroner, 2016). Comprehension of the individual stages and their interaction is crucial in understanding how granite-related mineralization forms, thereby allowing estimation of the mineralization potential of certain areas. Partial melting also influences the isotope systematics of melt and restite. Radiogenic and stable isotopes of magmatic rocks are commonly used to trace back the source of intrusions or to quantify mixing of magmas from different sources with distinct isotopic signatures (DePaolo and Wasserburg, 1979; Lesher, 1990; Chappell, 1996). These applications are based on the fundamental requirement that the isotopic signature in the melt reflects that of the bulk source from which it is derived. Different minerals in a protolith may have isotopic compositions of radiogenic isotopes that deviate from their whole rock signature (Ayres and Harris, 1997; Knesel and Davidson, 2002). In particular, old minerals with a distinct parent-to-daughter (P/D) ratio are expected to have a specific radiogenic isotope signature. As the partial melting reaction only involves selective phases in a protolith, the isotopic signature of the melt reflects that of the minerals involved in the melting reaction and, therefore, should be different from the bulk source signature. Similar considerations hold true for stable isotopes.

Description: Partielle Schmelzbildung ist ein zentraler Prozess für die geochemische Differentiation der Erdkruste (Vielzeuf et al., 1990). Die Umverteilung chemischer Elemente während der Schmelzbildung beeinflusst die Zusammensetzung der oberen und unteren Erdkruste entscheidend und stellt einen Mechanismus zur Konzentration und zum Transport chemischer Elemente dar. Das Verständnis der diversen Prozesse und der kontrollierenden Faktoren ist deshalb nicht nur von wissenschaftlichem Interesse sondern auch von ökonomischer Bedeutung um die Nachfrage für seltene Metalle zu decken. Die Umverteilung von Haupt- und Spurenelementen während des partiellen Aufschmelzens ist ein entscheidender Schritt für das Verständnis wie sich granitgebundene Lagerstätten bilden (Hedenquist and Lowenstern, 1994). Die Schmelzbildung und die Mobilisierung von Erz-Elementen (z. B. Sn, W, Nb, Ta) in die Schmelze hängt von der Zusammensetzung der sedimentären Ausgangsgesteine und den Schmelzbedingungen ab. Verschiedene Ausgangsgesteine haben aufgrund ihrer Ablagerungs- und Verwitterungsgeschichte unterschiedliche Zusammensetzungen. Dieses spezifische geochemische „Gedächtnis“ resultiert in unterschiedlichen Mineralparagenesen und Schmelzreaktionen in verschiedenen Ausgangsgesteinen während der prograden Metamorphose. (Brown and Fyfe, 1970; Thompson, 1982; Vielzeuf and Holloway, 1988). Diese Faktoren haben nicht nur einen wichtigen Einfluss auf die Verteilung chemischer Elemente während der Schmelzbildung, sie beeinflussen auch das Volumen an Schmelze, die Extraktion der Schmelze aus dem Ausgangsgestein und deren Aufstieg durch die Erdkruste (Le Breton and Thompson, 1988). Auf einer grösseren Skala unterliegen die Verteilung der Ausgangsgesteine und deren chemische Alteration (Verwitterung), die prograde Metamorphose mit partieller Schmelzbildung, Schmelzextraktion und die Platznahme granitischer Intrusionen einer plattentektonischen Kontrolle. Das Verständnis der einzelnen Schritte und deren Wechselwirkungen ist entscheidend um zu verstehen wie granitgebunden Lagerstätten entstehen und erlaubt es, das Mineralisierungspotential bestimmter Gebiete abzuschätzen. Partielles Aufschmelzen beeinflusst auch die Isotopensystematik der Schmelze und des Restites. Die Zusammensetzungen radiogener und stabiler Isotopen von magmatischen Gesteinen werden im Allgemeinen dazu verwendet um deren Ursprungsgesteine zu identifizieren oder um Mischungsprozesses von Magmen unterschiedlichen Ursprunges zu quantifizieren (DePaolo and Wasserburg, 1979; Lesher, 1990; Chappell, 1996). Diese Anwendungen basieren auf der fundamentalen Annahme, dass die Isotopenzusammensetzung der Schmelze derjenigen des Ausgangsgesteines entspricht. Unterschiedliche Minerale in einem Gestein können unterschiedliche, vom Gesamtgestein abweichende, Isotopenzusammensetzungen haben (Ayres and Harris, 1997; Knesel and Davidson, 2002). Insbesondere für alte Minerale, mit einem unterschiedlichen Mutter-Tochter Nuklidverhältnis, ist eine spezifische Isotopenzusammensetzung zu erwarten. Da im partiellen Schmelzprozess nur bestimmte Minerale eines Gesteines involviert sind, entspricht die Isotopenzusammensetzung der Schmelze derjenigen der Minerale welche an der Schmelzreaktion teilnehmen. Daher sollte die Isotopenzusammensetzung der Schmelze von derjenigen des Ursprungsgesteines abweichen. Ähnliche Überlegungen treffen auch für stabile Isotopen zu.

Description: SICOR is a sensor independent atmospheric correction of optical Earth Observation (EO) data from both multispectral and hyperspectral instruments. Currently, SICOR can be applied to Sentinel-2 and EnMAP data but the implementation of additional space- and airborne sensors is under development.

Description: CLI software for applying Fully Connected Neural Network (FNN) algorithm for LUCAS-2015

Description: Geochemical processes such as mineral dissolution and precipitation alter the microstructure of rocks, and thereby affect their hydraulic and mechanical behaviour. Quantifying these property changes and considering them in reservoir simulations is essential for a sustainable utilisation of the geological subsurface. Due to the lack of alternatives, analytical methods and empirical relations are currently applied to estimate evolving hydraulic and mechanical rock properties associated with chemical reactions. However, the predictive capabilities of analytical approaches remain limited, since they assume idealised microstructures, and thus are not able to reflect property evolution for dynamic processes. Hence, aim of the present thesis is to improve the prediction of permeability and stiffness changes resulting from pore space alterations of reservoir sandstones. A detailed representation of rock microstructure, including the morphology and connectivity of pores, is essential to accurately determine physical rock properties. For that purpose, three-dimensional pore-scale models of typical reservoir sandstones, obtained from highly resolved micro-computed tomography (micro-CT), are used to numerically calculate permeability and stiffness. In order to adequately depict characteristic distributions of secondary minerals, the virtual samples are systematically altered and resulting trends among the geometric, hydraulic, and mechanical rock properties are quantified. It is demonstrated that the geochemical reaction regime controls the location of mineral precipitation within the pore space, and thereby crucially affects the permeability evolution. This emphasises the requirement of determining distinctive porosity-permeability relationships by means of digital pore-scale models. By contrast, a substantial impact of spatial alterations patterns on the stiffness evolution of reservoir sandstones are only observed in case of certain microstructures, such as highly porous granular rocks or sandstones comprising framework-supporting cementations. In order to construct synthetic granular samples a process-based approach is proposed including grain deposition and diagenetic cementation. It is demonstrated that the generated samples reliably represent the microstructural complexity of natural sandstones. Thereby, general limitations of imaging techniques can be overcome and various realisations of granular rocks can be flexibly produced. These can be further altered by virtual experiments, offering a fast and cost-effective way to examine the impact of precipitation, dissolution or fracturing on various petrophysical correlations. The presented research work provides methodological principles to quantify trends in permeability and stiffness resulting from geochemical processes. The calculated physical property relations are directly linked to pore-scale alterations, and thus have a higher accuracy than commonly applied analytical approaches. This will considerably improve the predictive capabilities of reservoir models, and is further relevant to assess and reduce potential risks, such as productivity or injectivity losses as well as reservoir compaction or fault reactivation. Hence, the proposed method is of paramount importance for a wide range of natural and engineered subsurface applications, including geothermal energy systems, hydrocarbon reservoirs, CO2 and energy storage as well as hydrothermal deposit exploration

Description: Geochemische Lösungs- und Fällungsprozesse verändern die Struktur des Porenraums und können dadurch die hydraulischen und mechanischen Gesteinseigenschaften erheblich beeinflussen. Die Quantifizierung dieser Parameteränderung und ihre Berücksichtigung in Reservoirmodellen ist entscheidend für eine nachhaltige Nutzung des geologischen Untergrunds. Aufgrund fehlender Alternativen werden dafür bisher analytische Methoden genutzt. Da diese Ansätze eine idealisierte Mikrostruktur annehmen, können insbesondere Änderungen der Gesteinseigenschaften infolge von dynamischen Prozessen nicht zuverlässig abgebildet werden. Ziel der vorliegenden Doktorarbeit ist es deshalb, die Entwicklung von Gesteinspermeabilitäten und -steifigkeiten aufgrund von Porenraumveränderungen genauer vorherzusagen. Für die möglichst exakte Bestimmung physikalischer Gesteinsparameter ist eine detaillierte Darstellung der Mikrostruktur notwendig. Basierend auf mikro-computertomographischen Scans werden daher hochaufgelöste, dreidimensionale Modelle typischer Reservoirsandsteine erstellt und Gesteinspermeabilität und -steifigkeit numerisch berechnet. Um charakteristische Verteilungen von Sekundärmineralen abzubilden, wird der Porenraum dieser virtuellen Sandsteinproben systematisch verändert und die resultierenden Auswirkungen auf die granulometrischen, hydraulischen und elastischen Gesteinseigenschaften bestimmt. Die Ergebnisse zeigen deutlich, dass charakteristische Fällungsmuster unterschiedlicher geochemischer Reaktionsregime die Permeabilität erheblich beeinflussen. Folglich ist die Nutzung von porenskaligen Modellen zur Bestimmung der Porosität-Permeabilitätsbeziehungen unbedingt notwendig. Im Gegensatz dazu ist die Verteilung von Sekundärmineralen für die Gesteinssteifigkeit nur bei bestimmten Mikrostrukturen von Bedeutung, hierzu zählen hochporöse Sandsteine oder solche mit Korngerüst-stützenden Zementierungen. In der Arbeit wird außerdem ein Ansatz zur Konstruktion granularer Gesteine vorgestellt, welcher sowohl die Kornsedimentation als auch die diagenetische Verfestigung umfasst. Es wird gezeigt, dass die synthetischen Proben die mikrostrukturelle Komplexität natürlicher Reservoirsandsteine gut abbilden. Dadurch können generelle Limitationen von bildgebenden Verfahren überwunden und unterschiedlichste virtuelle Repräsentationen von granularen Gesteinen generiert werden. Die synthetischen Proben können zukünftig in virtuellen Experimenten verwendet werden, um die Auswirkungen von Lösungs- und Fällungsreaktionen auf verschiedene petrophysikalische Korrelationen zu untersuchen. Die vorgestellte Arbeit liefert methodische Grundlagen zur Quantifizierung von Permeabilitäts- und Steifigkeitsänderungen infolge geochemischer Prozesse. Die berechneten petrophysikalischen Beziehungen basieren direkt auf mikrostrukturellen Veränderungen des Porenraums. Daher bieten sie eine genauere Vorhersage der Gesteinseigenschaften als herkömmliche analytische Methoden, wodurch sich die Aussagekraft von Reservoirmodellen erheblich verbessert. Somit können Risiken, wie Produktivitäts- oder Injektivitätsverluste sowie Reservoirkompaktion oder Störungsreaktivierung, verringert werden. Die präsentierten Ergebnisse sind daher relevant für verschiedenste Bereiche der geologischen Untergrundnutzung wie CO2- oder Energiespeicherung, Geothermie, Kohlenwasserstoffgewinnung sowie die Erkundung hydrothermaler Lagerstätten.

Description: Due to the major role of greenhouse gas emissions in global climate change, the development of non-fossil energy technologies is essential. Deep geothermal energy represents such an alternative, which offers promising properties such as a high base load capability and a large untapped potential. The present work addresses barite precipitation within geothermal systems and the associated reduction in rock permeability, which is a major obstacle to maintaining high efficiency. In this context, hydro-geochemical models are essential to quantify and predict the effects of precipitation on the efficiency of a system. The objective of the present work is to quantify the induced injectivity loss using numerical and analytical reactive transport simulations. For the calculations, the fractured-porous reservoirs of the German geothermal regions North German Basin (NGB) and Upper Rhine Graben (URG) are considered. Similar depth-dependent precipitation potentials could be determined for both investigated regions (2.8-20.2 g/m3 fluid). However, the reservoir simulations indicate that the injectivity loss due to barite deposition in the NGB is significant (1.8%-6.4% per year) and the longevity of the system is affected as a result; this is especially true for deeper reservoirs (3000 m). In contrast, simulations of URG sites indicate a minor role of barite (〈 0.1%-1.2% injectivity loss per year). The key differences between the investigated regions are reservoir thicknesses and the presence of fractures in the rock, as well as the ionic strength of the fluids. The URG generally has fractured-porous reservoirs with much higher thicknesses, resulting in a greater distribution of precipitates in the subsurface. Furthermore, ionic strengths are higher in the NGB, which accelerates barite precipitation, causing it to occur more concentrated around the wellbore. The more concentrated the precipitates occur around the wellbore, the higher the injectivity loss. In this work, a workflow was developed within which numerical and analytical models can be used to estimate and quantify the risk of barite precipitation within the reservoir of geothermal systems. A key element is a newly developed analytical scaling score that provides a reliable estimate of induced injectivity loss. The key advantage of the presented approach compared to fully coupled reservoir simulations is its simplicity, which makes it more accessible to plant operators and decision makers. Thus, in particular, the scaling score can find wide application within geothermal energy, e.g., in the search for potential plant sites and the estimation of long-term efficiency.

Description: Aufgrund der tragenden Rolle der Treibhausgasemissionen für den globalen Klimawandel ist die Entwicklung von nicht-fossilen Energietechnologien essenziell. Die Tiefengeothermie stellt eine solche Alternative dar, welche vielversprechende Eigenschaften wie eine hohe Grundlastfähigkeit und ein großes ungenutztes Potenzial bietet. Die vorliegende Arbeit befasst sich mit Barytausfällungen inner- halb geothermaler Systeme und der damit einhergehenden Verringerung der Gesteinsdurchlässigkeit, welche ein Haupthindernis für die Aufrechterhaltung einer hohen Effizienz darstellen. Dabei sind hydro-geochemische Modelle unerlässlich, um die Auswirkungen von Ausfällungen auf die Effizienz eines Systems zu quantifizieren und vorherzusagen. Ziel der vorliegenden Arbeit ist es, mittels numerischer und analytischer reaktiver Transportsimulationen, den induzierten Injektivitätsverlust zu quantifizieren. Für die Berechnungen werden die klüftig-porösen Reservoire der deutschen Geothermieregionen Norddeutsches Becken (NDB) und Oberrheingraben (ORG) betrachtet. Für beide untersuchte Regionen konnte ein ähnliches, tiefenabhängiges Fällungspotenzial bestimmt werden (2,8–20,2 g/m3 Fluid). Die Reservoirsimulationen zeigen jedoch, dass der Injektivitätsverlust aufgrund von Barytablagerungen im NDB erheblich ist (1,8%–6,4% pro Jahr) und die Langlebigkeit der Anlage dadurch beeinträchtigt wird, dies gilt insbesondere für tiefere Reservoire (3000 m). Im Gegensatz dazu deuten die Simulationen der ORG-Standorte auf eine untergeordnete Rolle von Baryt hin (〈 0,1%–1,2% Injektivitätsverlust pro Jahr). Die entscheidenden Unterschiede zwischen den untersuchten Regionen sind die Reservoirmächtigkeiten und das Vorhandensein von Rissen im Gestein sowie die Ionenstärke der Fluide. Der ORG weist in der Regel klüftig-poröse Reservoire mit deutlich höheren Mächtigkeiten auf, was zu einer größeren Verteilung der Präzipitate im Untergrund führt. Weiterhin sind die Ionenstärken im NDB höher, was die Barytausfällung beschleunigt und diese dadurch konzentrierter um das Bohrloch herum entstehen lässt. Je konzentrierter die Präzipitate um die Bohrung herum auftreten, desto höher ist der Injektivitätsverlust. In dieser Arbeit wurde ein Workflow erarbeitet, innerhalb dessen mittels numerischer und analytischer Modelle das Risiko von Barytausfällungen innerhalb des Reservoirs geothermischer Systeme abgeschätzt und quantifiziert werden kann. Ein zentrales Element ist ein neu entwickelter, analytischer Scaling-Score, der eine zuverlässige Schätzung des induzierten Injektivitätsverlustes ermöglicht. Der entscheidende Vorteil des präsentierten Ansatzes im Vergleich zu voll-gekoppelten Reservoirsimulationen liegt in ihrer Einfachheit, die sie für Anlagenbetreiber und Entscheidungsträger zugänglicher macht. Somit kann insbesondere der Scaling-Score eine breite Anwendung innerhalb der Geothermie finden, z.B. bei der Suche nach potenziellen Anlagenstandorten und der Abschätzung der langfristigen Effizienz.