CN106712058B - The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation - Google Patents
The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation Download PDFInfo
- Publication number
- CN106712058B CN106712058B CN201710055248.1A CN201710055248A CN106712058B CN 106712058 B CN106712058 B CN 106712058B CN 201710055248 A CN201710055248 A CN 201710055248A CN 106712058 B CN106712058 B CN 106712058B
- Authority
- CN
- China
- Prior art keywords
- double
- wind power
- power plant
- fed
- fed blower
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000010248 power generation Methods 0.000 claims abstract description 28
- 238000011156 evaluation Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 238000004891 communication Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 230000005611 electricity Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000007665 sagging Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract 1
- 230000004044 response Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- H02J3/386—
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Abstract
The present invention relates to the control method for coordinating that a kind of double-fed blower wind power plant participates in electric system primary frequency modulation, belong to Operation of Electric Systems and control technology field.This method comprises: establishing the energy state evaluation index and power generation level evaluation index of blower；Set the communication coefficient matrix between blower；Wind power plant cooperation layer measurement frequency simultaneously generates Power Output for Wind Power Field reference value；Adjacent blower exchanges power generation level evaluation index, carries out Newton iteration, the reference value that tracking cooperation layer generates by outer ring alternative manner；In every step iteration, each blower is submitted necessary information by the approximate Newton direction of the inner ring recurrence calculation step for outer ring iterative.The kinetic energy that the method for the present invention makes full use of double-fed blower to store, the coordinated control system of design reduces the computation burden of cooperation layer, and pass through the Quadratic Convergence of distributed Newton method, realize the quick primary frequency modulation of wind power plant entirety, simultaneously, the reasonable distribution that can be realized power between blower, guarantees the safe operation of blower.
Description
Technical field
The present invention relates to the control method for coordinating that a kind of double-fed blower wind power plant participates in electric system primary frequency modulation, belong to electricity
Force system operation and control technology field.
Background technique
Environmental problem drives the innovation and change of electric system, and Renewable Energy Development power generation can optimize energy knot
Structure pushes energy-saving and emission-reduction and realizes sustainable economic development.As one of most representative renewable energy power generation mode, wind
Power, which generates electricity, is widely regarded as effective alternative solution of traditional thermal power generation.The fast development of wind-power electricity generation is gradually changing electricity
Traditional method of operation of Force system.Most obvious problem caused by extensive high permeability wind-power electricity generation accesses be due to wind energy from
The randomness and fluctuation of body will bring serious power swing to electric system.And since large-scale wind power generation is using different
The form of step motor, to reduce the inertia of system, further results in so that its rotor speed is mutually decoupled with system frequency
The fluctuation of frequency.Therefore, in wind-powered electricity generation high permeability area, if the wind-driven generator of large-scale wind power field is all made of maximum power
Tracing control mode will generate huge impact to electric system, or even influence safe and stable operation.
Currently, most of wind power plants use double-fed wind turbine power generation, and the kinetic energy stored in double-fed fan blade and flexible
Control characteristic provides possibility to participate in frequency adjusting.It, can be fast by adjusting the power electronic equipment of double-fed fan rotor side
Speed changes its active output, provides quick primary frequency modulation service using the kinetic energy of blower itself storage.However, most of at present
Double-fed blower participates in the method for primary frequency modulation all only using an individual blower as the unit for participating in frequency modulation, and in electric system
Middle wind power plant answers the frequency of Whole Response system to change, especially in the environment of electricity market wind power plant entirety output power
Even more its standard examined and settled accounts as participant in the market.If but carrying out the coordination that wind power plant participates in power system frequency control
Control generally need to belong to this type using the optimal control method of centralization, current existing method both domestic and external.But due to
Large-scale wind power field blower is large number of, and farther out, centerized fusion needs to build complicated communication network, and serious for geographical distribution
Dependent on wind power plant Centralized Controller.If wind power plant Centralized Controller breaks down, whole system will be at a standstill, therefore reliable
Property it is very low, and concentrate model maintenance and optimization calculate will also take a substantial amount of time, be not appropriate for response speed require very
High primary frequency modulation.
Summary of the invention
The purpose of the present invention is to propose to the coordinated control sides that a kind of double-fed blower wind power plant participates in electric system primary frequency modulation
Method makes it that the frequency of similar synchronous generator be presented can meet the target that wind power plant integrally participates in the adjusting of one secondary frequencies of system
Response characteristic, and can according to the regulating power of different double-fed blowers between double-fed blower reasonable distribution power, guarantee double-fed wind
The safe operation of machine improves the dynamic property of wind power plant primary frequency modulation, and method realizes that conveniently, control is simple.
Double-fed blower wind power plant proposed by the present invention participates in the control method for coordinating of electric system primary frequency modulation, including following
Step:
(1) the energy state evaluation index SOE of double-fed fan blade store kinetic energy in wind power plant is establishedi:
In above formula, SOEiIndicate the energy state evaluation index of i-th double-fed fan blade store kinetic energy in wind power plant,
ωt,iIndicate the revolving speed per unit value of i-th double-fed fan blade in wind power plant, ωtmin,iIndicate i-th double-fed blower in wind power plant
Guarantee the lower limit of the blade rotational speed per unit value of normal power generation, ωtmax,iIndicate that i-th double-fed blower guarantees normal hair in wind power plant
The upper limit of the blade rotational speed per unit value of electricity, i are between 1 to the arbitrary integer between m, and m is the quantity of double-fed blower in wind power plant；
According to the energy state evaluation index SOE of double-fed fan blade store kinetic energy in wind power plant obtained abovei, establish
Double-fed wind turbine power generation assessment of levels index x in wind power planti:
In above formula, Pe,iIndicate the electromagnetic power per unit value of i-th double-fed blower output in wind power plant, Pm,iIndicate wind power plant
In i-th double-fed blower capture mechanical output per unit value；
(2) correspondence set in wind power plant between all double-fed blowers is connection in topology, establishes double-fed blower
Between communication coefficient matrix A=[μij]m×m, element μ in matrixijBetween i-th double-fed blower and jth platform double-fed blower
Communication coefficient:
In above formula, NiIndicate the indexed set of all double-fed blowers for having correspondence with i-th double-fed blower, niAnd nj
Respectively indicate the number for the double-fed blower for having correspondence with i-th and jth platform double-fed blower；
Establish exchange coefficient matrix B=[B between double-fed blowerij]m×m:
B=I-2diag (A)+A
In above formula, I indicate m rank unit matrix, diag (A) indicate be made of the diagonal entry of matrix A to angular moment
Battle array；
(3) wind power plant starts the primary frequency modulation coordinated control to electric system, comprising the following steps:
(3-1) is according to set time period TfThe frequency f for measuring wind power plant and electric system grid entry point, is calculated using following formula
Obtain Power Output for Wind Power Field reference value Pref:
Pref=Pref *+Pf *+Pin *
In above formula, Pref *Indicate Power Output for Wind Power Field instruction value under the stable state that higher level's electric power system dispatching center provides,
Pf *And Pin *The correction amount that is respectively determined by the frequency departure of wind power plant and electric system grid entry point and by wind power plant and power train
The correction amount that the frequency change rate for grid entry point of uniting determines, calculation formula are as follows:
Pf *=Kf(f*-f)
In above formula, KfAnd KinThe respectively sagging coefficient of active power of wind power field-frequency and damped coefficient, KfValue model
It encloses for 0.1~10, KinValue range be 1~20, f*For the rated frequency of electric system；
(3-2) according to setting period of time TPMeasure the reality output active power of wind power plant and electric system grid entry point
Per unit value PPCC, and it is calculate by the following formula wind power plant output work in the reality output active power per unit value and above-mentioned steps (3-1)
Rate reference value PrefDeviation delta PPCC:
ΔPPCC=PPCC-Pref
And by above-mentioned deviation delta PPCCOccurred in a manner of broadcast communication to double-fed blower all in wind power plant；
(4) the double-fed blower of wind power plant starts the primary frequency modulation coordinated control to electric system, comprising the following steps:
When (4-1) is initialized, the double-fed blower starting of wind power plant is set to the primary frequency modulation coordinated control of electric system
Step number k=0 is originated, and makes the clock t for measuring wind power plant and the grid-connected dot frequency of electric systemfWith measurement wind power plant and electric system
The clock t of grid entry point reality output active power per unit valuePTwo clock timing since 0 moment；
The method that (4-2) utilizes above-mentioned steps (1) calculates in wind power plant double-fed blower the of primary frequency modulation coordinated control
The power generation level evaluation index x of k stepi(k):
All double-fed blowers in (4-3) wind power plant and there are other double-fed blowers of correspondence to intercourse by above-mentioned step
Suddenly the double-fed wind turbine power generation assessment of levels standard x for the kth step that (4-2) is calculatedi(k), i.e. the hair that i-th double-fed blower is walked in kth
Electric assessment of levels index is sent to all jth platform double-fed blowers directly connected in communication, wherein j ∈ Ni, and i-th
Platform double-fed blower is collected all jth platform double-fed blowers directly connected in communication and is referred in the power generation level evaluation that kth walks
Mark xi(k)；
Each double-fed blower in (4-4) wind power plant calculates parameter D required for kth step iterates to calculateii,kAnd gi,k:
In above formula, factor alpha is positive real number less than 1, and in one embodiment of the present of invention, the value of α is 0.5.
In above formula, K is weight coefficient, according to wind power plant capacity or the positive real number chosen by experiment, value range are as follows:
0.1~10；
(4-5) utilizes recurrence method, calculates Newton direction recursion of each double-fed blower when kth walks iteration in wind power plant
Value, specific recursive process are as follows:
When (4-5-1) is initialized, if recursion step t=0, and t=T when recursion step terminates is set, T takes 3~5, recursion
Process starts；
(4-5-2) calculates each respective Newton direction recursion initial value of double-fed blower in wind power plant, wherein i-th double-fed wind
The Newton direction recursion initial value of machine is di,k (0):
di,k (0)=Dii,k -1gi,k；
(4-5-3) judges t, if t<T, carries out step (4-5-4), if t>=T, terminates recursion, carries out step
(4-6)；
I-th double-fed blower is walked the Newton direction recursion value d of recursion by (4-5-4) in ti,k (t), it is all to be sent to other
There is the double-fed blower of correspondence, while i-th double-fed blower receives all other double-fed blowers for having correspondence and walk in t
Newton direction recursion value in recursion；
In (4-5-5) wind power plant each double-fed blower according to receive other have the newton side of the double-fed blower of correspondence
To recursion value, the Newton direction recursion value d of+1 step of kth is calculatedi,k (t+1), by taking i-th double-fed blower as an example, calculation formula is as follows:
In above formula, BijThe element in exchange coefficient matrix B between double-fed blower in above-mentioned steps (2)；
(4-5-6) makes t=t+1, return step (4-5-3)；
Each double-fed blower executes iteration in (4-6) wind power plant, calculates i-th double-fed blower of+1 step of kth using following formula
Power generation level evaluation index xi(k+1):
xi(k+1)=xi(k)-εdi,k (T)
In above formula, ε is Newton iteration step-length, is chosen according to the actual capacity of double-fed blower, or adjust value by experiment
Size, to obtain optimum efficiency, value range are as follows: 0.02~1；
Electromagnetic power of each double-fed blower according to new power generation level evaluation index adjustment output, benefit in (4-7) wind power plant
The electromagnetic power per unit value P of i-th double-fed blower output of+1 step of kth is calculated with following formulae,i(k):
Pe,i(k)=SOEi·xi(k)+Pm,i(k)
In above formula, Pm,i(k) the mechanical output per unit value acquired for i-th double-fed blower in kth successive step；
The clock t of (4-8) to measurement wind power plant and electric system grid entry point reality output active power per unit valuePSentenced
It is disconnected, if tP≥TP, make k=k+1, tPRestart timing, and enter step (4-9), if tP<TP, re-execute this step (4-
8)；
The clock t of (4-9) to measurement wind power plant and the grid-connected dot frequency of electric systemfJudged, if tf≥Tf, then return
Step (4-1), if tf<Tf, then return step (4-2).
Double-fed blower wind power plant proposed by the present invention participates in the control method for coordinating of electric system primary frequency modulation, advantage
It is:
1, the method for the present invention establishes the coordinated control frame that wind power plant integrally participates in primary frequency modulation, in blower side group in dilute
The point-to-point control for dredging communication network greatly improves the response speed of wind power plant control, and wind field side cooperation layer only needs to carry out
It measures and information broadcast, the centralized optimization for not needing to carry out complicated model maintenance calculates, burden substantially reduces.
2, this method proposes one based on distributed Newton method with second order convergence characteristic by matrix decomposition technology
Secondary frequency modulation control, this method have it is insensitive to measurement error, can effectively avoid contention control and over control, improve it is primary
The advantages that frequency modulation dynamic characteristic, compared to the distributed frequency modulation algorithm of linear convergence, method convergence proposed by the present invention faster, more
It is suitble to the primary frequency modulation more demanding to response speed.
It 3, can be with so as to easily measure in double-fed fan blade method introduces the concept of energy state index
The kinetic energy of exploitation, enables this method reasonably distribution power between double-fed blower under the premise of not losing safety, and
The wind energy in frequency modulation can be reduced using energy state index reasonable distribution power to lose.
4, the invention proposes a kind of double-fed blower wind power plant participate in electric system primary frequency modulation control method for coordinating,
Wind farm side need to only install tuning controller, and function is only limitted to monitor the frequency and output power of wind farm grid-connected point, generate
Power Output for Wind Power Field reference value, and information is sent to a blower in a broadcast manner.Then by distributed logical between blower
Letter and iteration control respond so that wind field tracks output power reference value.Method proposed by the present invention not only makes wind power plant to being
System is presented whole frequency response characteristic, and can between blower reasonable distribution power, and improve frequency control dynamic characteristic.
Detailed description of the invention
Fig. 1 is primary frequency modulation coordinated control flow diagram of the wind power plant involved in the method for the present invention to electric system.
Fig. 2 is that the generation wind power plant of wind power plant primary frequency modulation coordinated control involved in the method for the present invention is grid-connected to electric system
The control block diagram of point output power reference value.
Fig. 3 is primary frequency modulation coordinated control flow diagram of the double-fed blower that is related to of the method for the present invention to electric system.
Specific embodiment
Double-fed blower wind power plant proposed by the present invention participates in the control method for coordinating of electric system primary frequency modulation, including two
Control flow, i.e. wind power plant tuning controller control flow, double-fed blower distributed AC servo system process, two processes are according to certain
Loop cycle executes, and the information of station layer is periodically sent to double-fed by way of broadcast communication by wind power plant tuning controller
Controller of fan, this method specifically includes the following steps:
(1) the energy state evaluation index SOE of double-fed fan blade store kinetic energy in wind power plant is establishedi:
In above formula, SOEiIndicate the energy state evaluation index of i-th double-fed fan blade store kinetic energy in wind power plant,
ωt,iIndicate the revolving speed per unit value of i-th double-fed fan blade in wind power plant, ωtmin,iIndicate i-th double-fed blower in wind power plant
Guarantee the lower limit of the blade rotational speed per unit value of normal power generation, ωtmax,iIndicate that i-th double-fed blower guarantees normal hair in wind power plant
The upper limit of the blade rotational speed per unit value of electricity, i are between 1 to the arbitrary integer between m, and m is the quantity of double-fed blower in wind power plant；
According to the energy state evaluation index SOE of double-fed fan blade store kinetic energy in wind power plant obtained abovei, establish in wind power plant
Double-fed wind turbine power generation assessment of levels index xi:
In above formula, Pe,iIndicate the electromagnetic power per unit value of i-th double-fed blower output in wind power plant, Pm,iIndicate wind power plant
In i-th double-fed blower capture mechanical output per unit value；
(2) correspondence set in wind power plant between all double-fed blowers is connection in topology, establishes double-fed blower
Between communication coefficient matrix A=[μij]m×m, element μ in matrixijBetween i-th double-fed blower and jth platform double-fed blower
Communication coefficient:
In above formula, NiIndicate the indexed set of all double-fed blowers for having correspondence with i-th double-fed blower, niAnd nj
Respectively indicate the number for the double-fed blower for having correspondence with i-th and jth platform double-fed blower；
Establish exchange coefficient matrix B=[B between double-fed blowerij]m×m:
B=I-2diag (A)+A
In above formula, I indicate m rank unit matrix, diag (A) indicate be made of the diagonal entry of matrix A to angular moment
Battle array；
(3) wind power plant starts primary frequency modulation coordinated control to electric system, flow diagram as shown in Figure 1, include with
Lower step:
(3-1) is according to set time period Tf(recommended value 0.5 second, can be adjusted according to wind power plant suitable scale) measurement wind-powered electricity generation
The frequency f of field and electric system grid entry point, is calculated Power Output for Wind Power Field reference value P using following formularef:
Pref=Pref *+Pf *+Pin *
In above formula, Pref *Indicate Power Output for Wind Power Field instruction value under the stable state that higher level's electric power system dispatching center provides,
Pf *And Pin *The correction amount that is respectively determined by the frequency departure of wind power plant and electric system grid entry point and by wind power plant and power train
The correction amount that the frequency change rate for grid entry point of uniting determines, calculation formula are as follows:
Pf *=Kf(f*-f)
In above formula, KfAnd KinThe respectively sagging coefficient of active power of wind power field-frequency and damped coefficient, KfValue model
It encloses for 0.1~10, KinValue range be 1~20, f*For the rated frequency of electric system, in one embodiment of the present of invention,
KfAnd KinRespectively 0.8 and 2；
(3-2) according to setting period of time TPIt measures within (recommended value 0.08 second, can be adjusted according to wind power plant suitable scale)
The reality output active power per unit value P of wind power plant and electric system grid entry pointPCC, and be calculate by the following formula the reality output and have
Function power per unit value and Power Output for Wind Power Field reference value P in above-mentioned steps (3-1)refDeviation delta PPCC:
ΔPPCC=PPCC-Pref
And by above-mentioned deviation delta PPCCOccurred in a manner of broadcast communication to double-fed blower all in wind power plant；
(4) the double-fed blower of wind power plant starts the primary frequency modulation coordinated control to electric system, flow diagram such as Fig. 3 institute
Show, comprising the following steps:
When (4-1) is initialized, the double-fed blower starting of wind power plant is set to the primary frequency modulation coordinated control of electric system
Step number k=0 is originated, and makes tfAnd tPTwo clock timing since 0 moment；
The method that (4-2) utilizes above-mentioned steps (1) calculates in wind power plant double-fed blower the of primary frequency modulation coordinated control
The power generation level evaluation index x of k stepi(k):
All double-fed blowers in (4-3) wind power plant and there are other double-fed blowers of correspondence to intercourse by above-mentioned step
Suddenly the double-fed wind turbine power generation assessment of levels standard x for the kth step that (4-2) is calculatedi(k), i.e. the hair that i-th double-fed blower is walked in kth
Electric assessment of levels index is sent to all jth platform double-fed blowers directly connected in communication, wherein j ∈ Ni, and i-th
Platform double-fed blower is collected all jth platform double-fed blowers directly connected in communication and is referred in the power generation level evaluation that kth walks
Mark xi(k)；
Each double-fed blower in (4-4) wind power plant calculates parameter D required for kth step iterates to calculateIi, kAnd gI, k:
In above formula, factor alpha is positive real number less than 1, and in one embodiment of the present of invention, the value of α is 0.5.
In above formula, K is weight coefficient, according to wind power plant capacity or the positive real number chosen by experiment, value range are as follows:
0.1~10；
(4-5) utilizes recurrence method, calculates Newton direction recursion of each double-fed blower when kth walks iteration in wind power plant
Value, specific recursive process are as follows:
When (4-5-1) is initialized, if recursion step t=0, and t=T when recursion step terminates is set, T takes 3~5, recursion
Process starts；
(4-5-2) calculates each respective Newton direction recursion initial value of double-fed blower in wind power plant, wherein i-th double-fed wind
The Newton direction recursion initial value of machine is di,k (0):
di,k (0)=Dii,k -1gi,k；
(4-5-3) judges t, if t<T, carries out step (4-5-4), if t>=T, terminates recursion, carries out step
(4-6)；
I-th double-fed blower is walked the Newton direction recursion value d of recursion by (4-5-4) in ti,k (t), it is all to be sent to other
There is the double-fed blower of correspondence, while i-th double-fed blower receives all other double-fed blowers for having correspondence and walk in t
Newton direction recursion value in recursion；
In (4-5-5) wind power plant each double-fed blower according to receive other have the newton side of the double-fed blower of correspondence
To recursion value, the Newton direction recursion value d of+1 step of kth is calculatedi,k (t+1), by taking i-th double-fed blower as an example, calculation formula is as follows:
In above formula, BijThe element in exchange coefficient matrix B between double-fed blower in above-mentioned steps (2)；
(4-5-6) makes t=t+1, return step (4-5-3)；
Each double-fed blower executes iteration in (4-6) wind power plant, calculates i-th double-fed blower of+1 step of kth using following formula
Power generation level evaluation index xi(k+1):
xi(k+1)=xi(k)-εdi,k (T)
In above formula, ε is Newton iteration step-length, is chosen according to the actual capacity of double-fed blower, or adjust value by experiment
Size, to obtain optimum efficiency, value range are as follows: 0.02~1；
Electromagnetic power of each double-fed blower according to new power generation level evaluation index adjustment output, benefit in (4-7) wind power plant
The electromagnetic power per unit value P of i-th double-fed blower output of+1 step of kth is calculated with following formulae,i(k):
Pe,i(k)=SOEi·xi(k)+Pm,i(k)
In above formula, Pm,i(k) the mechanical output per unit value acquired for i-th double-fed blower in kth successive step；
The time t of (4-8) to measurement wind power plant and electric system grid entry point reality output active power per unit valuePSentenced
It is disconnected, if tP≥TP, make k=k+1, tPRestart timing, and enter step (4-9), if tP<TP, re-execute this step (4-
8)；
The time t of (4-9) to measurement wind power plant and the grid-connected dot frequency of electric systemfJudged, if tf≥Tf, then return
Step (4-1), if tf<Tf, then return step (4-2).
Claims (1)
1. the control method for coordinating that a kind of double-fed blower wind power plant participates in electric system primary frequency modulation, it is characterised in that this method packet
Include following steps:
(1) the energy state evaluation index SOE of double-fed fan blade store kinetic energy in wind power plant is establishedi:
In above formula, SOEiIndicate the energy state evaluation index of i-th double-fed fan blade store kinetic energy in wind power plant, ωt,iTable
Show the revolving speed per unit value of i-th double-fed fan blade in wind power plant, ωtmin,iIndicate that i-th double-fed blower guarantees just in wind power plant
The lower limit of the blade rotational speed per unit value often to generate electricity, ωtmax,iIndicate the leaf of i-th double-fed blower guarantee normal power generation in wind power plant
The upper limit of piece revolving speed per unit value, i are between 1 to the arbitrary integer between m, and m is the quantity of double-fed blower in wind power plant；
According to the energy state evaluation index SOE of double-fed fan blade store kinetic energy in obtained wind power planti, establish in wind power plant
Double-fed wind turbine power generation assessment of levels index xi:
In above formula, Pe,iIndicate the electromagnetic power per unit value of i-th double-fed blower output in wind power plant, Pm,iIt indicates the in wind power plant
The mechanical output per unit value of i platform double-fed blower capture；
(2) correspondence set in wind power plant between all double-fed blowers is connection in topology, is established between double-fed blower
Communication coefficient matrix A=[μij]m×m, element μ in matrixijIt is logical between i-th double-fed blower and jth platform double-fed blower
Believe coefficient:
In above formula, NiIndicate the indexed set of all double-fed blowers for having correspondence with i-th double-fed blower, niAnd njRespectively
Indicate the number for the double-fed blower for having correspondence with i-th and jth platform double-fed blower；
Establish exchange coefficient matrix B=[B between double-fed blowerij]m×m:
B=I-2diag (A)+A
In above formula, I indicates that the unit matrix of m rank, diag (A) indicate the diagonal matrix being made of the diagonal entry of matrix A；
(3) wind power plant starts the primary frequency modulation coordinated control to electric system, comprising the following steps:
(3-1) is according to set time period TfThe frequency f for measuring wind power plant and electric system grid entry point, is calculated using following formula
Power Output for Wind Power Field reference value Pref:
Pref=Pref *+Pf *+Pin *
In above formula, Pref *Indicate Power Output for Wind Power Field instruction value under the stable state that higher level's electric power system dispatching center provides, Pf *With
Pin *The correction amount that is respectively determined by the frequency departure of wind power plant and electric system grid entry point and by wind power plant with electric system simultaneously
The correction amount that the frequency change rate of site determines, calculation formula are as follows:
Pf *=Kf(f*-f)
In above formula, KfAnd KinThe respectively sagging coefficient of active power of wind power field-frequency and damped coefficient, KfValue range be
0.1~10, KinValue range be 1~20, f*For the rated frequency of electric system；
(3-2) according to setting period of time TPMeasure the reality output active power mark of wind power plant and electric system grid entry point
Value PPCC, and it is calculate by the following formula Power Output for Wind Power Field reference value in the reality output active power per unit value and step (3-1)
PrefDeviation delta PPCC:
ΔPPCC=PPCC-Pref
And by above-mentioned deviation delta PPCCOccurred in a manner of broadcast communication to double-fed blower all in wind power plant；
(4) the double-fed blower of wind power plant starts the primary frequency modulation coordinated control to electric system, comprising the following steps:
When (4-1) is initialized, starting of the double-fed blower starting of wind power plant to the primary frequency modulation coordinated control of electric system is set
Step number k=0, and make the clock t for measuring wind power plant and the grid-connected dot frequency of electric systemfIt is grid-connected with measurement wind power plant and electric system
The clock t of point reality output active power per unit valuePTwo clock timing since 0 moment；
The method that (4-2) utilizes step (1) calculates the hair that double-fed blower is walked in the kth of primary frequency modulation coordinated control in wind power plant
Electric assessment of levels index xi(k):
All double-fed blowers in (4-3) wind power plant and there are other double-fed blowers of correspondence to intercourse by step (4-2)
The double-fed wind turbine power generation assessment of levels standard x of the kth step of calculatingi(k), i.e., i-th double-fed blower is commented in the power generation level that kth walks
Valence index is sent to all jth platform double-fed blowers directly connected in communication, wherein j ∈ Ni, and i-th double-fed wind
Machine collects the power generation level evaluation index x that all jth platform double-fed blowers directly connected in communication are walked in kthi(k)；
Each double-fed blower in (4-4) wind power plant calculates parameter D required for kth step iterates to calculateii,kAnd gi,k:
In above formula, factor alpha is positive real number less than 1, and in one embodiment of the present of invention, the value of α is 0.5；
In above formula, K is weight coefficient, is chosen according to wind power plant capacity or by experiment, and the K is positive real number, value range are as follows:
0.1~10；
(4-5) utilizes recurrence method, calculates Newton direction recursion value of each double-fed blower when kth walks iteration in wind power plant, tool
Body recursive process are as follows:
When (4-5-1) is initialized, if recursion step t=0, and t=T when recursion step terminates is set, T takes 3~5, recursive process
Start；
(4-5-2) calculates each respective Newton direction recursion initial value of double-fed blower in wind power plant, wherein i-th double-fed blower
Newton direction recursion initial value is di,k (0):
di,k (0)=Dii,k -1gi,k；
(4-5-3) judges t, if t < T, carries out step (4-5-4), if t >=T, terminates recursion, carries out step (4-
6)；
I-th double-fed blower is walked the Newton direction recursion value d of recursion by (4-5-4) in ti,k (t), be sent to other it is all have it is logical
The double-fed blower of gateway system, while i-th double-fed blower receives all other double-fed blowers for having correspondence and walks recursion in t
In Newton direction recursion value；
In (4-5-5) wind power plant each double-fed blower according to receive other there is the Newton direction of the double-fed blower of correspondence to pass
Value is pushed away, the Newton direction recursion value d of+1 step of kth is calculatedi,k (t+1), by taking i-th double-fed blower as an example, calculation formula is as follows:
In above formula, BijThe element in exchange coefficient matrix B between double-fed blower in step (2)；
(4-5-6) makes t=t+1, return step (4-5-3)；
Each double-fed blower executes iteration in (4-6) wind power plant, and the power generation of i-th double-fed blower of+1 step of kth is calculated using following formula
Assessment of levels index xi(k+1):
xi(k+1)=xi(k)-εdi,k (T)
In above formula, ε is Newton iteration step-length, is chosen according to the actual capacity of double-fed blower, or adjust value size by experiment,
Value range are as follows: 0.02~1；
The electromagnetic power that each double-fed blower is exported according to new power generation level evaluation index adjustment in (4-7) wind power plant, under utilization
Formula calculates the electromagnetic power per unit value P of i-th double-fed blower output of+1 step of kthe,i(k):
Pe,i(k)=SOEi·xi(k)+Pm,i(k)
In above formula, Pm,i(k) the mechanical output per unit value acquired for i-th double-fed blower in kth successive step；
The clock t of (4-8) to measurement wind power plant and electric system grid entry point reality output active power per unit valuePJudged, if
tP≥TP, make k=k+1, tPRestart timing, and enter step (4-9), if tP< TP, re-execute this step (4-8)；
The clock t of (4-9) to measurement wind power plant and the grid-connected dot frequency of electric systemfJudged, if tf≥Tf, then return step
(4-1), if tf< Tf, then return step (4-2).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710055248.1A CN106712058B (en) | 2017-01-24 | 2017-01-24 | The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710055248.1A CN106712058B (en) | 2017-01-24 | 2017-01-24 | The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106712058A CN106712058A (en) | 2017-05-24 |
| CN106712058B true CN106712058B (en) | 2019-05-21 |
Family
ID=58909668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710055248.1A Active CN106712058B (en) | 2017-01-24 | 2017-01-24 | The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106712058B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107171368B (en) * | 2017-07-19 | 2020-12-08 | 国家电网公司 | Wind power generation primary frequency modulation function implementation method based on wind power plant power control |
| CN108110803B (en) * | 2018-01-23 | 2020-07-28 | 西南交通大学 | Method for coordinating and controlling double-fed fan auxiliary synchronous generator to participate in secondary frequency modulation of power grid |
| CN109038614B (en) * | 2018-07-12 | 2022-10-04 | 中国电力科学研究院有限公司 | Control method and system for photovoltaic power station to participate in system frequency modulation |
| CN110401222B (en) * | 2019-07-30 | 2020-10-27 | 北京理工大学 | Comprehensive control method and system for wind generating set participating in system frequency modulation |
| CN110783953B (en) * | 2019-10-31 | 2023-11-14 | 特变电工西安电气科技有限公司 | Quick frequency-modulation response communication method for string-type inverter photovoltaic power station |
| CN111523818B (en) * | 2020-04-28 | 2022-12-09 | 国网江西省电力有限公司电力科学研究院 | Thermal power generating unit primary frequency modulation performance evaluation method based on test information |
| CN114256857B (en) * | 2021-12-07 | 2024-03-26 | 国网河南省电力公司电力科学研究院 | Control method and system for improving frequency modulation capability of doubly-fed wind turbine generator |
| CN114726005B (en) * | 2021-12-31 | 2023-06-16 | 国网安徽省电力有限公司安庆供电公司 | Wind power plant double-layer frequency control method considering fan power optimization distribution |
| CN116014819B (en) * | 2023-03-23 | 2023-06-06 | 湖南大学 | Active power storage limit quantification method and system of permanent magnet wind generating set |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103795089A (en) * | 2013-12-12 | 2014-05-14 | 国家电网公司 | Over-speed and variable pitch coordination-based variable speed wind turbine generator set primary frequency modulation method |
| CN105207198A (en) * | 2015-09-06 | 2015-12-30 | 国家电网公司 | Primary frequency modulation control system based on frequency droop regulation method and method thereof |
| CN105826936A (en) * | 2016-05-06 | 2016-08-03 | 上海明华电力技术工程有限公司 | Thermal power generating unit intelligent primary frequency modulation control method aiming at large frequency-difference of power grid |
-
2017
- 2017-01-24 CN CN201710055248.1A patent/CN106712058B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103795089A (en) * | 2013-12-12 | 2014-05-14 | 国家电网公司 | Over-speed and variable pitch coordination-based variable speed wind turbine generator set primary frequency modulation method |
| CN105207198A (en) * | 2015-09-06 | 2015-12-30 | 国家电网公司 | Primary frequency modulation control system based on frequency droop regulation method and method thereof |
| CN105826936A (en) * | 2016-05-06 | 2016-08-03 | 上海明华电力技术工程有限公司 | Thermal power generating unit intelligent primary frequency modulation control method aiming at large frequency-difference of power grid |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106712058A (en) | 2017-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106712058B (en) | The control method for coordinating of double-fed blower wind power plant participation electric system primary frequency modulation | |
| CN102738809B (en) | Optimized control method for wind power field reactive power compensation capacity considering wind power distribution rule | |
| CN107565611B (en) | A kind of method of wind power plant inertia frequency modulation | |
| CN105140936B (en) | A kind of micro-capacitance sensor frequency modulation control method for minimizing cost of electricity-generating | |
| CN107069829B (en) | Control system, method and application of station-level virtual synchronous machine | |
| CN104578086B (en) | A kind of reactive voltage control method of wind-powered electricity generation and photovoltaic generation access power grid | |
| CN102684201B (en) | Voltage threshold probability-based reactive power optimizing method for grid containing wind power plant | |
| CN103605360B (en) | A kind of test system and method for wind farm power control strategy | |
| CN105262098B (en) | The quick automatic voltage control method assessed based on the fluctuation of wind power plant generated output | |
| CN102354992A (en) | Reactive-power control method of wind power field | |
| CN106160000B (en) | The wind-powered electricity generation that a kind of point of tune-province's tune two-stage is coordinated issues additional control method | |
| CN103762919B (en) | Output control device and method for direct-drive aerogenerator low voltage crossing | |
| CN107482649A (en) | A kind of two domain interacted system LOAD FREQUENCY control methods based on frequency dividing control | |
| CN109494811A (en) | A kind of wind power plant unit participates in the Poewr control method and system of frequency modulation and voltage modulation | |
| CN106058937A (en) | Generalized load modeling method for distribution network containing direct-drive permanent magnet wind farm | |
| CN108539760A (en) | A kind of double-fed induction Wind turbines frequency modulation PID control method based on group's grey wolf optimization algorithm | |
| CN110206686A (en) | A kind of adaptive maximum power tracking and controlling method for wind power generating set | |
| CN109842151A (en) | A kind of control method of hydroenergy storage station wide area power-generating control system | |
| Guo et al. | Double-layer feedback control method for synchronized frequency regulation of PMSG-based wind farm | |
| CN105470978B (en) | A kind of static passive compensation device control method for coordinating in groups | |
| CN108365627A (en) | A kind of wind storage isolated network power supply system control method for coordinating based on flexible coordinating factor | |
| CN105720573A (en) | Actually-measured data based modeling method for active power and reactive power control system of wind-light power storage station | |
| CN110854907A (en) | Collaborative optimization operation method and system for power distribution network wind power plant under communication fault | |
| CN115395587A (en) | Multi-source collaborative offshore wind farm reactive voltage control method | |
| CN111049178A (en) | Method for analyzing stability control of direct-drive permanent magnet wind turbine generator through VSC-HVDC grid connection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |