Noise Suppression Using a Near-Source Wavelet

Geophysics ◽  
2021 ◽  
pp. 1-56
Author(s):  
Afshin Aghayan ◽  
Priyank Jaiswal

Denoising becomes a non-trivial task when noise and signal overlap in multiple domains such as time, frequency, and velocity. Fortunately, signal and noise waveforms in general tend to remain morphologically different. This paper shows how morphological differences can be used to separate body-wave signals from other waveforms such as ground roll and cultural noise. The key was finding a wavelet that was a close approximation of the true source signature (SS) and remained uncontaminated by the Green’s function in any significant manner. An inverse filter designed using such a wavelet selectively compressed the body waves which was then extracted using median and low-pass filters. The overall phenomenon is explained with a synthetic example. The idea is also tested on a land dataset that was generated using a large weight drop source where a wavelet recorded ∼3 m from the source location fulfilled the criteria set in the proposed method. Results suggest that the incremental effort of recording an extra trace close to the source location during acquisition may provide previously unavailable denoising opportunities during processing although the trace itself may be redundant for imaging.

Geophysics ◽  
2008 ◽  
Vol 73 (4) ◽  
pp. D17-D33 ◽  
Author(s):  
Bing Zhou ◽  
Stewart Greenhalgh ◽  
Alan Green

Crosshole seismic tomography often is applied to image the velocity structure of an interwell medium. If the rocks are anisotropic, the tomographic technique must be adapted to the complex situation; otherwise, it leads to a false interpretation. We propose a nonlinear kinematic inversion method for crosshole seismic tomography in composite transversely isotropic media with known dipping symmetry axes. This method is based on a new version of the first-order traveltime perturbation equation. It directly uses the derivative of the phase velocity rather than the eigenvectors of the body-wave modes to overcome the singularity problem for application to the two quasi-shear waves. We applied an iterative nonlinear solver incorporating our kinematic ray-tracing scheme and directly compute the Jacobian matrix in an arbitrary reference medium. This reconstructs the five elastic moduli or Thomsen parameters from the first-arrival traveltimes of the three seismic body waves (qP, qSV, qSH) in strongly and weakly anisotropic media. We conducted three synthetic experiments that involve determining anisotropic parameters for a homogeneous rock, reconstructing a fault embedded in a strongly anisotropic background, and imaging a complicated four-layer model containing a small channel and a buried dipping interface. We compared results of our nonlinear inversion method with isotropic tomography and the traditional linear anisotropic inversion scheme, which showed the capability and superiority of the new scheme for crosshole tomographic imaging.


1999 ◽  
Vol 202 (23) ◽  
pp. 3423-3430 ◽  
Author(s):  
J.J. Videler ◽  
U.K. Muller ◽  
E.J. Stamhuis

Vertebrates swimming with undulations of the body and tail have inflection points where the curvature of the body changes from concave to convex or vice versa. These inflection points travel down the body at the speed of the running wave of bending. In movements with increasing amplitudes, the body rotates around the inflection points, inducing semicircular flows in the adjacent water on both sides of the body that together form proto-vortices. Like the inflection points, the proto-vortices travel towards the end of the tail. In the experiments described here, the phase relationship between the tailbeat cycle and the inflection point cycle can be used as a first approximation of the phase between the proto-vortex and the tailbeat cycle. Proto-vortices are shed at the tail as body vortices at roughly the same time as the inflection points reach the tail tip. Thus, the phase between proto-vortex shedding and tailbeat cycle determines the interaction between body and tail vortices, which are shed when the tail changes direction. The shape of the body wave is under the control of the fish and determines the position of vortex shedding relative to the mean path of motion. This, in turn, determines whether and how the body vortex interacts with the tail vortex. The shape of the wake and the contribution of the body to thrust depend on this interaction between body vortex and tail vortex. So far, we have been able to describe two types of wake. One has two vortices per tailbeat where each vortex consists of a tail vortex enhanced by a body vortex. A second, more variable, type of wake has four vortices per tailbeat: two tail vortices and two body vortices shed from the tail tip while it is moving from one extreme position to the next. The function of the second type is still enigmatic.


