TY - JOUR A1 - De Lucia, Marco A1 - Kühn, Michael T1 - DecTree v1.0-chemistry speedup in reactive transport simulations BT - purely data-driven and physics-based surrogates JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - The computational costs associated with coupled reactive transport simulations are mostly due to the chemical subsystem: replacing it with a pre-trained statistical surrogate is a promising strategy to achieve decisive speedups at the price of small accuracy losses and thus to extend the scale of problems which can be handled. We introduce a hierarchical coupling scheme in which "full-physics" equation-based geochemical simulations are partially replaced by surrogates. Errors in mass balance resulting from multivariate surrogate predictions effectively assess the accuracy of multivariate regressions at runtime: inaccurate surrogate predictions are rejected and the more expensive equation-based simulations are run instead. Gradient boosting regressors such as XGBoost, not requiring data standardization and being able to handle Tweedie distributions, proved to be a suitable emulator. Finally, we devise a surrogate approach based on geochemical knowledge, which overcomes the issue of robustness when encountering previously unseen data and which can serve as a basis for further development of hybrid physics-AI modelling. Y1 - 2021 U6 - https://doi.org/10.5194/gmd-14-4713-2021 SN - 1991-959X SN - 1991-9603 VL - 14 IS - 7 SP - 4713 EP - 4730 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Hennig, Theresa A1 - Kühn, Michael T1 - Surrogate model for multi-component diffusion of Uranium through Opalinus Clay on the host rock scale JF - Applied Sciences : open access journal N2 - Multi-component (MC) diffusion simulations enable a process based and more precise approach to calculate transport and sorption compared to the commonly used single-component (SC) models following Fick's law. The MC approach takes into account the interaction of chemical species in the porewater with the diffuse double layer (DDL) adhering clay mineral surfaces. We studied the shaly, sandy and carbonate-rich facies of the Opalinus Clay. High clay contents dominate diffusion and sorption of uranium. The MC simulations show shorter diffusion lengths than the SC models due to anion exclusion from the DDL. This hampers diffusion of the predominant species CaUO2(CO3)32-. On the one side, species concentrations and ionic strengths of the porewater and on the other side surface charge of the clay minerals control the composition and behaviour of the DDL. For some instances, it amplifies the diffusion of uranium. We developed a workflow to transfer computationally intensive MC simulations to SC models via calibrated effective diffusion and distribution coefficients. Simulations for one million years depict maximum uranium diffusion lengths between 10 m and 35 m. With respect to the minimum requirement of a thickness of 100 m, the Opalinus Clay seems to be a suitable host rock for nuclear waste repositories. KW - facies KW - uranium speciation KW - sorption KW - reactive transport KW - heterogeneity KW - PHREEQC KW - Mont Terri KW - repository far-field Y1 - 2021 U6 - https://doi.org/10.3390/app11020786 SN - 2076-3417 VL - 11 IS - 2 PB - MDPI CY - Basel ER - TY - JOUR A1 - Tranter, Morgan Alan A1 - De Lucia, Marco A1 - Kühn, Michael T1 - Numerical investigation of barite scaling kinetics in fractures JF - Geothermics : an international journal of geothermal research and its applications N2 - Barite stands out as one of the most ubiquitous scaling agents in deep geothermal systems, responsible for irreversible efficiency loss. Due to complex parameter interplay, it is imperative to utilise numerical simulations to investigate temporal and spatial precipitation effects. A one-dimensional reactive transport model is set up with heterogeneous nucleation and crystal growth kinetics. In line with geothermal systems in the North German Basin, the following parameters are considered in a sensitivity analysis: temperature (25 to 150 degrees C), pore pressure (10 to 50 MPa), fracture aperture (10(-4) to 10(-2) m), flow velocity (10(-3) to 10(0) m s(-1)), molar volume (50.3 to 55.6 cm(3) mol(-1)), contact angle for heterogeneous nucleation (0 degrees to 180 degrees), interfacial tension (0.07 to 0.134 J m(-2)), salinity (0.1 to 1.5 mol kgw(-1) NaCl), pH (5 to 7), and supersaturation ratio (1 to 30). Nucleation and consequently crystal growth can only begin if the threshold supersaturation is exceeded, therefore contact angle and interfacial tension are the most sensitive in terms of precipitation kinetics. If nucleation has occurred, crystal growth becomes the dominant process, which is mainly controlled by fracture aperture. Results show that fracture sealing takes place within months (median 33 days) and the affected range can be on the order of tens of metres (median 10 m). The presented models suggest that barite scaling must be recognised as a serious threat if the supersaturation threshold is exceeded, in which case, large fracture apertures could help to minimise kinetic rates. The models further are of use for adjusting the fluid injection temperature. KW - Geothermal KW - Fracture sealing KW - Reactive transport KW - Geochemical modelling KW - Crystal nucleation KW - PHREEQC Y1 - 2021 SN - 0375-6505 SN - 1879-3576 VL - 91 PB - Elsevier CY - Amsterdam [u.a.] ER - TY - JOUR A1 - Wetzel, Maria A1 - Kempka, Thomas A1 - Kühn, Michael T1 - Diagenetic trends of synthetic reservoir sandstone properties assessed by digital rock physics JF - Minerals N2 - Quantifying interactions and dependencies among geometric, hydraulic and mechanical properties of reservoir sandstones is of particular importance for the exploration and utilisation of the geological subsurface and can be assessed by synthetic sandstones comprising the microstructural complexity of natural rocks. In the present study, three highly resolved samples of the Fontainebleau, Berea and Bentheim sandstones are generated by means of a process-based approach, which combines the gravity-driven deposition of irregularly shaped grains and their diagenetic cementation by three different schemes. The resulting evolution in porosity, permeability and rock stiffness is examined and compared to the respective micro-computer tomographic (micro-CT) scans. The grain contact-preferential scheme implies a progressive clogging of small throats and consequently produces considerably less connected and stiffer samples than the two other schemes. By contrast, uniform quartz overgrowth continuously alters the pore space and leads to the lowest elastic properties. The proposed stress-dependent cementation scheme combines both approaches of contact-cement and quartz overgrowth, resulting in granulometric, hydraulic and elastic properties equivalent to those of the respective micro-CT scans, where bulk moduli slightly deviate by 0.8%, 4.9% and 2.5% for the Fontainebleau, Berea and Bentheim sandstone, respectively. The synthetic samples can be further altered to examine the impact of mineral dissolution or precipitation as well as fracturing on various petrophysical correlations, which is of particular relevance for numerous aspects of a sustainable subsurface utilisation. KW - digital core reconstruction KW - micro-CT scan KW - pore-scale KW - cementation KW - permeability-porosity relationship KW - elastic rock properties KW - numerical KW - simulation Y1 - 2021 U6 - https://doi.org/10.3390/min11020151 SN - 2075-163X VL - 11 IS - 2 PB - MDPI CY - Basel ER - TY - JOUR A1 - Steding, Svenja A1 - Kempka, Thomas A1 - Kühn, Michael T1 - How insoluble inclusions and intersecting layers affect the leaching process within potash seams JF - Applied Sciences : open access journal N2 - Potash seams are a valuable resource containing several economically interesting, but also highly soluble minerals. In the presence of water, uncontrolled leaching can occur, endangering subsurface mining operations. In the present study, the influence of insoluble inclusions and intersecting layers on leaching zone evolution was examined by means of a reactive transport model. For that purpose, a scenario analysis was carried out, considering different rock distributions within a carnallite-bearing potash seam. The results show that reaction-dominated systems are not affected by heterogeneities at all, whereas transport-dominated systems exhibit a faster advance in homogeneous rock compositions. In return, the ratio of permeated rock in vertical direction is higher in heterogeneous systems. Literature data indicate that most natural potash systems are transport-dominated. Accordingly, insoluble inclusions and intersecting layers can usually be seen as beneficial with regard to reducing hazard potential as long as the mechanical stability of leaching zones is maintained. Thereby, the distribution of insoluble areas is of minor impact unless an inclined, intersecting layer occurs that accelerates leaching zone growth in one direction. Moreover, it is found that the saturation dependency of dissolution rates increases the growth rate in the long term, and therefore must be considered in risk assessments. KW - salt dissolution KW - reactive transport KW - heterogeneity KW - density-driven KW - convection KW - PHREEQC KW - porous media Y1 - 2021 U6 - https://doi.org/10.3390/app11199314 SN - 2076-3417 VL - 11 IS - 19 PB - MDPI CY - Basel ER - TY - JOUR A1 - Hennig, Theresa A1 - Kühn, Michael T1 - Potential uranium migration within the geochemical gradient of the opalinus clay system at the Mont Terri JF - Minerals N2 - Transport properties of potential host rocks for nuclear waste disposal are typically determined in laboratory or in-situ experiments under geochemically controlled and constant conditions. Such a homogeneous assumption is no longer applicable on the host rock scale as can be seen from the pore water profiles of the potential host rock Opalinus Clay at Mont Terri (Switzerland). The embedding aquifers are the hydro-geological boundaries, that established gradients in the 210 m thick low permeable section through diffusive exchange over millions of years. Present-day pore water profiles were confirmed by a data-driven as well as by a conceptual scenario. Based on the modelled profiles, the influence of the geochemical gradient on uranium migration was quantified by comparing the distances after one million years with results of common homogeneous models. Considering the heterogeneous system, uranium migrated up to 24 m farther through the formation depending on the source term position within the gradient and on the partial pressure of carbon dioxide pCO2 of the system. Migration lengths were almost equal for single- and multicomponent diffusion. Differences can predominantly be attributed to changes in the sorption capacity, whereby pCO2 governs how strong uranium migration is affected by the geochemical gradient. Thus, the governing parameters for uranium migration in the Opalinus Clay can be ordered in descending priority: pCO2, geochemical gradients, mineralogical heterogeneity.

KW - PHREEQC KW - reactive transport KW - sorption KW - diffusion KW - repository far-field KW - hydro-geological system Y1 - 2021 U6 - https://doi.org/10.3390/min11101087 SN - 2075-163X VL - 11 IS - 10 PB - MDPI CY - Basel ER - TY - JOUR A1 - Steding, Svenja A1 - Kempka, Thomas A1 - Zirkler, Axel A1 - Kühn, Michael T1 - Spatial and temporal evolution of leaching zones within potash seams reproduced by reactive transport simulations JF - Water / Molecular Diversity Preservation International (MDPI) N2 - Leaching zones within potash seams generally represent a significant risk to subsurface mining operations and the construction of technical caverns in salt rocks, but their temporal and spatial formation has been investigated only rudimentarily to date. To the knowledge of the authors, current reactive transport simulation implementations are not capable to address hydraulic-chemical interactions within potash salt. For this reason, a reactive transport model has been developed and complemented by an innovative approach to calculate the interchange of minerals and solution at the water-rock interface. Using this model, a scenario analysis was carried out based on a carnallite-bearing potash seam. The results show that the evolution of leaching zones depends on the mineral composition and dissolution rate of the original salt rock, and that the formation can be classified by the dimensionless parameters of Peclet (Pe) and Damkohler (Da). For Pe > 2 and Da > 1, a funnel-shaped leaching zone is formed, otherwise the dissolution front is planar. Additionally, Da > 1 results in the formation of a sylvinitic zone and a flow barrier. Most scenarios represent hybrid forms of these cases. The simulated shapes and mineralogies are confirmed by literature data and can be used to assess the hazard potential. KW - carnallite KW - water rock interaction KW - density-driven flow KW - PHREEQC KW - Pitzer KW - equations Y1 - 2021 U6 - https://doi.org/10.3390/w13020168 SN - 2073-4441 VL - 13 IS - 2 PB - Molecular Diversity Preservation International CY - Basel ER - TY - JOUR A1 - Tranter, Morgan Alan A1 - De Lucia, Marco A1 - Kühn, Michael T1 - Barite scaling potential modelled for fractured-porous geothermal reservoirs JF - Minerals N2 - Barite scalings are a common cause of permanent formation damage to deep geothermal reservoirs. Well injectivity can be impaired because the ooling of saline fluids reduces the solubility of barite, and the continuous re-injection of supersaturated fluids forces barite to precipitate in the host rock. Stimulated reservoirs in the Upper Rhine Graben often have multiple relevant flow paths in the porous matrix and fracture zones, sometimes spanning multiple stratigraphical units to achieve the economically necessary injectivity. While the influence of barite scaling on injectivity has been investigated for purely porous media, the role of fractures within reservoirs consisting of both fractured and porous sections is still not well understood. Here, we present hydro-chemical simulations of a dual-layer geothermal reservoir to study the long-term impact of barite scale formation on well injectivity. Our results show that, compared to purely porous reservoirs, fractured porous reservoirs have a significantly reduced scaling risk by up to 50%, depending on the flow rate ratio of fractures. Injectivity loss is doubled, however, if the amount of active fractures is increased by one order of magnitude, while the mean fracture aperture is decreased, provided the fractured aquifer dictates the injection rate. We conclude that fractured, and especially hydraulically stimulated, reservoirs are generally less affected by barite scaling and that large, but few, fractures are favourable. We present a scaling score for fractured-porous reservoirs, which is composed of easily derivable quantities such as the radial equilibrium length and precipitation potential. This score is suggested for use approximating the scaling potential and its impact on injectivity of a fractured-porous reservoir for geothermal exploitation. KW - reactive transport KW - radial flow KW - geothermal energy KW - injectivity KW - phreeqc KW - formation damage Y1 - 2021 U6 - https://doi.org/10.3390/min11111198 SN - 2075-163X VL - 11 IS - 11 PB - MDPI CY - Basel ER - TY - JOUR A1 - De Lucia, Marco A1 - Kühn, Michael A1 - Lindemann, Alexander A1 - Lübke, Max A1 - Schnor, Bettina T1 - POET (v0.1): speedup of many-core parallel reactive transport simulations with fast DHT lookups JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - Coupled reactive transport simulations are extremely demanding in terms of required computational power, which hampers their application and leads to coarsened and oversimplified domains. The chemical sub-process represents the major bottleneck: its acceleration is an urgent challenge which gathers increasing interdisciplinary interest along with pressing requirements for subsurface utilization such as spent nuclear fuel storage, geothermal energy and CO2 storage. In this context we developed POET (POtsdam rEactive Transport), a research parallel reactive transport simulator integrating algorithmic improvements which decisively speed up coupled simulations. In particular, POET is designed with a master/worker architecture, which ensures computational efficiency in both multicore and cluster compute environments. POET does not rely on contiguous grid partitions for the parallelization of chemistry but forms work packages composed of grid cells distant from each other. Such scattering prevents particularly expensive geochemical simulations, usually concentrated in the vicinity of a reactive front, from generating load imbalance between the available CPUs (central processing units), as is often the case with classical partitions. Furthermore, POET leverages an original implementation of the distributed hash table (DHT) mechanism to cache the results of geochemical simulations for further reuse in subsequent time steps during the coupled simulation. The caching is hence particularly advantageous for initially chemically homogeneous simulations and for smooth reaction fronts. We tune the rounding employed in the DHT on a 2D benchmark to validate the caching approach, and we evaluate the performance gain of POET's master/worker architecture and the DHT speedup on a 3D benchmark comprising around 650 000 grid elements. The runtime for 200 coupling iterations, corresponding to 960 simulation days, reduced from about 24 h on 11 workers to 29 min on 719 workers. Activating the DHT reduces the runtime further to 2 h and 8 min respectively. Only with these kinds of reduced hardware requirements and computational costs is it possible to realistically perform the longterm complex reactive transport simulations, as well as perform the uncertainty analyses required by pressing societal challenges connected with subsurface utilization. Y1 - 2021 U6 - https://doi.org/10.5194/gmd-14-7391-2021 SN - 1991-959X SN - 1991-9603 VL - 14 IS - 12 SP - 7391 EP - 7409 PB - Copernicus CY - Göttingen ER -