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

• Journal of Industrial Technology
• /
• v.23 no.A
• /
• pp.63-68
• /
• 2003

Non-destructive method by tomography for safety diagnosis of civil engineering and building structures is tried. There are traveltime tomography that uses traveltime and fullwave tomography that uses the initial shock and seismic amplitude. But these methods have difficulty and weak points in accuracy and selection of initial value. In this study, corrected inversion method which is able to solve the two difficulty and this method is applied to theoretical pier model to calculate tomography.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

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

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

In Korean geology that crystalline rock is dominant, the properties of subsurface including the anisotropy are distributed complexly and changed abruptly. Because of such geological environments, cross-hole seismic traveltime tomography is widely used to obtain the high resolution image of the subsurface for the engineering purposes in the geotechnical sites. However, because the cross-hole tomography has a wide propagation angle coverage relatively, its data tend to include the seismic velocity anisotropy comparing with the surface seismic methods. It can cause the misinterpretation that the cross-hole seismic data including the anisotropic effects are analyzed and treated with the general processing techniques assuming the isotropy. Therefore, we need to consider the seismic anisotropy in cross-hole seismic traveltime tomography. The seismic anisotropic tomography algorithm, which is developed for evaluation of the velocity anisotropy, includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudo-beta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage.

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

Refraction tomography requires an algorithm for efficiently computing the traveltimes and their $Fr\'{e}chet$ derivatives. We have attempted to solve the damped wave field using the frequency domain finite element model ing and then invoked the reciprocity theorem to calculate the $Fr\'{e}chet$ derivative of the traveltime with respect to the subsurface parameter. Then, we used a damped least square method to invert the traveltimes of the Marmousi 2 model. Numerical tests demonstrate that the refraction tomography with large aperture data can be used to estimate the smooth velocity model for the prestack depth migration.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11a
• /
• pp.577-584
• /
• 2000

Seismic refraction survey was peformed for 10 lines along NE-SW and NW-SE directions above Nampoong gallery at Makyo-ri, Dogye, Samcheok, Kangwon-do. 48 geophones were laid in line with the interval of 1m, and a 5Kg hammer was used as a source at 5 points for each line. Data processing was done using reciprocal time method, GRM, and traveltime tomography which utilizes wavefront expansion method for forward process and SIRT for inversion. The result shows that the first layer has its lower boundary between 3.49m and 8.88m. The P-wave velocity of the first and the second layer were estimated as 270 360m/s and 1550 1940m/s respectively. When the boundary of the first and second layer is smooth enough and the velocity difference is large enough, GRM has little advantage over reciprocal time method. The result of reciprocal method and traveltime tomography shows consistency. The northeast part of the boundary has syncline structure, which is similar to the topography above. This implies that the collapse of the cavities of Nampoong gallery result in the subsidence of the ground surface. The subsidence is in progress across the Youngdong railroad, therefore a proper reinforcement work is required.

• Tunnel and Underground Space
• /
• v.18 no.4
• /
• pp.263-271
• /
• 2008

3D seismic tomography technique was assessed for applicability of developed 3D tomography algorithm based on Fresnel volume in the dam-planned area. Reconstructed 3D tomogram based on Fresnel volume and Fast Marching Method(FMM) reveals similar velocity structure to the other geotechnical survey results. With the correlation analysis between RMR data and seismic velocity information, it could provide reliable information of rock mass rate. The applicability of 3D seismic tomography was verified in this study. It would be expected to apply 3D tomography with new developed first arrival calculation and inversion algorithm to the engineering field economically.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.

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