• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Journal of the Korean Geophysical Society
• /
• v.3 no.1
• /
• pp.37-48
• /
• 2000

In order to investigate in-line ground motions caused by earthquakes, we examine the multicomponent complex trace analysis method (MCTAM) for the synthetic data and apply it to real earthquake data. An experimental result for synthetic data gives correct information on the arrival times, duration of individual phases, and approaching angles for body waves. Rayleigh waves are also easily identified with the MCTAM. A deep earthquake with magnitude of 7.3 was chosen to test various polarization attributes of ground motions. For P waves, instantaneous phase difference between the vertical and the in-line horizontal components ${\phi}(t)$, instantaneous reciprocal ellipticity ${\rho}(t)$, and approaching angle ${\tau}(t)$ are computed to be ${\pm}180^{\circ},\;0{\sim}0.25,\;and\;-30^{\circ}{\sim}-45^{\circ}$, respectively. For S waves, ${\phi}(t)$ tends to vary while ${\rho}(t)$ have values of $0{\sim}0.3\;and\;{\tau}(t)$ remains near vertical, respectively. A relatively low frequency signal registered just prior to the S wave event is interpreted as a P-wave phase based on its polarization characteristics. Velocities of P and S waves are computed to be 8.633 km/s and 4.762 km/s, and their raypath parameters 0.074 s/km and 0.197 s/km. Dynamic Poisson's ratio is obtained as 0.281 from the velocities of P and S waves.