Vertical radar profiling: Combined analysis of traveltimes, amplitudes, and reflections

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. H23-H35 ◽  
Author(s):  
Jens Tronicke ◽  
Göran Hamann

Vertical radar profiling (VRP) is a single-borehole geophysical technique, in which the receiver antenna is located within a borehole and the transmitter antenna is placed at one or various offsets from the borehole. Today, VRP surveying is primarily used to derive 1D velocity models by inverting the arrival times of direct waves. Using field data collected at a well-constrained test site in Germany, we evaluated a VRP workflow relying on the analysis of direct-arrival traveltimes and amplitudes as well as on imaging reflection events. To invert our VRP traveltime data, we used a global inversion strategy resulting in an ensemble of acceptable velocity models, and thus, it allowed us to appraise uncertainty issues in the estimated velocities as well as in porosity models derived via petrophysical translations. In addition to traveltime inversion, the analysis of direct-wave amplitudes and reflection events provided further valuable information regarding subsurface properties and architecture. The used VRP amplitude preprocessing and inversion procedures were adapted from ray-based crosshole ground-penetrating radar (GPR) attenuation tomography and resulted in an attenuation model, which can be used to estimate variations in electrical resistivity. Our VRP reflection imaging approach relied on corridor stacking, which is a well-established processing sequence in vertical seismic profiling. The resulting reflection image outlines bounding layers and can be directly compared to surface-based GPR reflection profiling. Our results of the combined analysis of VRP, traveltimes, amplitudes, and reflections were consistent with independent core and borehole logs as well as GPR reflection profiles, which enabled us to derive a detailed hydrostratigraphic model as needed, for example, to understand and model groundwater flow and transport.

Geophysics ◽  
1992 ◽  
Vol 57 (1) ◽  
pp. 9-14 ◽  
Author(s):  
Gérard C. Herman

A nonlinear inversion method is presented, especially suited for the determination of global velocity models. In a certain sense, it can be considered as a generalization of methods based on traveltimes of reflections, with the requirement of accurately having to determine traveltimes replaced by the (less stringent and less subjective) requirement of having to define time windows around main reflections (or composite reflections) of interest. It is based on an error norm, related to the phase of the wavefield, which is directly computed from wavefield measurements. Therefore, the cumbersome step of interpreting arrivals and measuring arrival times is avoided. The method is applied to the reconstruction of a depth‐dependent global velocity model from a set of plane‐wave responses and is compared to other methods. Despite the fact that the new error norm only makes use of data having a temporal bandwidth of a few Hz, its behavior is very similar to the behavior of the error norm used in traveltime inversion.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. H61-H69
Author(s):  
Niklas Allroggen ◽  
Stéphane Garambois ◽  
Guy Sénéchal ◽  
Dominique Rousset ◽  
Jens Tronicke

Crosshole ground-penetrating radar (GPR) is applied in areas that require a very detailed subsurface characterization. Analysis of such data typically relies on tomographic inversion approaches providing an image of subsurface parameters. We have developed an approach for processing the reflected energy in crosshole GPR data and applied it on GPR data acquired in different sedimentary settings. Our approach includes muting of the first arrivals, separating the up- and the downgoing wavefield components, and backpropagating the reflected energy by a generalized Kirchhoff migration scheme. We obtain a reflection image that contains information on the location of electromagnetic property contrasts, thus outlining subsurface architecture in the interborehole plane. In combination with velocity models derived from different tomographic approaches, these images allow for a more detailed interpretation of subsurface structures without the need to acquire additional field data. In particular, a combined interpretation of the reflection image and the tomographic velocity model improves the ability to locate layer boundaries and to distinguish different subsurface units. To support our interpretations of our field data examples, we compare our crosshole reflection results with independent information, including borehole logs and surface GPR data.


Geophysics ◽  
2021 ◽  
pp. 1-145
Author(s):  
Zhiming Ren ◽  
Qianzong Bao ◽  
Bingluo Gu

Full waveform inversion (FWI) suffers from the local minima problem and requires a sufficiently accurate starting model to converge to the correct solution. Wave-equation traveltime inversion (WETI) is an effective tool to retrieve the long-wavelength components of the velocity model. We develop a joint diving/direct and reflected wave WETI (JDRWETI) method to build the P- and S-wave velocity macromodels. We estimate the traveltime shifts of seismic events (diving/direct waves, PP and PS reflections) through the dynamic warping scheme and construct a misfit function using both the time shifts of diving/direct and reflected waves. We derive the adjoint wave equations and the gradients with respect to the background models based on the joint misfit function. We apply the kernel decomposition scheme to extract the kernel of the diving/direct wave and the tomography kernels of PP and PS reflections. For an explosive source, the kernels of diving/direct wave and PP reflections and the kernel of PS reflections are used to compute the P- and S-wave gradients of the background models, respectively. We implement JDRWETI by a two-stage inversion workflow: first invert the P- and S-wave velocity models using the P-wave gradients and then improve the S-wave velocity model using the S-wave gradients. Numerical tests on synthetic and field datasets reveal that the JDRWETI method successfully recovers the long-wavelength components of P- and S-wave velocity models, which can be used for an initial model for the subsequent elastic FWI. Moreover, the proposed JDRWETI method prevails over the existing reflection WETI method and the cascaded diving/direct and reflected wave WETI method, especially when large velocity errors are present in the shallow part of the starting models. The JDRWETI method with the two-stage inversion workflow can give rise to reasonable inversion results even for the model with different P- and S-wave velocity structures.


2021 ◽  
Author(s):  
Akinniyi Akinsunmade

AbstractSpatial and depth variability of soil characteristics greatly influence its optimum utilization and management. Concealing nature of soil subsurface horizons has made the traditional soil investigations which rely on point information less reliable. In this study, an alternative use of ground penetrating radar (GPR)—a near-surface geophysical survey method—was tested to address the shortcomings. The focus of the study was on assessment of characteristics variability of soil layers at a test site and evaluation of effects of compaction caused by machinery traffics on soil. GPR methods utilize electromagnetic energy in the frequency range of 10 MHz and 3.0 GHz. Fourteen profiles GPR data were acquired at the test site-a farmland in Krakow, Poland. Compaction on parts of the soil was induced using tractor movements (simulating traffic effects) at different passes. Data were processed using basic filtering algorithms and attributes computations executed in Reflexw software. Attempt made in the study was on use of GPR geophysical technique for soil assessment. The method allows delineation of the soil horizons which depicts characteristic depth changes and spatial variability within the horizons. Moreover, traffic effects that caused compaction on parts of the soil horizons were discernable from the GPR profile sections. Thus, similar densification like hardpan that may develop in natural setting can be investigated using the method. The results have shown the suitability of the method for quick, noninvasive and continuous soil investigation that may also allow assessment of temporal soil changes via repeated measurement.


2021 ◽  
Vol 225 (2) ◽  
pp. 1020-1031
Author(s):  
Huachen Yang ◽  
Jianzhong Zhang ◽  
Kai Ren ◽  
Changbo Wang

SUMMARY A non-iterative first-arrival traveltime inversion method (NFTI) is proposed for building smooth velocity models using seismic diving waves observed on irregular surface. The new ray and traveltime equations of diving waves propagating in smooth media with undulant observation surface are deduced. According to the proposed ray and traveltime equations, an analytical formula for determining the location of the diving-wave turning points is then derived. Taking the influence of rough topography on first-arrival traveltimes into account, the new equations for calculating the velocities at turning points are established. Based on these equations, a method is proposed to construct subsurface velocity models from the observation surface downward to the bottom using the first-arrival traveltimes in common offset gathers. Tests on smooth velocity models with rugged topography verify the validity of the established equations, and the superiority of the proposed NFTI. The limitation of the proposed method is shown by an abruptly-varying velocity model example. Finally, the NFTI is applied to solve the static correction problem of the field seismic data acquired in a mountain area in the western China. The results confirm the effectivity of the proposed NFTI.


2004 ◽  
Vol 36 (3) ◽  
pp. 1396 ◽  
Author(s):  
O. C. Galanis ◽  
C. B. Papazachos ◽  
P. M. Hatzidimitriou ◽  
E. M. Scordilis

In the past years there has been a growing demand for precise earthquake locations for seismotectonic and seismic hazard studies. Recently this has become possible because of the development of sophisticated location algorithms, as well as hardware resources. This is expected to lead to a better insight of seismicity in the near future. A well-known technique, which has been recently used for relocating earthquake data sets is the double difference algorithm. In its original implementation it makes use of a one-dimensional ray tracing routine to calculate seismic wave travel times. We have modified the implementation of the algorithm by incorporating a three-dimensional velocity model and ray tracing in order to relocate a set of earthquakes in the area of the Mygdonia Basin (Northern Greece). This area is covered by a permanent regional network and occasionally by temporary local networks. The velocity structure is very well known, as the Mygdonia Basin has been used as an international test site for seismological studies since 1993, which makes it an appropriate location for evaluating earthquake location algorithms, with the quality of the results depending only on the quality of the data and the algorithm itself. The new earthquake locations reveal details of the area's seismotectonic structure, which are blurred, if not misleading, when resolved by standard (routine) location algorithms.


2018 ◽  
Vol 23 (4) ◽  
pp. 489-496
Author(s):  
J. David Redman ◽  
A. Peter Annan ◽  
Nectaria Diamanti

Bulk electrical properties of media are important inherently for ground penetrating radar (GPR) applications and for providing a means to determine indirectly other physical properties such as moisture content. We have developed a reflector whose reflectivity can be controlled electronically. This variable reflector controlled by a GPR provides an effective method to measure bulk electrical properties of media. For sample measurements, the GPR is placed on one side of a sample and the variable reflector on the opposite side. GPR trace data are then acquired with the reflector in an on-state and in the off-state. By differencing these measurements, we improve the ability to detect the specific reflection event from the variable reflector. This process removes both the direct wave and clutter from the trace data, improving the quality of the refection event and our ability to accurately pick its arrival time and amplitude. We describe the variable reflector, a prototype instrument based on the reflector and numerical modeling performed to understand its response. We also show the results of testing applications to the measurement of wood chip moisture content and monitoring of the electrical properties of concrete during the curing process.


2015 ◽  
Vol 19 (3) ◽  
pp. 1125-1139 ◽  
Author(s):  
P. Klenk ◽  
S. Jaumann ◽  
K. Roth

Abstract. High-resolution time-lapse ground-penetrating radar (GPR) observations of advancing and retreating water tables can yield a wealth of information about near-surface water content dynamics. In this study, we present and analyze a series of imbibition, drainage and infiltration experiments that have been carried out at our artificial ASSESS test site and observed with surface-based GPR. The test site features a complicated but known subsurface architecture constructed with three different kinds of sand. It allows the study of soil water dynamics with GPR under a wide range of different conditions. Here, we assess in particular (i) the feasibility of monitoring the dynamic shape of the capillary fringe reflection and (ii) the relative precision of monitoring soil water dynamics averaged over the whole vertical extent by evaluating the bottom reflection. The phenomenology of the GPR response of a dynamically changing capillary fringe is developed from a soil physical point of view. We then explain experimentally observed phenomena based on numerical simulations of both the water content dynamics and the expected GPR response.


Sign in / Sign up

Export Citation Format

Share Document