Preconditioned Conjugate Gradient Algorithm-Based 2D Waveform Inversion for Shallow-Surface Site Characterization

Abstract

The full-waveform inversion (FWI) of a Love wave has become a powerful tool for shallow-surface site characterization. In classic conjugate gradient algorithm- (CG) based FWI, the energy distribution of the gradient calculated with the adjoint state method does not scale with increasing depth, resulting in diminished illumination capability and insufficient model updating. The inverse Hessian matrix (HM) can be used as a preprocessing operator to balance, filter, and regularize the gradient to strengthen the model illumination capabilities at depth and improve the inversion accuracy. However, the explicit calculation of the HM is unacceptable due to its large dimension in FWI. In this paper, we present a new method for obtaining the inverse HM of the Love wave FWI by referring to HM determination in inverse scattering theory to achieve a preconditioned gradient, and the preconditioned CG (PCG) is developed. This method uses the Love wave wavefield stress components to construct a pseudo-HM to avoid the huge calculation cost. It can effectively alleviate the influence of nonuniform coverage from source to receiver, including double scattering, transmission, and geometric diffusion, thus improving the inversion result. The superiority of the proposed algorithm is verified with two synthetic tests. The inversion results indicate that the PCG significantly improves the imaging accuracy of deep media, accelerates the convergence rate, and has strong antinoise ability, which can be attributed to the use of the pseudo-HM.