Machine learning (ML) algorithms, such as principal component analysis, independent component analysis, self-organizing maps, and artificial neural networks, have been used by geoscientists to not only accelerate the interpretation of their data, but also to provide a more quantitative estimate of the likelihood that any voxel belongs to a given facies. Identifying the best combination of attributes needed to perform either supervised or unsupervised ML tasks continues to be the most-asked question by interpreters. In the past decades, stepwise regression and genetic algorithms have been used together with supervised learning algorithms to select the best number and combination of attributes. For reasons of computational efficiency, these techniques do not test all of the seismic attribute combinations, potentially leading to a suboptimal classification. In this study, we have developed an exhaustive probabilistic neural network (PNN) algorithm that exploits the PNN’s capacity in exploring nonlinear relationships to obtain the optimal attribute subset that best differentiates target seismic facies of interest. We determine the efficacy of our proposed workflow in differentiating salt from nonsalt seismic facies in a Eugene Island seismic survey, offshore Louisiana. We find that from seven input candidate attributes, the exhaustive PNN is capable of removing irrelevant attributes by selecting a smaller subset of four seismic attributes. The enhanced classification using fewer attributes also reduces the computational cost. We then use the resulting facies probability volumes to construct the 3D distribution of the salt diapir geobodies embedded in a stratigraphic matrix.