Linear Programming Coordination for Overcurrent Relay in Electrical Distribution Systems with Distributed Generation

Abstract

Electric power distribution networks are generally radial in nature, with unidirectional power flows transmitted from the highest voltage levels to the consumption levels. The protection system in these distribution networks is relatively simple and consists mainly of fuses, reclosers (RC) and overcurrent relays (OCRs). The installation of distributed generation (DG) in a network causes coordination problems between these devices, because the power flows are no longer unidirectional and can flow upstream to the substation. For this reason, the work proposed here analyzes the most significant impacts that DG has on the protection devices and proposes an adjustment method for the OCRs based on linear programming (LP) techniques with the aim of improving their response time to the different faults that may occur in the main feeder of the network. The distribution system selected for the study is the IEEE 34 bus system using DIgSILENT 14.1 software for its modeling and Matlab for the adjustment of the overcurrent devices. Results indicate that better coordination between protection devices are achieved if LP is used.

1. Introduction

The current philosophy in the planning, management and control of a radial power distribution network is based on the assumption of the existence of unidirectional power flows, which is transmitted from the highest transport voltage levels to the distribution levels. We can assume that short-circuit currents behave similarly. These assumptions allow for the implementation of relatively simple and inexpensive protection schemes with which a selective operation of the protection system is achieved [1,2]. According to the principles of selectivity [3], only the protection device closest to the fault should operate to clear the fault, leaving the rest of the network energized.

The installation of DGs at medium and low voltage levels changes this fundamental basis. Power flows and short-circuit currents can now have different upstream directions and values [4]. As a consequence, the initial schemes implemented (e.g., main feeder protection) may no longer work or be less effective [1,5].

With the presence of DG the distribution system becomes more active, i.e., both load and generation significantly affect the state of the network [6]. The effects of DG on distribution network devices have been discussed in different literature [4,7]. To reduce them, several “mitigation methods” have been proposed, thus, several authors recommend acting directly on the DG [3,8], e.g., disconnecting it just before fault detection or limiting its power and therefore its effect on the protection system [9], others recommend reconfiguring the network topology [10,11], or the use of fault current limiters (FLCs) [12,13].

In recent years, mitigation methods based on coordination algorithms have been developed using linear programming (LP), non-linear programming (NPL), or genetic algorithms (GA) [14,15,16,17] that either modify the coordination interval time (CTI) between devices, determine the optimal location of the DG [15], or in combination, without compromising the protection system.

Although these methods can be effective, they have some disadvantages and limitations. Therefore, the disconnection of the DG immediately before fault detection proposed in [3,18], can cause asynchronous reconnections and cause severe damage to both the DG and the distribution network, especially at high penetration levels [3,10]. The limitation of the DG capacity proposed in [9] is also undesirable, since it also limits the penetration level and therefore the advantages of its installation close to the consumption points, (loss reduction, improvement of supply quality, etc.). The network reconfiguration suggested in [10,11] is costly and sometimes unfeasible. Finally, the use of current limiters [12,14] entails an additional installation cost that can be prohibitive.

Regarding the use of NLP methods for the determination of the coordination index mentioned in [15] can be very complex, particularly if the number of protection devices and DG’s is high, since the calculation time required to determine their operating parameters will also be high [19,20,21], on the other hand, not always the optimal location point of DG is viability from the perspective of network operation, e.g., availability of power evacuation at the point of common connection (PCC).

The investigations described in [14,16,17], use GA-based optimization algorithms to determine optimal coordination between OCR and distance relays [14]. In [16] a hybrid GA and NLP algorithm is developed to choose the optimal value of DIAL or time multiplier setting (TMS) of all OCR, the use of a continuous genetic algorithm (GAC) faster than GA for the same purpose is developed in [17]. However, all the optimization algorithms described do not consider DG.

In the present work:

• the most significant impacts that DG causes on protection devices are analyzed,
• an adjustment method for OCRs based on LP techniques is proposed,
• the analysis is carried out by using the IEEE 34-node test feeder system modeled in DIgSILENT PowerFactory.

As far as the authors know, it seems that there is a lack of similar studies in the available literature.

2. Effect Caused by the DG on Protective Devices

2.1. Loss of Sensitivity

The problem related to the loss of sensitivity to main feeder protection is often referred to as protection blinding [8]. Installation of DGs in the distribution system can reduce the value of the short-circuit current detected by the protection of the main substation and affect the response time of the circuit breaker, which will depend, to a large extent, on the size and location of the DG within the distribution network.

The loss of sensitivity of the protection device can be analyzed as a function of where the short circuit is located in relation to the location of the DG, distinguishing between short circuits located upstream of the DG and short circuits located downstream of it. Figure 1 shows the topology of a radial distribution network where S.E is the substation, R is the main feeder protection, RC is the automatic recloser, and F is a protection fuse for the branch line.

For faults located upstream (between the substation—S.E, and the DG), the contribution current of the S.E to the short-circuit is independent of the size of the DG. A three-phase short-circuit is also represented upstream of the DG, for different penetration levels: 17%, 33%, and 50%. In the operating characteristic curves of the time overcurrent (t/i) devices shown in Figure 2, this effect is observed, as the current detected by the relay remains constant at 1.59 kA, regardless of the of penetration level of the DG.

In contrast, for faults located downstream of the substation and the DG (as shown in Figure 3), it is observed that the current contribution of the DG to the short circuit causes the current measured by the main relay, R, of the S.E to decrease from 1.21 with 17% penetration to 1.14 kA with 50% penetration level of the DG, which represents a 6.1% loss in the sensitivity of the protection device. This effect will cause a delay in the operating time of the protection device, as shown in Figure 4. The extreme case of this effect occurs when the main relay does not detect the fault current, as shown in Figure 5. The maximum value at which the relay stops detecting the fault current appears when the DG is installed close to the substation and under fault conditions at the main feeder [12].

2.2. Loss of Coordination

Under normal operating conditions, protection devices are coordinated in such a way that, in the event of a fault in the network, the main protection acts before backup protection. A large majority of faults occurring in the network are of a temporary nature, the purpose of the recloser, RC, is to try to clear these faults and on the one hand avoid an unnecessary interruption of the electrical service, and on the other hand to safeguard the protection fuse of the branch line under its supervision. Depending on the initial settings of the coordination schemes, and taking into account the size, location, and type of DG installed, it may happen that a temporary fault in a branch line causes the loss of coordination between fuse and recloser, affecting their coordination time due to the additional current with which the DG contributes to the short-circuit. To observe this impact, the situation of temporary fault, I_CC, in a branch line (lateral) is shown in Figure 6.

The characteristic t/i trip curves are represented in Figure 7 where the impact called (fuse nuisance blowing) is shown. We can observe how a temporary fault in a branch line causes a loss of coordination between the fuse and the recloser. Thus, the fuse detects a fault current of 0.25 kA and trips in 0.15 s, while the recloser is crossed by a current of 0.650 kA, which would cause it to trip in 0.32 s. This time is higher than the fuse tripping.

3. Coordination with Linear Programming

OCRs are currently the most widely used protection devices [13,22], this device generally being used as backup protection. However, in some situations, it may be the only protection. The OCR should act as primary protection by removing faults that occur in the area under its supervision. Only in the event that the primary protection fails to operate should the backup protection initiate tripping. A typical electrical distribution network may consist of hundreds of equipment protection relays. Each relay in the system must be coordinated with the adjacent equipment protection relay. If backup protection is not well coordinated, coordination failure may occur, and therefore, coordination of backup protections is one of the main concerns in the network protection system [13,23]. Generally, protection coordination can be performed by topology [24,25], by optimization methods [26,27], or by expert methods [28]. Topological analysis is used for relay tuning in multiterminal networks, graph theory, and functional approximation techniques are employed to provide the best solution that does not necessarily have to be the optimal one. In optimization methods, some researchers [26,29] use non-linear programming techniques to determine the optimal relay settings, subject to constraints due to coordination and limits of the relay settings themselves. In reference [30], the big M method is proposed to find the optimal value of the OCR time multiplier setting (TMS) in which the stated values of the plug setting are assumed to be known and fixed.

The problem of OCR coordination in the distribution system with DG presence can thus be defined as an optimization problem with constraints. The objective is to minimize the operating time of the relay closest to the place where the fault has occurred. The constraints imposed are due to limits on the operating time of the relay, coordination criteria, and relay characteristics. In this work, the OCR coordination problem in radial distribution systems is formulated as a linear programming (LP) problem with constraints, where the sum of the operating times of the relays in the system for different minimum fault points:

$$\min \frac{\alpha }{{\left(k_{i,j}\right)}^{\beta }-1} {{\displaystyle \sum }}_{i=1}^n\left({{\displaystyle \sum }}_{j=1}^m{\mathrm{TMS}}_{i,j}-{\mathrm{TMS}}_{i+1,j+1}\right),$$

(1)

where α and β are OCR shape constants according to

Table 1. Form constants for the exponential equation IEC [31].

Type of Curve	α	β
Standard inverse time	0.14	0.02
Very inverse	13.5	1
Extremely inverse	80	2
Long inverse time	120	1

, k_i,j is the ratio between the short-circuit current, I_CC, and the relay setting current, I_ar, regarding relay i with fault in section j, and TMS is the time multiplier setting.

The proposed optimization problem is subjected to the following group of constraints:

• Coordination criterion—The protection coordination criterion establishes the minimum time that must elapse between the operation of the primary protection and the operation of the backup protection. The fault is detected simultaneously by the primary protection, PP, and the secondary protection, PS. To avoid erroneous operation, the PS will only have to operate in case the PP fails. If we define R_i as the primary fault protection at a certain point, j, and R_i+1 as the secondary or backup protection for the same fault. The constraint condition for coordination criteria is:

  $$t_{i+1,j}-t_{i,j}>\mathsf{\Delta }t,$$

  (2)

  where, t_i+1,j is the operation time of the back-up, PS, operation for fault at point j, t_i,j is the operation time of PP operation for the same fault, and ∆t is the coordination time interval (CTI).
• Relay operating time limits—The time taken by a relay to detect and isolate a fault produced in its zone of influence must be bounded. This is the constraint imposed by the operating time of the relays:

  $$t_{i,min}\le t_{i,j}\le t_{i,max}$$

  (3)

  where, t_i,min and t_i,max are the minimum and maximum operating times of relay i (at any point).
• Time multiplier setting, TMS, limits—The operating time of a relay is directly proportional to the TMS. Therefore:

  $${\mathrm{TMS}}_{i,\min }\le {\mathrm{TMS}}_{i,j}\le {\mathrm{TMS}}_{i,\max },$$

  (4)

  where TMS_i,min is the minimum TMS value for relay i and, TMS_i,max is maximum value of TMS for relay i. The TMS values usually taken are 0.025–1.2, respectively [32].
• Relay operating characteristics—To extend the work region of the overcurrent relay specified by the standard [33], the relay’s parameters are optimized considering the maximum value of the relay and the minimum value of the relay, as expressed respectively. We consider all relays with the same characteristics:

  $$t_{op}=\frac{\alpha }{{\mathrm{PSM}}^{\beta }-1}\mathrm{TMS},$$

  (5)

  where t_op is the relay operating time and PSM is the plug setting multiplier (k_i,j in Equation (1)). For inverse characteristic curves it is usually α = 0.02 and β = 0.14 [32,33].

The relay setting current, I_ar, is determined from the system requirements. Thus, the above equation can be rewritten as:

$$t_{op}=a\times \mathrm{TMS},$$

(6)

where:

$$a=\frac{\alpha }{{\mathrm{PMS}}^{\beta }-1}.$$

(7)

4. Case Study Results and Discussion

The distribution system selected for the simulations is the one taken from the test feeders of the Distribution System Analysis Subcommittee of the Institute of Electrical and Electronics Engineers (IEEE) [34], shown in Figure 8. The main characteristics of the distribution network on which the short circuits were simulated are the following:

• The voltage across the distribution network is 15 kV.
• The short-circuit power of the external network will be 100 MVA.
• The power demand of the electrical loads connected to the network is 12 MW and load flow analysis they have been considered balanced.
• Overhead network with double circuit with LA-110 conductor. The electrical characteristics of this conductor are shown in [35].

Three-phase and single-phase short circuits were simulated on the distribution network, at busbars 2, 7, and 21, as shown in Figure 8, and with different DG penetration levels ranging from 33% to 50%. A first distributed generation unit, DG1, was installed at busbar 20, whereas a second one, DG2, was installed at busbar 22.

The DG protection device should be properly coordinated with the other protection devices in the distribution network. When a fault occurs on the line to which the DG unit is connected, the DG unit must be disconnected before the protection device on that line is disconnected. This ensures that the DG does not cause disturbances in the operation of the line protections. Furthermore, when a fault occurs on a branch line, the DG must remain in service and will not disconnect before the fault line does [8]. According to this, the coordination criterion shown in

Table 2. Coordination criteria according fault location.

Bus (Node)	Location Faults	Main Protection	Backup Protection
2	F2	OCR2	OCR1
7	F7	RC	OCR1
21	F21	RC	OCR3

is proposed.

The transformation ratio of the current transformers used in the test is 400/5, chosen according to the full-load currents of the system and the short-circuit currents detected by each of the protection devices according to the location of the fault (shown in

Table 3. Short circuit (expressed in amperes, A) detected by the protection devices, for different grounding resistor values, R_f, expressed in ohm.

			F2			F7			F21
			Penetration Level DG			Penetration Level DG			Penetration Level DG
Fault Type			17%	33%	50%	17%	33%	50%	17%	33%	50%
Three-phase	IOCR1		1568	1568	1568	1193	1116	1104	321	325	328
	IOCR2(DG)		625	894	1077	570	778	919	392	407	417
	IOCR3		-	-	-	-	-	-	361	368	373
	IRC		-	-	-	1717	1922	2052	667	681	689
Single-phase	Rf = 0	IOCR1	539	608	665	505	556	575	404	397	395
		IOCR2(DG)	710	1076	1351	670	960	1165	480	465	474
		IOCR3	30.8	30	29	31	30	29	230	238	241
		IRC	644	627	619	1078	1445	1662	728	746	754
	Rf = 5	IOCR1	564	604	623	517	533	532	398	195	393
		IOCR2(DG)	736	943	1059	669	827	918	446	461	469
		IOCR3	32.45	33	33	31.6	32	32	211	216	219
		IRC	677	697	706	1060	1263	1359	733	748	754
	Rf = 10	IOCR1	523	535	535	490	488	462	396	393	464
		IOCR2(DG)	664	677	820	616	703	751	443	456	391
		IOCR3	33.3	33	33.9	32	33	33	194	199	200
		IRC	695	705	708	1000	1103	1146	735	747	753
	Rf = 15	IOCR1	488	491	488	466	459	452	394	391	389
		IOCR2	605	667	699	572	628	658	440	452	459
		IOCR3	33	33	33.6	32	33	33	180	184	195
		IRC	698	703	699	944	1004	1028	735	746	751
	Rf = 20	IOCR1	464	464	461	448	440	434	392	448	387
		IOCR2(DG)	563	607	629	540	580	602	437	389	454
		IOCR3	33	33.5	33	33	33	33	168	171	172
		IRC	697	699	699	900	040	956	734	744	748
	Rf = 25	IOCR1	447	446	443	435	427	422	390	387	385
		IOCR2(DG)	534	567	598	517	548	565	434	444	456
		IOCR3	33	33.3	33.3	33	33	33	158	160	161
		IRC	696	697	696	869	897	908	733	742	745

).

4.1. Objective Function According to Fault Location

According to the location of the fault, and based on the protection criteria included in Table 2, three objective functions to minimize can be obtained, each of them will result in a different adjustment, an adaptive adjustment.

The scenario analyzed is the three-phase short-circuit scenario with a DG1 penetration level of 17%, the results obtained for the other scenarios can be found in

Table A3. Operation time PL coordination versus operation time from classical coordination.

			F2						F7						F21
Type of Fault			top 17%		top 33%		top 50%		top 17%		top 33%		top 50%		top 17%		top 33%		top 50%
	Type Coordination		C	PL	C	PL	C	PL	C	PL	C	PL	C	PL	C	PL	C	PL	C	PL
Three-phase	IOCR1		0.42	0.39	0.42	0.39	0.42	0.39	-	-	-	-	-	-	-	-	-	-	-	-
	IOCR2(GD)		0.20	0.20	0.19	0.10	0.15	0.08	0.20	0.13	0.20	0.10	0.20	0.10	-	-		-	-	-
	IOCR3		-	-	-	-	-	-	-	-	-	-	-	-	0.42	0.38	0.42	0.40	0.42	0.40
	IRC		-	-	-	-	-	-	0.45	0.42	0.45	0.38	0.45	0.40	0.12	0.08	0.12	0.10	0.12	0.10
Single-phase	Rf = 0	IOCR1	0.42	0.39	0.43	0.37	0.42	0.39	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.13	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	0.40	0.39	0.14	0.38	0.14	0.38
		IRC	-	-	-	-	-	-	0.40	0.39	0.40	0.33	0.40	0.31	0.40	0.10	0.40	0.09	0.40	0.09
	Rf = 5	IOCR1	0.41	0.39	0.42	0.37	0.41	0.37	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.10	0.14	0.08	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	0.39	0.43	0.45	0.49	0.45	0.48
		IRC	-	-	-	-	-	-	0.40	0.39	0.40	0.3	0.40	0.40	0.40	0.09	0.40	0.09	0.40	0.09
	Rf = 10	IOCR1	0.41	0.37	0.42	0.40	0.41	0.40	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.11	0.14	0.10	0.14	0.09	0.14	0.09	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	-	0.39	-	0.46	-	0.45
		IRC	-	-	-	-	-	-	0.40	0.38	0.40	0.3	0.40	0.3	0.40	0.10	0.4	0.09	0.40	0.09
	Rf = 15	IOCR1	0.40	0.38	0.41	0.38	0.40	0.38	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.09	0.14	0.10	0.14	0.10	0.14	0.10	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	0.45	0.39	-	0.42	-	0.39
		IRC	-	-	-	-	-	-	0.40	0.39	0.40	0.3	0.40	0.3	0.40	0.09	-	0.09	-	0.09
	Rf = 20	IOCR1	0.40	0.39	0.41	0.39	0.40	0.39	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.11	0.14	0.10	0.14	0.10	0.14	0.10	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	-	0.39	-	0.46	-	0.45
		IRC	-	-	-	-	-	-	0.40	0.39	0.40	0.3	0.40	0.03	0.40	0.09	0.40	0.09	0.40	0.09
	Rf = 25	IOCR1	0.42	0.40	0.40	0.40	0.40	0.40	-	-	-	-	-	-	-	-	-	-	-	-
		IOCR2(DG)	0.14	0.10	0.14	0.10	0.14	0.11	0.14	0.10	0.14	0.11	0.14	0.11	-	-	-	-	-	-
		IOCR3	-	-	-	-	-	-	-	-	-	-	-	-	-	0.40	-	0.50	-	0.40
		IRC	-	-	-	-	-	-	0.40	0.37	0.40	0.3	0.40	0.03	0.40	0.09	0.40	0.09	0.40	0.09

C—classical coordination; PL—lineal programming coordination.

in the Appendix A.

Fault 2: OCR₂ with OCR₁

$$\min z={\displaystyle \sum _{i=1}^m{a}_{i,k}{(\mathrm{TMS})}_i=2.27 {\mathrm{TMS}}_1+3.33 {\mathrm{TMS}}_2}$$

(8)

Fault 7: Recloser with OCR₂

$$\min z={\displaystyle \sum _{i=1}^m{a}_{i,k}{(\mathrm{TMS})}_i=3.49 {\mathrm{TMS}}_2+2.21 {\mathrm{TMS}}_{RC2}}$$

(9)

Fault 21: OCR₃ with recloser

$$\min z={\displaystyle \sum _{i=1}^m{a}_{i,k}{(\mathrm{TMS})}_i=4.57 {\mathrm{TMS}}_3+3.23 {\mathrm{TMS}}_{RC}}$$

(10)

Constraints:

• The minimum operating time of each relay is 0.1 s.
• The normal range of TMS from 0.025 to 1.2.
• Typical CTI setting is 0.3 s.

Coordination criterion

Fault 2: OCR₂ with OCR₁:

2.27TMS₁ − 3.33TMS₂ ≥ 0.3

(11)

Fault 7: Recloser with OCR₂:

3.49TMS₂ − 2.21TMS_RC ≥ 0.3

(12)

Fault 21: OCR₃ with Recloser:

4.57TMS₃ − 3.23TMS_RC ≥ 0.3

(13)

Relay operating time limits

Fault 2: OCR₂ with OCR₁:

2.27TMS₁ ≥ 0.1 TMS₁ ≤ 0.044
   3.33TMS₂ ≥ 0.1 TMS₂ ≤ 0.030

(14)

Fault 7: Recloser with OCR₂:

3.49TMS₂ ≥ 0.1 TMS₂ ≤ 0.038.
  2.21TMS_RC ≥ 0.1 TMS_RC ≤ 0.045

(15)

Fault 21: OCR₃ with recloser:

4.57TMS₃ ≥ 0.1 TMS₃ ≤ 0.021.
   3.23TMS_RC ≥ 0.1 TMS_RC ≤ 0.030

(16)

The normal range of TMS is from 0.025 to 1.2:

TMS₁ ≥ 0.025
   TMS₂ ≥ 0.025
   TMS₃ ≥ 0.025
   TMS_RC ≥ 0.025

(17)

4.2. Relay Operating Time Calculation

Using linear programming, the TMS value was calculated as a function of the fault location point.

Substituting the TMS value and the obtained value of the relay constant in Equation (4), we obtain the operating time of each protection relay.

Fault 2: OCR₂ with OCR₁:

OCR1 t_op1 = a₁(TMS₁) = 2.27·0.17 = 0.39 s
   OCR2 t_op2 = a₂(TMS₂) = 3.33·0.030 = 0.10 s

(18)

Fault 7: Recloser with OCR₂:

OCR2 t_op2 = a₂(TMS₂) = 3.49·0.038 = 0.13 s
   RECLOSER t_opRC = a_RC(TMS_RC) = 2.21·0.11 = 0.24 s

(19)

Fault 21: OCR₃ with Recloser

OCR3 t_op3 = a₃(TMS₃) = 4.57·0.065 = 0.30 s
   RECLOSER t_opRC = a_RC(TMS_RC) = 3.23·0.025 = 0.1 s

(20)

The values of constants a and TMS and relay operation times top for the rest of the scenarios are shown in Appendix A 

Table A1. Values of constant relay (a_i,j) for the different scenarios analyzed.

			F2			F7			F21
Type of Fault			a17%	a33%	a50%	a17%	a33%	a50%	a17%	a33%	a50%
Three-phase	IOCR1		2.27	2.27	2.27	2.52	2.58	2.59	4.96	4.92	4.89
	IOCR2(DG)		3.33	2.83	2.62	3.49	3.00	2.79	4.33	4.23	4.17
	IOCR3		-	-	-	-	-	-	4.57	4.51	4.47
	IRC		-	-	-	2.21	2.13	2.08	3.23	3.19	3.18
Single-phase	Rf = 0	IOCR1	3.59	3.38	3.23	3.72	3.54	3.47	4.25	4.30	4.31
		IOCR2(DG)	3.13	2.62	2.4	3.22	2.74	2.54	3.83	3.90	3.86
		IOCR3		-	-	-	-	-	6.55	6.35	6.27
		IRC	3.28	3.33	3.35	2.62	2.34	2.23	3.10	3.06	3.05
	Rf = 5	IOCR1	3.51	3.39	3.34	3.68	3.63	3.62	4.29	7.78	4.32
		IOCR2(DG)	3.08	2.76	2.64	3.22	2.92	2.79	4.00	3.92	3.88
		IOCR3	-	-	-	-	-	-	7.14	6.97	6.88
		IRC	3.2	3.16	3.14	2.63	2.46	2.40	3.09	3.06	3.05
	Rf = 10	IOCR1	3.65	3.61	3.61	3.79	3.8	3.92	4.30	4.32	3.91
		IOCR2(DG)	3.23	3.2	2.93	3.35	3.15	3.05	4.02	3.95	4.34
		IOCR3	-	-	-	-	-	-	7.83	7.61	7.56
		IRC	3.16	3.14	3.14	2.7	2.59	2.56	3.08	3.06	3.05
	Rf = 15	IOCR1	3.8	3.78	3.8	3.9	3.93	3.97	4.32	4.34	4.35
		IOCR2(DG)	3.39	3.23	3.15	3.48	3.32	3.25	4.03	3.97	3.93
		IOCR3	-	-	-	-	-	-	8.56	8.33	7.78
		IRC	3.16	3.15	3.15	2.76	2.69	2.67	3.08	3.06	3.05
	Rf = 20	IOCR1	3.91	3.91	3.92	3.99	4.03	4.06	4.33	3.99	4.37
		IOCR2(DG)	3.51	3.38	3.32	3.59	3.46	3.39	4.05	4.35	3.96
		IOCR3	-	-	-	-	-	-	9.36	9.14	9.07
		IRC	3.16	3.15	3.15	2.82	2.77	2.75	3.08	3.06	3.06
	Rf = 25	IOCR1	3.99	4.00	4.02	4.06	4.11	4.13	4.34	4.37	4.38
		IOCR2(DG)	3.61	3.5	3.41	3.68	3.56	3.51	4.06	4.01	3.95
		IOCR3	-	-	-	-	-	-	10.2	10.0	9.93
		IRC	3.16	3.16	3.16	2.86	2.82	2.81	3.09	3.07	3.06

and

Table A2. Relay TMS values for the different scenarios analyzed.

			F2			F7			F21
Type of Fault			TMS17%	TMS33%	TMS50%	TMS17%	TMS33%	TMS50%	TMS17%	TMS33%	TMS50%
Three-phase	IOCR1		0.17	0.17	0.17	-	-	-	-	-.	-
	IOCR2(DG)		0.03	0.04	0.04	0.04	0.03	0.04	-	-	-
	IOCR3		-	-	-	-	-	-	0.08	0.09	0.09
	IRC		-	-	-	0.19	0.18	0.19	0.03	0.03	0.03
Single-phase	Rf = 0	IOCR1	0.11	0.11	0.12	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.04	0.04	0.03	0.04	-	-	-	-
		IOCR3	-	-	-	-	-	0.06	0.06	0.06	0.06
		IRC	-	-	-	0.15	0.14	0.03	0.03	0.03	0.03
	Rf = 5	IOCR1	0.11	0.11	0.11	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.04	0.04	0.03	0.03	-	-	-	-
		IOCR3	-	-	-	-		0.06	0.06	0.07	0.07
		IRC	-	-	-	0.15	0.01	0.03	0.03	0.03	0.03
	Rf = 10	IOCR1	0.10	0.11	0.11	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.03	0.04	0.03	0.03	-	-	-	-
		IOCR3	-	-	-	-	-	0.05	0.05	0.06	0.06
		IRC	-	-	-	0.14	0.01	0.03	0.03	0.03	0.03
	Rf = 15	IOCR1	0.10	0.10	0.10	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.03	0.03	0.03	-	-	-	-	-
		IOCR3	-	-	-	-	-	0.05	0.05	0.05	0.05
		IRC	-	-	-	0.14		0.03	0.03	0.03	0.03
	Rf = 20	IOCR1	0.10	0.10	0.10	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.03	0.03	0.03	-	-	-	-	-
		IOCR3	-	-	-	-	-	0.04	0.04	0.05	0.05
		IRC	-	-	-	0.14		0.03	0.03	0.03	0.03
	Rf = 25	IOCR1	0.10	0.10	0.10	-	-	-	-	-	-
		IOCR2(DG)	0.03	0.03	0.03	0.03	-	-	-	-	-
		IOCR3	-	-	-	-	-	0.04	0.04	0.05	0.04
		IRC	-	-	-	0.13	-	0.03	0.03	0.03	0.03

. Table A3, shows a comparison of the operation times between the coordination performed in a classical way and with the calculation of adjustments carried out by using LP. Implementing the results obtained in the simulation software, we obtain the following protection schemes shown below.

The classical coordination plots versus linear programming coordination in the case of three-phase short-circuits in the locations described with 17% DG1, are shown in Figure 9, Figure 10 and Figure 11. It can be seen that the response of the OCR protection devices is optimized using the method proposed in this work.

5. Conclusions

In this work, a unified protection system that offers double functionality is proposed. The system was proven to be capable of both optimizing the settings of the relays for each network situation.

The performance of the relay has been tested for different DG penetration levels, different types of short-circuit faults (single-phase to ground and three-phase), and different fault zones (faults in the main feeder and the laterals).

The proposed scheme can be implemented in distribution systems using the communication standard IEC61850 [31].

In all the cases studied, the system has proven capable of adapting, in real time, the performance parameters of the relays that protect the faulted line section depending on the operating conditions of the system at the time of the fault.

The proposed protection system has been applied to a distribution system in the presence of DG. The results demonstrate that the protection coordination proposed in this paper is capable of reducing relay operating times by more than 80% compared to classical coordination, which produces a more reliable operation of a network with DG.