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

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

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.101-107
• /
• 2003

3D seismic travel-time tomography algorithm baled on Fresnel volume was developed and its feasibility was investigated by the numerical experiments. To testify the field applicability of the developed algorithm, frequency characteristics and way coverage of the crossholel seismic raw data were investigated and 3D velocity tomogram cube with about 8m spatial resolution was obtained. When compared this 3D velocity cube with the conventional 2D ray tomogram, two results were matched well. We concluded that 3D seismic tomography algorithm developed in this study has enough potential to the field application.

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

3D travel time tomography was conducted to characterize the subsuface structure in the valley area. In this study, an area($200m{\times}200m$), where borehole informations were available to aid in the interpretation, was covered with wide source/receiver coverage. In data acquisition, both hole to hole and reverse VSP array was employed. For the inversion, 3D seismic traveltime tomography algorithm based on Fresnel volume was implemented. When compared 3D velocity cube with the geological survey and drilling logs, both results were matched well. From this, we concluded that 3D seismic travel time tomography has enough potential to the field application.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.85-98
• /
• 2009

The tau-p inversion algorithm is widely employed to generate starting models with most computer programs, which implement refraction tomography. This algorithm emphasises the vertical resolution of many layers, and as a result, it frequently fails to detect even large lateral variations in seismic velocities, such as the decreases which are indicative of shear zones. This study demonstrates the failure of the tau-p inversion algorithm to detect or define a major shear zone which is 50m or 10 stations wide. Furthermore, the majority of refraction tomography programs parameterise the seismic velocities within each layer with vertical velocity gradients. By contrast, the Generalized Reciprocal Method (GRM) inversion algorithms emphasise the lateral resolution of individual layers. This study demonstrates the successful detection and definition of the 50m wide shear zone with the GRM inversion algorithms. The existence of the shear zone is confirmed by a 2D analysis of the head wave amplitudes and by numerous closely spaced orthogonal seismic profiles carried out as part of a later 3D refraction investigation. Furthermore, an analysis of the shot record amplitudes indicates that a reversal in the seismic velocities, rather than vertical velocity gradients, occurs in the weathered layers. The major conclusion reached in this study is that while all seismic refraction operations should aim to provide as accurate depth estimates as is practical, those which emphasise the lateral resolution of individual layers generate more useful results for geotechnical and environmental applications. The advantages of the improved lateral resolution are obtained with 2D traverses in which the structural features can be recognised from the magnitudes of the variations in the seismic velocities. Furthermore, the spatial patterns obtained with 3D investigations facilitate the recognition of structural features such as faults which do not display any intrinsic variation or 'signature' in seismic velocities.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2007.03a
• /
• pp.155-177
• /
• 2007

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.109-115
• /
• 2008

It is very important to estimate the physical properties of survey area and delineate the geological basement in marine site survey for the design of offshore structures. For the purpose of providing high quality data by means of engineering site survey, it is necessary to apply several survey techniques and carry out the integrated interpretation to each other. In this study, we applied single channel seismic reflection method and OBC (Ocean Bottom Cable) type seismic refraction method at shallow marine. We used a dual boomer-single channel streamer as a source-receiver in seismic reflection survey and airgun source-the developed OBC type streamer in seismic refraction survey. We made 24 channels OBC type streamer which has 4m channel interval and each channel is composed of single hydrophone and preamplifier. We tested the field applicability of the proposed method and applied the typical seismic data processing methods to the obtained reflection data in order to enhance the data quality and image resolution. In order to estimate the geological velocity distribution from refraction data, seismic refraction tomography technique was applied. Therefore, we could successfully perform time-depth conversion using the velocity information as an integrated interpretation. The proposed method could provide reliable geologic information such as sediment layer thickness and 3D basement depth map.

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