TY - JOUR

T1 - Conversion of Induced Polarization Data and Their Uncertainty from Time Domain to Frequency Domain Using Debye Decomposition

AU - Hase, Joost

AU - Gurin, Grigory

AU - Titov, Konstantin

AU - Kemna, Andreas

PY - 2023/7/17

Y1 - 2023/7/17

N2 - The time-domain (TD) induced polarization (IP) method is used as an extension of direct current (DC) resistivity measurements to capture information on the ability of the subsurface to develop electrical polarization. In the TD, the transient voltage decay is measured after the termination of the current injection. To invert tomographic TD IP data sets into frequency-domain (FD) models of complex electrical resistivity, a suitable approach for converting TD IP transients and their corresponding uncertainties into the FD is essential. To apply existing FD inversion algorithms to TD IP measurements, a conversion scheme must transform the measured decay curves into FD impedances and also propagate the corresponding measurement uncertainty from the TD to the FD. Here, we present such an approach based on a Debye decomposition (DD) of the decay curve into a relaxation-time distribution and the calculation of the equivalent spectrum. The corresponding FD data error can be obtained by applying error propagation through all of these steps. To accomplish the DD we implement a non-linear Gauss–Newton inversion scheme. We test the conversion scheme in a synthetic study and demonstrate its application to field data on a tomographic TD IP data set measured on the Maletoyvaemskoie ore field (Kamchatka, Russia). The proposed conversion scheme yields accurate impedance data for relaxation processes, which are resolved by the TD measurements. The error propagation scheme provides a reasonable FD uncertainty estimate, as confirmed by a Monte Carlo analysis of the underlying parameter distributions.

AB - The time-domain (TD) induced polarization (IP) method is used as an extension of direct current (DC) resistivity measurements to capture information on the ability of the subsurface to develop electrical polarization. In the TD, the transient voltage decay is measured after the termination of the current injection. To invert tomographic TD IP data sets into frequency-domain (FD) models of complex electrical resistivity, a suitable approach for converting TD IP transients and their corresponding uncertainties into the FD is essential. To apply existing FD inversion algorithms to TD IP measurements, a conversion scheme must transform the measured decay curves into FD impedances and also propagate the corresponding measurement uncertainty from the TD to the FD. Here, we present such an approach based on a Debye decomposition (DD) of the decay curve into a relaxation-time distribution and the calculation of the equivalent spectrum. The corresponding FD data error can be obtained by applying error propagation through all of these steps. To accomplish the DD we implement a non-linear Gauss–Newton inversion scheme. We test the conversion scheme in a synthetic study and demonstrate its application to field data on a tomographic TD IP data set measured on the Maletoyvaemskoie ore field (Kamchatka, Russia). The proposed conversion scheme yields accurate impedance data for relaxation processes, which are resolved by the TD measurements. The error propagation scheme provides a reasonable FD uncertainty estimate, as confirmed by a Monte Carlo analysis of the underlying parameter distributions.

KW - Debye decomposition

KW - complex resistivity tomography

KW - induced polarization

UR - https://www.mendeley.com/catalogue/1fb82ea4-bd3b-3214-aa77-efe41c86b42f/

U2 - 10.3390/min13070955

DO - 10.3390/min13070955

M3 - Article

VL - 13

JO - Minerals

JF - Minerals

SN - 2075-163X

IS - 7

M1 - 955

ER -