TY - THES A1 - Tranter, Morgan Alan T1 - Numerical quantification of barite reservoir scaling and the resulting injectivity loss in geothermal systems N2 - 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. N2 - 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. KW - geothermal energy KW - formation damage KW - reactive transport KW - radial flow KW - barite KW - Geothermie KW - radiale Strömung KW - Baryt KW - reaktiver Transport KW - Formationsschaden Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-561139 ER - TY - THES A1 - Steding, Svenja T1 - Geochemical and Hydraulic Modeling of Cavernous Structures in Potash Seams T1 - Geochemische und hydraulische Modellierung kavernöser Strukturen in Kaliflözen N2 - Salt deposits offer a variety of usage types. These include the mining of rock salt and potash salt as important raw materials, the storage of energy in man-made underground caverns, and the disposal of hazardous substances in former mines. The most serious risk with any of these usage types comes from the contact with groundwater or surface water. It causes an uncontrolled dissolution of salt rock, which in the worst case can result in the flooding or collapse of underground facilities. Especially along potash seams, cavernous structures can spread quickly, because potash salts show a much higher solubility than rock salt. However, as their chemical behavior is quite complex, previous models do not account for these highly soluble interlayers. Therefore, the objective of the present thesis is to describe the evolution of cavernous structures along potash seams in space and time in order to improve hazard mitigation during the utilization of salt deposits. The formation of cavernous structures represents an interplay of chemical and hydraulic processes. Hence, the first step is to systematically investigate the dissolution and precipitation reactions that occur when water and potash salt come into contact. For this purpose, a geochemical reaction model is used. The results show that the minerals are only partially dissolved, resulting in a porous sponge like structure. With the saturation of the solution increasing, various secondary minerals are formed, whose number and type depend on the original rock composition. Field data confirm a correlation between the degree of saturation and the distance from the center of the cavern, where solution is entering. Subsequently, the reaction model is coupled with a flow and transport code and supplemented by a novel approach called ‘interchange’. The latter enables the exchange of solution and rock between areas of different porosity and mineralogy, and thus ultimately the growth of the cavernous structure. By means of several scenario analyses, cavern shape, growth rate and mineralogy are systematically investigated, taking also heterogeneous potash seams into account. The results show that basically four different cases can be distinguished, with mixed forms being a frequent occurrence in nature. The classification scheme is based on the dimensionless numbers Péclet and Damköhler, and allows for a first assessment of the hazard potential. In future, the model can be applied to any field case, using measurement data for calibration. The presented research work provides a reactive transport model that is able to spatially and temporally characterize the propagation of cavernous structures along potash seams for the first time. Furthermore, it allows to determine thickness and composition of transition zones between cavern center and unaffected salt rock. The latter is particularly important in potash mining, so that natural cavernous structures can be located at an early stage and the risk of mine flooding can thus be reduced. The models may also contribute to an improved hazard prevention in the construction of storage caverns and the disposal of hazardous waste in salt deposits. Predictions regarding the characteristics and evolution of cavernous structures enable a better assessment of potential hazards, such as integrity or stability loss, as well as of suitable mitigation measures. N2 - Salzlagerstätten bieten eine Vielzahl an Nutzungsmöglichkeiten. Diese umfassen den Abbau von Steinsalz und Kalisalz als wichtige Rohstoffe, die Speicherung von Energie in künstlich erzeugten Hohlräumen, sowie die Entsorgung gefährlicher Substanzen in stillgelegten Bergwerken. Die größte Gefahr bei jeder dieser Nutzungsarten ist der Kontakt mit Grund- oder Oberflächenwasser. Er bewirkt eine unkontrollierte Lösung des Salzgesteins, was im schlimmsten Fall zur Flutung oder zum Einsturz unterirdischer Infrastrukturen führt. Insbesondere entlang von Kaliflözen können sich kavernöse Strukturen schnell ausbreiten, da Kalisalze eine wesentlich höhere Löslichkeit besitzen als Steinsalz. Ihr chemisches Verhalten ist jedoch komplex, weshalb bisherige Modelle diese hochlöslichen Zwischenschichten vernachlässigen. Ziel der vorliegenden Doktorarbeit ist es daher, die Ausbreitung kavernöser Strukturen entlang von Kaliflözen räumlich und zeitlich zu beschreiben und damit die Möglichkeiten zur Gefahrenprävention bei der Nutzung von Salzlagerstätten zu verbessern. Die Bildung kavernöser Strukturen ist ein Zusammenspiel chemischer und hydraulischer Prozesse. Zunächst wird daher mithilfe eines geochemischen Reaktionsmodells systematisch untersucht, welche Lösungs- und Fällungsreaktionen beim Kontakt von Wasser und Kalisalz auftreten. Die Ergebnisse zeigen, dass nur ein Teil der Minerale gelöst wird, wodurch sich eine poröse, schwammartige Struktur bildet. Mit zunehmender Aufsättigung der Lösung treten verschiedene Sekundärminerale auf, deren Anzahl und Art vom Ausgangsgestein abhängen. Felddaten belegen dabei eine Korrelation zwischen Sättigungsgrad und Abstand vom Kavernenzentrum, wo die Lösung ein- und austritt. Anschließend wird das Reaktionsmodell mit einem Strömungs- und Transportcode gekoppelt und um einen neuartigen Ansatz namens "interchange" ergänzt. Dieser ermöglicht den Austausch von Lösung und Gestein zwischen Bereichen unterschiedlicher Porosität und Mineralogie, und damit letztlich das Wachstum der kavernösen Struktur. In mehreren Szenarienanalysen werden Kavernenform, Ausbreitungsgeschwindigkeit und Mineralogie systematisch untersucht und dabei auch heterogene Kaliflöze betrachtet. Die Ergebnisse zeigen, dass grundsätzlich vier Fälle zu unterscheiden sind, wobei in der Natur häufig Mischformen auftreten. Die Klassifizierung erfolgt auf Basis der dimensionslosen Kennzahlen Péclet und Damköhler und ermöglicht eine erste Abschätzung des Gefahrenpotentials. In Zukunft kann das Modell auf beliebige Feldbeispiele angewandt und mithilfe von Messdaten kalibriert werden. Die vorliegende Arbeit liefert ein reaktives Transportmodell, mit dem die Ausbreitung kavernöser Strukturen entlang von Kaliflözen erstmals räumlich und zeitlich beschrieben werden kann. Auch Mächtigkeit und Zusammensetzung der Übergangszone zwischen Kavernenzentrum und unberührtem Salzgestein können damit bestimmt werden. Letzteres ist insbesondere im Kalibergbau von Bedeutung, um natürliche kavernöse Strukturen rechtzeitig zu lokalisieren und damit das Risiko für eine Flutung von Bergwerken zu verringern. Auch bei der Herstellung von Speicherkavernen und der Einlagerung gefährlicher Substanzen im Salzgestein können die Modelle zu einer besseren Gefahrenprävention beitragen. Sie ermöglichen Prognosen über Beschaffenheit und Ausbreitungsverhalten kavernöser Strukturen, wodurch sowohl potentielle Gefahren, wie der Verlust von Dichtigkeit oder Stabilität, als auch geeignete Gegenmaßnahmen besser abschätzbar werden. KW - reactive transport KW - reaktiver Transport KW - salt rock KW - Salzgestein KW - water rock interactions KW - Wasser-Gesteins-Wechselwirkungen KW - density-driven flow KW - dichtegetriebene Strömung KW - PHREEQC Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-548182 ER -