A wind farm wake evaluation method, device, equipment and storage medium

Abstract

The application discloses a kind of wind farm wake evaluation method, device, equipment and storage medium, method includes: constructing physical model and generating coarse grid and fine grid;First calculation result and first flow field are obtained by CFD solution to coarse grid, and first flow field is as the initial flow field of fine grid;Actuator disc parameter is determined according to first calculation result;According to boundary condition, initial flow field and actuator disc parameter, the CFD solution based on actuator disc is carried out to fine grid to obtain second calculation result;Second machine site wind direction is constantly corrected according to second calculation result in second machine site wind direction Actuator disc orientation, until the difference between actuator disc orientation and second machine site wind direction is less than preset threshold, stop correction and obtain the wake parameter in each machine site and preset downstream range.The application corrects actuator disc orientation according to machine site wind direction, can more truly consider the influence of complex terrain and upstream wind turbine wake to downstream wind turbine, improve the accuracy and effectiveness of wind farm wake evaluation.

Description

Technical Field

[0001] This invention relates to the field of wind resource assessment technology, and in particular to a method, apparatus, equipment and storage medium for assessing wind farm wake. Background Technology

[0002] Wind farm wake assessment is a crucial component of wind resource assessment, significantly impacting turbine selection, layout, and overall investment efficiency. Currently, commonly used wake assessment methods can be broadly categorized into two types: linear engineering wake models and nonlinear computational fluid dynamics (CFD) models. Engineering wake models offer high solution speeds due to their computational simplicity, but the results heavily rely on empirical parameters, leading to poor accuracy and adaptability to terrain, especially complex mountainous terrain. In recent years, actuated disk-based CFD methods have gained widespread application in wind farm wake assessment due to their high computational accuracy, accurate representation of physical phenomena, and consideration for computational resources.

[0003] To implement CFD calculations based on actuator disks, it is necessary to configure actuator disk source terms and determine the regions where these source terms are added based on wind turbine parameters. In traditional actuator disk-based CFD calculations, the actuator disk orientation is directly specified based on the sector direction. In flat terrain, the wind direction at each turbine location is generally consistent with the sector direction, making this method suitable for flat terrain. However, in complex terrain, the wind direction at the turbine location is often deflected by the terrain, and the yaw strategy ensures that the wind turbine remains aligned with the wind. In this case, specifying the actuator disk orientation based on the sector direction will lead to errors in the magnitude and direction of the calculated forces, failing to fully account for the wake effect, thus reducing the accuracy and effectiveness of wind farm wake assessment. Summary of the Invention

[0004] In view of this, the purpose of the present invention is to provide a method, apparatus, equipment and storage medium for wind farm wake assessment, which solves the problem of inaccurate wind farm wake assessment in the prior art under complex terrain conditions.

[0005] To address the aforementioned technical problems, this invention provides a method for assessing wind farm wake, comprising:

[0006] Acquire measured topographic maps, satellite maps, projection systems, and camera location information, construct a physical model of the target terrain, and generate coarse and fine grids;

[0007] Obtain boundary conditions, and perform CFD solution on the coarse mesh to obtain the first calculation result and the first flow field, and use the first flow field as the initial flow field of the fine mesh;

[0008] Based on the first calculation result, the actuator disk parameters are determined, and the fine mesh is topologically mapped and the source terms are configured.

[0009] Based on the boundary conditions, the initial flow field, and the actuation disk parameters, the fine mesh is solved using CFD based on the actuation disk to obtain the second calculation result and the second flow field.

[0010] Based on the location information, the wind direction of the second location is extracted from the second calculation result;

[0011] The orientation of the actuator disk is continuously corrected based on the wind direction at the second machine point. When the difference between the orientation of the actuator disk and the wind direction at the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within a preset downstream range are obtained.

[0012] Optionally, obtaining the boundary conditions includes:

[0013] The boundary conditions are calculated based on roughness, thermal stability, and reference wind speed.

[0014] The roughness is obtained based on the roughness map data;

[0015] The inlet boundary conditions are obtained by solving the wind profile equation based on the Monin-Obukhov similarity theory.

[0016] Optionally, determining the actuator disk parameters based on the first calculation result and performing topology mapping and source term configuration on the fine mesh includes:

[0017] Based on the machine location information, the wind direction of the first machine location is extracted from the first calculation result, and the orientation of the actuation disk is determined based on the wind direction of the first machine location;

[0018] The configuration area of ​​the source item is determined by performing topological mapping on the fine mesh according to the wind turbine size;

[0019] The configuration parameters of the source term are determined based on the model curve, the fan size, and the wind direction at the first fan location.

[0020] Optionally, the step of performing CFD solution based on the actuated disk on the fine mesh to obtain the second calculation result and the second flow field includes:

[0021] The second calculation result is obtained by performing CFD solution based on the actuated disk on the fine mesh;

[0022] The actuator disk source terms are calculated by CFD solution based on the actuator disk parameters to obtain the second flow field.

[0023] The formula for calculating the actuation disk source term is as follows:

[0024]

[0025]

[0026] Among them, S u The actuation disk source term is represented by ρ, air density by Δx, actuation disk thickness by u1, and inflow velocity by C. x Let represent the induction coefficient, and 'a' represent the induction factor. The induction factor 'a' is expressed as:

[0027]

[0028] Where C P C represents the power factor. T This represents the thrust coefficient.

