Wake Parameter Determination Method, Device, Equipment and Computer Readable Storage Medium

Abstract

Embodiments of the present application provide a wake parameter determination method, apparatus, device, and computer-readable storage medium. The wake parameter determination method includes: obtaining wind farm data, operating data of a wind turbine generator set, and an actuator disk model of the impeller of the wind turbine generator set, where the operating data of the wind turbine generator set includes a first reference wind speed and a second reference wind speed; determining free stream data based on the wind farm data, operating data, and the actuator disk model, and determining the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model; and determining the wake parameters of the wind turbine generator set based on the actuator disk model according to the wind acceleration factor and the relative wind direction deflection angle. The wake parameter determination method for the wind turbine generator set provided by the present application can ensure better calculation accuracy, reduce the calculation amount, the result is closer to the actual operating result, and can better evaluate the influence of the wake existing in a large-scale wind farm on the power generation loss.

Description

Technical Field

[0001] The present application relates to the technical field of wind power generation. Specifically, the present application relates to a method, device, equipment and computer-readable storage medium for determining wake parameters. Background Art

[0002] A wake refers to the turbulent vortex flow behind a moving object or downstream of an object, also known as a wake. For a large-scale wind farm, the downstream wind turbines may be in the wake of the upstream wind turbines, which will greatly affect the power generation performance and safety of the downstream wind turbines. Therefore, accurately evaluating the wake impact of a wind farm can, on the one hand, accurately calculate the grid-connected power of wind turbine generators to ensure power generation benefits, and on the other hand, avoid the threat of inefficient assets and the safety of the turbines caused by too low downstream wind speed due to the wake.

[0003] Currently, the actuator disk model is considered an effective wake impact evaluation model. It replaces the wind wheel with an actuator disk and adds a volume momentum source to the NS equation to simulate the effect of the blades virtually, thus avoiding the simulation of actual blades and being a method for evaluating the wake impact of wind turbines. However, when the existing actuator disk model calculates the wake of a wind turbine, it requires inlet conditions of full wind speed and full wind direction for calculation, with a long calculation time and low efficiency, which cannot meet the needs of daily wind resource assessment services. Summary of the Invention

[0004] In view of the shortcomings of the existing methods, the present application provides a method, device, equipment and computer-readable storage medium for determining wake parameters to solve the technical problem of inaccurate evaluation of wake losses in a wind farm existing in the prior art.

[0005] In a first aspect, an embodiment of the present application provides a method for determining wake parameters of a wind turbine generator set, including:

[0006] Obtain wind farm data, operation data of the wind turbine generator set, and an actuator disk model of the impeller of the wind turbine generator set. The operation data of the wind turbine generator set includes a first reference wind speed and a second reference wind speed;

[0007] Determine free stream data based on the wind farm data, operation data, and the actuator disk model, and determine the wind acceleration factor and relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model;

[0008] Determine the wake parameters of the wind turbine generator set based on the actuator disk model according to the wind acceleration factor and the relative wind direction deflection angle.

[0009] In some implementations of the first aspect, the free stream data includes free stream wind speed data; determining the wind acceleration factor and relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model includes:

[0010] The actuator disk simulation data is determined based on the wind farm data, operation data, and actuator disk model;

[0011] According to a plurality of preset wind directions, the computational domain corresponding to the impeller in the actuator disk model is divided into a plurality of sectors;

[0012] According to the actuator disk simulation data, free stream wind speed data, and sectors, determine the wind acceleration factor of each sector and determine the relative wind direction deflection angle of each sector.

[0013] Combined with the first aspect and the above implementation, in some implementations of the first aspect, determining the wind acceleration factor of each sector includes:

[0014] Determine the first calibration wind speed and the second calibration wind speed according to the first reference wind speed, the second reference wind speed, and the free stream wind speed data, and divide the preset wind speed interval into three wind speed sub-intervals;

[0015] According to the wind speed sub-intervals and the free stream wind speed data of each sector, determine the actuator disk wind acceleration factor of each sector.

[0016] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the wind speed sub-intervals include a first wind speed sub-interval, a second wind speed sub-interval, and a third wind speed sub-interval. The first wind speed sub-interval is the wind speed interval less than the first calibration wind speed, the second wind speed sub-interval is the wind speed interval greater than or equal to the first calibration wind speed and less than or equal to the second calibration wind speed, and the third wind speed sub-interval is the wind speed interval greater than the third calibration wind speed;

[0017] And, according to the wind speed sub-intervals and the free stream wind speed data of each sector, determining the actuator disk wind acceleration factor of each sector includes:

[0018] If the free stream wind speed data of a sector is within the first wind speed sub-interval, determine the wind acceleration factor of the sector as the first actuator disk wind acceleration factor;

[0019] If the free stream wind speed data of a sector is within the third wind speed sub-interval, determine the wind acceleration factor of the sector as the third actuator disk wind acceleration factor;

[0020] If the free stream wind speed data of a sector is within the second wind speed sub-interval, determine the wind acceleration factor of the sector as the second actuator disk wind acceleration factor according to the first calibration wind speed, the second calibration wind speed, the first actuator disk wind acceleration factor, and the third actuator disk wind acceleration factor.

[0021] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the operation data of the wind turbine includes the measured data of the anemometer tower; determining the relative wind direction deflection angle of each sector includes:

[0022] Based on the simulation data of the actuator disk and the measured data of the wind measurement tower, determine the wind direction data of the wind measurement tower, the wind direction data of the machine position point, and the inlet wind direction data of each sector;

[0023] Based on the wind direction data of the wind measurement tower and the inlet wind direction data, determine the wind direction deflection angle of the wind measurement tower;

[0024] Based on the first linear interpolation function between the wind direction deflection angle and the inlet wind direction data of the wind measurement tower in each sector, determine the interpolated wind direction data of the wind measurement tower, and obtain the first functional relationship between the wind direction deflection angle of the wind measurement tower and the interpolated wind direction data of the wind measurement tower;

[0025] Based on the wind direction data of the machine position point and the inlet wind direction data, determine the wind direction deflection angle of the machine position point;

[0026] Based on the wind direction deflection angle of the machine position point and the inlet wind direction data, determine the interpolated wind direction data of the machine position point, and obtain the second functional relationship between the inlet wind direction data and the interpolated wind direction data of the machine position point;

[0027] Based on the first functional relationship and the second functional relationship, determine the relative wind direction deflection angle of each sector.

[0028] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the measured data of the wind measurement tower includes wind speed and wind direction; after determining the wind acceleration factor of each sector and the relative wind direction deflection angle of each sector, it further includes:

[0029] According to the preset wind speed interval and the preset wind direction interval, determine a number of wind speed zones and a number of wind direction zones;

[0030] Generate a wind frequency matrix according to the partition frequencies corresponding to the wind speed zones and the wind direction zones;

[0031] According to the wind acceleration factor and the relative wind direction deflection angle corresponding to each element in the wind frequency matrix, determine the wind speed value and the wind direction value at the machine position corresponding to each element in the wind frequency matrix;

[0032] Round the wind speed value and the wind direction value at the machine position to determine the actuator disk wind frequency matrix of each machine position point in each sector;

[0033] According to the actuator disk wind frequency matrix of each machine position point, determine the actuator disk power generation data.

[0034] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the measured data of the wind measurement tower is obtained from a number of wind measurement towers; determining the actuator disk wind frequency matrix of each machine position point in each sector includes:

[0035] Based on the actuator disk wind frequency matrices of each machine position point corresponding to a number of wind measurement towers, determine the wind frequency matrix of any one of the machine position points by weighting.

[0036] Combining the first aspect and the above implementation manners, in some implementation manners of the first aspect, the operation data of the wind turbine generator set further includes the unit power data of the wind turbine generator set; according to the actuation disk wind frequency matrix of each machine position point, the actuation disk power generation data is determined, including:

[0037] According to the unit power data, determine the power value at the first preset wind speed and the power value at the second preset wind speed in adjacent wind speed segments of any one machine position point;

[0038] According to the power value at the first preset wind speed and the power value at the second preset wind speed, obtain the average power of any one machine position;

[0039] According to the actuation disk wind frequency matrix and the average power of any one machine position point, determine the actuation disk power generation data.

[0040] Combining the first aspect and the above implementation manners, in some implementation manners of the first aspect, the actuation disk simulation data determined according to the wind farm data, the operation data and the actuation disk model includes:

[0041] Construct an actuation disk three-dimensional space grid according to the operation data of the wind turbine generator set;

[0042] Determine the boundary conditions according to the wind farm data, the first reference wind speed and the second reference wind speed, and perform initialization processing on the actuation disk three-dimensional space grid to obtain an initialized actuation disk three-dimensional space grid;

[0043] According to the computational fluid dynamics simulation determination method, perform iterative determination processing on the initialized actuation disk three-dimensional space grid to obtain the converged actuation disk simulation data.

[0044] Combining the first aspect and the above implementation manners, in some implementation manners of the first aspect, the wind farm data includes the terrain data of the wind turbine generator set, the machine position point coordinate data of the wind turbine generator set and the anemometer tower coordinate data;

[0045] Constructing an actuation disk three-dimensional space grid according to the operation data of the wind turbine generator set includes:

[0046] Generate a surface grid according to the terrain data of the wind turbine generator set;

[0047] According to the machine position point coordinate data and the anemometer tower coordinate data of the wind turbine generator set, increase the grid density of the surface grid to determine the encrypted surface grid;

[0048] Perform surface normal stretching on the encrypted surface grid to determine the actuation disk three-dimensional space grid.

