Aggregate thermostatically controlled loads (TCLs) are good candidates for providing load following services in power systems. This paper is concerned with the modeling, evaluation and control problems of a population of heterogeneous TCLs. Specifically, the heterogeneous population is divided into multiple homogeneous clusters and each cluster, i.e., TCL aggregator, is modeled by an approximated three-input single-output (TISO) state space model. Here, the aggregators serve as a bridge connecting the load utility and the terminal TCLs, which have their own decision makers and are responsible for aggregate estimation and command issuing. And aggregate evaluation is carried out for the aggregator so as to provide the aggregate regulation capacities and ramping rates, which is useful for setting of the reference power trajectory. Based on the established control model, we furthermore propose a hierarchical centralized control algorithm for a bus load utility to regulate all TCLs inside it so as to provide load following service, while not affecting the customers' comfort levels. Finally, simulation results with respect to a common bus load are provided to demonstrate the effectiveness of the proposed aggregate modeling and the centralized load following strategy.