A wind turbine incoming flow wind speed acquisition method, device, equipment and storage medium

By measuring the wind speed and direction angle in front of the wind turbine using lidar, calculating the standard deviation and function coefficients of the dimensionless wind speed loss, and reconstructing the velocity of the incoming flow field in front of the wind turbine, the problem of obtaining the incoming wind speed at different positions in front of the wind turbine is solved, thus improving the accuracy of wind turbine control and load assessment.

CN115220065BActive Publication Date: 2025-11-07华能通渭风电有限责任公司
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202210832840.9
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2022-07-15
Publication Date
2025-11-07
Estimated Expiration
2042-07-15

AI Technical Summary

Technical Problem

Existing technologies make it difficult to accurately obtain the incoming wind speed at different locations in front of the wind turbine, which affects the wind turbine's control strategy and load assessment.

Method used

By measuring the horizontal wind speed and wind direction angle at various set locations directly in front of the wind turbine using lidar, the standard deviation and function coefficients of the dimensionless wind speed loss are calculated. The incoming flow velocity reconstruction function of the blocking effect zone in front of the wind turbine is then reconstructed to calculate the incoming wind speed at any location.

Benefits of technology

It enables accurate acquisition of the incoming wind speed at any location within the blocking effect zone in front of the wind turbine, supporting the calculation of wind turbine blade load changes and unit feedforward control.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN115220065B_ABST
    Figure CN115220065B_ABST
Patent Text Reader

Abstract

The application discloses a wind turbine incoming flow wind speed acquisition method, device, equipment and readable storage medium. The method comprises the following steps: acquiring the horizontal wind speed and wind direction angle of different height layers at each first set position and second set position in front of the wind turbine; calculating the standard deviation of the dimensionless wind speed loss of each first set position; calculating the function coefficient according to the standard deviation of the dimensionless wind speed loss and the function relationship between the standard deviation and the first set position; and reconstructing the incoming flow field speed reconstruction function of the front blocking effect zone of the wind turbine according to the function relationship, the wind shear exponent of the far incoming flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the incoming flow wind speed in front of the wind turbine and the far incoming flow wind speed, so as to calculate the incoming flow wind speed at any position in front of the wind turbine. The technical scheme disclosed by the application realizes the calculation of the incoming flow wind speed at any position in the front blocking effect zone of the wind turbine by reconstructing the incoming flow field speed reconstruction function of the front blocking effect zone of the wind turbine.
Need to check novelty before this filing date? Find Prior Art

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of wind power generators, in particular to a wind turbine incoming flow wind speed acquisition method, device, equipment and readable storage medium. BACKGROUND

[0002] A wind power generator converts the kinetic energy of wind into electric energy by rotating a wind wheel to drive a generator. Therefore, wind speed is crucial to the power generation of a wind power generator.

[0003] Currently, in a wind farm, the wind speed in front of a wind power generator is measured by a few wind measurement towers, and the data of 2-3 height layers are measured. However, the wind measurement towers are generally far away from the wind power generator, and the nearest distance is generally not less than 2.5 times the diameter of the wind wheel. In addition, the highest measurement height of the wind measurement tower is generally the same as the hub center of the wind turbine. Therefore, the limited number of height layers cannot obtain the wind speed at any spatial position in front of the wind turbine, which will affect the control strategy research and load evaluation of the wind turbine.

[0004] In summary, how to obtain the incoming flow wind speed at different positions in front of the wind turbine is a technical problem to be solved by the technical personnel in the field. SUMMARY

[0005] Therefore, the purpose of the present application is to provide a wind turbine incoming flow wind speed acquisition method, device, equipment and readable storage medium for obtaining the incoming flow wind speed at different positions in front of the wind turbine.

[0006] In order to achieve the above purpose, the present application provides the following technical scheme:

[0007] A wind turbine incoming flow wind speed acquisition method, comprising:

[0008] acquiring the horizontal wind speed and wind direction angle of each first set position in front of the wind turbine measured by a laser radar, a second set position as a far incoming flow wind speed position, and different height layers;

[0009] calculating the standard deviation of the dimensionless wind speed loss of each first set position according to the horizontal wind speed and wind direction angle of different height layers at the second set position and each first set position;

[0010] calculating a function coefficient according to the standard deviation of the dimensionless wind speed loss of each first set position and the functional relationship between the first set position and the function coefficient;

[0011] According to the function relationship of the calculated function coefficient, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow speed in front of the wind turbine and the far flow speed, a flow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed to calculate the flow speed at any position in the front blocking effect zone of the wind turbine by using the flow field speed reconstruction function.

[0012] Preferably, the standard deviation of the dimensionless wind speed loss of each first setting position is calculated according to the horizontal wind speed and the wind direction angle of different height layers at the second setting position and each first setting position, comprising:

[0013] The axial wind speed at different height layers at the second setting position is calculated by using the horizontal wind speed and the wind direction angle of different height layers at the second setting position.

[0014] The axial wind speed at different height layers at each first setting position is calculated by using the horizontal wind speed and the wind direction angle of different height layers at each first setting position.

[0015] The wind speed loss at different height layers at each first setting position is calculated by using the axial wind speed at different height layers at the second setting position and the axial wind speed at different height layers at each first setting position.

[0016] The dimensionless wind speed loss at different height layers at each first setting position is obtained by dimensionless processing the wind speed loss at different height layers at each first setting position by using the axial wind speed at the hub center height layer at the second setting position.

[0017] The standard deviation of the dimensionless wind speed loss of each first setting position is calculated according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each first setting position.

[0018] Preferably, the standard deviation of the dimensionless wind speed loss of each first setting position is calculated according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each first setting position, comprising:

[0019] According to The standard deviation of the dimensionless wind speed loss of each first setting position is calculated by using the least square method fitting.

[0020] Wherein, Δu (x,0,z) is the wind speed loss at z height layer at x position, u (SSP,0,H) is the axial wind speed at the hub center height layer at the second setting position, D is the diameter of the wind wheel, H is the hub center height, is the dimensionless wind speed loss at z height layer at x position, CT is the thrust coefficient of the wind turbine, σ x is the standard deviation of the dimensionless wind speed deficit at the x position.

[0021] Preferably, the function coefficients are calculated according to the standard deviation of the dimensionless wind speed deficit at each of the first set positions and a functional relationship between the standard deviation and the first set positions, including:

[0022] the standard deviation of the dimensionless wind speed deficit at each of the first set positions and a functional relationship The function coefficients λ and η are calculated by least square fitting.

[0023] Preferably, the wind shear exponent of the far flow of the wind turbine is obtained, including:

[0024] The wind shear exponent of the far flow of the wind turbine is obtained by: wherein α is the wind shear exponent of the far flow of the wind turbine, u (SSP,0,z) is the axial wind speed at the z height layer at the second set position, u (SSP,0,ref) is the axial wind speed at the h (SSP,0,ref) height layer at the second set position.

[0025] Preferably, the inflow field velocity reconstruction function of the front blockage effect zone of the wind turbine is reconstructed according to the functional relationship of the calculated function coefficients, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed deficit, and the relationship between the inflow wind speed directly in front of the wind turbine and the far flow wind speed, including:

[0026] According to and the relationship between the inflow wind speed directly in front of the wind turbine and the far flow wind speed u (x,0,z) = u (SSP,0,z) - Δu (x,0,z) , the inflow field velocity reconstruction function of the front blockage effect zone of the wind turbine is: ; wherein u (x,y,z) is the inflow wind speed at (x, y, z) in the front blockage effect zone of the wind turbine.

[0027] A wind turbine inflow wind speed obtaining device, comprising:

[0028] An obtaining module, configured to obtain the horizontal wind speed and wind direction angle at each first set position directly in front of the wind turbine and at a second set position as a far flow wind speed position, measured by a laser radar.

[0029] A first calculating module, configured to calculate the standard deviation of the dimensionless wind speed deficit at each of the first set positions according to the horizontal wind speed and wind direction angle at different height layers at the second set position and at each of the first set positions.

[0030] a second calculation module, configured to calculate a function coefficient according to a standard deviation of the dimensionless wind speed loss of each of the first set positions and a function relationship between the first set positions;

[0031] a obtaining module, configured to obtain a flow field velocity reconstruction function of a front blockage effect zone of the wind turbine according to the function relationship of the calculated function coefficient, a wind shear exponent of a far flow of the wind turbine, a Gaussian function satisfied by the dimensionless wind speed loss, and a relationship between a flow speed in front of the wind turbine and the far flow speed of the wind turbine, so as to calculate the flow speed at any position in the front blockage effect zone of the wind turbine by using the flow field velocity reconstruction function.

[0032] Preferably, the first calculation module comprises:

[0033] a first calculation unit, configured to calculate the axial wind speed at different height layers at the second set position by using the horizontal wind speed and the wind direction angle at different height layers at the second set position;

[0034] a second calculation unit, configured to correspondingly calculate the axial wind speed at different height layers at each of the first set positions by using the horizontal wind speed and the wind direction angle at different height layers at each of the first set positions;

[0035] a third calculation unit, configured to correspondingly calculate the wind speed loss at different height layers at each of the first set positions by using the axial wind speed at different height layers at the second set position and the axial wind speed at different height layers at each of the first set positions;

[0036] a first obtaining unit, configured to obtain the dimensionless wind speed loss at different height layers at each of the first set positions by using the axial wind speed at the hub center height layer at the second set position to perform dimensionless processing on the wind speed loss at different height layers at each of the first set positions;

[0037] a fourth calculation unit, configured to calculate the standard deviation of the dimensionless wind speed loss of each of the first set positions according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each of the first set positions.

[0038] A wind turbine flow speed acquisition device, comprising:

[0039] a memory, configured to store a computer program;

[0040] a processor, configured to execute the computer program to realize the steps of the wind turbine flow speed acquisition method according to any one of the above.

[0041] A readable storage medium, in which a computer program is stored, the computer program, when executed by a processor, implements the steps of the wind turbine inflow wind speed acquisition method of any one of the above.

[0042] The application provides a wind turbine inflow wind speed acquisition method, device, equipment and readable storage medium, wherein the method comprises the following steps: acquiring horizontal wind speeds and wind direction angles of different height layers at each first set position in front of a wind turbine and at a second set position serving as a far-inflow wind speed position measured by a laser radar; calculating standard deviations of dimensionless wind speed deficits of each first set position according to the horizontal wind speeds and wind direction angles of different height layers at the second set position and at each first set position; calculating a function coefficient according to a function relationship between the standard deviations of the dimensionless wind speed deficits of each first set position and the first set positions; and reconstructing an inflow field speed reconstruction function of a front blocking effect zone of the wind turbine according to a function relationship of the calculated function coefficient, a wind shear exponent of far-inflow of the wind turbine, a Gaussian function satisfied by the dimensionless wind speed deficit, and a relationship between the inflow wind speed in front of the wind turbine and the far-inflow wind speed, so as to calculate the inflow wind speed at any position in the front blocking effect zone of the wind turbine by using the inflow field speed reconstruction function.

[0043] The above technical solution disclosed by the application uses the laser radar to measure the horizontal wind speeds and wind direction angles of different height layers at each first set position in front of the wind turbine and at the second set position serving as the far-inflow wind speed position, then calculates the standard deviations of the dimensionless wind speed deficits of each first set position according to the horizontal wind speeds and wind direction angles of different height layers at the second set position and at each first set position, and calculates the function coefficient of the function relationship between the standard deviations of the dimensionless wind speed deficits of each first set position and the first set positions according to the calculated standard deviations of the dimensionless wind speed deficits of each first set position, then reconstructs the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine according to the function relationship of the calculated function coefficient, the wind shear exponent of far-inflow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed deficit, and the relationship between the inflow wind speed in front of the wind turbine and the far-inflow wind speed, so that the inflow wind speed at any position in the front blocking effect zone of the wind turbine can be calculated by using the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine, thereby enabling relevant personnel to better understand the blocking effect zone of the wind turbine, and providing favorable support for calculating the load changes of the wind turbine blade along different spanwise directions, and providing an effective reference for the front feed control of the unit. BRIEF DESCRIPTION OF DRAWINGS

[0044] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the accompanying drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced. Obviously, the accompanying drawings in the following description only represent some of the embodiments of the present application, and for those skilled in the art, other drawings can be obtained without creative labor on the basis of the provided drawings.

[0045] Figure 1 A flow chart of a wind turbine incoming flow wind speed acquisition method provided by an embodiment of the present application;

[0046] Figure 2 A structural schematic diagram of a wind turbine incoming flow wind speed acquisition device provided by an embodiment of the present application;

[0047] Figure 3 A structural schematic diagram of a wind turbine incoming flow wind speed acquisition device provided by an embodiment of the present application. DETAILED DESCRIPTION

[0048] The existing wind turbine nacelle wind speed instrument is a single-point speed measurement at the position of the top of the nacelle behind the wind turbine wind wheel, and it is difficult to determine the wind speed in front of the wind turbine by using the data. At present, in a wind farm, the wind speed in front of the wind turbine is generally measured by a few wind measurement towers, and the real-time wind profile change can be obtained by the data of 2-3 height layers. However, the wind measurement tower is generally far away from the wind turbine, and the nearest distance is generally not less than 2.5 times the diameter of the wind wheel, so the wind speed flow state between the wind measurement tower and the wind turbine is unknown, and the closer to the wind turbine, the greater the blocking effect of the wind speed by the wind turbine. The widely used analytical expression speed prediction method based on the "vortex sheet theory" only establishes a model at one height (along the hub center height of the wind turbine) to obtain the wind speed in front of the wind turbine, and the wind speed at other heights or at any plane position in front of the wind turbine cannot be obtained.

[0049] Therefore, the present application provides a wind turbine incoming flow wind speed acquisition method, device, equipment and readable storage medium, which are used to acquire the incoming flow wind speed at different positions in front of the wind turbine.

[0050] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0051] Reference is made to Figure 1A flow chart of a wind turbine incoming flow wind speed acquisition method is shown, and the wind turbine incoming flow wind speed acquisition method can include the following steps.

[0052] S11: acquiring the horizontal wind speed and wind direction angle of each first set position in front of the wind turbine measured by the laser radar, as the second set position of the far incoming flow wind speed, the horizontal wind speed and wind direction angle of different height layers.

[0053] In the present application, first, the ground laser radar wind measurement technology can be used to measure the horizontal wind speed and wind direction angle of each first set position in front of the wind turbine, and the wind speed and wind direction angle of different height layers at the second set position in front of the wind turbine as the far incoming flow wind speed position. The ground laser radar can measure the horizontal wind speed of multiple height layers (usually 20 layers) vertically upward at one position, and the measurement height range is generally 40-300 m, which covers the top and bottom positions of the tip of the large wind turbine. Considering that the laser radar is expensive and mobile equipment, therefore, in a wind farm, a large number of equipment is not required, only 1-2 sets are required, and the laser radar only needs to be moved for specific measurement positions, so as to reduce the cost.

[0054] For the same wind turbine, when a laser radar is deployed to measure the horizontal wind speed and wind direction angle, the laser radar is first used to measure the horizontal wind speed and wind direction angle of each height layer at the first first set position. Then, the laser radar is moved to the second first set position in front of the wind turbine to measure the horizontal wind speed and wind direction angle of different height layers at the second first set position, and then the laser radar is moved to the third first set position in front of the wind turbine, and so on, until the measurement of the horizontal wind speed and wind direction angle of different height layers at the last first set position in front of the wind turbine is completed. Finally, the laser radar can be moved to the second set position in front of the wind turbine to measure the horizontal wind speed and wind direction angle of different height layers at the second set position. Of course, the measurement of the horizontal wind speed and wind direction angle of different height layers at the second set position can be performed before or after the measurement of any first set position, which is not limited in the present application.

[0055] For the same wind turbine, when two lidars are deployed, one of the lidars is used to measure the horizontal wind speed and wind direction angle of different height layers at the second set position, and the other is used to measure the horizontal wind speed and wind direction angle of different height layers at each first set position in turn by moving (i.e., after measuring the horizontal wind speed and wind direction angle of different height layers at the first first set position, moving to the next first set position for measurement, until the measurement of the last first set position is completed), and the two lidars can measure at the same time, that is, when the first lidar measures the horizontal wind speed and wind direction angle of different height layers at the second set position, the second lidar can measure the horizontal wind speed and wind direction angle of the corresponding height layer at the first set position being measured.

[0056] It should be noted that in order to reduce the measurement error and improve the accuracy of the wind turbine incoming flow wind speed, two lidars are preferably deployed to measure at the same time. In addition, the lidar can measure the horizontal wind speed and wind direction angle at the same time. Moreover, in order to improve the accuracy of the wind turbine incoming flow wind speed, the measurement time of the lidar at each height layer at each set position covers at least all operating states of the unit (from the start of the cut-in wind speed to the stop of the cut-out wind speed), and the lidar can measure the horizontal wind speed and wind direction angle of different height layers at the current position every second (of course, the measurement time interval can be set according to actual needs), and the horizontal wind speed and wind direction angle of different height layers at the current position output by the lidar are the horizontal wind speed and wind direction angle obtained by averaging the measured horizontal wind speed and wind direction angle of different height layers at the current position (of course, the lidar can also output all the measured horizontal wind speed and wind direction angle of different height layers at the current set position, and then the wind turbine incoming flow wind speed acquisition device can calculate the average to obtain the horizontal wind speed and wind direction angle of different height layers at the current set position).

[0057] The first set position mentioned above can be 0.5D, 1D, 1.5D, 2D, 2.5D, and the second set position can be 3D, that is, the lidar can measure the horizontal wind speed and wind direction angle of different height layers at 0.5D, 1D, 1.5D, 2D, 2.5D, and 3D positions in front of the wind turbine, and D is the diameter of the wind wheel. Of course, the first set position and the second set position can be adjusted according to actual needs, but the maximum value of the first set position is less than the second set position, and the second set position is preferably not more than 4D, and the first set position and the second set position are preferably uniformly distributed in front of the wind turbine.

[0058] On the basis, the wind turbine incoming flow wind speed acquisition device can acquire the horizontal wind speed and wind direction angle of each first setting position at different height layers in front of the wind turbine measured by the laser radar, and the horizontal wind speed and wind direction angle of each first setting position at different height layers in front of the wind turbine as the far incoming flow wind speed position.

[0059] S12: According to the horizontal wind speed and wind direction angle of each first setting position at different height layers and the second setting position, the standard deviation of the dimensionless wind speed loss of each first setting position is calculated.

[0060] In the wake of the wind turbine, the axial wind speed loss satisfies the Gaussian self-similarity, so from the flow condition of the wind turbine, the wind speed decay area with blocking effect in front of the wind turbine is more suitable as the research object of the speed loss. At the same time, the speed loss in front of the wind turbine also has self-similarity. The present application assumes that the wind speed loss of each position in the area with blocking effect in front of the wind turbine satisfies the Gaussian distribution.

[0061] On the basis of step S11, the wind speed loss of each first setting position can be calculated according to the horizontal wind speed and wind direction angle of each first setting position at different height layers and the second setting position, and the far incoming flow wind speed is dimensionless. Since the laser radar cannot measure the wind speed of each height layer of multiple setting positions at the same time, the present application calculates the standard deviation of the dimensionless wind speed loss of each first setting position based on the Gaussian distribution satisfied by the wind speed loss according to the horizontal wind speed and wind direction angle of each first setting position at different height layers and the second setting position measured by the laser radar, so as to reconstruct the wind turbine incoming flow field speed reconstruction function based on the standard deviation of the dimensionless wind speed loss of each first setting position.

[0062] It should be noted that in front of the wind turbine, when the distance from the wind turbine is from far to near, the wind speed is also from small to large due to the blocking effect of the wind turbine. Therefore, when the distance from the wind turbine is from far to near, the standard deviation of the dimensionless wind speed loss decreases from large to small.

[0063] S13: According to the standard deviation of the dimensionless wind speed loss of each first setting position and the functional relationship between the standard deviation of the dimensionless wind speed loss of each first setting position and the first setting position, the function coefficient is calculated.

[0064] After calculating the standard deviation of the dimensionless wind speed loss of each first setting position, the function coefficient in the functional relationship can be calculated according to the standard deviation of the dimensionless wind speed loss of each first setting position and the functional relationship between the standard deviation of the dimensionless wind speed loss of each first setting position and the first setting position (i.e. the functional relationship between the standard deviation of the dimensionless wind speed loss and the distance of the position in front of the wind turbine), so as to reconstruct the incoming flow field speed reconstruction function of the wind turbine according to the functional relationship of the calculated function coefficient.

[0065] S14: According to the function relationship of the calculated function coefficient, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow speed in front of the wind turbine and the far flow speed, the flow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed to calculate the flow speed at any position in the front blocking effect zone of the wind turbine.

[0066] On the basis of step S13, according to the function relationship of the calculated function coefficient, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow speed in front of the wind turbine and the far flow speed (the relationship is that the flow speed in front of the wind turbine at the x position of the z height layer is the difference between the far flow speed (the wind speed at the z height of the second set position) and the wind speed loss at the x position of the z height layer), the relationship between the flow speed in front of the wind turbine and the far flow speed is extended and reconstructed in the two-dimensional yz plane at the x position in front of the wind turbine, and the relationship between the flow speed in front of the wind turbine and the far flow speed is extended and reconstructed at any x position (x<second set position) in front of the wind turbine by the function relationship corresponding to the standard deviation of the dimensionless wind speed loss, thereby forming a three-dimensional space flow field function of the front blocking effect zone of the wind turbine, i.e. reconstructing the flow speed field reconstruction function of the front blocking effect zone of the wind turbine, so that the flow speed at any position in the front blocking effect zone of the wind turbine can be calculated by using the reconstructed flow speed field reconstruction function of the front blocking effect zone of the wind turbine. Specifically, the (x, y, z) position in the front blocking effect zone of the wind turbine to be calculated can be brought into the reconstructed flow speed field reconstruction function of the front blocking effect zone of the wind turbine to obtain the flow speed at the (x, y, z) position.

[0067] Through the above-mentioned manner, the flow speed at any position in front of the wind turbine is effectively simulated by the wind speed at a limited number of discrete point positions in front of the wind turbine.

[0068] The above technical solution disclosed in the application utilizes the laser radar to measure the horizontal wind speed and the wind direction angle of different height layers at each first set position in front of the wind turbine and at the second set position as the far flow wind speed position, then calculates the standard deviation of the dimensionless wind speed loss of each first set position according to the horizontal wind speed and the wind direction angle of different height layers at the second set position and at each first set position, calculates the function coefficient of the function relationship between the standard deviation of the dimensionless wind speed loss and the first set position according to the calculated standard deviation of the dimensionless wind speed loss of each first set position, and then reconstructs the flow field speed reconstruction function of the front blocking effect zone of the wind turbine according to the function relationship of the calculated function coefficient, the obtained wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow wind speed in front of the wind turbine and the far flow wind speed, so that the flow wind speed at any position in the front blocking effect zone of the wind turbine can be calculated by using the flow field speed reconstruction function of the front blocking effect zone of the wind turbine, thereby enabling relevant personnel to better understand the blocking effect zone of the wind turbine and provide favorable support for calculating the load change of the blade of the wind turbine along different spanwise directions, and providing an effective reference for the front feed control of the unit.

[0069] The wind turbine flow wind speed acquisition method provided by the embodiment of the application calculates the standard deviation of the dimensionless wind speed loss of each first set position according to the horizontal wind speed and the wind direction angle of different height layers at the second set position and at each first set position, which can include:

[0070] The horizontal wind speed and the wind direction angle of different height layers at the second set position are utilized to calculate the axial wind speed at different height layers at the second set position.

[0071] The horizontal wind speed and the wind direction angle of different height layers at each first set position are utilized to correspondingly calculate the axial wind speed at different height layers at each first set position.

[0072] The axial wind speed at different height layers at the second set position and the axial wind speed at different height layers at each first set position are utilized to correspondingly calculate the wind speed loss at different height layers at each first set position.

[0073] The axial wind speed at the hub center height layer at the second set position is utilized to perform dimensionless processing on the wind speed loss at different height layers at each first set position, to obtain the dimensionless wind speed loss at different height layers at each first set position.

[0074] The standard deviation of the dimensionless wind speed loss of each first set position is calculated according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each first set position.

[0075] In this application, the specific process for calculating the standard deviation of the dimensionless wind speed loss at each of the first set locations based on the horizontal wind speed and wind direction angle at different heights at the second set location and each of the first set locations is as follows:

[0076] 1) Utilizing the horizontal wind speed and wind direction angle at different height levels at the second predetermined location, using u (SSP,0,z) =U (SSP,0,z) ·cos(θ (SSP,0,z) Calculate the axial wind speed u at the second predetermined position directly in front of the wind turbine, at a height of z. (SSP,0,z) (This wind speed is taken as the incoming wind speed from a distance); where SSP is the second set position, U (SSP,0,z) θ is the horizontal wind speed at the z-height level of the second predetermined location measured by lidar. (SSP,0,z) The wind direction angle at the second set location (Z-height layer) measured by the lidar is the difference between the horizontal wind direction relative to true north and the azimuth angle of the wind turbine, z∈(0,300), and z is a custom value based on the center height of the wind turbine hub and the diameter of the rotor. y=0 indicates that the wind turbine is directly in front. Taking the second set location (SSP) as 3D as an example, this is where u is used. (3D,0,z) =U (3D,0,z) ·cos(θ (3D,0,z) Calculate the axial wind speed at the z-height layer at the 3D location.

[0077] 2) Utilizing the horizontal wind speed and wind direction angle at each first predetermined location x at a height z, using u (x,0,z) =U (x,0,z) ·cos(θ (x,0,z) The axial wind speed u at the first set position x, z-height level directly in front of the wind turbine is calculated accordingly. (x,0,z) Among them, U (x,0,z) Let θ be the horizontal wind speed at the first predetermined location x at the z-height level measured by the lidar. (x,0,z) The wind direction angle at the z-height layer of the first set location x, as measured by the lidar (i.e., the difference between the horizontal wind direction relative to true north and the azimuth angle of the wind turbine); to save costs, x = 2.5D, 2D, 1.5D, 1D, 0.5D are taken;

[0078] 3) Utilizing the axial wind speeds at different height levels at the second set position and at different height levels at the first set position, Δu is used. (x,0,z) =u (SSP,0,z) -u (x,0,z) Calculate the wind speed loss Δu at the z-height layer of the first set location x. (x,0,z) Taking the second set position SSP as 3D as an example, then Δu is used at this time. (x,0,z) =u (3D,0,z) -u (x,0,z) Calculate the wind speed loss Δu at the z-height layer of the first set location x.(x,0,z) .

[0079] 4) using the axial wind speed u at the hub center height layer at the second set position (SSP,0,H) The wind speed deficit at different height layers at each first set position is dimensionless, and the dimensionless wind speed deficit at different height layers at each first set position is obtained (specifically, the axial wind speed u at the hub center height layer at the second set position is used (SSP,0,H) The wind speed deficit at different height layers at each first set position is dimensionless, and the dimensionless wind speed deficit at different height layers at each first set position is obtained (specifically, the axial wind speed u at the hub center height layer at the second set position is used wherein, is the dimensionless wind speed deficit at the z height layer at the first set position x, A x The maximum normalized wind speed deficit at the hub center height of the x position is represented by H, and the relative change of the wind speed at the hub center height H is (Through the formula, the maximum normalized wind speed deficit at the hub center height of each first set position can be obtained by laser radar measurement), and x is the standard deviation of the dimensionless wind speed deficit at the first set position x, and H is the hub center height.

[0080] 5) On the basis of step 4), according to the Gaussian function satisfied by the dimensionless wind speed deficit at different height layers at each first set position The standard deviation σ of the dimensionless wind speed deficit at the first set position x is calculated by least square fitting x ; Taking the 3D at the second set position SSP as an example, at this time, according to the Gaussian function satisfied by the dimensionless wind speed deficit at different height layers at each first set position The standard deviation σ of the dimensionless wind speed deficit at the first set position x is calculated by least square fitting x , u (3D,0,H) is the axial wind speed at the hub center height H at the 3D position directly in front of the wind turbine.

[0081] Through the above manner, the standard deviation σ of the dimensionless wind speed deficit at each first set position can be accurately calculated x , in order to improve the accuracy of the inflow field speed reconstruction function reconstruction of the front blocking effect area of the wind turbine, so as to improve the accuracy of the inflow wind speed obtained at any position in the front blocking effect area of the wind turbine.

[0082] This application provides a method for obtaining the incoming wind speed of a wind turbine. Based on the Gaussian function satisfied by the dimensionless wind speed loss at different height levels at each first predetermined location, the method calculates the standard deviation of the dimensionless wind speed loss at each first predetermined location. This method may include:

[0083] according to The standard deviation of the dimensionless wind speed loss at each first set position is calculated using the least squares method.

[0084] Where, Δu (x,0,z) For the wind speed deficit at position x and height z, u (SSP,0,H) Here, D represents the axial wind speed at the hub center height level of the second designated location, where D is the rotor diameter and H is the hub center height. For the dimensionless wind speed loss at position x and height z, C T σ is the thrust coefficient of the wind turbine. x Let x be the standard deviation of the dimensionless wind speed loss at position x.

[0085] In this application, considering that lidar cannot simultaneously measure wind speeds at multiple locations and different height levels, a simple eddy current model is added to account for the maximum normalized wind speed loss at the hub center height. To Transform into (i.e., the Gaussian function satisfied by the dimensionless wind speed loss), where a is the axial induction factor, and is related to the wind turbine thrust coefficient C. T (Wind turbine thrust coefficient C) T (Originally obtained from the design curve by the wind turbine manufacturer) Then, based on the transformation obtained The standard deviation of the dimensionless wind speed loss at each of the first set locations is calculated using the least squares method. Taking the second set location SSP as 3D as an example, the simple eddy current model is as follows: After transformation, the Gaussian function satisfied by the dimensionless velocity deficit is:

[0086] The above method allows for the elimination of the need for prior procedures. Calculate A x Then, the standard deviation of the dimensionless wind speed loss is calculated, but only after determining C. T After D, the standard deviation of the dimensionless wind speed loss can be calculated using the least squares method based on the transformed formula above, thereby improving the ease of calculating the standard deviation of the dimensionless wind speed loss.

[0087] The wind turbine incoming flow wind speed acquisition method provided in the embodiment of the present application can comprise the following steps:

[0088] According to the standard deviation of the dimensionless wind speed loss of each first setting position and the functional relationship between the first setting position The least square method is used to fit the calculation function coefficients λ and η.

[0089] In the present application, the functional relationship between the standard deviation of the dimensionless wind speed loss of the first setting position and the first setting position can be obtained according to the half-wave width calculation formula in the related literature Then, the function coefficient λ and the function coefficient η in the functional relationship can be obtained by using the least square method fitting according to the standard deviation of the dimensionless wind speed loss of each first setting position and the above functional relationship. After the least square method fitting, the function coefficients λ and η are both constants, and λ is related to the rotor diameter, specifically, λ = m (1 + D) n , m and n are coefficients, and m and n can be fitted.

[0090] Through the above process, the function coefficients of the functional relationship can be accurately obtained, so as to reconstruct the incoming flow field speed reconstruction function of the front blocking effect zone of the wind turbine according to the functional relationship of the calculated function coefficients.

[0091] The wind turbine incoming flow wind speed acquisition method provided in the embodiment of the present application can comprise the following steps:

[0092] Obtaining Wherein, α is the wind shear index of the far incoming flow of the wind turbine, u (SSP,0,z) is the axial wind speed at the z height layer at the second setting position, u (SSP,0,ref) is the axial wind speed at the h (SSP,0,ref) height layer at the second setting position.

[0093] In the present application, the wind shear index of the far incoming flow of the wind turbine is obtained, and the wind shear index is fitted by the incoming flow data measured by the laser radar at the second setting position (far incoming flow) in front of the wind turbine, that is, the wind shear index formula is Wherein, u (SSP,0,z) is the axial wind speed at the z height layer at the second setting position, u (SSP,0,ref) is the axial wind speed at the reference height h (SSP,0,ref) at the second setting position, and h (SSP,0,z) represents the z height at the second setting position. Taking the second setting position SSP as 3D for example, the wind shear index formula is

[0094] By obtaining the wind shear index, the relationship between the inflow wind speed in front of the wind turbine and the inflow wind speed far away from the wind turbine (the inflow wind speed at the z height layer at the x position in front of the wind turbine is equal to the difference between the wind speed at the z height at the second set position far away from the wind turbine and the wind speed loss at the z height layer at the x position) can be transformed, so that the z height in the relationship between the inflow wind speed in front of the wind turbine and the inflow wind speed far away from the wind turbine can be transformed into a relationship related to the hub center height and the axial wind speed at the hub center position (that is, h (SSP,0,ref) = h (SSP,0,H) , u (SSP,0,ref) = u (SSP,0,H) ), so as to facilitate the convenience of the expression of the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine and the convenience of the inflow wind speed at any position in the front blocking effect zone of the wind turbine.

[0095] The wind turbine inflow wind speed acquisition method provided in the embodiments of the present application reconstructs the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine according to the function relationship of the calculated function coefficient, the wind shear index of the far inflow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed, which can include:

[0096] According to the function relationship of the calculated function coefficient, the wind shear index of the far inflow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed, the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed. The relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed is u (x,0,z) = u (SSP,0,z) - Δu (x,0,z) , and the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine is: ; wherein u (x,y,z) is the inflow wind speed at (x, y, z) in the front blocking effect zone of the wind turbine.

[0097] In the present application, when the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed according to the function relationship of the calculated function coefficient, the wind shear index of the far inflow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed, first, the relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed u (x,0,z) = u (SSP,0,z) - Δu (x,0,z) (specifically, the inflow wind speed at the x position in front of the wind turbine is the far inflow wind speed and the wind speed loss at the x position) and can be obtained Then, based on the laser radar and the fitted results, at the x position in front of the wind turbine, Reconstruction is performed on the two-dimensional yz plane, and the standard deviation function of the dimensionless wind speed loss is used to extend and reconstruct it at any x position (x < SSP) in front of the wind turbine, that is, according to and (z-H) is changed to sqrt((z-H 2 +y 2 ), so as to obtain the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine . Among them, taking the second set position SSP as an example, the 3D is obtained according to u (x,0,z) =u (3D,0,z) -Δu (x,0,z) , that is, Then, based on the laser radar and the fitting results, the Reconstruction is performed on the two-dimensional yz plane, and the standard deviation function of the dimensionless wind speed loss is used to extend and reconstruct it at any x position (x < SSP) in front of the wind turbine, so as to obtain the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine

[0098] The above-mentioned method can obtain the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine, so as to accurately calculate the inflow wind speed u (x,y,z) at (x, y, z) in the front blocking effect zone of the wind turbine according to the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine.

[0099] The embodiment of the application also provides a wind turbine inflow wind speed acquisition device, which is shown in Figure 2 which shows a structure schematic diagram of a wind turbine inflow wind speed acquisition device provided by the embodiment of the application, and can include:

[0100] The acquisition module 21 is configured to acquire the horizontal wind speed and wind direction angle of each first set position in front of the wind turbine measured by the laser radar, as the second set position of the far inflow wind speed position, and different height layers.

[0101] The first calculation module 22 is configured to calculate the standard deviation of the dimensionless wind speed loss of each first set position according to the horizontal wind speed and wind direction angle of different height layers at the second set position and each first set position.

[0102] The second calculation module 23 is configured to calculate the function coefficient according to the standard deviation of the dimensionless wind speed loss of each first set position and the function relationship between the first set position.

[0103] The obtaining module 24 is configured to obtain a flow field velocity reconstruction function of the front blockage effect zone of the wind turbine according to the function relationship of the calculated function coefficients, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow speed in front of the wind turbine and the far flow speed of the wind turbine, so as to calculate the flow speed at any position in the front blockage effect zone of the wind turbine by using the flow field velocity reconstruction function.

[0104] The wind turbine flow speed obtaining device provided in the embodiment of the application comprises a first calculation module 22, which can comprise:

[0105] The first calculation unit is configured to calculate the axial wind speed at different height layers at the second setting position by using the horizontal wind speed and the wind direction angle at different height layers at the second setting position.

[0106] The second calculation unit is configured to calculate the axial wind speed at different height layers at each first setting position by using the horizontal wind speed and the wind direction angle at different height layers at each first setting position.

[0107] The third calculation unit is configured to calculate the wind speed loss at different height layers at each first setting position by using the axial wind speed at different height layers at the second setting position and the axial wind speed at different height layers at each first setting position.

[0108] The first obtaining unit is configured to perform dimensionless processing on the wind speed loss at different height layers at each first setting position by using the axial wind speed at the hub center height layer at the second setting position, to obtain the dimensionless wind speed loss at different height layers at each first setting position.

[0109] The fourth calculation unit is configured to calculate the standard deviation of the dimensionless wind speed loss at each first setting position according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each first setting position.

[0110] The wind turbine flow speed obtaining device provided in the embodiment of the application comprises a fourth calculation unit, which can comprise:

[0111] The calculation subunit is configured to calculate the standard deviation of the dimensionless wind speed loss at each first setting position according to The standard deviation of the dimensionless wind speed loss at each first setting position is calculated by using the least square method fitting;

[0112] wherein, Δu (x,0,z) is the wind speed loss at the z height layer at the x position, u (SSP,0,H) is the axial wind speed at the hub center height layer at the second setting position, D is the diameter of the wind wheel, H is the hub center height, is the dimensionless wind speed loss at the z height layer at the x position, C T is the thrust coefficient of the wind turbine, and σ xThe standard deviation of the dimensionless wind speed deficit for the x position.

[0113] The second calculation module 23 can include:

[0114] The fifth calculation unit is configured to calculate the function coefficients λ and η according to the standard deviation of the dimensionless wind speed deficit for each first setting position and the function relationship The least square method is used to fit the calculation function coefficients λ and η.

[0115] The wind shear index acquisition module for acquiring the wind shear index of the far flow of the wind turbine can include:

[0116] The wind shear index acquisition unit is configured to acquire Wherein, α is the wind shear index of the far flow of the wind turbine, u (SSP,0,z) is the axial wind speed at the z height layer at the second setting position, u (SSP,0,ref) is the axial wind speed at the h (SSP,0,ref) height layer at the second setting position.

[0117] The obtaining module 24 can include:

[0118] The obtaining unit is configured to obtain and the relationship between the flow speed in front of the wind turbine and the far flow speed u (x,0,z) = u (SSP,0,z) - Δu (x,0,z) The flow field speed reconstruction function of the front blocking effect area of the wind turbine is obtained as follows: ; wherein, u (x,y,z) is the flow speed at (x, y, z) in the front blocking effect area of the wind turbine.

[0119] The embodiment of the present application further provides a wind turbine flow speed acquisition device, referring to Figure 3 which shows a structural schematic diagram of a wind turbine flow speed acquisition device provided by the embodiment of the present application, and can include:

[0120] The memory is configured to store the computer program.

[0121] The processor is configured to realize the following steps when executing the computer program stored in the memory.

[0122] The horizontal wind speed and the wind direction angle of each first setting position in front of the wind turbine and a second setting position as a far flow wind speed position, and different height layers are acquired by a laser radar; the standard deviation of the dimensionless wind speed loss of each first setting position is calculated according to the horizontal wind speed and the wind direction angle of the second setting position and different height layers of each first setting position; the function coefficient is calculated according to the standard deviation of the dimensionless wind speed loss of each first setting position and a functional relationship between the first setting position; the flow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed according to the functional relationship of the calculated function coefficient, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow wind speed in front of the wind turbine and the far flow wind speed, so as to calculate the flow wind speed at any position in the front blocking effect zone of the wind turbine by using the flow field speed reconstruction function.

[0123] The embodiment of the present application further provides a readable storage medium, and the readable storage medium stores a computer program.

[0124] The horizontal wind speed and the wind direction angle of each first setting position in front of the wind turbine and a second setting position as a far flow wind speed position, and different height layers are acquired by a laser radar; the standard deviation of the dimensionless wind speed loss of each first setting position is calculated according to the horizontal wind speed and the wind direction angle of the second setting position and different height layers of each first setting position; the function coefficient is calculated according to the standard deviation of the dimensionless wind speed loss of each first setting position and a functional relationship between the first setting position; the flow field speed reconstruction function of the front blocking effect zone of the wind turbine is reconstructed according to the functional relationship of the calculated function coefficient, the wind shear exponent of the far flow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and the relationship between the flow wind speed in front of the wind turbine and the far flow wind speed, so as to calculate the flow wind speed at any position in the front blocking effect zone of the wind turbine by using the flow field speed reconstruction function.

[0125] The readable storage medium can include a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various media that can store program codes.

[0126] The description of the related parts in the wind turbine flow wind speed acquisition device, the equipment and the readable storage medium provided by the present application can refer to the detailed description of the corresponding parts in the wind turbine flow wind speed acquisition method provided by the embodiment of the present application, and will not be repeated here.

[0127] It is to be understood that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting; it is not intended to exclude myriad other embodiments of the present application that other present or future technologies can enable. One of ordinary skill in the art will recognize that the functions explained herein can be combined in alternate implementations or partitioned differently without losing the intended benefits associated with the present application. One of ordinary skill in the art will also recognize that a plurality of hardware and software based devices and / or components can be used in carrying out the teachings of the present application. In these respects, the application specification is to be considered only as illustrative and not restrictive.

[0128] The above description of disclosed embodiments provides enough information to enable those with skill in the art to make and use the application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein can be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for obtaining the wind speed of the incoming flow of a wind turbine, characterized in that, The method comprises the following steps: acquiring horizontal wind speed and wind direction angle at each first set position in front of a wind turbine and at a second set position as a far flow wind speed position measured by a laser radar; calculating standard deviation of dimensionless wind speed loss at each first set position according to horizontal wind speed and wind direction angle at different height layers at the second set position and at each first set position; calculating a function coefficient according to standard deviation of dimensionless wind speed loss at each first set position and a function relationship between the first set position; reconstructing a flow field speed reconstruction function of a front blockage effect zone of the wind turbine according to a function relationship of the calculated function coefficient, a wind shear exponent of far flow of the wind turbine, a Gaussian function satisfied by dimensionless wind speed loss, and a relationship between flow wind speed in front of the wind turbine and far flow wind speed, so as to calculate flow wind speed at any position in the front blockage effect zone of the wind turbine by using the flow field speed reconstruction function; wherein the calculating standard deviation of dimensionless wind speed loss at each first set position according to horizontal wind speed and wind direction angle at different height layers at the second set position and at each first set position comprises: calculating axial wind speed at different height layers at the second set position by using horizontal wind speed and wind direction angle at different height layers at the second set position; correspondingly calculating axial wind speed at different height layers at each first set position by using horizontal wind speed and wind direction angle at different height layers at each first set position; correspondingly calculating wind speed loss at different height layers at each first set position by using axial wind speed at different height layers at the second set position and axial wind speed at different height layers at each first set position; performing dimensionless processing on wind speed loss at different height layers at each first set position by using axial wind speed at a hub center height layer at the second set position, to obtain dimensionless wind speed loss at different height layers at each first set position; calculating standard deviation of dimensionless wind speed loss at each first set position according to a Gaussian function satisfied by dimensionless wind speed loss at different height layers at each first set position.

2. The method of claim 1, wherein The calculating standard deviation of dimensionless wind speed loss at each first set position according to a Gaussian function satisfied by dimensionless wind speed loss at different height layers at each first set position comprises: According to , the standard deviation of the dimensionless wind speed deficit for each of the first set positions is calculated using a least squares fit; wherein is the wind speed deficit at the z-height level at the x-position, is the axial wind speed at the hub center height level at the second set position, D is the rotor diameter, and H is the hub center height, is the dimensionless wind speed deficit at the z-height level at the x-position, is the thrust coefficient of the wind turbine, is the standard deviation of the dimensionless wind speed deficit at the x-position.

3. The method of claim 2, wherein, The calculating a function coefficient according to standard deviation of dimensionless wind speed loss at each first set position and a function relationship between the first set position comprises: standard deviation of the dimensionless wind speed deficit and the functional relationship according to the respective first set position , the function coefficients are calculated by least square fitting and .

4. The method of claim 2, wherein acquiring a wind shear exponent of far flow of the wind turbine comprises: acquiring ; wherein is the wind shear exponent of the wind coming from a distance to the wind turbine, is the axial wind speed at the height layer z at the second set position, is the axial wind speed at the height layer z at the second set position, is the axial wind speed at the height layer z at the second set position, represents the height at the second set position, is the height at the second set position.

5. The method of wind speed acquisition from the incoming flow of a wind turbine according to claim 4, characterized in that, The reconstructing a flow field speed reconstruction function of a front blockage effect zone of the wind turbine according to a function relationship of the calculated function coefficient, a wind shear exponent of far flow of the wind turbine, a Gaussian function satisfied by dimensionless wind speed loss, and a relationship between flow wind speed in front of the wind turbine and far flow wind speed comprises: According to , , and the relationship between the inflow wind speed in front of the wind turbine and the far-field inflow wind speed , the inflow field speed reconstruction function of the front blocking effect zone of the wind turbine is obtained as: ; wherein, is the inflow wind speed at (x, y, z) in the front blocking effect zone of the wind turbine, is the inflow wind speed at the x position in front of the wind turbine.

6. A wind turbine incoming flow wind speed acquisition device, characterized in that, The method comprises the following steps: an acquisition module is configured to acquire horizontal wind speed and wind direction angle at each first set position in front of a wind turbine and at a second set position as a far flow wind speed position measured by a laser radar; The first calculation module is configured to calculate the standard deviation of the dimensionless wind speed loss of each first setting position according to the horizontal wind speed and the wind direction angle of different height layers at the second setting position and at each first setting position; The second calculation module is configured to calculate the function coefficient according to the standard deviation of the dimensionless wind speed loss of each first setting position and a function relationship between the first setting position and the function coefficient; The obtaining module is configured to obtain the inflow field speed reconstruction function of the front blockage effect zone of the wind turbine according to the function relationship of the calculated function coefficient, the wind shear exponent of the far inflow of the wind turbine, the Gaussian function satisfied by the dimensionless wind speed loss, and a relationship between the inflow wind speed in front of the wind turbine and the far inflow wind speed, so as to calculate the inflow wind speed at any position in the front blockage effect zone of the wind turbine by using the inflow field speed reconstruction function. The first calculation module comprises: The first calculation unit is configured to calculate the axial wind speed at different height layers at the second setting position by using the horizontal wind speed and the wind direction angle of different height layers at the second setting position; The second calculation unit is configured to calculate the axial wind speed at different height layers at each first setting position by using the horizontal wind speed and the wind direction angle of different height layers at each first setting position; The third calculation unit is configured to calculate the wind speed loss at different height layers at each first setting position by using the axial wind speed at different height layers at the second setting position and the axial wind speed at different height layers at each first setting position; The first obtaining unit is configured to perform dimensionless processing on the wind speed loss at different height layers at each first setting position by using the axial wind speed at the hub center height layer at the second setting position, to obtain the dimensionless wind speed loss at different height layers at each first setting position; The fourth calculation unit is configured to calculate the standard deviation of the dimensionless wind speed loss of each first setting position according to the Gaussian function satisfied by the dimensionless wind speed loss at different height layers at each first setting position.

7. A wind turbine incoming flow wind speed acquisition device, characterized in that The memory is configured to store a computer program; The processor is configured to implement the steps of the wind turbine inflow wind speed acquisition method according to any one of claims 1 to 5 when executing the computer program. The readable storage medium stores a computer program, and the computer program is executed by the processor to implement the steps of the wind turbine inflow wind speed acquisition method according to any one of claims 1 to 5.

8. A readable storage medium, characterized by, ​

Citation Information

Patent Citations

  • Method for calculating wind turbine wake flow by using two-dimensional Frandsen wake flow model

    CN112906321A

  • Three-dimensional wake flow wind speed distribution calculation method suitable for various terrains

    CN114398843A