• 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.9 no.3
• /
• pp.241-249
• /
• 2006

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 commonly used, however, may not give sufficiently precise computational results of traveltime data 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. Considering the complex geology of Korea, 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 performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• 한국지구물리탐사학회:학술대회논문집
• /
• 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 KSEEG Conference
• /
• 2001.04a
• /
• pp.115-118
• /
• 2001

Suppressed Wave equation Estimation of Traveltime (SWEET, Shin et, al., 2000) 알고리즘을 이용하여 등방성 및 이방성 탄성 매질에 대한 ched 주시와 진폭을 계산하는 방법을 개발하였다. SWEET 알고리즘을 2차원 등방성 및 이방성 탄성 매질의 속도 구조에 적용하여 P파의 초동 주시와 진폭을 계산 할 수 있었다. 본 논문에서는 간단한 등방성 균질 탄성 모형, 복잡한 등방성 탄성 모형 및 간단한 이방성 균질 모형에 대한 수치 계산 결과를 보여 줄 것이다.

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

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

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

The classic graphical method to determine the epicenter uses differences between the arrival times of P and S waves at each station. In this research, a robust approach is proposed, which provides a fast and intuitive estimation of earthquake epicenters. This method uses an empirical relationship between epicentral distance and traveltime of the first arrival P phase of local or regional earthquake. The relationship enables us to estimate epicentral distances and draw epicentral circles from each station with P-traveltimes counted from a probable origin time. As the assigned time is getting close to the origin time of the earthquake, epicentral circles begin to intersect each other at a possible location of the epicenter. Then the possibility of the epicenter can be expressed by a function of the time and the space. We choose the location which gives the minimum standard deviation of the origin time as an estimated epicenter. In this research, 918 P arrival times from 84 events occurring from 2005 to 2006 listed in the KMA earthquake catalog are used to determine the empirical P-traveltime function of epicentral distances.

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

• Proceedings of the Korean Geotechical Society Conference
• /
• 2002.03a
• /
• pp.695-702
• /
• 2002

Groundwater flow is primarily influenced by the presence of fractures, functioning as conduits. To block the flow, grouting operation is commonly used. Thereby the fractures are then expected to be sealed, which will add to enhance the shear strength in rock. This far, regarding the assessment of grouting efficiency, however, there's been a considerable uncertainty That is, several geophysical methods of high resolution such as tomography, S-wave logging have produced a significant amount of measurable response caused by grouting, but they can inevitably be used only for the qualitative assessment. Thus, this paper deals with an accurate quantitative assessment about the grouting result. In this, a new strategy is introduced, based mainly on evaluating the opening of fractures. For fracture-opening investigation purposes, borehole Televiewer has already proven to be an excellent logging technique that produces both amplitude image and traveltime image. As well known, the traveltime image can be converted to a high precision 3D caliper log with max. 288 arms, which allows to observe the opening of fractures. To evaluate the fracture opening from the traveltime image, an algorithm of practical use was developed, in which image correction due to the borehole deviation, feature discrimination of wall roughness from fractures, automatic evaluation procedure etc. were considered. Field examples are shown to confirm the efficiency of the suggested method.

• Proceedings of the KSEEG Conference
• /
• 2001.04a
• /
• pp.107-110
• /
• 2001