  Mechanistic induced polarization (IP) models describe the relationships
between the intrinsic properties of geomaterials and their frequency-dependent
complex conductivity spectra. However, the uncertainties associated with
estimating petrophysical properties from IP data are still poorly understood.
Therefore, practitioners rarely use mechanistic models to interpret actual IP
data. We propose a framework for critically assessing any IP model's
sensitivity and parameter estimation limitations. The framework consists of a
conditional variational autoencoder (CVAE), an unsupervised Bayesian neural
network specializing in data dimension reduction and generative modeling. We
train the CVAE on the IP signatures of synthetic mixtures of metallic mineral
inclusions in electrolyte-filled host geomaterials and describe the effect of
data transformations on the model. First, the CVAE's Jacobian reveals the
relative importance of each petrophysical property for generating spectral IP
data. The most critical parameters are the conductivity of the host, the
volumetric content of the inclusions, the characteristic length of the
inclusions, and the permittivity of the host. The inclusions' diffusion
coefficient, permittivity, and conductivity, as well as the host's diffusion
coefficient, only have marginal importance for generative IP modeling. A
parameter estimation experiment yields the standardized accuracy of
petrophysical properties using various model constraints scenarios and
corroborates the sensitivity analysis results. Finally, we visualize the
effects of data transformations and model constraints on the petrophysical
parameter space. We conclude that a common logarithm data transformation yields
optimal parameter estimation results and that constraining the electrochemical
properties of the geomaterial improves estimates of the characteristic length
of its metallic inclusions and vice versa.
