Abstract:Most finite-frequency traveltime tomography methods are based on the Born approximation, which requires that the scale of the velocity heterogeneity and the magnitude of the velocity perturbation should be small enough to satisfy the Born approximation. On the contrary, the Rytov approximation works well for large-scale velocity heterogeneity. Typically, the Rytov-approximation-based finite-frequency traveltime sensitivity kernel (Rytov-FFTSK) can be obtained by integrating the phase-delay sensitivity kernels with a normalized weighting function, in which the calculation of sensitivity kernels requires the numerical solution of Green's function. However, solving the Green's function explicitly is quite computationally demanding, especially for 3D problems. To avoid explicit calculation of the Green's function, we show that the Rytov-FFTSK can be obtained by crosscorrelating a forward-propagated incident wavefield and reverse-propagated adjoint wavefield in the time domain. In addition, we find that the action of the Rytov-FFTSK on a model-space vector, e.g., the product of the sensitivity kernel and a vector, can be computed by calculating the inner product of two time-domain fields. Consequently, the Hessian-vector product can be computed in a matrix-free fashion (i.e., first calculate the product of the sensitivity kernel and a model-space vector and then calculate the product of the transposed sensitivity kernel and a data-space vector), without forming the Hessian matrix or the sensitivity kernels explicitly. We solve the traveltime inverse problem with the Gauss-Newton method, in which the Gauss-Newton equation is approximately solved by the conjugate gradient using our matrix-free Hessian-vector product method. An example with a perfect acquisition geometry found that our Rytov-approximation-based traveltime inversion method can produce a high-quality inversion result with a very fast convergence rate. An overthrust synthetic data test demonstrates that large- to intermediate-scale model perturbations can be recovered by diving waves if long-offset acquisition is available.

Crosscorrelation-based Dynamic Time Warping and Its Application in Wave Equation Reflection Traveltime Inversion

Dynamic Time Warping As an Alternative to Windowed Cross Correlation in Seismological Applications

Multilevel Dynamic Time Warping: A Parameter-Light Method for Fast Time Series Classification

Improvement and Application of Hale's Dynamic Time Warping Algorithm

Time Reversal Based Channel Tracking for Underwater Acoustic Communications

Elastic Wave Equation Traveltime Inversion With Dynamic Time Warping Based on High-Speed Train Seismic Data

Suppression of Interference by Direct Arrival Wave Using Dynamic Time Warping Method

Time-Lapse Full-Waveform Inversion Using Cross-Correlation-Based Dynamic Time Warping

ESDTW: Extrema-based shape dynamic time warping

Wave equation traveltime tomography using Rytov approximation

Comparison between Time Shifting Deviation and Cross-correlation Methods

Rytov-approximation-based Wave-Equation Traveltime Tomography

Using Multi-Dimensional Dynamic Time Warping to Identify Time-Varying Lead-Lag Relationships

Affine and Regional Dynamic Time Warpng

Wave-equation Traveltime Slope Inversion by Combining Finite Difference and Cross-Correlation Methods

Asymptotic Dynamic Time Warping Calculation with Utilizing Value Repetition

Exact Dynamic Time Warping calculation for weak sparse time series

Data-domain wave equation reflection traveltime tomography

Seismic data extrapolation based on multi-scale dynamic time warping

A DTW Distance-Based Seismic Waveform Clustering Method for Layers of Varying Thickness

Robust Full-Waveform Inversion Based on Automatic Differentiation and Differentiable Dynamic Time Warping