[0049] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the operating data of the wind turbine includes the roughness map data corresponding to the wind turbine; obtaining the initialized actuator disk three-dimensional space grid includes:

[0050] According to the roughness map data, add the surface roughness value to the grid center point of the actuator disk three-dimensional space grid, and set the inlet boundary conditions to obtain the initialized actuator disk three-dimensional space grid respectively.

[0051] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the operating data of the wind turbine further includes reference height data and a plurality of inlet reference wind speed data; setting the inlet boundary conditions includes:

[0052] According to the wind profile equation, the turbulent kinetic energy equation, the turbulent dissipation equation, the surface roughness value, the reference height data, the first reference wind speed and the second reference wind speed, determine the inlet boundary conditions.

[0053] Combined with the first aspect and the above implementation, in some implementations of the first aspect, the software toolkit used in the computational fluid dynamics simulation determination method includes an open-source field operation and manipulation software toolkit;

[0054] According to the computational fluid dynamics simulation determination method, the steps of performing iterative determination processing on the initialized actuator disk three-dimensional space grid to obtain the converged actuator disk simulation data include: respectively performing iterative determination on the initialized actuator disk three-dimensional space grid to convergence according to the open-source field operation and manipulation software toolkit and the turbulence model to obtain the converged actuator disk simulation data;

[0055] The converged actuator disk simulation data includes wind speed vector data and turbulent kinetic energy data.

[0056] Combined with the first aspect and the above implementation, in some implementations of the first aspect, determining the free stream data includes:

[0057] According to the free stream model and the operating data of the wind turbine, determine the free stream simulation data;

[0058] According to the free stream simulation data, determine the free stream data including the free stream wind speed data.

[0059] Combined with the first aspect and the above implementation, in some implementations of the first aspect, determining the free stream simulation data includes:

[0060] Construct a three-dimensional space grid according to the operating data of the wind turbine;

[0061] Determine boundary conditions based on the operating data of the wind turbine generator set, and perform initialization processing on the three-dimensional space grid to obtain an initialized three-dimensional space grid;

[0062] According to the computational fluid dynamics simulation determination method, perform iterative determination processing on the initialized three-dimensional space grid to obtain the converged free flow simulation data.

[0063] Combined with the first aspect and the above implementation manners, in some implementation manners of the first aspect, the operating data of the wind turbine generator set includes the unit power data of the wind turbine generator set and the measured data of the anemometer tower;

[0064] According to the free flow simulation data, determine the free flow data including the free flow wind speed data, including:

[0065] According to the free flow simulation data and the measured data of the anemometer tower, determine the free flow wind acceleration factor and the free flow relative wind direction deflection angle;

[0066] According to the measured data of the anemometer tower, the free flow wind acceleration factor and the free flow relative wind direction deflection angle, determine the free flow wind frequency matrix of each machine position point;

[0067] According to the free flow wind frequency matrix of each machine position point, determine the free flow wind speed data.

[0068] In a second aspect, an embodiment of the present application provides a wake parameter determination device for a wind turbine generator set, including:

[0069] An acquisition module, configured to acquire wind farm data, the operating data of the wind turbine generator set, and the actuator disk model of the impeller of the wind turbine generator set, where the operating data of the wind turbine generator set includes a first reference wind speed and a second reference wind speed;

[0070] A simulation module, configured to determine free flow data according to the wind farm, the operating data, and the actuator disk model, and determine the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model;

[0071] A determination module, configured to determine the wake parameters of the wind turbine generator set based on the actuator disk model according to the wind acceleration factor and the relative wind direction deflection angle.

[0072] In a third aspect, an embodiment of the present application provides a wake parameter determination device for a wind turbine generator set, including:

[0073] A processor;

[0074] A memory, electrically connected to the processor;

[0075] At least one program, stored in a memory and configured to be executed by a processor, the at least one program being configured to: implement the method for determining wake parameters of a wind turbine as described in the embodiment of the first aspect of the present application.

[0076] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, which stores at least one segment of program or code set, and the at least one segment of program or code set is loaded and executed by a processor to implement the method for determining wake parameters of a wind turbine as described in the embodiment of the first aspect of the present application.

[0077] The beneficial technical effects brought by the technical solution provided by the embodiment of the present application are:

[0078] The method for determining wake parameters of a wind turbine provided by the present application regards the impeller movement trajectory of the wind turbine as an actuator disk model. In the actuator disk model, only two reference wind speeds need to be selected, and based on the wind farm data, operation data, and data obtained from the free flow model, the wind acceleration factor and relative wind direction deflection angle of each sector of the impeller wake in the actuator disk model are calculated to determine the wake parameters of the wind turbine, which can ensure better calculation accuracy, while reducing the calculation amount. The calculation result is closer to the actual operation result, and can better evaluate the influence of the wake existing in a large-scale wind farm on the power generation loss.

[0079] The additional aspects and advantages of the present application will be partially given in the following description, which will become obvious from the following description, or can be understood through the practice of the present application. Description of the Drawings

[0080] The above-mentioned and / or additional aspects and advantages of the present application will become obvious and easy to understand from the following description of the embodiments in conjunction with the drawings, where:

[0081] Figure 1 is a schematic flowchart of a method for determining wake parameters of a wind turbine provided by an embodiment of the present application;

[0082] Figure 2 is a schematic flowchart of a method for determining the wind acceleration factor and relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model provided by an embodiment of the present application;

[0083] Figure 3 is a schematic flowchart of a method for determining the actuator disk simulation data according to the wind farm data, operation data, and actuator disk model provided by an embodiment of the present application;

[0084] Figure 4 is a schematic flowchart of another method for determining wake parameters of a wind turbine provided by an embodiment of the present application;

[0085] Figure 5 Schematic structural framework diagram of a wake parameter determination device for a wind turbine provided by an embodiment of the present application;

[0086] Figure 6 Schematic structural framework diagram of a wake parameter determination device for a wind turbine provided by an embodiment of the present application. Specific implementation manners

[0087] The present application will be described in detail below. Examples of embodiments of the present application are shown in the accompanying drawings, where the same or similar reference numerals represent the same or similar components or components with the same or similar functions throughout. In addition, if the detailed description of the known technology is unnecessary for showing the features of the present application, it will be omitted. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present application and cannot be construed as a limitation of the present application.

[0088] Those skilled in the art of the present technology can understand that, unless otherwise defined, all terms (including technical terms and scientific terms) used herein have the same meaning as the general understanding of those of ordinary skill in the art to which the present application belongs. It should also be understood that terms such as those defined in a general dictionary should be understood to have a meaning consistent with the meaning in the context of the prior art, and will not be interpreted with an idealized or overly formal meaning unless specifically defined as here.

[0089] Those skilled in the art of the present technology can understand that, unless specifically stated, the singular forms "a", "an", "the" and "said" used herein may also include the plural forms. It should be further understood that the term "comprising" used in the specification of the present application means the presence of the described features, integers, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or their groups. It should be understood that the phrase "and / or" used herein includes all or any unit and all combinations of one or more of the associated listed items.

[0090] The technical solution of the present application and how the technical solution of the present application solves the above technical problems will be described in detail below with specific embodiments.

[0091] In the first aspect, an embodiment of the present application provides a method for determining wake parameters of a wind turbine, as Figure 1 shown, including the following steps:

[0092] S100: Obtain wind farm data, operating data of the wind turbine, and an actuator disk model of the impeller of the wind turbine. The operating data of the wind turbine includes a first reference wind speed and a second reference wind speed.

[0093] S200: Determine the free - stream data based on the wind farm data, operation data, and actuator disk model, and determine the wind acceleration factor and relative wind direction deflection angle for each sector of the impeller wake based on the actuator disk model.

[0094] S300: Determine the wake parameters of the wind turbine based on the actuator disk model according to the wind acceleration factor and relative wind direction deflection angle.

[0095] The actuator disk model mentioned in the above - mentioned method for determining the wake parameters of the wind turbine means that the trajectory shape swept by the impeller of the wind turbine is circular, and the impeller part is simplified into a cylindrical thin disk with a certain thickness, which is the actuator disk. During the rotation of the impeller, by absorbing the kinetic energy of the incoming wind and converting it into its own rotational mechanical energy, it will cause sudden changes in the wind speed and pressure of the incoming wind. A volume momentum source is added to the NS equations (Navier - Stokes equations) to simulate the effect of the blades virtually, avoiding the complex simulation of the actual blades.

[0096] The method for determining the wake parameters of the wind turbine provided in this application regards the impeller motion trajectory of the wind turbine as an actuator disk model, selects two reference wind speeds in the actuator disk model, and calculates the wind acceleration factor and relative wind direction deflection angle for each sector of the impeller wake according to the wind farm data, operation data, and data obtained from the free - stream model, thereby determining the wake parameters of the wind turbine. This can ensure better calculation accuracy, while reducing the calculation amount. The calculation result is closer to the actual operation result, and can better evaluate the influence of the wake existing in a large - scale wind farm on the power generation loss.

[0097] For S100, through detection devices or sensors in a wind farm including multiple wind turbines, such as a wind measurement tower, fan machinery, and electrical energy detection equipment, etc., obtain the wind farm data and the operation data of the wind turbine. Abstract the impeller of the wind turbine as an actuator disk model and perform subsequent calculations. To ensure calculation accuracy and at the same time not make the calculation time too long, two reference wind speeds in the operation data of the wind turbine are selected. Specifically, according to the power curve of the wind turbine, the wind speeds before and after the rated wind speed are selected, that is, a certain wind speed before the wind turbine reaches the rated power and a certain wind speed after reaching the rated power. The specific wind speed selection can be determined according to long - term calculation and evaluation time.

