A New Protection Scheme For Feeders of Microgrids With Inverter-Based Resources
A New Protection Scheme For Feeders of Microgrids With Inverter-Based Resources
A R T I C L E I N F O A B S T R A C T
Keywords: Due to the unavailability of a suitable relay for microgrid protection, various utilities are fitting available IEDs for
Microgrid protection the protection of microgrids. At present, microgrid protection is achieved using a combination of conventional
Distributed energy resources numerical relays. These numerical relays are not suitable for all kinds of microgrid architectures and do not
Digital relays
provide complete protection with inverter-based generators. Since these relays were designed considering the
Phasor measurement unit
fault characteristics of synchronous generators, they fail to respond to the fault characterized by inverter-based
generators. This paper proposes a new protection technique that is independent of the type of generating sources,
control philosophy of inverters as well as microgrid architecture. Simulations are performed using PSCAD/
EMTDC and performance of the protection scheme is also evaluated in real-time using RTDS. Pertinent results are
presented which demonstrate the effectiveness of the proposed scheme.
* Corresponding author.
E-mail address: [email protected] (T.S. Sidhu).
https://doi.org/10.1016/j.epsr.2023.109632
Received 9 December 2022; Received in revised form 10 March 2023; Accepted 21 June 2023
Available online 16 August 2023
0378-7796/© 2023 Elsevier B.V. All rights reserved.
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
2
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Fig. 6. Discrepant impedances for feeder between bus-4 and bus- 5 (a) before
the fault F1 and (b) during the fault at F-1(1-Φfault at 30% of line length). Fig. 7. Discrepant impedances for healthy feeders during a fault at F-1.
other electrical parameters are derived from these values. The proposed 2.1. No fault condition
scheme develops a new fault detection philosophy by estimating the
positive sequence ‘discrepant impedance’ of the feeder. The voltage and To explain the calculation of discrepant impedance, a two bus
current samples are obtained from two ends of the feeder. From these microgrid with both ends connected with inverter-interfaced distributed
sampled values, positive sequence ‘discrepant impedance’ is calculated generators is shown in Fig. 1. For a shunt fault at f at a distance x from
by the relay. The detailed theory and basis are described in the following Bus A, let the positive sequence voltages and currents at Bus-A and Bus-B
sections. be VA1, IA1, VB1, and IB1 respectively. Z1 is the positive-sequence
impedance of the feeder between Bus-A and Bus-B. The voltage VB1 at
Bus-B when current IA1 is flowing from Bus-A to Bus-B is: Fig. 1.
Sequence
3
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Then,
̅→ ̅→
→ VB1 − VA1
− Z1 = ̅→ (4)
IB1
-Z1 is the positive-sequence impedance of the feeder as estimated
form Bus-B using positive-sequence voltages from Bus-A and Bus-B, and
positive-sequence current measured at Bus-B.
Fig. 8. Discrepant Impedances during a fault F-2 (L-L fault at 60% of line Under no fault condition, IA1 is equal to IB1. Therefore, during no
length) on feeder 4–5. fault condition it can be observed from (2) and (4), the summation of the
̅→ ̅→
impedances Z1 and − Z1 is zero. This summation of impedances is
termed as ‘discrepant impedance’. This impedance signifies the
discrepancy between the feeder impedances when estimated from one
end of the line and the other end of the line.
̅→ ̅→ ̅→
VB1 − VA1 → → IA1 →
̅→ = − Z1 + Z1 x − ̅→Z1 x (8)
IB1 IB1
Where, Vf1 is positive sequence component of fault voltage at the
fault point f and Zf is a combination of negative sequence, zero sequence
or fault impedance depending upon the type of fault. Eqs. (7) and (8)
Fig. 10. Discrepant impedances during a fault at F-3(3-Φfault) on feeder 4–5. represent positive sequence impedance of the feeder when estimated
from Bus-A and Bus-B respectively. It is clear that the discrepant
impedance during a fault will not be zero as is the case during a no-fault
condition. Its exact value will depend on the fault location, currents and
voltages during the fault.
The R-X plot for the discrepant impedance estimated by the relay is
shown in Fig. 3. Under no-fault condition, the circle shown is Fig. 3 shall
theoretically have zero radius and lie at the origin. However, when the
errors encountered during measurement of currents and voltages are
accounted, therefore, the circle of no-fault region will have a small
radius. For the studies reported in this paper, the radius of the no-fault
Fig. 11. Discrepant impedances for healthy feeders during a fault at F-3(3- region circle was selected to be 3% of the feeder impedance. In field
Φ fault). applications, a setting close to 3–5% should suffice depending on the
measurement accuracy of instrument transformers and required sensi
→ → → →
V B1 = V A1 − I A1 Z1 (1) tivity of the relay. The relay will issue a trip command to the breaker
after five consecutive values of discrepant impedances are in fault region
Then, as depicted in Fig. 3 trip logic.
→ →
→ V A1 − V B1
Z1 = → (2) 3. Test results
I A1
The IEEE-9 bus test system with the combination of grid forming and
Z1 is the positive-sequence impedance of the feeder as estimated from
grid following inverters is developed in PSCAD. The model developed in
Bus-A using positive-sequence voltages from Bus-A and Bus-B, and
reference [8] with modifications is used as a test bed for this paper. Load
positive-sequence current measured at Bus-A.
data, generation data and line parameters are provided in Table-1 and
4
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Fig. 12. (a) Discrepant impedances for feeder between bus-4 and bus- 5 during a high-resistance a fault at F-1 and (b) discrepant impedances of the healthy feeder.
Fig. 13. (a) Discrepant impedances for feeder between bus-4 and bus- 5 during a high-resistance a fault at F-2 and (b) discrepant impedances of the healthy feeders.
Table-2 of Appendix. Fig. 4 describes the network architecture of ends of the line and discrepant impedances are estimated as discussed in
microgrid with three inverters, feeders and loads. Two inverters are Section 2. The proposed protection philosophy is validated for different
operating as grid-forming inverters and one inverter is operating as cases including different types of faults, fault locations, and type of
grid-following inverter. Each inverter is rated at 200 MVA. Solar and inverter present in the network both in gird-connected and islanded
wind resources were achieved using grid-following inverters whereas a modes.
battery source is connected via grid forming inverter. Loads connected
at buses are constant impedance loads. The lines are modeled using a pi
model. Positive-sequence voltages and currents are estimated at both
5
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Fig. 15. Single line diagram of a medium voltage microgrid modeled for real time testing.
Fig. 16. Discrepant impedances of line 3–7 and line 2–6 during LG fault (Grid
connected mode).
6
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Fig. 19. Differential relay (alpha plane) response for fault F-2 in Fig. 4.
Fig. 20. Current for a L-G fault on line 3–7 during grid-connected mode.
7
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
Fig. 21. Current for a L-G fault on line 3–7 during islanded mode.
Fig. 22. Discrepant Impedances for a L-G fault on line 3–7 during a) Islanded mode and b) grid connected mode.
Table 2
Line parameters of the test system in Fig. 4.
Sequence R (ohm/m) XL((ohm/m) XC(ohm/m)
4
Positive 1.07 × 10− 4.27 × 10− 4 2.5448 × 106
4
Zero 5.35 × 10− 1.153 × 10− 3 4.1642 × 106
Table 3
Line parameters of the test system in Fig. 15.
Sequence R (ohm/m) XL((ohm/m) XC(ohm/m)
− 4 − 4
Positive 1.73 × 10 4.317 × 10 3.626 × 106
Zero 3.5 × 10− 4 1.79 × 10− 3 8.846 × 106
8
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
and B-7 are provided to the relay. The discrepant impedances are exchange bewteen relays, PMUs and merging units can be wired or
computed by the relay. wireless depending on the geographical location of the microgrid.
It can be observed from Fig. 16 that discrepant impedances for line
3–7 before the fault are in the no-fault region and once the fault occurs, 7. Conclusions
the values are in the fault region. The relay computes discrepant im
pedances continuously. Once fault strikes on the feeder 3–7, the This paper presents a new microgrid feeder protection technique that
discrepant impedances are in the fault region and the relay issues a trip uses positive sequence voltages and currents from both ends of the
signal to the breaker after 1/4 cycle and once the relay issues a trip, the feeder. The technique estimates the discrepant impedance of the feeder
breaker clears the fault in 6 cycles as shown in Fig. 17. However, as seen under protection. The proposed protection scheme has been tested using
in Fig. 16, the discrepant impedances of adjacent healthy line 2–6 is not PSCAD simulations as well as via real-time implementation. Results
affected and are in the no-fault region before the fault and after fault demonstrated the effectiveness of the technique in detecting and
clearance. clearing faults in a variety of situations i.e. islanded mode, grid-
connected mode and for various faults including high-resistance faults.
4.1.2. Islanded mode of operation The results prove that the proposed protection scheme is effective irre
It can be observed from Fig. 18 that the discrepant impedances for spective of the type of inverter, microgrid topology, inverter control
line 3–7 computed during the fault are in the fault region, while philosophy, level of inverter-based generation and mode of operation.
discrepant impedance on the healthy feeder 2–6 are in the no-fault re The proposed scheme is fast to detect all types of faults in less than ½
gion before, during and after the fault is cleared. cycle. The technique can be implemented as a centralized, decentralized
protection scheme or as a hybrid protection scheme.
5. Comparison with differential protection
Author credits
In industry, two types of differential protection are commonly used
a) current differential using an alpha plane characteristic [30] and b) Tarlochan Sidhu: Conceptualization, Funding acquisition, Project
percentage differential relay. As demonstrated in this section, both dif administration, Writing – review/editing, Supervision
ferential protection methods do not protect a microgrid in the islanded Jigyesh Sharma: Simulations, Writing – Initial draft, Data curation,
mode of operation. On the other hand, the protection technique pro Formal analysis , Validation
posed in this paper works properly in both islanded and grid connected
modes. To verify this problem associated with line differential relay with
an alpha plane characteristic, a L-L fault at F-2 is simulated on a test Declaration of Competing Interest
microgrid shown in Fig. 4, and the response of the relay is shown in Fig
19. The authors declare that they have no known competing financial
The trajectory of the operating quantity doesn’t go into the trip zone, interests or personal relationships that could have appeared to influence
therefore relay does not operate for this fault. As shown in Fig. 8. the work reported in this paper.
discrepant impedances for the same fault fall in the fault region and the
proposed relay will trip. To demonstrate the problem with the per Data availability
centage differential relay, a LG fault is simulated on line 3–7 of a
microgrid shown in Fig. 15 when DG capacity is 2 MW and load L-6 is The authors are unable or have chosen not to specify which data has
removed. been used.
During grid-connected mode, the fault current is Imax= 2.33 kA and
during islanded mode, it is Imin=0.1 kA which can be observed from
Fig. 20 and Fig. 21 respectively. The setting of the differential relay is Appendix:. Test System Data
chosen to be 0.065 times the maximum fault current [31]. Therefore, the
differential relay is set at, 0.065x2.33kA which is 0.151 kA. The set value Table 1, 2 and 3
is higher than the minimum fault current observed when the microgrid
is in the islanded mode of operation. Therefore, the differential will not References
operate in islanded mode due to the dynamic nature of the fault current.
The discrepant impedances in islanded mode and grid connected modes [1] S. Parhizi, H. Lotfi, A. Khodaei, S. Bahramirad, State of the art in research on
are shown Fig. 22(a) and 22(b) respectively. The discrepant impedances microgrids: a review, IEEE Access 3 (2015) 890–925, https://doi.org/10.1109/
ACCESS.2015.2443119.
in both cases are in the fault region indicating correct operation of the [2] Pan Nianan and Wang Luyang, PSCAD studies demonstrate grid forming inverters
proposed relaying scheme. can improve weak grid of Australia. [Online]. Available:https://en.sungrowpower.
com/upload/file/20211201/Grid%20Forming%20Whitepaper.pdf.
[3] R.W. Kenyon, A. Sajadi, A. Hoke, and B.-M. Hodge, “Open-source PSCAD grid-
6. Practical implementation following and grid-forming inverters and a benchmark for zero-inertia power
system simulations,” in 2021 IEEE Kansas Power and Energy Conference (KPEC),
The proposed protection scheme can be realized either as a Manhattan, KS, USA, Apr. 2021, pp. 1–6. 10.1109/KPEC51835.2021.9446243.
[4] A. Srivastava, R. Mohanty, M.A.F. Ghazvini, L.A. Tuan, D. Steen, and O. Carlson, “A
centralized protection scheme or as a decentralized protection scheme as review on challenges and solutions in microgrid protection,” in 2021 IEEE Madrid
depicted in Fig. 23. In centralized protection scheme, the data from PowerTech, Madrid, Spain, Jun. 2021, pp. 1–6. 10.1109/PowerTech46648.2021
PMUs located at the ends of the line are transferred to a central unit .9495090.
[5] S.S. Rath, G. Panda, P.K. Ray, and A. Mohanty, “A comprehensive review on
where protection algorithm is implemented for all lines in the microgrid. microgrid protection: issues and challenges,” in 2020 3rd International Conference
In small microgrids, PMUs can be replaced with merging units, voltage on Energy, Power and Environment: Towards Clean Energy Technologies, Shillong,
and current samples can be transferred via IEC 61850 protocol to the Meghalaya, India, Mar. 2021, pp. 1–6. 10.1109/ICEPE50861.2021.9404520.
[6] M.W. Altaf, M.T. Arif, S. Saha, S.N. Islam, M.E. Haque, and Amt. Oo, “Renewable
central unit. In decentralized protection scheme the voltage and current
energy integration challenge on power system protection and its mitigation for
samples are exchanged between relays and each relay implements the reliable operation,” in IECON 2020 The 46th Annual Conference of the IEEE
proposed relay algorithm and scheme. A hybrid architecture of Industrial Electronics Society, Singapore, Oct. 2020, pp. 1917–1922. 10.110
centralized and decentralized can also be impleted provided that relays 9/IECON43393.2020.9254380.
[7] Central Power Research Institute, Protection challenges under bulk penetration of
are intellegient with features of phasor measurement and handling IEC renewable energy resources in power systems: a review, CSEE JPES 3 (4) (2017)
61,850–9–2 sampled value communication. The medium used for data 365–379, https://doi.org/10.17775/CSEEJPES.2017.00030. Dec.
9
J. Sharma and T.S. Sidhu Electric Power Systems Research 224 (2023) 109632
[8] J. Shiles et al., “Microgrid protection: an overview of protection strategies in North [20] T. Patel, S. Brahma, J. Hernandez-Alvidrez, and M.J. Reno, “Adaptive protection
American microgrid projects,” in 2017 IEEE Power and Energy Society General scheme for a real-world microgrid with 100% inverter-based resources,” in 2020
Meeting, Chicago, IL, Jul. 2017, pp. 1–5. 10.1109/PESGM.2017.8274519. IEEE Kansas Power and Energy Conference (KPEC), Manhattan, KS, USA, Jul. 2020,
[9] M.N. Alam, S. Chakrabarti, A.K. Pradhan, Protection of networked microgrids pp. 1–6. 10.1109/KPEC47870.2020.9167527.
using relays with multiple setting groups, IEEE Trans. Ind. Inf. 18 (6) (2022) [21] T.E. Sati, M.A. Azzouz, Optimal protection coordination for inverter dominated
3713–3723, https://doi.org/10.1109/TII.2021.3120151. Jun. islanded microgrids considering N-1 contingency, IEEE Trans. Power Delivery 37
[10] T.S. Sidhu and D. Bejmert, “Short-circuit current contribution from large scale PV (3) (2022) 2256–2267, https://doi.org/10.1109/TPWRD.2021.3108760. Jun.
power plant in the context of distribution power system protection performance,” [22] S. Gopalan, V. Sreeram, H. Iu, and Y. Mishra, “An improved protection strategy for
in IET Conference on Renewable Power Generation (RPG 2011), Edinburgh, UK, microgrids,” in IEEE PES ISGT Europe 2013, Lyngby, Denmark, Oct. 2013, pp. 1–5.
2011, pp. 134. doi: 10.1049/cp.2011.0153. 10.1109/ISGTEurope.2013.6695477.
[11] T. Neumann, I. Erlich, Short circuit current contribution of a photovoltaic power [23] M.A. Haj-ahmed, M.S. Illindala, The influence of inverter-based DGs and their
plant, IFAC Proceed. Volumes 45 (21) (2012) 343–348, https://doi.org/10.3182/ controllers on distribution network protection, IEEE Trans. Ind. Applicat. 50 (4)
20120902-4-FR-2032.00061. (2014) 2928–2937, https://doi.org/10.1109/TIA.2013.2297452. Jul.
[12] J. Li, X. Zhang, W. Yao, J. Xu, X. Shi, and F. Ren, “Inverse time differential current [24] C.A. Plet, M. Graovac, T.C. Green, and R. Iravani, “Fault response of grid-connected
protection for microgrid,” in 2022 IEEE 5th International Electrical and Energy inverter dominated networks,” in IEEE PES General Meeting, Minneapolis, MN,
Conference (CIEEC), Nangjing, China, May 2022, pp. 590–595. 10.1109/CIEE Jul. 2010, pp. 1–8. 10.1109/PES.2010.5589981.the.
C54735.2022.9845849. [25] M.J. Reno, S. Brahma, A. Bidram, M.E. Ropp, Influence of inverter-based resources
[13] Power System Relaying and Control CommitteeWorking Group C30, “Microgrid on microgrid protection: part 1: microgrids in radial distribution systems, IEEE
protection system,” IEEE, Technical report PES-TR71, Jul. 2019. [Online]. Power and Energy Mag. 19 (3) (2021) 36–46, https://doi.org/10.1109/
Available: https://resourcecenter.ieee-pes.org/publications/technical-reports/ MPE.2021.3057951. May.
PES_TP_TR71_PSRC_microgrid_082019.html. [26] A. Hooshyar, R. Iravani, A new directional element for microgrid protection, IEEE
[14] Z. Shuai, C. Shen, X. Yin, X. Liu, Z.J. Shen, Fault analysis of inverter-interfaced Trans. Smart Grid 9 (6) (2018) 6862–6876, https://doi.org/10.1109/
distributed generators with different control schemes, IEEE Trans. Power Delivery TSG.2017.2727400. Nov.
33 (3) (2018) 1223–1235, https://doi.org/10.1109/TPWRD.2017.2717388. Jun. [27] Behnam Mahamedi, John Edward Fletcher, Trends in the protection of inverter-
[15] E. Sortomme, J. Ren, and S.S. Venkata, “A differential zone protection scheme for based microgrids, IET Gener. Transm. Distrib. 13 (20) (2019) 4511–4522. Iss.
microgrids,” in 2013 IEEE Power and Energy Society General Meeting, Vancouver, [28] A protection method for inverter-based microgrid using current-only polarity
BC, 2013, pp. 1–5. 10.1109/PESMG.2013.6672113. comparison, J. Modern Power Syst. Clean Energy 8 (3) (2020) 446–453, https://
[16] T.S. Ustun, C. Ozansoy, and A. Zayegh, “Differential protection of microgrids with doi.org/10.35833/MPCE.2018.000722. May.
central protection unit support,” in IEEE 2013 Tencon - Spring, Sydney, Australia, [29] O. Nzimako and A. Rajapakse, "Real time simulation of a microgrid with multiple
Apr. 2013, pp. 15–19. 10.1109/TENCONSpring.2013.6584408. distributed energy resources," 2016 International Conference on Cogeneration,
[17] E. Casagrande, W.L. Woon, H.H. Zeineldin, D. Svetinovic, A differential sequence Small Power Plants and District Energy (ICUE), Bangkok, Thailand, 2016, pp. 1–6,
component protection scheme for microgrids with inverter-based distributed doi: 10.1109/COGEN.2016.7728945.
generators, IEEE Trans. Smart Grid 5 (1) (2014) 29–37, https://doi.org/10.1109/ [30] G. Benmouyal, The trajectories of line current differential faults in the alpha plane.
TSG.2013.2251017. Jan. Line Current Differential Protection: A Collection of Technical Papers Representing
[18] R. Lasseter and H. Nikkhajoei, “Microgrid protection,” p. 6. Modern Solutions, Schweitzer Engineering Laboratories, Inc, 2014.
[19] B. Liao, J. Cheng, and G. Ren, “Microgrid adaptive current instantaneous trip [31] Shi Chen, et al., Sequence-component-based current differential protection for
protection,” in 2019 IEEE Innovative Smart Grid Technologies - Asia (ISGT Asia), transmission lines connected with IIGs, IET Gener. Transm. Distrib. 12 (12) (2018)
Chengdu, China, May 2019, pp. 2074–2078. 10.1109/ISGT-Asia.2019.8881166. 3086–3096.
10