Robust full waveform inversion of surface waves

Masoni, I.; Brossier, Romain; Virieux, Jean; Boëlle, Jean-Luc

doi:10.1190/segam2014-1077.1

Cited by 12 publications

(8 citation statements)

References 7 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Reducing the nonlinearity of waveform inversion problems involves a judicious choice of the objective function (Luo and Schuster, 1991;Gee and Jordan, 1992;Dahlen and Baig, 2002;Fichtner et al, 2008;Bozdağ et al, 2011;Rickers et al, 2012;Masoni et al, 2014). Bozdağ et al (2011), in particular, discuss several misfit functions in seismic tomography.…”

Section: Surface Waves: a New Approachmentioning

confidence: 99%

Multiscale adjoint waveform tomography for surface and body waves

2015

View full text Add to dashboard Cite

International audienceWe have developed a wavelet-multiscale adjoint scheme for the elastic full-waveform inversion of seismic data, including body waves (BWs) and surface waves (SWs). We start the inversion on the SW portion of the seismograms. To avoid cycle skipping and reduce the dependence on the initial model of these dispersive waves, we commence by minimizing an envelope-based misfit function. Subsequently, we proceed to the minimization of a waveform-difference (WD) metric applied to the SWs only. After that, we fit BWs and SWs indiscriminately using a WD misfit metric. In each of these three steps, we guide the iterative inversion through a sequence of nested subspace projections in a wavelet basis. SW analysis preserves a wealth of near-surface features that would be lost in conventional BW tomography. We used a toy model to illustrate the dispersive and cycle-skipping behavior of the SWs, and to introduce the two ways by which we combat the nonlinearity of waveform inversions involving SWs. The first is the wavelet-based multiscale character of the method, and the second the envelope-based misfit function. Next, we used an industry synthetic model to perform realistic numerical experiments to further develop a strategy for SW and joint SW as well as BW tomography. The effect of incorrect density information on wave-speed inversions was also evaluated. We ultimately formalize a flexible scheme for full-waveform inversion based on adjoint methods that includes BWs and SWs, and also considers P- and S-wave speeds, as well as density. Our method is applicable to waveform inversion in exploration geophysics, geotechnical engineering, regional, and global seismology.Read More: http://library.seg.org/doi/10.1190/geo2014-0461.

show abstract

Section: Surface Waves: a New Approachmentioning

confidence: 99%

Multiscale adjoint waveform tomography for surface and body waves

2015

View full text Add to dashboard Cite

show abstract

“…(), Masoni et al . (), Masoni et al . () and Pérez Solano, Donno and Chauris () increased the robustness of 2D surface wave inversion by defining new cost functions in different data domains.…”

Section: Introductionmentioning

confidence: 94%

“…However, the extension of this concept to 2D surface wave inversion using local optimisation methods is a challenging task, due to potential cycle-skipping problems in case of the strongly dispersive surface wavefield. Masoni et al (2013), Masoni et al (2014a), Masoni et al (2014b) and Pérez Solano, Donno and Chauris (2014) increased the robustness of 2D surface wave inversion by defining new cost functions in different data domains. Apart from methodical developments and synthetic inversion tests, the number of 2D surface wave field data applications are currently very limited.…”

mentioning

confidence: 99%

Full waveform inversion of SH‐ and Love‐wave data in near‐surface prospecting

Dokter

Köhn

Wilken

et al. 2017

Geophysical Prospecting

View full text Add to dashboard Cite

We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village Čachtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls.

show abstract

“…Masoni et al . (, a,b) and Pérez Solano, Donno, and Chauris () increased the robustness of 2D surface‐wave inversion by defining new cost functions in different data domains. A 2D time‐domain Gauss–Newton Rayleigh‐wave FWI was developed by Tran and McVay () for geotechnical site characterization, such as the detection of sinkholes (Tran et al .…”

Section: Introductionmentioning

confidence: 99%

“…The extension to the 2D surface-wave inversion using local optimisation methods is a challenging task, due to potential cycle-skipping problems in case of the strongly dispersive surface wave field in the time domain. Masoni et al (2013Masoni et al ( , 2014a and Pérez Solano, Donno, and Chauris (2014) increased the robustness of 2D surface-wave inversion by defining new cost functions in different data domains. A 2D time-domain Gauss-Newton Rayleigh-wave FWI was developed by Tran and McVay (2012) for geotechnical site charac-the frequency content of the recorded data, are taken into account.…”

mentioning

confidence: 99%

Application of 2D elastic Rayleigh waveform inversion to ultrasonic laboratory and field data

Köhn

Meier

Fehr

et al. 2016

Near Surface Geophysics

View full text Add to dashboard Cite

In addition to geophysical applications from the near surface to a global scale, seismic full‐waveform inversion can be applied to ultrasonic data on the centimetre and decimetre scales for nondestructive testing of pavements, facades, plaster, sculptures, and load‐bearing structures such as pillars or core samples from boreholes, which can consist of geo‐materials and non‐geomaterials. Classical non‐destructive testing approaches are based on the inversion of body‐wave travel times to deduce P‐wave velocity models. In contrast, surface waves (especially Rayleigh waves) are well suited to quantify surficial alterations of material properties, e.g., due to weathering. Furthermore, ultrasonic measurements of test samples with known material parameters close the gap between synthetic tests or analytical solutions and field data applications to estimate the accuracy of seismic modelling and inversion codes. Due to the scale invariance of the problem, the full‐waveform‐inversion approaches developed on the ultrasonic scale are also applicable to larger scale geophysical problems. In this paper, we demonstrate the potential of two‐dimensional Rayleigh waveform inversion on the ultrasonic scale using two data examples acquired with a single‐fold, low‐coverage acquisition geometry. For a simple homogeneous Plexiglas block in a controlled laboratory environment, we discuss the accuracy of visco‐elastic Rayleigh wave modelling, as well as the sensitivity of two‐dimensional elastic full‐waveform inversion with respect to small data residuals. While the inversion is restricted to the recovery of the S‐wave velocity model, a passive visco‐elastic modelling approach with a homogeneous average Qs‐model is required in order to describe amplitude loss and dispersion of the Rayleigh wave. The applicability of this approach under field conditions is illustrated for an ultrasonic data example from the weathered sandstone facade of the Porta Nigra in Trier (Germany). In addition to a random medium resolution analysis of the full‐waveform‐inversion result, we particularly emphasise the importance of lateral model smoothing to mitigate small‐scale inversion artefacts and to avoid erroneous interpretations. The estimated two‐dimensional S‐wave velocity anomalies correlate well with prominent surficial weathering effects in the upper about 3 cm.

show abstract

Robust full waveform inversion of surface waves

Cited by 12 publications

References 7 publications

Multiscale adjoint waveform tomography for surface and body waves

Multiscale adjoint waveform tomography for surface and body waves

Full waveform inversion of SH‐ and Love‐wave data in near‐surface prospecting

Application of 2D elastic Rayleigh waveform inversion to ultrasonic laboratory and field data

Contact Info

Product

Resources

About