An Algorithm for Recognition of Fault Conditions in the Utility Grid with Renewable Energy Penetration

: Penetration level of renewable energy (RE) in the utility grid is continuously increasing to minimize the environmental concerns, risk of energy security, and depletion of fossil fuels. The uncertain nature and availability of RE power for a short duration have created problems related to the protection, grid security, power reliability, and power quality. Further, integration of RE sources near the load centers has also pronounced the protection issues, such as false tripping, delayed tripping, etc. Hence, this paper introduces a hybrid grid protection scheme (HGPS) for the protection of the grid with RE integration. This combines the merits of the Stockwell Transform, Hilbert Transform, and Alienation Coefﬁcient to improve performance of the protection scheme. The Stockwell Transform-based Median and Summation Index (SMSI) utilizing current signals, Hilbert Transform-based derivative index (HDI) utilizing voltage signals, and Alienation Coefﬁcient index (ACI) utilizing voltage signals were used to compute a proposed Stockwell Transform-, Hilbert Transform-, and Alienation-based fault index (SAHFI). This SAHFI was used to recognize the fault conditions. The fault conditions were categorized using the number of faulty phases and the proposed Stockwell Transform and Hilbert Transform-based ground fault index (SHGFI) utilizing zero sequence currents. The fault conditions, such as phase and ground (PGF), any two phases (TPF), any two phases and ground (TPGF), all three phases (ATPF), and all three phases and ground (ATPGF), were recognized effectively, using the proposed SAHFI. The proposed method has the following merits: performance is least affected by the noise, it is effective in recognizing fault conditions in minimum time, and it is also effective in recognizing the fault conditions in different scenarios of the grid. Performance of the proposed approach was found to be superior compared to the discrete wavelet transform (DWT)-based method reported in the literature. The study was performed using the hybrid grid test system realized by integrating wind and solar photovoltaic (PV) plants to the IEEE-13 nodes network in MATLAB software.


Introduction
Capability of the renewable energy (RE) sources to meet the growing demand of electrical power has led to the increased penetration level of these sources in the network of utilities in recent years [1]. This has presented technical challenges to the power system operators and planners in terms of frequency variations, voltage regulations, efficiency, power quality (PQ) issues, grid stability issues, power reliability, and protection issues [2]. The protection challenges are due to the conversion of passive traditional feeders into active ones by the use of RE sources near load centers. Protection challenges are also observed due to reduction in the number of radial feeders and increased use of multi-tapped transmission and distribution lines for integrating RE sources and simultaneously feeding the loads [3]. A detailed study of the protection schemes and protection challenges in the power system with high penetration level of RE sources is reported in Reference [4][5][6]. The central part of the protection scheme is the detection and classification of fault conditions. Signal processing, artificial intelligence, and machine learning techniques play important role in detection and classification of the fault conditions [7]. An algorithm for detection and classification of faults in the presence of RE generation using combined use of wavelet singular entropy (WSE) and fuzzy logic (FL) is reported in Reference [8]. This is effective in identifying the faults in a time duration of 10 ms after fault inception and differentiate the fault conditions from islanding events in the presence of RE sources. A fault recognition method using a fault index (FI) computed using the detailed coefficients obtained by the decomposition of current signals is proposed in Reference [9]. Application of this method for recognition of faults in the presence of wind generation by the use of FI computed using the detail coefficient at first decomposition level of current signals using DWT is reported in Reference [10]. Application of the approach for recognition of faults in the presence of solar energy is reported in Reference [11]. Performance of this DWT-based method is affected by noise, fault recognition time is high, and range of variation of peak values of FI is large, which may generate false tripping commands during the operational events. A method for identification of faults in a utility grid when a solar photovoltaic (PV) plant is integrated and based on the use of wavelet multi-resolution singular spectrum entropy (WMRSSE) and Support Vector Machine (SVM) is reported in Reference [12]. This method is effective in recognizing the different fault conditions at fast rate in the presence of solar PV generation. In Reference [13], authors proposed a method using the Machine Learning approach for tracing of the fault current for recognizing the fault conditions in distribution systems in the presence of distributed energy resources (DER). This overcomes the demerits of the artificial intelligent (AI)and SVM-based methods. Chen et al. [14] introduced an effective scheme for protection of a micro-grid in grid-connected mode by recognizing the fault conditions using wavelet energy and fuzzy neural network. In Reference [15], authors introduced a fuzzy-supported intelligent technique to detect and classify the fault conditions on distribution power lines where RE sources are integrated. This method is effective in recognizing the fault conditions in a time duration between quarter cycle to one cycle. A multi-class adaptive neuro-fuzzy classifier (MC-NFC) is introduced in Reference [16] to recognize the fault conditions in the presence of solar PV generation. It is established that MC-NFC has better fault classification accuracy compared to the artificial neural network (ANN) classifier. A protection scheme using the alienation coefficient and Wigner distribution function is proposed for identification of faults based on the features of current and voltage signals in the presence of solar energy [17] and both solar and wind energy [18]. This algorithm is effective to recognize faults with high penetration level of RE even in the presence of high noise level in a duration less than quarter cycle period. This article presents a method for recognition of fault conditions in the presence of renewable energy generation in addition to the conventional generation and loads. This method overcomes the demerits of the above discussed protection approaches. Hence, the main contribution of this article can be summarized by the following points: • This article introduced a fault recognition approach that can be used for the protection of the grid with RE penetration. •  This protection scheme uses a fault index based on the features of voltage and currents extracted  using the Stockwell Transform, Hilbert Transform, and Alienation Coefficient.  • This approach has the merits that its performance is least affected by noise, recognizes fault conditions in minimum time, effectively recognizes the fault conditions in different scenarios of the grid, and is free from generation of the false tripping command. • Proposed approach effectively identify the faults using the number of faulty phases and the proposed Stockwell and Hilbert Transform-based ground fault index.
The article is structured in nine sections. The first section details the introductory part, including research gaps and contribution. The second section is related to the description of hybrid grid test system. In Section 3, protection algorithm is discussed and different proposed indexes are described. Results of simulation for recognition of fault conditions are detailed under Section 4. The fifth section describes the categorization of fault conditions, while the sixth section describes the results for various cases of study. The application of the proposed protection approach to differentiate switching transients is discussed under Section 7. A comparative study is incorporated under Section 8. Lastly, the conclusions are discussed under Section 9.

