PEAK SHAVING CONTROL METHOD FOR EMERGENT SOURCE-GRID COORDINATION IN CASE OF FAULTY SENDING-END POWER GRID

Abstract

A peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid. The method includes: S1: evaluating dispatchability of a cluster virtual wind power unit; S2: developing a method for calculating a dispatchability index of the cluster virtual wind power unit; S3: analyzing a source-load peak-shaving resource strategy; and S4: distributing a control strategy for tie-line peak shaving. The present disclosure has the following beneficial effects: In the present disclosure, real-time dispatchability of wind power participating in real-time power balance is first analyzed, specific evaluation indexes and calculation methods are provided, and calculation examples are given for verification. Then, an optimized real-time dispatch strategy is provided based on demand-side response resources, and DC and AC tie-lines are coordinated for operation. The peak shaving control method for emergent source-grid coordination when a sending-end power grid is faulty can ensure normal operation.

Claims

1. A peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid, comprising the following steps: S1: evaluating dispatchability of a cluster virtual wind power unit; S2: developing a method for calculating a dispatchability index of the cluster virtual wind power unit; S3: analyzing a source-load peak-shaving resource strategy; and S4: distributing a control strategy for tie-line peak shaving.

2. The peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid according to claim 1, wherein the evaluation method in step S1 comprises: a1: conducting statistical analysis on probability distributions of wind power characteristics and a time-varying characteristic of a correlation coefficient, wherein the wind power characteristics concerned comprise fluctuation and uncertainty, and fluctuation ΔP.sup.VPG(t) and uncertainty ε.sup.VPG(t) of a wind cluster virtual power generation unit (WVPG) may be expressed as: $\begin{array}{cc}\Delta P^{VPG}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}\Delta P_i^W\left(t\right)& \left(1\right)\\ {\mathrm{.Math.}}^{VPG}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}{\mathrm{.Math.}}_i^W\left(t\right)& \left(2\right)\end{array}$ wherein n.sub.k denotes the number of wind farms in the WVPG, t denotes time, ΔP.sub.i.sup.W(t) denotes fluctuating power of an ith wind farm, and ε.sub.i.sup.W(t) denotes uncertainty of the ith wind farm; a2: establishing a model set based on a marginal probability distribution and an offline condition of a rank correlation coefficient matrix, comprising the following steps: a2.1: establishing an offline model set of conditional probability distributions of fluctuation ΔP.sub.i.sup.W(t) and uncertainty ε.sub.i.sup.W(t) of each wind farm in the WVPG; a2.2: establishing an offline model set of conditional rank correlation coefficient matrixes .sup.VPG, .sub.ΔP.sup.VPG of a cumulative probability distribution function of fluctuation ΔP.sub.i.sup.W(t) and uncertainty ε.sub.i.sup.W(t) among a plurality of wind farms in the WVPG; and a2.3: establishing in a rolling manner, by using a Monte Carlo simulation method based on a genetic algorithm, conditional probability distributions of the WVPG and each wind farm in the WVPG at each moment under a joint constraint of a marginal conditional probability distribution of each wind farm and conditional rank correlation coefficient matrixes of the plurality of wind farms; and a3: conducting online sampling and calculation based on a time-varying probability distribution of the WVPG, comprising the following steps: a3.1: conducting initialization, and randomly generating S n.sub.s×n.sub.k-order matrixes as an initial population of the genetic algorithm, wherein n.sub.s denotes a sampling scale, n.sub.k denotes the number of wind farms in the WVPG, and matrix elements are random numbers evenly distributed within [0,1]; a3.2: selecting a rank correlation coefficient matrix from the offline model set; first predicting a wind direction based on wind farms i and j at moment t, and calculating a wind direction difference coefficient; setting (,) as a data pair comprising elements in ith and jth columns in an mth random matrix, wherein m∈(1, 2, . . . , S); grouping data pairs in (,) into several subsets; setting .sub.ε,ij,c as a rank correlation coefficient calculated based on elements in a cth subset of wind farms i and j in the mth random matrix; setting .sub.εij,c as a rank correlation coefficient with the same value as condition set .sub.ε,ij,c in the offline model set, wherein an objective function of an mth matrix A.sub.m may be set as: $\begin{array}{cc}G\left(A_m\right)=\underset{i=1}{\overset{n_k}{.Math.}}\underset{j=1\&i\ne j}{\overset{n_k}{.Math.}}\underset{c=1}{\overset{n_c}{.Math.}}{\left({\widetilde{ℛ}}_{\mathrm{.Math.},\mathrm{ij},c}-{ℛ}_{\mathrm{.Math.},\mathrm{ij},c}\right)}^2& \left(3\right)\end{array}$ wherein n.sub.c denotes the number of subsets; a3.3: solving formula (3) by using the genetic algorithm, wherein elements in the S random matrixes are changed through genetic and mutation operations to allow a rank correlation coefficient matrix of an optimal matrix to approximate to target value .sub.ε,ij,c, that is, to minimize G(A.sub.m) in formula (3), wherein the optimal matrix allows G(A.sub.m) to be minimized in each iteration; a3.4: conducting convergence judgment; setting A.sub.k* as an optimal matrix of a kth generation of the genetic algorithm and A.sub.0* as an initial optimal matrix, wherein in this case, an improvement of the optimal matrix of the kth generation relative to an initial condition is G(A.sub.k*)−G(A.sub.0*); in the kth generation, if |G(A.sub.k*)−G(A.sub.k-1*)|≤σ|G(A.sub.k*)−G(A.sub.0*)|, genetic algorithm iteration stops, wherein σ denotes a convergence coefficient of the algorithm; otherwise, the process goes back to a3.2; a3.5: selecting cumulative distribution function (t) corresponding to wind farm i from the offline model set based on a value of an uncertainty condition set of wind farm i at moment t; setting a.sub.i,g* g=1, 2, . . . , n.sub.s as a gth element in an ith column of the optimal matrix in a3.4, wherein in this case, a gth sample of wind farm i is {tilde over (e)}.sub.i,g.sup.W(t)=(a.sub.i,g*) and a gth sample of the WVPG is ${\tilde{e}}_g^{\mathrm{VPG}}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}{\tilde{e}}_{i,g}^W\left(t\right);$ a3.6: checking sample time correlation; and a3.7: collecting statistics about frequency distributions of uncertainty of each wind farm and the WVPG at moment t based on samples, which is an estimate of a discrete probability distribution when the sampling scale n.sub.s is large enough, wherein in this way, discrete estimates .sub.ε|i,j.sup.W and of time-varying probability distributions of uncertainty of wind farm i and the WVPG at moment t can be obtained.

