• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.43-48
• /
• 1998

Recently seismic tomography has been widely used to visualize subsurface structure for resource explorations and construction site evaluation. We studied a way to include fresnel zone concept in the conventional ray-based traveltime tomography. The algorithm developed uses the same order of computing time as the conventional traveltime to mography but incorporates the rigorous wavepath concept of wave-equation tomography. Some experiments to synthetic and real data show reasonable results compared to conventional ray-based traveltime tomography.

• Economic and Environmental Geology
• /
• v.40 no.6
• /
• pp.739-749
• /
• 2007

The seismic velocity at the formation varies widely with physical properties in the layers. These features on seismic shot gathers are not capable of reproducing normally by numerical modeling of homogeneous medium, so that we need that of random inhomogeneous medium instead. In this study, we conducted Gaussian autocorrelation function (ACF), exponential autocorrelation function and von Karman autocorrelation function for getting inhomogeneous velocity model and applied a simple geological model. According to the results, von Karman autocorrelation function showed short wavelength to the inhomogeneous velocity medium. For numerical modeling for a gas hydrate, we determined a geological model based on field data set gathered in the East sea. The numerical modeling results showed that the von Karman autocorrelation function could properly describe scattering phenomena in the gas hydrate velocity model which contains an inhomogeneous layer. Besides, bottom-simulating-reflectors and scattered waves which appear at seismic shot gather of the field data showed properly in the inhomogeneous numerical modeling.

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

• Journal of the Korean Geophysical Society
• /
• v.2 no.4
• /
• pp.251-258
• /
• 1999

A high-resolution seismic survey was conducted on the abandoned tidal flat in Shihwa Lake on the west coast of Korea. A portable vibrator was used as a seismic source and 217 shot gathers with 48-channel system were acquired. F-k filtering, residual static corrections and post-stack frequency filtering are found to be useful for signal enhancement. The overburden is divided into two seismic depositional units. Unit I is deposited in tidal environment characterized by parallel and high continuity reflection pattern. This unit comprises a dry layer (Unit Ia) and a wet layer (Unit Ib) having averagely 5 and 15 meters thick, respectively. Unit II unconformably overlain by Unit I exhibits discontinuity and hummocky reflection pattern, indicating complex channel-fill sediments in estuary. The maximum thickness of this unit is approximately 20 meter. Acoustic basement is considered as Precambrian granitic gneiss occurred in the surrounding area.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.203-213
• /
• 2022

To understand an underlying geological structure via seismic imaging, the velocity information of the subsurface medium is crucial. Although the full-waveform inversion (FWI) method is considered useful for estimating subsurface velocity models, 3D FWI needs a lot-of computing power and time. Herein, we compare the calculation efficiency and accuracy of frequency-domain 2D and 3D acoustic FWIs. Thereafter, we demonstrate that the artifacts from 2D approximation can be partially suppressed via frequency-domain 2D FWI by employing diffraction angle filtering (DAF). By applying DAF, which employs only big reflection angle components, the impact of noise and out-of-plane reflections can be reduced. Additionally, it is anticipated that the DAF can create long-wavelength velocity structures for 3D FWI and migration.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.361-367
• /
• 2008

When constructing a synthetic seismogram in the earthquake study or in seismic data interpretation by using a ray-tracing technique, the most troublesome and error-prone task is to define a suite of ray codes for the corresponding rays to trace in advance. An infinite number of rays exist for any arbitrarily located source and receiver in a medium. Missing certain important rays or an inappropriate selection of ray codes in tracing rays may result in wrong interpretation of the earthquake record or seismogram. Automatic ray code generation could be able to eliminate those problems. In this study we have developed an efficient algorithm with which one can generate systematically all the ray codes for the source(s) and receiver(s) arbitrarily located in a model. The result of this work could be used not only in analysing multiples in seismic data processing and interpretation, but also in coda wave study, study on the amplification effects in a basin and phase identification of the waves multiply reflected/refracted in earthquake study.

• Geophysics and Geophysical Exploration
• /
• v.25 no.1
• /
• pp.14-25
• /
• 2022

The passive surface wave method measures seismic signals from ambient noises or vibrations of natural phenomena without using an artificial source. Since passive sources are usually in lower frequencies than artificial ones being able to ensure the information on deeper geological structures, the passive surface wave method can investigate deeper geological structures. In the passive method, frequency dispersion curves are obtained after data acquisition, and the dispersion curves are analyzed by assuming 1D-layered earth, which is like the method of active surface wave survey. However, when computing dispersion curves, the passive method first obtains and analyzes coherence curves of received signals from a set of receivers based on spatial autocorrelation. In this review, we explain how passive surface wave methods measure signals, and make data processing and interpretation, before analyzing field application cases.

• Geophysics and Geophysical Exploration
• /
• v.19 no.4
• /
• pp.220-227
• /
• 2016

Marine seismic data have not only primary signals from subsurface but also ghost signals reflected from the sea surface. The ghost decreases temporal resolution of seismic data because it attenuates specific frequency components. For eliminating the ghost signals effectively, the exact ghost delaytimes and reflection coefficients are required. Because of undulation of the sea surface and vertical movements of airguns and streamers, the ghost delaytime varies spatially and randomly while acquiring seismic data. The reflection coefficient is a function of frequency, incidence angle of plane-wave and the sea state. In order to estimate the proper ghost delaytimes considering these characteristics, we compared the ghost delaytimes estimated with L-1 norm, L-2 norm and kurtosis of the deghosted trace and its autocorrelation on synthetic data. L-1 norm of autocorrelation showed a minimal error and the reflection coefficient was calculated using Kirchhoff approximation equation which can handle the effect of wave height. We applied the estimated ghost delaytimes and the calculated reflection coefficients to remove the source and receiver ghost effects. By removing ghost signals, we reconstructed the frequency components attenuated near the notch frequency and produced the migrated stack section with enhanced temporal resolution.

• The Journal of Engineering Geology
• /
• v.26 no.3
• /
• pp.339-352
• /
• 2016

The best way of investigating the physical and mechanical properties of subsurface materials is the combined interpretation of data from borehole geophysical surveys and geotechnical experiments with rock samples. In this study two surface seismic surveys with refraction and surface-wave method are alternatively conducted for downhole seismic surveys in test site for geological field trip in Jeungpyung, Chungbuk. P- and S-wave velocity structures are delineated by refraction and MASW (multichannel analysis of shear waves) methods, respectively. Possion's ratio section, reconstructed from P- and S-wave velocities, is correlated to the outcrop geological features consisting of reddish sedimentary rock, gray volcanic rock, and joints/fractures. In addition, rock samples representative for reddish sedimentary and gray volcanic features are geotechnically analyzed to provide physical, mechanical properties, and elastic modulus. Dynamic elastic moduli estimated from geophysical data is found to be higher than the one from geotechnical data. Reddish sedimentary rock characterized with low porosity and moisture content corresponds to the zone of low electrical resistivities and their small variations in the resistivity sections between the rainy and dry days. This trend suggests that the weathered gray volcanic rock and the nearby fractures with higher low porosity and moisture content are interpreted to be good carrier especially in rainy season.

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