• Journal of the Korean Geophysical Society
• /
• v.2 no.3
• /
• pp.179-190
• /
• 1999

The overall structural framework was studied using the regional 2D seismic data, followed by the sequence stratigraphic study on the 3D seismic and well- log data in the margin of the South Con Son basin of the South China Sea. This research contributes to delineate depositional stratigraphy, depositional environment and geologic history in the 3D seismic area of highly complicated faulting. Eight Miocene sequences were indicated on the 3D seismic and well-log data, in which the structural maps of each sequence boundary and the isochron maps for the corresponding sequence were made. The seismic facies were analyzed for each sequence volume and sequence boundary surface. The 3D seismic area is characterized by coal beds deposited in the transgression environment (transgression systems tract) and channel distributions just above the sequence boundaries. During the Early Miocene, the coals and thick shales deposited in the mangrove swamp representing the lower coastal plain environment. During the Mid to Late Miocene, thick clastic sediments deposited in the coastal to shallow shelf by regional subsidence and marine transgression. The isochron maps and structural patterns indicate that the sediments were transported from west to east or from northwest to southeast.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.13 no.6
• /
• pp.490-500
• /
• 1988

In recent years the acquisition, processing and interpretation of three dimensional seisimic data, for the purpose of locating gas and reservoirs, have become practical. This paper exlores one way in which the volume data can be searched and visualized, which may aid the interpreter. The illusion of looking at a three dimensional volume can be obrained by fusing a stereoscopic pair of pictures. Each picture can be made by projecting each data point of the volume into a plane from a point where the eye is placed. The data valuse along any projection line can be summed to form the picture, or only a segment along the line can be selected. By selective projection, the volume can be searched and obscuring layers removed. The stereoscopic pictures show the physical models in there ture spatial positions. Projection of the envelope function of the seismic traces is shown to give improved depth perception compared with projection of the position amplitudes.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.24-30
• /
• 2010

This study utilises repeated numerical tests to understand the effects of variable near-surface conditions on time-lapse seismic surveys. The numerical tests were aimed at reproducing the significant scattering observed in field experiments conducted at the Naylor site in the Otway Basin for the purpose of $CO_2$ sequestration. In particular, the variation of elastic properties of both the top soil and the deeper rugose clay/limestone interface as a function of varying water saturation were investigated. Such tests simulate the measurements conducted in dry and wet seasons and to evaluate the contribution of these seasonal variations to seismic measurements in terms of non-repeatability. Full elastic pre-stack modelling experiments were carried out to quantify these effects and evaluate their individual contributions. The results show that the relatively simple scattering effects of the corrugated near-surface clay/limestone interface can have a profound effect on time-lapse surveys. The experiments also show that the changes in top soil saturation could potentially affect seismic signature even more than the corrugated deeper surface. Overall agreement between numerically predicted and in situ measured normalised root-mean-square (NRMS) differences between repeated (time-lapse) 2D seismic surveys warrant further investigation. Future field studies will include in situ measurements of the elastic properties of the weathered zone through the use of 'micro Vertical Seismic Profiling (VSP)' arrays and very dense refraction surveys. The results of this work may impact on other areas not associated with $CO_2$ sequestration, such as imaging oil production over areas where producing fields suffer from a karstic topography, such as in the Middle East and Australia.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.279-285
• /
• 2008

Korea Institute of Geoscience & Mineral Resources (KIGAM) carried out 2 dimensional multi-channel seismic surveys for Domi-Basin of east-southern part of Jeju Island, South Sea, Korea in 2007. The purpose of this survey is to investigate the structure of acoustic basement and the potential of energy resources in the Korean shelf. It is essential to produce fine stack and migration section to understand the structure of basement. However a basement can not be clearly defined where multiples exist between sea surface and seafloor. This study aimed at designing the optimal data processing parameter, especially to eliminate the peg-leg multiples. Main data processing procedure is composed of minimum phase predictive deconvolution, velocity analysis and Radon filter. We tested the efficiency of processing parameter from stack sections of each step. Our results confirmed that processing parameters are suitable for the seismic data of Domi-Basin.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.13-21
• /
• 2020

