Search Results

Now showing 1 - 8 of 8
  • Item
    Modelling the Surface Heat Flow Distribution in the Area of Brandenburg (Northern Germany)
    (Amsterdam [u.a.] : Elsevier, 2013) Cacace, Mauro; Scheck-Wenderoth, Magdalena; Noack, Vera; Cherubini, Yvonne; Schellschmidt, Rüdiger; Kühn, Michael; Juhlin, Christopher; Held, Hermann; Bruckman, Viktor; Tambach, Tim; Kempka, Thomas
    A lithosphere scale geological model has been used to determine the surface heat flow component due to conductive heat transport for the area of Brandenburg. The modelling results have been constrained by a direct comparison with available heat flow measurements. The calculated heat flow captures the regional trend in the surface heat flow distribution which can be related to existing thermal conductivity variations between the different sedimentary units. An additional advective component due to topography induced regional flow and focused flow within major fault zones should be considered to explain the spatial variation observed in the surface heat flow.
  • Item
    Effectiveness of acidizing geothermal wells in the South German Molasse Basin
    (Göttingen : Copernicus Publ., 2013) Schumacher, S.; Schulz, R.
    In Germany, many hydro-geothermal plants have been constructed in recent years, primarily in the region of Munich. As the host formation here mainly consists of carbonates, nearly all recently drilled wells have been acidized in order to improve the well yield. In this study, the effectiveness of these acid treatments is analyzed with respect to the amount of acid used and the number of acid treatments carried out per well. The results show that the first acid treatment has the largest effect, while subsequent acidizing improves the well only marginally. Data also indicate that continued acidizing can lead to degradation of the well. These findings may not only be important for geothermal installations in Germany but also for projects, for example, in Austria, France or China where geothermal energy is produced from carbonate formations as well.
  • Item
    Luminescence dating of volcanogenic outburst flood sediments from Aso volcano and tephric loess deposits, southwest Japan
    (Warsaw : De Gruyter, 2013) Tsukamoto, Sumiko; Kataoka, Kyoko; Miyabuchi, Yasuo
    Luminescence dating has been applied to volcanogenic outburst flood sediments (Takuma gravel bed) from Aso volcano, Japan, and tephric loess deposits overlying the gravel bed. The poly-mineral fine grains (4-11 μm) from loess deposits were measured with pulsed optically stimulated luminescence (pulsed OSL) and post-IR infrared stimulated luminescence (pIRIR) methods, whereas the Takuma gravel bed containing no quartz, was measured with IRSL and pIRIR methods using sand sized (150-200 μm) plagioclase. The loess deposits date back at least to ∼50 ka by consistent IRSL, pIRIR and pulsed OSL ages from the lowermost part of the loess deposits from one section. The ages obtained from the bottom part of the other loess section are not consistent each other. However, we consider that the pIRIR age (72±6 ka) which showed negligible anomalous fading is most reliable, and regard as a preliminary minimum age of the Takuma gravel bed. The equivalent doses (De) for the plagioclase from the Takuma gravel bed have a narrow distribution and the weighted mean of the three samples yield an age of 89±3 ka. This age is in agreement with the last caldera-forming eruption of Aso volcano (∼87 ka) and it is likely that the pIRIR signal has not been bleached before the deposi-tion. IRSL dating without applying pIRIR using small aliquots was also conducted, however, the IRSL signal shows no clear evidence of an additional bleaching during the event of outburst flood from the caldera lake. © 2013 Silesian University of Technology, Gliwice, Poland. All rights reserved.
  • 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.
  • 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.
  • Item
    Three-dimensional monitoring of soil water content in a maize field using Electrical Resistivity Tomography
    (Munich : EGU, 2013) Beff, L.; Günther, T.; Vandoorne, B.; Couvreur, V.; Javaux, M.
    A good understanding of the soil water content (SWC) distribution at the field scale is essential to improve the management of water, soil and crops. Recent studies proved that Electrical Resistivity Tomography (ERT) opens interesting perspectives in the determination of the SWC distribution in 3 dimensions (3-D). This study was conducted (i) to check and validate how ERT is able to monitor SWC distribution in a maize field during the late growing season; and (ii) to investigate how maize plants and rainfall affect the dynamics of SWC distribution. Time Domain Reflectometry (TDR) measurements were used to validate ERT-inverted SWC values. Evolution of water mass balance was also calculated to check whether ERT was capable of giving a reliable estimate of soil water stock evolution. It is observed that ERT was able to give the same average SWC as TDR (R2 = 0.98). In addition, ERT gives better estimates of the water stock than TDR thanks to its higher spatial resolution. The high resolution of ERT measurements also allows for the discrimination of SWC heterogeneities. The SWC distribution showed that alternation of maize rows and inter-rows was the main influencing factor of the SWC distribution. The drying patterns were linked to the root profiles, with drier zones under the maize rows. During short periods, with negligible rainfall, the SWC decrease took place mainly in the two upper soil horizons and in the inter-row area. In contrast, rainfall increased the SWC mostly under the maize rows and in the upper soil layer. Nevertheless, the total amount of rainfall during the growing season was not sufficient to modify the SWC patterns induced by the maize rows. During the experimental time, there was hardly any SWC redistribution from maize rows to inter-rows. Yet, lateral redistribution from inter-rows to maize rows induced by potential gradient generates SWC decrease in the inter-row area and in the deeper soil horizons.
  • Item
    Seismic imaging of sandbox experiments – laboratory hardware setup and first reflection seismic sections
    (Göttingen : Copernicus Publ., 2013) Krawczyk, C.M.; Buddensiek, M.-L.; Oncken, O.; Kukowski, N.
    With the study and technical development introduced here, we combine analogue sandbox simulation techniques with seismic physical modelling of sandbox models. For that purpose, we designed and developed a new mini-seismic facility for laboratory use, comprising a seismic tank, a PC-driven control unit, a positioning system, and piezoelectric transducers used here for the first time in an array mode. To assess the possibilities and limits of seismic imaging of small-scale structures in sandbox models, different geometry setups were tested in the first 2-D experiments that also tested the proper functioning of the device and studied the seismo-elastic properties of the granular media used. Simple two-layer models of different materials and layer thicknesses as well as a more complex model comprising channels and shear zones were tested using different acquisition geometries and signal properties. We suggest using well sorted and well rounded grains with little surface roughness (glass beads). Source receiver-offsets less than 14 cm for imaging structures as small as 2.0–1.5 mm size have proven feasible. This is the best compromise between wide beam and high energy output, and is applicable with a consistent waveform. Resolution of the interfaces of layers of granular materials depends on the interface preparation rather than on the material itself. Flat grading of interfaces and powder coverage yields the clearest interface reflections. Finally, sandbox seismic sections provide images of high quality showing constant thickness layers as well as predefined channel structures and indications of the fault traces from shear zones. Since these were artificially introduced in our test models, they can be regarded as zones of disturbance rather than tectonic shear zones characterized by decompaction. The multiple-offset surveying introduced here, improves the quality with respect to S / N ratio and source signature even more; the maximum depth penetration in glass-bead layers thereby amounts to 5 cm. Thus, the presented mini-seismic device is already able to resolve structures within simple models of saturated porous media, so that multiple-offset seismic imaging of shallow sandbox models, that are structurally evolving, is generally feasible.
  • 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.