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

Robust interval wind power dispatching method of considering frequency response characteristic of power system Download PDF

Info

Publication number
CN105870962A
CN105870962A CN201610195678.9A CN201610195678A CN105870962A CN 105870962 A CN105870962 A CN 105870962A CN 201610195678 A CN201610195678 A CN 201610195678A CN 105870962 A CN105870962 A CN 105870962A
Authority
CN
China
Prior art keywords
wind
sigma
overbar
response characteristic
frequency response
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610195678.9A
Other languages
Chinese (zh)
Other versions
CN105870962B (en
Inventor
陈建华
郭子明
张昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Jibei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jibei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Jibei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610195678.9A priority Critical patent/CN105870962B/en
Publication of CN105870962A publication Critical patent/CN105870962A/en
Application granted granted Critical
Publication of CN105870962B publication Critical patent/CN105870962B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Wind Motors (AREA)
  • Control Of Eletrric Generators (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

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 d ω d t = T m - T e ;
When wind power output or workload demand occur to change suddenly, by the electronic torque T of reflection to electromotor outputeChange Change, cause machine torque TmWith electronic torque TeDo 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 Tm0、Te0, then above formula can further indicate that For:
J d ω d t = T m 0 + ΔT m - T e 0 - ΔT e = ΔT m - ΔT e = ΔP m - Δ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 ω0 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;TGFor servo time constant;TRHFor the reheater time Constant;FHPThe proportionality coefficient of steam turbine general power is accounted for for high pressure turbine stage power;TCHBe 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:
ΔPe=Δ PL+DΔω;
Wherein, Δ PLFor 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:
Δ ω ΔP L = - 1 2 H s + D + 1 + F H P T R H s R ( 1 + sT G ) ( 1 + sT C H ) ( 1 + sT R H ) ; - - - ( 1 )
Wherein,For power system frequency response characteristic data, R is the speed variation of speed regulator, TGFor servo time Constant, TRHFor reheater time constant, FHPThe proportionality coefficient of steam turbine general power, T is accounted for for high pressure turbine stage powerCHIt 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, ΔPLFor the loaded portion response characteristic data insensitive to frequency;
Wherein, speed regulator response characteristic is by TRHImpact maximum, its usual value is 6~12s;And TG、TCHValue 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:
Δ ω ΔP L = - 1 2 H s + D + Σ i = 1 N 1 + F H P i T R H i s R i ( 1 + sT R H i ) ; - - - ( 2 )
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 TRHSensitivity 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:
f ( p i , t ) = min ( Σ t = t 0 + 1 T Σ i ∈ G c o n ( a i p i t 2 + b i p i t + c i ) + Σ j ∈ G w i n d Σ t = t 0 + 1 T λ j ( p ‾ j t w - p j t w , max ) ) ; - - - ( 5 )
Wherein: pi,t、pitIt is the plan of exerting oneself in the t period of i-th fired power generating unit, ai、bi、ciFor 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;GwindFor 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 j t w , m a x ≤ p ‾ j t w ; - - - ( 7 )
p j t w , m i n ≤ p ‾ j t w ; - - - ( 8 )
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:
u t = min p j t w , 1 ( Σ i ∈ G c o n p i , t + Σ i ∈ G c o n R i t u + Σ j ∈ G w i n d p j t w , 3 - D t ) ≥ 0 ; - - - ( 9 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 2 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 9 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m , ∀ j ∈ G w i n d ; - - - ( 9 - 3 ) n j t k ∈ { 0 , 1 }
d t = min p j t w , 2 ( D t - Σ i ∈ G c o n p i , t + Σ i ∈ G c o n R i t d - Σ j ∈ G w i n d p j t w , 4 ) ≥ 0 ; - - - ( 10 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 3 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 10 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m , ∀ j ∈ G w i n d ; - - - ( 10 - 3 ) n j t k ∈ { 0 , 1 }
Wherein, ut、dtIt is respectively the minimum upper and lower rotation of system of t period for nargin;pi,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;DtIt it is the system load demand of t period;GconFor conventional rack set.
4, the transmission section security constraint under the most severe scene represents with following formula:
max p j t w , 3 ( Σ i ∈ G c o n k l i p i t + Σ j ∈ G w i n d k l j p j t w , 5 ) ≤ T l ‾ ; - - - ( 11 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 4 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 11 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m ; - - - ( 11 - 3 ) n j t k ∈ { 0 , 1 }
min p j t w , 4 ( Σ i ∈ G c o n k l i p i t + Σ j ∈ G w i n d k l j p j t w , 6 ) ≥ T l ‾ ; - - - ( 12 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 5 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 12 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m ; - - - ( 12 - 3 ) n j t k ∈ { 0 , 1 }
Wherein, l is section numbering, l=1,2 ..., L, L are total section number;kliIt 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 i t u ≤ min ( p ‾ i - p i t , p i , t - 1 + Δp u , i T ′ - p i t , Δp u , i T ′ ) ; - - - ( 13 )
R i t d ≤ min ( p i t - p ‾ i , p i t - p i , t - 1 + Δp d , i T ′ , Δp d , i T ′ ) ; - - - ( 14 )
Wherein, T ' is the sampling interval.
6, the limit value constraint of exerting oneself of conventional power unit represents with following formula:
p ‾ i ≤ p i , t ≤ p ‾ i ; - - - ( 15 )
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:
pi,t-1-Δpd,iT≤pit≤pi,t-1+Δpu,iT;(16)
Wherein, Δ pui、ΔpdiIt 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 pitAnd 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:
Δ ω ΔP L = - 1 2 H s + D + 1 + F H P T R H s R ( 1 + sT G ) ( 1 + sT C H ) ( 1 + sT R H ) ; - - - ( 1 )
Wherein,For power system frequency response characteristic data, R is the speed variation of speed regulator, TGNormal for servo time Number, TRHFor reheater time constant, FHPThe proportionality coefficient of steam turbine general power, T is accounted for for high pressure turbine stage powerCHIt 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, Δ PLFor 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:
Δ ω ΔP L = - 1 2 H s + D + Σ i = 1 N 1 + F H P i T R H i s R i ( 1 + sT R H i ) ; - - - ( 2 )
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 TRHSensitivity 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:
f ( p i , t ) = min ( Σ t = t 0 + 1 T Σ i ∈ G c o n ( a i p i t 2 + b i p i t + c i ) + Σ j ∈ G w i n d Σ t = t 0 + 1 T λ j ( p ‾ j t w - p j w w , max ) ) ; - - - ( 5 )
Wherein: pi,t、pitIt is the plan of exerting oneself in the t period of i-th fired power generating unit, ai、bi、ciFor 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; GwindFor 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 j t w , m a x ≤ p ‾ j t w ; - - - ( 7 )
p j t w , m i n ≤ p ‾ j t w ; - - - ( 8 )
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:
u t = m i n p j t w , 1 ( Σ i ∈ G c o n p i , t + Σ i ∈ G c o n R i t u + Σ j ∈ G w i n d p j t w , 3 - D t ) ≥ 0 ; - - - ( 9 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 2 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ; - - - ( 9 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m , ∀ j ∈ G w i n d ; - - - ( 9 - 3 ) n j t k ∈ { 0 , 1 }
d t = m i n p j t w , 2 ( D t - Σ i ∈ G c o n p i , t + Σ i ∈ G c o n R i t d - Σ j ∈ G w i n d p j t w , 4 ) ≥ 0 ; - - - ( 10 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 3 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 10 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m , ∀ j ∈ G w i n d ; - - - ( 10 - 3 ) n j t k ∈ { 0 , 1 }
Wherein, ut、dtIt is respectively the minimum upper and lower rotation of system of t period for nargin;pi,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; DtIt it is the system load demand of t period;GconFor 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:
max p j t w , 3 ( Σ i ∈ G c o n k l i p i t + Σ j ∈ G w i n d k l j p j t w , 5 ) ≥ T l ‾ ; - - - ( 11 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 4 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 11 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m ; - - - ( 11 - 3 ) n j t k ∈ { 0 , 1 }
m i n p j t w , 4 ( Σ i ∈ G c o n k l i p i t + Σ j ∈ G w i n d k l j p j t w , 6 ) ≥ T l ‾ ; - - - ( 12 - 1 ) s . t . Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , min ≤ p j t w , 5 ≤ Σ k = 1 L n j t k p ‾ j k p ‾ j p j t w , max ; - - - ( 12 - 2 ) Σ j ∈ G w i n d Σ k = 1 L n j t k ≥ N G w i n d - m ; - - - ( 12 - 3 ) n j t k ∈ { 0 , 1 }
Wherein, l is section numbering, l=1,2 ..., L, L are total section number;kliIt 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 i t u ≤ m i n ( p ‾ i - p i t , p i , t - 1 + Δp u , i T ′ - p i t , Δp u , i T ′ ) ; - - - ( 13 )
R i t d ≤ m i n ( p i t - p ‾ i , p i t - p i , t - 1 + Δp d , i T ′ , Δp d , i T ′ ) ; - - - ( 14 )
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:
p ‾ i ≤ p i , t ≤ p ‾ i ; - - - ( 15 )
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:
pi,t-1-Δpd,iT≤pit≤pi,t-1+Δpu,iT; (16)
Wherein, Δ pui、ΔpdiIt is respectively conventional power unit climbing rate upwardly or downwardly.
CN201610195678.9A 2016-03-31 2016-03-31 A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic Expired - Fee Related CN105870962B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610195678.9A CN105870962B (en) 2016-03-31 2016-03-31 A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610195678.9A CN105870962B (en) 2016-03-31 2016-03-31 A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic

Publications (2)

Publication Number Publication Date
CN105870962A true CN105870962A (en) 2016-08-17
CN105870962B CN105870962B (en) 2019-03-12

Family

ID=56626699

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610195678.9A Expired - Fee Related CN105870962B (en) 2016-03-31 2016-03-31 A kind of Robust Interval wind-powered electricity generation dispatching method considering power system frequency response characteristic

Country Status (1)

Country Link
CN (1) CN105870962B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203855A (en) * 2017-08-03 2017-09-26 国网江苏省电力公司宿迁供电公司 The robust bi-level optimization model and conversion equivalent method of the Real-Time Scheduling containing wind power system
CN109980660A (en) * 2019-04-17 2019-07-05 陈建华 The Robust Interval wind-powered electricity generation dispatching method of frequency response characteristic is improved based on generator
CN110492531A (en) * 2019-08-06 2019-11-22 清华大学 Consider the electric power system dispatching operation method and system of synchronous rotary inertia level
CN110829487A (en) * 2019-10-22 2020-02-21 电子科技大学 Dynamic frequency prediction method for power system
CN111740416A (en) * 2020-07-08 2020-10-02 国网河南省电力公司经济技术研究院 Target cascade analysis-based receiving-end power grid direct-current maximum feed-in quantity evaluation method
WO2021109579A1 (en) * 2019-12-04 2021-06-10 国网青海省电力公司电力科学研究院 Unit combination method taking wind turbine generator security domain with frequency response control into consideration
CN113078639A (en) * 2021-04-14 2021-07-06 上海电力大学 High-proportion wind power system optimal scheduling method considering air conditioner load and frequency response

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104753086A (en) * 2015-04-09 2015-07-01 广西大学 Minimum spinning reserve scheduling method meeting requirements of robust running of various scenes
CN104868480A (en) * 2015-05-06 2015-08-26 清华大学 Active power and reactive power-coordinated active power distribution network robust voltage control method
CN104935017A (en) * 2015-05-19 2015-09-23 同济大学 Wind turbine generating unit and thermal power generating unit combination method based on light robust optimization model

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104753086A (en) * 2015-04-09 2015-07-01 广西大学 Minimum spinning reserve scheduling method meeting requirements of robust running of various scenes
CN104868480A (en) * 2015-05-06 2015-08-26 清华大学 Active power and reactive power-coordinated active power distribution network robust voltage control method
CN104935017A (en) * 2015-05-19 2015-09-23 同济大学 Wind turbine generating unit and thermal power generating unit combination method based on light robust optimization model

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
付娅: ""含大规模风电的互连系统的负荷频率控制"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107203855A (en) * 2017-08-03 2017-09-26 国网江苏省电力公司宿迁供电公司 The robust bi-level optimization model and conversion equivalent method of the Real-Time Scheduling containing wind power system
CN109980660A (en) * 2019-04-17 2019-07-05 陈建华 The Robust Interval wind-powered electricity generation dispatching method of frequency response characteristic is improved based on generator
CN110492531A (en) * 2019-08-06 2019-11-22 清华大学 Consider the electric power system dispatching operation method and system of synchronous rotary inertia level
CN110829487A (en) * 2019-10-22 2020-02-21 电子科技大学 Dynamic frequency prediction method for power system
WO2021109579A1 (en) * 2019-12-04 2021-06-10 国网青海省电力公司电力科学研究院 Unit combination method taking wind turbine generator security domain with frequency response control into consideration
US11777313B2 (en) 2019-12-04 2023-10-03 State Grid Qinghai Electric Power Research Institute Unit commitment method considering security region of wind turbine generator with frequency response control
CN111740416A (en) * 2020-07-08 2020-10-02 国网河南省电力公司经济技术研究院 Target cascade analysis-based receiving-end power grid direct-current maximum feed-in quantity evaluation method
CN113078639A (en) * 2021-04-14 2021-07-06 上海电力大学 High-proportion wind power system optimal scheduling method considering air conditioner load and frequency response

Also Published As

Publication number Publication date
CN105870962B (en) 2019-03-12

Similar Documents

Publication Publication Date Title
CN105870962A (en) Robust interval wind power dispatching method of considering frequency response characteristic of power system
Rezkalla et al. Electric power system inertia: requirements, challenges and solutions
Fischer et al. Operational experiences with inertial response provided by type 4 wind turbines
Ullah et al. Temporary primary frequency control support by variable speed wind turbines—Potential and applications
CN104242356B (en) Consider Robust Interval wind-powered electricity generation dispatching method and the device of wind energy turbine set collection cable malfunction
CN104600742B (en) Method for compensating wind power plant virtual inertia by utilizing energy accumulation device
Gonzalez-Longatt Effects of the synthetic inertia from wind power on the total system inertia: simulation study
Erlich et al. Impact of large wind power generation on frequency stability
Sedaghat et al. Design of a multilevel control strategy for integration of stand-alone wind/diesel system
Gonzalez-Longatt Impact of emulated inertia from wind power on under-frequency protection schemes of future power systems
Bhatt et al. Dynamic participation of doubly fed induction generator in automatic generation control
CN102545268B (en) Large grid active power real-time control method in restricted wind power state
Pulgar Painemal Wind farm model for power system stability analysis
CN106602606A (en) Comprehensive grid frequency modulation control method with consideration of wind power injection
Wang et al. Utilisation of kinetic energy from wind turbine for grid connections: a review paper
CN109728578A (en) Electric system stochastic and dynamic Unit Combination method based on Newton Algorithm quantile
Jietan et al. Dynamic frequency support by photovoltaic generation with “synthetic” inertia and frequency droop control
CN109599865B (en) Method for determining new energy ratio in power system
Krpan et al. Inertial and primary frequency response model of variable‐speed wind turbines
CN114465246B (en) Unit combination optimization method considering speed regulator dead zone and amplitude limiting link
Kamel Employing two novel mechanical fault ride through controllers for keeping stability of fixed speed wind generation systems hosted by standalone micro-grid
CN107346889A (en) Consider that the load of a frequency modulation frequency modulation and minimum frequency deviation cuts down Optimized model construction method
Wang et al. Participation in primary frequency regulation of wind turbines using hybrid control method
Vidyanandan et al. Issues in the grid frequency regulation with increased penetration of wind energy systems
CN109980660A (en) The Robust Interval wind-powered electricity generation dispatching method of frequency response characteristic is improved based on generator

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190312

Termination date: 20200331

CF01 Termination of patent right due to non-payment of annual fee