• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.09a
• /
• pp.3-7
• /
• 2003

The velocity structure beneath the CHNB broadband station is determined by receiver function analysis using by from teleseismic P waveforms. The detailed broadband receiver functions are obtained by stacking method for source-equalized vertical, radial and tangential components of teleseismic P waveforms. A time domain inversion uses the stacked radial receiver function to determine vertical P wave velocity structure beneath the station. The crustal velocity structures beneath the stations are estimated using the receiver function inversion method in the case at the crustal model parameterized by many thin, flat-lying, homogeneous layers. Events divide into 4 groups. four azimuths corresponding to events in group a(southwest), b(south), c(southeast), d(northeast). The result of crust at model inversion shows the crustal velocity structure beneath the CHNB station varies smoothly with increasing depth. The conard discontinuity lies around 18 km and moho discontinuity lies range from 30 to 34 km.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.173-182
• /
• 2009

In April 2008, KIGAM carried out an ocean-bottom seismometer (OBS) survey in the central Ulleung Basin where strong bottom simulating reflectors (BSRs) were revealed from previous surveys and some gas-hydrate samples were retrieved by direct sampling. The purpose of this survey is to estimate the velocity structure near the BSR in the gas hydrate prospect area using wide-angle seismic data recorded on the ocean-bottom seismometers. Along with the OBS survey, a 2-D seismic survey was performed whereby stratigraphic and preliminary velocity information was obtained. Two methods were applied to wide-angle data for estimating P wave velocity; one is velocity analysis in the $\tau$-p domain and the other is seismic traveltime inversion. A 1-D interval velocity profile was obtained by the first method, which was refined to layered velocity structure by the latter method. A layer stripping method was adopted for modeling and inversion. All velocity profiles at each OBS site clearly show velocity reversal at BSR depths due to the presence of gas hydrates. In addition, we could confirm high velocity in the column/chimney structure.

• Economic and Environmental Geology
• /
• v.48 no.1
• /
• pp.25-39
• /
• 2015

$CO_2$ geological storage plays an important role in reduction of greenhouse gas emissions, but there is a lack of research for CCS demonstration. To achieve the goal of CCS, storing $CO_2$ safely and permanently in underground geological formations, it is essential to understand the characteristics of them, such as total storage capacity, stability, etc. and establish an injection strategy. We perform the impedance inversion for the seismic data acquired from the Ulleung Basin in 2012. To review the possibility of $CO_2$ storage, we also construct porosity models and extract attributes of the prospects from the seismic data. To improve the quality of seismic data, amplitude preserved processing methods, SWD(Shallow Water Demultiple), SRME(Surface Related Multiple Elimination) and Radon Demultiple, are applied. Three well log data are also analysed, and the log correlations of each well are 0.648, 0.574 and 0.342, respectively. All wells are used in building the low-frequency model to generate more robust initial model. Simultaneous pre-stack inversion is performed on all of the 2D profiles and inverted P-impedance, S-impedance and Vp/Vs ratio are generated from the inversion process. With the porosity profiles generated from the seismic inversion process, the porous and non-porous zones can be identified for the purpose of the $CO_2$ sequestration initiative. More detailed characterization of the geological storage and the simulation of $CO_2$ migration might be an essential for the CCS demonstration.

• 한국신재생에너지학회:학술대회논문집
• /
• 2008.05a
• /
• pp.606-609
• /
• 2008

Deconvolution is one of the most used techniques for processing seismic reflection data. It is applied to improve temporal resolution by wavelet shaping and removal of short period reverberations. Several deconvolution algorithms such as predicted, spike, minimum entropy deconvolution and so on has been proposed to obtain such above purposes. Among of them, $\iota_1$ norm proposed by Taylor et al., (1979) and used to compared to minimum entropy deconvolution by Sacchi et al., (1994) has given some advantages on time computing and high efficiency. Theoritically, the deconvolution can be considered as inversion technique to invert the single seismic trace to the reflectivity, but it has not been successfully adopted due to noisy signals of the real data set and unknown source wavelet. After stacking, the seismic traces are moved to zero offset, thus each seismic traces now can be a single trace that is created by convolving the seismic source wavelet and reflectivity. In this paper, the fundamental of $\iota_1$ norm deconvolution method will be introduced. The method will be tested by synthetic data and applied to improve the stacked section of gas hydrate.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.303-306
• /
• 2007

After Kyoto protocol took effect, many countries are making efforts to reduce $CO_2$ one of effective which is geosequestration. But a percentage of geosequestration in total research budget is very small and the priority order of research also is receded in Korea. As one of efforts to activate the research on geosequestration in field of geophysics, we proposed the plan to build up a measurement system for rock physical properties monitoring during $CO_2$ injection which will function as original technology. The system consists of two part, one of which is a data acquisition system based on seismic and complex resistivity tomographic measurement and the other of which is a tri-axial compressive system to realize the in-situ condition. And also developments of various inversion algorithms are proposed to interpret data qualitatively such as a inversion algorithm for confined cylindrical boundary, a joint inversion algorithm and a 4-D inversion algorithm.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2011.10a
• /
• pp.49-53
• /
• 2011

• Tunnel and Underground Space
• /
• v.18 no.4
• /
• pp.263-271
• /
• 2008

3D seismic tomography technique was assessed for applicability of developed 3D tomography algorithm based on Fresnel volume in the dam-planned area. Reconstructed 3D tomogram based on Fresnel volume and Fast Marching Method(FMM) reveals similar velocity structure to the other geotechnical survey results. With the correlation analysis between RMR data and seismic velocity information, it could provide reliable information of rock mass rate. The applicability of 3D seismic tomography was verified in this study. It would be expected to apply 3D tomography with new developed first arrival calculation and inversion algorithm to the engineering field economically.

• The Journal of Engineering Geology
• /
• v.3 no.3
• /
• pp.231-239
• /
• 1993

This paper presents a result on the determination of velocity distribution by a tomographic inversion of crosshole seismic traveltimes in transversely isotropic(aniso tropic) media. The crosshole traveltimes used in this study are synthetic ones computed by ray tracing for some models having isotropic and transversely isotropic velocity distributions. The traveltimes are inverted by a general ART and ansotropic ART which considers the transversely isotropic effect during inver sion. The aniotropic ART gives accurate velodty distributions of transversely isotropic and isotropic models, while the isotropic ART determines accurate velocities only for the isotropic model but inaccurate for the transversely isotropic one. Therefore, the anisotropic ART may be used in case where no information is known on the isotropy or transverse isotropy of a survey area.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2005.03a
• /
• pp.18-25
• /
• 2005

Stochastic ground-motion method is adopted to simulate horizontal PGA values for the offshore Uljin earthquake recorded at nationwide seismic stations. For this purpose, the Fourier spectra are calculated at every stations based on comprehensive results of wave propagation and site effect which were previously revealed through inversion process applied to large accumulated spectral D/B. In addition, the apparent source spectrum of the offshore Uljin earthquake is estimated by removing the path and site response from the observed spectra. The distance dependent time-duration model is revised by iteratively fitting the PGA values generated by using the raw spectra data to the observed PGA data. The stochastic ground-motion method predicts the observed PGA values within a error of ${\sigma}_{log10}=0.1$. Transfer functions of a site relative to another site are estimated based on the error residual of the inversion results and used to convert PGA values at multiple stations to expected PGA values at a reference station of TJN. The converted PGA values can be used as basic data to evaluate the ground-motion attenuation relations developed for seismic hazard analysis that concerns the large damaging earthquakes.

• Economic and Environmental Geology
• /
• v.28 no.5
• /
• pp.481-485
• /
• 1995

This paper presents a new approaching method to determine the velocity and geometry of shallow subsurface from seismic refraction events. After picking the first breaks from seismic refraction data, we assume that field refraction seismogram can be replaced by the unit delta function having time shift of first break. Time curves are generated by shooting ray tracing. The partial derivatives seismogram for a damped least squares method is computed analytically at each step of the forward ray tracing. The technique is successfully tested on synthetic and real data. It has the advantage of real full waveform inversion, which is robust at low frequency band even if the initial guess is far from the true model.