Abstract:It is important to reconstruct both seismic velocities and anisotropy parameters by anisotropic full waveform inversion (FWI). But the trade-off between different parameter classes increases its nonlinearity severely. The significant difference of magnitude orders between seismic velocities and anisotropy parameters makes the inversion problem poorly conditioned. Judiciously incorporating the inverse Hessian operator within the FWI scheme can help to alleviate the negative impacts of these two problems on FWI. The truncated Newton (TN) method can better estimate the inverse Hessian operator. Therefore, we use the TN method to tackle these two problems and simultaneously reconstruct both seismic velocities and anisotropy parameters for acoustic VTI media.The TN method has two advantages. On one hand, it does not need explicit Hessian matrix, which is prohibitive to build due to high computational cost and large memory requirements. To get the updated direction, the TN method solves Newton equation at each time of iteration with the conjugate gradient method (CG), which only uses Hessian-vector products. We compute these products by the finite difference method. On the other hand, the TN method can adjust flexibly the updated direction when the Hessian matrix is non-positive definite. The CG method at each step of the TN iteration is terminated as soon as a direction of negative curvature is generated. With this strategy the TN method may perform better than Limited-memory Broyden-Fletcher-Goldfarb-Shanno method (L-BFGS) for ill-conditioned problems.We test the inversion capability of the TN method with two numerical experiments. In each experiment, the NMO velocity and anisotropy parameter eta are inverted simultaneously. The L-BFGS method is also tested for comparison. The first experiment is on a simple 2D VTI homogenous model with three different anomalies inside. The inversion is performed in the time domain with a full acquisition system. Inversion results show that the TN method can balance the updates of the two different parameters, and reconstruct more accurate models. In contrast, the L-BFGS method overestimates parameter eta but underestimates NMO velocity. The second experiment is on a more complex 2D model with a surface acquisition system. The two methods perform the same as they do in the first experiment.The reason for their different performances may be found from the Hessian analysis part that we have finished before the two experiments. We design a small size (51 X 51) inclusion model and compute explicitly the Hessian matrix of the background initial model. The analyses show that the Hessian matrix has several negative eigenvalues and a large condition number (ill-conditioned matrix), which may result in the failure of the L-BFGS method. The TN method better overcomes these difficulties. Hence, in the context of simultaneous inversion of two parameters for acoustic VTI media, the TN method can account for more accurate inverse Hessian operator than the L-BFGS method and get more accurate inversion results. It is a promising method for multi-parameter FWI.

A New Method for Dynamic Inversion of Anisotropic Equations

Full Waveform Inversion Based on Time-Integral-damping Wavefield

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

A New Inversion Method of Sedimentary Strata for Deriving Double Parameters and Its Applications

On the method of multidimensional Q-factor and velocity inversion with VSP data

A Joint Inversion of Sources and Seismic Waveforms for Velocity Distribution: 1-D and 2-D Examples

An efficient elastic full‐waveform inversion of multiple parameters with ocean‐bottom seismometer data

Monochromatic wave-equation traveltime inversion in vertical transverse isotropic media: Theory and application

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

Multi-Parameter Full Waveform Inversion for Acoustic VTI Media Using the Truncated Newton Method

2-D acoustic wave equation inversion combining plane-wave seismograms with well logs

JOINT INVERSION OF RAYLEIGH WAVE VERTICAL‐HORIZONTAL AMPLITUDE RATIOS AND DISPERSION BASED ON THE NEIGHBORHOOD ALGORITHM AND ITS APPLICATION

Efficient scattering approach to seismic full-waveform inversion in anisotropic elastic media with variable density

Optimal Inversion Method Based on Joint Waveform Inversion and Least Squares Reverse Time Migration

2D Frequency-Domain Elastic Full-Waveform Inversion Using the Block-Diagonal Pseudo-Hessian Approximation

Joint full-waveform inversion of towed-streamer and ocean-bottom seismic data with amplitude weighting

Elastic multi-parameter full-waveform inversion for anisotropic media

Dynamic inversion method based on the time-staggered stereo-modeling scheme and its acceleration

First-arrival multi-information joint inversion method based on wave-equation

Inversion of Rayleigh-wave Ellipticity for Shallow Strong-Velocity-contrast Structures

A Multi-parameter Full Waveform Inversion Strategy in Acoustic Media