[0029] Optionally, the step of continuously correcting the orientation of the actuator disk based on the wind direction at the second machine location, and stopping the correction when the difference between the orientation of the actuator disk and the wind direction at the second machine location is less than a preset threshold, and obtaining the wake parameters of each machine location and the wake parameters within a preset downstream range, includes:

[0030] If the difference between the orientation of the actuator disk and the wind direction of the second machine point is greater than the preset threshold, then the wind direction of the second machine point is taken as the orientation of the actuator disk, and the steps of performing topology mapping and source term configuration on the fine mesh are executed.

[0031] If the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than the preset threshold, the step of stopping correction and obtaining the wake parameters of each machine point and the wake parameters within the preset downstream range is executed.

[0032] Optionally, the step of acquiring measured topographic maps, satellite maps, projection systems, and camera location information, constructing a physical model of the target terrain, and generating coarse and fine grids includes:

[0033] The measured topographic map is processed according to the projection system and the satellite map to obtain the physical model;

[0034] The coarse mesh is obtained using the physical model and meshing tools described above;

[0035] The entire field mesh is refined based on the coarse mesh, and the fine mesh is obtained by locally refining the mesh at the machine points according to the machine point information and the mesh tool.

[0036] Optionally, the first flow field, as the initial flow field of the fine mesh, includes:

[0037] The first flow field is mapped onto the fine mesh using a three-dimensional interpolation method, so that the first flow field serves as the initial flow field of the fine mesh.

[0038] The present invention also provides a wind farm wake assessment device, comprising:

[0039] The acquisition and construction module is used to acquire measured topographic maps, satellite maps, projection systems and camera position information, construct a physical model of the target terrain, and generate coarse and fine grids;

[0040] The acquisition and solution module is used to acquire boundary conditions and perform CFD solution on the coarse mesh to obtain a first calculation result and a first flow field, and use the first flow field as the initial flow field of the fine mesh;

[0041] The determination module is used to determine the actuator disk parameters based on the first calculation result, and to perform topology mapping and source term configuration on the fine mesh;

[0042] The solver module is used to perform CFD solution based on the actuation disk on the fine mesh according to the boundary conditions, the initial flow field and the actuation disk parameters to obtain the second calculation result and the second flow field;

[0043] The extraction module is used to extract the wind direction of the second machine location from the second calculation result based on the machine location information;

[0044] The correction module is used to continuously correct the orientation of the actuator disk according to the wind direction of the second machine point. When the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within a preset downstream range are obtained.

[0045] The present invention also provides a wind farm wake assessment device, comprising:

[0046] Memory, used to store computer programs;

[0047] A processor is used to implement the steps of the wind farm wake assessment method described above when executing the computer program.

[0048] The present invention also provides a storage medium storing a computer program, which, when executed by a processor, implements the steps of the wind farm wake assessment method described above.

[0049] As can be seen, this invention acquires measured topographic maps, satellite maps, projection systems, and aircraft location information, constructs a physical model of the target terrain, and generates coarse and fine grids; acquires boundary conditions, and performs CFD solution on the coarse grid to obtain a first calculation result and a first flow field, using the first flow field as the initial flow field for the fine grid; determines the actuator disk parameters based on the first calculation result, and performs topology mapping and source term configuration on the fine grid; performs CFD solution based on the actuator disk on the fine grid based on the boundary conditions, initial flow field, and actuator disk parameters to obtain a second calculation result and a second flow field; extracts the second aircraft location wind direction from the second calculation result based on the aircraft location information; continuously corrects the actuator disk orientation based on the second aircraft location wind direction, and stops correction when the difference between the actuator disk orientation and the second aircraft location wind direction is less than a preset threshold, obtaining the wake parameters of each aircraft location and the wake parameters within a preset downstream range. This invention corrects the orientation of the actuator disk based on the wind direction at the turbine location, which can more realistically take into account the impact of complex terrain and the wake of upstream wind turbines on downstream wind turbines, improve the accuracy and effectiveness of wind farm wake assessment, overcome the problems of poor adaptability to complex terrain and low computational efficiency in wind farm wake assessment, effectively realize the refined assessment of wind farm wake, and provide a basis for wind farm power generation assessment and turbulence assessment, which is conducive to improving the reliability of wind resource assessment.

[0050] In addition, the present invention also provides a wind farm wake assessment device, equipment and storage medium, which also have the above-mentioned beneficial effects. Attached Figure Description

[0051] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0052] Figure 1 A flowchart of a wind farm wake assessment method provided in an embodiment of the present invention;

[0053] Figure 2 A flowchart illustrating a wind farm wake assessment method provided in an embodiment of the present invention;

[0054] Figure 3 This is a velocity distribution diagram of the wake region at location P2 provided in an embodiment of the present invention;

[0055] Figure 4 This is a schematic diagram of the structure of a wind farm wake assessment device provided in an embodiment of the present invention;

[0056] Figure 5This is a schematic diagram of a wind farm wake assessment device provided in an embodiment of the present invention. Detailed Implementation

[0057] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0058] Please refer to Figure 1 , Figure 1 A flowchart illustrating a wind farm wake assessment method provided in this embodiment of the invention. The method may include:

[0059] S101: Acquire measured topographic maps, satellite maps, projection systems, and camera location information, construct a physical model of the target terrain, and generate coarse and fine grids.

[0060] The execution subject in this embodiment is a terminal. This embodiment does not limit the type of terminal, as long as it can perform the wind farm wake assessment method described in this embodiment. For example, it can be a general-purpose terminal or a dedicated terminal. Based on the measured topographic map, projection system, satellite map, and machine location information, a physical model is constructed. Two sets of meshes are generated using the physical model, including a coarse mesh and a fine mesh.

[0061] Furthermore, to ensure the accuracy of fine grid division, the above-mentioned acquisition of measured topographic maps, satellite maps, projection systems, and camera position information, construction of a physical model of the target terrain, and generation of coarse and fine grids may include the following steps, specifically:

[0062] Step 61: Process the measured topographic map according to the projection system and satellite map to obtain a physical model.

[0063] Step 62: Obtain a coarse mesh using the physical model and meshing tools.

[0064] Step 63: Based on the coarse grid, refine the entire field grid, and locally densify the grid at the machine location points according to the machine location information and the grid tool to obtain a fine grid.

[0065] The obtained measured topographic map is processed based on the obtained projection system and satellite map to obtain a physical model. A coarse grid is generated on the physical model using a meshing tool. Based on the generated coarse grid, the entire field grid is refined. Furthermore, based on the location information of the aircraft, the meshing tool is used to locally densify the grid near the aircraft location to obtain a fine grid.

[0066] S102: Obtain boundary conditions and perform CFD solution on the coarse mesh to obtain the first calculation result and the first flow field, and use the first flow field as the initial flow field of the fine mesh.

[0067] Based on the obtained boundary conditions, CFD is performed on the coarse grid, the calculation results at the machine location are extracted as the first calculation result, and the terrain flow field under the condition of no wind turbine is obtained as the first flow field. The first flow field (the terrain flow field under the condition of no wind turbine) is used as the initial flow field of the fine grid.

[0068] Furthermore, to improve the accuracy of the calculation results, obtaining the boundary conditions described above may include the following steps, specifically:

[0069] Step 21: Calculate the boundary conditions based on roughness, thermal stability and reference wind speed.

[0070] Step 22: Roughness is obtained from roughness map data.

[0071] Step 23: The inlet boundary conditions are obtained by solving the wind profile equation based on the Monin-Obukhov similarity theory.

[0072] This embodiment does not limit the number or content of boundary conditions. For example, boundary conditions may include wall boundary conditions; or they may include symmetrical boundary conditions; or they may include inlet boundary conditions. Boundary conditions can be any of the above conditions or any combination of the above conditions. The boundary conditions are calculated based on information such as roughness, thermal stability, and reference wind speed, wherein the roughness is obtained based on roughness map data.

[0073] In this embodiment, the boundary conditions may include inlet boundary conditions. This embodiment does not limit the method of obtaining the inlet boundary. In this embodiment, the inlet boundary conditions are obtained by solving the wind profile equation based on the Monin-Obukhov similarity theory (the Monin-Obukhov similarity theory describes the dimensionless average flow and average temperature in the surface layer under non-neutral conditions as functions of dimensionless height parameters). The calculation formula is as follows:

[0074]

[0075]

[0076] Where U(z) represents the velocity distribution with height, z represents the height above the ground, and u * κ represents the frictional wind speed, z0 represents the von Kármán constant, L represents the roughness, and L represents the Obuhoff length.

[0077] S103: Determine the actuator disk parameters based on the first calculation result, and perform topology mapping and source term configuration on the fine mesh.

[0078] The first calculation result (the calculation result of the turbine location under the condition of no wind turbine) is used as the actuation disk parameter to perform topology mapping and resource configuration on the fine grid respectively.

[0079] Furthermore, accurately performing topology mapping and source term configuration on the mesh, as described above regarding determining the actuator disk parameters based on the first calculation result and performing topology mapping and source term configuration on the fine mesh, may include the following steps:

[0080] Step 31: Extract the wind direction of the first machine location from the first calculation result based on the machine location information, and determine the orientation of the actuation disk based on the wind direction of the first machine location.

[0081] Step 32: Perform topology mapping on the fine mesh according to the wind turbine size to determine the configuration area of ​​the source item.

[0082] Step 33: Determine the configuration parameters of the source term based on the model curve, fan size, and wind direction at the first fan location.

[0083] The physical parameters of the turbine location are extracted from the first flow field (the topographic flow field under turbine-free conditions) based on the turbine location information. These parameters may include wind direction and wind speed. Specifically, the wind direction extracted from the first calculation result (the calculation result of the turbine location under turbine-free conditions) is used to determine the initial actuation disk orientation. When performing topology mapping on the fine mesh, the configuration area of ​​the source term can be determined based on the turbine size. The configuration parameters of the source term are determined based on the turbine curve, turbine size, and the wind direction of the first turbine location, and then configured accordingly.

[0084] S104: Based on the boundary conditions, initial flow field, and actuation disk parameters, perform CFD solution based on the actuation disk on the fine mesh to obtain the second calculation result and the second flow field.

[0085] Based on the boundary conditions, initial flow field, and actuation disk parameters, CFD calculations based on the actuation disk are performed on the fine mesh to obtain the second calculation result (the calculation result of the turbine location under the condition of having a fan) and the second flow field (the terrain flow field under the condition of having a fan).

