A Bayesian multivariate model using Hamiltonian Monte Carlo inference to estimate total organic carbon content in shale

The prediction of total organic carbon (TOC) content using geophysical logs is one of the key steps in shale reservoir characterization. Various empirical relations have previously been used for the estimation of TOC content from well-logs; however, uncertainty quantification in the model estimation is often ignored while performing TOC estimation in a deterministic framework. We introduce the problem of TOC estimation in a Bayesian setting with the goal of enhancing the TOC content prediction together with the quantification of the uncertainty in the model prediction. To signify the uncertainty, we draw random samples of model parameters from the posterior distribution by realizing multidimensional stochastic processes within the Hamiltonian Monte Carlo algorithm. The posterior model for the variables that influence TOC estimation is conditioned on the available well-log observations and is further defined by a priori and likelihood distributions. We demonstrate examples of applications of this approach to estimate the TOC content on two real field data sets from the well-known Devonian Duvernay shale of Western Canada and the Silurian shale of the Ahnet Basin. The accuracy in the estimation is arbitrated by comparing the prediction results with those obtained using the two most widely used empirical models. Corroborating the results by the laboratory-measured TOC contents demonstrate that the Bayesian approach offers a more reliable and better confidence in predictions when compared with the empirical models, as it provides additional information on the prediction uncertainty. Finally, the implications of the present approach are derived in terms of depositional environments to characterize the high TOC content zone in the studied organic shale formations.