Proposed Test System
A hybrid grid test system (HGTS) is realized using the standard distribution network of IEEE-13 nodes and described in Figure 1. This HGTS has capacity of 5 MVA at frequency of 60 Hz and voltage levels of 0.48 kV and 4.16 kV. The node 634 is rated at voltage level of 0.48 kV and the rest of the nodes are rated at 4.16 kV. This test system is integrated to the utility grid rated at voltage level of 115 kV using a grid integration transformer (GIT) at node 650, which is considered as grid integration point (GIP). Loads with both the balanced and unbalanced nature are integrated to this grid [19,20]. A solar PV plant (SPP) of capacity 1 MW and output voltage level of 260 V is integrated on node 680 using a transformer, which is designated as SPPT. The SPP is designed using the single diode equivalent circuit of the solar cell using the modeling and data reported in Reference [21][22][23]. A wind power plant (WPP) of capacity 1.5 MW and output voltage level of 575 V is integrated on node 680 using a transformer, which is designated as WPPT. The WPP is designed using the modeling and data reported in Reference [24,25]. The node 680 is considered as the point of common coupling (PCC) for integration of RE plants. The current and voltage measurements are performed at node 650, where the proposed protection scheme is proposed to be installed. Hence, this node is also considered as hybrid grid protection relay location node (HGPRLN) in addition to the GIP. The fault conditions are realized on the node 675, which is considered as fault condition node (FCN). A distribution transformer (DT) is connected between the nodes 634 and 633 to maintain these nodes at different voltage levels. A description of the transformers used in the proposed HGTS is provided in Table 1, where the ratings, voltage levels, resistance, and reactance of the first winding (FW) and second winding (SW) of the transformers are provided. The load data, capacitor data, and feeder data reported in Reference [26,27] are used in the proposed study.

Protection Algorithm
The hybrid grid protection scheme (HGPS) for protection of the grid with RE integration introduced in this study is designed with five stages. The first stage calculates the Stockwell Transform based Median and Summation Index (SMSI) using the current signals. The second stage calculates the Hilbert Transform based derivative index (HDI) using the voltage signals. In third stage, Alienation Coefficient index (ACI) is computed using the voltage signals. In the fourth stage, the proposed Stockwell Transform-, Hilbert Transform-, and Alienation-based fault index (SAHFI) is computed using the SMSI, HDI, and ACI. In the fifth stage, different fault conditions are categorized using the number of faulty phases and the proposed Stockwell Transform-and Hilbert Transform-based ground fault index (SHGFI). All the stages of the proposed approach for fault recognition are detailed in Figure 2. A detailed description of the proposed indexes is provided in the following subsections.

Stockwell Transform-Based Median and Summation Index
The current signals are recorded on the bus 650, which is designated as the hybrid grid protection relay location node (HGPRLN). These current signals are processed using the Stockwell Transform (ST) to propose a Stockwell Transform based Median and Summation Index (SMSI) for recognition of fault conditions on the grid. The ST can be evaluated either by the use of the continuous WT (CWT) or the short time Fourier transform (STFT). The ST is effective for extraction of the information hidden in the phase spectrum and amplitude of the current signal. This information is obtained in the form of a complex matrix designated as STM. Each row of this matrix indicates definite frequency and each column indicates the definite time. High time resolution at a high frequency and a low time resolution at a low frequency is achieved using the ST. This has the additional advantage that there is minimum effect of noise on its performance and effectively identify the distortions of the waveform [28]. The current signal (i(t)) is decomposed using the short time Fourier transform (STFT), and output (STI) is computed as described below [29].
where t and f are, respectively, the spectral localization time and Fourier frequency, and g(t) represents the window function. The following Gaussian function (g(t)) is used in the above equation to evaluate the Stockwell Transform-based decomposition of i(t) [30].
Hence, the output matrix (which is complex in nature) of the Stockwell Transform-based decomposition of current signals (STM) can be given by the following relation [31].
This complex values STM matrix is used to extract the information of frequency and amplitude of the current signals. First, this complex values matrix is converted in a matrix with absolute values (AVM) using the following command.
Median of the matrix AVM is evaluated using the below detailed command.
The summation of each column of the AVM is evaluated using the below detailed command.
The MV and SV are multiplied element by element to compute the Stockwell Transform based median and summation index (SMSI) using the below detailed command.

Hilbert Transform-Based Derivative Index
Hilbert Transform (HT) is a signal processing approach, which can be implemented for the purpose of protection schemes. This is effective for spectral analysis and provide description of time-frequency-energy of the voltage time series data. The HT is effective for computation of the momentary frequencies, as well as amplitudes, which can be used for description of the voltage signals. Hence, HT can be effectively used for recognition of the fault conditions in the hybrid grid with RE penetration. The HT supported decomposition of voltage signals can be computed using the below mentioned relation [32].
Here, PV: Cauchy's principle value integral, t: time, and τ: time period. The voltage signal is used to compute Hilbert Transform-based derivative index (HDI). Voltage signal v(t) is decomposed using the Hilbert Transform to evaluate output matrix with absolute values (HTI) using the following command. HTI = abs(hilbert(v(t))).
The HTI is differentiated with respect to time for computing the HDI as detailed below.
This method has the merit that it gives physical instantaneous frequencies for voltage signal and minimum time for identification of the fault condition.

Alienation Coefficient Index
The sample supported alienation coefficient of voltage signals is computed for recognition of the faults and assigned the name Alienation Coefficient index (ACI). The correlation coefficient (CRC) between two voltage variables v 1 and v 2 is computed using the following relation [18].
Here, symbol representation, Ns: sample numbers in a cycle (Ns = 32 with sampling frequency of 1.92 kHz), v 1 : samples of voltage at time t 0 , v 2 : samples of voltage at time −T + t 0 , and T: time period. The ACI is computed using CRC by the following relation.
The moving window approach is used to compute ACI for a period of the quarter cycle. Here, samples of the quarter cycle considered are compared with the samples of the corresponding quarter cycle of past cycle, and the window is moved in forward direction with one sample step every time. This process is repeated. ACI helps to reduce the time of fault recognition due to the feature of sharp change at fault incidence moment, resulting in the fast protection scheme.

Stockwell, Alienation, and Hilbert Fault Index
The SMSI, HDI, and ACI are used to compute the Stockwell-, Hilbert-, and Alienation-based fault index (SAHFI), which is used to identify the different fault conditions. The SAHFI is evaluated using the multiplication of the SMSI, ACI, and HDI as detailed below.
Here, W: weight factor, which is taken as 100. This WF may be set depending on the complexity of the network and penetration level of the renewable energy at the time of installation of protection scheme. A threshold value of 2500 is selected for the SAHFI and designated as hybrid fault index threshold (HFITH). This algorithm has the merits that its performance is least affected by noise, recognizes fault conditions in minimum time, reduces the time of fault recognition, effectively recognizes the fault conditions in different scenarios of the grid, and is free from generation of false tripping command.

Stockwell Transform-and Hilbert Transform-Based Ground Fault Index
A Stockwell Transform and Hilbert Transform based ground fault index (SHGFI) is introduced to categorize the phase and ground (PGF) fault and any two phases (TPF) fault. This is also used to categorize all three phases and ground (ATPGF) fault and all three phases (ATPF) fault. Following procedure is used to find the SHGFI:-

•
Calculate the zero sequence current (I 0 ) using the currents I 1 , I 2 , and I 3 associated with phases A, B, and C, respectively, using the following relation.
• Calculate the SMSI 0 using the mathematical formulation detailed in Section 3.1 by processing the zero sequence current. • Calculate the HDI 0 using the mathematical formulation detailed in Section 3.2 by processing the zero sequence current. • Calculate the SHGFI by multiplying the SMSI 0 and HDI 0 as detailed below:- Here, W 0 is weight factor for the SHGFI, which is considered as 10 3 in this study. The threshold for the SHGFI (SHGFITH) is taken as 10 to detect the involvement of ground during the faulty condition.

Recognition of Fault Conditions: Results of Simulation
Results for implementation of the algorithm to identify and classify the fault conditions for protection of the grid network in the presence of RE are discussed in this section. The fault conditions, which include a phase and ground (PGF), any two phases (TPF), any two phases and ground (TPGF), all three phases (ATPF), and all three phases and ground (ATPGF), are analyzed. Results to classify the fault conditions are also discussed. The test bus 675 is considered as fault condition node (FCN) and voltage and currents are measured on test bus 650, which is considered as hybrid grid protection relay location node (HGPRLN).

Fault Condition Involving Ground and Phase A
A fault condition involving ground and phase A (PGF) is executed at 0.1 s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using the ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for PGF condition is computed using the SMSI, ACI, and HDI. The current signals, voltage signals, SMSI, ACI, HDI, and SAHFI for a PGF on phase A are described in Figure 3 in respective order. Figure 3c details that, due to the incidence of PGF on phase A, increase in the SMSI for phase A is high compared to the SMSI of phases B and C. This index has low values during the healthy condition. Figure 3d details that, due to the incidence of PGF on phase A, the ACI increases corresponding to all phases. This index also has low values during the healthy condition and increases sharply due to incidence of fault condition. Figure 3e details that, due to the incidence of PGF on phase A, increase in the HDI for phase A is high compared to the HDI of phases B and C. This index also has low values during the healthy condition. Figure 3f details that, due to the incidence of PGF on phase A, the proposed SAHFI has peak magnitude above the HFITH (2500) for phase A, whereas this index for phases B and C is below the HFITH. Hence, the condition of the PGF on phase A is effectively recognized using the proposed algorithm. This algorithm is also effective to recognize the PGF fault condition on phases B and C. A high resolution version of the SAHFI at the moment of incidence of PGF condition (phase A) is shown in Figure 4. The SAHFI associated with the faulty phase A achieves the value equal to HFITH (2500) at 0.10015 s after a time interval of 0.15 ms after fault incidence. Hence, the PGF event is detected in a time duration of (1/110)th fraction of a time of one cycle.

