[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Assessment of Stakeholder Benefits from Participating in Community-Shared Solar Photovoltaics Through Monthly Renting and Load Management in South Korea
Previous Article in Journal
Environmental-, Social-, and Governance-Oriented Pension Funds for Young Contributors: A Win–Win Option
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fault-Tolerant Control Implemented for Sustainable Active and Reactive Regulation of a Wind Energy Generation System

by
Adolfo R. Lopez
1,
Jesse Y. Rumbo-Morales
2,
Gerardo Ortiz-Torres
2,
Jesus E. Valdez-Resendiz
3,*,
Gerardo Vazquez
1 and
Julio C. Rosas-Caro
4
1
Laboratory of Electrical and Power Electronics, Tecnologico Nacional de Mexico/ITS de Irapuato, Irapuato 36821, Mexico
2
Centro Universitario de los Valles, Universidad de Guadalajara, Carretera Guadalajara-Ameca Km 45.5, Ameca 46600, Mexico
3
School of Engineering and Sciences, Tecnologico de Monterrey, Monterrey 64700, Mexico
4
Facultad de Ingenieria, Universidad Panamericana, Zapopan 45010, Mexico
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(24), 10875; https://doi.org/10.3390/su162410875
Submission received: 30 October 2024 / Revised: 21 November 2024 / Accepted: 26 November 2024 / Published: 12 December 2024
(This article belongs to the Special Issue Power Electronics on Recent Sustainable Energy Conversion Systems)

Abstract

:
This paper presents the design of a fault-tolerant control system based on fault estimation, aimed at enhancing the sustainability and efficiency of a wind energy conversion system using a doubly-fed induction generator. The control architecture comprises a rotor-side converter (RSC) and a grid-side converter (GSC). The RSC is responsible for regulating both active and reactive power, and its model incorporates two linear subsystem representations. A fault-tolerant control (FTC) scheme is developed using a state-feedback controller; this controller is applied to regulate stator and rotor currents. Additionally, for comparison purposes, Proportional–Integral (PI) and Sliding-Mode Controllers (SMCs) are designed to analyze the performance of each controller. Furthermore, a proportional integral observer is employed in the proposed fault-tolerant scheme for actuator fault estimation. Fault detection is achieved by comparing the fault estimation signal with a predefined threshold. The main contribution of this work is the design and validation of a comprehensive active FTC scheme that enhances system reliability and sustainability. It also includes a performance analysis comparing three controllers (PI, SMC, and state-feedback) applied to the RSC. These controllers are evaluated for their ability to regulate active and reactive power in a wind energy conversion system under conditions of non-constant actuator faults.

1. Introduction

Nowadays, the use of renewable energies has made a significant advance in recent years due to the high cost and negative impact of typical fossil fuels. The total electricity generation in the world increased by 2.3% in 2022, reaching 29,165 TWh; renewable energy sources made up 92% of this electricity increase [1].
Wind Energy Conversion Systems (WECS) involve converting the kinetic energy of the wind into electrical energy with the help of a generator. In recent years, technological advances have helped to develop more economical wind generation systems of great power and higher efficiency. These WECS are divided into three stages: the aerodynamic stage, the mechanical stage, and finally, the electrical stage (see Figure 1).
The generators commonly employed in Wind Energy Conversion Systems (WECS) can be broadly classified into two main categories: synchronous generators and asynchronous generators. The asynchronous generators extensively utilized in these systems include the Squirrel Cage Induction Generators (SCIG), the Wound Rotor Induction Generators (WRIG), and the Doubly-Fed Induction Generators (DFIG) [2].

1.1. Case Study

The case study of this paper focuses on a WECS based on a DFIG; this generator is connected to a multiplier gearbox attached to the rotation propeller shaft. The advantage of the DFIG is that both the rotor and stator windings can be accessed; in this case, the stator is fed directly to the network while the rotor is connected to a back-to-back converter. This converter consists of a Rotor-Side Converter (RSC), a Grid-Side Converter (GSC), and a DC bus. With this configuration, the converter only needs to manipulate the slip power, which oscillates around 30% of the nominal power of the machine [3,4]. Figure 2 shows the described WECS as well as the location of the RSC and the GSC.
Regarding control, the RSC is pivotal in overseeing and managing the active and reactive power transfer between the stator of the DFIG and the grid. This control is effectively executed by manipulating the rotor currents. Simultaneously, the GSC regulates the DC bus. Furthermore, the GSC controls the reactive power exchanged between the converter and the grid. This dual functionality underscores the importance of the GSC in maintaining stability and efficient power transfer within the system [5,6].
The control of the RSC involves meticulously computing the rotor references and generating the stator currents necessary to achieve the targeted levels of active and reactive power. The derivation of these references is facilitated through algebraic relations that intricately incorporate power variables, along with the rotor and stator currents of the DFIG. This methodical approach ensures a precise and coordinated control mechanism, where the algebraic relationships form the basis for determining the requisite parameters to optimize the performance of the DFIG [5,7]. The active and reactive power control of the RSC in WECS has been widely used in the literature [6,8,9,10], and good performance of such a controller is needed to obtain the needed value.
Given these characteristics and the increase in performance requirements at different operation conditions, when connected to the grid, WECS has become progressively vulnerable to faults. These faults are usually voltage harmonics, which in turn can generate voltage ripple, power or torque pulsations, overheating, and increased losses in the stator windings. This inevitably influences process stability, reliability, and safety dynamics [11].

1.2. Literature Review

When discussing the Fault Ride-Through (FRT) mechanism implemented in WECS based on DFIG, in [12], the authors implement a series converter on the stator side during a grid fault, while other papers demonstrate the use of linear control techniques to achieve optimal performance [13,14].
Fault-Tolerant Control (FTC) methods can be considered to identify malfunctions at any time and improve reliability and safety in a process. The FTC techniques are classified into two types [15]: passive and active. On the one hand, in the passive techniques, a fixed controller is designed to tolerate changes in the plant dynamics, and then, the stabilized system could satisfy its goals under all faulty conditions. This approach needs neither Fault Diagnosis (FD) schemes nor controller reconfiguration, but it has limited capabilities because it uses fixed parameters. On the other hand, active techniques adapt or reconfigure the controller parameters based on information from the FD system. This ensures the system’s stability and maintains acceptable performance despite the fault. According to [10], different names can be used to distinguish diagnostic steps: fault detection, fault isolation, and fault identification (also known as fault estimation).
In the context of fault detection and isolation applied to WECS, ref. [16] proposed a diagnostic approach for multiple open-circuit faults in the back-to-back converters of turbines. This diagnostic was implemented by using Permanent Magnet Synchronous Generators (PMSGs). The proposed fault detection method uses the information contained in the instantaneous amplitude (IA) in both the PMSG- and grid-side converters. This method enables accurate distinction between different states of any switch without being affected by speed variations. Simulation and experimental results showed that the proposed strategy detects and locates multiple faults.
In [17], a fault detection methodology of the power converter within a wind turbine chain equipped with a DFIG is presented. In this case, a short-circuit failure in one of the power semiconductor switches is considered and the fault detection method is obtained by the identification of the α - β signals of the RSC. Simulation results have been used to analyze the fault detection and identification method. Additionally, a passive FTC system has been proposed to ensure continuity of service. Fault diagnosis method-based detection for sensors, actuators, and system faults in a wind turbine machine is presented in [18]. Physical and analytical redundancy of sensors and actuators is used to achieve the proposed detection method. The simulation results show the efficiency of the diagnosis scheme.
While a state-of-the-art review on the fault diagnosis of wind power converters is presented in [19], this paper studies the typical fault modes of wind power converters, including short-circuit and open-circuit faults in power switches. Reliability and robustness performance are discussed.
Recent works on FTC design for WECS subjected to faults can be found in [20,21,22,23]. In [20], a nonlinear fault-tolerant approach is presented using a fractional-order non-singular terminal sliding mode controller. This method is designed to maximize power extraction in wind energy conversion systems operating below their rated wind speeds, even in the presence of actuator faults. The sliding surface is proposed to ensure fast finite-time convergence, and the fractional calculus enhances the convergence speed of system states without chattering. Numerical simulation is carried out in the WECS system in both fault-free and faulty scenarios.
A passive FTC system using Fuzzy-PID (Fuzzy-Proportional–Integral Derivative) controller for controlling the DFIG rotor-side converter in fault conditions is presented in [21]. The controller is added only in fault conditions and it reduces the current of the rotor circuit with a proper control of the RSC. The proposed fault control system controls the amplitude of the rotor current when it enters the circuit and prevents it from getting too high. The simulation results guarantee the effectiveness of the proposed method.
In [22] an optimal fault-tolerant pitch control strategy is addressed to adjust the pitch angle of the wind turbine system. The proposed scheme uses a fractional-calculus extended memory of pitch angles with a fractional-order PID controller to enhance the performance of the wind turbine. Fault-tolerant load frequency control design of interconnected wind power systems is investigated in [23]. The FTC is designed with the presence of known and unknown actuator faults to improve the frequency stability of the DFIG-based interconnected wind power systems, including the load disturbances and wind speed fluctuations.
The importance of the FTC in WECS is due to its efficiency. This is because the failure of the components of the WECS leads to low performance. In this case, the FTC keeps an optimal performance by dynamically adjusting the control strategy.
Another topic of importance of the FTC in WECS is security; this is because a failure in WECS could be catastrophic. In this case, the FTC improves security by detecting and mitigating the failures before it becomes a critical problem [24].
One gap in the current research of FTC in WECS is the lack of complex failure scenarios. The research of FTC in WECS focuses on simple or individual failures. More research needs to be carried out under complex circumstances [24,25].
Previous studies on active and reactive power control in wind energy conversion systems, incorporating fault detection, fault estimation, and fault-tolerant control, have seldom addressed non-constant faults in the actuator. In this case, this paper will exclusively focus on the RSC as the actuator. During a fault in the RSC, the DFIG cannot capture the maximum wind power, resulting in inaccurate measurements of active and reactive power. In general, an adequate function of the controller ensures the stability of the overall system [26,27]. Therefore, the motivation of this paper is to consider the RSC actuator partial faults with variable dynamics.
Additionally, this paper presents the design of a feedback controller and an FTC strategy using a proportional–integral observer for detection, estimation, accommodation, and reconfiguration of the fault. The main contribution of this work is the design and validation of a complete active FTC scheme and the performance analysis between three controllers (PI, SMC, and state-feedback) applied to the rotor-side converter for regulating the active and reactive power of the wind energy conversion system in the presence of non-constant actuator faults. The proposed scheme is easy to implement and also offers relevant characteristics compared to the robust controllers presented in this work (PI and SMC); our proposed FTC system provides information about the time of the fault occurrence and its dynamics, and it is possible to accommodate the fault by adding a signal to the nominal control law. The structure of the article is as follows. Section 2 presents the control model of the rotor-side converter, as well as the PI, the sliding-mode controller and the nominal pole placement controller. Section 3 validates the FTC scheme in simulation and the main graphs illustrating the performance of the system. Finally, the discussion and the conclusions are presented in Section 4 and Section 5, respectively.

2. Materials and Methods

2.1. Rotor-Side Converter Control Model

This section shows the control model of the rotor-side converter, which is essential to understanding the design of the different controllers shown in Section 2.2.3.
The model of the DFIG is the same as that of an SCIG type; however, in the case of the DFIG, the rotor is made up of a triphasic wound, which can be accessed through slip rings. The three-phase model of the DFIG is given in [28], which for the sake of brevity is not presented in this paper. This three-phase model can be rewritten in the reference frame d q 0 , which is the product of the transformation K T ( Θ s ) : f a b c f d q 0 of the three-phase model.
The d q 0 reference frame is a mathematical transformation used to represent the three-phase quantities (e.g., voltages or currents) in a rotating reference frame. This transformation simplifies analysis and control by converting three-phase AC signals ( a b c ) into two orthogonal components ( d q ), and a zero-sequence component. This approach is particularly useful in control systems, as it reduces the complexity of the calculations required for the DFIG.
This reference frame is orthogonal and rotating, with angular displacement Θ s . Assuming a balanced three-phase system, the zero-sequence component vanishes, obtaining the following d q model:
v d s = R s i d s ω λ q s + λ ˙ d s ,
v q s = R s i q s + ω λ d s + λ ˙ q s ,
v d r = R r i d r ( ω ω r ) λ q r + λ ˙ d r ,
v q r = R r i q r + ( ω ω r ) λ d r + λ ˙ q r ,
λ d s = L l s + 3 2 L m s i d s + 3 2 L m s i d r ,
λ q s = L l s + 3 2 L m s i q s + 3 2 L m s i q r ,
λ d r = L l r + 3 2 L m s i d r + 3 2 L m s i d s ,
λ q r = L l r + 3 2 L m s i q r + 3 2 L m s i q s ,
where i d , q s and i d , q r are the stator and rotor currents, while V d , q s and V d , q r are the stator and rotor voltages. λ d , q s and λ d , q r represent the rotor flux links; R r and R s are the resistances of the rotor and the stator. L l r is the rotor inductance, L l s is the stator inductance, L m s is the mutual inductance. In addition, ω is the angular velocity of the arbitrary reference frame, and ω r is the angular velocity of the rotor.
As mentioned before, the control of the RSC involves controlling the current that flows through the stator of the DFIG. Because this converter is connected to the rotor terminals, the stator current is controlled indirectly through the rotor currents.
This RSC control involves calculating the current references in the rotor; those references are obtained from algebraic relationships that involve the power variables, as well as the stator and rotor currents [5]. In this case, it was considered to orient the DFIG variables to the stator flow (Figure 3).
Due to this orientation, the equations of λ d s and λ q s are modified, obtaining the following:
λ d s = L s i d s + L M i d r = λ s = L M i m s ,
λ q s = L s i q s + L M i q r = 0 ,
where L s = L l s + 3 2 L m s ,   L M = 3 2 L m s . Clearing i d s and i q s from (9) and (10), and substituting in the equations of λ d r and λ q r , described in (7) and (8), respectively, the following is obtained:
V d r = R r i d r + σ L r d i d r d t c 1 ,
V q r = R r i q r + σ L r d i q r d t + c 2 ,
where
σ = 1 L M 2 L s L r , L r = L r + 3 2 L m s ,
and c 1 and c 2 are compensation terms, defined as:
c 1 = ( ω ω r ) ( L M i m s + σ L r i d r ) ,
c 2 = ( ω ω r ) σ L r i q r .
To design effective controllers for the linear subsystems (11) and (12), these are represented as transfer functions and as a continuous state-space model. Applying the Laplace transform on (11) and (12), without considering the compensation terms, the following are obtained [5]:
I d r ( s ) V d r ( s ) = 1 σ L r s + R r , I q r ( s ) V q r ( s ) = 1 σ L r s + R r .
Also, the nominal state-space model of the previous subsystems (11) and (12) can be represented by:
x ˙ = A x + B u + P , y = C x ,
with
A = R r σ L r 0 0 R r σ L r , B = 1 0 0 1 , P = c 1 σ L r c 2 σ L r , C = 1 σ L r 0 0 1 σ L r ,
where x = [ i d r , i q r ] is the state vector, u = [ V d r , V q r ] is the input vector of the system and y denotes the measurable output vector. In this representation, the compensation terms (14) and (15) are known perturbation vector P.

2.2. Rotor-Side Converter Controller Design

2.2.1. PI Controller

From (16), it can be observed that V d r and V q r depend on i q r , and therefore, the design of the linear controller of the RSC is obtained as [5,7]. The RSC’s linear Proportional–Integral (PI) controllers are designed, obtaining the control loop with the following transfer function:
G c ( s ) = K p s + K i s 2 ( σ L r ) + s ( K p + R r ) + K i .
Figure 4 shows the control scheme, where the conversion of the rotor current from the a b c scheme to the d q reference frame can be observed.
The reference values of P s * and Q s * come from the expressions in the d q reference frame that describe the active and reactive power at the DFIG stator terminal; these are:
P s = 3 2 V d s i d s + V q s i q s , Q s = 3 2 V q s i d s V d s i q s .
Because the stator flux is aligned with the d q reference, V d s = 0 and λ q s = 0 . R s is very small and therefore considered zero; so, V d s and V q s from (1) and (2) are simplified as follows:
V d s = 0 , V q s = ω λ d s = ω L M i m s .
From (9) and (10), the expressions that describe i d s and i q s are obtained:
i d s = L M L S ( i m s i d r ) , i q s = L M L S i q r .
Substituting (20) and (21) in (19), the equations are obtained with which it is possible to calculate reference currents i d r * and i q r * from the reference powers P s * and Q s * ; these are shown below [29]:
i d r * = i m s 2 3 Q s * L s v q r L M , i q r * = 2 3 P s * L s v q r L M .
In this case, a speed of 4200 rpm with an active power reference of 100 W and reactive power of 50 V A R was considered in the WECS. It is important to notice that in this paper, the GSC control, as well as the equations that define the value of the DC link capacitor, will not be aboard.

2.2.2. Sliding-Mode Controller

The basic idea behind SMC is to choose a sliding surface along which the system can be slid to its desired value (see Figure 5). For simplicity, this subsection will only explain the sliding-mode controller applied to the system (11), but the same analysis can be considered for the system (12). Consider the error expressed as e s = P s * ( 1 / σ L r ) i d r . In order to make the error and derivative of error equal to zero, the following sliding-mode surface function is defined:
s = γ e s ,
where γ > 0 . Then, the time derivative of the sliding surface is
s ˙ = γ e ˙ s = γ P ˙ s * γ σ L r i ˙ d r ,
and supposing that the reference value is constant, then by using (11), Equation (24) is rewritten as:
s ˙ = γ σ L r R r σ L r i d r + V d r + c 1 σ L r .
Now, the reaching law is given by [30]:
s ˙ = K s sgn ( s ) ϵ | s | b sgn ( | s | 1 ) s ,
where K s and ϵ are positive parameters, with 0 < b < 1 . The stability condition can be consulted in [30]. By using (25) and (26), and solving for the input V d r , the control law is defined as:
V d r = σ L r γ K s sgn ( s ) ϵ | s | b sgn ( | s | 1 ) s γ R r σ 2 L r 2 i d r + γ c 1 σ 2 L r 2 .

2.2.3. Nominal State-Feedback Controller

In this subsection, a feedback controller is designed for the RSC system, such that the stator and rotor currents tend to the desired reference, as shown in Figure 6. The state-feedback controller is a pole placement technique known as tracking controller as it requires output to follow the input command signal. The integral action reduces the finite error to zero [31]. This controller will be the nominal controller for the active FTC. When a partial fault is detected by the observer, then an alarm is activated and the observer will estimate the magnitude of the fault. Subsequently, a new control signal will be added to the state-feedback control law in order to compensate the fault. The integrator comparator is used to reach the desired time-varying values by defining the following error dynamics:
e ˙ = d c C x ,
where the error dynamics vector is e ˙ = [ e ˙ 1 , e ˙ 2 ] , and the desired reference vector is represented by d c = [ P s * , Q s * ] . By using (17), an augmented system x c = [ x , e ] can be proposed as [32]:
x ˙ c = A ¯ c x c + B ¯ c u + P ¯ + R ¯ d c , y c = C ¯ c x c ,
with
A ¯ c = A 0 2 × 2 C 0 2 × 2 , B ¯ c = B 0 2 × 2 , P ¯ c = P 0 2 × 1 , R ¯ c = 0 2 × 1 1 2 × 1 , C ¯ c = C 0 2 × 2 .
If the pair ( A ¯ c , B ¯ c ) is controllable, then a nominal control law u can be obtained as follows [33]:
u = K ¯ c x c = K 1 K 2 x e B c + P c ,
where the State-Feedback (SF) gain matrices are defined as K 1 and K 2 , and B c + is the pseudoinverse of the matrix B c . The term B c + P c is proposed to eliminate the effect of the vector P c on the controller design. Then, the problem is reduced to determine optimal values of the controller gains.
Then, synthesizing (30) in the linear system (29) and assuming that d c = 0 gives the closed-loop system represented by:
x ˙ c = ( A ¯ c + B ¯ c K ¯ c ) x c .
Global asymptotic convergence of the tracking controller is guaranteed by the following Theorem 1.
Theorem 1.
Given the system (29) and the comparator integrator control law (30), the closed-loop system (31) is asymptotically stable if there exist matrices X = X > 0 and M, such that
X A ¯ c + A ¯ c X B ¯ c M M B ¯ c < 0 .
The gain matrix for the controller is computed by K ¯ c = M X 1 .
The proof can be consulted in [34].
Usually, the Linear Matrix Inequality (LMI) (32) produces slow dynamics. This problem was solved using α -stability, as defined in [35], with the eigenvalues of ( A ¯ c B ¯ c K ¯ c ) located in the desired region, expressed as:
λ ( A ¯ c B ¯ c K ¯ c ) D ,
where D is the α -stability. Then, to avoid slow dynamics, (32) is replaced with (34). Therefore, the following corollary of Theorem 1 is obtained:
Corollary 1.
Given the system (29) and the comparator integrator control law (30), the closed-loop error system (31) is asymptotically stable if there exist a positive scalar α, matrices X = X > 0 and M, such that the condition is met.
X A ¯ c + A ¯ c X B ¯ c M M B ¯ c + 2 α X < 0 ,
and condition (33) holds. The gain matrix of the controller is computed by K ¯ c = M X 1 .

2.3. Active Fault-Tolerant Control Design

The linear state-space representation of the RSC system (17) is used to design and validate the feasibility of a fault-tolerant control (FTC) algorithm in simulation. The nominal pole placement SF controller is considered to validate the proposed FTC algorithm. A Fault Accommodation (FA) scheme is proposed to achieve this goal. In Figure 7, a general FA structure is depicted. Notice here that the scheme is composed of three main blocks: (1) the RSC system, (2) the nominal controller, and (3) the active fault-tolerant control system constructed by the fault detection (FD), estimation (FE), and accommodation subsystems (FA). This scheme is based on our previous work [36].
Only the nominal SF controller is applied when the RSC is in a fault-free case. Once the actuator fault occurs, the FE signal generated by a linear proportional–integral observer (PIO) detects and accommodates the actuator fault. The actuator fault is detected when the FE signal has a more precise value than a predefined threshold. Finally, the FTC law is computed using the FE signal when the fault is detected (faulty case).

2.3.1. Actuator Fault Estimation System

Following some ideas from [32], an observer has been designed for fault estimation. An actuator fault can be modeled as an additive external signal: u f = u + f , where f = [ f 1 , f 2 ] = θ u means the actuator fault vector that affects the RSC system and θ = diag ( θ 1 , θ 2 ) is the actuator loss of effectiveness (LoE). The fault signal is assumed to be constant over time or at least slowly varying, such that f ˙ = 0 . The value of θ i , with i = 1 and 2, indicates (i) a total actuator fault of i-th actuator when θ i = 1 , (ii) the i-th actuator is healthy when θ i = 0 , and (iii) LoE of the actuator when θ i = [ 0 , 1 ] . In this paper, only the case of the actuator LoE presented in the RSC system is considered.
A linear proportional–integral observer will now be applied to the RSC linear system to estimate the actuator fault vector f. The linear system (17) is rewritten in the augmented faulty linear form x a = [ x , f ] , as follows [32]:
x ˙ a = A ˘ a x a + B ˘ a u + P ˘ a , y a = C ˘ a x a ,
with
A ˘ a = A B 0 2 × 2 0 2 × 2 , B ˘ a = B 0 2 × 2 , P ˘ a = P 0 2 × 1 , C ˘ a = C 0 2 × 2 .
A linear PIO is written for the augmented system (35), as follows:
x ^ ˙ a = A ˘ a x ^ a + B ˘ a u + L ˘ P I e y + P ˘ a , y ^ a = C ˘ a x ^ a ,
where e y = y a y ^ a is the estimation error vector, and L ˘ P I = [ L P , L I ] represents the observer gain to be computed for estimate f and x a . The dynamics of the state estimation error e a = x a x ^ a between the augmented system (35) and the augmented PIO (36) are given by:
e ˙ a = ( A ˘ a L ˘ P I C ˘ a ) e a .
Sufficient conditions to guarantee asymptotic stability of the estimation error (37) are given through the following theorem:
Theorem 2.
Consider the linear system (35) and assume that the pair ( A ˘ a , C ˘ a ) is observable. The state and the fault estimation error dynamics (37) are asymptotically stable if there exists S = S and N, such that the following holds:
S A ˘ a + A ˘ a S N C ˘ a C ˘ a N < 0 .
The gain of the observer is computed by L ˘ P I = S 1 N .
The proof can be seen in [32].

2.3.2. Actuator Fault Detection System

The actuator additive fault estimation signal vector is used to detect the actuator fault at any time, as follows:
| f i ^ | μ i in   faulty   case ( Alarm = 1 ) , | f i ^ | < μ i in   fault-free   case ( Alarm = 0 ) ,
where μ i , with i = 1 and 2, is a constant threshold, chosen according to simulation results. In a fault-free case, the estimated value | f i ^ | is close to zero, while in a fault case, the estimated value has a greater value than the threshold to indicate the occurrence of a fault. If the fault estimation value is greater than the threshold, then it is considered a faulty case and the alarm indicator is one.

2.3.3. Fault Accommodation Control Law

When the partial fault is detected and estimated correctly, a new control law is added to the nominal controller to compensate for the effect of the fault on the system. The fault estimation signal is used with the nominal controller to ensure the tracking trajectory performance of the faulty system to the reference. Hence, the FA control law is expressed as u f = u f f ^ , where the first part of the equation is the input with actuator fault and the second part is the additive fault estimation to be added in order to accommodate the fault; see Figure 7.

3. Results

In this section, simulation implementation of the proposed fault-tolerant strategy is presented and validated in Matlab/Simulink 2023B (see Figure 8). Table 1 shows the parameters of the DFIG considered in this paper. Recalling from Section 2, the velocity considered was 4200 rpm, with an active power of 100 W and a reactive power of 50 V A R .
The PI controller gains are K p = 0.41 and K i = 68.54 . The sliding-mode parameters are established as γ = 0.001 , K s = 0.008 , b = 0.08 and ϵ = 100 . The synthesis of the pole placement state-feedback controller and the observer gains have been solved with Yalmip Toolbox [37]. The controller gain is computed by solving Corollary 1 with α = 3 to avoid slow dynamics. Thus, the following controller gain is obtained:
K ¯ c = 1.11 0 237.18 0 0 1.11 0 237.18 .
The observer gain is computed by solving Theorem 2, which relies on:
L ˘ P I = 1.87 0.01 4.78 0.02 0.01 1.87 0.02 4.78 .
For testing the effectiveness of the proposed scheme, two scenarios were considered.

3.1. Scenario 1—Fault-Free

Figure 9 and Figure 10 compare the PI, the sliding-mode and the nominal state-feedback controller applied to the RSC system in a fault-free case. It is easy to see that the three controllers fulfill their objective of regulating the stator and rotor currents. However, the nominal SF controller and the SMC provide fewer tracking errors between the output variables and the references (see Figure 9). The control law signal applied to the RSC system is presented in Figure 10. This figure shows that the sliding-mode control produces a chattering signal. This is a disadvantage compared to the feedback controller, as chattering can lead to wear on the actuators.

3.2. Scenario 2—Non-Constant Actuator Faults

In the second scenario non-constant additive faults are considered, applied as (40) (see Figure 11). In addition, as seen at the top Figure 11, a smaller error is achieved between the real fault and the estimate by using the PI observer. The bottom of Figure 11 shows the actuator fault detection with the alarm indicator.
f 1 = 0 , t < 0.032 s , sin ( 150 t ) , 0.032 s t 0.1 s , 78 t 10.14 , 0.13 s t 0.18 s , 0 , t 0.032 s . f 2 = 0 , t < 0.032 s , 60 t 7.8 , 0.032 s t 0.1 s , 1.6 sin ( 150 t ) + 1.6 , 0.13 s t 0.18 s , 0 , t 0.032 s .
These non-constant faults are applied to the RSC system in order to evaluate the performance of the nominal PI, the nominal SMC and the state-feedback controller adding the proposed FTC scheme. From Figure 12 and Figure 13, a comparison between the faulty RSC system with proportional–integral controller, with the sliding-mode controller, and with the state-feedback controller adding the FTC strategy shows that the proposed scheme is superior to regulating the stator and rotor currents, even with the presence of non-constant actuator faults.
As displayed in Figure 12, the proposed strategy significantly reduces the error between the references and the stator and rotor currents. The actuator fault detection is achieved after approximately 0.0016 s of its occurrence, and consequently, the fault accommodation is generated to compensate for the actuator faults. This reduction in error between the references allows for the achievement of the proposed active and reactive power in the RSC. In addition, Figure 12 shows that the nominal PI controller has bigger error between the desired and the real currents.

4. Discussion

Although the proposed controller and the sliding-mode controller are capable of regulating the faulty RSC system, the proposed FTC strategy can stabilize the system with less error than the SMC, as revealed in Figure 12 and Figure 14, where | e d r | = | P s * i d r | , and | e q r | = | Q s * i q r | . To compare the controller performances, the Integral of Absolute Error (IAE) and the Integral Time Absolute Error (ITAE) are obtained for the tracking error. Table 2 shows that the minimum values of the IAE and the ITAE are obtained for the state-feedback controller with the proposed fault-tolerant control strategy. These results show that the SF controller presents better performance for regulating the stator and rotor currents, even with the presence of non-constant actuator faults.
Figure 13 highlights the performance differences among the PI, SMC, and state-feedback with FTC controllers under non-constant actuator fault conditions. It can be observed that the state-feedback controller with FTC exhibits a rapid and stable response, maintaining control stability despite fault occurrences. In contrast, the SMC controller generates significant chattering, which may negatively impact the actuator’s lifespan and system stability over time. The PI controller, while stable, responds slower than the proposed state-feedback with FTC. Overall, the proposed fault-tolerant controller based on the state-feedback effectively mitigates the effects of faults, showcasing superior performance compared to the other controllers. The proposed controller can promptly command actuators to operate as expected, faster than the PI and the sliding mode controller, effectively reducing non-constant actuator faults.
The results show that the PI and sliding mode controllers are trying to compensate for the faults, but their performance is degraded due to the effect of the actuator fault. Even though the control regulation is robust against faults, these two controllers do not use the fault estimation system, so they cannot bring information about when the fault occurs or the dynamics of the fault. On the other hand, the proposed FTC scheme comprises the fault estimation system and the state-feedback controller. The fault estimation system is designed separately for the controller by considering the separation principle. The proposed scheme is robust against faults and can generate information about the time of the fault and the dynamics of the fault. With these data, fault compensation is achieved.

5. Conclusions

A fault-tolerant control strategy with a fault estimation and compensation for an RSC system of a WECS was proposed in this paper. Three controllers were designed to evaluate the system’s performance under actuator faults: PI, SMC and pole placement state-feedback controller with the proposed FTC system. Fault detection was performed by comparing the fault estimation signal with a predefined threshold. Upon occurrence of actuator faults, the FE signal generates a new control law. The actuator faults are accommodated, and then the performance of the system is improved. The effectiveness and performance of the proposed method and its feasibility in designing fault-tolerant control, fault detection, and identification have been demonstrated in simulation. The results indicate that the SF controller significantly reduces the error between the references and the stator and rotor currents. Additionally, the performance indexes IAE and ITAE applied to the PI, SMC and SF controller are obtained. These results highlight that the SF controller outperforms others in regulating the currents of the RSC system, even in the presence of non-constant actuator faults. Future work will consider the design of sensor fault diagnosis and sensor FTC for the RSC system.
Future work will focus on extending the proposed FTC scheme to address other fault types, such as sensor and electrical faults. It will also explore its application to different system configurations, including various wind turbine models and hybrid energy systems. Key challenges will involve managing computational complexity, ensuring scalability, and enhancing robustness to uncertainties. Additionally, the integration of machine learning and advanced sensor technologies will be explored to improve fault prediction and diagnosis, thereby contributing to more reliable and efficient wind energy systems.

Author Contributions

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

Funding

This research was funded by Consejo Nacional de Humanidades, Ciencias y Tecnologías, under grant CF-2023-G-1344. Desarrollo de nuevos enfoques para análisis, modelado y control de convertidores electrónicos.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. REN21 Members. Renewables 2023 Global Status Report; REN21: Paris, France, 2023; p. 49. [Google Scholar]
  2. Ansari, A.A.; Dyanamina, G.; Ansari, A.A. Decoupled Control of Rotor Side Power Electronic Converter for Grid Connected DFIG Based Wind Energy System. In Proceedings of the 2023 IEEE International Students’ Conference on Electrical, Electronics and Computer Science (SCEECS), Bhopal, India, 18–19 February 2023; pp. 1–6. [Google Scholar] [CrossRef]
  3. Garduño, R.; Rodríguez, A.A.; Sánchez, M.; Martínez, M.A.; Hernández, M.A.; Alcaide, I.X.; Re Lopez y Roberto Hernandez Gonzalez, V.G. Investigación, desarrollo e inovacion tecnológica de sistemas de control de aerogeneradores. Divulgación IIE 2012, 47–55. [Google Scholar]
  4. Hannan, M.A.; Al-Shetwi, A.Q.; Mollik, M.S.; Ker, P.J.; Mannan, M.; Mansor, M.; Al-Masri, H.M.K.; Mahlia, T.M.I. Wind Energy Conversions, Controls, and Applications: A Review for Sustainable Technologies and Directions. Sustainability 2023, 15, 3986. [Google Scholar] [CrossRef]
  5. Calderón, G.; Mina, J.; López, A. Modelado y simulación de un Sistema de Conversión de Energía Eólica de velocidad variable interconectado a la red eléctrica. In Proceedings of the XVI Congreso Latinoamericano de Control Automático, Cancún, Mexico, 14–17 October 2014; pp. 1012–1017. [Google Scholar]
  6. Pande, J.; Nasikkar, P.; Kotecha, K.; Varadarajan, V. A Review of Maximum Power Point Tracking Algorithms for Wind Energy Conversion Systems. J. Mar. Sci. Eng. 2021, 9, 1187. [Google Scholar] [CrossRef]
  7. Pena, R.; Clare, J.; Asher, G. Doubly fed induction generator using back-to-back PWM converters and its application to variable-speed wind-energy generation. IEE Proc.—Electr. Power Appl. 1996, 143, 231. [Google Scholar] [CrossRef]
  8. Ackermann, T. Wind Power in Power Systems; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef]
  9. Aguilar, O.; Tapia, R.; Valderrabano, A.; Minor, H. Design and performance comparison of PI and adaptive current controllers for a WECS. IEEE Lat. Am. Trans. 2015, 13, 1361–1368. [Google Scholar] [CrossRef]
  10. Mohd Zin, A.A.B.; Pesaran, M.H.A.; Khairuddin, A.B.; Jahanshaloo, L.; Shariati, O. An overview on doubly fed induction generators controls and contributions to wind based electricity generation. Renew. Sustain. Energy Rev. 2013, 27, 692–708. [Google Scholar] [CrossRef]
  11. Lopez-Nuñez, A.R.; Mina, J.D.; Aguayo, J.; Calderón, G. Proportional Integral Resonant controller for current harmonics mitigation in a wind energy conversion system. In Proceedings of the 2016 13th International Conference on Power Electronics (CIEP), Guanajuato, Mexico, 20–23 June 2016; pp. 232–237. [Google Scholar] [CrossRef]
  12. Abdel-Baqi, O.; Nasiri, A. A dynamic LVRT solution for doubly fed induction generators. IEEE Trans. Power Electron. 2010, 25, 193–196. [Google Scholar] [CrossRef]
  13. Bharti, O.P.; Saket, R.K.; Nagar, S.K. Design of PI controller for doubly fed induction generator using static output feedback. In Proceedings of the 2015 39th National Systems Conference, NSC 2015, Greater Noida, India, 14–16 December; pp. 2–7. [CrossRef]
  14. Ghafouri, M.; Karaagac, U.; Karimi, H.; Jensen, S.; Mahseredjian, J.; Faried, S.O. An LQR Controller for Damping of Subsynchronous Interaction in DFIG-Based Wind Farms. IEEE Trans. Power Syst. 2017, 32, 4934–4942. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Jiang, J. Bibliographical review on reconfigurable fault-tolerant control systems. Annu. Rev. Control 2008, 32, 229–252. [Google Scholar] [CrossRef]
  16. Xu, S.; Tao, S.; Zheng, W.; Chai, Y.; Ma, M.; Ding, L. Multiple Open-Circuit Fault Diagnosis for Back-to-Back Converter of PMSG Wind Generation System Based on Instantaneous Amplitude Estimation. IEEE Trans. Instrum. Meas. 2021, 70, 1–13. [Google Scholar] [CrossRef]
  17. Zhour, A.; Sihem, G.; Djalel, D. A New Approach to Fault Detection in the Power Converter in Wind Turbine Conversion Systems. Int. J. Robot. Control Syst. 2021, 1, 428–439. [Google Scholar] [CrossRef]
  18. Saci, A.; Cherroun, L.; Hafaifa, A.; Mansour, O. Effective fault diagnosis method for the pitch system, the drive train, and the generator with converter in a wind turbine system. Electr. Eng. 2022, 104, 1967–1983. [Google Scholar] [CrossRef]
  19. Liang, J.; Zhang, K.; Al-Durra, A.; Muyeen, S.; Zhou, D. A state-of-the-art review on wind power converter fault diagnosis. Energy Rep. 2022, 8, 5341–5369. [Google Scholar] [CrossRef]
  20. Mousavi, Y.; Bevan, G.; Küçükdemiral, I.B.; Fekih, A. Maximum power extraction from wind turbines using a fault-tolerant fractional-order nonsingular terminal sliding mode controller. Energies 2021, 14, 5887. [Google Scholar] [CrossRef]
  21. Moradi, A.; Hashemi, S.M.R. Introducing a new control system based on Fuzzy-PID for controlling DFIG rotor side converter in fault conditions. Int. J. Nonlinear Anal. Appl. 2022, 13, 2783–2799. [Google Scholar]
  22. Mousavi, Y.; Bevan, G.; Kucukdemiral, I.B. Fault-tolerant optimal pitch control of wind turbines using dynamic weighted parallel firefly algorithm. ISA Trans. 2022, 128, 301–317. [Google Scholar] [CrossRef]
  23. Kuppusamy, S.; Joo, Y.H.; Kim, H.S. Fault-tolerant load frequency control for DFIG-based interconnected wind power systems. Inf. Sci. 2022, 582, 73–88. [Google Scholar] [CrossRef]
  24. Fekih, A.; Habibi, H.; Simani, S. Fault Diagnosis and Fault Tolerant Control of Wind Turbines: An Overview. Energies 2022, 15, 7186. [Google Scholar] [CrossRef]
  25. Chen, Z.; Yin, M.; Zou, Y.; Meng, K.; Dong, Z. Maximum Wind Energy Extraction for Variable Speed Wind Turbines with Slow Dynamic Behavior. IEEE Trans. Power Syst. 2017, 4, 3321–3322. [Google Scholar] [CrossRef]
  26. Kong, X.; Zhang, Z.; Yin, X.; Wen, M. Study of fault current characteristics of the DFIG considering dynamic response of the RSC. IEEE Trans. Energy Convers. 2014, 29, 278–287. [Google Scholar] [CrossRef]
  27. Patel, R.; Hafiz, F.; Swain, A.; Ukil, A. Nonlinear rotor side converter control of DFIG based wind energy system. Electr. Power Syst. Res. 2021, 198, 107358. [Google Scholar] [CrossRef]
  28. Krause, P.C.; Wasynczuk, O.; Sudhoff, S.D. Analysis of Electric Machinery and Drive Systems, 2nd edi ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2002. [Google Scholar] [CrossRef]
  29. Calderón, G.; Mina, J.; Rosas-Caro, J.C.; Madrigal, M.; Claudio, A.; López, A. Simulation and Comparative Analysis of DFIG-based WECS Using Stator Voltage and Stator Flux Reference Frames. IEEE Lat. Am. Trans. 2017, 15, 1052–1059. [Google Scholar] [CrossRef]
  30. Wang, Y.; Feng, Y.; Zhang, X.; Liang, J. A new reaching law for antidisturbance sliding-mode control of PMSM speed regulation system. IEEE Trans. Power Electron. 2019, 35, 4117–4126. [Google Scholar] [CrossRef]
  31. Nise, N.S. Control Systems Engineering; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2020. [Google Scholar]
  32. Rumbo-Morales, J.Y.; Ortiz-Torres, G.; Sarmiento-Bustos, E.; Rosales, A.M.; Calixto-Rodriguez, M.; Sorcia-Vázquez, F.D.; Pérez-Vidal, A.F.; Rodríguez-Cerda, J.C. Purification and production of bio-ethanol through the control of a pressure swing adsorption plant. Energy 2024, 288, 129853. [Google Scholar] [CrossRef]
  33. Sanchez, R.O.; Rumbo Morales, J.Y.; Ortiz Torres, G.; Pérez Vidal, A.F.; Valdez Resendiz, J.E.; Sorcia Vázquez, F.d.J.; Nava, N.V. Discrete State-Feedback Control Design with D-Stability and Genetic Algorithm for LED Driver Using a Buck Converter. Int. Trans. Electr. Energy Syst. 2022, 2022, 8165149. [Google Scholar] [CrossRef]
  34. Lendek, Z.; Guerra, T.M.; Babuska, R.; De Schutter, B. Stability Analysis and Nonlinear Observer Design Using Takagi-Sugeno Fuzzy Models; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  35. Duan, G.R.; Yu, H.H. LMIs in Control Systems: Analysis, Design and Applications; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  36. Ortiz Torres, G.; Rumbo Morales, J.Y.; Ramos Martinez, M.; Valdez-Martínez, J.S.; Calixto-Rodriguez, M.; Sarmiento-Bustos, E.; Torres Cantero, C.A.; Buenabad-Arias, H.M. Active Fault-Tolerant Control Applied to a pressure swing adsorption process for the production of Bio-hydrogen. Mathematics 2023, 11, 1129. [Google Scholar] [CrossRef]
  37. Lofberg, J. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the 2004 IEEE International Conference on Robotics and Automation (IEEE Cat. No. 04CH37508), Taipei, Taiwan, 2–4 September 2004; pp. 284–289. [Google Scholar]
