• Economic and Environmental Geology
• /
• v.29 no.3
• /
• pp.345-355
• /
• 1996

The effectiveness of various electrode configurations in horizontal mappings and 1-D inversions of vertical sounding data for delineating 2-D structures was studied. Apparent resistivity values of three point, dipole-dipole, Wenner, and Schlumberger mappings were simulated for such structures as vertical dyke, tabular prism, buried vertical fault, ramp and complex structure by finite difference method (FDM) and they were compared with each other. Also 2-D cross sections for three structures obtained by interpolation of 1-D inverted sounding data in terms of three layers were compared for Schlumberger and Wenner arrays. On these cross sections, horizontal and vertical resistivity interfaces of the 2-D structures are revealed relatively clearly. Apparent resistivity curves of Schlumberger mapping show vertical resistivity discontinuities very well. On the whole, Schlumberger array is superior to the other arrays in electric sounding as well as mapping. This study clearly indicates that interpretations of 2-D structures based on 1-D inversion are possible.

• The Journal of Engineering Geology
• /
• v.7 no.2
• /
• pp.113-126
• /
• 1997

Electrical resistivity methods of dipole - dipole array profiling and Schiumberger array sounding were tested on a segment of the Woraksan granitic batholith for the research into the imaging of irregular attitudes of fracture zones in the crystaaline rock in terms of processing and interpretation schemes. By the dipole - dipole array method, inhomogeneities such as small scale of fracture zones were properly delineated down at some depth even within hard rock environment. Fracture zones were interpreted to be at the boundaries between the high amplitude zone and very low amplitude zone in the resistivity plot and they were also successfully outlined in two - dimensional layer and pseudo - three - dimensional volume constructed by the incorporation of vertical sounding data. The surface location of the fracture zones was correlated by the zero - crossing point in the VLF(very low frequency) electromagnetic data. Pseudo - three - dimensional attitudes of fracture zones were efficiently illuminated by optimum projection angle. The mean of bulk resistivity for the Woraksan granite and the near fracture zones is estimated to be approximately of 4,000 ohm - m which is much higher than the value of 700 ohm - m for the Rwachunri limesilicate environment. This difference is due to both the rock type, i.e., biotite granite vs limesilicate, and the occurrence of secondary openings of fold and fault associated with the intrusion of granite. In this study statistical analyses on the resistivity color plot were performed in terms of three representative statistical moments, i.e., standard deviation, skewness, and kurtosis. The fracture zones in the standard deviation plot were characterized by the higher value, compared to the value of homogeneous portion. The upper boundary of the high resistivity zone was also successfully delineated in the skewness and kurtosis plots.

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

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.104-106
• /
• 2002

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

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

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.89-97
• /
• 2007

A three-dimensional (3D) magnetotelluric (MT) survey has been carried out to delineate subsurface structures and possible fractures, for development of low-temperature geothermal resources in Pohang, Korea. Quite good quality MT data could be obtained throughout the survey region by locating the remote reference in Kyushu, Japan, which is ${\sim}480\;km$ from the centre of the field site. 3D modelling and inversion are performed taking into account the sea effect in MT measurements near the seashore. The nearby sea in the Pohang area affects MT data at frequencies below $1\;Hz{\sim}0.2\;Hz$, depending on the distance from the seashore. The most severe sea effects were observed in the south-east parts of the survey area, closer to Youngil Bay. 3D inversion with and without the seawater constraint showed very similar results at shallow depths, roughly down to 2 km. At greater depths, however, a strong sea effect seems to form a fictitious conductive structure in ordinary 3D inversion, especially in the south-eastern part of the survey region. Comparison between drilling results and the resistivity profiles from inversions showed that five layered structures can be distinguished the subsurface beneath the target area. They are: (a) semi-consolidated mudstones with resistivity less than $10\;{\Omega}m$, which are ${\sim}300\;m$ thick in the northern part and ${\sim}600\;m$ thick in the southern part of the survey area; (b) occasional occurrence of trachybasalt and lapilli tuff within the mudstone layer has resistivity of a few tens of${\Omega}m$, (c) intrusive rhyolite ${\sim}400\;m$ thick has resistivity of several hundreds of ${\Omega}m$, (d) alternating sandstone and mudstone down to 1.5 km depth shows resistivity of ${\sim}100\;{\Omega}m$, (e) a conductive structure was found at a depth of ${\sim}3\;km$, but more geological and geophysical study should be carried out to identify this structure.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.209-215
• /
• 2014

In order to develop a new mineral exploration technique, a study was carried out about the potential mapping of Moisan area using SIP (Spectral Induced Polarization) data. The SIP inversion results were classified according to the geological regions, and the distribution characteristics of resistivity and phase values of SIP data were analyzed at the ore region. Based on the characteristics of SIP of ore bodies, we performed 3D potential mapping of Moisan area. The analyzed potential map was verified using that the locations and patterns of high potential regions of the results are well matched with those of the known ore bodies. If we get the higher spatial resolution SIP data, the potential mapping technique using SIP data can be effectively applied to the estimation of mining deposit.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.33-40
• /
• 2009

In field surveys using the dipole-dipole electrical resistivity method, we often encounter negative apparent resistivity. The term 'negative apparent resistivity' refers to apparent resistivity values with the opposite sign to surrounding data in a pseudosection. Because these negative apparent resistivity values have been regarded as measurement errors, we have discarded the negative apparent resistivity data. Some people have even used negative apparent resistivity data in an inversion process, by taking absolute values of the data. Our field experiments lead us to believe that the main cause for negative apparent resistivity is neither measurement errors nor the influence of self potentials. Furthermore, we also believe that it is not caused by the effects of induced polarization. One possible cause for negative apparent resistivity is the subsurface geological structure. In this study, we provide some numerical examples showing that negative apparent resistivity can arise from geological structures. In numerical examples, we simulate field data using a 3D numerical modelling algorithm, and then extract 2D sections. Our numerical experiments demonstrate that the negative apparent resistivity can be caused by geological structures modelled by U-shaped and crescent-shaped conductive models. Negative apparent resistivity usually occurs when potentials increase with distance from the current electrodes. By plotting the voltage-electrode position curves, we could confirm that when the voltage curves intersect each other, negative apparent resistivity appears. These numerical examples suggest that when we observe negative apparent resistivity in field surveys, we should consider the possibility that the negative apparent resistivity has been caused by geological structure.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2011.10a
• /
• pp.95-96
• /
• 2011