Experimental and calculation method of power loss of wind turbine blade icing

Abstract

The present application relates to the technical field of fan blade icing prediction, and discloses a fan blade icing and power loss experiment and calculation method to establish a quantitative relationship model of fan blade icing power change and realize fan blade icing power loss calculation. The method comprises the following steps: testing a fan model in an artificial climate chamber to obtain acquisition data under different working conditions, wherein the acquisition data comprises measured current and voltage values under different icing thicknesses and different wind speed conditions; then obtaining fan total output, wind energy utilization coefficient and power loss function under different icing thicknesses and different wind speed conditions based on the acquisition data, and further obtaining power loss corresponding to different wind speeds and different icing thicknesses.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of wind turbine blade icing prediction, and particularly relates to a wind turbine blade icing and power loss experiment and calculation method. BACKGROUND

[0002] With the rapid increase of wind power clean energy installation and the gradual increase of wind power installation scale, the influence of wind power grid connection on the safety and stability of power systems gradually highlights. The large-scale construction of wind farms has brought challenges to the safe and stable operation of power grids. In recent years, there have been operation accidents of large-scale wind power system chain trip at home and abroad.

[0003] At present, for the problem of wind turbine icing, monitoring and early warning technology research is mainly carried out, and the power loss caused by wind turbine blade icing is not considered, which needs to be improved. SUMMARY

[0004] The present application aims to disclose a wind turbine blade icing and power loss experiment and calculation method to establish a quantitative relationship model of wind turbine blade icing power change and realize wind turbine blade icing power loss calculation.

[0005] In order to achieve the above purpose, the method of the present application comprises:

[0006] (1) Prepare basic experimental conditions. According to the requirements, select an artificial climate laboratory with temperature regulation, precipitation regulation and wind speed regulation; prepare a small wind turbine model and its power generation equipment; prepare an ammeter, a data acquisition card, a multimeter, a thermometer, a barometer, a vernier caliper and other equipment.

[0007] (2) Place the small wind turbine model and its power generation equipment in the artificial climate chamber simulation environment, combine the three-phase delta circuit according to the wind turbine output curve, set the resistance between each two phases to R ohms, and the equivalent Y-type circuit resistance value is R / 3 ohms, and connect the data acquisition terminals of the data acquisition card to the two ends of the resistance on phase A.

[0008] (3) Fix the vertical rod on the wind turbine base, and mark different positions of a wind turbine blade at the same time for recording the ice thickness at different positions of the wind turbine blade.

[0009] (4) Turn on the artificial climate chamber refrigeration system, set the environmental temperature to a constant value, and when the environmental temperature decreases to the value, turn on the simulated wind system, set the wind speed to v0 m / s, close the wind turbine connection circuit, turn on the data acquisition card driver program, and start receiving the three-phase circuit voltage and current waveform.

[0010] (5) After running for t minutes, turn off the simulated wind system and the liquid outlet system, and put on protective equipment. Enter the artificial climate chamber and measure and record the ice thickness on the marked position on the fan blade (including the blade thickness) using a vernier caliper: save the voltage and current waveforms recorded on the data acquisition card.

[0011] (6) Without deicing, stand up the fan again, and turn on the data acquisition program and the artificial climate chamber simulated wind system. Set multiple wind speeds at intervals of Δv, and repeat steps (5)-(6);

[0012] (7) Extend the running time by different amounts, and repeat steps (5)-(7);

[0013] (8) After the experiment is completed, save the recorded data, disconnect all power supplies, and organize the experimental instruments and restore the test site.

[0014] (9) According to the current and voltage values measured by the data acquisition card under different ice thicknesses and different wind speeds, calculate the total output power of the fan:

[0015]

[0016] (10) Further, the wind energy utilization coefficient of the fan blade under different ice thicknesses and different wind speeds can be calculated:

[0017]

[0018] wherein, C p(v,I) is the wind energy utilization coefficient when the wind speed is v and the ice thickness is I, P w is the output power, r is the wind wheel sweep radius, v is the incoming wind speed, and π is the circular constant;

[0019] ρ is the air density, which is calculated as follows:

[0020]

[0021] wherein, p is the air pressure, J is the gas constant, and R is the gas constant 287 J / kg.K.

[0022] (11) Based on the power curve and the wind energy utilization coefficient C p(v,0) without ice, the power loss coefficient under different wind speeds v and ice thicknesses I is calculated:

[0023]

[0024] wherein, ΔC p(v,I) is the power loss coefficient when the wind speed is v and the ice thickness is I; C p(v,I) is the wind energy utilization coefficient when the wind speed is v and the ice thickness is I, and Cp(v,0) is the wind energy utilization coefficient when the wind speed is v and the icing thickness is 0.

[0025] (12) the power loss when the wind speed is v and the icing thickness is I is calculated as follows:

[0026] ΔP w(v,I) = ΔC p(v,I) × P w(v,0)

[0027] wherein P w(v,0) is the power output value of the fan when the fan is not iced and the wind speed is v, which can be obtained from the fan's out-of-field power curve.

[0028] The present application has the following beneficial effects:

[0029] The fan model is tested in an artificial climate chamber to obtain collected data under different working conditions, which includes the measured current and voltage values under different icing thicknesses and different wind speeds; then the fan total output, the wind energy utilization coefficient under different icing thicknesses and different wind speeds, and the power loss function under the fan blade conditions are obtained based on the collected data, and the power loss corresponding to different wind speeds and different icing thicknesses is obtained. It is simple and easy to implement, and the collected data is reliable. The quantitative relationship model of the power change of the fan blade icing can be conveniently and reliably established, and the power loss calculation of the fan blade icing is realized.

[0030] The present application will be further described in detail below with reference to the accompanying drawings. BRIEF DESCRIPTION OF DRAWINGS

[0031] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application, and are incorporated in and constitute a part of this application. The embodiments of these drawings illustrate the preferred embodiments of the present application and, together with their description, serve to explain the application. In the drawings:

[0032] Figure 1 is the wind speed-power curve diagram under different icing thicknesses calculated by the present application. DETAILED DESCRIPTION

[0033] The embodiments of the present application will be described in detail below with reference to the accompanying drawings, but the present application can be implemented in various different ways limited and covered by the claims.

[0034] Example 1

[0035] Step 1, prepare basic experimental conditions. According to the demand to choose with temperature regulation, precipitation regulation, wind speed regulation artificial climate laboratory; Prepare NF-200TS small fan model and its power generation equipment with rated power of 200w, rated voltage of 24v, starting wind speed of 2m / s; Prepare cl011mm ammeter, data acquisition card, FLUKE17B multimeter, IES-1310 thermometer, Pdca9000 barometer, YCDCA-100 vernier caliper, etc.

[0036] Step 2, place the small fan model and its power generation equipment in the artificial climate chamber simulation environment, according to the fan output curve, combine three-phase delta circuit, set the resistance between each two phase to 180 ohm, the equivalent Y type circuit resistance value is 60 ohm, and connect AI0 and AI1 data acquisition terminals of data acquisition card respectively at both ends of A phase resistance.

[0037] Step 3, fix the vertical rod on the fan base, and mark different positions of the fan blade at the same time, which is used to record the ice thickness at different positions of the fan blade.

[0038] Step 4, open the artificial climate chamber refrigeration system, set the environment temperature to-4℃; When the environment temperature reaches-4℃, open the simulation wind system, set the wind speed to 3m / s, close the fan connection circuit, open the data acquisition card driver program, and start receiving three-phase circuit voltage and current waveform.

[0039] Step 5, after running for 10 minutes, close the simulation wind system and liquid outlet system, wear protective equipment, enter the artificial climate chamber, measure and record the ice thickness at the marked position of the fan blade (including blade thickness) using vernier caliper: save the voltage and current waveform recorded on the data acquisition card.

[0040] Step 6, re-erect the fan, and open the data acquisition program and artificial climate chamber simulation wind system, set to 4m / s, 5m / s, 6m / s……10m / s, repeat steps (5)-(6).

[0041] Step 7, extend the running time for different time, repeat steps (5)-(7).

[0042] Step 8, after the experiment is completed, save the recorded data, disconnect all power supply, arrange the experimental instruments, and restore the test site.

[0043] Step 9, calculate the total output of the fan according to the current and voltage values measured by the data acquisition card under different ice thickness and different wind speed conditions.

[0044]

[0045] P = U * I wFor the total output of the fan, U 有效 For the effective value of the A-phase voltage where the data acquisition card is located, I 有效 For the effective value of the A-phase current where the data acquisition card is located.

[0046] Step 10, calculate the wind energy utilization coefficient under the condition of different icing thickness and different wind speed of the fan blade:

[0047]

[0048] Wherein, C p(v,I) is the wind energy utilization coefficient when the wind speed is v and the icing thickness is I, P w is the output power, r is the wind wheel sweep radius, v is the incoming wind speed, π is the circular constant; ρ is the air density, which is calculated as follows:

[0049]

[0050] Wherein, p is the air pressure, J is the gas constant, usually 287 J / kg·K.

[0051] Step 11, based on the power curve and its wind energy utilization coefficient C p(v,0) under the condition of no icing, the power loss coefficient under the condition of different wind speed v and icing thickness I is calculated:

[0052]

[0053] Wherein, ΔC p(v,I) is the power loss coefficient when the wind speed is v and the icing thickness is I; C p(v,I) is the wind energy utilization coefficient when the wind speed is v and the icing thickness is I, C p(v,0) is the wind energy utilization coefficient when the wind speed is v and the icing thickness is 0 (i.e. no icing).

[0054] Step 12, thus the power loss when the wind speed is v and the icing thickness is I can be calculated:

[0055] ΔP w(v,I) = ΔC p(v,I) × P w(v,0)

[0056] Wherein, P w(v,0) is the power output value of the fan when it is not iced and the wind speed is v, which can be obtained from the field power curve of the fan.

[0057] The wind speed-power curve under different icing thickness calculated is shown in Figure 1 .

[0058] The above merely provides the preferred embodiments of the present application, and is not used to limit the present application. For those skilled in the art, the present application can have various modifications and changes. Any modifications, equivalent replacements, improvements, etc. made within the principles and technical scope of the present application shall fall into the scope of the present application.

Claims

1. A method for testing and calculating the icing and power loss of a wind turbine blade, characterized in that, Comprise: (1) Prepare basic experimental conditions; (2) The small fan model and its power generation equipment are placed in the artificial climate chamber simulation environment, and the three-phase delta circuit is combined according to the fan output curve, the resistance between each two phases is set to R ohms, and the equivalent Y type circuit resistance value is R / 3 ohms, and the data acquisition terminals of the data acquisition card are connected to the resistance on the A phase; (3) The vertical rod is fixed on the fan base, and the different positions of the fan blade are marked for recording the ice thickness of the different positions of the fan blade; (4) Open the artificial climate chamber refrigeration system, set the environmental temperature to a constant value, and when the environmental temperature decreases to the value, open the simulation wind system, set the wind speed to v0 m / s, close the fan connection circuit, open the data acquisition card driver program, and start receiving the three-phase circuit voltage and current waveform; (5) After running for t minutes, close the simulation wind system and liquid outlet system, wear protective equipment, enter the artificial climate chamber, measure and record the ice thickness of the marked position of the fan blade including the blade thickness using the vernier caliper: save the voltage and current waveform recorded on the data acquisition card; (6) Without deicing, the fan is erected again, and the data acquisition program and the artificial climate chamber simulation wind system are opened, set to multiple wind speeds with interval Δv, repeat steps (5)-(6); (7) The running time is extended for different lengths of time, and steps (5)-(7) are repeated; (8) After the experiment is completed, save the recorded data, disconnect all power supplies, arrange the experimental instruments, and restore the test site; (9) According to the current and voltage values measured by the data acquisition card under different ice thicknesses and different wind speeds, the total output of the fan is calculated: (10) Calculate the wind energy utilization coefficient of the fan blade under different ice thicknesses and different wind speeds: wherein C p(v,I) is the wind energy utilization coefficient for a wind speed of v and an ice thickness of I, P w is the output power, r is the wind wheel sweep radius, v is the incoming wind speed, π is the circle constant; and ρ is the air density, which is calculated as follows: Where, p is the air pressure, and J is the gas constant; (11) Based on the power curve when not iced and its wind energy utilization coefficient C p(v,0) The power loss coefficient under different wind speeds v and ice thickness I is calculated: Wherein, ΔC p(v,I) is the power loss coefficient when the wind speed v, the ice thickness I; C p(v,I) is the wind energy utilization coefficient when the wind speed v, the ice thickness I; C p(v,0) is the wind energy utilization coefficient when the wind speed v, the ice thickness 0; (12) Calculate the power loss when the wind speed is v and the ice thickness is I: ΔP w(v,I) = ΔC p(v,I) × P w(v,0) P w(v,0) is the power output value of the fan when it is not iced and the wind speed is v.

Images