TY - JOUR A1 - Wetzel, Maria A1 - Kempka, Thomas A1 - Kühn, Michael T1 - Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations JF - Materials N2 - The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution. KW - digital rock physics KW - micro-CT KW - elastic properties KW - numerical simulation KW - chemical-mechanical interaction KW - Code_Aster KW - composite properties Y1 - 2018 U6 - https://doi.org/10.3390/ma11040542 SN - 1996-1944 VL - 11 IS - 4 PB - MDPI CY - Basel ER - TY - GEN A1 - Wetzel, Maria A1 - Kempka, Thomas A1 - Kühn, Michael T1 - Quantifying rock weakening due to decreasing calcite mineral content by numerical simulations T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - The quantification of changes in geomechanical properties due to chemical reactions is of paramount importance for geological subsurface utilisation, since mineral dissolution generally reduces rock stiffness. In the present study, the effective elastic moduli of two digital rock samples, the Fontainebleau and Bentheim sandstones, are numerically determined based on micro-CT images. Reduction in rock stiffness due to the dissolution of 10% calcite cement by volume out of the pore network is quantified for three synthetic spatial calcite distributions (coating, partial filling and random) using representative sub-cubes derived from the digital rock samples. Due to the reduced calcite content, bulk and shear moduli decrease by 34% and 38% in maximum, respectively. Total porosity is clearly the dominant parameter, while spatial calcite distribution has a minor impact, except for a randomly chosen cement distribution within the pore network. Moreover, applying an initial stiffness reduced by 47% for the calcite cement results only in a slightly weaker mechanical behaviour. Using the quantitative approach introduced here substantially improves the accuracy of predictions in elastic rock properties compared to general analytical methods, and further enables quantification of uncertainties related to spatial variations in porosity and mineral distribution. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1092 KW - digital rock physics KW - micro-CT KW - elastic properties KW - numerical simulation KW - chemical-mechanical interaction KW - Code_Aster KW - composite properties Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-473089 SN - 1866-8372 IS - 1092 ER - TY - THES A1 - Otto, Christopher T1 - Numerical analysis of thermal, hydraulic and mechanical processes in the near- and far-field of underground coal gasification reactors T1 - Numerische Analyse von thermalen, hydraulischen und mechanischen Prozessen im Nah- und Fernfeld von Untertage-Kohlevergasungs-Reaktoren N2 - Die Untertagevergasung von Kohle (UTV) ermöglicht die Erschließung konventionell nicht förderbarer Kohleressourcen und bietet dadurch Potenzial zur Erhöhung der weltweiten Kohlereserven. Bei der in-situ Kohleumwandlung entsteht ein hochkalorisches Synthesegas, das elektrifiziert oder zur Gewinnung chemischer Rohstoffe und synthetischer Kraftstoffe eingesetzt werden kann. Neben den wirtschaftlichen Möglichkeiten, bestehen jedoch auch standort-spezifische Umweltgefährdungspotentiale durch Subsidenz und Schadstoffmigration von UTV-Rückständen in nutzbare Grundwasserleiter. Eine nachhaltige und effiziente UTV erfordert ein umfangreiches Verständnis der thermisch, hydraulisch und mechanisch gekoppelten Prozesse im UTV-Reaktornahbereich. Aufgrund der hohen Investitionskosten von UTV-Pilotanlagen, sind numerische Simulationen gekoppelter Prozesse von entscheidender Bedeutung für die Bewertung möglicher UTV-Umweltauswirkungen. Im Rahmen dieser Arbeit wird die UTV-induzierte Permeabilitätsveränderung, Erzeugung möglicher hydraulischer Kurzschlüsse benachbarter Reaktoren und Dynamik nicht-isothermer Multiphasenflüsse mit gekoppelten Simulationen analysiert. Die Simulationsergebnisse zeigen, dass eine Implementierung temperaturabhängiger thermo-mechanischer Gesteinsparameter nur für Untersuchungen von Permeabilitäts-änderungen im Reaktornachbereich notwendig ist. Die Ergebnisse erlauben somit eine recheneffiziente Realisierung von komplexen thermo-mechanisch gekoppelten Simulations-studien regionalskaliger Modelle mit konstanten Gesteinsparametern, bei nahezu gleichbleibender Ergebnisgenauigkeit, die zur Bewertung von UTV-Umweltgefährdungs-potenzialen beitragen. Simulationen zur Ausbildung hydraulischer Kurzschlüsse zwischen einzelnen UTV-Reaktoren auf regionaler Skala, verdeutlichen die Relevanz von geologischen Störungen an einem UTV-Standort, da diese durch Reaktivierung hydraulische Verbindungen induzieren und somit einen effizienten und nachhaltigen UTV-Betrieb negativ beeinträchtigen können. In diesem Zusammenhang kommt der Ausbildung einer Wasserdampfphase, der sogenannte „steam jacket“, im Hochtemperaturnahbereich von UTV-Reaktoren, als potenzielle Barriere zur Vermeidung von UTV-Schadstoffaustritten und zur potenziellen Minimierung von Energieverlusten eine entscheidende Bedeutung zu. Diese steam jackets entstehen durch evaporiertes Formationswasser und sind komplexe nicht-isotherme Multiphasenfluss-Phänomene. Für ein verbessertes Prozessverständnis dieser Multiphasenflüsse, wurde ein neuartiges Modellkonzept entwickelt, welches, validiert gegen Feldversuchsdaten, erstmals sowohl eine Quantifizierung als auch Prognose von Wasserflussraten in und aus einem UTV-Reaktor erlaubt. Die Ergebnisse der vorgelegten Doktorarbeit bilden eine wichtige Grundlage für eine erfolgreiche Integration gekoppelter thermo-hydro-mechanischer Simulationen in weiterführende Studien. Vor dem Hintergrund hoher UTV-Umweltgefährdungspotentiale, können diese zur verbesserten Bewertung und Minderung von UTV-Umweltauswirkungen beitragen, sowie die UTV-Effizienz nachhaltig optimieren. N2 - Underground coal gasification (UCG) has the potential to increase worldwide coal reserves by developing coal resources, currently not economically extractable by conventional mining methods. For that purpose, coal is combusted in situ to produce a high-calorific synthesis gas with different end-use options, including electricity generation as well as production of fuels and chemical feedstock. Apart from the high economic potentials, UCG may induce site‐specific environmental impacts, including ground surface subsidence and pollutant migration of UCG by-products into shallow freshwater aquifers. Sustainable and efficient UCG operation requires a thorough understanding of the coupled thermal, hydraulic and mechanical processes, occurring in the UCG reactor vicinity. The development and infrastructure costs of UCG trials are very high; therefore, numerical simulations of coupled processes in UCG are essential for the assessment of potential environmental impacts. Therefore, the aim of the present study is to assess UCG-induced permeability changes, potential hydraulic short circuit formation and non-isothermal multiphase fluid flow dynamics by means of coupled numerical simulations. Simulation results on permeability changes in the UCG reactor vicinity demonstrate that temperature-dependent thermo-mechanical parameters have to be considered in near-field assessments, only. Hence, far-field simulations do not become inaccurate, but benefit from increased computational efficiency when thermo-mechanical parameters are maintained constant. Simulations on potential hydraulic short circuit formation between single UCG reactors at regional-scale emphasize that geologic faults may induce hydraulic connections, and thus compromise efficient UCG operation. In this context, the steam jacket surrounding high-temperature UCG reactors plays a vital role in avoiding UCG by-products escaping into freshwater aquifers and in minimizing energy consumption by formation fluid evaporation. A steam jacket emerges in the close reactor vicinity due to phase transition of formation water and is a non-isothermal flow phenomenon. Considering this complex multiphase flow behavior, an innovative conceptual modeling approach, validated against field data, enables the quantification and prediction of UCG reactor water balances. The findings of this doctoral thesis provide an important basis for integration of thermo-hydro-mechanical simulations in UCG, required for the assessment and mitigation of its potential environmental impacts as well as optimization of its efficiency. KW - underground coal gasification KW - numerical simulation KW - non-isothermal multiphase flow KW - thermo-mechanical modeling KW - Untertage-Kohlevergasung KW - numerische Simulation KW - nichtisothermer Mehrphasenfluss KW - thermo-mechanische Modellierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-404625 ER - TY - JOUR A1 - Li, Zhen A1 - Spangenberg, Erik A1 - Schicks, Judith Maria A1 - Kempka, Thomas T1 - Numerical simulation of hydrate formation in the LArge-Scale Reservoir Simulator (LARS) JF - Energies : open-access journal of related scientific research, technology development and studies in policy and management N2 - The LArge-scale Reservoir Simulator (LARS) has been previously developed to study hydrate dissociation in hydrate-bearing systems under in-situ conditions. In the present study, a numerical framework of equations of state describing hydrate formation at equilibrium conditions has been elaborated and integrated with a numerical flow and transport simulator to investigate a multi-stage hydrate formation experiment undertaken in LARS. A verification of the implemented modeling framework has been carried out by benchmarking it against another established numerical code. Three-dimensional (3D) model calibration has been performed based on laboratory data available from temperature sensors, fluid sampling, and electrical resistivity tomography. The simulation results demonstrate that temperature profiles, spatial hydrate distribution, and bulk hydrate saturation are consistent with the observations. Furthermore, our numerical framework can be applied to calibrate geophysical measurements, optimize post-processing workflows for monitoring data, improve the design of hydrate formation experiments, and investigate the temporal evolution of sub-permafrost methane hydrate reservoirs. KW - methane hydrate KW - temperature sensor KW - electrical resistivity tomography KW - hydrate formation KW - numerical simulation Y1 - 2022 U6 - https://doi.org/10.3390/en15061974 SN - 1996-1073 VL - 15 IS - 6 PB - MDPI CY - Basel ER - TY - THES A1 - Li, Zhen T1 - Formation of Sub-Permafrost Methane Hydrate Reproduced by Numerical Modeling BT - Insights from LArge-scale Reservoir Simulator (LARS) to Mallik Site, Canadian Arctic N2 - Natural gas hydrates are ice-like crystalline compounds containing water cavities that trap natural gas molecules like methane (CH4), which is a potent greenhouse gas with high energy density. The Mallik site at the Mackenzie Delta in the Canadian Arctic contains a large volume of technically recoverable CH4 hydrate beneath the base of the permafrost. Understanding how the sub-permafrost hydrate is distributed can aid in searching for the ideal locations for deploying CH4 production wells to develop the hydrate as a cleaner alternative to crude oil or coal. Globally, atmospheric warming driving permafrost thaw results in sub-permafrost hydrate dissociation, releasing CH4 into the atmosphere to intensify global warming. It is therefore crucial to evaluate the potential risk of hydrate dissociation due to permafrost degradation. To quantitatively predict hydrate distribution and volume in complex sub-permafrost environments, a numerical framework was developed to simulate sub-permafrost hydrate formation by coupling the equilibrium CH4-hydrate formation approach with a fluid flow and transport simulator (TRANSPORTSE). In addition, integrating the equations of state describing ice melting and forming with TRANSPORTSE enabled this framework to simulate the permafrost evolution during the sub-permafrost hydrate formation. A modified sub-permafrost hydrate formation mechanism for the Mallik site is presented in this study. According to this mechanism, the CH4-rich fluids have been vertically transported since the Late Pleistocene from deep overpressurized zones via geologic fault networks to form the observed hydrate deposits in the Kugmallit–Mackenzie Bay Sequences. The established numerical framework was verified by a benchmark of hydrate formation via dissolved methane. Model calibration was performed based on laboratory data measured during a multi-stage hydrate formation experiment undertaken in the LArge scale Reservoir Simulator (LARS). As the temporal and spatial evolution of simulated and observed hydrate saturation matched well, the LARS model was therefore validated. This laboratory-scale model was then upscaled to a field-scale 2D model generated from a seismic transect across the Mallik site. The simulation confirmed the feasibility of the introduced sub-permafrost hydrate formation mechanism by demonstrating consistency with field observations. The 2D model was extended to the first 3D model of the Mallik site by using well-logs and seismic profiles, to investigate the geologic controls on the spatial hydrate distribution. An assessment of this simulation revealed the hydraulic contribution of each geological element, including relevant fault networks and sedimentary sequences. Based on the simulation results, the observed heterogeneous distribution of sub-permafrost hydrate resulted from the combined factors of the source-gas generation rate, subsurface temperature, and the permeability of geologic elements. Analysis of the results revealed that the Mallik permafrost was heated by 0.8–1.3 °C, induced by the global temperature increase of 0.44 °C and accelerated by Arctic amplification from the early 1970s to the mid-2000s. This study presents a numerical framework that can be applied to study the formation of the permafrost-hydrate system from laboratory to field scales, across timescales ranging from hours to millions of years. Overall, these simulations deepen the knowledge about the dominant factors controlling the spatial hydrate distribution in sub-permafrost environments with heterogeneous geologic elements. The framework can support improving the design of hydrate formation experiments and provide valuable contributions to future industrial hydrate exploration and exploitation activities. N2 - Gashydrate sind eisähnliche kristalline Verbindungen, die Moleküle wie Methan (CH4) in Hohlräumen einschließen. Die Mallik-Lagerstätte im Mackenzie-Delta in der kanadischen Arktis enthält ein großes Volumen an technisch förderbarem CH4-Hydrat unter dem Permafrostboden. Das Verständnis, wie die Hydrate verteilt sind, kann bei der Suche nach idealen Standorten für Förderbohrungen zu ihrer Erschließung als saubere Alternative zu Erdöl oder Kohle helfen. Weltweit führt die Erwärmung der Atmosphäre zum Auftauen des Permafrosts und zur Zersetzung der Hydrate, wodurch CH4 in die Atmosphäre freigesetzt und die globale Erwärmung verstärkt wird. Es ist also entscheidend, das potenzielle Risiko der Hydratauflösung aufgrund der Permafrostdegradation zu bewerten. Um die Verteilung und das Volumen von Hydraten in komplexen Sub-Permafrost Umgebungen quantitativ vorherzusagen, wurde ein numerischer Ansatz zur Simulation entwickelt. Hierzu wurde der Gleichgewichtsansatz für die CH4-Hydratbildung mit einem Strömungs- und Transportsimulator (TRANSPORTSE) kombiniert. Die zusätzliche Integrierung der Zustandsgleichungen, die das Schmelzen und die Bildung von Eis beschreiben, ermöglichte die Simulation der Permafrostentwicklung während der Hydratbildung. Für den Standort Mallik wird ein modifizierter Bildungsmechanismus in dieser Studie beschrieben. Demzufolge wurden die CH4-reichen Fluide seit dem späten Pleistozän aus tiefen Überdruckszonen vertikal über geologische Verwerfungssysteme transportiert, und haben die Hydratvorkommen gebildet. Der numerische Ansatz wurde anhand eines Benchmarks zur Hydratbildung verifiziert. Messdaten eines mehrstufigen Hydratbildungsexperiments im LArge scale Reservoir Simulator (LARS) dienten zur Kalibrierung. Basierend auf der guten Übereinstimmung zwischen der simulierten und beobachteten Hydratsättigung, wurde das LARS-Modell validiert. Im Anschluss erfolgte die Übertragung auf ein 2D-Modell im Feldmaßstab, das mithilfe einer seismischen Transekte durch den Mallik-Standort erstellt wurde. Die Übereinstimmung mit den Feldbeobachtungen bestätigte den beschriebenen Mechanismus zur Hydratbildung unterhalb des Permafrosts. Das 2D-Modell wurde basierend auf Bohrlochprotokollen und seismischen Profilen zum ersten 3D-Modell des Mallik-Standorts erweitert, um die geologischen Einflüsse auf die Hydratverteilung zu untersuchen. Die Auswertung verdeutlichte den Beitrag jedes geologischen Elements zum hydraulischen System, einschließlich relevanter Verwerfungssysteme und sedimentärer Abfolgen. Die beobachtete heterogene räumliche Verteilung der Hydrate ist auf die Gasproduktionsrate der Quelle, die Untergrundtemperatur und die Durchlässigkeit der geologischen Einheiten zurückzuführen. Die Analyse der Ergebnisse ergab, dass der Mallik-Permafrost um 0,8–1,3 °C erwärmt wurde, was durch den globalen Temperaturanstieg von 0,44 °C verursacht und durch die sogenannte polare Verstärkung seit Anfang der 1970er bis Mitte der 2000er Jahre beschleunigt wurde. Der in dieser Studie entwickelte numerische Ansatz zur Bildung von Permafrost-Hydrat-Systemen kann vom Labor- bis zum Feldmaßstab und über Zeitskalen von Stunden bis zu Millionen von Jahren angewendet werden. Mit den Simulationen konnten die dominierenden Faktoren identifiziert werden, welche die räumliche Hydratverteilung in Umgebungen mit heterogenen geologischen Strukturen steuern. Der Ansatz kann die Planung von Hydratbildungsexperimenten verbessern und einen wertvollen Beitrag für zukünftige industrielle Hydraterkundungen- und -erschließungen leisten. T2 - Bildung von Sub-Permafrost-Methanhydraten dargestellt durch numerische Modellierung: Erkenntnisse aus dem LArge-sclae Reservoir Simulator (LARS), angewandt auf den Mallik-Standort, Kanadische Arktis KW - methane hydrate KW - geologic fault KW - numerical simulation KW - hydrate formation KW - climate change KW - Mackenzie Delta KW - Mackenzie-Delta KW - Klimawandel KW - geologische Verwerfung KW - Hydratbildung KW - Methanhydrat KW - numerische Simulation Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-603302 ER - TY - GEN A1 - Kühn, Michael A1 - Kempka, Thomas A1 - de Lucia, Marco A1 - Scheck-Wenderoth, Magdalena T1 - Dissolved CO2 storage in geological formations with low pressure, low risk and large capacities T2 - Energy procedia N2 - Geological CO2 storage is a mitigation technology to reduce CO2 emissions from fossil fuel combustion. However, major concerns are the pressure increase and saltwater displacement in the mainly targeted deep groundwater aquifers due to injection of supercritical CO2. The suggested solution is storage of CO2 exclusively in the dissolved state. In our exemplary regional case study of the North East German Basin based on a highly resolved temperature and pressure distribution model and a newly developed reactive transport coupling, we have quantified that 4.7 Gt of CO2 can be stored in solution compared to 1.5 Gt in the supercritical state. KW - carbon dioxide KW - dissolved KW - storage capacity KW - numerical simulation KW - saline aquifer KW - Buntsandstein Y1 - 2017 U6 - https://doi.org/10.1016/j.egypro.2017.03.1607 SN - 1876-6102 VL - 114 SP - 4722 EP - 4727 PB - Elsevier CY - Amsterdam ER - TY - THES A1 - Jappsen, Anne-Katharina T1 - Present and early star formation : a study on rotational and thermal properties T1 - Sternentstehung im heutigen und im frühen Universum : eine Studie über Rotationseigenschaften und thermisches Verhalten N2 - We investigate the rotational and thermal properties of star-forming molecular clouds using hydrodynamic simulations. Stars form from molecular cloud cores by gravoturbulent fragmentation. Understanding the angular momentum and the thermal evolution of cloud cores thus plays a fundamental role in completing the theoretical picture of star formation. This is true not only for current star formation as observed in regions like the Orion nebula or the ρ-Ophiuchi molecular cloud but also for the formation of stars of the first or second generation in the universe. In this thesis we show how the angular momentum of prestellar and protostellar cores evolves and compare our results with observed quantities. The specific angular momentum of prestellar cores in our models agree remarkably well with observations of cloud cores. Some prestellar cores go into collapse to build up stars and stellar systems. The resulting protostellar objects have specific angular momenta that fall into the range of observed binaries. We find that collapse induced by gravoturbulent fragmentation is accompanied by a substantial loss of specific angular momentum. This eases the "angular momentum problem" in star formation even in the absence of magnetic fields. The distribution of stellar masses at birth (the initial mass function, IMF) is another aspect that any theory of star formation must explain. We focus on the influence of the thermodynamic properties of star-forming gas and address this issue by studying the effects of a piecewise polytropic equation of state on the formation of stellar clusters. We increase the polytropic exponent γ from a value below unity to a value above unity at a certain critical density. The change of the thermodynamic state at the critical density selects a characteristic mass scale for fragmentation, which we relate to the peak of the IMF observed in the solar neighborhood. Our investigation generally supports the idea that the distribution of stellar masses depends mainly on the thermodynamic state of the gas. A common assumption is that the chemical evolution of the star-forming gas can be decoupled from its dynamical evolution, with the former never affecting the latter. Although justified in some circumstances, this assumption is not true in every case. In particular, in low-metallicity gas the timescales for reaching the chemical equilibrium are comparable or larger than the dynamical timescales. In this thesis we take a first approach to combine a chemical network with a hydrodynamical code in order to study the influence of low levels of metal enrichment on the cooling and collapse of ionized gas in small protogalactic halos. Our initial conditions represent protogalaxies forming within a fossil HII region -- a previously ionized HII region which has not yet had time to cool and recombine. We show that in these regions, H2 is the dominant and most effective coolant, and that it is the amount of H2 formed that controls whether or not the gas can collapse and form stars. For metallicities Z <= 10-3 Zsun, metal line cooling alters the density and temperature evolution of the gas by less than 1% compared to the metal-free case at densities below 1 cm-3 and temperatures above 2000 K. We also find that an external ultraviolet background delays or suppresses the cooling and collapse of the gas regardless of whether it is metal-enriched or not. Finally, we study the dependence of this process on redshift and mass of the dark matter halo. N2 - Sterne sind fundamentale Bestandteile des Kosmos. Sie entstehen im Inneren von turbulenten Molekülwolken, die aus molekularem Wasserstoffgas und Staub bestehen. Durch konvergente Strömungen in der turbulenten Wolke bilden sich lokale Dichtemaxima, die kollabieren, falls die zum Zentrum der Wolke gerichtete Schwerkraft über die nach außen gerichteten Druckkräfte dominiert. Dies ist der Fall, wenn die Masse des Gases einen kritischen Wert überschreitet, der Jeansmasse genannt wird. Die Jeansmasse hängt von der Dichte und der Temperatur des Gases ab und fällt im isothermen Fall mit steigender Dichte stetig ab, so dass während des Kontraktionsprozesses immer kleinere Teilmassen instabil werden. Es kommt zur Fragmentierung der Molekülwolke zu protostellaren Kernen, den direkten Vorläufern der Sterne. In der vorliegenden Arbeit werden die zeitliche Entwicklung des Drehimpulses der protostellaren Kerne und der Einfluss der thermischen Eigenschaften des Gases mit Hilfe von dreidimensionalen hydrodynamischen Simulationen untersucht. Hierbei konzentrieren wir uns auf zwei fundamentale Probleme, die jede Theorie der Sternentstehung lösen muss: das "Drehimpulsproblem" und die Massenverteilung der Sterne (IMF). Die thermischen Eigenschaften des Gases sind nicht nur von Bedeutung für die derzeitige Sternentstehung in beobachtbaren Regionen wie z.B. der Orionnebel oder die ρ-Ophiuchi Molekülwolke, sondern auch für die Entstehung von Sternen der ersten und zweiten Generation im frühen Universum. Wir betrachten die Entwicklung des spezifischen Drehimpulses von protostellaren Kernen und vergleichen unsere Resultate mit beobachteten Werten. Wir finden eine gute Übereinstimmung zwischen den spezifischen Drehimpulsen der protostellaren Kerne in unserem Model und denen der beobachteten Kerne in Molekülwolken. In unseren Simulationen geht der gravitative Kollaps mit einem Verlust an spezifischem Drehimpuls einher. Somit kann das Drehimpulsproblem der Sternentstehung auch ohne Betrachtung der Magnetfelder entschärft werden. Ein weiterer Schwerpunkt der Arbeit ist die Untersuchung des Einflusses der thermodynamischen Eigenschaften des Gases auf die Massenverteilung der Sterne, die aus diesem Gas entstehen. Wir verwenden eine stückweise polytrope Zustandgleichung, die die Temperatur-Dichte-Beziehung genauer beschreibt. Wir zeigen, dass Veränderungen in der Zustandgleichung bei einer bestimmten Dichte einen direkten Einfluss auf die charakteristische Massenskala der Fragmentierung haben und somit den Scheitelpunkt der Sternmassenverteilung in der solaren Umgebung bestimmen. Des Weiteren sind die thermodynamischen Eigenschaften des Gases auch für die Sternentstehung im frühen Universum von Bedeutung. Das primordiale Gas, aus dem die ersten Sterne gebildet wurden, enthält keine Metalle (Elemente schwerer als H oder He), da diese erst durch Kernreaktionen in Sternen gebildet werden. In dieser Arbeit untersuchen wir den Einfluss einer geringen Metallizität auf das Kühlungs- und Kollapsverhalten von Gas, aus welchem die zweite Generation von Sternen entstanden ist. Dieses Gas ist anfänglich heiß und ionisiert und befindet sich in kleinen protogalaktischen Halos aus dunkler Materie. Unsere hydrodynamischen Simulationen, die auch ein adäquates chemisches Netzwerk beinhalten, zeigen, dass die Temperatur- und Dichteentwicklung des Gases während der Anfangsphase des Kollapses durch eine geringe Metallizität im Gas kaum beeinflusst wird. Wir stellen weiterhin fest, dass externe ultraviolette Strahlung den Kühlprozess des Gases ohne Metallizität und des Gases mit geringer Metallizität gleichermaßen verzögert oder sogar verhindert. Außerdem untersuchen wir den Einfluss der Rotverschiebung und der Masse des Halos aus dunkler Materie auf die Kühlung und den Kollaps des Gases. KW - Sternentstehung KW - Astrophysik KW - Turbulenz KW - Hydrodynamisches Modell KW - Gravitationskollaps KW - star formation KW - astrophysics KW - turbulence KW - numerical simulation KW - hydrodynamical model KW - self-gravity Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-7591 ER - TY - THES A1 - Codeco, Marta Sofia Ferreira T1 - Constraining the hydrology at Minas da Panasqueira W-Sn-Cu deposit, Portugal T1 - Bestimmung der Hydrologie der W-Sn-Cu Lagerstätte in der Miene Panasqueira, Portugal T1 - Considerações sobre a hidrologia do depósito de W-Sn-Cu da Panasqueira, Portugal N2 - This dissertation combines field and geochemical observations and analyses with numerical modeling to understand the formation of vein-hosted Sn-W ore in the Panasqueira deposit of Portugal, which is among the ten largest worldwide. The deposit is located above a granite body that is altered by magmatic-hydrothermal fluids in its upper part (greisen). These fluids are thought to be the source of metals, but that was still under debate. The goal of this study is to determine the composition and temperature of hydrothermal fluids at Panasqueira, and with that information to construct a numerical model of the hydrothermal system. The focus is on analysis of the minerals tourmaline and white mica, which formed during mineralization and are widespread throughout the deposit. Tourmaline occurs mainly in alteration zones around mineralized veins and is less abundant in the vein margins. White mica is more widespread. It is abundant in vein margins as well as alteration zones, and also occurs in the granite greisen. The laboratory work involved in-situ microanalysis of major- and trace elements in tourmaline and white mica, and boron-isotope analysis in both minerals by secondary ion mass spectrometry (SIMS). The boron-isotope composition of tourmaline and white mica suggests a magmatic source. Comparison of hydrothermally-altered and unaltered rocks from drill cores shows that the ore metals (W, Sn, Cu, and Zn) and As, F, Li, Rb, and Cs were introduced during the alteration. Most of these elements are also enriched in tourmaline and mica, which confirms their potential value as exploration guides to Sn-W ores elsewhere. The thermal evolution of the hydrothermal system was estimated by B-isotope exchange thermometry and the Ti-in-quartz method. Both methods yielded similar temperatures for the early hydrothermal phase: 430° to 460°C for B-isotopes and 503° ± 24°C for Ti-in-quartz. Mineral pairs from a late fault zone yield significantly lower median temperatures of 250°C. The combined results of thermometry with variations in chemical and B-isotope composition of tourmaline and mica suggest that a similar magmatic-hydrothermal fluid was active at all stages of mineralization. Mineralization in the late stage shows the same B-isotope composition as in the main stage despite a ca. 250°C cooling, which supports a multiple injection model of magmatic-hydrothermal fluids. Two-dimensional numerical simulations of convection in a multiphase NaCl hydrothermal system were conducted: (a) in order to test a new approach (lower dimensional elements) for flow through fractures and faults and (b) in order to identify conditions for horizontal fluid flow as observed in the flat-lying veins at Panasqueira. The results show that fluid flow over an intrusion (heat and fluid source) develops a horizontal component if there is sufficient fracture connectivity. Late, steep fault zones have been identified in the deposit area, which locally contain low-temperature Zn-Pb mineralization. The model results confirm that the presence of subvertical faults with enhanced permeability play a crucial role in the ascent of magmatic fluids to the surface and the recharge of meteoric waters. Finally, our model results suggest that recharge of meteoric fluids and mixing processes may be important at later stages, while flow of magmatic fluids dominate the early stages of the hydrothermal fluid circulation. N2 - In dieser Dissertation werden Feldbeobachtungen und geochemische Analysen mit numerischer Modellierung kombiniert, um die Bildung von Sn-W-Cu- Mineralisation in der Erzlagerstätte Panasqueira in Portugal zu verstehen. Panasqueira gehört zu den 10 größten Sn-W Lagerstätten weltweit, sie befindet sich oberhalb eines Granitkörpers, der im oberen Bereich durch magmatisch-hydrothermale Fluide alteriert ist (Greisenbildung). Es wird postuliert, dass magmatisch Fluide die Quelle für Metalle sind, das wurde aber bislang nicht eindeutig bestätigt. Das Ziel dieser Arbeiten ist es, die Zusammensetzung und Temperatur der hydrothermalen Fluide in Panasqueira zu bestimmen und mit diesen Informationen ein numerisches Modell des hydrothermalen Systems zu erstellen. Der Schwerpunkt liegt auf der Untersuchung von Turmalin und Hellglimmer, welche bei der Mineralisation gebildet wurden und in der gesamten Lagerstätte weit verbreitet sind. Turmalin kommt hauptsächlich in Alterationszonen um die vererzten Quarzgänge vor, sowie weniger häufig im Randbereich der Gänge. Hellglimmer dagegen ist stärker verbreitet. Es kommt sowohl in Quarzgangrändern und Alterationszonen vor als auch im Greisenkörper. Die Laborarbeiten umfassten in-situ Mikroanalytik der Haupt- und Spurenelementgehalte von Turmalin und Hellglimmer sowie die Analyse der Bor-isotopen in beiden Mineralen mittels Sekundärionen-Massenspektrometrie (SIMS). Die Bor-Isotopenzusammensetzung von Turmalin und Hellglimmer deuten auf eine magmatische Quelle hin. Der Vergleich von hydrothermal-überprägten mit unveränderten Gesteinsproben aus Bohrkernen zeigt, dass die Erzmetalle (W, Sn, Cu, Zn) sowie As, F, Li, Rb und Cs während der Alteration hinzugefügt wurden. Die meisten dieser Elemente sind auch in Turmalin und Glimmer angereichert, womit ihre Nützlichkeit als Explorationshilfe für Sn-W-Erze in anderen Gebieten bestätigt wird. Die thermische Entwicklung des Hydrothermalsystems wurde durch B-Isotopenaustausch-Thermometrie sowie durch die Ti-in-Quarz- Methode bestimmt. Beide Methoden ergaben für die frühe Hydrothermalphase ähnliche Temperaturen: 430° - 460°C für B-Isotope und 503° ± 24°C für Ti-in-Quarz. Mineral Paare aus einer späten Verwerfungszone ergaben deutlich niedrigere B-Isotopentemperaturen von durchschnittlich 250°C. Die Kombination der Thermometrie mit den chemischen und B-Isotopenvariationen in Turmalin und Glimmer deutet darauf hin, dass ein ähnliches magmatisch-hydrothermales Fluid in allen Mineralisierungsstufen beteiligt war. Die Mineralisierung im späten Stadium zeigt dieselbe B-Isotopenzusammensetzung wie die Hauptphase trotz der Abkühlung um ca. 250°C, was ein Mehrfachinjektionsmodell des magmatisch-hydrothermalen Fluids unterstützt. Zwei-dimensionale numerische Simulationen der Konvektion in einem mehrphasen NaCl System wurden durchgeführt um: a) eine neue Methode (lower dimensional elements) für hydrothermales Fließen durch Brüche und Störungszonen zu testen und b) die Voraussetzungen für die in Panasqueira dominierende horizontale Fluidbewegung in den flach liegenden Gängen zu identifizieren. Die Ergebnisse zeigen, dass Fluidströmungen immer dann eine starke horizontale Komponente haben wenn ausreichende Bruchverbindungen im Gestein vorhanden sind. Späte, steile Bruchzonen sind in der Umgebung der Lagerstätte identifiziert worden, welche lokal niedrig-temperierte Zn-Pb Mineralisierungen führen. Die Modellergebnisse bestätigen, dass das Vorhandensein subvertikaler Störungszonen mit höherer Permeabilität eine entscheidende Rolle für den Aufstieg magmatischer Fluide zur Oberfläche und das Eindringen von meteorischen Fluiden spielen. Schließlich schlagen unsere Simulationsergebnisse vor, dass das Eindringen meteorischer Fluide und Mischungsprozesse in späteren Phasen der hydrothermalen Zirkulation wichtig sind, während magmatische Fluide in frühen Phasen dominieren. N2 - O estudo apresentado na presente dissertação combina análises e observações campo e geoquímica (e.g. multi-elementar, mineral, isotópica) com modelação numérica por forma a compreender a evolução do sistema hidrotermal da Panasqueira e a formação dos filões sub-horizontais. O jazigo filoniano da Panasqueira, localizado em Portugal, encontra-se entre os dez maiores depósitos do tipo a nível mundial e é o maior produtor de W na União Europeia. O depósito desenvolveu-se a topo de um granito Varisco do tipo S, cuja cúpula se encontra greisenizada devido a circulação de fluidos de carácter magmato-hidrotermal. Pensa-se que estes fluidos sejam a fonte dos metais para a génese do jazigo, contudo esta questão tem constituído matéria de grande debate. Por forma, a compreender a evolução hidrotermal e construir um modelo numérico capaz de simular a hidrologia do sistema hidrotermal da Panasqueira, este trabalho envolveu a determinação da composição e temperatura dos fluidos hidrotermais. Para o efeito, este estudo concentrou-se na caracterização geoquímica e isotópica de turmalina e mica branca, as quais se formaram durante os processos iniciais de alteração hidrotermal e/ou mineralização, ocorrendo em diversos contextos. A turmalina ocorre essencialmente nos halos de alteração hidrotermal que encerram os veios mineralizados e está predominantemente associada aos estádios pré-mineralização. Em contraste, a mica branca ocorre em diversos contextos: greisen, filões (salbandas), e halos de alteração hidrotermal, estando associada quer aos estádios precoces, quer aos estádios principais de mineralização. O trabalho laboratorial envolveu análises in-situ de elementos maiores e traço em turmalina e mica, e análises isotópicas (isótopos de boro) em ambas as fases minerais através de espectrometria de massa por iões secundários (SIMS). As composições isotópicas da turmalina e mica branca sugerem uma fonte magmática para os fluidos hidrotermais. A comparação dos dados de litogeoquímica dos metassedimentos alterados e não alterados mostra que os metais (W, Sn, Cu e Zn), assim como em As, F, Li, Rb e Cs foram introduzidos durante o processo de alteração hidrotermal. Parte substancial destes elementos encontram-se também enriquecidos na mica e turmalina, o que confirma o seu potencial valor como vectores de prospecção mineral para os depósitos de W-Sn. A evolução térmica do sistema hidrotermal da Panasqueira foi estimada utilizando geotermómetros minerais. O geotermómetro do quartzo (Ti-in-quartz) indica temperaturas de 503° ± 24°C para a alteração precoce das rochas encaixantes, o que é consistente com as temperaturas médias de 430° a 460°C} obtidas através da geotermometria isotópica de boro em turmalina e mica branca nas salbandas micáceas. As zonas de falha estudadas através da utilização de pares-minerais indicam temperaturas médias substancialmente mais baixas (250°C). A combinação dos estudos de geotermometria mineral com as variações químicas e isotópicas obtidas para a turmalina e mica sugerem que um fluido magmático-hidrotermal relativamente homogéneo esteve activo durante todos os estádios de mineralização. Durante estádios tardios, a os fluidos mineralizantes possuem as mesmas composições isotópicas obtidas para os estádios principais, embora que registando um arrefecimento de ca. 250°C, o que suporta um modelo dinâmico com múltiplas injecções de fluidos magmático-hidrotermais. Simulações numéricas bidimensionais da convecção num sistema hidrotermal multifásico salino foram conduzidas: (i) para testar uma nova metodologia (“lower dimensional elements”) capaz de traduzir o fluxo de fluidos através de fracturas e falhas e, (ii) para identificar as condições do fluxo horizontal observado nos filões sub-horizontais da Panasqueira. Os resultados mostram que o escoamento dos fluidos em associação com uma intrusão (fonte de calor e fluidos) desenvolve uma componente horizontal, desde que haja conectividade suficiente. Falhas tardias inclinadas identificadas na área contem localmente mineralização de Zn e Pb de baixa temperatura. Os resultados dos modelos numéricos confirmam que a presença de falhas sub-verticais de permeabilidade acrescida tem um papel crucial na ascensão de fluidos magmáticos até a superfície e na infiltração de águas meteóricas. Por fim, os resultados das simulações sugerem que a infiltração de águas meteóricas e processos de mistura de fluidos possam ser importantes durantes os estádios tardios, enquanto os fluidos de carácter magmático dominam os estádios iniciais da circulação hidrotermal dos fluidos. KW - Panasqueira KW - tourmaline KW - muscovite KW - Boron isotopes KW - tungsten-tin deposits KW - magmatic-hydrothermal systems KW - LA-ICP-MS KW - SIMS KW - numerical simulation KW - fluid flow KW - fracture-controlled KW - alteration geochemistry KW - Panasqueira KW - Turmalin KW - Muscovit KW - Bor-isotopen KW - Wofram-Zinn Lagerstätte KW - magmatisch-hydrothermale Systeme KW - LA-ICP-MS KW - SIMS KW - numerische Modellierung KW - Fluid-strömungen KW - strukturelle Kontrolle KW - Alterationsgeochemie KW - Panasqueira KW - turmalina KW - muscovite KW - Isótopos de Boro KW - depósitos de volfrâmio-estanho KW - sistemas magmático-hidrotermais KW - LA-ICP-MS KW - SIMS KW - simulações numéricas KW - fluxo de fluidos KW - controlo estrutural KW - geoquímica da alteração hidrotermal Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-429752 ER -