Large wind farm flow field sub-region modeling calculation method, device and equipment and storage medium

Abstract

The application provides a large wind farm flow field sub-region modeling calculation method, device and equipment and storage medium, relates to the wind power generation field, and aims to solve the problem of large wind farm flow field calculation result error. The method comprises the following steps: dividing a large wind farm into different regions; the regions comprise: a free incoming flow region, a full-field blocking region, an inlet development region, a first transition region, a fully developed region, a second transition region and a group wake region; modeling and calculating the regions respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point comprises: the first row of wind turbine point wind speed, the inlet development region wind turbine point wind speed, the first transition region wind turbine point wind speed, the fully developed region wind turbine point wind speed and the second transition region wind turbine point wind speed; according to the wind turbine power curve and the wind speed of each wind turbine point, the power generation power of each wind turbine is obtained; the power generation power of each wind turbine is multiplied by the corresponding time to obtain the total power generation of the wind farm.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to the field of wind power generation, in particular, to a large wind farm flow field sub-region modeling calculation method, device, equipment and storage medium. BACKGROUND

[0002] China's wind power development is rapid, and large wind farm construction is increasing. Large wind farm area is large, the number of wind turbines in the field is large and densely arranged, which leads to extremely complex internal flow field evolution and development, specifically, different regions present different flow characteristics.

[0003] However, the engineering wake model widely used in the related art cannot accurately describe the flow characteristics of different regions of a large wind farm, resulting in large errors in the flow field calculation results of the large wind farm, and further leading to large deviations between the actual situation and the power generation calculation results. SUMMARY

[0004] Embodiments of the present application provide a large wind farm flow field sub-region modeling calculation method, device, equipment and storage medium, aiming to solve the problem that the engineering wake model widely used in the related art cannot accurately describe the flow characteristics of different regions of a large wind farm, resulting in large errors in the flow field calculation results of the large wind farm.

[0005] The first aspect of the embodiments of the present application provides a large wind farm flow field sub-region modeling calculation method, comprising:

[0006] The large wind farm is divided into different regions; the regions include: a free flow region, a full-field blocked region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a group wake region;

[0007] Modeling and calculating the regions respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point includes: the first row of wind turbine point wind speed, the inlet development region wind turbine point wind speed, the first transition region wind turbine point wind speed, the fully developed region wind turbine point wind speed, and the second transition region wind turbine point wind speed;

[0008] According to the wind turbine power curve and the wind speed of each wind turbine point, the power generation power of each wind turbine is obtained; the power generation power of each wind turbine is multiplied by the corresponding time to obtain the total power generation of the wind farm.

[0009] In an optional implementation, the modeling and calculating the regions respectively includes:

[0010] Modeling and calculating the full-field blocked region;

[0011] The first row of wind turbine point wind speed is calculated according to the following formula:

[0012]

[0013] wherein, U1 represents the first row of wind turbine point speed; U0 represents the free flow velocity of the free flow area; △u x represents the average distance between the adjacent wind turbines in the flow direction; △x y represents the average distance between the adjacent wind turbines in the spanwise direction; D represents the average impeller diameter of the wind turbine; N rows represents the row number of the wind turbine in the wind farm.

[0014] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0015] modeling and calculating the inlet development region;

[0016] calculating the wind turbine point speed of the inlet development region according to the following formula:

[0017]

[0018] wherein, U 入口,i represents the wind speed of the i th wind turbine point in the inlet development region, i.e. the wind turbine point speed of the inlet development region; U1 represents the first row of wind turbine point speed; △u ij represents the speed loss caused by the wake effect of the j th wind turbine at the i th wind turbine point; 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 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.

[0019] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0020] modeling and calculating the fully developed region; the wind turbine point speed of the fully developed region is constant;

[0021] calculating the wind turbine point speed of the fully developed region according to the following formula:

[0022]

[0023]

[0024] wherein, U3 represents the wind turbine point speed of the fully developed region; U0 represents the free flow velocity of the free flow area; 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; zh H represents the hub height of the wind turbine; D represents the average rotor diameter of the wind turbine; v represents the wind speed at the hub height of the wind turbine; x represents the streamwise distance. w * represents the equivalent eddy viscosity coefficient.

[0025] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0026] modeling and calculating the first transition region;

[0027] calculating the wind speed at the wind turbine point of the first transition region according to the following formula:

[0028]

[0029] wherein, U 过渡1 represents the wind speed at the wind turbine point of the first transition region; U3 represents the wind speed at the wind turbine point of the fully developed region; U2 represents the wind speed at the point corresponding to P2; P2 represents the lower limit of the range of the first transition region; P3 represents the upper limit of the range of the first transition region; and x represents the streamwise distance.

[0030] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0031] modeling and calculating the second transition region and the wake region;

[0032] calculating the wind speed at the wind turbine point of the second transition region and the wind speed at the point of the wake region according to the following formula:

[0033]

[0034] wherein, U 过渡2 represents the wind speed at the wind turbine point of the second transition region; U 群尾流 represents the wind speed at the point of the wake region; U6 represents the outlet velocity of the wake region; U4 represents the wind speed at the wind turbine point of the last row included in the fully developed region; 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; and x represents the streamwise distance.

[0035] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0036] modeling and calculating the second transition region;

[0037] calculating the wind speed at the wind turbine point of the second transition region according to the following formula:

[0038]

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

[0040] In an alternative embodiment, the modeling and calculating of the regions respectively comprises:

[0041] modeling and calculating the wake zone;

[0042] The point speed of the second transition zone and the point speed of the wake zone are calculated according to the following formula:

[0043]

[0044] wherein U 群尾流 represents the point speed of the wake zone; U5 represents the point speed corresponding to P5; P5 represents the upper limit of the second transition zone; U6 represents the outlet speed of the wake zone; P6 represents the upper limit of the wake zone; and x represents the flow distance.

[0045] In an alternative embodiment, the division of the large wind farm into different regions comprises:

[0046] The position of the first row of wind turbines in the wind farm is taken as the origin; the origin represents x=0, corresponding to P1, wherein x represents the flow distance;

[0047] The range of the free flow zone is determined to be x

[0048] The range of the full-field blockage zone is determined to be P0

[0049] The range of the inlet development zone is determined to be P1 ABL wherein a represents an empirical parameter; and H ABL represents the atmospheric boundary layer height above the wind farm;

[0050] The range of the first transition zone is determined to be P2

[0051] The range of the fully developed zone is determined to be P3

[0052] determining that the range of the second transition zone is P4

[0053] determining that the range of the group wake zone is P5

[0054] The embodiment of the application provides a large wind farm flow field sub-region modeling and calculating device, which comprises:

[0055] The partition module is used for dividing the large wind farm into different regions; the regions comprise a free incoming flow zone, a full-field blocking zone, an inlet development zone, a first transition zone, a fully developed zone, a second transition zone and a group wake zone.

[0056] The first calculating module is used for modeling and calculating the regions respectively, so as to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point comprises a first-row wind turbine point wind speed, an inlet development zone wind turbine point wind speed, a first transition zone wind turbine point wind speed, a fully developed zone wind turbine point wind speed and a second transition zone wind turbine point wind speed.

[0057] The second calculating module is used for obtaining the power generation of each wind turbine according to the wind turbine power curve and the wind speed of each wind turbine point; and multiplying the power generation of each wind turbine by the corresponding time, so as to obtain the total power generation of the wind farm.

[0058] In the third aspect of the embodiment of the application, an electronic device is provided, which comprises a memory, a processor and a computer program stored in the memory; the processor executes the computer program to realize the steps of the large wind farm flow field sub-region modeling and calculating method in the first aspect.

[0059] In the fourth aspect of the embodiment of the application, a computer readable storage medium is provided, which stores a computer program; the computer program is executed by a processor to realize the steps of the large wind farm flow field sub-region modeling and calculating method in the first aspect.

[0060] In the embodiment of the disclosure, the large wind farm is divided into different regions, the regions comprise a free incoming flow zone, a full-field blocking zone, an inlet development zone, a first transition zone, a fully developed zone, a second transition zone and a group wake zone; the regions are modeled and calculated respectively, so as 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; and the power generation of each wind turbine is multiplied by the corresponding time, so as to obtain the total power generation of the wind farm. The flow characteristics of different regions of the large wind farm are modeled and calculated, the accuracy of the large wind farm flow field and the power generation calculation is improved, and the investment risk of the wind farm is reduced. BRIEF DESCRIPTION OF DRAWINGS

[0061] To more clearly illustrate the technical solutions of the embodiments of this application, the drawings used in the description of the embodiments of this application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0062] Figure 1 This is a flowchart illustrating the steps of a method for regional modeling and calculation of flow field in a large wind farm, as proposed in an embodiment of this application.

[0063] Figure 2 This is a schematic diagram of the method for regional modeling and calculation of flow field in a large wind farm, as proposed in an embodiment of this application.

[0064] Figure 3 This is a schematic diagram of the wind speed distribution at each row of wind turbines in a large wind farm, based on the regional modeling and calculation method for the flow field in a large wind farm proposed in an embodiment of this application.

[0065] Figure 4 This is a schematic diagram of an electronic device according to this application. Detailed Implementation

[0066] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0067] In the accompanying drawings, the size of constituent elements, the thickness of layers, or areas may sometimes be exaggerated for clarity. Therefore, any implementation of this disclosure is not necessarily limited to the dimensions shown in the drawings, and the shapes and sizes of the components in the drawings do not reflect true proportions. Furthermore, the drawings schematically illustrate ideal examples, and any implementation of this disclosure is not limited to the shapes or values ​​shown in the drawings.

[0068] The wind power in China develops rapidly, and the newly added installed capacity is increasing. A large number of large-scale wind farms have been established or planned. The large-scale wind farm has a large regional range, and the number of wind turbines in the field is large and densely arranged, which leads to the extremely complex evolution and development of the internal flow field. The flow field development and evolution of the wind farm presents different flow characteristics in different regions. According to the flow characteristics, the wind farm can be divided into a free flow area, a full-field blocking area, an inlet area, a transition area, a fully developed area, an outlet area, and a group wake area. However, the widely used engineering wake model in the related art cannot accurately describe the flow characteristics of different regions of the large-scale wind farm, resulting in a large error in the calculation results of the flow field of the large-scale wind farm, and further leading to a large deviation between the actual situation and the calculation results of the power generation capacity. Therefore, it is necessary to perform regional modeling according to the flow characteristics of different regions of the large-scale wind farm, so as to improve the accuracy of the calculation of the flow field and the power generation capacity of the wind farm.

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

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

[0071] According to the flow characteristics of different regions of the large-scale wind farm, the wind farm is divided into a free flow area, a full-field blocking area, an inlet development area, a first transition area, a fully developed area, a second transition area, and a group wake area. Taking the first row of wind turbines of the wind farm as the origin, the regions are divided according to the flow characteristics of different regions of the large-scale wind farm, as Figure 2 indicated. The range of each region is related to 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 more accurately mastered, which provides strong support for the accurate calculation of the flow field and the power generation capacity of the large-scale wind farm.

[0072] Further, the large-scale wind farm is divided into different regions, including:

[0073] The position of the first row of wind turbines of the wind farm is taken as the origin; the origin represents x=0, corresponding to P1, wherein x represents the flow distance;

[0074] The range of the free flow area is determined to be x

[0075] determining that the full-field blockage region ranges from P0 to P1; the full-field blockage region ranges by 10D, that is, P1-P0=10D, wherein D represents the average impeller diameter of the wind turbine;

[0076] determining that the inlet development region ranges from P1 to P2; the inlet development region ranges by (P2-P1)α=H ABL , wherein α represents an empirical parameter; H ABL represents the atmospheric boundary layer height above the wind farm;

[0077] determining that the first transition region ranges from P2 to P3; the first transition region ranges by the distance between three rows of wind turbines;

[0078] determining that the fully developed region ranges from P3 to P4; the fully developed region ranges by the distance between the end of the first transition region and the third-to-last row of wind turbines;

[0079] determining that the second transition region ranges from P4 to P5; the second transition region ranges by the distance between the third-to-last row of wind turbines and the last row of wind turbines;

[0080] determining that the wake region of the group ranges from P5 to P6; the wake region of the group ranges by the entire downstream region of the wind farm.

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

[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 free flow region ranges from the region upstream of the wind farm x

[0083] Full-field blockage region: the flow in this region is affected by the full-field blockage effect of the wind farm. The full-field blockage effect refers to the phenomenon that the incoming flow wind speed is reduced due to the action of the wind turbine itself and the interaction between the wind turbines in the wind farm during operation. The full-field blockage region ranges from P0 to P1 upstream of the wind farm, and the specific range of the region is related to the scale of the wind farm, the layout of the wind turbines, the atmospheric environment, and the characteristics of the wind turbines. The influence range of the full-field blockage region can reach 10D (D is the average impeller diameter of the wind turbine), that is, P1-P0 can reach 10D.

[0084] Inlet development zone: the flow in this zone is affected by the wake effect of each wind turbine in the zone, and the wake of adjacent wind turbines is superimposed. The wake effect refers to the wind turbine obtaining energy from the wind while forming a wake zone with reduced wind speed downstream. The range of the inlet development zone is P1 < x < P2, and the specific range is related to the scale of the wind farm, the layout of the wind turbines, the atmospheric environment, and the characteristics of the wind turbines. Generally, it can be expressed by the following formula (1):

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

[0086] wherein a represents an empirical parameter, usually 0.8; H ABL represents the height of the atmospheric boundary layer above the wind farm.

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

[0088] Fully developed zone: the wind farm interacts with the atmospheric boundary layer, so that the flow in this zone is fully developed. The range of the fully developed zone is P3 < x < P4, and the range of the fully developed zone is set as the distance between the end of the transition zone and the third last row of wind turbines.

[0089] Second transition zone: the transition zone between the fully developed zone and the group wake zone. The range of the second transition zone is P4 < x < P5, and the range of the second transition zone is set as the distance between three rows of wind turbines, i.e. the range between the third last row of wind turbines and the last row of wind turbines.

[0090] Group wake zone: the downstream region of the entire wind farm, and the flow in this zone is affected by the group wake of the entire wind farm. The range of the group wake zone is P5 < x < P6, and the specific range is related to the scale of the wind farm, the layout of the wind turbines, the atmospheric environment, and the characteristics of the wind turbines. The influence range of the group wake zone can reach 50 km, i.e. P6-P5 can reach 50 km.

[0091] Step S12: modeling and calculating the regions respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point includes the first row of wind turbine point wind speed, the inlet development zone wind turbine point wind speed, the first transition zone wind turbine point wind speed, the fully developed zone wind turbine point wind speed, and the second transition zone wind turbine point wind speed.

[0092] The different regions are respectively modeled and calculated to obtain the wind speed of each wind turbine point in the wind farm, including the first row of wind turbine point wind speed, the inlet development zone wind turbine point wind speed, the first transition zone wind turbine point wind speed, the fully developed zone wind turbine point wind speed and the second transition zone wind turbine point wind speed. Different methods are used to calculate the flow field of each region. When calculating the flow field of a certain region, there are many choices for the calculation model, such as Jensen model, Frandsen model, Gaussian model, etc. for the inlet development zone model. 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 is helpful to improve the accuracy of large-scale wind farm flow field and power generation calculation.

[0093] Further, the full-field blockage region is modeled and calculated;

[0094] The first row of wind turbine point wind speed is calculated according to the following formula:

[0095]

[0096] Wherein, U1 represents the first row of wind turbine point wind speed; U0 represents the free flow speed of the free flow area; △ x represents the average distance between adjacent wind turbines in the flow direction; △ y represents the average distance between adjacent wind turbines in the spanwise direction; D represents the average impeller diameter of the wind turbine; N rows represents the row number of wind turbines in the wind farm.

[0097] As Figure 2 shown, in the full-field blockage region, the range is P0

[0098]

[0099] Wherein, U1 represents the first row of wind turbine point wind speed; U0 represents the free flow speed of the free flow area; △ x represents the average distance between adjacent wind turbines in the flow direction; △ y represents the average distance between adjacent wind turbines in the spanwise direction; D represents the average impeller diameter of the wind turbine; N rows represents the row number of wind turbines in the wind farm.

[0100] The free flow velocity of a wind farm can be measured by arranging wind speed measuring devices, such as meteorological towers or wind measuring poles, in the wind farm. 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 flow velocity of the free flow region can be obtained.

[0101] The flow direction generally refers to the main direction of air or gas flowing through the fan group, and the average distance between adjacent fans in the flow direction refers to the average distance between fans arranged continuously in this direction.

[0102] The spanwise direction generally refers to the direction perpendicular to the flow direction, that is, the transverse width direction of the fan array. The average distance between adjacent fans in the spanwise direction refers to the average distance between fans arranged continuously in this direction.

[0103] Further, the modeling and calculation of the regions respectively include:

[0104] Modeling and calculating the inlet development region;

[0105] The wind speed at the fan point of the inlet development region is calculated according to the following formula:

[0106]

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

[0108] As Figure 2 shown in the inlet development region, the range is P1

[0109]

[0110] wherein, U 入口,iU1+∑j=1j=i-1(△u)j ij Cj,i T,j Cj w Dj,i j Dj

[0111] The velocity deficit is the reduction of wind speed due to the rotation of wind turbines and the 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 the design of the wind turbine, the shape of the blades, wind speed, and other factors. The wake diffusion coefficient is an empirical parameter that describes the rate at which the wake spreads in space. Different wind farm conditions, terrain, and wind turbine layouts can result in different wake diffusion coefficients.

[0112] Further, the modeling and calculating of the regions respectively comprises:

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

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

[0115]

[0116]

[0117] wherein U3 represents the wind speed at the wind turbine point in the fully developed region; U0 represents the free flow speed in the free flow region; H ABL H represents the height of the atmospheric boundary layer above the wind farm; z 0,lo z represents the surface roughness; z 0,hi z represents the equivalent roughness of the wind farm; z h D represents the average rotor diameter of the wind turbine; v w * v represents the equivalent eddy viscosity coefficient.

[0118] As Figure 2As shown, in the fully developed region, the range is P3 < x < P4, the velocity decreases from U3 to U4 along the flow direction, wherein U3 is the fan point wind speed of the fully developed region; U4 is the last fan point wind speed contained in the fully developed region. The flow field in the fully developed region reaches the fully developed state, so the flow field in this region is uniformly distributed, and the velocity of each fan point is a constant value (U3 = U4), which can be obtained by the equivalent roughness model, and is calculated according to the following formulas (5) and (6):

[0119]

[0120] Wherein, U3 represents the fan point wind speed of the fully developed region; U0 represents the free flow velocity of the free flow 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 rotor diameter of the fan; v w * represents the equivalent eddy viscosity coefficient.

[0121] The flow field in the fully developed region reaches the fully developed state, so the flow field in this region is uniformly distributed, and the velocity of each fan point is a constant value, so U3 can represent the fan point wind speed of the fully developed region. The atmospheric boundary layer height refers to the height of the near-surface air layer formed due to the friction of the ground, which is generally between 400-1000 meters, but the specific height will be affected by many factors such as topography, landform, climate, etc. The surface roughness refers to the size of the frictional resistance of the air flow generated by the ground surface, which is mainly determined by the size and distribution of the surface roughness elements. The equivalent roughness of the wind farm refers to the equivalent frictional resistance of the air flow generated by the entire wind farm 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 rotating center of the wind turbine blade (i.e. the hub) to the ground, and the standard height of the wind turbine hub is generally around 80-120 meters.

[0122] Further, the modeling and calculation of the regions respectively include:

[0123] modeling and calculating the first transition region;

[0124] calculating the fan point wind speed of the first transition region according to the following formula:

[0125]

[0126] Wherein, 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 speed decreases from U2 to U3 along the flow direction. Among them, 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 speed U of each fan position in this area 过渡1 can be calculated according to the following formula (7):

[0128]

[0129] Among them, 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 group wake zone;

[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 group wake zone:

[0133]

[0134] Among them, 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 group wake zone; U6 represents the outlet speed of the group wake zone; 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 group wake zone; x represents the flow direction distance.

[0135] 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 position corresponding to P5, usually taking the wind speed at the last row of fan positions included in the second transition zone. In the group wake zone, the range is P5 < x < P6, and the speed further increases from U5 to U6 along the flow direction. U6 is the outlet speed of the group wake zone. When U5 cannot be determined, the speed of the second transition zone and the group wake zone can be calculated according to the following formula (8):

[0136]

[0137] Among them, 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 group wake zone; U6 represents the outlet speed of the group wake zone; U4 represents the wind speed at the position 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; P6 represents the upper limit of the range of the group wake zone; x represents the flow distance.

[0138] The free incoming flow velocity of the wind farm can be obtained by arranging 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 wind speed meters, wind speed gauges), 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 of the regions includes:

[0140] Modeling and calculating the second transition zone;

[0141] According to the following formula, calculate the wind speed at the fan position in the second transition zone:

[0142]

[0143] Among them, U 过渡2 represents the wind speed at the fan position in the second transition zone; U5 represents the wind speed at the position corresponding to P5; P5 represents the upper limit of the range of the second transition zone; U4 represents the wind speed at the position 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.

[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 position corresponding to P5, and usually the wind speed at the position of the last row of fans included in the second transition zone is taken. When U5 can be determined by calculation, the speed U 过渡2 at each fan position 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 position in the second transition zone; U5 represents the wind speed at the position corresponding to P5; U4 represents the wind speed at the position 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; P5 represents the upper limit of the range of the second transition zone; x represents the flow distance.

[0147] Optionally, the modeling and calculating of the regions respectively comprises:

[0148] modeling and calculating the wake region of the group;

[0149] The fan point wind speed of the second transition region and the point wind speed of the wake region of the group are calculated according to the following formula:

[0150]

[0151] wherein, U 群尾流 represents the point wind speed of the wake region of the group; U5 represents the point wind speed 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 distance in the flow direction.

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

[0153]

[0154] wherein, U 群尾流 represents the point wind speed of the wake region of the group; U6 represents the outlet speed of the wake region of the group; U5 represents the point wind speed 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 distance in the flow direction.

[0155] Step S13: According to the wind turbine power curve and the wind speed of each fan point, the power generation power of each fan is obtained; the power generation power of each fan is multiplied by the corresponding time to obtain the total power generation of the wind farm.

[0156] Each wind turbine has a specific power curve that describes the power generation capability of the wind turbine at different wind speeds. This curve is usually provided by the wind turbine manufacturer and is based on the design parameters and test data of the wind turbine. For each fan, according to its current wind speed value, the corresponding power generation power on the wind turbine power curve is found. This usually involves interpolation calculation, because the actual wind speed may not exactly correspond to a certain point on the power curve. Then the corresponding time used to calculate the total power generation is determined, which can be hours, days, months or years, etc., and the choice of time depends on the specific needs of the analysis and the availability of data. Finally, the power generation power of each fan is multiplied by the corresponding time to obtain the total power generation of the wind farm. The power generation capacity of the wind farm is quantified, which provides the 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 has 320 5.5MW wind turbines installed, the wind turbines are regularly arranged in the wind farm, there are 16 rows along the main wind direction (flow direction) and 20 columns in the span direction, the average rotor diameter of the wind turbines is D = 158m, the hub height of the wind turbine is z h = 110m, the average distance between adjacent wind turbines in the flow direction is △ x = 1100m, the average distance between adjacent wind turbines in the span direction is △ y = 450m, the thrust coefficient of the wind turbine in the wind speed range of the case is C T = 0.8, the atmospheric boundary layer height above the wind farm is H ABL = 1000m, the ground roughness is z 0,lo = 0.001m, the equivalent roughness of the wind farm is z 0,hi = 2.5m, the free flow wind speed is U0 = 8m / s, and the wake diffusion coefficient is k w = 0.05.

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

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

[0161] Region Range (m) Fan (discharge) Free-stream region x<-1580 Full-field blockage region -1580<x<0 1 Inlet development region 0<x<5623 2～6 First transition region 5623<x<8800 7～9 Fully developed region 8800<x<13200 10～13 Second transition region 13200<x<16500 14～16 Wake region 16500<x<66500

[0162] According to step S12, the different regions are respectively modeled and calculated to obtain the wind speed of each wind turbine point in the wind farm, and the results are shown in Table 2.

[0163] Table 2. Wind speed results of each wind turbine point

[0164]

[0165]

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

[0167] Based on the same inventive concept, the embodiments of the present application disclose a large wind farm flow field sub-region modeling and calculating device, which comprises:

[0168] The partition module is used for dividing the large wind farm into different regions; the regions include: a free flow region, a full-field blocking region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a group wake region.

[0169] The first calculation module is configured to perform modeling calculation on the regions respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point includes the first-row wind turbine point wind speed, the inlet development zone wind turbine point wind speed, the first transition zone wind turbine point wind speed, the fully developed zone wind turbine point wind speed and the second transition zone wind turbine point wind speed.

[0170] The second calculation module is configured to obtain the power generation of each wind turbine according to the wind turbine power curve and the wind speed of each wind turbine point, and multiply the power generation of each wind turbine by the corresponding time to obtain the total power generation of the wind farm.

[0171] In an optional implementation, the partition module is specifically configured to:

[0172] The position of the first-row wind turbine of the wind farm is taken as the origin; the origin represents x = 0, corresponding to P1, wherein x represents the flow distance;

[0173] The range of the free flow zone is determined as x < P0; the range of the free flow zone is the entire upstream region of the wind farm;

[0174] The range of the full-field blockage zone is determined as P0 < x < P1; the range of the full-field blockage zone is 10D, that is, P1-P0 = 10D, wherein D represents the average impeller diameter of the wind turbine;

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

[0176] The range of the first transition zone is determined as P2 < x < P3; the range of the first transition zone is the distance between three rows of wind turbines;

[0177] The range of the fully developed zone is determined as P3 < x < P4; the range of the fully developed zone is the distance between the end of the first transition zone and the third-to-last row of wind turbines;

[0178] The range of the second transition zone is determined as P4 < x < P5; the range of the second transition zone is the distance between the third-to-last row of wind turbines and the last row of wind turbines;

[0179] The range of the group wake zone is determined as P5 < x < P6; the range of the group wake zone is the entire downstream region of the wind farm.

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

[0181] modeling calculation is performed on the full-field blockage region;

[0182] The first exhaust fan point wind speed is calculated according to the following formula:

[0183]

[0184] wherein U1 represents the first exhaust fan point wind speed; U0 represents the free flow speed of the free flow region; Δu x represents the average distance between adjacent fans in the flow direction; Δx 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 row number of the wind turbine in the wind farm.

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

[0186] modeling calculation is performed on the inlet development region;

[0187] The fan point wind speed of the inlet development region is calculated according to the following formula:

[0188]

[0189] wherein U 入口,i represents the wind speed of the i-th fan point in the inlet development region, i.e. the fan point wind speed of the inlet development region; U1 represents the first exhaust fan point wind speed; Δu ij represents the speed loss caused by the wake effect of the j-th wind turbine at the i-th fan point; 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 fan; Δy represents the spanwise distance between the j-th wind turbine and the i-th fan.

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

[0191] modeling calculation is performed on the fully developed region; the fan point wind speed of the fully developed region is constant;

[0192] The fan point wind speed of the fully developed region is calculated according to the following formula:

[0193]

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

[0195] In one alternative implementation, the first computing module is specifically used for:

[0196] Modeling and calculations are performed on the first transition region;

[0197] Calculate the wind speed at the fan location in the first transition zone using the following formula:

[0198]

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

[0200] In one alternative implementation, the first computing module is specifically used for:

[0201] Modeling and calculation are performed on the second transition region and the group wake region;

[0202] Calculate the wind speed at the fan location in the second transition zone and the wind speed at the fan wake zone location using the following formula:

[0203]

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

[0205] In one alternative implementation, the first computing module is specifically used for:

[0206] Modeling and calculations are performed on the second transition region;

[0207] Calculate the wind speed at the fan location in the second transition zone using the following formula:

[0208]

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

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

[0211] model the wake zone;

[0212] calculate the point speed of the second transition zone and the point speed of the wake zone according to the following formula:

[0213]

[0214] wherein U 群尾流 represents the point speed of the wake zone; U5 represents the point speed corresponding to P5; P5 represents the upper limit of the second transition zone; U6 represents the outlet speed of the wake zone; P6 represents the upper limit of the wake zone; and x represents the streamwise distance.

[0215] The large-scale wind farm flow field sub-region modeling and calculation method provided by the embodiments of the present disclosure divides a large-scale wind farm into different regions, including a free incoming flow region, a full-field blocked region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a wake region. The regions are modeled and calculated respectively to obtain the wind speed of each wind turbine point in the wind farm. The power generation power of each wind turbine is obtained according to the wind turbine power curve and the wind speed of each wind turbine point. The total power generation of the wind farm is obtained by multiplying the power generation power of each wind turbine by the corresponding time. The flow characteristics of different regions of the large-scale wind farm are modeled and calculated, which improves the accuracy of the flow field and power generation calculation of the large-scale wind farm and reduces the investment risk of the wind farm.

[0216] It should be noted that, for the method embodiments, in order to simply describe, they are all described as a series of action combinations, but those skilled in the art should know that the sequence of the actions described is not used to limit the present disclosure, because according to the embodiments of the present disclosure, certain steps can be performed in other sequences or at the same time. Secondly, those skilled in the art should know that the embodiments described in the specification all belong to preferred embodiments, and the actions involved are not necessarily the necessary actions of the present disclosure.

[0217] The embodiments of the present disclosure also provide an electronic device, which refers to Figure 4 , Figure 4 is a schematic diagram of the electronic device shown in the embodiments of the present disclosure. As Figure 4As shown, the electronic device 100 includes a memory 110 and a processor 120, the memory 110 and the processor 120 are in communication connection through a bus, and the memory 110 stores a computer program, the computer program can run on the processor 120, and then realize the steps in the large wind farm flow field sub-regional modeling calculation method disclosed by the embodiments of the disclosure.

[0218] The embodiments of the disclosure also provide a computer readable storage medium, when 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 wind farm flow field sub-regional modeling calculation method disclosed by the embodiments of the disclosure.

[0219] Each of the embodiments in the specification is described in a progressive manner, and each embodiment focuses on the difference from other embodiments. The same and similar parts between the embodiments can be referred to each other.

[0220] The embodiments of the present application are described with reference to the flowcharts and / or block diagrams of the methods, devices, electronic devices and computer program products according to the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing terminal equipment to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal equipment realize the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one block or multiple blocks.

[0221] These computer program instructions can also be stored in a computer readable storage medium which can guide the computer or other programmable data processing terminal equipment to work in a specific way, so that the instructions stored in the computer readable storage medium generate a product including instruction devices, which realize the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one block or multiple blocks.

[0222] These computer program instructions can also be loaded into a computer or other programmable data processing terminal equipment, so that a series of operation steps are performed on the computer or other programmable terminal equipment to generate a computer implemented process, so that the instructions executed on the computer or other programmable terminal equipment provide steps for realizing the functions specified in the flowcharts and / or block diagrams. Figure 1 The functions specified in one flow or multiple flows and / or blocks Figure 1 The functions specified in one block or multiple blocks.

[0223] Although preferred embodiments of the application have been described in detail, those skilled in the art will appreciate that various modifications and alterations can be made to the embodiments without departing from the scope of the application. Accordingly, the appended claims are intended to encompass all such modifications and alterations. In particular, the application is intended to cover any and all combinations of the features set forth in the claims, and any and all combinations and sub-combinations of the features set forth in the detailed description above.

[0224] Finally, it should be noted that the terms "first", "second", and the like, herein do not denote any order, quantity, combination, or arrangement, but are used to distinguish one element from another, and are not intended to mean or imply that a specific order, quantity, combination, or arrangement is preferred, required, or necessary. Also, the use of "including", "containing", or "comprising" and variations thereof, does not imply that the listed material or step is the only one of its type to be encompassed. The terms "including", "containing", and "comprising" are used herein to mean that other elements can also be present. The use of "including" or "comprising" does not exclude the presence of elements other than those listed in the reference. The use of "including" or "comprising" does not exclude the presence of additional elements of the same type as those listed in the reference. The use of "including" or "comprising" does not exclude the presence of additional elements of the same type as those listed in the reference.

[0225] The above describes in detail the large wind farm flow field sub-region modeling calculation method, device, equipment and storage medium provided by the application. The principles and implementation manners of the application are described by using specific examples. The above embodiment descriptions are only used to help understand the method of the application and its core idea. Meanwhile, for those skilled in the art, according to the idea of the application, the specific implementation manners and application ranges can be changed. In summary, the content of the specification should not be understood as a limitation of the application.

Claims

1. A large wind farm flow field sub-region modeling calculation method, characterized in that, The method comprises the following steps: dividing a large wind farm into different regions; the regions include: a free 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; modeling and calculating the regions respectively to obtain the wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point includes: the first row of wind turbine point wind speed, the inlet development region wind turbine point wind speed, the first transition region wind turbine point wind speed, the fully developed region wind turbine point wind speed and the second transition region wind turbine point wind speed; obtaining the power generation of each wind turbine according to the wind turbine power curve and the wind speed of each wind turbine point; multiplying the power generation of each wind turbine by the corresponding time to obtain the total power generation of the wind farm; the modeling and calculating of the regions respectively include: modeling and calculating the full field blockage region; calculating the first row of wind turbine point wind speed according to the following formula: ， Wherein, U1 represents the first exhaust fan point position wind speed; U0 represents the free flow speed of the free flow area; △ x represents the average distance between the adjacent fans in the flow direction; △ y represents the average distance between the adjacent fans in the spanwise direction; D represents the average impeller diameter of the fan; N rows represents the row number of the wind turbine in the wind farm; modeling and calculating the inlet development region; calculating the inlet development region wind turbine point wind speed according to the following formula: ， ， wherein, U 入口,i represents the wind speed at the i-th fan point in the inlet development zone, i.e., the wind speed at the fan point in the inlet development zone; U1 represents the wind speed at the first row of fan points; Δu ij represents the speed loss caused by the j-th fan's wake effect at the i-th fan point; C T,j represents the thrust coefficient of the j-th fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; Δx represents the streamwise distance between the j-th fan and the i-th fan; and Δy represents the spanwise distance between the j-th fan and the i-th fan. modeling and calculating the fully developed region; the fully developed region wind turbine point wind speed is constant; calculating the fully developed region wind turbine point wind speed according to the following formula: ， ， where U3 represents the wind speed at the point of the wind turbine in the fully developed zone; U0 represents the free-stream velocity in the free-stream zone; H ABL represents the height of the atmospheric boundary layer 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 rotor diameter of the wind turbine; represents the equivalent eddy viscosity coefficient; modeling and calculating the first transition region; calculating the first transition region wind turbine point wind speed according to the following formula: ， wherein U 过渡1 U1 represents the fan point wind speed of the first transition region; U3 represents the fan point wind speed of the fully developed region; U2 represents the point wind speed corresponding to P2; P2 represents the lower limit of the range of the first transition region; P3 represents the upper limit of the range of the first transition region; and x represents the streamwise distance.

2. The method according to claim 1, wherein, the modeling and calculating of the regions respectively include: modeling and calculating the second transition region and the group wake region; calculating the second transition region wind turbine point wind speed and the group wake region point wind speed according to the following formula: ， wherein U 过渡2 represents the point speed of the second transition zone; U 群尾流 represents the point speed of the group wake zone; U6 represents the exit velocity of the group wake zone; U4 represents the point speed 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; and x represents the streamwise distance.

3. The method of claim 1, wherein, the modeling and calculating of the regions respectively include: modeling and calculating the second transition region; calculating the second transition region wind turbine point wind speed according to the following formula: ， wherein U 过渡2 U5 represents the point speed corresponding to P5; P5 represents the upper limit of the second transition zone; U4 represents the point speed of the last row of fans included in the fully developed zone; P4 represents the lower limit of the second transition zone; and x represents the streamwise distance.

4. The method according to claim 3, wherein, the modeling and calculating of the regions respectively include: modeling and calculating the group wake region; calculating the second transition region wind turbine point wind speed and the group wake region point wind speed according to the following formula: ， wherein U 群尾流 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 flow zone; P6 represents the upper limit of the range of the group wake flow zone; and x represents the flow distance.

5. The method according to any one of claims 1-4, wherein, the dividing of the large 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, wherein x represents the flow distance; determining that the range of the free flow region is x determining that the range of the full field blockage region is P0 determining a range of the inlet development zone as P1 < x < P2; the range of the inlet development zone is wherein a is expressed as an empirical parameter; H ABL denotes the atmospheric boundary layer height above the wind farm; determining that the range of the first transition region is P2 determining that the range of the fully developed region is P3 determining that the range of the second transition region is P4 determining that the range of the group wake region is P5 6. A large wind farm flow field sub-region modeling and computing device, characterized in that, the device includes: The partition module is used for dividing a large wind farm into different regions; the regions include: a free incoming flow region, a full-field blocking region, an inlet development region, a first transition region, a fully developed region, a second transition region, and a group wake region. The first calculation module is used for modeling calculation on the regions respectively, to obtain a wind speed of each wind turbine point in the wind farm; the wind speed of the wind turbine point includes: a first-row wind turbine point wind speed, an inlet development region wind turbine point wind speed, a first transition region wind turbine point wind speed, a fully developed region wind turbine point wind speed, and a second transition region wind turbine point wind speed. The second calculation module is used for obtaining a power generation of each wind turbine according to a wind turbine power curve and the wind speed of each wind turbine point; multiplying the power generation of each wind turbine by a corresponding time to obtain a total power generation of the wind farm. The first calculation module is specifically used for: modeling calculation on the full-field blocking region; calculating the first-row wind turbine point wind speed according to the following formula: ， wherein U1 represents the first row of fan point wind speed; U0 represents the free flow speed of 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 row number of the wind turbine of the wind farm; modeling calculation on the inlet development region; calculating the inlet development region wind turbine point wind speed according to the following formula: ， ， wherein, U 入口,i represents the wind speed at the i-th fan point in the inlet development zone, i.e., the wind speed at the fan point in the inlet development zone; U1represents the first row of fan point wind speed; △u ij represents the speed loss caused by the j-th fan wake effect at the i-th fan point; C T,j represents the thrust coefficient of the j-th fan; k w represents the wake diffusion coefficient; D j represents the impeller diameter of the j-th fan; △x represents the streamwise distance between the j-th fan and the i-th fan; and △y represents the spanwise distance between the j-th fan and the i-th fan. modeling calculation on the fully developed region; the fully developed region wind turbine point wind speed is constant; calculating the fully developed region wind turbine point wind speed according to the following formula: ， ， where U3 represents the wind speed at the point of the wind turbine in the fully developed zone; U0 represents the free-stream velocity in the free-stream zone; 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 rotor diameter of the wind turbine; represents the equivalent eddy viscosity coefficient; modeling calculation on the first transition region; calculating the first transition region wind turbine point wind speed according to the following formula: ， wherein U 过渡1 U1 represents the fan point wind speed of the first transition region; U3 represents the fan point wind speed of the fully developed region; U2 represents the point wind speed corresponding to P2; P2 represents the lower limit of the range of the first transition region; P3 represents the upper limit of the range of the first transition region; and x represents the streamwise distance.

7. An electronic device, comprising: The device comprises: a processor, a memory, and a computer program stored on the memory and executable on the processor; when the processor executes the program, the large wind farm flow field partition modeling calculation method in any one of claims 1 to 5 is implemented.

8. 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 wind farm flow field partition modeling calculation method in any one of claims 1 to 5.

Images