Bayesian decision theory allows for expressing the probability of a particular class given an observed x as
P(Cj│x) = P(x,Cj)/P(x) = (P(x│Cj) P(Cj))/P(x)
where P (x, Cj) is the joint probability of x and Cj, P (x | Cj) denotes the conditional probability of x given Cj. In order to establish a lithoclass x (i.e. gas sand, oil sand, brine sand, shale, etc.), a rock diagnostic is performed where, using the available well log data, a set of crossplots is generated utilizing the inversion output attributes (i.e. AI, SI, Vp/Vs, Poisson’s ratio, Lambda-rho, Mu-rho, etc.) and petrophysical logs (i.e. VCL, PHIT, Sw, etc.). The objective is to determine which elastic attributes best separate out the desired lithoclasses in the crossplot analysis. Last, probability density functions (pdfs) are generated from the well log data and are applied to the selected seismic inversion output cubes to produce class and probability cubes. At each sample, there is an associated probability which represents the likelihood of each class defined based on the cross plot analysis explained above. Several examples as well as a more comprehensive explanation of Bayesian statistics is given in the book Quantitative Seismic Interpretation by Avseth et al. (2005).
The next figure shows a map view of pay sand probability horizon slice. The result is part of a reservoir characterization study in the Burgos Basin, which includes a spatial lithology distribution based on lithoclasses generated using rock physics and Bayesian statistics. The lithology prediction was preceded by prestack simultaneous AVO inversion. The study was carried out to determine the lateral extent of the producing gas-saturated sand which has been penetrated twice successfully (Wells A and B) based solely on seismic amplitudes. Well C was a dry hole, which inspired the operator to find a different methodology to accurately characterize the reservoir sand.
The result may be interpreted as: given the information known from wells A and B, there is an 88.0% probability that hydrocarbons exist at Well A, which is where we know the well produced gas; the same goes for Well B, which shows a 61.5% probability. In Well C, we know there are no hydrocarbons and has a 47.6% probability. One may assign a conservative threshold of ~55% – 60% to determine whether gas saturated sand is present or not.