Learning relaxation time distributions from spectral induced polarization data with a complex-valued variational autoencoder

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Bérubé, Charles L., Gagnon, Sébastien, Nagasingha, Lahiru M. A., Gagnon, Jean-Luc, Kenko, E. Rachel, Ghanati, Reza, Baron, Frédérique
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913017986285568
author Bérubé, Charles L.
Gagnon, Sébastien
Nagasingha, Lahiru M. A.
Gagnon, Jean-Luc
Kenko, E. Rachel
Ghanati, Reza
Baron, Frédérique
author_facet Bérubé, Charles L.
Gagnon, Sébastien
Nagasingha, Lahiru M. A.
Gagnon, Jean-Luc
Kenko, E. Rachel
Ghanati, Reza
Baron, Frédérique
contents Spectral induced polarization (SIP) is a geophysical method used to characterize subsurface materials. It measures the frequency-dependent complex resistivity of rocks and soils through the application of a small alternating current in the subsurface or in laboratory samples. Debye decomposition (DD) is a standard method for analyzing and interpreting SIP data, as it allows estimation of the relaxation time distribution (RTD) of geomaterials. However, conventional DD approaches treat measurements independently, work in real-valued spaces despite the complex-valued nature of SIP data, and provide limited uncertainty quantification. These limitations reduce the effectiveness of conventional DD on heterogeneous datasets. We reformulate DD as an unsupervised machine learning problem and introduce a conditional variational autoencoder (CVAE) that learns a shared mapping from resistivity spectra to continuous RTDs. The model is validated on a dataset comprising 140 laboratory and field SIP measurements of granular mixtures, mineralized rocks, and cementitious materials. The CVAE operates in complex-valued data space and achieves reconstruction errors of 0.45 % and 0.24 % for the imaginary and phase components of resistivity, respectively, with statistically significant improvements over conventional methods (p-values of 4x10^-6 and 2x10^-3). The inferred RTDs are stable and physically consistent, and their total chargeability and mean relaxation time correlate with polarizable grain content and grain size, respectively, with coefficients of determination up to 0.98. An additional contribution of the proposed method is the learned latent representation, which organizes SIP spectra into a structured space. Unsupervised clustering in a two-dimensional projection of this space improves the Davies--Bouldin index by nearly a factor of three relative to conventional RTD parameters.
format Preprint
id arxiv_https___arxiv_org_abs_2603_13973
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Learning relaxation time distributions from spectral induced polarization data with a complex-valued variational autoencoder
Bérubé, Charles L.
Gagnon, Sébastien
Nagasingha, Lahiru M. A.
Gagnon, Jean-Luc
Kenko, E. Rachel
Ghanati, Reza
Baron, Frédérique
Geophysics
Spectral induced polarization (SIP) is a geophysical method used to characterize subsurface materials. It measures the frequency-dependent complex resistivity of rocks and soils through the application of a small alternating current in the subsurface or in laboratory samples. Debye decomposition (DD) is a standard method for analyzing and interpreting SIP data, as it allows estimation of the relaxation time distribution (RTD) of geomaterials. However, conventional DD approaches treat measurements independently, work in real-valued spaces despite the complex-valued nature of SIP data, and provide limited uncertainty quantification. These limitations reduce the effectiveness of conventional DD on heterogeneous datasets. We reformulate DD as an unsupervised machine learning problem and introduce a conditional variational autoencoder (CVAE) that learns a shared mapping from resistivity spectra to continuous RTDs. The model is validated on a dataset comprising 140 laboratory and field SIP measurements of granular mixtures, mineralized rocks, and cementitious materials. The CVAE operates in complex-valued data space and achieves reconstruction errors of 0.45 % and 0.24 % for the imaginary and phase components of resistivity, respectively, with statistically significant improvements over conventional methods (p-values of 4x10^-6 and 2x10^-3). The inferred RTDs are stable and physically consistent, and their total chargeability and mean relaxation time correlate with polarizable grain content and grain size, respectively, with coefficients of determination up to 0.98. An additional contribution of the proposed method is the learned latent representation, which organizes SIP spectra into a structured space. Unsupervised clustering in a two-dimensional projection of this space improves the Davies--Bouldin index by nearly a factor of three relative to conventional RTD parameters.
title Learning relaxation time distributions from spectral induced polarization data with a complex-valued variational autoencoder
topic Geophysics
url https://arxiv.org/abs/2603.13973