• 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.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.6-17
• /
• 2013

Multicomponent seismic data including both P- and PS-waves have advantages in discriminating the type of pore fluid, characterizing the lithologic attributes and producing the high resolution image. However, multicomponent seismic data recorded at the vertical and horizontal component receivers contain both P- and PS-waves which have different features, simultaneously. Therefore, the wavefield separation of P- and PS-waves as a preprocessing is inevitable in order to use the multicomponent seismic data successfully. In this study, we analyzed the previous study of the wavefield separation method suggested by Jeong and Byun in 2011, where the approximated reflection angle calculated only from one refernce depth is used in rotation transform, and showed its limitation for seismic data containing various reflected events from the multi-layered structure. In order to overcome its limitation, we suggested a new effective wavefield separation method of P- and PS-waves. In new method, we calculate the reflection angles with various reference depths and apply rotation transforms to the data with those reflection angles. Then we stack all results to obtain the final separated data. To verify our new method, we applied it to the synthetic data sets from a multi-layered model, a fault model, and the Marmousi-2 model. The results showed that the proposed method separated successfully P- and PS-reflection events from the multicomponent data from mild dipping layered model as long as the dip is not too steep.

• Geophysics and Geophysical Exploration
• /
• v.21 no.2
• /
• pp.112-124
• /
• 2018

Ocean-bottom seismic survey is a seismic acquisition technique which measures data by installing 4-component receiver on the sea floor. It can produce more improved data in quality than any other acquisition techniques. In the ocean-bottom seismic survey, however, the number of receivers is limited due to high cost. Since only a small number of receivers are used for acquisition, ocean-bottom seismic data may suffer from discontinuities of events over traces, which can result in spatial aliasing. In this paper, we implemented Kirchhoff migration using mirror-imaging algorithm to improve the quality of ocean-bottom seismic image. In order to implement the mirror imaging algorithm, the seismograms should be separated into up-going and down-going wavefields and the down-going wavefield should be used for migration. In this paper, we use the P-Z summation method to separate the wavefield. Numerical examples show that the migration results using mirror imaging algorithm have wider illumination than the conventional migration, especially in the shallow layers.

• The Journal of the Acoustical Society of Korea
• /
• v.15 no.4E
• /
• pp.9-13
• /
• 1996

We developed an iterative algorithm which could improve the resolution of reconstructed tomograms having random attenuation patterns and analyzed the limitation of this algorithm. The simple back-and forth propagation algorithm has depth resolution about four wavelengths. An iterative algorithm, based on back-and-forth propagation, can be used to improve the resolution of reconstructed tomograms. We analyzed the wavefield for multi-layered specimen and programmed iterative algorithm using Clanguage. Simulation results show that the images get clearer as the number of iterations increases. Also, unambiguous images can be reconstructed using this algorithm even when the layer separation is only two wavelengths. However, this iteration algorithm comes up with an incorrect solution for the number of projections less than five.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.234-241
• /
• 2011

Recently, exploration and development of oil sands are thriving due to high oil price. Because oil sands reservoir usually exists as a thin layer, multicomponent VSP, which has the advantage of the high-resolution around the borehole, is more effective than surface seismic survey in exploring oil sand reservoir. In addition, prestack phase-screen migration is effective for multicomponent seismic data because it is based on an one-way wave equation. In this study, we examined the applicability of the prestack phase-screen migration for multicomponent RVSP data to image the thin oil sand reservoir. As a preprocessing tool, we presented a method for separating P-wave and PS-wave from multicomponent RVSP data by using incidence angle and rotation matrix. To verify it, we have applied the developed wavefield separation method to synthetic data obtained from the velocity model including a horizontal layer and dipping layers. Also, we compared the migrated image by using P-wave with that by using PS-wave. As a result, the PS-wave migrated image has higher resolution and wide coverage than P-wave migrated image. Finally, we have applied the prestack phase-screen migration to the synthetic data from the velocity model simulating oil sand reservoir in Canada. The results show that the PS-wave migrated image describe the top and bottom boundaries of the thin oil sand reservoir more clearly than the P-wave migrated image.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.255-262
• /
• 2009

Conventionally, vertical-seismic-profiling (VSP) survey that provides high-resolution information has mainly performed to obtain the exact depth of the gas hydrate-bearing sediment, which is one of the key factors in the development of the gas hydrate. In this study, we extracted interval velocities and created corridor stacks from the first domestic zero-offset VSP data, which were acquired with three component receivers at UBGH09 borehole in Ulleung Basin where gas hydrate exists. Then we compared the corridor stacks with a CMP stacked section from surface seismic data. First of all, we converted the signals recorded with three component receivers to true vertical and horizontal components by phase rotation, and divided the data into direct waves and reflected waves by wavefield separation processing. The trend of the interval velocity extracted from the zero-offset VSP was similar to that of the sonic log obtained at the same borehole. Because the interval velocity of the gas hydrate-bearing sediment above the BSR was high, and it decreased suddenly through the BSR, we could infer that free gas is accumulated below the BSR. The results of comparing the corridor stacks to the CMP stacked section of the surface seismic data showed that most reflection events agreed well with those in the surface CMP stacked section and that the phase-rotated VSP data corresponded better with the surface seismic data than the VSP data without phase rotation. In addition, by comparing a corridor stack produced from the transverse component with the CMP stacked section of the surface seismic data, we could identify PS mode-converted reflections in the CMP stacked section.