CN115496391A - Method and device for calculating blade ice melting time of wind driven generator and storage medium - Google Patents
Method and device for calculating blade ice melting time of wind driven generator and storage medium Download PDFInfo
- Publication number
- CN115496391A CN115496391A CN202211226474.9A CN202211226474A CN115496391A CN 115496391 A CN115496391 A CN 115496391A CN 202211226474 A CN202211226474 A CN 202211226474A CN 115496391 A CN115496391 A CN 115496391A
- Authority
- CN
- China
- Prior art keywords
- blade
- ice
- icing
- parameter
- calculating
- 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.)
- Pending
Links
- 238000002844 melting Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 68
- 230000008018 melting Effects 0.000 title claims abstract description 62
- 230000008569 process Effects 0.000 claims abstract description 35
- 238000010521 absorption reaction Methods 0.000 claims abstract description 22
- 238000002474 experimental method Methods 0.000 claims abstract description 21
- 238000010309 melting process Methods 0.000 claims abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 161
- 238000000576 coating method Methods 0.000 claims description 30
- 239000011248 coating agent Substances 0.000 claims description 28
- 238000004364 calculation method Methods 0.000 claims description 22
- 101100425538 Pseudomonas aeruginosa (strain UCBPP-PA14) tis1 gene Proteins 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- 238000012423 maintenance Methods 0.000 abstract description 6
- 230000005611 electricity Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06311—Scheduling, planning or task assignment for a person or group
- G06Q10/063114—Status monitoring or status determination for a person or group
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
- 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/72—Wind turbines with rotation axis in wind direction
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Wind Motors (AREA)
Abstract
The disclosure provides a method, a device and a storage medium for calculating the ice melting time of a blade of a wind driven generator, wherein the method comprises the following steps: according to a first meteorological parameter, a blade parameter and a first experiment coefficient in an icing period, calculating a first heat quantity released in the icing process of the surface of a wind driven generator blade in a unit area, from the starting time of an ice melting period, calculating a second heat quantity absorbed in the ice melting process of the surface of the wind driven generator blade in the unit area according to a second meteorological parameter, the blade parameter and a second experiment coefficient in the ice melting period, determining the ending time of the ice melting period under the condition that the second heat quantity is equal to the first heat quantity, calculating the heat release quantity when the blade is iced and the heat absorption quantity when the blade is iced according to the real-time collected meteorological parameters, the blade parameters and the experiment coefficients, and deducing the ending time of the ice melting based on a thermodynamic balance principle, so that operation and maintenance personnel can start up the machine to generate electricity in time according to the ending time, and the safe and economic operation of a wind power plant unit and a power grid is facilitated to be guaranteed.
Description
Technical Field
The disclosure relates to the technical field of wind power generation, in particular to a method and a device for calculating the de-icing time of a blade of a wind driven generator and a storage medium.
Background
In winter, high mountain wind farms face severe blade icing problems. The icing of the blades of the wind driven generator can cause the change of the aerodynamic performance of the blades, the output power and the safe operation of a unit are seriously influenced, the output power of the whole wind driven generator unit is lost, and the problems of potential safety hazards such as blade stall, mechanical failure, icing falling and the like are easily caused. When a large amount of ice is coated on the surface of the fan blade, the unbalanced load of the unit is increased due to different ice loads on each blade, so that the service life of parts of the fan is shortened, and the unit is greatly damaged; due to different icing thicknesses, the original pneumatic appearance of the blade is changed, and the wind energy utilization coefficient of the wind turbine generator is reduced, so that the generated power is reduced, and the wind turbine generator cannot be normally started in serious cases; after ice coating, if the fan continues to operate, pieces of the thrown ice layer or falling ice cubes may injure the fan itself and people or things nearby. Therefore, the accurate state prediction of the icing and deicing of the blades of the wind driven generator has very important practical significance on the safety, reliability and economy of the wind turbine generator.
Disclosure of Invention
The disclosure provides a method and a device for calculating the blade ice melting time of a wind driven generator and a storage medium, and aims to solve at least one of technical problems in the related art to a certain extent.
The embodiment of the first aspect of the disclosure provides a method for calculating the blade ice melting time of a wind driven generator, which includes: calculating a first heat quantity released in the icing process of the surface of the blade of the wind driven generator in unit area according to a first meteorological parameter, the blade parameter and a first experiment coefficient in the icing period; calculating a second heat absorbed in the process of melting ice on the surface of the blade of the wind driven generator in unit area from the start time of the ice melting period according to a second meteorological parameter, the blade parameter and a second experimental coefficient in the ice melting period; and determining an end time of the ice-melt cycle if the second amount of heat equals the first amount of heat.
An embodiment of a second aspect of the present disclosure provides a device for calculating a blade ice melting time of a wind turbine, including: the first calculation module is used for calculating first heat released in the icing process of the surface of the wind driven generator blade in unit area according to a first meteorological parameter, a blade parameter and a first experiment coefficient in the icing period; the second calculation module is used for calculating second heat absorbed in the ice melting process of the surface of the blade of the wind driven generator in unit area from the start time of the ice melting period according to a second meteorological parameter, the blade parameter and a second experimental coefficient in the ice melting period; and a determining module for determining an end time of the ice-melting period when the second amount of heat is equal to the first amount of heat.
An embodiment of a third aspect of the present disclosure provides a computer device, including: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of wind turbine blade ice melting time calculation of an embodiment of the present disclosure.
A fourth aspect of the present disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the method for calculating the blade ice melting time of a wind turbine generator disclosed in the embodiments of the present disclosure.
In the embodiment, the first heat released in the process of icing the surface of the blade of the wind driven generator in unit area is calculated according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, the second heat absorbed in the process of deicing the surface of the blade of the wind driven generator in unit area is calculated according to the second meteorological parameter, the blade parameter and the second experimental coefficient in the deicing period from the starting time of the deicing period, the ending time of the deicing period is determined under the condition that the second heat is equal to the first heat, the heat release quantity when the blade is iced and the heat absorption quantity when the blade is deicing can be calculated according to the meteorological parameter, the blade parameter and the experimental coefficient which are collected in real time, and the ending time of the deicing is deduced based on the thermodynamic balance principle, so that operation and maintenance personnel can start the generator to generate electricity in time according to the ending time, and the safe and economic operation of a wind power plant unit and a power grid is facilitated to be guaranteed.
Additional aspects and advantages of the disclosure will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the disclosure.
Drawings
The above and/or additional aspects and advantages of the present disclosure will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a schematic flow chart diagram of a method for calculating blade ice melting time of a wind turbine provided according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram illustrating heat absorption and heat release relationships during blade icing and de-icing of a wind turbine provided according to an embodiment of the present disclosure;
FIG. 3 is a schematic flow chart diagram of a method for calculating blade ice melting time of a wind turbine generator according to another embodiment of the present disclosure;
FIG. 4 is a schematic view of an overall calculation flow of blade ice melting time of a wind turbine provided in an embodiment of the present disclosure;
FIG. 5 is a schematic diagram of a wind turbine blade ice-melt time calculation apparatus provided in accordance with another embodiment of the present disclosure;
FIG. 6 illustrates a block diagram of an exemplary computer device suitable for use in implementing embodiments of the present disclosure.
Detailed Description
Reference will now be made in detail to the embodiments of the present disclosure, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the drawings are exemplary only for the purpose of illustrating the present disclosure and should not be construed as limiting the same. On the contrary, the embodiments of the disclosure include all changes, modifications and equivalents coming within the spirit and terms of the claims appended hereto.
It should be noted that an execution subject of the method for calculating the blade-melting time of the wind turbine generator according to this embodiment may be a device for calculating the blade-melting time of the wind turbine generator, where the device may be implemented in a software and/or hardware manner, and the device may be configured in an electronic device, and the electronic device may include, but is not limited to, a terminal, a server, and the like.
Fig. 1 is a schematic flow chart of a method for calculating blade ice melting time of a wind turbine generator provided according to an embodiment of the present disclosure, as shown in fig. 1, the method includes:
s101: and calculating a first heat quantity released in the process of icing the surface of the blade of the wind driven generator in unit area according to the first meteorological parameter, the blade parameter and the first experiment coefficient in the icing period.
Wherein, the time period of the icing process (also called icing) on the surface of the wind turbine blade can be called an icing period, which can be T 1 [t1，t2]T1 represents the start time of ice coating, and may be the time when the ambient temperature is 0 ℃, t2 represents the end time of ice coating, and may be the time when the ambient temperature returns to 0 ℃ again from 0 ℃ or lower after t 1.
The first meteorological parameter is a meteorological parameter acquired in real time in an icing period, for example, when the icing condition is reached, the meteorological parameter at the height of the wind driven generator blade can be acquired to obtain the first meteorological parameter. Wherein the first meteorological parameter comprises, for example, an air relative humidity value R H1 Real-time wind speed v w1 Solar irradiance TIS1, water vapor content in air M v Temperature, and any other weather parameters as may be desired, without limitation.
It will be appreciated that the icing cycle is a longer time period and therefore the first weather parameter may be different at different points in time during the icing cycle. In this case, embodiments of the present disclosure may determine a plurality of time periods within an icing cycle, which may be referred to as a first time period t s1 A first time period t s1 In the case of sufficiently small values, the gas phase relative humidity value R in this time period can be considered H1 Real-time wind speed v w1 Blade surface air flow velocity v 0 The first weather parameters, such as solar irradiance TIS1, are constant, that is, embodiments of the present disclosure may operate at each first time period t s1 First meteorological parameters are acquired in real time.
And blade parameters, e.g. including linear speed v of blade rotation c1 Windward area S of blade m And any other possible blade-related parameters, without limitation. According to the embodiment, the first meteorological parameter can be acquired, the rotating speed data of the blades can be read from an SCADA (supervisory control and data acquisition) system of the wind power plant, and the rotating linear speed v of the blades can be determined by combining the wing profile data of the blades of the unit c1 Windward area S of blade m . It will be appreciated that the first time periods t are different s1 The blade parameters within may be the same, or different.
The first experimental coefficient is parameter or common sense data obtained by experiment, such as water vapor collection coefficient c W The retention amount of water vapor in the air on the surface of the blade can be determined; density of ice ρ i (ii) a At the same temperature, the water vapor reaches the absolute humidity rho of the saturated gas W (ii) a Coefficient of influence of solar irradiance on icing
And any other possible experimental coefficients, without limitation.
As is well known, the process of water icing requires heat to be released, and the embodiment of the disclosure can calculate the heat released in the process of icing the surface of the blade of the wind driven generator per unit area according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, wherein the heat can be called as first heat, and Q can be used for calculating the heat Put And (4) showing. FIG. 2 is a schematic diagram of heat absorption and heat release relationships during blade icing and de-icing of a wind turbine provided according to an embodiment of the disclosure, as shown in FIG. 2, Q Put Is the icing period [ t1, t2 ]]And during the period, the total heat released in the process of coating ice on the surface of the blade of the wind driven generator per unit area.
In some embodiments, before calculating the first heat amount released in the process of icing on the surface of the blade, meteorological prediction data of a predetermined time period may be further obtained, and whether the predetermined time period meets an icing condition is determined according to the meteorological prediction data, that is, whether the current environment meets the icing condition is determined, and the process of calculating the icing time is required to be performed when the icing condition is met.
Specifically, the embodiment of the present disclosure may obtain weather prediction data of a predetermined time period, for example, obtain weather information periodically collected and published by a national weather department, including weather prediction data of 24 hours in the future and longer periods, including weather prediction data of temperature, humidity, wind speed, and the like, and determine whether ice-covered weather conditions are met for 24 hours in the future and longer periods according to conditions such as soft rime formation and soft rime formation. Wherein the judgment is based on the actual temperature T being less than 0 deg.C and the humidity R H Greater than R 0 (different wind speeds, R 0 Different values) were considered to have reached icing conditions. If the meteorological condition reaches the blade surface icing condition in the future 24 hours, the blade icing and de-icing time predicting and calculating system (namely, the wind power plant calculating system executing the wind driven generator blade de-icing time calculating method) executes the icing and de-icing related programs; otherwise, the system continues to poll and judge the meteorological data at regular time.
When the predicted meteorological conditions reach the icing conditions on the blade surface, the system collects meteorological parameters such as air relative humidity value, real-time wind speed, solar irradiance, water vapor content in the air, temperature and the like at the height of the blade of the wind driven generator, judges whether the icing conditions are met or not in real time, and takes the meteorological parameters collected in real time as the first meteorological parameters and executes the calculation method of the ice melting time of the blade of the wind driven generator after the icing conditions are reached.
S102: and calculating the second heat absorbed in the process of melting ice on the surface of the blade of the wind driven generator in unit area according to the second meteorological parameter, the blade parameter and the second experimental coefficient in the ice melting period from the starting time of the ice melting period.
The time period of the ice melting process on the surface of the wind turbine blade can be referred to as an ice melting period. As shown in FIG. 2, the ice-melt period may be T 2 [t3，t4]T3 represents the start time of ice melting, mayTaking the moment when the ambient temperature is 0 ℃ (i.e. the ice coating end time t2 is the same as the ice melting start time t 3), t4 represents the end time of ice melting, namely: the time required by the embodiments of the present disclosure.
The second meteorological parameter is a meteorological parameter collected in real time from the start time of the ice melting cycle, namely: meteorological parameters acquired in real time starting from t3, including, for example, the relative humidity value R of the air H2 Blade surface air flow velocity v 1 (similarly, the solar irradiance TIS2, the difference Δ t between the air temperature and the ice temperature, and any other possible parameters may be calculated according to the real-time wind speed and the linear speed of blade rotation acquired in the ice melting period, which is not limited in this respect.
It will be appreciated that the ice-melt period is also a longer period of time, and therefore the second meteorological parameters may be different at different points in time within the ice-melt period. In this case, embodiments of the present disclosure may determine multiple time periods within the ice melt cycle, which may be referred to as a second time period t s2 A second time period t s2 In the case of sufficiently small values, the relative humidity value R of the air in this time period can be considered H2 Blade surface air flow velocity v 1 The solar irradiance TIS2, the air temperature and ice temperature difference Δ t, etc. are constant, that is, the disclosed embodiments may begin at t3 and begin at every second time period t s2 And acquiring a second meteorological parameter in real time.
And blade parameters, e.g. including the frontal area S of the blade m And any other possible blade-related parameters, without limitation. Similarly, in the embodiment, the rotation speed data of the blade can be read from the SCADA system of the wind power plant while the second meteorological parameter is acquired, and the windward area S of the blade is determined by combining the airfoil profile data of the blade of the unit m . It will be appreciated that the blade parameters may be the same, or different, during different second time periods.
The second experiment coefficient is parameter or common knowledge data obtained through experiments, for example, the second experiment coefficient comprises a correction coefficient beta, an ice surface heat exchange coefficient h and a solar irradiance influence coefficient during ice melting
And any other possible experimental coefficients, without limitation.
As is well known, the ice melting process needs to absorb heat, and the embodiment of the present disclosure may cumulatively calculate, from the ice melting start time t3, the heat that has been absorbed in the ice melting process on the surface of the wind turbine blade per unit area according to the second meteorological parameter, the blade parameter and the second experimental coefficient collected in real time, where the heat may be referred to as a second heat, which may be Q Suction device Shown, as in FIG. 2, Q Suction device Is the ice melting period [ t3, t4 ]]In the period, the heat absorbed in the process of melting ice on the surface of the blade of the wind driven generator in unit area, namely, the second heat Q along with the lapse of the ice melting time Suction device And also increases.
S103: in the event that the second amount of heat equals the first amount of heat, an end time of the ice melt cycle is determined.
The embodiment of the disclosure can monitor the second heat Q in real time Suction device Until the first heat quantity Q Placing the Are equal. At the second heat quantity Q, as shown in FIG. 2 Suction device Is equal to the first heat quantity Q Put In the case of (a): and if the heat absorbed and released by the ice melting and the ice coating is the same, the ice coating on the surface of the blade of the wind driven generator is completely melted, the moment is the ending time t4 of the ice melting period, and the blade of the wind driven generator can be started.
In the embodiment, the first heat released in the process of icing the surface of the blade of the wind driven generator in unit area is calculated according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, the second heat absorbed in the process of deicing the surface of the blade of the wind driven generator in unit area is calculated according to the second meteorological parameter, the blade parameter and the second experimental coefficient in the deicing period from the starting time of the deicing period, the ending time of the deicing period is determined under the condition that the second heat is equal to the first heat, the heat release quantity when the blade is iced and the heat absorption quantity when the blade is deicing can be calculated according to the meteorological parameter, the blade parameter and the experimental coefficient which are collected in real time, and the ending time of the deicing is deduced based on the thermodynamic balance principle, so that operation and maintenance personnel can start the generator to generate electricity in time according to the ending time, and the safe and economic operation of a wind power plant unit and a power grid is facilitated to be guaranteed.
FIG. 3 is a schematic flow chart of a method for calculating blade ice melting time of a wind turbine generator according to another embodiment of the present disclosure, as shown in FIG. 3, the method includes:
s301: and calculating the icing mass of the icing on the surface of the wind driven generator blade in unit area according to the first meteorological parameter, the blade parameter and the first experiment coefficient.
According to the embodiment of the disclosure, the first heat Q released in the ice coating period is calculated Placing the In the process, first, a first meteorological parameter (e.g., air relative humidity value R) may be determined H1 Real-time wind speed v w1 Solar irradiance TIS1, water vapor content in air M v ) Blade parameter (e.g. linear velocity v of blade rotation) c1 Windward area S of blade m ) And a first experimental coefficient (e.g., water vapor collection coefficient c) W Density of ice ρ i The water vapor reaches the absolute humidity rho of the saturated gas at the same temperature W Solar irradiance influence coefficient during ice coating
) Calculating the icing mass of the icing on the surface of the wind turbine blade per unit area, which can be M General assembly And (4) showing.
In some embodiments, the icing thickness of the icing on the surface of the wind turbine blade per unit area may be calculated according to the first meteorological parameter, the blade parameter and the first experimental coefficient, that is: total thickness of ice coating, which may be delta General (1) Represents; further, according to the thickness of ice coating delta General assembly S per unit area and ice density ρ i Calculating the icing mass M General (1) The calculation formula is as follows:
M general (1) ＝δ General assembly *S*ρ i
Some embodiments, the ice coating thickness delta is calculated General assembly In the process of (2), firstly, a plurality of first time periods t in the icing period can be determined s1 To (1)Calculating weather parameters, blade parameters and first experiment coefficients, and calculating the surface of the wind driven generator blade in unit area in a plurality of first time periods t s1 Multiple ice layer thicknesses of the ice coating, namely: each first time period t s1 The inner ice layer thickness, which can be represented by δ, wherein the calculation formula of the ice layer thickness δ is as follows:
wherein, t s1 Representing a first time period duration; rho i Represents the density of ice; c. C W Representing the water vapor collection coefficient; α represents a correction coefficient; m v Representing the moisture content of the air during a first time period; r H1 Representing a relative humidity value of the air during a first time period; v. of 0 Representing the blade surface air flow rate during a first time period,
v w1 representing real-time wind speed, v, over a first period of time c1 Representing the linear speed of blade rotation during a first time period; TIS1 represents solar irradiance during a first time period;
representing solar irradiance influence coefficient during ice coating; s m Representing the windward area of the blade; rho W Indicating that the water vapor reaches the absolute humidity of the saturated gas at the same temperature. Wherein R is H1 、v w1 TIS1 belongs to a first meteorological parameter, v c1 、S m Belonging to blade parameters, ρ i 、α、ρ W 、c W 、
Belonging to the first experimental coefficient.
Further, the ice layer thicknesses for each first time period are summed to obtain a total ice coating thickness δ General assembly The calculation formula is as follows:
where ε represents an empirical correction factor, which can be obtained from a large number of experimental data.
S302: the first heat is calculated based on the mass of the ice coating, the specific heat capacity of the ice, and the amount of temperature change at which water vapor condenses into ice from 0 ℃.
The above determination of icing Mass M General (1) Then, further, the ice coating quality M can be determined General assembly The specific heat capacity C of the ice and the temperature change delta t of the water vapor condensed into the ice from 0 ℃ are calculated, and the first heat Q is calculated Put The calculation formula is as follows:
Q placing the ＝M General assembly *C*Δt
S303: and calculating the instantaneous heat absorption capacity of the blade surface of the wind driven generator in unit area in the ice melting process of the second time period according to the second meteorological parameters, the blade parameters and the second experimental coefficient of the second time period which has passed in the ice melting cycle.
The second heat Q absorbed by the embodiment of the disclosure in the process of calculating the ice melting process of the surface of the wind driven generator blade in unit area Suction device In the process of (2), firstly, the ice melting period [ t3, t4 ] can be determined]Within the second time period t s2 The second meteorological parameters, the blade parameters and the second experimental coefficients, and the instantaneous heat absorption capacity in the ice melting process of each second time period is calculated, wherein the instantaneous heat absorption capacity can be Q Instantaneous suction Expressed, the calculation formula is as follows:
wherein Q Instantaneous suction Represents the instantaneous heat absorption; β represents a correction coefficient; r is H2 Representing the air-to-humidity value during the second time period; v. of 1 Representing a blade surface air flow rate for a second time period; TIS2 represents solar irradiance during a second time period; s. the m Representing the windward area of the blade; t is t s2 Representing a second time period duration; h represents the heat exchange coefficient of the ice surface; Δ t represents the air temperature and ice in the second time periodA temperature difference;
representing the solar irradiance influence coefficient during ice melting, wherein R H2 、v 1 TIS2, deltat belong to a second meteorological parameter, S m Belongs to the blade parameters of beta, h,
Belonging to the second experimental coefficient.
S304: and accumulating the instantaneous heat absorption quantity to obtain a second heat quantity.
Namely: the second time period t will have elapsed s2 Instantaneous heat absorption quantity Q Instantaneous suction Accumulating to obtain the second heat Q absorbed currently Suction device Up to accumulated Q Suction device Is equal to the first heat quantity Q Placing the The calculation formula is as follows:
therefore, the ice melting and ice coating periods can be divided into a plurality of time periods respectively to calculate the first heat and the second heat, and therefore the accuracy of the calculation result can be improved.
S305: in the event that the second amount of heat equals the first amount of heat, an end time of the ice melt cycle is determined.
For specific description of S305, refer to the foregoing embodiments, which are not described herein again.
In some embodiments, the calculation result can be compared with the actual on-site ice melting time, and the beta value is corrected, so that the calculation result is more accurate.
In the embodiment, the first heat released in the icing process of the blade surface of the wind driven generator in unit area is calculated according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, the second heat absorbed in the deicing process of the blade surface of the wind driven generator in unit area is calculated according to the second meteorological parameter, the blade parameter and the second experimental coefficient in the deicing period from the starting time of the deicing period, the ending time of the deicing period is determined under the condition that the second heat is equal to the first heat, the heat release amount when the blade is iced and the heat absorption amount when the blade is iced can be calculated according to the meteorological parameters, the blade parameters and the experimental coefficients collected in real time, and the ending time of the deicing is deduced based on the thermodynamic equilibrium principle, so that operation and maintenance personnel can start the generator generation in time according to the ending time, and the safe and economic operation of a wind power plant unit and a power grid is facilitated to be ensured. In addition, the ice melting and ice coating periods can be divided into a plurality of time periods respectively to calculate the first heat and the second heat, so that the accuracy of the calculation result can be improved.
In a specific example, fig. 4 is a schematic view of an overall calculation flow of blade ice melting time of a wind turbine provided according to an embodiment of the present disclosure, and as shown in fig. 4, a calculation process of blade ice melting time of a wind turbine is as follows:
acquiring meteorological prediction data; according to the acquired weather prediction data, determine whether icing weather conditions are met for 24 hours in the future? If so, starting a blade icing and de-icing time prediction and calculation system, acquiring wind speed, temperature, humidity and solar irradiance data at the height of the fan impeller in real time, reading impeller rotating speed data from an SCADA (supervisory control and data acquisition) system, and acquiring unit blade profile data, otherwise, circularly acquiring meteorological data; further, determine in real time whether the weather conditions for icing are present? If so, calculating the accumulated thickness of the ice coating on the surface of the blade in unit area in the period T, further calculating the ice coating quality on the surface of the blade in unit area in the period T, and further calculating the heat Q released in the ice coating process on the surface of the blade in unit area in the time period T1 when the environmental temperature is lower than 0 DEG C Put Further establishing a function relation between the heat absorbed in the ice coating process on the unit area of the blade surface in the time period T2 when the ambient temperature is higher than 0 ℃ and the specific heat, the air humidity, the temperature, the air flow velocity flowing through the blade surface, the solar irradiance and the accumulated time of the ice, and calculating the total heat absorption capacity Q at a certain moment Suction device (ii) a Further calculating and correcting ice melting time in real time according to a heat balance principle; and further, after the natural ice melting is finished, operation and maintenance personnel are reminded to start the wind turbine in time.
FIG. 5 is a schematic diagram of a device for calculating blade ice melting time of a wind turbine provided according to another embodiment of the present disclosure. As shown in FIG. 5, the wind turbine blade ice melting time calculation device 50 includes:
the first calculating module 501 is used for calculating a first heat quantity released in the icing process of the surface of the wind driven generator blade in unit area according to a first meteorological parameter, a blade parameter and a first experiment coefficient in the icing period;
the second calculating module 502 is configured to calculate, from the start time of the ice-melting period, a second heat quantity that has been absorbed in the ice-melting process of the blade surface of the wind turbine generator in a unit area according to a second meteorological parameter, the blade parameter, and a second experimental coefficient in the ice-melting period; and
a determining module 503, configured to determine the end time of the ice-melting period when the second heat amount is equal to the first heat amount.
In some embodiments, the first calculating module 501 is specifically configured to: calculating the icing mass of the icing on the surface of the wind driven generator blade in unit area according to the first weather parameter, the blade parameter and the first experiment coefficient; and calculating the first heat according to the ice coating quality, the specific heat capacity of the ice and the temperature variation of the water vapor condensed into the ice from 0 ℃.
In some embodiments, the first calculating module 501 is specifically configured to: calculating the icing thickness of the icing on the surface of the wind driven generator blade in unit area according to the first meteorological parameter, the blade parameter and the first experiment coefficient; and calculating the icing mass according to the icing thickness, the unit area and the ice density.
In some embodiments, the first calculating module 501 is specifically configured to: according to the first meteorological parameters, the blade parameters and the first experiment coefficients of a plurality of first time periods in the icing period, calculating a plurality of ice layer thicknesses of the icing of the surface of the wind driven generator blade in unit area in the plurality of first time periods; and accumulating the thicknesses of the plurality of ice layers to obtain the thickness of the ice coating.
Some embodiments, wherein the ice layer thickness is calculated for each first time period according to the following formula:
where δ represents the ice layer thickness, t s1 Representing a first time period duration; ρ is a unit of a gradient i Represents the density of ice; c. C W Representing the water vapor collection coefficient; α represents a correction coefficient; m v Representing the moisture content of the air during a first time period; r H1 Representing a relative humidity value of the air during a first time period; v. of 0 Representing the blade surface air flow rate during a first time period,
v w1 representing real-time wind speed, v, over a first period of time c1 Representing the linear speed of blade rotation during a first time period; TIS1 represents solar irradiance during a first time period;
representing solar irradiance influence coefficient during ice coating; s m Representing the windward area of the blade; rho W Indicating that the water vapor reaches the absolute humidity of the saturated gas at the same temperature; wherein R is H1 、v w1 TIS1 belongs to a first meteorological parameter, v c1 、S m Belonging to blade parameters, ρ i 、α、ρ W 、c W 、
Belonging to the first experimental coefficient.
In some embodiments, the second calculating module 502 is specifically configured to: calculating the instantaneous heat absorption capacity of the blade surface of the wind driven generator in unit area in the ice melting process of the second time period according to the second meteorological parameters, the blade parameters and the second experimental coefficient of the second time period which has passed in the ice melting period; and accumulating the instantaneous heat absorption to obtain a second heat.
Some embodiments, wherein the instantaneous heat absorption for each second time period is calculated according to the following equation:
wherein Q is Suction device Represents the instantaneous heat absorption; β represents a correction coefficient; r H2 Representing the air-to-humidity value during the second time period; v. of 1 Representing a blade surface air flow rate for a second time period; TIS2 represents solar irradiance during the second time period; s m Representing the windward area of the blade; t is t s2 Representing a second time period duration; h represents the heat exchange coefficient of the ice surface; Δ t represents the difference between the air temperature and the ice temperature over the second time period;
represents the solar irradiance influence coefficient during ice melting, wherein R H2 、v 1 TIS2, deltat belong to a second meteorological parameter, S m Belongs to the parameters of blade, beta, h,
Belonging to the second experimental coefficient.
In the embodiment, the first heat released in the icing process of the blade surface of the wind driven generator in unit area is calculated according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, the second heat absorbed in the deicing process of the blade surface of the wind driven generator in unit area is calculated according to the second meteorological parameter, the blade parameter and the second experimental coefficient in the deicing period from the starting time of the deicing period, the ending time of the deicing period is determined under the condition that the second heat is equal to the first heat, the heat release amount when the blade is iced and the heat absorption amount when the blade is iced can be calculated according to the meteorological parameters, the blade parameters and the experimental coefficients collected in real time, and the ending time of the deicing is deduced based on the thermodynamic equilibrium principle, so that operation and maintenance personnel can start the generator generation in time according to the ending time, and the safe and economic operation of a wind power plant unit and a power grid is facilitated to be ensured.
The present disclosure also provides a computer device, a readable storage medium, and a computer program product according to embodiments of the present disclosure.
In order to implement the foregoing embodiments, the present disclosure further provides a computer program product, which when executed by an instruction processor in the computer program product, performs the method for calculating the blade ice melting time of the wind turbine generator as set forth in the foregoing embodiments of the present disclosure.
FIG. 6 illustrates a block diagram of an exemplary computer device suitable for use in implementing embodiments of the present disclosure. The computer device 12 shown in FIG. 6 is only one example and should not impose any limitations on the functionality or scope of use of embodiments of the present disclosure.
As shown in FIG. 6, computer device 12 is in the form of a general purpose computing device. The components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that couples various system components including the system memory 28 and the processing unit 16.
Computer device 12 typically includes a variety of computer system readable media. Such media may be any available media that is accessible by computer device 12 and includes both volatile and nonvolatile media, removable and non-removable media.
Although not shown in FIG. 6, a disk drive for reading from and writing to a removable, nonvolatile magnetic disk (e.g., a "floppy disk") and an optical disk drive for reading from or writing to a removable, nonvolatile optical disk (e.g., a Compact disk Read Only Memory (CD-ROM), a Digital versatile disk Read Only Memory (DVD-ROM), or other optical media) may be provided. In these cases, each drive may be connected to bus 18 by one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to carry out the functions of embodiments of the disclosure.
A program/utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28, such program modules 42 including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or some combination of which may comprise an implementation of a network environment. Program modules 42 generally carry out the functions and/or methodologies of the embodiments described in this disclosure.
Computer device 12 may also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), with one or more devices that enable a user to interact with computer device 12, and/or with any devices (e.g., network card, modem, etc.) that enable computer device 12 to communicate with one or more other computing devices. Such communication may be through an input/output (I/O) interface 22. Moreover, computer device 12 may also communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), and/or a public Network such as the Internet) via Network adapter 20. As shown, the network adapter 20 communicates with the other modules of the computer device 12 over the bus 18. It should be understood that although not shown in the figures, other hardware and/or software modules may be used in conjunction with computer device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, among others.
The processing unit 16 executes various functional applications by running programs stored in the system memory 28, for example, implementing the wind turbine blade melting time calculation method mentioned in the foregoing embodiments.
Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This disclosure is intended to cover any variations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice in the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.
It will be understood that the present disclosure is not limited to the precise arrangements described above and shown in the drawings and that various modifications and changes may be made without departing from the scope thereof. The scope of the present disclosure is limited only by the appended claims.
It should be noted that, in the description of the present disclosure, the terms "first", "second", etc. are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. In addition, in the description of the present disclosure, the meaning of "a plurality" is two or more unless otherwise specified.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and the scope of the preferred embodiments of the present disclosure includes other implementations in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present disclosure.
It should be understood that portions of the present disclosure may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, any one or combination of the following technologies, which are well known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.
It will be understood by those skilled in the art that all or part of the steps carried by the method for implementing the above embodiments may be implemented by hardware related to instructions of a program, which may be stored in a computer readable storage medium, and when the program is executed, the program includes one or a combination of the steps of the method embodiments.
In addition, functional units in the embodiments of the present disclosure may be integrated into one processing module, or each unit may exist alone physically, or two or more units are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode. The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may also be stored in a computer readable storage medium.
The storage medium mentioned above may be a read-only memory, a magnetic or optical disk, etc.
In the description of the present specification, reference to the description of "one embodiment," "some embodiments," "an example," "a specific example," or "some examples" or the like means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present disclosure. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present disclosure have been shown and described above, it will be understood that the above embodiments are exemplary and not to be construed as limiting the present disclosure, and that changes, modifications, substitutions and alterations may be made to the above embodiments by those of ordinary skill in the art within the scope of the present disclosure.
Claims (10)
1. A method for calculating the blade ice melting time of a wind driven generator is characterized by comprising the following steps:
calculating a first heat released in the process of icing the surface of the wind driven generator blade in unit area according to a first meteorological parameter, a blade parameter and a first experiment coefficient in an icing period;
calculating a second heat absorbed in the unit area of the blade surface of the wind driven generator in the ice melting process according to a second meteorological parameter, the blade parameter and a second experimental coefficient in the ice melting period from the starting time of the ice melting period; and
determining an end time of the ice-melt period if the second amount of heat is equal to the first amount of heat.
2. The method according to claim 1, wherein calculating the first heat released during the icing process on the surface of the wind turbine blade per unit area according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing cycle comprises:
calculating the icing mass of the icing on the surface of the wind driven generator blade in the unit area according to the first meteorological parameter, the blade parameter and the first experiment coefficient; and
and calculating the first heat according to the ice coating quality, the specific heat capacity of the ice and the temperature variation of the water vapor condensed into the ice from 0 ℃.
3. The method according to claim 2, wherein calculating the icing mass of the unit area of the blade surface icing of the wind turbine according to the first meteorological parameter, the blade parameter and the first experimental coefficient comprises:
calculating the icing thickness of the icing on the surface of the wind driven generator blade in the unit area according to the first meteorological parameter, the blade parameter and the first experiment coefficient; and
and calculating the icing mass according to the icing thickness, the unit area and the ice density.
4. The method according to claim 3, wherein calculating the thickness of the ice coating on the surface of the wind turbine blade per unit area according to the first meteorological parameter, the blade parameter and the first experimental coefficient comprises:
according to the first meteorological parameters, the blade parameters and the first experiment coefficients of a plurality of first time periods in the icing period, calculating a plurality of ice layer thicknesses of the icing on the surface of the wind driven generator blade in the unit area in the plurality of first time periods; and
and accumulating the thicknesses of the plurality of ice layers to obtain the thickness of the ice coating.
5. The method of claim 4, wherein the ice layer thickness for each first time period is calculated according to the following equation:
wherein δ represents the ice layer thickness, t s1 Representing the first time period duration; rho i Represents the density of ice; c. C W Representing the water vapor collection coefficient; α represents a correction coefficient; m v Representing the water vapor content of the air during the first time period; r H1 Is shown inThe relative humidity of the air in the first time period; v. of 0 Representing a blade surface air flow rate during said first period of time,
v w1 representing real-time wind speed, v, over said first period of time c1 Representing a linear blade rotation speed during the first time period; TIS1 represents solar irradiance over the first time period;
representing solar irradiance influence coefficient during ice coating; s m Representing the windward area of the blade; rho W Indicating that the water vapor reaches the absolute humidity of the saturated gas at the same temperature; wherein R is H1 、v w1 TIS1 belongs to said first meteorological parameter, v c1 、S m Belonging to said blade parameter, p i 、α、ρ W 、c W 、
Belonging to said first experimental coefficient.
6. The method of claim 1, wherein calculating a second amount of heat that has been absorbed during the process of melting ice on the surface of the blade of the wind turbine per unit area according to a second meteorological parameter, the blade parameter, and a second experimental coefficient during the ice melting cycle comprises:
calculating the instantaneous heat absorption capacity of the blade surface of the wind driven generator in the unit area in the ice melting process of the second time period according to the second meteorological parameter, the blade parameter and the second experimental coefficient of the second time period which has passed in the ice melting period; and
and accumulating the instantaneous heat absorption quantity to obtain the second heat quantity.
7. The method of claim 6, wherein the instantaneous heat absorption for each second time period is calculated according to the following equation:
wherein Q is Instantaneous suction Representing the instantaneous heat absorption; β represents a correction coefficient; r is H2 Representing the air-to-humidity value during the second time period; v. of 1 Representing a blade surface air flow rate during the second period of time; TIS2 represents solar irradiance during the second time period; s. the m Representing the windward area of the blade; t is t s2 Representing the second time period duration; h represents the heat exchange coefficient of the ice surface; Δ t represents the difference between the air temperature and the ice temperature over the second time period;
representing the solar irradiance influence coefficient during ice melting, wherein R H2 、v 1 TIS2, deltat belong to the second meteorological parameter, S m Belongs to the parameters of blade, beta, h,
Belonging to said second experimental coefficient.
8. The method according to claim 1, wherein before calculating the first heat released in the process of icing the surface of the wind turbine blade per unit area according to the first meteorological parameter, the blade parameter and the first experimental coefficient in the icing period, the method further comprises the following steps:
acquiring meteorological forecast data of a preset time period; and
and judging whether the preset time period meets the icing condition or not according to the meteorological forecasting data.
9. A device for calculating the blade ice melting time of a wind driven generator is characterized by comprising:
the first calculating module is used for calculating first heat released in the icing process of the surface of the wind driven generator blade in unit area according to a first meteorological parameter, a blade parameter and a first experiment coefficient in the icing period;
the second calculation module is used for calculating second heat absorbed in the process of deicing the surface of the blade of the wind driven generator in the unit area according to a second meteorological parameter, the blade parameter and a second experimental coefficient in the deicing cycle from the starting time of the deicing cycle; and
a determining module for determining an end time of the ice-melt cycle if the second amount of heat is equal to the first amount of heat.
10. A non-transitory computer readable storage medium having stored thereon computer instructions for causing the computer to perform the method of any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211226474.9A CN115496391A (en) | 2022-10-09 | 2022-10-09 | Method and device for calculating blade ice melting time of wind driven generator and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211226474.9A CN115496391A (en) | 2022-10-09 | 2022-10-09 | Method and device for calculating blade ice melting time of wind driven generator and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115496391A true CN115496391A (en) | 2022-12-20 |
Family
ID=84472888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211226474.9A Pending CN115496391A (en) | 2022-10-09 | 2022-10-09 | Method and device for calculating blade ice melting time of wind driven generator and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115496391A (en) |
-
2022
- 2022-10-09 CN CN202211226474.9A patent/CN115496391A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gao et al. | A field study of ice accretion and its effects on the power production of utility-scale wind turbines | |
| KR101581102B1 (en) | Method for operating a wind energy plant in icing conditions | |
| US20120226485A1 (en) | Methods for predicting the formation of wind turbine blade ice | |
| Skrimpas et al. | Detection of icing on wind turbine blades by means of vibration and power curve analysis | |
| US9822762B2 (en) | System and method for operating a wind turbine | |
| EP3139038A1 (en) | A method for estimating the surface condition of a rotating blade | |
| Dierer et al. | Wind turbines in icing conditions: performance and prediction | |
| CN102418658A (en) | Method of and device for determining a mass condition of a rotor of a wind turbine, and method of operating a wind turbine | |
| CN114439706A (en) | Method for predicting icing state of fan blade of wind power plant | |
| WO2012149984A1 (en) | System and method for operating a wind turbine using an adaptive speed reference | |
| CN116662748A (en) | Prediction method and prediction device for fan icing thickness, storage medium and processor | |
| Seifert et al. | Aerodynamics of iced airfoils and their influence on loads and power production | |
| CN108062722B (en) | Mechanical power calculation method of mountain wind power plant model fan based on wind speed variation coefficient | |
| CN112682276A (en) | Fan blade icing state prediction method and device, medium and electronic equipment | |
| CN115585107A (en) | Method and system for predicting thickness of frozen fog ice coating of blade of wind generating set | |
| CN105138845B (en) | The method for obtaining wind-driven generator air speed value | |
| Botta et al. | Ice accretion at acqua spruzza and its effects on wind turbine operation and loss of energy production | |
| CN116702957A (en) | New energy power prediction method, equipment and storage medium for extreme weather | |
| Homola et al. | Modelling of ice induced power losses and comparison with observations | |
| CN115496391A (en) | Method and device for calculating blade ice melting time of wind driven generator and storage medium | |
| TW201934870A (en) | Wind power generation device and wind power generation system | |
| CN116760009A (en) | Wind power icing shutdown off-grid risk prediction method, wind power icing shutdown off-grid risk prediction device and processor | |
| CN116362382A (en) | Short-term power prediction method and system based on icing state of wind power plant | |
| Haaland | Estimating production loss due to icing on wind turbines | |
| Battisti | Icing on Wind Turbines |
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 |