TY - JOUR A1 - Hetenyi, Gyorgy A1 - Molinari, Irene A1 - Clinton, John A1 - Bokelmann, Gotz A1 - Bondar, Istvan A1 - Crawford, Wayne C. A1 - Dessa, Jean-Xavier A1 - Doubre, Cecile A1 - Friederich, Wolfgang A1 - Fuchs, Florian A1 - Giardini, Domenico A1 - Graczer, Zoltan A1 - Handy, Mark R. A1 - Herak, Marijan A1 - Jia, Yan A1 - Kissling, Edi A1 - Kopp, Heidrun A1 - Korn, Michael A1 - Margheriti, Lucia A1 - Meier, Thomas A1 - Mucciarelli, Marco A1 - Paul, Anne A1 - Pesaresi, Damiano A1 - Piromallo, Claudia A1 - Plenefisch, Thomas A1 - Plomerova, Jaroslava A1 - Ritter, Joachim A1 - Rumpker, Georg A1 - Sipka, Vesna A1 - Spallarossa, Daniele A1 - Thomas, Christine A1 - Tilmann, Frederik A1 - Wassermann, Joachim A1 - Weber, Michael A1 - Weber, Zoltan A1 - Wesztergom, Viktor A1 - Zivcic, Mladen A1 - Abreu, Rafael A1 - Allegretti, Ivo A1 - Apoloner, Maria-Theresia A1 - Aubert, Coralie A1 - Besancon, Simon A1 - de Berc, Maxime Bes A1 - Brunel, Didier A1 - Capello, Marco A1 - Carman, Martina A1 - Cavaliere, Adriano A1 - Cheze, Jerome A1 - Chiarabba, Claudio A1 - Cougoulat, Glenn A1 - Cristiano, Luigia A1 - Czifra, Tibor A1 - Danesi, Stefania A1 - Daniel, Romuald A1 - Dannowski, Anke A1 - Dasovic, Iva A1 - Deschamps, Anne A1 - Egdorf, Sven A1 - Fiket, Tomislav A1 - Fischer, Kasper A1 - Funke, Sigward A1 - Govoni, Aladino A1 - Groschl, Gidera A1 - Heimers, Stefan A1 - Heit, Ben A1 - Herak, Davorka A1 - Huber, Johann A1 - Jaric, Dejan A1 - Jedlicka, Petr A1 - Jund, Helene A1 - Klingen, Stefan A1 - Klotz, Bernhard A1 - Kolinsky, Petr A1 - Kotek, Josef A1 - Kuhne, Lothar A1 - Kuk, Kreso A1 - Lange, Dietrich A1 - Loos, Jurgen A1 - Lovati, Sara A1 - Malengros, Deny A1 - Maron, Christophe A1 - Martin, Xavier A1 - Massa, Marco A1 - Mazzarini, Francesco A1 - Metral, Laurent A1 - Moretti, Milena A1 - Munzarova, Helena A1 - Nardi, Anna A1 - Pahor, Jurij A1 - Pequegnat, Catherine A1 - Petersen, Florian A1 - Piccinini, Davide A1 - Pondrelli, Silvia A1 - Prevolnik, Snjezan A1 - Racine, Roman A1 - Regnier, Marc A1 - Reiss, Miriam A1 - Salimbeni, Simone A1 - Santulin, Marco A1 - Scherer, Werner A1 - Schippkus, Sven A1 - Schulte-Kortnack, Detlef A1 - Solarino, Stefano A1 - Spieker, Kathrin A1 - Stipcevic, Josip A1 - Strollo, Angelo A1 - Sule, Balint A1 - Szanyi, Gyongyver A1 - Szucs, Eszter A1 - Thorwart, Martin A1 - Ueding, Stefan A1 - Vallocchia, Massimiliano A1 - Vecsey, Ludek A1 - Voigt, Rene A1 - Weidle, Christian A1 - Weyland, Gauthier A1 - Wiemer, Stefan A1 - Wolf, Felix A1 - Wolyniec, David A1 - Zieke, Thomas T1 - The AlpArray seismic network BT - a large-scale european experiment to image the alpine orogen JF - Surveys in Geophysics N2 - The AlpArray programme is a multinational, European consortium to advance our understanding of orogenesis and its relationship to mantle dynamics, plate reorganizations, surface processes and seismic hazard in the Alps-Apennines-Carpathians-Dinarides orogenic system. The AlpArray Seismic Network has been deployed with contributions from 36 institutions from 11 countries to map physical properties of the lithosphere and asthenosphere in 3D and thus to obtain new, high-resolution geophysical images of structures from the surface down to the base of the mantle transition zone. With over 600 broadband stations operated for 2 years, this seismic experiment is one of the largest simultaneously operated seismological networks in the academic domain, employing hexagonal coverage with station spacing at less than 52 km. This dense and regularly spaced experiment is made possible by the coordinated coeval deployment of temporary stations from numerous national pools, including ocean-bottom seismometers, which were funded by different national agencies. They combine with permanent networks, which also required the cooperation of many different operators. Together these stations ultimately fill coverage gaps. Following a short overview of previous large-scale seismological experiments in the Alpine region, we here present the goals, construction, deployment, characteristics and data management of the AlpArray Seismic Network, which will provide data that is expected to be unprecedented in quality to image the complex Alpine mountains at depth. KW - Seismology KW - Alps KW - Seismic network KW - Geodynamics KW - Seismic imaging KW - Mountain building Y1 - 2018 U6 - https://doi.org/10.1007/s10712-018-9472-4 SN - 0169-3298 SN - 1573-0956 VL - 39 IS - 5 SP - 1009 EP - 1033 PB - Springer CY - Dordrecht ER - 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 - 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 - Tranter, Morgan Alan A1 - De Lucia, Marco A1 - Wolfgramm, Markus A1 - Kühn, Michael T1 - Barite scale formation and injectivity loss models for geothermal systems JF - Water N2 - Barite scales in geothermal installations are a highly unwanted effect of circulating deep saline fluids. They build up in the reservoir if supersaturated fluids are re-injected, leading to irreversible loss of injectivity. A model is presented for calculating the total expected barite precipitation. To determine the related injectivity decline over time, the spatial precipitation distribution in the subsurface near the injection well is assessed by modelling barite growth kinetics in a radially diverging Darcy flow domain. Flow and reservoir properties as well as fluid chemistry are chosen to represent reservoirs subject to geothermal exploration located in the North German Basin (NGB) and the Upper Rhine Graben (URG) in Germany. Fluids encountered at similar depths are hotter in the URG, while they are more saline in the NGB. The associated scaling amount normalised to flow rate is similar for both regions. The predicted injectivity decline after 10 years, on the other hand, is far greater for the NGB (64%) compared to the URG (24%), due to the temperature- and salinity-dependent precipitation rate. The systems in the NGB are at higher risk. Finally, a lightweight score is developed for approximating the injectivity loss using the Damkohler number, flow rate and total barite scaling potential. This formula can be easily applied to geothermal installations without running complex reactive transport simulations. KW - reactive transport KW - radial flow KW - geothermal energy KW - scaling KW - phreeqc KW - formation damage Y1 - 2020 U6 - https://doi.org/10.3390/w12113078 SN - 2073-4441 VL - 12 IS - 11 PB - MDPI 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 - 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 -