Waveform inversion with structural regularizing constraint based on gradient decomposition

Abstract

Abstract Full waveform inversion (FWI) can simultaneously update low-to-medium wavenumber velocity components and high-wavenumber velocity components. However, if seismic data lack large-offset data and effective low-frequency components, FWI updates will be dominated by high-wavenumber velocity perturbation. Meanwhile, providing that the initial model is inaccurate, inversion will have the problem of local minima. In this study, FWI is developed with structural regularizing constraint based on gradient decomposition (RGDFWI). By correlating the separated forward wavefield and backward wavefield with specific propagating direction, FWI gradient is decomposed into tomography-mode gradient and migration-mode gradient. We propose an optimized strategy taking full advantage of the two modes of FWI gradient. On the one hand, we use tomography-mode gradient to enhance low-to-medium wavenumber updates. On the other hand, we use migration-mode gradient to apply structural regularizing constraint by estimating structure dip and adding sparsity constraint in Seislet domain. During the inversion process, high-wavenumber structural information constrains and guides low-wavenumber model updates. The results of two numerical tests, Marmousi model test and Overthrust model test, validate the optimized strategy, which can produce a better initial velocity model for FWI. The inversion finally generates a high-precision and high-resolution velocity model.