[0098] For S200, based on the data obtained from S100, through free - flow modeling and simulation comprehensive calculation, free - flow data is prepared. Then, through actuator - disk modeling and simulation comprehensive calculation, various parameters for calculating the wake parameters of the wind turbine are determined, such as the wind acceleration factor and relative wind - direction deflection angle in each sector of the impeller wake. For each sector of the impeller wake involved in S200, the determined computational domain can be evenly divided into 16 equal parts, and each equal part is a sector of the impeller wake.

[0099] For S300, the final stage of actuator - disk modeling and simulation comprehensive calculation relies on the data calculated previously or obtained directly, such as the wind acceleration factor and relative wind - direction deflection angle, etc., to determine the wake parameters of the wind turbine. The wake parameters of the wind turbine include the actuator - disk wake wind speed and actuator - disk wake power generation, etc.

[0100] Optionally, in an implementation manner of the first - aspect embodiment of the present application, the free - flow data includes free - flow wind - speed data. As Figure 2 shown, determining the wind acceleration factor and relative wind - direction deflection angle in each sector of the impeller wake based on the actuator - disk model in S200 specifically includes:

[0101] S210: Determine the actuator - disk simulation data according to the wind - farm data, operation data, and actuator - disk model.

[0102] S220: Divide the actuator - disk corresponding to the impeller in the actuator - disk model into several sectors according to a number of preset wind directions.

[0103] S230: Determine the wind acceleration factor in each sector and the relative wind - direction deflection angle in each sector according to the actuator - disk simulation data, free - flow wind - speed data, and sectors.

[0104] According to the content described above, to determine the wake parameters of the wind turbine, it is necessary to first determine the wind acceleration factor and relative wind - direction deflection angle in each sector of the actuator - disk model. The determination of these two parameters requires obtaining the actuator - disk simulation data, free - flow wind - speed data, and the divided sectors. The division of the sectors is as described above, and the computational domain in the actuator - disk model corresponding to the impeller trajectory is evenly divided into multiple sub - regions, such as 16 sub - regions, that is, 16 sectors. The actuator - disk simulation data is calculated according to the wind - farm data, the operation data of the wind farm, and the selected actuator - disk model, and the specific calculation content will be described in detail later. The free - flow wind - speed data is a type of free - flow data, which is determined according to the free - flow model, the wind - farm data, the operation data of the wind farm, etc. The free - flow model is an ideal model that does not consider the influence of the wake in the wind turbine.

[0105] Optionally, combined with the above - mentioned implementation manner of the present application, in some implementation manners of the first - aspect embodiment of the present application, asFigure 3 As shown, S210 determines the actuator disk simulation data based on the wind farm data, operation data, and actuator disk model, including:

[0106] S211: Construct a three-dimensional space grid of the actuator disk according to the operation data of the wind turbine generator set.

[0107] S212: Determine the boundary conditions according to the wind farm data, the first reference wind speed, and the second reference wind speed, and perform initialization processing on the three-dimensional space grid of the actuator disk to obtain the initialized three-dimensional space grid of the actuator disk.

[0108] S213: Perform iterative determination processing on the initialized three-dimensional space grid of the actuator disk according to the computational fluid dynamics simulation determination method to obtain the converged actuator disk simulation data.

[0109] Optionally, in some implementation manners of the first aspect embodiment of the present application, the wind farm data includes the terrain data of the wind turbine generator set, the machine site coordinate data of the wind turbine generator set, and the anemometer tower coordinate data. In the foregoing steps, constructing a three-dimensional space grid of the actuator disk according to the operation data of the wind turbine generator set specifically includes:

[0110] Generate a surface grid according to the terrain data of the wind turbine generator set. Then, according to the machine site coordinate data of the wind turbine generator set and the anemometer tower coordinate data, increase the grid density of the surface grid to determine the encrypted surface grid. Then, perform surface normal stretching on the encrypted surface grid to determine the three-dimensional space grid of the actuator disk.

[0111] For example, process the terrain data in the original map format to determine the center point and calculation range of the computational domain, and then generate three-dimensional grids according to the number of sectors required for calculation. The steps of generating the three-dimensional grids are as follows:

[0112] Step 1: Use SnappyHexMesh in OpenFOAM to generate a surface grid of the map surface.

[0113] Step 2: Perform additional grid encryption at the anemometer tower and the machine site of the wind turbine generator set to make the parameters of the surface grid at these positions more comprehensive and simulate the actual situation at the anemometer tower and the machine site more finely, and better fit the real terrain data. Among them, an actuator disk grid is generated at the impeller of the machine site using the actuator disk model.

[0114] Step 3: Vertically stretch the surface grid to generate a volume grid and form the entire computational domain.

[0115] Optionally, in combination with the above implementation, in some implementations of the first aspect embodiment of the present application, the operation data of the wind turbine includes the roughness map data corresponding to the wind turbine; obtaining the initialized actuator disk three-dimensional space grid includes: according to the roughness map data, adding the surface roughness value to the grid center point of the actuator disk three-dimensional space grid, and setting the inlet boundary conditions to obtain the initialized actuator disk three-dimensional space grid respectively. For example, adding the surface roughness value z0 to the grid center point of the actuator disk three-dimensional space grid and setting the inlet boundary conditions, that is, setting the atmospheric boundary layer conditions. Optionally, the operation data of the wind turbine further includes the reference height data and several inlet reference wind speed data. Among them, setting the inlet boundary conditions specifically includes:

[0116] Determine the inlet boundary conditions according to the wind profile equation, turbulent kinetic energy equation, turbulent dissipation equation, surface roughness value, reference height data, first reference wind speed and second reference wind speed.

[0117] Optionally, the wind profile equation is Formula 1, the turbulent kinetic energy equation is Formula 2, and the turbulent dissipation equation is Formula 3, which are specifically as follows:

[0118]

[0119]

[0120]

[0121] In the above formulas, u * is the surface friction velocity of the atmospheric boundary layer. Under neutral atmospheric conditions, the inlet reference wind speed is v = 10 m / s, and the reference height is the height from the ground z = 500 m. Substituting into Formula 1, u * can be obtained. Among them, k0 is the turbulent kinetic energy of the inlet initial condition, κ is the von Kármán constant, which can be taken as 0.41, C μ is the turbulent model constant, usually can be taken as 0.09, and ε is the turbulent dissipation rate.

[0122] Optionally, the software tool package adopted by the computational fluid dynamics simulation determination method includes an open-source field operation and operation software package. According to the computational fluid dynamics simulation determination method, the steps of performing iterative determination processing on the initialized actuator disk three-dimensional space grid to obtain the converged actuator disk simulation data include: respectively iteratively determining the initialized actuator disk three-dimensional space grid to convergence according to the open-source field operation and operation software package and the turbulence model to obtain the converged actuator disk simulation data. The converged actuator disk simulation data includes wind speed vector data and turbulent kinetic energy data.

[0123] The method for determining computational fluid dynamics simulation is CFD (Computational Fluid Dynamics). The software toolkits that can be used in this method include the open-source field operation and manipulation software package, namely OpenFOAM (Field Operation and Manipulation). For each set of initialized three-dimensional grids of the actuator disk, CFD simulation calculations are performed using OpenFOAM. Among them, the standard k-ε two-equation turbulence model is selected for the turbulence model, that is, the aforementioned formulas 2 and 3 are used. The above toolkits and CFD simulation calculations are continuously used, and specific iterative calculations are performed until convergence to obtain the converged simulation data of the actuator disk.

[0124] Optionally, according to the aforementioned iterative calculation results, the simulation data of the actuator disk at the anemometer towers in each sector and each machine position of the wind turbine are extracted. These simulation data of the actuator disk can satisfy the wind speed vectors and turbulent kinetic energies in all hub height ranges of existing products. The wind direction and the inflow angle can be calculated based on the three velocity components in space. Among them, the horizontal wind speed can be calculated by formula 4, and the formula for calculating turbulence is formula 5. In addition, the calculation formula for turbulent kinetic energy is formula 6. The data used in formula 6 is obtained after iterative calculation and convergence according to the aforementioned inlet initial conditions such as k0 and other formulas. Specifically, formulas 4 to 6 are as follows:

[0125]

[0126]

[0127] k = (σ x 2 + σ y 2 + σ z 2 ) / 2 …………………… formula (6);

[0128] In the above formulas 4 to 6, the turbulence in the horizontal direction is calculated. Therefore When the flow of the wind current is isotropic, then Then the turbulence v x refers to the wind speed in the X-axis direction in the space coordinate system, v y refers to the wind speed in the Y-axis direction, v z refers to the wind speed in the Z-axis direction. The safety of the blades in the wind turbine can be characterized by wake parameters such as turbulent kinetic energy, turbulent dissipation rate, and the magnitude of turbulence.

[0129] Optionally, in combination with the above implementation manner, in some implementation manners of the first aspect embodiment of the present application, determining the wind acceleration factor of each sector includes:

[0130] Determine the first calibrated wind speed and the second calibrated wind speed based on the first reference wind speed, the second reference wind speed, and the free-stream wind speed data, and divide the preset wind speed interval into three wind speed sub-intervals. Determine the actuator disk wind acceleration factor for each sector according to the wind speed sub-intervals and the free-stream wind speed data of each sector.

