Next Article in Journal
Fast Constant-Time Modular Inversion over Fp Resistant to Simple Power Analysis Attacks for IoT Applications
Next Article in Special Issue
Switched-Delay Smith Predictor for the Control of Plants with Response-Delay Asymmetry
Previous Article in Journal
Two-Stage Intelligent Model for Detecting Malicious DDoS Behavior
Previous Article in Special Issue
Multi-Time Scale Trading Simulation of Source Grid Load Storage Based on Continuous Trading Mechanism for China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensor Fault-Tolerant Control of Microgrid Using Robust Sliding-Mode Observer

1
Department of Electrical Engineering, FET, International Islamic University, Islamabad 44000, Pakistan
2
Research and Innovation Centre of Excellence (KIOS CoE), University of Cyprus, P.O. Box 20537, Nicosia 1678, Cyprus
3
Centre of Renewable Energy, Government Advance Technical Training Centre, Peshawar 25100, Pakistan
4
Department of Electrical Engineering, University of Engineering and Technology, Mardan 23200, Pakistan
5
Computer Sciences Program, Turabah University College, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
6
Department of Electrical Engineering, University of Management and Technology Lahore, Sialkot Campus, Sialkot 51310, Pakistan
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(7), 2524; https://doi.org/10.3390/s22072524
Submission received: 12 December 2021 / Revised: 19 March 2022 / Accepted: 21 March 2022 / Published: 25 March 2022
(This article belongs to the Special Issue Nonlinear Control with Applications to Energy Systems)

Abstract

:
This work investigates sensor fault diagnostics and fault-tolerant control for a voltage source converter based microgrid (model) using a sliding-mode observer. It aims to provide a diagnosis of multiple faults (i.e., magnitude, phase, and harmonics) occurring simultaneously or individually in current/potential transformers. A modified algorithm based on convex optimization is used to determine the gains of the sliding-mode observer, which utilizes the feasibility optimization or trace minimization of a Ricatti equation-based modification of H-Infinity ( H ) constrained linear matrix inequalities. The fault and disturbance estimation method is modified and improved with some corrections in previous works. The stability and finite-time reachability of the observers are also presented for the considered faulty and perturbed microgrid system. A proportional-integral (PI) based control is utilized for the conventional regulations required for frequency and voltage sags occurring in a microgrid. However, the same control block features fault-tolerant control (FTC) functionality. It is attained by incorporating a sliding-mode observer to reconstruct the faults of sensors (transformers), which are fed to the control block after correction. Simulation-based analysis is performed by presenting the results of state/output estimation, state/output estimation errors, fault reconstruction, estimated disturbances, and fault-tolerant control performance. Simulations are performed for sinusoidal, constant, linearly increasing, intermittent, sawtooth, and random sort of often occurring sensor faults. However, this paper includes results for the sinusoidal nature voltage/current sensor (transformer) fault and a linearly increasing type of fault, whereas the remaining results are part of the supplementary data file. The comparison analysis is performed in terms of observer gains being estimated by previously used techniques as compared to the proposed modified approach. It also includes the comparison of the voltage-frequency control implemented with and without the incorporation of the used observer based fault estimation and corrections, in the control block. The faults here are considered for voltage/current sensor transformers, but the approach works for a wide range of sensors.

1. Introduction

Distributed electrical power generation by various renewable and non-conventional micro-sources provides diverse options for configuring modern electrical grids. A microgrid (MG) consists of a cluster of loads and distributed generators (DG) that operate as a single controllable system [1]. The interconnection of DGs to the utility grid through power electronic converters has raised concerns about safe operation and protection of equipment, as power management is required at micro levels as well as for large unified networks. MGs may operate in both autonomous and grid-connected modes, whereas stability and effective control are issues of concern for both modes of operation of dynamic microgrids [2]. A properly simulated and experimentally tested small-signal linearized MG model is required for reliable testing of the observer-based FTC techniques used in this work. The FTC approach enhances the safety and reliability of the entire system by generating alarms for fault magnitudes within certain thresholds. Moreover, it will maintain the operation close to normal instead of immediate shutdowns and prevent possible losses, along with providing maintenance opportunities by shifting the power switches to redundant generator sources using the hierarchical MG control strategy [3,4]. This work is specifically focused on the FTC of microgrids, particularly with reference to sensor faults, as the control of power systems relies on sensor measurements. If the sensors fail, have bad or broken connectivity, cause communication malfunction, etc., they will not only fail the control mechanism but also the power system. In general terms, the reliability enhancement, local voltage support, correction of voltage/frequency sags, and uninterruptible power supply (UPS) are major issues of focus in this work [5]. The small-signal modeling of MGs can be studied in detail in [4,6,7,8,9,10,11,12,13].
Sliding-mode observer (SMO), is different from ordinary Luenberger observers because of the nonlinear switch discontinuous term, which is properly processed for fault diagnostics (detection, isolation, and estimation). They are used for both linear/nonlinear systems with additive faults and are capable of handling a large class of perturbations, parametric variations, uncertainties, and unmodeled dynamics. Moreover, it provides compensation for observer mismatches while replicating the system, along with ensuring stability and reachability in finite time [14,15].
Walcot and Zak [16] used feedback output error linearly in observers, along with a Lyapunov-based stability analysis. The faults can be determined from the deviation of the system trajectory from the sliding surface, but it cannot determine which sensor or actuator is faulty. The basic requirement of fault detection is to detect faults, along with location identification, as mentioned in [17,18]. The difference between the system and observer outputs generates residual signals, which are processed further with static and dynamic thresholds for better diagnostics, reconfiguration of the system/observer, and fault reconstruction. Residual based methods were reviewed by [19,20]. SMO-based methods have also been proposed by [21,22]. Utkin [23] in his proposed observer used a discontinuous term with a suitably scaled gain, which requires software to solve the synthesis problem.
Edwards [24] extended the work by using both discontinuous and linear terms as feedback in SMO, which is transformed into canonical forms by suitable transformations constrained by conditions on invariant zeros of the system. The methods for computing the required gains of the observer are given, which provides an explicit solution, but all degrees of freedom are not exploited. The work was carried out by [25] using an output term injected with an appropriate gain for fault magnitude determination in the localized sensor/actuator while maintaining the sliding motion. Tan [26] used sliding-mode observer (SMO) with linear matrix inequality (LMI) based convex optimization algorithms to exploit more degrees of freedom for observer gains and establish a relationship between the sub-optimal observer and the linear part of the SMO using the linear quadratic Gaussian (LQG) theory. The paper also presented a modified algorithm for pole placement and clustering to place eigenvalues of the linear SMO part in a certain region to attain improved observer performance. The algorithms also involve some design matrices inherent in procedures that can modify the dynamics of the sliding motion. Tan [27] applied the SMO theory to find additive sensor faults by transforming them suitably to act as actuator faults. This work was further extended by [28] for localized additive sensor/actuator fault reconstruction with additive uncertainty in the system. It uses the minimization of the L 2 gain of uncertainty and reconstructed fault, scaling of the output injection signal, and application of LMIs for more design freedom. Although the sliding mode is often not retained in most cases in a very ideal sense in the presence of faults, the estimation is shown to be made possible in these works.
Ref. [29] also used the idea of correcting the sensor/actuator faults by using the estimated fault, but the technique could not be termed as strictly active or passive. An observer-based approach for online fault estimation was used to correct the sensor faults but was not active in the sense of updating the control law, which was investigated by Edwards in [30]. Yan [31] imposed some conditions on the fault and disturbance distribution matrices to achieve a specific fault reconstruction. Aldeen [32] estimated the states, faults, and disturbances/unknown inputs for nonlinear systems by designing an SMO.
Dhahri [33] initially presented a fault reconstruction methodology for matched faults and uncertainties, which is not very practical because of the very nature of uncertainties; hence, the method is not valid for a good range of practical systems. They further extended the work by determining the SMO gains by feasibility optimization of LMIs attained through Lyapunov stability criteria. It minimizes the H criteria (ratio of residual to disturbance), and the reduced effect of disturbance on fault reconstruction improves the fault estimation, along with validation of the work for unmatched uncertainties/ faults. This study proved the stability and reachability of the sliding-mode observer and evaluated the reachability time by using a suitable mathematical estimation.
Shi [34] presents H based FTC for actuator and sensor faults occurring in a wind energy system. It considers variable wind speed being modeled with stochastic affine models and uses linear quadratic regulator (LQR)-based state feedback control. This paper discusses various intelligent methods for designing and controlling small-and large-scale MGs operating in parallel or in isolation. The increasing complexity of MG systems with various combinations of active and passive elements requires individual converter units to act in a fault-tolerant manner for better efficiency and reliability in generation, distribution, and transmission [35].
In [36], SMO was applied to sensor faults, that is, for current transformer (C.T) magnitude faults and used an estimated fault for correction of sensor (C.T) fault. It also uses feedback pole placement control to track the real and complex powers, and in this way, a compensating sensor (C.T) fault without even replacement or alarm, which can be better termed as a method to design a software-based sensor. An adaptive-fuzzy-PID robust sliding control is proposed for an uncertain class of nonlinear systems, particularly focusing on aircraft flight control, is referred to in [37].
The work in [38] designed a sliding control with adaptive gains and an integral surface for robustness against actuator faults of a wind turbine system to maintain the rotor speed, with simulation results. Recent approaches are moving towards consensus/distributed controls of MGs along with trending applications of fuzzy control. A distributed adaptive fuzzy-based control is presented for large-scale systems with sensor and actuator faults [39].
Kangdi Lu [40] proposed robust PI control for load-frequency control (LFC) of multi area interconnected scenerio, where its parameters are estimated by constrained population extremal optimization. It uses H determined from LMIs as one constraint whereas the integral based time absolute error as another constraint. The papers showed effectiveness of scheme by its comparison with other PI control methods and an optimized model predictive control (MPC). Kang-Di Lu [41] developed adaptive, resilient event triggered PI-based LFC which works with energy constraint. Lyapunov theory to derive criterion of stability which is tested for various case studies, and is proved to be defeating denial of service attacks and reducing the communication burdens.
Damiano Rotando [42] developed LMI-based observer/controller to be employed for dynamic cost estimation and control for a saturated-actuator class of non-linear systems to maintain the cost below an upper bound. The scheme is tested by simulation results on rotational single-arm-inverted pendulum.
This study focuses on designing a sensor FTC mechanism applicable at both the primary and secondary levels with reference to hierarchical MG control [3] strategy. The FTC strategy uses the robust sliding-mode observer (RSMO) theory to be utilized for the detection and estimation of sensor faults of a distributed MG unit operating in autonomous mode of operation, using a reliably simulated and experimentally verified mathematical model of the MG. The considered system has one current transformer (C.T) and potential transformer (P.T), but it is easily scalable for multiple and redundant sensor sources (transformers), which can be utilized for enhanced security and reliability. The online fault estimation-based fault tolerance approach is quite general and can be extended to many types of systems by choosing the suitable observer/controller parameters discussed in this paper. The main contributions of this study are:.
  • We have presented an improved method, as it is capable of dealing with either one or both faulty sensors (transformers) with sinusoidal additive faults (composite of magnitude, phase, and harmonics). Earlier work was performed only for magnitude faults in the current transformer in [36].
  • Algorithm for determination of SMO gains is modified and improved by using Ricatti-equation-based modification of H enhanced LMIs, i.e., the combination of approaches represented by Edwards [26,33], which enhances robustness of SMO in fault diagnosis and reconstruction. The method for the mentioned approach is discussed in Theorem 2. The comparisons for results produced with SMO gains determined by proposed technique to those determined by earlier/base works (being used in this study), are also presented.
  • Fault estimation procedure is modified and contains some factors which were mistakenly missing in reduced order state error equations of earlier works such as [30]. The method also incorporates numerical solution of the DE for reduced order state error estimation. The mentioned approach is worked in Corollary 1.
  • Disturbances being unknown inputs are also accurately estimated for all fault types considered in the system. The modified fault estimation method gives very accurate disturbance/unknown inputs estimation as shown in Corollary 1.
  • Lyapunov stability analysis of SMO with considered additive faults and disturbances for the considered microgrid system is presented, as shown in Proposition 1.
  • Reachability of sliding mode in finite time is proved along with determination of reachability time, for the considered faulty and perturbed microgrid system, as shown in Theorem 3.
  • Along with voltage-frequency regulation using PI control and SVPWM in earlier studies, FTC is an additional feature in the same control mechanism by using the corrections in faulty sensor outputs achieved through the estimated faults and disturbances by using SMO.
  • The system is simulated for various fault types of practical importance particularly with reference to C.T/P.T and generally for any type of sensors, i.e., sinusoidal, linearly increasing, constant, square pulse, and random type faults with additive sinusoidal harmonics in the form of disturbances.
This paper is organized with an introduction, literature review, and brief discussion of the proposed scheme in Section 1, current/voltage transformer faults model and MG system model in Section 2, preliminaries of SMO based detection and estimation theory in Section 3, determination of H optimized SMO gains in Section 4, reduced-order sliding motion dynamics for fault estimation along with stability and reachability proofs in Section 5, voltage-frequency regulation enhanced with FTC in Section 6, results and discussions on various fault cases in Section 7, and concluding remarks and some possible future directions in Section 8. An Appendix A section is also added after the mention of funding and references. A simplified block diagram illustration of this study is presented in Figure 1.

2. System Model

2.1. Modelling of Current and Voltage Transformer Faults

The actual problem in C.T/P.T is saturation, which causes an increased core magnetization current and reduces secondary current compared to that required by relays and switches. The CT saturation on the primary side is caused by DC offset of the fault current preceding to the additional effects caused by the remnant flux of fault. The relay coils and even wires causing the cumulative impedance burden of the secondary side also plays its role in CT saturation. Hence, protection schemes require several C.Ts connected in parallel [43,44,45]. The inaccuracies of P.T/C.T, their output will result in incorrect signals generated by space-vector-pulse-width-modulation (SVPWM) for voltage source converter (VSC). The VSC will generate consequently the inaccurate voltages/currents, which will not track the required reactive and active powers. The errors of P.T/C.T are sinusoidal in nature inclusive of magnitude, harmonics and phase, whereas the previous studies have considered only magnitude faults for the sake of simplicity [36]. The additive faults/disturbances are considered in the system, such that the sensor output is added up with different sinusoidal faults/disturbances, which compositely form the complete mathematical model of the faults actually occurring in the C.Ts/P.Ts.
f ( t ) + ξ ( t ) = f o * s i n ( ω 1 t + ϕ 1 ) + ξ o * s i n ( ω 2 t + ϕ 2 )

2.2. Mathematical Model of Microgrid System

A small-signal MG model is used in this study to determine the workability of the proposed FTC scheme. The MG model is properly simulated, experimentally verified, and able to be used as a block in large grid networks [6].
The mathematical model of a VSC-based microgrid includes non-linear equations; however, its linearized version at the operating point is used in this study. The microgrid model, as proposed by [6] includes inverter equations, voltage source representing grid, LCL filter, PLL, current controller, droop control equations, voltage controller, SVPWM, and reactive and real power calculations, all connected in a cascade to give a non-linear model that needs linearization at the operating point. The LCL filters are connected as passive filters with the VSI, to manage the voltage/current spikes. The line/stray inductances, capacitances, and resistances are also modeled as being considered with LCL filters, i.e., r c and r f are the inductive parasitic resistances, whereas the damping resistor R d is serially connected to the filter capacitor. However, ESR of the capacitor is not considered explicitly instead its lumped into R d . The C.T/P.T are mounted on LCL filter to read the instantaneous voltage/current readings, which are used to determine the instantaneous values of reactive and real ( Q c , P r ) powers. The detailed inverter model as provided in Simulink is considered in this work, i.e., without any major erroneous behavior, it is assumed that the voltage being commanded appears at the input of the filter inductor i.e., V i d q * = V i d q . The losses in the diodes and IGBTs are neglected. The LCL part of the system modeled with KVL/KCL along with average VSI model and grid source taken together is considered as the required part of the mathematical model to design the SMO. Therefore, if there are no sensor faults in (P.T/C.T) of the system, the actual system states/outputs and those of model are in compromise, showing no difference. The considered MG system model is a continuous time LTI system, which follows the separable principle, i.e., the observer and controller can work in combination, and the control action can be performed on the the unknown system states being determined by an observer. The Simulink based model is shown in Figure 2.
Definition 1
(Principle of Separation of Estimation and Control). Considering deterministic linear systems, if state observer based estimated states are used to design a state feedback control for a LTI system, then the observer and feedback controller in combination are stable.
The linearized state space MG model is given as:
i ˙ i i ˙ o v ˙ o = r f / L f 0 1 / L f 0 r c / L c 1 / L c 1 / C f r f R d / L f ( 1 / C f r c R d / L c ) ( R d / L f + R d / L c ) i i i o v o
+ 1 / L f 0 1 / L f v i + 0 1 / L c R d / L c v g
The simulink based microgrid schemetic diagram (SLD) is given in Figure 2. The abc-dq0 transformation is a combination of the Park and Clark transformation for three-phase current/voltage, which is defined for any signal s(t) as
s d s q s 0 2 / 3 c o s ( θ ) c o s ( θ 2 π 3 ) c o s ( θ + 2 π 3 ) s i n ( θ ) s i n ( θ 2 π 3 ) s i n ( θ + 2 π 3 ) 2 / 2 2 / 2 2 / 2 s a s b s c
The frequency of the system is determined using PLL while working with dq frame voltages. The grid side angle of voltage is measured by PLL in the considered system, which is then used for all abc-dq/ dq-abc transformations, needed by the system. The dq transformed system model in Equation (2) is
A s = r f L f W p L L 0 0 1 L f 0 w p L L r f L f 0 0 0 1 L f 0 0 r c L c w p L L 1 L c 0 0 0 w p L L r c L c 0 1 L c 1 C f r f R d L f w p L L R d 1 C f + r c R d L c w p L L R d ( w p L L + R d L f + R d L c ) 0 w p L L R d 1 C f r f R d L f w p L L R d 1 C f + r c R d L c 0 ( w p L L + R d L f + R d L c )
B s = 1 / L f 0 0 1 / L f 0 0 0 0 R d / L f 0 0 R d / L f , B g = 0 0 0 0 1 / L c 0 0 1 / L c R d / L c 0 0 R d / L c , C s = 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 , E s = D s 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
x s = [ I i d , I i q , I o d , I o q , V o d , V o q ] T , w = [ V g d , V g q ] T , u = [ V i d , V i q ] T , ξ = [ ξ i d , ξ i q , ξ v d , ξ i d ] T , f = [ f i d , f i q , f v d , f v q ] T ,
The dq-axis output current and voltage measurements are used to determine the instantaneous reactive power (Q c ) and active power ( P r ) generated by the inverter.
P r = 3 / 2 * ( V o q I o q + V o d I o d )
Q c = 3 / 2 * ( V o q I o d V o d I o q )
where I o d , I o q , V o d , V o q are the dq and components of sensor (P.T/C.T) output currents/voltages. Instantaneous powers are low pass filtered using the corner frequency ω c to obtain the output power filtered with high frequency effects. The single line diagram (SLD) of the considered VSI-based MG is given Figure 2.
The generalized state-space model of system is:
x s ˙ = A s x s + B s u + B g w
y s = C s x s + E s f + D s ξ
In general form, the dimensions of vectors are x s R n * 1 , x h R p * 1 , w R m * 1 , u R m * 1 , f R q * 1 , ξ R q * 1 , y s R p * 1 , whereas the dimensions of system matrix, matrix of grid side dynamics, actuator matrix, and output matrix respectively, are A s R n * n , B g R n * m , B s R n * m , and C s R p * n , i.e., [ 0 p * ( n p ) , I p * p ] ; whereas dimensions of disturbance and fault distribution matrices with full column and row rank are D s R q * q = I q * q , E s R q * q = I q * q respectively, where n p q .
The Euclidean norm of any vector y is defined by
y y T y
The boundedness of the disturbance and fault magnitudes in terms of the Eucledian norm is given by f α and ξ < ξ o , which are design requirements, along with the above-mentioned matrix dimensions.

2.3. Stable Filtering and Augmented System

A stable filter is used to make the sensor outputs least depending on faults and to magnify the small/ insignificant faults [30]. The system and stable filtered output state vectors are augmented, where the state vector also includes the unscaled actual output variables. In this way, the stable filtered scaled augmented states also provides the isolation of faulty sensors, that is required for proper diagnosis process. A stable filter in some cases may be is a positive definite (PD) scaled identity matrix or a stable matrix with eigenvalues in the left half plane or a high-frequency noise suppressing first-order filter for the output signal in some applications. The state space form of stable filter equation is
x h ˙ = A h y s A h x h
where the stable filter matrix A h R p * p
x h ˙ = A h x h + A h C s x s + A h E s f + A h D s ξ
The states and stable filtered outputs are augmented for easy handling in a compact form. The system in augmented form is
x ˙ s x ˙ h = A s 0 A h C s A h x s x h + B s B g 0 0 u w + 0 A h E s f + 0 A h D s ξ
x ˙ c = A c x c + B c u c + E c f + D c ξ
y c = C c x c
where y c R p * 1
y c = 0 I x s x h
where x c = x s x h u c = u w , A c = A s 0 A h C s A h , B c = B s B g 0 0 , E c = 0 A h E s and D c = 0 A h D s
Dimensions of vectors are x c R n c * 1 , u c R 2 m * 1 , x h R p * 1 , whereas the dimensions of augmented system matrices are A c R n c * n c , B c R n c * p c , C c R p c * n c , i.e., [ 0 p c * ( n c p c ) , I p c * p c ] , A h D s = D o R q * q , A h E s = E o R q * q , where D c R n c * q c and E c R n c * q c , are disturbance and fault distribution matrices, respectively, in augmented system such that D c = E c = 0 n * p I p * p and for the MG system under consideration p c = p = 4 , q c = q = 4 , n c = n + p = 10 .
Remark 1.
The invariant zeros of ( A c , E c , C c ) λ ( A s ) if r a n k ( C c E c ) q as shown in [30], so the system ( A c , E c , C c ) will be minimum phase if the open loop system is stable. In other words, if the system has less inputs than outputs, that is, q p , then it is expected that the system will not have any invariant zeros.

3. Fault Diagnosis SMO Filtering

In the literature the phrase ‘fault diagnosis’ refers to three main objectives: detection, isolation, and estimation of faults occurring in the system [46]. The fault detection is often not sufficient in most of the cases, instead it also requires isolation of the fault location along with estimation of faults so that the corrective mechanism can be managed to ensure protection of the systems. These objectives are achieved using the methods which can be categorized into four main classes: model-based, signal-based, observer-based and parameter-estimation-based approaches of fault diagnosis [47,48,49,50]. The observer-based approach for fault detection/ estimation is used in this study, and more particularly it uses SMO for these tasks.
In the next subsection, the basis of the SMOs used for fault detection, isolation, and reconstruction, as shown in earlier works, is briefly discussed.

3.1. Fault Detection (FD) SMO Filtering

