Bayesian neural network and Bayesian physics-informed neural network via variational inference for seismic petrophysical inversion

Abstract

Bayesian neural network via variational inference (BNN-VI) aims to estimate probability distribution of the training parameters of the neural network, which are difficult to compute with traditional methods, by proposing a family of densities and finding the candidates that are close to the target. Variational inference is a Bayesian method that provides a faster tool than Markov chain Monte Carlo sampling algorithms. However, to provide an accurate estimation of the most likely models, BNN-VI requires sufficient training data, which is not always available in geophysical inverse problems. We propose a method that combines a Bayesian approach with a physics-informed neural networks algorithms via variational inference, namely Bayesian physics-informed neural networks via variational inference (BPINN-VI) and apply it to a geophysical inverse problem for the estimation of petrophysical properties from seismic data. The method is based on a Bayesian network approach in which the objective function is defined by the Kullback–Leibler (KL) divergence. In addition, the physical relations between data (seismic measurements) and model variables (petrophysical properties) are embedded into neural networks through a physics-dependent loss term. We tested the proposed method on a pre-stack seismic dataset from the North Sea with limited direct measurements of the target petrophysical properties from well logs. The outcome is the most likely model of petrophysical properties. Compared to BNN-VI, the BPINN-VI results show higher correlation and R2 coefficients, and lower mean square error as well as lower uncertainty and higher lateral continuity.