[0131] Optionally, in combination with the above implementation, in one implementation of the above embodiment, the wind speed sub-intervals include a first wind speed sub-interval, a second wind speed sub-interval, and a third wind speed sub-interval. The first wind speed sub-interval is the wind speed interval less than the first calibrated wind speed. The second wind speed sub-interval is the wind speed interval greater than or equal to the first calibrated wind speed and less than or equal to the second calibrated wind speed. The third wind speed sub-interval is the wind speed interval greater than the third calibrated wind speed. And, determining the actuator disk wind acceleration factor for each sector according to the wind speed sub-intervals and the free-stream wind speed data of each sector specifically includes:

[0132] If the free-stream wind speed data of a sector is within the first wind speed sub-interval, determine the wind acceleration factor of the sector as the first actuator disk wind acceleration factor.

[0133] If the free-stream wind speed data of a sector is within the third wind speed sub-interval, determine the wind acceleration factor of the sector as the third actuator disk wind acceleration factor.

[0134] If the free-stream wind speed data of a sector is within the second wind speed sub-interval, determine the wind acceleration factor of the sector as the second actuator disk wind acceleration factor according to the first calibrated wind speed, the second calibrated wind speed, the first actuator disk wind acceleration factor, and the third actuator disk wind acceleration factor.

[0135] Specifically, the above implementation and / or embodiment are as follows: When the anemometer tower is in the free-stream condition (i.e., not affected by the wake), the free-stream wind speed of a certain sector is determined to be V mast (For example, the wind speed data calculated at a reference wind speed of 10 m / s without an actuator disk model, the same as the anemometer tower wind speed described below). The preset wind speed interval is determined according to the possible minimum wind speed and maximum wind speed in the wind farm. For example, the preset wind speed interval of 0 - 50 m / s is divided into 0 - V mast_v1 、V mast_v1 ~V mast_v2 and V mast_v2 ~50 three intervals. This is to convert the anemometer tower calculated wind speed with a reference wind speed of 10 m / s to the calculated wind speed corresponding to the first reference wind speed v1 and the second reference wind speed v2, that is, the first calibrated wind speed V mast_v1 and the second calibrated wind speed V mast_v2 .

[0136] When the anemometer tower wind speed V mast is within the first wind speed sub-interval, that is, 0 - V mast_v1When it is within the range, the wind acceleration factor at the upwind point i of the wind turbine generator is When the wind speed V of the anemometer tower mast Within the third wind speed sub-interval, that is, V mastv1 ~50 within the range, the wind acceleration factor at the point i of the machine is In the calculation formulas of the above two wind acceleration factors, Vi_v1 and Vi_v2 are the calculated wind speeds at the point i of the machine at the first reference wind speed and the calculated wind speeds at the point i of the machine at the second reference wind speed, respectively.

[0137] And when the wind speed V of the anemometer tower mast Within the second sub-interval, that is, V mast_v1 ~V mast_v2 Within the range, the wind acceleration factor at the point i of the machine is a_idir_vj = f(Vmast), where f(V mast ) represents constructing a linear interpolation function according to Vmast_v1 and V mast_v2 And a_idir_v1 and a_idir_v2, and the relationship between the wind acceleration factor of the actuator disk model machine point and the inlet wind direction is shown in Table 1.

[0138] Table 1 Relationship between inlet wind direction and wind acceleration factor at point i of the machine

[0139]

[0140] The specific determination methods of the above-mentioned first calibration wind speed and second calibration wind speed are specifically Formulas 7 and 8:

[0141] V mast_v1 =(V mast / 10)*v1…………………………Formula (7),

[0142] V mast_v2 =(V mast / 10)*v2…………………………Formula (8).

[0143] Optionally, in combination with the above implementation method, in some implementation methods of the first aspect embodiment of the present application, the operation data of the wind turbine generator includes the measured data of the anemometer tower. Determining the relative wind direction deflection angle of each sector specifically includes:

[0144] First, according to the actuator disk simulation data and the measured data of the anemometer tower, determine the anemometer tower wind direction data, the machine point wind direction data, and the inlet wind direction data of each sector.

[0145] Secondly, according to the anemometer tower wind direction data and the inlet wind direction data, determine the wind direction deflection angle of the anemometer tower.

[0146] Again, according to the first linear interpolation function between the wind direction deflection angles of the anemometer towers in each sector and the inlet wind direction data, determine the anemometer tower wind direction interpolation data, and obtain the first functional relationship between the wind direction deflection angle of the anemometer tower and the anemometer tower wind direction interpolation data.

[0147] Next, according to the wind direction data at the machine point and the inlet wind direction data, determine the wind direction deflection angle at the machine point.

[0148] Then, according to the wind direction deflection angle at the machine point and the inlet wind direction data, determine the wind direction interpolation data at the machine point, and obtain the second functional relationship between the inlet wind direction data and the wind direction interpolation data at the machine point.

[0149] Finally, according to the first functional relationship and the second functional relationship, determine the relative wind direction deflection angles of each sector.

[0150] The above process can be understood with reference to the following case:

[0151] After obtaining the actuator disk simulation data and the anemometer tower measured data, the actuator disk model can be divided into 16 sectors, thereby defining the inlet wind direction data of each sector, denoted by the symbol dir_n, and correspondingly obtaining the anemometer tower wind direction data and the wind direction data at the machine point. Define the wind direction deflection angles of the anemometer towers in 16 sectors as Delta_mast_dir_n. The wind direction deflection angle of the anemometer tower is the difference between the anemometer tower wind direction and the corresponding anemometer tower inlet wind direction. A linear interpolation function, that is, the first linear interpolation function, can be constructed between the wind direction deflection angle of the anemometer tower and the anemometer tower inlet wind direction. Interpolate to obtain the anemometer tower inlet wind direction Delta_mast_idir for each inlet wind direction idir (an integer value in the range of 0 to 359, that is, the inlet wind direction for each degree), and thus obtain the relationship between the inlet wind direction idir and the anemometer tower wind direction D_mast_idir, that is, D_mast_idir = Delta_mast_idir + idir, as shown in Table 2.

[0152] Table 2 Relationship between inlet wind direction and calculated wind direction of anemometer tower

[0153]

[0154] Similarly, in sector dir_n, the value of the wind direction deflection angle at the i-th turbine location of the anemometer tower and the wind turbine, that is, the relative wind direction deflection angle of the turbine location relative to the anemometer tower, is specifically Δi_dir_n = Δi_dir_n - Δmast_dir_n, which is the difference between the wind direction at the turbine location Δmast_dir_n and the corresponding inlet wind direction at the turbine location Δi_dir_n, and can also be written as the difference between the wind direction at the turbine location and the wind direction at the anemometer tower: Δi_mast_dir_n = Di_dir_n - D_mast_dir_n. A linear interpolation function, that is, the second functional relationship, can also be constructed between Δi_mast_dir_n and the inlet wind direction dir_n of each sector to interpolate and obtain the relative deflection angle Δi_mast_idir between the i-th turbine location and the anemometer tower for each idir (an integer value in the range of 0 to 359). Thus, the relationship between idir and Δi_mast_idir is obtained, as shown in Table 3.

[0155] Table 3 Relationship between inlet wind direction and relative deflection angle of the i-th turbine location - anemometer tower

[0156] idir 0 1 2 … n … 359 Delta_i_mast_idir d0 d1 d2 … dn … d359

[0157] Optionally, in combination with the above implementation, in an implementation manner of the above embodiment, the measured data of the anemometer tower includes wind speed and wind direction. After the steps of determining the wind acceleration factor of each sector and determining the relative wind direction deflection angle of each sector in the above steps, the following steps are further included:

[0158] According to the preset wind speed interval and preset wind direction interval, a number of wind speed zones and a number of wind direction zones are determined. According to the partition frequencies corresponding to the wind speed zones and wind direction zones, a wind frequency matrix is generated. According to the wind acceleration factor and relative wind direction deflection angle corresponding to each element in the wind frequency matrix, the wind speed value and wind direction value at the turbine location corresponding to each element in the wind frequency matrix are determined. The wind speed value and wind direction value at the turbine location are rounded to determine the actuation disk wind frequency matrix of each turbine location in each sector. According to the actuation disk wind frequency matrix of each turbine location, the actuation disk power generation data is determined.

[0159] The measured data of the anemometer tower is usually a time series (the measured data includes wind speed and wind direction) data. For example, the possible wind speed range of 0 to 50 m / s (meters per second) and the possible direction range of 0 to 359° (degrees) are partitioned according to a preset wind speed interval of 0.1 m / s and a preset wind direction interval of 1°, respectively. The partition frequencies of each partition for wind speed and wind direction are statistically calculated to generate a wind frequency matrix fvvdd_mast (the data volume is 500 * 360). Arbitrarily select a partition from fvvdd_mast. For example, if the wind speed is iv and the wind direction is id, look up the corresponding wind acceleration factor a_id_iv and the relative wind direction deflection angle Delta_i_mast_idir in the aforementioned Table 2 and Table 4 respectively. Then the wind speed value at a certain machine position is iv * a_id_iv, and the wind direction value is Delta_i_mast_idir + id. Round the wind speed value iv * a_id_iv and the wind direction value Delta_i_mast_idir + id respectively. After the operation, we get:

[0160] fvvdd_wake_i[m,n],

[0161] where m = int(iv * a_id_iv), n = int(Delta_i_mast_idir + id), that is:

[0162] fvvdd_mast[iv,id] = fvvdd_wake_i[m,n],

[0163] That is, the wind frequency value is transferred to a new wind speed - wind direction interval, and finally the wind frequency matrix fvvdd_wake_i of the actuator disk at the machine position is obtained.

