• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.149-157
• /
• 2021

Given that induced polarization (IP) and direct current (DC) resistivity surveys are similar in terms of data acquisition, most DC resistivity systems are equipped with a time-domain IP data acquisition function. In addition, the time-domain IP data include the DC resistivity values. As such, IP and DC resistivity data are intimately linked, and the inversion of IP data is a two-step process based on DC resistivity inversions. Nevertheless, IP surveys are rarely applied, in contrast to DC resistivity surveys, as proper inversion software is unavailable. In this study, through numerical modeling and inversion experiments, we analyze the problems with the conventional inverse mapping technique used to invert time-domain IP data. Furthermore, we propose a modified inverse mapping technique that can effectively suppress inversion artifacts. The performance of the technique is confirmed through inversions applied to synthetic IP data.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.139-148
• /
• 2021

Negative apparent chargeability data in time-domain induced polarization (IP) survey are very often considered as noise and rejected before the inversion stage. Negative IP data can, however, occur naturally as a consequence of the distribution of chargeable zones in the ground. In some cases, the negative apparent chargeability values may account for most of the data measured. Negative IP data are caused by the geometry of chargeable zones and electrode positions. Negative apparent chargeability data appear most frequently when a dipole-dipole array is used. In this study, the effect of negative apparent chargeability data on inversion results is analyzed through the numerical 2D time-domain IP modeling and nonlinear inversion. The results demonstrate clearly that negative apparent chargeability data have to be included in the inversion as they contain important information on the distribution of subsurface chargeability.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.202-208
• /
• 2021

Induced polarization (IP) was first published in a Korean academic journal in 1973, and it was soon applied to coal and metal ore exploration. Then, in universities and research institutes, IP modeling studies using the finite element approach and experimental studies on IP responses for artificial samples were conducted. In the mid-1980s, the spectral IP (SIP) measurement module was introduced to Korea, and physical scale modeling and inversion approaches were developed. Due to the decline of the mineral resource industry, this method was not actively applied. However, the SIP method was not applied In the 1990s, IP exploration was applied in the investigation of hydrothermal deposits of sulfide minerals and bentonite mineralization zones, as well as to areas where the groundwater was contaminated by intruding seawater. In the 2000s, three-dimensional inversion of the IP approach was developed, and high-precision geophysical exploration was required to secure domestic and overseas mineral resources, so SIP experiments on rock samples and approaches for field exploration were developed. The SIP approach was proven useful for the exploration of metal deposits containing sulfide minerals by applying it to explore the mineralization zone of gold-silver deposits in the Haenam region. The IP method is considered to be effective in exploring critical minerals (lithium, cobalt, and nickel) in high-tech industries. It also is expected to be useful for environmental and geotechnical investigations.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.383-389
• /
• 2003

Clay minerals show a distinct induced-polarization phenomenon, which is one of the most important factors for predicting groundwater flow and contaminant transport. This paper presents a step-by-step process to estimate Cole-Cole parameters from spectral induced-polarization (IP) data measured on the surface of three-dimensional earth. First, the inversion of low-frequency resistivity survey data is made to identify the dc resistivity ${\rho}_dc$ of a volume having IP effects. The other parameters, chargeability m, time constant $\tau$, and frequency dependence c, are sought for the polarizable volume. Next, using multi-frequency data, c can be obtained as high or low asymptotes of the slope of log phase vs. log frequency. Further, for low m, intrinsic $\tau$ is approximated by apparent one, ${\tau}_a$, which is derived from the relation ${{\omega}{\tau}}_a$=1 at an angular frequency $\omega$, where the imaginary component of spectral IP data has an extreme value. Finally, to obtain intrinsic m a two-step linearized procedure has been derived. For a body of given $\tau$ and c, forward modeling with a progression of m values yields a plot of observed vs. intrinsic imaginary components for a frequency. Since this plot is essentially linear, to extract the intrinsic imaginary component is quite simple with an observed value. Using the plot of intrinsic imaginary component vs. m, intrinsic m is determined. We present a synthetic example to illustrate that the Cole-Cole parameters can be recovered from spectral IP data.

• Geophysics and Geophysical Exploration
• /
• v.20 no.2
• /
• pp.100-113
• /
• 2017

Induced polarization (IP) method is based on the measurement of a polarization effect known as overvoltage of the ground. IP techniques have been usually used to find mineral deposits, however, nowadays widely applied to hydrogeological investigations, surveys of groundwater pollution and foundation studies on construction sites. IP surveys can be classified by its source type, i.e., time-domain IP estimating chargeability, frequency-domain IP measuring frequency effect (FE), and complex resistivity (CR) and spectral IP (SIP) measuring complex resistivity. Recently, electromagnetic-based IP has been studied to avoid the requirement for spike electrodes to be placed in the ground. In order to understand IP methods in this study, we: 1) classify IP surveys by source type and measured data and illustrate their basic theories, 2) describe historical development of each IP forward modeling and inversion algorithm, and finally 3) introduce various case studies of IP measurements.