3. The peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid according to claim 1, wherein the calculation method in step S2 comprises: b1: calculating available power, wherein the available power is a possible value of actual output of the WVPG at moment t; available power of a thermal power unit is any point value between minimum stable combustion output and installed capacity, and a range of the available power remains unchanged at different moments when the thermal power unit operates online; available power of the WVPG complies with a probability distribution, and a probability distribution of the available power at different moments varies with a predicted power value and an uncertainty probability distribution; the available power of the WVPG at moment t is defined as formula (4):
P.sub.a.sup.VPG(t)=p.sub.f.sup.VPG(t)+ε.sup.VPG(t)  (4) wherein P.sub.a.sup.VPG(t) denotes the available power (in the unit of MW) of the WVPG at moment t and is a random variable, and a corresponding time-varying probability distribution is ; p.sub.f.sup.VPG(t) denotes a predicted power value (MW) of the WVPG at moment t; ε.sup.VPG(t) denotes a random variable (MW) of output uncertainty of the WVPG at moment t, and a corresponding time-varying probability distribution is ; the probability distribution of the available power is calculated as follows: (1) conducting a translation operation on a multi-state probability distribution of uncertainty of the WVPG at moment t, wherein a translation length is a predicted value p.sub.f.sup.VPG(t); respectively setting {tilde over (ε)}.sub.i,{tilde over (ε)}.sub.i e as values of an ith state of before and after the translation, and {tilde over (P)}.sub.ε,i,{tilde over (P)}.sub.ε,i as corresponding state probabilities, wherein in this case, $\begin{array}{cc}G\left(A_m\right)=\underset{i=1}{\overset{n_k}{.Math.}}\underset{j=1\&i\ne j}{\overset{n_k}{.Math.}}\underset{c=1}{\overset{n_c}{.Math.}}{\left({\widetilde{ℛ}}_{\mathrm{.Math.},\mathrm{ij},c}-{ℛ}_{\mathrm{.Math.},\mathrm{ij},c}\right)}^2;& \left(3\right)\end{array}$ (2) conducting downward truncation on the probability distribution after the translation from a point valued 0 based on a fact that the minimum available power of the WVPG is 0, wherein if values of state 1 to state a of the probability distribution after the translation is less than 0, probability {tilde over (P)}.sub.ε0″ corresponding to state 0 after the downward truncation is: $\begin{array}{cc}{\tilde{P}}_{\mathrm{.Math.}0}^{″}={\tilde{P}}_{\mathrm{.Math.}0}^{\prime }+\underset{i=1}{\overset{a}{.Math.}}{\tilde{P}}_{\mathrm{.Math.},i}^{\prime }& \left(5\right)\end{array}$ wherein {tilde over (P)}.sub.ε0′ denotes a probability of a multi-state distribution valued 0 before the truncation operation is performed at the point valued 0; if the distribution has no state valued 0 before the truncation, {tilde over (P)}.sub.ε0′=0; b2: calculating available power limits, wherein the available power limits are the maximum and minimum available power of the WVPG at moment t; an upper limit of output of the thermal power unit is the installed capacity, and a lower limit of the output is the minimum stable combustion output, and the lower and upper limits are fixed values at all moments; the available power limits of the WVPG are upper and lower boundary values of the probability distribution at a certain confidence level, and vary with the probability distribution of the available power; the upper limit of the available power of the WVPG at moment t is defined as formula (6) based on the value-at-risk theory:
p.sub.max.sup.VPG=min[a:P{P.sub.a.sup.VPG(t)≤a}≥β]  (6) wherein p.sub.max.sup.VPG(t) denotes the upper limit (MW) of the available power of the WVPG at moment t; β denotes a confidence probability; p.sub.max.sup.VPG(t) means that it is assured with a probability higher than β that the available power of the WVPG at moment t is less than or equal to a, and a minimum value of a is the upper limit of output; the lower limit of the available power of the WVPG at moment t is defined as formula (7):
p.sub.min,k.sup.VPG(t)=max[b:P{p.sub.b,k.sup.VPG(t)≥b}≥β]  (7) Wherein p.sub.min,k.sup.VPG(t) denotes the lower limit (MW) of the available power of the WVPG at moment t; p.sub.min,k.sup.VPG(t) means that it is assured with a probability higher than β that the available power of the WVPG at moment t is greater than or equal to b, and a maximum value of b is the lower limit of the available power; limit indexes of the available power are calculated as follows: separately setting {tilde over (p)}.sub.a|t,i, {tilde over (P)}.sub.a|t,i as an ith state value and state probability of multi-state probability distribution , wherein using the upper limit of the available power as an example, if a cumulative probability corresponding to a bth state value of satisfies formula (8): $\begin{array}{cc}\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a.Math.t,i}\le \beta <\underset{i=1}{\overset{b+1}{.Math.}}{\tilde{P}}_{a.Math.t,i}& \left(8\right)\end{array}$ the upper limit of the available power is defined as formula (9): $\begin{array}{cc}{p}_{\max }^{\mathrm{VPG}}\left(t\right)={\tilde{p}}_{a.Math.t,b}+\left({\tilde{p}}_{a.Math.t,b+1}-{\tilde{p}}_{a.Math.t,b}\right)\frac{\left(\beta -\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a.Math.t,i}\right)}{{\tilde{P}}_{a.Math.t,b+1}}& \left(9\right)\end{array}$ wherein {tilde over (p)}.sub.a|t,b+1−{tilde over (p)}.sub.a|t,b denotes a difference between values of a (b+1)th state and the bth state, and $\left(\beta -\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a.Math.t,i}\right)/{\tilde{P}}_{a.Math.t,b+1}$ denotes a ratio of a difference between the confidence probability β and a cumulative probability of the bth state to a difference between a cumulative probability of the (b+1)th state and that of the bth state; b3: calculating an output adjustment range, wherein the output adjustment range includes output up/down-adjustment limits of the WVPG at moment t; an output adjustment range of the thermal power unit is determined by the upper and lower limits of the available power; an output adjustment range of the WVPG is determined by both the upper and lower limits of the available power and a control method; an output up-adjustment limit p.sub.r,max.sup.VPG(t) (MW) of the WVPG at moment t is equal to the upper limit p.sub.max.sup.VPG(t) of the available power; the output down-adjustment limit is determined in the following two cases: when generator tripping is allowed, the output of the WVPG can be adjusted down to 0, and in this case, p.sub.r,min.sup.VPG(t)=0; when only pitch angle adjustment is allowed, a power down-adjustment limit (MW) of the wind power unit is set to be k % of the available power, and in this case, the output down-adjustment limit may be expressed as:
p.sub.r,min.sup.VPG(t)=k % E(P.sub.a.sup.VPG(t))  (10); and b4: calculating a ramp rate, wherein a ramp rate of the thermal power unit is a fixed value, whereas a ramp rate of the WVPG is determined by output fluctuation of the WVPG and is an upper/lower boundary value of fluctuation probability distribution with a certain confidence probability; in addition, when determining a ramp rate limit of the WVPG in a dispatch model, it is further necessary to consider a restriction imposed by a dispatch procedure on a wind power fluctuation range allowed per unit time; first, a ramp-up rate limit r.sub.u1.sup.VPG(t) (MW/min) of the WVPG determined by output fluctuation at moment t is defined as shown in formula (11):
r.sub.u1.sup.VPG(t)=[a:P{ΔP.sup.VPG(t)≤a}≥β]  (11) wherein ΔP.sup.VPG(t) denotes a random variable of output fluctuation of the WVPG at moment t; β denotes the confidence probability; r.sub.u1.sup.VPG(t) denotes the ramp-up rate limit of the WVPG determined by fluctuation at moment t, and is a positive value; formula (11) means that it is assured with a probability higher than β that a fluctuation value is less than or equal to a, and a minimum value of a is the ramp-up rate limit determined by fluctuation; a ramp-down rate limit r.sub.d1.sup.VPG(t) (MW/min) of the WVPG determined by output fluctuation at moment t is defined as formula (12):
r.sub.d1.sup.VPG=max[b: P{ΔP.sup.VPG(t)≥b}≥β]  (12) wherein r.sub.d1.sup.VPG(t) denotes the ramp-down rate limit of the WVPG determined by output fluctuation at moment t, and is a negative value; formula (12) means that it is assured with a probability higher than β that a fluctuation value is greater than or equal to b, and a maximum value of b is the ramp-down rate limit determined by fluctuation; considering both the ramp rate limits determined by output fluctuation and the restriction imposed by the dispatch procedure on the wind power fluctuation range allowed per unit time, the ramp-up and ramp-down rate limits of the WVPG in the dispatch model may be determined respectively by using formula (13) and formula (14); the ramp-up rate limit of the WVPG is calculated as follows:
r.sub.u.sup.VPG(t)=min[r.sub.u1.sup.VPG(t),r.sub.u2.sup.VPG(t)]  (13) wherein r.sub.u2.sup.VPG(t) denotes the ramp-up rate limit (MW/min) of the WVPG determined by the dispatch procedure at moment t, and r.sub.u.sup.VPG(t) denotes the ramp-up rate limit (MW/min) of the WVPG at moment t; when r.sub.u1.sup.VPG(t)<r.sub.u2.sup.VPG(t), indicates that the ramp-up rate limit determined by fluctuation is lower than the restriction imposed by the dispatch procedure from moment t−1 to moment t; the ramp-up rate limit needs to be set to r.sub.u.sup.VPG(t)=r.sub.u1.sup.VPG(t) to avoid overestimating a power rise rate; when r.sub.u1.sup.VPG(t)>r.sub.u2.sup.VPG(t), it indicates that the ramp-up rate limit determined by fluctuation is higher than the restriction imposed by the dispatch procedure from moment t−1 to moment t; the ramp-up rate limit needs to be set to r.sub.u.sup.VPG(t)=r.sub.u2.sup.VPG(t) to prevent power from rising too fast; the ramp-down rate limit of the WVPG is calculated as follows:
r.sub.d.sup.VPG(t)=min[r.sub.d1.sup.VPG(t),r.sub.d2.sup.VPG(t)]  (14) wherein r.sub.d2.sup.VPG(t) denotes the ramp-down rate limit (MW/min) of the WVPG determined by the dispatch procedure at moment t, and is a negative value; r.sub.d.sup.VPG(t) denotes the ramp-down rate limit (MW/min) of the WVPG at moment t, and is a negative value; if r.sub.d1.sup.VPG(t)<r.sub.d2.sup.VPG(t), it indicates that output determined by fluctuation drops too fast from moment t−1 to moment t, and the ramp-down rate limit needs to be set to r.sub.d.sup.VPG(t)=r.sub.d1.sup.VPG(t) to avoid underestimating the ramp-down rate; if r.sub.d.sup.VPG(t)>r.sub.d2.sup.VPG(t), it indicates that a downward fluctuation rate of output is less than a requirement of the dispatch procedure from moment t−1 to moment t; the ramp-down rate may be increased through a control method such as wind curtailment if necessary, and the ramp-down rate limit is set to r.sub.d.sup.VPG(t)=r.sub.d2.sup.VPG(t); and a ramp rate index is calculated as follows: conducting quantile operation on the multi-state probability distribution of fluctuation of the WVPG at moment t to obtain the ramp-up and ramp-down rate limits r.sub.u1.sup.VPG(t),r.sub.d1.sup.VPG(t) determined by output fluctuation; and then obtaining values of the ramp rate limits of the WVPG at the current moment by using formula (13) and formula (14).

