CN105870962A - Robust interval wind power dispatching method of considering frequency response characteristic of power system - Google Patents

Abstract

The invention discloses a robust interval wind power dispatching method of considering the frequency response characteristic of a power system. The method comprises the following steps: calculating and obtaining frequency response characteristic data of the power system according to the frequency response characteristic data of a generator and the frequency response characteristic data of a load; calculating and obtaining target data of the minimum system operation cost; building constraint conditions of a robust interval wind power dispatching model according to the frequency response characteristic data of the power system and the target data of the minimum system operation cost; and carrying out robust interval wind power dispatching on a research object according to the constraint conditions. According to the method, the frequency response characteristic of the system in a wind power fluctuation state can be analyzed by building the quantitative relationship between a short-term wind power prediction error and the power fluctuation magnitude of the system; the shortage of considering the frequency response characteristic of the system in the wind-power active power dispatching process is compensated; the frequency fluctuation range of the system can meet the operating requirements under all possible wind power output change conditions; and the operation cost of the system is the minimum.

Description

A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic

Technical field

The present invention relates to Operation of Electric Systems and control technical field, espespecially a kind of consideration power system frequency response characteristic Robust Interval wind-powered electricity generation dispatching method.

Background technology

Wind power output has stronger randomness, intermittence and undulatory property, and the concentration of large-scale wind power accesses to power system Normal active balance and FREQUENCY CONTROL bring bigger challenge.Suddenly change or the disturbance feelings such as blower fan line tripping operation in wind power output Under condition, moment relatively large unbalance amount will be produced between system generation load, cause the big ups and downs of frequency.And by wind energy itself And the impact of asynchronous variable-speed generator characteristic, Wind turbines cannot provide inertia to ring for system as traditional synchronous generator Answer and primary frequency modulation service, be further exacerbated by the difficulty that system frequency controls.Meanwhile, along with the large-scale grid connection of wind-powered electricity generation and excellent First dispatching, the tradition grid-connected quantity of fired power generating unit, in progressively downward trend, causes system inertia level to be gradually lowered, makes in turn In the case of disturbance, system frequency fluctuation increases further.In the U.S., FERC and ERCOT has begun working on large-scale wind power and accesses Impact on system frequency response, and using system frequency response characteristics as regenerative resource grid connection capacity assessment factor it One.Existing result of study also indicates that, the grid-connected inertia that will deteriorate system of large-scale wind power and a frequency response, jeopardizes and is The safety of system, needs system to provide more effective wind-powered electricity generation control strategy and more inertia support.

The existing more research of system reserve optimization problem under accessing about large-scale wind power at present, but about extensive wind System frequency Study on Safety Problem under being electrically accessed is less.Due to spinning reserve it is generally required to the relatively long time cycle (5～ 10min) could activate to balance generation load supply and demand deviation, occur in the sudden change of unexpected disturbance such as wind power output or the tripping operation of blower fan line Time, system inertia can only be relied on or limit wind power output in advance to reduce instantaneous power amount of unbalance size to stop system frequency Fall.Therefore, needing a kind of wind-powered electricity generation dispatching method badly, in order to study, the most how to control wind power output real-time to meet system Running the requirement to frequency security, in other words under current system inertia and wind-powered electricity generation forecast error level, maximum receiving is how many Wind-powered electricity generation can still meet system frequency surge requirements.

Summary of the invention

In the case of the present invention is directed to large-scale wind power concentration access, owing to system inertia is not enough and wind power output fluctuates suddenly The system frequency safety problem caused, proposes a kind of it can be considered that the Robust Interval wind-powered electricity generation of electromotor and LOAD FREQUENCY response characteristic Rolling scheduling model, is used for carrying out wind-powered electricity generation scheduling.

Concrete, the Robust Interval wind-powered electricity generation dispatching method of the consideration power system frequency response characteristic that the present invention proposes, bag Include: step 1, according to generator frequency response characteristic data and LOAD FREQUENCY response characteristic data, calculate and obtain power system frequency Rate response characteristic data；Step 2, calculates the target data that system operation cost is minimum；Step 3, rings according to power system frequency Answer performance data and the target data of system operation cost minimum, set up the constraints of Robust Interval wind-powered electricity generation scheduling model；Step Rapid 4, according to constraints, object of study is carried out the wind-powered electricity generation scheduling of Robust Interval.

The Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic that the present invention proposes, it is possible to set up short Phase wind-powered electricity generation forecast error and the quantitative relationship of system frequency fluctuation size, to the system frequency response characteristics under wind-powered electricity generation fluctuation status It is analyzed, the deficiency during making up wind-powered electricity generation active power dispatch, system frequency response characteristics considered.And enable in institute Under possible wind power output situation of change, the frequency fluctuation scope of system is satisfied by service requirement, and system operation cost Minimum.

Accompanying drawing explanation

Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not Constitute limitation of the invention.In the accompanying drawings:

Fig. 1 is the Robust Interval wind-powered electricity generation dispatching method stream considering power system frequency response characteristic of one embodiment of the invention Cheng Tu.

Fig. 2 is generator speed and the transmission function schematic diagram of torque relationship of one embodiment of the invention.

Fig. 3 is speed regulator speed and the transmission function schematic diagram of torque relationship of one embodiment of the invention.

Fig. 4 is electromotor and the collective frequency response characteristic schematic diagram of load of one embodiment of the invention.

Fig. 5 is electromotor and the governor parameter value table of one embodiment of the invention.

Fig. 6 is the Robust Interval scheduling result schematic diagram considering system frequency response characteristics of one embodiment of the invention.

Fig. 7 be one embodiment of the invention disturbance under system frequency wave process contrast schematic diagram.

Fig. 8 is the maximum allowable interval contrast schematic diagram of exerting oneself of the wind-powered electricity generation under the different Γ values of one embodiment of the invention.

Fig. 9 is the system operation cost table under the different Γ values of one embodiment of the invention.

Detailed description of the invention

Hereinafter coordinate diagram and presently preferred embodiments of the present invention, the present invention is expanded on further for reaching predetermined goal of the invention institute The technological means taked.

Fig. 1 is the Robust Interval wind-powered electricity generation dispatching method stream considering power system frequency response characteristic of one embodiment of the invention Cheng Tu.Wherein, the method includes:

Step S101, according to generator frequency response characteristic data and LOAD FREQUENCY response characteristic data, calculates and obtains electricity Force system frequency response characteristic data；

Step S102, calculates and obtains the target data that system operation cost is minimum；

Step S103, according to the target data that power system frequency response characteristic data and system operation cost are minimum, builds The constraints of vertical Robust Interval wind-powered electricity generation scheduling model；

Step S104, according to constraints, carries out the wind-powered electricity generation scheduling of Robust Interval to object of study.

Concrete, in step 1, it is necessary first to obtain generator frequency response characteristic, LOAD FREQUENCY response characteristic.

The acquisition methods of generator frequency response characteristic is as follows:

In the case of system suffers unexpected disturbance, act on the various machine torques on generator rotor shaft and electromagnetic torque Imbalance occur, system frequency starts to fluctuate with given pace, and this can be obtained by the generator amature equation of motion:

$J\frac{d\mathrm{\ω}}{dt}=T_m-T_e;$

When wind power output or workload demand occur to change suddenly, by the electronic torque T of reflection to electromotor output_eChange Change, cause machine torque T_mWith electronic torque T_eDo not mate, cause the change of generator amature movement velocity in turn.Assuming that Initial mechanical torque and electromagnetic torque under systematic steady state ruuning situation are respectively T_m0、T_e0, then above formula can further indicate that For:

$J\frac{d\mathrm{\ω}}{dt}=T_{m0}+{\mathrm{\ΔT}}_m-T_{e0}-{\mathrm{\ΔT}}_e={\mathrm{\ΔT}}_m-{\mathrm{\ΔT}}_e={\mathrm{\ΔP}}_m-{\mathrm{\ΔP}}_e;$

Wherein, in above formula, physical quantity is perunit value form.This formula can pass through the transmission function representation shown in Fig. 2, Wherein, H is the machinery inertial time constant of rotor.

Further, it is considered to the control action of speed regulator, speed regulator measures rotor speed ω and in real time with synchronous rotational speed ω₀ Comparing, velocity deviation forms control signal after being integrated amplification, is used for regulating the main valve for steam channel of steam turbine Or the gate of the hydraulic turbine, and then realize the control action to generator mechanical watt level.This regulation process can pass through Fig. 3 Shown transmission function representation.Wherein, R is the speed variation of speed regulator；T_GFor servo time constant；T_RHFor the reheater time Constant；F_HPThe proportionality coefficient of steam turbine general power is accounted for for high pressure turbine stage power；T_CHBe main enter steam space and steam chest time normal Number.

The acquisition methods of the frequency response characteristic of load is as follows:

