CN105909470B - The adaptive maximum power tracking and controlling method of wind power generating set - Google Patents

The adaptive maximum power tracking and controlling method of wind power generating set Download PDF

Info

Publication number
CN105909470B
CN105909470B CN201610235513.XA CN201610235513A CN105909470B CN 105909470 B CN105909470 B CN 105909470B CN 201610235513 A CN201610235513 A CN 201610235513A CN 105909470 B CN105909470 B CN 105909470B
Authority
CN
China
Prior art keywords
wind
wind speed
speed
value
generating set
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
Application number
CN201610235513.XA
Other languages
Chinese (zh)
Other versions
CN105909470A (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.)
SHANGHAI CHINESE CAR RIBERD INTELLIGENT SYSTEM Co.,Ltd.
Original Assignee
Shanghai Ruibode Intelligent System Ltd By Share 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 Shanghai Ruibode Intelligent System Ltd By Share Ltd filed Critical Shanghai Ruibode Intelligent System Ltd By Share Ltd
Priority to CN201610235513.XA priority Critical patent/CN105909470B/en
Publication of CN105909470A publication Critical patent/CN105909470A/en
Application granted granted Critical
Publication of CN105909470B publication Critical patent/CN105909470B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/06—Controlling wind motors  the wind motors having rotation axis substantially perpendicular to the air flow entering the rotor
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00—Control
    • F05B2270/10—Purpose of the control system
    • F05B2270/103—Purpose of the control system to affect the output of the engine
    • F05B2270/1032—Torque
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00—Control
    • F05B2270/10—Purpose of the control system
    • F05B2270/103—Purpose of the control system to affect the output of the engine
    • F05B2270/1033—Power (if explicitly mentioned)
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00—Control
    • F05B2270/30—Control parameters, e.g. input parameters
    • F05B2270/32—Wind speeds
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00—Control
    • F05B2270/30—Control parameters, e.g. input parameters
    • F05B2270/335—Output power or torque
    • 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/74—Wind turbines with rotation axis perpendicular to the wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Development (AREA)
  • Sustainable Energy (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Wind Motors (AREA)

Abstract

A kind of adaptive maximum power tracking and controlling method of wind power generating set, wind power generating set to be measured is set first and refers to tip-speed ratio and reference power coefficient, to obtain reference rotation velocity torque aircraft pursuit course and its deformation curve according to reference to tip-speed ratio and reference power coefficient；Then the interval division from incision wind speed to rated wind speed at small wind speed interval, in every each small wind speed interval, wind power generating set tracking is often run and is measured to result datas such as output powers in the period after a cycle and statistical average by deformation coefficient corresponding deformation curve operation.To which in continuous service, using deformation coefficient as optimized variable, output power is optimization aim, and optimal deformation coefficient is found by optimal method so that wind driven generator output power is maximum.Finally, after wind power generating set long-play, the deformation curve segmentation connection that each wind speed section deformation coefficient determines is obtained with the optimized rotating speed torque aircraft pursuit course in whole service wind speed range.

Description

The adaptive maximum power tracking and controlling method of wind power generating set
Technical field
Measure and calculate the present invention relates to physical field more particularly to wind-force, especially a kind of wind power generating set from Adapt to maximum power tracking and controlling method.
Background technology
H-type darrieus perpendicular shaft configuration form, wind wheel part-structure such as Fig. 1 and Fig. 2 institutes may be used in wind-driven generator Show.It is 2~4 that it, which uses H-type fixed-wing prismatic blade 2,2 numbers of blade,；Its each blade 2 and the angle of attack of direction of rotation are to fix Value α.3 connecting rod of horizontal cross is changed into the aerodynamic force on H-type prismatic blade to act on the generator amature of wind wheel center and turn Square.The ratio for radially projecting area of the blade 3 on blade geometry center rotation swept surface being accounted for entire swept surface is defined as solidity SD, solidity SD≤0.3 of general wind power generating set.The radius for scanning circle (cylinder) of wind wheel blade center line 4 is R, blade Height is H, and the planimetric area of swept surface is：A=2RH
The ratio of wind generator set blade linear velocity and wind speed is defined as tip-speed ratio TSR；Wind power generating set is transported Sector-style speed is ranging from from incision wind speed to rated wind speed, then arrives cut-out wind speed, such as typical data has：Incision wind speed is 3m/s, Rated wind speed is 12m/s, cut-out wind speed 20m/s.Wind power generating set at any constant wind speed V under cut-out wind speed, Relation curve between wind speed round and output machine torque has following features：
1. at a certain tip-speed ratio TSRTmax, i.e., when wind speed round is NTmax=9.55V/RTSRTmax, wind wheel Output machine torque reaches maximum Tmax, which is known as maximum torque point；
2. in tip-speed ratio between 0~TSRTmax, i.e., wind speed round is located at 0~NTmax, and wind wheel exports machine torque There are a minimum point (TSRTmin, Tmin), general point Tmin ＜ Tmax/4；The corresponding wind speed rounds of TSRTmin are NTmin=9.55V/RTSRTmin；
It is N0 3. rotating speed reaches highest when wind wheel zero load, N0 is more than NPmax；
4. in a certain tip-speed ratio TSRPmax, i.e., wind speed round is NPmax=9.55V/RTSRPmax, wind wheel output Power reaches maximum Pmax；NPmax ＞ NTmax, the operating point are known as maximum power point, and tip-speed ratio at this time is properly termed as optimal Tip-speed ratio；
5. between NTmin and NTmax, wind wheel output torque is increased monotonically；
6. between NTmax and N0, wind wheel output torque dullness reduces；
The relation curve schematic diagram of above-mentioned rotating speed N and mechanical output T is shown in Fig. 3, referred to as speed torque characteristic curve.
When constant wind speed takes different value, one group of torque and rotation speed characteristic can be obtained.Such as shown in Fig. 4, expression One group of speed torque characteristic curve when be different wind speed being respectively V1, V2, V3, V4 (V1 ＜ V2 ＜ V3 ＜ V4).Each constant wind The maximum power point on speed torque characteristic curve under speed links up, so that it may to obtain optimized rotating speed torque aircraft pursuit course. Under rated wind speed, if the generator electromagnetic braking torque energy of wind power generating set tracks rapidly the real-time rotating speed measured in real time Torque functional value on the optimized rotating speed torque aircraft pursuit course, then the wind power generating set can export machine as much as possible Tool power, and then electrical power is sent out as much as possible.Assuming that V3 is rated wind speed, then the optimized rotating speed torque under rated power Aircraft pursuit course is as shown in Fig. 4 heavy line parts.On rated wind speed, if continuing to track optimized rotating speed torque aircraft pursuit course, The electromagnetic braking torque that generator provides can then be resulted in the need for and be more than its design maximum, cause wind power generating set uncontrolled. Therefore, optimized rotating speed torque aircraft pursuit course is no longer tracked on rated wind speed, but by keeping or reducing wind-power electricity generation The rotating speed of unit so that wind power generating set is not operate near optimal tip-speed ratio, and blade is in stall fortune in the most of the time Row state, the energy that wind power generating set is obtained from wind reduce, can be using speed torque characteristic under constant wind speed in the slow-speed of revolution The characteristics of lower torque drop so that the torque of wind power generating set can be maintained at lower value when wind speed is larger, and should be compared with Low value is less than the design maximum torque value of generator, to ensure that wind power generating set can cut out wind on rated wind speed Safe operation under speed.
In order under rated wind speed wind power generating set run on the state of maximum power output, need to obtain most Excellent rotational speed and torque aircraft pursuit course；And according to above-mentioned analysis optimized rotating speed torque aircraft pursuit course in order to obtain, it is necessary to know each Speed torque characteristic curve under constant wind speed.
In fact, after a vertical axis aerogenerator group design and manufacture are completed, turn under each of which constant wind speed Fast torque characteristics and optimized rotating speed torque aircraft pursuit course are the inherent characteristics of physical presence, only cannot be only by design and meter It calculates accurately and really obtains, need to obtain by practical measurement and verify design value.And in practice, in the natural environment There is no constant wind speed and direction, cannot directly measure the speed torque characteristic under constant wind speed yet.It, can be in wind-tunnel environment Artificial manufacture constant wind speed, but influenced by wind-tunnel size and bolus effect, larger vertical axis aerogenerator group Accurate result can not be obtained by wind tunnel test.
The infield of the same model of wind power generating set is usually not secured to a place, and environment is also with the time Variation, therefore the speed torque characteristic curve and optimized rotating speed torque aircraft pursuit course that are obtained in specific place specific environment are by difference The influence of weather and landform can change, therefore track a fixed optimized rotating speed torque aircraft pursuit course and may result in this Model wind power generating set different from experimental field with cannot maximally send out electric energy under Different climate environmental condition.
Therefore it needs to find a kind of method, when enabling to wind power generating set that can run in naturally practical wind field certainly It adaptively tracks its optimized rotating speed torque aircraft pursuit course and obtains the optimized rotating speed torque aircraft pursuit course.
Invention content
The present invention is to provide for a kind of maximum power tracking and controlling method of wind power generating set so that wind power generating set It can be tracked to runtime adaptable in naturally practical wind field and obtain optimized rotating speed torque aircraft pursuit course.
The technical scheme is that using following steps method：
Step (1)
Wind power generating set needs to optimize the wind speed section of output power generally from incision wind speed V1 to rated wind speed V2.
The wind speed minizone of 0.5m/s width is divided into from total wind speed interval of incision wind speed V1 rated wind speeds V2, those Wind speed minizone label from low to high.
In this way, IMAX 0.5m/s width can be divided into from incision wind speed V1 to wind speed interval total rated wind speed V2 Wind speed minizone, IMAX=int (2V2+0.5)-int (2V1+0.5)+1, int () are bracket function；
After wind power generating set runs a cycle T1, if the mean wind speed in its T1 period is VSX, then the operation is all Phase, corresponding wind speed minizone label IX was obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
In order to store the data for corresponding to each wind speed minizone, the array that following width is IMAX is opened up；
KX () --- storage corresponds to the deformation coefficient of each wind speed minizone, initial value 1.0；
VSX () --- storage corresponds to the mean wind speed value of each wind speed minizone, initial value 0；
PJSXZ () --- storage corresponds to the net power output value after the conversion of each wind speed minizone, initial value 0；
NSXZ () --- storage corresponds to the reduced value of the mean speed value of each wind speed minizone, initial value 0；
NN () --- storage corresponds to the optimization number of each wind speed minizone, initial value 0；
FB () --- storage corresponds to the optimization perturbation direction mark of each wind speed minizone, initial value 1；
Initial deformation COEFFICIENT K XP=1.0 is set, and goes to next step.
Step (2)
It is according to the deformation curve for being worth to reference rotation velocity torque aircraft pursuit course of KXP：
T=CK (KXPN)2 (1)
Wherein coefficient
CK=6.716 × 10-3·A·(R/TSR)3·CP (2)
R scans radius surface for wind wheel blade；
Projected areas of the A for wind wheel blade swept surface on carrying out wind direction vertical plane；
TSR is with reference to tip-speed ratio, and value range is：2.0~3.0
CP is reference power coefficient, and value range is：0.1~0.25
Tip-speed ratio refers to the ratio of linear velocity and wind speed on the outside of wind-driven generator blade of wind-driven generator；Power coefficient refers to wind The ratio that the net electromotive power output of power generator group passes through the wind energy in rotor swept face with free stream.Net electromotive power output, letter It refers to the electric energy of charge controller output end for off-network type wind generator system to claim output power；Grid-connected wind is sent out Refer to the electric energy that current transformer output end is connected to the grid for electric system；
The derivation of formula (1) is：
At a certain wind speed V, when tip-speed ratio is TSR, the rotating speed of wind-driven generator group wind-wheel is set as N, then
N=VTSR/R9.55 (3)
V=NR/ (9.55TSR) (4)
At a certain wind speed V, when power coefficient is CP, the output power of wind power generating set is：
P=0.5 × 1.225AV3·CP (5)
Its corresponding torque is：T=P/N9.55 (6)
Formula (5) is substituted into according to formula (4), formula (5) is then substituted into (6) again, can be obtained：
T=6.716 × 10-3A·(R/TSR)3·CP·N3/ N=CKN2, (7)
Wherein CK=6.716 × 10-3·A·(R/TSR)3·CP (2)
The represented of formula (2.7) is a function, is obtained by being multiplied by a COEFFICIENT K XP to independent variable N with reference to turning The deformation curve of fast torque aircraft pursuit course, the deformation curve are expressed as
T=CK (KXPN)2 (1)
KXP is a variable, referred to as deformation coefficient, and value range is 0.5~2.0；
It is the derivation of formula (1) above.
When wind power generating set is run, the wind speed round N of wind power generating set is measured in real time, is then obtained by formula (1) To torque instruction value T, then controls wind power generating set reality output torque TL and track torque instruction value T operations.
Go to next step；
Step (3)
Wind power generating set by 2 periodic measurement of cycle T and obtains the T2 by step (2) method control operation T1 periods The real-time air speed value in period, real-time tachometer value and real-time output power value, T1=mT2, m are no less than 50 positive integer, T1> >T2；Wind-driven generator runs the T1 periods, obtains the real-time air speed value in m T2 period, real-time tachometer value and real-time output power value Average value VSX, NSX and PSX, and obtain the real-time tachometer value NS in first T2 period and the T1 periods in the T1 periods simultaneously The real-time tachometer value NE in the last one interior T2 period；
Then, it is calculate by the following formula to obtain the kinetic energy delivered power PDSX in the T1 periods：
PDSX=0.5J ((NS/9.55)2-(NE/9.55)2)
Wherein J is the rotary inertia of wind wheel rotating part；
The output power value discharged by wind wheel rotating part kinetic energy change in the T1 periods, is defined as kinetic energy delivered power PDSX。
It is defined on the average net power output PJSX that wind wheel in T1 is absorbed from wind；
It is calculate by the following formula to obtain the T again1The average net power output PJSX that wind wheel in period is absorbed from wind：
PJSX=PSX-PDSX
It is calculate by the following formula to obtain the T1The wind speed intermediate value in the affiliated 0.5m/s wind speed width sections period mean wind speed VSX VS：
VS=int (2VSX+0.5)/2
Int () is bracket function；
PJSX is converted by following formula by average net power output PJSXZP corresponding to wind speed intermediate value VS again：
PJSXZP=PJSX (VS/VSX)3
NSX is converted average net power output NSXZP corresponding to wind speed intermediate value VS by following formula：
NSXZP=NSX (VS/VSX)
For storing the array label IX for representing different wind speed intervals, can be obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
Go to next step；
Step (4)
Step (5)
Step (6)
Step (7)
If the total continuous operating time of wind-driven generator is more than certain setting value, i.e., by long lasting for continuing to optimize, The deformation coefficient and operation mean speed of each wind speed minizone of optimization are stored in KX (IX) and NSXZ (IX), under TX (IX) is calculated in formula：
TX (IX)=CK (KX (IX) NSXZ (IX))2, IX=1,2 ..., IMAX
So as to obtain one group of data point (NSXZ (IX), TX (IX))
The broken line that the data point connects is optimized rotating speed torque aircraft pursuit course；
Go to step (2) circular flow.
The beneficial effects of the invention are as follows：The method of the present invention enables to wind power generating set can be in naturally practical wind field Track and obtain to middle runtime adaptable its optimized rotating speed torque aircraft pursuit course.
Description of the drawings
Fig. 1 is the vertical view of vertical axis aerogenerator group wind wheel part；
Fig. 2 is the side view of vertical axis aerogenerator group wind wheel part；
Fig. 3 is the speed torque characteristic curve graph under any constant wind speed；
Fig. 4 is speed torque characteristic curve and optimized rotating speed torque trace plot under different wind speed；
Fig. 5 is the schematic diagram that optimized rotating speed torque aircraft pursuit course is fitted with deformation curve；
Fig. 6 is reference rotation velocity torque trace plot；
Fig. 7 is the deformation curve of one group of reference rotation velocity torque curve for being fitted optimized rotating speed torque aircraft pursuit course.
Specific implementation mode
The adaptive maximum power tracking and controlling method specific implementation of wind power generating set of the present invention is described below：
The reference tip-speed ratio and reference power coefficient for setting wind power generating set to be measured, further according to reference to tip-speed ratio and reference Power coefficient determines reference rotation velocity torque aircraft pursuit course and its deformation curve.
The reference tip-speed ratio TSR and power coefficient CP of wind power generating set to be measured can be determined by following three kinds of methods：
(1) by numerical computations, wind speed V, rotating speed N and output work when output power maximum under rated wind speed is calculated Rate P, to obtain TSR=N/9.55*R/V and CP=P/ (0.5 × 1.225AV3)；
(2) using the optimal tip-speed ratio and optimal power coefficient of close serial wind-driven generator group rated point.
(3) TSR=3.0, CP=0.1 are directly set.
Rated wind speed V2 refers to the certain wind speed of wind power generating set maximum power output；Incision wind speed V1 refers to wind-force hair When motor group starts power generation, the minimum mean wind speed at hub height；Power coefficient CP refers to the net output of wind power generating set The ratio that electrical power passes through the wind energy in rotor swept face with free stream.Net electromotive power output, abbreviation output power, to off-network type Refer to the electric energy of charge controller output end for wind generator system；Refer to unsteady flow for the Shunt-connected Wind Power Generation System The electric energy that device output end is connected to the grid.
Wind power generating set is needed the wind speed section for optimizing the power of output be set as { V | V1 ＜ V ＜ V2 }, wherein V1 is general To cut wind speed, V2 is generally rated wind speed.R scans radius surface for wind wheel blade.A is the vertical throwing of wind wheel blade swept surface Shadow area；
The rotating speed that wind power generating set corresponds to TSR is set as N；
N=VTSR/R9.55
Corresponding output power is：
P=0.5 × 1.225AV3·CP
Corresponding torque is：
T=P/N9.55
To obtain
T=6.716 × 10-3A·(R/TSR)3·CP·N3/ N=CKN2；
Wherein CK=6.716 × 10-3A·(R/TSR)3·CP
T=CKN2Deformation curve race T=CK (KXPN)2, KXP is known as deformation coefficient, and KXP is a variable, Between value range is 0.5~2.0, when it is less than 1, indicate curve to N positive direction extensional deformations；When it is more than 1, indicate Curve is to the contraction distortion of origin direction.
From incision wind speed to the section of rated wind speed, wind speed interval, each wind speed interval are divided by the intervals 0.5m/s Intermediate value when 0.5m/s integral multiple.Such as incision wind speed is 3m/s, rated wind speed, which is the wind speed interval of 12m/s, to be divided into [2.75,3.25), [3.25,3.75), [3.75,4.25) ..., [11.75,12.25) each section.
When the entire interval division of working speed, at many minizones, then in each minizone, there are one certainly A KXP deformation curves and actual optimum rotational speed and torque aircraft pursuit course are closest.When section is sufficiently small, the KXP deformation curves and Actual optimum curve co-insides.In engineering, when rotating speed section corresponds to the wind speed interval of a 0.5m/s width, wind power generating set KXP deformation curves are tracked with the ability of its output power of actual optimum curve with regard to essentially the same.
Fig. 5 illustrate an actual optimum rotational speed and torque aircraft pursuit course in 3 each adjacent rotating speed sections respectively use KX1, For the deformed curve of KX2, KX3 come the situation being fitted, heavy line therein is exactly the curved section of fitting.
Concrete instance is as follows：
If the rated wind speed of a vertical axis aerogenerator group is 12m/s, the design optimal tip-speed ratio of rated point is 2.5, It is set as referring to tip-speed ratio, rotor diameter is 1.5 meters, and wind wheel blade length is 1.0 meters, and power coefficient is：0.15, rated power：
P=0.5 × 1.225 × 1.5 × 0.15 × 12 × 12 × 12=238W；
Rated speed is：
N=2.5 × 12/0.75 × 9.55=382r/min
Nominal torque is：T=P/N9.55=5.95Nm
Incision wind speed is 3m/s.Cut-out wind speed is assumed to be 20m/s.In 12m/s and 20m/s or more, wind power generating set is pressed According to the method operation of limitation power；In 20m/s or more, wind power generating set cuts out braking.It is specified in incision wind speed 3m/s or more Wind speed 12m/s is hereinafter, wind power generating set need to track the operation of optimized rotating speed torque aircraft pursuit course with Maximum Power Output.Therefore Have：
CK=6.716 × 10-3×1.5×(0.75/2.5)3× 0.15=4.08 × 10-5
It can obtain following initial speed torque aircraft pursuit course T=CKN2；
Incision rotating speed is N1=2.5 × 3/0.75 × 9.55=95.5r/min
Multiple points are chosen between incision rotating speed and rated speed, reference rotation velocity torque aircraft pursuit course data can be obtained, Such as table 1：
The reference power curve that it is indicated is as shown in Figure 6.
By being carried out to reference rotation velocity torque curveNThe scaling of axis direction, deformation curve cluster are T=CK (KXPN )2。
For example, when KXP is { 0.92, 0.9,1,1.1,1.12When each curve it is as shown in Figure 7.
At the beginning because optimized rotating speed torque aircraft pursuit course is unknown, wind power generating set can be by reference rotation velocity torque The KXP deformation curves of aircraft pursuit course track operation as the tracking optimized rotating speed torque aircraft pursuit course of wind power generating set.In order to The power characteristic of wind power generating set is obtained, according to the experiment of GB/T 18451.2-2012 power characteristic of wind driven generator set and GB/ The method that the small-sized vertical axis wind power generation groups of T29494-2013 provide, often runs 10min, just counts the operation knot in the 10min Fruit data.The 10min durations can be set as cycle T 1.These data have：Mean wind speed, mean speed, mean power.At this In, also measurement data is averaged and is counted by every 10min.In addition it needs to increase by two data of storage：Section starts rotating speed area Between terminate rotating speed, the two data be used for calculate the output work caused by wind wheel rotating part kinetic energy change in the sections 10min Rate changing value, referred to as kinetic energy delivered power.
Definition mean wind speed is VSX, mean speed NSX, average output power PSX, kinetic energy delivered power are PDSX.
From start to end for the 10min, since speed changes, wind power generating set turns kinetic energy delivered power PDSX The power of the release of dynamic part kinetic energy；If first group of pretreated rotating speed average value is NS when the 10min starts, at the end of it is last One group of pretreated wind speed average value is NE, then kinetic energy delivered power is：
PDSX=0.5J ((NS/9.55)2-(NE/9.55)2)
Wherein J is total inertia of wind wheel rotating part；
PJSX=PSX-PDSX is defined, PJSX indicates the power that the wind wheel in the 10min is absorbed from wind, referred to as average net defeated Go out power P JSX；
Because the mean wind speed of each 10min, drops into some wind speed interval, but needs not be equal in wind speed interval Value, so the power of 10min twice cannot be compared under the same benchmark.Therefore it needs average net power output PJSX conversions To wind speed interval intermediate value.The reduced value is known as PJSXZ, PJSXZ=PJSX (VS/VSX)3, wherein VS is wind speed interval intermediate value, VS =int (2VSX+0.5)/2.Similarly, corresponding mean speed NSX is also required to conversion and arrives wind speed intermediate value, reduced value NSXZ =NSXVS/VSX.Int () is bracket function.
With the energy of wind power generating set output power when intermediate value wind speed of the PJSXZ to assess the wind speed interval in the 10min Power, and object function as an optimization.
The deformation coefficient used in the 10min is optimized variable.
The optimization problem is the continuous service by wind-driven generator, and optimal optimization is found in each wind speed interval and is become Amount is so that its output power maximizes after the deformation coefficient operation that wind power generating set is determined by the optimized variable.
The optimization problem uses following steps：
Step (1)
The wind speed minizone of 0.5m/s width is divided into from total wind speed interval of incision wind speed V1 rated wind speeds V2, those Wind speed minizone label from low to high.
In this way, IMAX 0.5m/s width can be divided into from incision wind speed V1 to wind speed interval total rated wind speed V2 Wind speed minizone, IMAX=int (2V2+0.5)-int (2V1+0.5)+1, int () are bracket function.
After wind power generating set runs a cycle T1, if the mean wind speed in its T1 period is VSX, then the operation is all Phase, corresponding wind speed minizone label IX was obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
In order to store the data for corresponding to each wind speed minizone, the array that following width is IMAX is opened up；
KX () --- storage corresponds to the deformation coefficient of each wind speed minizone, initial value 1.0；
VSX () --- storage corresponds to the mean wind speed value of each wind speed minizone, initial value 0；
PJSXZ () --- storage corresponds to the net power output value after the conversion of each wind speed minizone, initial value 0；
NSXZ () --- storage corresponds to the reduced value of the mean speed value of each wind speed minizone, initial value 0；
NN () --- storage corresponds to the optimization number of each wind speed minizone, initial value 0；
FB () --- storage corresponds to the optimization perturbation direction mark of each wind speed minizone, initial value 1；
Initial deformation COEFFICIENT K XP=1.0 is set, and goes to next step.
Step (2)
It is according to the deformation curve for being worth to reference rotation velocity torque aircraft pursuit course of KXP：
T=CK (KXPN)2
Wherein coefficient
CK=6.716 × 10-3·A·(R/TSR)3·CP
R scans radius surface for wind wheel blade；
Projected areas of the A for wind wheel blade swept surface on carrying out wind direction vertical plane；
TSR is with reference to tip-speed ratio, and value range is：2.0~3.0
CP is reference power coefficient, and value range is：0.1~0.25
When wind power generating set is run, the wind speed round N of wind power generating set is measured in real time, then passes through deformation curve letter Number obtains torque instruction value T, then controls wind power generating set reality output torque TL and tracks torque instruction value T operations.
Go to next step.
Step (3)
Wind power generating set by 2 periodic measurement of cycle T and obtains the T2 by step (2) method control operation T1 periods The real-time air speed value in period, real-time tachometer value and real-time output power value, T1=mT2, m are no less than 50 positive integer, T1> >T2；Wind-driven generator runs the T1 periods, obtains the real-time air speed value in m T2 period, real-time tachometer value and real-time output power value Average value VSX, NSX and PSX, and obtain the real-time tachometer value NS in first T2 period and the T1 periods in the T1 periods simultaneously The real-time tachometer value NE in the last one interior T2 period.
Then, it is calculate by the following formula to obtain the T1Kinetic energy delivered power PDSX in period：
PDSX=0.5J ((NS/9.55)2-(NE/9.55)2)
Wherein J is the rotary inertia of wind wheel rotating part.
It is calculate by the following formula to obtain the T again1The average net power output PJSX that wind wheel in period is absorbed from wind：
PJSX=PSX-PDSX
It is calculate by the following formula to obtain the T1The wind speed intermediate value in the affiliated 0.5m/s wind speed width sections period mean wind speed VSX VS：
VS=int (2VSX+0.5)/2
Int () is bracket function.
PJSX is converted by following formula by average net power output PJSXZP corresponding to wind speed intermediate value VS again：
PJSXZP=PJSX (VS/VSX)3
NSX is converted average net power output NSXZP corresponding to wind speed intermediate value VS by following formula：
NSXZP=NSX (VS/VSX)
For storing the array label IX for representing different wind speed intervals, can be obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
Go to next step.
Step (4)
Step (5)
Step (6)
Step (7)
If the total continuous operating time of wind-driven generator is more than certain setting value, i.e., by long lasting for continuing to optimize, The deformation coefficient and operation mean speed of each wind speed minizone of optimization are stored in KX (IX) and NSXZ (IX), under TX (IX) is calculated in formula：
TX (IX)=CK (KX (IX) NSXZ (IX))2, IX=1,2 ..., IMAX
So as to obtain one group of data point (NSXZ (IX), TX (IX))
The broken line that the data point connects is optimized rotating speed torque aircraft pursuit course.
Go to step (2) circular flow.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to restrict the invention, it is clear that those skilled in the art Various changes and modifications can be made to the invention by member without departing from the spirit and scope of the present invention.If in this way, the present invention Within the scope of the claims of the present invention and its equivalent technology, then the present invention is also intended to include these these modifications and variations Including modification and variation.

Claims (1)

1. a kind of adaptive maximum power tracking and controlling method of wind power generating set, which is characterized in that including steps are as follows：
Step (1)
It is divided into the wind speed minizone of 0.5m/s width, the wind from total wind speed interval of incision wind speed V1, rated wind speed V2 Fast minizone label from low to high；
In this way, the wind speed of IMAX 0.5m/s width can be divided into from incision wind speed V1 to wind speed interval total rated wind speed V2 Minizone, IMAX=int (2V2+0.5)-int (2V1+0.5)+1, int () are bracket function；
After wind power generating set runs a cycle T1, if mean wind speed in its T1 period is VSX, then cycle of operation pair The wind speed minizone label IX answered is obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
In order to store the data for corresponding to each wind speed minizone, the array that following width is IMAX is opened up；
KX () --- storage corresponds to the deformation coefficient of each wind speed minizone, initial value 1.0；
VSX () --- storage corresponds to the mean wind speed value of each wind speed minizone, initial value 0；
PJSXZ () --- storage corresponds to the net power output value after the conversion of each wind speed minizone, initial value 0；
NSXZ () --- storage corresponds to the reduced value of the mean speed value of each wind speed minizone, initial value 0；
NN () --- storage corresponds to the optimization number of each wind speed minizone, initial value 0；
FB () --- storage corresponds to the optimization perturbation direction mark of each wind speed minizone, initial value 1；
Initial deformation COEFFICIENT K XP=1.0 is set, and goes to next step；
Step (2)
It is according to the deformation curve for being worth to reference rotation velocity torque aircraft pursuit course of KXP：
T=CK (KXPN)2
Wherein coefficient
CK=6.716 × 10-3·A·(R/TSR)3·CP；
R scans radius surface for wind wheel blade；Projected areas of the A for wind wheel blade swept surface on carrying out wind direction vertical plane；
TSR is with reference to tip-speed ratio, and value range is：2.0~3.0；CP is reference power coefficient, and value range is：0.1~ 0.25；
When wind power generating set is run, the wind speed round N of wind power generating set is measured in real time, is then obtained by above-mentioned deformation curve To torque instruction value T, then controls wind power generating set reality output torque TL and track torque instruction value T operations, go to down One step；
Step (3)
Wind power generating set by 2 periodic measurement of cycle T and obtains the T2 periods by step (2) method control operation T1 periods Real-time air speed value, real-time tachometer value and real-time output power value, T1=mT2, m are larger positive integer, T1>>T2；Wind-force Generator runs the T1 periods, obtains the real-time air speed value in m T2 period, the average value of real-time tachometer value and real-time output power value VSX, NSX and PSX, and last is obtained in the T1 periods in the real-time tachometer value NS in first T2 period and the T1 periods simultaneously The real-time tachometer value NE in a T2 periods；
Then, it is calculate by the following formula to obtain the T1Kinetic energy delivered power PDSX in period：
PDSX=0.5J ((NS/9.55)2-(NE/9.55)2)
Wherein J is the rotary inertia of wind wheel rotating part；
It is calculate by the following formula to obtain the T again1The average net power output PJSX that wind wheel in period is absorbed from wind：
PJSX=PSX-PDSX
It is calculate by the following formula to obtain the T1The wind speed intermediate value VS in the affiliated 0.5m/s wind speed width sections period mean wind speed VSX：
VS=int (2VSX+0.5)/2
Int () is bracket function；
PJSX is converted by following formula by average net power output PJSXZP corresponding to wind speed intermediate value VS again：
PJSXZP=PJSX (VS/VSX)3
NSX is converted average net power output NSXZP corresponding to wind speed intermediate value VS by following formula：
NSXZP=NSX (VS/VSX)
For storing the array label IX for representing different wind speed intervals, can be obtained by following formula：
IX=int (2VSX+0.5)-int (2V1+0.5)+1
Go to next step；
Step (4)
Go to next step；
Step (5)
Go to next step；
Step (6)
Go to next step；
Step (7)
If the total continuous operating time of wind-driven generator is more than certain setting value, i.e., by long lasting for continuing to optimize, KX (IX) and in NSXZ (IX) it is stored with the deformation coefficient and operation mean speed of each wind speed minizone of optimization, according to the following formula TX (IX) is calculated：
TX (IX)=CK (KX (IX) NSXZ (IX))2, IX=1,2 ..., IMAX
So as to obtain one group of data point (NSXZ (IX), TX (IX))
The broken line that the data point connects is optimized rotating speed torque aircraft pursuit course；
Go to step (2) circular flow.
CN201610235513.XA 2016-04-15 2016-04-15 The adaptive maximum power tracking and controlling method of wind power generating set Active CN105909470B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610235513.XA CN105909470B (en) 2016-04-15 2016-04-15 The adaptive maximum power tracking and controlling method of wind power generating set

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610235513.XA CN105909470B (en) 2016-04-15 2016-04-15 The adaptive maximum power tracking and controlling method of wind power generating set

Publications (2)

Publication Number Publication Date
CN105909470A CN105909470A (en) 2016-08-31
CN105909470B true CN105909470B (en) 2018-08-24

Family

ID=56747096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610235513.XA Active CN105909470B (en) 2016-04-15 2016-04-15 The adaptive maximum power tracking and controlling method of wind power generating set

Country Status (1)

Country Link
CN (1) CN105909470B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107100795B (en) * 2017-07-05 2019-04-09 东方电气自动控制工程有限公司 A kind of low wind speed apparatus for lower wind generating set MPPT self-adaptation control method
CN109931230B (en) * 2017-12-19 2020-02-28 北京金风科创风电设备有限公司 Method and device for detecting active power of wind generating set
CN110067708B (en) * 2019-05-13 2020-08-25 北京天泽智云科技有限公司 Method for identifying yaw wind disharmony by using power curve
CN113007012B (en) * 2019-12-19 2022-09-23 新疆金风科技股份有限公司 Torque control coefficient optimizing method and device and wind generating set
CN113090453B (en) * 2019-12-23 2023-03-03 新疆金风科技股份有限公司 Control method and device of wind generating set and wind generating set
CN111336062B (en) * 2020-03-05 2021-11-09 中国大唐集团科学技术研究院有限公司华中电力试验研究院 Wind generating set maximum wind energy capture method based on measured wind speed
CN114320741A (en) * 2020-09-30 2022-04-12 新疆金风科技股份有限公司 Power control method and equipment of wind generating set
CN114673630A (en) * 2020-12-24 2022-06-28 新疆金风科技股份有限公司 Method and device for determining speed ratio of blade tip of wind turbine generator and main controller
CN116412074B (en) * 2023-03-07 2024-02-23 浙江大学 Wind driven generator torque control method and device

Also Published As

Publication number Publication date
CN105909470A (en) 2016-08-31

Similar Documents

Publication Publication Date Title
CN105909470B (en) The adaptive maximum power tracking and controlling method of wind power generating set
CN101943116B (en) Method and system for noise controlled operation of a wind turbine
EP2556249B1 (en) A wind turbine
CA2764612C (en) Methods and systems for determining a pitch angle offset signal and for controlling a rotor frequency of a rotor of a wind turbine for speed avoidance control
EP1906192B1 (en) Apparatus for evaluating sensors and/or for controlling operation of an apparatus that includes a sensor
US8022566B2 (en) Methods and systems for operating a wind turbine
EP2306003B1 (en) System and methods for controlling a wind turbine
US7352075B2 (en) Methods and apparatus for controlling rotational speed of a rotor
EP2556248B1 (en) A wind turbine
CA2810157C (en) Method of rotor-stall prevention in wind turbines
EP2273109A1 (en) Wind turbine acoustic emission control system and method
EP2757252A1 (en) Method of operating a wind turbine
CA2681784C (en) A speed control for wind turbines
CN102156044B (en) Model selection method of wind turbine simulator applicable to testing of direct driving type wind generating set
CN103375332B (en) Speed-changing oar-changing wind power generating set optimum moment of resistance dynamic optimization method
CN109882355A (en) The dynamic matching method and system of low wind speed apparatus for lower wind generating efficiency in a kind of promotion
CN109779837B (en) Yaw alignment correction method for wind generating set
CN116258393A (en) Wind power generation safety checking method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201821 Shanghai Jiading District City Yongsheng Road No. 2201, No. 655, building 5, Zhaoxian Road area, two layer 1 layer A

Applicant after: Shanghai ruibode intelligent system Limited by Share Ltd

Address before: 200335 Shanghai city Changning District road of 8 Building No. 3

Applicant before: Shanghai Ruibode Intelligent System Sci. & Tech. Co., Ltd.

GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 201821 Shanghai Jiading District City Yongsheng Road No. 2201, No. 655, building 5, Zhaoxian Road area, two layer 1 layer A

Patentee after: SHANGHAI CHINESE CAR RIBERD INTELLIGENT SYSTEM Co.,Ltd.

Address before: 201821 Shanghai Jiading District City Yongsheng Road No. 2201, No. 655, building 5, Zhaoxian Road area, two layer 1 layer A

Patentee before: SHANGHAI RUIBODE INTELLIGENT SYSTEMS Co.,Ltd.