This invention discloses a method for managing a power load of a grid includes performing a statistic-based distribution estimation of a composite load model using static and dynamic models with Gibbs sampling, and deriving a distribution estimation of load model coefficients.