Variational Bayesian inference for the probabilistic model of power load