Retrieving reflection arrivals from passive seismic data using Radon correlation 

Author(s):  
Diako Hariri Naghadeh ◽  
Christopher J Bean ◽  
Patrick Smith ◽  
Sergei Lebedev ◽  
Huda Mohamed

<p>Since explosive and impulsive seismic sources such as dynamite, air guns, gas guns, or even vibroseis can have a big impact on the environment, some companies have decided to record ambient seismic noise and use it to estimate the physical properties of the subsurface. Big challenges arise when the aim is extracting body-waves from recorded passive signals, especially in the presence of strong surface waves. In passive seismic signals, such body-waves are usually weak in comparison to surface waves which are much more prominent. To understand the characteristics of passive signals and the effect of natural source locations, three simple synthetic models were created. To extract body-waves from simulated passive signals we propose and test a Radon-correlation method. This is a time-spatial correlation of amplitudes with a train of time-shifted Dirac delta functions through different hyperbolic paths. It is tested on a two-layer horizontal model, three-layer model which includes a dipping layer (with and without lateral heterogeneity) and also on synthetic Marmousi model data sets. Synthetic tests show that the introduced method is able to reconstruct reflection events at the correct time-offset positions which are hidden in results obtained by the general cross-correlation method. Also, a depth migrated section shows a good match between imaged-horizons and the true model. It is possible to generate off-end virtual gathers by applying the method to a linear array of receivers and to construct a velocity model by semblance velocity analysis of individually extracted gathers.</p>

2021 ◽  
Vol 18 (2) ◽  
pp. 1-15
Author(s):  
Diako Hariri Naghadeh ◽  
Christopher J Bean ◽  
Florent Brenguier ◽  
Patrick J Smith

Abstract Since explosive and impulsive seismic sources such as dynamite, air guns, gas guns or even vibroseis can have a big impact on the environment, some companies have decided to record ambient seismic noise and use it to estimate the physical properties of the subsurface. Big challenges arise when the aim is extracting body waves from recorded passive signals, especially in the presence of strong surface waves. In passive seismic signals, such body waves are usually weak in comparison to surface waves that are much more prominent. To understand the characteristics of passive signals and the effect of natural source locations, three simple synthetic models were created. To extract body waves from simulated passive signals we propose and test a Radon-correlation method. This is a time-spatial correlation of amplitudes with a train of time-shifted Dirac delta functions through different hyperbolic paths. It is tested on a two-layer horizontal model, a three-layer model that includes a dipping layer (with and without lateral heterogeneity) and also on synthetic Marmousi model data sets. Synthetic tests show that the introduced method is able to reconstruct reflection events at the correct time-offset positions that are hidden in results obtained by the general cross-correlation method. Also, a depth migrated section shows a good match between imaged horizons and the true model. It is possible to generate off-end virtual gathers by applying the method to a linear array of receivers and to construct a velocity model by semblance velocity analysis of individually extracted gathers.


Geophysics ◽  
2017 ◽  
Vol 82 (3) ◽  
pp. S197-S205 ◽  
Author(s):  
Zhaolun Liu ◽  
Abdullah AlTheyab ◽  
Sherif M. Hanafy ◽  
Gerard Schuster

We have developed a methodology for detecting the presence of near-surface heterogeneities by naturally migrating backscattered surface waves in controlled-source data. The near-surface heterogeneities must be located within a depth of approximately one-third the dominant wavelength [Formula: see text] of the strong surface-wave arrivals. This natural migration method does not require knowledge of the near-surface phase-velocity distribution because it uses the recorded data to approximate the Green’s functions for migration. Prior to migration, the backscattered data are separated from the original records, and the band-passed filtered data are migrated to give an estimate of the migration image at a depth of approximately one-third [Formula: see text]. Each band-passed data set gives a migration image at a different depth. Results with synthetic data and field data recorded over known faults validate the effectiveness of this method. Migrating the surface waves in recorded 2D and 3D data sets accurately reveals the locations of known faults. The limitation of this method is that it requires a dense array of receivers with a geophone interval less than approximately one-half [Formula: see text].


Author(s):  
Julie Albaric ◽  
Daniela Kühn ◽  
Matthias Ohrnberger ◽  
Nadège Langet ◽  
Dave Harris ◽  
...  

