High resolution multi-facies realizations of sedimentary reservoir and aquifer analogs

Combined effects of aquifer heterogeneity and subsurface dam on nitrate contamination in coastal aquifers

Subsurface dams are effective for seawater intrusion mitigation, yet they can cause upstream nitrate accumulation. This research examines the interplay between subsurface dam construction and aquifer layering on nitrate pollution in coastal settings, employing numerical models to simulate density-driven flow and reactive transport. The study reveals that while subsurface dams are adept at curbing seawater intrusion, they inadvertently broaden the nitrate accumulation zone, especially when a low-K layer is present. Heterogeneous aquifers see more pronounced nitrate accumulation from subsurface dams. This effect is pronounced as it influences dissolved organic carbon dynamics, with a notable retreat inland correlating with the expansion of the nitrate pollution plume. A critical finding is that controlling seawater intrusion via dam height adjustment within the Effective Damming Region effectively reduces nitrate levels and bolsters freshwater output. However, exceeding the critical threshold-where the dam surpasses the low-K layer's bottom-results in a substantial shift in nitrate concentration, underscoring the need for precise dam height calibration to avoid aggravating nitrate pollution. This study's innovative contribution lies in its quantification of the nuanced effects of subsurface dams in stratified aquifers, providing an empirical basis for dam design that considers the layered complexities of coastal aquifers. The insights offer a valuable framework for managing nitrate contamination, thus informing sustainable coastal groundwater management and protection strategies.

GEOENT: A Toolbox for Calculating Directional Geological Entropy

Geological entropy is based on Shannon information entropy and measures order in the structure of a spatial random variable. Metrics have been defined to quantify geological entropy in multidimensional (2D and 3D) heterogeneous systems, for instance, porous and fractured geological media. This study introduces GEOENT, a toolbox that can efficiently be used to calculate geological entropy metrics for any kind of input-gridded field. Additionally, the definition of geological entropy metrics is updated to consider anisotropy in the structure of the heterogeneous system. Directional entrograms provide more accurate descriptions of spatial order over different Cartesian directions. This study presents the development of the geological entropy metrics, a description of the toolbox, and examples of its applications in different datasets, including 2D and 3D gridded fields, representing a variety of heterogeneous environments at different scales, from pore-scale microtomography (μCT) images to aquifer analogues.

Comprehensive geological dataset describing a crystalline rock mass for hydraulic stimulation experiments

High-resolution 3D geological models are crucial for underground development projects and corresponding numerical simulations with applications in e.g., tunneling, hydrocarbon exploration, geothermal exploitation and mining. Most geological models are based on sparse geological data sampled pointwise or along lines (e.g., boreholes), leading to oversimplified model geometries. In the framework of a hydraulic stimulation experiment in crystalline rock at the Grimsel Test Site, we collected geological data in 15 boreholes using a variety of methods to characterize a decameter-scale rock volume. The experiment aims to identify and understand relevant thermo-hydro-mechanical-seismic coupled rock mass responses during high-pressure fluid injections. Prior to fluid injections, we characterized the rock mass using geological, hydraulic and geophysical prospecting. The combination of methods allowed for compilation of a deterministic 3D geological analog that includes five shear zones, fracture density information and fracture locations. The model may serve as a decameter-scale analog of crystalline basement rocks, which are often targeted for enhanced geothermal systems. In this contribution, we summarize the geological data and the resulting geological interpretation.

Combining Hydraulic Head Analysis with Airborne Electromagnetics to Detect and Map Impermeable Aquifer Boundaries

Impermeable aquifer boundaries affect the flow of groundwater, transport of contaminants, and the drawdown of water levels in response to pumping. Hydraulic methods can detect the presence of such boundaries, but these methods are not suited for mapping complex, 3D geological bodies. Airborne electromagnetic (AEM) methods produce 3D geophysical images of the subsurface at depths relevant to most groundwater investigations. Interpreting a geophysical model requires supporting information, and hydraulic heads offer the most direct means of assessing the hydrostratigraphic function of interpreted geological units. This paper presents three examples of combined hydraulic and AEM analysis of impermeable boundaries in glacial deposits of eastern Nebraska, USA. Impermeable boundaries were detected in a long-term hydrograph from an observation well, a short-duration pumping test, and a water table map. AEM methods, including frequency-domain and time-domain AEM, successfully imaged the impermeable boundaries, providing additional details about the lateral extent of the geological bodies. Hydraulic head analysis can be used to verify the hydrostratigraphic interpretation of AEM, aid in the correlation of boundaries through areas of noisy AEM data, and inform the design of AEM surveys at local to regional scales.

A method to generate small-scale, high-resolution sedimentary bedform architecture models representing realistic geologic facies

Small-scale (mm to m) sedimentary structures (e.g. ripple lamination, cross-bedding) have received a great deal of attention in sedimentary geology. The influence of depositional heterogeneity on subsurface fluid flow is now widely recognized, but incorporating these features in physically-rational bedform models at various scales remains problematic. The current investigation expands the capability of an existing set of open-source codes, allowing generation of high-resolution 3D bedform architecture models. The implemented modifications enable the generation of 3D digital models consisting of laminae and matrix (binary field) with characteristic depositional architecture. The binary model is then populated with petrophysical properties using a textural approach for additional analysis such as statistical characterization, property upscaling, and single and multiphase fluid flow simulation. One example binary model with corresponding threshold capillary pressure field and the scripts used to generate them are provided, but the approach can be used to generate dozens of previously documented common facies models and a variety of property assignments. An application using the example model is presented simulating buoyant fluid (CO₂) migration and resulting saturation distribution.