skip to main content

Title: Decoupled Fréchet kernels based on a fractional viscoacoustic wave equation
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
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Page Range / eLocation ID:
T61 to T70
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    To complement velocity distributions, seismic attenuation provides additional important information on fluid properties of hydrocarbon reservoirs in exploration seismology, as well as temperature distributions, partial melting, and water content within the crust and mantle in earthquake seismology. Full waveform inversion (FWI), as one of the state‐of‐the‐art seismic imaging techniques, can produce high‐resolution constraints for subsurface (an)elastic parameters by minimizing the difference between observed and predicted seismograms. Traditional waveform inversion for attenuation is commonly based on the standard‐linear‐solid (SLS) wave equation, in which case the quality factor (Q) has to be converted to stress and strain relaxation times. When using multiple attenuation mechanisms in the SLS method, it is difficult to directly estimate these relaxation time parameters. Based on a time domain complex‐valued viscoacoustic wave equation, we present an FWI framework for simultaneously estimating subsurfacePwave velocity and attenuation distributions. BecauseQis explicitly incorporated into the viscoacoustic wave equation, we directly derivePwave velocity andQsensitivity kernels using the adjoint‐state method and simultaneously estimate their subsurface distributions. By analyzing the Gauss‐Newton Hessian, we observe strong interparameter crosstalk, especially the leakage from velocity toQ. We approximate the Hessian inverse using a preconditioned L‐BFGS method in viscoacoustic FWI, which enables us to successfully reduce interparameter crosstalk and produce accurate velocity and attenuation models. Numerical examples demonstrate the feasibility and robustness of the proposed method for simultaneously mapping complex velocity andQdistributions in the subsurface.

    more » « less

    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. Abstract

    We present the first continental‐scale seismic model of the lithosphere and underlying mantle beneath Southeast Asia obtained from adjoint waveform tomography (often referred to as full‐waveform inversion or FWI), using seismic data filtered at periods from 20 to 150 s. Based on >3,000 hr of analyzed waveform data gathered from ∼13,000 unique source‐receiver pairs, we image isotropicP‐wave velocity, radially anisotropicS‐wave velocity and density via an iterative non‐linear inversion that begins from a 1‐D reference model. At each iteration, the full 3‐D wavefield is determined through an anelastic Earth, accommodating effects of topography, bathymetry and ocean load. Our data selection aims to maximize sensitivity to deep structure by accounting for body wave arrivals separately.SASSY21, our final model after 87 iterations across seven period bands, is able to explain true‐amplitude data from events and receivers not included in the inversion. The trade‐off between inversion parameters is estimated through an analysis of the Hessian‐vector product.SASSY21reveals detailed anomalies down to the mantle transition zone, including multiple subduction zones. The most prominent feature is the (Indo‐)Australian plate descending beneath Indonesia, which is imaged as one continuous slab along the 180° curvature of the Banda Arc. The tomography confirms the existence of a hole in the slab beneath Mount Tambora and locates a highS‐wave velocity zone beneath northern Borneo that may be associated with subduction termination in the mid‐late Miocene. A previously undiscovered feature beneath the east coast of Borneo is also revealed, which may be a signature of post‐subduction processes, delamination or underthrusting from the formation of Sulawesi.

    more » « less

    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

    A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth’s viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

    more » « less