Fault diagnostics is performed using the SMO, that is, detection and isolation of faults by estimation of states/outputs for the MG application under consideration. In literal sense the isolation functions localization of the faulty sensors of the system; however, the fault estimation SMO is functioning for both the isolation/ estimation of faults. The first-order SMO as proposed by [23,24,30], gives the estimated states of the considered MG system.
x o ˙ = A c x o + B c u c + G o e o + G m ψ
where x o = x s * x h * R n c * 1 is the estimated state vector which augments estimated system states x s * R n * 1 and estimated stable filtered output states x h * R p * 1 . The matrix G o R n c * p c is the Luenberger gain of the output estimation error e o R p * 1 term in SMO, which also ensures the stability of the term ( A o = A c G o C c ) , whereas G m R n c * p c is the gain of the (switching) discontinuous term ( ψ ) of SMO. Both the gains G o and G m are required to be determined for SMO. The proposed form of the discontinuous term ( ψ ) is ( ψ = γ P o e o P o e o ) , where the constant gain factor ( γ ) is being appropriately chosen depending on the application under consideration. The proposed form of the discontinuous term gain G m = L T T T T , where the orthogonal matrix T R q * q is determined by using QR factorization. The matrices P o and L are sub parts of a PD Lyapunov matrix P > 0 , which is proposed to be of the form P = P 1 P 1 L L T P 1 T T P o T + L T P 1 L > 0 , where the matrices P R n c * n c , P 1 R n * n , P o R p * p , T R p * p , L R n * p are to be determined ([26]). The Luenberger gain G o is also determined from the Lyapunov matrix P, which is explained in detail in the next section.
The output estimated states of the system are given by
y o = C c x o
The output estimation error of the system is given by
e o = y o y c
and
A c = A 11 A 12 A 21 A 22 , x o = x s o x h o
As e o is not the actual output error, instead it is the (scaled) stable filtered output error, so it is suggested empirically to use the form of ψ in Equation (13) instead of its normalized form for the considered MG application.
ψ = γ * P o e o
Remark 2.
The observer as mentioned in Equation (14) is insensitive to faults (f) completely if:
1.
R a n k ( C c E c ) = q
2.
The System triple ( A c , E c , C c ) has the invariant zeros which lie in Left Half Plane (LHP) [30].

3.2. State Estimation Error SMO Filtering

The state estimation error is determined by taking the difference of SMO-based estimated states in Equation (14) and system states (being determined) from the mathematical model in (11), being discussed in the previous sections. The difference between Equations (11) and (14) gives the state estimation error.
e d = e = x o x c
The state estimation error SMO is
e ˙ = A c e E c f D c ξ G o e o + G m ψ
e ˙ = ( A c G o C c ) e E c f D c ξ + G m ψ
where A o = A c G o C c R n c * n c , e R n c * 1
The residual signal (i.e. the output estimation error) in terms of the augmented error vector e is given as
e o = r ( t ) = C c ( x c x o ) = C c e
Equation (20) is also a standard form of SMO [23], which is applied to the error system for augmented state error estimation. The state error is further also used for fault/disturbance estimations by attaining the sliding mode, where the state error surface is the considered sliding surface. The next section discusses the stability analysis and determination of gains of the fault detection/estimation SMO observers proposed for the MG system.

4. Determination of Sliding Mode Observer Gains through Stability Analysis

A lemma for the existence of the sliding mode is given, before the Lyapunov-based stability analysis of the proposed fault detection/ estimation SMOs for the considered MG system.
Lemma 1.
If σ ( e ) defines the sliding surface, then for the Lyapunov function
V = e T P e = e T P 1 / 2 P 1 / 2 e = P e 2
which implies that V = P e , σ ( e ) = P e , and σ 2 = P e 2 defines the distance from the sliding surface σ ( e ) = 0 . The sliding surface is reached if
d V d t = d V d σ d σ d t = σ T σ ˙ < 0
in the neighborhood of surface σ ( e ) = 0 and
σ ˙ = d σ d e o e ˙
Remark 3.
The Lemma discussed above can be studied in detail in [50], as it is the pivotal concept of SMO based FTC techniques used in this work.
Proposition 1.
If ( G o ) is SMO gain for the output error estimation term ( e o ) , G m the SMO gain of the discontinuous control term ( ψ ) is proposed to be of the form G m = L T T T T , the constant gain ( γ ) in the ( ψ ) term is constrained as ( γ η o E o α ) where ( η > 0 ) and P is a PD matrix, i.e., ( P > 0 ) of the form
P = P 1 P 1 L L T P 1 T T P o T + L T P 1 L > 0
which satisfies ( P A o + A o T P < 0 ) , then the estimation error e ( t ) stays bounded and hence asymptotically stable.
Proof. 
Let V define the Lyapunov function for an augmented error system. The stability of the equilibrium requires the Lyapunov function to be positive definite, and its time derivative to be negative or semi-negative-definite.
V ( e ) = e T P e
V ˙ ( e ) = e ˙ T P e + e T P e ˙
Using e ˙ from Equation (19)
V ˙ ( e ) = e T ( A c G o C c ) T P e + e T P ( A c G o C c ) e 2 e T P E c f 2 e T P D c ξ + 2 e T P G m ψ
Using the definition of G m as given in the statement of the proposition, and dropping the negative definite Lyapunov term because Z = A o T P + P A o < 0 , where A o = A c G o C c , the remaining terms are still negative, and the time derivative of the Lyapunov function becomes an inequality, which is always easier to handle in terms of parametric independence.
V ˙ ( e ) 2 e T P E c f 2 e T P D c ξ + 2 e T P G m ψ
Using the definitions of ( ψ , G m , P G o = C o T P o , P E c = C o T P o E o , and e T C o T = e o ) ,
V ˙ ( e ) 2 e o T C o T P o E o f 2 e o T C o T P o D o ξ 2 γ P o e o
Taking norm and upper bounds f < α , ξ < ξ o
V ˙ ( e ) 2 e o [ γ P o C o + P o C o E o α + P o C o D o ξ o ]
Let if γ P o C o + P o C o E o α + P o C o D o ξ o η o P o C o
γ η o E o α D o ξ o
Using the constraint on γ from Equation (30) in Equation (29)
V ˙ ( e ) 2 η o e o
Because the Thau inequality ([51]) for the Lyapunov equation for any PD matrix P o is defined for any vector x as [33,50].
x T P o 1 x λ m i n ( P o 1 ) x 2 2
The Lyapunov function in terms of the error function can be represented in the inequality form as: ( e T Z e > λ m i n ( Z ) e 2 ) . Using the Thau inequality and Equation (31), the inequality version of Equation (26) becomes:
V ˙ ( e ) λ m i n ( Z ) e 2 2 η o e o
which proves stability (i.e., negative definiteness) of time derivative of Lyapunov function.
Since x > λ m i n ( P 1 ) x 2 and can be seen in [33,50].
V ˙ ( e ) λ m i n ( Z ) e 2 2 η o λ m i n ( P o 1 ) P o e o
V ˙ ( e ) λ m i n ( Z ) λ m i n ( P 1 ) e 2 2 η o λ m i n ( P o 1 ) P o e o
The Equation (31) and its representation in terms of (34) show the asymptotic stability of the considered observer system.
Remark 4.
Proposition 1 is analytical proof, particularly for the considered microgrid system with faulty and perturbed sensors using the same steps as performed by Yuri and Edwards in Proposition 3.1 [50], whereas the Thau inequality and observer can be studied in more detail in [51].
Theorem 1.
Let Z = A o P + P A o < 0 such that P > 0 and the constant gain of the discontinuous switch term γ is constrained by γ η o E o α , then the augmented error system (state error and stable filtered output error) dynamics defined by Equation (19) remains bounded, such that the error magnitude remains bounded within the set X = { e | e 2 ξ o P D c λ m i n ( Z ) } , and the Lyapunov function in vector form gives the constraint in the form of LMI L i 1 = A o T P + P ( A o ) P D c D c T P 0 < 0 , where A o = A c G o C c , the LMI is further modified by general algebraic Ricatti equation with additional control parameters to give LMI L i 2 = A c T P + P A C c T F 1 1 C c Y Y T F 1 1 < 0 , and the iterative feasibility optimization or minimization of linear objective (i.e., trace) under LMI constraint L i 2 gives the SMO gain for output estimation error term to be G o = P 1 C T F 1 1 , such that Y is constrained to Y = C c T F 1 1 T > 0 , where ( F 1 > 0 ) .
Proof. 
Considering the Lyapunov function from Equation (26)
V ˙ ( e ) = e T ( A c G o C c ) P e + e T P ( A c G o C c ) e 2 e T P E c f 2 e T P D c ξ + 2 e T P G m ψ
Considering switch and fault terms from above equation, i.e., Equation (26)
Let
T s f = 2 e T P E c f + 2 e T P G m ψ
T s f = 2 e T P E c f + 2 e T P G m ( γ P o e o P o e o )
(Since P G o = C o T P o and P E c = C o T P o E o )
T s f = 2 e T C o T P o E o f 2 γ C o P o e o
Taking the norm, we use the Cauchy–Schwartz inequality and bounded fault with f α
T s f = 2 e o ( C o P o E o α + C o P o γ )
Because it is desired that the above term stays more negative T s f < 0 , for the negative definiteness of the time derivative of Lyapunov function in Equation (26), which is satisfied if
P o C o E o α + P o C o γ > η o P o c o γ η o E o α
Furthermore, since ( A o = A c G o C c ) , and if ( Z = A o P + P A o < 0 ) where ( P > 0 ) , so dropping the negative T s f terms, the Lyapunov function is still negative
V ˙ ( e ) e T ( A c G o C c ) P e + e T P ( A c G o C c ) e 2 e T P D c ξ
To show stability in terms of viable set form
V ˙ ( e ) λ m i n ( Z ) e 2 2 ξ o P D c e
V ˙ ( e ) e ( λ m i n ( Z ) e 2 2 ξ o P D c )
If ( e 2 > 2 ξ o P D c λ m i n ( Z ) ) , we obtain V ˙ < 0 , which ensures that the error magnitude remains bounded for the set
X = { e | e 2 ξ o P D c λ m i n ( Z ) }
Using Equation (36) and expressing in terms of the vector quadratic function, LMIs can be determined.
V ˙ ( e ) e T ξ T A o P + P A o P D c D c T P 0 e ξ
The constraint is true if the matrix in the vector quadratic function is negative definite.
L i 1 = A o P + P A o P D c D c T P 0 < 0
The LMI in Equation (40) can be solved for the trace minimization-based algorithm given by [26] to determine the optimized gains of the SMO. However, for the LMIs to be optimized with trace minimization without the H constraint for the determination of traditional SMO gains, it needs some modification. Applying the Schur complement to the LMI in Equation (40)
( A c G o C c ) T P + P ( A c G o C c ) 0
Because the LMI is not feasible in the actual form mentioned by L i 1 in Equation (40) with a trace-minimization-based optimization algorithm, some modifications are required. By adding/subtracting the term ( C c T F 1 1 C c ) and two more Lyapunov stable terms in above Lyapunov equation to form a balanced an algebraic Ricatti equation, i.e., ( Y F 1 Y T ) and ( P W P ) in the inequality, where ( Y = P G o )
A c T P + P A c Y C c C c T Y T + C c T F 1 1 C c C c F 1 1 C c + Y F 1 Y T + P W 1 P 0
A c T P + P A c C c T ( Y T F 1 1 C c ) Y ( C c F 1 Y T ) C c T F 1 1 C Y F 1 Y T + P W 1 P 0
Using/constraining C c F 1 Y T = 0 , which is equivalent to Y T F 1 1 C c = 0 gives the Luenberger SMO gain
G o = P 1 C c T F 1 1
Parameters such as the Lyapunov matrix P and matrix F are missing and need to be determined. Using the constraint C c F 1 Y T = 0 , Equation (43) becomes:
A c T P + P A c C c T F 1 1 C c + P W 1 P 0
Applying Schur complement on LMI in Equation (45), and the constraint F 1 1 C c = Y T
L i 2 = A c P + P A c C c T Y T P P W < 0
Remark 5.
The LMI L i 2 is optimized by using iterative convex feasibility optimization or alternatively trace minimization, as explained by an algorithm in ([26]), to determine the unknown parameters such as Lyapunov matrix P and matrix F, which are used to determine Luenberger gain G o of SMO. Moreover, the discontinuous term gain, G m = L T T T T as discussed in Proposition 1 can also be determined, as the parameter matrix L is also determined from the part P 12 of Lyapunov matrix P by the relation ( L = P 11 1 P 12 ) , where the orthogonal matrix T is can be determined by QR factorization.
The next subsection discusses the H enhanced trace minimization of LMIs to attain the robustness of the SMO for fault estimation with disturbance rejection.

H Optimized Robust Sliding Mode Observer Gains (Using LMIs)

H is a robust control criterion that may have several meanings with reference to context; however, in this work, it is desired that the fault detection and estimation be ensured with robustness against disturbances, whereas the H gain will ensure an upper bound of the disturbance voltage/current that will be suppressed to ensure the fault detection/estimation task. The method of incorporating H criteria may also have different approaches; however, this study uses a game-theoretic basis [46]. The criteria are incorporated in the Lyapunov function to derive the LMIs that are convex optimized (using the LMI-optimization toolbox in Matlab) for determination of the SMO gains.
Definition 2.
H based disturbance attenuation problem can be put into formulation as a two-player-zero-sum differential game with disturbance ξ and Luenberger observer gain G o are two players, respectively, where the SMO gains G o , G m are required to be designed to minimize the game /cost functional.
Lemma 2.
The maximum robustness attenuation problem
H = s u p [ ξ 0 , f = 0 ] r H 2 , [ 0 , t 1 ] ξ 2 , [ 0 , t 1 ] μ
for the system given in Equation (7), is satisfied if the cost functional
f H ( G o , G m , ξ ) = ( r H T r H μ ξ T ξ ) d t , ( f o r ) f = 0 .
is smaller than or equal to zero for any possible disturbance, where ( μ ) is the ( H ) parameter and ( r H = H r ( t ) = H e o = H C c e ) is the output estimation error. The problem is now viewed as a two-player-zero-sum differential game with the above defined cost function in Equation (48), where the minimizing player minimizes the functional through ( G o , G m ) and the maximizing player tries to maximize the function through ξ. The details can be seen in Problem 3.1 in [46].
Theorem 2.
Let Z = A o P + P A o < 0 , where P > 0 , and if ( H > 0 ) , then the augmented error system dynamics defined by Equation (19) remain asymptotically stable, and if the Lyapunov function in vector form is enhanced with H constraint, the disturbance attenuation gives the constraint in the form of LMI
L i 3 = A o P + P A o + C c T H T H C c P D c D c T P μ I < 0 ,
which is further modified by Riccati equation with additional control parameters to give constraint LMI
L i 4 = A c T P + P A c C c T Y T P D c D c T P μ 1 < 0 ,
and the iterative minimization of linear objective (i.e., trace) under LMI constraint L i 4 gives the Luenberger gain of output error estimation term as G o = P 1 C T F 1 , where Y = P G o > 0 is constrained to Y = C c T F > 0 and W = P D ¯ P .
Proof. 
Because V is a PD Lyapunov function, the negative definiteness of its derivative determines the stability of the system. If the H criterion is used to enhance the robustness of the traditional SMO, the criterion in terms of the norm is defined by
r H ( t ) μ ξ ( t )
Because the LMI solver constraints for feasible optimization require the residual to be of the form
r H ( t ) H r ( t )
where the output estimation error, i.e., residual is defined by
r ( t ) = e o = y o y c = C c e
and H is a scaling matrix. The Lyapunov function in inequality form in Equation (36) is added with the H criterion according to Lemma 2 as follows:
V ˙ + r H T r H μ ξ T ξ 0
V ˙ ( e ) e T [ ( A c G o C c ) P + P ( A c G o C c ) ] e + e T C c T H T H C c e 2 e T P D c ξ μ ξ T ξ
where
H T H = F = 0 F
The matrix F R p * p is a sub part of matrix H, which provides a control on the residual signal for stability purposes. However, as discussed above, in the definition of r H ( t ) , it is suggested by the LMI solver to meet the conditions of feasible optimization.
and if
H = H 1 H 2 , r H T r H = e T C c T H T H C c e = e T C c F C c e
where H R n c * p c , H 1 R n * p , H 2 = F R p * p
V ˙ ( e ) e T ξ T A o P + P A o + C c T F C a P D c D c T P μ I e ξ
The LMI, which can be solved for feasibility or trace minimization-based optimization, to determine the optimized parameters of the matrix is
L i 3 = A o P + P A o + C c T F C c P D c D c T P μ I < 0
The Schur complement can be applied to convert any bi-linearity to linearity. Applying Schur complement on LMI in Equation (56)
( A c G o C c ) T P + P ( A c G o C c ) + C c T F 1 C c + P D ( μ 1 I ) D T P 0
Using ( Y = P G o ) and adding subtracting the term ( Y F 1 Y T )
A c T P + P A C c ( Y C c T F ) Y T ( C c T Y F 1 ) + Y F 1 Y T + P D ¯ P 0
where D ¯ = D c ( μ I ) D c T
Setting ( Y C c T F = 0 ) which is equivalent to ( C c T Y F 1 = 0 ) gives
G o = P 1 C c T F 1
and the inequality reduces to
A c T P + P A c Y F 1 Y T + W 0
where W = P D ¯ P
Using Schur complement on inequality and C T = Y F 1
L i 4 = A c T P + P A c C T Y T P D c D c T P μ I < 0
Remark 6.
Trace minimization based optimization of Ricatti equation-motivated modification of LMIs, for determination of SMO gains is presented by [26], and H based feasibility optimization of LMIs for determination of SMO gains is used by [33]. However, this work combines the application of both, i.e., to carry out the feasibility or the linear objective (i.e., trace) minimization-based iterative convex optimization of Ricatti equation-based modification of H enhanced LMIs, which can work or may be tested on generally any system, and particularly here being applied on MG under consideration.

5. Reduced Order Sliding Motion