Figure 1. Stages of a wind energy conversion system.
Figure 1. Stages of a wind energy conversion system.
Sustainability 16 10875 g001
Figure 2. WECS based on DFIG and back-to-back converter.
Figure 2. WECS based on DFIG and back-to-back converter.
Sustainability 16 10875 g002
Figure 3. Vector diagram of the DFIG variables oriented to the stator flux.
Figure 3. Vector diagram of the DFIG variables oriented to the stator flux.
Sustainability 16 10875 g003
Figure 4. RSC PI control scheme.
Figure 4. RSC PI control scheme.
Sustainability 16 10875 g004
Figure 5. RSC sliding-mode control scheme.
Figure 5. RSC sliding-mode control scheme.
Sustainability 16 10875 g005
Figure 6. RSC state-feedback control scheme.
Figure 6. RSC state-feedback control scheme.
Sustainability 16 10875 g006
Figure 7. FTC scheme applied to the RSC system, with (1) the RSC system, (2) the nominal controller, and (3) the active fault-tolerant control system.
Figure 7. FTC scheme applied to the RSC system, with (1) the RSC system, (2) the nominal controller, and (3) the active fault-tolerant control system.
Sustainability 16 10875 g007
Figure 8. Simulation model of the FTC scheme with the state-feedback controller, developed in Matlab/Simulink 2023B, with (1) the RSC system, (2) the nominal controller, and (3) the active fault-tolerant control system.
Figure 8. Simulation model of the FTC scheme with the state-feedback controller, developed in Matlab/Simulink 2023B, with (1) the RSC system, (2) the nominal controller, and (3) the active fault-tolerant control system.
Sustainability 16 10875 g008
Figure 9. Scenario 1—comparison between PI, SMC and nominal state-feedback controller.
Figure 9. Scenario 1—comparison between PI, SMC and nominal state-feedback controller.
Sustainability 16 10875 g009
Figure 10. Scenario 1—input signal comparison between PI, SMC and nominal state-feedback controller.
Figure 10. Scenario 1—input signal comparison between PI, SMC and nominal state-feedback controller.
Sustainability 16 10875 g010
Figure 11. Scenario 2—additive fault estimation.
Figure 11. Scenario 2—additive fault estimation.
Sustainability 16 10875 g011
Figure 12. Scenario 2—comparison between PI, SMC and state-feedback controller with FTC.
Figure 12. Scenario 2—comparison between PI, SMC and state-feedback controller with FTC.
Sustainability 16 10875 g012
Figure 13. Scenario 2—input signal comparison between PI, SMC and state-feedback with FTC.
Figure 13. Scenario 2—input signal comparison between PI, SMC and state-feedback with FTC.
Sustainability 16 10875 g013
Figure 14. Scenario 2—comparison of absolute tracking error between the SMC and the state-feedback with FTC.
Figure 14. Scenario 2—comparison of absolute tracking error between the SMC and the state-feedback with FTC.
Sustainability 16 10875 g014
Table 1. DFIG parameters.
Table 1. DFIG parameters.
ParameterSymbolValueUnit
Rated DFIG PowerP372W
Number of polesp2-
Stator voltage (60 Hz)- 34.29 V
Rotor resistance R r 0.312 Ω
Stator resistance R s 0.343 Ω
Rotor inductance L l r 1.198 mH
Stator inductance L l s 1.198 mH
Mutual inductance L m s 38.62 mH
Table 2. Performance indexes of the PI, SMC and state-feedback controller with FTC applied to the faulty RSC system.
Table 2. Performance indexes of the PI, SMC and state-feedback controller with FTC applied to the faulty RSC system.
ErrorIAEITAE
PISMCSF with FTCPISMCSF with FTC
e d r = P s * i d r 0.29640.0047430.0008900.022810.0006428.411 × 10 6
e q r = Q s * i q r 0.15140.0043940.0005730.020450.0005577.147 × 10 6
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

