This paper proposes a geostatistical approach for geological modelling and for validating an interpreted geological model, by identifying the areas of an ore deposit with a high probability of being misinterpreted, based on quantitative coregionalised covariates correlated with the geological categories. This proposal is presented through a case study of an iron ore deposit at a stage where the only available data are from exploration drill holes. This study consists of jointly simulating the quantitative covariates with no previous geological domaining. A change of variables is used to account for stoichiometric closure, followed by projection pursuit multivariate transformation, multivariate Gaussian simulation, and conditioning to the drill hole data. Subsequently, a decision tree classification algorithm is used to convert the simulated values into a geological category for each target block and realisation. The determination of the prior (ignoring drill hole data) and posterior (conditioned to drill hole data) probabilities of categories provides a means of identifying the blocks for which the interpreted category disagrees with the simulated quantitative covariates.