Accurate prediction of the structure-dependent interphase drag coefficient in a computational cell is of importance for the simulation of heterogeneous gas-solid flow. However, the cluster diameter in the cell was usually correlated with macroscopic hydrodynamic parameters or simply approximated as a constant value, though there is a cluster size distribution (CSD) at the sub-grid scale. Based on the energy minimization multiscale (EMMS) model and the population balance model (PBM), a CFD-PBM-EMMS integrated model was proposed to account for the effect of CSD on the interphase drag within the computational grid in this article, in which the PBM was used to describe the spatio-temporal evolution of CSD by introducing an EMMS-based cluster growth rate model. The CFD-PBM-EMMS model was validated by the experiments in a pilot-scale circulating fluidized bed riser of Geldart A and/or B particles. Comparing with the two-fluid model using the EMMS drag or homogeneous drag law, the CFD-PBM-EMMS integrated model shows the best agreement with the experiments under various operating conditions. Much effort is being devoted to incorporating the coalescence and breakage kinetics of clusters to the PBM further to improve the simulation accuracy of the CFD-PBM-EMMS integrated model, especially for dense gas-solid flow.