CFD-based wind resource evaluation and generating capacity prediction method and system

Abstract

The invention relates to the technical field of wind field simulation, in particular to a CFD-based wind resource evaluation and generating capacity prediction method and system. The method comprises the steps of obtaining wind regime data of a wind measurement point and all result points in a group of wind direction sectors through directional calculation, obtaining wind measurement data of the wind measurement point in any wind direction, calculating weights of two wind direction sectors adjacent to the wind direction of the wind measurement data in a gradual working condition, and establishing a wind regime mapping relation between the wind measurement point and the result points, and calculating to obtain equivalent wind data of the result points when the wind measurement points are in any wind direction, analyzing the equivalent wind data of all the result points, evaluating wind resources of the target wind field, and predicting the generating capacity of the target wind field. According to the method, more accurate, efficient and reliable data support can be provided for wind resource assessment and power generation prediction through the equivalent wind data at all result points.

Description

Technical Field

[0001] The present invention relates to the technical field of wind farm simulation, and in particular to a CFD-based wind resource assessment and power generation prediction method and system. Background Art

[0002] As global warming continues, efforts to develop new energy sources, including wind power, have intensified. Therefore, how to capture wind energy economically and efficiently has become a long-term research topic in the wind energy industry, especially how to accurately evaluate wind resources and calculate power generation.

[0003] At present, the core content of wind farm wind resource assessment and power generation calculation is divided into two steps, including directional calculation and comprehensive calculation. Directional calculation is to establish the wind parameter relationship between the wind measurement point and the wind farm (wind turbine); comprehensive calculation is to use the directional calculation results, combined with the wake model, wind farm layout, wind measurement data, and wind turbine characteristic curves, to perform wind resource assessment and power generation calculation without or with wake consideration. However, in the above-mentioned wind farm resource assessment and power generation calculation, the wind measurement data obtained by equipment such as wind towers is difficult to cover the entire wind farm due to the cost and limited data points, and cannot provide more accurate, efficient and reliable data support for wind resource assessment and power generation prediction. Summary of the invention

[0004] The technical problem to be solved by the embodiments of the present invention is to provide a CFD-based wind resource assessment and power generation prediction method and system to solve the problem in the prior art that it is difficult to obtain wind measurement data to cover the entire wind farm, resulting in the inability to provide more accurate, efficient and reliable data support for wind resource assessment and power generation prediction.

[0005] The present invention discloses a wind resource assessment and power generation prediction method based on CFD, comprising: Establishing a CFD model of the target wind field, and obtaining wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation of the CFD model; Obtaining wind measurement data of the wind measurement point at any wind direction, and determining two wind direction sectors adjacent to the wind direction of the wind measurement data; Calculate the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data corresponding to the wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data; Establish a wind condition mapping relationship between the wind measurement point and the result point, combine the acquired wind measurement data, the calculated weight, and the wind condition data corresponding to the wind measurement point and the result point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, and calculate and obtain the equivalent wind data of the result point at any wind direction of the wind measurement point; Analyze the equivalent wind data of all the result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

[0006] Optionally, obtaining wind condition data of a wind measurement point and all result points within a set of wind direction sectors by directional calculation includes: A group of wind direction sectors having the same wind direction angle division range are established for the wind measurement point and each of the result points; Obtaining the wind resources and power generation of the CFD model, and obtaining the directional calculation results of the wind measurement point and each of the result points in a set of wind direction sectors through directional calculation, wherein the directional calculation results include wind speed components and turbulent kinetic energy; According to the directional calculation results of the wind measuring point and each of the result points, the wind condition data of the wind measuring point and each of the result points in each of the wind direction sectors are obtained by calculation, and the wind condition data include wind acceleration factor, inflow angle, wind direction deflection angle and calculated turbulence intensity.

[0007] Optionally, the calculating the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition includes: The wind direction angle is extracted from the wind measurement data. According to the wind direction angle of the wind measurement data and the wind direction deflection angle corresponding to the wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data are obtained by interpolation calculation. The weight calculation formula is:

[0008]

[0009] In the formula, Represents the weight of the first wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Represents the weight of the second wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates wind direction of wind measurement data wind direction angle.

[0010] Optionally, the equivalent wind data obtained by calculation includes the wind direction angle of the result point, and the calculation method of the wind direction angle of the result point includes: According to the wind direction angles of the two wind direction sectors adjacent to the wind direction of the wind measurement data and the wind direction deflection angle corresponding to the result point, an equivalent wind direction angle calculation formula for the result point under any wind direction is established. The equivalent wind direction angle calculation formula is:

[0011] In the formula, Indicates any wind direction at the wind measurement point i The equivalent wind direction angle of the result point is Indicates the wind direction deflection angle of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data.

[0012] Optionally, the equivalent wind data obtained by calculation include wind speed and turbulence intensity at the result point; The calculation method of the wind speed at the result point includes: The wind speed is extracted from the wind measurement data, and when any wind direction of the wind measurement point is established, the mapping relationship between the wind speed of the wind measurement data of the wind measurement point and the directional wind acceleration factor is expressed as follows:

[0013] When establishing the mapping relationship between the wind speed at the result point and the directional wind acceleration factor at any wind direction of the wind measurement point, the relationship expression is:

[0014] By combining the above two relational expressions, a calculation formula for the equivalent wind speed of the result point at any wind direction of the wind measurement point is established. The equivalent wind speed calculation formula is:

[0015] In the formula, Indicates the equivalent wind speed of the result point at any wind direction of the wind measurement point. Indicates the wind speed at any wind direction at the wind measuring point. Indicates the wind acceleration factor of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data; Indicates the wind acceleration factor of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Positive correlation; The calculation method of the turbulence intensity at the result point includes: Obtaining the measured environment turbulence intensity corresponding to the wind direction of the wind measurement data, and the interpolated calculated turbulence intensity of the wind measurement point in two adjacent wind direction sectors when the wind measurement point has any wind direction, and calculating the intensity difference between the measured environment turbulence intensity and the interpolated calculated turbulence intensity corresponding to the wind measurement point; When obtaining any wind direction at the wind measuring point, the turbulence intensity is calculated by interpolating the result point in two adjacent wind direction sectors; According to the wind acceleration factors of the wind measuring point and the result point in each wind direction sector, a ratio relationship between the wind measuring point and the result point corresponding to the wind acceleration factors when the wind measuring point has any wind direction is established, and the ratio relationship expression is:

