Next Article in Journal
Contribution to Active Thermal Protection Research—Part 1 Analysis of Energy Functions by Parametric Study
Previous Article in Journal
A Feasibility Analysis of a Solar Power Plant with Direct Steam Generation System in Sonora, Mexico
Previous Article in Special Issue
Optimal Cooperative Power Management Framework for Smart Buildings Using Bidirectional Electric Vehicle Modes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhanced Primary Frequency Control Using Model Predictive Control in Large-Islanded Power Grids with High Penetration of DFIG-Based Wind Farm

1
Engineering for Smart and Sustainable Systems Research Center, Mohammadia School of Engineers, Mohammed V University in Rabat, Rabat 10090, Morocco
2
Laboratory of Innovative Technologies, National School of Applied Sciences of Tangier, Abdelmalek Essaadi University, Tetouan 93000, Morocco
3
Department of Electrical and Control Engineering, Mokpo National University, Mokpo 58554, Republic of Korea
*
Author to whom correspondence should be addressed.
Energies 2023, 16(11), 4389; https://doi.org/10.3390/en16114389
Submission received: 24 April 2023 / Revised: 17 May 2023 / Accepted: 21 May 2023 / Published: 29 May 2023
(This article belongs to the Special Issue Towards Intelligent, Reliable and Flexible Stand-Alone Microgrids)

Abstract

:
A new primary frequency controller in power grids undergoing massive wind power penetration is the focus of this paper. The inescapable problem in a largely wind penetrated power grid is to ensure the maintenance of its frequency in the nominal band prescribed by the power system operator (PSO). However, with the massive arrival of wind farms with conventional control schemes, the operation of maintaining and restoring the frequency to the regulatory regimes remains very complicated. In order to overcome the above problem, this paper proposes a new strategy for primary frequency control in power grids using model predictive control (MPC) for a multi-cluster doubly-fed induction generator (DFIG)-based wind farm (WF), with a main objective of reducing the frequency nadir (FN), eliminating the second frequency dip (SFD), and providing the optimal support during wind speed variations. In this approach, a rolling prediction and optimization control strategy is developed based on the dynamic power system model to ideally predict the additional power to be provided. Moreover, in order to avoid second frequency dips, the wind turbines (WTs) are not allocated to extract additional power from the grid during the frequency event, the rotor speeds are not recovered to the maximum power point tracking (MPPT) operating points during the primary frequency control. The performance of the proposed controller was evaluated using a two-zone electrical system in MATLAB/Simulink®. The obtained results disclose that the frequency nadir is enhanced with more than 6.1% compared to the conventional schemes. In addition, the frequency response settling time has been improved with more than 10.51 s.

1. Introduction

Over the past few decades, a massive influx of large wind farms has taken place all over the world, and huge offshore wind farms have been especially integrated into power grids [1]. This introduces more and more adverse effects on the stability of the power system. One of the toughest issues faced is the impact of wind power generation on frequency control.

1.1. Literature Review

Power frequency control for wind farms so that they can be connected to the grid is a fundamental requirement in most countries around the world [2,3]. In the United Kingdom, the PSO requires a frequency droop in the 3 and 5% band for primary frequency control, while in Canada, a minimum 10% increase in active power during frequency events is required by the Independent Electricity System Operator. In Denmark, wind farms are required to adjust their output in the range of 10 to 100% during power events. EirGrid in Ireland requires a 15% reduction in active power during over-frequency events [2].
The literature review shows that a variety of research works have suggested a multitude of control strategies emulating the natural inertia of synchronous generators in WTs in order to mitigate frequency fluctuations during the most severe power system events. Indeed, studies in [4,5] have pointed out the rate of change of frequency (ROCOF) is significantly reduced by WT providing emulated inertial response. It is also shown that inertia emulation in WTs and droop control result in lower FN and faster recovery of the established frequency regime [6,7]. A combined control action of inertia emulation and the droop control method significantly improves the frequency response [8,9].
To achieve the contribution to frequency control, WTs increase or decrease their power outputs. This is achieved by operating the wind turbine in a deloaded operating point by reducing its pitch angle [3,7,8,9,10]. Meanwhile, this technique is recommended only for small and medium-sized wind turbines since the actuator time constant becomes large with increasing generator size, which adversely affects the frequency response. Thus, other studies have proposed to operate WTs at both suboptimal speed points [11,12]. Practically, frequency controllers based on wind turbine speed control have been shown to provide faster responses to system frequency [13]. However, all of these techniques introduce considerable loss given that they continuously operate generators at suboptimal points even before the frequency event. This means that the cash flow of the PSO is severely impaired.
In the same impetus, classical proportional integral (PID) and fractional order control have been proposed for angle of attack tuning and to improve the frequency response [14,15]. Fuzzy logic has also been used to solve the primary frequency tuning problem [16,17], but the results are highly dependent on the controller parameter settings. In addition, the studies in [18,19] develop a model predictive scheme to control the active power of the WF participating in primary frequency control. The development has been extended to a highly simplified power grid model that does not consider all elements of the power system. In the same perspective, interesting results have been obtained in [20,21], these studies used metaheuristic algorithms to build an optimized virtual rotor controller for energy systems connected to islanded microgrids. Table 1 synthesizes the strengths and weaknesses of the proposed control strategies in the literature.

1.2. Research Motivations

To design an optimal controller to addresses the drawbacks discussed in the literature review, it is legitimate to think of a control law producing an optimal output from a minimization of the quadratic cost function based on the dynamic model of the multi-zone power system and with various clusters of WTs operating at different wind speeds. The MPC technique considers the constraints on the output frequency deviation and load variation to compute the optimal control signal, as well as the constraints on the control signal. Indeed, a fast control strategy to increase the FN near the settlement frequency while ensuring the stable operation of the DFIG-based WTS under varying wind and load conditions is the objective of this paper. This objective is achieved through a novel MPC-based control scheme in which speeds are not rolled back to MPPT operating points during the primary frequency event support. The MPC controller measures the frequency deviation Δ f as an input and calculates the additional power based on the spatial state model of the power system frequency response. When the frequency stabilizes, Δ f does not cancel, the controller prevents the speed from recovering. In fact, the DFIG speed automatically recovers when the frequency is in its nominal band during the secondary frequency control; therefore, the MPC-based frequency control releases more deceleration kinetic energy compared to the generic control strategies. In addition, the FN in the proposed MPC control scheme is increased to be close to the establishment frequency without worrying about the SFD.

1.3. Contribution and Paper Organization

The main contributions of the proposed MPC-based fast frequency control are summarized as follows:
  • The FN is improved with more than 6.1% to be close to the settling frequency under all clusters wind and overload conditions, without SFD. This has never been achieved with existing conventional schemes, to our knowledge.
  • In order to meet the PSO requirements, the stabilization time to reach the new steady state frequency is significantly improved. It has been improved by more than 10.51 s compared to conventional schemes.
  • No electrical energy is required by the WTs to return the rotational speeds to the MPPT operating points during the frequency event, which positively improves the behavior of the power system with the proposed controller.
The remainder of this paper is structured as follows: Section 2 illustrates the theoretical modeling of the WTs, Section 3 details the proposed controller, In Section 4, the obtained results are discussed. Finally, Section 5 summarizes the conclusion.

2. Model of a DFIG-Based WT

2.1. Wind Turbine Model

