Coordinated optimization peak shaving method for plurality of power supplies based on fluctuation characteristics of renewable energy

Abstract

Disclosed is a coordinated peak shaving optimization method for a plurality of power supplies based on a fluctuation characteristic of a renewable energy source, including the following steps: s1: determining a weekly generated electricity quantity of hydropower based on an available capacity and a storage capacity of an electricity quantity; s2: predicting a renewable energy power generation curve and a load curve of a system weekly; s3: determining a start point of peak shaving of the hydropower based on an external transmission curve, the renewable energy generation curve, and the load curve of the system and a generating capacity of the hydropower; s4: determining a weekly peak shaving demand of the system; and s5: establishing an optimization model with a maximum peak shaving demand. The present disclosure proposes a reasonable arrangement for peak shaving, so as to resolve an accommodation problem caused by large-scale access of a renewable energy source.

Claims

1. A coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy, comprising the following steps: s1: determining a weekly generated electricity quantity of hydropower based on an available capacity and a storage capacity of an electricity quantity; s2: predicting a renewable energy power generation curve and a load curve of a system weekly; s3: determining a start point of peak shaving of the hydropower based on an external transmission curve, the renewable energy power generation curve, and the load curve of the system and a generating capacity of the hydropower; s4: determining a weekly peak shaving demand of the system; s5: establishing an optimization model with a maximum peak shaving demand; s6: calculating a system optimization model based on a branch and bound algorithm to determine a unit combination scheme of the system; and s7: calculating a peak shaving capability of the system and a margin of the peak shaving capability in a peak period of wind power generation, and controlling a thermal power generator to start up or shut down based on the unit combination scheme of the system.

2. The coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy according to claim 1, wherein the determining a start point of peak shaving of the hydropower in step s3 comprises the following steps: a1: determining a basic operating power P.sub.B of the hydropower, and determining, based on a capacity limit of a reservoir, that a limit of a planned weekly electricity quantity generated is [E.sub.min E.sub.max], wherein a difference between E.sub.min and E.sub.max is controlled to be between 2% and 5%; when the P.sub.B is calculated, a generated electricity quantity of the hydropower is calculated according to a formula E.sub.H=(E.sub.min+E.sub.max)/2 (1); a weekly load curve in a regional power grid is denoted as P.sub.l(t), an external transmission power curve is denoted as P.sub.T(t), renewable energy power generation is denoted as P.sub.R(t); and in this case, an equivalent load of the system is: $\begin{array}{cc}{P}_l^E\left(t\right)=P_l\left(t\right)+P_T\left(t\right)-P_R\left(t\right),& \left(2\right)\end{array}$ the generated electricity quantity of the hydropower meets: $\begin{array}{cc}\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}=E_H,& \left(3\right)\end{array}$ and a following constraint is met: $\begin{array}{cc}\max \left\{P_l^E\left(t\right)\right\}-P_B{R}_H^{\max };& \left(4\right)\end{array}$ a2: performing initialization based on historical data and operating experience to make $P_B=P_B^0,$ and setting an error parameter ; a3: testing whether an inequality $\begin{array}{cc}\max \left\{P_l^E\left(t\right)\right\}-P_B{R}_H^{\max }& \left(4\right)\end{array}$ true, and if the inequality does not hold true, increasing P.sub.B until the inequality holds true; a4: calculating an integral $\underset{0}{\overset{T}{}}\left[P_l^F\left(t\right)-P_B\right]dt$ according to a formula $\begin{array}{cc}\stackrel{T}{\underset{0}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}=E_H;\mathrm{and}\mathrm{if}& \left(3\right)\end{array}$ $\begin{array}{cc}.Math.\stackrel{T}{\underset{0}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}-E_H.Math.>,\mathrm{and}& \left(5\right)\end{array}$ $\begin{array}{cc}\stackrel{T}{\underset{0}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}>E_H,& \left(6\right)\end{array}$ setting P.sub.B=P.sub.BP (7), and on the contrary, setting P.sub.B=P.sub.B+P (8); and a5: repeating step a4 until $\begin{array}{cc}.Math.\stackrel{T}{\underset{0}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}-E_H.Math.,& \left(9\right)\end{array}$ and recording a current value of P.sub.B.

3. The coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy according to claim 1, wherein an equation condition of an active power balance satisfied by active power outputs, loads, and active power losses of all generators in a power grid at any time for the peak shaving demand in step s4 is as follows: $\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{Gi}-\underset{i=1}{\overset{m}{.Math.}}{P}_{Li}-P_{}=0& \left(10\right)\end{array}$ wherein when a startup mode of thermal power is determined, in order to ensure a maximum power output of the thermal power and reliable power supply of a load when the hydropower is not started, a power P.sub.B of the start point of the peak shaving of the hydropower is used to replace a total equivalent load $\underset{i=1}{\overset{m}{.Math.}}{P}_{Li},$ wherein i represents a total quantity of generators, m represents a total quantity of loads, and P.sub.Li represents an active power of an i.sup.th load, $\underset{i=1}{\overset{n}{.Math.}}{P}_{Gi}$ represents a sum of power outputs of all the generators in the power grid, and a power output P.sub.Gi of each generator meets upper and lower limit constraints shown in an equation (11): $\begin{array}{cc}{P}_{\mathrm{Gi}}^{\min }{P}_{\mathrm{Gi}}{P}_{\mathrm{Gi}}^{\max }& \left(11\right)\end{array}$ wherein P.sub. represents an active power loss of the power grid; and $P_{\mathrm{Gi}}^{\min }$ and $P_{\mathrm{Gi}}^{\max }$ respectively represent lower and upper limits of a power output of an i.sup.th generator; when an equivalent load of the system is reduced from P.sub.B to P.sub.min in a peak period of renewable energy power generation, the active power balance is written as follows: $\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{Gi}^{}-\left(P_B-P_L\right)-P_{}=0& \left(12\right)\end{array}$ wherein P.sub.L represents a variable of the equivalent load, which is a positive value; and P.sub.Gi represents a changed power of the i.sup.th generator; a following formula is obtained by subtracting the equation (10) from an equation (12): $\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{Gi}^{}-\underset{i=1}{\overset{n}{.Math.}}{P}_{Gi}=-P_L=-R_c;& \left(13\right)\end{array}$ and in this case, P.sub.L represents the variable of the equivalent load, and P.sub.L is always a negative value; R.sub.c is defined as a negative peak shaving capacity of a conventional generator; and in the formula (13), a negative peak shaving capability of the conventional generator is a capability of the conventional generator in reducing an active power output when an equivalent active load in the power grid decreases from the P.sub.B to the P.sub.min, wherein the negative peak shaving capacity is equal to a difference between an output after adjustment and an output before the adjustment, which is always negative.

4. The coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy according to claim 1, wherein the establishing an optimization model with a maximum peak shaving demand in step s5 comprises the following steps: b1: setting a maximum negative peak shaving capacity of a thermal power generator, wherein the maximum negative peak shaving capacity of the thermal power generator depends on a given lower limit parameter $\underset{i=1}{\stackrel{n.Math.}{.Math.}}{P}_{\mathrm{Gimin}}^{}$ of a power output of a generator combination; and in a power grid, there are a plurality of generator combinations meeting a given load requirement, and different generator combinations have different output lower limits; b2: finding, through optimization, a minimum value of the $\underset{i=1}{\stackrel{n.Math.}{.Math.}}{P}_{\mathrm{Gimin}}^{}$ meeting a load demand, so as to provide a negative peak shaving capacity and increase an accommodation proportion of renewable energy; b3: determining that there are N generators participating in the peak shaving in a regional power grid, and a power output of an i.sup.th generator is P.sub.Gi whose value range is defined as [P.sub.G min P.sub.G max], wherein P.sub.G min represents a lower limit of the power output of the i.sup.th generator, and P.sub.G max represents an upper limit of the power output of the i.sup.th generator; defining an N-dimensional vector to represent a generator combination, namely, C=[c.sub.1, c.sub.2, . . . c.sub.i, . . . c.sub.n], wherein C.sub.i represents a state of the i.sup.th generator, and the corresponding generator does not participate in power generation when the state of the i.sup.th generator is 0, or has provided a certain power output to the power grid when the state of the i.sup.th generator is 1; defining a vector P.sub.G min=[P.sub.G1 min, . . . P.sub.Gi min, . . . P.sub.Gn min] and a vector P.sub.G max=[P.sub.G1 max, . . . P.sub.Gi max, . . . P.sub.Gn max] and then obtaining $\begin{array}{cc}{R}_c^{\max }\left(C\right)=P_B+P_{}-C^T{P}_{G\min };& \left(14\right)\end{array}$ and using to represent a value space of C, wherein in this space, a quantity of states of the generator combination is 2.sup.n-1, and a purpose of the generator combination is to search for the generator state vector C in the state space to maximize a peak shaving capacity, so a corresponding objective function is written as a formula $\begin{array}{cc}\max R_c^{\max }\left(C\right)=P_B+P_{}-C^T{P}_{G\min };& \left(15\right)\end{array}$ b4: performing summarization according to step b3 to obtain $\begin{array}{cc}\{\begin{array}{cc}{C}^T{P}_{G\min }{P}_B+P_{}{C}^T{P}_{G\max }& \left(a\right)\\ P_{\mathrm{Gmin}}{P}_{Gi}{P}_{\mathrm{Gmax}}& \left(b\right)\\ c_i\left[0,1\right]& \left(c\right)\end{array};& \left(16\right)\end{array}$ b5: the equivalent load changes within a range of L=[P.sub.min P.sub.B] when the start point P.sub.B of the peak shaving of the hydropower is met, and when the equivalent load of the system continues to decrease, a limit of the negative peak shaving capacity of the system continues to decrease, for ensuring adequacy of a system capacity in the case of a minimum load, namely, P.sub.minC.sup.TP.sub.G minR.sub.th (17), wherein R.sub.th represents a minimum peak shaving demand based on a load and a renewable energy prediction error; and b6: obtaining a final optimization model of a thermal power generator combination based on the peak shaving of the hydropower, namely, $\begin{array}{cc}\{\begin{array}{cc}\max R_c^{\max }\left(C\right)=P_B^{}-C^T{P}_{\mathrm{Gmin}}& \left(a\right)\\ s.t.C^T{P}_{\mathrm{Gmin}}{P}_B+P_{}{C}^T{P}_{\mathrm{Gmax}}& \left(b\right)\\ P_{\mathrm{Gmin}}{P}_{Gi}{P}_{\mathrm{Gmax}}& \left(c\right)\\ P_{\min }-C^T{P}_{\mathrm{Gmin}}{R}_{\mathrm{th}}& \left(d\right)\\ c_i\left[0,1\right]& \left(e\right)\end{array}.& \left(18\right)\end{array}$

5. The coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy according to claim 1, wherein the branch and bound algorithm in step s6 is as follows: C1: placing an original problem P0 into a to-be-resolved problem set L, and setting a target value z*= and a solution variable x*=; C2: determining whether the L is null, wherein if the L is null, the current operation is stopped; or if the L is not null, ILP(k) is selected from the L according to a certain policy, and the ILP(k) is deleted from the L; C3: solving a linear relaxation problem LP(k) of the ILP(k), wherein if the LP(k) has no feasible solution, step C2 is performed; otherwise, z.sub.LP.sub.k is set to represent an objective function value of the LP(k), wherein x.sub.LP.sub.k represents the corresponding solution of the LP(k); and C4: if z.sub.LP.sub.kz*, performing step C2, wherein if x.sub.LP.sub.k does not meet an integer constraint, step C3 is performed; otherwise, z*=z.sub.LP.sub.k and x*=x.sub.LP.sub.k are set, and step C2 is performed.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic flowchart of a coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy.

DETAILED DESCRIPTION

(2) The following clearly and completely describes the technical solutions in the embodiments of the present disclosure with reference to the embodiments of the present disclosure. Apparently, the described embodiments are merely a part rather than all of the embodiments of the present disclosure. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

(3) As shown in FIG. 1, the present disclosure provides a coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy.

(4) The coordinated optimization peak shaving method includes the following steps:

(5) s1: Determine a weekly generated electricity quantity of hydropower based on an available capacity and a storage capacity of an electricity quantity.

(6) s2: Predict a renewable energy power generation curve and a load curve of a system weekly.

(7) s3: Determine a start point of peak shaving of the hydropower based on an external transmission curve, the renewable energy power generation curve, and the load curve of the system and a generating capacity of the hydropower.

(8) s4: Determine a weekly peak shaving demand of the system.

(9) s5: Establish an optimization model with a maximum peak shaving demand.

(10) s6: Calculate a system optimization model based on a branch and bound algorithm, and determine a unit combination scheme of the system.

(11) s7: Calculate a peak shaving capability of the system and a margin of the peak shaving capability in a peak period of wind power generation.

(12) The determining a start point of peak shaving of the hydropower in step s3 includes the following steps: a1: Determine a basic operating power P.sub.B of the hydropower, and determine, based on a capacity limit of a reservoir, that a limit of a planned weekly electricity quantity generated is [E.sub.min E.sub.max], where a difference between E.sub.min and E.sub.max is controlled to be between 2% and 5%; when the P.sub.B is calculated, a generated electricity quantity of the hydropower is calculated according to a formula E.sub.H=(E.sub.min+E.sub.max)/2 (1); a weekly load curve in a regional power grid is denoted as P.sub.l(t), an external transmission power curve is denoted as P.sub.T(t), renewable energy power generation is denoted as P.sub.R(t); and in this case, an equivalent load of the system is:

(13) $\begin{array}{cc}{P}_l^E\left(t\right)=P_l\left(t\right)+P_T\left(t\right)-P_R\left(t\right),& \left(2\right)\end{array}$
the generated electricity quantity of the hydropower meets:

(14) $\begin{array}{cc}\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}=E_H,& \left(3\right)\end{array}$
and a following constraint is met:

(15) $\begin{array}{cc}\max \left\{P_l^E\left(t\right)\right\}-P_B{R}_H^{\max };& \left(4\right)\end{array}$ a2: Perform initialization based on historical data and operating experience to make

(16) $P_B=P_B^0,$
and set an error parameter ; a3: Test whether an inequality

(17) $\begin{array}{cc}\max \left\{P_l^E\left(t\right)\right\}-P_B{R}_H^{\max }& \left(4\right)\end{array}$
holds true, and if the inequality does not hold true, increase P.sub.B until the inequality holds true. a4: Calculate an integral

(18) 0$\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}$
according to a formula

(19) $\begin{array}{cc}\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}=E_H;\mathrm{and}\mathrm{if}& \left(3\right)\end{array}$$\begin{array}{cc}.Math.\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}-E_H.Math.>,\mathrm{and}& \left(5\right)\end{array}$$\begin{array}{cc}\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}>E_H,& \left(6\right)\end{array}$
set P.sub.B=P.sub.BP (7), and on the contrary, set P.sub.B=P.sub.B+P (8). a5: Repeat step a4 until

(20) $\begin{array}{cc}.Math.\underset{0}{\overset{T}{}}\left[P_l^E\left(t\right)-P_B\right]\mathrm{dt}-E_H.Math.,& \left(9\right)\end{array}$
and record a current value of P.sub.B.

(21) An equation condition of an active power balance satisfied by active power outputs, loads, and active power losses of all generators in a power grid at any time for the peak shaving demand in step s4 is as follows:

(22) $\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{\mathrm{Gi}}-\underset{i=1}{\overset{m}{.Math.}}{P}_{\mathrm{Li}}-P_{.Math.}=0& \left(10\right)\end{array}$ where when a startup mode of thermal power is determined, in order to ensure a maximum power output of the thermal power and reliable power supply of a load when the hydropower is not started, a power P.sub.B of the start point of the peak shaving of the hydropower is used to replace a total equivalent load

(23) $\underset{i=1}{\overset{m}{.Math.}}{P}_{\mathrm{Li}},$
where i represents a total quantity of generators, m represents a total quantity of loads, and P.sub.Li represents an active power of an i.sup.th load,

(24) $\underset{i=1}{\overset{n}{.Math.}}{P}_{\mathrm{Gi}}$
represents a sum of power outputs of all the generators in the power grid, and a power output P.sub.Gi of each generator must meet upper and lower limit constraints shown in an equation (11):

(25) $\begin{array}{cc}{P}_{\mathrm{Gi}}^{\min }{P}_{\mathrm{Gi}}{P}_{\mathrm{Gi}}^{\max }& \left(11\right)\end{array}$ where P.sub. represents an active power loss of the power grid;

(26) $P_{Gi}^{\min }$
and

(27) $P_{Gi}^{\max }$
respectively represent lower and upper limits of a power output of an i.sup.th generator; when an equivalent load of the system is reduced from P.sub.B to P.sub.min in a peak period of renewable energy power generation, the active power balance is written as follows:

(28) $\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{\mathrm{Gi}}^{}-\left(P_B-P_L\right)-P_{.Math.}=0& \left(12\right)\end{array}$ where P.sub.L represents a variable of the equivalent load, which is a positive value; P.sub.B represents the power of the start point of the peak shaving of the hydropower, in other words, the basic operating power of the hydropower; and P.sub.Gi represents a changed power of the i.sup.th generator; a following formula is obtained by subtracting the equation (10) from an equation (12):

(29) 0$\begin{array}{cc}\underset{i=1}{\overset{n}{.Math.}}{P}_{\mathrm{Gi}}^{}-\underset{i=1}{\overset{n}{.Math.}}{P}_{\mathrm{Gi}}=-P_L=-R_c;& \left(13\right)\end{array}$
and in this case, P.sub.L represents the variable of the equivalent load, and P.sub.L is always a negative value; R.sub.c is defined as a negative peak shaving capacity of a conventional generator; and in the formula (13), a negative peak shaving capability of the conventional generator is a capability of the conventional generator in reducing an active power output when an equivalent active load in the power grid decreases from the P.sub.B to the P.sub.min, where the negative peak shaving capacity is equal to a difference between an output after adjustment and an output before the adjustment, which is always negative.

(30) The establishing an optimization model with a maximum peak shaving demand in step s5 includes the following steps:

(31) b1: Set a maximum negative peak shaving capacity of a thermal power generator, where the maximum negative peak shaving capacity of the thermal power generator depends on a given lower limit parameter

(32) $\underset{i=1}{\overset{n.Math.}{.Math.}}{P}_{\mathrm{Gimin}}^{}$
of a power output of a generator combination; and in a power grid, there are a plurality of generator combinations meeting a given load requirement, and different generator combinations have different output lower limits.

(33) b2: Find, through optimization, a minimum value of the

(34) $\underset{i=1}{\overset{n.Math.}{.Math.}}{P}_{\mathrm{Gimin}}^{}$
meeting a load demand, so as to provide a negative peak shaving capacity and increase an accommodation proportion of a renewable energy source.

(35) b3: Determine that there are N generators participating in the peak shaving in a regional power grid, and the power output of the i.sup.th generator is P.sub.Gi whose value range is defined as [P.sub.G min P.sub.G max] where P.sub.G min represents the lower limit of the power output of the i.sup.th generator, and P.sub.G max represents the upper limit of the power output of the i.sup.th generator; define an N-dimensional vector to represent a generator combination, namely, C=[c.sub.1, c.sub.2, . . . c.sub.i, . . . c.sub.n], where c.sub.i represents a state of the i.sup.th generator, and the corresponding generator does not participate in power generation when the state of the i.sup.th generator is 0, or has provided a certain power output to the power grid when the state of the i.sup.th generator is 1; define a vector P.sub.G min=[P.sub.G1 min, . . . P.sub.Gi min, . . . P.sub.Gn min] and a vector P.sub.G max=[P.sub.G1 max, . . . P.sub.Gi max, . . . P.sub.Gn max], and then obtain

(36) $\begin{array}{cc}{R}_c^{\max }\left(C\right)=P_B+P_{}-C^T{P}_{G\min };& \left(14\right)\end{array}$
and use to represent a value space of C, where in this space, a quantity of states of the generator combination is 2.sup.n-1, and a purpose of the generator combination is to search for the generator state vector C in the state space to maximize a peak shaving capacity, so a corresponding objective function is written as a formula

(37) $\begin{array}{cc}\max R_c^{\max }\left(C\right)=P_B+P_{}-C^T{P}_{G\min }.& \left(15\right)\end{array}$

(38) b4: Perform summarization according to the b3 to obtain

