skip to main content


Title: Rayleigh-wave multicomponent crosscorrelation-based source strength distribution inversions. Part 2: a workflow for field seismic data
SUMMARY Estimation of ambient seismic source distributions (e.g. location and strength) can aid studies of seismic source mechanisms and subsurface structure investigations. One can invert for the ambient seismic (noise) source distribution by applying full-waveform inversion (FWI) theory to seismic (noise) crosscorrelations. This estimation method is especially applicable for seismic recordings without obvious body-wave arrivals. Data pre-processing procedures are needed before the inversion, but some pre-processing procedures commonly used in ambient noise tomography can bias the ambient (noise) source distribution estimation and should not be used in FWI. Taking this into account, we propose a complete workflow from the raw seismic noise recording through pre-processing procedures to the inversion. We present the workflow with a field data example in Hartoušov, Czech Republic, where the seismic sources are CO2 degassing areas at Earth’s surface (i.e. a fumarole or mofette). We discuss factors in the processing and inversion that can bias the estimations, such as inaccurate velocity model, anelasticity and array sensitivity. The proposed workflow can work for multicomponent data across different scales of field data.  more » « less
Award ID(s):
1643795
NSF-PAR ID:
10203459
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geophysical Journal International
Volume:
222
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
2084 to 2101
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Both natural and anthropogenic seismic sources generate so‐called ambient seismic waves. One in turn can use ambient seismic waves to estimate these source distributions and study source characteristics, for instance the source mechanism. A commonly used estimation method is called matched field processing (MFP), and the MFP results are inherently smeared by the array geometry. Another approach to estimate ambient seismic sources is to apply full waveform inversion (FWI) to the crosscorrelations of ambient seismic wave recordings. Both methods have pros and cons, but the model resolution and uncertainty in these two methods are important for the interpretation. Unfortunately, this topic has attracted little attention in the past. We propose to estimate both the model resolution matrix and model covariance matrix of the inversion using singular value decomposition. We demonstrate our estimates using two examples, one of which is an actual field array geometry. We quantitatively compare the model resolution of the two methods and discuss the model null space. We demonstrate that FWI has superior resolution with enough independent data and should be used when computational resources permit.

     
    more » « less
  2. Seismic imaging and monitoring of the near-surface structure are crucial for the sustainable development of urban areas. However, standard seismic surveys based on cabled or autonomous geophone arrays are expensive and hard to adapt to noisy metropolitan environments. Distributed acoustic sensing (DAS) with pre-existing telecom fiber optic cables, together with seismic ambient noise interferometry, have the potential to fulfill this gap. However, a detailed noise wavefield characterization is needed before retrievingcoherent waves from chaotic noise sources. We analyze local seismic ambient noise by tracking five-month changes in signal-to-noise ratio (SNR) of Rayleigh surface wave estimated from traffic noise recorded by DAS along the straight university campus busy road. We apply the seismic interferometry method to the 800 m long part of the Penn State Fiber-Optic For Environment Sensing (FORESEE) array. We evaluate the 160 virtual shot gathers (VSGs) by determining the SNR using the slant-stack technique. We observe strong SNR variations in time and space. We notice higher SNR for virtual source points close to road obstacles. The spatial noise distribution confirms that noise energy focuses mainly on bumps and utility holes. We also see the destructive impact of precipitation, pedestrian traffic, and traffic along main intersections on VSGs. A similar processing workflow can be applied to various straight roadside fiber optic arrays in metropolitan areas.

     
    more » « less
  3. SUMMARY

    Ambient noise tomography is a well-established tomographic imaging technique but the effect that spatially variable noise sources have on the measurements remains challenging to account for. Full waveform ambient noise inversion has emerged recently as a promising solution but is computationally challenging since even distant noise sources can have an influence on the interstation correlation functions and therefore requires a prohibitively large numerical domain, beyond that of the tomographic region of interest. We investigate a new strategy that allows us to reduce the simulation domain while still being able to account for distant contributions. To allow nearby numerical sources to account for distant true sources, we introduce correlated sources and generate a time-dependent effective source distribution at the boundary of a small region of interest that excites the correlation wavefield of a larger domain. In a series of 2-D numerical simulations, we demonstrate that the proposed methodology with correlated sources is able to successfully represent a far-field source that is simultaneously present with nearby sources and the methodology also successfully results in a robustly estimated noise source distribution. Furthermore, we show how beamforming results can be used as prior information regarding the azimuthal variation of the ambient noise sources in helping determine the far-field noise distribution. These experiments provide insight into how to reduce the computational cost needed to perform full waveform ambient noise inversion, which is key to turning it into a viable tomographic technique. In addition, the presented experiments may help reduce source-induced bias in time-dependent monitoring applications.

     
    more » « less
  4. SUMMARY

    Infrasound sensors are deployed in a variety of spatial configurations and scales for geophysical monitoring, including networks of single sensors and networks of multisensor infrasound arrays. Infrasound signal detection strategies exploiting these data commonly make use of intersensor correlation and coherence (array processing, multichannel correlation); network-based tracking of signal features (e.g. reverse time migration); or a combination of these such as backazimuth cross-bearings for multiple arrays. Single-sensor trace-based denoising techniques offer significant potential to improve all of these various infrasound data processing strategies, but have not previously been investigated in detail. Single-sensor denoising represents a pre-processing step that could reduce the effects of ambient infrasound and wind noise in infrasound signal association and location workflows. We systematically investigate the utility of a range of single-sensor denoising methods for infrasound data processing, including noise gating, non-negative matrix factorization, and data-adaptive Wiener filtering. For the data testbed, we use the relatively dense regional infrasound network in Alaska, which records a high rate of volcanic eruptions with signals varying in power, duration, and waveform and spectral character. We primarily use data from the 2016–2017 Bogoslof volcanic eruption, which included multiple explosions, and synthetics. The Bogoslof volcanic sequence provides an opportunity to investigate regional infrasound detection, association, and location for a set of real sources with varying source spectra subject to anisotropic atmospheric propagation and varying noise levels (both incoherent wind noise and coherent ambient infrasound, primarily microbaroms). We illustrate the advantages and disadvantages of the different denoising methods in categories such as event detection, waveform distortion, the need for manual data labelling, and computational cost. For all approaches, denoising generally performs better for signals with higher signal-to-noise ratios and with less spectral and temporal overlap between signals and noise. Microbaroms are the most globally pervasive and repetitive coherent ambient infrasound noise source, with such noise often referred to as clutter or interference. We find that denoising offers significant potential for microbarom clutter reduction. Single-channel denoising of microbaroms prior to standard array processing enhances both the quantity and bandwidth of detectable volcanic events. We find that reduction of incoherent wind noise is more challenging using the denoising methods we investigate; thus, station hardware (wind noise reduction systems) and site selection remain critical and cannot be replaced by currently available digital denoising methodologies. Overall, we find that adding single-channel denoising as a component in the processing workflow can benefit a variety of infrasound signal detection, association, and location schemes. The denoising methods can also isolate the noise itself, with utility in statistically characterizing ambient infrasound noise.

     
    more » « less
  5. We develop a new 3D ambient noise tomography (3D ANT) method for geotechnical site characterization. It requires recording ambient noise fields using a 2D surface array of geophones, from which experimental crosscorrelation functions (CCFs) are then extracted and directly inverted to obtain an S-wave velocity ([Formula: see text]) structure. The method consists of a forward simulation using 3D P-SV elastic wave equations to compute the synthetic CCF and an adjoint-state inversion to match the synthetic CCFs to the experimental CCFs for extraction of [Formula: see text]. The main advantage of the presented method, as compared with conventional passive-source seismic methods using characteristics of Green’s function (GF), is that it does not require equal energy on both sides of each receiver pair or far-field wavefields to retrieve the true GF. Instead, the source power spectrum density is inverted during the analysis and incorporated into the forward simulation of the synthetic CCFs to account for source energy distribution. After testing on synthetic data, the 3D ANT method is applied to 3 h of ambient noise recordings at the Garner Valley Downhole Array (GVDA) site in California, using a surface array of 196 geophones placed on a 14 × 14 grid with 5 m spacing. The inverted 3D [Formula: see text] model is found to be consistent with previous invasive and noninvasive geotechnical characterization efforts at the GVDA site. 
    more » « less