Angle-domain generalized Radon transform for elastic multiparameter inverse scattering inversion

Abstract

Linearized algorithms based on the Born approximation are well-known and popular techniques for quantitative seismic imaging and inversion. However, linearization methods usually suffer from some significant problems, such as the computational cost for the required number of iterations, requirement for background models, and uncertain and unstable multiparameter extraction, which make the methods difficult to implement in practical applications. To avoid these problems, we have developed an angle-domain generalized Radon transform (AD-GRT) inversion in 2D elastic isotropic media. This AD-GRT is an approximate transform between the seismic data and an angle-domain model, which acts as a scattering function, and the seismic data can be reconstructed accurately, even when the background models are incorrect. The density and Lamé moduli perturbation parameters can be extracted stably from the inverted angle-domain scattering function. Deconvolution of the source wavelet is taken into account to remove the effect of the wavelet and improve the resolution and accuracy of the inversion results. The derived AD-GRT inversion is noniterative and is as efficient as the traditional elastic GRT method. The additional dimension of the angle domain has little effect on the computational cost of the AD-GRT, as opposed to other extended-domain inversion/migration methods. Our method also can be used to solve nonlinear Born inversion problems using iteration, which can significantly improve their convergence rate. Three numerical examples illustrate that the angle-domain scattering function inversion, data reconstruction, and multiparameter extraction using the presented AD-GRT inversion are effective.