[0016] According to the intensity difference, the interpolated calculated turbulence intensity corresponding to the result point, and the ratio relationship, an equivalent turbulence intensity calculation formula for the result point under any wind direction is established. The equivalent turbulence intensity calculation formula is:

[0017] In the formula, It represents the equivalent turbulence intensity of the result point at any wind direction of the wind measurement point. Indicates the interpolation calculation turbulence intensity corresponding to the result point, represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity corresponding to the wind measurement point, It indicates the wind acceleration factor at the wind measuring point at any wind direction. Indicates the wind acceleration factor of the result point at any wind direction of the wind measurement point. It indicates the ratio of the wind acceleration factors corresponding to the wind measurement point and the result point.

[0018] Optionally, the calculating the intensity difference between the measured environment turbulence intensity and the interpolated calculated turbulence intensity corresponding to the wind measurement point includes: The wind speed and the wind speed standard deviation are extracted from the wind measurement data, and the measured environmental turbulence intensity corresponding to the wind direction of the wind measurement data is calculated based on the wind speed and the wind speed standard deviation of the wind measurement data. The calculation formula for the measured environmental turbulence intensity is:

[0019] In the formula, Indicates the measured environmental turbulence intensity of the wind direction corresponding to the wind measurement point. Indicates the standard deviation of wind speed at the wind measurement point corresponding to the wind direction of the wind measurement data, Indicates the wind speed at the wind measuring point corresponding to the wind direction of the wind measurement data; The wind speed component of the wind measuring point is extracted from the directional calculation result, and the turbulence intensity of the wind measuring point in two adjacent wind direction sectors is calculated when obtaining any wind direction of the wind measuring point. The calculation formula for the turbulence intensity of the wind measuring point is:

[0020] In the formula, Indicates that the wind measuring point is in the adjacent wind direction sector when the wind direction is arbitrary i The calculated turbulence intensity is Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The turbulent kinetic energy corresponding to the wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from west to east at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from south to north at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from bottom to top corresponding to the wind measuring point; According to the calculated turbulence intensity of the wind measuring point in two adjacent wind direction sectors, the interpolated calculated turbulence intensity of the wind measuring point in two adjacent wind direction sectors is obtained by interpolation calculation. The calculation formula of the interpolated calculated turbulence intensity is:

[0021] In the formula, It indicates the interpolation calculation of turbulence intensity at the wind measuring point in two adjacent wind direction sectors. It represents the calculated turbulence intensity of the wind measuring point in the first adjacent wind direction sector when the wind measuring point has any wind direction. It represents the calculated turbulence intensity of the wind measuring point in the adjacent second wind direction sector when the wind direction is arbitrary at the wind measuring point; The intensity difference is calculated according to the measured environment turbulence intensity obtained by calculation and the interpolation calculated turbulence intensity. The calculation formula of the intensity difference is:

[0022] In the formula, Represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity.

[0023] Optionally, a method for obtaining different types of equivalent wind data according to the type of the wind measurement data is also included: Based on data collection, the time series wind measurement data at the wind measurement point is acquired, and the time series wind measurement data is calculated through equivalent wind data to obtain the time series equivalent wind data of the result point at the wind measurement point in any wind direction; Acquire TAB format wind measurement data based on file upload, calculate the TAB format wind measurement data through equivalent wind data, and obtain TAB format equivalent wind data of the result point at any wind direction of the wind measurement point; If the TAB format wind measurement data uses an intermediate value to represent the wind speed bin interval, the bin width of each wind speed bin interval in the TAB format wind measurement data is increased by a preset multiple, and then the processed TAB format wind measurement data is used to calculate the equivalent wind data; If the wind measurement data in the TAB format represents the wind speed bin interval with an upper limit value, the wind measurement data in the TAB format is used for equivalent wind data calculation; Compare and analyze the time series equivalent wind data and TAB format equivalent wind data of all the result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

[0024] Optionally, the method further includes calculating and obtaining the single-tower equivalent wind data of the result point according to different types of wind measurement data: When the wind measurement data is time series wind measurement data, the wind measurement data is calculated through equivalent wind data to obtain the single tower time series equivalent wind data of the result point when the wind measurement point is in any wind direction; When the wind measurement data is in TAB format, the wind speed range and wind direction sector are added to the wind measurement data to obtain TAB format subdivided wind measurement data, the TAB format subdivided wind measurement data is calculated through equivalent wind data to obtain the TAB format subdivided equivalent wind data of the result point, and then all the TAB format subdivided equivalent wind data of the result points are aggregated to obtain the single-tower TAB format equivalent wind data of the result point when the wind measurement point is in any wind direction.

[0025] Optionally, the method further includes obtaining multi-tower equivalent wind data of the result point according to the single-tower equivalent wind data of the result point: According to the single-tower equivalent wind data of the result point, a multi-tower weighted comprehensive calculation is performed to obtain the multi-tower equivalent wind data of the result point. The calculation formula of the multi-tower equivalent wind data of the result point is:

[0026]

[0027] In the formula, for Total number of wind towers, For the The confidence factor of the wind tower, The result point is based on The wind speed components and turbulence intensity calculated by the wind towers are Click to the result The distance between the wind towers, For the The weighted average coefficient of the wind towers is Multi-tower equivalent wind data for the result point.