[0086] Furthermore, to improve the computational efficiency of the actuation disk source terms, the above-mentioned CFD solution based on the actuation disk for the fine mesh to obtain the second calculation result and the second flow field may include the following steps:

[0087] The second calculation result is obtained by performing CFD solution based on the actuated disk on the fine mesh;

[0088] The actuation disk source terms are calculated by performing CFD solution based on the actuation disk on the fine mesh according to the actuation disk parameters to obtain the second flow field;

[0089] The formula for calculating the actuator disk source term is as follows:

[0090]

[0091]

[0092] Among them, S u The actuation disk source term is represented by ρ, air density by Δx, actuation disk thickness by u1, and inflow velocity by C. x Let represent the induction coefficient, and 'a' represent the induction factor. The induction factor 'a' is expressed as:

[0093]

[0094] Where C P C represents the power factor. T This represents the thrust coefficient.

[0095] The model curves include power curves and thrust curves. The power coefficient and thrust coefficient are obtained from the model curves to obtain the induction factor. Then, based on the actuator disk parameters, the fine mesh is calculated using CFD based on the actuator disk to obtain the actuator disk source term, thereby obtaining the second flow field (the terrain flow field under the condition of having a wind turbine).

[0096] S105: Extract the wind direction of the second machine location from the second calculation result based on the machine location information.

[0097] The physical parameters, including wind direction and wind speed, are obtained from the second calculation result (the calculation result under the condition of having a wind turbine) based on the turbine location information. The wind direction is then determined as the wind direction at the second turbine location.

[0098] S106: The direction of the actuator disk is continuously corrected according to the wind direction of the second machine point. When the difference between the direction of the actuator disk and the wind direction of the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within the preset downstream range are obtained.

[0099] The actuator disk orientation is continuously corrected based on the wind direction at the second machine location until the difference between the actuator disk orientation and the wind direction at the second machine location is less than a preset threshold. At this point, the wake parameters for each machine location and the wake parameters within a preset downstream range are obtained. This embodiment does not limit the preset threshold and preset downstream range; users can set them according to their actual needs.

[0100] Furthermore, to ensure the accuracy of the correction results, the aforementioned continuous correction of the actuator disk orientation based on the wind direction at the second machine location, stopping the correction when the difference between the actuator disk orientation and the wind direction at the second machine location is less than a preset threshold, and obtaining the wake parameters for each machine location and the wake parameters within a preset downstream range, may include the following steps, specifically:

[0101] Step 51: If the difference between the orientation of the actuator disk and the wind direction of the second machine point is greater than a preset threshold, then the wind direction of the second machine point is taken as the orientation of the actuator disk, and the steps of topology mapping and source term configuration for the fine grid are executed.

[0102] Step 52: If the difference between the direction of the actuator disk and the wind direction of the second machine point is less than the preset threshold, execute the step of stopping correction and obtaining the wake parameters of each machine point and the wake parameters within the preset downstream range.

[0103] When the difference between the wind direction at the second turbine location and the orientation of the actuator disk is greater than a preset threshold, the wind direction at the second turbine location is used as the orientation of the actuator disk, overwriting the previous orientation value of the actuator disk, and the topology mapping and source term configuration of the fine grid and subsequent steps in S103 are continued; when the difference between the wind direction at the second turbine location and the orientation of the actuator disk is greater than or less than a preset threshold, the correction is stopped, and the wake parameters of each turbine location and the wake parameters within a preset downstream range are obtained for subsequent comprehensive calculations. This embodiment does not limit the subsequent actions, for example, it can be used for power generation calculation; or it can be used for turbulence calculation.

[0104] Furthermore, based on any of the above embodiments, to ensure the accuracy of the initial flow field and the convergence of subsequent calculations based on the actuated disk CFD, the first flow field, as the initial flow field of the fine mesh, may include:

[0105] A three-dimensional interpolation method is used to map the first flow field onto a fine mesh so that the first flow field serves as the initial flow field for the fine mesh.

[0106] This embodiment does not limit the method of using the first flow field as the initial flow field of the fine mesh. For example, a mesh point interpolation method can be used; or a linear interpolation method, i.e., three-dimensional interpolation, can be used; or a nearest-point method can be used. This embodiment can use a three-dimensional interpolation mapping method to map the first flow field to the system mesh to obtain a better initial flow field.

[0107] The wind farm wake assessment method provided in this invention involves acquiring measured topographic maps, satellite maps, projection systems, and turbine location information, constructing a physical model of the target terrain, and generating coarse and fine grids. Boundary conditions are obtained, and CFD is performed on the coarse grid to obtain a first calculation result and a first flow field, which is then used as the initial flow field for the fine grid. Actuator disk parameters are determined based on the first calculation result, and topology mapping and source term configuration are performed on the fine grid. Based on the boundary conditions, initial flow field, and actuator disk parameters, actuator disk-based CFD is performed on the fine grid to obtain a second calculation result and a second flow field. The wind direction of the second turbine location is extracted from the second calculation result based on the turbine location information. The actuator disk orientation is continuously corrected based on the wind direction of the second turbine location. When the difference between the actuator disk orientation and the wind direction of the second turbine location is less than a preset threshold, the correction is stopped, and wake parameters for each turbine location and wake parameters within a preset downstream range are obtained. This invention corrects the orientation of the actuator disk based on the wind direction at the turbine location, which can more realistically take into account the impact of complex terrain and the wake of upstream wind turbines on downstream wind turbines, improve the accuracy and effectiveness of wind farm wake assessment, overcome the problems of poor adaptability to complex terrain and low computational efficiency in wind farm wake assessment, effectively realize the refined assessment of wind farm wake, and provide a basis for wind farm power generation assessment and turbulence assessment, which is conducive to improving the reliability of wind resource assessment. Furthermore, a three-dimensional interpolation method is used to map the first flow field onto a fine mesh to obtain a better initial flow field and accelerate the convergence of the actuator disk-based CFD calculation. Also, when the wind direction at the second turbine location exceeds a preset threshold relative to the actuator disk orientation, the wind direction at the second turbine location is set as the current actuator disk orientation, and step S103 is executed to continuously correct the actuator disk orientation, ensuring the accuracy of the orientation correction. Furthermore, fewer parameters are used to calculate the actuator disk source terms, improving computational efficiency. Moreover, the actuator disk orientation, the configuration area of ​​the source terms, and the configuration parameters of the source terms are accurately determined based on the turbine size, model curve, and wind direction at the first turbine location. Furthermore, the inlet boundary conditions are obtained using the wind profile equation, improving the accuracy and efficiency of the calculated boundary conditions. Finally, a fine mesh is obtained using a coarse mesh, mesh tools, and turbine location information, achieving accurate mesh generation.

[0108] For a clearer understanding of this invention, please refer to the following details. Figure 2 , Figure 2The flowchart of a wind farm wake assessment method provided in this embodiment of the invention includes: processing a measured topographic map based on a projection system and satellite map to obtain a physical model; generating a coarse grid using a meshing tool; locally refining the mesh at the turbine locations using the meshing tool and turbine location information to generate a fine grid; performing CFD-based calculations on the coarse grid to obtain a first calculation result and a first flow field; mapping the first flow field using a three-dimensional interpolation method to obtain the initial flow field of the fine grid; post-processing the first calculation result and the first flow field obtained from the CFD calculation based on the turbine location information to obtain the wind direction at the first turbine location as the initial actuator disk orientation; and determining the region for source term configuration using topology mapping in conjunction with the turbine size; and determining the configuration parameters of the source term by combining the turbine curve, turbine size, and wind direction at the first turbine location, thereby completing the source term configuration.

[0109] Based on the determined initial flow field, source term configuration parameters, and boundary conditions, a fine-grid CFD solution based on an actuator disk is used to obtain the wind direction at the second turbine location. It is then determined whether the wind direction at the second turbine location is consistent with the orientation of the actuator disk (whether the difference is within a preset threshold range). If so, the wake parameters of each turbine location and the preset downstream range are obtained, including wind speed, wind direction, and pressure. If not, the process continues to execute the region for determining the source term configuration using topology mapping in conjunction with the turbine size and subsequent steps until the wind direction at the second turbine location is consistent with the orientation of the actuator disk (the difference is within a preset threshold range), and the wake parameters of each turbine location and the preset downstream range are obtained.

[0110] Experiments were conducted using this embodiment and the conventional method, and the results are compared below:

[0111] Taking the topography of an actual wind farm in South my country as the research object, the wake calculation results were obtained by using the traditional method and the modified method of the present invention, respectively. The traditional method directly uses the sector direction as the orientation of the actuator disk at each turbine location, and uses the same set of fine grids as in this embodiment.

[0112] Taking the 0° sector as an example, Table 1 shows the calculation results obtained at the machine point P2 using the traditional method and the method of this embodiment; Figure 3 This is a velocity distribution diagram of the wake region at location P2 provided in an embodiment of the present invention. (From Table 1 and...) Figure 3It can be seen that due to the influence of complex terrain, there is a significant difference between the wind direction and the sector direction at the turbine site, with a wind direction difference of up to 15°. In this case, traditional methods would lead to a deviation between the actuator disk orientation and the actual wind direction, while the correction method can effectively align the actuator disk orientation with the actual wind direction (normally, the actuator disk orientation is the windward direction). The actuator disk orientation causes differences in wind speed and direction at the turbine site (wind speed difference can reach 6%, wind direction difference is about 1.35°), which directly affects the power generation calculation at this turbine site, potentially causing an error of about 10% in power generation assessment. Simultaneously, the velocity distribution downstream of the turbine site is also affected; the deflection of the actuator disk orientation causes the force direction to tilt, thus underestimating wake losses. It is worth noting that the correction method in this embodiment only requires two corrections to make the actuator disk orientation basically consistent with the actual wind direction, and the computational efficiency is improved by about 9%. Experiments show that the flow field mapping method can effectively accelerate computational convergence. Therefore, the correction method of this invention demonstrates significant advantages in both the accuracy of the results and computational efficiency.

[0113] Table 1

[0114]

[0115]

[0116] The following describes the wind farm wake assessment device provided in the embodiments of the present invention. The wind farm wake assessment device described below and the wind farm wake assessment method described above can be referred to each other.

[0117] Please refer to the details. Figure 4 , Figure 4 A schematic diagram of a wind farm wake assessment device provided in an embodiment of the present invention may include:

[0118] The acquisition and construction module 100 is used to acquire measured topographic maps, satellite maps, projection systems and camera position information, construct a physical model of the target terrain, and generate coarse and fine grids;

