• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.83-86
• /
• 2006

From numerical modeling test $-{\rho}a$ by one steel casing borehole near resistivity survey line can be acquired. Negative apparent resistivities even in the flat area are surely subsurface information. Inversion technique for those need to be developed in the near future.

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

Due to the recent development of instruments which enabled the measurements of subtle IP effect in the ground and analysis algorithms, complex resistivity (CR) method was expanding its application to various field. In this study, we applied the CR method to the test site where the ground reinforcement had been done by injecting the cement mortar for investigating the effect of ground reinforcement. For this site, resistivity monitoring and tomography survey was carried out while the ground reinforcement had been made by the grouting. From the result, the anomalous region that was shown on the result of resistivity 4D monitoring was coincident with those of phase section in the CR method, because the cement grouting material had the strong IP effects. It might be expected that the CR method should be very powerful surveying tool for the similar purpose.

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

We have measured resistivities for undisturbed soil samples collected from dredged bank with sea sand, and analyzed with one-dimensional inversion results from small-loop electromagnetic survey data. From the relationship between the two resistivities, it appeared that calculated resistivities were remarkably consistent with measured resistivities. Correlation relationships between resistivity values and cone resistance were analyzed after comparing inversion results with cone resistance. It turns out that the region with below 1 ohm-m is correspondent to that of with less than $50\;kgf/cm^2$ in dredged bank with sea sand. From the study result, resistivity monitoring of small-loop EM periodically is proved to be more effective to maintain the stability of embankment dike.

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

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2003.09a
• /
• pp.112-115
• /
• 2003

쓰레기매립장의 침출수 누출 탐지를 위한 물리탐사법은 침출수의 전기전도도를 대상으로 한 전기탐사법이 주로 적용되어 왔다. 이 연구는 제주도에 위치한 비위생매립지를 대상으로 주변 수리지질환경 파악을 위한 전기비저항 탐사와 더불어 매립장 경계부와 토양층을 통한 침출수 누출 영역을 효율적으로 규명하기 위하여 다중주파수를 이용한 소형루프 전자탐사를 병행하여 수행하였으며, 아울러 오염지역의 천부 지하수 유동 방향을 추정하기 위하여 침출수의 유동에 의해 발생될 수 있는 자연전위를 모니터링 하였다. 전기비저항 탐사 및 자연전위 모니터링 결과 매립지 주변을 통한 침출수의 누출 가능성은 거의 없는 것으로 나타났으며, 이는 주변 하천에서 정기적으로 실시한 수질 분석 결과와 일치된다. 또한 소형루프 전자탐사 자료에 대하여 공간 필터링 및 1차원 역산법을 적용한 결과 매립장 경계부의 위치 및 매립장 내부의 심도별 매립물에 의한 저비저항 이상대를 효과적으로 도출 할 수 있었다.

• Journal of the Korean Geophysical Society
• /
• v.2 no.3
• /
• pp.217-224
• /
• 1999

Although faults can be found by geological surveys, the surface traces of faults are not easily discovered by traditional geological surveys due to alluvia. In and around faults and fracture zones, the electrical resistivity appears to be lower than that of the surroundings due to the content of groundwater and clay minerals. Therefore, electrical resistivity surveys are effective to search buried faults and fracture zones. The dipole-dipole array electrical resistivity surveys, which could show the two dimensional subsurface electrical resistivity structure, were carried out in two areas, Yongdang-ri, Woongsang-eup, Yangsan-si, Kyungsangnam-do and Malbang-ri, Woedong-eup, Kyungju-si, Kyungsangpook-do. The one was next to the Dongrae Fault and the other near the Ulsan Fault was close to the region in which debatable quaternary fault traces had been found recently. From each measured data set, the electrical resistivity cross-section was obtained using the inversion program the reliability of which was analyzed using analytic solutions. A low resistivity zone was found in the inverted cross-section from the Yongdang-ri area survey data, and two low resistivity zones were found in that from the Malbang-ri area survey data. They were almost vertical and were 15∼20 m wide. Accounting the shape and the very low resistivity values of those zones (<100 Ωm)in the inverted section, they were interpreted as fracture zones although they should be proven by trenching. The reliability of the interpretation might be improved by adding some more parallel resistivity survey lines and interpreting the results in 3 and/or adding other geophysical survey.

• Journal of the Korean earth science society
• /
• v.29 no.6
• /
• pp.447-453
• /
• 2008

Recently, the electrical DC survey has frequently been performed to assess the safety of embankment. This study showed that the damaged section of embankment could be appropriately detected by the survey only when the three dimensional effect was correctly considered. The shape of the three dimensional embankment was numerically implemented, and a proper modeling was performed to confirm the effect by analyzing the apparent and inverted true resistivity. Then, the field work was carried out. The three dimensional interpretation distinguished the erroneous weak zones from the geometrical artifact, and the embankment was ensured as safe both by the additional survey performed in rainy season and the partial excavation for direct observation.

• Proceedings of the KIEE Conference
• /
• 1997.07a
• /
• pp.3-5
• /
• 1997

지구 물리학이나 의공학 분야등에서 이용되왔던 전기비저항 탐사법은 관심 영역에 전류 입력을 가한 후, 그에 대한 전압 응답을 측정하여 관심 영역 내의 전기비저항 분포를 규명하는 방법으로서 역해석 문제의 범주에 포함된다. 따라서 일반적인 역해석 문제가 지니고 있는 해의 존재성, 유일성, 그리고 측정 데이터에 대한 해의 연속적 의존성이라는 기본적 문제들을 가지게된다. 이러한 역해석 문제의 해결에는 정확한 정해석 풀이법과 효율적인 역해석 방법이 요구되어진다. 본 논문에서는 정해석 방법으로 유한요소법을, 역해석 방법으로는 전체 최적점을 발견할 가능성이 높은 유전 알고리즘을 최적화 방법으로 사용하였다. 기존의 역해석 문제의 해결책으로 제시되어왔던 기울기 방법에 기반한 결정론적 최적화 알고리즘들이 지니고 있는 국소해로의 수렴, 즉 단순한 전기비저항 분포의 불연속성 확인이라는 한정된 정보의 획득을 넘어서 실제 전기비저항 분포와 가장 가까운 분포는 전체 최적점 근처에서 발견될 수 있음을 보이고자 한다. 이러한 전기비저항 분포의 역해석적인 규명을 간단한 2차원 수치해석문제를 풀어보므로서 확인해본다.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.262-271
• /
• 2002

This study presents a practical procedure for the Bayesian inversion of geophysical data. We have applied geostatistical techniques for the acquisition of prior model information, then the Markov Chain Monte Carlo (MCMC) method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter.