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

• 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.3 no.1
• /
• pp.7-12
• /
• 2000

The paper discusses a data processing methodology that derives a three dimensional porosity volume information from the 3-D seismic dataset. The methodology consists of preprocessing and inversion procedures. The purpose of the preprocessing is balancing the amplitudes of seismic traces by using reflectivity series derived from sonic and density logs. There are eight sonic logs are available in the study area; therefore, we can compute only 8 balance functions. The balance function for every seismic trace was derived from these 8 balance functions by kriging. In order to derive a wide-band acoustic impedance --similar to the one can be derived from a sonic log- from a band-limited reflection seismogram, we need to recover missing low- and high-frequency information of the seismic trace. For that Purpose we use the autoregressive method.

• 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.2
• /
• pp.153-160
• /
• 2008

Seismic modeling is used to simulate wave propagation in the earth. Although the earth's subsurface is usually semi-infinite, we cannot handle the semi-infinite model in seismic modeling because of limited computational resources. For this reason, we usually assume a finite-sized model in seismic modeling. In that case, we need to eliminate the edge reflections arising from the artificial boundaries introducing a proper boundary condition. In this study, we changed three kinds of boundary conditions (sponge boundary condition, Clayton and Engquist's absorbing boundary condition, and Higdon's transparent boundary condition) so that they can be applied in elastic wave modeling for anisotropic media. We then apply them to several models whose Poisson's ratios are different. Clayton and Engquist's absorbing boundary condition is unstable in both isotropic and anisotropic media, when Poisson's ratio is large. This indicates that the absorbing boundary condition can be applied in anisotropic media restrictively. Although the sponge boundary condition yields good results for both isotropic and anisotropic media, it requires too much computational memory and time. On the other hand, Higdon's transparent boundary condition is not only inexpensive, but also reduce reflections over a wide range of incident angles. We think that Higdon's transparent boundary condition can be a method of choice for anisotropic media, where Poisson's ratio is large.

• 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.11 no.1
• /
• pp.1-14
• /
• 2008

We present a method for interpreting seismic records with arrivals and waveforms having characteristics which could be generated by extremely inhomogeneous velocity structures, such as non-typical oceanic crust, decollement at subduction zones, and seamounts in oceanic regions, by comparing them with synthetic waveforms. Recent extensive refraction and wide-angle reflection surveys in oceanic regions have provided us with a huge number of high-resolution and high-quality seismic records containing characteristic arrivals and waveforms, besides first arrivals and major reflected phases such as PmP. Some characteristic waveforms, with significant later reflected phases or anomalous amplitude decay with offset distance, are difficult to interpret using only a conventional interpretation method such as the traveltime tomographic inversion method. We find the best process for investigating such characteristic phases is to use an interactive interpretation method to compare observed data with synthetic waveforms, and calculate raypaths and traveltimes. This approach enables us to construct a reasonable structural model that includes all of the major characteristics of the observed waveforms. We present results here with some actual observed examples that might be of great help in the interpretation of such problematic phases. Our approach to the analysis of waveform characteristics is endorsed as an innovative method for constructing high-resolution and high-quality crustal structure models, not only in oceanic regions, but also in the continental regions.

• Journal of Industrial Technology
• /
• v.20 no.A
• /
• pp.303-309
• /
• 2000

Seismic reflection method is applied to detect shallow location of limestone in Chechon area. The data using hammer source is compared with that of weight drop. Small size hammer and weight-drop are used as energy source and 100Hz geophones are used for data aquisition. Data processing is conducted utilizing the available processing technique of "Geobit", which is seismic data processing software developed by KIGAM. The result of above data processing, the velocity of topsoil layer is 1,250m/sec. The velocity of this area is higher than other area because loading trucks pass this area and make this layer compact. And in limestone area, hammer is proposed to energy source instead of weight drop because the energy propagates the layer very well.

• Geophysics and Geophysical Exploration
• /
• v.19 no.2
• /
• pp.76-83
• /
• 2016

Migration image can be distorted due to reflected waves in the source and receiver wavefields when discontinuities of input velocity model exist in seismic imaging. To remove reflected waves coming from layer interfaces, it is a common practice to smooth the velocity model for migration. If the velocity model is smoothed, however, the subsurface image can be distorted because the velocity changes around interfaces. In this paper, we attempt to minimize the distortion by reducing reflection energy in the source and receiver wavefields through acoustic impedance homogenization. To make acoustic impedance constant, we define fake density model and use it for migration. When the acoustic impedance is constant over all layers, the reflection coefficient at normal incidence becomes zero and the minimized reflection energy results in the improvement of migration result. To verify our algorithm, we implement the reverse-time migration using cell-based finite-difference method. Through numerical examples, we can note that the migration image is improved at the layer interfaces with high velocity contrast, and it shows the marked improvement particularly in the shallow part.

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