Design, Evaluation and Implementation of an Islanding Detection Method for a Micro-grid

: Correct and fast detection of a micro-grid (MG) islanding is essential to the MG since operation, control, and protection of the MG depend on its operating mode i.e., an interconnected mode or islanding mode. This study describes the design, evaluation and implementation of an islanding detection method for an MG, which includes a natural gas-ﬁred generator, a doubly fed induction generator type wind generator, a photovoltaic generator, and some associated local loads. The proposed method is based on the instantaneous active and reactive powers at the point of common coupling (PCC) of the MG. During the islanding mode, the instantaneous active and reactive powers at the PCC are constants, which depend on the voltage of the PCC and the impedance of the dedicated line. The performance of the proposed method is veriﬁed under various scenarios including islanding conditions for the different outputs of the MG, and fault conditions by varying the position, type, inception angle and resistance of the fault, using the PSCAD/EMTDC simulator. This paper also concludes by implementing proposed method into a TMS320C6701 digital signal processor. The results indicate that the proposed method successfully detects islanding for the MG in islanding conditions, and remains stable in fault conditions.


Introduction
Distributed generations (DGs) are small-scale power generations, which are considered to be a promising approach as a solution of economic and environmental issues for conventional power systems [1][2][3]. DGs provide economic advantages to reduce the amount of energy lost in transmitting electricity, and the number and size of power distribution lines, since the electricity can be generated near where it is consumed. In addition, the integration of distributed renewable energy sources i.e., wind energy, solar energy, bio-energy, hydraulic energy and so on, into the traditional electric power system, can reduce the emission of the greenhouse gases, and thus solve the environmental problems. However, the limitation of integrating DGs is that, they must be obedient to the operating conditions of the interconnected distribution networks. When the grid breaks the connection between the distribution network and DGs, DGs must detect the islanding condition and immediately stop producing power, which is called anti-islanding. Otherwise, grid operators may not realize that a circuit is still powered by DGs, and automatic re-connection of devices may be prevented. This will reduce the capacity factor of DGs and increase the power lost in transmitting electricity, since the loads which were powered by DGs would be powered by other remote power sources.
impedance of the fault, are analyzed using the PSCAD/EMTDC generated data. Lastly, this paper concludes by implementing the proposed method into a TMS320C6701 digital signal processor (DSP).

