Scattering Gaussian wave packet modeling and imaging method in anisotropic media based on Gabor function

Wave anisotropy exists widely in the earth's medium, which provides new information for the study of the earth' s internal structures and components. Nowadays, the anisotropy of media plays a more and more important role in seismic wave propagation, imaging and inversion. Gaussian wave packet/beam forward and inversion methods have made great progress in seismic background velocity modeling and imaging in recent years because of their high efficiency and simplicity. Therefore, under the pseudo acoustic approximation, this paper develops a scattering Gaussian wave packet modeling and imaging method for anisotropic media based on Gabor framework. Different from the traditional grid parameter perturbation method, the characteristics of spatially distributed localized Gabor basis functions are used to express the parameter perturbation of underground media. Then the mapping relationship between Gabor function in model domain and Gaussian wave packet in data domain is established. This method not only inherits the advantages and characteristics of isotropic Gaussian wave packet, such as (1) the scattering Gaussian wave packet modeling based on Gabor framework has the characteristics of high computational efficiency and small memory occupation; (2) it can deal with multi wave arrivals and can be applied to the imaging of complex structures (salt dome, gas cloud boundary) and middle or deep structures; (3) Gabor function can be used to sparsely express the underground structures, which is suitable for target oriented flexible imaging; in addition, (4) compared with the isotropic Gaussian wave packet imaging method, considering the anisotropy of the medium, the results in this paper are more accurate and the energy of the events of the imaging section is more focused. The proposed method is tested by numerical experiments. The experimental results confirm the effectiveness, efficiency and flexibility of the method.