Title: Source Independent Velocity Recovery Using Imaginary FWI
We consider the problem of velocity inversion/calibration in passive survey, where the seismic source is also an unknown. In earthquake detection or microseismic localization, the major task is to reconstruct the passive seismic sources, but due to the source-velocity coupling, source reconstructions are inherently affected by inaccurate knowledge of the velocity, bringing the need of velocity calibration. We propose a source independent velocity calibration method that recovers the velocity without the source information, thus providing a better ground for source inversion. Unlike existing methods that assume sources to be a linear combination of separated point sources, the proposed method allows sources to lie on a line singularity (representing rock cracks), as long as the activation time is relatively brief. The proposed approach is based on the observation that the spatial distribution of the source is separable from the velocity model after a proper Helmholtz domain projection.  more » « less
Award ID(s):
1909523 2006881 1715178
PAR ID:
10308763
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
82nd EAGE Annual Conference & Exhibition
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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
  2. Abstract Environmental seismology is a rapidly growing field that is providing new insights into many surficial and near-surface Earth processes, largely through passive detection of events (e.g., landslides and sediment movement) and subsurface velocity changes that reflect changes in saturation (e.g., water table changes). Recent advances in instrumentation (e.g., distributed acoustic sensing) and analysis techniques (e.g., machine learning and artificial intelligence) are expanding the use cases and data volumes of passive environmental seismology. However, despite these successes, passive environmental seismology has two disadvantages that may limit the usefulness of the data. First, passive-source studies that rely on surface-wave observations usually produce S-wave velocity without accompanying P-wave velocity measurements, thus incompletely characterizing subsurface elastic properties. Second, passive-source studies often do not provide subsurface structure at the resolution necessary to distinguish the near-surface units that host and control environmental processes (e.g., soil, saprolite, and fractured bedrock). We suggest that these limitations can be surmounted by incorporating active-source seismic imaging of the critical zone (CZ). Active-source data records the entire wavefield, thus providing detailed P- and S-velocity structure that can contextualize passive-source measurements. Active-source CZ seismology is a rapidly changing field, with notable recent advances in instrumentation (e.g., small, autonomous nodes) and analysis techniques (e.g., wider adoption of full-waveform inversion). As a result, new types of surveys (e.g., three-component, 3D, and time-lapse) are imaging the CZ in novel ways and with unprecedented resolution. Here, we review the state of the art in active-source CZ seismology and provide several examples demonstrating that active-source environmental seismology can achieve meter-scale resolution of CZ structure and processes. 
    more » « less
  3. ABSTRACT Distributed acoustic sensing (DAS) technology is an emerging field of seismic sensing that enables recording ambient noise seismic data along the entire length of a fiber-optic cable at meter-scale resolution. Such a dense spatial resolution of recordings over long distances has not been possible using traditional methods because of limited hardware resources and logistical concerns in an urban environment. The low spatial resolution of traditional passive seismic acquisition techniques has limited the accuracy of the previously generated velocity profiles in many important urban regions, including the Reno-area basin, to the top 100 m of the underlying subsurface. Applying the method of seismic interferometry to ambient noise strain rate data obtained from a dark-fiber cable allows for generating noise cross correlations, which can be used to infer shallow and deep subsurface properties and basin geometry. We gathered DAS ambient noise seismic data for this study using a 12 km portion of a dark-fiber line in Reno, Nevada. We used gathered data to generate and invert dispersion curves to estimate the near-surface shear-wave velocity structure. Comparing the generated velocity profiles with previous regional studies shows good agreement in determining the average depth to bedrock and velocity variations in the analyzed domain. A synthetic experiment is also performed to verify the proposed framework further and better understand the effect of the infrastructural cover along the cable. The results obtained from this research provide insight into the application of DAS using dark-fiber lines in subsurface characterization in urban environments. It also discusses the potential effects of the conduit that covers such permanent fiber installations on the produced inversion results. 
    more » « less
  4. 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
  5. The likelihood and severity of mining-induced seismicity is dependent on the stress magnitude and distribution within the rockmass. Both geomechanical modelling and seismic imaging provide volumetric information about stress redistribution within a rock mass; however, both require some degree of calibration and/or uncertainly quantification. Typically, geomechanical models are calibrated with point-location geotechnical measurements; however, improved model calibration could be obtained using volumetric information. Seismic imaging, through passive seismic tomography, provides information about stress distribution throughout a rockmass, but some portions of the volume have less reliable results than others. Often, to determine which portions of the volume have reliable results, one starts with a checkerboard model, where synthetic rays are propagated through a model with alternating high and low velocities. The results are then inverted and compared to the initial, synthetic velocity model to determine which portions of the volume can be reliably trusted. The problem with this approach, however, is that it is arbitrary and unrealistically represents the actual rockmass; we suggest that stress information from a geomechanical model could be used instead. This paper presents an approach to use the results from both geomechanical modelling and passive seismic tomography to jointly calibrate and quantify uncertainty in both methods. This approach will benefit the geomechanical modelling effort by not only providing a method for calibration but also providing data about the actual stress distribution within the rockmass, which is affected by geologic inhomogeneity that may not be properly accounted for in the geomechanical model. 
    more » « less