[0028] The present invention also discloses an evaluation and prediction system, which adopts the above-mentioned CFD-based wind resource evaluation and power generation prediction method, and the system includes: Directional calculation module, used to establish the CFD model of the target wind field, and obtain the wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation; A data processing module is used to obtain wind measurement data of a wind measurement point at any wind direction and determine two wind direction sectors adjacent to the wind direction of the wind measurement data; A weight calculation module, used to calculate the weights of two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data of the two wind direction sectors adjacent to the wind direction of the wind measurement data; The equivalent wind data calculation module is used to establish the wind condition mapping relationship between the wind measurement point and the result point, and calculate the equivalent wind data of the result point at any wind direction of the wind measurement point by combining the acquired wind measurement data, the calculated weight, and the wind condition data of the two wind direction sectors adjacent to the wind measurement data. The evaluation and prediction module is used to analyze the equivalent wind data of all result points, and to conduct wind resource evaluation and power generation prediction for the target wind farm.

[0029] Compared with the prior art, the CFD-based wind resource assessment and power generation prediction method and system provided by the embodiment of the present invention have the following beneficial effects: Directional calculation of the CFD model of the target wind farm provides wind condition data in a set of wind direction sectors for the analysis objects including wind measurement points and all result points in the calculation domain. By calculating the weights of two wind direction sectors adjacent to the wind measurement data under each working condition, and using the calculation of equivalent wind data of the result point, the wind measurement data of the wind measurement point can be mapped to the result point, so that the result point has data in the same form as the wind measurement data, so as to obtain equivalent wind data at all result points. Therefore, based on the equivalent wind data at all result points, more accurate, efficient and reliable data support can be provided for wind resource assessment and power generation prediction. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments, in which: Figure 1 A schematic block diagram of the steps of a CFD-based wind resource assessment and power generation prediction method provided in an embodiment of the present invention; Figure 2 A relationship diagram between a result point under any wind direction and the directional results of the corresponding wind measurement points in the adjacent wind direction sectors provided by an embodiment of the present invention; Figure 3 A schematic block diagram of a process for calculating equivalent wind data of a single wind tower and multiple wind towers according to time series data and TAB format data provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0031] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. Now, in conjunction with the accompanying drawings, the preferred embodiments of the present invention are described in detail.

[0032] The present invention discloses a wind resource assessment and power generation prediction method based on CFD, such as Figure 1 As shown, including: S1. Establish a CFD model of the target wind field, and obtain wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation of the CFD model; S2. Obtain wind measurement data of the wind measurement point at any wind direction, and determine two wind direction sectors adjacent to the wind direction of the wind measurement data; S3, calculating the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data corresponding to the wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data; S4, establishing a wind condition mapping relationship between the wind measurement point and the result point, combining the acquired wind measurement data, the calculated weight, and the wind condition data corresponding to the wind measurement point and the result point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, and calculating and obtaining the equivalent wind data of the result point at any wind direction of the wind measurement point; S5. Analyze the equivalent wind data of all result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

[0033] Through the implementation of the above-mentioned wind resource assessment and power generation prediction method embodiment, computational fluid dynamics (CFD) is a science that studies fluid flow through numerical methods and computer simulations. It can simulate the discrete distribution of fluids in a continuous area by solving differential equations of fluid motion, thereby approximately simulating the actual fluid flow situation. Based on CFD wind field modeling, a CFD model of the target wind field can be established by obtaining some terrain parameters of the target wind field, and the CFD model of the target wind field can be used for directional calculation to provide wind condition data of the analysis object including wind measurement points and all result points in the calculation domain within a set of wind direction sectors, that is, the directional results of wind measurement points and the directional results of result points. By calculating the weights of two wind direction sectors adjacent to the wind direction of the wind data under each working condition, and using the calculation of equivalent wind data of the result point, the wind measurement data of the wind measurement point can be mapped to the result point, so that the result point has data in the same form as the wind measurement data, so as to integrate the equivalent wind data at all result points, cover the entire target wind farm, and estimate the average wind force received by the target wind farm based on the wind direction, wind speed, turbulence intensity and other data contained in the equivalent wind data of the result point, and the equivalent wind speed can be used to predict the power generation of the target wind farm. Then, by collecting and analyzing these data, more accurate, efficient and reliable data support is provided for wind resource assessment and power generation prediction.

[0034] Furthermore, the wind condition data of the wind measurement point and all result points within a set of wind direction sectors are obtained through directional calculation, including: A group of wind direction sectors with the same wind direction angle division range are established for the wind measurement point and each result point; Obtain the wind resources and power generation of the CFD model, and obtain the directional calculation results of the wind measurement point and each result point in a set of wind direction sectors through directional calculation. The directional calculation results include wind speed components and turbulent kinetic energy; According to the directional calculation results of the wind measuring point and each result point, the wind condition data of the wind measuring point and each result point in each wind direction sector are obtained by calculation. The wind condition data include wind acceleration factor, inflow angle, wind direction deflection angle and calculated turbulence intensity.

[0035] Through the implementation of the above-mentioned wind resource assessment and power generation prediction method embodiment, since the wind direction angle is usually measured clockwise from the north direction, ranging from 0° (due north) to 360°, the wind direction angle can be divided according to a certain interval, so that each wind direction sector in a group of wind direction sectors corresponds to a wind direction angle and is included in the range of 0°-360°, as shown in the wind direction sector division table: Wind direction sector division table

[0036] As shown in the wind direction sector division table, the embodiment of the present invention preferably divides the range of 0°-360° into 16 wind direction sectors, and for the convenience of interpolation calculation, on the basis of the 16 wind direction sectors, an additional directional calculation result of 360° in a specific direction is added, which is consistent with 0°. Among them, in the wind direction sector division table, each wind direction sector corresponds to a box for statistics of wind condition data such as wind acceleration factor S, inflow angle In, wind direction deflection angle Dev, and calculated turbulence intensity TI at each point, that is, the wind measurement point and each result point are statistically analyzed for wind condition data through the above-mentioned wind direction sector division table. Here, only a group of wind direction sector divisions are shown in a tabular form, and the statistics of specific data are not involved.