The output power of a WT is expressed as follows [29]:
P m = 1 2 ρ A C p ( λ , β ) V ω 3
where ρ is the air density (in kg/m 3 ) and A is the section formulated by the turbine blades (in m 2 ). The power conversion coefficient, denoted by C p , is a function of two variables: the turbine pitch angle, measured in degrees, and the tip–speed ratio (TSR). The expression for λ is given by:
λ = ω r R V ω
The parameters ω r and R represent, respectively, the rotational speed of the shaft (in rad/s), and the radius formed by the turbine blades (in m). The power conversion coefficient C p is a nonlinear function of the turbine pitch angle (in degrees) and the tip–speed ratio λ , which is defined as the ratio of the speed of the blade tips to the speed of the wind that drives them. It is noteworthy that the coefficient C p has a specific maximum value at a particular tip–speed ratio λ O p t .
To characterize the performance of a wind turbine, a generalized equation is employed to model the power conversion coefficient C p ( λ , β ) . This equation is expressed as follows:
C p ( λ , β ) = c 1 c 2 λ i c 3 β c 4 e c 5 λ i + c 6 λ
1 λ i = 1 λ + 0.08 β 0.035 β 3 + 1
where c 1 , c 2 , c 3 , c 4 , c 5 , c 6 represent the WT’s manufacturer aerodynamic power coefficients.
In order to extract the optimal power that can be served by the WTs for given wind conditions, the speed controller operates with the MPPT characteristic depicted in Figure 1, which can be expressed by Equation (4), where K O p t is the optimal coefficient for the MPPT, whose value is given in Table A3.
According to the characteristic in Figure 1, the segment A B aligns with the starting zone, and the turbine is engaged in this zone. In the segment B C , which is called the optimal zone, the controller adjusts the rotational speed by regulating DFIG rotor current, and then modifies the power extracted from the turbine. Finally, in the C D segment, the speed is maintained nearly unchanged until the power attains its rated value. Beyond point D, the generator can no longer increase its power for protection reasons, and the pitch angle controller performs control maneuvers, which saturates the power output.
P m = K O p t ω 0 3 ω 0 ω r m i n ω r ω r m i n , ω r m i n ω r < ω 0 K O p t ω r 3 , ω 0 ω r ω 1 P m a x K O p t ω 1 3 ω r m a x ω 1 ω r ω r m a x , ω 1 < ω r ω r m a x

2.2. Power Reserve Deployment Requirements for Primary Frequency Control

Figure 2 portrays the basis of power system frequency behavior following a sudden loss of energy production or following an overloading event; it also illustrates the involved energy reserve levels. The grid frequency begins to drop, and only the instantaneous inertia reserves of the synchronous generators can renew the power balance in the system. Once a specific activation frequency is attained, the primary power reserves are triggered, this escalates the mechanical power of the generators until a different balance between the required electrical power and the generated mechanical power is reached. At this point, the frequency stops falling, it stabilizes on a new threshold below its nominal value.
Moreover, the FN and settling frequency are greatly determined by the slope of the frequency droop characteristic [35,36], the slower the response of the synchronous generator controller, the lower the FN.
Figure 2. Concepts definition (The designated values for time frames and frequencies comply with the recommendations of ENTSO-E) [37].
Figure 2. Concepts definition (The designated values for time frames and frequencies comply with the recommendations of ENTSO-E) [37].
Energies 16 04389 g002
The deployment of the secondary reserves rectifies the frequency levels to standard operational parameters, which triggers the deactivation of the primary reserves. The secondary reserves are implemented until they are wholly taken over by the tertiary reserves.

2.3. Emulation of Inertial Response in a DFIG-Based WT

It is possible that a DFIG-based WTS emulates the inertial response of traditional synchronous generators, since each rotating mass generates kinetic energy as it decelerates. To reach this target, it is essential to make some modifications to control scheme. In the aerodynamic characteristic of the DFIG-based WTS depicted in Figure 1, if the WT is only governed by the MPPT control curve and then an underfrequency event happens (e.g., due to a sudden overloading in the power grid or by loss of a generation unit); the conventional SGs are then called to raise their power production by making use of the primary energy reserves available—this will result in reducing the imbalance between the requested active power and the produced power.
The PSO requirements require WTS to participate in supporting grid contingencies; they are required to increase their output active power by moving immediately from the MPPT operating point located at E, to somewhere around point F in the left shifted MPPT as shown in Figure 1. Operating at point F will support the demand in active power and reduce the suffered imbalance. The active power demanded in F exceeds the mechanical power that can be delivered by the WT for the actual wind conditions; this will trigger rotor deceleration to the equilibrium point G ( ω r 0 ω r 1 ) . Throughout the (EFG) operation process, the DFIG-WT has been giving away a portion of its stored kinetic energy to reduce the FN. Although, going forward, it is working at a sub-optimal point G; the process of recovering the active power reference signal to its maximum point P E will occur during the secondary frequency regulation, passing by P G . Likewise, during the over-frequency event, the DFIG-WT inertial response describes the path (EHII’), except that this time, the rotor speed is restored to the MPPT point E directly at the end of the frequency disturbance.
It is infeasible to limit the deviation between the FN and the settling frequency if only the MPPT controller is operated. This is due to the fact that the speed recovery inevitably consumes additional energy from power grid to accelerate the turbine into MPPT point. In the MPC-based frequency control strategy, the FN is raised closer to the settling frequency since no speed recovery is allocated.
In the additional Δ f loop proposed in [38,39], the frequency responses remain limited. For good results, the MPPT and Δ f components are combined in the MPC-based frequency control strategy. Indeed, with the MPPT curve, the power reference will moderately decrease with the DFIG-WT deceleration, making the DFIG-WT behavior sufficiently stable. As depicted in Figure 1, during frequency support, the OPPT power curve used for limitation decreases with ω r from point F to point G. When supporting the frequency in the FG segment, the controller inherits MPPT proprieties, it has the privilege to adapt its operating point if the wind conditions change.
The analysis conducted above shows that the proposed controller allows the WTS power to decrease and stabilize smoothly at the power point G without introducing oscillations in the power system frequency, which cannot be provided by conventional SGs with highly oscillating characteristics.

3. Model Predictive Control for DFIG Supporting Frequency Control

3.1. DFIG Model

In stator voltage oriented control, the d axis of the d q reference frame is selected to carry the mains voltage vector v s , and the rotating reference frame d q , rotates at the synchronous speed ω s of the stator voltage. The DFIG dynamic equations are given by [40,41]:
v s k ( t ) = R s i s k ( t ) + d d t ϕ s k ( t ) + ω s ( t ) J ϕ s k ( t ) v r k ( t ) = R r i r k ( t ) + d d t ϕ r k ( t ) + ω s l ( t ) J ϕ r k ( t )
where, J = 0 1 1 0 , ω s l ( t ) = ω s ( t ) ω r ( t ) is the angular sliding pulsation, the vector x k ( t ) denotes the d q attributes of each variable, it is given by x k = ( x d , x q ) T .
The stator and rotor flux linkages are represented as [40,41]:
ϕ s k ( t ) = L s i s k ( t ) + L m i r k ( t ) ϕ r k ( t ) = L r i r k ( t ) + L m i s k ( t )
The magnitude of the stator voltage d-axis component v s d is identical to v s while the q-axis element v s q is constantly maintained at zero, v s d ( t ) = | | v s ( t ) | | and v s q ( t ) = 0 . The replacement of the value of ϕ r k ( t ) from Equation (6) to Equation (4) results in the expression of the rotor voltages as v r k ( t ) , as shown in Equation (7).
v r k ( t ) = R r i r k ( t ) + σ L r d d t i r k ( t ) + L m L s d d t ϕ s k ( t ) + ω s l ( t ) σ L r J i r k ( t ) + ω s l ( t ) L m L s J i s k ( t )
where σ = 1 L m 2 L s L r represents the DFIG leakage coefficient.
The expressions for the derivatives of the rotor currents are obtained in Equation (8). This is achieved with the aim of achieving current control.
d d t i r d ( t ) = 1 σ L s L r R r L s i r d ( t ) + ( ω s l ( t ) L s L r ω s ( t ) L m 2 ) i r q ( t ) + R s L m i s d ( t ) ω r ( t ) L m L s i s q ( t ) + L s v r d ( t ) L m v s d ( t ) d d t i r q ( t ) = 1 σ L s L r R r L s i r q ( t ) + ( ω s l ( t ) L s L r ω s ( t ) L m 2 ) i r d ( t ) + R s L m i s q ( t ) + ω r ( t ) L m L s i s d ( t ) + L s v r q ( t ) L m v s q ( t )
The DFIG’s rotor currents can be acquired directly from Equation (9) [42].
i r d ( t ) = 2 3 L s L m v s d ( t ) P s ( t ) i r q ( t ) = 2 3 L s L m v s d ( t ) Q s ( t ) v s d ( t ) ω s ( t ) L m
. where, P s ( t ) and Q s ( t ) denote the DFIG’s active and reactive powers, respectively.