The error system augments state/ stable filtered output error vectors that is,
e s ˙ e o ˙ = A 11 A 12 A 21 A 22 e s e o G 1 G 2 e o 0 E o f 0 D o ξ L T T T T ψ
where e = e s e o ;   G o = G 1 G 2 ;   E c = 0 E o ;   G m = L T T T T
The generalized dimensions of matrices and vectors are e o R p * 1 , e s R n * 1 , ξ R q * 1 , ψ R q * 1 , f R q * 1 , A 11 R n * n , A 12 R n * p , A 21 R p * n , A 22 R p * p , E o R q * q , D o R q * q , G o R n c * p , G 1 R n * p , G 2 R p * p , L R n * p , T R q * q whereas for the MG system considered p = 4 ,   n = 6 ,   m c = m = 2 ,   m = 2 ,   m c = m = 2 ,   n c = n + p = 10 ,   p c = p = 4 ,   q c = q = 4 .
A transformation is used to make the state error part of the error system least dependent on faults to retain the sliding motion, along with the attainment of reduced-order sliding motion for fault estimation. Ideally, the fault term in the state error part of the system and the state error term in the output error part should be eliminated; however, such a transformation cannot be designed. Thus, practically and very strictly, the sliding motion may not be retained owing to faults and disturbances, but the transformation not only provides control over the state error part through the gain matrix ( L ) but also reduces the effect of the state error part in the fault term to a negligible order in terms of magnitude for the considered MG system. However, the reduced-order state estimation error increases with the passage of time, which can be controlled by scalar ( δ ) and ( β ) parameters. Applying Transformation
T L = I n p L 0 T
where its inverse, which explains the form of the discontinuous term gain of the SMO, is given by
T L 1 = I n p L T T 0 T T
The transformation obtains the system to attain a reduced-order sliding surface along with fault estimation. The transformation and its inverse also explain the sense of the considered form of the switching term gain G m in the SMO. The orthogonal matrix T was determined using the QR transformation method.
A ¯ = T L A T L 1 = A ¯ 11 A ¯ 12 A ¯ 21 A ¯ 22
where
A ¯ 11 = A 11 + L A 21 , A ¯ 21 = T A 21
A ¯ 12 = L A 22 ) T T ( A 11 + L A 21 ) L T T + ( A 12
A ¯ 22 = T A 22 T T T A 21 L T T
T L G o = G 1 + L G 2 T G 2 , T L E c = L E o T E o ,
T L D c = L D o T D o , T L G m = 0 I , C a T L = 0 T
e s ˙ e o ˙ = A ¯ 11 A ¯ 12 A ¯ 21 A ¯ 22 e s e o G 1 + L G 2 T G 2 0 e o L E o T E o f L D o T D o ξ 0 I ψ
e s ˙ = A ¯ 11 e s + A ¯ 12 e o ( G 1 + L G 2 ) e o L E o f L D o ξ
e o ˙ = A ¯ 21 e s + A ¯ 22 e o + T G 2 e o T E o f T D o ξ + ψ
A reduced-order error system is needed for fault and disturbance estimation if the observer mechanism remains stable along with reachability of the sliding surface in the finite time. The scalar gain ( γ ) to be used with he output estimation error term ( e o ) in the discontinuous control law defined by term ( ψ ) is also determined from the analytical work of the stability analysis.

5.1. Reachability and Stability Analysis

The complete stability analysis of SMOs used for fault detection and estimation is shown in terms of the estimation error as a sliding surface. The error vector comprises of the state estimation error and the stable filtered output estimation error in the augmented vector form. A transformation is used in the above section for reduced-order error dynamics, which is particularly needed for fault and disturbance estimations. Therefore, considering the state error stability from the earlier section in Proposition 1, the Lyapunov stability of only the stable filtered output estimation error part is performed again to show the reachability of the sliding surface and hence ensures the estimation of the fault in finite time. The reachability is directly related to stability, and its analytical work is also used to determine the sliding surface reachability time. Before detailed analytical work, a Lemma for the stability and reachability conditions is provided.
Lemma 3.
If the sliding mode for the surface σ ( e ) is attained in finite time, that is, σ ( e ) = 0 , then the time derivative of the Lyapunov function d V / d t may follow the inequality
d V d t k ( V ) z
must be bounded more strongly, somewhat away from zero. That is, the attraction to the sliding mode will only be asymptotic, if it vanishes too quickly. The solution of inequality in Equation (65) is given by
2 V ( t ) k V ( t ) + 2 V 0
Because V 0 , then the inequality in Equation (71) states that V must reach V = 0 in finite time. In addition, because V is proportional to the Euclidean norm of the sliding surface σ ( e o ) , that is, σ ( e o ) , the rate of reaching the sliding surface is bounded away from zero [24,50].
Theorem 3.
If the augmented state estimation error system defined by Equation (66) is transformed by matrix T L = I n p L 0 T to induce reduced order sliding motion on estimation SMO,
and if the constraint on Lyapunov design matrix is,
Q = P o T ( A 22 G 2 ) + ( A 22 T G 2 T ) T T P o < 0 ) ,
and if σ ( e o ) = { e o : C e o = 0 } governs the reduced order sliding motion,
such that the magnitudes of fault/ disturbance are bounded, i.e.,
( f < α ) and ( ξ ξ o ) ,
and if the gain factor ( γ ) for output error estimation term is bounded by
γ T A ¯ 21 e s T E o α T D o ξ o + η ,
then the fault detector/ estimator SMOs utilized for MG system are ensured to be stable in terms of Lyapunov criteria along with the reachability of the sliding motion in the finite time. The finite time required to reach/hit the sliding surface is given by:
T R e o T P o e o λ m i n ( P o 1 )
Proof. 
Let the Lyapunov function for the error system be V e , where the error vector is augmented by the state estimation error and stable filtered output error for the considered system. Thus, the complete Lyapunov function can be represented as:
V e = V e s + V e o = e s T P 1 e s + e o T P o e o
Taking the time derivative of the Lyapunov function to prove the stability of the error system based on its negative definiteness.
V e ˙ = e s T ˙ P 1 e s + e s T P 1 e s ˙ + e o T ˙ P o e o + e o T P o e o ˙
Substituting the values for e s ˙ from Equation (63)
V e s ˙ = e s T A ¯ 11 T P 1 e s + e o T A ¯ 12 T P 1 e s f T E o T L T P 1 e s ξ T D o T L T P 1 e s e o T G 1 T P 1 e s e o T G 2 T L T P 1 e s + e s T P 1 A ¯ 11 e s + e s T P 1 A ¯ 12 e o e s T P 1 L E o f e s T P 1 L D o ξ e s T P 1 G 1 e o e s T P 1 L G 2 e o
The stability of the sliding motion in reduced-order is proved by considering the Lyapunov function on the output part of the error vector only.
V e o ˙ = e o ˙ T P e o + e o T P e o ˙
V e o ˙ = e s T A ¯ 21 T T T P o e o + e o T A ¯ 22 T T T P o e o + e o T G 2 T T T P o e o f T E o T T T P e o ξ T D o T T T P o e o + ψ T P o e o + e o T P o T A ¯ 22 e o + e o T P o T G 2 e o + e o T P o T G 2 e o + e o T P o T E o f e o T P o T D o ξ + e o T P o ψ
V e o ˙ = 2 e o T P o T T A ¯ 22 e o + 2 e o T P o T T G 2 e o + 2 e o T P o T T G 2 e o + 2 e o T P o E o f 2 e o T P o D o ξ + 2 e o T P o ψ
V e o ˙ = e o T [ P o T ( A ¯ 22 G 2 ) + ( A ¯ 22 T G 2 T ) T T P o ] e o + 2 e o T P o T A ¯ 21 e s 2 e o P o T E o f 2 e o T P o T D o ξ + 2 e o T P o ψ
If designed suitably
Q = P o T ( A 22 G 2 ) + ( A 22 T G 2 T ) T T P o < 0
Dropping the negative definite term, the inequality version of Equation (78) is
V e o ˙ < 2 e o T P o T A ¯ 21 e s 2 e o P o T E o f 2 e o T P o T D o ξ + 2 e o T P o ψ
Taking norm of above equation, and using the Cauchy–Schwartz inequality
V e o ˙ 2 e o P o T A ¯ 21 e s 2 e o P o T E o f 2 e o P o T D o ξ 2 γ P o e o
V e o ˙ 2 P o e o [ T A 21 e s T E o f T D o ξ o γ ]
Since Furthermore, since magnitudes of fault and disturbance are bounded: f ( t ) < α and ξ ( t ) ξ o <
V e o ˙ 2 P o e o [ T A ¯ 21 e s T E o α T D o ξ o γ ]
The derivative of Lyapunov function is negative definite if
γ > T A ¯ 21 e s T E o α T D o ξ o
where η is a positive constant that ensures Lyapunov stability in terms of V e o ˙ . Using the bound for γ from Equation (84) in Equation (83),
V e o ˙ 2 η P o e o
If P is a positive definite matrix, then x 2 2 > λ m i n ( P 1 ) x 2 2 ([51])
e o 2 = e o T e o = e o T P o 1 / 2 P o 1 P o 1 / 2 e o λ m i n ( P o 1 e o 2
e o λ m i n ( P o 1 ) V e o
where P o 1 / 2 e o = V e o .
Using Lemma 3, Equation (85) shows the stability and finite-time reachability of the sliding motion for the sliding-mode observers. The finite-time reachability governs and ensures the real-time operation of the proposed observers for fault detection and estimation, that is, the sliding surface is ensured to be reached/hit in a finite time. Using Equations (79), (80), (84), and (85) in the Lyapunov function in Equation (78), in the inequality form yields:
V e o ˙ λ m i n ( Q ) e o 2 2 η P o e o
Using Equation (86)
V e o ˙ λ m i n ( Q ) λ m i n ( P o 1 ) V e o 2 η λ m i n ( P o 1 ) V e o
V e o ˙ λ m i n ( Q ) λ m i n ( P o 1 ) V e o 2 η λ m i n ( P o 1 ) V e o
For the DE of the form ( x ˙ ( t ) a x ( t ) b x ( t ) ) ([33]), the time required for x ( t ) to move from x 0 = x ( t = 0 ) to x 1 = x ( t = 1 ) is
T R 2 a l n ( a x o + b a x 1 + b ) 2 b ( x o x 1 )
The reachability time ( T R ) is given by
T R 1 η ( V o ( e o ( t o ) ) V o ( e o ( t o + T R ) ) λ m i n ( P o 1 ) ) e o T P o e o λ m i n ( P o 1 ) .
where
V o ( e o ( t o + T R ) ) = e o T ( t o + T R ) P o e o ( t o + T R ) = 0
Because the stability and finite-time reachability of the sliding-mode observers is shown and the gain of the output error injection term is determined, the estimation of faults and disturbances can proceed with the results determined in previous sections.

5.2. Reconstructed Fault and Estimated Disturbance

This subsection describes the analytical work of faults and disturbance reconstruction based on the results obtained in the previous sections.
Corollary 1.
Considering the MG system/ observer in Equations (11) and (14) respectively, if the SMO gains are determined by solving the LMI optimization problem described in Theorem 2, the transformed error system is defined by Equations (68) and (69), and if (γ) being constant gain of output error estimation term is constrained to satisfy (84) to ensure the stability of observers, i.e.,
γ T A ¯ 21 e s T E o α T D o ξ o + η ,
then sliding surface attained in reduced order gives the estimated sensor fault ( f * ) / disturbance ( ξ * ) to be:
f * = β E o 1 T 1 ψ e q
ξ * = f f * E o 1 T 1 A 21 e s
where β serves as a scaling constant.
Proof. 
When observers enter the sliding mode, that is, the sliding surface is reached, the condition in terms of the error surface ideally is ( C e = 0 ) ; however, for the considered perturbed and faulty sensor MG systems, the sliding mode/surface is reached only for the output error term ( e o = e o ˙ = 0 ) and not for the state estimation error term ( e s ). However, if there is no disturbance term, then ( e s 0 ; e s ˙ 0 ) are also approached in finite time.
e s ˙ = A ¯ 11 e s L E o f L D o ξ + T ψ e q
0 = A ¯ 11 e s T E o f T D o ξ + T ψ e q
From above Equation (92)
f = E o 1 T 1 A 21 e s E o 1 D o ξ + E o 1 T 1 ψ e q
Using the fault estimate Equation (93) in the state error Equation (94) yields the reduced-order state error equation
e s ˙ = ( A ¯ 11 L T 1 A ¯ 21 ) e s L T 1 ψ e q
This DE needs to be solved linearly at every instant of time to be back-substituted in the fault estimation Equation (93) to obtain a closer fault estimate. As evident from Equation (94), the solution to DE also requires a run time evaluation of factor ( ψ e q ) , because the aforementioned approach is necessary for disturbance estimation. Considering the terms in the fault estimate Equation (93), because the state error term e s in the fault equation cannot be removed or isolated by the transformation ( T L ) applied to system (62), but its magnitude becomes insignificant for the considered MG application, whereas the disturbance term in the equation is known to be an undesired and diverging term, the reconstructed fault should be approximated equivalently by ( ψ ) term only, which is based on the output estimation error ( e o ) term, design matrices ( E o , T ) , and gamma ( γ ) parameter.
So from Equation (93) the reconstructed fault is
f * = E o 1 T 1 ψ e q
The equivalent switched control in sliding mode is ψ e q is
ψ e q = γ P o e o P o e o + δ
where the term ψ e q is meant to maintain the sliding motion in reduced-order motion in the presence of disturbances and faults, and δ is a constant control parameter for this purpose. □
Remark 7. ( γ ) is constant gain factor of discontinuous control term ( ψ e q ) whereas ( δ ) is its additional constant control parameter. The ( ψ e q ) term controls the gradually increasing reduced-order state error ( e s ) and the value of ( γ ) being in a closed loop, which, if not controlled properly, causes a rapidly increasing reduced-order state error ( e s ) and hence the total divergence of the whole process. It is empirically suggested that, for the considered MG application, the reconstructed fault is
f * = β E o 1 T 1 ψ e q
where β is a constant magnitude scaling parameter that is required to adjust the attenuation achieved by stable filtering on the output error term.
Because E o = D o E o 1 D o = I d e n t i t y , using Equation (97) in Equation (93), the estimated disturbance in the system is given by
ξ = f f * E o 1 T 1 A 21 e s
Remark 8.
The values of parameters (δ, β, γ) are considered in this paper intuitively by hit and trial, however setting the parameters sub optimally can be addressed in future work, which would also provide generic applicability of observer based fault-tolerant approach for many other applications with different fault/disturbance magnitudes and sensitivity requirements.
Remark 9.
The presence of disturbances and faults result that the sliding motion is not attained in a strict sense; however the effect of disturbances and faults on fault estimation can be minimized using ( H ) norm as presented in Theorem 2 which is ratio of L 2 norms of residual ( r ( t ) = e o ( t ) ) and disturbance (ξ), respectively.
Remark 10.
The state error term ( e s ) in output error term and fault term in state error term are undesired in terms of fault estimation. Although the transformation ( T L ) cannot ideally achieve the desired elimination of undesired terms, it provides control through the gain matrix ( L ) , which is determined by the LMIs.
Remark 11.
Moreover the transformation ( T L ) reduces the effect of state error part ( e s ) in estimated fault to quite negligible order in terms of magnitude for the considered system which is almost zero in the case of micro-grids. The results were verified through simulations.

6. Fault Tolerant Control

The estimated faults in C.T/P.T are used for corrections in the faulty sensors’ (transformers) readings to avoid the wrong current/voltage readings to calculate the instantaneous powers and correspondingly generation of wrong PWM signals. The corrected dq-currents/voltages of current/potential transformers are given by the vector
Y s c = Y s E s * f * D s * ξ *
The corrected sensor voltages/currents are given to voltage-frequency control block for proper tracking of reference reactive and real powers.

P ω and Q V Control Scheme

The control block mainly utilizes PI control-based current/voltage controllers and uses droop control and a PLL block to regulate real and complex power in relation to frequency and quadrature voltage, respectively. In other words it is said that, the same P f and Q V control mechanisms are enhanced as FTC by using observer-based fault estimates for the correction of sensor outputs. The control block is not discussed in detail except for a brief continuation, as it is taken from the work of [6], from which the model of the MG is considered in Section 2. All the system parameters and PI constants are considered to be the same as those given in the aforementioned paper and detailed PhD thesis.
DROOP CONTROL
The voltage and frequency references for the standalone mode of operation of the microgrids must be internally generated.
ω * = ω n m P r
V * = V o q n n Q c
where P r = 3 / 2 * ( V o q I o q + V o d I o d ) and Q c = 3 / 2 * ( V o q I o d ˘ V o d I o q ) Because the values are accepted directly from the grid and low-pass filtered output power is processed, the high-frequency current and voltage spikes must be neglected.
P = ω c , P L L S + ω c , P L L P r , Q = ω c , P L L S + ω c , P L L Q r , where ω c , P L L = 7853.98 rad/s
m = ω 1 ω 2 P 1 P 2 ; n = V 1 V 2 Q 1 Q 2
The value of m = n = 1 / 1000 is used for the simulation, which can be varied for different grid ratings and other applications.
PLL BLOCK
The block computes the grid angle by using the phase-locked loop (PLL) at the instant that is then used for all dq conversions, SVPWM block, and other required computations.
V o d f ˙ = ω c , P L L V o d ω c , P L L V o d f
θ ˙ = ω P L L θ = ω P L L d t
For f = 60 Hz ω P L L = 377 rad/s and ω c , P L L = 7853.98 ,
w P L L = 377 k p P L L V o d f k i P L L V o d f d t
The values of constants and gains used in simulations are: k p = 0.25 , k i = 2 , w c = 377 rad/s.
VOLTAGE CONTROLLER
The PI-based control equations managing integral-based control serve as a voltage controller and generate the corresponding references of the dq components of the current.
i i d * = k i v d ( w P L L w * ) d t + k p v d ( w P L L w * )
i i q * = k i v q ( V o q * V o q ) d t + k p v q ( ( V o q * V o q ) )
The values of constants and gains used in simulations are: k i v d = 0.25 ,   k i v q = 0.25 ,   k p v d = 0.5 ,   k p v q = 0.5 .
CURRENT CONTROLLER
The block uses the desired current values generated by the above block and generates the corresponding voltage values.
V i d * = w n L f i i q + k i c d ( i i d * i i d ) d t + k p c d ( i i d * i i d )
V i q * = w n L f i i d + k i c q ( i i d * i i d ) d t + k p c q ( i i q * i i q )
The values of constants and gains used in simulations are w n = 377 ,   k p c d = 1 , k p c q = 1 , k i c d = 100 ,   k i c q = 100 .
SVPWM CONTROL
The SVPWM control uses a 10,000 Hz frequency for the considered MG operation. The desired voltage values generated by current controllers are transformed back to a three-phase voltage representation and provided to the SVPWM block, which manages the switching sequences and times of inverter switches to generate the desired values of current and voltage by VSC. These current/voltage values correspondingly manage the real power vs. frequency and complex power vs. voltage sags. The converter model used was a detailed model provided in MATLAB (Simulink) instead of the average model.
The complete procedure for the proposed observer-based FTC approach is explained in Algorithm 1.
Algorithm 1. Procedure Algorithm
ine Inputs: w, u, ξ , f
Outputs: Matrices P, G m , G o , y s , x s , x o , e, γ , Ψ e q , ξ * , f * , e s , y o , e o
while (Simulation Time)
START:
         1: Linearized system model in Equation (2) is d-q-0 transformed using Equation (3);
         2: Acquisition of grid/converter voltages (u, w) from microgrid in Figure 2 are given to system model ( x s ˙ = A s x s + B s u + B g w ), Equation (6)
         3: Sensor/C.T/P.T faults/disturbances ( f ( t ) = f o * s i n ( ω 1 t + ϕ 1 ) ), ( ξ o * s i n ( ω 2 t + ϕ 2 ) ) are generated to be added in system output ( y s = C s x s + E s f + D s ξ ), Equation (7)
         4: Faulty system output is stable filtered using ( x h ˙ = A h y s A h x h )  Equation (8);
         5: Augment system states with stable filtered outputs, Equations (10) and (11);
         6: Determine the SMO gains G o , G m ) using convex (feasibility and trace minimization) optimization performed on system parameters according to methods explained in Theorems 1 and 2;
< (i) Perform convex feasibility optimization on the LMI constrained with Lyapunov stability given in Equation (46) (to be used for trace optimization), (ii) Lyapunov stability and H constrained LMI in Equation (56) (to be feasibility optimized), and (iii) Ricatti equation-based modification of H constrained LMI in Equation (61) (to be feasibility optimized), to find the respective optimal gains G o * , G m . >
         7: Augmented system state is given to state estimator (SMO) ( x o ˙ = A c x o + B c u c + G o e o + G m ψ ), Equation (14);
         8: Difference between system and observer states/outputs ( e d = e = x o x c ) as defined in Equation (18) gives Augmented state estimation error;
         9: Augmented system state is given to state error estimator SMO ( e ˙ = ( A c G o C c ) e E c f D c ξ + G m ψ ), Equation (19) or (20);
         10: Stable filtered output estimation error ( e o = y o y c ) in Equation (16) is fed to state estimator SMO in step 6 and the state error estimator SMO in step 9 to attain the sliding mode;
         11: Determine the gains ( γ > T A ¯ 21 e s T E o α T D o ξ o ) and the suitable value of Delta, Equation (84);
         12: Determine the gain ψ e q = γ P o e o P o e o + δ as in Equation (96);
         13: Determine the gain reduced order state error e s )) by Simulink based numerical solution of ( e s ˙ = ( A ¯ 11 L T 1 A ¯ 21 ) e s L T 1 ψ e q ) as in Equation (94);
         14: Compute the estimated fault by ( f * = β E o 1 T 1 ψ e q ), Equation (97)
         15: Compute the estimated disturbance (unknown input) by ( ξ = f f * E o 1 T 1 A 21 e s ), Equation (98);
         16: The erroneous sensor outputs are corrected using estimated faults/disturbances i.e. ( Y s c = Y s E s * f * D s * ξ * ), Equation (99).
         17: The corrected sensor output values are fed to PI/ Droop-based voltage/current/complex power/real power control block, Figure 3;
         18: Repeat Step 6 and on-wards for the feasibility optimization of H constrained LMI in Equation (56) and the Ricatti equation based modification of H constrained LMI in Equation (61).
END (while)

7. Results and Discussions

According to the literature [46,47,48], the most commonly occurring faults in sensors are impulsive (incipient), occurring at regular intervals (intermittent), linearly increasing, constant measurement error, and random faults. However, for current and potential transformers in saturation, the faults occurring are mainly sinusoidal magnitude, phase, and harmonic faults. The simulations are performed for state/output estimation, stable filtered state/output estimation, state estimation error, stable filtered output estimation error, state estimation filter error, reduced order state estimation error, disturbance estimation, fault reconstruction, and fault tolerant control (FTC) performance. All these are performed for six fault cases: constant, ramp, sawtooth, square sinusoidal, and random types of additive faults along with additive sinusoidal disturbances of the first, second, and third harmonics. The same types of faults are considered for both voltage and current sensors simultaneously with sinusoidal additive disturbances of different and same frequencies; however, the proposed FD and FTC mechanisms are quite robust against the cross options as well. To avoid unnecessary details and length of paper, the reduced order errors and stable filtered states/outputs are not given for any case, whereas the FTC performance is given for worst-case sinusoidal faults (among the considered ones). The behavior of the system and the results for all simulations were consistent, and no statistical analysis was required for a deterministic system and simulation platform.
After detailed testing of the system, the results are quite good in terms of accuracy, except for the fault and disturbance signals in the near frequency range of 60o, which produced scaled current fault estimation in some cases and delayed current fault estimation results in some cases; however, the voltage faults are accurately estimated in all cases. Moreover, there is an occurrence of time delay problem, which needs discrete time compensation for the phase. Some fault estimation errors and less accurate FTC performance of current is due to the very reason. However the problem is not corrected here, instead will be considered for future works.
A Simulink-based detailed three-phased inverter model was considered in the simulation. The grid/ DC source voltage both operate at 600 V. the phase of the grid voltage is used for PLL block and all abc-dq transformations; SVPWM operates at a frequency of 10,000 Hz with sampling time of 0.0001. The discontinuities are caused by greater sampling times, which can be reduced to improve accuracy at the cost of lesser ability of online working due to the increased response time. Because the continuous time simulations move at very low speeds, which are not viable for real-time online performance; therefore, fixed-step solvers are used for simulation in Simulink (MATLAB) with single task handling to avoid complexities with very minor compromises on accuracy.
Regarding some other simulation constants, since the residual magnitude if considered peak to peak is 0.2 for the (dq) currents/voltages for the considered time of simulation; whereas the ( m a x ( ξ ) = 6.32 ; so the H norm μ in this case has the numerical value 0.01.
The value of η is a small positive constant considered η = 10 to ensure the constraint in inequality (90), whereas the upper bounds for the current and voltage are α = 1 A / 10 V , respectively. However, for the worst cases and increased magnitudes, it is considered up to 10 A/100 V for I/V, whereas the upper bounds of ξ are normally considered as 0.2 A/2 V for I/V, and for the worst case with increased magnitudes in simulations 3 A/10 V for I/V. The microgrid system model parameters are listed in Table 1, and the controller parameters are listed in Table 2.

7.1. Case I—Sinusoidal Faults and Disturbances/Worst Case Scenario (among the Considered Ones)

Simulations were performed for randomly considered different frequencies and phases for both faults and disturbances and quite high fault/disturbance magnitudes, and they provided quite good performance. The voltage and current fault magnitudes are 100 V and 10 A, the frequencies are 120 Hz and 180 Hz, and the phases are 75 and 240 , respectively, whereas the disturbance magnitudes are 10 V and 3 A, the frequencies are 75 Hz and 300 Hz, and the phases are 45 and 310 , respectively. The results are shown for the state/output estimation in Figure 4, Figure 5 and Figure 6, the state/output estimation errors in Figure 7, which are performed with modified H constrained SMO gains only. The results in Figure 8 presents voltage fault reconstruction, whereas the current fault reconstruction is presented in Figure 9, voltage fault estimation error in Figure 10, current fault estimation error in Figure 11, current disturbance estimation in Figure 12, FTC performance for q-component of output current in Figure 13, FTC performance for d-component of output voltage in Figure 14 and current and voltage faults reconstruction for ramp/linearly increasing faults in Figure 15. The above mentioned results are compared for SMO gains which optimized for Trace minimization algorithm proposed by Edwards in [25], feasibility optimization of H-infinity enhanced SMO, and Riccatti equation based modification of H-infinity enhanced feasibility optimized gains. Voltage fault reconstruction has a deformation of the lesser grade; however, the current fault reconstruction and fault estimation error increases due to time delay issue, which can be dealt with separately in future works.
An important Remark on Results: The trace minimization of Ricatti-equation-based modification of H -enhanced LMI may not be feasible for some forms of LMIs; however the feasibility optimization of the LMI works well nearly in all cases. The minimization of trace is a linear objective for convex optimization, which can work well on some forms of LMIs or some systems depending upon its parameters. However, the trace minimization of Ricatti equation based modification of LMIs not enhanced with H criterion presented by Edwards in [25] also works well nearly in all cases, but its not giving the accuracy of results in the considered microgrid system, which is included in comparisons. This particular aspect still requires much of rigorous testing to comment in a more deterministic way. The feasibility optimized H -enhanced-LMIs also have a phase lag effect particularly in current signal, which needs to be compensated, but compensation of phase is not considered in this study.
Figure 4, Figure 5 and Figure 6 are showing the state/output estimation using SMO gains determined by the proposed Riccati equation based modification of H enhanced gains. Instead of all states, the sate estimation performance results are given for estimation of input current (d-component), output current (d-component), Iod, and output voltage )q-component). The results are apparently quite close in terms of estimation except some scaling which can be adjusted further as well.
Figure 7 is showing the state/output estimation error using the modified method of determination of SMO gains. The state estimation error can be considered as a performance index for working of state/output estimator SMO, which also acts as fault detection filter. In the fault magnitude of 10 A/100 V the index stays within a controlled range effectively 1 V/0.5 A, with a negligible small mean value.
Figure 8 and Figure 9 are showing the voltage and current faults estimation/reconstruction using SMOs using the gains optimized with Trace optimization, feasibility optimization of H constrained LMIs, and trace/feasibility optimization of Ricatti equation based modification of H constrained LMIs. The results of the later two as given in proposed technique in Theorem 2 are quite better in terms of accuracy of reconstruction.
Figure 10 and Figure 11 are showing voltage and current fault estimation errors, which can be considered as performance index of working of fault estimations and ultimately the FTC, whose accuracy depends. The results are compared for voltage/current fault estimations performed with SMOs using the gains optimized with Trace optimization, feasibility optimization of H constrained LMIs and trace/feasibility optimization of Ricatti equation based modification of H constrained LMIs. The results of the later two are nearly comparable.
Figure 12 is showing the current disturbances estimation/reconstruction using SMO using gains with the above three methods to lead towards a comparison analysis with the proposed work. The results for all three methods are in compromise and no one method can apparently be said to be better than the others. This accuracy is achieved due to the modified and corrected disturbance estimation procedure mentioned in Corollary 1.
Figure 13 and Figure 14 give the FTC analysis shown for q-components of output voltage. The faulty q-components being corrected by SM observer based fault estimation. The fault estimate and correction using the above mentioned techniques of SMO is taken as comparison, because the accuracy of fault/disturbance estimation ensures the reliable FTC working. The FTC achieved through proposed SMO observer is better as shown by graphical result.

