Fan regular arrangement method and device considering starting phase and row of dislocations

Abstract

This application discloses a method and apparatus for wind turbine layout considering starting phase and row-column misalignment. The method includes the following steps: acquiring calculation data, acquiring the deployable area, initializing rotation angle, constructing and mapping a coordinate system, traversing parameter combinations, generating a candidate turbine location set, determining the layout scheme, reversing the rotation and mapping back to the global coordinate system, evaluating and checking the wake and outputting the results. Compared with existing technologies, the wind turbine regular layout method of this application introduces a starting point phase: by translating the array starting point within the basic unit, the system covers the "landing point difference" of the turbine positions caused by different starting points, significantly increasing the number of feasible solutions. By defining row and column offsets by g(j,ω), parallelogram arrays can be formed, and various regular array forms such as alternating odd and even offsets and periodic offsets can be uniformly described. At the same time, this application is adapted to complex boundary and multi-restricted area sites: by generating candidates and then performing constraint judgment and screening, it is applicable to complex geometric scenarios such as multi-boundary, concave polygons, and restricted areas with holes. Furthermore, it is configured with discrete search to adapt to different precisions: the range and step size of column spacing, row spacing, row and column offset parameters, and phase parameters can be configured by users or empirical rules, which is convenient for the needs and computing power control at different stages of the project. Moreover, it can be combined with the engineering wake model, based on the wind rose frequency or wind speed and direction time series weighting, to directly output the optimal solution for AEP / wake loss.

Description

Technical Field

[0001] This application relates to wind farm design, and in particular to a method and apparatus for regular wind turbine arrangement that takes into account starting phase and staggered rows and columns. Background Technology

[0002] Wind turbine layout is a critical aspect of wind farm design, directly impacting wake losses, annual power generation (AEP), project investment costs, and operation and maintenance efficiency. In engineering practice, regular row-and-column layouts, parallelogram / staggered layouts, and experience-based manual layouts are commonly used. However, existing layout methods typically have the following shortcomings: (1) Insufficient controllability and exhaustiveness: Regular arrays often require manual adjustment of the starting position and arrangement direction, making it difficult to systematically enumerate and cover a large number of feasible solutions under controllable computational load; (2) The starting point phase factor was not systematically incorporated: Under the same column spacing and row spacing, the translation of the array starting point in a basic unit (starting point phase) will significantly change the number and spatial distribution of positions falling within the boundary, especially when there is a concave polygon boundary or a restricted area; if the starting point is fixed, it is easy to miss feasible solutions and affect the evaluation results; (3) Insufficient expressive power of staggered mechanism: Common staggered mechanisms often use simple rules such as fixed offset or alternation of odd and even, which are difficult to uniformly describe different offset amplitudes, different periods and different staggered forms, and are difficult to meet the needs of engineering to explore and optimize more feasible solutions; (4) Insufficient representativeness of a single arrangement direction: Generating a regular array based solely on the dominant wind direction or a single direction may not cover candidate layouts across the entire circumference, thus missing out on better options with less wake loss or higher AEP.

[0003] Therefore, existing technologies urgently need a wind turbine layout optimization method that can systematically generate regular array candidates within a configurable parameter range, explicitly introduce starting phase and staggered mechanism, and combine wake model for weighted evaluation under multiple wind conditions. Summary of the Invention

[0004] The purpose of this application is to provide a wind turbine layout optimization method that can generate a regular array of candidates within a configurable parameter range, explicitly introduce a starting phase and stagger mechanism, and combine a wake model for weighted evaluation under multiple wind conditions.

[0005] Specifically, this application discloses a wind turbine arrangement method that considers the starting phase and row-column misalignment, characterized in that the arrangement method includes the following steps: S1 obtains calculation data Obtain the data required for calculating the wind turbine layout rules. This data includes: site boundaries, wind turbine diameter, minimum turbine spacing constraints, boundary setback constraints, restricted area constraints, unit quantity constraints, and wind condition and turbine model data. S2 obtains the deployable area. The site boundary is shrunk inward according to the setback constraint, and the restricted area is deducted to obtain the deployable area.

[0006] S3 initialization Initialize rotation angle make Set the angle step size and determine the wind direction acquisition method; S4 constructs and maps the coordinate system. At the current rotation angle A local coordinate system is constructed under the corresponding target wind direction, so that the local y-axis is consistent with the downwind direction corresponding to the target wind direction, and the deployable area / restricted area is mapped to the local coordinate system. After mapping to the local coordinate system, the local bounding rectangle of the deployable area is calculated, and reference anchor points are set for the starting positioning of rule traversal.

[0007] The target wind direction is used to determine the reference wind direction of the downwind axis in the local coordinate system. In the local coordinate system, the target wind direction corresponds to the downwind direction, i.e., the y-axis, and its vertical direction corresponds to the crosswind direction, i.e., the x-axis.

[0008] S5 traversal parameter combinations Rotation angle in step S4 The parameter combinations are iterated through to form multiple candidate parameter combinations, including: column spacing, row spacing, row-to-column offset parameters, and phase parameters; wherein: Column spacing represents the distance between adjacent columns in the crosswind direction (local x-axis), and row spacing represents the distance between adjacent rows in the downwind direction (local y-axis). The range of values ​​for row spacing and column spacing can be input by the user or preset based on engineering experience, and the range is not fixed. The phase parameter is expressed as: This is used to describe the translational position of the starting point of the regular array within the basic unit. To ensure the periodic equivalence of the phase, it can take any real number and be normalized to a value by modulo. That is, adopt and The equivalent phase is obtained. Therefore, the starting phase offset vector p is defined as follows: ; The row offset parameter is used to control the degree of offset between adjacent rows in the crosswind direction, and the horizontal offset of the j-th row is defined by the offset function: Where: j is the row index, function For j about The function, For row and column shifting parameters. To achieve periodic equivalence and equivalent deduplication of the shifted column pattern, It can take any real number and pass through Normalization to ; S6 generates a set of candidate aircraft positions. Using the determined reference anchor point, and combining the starting phase offset vector P with the row and column offsets, a set of candidate aircraft positions is generated according to the following recursive rule: in i , j Integer index, As a reference anchor point, P The starting phase offset vector, For row and column offset, The column spacing is... The line spacing is [the specified value]. S7 Determines the layout plan Constraint determination is performed on the candidate aircraft position set in step S6 to obtain a feasible aircraft position set that meets the constraints, and an arrangement scheme is formed based on the target number of aircraft units.

[0009] S8 reverse rotation mapped back to the global coordinate system The feasible machine position coordinates after constraint screening are mapped back to the global coordinate system through reverse rotation, resulting in a rotation angle with the original coordinate system. The actual layout plan.

[0010] S9 Wake Stream Evaluation A wake evaluation model is used to perform a weighted evaluation of the layout scheme under multiple wind conditions. The wake evaluation module is a replaceable engineering evaluation unit, configured to calculate the field-level power output based on the layout scheme, wind data, and turbine information, and further obtain the total power generation average efficiency (AEP), wake loss, or a combination thereof as evaluation indicators. The wind data can come from wind roses (discrete wind direction sectors and their frequencies, representative wind speeds) or wind speed / direction time series; the evaluation indicators can be maximizing AEP, minimizing wake loss, or a weighted combination of both. After the evaluation is completed, the evaluation values ​​of the current parameter set and the corresponding scheme are compared with the historical best records, and the best records are updated.

[0011] S10 checks and outputs. First, confirm the current rotation angle. Check if all parameter combinations have been traversed. If not, return to step S6 to continue evaluating the untraversed combinations. If the current rotation angle If all parameter combinations have been traversed, then determine the current rotation angle. Has it exceeded 360°? If not, then perform a rotation angle check. and restart the calculation from step S4; where The angle step size set in step 3 and when If any value exceeds 360 degrees, the calculation stops and the globally optimal solution and parameters recorded in the wake evaluation module are output.

[0012] In a preferred embodiment, the values ​​of the column spacing and row spacing are set with upper and lower limits according to the fan diameter.

[0013] In a preferred embodiment, the staggered function includes the following form: Error-free column array: ; Parallelogram array: ,Right now ; Alternating odd-even staggered array: Adjacent rows switch between two preset offset states, such as odd rows. Even-numbered behaviors ;and Periodic misalignment: Repeating a row sequence with a period of length T.

[0014] In a preferred embodiment, the constraint determination in step S7 includes: boundary inclusion determination: whether the candidate point falls within the arrangeable area, and can further determine whether the minimum distance from the candidate point to the boundary meets the setback requirement; Restricted Area Exclusion Criteria: Whether a candidate point falls within a restricted area or is less than a preset safe distance from the restricted area boundary; and Minimum fan spacing determination: Whether the distance between any two fans meets the minimum fan spacing constraint, such as not being less than a multiple of the fan diameter.

[0015] In a preferred embodiment, the upper and lower limits of the row and column indexes in step S5 are extended with four-way redundancy that can include negative indexes. That is, based on the coverage of the outer rectangle, a certain index margin is extended to the four sides, and then the pruning is performed by boundary determination in the subsequent process.

[0016] In a preferred embodiment, step S7 further introduces a unit number redundancy parameter k, relaxing the target unit number N from "equal to N" to "falls into" The range will be defined, and solutions falling within the range will be retained for subsequent evaluation.

[0017] In a preferred embodiment, the wake evaluation module can be implemented using a Jensen model, a Gaussian model, or a TurbOPark model; it can also be implemented using existing wake simulation software / libraries.

[0018] In a preferred embodiment, the reference anchor point is located at the lower left corner of the local bounding rectangle.

[0019] In a preferred embodiment, the rotation mapping base point in the mapping process of step S4 is predetermined in the global coordinate system and is independent of the rotation angle.

[0020] The second aspect of this application discloses an arrangement device applicable to the above-mentioned wind turbine rules. The device includes a processor and a memory. The memory stores a computer program. When the computer program is run on the processor, it causes the processor to execute the steps of the method to generate a wind turbine arrangement scheme, determine constraints and perform wake weighted evaluation, and output optimization results.

[0021] The main advantages of this invention are: (1) Introducing the starting phase: through By shifting the array starting point within the basic unit, the system covers the "landing point difference" of the machine position caused by different starting points, significantly increasing the number of feasible solutions.

[0022] (2) Unify and expand the staggered form: through Define row and column offsets to form parallelogram arrays, and uniformly describe various regular array forms such as alternating odd and even offsets and periodic offsets; (3) Adapting to complex boundary and multi-restricted area sites: After candidate generation, constraint judgment and screening are performed, which is suitable for complex geometric scenes such as multi-boundary, concave polygon, and restricted areas with holes; (4) Configure discrete search to adapt to different precisions: The range and step size of column spacing, row spacing, row staggered column parameters and phase parameters can be configured by the user or by experience rules, which is convenient for different stages of the project and control of computing power. (5) Engineering evaluation: It can be combined with the engineering wake model and directly output the optimal solution for AEP / wake loss based on the time series weighting of wind rose frequency or wind speed and direction.

[0023] The specification of this application contains numerous technical features distributed across various technical solutions. Listing all possible combinations of these technical features (i.e., technical solutions) would make the specification excessively lengthy. To avoid this problem, all technical features disclosed in the above-described invention, in the various embodiments and examples below, and in the accompanying drawings should be considered as already described in this specification, unless such a combination of technical features is technically infeasible. For example, one example discloses feature A+B+C, and another example discloses feature A+B+D+E. Features C and D are equivalent technical means that serve the same function, and technically only one needs to be used; they cannot be used simultaneously. Feature E can technically be combined with feature C. Therefore, the solution A+B+C+D should not be considered as already described due to technical infeasibility, while the solution A+B+C+E should be considered as already described. Attached Figure Description

[0024] Figure 1 This is a flowchart of a fan arrangement method according to an embodiment of the present invention; Figure 2 This refers to the site location and local coordinate system under the global coordinate system (X, Y) defined in an embodiment of the present invention; Figure 3 This is a schematic diagram of the circumscribed rectangle, anchor points, and rule traversal under a local coordinate system (mapped and aligned) defined in an embodiment of the present invention; Figure 4 This is a schematic diagram of the starting vector P according to an embodiment of the present invention; Figure 5 This is a schematic diagram of the arrangement when the staggered parameter is 0 according to an embodiment of the present invention.

[0025] Figure 6 This is a schematic diagram of a parallelogram array arrangement of the staggered function according to an embodiment of the present invention; Figure 7 This is a schematic diagram of the staggered function in one embodiment of the present invention when it is arranged in an odd-even staggered manner; Figure 8 This is a structural block diagram of a fan regular arrangement device according to an embodiment of the present invention. Detailed Implementation

[0026] Through meticulous and in-depth research, the inventors of this invention have developed for the first time a wind turbine regular arrangement method that considers the starting phase and row-column misalignment. Compared with existing technologies, the wind turbine regular arrangement method of this application introduces the starting phase: through... By shifting the array starting point within the basic unit, the system covers the "landing point difference" of the machine positions caused by different starting points, significantly increasing the number of feasible solutions, and through... By defining row and column offsets, parallelogram arrays can be formed, and various regular array forms such as alternating odd and even offsets and periodic offsets can be uniformly described. At the same time, this application is adapted to complex boundary and multi-restricted area sites: after candidate generation, constraint judgment and screening are performed, which is applicable to complex geometric scenarios such as multi-boundary, concave polygons, and restricted areas with holes. It also configures discrete search to adapt to different precisions: the range and step size of column spacing, row spacing, row and column offset parameters and phase parameters can be configured by users or empirical rules, which is convenient for different stages of engineering needs and computing power control. Furthermore, it can be combined with engineering wake models, based on wind rose frequency or wind speed and direction time series weighting, to directly output the optimal solution for AEP / wake loss.

[0027] The present invention will now be described in detail with reference to specific embodiments. These embodiments will help those skilled in the art to further understand the invention, but do not limit the invention in any way. It should be understood that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention. These modifications and improvements all fall within the scope of protection of the present invention.

[0028] Example 1 Specifically, the first embodiment of the present invention is, for example... Figure 1-8 As shown, this embodiment discloses a method and apparatus for regular wind turbine arrangement considering starting phase and row / column misalignment. The wind turbine regular arrangement method is as follows: Figure 1 As shown, it includes the following steps: S1 obtains calculation data Obtain the data required for calculating the wind turbine layout rules. This data includes: site boundaries, wind turbine diameter, minimum turbine spacing constraints, boundary setback constraints, restricted area constraints, unit quantity constraints, and wind condition and turbine model data. Wind conditions can be wind speed and direction frequency statistics (total duration must be provided, e.g., 8760 hours per year) or wind speed / direction time series. The turbine model data includes at least power curves and necessary unit parameters.

[0029] S2 obtains the deployable area. The site boundary is shrunk inward according to the setback constraint, and the restricted area is deducted to obtain the deployable area.

[0030] Optionally, in one embodiment, the setback distance of the wind farm boundary is typically not less than 0.5D, while also meeting the principle of saving sea resources.

[0031] S3 initialization Initialize rotation angle Set angle step The wind direction acquisition method is determined; the target wind direction is used to determine the reference wind direction along the downwind axis of the local coordinate system. In the local coordinate system, the target wind direction corresponds to the downwind direction (y-axis), and its vertical direction corresponds to the crosswind direction (x-axis). When traversal is used, the rotation angle is calculated according to... Generate and iterate from 0° to the maximum angle less than 360°.

[0032] S4 constructs and maps the coordinate system. like Figure 2 As shown, at the current rotation angle A local coordinate system is constructed under the corresponding target wind direction, aligning the local y-axis with the downwind direction corresponding to that target wind direction, and the deployable area / restricted area is mapped to the local coordinate system; after mapping to the local coordinate system, as follows: Figure 3The calculation shows the local bounding rectangle of the deployable area, and a reference anchor point is set for the starting positioning of the rule traversal. In this embodiment, the reference anchor point is set at the lower left corner of the bounding rectangle.

[0033] The target wind direction is used to determine the reference wind direction of the downwind axis in the local coordinate system. In the local coordinate system, the target wind direction corresponds to the downwind direction, i.e., the y-axis, and its vertical direction corresponds to the crosswind direction, i.e., the x-axis.

[0034] Optionally, in one embodiment, the rotation mapping base point O of the mapping is predetermined in the global coordinate system and is related to the rotation angle. Irrelevant; for example, O can be the center of the globally bounding rectangle of the area to be crammed, or it can be specified by the user. The area is then oriented around O by angle. Perform coordinate mapping to the local coordinate system. For ease of understanding, this mapping is equivalent to "rotating the coordinate system clockwise", or it can be understood as rotating the site in the opposite direction by the same angle in the plane, so that the local coordinate axes are aligned with the global axes.

[0035] S5 traversal parameter combinations Rotation angle in step S4 The parameter combinations are iterated through to form multiple candidate parameter combinations, including: column spacing, row spacing, row-to-column offset parameters, and phase parameters; wherein: Column spacing represents the distance between adjacent columns in the crosswind direction (local x-axis), and row spacing represents the distance between adjacent rows in the downwind direction (local y-axis). The range of values ​​for row spacing and column spacing can be input by the user or preset based on engineering experience, and the range is not fixed. The phase parameters are as follows Figure 4 As shown, it is represented as This is used to describe the translational position of the starting point of the regular array within the basic unit. To ensure the periodic equivalence of the phase, it can take any real number and be normalized to a value by modulo. That is, adopt and The equivalent phase is obtained. Therefore, the starting phase offset vector P is defined as follows: ; The row offset parameter is used to control the degree of offset between adjacent rows in the crosswind direction, and the horizontal offset of the j-th row is defined by the offset function: Where: j is the row index, function For j about The function, For row and column shifting parameters. To achieve periodic equivalence and equivalent deduplication of the shifted column pattern, It can take any real number and pass through Normalization to ; Optionally, in one embodiment, the staggered parameters include: like Figure 5 The error-free column array shown: ; Parallelogram array: ,Right now ; like Figure 7 The alternating odd-even staggered array shown: adjacent rows switch between two preset offset states, such as odd rows. Even-numbered behaviors ;and Periodic misalignment: Repeating a row sequence with a period of length T.

[0036] S6 generates a set of candidate aircraft positions. Using the determined reference anchor point, and combining the starting phase offset vector P with the row and column offsets, a set of candidate aircraft positions is generated according to the following recursive rule: in i , j Integer index, As a reference anchor point, P The starting phase offset vector, For row and column offset, The column spacing is... The line spacing is [value].

[0037] Alternatively, in one embodiment, as Figure 5 As shown, the upper and lower limits of the row and column indexes are extended with four-way redundancy that can include negative indexes. That is, based on the coverage of the outer rectangle, a certain index margin is extended in all directions, and then pruning is performed by boundary determination in the subsequent process.

[0038] S7 Determines the layout plan Constraint determination is performed on the candidate aircraft position set in step S6 to obtain a feasible aircraft position set that meets the constraints, and an arrangement scheme is formed based on the target number of aircraft units.

[0039] The constraint determination includes: boundary inclusion determination: whether the candidate point falls within the arrangeable area, and can further determine whether the minimum distance from the candidate point to the boundary meets the setback requirement; Restricted Area Exclusion Criteria: Whether a candidate point falls within a restricted area or is less than a preset safe distance from the restricted area boundary; and Minimum fan spacing determination: Whether the distance between any two fans meets the minimum fan spacing constraint, such as not being less than a multiple of the fan diameter.

[0040] S8 reverse rotation mapped back to the global coordinate system The feasible machine position coordinates after constraint screening are mapped back to the global coordinate system through reverse rotation, resulting in a rotation angle with the original coordinate system. The actual layout plan.

[0041] S9 Wake Stream Evaluation A wake evaluation model is used to perform a weighted evaluation of the layout scheme under multiple wind conditions. The wake evaluation module is a replaceable engineering evaluation unit, configured to calculate the field-level power output based on the layout scheme, wind data, and turbine information, and further obtain the total power generation average efficiency (AEP), wake loss, or a combination thereof as evaluation indicators. The wind data can come from wind roses (discrete wind direction sectors and their frequencies, representative wind speeds) or wind speed / direction time series; the evaluation indicators can be maximizing AEP, minimizing wake loss, or a weighted combination of both. After the evaluation is completed, the evaluation values ​​of the current parameter set and the corresponding scheme are compared with the historical best records, and the best records are updated.

[0042] The AEP described in this embodiment is based on wind rose weighting. The commonly used annual power generation AEP in engineering is as follows: in Let Q be the frequency of the Qth wind condition, and the sum of the frequencies of all wind conditions is 1. This represents the field-level average power under this wind condition; Or annual wake loss: in, The reference power for a certain wind condition without wake is given.

[0043] The "target wind direction" is used to determine the downwind axis direction in the local coordinate system: in the local coordinate system, the target wind direction corresponds to the y-axis direction, and its perpendicular direction is the local x-axis direction. In this embodiment, a rotation angle is introduced. This is used to align the y-axis of the local coordinate system with the downwind direction corresponding to the target wind direction. The rotation mapping base point O is predetermined in the global coordinate system and is aligned with... Regardless of the context, the center of the globally bounding rectangle of the arbitrarily sized area can be taken, or it can be specified by the user. The arbitrarily sized area and the restricted area are rotated clockwise around the O. Perform coordinate mapping to the local coordinate system S10 checks and outputs. First, confirm the current rotation angle. Check if all parameter combinations have been traversed. If not, return to step S6 to continue evaluating the untraversed combinations. If the current rotation angle If all parameter combinations have been traversed, then determine the current rotation angle. Has it exceeded 360°? If not, then perform a rotation angle check. and restart the calculation from step S4; where The angle step size set in step 3 and when If any value exceeds 360 degrees, the calculation stops and the globally optimal solution and parameters recorded in the wake evaluation module are output.

[0044] Optionally, in one embodiment, the values ​​of the column spacing and row spacing are set with upper and lower limits according to the fan diameter.

[0045] Optionally, in one embodiment, a unit number redundancy parameter k can be further introduced in step S7, relaxing the target unit number N from "equal to N" to "falling into". The range will be defined, and solutions falling within the range will be retained for subsequent evaluation.

[0046] Optionally, in one embodiment, the wake evaluation module can be implemented using a Jensen model, a Gaussian model, or a TurbOPark model; it can also be implemented using existing wake simulation software / libraries.

[0047] Optionally, in one embodiment, the reference anchor point is located at the lower left corner of the local circumscribed rectangle.

[0048] This embodiment also discloses a method such as Figure 8 The arrangement device shown is suitable for the above-mentioned wind turbine rules. The device includes a processor and a memory. The memory stores a computer program. When the computer program runs on the processor, it causes the processor to execute the steps of the method to generate wind turbine arrangement schemes, determine constraints, perform wake weighted evaluation, and output optimization results. Furthermore, it should be noted that in this patent application, relational terms such as "first" and "second" are used merely to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one" does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element. In this patent application, if it refers to performing an action according to a certain element, it means performing the action at least according to that element, including two cases: performing the action only according to that element, and performing the action according to that element and other elements. Expressions such as "multiple," "repeatedly," and "various" include two, two times, two kinds, and more than two, more than two times, and more than two kinds.

Claims