3.2. RL Filter and Inverter Models

In the d q rotating frame, the RL filter model can be written as follows [42]:
v s k ( t ) = R f i f k ( t ) + L f d d t i f k ( t ) + ω s ( t ) L f J i f k ( t ) + v f k ( t )
The expressions for the derivatives of the filter current are derived in Equation (11).
d d t i f d ( t ) = R f L f i f d ( t ) + ω s ( t ) i f q ( t ) + 1 L f ( v s d ( t ) v f d ( t ) ) d d t i f q ( t ) = R f L f i f d ( t ) ω s ( t ) i f d ( t ) + 1 L f ( v s q ( t ) v f q ( t ) )
.
In the insulated-gate bipolar transistor (IGBT)-based three-phase voltage inverter illustrated in Figure 3, current can flow in both directions through six bi-directional switches (that can be controlled to be on or off) and anti-parallel diodes. All components are treated as ideal switches.
The switch states allow for expressing the voltages as:
v a n v b n v c n = v d c 3 2 1 1 1 2 1 1 1 2 S a S b S c
By coding the possible switching of the IGBT switches using three states ( S a , S b , S c ) , it becomes feasible to determine the voltage vector obtained in the stationary frame ( α , β ) for a particular switch combination. This information is then used to construct the switching polygon as depicted in Figure 4.

3.3. Finite Control Set-Model Predictive Control (FCS-MPC) for DFIG-Based Wind Turbine

The direct Euler in Equation (13) discrete-time model is employed for predicting currents in the future sampling periods [42,43]. The DFIG P + P 2 FCS-MPC discrete model at the first future sampling period can be expressed as illustrated in Equations (14) and (15) [42]:
d d t x ( t ) t = t k A x ( t k ) + B u ( t k ) = x ( t k + 1 ) x ( t k ) T s
i r d ( t k + 1 ) i r q ( t k + 1 ) = 1 0 0 1 i r d ( t k ) i r q ( t k ) + T s R r σ L s L r ω s l ( t k ) L s L r ω s ( t k ) L m 2 L s L r ω s l ( t k ) L s L r ω s ( t k ) L m 2 L s L r R r σ L s L r i r d ( t k ) i r q ( t k ) + T s σ L s L r R s L m ω r ( t k ) L m L s ω r ( t k ) L m L s R s L m + T s σ L r v r d ( t k ) v r q ( t k ) T s L m σ L r v s d ( t k ) v s q ( t k )
i f d ( t k + 1 ) i f q ( t k + 1 ) = 1 0 0 1 i f d ( t k ) i f q ( t k ) + T s R f L f ω s ( t k ) ω s ( t k + 1 ) R f L f i f d ( t k ) i f q ( t k ) + T s L f v s d ( t k ) v s q ( t k ) T s L f v f d ( t k ) v f q ( t k )
The model expressed by Equations (14) and (15) predicts the behavior of the DFIG rotor currents at a single sampling time, these values are reproduced in the discrete model to obtain the responses of the currents at two future sampling periods, the final model is expressed by Equations (16) and (17) [42].
i r d ( t k + 2 ) i r q ( t k + 2 ) = 1 0 0 1 i r d ( t k + 1 ) i r q ( t k + 1 ) + T s R r σ L s L r ω s l ( t k + 1 ) L s L r ω s ( t k + 1 ) L m 2 L s L r ω s l ( t k + 1 ) L s L r ω s ( t k + 1 ) L m 2 L s L r R r σ L s L r i r d ( t k + 1 ) i r q ( t k + 1 ) + T s σ L s L r R s L m ω r ( t k + 1 ) L m L s ω r ( t k + 1 ) L m L s R s L m + T s σ L r v r d ( t k + 1 ) v r q ( t k + 1 ) T s L m σ L r v s d ( t k + 1 ) v s q ( t k + 1 )
i f d ( t k + 2 ) i f q ( t k + 2 ) = 1 0 0 1 i f d ( t k + 1 ) i f q ( t k + 1 ) + T s R f L f ω s ( t k + 1 ) ω s ( t k + 1 ) R f L f i f d ( t k + 1 ) i f q ( t k + 1 ) + T s L f v s d ( t k + 1 ) v s q ( t k + 1 ) T s L f v f d ( t k + 1 ) v f q ( t k + 1 )
Extrapolations to a single and two sampling steps use Equations (18) and (19) [42].
i r e f k ( t k + 1 ) = 3 i r e f k ( t k ) 3 i r e f k ( t k 1 ) + i r e f k ( t k 2 )
i r e f k ( t k + 2 ) = 6 i r e f k ( t k ) 8 i r e f k ( t k 1 ) + 3 i r e f k ( t k 2 )
The best switching state is obtained by minimizing the objective function given in Equation (20). The first part of g represents the currents absolute error, while the second eliminates switching states resulting in exceeding the rated current of DFIG. The third term reduces the switching frequency of power converters. The term λ s w is set to 3.8 × 10 7 for the RSC and λ s w = 1.5 × 10 7 for the GSC [42].
Algorithm 1 and Figure 5 depict the proposed control scheme.
g R S C , G S C = | i r e f d ( t k + 2 ) i p d ( t k + 2 ) | + | i r e f q ( t k + 2 ) i p q ( t k + 2 ) | + 0 if i p d ( t k + 2 ) 2 + i p q ( t k + 2 ) 2 i m a x + if i p d ( t k + 2 ) 2 + i p q ( t k + 2 ) 2 > i m a x           + λ s w x = a , b , c ( S x p ( t k ) S x o p t ( t k ) )
Algorithm 1 FCS-MPC Rotor Currents Regulation Algorithm
  • MEASUREMENT ( x r k * ( t k ) , x f k * ( t k ) )
  • EXTRAPOLATION ( x r k * ( t k + 2 ) , x f k * ( t k + 2 ) )
  •      ( x r k * ( t k + 2 ) , x f k * ( t k + 2 ) ) Equation ( 19 ) .
  • MEASUREMENT ( x r k ( t k ) , x f k ( t k ) )
  • PREDICTION ( x r kp ( t k + 2 ) , x f kp ( t k + 2 ) )
  •     for i = 1,...,7
  • ( x r kp ( t k + 2 ) | x f kp ( t k + 2 ) ) (Equation (16) | Equation (17)).
  •      ( g i RSC | g i GSC ) (Equation (20)).
  •     end
  •      f sw RSC ( t k ) m i n g i RSC i = 1 , 2 , . . . , 7
  •      f sw GSC ( t k ) m i n g i GSC i = 1 , 2 , . . . , 7
  • RETURN f sw RSC ( t k ) , f sw GSC ( t k ) .

3.4. Continuous Control Set-Model Predictive Control (CCS-MPC) for DFIG-Based Wind Turbine Participation in Frequency Support

