• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.111-123
• /
• 2007

Unlike the conventional reverse time migration method which is implemented by simply extrapolating wavefield in reverse time, this paper presents a derivation of another reverse time migration operator as the exact adjoint of the presumed forward wavefield extrapolation operator. The adjoint operator is obtained by formulating the forward time extrapolation operator in an explicit matrix equation form and then taking the adjoint to this matrix equation followed by determining the corresponding operator. The reverse time migration operator as the exact adjoint to the implied forward operator can be used not only as a migration algorithm but also as an adjoint operator which is required in the imaging through an inversion such as least-squares migration.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.70-75
• /
• 2000

In the seismic migration, Kirchhoff and reverse time migration are used in general. In the reverse time migration using wave equation, two-way and one-way wave equation are applied. The approach of one-way wave equation uses approximately computed downward continuation extrapolator, it need tess amounts of calculations and core memory in compared to that of two-way wave equation. In this paper, we applied one-way wave equation to pre-stack reverse time migration. In the frequency-space domain, forward propagation of source wavefield and back propagration of measured wavefield were executed by using monochromatic one-way wave equation, and zero-lag cross correlation of two wavefield resulted in the image of subsurface. We had implemented prestack migration on a massively parallel processors (MPP) CRAYT3E, and knew the algorithm studied here is efficiently applied to the prestck migration due to its suitability for parallelization.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.315-322
• /
• 2010

The phase-screen and the split-step Fourier migrations, which are implemented in both the frequency-wavenumber and frequency-space domains by using one-way scalar wave equation, allow imaging in laterally heterogeneous media with less computing time and efficiency. The generalized-screen migration employs the series expansion of the exponential, unlike the phase-screen and the split-step Fourier migrations which assume the vertical propagation in frequency-wavenumber domain. In addition, since the generalized-screen migration generalizes the series expansion of the vertical slowness, it can utilize higher-order terms of that series expansion. As a result, the generalized-screen migration has higher accuracy in computing the propagation with wide angles than the phase-screen and split-step Fourier migrations for media with large and rapid lateral velocity variations. In this study, we developed a 2D prestack generalized-screen migration module for imaging a complex subsurface efficiently, which includes various dips and large lateral variations. We compared the generalized-screen propagator with the phase-screen propagator for a constant perturbation model and the SEG/EAGE salt dome model. The generalized-screen propagator was more accurate than the phase-screen propagator in computing the propagation with wide angles. Furthermore, the more the higher-order terms were added for the generalized-screen propagator, the more the accuracy was increased. Finally, we compared the results of the generalizedscreen migration with those of the phase-screen migration for a model which included various dips and large lateral velocity variations and the synthetic data of the SEG/EAGE salt dome model. In the generalized-screen migration section, reflectors were positioned more accurately than in the phase-screen migration section.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.191-202
• /
• 2011

Reserve-time migration (RTM) using a two-way wave equation is one of the most accurate migration techniques. RTM has been conducted by assuming that subsurface media are isotropic. However, anisotropic media are commonly encountered in reality. Conventional isotropic RTM may yield inaccurate results for anisotropic media. In this paper, we develop RTM algorithms for vertical transversely isotropic media (VTI) and tilted transversely isotropic media (TTI). For this, the pseudo-acoustic wave equations are used. The modeling algorithms are based on the high-order finite-difference method (FDM). The RTM algorithms are composed using the cross-correlation imaging condition or the imaging condition using virtual sources. By applying the developed RTM algorithms to the Hess VTI and BP TTI models, we could obtain better images than those obtained by the conventional isotropic RTM.

• Geophysics and Geophysical Exploration
• /
• v.6 no.3
• /
• pp.153-159
• /
• 2003