Fault Condition Involving Two Phases
A two phases fault (TPF) condition involving phases A and B is executed at 0.1 s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using the ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for TPF condition is computed using the SMSI, ACI, and HDI. The current signals, voltage signals, SMSI, ACI, HDI, and SAHFI for a TPF condition involving phases A and B are described in Figure 5 in respective order. Figure 5c details that, due to the incidence of the TPF on phases A and B, increase in the SMSI for phases A and B is high, and this SMSI is low for phase C. This index has low values during the healthy condition. Figure 5d details that, due to the incidence of TPF condition on phases A and B, the ACI increases corresponding to all phases. This index also has low values during the healthy condition and increase sharply due to incidence of TPF fault condition. Figure 5e details that, due to the incidence of TPF on phases A and B, increase in the HDI for phases A and B is high compared to the HDI of phase C. This index also has low values during the healthy condition. Figure 5f details that, due to the incidence of TPF on phases A and B, the proposed SAHFI has peak magnitude above the HFITH (2500) for phases A and B, whereas this index for phase C is below the HFITH. Hence, the condition of the TPF on phases A and B is effectively recognized using the proposed algorithm. This algorithm is also effective to recognize the TPF fault condition involving other combinations of phases, such as B and C and C and A. A high resolution version of the SAHFI at the moment of incidence of TPF condition on phases A and B is shown in Figure 6. The SAHFIs associated with the faulty phases A and B achieve the value equal to HFITH (2500) at 0.10002 s in a time interval of 0.02 ms after fault incidence. Hence, the PGF event is detected in a time duration of (1/830)th fraction of a time period of one cycle.

Fault Condition Involving Two Phases and Ground
A two phases fault condition involving phases A and B and ground (TPGF) is executed at 0.1 s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using the ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for TPF condition is computed using the SMSI, ACI, and HDI. The current signals, voltage signals, SMSI, ACI, HDI, and SAHFI for a TPGF condition involving phases A and B are described in Figure 7 in respective order. Figure 7c details that, due to the incidence of TPGF on phases A and B, increase in the SMSI for phases A and B is high and this SMSI is low for phase C. This index has low values during the healthy condition. Figure 7d details that, due to the incidence of TPGF condition on phases A and B, the ACI increases corresponding to all phases. This index also has low values during the healthy condition and increase sharply due to incidence of TPGF fault condition. Figure 7e details that, due to the incidence of TPGF on phases A and B, increase in the HDI for phases A and B is high compared to the HDI of phase C. This index also has low values during the healthy condition. Figure 7f details that, due to the incidence of TPGF on phases A and B, the proposed SAHFI has peak magnitude above the HFITH (2500) for phases A and B, whereas this index for phase C is below the HFITH. Hence, the condition of the TPGF on phases A and B is effectively recognized using the proposed algorithm. This algorithm is also effective to recognize the TPGF fault condition involving other combinations of phases, such as B and C and C and A. A high resolution version of the SAHFI at the moment of incidence of TPGF condition on phases A and B is shown in Figure 8. The SAHFI associated with the faulty phase A achieves the value equal to HFITH (2500) at 0.10005 s after a time interval of 0.5 ms after fault incidence. This SAHFI associated with the faulty phase B achieves the value equal to HFITH (2500) at 0.10001 s after a time interval of 0.1 ms after fault incidence. Hence, the TPGF event is detected in a time duration of (1/330)th fraction of a time duration of one cycle.

Fault Involving All the Phases
An all the three phases fault (ATPF) condition involving all phases is executed at 0.1 s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using the ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for ATPF condition is computed using the SMSI, ACI, and HDI. The current signals, voltage signals, SMSI, ACI, HDI, and SAHFI for a TPF condition involving all the phases are described in Figure 9 in respective order. Figure 9c details that, due to the incidence of ATPF on all phases, increase in the SMSI for all phases is high. This index has low values during the healthy condition. Figure 9d details that, due to the incidence of ATPF condition on all phases, the ACI increases corresponding to all phases. This index also has low values during the healthy condition and increase sharply due to incidence of ATPF fault condition. Figure 9e details that, due to the incidence of ATPF on all phases, increase in the HDI for all phases is high. This index also has low values during the healthy condition. Figure 9f details that, due to the incidence of ATPF on all phases, the proposed ATPF has peak magnitude above the HFITH (2500) for all the phases. Hence, the condition of the ATPF on all the phases is effectively recognized using the proposed algorithm. A high resolution version of the SAHFI at the moment of incidence of ATPF condition on all the phases is shown in Figure 10. The SAHFI associated with phase A achieves the value equal to HFITH (2500) at 0.10008 s after a time interval of 0.08 ms after fault incidence. This SAHFI associated with the faulty phase B achieves the value equal to HFITH (2500) at 0.10001 s after a time interval of 0.001 ms after fault incidence. Similarly, this SAHFI associated with the faulty phase C achieves the value equal to HFITH (2500) at 0.10015 s after a time interval of 0.15 ms after fault incidence. Hence, the ATPF event is detected in a time duration of (1/110)th fraction of a time duration of one cycle.

Fault Involving All the Phases and Ground
An all the three phases to ground fault (ATPGF) condition involving all phases is executed at 0.1s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using the ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for ATPGF condition is computed using the SMSI, ACI, and HDI. The current signals, voltage signals, SMSI, ACI, HDI, and SAHFI for an ATPGF condition are described in Figure 11 in respective order. Figure 11c details that, due to the incidence of ATPGF on all the phases, increase in the SMSI for all phases is high. This index has low values during the healthy condition. Figure 11d details that, due to the incidence of ATPGF condition on all the phases, the ACI increases corresponding to all the phases. This index also has low values during the healthy condition and increase sharply due to incidence of ATPGF fault condition. Figure 11e details that, due to the incidence of ATPGF on all the phases, increase in the HDI for all phases is high. This index also has low values during the healthy condition. Figure 11f details that, due to the incidence of ATPGF on all the phases, the proposed ATPGF has peak magnitude above the HFITH (2500) for all the phases. Hence, the condition of the ATPGF on all the phases is effectively recognized using the proposed algorithm.
A high resolution version of the SAHFI at the moment of incidence of ATPGF condition on all the phases is shown in Figure 12. The SAHFI associated with phase A achieves the value equal to HFITH (2500) at 0.10008 s after a time interval of 0.08 ms after fault incidence. This SAHFI associated with the faulty phase B achieves the value equal to HFITH (2500) at 0.10001 s after a time interval of 0.001 ms after fault incidence. Similarly, this SAHFI associated with the faulty phase C achieves the value equal to HFITH (2500) at 0.10015 s after a time interval of 0.15 ms after fault incidence. Hence, an ATPGF event is detected in a time duration of (1/110)th fraction of a time interval of one cycle.

Categorization of Fault Conditions
The fault conditions are categorized using the number of phases involved in the fault. The phases involved in the fault condition are recognized using the proposed SAHFI for a phase and comparing with the HFITH (2500). A faulty phase has a peak value of SAHFI higher compared to the HFITH, whereas this index has a peak value below HFITH for healthy phases. If only one phase is faulty then PGF fault is present. Two phases will be faulty during the conditions of TPF or the TPGF faults. All the three phases will be faulty when either the ATPF or the ATPGF fault is present.
There are two faulty phases during the events of TPF and TPGF. A Hilbert Transform and Stockwell Transform-based ground fault index (SHGFI) is introduced to categorize the TPF and TPGF faults, which is detailed in Figure 13. This is observed that for the event of TPF, the SHGFI has low values near to zero, whereas this index has a high magnitude peak for the TPGF faulty event at the time of fault incidence, which crosses the threshold value of SHGFI (SHGFITH). Here, SHGFITH is considered equal to 10. Therefore, it is concluded that the proposed method is effective in the categorization of all the types of investigated fault conditions. Similarly, the ATPF and ATPGF fault events are also successfully identified using the SHGFI.

Case Studies: Implementation of Protection Algorithm
The protection approach is validated for different fault conditions, such as different values of impedance during fault condition, fault incident angle (FIA), different FCN, and noisy condition. Results for the PGF condition are elaborated in this section. However, the protection approach is equally effective for the conditions of TPF, TPGF, ATPF, and ATPGF.

Impact of Fault Impedance
The algorithm is tested for different values of the impedance involved during the fault condition. The algorithm is tested for fault impedance values of 0 Ω, 2 Ω, 4 Ω, 6 Ω, 8 Ω, and 10 Ω during the PGF fault condition and peak magnitudes of respective SAHFI are tabulated in Table 2. From this it is inferred that SAHFI has a peak magnitude higher than the HFITH (2500) for phase A and lower for phases B and C for all the values of investigated fault impedances. Therefore, the PGF fault event is effectively recognized for the fault impedance values up to 10 Ω. The algorithm is also effective in recognizing the TPF, TPGF, ATPF, and ATPGF fault conditions for the fault impedance values up to 10 Ω.

