The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves

Geophysics ◽  
2014 ◽  
Vol 79 (6) ◽  
pp. R247-R261 ◽  
Author(s):  
Lisa Groos ◽  
Martin Schäfer ◽  
Thomas Forbriger ◽  
Thomas Bohlen

Full-waveform inversion (FWI) of Rayleigh waves is attractive for shallow geotechnical investigations due to the high sensitivity of Rayleigh waves to the S-wave velocity structure of the subsurface. In shallow-seismic field data, the effects of anelastic damping are significant. Dissipation results in a low-pass effect as well as frequency-dependent decay with offset. We found this by comparing recorded waveforms with elastic and viscoelastic wave simulation. The effects of anelastic damping must be considered in FWI of shallow-seismic Rayleigh waves. FWI using elastic simulation of wave propagation failed in synthetic inversion tests in which we tried to reconstruct the S-wave velocity in a viscoelastic model. To overcome this, [Formula: see text]-values can be estimated from the recordings to quantify viscoelasticity. Waveform simulation in the FWI then uses these a priori values when inferring seismic velocities and density. A source-wavelet correction, which is inevitable in FWI of field data, can compensate a significant fraction of the residuals between elastically and viscoelastically simulated data by narrowing the signals’ bandwidth. This way, elastic simulation becomes applicable in FWI of data from anelastic media. This approach, however, was not able to produce a frequency-dependent amplitude decay with offset. Reconstruction, therefore, was more accurate when using appropriate viscoelastic modeling in FWI of shallow-seismic Rayleigh waves. We found this by synthetic inversion tests using elastic forward simulation as well as viscoelastic simulation with different a priori values for [Formula: see text].

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. R109-R117 ◽  
Author(s):  
Lisa Groos ◽  
Martin Schäfer ◽  
Thomas Forbriger ◽  
Thomas Bohlen

The S-wave velocity of the shallow subsurface can be inferred from shallow-seismic Rayleigh waves. Traditionally, the dispersion curves of the Rayleigh waves are inverted to obtain the (local) S-wave velocity as a function of depth. Two-dimensional elastic full-waveform inversion (FWI) has the potential to also infer lateral variations. We have developed a novel workflow for the application of 2D elastic FWI to recorded surface waves. During the preprocessing, we apply a line-source simulation (spreading correction) and perform an a priori estimation of the attenuation of waves. The iterative multiscale 2D elastic FWI workflow consists of the preconditioning of the gradients in the vicinity of the sources and a source-wavelet correction. The misfit is defined by the least-squares norm of normalized wavefields. We apply our workflow to a field data set that has been acquired on a predominantly depth-dependent velocity structure, and we compare the reconstructed S-wave velocity model with the result obtained by a 1D inversion based on wavefield spectra (Fourier-Bessel expansion coefficients). The 2D S-wave velocity model obtained by FWI shows an overall depth dependency that agrees well with the 1D inversion result. Both models can explain the main characteristics of the recorded seismograms. The small lateral variations in S-wave velocity introduced by FWI additionally explain the lateral changes of the recorded Rayleigh waves. The comparison thus verifies the applicability of our 2D FWI workflow and confirms the potential of FWI to reconstruct shallow small-scale lateral changes of S-wave velocity.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. R105-R115 ◽  
Author(s):  
Edgar Manukyan ◽  
Hansruedi Maurer ◽  
André Nuber

Seismic full-waveform inversion (FWI) is potentially a powerful method for obtaining high-resolution subsurface images, but the results are often distorted by nonlinear effects and parameter trade-offs. Such distortions can be particularly severe in the case of multiparameter FWI, such as elastic FWI, in which inversion is performed simultaneously for P- and S-wave velocities and density. The problem can be alleviated by adding constraints in the form of plausible a priori information. A usually well-justified constraint includes the structural similarity of different model parameters; i.e., an anomalous body likely exhibits variations in all elastic properties, although their magnitudes may be different. To consider such types of a priori information, we have developed a structurally constrained elastic FWI, which is based on minimization of the cross products of gradients of different model parameters. Our synthetic 2D experiments show that structurally constrained FWI can significantly improve model reconstruction. It is also demonstrated that our approach still leads to improved results, even when the structural similarity between the individual parameter types is not exactly met. Inversions of field data show that in comparison to conventional FWI, structurally constrained FWI is able to match the field data equally well while requiring less structural complexity of the subsurface.


2021 ◽  
Author(s):  
Tan Qin ◽  
Thomas Bohlen ◽  
Yudi Pan

<p>Shallow-seismic surface wave and ground penetrating radar (GPR) are employed in a wide range of engineering and geosciences applications. Full-waveform inversion (FWI) of either seismic or multi-offset GPR data are able to provide high-resolution subsurface characterization and have received particular attention in the past decade. Those two geophysical methods are involved in the increasing requirements of comprehensive site and material investigations. However, it is still challenging to provide an effective integration between seismic data and electromagnetic data. In this paper, we investigated the joint petrophysical inversion (JPI) of shallow-seismic and multi-offset GPR data for more consistent imaging of near surface. As a bridge between the seismic parameters (P-wave velocity, S-wave velocity, and density) and GPR parameters (relative dielectric permittivity and electric conductivity), the petrophysical relationships with the parameters namely porosity and saturation are employed to link two data sets. We first did a sensitivity analysis of the petrophysical parameters to the seismic and GPR parameters and then determined an efficient integration of using shallow-seismic FWI to update porosity and GPR FWI to update saturation, respectively. A comparison of several parameterisation combinations shows that the seismic velocity parameterisation in shallow-seismic FWI and a modified logarithm parameterisation in GPR FWI works well in reconstructing reliable S-wave velocity and relative dielectric permittivity models, respectively. With the help from the petrophysical links, we realized JPI by transforming those well recovered parameters to the petrophysical parameters and then to other seismic and GPR parameters. A synthetic test indicates that, compared with the individual petrophysical inversion and individual FWI, JPI outperforms in simultaneously reconstructing all seismic, GPR, and petrophysical parameters with higher resolution and improved details. It is proved that JPI would be a potential data integration approach for the shallow subsurface investigation.</p>


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B95-B105 ◽  
Author(s):  
Yao Wang ◽  
Richard D. Miller ◽  
Shelby L. Peterie ◽  
Steven D. Sloan ◽  
Mark L. Moran ◽  
...  

We have applied time domain 2D full-waveform inversion (FWI) to detect a known 10 m deep wood-framed tunnel at Yuma Proving Ground, Arizona. The acquired seismic data consist of a series of 2D survey lines that are perpendicular to the long axis of the tunnel. With the use of an initial model estimated from surface wave methods, a void-detection-oriented FWI workflow was applied. A straightforward [Formula: see text] quotient masking method was used to reduce the inversion artifacts and improve confidence in identifying anomalies that possess a high [Formula: see text] ratio. Using near-surface FWI, [Formula: see text] and [Formula: see text] velocity profiles were obtained with void anomalies that are easily interpreted. The inverted velocity profiles depict the tunnel as a low-velocity anomaly at the correct location and depth. A comparison of the observed and simulated waveforms demonstrates the reliability of inverted models. Because the known tunnel has a uniform shape and for our purposes an infinite length, we apply 1D interpolation to the inverted [Formula: see text] profiles to generate a pseudo 3D (2.5D) volume. Based on this research, we conclude the following: (1) FWI is effective in near-surface tunnel detection when high resolution is necessary. (2) Surface-wave methods can provide accurate initial S-wave velocity [Formula: see text] models for near-surface 2D FWI.


Sign in / Sign up

Export Citation Format

Share Document