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

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

In the geophysical monitoring to understand the change of subsurface material properties with time, the time-invariant static subsurface model is commonly adopted to reconstruct a time-lapse image. This assumption of static model, however, can be invalid particularly when fluid migrates very quickly in highly permeable medium in the brine injection experiment. In such case, the resultant subsurface images may be severely distorted. In order to alleviate this problem, we develop a new least-squares inversion algorithm under the assumption that the subsurface model will change continuously in time. Instead of sampling a time-space model into numerous space models with a regular time interval, a few reference models in space domain at different times pre-selected are used to describe the subsurface structure continuously changing in time; the material property at a certain space coordinate are assumed to change linearly in time. Consequently, finding a space-time model can be simplified into obtaining several reference space models. In order to stabilize iterative inversion and to calculate meaningful subsurface images varying with time, the regularization along time axis is introduced assuming that the subsurface model will not change significantly during the data acquisition. The performance of the proposed algorithm is demonstrated by the numerical experiments using the synthetic data of crosshole dc resistivity tomography.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2000.09a
• /
• pp.112-131
• /
• 2000

In spite of its many strong points, pole-pole electrode configuration is not often used for ground electrical resistivity survey except for bore hole survey as normal logging and for archaeological investigation. Above all, poor spatial resolution of pole-pole survey may be responsible for this. But recent experiences so far gained by the present authers lead them to think that pole-pole survey can be at least a viable means of reconnaissance survey in near-surface conductive environment and an effective interpretational scheme may augment its resolution. As well known, a response of any other electrode configuration is a linear combination of pole-pole responses. Based on this principle of linear superposition and the principle of reciprocity, the other 'responses' can be derived with simple additions and subtractions of pole-pole responses. Though such responses are not always correct due to the adverse effects of noises, combined with the potential decay curves, they can be helpful to interprete better pole-pole survey data especially in connection of the resolution. This can be comparable to the use of the first or second derivatives of gravity and magnetic intensity surveys.

• Journal of the Korean Geophysical Society
• /
• v.6 no.1
• /
• pp.39-46
• /
• 2003

Twin stone pagodas of the ruins of Kamunsa temple at Kyongju city, Kyungsangbukdo were believed to be built in 682 during the Unified Shilla Kingdom. The 13.4-m-high granodiolite pagodas with the base of 6.78 m x 4.4 m are the largest three-story stone pagoda in Korea. The western pagoda which was re-organized in 1959 is observed to be on the process of severe weathering. Also, some stone contacts are represented by the shape of sharp chevron, which is probably caused by the uneven loading due to the structural unbalance. For the structure-safety diagnosis of the western pagoda, it is necessary to understand its site characteristics and surrounding subsurface environment. Combined geophysical survey such as seismic and resistivity methods was carried out around the western pagoda. The range of 55∼350 Ωm is shown around the pagoda from the electrical resistivity mapping by the Wenner method. The higher resistivities occur the southwestern area, while the lower (<100 Ωm) values indicating the weaker subsurface appear to be on the northeastern area. This result coincides with the measurement of a leaning angle of the pagoda. Along 6 seismic lines, about 3-m-thick uppermost section around the pagoda shows the P-wave velocity of 200∼700 m/s from the refraction survey. Based on the integrated geophysical survey, the foundation of the pagoda is estimated to be in the form of 11-m-side square down to the depth of 3 m.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.56-63
• /
• 2002

Tracing the history of study, problems of seawater intrusion are commonly investigated with electrical techniques because seawater saturated zone is indicative of the low resistivity anomaly. There we, however, silt and mud layers in the western and southern coastal areas of Korea, so we may make a mistake in case we determine seawater intrusion only with resistivity survey. Hence, reference IP survey was carried out in Kimje, Jeollabuk-Do and Youngkwang, Jeollanam-Do in order to decide whether or not the area is under the influence of seawater intrusion. With the use of a electric field cable to minimize EM coupling, we obtained more accurate results by appling reference If technique measuring simultaneously wavelet of current as well as potential. With the aid of reference IP technique, it is possible that we can exactly evaluate seawater intrusion by discriminating seawater saturated area (no IP effect) from very highly conductive layer composed of clay mineral (high IP effect).

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.309-315
• /
• 2002

