More Like this

1. Eulerian partial-differential-equation methods for complex-valued eikonals in attenuating media

   https://doi.org/10.1190/GEO2020-0659.1  

   Hu, Jiangtao; Qian, Jianliang; Song, Jian; Ouyang, Min; Cao, Junxing; Leung, Shingyu (July 2021, GEOPHYSICS)

   null (Ed.)

   Seismic waves in earth media usually undergo attenuation, causing energy losses and phase distortions. In the regime of high-frequency asymptotics, a complex-valued eikonal is an essential ingredient for describing wave propagation in attenuating media, where the real and imaginary parts of the eikonal function capture dispersion effects and amplitude attenuation of seismic waves, respectively. Conventionally, such a complex-valued eikonal is mainly computed either by tracing rays exactly in complex space or by tracing rays approximately in real space so that the resulting eikonal is distributed irregularly in real space. However, seismic data processing methods, such as prestack depth migration and tomography, usually require uniformly distributed complex-valued eikonals. Therefore, we have developed a unified framework to Eulerianize several popular approximate real-space ray-tracing methods for complex-valued eikonals so that the real and imaginary parts of the eikonal function satisfy the classic real-space eikonal equation and a novel real-space advection equation, respectively, and we dub the resulting method the Eulerian partial-differential-equation method. We further develop highly efficient high-order methods to solve these two equations by using the factorization idea and the Lax-Friedrichs weighted essentially nonoscillatory schemes. Numerical examples demonstrate that our method yields highly accurate complex-valued eikonals, analogous to those from ray-tracing methods. Our methods can be useful for migration and tomography in attenuating media. 

   more » « less
2. 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
3. Liouville partial-differential-equation methods for computing 2D complex multivalued eikonals in attenuating media

   https://doi.org/10.1190/GEO2021-0113.1  

   Leung, Shingyu; Hu, Jiangtao; Qian, Jianliang (March 2022, GEOPHYSICS)

   We have developed a Liouville partial-differential-equation (PDE)-based method for computing complex-valued eikonals in real phase space in the multivalued sense in attenuating media with frequency-independent qualify factors, where the new method computes the real and imaginary parts of the complex-valued eikonal in two steps by solving Liouville equations in real phase space. Because the earth is composed of attenuating materials, seismic waves usually attenuate so that seismic data processing calls for properly treating the resulting energy losses and phase distortions of wave propagation. In the regime of high-frequency asymptotics, the complex-valued eikonal is one essential ingredient for describing wave propagation in attenuating media because this unique quantity summarizes two wave properties into one function: Its real part describes the wave kinematics and its imaginary part captures the effects of phase dispersion and amplitude attenuation. Because some popular ordinary-differential-equation (ODE)-based ray-tracing methods for computing complex-valued eikonals in real space distribute the eikonal function irregularly in real space, we are motivated to develop PDE-based Eulerian methods for computing such complex-valued eikonals in real space on regular meshes. Therefore, we solved novel paraxial Liouville PDEs in real phase space so that we can compute the real and imaginary parts of the complex-valued eikonal in the multivalued sense on regular meshes. We call the resulting method the Liouville PDE method for complex-valued multivalued eikonals in attenuating media; moreover, this new method provides a unified framework for Eulerianizing several popular approximate real-space ray-tracing methods for complex-valued eikonals, such as viscoacoustic ray tracing, real viscoelastic ray tracing, and real elastic ray tracing. In addition, we also provide Liouville PDE formulations for computing multivalued ray amplitudes in a weakly viscoacoustic medium. Numerical examples, including a synthetic gas-cloud model, demonstrate that our methods yield highly accurate complex-valued eikonals in the multivalued sense. 

   more » « less
4. EikoNet: Solving the Eikonal Equation With Deep Neural Networks

   https://doi.org/10.1109/tgrs.2020.3039165  

   Smith, Jonathan D.; Azizzadenesheli, Kamyar; Ross, Zachary E. (December 2020, IEEE Transactions on Geoscience and Remote Sensing)

   null (Ed.)

   The recent deep learning revolution has created enormous opportunities for accelerating compute capabilities in the context of physics-based simulations. In this article, we propose EikoNet, a deep learning approach to solving the Eikonal equation, which characterizes the first-arrival-time field in heterogeneous 3-D velocity structures. Our grid-free approach allows for rapid determination of the travel time between any two points within a continuous 3-D domain. These travel time solutions are allowed to violate the differential equation--which casts the problem as one of optimization--with the goal of finding network parameters that minimize the degree to which the equation is violated. In doing so, the method exploits the differentiability of neural networks to calculate the spatial gradients analytically, meaning that the network can be trained on its own without ever needing solutions from a finite-difference algorithm. EikoNet is rigorously tested on several velocity models and sampling methods to demonstrate robustness and versatility. Training and inference are highly parallelized, making the approach well-suited for GPUs. EikoNet has low memory overhead and further avoids the need for travel-time lookup tables. The developed approach has important applications to earthquake hypocenter inversion, ray multipathing, and tomographic modeling, as well as to other fields beyond seismology where ray tracing is essential. 

   more » « less
5. Heterogeneous mantle effects on the behaviour of SmKS waves and outermost core imaging

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

   Frost, Daniel_A; Garnero, Edward_J; Creasy, Neala; Wolf, Jonathan; Bozdağ, Ebru; Long, Maureen_D; Aderoju, Adeolu; Vite, Reynaldo (April 2024, Geophysical Journal International)

   SUMMARY Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure. 

   more » « less