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

• 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.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.5 no.1
• /
• pp.33-45
• /
• 2002

It was necessary to devise new techniques to overcome and enhance the resolution limits of traveltime tomography. Waveform inversion has been one of the methods for giving very high resolution result. High resolution image could be acquired because waveform inversion used not only phase but amplitude. But waveform inversion was much time consuming Job because forward and backward modeling was needed at each iteration step. Velocity-stress method was used for effective modeling. Resolution limits of imaging methods such as travel time inversion, acoustic and elastic waveform inversion were investigated with numerical models. it was investigated that Resolution limit of waveform inversion was similar tn resolution limit of migration derived by Schuster. Horizontal resolution limit could be improved with increased coverage by adding VSP data in cross hole that had insufficient coverage. Also, waveform inversion was applied to realistic models to evaluate applicability and using initial guess of travel time tomograms to reduce non-linearity of waveform inversion showed that the better reconstructed image could be acquired.

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

• Tunnel and Underground Space
• /
• v.11 no.2
• /
• pp.120-121
• /
• 2001

The purpose of this study is to investigate concrete flaw using seismic first arrival and various inversion method. Seismic wave propagation was calculated using finite element method in theoretical modelling and tomogram was made using various inversion methods in theoretical and experimental modelling. Five steps of seismic first arrival were selected from FEM results and these data were used to calculate seismic velocity section. According to the results, exact seismic first arrival picking method was proposed and experimental modelling was conducted.

• The Journal of Engineering Geology
• /
• v.3 no.3
• /
• pp.231-239
• /
• 1993

This paper presents a result on the determination of velocity distribution by a tomographic inversion of crosshole seismic traveltimes in transversely isotropic(aniso tropic) media. The crosshole traveltimes used in this study are synthetic ones computed by ray tracing for some models having isotropic and transversely isotropic velocity distributions. The traveltimes are inverted by a general ART and ansotropic ART which considers the transversely isotropic effect during inver sion. The aniotropic ART gives accurate velodty distributions of transversely isotropic and isotropic models, while the isotropic ART determines accurate velocities only for the isotropic model but inaccurate for the transversely isotropic one. Therefore, the anisotropic ART may be used in case where no information is known on the isotropy or transverse isotropy of a survey area.

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

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.393-398
• /
• 2007

Seismic traveltime tomography technique was conducted at a site composed of black shale. It is well known that black shale has strong anisotropic property. Therefore, the anisotropic property of black shale has to be considered to obtain the appropriate subsurface velocity model by an inversion process. To estimate the anisotropic constant of the velocity of the black shale in the survey area, the relation between the velocity, which is calculated by the straight ray path and the first arrival time, and the angle of the ray propagation was examined. The elliptically shaped relation was found and it reveals that the black shale contains the anisotropic property of velocity. It was also noticed that the horizontal velocity is faster than the vertical velocity. When the estimated anisotropic constant was applied in the process of the velocity inversion for three sets of field data, we could obtain the appropriate velocity structures of the site that is consistent with the result of the geological survey.

• Geophysics and Geophysical Exploration
• /
• v.20 no.1
• /
• pp.1-11
• /
• 2017

In order to improve the accuracy of microseismic epicenter location through the inversion techniques using P and S wave first arrivals, field experiments of microseismic monitoring were performed using borehole 3-component geophones. The direction of epicenter was estimated from the hodograms of P-wave first arrivals through the weight drop experiments in which the $\times$ component of 3-component geophone was aligned to the magnetic north. The picking of S wave first arrival was possible in the polarization filtered data even if S waves are difficult to be identified in raw data. The inversion technique using only P wave first arrival times can often converge to the local minimum when the initial values for epicenter are largely apart from the true epicenter, so that the correct solution can not be found. To solve this problem, the epicenter determination method using differences between P and S wave arrival times was used to estimate proper initial values of epicenter. The inversion result using only P-wave first arrival times which started from the estimated initial values showed the improved accuracy of the epicenter location.