• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2006.03a
• /
• pp.25-32
• /
• 2006

We present a imaging method of seismic sources by time reversal propagation of seismic waves. Time-reversal wave propagation is actively used in medical imaging, non destructive testing and waveform tomography. Time-reversal wave propagation is based on the time-reversal invariance and the spatial reciprocity of the wave equation. A signal is recorded by an array of receivers, time-reversed and then back-propagated into the medium. The time-reversed signal propagates back into the same medium and the energy refocuses back at the source location. The increasing power of computers and numerical methods makes it possible to simulate more accurately the propagation of seismic waves in heterogenous media. In this work, a staggered-grid finite-difference solution of the elastic wave equation is employed for the wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• Geomechanics and Engineering
• /
• v.8 no.3
• /
• pp.377-390
• /
• 2015

The governing equations of wave propagation in one dimension of elastic continuum materials are investigated by taking the influence of the initial stress into account. After a short review of the theory of elastic wave propagation in a rock mass with an initial stress, results indicate that the initial stress differentially influences P-wave and S-wave propagation. For example, when the initial stress is homogeneous, for the P-wave, the initial stress only affects the magnitude of the elastic coefficients, but for the S-wave, the initial stress not only influences the elastic coefficients but also changes the governing equation of wave propagation. In addition, the P-wave and S-wave velocities were measured for granite samples at a low initial stress state; the results indicate that the seismic velocities increase with the initial stress. The analysis of the previous data of seismic velocities and elastic coefficients in rocks under ultra-high hydrostatic initial stress are also investigated.

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

There have been several studies about empirical relation between seismic source parameters(e.g., focal mechanisms, depths, magnitudes, etc.) and T-wave observation. In order to delineate the relation, numerical and theoretical approaches to figure out T-wave excitation mechanism are required. In an attempt to investigate source radiation and wave scattering effects in the oceanic crust on T-wave envelopes, we perform three-dimensional numerical modeling to synthesize T-wave envelopes. We first calculate seismic P- and SV-wave energy on the seafloor using the Direct Simulation Monte Carlo based on the Radiative Transfer Theory, which enables us to take into account both realistic seismic source parameters and wave scattering in heterogeneous media, and then estimate excited T-wave energy by normal mode computation. The numerical simulation has been carried out considering the following different conditions: source types (strike and normal faults), source depths (shallow and deep), and wave propagation through homogeneous and heterogeneous Earth media. From the results of numerical modeling, we confirmed that T-wave envelopes vary according to spatial seismic energy distributions on the seafloor for the various input parameters. Furthermore, the synthesized T-wave envelopes show directional patterns due to anisotropic source radiation, and the slope change of T-wave envelopes caused by focal depth. Seismic wave scattering in the oceanic crust is likely to control the shape of envelopes.

• Journal of the Korean earth science society
• /
• v.27 no.5
• /
• pp.539-547
• /
• 2006

The activity of the seismic wave propagation around the cavity is investigated for the exact inversion of the crosshole tomography data, in order to understand the possibility of the existence inside the underground cavity. It is found that the adequate frequency range for the tunnel investigation is about 2 kHz to 5 kHz, and the grid space should be set up to 1/10 of the wavelength. The propagation of the seismic wave near the cavity may go through or detour the cavity according to the seismic velocity inside the cavity. The detouring wave propagates with the seismic velocity of mother rock in spite of the velocity of inside of the cavity. The smaller the velocity difference is between the mother rock and cavity, the more frequent penetration of the seismic wave through the cavity appears.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.81-88
• /
• 2000

The elastic wave equation is solved using the finite-difference method in 3D space to simulate the seismic wave propagation. It is based on the velocity-stress formulation of the equation of motion on a staggered grid. The nonreflecting boundary conditions are used to attenuate the wave field close to the numerical boundary. To satisfy the stress-free conditions at the free-surface boundary, a new formulation combining the zero-stress formalism with the vacuum one is applied. The effective media parameters are employed to satisfy the traction continuity condition across the media interface. With use of the moment-tensor components, the wide range of source mechanism parameters can be specified. The numerical experiments are carried out in order to test the applicability and accuracy of this scheme and to understand the fundamental features of the wave propagation under the generalized elastic media structure. Computational results show that the scheme is sufficiently accurate for modeling wave propagation in 3D elastic media and generates all the possible phases appropriately in under the given heterogeneous velocity structure. Also the characteristics of the ground motion in an sedimentary basin such as the amplification, trapping, and focusing of the elastic wave energy are well represented. These results demonstrate the use of this simulation method will be helpful for modeling the ground motion of seismological and engineering purpose like earthquake hazard assessment, seismic design, city planning, and etc..

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.181-186
• /
• 2006

