Search Results

Now showing 1 - 10 of 19
  • Item
    Rearrangement of stresses in fault zones - detecting major issues of coupled hydraulic–mechanical processes with relevance to geothermal applications
    (Göttingen : Copernicus Publ., 2014) Ziefle, G.
    The South German Molasse Basin provides favourable conditions for geothermal plants. Nevertheless, micro-seismic events occur in the vicinity of the geothermal Unterhaching Gt2 well and seem to be caused by the geothermal plant. The injection and production are located in an existing fault system. The majority of seismic events takes place at a horizontal distance of 500 m or less of the borehole. However, none of the seismic events are located in the injection reservoir but in fact at a significantly greater depth. A deeper process understanding of the interacting thermal–hydraulic–mechanical effects in the vicinity of the well is desired. This article presents a significantly simplified 2-D model, investigating interactions of the stress field in the vicinity of the geothermal well and movements in the fault system. This might be of special interest, as the operation of the geothermal plant might lead to changes in the material and fracture properties on the one hand and in the equilibrium state on the other. A detailed description of the model, as well as various parameter studies, is presented. It can be seen that boundary conditions such as direction of the stress field in relation to the fault system, geometry of the fault system and parameters of the fractures have a significant influence on stresses in the proximity of the geothermal well. A variation in the spatial stress field in some parts of the fault system is to be expected. For the chosen assumptions the dimension of this variation is about 25% of the assumed stresses. Future work on this model might focus on the characteristics of the fault system, as well as on the influence of the coupled thermal–hydraulic–mechanical effects.
  • 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
    A remote-control datalogger for large-scale resistivity surveys and robust processing of its signals using a software lock-in approach
    (Göttingen : Copernicus Publ., 2018) Oppermann, Frank; Günther, Thomas
    We present a new versatile datalogger that can be used for a wide range of possible applications in geosciences. It is adjustable in signal strength and sampling frequency, battery saving and can remotely be controlled over a Global System for Mobile Communication (GSM) connection so that it saves running costs, particularly in monitoring experiments. The internet connection allows for checking functionality, controlling schedules and optimizing pre-amplification. We mainly use it for large-scale electrical resistivity tomography (ERT), where it independently registers voltage time series on three channels, while a square-wave current is injected. For the analysis of this time series we present a new approach that is based on the lock-in (LI) method, mainly known from electronic circuits. The method searches the working point (phase) using three different functions based on a mask signal, and determines the amplitude using a direct current (DC) correlation function. We use synthetic data with different types of noise to compare the new method with existing approaches, i.e. selective stacking and a modified fast Fourier transformation (FFT)-based approach that assumes a 1∕f noise characteristics. All methods give comparable results, but the LI is better than the well-established stacking method. The FFT approach can be even better but only if the noise strictly follows the assumed characteristics. If overshoots are present in the data, which is typical in the field, FFT performs worse even with good data, which is why we conclude that the new LI approach is the most robust solution. This is also proved by a field data set from a long 2-D ERT profile.
  • Item
    Large-scale electrical resistivity tomography in the Cheb Basin (Eger Rift) at an International Continental Drilling Program (ICDP) monitoring site to image fluid-related structures
    (Göttingen : Copernicus Publ., 2019) Nickschick, Tobias; Flechsig, Christina; Mrlina, Jan; Oppermann, Frank; Löbig, Felix; Günther, Thomas
    The Cheb Basin, a region of ongoing swarm earthquake activity in the western Czech Republic, is characterized by intense carbon dioxide degassing along two known fault zones – the N–S-striking Počatky–Plesná fault zone (PPZ) and the NW–SE-striking Mariánské Lázně fault zone (MLF). The fluid pathways for the ascending CO2 of mantle origin are one of the subjects of the International Continental Scientific Drilling Program (ICDP) project “Drilling the Eger Rift” in which several geophysical surveys are currently being carried out in this area to image the topmost hundreds of meters to assess the structural situation, as existing boreholes are not sufficiently deep to characterize it. As electrical resistivity is a sensitive parameter to the presence of conductive rock fractions as liquid fluids, clay minerals, and also metallic components, a large-scale dipole–dipole experiment using a special type of electric resistivity tomography (ERT) was carried out in June 2017 in order to image fluid-relevant structures. We used permanently placed data loggers for voltage measurements in conjunction with moving high-power current sources to generate sufficiently strong signals that could be detected all along the 6.5 km long profile with 100 and 150 m dipole spacings. After extensive processing of time series for voltage and current using a selective stacking approach, the pseudo-section is inverted, which results in a resistivity model that allows for reliable interpretations depths of up than 1000 m. The subsurface resistivity image reveals the deposition and transition of the overlying Neogene Vildštejn and Cypris formations, but it also shows a very conductive basement of phyllites and granites that can be attributed to high salinity or rock alteration by these fluids in the tectonically stressed basement. Distinct, narrow pathways for CO2 ascent are not observed with this kind of setup, which hints at wide degassing structures over several kilometers within the crust instead. We also observed gravity and GPS data along this profile in order to constrain ERT results. A gravity anomaly of ca. −9 mGal marks the deepest part of the Cheb Basin where the ERT profile indicates a large accumulation of conductive rocks, indicating a very deep weathering or alteration of the phyllitic basement due to the ascent of magmatic fluids such as CO2. We propose a conceptual model in which certain lithologic layers act as caps for the ascending fluids based on stratigraphic records and our results from this experiment, providing a basis for future drillings in the area aimed at studying and monitoring fluids.
  • Item
    Classification and quantification of pore shapes in sandstone reservoir rocks with 3-D X-ray micro-computed tomography
    (Göttingen : Copernicus Publ., 2016) Schmitt, Mayka; Halisch, Matthias; Müller, Cornelia; Fernandes, Celso Peres
    Recent years have seen a growing interest in the characterization of the pore morphologies of reservoir rocks and how the spatial organization of pore traits affects the macro behavior of rock–fluid systems. With the availability of 3-D high-resolution imaging, such as x-ray micro-computed tomography (µ-CT), the detailed quantification of particle shapes has been facilitated by progress in computer science. Here, we show how the shapes of irregular rock particles (pores) can be classified and quantified based on binary 3-D images. The methodology requires the measurement of basic 3-D particle descriptors (length, width, and thickness) and a shape classification that involves the similarity of artificial objects, which is based on main pore network detachments and 3-D sample sizes. Two main pore components were identified from the analyzed volumes: pore networks and residual pore ganglia. A watershed algorithm was applied to preserve the pore morphology after separating the main pore networks, which is essential for the pore shape characterization. The results were validated for three sandstones (S1, S2, and S3) from distinct reservoirs, and most of the pore shapes were found to be plate- and cube-like, ranging from 39.49 to 50.94 % and from 58.80 to 45.18 % when the Feret caliper descriptor was investigated in a 10003 voxel volume. Furthermore, this study generalizes a practical way to correlate specific particle shapes, such as rods, blades, cuboids, plates, and cubes to characterize asymmetric particles of any material type with 3-D image analysis.
  • Item
    Finite-difference modelling to evaluate seismic P-wave and shear-wave field data
    (Göttingen : Copernicus Publ., 2015) Burschil, T.; Beilecke, T.; Krawczyk, C.M.
    High-resolution reflection seismic methods are an established non-destructive tool for engineering tasks. In the near surface, shear-wave reflection seismic measurements usually offer a higher spatial resolution in the same effective signal frequency spectrum than P-wave data, but data quality varies more strongly. To discuss the causes of these differences, we investigated a P-wave and a SH-wave seismic reflection profile measured at the same location on the island of Föhr, Germany and applied seismic reflection processing to the field data as well as finite-difference modelling of the seismic wave field. The simulations calculated were adapted to the acquisition field geometry, comprising 2 m receiver distance (1 m for SH wave) and 4 m shot distance along the 1.5 km long P-wave and 800 m long SH-wave profiles. A Ricker wavelet and the use of absorbing frames were first-order model parameters. The petrophysical parameters to populate the structural models down to 400 m depth were taken from borehole data, VSP (vertical seismic profile) measurements and cross-plot relations. The simulation of the P-wave wave-field was based on interpretation of the P-wave depth section that included a priori information from boreholes and airborne electromagnetics. Velocities for 14 layers in the model were derived from the analysis of five nearby VSPs (vP =1600–2300 m s-1). Synthetic shot data were compared with the field data and seismic sections were created. Major features like direct wave and reflections are imaged. We reproduce the mayor reflectors in the depth section of the field data, e.g. a prominent till layer and several deep reflectors. The SH-wave model was adapted accordingly but only led to minor correlation with the field data and produced a higher signal-to-noise ratio. Therefore, we suggest to consider for future simulations additional features like intrinsic damping, thin layering, or a near-surface weathering layer. These may lead to a better understanding of key parameters determining the data quality of near-surface shear-wave seismic measurements.
  • Item
    A numerical sensitivity study of how permeability, porosity, geological structure, and hydraulic gradient control the lifetime of a geothermal reservoir
    (Göttingen : Copernicus Publ., 2019) Bauer, Johanna F.; Krumbholz, Michael; Luijendijk, Elco; Tanner, David C.
    Geothermal energy is an important and sustainable resource that has more potential than is currently utilized. Whether or not a deep geothermal resource can be exploited, mostly depends on, besides temperature, the utilizable reservoir volume over time, which in turn largely depends on petrophysical parameters. We show, using over 1000 (n=1027) 4-D finite-element models of a simple geothermal doublet, that the lifetime of a reservoir is a complex function of its geological parameters, their heterogeneity, and the background hydraulic gradient (BHG). In our models, we test the effects of porosity, permeability, and BHG in an isotropic medium. Furthermore, we simulate the effect of permeability contrast and anisotropy induced by layering, fractures, and a fault. We quantify the lifetime of the reservoir by measuring the time to thermal breakthrough, i.e. how many years pass before the temperature of the produced fluid falls below the 100 ∘C threshold. The results of our sensitivity study attest to the positive effect of high porosity; however, high permeability and BHG can combine to outperform the former. Particular configurations of all the parameters can cause either early thermal breakthrough or extreme longevity of the reservoir. For example, the presence of high-permeability fractures, e.g. in a fault damage zone, can provide initially high yields, but it channels fluid flow and therefore dramatically restricts the exploitable reservoir volume. We demonstrate that the magnitude and orientation of the BHG, provided permeability is sufficiently high, are the prime parameters that affect the lifetime of a reservoir. Our numerical experiments show also that BHGs (low and high) can be outperformed by comparatively small variations in permeability contrast (103) and fracture-induced permeability anisotropy (101) that thus strongly affect the performance of geothermal reservoirs.
  • Item
    X-ray computed tomography investigation of structures in Opalinus Clay fromlarge-scale to small-scale after mechanical testing
    (Göttingen : Copernicus Publ., 2016) Kaufhold, Annette; Halisch, Matthias; Zacher, Gerhard; Kaufhold, Stephan
    In the past years X-ray computed tomography (CT) has became more and more common for geoscientific applications and is used from the µm-scale (e.g. for investigations of microfossils or pore-scale structures) up to the dm-scale (full drill cores or soil columns). In this paper we present results from CT imaging and mineralogical investigations of an Opalinus Clay core on different scales and different regions of interest, emphasizing especially the 3-D evaluation and distribution of cracks and their impact on mechanical testing of such material. Enhanced knowledge of the testing behaviour of the Opalinus Clay is of great interest, especially since this material is considered for a long-term radioactive waste disposal and storage facility in Switzerland. Hence, results are compared regarding the mineral (i.e. phase) contrast resolution, the spatial resolution, and the overall scanning speed. With this extensive interdisciplinary scale-down approach it has been possible to characterize the general fracture propagation in comparison to mineralogical and textural features of the Opalinus Clay. Additionally, and as far as we know, a so-called mylonitic zone, located at an intersect of two main fractures, has been observed for the first time for an experimentally deformed Opalinus sample. The multi-scale results are in good accordance to data from naturally deformed Opalinus Clay samples, which enables us to perform systematical research under controlled laboratory conditions. Accompanying 3-D imaging greatly enhances the capability of data interpretation and assessment of such a material.
  • Item
    Asymmetry of high-velocity lower crust on the South Atlantic rifted margins and implications for the interplay of magmatism and tectonics in continental breakup
    (Göttingen : Copernicus Publ., 2014) Becker, K.; Franke, D.; Trumbull, R.; Schnabel, M.; Heyde, I.; Schreckenberger, B.; Koopmann, H.; Bauer, K.; Jokat, W.; Krawczyk, C.M.
    High-velocity lower crust (HVLC) and seaward-dipping reflector (SDR) sequences are typical features of volcanic rifted margins. However, the nature and origin of HVLC is under discussion. Here we provide a comprehensive analysis of deep crustal structures in the southern segment of the South Atlantic and an assessment of HVLC along the margins. Two new seismic refraction lines off South America fill a gap in the data coverage and together with five existing velocity models allow for a detailed investigation of the lower crustal properties on both margins. An important finding is the major asymmetry in volumes of HVLC on the conjugate margins. The seismic refraction lines across the South African margin reveal cross-sectional areas of HVLC 4 times larger than at the South American margin, a finding that is opposite to the asymmetric distribution of the flood basalts in the Paraná–Etendeka Large Igneous Province. Also, the position of the HVLC with respect to the SDR sequences varies consistently along both margins. Close to the Falkland–Agulhas Fracture Zone in the south, a small body of HVLC is not accompanied by SDRs. In the central portion of both margins, the HVLC is below the inner SDR wedges while in the northern area, closer to the Rio Grande Rise-Walvis Ridge, large volumes of HVLC extend far seaward of the inner SDRs. This challenges the concept of a simple extrusive/intrusive relationship between SDR sequences and HVLC, and it provides evidence for formation of the HVLC at different times during the rifting and breakup process. We suggest that the drastically different HVLC volumes are caused by asymmetric rifting in a simple-shear-dominated extension.
  • Item
    High-resolution shear-wave seismic reflection as a tool to image near-surface subrosion structures – a case study in Bad Frankenhausen, Germany
    (Göttingen : Copernicus Publ., 2016) Wadas, Sonja H.; Polom, Ulrich; Krawczyk, Charlotte M.
    Subrosion is the subsurface leaching of soluble rocks that results in the formation of depression and collapse structures. This global phenomenon is a geohazard in urban areas. To study near-surface subrosion structures, four shear-wave seismic reflection profiles, with a total length of ca. 332 m, were carried out around the famous leaning church tower of Bad Frankenhausen in northern Thuringia, Germany, which shows an inclination of 4.93° from the vertical. Most of the geological underground of Thuringia is characterized by soluble Permian deposits, and the Kyffhäuser Southern Margin Fault is assumed to be a main pathway for water to leach the evaporite. The seismic profiles were acquired with the horizontal micro-vibrator ELVIS, developed at Leibniz Institute for Applied Geophysics (LIAG), and a 72 m long landstreamer equipped with 72 horizontal geophones. The high-resolution seismic sections show subrosion-induced structures to a depth of ca. 100 m and reveal five features associated with the leaching of Permian deposits: (1) lateral and vertical varying reflection patterns caused by strongly heterogeneous strata, (2) discontinuous reflectors, small offsets, and faults, which show the underground is heavily fractured, (3) formation of depression structures in the near-surface, (4) diffractions in the unmigrated seismic sections that indicate increased scattering of the seismic waves, and (5) varying seismic velocities and low-velocity zones that are presumably caused by fractures and upward-migrating cavities. A previously undiscovered southward-dipping listric normal fault was also found, to the north of the church. It probably serves as a pathway for water to leach the Permian formations below the church and causes the tilting of the church tower. This case study shows the potential of horizontal shear-wave seismic reflection to image near-surface subrosion structures in an urban environment with a horizontal resolution of less than 1 m in the uppermost 10–15 m.