A Novel Approach for Secure Hybrid Islanding Detection Considering the Dynamic Behavior of Power and Load in Electrical Distribution Networks

: In the arena of modern electrical power distribution systems, distributed generators (DGs) are emerging as a manifestation of electric power personalization. Even though DGs have various advantages, unintentional islanding phenomena caused by DGs during abnormal grid operations can damage equipment connected to the grid. Therefore, islanding detection mechanisms are essential for DGs in grid-connected mode to disconnect the DG from the grid in case of grid abnormalities by obeying to speciﬁc grid codes. In this regard, a novel approach to develop a secure hybrid islanding detection method (IDM) is presented in this paper. The proposed hybrid IDM is developed by combining two passive IDMs known as rate of change of active power and rate of change of reactive power with an active IDM called load connecting strategy. An 11 kV Malaysian distribution system integrated with three types of DGs, namely synchronous generator, photovoltaic, and biomass, has been chosen as a testbed for the veriﬁcation of the proposed hybrid IDM. Seven different case studies have been conducted in the PSCAD/EMTDC platform to validate the performance of the proposed IDM for islanding and non-islanding events. The simulation results conﬁrm that the proposed IDM can detect islanding within 0.09 s, which is within 2 s complying with IEEE and IEC standards. Further, a comparative study based on the detection time and non-detection zone has been carried out, which has conﬁrmed that the proposed IDM demonstrates better performance compared to the previously developed hybrid IDMs.


Introduction
All over the world in recent times, the use of distributed generators (DGs) has increased in distribution networks due to their advantages such as cost-effective solutions to increasing load demand, peak shaving, improved power quality and reduced power losses to enhance the distribution system's reliability and performance [1,2]. Though DG has become an essential part of modern power system networks, during grid abnormalities (undervoltage and short circuit events), a DG needs to be detached from the grid for protection purposes by satisfying a set of grid standards produced by IEEE and IEC [3][4][5]. If the DG does not disconnect from the grid during unprecedented grid abnormality, then this phenomenon may lead to the damage of equipment connected to the grid and becomes a threat to the lives of distribution system line workers [6][7][8]. Therefore, to ensure the protection of DGs against unintentional islanding, islanding should take place to detach the DG from the grid within a specified time (seconds) using an autonomous disconnecting technique known as the islanding detection method (IDM). As per IEEE standard 1547, once islanding occurs, within 100 cycles or 2 s, it is mandatory for the DG to be detached from Maintenance and implementation cost are very high. [14,15] Phasor Measuring Units Detection process does not require any extra device, so it is easy to implement Shows less robustness to handle different types of signals across the network. [16,17] Transfer Trip It is a very simple concept to implement with a very small NDZ.
Maintenance and implementation cost are very high. [18] Active IDMs

Active and Reactive Power Injection
The detection accuracy is high due to the injection of powers.
The voltage at distribution side rises which is a concern. [19] Active Frequency Drift Balanced islanding conditions and small NDZ can be achieved. Power quality degrades. [20] Impedance Measurement The method operates well because of the absence of NDZ This method is not suitable for parallel inverter connection. [21] Harmonic Signal Injection During islanding, power balance can be achieved among generator and demand.
Detection time is high. [22] Slip Mode Frequency Shift This method has small NDZ. Inaccurate in measurement due to the presence of phase shift parabola [23] Sandia Frequency Shift (SFS) Implementation is easy due to having very small NDZ.
Power system stability and quality are the concerns. [24] Sandia Voltage Shift Islanding detection speed is fast. Power quality and transient response of the system get affected. [25]

Method Advantage Disadvantage References
Frequency Jump Effective for non-parallel multi DGs. Less efficient for parallel DGs. [26] Virtual Capacitor and Inductor Harmonics are lower at the output. Power quality degrades. [27] Passive IDMs Rate of change of Power (ROCOP) Suitable for large power mismatch. Selection of threshold values are difficult. [28,29] Rate of change of Frequency (ROCOF) Islanding detection speed is fast.
For small DGs, threshold values can be chosen accurately but for medium and large DGs it is difficult. [30] Over/Under Voltage and Frequency Low cost and implementation is easy. Due to large NDZ, detection time is long. [31] Change of Impedance Suitable for small, medium and large DGs with large power mismatch.
Initialization of unwanted tripping is a concern. [32] Voltage Unbalance It can easily identify unbalance in the 3-phase system.
For a single-phase system, it is not suitable. [33] ROCOF over ROCOP Small power mismatch can be detected between load and DG Threshold selection can cause incorrect detection. [34] Phase Jump Implementation is easy. When the DG meets local demand, it fails to detect islanding condition. [35] Hybrid IDMs

Voltage and Reactive Power Shift
Fault tolerant capacity and robustness of the system improved.
Power system stability and quality are the concerns. [36] SFS and Q-f Based Scheme Voltage regulation, and power factor improved.
Selection of threshold values are difficult and power quality degrade. [37] Positive Feedback and voltage unbalance Tripping rate and false detection can be reduced.
For a single-phase system, it is not suitable. [38] SFS and ROCOF Suitable for multi-DG system along with high accuracy and fast detection.
Sometimes allocation of trip boundary is tricky. [39] Rate of change of reactive power (ROCORP) and load connecting strategy (LCS) It has fast detection speed.
Power system stability and quality are the concerns. [40] ROCOF over ROCORP Fast detection speed with high accuracy.
Selection of threshold values are difficult. [41] Combined rate of change of voltage (ROCOV) and ROCORP Small mismatch in power between DG and load can be easily detected.
Power system stability and quality are the concerns. [42] Rate of change of regulator voltage (ROCRV) over ROCORP It can easily detect small mismatch in power among DGs and load.
Selection of threshold values are difficult. [43] Intelligent IDMs

Fuzzy Logic
Suitable for multi-inverter-based DG system. It has good accuracy.
The results are dependent on a set of predefined rulesets. [44] Adaptive network-based fuzzy inference system (ANFIS) [45]  Information can be lost due to less interpretation of independent variables. [54]

