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

• Economic and Environmental Geology
• /
• v.39 no.6 s.181
• /
• pp.711-717
• /
• 2006

In order to study gas hydrate, potential future energy resources, Korea Institute of Geoscience and Mineral Resources has conducted seismic reflection survey in the East Sea since 1997. one of evidence for presence of gas hydrate in seismic reflection data is a bottom simulating reflector (BSR). The BSR occurs at the interface between overlaying higher velocity, hydrate-bearing sediment and underlying lower velocity, free gas-bearing sediment. That is often characterized by large reflection coefficient and reflection polarity reverse to that of seafloor reflection. In order to apply depth migration to seismic reflection data. we need high performance computers and a parallelizing technique because of huge data volume and computation. Phase shift plus interpolation (PSPI) is a useful method for migration due to less computing time and computational efficiency. PSPI is intrinsically parallelizing characteristic in the frequency domain. We conducted conventional data processing for the gas hydrate data of the Ease Sea and then applied prestack depth migration using message-passing-interface PSPI (MPI_PSPI) that was parallelized by MPI local-area-multi-computer (MPI_LAM). Velocity model was made using the stack velocities after we had picked horizons on the stack image with in-house processing tool, Geobit. We could find the BSRs on the migrated stack section were about at SP 3555-4162 and two way travel time around 2,950 ms in time domain. In depth domain such BSRs appear at 6-17 km distance and 2.1 km depth from the seafloor. Since energy concentrated subsurface was well imaged we have to choose acquisition parameters suited for transmitting seismic energy to target area.

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

• 한국신재생에너지학회:학술대회논문집
• /
• 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.

• 한국신재생에너지학회:학술대회논문집
• /
• 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.5 no.3
• /
• pp.185-192
• /
• 2002

We made a study of 3-D migration which could precisely image data of GPR (Ground Penetrating Radar) applied to NDT (Non-Destructive Test) field for the inspection of structural safety. In this study, we obtained 3-D migrated images of important targets in structuresurvey (e.g. steel pipes, cracks) by using 3-D Kirchhoff prestack depth migration scheme developed for seismic data processing. For a concrete model consisting of steel pipe and void, the targets have been well defined with opposite amplitude according to the parameters of the targets. And migrated images using Parallel-Broadside array (XX configuration) have shown higher resolution than those using Perpendicular-Broadside array (YY configuration) when steel pipes had different sizes. Therefore, it is required to analyze the migrated image of XX configuration as well as that of general YY configuration in order to get more accurate information. As the last stage, we chose a model including two steel pipes which cross each other. The upper pipe has been resolved clearly but the lower has been imaged bigger than the model size due to the high conductivity of the upper steel.

• Geophysics and Geophysical Exploration
• /
• v.22 no.4
• /
• pp.225-238
• /
• 2019

Migration velocity analysis (MVA) for creating optimum depth-domain velocities in seismic imaging was applied to marine long-offset multi-channel data, and the effectiveness of the MVA approach was demonstrated by the combinations of conventional data processing procedures. The time-domain images generated by conventional time-processing scheme has been considered to be sufficient so far for the seismic stratigraphic interpretation. However, when the purpose of the seismic imaging moves to the hydrocarbon exploration, especially in the geologic modeling of the oil and gas play or lead area, drilling prognosis, in-place hydrocarbon volume estimation, the seismic images should be converted into depth domain or depth processing should be applied in the processing phase. CMP-based velocity analysis, which is mainly based on several approximations in the data domain, inherently contains errors and thus has high uncertainties. On the other hand, the MVA provides efficient and somewhat real-scale (in depth) images even if there are no logging data available. In this study, marine long-offset multi-channel seismic data were optimally processed in time domain to establish the most qualified dataset for the usage of the iterative MVA. Then, the depth-domain velocity profile was updated several times and the final velocity-in-depth was used for generating depth images (CRP gather and stack) and compared with the images obtained from the velocity-in-time. From the results, we were able to confirm the depth-domain results are more reasonable than the time-domain results. The spurious local minima, which can be occurred during the implementation of full waveform inversion, can be reduced when the result of MVA is used as an initial velocity model.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.18-22
• /
• 2002

As a single arrival traveltime, maximum energy arrival traveltime has been known as the most proper operator for Kirchhoff migration. In case of the model having the simple structure, both the first arrival traveltime and the maximum energy arrival traveltime can be used as the correct operators for Kirchhoff migration. However for some model having the complex and high velocity contrast structure, the migration using the first arrival traveltime can't give the correct depth section. That is, traveltime to be required in Kirchhoff migration is the maximum energy traveltime, but, needs considerably more calculation time than that of first arrival. In this paper, we propose the method for calculating the traveltime approximated to the maximum energy arrival using one-way wave equation. After defining the WAS(Wrap Around Suppression) factor to be used for calculating the first arrival traveltime using one-way wave equation as the function of lateral grid interval and depth and considering the delay time of source wavelet. we calculate the traveltime approximated to the maximum energy arrival. to verify the validity of this traveltime, we applied this to the migraion for simple structure and complex structure and compared the depth section with that obtained by using the first arrival traveltime.

• The Korean Journal of Petroleum Geology
• /
• v.11 no.1 s.12
• /
• pp.9-17
• /
• 2005

In the main target area of the block II, Targe-scale faults occur below the unconformity developed around 1 km in depth. The contrast of seismic velocity around the unconformity is generally so large that the strong multiples and the radical velocity variation would deteriorate the quality of migrated section due to serious distortion. More than 15 kinds of data processing techniques have been applied to improve the image resolution for the structures farmed from this active crustal activity. The bad and noisy traces were edited on the common shot gathers in the first step to get rid of acquisition problems which could take place from unfavorable conditions such as climatic change during data acquisition. Correction of amplitude attenuation caused from spherical divergence and inelastic attenuation has been also applied. Mild F/K filter was used to attenuate coherent noise such as guided waves and side scatters. Predictive deconvolution has been applied before stacking to remove peg-leg multiples and water reverberations. The velocity analysis process was conducted at every 2 km interval to analyze migration velocity, and it was iterated to get the high fidelity image. The strum noise caused from streamer was completely removed by applying predictive deconvolution in time space and ${\tau}-P$ domain. Residual multiples caused from thin layer or water bottom were eliminated through parabolic radon transform demultiple process. The migration using curved ray Kirchhoff-style algorithm has been applied to stack data. The velocity obtained after several iteration approach for MVA (migration velocity analysis) was used instead or DMO for the migration velocity. Using various testing methods, optimum seismic processing parameter can be obtained for structural and stratigraphic interpretation in the Block II, Yellow Sea Basin.