Refine
Year of publication
Document Type
- Article (49)
- Postprint (3)
- Monograph/Edited Volume (1)
- Conference Proceeding (1)
- Other (1)
- Part of Periodical (1)
Keywords
- Ground-penetrating radar (4)
- Electromagnetics (2)
- Global inversion (2)
- Ground penetrating radar (2)
- ground-penetrating radar (2)
- preferential flow (2)
- tracer (2)
- 3D modeling (1)
- Applied geophysics (1)
- Aquifer (1)
As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.
The in-phase response collected by portable loop-loop electromagnetic induction (EMI) sensors operating at low and moderate induction numbers (<= 1) is typically used for sensing the magnetic permeability (or susceptibility) of the subsurface. This is due to the fact that the in-phase response contains a small induction fraction and a preponderant induced magnetization fraction. The magnetization fraction follows the magneto-static equations similarly to the magnetic method but with an active magnetic source. The use of an active source offers the possibility to collect data with several loop-loop configurations, which illuminate the subsurface with different sensitivity patterns. Such multiconfiguration soundings thereby allows the imaging of subsurface magnetic permeability/susceptibility variations through an inversion procedure. This method is not affected by the remnant magnetization and theoretically overcomes the classical depth ambiguity generally encountered with passive geomagnetic data. To invert multiconfiguration in-phase data sets, we propose a novel methodology based on a full-grid 3-D multichannel deconvolution (MCD) procedure. This method allows us to invert large data sets (e.g. consisting of more than a hundred thousand of data points) for a dense voxel-based 3-D model of magnetic susceptibility subject to smoothness constraints. In this study, we first present and discuss synthetic examples of our imaging procedure, which aim at simulating realistic conditions. Finally, we demonstrate the applicability of our method to field data collected across an archaeological site in Auvergne (France) to image the foundations of a Gallo-Roman villa built with basalt rock material. Our synthetic and field data examples demonstrate the potential of the proposed inversion procedure offering new and complementary ways to interpret data sets collected with modern EMI instruments.
A comprehensive workflow to analyze ensembles of globally inverted 2D electrical resistivity models
(2022)
Electrical resistivity tomography (ERT) aims at imaging the subsurface resistivity distribution and provides valuable information for different geological, engineering, and hydrological applications. To obtain a subsurface resistivity model from measured apparent resistivities, stochastic or deterministic inversion procedures may be employed. Typically, the inversion of ERT data results in non-unique solutions; i.e., an ensemble of different models explains the measured data equally well. In this study, we perform inference analysis of model ensembles generated using a well-established global inversion approach to assess uncertainties related to the nonuniqueness of the inverse problem. Our interpretation strategy starts by establishing model selection criteria based on different statistical descriptors calculated from the data residuals. Then, we perform cluster analysis considering the inverted resistivity models and the corresponding data residuals. Finally, we evaluate model uncertainties and residual distributions for each cluster. To illustrate the potential of our approach, we use a particle swarm optimization (PSO) algorithm to obtain an ensemble of 2D layer-based resistivity models from a synthetic data example and a field data set collected in Loon-Plage, France. Our strategy performs well for both synthetic and field data and allows us to extract different plausible model scenarios with their associated uncertainties and data residual distributions. Although we demonstrate our workflow using 2D ERT data and a PSObased inversion approach, the proposed strategy is general and can be adapted to analyze model ensembles generated from other kinds of geophysical data and using different global inversion approaches.
Modern mobile devices (i.e. smartphones and tablet computers) are widespread, everyday tools, which are equipped with a variety of sensors including three-axis magnetometers. Here, we investigate the feasibility and the potential of using such mobile devices to mimic geophysical experiments in the classroom in a table-top setup. We focus on magnetic surveying and present a basic setup of a table-top experiment for collecting three-component magnetic data across well-defined source bodies and structures. Our results demonstrate that the quality of the recorded data is sufficient to address a number of important basic concepts in the magnetic method. The shown examples cover the analysis of magnetic data recorded across different kinds of dipole sources, thus illustrating the complexity of magnetic anomalies. In addition, we analyze the horizontal resolution capabilities using a pair of dipole sources placed at different horizontal distances to each other. Furthermore, we demonstrate that magnetic data recorded with a mobile device can even be used to introduce filtering, transformation, and inversion approaches as they are typically used when processing magnetic data sets recorded for real-world field applications. Thus, we conclude that such table-top experiments represent an easy-to-implement experimental procedure (as student exercise or classroom demonstration) and can provide first hands-on experience in the basic principles of magnetic surveying including the fundamentals of data acquisition, analysis and processing, as well as data evaluation and interpretation.
Direct current systems employing a kinematic surveying strategy allow to analyze the electrical resistivity of the subsurface for large areas (i.e., several hectares). Typical applications are found in precision agriculture, archaeological prospecting and soil sciences. With the typical survey setting, the collected data sets are often characterized by a rather high level of noise and a rather coarse lateral sampling compared to data acquired with fixed electrodes. We therefore present an efficient one-dimensional inversion approach in which we put special attention on modeling the effects of noise. We apply this method to data recorded with a five-offset equatorial dipole-dipole system employing rolling electrodes. By performing several synthetic tests with realistic noise levels, we found that the considered five-configuration soundings allow for a reliable imaging of two-layer cases in the uppermost two meters of the subsurface, where the subsurface can be assumed to follow a horizontally layered geometry within 3 m around the system. By analyzing the corresponding sensitivity functions, we also show that the equatorial dipole-dipole array is relatively well suited for a 1D inversion approach compared to standard in-line electrode arrays. To illustrate this aspect, we show that our method can provide results similar to those obtained with a 2D Wenner imaging procedure for data recorded across a well-constrained 2D target. We finally apply our method to a large five-offset data set acquired in an agricultural study. The final pseudo-3D model of electrical resistivity is in accordance with borehole data available for the surveyed area. Our results demonstrate the applicability and the versatility of the presented inversion approach for large-scale data sets as they are typically collected with such rolling electrode systems. (C) 2017 Elsevier B.V. All rights reserved.
We present an algorithm that performs sequentially one-dimensional inversion of subsurface magnetic permeability and electrical conductivity by using multi-configuration electromagnetic induction sensor data. The presented method is based on the conversion of the in-phase and out-of-phase data into effective magnetic permeability and electrical conductivity of the equivalent homogeneous half-space. In the case of small-offset systems, such as portable electromagnetic induction sensors, for which in-phase and out-of-phase data are moderately coupled, the effective half-space magnetic permeability and electrical conductivity can be inverted sequentially within an iterative scheme. We test and evaluate the proposed inversion strategy using synthetic and field examples. First, we apply it to synthetic data for some highly magnetic environments. Then, the method is tested on real field data acquired in a basaltic environment to image a formation of archaeological interest. These examples demonstrate that a joint interpretation of in-phase and out-of-phase data leads to a better characterisation of the subsurface in magnetic environments such as volcanic areas.
In an attempt to map the shallow geometry of the Maleme Fault Zone (North Island, New Zealand) and estimate vertical displacements of selected fault strands, we have collected 2D and 3D georadar data using 100 MHz antennae. The 2D data consisted of three parallel georadar lines recorded perpendicular to the axis of the well-defined graben of the Maleme Fault Zone. These similar to 160 in long lines, which were 7.5 m apart, crossed several fault strands on either side of the graben axis. The processed georadar sections revealed two prominent parallel reflections that originated from the boundaries of Late Pleistocene lacustrine and tephra deposits. Distinct vertical offsets of these reflections allowed us to estimate displacernents at individual fault strands across the entire inner graben. The total displacements represented by these offsets was similar to 10-20% greater than that inferred from geomorphological studies, thus demonstrating the limitations of surface observations for determining cumulative fault movements. The 3D georadar data set, recorded across an area of similar to 70x similar to 20 in to one side of the graben axis, provided key details on individual fault strands. For the 3D visualization of fault-related structures, various spatial attribute analyses based on the cosine of the instantaneous phase proved to be useful
Inversions of an individual geophysical data set can be highly nonunique, and it is generally difficult to determine petrophysical parameters from geophysical data. We show that both issues can be addressed by adopting a statistical multiparameter approach that requires the acquisition, processing, and separate inversion of two or more types of geophysical data. To combine information contained in the physical-property models that result from inverting the individual data sets and to estimate the spatial distribution of petrophysical parameters in regions where they are known at only a few locations. we demonstrate the potential of the fuzzy c-means (FCM) clustering technique. After testing this new approach on synthetic data, we apply it to limited crosshole georadar, crosshole seismic, gamma-log, and slug-test data acquired within a shallow alluvial aquifer. The derived multiparameter model effectively outlines the major sedimentary units observed in numerous boreholes and provides plausible estimates for the spatial distributions of gamma-ray emitters and hydraulic conductivity
Interdisziplinäres Zentrum für Musterdynamik und Angewandte Fernerkundung Workshop vom 9. - 10. Februar 2006
Near-surface geophysical techniques are extensively used in a variety of engineering, environmental, geologic, and hydrologic applications. While many of these applications ask for detailed, quantitative models of selected material properties, geophysical data are increasingly used to estimate such target properties. Typically, this estimation procedure relies on a two-step workflow including (1) the inversion of geophysical data and (2) the petrophysical translation of the inverted parameter models into the target properties. Standard deterministic implementations of such a quantitative interpretation result in a single best-estimate model, often without considering and propagating the uncertainties related to the two steps. We address this problem by using a rather novel, particle-swarm-based global joint strategy for data inversion and by implementing Monte Carlo procedures for petrophysical property estimation. We apply our proposed workflow to crosshole ground-penetrating radar, P-, and S-wave data sets collected at a well-constrained test site for a detailed geotechnical characterization of unconsolidated sands. For joint traveltime inversion, the chosen global approach results in ensembles of acceptable velocity models, which are analyzed to appraise inversion-related uncertainties. Subsequently, the entire ensembles of inverted velocity models are considered to estimate selected petrophysical properties including porosity, bulk density, and elastic moduli via well-established petrophysical relations implemented in a Monte Carlo framework. Our results illustrate the potential benefit of such an advanced interpretation strategy; i.e., the proposed workflow allows to study how uncertainties propagate into the finally estimated property models, while concurrently we are able to study the impact of uncertainties in the used petrophysical relations (e.g., the influence of uncertain, user-specified parameters). We conclude that such statistical approaches for the quantitative interpretation of geophysical data can be easily extended and adapted to other applications and geophysical methods and might be an important step toward increasing the popularity and acceptance of geophysical tools in engineering practice.
Foreword
(2019)
Ground-penetrating radar is widely used to provide highly resolved images of subsurface sedimentary structures, with implications for processes active in the vadose zone. Frequently overlooked among these structures are tunnels excavated by fossorial animals (e.g., moles). We present two repeated ground-penetrating radar surveys performed a year apart in 2016 and 2017. Careful three-dimensional data processing reveals, in each data set, a pattern of elongated structures that are interpreted as a subsurface mole tunnel network. Our data demonstrate the ability of three-dimensional ground-penetrating radar imaging to non-invasively delineate the small animal tunnels (similar to 5 cm diameter) at a higher spatial and geolocation resolution than has previously been achieved. In turn, this makes repeated surveys and, therefore, long-term monitoring possible. Our results offer valuable insight into the understanding of the near-surface and showcase a potential new application for a geophysical method as well as a non-invasive method of ecological surveying.
Ground-penetrating radar (GPR) is an established geophysical tool to explore a wide range of near-surface environments. Today, the use of synthetic GPR data is largely limited to 2D because 3D modeling is computationally more expensive. In fact, only recent developments of modeling tools and powerful hardware allow for a time-efficient computation of extensive 3D data sets. Thus, 3D subsurface models and resulting GPR data sets, which are of great interest to develop and evaluate novel approaches in data analysis and interpretation, have not been made publicly available up to now. <br /> We use a published hydrofacies data set of an aquifer-analog study within fluvio-glacial deposits to infer a realistic 3D porosity model showing heterogeneities at multiple spatial scales. Assuming fresh-water saturated sediments, we generate synthetic 3D GPR data across this model using novel GPU-acceleration included in the open-source software gprMax. We present a numerical approach to examine 3D wave-propagation effects in modeled GPR data. Using the results of this examination study, we conduct a spatial model decomposition to enable a computationally efficient 3D simulation of a typical GPR reflection data set across the entire model surface. We process the resulting GPR data set using a standard 3D structural imaging sequence and compare the results to selected input data to demonstrate the feasibility and potential of the presented modeling studies. We conclude on conceivable applications of our 3D GPR reflection data set and the underlying porosity model, which are both publicly available and, thus, can support future methodological developments in GPR and other near-surface geophysical techniques.
Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.
Ground-penetrating radar (GPR) is a standard geophysical technique used to image near-surface structures in sedimentary environments. In such environments, GPR data acquisition and processing are increasingly following 3D strategies. However, the processed GPR data volumes are typically still interpreted using selected 2D slices and manual concepts such as GPR facies analyses. In seismic volume interpretation, the application of (semi-)automated and reproducible approaches such as 3D attribute analyses as well as the production of attribute-based facies models are common practices today. In contrast, the field of 3D GPR attribute analyses and corresponding facies models is largely untapped. We have developed and applied a workflow to produce 3D attribute-based GPR facies models comprising the dominant sedimentary reflection patterns in a GPR volume, which images complex sandy structures on the dune island of Spiekeroog (Northern Germany). After presenting our field site and details regarding our data acquisition and processing, we calculate and filter 3D texture attributes to generate a database comprising the dominant texture features of our GPR data. Then, we perform a dimensionality reduction of this database to obtain meta texture attributes, which we analyze and integrate using composite imaging and (also considering additional geometric information) fuzzy c-means cluster analysis resulting in a classified GPR facies model. Considering our facies model and a corresponding GPR facies chart, we interpret our GPR data set in terms of near-surface sedimentary units, the corresponding depositional environments, and the recent formation history at our field site. Thus, we demonstrate the potential of our workflow, which represents a novel and clear strategy to perform a more objective and consistent interpretation of 3D GPR data collected across different sedimentary environments.
Hydrogeophysik : Erkundungen und Sicherung der Ressource Wasser : Antrittsvorlesung 2006-06-01
(2006)
Die weltweite Wasserversorgung basiert zu einem überwiegenden Teil auf Grundwasser. Die Erkundung, der Schutz, die nachhaltige Nutzung sowie die eventuelle Sanierung dieser Grundwasserressourcen sind somit global von fundamentalem gesellschaftlichem Interesse. Bei vielen dieser grundwasserbezogenen Fragestellungen ist häufig eine effiziente und detaillierte Charakterisierung des Untergrundes notwendig. Geophysikalische Messverfahren liefern Abbilder der physikalischen Eigenschaften, wie beispielsweise des elektrischen Widerstandes, die wichtige Informationen über den geometrischen und stofflichen Aufbau des verborgenen Untergrundes liefern. In der Vorlesung wird gezeigt, wie die Verfahren der Angewandten Geophysik auf Fragestellungen hinsichtlich der Präsenz, Ausbreitung und Qualität der Ressource Grundwasser eingesetzt werden können. Darüber hinaus werden aktuelle Forschungsthemen und offene Fragen angesprochen.
There are far-reaching conceptual similarities between bi-static surface georadar and post-stack, "zero-offset" seismic reflection data, which is expressed in largely identical processing flows. One important difference is, however, that standard deconvolution algorithms routinely used to enhance the vertical resolution of seismic data are notoriously problematic or even detrimental to the overall signal quality when applied to surface georadar data. We have explored various options for alleviating this problem and have tested them on a geologically well-constrained surface georadar dataset. Standard stochastic and direct deterministic deconvolution approaches proved to be largely unsatisfactory. While least-squares-type deterministic deconvolution showed some promise, the inherent uncertainties involved in estimating the source wavelet introduced some artificial "ringiness". In contrast, we found spectral balancing approaches to be effective, practical and robust means for enhancing the vertical resolution of surface georadar data, particularly, but not exclusively, in the uppermost part of the georadar section, which is notoriously plagued by the interference of the direct air- and groundwaves. For the data considered in this study, it can be argued that band- limited spectral blueing may provide somewhat better results than standard band-limited spectral whitening, particularly in the uppermost part of the section affected by the interference of the air- and groundwaves. Interestingly, this finding is consistent with the fact that the amplitude spectrum resulting from least-squares-type deterministic deconvolution is characterized by a systematic enhancement of higher frequencies at the expense of lower frequencies and hence is blue rather than white. It is also consistent with increasing evidence that spectral "blueness" is a seemingly universal, albeit enigmatic, property of the distribution of reflection coefficients in the Earth. Our results therefore indicate that spectral balancing techniques in general and spectral blueing in particular represent simple, yet effective means of enhancing the vertical resolution of surface georadar data and, in many cases, could turn out to be a preferable alternative to standard deconvolution approaches.
We present cross-hole P- and S-wave seismic experiments that have been performed along a similar to 100 m long transect for the detailed characterization of a contaminated sedimentary site (Bitterfeld research test site, Germany). We invert the corresponding first break arrival times for the P- and S-wave velocity structure and compare two different strategies to interpret these models in terms of pertinent lithological and geotechnical parameter variations. The first (common) approach is based on directly translating the tomographic velocity models into the parameters of interest (e.g., elastic moduli). The second (zonal) approach first reduces the tomographic parameter information to a limited number of characteristic velocity combinations via k-means cluster analysis. Then, for each zone (cluster) further parameters including uncertainties can be estimated. In the presented case study, Our results indicate that the zonal approach provides an effective means for the integrated interpretation of different co-located data.
In many hydrological applications, ground-wave velocity measurements are increasingly used to map and monitor shallow soil water content. In this study, we propose an automated spectral velocity analysis method to determine the direct ground-wave (DGW) velocity from common midpoint (CMP) or multi-offset ground-penetrating radar (GPR) data. The method introduced in this paper is a variation of the well-known spectral velocity analysis for seismic and GPR reflection events where velocity spectra are computed using different coherency measures along hyperbolas following the normal moveout model. Here, the unnormalized cross-correlation is computed between waveforms across data gathers that are corrected with a linear moveout equation using a predefined range of velocities. Peaks in the resulting velocity spectra identify linear events in the GPR data gathers like DGW events and allow for estimating the corresponding velocities. In addition to obtaining a DGW velocity measurement, we propose a robust method to estimate the associated velocity uncertainties based on the width of the peak in the calculated velocity spectrum. Our proposed method is tested on synthetic data examples to evaluate the influence of subsurface velocity, surveying geometry and signal frequency on the accuracy of estimated ground-wave velocities. In addition, we investigate the influence of such velocity uncertainties on subsequent soil water content estimates using an established petrophysical relationship. Furthermore, we apply our approach to analyse field data, which were collected across a test site in Canada to monitor a wide range of seasonal soil moisture variations. A comparison between our spectral velocity estimates and results derived from manually picked ground-wave arrivals shows good agreement, which illustrates that our spectral velocity analysis is a feasible tool to analyse DGW arrivals in multi-offset GPR data gathers in an objective and more automated manner.