[0119] The acquisition and solution module 200 is used to acquire boundary conditions and perform CFD solution on the coarse mesh to obtain the first calculation result and the first flow field, and use the first flow field as the initial flow field of the fine mesh;

[0120] The determination module 300 is used to determine the actuator disk parameters based on the first calculation result, and to perform topology mapping and source term configuration on the fine mesh.

[0121] The solver module 400 is used to perform CFD solution based on the actuation disk on the fine mesh according to the boundary conditions, initial flow field and actuation disk parameters to obtain the second calculation result and the second flow field;

[0122] The extraction module 500 is used to extract the wind direction of the second machine location from the second calculation result based on the machine location information.

[0123] The correction module 600 is used to continuously correct the orientation of the actuator disk according to the wind direction of the second machine point. When the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within the preset downstream range are obtained.

[0124] Furthermore, based on the above embodiments, the acquisition and solution module 200 may include:

[0125] The first calculation unit is used to calculate the boundary conditions based on roughness, thermal stability and reference wind speed;

[0126] The first obtaining unit is used to obtain the roughness based on the roughness map data;

[0127] The first solution element, the inlet boundary conditions used in the boundary conditions, are obtained by solving the wind profile equation based on the Monin-Obukhov similarity theory.

[0128] Furthermore, based on the above embodiments, the determining module 300 may include:

[0129] The first determining unit is used to extract the wind direction of the first machine location from the first calculation result based on the machine location information, and to determine the orientation of the actuation disk based on the wind direction of the first machine location.

[0130] The second determining unit is used to perform topological mapping on the fine mesh according to the wind turbine size to determine the configuration area of ​​the source item;

[0131] The third determining unit is used to determine the configuration parameters of the source term based on the model curve, the fan size, and the wind direction at the first fan location.

[0132] Furthermore, based on the above embodiments, the solving module 400 may include:

[0133] The solver module is used to perform CFD solving based on the actuated disk on the fine mesh to obtain the second calculation result;

[0134] The actuation disk source term calculation module is used to calculate the actuation disk source term based on the actuation disk CFD solution of the fine mesh according to the actuation disk parameters, so as to obtain the second flow field;

[0135] The formula for calculating the actuator disk source term is as follows:

[0136]

[0137]

[0138] Among them, Su The actuation disk source term is represented by ρ, air density by Δx, actuation disk thickness by u1, and inflow velocity by C. x Let represent the induction coefficient, and 'a' represent the induction factor. The induction factor 'a' is expressed as:

[0139]

[0140] Where C P C represents the power factor. T This represents the thrust coefficient.

[0141] Furthermore, based on the above embodiments, the correction module 600 may include:

[0142] The first execution unit is used to perform the steps of topology mapping and source term configuration on the fine grid if the difference between the orientation of the actuator disk and the wind direction of the second machine point is greater than a preset threshold.

[0143] The second execution unit is used to execute the steps of stopping correction and obtaining the wake parameters of each machine point and the wake parameters within a preset downstream range if the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than a preset threshold.

[0144] Furthermore, based on the above embodiments, the acquisition and construction module 100 may include:

[0145] The physical model building unit is used to process the measured topographic map based on the projection system and satellite map to obtain a physical model;

[0146] Coarse mesh definition elements are used to obtain coarse meshes using physical models and meshing tools;

[0147] The fine mesh determination element is used to refine the entire field mesh based on the coarse mesh, and to locally refine the mesh at the machine point based on the machine point information and meshing tools to obtain the fine mesh.

[0148] Furthermore, based on any of the above embodiments, the acquisition and solution module 200 may include:

[0149] The mapping element is used to map the first flow field to a fine mesh using a three-dimensional interpolation method, so that the first flow field serves as the initial flow field for the fine mesh.

[0150] It should be noted that the order of the modules and units in the aforementioned wind farm wake assessment device can be changed without affecting the logic.

[0151] This invention provides a wind farm wake assessment device. The device comprises an acquisition and construction module 100, which acquires measured topographic maps, satellite maps, projection systems, and machine location information, constructs a physical model of the target terrain, and generates coarse and fine meshes. An acquisition and solution module 200 acquires boundary conditions and performs CFD solving on the coarse mesh to obtain a first calculation result and a first flow field, using the first flow field as the initial flow field for the fine mesh. A determination module 300 determines the actuator disk parameters based on the first calculation result and performs topological mapping and source mapping on the fine mesh. The system comprises the following components: a solver module 400, which performs CFD-based solution on a fine mesh to obtain a second calculation result and a second flow field based on boundary conditions, initial flow field, and actuator disk parameters; an extraction module 500, which extracts the wind direction of the second turbine location from the second calculation result based on turbine location information; and a correction module 600, which continuously corrects the actuator disk orientation based on the wind direction of the second turbine location. When the difference between the actuator disk orientation and the wind direction of the second turbine location is less than a preset threshold, the correction stops, and the wake parameters of each turbine location and the wake parameters within a preset downstream range are obtained. This device corrects the actuator disk orientation based on the wind direction of the turbine location, which more realistically considers the impact of complex terrain and upstream turbine wake on downstream turbines, improving the accuracy and effectiveness of wind farm wake assessment. It overcomes the problems of poor adaptability to complex terrain and low computational efficiency in wind farm wake assessment, effectively realizing refined assessment of wind farm wake and providing a basis for wind farm power generation assessment and turbulence assessment, thus improving the reliability of wind resource assessment. Furthermore, a three-dimensional interpolation method is used to map the first flow field onto a fine mesh to obtain a better initial flow field and accelerate the convergence of the actuator disk-based CFD calculation. Also, when the wind direction at the second turbine location exceeds a preset threshold relative to the actuator disk orientation, the wind direction at the second turbine location is set as the current actuator disk orientation, and step S103 is executed to continuously correct the actuator disk orientation, ensuring the accuracy of the orientation correction. Furthermore, fewer parameters are used to calculate the actuator disk source terms, improving computational efficiency. Moreover, the actuator disk orientation, the configuration area of ​​the source terms, and the configuration parameters of the source terms are accurately determined based on the turbine size, model curve, and wind direction at the first turbine location. Furthermore, the inlet boundary conditions are obtained using the wind profile equation, improving the accuracy and efficiency of the calculated boundary conditions. Finally, a fine mesh is obtained using a coarse mesh, mesh tools, and turbine location information, achieving accurate mesh generation.