Gauss-Newton Algorithm
Due to tidy error estimates, the accuracy is high.
High computation time and implementation cost is also high. [55] Phaselet Algorithm Fast estimation can be achieved by calculating the phasor of variable data.
Due to variable window size unwanted classification occurs during transients. [56] To overcome the problems associated with remote IDMs, local IDMs are introduced, which are categorized into active and passive methods. There are several active IDMs proposed in the literature [19][20][21][22][23][24][25][26][27], which are shown in Table 1. Though for active IDMs, non-detection zones (NDZ) are very small, the power quality of the network is affected due to the injection of disturbance into the system, and their islanding detection speed is slower than passive IDMs [20,[22][23][24][25]. On the other hand, passive IDM is one of the most cost-effective IDMs as the different operating parameters (power, frequency, voltage, current, harmonic distortion, etc.) are measured at the point of common coupling (PCC) to identify the abnormalities in the system [31]. These methods can be applied to both meshed and radial distribution networks consisting of DGs because no disturbance is created by them in the network and they have a faster detection speed than active ones [35]. The shortcoming of these methods is that they have a large NDZ range, in which it is not possible to detect islanding, because any islanding occurrence disturbed signal cannot be recognized by passive IDMs in that range [31].
By combining different active and passive IDMs with each other or with signal processing and artificial intelligence-based IDMs, another method known as hybrid IDMs has been developed to address the problems of passive and active IDMs. These methods have small NDZ and, due to combined characteristics, they can improve the robustness and fault tolerant nature of the distribution network [36]. In addition, this IDM has high accuracy and a fast detection time [39]. However, power quality and system stability remain an issue for this IDM [36,39]. Further, due to the combined feature, the islanding detection time of hybrid IDMs is higher than independent passive and active IDMs [40,42].
All the IDMs discussed earlier either depend on the sensors' pre-defined data limits or on the real-time data, even though most of the time they are accurate, but at certain anomalies, the pre-set data may not be sufficient, and this may lead to an unintended trip [6]. Therefore, advanced IDMs, namely intelligent and signal-processing IDMs, were introduced to overcome the issues of conventional IDMs. The advantages of intelligent IDMs over signal processing IDMs are high precision, good reliability, and a non-requirement of threshold value selection [49]. Even though intelligent IDMs can be easily implemented in an SG-based distribution system, in an inverter-based distribution system, they may not successfully detect islanding [44]. On the other hand, signal processing-based IDMs have a smaller NDZ compared to intelligent IDMs. The common drawbacks of both these methods are high computation and implementations are difficult [48,51].
It can be summarized from the earlier discussion that various methods of IDMs have been developed for detecting islanding efficiently and cost effectively. Nevertheless, these methods still have some limitations, such as complications in implementing the algorithm of IDM, a higher cost associated during the implementation of IDM, higher computation and detection time of the IDM algorithm, higher cost for maintenance and lack of implementation by considering a real-world power system network. As a result, robust IDM development with a simple algorithm and faster response is necessary, which could detect islanding within very short time. In addition, developed IDMs also have less implementation complexity for a practical power system network with a lower NDZ range. In [57], a hybrid method by combining ROCOP and LCS is proposed, which shows a faster detection time. However, the analysis of NDZ for IDM was not considered is this study. This paper aims to present a novel hybrid IDM to detect islanding conditions of DGs, namely a synchronous generator (SG), photovoltaic (PV) generator, and biomass generator. This paper considers a low-voltage distribution system of Malaysia with an 11 kV distribution voltage to evaluate the performances of the proposed hybrid IDM. A combination of two passive IDMs, namely ROCOAP and ROCORP, along with an active LCS method is considered in the proposed novel hybrid IDM. LCS works as active method because LCS checks the necessity of islanding by injecting extra amounts of load to the distribution system. Both the passive IDMs of ROCOAP and ROCORP are considered due the sensitivity of changing the reactive power in the system, as the reactive power changes might cause false islanding detection. Sometimes, IDM fails to detect islanding and non-islanding mode with the small changes in the system's reactive power demand, which might cause inappropriate disconnection of DGs from the distribution network. For this reason, the proposed IDM considers the combination of ROCOAP and ROCORP for proper distinguishing between the necessity of the true islanding state and the false islanding state. The proposed hybrid IDM has a fast islanding detection response as it contains the high-speed passive IDMs. This IDM, with a small detection time and small NDZ, is satisfactory for both multiple DGs and single synchronous DG. In addition, the proposed IDM can also avert short circuit fall. On top of this, the proposed IDM is simple and can be easily implemented in real-world distribution networks in comparison with common passive and active IDMs.
The main contributions of this paper are as follows: A hybrid IDM has been modelled and developed based on the combination of two passive IDMs, namely ROCORP and ROCOAP, and LCS as an active IDM to detect islanding phenomena at PCC. The proposed IDM's performance has been validated in a PSCAD environment for various cases, such as islanding, fault analysis, quality factor, load variations, DG tripping, power mispatch and NDZ range. The proposed IDM is applied at a PCC between DGs and an existing 11 kV Malaysian distribution network, which are modelled using the modules available from the PSCAD library.
Finally, a comparative study has been conducted based on islanding detection time and NDZ to prove the better performance and effectiveness of the proposed IDM.
The organization of the paper is as follows: the proposed IDM's detailed modelling is presented in Section 2. In Section 3, the proposed IDM's NDZ range is presented. The details of the testbed are presented in Section 4. Section 5 presents the simulation results and analysis of the study along with the comparative study. Finally, in Section 6, the outcome of the research is summarized.

Modelling of Proposed Hybrid Islanding Detection Strategy
Active power and reactive power of DGs are represented in Equations (1) and (2), respectively. Here, the active power and reactive power of the grid, load and DGs are considered.
In Equation (1), P DG is active power at the DG side, R L is load resistance, I L is load current, R DG is DG resistance, I DG is the current of the DG, and ∆P is active power difference of load and DG. In Equation (2), Q DG is reactive power at the DG side, X L is load reactance, I L is load current, X DG is reactance of DG, I DG is the current of DG, and ∆Q is the reactive power difference of load and the DG.
This new hybrid method combines power change variations with the connected load variations to obtain the final result, which triggers the output signal for the breakers.

Active Power Method
One of the examples of passive IDM is the active power method. In this case, the power of a three-phase system is measured at PCC and is given by Equation (3).
The phase values of the output voltages are V a , V b , and V c and the phase values of the output currents are I a , I b , and I c . The system will not be islanded when the values of ROCOAP are lower than the threshold point, and at that time, the DGs will be connected with the grid. However, the system will operate in islanding mode while the values of ROCOAP are more than the threshold point. The values of ROCOAP can be obtained by differentiating Equation (3). During the islanding operation mode, the change in the ROCOAP value is significant. IDM modules send signals to DG breakers for disconnecting the DGs by observing the ROCOAP value. When the load's required power matches the supplied power by the DG, fewer oscillations are observed in the output, and for this, the IDM sometimes gives the wrong results. Hence, investigating this phenomenon is necessary to improve the IDM.

Reactive Power Method
The reactive power method is another passive IDM which is like the active power method, except for the sensitivity of ROCORP. The sensitivity of ROCORP is higher than the sensitivity of ROCOAP. System reactive power is monitored continuously for ROCORP. The reactive power equation of ROCORP at the PCC is given by Equation (4).
The values of ROCORP can be obtained at PCC by differentiating Equation (4). While the system is connected to the grid, the oscillation does not cross the threshold point, but it crosses the threshold during islanding, and this is similar to ROCOAP.

IDM Final Stage including ROCOAP, ROCORP and LCS
Three different IDMs are considered to develop the hybrid IDM in this paper. The step-by-step process of implementing these three IDMs is depicted in Figure 1. From the figure, it can be seen that in the beginning, the module first considers active and reactive power responses from the DG side and compares the dataset with the margin limit. If the data cross the margin limit, then the module sends a positive signal to the breakers to trip. The breakers will receive the positive signal when Equations (9) and (14) are veracious. First stage: Given equations are the steps for the flowchart algorithm. Based on these equations, the module will work. p = p measured i , p measured i+1 , p measured i+2 · · · · · · · · · · · · · · · p island (5) q = q measured i , q measured i+1 , q measured i+2 · · · · · · · · · · · · · · · q island (6) First step: dp = dp measured i , dq = dq measured i dp measured i dt > dp dt Island If Equation (7) is not true, then dp = dp measured i+1 dp measured i+1 dt > dp dt Island If Equation (8) is true, then If Equation (9) is true, then the system is islanded. If it is not true, then the LCS will be initiated.
At a specific time, (dp/dt) measured and (dq/dt) measured are the values of ROCOAP and ROCORP, respectively, received by the module. (dp/dt) Island and (dq/dt) Island are the margin values which specify the data crossing limit for the islanding scenario. If reactive power is not more than the reference threshold value, then in the second stage, the LCS is connected for further clarification of islanding and non-islanding events, such as when there is an increase in total load, but the grid is still connected, or any induction motor has started or any capacitor bank is connected. When the power imbalance between the generation and demand is small, the LCS is connected. To accelerate the power change, a suitable value of R-L load is connected. The LCS will be activated when Equation (9) is not satisfied.
The verification of the whole process has been carried out in a PSCAD/EMTDC simulation environment while the IDM module is developed using the FORTRAN language. The required real and reactive power of the load can be calculated using Equations (15) and (16). Islanding will be detected by the module based on the comparison of active and reactive power responses and LCS together.
where P load and Q load are the load's real and reactive powers, respectively. ∆P and ∆Q are the changes in real and reactive power, respectively. P and Q are the total active and reactive powers of the DGs, respectively. The amplitude of the RLC load, including phase angle (φ LOAD ), resonant frequency (ƒ 0 ) and the quality factor (Q f ), are given by Equations (17)- (20), respectively.

Non-Detection Zone (NDZ) of the Proposed Method
In dictating the precision and potency of the IDM, the NDZ plays an indispensable role. The NDZ represents the area or zone or the region in which the IDM will not be proficient enough to perceive islanding. It can be represented as a gap between the power supplied by the grid and the required resistive, capacitive or inductive load by consumers. The smaller the region of the NDZ, the more precise the IDM is.
The active power for NDZ can be calculated using Equation (21) [40]: where ∆P = change in active power, V = estimated voltage, I = estimated current and ∆V = difference between voltage upper limit and lower limit. The reactive power equation for NDZ is given by (22) [40]: where ∆Q = change in reactive power, V = estimated voltage, f n = nominal frequency, ∆f = difference between frequency upper limit and frequency lower limit and Q f = quality factor. Equations (21) and (22) can also be represented as Equations (25) and (26), respectively. (26) where f = estimated frequency, V = estimated voltage, P = active power, Q = reactive power, Q f = quality factor, ∆P = change in active power, ∆Q = change in reactive power, V max , V min = upper limit and lower limit of voltage and f max , f min = upper and lower limit of frequency.
For this study, the maximum and minimum values of the rated voltage are V max = 115% and V min = 85%, and maximum and minimum values of the frequency are f max = 51.5 Hz and f min = 48.5 Hz, and a quality factor Q f = 2.5 is considered according to the IEC 62116 standard [5].
According to IEC 62116, the NDZ region for an active and reactive power imbalance of the distribution system should be in between (0.3841 MW to −0.2439 MW) and (0.1435 MVar and −0.1570 MVar), respectively [5]. The NDZ region for the proposed hybrid IDM is depicted in Figure 2. From Figure 2, it is observed that the proposed technique works efficiently when the real power imbalance is between 0.043 MW and −0.05 MW and the reactive power imbalance is between 0.043 MVar and −0.045 MVar. Thus, the proposed IDM increases the accuracy and also reduces the NDZ.

Testbed under Study
The effectiveness of the proposed IDM mode is tested in this study by developing an 11 kV distribution system that is symmetrical with the Malaysian distribution system. The developed system consists of a synchronous generator (SG), a photovoltaic (PV) generation system, a biomass generator, 32 buses including DG buses, and 29 loads, including LCS. The testbed along with the DGs is modelled in the PSCAD environment. The modelling of the SG, PV and biomass systems is adopted from [40,58,59]. Figure 3 shows the test frame satisfying the IDM operations as per the IEEE 1547 standard recommendations and Figure 4 shows the developed system as per the IEEE 1547 standard. It can be seen from Figure 4 that the grid and LCS breakers can initiate islanding and LCS, respectively, in the testbed system. In addition, the connected load can be designed as various types, such as parallel RLC load, as the IDM performance depends on the load type. It can be concluded from Figure 4 that the developed system considers all the possible considerations mentioned in the IEEE 1547 standard.
The values of R, L and C loads can be calculated using Equations (29)-(31) when the power factor is equal to 1 [40]: where Q f = quality factor.
Values of loads can be calculated by using Equations (32)- (35) when the power factor is not equal to 1 [40]: where X L and X c are the inductive and capacitive reactances. Tables A1 and A2 (Appendix A) exhibit the SG, transformer, PV, biomass and load data of an 11 KV Malaysian distribution system. The operating voltage of the SG and biomass is 3.3 KV, and for PV it is 0.23 kV. A step-up transformer rated at 2 MVA is connected to the DG units to step the voltage up to 11 kV. Three DG units supply 3.6 MW and 1.3 MVar power to the load.

Simulation Results
Various cases including voltage and frequency variations have been considered to verify the performance of the proposed IDM on the PSCAD platform in an 11 kV, 27 bus Malaysian distribution system. Seven case studies are considered to examine the performances of the module under islanding and non-islanding cases to see whether the module can identify islanding cases as islanding and non-islanding cases as non-islanding, or mistakenly take a non-islanding as islanding. The threshold values were set according to the distribution system and DG responses. The module compares the instantaneous result with the given threshold to check the difference that occurs when the system enters the islanding mode, or in the case of other phenomena.

Case 1: Grid Supply Disconnected for Intentional Islanding Operation
Intentional islanding of the grid was performed by disconnecting the grid circuit breaker (CB) to check the module's functionality against the islanding detection, which is shown Figure 5. It can be seen from Figure 5 that the DG with the distribution network falls under the islanding condition due to the disconnection of the grid at t = 3 s. Figure 5 also shows that at t = 3.04 s, dp/dt reaches 4 MW/s, crossing the margin limit of 0.8 MW/s. Now, the module will measure ROCORP according to the flowchart in Figure 1, and the ROCORP DGs at PCC are shown in Figure 6. It can be seen from Figure 6 that the value of dq/dt crosses the margin of 0.019 MVar/s at t = 3.03 s. Therefore, the signal for disconnection of the DGs are sent to all DG breakers, as at t = 3.07, the islanding detection conditions are true. The real power and reactive power of DGs are depicted in Figure 7 at the time of islanding detection.

Case 2: Varying Quality Factors
The performance of the developed IDM was carried out by varying the quality factors (Q f ) according to IEEE Std. 929, IEEE Std. 1547, and UL1741 Std. Various standards encourage different values of quality factors; for example, in the IEEE Std. 1547, the value of Q f = 1, while in the IEEE Std. 929, the value of Q f is ≤2.5, and in the UL1741 Std. the recommended value of Q f is ≤1.8. All these three standards' recommendations are considered while choosing the Q f values to evaluate the developed IDM performance. The corresponding values of R, L, and C for different Q f values are shown in Table 2. As shown in Figure 8, the value of ROCOAP drastically increases at t = 3 s, when DGs entered the islanding mode. After that, when the system checks the ROCORP, from Figure 9 it can be seen that it also crosses the threshold, thus detecting islanding within 0.08 s. Figure 10 represents PCC reactive and real power responses.

Case 3: Initiating the Connection of Varying Reactive Power
For improving the voltage sag and power factor, capacitor banks are connected in parallel with the load, and during this time, the passive parameters of the system change.
To check the module performance during the changes in the passive parameters, in this case, different ratings of capacitor banks (1 MVar, 1.5 MVar and 0.5 MVar) are connected at t = 3 s. As can be seen from Figure 11, after connecting the capacitor banks, the ROCOAP value does not cross the maximum margin of 0.7 MW/s. As a result, according to the algorithm, the islanding is correctly rejected for ROCOAP. However, as can be seen from Figure 12, even though ROCORP crosses the threshold, the case is detected as non-islanding by the module, because, according to the algorithm, if ROCOAP does not cross the threshold limit, then the module will not check for the ROCORP condition. This indicates that the proposed IDM is also able to detect non-islanding cases accurately. The proposed IDM results for the selected Q c are presented in Table 3, and reactive and real power for different Q c values are shown in Figure 13.

Case 4: Fault Analysis
In this section, the performance of the proposed IDM's performance is validated by applying different types of faults. The threshold limits of ROCOAP and ROCORP are set to 0.50 MW/s and 0.1 MVar/s, respectively, during fault conditions. For this case, the load values of active and reactive powers are considered to be 1.1 MW and 0.79 MVar, respectively. The ROCOAP threshold is set to 1 MW/s after connecting the LCS because of the addition of new loads. Under different fault conditions, the detection performances of the proposed IDM are presented in Table 4. Different types of faults (L-G, L-LG, L-L-L-G, and L-L) are applied at the bus 20 Malaysian distribution network at t = 3.02 s. The ROCOAP and ROCORP responses at PCC are shown in Figures 14-17 for fault resistances of 0.01 (Ω) and 0.02 (Ω). From the figures, it can be seen that at t = 3.06 s, the ROCOAP value is more than the threshold 0.50 MW/s. In contrast, the ROCORP value does not exceed the set value of 0.1 MVar/s, because due to the applied fault, the threshold values are changed. Now, according to the proposed algorithm, the LCS started operating and the ROCOAP response during LCS operation is presented in Figures 18 and 19 for fault resistances of 0.01 (Ω) and 0.02 (Ω), respectively. From the figures, it is observed that the ROCOAP value does not exceed the threshold limit of 1 MW/s after connecting the LCS. Since ROCOAP has not exceeded the minimum acceptable limit within the specified period, the system continues to operate by rejecting islanding. The real and reactive power responses of DG units during islanding are depicted in Figure 20.

Case 5: Starting of Induction Motor
To show that, at the time of induction, the motor starting the module properly perceives the islanding state, in this section, the performance of the proposed IDM is validated by connecting the induction motor at t = 3 s to the network. The ROCOAP and ROCORP threshold values are set to 0.4 MW/s and 0.39 MVar/s, respectively, for this case. The load values of the active and reactive powers are set to 1 MW and 0.8 MVar, respectively.
According to Figure 21, ROCOAP has crossed the threshold of 0.4 MW/s at t = 3 s, whereas the ROCORP value has not crossed the threshold of 0.39 MW/s, which is shown in Figure 22. As a result, the module initiates LCS, and after connecting LCS, it can be seen from Figure 23 that the ROCOAP value has exceeded the threshold limit. On the other hand, after connecting the LCS, ROCORP has not exceeded the threshold of 0.39 MVar/s according to Figure 24, which confirms that the case has accurate non-islanding case detection. Therefore, it is validated that by using the proposed algorithm, non-islanding cases can be detected, and the system continues to work in grid-connected mode with a lagging power factor. DG units' power (real and reactive) responses during islanding are depicted in Figure 25.

Case 6: Zero Power Mismatch
The performances of the proposed IDM are validated during zero active, reactive and total power mismatches in this case study.

Zero Active Power Mismatch
To create a zero active power mismatch scenario, the real power flow between the entire DG systems and the grid is maintained to a zero value at t = 3 s by disconnecting the grid. Figure 26 shows that at t = 3.04 s, the ROCOAP value exceeds the threshold of 0.7 MW/s. As a result, the event is detected as an islanding condition by the module within a three-cycle time period.

Zero Reactive Power Mismatch
To create a zero active power mismatch scenario, the reactive power flow between the entire DG system and the grid is maintained to a zero value at t = 3 s by disconnecting the grid. Figure 27 shows that at t = 3.04 s, the ROCORP value is 1.2 MVar/s, which exceeds the threshold of 0.08 MVar/s. As a result, islanding conditions are detected accurately by the module and the signals to all DG breakers are sent for disconnection at t = 3.04 s.

Case 7: DG Tripping
In the DG-connected distribution system, the tripping of a DG is a very serious issue which can lead to false islanding detection. To verify the performance of the proposed IDM for this case, the PV and biomass are tripped at t = 3 s. For PV and biomass, it is an islanding condition; however, SG is not islanded. From Figure 29, it can be seen that the ROCOAP value increases suddenly up to 0.49 MW/s at t = 3.017 s, but still below the threshold of 0.8 MW/s. As a result, the case is identified as a precise non-islanding case. On the other hand, according to Figure 30, the ROCORP value at t = 3.02 s is found to have exceeded the threshold of 0.019 MVar/s. Even though ROCORP crosses the threshold, the case is detected as non-islanding by the module, because, according to the algorithm, if ROCOAP does not cross the threshold limit, the module will not check for ROCORP. This indicates that the proposed IDM is able to detect both islanding and non-islanding cases accurately.

Case 7: DG Tripping
In the DG-connected distribution system, the tripping of a DG is a very serious issue which can lead to false islanding detection. To verify the performance of the proposed IDM for this case, the PV and biomass are tripped at t = 3 s. For PV and biomass, it is an islanding condition; however, SG is not islanded. From Figure 29, it can be seen that the ROCOAP value increases suddenly up to 0.49 MW/s at t = 3.017 s, but still below the threshold of 0.8 MW/s. As a result, the case is identified as a precise non-islanding case. On the other hand, according to Figure 30, the ROCORP value at t = 3.02 s is found to have exceeded the threshold of 0.019 MVar/s. Even though ROCORP crosses the threshold, the case is detected as non-islanding by the module, because, according to the algorithm, if ROCOAP does not cross the threshold limit, the module will not check for ROCORP. This indicates that the proposed IDM is able to detect both islanding and non-islanding cases accurately.

Based on Islanding Detection Time
In this section, the proposed IDM performance is compared with the other hybrid methods in [36,[41][42][43] based on the detection time, which is presented in Table 5. From the table, it can be seen that the combined voltage and reactive power shift IDM [36] took 160 ms to detect islanding. On the other hand, the ROCOF and ROCORP combined IDM [41], similar to the proposed method, has a detection time of 200 ms, which is higher than the proposed method. The combination of a ROCOV and ROCORP-based hybrid IDM [42] has a detection time of 250 ms, which is higher than both the proposed method and the method in [41]. Finally, it is observed that the ROCORV and ROCORP combined hybrid IDM [43] has a longer detection time (640 ms) than previous methods. From the above discussion, it is clear that even though all the hybrid IDMs were able to detect the islanding within 2 s, which satisfies the IEEE 1547 standards, the proposed hybrid IDM's detection time is faster than others.    [36,[41][42][43] based on the NDZ that is presented in Table 6. The comparison results of different IDMs based on only NDZ is presented in Figure 31. From Figure 31, it can be seen that the proposed IDM has a smaller NDZ in comparison with other IDMs, and islanding is recognized satisfactorily using the proposed hybrid IDM. Further, the proposed IDM is not negatively influenced by any kinds of faults.

Conclusions
This work has presented a novel hybrid IDM based on the estimation of power changes and load fluctuation in the distribution system. The performance of the proposed IDM has validated for different islanding and non-islanding cases in the PSCAD/EMTDC platform, whereas a testbed Malaysian distribution network (11 kV) has also been considered. From the simulation results, the efficacy of the proposed IDM is perceived because within 0.09 s, the proposed IDM was able to detect islanding, which is lower than 2 s, which is in line with the IEEE standard 1547. In addition, for a distribution network consisting of multiple and different DGs (SG, PV and biomass), the proposed IDM has also shown excellent performance by successfully differentiating between islanding and non-islanding cases.
A comparative study based on the detection period between the proposed IDM and three different hybrid detection methods to prove the better performance of the proposed method has been carried out in this paper. It can be seen that the detection time taken by the proposed IDM is 0.09 s, which is 0.07 s, 0.11 s, 0.16 s and 0.55 s less than the other hybrid IDMs, namely the combined voltage and reactive power shift IDM, ROCOF and ROCORP IDM, ROCOV and ROCORP IDM, and ROCORV and ROCORP IDM, respectively. Furthermore, it is observed that the proposed IDM has a smaller NDZ compared to other available hybrid IDMs.
In this work, through simulation only the performance of the proposed IDM has been validated in an existing Malaysian network. Therefore, to prove its effectivity also in real-world conditions in the future, the performance of the proposed IDM module will be verified through either hardware-in-loop simulations or hardware implementation.
Author Contributions: Conceptualization, methodology, and simulation, U.K.J. and S.A.; manuscript writing and editing, U.K.J., S.A. and T.A. All authors have read and agreed to the published version of the manuscript.

Funding:
The authors would like to thank University of Malaya, Malaysia and Chittagong University of Engineering and Technology, Chittagong, Bangladesh for providing research facilities.

Data Availability Statement:
The data used in this research is enclosed within the manuscript. No external data sources are used in this research.

Conflicts of Interest:
The authors declare no conflict of interest.