• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.18-22
• /
• 2002

As a single arrival traveltime, maximum energy arrival traveltime has been known as the most proper operator for Kirchhoff migration. In case of the model having the simple structure, both the first arrival traveltime and the maximum energy arrival traveltime can be used as the correct operators for Kirchhoff migration. However for some model having the complex and high velocity contrast structure, the migration using the first arrival traveltime can't give the correct depth section. That is, traveltime to be required in Kirchhoff migration is the maximum energy traveltime, but, needs considerably more calculation time than that of first arrival. In this paper, we propose the method for calculating the traveltime approximated to the maximum energy arrival using one-way wave equation. After defining the WAS(Wrap Around Suppression) factor to be used for calculating the first arrival traveltime using one-way wave equation as the function of lateral grid interval and depth and considering the delay time of source wavelet. we calculate the traveltime approximated to the maximum energy arrival. to verify the validity of this traveltime, we applied this to the migraion for simple structure and complex structure and compared the depth section with that obtained by using the first arrival traveltime.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.402-409
• /
• 2003

Inversion of traveltime requires an efficient algorithm for computing the traveltime as well as its $Frech\hat{e}t$ derivative. We compute the traveltime of the head waves using the damped wave solution in the Laplace domain and then present a new algorithm for calculating the $Frech\hat{e}t$ derivative of the head wave traveltimes by exploiting the numerical structure of the finite element method, the modem sparse matrix technology, and SWEET algorithm developed recently. Then, we use a properly regularized steepest descent method to invert the traveltime of the Marmousi-2 model. Through our numerical tests, we will demonstrate that the refraction tomography with large aperture data can be used to construct the initial velocity model for the prestack depth migration.

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

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms for the travel time computation commonly used, however, may not give sufficiently precise results particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. We assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution.

• The Journal of Engineering Geology
• /
• v.17 no.4
• /
• pp.633-642
• /
• 2007

Traveltime tomogram is generally used for interpretation of seismic tunnel data. In the field data, the first arrival traveltime is less dispersive with increasing source-receiver seperation compared to theoretical model data. So the result of calculation can be serious despite of small errors such as traveltime picking. In this study, amplitude method and error tomogram method are tried to overcome these problems. This method will help the interpretation of the data from the underground tunnel.

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

• Geophysics and Geophysical Exploration
• /
• v.7 no.4
• /
• pp.234-244
• /
• 2004

In order to calculate 3-dimensional wavefield using finite-element method in frequency domain, we must factor so huge sparse impedance matrix. Because of difficulties of handling of this huge impedance matrix, 3-dimensional wave equation modeling is conducted mainly in time domain. In this study, we simulate the 3-D wavefield using finite-element method in Laplace domain by combining high-performance parallel finite-element solver and SWEET (Suppressed Wave Equation Estimation of Traveltime) algorithm which can calculate the traveltime and the amplitude. To verify this combination, we applied it to the SEG/EAGE 3D salt model in serial and parallel computing environments.

• Tunnel and Underground Space
• /
• v.18 no.3
• /
• pp.202-213
• /
• 2008

In this study, theoretical approach of 3D seismic traveltime tomography was investigated. To guarantee the successful field application of 3D tomography, appropriate control of problem associated with blind zone is pre-requisite. To overcome the velocity distortion of the reconstructed tomogram due to insufficient source-receiver array coverage, the algorithm of 3D seismic traveltime tomography based on the Fresnel volume was developed as a technique of ray-path broadening. For the successful reconstruction of velocity cube, 3D traveltime algorithm was explored and employed on the basis of 2nd order Fast Marching Method(FMM), resulting in improvement of precision and accuracy. To prove the validity and field application of this algorithm, two numerical experiments were performed for globular and layered models. The algorithm was also found to be successfully applicable to field data.

• Economic and Environmental Geology
• /
• v.25 no.1
• /
• pp.79-85
• /
• 1992

An Interpretation technique is presented to determine strike, dip, velocity and depth of 3-D planar layers from refraction or reflection traveltime curve. This interpretation technique determines the direction of emerging ray from the slope of the traveltime curve and traces the emerging ray to the refractor or reflector. The ray direction in the last layer is used to decide the normal vector to the refractor or reflector from whick its dip, strike and velocity are calculated. The vertical depth to the refractor or reflector is computed by using the intercept or zero-offset time and the ray direction in each layer. Some tests on the interpretation method are performed for the sysnthetic traveltimes generated in 3-D model layers and show that the paramerters of the model layers are accurately determined.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.46-55
• /
• 2002

We investigated the distorting factors of velocity structure reconstructed by traveltime inversion. The set of models that fit the data in a numerical sense usually contains unrealistic models. Reconstructed velocity structure was enhanced because unreasonable models were eliminated by defining constraint of variable grid using a priori information. To correct time delay of source explosion, which distorts traveltime tomograms, terms for correction of time delay was formulated into equation of travel time tomography.

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