Optimal Power Reserve of a Wind Turbine System Participating in Primary Frequency Control

Participation of a wind turbine (WT) in primary frequency control (PFC) requires reserving some active power. The reserved power can be used to support the grid frequency. To maintain the required amount of reserve power, the WT is de-loaded to operate under its maximum power. The objective of this article is to design a control method for a WT system to maintain the reserved power of the WT, by controlling both pitch angle and rotor speed simultaneously in order to optimize the operation of the WT system. The pitch angle is obtained such that the stator current of the permanent magnet synchronous generator (PMSG) is reduced. Therefore, the resistive losses in the machine and the conduction losses of the converter are minimized. To avoid an excessive number of pitch motor operations, the wind forecast is implemented in order to predict consistent pitch angle valid for longer timeframe. Then, the selected pitch angle and the known curtailed power are used to find the optimal rotor speed by applying a nonlinear equation solver. To validate the proposed de-loading approach and control method, a detailed WT system is modeled in Matlab/Simulink (The Mathworks, Natick, MA, USA, 2017). Then, the proposed control scheme is validated using hardware-in-the-loop and real time simulation built in Opal-RT (10.4.14, Opal-RT Inc., Montreal, PQ, Canada).


Introduction
In power systems with less conventional generators (mainly synchronous), several functions must be added to the renewable resources to compensate for the grid requirements, such as frequency response, negative sequence, harmonics mitigation, and resynchronization and black-start ability [1].
For adjusting the frequency of a grid during frequency deviation, there are different control strategies and approaches that have been implemented and presented in the literature.These control methods can be categorized based on duration and capability with certain features and limitations [2][3][4].
The participation of wind turbines (WTs) in primary frequency control (PFC) has been discussed widely in several publications.Controllers based on rotor speed regulation or pitch angle regulation to de-load WT for PFC have been discussed in [5][6][7][8][9][10][11].
In [12,13], a coordination control was proposed between inertial, rotor speed, and pitch angle for WT, based on a doubly-fed induction generator (DFIG).This approach depends on a new classification of wind speed to select among the proposed controllers.Storing kinetic energy in the rotor shaft to be used during the PFC stage has been proposed [14,15].An adjustable droop controller was proposed to improve PFC in WT, based on a DFIG in [7].

Wind Turbine System
The system considered in this paper is direct-drive, based on PMSG.The WT system is decoupled from the power system using a back-to-back power converter.The model of the WT and PMSG used for this study are presented in [22][23][24][25].Detailed discussion of the model will be elaborated on in the following subsections.

Understanding the Wind Power Curve for Primary Frequency Control
To understand the relationship between pitch angle and rotor speed with the response of the torque and stator current of the machine, the model of the wind turbine and PMSG should be presented first.
The wind turbine model is a function of two main variables, the pitch angle (β) and the rotor speed (ω r ), as given by [25]: where A W is the swept area of the WT and the power coefficient C P (λ, β) is expressed as The coefficients c 1 to c 6 were taken from [25] (c 1 = 05176, c 2 = 116, c 3 = 0.4, c 4 = 5, c 5 = 21, c 6 = 0.0068), and the tip speed ratio is defined as Here, R w represents the radius of the wind turbine, and V W represents the wind speed.
From the curve characteristics of the mechanical power Equation (1) and the power coefficient of Equation (2), the maximum power is achieved at the zero pitch angle and optimum tip speed ratio for underrated wind speed conditions.Figure 1 shows the characteristics of the wind turbine power for different pitch angles and different rotor speeds for a specific wind turbine.As the pitch angle increases, the power produced by the WT decreases.This helps identify the operating point of the WT system for de-loading purposes.In addition, the rotor speed can be used to reach a specific production of power.To implement de-loading to a WT system, the power is curtailed by a certain amount, (∆P), using either the pitch angle or rotor speed controller.In Figure 1, the power is curtailed to (P curt.), providing a reserve power of ∆P.
amount, (Δ ), using either the pitch angle or rotor speed controller.In Figure 1, the power is curtailed to ( . ), providing a reserve power of Δ .The power curve can be represented in terms of rotor speed for different pitch angles for several wind speeds.This will allow us to select an optimal operating point of rotor speed and pitch angle for the purpose of PFC. Figure 1 represents a three-dimensional plot of the power curve.During regular operation of WT, the maximum power is tracked.However, for PFC, the sub-optimal operating point is selected to reserve the required power in the WT system.The output power is plotted using a flat surface, where the intersections of the two surfaces demonstrate all possible sets of pitch angles and rotor speed values that can be used.