Currently, the WTS inertial response is obtained by introducing an additional power control loop. During frequency events, this loop increases the WT output power to meet power imbalance. The traditional controller consists of two separate loops, an artificial inertial loop, proposed in [44], its output control signal is proportional to the frequency derivative. The second loop is a droop-based characteristic, its control sets are proportional to the frequency deviation [45]. The control scheme is shown in Figure 6 [46]. K i n is a coefficient that scales the frequency derivative, whereas K D scales the frequency deviation. A high-pass filter on the droop loop is used in order to eliminate continuous activation of the control system by permanent frequency deviation and transient variations. Thus, a low-pass filter is placed on the frequency differentiator loop with the aim of ensuring that the frequency measurement noise does not affect the system. The requested power is regulated by the rotor side converter (RSC) controller by adjusting the rotor current i r d ( t k ) . It is supplied as a deceleration kinetic energy. In this way, the WT emulates the behavior of the traditional synchronous generator and virtually reinforce the total power system inertia.
In order to overcome all the mentioned limitations, the MPC uses the power system dynamics shown by Equation (21). The frequency deviation represents the frequency controller input, while the additional active power reference is its output. The model of frequency response illustrated in Figure 7 is composed of an active power response unit based on MPC which is added to a rotor speed response unit according to traditional MPPT.
Δ f ( t k + 1 ) = D 2 H Δ f ( t k ) + 1 2 H Δ P e ( t k ) 0 < Δ P e ( t k ) P O P P T P M P P T
The power system dynamic model is a sole input and a sole output, according to Equation (22) [47].
x m ( t k + 1 ) = A m x m ( t k ) + B m u ( t k ) y ( t k ) = C m x m ( t k ) + D m u ( t k )
where u ( t k ) is the control variable, y ( t k ) is the system output, and x m ( t k ) is the state variable.
Subtracting both sides of Equation (22) leads to;
x m ( t k + 1 ) x m ( t k ) = A m ( x m ( t k ) x m ( t k 1 ) ) + B m ( u ( t k ) u ( t k 1 ) )
The state and control variable differences are designated by
Δ x m ( t k + 1 ) = x m ( t k + 1 ) x m ( t k ) ; Δ x m ( t k ) = x m ( t k ) x m ( t k 1 ) ; u ( t k ) = u ( t k ) u ( t k 1 )
The state space model receives as input Δ u ( t k ) . In order to attach Δ x m ( t k ) to the output y ( t k ) , a new vector containing state variables is defined as x ( t k ) = Δ x m ( t k ) T y ( t k ) T
y ( t k + 1 ) y ( t k ) = C m ( x m ( t k ) x m ( t k 1 ) ) = C m Δ x m ( t k + 1 ) = C m A m Δ x m ( t k ) + C m B m Δ u ( t k )
This induces a new state space model, described in Equation (26).
Δ x m ( t k + 1 ) y ( t k + 1 ) x ( t k + 1 ) = A m o m T C m A m C m B m A · Δ x m ( t k ) y ( t k ) x ( t k ) + B m C m B m B · Δ u ( t k ) y ( t k ) = o m 1 C · Δ x m ( t k ) y ( t k )
where, o m = 0 0 . . . 0 n , n represents the state vector size. The augmented model (A, B, C) is formulated based on the system parameters. At sampling time k i , k i > 0 the state variable x ( t k i ) is directly measured, the current frequency information is provided by the state variable x ( t k i ) . The future control trajectory is denoted by
Δ u ( t k i ) , Δ u ( t k i + 1 ) , Δ u ( t k i + 2 ) , . . . , Δ u ( t k i + N c 1 )
where N C denotes the control horizon, and using the available x ( t k i ) , the state variables are predicted for a prediction horizon of N P . N P denotes the length of the optimization window. The incoming state variables are denoted by
x ( t k i + 1 | t k i ) , x ( t k i + 2 | t k i ) , . . . , x ( t k i + m | t k i ) , . . . , x ( t k i + N p | t k i )
where x ( t k i + m | t k i ) corresponds to the predicted state variable at t k i + m with the available x ( t k i ) information. Based on (A, B, C) the future control parameters and predicted states, the incoming state variables are computed sequentially in a stepwise manner.
x ( t k i + 1 | t k i ) = A x ( t k i ) + B Δ u ( t k i ) x ( t k i + 2 | t k i ) = A x ( t k i + 1 | t k i ) + B Δ u ( t k i + 1 ) . . . x ( t k i + N P | t k i ) = A N P x ( t k i ) + A N P 1 B Δ u ( t k i ) + A N P 2 B Δ u ( t k i + 1 ) + B Δ u ( t k i + 1 ) + , . . . , + A N P N C B Δ u ( t k i + N C 1 )
y ( t k i + 1 | t k i ) = CA x ( t k i ) + CB Δ u ( t k i ) y ( t k i + 2 | t k i ) = CA 2 x ( t k i ) + CAB Δ u ( t k i ) + CB Δ u ( t k i + 1 ) y ( t k i + 3 | t k i ) = CA 3 x ( t k i ) + CA 2 B Δ u ( t k i ) + CAB Δ u ( t k i + 1 ) + CB Δ u ( t k i + 1 ) . . . y ( t k i + N P | t k i ) = CA N P x ( t k i ) + CA N P 1 B Δ u ( t k i ) + CA N P 2 B Δ u ( t k i + 1 ) + B Δ u ( t k i + 1 ) + , . . . , + CA N P N C B Δ u ( t k i + N C 1 )
x ( t k i + m ) are formulated as information based on x ( t k i ) and Δ u ( t k i + j ) , where j = 0 , 1 , . . . N C 1 . Y and Δ U are defined as follows
Y = y ( t k i + 1 | t k i ) , y ( t k i + 2 | t k i ) , . . . , y ( t k i + N P | t k i ) T Δ U = Δ u ( t k i + 1 | t k i ) , Δ u ( t k i + 2 | t k i ) , . . . , Δ u ( t k i + N P | t k i ) T
Y and Δ U are N P and N C dimension. The final compact matrix form of the state model is
Y = F x ( t k i ) + ϕ Δ U
where;
F = CA CA 2 CA 3 . . . CA N P ϕ = CB 0 0 CAB CB 0 CA 2 B CAB 0 CA N P 1 B CA N P 2 B CA N P N C B
At t k i and a given reference r ( t k i ) , the MPC reduces the discrepancy between the predicted output and the reference. This is achieved by selecting the best Δ U control vector that minimizes the error function. The data vector containing the setpoint is
R s T = 1 1 . . . 1 N P r ( t k i ) = R ¯ r ( t k i )
The controller cost function J is defined by Equation (32).
J = ( R s Y ) T ( R s Y ) + Δ U T R ¯ Δ U
The first term consists of the squared error minimization, the significance of the value of Δ U is captured by the second term. R ¯ is a diagonal matrix as R ¯ = r ω I N C x N C . where r ω represents the optimization parameter. It ensures the desired closed loop performance. If r ω = 0 , the controller does not pay special attention to the output value of Δ U and only takes care of minimizing the error ( R s Y ) T ( R s Y ) . In the case of a large r ω , the Δ U values are smaller. To find the optimal Δ U that minimizes J, using Equation (30), J becomes
J = ( R s F x ( t k i ) ) T ( R s F x ( t k i ) ) 2 Δ U T ϕ T ( R s F x ( t k i ) ) R ¯ Δ U + Δ U T ( ϕ T ϕ + R ¯ ) Δ U
The optimal solution for the control signal expressed by Equation (35) is obtained at the point where J Δ U = 0 .
J Δ U = 2 ϕ T ( R s F x ( t k i ) ) + 2 ( ϕ T ϕ + R ¯ ) Δ U = 0
Δ U = ( ϕ T ϕ + R ¯ ) 1 ϕ T ( R s F x ( t k i ) )
With the existence and uniqueness of the Hessian matrix ( ϕ T ϕ + R ¯ ) 1 , the control signal can be calculated. Its optimal solution is related to the reference signal r ( t k i ) and the state variable x ( t k i ) as follows
Δ U = ( ϕ T ϕ + R ¯ ) 1 ϕ T ( R ¯ r ( t k i ) F x ( t k i ) )
Algorithm 2 illustrates the proposed control scheme for frequency support.
Since we are dealing with a controller that will have to equip each WT, the number of clusters formulating the wind farm are not affecting the performance of the proposed controller. However, the computational time is strongly impacted by the parameters N p and N c determining the control horizon chosen by the user.
Algorithm 2 CCS-MPC Frequency Control Algorithm
  • SET ( A m , B m , C m , N P , N C ) and MEASURE ( r ( t k i ) )
  • FORMULATION ( F , ϕ , R s )
  •     for i = 1,...,  N P
  •         R s 1 1 . . . 1 N P r ( t k i )
  •        for j = 1,..., N C
  •         ( F , ϕ ) Equation ( 31 ) .
  •        end
  •     end
  • MINIMIZATION ( J )
  •      Δ U ( ϕ T ϕ + R ¯ ) 1 ϕ T ( R ¯ r ( t k i ) F x ( t k i ) ) Equation ( 36 ) .
  •      Δ U ( t k i ) 1 0 . . . 0 N C Δ U .
  • RETURN U ( t k i ) Δ U ( t k i ) + U ( t k i 1 ) .

4. Case Studies