1. A method for regular wind turbine arrangement considering starting phase and row / column misalignment, characterized in that, The arrangement method includes the following steps: S1 obtains calculation data Obtain the data required for calculating the wind turbine layout rules. This data includes: site boundaries, wind turbine diameter, minimum turbine spacing constraints, boundary setback constraints, restricted area constraints, unit quantity constraints, and wind condition and turbine model data. S2 obtains the deployable area. The site boundary is shrunk inward according to the setback constraint, and the restricted area is deducted to obtain the deployable area. S3 Rotation Angle Initialization Initialize rotation angle make Set the angle step size and determine the wind direction acquisition method; S4 constructs and maps the coordinate system. At the current rotation angle A local coordinate system is constructed under the corresponding target wind direction, so that the local y-axis is consistent with the downwind direction corresponding to the target wind direction, and the deployable area / restricted area is mapped to the local coordinate system. After mapping to the local coordinate system, the local bounding rectangle of the deployable area is calculated, and reference anchor points are set for the starting positioning of rule traversal. The target wind direction is used to determine the reference wind direction of the downwind axis in the local coordinate system. In the local coordinate system, the target wind direction corresponds to the downwind direction, i.e., the y-axis, and its vertical direction corresponds to the crosswind direction, i.e., the x-axis. S5 traversal parameter combinations The rotation angle determined in step S4 The parameter combinations are iterated through to form multiple candidate parameter combinations, including: column spacing, row spacing, row-to-column offset parameters, and phase parameters; wherein: The column spacing represents the crosswind direction, i.e., the distance between adjacent columns on the local x-axis, and the row spacing represents the downwind direction, i.e., the distance between adjacent rows on the local y-axis. The range of values ​​for the row spacing and column spacing can be input by the user or preset based on engineering experience, and the range is not fixed. The phase parameter is expressed as: This is used to describe the translational position of the starting point of the regular array within the basic unit. To ensure the periodic equivalence of the phase, it can take any real number and be normalized to a value by modulo. That is, adopt and The equivalent phase is obtained. Therefore, the starting phase offset vector p is defined as follows: ； The row offset parameter is used to control the degree of offset between adjacent rows in the crosswind direction, and the horizontal offset of the j-th row is defined by the offset function: Where: j is the row index, function For j about The function, For row and column shifting parameters. To achieve periodic equivalence and equivalent deduplication of the shifted column pattern, It can take any real number and pass through Normalization to ; S6 generates a set of candidate aircraft positions. Using the lower left corner of the local bounding rectangle as the reference anchor point, and combining the starting phase offset vector P with the row and column offsets, a set of candidate aircraft positions is generated according to the following recursive rules: in i , j Integer index, As a reference anchor point, P The starting phase offset vector, For row and column offset, The column spacing is... The line spacing is the specified line spacing. S7 Determines the layout plan Constraint determination is performed on the candidate aircraft position set in step S6 to obtain a feasible aircraft position set that meets the constraints, and an arrangement scheme is formed based on the target number of aircraft units. S8 reverse rotation mapped back to the global coordinate system The feasible machine position coordinates after constraint screening are mapped back to the global coordinate system through reverse rotation, resulting in a rotation angle with the original coordinate system. The actual layout plan. S9 Wake Stream Evaluation The wake evaluation model is used to perform a weighted evaluation of the layout scheme under multiple wind conditions. The wake evaluation module is a replaceable engineering evaluation unit and is configured to calculate the field-level power output based on the layout scheme, wind condition data and turbine information, and further obtain the total power generation AEP, wake loss or their combination evaluation index. After the evaluation is completed, the evaluation values ​​of the current parameter group and the corresponding scheme are compared with the historical best record and the best record is updated. S10 checks and outputs. First, confirm the current rotation angle. Check if all parameter combinations have been traversed. If not, return to step S6 to continue evaluating the untraversed combinations. If the current rotation angle If all parameter combinations have been traversed, then determine the current rotation angle. Has the rotation angle exceeded 360°? If not, adjust the current rotation angle. and restart the calculation from step S4; where The angle step size set in step 3 and when If the angle exceeds 360 degrees, the calculation should also be stopped and the next judgment should be made. If the current rotation angle If 360 degrees have been reached, the globally optimal solution and parameters recorded in the wake evaluation module will be output.

2. The arrangement method according to claim 1, characterized in that, The range of values ​​for column spacing and row spacing is set with upper and lower limits according to the fan diameter.

3. The arrangement method according to claim 1, characterized in that, The row offset parameters are configured to control the degree of offset between adjacent rows in the crosswind direction, and the lateral offset of the corresponding row number is defined by the offset function; the offset function includes the following forms: no offset array, parallelogram array, alternating odd and even offset array and / or periodic offset.

4. The arrangement method according to claim 1, characterized in that, The constraint determination in step S7 includes: boundary inclusion determination, restricted area exclusion determination, and minimum machine distance determination.

5. The arrangement method according to claim 1, characterized in that, In step S5, the upper and lower limits of the row and column indexes are extended with four-way redundancy that can include negative indexes. That is, based on the coverage of the outer rectangle, a certain index margin is extended to the four sides, and then the pruning is performed by boundary determination in the subsequent process.

6. The arrangement method according to claim 1, characterized in that, In step S7, a unit number redundancy parameter k is further introduced, relaxing the target unit number N from "equal to N" to "falling into". The range will be defined, and solutions falling within the range will be retained for subsequent evaluation.

7. The arrangement method according to claim 1, characterized in that, The wake evaluation module can be implemented using the Jensen model, Gaussian model, or TurbOPark model; it can also be implemented by calling existing wake simulation software / libraries.

8. The arrangement method according to claim 1, characterized in that, The reference anchor point is set at the lower left corner of the local bounding rectangle.

9. The arrangement method according to claim 1, characterized in that, The rotation mapping base point in the mapping process of step S4 is predetermined in the global coordinate system and is independent of the rotation angle.

10. A fan arrangement device suitable for any one of claims 1-9, characterized in that, The device includes a processor and a memory. The memory stores a computer program. When the computer program is run on the processor, it causes the processor to execute the steps of the method to generate a wind turbine layout scheme, determine constraints and perform wake weighted evaluation, and output optimization results.

Images