Machine Losses and Rotor Speed
The machine considered in this paper is based on PMSG.The detailed model of the machine is presented in [23].The equation of the electromechanical torque is given as Assuming a round rotor type (i.e., = ), the quadrature current can be given as where is the number of poles of the machine and is the flux linkage.and represent the quadrature and direct inductance, respectively, and and represent the quadrature and direct currents of the d-q reference frame of the PMSG.
At constant power production, a higher rotor speed means lower torque.From Equation (5), it is clear that is proportional to the torque (i.e., lower torque results in lower current).This contributes to the resistive losses of the PMSG.The resistive (copper) losses in the stator equivalent circuit of the PMSG can be expressed as where represents the stator resistance of the PMSG.The copper losses are directly related to the current of the quadrature axis and the direct axis.Thus, by controlling the quadrature current of the machine ( = 0), the resistive losses can be reduced by involving pitch angle control for PFC.The power curve can be represented in terms of rotor speed for different pitch angles for several wind speeds.This will allow us to select an optimal operating point of rotor speed and pitch angle for the purpose of PFC. Figure 1 represents a three-dimensional plot of the power curve.During regular operation of WT, the maximum power is tracked.However, for PFC, the sub-optimal operating point is selected to reserve the required power in the WT system.The output power is plotted using a flat surface, where the intersections of the two surfaces demonstrate all possible sets of pitch angles and rotor speed values that can be used.

Machine Losses and Rotor Speed
The machine considered in this paper is based on PMSG.The detailed model of the machine is presented in [23].The equation of the electromechanical torque is given as Assuming a round rotor type (i.e., L q = L d ), the quadrature current can be given as where P is the number of poles of the machine and ψ m is the flux linkage.L q and L d represent the quadrature and direct inductance, respectively, and I q and I d represent the quadrature and direct currents of the d-q reference frame of the PMSG.At constant power production, a higher rotor speed means lower torque.From Equation (5), it is clear that I q is proportional to the torque (i.e., lower torque results in lower current).This contributes to the resistive losses of the PMSG.The resistive (copper) losses in the stator equivalent circuit of the PMSG can be expressed as where R S represents the stator resistance of the PMSG.The copper losses are directly related to the current of the quadrature axis and the direct axis.Thus, by controlling the quadrature current of the machine (I d = 0), the resistive losses can be reduced by involving pitch angle control for PFC.

Conduction Losses of Insulated Gate Bipolar Transistor Switches
The power converter used is based on insulated gate bipolar transistor (IGBT) switches.Power losses of such switches are conduction losses and switching losses.The conduction losses occur when the switch is turned on and conducts current.
To evaluate the conduction losses of the IGBT, a voltage drop and a small resistor are added into the series with an ideal switch.The simple expression of the conduction losses is given as where I s is the current crossing the IGBT switch during conducting, and V on and R on represent the voltage drop and the resistance of the switch, respectively The conduction losses of an IGBT depends on the temperature, which causes variation of the voltage drop across the switch.Also, the loss depends on the current passing through the switch.Therefore, decreasing the current through the switch will reduce the conduction losses.The conduction losses of the semiconductor devices can be evaluated using information from the data sheet from the manufacturer, considering the thermal model of the switch.

Wind Data and Speed Prediction
Using wind speed sensor provides a reference for the power tracking controller.With the availability of wind data, a wind forecast can be implemented to predict the wind speed.This allows us to estimate the behavior of the turbine and select the optimal operating point, ensuring control stability.Using a Light Detection and Ranging (LIDAR) system improves the accuracy of the measurements of wind speed.The National Renewable Energy Lab (NREL) investigated the use of LIDAR in feedforward control in WT systems [26].
Knowledge of future wind speed can be utilized in the algorithm to enhance wind turbine operation.In the literature, there are several approaches to estimate the average of the wind speed [27].Here, we implement the auto-regression model (AR) to generate and predict the mean of wind speed as where y t and y t−i represent the current and previous data (wind speed), respectively, Φ i is the autoregressive parameter of the model, and t is a Gaussian noise.The predicted wind speed is used to predict the future set points that required to maintain the reserved power for next timeframe (t + 1), as will be discussed in the next section.
Maximizing power using a speed controller is preferable during underrated wind speeds.However, in some cases the pitch controller can be involved to ensure smooth power production near rated wind speeds.In addition, combining a pitch angle controller with a speed controller can improve the efficiency of the machine, as will be discussed in the following section.

Proposed De-Loading Approach
The WT system considered in this paper is a direct drive configuration, where the power converter is directly connected to the machine's stator winding.In such a system, the rotor speed is fully controlled by the converter, while the blades are controlled by a pitch angle control system.
The reserve power for PFC is achieved by controlling both pitch angle and rotor speed simultaneously, in order to optimize the operation of the WT system.The goal is to minimize the power losses in the PMSG given in ( 6) and (7).This is accomplished by involving a pitch angle controller during underrated wind speeds.
Initially, with the prediction of wind speed, the pitch angle is obtained such that the rotor speed is increased to reduce the power losses in the stator of the PMSG.Predicted wind speed is averaged to avoid blade fluctuation and an excessive number of pitch motor operations.Therefore, averaged pitch angle is selected to ensure consistent operation of the WT, while minimizing resistive losses of the machine.
First, the measured wind speed is averaged and used to calculate the maximum power that can be achieved by the turbine.The wind forecast is needed to estimate the maximum possible power in advance (t + k) to help reduce the number of pitch operations.For a given specific reserve power, the output power is obtained and used to determine the required pitch angle instantaneously.
The optimum operating point of pitch angle and rotor speed is determined using the steps as follows: 1.
Calculate the current and future maximum available power as 2.
Calculate the output power as where P PFC is the output power of the droop curve and ∆P is the reserved power.

3.
Set β = 0 and solve the nonlinear power equation for ω r as: Increment β and repeat Step 3 until β max .

5.
Select an optimal solution such that: All solutions found by the solver in Step 3 provide the same output power as demonstrated in Figure 2.However, the optimal solution, S * = {β * , ω * r }, is selected such that pitch angle is as maximum as possible.As a result, the required rotor speed of the WT is increased and obtained by the nonlinear solver in Step 3.This will minimize the resistive power in the machine as in Equation ( 6) and the conduction losses in Equation (7).The algorithm must consider operation limits, such as Figure 3 summarizes the algorithm used to generate the optimum operating point of pitch angle and rotor speed.Note that the output of the droop curve can be involved in the equation, in order to define the rotor speed for the amount of power necessary to be injected during a frequency drop.Otherwise, the PFC power is considered to be zero.
All solutions found by the solver in Step 3 provide the same output power as demonstrated in Figure 2.However, the optimal solution, * = { * , * } , is selected such that pitch angle is as maximum as possible.As a result, the required rotor speed of the WT is increased and obtained by the nonlinear solver in Step 3.This will minimize the resistive power in the machine as in Equation ( 6) and the conduction losses in Equation (7).The algorithm must consider operation limits, such as  Figure 3 summarizes the algorithm used to generate the optimum operating point of pitch angle and rotor speed.Note that the output of the droop curve can be involved in the equation, in order to define the rotor speed for the amount of power necessary to be injected during a frequency drop.Otherwise, the PFC power is considered to be zero.
The proposed algorithm with the control loop is shown in Figure 4.The pitch controller is based on a conventional method, where the pitch angle is controlled considering the pitch servo transfer function [28].The rotor speed is controlled using a designed proportional integral (PI) controller with an inner loop to control the current.
The controller must ensure a smooth and reliable response to frequency drop.Therefore, the controller is required to adjust the frequency variation by considering the output of a droop curve.The proposed algorithm with the control loop is shown in Figure 4.The pitch controller is based on a conventional method, where the pitch angle is controlled considering the pitch servo transfer function [28].The rotor speed is controlled using a designed proportional integral (PI) controller with an inner loop to control the current.
The controller must ensure a smooth and reliable response to frequency drop.Therefore, the controller is required to adjust the frequency variation by considering the output of a droop curve.

Simulation and Hardware-in-the-Loop Results
The parameters of the system considered in this study including WT, PMSG are presented in Table 1.The wind turbine has a nominal rotor speed of 25.88 rpm at a 12 m/s rated wind speed.The rated power of the machine and WT is 2 MW.

Simulation and Hardware-in-the-Loop Results
The parameters of the system considered in this study including WT, PMSG are presented in Table 1.The wind turbine has a nominal rotor speed of 25.88 rpm at a 12 m/s rated wind speed.The rated power of the machine and WT is 2 MW.

Simulation Studies
For simulation study, different scenarios were implemented to achieve de-loading.The first scenario was based on varying the rotor speed only, while keeping pitch angle constant.The second scenario was performed using the proposed method, but allowing pitch angle variation without wind speed forecasting.In the last scenario, the proposed de-loading method was implemented to obtain an optimal pitch angle that is valid for a longer time frame, and updating the rotor speed instantaneously.
The simulation experiment was repeated again to perform de-loading for PFC purposes.In this scenario, the pitch angle was kept constant (i.e., 0 degrees), while the de-loading was performed using a rotor speed controller.The maximum power point tracking (MPPT) was modified to track the output power by subtracting the maximum power by the required power to be reserved.Then, the reference rotor speed was calculated to provide the required reserve power.
The proposed algorithm was used to de-load the WT, in order to maintain an adequate amount of power for PFC use.In this scenario, the wind profile was assumed to be predicted by wind forecast.The averaged estimated wind speed was used to find the optimum pitch angle for a longer period.Then, the current and future wind speed is used to estimate the optimal possible pitch angle for that time frame.The rotor speed was extracted using the nonlinear relation, using Equation (10).

Long-Term (One Day) Simulation
First, a simulation of one day of wind profile was performed.In this study, simplified model of the WT and PMSG were used to expedite the simulation.The wind profile used in this study is shown in Figure 5. Figure 6 shows the maximum power available and the de-loaded power.The dotted line represents the WT power achieved by rotor speed only.The WT power obtained by coordinating between rotor speed and pitch angle is presented by the dashed line.It is clear that both approaches provide the same amount of reserve power (∆P).In this study, the reserve is fixed at 0.3 MW for all wind profiles.time frame.The rotor speed was extracted using the nonlinear relation, using Equation (10).
4.1.1.Long-Term (One Day) Simulation First, a simulation of one day of wind profile was performed.In this study, simplified model of the WT and PMSG were used to expedite the simulation.The wind profile used in this study is shown in Figure 5. Figure 6 shows the maximum power available and the de-loaded power.The dotted line represents the WT power achieved by rotor speed only.The WT power obtained by coordinating between rotor speed and pitch angle is presented by the dashed line.It is clear that both approaches provide the same amount of reserve power (Δ ).In this study, the reserve is fixed at 0.3 MW for all wind profiles.Although every method gives the same reserve, the rotor speed is different, as shown in Figure 7.The rotor speed is increased by involving the pitch angle.This is because of the higher pitch angle that results in lower torque of the WT.Therefore, to maintain the same reserve, the rotor speed is increased.As a result, the quadrature current of the PMSG ( ) is reduced, as shown in Figure 8.In this simulation, we controlled the current of the direct axis to be zero.Therefore, the only current that appears in the losses is the current of the q-axis.From Figure 8, it is clear that with a 0 pitch angle, the current is higher for whole period considered in this simulation.Current produced by the proposed approach reduces the power losses in the stator of the machine and the conduction losses of the switch, as given in ( 6) and (7).In Figure 9, the pitch angle is represented using the dashed line.Although every method gives the same reserve, the rotor speed is different, as shown in Figure 7.The rotor speed is increased by involving the pitch angle.This is because of the higher pitch angle that results in lower torque of the WT.Therefore, to maintain the same reserve, the rotor speed is increased.As a result, the quadrature current of the PMSG (I q ) is reduced, as shown in Figure 8.In this simulation, we controlled the current of the direct axis to be zero.Therefore, the only current that appears in the losses is the current of the q-axis.From Figure 8, it is clear that with a 0 o pitch angle, the current is higher for whole period considered in this simulation.Current produced by the proposed approach reduces the power losses in the stator of the machine and the conduction losses of the switch, as given in ( 6) and (7).In Figure 9, the pitch angle is represented using the dashed line.
The proposed strategy has the ability to utilize a wind speed forecast and predict the optimal solution.In this case, one can combine the optimal operating point of the current and future state in order to minimize the operation of the pitch angle.For the same wind profile, one model included the wind forecast and the other model did not consider the wind prediction in the decision.The pitch angle changed 82 times for the model without a forecast.In contrast, in the case where there was a wind forecast, the pitch angle changed 64 times, as shown in Figure 10.In both cases, the power produced by the WT was the same.
increased.As a result, the quadrature current of the PMSG ( ) is reduced, as shown in Figure 8.In this simulation, we controlled the current of the direct axis to be zero.Therefore, the only current that appears in the losses is the current of the q-axis.From Figure 8, it is clear that with a 0 pitch angle, the current is higher for whole period considered in this simulation.Current produced by the proposed approach reduces the power losses in the stator of the machine and the conduction losses of the switch, as given in ( 6) and (7).In Figure 9, the pitch angle is represented using the dashed line.appears in the losses is the current of the q-axis.From Figure 8, it is clear that with a 0 pitch angle, the current is higher for whole period considered in this simulation.Current produced by the proposed approach reduces the power losses in the stator of the machine and the conduction losses of the switch, as given in ( 6) and (7).In Figure 9, the pitch angle is represented using the dashed line.The proposed strategy has the ability to utilize a wind speed forecast and predict the optimal solution.In this case, one can combine the optimal operating point of the current and future state in order to minimize the operation of the pitch angle.For the same wind profile, one model included the wind forecast and the other model did not consider the wind prediction in the decision.The pitch angle changed 82 times for the model without a forecast.In contrast, in the case where there was a wind forecast, the pitch angle changed 64 times, as shown in Figure 10.In both cases, the power produced by the WT was the same.solution.In this case, one can combine the optimal operating point of the current and future state in order to minimize the operation of the pitch angle.For the same wind profile, one model included the wind forecast and the other model did not consider the wind prediction in the decision.The pitch angle changed 82 times for the model without a forecast.In contrast, in the case where there was a wind forecast, the pitch angle changed 64 times, as shown in Figure 10.In both cases, the power produced by the WT was the same.

Short-Term Simulation
To test the wind system with the designed controller, a study was performed to validate the proposed method, using a detailed model of the whole system.The proposed WT system and control scheme shown in Figure 4 was modeled in Matlab/Simulink.
The aim of this short-term simulation was to validate the performance of the proposed algorithm, and examine its accuracy and response to wind variation.Therefore, the approach should give the reference for the pitch angle and rotor speed as quickly as possible.
In this simulation study, the wind speed was scaled to be in seconds instead of minutes.The wind turbine system was tested with the wind profile shown in Figure 11.Like the one-day simulation, in order to demonstrate the performance of the proposed de-loading algorithm, two different approaches were performed, with different methods of de-loading.First, the rotor speed only was used as a tool to maintain the reserved power.Then, the proposed de-loading approach was implemented using the same wind profile for the purpose of comparison.

Short-Term Simulation
To test the wind system with the designed controller, a study was performed to validate the proposed method, using a detailed model of the whole system.The proposed WT system and control scheme shown in Figure 4 was modeled in Matlab/Simulink.
The aim of this short-term simulation was to validate the performance of the proposed algorithm, and examine its accuracy and response to wind variation.Therefore, the approach should give the reference for the pitch angle and rotor speed as quickly as possible.
In this simulation study, the wind speed was scaled to be in seconds instead of minutes.The wind turbine system was tested with the wind profile shown in Figure 11.Like the one-day simulation, in order to demonstrate the performance of the proposed de-loading algorithm, two different approaches were performed, with different methods of de-loading.First, the rotor speed only was used as a tool to maintain the reserved power.Then, the proposed de-loading approach was implemented using the same wind profile for the purpose of comparison.Initially, the MPPT was implemented to be used as a baseline.Then two scenarios discussed in the previous section were examined.Figure 12 shows the maximum and de-loaded power for the two methods.The rotor speed for the different methods are demonstrated in Figure 13.Both approaches provide the same amount of reserve power (Δ ).As a one-day simulation, the reserve was fixed at 0.3 MW for all wind speeds.For the rotor speed method, the rotor speed was reduced compared to the other method.This low speed results in higher stator current, as shown in Figure 14.The participation of the pitch angle controller is demonstrated in Figure 15.Initially, the MPPT was implemented to be used as a baseline.Then two scenarios discussed in the previous section were examined.Figure 12 shows the maximum and de-loaded power for the two methods.The rotor speed for the different methods are demonstrated in Figure 13.Both approaches provide the same amount of reserve power (∆P).As a one-day simulation, the reserve was fixed at 0.3 MW for all wind speeds.For the rotor speed method, the rotor speed was reduced compared to the other method.This low speed results in higher stator current, as shown in Figure 14.The participation of the pitch angle controller is demonstrated in Figure 15.Initially, the MPPT was implemented to be used as a baseline.Then two scenarios discussed in the previous section were examined.Figure 12 shows the maximum and de-loaded power for the two methods.The rotor speed for the different methods are demonstrated in Figure 13.Both approaches provide the same amount of reserve power (Δ ).As a one-day simulation, the reserve was fixed at 0.3 MW for all wind speeds.For the rotor speed method, the rotor speed was reduced compared to the other method.This low speed results in higher stator current, as shown in Figure 14.The participation of the pitch angle controller is demonstrated in Figure 15.Initially, the MPPT was implemented to be used as a baseline.Then two scenarios discussed in the previous section were examined.Figure 12 shows the maximum and de-loaded power for the two methods.The rotor speed for the different methods are demonstrated in Figure 13.Both approaches provide the same amount of reserve power (Δ ).As a one-day simulation, the reserve was fixed at 0.3 MW for all wind speeds.For the rotor speed method, the rotor speed was reduced compared to the other method.This low speed results in higher stator current, as shown in Figure 14.The participation of the pitch angle controller is demonstrated in Figure 15.The high current on the stator causes an increase in the resistive losses and conduction losses, as demonstrated in Table 2.The table presents the averaged resistive losses and conduction losses for the whole period (240 s).The resistive losses are calculated using Equation ( 6), where the d-axis current is controlled to be fixed at 0 A. The IGBT characteristics were taken from the ASEA Brown Boveri (ABB) data sheet that matches the rating of the WT system, (5SNA 3600E170300 HiPak IGBT module).The conduction losses were calculated using web-based software (SEMIS) developed by Plexim (Zürich, Switzerland) for ABB data-sheet characteristics, in order to consider the detailed thermal module of the converter.

Real-Time and Hardware-in-Loop Results
To validate the response and the performance of the proposed approach with the designed controllers, hardware-in-the-loop (HIL) was implemented, as shown in Figure 16.The HIL study provides a real-time benchmark to test the proposed method and the designed controllers.The WT system is modeled and compiled inside an Opal-RT real time simulator.Then, C-code was generated for the controller of the machine side converter, and was implemented inside a Texas Instruments (TI) Digital signal processor (DSP), (TI TMS230F28335).The high current on the stator causes an increase in the resistive losses and conduction losses, as demonstrated in Table 2.The table presents the averaged resistive losses and conduction losses for the whole period (240 s).The resistive losses are calculated using Equation ( 6), where the d-axis current is controlled to be fixed at 0 A. The IGBT characteristics were taken from the ASEA Brown Boveri (ABB) data sheet that matches the rating of the WT system, (5SNA 3600E170300 HiPak IGBT module).The conduction losses were calculated using web-based software (SEMIS) developed by Plexim (Zürich, Switzerland) for ABB data-sheet characteristics, in order to consider the detailed thermal module of the converter.To validate the response and the performance of the proposed approach with the designed controllers, hardware-in-the-loop (HIL) was implemented, as shown in Figure 16.The HIL study provides a real-time benchmark to test the proposed method and the designed controllers.The WT system is modeled and compiled inside an Opal-RT real time simulator.Then, C-code was generated for the controller of the machine side converter, and was implemented inside a Texas Instruments (TI) Digital signal processor (DSP), (TI TMS230F28335).The controller (DSP) receives the rotor speed and the machine's current as an analog signal (fixed-sample) from Opal-RT.Then, the controller processes the signals and generates pulse width modulation (PWM) pulses to control the switches of the converter.The sampling frequency of the DSP is selected to be 9 kHz.The s-domain PI controllers were discretized using the backward Euler conversion method.The inner loop (current) controller is ensured to be faster than the outer loop The controller (DSP) receives the rotor speed and the machine's current as an analog signal (fixed-sample) from Opal-RT.Then, the controller processes the signals and generates pulse width modulation (PWM) pulses to control the switches of the converter.The sampling frequency of the DSP is selected to be 9 kHz.The s-domain PI controllers were discretized using the backward Euler conversion method.The inner loop (current) controller is ensured to be faster than the outer loop (speed) controller.
In this study, the controller and proposed method was tested with the wind speed profile shown in Figure 17.Like the simulation studies, two different approaches for de-loading were implemented, using the same wind profile for comparison.A baseline of MPPT was tested first, and the maximum power was achieved.Then, the de-loading was implemented using the proposed method summarized in Figure 3, using the optimal pitch angle and rotor speed set.Then, de-loading was achieved using rotor speed only.A fixed reserve power of 0.3 MW was set for both approaches.
The controller (DSP) receives the rotor speed and the machine's current as an analog signal (fixed-sample) from Opal-RT.Then, the controller processes the signals and generates pulse width modulation (PWM) pulses to control the switches of the converter.The sampling frequency of the DSP is selected to be 9 kHz.The s-domain PI controllers were discretized using the backward Euler conversion method.The inner loop (current) controller is ensured to be faster than the outer loop (speed) controller.
In this study, the controller and proposed method was tested with the wind speed profile shown in Figure 17.Like the simulation studies, two different approaches for de-loading were implemented, using the same wind profile for comparison.A baseline of MPPT was tested first, and the maximum power was achieved.Then, the de-loading was implemented using the proposed method summarized in Figure 3, using the optimal pitch angle and rotor speed set.Then, de-loading was achieved using rotor speed only.A fixed reserve power of 0.3 MW was set for both approaches.
For consistency, the same wind speed and time frame were used for all scenarios.First, the controller received the reference to perform the MPPT.Then, the same controller was fed with a rotor reference to de-load the WT system.Finally, the reference generated by the proposed method was used to obtain the reference set points (rotor speed and pitch angle).The results obtained from real time simulator and HIL are demonstrated in Figures 18-20.For the proposed experimental method, the three-phase current for machine is shown in Figure 18.The rotor speed (measured and reference) and q-axis current for all cases are shown in Figure 19.It is clear For consistency, the same wind speed and time frame were used for all scenarios.First, the controller received the reference to perform the MPPT.Then, the same controller was fed with a rotor reference to de-load the WT system.Finally, the reference generated by the proposed method was used to obtain the reference set points (rotor speed and pitch angle).
The results obtained from real time simulator and HIL are demonstrated in Figures 18-20.For the proposed experimental method, the three-phase current for machine is shown in Figure 18.The rotor speed (measured and reference) and q-axis current for all cases are shown in Figure 19.It is clear that the proposed method reduces the q-axis current.As a result, the power loss is decreased.The rotor speed in the plot is in (rad/s), and the gain is 1.For instant, the rotor speed at the first 5-s time slot was 2.3 (rad/s) for the MPPT case.The rotor speed that resulted from the proposed method is 2.2 (rad/s), whereas, the rotor speed is 1.6 (rad/s) when the pitch angle is kept at 0.
The gain for the q-axis current is 2000.Therefore, the current at the MPPT case is 3400 A for the first 5-s time slot.For the proposed and fixed pitch angle methods, the currents are 2600 A and 3600 A, respectively.Figure 20 shows the power generated from the WT system.Both approaches provide the required reserve power, as can be seen from the plot.The proposed method involves a pitch angle controller, as shown in Figure 20b.
(rad/s), whereas, the rotor speed is 1.6 (rad/s) when the pitch angle is kept at 0.
The gain for the q-axis current is 2000.Therefore, the current at the MPPT case is 3400 A for the first 5-s time slot.For the proposed and fixed pitch angle methods, the currents are 2600 A and 3600 A, respectively.Figure 20 shows the power generated from the WT system.Both approaches provide the required reserve power, as can be seen from the plot.The proposed method involves a pitch angle controller, as shown in Figure 20b.(rad/s), whereas, the rotor speed is 1.6 (rad/s) when the pitch angle is kept at 0. The gain for the q-axis current is 2000.Therefore, the current at the MPPT case is 3400 A for the first 5-s time slot.For the proposed and fixed pitch angle methods, the currents are 2600 A and 3600 A, respectively.Figure 20 shows the power generated from the WT system.Both approaches provide the required reserve power, as can be seen from the plot.The proposed method involves a pitch angle controller, as shown in Figure 20b.(rad/s), whereas, the rotor speed is 1.6 (rad/s) when the pitch angle is kept at 0. The gain for the q-axis current is 2000.Therefore, the current at the MPPT case is 3400 A for the first 5-s time slot.For the proposed and fixed pitch angle methods, the currents are 2600 A and 3600 A, respectively.Figure 20 shows the power generated from the WT system.Both approaches provide the required reserve power, as can be seen from the plot.The proposed method involves a pitch angle controller, as shown in Figure 20b.

Conclusions
This paper is aimed to design an optimal control method for a WT system intended to provide PFC to support power system stability.The proposed method helps maintain reserve power to be used to support the stability of the power system, while reducing power losses in the WT system.The reserved power of the WT is achieved by controlling both pitch angle and rotor speed simultaneously.The main contribution of the proposed method is to reduce the current of the PMSG, and as a result, the copper losses and conduction losses are reduced.Comprehensive case studies were implemented in simulation, as well as HIL with a real-time simulator.The results showed good performance of the proposed algorithm.The method is able to maintain the required reserve power for different wind speeds.The copper losses were calculated and compared with conventional method to de-load the power of WT.In addition, the conduction losses were evaluated using web-based software developed for ABB data sheet characteristics.The copper and conduction losses were reduced significantly compared to the conventional rotor speed control method.De-loading of a WT system under its maximum capacity means that some power is not being used.Therefore, implementing this optimal strategy to a large wind farm can improve the overall power generation by reducing a significant amount of power loss.

Figure 1 .
Figure 1.Power curve and output power versus rotor speed and pitch angle.

Figure 1 .
Figure 1.Power curve and output power versus rotor speed and pitch angle.

Figure 2 .
Figure 2. Demonstrating all possible operating points.

Figure 2 .
Figure 2. Demonstrating all possible operating points.

Figure 4 .
Figure 4. Wind turbine system with a control loop, including the proposed de-loading method.

Figure 4 .
Figure 4. Wind turbine system with a control loop, including the proposed de-loading method.

Figure 4 .
Figure 4. Wind turbine system with a control loop, including the proposed de-loading method.

Figure 6 .
Figure 6.Mechanical power for the maximum power point tracking (MPPT) and de-loaded scenarios.

Figure 10 .
Figure 10.The pitch angle with a wind forecast and without a wind forecast.

Figure 10 .
Figure 10.The pitch angle with a wind forecast and without a wind forecast.

Figure 12 .
Figure 12.Mechanical power for the MPPT and de-loaded scenarios (short-term simulation).

Figure 12 .
Figure 12.Mechanical power for the MPPT and de-loaded scenarios (short-term simulation).

Figure 12 .
Figure 12.Mechanical power for the MPPT and de-loaded scenarios (short-term simulation).

Figure 19 .
Figure19.Rotor speed and q-axis current (I q gain is 2000): (a) maximum power, (b) de-loading using proposed method, and (c) de-loading by controlling rotor speed only.

Table 2 .
Power losses for different scenarios.

Table 2 .
Power losses for different scenarios.