Using reciprocity theorem, one direction aperture of seismogram can be extended to full aperture seismogram. Modified seismogram is applied to reverse time migration only to acquire improved migration image. In this paper, we tested reverse time migration with the Marmousi velocity data to examine efficiency of modified seismogram. And linearly increasing velocity model is selected and examined for a case where velocity data is insufficient. When true velocity is applied, using modified seismogram enhances the reverse time migration image more than using original seismogram. In the case of using linearly increasing velocity model, migration image is distorted. So low frequency source is brought in migration process. Reverse time migration image with low frequency source and linearly increasing velocity model is improved when modified seismogram is used. From the result of study, seismogram modification by reciprocity theorem is useful and migration image can be enhanced.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.150-152
• /
• 2002

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.3
• /
• pp.247-258
• /
• 2004

Two 3-dimensional data processing techniques to predict the fractured zone ahead of a tunnel face by the tunnel seismic survey were proposed so that the geometric formation of the fractured zone could be estimated. The first 3-dimensional data processing technique was developed based on the principle of ellipsoid, The input data needed for the 3D migration can be obtained from the 2-dimensional tunnel seismic prediction (TSP) test where the TSP test should be performed in each sidewall of a tunnel. The second 3-dimensional migration technique that was developed based on the concept of wave travel plane was proposed. This technique can be applied when the TSP is operated with sources in one sidewall of a tunnel while the receivers are installed in both sidewalls. New migration technique was applied to an in-situ tunnelling site. The 3-dimensional migration was performed using measured TSP data and its results were compared with the geological investigation results that were monitored during tunnel construction. This comparison revealed that the proposed migration technique could reconstruct the discontinuity planes reasonably well.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2018.05a
• /
• pp.96-96
• /
• 2018

최근 기후변동성으로 유발되는 불안정한 기상상태를 효과적으로 관측하고자 기상레이더가 도입되고 있다. 기상레이더는 경험식으로 산정된 Z-R 관계식을 통하여 레이더 강수량을 제시하게 된다. 이 과정에서 레이더 강수량은 필연적으로 실제 지상에 도달하는 강수량과는 정량적으로 오차가 발생하게 된다. 레이더 강수량에 포함된 오차는 다양한 원인으로 발생하게 되므로 레이더 강수량의 오차 성분을 규명하는 것은 레이더 강수량 활용을 위하여 필수적으로 선행되어야 한다. 본 연구는 지상강수량과 레이더 강수량의 편의를 보정하기 위한 확률통계학적 방법론을 개발하였다. 레이더 강수량의 편의오차를 보정하기 위하여 수문통계학에서 널리 활용되고 있는 계층적 Bayesian 구조를 기반으로 하였으며 자료통합(data pooling) 기법을 이용하여 편의보정 매개변수 추정과정의 불확실성 추정 효율성을 증대시켰다. 본 연구를 통하여 개발된 레이더 강수량 편의보정기법은 계층적 Bayesian 구조를 도입함으로써 편의보정 매개계수의 불확실성을 정량적으로 제시하였으며 유역 단위의 강수상관성을 현실적으로 복원하는 것을 확인하였다. 따라서 본 연구에서 제안하는 편의보정기법은 편의보정 과정에서 발생할 수 있는 매개변수의 불확실성 및 레이더 강수량의 오차구조를 정량적으로 규명하여 고해상도의 강수정보를 생산함으로써 고도화된 수문해석을 가능케 할 것으로 판단된다.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.102-107
• /
• 2006

Recently, rapid developments in computer hardware have enabled reverse-time migration to be applied to various production imaging problems. As a wave-equation technique using the two-way wave equation, reverse-time migration can handle not only multi-path arrivals but also steep dips and overturned reflections. However, reverse-time migration causes unwanted artefacts, which arise from the two-way characteristics of the hyperbolic wave equation. Zero-lag cross correlation with diving waves, head waves and back-scattered waves result in spurious artefacts. These strong artefacts have the common feature that the correlating forward and backward wavefields propagate in almost the opposite direction to each other at each correlation point. This is because the ray paths of the forward and backward wavefields are almost identical. In this paper, we present several tactics to avoid artefacts in shot-domain reverse-time migration. Simple muting of a shot gather before migration, or wavefront migration which performs correlation only within a time window following first arriving travel times, are useful in suppressing artefacts. Calculating the wave propagation direction from the Poynting vector gives rise to a new imaging condition, which can eliminate strong artefacts and can produce common image gathers in the reflection angle domain.

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