• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.4
• /
• pp.361-371
• /
• 2009

Resistivity tomography is very effective geophysical method to find out the resistivity distribution and its change in time around a tunnel. Thus, the resistivity tomogram can provide helpful information which is necessary for the effective maintenance of the tunnel. However, an air filled tunnel severely distorts tomography data, especially when the current or potential electrode is placed near the tunnel. Moreover, the distortion can often lead to misinterpretation of tomography monitoring data. To solve these problem, we developed a resistivity modeling and time-lapse inversion program which include a tunnel. In this study, using the developed program we assured that the inversion including a tunnel gives much more accurate image around a tunnel, compared with the conventional tomogram where the tunnel is not included. We also confirmed that the time-lapse inversion of resistivity monitoring data defines well resistivity changed areas around a tunnel in time.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.69-72
• /
• 2008

EZTOMO is a crosswell seismic tomography software system. The system has capability of event picking, raytracing, inversion, error analysis, and visualization of the processing results. Waveform of the first arrival signal has been utilized to select the event of the first motion, and uncertainty measured in estimation of the first breaks has been utilized to improve the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.1 no.3
• /
• pp.155-160
• /
• 1998

EM tomography technique has been developed. The algorithm used the extended Born approximations for forward modeling and reconstructed a conductivity image by a smoothness constraint least squares inversion method. Observed data, the vertical components of secondary magnetic fields, were simulated with the 3-D integral equation code. The results showed that the location of anomalous body could be imaged very well, but conductivity of the body was lower than real one and the vertical resolution was much higher than the horizontal resolution.

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

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.

• 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.8 no.2
• /
• pp.129-136
• /
• 2005

The damped least-squares inversion has become a most popular method in finding the solution in geophysical problems. Generally, the least-squares inversion is to minimize the object function which consists of data misfits and model constraints. Although both the data misfit and the model constraint take an important part in the least-squares inversion, most of the studies are concentrated on what kind of model constraint is imposed and how to select an optimum regularization parameter. Despite that each datum is recommended to be weighted according to its uncertainty or error in the data acquisition, the uncertainty is usually not available. Thus, the data weighting matrix is inevitably regarded as the identity matrix in the inversion. We present a new inversion scheme, in which the data weighting matrix is automatically obtained from the analysis of the data resolution matrix and its spread function. This approach, named 'extended active constraint balancing (EACB)', assigns a great weighting on the datum having a high resolution and vice versa. We demonstrate that by applying EACB to a two-dimensional resistivity tomography problem, the EACB approach helps to enhance both the resolution and the stability of the inversion process.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.229-240
• /
• 2007

In this study, seismic anisotropic tomography algorithm was developed for imaging the seismic velocity anisotropy of the subsurface. This algorithm includes several inversion schemes in order to make the inversion process stable and robust. First of all, the set of the inversion parameters is limited to one slowness, two ratios of slowness and one direction of the anisotropy symmetric axis. The ranges of the inversion parameters are localized by the pseudobeta transform to obtain the reasonable inversion results and the inversion constraints are controlled efficiently by ACB(Active Constraint Balancing) method. Especially, the inversion using the Fresnel volume is applied to the anisotropic tomography and it can make the anisotropic tomography more stable than ray tomography as it widens the propagation angle coverage. The algorithm of anisotropic tomography is verified through the numerical experiments. And, it is applied to the real field data measured at limestone region and the results are discussed with the drill log and geological survey data. The anisotropic tomography algorithm will be able to provide the useful tool to evaluate and understand the geological structure of the subsurface more reasonably with the anisotropic characteristics.

• 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.6-17
• /
• 2002

Refraction tomography was developed to interpret subsurface velocity structure easily in topographic conditions. It was applied to synthetic refraction data to find the factors for optimization of applicability of refraction tomography such as configuration of profiling and its length, spacing of geophones and sources and topographic conditions. Also, low velocity layer near VSP hole could be detected by joint inversion with refraction and VSP data. Continuity of subsurface velocity structure in two different spread lines for area of house land development was good in case of applying our algorithm and velocity structure was classified quantitatively to evaluate rippability for engineering works.