Evaluation of the Behavior of a 3 MW Wind Warm with a Permanent Magnet Synchronous Generator under the Wind Conditions of Benin ()
1. Introduction
From 2000 to 2020, the global wind energy market grew from 17.4 GW to 743 GW, with an average annual growth rate of about 12.6%, driven by the development of variable-speed technologies. With the increasing demand for electricity, projections indicate that wind energy will reach 29.1% in 2030 and 34.2% in 2050 [1] [2]. Benin, characterized by an electricity deficit and dependence on increasingly scarce fossil energy sources, often faces instability issues within its power grid. To address this energy challenge, Benin relies on its photovoltaic potential and aligns its energy policy with the renewable energy production projects promoted by ECOWAS for the 2019-2033 period. To maintain the balance between supply and demand, the construction and connection of new production units using inertial renewable sources such as wind turbines should be considered. Increasing production capacity contributes to the stability of the transmission network. However, the integration of wind energy into the national energy mix remains conditional upon favorable wind conditions and the advancement of wind power technologies [3]-[5].
Variable-speed wind turbines offer the ability to capture the maximum possible energy from wind, which is inherently random and intermittent [6]. They have accelerated the development of large-scale wind turbine systems capable of contributing to the stability of transmission networks, thus justifying their use in large industrial applications. These systems are based on asynchronous and synchronous machines. Unlike wind turbines using asynchronous generators, those based on synchronous generators with a high number of pole pairs and permanent magnets (PMSG) operating at low speed do not require slip rings, brushes, or gearboxes. This significantly reduces maintenance costs, especially for offshore wind projects [7] [8]. Wind turbines operate in four distinct zones, but mainly in two: the partial-load zone for low wind speeds and the full-load zone for medium and high wind speeds.
To extract the maximum power during partial-load operation, wind turbines use Maximum Power Point Tracking (MPPT) techniques, which are divided into two categories. The speed-sensorless MPPT technique, based on wind speed estimation, allows indirect adjustment of the turbine’s mechanical rotational speed to its optimal value, enabling maximum power extraction for any wind speed below
VWN [9] [10]. However, this estimation-based strategy depends on the specific characteristics of the wind turbine, even though it remains simple, fast, independent of direct wind measurement, robust, and easy to implement [11]. It is widely used in direct-drive wind turbines to facilitate the development of vector control strategies for electrical quantities such as current or power [12] [13].
In contrast, during full-load operation, due to the mechanical load constraints of the turbine and the operating limits of the power converter in a PMSG designed for 100% of the rated power, the wind turbine system must be protected by a power-limitation mechanism [14]. Effective pitch angle control enables this limitation even under highly turbulent and strong wind conditions. A first-order control system is sufficient to achieve proper regulation of the blade pitch orientation.
Furthermore, any control strategy for the PMSG must contribute to improving its performance and, consequently, the overall performance of the entire system [15] [16]. Moreover, the E-115 3 MW variable-speed wind turbines supplied by ENERCON, featuring direct-drive systems and active pitch control, are of great interest for wind sites due to their high-power coefficient [17] [18].
In the context of harnessing Benin’s wind potential, the question arises of how to ensure the optimal operation of a 3 MW wind turbine under fluctuating wind conditions at a Beninese site and which wind technology is most suitable. The objective of this study is to evaluate the behavior of a 3 MW wind turbine with a Permanent Magnet Synchronous Generator (PMSG) under the wind conditions of Benin. To achieve this objective, we formulated two hypotheses. Modeling the wind turbine system allows for the evaluation of its behavior under variable wind conditions. The implemented control laws improve the overall performance of the wind turbine system.
2. Materials and Methods
2.1. Presentation of Research Diagram
Figure 1 presents the research diagram of the study. The process begins with the characterization of the E-115/3 MW wind turbine and the wind profile at a site in Benin. Subsequently, the wind profile, the wind turbine, and its shaft are modeled, along with the PMSG. Based on these models, a simple MPPT control strategy is developed, followed by the pitch angle regulation algorithm and the vector control of the PMSG. The capability of the wind turbine system and its control to deliver power under low, medium, and high wind conditions is then evaluated. The robustness of this controller will be assessed under realistic and strong wind conditions at a wind site in Benin.
Figure 1. Research diagram.
2.2. Presentation of the Wind Turbine Characteristics and Wind
Profile Data
Direct-drive technologies (without gearboxes) have the advantage of contributing to grid stability while reducing mechanical wear on components [16]. The EP5-type E115-3 MW direct-drive wind turbines with ring-type permanent magnet synchronous generators are an improved version of the EP3 turbines, which are directly coupled to a multipolar synchronous generator with separate excitation. The technical data of the E115-3 MW turbine are presented in Table 1.
Due to its proximity (38 km) to the 161 kV WAPP network in Sakété, an interconnection point, and its wind profile suitable for variable-speed wind turbines, the Sèkandji site has been selected for the study. Table 2 provides the parameters used for the wind speed simulation.
After modeling the turbine and the control laws, the technical data of the wind turbine and the wind profile will be used for simulation and validation of the results using Matlab/Simulink software.
Table 1. Technical data of the E-115/3 MW [17].
|
Technical Data of the Wind Turbine |
Characteristics |
Values |
Characteristics of the 115 m/3 MW wind
turbine |
Nominal Power |
3000 kW |
Starting speed |
2 m/s |
Nominal Wind speed |
11.5 m/s |
Wind speed at disconnection |
25 m/s |
Rotor blade diameter |
115.7 m |
Swept area of the roto |
10515.5 m2 |
Number of blades |
3 |
Rotational speed, max: |
12.8 U/min |
Tip speed |
78 m/s |
Hub height |
92 m |
Rotor rotational speed |
4 - 12.4 tours/min |
Reference pitch angle |
1˚ |
Table 2. Simulation parameters for the realistic wind speed at the Sèkandji site.
|
Characteristic data of the wind profile |
Parameters |
Value |
Unit |
Parameters of the realistic wind model |
Wind characteristics |
Average wind speed |
7.9544 |
m/s |
Scaling factor |
9.4669 |
m/s |
Shape factor |
1.9065 |
|
Average of variations |
0.1 |
m/s |
Data of the filter and the number-generating function |
Filter time constant
|
0.5 |
s |
Sampling period |
0.1 |
s |
2.3. Wind Modeling
A realistic wind at a given site is composed of two components: the mean wind
and the fluctuating component
, which represents variations caused by turbulent winds. Its profile is generated using Equation (2) [19].
(1)
(2)
where
is the magnitude of the function rand(t), which generates a sequence of random numbers between 0 and 1 following a uniform distribution and conforming to the Weibull distribution.
is the shape factor and
is the scale factor calculated from the Weibull breakdown at hub height.
represents the average wind variations at the site.
The very high-frequency components can be removed from the wind profile by filtering the wind speed generating function using Equation (3), which allows taking into account the mean velocity gradient across the entire surface of the turbine rotor [19].
(3)
2.4. Modeling of the Wind Turbine Dynamics
The modeling of the wind turbine system concerns the wind motion, the aerodynamic conversion of the turbine, and the shaft dynamics.
The wind turbine captures a portion of the power contained as kinetic energy in the wind passing through the area swept by the turbine. Betz’s theory, which allows the evaluation of the aerodynamic power
delivered by the wind turbine, is represented by Equation (4).
(4)
where
is the upstream wind speed (m/s).
is the power coefficient of the wind turbine (theoretical maximum value 0.593) [20].
is the air density (1.225 kg∙m−3).
is the area swept by the rotor blades (m2) with radius
(m).
The tip speed ratio, denoted by
, characterizes the ratio between the blade peripheral speed
and the wind speed
[20] and is given by Equation (5).
(5)
With
being the angular rotation speed of the blades in rd/s.
The aerodynamic efficiency
as a function of the pitch angle
and the tip speed ratio
is generally defined by empirical formulas. However, there are formulas implemented in MATLAB for large three-bladed turbines developed in [21] and [22], presented in Equation (6) and Equation (7).
(6)
(7)
is a coefficient used to adjust
for the aerodynamic power at the rated point provided by the manufacturer.
The wind torque
developed by a wind turbine rotating at the angular speed
is described by Equation (8) and Equation (9)
(8)
(9)
The generator shaft, mechanically linked to the rotor blades, converts mechanical energy into electrical energy. The mechanical torque seen by the directly coupled synchronous machine is different from that developed by the turbine. Indeed, the intrinsic viscous friction and the inertia of the wind turbine generate dynamics that oppose the turbine structure. This can be expressed mathematically by Equation (10).
(10)
With
: Electromagnetic torque produced by the generator.
: Viscous friction coefficient.
: Mechanical rotational speed of the generator.
: Turbine moment of inertia (kg∙m2).
The machine is modeled in the Park reference frame. Assuming that the rotor flux vector coincides with the “d” axis of the synchronous (d,q) frame and neglecting the homopolar component [16], the expressions for the output current vector components are given by Equation (11), and the electromagnetic torque by Equation (12)
(11)
(12)
Synchronous machines with a large number of poles have smooth poles, and the assumption given by relation (13).
(13)
Thus, the expression of the electromagnetic torque becomes (14)
(14)
where
is the number of pole pairs of the synchronous generator.
The angular frequency of the sinusoidal quantities is defined by (15)
(15)
2.5. Development of Control Lows
The control laws concern the strategies for controlling the dynamics of the wind turbine and the electrical conversion system to optimally exploit the potential of the wind site.
The MPPT technique without speed control based on wind speed estimation allows indirect adjustment of the mechanical rotation speed of the turbine to its optimal value, as indicated by the expression (16)
(16)
where
(17)
The proposed pitch controller is based on a proportional controller that amplifies the rotational speed error over a wide wind speed range from 1 m/s to 25 m/s. A memory placed after the controller allows the regulator to track the wind dynamics and adjust the pitch angle to optimize and limit the power to its rated value. Its algorithm is provided by the block diagram in Figure 2.
In vector control, the stator flux is oriented so as to independently control the electromagnetic torque, similar to a DC machine, by decoupling the (dq) components of the stator current [23]. The commonly used control consists of setting the direct-axis current to zero [24]. Thus, with the stator-induced flux
kept constant, the q-axis current provides the reference value of the electromagnetic torque according to Equation (18) [24]-[28].
Figure 2. Block diagram of the pitch controller.
(18)
The reference quadrature current is given by Equation (19).
(19)
The control voltages
and
are generated by compensating the (dq) coupling components of the voltages according to Equations (20) and (21).
(20)
(21)
The expressions in Equation (21) allow generating the reference voltages available at the output of the PI controllers of the direct and quadrature current regulation loops.
The PI controllers are synthesized using the pole placement technique for the closed-loop dynamics. The open-loop transfer function is presented in Equation (22) [29].
(22)
Using the pole compensation method, we obtain Equation (23) which gives the time constant of the original system.
(23)
The open-loop transfer function thus reduces to (24).
(24)
In closed loop, the expressions in Equation (25) are used, with
being the time constant of the reduced closed-loop system.
(25)
A response time
of the closed-loop system to reach 95% of the setpoint allows calculating the controller parameters, as indicated by the expressions in (26).
(26)
3. Results and Discussion
Three cases are simulated. In the first case, we set
with
and
for a nominal and ideal wind speed of 11.5 m/s. In the second case, the pitch angle is regulated by the controller we proposed. The third case considers three realistic wind speeds of 9 m/s, 11.5 m/s, and 22 m/s applied to the turbine at regular time intervals of 100 s. We focused on the evolution of the turbine rotational speed, the electromagnetic torque and power, the Id and Iq currents, and the AC voltages and currents.
3.1. Simulation Case 1
The curves in Figure 3 show the evolution of the rotational speed, electromagnetic torque, and electromagnetic power for a nominal wind speed of 11.5 m/s.
Figure 3. Evolution of the dynamic parameters with
for ideal wind.
It can be observed that the nominal values of the turbine speed and power are exceeded, reaching 1.5 rad/s for speed and over 4 MW for power. This phenomenon can be controlled by a pitch angle controller.
3.2. Simulation Case 2
The curves in Figure 4 show the evolution of the rotational speed, electromagnetic torque, and electromagnetic power for an ideal nominal wind speed of 11.5 m/s in the presence of the pitch controller.
Figure 4. Evolution of the dynamic parameters with pitch control for ideal wind.
In Figure 4, It can be observed that the nominal values of the turbine speed and power are kept within their limits: 1.348 rad/s for speed and 3 MW for power. The pitch controller thus not only optimizes the extracted power but also maintains the power at its nominal value. By keeping
for low and medium wind speeds, the power coefficient can be adjusted to achieve an aerodynamic power of 3 MW at 11.5 m/s. However, this approach does not allow optimizing the electrical power, which increases with aerodynamic efficiency at low values of β. This aspect is managed by the pitch angle controller. Nevertheless, it is necessary to evaluate the controller’s performance in limiting the power under strong and realistic wind dynamics.
3.3. Simulation Case 3
The curves in Figure 5 show the evolution of the rotational speed, electromagnetic torque, and electromagnetic power for three realistic wind speeds—9 m/s, 11.5 m/s, and 22 m/s—in the presence of the pitch controller over a period of 300 s.
Figure 5. Evolution of the dynamic parameters with pitch control for realistic wind.
In Figure 5, it can be observed that the nominal values of the turbine speed and power are maintained within their limits: 1.348 rad/s for speed and 3 MW for power. The pitch controller thus optimizes and limits the extracted power over a wide range of wind speeds, demonstrating the robustness of the pitch controller. This performance is achieved with a proportional gain of 26.765. However, the high gain value may alter the gain margin of the timing angle. This can be corrected by introducing a primary system or by using programmed gain correctors.
Figure 6. Evolution of the Id current and its reference in the PMSG.
Figure 6 shows that the Id current in the GSAP closely follows its zero-reference value, with peaks appearing at moments of sudden wind speed changes.
In Figure 7, the Iq current also follows its reference, displaying a pattern similar to the rotational speed. Indeed, the Iq current is proportional to the square of the speed due to the MPPT algorithm of the electromagnetic torque applied to the turbine. The Id and Iq currents evolving near their references demonstrate that the decoupling is well achieved, thus ensuring proper operation of the vector control.
Figure 7. Evolution of the Iq current and its reference in the PMSG.
Figure 8. Evolution of the three-phase currents Iabc in the PMSG.
Figure 8 and Figure 9 respectively show that the three-phase generator currents and the three-phase voltages follow the speed variations under low wind conditions. However, when the wind speed exceeds the nominal value, the pitch control comes into action and limits the electrical power to its nominal value. Consequently, the GSAP currents and voltages at variable frequency, decoupled from the grid, are limited to their maximum values.
Figure 9. Zoom on the three-phase currents Iabc in the PMSG.
Figure 10 and Figure 11, which zoom in on the phase currents and phase voltages, illustrate that these quantities are sinusoidal.
Figure 10. Evolution of the three-phase voltages Vabc at the PMSG output.
Figure 11. Zoom on the three-phase currents Vabc in the PMSG.
4. Conclusions
In the context of promoting wind energy development in Benin, we modeled the components of the energy conversion chain of a variable-speed wind turbine of the E115/3M type in order to evaluate its behavior under fluctuating wind conditions. The electrical converter is a high-pole permanent magnet generator. The control of the electromagnetic torque using an MPPT algorithm without speed regulation allowed extracting the maximum electrical power at low wind speeds. For realistic medium and high wind speeds, the action of the blade pitch regulation system, based on a proportional controller with memory, limited the power to 3 MW. The ability of the pitch controller to maintain the power at its nominal value despite wind speed fluctuations demonstrates its robustness, speed, and precision, ensuring stable and safe operation of the wind turbine. The pitch controller not only efficiently limits the power at full load and optimizes power at partial load, but it also has the advantage of being simple and less computationally demanding. Its simplicity facilitates integration into an optimization algorithm, particularly for stability studies.
However, the high gain value may alter the gain margin of the pitch angle. This can be corrected by introducing a first-order system or by using controllers with scheduled gains.
The synchronous generator is modeled in the two-phase reference frame. Vector control is used to decouple the voltages and align the flux along the quadrature axis. Validation tests showed that the Id and Iq currents evolve near their reference values, demonstrating that the decoupling is properly achieved and thus ensuring the correct operation of the vector control and the entire system.
Power optimization and safe operation at nominal power facilitate the integration of wind energy into the electricity mix of countries like Benin.