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

• Journal of the Korean Geophysical Society
• /
• v.8 no.2
• /
• pp.67-74
• /
• 2005

Two-dimensional S-wave velocity sections from SH-wave refraction tomography and surface wave dispersions were obtained by inverting traveltimes of first arrivals and surface wave dispersions, respectively. For the purpose of comparison, a P-wave velocity tomogram was also obtained from a P-wave refraction profiling. P and Rayleigh waves generated by vertical blows on a plate with a sledgehammer were received by 100- and 4.5-Hz geophones, respectively. SH-waves generated by horizontal blows on both sides of a 50 kg timber were received by 8 Hz horizontal geophones. The shear-wave signals were enhanced subtracting data of left-side blows from ones of the right-side blows. Shear-wave velocities from tomography inversion of first-arrival times were compared with ones from inverting dispersion curves of Rayleigh waves. Although the two velocity sections look similar to each other in general, the one from the surface waves tends to have lower velocities. First arrival picking of SH waves is troublesome since P and PS-converted waves arrive earlier than SH waves. Application of the surface wave method, on the other hand, is limited where lateral variation of subsurface tructures is not mild.

• Journal of the Korean earth science society
• /
• v.27 no.4
• /
• pp.425-432
• /
• 2006

A seismic tomography using refraction waves is applied to provide information on depth of basement rocks and leachate distribution of the Noeun waste landfill site for the stage of preliminary environmental survey. This method is generally applied to civil and environmental areas. Three lines, apparently perpendicular to the potential leachate flow direction in this site, were installed to investigate the waste landfill site in pseudo three dimensional geometry. The results show that the site is composed of 3 layers and depth of basement becomes shallower at the upstream area of the landfill site than that of the downstream area. Moreover, some parts of the second layer and the basement at the down stream area are partially infiltrated by the leachate, probably related to the disturbed distribution of the different velocity materials within the second layer. In Conclusion, refraction wave tomography is found to be one of the most efficient way to investigate waste landfill site.

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

• Proceedings of the KSR Conference
• /
• 2005.05a
• /
• pp.610-614
• /
• 2005

Rock mass classification by using electrical resistivity tomography(ERT) method is widely performed for the determination of rock support type in tunnel design. In the region of high electrical noise level, however, the result of the ERT will have many erroneous features. In this study, the back ground electrical noise had been measured to find out the reason why the results of ERT in this area did not agree to the expected geology confirmed by boreholes. In order to overcome this limitation of ERT, a hybrid surface-borehole array seismic refraction tomography had been followed. Using this technique, we could get P-wave velocity section including the depth level of tunnel. The comparison of the P-wave velocity and RMR shows fairly good statistical relationship to make it possible to set up the rock mass classification for the entire tunnel line.

• Journal of the Korean Institute of Telematics and Electronics B
• /
• v.28B no.5
• /
• pp.335-342
• /
• 1991

A real-tme nonlinear parameter tomography is pumping wave method. This tomorgraphy has a merit which requires no 180$^{\circ}$ projection datum, while the ray-bending effect is remrkably remained on the reconstructed image. In this paper we intend to compensate this ray-bending effect using the perturbation method. Impulsive pumping wave makes derived compensative term simple form, nad the compensative image is easily obtained. We perform computer simulation to confirm the improvement of corrected imate.

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

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.19-24
• /
• 1998

This paper deals with tomographic travel time inversion of three dimensional seismic refraction survey conducted over a dipping interface. The slowness, and thus velocity as its reciprocal, distribution on the subsurface interface is to be determined applying an ART with under-relaxtion parameter. The models chosen are realistic, i.e., most likely to be met in engineering seismics, and the interface includes anomalous zones. It is found that, generally speaking, the inversion could be misleading or meaningless without the correction of the dip of the interface. This is rather surprising when we recall that usual assumption for the interpretation of refraction seismics data is the horizontal attitude of structures within the limit of $15^{\circ}$ dip or so. To make the present method tenable for a new means of routine seismics, some practical ways of identifying head wave arrivals are to be devised.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.05a
• /
• pp.287-291
• /
• 2005

Two-dimensional velocity tomograms of P- and S-waves were obtained by inverting traveltimes of first arrivals. The two sections of shear-wave velocity show similar features as a whole, with smaller values on the section from surface wave dispersions. Difficulties in picking SH-wave phases due to noise and later arrivals than P waves and PS converted waves are experienced. In addition, a flat layer model based on the surface wave inversion prohibits applications of the method where sgear wave velocities vary strongly in the lateral direction.