4. The peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid according to claim 1, wherein the analysis method in step S3 comprises: c1. variables c1.1. integer variables 1) v(i,t): a variable that denotes an on/off state of unit i at moment t, wherein v(i,t)=1 indicates that unit i is on at moment t, and v(i,t)=0 indicates that unit i is off at moment t; 2) y(i,t): a variable that denotes a conversion state of unit i at moment t, wherein y(i,t)=1 indicates that unit i changes from the off state to the on state at moment t, and y(i,t)=0 indicates that the unit is in another state; 3) z(i,t): a variable that denotes the conversion state of unit i at moment t, wherein z(i,t)=1 indicates that unit i changes from the on state to the off state at moment t, and z(i,t)=0 indicates that the unit is in another state; 4) ω(i,t, j): a variable that denotes a downtime state of the unit, wherein ω(i,t,j)=1 indicates that unit i has been off for j hours at moment t, and ω(i,t,j)=0 indicates that the unit is in another state; and Tr(i,t): a variable valued 0 or 1 added to denote the startup time of a translatable load, wherein Tr(i,t)=1 indicates that translatable load i starts to operate at moment t, and Tr(i,t)=0 indicates that translatable load i does not start to operate at moment t; c1.2. continuous variable p(i,t): output of unit i at moment t; c2. objective function $\begin{array}{cc}\underset{p\left(i,t\right),y\left(i,t\right)}{\mathrm{Min}}\underset{t=1}{\overset{\mathrm{TT}}{.Math.}}\underset{i=1}{\overset{\mathrm{NG}}{.Math.}}\left[c\left(i,t\right)+b\left(i,t\right)\right]& \left(15\right)\end{array}$
wherein c(i,t)=A(i)×p(i,t).sup.2+B(i)×p(i,t)+C(i)×v(i,t)  (16) wherein TT denotes total operation duration; NG denotes the number of units; A(i), B(i), C(i) respectively denote a quadratic coefficient, a linear coefficient, and a constant term in a quadratic function of power generation costs and output power of unit i; c(i,t) denotes power generation costs of unit i at moment t; b(i,t) denotes startup costs of unit i at moment t; calculation of b(i,t) is related to downtime of the unit before startup, and detailed calculation is described in constraints; shutdown costs is usually not considered because it has little impact on total costs; c3. constraints c3.1. power balance constraint $\begin{array}{cc}\underset{i}{\overset{NG}{.Math.}}p\left(i,t\right)=D_0\left(t\right)+\underset{i=1}{\overset{NJ}{.Math.}}PD\left(i,t\right)\forall t\in T& \left(17\right)\end{array}$ wherein PD(i,t) denotes a load value of translatable load i at moment t; D.sub.0(t) denotes a fixed system load at moment t; NI denotes the number of translatable loads; the sum of the fixed loads and the translatable load after translation matches the total output of the unit based on this constraint; c3.2. constraint on upper and lower output limits of the unit
p(i,t)×v(i,t)≤p(i,t)t)×v(i,t)∀i∈G,t∈T  (18) wherein p(i,t) denotes the lower output limit of unit i at moment t; p(i,t) denotes the upper output limit of unit i at moment t; G denotes a set of units; this constraint ensures that the units can operate in an output range for normal operation; c3.3. constraints on startup/shutdown auxiliary variable
y(i,t)+z(i,t)≤1 ∀i∈G,t∈T  (19)
y(i,t)−z(i,t)=v(i,t)−v(i,t−1)∀i∈G,t∈T  (20) wherein formula (19) imposes a restriction that prevents a unit from being started and shut down at the same time, and formula (20) restricts consistency between startup and shutdown actions and an operating state of the unit; c3.4. constraints on minimum runtime $\begin{array}{cc}\underset{k=1}{\overset{G\left(i\right)}{.Math.}}\left[1-y\left(i,k\right)\right]=0& \left(21\right)\\ \underset{k=t}{\overset{t+T_{\mathrm{up}}\left(i\right)-1}{.Math.}}v\left(i,k\right)\ge T_{up}\left(i\right)\times y\left(i,t\right)\forall t=G\left(i\right)+1LTT-T_{up}\left(i\right)+1& \left(22\right)\\ \underset{k=t}{\overset{TT}{.Math.}}\left[v\left(i,k\right)-y\left(i,t\right)\right]\ge 0\forall t=TT-T_{up}\left(i\right)+2\mathrm{LTT}& \left(23\right)\end{array}$ wherein G(i)=Min[TT,(T.sub.up(i)−T.sub.0(i))×U.sub.0(i)] denotes minimum runtime initially required by unit i after operation starts; T.sub.up(i) denotes minimum startup runtime of unit i; U.sub.0(i) denotes an initial operating state of unit i, and if U.sub.0(i)=1, unit i is initially in operation, and if U.sub.0(i)=0, unit i is initially in shutdown state; T.sub.0(i) denotes initial runtime or downtime of unit i, and is a positive number when it denotes the initial runtime and a negative number when it denotes the initial downtime; impact of an operating state at the end of a previous day on start of operation on a new day is mainly considered in formulas (22) and (23); if the unit is initially in operation, formula (22) indicates that the unit requires at least a runtime constraint; formula (23) denotes a minimum runtime constraint in last T.sub.up(i)−1 hours of runtime; if the unit starts in this time period, it needs to operate until the end of the day and cannot be shut down again; c3.5. constraints on minimum downtime $\begin{array}{cc}\underset{k=1}{\overset{L\left(i\right)}{.Math.}}v\left(i,k\right)=0& \left(24\right)\\ \underset{k=t}{\overset{t+T_{\mathrm{down}}\left(i\right)-1}{.Math.}}\left[1-v\left(i,k\right)\right]\ge T_{down}^{}\left(i\right)\times z\left(i,t\right)\forall t=L\left(i\right)+1\mathrm{.Math.}\mathrm{TT}-T_{down}\left(i\right)+1& \left(25\right)\\ \underset{k=t}{\overset{TT}{.Math.}}\left[1-v\left(i,k\right)-z\left(i,t\right)\right]\ge 0\forall t=TT-T_{down}\left(i\right)+2\mathrm{.Math.}\mathrm{TT}& \left(26\right)\end{array}$ wherein L(i)=Min[TT,(T.sub.down(i)+T.sub.0(i))×(1−U.sub.0(i))], and it denotes minimum downtime initially required by unit i after operation starts; T.sub.down(i) denotes minimum downtime of unit i; impact of an operating state at the end of a previous day on start of operation on a new day is mainly considered in formulas (25) and (26); if the unit is initially in shutdown state, it indicates that the unit requires at least a downtime constraint; formula (25) denotes a minimum downtime constraint each time unit i is shut down; formula (26) denotes a minimum downtime constraint in last T.sub.down(i)−1 hours of runtime; if the unit is shut down in the time period, it needs to be in shutdown state until the end of the day and cannot be started again; c3.6. ramp rate constraints
p(i,t)−p(i,t−1)≤RATEUP(i)∀i∈G,t=2,3 . . . TT  (27)
p(i,t−1)−p(i,t)≤RATEDN(i)∀i∈G,t=2,3 . . . TT  (28) wherein RATEUP(i) denotes a ramp-up rate of unit i, that is, a rate at which output of unit i increases; RATEDN(i) denotes a ramp-down rate of unit i, that is, a rate at which output of unit i decreases; an output change between two adjacent time periods is limited because the output of the unit is adjusted at a limited rate; and the ramp rate constraint represents this limitation; c3.7. downtime counter constraints
S(i,t)≤S(i,t−1)+1 ∀t∈T
S(i,t)+(TT+1)×v(i,t)≥S(i,t−1)+1 ∀t∈T
S(i,t)−TT×[1−v(i,t)]≤0 ∀t∈T
S(i,t)≥0 ∀t∈T  (29) wherein S(i,t) is a downtime count variable that denotes accumulated consecutive downtime of unit i at moment t; it is necessary to use the constraint to record downtime because startup costs in the model are segmented and related to the downtime of the unit; a function implemented by the four constraints is as follows: when v(i,t)=0, S(i,t)=S(i, t−1)+1, and when v(i,t)=1, S(i,t)=0, thereby implementing accumulation of consecutive downtime; c3.8. constraints on segmented startup costs $\begin{array}{cc}b\left(i,t\right)=\underset{j=1}{\overset{ND}{.Math.}}\mathrm{UP}\cos t\left(i,j\right)\times \omega \left(i,t,j\right)\forall t\in T\underset{j=1}{\overset{ND}{.Math.}}\omega \left(i,t,j\right)=y\left(i,t\right)\forall t\in TS\left(i,t-1\right)=\underset{j=1}{\overset{ND-1}{.Math.}}j\times \omega \left(i,t,j\right)+m\left(i,t\right),\forall t\in Tm\left(i,t\right)\le TT\times \left[\omega \left(i,t,ND\right)-y\left(i,t\right)+1\right]\forall t\in Tm\left(i,t\right)\ge \mathrm{ND}\times \omega \left(i,t,\mathrm{ND}\right)\forall t\in T& \left(30\right)\end{array}$ wherein in these constraints, segmented startup costs b(i,t) are calculated by using a counter of consecutive downtime; UP cos t(i,j) denotes startup costs after unit i is shut down for j hours; longer downtime usually indicates higher startup costs, and the startup costs no longer change after the unit is shut down for a certain time period; ND denotes a time threshold that finally allows the startup costs to stop changing in segmented startup costs calculation; ω(i,t,j) is a variable valued 0 or 1 that denotes a downtime state of the unit, wherein ω(i,t,j)=1 indicates that unit i has been shut down for j hours at moment t, and ω(i,t,j)=0 indicates that the unit is in another state; in this constraint, for each time period of each unit, at most one ω(i,t,j) can be set to 1; in formula (30), dummy variable m(i,t) is used to associate consecutive downtime count S(i,t) with ω(i,t,j); c3.9. constraint on period-based load translation $\begin{array}{cc}PD\left(i,t\right)=\underset{j=1}{\overset{T_0}{.Math.}}\left[Tr\left(i,t-j+1\right)\times D_{\mathrm{tr}}\left(i,j\right)\right]\forall i\in NI,t\in NT\left(i\right)\underset{t=1}{\overset{TT}{.Math.}}Tr\left(i,t\right)=1,\forall i\in NI\underset{t=TT-NT\left(i\right)+1}{\overset{TT}{.Math.}}Tr\left(i,t\right)=0,\forall i\in NI& \left(31\right)\end{array}$ wherein T.sub.0=min[t,NT(i)]; D.sub.tr(i,j) denotes a load value of a jth segment in a load sequence of translatable load i; NT(i) denotes a time length of the sequence of translatable load i; and NI denotes a set of sequence numbers of translatable loads.

