Abstract We propose a simple method to estimate facies probabilities based on statistical analysis of multidimensional crossplots of seismic attributes. Log scale facies flags related to thick sand bodies are created using petrophysics and rock physics analyses. These flags are then used to color crossplots of seismic scale attributes derived from AVO inversion of PP data (V P , V S , and density) and inversion of poststack fast and slow PS components of a 3C-3D survey. We show that by using these five seismic attributes and facies flags (like a colored five-dimensional crossplot), we can estimate the probability of thick sand bodies much better than when we crossplot two attributes at a time. Unlike commonly used approaches to map facies or lithologies from seismic data based on selecting regions in seismic attribute crossplots, our approach accounts properly for overlap among different facies and quantifies the probability of their occurrence.