[0164] Optionally, in combination with the above implementation method, in an implementation manner of the above - mentioned embodiment, the measured data of the anemometer tower is obtained from several anemometer towers; determining the wind frequency matrix of the actuator disk at each machine position in each sector of the wind turbine includes: weighted - determining the wind frequency matrix of any one of the machine positions based on the wind frequency matrices of the actuator disks at the corresponding machine positions of several anemometer towers.

[0165] If several anemometer towers around a certain wind turbine participate in the calculation, then the wind frequency matrix of the i - th machine position of the wind turbine can be obtained by weighted operation of the wind frequency matrices respectively deduced by several anemometer towers. This data is more comprehensive and can reflect a more actual situation. If fvvdd_i_j represents that the wind frequency matrix of the i - th machine position is deduced by the j - th anemometer tower, then the wind frequency matrix fvvdd_i of the i - th machine position considering the influence of l anemometer towers is calculated by Formula 9:

[0166]

[0167] In Formula 9, a_i_j represents the weighting factor of the j-th anemometer tower, and the calculation formula is where c j is the confidence coefficient, with a default value of 1, and d j is the straight-line distance between the i-th turbine site and the j-th anemometer tower. ε is taken as 1 to avoid a zero denominator.

[0168] Optionally, in combination with the above implementation method, in another implementation manner of the above embodiment, the operation data of the wind turbine generator set further includes the unit power data of the wind turbine generator set; according to the actuator disk wind frequency matrix of each turbine site, the actuator disk power generation data is determined, including:

[0169] According to the unit power data, determine the power value at the first preset wind speed and the power value at the second preset wind speed in the adjacent wind speed sections of any one turbine site. According to the power value at the first preset wind speed and the power value at the second preset wind speed, obtain the average power of any one turbine position. According to the actuator disk wind frequency matrix and the average power of any one turbine site, determine the actuator disk power generation data.

[0170] Specifically, the calculation formula of the actuator disk wake wind speed in the actuator disk wake parameters can adopt Formula 10:

[0171]

[0172] The actuator disk wake power generation P in the actuator disk wake parameters wake The calculation formula is:

[0173]

[0174] where p(t / 10) is the power value at the t / 10 wind speed, p((t + 1) / 10) is the power value at the (t + 1) / 10 wind speed, the value of t ranges from 0 to 500, which is an index number, and this index number uses integers, while the wind speed may have decimal places, and the wind speed range is 0 to 50 m / s (meters per second), corresponding one-to-one with the index number t. In Formulas 10 and 11, fvvdd_wake[t,:] refers to the actuator disk wind frequency matrix in the t wind speed section, refers to the sum of the actuator disk wind frequency matrices of the entire sector in the t wind speed section.

[0175] Based on the foregoing calculation data and the above formulas, the actuator disk wake parameters that are more in line with the actual power generation data can be obtained.

[0176] Optionally, in combination with the above implementation method, in an implementation manner of the above embodiment, determining the free stream data includes: determining the free stream simulation data according to the free stream model and the operation data of the wind turbine generator set; determining the free stream data including the free stream wind speed data according to the free stream simulation data.

[0177] During the process of calculating the wake parameters of the actuator disk, the free stream wind speed data is required. The free stream wind speed data is corrected using the actuator disk model to obtain more realistic power generation data. The acquisition of the free stream data is based on the free stream model and the operating data of the wind turbine generator set, and through appropriate CFD simulation calculations, the free stream simulation data is first determined.

[0178] Optionally, in the above implementation, determining the free stream simulation data further includes:

[0179] Construct a three-dimensional space grid according to the operating data of the wind turbine generator set. Determine the boundary conditions according to the operating data of the wind turbine generator set, and perform initialization processing on the three-dimensional space grid to obtain an initialized three-dimensional space grid. According to the computational fluid dynamics simulation determination method, perform iterative determination processing on the initialized three-dimensional space grid to obtain the converged free stream simulation data.

[0180] Generate a surface grid according to the terrain data of the wind turbine generator set. Then, according to the machine position point coordinate data and the anemometer tower coordinate data of the wind turbine generator set, increase the grid density of the surface grid to determine the encrypted surface grid. Then stretch the encrypted surface grid in the surface normal direction to determine the free stream three-dimensional space grid. For example, process the terrain data in the original map format to determine the center point and calculation range of the computational domain, and then generate three-dimensional grids according to the number of sectors required for the calculation.

[0181] Add a low roughness value z0 to the three-dimensional space grid of the free stream, and assign the inlet boundary conditions, specifically including: determining the inlet boundary conditions according to the wind profile equation, turbulent kinetic energy equation, turbulent dissipation equation, surface roughness value, reference height data, first reference wind speed, and second reference wind speed. This process can refer to the previous processing method for initializing the three-dimensional space grid of the actuator disk.

[0182] The steps of using the computational fluid dynamics simulation determination method including the open source field operation and operation software package to perform iterative determination processing on the initialized three-dimensional space grid to obtain the converged free stream simulation data include: respectively iteratively determining the initialized three-dimensional space grid to convergence according to the open source field operation and operation software package and the turbulence model to obtain the converged free stream simulation data. The converged free stream simulation data includes wind speed vector data and turbulent kinetic energy data. Among them, the standard k-ε two-equation turbulence model is selected for the turbulence model. A more detailed process can refer to the previous convergence calculation processing of the actuator disk simulation data.

[0183] Optionally, in combination with the above implementation, in some implementations of the first aspect embodiment of the present application, the operation data of the wind turbine includes the unit power data of the wind turbine and the measured data of the anemometry tower. According to the free flow simulation data, the free flow data including the free flow wind speed data is determined, specifically including:

[0184] According to the free flow simulation data and the measured data of the anemometry tower, the free flow wind acceleration factor and the free flow relative wind direction deflection angle are determined.

[0185] According to the measured data of the anemometry tower, the free flow wind acceleration factor and the free flow relative wind direction deflection angle, the free flow wind frequency matrix of each machine position point is determined.

[0186] According to the free flow wind frequency matrix of each machine position point, the free flow wind speed data is determined.

[0187] In the above process, the specific process of determining the free flow wind acceleration factor is as follows:

[0188] At the inlet wind direction dir, the inlet wind direction dir_n of each sector, the inlet wind direction of each degree of wind direction is idir, and the horizontal wind speed V at the machine position point 10_i (the wind speed calculated when the inlet reference wind speed of the i-th machine position is 10 m / s), and the simulated wind speed V under the free flow of the anemometry tower mast (the wind speed calculated when the inlet reference wind speed is 10 m / s), and the free flow wind acceleration factor V of the fan in 16 sectors can be obtained 10_i / V mast . Construct a linear interpolation function with it and the inlet wind direction sequence, and interpolate to obtain the wind acceleration factor a_idir of each idir (i is an integer value in the range of 0 to 359), as shown in Table 4 below:

[0189] Table 4 Relationship between inlet wind direction and wind acceleration factor of machine position point i

[0190] idir 0 1 2 3 … n … 359 a_idir a0 a1 a2 a3 … an … a359

[0191] The determination method of the free flow relative direction deflection angle is the same as the determination method of the relative wind direction deflection angle of each sector in the actuator disk model calculation described above, and will not be elaborated here.

[0192] The measured data of a wind measurement tower is usually a time series (the measured data includes wind speed and wind direction) data. For example, the possible wind speed range of 0 - 50 m / s and the possible direction range of 0 - 359° are partitioned according to a preset wind speed interval of 0.1 m / s and a preset wind direction interval of 1° respectively. The partition frequencies of each partition for wind speed and wind direction are statistically calculated to generate a wind frequency matrix fvvdd_mast (the data volume is 500 * 360). Arbitrarily select a partition from fvvdd_mast. For example, when the wind speed is iv and the wind direction is id, look up the aforementioned Table 1 and Table 3 respectively to obtain the corresponding wind acceleration factor a_id_ and the relative wind direction deflection angle Delta_i_mast_idir. Then the wind speed value at a certain turbine position is iv * a_id, and the wind direction value is Delta_i_mast_idir + id. Perform rounding operations on the wind speed value iv * a_id and the wind direction value Delta_i_mast_idir + id respectively. After the operation, we get:

[0193] fvvdd_i[m,n], where m = int(iv * a_id), n = int(Delta_i_mast_idir + id), that is, fvvdd_mast[iv,id] = fvvdd_i[m,n], which means the wind frequency value is transferred to a new wind speed and wind direction interval, and finally the wind frequency matrix fvvdd_i of the free flow at the turbine position is obtained.

[0194] If there are several wind measurement towers around a certain wind turbine participating in the calculation, then the wind frequency matrix of the i-th turbine position of this wind turbine can be obtained by performing a weighted operation on the wind frequency matrices respectively deduced by several wind measurement towers. This data is more comprehensive and can reflect a more actual situation. fvvdd_i_j represents that the wind frequency matrix of the i-th turbine position is deduced by the j-th wind measurement tower. Then the calculation formula for the wind frequency matrix of the i-th turbine position considering the influence of l wind measurement towers adopts Formula 9 described above.

[0195] Subsequently, the free flow wind speed and the free flow power generation are calculated according to the free flow wind speed calculation formula (Formula 12) and the free flow power generation calculation formula (Formula 13) respectively.

[0196]

[0197]

[0198] Where p(i / 10) is the power value at the wind speed of i / 10, and p((i + 1) / 10) is the power value at the wind speed of (i + 1) / 10. fvvdd_free[t,:] refers to the free flow wind frequency matrix in the t wind speed section, which refers to the sum of the free flow wind frequency matrices of the entire sector in the t wind speed section.

[0199] To comprehensively understand the method concept of this application, the following provides a general description with the accompanying drawings:

[0200] As Figure 4 shown, in an embodiment of this application, the method for determining the wake parameters of a wind turbine mainly consists of two parts: free-stream modeling and simulation comprehensive calculation, and actuator disk modeling and simulation comprehensive calculation. Among them, the free-stream modeling and simulation comprehensive calculation includes three processes: free-stream model modeling, simulation calculation, and comprehensive calculation. The actuator disk modeling and simulation comprehensive calculation also includes three processes: actuator disk modeling, simulation calculation, and comprehensive calculation. And in the process of actuator disk modeling and simulation calculation, the result files obtained from the free-stream modeling and simulation comprehensive calculation are required, such as free-stream wind speed data.

[0201] As Figure 4 shown, the actuator disk modeling and simulation comprehensive calculation first determines the three-dimensional spatial grid of the actuator disk (corresponding to the three-dimensional actuator disk spatial grid in the figure) based on data such as terrain data, machine location coordinates, and wind measurement tower coordinates, in combination with the actuator disk model. Then, in combination with the roughness map, boundary conditions are set and full-field initialization is performed on the three-dimensional spatial grid of the actuator disk, and fluid simulation calculations are performed on all divided sectors, that is, full-sector CFD simulation calculations, to obtain convergent results. Then, the final comprehensive calculation is performed, that is, the actuator disk wake calculation. This calculation process first includes the calculation of the wind acceleration factor, and secondly includes the calculation of the relative wind direction deflection angle. A wake wind frequency matrix is generated based on the previous calculation results. When multiple wind measurement towers are involved, multi-tower calculation is also introduced. Finally, the post-wake wind speed, power generation, and wake reduction at the machine location are calculated, and the calculation results are closer to the actual operating results, and can better evaluate the impact of the wake existing in a large wind farm on power generation loss.

[0202] In the whole calculation process, the free-stream modeling and simulation comprehensive calculation is generally the same as the actuator disk modeling and simulation calculation. However, the calculation models used are different, and there are also significant differences between the two in the calculation of the wind acceleration factor. For detailed information, please refer to the corresponding parts in the previous text. In addition, data contents such as terrain data, machine location coordinates, wind measurement tower coordinates, roughness map, unit power curve, and measured data of the wind measurement tower are used in the process of free-stream modeling and simulation comprehensive calculation and actuator disk modeling and simulation calculation.

[0203] Based on the same inventive concept, an embodiment of the second aspect of this application provides a device 10 for determining the wake parameters of a wind turbine, as Figure 5 shown, including: an acquisition module 11, a simulation module 12, and a determination module 13.

[0204] Among them, the acquisition module 11 is used to acquire wind farm data, the operation data of the wind turbine generator, and the actuator disk model of the impeller of the wind turbine generator. The operation data of the wind turbine generator includes a first reference wind speed and a second reference wind speed.

[0205] The simulation module 12 is used to determine the free flow data according to the wind farm, the operation data, and the actuator disk model, and to determine the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model.

[0206] The determination module 13 is used to determine the wake parameters of the wind turbine generator based on the actuator disk model according to the wind acceleration factor and the relative wind direction deflection angle.

[0207] The wake parameter determination device of the wind turbine generator provided by the present application regards the movement trajectory of the impeller of the wind turbine as an actuator disk model, selects two reference wind speeds in the actuator disk model, and calculates the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake of the actuator disk model according to the wind farm data, the operation data, and the data obtained from the free flow model, and determines the wake parameters of the wind turbine generator, which can ensure better calculation accuracy, reduce the calculation amount at the same time, the calculation result is closer to the actual operation result, and can better evaluate the influence of the wake existing in the large wind farm on the power generation loss.

[0208] Optionally, the free flow data includes free flow wind speed data; the simulation module 12 determines the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model, including:

[0209] Determine the actuator disk simulation data according to the wind farm data, the operation data, and the actuator disk model. Divide the calculation domain corresponding to the impeller in the actuator disk model into several sectors according to several preset wind directions. Determine the wind acceleration factor of each sector according to the actuator disk simulation data, the free flow wind speed data, and the sectors, and determine the relative wind direction deflection angle of each sector.

[0210] Optionally, the simulation module 12 determines the wind acceleration factor of each sector, including: determining a first calibration wind speed and a second calibration wind speed according to the first reference wind speed, the second reference wind speed, and the free flow wind speed data, and dividing the preset wind speed interval into three wind speed sub-intervals; determining the actuator disk wind acceleration factor of each sector according to the wind speed sub-intervals and the free flow wind speed data of each sector.

[0211] Optionally, the wind speed sub - intervals include a first wind speed sub - interval, a second wind speed sub - interval, and a third wind speed sub - interval. The first wind speed sub - interval is the wind speed interval less than the first calibrated wind speed. The second wind speed sub - interval is the wind speed interval greater than or equal to the first calibrated wind speed and less than or equal to the second calibrated wind speed. The third wind speed sub - interval is the wind speed interval greater than the third calibrated wind speed. And, the simulation module 12 determines the actuator disk wind acceleration factor for each sector according to the wind speed sub - intervals and the free - stream wind speed data of each sector, including: if the free - stream wind speed data of a sector is within the first wind speed sub - interval, determining the wind acceleration factor of the sector as the first actuator disk wind acceleration factor; if the free - stream wind speed data of a sector is within the third wind speed sub - interval, determining the wind acceleration factor of the sector as the third actuator disk wind acceleration factor; if the free - stream wind speed data of a sector is within the second wind speed sub - interval, determining the wind acceleration factor of the sector as the second actuator disk wind acceleration factor according to the first calibrated wind speed, the second calibrated wind speed, the first actuator disk wind acceleration factor, and the third actuator disk wind acceleration factor.

[0212] Optionally, the operating data of the wind turbine includes the measured data of the anemometer tower. The simulation module 12 determines the relative wind direction deflection angle of each sector, including:

[0213] Determining the anemometer tower wind direction data, the machine location wind direction data, and the inlet wind direction data of each sector according to the actuator disk simulation data and the measured data of the anemometer tower.

[0214] Determining the wind direction deflection angle of the anemometer tower according to the anemometer tower wind direction data and the inlet wind direction data;

[0215] Determining the anemometer tower wind direction interpolation data according to the first linear interpolation function between the wind direction deflection angle of the anemometer tower of each sector and the inlet wind direction data, and obtaining the first functional relationship between the wind direction deflection angle of the anemometer tower and the anemometer tower wind direction interpolation data.

[0216] Determining the wind direction deflection angle of the machine location according to the machine location wind direction data and the inlet wind direction data.

[0217] Determining the machine location wind direction interpolation data according to the wind direction deflection angle of the machine location and the inlet wind direction data, and obtaining the second functional relationship between the inlet wind direction data and the machine location wind direction interpolation data.

[0218] Determining the relative wind direction deflection angle of each sector according to the first functional relationship and the second functional relationship.

[0219] Optionally, the measured data of the anemometer tower includes wind speed and wind direction. After the simulation module 12 determines the wind acceleration factor of each sector and determines the relative wind direction deflection angle of each sector, it further includes:

[0220] Determine a number of wind speed zones and a number of wind direction zones according to a preset wind speed interval and a preset wind direction interval. Generate a wind frequency matrix based on the partition frequencies corresponding to the wind speed zones and the wind direction zones. Determine the wind speed value and the wind direction value at each machine position corresponding to each element in the wind frequency matrix according to the wind acceleration factor and the relative wind direction deflection angle corresponding to each element in the wind frequency matrix. Round the wind speed value and the wind direction value at each machine position to determine the actuation disk wind frequency matrix of each machine position in each sector. Determine the actuation disk power generation data according to the actuation disk wind frequency matrix of each machine position.

[0221] Optionally, the measured data of the wind measurement tower is obtained from a number of wind measurement towers; the simulation module 12 determines the actuation disk wind frequency matrix of each machine position in each sector, including: determining the wind frequency matrix of any one of the machine positions among the machine positions by weighting according to the actuation disk wind frequency matrix of each machine position corresponding to the number of wind measurement towers.

[0222] Optionally, the operation data of the wind turbine generator set further includes the unit power data of the wind turbine generator set; the determination module 13 determines the actuation disk power generation data according to the actuation disk wind frequency matrix of each machine position, including: determining the power value at the first preset wind speed and the power value at the second preset wind speed of adjacent wind speed segments of any one of the machine positions according to the unit power data. Obtain the average power value of any one of the machine positions according to the power value at the first preset wind speed and the power value at the second preset wind speed. Determine the actuation disk power generation data according to the actuation disk wind frequency matrix and the average power value of any one of the machine positions.

[0223] Optionally, the simulation module 12 determines the actuation disk simulation data according to the wind farm data, the operation data and the actuation disk model, including: constructing an actuation disk three-dimensional space grid according to the operation data of the wind turbine generator set. Determine the boundary conditions according to the wind farm data, the first reference wind speed and the second reference wind speed, and perform an initialization process on the actuation disk three-dimensional space grid to obtain an initialized actuation disk three-dimensional space grid. Perform an iterative determination process on the initialized actuation disk three-dimensional space grid according to the computational fluid dynamics simulation determination method to obtain the converged actuation disk simulation data.

[0224] Optionally, the wind farm data includes the terrain data of the wind turbine generator set, the machine position coordinate data of the wind turbine generator set and the wind measurement tower coordinate data. Constructing an actuation disk three-dimensional space grid according to the operation data of the wind turbine generator set includes: generating a surface grid according to the terrain data of the wind turbine generator set. According to the machine position coordinate data and the wind measurement tower coordinate data of the wind turbine generator set, increase the grid density of the surface grid to determine the encrypted surface grid. Stretch the encrypted surface grid in the surface normal direction to determine the actuation disk three-dimensional space grid.

[0225] Optionally, the operation data of the wind turbine includes the roughness map data corresponding to the wind turbine. The obtaining module 11 obtains the initialized actuator disk three-dimensional space grid, including:

[0226] According to the roughness map data, add the surface roughness value to the grid center point of the actuator disk three-dimensional space grid, and set the inlet boundary condition to obtain the initialized actuator disk three-dimensional space grid respectively.

[0227] Optionally, the operation data of the wind turbine further includes reference height data and several inlet reference wind speed data; the obtaining module 11 sets the inlet boundary condition, including: determining the inlet boundary condition according to the wind profile equation, turbulent kinetic energy equation, turbulent dissipation equation, surface roughness value, reference height data, first reference wind speed and second reference wind speed.

[0228] Optionally, the software toolkit used in the computational fluid dynamics simulation determination method includes an open-source field operation and operation software package. The simulation module 12 performs iterative determination processing on the initialized actuator disk three-dimensional space grid according to the computational fluid dynamics simulation determination method to obtain the converged actuator disk simulation data, including: respectively performing iterative determination on the initialized actuator disk three-dimensional space grid until convergence according to the open-source field operation and operation software package and the turbulence model to obtain the converged actuator disk simulation data. The converged actuator disk simulation data includes wind speed vector data and turbulent kinetic energy data.

[0229] Optionally, the simulation module 12 determines the free stream data, including: determining the free stream simulation data according to the free stream model and the operation data of the wind turbine. According to the free stream simulation data, determine the free stream data including the free stream wind speed data.

[0230] Optionally, the simulation module 12 determines the free stream simulation data, including: constructing a three-dimensional space grid according to the operation data of the wind turbine; determining the boundary condition according to the operation data of the wind turbine and performing initialization processing on the three-dimensional space grid to obtain the initialized three-dimensional space grid; performing iterative determination processing on the initialized three-dimensional space grid according to the computational fluid dynamics simulation determination method to obtain the converged free stream simulation data.

[0231] Optionally, the operation data of the wind turbine includes the unit power data of the wind turbine and the measured data of the anemometer tower. The simulation module 12 determines the free stream data including the free stream wind speed data according to the free stream simulation data, including: determining the free stream wind acceleration factor and the free stream relative wind direction deflection angle according to the free stream simulation data and the measured data of the anemometer tower. According to the measured data of the anemometer tower, the free stream wind acceleration factor and the free stream relative wind direction deflection angle, determine the free stream wind frequency matrix of each machine position point. According to the free stream wind frequency matrix of each machine position point, determine the free stream wind speed data.

[0232] Based on the same inventive concept, an embodiment of the third aspect of the present application provides a wake parameter determination device for a wind turbine generator set, and the device includes:

[0233] a processor and a memory; the memory is electrically connected to the processor;

[0234] at least one program, which is stored in the memory and configured to be executed by the processor, and the at least one program is configured to: implement the wake parameter determination method for the wind turbine generator set described in the embodiment of the first aspect of the present application.

[0235] Those skilled in the art of the present technology can understand that the electronic device provided in the embodiment of the present application can be specially designed and manufactured for the required purpose, or can also include known devices in a general-purpose computer. These devices have computer programs stored therein, and these computer programs are selectively activated or reconstructed. Such computer programs can be stored in a device (such as a computer) readable medium or stored in any type of medium suitable for storing electronic instructions and coupled to the bus respectively.

[0236] Compared with the prior art, it can be realized that: by regarding the impeller movement trajectory of the wind turbine as an actuator disk model, two reference wind speeds are selected in the actuator disk model, and according to the wind farm data, operation data and the data obtained from the free flow model, the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake of the actuator disk model are calculated to determine the wake parameters of the wind turbine generator set, which can ensure better calculation accuracy, while reducing the calculation amount, the calculation result is closer to the actual operation result, and it can better evaluate the influence of the wake existing in a large-scale wind farm on the power generation loss.

[0237] The present application provides a wake parameter determination device in an alternative embodiment, as Figure 6 shown Figure 6 The wake parameter determination device 1000 shown includes a processor 1001 and a memory 1003. Among them, the processor 1001 and the memory 1003 are electrically connected, such as connected through a bus 1002.

[0238] The processor 1001 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field-Programmable Gate Array), or other programmable logic devices, transistor logic devices, hardware components, or any combination thereof. It can implement or execute various exemplary logical blocks, modules, and circuits described in connection with the disclosure of this application. The processor 1001 may also be a combination that implements computing functions, such as a combination of one or more microprocessors, a combination of a DSP and a microprocessor, etc.

[0239] The bus 1002 may include a path for transmitting information between the above components. The bus 1002 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus, etc. The bus 1002 may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 6 only a thick line is used to represent it in the figure, but it does not mean that there is only one bus or one type of bus.

[0240] The memory 1003 may be a ROM (Read-Only Memory) or other type of static storage device that can store static information and instructions, a RAM (random access memory), or other type of dynamic storage device that can store information and instructions. It may also be an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read-Only Memory), or other optical disc storage, optical disc storage (including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), magnetic storage media, or other magnetic storage devices, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto.

[0241] Optionally, the wake parameter determination device 1000 may further include a transceiver 1004. The transceiver 1004 can be used for receiving and transmitting signals. The transceiver 1004 may allow the wake parameter determination device 1000 to communicate with other devices wirelessly or wiredly to exchange data. It should be noted that in practical applications, the transceiver 1004 is not limited to one.

[0242] Optionally, the wake parameter determination device 1000 may further include an input unit 1005. The input unit 1005 can be used for receiving input digital, character, image, and / or sound information, or generating key signal inputs related to the user settings and function controls of the wake parameter determination device 1000. The input unit 1005 may include, but is not limited to, one or more of a touch screen, a physical keyboard, function keys (such as volume control keys, switch keys, etc.), a trackball, a mouse, a joystick, a photographing device, a pickup, etc.

[0243] Optionally, the wake parameter determination device 1000 may further include an output unit 1006. The output unit 1006 can be used for outputting or presenting the information processed by the processor 1001. The output unit 1006 may include, but is not limited to, one or more of a display device, a speaker, a vibration device, etc.

[0244] Although Figure 6 the wake parameter determination device 1000 with various devices is shown, it should be understood that it is not required to implement or have all the shown devices. More or fewer devices may be alternatively implemented or had.

[0245] Optionally, the memory 1003 is used for storing the application program code for executing the solution of this application, and is controlled by the processor 1001 to execute. The processor 1001 is used for executing the application program code stored in the memory 1003 to implement any one of the wake parameter determination methods of the wind turbine generator set provided in the embodiments of this application.

[0246] Based on the same inventive concept, the embodiments of this application provide a computer-readable storage medium. The computer-readable storage medium stores at least one segment of program or code set, and the at least one segment of program or code set is loaded and executed by a processor to implement the wake parameter determination method of the wind turbine generator set provided in the first aspect embodiment of this application.

[0247] The computer-readable storage medium provided by the embodiments of the present application, by executing the method for determining the wake parameters of a wind turbine provided by the embodiments of the present application, regards the impeller movement trajectory of the wind turbine as an actuator disk model, selects two reference wind speeds in the actuator disk model, and calculates the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake in the actuator disk model according to the wind farm data, operation data, and data obtained from the free flow model, and determines the wake parameters of the wind turbine, which can ensure better calculation accuracy, while reducing the calculation amount, the calculation result is closer to the actual operation result, and can better evaluate the influence of the wake existing in a large wind farm on the power generation loss.

[0248] Those skilled in the art of the present technology can understand that the various operations, methods, steps, measures, and solutions in the processes discussed in the present application can be alternated, changed, combined, or deleted. Further, other steps, measures, and solutions in the various operations, methods, and processes discussed in the present application can also be alternated, changed, rearranged, decomposed, combined, or deleted. Further, the steps, measures, and solutions in the prior art that are the same as those disclosed in the present application can also be alternated, changed, rearranged, decomposed, combined, or deleted.

[0249] The terms "first" and "second" are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include one or more of such features. In the description of the present application, unless otherwise stated, the meaning of "a plurality" is two or more.

[0250] It should be understood that although the steps in the flowchart of the accompanying drawings are shown in sequence according to the indication of the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless there is a clear indication in this article, the execution of these steps is not strictly limited in order, and they can be executed in other orders. Moreover, at least a part of the steps in the flowchart of the accompanying drawings may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be executed alternately or alternately with at least a part of other steps or sub-steps or stages of other steps.

[0251] The above are only some embodiments of the present application. It should be noted that for those of ordinary skill in the art of the present technology, without departing from the principle of the present application, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present application.

Claims

1. A method for determining wake parameters of a wind turbine generator, characterized in that Including: Obtaining wind farm data, operating data of a wind turbine generator set, and an actuator disk model of an impeller of the wind turbine generator set, wherein the operating data of the wind turbine generator set includes a first reference wind speed and a second reference wind speed; Determining free flow data according to the wind farm data, the operating data, and the actuator disk model, and determining a wind acceleration factor and a relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model; Determining wake parameters of the wind turbine generator set based on the actuator disk model according to the wind acceleration factor and the relative wind direction deflection angle; The free flow data includes free flow wind speed data; The determining the wind acceleration factor and the relative wind direction deflection angle of each sector of the impeller wake based on the actuator disk model includes: Determining actuator disk simulation data according to the wind farm data, the operating data, and the actuator disk model; Dividing a calculation domain corresponding to the impeller in the actuator disk model into a plurality of sectors according to a plurality of preset wind directions; Determining the wind acceleration factor of each sector and determining the relative wind direction deflection angle of each sector according to the actuator disk simulation data, the free flow wind speed data, and the sectors; The determining the free flow data includes: Determining free flow simulation data according to a free flow model and the operating data of the wind turbine generator set; Determining the free flow data including free flow wind speed data according to the free flow simulation data; The operating data of the wind turbine generator set includes unit power data of the wind turbine generator set and measured data of a wind measurement tower; The determining the free flow data including free flow wind speed data according to the free flow simulation data includes: Determining a free flow wind acceleration factor and a free flow relative wind direction deflection angle according to the free flow simulation data and the measured data of the wind measurement tower; Determining a free flow wind frequency matrix of each machine position point according to the measured data of the wind measurement tower, the free flow wind acceleration factor, and the free flow relative wind direction deflection angle; Determining free flow wind speed data according to the free flow wind frequency matrix of each machine position point; 2. The method for determining the wake parameters of a wind turbine generator according to claim 1, wherein The determining the wind acceleration factor of each sector includes: Determining a first calibration wind speed and a second calibration wind speed according to the first reference wind speed, the second reference wind speed, and the free flow wind speed data, and dividing a preset wind speed interval into three wind speed sub-intervals; Determining the actuator disk wind acceleration factor of each sector according to the wind speed sub-intervals and the free flow wind speed data of each sector; 3. The method for determining the wake parameters of a wind turbine generator according to claim 2, wherein The wind speed sub-intervals include a first wind speed sub-interval, a second wind speed sub-interval, and a third wind speed sub-interval. The first wind speed sub-interval is a wind speed interval less than the first calibration wind speed, the second wind speed sub-interval is a wind speed interval greater than or equal to the first calibration wind speed and less than or equal to the second calibration wind speed, and the third wind speed sub-interval is a wind speed interval greater than a third calibration wind speed; And, the determining the actuator disk wind acceleration factor of each sector according to the wind speed sub-intervals and the free flow wind speed data of each sector includes: If the free flow wind speed data of a sector is within the first wind speed sub-interval, determining the wind acceleration factor of the sector as a first actuator disk wind acceleration factor; If the free - stream wind speed data of a sector is within the third wind speed sub - interval, determine that the wind acceleration factor of the sector is the third actuator disk wind acceleration factor; If the free - stream wind speed data of a sector is within the second wind speed sub - interval, determine that the wind acceleration factor of the sector is the second actuator disk wind acceleration factor according to the first calibrated wind speed, the second calibrated wind speed, the first actuator disk wind acceleration factor, and the third actuator disk wind acceleration factor.

4. The method for determining the wake parameters of a wind turbine generator according to claim 1, wherein The operation data of the wind turbine generator set includes the measured data of the anemometer tower; determining the relative wind direction deflection angle of each sector includes: Determine the anemometer tower wind direction data, the machine position point wind direction data, and the inlet wind direction data of each sector according to the actuator disk simulation data and the measured data of the anemometer tower; Determine the wind direction deflection angle of the anemometer tower according to the anemometer tower wind direction data and the inlet wind direction data; Determine the anemometer tower wind direction interpolation data according to the first linear interpolation function between the wind direction deflection angle of the anemometer tower of each sector and the inlet wind direction data, and obtain the first functional relationship between the wind direction deflection angle of the anemometer tower and the anemometer tower wind direction interpolation data; Determine the wind direction deflection angle of the machine position point according to the machine position point wind direction data and the inlet wind direction data; Determine the machine position point wind direction interpolation data according to the wind direction deflection angle of the machine position point and the inlet wind direction data, and obtain the second functional relationship between the inlet wind direction data and the machine position point wind direction interpolation data; Determine the relative wind direction deflection angle of each sector according to the first functional relationship and the second functional relationship.

5. The method for determining the wake parameters of a wind turbine according to claim 4, characterized in that, The measured data of the anemometer tower includes wind speed and wind direction; after the steps of determining the wind acceleration factor of each sector and determining the relative wind direction deflection angle of each sector, it further includes: Determine a number of wind speed sub - regions and a number of wind direction sub - regions according to the preset wind speed interval and the preset wind direction interval; Generate a wind frequency matrix according to the partition frequencies corresponding to the wind speed sub - regions and the wind direction sub - regions; Determine the wind speed value and the wind direction value at the machine position corresponding to each element in the wind frequency matrix according to the wind acceleration factor and the relative wind direction deflection angle corresponding to each element in the wind frequency matrix; Round the wind speed value and the wind direction value at the machine position to determine the actuator disk wind frequency matrix of each machine position point in each sector; Determine the actuator disk power generation data according to the actuator disk wind frequency matrix of each machine position point.

6. The method for determining the wake parameters of a wind turbine according to claim 5, wherein The measured data of the anemometer tower is obtained from a number of anemometer towers; determining the actuator disk wind frequency matrix of each machine position point in each sector includes: Weightedly determine the wind frequency matrix of any one machine position point among each machine position point according to the actuator disk wind frequency matrix of each machine position point corresponding to a number of anemometer towers.

7. The method for determining the wake parameters of a wind turbine according to claim 5, characterized in that, The operation data of the wind turbine generator set further includes the unit power data of the wind turbine generator set; The determining the actuator disk power generation data according to the actuator disk wind frequency matrix of each machine position point includes: Determine the power value at the first preset wind speed and the power value at the second preset wind speed of adjacent wind speed segments of any one machine position point according to the unit power data; Obtain the average power value of any wind turbine position based on the power value at the first preset wind speed and the power value at the second preset wind speed; Determine the power generation data of the actuator disk based on the actuator disk wind frequency matrix and the average power value of any wind turbine position.

8. The method for determining the wake parameters of a wind turbine generator according to claim 1, characterized in that, Determine the actuator disk simulation data based on the wind farm data, the operation data, and the actuator disk model, including: Construct a three-dimensional spatial grid of the actuator disk according to the operation data of the wind turbine generator; Determine the boundary conditions based on the wind farm data, the first reference wind speed, and the second reference wind speed, and perform initialization processing on the three-dimensional spatial grid of the actuator disk to obtain an initialized three-dimensional spatial grid of the actuator disk; Perform iterative determination processing on the initialized three-dimensional spatial grid of the actuator disk according to the computational fluid dynamics simulation determination method to obtain the converged actuator disk simulation data.

9. The method for determining the wake parameters of a wind turbine according to claim 8, characterized in that, The wind farm data includes the terrain data of the wind turbine generator, the coordinates data of the wind turbine positions, and the coordinates data of the wind measurement towers; Construct a three-dimensional spatial grid of the actuator disk according to the operation data of the wind turbine generator, including: Generate a surface grid according to the terrain data of the wind turbine generator; Increase the grid density of the surface grid according to the coordinates data of the wind turbine positions and the coordinates data of the wind measurement towers to determine the encrypted surface grid; Perform surface normal stretching on the encrypted surface grid to determine the three-dimensional spatial grid of the actuator disk.

10. The method for determining the wake parameters of a wind turbine according to claim 8, characterized in that, The operation data of the wind turbine generator includes the roughness map data corresponding to the wind turbine generator; The step of obtaining the initialized three-dimensional spatial grid of the actuator disk includes: Add the surface roughness value to the grid center point of the three-dimensional spatial grid of the actuator disk according to the roughness map data, and set the inlet boundary conditions to obtain the initialized three-dimensional spatial grid of the actuator disk respectively.

11. The method for determining the wake parameters of a wind turbine according to claim 10, characterized in that, The operation data of the wind turbine generator further includes reference height data and a plurality of inlet reference wind speed data; The setting of the inlet boundary conditions includes: Determine the inlet boundary conditions according to the wind profile equation, the turbulent kinetic energy equation, the turbulent dissipation equation, the surface roughness value, the reference height data, the first reference wind speed, and the second reference wind speed.

12. The method for determining the wake parameters of a wind turbine according to claim 8, characterized in that, The software toolkit used in the computational fluid dynamics simulation determination method includes an open-source field operation and operation software toolkit; The step of performing iterative determination processing on the initialized three-dimensional spatial grid of the actuator disk according to the computational fluid dynamics simulation determination method to obtain the converged actuator disk simulation data includes: iteratively determining the initialized three-dimensional spatial grid of the actuator disk to convergence respectively according to the open-source field operation and operation software toolkit and the turbulence model to obtain the converged actuator disk simulation data; The converged actuator disk simulation data includes wind speed vector data and turbulent kinetic energy data.

13. The method for determining the wake parameters of a wind turbine according to claim 1, characterized in that, The determination of the free flow simulation data includes: Construct a three-dimensional spatial grid according to the operation data of the wind turbine generator; Determine the boundary conditions according to the operation data of the wind turbine generator, and perform initialization processing on the three-dimensional spatial grid to obtain an initialized three-dimensional spatial grid; According to the computational fluid dynamics simulation determination method, perform iterative determination processing on the initialized three-dimensional space grid to obtain the converged free flow simulation data.

14. A wake parameter determination device for a wind turbine, characterized in that a processor; a memory electrically connected to the processor; at least one program stored in the memory and configured to be executed by the processor, the at least one program being configured to: implement the wake parameter determination method for a wind turbine as described in any one of claims 1 to 13.

15. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores at least one program or code set, and the at least one program or the code set is loaded and executed by a processor to implement the wake parameter determination method for a wind turbine as described in any one of claims 1 to 13.

Images