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

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.163-166
• /
• 2002

Initial velocity model close to real velocity structure of the subsurface governs the quality of image of prestack depth migration. Geophysicists employ velocity estimation tools such as velocity analysis (curvature method, coherency inversion), tomography and waveform inversion. We present a reflection tomography that parameterizes the subsurface into the movable blocks. By interpreting the depth-migrated section or stacked section, we can design an initial constant velocity model having only impedance boundaries. We use shooting-raytracing method that allows us to calculate the Jacobian-matrix efficiently.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.314-321
• /
• 2007

Recently, multi-dimensional joint inversion of geophysical data based on fundamentally different physical properties is being actively studied. Joint inversion can provide a way to obtaining much more accurate image of the subsurface structure. Through the joint inversion, furthermore, it is possible to directly estimate non-geophysical material properties from geophysical measurements. In this study, we developed a new algorithm for jointly inverting dc resistivity and seismic traveltime data based on the multiple constraints: (1) structural similarity based on cross-gradient, (2) correlation between two different material properties, and (3) a priori information on the material property distribution. Through the numerical experiments of surface dc resistivity and seismic refraction surveys, the performance of the proposed algorithm was demonstrated and the effects of different regularizations were analyzed. In particular, we showed that the hidden layer problem in the seismic refraction method due to an inter-bedded low velocity layer can be solved by the joint inversion when appropriate constraints are applied.

• Journal of the Korean Geophysical Society
• /
• v.4 no.1
• /
• pp.21-33
• /
• 2001

Simultaneous inversion of first-arrivals of local earthquakes recorded by the Korea Meteorological Administration(KMA) seismograph network from 1991 to 1998 is made to derive 1D crustal velocity structure of the Korean peninsula. Twenty-nine events with 178 observations are used in the inversion. Average crustal P-wave velocity turns out to be about 6.3 km/sec, and crustal thickness and upper mantle P-wave velocity are estimated as 33 km and 7.9 km/sec, respectively. Results of inversion indicate the possibility of the low velocity layer in the lower crust. Joint inversion is applied to estimate hypocenters, station delays, and velocities simultaneously. Relative station corrections for 11 stations range from zero to about 1.2 sec. Analysis of the synthetic data shows that estimates of hypocenter locations and station corrections as well as averaged crustal structure are reliable for the given data set..

• Journal of the Korean earth science society
• /
• v.22 no.1
• /
• pp.20-29
• /
• 2001

P wave traveltime delays has been measured and a preliminary tomouaphic inversion has been performed using the data collected from the network deployed by Korea National University of Education in the southcentral part Korea. The maximum variation in relative traveltime residuals is almost 0.7 seconds. A large azimuthal variation in traveltime residuals is observed at the stations in the northwestern part of the study area. This might indicate the existence of lateral velocity heterogeneity beneath the study area. Although, the inverted tomographic image seems to be correlated with the tectonic boundary between Kyunggi Massif and Okcheon Belt, the confidence level is presently low. We should place on the importance of this study for finding preferable inversion parameters and predicting probable result. Better tomographic image of the study area can be obtained in the near future when the data are accumulated.

• Geophysics and Geophysical Exploration
• /
• v.6 no.2
• /
• pp.95-100
• /
• 2003

In this study, we developed reflection tomography inversion algorithm using Straight Ray Technique (SRT) which can calculate travel time easily and fast for complex geological structure. The inversion process begins by setting the initial velocity model as a constant velocity model that hat only impedance boundaries. The inversion process searches a layer-interface structure model that is able to explain the given data satisfactorily by inverting to minimize data misfit. For getting optimal solution, we used Gauss-Newton method that needed constructing the approximate Hessian matrix. We also applied the Marquart-Levenberg regularization method to this inversion process to prevent solution diverging. The ability of the method to resolve typical target structures was tested in a synthetic salt dome inversion. Using the inverted velocity model, we obtained the migration image close to that of the true velocity model.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.140-145
• /
• 2011

In order to analyze the influence of the noise on a cross-well traveltime tomography to detect a small scale low velocity body in a homogeneous medium, the first arrival travel times were computed one a tunnel model by a finite-difference ray tracing scheme. Three different types and four different intensity levels of white noises were added to the computed first arrival travel times, and velocity tomograms were constructed using an iterative inversion method (SIRT). Tomograms with the noise intensity up to 10% of the maximum traveltime delay in the tunnel model, showed the exact location of the tunnel. However, the velocity shown at the tunnel location was not close to air velocity but only slightly less than the velocity of the background medium. The additive random noise showed significantly less degree of influence on the resulting tomogram than the source- and receiver consistent noise.

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

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.49-56
• /
• 1998

Seismic tomography has been widely used as high resolution subsurface imaging techniques in engineering applications. Although most of the techniques have been using travel time inversion, waveform method is being driven forward owing to the progress of computational environments. Although full-waveform inversion method has been known as the best method in terms of model resolving power without high-frequency restriction and weak scattering approximation, it has practical disadvantage that it is apt to get stuck in local minimum if the initial guess is far from the actual model and it consumes so much time to calculate. In this study, 2-D full-waveform inversion algorithm in acoustic medium is developed, which uses result of traveltime tomography as initial model. From the application on synthetic data, it is proved that this approach can efficiently reduce the problem of conventional approaches: our algorithm shows much faster convergence rate and improvement of model resolution. Result of application on physical modeling data also shows much improvement. It is expected that this algorithm can be applicable to real data.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.142-152
• /
• 2008

SPT-Uphole method was introduced for the evaluation of near subsurface shear wave velocity (Vs) profile. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in geotechnical site investigation was used as a source and several surface geophones in line were used as receivers. 1D shearwave velocity profile can be obtained in the manner of downhole method, Vs distribution map which is the triangular shape around the boring point can be developed by tomography inversion. To obtain the exact travel time information of shear wave component, a procedure using the magnitude summation of vertical and horizontal components was used based on the evaluation of particle motion at the surface. It was verified that proposed method could give reliable Vs distribution map through the numerical study using the FEM (Finite Element Method) model. Finally, SPT-Uphole method was performed and the feasibility of proposed method was verified in the field.