Assessing the saltwater remediation potential of a three-dimensional,heterogeneous, coastal aquifer system

This paper evaluates the remediation potential of a salinized coastal aquifer by utilizing a scenario simulation. Therefore, the numerical model OpenGeoSys is first validated against analytical and experimental data to represent transient groundwater level development and variable density saline intrusion. Afterwards, a regional scale model with a three-dimensional, heterogeneous hydrogeology is calibrated for a transient state and used to simulate a best-case scenario. Water balances are evaluated in both the transient calibration and scenario run. Visualization techniques help to assess the complex model output providing valuable insight in the occurring density-driven flow processes. Furthermore, modeling and visualization results give information on the time scale for remediation activities and, due to limitations in data quality and quantity reveal potential for model improvement.     

Visualisation strategies for environmental modelling data

We present a framework that allows users to apply a number of strategies to view and modify a wide range of environmental data sets for the modelling of natural phenomena. These data sets can be concurrently visualised to find inconsistencies or artefacts. This ensures at an early stage that models set up for the simulation of hydrological or thermal processes will not give implausible results due to complications based on input data. A number of generally applicable visualisation techniques are provided by our framework to help researchers detect potential problems. We also propose a number of mapping algorithms for the integration of multiple data sets to resolve some of the most common issues. Techniques for the presentation of input- and modelling data in combination with simulation results are proposed with respect to the benefits of visualisation of environmental data within specialised environments. The complete workflow from input data to presentation is demonstrated based on a case study in Central Germany. We identify typical problems, propose approaches for a suitable data integration for this case study and compare results of the original and modified data sets.     

Modeling and simulation of two-phase two-component flow with disappearing nonwetting phase

Carbon capture and storage is a recently discussed new technology, aimed at allowing an ongoing use of fossil fuels while preventing the produced CO₂ to be released to the atmosphere. CCS can be modeled with two components (water and CO₂) in two phases (liquid and CO₂). To simulate the process, a multiphase flow equation with equilibrium phase exchange is used. One of the big problems arising in two-phase two-component flow simulations is the disappearance of the nonwetting phase, which leads to a degeneration of the equations satisfied by the saturation. A standard choice of primary variables, which is the pressure of one phase and the saturation of the other phase, cannot be applied here. We developed a new approach using the pressure of the nonwetting phase and the capillary pressure as primary variables. One important advantage of this approach is the fact that we have only one set of primary variables that can be used for the biphasic as well as the monophasic case. We implemented this new choice of primary variables in the DUNE simulation framework and present numerical results for some test cases.     

Visual data exploration for hydrological analysis

Hydrological research projects for integrated water resources management such as the IWAS initiative often accumulate large amounts of heterogeneous data from different sources. Given the number of partners taking part in such projects, surveying and accessing the available data sets, as well as searching for a defined subset, becomes increasingly difficult. We propose an integrated approach for a system combining visual data management and numerical simulation which allows to survey and select data sets based on keywords such as a region of interest or given indicators. An adequate 3D visualisation of such subsets helps to convey information and significantly supports the assessment of relations between different types of data. Furthermore, the interface between the visual data management system and finite element codes allows for the straightforward integration of information into the numerical simulation process and the subsequent visualisation of simulation results in a geographical context. We demonstrate typical workflows for integration and processing within the system based on data from the IWAS model region in Saudi Arabia and the TERENO Bode Observatory in the Harz Mountains in Germany. In addition, we present examples for data import and export based on established standard file formats.     

Thermo-hydro-mechanical modeling of carbon dioxide injection for enhanced gas-recovery (CO2-EGR): a benchmarking study for code comparison

The objective of this paper was to investigate the THM-coupled responses of the storage formation and caprock, induced by gas production, CO₂-EGR (enhanced gas recovery), and CO₂-storage. A generic 3D planer model (20,000?×?3,000?×?100?m, consisting of 1,200?m overburden, 100?m caprock, 200?m gas reservoir, and 1,500?m base rock) is adopted for the simulation process using the integrated code TOUGH2/EOS7C-FLAC3D and the multi-purpose simulator OpenGeoSys. Both simulators agree that the CO₂-EGR phase under a balanced injection rate (31,500?tons/year) will cause almost no change in the reservoir pressure. The gas recovery rate increases 1.4?% in the 5-year CO₂-EGR phase, and a better EGR effect could be achieved by increasing the distance between injection and production wells (e.g., 5.83?% for 5?km distance, instead of 1.2?km in this study). Under the considered conditions there is no evidence of plastic deformation and both reservoir and caprock behave elastically at all operation stages. The stress path could be predicted analytically and the results show that the isotropic and extensional stress regime will switch to the compressional stress regime, when the pore pressure rises to a specific level. Both simulators agree regarding modification of the reservoir stress state. With further CO₂-injection tension failure in reservoir could occur, but shear failure will never happen under these conditions. Using TOUGH-FLAC, a scenario case is also analyzed with the assumption that the reservoir is naturally fractured. The specific analysis shows that the maximal storage pressure is 13.6?MPa which is determined by the penetration criterion of the caprock.     

Sedimentology and geochemistry of an exceptionally preserved last interglacial sapropel S5 in the Levantine Basin (Mediterranean Sea)

A combined study of lithological, geochemical and physical sediment properties is reported from a completely laminated S5 sapropel, recovered in three gravity cores (M40-4 SL67, M51-3 SL103, M51-3 SL104) from the Pliny Trench region of the eastern Mediterranean. The thickness of the studied sapropel S5 varies between 85 and 91 cm and tops most S5-sapropels in the Mediterranean. Based on optical features like color and thickness of laminae, the sapropels were subdivided into thirteen distinct lithostratigraphic zones. These zones, as well as the finer layering pattern within them, could be followed exactly among the three cores, indicating that the processes responsible for this variation acted at least on a regional scale. The sapropel sediment is characterized by exceptionally high porosity, which is strongly correlated with Si/Ca. This relationship implies that the sapropel is in essence an organic-matter rich diatomite and its exceptional thickness can be explained by preservation of diatoms forming a loosely packed sediment fabric. Compared to other S5 sapropels, the preservation of diatoms has apparently led to a twofold increase in the thickness of the sapropel layer. Relative abundances of 10 elements were determined at ultra-high resolution (0.2 mm) by XRF-scanner over the complete length of each sapropel including several cm of enclosing marl. An analysis of the chemical data indicates that the lowermost 13 cm of the sapropel is chemically more similar to the underlying marl and that the sediment chemistry shows different signals at different scales. The strongest pattern is the contrast between the sapropel and the surrounding marl, which is accentuated in elements indicative for redox conditions as well as terrigenous sediment input and productivity. Within the sapropel, a mm- to cm-scale layering is observed. The abundances of many elements are systematically linked to the pattern of these layers, indicating a common origin, related to productivity and/or terrigenous sediment and/or redox conditions. This pattern indicates a link to a regional climatic process, making the S5 sapropel horizon in M40-4 SL67, M51-3 SL103 and M51-3 SL104 a potential high-resolution archive of climatic variability during the last interglacial in the Mediterranean Sea and its adjacent landmasses.