Application of the Truncated Plurigaussian Method to delineate hydrofacies distribution in heterogeneous aquifers