The acquisition and processing of long-offset data are essential for imaging deep geological structures in marine seismic surveys. It is challenging to derive an accurate subsurface image by employing conventional data processing to long-offset data owing to the normal moveout (NMO) stretch and non-hyperbolic moveout phenomena induced by seismic anisotropy. In 2017, the Korea Institute of Geoscience and Mineral Resources conducted a simultaneous two-dimensional multichannel streamer and ocean-bottom seismic survey using a 5.7-km streamer and an ocean-bottom seismometer to identify the deep geological structure of the Ulleung Basin. Herein, the actual geological subsurface structure was obtained via the sequential iterative updating of the velocity and anisotropic parameters of the long-offset data obtained using a multichannel streamer, and anisotropic prestack Kirchhoff migration was performed using the updated velocity and anisotropic parameters as input parameters. As a result, the reflection energy in the long-offset traces, which showed non-hyperbolic moveout owing to seismic anisotropy, was well aligned horizontally and NMO stretches were also reduced. Thus, a more precise and accurate migrated image was obtained, minimizing the distortion of reflectors and mispositioned reflection energy.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.1-18
• /
• 2007

Non-linear elastic wavefield inversion is a powerful method for estimating elastic parameters for physical constraints that determine subsurface rock and properties. Here, I introduce six elastic-wave velocity models by reconstructing elastic-wave velocity variations from real data and a 2D elastic-wave velocity model. Reflection seismic data information is often decoupled into short and long wavelength components. The local search method has difficulty in estimating the longer wavelength velocity if the starting model is far from the true model, and source frequencies are then changed from lower to higher bands (as in the 'frequency-cascade scheme') to estimate model elastic parameters. Elastic parameters are inverted at each inversion step ('simultaneous mode') with a starting model of linear P- and S-wave velocity trends with depth. Elastic parameters are also derived by inversion in three other modes - using a P- and S-wave velocity basis $('V_P\;V_S\;mode')$; P-impedance and Poisson's ratio basis $('I_P\;Poisson\;mode')$; and P- and S-impedance $('I_P\;I_S\;mode')$. Density values are updated at each elastic inversion step under three assumptions in each mode. By evaluating the accuracy of the inversion for each parameter set for elastic models, it can be concluded that there is no specific difference between the inversion results for the $V_P\;V_S$ mode and the $I_P$ Poisson mode. The same conclusion is expected for the $I_P\;I_S$ mode, too. This gives us a sound basis for full wavelength elastic wavefield inversion.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.184-196
• /
• 2008

We perform MT soundings in Bendigo, the northern part of Victoria, Australia, to investigate the deep subsurface geologic structure. The primary purpose of this survey is to figure out whether the discontinuity such as faults extends northward. The time series of MT signal were measured over 11 days at 71 measurement stations together with at remote reference, which help enhance the quality of impedance estimation and its interpretation. The impedances are estimated by robust processing using remote reference technique and then inverted with 2D MT 2D inversion. We can see that known faults are clearly imaged in MT 2D inversion. Comparing resistivity images from MT 2D inversion with interpreted boundary from reflection seismic exploration, two interpretations match well each other.

• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.78-88
• /
• 2021

Owing to the abnormal weather conditions due to global warming, carbon capture and storage (CCS) technology has attracted global attention as a countermeasure to reduce CO₂ emissions. In the Pohang CCS demonstration project in South Korea, 100 tons of CO₂ were successfully injected into the subsurface CO₂ storage in early 2017. However, after the 2017 Pohang earthquake, the Pohang CCS demonstration project was suspended due to an increase in social concerns about the safety of the CCS project. In this study, to reconfirm the structural suitability of the CO₂ storage site in the Pohang Basin, we employed seismic imaging based on reverse-time migration (RTM) to analyze small-scale ocean-bottom seismic data, which have not been utilized in previous studies. Compared with seismic images using marine streamer data, the continuity of subsurface layers in the RTM image using the ocean-bottom seismic data is improved. Based on the obtained subsurface image, we discuss the structural suitability of the Pohang CO₂ storage site.