Conventional electromagnetic (EM) method using small loops as a source and receiver has been used in detection of conductive buried objects like a metal detector or in qualitative estimation of the subsurface conductivity variation. Recently, however, since detection of buried objects and imaging of the subsurface conductivity distribution in a relatively conductive area are in a high demand for environmental and engineering purposes, the quantitative interpretation technique of EM data is actively studied. In this regard, we introduce a brief principle of EM survey and show an example of the detection of buried conductive material and imaging of the subsurface conductivity distribution based on data measured at a test survey area. Through this study, we show that multi-frequency EM surveys using small loops may be a good solution to give quick and detail information of subsurface in a conductive survey area.

• The Journal of Engineering Geology
• /
• v.18 no.2
• /
• pp.127-136
• /
• 2008

Geophysical survey for establishing a wide site for the distribution of water content, wetting front infiltration due to the rainfall, and distribution of groundwater level has been performed by using 8round penetration radar (GPR) method, electrical resistivity method, and so on. On the other hand, a narrow area survey was performed to use a permittivity method such as time domain reflectometry, frequency domain reflectometry, and amplitude domain reflectometry methods for estimating volumetric water content, soil density, and concentration of contaminant in surface and subsurface. The permittivity methods establish more corrective physical parameters than different found survey technologies mentioned above. In this study for establishment of infiltration behaviors for wetting front in the unsaturated soil caused by an artificial rainfall, soil physical parameters for volumetric water content, pore water pressure, and pore air pressure were measured by FDR measurement device and pore water pressure meter which are installed in the unsaturated weathered granite soil with different depths. Consequently, the authors were proposed to a new establishment method for analyzing the variations of volumetric water content and wetting front infiltration from the responses of infiltrating pore water in the unsaturated soil.

• Geophysics and Geophysical Exploration
• /
• v.19 no.4
• /
• pp.212-219
• /
• 2016

The development of more advanced geophysical exploration techniques is necessary because the orebodies as yet discovered are increasingly changing in characteristics from shallow/high-grade to deep/low-grade. In this work, laboratory measurement of physical properties of rock samples and a field survey and interpretation of spectral induced polarization (SIP) have been conducted in a skarn metallic deposit, Gagok mine. The purpose of this study is that the applicability of SIP in the exploration of skarn metallic deposits is verified by the comprehensive interpretation between SIP characteristics of rocks obtained from the laboratory measurements and inverted survey results from the field data. In order to understand the SIP characteristics of each lithology, the data of eighty nine rock samples utilized in the previous studies were revaluated. The field survey was carried out using frequency of 0.25 Hz along a survey line designed for intersecting lithological boundaries and evaluating mineralized zones. The mineralized rocks were more conductive (low-resistivity) and capacitive (high-chargeability or strong-phase) than other rocks. Thus, SIP can be one of the very useful tools for the mineral exploration of the skarn deposits.

• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.73-77
• /
• 2021

Induced polarization (IP) is useful for mineral exploration and hydrogeological studies by visualizing the electrochemical reactions at the interface between polarized minerals and groundwater. Frequency-domain IP (FDIP) is not actively applied to field surveys because it takes longer to acquire data, despite its higher data quality than conventional time-domain IP. However, data quality is more important in current mineral exploration as the targets gradually shift to deep or low-grade ore bodies. In addition, the measurement time reduced by automated instrumentation increases the potential for FDIP field applications. Therefore, we demonstrate that FDIP can detect mineral exploration targets by performing geophysical logging in the boreholes of a skarn deposit, in South Korea. Alternating current (AC) resistivity, percent frequency effect (PFE) and metal factor (MF) were calculated from impedance values obtained at two different frequencies. Skarn zones containing magnetite or pyrite showed relatively low AC resistivity, high PFE, and high MF compared to other zones. Therefore, FDIP surveys are considered to be useful for mineral exploration.