An imaging method of seismic sources using time-reversal wave propagation is presented. The method is based on the time-reversal invariance and the spatial reciprocity of the wave equation. Time-reversal wave propagation has been used to image anomalous features of a midium in medical imaging, non destructive testing and waveform tomography. Seismogram is the record whose energy is propagated from the seismic source. If time-reversed seismogram propagates back into the medium, seismic energy is concentrated at the origin time of the event and at the source location. In this work, a staggered-grid finite-difference method of the elastic wave equation is parallelized for 3-D wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.291-296
• /
• 2006

The aspect of wave propagation around cavity was investigated for the exact inversion of crosshole tomography data in order to understand the possibility of the existence of underground cavity. We found that the adequate frequency range for the tunnel investigation was about 2kHz to 5kHz, and the grid space was set up to 1/10 length of wavelength. The propagation of the seismic wave near the cavity may go through or detour the cavity according to the seismic velocity of inside of cavity. The detouring wave propagates with the seismic velocity of mother rock in spite of the velocity of inside of cavity. The smaller the velocity difference between the mother rock and cavity, the more frequent penetration of the seismic wave through the cavity was appeared.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1999.10a
• /
• pp.91-98
• /
• 1999

In this study a modified ground strain model is developed for an equivalent earthquake load and is applied to the seismic analysis of buried pipelines, The ground strain can be obtained using the ratio of a maximum ground velocity to a wave propagation velocity. To reflect soil conditions and seismic characteristics the wave propagation velocity is evaluated by a proposed dispersion curve based on wave energy distribution. In order to verify the procedures the observed earthquake data and the results of this study are compared. For the application of an equivalent earthquake load to the seismic analysis the buried pipelines are modeled using the beam theory. the results of the analyses are compared with those of a dynamic analysis code and those obtained from the response displacement method. Finally various parametric studies considering different soil conditions and seismic loads are examined.

• Journal of Industrial Technology
• /
• v.27 no.A
• /
• pp.9-14
• /
• 2007

This study was designated to clarify the aspect of the wave propagation around the cavity. The change of traveltime and amplitude of the seismic wave was observed according to the various wave velocities of the cavity. The seismic wave detour or penetrate the cavity depending on the seismic velocity of the in-filled material. Generally, seismic wave detours toward high velocity zone around the cavity, and when the velocity of the cavity material reaches to 80 % of the base rock, the wave penetrates the cavity. The traveltime of the detouring seismic wave is not sensitive to the change of the cavity velocity, but as the velocity of the cavity increases, the fall of the amplitude was reduced. The penetrating wave showed the steeply increasing amplitude due to the reiteration of the detouring wave.

• Geomechanics and Engineering
• /
• v.12 no.1
• /
• pp.161-183
• /
• 2017

This paper investigates the damage identification of the concrete pile element through axial wave propagation technique using computational and experimental studies. Now-a-days, concrete pile foundations are often common in all engineering structures and their safety is significant for preventing the failure. Damage detection and estimation in a sub-structure is challenging as the visual picture of the sub-structure and its condition is not well known and the state of the structure or foundation can be inferred only through its static and dynamic response. The concept of wave propagation involves dynamic impedance and whenever a wave encounters a changing impedance (due to loss of stiffness), a reflecting wave is generated with the total strain energy forked as reflected as well as refracted portions. Among many frequency domain methods, the Spectral Finite Element method (SFEM) has been found suitable for analysis of wave propagation in real engineering structures as the formulation is based on dynamic equilibrium under harmonic steady state excitation. The feasibility of the axial wave propagation technique is studied through numerical simulations using Elementary rod theory and higher order Love rod theory under SFEM and ABAQUS dynamic explicit analysis with experimental validation exercise. Towards simulating the damage scenario in a pile element, dis-continuity (impedance mismatch) is induced by varying its cross-sectional area along its length. Both experimental and computational investigations are performed under pulse-echo and pitch-catch configuration methods. Analytical and experimental results are in good agreement.