1976 ◽  
Vol 66 (5) ◽  
pp. 1485-1499 ◽  
Author(s):  
L. J. Burdick ◽  
George R. Mellman

abstract The generalized linear inverse technique has been adapted to the problem of determining an earthquake source model from body-wave data. The technique has been successfully applied to the Borrego Mountain earthquake of April 9, 1968. Synthetic seismograms computed from the resulting model match in close detail the first 25 sec of long-period seismograms from a wide range of azimuths. The main shock source-time function has been determined by a new simultaneous short period-long period deconvolution technique as well as by the inversion technique. The duration and shape of this time function indicate that most of the body-wave energy was radiated from a surface with effective radius of only 8 km. This is much smaller than the total surface rupture length or the length of the aftershock zone. Along with the moment determination of Mo = 11.2 ×1025 dyne-cm, this radius implies a high stress drop of about 96 bars. Evidence in the amplitude data indicates that the polarization angle of shear waves is very sensitive to lateral structure.


1980 ◽  
Vol 70 (3) ◽  
pp. 791-808
Author(s):  
Anne Suteau ◽  
Louis Martel

abstract The transmitted field due to surface waves incident on a local irregularity of a plane-layered medium has been studied. A perturbation method to the first order and the Born approximation can be used if the variations in the thickness of the layers are sufficiently smooth and the wavelengths are long when compared to the size of the irregularities. The spectrum of the perturbed part of the displacement field at the surface is a sum over the surface-wave modes for the regular medium, with an additional term involving the scattered body waves. Numerical computations have been performed for structures composed of a layer overlying a half-space. The contribution of the various modes to the transmitted Love or Rayleigh fields has been studied for several structures. A general method has been obtained to analyze the effect of a complex structure as the superposition of the fields due to simpler ones. When the layer thickness is kept unchanged, the incident mode is not perturbed to the first order. Synthetic seismograms, computed at stations sufficiently close to the irregular region, show how the perturbation of the signal depends on distance. A comparison has been made for Love waves with a finite element method. Both methods give very similar results when the stations are not too close to the irregularities so that the body-wave contribution is negligible. The local phase velocity shows departures from the curves for a regular model.


1984 ◽  
Vol 74 (2) ◽  
pp. 361-376
Author(s):  
John Boatwright ◽  
Jon B. Fletcher

Abstract Seventy-three digitally recorded body waves from nine multiply recorded small earthquakes in Monticello, South Carolina, are analyzed to estimate the energy radiated in P and S waves. Assuming Qα = Qβ = 300, the body-wave spectra are corrected for attenuation in the frequency domain, and the velocity power spectra are integrated over frequency to estimate the radiated energy flux. Focal mechanisms determined for the events by fitting the observed displacement pulse areas are used to correct for the radiation patterns. Averaging the results from the nine events gives 27.3 ± 3.3 for the ratio of the S-wave energy to the P-wave energy using 0.5 〈Fi〉 as a lower bound for the radiation pattern corrections, and 23.7 ± 3.0 using no correction for the focal mechanisms. The average shift between the P-wave corner frequency and the S-wave corner frequency, 1.24 ± 0.22, gives the ratio 13.7 ± 7.3. The substantially higher values obtained from the integral technique implies that the P waves in this data set are depleted in energy relative to the S waves. Cursory inspection of the body-wave arrivals suggests that this enervation results from an anomalous site response at two of the stations. Using the ratio of the P-wave moments to the S-wave moments to correct the two integral estimates gives 16.7 and 14.4 for the ratio of the S-wave energy to the P-wave energy.


1982 ◽  
Vol 72 (4) ◽  
pp. 1093-1109
Author(s):  
Jeffrey W. Given ◽  
Terry C. Wallace ◽  
Hiroo Kanamori

abstract The source mechanisms of the three largest events of the 1980 Mammoth Lakes earthquake sequence have been determined using surface waves recorded on the global digital seismograph network and the long-period body waves recorded on the WWSSN network. Although the fault-plane solutions from local data (Cramer and Toppozada, 1980; Ryall and Ryall, 1981) suggest nearly pure left-lateral strike-slip on north-south planes, the teleseismic waveforms require a mechanism with oblique slip. The first event (25 May 1980, 16h 33m 44s) has a mechanism with a strike of N12°E, dip of 50°E, and a rake of −35°. The second event (27 May 19h 44m 51s) has a mechanism with a strike of N15°E, dip of 50°, and a slip of −11°. The third event (27 May, 14h 50m 57s) has a mechanism with a strike of N22°E, dip of 50°, and a rake of −28°. The first event is the largest and has a moment of 2.9 × 1025 dyne-cm. The second and third events have moments of 1.3 and 1.1 × 1025 dyne-cm, respectively. The body- and surface-wave moments for the first and third events agree closely while for the second event the body-wave moment (approximately 0.6 × 1025 dyne-cm) is almost a factor of 3 smaller than the surface-wave moment. The principal axes of extension of all three events is in the approximate direction of N65°E which agrees with the structural trends apparent along the eastern front of the Sierra Nevada.


1978 ◽  
Vol 68 (5) ◽  
pp. 1281-1292
Author(s):  
John E. Ebel ◽  
L. J. Burdick ◽  
Gordon S. Stewart

abstract The El Golfo earthquake of August 7, 1966 (mb = 6.3, MS = 6.3) occurred near the mouth of the Colorado River at the northern end of the Gulf of California. Synthetic seismograms for this event were computed for both the body waves and the surface waves to determine the source parameters of the earthquake. The body-wave model indicated the source was a right lateral, strike-slip source with a depth of 10 km and a far-field time function 4 sec in duration. The body-wave moment was computed to be 5.0 × 1025 dyne-cm. The surface-wave radiation pattern was found to be consistent with that of the body waves with a surface-wave moment of 6.5 × 1025 dyne-cm. The agreement of the two different moments indicates that the earthquake had a simple source about 4 sec long. A comparison of this earthquake source with the Borrego Mountain and Truckee events demonstrates that all three of these earthquakes behaved as high stress-drop events. El Golfo was shown to be different from the low stress-drop, plate-boundary events which were located on the Gibbs fracture zone in 1967 and 1974.


2019 ◽  
Vol 110 (1) ◽  
pp. 110-126
Author(s):  
Leiph Preston ◽  
Christian Poppeliers ◽  
David J. Schodt

ABSTRACT As a part of the series of Source Physics Experiments (SPE) conducted on the Nevada National Security Site in southern Nevada, we have developed a local-to-regional scale seismic velocity model of the site and surrounding area. Accurate earth models are critical for modeling sources like the SPE to investigate the role of earth structure on the propagation and scattering of seismic waves. We combine seismic body waves, surface waves, and gravity data in a joint inversion procedure to solve for the optimal 3D seismic compressional and shear-wave velocity structures and earthquake locations subject to model smoothness constraints. Earthquakes, which are relocated as part of the inversion, provide P- and S-body-wave absolute and differential travel times. Active source experiments in the region augment this dataset with P-body-wave absolute times and surface-wave dispersion data. Dense ground-based gravity observations and surface-wave dispersion derived from ambient noise in the region fill in many areas where body-wave data are sparse. In general, the top 1–2 km of the surface is relatively poorly sampled by the body waves alone. However, the addition of gravity and surface waves to the body-wave dataset greatly enhances structural resolvability in the near surface. We discuss the methodology we developed for simultaneous inversion of these disparate data types and briefly describe results of the inversion in the context of previous work in the region.


2020 ◽  
Author(s):  
Boris Boullenger ◽  
Merijn de Bakker ◽  
Arie Verdel ◽  
Stefan Carpentier

