More Like this

1. New likelihood functions and level-set prior for Bayesian full-waveform inversion

   https://doi.org/10.1190/segam2020-3428223.1  

   Dunlop, Matt ; Yang, Yunan ( September 2020 , SEG Technical Program Expanded Abstracts 2020)

   null (Ed.)

   Seismic full-waveform inversion aims to reconstruct subsurface medium parameters from recorded seismic data. It is solved as a constrained optimization problem in the deterministic approach. Many different objective functions have been proposed to tackle the nonconvexity that originated from the cycle-skipping issues. The analogy between objective functions in the deterministic inversion and likelihood functions in Bayesian inversion motivates us to analyze the noise model each objective function accounts for under the Bayesian inference setting. We also show the existence and wellposedness of their corresponding posterior measures. In particular, the theorem shows that theWasserstein-type likelihood offers better stability with respect to the noise in the recorded data. Together with an application of the level-set prior, we demonstrate by numerical examples the successful reconstruction from Bayesian full-waveform inversion under the proper choices of the likelihood function and the prior distribution. 

   more » « less
2. 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)

   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
3. An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

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

   Puel, S. ; Becker, T. W. ; Villa, U. ; Ghattas, O. ; Liu, D. ( November 2023 , Geophysical Journal International)

   Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

    

   more » « less
4. 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)

   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
5. Using convolutional neural networks to develop starting models for near-surface 2-D full waveform inversion

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

   Vantassel, Joseph P. ; Kumar, Krishna ; Cox, Brady R. ( May 2022 , Geophysical Journal International)

   Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamentally change near-surface (<30 m) site characterization by enabling the recovery of high-resolution (metre-scale) 2-D/3-D maps of subsurface elastic material properties. Yet, FWI results are quite sensitive to their starting model due to their dependence on local-search optimization techniques and inversion non-uniqueness. Starting model dependence is particularly problematic for near-surface FWI due to the complexity of the recorded seismic wavefield (e.g. dominant surface waves intermixed with body waves) and the potential for significant spatial variability over short distances. In response, convolutional neural networks (CNNs) are investigated as a potential tool for developing starting models for near-surface 2-D elastic FWI. Specifically, 100 000 subsurface models were generated to be representative of a classic near-surface geophysics problem; namely, imaging a two-layer, undulating, soil-over-bedrock interface. A CNN has been developed from these synthetic models that is capable of transforming an experimental wavefield acquired using a seismic source located at the centre of a linear array of 24 closely spaced surface sensors directly into a robust starting model for FWI. The CNN approach was able to produce 2-D starting models with seismic image misfits that were significantly less than the misfits from other common starting model approaches, and in many cases even less than the misfits obtained by FWI with inferior starting models. The ability of the CNN to generalize outside its two-layered training set was assessed using a more complex, three-layered, soil-over-bedrock formation. While the predictive ability of the CNN was slightly reduced for this more complex case, it was still able to achieve seismic image and waveform misfits that were comparable to other commonly used starting models, despite not being trained on any three-layered models. As such, CNNs show great potential as tools for rapidly developing robust, site-specific starting models for near-surface elastic FWI.

    

   more » « less