Large-scale wind power plant flow field regional modeling calculation method, device and equipment and storage medium

Abstract

The invention provides a large-scale wind power plant flow field regional modeling calculation method, device and equipment and a storage medium, relates to the field of wind power generation, and aims to solve the problem that a large-scale wind power plant flow field calculation result is large in error. The method comprises the following steps: dividing a large-scale wind power plant into different areas; the areas comprise a free incoming flow area, a full-field blocking area, an entrance development area, a first transition area, a full development area, a second transition area and a group wake flow area; respectively carrying out modeling calculation on the areas to obtain the wind speed of each fan point in the wind power plant; the wind speed of the fan point location comprises a first exhaust fan point location wind speed, an inlet development area fan point location wind speed, a first transition area fan point location wind speed, a full development area fan point location wind speed and a second transition area fan point location wind speed; according to the wind turbine power curve and the wind speed of each wind turbine point location, the generated power of each wind turbine is obtained; and multiplying the generated power of each fan by the corresponding time to obtain the total generating capacity of the wind power plant.

Description

Technical Field

[0001] The embodiments of the present application relate to the field of wind power generation, and more specifically, to a method, device, equipment and storage medium for regional modeling and calculation of flow fields in large wind farms. Background Art

[0002] Guided by the dual carbon goals, my country's wind power has developed rapidly, and the construction of large-scale wind farms has increased. Large-scale wind farms have large areas, a large number of wind turbines and dense layout, which leads to extremely complex evolution of the internal flow field, which is manifested in different flow characteristics in different areas.

[0003] However, the engineering wake model widely used in related technologies cannot accurately describe the flow characteristics of different areas in large wind farms, resulting in large errors in the flow field calculation results of large wind farms, which in turn leads to large deviations between the power generation measurement results and the actual situation. Summary of the invention

[0004] The embodiments of the present application provide a method, device, equipment and storage medium for regional modeling and calculation of the flow field of a large wind farm, aiming to solve the problem that the engineering wake model widely used in the relevant technology cannot accurately describe the flow characteristics of different areas of a large wind farm, resulting in large errors in the calculation results of the flow field of a large wind farm.

[0005] A first aspect of an embodiment of the present application provides a method for calculating flow field regional modeling in a large wind farm, comprising:

[0006] Dividing a large wind farm into different areas; the areas include: free inflow area, full-field blocking area, inlet development area, first transition area, full development area, second transition area and group tail area;

[0007] Modeling and calculating the areas respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speeds of the wind turbine points include: the wind speed of the first row of wind turbines, the wind speed of the wind turbine points in the entrance development area, the wind speed of the wind turbine points in the first transition area, the wind speed of the wind turbine points in the fully developed area, and the wind speed of the wind turbine points in the second transition area;

[0008] The power generation of each wind turbine is obtained according to the wind turbine power curve and the wind speed at each wind turbine point; the total power generation of the wind farm is obtained by multiplying the power generation of each wind turbine by the corresponding time.

[0009] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0010] Modeling and calculating the full-field blocking area;

[0011] According to the following formula, the wind speed of the first row of fans is calculated:

[0012]

[0013] Among them, U1 represents the wind speed at the first row of fans; U0 represents the free flow velocity in the free flow area; x Indicates the average distance between adjacent fans in the flow direction; y represents the average distance between adjacent fans in the span direction; D represents the average impeller diameter of the fan; N rows Indicates the number of wind turbine rows in a wind farm.

[0014] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0015] Modeling and calculating the entrance development area;

[0016] The wind speed at the wind turbine point in the inlet development area is calculated according to the following formula:

[0017]

[0018] Among them, U 入口,i It represents the wind speed of the i-th fan point in the entrance development area, that is, the wind speed of the fan point in the entrance development area; U1 represents the wind speed of the first row of fans; △u ij represents the speed loss caused by the wake effect of the jth wind turbine at the i-th wind turbine location; C T,j represents the thrust coefficient of the jth fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; △x represents the flow distance between the j-th fan and the ith fan; △y represents the span distance between the j-th fan and the ith fan.

[0019] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0020] Modeling and calculating the fully developed area; the wind speed at the wind turbine point in the fully developed area is constant;

[0021] The wind speed at the wind turbine point in the fully developed area is calculated according to the following formula:

[0022]

[0023]

[0024] Among them, U3 represents the wind speed at the fan point in the fully developed area; U0 represents the free flow velocity in the free flow area; H ABL represents the height of the atmospheric boundary layer above the wind farm; z 0,lo Indicates the surface roughness; z 0,hirepresents the equivalent roughness of the wind farm; z h represents the hub height of the wind turbine; D represents the average impeller diameter of the wind turbine; ν w * represents the equivalent eddy viscosity coefficient.

[0025] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0026] Performing modeling calculation on the first transition zone;

[0027] The wind speed at the fan point in the first transition zone is calculated according to the following formula:

[0028]

[0029] Among them, U 过渡1 Represents the wind speed of the fan point in the first transition zone; U3 represents the wind speed of the fan point in the fully developed zone; U2 represents the wind speed of the point corresponding to P2; P2 represents the lower limit of the range of the first transition zone; P3 represents the upper limit of the range of the first transition zone; x represents the flow distance.

[0030] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0031] Performing modeling calculation on the second transition zone and the group wake zone;

[0032] The wind speed at the wind turbine point in the second transition zone and the wind speed at the group wake zone are calculated according to the following formula:

[0033]

[0034] Among them, U 过渡2 Indicates the wind speed at the fan point in the second transition zone; U 群尾流 It represents the wind speed at the point in the group wake zone; U6 represents the exit speed of the group wake zone; U4 represents the wind speed at the point in the last row of wind turbines included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; P6 represents the upper limit of the range of the group wake zone; x represents the flow distance.

[0035] In an optional implementation, the performing modeling calculations on the regions respectively includes:

[0036] Performing modeling calculation on the second transition zone;

[0037] The wind speed at the fan point in the second transition zone is calculated according to the following formula:

[0038]

[0039] Among them, U 过渡2U5 represents the wind speed at the point corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U4 represents the wind speed at the point of the last row of turbines included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; x represents the flow direction distance.

[0040] In an alternative embodiment, the separately modeling and calculating the regions includes:

[0041] Modeling and calculating the wake region of the group;

[0042] According to the following formula, calculate the wind speed at the point of the second transition zone turbines and the wind speed at the point of the wake region of the group:

[0043]

[0044] where, U 群尾流 represents the wind speed at the point of the wake region of the group; U5 represents the wind speed at the point corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U6 represents the outlet velocity of the wake region of the group; P6 represents the upper limit of the range of the wake region of the group; x represents the flow direction distance.

[0045] In an alternative embodiment, the dividing the large-scale wind farm into different regions includes:

[0046] Taking the position of the first row of turbines in the wind farm as the origin; the origin represents x = 0, corresponding to P1, where x represents the flow direction distance;

[0047] Determining the range of the free stream region as x < P0; the range of the free stream region is the entire upstream region of the wind farm;

[0048] Determining the range of the full-field blockage region as P0 < x < P1; the range of the full-field blockage region is 10D, that is, P1 - P0 = 10D, where D represents the average impeller diameter of the turbines;

[0049] Determining the range of the inlet development region as P1 < x < P2; the range of the inlet development region is (P2 - P1)α = H ABL , where α represents an empirical parameter; H ABL represents the height of the atmospheric boundary layer above the wind farm;

[0050] Determining the range of the first transition zone as P2 < x < P3; the range of the first transition zone is the distance between three rows of turbines;

[0051] Determining the range of the fully developed region as P3 < x < P4; the range of the fully developed region is the distance from the end of the range of the first transition zone to the distance between the third last row of turbines;

[0052] It is determined that the range of the second transition region is P4 < x < P5; the range of the second transition region is the distance between the penultimate row of wind turbines and the last row of wind turbines.

[0053] It is determined that the range of the wake region of the group is P5 < x < P6; the range of the wake region of the group is the downstream area of the entire wind farm.

[0054] In a second aspect of the embodiments of the present application, a large-scale wind farm flow field sub-region modeling and calculation device is provided, including:

[0055] A zoning module: used to divide a large-scale wind farm into different regions; the regions include: a free incoming flow region, a full-field blockage region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a wake region of the group.

[0056] A first calculation module: used to perform modeling calculations on the regions respectively to obtain the wind speed at each wind turbine point in the wind farm; the wind speed at the wind turbine point includes: the wind speed at the first row of wind turbine points, the wind speed at the wind turbine points in the inlet development region, the wind speed at the wind turbine points in the first transition region, the wind speed at the wind turbine points in the fully developed region, and the wind speed at the wind turbine points in the second transition region.

[0057] A second calculation module: used to obtain the power generation power of each wind turbine according to the wind turbine power curve and the wind speed at each wind turbine point; multiply the power generation power of each wind turbine by the corresponding time to obtain the total power generation of the wind farm.

[0058] In a third aspect of the embodiments of the present application, an electronic device is provided, including: a memory, a processor, and a computer program stored on the memory, and the processor executes the computer program to implement the steps of the large-scale wind farm flow field sub-region modeling and calculation method described in the first aspect.

[0059] In a fourth aspect of the embodiments of the present application, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the large-scale wind farm flow field sub-region modeling and calculation method described in the first aspect are implemented.

[0060] In the embodiments of the present disclosure, by dividing a large-scale wind farm into different regions, the regions include: a free incoming flow region, a full-field blockage region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a wake region of the group; performing modeling calculations on the regions respectively to obtain the wind speed at each wind turbine point in the wind farm; obtaining the power generation power of each wind turbine according to the wind turbine power curve and the wind speed at each wind turbine point; multiplying the power generation power of each wind turbine by the corresponding time to obtain the total power generation of the wind farm. By performing sub-region modeling calculations for the flow characteristics of different regions of a large-scale wind farm, the accuracy of calculating the flow field and power generation of the large-scale wind farm is improved, and the investment risk of the wind farm is reduced. Description of the Drawings

[0061] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the description of the embodiments of the present application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0062] Figure 1 This is a flowchart of the steps of a large-scale wind farm flow field regional modeling calculation method proposed in an embodiment of the present application;

[0063] Figure 2 It is a schematic diagram of the regional modeling and calculation of the flow field of a large wind farm proposed in an embodiment of the present application;

[0064] Figure 3 It is a schematic diagram of the wind speed distribution of each row of wind turbines in a large wind farm according to the regional modeling calculation method of the flow field of a large wind farm proposed in one embodiment of the present application;

[0065] Figure 4 is a schematic diagram of an electronic device in the present application. DETAILED DESCRIPTION

[0066] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0067] In the drawings, the size of the constituent elements, the thickness of the layer or the area may be exaggerated for the sake of clarity. Therefore, any implementation of the present disclosure is not necessarily limited to the size shown in the drawings, and the shapes and sizes of the components in the drawings do not reflect the true proportions. In addition, the drawings schematically show ideal examples, and any implementation of the present disclosure is not limited to the shapes or values ​​shown in the drawings.

[0068] Under the guidance of the dual-carbon goal, the wind power in China has developed rapidly, with the newly installed capacity increasing continuously. A large number of large-scale wind farms have been established or planned. The area of a large-scale wind farm is relatively large, and there are numerous and densely arranged wind turbines in the farm, resulting in extremely complex evolution of the internal flow field. Specifically, the flow field evolution in the wind farm shows different flow characteristics in different regions. According to the flow characteristics, the wind farm can be divided into free incoming flow area, full-field blockage area, inlet area, transition area, fully developed area, outlet area, group wake area, etc. However, the engineering wake models widely used in related technologies cannot accurately describe the flow characteristics of different regions in a large-scale wind farm, resulting in large errors in the calculation results of the flow field of the large-scale wind farm, and further leading to large deviations between the calculated power generation and the actual situation. Therefore, it is necessary to build sub-regional models for different regions in a large-scale wind farm to improve the accuracy of calculating the flow field and power generation of the wind farm.

[0069] Referring to Figure 1 , Figure 1 is the step flowchart of the sub-regional modeling calculation method for the flow field of a large-scale wind farm proposed in Embodiment 1 of the present application. As Figure 1 shown, the method includes the following steps S11 to S13:

[0070] Step S11: Divide the large-scale wind farm into different regions; the regions include: free incoming flow area, full-field blockage area, inlet development area, first transition area, fully developed area, second transition area, and group wake area.