Abstract We analyze data from passive and active seismic experiments conducted in the Adventdalen valley of Svalbard in the Norwegian Arctic. Our objective is to characterize the ambient wavefield of the region and to investigate permafrost dynamics through estimates of seismic velocity variations. We are motivated by a need for early geophysical detection of potentially hazardous changes to permafrost stability. We draw upon several data sources to constrain various aspects of seismic wave propagation in Adventdalen. We use f-k analysis of five years of continuous data from the Spitsbergen seismic array (SPITS) to demonstrate that ambient seismic noise on Svalbard consists of continuously present body waves and intermittent surface waves appearing at regular intervals. A change in wavefield direction accompanies the sudden onset of surface waves when the average temperature rises above the freezing point, suggesting a cryogenic origin. This hypothesis is supported further by our analysis of records from a temporary broadband network, which indicates that the background wavefield is dominated by icequakes. Synthetic Green’s functions calculated from a 3D velocity model match well with empirical Green’s functions constructed from the recorded ambient seismic noise. We use a shallow shear-wave velocity model, obtained from active seismic measurements, to estimate the maximum depth of Rayleigh wave sensitivity to changes in shear velocity to be in the 50–100 m range. We extract seasonal variations in seismic velocities from ambient noise cross-correlation functions computed over three years of SPITS data. We attribute relative velocity variations to changes in the ice content of the shallow (2–4 m depth) permafrost, which is sensitive to seasonal temperature changes. A linear decreasing trend in seismic velocity is observed over the years, most likely due to permafrost warming.


Geophysics ◽  
2011 ◽  
Vol 76 (6) ◽  
pp. V115-V128 ◽  
Author(s):  
Ning Wu ◽  
Yue Li ◽  
Baojun Yang

To remove surface waves from seismic records while preserving other seismic events of interest, we introduced a transform and a filter based on recent developments in image processing. The transform can be seen as a weighted Radon transform, in particular along linear trajectories. The weights in the transform are data dependent and designed to introduce large amplitude differences between surface waves and other events such that surface waves could be separated by a simple amplitude threshold. This is a key property of the filter and distinguishes this approach from others, such as conventional ones that use information on moveout ranges to apply a mask in the transform domain. Initial experiments with synthetic records and field data have demonstrated that, with the appropriate parameters, the proposed trace transform filter performs better both in terms of surface wave attenuation and reflected signal preservation than the conventional methods. Further experiments on larger data sets are needed to fully assess the method.


Geophysics ◽  
2021 ◽  
pp. 1-84
Author(s):  
Chunying Yang ◽  
Wenchuang Wang

Irregular acquisition geometry causes discontinuities in the appearance of surface wave events, and a large offset causes seismic records to appear as aliased surface waves. The conventional method of sampling data affects the accuracy of the dispersion spectrum and reduces the resolution of surface waves. At the same time, ”mode kissing” of the low-velocity layer and inhomogeneous scatterers requires a high-resolution method for calculating surface wave dispersion. This study tested the use of the multiple signal classification (MUSIC) algorithm in 3D multichannel and aliased wavefield separation. Azimuthal MUSIC is a useful method to estimate the phase velocity spectrum of aliased surface wave data, and it represent the dispersion spectra of low-velocity and inhomogeneous models. The results of this study demonstrate that mode-kissing affects dispersion imaging, and inhomogeneous scatterers change the direction of surface-wave propagation. Surface waves generated from the new propagation directions are also dispersive. The scattered surface wave has a new dispersion pattern different to that of the entire record. Diagonal loading was introduced to improve the robustness of azimuthal MUSIC, and numerical experiments demonstrate the resultant effectiveness of imaging aliasing surface waves. A phase-matched filter was applied to the results of azimuthal MUSIC, and phase iterations were unwrapped in a fast and stable manner. Aliased surface waves and body waves were separated during this process. Overall, field data demonstrate that azimuthal MUSIC and phase-matched filters can successfully separate aliased surface waves.


1995 ◽  
Vol 85 (3) ◽  
pp. 716-735 ◽  
Author(s):  
John F. Cassidy ◽  
Garry C. Rogers