The Kundurs two-areas power grid model illustrated in Figure 7 [48] is used to demonstrate the effectiveness of the proposed MPC-based frequency control strategy. Ref. [48] details the synchronous generators (SGs) and automatic voltage regulators (AVRs) parameters. Furthermore, Figure 8 and Table 2 highlight the steam turbine governor model and its relevant constants, respectively [48,49]. SGs use the Type PSS4C multi-band power system stabilizer with the detailed model and parameters in [50]. Parameters of the synchronous generators and loads are illustrated in Table A1 and Table A2.
The DFIG-based WTS parameters are illustrated in Table A3. The WF consists of three aggregated clusters consisting of 100 DFIG-based WTs of 1.5 MW each. A 20/230 kV − 900 MVA transformer station and a 10 km long 230 kV transmission line connect the WF to the power grid.
Figure 1 shows the WT power characteristics, its input wind speed is 6 m/s, its rated speed is 12 m/s, the operating window of the DFIG rotor speed is between 0.7 p.u and 1.2 p.u.
Nevertheless, it is reminded that the contribution of this paper aims only at primary frequency control, no secondary frequency control means have been designed for SGs.
The performance of the proposed controller was evaluated under different wind with load variation conditions. The frequency controller was tested over the entire operating band of the DFIG to determine its stability. Three cases summarize this: operation with wind speed of 8 m/s at Cluster 1, 9 m/s at Cluster 2, 11 m/s, at Cluster 3 and 160 MW overload in the power grid. Furthermore, the controller was also tested during permanent wind variations. The efficiency of the proposed MPC controller is compared with the generic combined inertial and droop control, in which R is set to 5 % [26], and K i n is set to 2 * H [1]. In addition, it is compared with the sectional droop controller, in which the droop is dynamically adjusted according to Figure 9.

4.1. Case 1: Wind Speed of 8 m/s at Cluster 1, 9 m/s at Cluster 2, 11 m/s at Cluster 3 and 160 MW Overload in the Power Grid

Figure 10 and Figure 11 outline the obtained results during the simulation of case study 1.
Figure 10a show that with the proposed control scheme, the FN is arrested closely to the settling frequency, the obtained FN for MPPT operation, scheme #1, scheme #2, and the proposed scheme are evaluated to 49.8503, 49.8401, 49.8367, and 49.8324 Hz, respectively. Δ f is increased from −0.1676 to −0.1497 Hz as depicted in Figure 10b. The overproduction ensured by scheme #1, scheme #2 are almost the same during the initial stage, the overproduction is ensured for about 40s, then decreases rapidly with the frequency settling until reaching a much lower level and leave for the SGs that take over the excess demand. However, in the proposed controller, the overproduction slowly increases to reach its maximum level after about 8s, then slowly decreases and causes a deceleration of the DFIG-Based WTS, which restores more kinetic energy compared to conventional schemes.
In scheme #1 and scheme #2, the rapid cancellation of overproduction lead to the rotor speed recovery process to the MPPT points (Figure 10d and Figure 11a), which are ω r = 0.8 p.u in Cluster 1, ω r = 0.9 p.u in Cluster 2, and ω r = 1.1 p.u in Cluster 3. In the proposed scheme, this recovery is not allocated, the rotor speed converges to ω r = 0.7088 p.u in Cluster 1, ω r = 0.7539 p.u in Cluster 2, ω r = 0.9690 p.u in Cluster 3, and stabilizes at these points until the frequency is restored to its nominal band during the secondary control. It is this deceleration that allowed us to restore enough kinetic energy to the power grid and arrest the FN at the value closest to settling frequency which amounts to 49.8761 Hz. It can also be observed in Figure 10a,b that, with the proposed scheme, no oscillations were experienced in the power system frequency and the WTs output power. While with scheme #1 and scheme #2, the output power and frequency have remarkable oscillations and take longer to achieve the new steady state.
Figure 10b also illustrates that the WT rotor speed does not oscillate, which means the alleviation of mechanical stresses exerted on the WTs shafts. This is significantly important for power systems experiencing high WFs penetration with very large capacities. Figure 11a,b show that the stabilized operation point under the proposed controller is slightly lower than that of MPPT operation, as with the proposed controller, the WTs continue to operate at the OPPT point. However, the settling frequency difference between the two points is only 0.0064 Hz, which is very minimal. Meanwhile, the settling time for reaching the new frequency steady state is significantly improved by 10.51 s by the proposed controller compared to the conventional control, the settling time of scheme #1 and scheme #2 and the proposed scheme are evaluated to 43.5954, 42.3113, and 31.7937 s, respectively. Obtained results are summarized in Table 3 and Figure 11c,d.

4.2. Case 2: Variable Wind Speeds and Fixed Load

The obtained results during the simulation of case study 2 are shown in Figure 12.
Remarkable fluctuations are caused by changes in weather conditions, especially wind speed. In order to evaluate the robustness to these data, a variation of ±0.5 m/s was introduced in the different sites as shown in Figure 12a. It can be observed that the proposed MPC predictive controller damped the oscillation of the system frequency considerably compared to the traditional controllers. This is due to the fact that the output of the MPC controller is optimally calculated following the actual deviation caused by the change in wind speed while the traditional controllers rely on well-defined droop and coefficients beforehand. Furthermore, the results show that traditional controllers have more irregular frequency responses. The performance of the different controllers was evaluated using a set of performance indicators as expressed by Equation (37) [51].
Obtained results are shown in the Figure 12 and Table 4.
I S E = 0 ( f f n o m ) 2 d t I T S E = 0 t · ( f f n o m ) 2 d t I A E = 0 | ( f f n o m ) | d t I T A E = 0 t · | ( f f n o m ) | d t

5. Conclusions

This paper proposes a fast MPC-based frequency support controller using the the WTs rotational stored kinetic energy in multi-clusters wind farm. The suggested scheme raises the FN closely to the settling frequency, which reduces the mechanical stresses on the WTs during frequency events. The key novelties of the proposed scheme are listed as follows:
  • The frequency Nadir enhancement: The FN is arrested closely to the settling frequency. Indeed, the FN is improved by more than 6.1% compared to conventional schemes.
  • Eliminating SFDs through the WTS’s rotor speed recovery approach: WTS rotor speed is not designed to revert back to the MPPT point during the main frequency support. This means that the WT does not consume any additional power from the grid for speed recovery to MPPT during the frequency event. Therefore, no SDF will occur.
  • Improving the frequency settling time: the fast MPC-based controller has enhanced the settling time for reaching the new frequency steady state by more than 10.51 s, compared to conventional schemes.
  • Improving the frequency response during transient wind conditions and fluctuations.
Regarding the difference in wind conditions between sites on controller performance, the stability of the WTs is ensured for all the operating band.

Author Contributions

Conceptualization, Y.A.A. and M.O.; methodology, Y.A.A. and M.O.; software, Y.A.A.; validation, Y.A.A., M.O., S.-H.L. and Z.C.; formal analysis, Y.A.A. and Z.C.; investigation, S.-H.L. and Z.C.; resources, M.O., Z.C. and S.-H.L.; writing—original draft preparation, Y.A.A. and M.O.; writing—review and editing, Y.A.A., M.O., Z.C. and S.-H.L.; supervision, M.O.; project administration, M.O.; funding acquisition S.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the Research Centre for Engineering in Smart and Sustainable Systems of Mohammed V University in Rabat, Morocco. It is was also supported by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) Grant funded by the Korea Government (MOTIE) in part (Development of DC Power Trade Platform System in Public Community Which is Connected by EV-Renewable Based on Block Chain Technology) under Grant 20192010107050, and in part (20223A10100030, Development of Synchronous Condenser Model and Power System Inertia Operating Technology).

Data Availability Statement

The data used in this research are available in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CCS-MPCContinuous Control Set Model Predictive Control.
DFIGDoubly-Fed Induction Generator.
FCS-MPCFinite Control Set Model Predictive Control.
FNFrequency Nadir.
GSCGrid Side Converter.
IGBTInsulated-Gate Bipolar Transistor.
MPPTMaximum Power Point Tracking.
OPPTOptimal Power Point Tracking.
PSOPower System Operator.
ROCOFRate Of Change Of Frequency.
RSCRotor Side Converter.
SFDSecond Frequency Dip.
WFWind Farm.
WTWind Turbine.
f , f n o m Actual and nominal system frequency.
Δ f Frequency deviation.
d f / d t Frequency derivative.
P m WT mechanical output power.
v ω WT input Wind speed.
β Pitch angle.
ω r DFIG rotor speed.
ω r m i n Cut-in rotor speed.
ω r m a x Rated rotor speed.
K O p t Optimal coefficient for MPPT control.
K i n , K D ROCOF and droop loop gains.
HWT inertia constant.
v s k , i s k dq-axis stator voltage and current components.
v r k , i r k dq-axis rotor voltage and current components.
v f k , i f k dq-axis filter voltage and current components.
t k Current sampling time.

Appendix A

Table A1. Operating parameters of synchronous generator.
Table A1. Operating parameters of synchronous generator.
ParameterValue
Rated Power P e (MVA) 900.00
D-Axis Synchronous Reactance X d (p.u) 1.8000
Q-Axis Synchronous Reactance X q (p.u) 1.7000
Stator Resistance R a (p.u) 0.0025
Stator Leakage Inductance X l (p.u) 0.2000
D-Axis Transient Reactance X d (p.u) 0.3000
Q-Axis Transient Reactance X q (p.u) 0.5500
D-Axis Subtransient Reactance X d (p.u) 0.2500
Q-Axis Subtransient Reactance X q (p.u) 0.2500
D-Axis Transient Open Circuit Time Constant T d 0 (s) 8.0000
Q-Axis Transient Open Circuit Time Constant T q 0 (s) 0.4000
D-Axis Subtransient Open Circuit Time Constant T d 0 (s) 0.0300
Q-Axis Subtransient Open Circuit Time Constant T q 0 (s) 0.0500
Total inertia constant H (s) 6.5000

Appendix B

Table A2. Loads and shunt capacitance connected in Kundur’s AC grid.
Table A2. Loads and shunt capacitance connected in Kundur’s AC grid.
BusLoadShunt Capacitance
Bus B 7 ( 900 + j 100 ) MVA200 MVAr
Bus B 8 ( 900 + j 100 ) MVA350 MVAr

Appendix C

Table A3. Parameters of DFIG-based WTS.
Table A3. Parameters of DFIG-based WTS.
ParameterValue
Rated Power P e (MVA) 1.50000
Rated Stator Voltage v s (V) 575.000
Rated Frequency f (Hz) 50.0000
Stator Resistance R s (p.u) 0.00706
Rotor Resistance R r (p.u) 0.00500
Stator Inductance L s (p.u) 0.17100
Rotor Inductance L r (p.u) 0.15600
Mutual Inductance L m (p.u) 2.90000
Pole Pairs p3
Leakage constant σ (p.u) 1 L m 2 L s L r
Converters rated power (p.u) 0.5
Total inertia constant H (s) 5.04000
ω r m i n (p.u) 0.70000
ω r m a x (p.u) 1.20000
[ β m i n , β m a x ][ 0.00000 , 45.0000 ]
Maximum d β d t ( /s) 2.00000
[ c 1 , c 2 , c 3 , c 4 , c 5 , c 6 ][ 0.5176 , 116, 0.4 , 5, 21, 0.0068 ]
v n o m (m/s) at P e = 0.73 (p.u) 12.0000
K O p t 0.42250
K i n 2 H
K D 1 R ( % )
Pitch controller gain K p 500

References

  1. Yingcheng, X.; Nengling, T. System frequency regulation in doubly fed induction generators. Int. J. Electr. Power Energy Syst. 2012, 43, 977–983. [Google Scholar] [CrossRef]
  2. Datta, U.; Shi, J.; Kalam, A. Primary frequency control of a microgrid with integrated dynamic sectional droop and fuzzy based pitch angle control. Int. J. Electr. Power Energy Syst. 2019, 111, 248–259. [Google Scholar] [CrossRef]
  3. Tian, X.; Wang, W.; Chi, Y.; Li, Y.; Liu, C. Virtual inertia optimisation control of DFIG and assessment of equivalent inertia time constant of power grid. IET Renew. Power Gener. 2018, 12, 1733–1740. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Melin, A.M.; Djouadi, S.M.; Olama, M.M.; Tomsovic, K. Provision for Guaranteed Inertial Response in Diesel-Wind Systems via Model Reference Control. IEEE Trans. Power Syst. 2018, 33, 6557–6568. [Google Scholar] [CrossRef]
  5. Liu, F.; Liu, Z.; Mei, S.; Wei, W.; Yao, Y. ESO-Based Inertia Emulation and Rotor Speed Recovery Control for DFIGs. IEEE Trans. Energy Convers. 2017, 32, 1209–1219. [Google Scholar] [CrossRef]
  6. Van de Vyver, J.; De Kooning, J.D.M.; Meersman, B.; Vandevelde, L.; Vandoorn, T.L. Droop Control as an Alternative Inertial Response Strategy for the Synthetic Inertia on Wind Turbines. IEEE Trans. Power Syst. 2016, 31, 1129–1138. [Google Scholar] [CrossRef]
  7. Tan, Y.; Meegahapola, L.; Muttaqi, K.M. A Suboptimal Power-Point-Tracking-Based Primary Frequency Response Strategy for DFIGs in Hybrid Remote Area Power Supply Systems. IEEE Trans. Energy Convers. 2016, 31, 93–105. [Google Scholar] [CrossRef]
  8. Wang, Z.; Wu, W. Coordinated Control Method for DFIG-Based Wind Farm to Provide Primary Frequency Regulation Service. IEEE Trans. Power Syst. 2018, 33, 2644–2659. [Google Scholar] [CrossRef]
  9. Ghosh, S.; Kamalasadan, S.; Senroy, N.; Enslin, J. Doubly Fed Induction Generator (DFIG)-Based Wind Farm Control Framework for Primary Frequency and Inertial Response Application. IEEE Trans. Power Syst. 2016, 31, 1861–1871. [Google Scholar] [CrossRef]
  10. Wang, S.; Hu, J.; Wang, S.; Tang, H.; Chi, Y. Comparative study on primary frequency control schemes for variable-speed wind turbines. J. Eng. 2017, 2017, 1332–1337. [Google Scholar] [CrossRef]
  11. Zhang, X.; Zha, X.; Yue, S.; Chen, Y. A Frequency Regulation Strategy for Wind Power Based on Limited Over-Speed De-Loading Curve Partitioning. IEEE Access 2018, 6, 22938–22951. [Google Scholar] [CrossRef]
  12. Gomez, L.A.G.; Lourenco, L.F.N.; Grilo, A.P.; Salles, M.B.C.; Meegahapola, L.; Filho, A.J.S. Primary Frequency Response of Microgrid Using Doubly Fed Induction Generator With Finite Control Set Model Predictive Control Plus Droop Control and Storage System. IEEE Access 2020, 8, 189298–189312. [Google Scholar] [CrossRef]
  13. Boyle, J.; Littler, T.; Muyeen, S.; Foley, A.M. An alternative frequency-droop scheme for wind turbines that provide primary frequency regulation via rotor speed control. Int. J. Electr. Power Energy Syst. 2021, 133, 107219. [Google Scholar] [CrossRef]
  14. Mahvash, H.; Taher, S.A.; Rahimi, M.; Shahidehpour, M. Enhancement of DFIG performance at high wind speed using fractional order PI controller in pitch compensation loop. Int. J. Electr. Power Energy Syst. 2019, 104, 259–268. [Google Scholar] [CrossRef]
  15. Ibrahim, A.O.; Nguyen, T.H.; Lee, D.C.; Kim, S.C. A Fault Ride-Through Technique of DFIG Wind Turbine Systems Using Dynamic Voltage Restorers. IEEE Trans. Energy Convers. 2011, 26, 871–882. [Google Scholar] [CrossRef]
  16. Zhang, S.; Mishra, Y.; Shahidehpour, M. Fuzzy-Logic Based Frequency Controller for Wind Farms Augmented With Energy Storage Systems. IEEE Trans. Power Syst. 2016, 31, 1595–1603. [Google Scholar] [CrossRef]
  17. Miao, L.; Wen, J.; Xie, H.; Yue, C.; Lee, W.J. Coordinated Control Strategy of Wind Turbine Generator and Energy Storage Equipment for Frequency Support. IEEE Trans. Ind. Appl. 2015, 51, 2732–2742. [Google Scholar] [CrossRef]
  18. Mohamed, T.H.; Morel, J.; Bevrani, H.; Hiyama, T. Model predictive based load frequency control_design concerning wind turbines. Int. J. Electr. Power Energy Syst. 2012, 43, 859–867. [Google Scholar] [CrossRef]
  19. Liu, J.; Yao, Q.; Hu, Y. Model predictive control for load frequency of hybrid power system with wind power and thermal power. Energy 2019, 172, 555–565. [Google Scholar] [CrossRef]
  20. Abubakr, H.; Vasquez, J.C.; Hassan Mohamed, T.; Guerrero, J.M. The concept of direct adaptive control for improving voltage and frequency regulation loops in several power system applications. Int. J. Electr. Power Energy Syst. 2022, 140, 108068. [Google Scholar] [CrossRef]
  21. Abubakr, H.; Guerrero, J.M.; Vasquez, J.C.; Mohamed, T.H.; Mahmoud, K.; Darwish, M.M.F.; Dahab, Y.A. Adaptive LFC Incorporating Modified Virtual Rotor to Regulate Frequency and Tie-Line Power Flow in Multi-Area Microgrids. IEEE Access 2022, 10, 33248–33268. [Google Scholar] [CrossRef]
  22. Wu, Z.; Gao, W.; Gao, T.; Yan, W.; Zhang, H.; Yan, S.; Wang, X. State-of-the-art review on frequency response of wind power plants in power systems. J. Mod. Power Syst. Clean Energy 2018, 6, 1–16. [Google Scholar] [CrossRef]
  23. Dreidy, M.; Mokhlis, H.; Mekhilef, S. Inertia response and frequency control techniques for renewable energy sources: A review. Renew. Sustain. Energy Rev. 2017, 69, 144–155. [Google Scholar] [CrossRef]
  24. Hwang, M.; Muljadi, E.; Jang, G.; Kang, Y.C. Disturbance-Adaptive Short-Term Frequency Support of a DFIG Associated With the Variable Gain Based on the ROCOF and Rotor Speed. IEEE Trans. Power Syst. 2017, 32, 1873–1881. [Google Scholar] [CrossRef]
  25. Yang, D.; Gao, H.C.; Zhang, L.; Zheng, T.; Hua, L.; Zhang, X. Short-term frequency support of a doubly-fed induction generator based on an adaptive power reference function. Int. J. Electr. Power Energy Syst. 2020, 119, 105955. [Google Scholar] [CrossRef]
  26. Margaris, I.D.; Papathanassiou, S.A.; Hatziargyriou, N.D.; Hansen, A.D.; Sorensen, P. Frequency Control in Autonomous Power Systems With High Wind Power Penetration. IEEE Trans. Sustain. Energy 2012, 3, 189–199. [Google Scholar] [CrossRef]
  27. Liu, T.; Pan, W.; Quan, R.; Liu, M. A Variable Droop Frequency Control Strategy for Wind Farms that Considers Optimal Rotor Kinetic Energy. IEEE Access 2019, 7, 68636–68645. [Google Scholar] [CrossRef]
  28. Datta, U.; Kalam, A.; Shi, J. Frequency performance analysis of multi-gain droop controlled DFIG in an isolated microgrid using real-time digital simulator. Eng. Sci. Technol. Int. J. 2020, 23, 1028–1041. [Google Scholar] [CrossRef]
  29. Hwang, M.; Muljadi, E.; Park, J.W.; Sorensen, P.; Kang, Y.C. Dynamic Droop–Based Inertial Control of a Doubly-Fed Induction Generator. IEEE Trans. Sustain. Energy 2016, 7, 924–933. [Google Scholar] [CrossRef]
  30. Abazari, A.; Monsef, H.; Wu, B. Load frequency control by de-loaded wind farm using the optimal fuzzy-based PID droop controller. IET Renew. Power Gener. 2019, 13, 180–190. [Google Scholar] [CrossRef]
  31. Abazari, A.; Dozein, M.G.; Monsef, H. An optimal Fuzzy-logic based frequency control strategy in a high wind penetrated power system. J. Frankl. Inst. 2018, 355, 6262–6285. [Google Scholar] [CrossRef]
  32. Elyaalaoui, K.; Ouassaid, M.; Cherkaoui, M. Primary frequency control using hierarchal fuzzy logic for a wind farm based on SCIG connected to electrical network. Sustain. Energy Grids Netw. 2018, 16, 188–195. [Google Scholar] [CrossRef]
  33. Sambariya, D.; Nath, V. Load Frequency Control Using Fuzzy Logic Based Controller for Multi-area Power System. Br. J. Math. Comput. Sci. 2016, 13, 1–19. [Google Scholar] [CrossRef] [PubMed]
  34. Obaid, Z.A.; Cipcigan, L.; Muhssin, M.T. Fuzzy hierarchal approach-based optimal frequency control in the Great Britain power system. Electr. Power Syst. Res. 2016, 141, 529–537. [Google Scholar] [CrossRef]
  35. Kayikci, M.; Milanovic, J. Dynamic Contribution of DFIG-Based Wind Plants to System Frequency Disturbances. IEEE Trans. Power Syst. 2009, 24, 859–867. [Google Scholar] [CrossRef]
  36. Sun, Y.z.; Zhang, Z.s.; Li, G.j.; Lin, J. Review on frequency control of power systems with wind power penetration. In Proceedings of the 2010 International Conference on Power System Technology, Hangzhou, China, 24–28 October 2010; pp. 1–8. [Google Scholar] [CrossRef]
  37. Díaz-González, F.; Hau, M.; Sumper, A.; Gomis-Bellmunt, O. Participation of wind power plants in system frequency control: Review of grid code requirements and control methods. Renew. Sustain. Energy Rev. 2014, 34, 551–564. [Google Scholar] [CrossRef]
  38. Keung, P.K.; Li, P.; Banakar, H.; Ooi, B.T. Kinetic Energy of Wind-Turbine Generators for System Frequency Support. IEEE Trans. Power Syst. 2009, 24, 279–287. [Google Scholar] [CrossRef]
  39. Yang, D.; Kim, J.; Kang, Y.C.; Muljadi, E.; Zhang, N.; Hong, J.; Song, S.H.; Zheng, T. Temporary Frequency Support of a DFIG for High Wind Power Penetration. IEEE Trans. Power Syst. 2018, 33, 3428–3437. [Google Scholar] [CrossRef]
  40. Ali, Y.A.; Ouassaid, M. Advanced Control Strategy of DFIG based Wind Turbine using combined Artificial Neural Network and PSO Algorithm. In Proceedings of the 2020 International Conference on Electrical and Information Technologies (ICEIT), Rabat, Morocco, 4–7 March 2020; pp. 1–7. [Google Scholar] [CrossRef]
  41. Blaabjerg, F.; Xu, D.; Chen, W.; Zhu, N. Advanced Control of Doubly Fed Induction Generator for Wind Power Systems; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2018. [Google Scholar] [CrossRef]
  42. Ali, Y.A.; Ouassaid, M.; Muljadi, E. Reduced Switching Frequency Finite Control Set Model Predictive Control (FCS-MPC) for DFIG. In Proceedings of the 2021 Innovations in Power and Advanced Computing Technologies (i-PACT), Kuala Lumpur, Malaysia, 27–29 November 2021; pp. 1–7. [Google Scholar] [CrossRef]
  43. Yaramasu, V.; Wu, B. Model Predictive Control of Wind Energy Conversion Systems, 1st ed.; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar] [CrossRef]
  44. Ekanayake, J.; Jenkins, N. Comparison of the Response of Doubly Fed and Fixed-Speed Induction Generator Wind Turbines to Changes in Network Frequency. IEEE Trans. Energy Convers. 2004, 19, 800–802. [Google Scholar] [CrossRef]
  45. Mauricio, J.; Marano, A.; Gomez-Exposito, A.; Martinez Ramos, J. Frequency Regulation Contribution Through Variable-Speed Wind Energy Conversion Systems. IEEE Trans. Power Syst. 2009, 24, 173–180. [Google Scholar] [CrossRef]
  46. Morren, J.; de Haan, S.; Kling, W.; Ferreira, J. Wind Turbines Emulating Inertia and Supporting Primary Frequency Control. IEEE Trans. Power Syst. 2006, 21, 433–434. [Google Scholar] [CrossRef]
  47. Model Predictive Control System Design and Implementation Using MATLAB®; Advances in Industrial Control; Springer: London, UK, 2009. [CrossRef]
  48. Kundur, P. Power System Stability and Control; McGraw Hill: New York, NY, USA, 1994. [Google Scholar]
  49. Vahidi, B.; Tavakoli, M.; Gawlik, W. Determining Parameters of Turbine’s Model Using Heat Balance Data of Steam Power Unit for Educational Purposes. IEEE Trans. Power Syst. 2007, 22, 1547–1553. [Google Scholar] [CrossRef]
  50. Report, IEEE Committee. Excitation System Models for Power System Stability Studies. IEEE Trans. Power Appar. Syst. 1981, PAS-100, 494–509. [Google Scholar] [CrossRef]
  51. Almabrok, A.; Psarakis, M.; Dounis, A. Fast Tuning of the PID Controller in An HVAC System Using the Big Bang–Big Crunch Algorithm and FPGA Technology. Algorithms 2018, 11, 146. [Google Scholar] [CrossRef]
Figure 1. WT power characteristics with β = 0° and at different wind speeds; bold red (MPPT curve), bold blue and black (Shifted optimal power point tracking (OPPT) curves).
Figure 1. WT power characteristics with β = 0° and at different wind speeds; bold red (MPPT curve), bold blue and black (Shifted optimal power point tracking (OPPT) curves).
Energies 16 04389 g001
Figure 3. Structure of the three-phase voltage source inverter.
Figure 3. Structure of the three-phase voltage source inverter.
Energies 16 04389 g003
Figure 4. Switching pattern of the three-phase voltage source inverter.
Figure 4. Switching pattern of the three-phase voltage source inverter.
Energies 16 04389 g004
Figure 5. DFIG-based WT overall control scheme.
Figure 5. DFIG-based WT overall control scheme.
Energies 16 04389 g005
Figure 6. Inertial controller. Upper section: MPPT control. Middle section: frequency droop control. Lower section: inertia emulation.
Figure 6. Inertial controller. Upper section: MPPT control. Middle section: frequency droop control. Lower section: inertia emulation.
Energies 16 04389 g006
Figure 7. Kundur’s AC grid simulation system.
Figure 7. Kundur’s AC grid simulation system.
Energies 16 04389 g007
Figure 8. IEEE-G tandem compound steam turbine governor model.
Figure 8. IEEE-G tandem compound steam turbine governor model.
Energies 16 04389 g008
Figure 9. Conventional and sectional droop control. Scheme #1 [26] and Scheme #2 [28].
Figure 9. Conventional and sectional droop control. Scheme #1 [26] and Scheme #2 [28].
Energies 16 04389 g009
Figure 10. Simulation results for case 1. (a) Frequency. (b) Frequency deviation. (c) Additional power. (d) Rotor speed.
Figure 10. Simulation results for case 1. (a) Frequency. (b) Frequency deviation. (c) Additional power. (d) Rotor speed.
Energies 16 04389 g010
Figure 11. Simulation results for case 1. (a) WTS power. (b) WTS power–speed trajectory. (c) Comparison results for frequency deviation. (d) Comparison results for settling time.
Figure 11. Simulation results for case 1. (a) WTS power. (b) WTS power–speed trajectory. (c) Comparison results for frequency deviation. (d) Comparison results for settling time.
Energies 16 04389 g011
Figure 12. Simulation results for case 2. (a) Wind variation. (b) Frequency. (c) ISE. (d) ITSE. (e) IAE. (f) ITAE.
Figure 12. Simulation results for case 2. (a) Wind variation. (b) Frequency. (c) ISE. (d) ITSE. (e) IAE. (f) ITAE.
Energies 16 04389 g012
Table 1. Control strategies for primary frequency control comparison.
Table 1. Control strategies for primary frequency control comparison.
ControlAdvantagesDisadvantagesSFD
Consideration
Ref
Inertial control
-
Easy to implement.
-
Low computation rate.
-
Dependence of the ROCOF and K i n .
-
Slow response.
No[22]
[23]
Droop control
-
Easy to implement.
-
Low computation rate.
-
Dependence of the Droop and K D .
-
Slow response.
No[7]
[24]
[25]
[26]
Combined control & Variable Droop
-
Easy to implement.
-
Low computation rate.
-
Dependence of the ROCOF and Droop.
-
Tuning of K D and K i n is mandatory.
No[12]
[27]
[28]
[29]
PID-Fuzzy
-
Relatively easy to implement.
-
Moderated computation rate.
-
Tuning of K P and K I is mandatory.
-
Hard definition of control rules.
No[30]
[31]
Fuzzy Logic
-
Relatively easy to implement.
-
Moderated computation rate.
-
Tuning of K P and K I is mandatory.
-
Hard definition of control rules.
No[17]
[32]
[33]
[34]
MPC Control
-
The consideration of control constraints.
-
The output is optimally calculated.
-
Advanced level of control difficulty.
-
Heavy computations are involved.
YesProposed
Table 2. Parameters of The IEEE-G tandem compound steam turbine governor model.
Table 2. Parameters of The IEEE-G tandem compound steam turbine governor model.
K V H P K H P K I P K L P
0 0.36 0.36 0.28
T S C ( s ) T R H 1 ( s ) T R H 2 ( s ) T C O ( s )
0.5 3.3 101
Table 3. Obtained results for case 1.
Table 3. Obtained results for case 1.
ParameterController Results
Proposed 49.8503
FrequencyScheme #2 49.8401
Nadir (Hz)Scheme #1 49.8367
MPPT 49.8324
MaximumProposed −0.1497
frequencyScheme #2 −0.1599
deviationScheme #1 −0.1633
(Hz)MPPT −0.1676
SettlingProposed 31.7937
timeScheme #2 42.3113
(s)Scheme #1 43.5954
SettlingProposed 49.8761
frequencyScheme #2 49.8836
(Hz)Scheme #1 49.8836
MPPT 49.8831
Cluster 10.7088
ProposedCluster 20.7539
Converged Cluster 30.9690
rotor speed (p.u) Cluster 10.8001
Scheme #1 & #2Cluster 20.9002
Cluster 31.1002
Table 4. Comparison of performance index for case 2.
Table 4. Comparison of performance index for case 2.
Control SchemeISEITSEIAEITAE
Proposed0.00280.02020.30226.7177
Scheme #10.00690.05020.41808.1983
Scheme #20.00700.05420.43228.2305
MPPT operation0.00950.05670.45319.2788
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ali, Y.A.; Ouassaid, M.; Cabrane, Z.; Lee, S.-H. Enhanced Primary Frequency Control Using Model Predictive Control in Large-Islanded Power Grids with High Penetration of DFIG-Based Wind Farm. Energies 2023, 16, 4389. https://doi.org/10.3390/en16114389

AMA Style

Ali YA, Ouassaid M, Cabrane Z, Lee S-H. Enhanced Primary Frequency Control Using Model Predictive Control in Large-Islanded Power Grids with High Penetration of DFIG-Based Wind Farm. Energies. 2023; 16(11):4389. https://doi.org/10.3390/en16114389

Chicago/Turabian Style

Ali, Youssef Ait, Mohammed Ouassaid, Zineb Cabrane, and Soo-Hyoung Lee. 2023. "Enhanced Primary Frequency Control Using Model Predictive Control in Large-Islanded Power Grids with High Penetration of DFIG-Based Wind Farm" Energies 16, no. 11: 4389. https://doi.org/10.3390/en16114389

APA Style

Ali, Y. A., Ouassaid, M., Cabrane, Z., & Lee, S. -H. (2023). Enhanced Primary Frequency Control Using Model Predictive Control in Large-Islanded Power Grids with High Penetration of DFIG-Based Wind Farm. Energies, 16(11), 4389. https://doi.org/10.3390/en16114389

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop