TY - THES A1 - Wriedt, Gunter T1 - Modelling of nitrogen transport and turnover during soil and groundwater passage in a small lowland catchment of Northern Germany N2 - Stoffumsatzreaktionen und hydraulische Prozesse im Boden und Grundwasser können in Tieflandeinzugsgebieten zu einer Nitratretention führen. Die Untersuchung dieser Prozesse in Raum und Zeit kann mit Hilfe geeigneter Modelle erfolgen. Ziele dieser Arbeit sind: i) die Entwicklung eines geeigneten Modellansatzes durch Kombination von Teilmodellen zur Simulation des N-Transportes im Boden und Grundwasser von Tieflandeinzugsgebieten und ii) die Untersuchung von Wechselwirkungen zwischen Gebietseigenschaften und N-Transport unter besonderer Berücksichtigung der potentiellen N-Zufuhr in die Oberflächengewässer. Der Modellansatz basiert auf der Kombination verschiedener Teilmodelle: das Bodenwasser- und -stickstoffmodell mRISK-N, das Grundwassermodell MODFLOW und das Stofftransportmodell RT3D. Zur Untersuchung der Wechselwirkungen mit den Gebietseigenschaften muss die Verteilung und Verfügbarkeit von Reaktionspartnern berücksichtigt werden. Dazu wurde ein Reaktionsmodul entwickelt, welches chemische Prozesse im Grundwasser simuliert. Hierzu gehören die Mineralisation organischer Substanz durch Sauerstoff, Nitrat und Sulfat sowie die Pyritoxidation durch Sauerstoff und Nitrat. Der Modellansatz wurde in verschiedenen Einzelstudien angewandt, wobei jeweils bestimmte Teilmodelle im Vordergrund stehen. Alle Modellstudien basieren auf Daten aus dem Schaugrabeneinzugsgebiet (ca. 25 km²), in der Nähe von Osterburg(Altmark) im Norden Sachsen-Anhalts. Die folgenden Einzelstudien wurden durchgeführt: i) Evaluation des Bodenmodells anhand von Lysimeterdaten, ii) Modellierung eines Tracerexperimentes im Feldmaßstab als eine erste Anwendung des Reaktionsmoduls, iii) Untersuchung hydraulisch-chemischer Wechselwirkungen an einem 2D-Grundwassertransekt, iv) Flächenverteilte Modellierung von Grundwasserneubildung und Bodenstickstoffaustrag im Untersuchungsgebiet als Eingangsdaten für nachfolgende Grundwassersimulationen, und v) Untersuchung der Ausbreitung von Nitrat im Grundwasser und des Durchbruchs in die Oberflächengewässer im Untersuchungsgebiet auf Basis einer 3D-Modellierung von Grundwasserströmung und reaktivem Stofftransport. Die Modellstudien zeigen, dass der Modellansatz geeignet ist, die Wechselwirkungen zwischen Stofftransport und –umsatz und den hydraulisch-chemischen Gebietseigenschaften zu modellieren. Die Ausbreitung von Nitrat im Sediment wird wesentlich von der Verfügbarkeit reaktiver Substanzen sowie der Verweilzeit im Grundwasserleiter bestimmt. Bei der Simulation des Untersuchungsgebietes wurde erst nach 70 Jahren eine der gegebenen Eintragssitutation entsprechende Nitratkonzentration im Grundwasserzustrom zum Grabensystem erreicht (konservativer Transport). Die Berücksichtigung von reaktivem Stofftransport führt zu einer deutlichen Reduktion der Nitratkonzentrationen. Die Modellergebnisse zeigen, dass der Grundwasserzustrom die beobachtete Nitratbelastung im Grabensystem nicht erklären kann, da der Großteil des Nitrates durch Denitrifikation verloren geht. Andere Quellen, wie direkte Einträge oder Dränagenzuflüsse müssen ebenfalls in Betracht gezogen werden. Die Prognosefähigkeit des Modells für das Untersuchungsgebiet wird durch die Datenunsicherheiten und die Schätzung der Modellparameter eingeschränkt. Dennoch ist der Modellansatz eine wertvolle Hilfe bei der Identifizierung von belastungsrelevanten Teilflächen (Stoffquellen und -senken) sowie bei der Modellierung der Auswirkungen von Managementmaßnahmen oder Landnutzungsveränderungen auf Grundlage von Szenario-Simulationen. Der Modellansatz unterstützt auch die Interpretation von Beobachtungsdaten, da so die lokalen Informationen in einen räumlichen und zeitlichen Zusammenhang gestellt werden können. N2 - Chemical transformations and hydraulic processes in soil and groundwater often lead to an apparent retention of nitrate in lowland catchments. Models are needed to evaluate the interaction of these processes in space and time. The objectives of this study are i) to develop a specific modelling approach by combining selected modelling tools simulating N-transport and turnover in soils and groundwater of lowland catchments, ii) to study interactions between catchment properties and nitrogen transport. Special attention was paid to potential N-loads to surface waters. The modelling approach combines various submodels for water flow and solute transport in soil and groundwater: The soil-water- and nitrogen-model mRISK-N, the groundwater flow model MODFLOW and the solute transport model RT3D. In order to investigate interactions of N-transport and catchment characteristics, the distribution and availability of reaction partners have to be taken into account. Therefore, a special reaction-module is developed, which simulates various chemical processes in groundwater, such as the degradation of organic matter by oxygen, nitrate, sulphate or pyrite oxidation by oxygen and nitrate. The model approach is applied to different simulation, focussing on specific submodels. All simulation studies are based on field data from the Schaugraben catchment, a pleistocene catchment of approximately 25 km², close to Osterburg(Altmark) in the North of Saxony-Anhalt. The following modelling studies have been carried out: i) evaluation of the soil-water- and nitrogen-model based on lysimeter data, ii) modelling of a field scale tracer experiment on nitrate transport and turnover in the groundwater as a first application of the reaction module, iii) evaluation of interactions between hydraulic and chemical aquifer properties in a two-dimensional groundwater transect, iv) modelling of distributed groundwater recharge and soil nitrogen leaching in the study area, to be used as input data for subsequent groundwater simulations, v) study of groundwater nitrate distribution and nitrate breakthrough to the surface water system in the Schaugraben catchment area and a subcatchment, using three-dimensional modelling of reactive groundwater transport. The various model applications prove the model to be capable of simulating interactions between transport, turnover and hydraulic and chemical catchment properties. The distribution of nitrate in the sediment and the resulting loads to surface waters are strongly affected by the amount of reactive substances and by the residence time within the aquifer. In the Schaugraben catchment simulations, it is found that a period of 70 years is needed to raise the average seepage concentrations of nitrate to a level corresponding to the given input situation, if no reactions are considered. Under reactive transport conditions, nitrate concentrations are reduced effectively. Simulation results show that groundwater exfiltration does not contribute considerably to the nitrate pollution of surface waters, as most nitrate entering soils and groundwater is lost by denitrification. Additional sources, such as direct inputs or tile drains have to be taken into account to explain surface water loads. The prognostic value of the models for the study site is limited by uncertainties of input data and estimation of model parameters. Nevertheless, the modelling approach is a useful aid for the identification of source and sink areas of nitrate pollution as well as the investigation of system response to management measures or landuse changes with scenario simulations. The modelling approach assists in the interpretation of observed data, as it allows to integrate local observations into a spatial and temporal framework. KW - Stickstoff KW - Nitrat KW - Modellierung KW - Grundwasser KW - Einzugsgebiet KW - Reaktiver Stofftransport KW - Denitrifikation KW - Nitrogen KW - Nitrate KW - modelling KW - groundwater KW - catchment KW - reactive transport KW - denitrification Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-0001307 ER - 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 -