Impact of FCN
Location of the fault at different buses of the test grid affects the SAHFI. Hence, the algorithm is tested by realizing the fault at different nodes. The peak magnitude of SAHFI during the condition of PGF at various nodes of the test grid in the presence of RE generation are tabulated in Table 3. From this it is inferred that SAHFI has peak magnitude higher than the HFITH (2500) for phase A and lower for phases B and C for all the locations of the faults. Therefore, the PGF fault event is effectively recognized when there is a fault incident at various nodes of the test grid. The algorithm is also effective in recognizing the TPF, TPGF, ATPF, and ATPGF fault conditions at different nodes.

Impact of Fault Incidence Angle
A fault condition of PGF on phase A is executed at 0.01 s with fault incident angle (FIA) of 90°of the waveform on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. The current signals with FIA of 90°are processed using ST to compute the SMSI. The voltage signals with FIA of 90°are processed using the Hilbert Transform to compute the HDI. The voltage signals with FIA of 90°are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for the PGF condition is computed using the SMSI, ACI, and HDI, which is illustrated in Figure 14. Figure 14 details that, due to the incidence of PGF condition on phase A, the proposed SAHFI has a peak magnitude above the HFITH (2500) for phase A, whereas this index for phases B and C is below the HFITH. Hence, the condition of the PGF on phase A with FIA of 90°is effectively recognized using the proposed algorithm. It is also observed that the voltage magnitude of the SAHFI has increased significantly compared to the noise-free environment. The peak magnitude of SAHFI during the condition of PGF with FIA of 0°, 45°, 90°, and 135°are tabulated in Table 4. From this it is inferred that SAHFI has peak magnitude higher than the HFITH (2500) for phase A and lower for phases B and C for all the FIAs. Therefore, the PGF fault event is effectively recognized for the different angles of incidence of faults.

Impact of Noise
A fault condition of PGF on phase A is executed at 0.1 s on FCN of the grid with RE sources. Signals of voltage and current are recorded at the HGPRLN. A noise of level 20 dB SNR is superimposed on both the voltage and current signals. The current signals with noise of 20 dB SNR are processed using ST to compute the SMSI. The voltage signals with noise of 20 dB SNR are processed using the Hilbert Transform to compute the HDI. The voltage signals with noise of 20 dB Signal to noise ratio (SNR) are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for PGF condition is computed using the SMSI, ACI, and HDI, which is illustrated in Figure 15. Figure 15 details that, due to the incidence of the condition on phase A, the proposed SAHFI has a peak magnitude above the HFITH (2500) for phase A, whereas this index for phases B and C is below the HFITH. Hence, the condition of the PGF on phase A with the superimposed noise of 20 dB SNR on both the voltage and current signals is effectively recognized using the proposed algorithm. It is also observed that the voltage magnitude of the SAHFI has increased significantly compared to the noise-free environment.

Performance of Protection Algorithm in the Presence of Switching Transients
It is expected that a protection scheme should not operate during the operational events, such as operation of a power plant, load, feeder, and capacitor banks. Hence, the section presents results for the operational conditions to establish that trip commands for the circuit breaker (CB) will not be generated during the operational events.

Tripping of Wind Plant
The wind power plant of capacity 1.5 MW integrated on bus 680 is tripped at 0.01 s. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for the event of wind plant tripping is computed using the SMSI, ACI, and HDI. The SMSI, ACI, HDI and SAHFI and described in Figure 16a-d in respective order. Figure 16a details that the magnitude of SMSI increase for all phases at instant of tripping of the wind plant. Figure 16b details that the magnitude of ACI increase for all phases at instant of tripping of the wind plant. Figure 16c details that the magnitude of HDI increase for all phases at instant of tripping of wind plant. Figure 16d details that the magnitude of the increase for all phases at instant of tripping of the wind plant, but the increase in the magnitude is much less compared to the SAHFITH (2500). Hence, the peak magnitude of the corresponding to all phases is below the SAHFITH. Therefore, a CB trip command will not be generated during this event of wind plant tripping, and the power system network will operate in normal operating condition.