The frequency response characteristic of load is relevant to its kind: as the resistor-type load of illumination and heating load etc, its electricity Power is unrelated with frequency；But as motor and the induction-motor load of pump class, the fluctuation with frequency is changed by its rotating speed, cause output Electrical power change therewith.Therefore, the response characteristic that frequency is changed by load can approximate representation be both the above form sum:

ΔP_e=Δ P_L+DΔω；

Wherein, Δ P_LFor the loaded portion insensitive to frequency；D Δ ω is the loaded portion to frequency sensitive；D is load Damped coefficient.

The acquisition methods of power system frequency response characteristic data is as follows:

The electromotor obtained from above-mentioned calculating and the comprehensive frequency of load collective frequency response characteristic, electromotor and load Rate response characteristic can pass through the transmission function representation shown in Fig. 4.

As shown in Figure 4, the formula of power system frequency response characteristic data is as follows:

$\frac{\mathrm{\Δ}\mathrm{\ω}}{{\mathrm{\ΔP}}_L}=-\frac{1}{2Hs+D+\frac{1+F_{HP}{T}_{RH}s}{R(1+{\mathrm{sT}}_G)(1+{\mathrm{sT}}_{CH})(1+{\mathrm{sT}}_{RH})}};---\left(1\right)$

Wherein,For power system frequency response characteristic data, R is the speed variation of speed regulator, T_GFor servo time Constant, T_RHFor reheater time constant, F_HPThe proportionality coefficient of steam turbine general power, T is accounted for for high pressure turbine stage power_CHIt is that master enters Steam space and the time constant of steam chest, H is the machinery inertial time constant of rotor, and D is load damped coefficient, and s is differential operator, ΔP_LFor the loaded portion response characteristic data insensitive to frequency；

Wherein, speed regulator response characteristic is by T_RHImpact maximum, its usual value is 6～12s；And T_G、T_CHValue 0.2～ About 0.3s, affects relatively small, can ignore.

Meanwhile, further formula (1) can be expanded to multi-machine power system, multi-machine power system frequency response characteristic data Formula as follows:

$\frac{\mathrm{\Δ}\mathrm{\ω}}{{\mathrm{\ΔP}}_L}=-\frac{1}{2Hs+D+\underset{i=1}{\overset{N}{\Σ}}\frac{1+F_{HPi}{T}_{RHi}s}{R_i(1+{\mathrm{sT}}_{RHi})}};---\left(2\right)$

Wherein, the number of units of conventional power generation usage unit during N is system.

From formula (2), the change of disturbance starting stage frequency depends mainly on the size of disturbance and system inertia size.Greatly In the case of scale wind power integration, disturbance generally refers to fluctuation and the blower fan line tripping fault suddenly of wind power output；And system is used to Property not only includes the inertia time constant of generator amature, also includes the speed regulating constant of reheater time constant, speed regulator And the damping constant of load.

Further, shown by research, owing to frequency is to T_RHSensitivity is less, for convenience of calculation, sets all reheating Device time constant is T, meanwhile,

Order

Then the time domain response tables of data of Δ ω is shown as:

Δ ω (t) derivation is obtained maximum is:

Wherein,

In step 2, calculate and obtain the target data that system operation cost is minimum, including:

Operating cost includes that the coal consumption cost of conventional fired power generating unit and minimum are abandoned the system under wind requires and abandoned eolian, its In, the formula of operating cost is as follows:

$\begin{array}{c}f\left(p_{i,t}\right)=\min \left(\underset{t=t_0+1}{\overset{T}{\Σ}}\underset{i\∈G_{con}}{\Σ}\right(a_i{p}_{it}^2+b_i{p}_{it}+c_i)+\\ \underset{j\∈G_{wind}}{\Σ}\underset{t=t_0+1}{\overset{T}{\Σ}}{\mathrm{\λ}}_j({\stackrel{\‾}{p}}_{jt}^w-p_{jt}^{w,\max }));\end{array}---\left(5\right)$

Wherein: p_i,t、p_itIt is the plan of exerting oneself in the t period of i-th fired power generating unit, a_i、b_i、c_iFor predetermined coefficient, For in scheduling process allow wind-powered electricity generation EIAJ,EIAJ is predicted for wind-powered electricity generation.

According to equal incremental principle, when meeting λ_jRealize minimum during ＞ 0 and abandon the purpose of wind, calculated by formula (5) and obtain system The target data that system operating cost is minimum.

In step 3, constraints includes:

1, the system frequency security constraint under the most severe scene represents with following formula:

Wherein, Γ ∈ [0,1], for the uncertainty factor of total wind power output；Be respectively due to wind energy with The output of wind electric field variable quantity that machine fluctuation and the tripping operation of blower fan line cause；System frequency maximum for allowing in scheduling process becomes Change amount；G_windFor Wind turbines set；Bar number for all integrated wind plant collection electric line；M is maximum allowable collection Electric line tripping strip number；Total installation of generating capacity for jth wind energy turbine set；For the dress of kth bar collection electric line in jth wind energy turbine set Machine capacity；For the kth article collection electric line of jth wind energy turbine set in the running status of t period: take 1 time properly functioning, accident 0 is taken during tripping operation.

Wherein, formula (6-1) is the bound constraint of single output of wind electric field change, wherein exerting oneself jth wind energy turbine set Bound in being expressed as field all collection electric line of running exert oneself the summation of bound, be approximately considered collection electric line and hold by installation Amount distribution is exerted oneself；

Formula (6-2) is to exert oneself excursion constraint due to the blower fan line wind field that causes of tripping operation；

Formula (6-3) is the constraint of collection electric line all to system maximum allowable tripping strip number；

Formula (6-4) is the span constraint to collection electric line running status；

Formula (6-5) is the constraint of the instantaneous maximum sudden change amount of output of wind electric field all to the whole network；

Wherein, the value of uncertainty factor Γ reflects the balance between system performance driving economy and safety, and Γ is more Big then solution is the most conservative, security of system is the highest and economy is the poorest；Otherwise it is the most optimistic.

When Γ=0, showing not consider the probabilistic impact of wind power output in scheduling process, model degradation is tradition Economic load dispatching model.

Γ=1 shows to consider all possible wind power output situation in scheduling process, now the most conservative by obtaining Scheduling result.

2, the constraint of wind power output bound represents with following formula:

When considering to abandon wind factor, the wind power output plan bound not higher than wind-powered electricity generation prediction allowed in scheduling process is exerted oneself Interval bound,

$p_{jt}^{w,max}\≤{\stackrel{\‾}{p}}_{jt}^w;---\left(7\right)$

$p_{jt}^{w,min}\≤{\underset{\‾}{p}}_{jt}^w;---\left(8\right)$

Wherein,Minimum load is predicted for wind-powered electricity generation；Minimum load is allowed for wind-powered electricity generation.

3, the standby margin constraints of the rotation of the system under the most severe scene represents with following formula:

$\left\{\begin{array}{c}{u}_t=\underset{p_{jt}^{w,1}}{\min }(\underset{i\∈G_{con}}{\Σ}{p}_{i,t}+\underset{i\∈G_{con}}{\Σ}{R}_{it}^u+\underset{j\∈G_{wind}}{\Σ}{p}_{jt}^{w,3}-D_t)\≥0;---(9-1)\\ \begin{array}{cc}s.t.& \underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,2}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(9-2)\\ & \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m,\∀j\∈G_{wind};---(9-3)\\ & n_{jt}^k\∈\{0,1\}\end{array}\end{array}\right.$

$\left\{\begin{array}{c}{d}_t=\underset{p_{jt}^{w,2}}{\min }(D_t-\underset{i\∈G_{con}}{\Σ}{p}_{i,t}+\underset{i\∈G_{con}}{\Σ}{R}_{it}^d-\underset{j\∈G_{wind}}{\Σ}{p}_{jt}^{w,4})\≥0;---(10-1)\\ \begin{array}{cc}s.t.& \underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,3}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(10-2)\\ & \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m,\∀j\∈G_{wind};---(10-3)\\ & n_{jt}^k\∈\{0,1\}\end{array}\end{array}\right.$

Wherein, u_t、d_tIt is respectively the minimum upper and lower rotation of system of t period for nargin；p_i,tAndIt is respectively i-th Platform fired power generating unit the t period exert oneself plan and upper and lower rotation for capacity；It is respectively the wind-powered electricity generation under two kinds of scenes Go out force value；D_tIt it is the system load demand of t period；G_conFor conventional rack set.

4, the transmission section security constraint under the most severe scene represents with following formula:

$\left\{\begin{array}{c}\underset{p_{jt}^{w,3}}{\max }(\underset{i\∈G_{con}}{\Σ}{k}_{li}{p}_{it}+\underset{j\∈G_{wind}}{\Σ}{k}_{lj}{p}_{jt}^{w,5})\≤\stackrel{\‾}{T_l};---(11-1)\\ s.t.\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,4}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(11-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m;---(11-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

$\left\{\begin{array}{c}\underset{p_{jt}^{w,4}}{\min }(\underset{i\∈G_{con}}{\Σ}{k}_{li}{p}_{it}+\underset{j\∈G_{wind}}{\Σ}{k}_{lj}{p}_{jt}^{w,6})\≥\underset{\‾}{T_l};---(12-1)\\ s.t.\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,5}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(12-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m;---(12-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

Wherein, l is section numbering, l=1,2 ..., L, L are total section number；k_liIt is i-th unit spirit to l section Sensitivity；It is respectively the wind power output under two kinds of scenes.

5, the rotation of conventional power unit represents with following formula for capacity-constrained:

$R_{it}^u\≤\min ({\stackrel{\‾}{p}}_i-p_{it},p_{i,t-1}+{\mathrm{\Δp}}_{u,i}{T}^{\′}-p_{it},{\mathrm{\Δp}}_{u,i}{T}^{\′});---\left(13\right)$

$R_{it}^d\≤\min (p_{it}-{\underset{\‾}{p}}_i,p_{it}-p_{i,t-1}+{\mathrm{\Δp}}_{d,i}{T}^{\′},{\mathrm{\Δp}}_{d,i}{T}^{\′});---\left(14\right)$

Wherein, T ' is the sampling interval.

6, the limit value constraint of exerting oneself of conventional power unit represents with following formula:

${\underset{\‾}{p}}_i\≤p_{i,t}\≤{\stackrel{\‾}{p}}_i;---\left(15\right)$

Wherein, p _iIt is respectively the bound of exerting oneself of conventional power unit.

7, the climbing rate constraint of conventional power unit represents with following formula:

p_i,t-1-Δp_d,iT≤p_it≤p_i,t-1+Δp_u,iT；(16)

Wherein, Δ p_ui、Δp_diIt is respectively conventional power unit climbing rate upwardly or downwardly.

Wind turbines is controlled to use Robust Interval control model by the present invention.On the one hand the selection in wind power output interval should expire Minimum in terms of foot economy abandons wind requirement, on the other hand should meet the system safety in operation in the case of the most severe wind power output Requirement.Set up following consideration electromotor and the Robust Interval wind-powered electricity generation scheduling model of LOAD FREQUENCY response characteristic accordingly.

The solution of Robust Interval wind-powered electricity generation scheduling model (formula (5)-formula (16)) is to try to achieve such a optimum thermoelectricity to exert oneself meter Draw value p_itAnd wind-powered electricity generation high output is intervalIt makes under this kind of plan arrangement, and system always has enough Nargin the system frequency fluctuation, standby deficiency and the section that cause due to wind-powered electricity generation forecast error or blower fan line fault trip with reply The security of system problems such as through-put power is out-of-limit, and system economy is optimum under this kind of plan arrangement.

It can be seen that Robust Interval wind-powered electricity generation scheduling model (formula (5)-formula (16)) is two-layer hybrid integer optimization problem, on There is coupling by wind-powered electricity generation high output interval variable in lower floor's Optimized model, it is impossible to direct solution.But under the feature of this model is Layer optimization problem is linear programming problem, and the object function of lower floor's optimization problem to participate in upper strata with the form of constraints excellent During change.Therefore, according to character and the strong dual principle of linear programming, the integer variable of lower floor's optimization problem can be expanded to Continuous variable on [0,1] interval, simultaneously can be by lower floor's optimization problem by its dual problem equivalencing, thus can be by former Problem is converted into the quadratic programming problem of monolayer, and can conveniently use prim al-dual interior point m ethod to solve.

Below in conjunction with step 4, according to constraints, a concrete object of study is carried out the wind-powered electricity generation scheduling of Robust Interval.

With IEEE RTS system as object of study, the #14 unit at #13 bus is replaced with an installed capacity is The #1 wind energy turbine set of 600MW, meanwhile, at #7 bus, one installed capacity of interpolation is the #2 wind energy turbine set of 350MW.Conventional generator is climbed Ratio of slope is taken as the 1% of rated capacity, and the sampling interval is 5min.Load damped coefficient D=1, T=8s, electromotor and speed regulator ginseng Number value refers to shown in Fig. 5.

1, to Γ value 0.4,Value 0.2Hz, now, it is considered to maximum with wind-powered electricity generation when not considering that system frequency retrains Allow to exert oneself interval as shown in Figure 6.

From fig. 6, it can be seen that when considering system frequency response characteristics, the maximum allowable interval of exerting oneself of wind-powered electricity generation is when major part Section is essentially identical with when not considering system frequency response characteristics, change greatly in part wind power output interval the period (period 3, time Section 12 decreases, and this is mainly for preventing unexpected wind power output fluctuation from causing system frequency to exceed normal consistency permission model Enclose.As a example by the period 12, under two kinds of wind power output projected conditions of Fig. 6, when wind power output occurs to fluctuate suddenly, it is considered to Do not consider system peak frequency wave process when system frequency retrains as it is shown in fig. 7, wherein emulation Period Length is taken as 20s, Simulation step length 0.1s.

As seen from Figure 7, when considering system frequency response, in the case of disturbance, system frequency fluctuation maximum is-0.2Hz, In disturbance, about about 10s occurring, system frequency reaches stationary value；And when not considering system frequency response characteristics, system under disturbance Frequency fluctuation maximum is up to-0.8Hz, beyond the mains frequency security requirement under normal operation.

2, rightValue 0.2Hz, uncertainty factor Γ value 0.2,0.5,0.9 respectively, now, wind-powered electricity generation is maximum allowable Exert oneself interval as shown in Figure 8, system operation cost is as shown in Figure 9.

As seen from Figure 8, along with the increase of Γ value, the maximum allowable interval upper limit of exerting oneself of wind-powered electricity generation is gradually lowered, and meanwhile, system is transported Row cost increases, and this is basically identical with the balance conclusion of safety to system performance driving economy with part 2.

The Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic that the present invention proposes, it is possible to set up short Phase wind-powered electricity generation forecast error and the quantitative relationship of system frequency fluctuation size, to the system frequency response characteristics under wind-powered electricity generation fluctuation status It is analyzed, the deficiency during making up wind-powered electricity generation active power dispatch, system frequency response characteristics considered.And enable in institute Under possible wind power output situation of change, the frequency fluctuation scope of system is satisfied by service requirement, and system operation cost Minimum.

Particular embodiments described above, has been carried out the purpose of the present invention, technical scheme and beneficial effect the most in detail Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the present invention, the guarantor being not intended to limit the present invention Protect scope, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, should be included in this Within the protection domain of invention.

Claims (12)

1. the Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic, it is characterised in that the method bag Include:

Step 1, according to generator frequency response characteristic data and LOAD FREQUENCY response characteristic data, calculates and obtains power system frequency Rate response characteristic data；

Step 2, calculates and obtains the target data that system operation cost is minimum；

Step 3, according to the target data that power system frequency response characteristic data and system operation cost are minimum, sets up robust district Between the constraints of wind-powered electricity generation scheduling model；

Step 4, according to constraints, carries out the wind-powered electricity generation scheduling of Robust Interval to object of study.

Method the most according to claim 1, it is characterised in that in step 1, according to generator frequency response characteristic data And LOAD FREQUENCY response characteristic data, calculate and obtain power system frequency response characteristic data, including:

The formula of power system frequency response characteristic data is as follows:

$\frac{\mathrm{\Δ}\mathrm{\ω}}{{\mathrm{\ΔP}}_L}=-\frac{1}{2Hs+D+\frac{1+F_{HP}{T}_{RH}s}{R(1+{\mathrm{sT}}_G)(1+{\mathrm{sT}}_{CH})(1+{\mathrm{sT}}_{RH})}};---\left(1\right)$

Wherein,For power system frequency response characteristic data, R is the speed variation of speed regulator, T_GNormal for servo time Number, T_RHFor reheater time constant, F_HPThe proportionality coefficient of steam turbine general power, T is accounted for for high pressure turbine stage power_CHIt is that master enters vapour Volume and the time constant of steam chest, H is the machinery inertial time constant of rotor, and D is load damped coefficient, and s is differential operator, Δ P_LFor the loaded portion response characteristic data insensitive to frequency；

Formula (1) expands to multi-machine power system, and the formula of multi-machine power system frequency response characteristic data is as follows:

$\frac{\mathrm{\Δ}\mathrm{\ω}}{{\mathrm{\ΔP}}_L}=-\frac{1}{2Hs+D+\underset{i=1}{\overset{N}{\Σ}}\frac{1+F_{HPi}{T}_{RHi}s}{R_i(1+{\mathrm{sT}}_{RHi})}};---\left(2\right)$

Wherein, the number of units of conventional power generation usage unit during N is system.

Method the most according to claim 2, it is characterised in that owing to frequency is to T_RHSensitivity is less, for convenience of calculation, if Fixed all reheater time constants are T, meanwhile,

Order

Then the time domain response tables of data of Δ ω is shown as:

Δ ω (t) derivation is obtained maximum is:

Wherein,

Method the most according to claim 3, it is characterised in that in step 2, calculates and obtains system operation cost minimum Target data, including:

Operating cost includes that the coal consumption cost of conventional fired power generating unit and minimum are abandoned the system under wind requires and abandoned eolian, wherein, and fortune The formula of row cost is as follows:

$\begin{array}{c}f\left(p_{i,t}\right)=\min \left(\underset{t=t_0+1}{\overset{T}{\Σ}}\underset{i\∈G_{con}}{\Σ}\right(a_i{p}_{it}^2+b_i{p}_{it}+c_i)+\\ \underset{j\∈G_{wind}}{\Σ}\underset{t=t_0+1}{\overset{T}{\Σ}}{\mathrm{\λ}}_j\left({\stackrel{\‾}{p}}_{jt}^w-p_{jw}^{w,\max }\right));\end{array}---\left(5\right)$

Wherein: p_i,t、p_itIt is the plan of exerting oneself in the t period of i-th fired power generating unit, a_i、b_i、c_iFor predetermined coefficient,For adjusting The wind-powered electricity generation EIAJ allowed during degree,EIAJ is predicted for wind-powered electricity generation；

According to equal incremental principle, when meeting λ_j> 0 time realize minimum and abandon the purpose of wind, calculate acquisition system by formula (5) and run The target data of cost minimization.

Method the most according to claim 4, it is characterised in that in step 3, according to power system frequency response characteristic number According to and the minimum target data of system operation cost, set up the constraints of Robust Interval wind-powered electricity generation scheduling model, wherein, retrain bar Part includes: the system rotation under the constraint of system frequency security constraint under the most severe scene, wind power output bound, the most severe scene Transmission section security constraint under standby margin constraints, the most severe scene, the standby capacity-constrained of rotation of conventional power unit, the going out of conventional power unit The constraint of power limit value, the climbing rate constraint of conventional power unit.

Method the most according to claim 5, it is characterised in that the system frequency security constraint under the most severe described scene is used Following formula represents:

Wherein, Γ ∈ [0,1], for the uncertainty factor of total wind power output；It is respectively due to wind energy random fluctuation And the output of wind electric field variable quantity that the tripping operation of blower fan line causes；For the system frequency maximum variable quantity allowed in scheduling process； G_windFor Wind turbines set；Bar number for all integrated wind plant collection electric line；M is maximum allowable collection electric line Tripping strip number；Total installation of generating capacity for jth wind energy turbine set；Hold for the installation of kth bar collection electric line in jth wind energy turbine set Amount；For the kth article collection electric line of jth wind energy turbine set in the running status of t period: take 1 time properly functioning, emergency stop valve trip Time take 0；

Wherein, formula (6-1) is the bound constraint of single output of wind electric field change, wherein exerting oneself jth wind energy turbine set up and down Limit in being expressed as field all collection electric line of running exert oneself the summation of bound, be approximately considered collection electric line and divide by installed capacity Allot power；

Formula (6-2) is to exert oneself excursion constraint due to the blower fan line wind field that causes of tripping operation；

Formula (6-3) is the constraint of collection electric line all to system maximum allowable tripping strip number；

Formula (6-4) is the span constraint to collection electric line running status；

Formula (6-5) is the constraint of the instantaneous maximum sudden change amount of output of wind electric field all to the whole network；

Wherein, the value of uncertainty factor Γ reflects the balance between system performance driving economy and safety, and Γ is the biggest then Xie Yue guards, and security of system is the highest and economy is the poorest；When Γ=0, show not consider wind power output in scheduling process Probabilistic impact, model degradation is traditional economic load dispatching model；Γ=1 show to consider in scheduling process all can The wind power output situation of energy, now will obtain the most conservative scheduling result.

Method the most according to claim 5, it is characterised in that described wind power output bound retrains with following formula table Show:

When considering to abandon wind factor, the wind power output plan bound not higher than wind-powered electricity generation allowed in scheduling process predicts interval of exerting oneself Bound,

$p_{jt}^{w,max}\≤{\stackrel{\‾}{p}}_{jt}^w;---\left(7\right)$

$p_{jt}^{w,min}\≤{\underset{\‾}{p}}_{jt}^w;---\left(8\right)$

Wherein,Minimum load is predicted for wind-powered electricity generation；Minimum load is allowed for wind-powered electricity generation.

Method the most according to claim 5, it is characterised in that the standby margin constraints of system rotation under the most severe described scene is used Following formula represents:

$\left\{\begin{array}{c}{u}_t=\underset{p_{jt}^{w,1}}{min}(\underset{i\∈G_{con}}{\Σ}{p}_{i,t}+\underset{i\∈G_{con}}{\Σ}{R}_{it}^u+\underset{j\∈G_{wind}}{\Σ}{p}_{jt}^{w,3}-D_t)\≥0;---(9-1)\\ \begin{array}{cc}s.t.& \underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,2}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\end{array};---(9-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m,\∀j\∈G_{wind};---(9-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

$\left\{\begin{array}{c}{d}_t=\underset{p_{jt}^{w,2}}{min}(D_t-\underset{i\∈G_{con}}{\Σ}{p}_{i,t}+\underset{i\∈G_{con}}{\Σ}{R}_{it}^d-\underset{j\∈G_{wind}}{\Σ}{p}_{jt}^{w,4})\≥0;---(10-1)\\ \begin{array}{cc}s.t.& \underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,3}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max }\end{array};---(10-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m,\∀j\∈G_{wind};---(10-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

Wherein, u_t、d_tIt is respectively the minimum upper and lower rotation of system of t period for nargin；p_i,tAndIt is respectively i-th thermoelectricity Unit the t period exert oneself plan and upper and lower rotation for capacity；It is respectively the wind power output value under two kinds of scenes； D_tIt it is the system load demand of t period；G_conFor conventional rack set.

Method the most according to claim 5, it is characterised in that the transmission section security constraint under the most severe described scene is used Following formula represents:

$\left\{\begin{array}{c}\underset{p_{jt}^{w,3}}{\max }(\underset{i\∈G_{con}}{\Σ}{k}_{li}{p}_{it}+\underset{j\∈G_{wind}}{\Σ}{k}_{lj}{p}_{jt}^{w,5})\≥\stackrel{\‾}{T_l};---(11-1)\\ s.t.\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,4}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(11-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m;---(11-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

$\left\{\begin{array}{c}\underset{p_{jt}^{w,4}}{min}(\underset{i\∈G_{con}}{\Σ}{k}_{li}{p}_{it}+\underset{j\∈G_{wind}}{\Σ}{k}_{lj}{p}_{jt}^{w,6})\≥\underset{\‾}{T_l};---(12-1)\\ s.t.\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\min }\≤p_{jt}^{w,5}\≤\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\frac{{\stackrel{\‾}{p}}_j^k}{{\stackrel{\‾}{p}}_j}{p}_{jt}^{w,\max };---(12-2)\\ \underset{j\∈G_{wind}}{\Σ}\underset{k=1}{\overset{L}{\Σ}}{n}_{jt}^k\≥N_{G_{wind}}-m;---(12-3)\\ n_{jt}^k\∈\{0,1\}\end{array}\right.$

Wherein, l is section numbering, l=1,2 ..., L, L are total section number；k_liIt is sensitive to l section of i-th unit Degree；It is respectively the wind power output under two kinds of scenes.

Method the most according to claim 5, it is characterised in that the rotation of described conventional power unit is for the following calculation of capacity-constrained Formula represents:

$R_{it}^u\≤min({\stackrel{\‾}{p}}_i-p_{it},p_{i,t-1}+{\mathrm{\Δp}}_{u,i}{T}^{\′}-p_{it},{\mathrm{\Δp}}_{u,i}{T}^{\′});---\left(13\right)$

$R_{it}^d\≤min(p_{it}-{\underset{\‾}{p}}_i,p_{it}-p_{i,t-1}+{\mathrm{\Δp}}_{d,i}{T}^{\′},{\mathrm{\Δp}}_{d,i}{T}^{\′});---\left(14\right)$

Wherein, T ' is the sampling interval.

11. methods according to claim 5, it is characterised in that following calculation is used in the limit value constraint of exerting oneself of described conventional power unit Formula represents:

${\underset{\‾}{p}}_i\≤p_{i,t}\≤{\stackrel{\‾}{p}}_i;---\left(15\right)$

Wherein, p _iIt is respectively the bound of exerting oneself of conventional power unit.

12. methods according to claim 5, it is characterised in that following formula is used in the climbing rate constraint of described conventional power unit Represent:

p_i,t-1-Δp_d,iT≤p_it≤p_i,t-1+Δp_u,iT； (16)

Wherein, Δ p_ui、Δp_diIt is respectively conventional power unit climbing rate upwardly or downwardly.