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.

Fast Multi-Trace Impedance Inversion Using Anisotropic Total P-Variation Regularization in the Frequency Domain

Multichannel impedance inversion in the frequency domain via anisotropic total variation with overlapping group sparsity regularization

Seismic Acoustic Impedance Inversion Based on Arctangent Total Variation and Hybrid Domain Constraints

Seismic Acoustic Impedance Inversion Using Mixed Second-Order Fractional ATpV Regularization.

Multi-trace Seismic Impedance Inversion with Structure-oriented Minimum Entropy Stabilizer

Data-driven Multichannel Seismic Impedance Inversion with Anisotropic Total Variation Regularization

Seismic Acoustic Impedance Inversion with Multi-Parameter Regularization

Relative acoustic impedance inversion via Toeplitz-Sparse Matrix Factorization

Multichannel Blind Acoustic Impedance Inversion Based on 2D-TV Regularization

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Structure-Guided Multiscale Impedance Inversion Based on Modified Total Variation Regularization

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

Sharp and Laterally Constrained Multitrace Impedance Inversion Based on Blocky Coordinate Descent

Impedance Inversion Based on L1 Norm Regularization

An Improved TV-Type Variational Regularization Method for Seismic Impedance Inversion

Seismic Acoustic Impedance Inversion Using Reweighted L1-Norm Sparse Constraint

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

Seismic Absolute Acoustic Impedance Inversion with L1 Norm Reflectivity Constraint and Combined First- and Second-Order Total Variation Regularizations

A New Approach for Blind Nonlinear Acoustic Impedance Inversion

Acoustic- and elastic-waveform inversion with total generalized p-variation regularization

Post‐stack impedance blocky inversion based on analytic solution of viscous acoustic wave equation