[0071] According to the flow characteristics of different regions in the large-scale wind farm, the wind farm is divided into free incoming flow area, full-field blockage area, inlet development area, first transition area, fully developed area, second transition area, and group wake area. Taking the first row of wind turbines in the wind farm as the origin, according to the flow characteristics of different regions in the large-scale wind farm, the region is divided, as Figure 2 shown. The range of each region is related to factors such as the scale of the wind farm, the layout of the wind turbines, the atmospheric environment, and the characteristics of the wind turbines. By dividing the large-scale wind farm into regions, the wind flow characteristics of different regions in the wind farm can be grasped more accurately, providing strong support for the accurate calculation of the flow field and power generation of the large-scale wind farm.

[0072] Furthermore, the dividing the large-scale wind farm into different regions includes:

[0073] Taking the position of the first row of wind turbines in the wind farm as the origin; the origin represents x = 0, corresponding to P1, where x represents the flow direction distance;

[0074] Determining the range of the free incoming flow area as x < P0; the range of the free incoming flow area is the entire upstream area of the wind farm;

[0075] Determine that the range of the full-field blocked area is P0 < x < P1; the range of the full-field blocked area is 10D, that is, P1 - P0 = 10D, where D represents the average impeller diameter of the fan;

[0076] Determine that the range of the inlet development area is P1 < x < P2; the range of the inlet development area is (P2 - P1)α = H ABL , where α represents an empirical parameter; H ABL represents the height of the atmospheric boundary layer above the wind farm;

[0077] Determine that the range of the first transition area is P2 < x < P3; the range of the first transition area is the distance between three rows of fans;

[0078] Determine that the range of the fully developed area is P3 < x < P4; the range of the fully developed area is the distance from the end of the range of the first transition area to the distance between the third row of fans from the end;

[0079] Determine that the range of the second transition area is P4 < x < P5; the range of the second transition area is the distance between the third row of fans from the end and the last row of fans;

[0080] Determine that the range of the wake cluster area is P5 < x < P6; the range of the wake cluster area is the downstream area of the entire wind farm.

[0081] Taking the first row of fans in the wind farm as the origin (the position of x = 0, corresponding to point P1, and x is the flow direction distance), according to the flow characteristics of different regions of a large wind farm, regional division is carried out, as Figure 2 shown. The range of each region is related to factors such as the scale of the wind farm, the layout of the fans, the atmospheric environment, and the characteristics of the fans.

[0082] Free flow region: The flow field in this region is not affected by the wind farm, and the corresponding wind speed is the free flow wind speed of the wind farm. The range of the free flow region is the region where x < P0 upstream of the wind farm.

[0083] Full-field blocked area: The flow in this region is affected by the full-field blocking effect of the wind farm. The full-field blocking effect refers to the phenomenon that during the operation of all wind turbines in the wind farm, due to the action of the wind turbines themselves and each other, the incoming wind speed decreases. The range of the full-field blocked area is the region where P0 < x < P1 upstream of the wind farm. The specific regional range is related to factors such as the scale of the wind farm, the layout of the fans, the atmospheric environment, and the characteristics of the fans. The influence range of the full-field blocked area can reach 10D (D is the average impeller diameter of the fan), that is, P1 - P0 can reach 10D.

[0084] Inlet development area: The flow in this area is affected by the wake effect generated by each fan in the area, and the wakes between adjacent fans will be superimposed. The wake effect refers to the formation of a wake area with a reduced wind speed downstream while the fan extracts energy from the wind. The range of the inlet development area is P1 < x < P2, and the specific area range is related to factors such as the scale of the wind farm, the layout of the fans, the atmospheric environment, and the characteristics of the fans. Generally, it can be expressed by the following formula (1):

[0085] (P2 - P1)α = H ABL (1)

[0086] Where α represents the empirical parameter, usually taken as 0.8; H ABL represents the height of the atmospheric boundary layer above the wind farm.

[0087] First transition area: The transition area between the inlet development area and the fully developed area. The range of the first transition area is P2 < x < P3, and the range of the first transition area is set as the distance between three rows of fans.

[0088] Fully developed area: The interaction between the wind farm and the atmospheric boundary layer makes the flow in this area present a fully developed state. The range of the fully developed area is P3 < x < P4, and the range of the fully developed area is set as the distance from the end of the transition area to the third row of fans from the bottom.

[0089] Second transition area: The transition area between the fully developed area and the group wake area. The range of the second transition area is P4 < x < P5, and the range of the second transition area is set as the interval between three rows of fans, that is, the range from the third row of fans from the bottom to the last row of fans.

[0090] Group wake area: The downstream area of the entire wind farm, and the flow in this area is affected by the group wake of the entire wind farm. The range of the group wake area is P5 < x < P6, and the specific area range is related to factors such as the scale of the wind farm, the layout of the fans, the atmospheric environment, and the characteristics of the fans. The influence range of the group wake area can reach 50 km, that is, P6 - P5 can reach 50 km.

[0091] Step S12: Model and calculate each of the above-mentioned areas respectively to obtain the wind speed at each fan location in the wind farm; the wind speed at the fan location includes: the wind speed at the first row of fan locations, the wind speed at the fan locations in the inlet development area, the wind speed at the fan locations in the first transition area, the wind speed at the fan locations in the fully developed area, and the wind speed at the fan locations in the second transition area.

[0092] Modeling and calculation are respectively carried out on different divided regions to obtain the wind speeds at each fan position in the wind farm. The wind speeds at the fan positions include: the wind speed at the first row of fan positions, the wind speed at the fan positions in the inlet development area, the wind speed at the fan positions in the first transition area, the wind speed at the fan positions in the fully developed area, and the wind speed at the fan positions in the second transition area. Different methods are used to calculate the flow field in each region. When specifically calculating the flow field in a certain region, there are multiple choices for the calculation model. For example, for the inlet development area model, there are Jensen model, Frandsen model, Gaussian model, etc. For different wind farm environments, the corresponding model can be freely and flexibly selected to replace the original model, thereby further improving the calculation accuracy. It helps to improve the accuracy of the flow field and power generation calculation of large-scale wind farms.

[0093] Further, modeling and calculation are carried out on the full-field blocked area;