Lopez, A.R.; Rumbo-Morales, J.Y.; Ortiz-Torres, G.; Valdez-Resendiz, J.E.; Vazquez, G.; Rosas-Caro, J.C. Fault-Tolerant Control Implemented for Sustainable Active and Reactive Regulation of a Wind Energy Generation System. Sustainability 2024, 16, 10875. https://doi.org/10.3390/su162410875

AMA Style

Lopez AR, Rumbo-Morales JY, Ortiz-Torres G, Valdez-Resendiz JE, Vazquez G, Rosas-Caro JC. Fault-Tolerant Control Implemented for Sustainable Active and Reactive Regulation of a Wind Energy Generation System. Sustainability. 2024; 16(24):10875. https://doi.org/10.3390/su162410875

Chicago/Turabian Style

Lopez, Adolfo R., Jesse Y. Rumbo-Morales, Gerardo Ortiz-Torres, Jesus E. Valdez-Resendiz, Gerardo Vazquez, and Julio C. Rosas-Caro. 2024. "Fault-Tolerant Control Implemented for Sustainable Active and Reactive Regulation of a Wind Energy Generation System" Sustainability 16, no. 24: 10875. https://doi.org/10.3390/su162410875

APA Style

Lopez, A. R., Rumbo-Morales, J. Y., Ortiz-Torres, G., Valdez-Resendiz, J. E., Vazquez, G., & Rosas-Caro, J. C. (2024). Fault-Tolerant Control Implemented for Sustainable Active and Reactive Regulation of a Wind Energy Generation System. Sustainability, 16(24), 10875. https://doi.org/10.3390/su162410875

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