Coupling hydrogeological with surface runoff model in a Poltva case study in Western Ukraine

This paper presents the hydrological coupling of the software framework OpenGeoSys (OGS) with the EPA Storm Water Management Model (SWMM). Conceptual models include the Saint Venant equation for river flow, the 2D Darcy equations for confined and unconfined groundwater flow, a two-way hydrological coupling flux in a compartment coupling approach (conductance concept), and Lagrangian particles for solute transport in the river course. A SWMM river–OGS aquifer inter-compartment coupling flux is examined for discharging groundwater in a systematic parameter sensitivity analysis. The parameter study involves a small perturbation (first-order) sensitivity analysis and is performed for a synthetic test example base-by-base through a comprehensive range of aquifer parametrizations. Through parametrization, the test cases enables to determine the leakance parameter for simulating streambed clogging and non-ocillatory river-aquifer water exchange rates with the sequential (partitioned) coupling scheme. The implementation is further tested with a hypothetical but realistic 1D river–2D aquifer model of the Poltva catchment, where discharging groundwater in the upland area affects the river–aquifer coupling fluxes downstream in the river course (propagating feedbacks). Groundwater contribution in the moving river water is numerically determined with Lagrangian particles. A numerical experiment demonstrates that the integrated river–aquifer model is a serviceable and realistic constituent in a complete compartment model of the Poltva catchment.     

Numerical analysis of thermal effects during carbon dioxide injection with enhanced gas recovery: a theoretical case study for the Altmark gas field

Prediction about reservoir temperature change during carbon dioxide injection requires consideration of all, often subtle, thermal effects. In particular, Joule?CThomson cooling (JTC) and the viscous heat dissipation (VHD) effect are factors that cause flowing fluid temperature to differ from the static formation temperature. In this work, warm-back behavior (thermal recovery after injection completed), as well as JTC and VHD effects, at a multi-layered depleted gas reservoir are demonstrated numerically. OpenGeoSys (OGS) is able to solve coupled partial differential equations for pressure, temperature and mole-fraction of each component of the mixture with a combination of monolithic and staggered approaches. The Galerkin finite element approach is adapted for space discretization of governing equations, whereas for temporal discretization, a generalized implicit single-step scheme is used. For numerical modeling of warm-back behavior, we chose a simplified test case of carbon dioxide injection. This test case is numerically solved by using OGS and FeFlow simulators independently. OGS differs from FeFlow in the capability of representing multi-componential effects on warm-back behavior. We verify both code results by showing the close comparison of shut-in temperature profiles along the injection well. As the JTC cooling rate is inversely proportional to the volumetric heat capacity of the solid matrix, the injection layers are cooled faster as compared to the non-injection layers. The shut-in temperature profiles are showing a significant change in reservoir temperature; hence it is important to account for thermal effects in injection monitoring.     

Modeling, parameterization and evaluation of monitoring methods for CO2 storage in deep saline formations: the CO2-MoPa project

Capture and geological sequestration of CO₂ from large industrial sources is considered a measure for reducing anthropogenic emissions of CO₂ and thus mitigating climate change. One of the main storage options proposed are deep saline formations, as they provide the largest potential storage capacities among the geologic options. A thorough assessment of this type of storage site therefore is required. The CO₂-MoPa project aims at contributing to the dimensioning of CO₂ storage projects and to evaluating monitoring methods for CO₂ injection by an integrated approach. For this, virtual, but realistic test sites are designed geometrically and fully parameterized. Numerical process models are developed and then used to simulate the effects of a CO₂ injection into the virtual test sites. Because the parameterization of the virtual sites is known completely, investigation as well as monitoring methods can be closely examined and evaluated by comparing the virtual monitoring result with the simulation. To this end, the monitoring or investigation method is also simulated, and the (virtual) measurements are recorded and evaluated like real data. Application to a synthetic site typical for the north German basin showed that pressure response has to be evaluated taking into account the layered structure of the storage system. Microgravimetric measurements are found to be promising for detecting the CO₂ phase distribution. A combination of seismic and geoelectric measurements can be used to constrain the CO₂ phase distribution for the anticline system used in the synthetic site.     

CLEAN: project overview on CO2 large-scale enhanced gas recovery in the Altmark natural gas field (Germany)

The joint research project CLEAN was conducted in the years 2008?C2011 by a German research and development (R&D) alliance of 16 partners from science and industry. The project was set-up as pilot project to investigate the processes relevant to enhanced gas recovery (EGR) by the injection of CO₂ into a subfield of the almost depleted Altmark natural gas field. Despite the setback that permission for active injection was not issued by the mining authority during the period of the project, important results fostering the understanding of processes linked with EGR were achieved. Work carried out led to a comprehensive evaluation of the EGR potential of the Altmark field and the Altensalzwedel subfield in particular. The calculated safety margins emphasize that technical well integrity of the 12 examined boreholes is given for EGR without a need for any further intervention. The laboratory and field tests confirm that the Altensalzwedel subfield is suitable for the injection of 100,000?t of CO₂. Numerical simulations provide sound predictions for the efficiency and safety of the EGR technology based on the CO₂ injection. The development and testing of different monitoring techniques facilitate an improved surveying of CO₂ storage sites in general. The CLEAN results provide the technological, logistic and conceptual prerequisites for implementing a CO₂-based EGR project in the Altmark and provide a benchmark for similar projects in the world.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

Chemical and structural relations of epitaxial xenotime and zircon substratum in sedimentary and hydrothermal environments: a TEM study

Xenotime overgrowths on detrital zircon in siliciclastic sediments have been reported in numerous studies. However, in natural samples, solid solution of zircon and xenotime is limited to near-end-member compositions. In order to characterize the interface region between both minerals and to draw inferences on the growth mechanisms of authigenic xenotime, we studied xenotime overgrowths on detrital zircon grains from two Phanerozoic sandstone samples with contrasting post-depositional histories. In one sample, the small (≤10 μm), pyramidal xenotime overgrowths are of diagenetic origin and grew without major discontinuity on the detrital zircon grain. The second sample shows up to >50-μm-wide, porous and inclusion-rich, hydrothermal xenotime overgrowths on detrital zircon, whereas the transition zone between both minerals is accompanied by large pore volume. Chemical compositions of the xenotime precipitates from the two samples differ particularly in Y, REE, Th and Sc concentrations, whereas high MREE availability in the diagenetic sample and the presence of Sc in the hydrothermal sample, respectively, appear to have promoted xenotime growth. Transmission electron microscopy on electron-transparent foils cut from the interface region shows that both the diagenetic xenotime and the hydrothermal xenotime are crystalline and grew in optical and crystallographic continuity to their detrital zircon substrata. Only a narrow transition zone (≤90 nm—diagenetic sample, 200–300 nm—hydrothermal sample) between zircon and xenotime is in part made up of nanometre-scale crystalline domains that are slightly distorted and may have formed from dissolution–re-precipitation processes at the zircon rim along with precipitation from the respective fluid.     

A data exploration framework for validation and setup of hydrological models

Over the course of hydrological research projects often a large number of heterogeneous data sets are acquired from sources as diverse as boreholes, gauging stations or satellite imagery. This data then need to be integrated into models for the simulation of hydrological processes. We propose a framework for exploration of geoscientific data and visually guided preparation of such models. Data sets from a large number of sources can be imported, combined and validated to avoid potential problems due to artefacts or inconsistencies between data sets in a subsequent simulation. Boundary conditions and domain discretisations for surface and subsurface models can be created and tested regarding criteria indicating possible numerical instabilities. All data sets including simulation results can be integrated into a user-controlled 3D scene and aspects of the data can be enhanced using a number of established visualisation techniques including thresholding and user-defined transfer functions. We present the application of this framework for the preparation of a model for simulation of groundwater flow in a river catchment in southwest Germany investigated in the scope of the WESS project.