Assessing the uncertainty in reservoir performance is a necessary step during the exploration phase. To examine the uncertainty in flow response, a large set of realizations must be processed. There are several stochastic geostatistical algorithms capable of simulating multiple equiprobable realizations. Although these can show us the possible realities highlighting the spatial uncertainty, their handling is time- and CPU-consuming during the later processes, such as flow simulations. Consequently, only a small number of realizations can be post-processed in industrial practice. The purpose of this work is to develop a method, which will reduce the huge number of realizations in a way that the remaining ones retain the spatial uncertainty of a reservoir’s flow behavior, as would be demonstrated by a larger set of realizations. To solve this problem, ranking methods can be applied. Traditional ranking techniques, such as probability selection, are highly dependent on the applied static properties. In this paper, an alternative selection method is parameterized for measuring the pairwise dissimilarity between geostatistical models, with a distance function based on the hydrodynamic properties of the hydrocarbon reservoirs. The effectivity of the method is highly dependent upon the selected criteria. Thus, the distance function refers to the flow responses and allows visualizing the space of uncertainty through multidimensional scaling. A kernel transformation of the MDS data set is required to obtain a feature space where the K-means algorithm can discover non-linear structures in the basic data set. The final step of the method is the selection of the Earth models closest to the cluster centers. This tool allows for the selection of a subset of representative realizations, containing similar properties to the larger set.