More Like this

1. Rayleigh-wave multicomponent crosscorrelation-based source strength distribution inversions. Part 2: a workflow for field seismic data

   https://doi.org/10.1093/gji/ggaa284  

   Xu, Zongbo; Mikesell, T Dylan; Umlauft, Josefine; Gribler, Gabriel (September 2020, Geophysical Journal International)

   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. Three-dimensional Gauss–Newton constant- Q viscoelastic full-waveform inversion of near-surface seismic wavefields

   https://doi.org/10.1093/gji/ggac287  

   Mirzanejad, Majid; Tran, Khiem T.; Wang, Yao (July 2022, Geophysical Journal International)

   SUMMARY Full-waveform inversion (FWI) methods rely on accurate numerical simulation of wave propagation in the analysed medium. Acoustic or elastic wave equations are often used to model seismic wave propagation. These types of simulations do not account for intrinsic attenuation effects due to material anelasticity, and thus correction techniques have been utilized in practice to partially compensate the anelasticity. These techniques often only consider the waveform amplitude correction based on averaging of overall amplitude response over the entire data set, and ignore the phase correction. Viscoelastic wave equations account for the anelastic response in both waveform amplitude and phase, and are therefore a more suitable alternative. In this study, we present a novel 3-D Gauss–Newton viscoelastic FWI (3-D GN-VFWI) method. To address the main challenge of the Gauss–Newton optimization, we develop formulas to compute the Jacobian efficiently by the convolution of virtual sources and backward wavefields. The virtual sources are obtained by directly differentiating the viscoelastic wave equations with respect to model parameters. In order to resolve complex 3-D structures with reasonable computational effort, a homogeneous attenuation (Q factor) is used throughout the analysis to model the anelastic effects. Synthetic and field experiments are performed to demonstrate the utility of the method. The synthetic results clearly demonstrate the ability of the method in characterizing a challenging velocity profile, including voids and reverse velocity layers. The field experimental results show that method successfully characterizes the complex substructure with two voids and undulating limestone bedrock, which are confirmed by invasive tests. Compared to 3-D elastic FWI results, the presented viscoelastic method produces more accurate results regarding depths of the voids and bedrock. This study suggests that the improvement of imaging accuracy would warrant the widespread use of viscoelastic wave equations in FWI problems. To our best knowledge, this is the first reported study on 3-D GN-VFWI at any scale. This study provides the new theory and formulation for the use of Gauss–Newton optimization on the 3-D viscoelastic problem. 

   more » « less
3. Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation

   https://doi.org/10.1190/geo2021-0248.1  

   Xing, Guangchi; Zhu, Tieyuan (January 2022, GEOPHYSICS)

   We have formulated the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we find that the adjoint wave propagator preserves the dispersion and compensates the amplitude, whereas the time-reversed adjoint wave propagator behaves identically to the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit [Formula: see text] parameterization, which avoids the implicit [Formula: see text] in the conventional viscoacoustic/viscoelastic full-waveform inversion ([Formula: see text]-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. We find that the lossless propagation kernel dominates the velocity kernel, whereas the dissipation kernel dominates the attenuation kernel over the dispersion kernel. After validating the Fréchet kernels using the finite-difference method, we conduct a numerical example to demonstrate the capability of the kernels to characterize the velocity and attenuation anomalies. The kernels of different misfit measurements are presented to investigate their different sensitivities. Our results suggest that, rather than the traveltime, the amplitude and the waveform kernels are more suitable to capture attenuation anomalies. These kernels lay the foundation for the multiparameter inversion with the fractional formulation, and the decoupled nature of them promotes our understanding of the significance of different physical processes in [Formula: see text]-FWI. 

   more » « less
4. Resolution and trade-offs in global anelastic full-waveform inversion

   https://doi.org/10.1093/gji/ggad462  

   Espindola-Carmona, Armando; Örsvuran, Rıdvan; Mai, P. Martin; Bozdağ, Ebru; Peter, Daniel B. (November 2023, Geophysical Journal International)

   SUMMARY Improving the resolution of seismic anelastic models is critical for a better understanding of the Earth’s subsurface structure and dynamics. Seismic attenuation plays a crucial role in estimating water content, partial melting and temperature variations in the Earth’s crust and mantle. However, compared to seismic wave-speed models, seismic attenuation tomography models tend to be less resolved. This is due to the complexity of amplitude measurements and the challenge of isolating the effect of attenuation in the data from other parameters. Physical dispersion caused by attenuation also affects seismic wave speeds, and neglecting scattering/defocusing effects in classical anelastic models can lead to biased results. To overcome these challenges, it is essential to account for the full 3-D complexity of seismic wave propagation. Although various synthetic tests have been conducted to validate anelastic full-waveform inversion (FWI), there is still a lack of understanding regarding the trade-off between elastic and anelastic parameters, as well as the variable influence of different parameter classes on the data. In this context, we present a synthetic study to explore different strategies for global anelastic inversions. To assess the resolution and sensitivity for different misfit functions, we first perform mono-parameter inversions by inverting only for attenuation. Then, to study trade-offs between parameters and resolution, we test two different inversion strategies (simultaneous and sequential) to jointly constrain the elastic and anelastic parameters. We found that a sequential inversion strategy performs better for imaging attenuation than a simultaneous inversion. We also demonstrate the dominance of seismic wave speeds over attenuation, underscoring the importance of determining a good approximation of the Hessian matrix and suitable damping factors for each parameter class. 

   more » « less
5. Sinkhole detection with 3D full seismic waveform tomography

   https://doi.org/10.1190/GEO2019-0490.1  

   Mirzanejad, Majid; Tran, Khiem T.; McVay, Michael; Horhota, David; Wasman, Scott J. (September 2020, GEOPHYSICS)

   Sinkhole collapse may result in significant property damage and even loss of life. Early detection of sinkhole attributes (buried voids, raveling zones) is critical to limit the cost of remediation. One of the most promising ways to obtain subsurface imaging is 3D seismic full-waveform inversion. For demonstration, a recently developed 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method is used to detect buried voids, raveling soils, and characterize variable subsurface soil/rock layering. It is based on a finite-difference solution of 3D elastic wave equations and Gauss-Newton optimization. The method is tested first on a data set constructed from the numerical simulation of a challenging synthetic model and subsequently on field data collected from two separate test sites in Florida. For the field tests, receivers and sources are placed in uniform 2D surface grids to acquire the seismic wavefields, which then are inverted to extract the 3D subsurface velocity structures. The inverted synthetic results suggest that the approach is viable for detecting voids and characterizing layering. The field seismic results reveal that the 3D waveform analysis identified a known manmade void (plastic culvert), unknown natural voids, raveling, as well as laterally variable soil/rock layering including rock pinnacles. The results are confirmed later by standard penetration tests, including depth to bedrock, two buried voids, and a raveling soil zone. Our study provides insight into the application of the 3D seismic FWI technique as a powerful tool in detecting shallow voids and other localized subsurface features. 

   more » « less