14 results
Search Results
Now showing 1 - 10 of 14
- ItemA 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, ThomasWe 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.
- ItemLarge-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, ThomasThe 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.
- ItemClassification 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 PeresRecent 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.
- ItemFinite-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.
- ItemA 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.
- ItemX-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, StephanIn 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.
- ItemHigh-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.
- ItemPore-scale tomography and imaging: applications, techniques and recommended practice(Göttingen : Copernicus Publ., 2016) Halisch, Matthias; Steeb, Holger; Henkel, Steven; Krawczyk, Charlotte M.[No abstract available]
- ItemElectric resistivity and seismic refraction tomography: a challenging joint underwater survey at Äspö Hard Rock Laboratory(Göttingen : Copernicus Publ., 2017) Ronczka, Mathias; Hellman, Kristofer; Günther, Thomas; Wisén, Roger; Dahlin, TorleifTunnelling below water passages is a challenging task in terms of planning, pre-investigation and construction. Fracture zones in the underlying bedrock lead to low rock quality and thus reduced stability. For natural reasons, they tend to be more frequent at water passages. Ground investigations that provide information on the subsurface are necessary prior to the construction phase, but these can be logistically difficult. Geophysics can help close the gaps between local point information by producing subsurface images. An approach that combines seismic refraction tomography and electrical resistivity tomography has been tested at the Äspö Hard Rock Laboratory (HRL). The aim was to detect fracture zones in a well-known but logistically challenging area from a measuring perspective. The presented surveys cover a water passage along part of a tunnel that connects surface facilities with an underground test laboratory. The tunnel is approximately 100 m below and 20 m east of the survey line and gives evidence for one major and several minor fracture zones. The geological and general test site conditions, e.g. with strong power line noise from the nearby nuclear power plant, are challenging for geophysical measurements. Co-located positions for seismic and ERT sensors and source positions are used on the 450 m underwater section of the 700 m profile. Because of a large transition zone that appeared in the ERT result and the missing coverage of the seismic data, fracture zones at the southern and northern parts of the underwater passage cannot be detected by separated inversion. Synthetic studies show that significant three-dimensional (3-D) artefacts occur in the ERT model that even exceed the positioning errors of underwater electrodes. The model coverage is closely connected to the resolution and can be used to display the model uncertainty by introducing thresholds to fade-out regions of medium and low resolution. A structural coupling cooperative inversion approach is able to image the northern fracture zone successfully. In addition, previously unknown sedimentary deposits with a significantly large thickness are detected in the otherwise unusually well-documented geological environment. The results significantly improve the imaging of some geologic features, which would have been undetected or misinterpreted otherwise, and combines the images by means of cluster analysis into a conceptual subsurface model.
- ItemDistinct 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, TorstenThe 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.