5. The peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid according to claim 1, wherein the distribution method in step S4 comprising: d1: establishing a general model of an optimal power distribution problem; d2: simplifying the general model of the optimal power distribution problem; d3: solving the simplified model of the optimal power distribution problem; and d4: conducting simulation analysis on a solved value.

Description

DETAILED DESCRIPTION

[0099] 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.

[0100] The present disclosure provides a peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid. The method includes the following steps:

[0101] S1: Evaluate dispatchability of a cluster virtual wind power unit.

[0102] S2: Develop a method for calculating a dispatchability index of the cluster virtual wind power unit.

[0103] S3: Analyze a source-load peak-shaving resource strategy.

[0104] S4: Distribute a control strategy for tie-line peak shaving.

[0105] The evaluation method in step S1 includes the following steps:

[0106] a1: Conduct statistical analysis on probability distributions of wind power characteristics and a time-varying characteristic of a correlation coefficient. The wind power characteristics concerned include fluctuation and uncertainty. Fluctuation ΔP.sup.VPG(t) and uncertainty ε.sup.VPG(t) of a wind cluster virtual power generation unit (WVPG) may be expressed as:

[00016]$\begin{array}{cc}\Delta P^{VPG}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}\Delta P_i^W\left(t\right)\mathrm{and}& \left(1\right)\\ {\mathrm{.Math.}}^{VPG}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}{E}_i^W\left(t\right).& \left(2\right)\end{array}$