[0152] The following describes the wind farm wake assessment device provided in the embodiments of the present invention. The wind farm wake assessment device described below and the wind farm wake assessment method described above can be referred to each other.

[0153] Please refer to Figure 5 , Figure 5 A schematic diagram of a wind farm wake assessment device provided in an embodiment of the present invention may include:

[0154] Memory 10 is used to store computer programs;

[0155] Processor 20 is used to execute computer programs to implement the wind farm wake assessment method described above.

[0156] The system includes a memory 10, a processor 20, a communication interface 31, and a communication bus 32. The memory 10, processor 20, and communication interface 31 communicate with each other through the communication bus 32.

[0157] In this embodiment of the invention, the memory 10 is used to store one or more programs. The programs may include program code, which includes computer operation instructions. In this embodiment, the memory 10 may store programs for implementing the following functions:

[0158] Acquire measured topographic maps, satellite maps, projection systems, and camera location information, construct a physical model of the target terrain, and generate coarse and fine grids;

[0159] Obtain the boundary conditions and perform CFD solution on the coarse mesh to obtain the first calculation result and the first flow field. Use the first flow field as the initial flow field of the fine mesh.

[0160] Based on the first calculation result, determine the actuator disk parameters, and perform topology mapping and source term configuration on the fine mesh;

[0161] Based on the boundary conditions, initial flow field, and actuator disk parameters, a second calculation result and a second flow field are obtained by performing CFD solution based on actuator disk on the fine mesh.

[0162] The wind direction of the second machine location is extracted from the second calculation result based on the machine location information;

[0163] The orientation of the actuator disk is continuously corrected based on the wind direction at the second machine point. When the difference between the orientation of the actuator disk and the wind direction at the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within a preset downstream range are obtained.

[0164] In one possible implementation, the memory 10 may include a program storage area and a data storage area, wherein the program storage area may store the operating system and applications required for at least one function; and the data storage area may store data created during use.

[0165] Furthermore, memory 10 may include read-only memory and random access memory, providing instructions and data to the processor. A portion of the memory may also include NVRAM. The memory stores operating systems and operating instructions, executable modules, or data structures, or subsets thereof, or extended sets thereof, wherein the operating instructions may include various operating instructions for implementing various operations. The operating system may include various system programs for implementing various basic tasks and handling hardware-based tasks.

[0166] Processor 20 can be a central processing unit (CPU), an application-specific integrated circuit, a digital signal processor, a field-programmable gate array, or other programmable logic device. Processor 20 can be a microprocessor or any conventional processor. Processor 20 can call programs stored in memory 10.

[0167] Communication interface 31 can be an interface for the communication module, used to connect with other devices or systems.

[0168] Of course, it should be noted that, Figure 5 The structure shown does not constitute a limitation on the wind farm wake assessment device in the embodiments of this application. In practical applications, the wind farm wake assessment device may include more than Figure 5 More or fewer components as shown, or combinations of certain components.

[0169] The storage medium provided in the embodiments of the present invention is described below. The storage medium described below can be referred to in correspondence with the wind farm wake assessment method described above.

[0170] The present invention also provides a storage medium storing a computer program, which, when executed by a processor, implements the steps of the wind farm wake assessment method described above.

[0171] The storage medium can include various media that can store program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0172] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to in the method section.