(39) $\begin{array}{cc}\{\begin{array}{cc}{C}^T{P}_{G\min }{P}_B+P_{}{C}^T{P}_{G\max }& \left(a\right)\\ P_{\mathrm{Gmin}}{P}_{Gi}{P}_{\mathrm{Gmax}}& \left(b\right)\\ c_i\left[0,1\right]& \left(c\right)\end{array}.& \left(16\right)\end{array}$

(40) b5: the equivalent load changes within a range of L=[P.sub.min P.sub.B] min when the start point P.sub.B of the peak shaving of the hydropower is met, and when the equivalent load of the system continues to decrease, a limit of the negative peak shaving capacity of the system continues to decrease, for ensuring adequacy of a system capacity in the case of a minimum load, namely,

(41) $\begin{array}{cc}{P}_{\min }-C^T{P}_{G\min }{R}_{\mathrm{th}},& \left(17\right)\end{array}$
where R.sub.th represents a minimum peak shaving demand based on a load and a renewable energy prediction error.

(42) b6: Obtain a final optimization model of a thermal power generator combination based on the peak shaving of the hydropower, namely,

(43) $\begin{array}{cc}\{\begin{array}{cc}\max R_c^{\max }\left(C\right)=P_B^{}-C^T{P}_{\mathrm{Gmin}}& \left(a\right)\\ s.t.C^T{P}_{\mathrm{Gmin}}{P}_B+P_{}{C}^T{P}_{\mathrm{Gmax}}& \left(b\right)\\ P_{\mathrm{Gmin}}{P}_{Gi}{P}_{\mathrm{Gmax}}& \left(c\right)\\ P_{\min }-C^T{P}_{\mathrm{Gmin}}{R}_{\mathrm{th}}& \left(d\right)\\ c_i\left[0,1\right]& \left(e\right)\end{array}.& \left(18\right)\end{array}$

(44) The branch and bound algorithm in step s6 is as follows:

(45) C1: Place an original problem P0 into a to-be-resolved problem set L, and set a target value z*= and a solution variable x*=.

(46) C2: Determine whether the L is null, where if the L is null, the current operation is stopped; or if the L is not null, ILP(k) is selected from the L according to a certain policy, and the ILP(k) is deleted from the L.

(47) C3: Solve a linear relaxation problem LP(k) of the ILP(k), where if the LP(k) has no feasible solution, step C2 is performed; otherwise, z.sub.LP.sub.k is set to represent an objective function value of the LP(k), wherein x.sub.LP.sub.k represents the corresponding solution of the LP(k).

(48) C4: If z.sub.LP.sub.kz*, perform step C2, where if x.sub.LP.sub.k does not meet an integer constraint, step C3 is performed; otherwise, z*=z.sub.LP.sub.k and z*=x.sub.LP.sub.k are set, and step C2 is performed.

(49) The present disclosure provides a coordinated optimization peak shaving method for a plurality of power supplies based on fluctuation characteristics of renewable energy. A working principle of the coordinated optimization peak shaving method is as follows: determining a weekly generated electricity quantity of hydropower based on an available capacity and a storage capacity of an electricity quantity; predicting a renewable energy power generation curve and a load curve of a system weekly; determining a start point of peak shaving of the hydropower based on an external transmission curve, the renewable energy power generation curve, and the load curve of the system and a generating capacity of the hydropower; determining a weekly peak shaving demand of the system, and establishing an optimization model with a maximum peak shaving demand; calculating a system optimization model based on a branch and bound algorithm, and determining a unit combination scheme of the system; and finally calculating a peak shaving capability of the system and a margin of the peak shaving capability in a peak period of wind power generation. In conclusion, the present disclosure proposes a reasonable arrangement for peak shaving, so as to resolve an accommodation problem caused by large-scale access of renewable energy. Further, in order to meet a peak shaving demand for accommodating the renewable energy, an optimization model of a startup and shutdown combination of a thermal power generator is proposed to meet a peak shaving demand of the system. A weekly startup mode of the thermal power generator can be determined by solving and establishing the model, which can ensure a peak shaving capability of the system while avoiding frequent startup and shutdown of the thermal power generator.

(50) It should be understood that in the description of the present disclosure, the indicated orientation or position relationships are based on the embodiments. They are merely intended to facilitate or simplify description of the present disclosure, rather than to indicate or imply that the mentioned apparatus or component must have a specific orientation and must be constructed and operated in a specific orientation. Therefore, these terms should not be construed as a limitation to the present disclosure.

(51) In the present disclosure, unless otherwise clearly specified, for example, a connection may be a fixed connection, a removable connection, or an integral connection; may be a mechanical connection or an electrical connection; may be a direct connection or an indirect connection by using an intermediate medium; or may be intercommunication between two components or interaction between two components. Unless otherwise clearly specified, a person of ordinary skill in the art can understand a specific meaning of the term in the present disclosure based on a specific situation.

(52) Although the embodiments of the present disclosure have been illustrated and described, it should be understood that those of ordinary skill in the art may make various changes, modifications, replacements and variations to the above embodiments without departing from the principle and spirit of the present disclosure, and the scope of the present disclosure is limited by the appended claims and their legal equivalents.