7.2. Case II—Ramp Faults and Sinusoidal Disturbances

Simulations are performed in this case for current/voltage injected linearly increasing faults with gradients 4 and 1, respectively, and third-harmonic sinusoidal disturbance, that is, 180 Hz for both V/I. The results for current/voltage fault reconstruction in Figure 15 are given.

8. Conclusions and Future Directions

The study uses SMO theory for detection, isolation, estimation, and correction of faults as a scheme of fault-tolerant control, for faults occurring in sensors (current/potential transformers) mounted on microgrid, by using aVSI-based microgrid model. The work is generally applicable for a wide range of sensor/actuator faults and systems. The estimation and correction of faults using SMOs is like providing software based sensors (transformers) replacing the real ones at the time of occurrence of faults. The SMO-based estimated faults are used to correct the faulty readings of the dq-currents/voltages, to be supplied to PI-based conventional voltage-frequency control block, to determine the actual instantaneous values of reactive and real powers, and hence providing the correct SVPWM pulses provided to VSI. The gains of SMOs are determined using convex optimization of Lyapunov stability ensured LMIs. The comparisons for results produced with SMO gains determined by proposed technique to those determined by earlier/base works used in this study, are also presented, which shows the improvements. The comparison of ordinary voltage-frequency control without the proposed SMO-observer-based corrective mechanism is also shown. The state and fault estimation errors are considered as indices to show the effectiveness of the proposed method. The finite-time reachability of detection/estimation SMOs are also presented to show the real-time applicability of the study in VSI-based microgrids. The work sums up various previous works along with modifications in fault estimation and LMI optimization algorithm for determination of SMO gains. The method is transformable to provide fault detection/estimation/tolerance in various types of systems by selecting the suitable parameters.The work is intended to be enhanced for networked microgrids using the approaches of distributed controls while managing the optimal power flow.

Author Contributions

Conceptualization, E.S., A.U.K. and A.W.; Data curation, E.S., A.S., G.H. and A.W.; Formal analysis, E.S., M.I., A.S., G.H., A.W., F.R.A. and Z.U.; Funding acquisition, G.H., A.W. and Z.U.; Investigation, E.S., M.I., A.S., A.W. and Zahid Ullah; Methodology, E.S., A.U.K., A.S., G.H. and A.W.; Project administration, A.U.K., G.H., A.W. and F.R.A.; Resources, E.S., A.U.K., M.I., A.W. and F.R.A.; Software, E.S., M.I., A.W. and F.R.A.; Supervision, A.U.K., G.H., A.W. and F.R.A.; Validation, E.S., M.I., A.W. and F.R.A.; Visualization, E.S., A.W., F.R.A. and Z.U.; Writing—original draft, E.S., A.S. and A.W.; Writing—review & editing, M.I., A.S., A.W., F.R.A. and Z.U. All authors have read and agreed to the published version of the manuscript.

Funding

The APC of this manuscript is funded by Taif University Researchers Supporting Project Number (TURSP-2020/331), Taif University, Taif, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data including Simulink-based simlation files, Matlab based optimization m-files, supplementary result file and some necessarily required documents will be uploaded as the article gets published.

Acknowledgments

The authors would like to acknowledge the support from Taif University Researchers Supporting Project Number (TURSP-2020/331), Taif University, Taif, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Sliding Mode ObserverSMOSliding Mode ControlSMC
Robust Sliding Mode ObserverRSMOLoad Frequency controlLFC
Infemum/Supremuminf/supDifferential EquationDE
Fault Detection/DiagnosticsFDFault EstimationFE
Un-interrupted Power SupplyUPSVoltage Source ConverterVSC
Current/Potential TransformerC.T/P.TPhase Locked LoopPLL
Inductor-Capacitor-InductorLCLWind Energy SystemWES
MicrogridMGLinear Quadratic GaussianLQG
Linear Quadratic RegulatorLQRGenetic AlgorithmGA
Particle Swarm OptimizationPSOPositive DefinitePD
Model Reference Adaptive ControlMRACPhoto-voltaicPV
Linear Matrix InequalityLMILinear Parameter VaryingLPV
Model Predictive ControlMPCFault-Tolerant ControlFTC
Left/Right Half PlaneLHP/RHPProportional-IntegralPI
Proportional-Integral-DifferentiatorPIDDistributed GeneratorsDGs
Phase Locked Loop Frequency ω p L L Fault/Disturbancef/ ξ
Discontinuous Term ψ Discontinuous gain γ
H-Infinity Coefficient μ Fault upper bound α
Disturbance upper bound ξ o Scaling factor β
Reachability Time T R H-Infinity Norm H
Minimum Eigen Value λ m i n Sliding surface σ ( . )
Space Vector Pulse width ModulationSVPWM

Appendix A

Proof for Proposed Form of Matrix ‘P’

For the form of the Lyapunov matrix P as used in [23,26], for the complete Lyapunov function in Equation (73), substituting complete forms from Equations (74) and (77) in (73)
V e ˙ = e s T e o T f T ξ T [ M ] e s e o f ξ
M 1 = A 11 ¯ T P 1 + P 1 A 11 ¯ P 1 A 12 ¯ + A 21 ¯ T T P o P 1 G 1 P 1 L G 2 P 1 L E o P 1 L D o A 12 ¯ T P 1 + P o T A 21 ¯ A 22 ¯ T T P o + P o T A 22 ¯ G 2 T T T P o P o T G 2 P o T E o P o T D o E o T L T P 1 E o T T T P o 0 0 D o T L T P 1 D o T T T P o 0 μ I ,
and Considering P = P 11 P 12 P 21 P 22 > 0 (and) A 11 ¯ A 12 ¯ A 21 ¯ A 22 ¯ in term ( A c G C c ) T P + P ( A G C ) + C c T C c and comparing with concerning part of matrix given in vector Lyapunov equation in (104) implies
P ¯ = P 1 0 0 T T P o T
Because the form of P ¯ is determined from the T L transformed system, applying the inverse transformation T L 1 P ¯ T L gives
P = P 1 P 1 L L P o T 0 P o + T T P 1 L
Applying the transformation T L T P T L on P ¯ gives:
P = P 1 P 1 L L T P 1 T T P o T + L T P L > 0

