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

• Proceedings of the KSEEG Conference
• /
• 2001.04a
• /
• pp.56-58
• /
• 2001

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.09a
• /
• pp.21-36
• /
• 2005

Although the main geology of Korea consists of granite and gneiss, it Is not uncommon to encounter anisotropy Phenomena in crosshole radar tomography even when the basement is crystalline rock. To solve the anisotropy Problem, we have developed and continuously upgraded an anisotropic inversion algorithm assuming a heterogeneous elliptic anisotropy to reconstruct three kinds of tomograms: tomograms of maximum and minimum velocities, and of the direction of the symmetry axis. In this paper, we discuss the developed algorithm and introduce some case histories on the application of anisotropic radar tomography in Korea. The first two case histories were conducted for the construction of infrastructure, and their main objective was to locate cavities in limestone. The last two were performed In a granite and gneiss area. The anisotropy in the granite area was caused by fine fissures aligned in the same direction, while that in the gneiss and limestone area by the alignment of the constituent minerals. Through these case histories we showed that the anisotropic characteristic itself gives us additional important information for understanding the internal status of basement rock. In particular, the anisotropy ratio defined by the normalized difference between maximum and minimum velocities as well as the direction of maximum velocity are helpful to interpret the borehole radar tomogram.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.7-17
• /
• 2011

We have developed an inversion code for three-dimensional (3D) resistivity tomography including the anisotropy effect. The algorithm is based on the finite element approximations for the forward modelling and Active Constraint Balancing method is adopted to enhance the resolving power of the smoothness constraint least-squares inversion. Using numerical experiments, we have shown that anisotropic inversion is viable to get an accurate image of the subsurface when the subsurface shows strong electrical anisotropy. Moreover, anisotropy can be used as additional information in the interpretation of subsurface. This algorithm was also applied to the field dataset acquired in the abandoned old mine area, where a high-rise apartment block has been built up over a mining tunnel. The main purpose of the investigation was to evaluate the safety analysis of the building due to old mining activities. Strong electrical anisotropy has been observed and it was proven to be caused by geological setting of the site. To handle the anisotropy problem, field data were inverted by a 3D anisotropic tomography algorithm and we could obtain 3D subsurface images, which matches well with geology mapping observations. The inversion results have been used to provide the subsurface model for the safety analysis in rock engineering and we could assure the residents that the apartment has no problem in its safety after the completion of investigation works.

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

• The Journal of Engineering Geology
• /
• v.25 no.2
• /
• pp.251-263
• /
• 2015

Geological mapping, borehole drilling, electrical resistivity, and seismic tomography surveys were conducted in order to map underground cavities and better understand the mechanisms driving subsidence in a limestone region in Korea. Limestone outcrops in the study area generally alternate between calcite-rich and calcite-poor rock. The results reveal that in areas experiencing subsidence, cavities occur mainly around soil-rock boundaries at depths of 7~14 m. These results are based on comparative analyses of electrical resistivity, seismic tomography, and borehole logging data. The volumes of the cavities are relatively small in a range of 558~835 ㎥ and they have a shape typical of suffosion sinkholes, which are typically found where sandy soils overlie bedrock cavities.

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

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2004.08a
• /
• pp.138-156
• /
• 2004

GPR(Ground Penetrating Radar) was used for imaging the interior of the historical masonry such as stone pagoda in order to provide the basic information of safely inspection. The scope of the imaging was restricted to the foundation part of stone pagoda that transferred the load of the pagoda to the ground. Kirchhoff migration and traveltime tomography was used for imaging the outer stone and the inside of stone pagoda, respectively. From the migrated images, we could measure the thickness and the shape of the boundaries of the outer stone in the foundation part. From the reconstructed tomograms for the physical model, we could get the GPR propagation velocity distribution and exactly find the position of the air in the model and calculate the average velocity with respect to the different filling materials. The properties and the shape of the interior materials of stone pagoda can be basic informations for the safety inspection.