Gaussian Process-based Approach for Bilevel Optimization in the Power System -- A Critical Load Restoration Case