References

  1. Olivares, D.E.; Mehrizi-Sani, A.; Etemadi, A.H.; Cañizares, C.A.; Iravani, R.; Kazerani, M.; Hajimiragha, A.H.; Gomis-Bellmunt, O.; Saeedifard, M.; Palma-Behnke, R.; et al. Trends in microgrid control. IEEE Trans. Smart Grid 2014, 5, 1905–1919. [Google Scholar]
  2. Rocabert, J.; Luna, A.; Blaabjerg, F.; Rodriguez, P. Control of power converters in AC microgrids. IEEE Trans. Power Electron. 2012, 27, 4734–4749. [Google Scholar]
  3. Vasquez, J.C.; Guerrero, J.M.; Miret, J.; Castilla, M.; De Vicuna, L.G. Hierarchical control of intelligent microgrids. IEEE Ind. Electron. Mag. 2010, 4, 23–29. [Google Scholar]
  4. Poon, J.; Konstantakopoulos, I.C.; Spanos, C.; Sanders, S.R. Real-time model-based fault diagnosis for switching power converters. In Proceedings of the 2015 IEEE Applied Power Electronics Conference and Exposition (APEC), Charlotte, NC, USA, 15–19 March 2015; pp. 358–364. [Google Scholar]
  5. Odgaard, P.F.; Aubrun, C.; Majanne, Y. Fault tolerant control of power grids. Int. J. Robust Nonlinear Control 2014, 24, 1281–1282. [Google Scholar]
  6. Rasheduzzaman, M.; Mueller, J.A.; Kimball, J.W. An accurate small-signal model of inverter-dominated islanded microgrids using dq reference frame. IEEE J. Emerg. Sel. Top. Power Electron. 2014, 2, 1070–1080. [Google Scholar]
  7. Hassan, M.; Abido, M. Optimal design of microgrids in autonomous and grid-connected modes using particle swarm optimization. IEEE Trans. Power Electron. 2010, 26, 755–769. [Google Scholar]
  8. Li, J.; Zhang, X.; Li, W. An efficient wind-photovoltaic hybrid generation system for DC micro-grid. In Proceedings of the 8th International Conference on Advances in Power System Control, Operation and Management (APSCOM 2009), London, UK, 8–11 November 2009. [Google Scholar]
  9. Kroutikova, N.; Hernandez-Aramburo, C.A.; Green, T.C. State-space model of grid-connected inverters under current control mode. IET Electr. Power Appl. 2007, 1, 329–338. [Google Scholar]
  10. Andrade, F.; Romeral, L.; Kampouropoulos, K.; Cusido, J. New mathematical model of an inverter-based generator for stability studies of microgrid systems. In Proceedings of the 4th International Conference on Power Engineering, Energy and Electrical Drives, Istanbul, Turkey, 13–17 May 2013; pp. 944–949. [Google Scholar]
  11. Park, J.D.; Candelaria, J. Fault detection and isolation in low-voltage DC-bus microgrid system. IEEE Trans. Power Deliv. 2013, 28, 779–787. [Google Scholar]
  12. Hu, K.; Liaw, C. On the flywheel/battery hybrid energy storage system for DC microgrid. In Proceedings of the 2013 1st International Future Energy Electronics Conference (IFEEC), Tainan, Taiwan, 3–6 November 2013; pp. 119–125. [Google Scholar]
  13. Johnson, B.; Lasseter, R. An industrial power distribution system featuring UPS properties. In Proceedings of the IEEE Power Electronics Specialist Conference—PESC’93, Seattle, WA, USA, 20–24 June 1993; pp. 759–765. [Google Scholar]
  14. Alwi, H.; Edwards, C.; Marcos, A. Fault reconstruction using a LPV sliding mode observer for a class of LPV systems. J. Frankl. Inst. 2012, 349, 510–530. [Google Scholar]
  15. Kalsi, K.; Hui, S.; Żak, S.H. Unknown input and sensor fault estimation using sliding-mode observers. In Proceedings of the 2011 American Control Conference, San Francisco, CA, USA, 29 June–1 July 2011; pp. 1364–1369. [Google Scholar]
  16. Walcott, B.; Zak, S. State observation of nonlinear uncertain dynamical systems. IEEE Trans. Autom. Control 1987, 32, 166–170. [Google Scholar]
  17. Garcia, E.A.; Frank, P.M. Deterministic nonlinear observer-based approaches to fault diagnosis: A survey. Control Eng. Pract. 1997, 5, 663–670. [Google Scholar]
  18. Chen, J.; Patton, R.J.; Chen, Z. Active fault-tolerant flight control systems design using the linear matrix inequality method. Trans. Inst. Meas. Control 1999, 21, 77–84. [Google Scholar]
  19. Arsan, M.; Mouyon, P.; Magni, J.F. Fault diagnosis in the presence of parametric variations. IFAC Proc. Vol. 1994, 27, 83–88. [Google Scholar]
  20. Zhang, Y.; Jiang, J. Active fault-tolerant control system against partial actuator failures. IEE Proc. Control Theory Appl. 2002, 149, 95–104. [Google Scholar]
  21. Sreedhar, R.; Fernandez, B.; Masada, G. Robust fault detection in nonlinear systems using sliding mode observers. In Proceedings of the IEEE International Conference on Control and Applications, Vancouver, BC, Canada, 13–16 September 1993; pp. 715–721. [Google Scholar]
  22. Hermans, F.; Zarrop, M. Sliding mode observers for robust sensor monitoring. IFAC Proc. Vol. 1996, 29, 6530–6535. [Google Scholar]
  23. Utkin, V.I. Sliding Modes in Optimization and Control Problems; Springer: Berlin/Heidelberg, Germany, 1992. [Google Scholar]
  24. Edwards, C.; Spurgeon, S.K. On the development of discontinuous observers. Int. J. Control 1994, 59, 1211–1229. [Google Scholar]
  25. Edwards, C.; Spurgeon, S.K.; Patton, R.J. Sliding mode observers for fault detection and isolation. Automatica 2000, 36, 541–553. [Google Scholar]
  26. Tan, C.P.; Edwards, C. An LMI approach for designing sliding mode observers. Int. J. Control 2001, 74, 1559–1568. [Google Scholar]
  27. Tan, C.P.; Edwards, C. Sliding mode observers for detection and reconstruction of sensor faults. Automatica 2002, 38, 1815–1821. [Google Scholar]
  28. Tan, C.P.; Edwards, C. Sliding mode observers for robust detection and reconstruction of actuator and sensor faults. Int. J. Robust Nonlinear Control IFAC-Affil. J. 2003, 13, 443–463. [Google Scholar]
  29. Shin, D.; Moon, G.; Kim, Y. Design of reconfigurable flight control system using adaptive sliding mode control: Actuator fault. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2005, 219, 321–328. [Google Scholar]
  30. Edwards, C.; Tan, C.P. Sensor fault tolerant control using sliding mode observers. Control Eng. Pract. 2006, 14, 897–908. [Google Scholar]
  31. Yan, X.G.; Edwards, C. Nonlinear robust fault reconstruction and estimation using a sliding mode observer. Automatica 2007, 43, 1605–1614. [Google Scholar]
  32. Aldeen, M.; Sharma, R. Estimation of states, faults and unknown disturbances in non-linear systems. Int. J. Control 2008, 81, 1195–1201. [Google Scholar]
  33. Dhahri, S.; Sellami, A. Robust sliding mode observer design for a class of uncertain linear systems with fault reconstruction synthesis. Int. J. Phys. Sci. 2012, 7, 1259–1269. [Google Scholar]
  34. Shi, Y.T.; Kou, Q.; Sun, D.H.; Li, Z.X.; Qiao, S.J.; Hou, Y.J. H fault tolerant control of WECS based on the PWA model. Energies 2014, 7, 1750–1769. [Google Scholar]
  35. Cecati, C.; Genduso, F.; Miceli, R.; Galluzzo, G.R. A suitable control technique for fault-tolerant converters in distributed generation. In Proceedings of the 2011 IEEE International Symposium on Industrial Electronics, Gdansk, Poland, 27–30 June 2011; pp. 107–112. [Google Scholar]
  36. Gholami, S.; Saha, S.; Aldeen, M. Sensor fault tolerant control of microgrid. In Proceedings of the 2016 IEEE Power and Energy Society General Meeting (PESGM), Boston, MA, USA, 17–21 July 2016; pp. 1–5. [Google Scholar]
  37. Van, M. An enhanced robust fault tolerant control based on an adaptive fuzzy PID-nonsingular fast terminal sliding mode control for uncertain nonlinear systems. IEEE/ASME Trans. Mechatron. 2018, 23, 1362–1371. [Google Scholar]
  38. Azizi, A.; Nourisola, H.; Shoja-Majidabad, S. Fault tolerant control of wind turbines with an adaptive output feedback sliding mode controller. Renew. energy 2019, 135, 55–65. [Google Scholar]
  39. Zhang, J.; Li, S.; Xiang, Z. Adaptive fuzzy finite-time fault-tolerant control for switched nonlinear large-scale systems with actuator and sensor faults. J. Frankl. Inst. 2020, 357, 11629–11644. [Google Scholar]
  40. Lu, K.; Zhou, W.; Zeng, G.; Zheng, Y. Constrained population extremal optimization-based robust load frequency control of multi-area interconnected power system. Int. J. Electr. Power Energy Syst. 2019, 105, 249–271. [Google Scholar]
  41. Lu, K.D.; Zeng, G.Q.; Luo, X.; Weng, J.; Zhang, Y.; Li, M. An adaptive resilient load frequency controller for smart grids with DoS attacks. IEEE Trans. Veh. Technol. 2020, 69, 4689–4699. [Google Scholar]
  42. Rotondo, D.; Buciakowski, M. Guaranteed cost estimation and control for a class of nonlinear systems subject to actuator saturation. Eur. J. Control 2021, 61, 119–132. [Google Scholar]
  43. Prochazka, M. Modeling of Current Transformers Under Saturation Conditions. Adv. Electr. Electron. Eng. 2011, 5, 94–97. [Google Scholar]
  44. Matussek, R.; Dzienis, C.; Blumschein, J.; Schulte, H. Current transformer model with hysteresis for improving the protection response in electrical transmission systems. J. Phys. Conf. Ser. 2014, 570, 062001. [Google Scholar]
  45. Muthumuni, D.; Ruchkall, L.; Jayasinghe, R. Modeling current transformer (CT) saturation for detailed protection studies. Pulse Manit. Hvdc Res. Cent. J. 2011, 1–5. [Google Scholar]
  46. Abid, M. Fault Detection in Nonlinear Systems: An Observer-Based Approach; Universitat Duisburg-Essen: Duisburg, Germany, 2010. [Google Scholar]
  47. Chen, W. Fault Detection and Isolation in Nonlinear Systems: Observer and Energy-Balance Based Approaches; Südwestdeutscher Verlag für Hochschulschriften: Duisberg, Germany , 2012. [Google Scholar]
  48. Qayyum Khan, A. Observer-Based Fault Detection in Nonlinear Systems. Ph.D. Thesis, Institute of Automatic Control and Complex Systems (AKS), Duisburg, Germany, 2010. [Google Scholar]
  49. Raza, M.T. Fault Diagnosis and Fault Tolerant Control of Switched Dynamical Systems. Ph.D. Thesis, Pakistan Institute of Engineering and Applied Sciences, Nilore, Islamabad, Pakistan, 2016. [Google Scholar]
  50. Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. Sliding Mode Control and Observation; Springer: Berlin/Heidelberg, Germany, 2014; Volume 10. [Google Scholar]
  51. Starkov, K.E.; Coria, L.N.; Aguilar, L.T. On synchronization of chaotic systems based on the Thau observer design. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 17–25. [Google Scholar]
Figure 1. Sensor (C.T/P.T) fault diagnosis and fault-tolerant control approach.
Figure 1. Sensor (C.T/P.T) fault diagnosis and fault-tolerant control approach.
Sensors 22 02524 g001
Figure 2. SLD of microgrid.
Figure 2. SLD of microgrid.
Sensors 22 02524 g002
Figure 3. PI and droop control for voltage-frequency regulation.
Figure 3. PI and droop control for voltage-frequency regulation.
Sensors 22 02524 g003
Figure 4. Fault diagnostic observer based input current estimation (d-component).
Figure 4. Fault diagnostic observer based input current estimation (d-component).
Sensors 22 02524 g004
Figure 5. Fault diagnostic observer based output current estimation (d-component).
Figure 5. Fault diagnostic observer based output current estimation (d-component).
Sensors 22 02524 g005
Figure 6. Fault diagnostic observer based state output voltage estimation (q-component).
Figure 6. Fault diagnostic observer based state output voltage estimation (q-component).
Sensors 22 02524 g006
Figure 7. Fault diagnostic observer based state estimation error (dq-components).
Figure 7. Fault diagnostic observer based state estimation error (dq-components).
Sensors 22 02524 g007
Figure 8. Voltage (dq) fault reconstruction.
Figure 8. Voltage (dq) fault reconstruction.
Sensors 22 02524 g008
Figure 9. Current (dq) fault reconstruction.
Figure 9. Current (dq) fault reconstruction.
Sensors 22 02524 g009
Figure 10. Voltage (dq) fault estimation error.
Figure 10. Voltage (dq) fault estimation error.
Sensors 22 02524 g010
Figure 11. Current (dq) fault estimation error.
Figure 11. Current (dq) fault estimation error.
Sensors 22 02524 g011
Figure 12. Current (dq) disturbance estimation.
Figure 12. Current (dq) disturbance estimation.
Sensors 22 02524 g012
Figure 13. FTC performance for q-component of output current.
Figure 13. FTC performance for q-component of output current.
Sensors 22 02524 g013
Figure 14. FTC performance for d-component of output voltage.
Figure 14. FTC performance for d-component of output voltage.
Sensors 22 02524 g014
Figure 15. Current (dq) and voltage (dq) fault reconstruction (ramp fault/Case-II).
Figure 15. Current (dq) and voltage (dq) fault reconstruction (ramp fault/Case-II).
Sensors 22 02524 g015
Table 1. Microgrid system parameters.
Table 1. Microgrid system parameters.
ParameterValueParameterValue
V d c 600 V V g 600 V
L f 1 4.20 mH L f 2 4.20 mH
L c 1 0.50 mH L c 2 0.50 mH
C f 1 15 μ F C f 2 15 μ F
R d 1 2.025 Ω R d 2 2.025 Ω
r f 1 0.50 Ω r f 2 0.50 Ω
r f 1 0.09 Ω r f 2 0.09 Ω
ω c 50.26 rad/s ω n 377 rad/s
V o q n 85 V m , n 1 / 1000
θ g r i d 60 L f 3 4.20 mH
L c 3 0.50 mH C f 3 15 μ F
R d 3 2.025 Ω r f 3 0.50 Ω
r f 3 0.09 Ω ω P L L 377 rad/s
ine ω c , P L L 7853.98 rad/s
Table 2. Controller gains/parameters.
Table 2. Controller gains/parameters.
PI GainsParameterValue
Voltage k p v d , k p v q 0.5
Controllers k i v d , k i v q 25
Current k p c d , k p c q 1
Controllers k i c d , k i c q 100
PLL k p P L L 0.25
Controller k i P L L 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shahzad, E.; Khan, A.U.; Iqbal, M.; Saeed, A.; Hafeez, G.; Waseem, A.; Albogamy, F.R.; Ullah, Z. Sensor Fault-Tolerant Control of Microgrid Using Robust Sliding-Mode Observer. Sensors 2022, 22, 2524. https://doi.org/10.3390/s22072524

AMA Style

Shahzad E, Khan AU, Iqbal M, Saeed A, Hafeez G, Waseem A, Albogamy FR, Ullah Z. Sensor Fault-Tolerant Control of Microgrid Using Robust Sliding-Mode Observer. Sensors. 2022; 22(7):2524. https://doi.org/10.3390/s22072524

Chicago/Turabian Style

Shahzad, Ebrahim, Adnan Umar Khan, Muhammad Iqbal, Ahmad Saeed, Ghulam Hafeez, Athar Waseem, Fahad R. Albogamy, and Zahid Ullah. 2022. "Sensor Fault-Tolerant Control of Microgrid Using Robust Sliding-Mode Observer" Sensors 22, no. 7: 2524. https://doi.org/10.3390/s22072524

APA Style

Shahzad, E., Khan, A. U., Iqbal, M., Saeed, A., Hafeez, G., Waseem, A., Albogamy, F. R., & Ullah, Z. (2022). Sensor Fault-Tolerant Control of Microgrid Using Robust Sliding-Mode Observer. Sensors, 22(7), 2524. https://doi.org/10.3390/s22072524

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