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

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

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.233-238
• /
• 2009

An inversion method regulated by the error in the measurement of the first arrival time was developed, and we conducted a feasibility study by applying the method to a real cross-well seismic data. The inversion is a two-step regulation process; 1) derive the measurement error bound based on the resolution of the velocity image want to derive, and exclude the records whose picking error is larger than the error bound, 2) set the travel time residual to zero if the residual is less than the measurement error. This process prevents the trivial residuals are accumulated and contribute to the velocity-model update. Comparison of two velocity images, one by using all records and another by using the regulate inversion method, shows that the later velocity image exhibits less numerical artefacts, and it also indicates that, according to the Fermat's principle, the latter image is a more feasible velocity model.

• 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.5 no.3
• /
• pp.150-156
• /
• 2002

In this study, a software for crosswell seismic tomography is developed. The software consists of first arrival picking and adjusting module, crosswell traveltime tomography module, and imaging module. This software allows saying the picked first arrival times into the header of seismic data, and using this data directly to the input of crosswell seismic tomography. With an imaging module, velocity structures and ray path can be imaged directly from the output of the tomography module. Because it is developed on the basis of the SU under the Linux and the GUI environment for user, this software can be carried out directly the first arrival picking, inversion and tomogram for crosswell tomography data in the field. Therefore, this software can be improved the applicability of site investigation by tomography method.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.167-173
• /
• 1999

Since the resolution of the 2-D hole-to-hole seismic traveltime tomography is affected by the limited ray transmission angle, various methods were used to improve the resolution. Linear traveltime interpolation(LTI) ray tracing method was chosen for forward-modeling method. Inversion results using the LTI method were compared with those using the other ray tracing methods. As an inversion algorithm, SIRT method was used. In the iterative non-linear inversion method, the cost of ray tracing is quite expensive. To reduce the cost, each raypath was stored and the inversion was performed from this information. Using the proposed method, fast convergence was achieved. Inversion results are likely to be affected by the initial velocity guess, especially when the ray transmission angle was limited. To provide a good initial guess for the inversion, generalized regression neural network(GRNN) method was used. When the transmitted raypath angle is not limited or the geological model is very complex, the inversion results are not affected by initial velocity model very much. Since the raypath angles, however, are limited in most geophysical tomographic problems, the enhancement of resolution in tomography can be achieved by providing a proper initial velocity model by another inversion algorithm such as GRNN.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.330-335
• /
• 2010

In order to estimate the confidence level of the velocity distribution shown in a velocity image reconstructed from a travel-time tomography, the ray coverage and the inversion characteristics of the system matrix were investigated. The targets of the analysis is the first arrival travel-time, the raypath information, and the resulting velocity model. The ray coverage, degree of ray and model coupling, was estimated by the number of rays and total ray length in a velocity grid, and information regarding the resolution and uncertainties involved in the reconstructed velocity model was derived from the results of the SVD analysis of the system matrix that relates the data space (first arrival travel times) to the model space (velocity distribution in tomogram).

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.157-163
• /
• 2004

In this study, a travel-time tomography software was developed under the MS Windows system and GUI environment for user. The software supports following features: (1) supporting various data input format (2) flexible treatment of shot and receiver coordinate coding (3) flexible first arrival picking and modification (4) easy modification of intermediate tomogram. It is expected that the effort of the user can be minimized in each data processing step.

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