Niklas Linde

Improved hydrogeophysical characterization using joint inversion of cross-hole electrical resistance and ground-penetrating radar traveltime data

Appropriate regularizations of geophysical inverse problems and joint inversion of different data types improve geophysical models and increase their usefulness in hydrogeological studies. We have developed an efficient method to calculate stochastic regularization operators for given geostatistical models. The method, which combines circulant embedding and the diagonalization theorem of circulant matrices, is applicable for stationary geostatistical models when the grid discretization, in each spatial direction, is uniform in the volume of interest. We also used a structural approach to jointly invert cross-hole electrical resistance and ground-penetrating radar traveltime data in three dimensions. The two models are coupled by assuming, at all points, that the cross product of the gradients of the two models is zero. No petrophysical relationship between electrical conductivity and relative permittivity is assumed but is instead obtained as a by-product of the inversion. The approach has been applied to data collected in a U.K. sandstone aquifer in order to improve characterization of the vadose zone hydrostratigraphy. By analyzing scatterplots of electrical conductivity versus relative permittivity together with petrophysical models a zonation could be obtained with corresponding estimates of the electrical formation factor, the water content, and the effective grain radius of the sediments. The approach provides greater insight into the hydrogeological characteristics of the subsurface than by using conventional geophysical inversion methods.

Streaming current generation in two‐phase flow conditions

Self-potential (SP) signals that are generated under two-phase flow conditions could be used to study vadose zone dynamics and to monitor petroleum production. These streaming-potentials may also act as an error source in SP monitoring of vulcanological activity and in magnetotelluric studies. We propose a two-phase flow SP theory that predicts streaming currents as a function of the pore water velocity, the excess of charge in the pore water, and the porosity. The source currents that create the SP signals are given by the divergence of the streaming currents, and contributions are likely to be located at infiltration fronts, at the water table, or at geological boundaries. Our theory was implemented in a hydrogeological modeling code to calculate the SP distribution during primary drainage. Forward and inverse modeling of a well-calibrated 1D drainage experiment suggest that our theory can predict streaming potentials in the vadose zone.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petrophysical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine theparameters of a Gaussian mixture model that defines zonal geometries ...

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform cross-hole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4–10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Critical steps for the continuing advancement of hydrogeophysics

Special hydrogeophysics issues published by hydrology and geophysics journals, special sessions and workshops at conferences, and an increasing number of short courses demonstrate the growing interest in the use of geophysics for hydrologic investigations. The formation of the hydrogeophysics technical subcommittee of AGUs Hydrology section adds further evidence of the recognized significance of this growing interdisciplinary field. Given the clear value of nondestructive and nonintrusive imaging for subsurface investigations, we believe the advances in the adoption of existing geophysical methods, the development of novel methods, and the merging of geophysical and other data made in hydrogeophysics could be applied to a wide range of geological, environmental, and engineering applications.

3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water

The hydrogeological properties and responses of a productive aquifer in northeastern Switzerland are investigated. For this purpose, 3D crosshole electrical resistivity tomography (ERT) is used to define the main lithological structures within the aquifer (through static inversion) and to monitor the water infiltration from an adjacent river. During precipitation events and subsequent river flooding, the river water resistivity increases. As a consequence, the electrical characteristics of the infiltrating water can be used as a natural tracer to delineate preferential flow paths and flow velocities. The focus is primarily on the experiment installation, data collection strategy, and the structural characterization of the site and a brief overview of the ERT monitoring results. The monitoring system comprises 18 boreholes each equipped with 10 electrodes straddling the entire thickness of the gravel aquifer. A multichannel resistivity system programmed to cycle through various four-point electrode configurations of the 180 electrodes in a rolling sequence allows for the measurement of approximately 15,500 apparent resistivity values every 7 h on a continuous basis. The 3D static ERT inversion of data acquired under stable hydrological conditions provides a base model for future time-lapse inversion studies and the means to investigate the resolving capability of our acquisition scheme. In particular, it enables definition of the main lithological structures within the aquifer. The final ERT static model delineates a relatively high-resistivity, low-porosity, intermediate-depth layer throughout the investigated aquifer volume that is consistent with results from well logging and seismic and radar tomography models. The next step will be to define and implement an appropriate time-lapse ERT inversion scheme using the river water as a natural tracer. The main challenge will be to separate the superposed time-varying effects of water table height, temperature, and salinity variations associated with the infiltrating water.

Joint inversion of crosshole radar and seismic traveltimes acquired at the South Oyster Bacterial Transport Site

The structural approach to joint inversion, entailing common boundaries or gradients, offers a flexible and effective way to invert diverse types of surface-based and/or crosshole geophysical data. The cross-gradients function has been introduced as a means to construct models in which spatial changes in two distinct physical-property models are parallel or antiparallel. Inversion methods that use such structural constraints also provide estimates of nonlinear and nonunique field-scale relationships between model parameters. Here, we jointly invert crosshole radar and seismic traveltimes for structurally similar models using an iterative nonlinear traveltime tomography algorithm. Application of the inversion scheme to synthetic data demonstrates that it better resolves lithologic boundaries than the individual inversions alone. Tests of the scheme on GPR and seismic data acquired within a shallow aquifer illustrate that the resultant models have improved correlations with flowmeter data in comparison with models based on individual inversions. The highest correlation with the flowmeter data is obtained when the joint inversion is combined with a stochastic regularization operator and the vertical integral scale is estimated from the flowmeter data. Point-spread functions show that the most significant resolution improvements offered by the joint inversion are in the horizontal direction.

Inferring transport characteristics in a fractured rock aquifer by combining single‐hole ground‐penetrating radar reflection monitoring and tracer test data

[1] Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Detection and localization of hydromechanical disturbances in a sandbox using the self‐potential method

Four sandbox experiments were performed to understand the self-potential response to hydro-mechanical disturbances in a water-infiltrated controlled sandbox. In the first two experiments, ∼0.5 mL of water was abruptly injected through a small capillary at a depth of 15 cm using a syringe impacted by a hammer stroke. In the second series of experiments, ∼0.5 mL of pore water was quickly pumped out of the tank, at the same depth, using a syringe. In both type of experiments, the resulting self-potential signals were measured using 32 sintered Ag/AgCl medical electrodes. In two experiments, these electrodes were located 3 cm below the top surface of the tank. In two other experiments, they were placed along a vertical section crossing the position of the capillary. These electrodes were connected to a voltmeter with a sensitivity of 0.1 μV and an acquisition frequency of 1.024 kHz. The injected/pumped volumes of water produced hydro-mechanical disturbances in the sandbox. In turn, these disturbances generated dipolar electrical anomalies of electrokinetic nature with an amplitude of few microvolts. The source function is the product of the dipolar Greens function by a source intensity function that depends solely on the product of the streaming potential coupling coefficient of the sand to the pore fluid overpressure with respect to the hydrostatic pressure. Numerical modeling using a finite element code was performed to solve the coupled hydro-mechanical problem and to determine the distribution of the resulting self-potential during the course of these experiments. We use 2D and 3D algorithms based on the cross-correlation method and wavelet analysis of potential fields to show that the source was a vertical dipole. These methods were also used to localize the position of the source of the hydromechanical disturbance from the self-potential signals recorded at the top surface of the tank. The position of the source agrees with the position of the inlet/outlet of the capillary showing the usefulness of these methods for application to active volcanoes.

Local earthquake (LE) tomography with joint inversion for P‐ and S‐wave velocities using structural constraints

Local earthquake (LE) tomography with joint inversion for P- and S-wave velocities using structural constraints