scholarly journals 3D wave-equation dispersion inversion of surface waves recorded on irregular topography

Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. R147-R161 ◽  
Author(s):  
Zhaolun Liu ◽  
Jing Li ◽  
Sherif M. Hanafy ◽  
Kai Lu ◽  
Gerard Schuster
Keyword(s):  
Wave Equation ◽  
Surface Wave ◽  
Surface Waves ◽  
Wave Velocity ◽  
Field Data ◽  
P Wave ◽  
Topographic Effects ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Irregular Topography

Irregular topography can cause strong scattering and defocusing of propagating surface waves, so it is important to account for such effects when inverting surface waves for shallow S-wave velocity structures. We have developed a 3D surface-wave dispersion inversion method that takes into account the topographic effects modeled by a 3D spectral element solver. The objective function is the frequency summation of the squared wavenumber differences [Formula: see text] along each azimuthal angle of the fundamental mode or higher-order modes of Rayleigh waves in each shot gather. The wavenumbers [Formula: see text] associated with the dispersion curves are calculated using the data recorded along the irregular free surface. Numerical tests on synthetic and field data demonstrate that 3D topographic wave equation dispersion inversion (TWD) can accurately invert for the S-wave velocity model from surface-wave data recorded on irregular topography. Field data tests for data recorded across an Arizona fault demonstrate that, for this example, the 2D TWD model can be as accurate as the 3D tomographic model. This suggests that in some cases, the 2D TWD inversion is preferred over 3D TWD because of its significant reduction in computational costs. Compared to the 3D P-wave velocity tomogram, the 3D S-wave tomogram agrees much more closely with the geologic model taken from the trench log. The agreement with the trench log is even better when the [Formula: see text] tomogram is computed, which reveals a sharp change in velocity across the fault. The localized velocity anomaly in the [Formula: see text] tomogram is in very good agreement with the well log. Our results suggest that integrating the [Formula: see text] and [Formula: see text] tomograms can sometimes give the most accurate estimates of the subsurface geology across normal faults.

Crustal structure of the Iceland region from dispersion of surface waves

1962 ◽  
Vol 52 (2) ◽  
pp. 359-388
Author(s):  
Eysteinn Tryggvason
Keyword(s):  
Surface Layer ◽  
Surface Wave ◽  
Surface Waves ◽  
Wave Velocity ◽  
Crustal Structure ◽  
Wave Dispersion ◽  
P Wave ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Mid Atlantic Ridge

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

Full-waveform matching of VP and VS models from surface waves

2019 ◽  
Vol 218 (3) ◽  
pp. 1873-1891 ◽  
Author(s):  
Farbod Khosro Anjom ◽  
Daniela Teodor ◽  
Cesare Comina ◽  
Romain Brossier ◽  
Jean Virieux ◽  
...  
Keyword(s):  
Surface Waves ◽  
Wave Velocity ◽  
Real Data ◽  
Test Site ◽  
P Wave ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Data Set ◽  
Near Surface ◽  
Velocity Models

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.

A comprehensive comparison between the refraction microtremor and seismic interferometry methods for phase-velocity estimation

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. EN99-EN108 ◽  
Author(s):  
Zongbo Xu ◽  
T. Dylan Mikesell ◽  
Jianghai Xia ◽  
Feng Cheng
Keyword(s):  
Surface Wave ◽  
Surface Waves ◽  
Wave Velocity ◽  
Velocity Estimation ◽  
Seismic Interferometry ◽  
S Wave ◽  
Noise Sources ◽  
Surface Wave Dispersion ◽  
Near Surface ◽  
Refraction Microtremor

Passive-source seismic-noise-based surface-wave methods are now routinely used to investigate the near-surface geology in urban environments. These methods estimate the S-wave velocity of the near surface, and two methods that use linear recording arrays are seismic interferometry (SI) and refraction microtremor (ReMi). These two methods process noise data differently and thus can yield different estimates of the surface-wave dispersion, the data used to estimate the S-wave velocity. We have systematically compared these two methods using synthetic data with different noise source distributions. We arrange sensors in a linear survey grid, which is conveniently used in urban investigations (e.g., along roads). We find that both methods fail to correctly determine the low-frequency dispersion characteristics when outline noise sources become stronger than inline noise sources. We also identify an artifact in the ReMi method and theoretically explain the origin of this artifact. We determine that SI combined with array-based analysis of surface waves is the more accurate method to estimate surface-wave phase velocities because SI separates surface waves propagating in different directions. Finally, we find a solution to eliminate the ReMi artifact that involves the combination of SI and the [Formula: see text]-[Formula: see text] transform, the array processing method that underlies the ReMi method.

scholarly journals Two-Dimensional Full-Waveform Joint Inversion of Surface Waves Using Phases and Z/H Ratios

2021 ◽  
Vol 11 (15) ◽  
pp. 6712
Author(s):  
Chao Zhang ◽  
Ting Lei ◽  
Yi Wang
Keyword(s):  
Surface Wave ◽  
Surface Waves ◽  
Wave Velocity ◽  
Joint Inversion ◽  
Waveform Inversion ◽  
Small Scale ◽  
Imaging Method ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Full Waveform

Surface-wave dispersion and the Z/H ratio are important parameters used to resolve the Earth’s structure, especially for S-wave velocity. Several previous studies have explored using joint inversion of these two datasets. However, all of these studies used a 1-D depth-sensitivity kernel, which lacks precision when the structure is laterally heterogeneous. Adjoint tomography (i.e., full-waveform inversion) is a state-of-the-art imaging method with a high resolution. It can obtain better-resolved lithospheric structures beyond the resolving ability of traditional ray-based travel-time tomography. In this study, we present a systematic investigation of the 2D sensitivities of the surface wave phase and Z/H ratio using the adjoint-state method. The forward-modeling experiments indicated that the 2D phase and Z/H ratio had different sensitivities to the S-wave velocity. Thus, a full-waveform joint-inversion scheme of surface waves with phases and a Z/H ratio was proposed to take advantage of their complementary sensitivities to the Earth’s structure. Both applications to synthetic data sets in large- and small-scale inversions demonstrated the advantage of the joint inversion over the individual inversions, allowing for the creation of a more unified S-wave velocity model. The proposed joint-inversion scheme offers a computationally efficient and inexpensive alternative to imaging fine-scale shallow structures beneath a 2D seismic array.

Shear Wave Velocity Analysis of 2-D Multichannel Analysis of Surface Wave (MASW) to investigate subsurface Fault of Alternative Bridge Construction in Kelok Sago Jambi

2020 ◽  
Vol 25 (1) ◽  
pp. 18-20
Author(s):  
Muhammad Farhan ◽  
◽  
Gunawan Handayani ◽  
Keyword(s):  
Surface Wave ◽  
Surface Waves ◽  
Wave Velocity ◽  
Urban Areas ◽  
Signal To Noise Ratio ◽  
Body Wave ◽  
Geophysical Methods ◽  
P Wave ◽  
S Wave ◽  
Multichannel Analysis

Every geotechnical measurement requires geophysical methods to classify soil types under the ground. S-wave velocity (Vs), P-wave velocity (Vp), and density (ρ), are the most important parameters in the classification of soils. There are various methods to determine Vs, one of them is P-S logging method. However, this method is less suitable to be applied in urban areas due to the difficulties of data acquisition and high expense in operational costs. In 1999, a seismic method uses surface waves to de-termine Vs profile with a higher signal to noise ratio which was known by the name of Multichannel Analysis of Surface Waves (MASW). A surface wave, especially Rayleigh wave, creeps slowly on the surface with a larger amplitude than a body wave. The wavelengths of the surface wave will disperse in the layers system i.e. the phase velocity of the surface waves is now func-tion of frequency. MASW 2-D method is used in this paper to determine subsoil properties and to identify the fault under the bridge abutments plan (abutment 1 and abutment 2) in Kelok Sago Jambi.

scholarly journals Joint Body‐ and Surface‐Wave Tomography of Yucca Flat, Nevada, Using a Novel Seismic Source

2019 ◽  
Vol 109 (5) ◽  
pp. 1922-1934 ◽  
Author(s):  
Liam D. Toney ◽  
Robert E. Abbott ◽  
Leiph A. Preston ◽  
David G. Tang ◽  
Tori Finlay ◽  
...  
Keyword(s):  
Surface Wave ◽  
Phase Velocity ◽  
Wave Velocity ◽  
Velocity Structure ◽  
Joint Inversion ◽  
P Wave ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Velocity Models ◽  
Phase Velocity Dispersion

Abstract In preparation for the next phase of the Source Physics Experiments, we acquired an active‐source seismic dataset along two transects totaling more than 30 km in length at Yucca Flat, Nevada, on the Nevada National Security Site. Yucca Flat is a sedimentary basin which has hosted more than 650 underground nuclear tests (UGTs). The survey source was a novel 13,000 kg modified industrial pile driver. This weight drop source proved to be broadband and repeatable, richer in low frequencies (1–3 Hz) than traditional vibrator sources and capable of producing peak particle velocities similar to those produced by a 50 kg explosive charge. In this study, we performed a joint inversion of P‐wave refraction travel times and Rayleigh‐wave phase‐velocity dispersion curves for the P‐ and S‐wave velocity structure of Yucca Flat. Phase‐velocity surface‐wave dispersion measurements were obtained via the refraction microtremor method on 1 km arrays, with 80% overlap. Our P‐wave velocity models verify and expand the current understanding of Yucca Flat’s subsurface geometry and bulk properties such as depth to Paleozoic basement and shallow alluvium velocity. Areas of disagreement between this study and the current geologic model of Yucca Flat (derived from borehole studies) generally correlate with areas of widely spaced borehole control points. This provides an opportunity to update the existing model, which is used for modeling groundwater flow and radionuclide transport. Scattering caused by UGT‐related high‐contrast velocity anomalies substantially reduced the number and frequency bandwidth of usable dispersion picks. The S‐wave velocity models presented in this study agree with existing basin‐wide studies of Yucca Flat, but are compromised by diminished surface‐wave coherence as a product of this scattering. As nuclear nonproliferation monitoring moves from teleseismic to regional or even local distances, such high‐frequency (>5  Hz) scattering could prove challenging when attempting to discriminate events in areas of previous testing.

Nonlinear two‐dimensional elastic inversion of multioffset seismic data

Geophysics ◽  
1987 ◽  
Vol 52 (9) ◽  
pp. 1211-1228 ◽  
Author(s):  
Peter Mora
Keyword(s):  
Wave Equation ◽  
Wave Velocity ◽  
A Priori ◽  
Gaussian Model ◽  
Wave Impedance ◽  
P Wave ◽  
S Wave ◽  
Elastic Wave Equation ◽  
Gradient Direction ◽  
P Wave Velocity

The treatment of multioffset seismic data as an acoustic wave field is becoming increasingly disturbing to many geophysicists who see a multitude of wave phenomena, such as amplitude‐offset variations and shearwave events, which can only be explained by using the more correct elastic wave equation. Not only are such phenomena ignored by acoustic theory, but they are also treated as undesirable noise when they should be used to provide extra information, such as S‐wave velocity, about the subsurface. The problems of using the conventional acoustic wave equation approach can be eliminated via an elastic approach. In this paper, equations have been derived to perform an inversion for P‐wave velocity, S‐wave velocity, and density as well as the P‐wave impedance, S‐wave impedance, and density. These are better resolved than the Lamé parameters. The inversion is based on nonlinear least squares and proceeds by iteratively updating the earth parameters until a good fit is achieved between the observed data and the modeled data corresponding to these earth parameters. The iterations are based on the preconditioned conjugate gradient algorithm. The fundamental requirement of such a least‐squares algorithm is the gradient direction which tells how to update the model parameters. The gradient direction can be derived directly from the wave equation and it may be computed by several wave propagations. Although in principle any scheme could be chosen to perform the wave propagations, the elastic finite‐ difference method is used because it directly simulates the elastic wave equation and can handle complex, and thus realistic, distributions of elastic parameters. This method of inversion is costly since it is similar to an iterative prestack shot‐profile migration. However, it has greater power than any migration since it solves for the P‐wave velocity, S‐wave velocity, and density and can handle very general situations including transmission problems. Three main weaknesses of this technique are that it requires fairly accurate a priori knowledge of the low‐ wavenumber velocity model, it assumes Gaussian model statistics, and it is very computer‐intensive. All these problems seem surmountable. The low‐wavenumber information can be obtained either by a prior tomographic step, by the conventional normal‐moveout method, by a priori knowledge and empirical relationships, or by adding an additional inversion step for low wavenumbers to each iteration. The Gaussian statistics can be altered by preconditioning the gradient direction, perhaps to make the solution blocky in appearance like well logs, or by using large model variances in the inversion to reduce the effect of the Gaussian model constraints. Moreover, with some improvements to the algorithm and more parallel computers, it is hoped the technique will soon become routinely feasible.

P- and S-wave velocity models from surface wave dispersion curves data transform

2017 ◽  
Author(s):  
Valentina Socco ◽  
Farbod Khosro Anjom ◽  
Cesare Comina ◽  
Daniela Teodor
Keyword(s):  
Surface Wave ◽  
Wave Velocity ◽  
Wave Dispersion ◽  
Dispersion Curves ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Velocity Models

Multi-Window Weighted Stacking of Surface-Wave Dispersion

Geophysics ◽  
2020 ◽  
pp. 1-53
Author(s):  
Sylvain Pasquet ◽  
Wei Wang ◽  
Po Chen ◽  
Brady A. Flinchum
Keyword(s):  
Surface Wave ◽  
Surface Waves ◽  
Real Data ◽  
Wave Dispersion ◽  
Lateral Resolution ◽  
Data Sets ◽  
S Wave ◽  
Surface Wave Dispersion ◽  
Shallow Subsurface ◽  
Depth Of Investigation

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.

Sign in / Sign up

Export Citation Format

Share Document