Multi-Window Weighted Stacking of Surface-Wave Dispersion

Geophysics ◽  
2020 ◽  
pp. 1-53
Author(s):  
Sylvain Pasquet ◽  
Wei Wang ◽  
Po Chen ◽  
Brady A. Flinchum

Surface wave (SW) methods are classically used to characterize shear (S-) wave velocities ( VS) of the shallow subsurface through the inversion of dispersion curves. When targeting 2D shallow structures with sharp lateral heterogeneity, windowing and stacking techniques can be implemented to provide a better description of VS lateral variations. These techniques, however, suffer from the trade-off between lateral resolution and depth of investigation, well-known when using multichannel analysis of surface waves (MASW). We propose a novel methodology aimed at enhancing both lateral resolution and depth of investigation of MASW results through the use of multi-window weighted stacking of surface waves (MW-WSSW). MW-WSSW consists in stacking dispersion images obtained from data segments of different sizes, with a wavelength-based weight that depends on the aperture of the data selection window. In that sense, MW-WSSW provides additional weight to short wavelengths in smaller windows so as to better inform shallow parts of the subsurface, and vice versa for deeper velocities. Using multiple windows improves the depth of investigation, while applying wavelength-based weights enhances shallow lateral resolution. MW-WSSW was implemented within the open-source package SWIP, and applied to the processing of synthetic and real data sets. In both cases we compared it with standard windowing and stacking procedures that are already implemented in SWIP. MW-WSSW provided convincing results with optimized lateral extent, improved shallow resolution, and increased depth of investigation.

1962 ◽  
Vol 52 (2) ◽  
pp. 359-388
Author(s):  
Eysteinn Tryggvason

ABSTRACT A number of Icelandic records of earthquakes originating in the Mid-Atlantic Seismic Belt between 52° and 70° N. lat. have been investigated. The surface waves on these records are chiefly in the period interval 3–10 sec, and are first mode Love-waves and Rayleigh-waves. The surface wave dispersion can be explained by a three-layered crustal structure as follows. A surface layer of S-wave velocity about 2.7 km/sec covering the whole region studied, a second layer of S-wave velocity about 3.6 km/sec covering Iceland and extending several hundred kilometers off the coasts and a third layer of S-wave velocity about 4.3 km/sec and P-wave velocity about 7.4 km/sec underlying the whole region. The thickness of the surface layer appears to be about 4 km on the Mid-Atlantic Ridge south of Iceland and in western Iceland, 3 km in central Iceland and 7 km northwest of Iceland. The second layer is apparently of similar thickness than the surface layer, while the third layer is thick; and the surface wave dispersion does not indicate any layer of higher wave velocity. This 7.4-layer is supposed to belong to the mantle, although its wave velocity is significantly lower than usually found in the upper mantle


2020 ◽  
Author(s):  
Shufan Hu ◽  
Yonghui Zhao ◽  
Wenda Bi ◽  
Ruiqing Shen ◽  
Bo Li ◽  
...  

<p>Ground penetrating radar (GPR) and Seismic Surface Wave methods (SWMs) are two nondestructive testing (NDT) methods commonly used in near-surface site investigations. These two methods investigate the media properties of subsurface based on different physical phenomena. GPR has a good resolvability to characterize the layered structure since the propagation of electromagnetic wave is sensitive to the change of electrical properties, while, the geometric dispersion of surface waves can be used to retrieve the variation of S-wave velocity (<em>V</em>s) with depth. In most situations, these two data sets are processed separately, and the results are later used for comprehensive interpretation. Constrained inversion, as a way to implement data fusion, can alleviate the non-uniqueness of the solution and produce more consistent information for the comprehensive site and material investigations.</p><p>We present an algorithm for the inversion of surface-wave dispersion curves with GPR interface constraints in 2D media. The reflection interfaces interpreted from the GPR profile are integrated into a cell- and boundary-based <em>V</em>s model. This implementation allows both vertical and lateral changes within each region while also allows sharp changes across the boundaries. In addition, our algorithm simultaneously inverts several dispersion curves extracted along the survey line using multi-size spatial windows, which mitigates the adverse effects of 1D assumption in traditional surface-wave dispersion inversion and improves the matching of GPR and SWMs in lateral variations. We use synthetic and field data sets to test the effectivity of the proposed method. Both results show the improved resolution of the <em>V</em>s model retrieved by the constrained inversion compared to the standard inversion.</p>


2013 ◽  
Vol 353-356 ◽  
pp. 1196-1202 ◽  
Author(s):  
Jian Qi Lu ◽  
Shan You Li ◽  
Wei Li

Surface wave dispersion imaging approach is crucial for multi-channel analysis of surface wave (MASW). Because the resolution of inversed S-wave velocity and thickness of a layer are directly subjected to the resolution of imaged dispersion curve. The τ-p transform approach is an efficient and commonly used approach for Rayleigh wave dispersion curve imaging. However, the conventional τ-p transform approach was severely affected by waves amplitude. So, the energy peaks of f-v spectrum were mainly gathered in a narrow frequency range. In order to remedy this shortage, an improved τ-p transform approach was proposed by this paper. Comparison has been made between phase shift and improved τ-p transform approaches using both synthetic and in situ tested data. Result shows that the dispersion image transformed from proposed approach is superior to that either from conventionally τ-p transform or from phase shift approaches.


2017 ◽  
Author(s):  
Valentina Socco ◽  
Farbod Khosro Anjom ◽  
Cesare Comina ◽  
Daniela Teodor

1967 ◽  
Vol 57 (5) ◽  
pp. 959-981
Author(s):  
Victor Gregson