<p>The theory of ambient seismic noise interferometry offers techniques to retrieve estimates of inter-receiver responses from continuously recorded ambient seismic noise. This is usually achieved by correlating and stacking successive noise panels over sufficiently long periods of time. If the noise panels contain significant body-wave energy, the stacked correlations expected to result in retrieved estimates of the body-wave responses, including reflections. Such application combined with a dense surface seismic array is promising for imaging the subsurface structures at lower cost and lower environmental impact as compared to with controlled seismic sources. Subsequently, this technique can be an alternative to active-source surveys in a range of challenging scenarios and locations, and can also be used to perform time-lapse subsurface characterization.</p><p>In this study, we apply seismic body-wave noise interferometry to 30-days of continuous records from a surface line of 31 receivers spaced by 25 meters in the South of the Netherlands with the aim to image subsurface reflectors, at depths from a few hundreds of meters to a few kilometers. As a first step, we compute stacked auto-correlations and compare the retrieved zero-offset section with a co-located stacked section from a past active reflection survey on the site.</p><p>Yet, the retrieval of reflectivity estimates relies on the identification and collection of a sufficient number of noise panels with recorded body waves that have travelled from the subsurface towards the array. Even in the case of favorable body-wave noise conditions, the panels are most often contaminated with stronger anthropogenic coherent seismic noise, mainly in the form of surface waves, which in turn prevents the stacked correlations to reveal reflectivity. Because of the limited effect of frequency filtering, the application of seismic body-wave noise interferometry requires in fact extensive effort to identify noise panels without prominent coherent noise from the surface activity. Typically, this leads to disregard a significant amount of actually useful data.</p><p>For this reason, we designed, trained and tested a deep convolutional neural network to perform this classification task more efficiently and facilitate the repetition of the retrieval method over long periods of time. We tested several supervised learning schemes to classify the panels, where two classes are defined, according to the presence or absence of prominent coherent noise. The retained classification models achieved close to 90% of prediction accuracy on the test set.</p><p>We used the trained classification models to correlate and stack panels which were predicted in the class with coherent noise absent. The resulting stacked correlations exhibit potential reflectors in a larger depth range than previously achieved. The results show the benefits of using machine learning to collect efficiently a maximum amount of favorable noise panels and a way forward to the upscaling of seismic body-wave noise interferometry for reflectivity imaging.</p>


2020 ◽  
Author(s):  
Hrvoje Tkalčić ◽  
Sheng Wang ◽  
Thanh Son Pham

<p>We have recently shown that all features in the earthquake-coda correlogram can be explained by the similarity of seismic phases that have a common slowness for the analysed receiver pair. This includes both the features that have their equivalents in the conventional traveltime stacks, but also those that were previously unexplained. Consequently, the information contained in the correlograms – cross-correlated ground-motion time-series in a two-dimensional representation – can be used to constrain Earth’s internal structure, however, that requires a proof of concept and further investigation into the origin of the correlation wavefield. We thus first decompose relevant correlogram features into discrete constituents with respect to their arrival times and we uniquely identify contributing seismic phases to each constituent. This confirms that the correlation wavefield does not arise due to the reconstruction of body waves between the two receivers (a.k.a. Green’s function) – instead, it is dominated by the interaction of various body waves, and its features are characterised by complex sensitivity kernels.</p><p>We demonstrate that the event locations relative to the receivers alter the similarities between the body waves, and may result in significant waveform distortions and inaccuracies in arrival-time predictions. We further show that the nature of source-mechanism and energy-release dynamics are the key influencers responsible for individual correlograms equal in quality to a stack of hundreds of correlograms. In other words, a single seismic event that meets a set of criteria in the presence of multiple receivers can completely `illuminate’ the Earth’s interior. Quantitative kernel-decomposition and identification of body-wave pairs that contribute to a given feature in the correlogram, along with informed choices of seismic events, thus makes the correlation-wavefield tomography and other applications fully feasible. This has the potential to change the course of global seismology in the coming decades.</p>


Sign in / Sign up

Export Citation Format

Share Document