The majority of the geostatistical realizations ranking methods disregard the production history in selection of realizations, due to its requirement to high simulation run time. They also ignore to consider the degree of linear relationship between the “ranks based on the ranking measure” and “ranks based on the performance parameter” in choosing the employed ranking measure. To address these concerns, we propose an uncertainty quantification workflow, which includes two sequential stages of history matching and realization selection. In the first stage, production data are incorporated in the uncertainty quantification procedure through a history matching process. A fast simulator is employed to find the realizations with consistent flow behavior with the production history data in shorter time, compared to a comprehensive simulator. The selected realizations are the input of the second stage of the workflow, which can be any type of realization selection method. In this study, we used the most convenient realization selection method, i.e. ranking of the realizations. To select the most efficient ranking measure, we investigated the degree of the linear correlation between the ranks based on the several ranking measures and the performance parameter. In addition, due to the shortcomings of the traditional ranking methods in uncertainty quantiles identification, a modified ranking method is introduced. This modification increases the certainty in the probability of the selected realizations. The obtained results on 3D close-to-real synthetic reservoir models revealed capability of the modified ranking method in more accurate quantification of the uncertainty in reservoir performance prediction.