CN105373858B - A wind farm active power optimization method based on wind speed time series decomposition - Google Patents
A wind farm active power optimization method based on wind speed time series decomposition Download PDFInfo
- Publication number
- CN105373858B CN105373858B CN201510843636.7A CN201510843636A CN105373858B CN 105373858 B CN105373858 B CN 105373858B CN 201510843636 A CN201510843636 A CN 201510843636A CN 105373858 B CN105373858 B CN 105373858B
- Authority
- CN
- China
- Prior art keywords
- unit
- wind
- active power
- wind speed
- value
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000000354 decomposition reaction Methods 0.000 title 1
- 230000000694 effects Effects 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 230000001939 inductive effect Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000005684 electric field Effects 0.000 description 6
- 230000005611 electricity Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a kind of active power of wind power field optimization methods decomposed based on wind speed timing, comprising the following steps: step 1: the forecasting wind speed Curve transform by wind power plant next control period is discrete wind series;Step 2: being considered as constant natural wind speed for the corresponding air speed value of discrete wind series each element, calculates separately under these constant natural wind speed effects, the active power reference value and its initial time of each unit;Step 3: by each unit active power reference value and its initial time, each unit active power reference value controlling curve under wind series effect is generated, each unit is run according to this controlling curve.The lower calculation amount of this optimization method, less hardware resource realize that active power of wind power field optimizes under Dynamic Wind Speed, and forming one kind has practical value scheme, promotes wind power plant generated energy, improves wind power plant gene-ration revenue.
Description
Technical field
The present invention relates to wind power plant control technology field, in particular to a kind of wind power plant decomposed based on wind speed timing is active
Power optimization method.
Background technique
In recent years, the Wind Power Generation Industry in China remains strong growth, with the development of wind-powered electricity generation, there is good wind energy
Land resources are fewer and fewer, further develop for support Wind Power Generation Industry, extensive development marine wind electric field becomes inexorable trend.When
Preceding marine wind electric field construction cost is higher, improves the generating efficiency of marine wind electric field, increases the economic benefit of wind-powered electricity generation owner, is sea
One of the critical issue of upper wind power plant operation.
Compared with landwid electric field, marine wind electric field is extremely violent by wake effect, and traditional separate unit Wind turbines are maximum
Wind energy extraction will lead to upwind unit and absorb biggish wind energy, and lower wind direction unit input wind speed reduces, to lose wind power plant
Active power, studies have shown that traditional single machine maximal wind-energy capture scheme, the marine wind electric field as caused by wake effect are active
Power loss may be up to 30%, although wake effect can be effectively reduced by wind power plant microcosmic structure, Wind turbines layout optimization
Influence, but due to the factors such as wind power plant place and construction cost restrict, the distance between Wind turbines are generally in impeller diameter
5 to 7 times between, influence of the wake effect to active power of wind power field is still obvious.
Domestic and international some scholars are based on wake effect and have carried out a series of researchs to active power of wind power field optimization, form
A variety of wake models, existing research concentrate on substantially wind speed it is constant when, influence of the wake effect to active power of wind power field and its
In terms of ameliorative way.But in wind speed dynamic change, the research of the active power optimization aspect of wind power plant is still immature, main
The reason is that Wind turbines have it is very strong non-linear, by traditional method of thinking, using the wind speed of longer period of time as one
Entirety is investigated, and active power of wind power field optimization method is established, and since the possibility combination of time dimension is extremely more, makes equation solution
Computationally intensive, corresponding control strategy realizes difficulty, is unfavorable for practical application.
Summary of the invention
That in order to solve the above technical problem, the present invention provides a kind of calculation amounts is small, at low cost, practical value is high based on wind
The active power of wind power field optimization method that fast timing is decomposed.
Technical proposal that the invention solves the above-mentioned problems is: a kind of active power of wind power field decomposed based on wind speed timing is excellent
Change method, comprising the following steps:
Step 1: the forecasting wind speed curve in wind power plant next control period is chronologically converted into the discrete of 30 elements
Wind series;
Step 2: being considered as constant natural wind speed for the corresponding air speed value of discrete wind series each element, is based on wake effect,
It calculates separately under each constant natural wind speed effect, the active power reference value of each unit;Meanwhile it being based on wake flow propagation delay,
It calculates separately under each constant natural wind speed effect, at the time of each unit active power reference value starts effect;
Step 3: by each unit active power reference value and its act at the time of, generate the effect of discrete wind series
Under each unit active power reference value controlling curve, each unit according to this controlling curve run.
The above-mentioned active power of wind power field optimization method decomposed based on wind speed timing, in the step 1, discrete wind speed sequence
Arrange the conversion formula of s (m) are as follows:
①
Wherein, m is natural number, and value range is 1≤m≤31, scIt (t) is t moment predicting wind speed of wind farm value, tNTable
Show that the initial time in next control period, T indicate one control period corresponding duration of wind power plant, TsIt is between wind speed discretization
Every the time, Ts=T/30.
The above-mentioned active power of wind power field optimization method decomposed based on wind speed timing, in the step 2, each unit has
Steps are as follows for the calculating of function value and power reference:
1) with the axial inducible factor a of unit iiFor variable, then the power coefficient of unit iAnd its thrust coefficientAre as follows:
②
③
2) m in discrete wind series s (m) takes occurrence n, enables n=1;
3) it sets wind power plant to be acted on by constant natural wind speed s (n), uses viThe input wind speed for indicating unit i, if v1It is upwind
Boundary unit, unit i+1 are the upward First units of unit i leeward, 4. 5. calculate the input wind speed of unit i+1 with formula by formula:
v1=s (n) is 4.
⑤
Wherein, k is wind power plant terrain rough factor, riFor unit i impeller radius, x is unit i hub centre and unit
The distance that i+1 hub centre projects on wind direction;
4) each unit active power is calculated as follows:
⑥
Wherein, PiIndicate the active power of unit i, vratedIndicate unit rated wind speed, vin、vcutThe incision of expression unit,
Cut-out wind speed;ρ indicates atmospheric density,Indicate the rated power of unit i;
5) it is as follows to establish active power of wind power field Optimized model:
⑦
Wherein, PallIndicate that active power of wind power field, i.e. the sum of each unit active power of wind power plant, N indicate machine in wind power plant
Group sum,It indicates through optimization aiThe formula in bracket is set to obtain maximum value;
6) solution procedure 5) given by equation, obtain the optimization solution of unit i inducible factor, and substituted into formula 2.-
5. obtaining unit i power coefficient, thrust coefficient, the specific value for inputting wind speed;
7) by unit i power coefficient value substitute into following formula, then obtain natural wind speed be s (n) when, the wattful power of unit i
Rate reference value
⑧
Make i value traverse each unit of wind power plant, then obtain natural wind speed be s (n) when, the active power reference of each unit
Value;
8) judge that the value size of n enables n=n+1 if n≤30, jump to step 3), otherwise terminate operation.
The above-mentioned active power of wind power field optimization method decomposed based on wind speed timing, in the step 2, each unit has
Calculation method at the time of function value and power reference starts effect is as follows: when the m in discrete wind series s (m) takes occurrence n, counting
Calculate the time delay D that wind speed s (n) travels to unit i from wind power plant upwind boundary uniti,n, then when wind speed is s (n), unit i's has
Function value and power referenceThe initiation momentFor
⑨。
The above-mentioned active power of wind power field optimization method decomposed based on wind speed timing, the active power reference value of each unit
Start at the time of effect in calculation method, wind speed s (n) travels to the time delay D of unit i from wind power plant upwind boundary uniti,n's
Calculation method is as follows:
1) wind is propagated by wind direction, on wind speed propagation path, seeks wind power plant upwind boundary unit impeller center and machine
The distance of subpoint of the group i impeller center on wind direction, is denoted as L;
2) D is calculated as followsi,n, and i value is made to traverse each unit of wind power plant, n value traverses each nature in 1≤n≤31
Number,
⑩。
The above-mentioned active power of wind power field optimization method decomposed based on wind speed timing, in the step 3, unit wattful power
The generating mode of rate reference value controlling curve is as follows: by i value, n value difference is identicalWithComposition control amount
It is meant thatThe active power reference value at moment, unit i isIt, will to determining i valueValue by n by it is small to
Longer spread,It arrivesBetween time, the active power reference value of unit i maintainsIt is constant, to obtain unit i
Active power reference value controlling curve;I value is set to traverse each unit of wind power plant, m value traverses each natural number in 1≤m≤31,
To obtain each unit active power reference value controlling curve of wind power plant.
The beneficial effects of the present invention are: the prediction of wind speed Curve transform in wind power plant next control period is by the present invention
Discrete wind series calculate constant nature after the corresponding air speed value of discrete wind series each element is considered as constant natural wind speed
Wind speed effect under each unit active power reference value and its start effect at the time of, and then generate wind series effect under each unit
Active power reference value controlling curve, each unit are run according to this controlling curve;It is the lower calculation amount of this optimization method, less
Hardware resource, realize that active power of wind power field optimizes under Dynamic Wind Speed, formed it is a kind of there is practical value scheme, promote wind-powered electricity generation
Field generated energy, improves wind power plant gene-ration revenue.
Detailed description of the invention
Fig. 1 is general flow chart of the invention.
Fig. 2 is the flow chart of step 3 of the present invention.
Specific embodiment
The present invention is further illustrated with reference to the accompanying drawings and examples.
As shown in Figure 1, 2, the present invention the following steps are included:
Step 1: the forecasting wind speed curve in wind power plant next control period is chronologically converted into the discrete of 30 elements
Wind series.
The conversion formula of discrete wind series s (m) are as follows:
①
Wherein, m is natural number, and value range is 1≤m≤31, scIt (t) is t moment predicting wind speed of wind farm value, tNTable
Show that the initial time in next control period, T indicate one control period corresponding duration of wind power plant, TsIt is between wind speed discretization
Every the time, Ts=T/30.
Step 2: being considered as constant natural wind speed for the corresponding air speed value of discrete wind series each element, is based on wake effect,
It calculates separately under each constant natural wind speed effect, the active power reference value of each unit;Meanwhile it being based on wake flow propagation delay,
It calculates separately under each constant natural wind speed effect, at the time of each unit active power reference value starts effect.
Steps are as follows for the calculating of the active power reference value of each unit:
1) with the axial inducible factor a of unit iiFor variable, then the power coefficient of unit iAnd its thrust coefficientAre as follows:
②
③
2) m in discrete wind series s (m) takes occurrence n, enables n=1;
3) it sets wind power plant to be acted on by constant natural wind speed s (n), uses viThe input wind speed for indicating unit i, if v1It is upwind
Boundary unit, unit i+1 are the upward First units of unit i leeward, 4. 5. calculate the input wind speed of unit i+1 with formula by formula:
v1=s (n) is 4.
⑤
Wherein, k is wind power plant terrain rough factor, riFor unit i impeller radius, x is unit i hub centre and unit
The distance that i+1 hub centre projects on wind direction;
4) each unit active power is calculated as follows:
⑥
Wherein, PiIndicate the active power of unit i, vratedIndicate unit rated wind speed, vin、vcutThe incision of expression unit,
Cut-out wind speed;ρ indicates atmospheric density,Indicate the rated power of unit i;
5) it is as follows to establish active power of wind power field Optimized model:
⑦
Wherein, PallIndicate that active power of wind power field, i.e. the sum of each unit active power of wind power plant, N indicate machine in wind power plant
Group sum,It indicates through optimization aiValue makes the formula in bracket obtain maximum value;
6) solution procedure 5) given by equation, obtain the optimization solution of unit i inducible factor, and substituted into formula 2.-
5. obtaining unit i power coefficient, thrust coefficient, the specific value for inputting wind speed;
7) by unit i power coefficient value substitute into following formula, then obtain natural wind speed be s (n) when, the wattful power of unit i
Rate reference value
⑧
Make i value traverse each unit of wind power plant, then obtain natural wind speed be s (n) when, the active power reference of each unit
Value;
8) judge that the value size of n enables n=n+1 if n≤30, jump to step 3), otherwise terminate the partial arithmetic.
It is as follows that the active power reference value of each unit starts calculation method at the time of effect:
When m in discrete wind series s (m) takes occurrence n, calculation of wind speed s (n) is passed from wind power plant upwind boundary unit
It is multicast to the time delay D of unit ii,n, then when wind speed is s (n), the active power reference value of unit iThe initiation momentFor
⑨
Wind speed s (n) travels to the time delay D of unit i from wind power plant upwind boundary uniti,nCalculation method it is as follows:
1) wind is propagated by wind direction, on wind speed propagation path, seeks wind power plant upwind boundary unit impeller center and machine
The distance of subpoint of the impeller center of group i on wind direction, is denoted as L;
2) D is calculated as followsi,n, and i value is made to traverse each unit of wind power plant, n value traverses each nature in 1≤n≤31
Number,
⑩
Step 3: by each unit active power reference value and its act on the moment, generate under the effect of discrete wind series
Each unit active power reference value controlling curve, each unit are run according to this controlling curve.
The generation method of unit active power reference value controlling curve is as follows: by i value, n value difference is identicalWithGroup
At control amountIt is meant thatThe active power reference value at moment, unit i isIt, will to determining i valueValue press the ascending arrangement of n,It arrivesBetween time, the active power reference value of unit i maintainsIt is constant, so that i value is traversed each unit of wind power plant, m value traverses each natural number in 1≤m≤31, to obtain wind-powered electricity generation
Each unit active power reference value controlling curve in field.It can specifically carry out as follows:
1) if wind power plant unit sum is N, the bivariate table of three N rows 30 column is constructed, TableT is respectively designated as,
TableP and TalblePandT;
2) the possibility value of i and n is traversed, it willFill out the i-th row n-th column in TableT;
3) the possibility value of i and n is traversed, it willFill out the i-th row n-th column in TableP;
4) i-th the n-th column element of row of TableT and i-th the n-th column element of row of TableP are complex asAnd by its
Write on the i-th row n-th column of TalblePandT;
5) using the time as abscissa, active power reference value is that ordinate constitutes coordinate system, takes the i-th of TalblePandT
Row data, by row number, sequence is successively fetched evidence from small to large, and is demarcated on coordinate system each Point, also,
It arrivesBetween time, the active power reference value of unit i is maintainedIt is constant, to obtain the wattful power of unit i
Rate reference value controlling curve;All possible values of i traversal, repeat the step, obtain the next control period active power of each unit
Reference value controlling curve.
Claims (6)
1. a kind of active power of wind power field optimization method decomposed based on wind speed timing, comprising the following steps:
Step 1: the forecasting wind speed curve in wind power plant next control period is chronologically converted to the discrete wind speed of 30 elements
Sequence;
Step 2: being considered as constant natural wind speed for the corresponding air speed value of discrete wind series each element, is based on wake effect, respectively
It calculates under each constant natural wind speed effect, the active power reference value of each unit;Meanwhile it being based on wake flow propagation delay, respectively
It calculates under each constant natural wind speed effect, at the time of each unit active power reference value starts effect;
Step 3: by each unit active power reference value and its act at the time of, generate each under the effect of discrete wind series
Unit active power reference value controlling curve, each unit are run according to this controlling curve.
2. the active power of wind power field optimization method according to claim 1 decomposed based on wind speed timing, it is characterised in that:
In the step 1, the conversion formula of discrete wind series s (m) are as follows:
①
Wherein, m is natural number, and value range is 1≤m≤31, scIt (t) is t moment predicting wind speed of wind farm value, tNUnder expression
The initial time in one control period, T indicate one control period corresponding duration of wind power plant, TsWhen being wind speed discretization interval
Between, Ts=T/30.
3. the active power of wind power field optimization method according to claim 2 decomposed based on wind speed timing, it is characterised in that:
In the step 2, steps are as follows for the calculating of the active power reference value of each unit:
1) with the axial inducible factor a of unit iiFor variable, then the power coefficient of unit iAnd its thrust coefficientAre as follows:
②
③
2) m in discrete wind series s (m) takes occurrence n, enables n=1;
3) it sets wind power plant to be acted on by constant natural wind speed s (n), uses viThe input wind speed for indicating unit i, if v1It is upwind boundary
Unit, unit i+1 are the upward First units of unit i leeward, 4. 5. calculate the input wind speed of unit i+1 with formula by formula:
v1=s (n) is 4.
⑤
Wherein, k is wind power plant terrain rough factor, riFor unit i impeller radius, x is that unit i hub centre and unit i+1 take turns
The distance that hub center projects on wind direction;
4) each unit active power is calculated as follows:
⑥
Wherein, PiIndicate the active power of unit i, vratedIndicate unit rated wind speed, vin、vcutIt indicates the incision of unit, cut out
Wind speed;ρ indicates atmospheric density,Indicate the rated power of unit i;
5) it is as follows to establish active power of wind power field Optimized model:
s.t.
v1=s (n)
⑦
Wherein, PallIndicate that active power of wind power field, i.e. the sum of each unit active power of wind power plant, N indicate that unit is total in wind power plant
Number,It indicates through optimization aiThe formula in bracket is set to obtain maximum value;
6) solution procedure 5) given by equation, obtain the optimization solution of unit i inducible factor, and substituted into formula 2. -5., obtain
Unit i power coefficient, thrust coefficient, the specific value for inputting wind speed out;
7) by unit i power coefficient value substitute into following formula, then obtain natural wind speed be s (n) when, unit i active power ginseng
Examine value
⑧
Make i value traverse each unit of wind power plant, then obtain natural wind speed be s (n) when, the active power reference value of each unit;
8) judge that the value size of n enables n=n+1 if n≤30, jump to step 3), otherwise terminate operation.
4. the active power of wind power field optimization method according to claim 3 decomposed based on wind speed timing, it is characterised in that:
In the step 2, the calculation method that the active power reference value of each unit starts at the time of acting on is as follows: discrete wind series s
(m) when the m in takes occurrence n, calculation of wind speed s (n) travels to the time delay D of unit i from wind power plant upwind boundary uniti,n, then
When wind speed is s (n), the active power reference value of unit iThe initiation momentFor
⑨。
5. the active power of wind power field optimization method according to claim 4 decomposed based on wind speed timing, it is characterised in that:
The active power reference value of each unit starts at the time of effect in calculation method, and wind speed s (n) is from wind power plant upwind boundary machine
Group travels to the time delay D of unit ii,nCalculation method it is as follows:
1) wind is propagated by wind direction, on wind speed propagation path, seeks wind power plant upwind boundary unit impeller center and unit i leaf
The distance of subpoint of the wheel center on wind direction, is denoted as L;
2) D is calculated as followsi,n, and i value is made to traverse each unit of wind power plant, n value traverses each natural number in 1≤n≤31,
⑩。
6. the active power of wind power field optimization method according to claim 5 decomposed based on wind speed timing, it is characterised in that:
In the step 3, the generating mode of unit active power reference value controlling curve is as follows: by i value, n value difference is identical
WithComposition control amountIt is meant thatThe active power reference value at moment, unit i isTo determining i
Value, willValue press the ascending arrangement of n,It arrivesBetween time, the active power reference value of unit i maintainsIt is constant, to obtain the active power reference value controlling curve of unit i;I value is set to traverse each unit of wind power plant, m takes
Each natural number in value traversal 1≤m≤31, to obtain each unit active power reference value controlling curve of wind power plant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510843636.7A CN105373858B (en) | 2015-11-26 | 2015-11-26 | A wind farm active power optimization method based on wind speed time series decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510843636.7A CN105373858B (en) | 2015-11-26 | 2015-11-26 | A wind farm active power optimization method based on wind speed time series decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105373858A CN105373858A (en) | 2016-03-02 |
CN105373858B true CN105373858B (en) | 2019-05-07 |
Family
ID=55376037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510843636.7A Active CN105373858B (en) | 2015-11-26 | 2015-11-26 | A wind farm active power optimization method based on wind speed time series decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105373858B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105868848B (en) * | 2016-03-25 | 2019-05-21 | 湘潭大学 | The wind power plant synergisting method of yaw and active power comprehensive coordination between a kind of unit |
CN106203695B (en) * | 2016-07-07 | 2020-01-14 | 华北电力大学 | Optimal scheduling method for reducing wake effect in wind power plant |
CN109409596B (en) * | 2018-10-22 | 2021-04-13 | 东软集团股份有限公司 | Processing method, device, equipment and computer-readable storage medium for predicting wind speed |
CN109274121B (en) * | 2018-11-15 | 2021-03-23 | 山东中车风电有限公司 | Wind power plant control parameter optimization method and system |
CN110321632B (en) * | 2019-07-02 | 2021-03-16 | 华北电力大学 | A method for calculating the equivalent roughness of fully developed wind farms |
CN111259552A (en) * | 2020-01-19 | 2020-06-09 | 中国电力科学研究院有限公司 | Large-scale wind power plant probability type time sequence output simulation method and device |
CN113708386B (en) * | 2021-08-26 | 2023-08-18 | 广东电网有限责任公司 | Method, device, equipment and medium for determining stability domain of offshore wind power network |
CN116123028A (en) * | 2022-12-20 | 2023-05-16 | 中国华能集团清洁能源技术研究院有限公司 | Wind farm-level MPPT predictive model control method and device |
CN118174362B (en) * | 2024-05-16 | 2024-07-30 | 国网上海市电力公司 | An island operation scheduling method considering the differences of distributed generation in different fault areas |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663251A (en) * | 2012-04-09 | 2012-09-12 | 华北电力大学 | Physical prediction method for wind power station power based on computational fluid mechanics model |
CN102767474A (en) * | 2012-08-06 | 2012-11-07 | 广东电网公司电力科学研究院 | Real-time assessment method and system for active power for wind power |
CN103745024A (en) * | 2013-11-26 | 2014-04-23 | 沈阳工业大学 | Power characteristic evaluation method for correcting tail wind speed of wind turbine generator based on three-dimensional wake flow model |
CN103942608A (en) * | 2014-03-28 | 2014-07-23 | 华北电力大学 | Optimized dispatching method for wind power farm based wake flow models |
CN104201686A (en) * | 2014-07-24 | 2014-12-10 | 清华大学 | Wind power plant reactive power control method by taking mountain earth wind resource characteristics into account |
-
2015
- 2015-11-26 CN CN201510843636.7A patent/CN105373858B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663251A (en) * | 2012-04-09 | 2012-09-12 | 华北电力大学 | Physical prediction method for wind power station power based on computational fluid mechanics model |
CN102767474A (en) * | 2012-08-06 | 2012-11-07 | 广东电网公司电力科学研究院 | Real-time assessment method and system for active power for wind power |
CN103745024A (en) * | 2013-11-26 | 2014-04-23 | 沈阳工业大学 | Power characteristic evaluation method for correcting tail wind speed of wind turbine generator based on three-dimensional wake flow model |
CN103942608A (en) * | 2014-03-28 | 2014-07-23 | 华北电力大学 | Optimized dispatching method for wind power farm based wake flow models |
CN104201686A (en) * | 2014-07-24 | 2014-12-10 | 清华大学 | Wind power plant reactive power control method by taking mountain earth wind resource characteristics into account |
Also Published As
Publication number | Publication date |
---|---|
CN105373858A (en) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105373858B (en) | A wind farm active power optimization method based on wind speed time series decomposition | |
Sun et al. | Investigation into spacing restriction and layout optimization of wind farm with multiple types of wind turbines | |
Barnes et al. | Determining the impact of VAWT farm configurations on power output | |
CN106979126B (en) | Wind power generating set high wind speed section effective wind speed estimation method based on SVR | |
Lee et al. | Performance study on a counter-rotating tidal current turbine by CFD and model experimentation | |
Mohamed et al. | Optimization of blade pitch angle of an axial turbine used for wave energy conversion | |
CN107240933B (en) | A wind-fire coordination rolling scheduling method considering wind power characteristics | |
Shaaban et al. | Optimization of H-Rotor Darrieus turbines' mutual interaction in staggered arrangements | |
WO2023087521A1 (en) | Wind power plant layout optimization method based on mathematical programming | |
CN106374465B (en) | Short-term wind power forecast method based on GSA-LSSVM model | |
CN101949363A (en) | Method for grouping wind generating sets by taking input wind speed and random fluctuation of wind direction of wind farm into consideration | |
CN105868848B (en) | The wind power plant synergisting method of yaw and active power comprehensive coordination between a kind of unit | |
CN103617308B (en) | A kind of construction method of wind power plant frequency domain equivalent model | |
CN107035618A (en) | A kind of wind power plant optimal control method based on wake effect | |
CN102254092B (en) | A dynamic equivalent method for large wind farms with doubly-fed wind turbines | |
Lanzilao et al. | Effects of self-induced gravity waves on finite wind-farm operations using a large-eddy simulation framework | |
Duan et al. | Wake topology and energy recovery in floating horizontal-axis wind turbines with harmonic surge motion | |
CN106529706A (en) | Support-vector-machine-regression-based method for predicting wind speed of wind power plant | |
CN103760779A (en) | Direct-driven wind generating set dynamic simulation system | |
CN115544884A (en) | A data-driven fast calculation method and system for large-scale wind farm wake | |
CN104675629A (en) | Maximum wind energy capturing method of variable-speed wind generating sets | |
CN102738792A (en) | Wind power predicting method | |
CN110571850A (en) | A method for predicting and correcting control of wind farm power fluctuation trajectory | |
CN102609590B (en) | Wind power farm group layout method | |
Hachmann et al. | Experimentally validated study of the impact of operating strategies on power efficiency of a turbine array in a bi-directional tidal channel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |