A method, system and device for determining wind conditions for wind turbines

Abstract

The present invention discloses a method, system, and device for determining wind conditions for a wind turbine generator set, which are applied to the field of wind power technology. The method comprises setting a plurality of nodes on the plane where the wind rotor of the wind turbine generator set is located; determining the actual wind speed history of each node, where the actual wind speed history is the corresponding relationship between the actual wind speed and time; determining the average wind speed distribution of each node based on the actual wind speed history of each node; and determining the simulated wind speed history of each node based on the actual wind speed history. The turbulent wind speed history represented by the actual wind speed history is relatively accurate, the average wind speed distribution represents the dynamic response characteristics under stable wind conditions at a complex site, and the turbulent wind speed history represented by the simulated wind speed history can be quickly generated in batches. Analyzing the dynamic response of a wind turbine generator set based on the actual wind speed history, the average wind speed distribution, or the simulated wind speed history is beneficial for guiding turbine design, improving turbine operation and maintenance, and ensuring turbine structural safety.

Description

Technical Field

[0001] The present invention relates to the field of wind power technology, and in particular to a method, system and device for determining wind conditions applicable to a wind turbine generator set. Background Art

[0002] In related technologies, wind turbine characteristics are obtained and then controlled to generate electricity based on these characteristics. However, wind conditions generated based on standards are difficult to reflect the wind characteristics of complex sites, resulting in units that have passed standard evaluations being difficult to adapt to complex sites, and prone to load overload, clearance risks, frequent shutdowns, and reduced power generation. In complex site environments such as high slopes, peaks, and valleys, there are significant differences or mutations in the wind speed spectrum in the horizontal direction or locally, which can easily lead to dynamic response characteristics of the unit that are beyond the standard evaluation, causing operation and maintenance accidents. Summary of the Invention

[0003] The object of the present invention is to provide a method, system and device for determining wind conditions of a wind turbine generator system, wherein the obtained data is more accurate and the determined dynamic response is more in line with the requirements.

[0004] To solve the above technical problems, the present invention provides a method for determining wind conditions for a wind turbine generator system, comprising:

[0005] A plurality of nodes are set on the plane where the wind rotor of the wind turbine is located;

[0006] Determining an actual wind speed time series for each of the nodes, wherein the actual wind speed time series is a correspondence between actual wind speed and time;

[0007] Determine the average wind speed distribution of each node according to the actual wind speed time history of each node, wherein the average wind speed distribution is the average value of the actual wind speed within a preset time;

[0008] Determining a simulated wind speed time history for each of the nodes based on the actual wind speed time history, wherein the simulated wind speed time history represents multiple wind speeds for each of the nodes;

[0009] The dynamic response of the wind turbine is determined according to the actual wind speed time history, the average wind speed distribution or the simulated wind speed time history, wherein the dynamic response includes the power generation of the wind turbine, component loads and the distance between the blades and the tower of the wind turbine.

[0010] On the other hand, multiple nodes are set on the plane where the wind turbine rotor is located, including:

[0011] The plane where the wind turbine rotor is located is divided into preset grids, and the vertices of each grid are used as the nodes.

[0012] On the other hand, dividing the grid on the plane where the wind rotor of the wind turbine generator set is located includes:

[0013] The number and area of ​​the grids are determined so that the locations where the wind speed changes on the plane are aligned with the nodes.

[0014] On the other hand, determining the actual wind speed time history of each of the nodes includes:

[0015] The actual wind speed time course V of each node is determined by wind measuring radar experiment (y j ,z j ,t,i);

[0016] Among them, y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0017] On the other hand, determining the actual wind speed time history of each of the nodes includes:

[0018] Setting the wind direction and plane of the plane, and establishing a fluid simulation model of the plane where the wind wheel is located;

[0019] Determine the actual wind speed time history V of each node according to the fluid simulation model experiment (y j ,z j ,t,i), where y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0020] On the other hand, determining the average wind speed distribution according to the actual wind speed time history of each of the nodes includes:

[0021] Determining an average wind speed of each of the nodes within a preset time period according to the wind speed time history;

[0022] Determine the average wind speed distribution matrix based on the average wind speed and the nodes The number of nodes in the horizontal direction is k, and the number of nodes in the vertical direction is h. is the average wind speed of the node at the hth row and the kth column in the i direction within the preset time, where i=1, 2, and 3 represent the longitudinal, horizontal, and vertical directions respectively.

[0023] On the other hand, determining the simulated wind speed time history of each of the nodes according to the actual wind speed time history includes:

[0024] The wind speeds of the actual wind speed time history matrix are converted from the time domain to the frequency domain to obtain the wind spectrum characteristic distribution matrix Among them S hk (w) is the wind power spectrum density value in direction i at the node in row h and column k within the preset time;

[0025] Determining the autospectral value of the wind speed spectrum density matrix according to the wind spectrum characteristic distribution matrix, and determining the cross-spectral value of the wind speed spectrum density matrix according to the spatial coherence function;

[0026] The wind speed spectral density matrix is ​​converted from the frequency domain to the time domain to obtain the simulated wind speed time history.

[0027] On the other hand, converting the wind speed spectral density matrix from the frequency domain to the time domain to obtain the simulated wind speed time history includes:

[0028] The wind speed spectrum density matrix S i (w) According to the first relation to decompose;

[0029] According to the second relationship, the wind speed spectrum density matrix is ​​converted from the frequency domain to the time domain to obtain the simulated wind speed time history

[0030] Among them, i=1, 2, 3 represent the longitudinal, horizontal, and vertical directions respectively, y j and z j represent the horizontal and vertical coordinates of the j-th point, respectively. represents the frequency increment, w up is the cutoff frequency, N is an integer greater than 1; Φ ml is a random phase angle, uniformly distributed in the interval [0, 2π]; θ jm (w ml ) is H jm The complex angle of (w), H jm (w) is the value of the jth row and mth column in the matrix Hi(w),

[0031] To solve the above technical problems, the present invention further provides a wind condition determination system applicable to a wind turbine generator system, comprising:

[0032] A setting unit, used for setting a plurality of nodes on the plane where the wind rotor of the wind turbine generator set is located;

[0033] a first determining unit, configured to determine an actual wind speed time series of each of the nodes, wherein the actual wind speed time series is a correspondence between actual wind speed and time;

[0034] A second determining unit is configured to determine an average wind speed distribution of each of the nodes based on the actual wind speed time history of each of the nodes, wherein the average wind speed distribution is an average value of the actual wind speed within a preset time;

[0035] a third determining unit, configured to determine a simulated wind speed time history of each of the nodes according to the average wind speed distribution, wherein the simulated wind speed time history represents a plurality of wind speeds of each of the nodes;

[0036] The fourth determining unit is configured to determine a dynamic response of the wind turbine generator set according to the actual wind speed time series, the average wind speed distribution or the simulated wind speed time series, wherein the dynamic response includes power generation of the wind turbine generator set and positions of components.

[0037] To solve the above technical problems, the present invention further provides a wind condition determination device applicable to a wind turbine generator system, comprising:

[0038] Memory for storing computer programs;

[0039] The processor is configured to implement the steps of the above-mentioned method for determining wind conditions applicable to a wind turbine generator system when executing the computer program.

[0040] The present invention discloses a method, system, and device for determining wind conditions applicable to a wind turbine generator set, which are applied to the field of wind power technology. The method comprises setting a plurality of nodes on the plane where the wind rotor of the wind turbine generator set is located; determining the actual wind speed history of each node, where the actual wind speed history is the corresponding relationship between the actual wind speed and time; determining the average wind speed distribution of each node based on the actual wind speed history of each node; and determining the simulated wind speed history of each node based on the actual wind speed history. The turbulent wind speed history represented by the actual wind speed history is relatively accurate, the average wind speed distribution represents the dynamic response characteristics under stable wind conditions at a complex site, and the turbulent wind speed history represented by the simulated wind speed history can be quickly generated in batches. Analyzing the dynamic response of a wind turbine generator set based on the actual wind speed history, the average wind speed distribution, or the simulated wind speed history is beneficial for guiding turbine design, improving turbine operation and maintenance, and ensuring turbine structural safety. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0042] Figure 1 A flow chart of a method for determining wind conditions applicable to a wind turbine generator system provided by the present invention;

[0043] Figure 2A schematic diagram of a wind rotor plane grid division provided by the present invention;

[0044] Figure 3 A schematic diagram of a typical average wind speed distribution at a wind rotor plane grid point provided by the present invention;

[0045] Figure 4 A three-dimensional wind speed time history diagram of a typical grid point on a wind rotor plane provided by the present invention;

[0046] Figure 5 A typical dynamic parameter time history diagram under complex site wind conditions provided by the present invention;

[0047] Figure 6 A schematic structural diagram of a wind condition determination system applicable to a wind turbine generator system provided by the present invention;

[0048] Figure 7 This is a structural schematic diagram of a wind condition determination device applicable to a wind turbine generator system provided by the present invention. DETAILED DESCRIPTION

[0049] The core of the present invention is to provide a method, system and device for determining wind conditions of wind turbines, so that the obtained data is more accurate and the determined dynamic response is more in line with the requirements.

[0050] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts shall fall within the scope of protection of the present invention.

[0051] Figure 1 The present invention provides a flow chart of a method for determining wind conditions for a wind turbine generator system, including:

[0052] S11: setting a plurality of nodes on the plane where the wind turbine rotor is located;

[0053] Taking into account that the wind rotor of the wind turbine is constantly rotating and the sizes of different parts of the wind rotor are different, the wind speed at each position of the plane where the wind rotor is located is different. In order to improve the accuracy of wind speed measurement, multiple nodes are set on the plane where the wind rotor is located. By determining the wind speed of each node, the wind speed of the plane where the wind rotor is located can be determined.

[0054] Specifically, by adjusting the positions of the nodes, the wind speed at each position on the plane where the wind rotor is located can be obtained more accurately.

[0055] S12: Determine the actual wind speed time history of each node, where the actual wind speed time history is the corresponding relationship between the actual wind speed and time;

[0056] Determine the actual wind speed of each node and obtain the corresponding relationship between the actual wind speed and time of each node within the preset time. The actual wind speed time series can be a matrix or other expression method, and this application does not make too many restrictions here.

[0057] Based on the actual wind speed, the actual wind speed can be further processed to obtain more accurate wind conditions.

[0058] Specifically, the method for determining the actual wind speed is not limited in this application.

[0059] S13: Determine the average wind speed distribution of each node based on the actual wind speed time history of each node, where the average wind speed distribution is the average value of the actual wind speed within a preset time;

[0060] After determining the actual wind speed time series, the actual wind speed time series within the preset time can be determined and its average value can be calculated. The calculated wind speed can represent the wind speed within the preset time. If the wind speed is relatively stable, the average wind speed can evaluate the dynamic response characteristics of the wind turbine under stable wind conditions at a complex site.

[0061] S14: determining a simulated wind speed history for each node according to the actual wind speed history, where the simulated wind speed history represents multiple wind speeds for each node;

[0062] Taking into account that under the condition of meeting the specified frequency, there may be many wind speeds at the same node, and the actual wind speed detected belongs to only one of them, the simulated wind speed time series determined according to the actual wind speed time series can determine multiple wind speeds at the same node, which makes it faster to calculate the dynamic response of the wind turbine.

[0063] S15: Determine a dynamic response of the wind turbine generator set according to the actual wind speed time history, the average wind speed distribution or the simulated wind speed time history, where the dynamic response includes power generation of the wind turbine generator set, component loads and a distance between blades and a tower of the wind turbine generator set.

[0064] The above-mentioned actual wind speed time series, average wind speed distribution and simulated wind speed time series all have their own advantages. The turbulent wind speed time series of the actual wind speed time series is relatively accurate. The average wind speed distribution represents the dynamic response characteristics under stable wind conditions at complex sites. The turbulent wind speed time series represented by the simulated wind speed time series can be quickly generated in batches. The three wind speeds all have their own advantages. One of the three wind speeds can be selected as the basis for calculating the dynamic response according to actual needs.

[0065] The present invention discloses a method, system, and device for determining wind conditions applicable to a wind turbine generator set, which are applied to the field of wind power technology. The method comprises setting a plurality of nodes on the plane where the wind rotor of the wind turbine generator set is located; determining the actual wind speed history of each node, where the actual wind speed history is the corresponding relationship between the actual wind speed and time; determining the average wind speed distribution of each node based on the actual wind speed history of each node; and determining the simulated wind speed history of each node based on the actual wind speed history. The turbulent wind speed history represented by the actual wind speed history is relatively accurate, the average wind speed distribution represents the dynamic response characteristics under stable wind conditions at a complex site, and the turbulent wind speed history represented by the simulated wind speed history can be quickly generated in batches. Analyzing the dynamic response of a wind turbine generator set based on the actual wind speed history, the average wind speed distribution, or the simulated wind speed history is beneficial for guiding turbine design, improving turbine operation and maintenance, and ensuring turbine structural safety.

[0066] Based on the above embodiment:

[0067] Figure 2 A schematic diagram of a wind rotor plane grid division provided by the present invention;

[0068] In some embodiments, a plurality of nodes are set on the plane where the wind turbine rotor is located, including:

[0069] The plane where the wind turbine rotor is located is divided into preset grids, and the vertices of each grid are used as nodes.

[0070] In some embodiments, dividing a grid on a plane where a wind turbine rotor is located includes:

[0071] Determine the number and area of ​​the grids so that the wind speed changes on the plane are aligned with the nodes.

[0072] Assume that the number of horizontal grid points is k and the number of vertical grid points is h. If the unit size is too large and the number of grid nodes is too small, the wind speed distribution characteristics of the wind turbine surface cannot be effectively captured; if the unit size is too small and the number of grid nodes is too large, the implementation cost will increase significantly.

[0073] For areas where the wind speed on the rotor plane changes dramatically, it is necessary to align the grid nodes or appropriately encrypt the grid to obtain higher capture accuracy to reflect the peak and valley changes in wind speed.

[0074] Specifically, based on the current mainstream wind rotor size, the unit size is generally 5 to 10 meters, and the number of grid nodes is generally 20 to 30.

[0075] Figure 3 A schematic diagram of a typical average wind speed distribution at a wind rotor plane grid point provided by the present invention;

[0076] like Figure 3As shown in the figure, the average wind speed distribution at the rotor surface for the complex site wind condition obtained by fluid simulation shows that the average wind speed increases with height, and the wind speed gradually increases from the -Y axis to the +Y axis. Therefore, the wind speed is the highest at the upper right corner of the rotor and the lowest at the lower left corner of the rotor.

[0077] In some embodiments, determining the actual wind speed time history of each node includes:

[0078] The actual wind speed time course V of each node is determined by wind radar experiment (y j ,z j ,t,i);

[0079] Among them, y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0080] In some embodiments, determining the actual wind speed time history of each node includes:

[0081] Set the wind direction and plane of the plane, and build a fluid simulation model of the plane where the wind wheel is located;

[0082] Determine the actual wind speed time history V at each node based on the fluid simulation model experiment (y j ,z j ,t,i), where y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0083] Complex site environments are typically found in areas with significant terrain variations, such as high slopes, valleys, and mountain peaks. Wind characteristics such as the spatial average wind speed distribution, turbulence intensity, and wind power spectrum at these sites differ significantly from regulatory standards, making it impossible to generate corresponding wind speed time histories based on these standards. Therefore, these characteristics must be determined through field measurements or simulations. Wind radars can scan the wind speed distribution along the wind turbine rotor plane in real time. Large grid sizes require the installation of two or more radars operating simultaneously. By establishing a three-dimensional model of the complex site and setting the inflow and boundary conditions, computational fluid dynamics (CFD) time-domain simulations can also be used to determine the wind speed time histories at each grid point.

[0084] Specific acquisition methods include but are not limited to the above two methods, and this application does not make too many restrictions here.

[0085] In some embodiments, determining the average wind speed distribution based on the actual wind speed time history of each node includes:

[0086] Determine the average wind speed of each node within a preset time based on the wind speed time history;

[0087] Determine the average wind speed distribution matrix based on the average wind speed and nodes The number of nodes in the horizontal direction is k, and the number of nodes in the vertical direction is h. is the average wind speed of the node at the hth row and the kth column in the i direction within the preset time, where i=1, 2, and 3 represent the longitudinal, horizontal, and vertical directions respectively.

[0088] The wind direction is divided into three directions: longitudinal, horizontal and vertical. The wind speed in each direction is calculated in turn. When i=1, the longitudinal wind speed of each node is determined, and the average longitudinal wind speed within the preset time is calculated. When i=2, the horizontal wind speed of each node is determined, and the average horizontal wind speed within the preset time is calculated. When i=3, the vertical wind speed of each node is determined, and the average vertical wind speed within the preset time is calculated. Finally, three average wind speed distribution matrices are formed: and It reflects the wind speed distribution characteristics of the wind wheel surface in a complex site environment.

[0089] Figure 4 A three-dimensional wind speed time history diagram of a typical grid point on a wind rotor plane provided by the present invention;

[0090] The wind speed time histories in the longitudinal (u), horizontal (v) and vertical (w) directions of typical points are given.

[0091] Figure 5 A typical dynamic parameter time history diagram under complex site wind conditions provided by the present invention;

[0092] The time history diagrams of typical parameters such as pitch angle, blade tip out-of-plane deformation, tower base My torque, and generator power of the MW unit under complex site wind conditions are given.

[0093] In some embodiments, determining a simulated wind speed time history for each node based on the actual wind speed time history includes:

[0094] Convert each wind speed of the actual wind speed time history matrix from the time domain to the frequency domain to obtain the wind spectrum characteristic distribution matrix Among them S hk (w) is the wind power spectrum density value in direction i at the node in row h and column k within the preset time;

[0095] The autospectral value of the wind speed spectrum density matrix is ​​determined according to the wind spectrum characteristic distribution matrix, and the cross-spectral value of the wind speed spectrum density matrix is ​​determined according to the spatial coherence function;

[0096] The wind speed spectral density matrix is ​​converted from the frequency domain to the time domain to obtain the simulated wind speed time history.

[0097] Based on the fast Fourier transform, the wind speed time history data of each point on the wind rotor plane is converted from the time domain to the frequency domain, and the wind spectrum characteristic distribution matrix of the wind rotor surface is obtained, that is, Convert to

[0098] In order to determine the wind speed spectral density matrix S i (w), S i (w) contains the self-spectral value and the cross-spectral value, and the wind spectrum characteristic distribution matrix can be directly used Alternatively, the autospectral values ​​at different grid points can be calculated using a typical wind spectrum. Typical wind spectra can include Kaimal, VonKarman, or Mann spectra. Cross-spectral values ​​between different grid points can be calculated based on a spatial coherence function, such as the Davenport function.

[0099] In addition, the simulated wind speed time history can be obtained by converting the standard wind spectrum formula from the frequency domain to the time domain.

[0100] In some embodiments, converting the wind speed spectral density matrix from the frequency domain to the time domain to obtain a simulated wind speed time history includes:

[0101] The wind speed spectral density matrix S i (w) According to the first relation to decompose;

[0102] According to the second relationship, the wind speed spectrum density matrix is ​​converted from the frequency domain to the time domain to obtain the simulated wind speed time history.

[0103] Among them, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively, y j and z j represent the horizontal and vertical coordinates of the j-th point, respectively. represents the frequency increment, w up is the cutoff frequency, N is an integer greater than 1; Φ ml is a random phase angle, uniformly distributed in the interval [0, 2π]; θ jm (w ml ) is H jm The complex angle of (w), H jm (w) is the value of the jth row and mth column in the matrix Hi(w),

[0104] Due to the existence of random phase angle Φ ml The calculated simulated wind speed time series can be used to obtain different wind speeds by adjusting the random phase angle, thereby achieving batch acquisition of the wind speed of the node in multiple directions.

[0105] The three wind speed matrices are written into a binary format wind file, which is then called by the wind turbine simulation software to simulate the wind turbine dynamic response under complex site wind conditions and guide the turbine design. The three wind speed matrices written are:

[0106] The first is the average wind speed distribution matrix It can simulate and evaluate the dynamic response characteristics of the unit under stable wind conditions at complex sites;

[0107] The second is the actual wind speed time series V obtained by fluid simulation or wind radar. experiment (y j ,z j ,t,i), the turbulent wind speed time history is relatively accurate;

[0108] The third is the simulated wind speed time history matrix V obtained based on wind spectrum simulation. simulation (y j ,z j ,t,i), the turbulent wind speed time history can be quickly generated in batches.

[0109] By evaluating and optimizing the overall machine control response, load safety, clearance risk, and electrical performance under complex site wind conditions, the adaptability and operation and maintenance safety of wind turbines in complex site environments can be improved.

[0110] Figure 6 A schematic diagram of a wind condition determination system for a wind turbine generator system provided by the present invention includes:

[0111] A setting unit 21 is used to set a plurality of nodes on the plane where the wind rotor of the wind turbine is located;

[0112] A first determining unit 22 is configured to determine an actual wind speed time series of each node, where the actual wind speed time series is a correspondence between actual wind speed and time;

[0113] A second determining unit 23 is configured to determine an average wind speed distribution of each node based on the actual wind speed time history of each node, where the average wind speed distribution is an average value of the actual wind speed within a preset time period;

[0114] A third determining unit 24 is configured to determine a simulated wind speed history for each node based on the actual wind speed history, wherein the simulated wind speed history represents multiple wind speeds for each node;

[0115] The fourth determination unit 25 is used to determine the dynamic response of the wind turbine according to the actual wind speed time history, the average wind speed distribution or the simulated wind speed time history, the dynamic response including the power generation of the wind turbine, the component load and the distance between the blades and the tower of the wind turbine.

[0116] The setting unit 21 is specifically configured to divide the plane where the wind turbine rotor is located into preset grids, and use the vertices of each grid as nodes.

[0117] The setting unit 21 is specifically used to determine the number and area of ​​the grids so that the wind speed changes on the plane are aligned with the nodes, and the vertices of each grid are used as nodes.

[0118] The first determining unit 22 is specifically configured to determine the actual wind speed time series V of each node through a wind measuring radar. experiment (y j ,z j ,t,i);

[0119] Among them, y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0120] A setting unit is used to set the wind direction and plane of the plane and establish a fluid simulation model of the plane where the wind wheel is located;

[0121] The first determining unit 22 is specifically configured to determine the actual wind speed time series V of each node according to the fluid simulation model. experiment (y j ,z j ,t,i), where y j and z j Represent the horizontal and vertical coordinates of the jth point respectively, j = 1, 2…, k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is time, i = 1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

[0122] a fifth determining unit, configured to determine an average wind speed of each node within a preset time period according to the wind speed time history;

[0123] The second determining unit 23 is specifically configured to determine the average wind speed distribution matrix according to the average wind speed and the nodes. The number of nodes in the horizontal direction is k, and the number of nodes in the vertical direction is h. is the average wind speed of the node at the hth row and the kth column in the i direction within the preset time, where i=1, 2, and 3 represent the longitudinal, horizontal, and vertical directions respectively.

[0124] The first conversion unit is used to convert each wind speed of the actual wind speed time series from the time domain to the frequency domain to obtain the wind spectrum characteristic distribution matrix Among them S hk (w) is the wind power spectrum density value in direction i at the node in row h and column k within the preset time;

[0125] a sixth determining unit, configured to determine an autospectral value of a wind speed spectrum density matrix according to the wind spectrum characteristic distribution matrix, and to determine a cross-spectral value of the wind speed spectrum density matrix according to a spatial coherence function;

[0126] The third determining unit 24 is specifically configured to convert the wind speed spectrum density matrix from the frequency domain to the time domain to obtain a simulated wind speed time history.

[0127] Decomposition unit, used to convert the wind speed spectrum density matrix S i (w) According to the first relation to decompose;

[0128] The third determining unit 24 is specifically configured to convert the wind speed spectrum density matrix from the frequency domain to the time domain according to the second relationship to obtain a simulated wind speed time history.

[0129] Among them, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively, y j and z j represent the horizontal and vertical coordinates of the j-th point, respectively. represents the frequency increment, w up is the cutoff frequency, N is an integer greater than 1; Φ ml is a random phase angle, uniformly distributed in the interval [0, 2π]; θ jm (w ml ) is H jm The complex angle of (w), H jm (w) is the value of the jth row and mth column in the matrix Hi(w),

[0130] Figure 7 This is a schematic structural diagram of a wind condition determination device for a wind turbine generator system provided by the present invention, the device comprising:

[0131] Memory 31, for storing computer programs;

[0132] The processor 32 is configured to implement the steps of the above-mentioned method for determining wind conditions applicable to a wind turbine when executing a computer program.

[0133] For an introduction to the wind condition determination device applicable to a wind turbine generator system provided in this application, please refer to the above embodiments and will not be repeated here.

[0134] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. Reference can be made to the common and similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the description is relatively simple, and the relevant parts can be referred to the method description.

[0135] It should also be noted that, in this specification, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variants thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus comprising the element.

[0136] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0137] The above description of the disclosed embodiments is intended to enable one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention is not limited to the embodiments shown herein but is intended to conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for determining wind conditions for a wind turbine generator system, characterized in that: include: A plurality of nodes are set on the plane where the wind rotor of the wind turbine is located; Determining an actual wind speed time series for each of the nodes, wherein the actual wind speed time series is a correspondence between actual wind speed and time; Determine the average wind speed distribution of each node according to the actual wind speed time history of each node, wherein the average wind speed distribution is the average value of the actual wind speed within a preset time; Determining a simulated wind speed time history for each of the nodes based on the actual wind speed time history, wherein the simulated wind speed time history represents multiple wind speeds for each of the nodes; determining a dynamic response of the wind turbine generator set according to the actual wind speed time history, the average wind speed distribution, or the simulated wind speed time history, wherein the dynamic response includes power generation of the wind turbine generator set, component loads, and a distance between blades and a tower of the wind turbine generator set; Determining a simulated wind speed time history of each of the nodes according to the actual wind speed time history includes: The wind speeds of the actual wind speed time series are converted from the time domain to the frequency domain to obtain the wind spectrum characteristic distribution matrix ,in is the wind power spectrum density value in direction i at the node in row h and column k within the preset time; Determining the autospectral value of the wind speed spectrum density matrix according to the wind spectrum characteristic distribution matrix, and determining the cross-spectral value of the wind speed spectrum density matrix according to the spatial coherence function; Converting the wind speed spectral density matrix from the frequency domain to the time domain to obtain the simulated wind speed time history; Converting the wind speed spectral density matrix from the frequency domain to the time domain to obtain the simulated wind speed time history includes: The wind speed spectral density matrix According to the first relation to decompose; According to the second relationship, the wind speed spectrum density matrix is ​​converted from the frequency domain to the time domain to obtain the simulated wind speed time history ; Among them, t is time, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively, y j and z j represent the horizontal and vertical coordinates of the j-th point, respectively. represents the frequency increment, w up is the cutoff circular frequency, N is an integer greater than 1; is a random phase angle in the interval Evenly distributed within; for The complex angle, is the matrix H i The value of row j and column m in (w) is , the number of nodes in the vertical direction is h.

2. The method for determining wind conditions for a wind turbine generator system according to claim 1, wherein: Multiple nodes are set on the plane where the wind turbine rotor is located, including: The plane where the wind turbine rotor is located is divided into a preset number of grids, and the vertices of each grid are used as the nodes.

3. The method for determining wind conditions for a wind turbine generator system according to claim 2, wherein: Dividing a grid on a plane where a wind rotor of the wind turbine generator set is located includes: The number and area of ​​the grids are determined so that the locations where the wind speed changes on the plane are aligned with the nodes.