Abstract On 6 April 1992, a magnitude 6.8 (MS) earthquake occurred in the triple-junction region at the northern end of the Cascadia subduction zone. This was the largest earthquake in at least 75 yr to occur along the 110-km-long Revere-Dellwood-Wilson (RDW) transform fault and the first large earthquake in this region recorded by modern broadband digital seismic networks. It thus provides an opportunity to examine the rupture process along a young (<2 Ma) oceanic transform fault and to gain better insight into the tectonics of this triple-junction region. We have investigated the source parameters and the rupture process of this earthquake by modeling broadband body waves and long-period surface waves and by accurately locating the mainshock and the first 10 days of aftershocks using a well-located “calibration” event recorded during an ocean-bottom seismometer survey. Analysis of P and SH waveforms reveals that this was a complex rupture sequence consisting of three strike-slip subevents in 12 sec. The initial rupture occurred 5 to 6 km to the SW of the seafloor trace of the RDW fault at 50.55° N, 130.46° W. The dominant subevent occurred 2 to 3 sec later and 4.3 km beneath the seafloor trace of the RDW fault, and a third subevent occurred 5 sec later, 18 km to the NNW, suggesting a northwestward propagating rupture. The aftershock sequence extended along a 60- to 70-km-long segment of the RDW fault, with the bulk of the activity concentrated ∼30 to 40 km to the NNW of the epicenter, consistent with this interpretation. The well-constrained mechanism of the initial rupture (strike/dip/slip 339°/90°/−168°) and of the largest aftershock (165°/80°/170°) are rotated 15° to 20° clockwise relative to the seafloor trace of the RDW fault but are parallel to the Pacific/North America relative plate motion vector. In contrast, the mechanisms of the dominant subevent (326°/87°/−172°), and the long-period solution derived from surface waves aligns with the RDW fault. This suggests that small earthquakes (M < 6) in this area occur along faults that are optimally aligned with respect to the regional stress field, whereas large earthquakes, involving tens of kilometers of rupture, activate the RDW fault. For the mainshock, we estimate a seismic moment (from surface waves) of 1.0 × 1026 dyne-cm, a stress drop of 60 bars, and an average slip of 1.2 m. This represents only 21 yr of strain accumulation, implying that there is either a significant amount of aseismic slip along the RDW fault or that much of the strain accumulation manifests itself as deformation within the Dellwood and Winona blocks or along the continental margin.


Author(s):  
L. Borcea ◽  
J. Garnier ◽  
K. Sølna

We derive a radiative transfer equation that accounts for coupling from surface waves to body waves and the other way around. The model is the acoustic wave equation in a two-dimensional waveguide with reflecting boundary. The waveguide has a thin, weakly randomly heterogeneous layer near the top surface, and a thick homogeneous layer beneath it. There are two types of modes that propagate along the axis of the waveguide: those that are almost trapped in the thin layer, and thus model surface waves, and those that penetrate deep in the waveguide, and thus model body waves. The remaining modes are evanescent waves. We introduce a mathematical theory of mode coupling induced by scattering in the thin layer, and derive a radiative transfer equation which quantifies the mean mode power exchange. We study the solution of this equation in the asymptotic limit of infinite width of the waveguide. The main result is a quantification of the rate of convergence of the mean mode powers toward equipartition.


1967 ◽  
Vol 57 (1) ◽  
pp. 55-81
Author(s):  
E. J. Douze

abstract This report consists of a summary of the studies conducted on the subject of short-period (6.0-0.3 sec period) noise over a period of approximately three years. Information from deep-hole and surface arrays was used in an attempt to determine the types of waves of which the noise is composed. The theoretical behavior of higher-mode Rayleigh waves and of body waves as measured by surface and deep-hole arrays is described. Both surface and body waves are shown to exist in the noise. Surface waves generally predominate at the longer periods (of the period range discussed) while body waves appear at the shorter periods at quiet sites. Not all the data could be interpreted to define the wave types present.


Author(s):  
Tianshi Liu ◽  
Haiming Zhang

The cross-correlations of ambient noise or earthquake codas are massively used in seismic tomography to measure the dispersion curves of surface waves and the travel times of body waves. Such measurements are based on the assumption that these kinematic parameters in the cross-correlations of noise coincide with those in Green's functions. However, the relation between the cross-correlations of noise and Green's functions deserves to be studied more precisely. In this paper, we use the asymptotic analysis to study the dispersion relations of surface waves and the travel times of body waves, and come to the conclusion that for the spherically symmetric Earth model, when the distribution of noise sources is laterally uniform, the dispersion relations of surface waves and the travel times of SH body-wave phases in noise correlations should be exactly the same as those in Green's functions.


Sign in / Sign up

Export Citation Format

Share Document