• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.164-170
• /
• 2021

Frequency-domain and time-domain induced polarization methods can provide spectral information about subsurface media. Analysis of spectral characteristics has been studied mainly in the frequency-domain, however, time-domain induced polarization research has recently become popular. In this study, assuming a homogeneous half-space model, an inversion method was developed to extract Cole-Cole parameters from the measured secondary potential or electrical resistivity. Since the Cole-Cole parameters of chargeability, time constant, and frequency index are not independent of each other, various problems, such as slow convergence rate, initial model problem, local minimum problem, and divergence, frequently occur when conventional nonlinear inversion is applied. In this study, we developed an effective inversion method using the initial model close to the true model by introducing a grid search method. Finally, the validity of the developed inversion method was verified using inversion experiments.

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

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.158-163
• /
• 2021

The frequency-domain induced polarization (IP) response based on Cole-Cole model is expressed as a simple equation in close form. However, it is difficult to compute the time-domain IP response based on Cole-Cole model or any other relaxation model because it cannot be written in closed form. In this study, using numerical experiments, we compared three numerical methods for calculating the time-domain IP response of the Cole-Cole model asymptotically: series expansion, digital linear filtering and Fourier transform. The series expansion method is inadequately accurate for certain time values and converges very slowly. A digital linear filter specially designed to calculate the time-domain IP response does not present the desired accuracy, especially at later times. The Fourier transform method can overcome the abovementioned problems and present the time-domain IP response with adequate accuracy for all time values, even though more computing time is required.

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

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.171-179
• /
• 2021

We outline a process for estimating Cole-Cole parameters from time-domain induced polarization (IP) data. The IP transients are all inverted to 2D Cole-Cole earth models that include resistivity, chargeability, relaxation time, and the frequency exponent. Our inversion algorithm consists of two stages. We first convert the measured voltage decay curves into time series of current-on time apparent resistivity to circumvent the negative chargeability problem. As a first step, a 4D inversion recovers the resistivity model at each time channel that increases monotonically with time. The desired intrinsic Cole-Cole parameters are then recovered by inverting the resistivity time series of each inversion block. In the second step, the Cole-Cole parameters can be estimated readily by setting the initial model close to the true value through a grid search method. Finally, through inversion procedures applied to synthetic data sets, we demonstrate that our algorithm can image the Cole-Cole earth models effectively.

• Current Optics and Photonics
• /
• v.5 no.1
• /
• pp.66-71
• /
• 2021

The resonance characteristics of H-shaped metamaterials, whose widths were varied while keeping the height constant, were investigated in the terahertz (THz) frequency range. The H-shaped metamaterials were numerically analyzed in two modes in which the polarization of the incident THz electric field was either parallel or perpendicular to the width of the H-shaped structure. The resonant frequency of the metamaterial changed stably in each mode, even if only the width of the H shape was changed. The resonant frequency of the metamaterial operating in the two modes increases without significant difference regardless of the polarization of the incident electromagnetic wave as the width of the H-shaped metamaterial increases. The electric field distribution and the surface current density induced in the metamaterial in the two modes were numerically analyzed by varying the structure ratio of the metamaterial. The numerical analysis clearly revealed the cause of the change in the resonance characteristics as the width of the H-shaped metamaterial changed. The efficacy of the numerical analysis was verified experimentally using the THz-TDS (time-domain spectroscopy) system. The experimental results are consistent with the simulations, clearly demonstrating the meaningfulness of the numerical analysis of the metamaterial. The analyzed resonance properties of the H-shaped metamaterial in the THz frequency range can be applied for designing THz-tunable metamaterials and improving the sensitivity of THz sensors.

• Economic and Environmental Geology
• /
• v.18 no.4
• /
• pp.301-308
• /
• 1985

I.P exploration is conducted mainly at the alluvium covered granodiorite stock of the Red Hill area at the Dongjeom coal mine, employing Canadian Mcphar equipments of variable frequencies-domain method along a total of 8.5 survey line. Mineralization zone is found by LP anomalies along the I.P the profile of frequency mode. Comparing with the past drilling data, the cause of each anomly is furthermore identified as copper bearing mineralized zones. As alteration and mineralization cover all over the Red-Hill altered grandiorite, copper bearing sulfide veinlets and strings filled out the fractured and altered zones.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.180-193
• /
• 2021

The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.