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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.231-236
• /
• 2007

Electrical resistivity, self-potential and time-domain induced polarization methods were conducted for study of electrical responses on vein-type sulfides ore, which is intruding limestone and dolomite of Ordovician, of Geumpung mine located in Dojeon-ri, Susan-myeon, Jecheon-si, Chungbuk. Sulfides bearing chalcopyrite, pyrrhotite and galena etc. are deposited in disseminated or vein-type. Good result that resistivity and self-potential surveys detect high grade-estimated mineralized zone located in upper part of existing low grade ore zone is acquisited and is to some extent consistent with induced polarization. Furthmore, a new mineralization zone directing EW is detected.

• Economic and Environmental Geology
• /
• v.19 no.3
• /
• pp.193-197
• /
• 1986

A simple and fast solution is derived to evaluate the effects of time-domain electromagnetic coupling in induced polarization surveys on a uniform earth. The simplified solution gives an explicit statement of the dependence of time-domain electromagnetic coupling on the model parameters, and yields sufficiently accurate results for most situations encountered in practice. The co-linear dipole-dipole and Wenner arrays are used as examples in this paper, but th numerical solution can be applied to any electrode configuration.

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

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

• Proceedings of the Korean Vacuum Society Conference
• /
• 2013.08a
• /
• pp.99-99
• /
• 2013

In this talk, I will present two research works in progress, which are: i) mapping of piezoelectric polarization and associated charge density distribution in the heteroepitaxial InGaN/GaN multi-quantum well (MQW) structure of a light emitting diode (LED) by using inline electron holography and ii) in-situ observation of the polarization switching process of an ferroelectric Pb(Zr1-x,Tix)O3 (PZT) thin film capacitor under an applied electric field in transmission electron microscope (TEM). In the first part, I will show that strain as well as total charge density distributions can be mapped quantitatively across all the functional layers constituting a LED, including n-type GaN, InGaN/GaN MQWs, and p-type GaN with sub-nm spatial resolution (~0.8 nm) by using inline electron holography. The experimentally obtained strain maps were verified by comparison with finite element method simulations and confirmed that not only InGaN QWs (2.5 nm in thickness) but also GaN QBs (10 nm in thickness) in the MQW structure are strained complementary to accommodate the lattice misfit strain. Because of this complementary strain of GaN QBs, the strain gradient and also (piezoelectric) polarization gradient across the MQW changes more steeply than expected, resulting in more polarization charge density at the MQW interfaces than the typically expected value from the spontaneous polarization mismatch alone. By quantitative and comparative analysis of the total charge density map with the polarization charge map, we can clarify what extent of the polarization charges are compensated by the electrons supplied from the n-doped GaN QBs. Comparison with the simulated energy band diagrams with various screening parameters show that only 60% of the net polarization charges are compensated by the electrons from the GaN QBs, which results in the internal field of ~2.0 MV cm-1 across each pair of GaN/InGaN of the MQW structure. In the second part of my talk, I will present in-situ observations of the polarization switching process of a planar Ni/PZT/SrRuO3 capacitor using TEM. We observed the preferential, but asymmetric, nucleation and forward growth of switched c-domains at the PZT/electrode interfaces arising from the built-in electric field beneath each interface. The subsequent sideways growth was inhibited by the depolarization field due to the imperfect charge compensation at the counter electrode and preexisting a-domain walls, leading to asymmetric switching. It was found that the preexisting a-domains split into fine a- and c-domains constituting a $90^{\circ}$ stripe domain pattern during the $180^{\circ}$ polarization switching process, revealing that these domains also actively participated in the out-of-plane polarization switching. The real-time observations uncovered the origin of the switching asymmetry and further clarified the importance of charged domain walls and the interfaces with electrodes in the ferroelectric switching processes.

• Tunnel and Underground Space
• /
• v.33 no.6
• /
• pp.483-494
• /
• 2023

Measurement of the physical properties of rocks or minerals is an important factor in determining the distribution of the underground medium as well as mineral resource investigations. Resistivity and induced polarization, which are widely used in Korea, are methods for measuring electrical properties, which are representative properties of obtaining subsurface information. In order to precisely analyze the exploration data obtained from various sites, it is important to accurately measure the material properties. Electrical properties of rock is measured using two-electrode or four-electrode method. Compared to the four-electrode method, the two-electrode method is generally used because it is very easy to contact the sample and the electrode, but there is a problem in that the impedance of the electrode and the sample is measured together. In this study, the time-domain the induced polarization effects were measured using the 2-electrode method and the 4-electrode method for artificial samples mixed with graphite and cement having induced polarization characteristics, and the results were compared. Although the 4-electrode method has difficulties in installing potential electrodes, it was confirmed that it is effective in measuring electrical properties because it can reduce the problem caused by the impedance of potential electrodes compared to the 2-electrode method.