Filters

2012 - 20196

More filters

20226
Reset filters

Settings

DDC: 550 × Subject: numerical model × WGL Institute: LIAG ×

Search Results

Now showing 1 - 6 of 6
  • Thumbnail Image
    Item
    Ground-penetrating radar insight into a coastal aquifer: the freshwater lens of Borkum Island
    (Munich : EGU, 2013) Igel, J.; Günther, T.; Kuntzer, M.
    Freshwater lenses, as important resource for drinking water, are sensitive to climate changes and sea level rise. To simulate this impact on the groundwater systems, hydraulic subsurface models have to be designed. Geophysical techniques can provide information for generating realistic models. The aim of our work is to show how ground-penetrating radar (GPR) investigations can contribute to such hydrological simulations. In the pilot area, Borkum island, GPR was used to map the shape of the groundwater table (GWT) and to characterise the aquifer. In total, 20 km of constant offset (CO) profiles were measured with centre frequencies of 80 and 200 MHz. Wave velocities were determined by common midpoint (CMP) measurements and vertical radar profiling (VRP) in a monitoring well. The 80 MHz CO data show a clear reflection at the groundwater table, whereas the reflection is weaker for the 200 MHz data. After correcting the GPR water tables for the capillary rise, they are in good accordance with the pressure heads of the observation wells in the area. In the centre of the island, the groundwater table is found up to 3.5 m above sea level, however it is lower towards the coastline and marshland. Some local depressions are observed in the region of dune valleys and around pumping stations of the local water supplier. GPR also reveals details within the sediments and highly-permeable aeolian sands can be distinguished from less-permeable marine sediments. Further, a silt loam layer below the water table could be mapped on a large area. The reflection characteristics indicates scattered erosion channels in this layer that cause it to be an aquitard with some leakage. GPR provides a high resolution map of the groundwater table and insight into the stratigraphy of the sediments and their hydraulic properties. This is valuable complementary information to the observation of sparsely distributed monitoring wells as input to hydraulic simulation.
  • Thumbnail Image
    Item
    Distinct element geomechanical modelling of the formation of sinkhole clusters within large-scale karstic depressions
    (Göttingen : Copernicus Publ., 2019) Al-Halbouni, Djamil; Holohan, Eoghan P.; Taheri, Abbas; Watson, Robert A.; Polom, Ulrich; Schöpfer, Martin P. J.; Emam, Sacha; Dahm, Torsten
    The 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvala-like depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea.
  • Thumbnail Image
    Item
    Development of a numerical workflow based on μ-CT imaging for the determination of capillary pressure–saturation-specific interfacial area relationship in 2-phase flow pore-scale porous-media systems: a case study on Heletz sandstone
    (Göttingen : Copernicus Publ., 2016) Peche, Aaron; Halisch, Matthias; Bogdan Tatomir, Alexandru; Sauter, Martin
    In this case study, we present the implementation of a finite element method (FEM)-based numerical pore-scale model that is able to track and quantify the propagating fluid–fluid interfacial area on highly complex micro-computed tomography (μ-CT)-obtained geometries. Special focus is drawn to the relationship between reservoir-specific capillary pressure (pc), wetting phase saturation (Sw) and interfacial area (awn). The basis of this approach is high-resolution μ-CT images representing the geometrical characteristics of a georeservoir sample. The successfully validated 2-phase flow model is based on the Navier–Stokes equations, including the surface tension force, in order to consider capillary effects for the computation of flow and the phase-field method for the emulation of a sharp fluid–fluid interface. In combination with specialized software packages, a complex high-resolution modelling domain can be obtained. A numerical workflow based on representative elementary volume (REV)-scale pore-size distributions is introduced. This workflow aims at the successive modification of model and model set-up for simulating, such as a type of 2-phase problem on asymmetric μ-CT-based model domains. The geometrical complexity is gradually increased, starting from idealized pore geometries until complex μ-CT-based pore network domains, whereas all domains represent geostatistics of the REV-scale core sample pore-size distribution. Finally, the model can be applied to a complex μ-CT-based model domain and the pc–Sw–awn relationship can be computed.
  • Thumbnail Image
    Item
    Spectral induced polarization measurements for predicting the hydraulic conductivity in sandy aquifers
    (Munich : EGU, 2013) Attwa, M.; Günther, T.
    Field and laboratory spectral induced polarization (SIP) measurements are integrated to characterize the hydrogeological conditions at the Schillerslage test site in Germany. The phase images are capable of monitoring thin peat layers within the sandy aquifers. However, the field results show limitations of decreasing resolution with depth. In comparison with the field inversion results, the SIP laboratory measurements show a certain shift in SIP response due to different compaction and sorting of the samples. The SIP data are analyzed to derive an empirical relationship for predicting the hydraulic conductivity (K). In particular, two significant but weak correlations between individual real resistivities (ρ') and relaxation times (τ), based on a Debye decomposition (DD) model, with measured K are found for the upper groundwater aquifer. The maximum relaxation time (τmax) and logarithmically weighted average relaxation time (τlw) show a better relation with K values than the median value τ50. A combined power law relation between individual ρ' and τ with K is developed with an expression of A · (ρ')B · (τlw)C, where A, B and C are determined using a least-squares fit between the measured and predicted K. The suggested approach with the calculated coefficients of the first aquifer is applied for the second. Results show good correlation with the measured K indicating that the derived relationship is superior to single phase angle models as Börner or Slater models.
  • Thumbnail Image
    Item
    Numerical modelling of climate change impacts on freshwater lenses on the North Sea Island of Borkum using hydrological and geophysical methods
    (Munich : EGU, 2012) Sulzbacher, H.; Wiederhold, H.; Siemon, B.; Grinat, M.; Igel, J.; Burschil, T.; Günther, T.; Hinsby, K.
    A numerical, density dependent groundwater model is set up for the North Sea Island of Borkum to estimate climate change impacts on coastal aquifers and especially the situation of barrier islands in the Wadden Sea. The database includes information from boreholes, a seismic survey, a helicopter-borne electromagnetic (HEM) survey, monitoring of the freshwater-saltwater boundary by vertical electrode chains in two boreholes, measurements of groundwater table, pumping and slug tests, as well as water samples. Based on a statistical analysis of borehole columns, seismic sections and HEM, a hydrogeological model is set up. The groundwater model is developed using the finite-element programme FEFLOW. The density dependent groundwater model is calibrated on the basis of hydraulic, hydrological and geophysical data, in particular spatial HEM and local monitoring data. Verification runs with the calibrated model show good agreement between measured and computed hydraulic heads. A good agreement is also obtained between measured and computed density or total dissolved solids data for both the entire freshwater lens on a large scale and in the area of the well fields on a small scale. For simulating future changes in this coastal groundwater system until the end of the current century, we use the climate scenario A2, specified by the Intergovernmental Panel on Climate Change and, in particular, the data for the German North Sea coast. Simulation runs show proceeding salinisation with time beneath the well fields of the two waterworks Waterdelle and Ostland. The modelling study shows that the spreading of well fields is an appropriate protection measure against excessive salinisation of the water supply until the end of the current century.
  • Thumbnail Image
    Item
    Two-dimensional numerical investigations on the termination of bilinear flow in fractures
    (Göttingen : Copernicus Publ., 2013) Ortiz R., A.E.; Jung, R.; Renner, J.
    Bilinear flow occurs when fluid is drained from a permeable matrix by producing it through an enclosed fracture of finite conductivity intersecting a well along its axis. The terminology reflects the combination of two approximately linear flow regimes: one in the matrix with flow essentially perpendicular to the fracture, and one along the fracture itself associated with the non-negligible pressure drop in it. We investigated the characteristics, in particular the termination, of bilinear flow by numerical modeling allowing for an examination of the entire flow field without prescribing the flow geometry in the matrix. Fracture storage capacity was neglected relying on previous findings that bilinear flow is associated with a quasi-steady flow in the fracture. Numerical results were generalized by dimensionless presentation. Definition of a dimensionless time that, other than in previous approaches, does not use geometrical parameters of the fracture permitted identifying the dimensionless well pressure for the infinitely long fracture as the master curve for type curves of all fractures with finite length from the beginning of bilinear flow up to fully developed radial flow. In log–log scale the master curve's logarithmic derivative initially follows a 1/4-slope straight line (characteristic for bilinear flow) and gradually bends into a horizontal line (characteristic for radial flow) for long times. During the bilinear flow period, isobars normalized to well pressure propagate with the fourth and second root of time in fracture and matrix, respectively. The width-to-length ratio of the pressure field increases proportional to the fourth root of time during the bilinear period, and starts to deviate from this relation close to the deviation of well pressure and its derivative from their fourth-root-of-time relations. At this time, isobars are already significantly inclined with respect to the fracture. The type curves of finite fractures all deviate counterclockwise from the master curve instead of clockwise or counterclockwise from the 1/4-slope straight line as previously proposed. The counterclockwise deviation from the master curve was identified as the arrival of a normalized isobar reflected at the fracture tip 16 times earlier. Nevertheless, two distinct regimes were found in regard to pressure at the fracture tip when bilinear flow ends. For dimensionless fracture conductivities TD < 1, a significant pressure increase is not observed at the fracture tip until bilinear flow is succeeded by radial flow at a fixed dimensionless time. For TD > 10, the pressure at the fracture tip has reached substantial fractions of the associated change in well pressure when the flow field transforms towards intermittent formation linear flow at times that scale inversely with the fourth power of dimensionless fracture conductivity. Our results suggest that semi-log plots of normalized well pressure provide a means for the determination of hydraulic parameters of fracture and matrix after shorter test duration than for conventional analysis.