abstract Elastic waves produced by an impact were recorded at the surface of a solid 12.0 inch diameter steel sphere coated with a 0.3 inch copper layer. Conventional modeling techniques employing both compressional and shear piezoelectric transducers were used to record elastic waves for one millisecond at various points around the great circle of the sphere. Body, PL, and surface waves were observed. Density, layer thickness, compressional and shear-wave velocities were measured so that accurate surface-wave dispersion curves could be computed. Surface-wave dispersion was measured as well as computed. Measured PL mode dispersion compared favorably with theoretical computations. In addition, dispersion curves for Rayleigh, Stoneley, and Love modes were computed. Measured surface-wave dispersion showed Rayleigh and Love modes were observed but not Stoneley modes. Measured dispersion compared favorably with theoretical computations. The curvature correction applied to dispersion calculations in a flat space has been estimated to correct dispersion values at long-wave lengths to about one per cent of correct dispersion in a spherical model. Measured dispersion compared with such flat space dispersion corrected for curvature proved accurate within one per cent at long wave lengths. Two sets of surface waves were observed. One set was associated with body waves radiating outward from impact. The other set was associated with body waves reflecting at the pole opposite impact. For each set of surface waves, measured dispersion compared favorably with computed dispersion.


Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. EN1-EN15 ◽  
Author(s):  
Paolo Bergamo ◽  
Ben Dashwood ◽  
Sebastian Uhlemann ◽  
Russell Swift ◽  
Jonathan E. Chambers ◽  
...  

The UK’s transportation network is supported by critical geotechnical assets (cuttings/embankments/dams) that require sustainable, cost-effective management, while maintaining an appropriate service level to meet social, economic, and environmental needs. Recent effects of extreme weather on these geotechnical assets have highlighted their vulnerability to climate variations. We have assessed the potential of surface wave data to portray the climate-related variations in mechanical properties of a clay-filled railway embankment. Seismic data were acquired bimonthly from July 2013 to November 2014 along the crest of a heritage railway embankment in southwest England. For each acquisition, the collected data were first processed to obtain a set of Rayleigh-wave dispersion and attenuation curves, referenced to the same spatial locations. These data were then analyzed to identify a coherent trend in their spatial and temporal variability. The relevance of the observed temporal variations was also verified with respect to the experimental data uncertainties. Finally, the surface wave dispersion data sets were inverted to reconstruct a time-lapse model of S-wave velocity for the embankment structure, using a least-squares laterally constrained inversion scheme. A key point of the inversion process was constituted by the estimation of a suitable initial model and the selection of adequate levels of spatial regularization. The initial model and the strength of spatial smoothing were then kept constant throughout the processing of all available data sets to ensure homogeneity of the procedure and comparability among the obtained [Formula: see text] sections. A continuous and coherent temporal pattern of surface wave data, and consequently of the reconstructed [Formula: see text] models, was identified. This pattern is related to the seasonal distribution of precipitation and soil water content measured on site.


2020 ◽  
Vol 224 (3) ◽  
pp. 2077-2099
Author(s):  
J K Magali ◽  
T Bodin ◽  
N Hedjazian ◽  
H Samuel ◽  
S Atkins

SUMMARY In the Earth’s upper mantle, seismic anisotropy mainly originates from the crystallographic preferred orientation (CPO) of olivine due to mantle deformation. Large-scale observation of anisotropy in surface wave tomography models provides unique constraints on present-day mantle flow. However, surface waves are not sensitive to the 21 coefficients of the elastic tensor, and therefore the complete anisotropic tensor cannot be resolved independently at every location. This large number of parameters may be reduced by imposing spatial smoothness and symmetry constraints to the elastic tensor. In this work, we propose to regularize the tomographic problem by using constraints from geodynamic modelling to reduce the number of model parameters. Instead of inverting for seismic velocities, we parametrize our inverse problem directly in terms of physical quantities governing mantle flow: a temperature field, and a temperature-dependent viscosity. The forward problem consists of three steps: (1) calculation of mantle flow induced by thermal anomalies, (2) calculation of the induced CPO and elastic properties using a micromechanical model, and (3) computation of azimuthally varying surface wave dispersion curves. We demonstrate how a fully nonlinear Bayesian inversion of surface wave dispersion curves can retrieve the temperature and viscosity fields, without having to explicitly parametrize the elastic tensor. Here, we consider simple flow models generated by spherical temperature anomalies. The results show that incorporating geodynamic constraints in surface wave inversion help to retrieve patterns of mantle deformation. The solution to our inversion problem is an ensemble of models (i.e. thermal structures) representing a posterior probability, therefore providing uncertainties for each model parameter.


2019 ◽  
Vol 218 (3) ◽  
pp. 1873-1891 ◽  
Author(s):  
Farbod Khosro Anjom ◽  
Daniela Teodor ◽  
Cesare Comina ◽  
Romain Brossier ◽  
Jean Virieux ◽  
...  

SUMMARY The analysis of surface wave dispersion curves (DCs) is widely used for near-surface S-wave velocity (VS) reconstruction. However, a comprehensive characterization of the near-surface requires also the estimation of P-wave velocity (VP). We focus on the estimation of both VS and VP models from surface waves using a direct data transform approach. We estimate a relationship between the wavelength of the fundamental mode of surface waves and the investigation depth and we use it to directly transform the DCs into VS and VP models in laterally varying sites. We apply the workflow to a real data set acquired on a known test site. The accuracy of such reconstruction is validated by a waveform comparison between field data and synthetic data obtained by performing elastic numerical simulations on the estimated VP and VS models. The uncertainties on the estimated velocity models are also computed.


Sign in / Sign up

Export Citation Format

Share Document