Joint FWI of Active Source Data and Passive Virtual Source Data Reconstructed Using an Improved Multidimensional Deconvolution

Traditional full waveform inversion (FWI) highly depends on sufficient low-frequency data or a good initial model. Passive seismic data contain rich low-frequency components, and passive seismic FWI using virtual source data by seismic interferometry (SI) is a promising method. However, the distribution of passive sources in the subsurface is always inhomogeneous, which will lead to artifacts in the reconstruction results by SI using cross-correlation (CC). SI by multidimensional deconvolution (MDD) can counteract the source inhomogeneity but requires the separation of the reference wavefields, which is difficult to achieve for noise source data. To mitigate this problem, we propose an improved SI method by linear Radon transform based multidimensional deconvolution (LRTMDD). The interferometric point-spread function can be extracted more accurately and efficiently in the linear Radon domain, thus improving the reconstruction results. The passive virtual source full waveform inversion (PVSFWI) based on LRTMDD is further proposed, which can effectively use the low-frequency information in the virtual source data to invert the macroscopic velocity structures even in the case of inhomogeneous source distributions, and without the need to estimate the virtual source wavelets. A joint multi-source FWI strategy is proposed to solve the problem of missing low-frequency data suffered by active source FWI. Numerical experiments on the Marmousi model and the SEG/EAGE overthrust model show that the proposed methods can fully combine the respective advantages of multisource seismic data to stably achieve high-accuracy velocity models in the case of inhomogeneous passive source distributions and the lack of low-frequency data in active seismic data.