Fan blade icing type detection method and system
By obtaining the ice covering weight on the rotating metal conductor, establishing an ice covering rate model and a thermal equilibrium model, the problem of difficulty in monitoring the type of ice covering of fan blades is solved, and data support for the anti-ice disaster reduction work of wind farms is achieved.
Patent Information
- Application Number
- CN202411519298.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-29
- Publication Date
- 2025-02-07
AI Technical Summary
Monitoring of the type of ice covering of fan blades is difficult, mainly because parameters such as liquid water content in the air and median diameter of water droplets are difficult to obtain under ice covering conditions.
By obtaining the time change curve of the ice-covered weight on the rotating metal conductor, establishing the conductor ice-covered rate model, inversely obtaining the ice-covered environmental parameters, and then building a thermal equilibrium model to solve the water droplet freezing coefficient on the surface of the conductor, judging the ice-covered type and calculating the ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-covered ice-cover thickness.
It has realized the real-time judgment of the type of ice covering of the fan blades under harsh conditions of ice covering, avoiding the impact of the rapid growth period of ice covering, and providing data support for the anti-ice disaster reduction work of wind farms.
Smart Images
Figure CN119397357A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of ice coating detection, and in particular to a method and system for detecting ice coating types on fan blades. Background Art
[0002] Icing on wind turbine blades has become a long-term safety issue threatening the operation of wind turbines. The hazards of wind turbine icing are closely related to the type of icing. In previous studies, icing types are generally distinguished according to the different icing densities, and are divided into rime, rime, mixed rime, and hoarfrost. When the ambient temperature is high (below zero degrees Celsius), the wind speed is low, and the liquid water content in the air is sufficient, the icing type is generally rime. Under this icing type, a transparent or milky white smooth ice layer is easily formed on the surface of the wind turbine blades. Due to the presence of water film, the density of rime is relatively large (greater than 0.8g / cm3), which can easily cause the wind turbine blades to overload and break. When the ambient temperature drops further, the icing type on the surface of the wind turbine blades will gradually change to mixed rime, hard rime, and soft rime according to the different freezing speeds of water droplets. When the freezing speed of water droplets increases, the water film flows slowly, the ice ridges decrease, and the icing type is often mixed rime; when the water film flow speed further decreases, the flow stagnates or even the water film disappears, the fan blade type changes to hard rime or soft rime. Compared with rain rime, the density of rime is smaller and it is not easy to form a serious overload, but the mixed rime of alternating rain rime and rime is also the focus of icing monitoring.
[0003] As a key equipment for renewable energy power generation, the operating status of wind turbines directly affects the safety and stability of the power grid. However, with the vigorous development of wind power generation in recent years, the harm of icing to wind turbine operation has become increasingly prominent. Accurately judging the development type of icing on wind turbine blades can assist operation and maintenance personnel in accurately carrying out wind turbine shutdown, anti-icing, and de-icing work. The icing process of wind turbine blades is complex. Various environmental parameters including wind speed, ambient temperature, liquid water content in the air, median diameter of water droplets, etc. will affect the type and rate of icing development. However, parameters such as liquid water content in the air and median diameter of water droplets under icing conditions are extremely difficult to obtain, which makes it difficult to monitor the type of icing on wind turbine blades.
[0004] Therefore, it is necessary to propose a high-efficiency and high-precision method for monitoring the type of wind turbine blade icing. This method can avoid the influence of severe icing conditions and judge the icing development type in real time during the rapid growth period of icing, thereby providing data support for wind farm anti-icing and disaster reduction work. Summary of the invention
[0005] In view of the above existing problems, the present invention is proposed.
[0006] Therefore, the present invention provides a method and system for detecting the type of icing on fan blades to solve the problem that parameters such as the liquid water content in the air and the median diameter of water droplets are extremely difficult to obtain under current icy conditions, which makes it difficult to monitor the type of icing on fan blades.
[0007] In order to solve the above technical problems, the present invention provides the following technical solutions:
[0008] In a first aspect, the present invention provides a method for detecting the type of ice coating on a wind turbine blade, comprising:
[0009] Acquiring parameters of an icing environment, wherein the icing environment is provided with a rotating metal conductor;
[0010] Obtaining the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor;
[0011] Converting the time variation curve into a conductor icing rate model to obtain liquid water parameters in an icing environment;
[0012] The water drop collision coefficient on the conductor surface is calculated based on the conductor ice coating rate model, and the conductor surface thermal balance model is constructed based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface;
[0013] The type of ice coating is determined based on the freezing coefficient, and the thickness of ice coating on the conductor surface is calculated.
[0014] As a preferred solution of the wind turbine blade icing type detection method of the present invention, the time variation curve is converted into a conductor icing rate model, and the conductor icing rate model constructs a function through wind speed, liquid water content in the air, median volume diameter of water droplets, ambient temperature and conductor length, which is expressed as:
[0015] dM i (Δt)=f(U,MVD,w,T,D i )(i=1,2)
[0016] Where U is the wind speed, MVD is the median volume diameter of the water droplet, w is the liquid water content, T is the ambient temperature, and D i is the conductor length, and i represents the number of conductors.
[0017] As a preferred solution of the method for detecting icing type of wind turbine blades described in the present invention, the liquid water parameters in the icing environment include the median volume diameter of water droplets and the liquid water content, which are obtained by reverse deduction of the conductor icing rate model.
[0018] As a preferred solution of the method for detecting the type of icing on wind turbine blades according to the present invention, the water drop collision coefficient on the conductor surface is calculated based on the conductor icing rate model, including:
[0019] The conductor icing rate model is converted into a function constructed by the conductor length, wind speed, liquid water content, water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface during the period; wherein the water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface are unknown;
[0020] The airflow velocity at different positions is calculated based on the time variation curve combined with the use of the circular conductor potential stream function;
[0021] The speed and trajectory of water droplets are obtained according to the drag force of airflows of different speeds on water droplets;
[0022] Based on the calculated parameters, the value of the water droplet collision coefficient on the rotating conductor surface can be obtained by reverse deduction.
[0023] As a preferred solution of the wind turbine blade icing type detection method of the present invention, a conductor surface thermal balance model is constructed based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface, including:
[0024] Solve the thermal equilibrium model including the freezing coefficient of water droplets on the rotating conductor surface;
[0025] The water drop freezing coefficient is expressed as:
[0026]
[0027] Among them, α 2 is the freezing coefficient of the conductor surface, M 冻结 M is the mass of ice frozen on the conductor surface per unit time, 碰撞捕获 is the mass of water droplets that collide with the conductor and are captured per unit time;
[0028] The heat balance model is expressed as:
[0029] q c +q e +q r +q l +q s =q f +q v +q k +q d
[0030] Among them, q s is the heat carried away by the unfrozen water droplets; q l The heat consumed to heat the water droplets to freezing temperature; q r is the long-wave radiation heat loss; q e is the heat loss due to evaporation and sublimation; q c represents the convective heat transfer loss with air;
[0031] q v represents the friction heating of the round conductor by the airflow; q d The heat released when ice at 0°C freezes to the steady-state temperature of the circular conductor surface; q k represents the collision kinetic energy of supercooled water; q f It is the latent heat released when water freezes.
[0032] As a preferred solution of the method for detecting ice coating type of wind turbine blades according to the present invention, the method of determining the ice coating type according to the freezing coefficient includes:
[0033] When the freezing coefficient is less than the first threshold, it is judged as wet growth ice;
[0034] When the freezing coefficient is greater than or equal to the first threshold, it is determined to be dry growth icing.
[0035] As a preferred solution of the method for detecting the type of ice coating on wind turbine blades according to the present invention, it further includes:
[0036] The conductor icing rate model and solution process are constructed repeatedly to continuously output the icing development type within the current time step;
[0037] The conductor diameter parameter is updated by the continuously obtained ice thickness.
[0038] In a second aspect, the present invention provides a wind turbine blade icing type detection system, comprising:
[0039] A first acquisition module is used to acquire parameters of an icing environment, wherein the icing environment is provided with a rotating metal conductor;
[0040] A second acquisition module is used to acquire the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor;
[0041] A first model solving module is used to convert the time variation curve into a conductor icing rate model to obtain liquid water parameters in an icing environment;
[0042] The second model solving module is used to calculate the water drop collision coefficient on the conductor surface based on the conductor ice coating rate model, and to construct a conductor surface thermal balance model based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface;
[0043] The judgment module is used to judge the type of ice coating according to the freezing coefficient and calculate the ice coating thickness on the conductor surface.
[0044] In a third aspect, the present invention provides a computing device, comprising:
[0045] Memory and processor;
[0046] The memory is used to store computer executable instructions, and the processor is used to execute the computer executable instructions. When the computer executable instructions are executed by the processor, the steps of the method for detecting the type of ice coating on wind turbine blades are implemented.
[0047] In a fourth aspect, the present invention provides a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the steps of the wind turbine blade icing type detection method.
[0048] Compared with the prior art, the beneficial effects of the present invention are as follows: by acquiring the time-varying data of the ice weight of two rotating conductors, an equation for the ice rate on the conductor surface is established, the ice environment parameters are obtained by inverse calculation, and then a heat balance equation is established to solve the conductor freezing coefficient, and the ice type of the wind turbine blades is judged accordingly, avoiding the influence of severe icing conditions, and the ice development type can be judged in real time during the rapid growth period of icing, thereby providing data support for the anti-icing and disaster reduction work of wind farms and bringing convenience to the technical reference for anti-icing and disaster reduction work. BRIEF DESCRIPTION OF THE DRAWINGS
[0049] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.
[0050] Figure 1 A schematic diagram of the overall process logic of a method for detecting ice coating types on wind turbine blades according to an embodiment of the present invention;
[0051] Figure 2 A specific process logic diagram of a method for detecting ice coating type of wind turbine blades according to an embodiment of the present invention;
[0052] Figure 3 A physical diagram of two rotating circular conductors of different diameters involved in a method for detecting ice coating types on wind turbine blades according to an embodiment of the present invention;
[0053] Figure 4 A schematic diagram of the collision coefficient value of a rotating conductor water droplet involved in a method for detecting ice coating types on wind turbine blades according to an embodiment of the present invention;
[0054] Figure 5 This is a physical diagram of the icing result states of a rotating conductor under different icing types in a method for detecting icing types of wind turbine blades according to an embodiment of the present invention. DETAILED DESCRIPTION
[0055] In order to make the above-mentioned purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work should fall within the scope of protection of the present invention.
[0056] Example 1, reference Figure 1-Figure 5 Tables 1 and 2 are an embodiment of the present invention, which provides a method for detecting ice coating types of wind turbine blades, including:
[0057] S100: Acquire parameters of an icing environment, where a rotating metal conductor is provided in the icing environment;
[0058] Specifically, the environmental parameters include wind speed and ambient temperature.
[0059] Specifically, two rotating metal conductors are provided in the ice-covered environment, and each metal conductor has the same length.
[0060] Specifically, the metal conductor can be made of iron, steel or other materials, preferably an aluminum conductor.
[0061] Exemplarily, as shown in Table 1, two rotating conductors are evenly arranged in an ice-covered environment.
[0062] Preferably, the length of the metal conductor is 30 cm-50 cm.
[0063] Preferably, Figure 3 As shown, two metal conductors can be evenly arranged in an ice-covered environment and rotated at a speed of 1 revolution per minute.
[0064] Table 1 Basic parameters of two rotating conductors
[0065]
[0066] Specifically, the wind speed and the ambient temperature are measured by arranging a meteorological sensor in parallel with the rotating conductor.
[0067] S200: Obtaining the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor;
[0068] Specifically, the weight of ice covering the two rotating conductors can be collected in real time by using a high-precision tension sensor.
[0069] It should be noted that the obtained time variation curve of the conductor ice coating weight is the input data of the conductor ice coating rate model function, which is used to solve the input of the model equation.
[0070] S300: converting the time variation curve into a conductor icing rate model to obtain liquid water parameters in the icing environment;
[0071] In the embodiment of the present application, the time variation curve is converted into a conductor icing rate model. The conductor icing rate model constructs a function through wind speed, liquid water content in the air, median volume diameter of water droplets, ambient temperature and conductor length, which is expressed as:
[0072] dM i (Δt)=f(U,MVD,w,T,D i )(i=1,2)
[0073] Where U is the wind speed, MVD is the median volume diameter of the water droplet, w is the liquid water content, T is the ambient temperature, and D i is the conductor length, and i represents the number of conductors.
[0074] Specifically, the conductor icing rate dM containing two unknown icing environment parameters, the liquid water content w in the air and the median volume diameter MVD of the water droplets, is i (Δt) equation, expressed as:
[0075]
[0076] In the embodiment of the present application, the liquid water parameters in the icing environment include the median volume diameter of the water droplets and the liquid water content, which are obtained by reversely solving the conductor icing rate model.
[0077] S400: calculating a water drop collision coefficient on the conductor surface based on a conductor ice coating rate model, and constructing a conductor surface thermal balance model based on the obtained collision coefficient to obtain a water drop freezing coefficient on the conductor surface;
[0078] In the embodiment of the present application, the water droplet collision coefficient on the conductor surface is calculated based on the conductor ice coating rate model, including:
[0079] The conductor icing rate model is converted into a function constructed by the conductor length, wind speed, liquid water content, water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface during the period; wherein the water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface are unknown;
[0080] The airflow velocity at different positions is calculated based on the time variation curve combined with the use of the circular conductor potential stream function;
[0081] The speed and trajectory of water droplets are obtained according to the drag force of airflows of different speeds on water droplets;
[0082] Based on the calculated parameters, the value of the water droplet collision coefficient on the rotating conductor surface can be obtained by reverse deduction.
[0083] Specifically, based on steps S300 and S400, the specific solution process of the f function is as follows:
[0084] dM i After (Δt) transformation, it can be expressed as:
[0085]
[0086] Among them, K and are two dimensionless parameters. K is the Stokes number of the inertia of the moving water droplets; Re is the Reynolds number of the relative motion of the water droplets in the airflow; Defined as the ratio of Re to K.
[0087]
[0088] Among them, R d is the radius of the water drop, m; ρ d and ρ a are the water drop and air densities, kg / m 3 ; u and v are the velocity vectors of gas and water droplet, respectively, m / s; U is the wind speed, m / s; R is the conductor radius, m; μ is the kinematic viscosity of air, m 2 / s.
[0089] It should be noted that u and v are unknown in formula (5). However, according to the uniform ice-covered morphology of the rotating conductor, it can maintain a circular outer contour with or without ice-covered conditions. The potential flow function of the circular conductor, such as formula (6), can be used to calculate the velocity u of the airflow at different positions, and then the velocity v and trajectory of the water droplet can be obtained according to the drag force of the airflow at different speeds on the water droplet.
[0090] Therefore, according to formulas (3) to (7), the water droplet collision coefficient α on the rotating conductor surface can be obtained: 1 The value of Figure 4 shown.
[0091] In the embodiment of the present application, a conductor surface thermal balance model is constructed based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface, including:
[0092] Solve the thermal equilibrium model including the freezing coefficient of water droplets on the rotating conductor surface;
[0093] The water drop freezing coefficient is expressed as:
[0094]
[0095] Among them, α 2 is the freezing coefficient of the conductor surface, M 冻结 M is the mass of ice frozen on the conductor surface per unit time, 碰撞捕获is the mass of water droplets that collide with the conductor and are captured per unit time;
[0096] The heat balance model is expressed as:
[0097] q c +q e +q r +q l +q s =q f +q v +q k +q d (9)
[0098] Among them, q s is the heat carried away by the unfrozen water droplets; q l The heat consumed to heat the water droplets to freezing temperature; q r is the long-wave radiation heat loss; q e is the heat loss due to evaporation and sublimation; q c represents the convective heat transfer loss with air;
[0099] q v represents the friction heating of the round conductor by the airflow; q d The heat released when ice at 0°C freezes to the steady-state temperature of the circular conductor surface; q k represents the collision kinetic energy of supercooled water; q f It is the latent heat released when water freezes.
[0100] Solving for α 2 The expression is:
[0101]
[0102] Among them, T s 、T s are the freezing point temperature and the ambient temperature, respectively, °C; h is the convective heat transfer coefficient; S b , D are the surface area and diameter of the rotating conductor respectively; x is the evaporation coefficient; emissivity ε=0.95; e(T) is the saturated water vapor pressure on the surface of the circular conductor at temperature T; δ R is the Stefan-Boltzmann constant; L f is the latent heat of melting of ice, J / kg; C w is the specific heat capacity of water, J / (kg.℃).
[0103] It should be noted that q in formula (9) f =α 2 ·M 碰撞捕获 ·C w (T s -T), after replacement, the other items can be solved by referring to ISO12494 standard. By solving, we can get α in formula (10)2 expression.
[0104] It can also be seen from this that in α 1 , α 2 Under known conditions, the f function of equation (1) can be solved, and equation (1) becomes a system of equations including two unknowns w and MVD, which can also be solved.
[0105] S500: Determine the type of ice coating according to the freezing coefficient, and calculate the thickness of ice coating on the conductor surface.
[0106] In the embodiment of the present application, judging the ice type according to the freezing coefficient includes:
[0107] When the freezing coefficient is less than the first threshold, it is judged as wet growth ice;
[0108] When the freezing coefficient is greater than or equal to the first threshold, it is determined to be dry growth icing.
[0109] For example, the step time is set to 1h, and then according to the ice mass of the two rotating conductors measured and the wind speed and ambient temperature measured by the meteorological sensor, the results of some ice types are shown in Table 2. Figure 5 The ice cover types shown are dry growth on the left and wet growth on the right.
[0110] Table 2 Results of icing types
[0111]
[0112] In the embodiment of the present application, it also includes:
[0113] S600: repeatedly construct the conductor icing rate model and solution process, and continuously output the icing development type within the current time step; update the conductor diameter parameter by the continuously obtained icing thickness, and judge the development of the icing type in real time.
[0114] In summary, icing on wind turbine blades seriously endangers the safe and stable operation of wind turbines, and different types of icing have different hazards. This method uses two rotating conductors to monitor the icing type of wind turbine blades in real time, providing a technical reference for wind turbine anti-icing and disaster reduction work.
[0115] The above is a schematic scheme of a method for detecting ice coating type of a fan blade in this embodiment. It should be noted that the technical scheme of the system for detecting ice coating type of a fan blade is the same as the technical scheme of the above method for detecting ice coating type of a fan blade. For details not described in detail in the technical scheme of the system for detecting ice coating type of a fan blade in this embodiment, please refer to the description of the technical scheme of the above method for detecting ice coating type of a fan blade.
[0116] In this embodiment, the wind turbine blade icing type detection system includes:
[0117] A first acquisition module is used to acquire parameters of an icing environment, where a rotating metal conductor is provided in the icing environment;
[0118] The second acquisition module is used to acquire the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor;
[0119] The first model solving module is used to convert the time variation curve into a conductor icing rate model to obtain liquid water parameters in an icing environment;
[0120] The second model solving module is used to calculate the water drop collision coefficient on the conductor surface based on the conductor ice coating rate model, and to construct a conductor surface thermal balance model based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface;
[0121] The judgment module is used to judge the type of ice coating according to the freezing coefficient and calculate the ice coating thickness on the conductor surface.
[0122] This embodiment further provides a computing device, which is applicable to the case of detecting the type of ice coating on wind turbine blades, and includes:
[0123] Memory and processor; the memory is used to store computer executable instructions, and the processor is used to execute computer executable instructions to implement the method for detecting the type of ice coating on wind turbine blades as proposed in the above embodiment.
[0124] This embodiment further provides a storage medium on which a computer program is stored. When the program is executed by a processor, the method for detecting the type of ice coating on wind turbine blades as proposed in the above embodiment is implemented.
[0125] The storage medium proposed in this embodiment and the method for detecting the type of icing on wind turbine blades proposed in the above embodiment belong to the same inventive concept. The technical details not fully described in this embodiment can be referred to the above embodiment, and this embodiment has the same beneficial effects as the above embodiment.
[0126] Through the above description of the implementation methods, the technicians in the relevant field can clearly understand that the present invention can be implemented by means of software and necessary general hardware, and of course can also be implemented by hardware. Based on such an understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as a computer floppy disk, read-only memory (ReadOnly, Memory, ROM), random access memory (Random Access Memory, RAM), flash memory (FLASH), hard disk or optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, server, or network device, etc.) to perform the methods of various embodiments of the present invention.
[0127] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit the present invention. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.
Claims
1. A method for detecting the type of ice coating on a wind turbine blade, characterized in that: include: Acquiring parameters of an icing environment, wherein the icing environment is provided with a rotating metal conductor; Obtaining the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor; Converting the time variation curve into a conductor icing rate model to obtain liquid water parameters in an icing environment; The water drop collision coefficient on the conductor surface is calculated based on the conductor ice coating rate model, and the conductor surface thermal balance model is constructed based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface; The type of ice coating is determined based on the freezing coefficient, and the thickness of ice coating on the conductor surface is calculated.
2. The method for detecting the type of ice coating on a wind turbine blade according to claim 1, characterized in that: The time variation curve is converted into a conductor icing rate model, and the conductor icing rate model constructs a function through wind speed, liquid water content in the air, median volume diameter of water droplets, ambient temperature and conductor length, which is expressed as: dM i (Δt)=f(U,MVD,w,T,D i )(i=1,2) Where U is the wind speed, MVD is the median volume diameter of the water droplet, w is the liquid water content, T is the ambient temperature, and D i is the conductor length, and i represents the number of conductors.
3. The method for detecting the type of ice coating on wind turbine blades according to claim 2, characterized in that: The liquid water parameters in the icing environment include the median volume diameter of the water droplets and the liquid water content, which are obtained by reversely solving the conductor icing rate model.
4. The method for detecting the type of ice coating on a wind turbine blade according to claim 3, characterized in that: The water droplet collision coefficient on the conductor surface is calculated based on the conductor icing rate model, including: The conductor icing rate model is converted into a function constructed by the conductor length, wind speed, liquid water content, water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface during the period; wherein the water drop collision coefficient on the rotating conductor surface and freezing coefficient on the circular conductor surface are unknown; The airflow velocity at different positions is calculated based on the time variation curve combined with the use of the circular conductor potential stream function; The speed and trajectory of water droplets are obtained according to the drag force of airflows of different speeds on water droplets; Based on the calculated parameters, the value of the water droplet collision coefficient on the rotating conductor surface can be obtained by reverse deduction.
5. The method for detecting the type of ice coating on wind turbine blades according to claim 4, characterized in that: Based on the obtained collision coefficient, a thermal balance model of the conductor surface is constructed to obtain the water drop freezing coefficient on the conductor surface, including: Solve the thermal equilibrium model including the freezing coefficient of water droplets on the rotating conductor surface; The water drop freezing coefficient is expressed as: Among them, α2 is the freezing coefficient of the conductor surface, M 冻结 M is the mass of ice frozen on the conductor surface per unit time, 碰撞捕获 is the mass of water droplets that collide with the conductor and are captured per unit time; The heat balance model is expressed as: q c +q e +q r +q l +q s =q f +q v +q k +q d Among them, q s is the heat carried away by the unfrozen water droplets; q l The heat consumed to heat the water droplets to freezing temperature; q r is the long-wave radiation heat loss; q e is the heat loss due to evaporation and sublimation; q c represents the convective heat loss with air; q v represents the friction heating of the round conductor by the airflow; q d The heat released when ice at 0°C freezes to the steady-state temperature of the circular conductor surface; q k represents the collision kinetic energy of supercooled water; q f It is the latent heat released when water freezes.
6. The method for detecting the type of ice coating on wind turbine blades according to claim 5, characterized in that: The step of judging the ice type according to the freezing coefficient includes: When the freezing coefficient is less than the first threshold, it is judged as wet growth ice; When the freezing coefficient is greater than or equal to the first threshold, it is determined to be dry growth icing.
7. The method for detecting the type of icing on wind turbine blades according to claim 5 or 6, characterized in that: Also includes: The conductor icing rate model and solution process are constructed repeatedly to continuously output the icing development type within the current time step; The conductor diameter parameter is updated by the continuously obtained ice thickness.
8. A system using the method for detecting ice coating type of wind turbine blades as claimed in any one of claims 1 to 7, characterized in that: include: A first acquisition module is used to acquire parameters of an icing environment, wherein the icing environment is provided with a rotating metal conductor; A second acquisition module is used to acquire the weight of ice on the metal conductor to obtain a time variation curve of the weight of ice on the conductor; A first model solving module is used to convert the time variation curve into a conductor icing rate model to obtain liquid water parameters in an icing environment; The second model solving module is used to calculate the water drop collision coefficient on the conductor surface based on the conductor ice coating rate model, and to construct a conductor surface thermal balance model based on the obtained collision coefficient to obtain the water drop freezing coefficient on the conductor surface; The judgment module is used to judge the type of ice coating according to the freezing coefficient and calculate the ice coating thickness on the conductor surface.
9. An electronic device, comprising: Memory and processor; The memory is used to store computer executable instructions, and the processor is used to execute the computer executable instructions. When the computer executable instructions are executed by the processor, the steps of the method for detecting the type of icing on wind turbine blades according to any one of claims 1 to 7 are implemented.
10. A computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the steps of the method for detecting ice coating type of wind turbine blades as claimed in any one of claims 1 to 7.