Numerical simulations of pure quasi-P-waves in orthorhombic anisotropic media

Abstract

The accurate simulation of anisotropic media is critical in seismic imaging and inversion. In recent years, some scholars have dedicated efforts to the study of precise elastic waves in anisotropic media; however, it is easy to separate P-wave and S-wave from elastic wave fields in isotropic media but difficult to separate them in anisotropic media. To address this issue, others have proposed pseudo-pure-wave equations based on the theory of wave-mode separation, but shear wave interference still exists. Therefore, we derived the first-order pure quasi-P-wave equation with no shear wave component in orthorhombic anisotropic media (ORT) which is common in the Earth’s crust and has very important research value. The presence of a pseudodifferential operator in the equation poses a challenge for solving. In order to solve the pure wave equation, we decomposed the original pseudodifferential operator into an elliptic differential operator and a scalar operator, both of which are easily solvable. In addition, we extended the equation from ORT media to tilted ORT (TORT) media. The example results indicate that our pure quasi-P-wave equation can yield a more stable and accurate P-wave field. The pure wave equation we propose can be applied in reverse time migration (RTM), the least squares RTM (LSRTM), and even the full waveform inversion (FWI).