A computer implemented method for identifying reservoir facies in a subsurface region includes obtaining a set of seismic data points of both petrophysical and geophysical parameters relating to the subsurface region, identifying one or more correlated clusters of petrophysical parameters, generating, from the one or more correlated clusters of petrophysical parameters, one or more corresponding multi-dimensional clusters of seismic data points, storing, in a facies database, a multi-dimensional cluster center point for at least one multi-dimensional clusters, and recursively splitting the multi-dimensional clusters into distinct sub-clusters of seismic data points corresponding to facies types.