[0107] In formulas (1) and (2), n.sub.k denotes the number of wind farms in the WVPG, t denotes time, Δ.sub.i.sup.W(t) denotes fluctuating power of an ith wind farm, and ε.sub.i.sup.W(t) denotes uncertainty of the ith wind farm.

[0108] a2: Establish a model set based on a marginal probability distribution and an offline condition of a rank correlation coefficient matrix, including the following steps:

[0109] a2.1: Establish an offline model set of conditional probability distributions of fluctuation ΔP.sub.i.sup.W(t) and uncertainty ε.sub.i.sup.W(t) of each wind farm in the WVPG.

[0110] a2.2: Establish an offline model set of conditional rank correlation coefficient matrixes .sub.ε.sup.VPG,.sub.ΔP.sup.VPG of a cumulative probability distribution function of fluctuation ΔP.sub.i.sup.W(t) and uncertainty ε.sub.i.sup.W(t) among a plurality of wind farms in the WVPG.

[0111] a2.3: Establish in a rolling manner, by using a Monte Carlo simulation method based on a genetic algorithm, conditional probability distributions of the WVPG and each wind farm in the WVPG at each moment under a joint constraint of a marginal conditional probability distribution of each wind farm and conditional rank correlation coefficient matrixes of the plurality of wind farms.

[0112] a3: Conduct online sampling and calculation based on a time-varying probability distribution of the WVPG, including the following steps:

[0113] a3.1: Conduct initialization, and randomly generate S n.sub.s×n.sub.k-order matrixes as an initial population of the genetic algorithm, where n.sub.s denotes a sampling scale, n.sub.k denotes the number of wind farms in the WVPG, and matrix elements are random numbers evenly distributed within [0,1].

[0114] a3.2: Select a rank correlation coefficient matrix from the offline model set. First predict a wind direction based on wind farms i and j at moment t. Calculate a wind direction difference coefficient. Set (, ) as a data pair including elements in ith and jth columns in an mth random matrix, where m∈(1,2, . . . , S). Group data pairs in (,) into several subsets. Set as a rank correlation coefficient calculated based on elements in a cth subset of wind farms i and j in the mth random matrix. Set as a rank correlation coefficient with the same value as condition set in the offline model set. An objective function of an mth matrix A.sub.m may be set as:

[00017]$\begin{array}{cc}G\left(A_m\right)=\underset{i=1}{\overset{n_k}{.Math.}}\underset{j=1\&i\ne j}{\overset{n_k}{.Math.}}\underset{c=1}{\overset{n_C}{.Math.}}{\left({\widetilde{ℛ}}_{\mathrm{.Math.},\mathrm{ij},c}-{ℛ}_{\mathrm{.Math.},\mathrm{ij},c}\right)}^2.& \left(3\right)\end{array}$

[0115] In formula (3), n.sub.c denotes the number of subsets.

[0116] a3.3: Solve formula (3) by using the genetic algorithm. Elements in the S random matrixes are changed through genetic and mutation operations to allow a rank correlation coefficient matrix of an optimal matrix (which allows G(A.sub.m) to be minimized in each iteration) to approximate to target value , that is, to minimize in formula (3).

[0117] a3.4: Conduct convergence judgment. Set A.sub.k* as an optimal matrix of a kth generation of the genetic algorithm and A.sub.0* as an initial optimal matrix. In this case, an improvement of the optimal matrix of the kth generation relative to an initial condition is G(A.sub.k*)−G(A.sub.0*). In the kth generation, if |G(A.sub.k*)−G(A.sub.k-1*)|≤σ|G(A.sub.k*)−G(A.sub.0*)|, genetic algorithm iteration stops, where σ denotes a convergence coefficient of the algorithm; otherwise, the process goes back to a3.2.

[0118] a3.5: Select cumulative distribution function (t) corresponding to wind farm i from the offline model set based on a value of an uncertainty condition set of wind farm i at moment t. Set a.sub.i,g* g=1, 2, . . . , n.sub.s as a gth element in an ith column of the optimal matrix in a3.4. In this case, a gth sample of wind farm i is {tilde over (e)}.sub.i,g.sup.W=(a.sub.i,g*), and a gth sample of the WVPG is

[00018]${\tilde{e}}_g^{VPG}\left(t\right)=\underset{i=1}{\overset{n_k}{.Math.}}{\tilde{e}}_{i,g}^W\left(t\right).$

[0119] a3.6: Check sample time correlation.

[0120] a3.7: Collect statistics about frequency distributions of uncertainty of each wind farm and the WVPG at moment t based on samples, which is an estimate of a discrete probability distribution when the sampling scale n.sub.s is large enough. In this way, discrete estimates and of time-varying probability distributions of uncertainty of wind farm i and the WVPG at moment t can be obtained.

[0121] The calculation method in step S2 includes the following steps:

[0122] b1: Calculate available power. The available power is a possible value of actual output of the WVPG at moment t. Available power of a thermal power unit is any point value between minimum stable combustion output and installed capacity, and a range of the available power remains unchanged at different moments when the thermal power unit operates online. Available power of the WVPG complies with a probability distribution, and a probability distribution of the available power at different moments varies with a predicted power value and an uncertainty probability distribution.

[0123] The available power of the WVPG at moment t is defined as formula (4):

P.sub.a.sup.VPG(t)=p.sub.f.sup.VPG(t)+ε.sup.VPG(t)  (4).

[0124] In formula (4), P.sub.a.sup.VPG(t) denotes the available power (in the unit of MW) of the WVPG at moment t and is a random variable, and a corresponding time-varying probability distribution is ; p.sub.f.sup.VPG(t) denotes a predicted power value (MW) of the WVPG at moment t; ε.sup.VPG(t) denotes a random variable (MW) of output uncertainty of the WVPG at moment t, and a corresponding time-varying probability distribution is .

[0125] The probability distribution of the available power is calculated as follows:

[0126] (1) Conduct a translation operation on a multi-state probability distribution of uncertainty of the WVPG at moment t, where a translation length is a predicted value p.sub.f.sup.VPG(t); respectively set {tilde over (ε)}.sub.i,{tilde over (ε)}.sub.i′ as values of an ith state of before and after the translation, and {tilde over (P)}.sub.ε,j,{tilde over (P)}.sub.ε,i′ as corresponding state probabilities. In this case,

[00019]$\begin{array}{cc}G\left(A_m\right)=\underset{i=1}{\overset{n_k}{.Math.}}\underset{j=1\&i\ne j}{\overset{n_k}{.Math.}}\underset{c=1}{\overset{n_C}{.Math.}}{\left({\widetilde{ℛ}}_{\mathrm{.Math.},\mathrm{ij},c}-{ℛ}_{\mathrm{.Math.},\mathrm{ij},c}\right)}^2.& \left(3\right)\end{array}$

[0127] (2) Conduct downward truncation on the probability distribution after the translation from a point valued 0 based on a fact that the minimum available power of the WVPG is 0. If values of state 1 to state a of the probability distribution after the translation is less than 0, probability {tilde over (P)}.sub.ε0″ corresponding to state 0 after the downward truncation is:

[00020]$\begin{array}{cc}{\tilde{P}}_{\mathrm{.Math.}0}^{″}={\tilde{P}}_{\mathrm{.Math.}0}^{\prime }+\underset{i=1}{\overset{a}{.Math.}}{\tilde{P}}_{\mathrm{.Math.},i}^{\prime }& \left(5\right)\end{array}$

[0128] In formula (5), {tilde over (P)}.sub.ε0′ denotes a probability of a multi-state distribution valued 0 before the truncation operation is performed at the point valued 0. If the distribution has no state valued 0 before the truncation, {tilde over (P)}.sub.ε0′=0.

[0129] b2: Calculate available power limits. The available power limits are the maximum and minimum available power of the WVPG at moment t. An upper limit of output of the thermal power unit is the installed capacity, and a lower limit of the output is the minimum stable combustion output. The lower and upper limits are fixed values at all moments. The available power limits of the WVPG are upper and lower boundary values of the probability distribution at a certain confidence level, and vary with the probability distribution of the available power.

[0130] The upper limit of the available power of the WVPG at moment t is defined as formula (6) based on the value-at-risk theory:

p.sub.max.sup.VPG(t)=min[a: P{P.sub.a.sup.VPG(t)≤a}≥β]  (6).

[0131] In formula (6), p.sub.max.sup.VPG(t) denotes the upper limit (MW) of the available power of the WVPG at moment t, and β denotes a confidence probability. p.sub.max.sup.VPG(t) means that it is assured with a probability higher than β that the available power of the WVPG at moment t is less than or equal to a, and a minimum value of a is the upper limit of output.

[0132] The lower limit of the available power of the WVPG at moment t is defined as formula (7):

p.sub.min,k.sup.VPG(t)=max[b: P{p.sub.b,k.sup.VPG(t)≥b}≥β]  (7).

[0133] In formula (7), p.sub.min,k.sup.VPG(t) denotes the lower limit (MW) of the available power of the WVPG at moment t. p.sub.min,k.sup.VPG(t) means that it is assured with a probability higher than β that the available power of the WVPG at moment t is greater than or equal to b, and a maximum value of b is the lower limit of the available power.

[0134] Limit indexes of the available power are calculated as follows:

[0135] Separately set {tilde over (p)}.sub.a|t,i,{tilde over (P)}.sub.a|t,i as an ith state value and state probability of multi-state probability distribution . Using the upper limit of the available power as an example, if a cumulative probability corresponding to a bth state value of satisfies formula (8):

[00021]$\begin{array}{cc}\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a|t,i}\le \beta <\underset{i=1}{\overset{b+1}{.Math.}}{\tilde{P}}_{a|t,i},& \left(8\right)\end{array}$

[0136] the upper limit of the available power is defined as formula (9):

[00022]$\begin{array}{cc}{p}_{\max }^{\mathrm{VPG}}\left(t\right)={\tilde{p}}_{a|t,b}+\left({\tilde{p}}_{a|t,b+1}-{\tilde{p}}_{a|t,b}\right)\frac{\left(\beta -\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a|t,i}\right)}{{\tilde{P}}_{a|t,b+1}}.& \left(9\right)\end{array}$

[0137] In formula (9), {tilde over (p)}.sub.a|t,b+1−{tilde over (p)}.sub.a|t,b denotes a difference between values of a (b+1)th state and the bth state,

[00023]$\left(\beta -\underset{i=1}{\overset{b}{.Math.}}{\tilde{P}}_{a|t,i}\right)/{\tilde{P}}_{a|t,b+1}$

and denotes a ratio of a difference between the confidence probability β and a cumulative probability of the bth state to a difference between a cumulative probability of the (b+1)th state and that of the bth state.

[0138] b3: Calculate an output adjustment range, where the output adjustment range includes output up/down-adjustment limits of the WVPG at moment t. An output adjustment range of the thermal power unit is determined by the upper and lower limits of the available power. An output adjustment range of the WVPG is determined by both the upper and lower limits of the available power and a control method.

[0139] An output up-adjustment limit p.sub.r,max.sup.VPG(t) of the WVPG at moment t is equal to the upper limit p.sub.max.sup.VPG(t) of the available power. The output down-adjustment limit is determined in the following two cases: When generator tripping is allowed, the output of the WVPG can be adjusted down to 0. In this case, p.sub.r,min.sup.VPG(t)=0. When only pitch angle adjustment is allowed, a power down-adjustment limit (MW) of the wind power unit is set to be k % of the available power. In this case, the output down-adjustment limit may be expressed as:

p.sub.r,min.sup.VPG(t)=k % E(P.sub.a.sup.VPG(t))  (10).

[0140] b4: Calculate a ramp rate. A ramp rate of the thermal power unit is a fixed value, whereas a ramp rate of the WVPG is determined by output fluctuation of the WVPG and is an upper/lower boundary value of fluctuation probability distribution with a certain confidence probability. In addition, when determining a ramp rate limit of the WVPG in a dispatch model, it is further necessary to consider a restriction imposed by a dispatch procedure on a wind power fluctuation range allowed per unit time.

[0141] First, a ramp-up rate limit r.sub.u1.sup.VPG(t) (MW/min) of the WVPG determined by output fluctuation at moment t is defined as shown in formula (11):

r.sub.u1.sup.VPG(t)=min[a:P{ΔP.sup.VPG(t)≤a}≥β]  (11).

[0142] In formula (11), ΔP.sup.VPG(t) denotes a random variable of output fluctuation of the WVPG at moment t; β denotes the confidence probability; r.sub.u1.sup.VPG(t) denotes the ramp-up rate limit of the WVPG determined by fluctuation at moment t, and is a positive value. Formula (11) means that it is assured with a probability higher than β that a fluctuation value is less than or equal to a, and a minimum value of a is the ramp-up rate limit determined by fluctuation.

[0143] A ramp-down rate limit rdlPG (t) (MW/min) of the WVPG determined by output fluctuation at moment t is defined as formula (12):

r.sub.d1.sup.VPG(t)=max[b: P{ΔP.sup.VPG(t)≥b}≥β]  (12).

[0144] In formula (12), r.sub.d1.sup.VPG(t) denotes the ramp-down rate limit of the WVPG determined by output fluctuation at moment t, and is a negative value. Formula (12) means that it is assured with a probability higher than β that a fluctuation value is greater than or equal to b, and a maximum value of b is the ramp-down rate limit determined by fluctuation.

[0145] Considering both the ramp rate limits determined by output fluctuation and the restriction imposed by the dispatch procedure on the wind power fluctuation range allowed per unit time, the ramp-up and ramp-down rate limits of the WVPG in the dispatch model may be determined respectively by using formula (13) and formula (14).

[0146] The ramp-up rate limit of the WVPG is calculated as follows:

r.sub.u.sup.VPG(t)=min[r.sub.u1.sup.VPG(t),r.sub.u2.sup.VPG(t)]  (13).

[0147] In formula (13), r.sub.u2.sup.VPG(t) denotes the ramp-up rate limit (MW/min) of the WVPG determined by the dispatch procedure at moment t, and r.sub.u.sup.VPG(t) denotes the ramp-up rate limit (MW/min) of the WVPG at moment t.

[0148] When r.sub.u1.sup.VPG(t)<r.sub.u2.sup.VPG(t), indicates that the ramp-up rate limit determined by fluctuation is lower than the restriction imposed by the dispatch procedure from moment t−1 to moment t. The ramp-up rate limit needs to be set to r.sub.u.sup.VPG(t)=r.sub.u1.sup.VPG(t) to avoid overestimating a power rise rate.

[0149] When r.sub.u1.sup.VPG(t)>r.sub.u2.sup.VPG(t), it indicates that the ramp-up rate limit determined by fluctuation is higher than the restriction imposed by the dispatch procedure from moment t−1 to moment t. The ramp-up rate limit needs to be set to r.sub.u.sup.VPG(t)=r.sub.u2.sup.VPG(t) to prevent power from rising too fast.

[0150] The ramp-down rate limit of the WVPG is calculated as follows:

r.sub.d.sup.VPG(t)=min[r.sub.d1.sup.VPG(t),r.sub.d2.sup.VPG(t)]  (14).

[0151] In formula (14), r.sub.d2.sup.VPG(t) denotes the ramp-down rate limit (MW/min) of the WVPG determined by the dispatch procedure at moment t, and is a negative value; r.sub.d.sup.VPG(t) denotes the ramp-down rate limit (MW/min) of the WVPG at moment t, and is a negative value.

[0152] If r.sub.d1.sup.VPG (t)<r.sub.d2.sup.VPG(t), indicates that output determined by fluctuation drops too fast from moment t−1 to moment t. The ramp-down rate limit needs to be set to r.sub.d.sup.VPG(t)=r.sub.d1.sup.VPG(t) to avoid underestimating the ramp-down rate.

[0153] If r.sub.d1.sup.VPG(t)>r.sub.d2.sup.VPG(t), itindicates that a downward fluctuation rate of output is less than a requirement of the dispatch procedure from moment t−1 to moment t. The ramp-down rate may be increased through a control method such as wind curtailment if necessary, and the ramp-down rate limit is set to r.sub.d.sup.VPG(t)=r.sub.d2.sup.VPG(t).

[0154] A ramp rate index is calculated as follows:

[0155] Conduct quantile operation on the multi-state probability distribution of fluctuation of the WVPG at moment t to obtain the ramp-up and ramp-down rate limits r.sub.u1.sup.VPG(t),r.sub.d1.sup.VPG(t) determined by output fluctuation. Then obtain values of the ramp rate limits of the WVPG at the current moment by using formula (13) and formula (14).

[0156] The analysis method in step S3 includes the following:

[0157] c1. variables

[0158] c1.1. integer variables [0159] 1) v(i,t): a variable that denotes an on/off state of unit i at moment t, where v(i,t)=1 0 t)=1 indicates that unit i is on at moment t, and v(i,t)=0 indicates that unit i is off at moment t; [0160] 2) y(i,t): a variable that denotes a conversion state of unit i at moment t, where y(i,t)=1 indicates that unit i changes from the off state to the on state at moment t, and y(i,t)=0 indicates that the unit is in another state; [0161] 3) z(i,t): a variable that denotes the conversion state of unit i at moment t, where z(i,t)=1 indicates that unit i changes from the on state to the off state at moment t, and z(i,t)=0 indicates that the unit is in another state; and [0162] 4)ω(i,tj): a variable that denotes a downtime state of the unit, where ω(i,t,j)=1 indicates that unit i has been off for j hours at moment t, and ω(i,t,j)=0 indicates that the unit is in another state.

[0163] Tr (i,t), a variable valued 0 or 1, further needs to be added to denote the startup time of a translatable load. Tr (i,t)=1 indicates that translatable load i starts to operate at moment t. Tr (i,t)=0 indicates that translatable load i does not start to operate at moment t.

[0164] c1.2. continuous variable

[0165] p(i,t): output of unit i at moment t.

[0166] c2. objective function

[00024]$\begin{array}{cc}\underset{p\left(i.t\right),y\left(i,t\right)}{\mathrm{Min}}{.Math.}_{t=1}^{TT}{.Math.}_{i=1}^{NG}\left[c\left(i,t\right)+b\left(i,t\right)\right],& \left(15\right)\end{array}$
where c(i,t)=A(i)×p(i,t).sup.2+B(i)×p(i,t)+C(i)×v(i,t)  (16).

[0167] In the formulas, TT denotes total operation duration; NG denotes the number of units; A(i), B(i), C(i) respectively denote a quadratic coefficient, a linear coefficient, and a constant term in a quadratic function of power generation costs and output power of unit i; c(i,t) denotes power generation costs of unit i at moment t; b(i,t) denotes startup costs of unit i at moment t; calculation of b(i,t) is related to downtime of the unit before startup, and detailed calculation is described in constraints; shutdown costs is usually not considered because it has little impact on total costs.

[0168] c3. constraints

[0169] c3.1. power balance constraint

[00025]$\begin{array}{cc}\underset{i}{\overset{NG}{.Math.}}p\left(i,t\right)=D_0\left(t\right)+\underset{i=1}{\overset{\mathrm{NI}}{.Math.}}PD\left(i,t\right)\forall t\in T.& \left(17\right)\end{array}$

[0170] In formula (17), PD(i,t) denotes a load value of translatable load i at moment t; D.sub.0(t) denotes a fixed system load at moment t; and NI denotes the number of translatable loads. The sum of the fixed loads and the translatable load after translation matches the total output of the unit based on this constraint.

[0171] c3.2. constraint on upper and lower output limits of the unit

p(i,t)×v(i,t)≤p(i,t)≤p(i,t)×v(i,t)∀i∈G,t∈T  (18).

[0172] p(i,t) denotes the lower output limit of unit i at moment t; p(i,t) denotes the upper output limit of unit i at moment t; G denotes a set of units. This constraint ensures that the units can operate in an output range for normal operation.

[0173] c3.3. constraints on startup/shutdown auxiliary variable

y(i,t)+z(i,t)≤1 ∀i∈G,t∈T  (19)

y(i,t)−z(i,t)=v(i,t)−v(i,t−1)∀i∈G,t∈T  (20).

[0174] Formula (19) imposes a restriction that prevents a unit from being started and shut down at the same time. Formula (20) restricts consistency between startup and shutdown actions and an operating state of the unit.

[0175] c3.4. constraints on minimum runtime

[00026]$\begin{array}{cc}\underset{k=1}{\overset{G\left(i\right)}{.Math.}}\left[1-v\left(i,k\right)\right]=0,& \left(21\right)\\ \underset{k=t}{\overset{t+T_{\mathrm{up}}\left(i\right)-1}{.Math.}}v\left(i,k\right)\ge T_{up}\left(i\right)\times y\left(i,t\right)\forall t=G\left(i\right)+1\mathrm{.Math.}\mathrm{TT}-T_{up}\left(i\right)+1& \left(22\right)\\ \underset{k=t}{\overset{TT}{.Math.}}\left[v\left(i,k\right)-y\left(i,t\right)\right]\ge 0\forall t=TT-T_{up}\left(i\right)+2\mathrm{.Math.}\mathrm{TT}.& \left(23\right)\end{array}$

[0176] In the formulas, G(i)=Min[TT,(T.sub.up(i)−T.sub.0(i))×U.sub.0(i)] denotes minimum runtime initially required by unit i after operation starts; T.sub.up(i) denotes minimum startup runtime of unit i; U.sub.0(i) denotes an initial operating state of unit I, and if U.sub.0(i)=1, unit i is initially in operation, and if U.sub.0(i)=0, unit i is initially in shutdown state; T.sup.0(i) denotes initial runtime or downtime of unit i, and is a positive number when it denotes the initial runtime and a negative number when it denotes the initial downtime. Impact of an operating state at the end of a previous day on start of operation on a new day is mainly considered in formulas (22) and (23). If the unit is initially in operation, formula (22) indicates that the unit requires at least a runtime constraint. Formula (23) denotes a minimum runtime constraint in last T.sub.up(i)−1 hours of runtime. If the unit starts in this time period, it needs to operate until the end of the day and cannot be shut down again.

[0177] c3.5. constraints on minimum downtime

[00027]$\begin{array}{cc}\underset{k=1}{\overset{L\left(i\right)}{.Math.}}v\left(i,k\right)=0,& \left(24\right)\\ \underset{k=t}{\overset{t+T_{\mathrm{down}}\left(i\right)-1}{.Math.}}\left[1-v\left(i,k\right)\right]\ge T_{\mathrm{down}}\left(i\right)\times z\left(i,t\right)\forall t=L\left(i\right)+1\mathrm{.Math.}\mathrm{TT}-T_{\mathrm{down}}\left(i\right)+1& \left(25\right)\\ \underset{k=t}{\overset{TT}{.Math.}}\left[1-v\left(i,k\right)-z\left(i,t\right)\right]\ge 0\forall t=TT-T_{\mathrm{down}}\left(i\right)+2\mathrm{.Math.}\mathrm{TT}.& \left(26\right)\end{array}$

[0178] In the formulas, L(i)=Min[TT,(T.sub.down(i)+T.sub.0(i))×(1−U.sub.0(i))], and it denotes minimum downtime initially required by unit i after operation starts; T.sub.down(i) denotes minimum downtime of unit i. Impact of an operating state at the end of a previous day on start of operation on a new day is mainly considered in formulas (25) and (26). If the unit is initially in shutdown state, it indicates that the unit requires at least a downtime constraint. Formula (25) denotes a minimum downtime constraint each time unit i is shut down. Formula (26) denotes a minimum downtime constraint in last T.sub.down(i)−1 hours of runtime. If the unit is shut down in the time period, it needs to be in shutdown state until the end of the day and cannot be started again.

[0179] c3.6. ramp rate constraints

p(i,t)−p(i,t−1)≤RATEUP(i)∀i∈G,t=2,3 . . . TT  (27)

p(i,t−1)−p(i,t)≤RATEDN(i)∀i∈G,t=2,3 . . . TT  (28).

[0180] RATEUP(i) denotes a ramp-up rate of unit i, that is, a rate at which output of unit i increases. RATEDN(i) denotes a ramp-down rate of unit i, that is, a rate at which output of unit i decreases. An output change between two adjacent time periods is limited because the output of the unit is adjusted at a limited rate. The ramp rate constraint represents this limitation.

[0181] c3.7. downtime counter constraints

S(i,t)≤S(i,t−1)+1 ∀t∈T

S(i,t)+(TT+1)×v(i,t)≥S(i,t−1)+1 ∀t∈T

S(i,t)−TT×[1−v(i,t)]≤0 ∀t∈T

S(i,t)≥0 ∀t∈T  (29).

[0182] S(i,t) is a downtime count variable that denotes accumulated consecutive downtime of unit i at moment t. It is necessary to use the constraint to record downtime because startup costs in the model are segmented and related to the downtime of the unit. A function implemented by the four constraints is as follows: when v(i,t)=0, S(i,t)=S(i, t−1)+1, and when v(i,t)=1 S(i,t)=0, thereby implementing accumulation of consecutive downtime.

[0183] c3.8. constraints on segmented startup costs

[00028]$\begin{array}{cc}b\left(i,t\right)=\underset{j=1}{\overset{ND}{.Math.}}\mathrm{UP}\mathrm{cost}\left(i,j\right)\times \omega \left(i,t,j\right)\forall t\in T\underset{j=1}{\overset{ND}{.Math.}}\omega \left(i,t,j\right)=y\left(i,t\right)\forall t\in TS\left(i,t-1\right)=\underset{j=1}{\overset{\mathrm{ND}-1}{.Math.}}j\times \omega \left(i,t,j\right)+m\left(i,t\right)\forall t\in Tm\left(i,t\right)\le TT\times \left[\omega \left(i,t,\mathrm{ND}\right)-y\left(i,t\right)+1\right]\forall t\in Tm\left(i,t\right)\ge \mathrm{ND}\times \omega \left(i,t,\mathrm{ND}\right)\forall t\in T.& \left(30\right)\end{array}$

[0184] In these constraints, segmented startup costs b(i,t) are calculated by using a counter of consecutive downtime. UP cos t(i,j) denotes startup costs after unit i is shut down for j hours. Longer downtime usually indicates higher startup costs. The startup costs no longer change after the unit is shut down for a certain time period. ND denotes a time threshold that finally allows the startup costs to stop changing in segmented startup costs calculation. ω(i,t,j) is a variable valued 0 or 1 that denotes a downtime state of the unit. ω(i,t,j)=1 indicates that unit i has been shut down for j hours at moment t, and ω(i,t,j)=0 indicates that the unit is in another state. In this constraint, for each time period of each unit, at most one ω(i,t,j) can be set to 1. In formula (30), dummy variable m(i,t) is used to associate consecutive downtime count S(i,t) with ω(i,t,j)

[0185] c3.9. constraint on period-based load translation

[00029]$\begin{array}{cc}PD\left(i,t\right)=\underset{j=1}{\overset{T_0}{.Math.}}\left[Tr\left(i,t-j+1\right)\times D_{\mathrm{tr}}\left(i,j\right)\right]\forall i\in NI,t\in NT\left(i\right)\underset{t=1}{\overset{TT}{.Math.}}Tr\left(i,t\right)=1,\forall i\in NI\underset{t=TT-NT\left(i\right)+1}{\overset{TT}{.Math.}}Tr\left(i,t\right)=0,\forall i\in \mathrm{NI}.& \left(31\right)\end{array}$

[0186] T.sub.0=min[t, NT(i)]. D.sub.tr(i,j) denotes a load value of a jth segment in a load sequence of translatable load i. NT(i) denotes a time length of the sequence of translatable load i. NI denotes a set of sequence numbers of translatable loads.

[0187] The distribution method in step S4 includes the following steps:

[0188] d1: Establish a general model of an optimal power distribution problem.

[0189] d2: Simplify the general model of the optimal power distribution problem.

[0190] d3: Solve the simplified model of the optimal power distribution problem.

[0191] d4: Conduct simulation analysis on a solved value.

[0192] In specific implementation, in the peak shaving control method for emergent source-grid coordination in case of a faulty sending-end power grid provided in the present disclosure, the dispatchability of the cluster virtual wind power unit is evaluated first, and the method for calculating the dispatchability index of the cluster wind power unit is developed. Then the source-load peak-shaving resource strategy is analyzed. Finally, the control strategy for tie-line peak shaving is distributed. In the present disclosure, real-time dispatchability of wind power participating in real-time power balance is first analyzed, specific evaluation indexes and calculation methods are provided, and calculation examples are given for verification. Then, an optimized real-time dispatch strategy is provided based on demand-side response resources, and DC and AC lines are coordinated for operation in terms of tie-lines. The peak shaving control method for emergent source-grid coordination when a sending-end power grid is faulty is provided to ensure normal operation.

[0193] 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, they should not be construed as a limitation to the present disclosure.

[0194] 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.

[0195] Although the embodiments of the present disclosure have been illustrated and described, it should be understood that a person 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.