Bayesian inference of petrophysical properties with generative spectral induced polarization models

TL;DR

This paper introduces a Bayesian framework using a CVAE to assess and improve the estimation of petrophysical properties from spectral IP data, addressing uncertainties in mechanistic models.

Contribution

It presents a novel Bayesian neural network approach for sensitivity analysis and parameter estimation in spectral IP models, enhancing interpretability and reliability.

Findings

Logarithmic data transformation improves parameter estimation accuracy.

Constraining electrochemical properties enhances estimates of inclusion length.

Sensitivity analysis identifies key petrophysical parameters affecting spectral IP data.

Abstract

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.