[0094] According to the following formula, calculate the wind speed at the first row of fan positions:

[0095]

[0096] where U1 represents the wind speed at the first row of fan positions; U0 represents the free flow velocity in the free flow area; △ x represents the average distance between adjacent fans in the flow direction; △ y represents the average distance between adjacent fans in the spanwise direction; D represents the average impeller diameter of the fan; N rows represents the number of rows of fans in the wind farm.

[0097] As Figure 2 shown, in the full-field blocked area, the range is P0 < x < P1, and the velocity decreases from U0 to U1 along the flow direction, where U0 is the free flow velocity in the free flow area; U1 is the actual inflow velocity of the wind farm affected by the full-field blocked area, that is, the wind speed at the first row of fan positions. According to the following formula (2), calculate the wind speed at the first row of fan positions:

[0098]

[0099] where U1 represents the wind speed at the first row of fan positions; U0 represents the free flow velocity in the free flow area; △ x represents the average distance between adjacent fans in the flow direction; △ y represents the average distance between adjacent fans in the spanwise direction; D represents the average impeller diameter of the fan; N rows represents the number of rows of fans in the wind farm.

[0100] The free incoming flow velocity of a wind farm can be obtained by deploying wind speed measurement devices in the wind farm, such as meteorological towers or wind measurement poles. These devices are usually equipped with various wind speed sensors (such as ultrasonic anemometers, anemometers), which can accurately measure the wind speed in the incoming flow. In this way, the free incoming flow velocity in the free incoming flow area can be obtained.

[0101] The flow direction usually refers to the main direction when air or gas flows through a group of wind turbines. The average distance between adjacent wind turbines in the flow direction refers to the average spacing between wind turbines arranged continuously along this direction.

[0102] The spanwise direction usually refers to the direction perpendicular to the flow direction, that is, the lateral width direction of the wind turbine array. The average distance between adjacent wind turbines in the spanwise direction refers to the average spacing between wind turbines arranged continuously along this direction.

[0103] Furthermore, the separately modeling and calculating of the regions includes:

[0104] Modeling and calculating the inlet development region;

[0105] According to the following formula, calculate the wind speed at the wind turbine position in the inlet development region:

[0106]

[0107] where, U 入口,i represents the wind speed at the i-th wind turbine position in the inlet development region, that is, the wind speed at the wind turbine position in the inlet development region; U1 represents the wind speed at the first row of wind turbine positions; △u ij represents the velocity deficit caused by the wake effect of the j-th wind turbine at the i-th wind turbine position; C T,j represents the thrust coefficient of the j-th wind turbine; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th wind turbine; △x represents the flow direction distance between the j-th wind turbine and the i-th wind turbine; △y represents the spanwise distance between the j-th wind turbine and the i-th wind turbine.

[0108] As Figure 2 shown, in the inlet development region, the range is P1 < x < P2, and the velocity decreases from U1 to U2 along the flow direction, where U1 is the wind speed at the first row of wind turbine positions; U2 is the wind speed at the position corresponding to P2, usually taking the wind speed at the last row of wind turbine positions included in the inlet development region. The wind speed at each wind turbine position in the region can be obtained by analyzing the wake model and the wake superposition model, and calculated according to the following formulas (3) and (4):

[0109]

[0110] where, U 入口,iIt represents the wind speed of the i-th fan point in the entrance development area, that is, the wind speed of the fan point in the entrance development area; U1 represents the wind speed of the first row of fans; △u ij represents the speed loss caused by the wake effect of the jth wind turbine at the i-th wind turbine location; C T,j represents the thrust coefficient of the jth fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; △x represents the flow distance between the j-th fan and the ith fan; △y represents the span distance between the j-th fan and the ith fan.

[0111] Speed ​​loss is the phenomenon of wind speed reduction caused by wind turbine rotation and wake effect. The thrust coefficient is a parameter that describes the efficiency of wind turbines in extracting energy from the wind. It is related to many factors such as wind turbine design, blade shape, wind speed, etc. The wake diffusion coefficient is an empirical parameter used to describe the rate at which the wake diffuses in space. Different wind field conditions, terrain and wind turbine layout may result in different wake diffusion coefficients.

[0112] Furthermore, the modeling and calculation of the regions respectively include:

[0113] Modeling and calculating the fully developed area; the wind speed at the wind turbine point in the fully developed area is constant;

[0114] The wind speed at the wind turbine point in the fully developed area is calculated according to the following formula:

[0115]

[0116]

[0117] Among them, U3 represents the wind speed at the fan point in the fully developed area; U0 represents the free flow velocity in the free flow area; H ABL represents the height of the atmospheric boundary layer above the wind farm; z 0,lo Indicates the surface roughness; z 0,hi represents the equivalent roughness of the wind farm; z h represents the hub height of the wind turbine; D represents the average impeller diameter of the wind turbine; ν w * represents the equivalent eddy viscosity coefficient.

[0118] like Figure 2As shown, in the fully developed region, where P3 < x < P4, the velocity decreases from U3 to U4 along the flow direction. Here, U3 is the wind speed at the fan position in the fully developed region; U4 is the wind speed at the last row of fan positions included in the fully developed region. Since the flow field in the fully developed region has reached a fully developed state, the flow field in this region is evenly distributed, and the velocity at each fan position is a constant value (U3 = U4), which can be obtained through the equivalent roughness model and calculated according to the following formulas (5) and (6):

[0119]

[0120] Among them, U3 represents the wind speed at the fan position in the fully developed region; U0 represents the free stream velocity in the free stream region; H ABL represents the atmospheric boundary layer height above the wind farm; z 0,lo represents the surface roughness; z 0,hi represents the equivalent roughness of the wind farm; z h represents the hub height of the wind turbine; D represents the average impeller diameter of the fan; ν w * represents the equivalent eddy viscosity coefficient.

