Method for real-time scheduling of multi-energy complementary micro-grids based on rollout algorithm

Abstract

The invention relates to a method for real-time scheduling of multi-energy complementary micro-grids based on a Rollout algorithm, which is technically characterized by comprising the following steps of: Step 1, setting up a moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs, and establishing constraint conditions for the real-time scheduling; Step 2, establishing a target function of the real-time scheduling; Step 3, dividing a single complete scheduling cycle into a plurality of scheduling intervals, and finding one basic feasible solution meeting the constraint conditions for the real-time scheduling based on a greedy algorithm; and Step 4, finding a solution to the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids by using the Rollout algorithm based on the basic feasible solution from Step 3. With the consideration of the fluctuations in the new-energy outputs, the present invention solves the problems of low speed and low efficiency of a traditional algorithm at the same time, enabling high-speed efficient multi-energy complementary micro-grid real-time scheduling.

Claims

1. A method for real-time scheduling of multi-energy complementary micro-grids based on a Rollout algorithm, characterized by comprising the following steps of: Step 1, setting up a moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs, and establishing constraint conditions for the real-time scheduling; Step 2, establishing a target function of the real-time scheduling for the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with the random new-energy outputs, with the goal of minimum operating cost of a micro-grid system in a moving-horizon Markov decision cycle; Step 3, dividing a single complete scheduling cycle into a plurality of scheduling intervals, and finding one basic feasible solution meeting the constraint conditions for the real-time scheduling based on a greedy algorithm; and Step 4, finding a solution to the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids by using the Rollout algorithm based on the basic feasible solution from Step 3; wherein the constraint conditions established for the real-time scheduling in Step 1 comprises: micro-grid electric equilibrium constraints, storage battery operating constraints, exchange electric power constraints for the micro-grids and a main grid, and electric power output constraints for combined heat and power equipment; the micro-grid electric equilibrium constraints are as follows: $p^G\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{p}_i^c\left(t\right)+p^B\left(t\right)+p^w\left(t\right)=p^D\left(t\right)$ in the formula, t is a time parameter; p.sup.G(t) is exchange electric power for the micro-grids and the main grid at a time t, which is positive during purchasing of electricity from the main grid and negative during selling of electricity to the main grid; N is the quantity of the combined heat and power equipment; p.sub.i.sup.c(t) is output electric power of the i th combined heat and power equipment at the time t; p.sup.B(t) is charging/discharging power of the storage battery at the time t, which is negative during charging and positive during discharging; p.sup.w(t) is generated output of wind power at the time t; and p.sup.D(t) is an electric load demand at the time t; the storage battery operating constraints are as follows: $\hspace{1em}\{\begin{array}{c}E\left(t+1\right)=E\left(t\right)-p^B\left(t\right).Math.\Delta T.Math.{\alpha }_c\\ E\left(t+1\right)=E\left(t\right)-p^B\left(t\right).Math.\Delta T\text{/}{\alpha }_d\\ \underset{\_}{E}\le E\left(t\right)\le \stackrel{\_}{E}\\ .Math.p^B\left(t\right).Math.\le \stackrel{\_}{p^B}\end{array}$ in the formulae, E(t) and E(t+1) are energy storage levels of the storage battery at the time t and a time t+1 respectively; Ē and Ē are upper and lower boundaries of the capacity of the storage battery respectively; ΔT is a time interval from the time t to the time t+1; α.sub.c and α.sub.d are charging and discharging efficiencies of the storage battery respectively; p.sup.B(t) is charging/discharging power of the storage battery at the time t, which is negative during charging and positive during discharging; and p.sup.B is an upper limit value of the charging/discharging power of the storage battery; the exchange electric power constraints for the micro-grids and the main grid are as follows: $\hspace{1em}\{\begin{array}{c}.Math.p^G\left(t\right).Math.\hspace{1em}\hspace{1em}\le \stackrel{\_}{p^G}\\ .Math.p^G\left(t\right)-p^G\left(t-1\right).Math.\le \delta \end{array}$ in the formulae, p.sup.G(t) and p.sup.G(t−1) are exchange electric power between the micro-grids and the main grid at the times t and t−1 respectively, which is positive during purchasing of electricity from the main grid and negative during selling of electricity to the main grid; p.sup.G is an upper limit value of the exchange electric power between the micro-grids and the main grid, and δ is an upper fluctuation limit of the exchange electric power between the micro-grids and the main grid; the electric power output constraints for the combined heat and power equipment are as follows: $\hspace{1em}\{\begin{array}{c}0\le p_i^c\left(t\right)\hspace{1em}\le \stackrel{\_}{p_i^c},i=1,2,\mathrm{.Math.},N\\ \underset{\_}{v_i}.Math.H_i\left(t\right)\le p_i^c\left(t\right)\le \stackrel{\_}{v_i}.Math.H_i\left(t\right)\end{array}$ in the formulae, p.sub.i.sup.c(t) is output electric power of the i th combined heat and power equipment at the time t; p.sub.i.sup.c is an upper limit of the output electric power of the ith combined heat and power equipment H.sub.i(t) is thermal power required to be supplied by the ith combined heat and power equipment at the time t; v.sub.i and v.sub.i are lower and upper limits of electric-thermal power conversion efficiency of the ith combined heat and power equipment.

2. The method for real-time scheduling of multi-energy complementary micro-grids based on the Rollout algorithm according to claim 1, wherein the establishing a target function of the real-time scheduling for the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs in Step 2 specifically comprises the following sub-steps of: first setting up an operating cost function of the micro-grid system at a single scheduling interval with the goal of minimum operating cost of the micro-grid system at the single scheduling interval, and then establishing a target function of the real-time scheduling with the goal of the minimum operating cost of the micro-grid system in the moving-horizon Markov decision cycle; the operating cost function of the micro-grid system at the single scheduling interval is as follows: $\begin{array}{cc}{c}_t\left(X\left(t\right),A\left(t\right)\right)=\lambda \left(t\right).Math.p^G\left(t\right).Math.\Delta T+c.Math.\underset{i=1}{\overset{N}{.Math.}}{F}_i^c\left(t\right)& \\ \mathrm{Wherein}& \\ \{\begin{array}{c}{F}_i^c\left(t\right)=a_i.Math.p_i^c\left(t\right)+b_i\\ X\left(t\right)=\left[E\left(t\right),p^G\left(t-1\right),p^w\left(t\right)\right]\\ A\left(t\right)=\left[p_i^c\left(t\right),p^G\left(t\right),p^B\left(t\right)\right]\end{array}& \end{array}$ in the formulae, X(t) is a state variable of the micro-grid system at the time t; A(t) is a control variable of the micro-grid system at the time t; c.sub.t(X(t),A(t)) is a function of system operating cost at the single scheduling interval; λ(t) is a grid electricity price at the time t; c is a price of natural gas; F.sub.i.sup.c(t) is a linear function between a gas consumption and an electric output of the i th combined heat and power equipment; and a.sub.i and b.sub.i are coefficients of the linear function between the gas consumption and the electric output of the i th combined heat and power equipment; the target function of the real-time scheduling is as follows: $\min J_t\left(X\left(t\right),A\left(t\right)\right)=c_t\left(X\left(t\right),A\left(t\right)\right)+E\left[\underset{t_1=t+1}{\overset{t+T-1}{.Math.}}{c}_{t_1}\left(X\left(t_1\right),A\left(t_1\right)\right)\right]$ in the formula, J.sub.i(X(t),A(t)) is a function of the operating cost of the micro-grid system in the moving-horizon Markov decision cycle.

3. The method for real-time scheduling of multi-energy complementary micro-grids based on the Rollout algorithm according to claim 1, wherein Step 3 specifically comprises the following sub-steps of: dividing a complete scheduling cycle into a plurality of scheduling intervals, finding a solution specific to a scheduling optimization problem in each of the scheduling intervals based on the greedy algorithm respectively, and finally synthesizing locally optimal solutions to respective scheduling intervals into one basic feasible solution across the complete scheduling interval.

4. The method for real-time scheduling of multi-energy complementary micro-grids based on the Rollout algorithm according to claim 2, characterized in that the finding a solution specific to a scheduling optimization problem in each of the scheduling intervals based on the greedy algorithm respectively in Step 3 specifically comprises the following sub-steps of: (1) according to the operating cost function of the micro-grid system in a single scheduling interval in Step 2, listing the target function and the constraint conditions as follows: $\min c_t\left(X\left(t\right),A\left(t\right)\right)=\lambda \left(t\right).Math.p^G\left(t\right).Math.\Delta T+\underset{i=1}{\overset{N}{.Math.}}c.Math.a_i.Math.p_i^c\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{b}_i$ the constraint conditions are as follows: $\begin{array}{cc}\{\begin{array}{c}{p}^G\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{p}_i^c\left(t\right)+p^B\left(t\right)+p^w\left(t\right)=p^D\left(t\right)\\ \underset{\_}{\mathrm{pb}}\left(t\right)\le p^B\left(t\right)\le \stackrel{\_}{\mathrm{pb}}\left(t\right)\\ \underset{\_}{\mathrm{pg}}\left(t\right)\le p^G\left(t\right)\le \stackrel{\_}{\mathrm{pg}}\left(t\right)\\ \underset{\_}{{\mathrm{pc}}_i}\left(t\right)\le p_i^c\left(t\right)\le \stackrel{\_}{{\mathrm{pc}}_i}\left(t\right)i=1,2,\mathrm{.Math.},N\end{array}& \\ \mathrm{wherein}& \\ \{\begin{array}{c}\underset{\_}{\mathrm{pb}}\left(t\right)=\max \left\{\left(E\left(t\right)-\stackrel{\_}{E}\right)\text{/}\left(\Delta T.Math.{\alpha }_c\right),-\stackrel{\_}{p^B}\right\}\\ \stackrel{\_}{\mathrm{pb}}\left(t\right)=\min \left\{\left(E\left(t\right)-\underset{\_}{E}\right).Math.{\alpha }_d\text{/}\Delta T,\stackrel{\_}{p^B}\right\}\\ \underset{\_}{\mathrm{pg}}\left(t\right)=\max \left\{\underset{\_}{p^G},p^G\left(t-1\right)-\delta \right\}\\ \stackrel{\_}{\mathrm{pg}}\left(t\right)=\min \left\{\stackrel{\_}{p^G},p^G\left(t-1\right)+\delta \right\}\\ \underset{\_}{{\mathrm{pc}}_i}\left(t\right)=\max \left\{0,v_i.Math.H_i\left(t\right)\right\}\\ \stackrel{\_}{{\mathrm{pc}}_i}\left(t\right)=\min \left\{\stackrel{\_}{p_i^c},\stackrel{\_}{v_i}.Math.H_i\left(t\right)\right\}\end{array}& \end{array}$ in the formulae, pb(t) and pb(t) are new lower and upper limits of the charging/discharging power of the storage battery during building of the basic feasible solution respectively; pg(t) and pg(t) are new lower and upper limits of the exchange power between the micro-grid and main grid during the building of the basic feasible solution respectively; pc.sub.i (t) and pc.sub.i (t) are new lower and upper limits of the electric output of the i th combined heat and power equipment during the building of the basic feasible solution respectively; (2) to obtain an optimized result in one of the scheduling intervals, ranking N+2 coefficients of the target function in an ascending order to obtain a scheduling sequence for respective decision variables, wherein the N+2 coefficients are λ(t).Math.ΔT, 0, c.Math.a.sub.i (i=1, 2, . . . , N), and the presence of the coefficient 0 is because the target function does not comprise p.sup.B; then forcing a load difference to be d(t)=p.sup.D(t)−p.sup.w(t); and finding a final optimized result according to the value of d(t) and the scheduling sequence of respective decision variables; and (3) finding solutions for other scheduling intervals in turn to obtain one basic feasible solution meeting the micro-grid operating demand finally.

5. The method for real-time scheduling of multi-energy complementary micro-grids based on the Rollout algorithm according to claim 1, wherein Step 4 specifically comprises the following sub-steps of: (1) supposing the basic feasible solution obtained from Step 3 to be π.sub.b=(A.sub.b,t+1, A.sub.b,t+2, . . . , A.sub.b,t+T−1) to obtain an approximate value of the target function for the minimum operating cost from the time t to the time t+T−1 during one real-time scheduling: $\stackrel{\_}{J_t}\left(X\left(t\right),A\left(t\right)\right)=c_t\left(X\left(t\right),A\left(t\right)\right)+E\left[\underset{t_1=t+1}{\overset{t+T-1}{.Math.}}{c}_{t_1}\left(X\left(t_1\right),A_{b,t_1}\right)\right]$ in the formula, J(X(t),A(t)) is a function of the operating cost of the micro-grid system in the moving-horizon Marcov decision cycle as solved with the Rollout algorithm; (2) forcing the approximate value to approach the minimum with the Rollout algorithm: $\stackrel{\_}{a_t}=\arg \underset{A\left(t\right)}{\min }\stackrel{\_}{J_t}\left(X\left(t\right),A\left(t\right)\right)$ (3) according to X(t) in one of the current scheduling intervals and in combination with a generated output value of wind power at the time t, finding X(t+1) in a next scheduling interval, finding the generated output value of the wind power at the time t+1 at the same time, then calculating X(t+2) in the next scheduling interval, and repeating the process until the whole scheduling cycle ends.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 shows a process flow diagram of the present invention.

DETAILED DESCRIPTION

(2) The embodiments of the present invention are further described in detail below with reference to the accompanying drawings:

(3) The present invention provides a method for real-time scheduling of multi-energy complementary micro-grids based on the Rollout algorithm, which not only takes the fluctuations in the new energy outputs into consideration, but also more effectively solves the problem on the moving-horizon scheduling of the multi-energy complementary micro-grids, solving the problem on the moving-horizon Marcov decision model with the Rollout algorithm. According to the method, at first, the moving-horizon Markov decision process model for multi-energy complementary micro-grid real-time scheduling with random new-energy output is set up, and the constraint conditions and the target function for the real-time scheduling are established; then, a complete scheduling cycle is divided into a plurality of scheduling intervals, and one basic feasible solution meeting the constraint conditions for the real-time scheduling is found based on the greedy algorithm; and finally, a solution to the moving-horizon Markov decision model for multi-energy complementary micro-grids is found by using the Rollout algorithm based on the basic feasible solution above.

(4) A method for real-time scheduling of multi-energy complementary micro-grids based on a Rollout algorithm, as shown in FIG. 1, comprises the following steps of:

(5) Step 1, setting up a moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs, and establishing constraint conditions for the real-time scheduling;

(6) where the constraint conditions established for the real-time scheduling in Step 1 comprises: micro-grid electric equilibrium constraints, storage battery operating constraints, exchange electric power constraints for the micro-grids and a main grid, and electric power output constraints for combined heat and power equipment;
the micro-grid electric equilibrium constraints are as follows:

(7) $p^G\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{p}_i^c\left(t\right)+p^B\left(t\right)+p^w\left(t\right)=p^D\left(t\right)$
in the formula, t is a time parameter; p.sup.G(t) is exchange electric power for the micro-grids and the main grid at a time t, which is positive during purchasing of electricity from the main grid and negative during selling of electricity to the main grid; N is the quantity of the combined heat and power equipment; p.sub.i.sup.c(t) is output electric power of the ith combined heat and power equipment at the time t; p.sup.B(1) is charging/discharging power of the storage battery at the time t, which is negative during charging and positive during discharging; p.sup.w(t) is generated output of wind power at the time t; and p.sup.D(t) is an electric load demand at the time t;
the storage battery operating constraints are as follows:

(8) $\hspace{1em}\{\begin{array}{c}E\left(t+1\right)=E\left(t\right)-p^B\left(t\right).Math.\Delta T.Math.{\alpha }_c\\ E\left(t+1\right)=E\left(t\right)-p^B\left(t\right).Math.\Delta T\text{/}{\alpha }_d\\ \underset{\_}{E}\le E\left(t\right)\le \stackrel{\_}{E}\\ .Math.p^B\left(t\right).Math.\le \stackrel{\_}{p^B}\end{array}$
in the formulae, E(t) and E(t+1) are energy storage levels of the storage battery at the time t and a time t+1 respectively; E and Ē are upper and lower boundaries of the capacity of the storage battery respectively; ΔT is a time interval from the time t to the time t+1; a.sub.c and a.sub.d are charging and discharging efficiencies of the storage battery respectively; p.sup.B(t) is charging/discharging power of the storage battery at the time t, which is negative during charging and positive during discharging; and p.sup.B is an upper limit value of the charging/discharging power of the storage battery;
the exchange electric power constraints for the micro-grids and the main grid are as follows:

(9) $\hspace{1em}\{\begin{array}{c}.Math.p^G\left(t\right).Math.\le \stackrel{\_}{p^G}\\ .Math.p^G\left(t\right)-p^G\left(t-1\right).Math.\le \delta \end{array}$
in the formulae, p.sup.G(t) and p.sup.G(t−1) are exchange electric power between the micro-grids and the main grid at the times t and t−1 respectively, which is positive during purchasing of electricity from the main grid and negative during selling of electricity to the main grid; p.sup.G is an upper limit value of the exchange electric power between the micro-grids and the main grid, and δ is an upper fluctuation limit of the exchange electric power between the micro-grids and the main grid;
the electric power output constraints for the combined heat and power equipment are as follows:

(10) $\hspace{1em}\{\begin{array}{c}0\le p_i^c\left(t\right)\le \stackrel{\_}{p_i^c},i=1,2,\mathrm{.Math.},N\\ \underset{\_}{v_i}.Math.H_i\left(t\right)\le p_i^c\left(t\right)\le \stackrel{\_}{v_i}.Math.H_i\left(t\right)\end{array}$
in the formulae, p.sub.i.sup.c(t) is output electric power of the ith combined heat and power equipment at the time t; p.sub.i.sup.c is an upper limit of the output electric power of the ith combined heat and power equipment; H.sub.i(t) is thermal power required to be supplied by the ith combined heat and power equipment at the time t; v.sub.i and v.sub.i are lower and upper limits of electric-thermal power conversion efficiency of the ith combined heat and power equipment.

(11) Step 2, establishing a target function of the real-time scheduling for the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs, with the goal of minimum operating cost of a micro-grid system in a moving-horizon Markov decision cycle;

(12) wherein when in a grid-connected state with the main grid, the multi-energy complementary micro-grids can exchange electricity with the main grid, with energy supply equipment comprising wind driven generators, combined heat and power (CHP) equipment, and storage batteries; and the target function for real-time scheduling is to achieve the minimum operating cost, including system electricity purchasing cost and fuel cost of the CHP equipment, for the micro-grid system.
the establishing a target function of the real-time scheduling for the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids with random new-energy outputs specifically comprises the following sub-steps of: first setting up an operating cost function of the micro-grid system at a single scheduling interval with the goal of minimum operating cost of the micro-grid system at the single scheduling interval, and then establishing a target function of the real-time scheduling with the goal of the minimum operating cost of the micro-grid system in the moving-horizon Markov decision cycle;
the operating cost function of the micro-grid system at the single scheduling interval is as follows:

(13) $\begin{array}{cc}{c}_t\left(X\left(t\right),A\left(t\right)\right)=\lambda \left(t\right).Math.p^G\left(t\right).Math.\Delta T+c.Math.\underset{i=1}{\overset{N}{.Math.}}{F}_i^c\left(t\right)& \\ \mathrm{Wherein}& \\ \{\begin{array}{c}{F}_i^c\left(t\right)=a_i.Math.p_i^c\left(t\right)+b_i\\ X\left(t\right)=\left[E\left(t\right),p^G\left(t-1\right),p^w\left(t\right)\right]\\ A\left(t\right)=\left[p_i^c\left(t\right),p^G\left(t\right),p^B\left(t\right)\right]\end{array}& \end{array}$
in the formulae, X(t) is a state variable of the micro-grid system at the time t; A(t) is a control variable of the micro-grid system at the time t; c.sub.t(X(t),A(t)) is a function of system operating cost at the single scheduling interval; λ(t) is a grid electricity price at the time t; c is a price of natural gas; F.sub.i.sup.c(t) is a linear function between a gas consumption and an electric output of the ith combined heat and power equipment; and a.sub.i and b.sub.i are coefficients of the linear function between the gas consumption and the electric output of the ith combined heat and power equipment;
the target function of the real-time scheduling is as follows:

(14) $\min J_t\left(X\left(t\right),A\left(t\right)\right)=c_t\left(X\left(t\right),A\left(t\right)\right)+E\left[\underset{t_1=t+1}{\overset{t+T-1}{.Math.}}{c}_{t_1}\left(X\left(t_1\right),A\left(t_1\right)\right)\right]$
in the formula, J.sub.t(X(t),A(t)) is a function of the operating cost of the micro-grid system in the moving-horizon Markov decision cycle;

(15) Step 3, dividing a single complete scheduling cycle into a plurality of scheduling intervals, and finding one basic feasible solution meeting the constraint conditions for the real-time scheduling based on a greedy algorithm;

(16) where Step 3 specifically comprises the following sub-steps of: dividing a complete scheduling cycle into a plurality of scheduling intervals, finding a solution specific to a scheduling optimization problem in each of the scheduling intervals based on the greedy algorithm respectively, and finally synthesizing locally optimal solutions to respective scheduling intervals into one basic feasible solution across the complete scheduling interval.
the finding a solution specific to a scheduling optimization problem in each of the scheduling intervals based on the greedy algorithm respectively in Step 3 specifically comprises the following sub-steps of: (1) according to the operating cost function of the micro-grid system in a single scheduling interval in Step 2, listing the target function and the constraint conditions as follows:

(17) $\min c_t\left(X\left(t\right),A\left(t\right)\right)=\lambda \left(t\right).Math.p^G\left(t\right).Math.\Delta T+\underset{i=1}{\overset{N}{.Math.}}c.Math.a_i.Math.p_i^c\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{b}_i$
the constraint conditions are as follows:

(18) $\begin{array}{cc}\{\begin{array}{c}{p}^G\left(t\right)+\underset{i=1}{\overset{N}{.Math.}}{p}_i^c\left(t\right)+p^B\left(t\right)+p^w\left(t\right)=p^D\left(t\right)\\ \underset{\_}{\mathrm{pb}}\left(t\right)\le p^B\left(t\right)\le \stackrel{\_}{\mathrm{pb}}\left(t\right)\\ \underset{\_}{\mathrm{pg}}\left(t\right)\le p^G\left(t\right)\le \stackrel{\_}{\mathrm{pg}}\left(t\right)\\ \underset{\_}{{\mathrm{pc}}_i}\left(t\right)\le p_i^c\left(t\right)\le {\mathrm{pc}}_i\left(t\right)i=1,2,\mathrm{.Math.},N\end{array}& \\ \mathrm{wherein}& \\ \{\begin{array}{c}\underset{\_}{\mathrm{pb}}\left(t\right)=\max \left\{\left(E\left(t\right)-\stackrel{\_}{E}\right)\text{/}\left(\Delta T.Math.{\alpha }_c\right),-\stackrel{\_}{p^B}\right\}\\ \stackrel{\_}{\mathrm{pb}}\left(t\right)=\min \left\{\left(E\left(t\right)-\underset{\_}{E}\right).Math.{\alpha }_d\text{/}\Delta T,\stackrel{\_}{p^B}\right\}\\ \underset{\_}{\mathrm{pg}}\left(t\right)=\max \left\{\underset{\_}{p^G},p^G\left(t-1\right)-\delta \right\}\\ \stackrel{\_}{\mathrm{pg}}\left(t\right)=\min \left\{\stackrel{\_}{p^G},p^G\left(t-1\right)+\delta \right\}\\ \underset{\_}{{\mathrm{pc}}_i}\left(t\right)=\max \left\{0,\underset{\_}{v_i}.Math.H_i\left(t\right)\right\}\\ \stackrel{\_}{{\mathrm{pc}}_i}\left(t\right)=\min \left\{\stackrel{\_}{p_i^c},\stackrel{\_}{v_i}.Math.H_i\left(t\right)\right\}\end{array}& \end{array}$
in the formulae, pb (t)and pb(t) are new lower and upper limits of the charging/discharging power of the storage battery during building of the basic feasible solution respectively; pg(t) and pg(t) are new lower and upper limits of the exchange power between the micro-grid and main grid during the building of the basic feasible solution respectively; pc.sub.i (t) and pc.sub.i (t) are new lower and upper limits of the electric output of the ith combined heat and power equipment during the building of the basic feasible solution respectively; (2) to obtain an optimized result in one of the scheduling intervals, ranking N+2 coefficients of the target function in an ascending order to obtain a scheduling sequence for respective decision variables, wherein the N+2 coefficients are λ(t).Math.ΔT, 0, c.Math.a.sub.i (i=1, 2, . . . , N), the presence of the coefficient 0 is because the target function does not comprise p.sup.B (t); then forcing a load difference to be d(t)=p.sup.D(t)−p.sup.w(t); and finding a final optimized result according to the value of d(t) and the scheduling sequence of respective decision variables; and (3) finding solutions for other scheduling intervals in turn to obtain one basic feasible solution meeting the micro-grid operating demand finally.

(19) Step 4, finding a solution to the moving-horizon Markov decision process model for the real-time scheduling of the multi-energy complementary micro-grids by using the Rollout algorithm based on the basic feasible solution from Step 3.

(20) Step 4 specifically comprises the following sub-steps of: (1) supposing the basic feasible solution obtained from Step 3 to be π.sub.b=(A.sub.b,t+1, A.sub.b,t+2, . . . , A.sub.b,t+T−1) to obtain an approximate value of the target function for the minimum operating cost from the time t to the time t+T−1 during one real-time scheduling:

(21) $\stackrel{\_}{J_t}\left(X\left(t\right),A\left(t\right)\right)=c_t\left(X\left(t\right),A\left(t\right)\right)+E\left[\underset{t_1=t+1}{\overset{t+T-1}{.Math.}}{c}_{t_1}\left(X\left(t_1\right),A_{b,t_1}\right)\right]$
in the formula, J.sub.t (X(t),A(0) is a function of the operating cost of the micro-grid system in the moving-horizon Marcov decision cycle as solved with the Rollout algorithm (2) forcing the approximate value to approach the minimum with the Rollout algorithm:

(22) 0$\stackrel{\_}{a_t}=\arg \underset{A\left(t\right)}{\min }\stackrel{\_}{J_t}\left(X\left(t\right),A\left(t\right)\right)$ (3) according to X(t) in one of the current scheduling intervals and in combination with a generated output value of wind power at the time t, finding X (t+1) in a next scheduling interval, finding the generated output value of the wind power at the time t+1 at the same time, then calculating X(t+2) in the next scheduling interval, and repeating the process until the whole scheduling cycle ends.

(23) It should be noted that the described embodiments of the present invention are for an illustrative purpose rather than a limiting purpose, and the present invention thus includes but not limited to the embodiments described in the Description of Preferred Embodiments. Any other embodiments obtained by those skilled in the art according to the technical solution of the present invention likewise fall within the protection scope of the present invention.