• The Journal of Engineering Geology
• /
• v.7 no.1
• /
• pp.27-36
• /
• 1997

This paper shows that crosshole seismic survey allows to detect even a small size of underground tunnel (about 2m$\times$2m). Such a small tunnel (e.g. infiltration tunnel) causes diffraction, as the seismic wave propagates, which results in distinctive variations of traveltime and amplitude of the first arrivals. This effect (or tunnel effect) is a typical indicator for the existence of tunnel and thereby an information about the tunnel location can be obtained. It was shown that the tunnel effect illustrated by numerical modeling (FDM) could be also observed in field measurements. The depth and shape of the tunnel were determined by a simplified processing method based on the use of amplitude variation of the first arrivals. The estimated location of the tunnel was well matched to that of the real tunnel.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.26-33
• /
• 2008

The analysis of wide-angle seismic reflection and refraction data plays an important role in lithospheric-scale crustal structure study. However, it is extremely difficult to develop an appropriate velocity structure model directly from the observed data, and we have to improve the structure model step by step, because the crustal structure analysis is an intrinsically non-linear problem. There are several subjective processes in wide-angle crustal structure modelling, such as phase identification and trial-and-error forward modelling. Because these subjective processes in wide-angle data analysis reduce the uniqueness and credibility of the resultant models, it is important to reduce subjectivity in the analysis procedure. From this point of view, we describe two software tools, PASTEUP and MODELING, to be used for developing crustal structure models. PASTEUP is an interactive application that facilitates the plotting of record sections, analysis of wide-angle seismic data, and picking of phases. PASTEUP is equipped with various filters and analysis functions to enhance signal-to-noise ratio and to help phase identification. MODELING is an interactive application for editing velocity models, and ray-tracing. Synthetic traveltimes computed by the MODELING application can be directly compared with the observed waveforms in the PASTEUP application. This reduces subjectivity in crustal structure modelling because traveltime picking, which is one of the most subjective process in the crustal structure analysis, is not required. MODELING can convert an editable layered structure model into two-way traveltimes which can be compared with time-sections of Multi Channel Seismic (MCS) reflection data. Direct comparison between the structure model of wide-angle data with the reflection data will give the model more credibility. In addition, both PASTEUP and MODELING are efficient tools for handling a large dataset. These software tools help us develop more plausible lithospheric-scale structure models using wide-angle seismic data.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.116-126
• /
• 2008

In gas hydrate survey, seismic amplitude and frequency characteristics play a very important role in determining whether gas hydrate exists. According to the variation of source frequency and scatterer size, we study seismic amplitude characteristics using elastic modeling applied at staggered grids. Generally speaking, scattering occurs in proportion to the square of source frequency and the scatterer volume, which has an effect on seismic amplitude. The higher source frequency is, the more scattering occurs in gas hydrate bearing zone. Therefore, BSR is hardly observed in high frequencies. On the other side, amplitude blanking zone and BSR is clearly observed in lower frequencies although the resolution is poor as a whole. Seismic reflections traveling through free-gas layer below gas hydrate bearing zone decay so severely a high frequency component that a low frequency term is dominant. Amplitude anomaly of BSR result from high acoustic impedance contrast due to free-gas, which is a very crucial factor to estimate gas hydrate bearing zone. Seismic frequency analysis is carried out using wavelet transform method that frequency component could be decomposed with time variation. In application of wavelet transform to the seismic physical experiments data, we can observe that reflections traveling through air layer, which corresponds to the free-gas layer, decay a high frequency component.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2011.10a
• /
• pp.122-122
• /
• 2011

• 한국지구물리탐사학회:학술대회논문집
• /
• 2010.10a
• /
• pp.79-82
• /
• 2010

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.61-69
• /
• 2000

The free-surface boundary conditions are persistent problem in elastic wave equation modeling by finite-difference method, which can be summarized with the degradation of the accuracy of the solution and limited stability range in Poisson's ratio. In this paper, we propose the mid-point averaging scheme as an alternative way of implementing the free-surface boundary conditions, and present the solution to Lamb's problem to verify our approach.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.183-191
• /
• 2009

Finite difference method using not general SSG (standard staggered grid) but RSG (rotated staggered grid) was applied to simulation of elastic wave propagation. Special free surface boundary condition such as imaging method is needed in finite difference method using SSG in elastic wave propagation. But free surface boundary condition in finite difference method using RSG is easily solved with adding air layer or vacuum layer. Recently PML (Perfectly Matched layer) is widely used to eliminate artificial reflection waves from finite boundary because of its' greate efficiency. Absorbing ability of CPML (convolutional Perfectly Matched Layer) that is more efficient than that of PML and CPML that don't use splitting of wave equation that should be adapted to PML was applied to FDM using RSG in this study. Frequency absorbing characteristic and energy absorbing ability in CPML layer were investigated and CPML eliminated artificial boundary waves very effectively in FDM using RSG in being compared with that of Cerjan's absorbing method. CPML method also diminished amplitude of waves in boundary layer of solid-liquid model very well.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.148-154
• /
• 2006

The effect of small-scale near surface inhomogeneity on Rayleigh wave propagation and dispersion has been investigated in this study using two-dimensional FEM elastic modeling. Various inhomogeneity models with a variety of geometrical shape and embedment depth which exist in homogeneous half-space and two-layered media are considered. Results show that any near surface inhomogeneity greater than one wavelength in terms of minimum wavelength of Rayleigh wave shows dispersion characteristics. Such dispersion effect become stronger as the dimensions of the inhomogeneity increase. The effect of horizontal dimension is more dominant factor governing the dispersion characteristics than vertical dimension. However, the dispersion effect can not be identifiable in seismogram if the horizontal dimension is not wide enough. Nonetheless, even in this case, the existence of inhomogeneity can be inferred by the reflection or transmission event of Rayleigh wave. The results can be expected to provide insights on the behavior of Rayleigh wave which may be helpful for designating field work or new processing scheme to detect near surface inhomogeneity by surface wave method.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.2
• /
• pp.5-17
• /
• 2007

Various field setup and filtering criteria have been suggested to avoid the near field effects in surface wave methods. Unlike other surface wave methods HWAW method uses the near field component positively. It is possible by using maximum energy point based on time-frequency map and inversion method to consider receiver locations from the source point and body wave component. To verify the HWAW method in the near field numerical study was performed and the wave propagation in the stratified soil media was simulated due to a surface point load. All of five representative soil models were used. The experimental dispersion curves, determined by HWAW method at the various receiver distances in the region of near field, all coincided well with the theoretical dispersion curves determined by 3D forward modeling (Kausel's method). Consequently, it was considered that the HWAW method can provide reliable $V_s$ profiles effectively in the near field.