Capacitor Bank and RL Load Operation
The load consisting of 843 kW and 462 kVAr and capacitor bank of 600 kVAr connected on bus 675 of the test system are simultaneously tripped at 0.1 s. Signals of voltage and current are recorded at the HGPRLN. The current signals are processed using ST to compute the SMSI. The voltage signals are processed using the Hilbert Transform to compute the HDI. The voltage signals are also processed using the Alienation Coefficient to compute the ACI. The SAHFI for the PGF condition is computed using the SMSI, ACI, and HDI. The SMSI, ACI, HDI, and SAHFI for event of simultaneous tripping of capacitor bank and load are described in Figure 17a-d in respective order. Figure 17a details that magnitude of the SMSI increase for all phases at instant of simultaneous tripping of capacitor bank and load. Figure 17b details that the magnitude of the ACI increases for all phases at instant of simultaneous tripping of capacitor bank and load. Figure 17c details that the magnitude of the HDI increases for all phases at instant of simultaneous tripping of capacitor bank and load. Figure 17d details that the magnitude of the SAHFI increases for all phases at instant of simultaneous tripping of capacitor bank and load, but the the increase in magnitude is much less compared to the SAHFITH (2500). Hence, the peak magnitude of the SAHFI corresponding to all phases is below the SAHFITH. Therefore, a CB trip command will not be generated during this event of simultaneous tripping of the capacitor bank and load and the power system network will operate in a normal operating condition.

Performance of Protection Scheme: A Comparative Study
The performance of the protection approach discussed in this article is compared with the DWT based approach reported in Reference [9][10][11]. The DWT based approach used a fault index (FI) computed using the detailed coefficients obtained by the decomposition of current signals and applied for the protection of the transmission line [9], power system with wind generation [10] and power system network with solar PV generation [11]. Performance of this DWT-based method is affected by noise, the fault recognition time is high, and the range of variation of peak values of FI is large, which may generate false tripping commands during the operational events. The FI computed using detailed coefficient at third decomposition level and used for protection of transmission line is shown in Figure 18. FI computed using detailed coefficient at first decomposition level and used for protection of power system with wind generation is shown in Figure 19. It is observed from Figures 18 and 19 that the fault recognition time by the use of DWT approach is high compared to the method proposed in this article. Further, the performance of the proposed approach is not affected by the noise level of 20 dB SNR, while the DWT based approach is effective only for the noise level up to 50 dB SNR. Further, the variations in the SAHFI proposed in this article are small from one fault event to another and during variations of parameters, such as FIA, fault impedance, and FCN. There is a wide gap between the SAHFI for the faulty events and operational events. Hence, the proposed method is free from the generation of false tripping command compared to the DWT approach. Therefore, the protection approach proposed in this article is more effective compared to the DWT based approach reported in Reference [9][10][11].

Conclusions
A protection scheme for the grid with RE penetration is proposed in this article. This is based on the indexes proposed, such as SMSI, HDI, and ACI. These indexes are used to compute the fault index SAHFI, which is used for the recognition of fault conditions. This is achieved by comparing the absolute magnitude of SAHFI with its threshold value HFITH, which is considered equal to 2500. If computed value of SAHFI for a phase is greater compared to HFITH, then fault is present. Fault conditions are categorized using the number of faulty phases and a proposed ground fault index SHGFI. Involvement of ground in the fault condition is identified by comparing the computed value of SHGFI with the threshold SHGFITH, which is considered equal to 10. If computed value of SHGFI is greater compared to SHGFITH, then fault is present. It is concluded that the fault conditions, such as phase and ground (PGF), any two phases (TPF), any two phases and ground (TPGF), all three phases (ATPF), and all three phases and ground (ATPGF), are recognized effectively using the proposed SAHFI and SHGFI. Further, the proposed method has the follwing merits: performance is least affected by noise, is effective in recognizing fault conditions in minimum time, and is effective in recognizing the fault conditions in different scenarios of the grid. This scheme is free from the generation of the false tripping indications. This approach effectively recognizes the fault conditions in such cases as variations of fault impedance, FIA, presence of noise, and different FCN. This approach effectively discriminates the faulty condition from the healthy condition. The performance of the proposed approach was found to be superior compared to the DWT-based method reported in literature. Thus, this study was successfully performed on the hybrid grid formed by integrating wind and solar generators to the IEEE-13 node test system using MATLAB software.