Ensemble-based analyses are useful to compare equiprobable scenarios of the reservoir models. However, they require a large suite of reservoir models to cover high uncertainty in heterogeneous and complex reservoir models. For stable convergence in ensemble Kalman filter (EnKF), increasing ensemble size can be one of the solutions, but it causes high computational cost in large-scale reservoir systems. In this paper, we propose a preprocessing of good initial model selection to reduce the ensemble size, and then, EnKF is utilized to predict production performances stochastically. In the model selection scheme, representative models are chosen by using principal component analysis (PCA) and clustering analysis. The dimension of initial models is reduced using PCA, and the reduced models are grouped by clustering. Then, we choose and simulate representative models from the cluster groups to compare errors of production predictions with historical observation data. One representative model with the minimum error is considered as the best model, and we use the ensemble members near the best model in the cluster plane for applying EnKF. We demonstrate the proposed scheme for two 3D models that EnKF provides reliable assimilation results with much reduced computation time.