4. The method for determining wind conditions for a wind turbine generator system according to claim 1, wherein: Determining the actual wind speed time history of each of the nodes, including: Determine the actual wind speed time history of each node by wind measuring radar ; Among them, y j and z j Represent the horizontal and vertical coordinates of the jth point, j=1, 2… ., k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is the time, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

5. The method for determining wind conditions for a wind turbine generator system according to claim 1, wherein: Determining the actual wind speed time history of each of the nodes, including: Setting the wind direction and plane of the plane, and establishing a fluid simulation model of the plane where the wind wheel is located; Determine the actual wind speed time history of each node according to the fluid simulation model , where y j and z j Represent the horizontal and vertical coordinates of the jth point, j=1, 2…., k*h, k is the number of nodes in the horizontal direction, h is the number of nodes in the vertical direction, t is the time, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively.

6. The method for determining wind conditions for a wind turbine generator system according to claim 1, wherein: Determining the average wind speed distribution according to the actual wind speed time history of each of the nodes includes: Determining an average wind speed of each of the nodes within a preset time period according to the wind speed time history; Determine the average wind speed distribution matrix based on the average wind speed and the nodes , where the number of nodes in the horizontal direction is k and the number of nodes in the vertical direction is h, is the average wind speed of the node at the hth row and kth column in the i direction within the preset time, where i=1, 2, and 3 represent the longitudinal, horizontal, and vertical directions respectively.

7. A wind condition determination system suitable for a wind turbine generator system, characterized in that: include: A setting unit, used for setting a plurality of nodes on the plane where the wind rotor of the wind turbine generator set is located; a first determining unit, configured to determine an actual wind speed time series of each of the nodes, wherein the actual wind speed time series is a correspondence between actual wind speed and time; A second determining unit is configured to determine an average wind speed distribution of each of the nodes based on the actual wind speed time history of each of the nodes, wherein the average wind speed distribution is an average value of the actual wind speed within a preset time; a third determining unit, configured to determine a simulated wind speed time history of each of the nodes according to the actual wind speed time history, wherein the simulated wind speed time history represents a plurality of wind speeds of each of the nodes; a fourth determining unit, configured to determine a dynamic response of the wind turbine generator set according to the actual wind speed time history, the average wind speed distribution, or the simulated wind speed time history, the dynamic response including power generation of the wind turbine generator set, component loads, and a distance between blades and a tower of the wind turbine generator set; Also includes: The first conversion unit is used to convert each wind speed of the actual wind speed time series from the time domain to the frequency domain to obtain a wind spectrum characteristic distribution matrix ,in is the wind power spectrum density value in direction i at the node in row h and column k within the preset time; a sixth determining unit, configured to determine an autospectral value of a wind speed spectral density matrix according to the wind spectrum characteristic distribution matrix, and determine a cross-spectral value of the wind speed spectral density matrix according to a spatial coherence function; The third determining unit is specifically configured to convert the wind speed spectral density matrix from the frequency domain to the time domain to obtain the simulated wind speed time history; Decomposition unit, used to decompose the wind speed spectrum density matrix According to the first relation to decompose; The third determining unit is specifically configured to convert the wind speed spectrum density matrix from the frequency domain to the time domain according to the second relationship to obtain the simulated wind speed time history. ; Among them, t is time, i=1, 2, 3 represent the longitudinal, horizontal and vertical directions respectively, y j and z j represent the horizontal and vertical coordinates of the j-th point, respectively. represents the frequency increment, w up is the cutoff circular frequency, N is an integer greater than 1; is a random phase angle in the interval Evenly distributed within; for The complex angle, Take the value of the jth row and mth column in the matrix Hi(w), , the number of nodes in the vertical direction is h.

8. A wind condition determination device suitable for a wind turbine generator system, characterized in that: include: memory for storing computer programs; A processor is configured to implement the steps of the method for determining wind conditions applicable to a wind turbine generator set as claimed in any one of claims 1 to 6 when executing the computer program.

Images