[0037] In addition, through directional calculation, the calculation results of the wind measurement point and each result point in a set of wind direction sectors are obtained respectively. The directional calculation results mainly include wind speed components (u, v, w), turbulent kinetic energy (k), atmospheric pressure (p), air density ( ) and other physical fields. Among them, the wind speed component usually refers to the projection of the wind speed in the direction of the wind direction angle, and the turbulent kinetic energy is a parameter that describes the turbulence intensity of the fluid. Based on the wind speed component and turbulent kinetic energy, the wind acceleration factor S, inflow angle In, wind direction deflection angle Dev, calculated turbulence intensity TI and other wind condition data can be calculated through the relevant definition formulas. Among them, the wind acceleration factor describes the degree of wind acceleration during the flow process, which is usually used to evaluate the potential of wind energy resources. The inflow angle refers to the angle between the wind direction and the plane of the wind turbine blade, which affects the capture efficiency of wind energy. The wind direction deflection angle describes the angle difference between the actual wind direction and the reference wind direction (usually the wind direction at the wind measurement point). The calculated turbulence intensity measures the intensity of turbulence and has an important impact on the stability and performance of the wind turbine. Therefore, the wind acceleration factor, inflow angle, wind direction deflection angle and calculated turbulence intensity of the wind measurement point and each result point in each wind direction sector are obtained by calculation, so as to provide accurate, efficient and reliable data for the subsequent equivalent wind data calculation.

[0038] As mentioned above, in step S2, while obtaining the wind measurement data, it is also necessary to prepare the analysis object, including obtaining the coordinates of the wind measurement point and each result point in a unified coordinate system. Since the wind acceleration factor of the wind measurement point and all the result points in a set of wind direction sectors is obtained by directional calculation based on the CFD model of the target wind field in step S1, S、 Inflow Angle In , wind direction deflection angle Dev , calculate turbulence intensity TI, Through the coordinates of the wind measurement points and the result points, the wind acceleration factor at any point can be quickly located and retrieved. S、 Inflow Angle In , wind direction deflection angle Dev , calculate turbulence intensity TIWind field information is provided to facilitate subsequent equivalent data calculations.

[0039] Furthermore, combined with Figure 2 As shown in the figure, the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition are calculated, including: The wind direction angle is extracted from the wind measurement data. According to the wind direction angle of the wind measurement data and the wind direction deflection angle of the corresponding wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data are obtained by interpolation calculation. The weight calculation formula is:

[0040]

[0041] In the formula, Represents the weight of the first wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Represents the weight of the second wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction angle of wind direction i in wind measurement data.

[0042] Through the implementation of the above-mentioned wind resource assessment and power generation prediction method embodiment, combined with the above-mentioned wind direction sector division table, for example, the wind direction angle of the wind direction of the wind measurement data If 50° is used, the weight of the 45° wind direction sector is calculated. and the weight of the 67.5° wind direction sector .

[0043] Furthermore, combined with Figure 2 As shown, the equivalent wind data obtained by calculation includes the wind direction angle of the result point, and the calculation method of the wind direction angle of the result point includes: According to the wind direction angles of the two wind direction sectors adjacent to the wind direction of the wind measurement data and the wind direction deflection angle of the corresponding result point, the equivalent wind direction angle calculation formula of the result point under any wind direction is established. The equivalent wind direction angle calculation formula is:

[0044] In the formula, Indicates any wind direction at the wind measurement point i The equivalent wind direction angle of the result point is Indicates the wind direction deflection angle of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data.

[0045] Furthermore, the equivalent wind data obtained by calculation include wind speed and turbulence intensity at the result point; The calculation methods of wind speed at the result point include: The wind speed is extracted from the wind measurement data, and the mapping relationship between the wind speed of the wind measurement data at the wind measurement point and the directional wind acceleration factor is established when any wind direction is set. The relationship expression is:

[0046] When establishing an arbitrary wind direction at a wind measurement point, the mapping relationship between the wind speed at the result point and the directional wind acceleration factor is expressed as follows:

[0047] Combining the above two relational expressions, the calculation formula for the equivalent wind speed of the result point at any wind direction at the wind measurement point is established. The calculation formula for the equivalent wind speed is:

[0048] In the formula, Indicates the equivalent wind speed of the result point at any wind direction of the wind measurement point. Indicates the wind speed at any wind direction at the wind measuring point. Indicates the wind acceleration factor of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data; Indicates the wind acceleration factor of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Positive correlation; The calculation methods of turbulence intensity at result points include: Obtain the measured environmental turbulence intensity corresponding to the wind direction of the wind measurement data, and the interpolated turbulence intensity of the wind measurement point in two adjacent wind direction sectors when the wind measurement point has any wind direction, and calculate the intensity difference between the measured environmental turbulence intensity and the interpolated turbulence intensity corresponding to the wind measurement point; When obtaining any wind direction at a wind measuring point, the turbulence intensity is calculated by interpolating the result point in two adjacent wind direction sectors; According to the wind acceleration factors of the wind measuring point and the result point in each wind direction sector, the ratio relationship between the wind measuring point and the corresponding wind acceleration factors of the result point is established when the wind measuring point has any wind direction. The ratio relationship expression is:

[0049] According to the intensity difference, the interpolation calculation of the turbulence intensity and the ratio relationship corresponding to the result point, the equivalent turbulence intensity calculation formula of the result point under any wind direction is established. The equivalent turbulence intensity calculation formula is:

[0050] In the formula, It represents the equivalent turbulence intensity of the result point at any wind direction of the wind measurement point. Indicates the interpolation calculation turbulence intensity corresponding to the result point, represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity corresponding to the wind measurement point, It indicates the wind acceleration factor at the wind measuring point at any wind direction. Indicates the wind acceleration factor of the result point at any wind direction of the wind measurement point. It indicates the ratio of the wind acceleration factors corresponding to the wind measurement point and the result point.

[0051] Further, the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity corresponding to the wind measurement point is calculated, including: The wind speed and wind speed standard deviation are extracted from the wind measurement data. The measured environmental turbulence intensity corresponding to the wind direction of the wind measurement data is calculated based on the wind speed and standard deviation of the wind measurement data. The calculation formula for the measured environmental turbulence intensity is:

[0052] In the formula, It indicates the measured environmental turbulence intensity at the wind measurement point corresponding to the wind direction of the wind measurement data. Indicates the standard deviation of wind speed at the wind measurement point corresponding to the wind direction of the wind measurement data. Indicates the wind speed at the wind measurement point corresponding to the wind direction of the wind measurement data; The wind speed component of the wind measuring point is extracted from the directional calculation results. When the wind direction of the wind measuring point is obtained, the turbulence intensity of the wind measuring point in the two adjacent wind direction sectors is calculated. The calculation formula for the turbulence intensity of the wind measuring point is:

[0053] In the formula, Indicates that the wind measuring point is in the adjacent wind direction sector when the wind direction is arbitrary The calculated turbulence intensity is Indicates the wind direction sector adjacent to the wind direction of the wind measurement data The turbulent kinetic energy corresponding to the wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data The wind speed component from west to east at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data The wind speed component from south to north at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data The wind speed component from bottom to top corresponding to the wind measuring point; According to the calculated turbulence intensity of the wind measuring point in two adjacent wind direction sectors, the interpolation calculation of the turbulence intensity of the wind measuring point in two adjacent wind direction sectors is obtained by interpolation calculation. The calculation formula of the interpolation calculation of turbulence intensity is:

[0054] In the formula, It indicates the interpolation calculation of turbulence intensity at the wind measuring point in two adjacent wind direction sectors. It represents the calculated turbulence intensity of the wind measuring point in the first adjacent wind direction sector when the wind measuring point has any wind direction. It represents the calculated turbulence intensity of the wind measuring point in the adjacent second wind direction sector when the wind direction is arbitrary at the wind measuring point; According to the calculated measured environmental turbulence intensity and the interpolated turbulence intensity, the intensity difference is calculated. The calculation formula of the intensity difference is:

[0055] In the formula, Represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity.

[0056] Combining the above calculation formulas, we can further obtain the equivalent turbulence intensity calculation formula:

[0057]

[0058]

[0059] Furthermore, combined with Figure 3 As shown, it also includes a method for obtaining different types of equivalent wind data according to the type of wind measurement data: Based on data collection, the time series wind measurement data at the wind measurement point is obtained, and the time series wind measurement data is calculated through equivalent wind data to obtain the time series equivalent wind data of the result point at any wind direction of the wind measurement point; Obtain TAB format wind measurement data based on file upload, calculate the TAB format wind measurement data through equivalent wind data, and obtain TAB format equivalent wind data of the result point at any wind direction; If the TAB format wind measurement data uses the middle value to represent the wind speed bin interval, the bin width of each wind speed bin interval in the TAB format wind measurement data is increased by a preset multiple, and then the equivalent wind data is calculated for the processed TAB format wind measurement data; If the wind measurement data in TAB format represents the wind speed bin interval with the upper limit value, the equivalent wind data calculation is performed on the wind measurement data in TAB format; Compare and analyze the time series equivalent wind data and TAB format equivalent wind data of all result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

[0060] Through the implementation of the above-mentioned wind resource assessment and power generation prediction method embodiment, the wind measurement data includes wind measurement data of a wind tower, mesoscale data, multi-mesoscale data, etc., and its types are divided into time series wind measurement data and TAB format wind measurement data.

[0061] Time series wind measurement data are represented by wind towers with timestamps, radar wind measurement data and mesoscale data. The wind towers involved in the comprehensive calculation need to use wind measurement data at the highest altitude. Its data variables are shown in the time series wind measurement data variable table, including time date&time, wind speed v, wind direction dir, turbulence intensity ti, standard deviation std, temperature tem, relative humidity rh, and air pressure p.

[0062] Time series wind data variable table

[0063] There are two types of TAB format wind measurement data: one is from time series data conversion, and the other is from file upload. There are also two types of TAB format wind measurement data based on file upload: (1) If the TAB format data comes from a type such as Windographer that uses the middle value to represent the wind speed bin interval, the wind speed used in average wind speed, energy density, annual power generation, etc. needs to be additionally processed by adding 0.5*bin. This is because the wind speed reported by software such as Windographer is an estimate of the middle value of each bin, while the actual average wind speed or energy density and other parameters may require more accurate bin boundary values; (2) If the TAB format data comes from a type of data such as MeteodynWT that uses the upper limit value to represent the wind speed bin interval, no additional processing is required for average wind speed, energy density, annual power generation, etc. This means that the center value of each wind speed bin (i.e. the middle value of the bin) can be used directly to calculate parameters such as average wind speed, energy density, and annual power generation. Because the wind speed reported by software such as MeteodynWT is the upper limit value of each bin, using the bin middle value directly is accurate enough.

[0064] As mentioned above, Windographer is a wind resource and visualization tool designed specifically for meteorological professionals. It is mainly used to import, analyze and visualize wind resource data measured by meteorological towers, sodas or lidars. Meteodyn WT is a software mainly used for wind resource assessment in complex terrain. Its core is a set of computational fluid dynamics (CFD) codes specifically used to study atmospheric airflow. It can significantly reduce the uncertainty in wind farm design caused by complex terrain (such as mountains, forests, mudflats, etc.) and improve the accuracy of wind energy resource assessment. ‌

[0065] For the TAB format wind measurement data uploaded based on files, for the wind measurement data files with timestamp, wind speed, wind direction, and wind speed standard deviation, without considering the data truncation error, it is approximately equivalent to 1 frequency distribution table + 1 turbulence intensity table + 1 turbulence intensity relative standard deviation table. The frequency distribution table shows the frequency distribution of wind speed and wind direction, which usually includes each bin of wind speed and the corresponding wind speed value, as well as the number of wind speed occurrences (frequency) in each wind direction sector. In the TAB format wind measurement data, wind speed and wind direction are recorded by timestamp, and the frequency distribution table can reflect the statistical characteristics of wind speed and wind direction; the turbulence intensity table records the turbulence intensity under different wind speed and wind direction conditions. Turbulence intensity is usually expressed as the standard deviation of wind speed, reflecting the volatility of wind speed. In the approximate conversion, this table will be filled based on the wind speed standard deviation in the wind measurement data, and each wind speed bin and wind direction sector corresponds to a turbulence intensity value; the turbulence intensity relative standard deviation table provides the relative standard deviation of turbulence intensity, that is, the ratio of the turbulence intensity standard deviation to the average turbulence intensity. The relative standard deviation is a relative measure that can help understand the degree of fluctuation of turbulence intensity relative to the average level. Through this approximation, complex TAB format wind measurement data can be simplified into a tabular form that is easier to analyze and compare.

[0066] Furthermore, it also includes a method for calculating and obtaining equivalent wind data of a single tower at a result point according to different types of wind measurement data: When the wind measurement data is time series wind measurement data, the wind measurement data is calculated through equivalent wind data to obtain the single tower time series equivalent wind data of the result point when the wind measurement point is in any wind direction, as shown in the single tower time series equivalent wind data table of the result point; Results point single tower time series equivalent wind data table

[0067] When the wind measurement data is in TAB format, the wind speed range and wind direction sector are added to the wind measurement data to obtain TAB format subdivided wind measurement data, and the TAB format subdivided wind measurement data is calculated through equivalent wind data to obtain the TAB format subdivided equivalent wind data of the result point, and then all the TAB format subdivided equivalent wind data of the result point are aggregated to obtain the single-tower TAB format equivalent wind data of the result point when the wind measurement point is in any wind direction.

[0068] Here, the wind speed is [7,8) and the wind direction is [45,67.5). Assume that the frequency is , the turbulence intensity is , wind speed V is at intervals of 0.2m / s, wind direction Dir is at intervals of 2.5°, the subdivided frequency matrix is ​​the average division of the original grid data, and the increased wind speed range and wind direction sectors are displayed through 1 frequency distribution table + 1 turbulence intensity table + 1 turbulence intensity relative standard deviation table, without involving specific data: Frequency distribution table

[0069] Turbulence Intensity Table

[0070] Table of relative standard deviation of turbulence intensity

[0071] Furthermore, it also includes a method for obtaining equivalent wind data of multi-tower result points according to the single-tower equivalent wind data of the result point: According to the single-tower equivalent wind data of the result point, multi-tower weighted comprehensive calculation is performed to obtain the equivalent wind data of the multi-tower result point. The calculation formula of the equivalent wind data of the multi-tower result point is:

[0072]

[0073] In the formula, for Total number of wind towers, For the The confidence factor of the wind tower, The result point is based on The wind speed components and turbulence intensity calculated by the wind towers are Click to the result i The distance between the wind towers, For the i The weighted average coefficient of the wind towers is Multi-tower equivalent wind data for the result points.

[0074] Among them, time series data are integrated from multiple towers, and the result vectors (such as wind speed) are weighted by vectors, and scalars (such as turbulence intensity) are weighted by scalars.

[0075] As mentioned above, when the wind measurement data is time series wind measurement data, the multi-tower weighted comprehensive calculation can be directly performed on the time series equivalent wind data of the single tower of the result point; when the wind measurement data is TAB format wind measurement data, considering that the TAB format wind measurement data is an approximation of the time series wind measurement data, the error is related to the wind speed range and the size of the wind direction sector division. In order to reduce this approximate error, it is often further refined on the basis of the TAB format wind measurement data, and the equivalent wind data formula is used to calculate the result point TAB format subdivision equivalent wind data, and the result point TAB format subdivision equivalent wind data is obtained by multi-tower weighted comprehensive calculation of the result point TAB format subdivision equivalent wind data, and finally aggregated into the result point multi-tower TAB format equivalent wind data in the general sense (reference Figure 3 ).

[0076] The embodiment of the present invention is based on a CFD-based wind resource assessment and power generation prediction method, including equivalent wind data based on time series data and equivalent wind data based on TAB format data for a single wind tower and multiple wind towers. The equivalent wind data calculation is essentially a set of relationships established through directional calculation, which provides information such as wind acceleration factor, wind direction deviation, inflow angle, and turbulence intensity in specific wind directions of each sector for all result points in the calculation domain (including point elements of all analysis objects including wind measurement points, etc.), and maps the wind measurement data information of the wind measurement points to the analysis objects, so that the result points in the analysis objects have data in the same form as the wind measurement data, so as to provide more accurate, efficient, and reliable data support for wind resource assessment and power generation prediction through the equivalent wind data at all result points.

[0077] Based on the above-mentioned wind resource assessment and power generation prediction method, the present invention also discloses an assessment and prediction system, including: Directional calculation module, used to establish the CFD model of the target wind field, and obtain the wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation; A data processing module is used to obtain wind measurement data of a wind measurement point at any wind direction and determine two wind direction sectors adjacent to the wind direction of the wind measurement data; A weight calculation module, used to calculate the weights of two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data of the two wind direction sectors adjacent to the wind direction of the wind measurement data; The equivalent wind data calculation module is used to establish the wind condition mapping relationship between the wind measurement point and the result point, and calculate the equivalent wind data of the result point at any wind direction of the wind measurement point by combining the acquired wind measurement data, the calculated weight, and the wind condition data of the two wind direction sectors adjacent to the wind measurement data. The evaluation and prediction module is used to analyze the equivalent wind data of all result points, conduct wind resource evaluation and power generation prediction for the target wind farm.

[0078] Based on the above-mentioned wind resource assessment and power generation prediction method, the present invention also discloses a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the above-mentioned wind resource assessment and power generation prediction method are implemented.

[0079] Based on the above-mentioned wind resource assessment and power generation prediction method, the present invention also discloses a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the above-mentioned wind resource assessment and power generation prediction method when executing the computer program.

[0080] The present invention is described in terms of flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to specific embodiments. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0081] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0082] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0083] It should be understood that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. For those skilled in the art, the technical solutions described in the above embodiments can be modified, or some of the technical features therein can be replaced by equivalents; and all these modifications and replacements should fall within the scope of protection of the claims attached to the present invention.

Claims

1. A wind resource assessment and power generation prediction method based on CFD, characterized in that: The wind resource assessment and power generation prediction method comprises: Establishing a CFD model of the target wind field, and obtaining wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation of the CFD model; Obtaining wind measurement data of the wind measurement point at any wind direction, and determining two wind direction sectors adjacent to the wind direction of the wind measurement data; Calculate the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data corresponding to the wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data; Establish a wind condition mapping relationship between the wind measurement point and the result point, combine the acquired wind measurement data, the calculated weight, and the wind condition data corresponding to the wind measurement point and the result point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, and calculate and obtain the equivalent wind data of the result point at any wind direction of the wind measurement point; Analyze the equivalent wind data of all the result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

2. The CFD-based wind resource assessment and power generation prediction method according to claim 1, characterized in that: The obtaining of wind condition data of the wind measuring point and all result points within a set of wind direction sectors by directional calculation of the CFD model includes: A group of wind direction sectors having the same wind direction angle division range are established for the wind measurement point and each of the result points; Obtaining the wind resources and power generation of the CFD model, and obtaining the directional calculation results of the wind measurement point and each of the result points in a set of wind direction sectors through directional calculation, wherein the directional calculation results include wind speed components and turbulent kinetic energy; According to the directional calculation results of the wind measuring point and each of the result points, the wind condition data of the wind measuring point and each of the result points in each of the wind direction sectors are obtained by calculation, and the wind condition data include wind acceleration factor, inflow angle, wind direction deflection angle and calculated turbulence intensity.

3. The CFD-based wind resource assessment and power generation prediction method according to claim 2, characterized in that: The calculating the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition includes: The wind direction angle is extracted from the wind measurement data. According to the wind direction angle of the wind measurement data and the wind direction deflection angle corresponding to the wind measurement point in the two wind direction sectors adjacent to the wind direction of the wind measurement data, the weights of the two wind direction sectors adjacent to the wind direction of the wind measurement data are obtained by interpolation calculation. The weight calculation formula is: In the formula, Represents the weight of the first wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Represents the weight of the second wind direction sector adjacent to the wind direction of the wind measurement data, Indicates the wind direction angle of the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Indicates wind direction of wind measurement data i wind direction angle.

4. The CFD-based wind resource assessment and power generation prediction method according to claim 3, characterized in that: The equivalent wind data obtained by calculation includes the wind direction angle of the result point, and the calculation method of the wind direction angle of the result point includes: According to the wind direction angles of the two wind direction sectors adjacent to the wind direction of the wind measurement data and the wind direction deflection angle corresponding to the result point, an equivalent wind direction angle calculation formula for the result point under any wind direction is established. The equivalent wind direction angle calculation formula is: In the formula, Indicates any wind direction at the wind measurement point i The equivalent wind direction angle of the result point is Indicates the wind direction deflection angle of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind direction deflection angle of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data.

5. The CFD-based wind resource assessment and power generation prediction method according to claim 3, characterized in that: The equivalent wind data obtained by calculation include the wind speed and turbulence intensity at the result point; The calculation method of the wind speed at the result point includes: The wind speed is extracted from the wind measurement data, and when any wind direction of the wind measurement point is established, the mapping relationship between the wind speed of the wind measurement data of the wind measurement point and the directional wind acceleration factor is expressed as follows: When establishing the mapping relationship between the wind speed at the result point and the directional wind acceleration factor at any wind direction of the wind measurement point, the relationship expression is: By combining the above two relational expressions, a calculation formula for the equivalent wind speed of the result point at any wind direction of the wind measurement point is established. The equivalent wind speed calculation formula is: In the formula, Indicates the equivalent wind speed of the result point at any wind direction of the wind measurement point. Indicates the wind speed at the wind measuring point when the wind direction is arbitrary. Indicates the wind acceleration factor of the wind measurement point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the wind measurement point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data; Indicates the wind acceleration factor of the result point corresponding to the first wind direction sector adjacent to the wind direction of the wind measurement data. Indicates the wind acceleration factor of the result point corresponding to the second wind direction sector adjacent to the wind direction of the wind measurement data. Positive correlation; The calculation method of the turbulence intensity at the result point includes: Obtaining the measured environment turbulence intensity corresponding to the wind direction of the wind measurement data, and the interpolated calculated turbulence intensity of the wind measurement point in two adjacent wind direction sectors when the wind measurement point has any wind direction, and calculating the intensity difference between the measured environment turbulence intensity and the interpolated calculated turbulence intensity corresponding to the wind measurement point; When obtaining any wind direction at the wind measuring point, the turbulence intensity is calculated by interpolating the result point in two adjacent wind direction sectors; According to the wind acceleration factors of the wind measuring point and the result point in each wind direction sector, a ratio relationship between the wind measuring point and the result point corresponding to the wind acceleration factors when the wind measuring point has any wind direction is established, and the ratio relationship expression is: According to the intensity difference, the interpolated calculated turbulence intensity corresponding to the result point, and the ratio relationship, an equivalent turbulence intensity calculation formula for the result point under any wind direction is established. The equivalent turbulence intensity calculation formula is: In the formula, It represents the equivalent turbulence intensity of the result point at any wind direction of the wind measurement point. Indicates the interpolation calculation turbulence intensity corresponding to the result point, represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity corresponding to the wind measurement point, It indicates the wind acceleration factor at the wind measuring point at any wind direction. Indicates the wind acceleration factor of the result point at any wind direction of the wind measurement point. It indicates the ratio of the wind acceleration factors corresponding to the wind measurement point and the result point.

6. The CFD-based wind resource assessment and power generation prediction method according to claim 5, characterized in that: The calculating the intensity difference between the measured environment turbulence intensity and the interpolated calculated turbulence intensity corresponding to the wind measurement point comprises: The wind speed and the wind speed standard deviation are extracted from the wind measurement data, and the measured environmental turbulence intensity corresponding to the wind direction of the wind measurement data is calculated based on the wind speed and the wind speed standard deviation of the wind measurement data. The calculation formula for the measured environmental turbulence intensity is: In the formula, Indicates the measured environmental turbulence intensity of the wind direction corresponding to the wind measurement point. Indicates the standard deviation of wind speed at the wind measurement point corresponding to the wind direction of the wind measurement data, Indicates the wind speed at the wind measuring point corresponding to the wind direction of the wind measurement data; The wind speed component of the wind measuring point is extracted from the directional calculation result, and the turbulence intensity of the wind measuring point in two adjacent wind direction sectors is calculated when obtaining any wind direction of the wind measuring point. The calculation formula for the turbulence intensity of the wind measuring point is: In the formula, Indicates that the wind measuring point is in the adjacent wind direction sector when the wind direction is arbitrary i The calculated turbulence intensity is Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The turbulent kinetic energy corresponding to the wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from west to east at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from south to north at the corresponding wind measuring point is: Indicates the wind direction sector adjacent to the wind direction of the wind measurement data i The wind speed component from bottom to top corresponding to the wind measuring point; According to the calculated turbulence intensity of the wind measuring point in two adjacent wind direction sectors, the interpolated calculated turbulence intensity of the wind measuring point in two adjacent wind direction sectors is obtained by interpolation calculation. The calculation formula of the interpolated calculated turbulence intensity is: In the formula, It indicates the interpolation calculation of turbulence intensity at the wind measuring point in two adjacent wind direction sectors. It represents the calculated turbulence intensity of the wind measuring point in the first adjacent wind direction sector when the wind measuring point has any wind direction. It represents the calculated turbulence intensity of the wind measuring point in the adjacent second wind direction sector when the wind direction is arbitrary at the wind measuring point; The intensity difference is calculated according to the measured environment turbulence intensity obtained by calculation and the interpolation calculated turbulence intensity. The calculation formula of the intensity difference is: In the formula, Represents the intensity difference between the measured ambient turbulence intensity and the interpolated turbulence intensity.

7. The CFD-based wind resource assessment and power generation prediction method according to any one of claims 1 to 6, characterized in that: The method also includes obtaining different types of equivalent wind data according to the type of wind measurement data: Based on data collection, the time series wind measurement data at the wind measurement point is acquired, and the time series wind measurement data is calculated through equivalent wind data to obtain the time series equivalent wind data of the result point at the wind measurement point in any wind direction; Acquire TAB format wind measurement data based on file upload, calculate the TAB format wind measurement data through equivalent wind data, and obtain TAB format equivalent wind data of the result point at any wind direction of the wind measurement point; If the TAB format wind measurement data uses an intermediate value to represent the wind speed bin interval, the bin width of each wind speed bin interval in the TAB format wind measurement data is increased by a preset multiple, and then the processed TAB format wind measurement data is used to calculate the equivalent wind data; If the wind measurement data in the TAB format represents the wind speed bin interval with an upper limit value, the wind measurement data in the TAB format is used for equivalent wind data calculation; Compare and analyze the time series equivalent wind data and TAB format equivalent wind data of all the result points, and conduct wind resource assessment and power generation forecast for the target wind farm.

8. The CFD-based wind resource assessment and power generation prediction method according to claim 7, characterized in that: It also includes a method for calculating and obtaining the single-tower equivalent wind data at the result point according to different types of wind measurement data: When the wind measurement data is time series wind measurement data, the wind measurement data is calculated through equivalent wind data to obtain the single tower time series equivalent wind data of the result point when the wind measurement point is in any wind direction; When the wind measurement data is in TAB format, the wind speed range and wind direction sector are added to the wind measurement data to obtain TAB format subdivided wind measurement data, the TAB format subdivided wind measurement data is calculated through equivalent wind data to obtain the TAB format subdivided equivalent wind data of the result point, and then all the TAB format subdivided equivalent wind data of the result points are aggregated to obtain the single-tower TAB format equivalent wind data of the result point when the wind measurement point is in any wind direction.

9. The CFD-based wind resource assessment and power generation prediction method according to claim 8, characterized in that: It also includes a method for obtaining multi-tower equivalent wind data at the result point according to the single-tower equivalent wind data at the result point: According to the single-tower equivalent wind data of the result point, a multi-tower weighted comprehensive calculation is performed to obtain the multi-tower equivalent wind data of the result point. The calculation formula of the multi-tower equivalent wind data of the result point is: In the formula, for Total number of wind towers, For the The confidence factor of the wind tower, The result point is based on The wind speed components and turbulence intensity calculated by the wind towers are Click to the result The distance between the wind towers, For the The weighted average coefficient of the wind towers is Multi-tower equivalent wind data for the result points.

10. An evaluation and prediction system, using the CFD-based wind resource evaluation and power generation prediction method according to any one of claims 1 to 9, the system comprising: Directional calculation module, used to establish the CFD model of the target wind field, and obtain the wind condition data of the wind measurement point and all result points within a set of wind direction sectors through directional calculation; A data processing module is used to obtain wind measurement data of a wind measurement point at any wind direction and determine two wind direction sectors adjacent to the wind direction of the wind measurement data; A weight calculation module, used to calculate the weights of two wind direction sectors adjacent to the wind direction of the wind measurement data under each working condition according to the wind condition data of the two wind direction sectors adjacent to the wind direction of the wind measurement data; The equivalent wind data calculation module is used to establish the wind condition mapping relationship between the wind measurement point and the result point, and calculate the equivalent wind data of the result point at any wind direction of the wind measurement point by combining the acquired wind measurement data, the calculated weight, and the wind condition data of the two wind direction sectors adjacent to the wind measurement data. The evaluation and prediction module is used to analyze the equivalent wind data of all result points, and to conduct wind resource evaluation and power generation prediction for the target wind farm.

Images