Block-diagonal Approximate Hessian for Preconditioning in Full Waveform Inversion

Author(s):  
N. Korta ◽  
A. Fichtner ◽  
V. Sallarčs
Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. R247-R259 ◽  
Author(s):  
Yuwei Wang ◽  
Liangguo Dong ◽  
Yuzhu Liu ◽  
Jizhong Yang

Elastic full-waveform inversion (EFWI) of multicomponent seismic data is a powerful tool for estimating the subsurface elastic parameters with high accuracy. However, the trade-offs between multiple parameters increase the nonlinearity of EFWI. Although the conventional diagonal-approximate Hessian matrix describes the illumination and limited bandwidth effects, it ignores the trade-off effects and decreases the convergence rate of EFWI. We have developed a block-diagonal pseudo-Hessian operator for 2D frequency-domain EFWI to take into account the approximate trade-offs among the P-wave (compressional-wave) velocity, S-wave (shear-wave) velocity, and density without extra computational costs on forward simulations. The Hessian matrix tends toward a block-diagonal matrix as the frequency grows to infinity; thus, the proposed block-diagonal pseudo-Hessian matrix is more accurate at higher frequencies. The inverse of the block-diagonal pseudo-Hessian matrix is used as a preconditioner for the nonlinear conjugate-gradient method to simultaneously reconstruct P- and S-wave velocities and density. This approach effectively mitigates the crosstalk artifacts by correcting the gradients from the trade-off effects and produces more rapid inversion convergence, which becomes more significant at higher frequencies. Synthetic experiments on an inclusion model and the elastic Marmousi2 model demonstrate its feasibility and validity in EFWI.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. R325-R337 ◽  
Author(s):  
Yuzhu Liu ◽  
Zheng Wu ◽  
Hao Kang ◽  
Jizhong Yang

The truncated Newton method uses information contained in the exact Hessian in full-waveform inversion (FWI). The exact Hessian physically contains information regarding doubly scattered waves, especially prismatic events. These waves are mainly caused by the scattering at steeply dipping structures, such as salt flanks and vertical or nearly vertical faults. We have systematically investigated the properties and applications of the exact Hessian. We begin by giving the formulas for computing each term in the exact Hessian and numerically analyzing their characteristics. We show that the second term in the exact Hessian may be comparable in magnitude to the first term. In particular, when there are apparent doubly scattered waves in the observed data, the influence of the second term may be dominant in the exact Hessian and the second term cannot be neglected. Next, we adopt a migration/demigration approach to compute the Gauss-Newton-descent direction and the Newton-descent direction using the approximate Hessian and the exact Hessian, respectively. In addition, we determine from the forward and the inverse perspectives that the second term in the exact Hessian not only contributes to the use of doubly scattered waves, but it also compensates for the use of single-scattering waves in FWI. Finally, we use three numerical examples to prove that by considering the second term in the exact Hessian, the role of prismatic waves in the observed data can be effectively revealed and steeply dipping structures can be reconstructed with higher accuracy.


2015 ◽  
Vol 6 (2) ◽  
pp. 5-16 ◽  
Author(s):  
Sergio Alberto Abreo Carrillo ◽  
Ana B. Ramirez ◽  
Oscar Reyes ◽  
David Leonardo Abreo-Carrillo ◽  
Herling González Alvarez

Sign in / Sign up

Export Citation Format

Share Document