CN113947037A - A kind of wind turbine load calculation method and device - Google Patents
A kind of wind turbine load calculation method and device Download PDFInfo
- Publication number
- CN113947037A CN113947037A CN202111028997.8A CN202111028997A CN113947037A CN 113947037 A CN113947037 A CN 113947037A CN 202111028997 A CN202111028997 A CN 202111028997A CN 113947037 A CN113947037 A CN 113947037A
- Authority
- CN
- China
- Prior art keywords
- turbulence intensity
- fitting
- wind speed
- data
- wind turbine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000013461 design Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 15
- 230000002159 abnormal effect Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Primary Health Care (AREA)
- Fluid Mechanics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a wind turbine load calculation method and a wind turbine load calculation device, wherein the method comprises the following steps: acquiring turbulence intensity data of each machine position point of the wind turbine generator at different wind speeds, and extracting enveloping turbulence intensity corresponding to different wind speeds; fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity; and calculating the load value of the wind turbine generator according to the fitting turbulence intensity, and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator. The method comprises the steps of obtaining turbulence intensity data of each machine position of the wind turbine generator at different wind speeds, extracting enveloping turbulence intensity, performing data fitting and calculating a load value of the wind turbine generator, replacing a traditional calculation method to perform safety recheck on each machine position of a wind field, and dividing a maximum rechecking turbulence intensity by using a reasonable wind speed interval, so that load calculation is performed, and loads are more reasonable.
Description
Technical Field
The invention relates to the technical field of wind power control, in particular to a wind turbine load calculation method and device.
Background
When the safety of the large-scale wind generating set is rechecked, the wind resource data of the wind power plant is needed to be used for carrying out load calculation. The formation of the wind resource is a result of the comprehensive action of various complex natural factors, and the evaluation of the turbulence intensity of the wind resource directly influences the magnitude of the ultimate load and the fatigue load of the wind driven generator. Therefore, the method is particularly important for analyzing and evaluating the wind conditions of the wind resources, and based on the prior art, after the wind measurement data is output from wind resource evaluation software, the load calculation is carried out by directly combining various factors such as wind resource turbulence intensity, annual average wind speed and air density of the wind power plant with IEC/GL specifications. The low probability wind speed of an individual sector in software is large, calculation is over conservative, factors such as turbulence data abnormity cause large or even unreasonable load calculation results, and judgment on the safety of the wind power plant is affected.
The existing safety rechecking method excessively depends on wind resource turbulence data output by software, in a wind resource report based on the software output, the wind measuring data mostly takes half a year or more, the probability of occurrence of high wind is low, particularly, the turbulence data of a part of sectors from high wind speed to back are nearly uniform, because the high wind speed and the wind frequency are low, the turbulence data can not be obtained in a short time. When the turbulence data under the wind speed cannot be obtained in the period measured by the wind measuring radar, the software automatically defaults that the turbulence intensity under the wind speed is consistent with the previous wind speed. And combining the wind frequency and the turbulence intensity of each sector, and partial data irrationality exists in the finally obtained effective turbulence intensity value.
Disclosure of Invention
The embodiment of the invention aims to provide a wind turbine load calculation method and device, which are used for extracting enveloping turbulence intensity by acquiring turbulence intensity data of each machine position of a wind turbine at different wind speeds, further performing data fitting and calculating a load value of the wind turbine, replacing the traditional calculation method to perform security recheck on each machine position of a wind field, and dividing a reasonable wind speed interval for maximum rechecking turbulence intensity, so that load calculation is performed, the load is more reasonable, and sudden change or abnormal load at a certain wind speed cannot be generated.
In order to solve the technical problem, a first aspect of the embodiments of the present invention provides a wind turbine load calculation method, including the following steps:
acquiring turbulence intensity data of each machine position point of the wind turbine generator at different wind speeds, and extracting enveloping turbulence intensity corresponding to the different wind speeds respectively;
fitting the turbulence intensity data in a preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity;
and calculating the load value of the wind turbine generator according to the fitting turbulence intensity, and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
Further, fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain a fitting turbulence intensity includes:
acquiring the enveloping turbulence intensity of a first preset wind speed interval, and fitting the turbulence intensity data of the first preset wind speed interval to obtain a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and a tower thereof; and/or
Acquiring the enveloping turbulence intensity of a second preset wind speed interval, and fitting the turbulence intensity data of the second preset wind speed interval to obtain a second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator; and/or
And acquiring the enveloping turbulence intensity of a third preset wind speed interval, and fitting the turbulence intensity data of the third preset wind speed interval to obtain a third fitting turbulence intensity for verifying the fatigue load of the blade.
Further, the first preset wind speed interval is Vave-2 Vave;
the second preset wind speed interval is Vrate-Vout;
the third preset wind speed interval is Vave-2 Vave;
wherein Vave is the annual average wind speed, Vrate is the rated wind speed of the fan, and Vout is the cut-out wind speed of the fan.
Further, after the fitting of the turbulence intensity data in the preset wind speed interval, the method further includes:
and removing the turbulence intensity corresponding to the point where the fitting curve is too high or too low.
Further, after extracting the enveloping turbulence intensity corresponding to the different wind speeds, the method further includes:
and judging whether the enveloping turbulence intensity is in a normal numerical range or not by combining the sector turbulence intensity of the wind turbine generator.
Accordingly, a first aspect of an embodiment of the present invention provides a wind turbine load calculation apparatus, including:
the data acquisition module is used for acquiring turbulence intensity data of each machine site of the wind turbine generator at different wind speeds and extracting enveloping turbulence intensity corresponding to the different wind speeds;
the data fitting module is used for fitting the turbulence intensity data in a preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity;
and the data judgment module is used for calculating the load value of the wind turbine generator according to the fitting turbulence intensity and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
Further, the data fitting module comprises:
the first data fitting unit is used for acquiring the enveloping turbulence intensity of a first preset wind speed interval, fitting the turbulence intensity data of the first preset wind speed interval, and obtaining a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and the tower thereof; and/or
The second data fitting unit is used for acquiring the enveloping turbulence intensity of a second preset wind speed interval, fitting the turbulence intensity data of the second preset wind speed interval, and obtaining second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator; and/or
And the third data fitting unit is used for acquiring the enveloping turbulence intensity of a third preset wind speed interval, fitting the turbulence intensity data of the third preset wind speed interval, and obtaining a third fitting turbulence intensity for verifying the fatigue load of the blade.
Further, the first preset wind speed interval is Vave-2 Vave;
the second preset wind speed interval is Vrate-Vout;
the third preset wind speed interval is Vave-2 Vave;
wherein Vave is the annual average wind speed, Vrate is the rated wind speed of the fan, and Vout is the cut-out wind speed of the fan.
Further, the data fitting module further comprises:
and the data processing unit is used for removing the turbulence intensity corresponding to the point where the fitting curve is too high or too low.
Further, the data acquisition module further comprises:
and the data judgment unit is used for judging whether the enveloping turbulence intensity is positioned in a normal numerical value interval or not by combining the sector turbulence intensity of the wind turbine generator.
Accordingly, a third aspect of the embodiments of the present invention further provides an electronic device, including: at least one processor; and a memory coupled to the at least one processor; the storage stores instructions executable by a processor, and the instructions are executed by the processor to enable at least one processor to execute the wind turbine load calculation method.
In addition, a fourth aspect of the embodiments of the present invention also provides a computer-readable storage medium, on which computer instructions are stored, and the instructions, when executed by a processor, implement the wind turbine load calculation method.
The technical scheme of the embodiment of the invention has the following beneficial technical effects:
the method comprises the steps of obtaining turbulence intensity data of each machine position of the wind turbine generator at different wind speeds, extracting enveloping turbulence intensity, performing data fitting and calculating a load value of the wind turbine generator, replacing a traditional calculation method to perform safety recheck on each machine position of a wind field, and dividing the maximum rechecking turbulence intensity by using a reasonable wind speed interval, so that load calculation is performed, the load is more reasonable, and sudden change or abnormal load at a certain wind speed cannot be generated.
Drawings
FIG. 1 is a flow chart of a wind turbine load calculation method provided by an embodiment of the invention;
FIG. 2 is a logic diagram of a wind turbine load calculation method provided by the embodiment of the invention;
fig. 3 is a schematic diagram of a turbulence intensity (m-4) fit provided by an embodiment of the present invention;
fig. 4 is a schematic diagram of a turbulence intensity (m-1) fit provided by an embodiment of the present invention;
fig. 5 is a schematic diagram of a turbulence intensity (m-10) fit provided by an embodiment of the present invention;
FIG. 6 is a block diagram of a wind turbine load calculating device provided by an embodiment of the invention;
FIG. 7 is a block diagram of a data fitting module provided by an embodiment of the present invention;
fig. 8 is a block diagram of a data determination module according to an embodiment of the present invention.
Reference numerals:
1. the device comprises a data acquisition module, 2, a data fitting module, 21, a first data fitting unit, 22, a second data fitting unit, 23, a third data fitting unit, 24, a data processing unit, 3, a data judgment module, 31 and a data judgment unit.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings in conjunction with the following detailed description. It should be understood that the description is intended to be exemplary only, and is not intended to limit the scope of the present invention. Moreover, in the following description, descriptions of well-known structures and techniques are omitted so as to not unnecessarily obscure the concepts of the present invention.
With the technical development and the improvement of power, the blades used by the large-scale wind generating set are longer and longer, so that the load of the set is more and more sensitive to the aerodynamic performance of wind resource turbulence intensity. Under the high wind speed turbulence intensity of the wind turbine generator, the limit load of the turbine and the components is greatly influenced, and the fatigue load of the turbine and the components is greatly influenced by the turbulence intensity near the rated wind speed.
At present, the specific process of the large-scale wind generating set when carrying out security recheck is as follows: obtaining turbulence intensity values under three different wind speeds, namely m is 1, m is 4 and m is 10, according to the external conditions of the wind resources output by the software; calculating the limit load of the wind driven generator set and the component by combining the maximum turbulence intensity of each machine position point in the category of m-1 with the IEC/GL specification; calculating the fatigue load of the wind generating set by combining the maximum turbulence intensity of each wind speed of each machine position point of which the m is 4 types with the IEC/GL specification; calculating the fatigue load of the blade of the wind driven generator by combining the maximum turbulence intensity of each wind speed of each machine position with the m-10 category and IEC/GL specifications; if the load calculation result is within the deviation range of 5% of the designed load of the wind driven generator, the wind power plant has no potential safety hazard; the load calculation result exceeds the design load of the wind driven generator, the deviation is more than 5%, the potential safety hazard of the wind power plant is large, and protective measures need to be taken in the running process of the unit or the wind driven generator is not put in the unit position.
Referring to fig. 1 and fig. 2, a first aspect of an embodiment of the present invention provides a wind turbine load calculation method, including the following steps:
s100, turbulence intensity data of each machine position point of the wind turbine generator at different wind speeds are obtained, and enveloping turbulence intensity corresponding to the different wind speeds is extracted.
S200, fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity.
And S300, calculating the load value of the wind turbine generator according to the fitting turbulence intensity, and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
Specifically, in step S200, fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain a fitting turbulence intensity, including:
s210, obtaining the enveloping turbulence intensity of the first preset wind speed interval, and fitting the turbulence intensity data of the first preset wind speed interval to obtain a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and the tower thereof.
Specifically, when fatigue loads of the unit and the tower are verified, the enveloping turbulence intensity when m is 4 is extracted by adopting turbulence data of each unit position when m is 4.
And judging the rationality of the envelope turbulence intensity by combining the turbulence intensity of the sector. The wind resource turbulence intensity is defined to know that the turbulence intensity is the ratio of the mean square error sigma of the difference between the instantaneous wind speed and the average wind speed to the average wind speed v, and the occurrence frequency of high wind speed in a anemometry period is low, which may cause that software outputs some effective turbulence data to be abnormal, and the selection of the enveloping turbulence intensity is influenced due to the excessively high or excessively low turbulence intensity.
According to the Weibull wind frequency distribution and the equivalent fatigue load calculation, the following results are obtained: in the Vave (annual average wind speed) and 2Vave wind speed interval, the wind frequency ratio is the highest, the fatigue load is greatly influenced, so that the turbulence data in the interval is taken as the main, the envelope turbulence intensity between the annual average wind speed Vave and 2Vave is taken by adopting a formula of sigma-I (0.75Vhub + b) to fit the turbulence data in the wind speed interval, and the points which are too high or too low on the fitting curve are removed.
Substituting fitting turbulence intensity, calculating fatigue loads of the unit and the tower according to IEC/GL specifications, and judging whether the wind field has potential safety hazards or not by comparing design loads.
And/or
S220, obtaining the enveloping turbulence intensity of the second preset wind speed interval, and fitting the turbulence intensity data of the second preset wind speed interval to obtain a second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator.
Specifically, when the limit loads of the unit and the blade are verified, the enveloping turbulence intensity when m is 1 is extracted by adopting turbulence data of each unit position when m is 1.
And judging the rationality of the envelope turbulence intensity by combining the turbulence intensity of the sector. The wind resource turbulence intensity is defined to know that the turbulence intensity is the ratio of the mean square error sigma of the difference between the instantaneous wind speed and the average wind speed to the average wind speed v, and the occurrence frequency of high wind speed in a anemometry period is low, which may cause that software outputs some effective turbulence data to be abnormal, and the selection of the enveloping turbulence intensity is influenced due to the excessively high or excessively low turbulence intensity.
In a range between a Vrate (rated wind speed of a fan) and a Vout (cut-out wind speed of the fan), extreme loads are easy to occur, so turbulence intensity data in the range between the Vrate and the Vout are selected, an envelope turbulence intensity between the Vrate and the Vout is adopted to fit the turbulence data in the range by adopting a sigma-I (0.75Vhub + b) formula, and points which are too high or too low on a fitting curve are removed.
Substituting fitting turbulence intensity, calculating limit loads of the unit and the blade according to IEC/GL specifications, and judging whether the wind field has potential safety hazards or not by comparing design loads.
And/or
And S230, acquiring the enveloping turbulence intensity of a third preset wind speed interval, and fitting the turbulence intensity data of the third preset wind speed interval to obtain a third fitting turbulence intensity for verifying the fatigue load of the blade.
Specifically, when verifying the blade fatigue load, the envelope turbulence intensity at m 10 is extracted using the turbulence data at each machine position at m 10.
And judging the rationality of the envelope turbulence intensity by combining the turbulence intensity of the sector. The wind resource turbulence intensity is defined to know that the turbulence intensity is the ratio of the mean square error sigma of the difference between the instantaneous wind speed and the average wind speed to the average wind speed v, and the occurrence frequency of high wind speed in a anemometry period is low, which may cause that software outputs some effective turbulence data to be abnormal, and the selection of the enveloping turbulence intensity is influenced due to the excessively high or excessively low turbulence intensity.
According to the Weibull wind frequency distribution and the equivalent fatigue load calculation, the following results are obtained: in the Vave (annual average wind speed) and 2Vave wind speed interval, the wind frequency ratio is the highest, the fatigue load is greatly influenced, so that the turbulence data in the interval is taken as the main, the envelope turbulence intensity between the annual average wind speed Vave and 2Vave is taken by adopting a formula of sigma-I (0.75Vhub + b) to fit the turbulence data in the wind speed interval, and the points which are too high or too low on the fitting curve are removed.
Substituting fitting turbulence intensity, calculating blade fatigue load according to IEC/GL specifications, and comparing design load to judge whether the wind field has potential safety hazards.
In order to make the objects, technical solutions and advantages of the present invention more clear, the present invention is further described in detail below with reference to a flow chart and a security review of wind resource data and model configuration of a certain wind farm. The wind field unit is a 2.2MW/131 unit, the total number of the seven unit sites is seven, the height of a tower is 80m, the cut-in wind speed is 3m/s, the cut-out wind speed is 20m/s, the annual average wind speed is 7m/s, and the rated wind speed of the unit is 9 m/s. It should be understood that the description is intended to be exemplary only, and is not intended to limit the scope of the present invention.
Step 1: and extracting the maximum turbulence intensity at the wind speed, namely the envelope turbulence intensity according to the turbulence intensity values at different wind speeds of the machine points output by the software.
Step 2: and when fatigue loads of the unit and the tower are verified, extracting the enveloping turbulence intensity when m is 4 by adopting turbulence data of each unit position when m is 4.
Table 1: turbulence intensity of each machine position when M is 4
And step 3: and judging the rationality of the envelope turbulence intensity by combining the turbulence intensity of the sector. The wind resource turbulence intensity is defined to know that the turbulence intensity is the ratio of the mean square error sigma of the difference between the instantaneous wind speed and the average wind speed to the average wind speed v, and the occurrence frequency of high wind speed in a anemometry period is low, which may cause that software outputs some effective turbulence data to be abnormal, and the selection of the enveloping turbulence intensity is influenced due to the excessively high or excessively low turbulence intensity.
And 4, step 4: according to the Weibull wind frequency distribution and the equivalent fatigue load calculation, the following results are obtained: in the wind speed interval of 7m/s (annual average wind speed) and 14m/s, the wind frequency ratio is the highest, the fatigue load is greatly influenced, so that the turbulence data in the interval are taken as the main, the envelope turbulence intensity between 7m/s and 14m/s of the annual average wind speed is taken by adopting a formula of sigma I (0.75Vhub + b) to fit the turbulence data in the wind speed interval, and the points which are too high or too low on the fitted curve are removed. As shown in fig. 3, the turbulence intensity for low wind speeds and cut-out wind speeds with less fatigue load effect can be removed and the turbulence intensity in the fitted 7-14m/s wind speed interval is selected as the calculated value.
And 5: substituting the fitting turbulence intensity to calculate the fatigue loads of the unit and the tower according to the IEC/GL standard, and comparing the design load to judge whether the wind field has potential safety hazards.
And 6, when the limit loads of the unit and the blade are verified, extracting the enveloping turbulence intensity when m is 1 by adopting turbulence data of each unit position when m is 1.
Table 2: turbulence intensity of time group when m is 1
And 7: after repeating the step 3, as shown in fig. 4, extreme loads are likely to occur in a wind speed interval between Vrate (rated wind speed of the fan) and Vout (cut-out wind speed of the fan), so turbulence intensity data in the wind speed interval between Vrate and Vout are selected, and the envelope turbulence intensity between Vrate and Vout is fitted to the turbulence data in the wind speed interval by using a formula of σ ═ I (0.75Vhub + b), so that points that are too high or too low on the fitted curve are removed.
And 8: substituting the fitting turbulence intensity to calculate the limit load of the unit and the blade according to the IEC/GL standard, and judging whether the wind field has potential safety hazards or not by comparing the design load.
And step 9: when the fatigue load of the blade is verified, the enveloping turbulence intensity when m is 10 is extracted by adopting turbulence data of each machine position when m is 10.
TABLE 3 time group turbulence intensity when m is 10
Step 10: and (5) repeating the steps 3, 4 and 5 to complete the safety rechecking of the wind power plant unit, the tower and the blade as shown in the figure 5.
Accordingly, referring to fig. 6, a first aspect of the embodiments of the present invention provides a wind turbine load calculating device, including: the device comprises a data acquisition module 1, a data fitting module 2 and a data judgment module 3.
The data acquisition module 1 is used for acquiring turbulence intensity data of each machine site of the wind turbine generator at different wind speeds and extracting enveloping turbulence intensity corresponding to different wind speeds; the data fitting module 2 is used for fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity; the data judgment module 3 is used for calculating the load value of the wind turbine generator according to the fitting turbulence intensity and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
Further, referring to fig. 7, the data fitting module 2 includes: a first data fitting unit 21, a second data fitting unit 22 and/or a third data fitting unit 23.
The first data fitting unit 21 is configured to obtain the enveloping turbulence intensity of the first preset wind speed interval, and fit the turbulence intensity data of the first preset wind speed interval to obtain a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and the tower thereof. And the second data fitting unit 22 is configured to obtain the enveloping turbulence intensity of the second preset wind speed interval, and fit the turbulence intensity data of the second preset wind speed interval to obtain a second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator. The third data fitting unit 23 is configured to obtain an enveloping turbulence intensity of a third preset wind speed interval, and fit turbulence intensity data of the third preset wind speed interval to obtain a third fitting turbulence intensity for verifying the fatigue load of the blade.
Further, the first preset wind speed interval is Vave-2 Vave; the second preset wind speed interval is Vrate-Vout; the third preset wind speed interval is Vave-2 Vave. Wherein Vave is the annual average wind speed, Vrate is the rated wind speed of the fan, and Vout is the cut-out wind speed of the fan.
Further, the data fitting module 2 further includes: a data processing unit 24. The data processing unit 24 is used to remove turbulence intensity corresponding to points where the fitted curve is too high or too low.
Further, referring to fig. 8, the data obtaining module 3 further includes: a data judging unit 31. The data judging unit 31 is configured to judge, in combination with the sector turbulence intensity of the wind turbine, whether the envelope turbulence intensity is within a normal numerical range.
Accordingly, a third aspect of the embodiments of the present invention further provides an electronic device, including: at least one processor; and a memory coupled to the at least one processor; the storage stores instructions executable by a processor, and the instructions are executed by the processor to enable at least one processor to execute the wind turbine load calculation method.
In addition, a fourth aspect of the embodiments of the present invention also provides a computer-readable storage medium, on which computer instructions are stored, and the instructions, when executed by a processor, implement the wind turbine load calculation method.
The embodiment of the invention aims to protect a wind turbine load calculation method and a wind turbine load calculation device, wherein the method comprises the following steps: acquiring turbulence intensity data of each machine position point of the wind turbine generator at different wind speeds, and extracting enveloping turbulence intensity corresponding to different wind speeds; fitting the turbulence intensity data in the preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity; and calculating the load value of the wind turbine generator according to the fitting turbulence intensity, and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator. The technical scheme has the following effects:
the method comprises the steps of obtaining turbulence intensity data of each machine position of the wind turbine generator at different wind speeds, extracting enveloping turbulence intensity, performing data fitting and calculating a load value of the wind turbine generator, replacing a traditional calculation method to perform safety recheck on each machine position of a wind field, and dividing the maximum rechecking turbulence intensity by using a reasonable wind speed interval, so that load calculation is performed, the load is more reasonable, and sudden change or abnormal load at a certain wind speed cannot be generated.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (10)
1. A wind turbine load calculation method is characterized by comprising the following steps:
acquiring turbulence intensity data of each machine position point of the wind turbine generator at different wind speeds, and extracting enveloping turbulence intensity corresponding to the different wind speeds respectively;
fitting the turbulence intensity data in a preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity;
and calculating the load value of the wind turbine generator according to the fitting turbulence intensity, and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
2. The wind turbine load calculation method according to claim 1, wherein fitting the turbulence intensity data in a preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain a fitting turbulence intensity comprises:
acquiring the enveloping turbulence intensity of a first preset wind speed interval, and fitting the turbulence intensity data of the first preset wind speed interval to obtain a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and a tower thereof; and/or
Acquiring the enveloping turbulence intensity of a second preset wind speed interval, and fitting the turbulence intensity data of the second preset wind speed interval to obtain a second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator; and/or
And acquiring the enveloping turbulence intensity of a third preset wind speed interval, and fitting the turbulence intensity data of the third preset wind speed interval to obtain a third fitting turbulence intensity for verifying the fatigue load of the blade.
3. The wind turbine load calculation method according to claim 2,
the first preset wind speed interval is Vave-2 Vave;
the second preset wind speed interval is Vrate-Vout;
the third preset wind speed interval is Vave-2 Vave;
wherein Vave is the annual average wind speed, Vrate is the rated wind speed of the fan, and Vout is the cut-out wind speed of the fan.
4. The wind turbine load calculation method according to any one of claims 2, wherein after fitting the turbulence intensity data in the preset wind speed interval, the method further comprises:
and removing the turbulence intensity corresponding to the point where the fitting curve is too high or too low.
5. The wind turbine load calculation method according to any one of claims 1 to 4, wherein after extracting the enveloping turbulence intensity corresponding to the different wind speeds, the method further comprises:
and judging whether the enveloping turbulence intensity is in a normal numerical range or not by combining the sector turbulence intensity of the wind turbine generator.
6. A wind turbine load calculation device, comprising:
the data acquisition module is used for acquiring turbulence intensity data of each machine site of the wind turbine generator at different wind speeds and extracting enveloping turbulence intensity corresponding to the different wind speeds;
the data fitting module is used for fitting the turbulence intensity data in a preset wind speed interval according to the enveloping turbulence intensity of the preset wind speed interval to obtain fitting turbulence intensity;
and the data judgment module is used for calculating the load value of the wind turbine generator according to the fitting turbulence intensity and judging whether the load value is in a normal range or not by combining the design load of the wind turbine generator.
7. The wind turbine load calculation device of claim 6, wherein the data fitting module comprises:
the first data fitting unit is used for acquiring the enveloping turbulence intensity of a first preset wind speed interval, fitting the turbulence intensity data of the first preset wind speed interval, and obtaining a first fitting turbulence intensity for verifying the fatigue load of the wind turbine generator and the tower thereof; and/or
The second data fitting unit is used for acquiring the enveloping turbulence intensity of a second preset wind speed interval, fitting the turbulence intensity data of the second preset wind speed interval, and obtaining second fitting turbulence intensity for verifying the wind turbine generator and the ultimate load of the blades of the wind turbine generator; and/or
And the third data fitting unit is used for acquiring the enveloping turbulence intensity of a third preset wind speed interval, fitting the turbulence intensity data of the third preset wind speed interval, and obtaining a third fitting turbulence intensity for verifying the fatigue load of the blade.
8. The wind turbine load calculation device of claim 7,
the first preset wind speed interval is Vave-2 Vave;
the second preset wind speed interval is Vrate-Vout;
the third preset wind speed interval is Vave-2 Vave;
wherein Vave is the annual average wind speed, Vrate is the rated wind speed of the fan, and Vout is the cut-out wind speed of the fan.
9. The wind turbine load calculation device of claim 7, wherein the data fitting module further comprises:
and the data processing unit is used for removing the turbulence intensity corresponding to the point where the fitting curve is too high or too low.
10. The wind turbine load calculation device according to any one of claims 6 to 9, wherein the data acquisition module further comprises:
and the data judgment unit is used for judging whether the enveloping turbulence intensity is positioned in a normal numerical value interval or not by combining the sector turbulence intensity of the wind turbine generator.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111028997.8A CN113947037B (en) | 2021-09-03 | 2021-09-03 | A method and device for calculating load of wind turbine generator set |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111028997.8A CN113947037B (en) | 2021-09-03 | 2021-09-03 | A method and device for calculating load of wind turbine generator set |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113947037A true CN113947037A (en) | 2022-01-18 |
| CN113947037B CN113947037B (en) | 2024-12-13 |
Family
ID=79327786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111028997.8A Active CN113947037B (en) | 2021-09-03 | 2021-09-03 | A method and device for calculating load of wind turbine generator set |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113947037B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114856931A (en) * | 2022-03-30 | 2022-08-05 | 许昌许继风电科技有限公司 | Method, system and electronic device for evaluating fatigue load of wind turbine |
| CN114912756A (en) * | 2022-04-13 | 2022-08-16 | 浙江运达风电股份有限公司 | Net-based wind generating set fatigue load fine evaluation method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760669A (en) * | 2016-02-15 | 2016-07-13 | 重庆三峡学院 | Wind power converter power module fault rate assessment method and system |
| CN112528504A (en) * | 2020-12-14 | 2021-03-19 | 中国船舶重工集团海装风电股份有限公司 | Wind turbine generator fatigue load calculation method based on turbulent flow distribution |
-
2021
- 2021-09-03 CN CN202111028997.8A patent/CN113947037B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105760669A (en) * | 2016-02-15 | 2016-07-13 | 重庆三峡学院 | Wind power converter power module fault rate assessment method and system |
| CN112528504A (en) * | 2020-12-14 | 2021-03-19 | 中国船舶重工集团海装风电股份有限公司 | Wind turbine generator fatigue load calculation method based on turbulent flow distribution |
Non-Patent Citations (3)
| Title |
|---|
| 付德义;薛扬;焦渤;边伟;李松迪;: "湍流强度对风电机组疲劳等效载荷的影响", 华北电力大学学报(自然科学版), no. 01, 30 January 2015 (2015-01-30), pages 45 - 50 * |
| 戴靠山 等: "风力发电机台风载荷响应分析", 《结构工程师》, vol. 31, no. 6, 31 December 2015 (2015-12-31), pages 98 - 106 * |
| 蒋红武 等: "风电场湍流强度对机组结构疲劳安全性的研究", 《应用能源技术》, no. 2, 29 February 2020 (2020-02-29), pages 45 - 50 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114856931A (en) * | 2022-03-30 | 2022-08-05 | 许昌许继风电科技有限公司 | Method, system and electronic device for evaluating fatigue load of wind turbine |
| CN114912756A (en) * | 2022-04-13 | 2022-08-16 | 浙江运达风电股份有限公司 | Net-based wind generating set fatigue load fine evaluation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113947037B (en) | 2024-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109296506B (en) | Vibration detection method, control method and device for wind turbine | |
| CN106815771B (en) | A long-term assessment method for wind farm loads | |
| CN103825521B (en) | Drive system and the discrimination method of generator parameter in a kind of double-fed fan motor unit | |
| CN113947037B (en) | A method and device for calculating load of wind turbine generator set | |
| CN108695862B (en) | Power grid inertia characteristic online evaluation method based on PMU measured data | |
| CN110390146B (en) | Method for calculating fatigue damage of welding seam of tower cylinder of wind turbine generator based on sector load and product | |
| CN110414102A (en) | Fatigue Analysis Method for Offshore Wind Power Structure | |
| CN105508149A (en) | Fault detection method and device for wind generating set | |
| CN110797976A (en) | A test platform and test method suitable for wind farm power controller evaluation | |
| CN113931806A (en) | Wind generating set and control method, controller and control system thereof | |
| CN104794287A (en) | Wind engineering wake flow calculation method | |
| CN113586336B (en) | Control method and control device of wind generating set and computer readable storage medium | |
| CN111287911A (en) | A kind of early warning method and system for wind turbine fatigue load | |
| CN107038264B (en) | Sector division method and system for wind turbine generator | |
| CN114778116A (en) | Method and system for early warning of wind turbine pitch bearing fault | |
| CN112287621B (en) | Wind turbine generator operating state threshold curve determining method, evaluation method and system | |
| CN117272869B (en) | Full wake analysis method considering characteristics of near wake and far wake of wind turbine | |
| CN113847212A (en) | A method for monitoring the natural frequency of wind turbine blades | |
| CN116644618B (en) | A fatigue life assessment method, system, equipment and medium for wind turbines | |
| CN112555101B (en) | Method and device for identifying impeller aerodynamic state of wind generating set | |
| CN115898787A (en) | Method and device for dynamically identifying static yaw error of wind turbine generator | |
| CN116776064A (en) | Method for calculating generating capacity of fan | |
| CN111985104B (en) | Wind driven generator equivalent fatigue load estimation method and device and computer readable storage medium | |
| CN110661261B (en) | PMU data processing method and device for low-frequency oscillation online analysis | |
| CN111211556A (en) | A reliability assessment method for distribution network power supply considering wind power |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |