Next Article in Journal
Parasitic Effects of PWM-VSI Control Leading to Torque Harmonics in AC Drives
Next Article in Special Issue
A Straightforward Approach to Site-Wide Assessment of Wind Turbine Tower Lifetime Extension Potential
Previous Article in Journal
Increasing Parameters of Diesel Engines by Their Transformation for Methanol Conversion Products
Previous Article in Special Issue
A Stochastic Petri Net Model for O&M Planning of Floating Offshore Wind Turbines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Control of a PMSG-Based Wind Turbine Generator Using Lyapunov Function

by
Roghayyeh Pourebrahim
1,
Amin Mohammadpour Shotorbani
1,2,
Fausto Pedro García Márquez
3,*,
Sajjad Tohidi
1 and
Behnam Mohammadi-Ivatloo
1
1
Department of Electrical and Computer Engineering, University of Tabriz, Tabriz 51666-16471, Iran
2
Faculty of Applied Science, The University of British Columbia, Vancouver, BC V6T 1Z4, Canada
3
Ingenium Research Group, Universidad Castilla-La Mancha, 13071 Ciudad Real, Spain
*
Author to whom correspondence should be addressed.
Energies 2021, 14(6), 1712; https://doi.org/10.3390/en14061712
Submission received: 9 February 2021 / Revised: 8 March 2021 / Accepted: 11 March 2021 / Published: 19 March 2021
(This article belongs to the Special Issue Lifetime Extension of Wind Turbines and Wind Farms)

Abstract

:
This paper proposes a robust finite-time controller (FTC) for a permanent magnet synchronous generator (PMSG)-based wind turbine generator (WTG). An adaptive observer is used for the rotor angle, rotor speed, and turbine torque estimations of the PMSG, thus eliminating the use of anemometers. The robustness of the proposed FTC regarding parameter uncertainty and the external weak power grid is analyzed. The impacts of the power grid short-circuit ratio (SCR) at the point of common coupling (PCC) on the conventional proportional-integral (PI) controller and the proposed FTC are discussed. Case studies illustrate that the proposed observer-based FTC is able to estimate the mechanical variables accurately and provides robust control for WTGs with parameter uncertainty and weak power grids.

1. Introduction

The control and integration of wind energy into a power grid are challenging due to geographical constraints and the intermittency of energy production. The maximum power point tracking (MPPT) control plays a vital role in maximizing the wind energy yield from wind turbines [1]. Wind-characteristic-based MPPT methods require online wind speed and mechanical torque measurements using mechanical sensors, e.g., an anemometer installed at a nacelle [2]. Sensorless controls of wind turbine generators (WTGs) can eliminate expensive mechanical sensors and improve the reliability of the WTGs [3,4,5]. Sensorless controls can also facilitate sensor fault detection and isolation when mechanical sensors are present [6].
Kalman filters can be used for rotor position and speed observations of WTGs [4]. However, Kalman filters are based on linearized models of the nonlinear WTG system [7]; thus, they are not suitable for WTGs at low speeds with small back electromotive forces (EMFs). Rotor position observation through high-frequency signal injection is also challenging at high speeds [8]. Artificial neural network (ANN)-based sensorless controls are proposed in [9,10], which requires extensive offline data collection and training to achieve nonlinear input–output mapping. In [11], the rotor position is estimated using a phase-locked loop (PLL). Back EMF-based PLL observers have been studied for sensorless control of permanent magnet synchronous motors (PMSMs) to estimate the position/speed of the rotor [12,13]. A PLL combined with synchronous frequency-extract filters was proposed in [12] to mitigate the estimated back EMF harmonic error in the sensorless control of surface PMSMs. In [13], a new finite-position set–PLL (FPS-PLL) based on model predictive control was investigated, where the rotor position was identified from the finite set of rotor positions. The consequences demonstrated excellent performance in estimating the position/speed of the rotor in the FPS-PLL compared to the conventional PLL, as well as a significant reduction in dynamic and steady-state fluctuations. Two computationally efficient finite-position set–PLL were proposed in [14] to diminish the computational effort of the FPS-PLL proposed in [13].
Sliding-mode observers for PMSM sensorless controls have been extensively surveyed [15,16,17,18,19,20]. A wind speed estimation and sensorless control for a surface-mounted permanent magnet synchronous generator (PMSG) using sliding mode control (SMC) based on a linear matrix inequality (LMI) was proposed in [15], in which the rotor speed is estimated according to the rate of change in the estimated rotor position. In [19] and [20], the estimations of the rotor position and the stator flux linkage were investigated using a quasi-sliding mode observer and a model reference adaptive system using an sliding mode observer (SMO). However, the turbine mechanical torque was not estimated in [11,19,20]. Additionally, the maximum back EMF and the sampling frequency impose constraints on the convergence of the quasi-sliding mode observer [19]. In [21,22], second-order sliding mode (SOSM) control for a doubly fed induction-generator-based wind turbine has been discussed. This control method improves the classical sliding mode (first-order) chattering problem. Moreover, the main features of SOSM control are robustness with respect to perturbations (grid fault) and chattering-free behavior. An adaptive Kalman-like observer was proposed in [23], which can provide accurate and reliable estimations of system states and parameters for a generic nonlinear system. It will be used in the paper for estimations of the rotor speed and the turbine mechanical torque of the WTG system under a wide range of operating conditions.
A linear controller for the WTG is the most widely used control method, e.g., proportional-integral (PI) controllers [24,25,26,27], the fuzzy propotional-integral-derivative (PID) [28], and the self-tuned adaptive PI controller [29]. PI controllers are very common in the industry due to their simple structure and easy implementation in various applications. Such controllers are designed for a single operating point. The operating point changes due to the disturbances and variation in parameters. In this case, the controller will not have a proper performance, and voltage and power oscillations will appear. In addition, such controllers are very sensitive to parameter variations. Hence, an accurate mathematical system model is required to determine the control parameters. However, accurate calculation of the system model is not possible due to unknown machine and environment parameters, core saturation, wind oscillations, temperature variations, and system disturbances [8,24,25].
A linear controller is fundamental to the stable and reliable operation of the WTG system. Connecting the WTG system to a weak power grid represents challenges and requires detailed stability analysis of the WTG and the power grid [30]. The inherent right half-plane (RHP) zero and pole in the converter’s transfer function are major concerns for the linear controller design of the PMSG-based WTG system when connected to a weak grid [31,32]. Accordingly, root-locus analysis shows that a pole pair enters the RHP and becomes unstable by decreasing the short-circuit ratio (SCR) of the grid [33]. The PLL dynamic also plays a critical role in connecting WTGs to a weak grid [34]. In [35], the power synchronization method was developed to connect an off-shore wind farm to a weak power grid. A three-level hierarchical control scheme was used in [36] to support the point of common coupling (PCC) voltage for connecting the WTG to a weak grid.
Linear controllers are designed under specific operating conditions and are influenced by large external perturbations and parameter uncertainties. Therefore, several nonlinear control techniques are proposed, including the intelligent fuzzy sliding-mode control [37], radial basis function network-based neural network control [38], and adaptive fuzzy control [39]. However, these artificial intelligence control methods require prior behavioral knowledge about the WTG system and extensive training data. Although the nonlinear controllers in [37,38,39,40] led to exponential convergence of the state trajectories, finite-time controller (FTC) architectures [41] can push the control system error trajectories to zero in a pre-defined time, thus resulting in fast convergence of system states, high control accuracy in the steady state, and excellent robustness against perturbations and uncertainties.
This paper proposes a robust FTC method for the PMSG-based WTG system. The nonlinear controls of the machine-side converter (MSC) and the grid-side converter (GSC) are based on the control Lyapunov function (CLF) and back-stepping control design techniques. The main contributions of the paper are summarized as follows:
  • A novel back-stepping finite-time controller is proposed for the PMSG-based wind generation system, and its stability and finite-time convergence are proved mathematically using the Lyapunov stability theorem.
  • The proposed multi-loop output-feedback nonlinear FTC improves the control robustness against parameter uncertainty by the proper tuning of its control gains. The robustness of the FTC is verified using the perturbed WTG system with 20% parameter variation and a weak power grid at the PCC.
  • The proposed FTC leads to fast convergence of the system states and small steady-state errors for normal and perturbed WTG systems with parameter uncertainties. Therefore, the FTC achieves faster maximum power point tracking (MPPT) and extracts more wind power compared to conventional linear and nonlinear controllers.
  • An adaptive Kalman-like observer is proposed for the WTG system to estimate the rotor position, rotor speed, and turbine mechanical torque. Thus, the proposed FTC achieves mechanical sensorless control of the WTG and increases the reliability of the WTG system.
  • Using the reactive power control loop of the proposed controller, the grid’s voltage stability in a weak grid with a low circuit ratio is investigated.

2. PMSG-Based Wind Turbine System Model

The PMSG-based WTG connected to the AC power grid at the PCC, i.e., E n , is shown in Figure 1. The PMSG is connected to the grid through a voltage source converter (VSC) consisting of an MSC and a GSC. The AC power grid is modeled as a voltage source ( E g ) behind a series branch ( r n , L n ). The mechanical power extracted from a wind turbine is calculated as
P m = 1 2 c p ρ A v w 3
where c p = c p ( λ ) , λ = R ω m / v w , R, ω m , ρ , v w and A = π R 2 are the turbine power coefficient, tip speed ratio, turbine blade, radius, turbine shaft speed, air density, wind speed, and rotor swept area, respectively.
The state–space model of the PMSG-based WTG system in the qd rotor reference frame is expressed as [40]
{ d ω r d t = 1 J ( F ω r + K M i q s T m ) d i q s d t = 1 L s ( r s i q s + L s p ω r i d s K M ω r + v q s ) d i d s d t = 1 L s ( r s i d s L s p ω r i d s + v d s ) d v d c 2 d t = 1 C ( E d n i d n + E q n i q n v d c i R ) d i d n d t = 1 L 0 ( E d n L 0 ω n i q n v d n ) d i q n d t = 1 L 0 ( E q n + L 0 ω n i d n v q n )
where T m is the turbine mechanical torque, K M = p 3 / 2 φ r is the flux generator constant, ϕr is the flux generator constant, and ir is the output DC current of the MSC. The variables and quantities of the WTG system used in the paper are listed in Table 1 [27].
The WTG system can be separated into two inter-connected subsystems according to Equation (2) as the WT-PMSG-MSC and the GSC-PCC grid. Each subsystem consists of two control loops, as depicted in Figure 1. The reference rotor speed ω r *   is determined by the MPPT block based on the PMSG stator output active power P s . The MSC’s d axis current reference i d s * is set to zero ( i d s * = 0 ) to attain the maximum torque per ampere [3]. The DC-link voltage reference v d c * is set to be constant to stabilize converter operation. The PCC reactive power reference Q * is be determined based on the PCC voltage regulation or constant power factor according to IEEE Standard 1547. The observer block in Figure 1 estimates the rotor speed, mechanical torque, and flux linkage using stator voltages and currents in the αβ stationary reference frame. The rotor position estimation θ ^ is calculated based on the estimated flux linkage and is used for a b c d q transformation. These nonlinear regulators and the observer blocks in Figure 1 will be introduced in Section 3, Section 4, Section 5 and Section 6.

3. Linear Control Scheme for the WTG System

The control loops for the MSC and the GSC using the PI control scheme in a cascaded structure is given in Figure 2. The outer control loops determine the dq axis current reference of each VSC, and the inner control loops regulate the dq axis currents of the MSC and the GSC to the references. The inner loops in Figure 2 are usually designed to be approximately 2–5 times faster than the outer loop to track the references in a suitable time. The PI controllers are tuned considering the stability and performance of the linearized model of the system at the operating points and the desired loop bandwidth. Various on-site PI controller tuning rules are proposed in the literature, e.g., the Cohen–Coon, Ziegler–Nichols, Tyreus–Luyben, Ciancone–Marline, internal model control, and auto-tune methods, with different objectives.
Tuning the PI controller using linear control theory requires a linear plant model. For stability considerations, the real parts of the poles of the plant at the operating point must have negative values. Further, increasing the loop bandwidth enhances the performance and provides a faster response. However, the loop bandwidth and response speed of the plant with a low-order linear controller is limited by control stability considerations. An appropriate gain and phase margin (i.e., typically 60 ° ) of the loop must be considered for adequate stability margins and as a trade-off between the performance and robustness of the plant. Since linear controllers are tuned based on specific operating points, the performance of the linear controller can be affected by the power grid strength at the PCC, e.g., a weak grid, and parameter uncertainty of the WTG.
Next, the analysis of the linearized model of the GSC at the PCC and the influence of a weak grid on linear PI controllers of the GSC are discussed. In a weak grid, the three phase voltages at the PCC are affected by the power flow generated by the WTG. Considering this impact on the linearized model of the feed-forward power control loop (FFPCL) [31], two main stability concerns can be identified, i.e., the existence of an RHP zero and an RHP pole in rectifier and inverter modes of operation, respectively [29,30,31]. The RHP zero and pole are presented briefly next.
Consider the transfer function between the PCC active ( Δ P n ) and VSC output voltage ( Δ v n ) which is part of the forward path of the FFPCL and is calculated as [31]
Δ P n 1 ( L o + L n ) s + R n [ ( L n i ¯ d n s + E ¯ d n ) Δ v n + ( L o i ¯ d n s E ¯ d n ) Δ E d g + ( ω n ( L o + L n ) E ¯ d n ) Δ i q n ]
where L n , R n , i ¯ d n , E ¯ d n , and s are grid inductance and resistance, the current and voltage of the PCC at the operating point, and the Laplace operator, respectively. As seen in Equation (3), an RHP zero exists in the FFPCL in a rectifier mode of operation (i.e., i ¯ d n < 0 ). This RHP zero causes transient undershoot, limits the use of high control gain due to instability, and restricts the closed-loop bandwidth.
Besides, the DC-link-linearized dynamic in the FFPCL is
Δ v d c v ¯ d c C s v ¯ d c 2 P ¯ n ( Δ P n + v ¯ d c Δ i R )
where v ¯ d c and P ¯ n are the DC-link voltage and the PCC active power at the operating point C and i R are the DC-link capacitance and the output DC current of the MSC shown in Figure 1. Considering the denominator of Equation (4), an RHP exists in the inverter mode of operation ( P ¯ n > 0 ), which limits the control gain due to instability [31].
Consequently, designing a linear PI controller based on the linearized model of the FFPCL would be insufficient for large variations in the operating point, which is further confirmed through root-locus analysis in [33]. Therefore, a nonlinear control law [42] with global stability will be desirable since it is independent of the operating point and compensates for the variation in operating conditions and grid strength.

4. Finite-Time Control Design for the PMSG-Based Wind Turbine Generator

In this section, the finite-time stability [41] and back-stepping control design [42] for a nonlinear dynamic system are reviewed first. Then, FTC laws for the MSC and the GSC are designed based on the CLF and back-stepping design techniques.

4.1. FTC and Back-Stepping Control Design

Assume a continuous function f : D n on D containing the origin (i.e., 0 D ), and a nonlinear system as
x ˙ = f ( x ) ,   f ( 0 ) = 0 ,   x n
Suppose V ( x , t ) : D is a Lyapunov function. Then, the origin of Equation (5) is a local finite-time stable equilibrium if
t V + ρ V α 0
where ρ > 0 and 0 < α < 1 are real numbers.
Moreover, the origin is a globally finite-time stable equilibrium if D = n and V is radially unbounded. The finite convergence time of the trajectories, i.e., T ( x 0 ) , satisfies [34].
  T ( x 0 ) V 1 α ( x 0 , 0 ) ρ ( 1 α )
where x 0 = x ( t 0 ) is the initial state.
The back-stepping control design technique is often used for a nonlinear system in the form of the following:
{ χ ˙ = h ( χ ) + g ( χ ) ξ   ξ ˙ = ν ( χ ) + u ( χ , ξ )
where χ = [ χ i ]   and ξ = [ ξ j ] are system state vectors; u is the control vector; and h , g , and ν are smooth vector functions.
The back-stepping design achieves the control of state χ through the control of state ξ since there is no direct control input u available for χ , as shown in Equation (8). Thus, by controlling ξ , the control input u indirectly controls the state variable χ in a back-stepping fashion.
Proposition 1.
Assume that the proposed virtual controller ξ = ξ F T * ( χ ) is a finite-time stabilizing control law for χ ˙ in Equation (8) and satisfies Equation (6) with a Lyapunov function. Then, the proposed control law (i.e., Equation (9)) stabilizes the origin of the system (i.e., Equation (8)) in finite time:
u F T ( χ , ξ F T * ) = ξ F T * χ ( h + g ξ ) V a χ g k ( ξ ξ F T * ) ν k t α ( ξ ξ F T * )
where k and k t are positive real diagonal matrices, i.e., k = d i a g ( k i ) and k t = d i a g ( k t i ) , and is the finite-time fractional term proposed as ( z ) = | sinh ( z ) | α s i g n ( z ) , z n ,   0 < α < 1 ,   α .
In the proposed controller, ξ F T * is designed as
ξ F T * ( χ ) = g 1 ( h K χ K t sign ( χ ) )
Proof. 
First, we show the finite-time stability of system χ ˙ with the proposed controller ξ F T * ( χ ) and then prove the finite-time stability of the system (i.e., Equation (8)) with the proposed back-stepping controller (i.e., Equation (9)).
The finite-time Lyapunov stability of system χ ˙ can be verified using the candidate control Lyapunov function (CLF). By substituting ξ F T * ( χ ) to the derivative of the CLF as V ˙ a ( χ ) = V a χ χ ˙ = χ T χ ˙ , we have
V ˙ a ( χ ) = χ T ( h + g ξ   ) = χ T ( K χ K t α ( χ )   ) = K χ T χ K t χ T | sinh ( χ ) | α sign ( χ )
Considering that | z | sinh ( | z | ) e | z | 1 for z , and thus | sinh ( χ ) | α | χ | α , one concludes
V ˙ a ( χ ) K χ T χ K t | χ T | | χ | α k i χ i 2 k t i | χ i | α + 1 V ˙ a 2 k V a 2 β k t V a β
where 0 < β = α + 1 2 < 1 is the fractional power and k and k t are the minimum eigenvalues of the gain matrices K and K t , respectively. The inequality V ˙ a ( χ ) 2 β k t V a β shows the finite-time stability of system χ ˙ with the proposed controller ξ F T * ( χ ) with respect to Equation (7).
Moreover, the finite-time stability of Equation (8) using the proposed controller (Equation (9)) can be proved by using the CLF V b = 1 2 ( ξ ξ F T * ) T ( ξ ξ F T * ) . The detailed proof and derivation of the proposed control laws are given in Appendix A. □
The FTC-based regulators for the MSC and GSC of the WTG are proposed using the back-stepping control design (i.e., Equation (9)). It is noticed that the dynamic equations of ω r and v d c 2 in the WTG model (i.e., Equation (2)) are not directly regulated by any control inputs. Thus, the back-stepping technique is used to design the control law. The q axis current of the MSC and the d axis current of the GSC are used to regulate ω r and v d c 2 . Comparing the state space model (i.e., Equation (2)) with Equation (8), we obtain
{ χ = [ ω r , v d c 2 ] T ξ = [ K M i q s / J , E d n i d n / C ] T u = [ v q s / L s , v d n / L 0 ]

4.2. Rotor Speed and Stator Current FTC Design for the MSC

The reference rotor speed ω r * is determined through the MPPT block in Figure 1. Based on the back-stepping design in Equation (9) and the virtual controller ξ 1 = K M i q s / J in Equation (10), the control inputs v q s * and v d s * are proposed as
ξ 1 * = 1 J ( F ω r T m ) + ω ˙ r * k 1 e 1 k t 1 sign ( e 1 )
v q s * = ( r s i q s + L s p ω r i d s K M ω r ) + J L s K M ( ξ ˙ 1 * e 1 k 2 e 2 k t 2 α ( e 2 ) )
v d s * = L s ( k 3 e 3 r s L s e 3 + p ω r i q s + k t 3 α ( e 3 ) )
where e 1 = ω r ω r * , e 2 = ξ 1 ξ 1 * , e 3 = i d s i d s * ; k t i and k i are real positive design constants for i = 1 ,   2 ,   3 ; and α ( e i ) = | sinh ( e i ) | α sign ( e i ) . Moreover, with respect to Equation (11), ξ ˙ 1 * is calculated in Appendix A.
The nonlinear control law (Equations (11)–(13)) realizes the ω r and i d s regulators of the MSC in Figure 1. The detailed derivation of Equations (11)–(13) is given in Appendix A.

4.3. DC-Link Voltage and Reactive Power FTC for the GSC

The DC-link voltage of the MSC needs to be regulated exactly to its constant reference value ( v d c * ) for stable operation of the WTG. The reactive power injected by the GSC is used to control the PCC voltage. Therefore, the dq currents of the GSC are regulated to achieve the DC voltage and the injected reactive power controls. Based on the back-stepping design in Equation (9) and the virtual controller ξ 2 = E d n i d n / C in Equation (10), the FTCs v d n * and v q n * are designed as
ξ 2 * = σ 1 k 4 e 4 k t 4 sign ( e 4 )
v d n * = C L 0 E d n ( σ 2 + e 4 + k 5 e 5 + k t 5 α ( e 5 ) )
v q n * = L 0 E d n ( ω n E q n i q n + Q ˙ n * ) + ω n L 0 i d n + E q n E d n v d n * L 0 E d n ( k 6 e 6 + k t 6 α ( e 6 ) )
where e 4 = v d c 2 ( v d c * ) 2 , e 5 = ξ 2 ξ 2 * , and e 6 = Q n Q n * . In Equations (14)–(16), k t and k i are real positive control design parameters, s i = s i g n ( e i ) for i = 4 ,   5 ,   6 , and σ 1 and σ 2 are
σ 1 = 1 C E q n i q n + 1 C v d c i R
σ 2 = k 4 e ˙ 4 + σ ˙ 1 + 1 C L 0 E d n 2 1 C ω n E d n i q n
The nonlinear control law (Equations (14)–(16)) realizes the v d c and Q regulators of the GSC in Figure 1. The detailed derivation of Equations (14)–(16) is not included in the paper for space considerations. It follows a similar derivation procedure as for the MSC. The detailed demonstration of the regulators in Figure 1 for the MSC and the GSC using the proposed nonlinear FTC are depicted in Figure 3. The variables in this diagram are defined in Equations (11)–(16).

4.4. Global Stability Proof of the Proposed FTC for the PMSG

Finite-time stability of the WTG system (Equation (2)) with the proposed FTC is proved using Lyapunov stability criteria. Given the Lyapunov function V t = 1 2 i = 1 6 e i 2 , we have
V ˙ t = i = 1 6 e i e ˙ i
where e i , i = 1 , , 6 are the errors of the states and have been defined previously. The error dynamics, i.e., e ˙ i , can be calculated using system dynamics (Equation (2)). Substituting the proposed FTC from Equations (11)–(16) into e ˙ i in Equation (19), we conclude
V ˙ t = i = 1 6 k i e i 2 i = 1 6 k t | e i | α + 1
Therefore, the WTG system is Lyapunov-stable since V ˙ t in Equation (20) is negative definite for positive errors.
In addition, considering i = 1 6 k i e i 2 0 in (20), we have
V ˙ t k t i = 1 6 | e i | α + 1 2 β k t ( 1 2 i = 1 6 e i 2 ) β
This indicates V ˙ t ρ m V t α with ρ m = 2 β k t and β = ( α + 1 ) / 2 , which proves the finite-time stability of Equation (2) with the proposed FTC law (Equations (11)–(16)) in accordance with Equation (6).
If the coefficients k t of the sign functions are set to zero, the control law (Equations (11)–(16)) becomes the conventional nonlinear exponential convergent controller (ECC) and does not satisfy finite-time convergence in Equation (6) anymore. The comparison of the conventional nonlinear ECC and the proposed FTC is performed in Section 7 for the WTG system.

5. Robustness of and Chattering in the Proposed FTC

Uncertainties due to modeling errors and non-ideal operating conditions are inevitable in practice. Therefore, a robust design of a control system against uncertainties such as parameter variations is important for practical applications. Robustness enhancement of a closed-loop system by the proposed FTC is investigated in this section.

5.1. Robustness against Parameter Uncertainty

To analyze the effect of perturbations on system stability, assume the parameter vector θ in Equation (2) as
θ = [ J , F , K M , T m ,   r s ,   L s ,   E q n ,   E d n , C ,   L 0 ] T
An uncertain parameter is formulated as
θ j + Δ θ j = ( 1 + ξ θ j ) θ j ,   0 ξ θ j 1
where θ j , Δ θ j = ξ θ j θ j , and ξ θ j are the nominal parameter, parameter variation, and the uncertainty factor, respectively.
Correspondingly, the dynamics of the perturbed system (Equation (2)) can be separated into a nominal part ( F N = [ f N i ] ,   i = 1 6 ) and a perturbed part ( F δ = [ f δ i ] ,   i = 1 6 ) as [35]
e ˙ i = f N i ( t ) + f δ i ( t )
It is note that the nominal part is identical to Equation (2), while the perturbed part is a time-variant function of the system states and nominal parameters, i.e., f δ i = f δ i ( θ , x , t ) .
Suppose the WTG system model (Equation (2)) is perturbed by parameter uncertainties in accordance with Equation (24), the perturbed part F δ is norm-bounded [42] as
F δ a + b e
where a and b are known real positive constants.
Taking the time derivative of the Lyapunov function V t = 1 2 i = 1 6 e i 2 for the perturbed system considering Equation (24), we have
V t t = i = 1 6 e i e ˙ i = i = 1 6 e i ( f N i + f δ i )
According to Equation (20), the time derivative of the Lyapunov function V t for the nominal system can be represented as
i = 1 6 e i f N i = i = 1 6 k i e i 2 i = 1 6 k t | e i | α + 1
Substituting Equations (25) and (27) into Equation (26) gives
V ˙ i = 1 6 ( k i e i 2 k t | e i | α + 1 + e i ( a + b e ) )
By choosing k t a and k i b , we obtain
V ˙ i = 1 6 ( k i b ) e i 2 k t | e i | α + 1 + a e i
which gives V ˙ 0 as long as k i b and a e i k t | e i | α + 1 0 . The first inequality constraint can be applied by choosing control gains k i large enough, k i b ; the second inequality constraint is always true for negative errors, e i 0 . For positive errors, e i 0 , the second inequality constraint requires a k t e i α . Assume that it is chosen to become robust for a desired range of error, i.e., | e i | c > 0 , where c is a real positive constant. Consequently, we must have a k t e i α when | e i | α c , which can be achieved by choosing k t large enough so that it satisfies k t a / c .
According to Equation (6), Equation (29) indicates the robust stability of the closed-loop system with the proposed FTC with k i b and k t a / c , for the error range | e i | α c , in the presence of perturbations. In other words, using sufficiently large control gains k i and k t , the closed-loop system (Equation (2)) with the proposed FTC is robust against the bounded uncertainties. The robustness of the proposed FTC for the WTG system will be verified in Section 7.

5.2. Continuous FTC with Chattering Elimination

The presence of the discontinuous s i g n function in the FTC (Equations (11)–(16)) leads to a chattering problem and produces undesirable oscillatory control outputs. In this paper, the s i g n function is approximated using a continuous inverse tangent function to eliminate chattering as
s i g n ( e ) tanh ( ε e )
where ε is a design-specific real constant.
It is noted that increasing the value of ε leads to better approximation of the s i g n function. However, a very large value of ε may increase the possibility of chattering occurrence. At the same time, ε should be sufficiently large to ensure the robustness of the FTC. When applying the approximation (Equation (30)) in the proposed FTC, the steady-state error of the perturbed system will be bounded (i.e., e ϵ ,   ϵ > 0   ) [42]. The boundary region D : e ϵ is determined by the parameter ε . The closed-loop trajectories of the perturbed system converge toward and enter the region D in a finite time and will remain inside. An appropriate ε is selected through numerical simulations in this paper.

6. State and Parameter Estimations Using an Adaptive Observer

To achieve sensorless control of the PMSG, an adaptive observer [23] is employed to estimate mechanical variables. Since the rotor position is not available, the stationary ( α , β ) reference frame is used for the PMSG model:
{ x ˙ 1 = A 1 x 1 + B 1 + [ 0 1 / J ] T T m x ˙ 2 = A 2 x 2 + B 2
where the matrices are defined as
{ A 1 = [ r s L s p ϕ β r L s 3 p ϕ β r 2 J F J ] , B 1 = [ u α s L s 3 p ϕ α r i β s 2 J ] , x 1 = [ i α s ω r ] A 2 = [ r s L s p ω r L s 0 0 0 p ω r 0 p ω r 0 ] , B 2 = 1 L s [ u β s 0 0 ] , x 2 = [ i β s ϕ α r ϕ β r ]
In Equation (32), u α s , u β s , i α s , i β s , ϕ α r , and ϕ β r are stator voltages, stator currents, and rotor fluxes, respectively, in the stationary reference frame. The adaptive observer and the adaptive law are formulated as Equation (33) and Equation (34), respectively [21]:
{ x ^ ˙ 1 = A 1 ( x ^ 2 ) + B 1 ( x ^ 2 ) + [ 0 1 / J ] T T ^ m + λ x ^ 1 [ Γ Γ T / μ + M 1 1 ] [ e α i 0 ] T x ^ ˙ 2 = A 2 ( ω ^ r ) x ^ 2 + B 2 + M 2 1 [ e β i 0 0 ] T
{ T ^ ˙ m = γ 1 λ e α i / μ μ ˙ = η 1 μ + γ 1 2 Γ ˙ = A 1 Γ + [ λ γ 1 ( M 1 1 ) 1 , 1 1 / J ] T   M ˙ 1 = η 0 M 1 A 1 T ( x ^ 2 ) M 1 M 1 A 1 ( x ^ 2 ) + δ Ι 2 T I 2   M ˙ 2 = η 2 M 2 A 2 T ( ω ^ r ) M 2 M 2 A 2 ( ω ^ r ) + δ Ι 3 T I 3
where η 0 ,   η 1 ,   η 2 ,   δ , and λ are real positive constants; μ ( t ) is a time-variable scalar; M 1 and M 2 are time-variant matrices; ( M 1 1 ) 1 , 1 is the (1,1) element of matrix M 1 1 ; e α i = i α s i ^ α s ; e β i = i β s i ^ β s ; Γ = [ γ 1 γ 2 ] T ; Ι 2 = [ 1 , 0 ] ; and Ι 3 = [ 1 , 0 , 0 ] .
In the dynamic subsystem x ˙ 1 of Equation (31), it is desired to estimate the unknown variable T m and the state vector x 1 simultaneously. Therefore, an adaptive Kalman-like observer is used to estimate x 1 and x 2 by x ^ ˙ 1 and x ^ ˙ 2 . The stability of the observer can be proved using the Lyapunov function [23]:
V o b s = e o 1 T M 1 e o 1 + e o 2 T M 2 e o 2 + μ e o 3 2
with e o 1 = x 1 x ^ 1 Γ e o 3 , e o 2 = x 2 x ^ 2 , and e o 3 = T m T ^ m .
It is assumed that under normal operation of the PMSG wind energy conversion (WEC) system, the data from both wind and torque sensors are used in the proposed controllers rather than the data from the observer. However, the data from the system/parameter observer can be used during sensor failure for a limited period of time, which can improve the reliability of the system. The use of observed data can also be effective if there is a malfunction in the wind and/or torque sensors, such as when the sensor output is out of range. However, although the stability of the observer system is proved by the Lyapunov stability criteria, the delay caused by the observer system has a negative impact on the stability and performance of the system [43], which requires a dedicated investigation.
Applying the chattering elimination in Equation (30) and the adaptive observer in Equation (33), the final output-feedback form of the proposed FTC can be derived from Equations (11)–(16). As an example, the control output v d s * of the proposed FTC is presented as
v d s * = L s [ k 3 e 3 r s L s e 3 + p ω ^ r i q s + 2 π k t 3 α ( e 3 ) + k t | sinh ( e 3 ) | α tanh ( ε e 3 ) ]
The other control outputs of the FTC in Equations (11)–(16) are not presented due to space considerations.

7. Case Study

A 5MW PMSG-based WTG is simulated using Simulink/SimPowerSystems to verify the proposed FTC and observer algorithms. The WTG parameters are listed in Table 1. Two nonlinear controllers, i.e., the ECC and the FTC, are implemented for the WTG. The ECC is derived by setting k t = 0 in the proposed FTC, as described in Section 4. The conventional PI controller is also implemented for the WTG system for benchmark comparison of the control performance to the ECC and FTC. The digital signal processor (DSP) system toolbox is used for the matrix operations of the presented adaptive observer. The initial values of the matrices in Equation (33) are selected as the unity matrix. The parameters of the MSC and GSC of the WTG, the proposed FTC, the adaptive observer, and the PI controller are given in Table 2 and Table 3. Control gains are selected considering the current limits and capacity of the MSC and GSC.
A variant wind profile, illustrated in Figure 4a, is selected for the evaluation of the proposed FTC. The dynamic of the wind profile consists of average value, ramps, gusts, and noise. The mechanical torque of the wind turbine is extracted by the MPPT and rotor speed control and is identical for the FTC and the ECC. The mechanical torque of the wind turbine is shown in Figure 4b.
For the MPPT, the extremum-seeking technique is used since it is independent of the wind turbine’s characteristics. Each harvestable power diagram for a specific wind speed has a peak value, i.e.,   ( ω r , P m a x ) . Optimum rotor speed reference tracking leads to the MPPT of the WTG. In this paper, a fouth-order polynomial of the stator active power P s in pu value is used to represent the turbine’s MPPT reference rotor speed as [40,44]
ω r * ( P s ) = 2.0768 P s 4 + 5.6965 P s 3 5.9254 P s 2 + 3.6040 P s + 0.5516

7.1. Rotor Speed Reference Tracking

Figure 5 shows the comparison of the rotor-speed-tracking performance between the proposed FTC and the ECC. The reference rotor speed for the maximum power extraction is determined by the MPPT algorithm. As shown in Figure 5, the FTC provides faster tracking to the reference rotor speed than the ECC does. The average values of the extracted real power using the FTC and the ECC are listed in Table 4. Two scenarios are compared for different controllers; measurement-based values are calculated using direct rotor speed and WT torque measurements, while observer-based values specify the impact of the adaptive observer transient response. As shown in Table 4, the extracted wind power using the FTC is higher than that using the ECC by 1% since the FTC tracks the optimum rotor speed faster than the ECC does. The relative impact of the transient caused by the adaptive observer on the sensorless MPPT is less than 0.01% compared to the measurement-based MPPT.
To compare the controller outputs of the proposed FTC, the ECC, and the PI controller, an index called control effort (CE) is defined using the mean value of the WTG’s control outputs as
C E = ( 1 τ 0 τ v d s 2 + v q s 2 + v d n 2 + v q n 2 ) 1 2
where τ is the duration of control action.
The average CE for the FTC, the ECC, and the PI controller is 5951, 5948, and 5952, respectively (Figure 5 and Table 4). The proposed FTC extracts more power compared to the ECC and the PI controller though using an identical CE, which relatively differs only by 0.05%.
Figure 6 shows the unextracted wind powers by the FTC and the ECC when the adaptive observer is used. As shown in Figure 6, the FTC presents a smaller unextracted wind power compared to the ECC and improves the MPPT algorithm by accurate and fast tracking of the reference rotor speed.
The GSC of the WTG regulates a constant DC-link voltage under the condition of varying wind power output, as depicted in Figure 7. The reactive power reference of the GSC regulates a unity power factor at the PCC. As shown in Figure 7, the DC-link voltage is effectively regulated by the FTC and the ECC under fluctuating wind power. Figure 8 illustrates the control outputs of the proposed FTC for the case study. It is observed in Figure 8 that the control signals are in the practical range without chattering and saturation.

7.2. Robustness of the Proposed Nonlinear Controls against Uncertainties

The robustness of the closed-loop WTG system using the proposed FTC is investigated under the perturbed condition in which the magnitudes of the WTG parameters increase by 20% from the nominal values. The performance of the proposed FTC is compared to the ECC and the PI controller under nominal and perturbed conditions (Figure 9).
CE values are calculated for the FTC, the ECC, and the PI controller under nominal and perturbed conditions. Figure 10 illustrates that the control effort of the proposed FTC is analogous to the ECC and the PI controller. The proposed FTC provides lower overshoot compared to the PI controller and the ECC, although it uses a similar average CE. Therefore, the FTC is robust against uncertainty and improves the transient response under the perturbed condition as well as the nominal condition. In addition, the Bode diagram of the PI controller’s outer control loops for rotor speed and DC-link voltage in Figure 11 illustrates appropriate tuning of the PI controller with high and sufficient stability margins and robustness.
As discussed in Section 4, the robustness of the proposal FTC is influenced by the control gains k i and k t . Three scenarios are created to analyze the impact of control gains. First, the response of the nominal WTG is evaluated for the FTC with nominal control gains given in Table 3. The simulated rotor speed response is denoted as FTC0 in Figure 12. Second, the same control gains k i and k t in Table 3 are used by the FTC and the ECC for the perturbed system. The resultant responses are denoted as FTC1 and ECC1 in Figure 12. Third, the control gains k i and k t for the FTC and the ECC are multiplied by a factor of 1.5 and denoted by FTC2 and ECC2. As depicted in Figure 12, the settling time and the steady-state error of the rotor speed are reduced by increasing the control gains k i and the finite-time gain k t under the perturbed condition. It is also shown in Figure 12 that the steady-state error and the settling time of the response with the proposed FTC are lower than the ECC’s. The relative tracking error by the proposed FTC is 0.12% for the perturbed system and is zero for the nominal system.
In Figure 13, the impact of increasing the control gains k i and k t on the settling time and the steady-state error of the closed-loop WTG’s response is illustrated for the perturbed system. For the label, FTC k i inc., used in Figure 13, the FTC’s gain k i increases from 1 to 1.5 and 2 pu, while k t is kept constant at 1 pu. Other labels in Figure 13 follow the similar notation as FTC k i inc. It is observed in Figure 13 that the FTC presents a shorter settling time and a smaller steady-state error compared to the ECC for various control gains k t and k i .

7.3. Implemmention of the Adaptive Observer

The implementation of the adaptive observer proposed in Section 5 is investigated for the FTC in this section. The estimation of the mechanical variables ω r and T m by the adaptive observer are shown in Figure 14 for the step change in wind speed. The results illustrate that the proposed observer accurately tracks the reference rotor speed and the mechanical torque in the steady state. As depicted in Figure 14b, the mechanical torque T m estimation has been designed to have a slower transient response in order to suppress any measurement noise presented in T m . The proposed adaptive observer generates an accurate estimation of T m in the steady state, which is essential for the controller stability. It is noted that the estimation error may lead to a minor loss of MPPT during the transient response of the observation.

7.4. Performance of the GSC FTC for a Weak Grid

The robustness of the proposed FTC, ECC, and PI controller was analyzed in Section 5. In this subsection, the impact of the different grid strength and injected reactive power on the variations in the PCC’s voltage is investigated. A step change is applied to the wind speed, which is identical to Figure 9. The performance and resiliency of the proposed FTC are compared to the PI controller for the GSC of the WTG connecting to a weak grid in Figure 15a and Figure 15b, respectively. Commonly, a grid with an SCR lower than 20 is considered a weak grid. Furthermore, the impact of the PCC reactive power control (or GSC q axis current control) is depicted in Figure 15 for a grid with SCRs 5 and 19.
When a step change in the active power of the MSC occurs, the power has to be injected into the grid through the GSC in order to regulate a constant DC-link voltage. Moreover, enabling the reactive power control at the PCC compensates the steady-state PCC voltage variation caused by the grid strength. The DC-link voltage and PCC reactive power are regulated using the GSC d and q axis currents, respectively.
The proposed PI controller is also effective for the weak power grid.
As demonstrated in Figure 15, the steady-state variation and the transient overshoot in the PCC’s voltage magnitude are lower with the proposed FTC (Figure 15a) compared to the conventional PI controller (Figure 15b) since the FTC provides an advanced transient in the GSC dq axis currents and eliminates overshoot.
Similar to the uncompensated case, the proposed FTC (Figure 15a) further improves the transient in the PCC’s voltage magnitude compared to the conventional PI controller (Figure 15b) when the PCC reactive power control is enabled (i.e., Q n * 0 ), as well as with unity power factor (i.e., Q n * = 0 ).

7.5. Discussion

The proposed FTC is developed based on the existing ECC [40] with minor modifications, i.e., adding a finite-time term to each dynamic of the control errors. It is a critical variation, as the system can move from an FTC to an ECC and vice versa easily, yet take the advantages of finite-time convergence. The performance of the proposed FTC is compared to the ECC proposed in [40] and the conventional PI control method. Although the performance of the proposed nonlinear FTC is comparable to nonlinear high-order sliding mode control, the proposed FTC has the advantage that it can be easily implemented by slightly modifying the conventional ECC, i.e., by adding finite-time terms with a fractional power. Since the proposed FTC is developed by adding simple terms to the conventional ECC, its implementation in practice should have the same level of difficulties.

8. Conclusions

This paper proposed a finite-time nonlinear controller for the PMSG-based WTG system. The nonlinear controls of the machine-side converter and the grid-side converter are based on the control Lyapunov function and back-stepping control design techniques. Rotor speed and current are regulated through the machine-side converter for MPPT, whereas the voltage of the DC link is controlled using the grid-side converter’s current regulators. The robustness of the proposed controller was studied for bounded uncertainties. The proposed multi-loop output-feedback nonlinear FTC improves the control robustness against a bounded uncertainty, such as parameter variations, only by the proper tuning of its control gains, as the most effective part in the proposed scheme responsible for robustness against parameter uncertainties. Minimum required control gains to provide robustness against bounded uncertainties are calculated based on the Lyapunov stability theorem. A 5MW WTG is used in the case study to verify the effectiveness of the proposed nonlinear finite-time control techniques. The performance of the proposed FTC is compared with the ECC and the PI controller for the WTG under nominal and perturbed conditions. Simulation studies verify that the proposed FTC achieves fast and accurate MPPT and improves the controller robustness for parameter uncertainties and for a weak power grid compared to the conventional ECC and PI controllers. It also enhances stability of the system. The robustness of the proposed FTC is verified using aperturbed WTG system with 20% parameter variation and a weak power grid at the PCC. Moreover, the outer reactive power control loop is used to control the AC voltage in a weak grid at the PCC with a low circuit ratio. In a weak grid, the proposed FTC improves the steady-state variation and the transient overshoot in the PCC’s voltage compared to the conventional PI controller. As future research work, the authors aim to address the impacts of the delay caused by the observer on the performance and stability of the PMSG WEC system with different control methods, using the Lyapunov–Krasovskii linear matrix inequality (LMI) method for time delay systems.

Author Contributions

Conceptualization, R.P., A.M.S., S.T., B.M.-I. and F.P.G.M.; methodology, R.P. and A.M.S.; software, R.P. and A.M.S.; validation, R.P., A.M.S. and B.M.-I.; resources, S.T., B.M.-I. and F.P.G.M.; writing—original draft preparation, R.P. and A.M.S.; writing—review and editing, R.P., A.M.S. and B.M.-I.; supervision, S.T., B.M.-I.; project administration, S.T., B.M.-I. and F.P.G.M.; funding F.P.G.M. All authors have read and agreed to the published version of the manuscript.

Funding

The work reported herewith was financially supported by the Dirección General de Universidades, Investigación e Innovación of Castilla-La Mancha under the Research Grant ProSeaWind project (SBPLY/19/180501/000102).

Acknowledgments

The authors would like to sincerely offer their special thanks to liwei Wang (School of Engineering, University of British Columbia Okanagan, Kelowna, BC, Canada) for his wholehearted encouragement and support during this research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Design and Stability Proof of the Proposed FTC

The complete control design approach for the WTG-MSC system, i.e., FTC in Equations (11)(13), is presented here. The FTC for the GSC is derived using a similar approach for the MSC as follows. Therefore, the FTC design approach for the GSC is not presented for briefness.
1. Designing ξ 1 * : Assuming a Lyapunov function V 1 = 1 2 e 1 2 , the time derivative of V 1 can be calculated as
V ˙ 1 = e 1 e ˙ 1 = e 1 (   ξ 1 1 J ( F ω r T m ) ω ˙ r * )
By choosing the virtual control law ξ 1 = ξ 1 * defined as Equation (11) and substituting into Equation (A1), we get V ˙ 1 = k 1 e 1 2 k t 1 | e 1 | , which is further simplified as Equation (A2), considering that k 1 e 1 2 0 .
V ˙ 1 k t 1 | e 1 | 2 k t 1 ( 1 2 e 1 2 ) 1 2 2 1 2 k t 1 V 1 1 2
Based on Equations (6) and (A2), ω ˙ r is finite-time-stable with the control law ξ 1 * .
To calculate the dynamics of e ˙ 1 under the proposed control law, we define e 2 = ξ 1 ξ 1 * and substitute it into e ˙ 1 as
e ˙ 1 = 1 J ( F ω r T m ) + ( ξ 1 * + e 2 ) ω ˙ r *
Substituting Equation (11) into Equation (A3) with further simplification gives
e ˙ 1 = k 1 e 1 k t 1 sign ( e 1 ) + e 2
2. Designing v q s * : By supposing a Lyapunov function V 2 = 1 2 i = 1 2 e i 2 , the time derivative of V 2 is calculated as
V ˙ 2 = e 1 e ˙ 1 + e 2 e ˙ 2
Substituting e ˙ 1 from Equation (A3) and e ˙ 2 = ξ ˙ 1 ξ ˙ 1 * into Equation (A5) yields
V ˙ 2 = e 1 ( k 1 e 1 k t 2 α ( e 2 ) + e 2 ) + e 2 ( ξ ˙ 1 ξ ˙ 1 * )
where ξ ˙ 1 and ξ ˙ 1 * are calculated as follows:
ξ ˙ 1 = K M J d i q s d t = K M J L s ( r s i q s + L s P ω r i d s K M ω r + v q s )
ξ ˙ 1 * = 1 J ( F ω ˙ r T ˙ m ) + ω ¨ r * + k 1 e ˙ 1 = 1 J [ F ( 1 J ( F ω r + K M i q s T m ) ) T ˙ m ] + k 1 ( k 1 e 1 + k t 1 sign ( e 1 ) e 2 )
where ω ˙ r and e ˙ 1 are calculated in Equations (2) and (A4), respectively.
Substituting ξ ˙ 1 from Equation (A7) and ξ ˙ 1 * from Equation (A8) into Equation (A6), we design the control law v q s * as Equation (12). Substituting the proposed FTC law v q s * from Equation (12) into V ˙ 2 Equation (A5) yields
V ˙ 2 = e 1 ( k 1 e 1 k t 1 sign ( e 1 ) + e 2 ) + e 2 ( e 1 k 1 e 2 k t 2 α ( e 2 ) )
which is further simplified as Equation (A10), considering that | z | sinh ( | z | ) for z , and thus | sinh ( z ) | α | z | α .
V ˙ 2 k 1 e 1 2 k t 1 | e 1 | k 2 e 2 2 k t 2 | e 2 | α + 1
According to the triangular inequality and (A10), we obtain:
V ˙ 2 2 k ( 1 2 i = 1 2 e i 2 ) 2 β k t 2 ( 1 2 i = 1 2 e i 2 ) β 2 k V 2 2 β k t V 2 β
V ˙ 2 2 β k t V 2 β
where 0 < β = α + 1 2 < 1 is the fractional power and k = min { k 1 ,   k 2 } and k t = min { k t 1 ,   k t 2 } are the control gains.
Based on Equations (6) and (A12), the subsystem [ e ˙ 1 ,   e ˙ 2 ] T is finite-time-stable with the proposed FTC law v q s * in Equation (15). The dynamic of e ˙ 2 with the proposed FTC law in Equation (15) is calculated by substituting ξ ˙ 1 in Equation (11) and v q s * in Equation (12) into e ˙ 2 as
e ˙ 2 = ξ ˙ 1 ξ ˙ 1 * = e 1 k 1 e 2 k t 2 α ( e 2 )
3. Designing v d s * : Defining a Lyapunov function V 3 = 1 2 e 3 2 and taking its time derivative considering i d s * = 0 yields
V ˙ 3 = e 3 e ˙ 3 = e 3 ( 1 L s ( r s i d s L s P ω r i d s + v d s ) )
According to Equation (A14), we define the FTC law v d s * as Equation (13). By substituting v d s * in Equation (13) into Equation (A14) and canceling similar terms, we have
V ˙ 3 k 3 e 3 2 k t 3 | e 3 | α + 1 k t 3 | e 3 | α + 1
According to triangular inequality and Equation (A15), we conclude
V ˙ 3 2 β k t 3 ( 1 2 e 3 2 ) β 2 β k t 3 V 3 β
Based on Equations (6) and (A16), d i d s d t is finite-time-stable with the proposed FTC law v d s * in Equation (13) into e ˙ 3 and some simple calculations as
e ˙ 3 = k 3 e 3 k t 3 α ( e 3 )
The dynamic of the closed-loop WTG system (Equation (2)) with the proposed FTC (Equations (11)–(16)) is represented as
[ e ˙ 1 e ˙ 2 e ˙ 3 ] = [ k 1 1 0 1 k 2 0 0 0 k 3 ] [ e 1 e 2 e 3 ] d i a g ( [ k t 1 k t 2 k t 3 ] ) [ sign ( e 1 ) α ( e 2 ) α ( e 3 ) ]
[ e ˙ 4 e ˙ 5 e ˙ 6 ] = [ k 4 1 0 1 k 5 0 0 0 k 6 ] [ e 4 e 5 e 6 ] d i a g ( [ k t 4 k t 5 k t 6 ] ) [ sign ( e 4 ) α ( e 5 ) α ( e 6 ) ]

References

  1. Chinchilla, M.; Arnaltes, S.; Burgos, J.C. Control of permanent-magnet generators applied to variable-speed wind-energy systems connected to the grid. IEEE Trans. Energy Convers. 2006, 21, 130–135. [Google Scholar] [CrossRef] [Green Version]
  2. Xia, Y.; Ahmed, K.H.; Williams, B.W. Wind turbine power coefficient analysis of a new maximum power point tracking technique. IEEE Trans. Ind. Electron. 2013, 60, 1122–1132. [Google Scholar] [CrossRef]
  3. Yang, X.; Gong, X.; Qiao, W. Mechanical sensorless maximum power tracking control for direct-drive PMSG wind turbines. In Proceedings of the 2010 IEEE Energy Conversion Congress and Exposition (ECCE), Atlanta, GA, USA, 12–16 September 2010. [Google Scholar]
  4. Rocha, R. A sensorless control for a variable speed wind turbine operating at partial load. Renew. Energy 2011, 36, 132–141. [Google Scholar] [CrossRef] [Green Version]
  5. Tan, K.; Islam, S. Optimum control strategies in energy conversion of PMSG wind turbine system without mechanical sensors. IEEE Trans. Energy Convers. 2008, 19, 392–399. [Google Scholar] [CrossRef]
  6. Blaabjerg, F.; Chen, Z.; Kjaer, S.B. Power electronics as efficient interface in dispersed power generation systems. IEEE Trans. Power Electron. 2004, 19, 1184–1194. [Google Scholar] [CrossRef]
  7. Rigatos, G.; Siano, P.; Zervos, N. Sensorless control of distributed power generators with the derivative-free nonlinear Kalman filter. IEEE Trans. Ind. Electron. 2014, 61, 6369–6382. [Google Scholar] [CrossRef]
  8. Pourebrahim, R.; Tohidi, S.; Younesi, A. Sensorless model adaptive control of DFIG by using High frequency signal injection and fuzzy logic control. Iran. J. Electr. Electron. Eng. 2018, 14, 11–21. [Google Scholar]
  9. Qiao, W.; Zhou, W.; Aller, J.M.; Harley, R.G. Wind Speed Estimation Based Sensorless Output Maximization Control for a Wind Turbine Driving a DFIG. IEEE Trans. Power Electron. 2008, 23, 1156–1169. [Google Scholar] [CrossRef]
  10. Qiao, W.; Harley, R.G.; Venayagamoorthy, G.K. Coordinated Reactive Power Control of a Large Wind Farm and a STATCOM Using Heuristic Dynamic Programming. IEEE Trans. Energy Convers. 2009, 24, 493–503. [Google Scholar] [CrossRef] [Green Version]
  11. Fatu, M.; Blaabjerg, F.; Boldea, I. Grid to standalone transition motion-sensorless dual-inverter control of PMSG with asymmetrical grid voltage sags and harmonics filtering. IEEE Trans. Power Electron. 2014, 29, 3463–3472. [Google Scholar] [CrossRef]
  12. Song, X.; Fang, J.; Han, B.; Zheng, S. Adaptive Compensation Method for High-Speed Surface PMSM Sensorless Drives of EMF-Based Position Estimation Error. IEEE Trans. Power Electron. 2016, 31, 1438–1449. [Google Scholar] [CrossRef]
  13. Abdelrahem, M.; Hackl, C.M.; Kennel, R. Finite Position Set-Phase Locked Loop for Sensorless Control of Direct-Driven Permanent-Magnet Synchronous Generators. IEEE Trans. Power Electron. 2018, 33, 3097–3105. [Google Scholar] [CrossRef]
  14. Abdelrahem, M.; Hackl, C.M.; Kenne, R.; Rodríguez, J. Computationally Efficient Finite-Position-Set-Phase-Locked Loop for Sensorless Control of PMSGs in Wind Turbine Applications. IEEE Trans. Power Electron. 2021, 36, 3007–3016. [Google Scholar] [CrossRef]
  15. Ali, M.A.S.; Mehmood, K.K.; Baloch, S.; Kim, C. Wind-Speed Estimation and Sensorless Control for SPMSG-Based WECS Using LMI-Based SMC. IEEE Access 2020, 8, 26524–26535. [Google Scholar] [CrossRef]
  16. Liang, D.; Li, J.; Qu, R.; Kong, W. Adaptive Second-Order Sliding-Mode Observer for PMSM Sensorless Control Considering VSI Nonlinearity. IEEE Trans. Power Electron. 2018, 33, 8994–9004. [Google Scholar] [CrossRef]
  17. Wang, Y.; Xu, Y.; Zou, J. Sliding-Mode Sensorless Control of PMSM with Inverter Nonlinearity Compensation. IEEE Trans. Power Electron. 2019, 34, 10206–10220. [Google Scholar] [CrossRef]
  18. Pan, Y.; Yang, C.; Pan, L.; Yu, H. Integral Sliding Mode Control: Performance, Modification, and Improvement. IEEE Trans. Ind. Inform. 2018, 14, 3087–3096. [Google Scholar] [CrossRef]
  19. Zhang, Z.; Zhao, Y.; Qiao, W.; Qu, L. A space-vector-modulated sensorless direct-torque control for direct-drive PMSG wind turbines. IEEE Trans. Ind. Appl. 2014, 50, 2331–2341. [Google Scholar] [CrossRef]
  20. Yan, J.; Lin, H.; Feng, Y.; Guo, X.; Huang, Y.; Zhu, Z.Q. Improved sliding mode model reference adaptive system speed observer for fuzzy control of direct-drive permanent magnet synchronous generator wind power generation system. IET Renew. Power Gener. 2013, 7, 28–35. [Google Scholar] [CrossRef]
  21. Beltran, B.; Benbouzid, M.E.H.; Ahmed-Ali, T. Second-Order Sliding Mode Control of a Doubly Fed Induction Generator Driven Wind Turbine. IEEE Trans. Energy Convers. 2012, 27, 261–269. [Google Scholar] [CrossRef] [Green Version]
  22. Benbouzid, M.E.H.; Beltran, B.; Amirat, Y.; Yao, G.; Han, J. Second-Order Sliding Mode Control for DFIG-Based Wind Turbines Fault Ride-Through Capability Enhancement. ISA Trans. 2014, 53, 827–833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Besancon, G.; Leon-Morales, J.D.; Huerta-Guevara, O. On adaptive observers for state affine systems. Int. J. Control 2006, 79, 581–591. [Google Scholar] [CrossRef]
  24. Haque, M.E.; Negnevitsky, M.; Muttaqi, K.M. A novel control strategy for a variable-speed wind turbine with a permanent-magnet synchronous generator. IEEE Trans. Ind. Appl. 2010, 46, 331–339. [Google Scholar] [CrossRef] [Green Version]
  25. Geng, H.; Xu, D. Stability analysis and improvements for variable-speed multipole permanent magnet synchronous generator-based wind energy conversion system. IEEE Trans. Sustain. Energy 2011, 2, 459–467. [Google Scholar] [CrossRef]
  26. Geng, H.; Yang, G.; Xu, D.; Wu, B. Unified Power Control for PMSG-Based WECS Operating Under Different Grid Conditions. IEEE Trans. Energy Convers. 2011, 26, 822–830. [Google Scholar] [CrossRef]
  27. Debnath, S.; Saeedifard, M. A new hybrid modular multilevel converter for grid connection of large wind turbines. IEEE Trans. Sustain. Energy 2013, 4, 1051–1064. [Google Scholar] [CrossRef]
  28. Li, S.; Haskew, T.A.; Xu, L. Conventional and novel control designs for direct driven PMSG wind turbines. Electr. Power Syst. Res. 2010, 80, 328–338. [Google Scholar] [CrossRef]
  29. Giraldo, E.; Garces, A. An adaptive control strategy for a wind energy conversion system based on PWM-CSC and PMSG. IEEE Trans. Power Syst. 2014, 29, 1446–1453. [Google Scholar] [CrossRef]
  30. Muljadi, E.; Butterfield, C.P.; Parsons, B.; Ellis, A. Effect of variable speed wind turbine generator on stability of a weak grid. IEEE Trans. Energy Convers. 2007, 22, 29–36. [Google Scholar] [CrossRef]
  31. Wang, W.; Beddard, A.; Barnes, M.; Marjanovic, O. Analysis of active power control for VSC–HVDC. IEEE Trans. Power Deliv. 2014, 29, 1978–1988. [Google Scholar] [CrossRef]
  32. Yuan, X.; Wang, F.; Boroyevich, D.; Li, Y.; Burgos, R. DC-link voltage control of a full power converter for wind generator operating in weak-grid systems. IEEE Trans. Power Electron. 2009, 24, 2178–2192. [Google Scholar] [CrossRef]
  33. Strachan, N.P.W.; Jovcic, D. Stability of a variable-speed permanent magnet wind generator with weak AC grids. IEEE Transactions on Power Delivery 2010, 25, 2779–2788. [Google Scholar] [CrossRef]
  34. Xi, X.; Geng, H.; Yang, G. Enhanced model of the doubly fed induction generator-based wind farm for small-signal stability studies of weak power system. IET Renew. Power Gener. 2014, 8, 765–774. [Google Scholar] [CrossRef]
  35. Mitra, P.; Zhang, L.; Harnefors, L. Offshore wind integration to a weak grid by VSC-HVDC links using power-synchronization control: A case study. IEEE Trans. Power Deliv. 2014, 29, 453–461. [Google Scholar] [CrossRef]
  36. Neris, A.S.; Vovos, N.A.; Giannakopoulos, G.B. A variable speed wind energy conversion scheme for connection to weak AC systems. IEEE Trans. Energy Convers. 1999, 14, 122–127. [Google Scholar] [CrossRef]
  37. Lin, W.M.; Hong, C.M.; Lee, M.R.; Huang, C.H.; Huang, C.C.; Wu, B.L. Fuzzy sliding mode-based control for PMSG maximum wind energy capture with compensated pitch angle. In Proceedings of the 2010 International Symposium on Computer, Communication, Control and Automation (3CA), Tainan, Taiwan, 5–7 May 2010; pp. 397–400. [Google Scholar]
  38. Lin, W.M.; Hong, C.M.; Chen, C.H. Neural-network-based MPPT control of a stand-alone hybrid power generation system. IEEE Trans. Power Electron. 2011, 26, 3571–3581. [Google Scholar] [CrossRef]
  39. Koutiva, X.I.; Vrionis, T.D.; Vovos, N.A.; Giannakopoulos, G.B. Optimal integration of an offshore wind farm to a weak AC grid. IEEE Trans. Power Deliv. 2006, 21, 987–994. [Google Scholar] [CrossRef]
  40. Magri, A.E.; Giri, F.; Besancon, G.; Fadili, A.E.; Dugard, L.; Chaoui, F.Z. Sensorless adaptive output feedback control of wind energy systems with PMS generators. Control Eng. Pract. 2013, 21, 530–543. [Google Scholar] [CrossRef] [Green Version]
  41. Shotorbani, A.M.; Mohammadi, B.; Wang, L.; Marzband, M.; Sabahi, M. Application of finite-time control Lyapunov function in low-power PMSG wind energy conversion systems for sensorless MPPT. Electr. Power Energy Syst. 2019, 106, 169–182. [Google Scholar] [CrossRef]
  42. Khalil, H.K. Nonlinear Systems, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2011. [Google Scholar]
  43. Sami, I.; Ullah, S.; Ullah, N.; Ro, J.S. Sensorless fractional order composite sliding mode control design for wind generation system. ISA Trans. 2020. [Google Scholar] [CrossRef] [PubMed]
  44. Gonzales, L.G.; Figueres, E.; Garcera, G.; Carranza, O. Maximum-power-point tracking with reduced mechanical stress applied to wind-energy-conversion-systems. Appl. Energy 2010, 87, 2304–2312. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Permanent magnet synchronous generator (PMSG)-based wind energy conversion system.
Figure 1. Permanent magnet synchronous generator (PMSG)-based wind energy conversion system.
Energies 14 01712 g001
Figure 2. PI control scheme for comparison with the proposed finite-time controller (FTC).
Figure 2. PI control scheme for comparison with the proposed finite-time controller (FTC).
Energies 14 01712 g002
Figure 3. Block diagram of the proposed FTC.
Figure 3. Block diagram of the proposed FTC.
Energies 14 01712 g003
Figure 4. Input and output of the wind turbine: (a) wind speed profile and (b) mechanical output torque of the wind turbine.
Figure 4. Input and output of the wind turbine: (a) wind speed profile and (b) mechanical output torque of the wind turbine.
Energies 14 01712 g004
Figure 5. Rotor speed tracking by the ECC and the FTC.
Figure 5. Rotor speed tracking by the ECC and the FTC.
Energies 14 01712 g005
Figure 6. The unextracted wind power due to the transients in reference tracking of the FTC and the ECC.
Figure 6. The unextracted wind power due to the transients in reference tracking of the FTC and the ECC.
Energies 14 01712 g006
Figure 7. DC-link voltage.
Figure 7. DC-link voltage.
Energies 14 01712 g007
Figure 8. Reference control outputs of the MSC and the GSC by the FTC.
Figure 8. Reference control outputs of the MSC and the GSC by the FTC.
Energies 14 01712 g008
Figure 9. Comparison of the proposed FTC to the ECC and the PI controller: rotor speed reference-tracking performance.
Figure 9. Comparison of the proposed FTC to the ECC and the PI controller: rotor speed reference-tracking performance.
Energies 14 01712 g009
Figure 10. Average control efforts of the FTC, the ECC, and the PI controller for nominal and perturbed systems.
Figure 10. Average control efforts of the FTC, the ECC, and the PI controller for nominal and perturbed systems.
Energies 14 01712 g010
Figure 11. Bode diagram of the closed-loop WTG system with the proposed PI controller.
Figure 11. Bode diagram of the closed-loop WTG system with the proposed PI controller.
Energies 14 01712 g011
Figure 12. Responses of FTCs and ECCs under nominal and perturbed conditions.
Figure 12. Responses of FTCs and ECCs under nominal and perturbed conditions.
Energies 14 01712 g012
Figure 13. Effects of an increase in k i and k t on responses of the perturbed WTG with the FTC and ECC: (a) settling time and (b) relative steady-state error.
Figure 13. Effects of an increase in k i and k t on responses of the perturbed WTG with the FTC and ECC: (a) settling time and (b) relative steady-state error.
Energies 14 01712 g013
Figure 14. Mechanical output from the adaptive observer: (a) rotor speed ω r and (b) mechanical torque T m .
Figure 14. Mechanical output from the adaptive observer: (a) rotor speed ω r and (b) mechanical torque T m .
Energies 14 01712 g014
Figure 15. Impact of the grid strength on the PCC voltage magnitude with different controllers: (a) proposed FTC and (b) conventional PI controller.
Figure 15. Impact of the grid strength on the PCC voltage magnitude with different controllers: (a) proposed FTC and (b) conventional PI controller.
Energies 14 01712 g015
Table 1. Parameter discription of the wind turbine generator (WTG) model.
Table 1. Parameter discription of the wind turbine generator (WTG) model.
SymbolQuantityValue
J Rotor inertia1.06 × 107 Nm/rad/s2
F Viscous coefficient1.417 Nm/rad/s
P Number of pole pairs145
K M Flux generator constant2.6354 × 103
r s PMSG stator resistance3 mΩ
L s PMSG stator inductance1 mH
E n Grid voltage4.16 kV
L 0 Output inductance1 mH
ω n Grid voltage frequency2π × 50 rad/s
Table 2. Parameters of the WTG used in the case study.
Table 2. Parameters of the WTG used in the case study.
QuantityValue
Turbine-rated torque3.226 × 106 Nm
Rated wind speed11.8 m/s
PMSG-rated speed1.55 rad/s
PMSG-rated power5 MVA
Table 3. Parameters of the controllers and the observers.
Table 3. Parameters of the controllers and the observers.
QuantityValue
Control   gains   k i ;   i = 1 6 2.7, 9.3 × 103 330, 4, 3.1 × 103, 2.1 × 103
Finite - time   gain   k t 1
Sign approximation gain ε 20
PI   controller   proportional   gains   k P _ i ;   i = 1 7 5 × 104, 1.5, 1.5, 5, 1.1, 1.36, 0.01
PI   controller   integral   gains   k I _ i ;   i = 1 7 2.5 × 105, 1.2 × 103, 1 × 103, 6 × 102, 8 × 102, 1.2 × 103, 3
Observer adaptive gains η 0 , η 1 , η 2 , δ 1.5 × 104, 2.3 × 104, 1.5 × 104, 1 × 10−6
Table 4. Mean available and injected active power (MW).
Table 4. Mean available and injected active power (MW).
ControllerFTCECCPI Controller
Maximum available 2.091
Measurement based1.8541.8351.817
Observer based1.8461.8241.806
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pourebrahim, R.; Shotorbani, A.M.; Márquez, F.P.G.; Tohidi, S.; Mohammadi-Ivatloo, B. Robust Control of a PMSG-Based Wind Turbine Generator Using Lyapunov Function. Energies 2021, 14, 1712. https://doi.org/10.3390/en14061712

AMA Style

Pourebrahim R, Shotorbani AM, Márquez FPG, Tohidi S, Mohammadi-Ivatloo B. Robust Control of a PMSG-Based Wind Turbine Generator Using Lyapunov Function. Energies. 2021; 14(6):1712. https://doi.org/10.3390/en14061712

Chicago/Turabian Style

Pourebrahim, Roghayyeh, Amin Mohammadpour Shotorbani, Fausto Pedro García Márquez, Sajjad Tohidi, and Behnam Mohammadi-Ivatloo. 2021. "Robust Control of a PMSG-Based Wind Turbine Generator Using Lyapunov Function" Energies 14, no. 6: 1712. https://doi.org/10.3390/en14061712

APA Style

Pourebrahim, R., Shotorbani, A. M., Márquez, F. P. G., Tohidi, S., & Mohammadi-Ivatloo, B. (2021). Robust Control of a PMSG-Based Wind Turbine Generator Using Lyapunov Function. Energies, 14(6), 1712. https://doi.org/10.3390/en14061712

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