With the penetration of renewable energy increasing, the power system requires higher flexibility of power regulation. Virtual power plants can aggregate distributed flexible loads to improve the utilization of distributed renewable energy. In this paper, a bi-level scheduling model for virtual power plants with a large number of distributed thermostatically controlled loads and intermittent renewable energy is established to reduce the net exchange power deviation caused by the forecast error of renewable energy. The upper level optimizes the exchange power curve and reduces the imbalance costs in intraday, while the lower level tracks the optimized power curve in real-time to complete the regulation target. Static and dynamic aggregation method reflecting the regulation characteristics of aggregated thermostatically controlled loads is proposed and applied in lower/upper level, respectively. In addition, a two-step simplified strategy is proposed to solve the mixed integer nonlinear programming in upper level. Simulation results show that the proposed method can reduce the maximum imbalance power, and it is not affected by parameters heterogeneity, which is suitable for virtual power plants with diversified users.