Calculation of the Instantaneous Active and Reactive Power
The three-phase instantaneous active power (p 3ph ) delivered from the MG to the dedicated line, which is defined and given in (1), can be calculated from the voltages (v a , v b and v c ) and currents (i a , i b and i c ) measured at the PCC.
The instantaneous reactive power could be calculated by taking the imaginary part of the complex power, whose real part represents the instantaneous active power [21,22]. However, this calculation of the instantaneous reactive power is invalid under unbalanced operating conditions [23,24]. To correctly calculate the instantaneous reactive power even under an unbalanced fault condition, the instantaneous reactive power (q 3ph ) delivered from the MG to the dedicated line is calculated from the voltages (v a ', v b ' and v c '), which respectively lag v a , v b and v c by a quarter of period, and currents (i a , i b and i c ) [25], as shown in (2).

Islanding Detection Method Based on the Instantaneous Active and Reactive Power at the PCC
In this subsection, the proposed islanding detection method for an MG based on p 3ph and q 3ph is described. When the MG is disconnected from the distributed network by opening the circuit breakers at the grid side, p 3ph and q 3ph will be constant, which depend on the voltages at the PCC as well as the impedance of the dedicated line. p 3ph will become almost zero because little resistance exists in the dedicated line. On the other hand, q 3ph has some value corresponding to the series inductance and the shunt capacitance of the dedicated line. Since the parameters of the dedicated line can be obtained, these constant active and reactive powers can be easily calculated. Therefore, if the calculated instantaneous active and reactive powers converge to the pre-calculated constant values, the islanding inception will be detected. Figure 1 shows the islanding detection region. The reference value for the complex power (S ref ) is given by where, V PCC and Z line * represent the rated line-line voltage of the PCC and the impedance of the dedicated line, respectively.

Calculation of the Instantaneous Active and Reactive Power
The three-phase instantaneous active power (p3ph) delivered from the MG to the dedicated line, which is defined and given in (1), can be calculated from the voltages (va, vb and vc) and currents (ia, ib and ic) measured at the PCC.
The instantaneous reactive power could be calculated by taking the imaginary part of the complex power, whose real part represents the instantaneous active power [21,22]. However, this calculation of the instantaneous reactive power is invalid under unbalanced operating conditions [23,24]. To correctly calculate the instantaneous reactive power even under an unbalanced fault condition, the instantaneous reactive power (q3ph) delivered from the MG to the dedicated line is calculated from the voltages (va', vb' and vc'), which respectively lag va, vb and vc by a quarter of period, and currents (ia, ib and ic) [25], as shown in (2).

Islanding Detection Method Based on the Instantaneous Active and Reactive Power at the PCC
In this subsection, the proposed islanding detection method for an MG based on p3ph and q3ph is described. When the MG is disconnected from the distributed network by opening the circuit breakers at the grid side, p3ph and q3ph will be constant, which depend on the voltages at the PCC as well as the impedance of the dedicated line. p3ph will become almost zero because little resistance exists in the dedicated line. On the other hand, q3ph has some value corresponding to the series inductance and the shunt capacitance of the dedicated line. Since the parameters of the dedicated line can be obtained, these constant active and reactive powers can be easily calculated. Therefore, if the calculated instantaneous active and reactive powers converge to the pre-calculated constant values, the islanding inception will be detected. Figure 1 shows the islanding detection region. The reference value for the complex power (Sref) is given by where, VPCC and Zline* represent the rated line-line voltage of the PCC and the impedance of the dedicated line, respectively. The criteria of the islanding detection are given by where k 1 and k 2 depend on the variation of the voltages at the PCC and the measurement ratio errors of the current transformer (CT) and the potential transformer (PT). The variation of the voltages at the PCC is ±20%, with the fully consideration of the voltage deviation both in the steady-state and the transient state after an islanding inception, and sufficient margin. Thus, k 1 and k 2 are set by where 1.0, 0.8, 1.2, ±3% and ±6% are respectively per unit value of the rated voltage, lower and upper limits of voltage variation, limits of CT and PT ratio errors [26,27]. In (6), the minimal magnitude of the complex power at the PCC in islanding conditions is calculated when the actual voltages at the PCC are only 80% of the rated voltages considering maximal voltage variation, and the measured currents and voltages at the PCC are respectively 97% and 94% of the actual currents and voltages considering the maximal ratio errors of CTs and PTs. Meanwhile, as shown in (7), the maximal magnitude of the complex power at the PCC in islanding conditions is calculated when the actual voltages at the PCC are 120% of the rated voltages, and the measured currents and voltages at the PCC are respectively 103% and 106% of the actual currents and voltages. k 3 is set to be 15 • , considering the limits of phase errors of both CTs and PTs, also with sufficient margin [26,27]. These three coefficients, k 1 , k 2 and k 3 , only depend on the limits of voltage variation and the limits of measurement errors of CTs and PTs defined in the IEC Standards [26,27]. Thus, when the proposed method is applied in another MG, k 1 , k 2 and k 3 will not be changed, and only S ref should be pre-calculated considering the parameters of the new dedicated line. As described above, the islanding will be detected when the trajectory of the point (p 3ph , q 3ph ) moves into the islanding detection region, which is used to prevent the mal-detection for the transient disturbances.

Case Studies
An MG, which is connected to a 22.9 kV, 60 Hz distribution network through a Y-Y transformer (22.9/6.6 kV) and a dedicated line (1 km), is shown in Figure 2. The dedicated line is modeled ACSR 58 mm 2 , whose series resistance, series inductance and shunt capacitance are 0.8316 Ω/km, 0.0022 H/km and 0.0021 µF/km, respectively. The MG is composed of three DGs i.e., a 2 MW natural gas-fired generator, a 2 MW DFIG type wind generator, and a 1 MW photovoltaic generator, and some associated local loads. The system is modeled using the PSCAD/EMTDC simulator, where the sampling rate is 32 samples/cycle. The signals of the currents and voltages measured at the PCC are passed through anti-aliasing RC low-pass filters with the cutoff frequency of 960 Hz, which is half the sampling frequency.
The performance of the proposed islanding detection method is verified under various islanding conditions for the different outputs of the MG, as well as various fault conditions varying the fault position, fault inception angle, fault type and fault impedance, as shown in Tables 1 and 2. When an islanding incepts, the proposed islanding detecting method should activate the relay as soon as possible; on the contrary, the proposed method should not activate the relay in fault conditions.    Case 8 Dedicated line 0° Single line-toground 0 Ω Figure 10 Case 9 Dedicated line 0° Double lineto-ground 0 Ω Figure 11 Case 10 Dedicated line 0° Line-to-line 0 Ω Figure 12 Case 11 Dedicated line 0° Single line-toground 1 Ω Figure 13 Case 12 Dedicated line 0° Single line-toground 5 Ω Figure 14

Islanding Conditions
The performance of the proposed islanding detection method is verified under various islanding conditions for the different outputs of the MG. The results of three cases, where the generating power of the MG is smaller than (Case 1), close to (Case 2) and larger than (Case 3) the local loads, are shown in     Case 8 Dedicated line 0 • Single line-to-ground 0 Ω Figure 10 Case 9 Dedicated line 0 • Double line-to-ground 0 Ω Figure 11 Case 10 Dedicated line 0 • Line-to-line 0 Ω Figure 12 Case 11 Dedicated line 0 • Single line-to-ground 1 Ω Figure 13 Case 12 Dedicated line 0 • Single line-to-ground 5 Ω Figure 14

Islanding Conditions
The performance of the proposed islanding detection method is verified under various islanding conditions for the different outputs of the MG. The results of three cases, where the generating power of the MG is smaller than (Case 1), close to (Case 2) and larger than (Case 3) the local loads, are shown in      Figure 3a indicates the voltages and currents measured at the PCC, where v a (solid), v b (dashed), and v c (dotted) are shown in the upper subfigure and i a (solid), i b (dashed), and i c (dotted) are shown in the lower subfigure. After the islanding inception, the voltages decrease slightly whilst the currents decrease to nearly zero. p 3ph and q 3ph are respectively shown in the upper subfigure and lower subfigure of Figure 3b. After islanding inception, p 3ph becomes almost zero after a slight fluctuation whilst q 3ph becomes a very small value directly. Both p 3ph and q 3ph change and become stable again in several milliseconds after the island incepts. In Figure 3c, the trajectory of the point (p 3ph , q 3ph ), which is shown with the marks of "o", is located in the third quadrant prior to the islanding inception. This is because the power is delivered from the grid to the MG prior to the islanding inception, which can also be confirmed in Figure 3b. When the islanding incepts, p 3ph is nearly zero and q 3ph has the negative value because of the characteristics of the dedicated line. Therefore, the trajectory moves to the islanding detection region. In Figure 3d, where "0" and "1" respectively mean mode of interconnection and islanding, the islanding detection signal is activated at 22.53 ms after the islanding inception. The results indicate that the proposed method can successfully and quickly detect the islanding operation in 1.5 cycles after the islanding inception. Figures 4 and 5 show the results for Case 2 and 3, where the generating power is close to and larger than the local loads of the MG, respectively. In both cases, an islanding incepts at 33.33 ms. In Case 2, the power transmitted between the MG and the grid prior to the islanding inception is nearly zero. Due to the little variation of the transmitted power in the dedicated line prior to and after the islanding inception, the magnitude and the phase angle of the voltage and the frequency measured at the PCC do not change significantly. However, the proposed method can discriminate the islanding inception from normal load variations by considering the trajectory of the point (p 3ph , q 3ph ). As seen in Figure 4d, the trajectory of the point moves into islanding detection region at 20.45 ms after islanding inception. In Case 3, the generating power of the MG is larger than the local loads of the MG. Thus, the point (p 3ph , q 3ph ) is in the fourth quadrant prior to the islanding inception. As expected, the trajectory of the point (p 3ph , q 3ph ) enters the islanding detection region from the fourth quadrant at 22.53 ms after the islanding inception, as shown in Figure 5c,d.
The results for Cases 1-3 indicate that the proposed method can successfully detect the islanding operation irrespectively of the relationship between the generating power of the MG and its local loads. In addition, the detection speed is much faster than that in [10], almost 1.5 cycle after the islanding incepts.

Fault Conditions
The performance of the proposed islanding detection method is also verified under various fault conditions varying the position and type of the fault. In addition, the qualitative analysis about the effect of the fault inception angle and fault resistance is also included in this subsection. In all these fault cases, the power transmitted between the MG and the grid prior to the fault inception is nearly zero, which is the same as that in Case 2. All faults occur at 33.33 ms, and the proposed method should not activate the islanding detection signal in fault conditions.

Faults with Different Position
In this subsection, three-phase faults, whose inception angles are all 0 • , with different position i.e., distribution line of the MG (Case 4) and dedicated line (Case 5) are considered. The results are shown in Figures 6 and 7, respectively. Figure 6 shows the results for Case 4. In this case, it is assumed that a three-phase fault occurs at the distribution line in the MG. As shown in Figure 6a, the voltages decrease whilst the currents increase sharply when the fault occurs. From Figure 6b,c, both p 3ph and q 3ph fluctuates after the fault incepts due to the large fault current. The point (p 3ph , q 3ph ) is near the origin prior to the fault inception, since the generating power of the MG is same as the local loads. However, when the fault occurs, the point moves far away the origin, since the fault currents are considerably large. The proposed method does not activate the islanding detection signal (Figure 6d). Figure 7 shows the results for Case 5, where a three-phase fault occurs at the dedicated line. The results are similar with those for Case 4. The trajectories of the point (p 3ph , q 3ph ) do not move into the islanding detection region, as shown in Figure 7c. Thus, the islanding detection signal is not activated.
Energies 2018, 11, x FOR PEER REVIEW 10 of 24 point moves far away the origin, since the fault currents are considerably large. The proposed method does not activate the islanding detection signal (Figure 6d). Figure 7 shows the results for Case 5, where a three-phase fault occurs at the dedicated line. The results are similar with those for Case 4. The trajectories of the point (p3ph, q3ph) do not move into the islanding detection region, as shown in Figure 7c. Thus, the islanding detection signal is not activated.
The results indicate that the proposed method does not activate the islanding detection signal no matter where a fault occurs. The results indicate that the proposed method does not activate the islanding detection signal no matter where a fault occurs.

Faults with Different Inception Angle
In this subsection, three different fault inception angles i.e., 0 • (Case 5), 45 • (Case 6) and 90 • (Case 7) are compared. The fault occurs at the dedicated line at 33.33 ms. The results of Cases 6 and 7 are shown in Figures 8 and 9. From the results and analysis in the previous subsection, the considerably large fault current is why the point (p 3ph , q 3ph ) moves far away the islanding detection region and consequently inactivates the islanding detection signal. In case of faults with different inception angle, even the magnitudes and waveforms of fault current in Cases 6 and 7 were different from those in Case 5, as seen in Figures 8a and 9a, the fault currents were still considerably large. Finally, the islanding detection signal could not be activated as shown in Figures 8d and 9d.
It can be concluded that the proposed islanding detection method can remain stable no matter when a fault occurs.

Faults with Different Type
In this subsection, four kinds of fault types i.e., single line-to-ground (SLG, Case 8), double line-to-ground (DLG, Case 9), line-to-line (LL, Case 10), and three-phase (3P, Case 5) are considered. The fault occurs at the dedicated line, which is same with that of Case 5. The results of Cases 8−10 are shown in Figures 10-12. signal could not be activated as shown in Figures 8d and 9d.
It can be concluded that the proposed islanding detection method can remain stable no matter when a fault occurs.

Faults with Different Type
In this subsection, four kinds of fault types i.e., single line-to-ground (SLG, Case 8), double line-to-ground (DLG, Case 9), line-to-line (LL, Case 10), and three-phase (3P, Case 5) are considered. The fault occurs at the dedicated line, which is same with that of Case 5. The results of Cases 8−10 are shown in Figures 10-12. Figure 10 shows the results for Case 8, where an A-phase SLG fault occurs. From Figure 10a, the A-phase voltage decreases significantly since the fault occurs very close to the PCC, whilst the voltages of other phases do not change much. The A-phase current increases significantly whilst the currents of other phases increase slightly, compared with the faulted phase, due to the zero-sequence component of the fault current. Due to the unbalanced three-phase voltages and currents, p3ph and q3ph fluctuate even when the transient state is finished, as seen in Figure 10b. Hence, the trajectory of the point (p3ph, q3ph) cannot remain stable at one point and move into the islanding detection region in Figure 10c. In addition, the islanding detection signal is not activated.  Figures 11b and 12b, and the trajectory of the point (p3ph, q3ph) does not move into the islanding detection region in Figures 11c and 12c. As expected, the islanding detection signal is not activated.
It can be concluded that the proposed islanding detection method can remain stable in fault conditions irrespective of the type of fault.

Faults with Different Fault Impedance
In this subsection, three SLG faults, whose inception angles are all 0°, with different fault impedance i.e., 0 Ω (Case 5), 1 Ω (Case 11) and 5 Ω (Case 12), are analyzed together. The results for Cases 11 and 12 are shown in Figures 13 and 14, respectively.
From the analysis in the previous subsection, when an unbalanced fault i.e., SLG, DLG and LL fault occurs, both p3ph and q3ph fluctuate and the trajectory of the point (p3ph, q3ph) cannot remain stable at one point and move into the islanding detection region. The similar results and analysis can be easily drawn even the fault resistance exists in these unbalanced fault conditions. In addition, fault resistance has no effect on the magnitude and waveform of the fault current in the case of a balanced fault (3P fault). Hence, it could be easily concluded that the fault resistance would not affect the performance of the proposed islanding detection method.
The results for all fault cases indicate that the proposed islanding detection method does not activate the islanding detection signal under various fault scenarios considering different position, inception angle, type and impedance of fault. Hence, the proposed islanding detection method can remain stable as expected in fault conditions, irrespectively of the position, type, inception angle and resistance of a fault.   Figure 10a, the A-phase voltage decreases significantly since the fault occurs very close to the PCC, whilst the voltages of other phases do not change much. The A-phase current increases significantly whilst the currents of other phases increase slightly, compared with the faulted phase, due to the zero-sequence component of the fault current. Due to the unbalanced three-phase voltages and currents, p 3ph and q 3ph fluctuate even when the transient state is finished, as seen in Figure 10b. Hence, the trajectory of the point (p 3ph , q 3ph ) cannot remain stable at one point and move into the islanding detection region in Figure 10c. In addition, the islanding detection signal is not activated. Figures 11 and 12 show the results for Cases 9 and 10, where a DLG fault and an LL fault occurs at the dedicated line, respectively. Similar with the results for Case 8, p 3ph and q 3ph fluctuate in Figures 11b and 12b, and the trajectory of the point (p 3ph , q 3ph ) does not move into the islanding detection region in Figures 11c and 12c. As expected, the islanding detection signal is not activated.
It can be concluded that the proposed islanding detection method can remain stable in fault conditions irrespective of the type of fault.

Faults with Different Fault Impedance
In this subsection, three SLG faults, whose inception angles are all 0 • , with different fault impedance i.e., 0 Ω (Case 5), 1 Ω (Case 11) and 5 Ω (Case 12), are analyzed together. The results for Cases 11 and 12 are shown in Figures 13 and 14, respectively.
From the analysis in the previous subsection, when an unbalanced fault i.e., SLG, DLG and LL fault occurs, both p 3ph and q 3ph fluctuate and the trajectory of the point (p 3ph , q 3ph ) cannot remain stable at one point and move into the islanding detection region. The similar results and analysis can be easily drawn even the fault resistance exists in these unbalanced fault conditions. In addition, fault resistance has no effect on the magnitude and waveform of the fault current in the case of a balanced fault (3P fault). Hence, it could be easily concluded that the fault resistance would not affect the performance of the proposed islanding detection method. The results for all fault cases indicate that the proposed islanding detection method does not activate the islanding detection signal under various fault scenarios considering different position, inception angle, type and impedance of fault. Hence, the proposed islanding detection method can remain stable as expected in fault conditions, irrespectively of the position, type, inception angle and resistance of a fault.

Hardware Implementation
Practically, noise signals are contained in real measured three-phase voltage and current signals. As the currents flowing through the PCC in an islanding condition are extremely small, the effect of noise signals of the measured voltages and currents on the performance of the proposed method might not be ignored. Therefore, to verify the performance of the proposed method when noise signals are contained, the method is tested under practical conditions and this section shows the results of hardware implementation of the method into a TMS320C6701 DSP. Figure 15 shows the configuration of hardware implementation. Three-phase voltages and currents generated by PSCAD/EMTDC simulator are converted into analog signals using PCI 1724 U board and then injected into the Intelligent Electronic Device (IED) based on a TMS320C6701 DSP. The signals are then passed through first-order RC filter (f c = 960 Hz) to the 16-bit A/D converters operating at a sampling rate of 32 s/c. All calculation and process of islanding detection are done in the IED. Figures 16 and 17 show the results of Case 1, in which islanding incepts at 33.33 ms and the islanding detection signal should be activated, and Case 4, in which a 3P fault occurs at 33.33 ms and the islanding detection signal should not be activated. As shown in Figure 16, the point (p 3ph , q 3ph ) cannot remain stable at one point even when the transient state is over. This is because p 3ph and q 3ph slightly fluctuate due to the noise signals in the voltages and currents. To prevent mal-operation due to these noise signals, the islanding detection region is appropriately expanded and set to be a circle.
The results indicate that the proposed method can successfully and fast detect the islanding inception at 17.71 ms after the islanding inception. In Figure 17, even noise signals are contained in real voltage and current signals, the proposed islanding detection method does not activate the islanding detection signal due to large p 3ph and q 3ph after the fault inception.

Hardware Implementation
Practically, noise signals are contained in real measured three-phase voltage and current signals. As the currents flowing through the PCC in an islanding condition are extremely small, the effect of noise signals of the measured voltages and currents on the performance of the proposed method might not be ignored. Therefore, to verify the performance of the proposed method when noise signals are contained, the method is tested under practical conditions and this section shows the results of hardware implementation of the method into a TMS320C6701 DSP. Figure 15 shows the configuration of hardware implementation. Three-phase voltages and currents generated by PSCAD/EMTDC simulator are converted into analog signals using PCI 1724 U board and then injected into the Intelligent Electronic Device (IED) based on a TMS320C6701 DSP. The signals are then passed through first-order RC filter (fc = 960 Hz) to the 16-bit A/D converters operating at a sampling rate of 32 s/c. All calculation and process of islanding detection are done in the IED. Figures 16 and 17 show the results of Case 1, in which islanding incepts at 33.33 ms and the islanding detection signal should be activated, and Case 4, in which a 3P fault occurs at 33.33 ms and the islanding detection signal should not be activated. As shown in Figure 16, the point (p3ph, q3ph) cannot remain stable at one point even when the transient state is over. This is because p3ph and q3ph slightly fluctuate due to the noise signals in the voltages and currents. To prevent mal-operation due to these noise signals, the islanding detection region is appropriately expanded and set to be a circle.
The results indicate that the proposed method can successfully and fast detect the islanding inception at 17.71 ms after the islanding inception. In Figure 17, even noise signals are contained in real voltage and current signals, the proposed islanding detection method does not activate the islanding detection signal due to large p3ph and q3ph after the fault inception.

Conclusions
This paper proposes an islanding detection method for the MG based on the instantaneous active and reactive power delivered from the MG to the dedicated line. The instantaneous active and

Conclusions
This paper proposes an islanding detection method for the MG based on the instantaneous active and reactive power delivered from the MG to the dedicated line. The instantaneous active and reactive power are calculated and used to monitor whether the islanding incepts or not. When the circuit breakers at the grid side open, the monitored instantaneous active power from the MG to the dedicated line is converged to nearly zero, whilst the instantaneous reactive power from the MG to the dedicated line has some small value according to the shunt capacitance and the series inductance of the dedicated line. Therefore, the trajectory of the point (p 3ph , q 3ph ) moves into the islanding detection region, and the islanding detection signal is consequently activated. On the contrary, the trajectory of the point (p 3ph , q 3ph ) would move to another point or fluctuates in a fault condition. The islanding detection region can be pre-defined considering the parameters of the dedicated line, variation of the voltage, and the possible measurement errors of CTs and PTs.
The performance of the proposed islanding detection method is verified under various islanding conditions for the different outputs of the MG, as well as various fault conditions varying the position, type, inception angle and resistance of the fault, with the PSCAD/EMTDC generated data. The results indicate that the proposed method can successfully and fast detect the islanding operation irrespective of the relationship between the generating power and the local loads. In addition, the proposed method does not mal-operate irrespectively of the position, type inception angle and resistance of the fault. A prototype relay based on the described scheme successfully detects islanding inception.
The proposed method can correctly and quickly detect the islanding inception. Consequently, the strategies of operation and control of the MG could be re-decided. In addition, the threshold values for the protection relays in the MG could be properly re-set according to the mode of islanding or not, to increase the reliability of the protection system of the MG.