[0173] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0174] Finally, it should be noted that in this document, relationships such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0175] The above provides a detailed description of the wind farm wake assessment method, apparatus, equipment, and storage medium provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for assessing the wake of a wind farm, characterized in that, include: Acquire measured topographic maps, satellite maps, projection systems, and camera location information, construct a physical model of the target terrain, and generate coarse and fine grids; Obtain boundary conditions, and perform CFD solution on the coarse mesh to obtain the first calculation result and the first flow field, and use the first flow field as the initial flow field of the fine mesh; Based on the first calculation result, the actuator disk parameters are determined, and the fine mesh is topologically mapped and the source terms are configured. Based on the boundary conditions, the initial flow field, and the actuation disk parameters, the fine mesh is solved using CFD based on the actuation disk to obtain the second calculation result and the second flow field. Based on the location information, the wind direction of the second location is extracted from the second calculation result; The orientation of the actuator disk is continuously corrected according to the wind direction at the second machine point. When the difference between the orientation of the actuator disk and the wind direction at the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within a preset downstream range are obtained. The step of performing CFD solution based on the actuated disk on the fine mesh to obtain the second calculation result and the second flow field includes: The second calculation result is obtained by performing CFD solution based on the actuated disk on the fine mesh; The actuator disk source terms are calculated by CFD solution based on the actuator disk parameters to obtain the second flow field. The formula for calculating the actuation disk source term is as follows: ； ； in, Indicates the source item of the actuator disk. Indicates air density, Indicates the thickness of the actuator disk. Indicates the incoming wind speed. Represents the induction coefficient. Represents the inducing factor; wherein the inducing factor Represented as: ； in Indicates the power factor. Indicates the thrust coefficient; The process of continuously correcting the orientation of the actuator disk based on the wind direction at the second machine location, and stopping the correction when the difference between the orientation of the actuator disk and the wind direction at the second machine location is less than a preset threshold, and obtaining the wake parameters of each machine location and the wake parameters within a preset downstream range, includes: If the difference between the orientation of the actuator disk and the wind direction of the second machine point is greater than the preset threshold, then the wind direction of the second machine point is taken as the orientation of the actuator disk, and the steps of performing topology mapping and source term configuration on the fine mesh are executed. If the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than the preset threshold, the step of stopping correction and obtaining the wake parameters of each machine point and the wake parameters within the preset downstream range is executed.

2. The wind farm wake assessment method according to claim 1, characterized in that, The acquisition of boundary conditions includes: The boundary conditions are calculated based on roughness, thermal stability, and reference wind speed. The roughness is obtained based on the roughness map data; The inlet boundary conditions are obtained by solving the wind profile equation based on the Monin-Obukhov similarity theory.

3. The wind farm wake assessment method according to claim 1, characterized in that, The step of determining the actuator disk parameters based on the first calculation result and performing topology mapping and source term configuration on the fine mesh includes: Based on the machine location information, the wind direction of the first machine location is extracted from the first calculation result, and the orientation of the actuation disk is determined based on the wind direction of the first machine location; The configuration area of ​​the source item is determined by performing topological mapping on the fine mesh according to the wind turbine size; The configuration parameters of the source term are determined based on the model curve, the fan size, and the wind direction at the first fan location.

4. The wind farm wake assessment method according to claim 1, characterized in that, The process of acquiring measured topographic maps, satellite maps, projection systems, and camera position information, constructing a physical model of the target terrain, and generating coarse and fine grids includes: The measured topographic map is processed according to the projection system and the satellite map to obtain the physical model; The coarse mesh is obtained using the physical model and meshing tools described above; The entire field mesh is refined based on the coarse mesh, and the fine mesh is obtained by locally refining the mesh at the machine points according to the machine point information and the mesh tool.

5. The wind farm wake assessment method according to any one of claims 1 to 4, characterized in that, The first flow field, serving as the initial flow field for the fine mesh, includes: The first flow field is mapped onto the fine mesh using a three-dimensional interpolation method, so that the first flow field serves as the initial flow field of the fine mesh.

6. A wind farm wake assessment device, characterized in that, include: The acquisition and construction module is used to acquire measured topographic maps, satellite maps, projection systems and camera position information, construct a physical model of the target terrain, and generate coarse and fine grids; The acquisition and solution module is used to acquire boundary conditions and perform CFD solution on the coarse mesh to obtain a first calculation result and a first flow field, and use the first flow field as the initial flow field of the fine mesh; The determination module is used to determine the actuator disk parameters based on the first calculation result, and to perform topology mapping and source term configuration on the fine mesh; The solver module is used to perform CFD solution based on the actuation disk on the fine mesh according to the boundary conditions, the initial flow field and the actuation disk parameters to obtain the second calculation result and the second flow field; The extraction module is used to extract the wind direction of the second machine location from the second calculation result based on the machine location information; The correction module is used to continuously correct the orientation of the actuator disk according to the wind direction of the second machine point. When the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than a preset threshold, the correction is stopped and the wake parameters of each machine point and the wake parameters within a preset downstream range are obtained. The solution module includes: The solving element is used to perform CFD solving based on the actuated disk on the fine mesh to obtain the second calculation result; The actuator disk source term calculation unit is used to perform CFD calculation based on the actuator disk on the fine mesh according to the actuator disk parameters to calculate the actuator disk source term, so as to obtain the second flow field; The formula for calculating the actuator disk source term is as follows: ； ； in, Indicates the source item of the actuator disk. Indicates air density, Indicates the thickness of the actuator disk. Indicates the incoming wind speed. Represents the induction coefficient. Represents the inducing factor; wherein the inducing factor Represented as: ； in Indicates the power factor. Indicates the thrust coefficient; The correction module includes: The first execution unit is used to perform the steps of topology mapping and source term configuration on the fine grid if the difference between the orientation of the actuator disk and the wind direction of the second machine point is greater than a preset threshold. The second execution unit is used to execute the steps of stopping correction and obtaining the wake parameters of each machine point and the wake parameters within a preset downstream range if the difference between the orientation of the actuator disk and the wind direction of the second machine point is less than a preset threshold.

7. A wind farm wake assessment device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the steps of the wind farm wake assessment method as described in any one of claims 1 to 5.

8. A storage medium, characterized in that, The storage medium stores a computer program, which, when executed by a processor, implements the steps of the wind farm wake assessment method as described in any one of claims 1 to 5.

Images