Yaw angle rapid optimization method and system for active yaw control

Abstract

The embodiment of the invention provides a yaw angle rapid optimization method and system for active yaw control, and the method comprises the steps: setting an initial yaw matrix formed by all fan yaw angles as a zero matrix, calculating the corresponding power of each fan traversing a preset fan yaw angle candidate set, and calculating the corresponding power of the fan yaw angle candidate set; selecting a first fan yaw angle corresponding to the maximum power of each fan to form a first yaw matrix of all fans; for each fan, after different offsets are added on the basis of the first fan yaw angle, the corresponding power is calculated, the offset value is in a preset offset angle candidate set, and a second fan yaw angle corresponding to the maximum power of each fan is selected to form a second yaw matrix of all the fans; and for each fan, different correction disturbance values are added on the basis of the second fan yaw angle, the power of the wind field is calculated, and a third fan yaw angle corresponding to the maximum power of each fan is selected to form a third yaw matrix of all the fans.

Description

Technical Field

[0001] This document relates to the field of wind power generation technology, and in particular to a method and system for rapid yaw angle optimization for active yaw control. Background Art

[0002] In a wind farm, the wake of upstream wind turbines can lead to a decrease in the aerodynamic performance and an increase in fatigue loads on downstream wind turbines. Active yaw control is a frequently used power augmentation method during the operation phase of a wind farm. Its advantage is that it does not require structural adjustments and is applicable to a wide range of situations. One of the key issues in active yaw is the optimization algorithm for the yaw angle of the wind turbine. Currently, commonly used active yaw control methods are mostly based on gradient optimization algorithms (such as genetic algorithms, gray wolf algorithms, etc.). These algorithms determine the optimal yaw angle for each wind turbine through iterative calculations. However, when the number of wind turbines in a wind farm is large, the time required for the optimization process increases significantly, making it difficult to meet the demand for real-time control of the yaw angle under conditions of multiple wind directions, multiple wind speeds, and multiple degrees of freedom. Summary of the Invention

[0003] One or more embodiments of this specification provide a method for rapid yaw angle optimization for active yaw control, including:

[0004] Rough calculation: Set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse the preset set of candidate wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the first yaw matrix of all wind turbines;

[0005] Accurate calculation: For each wind turbine, the corresponding power is calculated by adding different offsets to the first wind turbine yaw angle, where the offset value is taken from a preset offset angle candidate set. The second wind turbine yaw angle corresponding to the maximum power of each wind turbine is selected to form the second yaw matrix of all wind turbines;

[0006] Disturbance correction: For each wind turbine, different correction disturbance values ​​are added based on the yaw angle of the second wind turbine, the wind farm power is calculated, and the yaw angle of the third wind turbine corresponding to the maximum power of each wind turbine is selected to form the third yaw matrix of all wind turbines.

[0007] Furthermore, in the preset candidate set of wind turbine yaw angles, the wind turbine yaw angles are respectively set to -25°, -15°, -5°, 5°, 15° and 25°.

[0008] Furthermore, in the preset offset angle candidate set, the offset angle values ​​are -5°, -3°, -1°, 1°, 3° and 5° respectively.

[0009] Furthermore, the modified disturbance value is set to three random numbers between 0 and 1.

[0010] Furthermore, the specific method for calculating the power of a single wind turbine is as follows:

[0011] Sort the wind turbines by wind direction, add calculation grid points on the rotor plane according to the wind direction and the yaw angle of the wind turbine, calculate the approach turbulence and wind speed corresponding to each wind turbine, and obtain the wind speed at different grid points on the rotor plane;

[0012] The power of a single wind turbine is calculated based on the wind speed at different grid points on the rotor plane, as shown below:

[0013]

[0014] Among them, P i is the power of the i-th fan; ρ is the air density; A j is the representative area of ​​the wind wheel grid point; u j is the incoming wind speed at the j-th grid point; γ i is the yaw angle; α is a constant, which represents the exponential relationship between the power and the cosine of the yaw angle during yaw, and is generally taken as 2-3.

[0015] Furthermore, the specific calculation method of the wind speed at different grid points on the wind rotor plane is:

[0016] Set the turbulence and wind speed of the first wind turbine as the initial approach turbulence and wind speed. Then, traverse the second to the Nth wind turbines in sequence. In each traversal, calculate the additional turbulence and wind speed at the i-th wind turbine from the first i-1 upstream wind turbines and superimpose them. The impact of upstream wind turbines on the wind speed and turbulence of downstream wind turbines can be studied using different wake models, such as the Ishihara wake model. The following model can be used to superimpose the wakes of different wind turbines:

[0017]

[0018] Among them, u i and Ia i represents the wind speed and turbulence at fan i; u0 and Indicates the inflow wind speed and turbulence; u j represents the wind speed at fan j; u ij and ΔIa ij represents the wind speed and additional turbulence intensity of wind turbine j at the downstream wind turbine i; m is a constant, representing different wake superposition models.

[0019] Furthermore, the total wind farm power is calculated based on the power of each wind turbine, specifically:

[0020]

[0021] Among them, P T is the total wind farm power.

[0022] One or more embodiments of this specification provide a yaw angle rapid optimization system for active yaw control, including:

[0023] The first calculation module is used to set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse a preset set of candidate wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of the wind farm to form the first yaw matrix of all wind turbines;

[0024] A second calculation module is configured to calculate the corresponding power of each wind turbine by adding different offsets to the yaw angle of the first wind turbine, where the offset value is taken from a preset offset angle candidate set, and select the yaw angle of the second wind turbine corresponding to the maximum power of the wind farm to form a second yaw matrix for all wind turbines;

[0025] The third calculation module is used to add different correction disturbance values ​​to each wind turbine based on the yaw angle of the second wind turbine, calculate the wind farm power, and select the yaw angle of the third wind turbine corresponding to the maximum power of the wind farm to form the third yaw matrix of all wind turbines.

[0026] One or more embodiments of this specification provide an electronic device, including:

[0027] processor; and,

[0028] A memory is arranged to store computer executable instructions, which, when executed, cause the processor to implement the steps of the above-mentioned method for fast yaw angle optimization for active yaw control.

[0029] One or more embodiments of this specification provide a storage medium for storing computer-executable instructions, which, when executed, implement the steps of the above-mentioned method for rapid yaw angle optimization for active yaw control.

[0030] By adopting the embodiment of the present invention, the running time can be reduced by one order of magnitude while maintaining the same calculation results as the existing optimization algorithm, thereby significantly improving the calculation efficiency.

[0031] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are specifically listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] In order to more clearly illustrate one or more embodiments of this specification or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0033] Figure 1 A flowchart of a method for rapid yaw angle optimization for active yaw control provided in one or more embodiments of this specification;

[0034] Figure 2 A flow chart of a fast optimization algorithm for a yaw angle fast optimization method for active yaw control provided in one or more embodiments of this specification;

[0035] Figure 3 A schematic diagram of wind turbine layout and wind turbine numbering in a wind farm for a method for rapid yaw angle optimization for active yaw control provided in one or more embodiments of this specification;

[0036] Figure 4 A schematic diagram of wind rotor plane grid points for a method for rapid yaw angle optimization for active yaw control provided in one or more embodiments of this specification;

[0037] Figure 5 A schematic diagram of the composition of a yaw angle rapid optimization system for active yaw control provided in one or more embodiments of this specification;

[0038] Figure 6 A schematic diagram of the structure of an electronic device provided in one or more embodiments of this specification. DETAILED DESCRIPTION

[0039] In order to enable those skilled in the art to better understand the technical solutions in one or more embodiments of this specification, the technical solutions in one or more embodiments of this specification will be clearly and completely described below in conjunction with the drawings in one or more embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this document.

[0040] Method Example

[0041] According to an embodiment of the present invention, a method for rapid optimization of yaw angle for active yaw control is provided. Figure 1A flowchart of a method for rapid yaw angle optimization for active yaw control provided in one or more embodiments of this specification is shown in FIG. Figure 1 As shown, the yaw angle rapid optimization method for active yaw control according to an embodiment of the present invention specifically includes:

[0042] S1. Rough calculation: Set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse the preset set of candidate wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the first yaw matrix of all wind turbines.

[0043] First, in the wind farm, the wind turbines in the wind farm are numbered in a certain order according to the incoming wind direction. In this embodiment, the layout of the wind turbines in the wind farm is as follows: Figure 3 As shown in the figure, the incoming wind flows along the positive direction of the y-axis through a wind farm containing N wind turbines. The wind turbines are numbered from 1 to N according to the order in which the incoming wind flows through them, with 1 representing the upstreammost wind turbine and N representing the downstreammost wind turbine. When several wind turbines are arranged perpendicular to the incoming wind, they are numbered from left to right.

[0044] After the marking is completed, the yaw angle of the wind turbine is roughly calculated, such as Figure 2 As shown in the rough calculation process, the initial yaw matrix composed of the yaw angles of all wind turbines is set to zero matrix, and for the i-th wind turbine (1≤i≤N-1), the yaw angle γ i Take the value to traverse the candidate set, and then calculate the power corresponding to the i-th wind turbine at this time. The yaw angle with the largest power is recorded as the yaw angle of the first wind turbine γ ai Repeat the above process, after traversing the wind turbines (1≤i≤N-1), the yaw angles of all wind turbines are combined into the first yaw matrix, i.e., γ a Yaw matrix. In the preset candidate set of wind turbine yaw angles, the wind turbine yaw angles are set to -25°, -15°, -5°, 5°, 15°, and 25°.

[0045] S2. Accurate calculation: For each wind turbine, calculate the corresponding power by adding different offsets to the first wind turbine yaw angle, where the offset value is taken from a preset offset angle candidate set. The second wind turbine yaw angle corresponding to the maximum power of each wind turbine is selected to form the second yaw matrix of all wind turbines.

[0046] After roughly calculating the first yaw matrix, the first yaw angle in the matrix is ​​accurately calculated, such as Figure 2 As shown in the precise calculation process, for the i-th wind turbine (1≤i≤N-1), at the first wind turbine yaw angle γ a Based on this, let the yaw angle of the i-th wind turbine γ a Added offset Δγ i , Δγi Take the value to traverse the offset angle candidate set, and then calculate the power corresponding to the i-th wind turbine. The yaw angle when the power is maximum is recorded as the yaw angle of the second wind turbine γ bi Repeat the above process and traverse the wind turbines (1≤i≤N-1). The yaw angle combination of all wind turbines is the second yaw matrix, i.e., γ b Yaw matrix. In the preset offset angle candidate set, the offset angle values ​​are -5°, -3°, -1°, 1°, 3°, and 5°.

[0047] S3. Disturbance correction: For each wind turbine, add different correction disturbance values ​​based on the yaw angle of the second wind turbine, calculate the wind farm power, and select the third wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the third yaw matrix of all wind turbines.

[0048] After the accurate calculation is completed, the yaw angle is corrected again. For the i-th wind turbine (1≤i≤N-1), the yaw angle γ of the second wind turbine is b On the basis of , let the yaw angle of the i-th wind turbine add different correction disturbance values, calculate the power of the wind field under different yaw angles at this time, and the third yaw angle combination with the largest power is recorded as γ ci Repeat the above process, after traversing the wind turbines (1≤i≤N-1), the yaw angles of all wind turbines are combined into the third yaw matrix, namely γ c Yaw matrix, where the correction perturbation value is set to three random numbers between 0 and 1.

[0049] In the above three calculation processes, the specific calculation method for the power of a single fan is as follows:

[0050] First, sort the wind turbines according to the wind direction, and add calculation grid points on the wind rotor plane according to the wind direction and the yaw angle of the wind turbine, such as Figure 4 As shown, the total number of grid points is M. The turbulence and wind speed of the first wind turbine are set as the initial approach turbulence and wind speed. Then, the 2nd to Nth wind turbines are traversed in sequence. In each traversal, the additional turbulence and wind speed of the first i-1 wind turbines upstream at the i-th wind turbine are calculated and superimposed. The influence of the upstream wind turbine on the wind speed and turbulence of the downstream wind turbine can be studied using different wake models, such as the Ishihara wake model. The following model can be used to superimpose the wakes of different wind turbines:

[0051]

[0052] Among them, u i and Ia i represents the wind speed and turbulence at fan i; u0 and Indicates the inflow wind speed and turbulence; u j represents the wind speed at fan j; u ij and ΔIaij represents the wind speed and additional turbulence intensity of wind turbine j at the downstream wind turbine i; m is a constant, representing different wake superposition models, such as the linear superposition value is 1, and the square model addition value is 2.

[0053] According to the above wake model and wake superposition model, the wind speed at different grid points on the rotor plane can be obtained. The wind turbine power can be calculated based on the wind speed at different grid points on the rotor plane as shown below:

[0054]

[0055] Among them, P i is the power of the i-th wind turbine; ρ is the air density; Aj is the representative area of ​​the wind wheel grid point; u j is the incoming wind speed at the j-th grid point; γ i is the yaw angle; α is a constant, which represents the exponential relationship between the power and the cosine of the yaw angle during yaw, and is generally taken as 2-3.

[0056] The total power of the wind farm is calculated based on the power of a single wind turbine, specifically:

[0057]

[0058] Among them, P T is the total wind farm power.

[0059] The beneficial effects of the present invention are as follows:

[0060] By adopting the embodiment of the present invention, the running time can be reduced by one order of magnitude while maintaining the same calculation results as the existing optimization algorithm, thereby significantly improving the calculation efficiency.

[0061] System Example

[0062] According to an embodiment of the present invention, a yaw angle rapid optimization system for active yaw control is provided. Figure 5 A schematic diagram of a yaw angle rapid optimization system for active yaw control provided in one or more embodiments of this specification is shown in FIG. Figure 5 As shown, the yaw angle rapid optimization system for active yaw control according to an embodiment of the present invention specifically includes:

[0063] A first calculation module 50 is configured to set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse a preset set of candidate wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of the wind farm to form the first yaw matrix of all wind turbines;

[0064] A second calculation module 52 is configured to calculate the corresponding power of each wind turbine by adding different offsets to the first wind turbine yaw angle, wherein the offsets are taken from a preset offset angle candidate set, and select the second wind turbine yaw angle corresponding to the maximum power of the wind farm to form a second yaw matrix for all wind turbines;

[0065] The third calculation module 54 is used to add different correction disturbance values ​​to the yaw angle of the second wind turbine for each wind turbine, calculate the wind farm power, and select the third wind turbine yaw angle corresponding to the maximum wind farm power to form a third yaw matrix for all wind turbines.

[0066] The embodiment of the present invention is a system embodiment corresponding to the above-mentioned method embodiment. The specific operations of each module can be understood by referring to the description of the method embodiment, which will not be repeated here.

[0067] Device Example 1

[0068] An embodiment of the present invention provides an electronic device, such as Figure 6 As shown, it includes: a memory 60, a processor 62, and a computer program stored in the memory 60 and executable on the processor 62. When the computer program is executed by the processor 62, the following method steps are implemented:

[0069] S1. Rough calculation: Set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix. Calculate the corresponding power of each wind turbine by traversing a preset set of candidate wind turbine yaw angles. Select the first wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the first yaw matrix of all wind turbines.

[0070] S2. Accurate calculation: For each wind turbine, the corresponding power is calculated by adding different offsets to the first wind turbine yaw angle, where the offset values ​​are taken from a preset set of candidate offset angles. The second wind turbine yaw angle corresponding to the maximum power of each wind turbine is selected to form a second yaw matrix for all wind turbines;

[0071] S3. Disturbance correction: For each wind turbine, add different correction disturbance values ​​based on the yaw angle of the second wind turbine, calculate the wind farm power, and select the third wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the third yaw matrix of all wind turbines.

[0072] Device Example 2

[0073] An embodiment of the present invention provides a computer-readable storage medium having stored thereon a program for implementing information transmission. When the program is executed by the processor 62, the following method steps are implemented:

[0074] S1. Rough calculation: Set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix. Calculate the corresponding power of each wind turbine by traversing a preset set of candidate wind turbine yaw angles. Select the first wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the first yaw matrix of all wind turbines.

[0075] S2. Accurate calculation: For each wind turbine, the corresponding power is calculated by adding different offsets to the first wind turbine yaw angle, where the offset values ​​are taken from a preset set of candidate offset angles. The second wind turbine yaw angle corresponding to the maximum power of each wind turbine is selected to form a second yaw matrix for all wind turbines;

[0076] S3. Disturbance correction: For each wind turbine, add different correction disturbance values ​​based on the yaw angle of the second wind turbine, calculate the wind farm power, and select the third wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the third yaw matrix of all wind turbines.

[0077] The computer-readable storage medium in this embodiment includes, but is not limited to, ROM, RAM, magnetic disk, or optical disk.

[0078] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the above embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for rapid optimization of yaw angle for active yaw control, characterized in that: include: Rough calculation: Set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse the preset candidate set of wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of each wind turbine to form the first yaw matrix of all wind turbines; Accurate calculation: for each wind turbine, different offsets are added to the yaw angle of the first wind turbine to calculate the corresponding power, where the offset value is taken from a preset offset angle candidate set, and the second wind turbine yaw angle corresponding to the maximum power of each wind turbine is selected to form the second yaw matrix of all wind turbines; Disturbance correction: For each wind turbine, different correction disturbance values ​​are added based on the yaw angle of the second wind turbine, the wind farm power is calculated, and the yaw angle of the third wind turbine corresponding to the maximum power of each wind turbine is selected to form the third yaw matrix of all wind turbines.

2. The method according to claim 1, characterized in that In the preset candidate set of wind turbine yaw angles, the wind turbine yaw angles are respectively set to -25°, -15°, -5°, 5°, 15° and 25°.

3. The method according to claim 1, characterized in that In the preset offset angle candidate set, the offset angle values ​​are -5°, -3°, -1°, 1°, 3° and 5° respectively.

4. The method according to claim 1, characterized in that: The modified disturbance value is set to three random numbers between 0 and 1.

5. The method according to claim 1, characterized in that The specific method for calculating the power of a single fan is as follows: Sort the wind turbines by wind direction, add calculation grid points on the wind rotor plane according to the wind direction and the yaw angle of the wind turbine, calculate the approach turbulence and wind speed corresponding to each wind turbine, and obtain the wind speed at different grid points on the wind rotor plane; The wind turbine power is calculated based on the wind speed at different grid points on the rotor plane as shown below: Among them, P i : power of the i-th wind turbine; ρ: air density; Aj: representative area of ​​the wind wheel grid point; u j : wind speed at the jth grid point; γ i : Yaw angle; α: constant, which indicates the exponential relationship between power and cosine of yaw angle during yaw, and its value is usually 2-3.

6. The method according to claim 5, characterized in that The specific calculation method of wind speed at different grid points on the wind rotor plane is: Set the turbulence and wind speed of the first wind turbine as the initial approach turbulence and wind speed; then, traverse the second to Nth wind turbines in turn. In each traversal, calculate the additional turbulence and wind speed of the first i-1 wind turbines upstream at the i-th wind turbine and superimpose them; the influence of the upstream wind turbine on the wind speed and turbulence of the downstream wind turbine can be measured by different wake models such as the Ishihara wake model. The following model can be used for the superposition of wakes of different wind turbines: Among them, u i and Ia i represents the wind speed and turbulence at fan i; u0 and Indicates the inflow wind speed and turbulence; u j represents the wind speed at fan j; u ij and ΔIa ij represents the wind speed and additional turbulence intensity of wind turbine j at the downstream wind turbine i; m is a constant, representing different wake superposition models.

7. The method according to claim 5, characterized in that The total power of the wind farm is calculated based on the power of a single wind turbine, specifically: Among them, P T is the total wind farm power.

8. A yaw angle rapid optimization system for active yaw control, characterized in that: include: The first calculation module is used to set the initial yaw matrix formed by the yaw angles of all wind turbines to a zero matrix, traverse the preset candidate set of wind turbine yaw angles for each wind turbine to calculate the corresponding power, and select the first wind turbine yaw angle corresponding to the maximum power of the wind farm to form the first yaw matrix of all wind turbines; The second calculation module is used to calculate the corresponding power of each wind turbine after adding different offsets on the basis of the yaw angle of the first wind turbine, wherein the offset value is taken from a preset offset angle candidate set, and the yaw angle of the second wind turbine corresponding to the maximum power of the wind farm is selected to form a second yaw matrix of all wind turbines; The third calculation module is used to add different correction disturbance values ​​to each wind turbine based on the yaw angle of the second wind turbine, calculate the wind farm power, and select the yaw angle of the third wind turbine corresponding to the maximum power of the wind farm to form the third yaw matrix of all wind turbines.

9. An electronic device, characterized in that: include: processor; as well as, A memory arranged to store computer executable instructions, which, when executed, cause the processor to implement the steps of the method for fast yaw angle optimization for active yaw control according to any one of claims 1 to 7.

10. A storage medium, characterized in that: Used to store computer executable instructions, which, when executed, implement the steps of the yaw angle rapid optimization method for active yaw control according to any one of claims 1 to 7.

Images