A+ A A-

I. Introduction:

In this test system, both near-end and far-end 3ϕ fault locations are considered, as shown in the figure given below.
This network consists of 4 buses, 4 branches, and 8 directional overcurrent relays (DOCRs).
The objective function of this test system is mathematically expressed as follows:
min Z=mp=1Tnearpr,p+nq=1Tfarpr,q ...... (1)


where Tnearpr,p and Tfarpr,q are respectively the operating time of the primary relay at the near-end 3ϕ fault (at the pth location) and the far-end 3ϕ fault (at the qth location), which can be calculated as follows:
Tnearpr,p=TMSp[β(apPSp×bp)α1] ...... (2)

Tfarpr,q=TMSq[β(cpPSq×dp)α1] ...... (3)

where:
α and β are the time-current characteristic curve (TCCC) constants. For IDMT-based DOCRs, α and β are equal to 0.14 and 0.02, respectively. 
Note that the constants of Tfarpr in Eq.(3) are calculated in sequence of p instead of q.
Also, these p and q notations are i and j in the original references. This replacement is necessary to avoid any confusion with other systems because the notations i and j are assigned for primary and backup relays, respectively.
The number of selectivity constraints depends on the network where only 9 selectivity constraints out of 16 are considered. This results in having a total number of 33 constraints.
The selectivity constraint is mathematically expressed as follows:
TjkTikCTI ...... (4)

where Tjk and Tik are respectively the operating times of the jth backup and ith primary relays for a 3ϕ fault occurs at the kth location. They can be computed by the following equations:
Tjk=TMSj[β(epPSj×fp)α1] ...... (5)

Tik=TMSi[β(gpPSi×hp)α1] ...... (6)

For this test system, the primary operating time of each relay (Ti) should be bounded between the following lower and upper limits:
Tmin(=0.05)TTmax(=1.0) ...... (7)

Both plug-setting (PS) and time-multiplier setting (TMS) are considered as continuous independent variables. The CT rations (CTRs) for the relays R1 to R8, the listed primary/backup (P/B) relay pairs, and the corresponding 3ϕ fault currents available as {a,b,c,d,e,f,g,h} constants all are given below (click on them for bigger size):

II. Single-Line Diagram:

This single-line diagram was drawn by Ali R. Al-Roomi in Sept. 2013 and all the necessary data were coded in MATLAB m-files.

III. Files:

System DATA (MATLAB, m-file Format) [Download]
Results Tester (MATLAB, m-file Format) [Download]

IV. Citation Policy:

If you publish material based on databases obtained from this repository, then, in your acknowledgments, please note the assistance you received by using this repository. This will help others to obtain the same data sets and replicate your experiments. We suggest the following pseudo-APA reference format for referring to this repository:

Ali R. Al-Roomi (2015). Optimal Relay Coordination Test Systems Repository [https://www.al-roomi.org/coordination]. Halifax, Nova Scotia, Canada: Dalhousie University, Electrical and Computer Engineering.

Here is a BiBTeX citation as well:

@MISC{Al-Roomi2015,
author = {Ali R. Al-Roomi},
title = {{Optimal Relay Coordination Test Systems Repository}},
year = {2015},
address = {Halifax, Nova Scotia, Canada},
institution = {Dalhousie University, Electrical and Computer Engineering},
url = {https://www.al-roomi.org/coordination}
}

V. References:

[1] Ali R. Al-Roomi. Optimal Coordination of Power Protective Devices with Illustrative Examples. Hoboken, New Jersey: Wiley-IEEE Press, 2022.