• The Journal of Engineering Geology
• /
• v.17 no.4
• /
• pp.545-554
• /
• 2007

Landslides is mainly induced by a heavy rainfall, earthquake ground motion, and some other factors like soil mechanics, morphological-geological factors etc. Since the starting point of the failure seemed to be originated at a construction site in the study, it is meaningful to find out the relationship between the landslide and the construction. For this study, the slope failure factor was examined carefully to see that the original natural slope had vulnerability and that the complex ground had unstability changed by construction. A field survey was conducted on the original ground surface and filled-up ground. A laboratory test was also conducted to determine the geomechanical properties of soil samples. 2D and 3D limit equilibrium analysis with changing groundwater level were conducted at the failure depth using a seismic refraction survey. The result shows that the factor of safety is similar stability under all condition, but unstable under saturated condition.

• 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.1 no.1
• /
• pp.8-18
• /
• 1998

In exploration seismology, the Kirchhoff hyperbola has been successfully used to migrate reflection seismo-grams. The mathematical basis of Kirchhoff hyperbola has not been clearly defined and understood for the application of prestack or poststack migration. The travel time from the scatterer in the subsurface to the receivers (exploding reflector model) on the surface can be a kinematic approximation of Green's function when the source is excited at position of the scatterer. If we add the travel time from the source to the scatterer in the subsurface to the travel time of exploding reflector model, we can view this travel time as a kinematic approximation of the partial derivative wavefield with respect to the velocity or the density in the subsurface. The summation of reflection seismogram along the Kirchhoff hyperbola can be evaluated as an inner product between the partial derivative wavefield and the field reflection seismogram. In addition to this kinematic interpretation of Kirchhoff hyperbola, when we extend this concept to shallow refraction seismic data, the stacking of refraction data along the straight line can be interpreted as a measurement of an inner product between the first arrival waveform of the partial derivative wavefield and the field refraction data. We evaluated the Kirchhoff hyperbola and the straight line for stacking the refraction data in terms of the first arrival waveform of the partial derivative wavefield with respect to the velocity or the density in the subsurface. This evaluation provides a firm and solid basis for the conventional Kirchhoff migration and the straight line stacking of the refraction data.

• Geophysics and Geophysical Exploration
• /
• v.24 no.1
• /
• pp.6-19
• /
• 2021

Subsurface physical properties can be obtained and imaged by seismic exploration, and various algorithms have been developed for this purpose. In this regard, root mean square error (RMSE) has been widely used to quantitatively evaluate the accuracy of the developed algorithms. Although RMSE has the advantage of being numerically simple, it has limitations in assessing structural similarity. To supplement this, full-reference image quality assessment (FR-IQA) techniques, which reflect the human visual system, are being investigated. Therefore, we selected six FR-IQA techniques that could evaluate the obtained physical properties. In this paper, we used the full-waveform inversion, because the algorithm can provide the physical properties. The inversion results were applied to the six selected FR-IQA techniques using three benchmark models. Using salt models, it was confirmed that the inversion results were not satisfactory in some aspects, but the value of RMSE decreased. On the other hand, some FR-IQA techniques could definitely improve the evaluation.

• Journal of the Korean Geophysical Society
• /
• v.3 no.1
• /
• pp.57-66
• /
• 2000

The location and geometry of the Ulsan Fault play important roles in interpreting tectonic evolution of the southeastern part of the Korean Peninsula. Dipole-dipole electrical resistivity surveys and seismic refraction surveys were carried out in the Yaksoo area, Ulsan in order to measure the thickness of the alluvium covering the Ulsan Fault and to find associated fracture zones and possibly the location of its major fault plane. The collected data were analyzed and interpreted. Some results reported previously by others were also used in this interpretation. No low resistivity anomalies were found in the cross-sectional resistivity image of the survey line located in the east of the Dong River. In contrast, well-developed continuous low resistivity anomalies were detected in the west of the Dong River. This strongly suggests that the major fault plane of the Ulsan Fault is located under or in the west part of the Dong River. Two refraction boundaries corresponding to the underground water level and the bottom of the alluvium were found by refraction surveys carried out on the limited part of the east survey line. The thickness of the alluvium was found to be about 30 m. Small faults in the basement rock identified by reflection surveys were not detected by both resistivity and refraction seismic surveys. This might be explained by assuming that low resistivity anomaly is more closely related to the clay contents than the water contents. On the other hand, it may be resulted by the limited resolution of the resistivity and refraction surveys. Detailed study is required to clarify the reason. Resistivity survey is frequently considered to be a good exploration method to detect subsurface faults. However, it appears to be less useful than reflection seismic survey in this work. In dipole-dipole resistivity survey, the number of separation should be increased to survey deeper subsurface with the same resolution. However, signal to noise ratio decreases as the number of separation increases. In this survey area, the signal to noise ratio of up to sixteen separations was good enough based on the statistical properties of measurements.

• 한국터널공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.317-331
• /
• 2005

[ ${\circ}$ ] 지반 분포 상환 파악 : 도모그래피 기법 ${\circ}$ 지반의 불연속면 파악 : BHTV, BIPS 및 TSP ${\circ}$ 지반의 횡파 속도 파악, 내진 설계: 표면파 기법, 다운홀 및 크로스홀 ${\circ}$ 시추공의 심도 별 원위치 물성 산출: 밀도 검층, Suspension PS 검충

• Economic and Environmental Geology
• /
• v.30 no.2
• /
• pp.137-142
• /
• 1997

At the tailings dam of the disused Brukunga pyrite mine in South Australia, reaction of groundwater with the tailings causes the formation and discharge of sulphuric acid. There is a need to improve remediation efforts by decreasing groundwater flow through the tailings dam. Geophysical methods have been investigated to determine whether they can be used to characterise variations in depth to watertable and map preferred groundwater flow paths. Three methods were used: transient electromagnetic (TEM) soundings, direct current (DC) soundings and profiling, and self potential (SP) profiling. The profiling methods were used to map the areal extent of a given response, while soundings was used to determine the variation in response with depth. The results of the geophysical surveys show that the voltages measured with SP profiling are small and it is hard to determine any preferred channels of groundwater flow from SP data alone. Results obtained from TEM and DC soundings, show that the DC method is useful for determining layer boundaries at shallow depths (less than about 10 m), while the TEM method can resolve deeper structures. Joint use of TEM and DC data gives a more complete and accurate geoelectric section. The TEM and DC measurements have enabled accurate determination of depth to groundwater. For soundings centred at piezometers, this depth is consistent with the measured watertable level in the corresponding piezometer. A map of the watertable level produced from all the TEM and DC soundings at the site shows that the shallowest level is at a depth of about 1 m, and occurs at the southeast of the site, while the deepest watertable level (about 17 m) occurs at the northwest part of the site. The results indicate that a possible source of groundwater occurs at the southeast area of the dam, and the aquifer thickness varies between 6 and 13 m. A map of the variation of resistivity of the aquifer has also been produced from the TEM and DC data. This map shows that the least resistive (i.e., most conductive) section of the aquifer occurs in the northeast of the site, while the most resistive part of the aquifer occurs in the southeast. These results are interpreted to indicate a source of fresh (resistive) groundwater in the southeast of the site, with a possible further source of conductive groundwater in the northeast.

• Journal of the Korean Geophysical Society
• /
• v.5 no.2
• /
• pp.143-151
• /
• 2002

Bulku temple in the city of Kyungju, Korea, built in 791 and reconstructed in the 20th century, is the home of seven national treasures including two three-story stone pagodas, Dabotap (height 10.4m, width 7.4m, weight 123.2ton) and Seokgatap (height 10.8m, width 4.4m, weight 82.3 ton). An earlier archaeological investigation shows that stone pagodas have experienced severe weathering process which will threaten their stability. At the base part of Dabotap, an offset of the stone alignment is also observed. For the purpose of the structural safety diagnosis of two pagodas, we introduce the nondestructive geophysical methods. Site characteristics around the pagodas are determined by the measurement of multiple properties such as seismic velocity, resistivity, image of GPR(ground-penetrating radar). Near the pagodas, the occurrence of high resistivity (up to 2200 Ωm) is obvious whereas their outskirts have as low as 200 Ωm. For the velocity of the P wave, the site of Dabotap has the range of 500~800 m/s which is higher than counterpart of Seokgatap with the velocity of 300~500m/s, indicating the solider stability of Dabotap site. Consequently, in addition to GPR images, the foundation boundaries beneath each stone pagodas are revealed. The Dabotap site is in the form of an octagon having 6-m-long side with the depth of ~4m, whereas the Seokgatap site the 9m × 10m rectangle with the depth of 3m. These subsurface structures appear to reflect the original foundations constructed against the stone load of ~8 ton/㎡. At the subsurface beneath the northeast of each pagoda, low seismic velocity as well as low resistivity is prominent. It is interpreted to represent the weak underground condition.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.60-68
• /
• 2007

Conventional analysis of spectral induced polarization (SIP) data consists of measuring impedances over a range of frequencies, followed by spectral analysis to estimate spectral parameters. For the quantitative and accurate estimation of subsurface SIP parameter distribution, however, a sophisticated and stable inversion technique is required. In this study, we have developed a two-step inversion approach to obtain the two-dimensional distribution of SIP parameters. In the first inversion step, all the SIP data measured over a range of frequencies are simultaneously inverted, adopting cross regularisation of model complex resistivities at each frequency. The cross regularisation makes it possible to enhance the noise characteristics of the inversion by imposing a strong assumption, that complex resistivities should show similar characteristics over a range of frequencies. In numerical experiments, we could verify that our inversion approach successfully reduced inversion artefacts. As a second step, we have also developed an inversion algorithm to obtain SIP parameters based on the Cole-Cole model, in which frequency-dependent complex resistivities from the first step are inverted to obtain a two-dimensional distribution of SIP parameters. In numerical tests, the SIP parameter images showed a fairly good match with the exact model, which suggests that SIP imaging can provide a very useful subsurface image to complement resistivity.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.136-147
• /
• 2004

In this study, the characteristics of magnetotelluric (MT) responses due to a three-dimensional (3-D) body are analyzed with 3-D numerical modeling. The first model for the analysis consists of a single isolated conductive body embedded in a resistive homogeneous half-space. The second model has an additional conductive overburden while the other conditions remain the same as the first one. The analysis of apparent resistivities shows well that the 3-D effects are dominant over some frequency range for the first model. Two mechanisms, current channeling and induction, for secondary electric fields due to the conductive body are analyzed at various frequencies: at high frequencies induction is more dominant than channeling, while at low frequencies channeling is more dominant than induction. Tippers have a strong relation to the position of anomalous body and the real and imaginary parts of induction vector also indicate the position of anomalous body. off-line conductive anomaly sometimes causes severe problem in 2-D interpretation. In such case, induction vector analysis can give information on the existence and location of the anomalous body. Each parameter of the second model shows similar responses as those of the first model. The only difference is that the magnitude of all parameters is decreased and that the domain showing the 3-D effects becomes narrower. As shown in this study, the analysis of 3-D effects provides a useful and effective means to understand the 3-D subsurface structure and to interpret MT survey data.