Belief Reliability Modeling Method for Wind Farms Considering Two-Directional Rotor Equivalent Wind Speed

: Compared to conventional energy sources, wind power is a clean energy source with high intermittence and uncertainty. As a system that converts wind energy into electricity, wind farms inevitably face severe reliability issues. In this paper, based on reliability theory, a new reliability modeling method for wind farms is proposed. Firstly, a belief reliability model for wind farms is constructed. Then, a power generation model based on two-directional rotor equivalent wind speed is established to represent the wind farm performance in the belief reliability model. Finally, several numerical studies are conducted to verify the power generation model under different wind speeds and directions, to demonstrate the belief reliability model with different levels of uncertainty, and to compare the belief reliability considering two-directional rotor equivalent wind speed with other methods.


Introduction
The energy crisis and environmental pollution are major issues around the world, and accelerating the development of clean energy is the trend of sustainable development.As an important renewable energy, wind has been increasingly considered due to its clean, safe, and large reserves.Wind power generation has also received widespread attention from countries around the world.The Global Wind Energy Council (GWEC) released the 2023 Global Wind Energy Report [1], which states that in 2022, the world's newly installed wind power capacity was 77.6 GW, with a cumulative capacity of 906 GW, an increase of 9% compared to 2021.Meanwhile, the report predicts that the cumulative installed capacity of global wind power will increase to 1586 GW over the next five years (2023-2027) under existing policies, with a very optimistic outlook.
Reliability issues relating to wind farms (WFs) have also emerged.In contrast to traditional power generation methods under steady-state conditions, wind energy is characterized by fluctuation, variation and intermittency.This results in wind turbines operating under non-stationary loads in highly stochastic environments, which introduces significant uncertainty.Much research has begun to focus on the reliability of WFs.
Some studies [2][3][4] applied traditional reliability metrics to wind turbines (WTs), such as failure rate, Mean Time Between Failures (MTTF), Mean Time To Repair (MTTR), and availability to obtain WF reliability through analytical or simulation methods.The analytical methods [5,6] derive the analytical model of WF reliability through the probability expression of the state of WTs.The simulation methods [7] mainly adopt Monte Carlo simulation, which samples the wind speed and the states of WTs to statistically obtain the reliability results.However, the reliability of WFs is mainly evaluated by reliability indices directly related to the performance of power generation, such as Duration per Interruption (D) [8,9], Load Not Supplied per Interruption (LNSI) [8,9], Energy Not Supplied Interruption (ENSI) [8,9] and other deterministic indices, as well as expected energy Symmetry 2024, 16, 614 3 of 20 represent the wind speed at all positions within that height layer on the wind wheel rotation plane, without considering the variation of wind speed in the horizontal direction.However, under the influence of wake effects, the wind speed at different horizontal directions within the same height layer may be quite different.For example, there is a situation where the downwind wheel rotation plane is partially affected by the wake and partially unaffected.The influence of wake within the same height layer will decrease to zero as the distance from the upwind wake center increases, and the wind speed will correspondingly increase until no reduction occurs.Therefore, considering the variation of wind speed on the height and the width of the wind turbine, we propose a two-directional rotor equivalent wind speed (2D-REWS) model to calculate the wind generation for wind farm belief reliability modeling.
This paper is organized as follows.Section 2 provides some preliminaries on the belief reliability theory and the three-dimensional Jensen-Gaussian wake model.The wind farm belief reliability model and the power generation model based on 2D-REWS are introduced in Sections 3 and 4, respectively.In Section 5, two numerical examples are employed to demonstrate the proposed method.Section 6 concludes this paper.

Preliminaries
Some basic concepts of belief reliability theory and a three-dimensional Jensen-Gaussian (3DJG) wake model are introduced in this section, which are used to deduce the mathematical model of the proposed method in this paper.

Belief Reliability Theory
Belief reliability is a new reliability metric proposed by Kang [12], which builds a theoretical discourse of reliability expressed via a margin equation, degradation equation and metric equation based on some basic definitions and principles.Reliability is the ability of a product to perform a specified function within a specified time under specified conditions.To describe this meaning clearly, the theory of belief reliability introduces the concepts of the state variable and the feasibility domain.The state variable is a set of variables that describe the behavior of a product, and the feasibility domain is the space of values of the state variable that allow a product not to completely lose its function.If the state variable is in the feasible domain, it means that the product works.The state variables describe the behavior of the product, and the feasibility domain reflects the criterion of functional feasibility.The failure behavior of the product is actually caused by the failure of the functional behavior, which is the ultimate expression of the functional behavior, so the core state variable is the variable that can describe the functional behavior.Functionality reflects human needs for the product, which can be represented by one or more specific performance parameters, and whether the product can perform the specified function depends on how much performance margin is provided for the function.Therefore, the basic principles of reliability theory are as follows: Principle 1 Margin-based Reliable Principle: The performance margin determines how reliable the object is.
Principle 2 Eternal Degradation Principle: The object performance margin undergoes irreversible degradation along the time vector.
Principle 3 Uncertainty Principle: The object performance margin and its degradation process are uncertain.
The related definitions of belief reliability, performance margin, and performance threshold of requirements are as follows.
Definition 1. Belief reliability: Assume that the state variable of a system ξ is an uncertain random variable whose reliable domain is Ξ.The belief reliability of this system is defined as the chance that the state variable is in the reliable domain, that is Definition 2. Performance Threshold of Requirements: The requirements of the subject for the specified function of the object.

Definition 3. Performance Margin:
The performance margin represents the distance between the performance parameter and the performance threshold of requirements, which can be written as a margin equation M = m→ t (P, P th ), where P is the performance parameter vector, and P th is the performance threshold of the requirement vector.
Based on the reliability principle of margin and uncertainty, the performance margin determines how reliable the object is, and the object's performance margin is uncertain.Thus, when the state variable is referred to as the performance margin, the belief reliability function of the performance margin can be written as a metric equation: The performance parameters are derived from the physical laws of various fields, which can be written as an interdisciplinary equation: where X and Y are the vectors of internal and external parameters, respectively, and t is the time vector.

Three-Dimensional Jensen-Gaussian (3DJG) Wake Model
The 3DJG model [26] is a wake model considering the wind shear and turbulence intensity on wake recovery rate based on the Jensen model using a Gaussian distribution.Some important wind speed equations are as follows.
With the hub center of the wind turbine as the coordinate origin, the horizontal velocity distribution at arbitrary height z and at the radial position y at the horizontal level at the downwind position x in the wake region is where u ref is the incoming wind speed measured at a reference height z ref (m/s), z + z hub is the absolute height above ground (m), α is the wind shear index, σ y is decided by r y = 2.58σ y = k y x + r 0 (m), and k y represents the wake expansion coefficient in the horizontal plane.
ζ is the speed distribution in the vertical direction (m/s): where and where the axial induction factor is a = 1 − √ 1 − C T /2, C T is the thrust coefficient, the initial wake radius r 0 = R (1 − a)/(1 − 2a) (m), and the wake radius is r z (x) = 2.58σ z = k z x + r 0 (m).

Belief Reliability Model for Wind Farm
Reliability is defined as the ability to perform a specified function within a specified period of time under a specified condition.For wind farms, their function is to ensure that the power generation meets its demand and to avoid mismatch at the medium-to-long-term scale.Wind farm reliability refers to the ability of a wind farm to match its power generation with the demand.Based on belief reliability theory, we propose wind farm belief reliability, build the wind farm margin equation, and measure the uncertainty.

Description of Wind Farm
A wind farm, denoted by WF, is the system that arranges multiple wind turbines in an area with good wind energy resources to generate power in a certain layout.The research object of this paper focuses on the wind farm systems composed of multiple wind turbines, whose power generation is influenced by the layout, as well as external meteorological and geographical factors.The wind farm systems do not include the transmission electronic systems, the control subsystems, the network monitoring systems, and the energy storage devices.
The total number of wind turbine subsystems in the wind farm is N, and each wind turbine is a subsystem of WF, denoted by WT n , that satisfies WT n ⊂ WF,n = 1, 2, • • • N.

Wind Farm Layout
The layout is the spatial arrangement of wind turbines in a wind farm.For the layout of wind farm systems in this paper, there are two representation methods.One is the absolute position coordinate system, and it is used to input the layout of the wind farms, and the other is the relative position coordinate system, which is used for calculating the power of wind farms in different wind directions.The two-coordinate system sometimes requires mutual conversion; here, an example to illustrate the conversion process is shown.

Example 1. Coordinate systems of the wind farm and their conversion:
As shown inFigure 1, in this example, we adopt the right-hand coordinate system O − XYZ as the absolute position coordinate system of the wind farm, and the position coordinate of the wheel hub height of the wind turbine WT n is (x n , y n , z n ).The coordinates are in meters.
where the axial induction factor is , is the thrust coefficient, the initial wake radius (m), and the wake radius is (m).

Belief Reliability Model for Wind Farm
Reliability is defined as the ability to perform a specified function within a specified period of time under a specified condition.For wind farms, their function is to ensure that the power generation meets its demand and to avoid mismatch at the medium-to-longterm scale.Wind farm reliability refers to the ability of a wind farm to match its power generation with the demand.Based on belief reliability theory, we propose wind farm belief reliability, build the wind farm margin equation, and measure the uncertainty.

Description of Wind Farm
A wind farm, denoted by , is the system that arranges multiple wind turbines in an area with good wind energy resources to generate power in a certain layout.The research object of this paper focuses on the wind farm systems composed of multiple wind turbines, whose power generation is influenced by the layout, as well as external meteorological and geographical factors.The wind farm systems do not include the transmission electronic systems, the control subsystems, the network monitoring systems, and the energy storage devices.
The total number of wind turbine subsystems in the wind farm is , and each wind turbine is a subsystem of , denoted by , that satisfies .

Wind Farm Layout
The layout is the spatial arrangement of wind turbines in a wind farm.For the layout of wind farm systems in this paper, there are two representation methods.One is the absolute position coordinate system, and it is used to input the layout of the wind farms, and the other is the relative position coordinate system, which is used for calculating the power of wind farms in different wind directions.The two-coordinate system sometimes requires mutual conversion; here, an example to illustrate the conversion process is shown.

Example 1. Coordinate systems of the wind farm and their conversion:
As shown in Figure 1, in this example, we adopt the right-hand coordinate system as the absolute position coordinate system of the wind farm, and the position coordinate of the wheel hub height of the wind turbine is .The coordinates are in meters.Under the wind direction θ , a new right-handed coordinate system was established.The hub height center of WT q is taken as the coordinate origin, the downwind direction is designated as the positive direction of X q axis, and the Y q axis is designated as a rotation of 90 • counterclockwise perpendicular to the downwind direction.The coordinate system of the wind farm was called the relative position coordinate system marked as θ − q.
At this time, the angle between the X q axis and the X axis is θ q (counterclockwise is positive), then the coordinates of each WT n at the hub height under the absolute coordinate Symmetry 2024, 16, 614 6 of 20 system (x n , y n , z n ) are transformed to the coordinates x q n , y q n , z q n under the θ − q relative position coordinate system.The coordinates x q n , y q n , z q n are as follows: x q n = x n − x q cos θ q + y n − y q sin θ q , (9) y q n = y n − y q cos θ q − x n − x q sin θ q , (10)

Wind Farm Interference Relationship
The interference relationship of a wind farm refers to the influences formed by the wake effect between wind turbines within the wind farm area.To describe the interference relationship between wind turbines, three necessary definitions are given in this paper.Definition 4. Upwind Turbine: Let the relative position coordinate system established with any wind turbine numbered q.For any WT q n 1 and WT q n 2 (n 1 < n 2 ), if their coordinates satisfy x q n 1 ≤ x q n 2 , then the WT q n 1 is said to be the upwind turbine of WT q n 2 ; otherwise, the WT q n 2 is said to be the upwind turbine of WT q n 1 .
Definition 5.The Most Upwind Turbine: Let the relative position coordinate system established with any wind turbine numbered q, for all wind turbines WT q n , if (1) WT q s satisfies: ∀n, x q s ≤ x q n , then the WT q s is considered to be the most upstream turbines.(2) WT q s satisfies: s = min{S λ } and x q S λ ≤ x q n , when there are multiple turbines WT , the turbine WT s n with the smallest serial number is considered to be the most upstream turbine.Definition 6.The Bidirectional interference Relationship Matrix: The bidirectional interference relation matrix W of the wind farm is a matrix of N × N, and each element w s ij represents whether there is an interference relationship between wind turbine WT s i and WT s j in the coordinate system with the most upwind turbine WT s as the origin, as follows: where Γ s i represents the wind turbine set affected by wind turbine WT s i in the coordinate system with the most upwind turbine WT s as the origin.

Belief Reliability Analysis of Wind Farm
Different from traditional reliability modeling and analysis methods, the belief reliability analysis method uses a margin equation to describe the degree to which the functions and performance of wind farms meet the requirements, that is, the reliability degree, and the metric equation is used to describe the influence of uncertainty on the reliability degree of the system.The former is the part used to obtain the reliability certainty of wind farm systems, and the latter is the influence of uncertainty during the process of completing the function.Both jointly complete the analysis of wind farm reliability.Therefore, this section describes the belief reliability analysis method based on the margin equation and measurement equation of wind farms.

Wind Farm Margin Equation
According to the Margin-based Reliable Principle [12], the performance margin determines how reliable the object is.The normal status is that the margin of the system output performance is greater than zero.When the margin of performance is less than or equal to zero, the system is in the failure status.According to Definition 3, the system's performance margin is characterized by its performance parameters and the corresponding performance threshold.The system's performance, which includes aspects such as motor output torque, battery capacity, and gear ratio, is an external manifestation of its functions.In line with the definition of reliability, the system's performance parameters should be capable of describing the execution of the system's fixed functions.Performance thresholds, on the other hand, set limits on the performance parameters associated with these functions.For wind farms, the most important function is to convert wind energy into electricity, the corresponding performance parameter is the power generation, and the performance threshold is the demand for power generation.The performance margin of wind farm M WF s the distance between power generation and its requirement: where E WF is the power generation, and E th is the requirement of wind farm power generation.
In the context of a wind farm system, the normal operational state is characterized by power generation exceeding the demand for power.Obviously, the power generation belongs to the larger-the-better performance parameter, and the greater the power generation exceeds the power demand, the larger the performance margin.Therefore, the distance function can be used to describe the relatively simple Euclidean distance description that can describe the desired performance, so the performance margin of a wind farm can be described as follows: This equation can ensure that when the power generation is greater than the power generation demand value, the performance margin is greater than 0, and the greater the power generation, the greater the performance margin.

Wind Farm Uncertainty Measurement
In practice, wind farms face various uncertainties, which may render them unreliable despite having sufficient margins.Therefore, to enhance the reliability assessment of wind farms, it is imperative to first analyze and quantify these uncertainties.
External uncertainty is the main source of uncertainty in wind farms.The working environment of wind farms is complex and variable, thus there are fluctuations and changes in related parameters such as meteorological conditions and geographical environment.These factors can have an impact on the margin of wind farms.Among those external parameters Y, wind speed v, which serves as the input for wind farm power generation, is the main factor.Thus, in this paper, we consider the uncertainty of wind speed.
For a given wind distribution, wind speed s regarded as a multi- dimensional random variable and follows a multidimensional normal distribution: where

Wind Farm Metric Equation
Given that wind speed directly influences the power generation of wind farms, the uncertainty associated with wind speed propagates to the performance parameters of the wind farms, which, in turn, causes the performance margin to have a distribution.According to the Uncertainty Principle [12], the performance margin greater than 0 means that the system is reliable, and the reliability of the wind farm is expressed as the probability that the performance margin of the system is greater than 0. According to the Definition 1, the metric equation of wind farm belief reliability is where R B s belief reliability, and M WF s the wind farm performance margin considered uncertain.
The expression for reliability can be further obtained as follows: where Φ M WF (v) epresents the distribution function of the performance margin M WF , which is related to the distribution function of wind speed, and the specific relationship is determined by the functional relationship between wind speed and power generation, which is referred to in Section 4.
From Equation ( 16), it is evident that the metric equation of belief reliability can be used to evaluate a system's reliability by integrating the uncertainties inherent in the system into the system performance margins.This approach significantly deviates from traditional reliability methods, which necessitate the input of component failure rates or state transition probabilities.This unique feature endows the belief reliability assessment method proposed in this paper for wind farms with the capability to be applied throughout the entire life cycle of the system, eliminating the need to acquire statistical information, such as component failure rates.

Power Generation Model Based on Two-Directional Rotor Equivalent Wind Speed
Drawing upon belief reliability theory, the first step in modeling the reliability of a wind farm involves establishing interdisciplinary equations for the performance parameters.Over an extended time scale, power generation is typically used as the primary indicator of wind farm performance.To account for the variation of wind speed across both the height and width of the wind turbine, we propose a power generation model based on two-dimensional rotor equivalent wind speed (2D-REWS).

Two-Directional Rotor Equivalent Wind Speed (2D-REWS) Model
In recent years, with increased wind turbine diameter, the concept of rotor equivalent wind speed (REWS) considering wind shear has been proposed, and the IEC 61400-12-1 standard (2nd revised in 2017) [27] also calculates the equivalent wind speed of a wind turbine by measuring the wind speed at multiple heights within the wind turbine.This is carried out to evaluate the power of the wind turbine.The equivalent wind speed of a wind turbine is defined as follows: where i is the number of each layer divided by the height of the wind turbine, v i is the wind speed at different heights of the wind turbine (m/s), A i is the area of the wind turbine rotor at the corresponding height of v i (m 2 ), and A = πR 2 is the area of a wind turbine rotor (m 2 ).
Based on the equivalent wind speed of the wind turbine, this paper not only considers the variation of wind speed on the height of the wind turbine but also that on the width of the wind turbine.In combination with a three-dimensional wake model considering wind shear, a two-directional rotor equivalent wind speed (2D-REWS) and its calculation method for wind turbines affected by wake effects are proposed.The specific approach involves using the concept of double integration along the height (m) and width of the Symmetry 2024, 16, 614 9 of 20 wind turbine rotor (m), and the wind speed (m/s) at each position on the rotor of a wind turbine can be equivalent to the height of the hub: In practice, the above process can be discretized and divided into multiple grids along the height and width of the rotor of wind turbine, as shown in Figure 2. The coordinates of each grid's center point represent the wind speed within that grid, and then equivalent calculations can be performed.
Based on the equivalent wind speed of the wind turbine, this paper not only considers the variation of wind speed on the height of the wind turbine but also that on the width of the wind turbine.In combination with a three-dimensional wake model considering wind shear, a two-directional rotor equivalent wind speed (2D-REWS) and its calculation method for wind turbines affected by wake effects are proposed.The specific approach involves using the concept of double integration along the height (m) and width of the wind turbine rotor (m), and the wind speed (m/s) at each position on the rotor of a wind turbine can be equivalent to the height of the hub: (19) or (20) In practice, the above process can be discretized and divided into multiple grids along the height and width of the rotor of wind turbine, as shown in Figure 2. The coordinates of each grid's center point represent the wind speed within that grid, and then equivalent calculations can be performed.

Wind Farm Power Generation Affected by Multiple Wakes
In practice, wake effects have a significant impact on the output power of wind farms.According to research, the power generation loss caused by wake effects is generally 5% to 15%, sometimes even up to 30% to 40% [28,29].Furthermore, a wind farm usually comprises more than one wind turbine, with each downstream wind turbine potentially being affected by multiple upstream wind turbines.Therefore, it is essential to consider the influence of multiple wind turbines.Consequently, in accordance with the 2D-REWS model we proposed, we adopted a three-dimensional Jensen-Gaussian (3DJG) wake model and considered the impact of multiple wakes to construct a wind farm power generation model.

Wind Farm Power Generation Affected by Multiple Wakes
In practice, wake effects have a significant impact on the output power of wind farms.According to research, the power generation loss caused by wake effects is generally 5% to 15%, sometimes even up to 30% to 40% [28,29].Furthermore, a wind farm usually comprises more than one wind turbine, with each downstream wind turbine potentially being affected by multiple upstream wind turbines.Therefore, it is essential to consider the influence of multiple wind turbines.Consequently, in accordance with the 2D-REWS model we proposed, we adopted a three-dimensional Jensen-Gaussian (3DJG) wake model and considered the impact of multiple wakes to construct a wind farm power generation model.
Considering any position (x s , y s , z s ) in the wind farm, as shown in Figure 3, the yellow area represents the wake area of wind turbine WT j , and the cross-section at the downwind position x s is a circle with the radius r z j x s − x s j .When position (x s , y s , z s ) is on and inside the circle, it will be affected by wake, and the wind speed follows Equation (5).When position (x s , y s , z s ) is outside the circle, it will not be affected by wake, and the wind speed can be expressed as the logarithmic wind profile: where z s s is the vertical coordinates of WT s (m).
Symmetry 2024, 16, x FOR PEER REVIEW 10 of 21 Considering any position in the wind farm, as shown in Figure 3, the yellow area represents the wake area of wind turbine , and the cross-section at the downwind position is a circle with the radius .When position is on and inside the circle, it will be affected by wake, and the wind speed follows Equation ( 5).When position is outside the circle, it will not be affected by wake, and the wind speed can be expressed as the logarithmic wind profile: (21) where is the vertical coordinates of (m).Therefore, in the coordinate system with wind turbine as the origin, the wind speed at location that is only affected by can be expressed as a piecewise function: (22) To determine the wind speed of downwind turbines affected by a single upwind turbine, it is essential to construct an interference relationship matrix for the wind farms.This matrix will facilitate the determination of the interference relationship between wind turbines.The specific process is as follows: The positional relationship between the upwind turbine wake circle and the rotor circle of the downwind turbine is crucial when constructing the interference relationship matrix.This positional relationship can include five kinds: external separation, external tangent, intersection, internal tangent, and inclusive, as shown in Figure 4. Therefore, in the coordinate system with wind turbine WT s as the origin, the wind speed at location (x s , y s , z s ) that is only affected by WT j can be expressed as a piecewise function: To determine the wind speed of downwind turbines affected by a single upwind turbine, it is essential to construct an interference relationship matrix for the wind farms.This matrix will facilitate the determination of the interference relationship between wind turbines.The specific process is as follows: The positional relationship between the upwind turbine wake circle and the rotor circle of the downwind turbine is crucial when constructing the interference relationship matrix.This positional relationship can include five kinds: external separation, external tangent, intersection, internal tangent, and inclusive, as shown in Figure 4.The downwind turbine is influenced by the wake of the upwind turbine when the two circles are at the last four different positional relationships: external tangent, intersection, internal tangent, or inclusion.According to the positional relationship between the two circles, the set of wind turbines , which are affected by the upwind turbine , can be defined as follows: The downwind turbine is influenced by the wake of the upwind turbine when the two circles are at the last four different positional relationships: external tangent, intersection, internal tangent, or inclusion.According to the positional relationship between the two circles, the set of wind turbines Γ s j , which are affected by the upwind turbine WT j , can be defined as follows: where r z j x s i − x s j is the wake radius of upwind turbine WT j at the downwind position x s i , and R i is the wind turbine radius of downwind turbine WT i .Further, we adopt an energy balance (EB) method to calculate the wind speed under the interaction of multiple wind turbine wakes.In the coordinate system with wind turbine WT s as the origin, the wind speed of WT i affected by the wake of J wind turbines can be expressed as follows: where v 0 is the incoming wind speed (m/s), v s eq,j is the 2D-REWS of WT j (m/s), v s,j eq,i is the 2D-REWS of WT i only affected by WT j (m/s).
Then, the wind power curve of wind turbine WT n is as follows: where ρ is the air density (kg/m 3 ), C p,n is the power coefficient of WT n , v eq,n is the 2D-REWS of WT n affected by multiple wake (m/s), p r,n is the rated power of WT n (kW), v in,n is the cut-in wind speed of WT n (m/s), v r,n is the rated wind speed of WT n (m/s), and v out,n is the cut-out wind speed of WT n (m/s).
The power generation of a wind farm is the sum of the power generation of all wind turbines: where f v eq,n is the wind distribution of v eq,n , E WF is the power generation of a wind farm (kWh), and t is the power generation time (h).

Numerical Examples
Two sets of numerical examples are presented in this section to better illustrate the belief reliability modeling method we proposed.The first set of example shows the changes in wind farm power generation under different wind speeds and directions, and the second provides belief reliability results under conditions of varying degrees of wind speed uncertainty.All the results presented in this section were calculated in MATLAB R2021b.

Power of Wind Farm under Different Wind Speeds and Directions
To verify the power generation model of wind farms based on the 2D-REWS model, we performed calculations and analyses of wind farm power under varying wind speeds and directions.
In the numerical examples used in this section, the wind farm is composed of nine wind turbines, and the layout is given as shown in Figure 5.The parameters are shown in Table 1.
To verify the power generation model of wind farms based on the 2D-REWS model, we performed calculations and analyses of wind farm power under varying wind speeds and directions.
In the numerical examples used in this section, the wind farm is composed of nine wind turbines, and the layout is given as shown in Figure 5.The parameters are shown in Table 1.In addition, all wind turbines in the wind farm are of the same type, and the related parameters are shown in Table 2.In addition, all wind turbines in the wind farm are of the same type, and the related parameters are shown in Table 2.We consider the power of a wind farm under a single wind direction and different wind speeds.The wind direction is 0 • (along the positive x-axis), and the wind speed ranges from 1 m/s to 32 m/s with intervals of 0.5 m/s.Other environmental parameters include the following: the wind speed reference height z ref = 65 m, the air density ρ = 1.225 kg/m 3 , the turbulence intensity k y = k z = 0.1 and the wind shear index α = 0.15.The power of the wind farm is shown in Figure 6.The following was found: 1.
The power of wind farms under different wind speeds basically conforms to their variation law with wind speed; increases accordingly.When the wind speed reaches a of around 13.5 m/s, the power of the wind farm remains stable.When the wind speed continues to increase to a of 25 m/s or above, the power generated by the wind farm becomes 0 again.This is because the wind speed has reached the wind turbine cut-out wind speed, and in order to protect the equipment from damage, the wind turbines will stop generating power.2. The power output of wind farms exhibits noticeable differences when the influence of wake is considered.This is especially pronounced when the wind speed approaches the rated wind speed and when it exceeds the cut-out wind speed.
The reason for this difference can be as follows: When the wind speed increases to , the upwind reaches , while the downwind affected by the wake has not yet reached , only the upwind starts generating electricity, and other downwind have not yet started.As the wind When the wind speed is lower than a v in of 3.5 m/s, the power of the wind farm is 0. When the wind speed gradually increases from 3.5 m/s, the power of the wind farm increases accordingly.When the wind speed reaches a v r of around 13.5 m/s, the power of the wind farm remains stable.When the wind speed continues to increase to a v out of 25 m/s or above, the power generated by the wind farm becomes 0 again.This is because the wind speed has reached the wind turbine cut-out wind speed, and in order to protect the equipment from damage, the wind turbines will stop generating power.

2.
The power output of wind farms exhibits noticeable differences when the influence of wake is considered.This is especially pronounced when the wind speed approaches the rated wind speed and when it exceeds the cut-out wind speed.
The reason for this difference can be as follows: When the wind speed increases to v in , the upwind WT reaches v in , while the downwind WT affected by the wake has not yet reached v in , only the upwind WT starts generating electricity, and other downwind WT s have not yet started.As the wind speed continues to increase, the amount of WT s reaching v in increases accordingly, increasing the power of the wind farm.However, the power of downwind WT s is still lower than that without considering the wake effect.This is due to the fact that the wind speed of the downstream wind turbines is lower than the ambient wind speed due to the wake effect, resulting in lower power generation than without the wake effect.Therefore, when considering the wake effect, the power of the wind farm is always lower than when not considering the wake effect.
When the wind speed reaches v r , the upwind WT reaches v r , while the downwind WT s affected by the wake has not yet reached v r .As the wind speed continues to increase, all the downwind WT s reach v r , there is no difference in the power generated by the wind farm with or without considering wake effect.As the wind speed continues to increase, the downwind WT s under the influence of the wake can still reach above v r , and those downwind WT s can also generate at its rated power.
Similarly, when the wind speed reaches v out , the upwind WT reaches v out and stops running, while the downwind WT 7 and WT 9 affected by the wake have not yet reached v out and continued generating electricity.With the gradual increase in ambient wind speed, the wind speed of downwind WT 7 reaches the cut-out wind speed preferentially, while the wind speed of downwind WT 9 still does not reach the cut-out wind speed due to more influence from the wake effect.As the wind speed continues to increase, all the downwind WT s reach v r , there is no difference in the power generated by the wind farm with or without considering wake effect.As the wind speed continues to increase, the downwind WT s under the influence of the wake can still reach above v r , and those downwind WT s can also generate at its rated power.

Under Different Wind Directions
We consider the power of a wind farm under a single wind speed and different wind directions.The wind speed is 12 m/s, and the wind direction is in the range of 0 • to 360 • , with intervals of 45 • .The other environmental parameters are the same as those shown in 5.1.1.The power of the wind farm is shown in Figure 7.The following can be found:

Belief Reliability of Wind Farm with Wind Speed Uncertainty
To present the belief reliability results based on the modeling method we proposed, a conventional wind distribution is utilized, followed by a comparison of the outcomes under wind distributions with varying degrees of dispersion.

Wind Farm Belief Reliability under Classical Wind Distribution
We consider the annual power generation of a wind farm under multiple wind directions and multiple wind speeds, which is a classical wind distribution used in [30][31][32][33][34].The wind speeds are 8 m/s, 12 m/s and 17 m/s, the wind direction is in the range of 0° to 360° with intervals of 10°, and the frequency of each wind speed and direction are shown in There are significant differences in the power of the wind farm in different wind directions.Under the eight wind directions, the power of the wind farm in descending order is 90 The wake influence between wind turbines varies in different wind directions, as shown in Figure 8. Wind farms with similar wake influences have similar power.

Belief Reliability of Wind Farm with Wind Speed Uncertainty
To present the belief reliability results based on the modeling method we proposed, a conventional wind distribution is utilized, followed by a comparison of the outcomes under wind distributions with varying degrees of dispersion.At 90 • and 270 • , the wake influences of the wind farm are symmetrical, and the power generated by the wind farm is equal.
At 45 • and 225 • , WT 3 and WT 7 are symmetrical, WT 1 , WT 2 and WT 4 are approximately symmetrical because their positions are approximately in a straight line and the distance between WT 1 and WT 2 is approximately equal to the distance between WT 2 and WT 4 ; thus, the power of the wind farm is similar under the two wind directions 45 • and 225 • .
At 0 • and 180 • , although the influence relationships of WT s are symmetrical, the intensity of wake influence is different; thus, the power of the wind farm is different under the two wind directions.

Belief Reliability of Wind Farm with Wind Speed Uncertainty
To present the belief reliability results based on the modeling method we proposed, a conventional wind distribution is utilized, followed by a comparison of the outcomes under wind distributions with varying degrees of dispersion.

Wind Farm Belief Reliability under Classical Wind Distribution
We consider the annual power generation of a wind farm under multiple wind directions and multiple wind speeds, which is a classical wind distribution used in [30][31][32][33][34].The wind speeds are 8 m/s, 12 m/s and 17 m/s, the wind direction is in the range of 0 • to 360 • with intervals of 10 • , and the frequency of each wind speed and direction are shown in Figure 9.The other environmental parameters are the same as those shown in Section 5.1.1.For the uncertainty of wind speed, to simplify the calculation and reduce the sampling dimension while keeping the sum of the three wind speeds unchanged, we considered that the first two wind speeds follow a two-dimensional normal distribution: (27) where , .
Here, we conducted 200 simulation samples.In each sample, the sum of the three wind speeds remained 37 and the frequency remained unchanged.The requirement of wind farm power generation is as follows: = 1.0 × 10 11 kWh.The belief reliability of the wind farm is 0.69.

Reliability Comparison with Different Dispersion Degrees of Wind Distribution
To compare the wind farm belief reliability under different dispersion degrees of wind distribution, we calculated the results under four other dispersion degrees of wind following different two-dimensional normal distribution.Still, the expected value in the two-dimensional normal distribution is , and the covariance matrices are For the uncertainty of wind speed, to simplify the calculation and reduce the sampling dimension while keeping the sum of the three wind speeds unchanged, we considered that the first two wind speeds follow a two-dimensional normal distribution: where µ = (8, 12), Σ = 0.3 0 0 0.3 .
Here, we conducted 200 simulation samples.In each sample, the sum of the three wind speeds remained 37 and the frequency remained unchanged.The requirement of wind farm power generation is as follows: E th = 1.0 × 10 11 kWh.The belief reliability of the wind farm is 0.69.

Reliability Comparison with Different Dispersion Degrees of Wind Distribution
To compare the wind farm belief reliability under different dispersion degrees of wind distribution, we calculated the results under four other dispersion degrees of wind following different two-dimensional normal distribution.Still, the expected value in the two-dimensional normal distribution is µ = (8, 12), and the covariance matrices are  2.The belief reliability of a wind farm is positively correlated with the uncertainty of wind speed.
With the increased variance of wind speed, the belief reliability of wind farms decreases.This indicates that the greater the dispersion of wind speed, the greater the uncertainty of the performance margin, the lower belief reliability.
These trends regarding the relationship between wind speed uncertainty and both performance margin and belief reliability can reflect the actual patterns.2. The belief reliability of a wind farm is positively correlated with the uncertainty of wind speed.

Comparison of Belief Reliability Considering HHWS, REWS and 2D-REWS
With the increased variance of wind speed, the belief reliability of wind farms decreases.This indicates that the greater the dispersion of wind speed, the greater the uncertainty of the performance margin, the lower belief reliability.
These trends regarding the relationship between wind speed uncertainty and both performance margin and belief reliability can reflect the actual patterns.

Comparison of Belief Reliability Considering HHWS, REWS and 2D-REWS
First, to provide a more intuitive explanation of the difference between HHWS, REWS and 2D-REWS, we illustrated three representative situations.Figures 4 and 12 show As the variance of wind speed increases, both the dispersion of the wind speed and the performance margin also increase.Under a relatively dispersed wind distribution, although there are higher performance margins, there are also more performance margins less than 0, and this indicates that the uncertainty of wind farm performance margin is also greater.

2.
The belief reliability of a wind farm is positively correlated with the uncertainty of wind speed.
With the increased variance of wind speed, the belief reliability of wind farms decreases.This indicates that the greater the dispersion of wind speed, the greater the uncertainty of the performance margin, the lower belief reliability.
These trends regarding the relationship between wind speed uncertainty and both performance margin and belief reliability can reflect the actual patterns.

Comparison of Belief Reliability Considering HHWS, REWS and 2D-REWS
First, to provide a more intuitive explanation of the difference between HHWS, REWS and 2D-REWS, we illustrated three representative situations.Figures 4 and 12 show the three-positional relationship between the upwind turbine wake circle (the yellow one, abbreviated as C up in the following) and the rotor circle of the turbine (the transparent one, abbreviated as C down in the following).According to HHWS, the wind speed of the center of Cdown (denoted by in (a) of Situation 1/2/3 in Figure 12) represents the speed of the whole downwind turbine.When is outside Cup, the entire wind turbine will be considered unaffected by the wake, such as in the case of Situation 1 (a), where the wind speed is overestimated; when is inside Cup, the entire wind turbine will be considered affected by the wake, such as in the case of Situation 2 (a), where the wind speed is underestimated; and when the centers of Cup and Cdown overlap, the wake effect far from the center is lower than that near the center, such as in the case of Situation 3 (a), where the wind speed is underestimated.
According to REWS, the wind speed at different heights on the line perpendicular to the center of Cdown (denoted by the blue line in (b) of Situation 1/2/3 in Figure 12) represents the speed at different heights of the downwind turbine, and the equivalent wind speed is calculated using Equation (18).When the blue line is totally outside Cup, the entire wind turbine will be considered unaffected by the wake, such as in the case of Situation 1 (b), where the wind speed is overestimated; when the blue line is totally inside Cup, the entire wind turbine will be considered affected by the wake, such as in the case of Situation 2 According to HHWS, the wind speed of the center of C down (denoted by O i in (a) of Situation 1/2/3 in Figure 12) represents the speed of the whole downwind turbine.When O i is outside C up , the entire wind turbine will be considered unaffected by the wake, such as in the case of Situation 1 (a), where the wind speed is overestimated; when O i is inside C up , the entire wind turbine will be considered affected by the wake, such as in the case of Situation 2 (a), where the wind speed is underestimated; and when the centers of C up and C down overlap, the wake effect far from the center is lower than that near the center, such as in the case of Situation 3 (a), where the wind speed is underestimated.
According to REWS, the wind speed at different heights on the line perpendicular to the center of C down (denoted by the blue line in (b) of Situation 1/2/3 in Figure 12) represents the speed at different heights of the downwind turbine, and the equivalent wind speed is calculated using Equation (18).When the blue line is totally outside C up , the entire wind turbine will be considered unaffected by the wake, such as in the case of Situation 1 (b), where the wind speed is overestimated; when the blue line is totally inside C up , the entire wind turbine will be considered affected by the wake, such as in the case of Situation 2 (b), where the wind speed is underestimated; when the centers of C up and C down overlap, the wake effect far from the blue line is lower than that near the blue line, such as in the case of Situation 3 (b), where the wind speed is underestimated.
According to discretization 2D-REWS method, the wind speed of every intersection point of blue and green lines in C down are used to calculate the equivalent wind speed, and there is no situation where the wind speed is underestimated or overestimated based on the continuous 2D-REWS method (Equations ( 19) and (20)).
Then, to analyze the wind farm belief reliability when using different methods for representing the speed of the wind turbine, we calculated the belief reliability results considering HHWS, REWS and 2D-REWS of three different wind farm layouts.
We consider the reliability of the wind farm under a single wind speed and single wind direction.The wind speed is 12 m/s, the wind direction is 0 • .The relation between WTs is shown in the first subfigure in Figure 8.The rotor circles and wake circles of WT 2 , WT 7 , and WT 9 are all concentric circles.
Wind farm B is a control group whose layout is the same as that shown in Section 5.2.1.The layouts of wind farms A and C are changed.To ensure that the qualitative relationship between the upwind and downwind turbines remains unchanged and only the quantitative influence between WTs changes, the position of WT 9 is moved.The position of WT 9 in wind farm A is (2200,3400,0), and the position of WT 9 in wind farm C is (2200,3600,0).
The requirement of wind farm power generation is E th = 7.72 × 10 10 kWh.The other parameters are the same as those shown in Section 5.2.1.Here, we also conducted 200 simulation samples.The belief reliability of the wind farm is shown in Table 3.The following can be found.Wind farm layouts A and C are similar to Situation 1, shown in Figure 8; thus, the wind speeds are overestimated.Wind farm layout B is similar to Situation 3, shown in Figure 8; thus, the wind speeds is underestimated.The values of belief reliability considering HHWS and REWS of wind farm layout A and C are both higher than that considering 2D-REWS, and the values of belief reliability considering HHWS and REWS of wind farm layout B are both lower than that considering 2D-REWS.The situation where belief reliability is underestimated and overestimated is the same as wind speed.In practice, more deviations will inevitably be introduced in reliability results when using HHWS or REWS.Thus, the belief reliability model considering 2D-REWS is more reasonable.

Conclusions
This paper proposes a new belief reliability modeling method for wind farms considering two-directional rotor equivalent wind speeds.Several conclusions can be drawn from this paper: 1.
Only based on the performance margin and its threshold of wind farms, the proposed belief reliability modeling method allows for reliability analysis without relying on wind turbine failure information.

2.
A 2D-REWS model was introduced, which incorporates a power generation model to enable belief reliability analysis that accounts for variations in wind speed across different heights and widths of wind turbines.

3.
Numerical examples demonstrate that the variation of wind farm power with wind speed and direction is consistent with the actual patterns, the uncertainties in performance margin and belief reliability of a wind farm are both positively correlated with uncertainties in wind speed, and the belief reliability model considering 2D-REWS is more reasonable.
Future research directions may include incorporating factors such as the degradation of individual wind turbines and additional sources of uncertainty into the belief reliability model for wind farms.Additionally, we are interested in exploring layout optimization models based on the concept of belief reliability specifically tailored to enhance overall performance.

Figure 1 .
Figure 1.Coordinate systems of the wind farm layout.Figure 1. Coordinate systems of the wind farm layout.

Figure 1 .
Figure 1.Coordinate systems of the wind farm layout.Figure 1. Coordinate systems of the wind farm layout.

Figure 3 .
Figure 3. Diagram of wake effect in wind farm.

Figure 3 .
Figure 3. Diagram of wake effect in wind farm.

Symmetry 2024 , 21 Figure 4 .
Figure 4. Relation between the upwind turbine wake circle and the rotor circle of downwind turbine.

Figure 4 .
Figure 4. Relation between the upwind turbine wake circle and the rotor circle of downwind turbine.

Figure 6 .
Figure 6.The power of a wind farm under different wind speeds.

Figure 6 .
Figure 6.The power of a wind farm under different wind speeds.

Symmetry 2024 , 21 Figure 7 .
Figure 7.The power of the wind farm under different wind directions.

Figure 8 .
Figure 8.The wake influence between wind turbines under different wind directions.

Symmetry 2024 , 21 Figure 7 .
Figure 7.The power of the wind farm under different wind directions.

Figure 8 .
Figure 8.The wake influence between wind turbines under different wind directions.

igure 8 .
The wake influence between wind turbines under different wind directions.

21 Figure 10 .
Figures 10 and 11 present the frequency histogram of the performance margin and the belief reliability results, and the following can be found: 1.The uncertainty of the wind farm performance margin is positively correlated with the uncertainty of wind speed.Symmetry 2024, 16, x FOR PEER REVIEW 17 of 21

Figure 11 .
Figure 11.Belief reliability under different dispersion degrees of wind distribution.

Figure 11 .
Figure 11.Belief reliability under different dispersion degrees of wind distribution.

Figure 11 .
Figure 11.Belief reliability under different dispersion degrees of wind distribution.

Symmetry 2024 , 21 Figure 12 .
Figure 12.Diagram of different methods for representing the speed of the wind turbine.

3 Figure 12 .
Figure 12.Diagram of different methods for representing the speed of the wind turbine.

Author Contributions:
Conceptualization, S.L. and R.K.; methodology, S.L.; software, S.L.; validation, S.L., M.W. and T.Z.; formal analysis, S.L., M.W. and T.Z.; investigation, S.L.; resources, M.W. and R.K.; data curation, S.L.; writing-original draft preparation, S.L.; writing-review and editing, S.L. and T.Z.; visualization, S.L.; supervision, M.W. and R.K.; project administration, M.W. and R.K.; funding acquisition, R.K. and M.W.All authors have read and agreed to the published version of the manuscript.Funding: This research was funded by the National Natural Science Foundation of China, grant number 62073009 and the Stable Supporting Project of Science & Technology on Reliability & Environmental Engineering Laboratory, grant number WDZC20220102.

Table 1 .
Parameters of wind farm layout.

Table 1 .
Parameters of wind farm layout.
Wind Turbine Order Number nx n (m) y n (m) z n (m)

Table 2 .
Parameters of wind turbines.