Bayesian inference of petrophysical properties with generative spectral induced polarization models

Charles L. Bérubé,Frédérique Baron
DOI: https://doi.org/10.1190/geo2022-0495.1
2022-08-02
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.
Geophysics
What problem does this paper attempt to address?