[0121] Since the flow field in the fully developed region has reached a fully developed state, the flow field in this region is evenly distributed, and the velocity at each fan position is a constant value. Therefore, U3 can represent the wind speed at the fan position in the fully developed region. The atmospheric boundary layer height refers to the height of the near-surface air layer formed due to the ground friction in the atmosphere, generally between 400 and 1000 meters, but the specific height will be affected by various factors such as terrain, landform, and climate. The surface roughness refers to the magnitude of the frictional resistance generated by the ground surface to the air flow, which is mainly determined by the size and distribution of the surface roughness elements. The equivalent roughness of the wind farm refers to the magnitude of the equivalent frictional resistance generated by the entire wind farm to the air flow after considering the specific arrangement and interaction of the wind turbines in the wind farm. The hub height of the wind turbine refers to the vertical distance from the rotation center (i.e., the hub) of the wind turbine blade to the ground. Generally, the standard height of the wind turbine hub is about 80 - 120 meters.

[0122] Furthermore, the modeling and calculation of the region respectively include:

[0123] Modeling and calculating the first transition region;

[0124] Calculate the wind speed at the fan position in the first transition region according to the following formula:

[0125]

[0126] Among them, U 过渡1U1 represents the wind speed at the fan position in the first transition zone; U3 represents the wind speed at the fan position in the fully developed zone; U2 represents the wind speed at the position corresponding to P2; P2 represents the lower limit of the range of the first transition zone; P3 represents the upper limit of the range of the first transition zone; x represents the flow direction distance.

[0127] As Figure 2 shown, in the first transition zone, the range is P2 < x < P3, and the velocity decreases from U2 to U3 along the flow direction, where U2 is the wind speed at the position corresponding to P2, usually taking the wind speed at the last row of fan positions included in the inlet development zone; U3 is the wind speed at the fan position in the fully developed zone. The velocity U at each fan position in this area 过渡1 can be calculated according to the following formula (7):

[0128]

[0129] where U 过渡1 represents the wind speed at the fan position in the first transition zone; U3 represents the wind speed at the fan position in the fully developed zone; U2 represents the wind speed at the position corresponding to P2; P2 represents the lower limit of the range of the first transition zone; P3 represents the upper limit of the range of the first transition zone; x represents the flow direction distance.

[0130] Furthermore, the modeling and calculation of the regions respectively include:

[0131] Modeling and calculating the second transition zone and the wake region of the group;

[0132] According to the following formula, calculate the wind speed at the fan position in the second transition zone and the wind speed at the position in the wake region of the group:

[0133]

[0134] where U 过渡2 represents the wind speed at the fan position in the second transition zone; U 群尾流 represents the wind speed at the position in the wake region of the group; U6 represents the outlet velocity of the wake region of the group; U4 represents the wind speed at the last row of fan positions included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; P6 represents the upper limit of the range of the wake region of the group; x represents the flow direction distance.

[0135] As Figure 2 shown, in the second transition zone, the range is P4 < x < P5, and the velocity increases from U4 to U5 along the flow direction, where U5 is the wind speed at the position corresponding to P5, usually taking the wind speed at the last row of fan positions included in the second transition zone. In the wake region of the group, the range is P5 < x < P6, and the velocity further increases from U5 to U6 along the flow direction, where U6 is the outlet velocity of the wake region of the group. When U5 cannot be determined, the velocities in the second transition zone and the wake region of the group can be calculated according to the following formula (8):

[0136]

[0137] Among them, U 过渡2 represents the wind speed at the fan point in the second transition zone; U 群尾流 represents the wind speed at the point in the group wake zone; U6 represents the outlet speed of the group wake zone; U4 represents the wind speed at the fan point of the last row included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; P6 represents the upper limit of the range of the group wake zone; x represents the flow direction distance.

[0138] The free incoming flow velocity of the wind farm can be obtained by arranging wind speed measuring devices in the wind farm, such as meteorological towers or wind measurement poles. These devices are usually equipped with various wind speed sensors (such as ultrasonic anemometers, anemometers), which can accurately measure the wind speed in the incoming flow. In this way, the outlet speed of the group wake zone can be obtained. Usually, the outlet speed of the group wake zone is the same as the free incoming flow velocity of the free incoming flow zone, that is, U6 = U0.

[0139] Optionally, the separately modeling and calculating the regions includes:

[0140] Modeling and calculating the second transition zone;

[0141] Calculating the wind speed at the fan point in the second transition zone according to the following formula:

[0142]

[0143] Among them, U 过渡2 represents the wind speed at the fan point in the second transition zone; U5 represents the wind speed at the point corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U4 represents the wind speed at the fan point of the last row included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; x represents the flow direction distance.

[0144] As Figure 2 shown, in the second transition zone, the range is P4 < x < P5, and the speed increases from U4 to U5 along the flow direction. U5 is the wind speed at the point corresponding to P5, and usually takes the wind speed at the fan point of the last row included in the second transition zone. When U5 can be determined by calculation, the speed U 过渡2 at each fan point in the second transition zone can be calculated according to the following formula (9):

[0145]

[0146] Among them, U 过渡2 represents the wind speed at the fan point in the second transition zone; U5 represents the wind speed at the point corresponding to P5; U4 represents the wind speed at the fan point of the last row included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; P5 represents the upper limit of the range of the second transition zone; x represents the flow direction distance.

[0147] Optionally, the modeling calculations for the regions respectively include:

[0148] Performing a modeling calculation on the wake region of the group;

[0149] According to the following formula, calculating the wind speed at the fan position in the second transition region and the wind speed at the position in the wake region of the group:

[0150]

[0151] where U 群尾流 represents the wind speed at the position in the wake region of the group; U5 represents the wind speed at the position corresponding to P5; P5 represents the upper limit of the range of the second transition region; U6 represents the outlet speed of the wake region of the group; P6 represents the upper limit of the range of the wake region of the group; x represents the flow direction distance.

[0152] As Figure 2 shown, in the wake region of the group, the range is P5 < x < P6, and the speed further recovers and increases from U5 to U6 along the flow direction, where U6 is the outlet speed of the wake region of the group. When U5 can be determined by calculation, the speed in the wake region of the group can be calculated according to the following formula (10):

[0153]

[0154] where U 群尾流 represents the wind speed at the position in the wake region of the group; U6 represents the outlet speed of the wake region of the group; U5 represents the wind speed at the position corresponding to P5; P4 represents the lower limit of the range of the second transition region; P6 represents the upper limit of the range of the wake region of the group; x represents the flow direction distance.

[0155] Step S13: Obtain the power generation power of each wind turbine according to the power curve of the wind turbine and the wind speed at each fan position; multiply the power generation power of each wind turbine by the corresponding time to obtain the total power generation of the wind farm.

[0156] Each wind turbine has a specific power curve, which describes the power generation power that the wind turbine can generate at different wind speeds. Usually, this curve is provided by the wind turbine manufacturer and is drawn based on the design parameters and test data of the wind turbine. For each fan, according to its current wind speed value, find the corresponding power generation power on the power curve of the wind turbine. This usually involves interpolation calculation because the actual wind speed may not exactly correspond to a certain point on the power curve. Then determine the corresponding time used for calculating the total power generation, which can be hours, days, months, or years, etc. The choice of time depends on the specific requirements of the analysis and the availability of data. Finally, multiply the power generation power of each wind turbine by the corresponding time to obtain the total power generation of the wind farm. Quantify the power generation capacity of the wind farm to provide a basis for power dispatching and trading; by accurately calculating the power generation, the economic benefits of the wind farm can be evaluated and the investment risk of the wind farm can be reduced.

[0157] For ease of understanding, the embodiments of the present disclosure are described through specific application scenarios.

[0158] A large wind farm is equipped with 320 5.5MW wind turbines. The wind turbines are arranged regularly in the wind farm, with 16 rows along the main wind direction (flow direction) and 20 rows in the span direction. The average impeller diameter of the wind turbine is D = 158m, and the hub height of the wind turbine is z h =110m, average distance between adjacent fans in the flow direction△ x =1100m, the average distance between adjacent wind turbines in the span direction△ y =450m, the thrust coefficient of the fan in this case is C T =0.8, the atmospheric boundary layer height above the wind farm H ABL =1000m, surface roughness z 0,lo =0.001m, wind farm equivalent roughness z 0,hi =2.5m, free stream wind speed U0=8m / s, wake diffusion coefficient k w =0.05.

[0159] According to step S11, the large wind farm is divided into different areas, and the results are shown in Table 1.

[0160] Table 1. Results of wind farm division into different regions

[0161] area Range(m) Fan (exhaust) Free flow area x<-1580 Full-court blocking area -1580<x<0 1 Entrance Development Area 0<x<5623 2～6 First transition zone 5623<x<8800 7～9 Zone of Full Development 8800<x<13200 10～13 Second transition zone 13200<x<16500 14～16 Group wake area 16500<x<66500

[0162] According to step S12, modeling and calculation are performed on different areas respectively to obtain the wind speed of each wind turbine point in the wind farm. The results are shown in Table 2.

[0163] Table 2. Wind speed results for each fan location

[0164]

[0165]

[0166] At the same time, the wind speed distribution of each row of wind turbines in the wind farm can be obtained, such as Figure 3 Combined with the wind turbine power curve, the power generation of each wind turbine can be obtained, and the power generation of the entire wind farm can be calculated by multiplying it by the corresponding time.

[0167] Based on the same inventive concept, the embodiment of the present application discloses a large-scale wind farm flow field regional modeling and calculation device, which includes:

[0168] Zoning module: used to divide a large wind farm into different areas; the areas include: free inflow area, full-field blocking area, entrance development area, first transition area, full development area, second transition area and group tail area;

[0169] The first calculation module: used to perform modeling calculations on the respective regions to obtain the wind speeds at each fan position in the wind farm; the wind speeds at the fan positions include: the wind speeds at the first row of fan positions, the wind speeds at the fan positions in the inlet development zone, the wind speeds at the fan positions in the first transition zone, the wind speeds at the fan positions in the fully developed zone, and the wind speeds at the fan positions in the second transition zone;

[0170] The second calculation module: used to obtain the power generation of each fan according to the wind turbine power curve and the wind speed at each fan position; multiply the power generation of each fan by the corresponding time to obtain the total power generation of the wind farm.

[0171] In an alternative embodiment, the zoning module is specifically configured to:

[0172] Take the position of the first row of fans in the wind farm as the origin; the origin represents x = 0, corresponding to P1, where x represents the flow direction distance;

[0173] Determine that the range of the free stream region is x < P0; the range of the free stream region is the entire upstream region of the wind farm;

[0174] Determine that the range of the full-field blockage region is P0 < x < P1; the range of the full-field blockage region is 10D, that is, P1 - P0 = 10D, where D represents the average impeller diameter of the fan;

[0175] Determine that the range of the inlet development zone is P1 < x < P2; the range of the inlet development zone is (P2 - P1)α = H ABL , where α represents an empirical parameter; H ABL represents the height of the atmospheric boundary layer above the wind farm;

[0176] Determine that the range of the first transition zone is P2 < x < P3; the range of the first transition zone is the distance between three rows of fans;

[0177] Determine that the range of the fully developed zone is P3 < x < P4; the range of the fully developed zone is the distance from the end of the range of the first transition zone to the distance between the third row of fans from the end;

[0178] Determine that the range of the second transition zone is P4 < x < P5; the range of the second transition zone is the distance between the third row of fans from the end to the last row of fans;

[0179] Determine that the range of the wake interference region is P5 < x < P6; the range of the wake interference region is the entire downstream region of the wind farm.

[0180] In an alternative embodiment, the first calculation module is specifically configured to:

[0181] Modeling and calculating the full-field blocking area;

[0182] According to the following formula, the wind speed of the first row of fans is calculated:

[0183]

[0184] Among them, U1 represents the wind speed at the first row of fans; U0 represents the free flow velocity in the free flow area; x Indicates the average distance between adjacent fans in the flow direction; y represents the average distance between adjacent fans in the span direction; D represents the average impeller diameter of the fan; N rows Indicates the number of wind turbine rows in a wind farm.

[0185] In an optional implementation, the first calculation module is specifically configured to:

[0186] Modeling and calculating the entrance development area;

[0187] The wind speed at the wind turbine point in the inlet development area is calculated according to the following formula:

[0188]

[0189] Among them, U 入口,i It represents the wind speed of the i-th fan point in the entrance development area, that is, the wind speed of the fan point in the entrance development area; U1 represents the wind speed of the first row of fans; △u ij represents the speed loss caused by the wake effect of the jth wind turbine at the i-th wind turbine location; C T,j represents the thrust coefficient of the jth fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; △x represents the flow distance between the j-th fan and the ith fan; △y represents the span distance between the j-th fan and the ith fan.

[0190] In an optional implementation, the first calculation module is specifically configured to:

[0191] Modeling and calculating the fully developed area; the wind speed at the wind turbine point in the fully developed area is constant;

[0192] The wind speed at the wind turbine point in the fully developed area is calculated according to the following formula:

[0193]

[0194] Among them, U3 represents the wind speed at the fan point in the fully developed area; U0 represents the free flow speed in the free flow area; H ABL represents the height of the atmospheric boundary layer above the wind farm; z 0,loIndicates the surface roughness; z 0,hi represents the equivalent roughness of the wind farm; z h represents the hub height of the wind turbine; D represents the average impeller diameter of the wind turbine; ν w * represents the equivalent eddy viscosity coefficient.

[0195] In an optional implementation, the first calculation module is specifically configured to:

[0196] Performing modeling calculation on the first transition zone;

[0197] The wind speed at the fan point in the first transition zone is calculated according to the following formula:

[0198]

[0199] Among them, U 过渡1 Represents the wind speed of the fan point in the first transition zone; U3 represents the wind speed of the fan point in the fully developed zone; U2 represents the wind speed of the point corresponding to P2; P2 represents the lower limit of the range of the first transition zone; P3 represents the upper limit of the range of the first transition zone; x represents the flow distance.

[0200] In an optional implementation, the first calculation module is specifically configured to:

[0201] Performing modeling calculation on the second transition zone and the group wake zone;

[0202] The wind speed at the wind turbine point in the second transition zone and the wind speed at the group wake zone are calculated according to the following formula:

[0203]

[0204] Among them, U 过渡2 Indicates the wind speed at the fan point in the second transition zone; U 群尾流 It represents the point wind speed in the group wake area; U6 represents the exit speed in the group wake area; U4 represents the point wind speed of the last row of wind turbines in the fully developed area; P4 represents the lower limit of the second transition area; P6 represents the upper limit of the group wake area; x represents the flow distance.

[0205] In an optional implementation, the first calculation module is specifically configured to:

[0206] Performing modeling calculation on the second transition zone;

[0207] The wind speed at the fan point in the second transition zone is calculated according to the following formula:

[0208]

[0209] Among them, U 过渡2Represents the wind speed of the fan point in the second transition zone; U5 represents the wind speed of the point corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U4 represents the wind speed of the last row of fans included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; x represents the flow distance.

[0210] In an optional implementation, the first calculation module is specifically configured to:

[0211] Modeling and calculating the group wake area;

[0212] The wind speed at the wind turbine point in the second transition zone and the wind speed at the group wake zone are calculated according to the following formula:

[0213]

[0214] Among them, U 群尾流 Represents the point wind speed in the group wake zone; U5 represents the point wind speed corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U6 represents the exit speed of the group wake zone; P6 represents the upper limit of the range of the group wake zone; x represents the flow direction distance.

[0215] The method for regional modeling and calculation of the flow field of a large wind farm provided by the embodiment of the present disclosure is to divide the large wind farm into different areas, including: free flow area, full-field blocking area, entrance development area, first transition area, full development area, second transition area and group tail area; modeling and calculation are performed on each area to obtain the wind speed of each wind turbine point in the wind farm; the power generation of each wind turbine is obtained according to the wind turbine power curve and the wind speed of each wind turbine point; the power generation of each wind turbine is multiplied by the corresponding time to obtain the total power generation of the wind farm. Regional modeling and calculation are performed according to the flow characteristics of different areas of the large wind farm, so as to improve the accuracy of the calculation of the flow field and power generation of the large wind farm and reduce the investment risk of the wind farm.

[0216] It should be noted that, for the method embodiments, for the sake of simplicity, they are all described as a series of action combinations, but those skilled in the art should be aware that the embodiments of the present disclosure are not limited by the order of the actions described, because according to the embodiments of the present disclosure, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions involved are not necessarily required by the embodiments of the present disclosure.

[0217] The present disclosure also provides an electronic device, referring to Figure 4 , Figure 4 Schematic diagram of an electronic device according to an embodiment of the present disclosure. Figure 4As shown, the electronic device 100 includes: a memory 110 and a processor 120. The memory 110 and the processor 120 are connected via a bus communication. A computer program is stored in the memory 110. The computer program can be run on the processor 120 to implement the steps in the large wind farm flow field regional modeling calculation method disclosed in the embodiment of the present disclosure.

[0218] The embodiment of the present disclosure also provides a computer-readable storage medium. When the instructions in the computer-readable storage medium are executed by a processor of a computer device, the computer device can execute the steps in the large-scale wind farm flow field regional modeling calculation method disclosed in the embodiment of the present disclosure.

[0219] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0220] The embodiments of the present application are described with reference to the flowcharts and / or block diagrams of the methods, apparatuses, electronic devices, and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the 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 terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal 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.

[0221] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including 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.

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

[0223] Although the preferred embodiments of the present application have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the embodiments of the present application.

[0224] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of further restrictions, the elements defined by the sentence "including one..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

[0225] The above is a detailed introduction to the large-scale wind farm flow field regional modeling calculation method, device, equipment and storage medium provided by the present application. This article uses specific examples to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea; at the same time, for general technical personnel in this field, according to the idea of ​​the present application, there will be changes in the specific implementation method and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A method for modeling and calculating the flow field of a large wind farm in different regions, characterized in that: Including: Dividing a large-scale wind farm into different regions; The regions include: free incoming flow region, full-field blockage region, inlet development region, first transition region, fully developed region, second transition region, and wake region of the group; Performing modeling calculations on the regions respectively to obtain the wind speed at each wind turbine location in the wind farm; the wind speed at the wind turbine location includes: wind speed at the first row of wind turbine locations, wind speed at the inlet development region wind turbine locations, wind speed at the first transition region wind turbine locations, wind speed at the fully developed region wind turbine locations, and wind speed at the second transition region wind turbine locations; According to the wind turbine power curve and the wind speed at each wind turbine location, obtaining the power generation of each wind turbine; multiplying the power generation of each wind turbine by the corresponding time to obtain the total power generation of the wind farm.

2. The method for calculating flow field in different regions of a large wind farm according to claim 1 is characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the full-field blockage region; Calculating the wind speed at the first row of wind turbine locations according to the following formula: Among them, U1 represents the wind speed at the first row of fans; U0 represents the free flow velocity in the free flow area; x Indicates the average distance between adjacent fans in the flow direction; y represents the average distance between adjacent fans in the span direction; D represents the average impeller diameter of the fan; N rows Indicates the number of wind turbine rows in a wind farm.

3. The method for calculating flow field in different regions of a large wind farm according to claim 2 is characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the inlet development region; Calculating the wind speed at the inlet development region wind turbine locations according to the following formula: Among them, U 入口,i It represents the wind speed of the i-th fan point in the entrance development area, that is, the wind speed of the fan point in the entrance development area; U1 represents the wind speed of the first row of fans; △u ij represents the speed loss caused by the wake effect of the jth wind turbine at the i-th wind turbine location; C T,j represents the thrust coefficient of the jth fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; △x represents the flow distance between the j-th fan and the ith fan; △y represents the span distance between the j-th fan and the ith fan.

4. The method for calculating flow field regional modeling in a large wind farm according to claim 3 is characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the fully developed region; the wind speed at the fully developed region wind turbine locations is constant; Calculating the wind speed at the fully developed region wind turbine locations according to the following formula: Among them, U3 represents the wind speed at the fan point in the fully developed area; U0 represents the free flow speed in the free flow area; H ABL represents the height of the atmospheric boundary layer above the wind farm; z 0,lo Indicates the surface roughness; z 0,hi represents the equivalent roughness of the wind farm; z h represents the hub height of the wind turbine; D represents the average impeller diameter of the wind turbine; ν w * represents the equivalent eddy viscosity coefficient.

5. The method for calculating flow field regional modeling in a large wind farm according to claim 4, characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the first transition region; Calculating the wind speed at the first transition region wind turbine locations according to the following formula: Among them, U 过渡1 Represents the wind speed of the fan point in the first transition zone; U3 represents the wind speed of the fan point in the fully developed zone; U2 represents the wind speed of the point corresponding to P2; P2 represents the lower limit of the range of the first transition zone; P3 represents the upper limit of the range of the first transition zone; x represents the flow distance.

6. The method for calculating flow field in different regions of a large wind farm according to claim 5, characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the second transition region and the wake region of the group; Calculating the wind speed at the second transition region wind turbine locations and the wind speed at the wake region of the group locations according to the following formula: Among them, U 过渡2 Indicates the wind speed at the fan point in the second transition zone; U 群尾流 It represents the wind speed at the point in the group wake zone; U6 represents the exit speed of the group wake zone; U4 represents the wind speed at the point in the last row of wind turbines included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; P6 represents the upper limit of the range of the group wake zone; x represents the flow distance.

7. The method for calculating flow field in different regions of a large wind farm according to claim 5, characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the second transition region; Calculating the wind speed at the second transition region wind turbine locations according to the following formula: Among them, U 过渡2 Represents the wind speed of the fan point in the second transition zone; U5 represents the wind speed of the point corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U4 represents the wind speed of the last row of fans included in the fully developed zone; P4 represents the lower limit of the range of the second transition zone; x represents the flow distance.

8. The method for calculating flow field in different regions of a large wind farm according to claim 7, characterized in that: The performing modeling calculations on the regions respectively includes: Performing modeling calculations on the wake region of the group; Calculating the wind speed at the second transition region wind turbine locations and the wind speed at the wake region of the group locations according to the following formula: Among them, U 群尾流 Represents the point wind speed in the group wake zone; U5 represents the point wind speed corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U6 represents the exit speed of the group wake zone; P6 represents the upper limit of the range of the group wake zone; x represents the flow direction distance.

9. The method for calculating flow field regional modeling in a large wind farm according to any one of claims 2 to 8, characterized in that: The dividing the large-scale wind farm into different regions includes: Taking the position of the first row of wind turbines in the wind farm as the origin; the origin represents x = 0, corresponding to P1, where x represents the flow direction distance; Determining the range of the free incoming flow region as x < P0; the range of the free incoming flow region is the entire upstream region of the wind farm; Determining the range of the full-field blockage region as P0 < x < P1; the range of the full-field blockage region is 10D, that is, P1 - P0 = 10D, where D represents the average impeller diameter of the wind turbine; Determine that the range of the inlet development area is P1 < x < P2; the range of the inlet development area is (P2 - P1)α = H ABL , where α represents an empirical parameter; H ABL represents the height of the atmospheric boundary layer above the wind farm; Determining the range of the first transition region as P2 < x < P3; the range of the first transition region is the distance between three rows of wind turbines; Determining the range of the fully developed region as P3 < x < P4; the range of the fully developed region is the distance from the end of the range of the first transition region to the distance between the third last row of wind turbines; Determining the range of the second transition region as P4 < x < P5; the range of the second transition region is the distance between the third last row of wind turbines and the last row of wind turbines; It is determined that the range of the group wake region is P5 < x < P6; the range of the group wake region is the downstream area of the entire wind farm.

10. A large-scale wind farm flow field regional modeling and calculation device, characterized in that: The device includes: A zoning module: used to divide a large-scale wind farm into different regions; the regions include: a free incoming flow region, a full-field blockage region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a group wake region; A first calculation module: used to perform modeling calculations on the regions respectively to obtain the wind speeds at each fan position in the wind farm; the wind speeds at the fan positions include: the wind speeds at the first row of fan positions, the wind speeds at the fan positions in the inlet development region, the wind speeds at the fan positions in the first transition region, the wind speeds at the fan positions in the fully developed region, and the wind speeds at the fan positions in the second transition region; A second calculation module: used to obtain the power generation power of each fan according to the wind turbine power curve and the wind speed at each fan position; multiplying the power generation power of each fan by the corresponding time to obtain the total power generation of the wind farm.

11. An electronic device, characterized in that: It includes: A processor, a memory, and a computer program stored on the memory and executable on the processor. When the processor executes the program, it implements the large-scale wind farm flow field sub-region modeling calculation method described in any one of claims 1 to 9.

12. A computer-readable storage medium, characterized in that: When the instructions in the storage medium are executed by the processor of the electronic device, the electronic device can execute the large-scale wind farm flow field sub-region modeling calculation method described in any one of claims 1 to 9.

Images