• Economic and Environmental Geology
• /
• v.43 no.4
• /
• pp.349-358
• /
• 2010

For imaging complex subsurface structures such as salt dome, faults, thrust belt, and folds, seismic prestack reverse-time migration in depth domain is widely used, which is performed by the cross-correlation of shot-domain wavefield extrapolation with receiver-domain wavefield extrapolation. We apply the prestack reverse-time migration, which had been developed at KIGAM, to the seismic field data set of Block 5 in Jeju basin of Korea continental shelf in order to improve subsurface syncline stratigraphy image of the deep structures under the shot point 8km at the surface. We performed basic data processing for improving S/N ratio in the shot gathers, and constructed a velocity model from stack velocity which was calculated by the iterative velocity spectrum. The syncline structure of the stack image appears as disconnected interfaces due to the diffractions, but the result of the prestack migration shows that the syncline image is improved as seismic energy is concentrated on the geological interfaces.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.11a
• /
• pp.564-568
• /
• 2007

Gas hydrate has been attractive topic for two dedicates because it may cause the global warming, ocean hazards associated with the instability of marine slope due to the gas hydrate release as well as high potential of future energy resources. The study on gas hydrate in Ulleung basin has been performed since 1999 to explore the potential and distribution of gas hydrate offshore Korea. The numerous multi channel seismic data have been acquired and processed by Korea Institute of Geosciences and Mineral Resources (KIGAM). The results showed clearly the gas hydrate indicators such as pull up structure, bottom simulating reflector (BSR), seismic blanking zone. The prestack depth migration has been considered as fast and accurate technique to image the subsurface. In this paper, we will present both the conventional seismic data processing and apply Kirchhoff prestack depth migration for gas hydrate data set. The results will be applied for core sample collections and for proposal more detail 2D with long offset or 3D seismic exploration.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.10a
• /
• pp.203-207
• /
• 2008

Prestack depth migration is used to image for complex geological structure such as faults, folds, and subsalt. In this case, it is widely used the surface reflection data as a input data. However, the surface reflection data have intrinsic problems to image the subsalt and the salt flank due to the complex wavefields and multiples which come from overburden. For overcoming the structural defect of the surface reflection data in the imaging, I used the virtual sources in terms of seismic interferometry to image the subsurface and suppress the multiples using the velocity model of the lower part of the virtual sources. The results of the prestack depth migration using virtual source gathers and velocity model below receivers are similar geological interfaces to the results from shot gathers of the conventional ocean bottom seismic survey. And especially artificial interfaces by multiples were suppressed without applying any other data processing to eliminate multiples. This study results by numerical modeling can make a valuable imaging tool when it is applied to satisfied field data for specific condition.

• 한국신재생에너지학회:학술대회논문집
• /
• 2007.06a
• /
• pp.493-496
• /
• 2007

Korean Institute of Geosciences and Mineral Resources (KIGAM) has studied on gas hydrate in the Ulleung Basin, East sea of Korea since 1997. Most of all, a evidence for existence of gas hydrate, possible new energy resources, in seismic reflection data is bottom simulating reflection (BSR) which parallel to the sea bottom. Here we conducted the conventional data processing for gas hydrate data and Kirchhoff prestack depth migration. Kirchhoff migration is widely used for pre- and post-stack migration might be helpful to better image as well as to get the geological information. The processed stack image by GEOBIT showed some geological structures such as faults and shallow gas hydrate seeping area indicated by strong BSR. The BSR in the stack image showed at TWT 3.07s between shot gather No 3940 to No 4120. The estimated gas seeping area occurred at the shot point No 4187 to No 4203 and it seems to have some minor faults at shot point No 3735, 3791, 3947 and 4120. According to the result of depth migration, the BSR showed as 2.3km below the sea bottom.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.54-62
• /
• 1999

This paper describes a datuming scheme for a prestack dataset which uses wavefield depth extrapolation. The formulation of the prestack datuming algorithm is performed by finding the adjoint operator to the corresponding forward prestack wavefield extrapolation from a flat surface to an irregular surface. Here I used double-square-root (DSR) equation to extrapolate wavefield in prestack sense. This correspond to the forward model of the well known `survey sinking` prestack imaging algorithm.

• 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.4 no.3
• /
• pp.84-88
• /
• 2001

For the elastic migation, the velocity errors between the initial velocity model and true velocity model seriously affect the migrated images. The assumption of an initial velocity model, thus, is one of the critical factor for the successful migration. In case of applying the layered earth model as an initial velocity model, the layer boundary having large velocity contrast can not be defined well with conventional traveltime calculation algolithms and we have the difficulties for expressing the characteristics of the real subsurface. Smooth Background Model (SBM) we have applied as an initial velocity model in our study is characterized to be linearly varying the velocity with the depth, which can express the velocity variation in the subsurface properly. Thus it can properly be applied to traveltime calculation algolithms such as Vidale's method. In this study, Kirchhoff operator for prestack migration was used and the absolute amplitude obtained by modeling was applied as a weighted value to consider the true amplitude for initial model. Initial velocity model for migration was determined by using stacking velocity and we applied this model to real data.

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

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.638-642
• /
• 2005

For gas hydrate exploration, long offset multichannel seismic data acquired using by the 4km streamer length in Ulleung basin of the East Sea. The dataset was processed to define the BSRs (Bottom Simulating Reflectors) and to estimate the amount of gas hydrates. Confirmation of the presence of Bottom Simulating reflectors (BSR) and investigation of its physical properties from seismic section are important for gas hydrate detection. Specially, faster interval velocity overlying slower interval velocity indicates the likely presences of gas hydrate above BSR and free gas underneath BSR. In consequence, estimation of correct interval velocities and analysis of their spatial variations are critical processes for gas hydrate detection using seismic reflection data. Using Dix's equation, Root Mean Square (RMS) velocities can be converted into interval velocities. However, it is not a proper way to investigate interval velocities above and below BSR considering the fact that RMS velocities have poor resolution and correctness and the assumption that interval velocities increase along the depth. Therefore, we incorporated Migration Velocity Analysis (MVA) software produced by Landmark CO. to estimate correct interval velocities in detail. MVA is a process to yield velocities of sediments between layers using Common Mid Point (CMP) gathered seismic data. The CMP gathered data for MVA should be produced after basic processing steps to enhance the signal to noise ratio of the first reflections. Prestack depth migrated section is produced using interval velocities and interval velocities are key parameters governing qualities of prestack depth migration section. Correctness of interval velocities can be examined by the presence of Residual Move Out (RMO) on CMP gathered data. If there is no RMO, peaks of primary reflection events are flat in horizontal direction for all offsets of Common Reflection Point (CRP) gathers and it proves that prestack depth migration is done with correct velocity field. Used method in this study, Tomographic inversion needs two initial input data. One is the dataset obtained from the results of preprocessing by removing multiples and noise and stacked partially. The other is the depth domain velocity model build by smoothing and editing the interval velocity converted from RMS velocity. After the three times iteration of tomography inversion, Optimum interval velocity field can be fixed. The conclusion of this study as follow, the final Interval velocity around the BSR decreased to 1400 m/s from 2500 m/s abruptly. BSR is showed about 200m depth under the seabottom