# Modeling Electromechanical Overcurrent Relays Using Singular Value Decomposition

Authors: Feng-Jih Wu; Chih-Ju Chou; Ying Lu; Jarm-Long Chung
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

This paper presents a practical and effective novel approach to curve fit electromechanical (EM) overcurrent (OC) relay characteristics. Based on singular value decomposition (SVD), the curves are fitted with equation in state space under modal coordinates. The relationships between transfer function and Markov parameters are adopted in this research to represent the characteristic curves of EM OC relays. This study applies the proposed method to two EM OC relays: the GE IAC51 relay with moderately inverse-time characteristic and the ABB CO-8 relay with inverse-time characteristic. The maximum absolute values of errors of hundreds of sample points taken from four time dial settings (TDS) for each relay between the actual characteristic curves and the corresponding values from the curve-fitting equations are within the range of 10 milliseconds. Finally, this study compares the SVD with the adaptive network and fuzzy inference system (ANFIS) to demonstrate its accuracy and identification robustness.


---

## Body

## 1. Introduction

Power generation systems generally have few large generators connected directly to their subtransmission networks and distribution networks. Thus, the fault currents of buses do not differ much from those of transmission lines. This makes low-cost, reliable, and easily coordinated electromechanical (EM) overcurrent (OC) relays suitable for protection coordination relay in the subtransmission networks and distribution networks. Although some older relays have been replaced by new digital ones, there are still many EM OC relays in service.The operation principle of the EM OC relay is to introduce an electric current into the coil of an electromagnet to produce eddy currents with phase differences. This in turn generates induction torque on the rotation disc of the relay. The proper contact closing time can be set by adjusting the distance between the fixed and the movable contacts, achieving protection coordination between upstream and downstream. However, due to the mechanical nature of the relay, there are inertial and frictional effects. Therefore, unlike a digital relay [1], it is not possible to describe the characteristic curve of an EM OC relay using a single precise equation. Manufacturer manuals generally provide families of characteristic curves with different time dial (TDS) values. These are all piecewise nonlinear continuous smooth descending curves.Accurate representations of the inverse-time EM OC relay characteristics play an important role in the coordination of power network protection schemes. In the early days, researchers were interested in fitting EM OC relay characteristics curves [2] to facilitate protection coordination in conventional power systems. After the introduction of the digital OC relays, better curve-fitting of EM OC relay characteristic [3–6] is even more important for proper protection coordination. Most of the literature about curve fitting shows the absolute values of errors [7–9], while some show the averages of absolute values of errors [9, 10]. Only a few studies show the maximum absolute values of percentage errors [10], which are harder to curve fit, and no studies show the maximum absolute values of errors, which are the hardest to curve fit. For EM OC relays at small values of M (multiples of tap value current), for example, 1.3–3.5, the relay operating time changes nonlinearly and drastically, and only one study [10] shows the curve fitting results in this range of M values.This paper applies the Hankel matrices and the singular value decomposition (SVD) [11–13] to obtain the curve-fitting equation of the characteristic curves of two EM OC relays under state space with modal coordinates [14, 15]. To demonstrate the accuracy of the curve-fitting results, the current study not only shows all the maximum absolute values of errors, maximum absolute values of percentage errors, average of absolute values of errors, and average of absolute values of percentage errors, but also considers smaller values of M where the relay operating time changes nonlinearly and drastically. For even better accuracy, this paper reduces the maximum absolute values of errors to less than an alternating current cycle in the range of a few milliseconds (ms), as opposed to 3 cycles in [3] or 2 cycles in [8].The paper proposed a new application algorithm to fit the characteristic curves of the EM OC relays, using one unified equation to represent their characteristics. Finally, this study uses the SVD method to fit the characteristic curves of EM OC relays. Comparing the results with those obtained by [9] demonstrates the accuracy and identification robustness of the SVD method.The content of this paper is as follows: Section2 introduces the representations for the characteristic curves of EM OC relays, the mathematical derivation of the curve fitting of the characteristic curves of the EM OC relays using SVD method is outlined in Section 3, two cases are studied in Section 4 to verify the approach, the performance of the SVD and the well-behaved ANFIS algorithm are compared in Section 5, and the conclusion is in Section 6.


## 2. Models of the Electromechanical (EM) Overcurrent (OC) Relay Characteristic

The characteristic of an EM OC relay is determined by its magnetic circuit design, and the manufacturers provide the characteristics in the relay manuals with curves in a two-dimensional plot withM the abscissa and operating time the ordinate. Some typical models for the characteristic of an EM OC relay are as follows.


### 2.1. Exponential and Polynomial Forms

Various exponential and polynomial forms of equations are summarized and recommended by the IEEE Committee [3], for example, (2.1)–(2.5) below, for EM OC relay characteristic curve fitting. In some studies [7, 8, 10] that apply numerical methods to determine the best coefficients of the curve-fitting equations, the maximum absolute values of percentage errors are as large as 15% [10], so there is still much room for improvement. Consider
(2.1)t=Y×TDS(Mp-1)+Z=Y×TDS((i/In)p-1)+Z,(2.2)t=a0+a1×TDS+a2×TDS2+a3×TDS3+⋯,(2.3)t=b0+b1(M-1)+b2(M-1)2+b3(M-1)3+⋯,(2.4)t=c0+c1logM+c2(logM)2+c3(logM)3+⋯,(2.5)t=d0+d1×TDS+d2×TDS(M-1)2+d3×TDS2(M-1)+d4×TDS2(M-1)2+d5×TDS(M-1)3+d6×TDS2(M-1)4,(2.6)t=TDS×(F(Mp-1)+G)+H,
where t: relay operating time. TDS: time dial setting. i: fault current on the secondary side of the CT. In: current tap setting. M: multiples of tap value current, M=i/In. F,G,H,Y,Z,an,bn,cn,dn,p: constants.


### 2.2. Customized Characteristic Equation

A customized characteristic Equation (2.6) is obtained by modifying (1) in [16] for simulation. The IEEE normal standard inverse-time digital relay characteristic representation is obtained by letting K=0 [5], and the IEC normal standard inverse-time digital relay characteristic representation is obtained by letting both K=0 and B=0 [17].Take the characteristic curves of the ABB’s EM OC relay CO-8 as an example [18]. The recommended values of A, B, K, and p in (2.6) are 8.9341, 0.17966, 0.028, and 2.0938, respectively in [16]. Figure 1 shows the actual and the fitted characteristic curves with TDS settings of 0.5, 2, 5, and 10. The averages of absolute values of errors of the 488 sampling operating times for TDS settings 0.5, 2, 5, and 10 are 99.95, 189.18, 382.16, and 449.94 ms, respectively. The fitted curve differs considerably from the actual characteristic curve and cannot be used directly as a good replacement.Figure 1
The actual characteristic curves and the curves fitted by (2.6) for the CO-8 relay.Consider the case [19] in which the TDS settings in (2.6) are modified to 0.3, 1.5, 4, and 8.7, while their manufacture data counterparts remain the same as 0.5, 2, 5, and 10, respectively (expressed as Figure 2). The averages of absolute values of errors of the four sets of 488 sampling operating times taken between the modified fitted curve and the actual characteristic curve are 28.64, 32.54, 90.97, and 168.74 ms, respectively. Although this modified (2.6) has better accuracy, it is still not good enough to represent the actual characteristic curves provided by the manufacturer.Figure 2
The actual characteristic curves and the modified curves fitted by (2.6) for the CO-8 relay.


### 2.3. Data Base Method

The values ofM and the corresponding operating times are stored directly [3, 4]. This type of representation is commonly used today, but requires large data storage. Because the relay characteristics cannot be represented by an equation, interpolation method is usually applied to estimate data points not stored.


### 2.4. Artificial Intelligence Techniques

Researchers are applying more artificial neural network and fuzzy model techniques [7, 9] to optimal curve fitting. Among these approaches, the ANFIS algorithm developed by Geethanjali and Slochanal [9] shows promising results.


## 2.1. Exponential and Polynomial Forms

Various exponential and polynomial forms of equations are summarized and recommended by the IEEE Committee [3], for example, (2.1)–(2.5) below, for EM OC relay characteristic curve fitting. In some studies [7, 8, 10] that apply numerical methods to determine the best coefficients of the curve-fitting equations, the maximum absolute values of percentage errors are as large as 15% [10], so there is still much room for improvement. Consider
(2.1)t=Y×TDS(Mp-1)+Z=Y×TDS((i/In)p-1)+Z,(2.2)t=a0+a1×TDS+a2×TDS2+a3×TDS3+⋯,(2.3)t=b0+b1(M-1)+b2(M-1)2+b3(M-1)3+⋯,(2.4)t=c0+c1logM+c2(logM)2+c3(logM)3+⋯,(2.5)t=d0+d1×TDS+d2×TDS(M-1)2+d3×TDS2(M-1)+d4×TDS2(M-1)2+d5×TDS(M-1)3+d6×TDS2(M-1)4,(2.6)t=TDS×(F(Mp-1)+G)+H,
where t: relay operating time. TDS: time dial setting. i: fault current on the secondary side of the CT. In: current tap setting. M: multiples of tap value current, M=i/In. F,G,H,Y,Z,an,bn,cn,dn,p: constants.


## 2.2. Customized Characteristic Equation

A customized characteristic Equation (2.6) is obtained by modifying (1) in [16] for simulation. The IEEE normal standard inverse-time digital relay characteristic representation is obtained by letting K=0 [5], and the IEC normal standard inverse-time digital relay characteristic representation is obtained by letting both K=0 and B=0 [17].Take the characteristic curves of the ABB’s EM OC relay CO-8 as an example [18]. The recommended values of A, B, K, and p in (2.6) are 8.9341, 0.17966, 0.028, and 2.0938, respectively in [16]. Figure 1 shows the actual and the fitted characteristic curves with TDS settings of 0.5, 2, 5, and 10. The averages of absolute values of errors of the 488 sampling operating times for TDS settings 0.5, 2, 5, and 10 are 99.95, 189.18, 382.16, and 449.94 ms, respectively. The fitted curve differs considerably from the actual characteristic curve and cannot be used directly as a good replacement.Figure 1
The actual characteristic curves and the curves fitted by (2.6) for the CO-8 relay.Consider the case [19] in which the TDS settings in (2.6) are modified to 0.3, 1.5, 4, and 8.7, while their manufacture data counterparts remain the same as 0.5, 2, 5, and 10, respectively (expressed as Figure 2). The averages of absolute values of errors of the four sets of 488 sampling operating times taken between the modified fitted curve and the actual characteristic curve are 28.64, 32.54, 90.97, and 168.74 ms, respectively. Although this modified (2.6) has better accuracy, it is still not good enough to represent the actual characteristic curves provided by the manufacturer.Figure 2
The actual characteristic curves and the modified curves fitted by (2.6) for the CO-8 relay.


## 2.3. Data Base Method

The values ofM and the corresponding operating times are stored directly [3, 4]. This type of representation is commonly used today, but requires large data storage. Because the relay characteristics cannot be represented by an equation, interpolation method is usually applied to estimate data points not stored.


## 2.4. Artificial Intelligence Techniques

Researchers are applying more artificial neural network and fuzzy model techniques [7, 9] to optimal curve fitting. Among these approaches, the ANFIS algorithm developed by Geethanjali and Slochanal [9] shows promising results.


## 3. The Singular Value Decomposition (SVD) Method

Based on the concept of transfer function, this paper proposes an algorithm to represent the characteristic curves of EM OC relays to calculate Markov parameters [20, 21]. After the Hankel matrices [20] are constructed, they are decomposed by SVD and transformed to state space system under modal coordinates. Finally, the state space solution is found, and the fitting equation is obtained by transformation back to the continuous M-domain system [21]. The proposed algorithm processes the M-domain data directly and is thus an identification method in the M-domain.The SVD method is as follows [11–15, 20, 21].Step 1.
Find the estimated operating timet0 on the EM OC relay characteristic curve corresponding to a specific multiple of tap value current M- as in (3.1). Repeat this process to the right with incremental step ΔM to form the sample sequence tk, k=1,2,3,…, where
(3.1)M-=Mmin-ΔMMmin: minimum multiple of tap value current. ΔM: increment sampling step.Step 2.
Usetk, k=1,2,3,… to form the Hankel matrices H(0) and H(1) in (3.2) as follows:
(3.2)H(0)=[t1t2⋯tht2t3⋯th+1⋮⋮⋱⋮thth+1⋯t2h-1]h×h,H(1)=[t2t3⋯th+1t3t4⋯th+2⋮⋮⋱⋮th+1th+2⋯t2h]h×h.Step 3.
Apply SVD to the Hankel matrixH(0) to obtain matrices R, Σ, and S in (3.3)
(3.3)H(0)=RΣST.Step 4.
Determine the proper dimensionn for the modal coordination system in (3.3) and obtain matrices Rn, Σn, and Sn in (3.4), where n is also the number of fitting waveform components and its range is form 1 to the rank of Σ(3.4)RΣST=[RnR0][Σn0_0_Σ0][SnTS0T].Step 5.
Calculate the matricesA^, B^, and C^, which are the estimates of the matrices A, B, and C in the state space system (3.5), as (3.6) and (3.7) show
(3.5)M(k+1)=AM(k)+Bu(k)k=0,1,2,…,t(k)=CM(k)+Du(k),(3.6)A^=(Σn)-1/2RnTH(1)Sn(Σn)-1/2,(3.7)B^=(Σn)1/2SnTE1,C^=E1TRn(Σn)1/2,
where E1T is shown as (3.8)
(3.8)E1T=[10⋯0]
and the system matrix D is as shown in (3.9)
(3.9)D=t0.Step 6.
Transform the state space Equation (3.5) into modal coordinate system to find Λ, Bm, and Cm as (3.10), (3.12), (3.13), and (3.14) show
(3.10)Mm(k+1)=ΛMm(k)+Bmu(k),k=0,1,2,…,t(k)=CmMm(k)+Du(k),
where Mm(k) is as shown in (3.11)
(3.11)Mm(k)=Ψ-1x(k),(3.12)Λ=Ψ-1A^Ψ=[λ10_λ2⋱0_λn],(3.13)Bm=Ψ-1B^,(3.14)Cm=C^Ψ,
where λi: eigenvalues of A^, i=1,2,…,n. Ψ: matrix whose columns are the eigenvectors of A^.Step 7.
Obtain the unit impulse response sequence from (3.10), the state space equations under modal coordinates, as known system Markov parameters in (3.15) below:
(3.15)t(k)=t0,k=0=CAk-1B=CmΛk-1Bmk=1,2,3,…=∑i=1nciλik-1biStep 8.
Derive the equation of the fitted relay characteristic curve by transforming back to continuousM-domain system (3.16)
(3.16)t(M)=∑i=1n1Cie-αi(M-M-)+2∑i=1n2KiAri-fi(M-M-)cos(2πfi(M-M-)+φi)+∑i=(n2+1)n3KiAri-fi(M-M-)cos(2πfi(M-M-)+φi),
where t is the operating time with M as its variable. n1 is the number of the smooth waveform components. n2 is the number of the paired oscillation waveform components (thus the coefficient 2). n3 is the number of the independent oscillation components. (n3-n2) is the number of the unpaired oscillation waveform components. Ci,αi,Ki, and Ari are the constants. fi is the oscillation frequency, in Hertz. φi is the oscillation phase shift, in radian.
Equation (3.17) describes the relationship among n, n1, n2, and n3(3.17)n=n1+2n2+(n3-n2)=n1+n2+n3.The characteristic curves of an EM OC relay can be represented by a digital state space model, and the desired bound of the maximum absolute value of errors may be established by selecting an appropriate system model ordern. Equation (3.16) is an unified equation by SVD method to fit EM OC relays characteristic curves, all of which are piecewise nonlinear continuous smooth descending curves. Any of such characteristic curves can be fitted by simply changing the values of the parameters in the equation.


## 4. Cases Study

The following case study involves two EM OC relays: the GE moderately inverse-time relay IAC51 [22] and the ABB inverse-time relay CO-8 [18]. To be both accurate and reasonable, the constraint set for the fit is that the maximum absolute value of errors between all the fitted sampling points and the actual characteristic curves for each TDS be less than 10 milliseconds. The calculations in this study were made using MATLAB.


### 4.1. Case 1: GE IAC51 Relay

Four characteristic curves corresponding to TDS 1, 4, 7, and 10 were selected for curve fitting and 486 sample points of relay operating time are taken forM ranging from 1.5 to 50.0 with steps of 0.1. Table 1 lists the calculated mathematic parameters, where the selected number of fitted waveform components n is the minimum number of waveform needed to construct the TDS curve. Equation (3.16) can be rewritten as (4.1) to represent the characteristic curves corresponding to the four TDS as
(4.1)t(M)=∑i=15Cie-αi(M-1.4)+2∑i=14KiAri-fi(M-1.4)cos(2πfi(M-1.4)+φi).Table 1
The 26 parameters in (4.1) for the IAC51 relay.
TDS
1
4
7
10
n
4
9
12
13
C
1
4.7488
E
-
01
1.7515
E
+
00
3.4807
E
+
00
5.6144
E
+
00
α
1
2.4356
E
+
00
3.5435
E
+
00
3.5331
E
+
00
3.5255
E
+
00
C
2
3.3871
E
-
01
1.7317
E
+
00
3.4524
E
+
00
5.1574
E
+
00
α
2
5.9703
E
-
01
9.9106
E
-
01
9.5370
E
-
01
9.6362
E
-
01
C
3
2.4584
E
-
01
8.9496
E
-
01
1.2812
E
+
00
2.6120
E
+
00
α
3
4.0672
E
-
03
2.6039
E
-
01
5.5443
E
-
03
2.6972
E
-
01
C
4
1.7784
E
-
01
6.3240
E
-
01
1.4444
E
+
00
1.8188
E
+
00
α
4
1.2779
E
-
01
3.2351
E
-
03
1.0340
E
-
01
4.7116
E
-
03
C
5
0
4.7301
E
-
01
0
1.4377
E
+
00
α
5
0
6.2658
E
-
02
0
7.2246
E
-
02
K
1
0
3.9552
E
-
04
5.5417
E
-
01
3.4270
E
-
03
f
1
0
1.0527
E
-
01
1.7904
E
-
02
1.2421
E
-
01
A
r
1
0
1.3835
E
+
00
1.7813
E
+
05
1.9980
E
+
00
φ
1
0
-
1.9811
E
+
00
6.7820
E
-
01
1.4603
E
+
00
K
2
0
3.2805
E
-
02
8.4039
E
-
04
2.4583
E
-
03
f
2
0
4.4042
E
+
00
1.0384
E
-
01
2.2450
E
-
01
A
r
2
0
2.4185
E
+
00
1.5521
E
+
00
1.5293
E
+
00
φ
2
0
-
6.6692
E
+
00
4.0532
E
+
00
-
6.1765
E
+
00
K
3
0
0
1.7842
E
-
04
6.4674
E
-
04
f
3
0
0
1.7112
E
-
01
3.1642
E
-
01
A
r
3
0
0
1.1211
E
+
00
1.1454
E
+
00
φ
3
0
0
1.5805
E
+
00
1.7011
E
-
01
K
4
0
0
6.6388
E
-
02
1.2479
E
-
01
f
4
0
0
4.4591
E
+
00
4.5025
E
+
00
A
r
4
0
0
2.2498
E
+
00
2.4045
E
+
00
φ
4
0
0
-
4.9733
E
-
01
-
6.0588
E
-
01
TDS: time dial setting;n: number of fitted waveform components.The number of the smooth waveform componentsn1 was set to the maximum of the fitted result for the four TDS characteristic curves, and likewise for n2 and n3. The parameters of the waveform components not used in the fitted result were set to zero. This simple operating equation contains 26 parameters and includes n1=5 smooth waveform components, n2=4 paired oscillation waveform components and n3-n2=0 unpaired oscillation waveform components. Since the minimum value of the samples of the multiples of tap value current of the characteristic curve is 1.5 and the sampling step ΔM=0.1, M-=1.5-0.1=1.4.Table2 compares the 17 relay operating times and the corresponding SVD fitted values for each TDS. The average of absolute values of errors for each TDS is in the range 0.40–1.27 ms.Table 2
The curve-fitting errors of 17 actual operating times for each TDS of the IAC51 relay in milliseconds.
TDS
1
4
7
10
M
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
1.5
1115.0
1112.6
2.34
4738.9
4738.9
0.02
9067.1
9067.1
0.01
14322.0
14322.0
0.16
2
753.7
756.8
3.11
3009.2
3011.1
1.87
5700.4
5702.3
1.89
8960.3
8964.3
3.99
3
530.5
529.2
1.31
2009.6
2008.3
1.30
3772.1
3770.0
2.05
5900.0
5896.7
3.25
4
442.2
443.4
1.15
1616.0
1616.3
0.32
3013.9
3014.0
0.10
4698.1
4697.7
0.49
6
362.4
361.8
0.61
1266.8
1266.2
0.59
2342.1
2340.9
1.22
3632.3
3631.1
1.21
8
323.3
322.4
0.82
1094.7
1094.4
0.34
2010.9
2010.4
0.48
3107.9
3106.7
1.18
10
298.2
298.6
0.47
985.7
986.2
0.48
1801.4
1802.0
0.66
2776.2
2778.1
1.92
13
275.0
275.2
0.22
881.8
881.9
0.12
1601.5
1601.3
0.29
2457.8
2455.8
2.04
16
259.3
259.3
0.04
812.5
812.8
0.34
1468.8
1469.2
0.42
2249.4
2250.6
1.25
19
247.4
247.6
0.20
763.2
763.1
0.04
1373.9
1373.6
0.28
2099.8
2099.0
0.84
22
239.5
238.9
0.66
726.4
726.2
0.18
1302.1
1302.0
0.12
1985.9
1984.6
1.29
26
230.6
230.1
0.49
686.3
686.7
0.46
1224.8
1226.1
1.36
1864.8
1865.8
1.06
30
223.7
223.4
0.21
655.8
655.8
0.08
1166.7
1166.5
0.18
1772.3
1773.0
0.66
34
217.4
218.1
0.66
630.5
630.8
0.29
1118.8
1119.1
0.24
1695.6
1696.3
0.67
39
211.6
212.4
0.88
604.8
604.7
0.11
1070.0
1069.4
0.66
1619.2
1619.0
0.15
44
207.1
207.5
0.41
583.8
583.9
0.08
1029.2
1029.4
0.20
1554.7
1553.9
0.88
49
202.9
203.0
0.09
566.3
566.1
0.26
994.7
994.3
0.38
1500.4
1499.8
0.64
AV
=
0.80
AV
=
0.40
AV
=
0.62
AV
=
1.27
TDS: time dial setting;M: multiples of tap value current; top: actual operating time in milliseconds; SVD: fitted value by singular value decomposition; Err: absolute value of difference between top and SVD value; AV: average of absolute values of errors in milliseconds.Table3 summarizes the complete fitting results. The maximum absolute values of errors range from 3.62 to 8.73 ms, and all occur at smaller M values (1.8, 2.2, 2.8, and 2.2). The maximum absolute values of percentage errors range from 0.15 to 0.48, and mostly occur in the smaller one-third of the range of M (i.e., 41.4, 3.9, 3.9, and 12.1). The average of absolute values of errors range from 0.28 to 0.89 ms, and the average of absolute values of percentage errors ranges from 0.03 to 0.19.Table 3
Summary of the SVD results for the IAC51 relay.
TDS
Max_Err/M
Max_Err%/M
AV
AV%
1
4.05/1.8
0.48/41.4
0.53
0.19
4
3.62/2.2
0.16/3.9
0.28
0.03
7
5.08/2.8
0.15/3.9
0.50
0.03
10
8.73/2.2
0.20/12.1
0.89
0.04
TDS: time dial setting; Max_Err/M: maximum absolute values of errors in millisecond in multiples of tap value current; Max_Err%/M: maximum absolute values of percentage errors in multiples of tap value current; AV: average of absolute values of errors in milliseconds; AV%: averages of absolute values of percentage errors.Finally, Figure3 shows the samples of the actual relay operating times and the corresponding values from the curve-fitting equations obtained by SVD. This figure shows that two sets of values are so closely matched that they are virtually indistinguishable. This clearly demonstrates the accuracy and identification robustness of the SVD method.Figure 3
The actual characteristic curves and the curves fitted by (4.1) for the IAC51 relay.


### 4.2. Case 2: ABB CO-8 Relay

Four characteristic curves corresponding to TDS 0.5, 3, 6, and 9 were selected for curve fitting, and 488 sample points of relay operating time were taken forM ranging from 1.3 to 50.0 with steps of 0.1. Table 4 lists the calculated mathematic parameters. Equation (3.16) can be rewritten as (4.2) to represent the characteristic curves corresponding to the four TDS as
(4.2)t(M)=∑i=16Cie-αi(M-1.2)+2∑i=16KiAri-fi(M-1.2)cos(2πfi(M-1.2)+φi)+K7Ar7-f7(M-1.2)cos(2πf7(M-1.2)+φ7).Table 4
The 40 parameters in (4.2) for the CO-8 relay.
TDS
0.5
3
6
9
n
9
12
15
19
C
1
2.2260
E
+
00
2.5261
E
+
01
3.7055
E
+
01
5.8957
E
+
01
α
1
3.6161
E
+
00
3.7976
E
+
00
4.2328
E
+
00
3.5943
E
+
00
C
2
9.4399
E
-
01
4.8815
E
-
01
2.3176
E
+
01
2.9885
E
+
01
α
2
9.7475
E
-
01
3.7793
E
-
03
1.8071
E
+
00
1.5467
E
+
00
C
3
1.6429
E
-
01
5.1059
E
-
01
7.0000
E
+
00
9.1230
E
+
00
α
3
2.6196
E
-
01
1.1208
E
-
01
6.2423
E
-
01
6.0876
E
-
01
C
4
2.3958
E
-
02
1.9846
E
+
00
1.0226
E
+
00
1.6539
E
+
00
α
4
-
9.7055
E
-
03
4.1472
E
-
01
2.1111
E
-
01
2.3403
E
-
01
C
5
5.1724
E
-
02
0
9.8116
E
-
01
1.5020
E
+
00
α
5
1.9268
E
-
02
0
3.4244
E
-
03
4.0665
E
-
03
C
6
0
0
5.4251
E
-
01
1.0217
E
+
00
α
6
0
0
7.9622
E
-
02
7.9132
E
-
02
K
1
1.2639
E
-
04
3.1088
E
+
00
3.7043
E
+
00
2.9481
E
-
01
f
1
8.1404
E
-
02
6.7874
E
-
02
1.5581
E
+
00
1.0716
E
+
00
A
r
1
1.5162
E
+
00
7.4100
E
+
05
5.2041
E
+
02
1.3713
E
+
01
φ
1
-
3.4533
E
-
02
1.6373
E
-
01
1.2089
E
+
00
-
2.5188
E
+
00
K
2
8.8486
E
-
02
2.3338
E
-
02
6.1002
E
-
03
4.7967
E
-
03
f
2
3.3685
E
+
00
6.0254
E
-
01
2.8122
E
+
00
2.0997
E
+
00
A
r
2
8.2508
E
+
00
4.8391
E
+
00
1.1917
E
+
00
1.1848
E
+
00
φ
2
-
1.7509
E
+
00
-
2.5597
E
+
00
-
7.5985
E
+
00
7.1928
E
-
02
K
3
0
5.8643
E
+
00
2.1803
E
-
02
1.6860
E
-
02
f
3
0
2.8484
E
+
00
3.3973
E
+
00
2.4543
E
+
00
A
r
3
0
1.5179
E
+
02
1.3685
E
+
00
1.3082
E
+
00
φ
3
0
-
6.8235
E
+
00
3.3636
E
+
00
-
1.1421
E
+
00
K
4
0
1.5217
E
-
01
4.9041
E
-
01
3.8013
E
-
02
f
4
0
4.3811
E
+
00
4.2044
E
+
00
2.8364
E
+
00
A
r
4
0
1.7508
E
+
00
1.9500
E
+
00
1.3904
E
+
00
φ
4
0
4.6268
E
-
01
-
4.5837
E
+
00
-
1.5897
E
+
00
K
5
0
0
0
1.1831
E
-
01
f
5
0
0
0
3.3797
E
+
00
A
r
5
0
0
0
1.6111
E
+
00
φ
5
0
0
0
-
2.6524
E
+
00
K
6
0
0
0
7.6744
E
-
01
f
6
0
0
0
4.0923
E
+
00
A
r
6
0
0
0
1.9734
E
+
00
φ
6
0
0
0
-
3.8092
E
+
00
K
7
0
0
2.8143
E
+
00
2.9320
E
+
00
f
7
0
0
5.0000
E
+
00
5.0000
E
+
00
A
r
7
0
0
3.5972
E
+
00
2.7190
E
+
00
φ
7
0
0
-
3.1416
E
+
00
-
3.1416
E
+
00
Please refer to Table1 for the meanings of the terms.This simple operating equation contains 40 parameters and includesn1=6 smooth waveform components, n2=6 paired oscillation waveform components, and n3-n2=1 unpaired oscillation waveform components. Since the minimum value of the samples of the multiples of tap value current of the characteristic curve is 1.3 and the sampling step ΔM=0.1, M-=1.3-0.1=1.2.Table5 compares the 17 relay operating times and the corresponding SVD fitted values for each TDS. The average of absolute values of errors for each TDS ranges from 0.37 to 0.95 ms.Table 5
The curve-fitting errors of 17 actual operating times for each TDS of the CO-8 relay in milliseconds.
TDS
0.5
3
6
9
M
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
1.5
1682.2
1682.1
0.11
15254.0
15254.0
0.10
31907.0
31907.0
0.26
50681.0
50681.0
0.16
2
760.5
764.2
3.69
6182.1
6182.6
0.49
13232.0
13232.0
0.21
21256.0
21257.0
1.11
3
343.9
343.6
0.25
2577.6
2581.6
3.98
5334.0
5337.1
3.10
8449.4
8450.2
0.78
4
213.7
214.3
0.60
1578.5
1578.0
0.44
3341.9
3339.4
2.51
5222.6
5218.7
3.92
6
128.1
127.6
0.55
1004.0
1002.7
1.24
2062.6
2060.3
2.32
3221.3
3218.3
3.02
8
99.7
99.7
0.01
821.3
820.2
1.12
1618.9
1618.4
0.44
2542.7
2540.0
2.66
10
86.1
86.3
0.19
711.7
712.9
1.13
1407.3
1409.4
2.16
2210.1
2211.9
1.79
13
75.8
75.7
0.04
618.3
617.9
0.42
1244.0
1243.5
0.56
1944.8
1944.7
0.10
16
70.1
70.0
0.10
563.1
563.1
0.01
1145.0
1145.3
0.30
1783.3
1783.9
0.62
19
66.5
66.6
0.13
526.5
527.1
0.56
1078.3
1078.6
0.28
1672.8
1672.7
0.06
22
64.7
64.6
0.07
501.4
501.2
0.12
1030.0
1029.9
0.06
1589.8
1589.9
0.13
26
62.9
62.9
0.01
476.0
476.2
0.23
982.0
982.0
0.02
1506.4
1506.5
0.02
30
61.4
61.4
0.00
458.5
458.1
0.47
946.5
946.1
0.36
1442.9
1442.6
0.33
34
60.1
60.4
0.35
443.7
444.2
0.45
917.4
917.8
0.38
1390.5
1391.4
0.93
39
59.6
59.6
0.02
430.7
430.6
0.19
889.1
889.1
0.02
1339.7
1339.5
0.14
44
59.0
58.9
0.09
419.3
419.5
0.20
865.5
865.5
0.03
1297.0
1296.7
0.27
49
58.6
58.7
0.12
409.8
409.9
0.11
845.0
845.1
0.14
1260.1
1259.9
0.15
AV
=
0.37
AV
=
0.66
AV
=
0.77
AV
=
0.95
Please refer to Table2 for the meanings of the terms.Table6 summarizes the complete fitting results. The maximum absolute values of errors range from 3.69 to 8.22 ms, and all occur at smaller M values (i.e., 2.0, 3.9, 2.5, and 4.4). The maximum absolute values of percentage errors ranges from 0.18 to 0.73, and also occur at smaller M values. The average of absolute values of errors ranges from 0.17 to 0.58 ms, while the averages of absolute values of percentage errors ranges from 0.03 to 0.16.Table 6
Summary of the SVD results for the CO-8 relay.
TDS
Max_Err/M
Max_Err%/M
AV
AV%
0.5
3.69/2.0
0.73/3.9
0.17
0.16
3
8.22/3.9
0.50/3.9
0.42
0.05
6
7.72/2.5
0.29/5.2
0.53
0.03
9
6.08/4.4
0.18/6.9
0.58
0.03
Please refer to Table3 for the meanings of the terms.The absolute values of error of 488 sample points with TDS 3 CO8 relay are shown in Figure4. It can be seen that the absolute values of errors are within preset range, and the absolute values of errors at smaller values of M are larger, in accordance with what was mentioned in Section 1 that, for EM OC relays, curve fitting is harder at small values of M where the relay operating time changes nonlinearly and drastically.Figure 4
Absolute values of error forTDS=3 CO-8 relay.Finally, Figure5 shows the samples of the actual relay operating times and the corresponding values from the curve-fitting equations obtained by SVD. This figure shows that these two sets of values are also closely matched.Figure 5
The actual characteristic curves and the curves fitted by (4.2) for the CO-8 relay.


### 4.3. Result Analysis

The following paragraphs analyze and compare the SVD fitting results for the IAC51 and CO-8 EM OC relays.The maximum absolute values of errors of the hundreds of sample points calculated by SVD are within the millisecond range, which is accurate enough for practical purposes. However, more fitted waveform components can be selected if greater precision is desired. Although the maximum absolute values of errors and maximum absolute values of percentage errors primarily occur at smallM values where the relay operating time is highly nonlinear and changes rapidly, the curves fitted by SVD method still match the relay characteristic very well in the case study.The SVD method can be applied to different types of inverse-time EM OC relay. As shown in Table1 for the IAC51 relay, the numbers of fitted waveform components for the four characteristic curves with TDS 1, 4, 7, and 10 are 4, 9, 12, and 13, respectively. Table 4 shows the numbers of fitted waveform components for the four characteristic curves with TDS 0.5, 3, 6, and 9 are 9, 12, 15, and 19, respectively for the CO-8 relay. Relays with more inverse-time characteristics require more fitted waveform components to achieve the same range of maximum absolute values of errors. For a particular inverse-time EM OC relay, a larger number of fitted waveform components is also required to achieve the same range of maximum absolute values of errors for higher TDS. In any case, the fitted results for both types of relays are accurate enough for practical applications.SVD can be applied to achieve excellent protection coordination between widely used digital relays and conventional EM OC relays. The procedures may be outlined below using CO-8 relay as an example, with a coordination time interval of 0.3 second.(a) Use the CO-8 Relay as the Primary Relay
SetTDS=5 so that (t,M) = (1.098, 18). The setting of the backup digital relay should be (t,M)=(1.398,18) and can be easily accomplished by setting TDS=6.49. This makes it easy to set the digital relay using the curve-fitted equation of the EM OC relay’s characteristic curves with TDS=6 and TDS=7 as a base reference.(b) Use the Digital Relay to Achieve Adaptive Relay [1, 8] Protection Coordination
First, apply SVD to determine the settings of the CO-8 relay and then calculate the operating time for the digital relay by adding 0.3 s to, or subtracting 0.3 s from, that of the CO-8 relay if the digital relay is to function as a backup or primary relay, respectively.


## 4.1. Case 1: GE IAC51 Relay

Four characteristic curves corresponding to TDS 1, 4, 7, and 10 were selected for curve fitting and 486 sample points of relay operating time are taken forM ranging from 1.5 to 50.0 with steps of 0.1. Table 1 lists the calculated mathematic parameters, where the selected number of fitted waveform components n is the minimum number of waveform needed to construct the TDS curve. Equation (3.16) can be rewritten as (4.1) to represent the characteristic curves corresponding to the four TDS as
(4.1)t(M)=∑i=15Cie-αi(M-1.4)+2∑i=14KiAri-fi(M-1.4)cos(2πfi(M-1.4)+φi).Table 1
The 26 parameters in (4.1) for the IAC51 relay.
TDS
1
4
7
10
n
4
9
12
13
C
1
4.7488
E
-
01
1.7515
E
+
00
3.4807
E
+
00
5.6144
E
+
00
α
1
2.4356
E
+
00
3.5435
E
+
00
3.5331
E
+
00
3.5255
E
+
00
C
2
3.3871
E
-
01
1.7317
E
+
00
3.4524
E
+
00
5.1574
E
+
00
α
2
5.9703
E
-
01
9.9106
E
-
01
9.5370
E
-
01
9.6362
E
-
01
C
3
2.4584
E
-
01
8.9496
E
-
01
1.2812
E
+
00
2.6120
E
+
00
α
3
4.0672
E
-
03
2.6039
E
-
01
5.5443
E
-
03
2.6972
E
-
01
C
4
1.7784
E
-
01
6.3240
E
-
01
1.4444
E
+
00
1.8188
E
+
00
α
4
1.2779
E
-
01
3.2351
E
-
03
1.0340
E
-
01
4.7116
E
-
03
C
5
0
4.7301
E
-
01
0
1.4377
E
+
00
α
5
0
6.2658
E
-
02
0
7.2246
E
-
02
K
1
0
3.9552
E
-
04
5.5417
E
-
01
3.4270
E
-
03
f
1
0
1.0527
E
-
01
1.7904
E
-
02
1.2421
E
-
01
A
r
1
0
1.3835
E
+
00
1.7813
E
+
05
1.9980
E
+
00
φ
1
0
-
1.9811
E
+
00
6.7820
E
-
01
1.4603
E
+
00
K
2
0
3.2805
E
-
02
8.4039
E
-
04
2.4583
E
-
03
f
2
0
4.4042
E
+
00
1.0384
E
-
01
2.2450
E
-
01
A
r
2
0
2.4185
E
+
00
1.5521
E
+
00
1.5293
E
+
00
φ
2
0
-
6.6692
E
+
00
4.0532
E
+
00
-
6.1765
E
+
00
K
3
0
0
1.7842
E
-
04
6.4674
E
-
04
f
3
0
0
1.7112
E
-
01
3.1642
E
-
01
A
r
3
0
0
1.1211
E
+
00
1.1454
E
+
00
φ
3
0
0
1.5805
E
+
00
1.7011
E
-
01
K
4
0
0
6.6388
E
-
02
1.2479
E
-
01
f
4
0
0
4.4591
E
+
00
4.5025
E
+
00
A
r
4
0
0
2.2498
E
+
00
2.4045
E
+
00
φ
4
0
0
-
4.9733
E
-
01
-
6.0588
E
-
01
TDS: time dial setting;n: number of fitted waveform components.The number of the smooth waveform componentsn1 was set to the maximum of the fitted result for the four TDS characteristic curves, and likewise for n2 and n3. The parameters of the waveform components not used in the fitted result were set to zero. This simple operating equation contains 26 parameters and includes n1=5 smooth waveform components, n2=4 paired oscillation waveform components and n3-n2=0 unpaired oscillation waveform components. Since the minimum value of the samples of the multiples of tap value current of the characteristic curve is 1.5 and the sampling step ΔM=0.1, M-=1.5-0.1=1.4.Table2 compares the 17 relay operating times and the corresponding SVD fitted values for each TDS. The average of absolute values of errors for each TDS is in the range 0.40–1.27 ms.Table 2
The curve-fitting errors of 17 actual operating times for each TDS of the IAC51 relay in milliseconds.
TDS
1
4
7
10
M
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
1.5
1115.0
1112.6
2.34
4738.9
4738.9
0.02
9067.1
9067.1
0.01
14322.0
14322.0
0.16
2
753.7
756.8
3.11
3009.2
3011.1
1.87
5700.4
5702.3
1.89
8960.3
8964.3
3.99
3
530.5
529.2
1.31
2009.6
2008.3
1.30
3772.1
3770.0
2.05
5900.0
5896.7
3.25
4
442.2
443.4
1.15
1616.0
1616.3
0.32
3013.9
3014.0
0.10
4698.1
4697.7
0.49
6
362.4
361.8
0.61
1266.8
1266.2
0.59
2342.1
2340.9
1.22
3632.3
3631.1
1.21
8
323.3
322.4
0.82
1094.7
1094.4
0.34
2010.9
2010.4
0.48
3107.9
3106.7
1.18
10
298.2
298.6
0.47
985.7
986.2
0.48
1801.4
1802.0
0.66
2776.2
2778.1
1.92
13
275.0
275.2
0.22
881.8
881.9
0.12
1601.5
1601.3
0.29
2457.8
2455.8
2.04
16
259.3
259.3
0.04
812.5
812.8
0.34
1468.8
1469.2
0.42
2249.4
2250.6
1.25
19
247.4
247.6
0.20
763.2
763.1
0.04
1373.9
1373.6
0.28
2099.8
2099.0
0.84
22
239.5
238.9
0.66
726.4
726.2
0.18
1302.1
1302.0
0.12
1985.9
1984.6
1.29
26
230.6
230.1
0.49
686.3
686.7
0.46
1224.8
1226.1
1.36
1864.8
1865.8
1.06
30
223.7
223.4
0.21
655.8
655.8
0.08
1166.7
1166.5
0.18
1772.3
1773.0
0.66
34
217.4
218.1
0.66
630.5
630.8
0.29
1118.8
1119.1
0.24
1695.6
1696.3
0.67
39
211.6
212.4
0.88
604.8
604.7
0.11
1070.0
1069.4
0.66
1619.2
1619.0
0.15
44
207.1
207.5
0.41
583.8
583.9
0.08
1029.2
1029.4
0.20
1554.7
1553.9
0.88
49
202.9
203.0
0.09
566.3
566.1
0.26
994.7
994.3
0.38
1500.4
1499.8
0.64
AV
=
0.80
AV
=
0.40
AV
=
0.62
AV
=
1.27
TDS: time dial setting;M: multiples of tap value current; top: actual operating time in milliseconds; SVD: fitted value by singular value decomposition; Err: absolute value of difference between top and SVD value; AV: average of absolute values of errors in milliseconds.Table3 summarizes the complete fitting results. The maximum absolute values of errors range from 3.62 to 8.73 ms, and all occur at smaller M values (1.8, 2.2, 2.8, and 2.2). The maximum absolute values of percentage errors range from 0.15 to 0.48, and mostly occur in the smaller one-third of the range of M (i.e., 41.4, 3.9, 3.9, and 12.1). The average of absolute values of errors range from 0.28 to 0.89 ms, and the average of absolute values of percentage errors ranges from 0.03 to 0.19.Table 3
Summary of the SVD results for the IAC51 relay.
TDS
Max_Err/M
Max_Err%/M
AV
AV%
1
4.05/1.8
0.48/41.4
0.53
0.19
4
3.62/2.2
0.16/3.9
0.28
0.03
7
5.08/2.8
0.15/3.9
0.50
0.03
10
8.73/2.2
0.20/12.1
0.89
0.04
TDS: time dial setting; Max_Err/M: maximum absolute values of errors in millisecond in multiples of tap value current; Max_Err%/M: maximum absolute values of percentage errors in multiples of tap value current; AV: average of absolute values of errors in milliseconds; AV%: averages of absolute values of percentage errors.Finally, Figure3 shows the samples of the actual relay operating times and the corresponding values from the curve-fitting equations obtained by SVD. This figure shows that two sets of values are so closely matched that they are virtually indistinguishable. This clearly demonstrates the accuracy and identification robustness of the SVD method.Figure 3
The actual characteristic curves and the curves fitted by (4.1) for the IAC51 relay.


## 4.2. Case 2: ABB CO-8 Relay

Four characteristic curves corresponding to TDS 0.5, 3, 6, and 9 were selected for curve fitting, and 488 sample points of relay operating time were taken forM ranging from 1.3 to 50.0 with steps of 0.1. Table 4 lists the calculated mathematic parameters. Equation (3.16) can be rewritten as (4.2) to represent the characteristic curves corresponding to the four TDS as
(4.2)t(M)=∑i=16Cie-αi(M-1.2)+2∑i=16KiAri-fi(M-1.2)cos(2πfi(M-1.2)+φi)+K7Ar7-f7(M-1.2)cos(2πf7(M-1.2)+φ7).Table 4
The 40 parameters in (4.2) for the CO-8 relay.
TDS
0.5
3
6
9
n
9
12
15
19
C
1
2.2260
E
+
00
2.5261
E
+
01
3.7055
E
+
01
5.8957
E
+
01
α
1
3.6161
E
+
00
3.7976
E
+
00
4.2328
E
+
00
3.5943
E
+
00
C
2
9.4399
E
-
01
4.8815
E
-
01
2.3176
E
+
01
2.9885
E
+
01
α
2
9.7475
E
-
01
3.7793
E
-
03
1.8071
E
+
00
1.5467
E
+
00
C
3
1.6429
E
-
01
5.1059
E
-
01
7.0000
E
+
00
9.1230
E
+
00
α
3
2.6196
E
-
01
1.1208
E
-
01
6.2423
E
-
01
6.0876
E
-
01
C
4
2.3958
E
-
02
1.9846
E
+
00
1.0226
E
+
00
1.6539
E
+
00
α
4
-
9.7055
E
-
03
4.1472
E
-
01
2.1111
E
-
01
2.3403
E
-
01
C
5
5.1724
E
-
02
0
9.8116
E
-
01
1.5020
E
+
00
α
5
1.9268
E
-
02
0
3.4244
E
-
03
4.0665
E
-
03
C
6
0
0
5.4251
E
-
01
1.0217
E
+
00
α
6
0
0
7.9622
E
-
02
7.9132
E
-
02
K
1
1.2639
E
-
04
3.1088
E
+
00
3.7043
E
+
00
2.9481
E
-
01
f
1
8.1404
E
-
02
6.7874
E
-
02
1.5581
E
+
00
1.0716
E
+
00
A
r
1
1.5162
E
+
00
7.4100
E
+
05
5.2041
E
+
02
1.3713
E
+
01
φ
1
-
3.4533
E
-
02
1.6373
E
-
01
1.2089
E
+
00
-
2.5188
E
+
00
K
2
8.8486
E
-
02
2.3338
E
-
02
6.1002
E
-
03
4.7967
E
-
03
f
2
3.3685
E
+
00
6.0254
E
-
01
2.8122
E
+
00
2.0997
E
+
00
A
r
2
8.2508
E
+
00
4.8391
E
+
00
1.1917
E
+
00
1.1848
E
+
00
φ
2
-
1.7509
E
+
00
-
2.5597
E
+
00
-
7.5985
E
+
00
7.1928
E
-
02
K
3
0
5.8643
E
+
00
2.1803
E
-
02
1.6860
E
-
02
f
3
0
2.8484
E
+
00
3.3973
E
+
00
2.4543
E
+
00
A
r
3
0
1.5179
E
+
02
1.3685
E
+
00
1.3082
E
+
00
φ
3
0
-
6.8235
E
+
00
3.3636
E
+
00
-
1.1421
E
+
00
K
4
0
1.5217
E
-
01
4.9041
E
-
01
3.8013
E
-
02
f
4
0
4.3811
E
+
00
4.2044
E
+
00
2.8364
E
+
00
A
r
4
0
1.7508
E
+
00
1.9500
E
+
00
1.3904
E
+
00
φ
4
0
4.6268
E
-
01
-
4.5837
E
+
00
-
1.5897
E
+
00
K
5
0
0
0
1.1831
E
-
01
f
5
0
0
0
3.3797
E
+
00
A
r
5
0
0
0
1.6111
E
+
00
φ
5
0
0
0
-
2.6524
E
+
00
K
6
0
0
0
7.6744
E
-
01
f
6
0
0
0
4.0923
E
+
00
A
r
6
0
0
0
1.9734
E
+
00
φ
6
0
0
0
-
3.8092
E
+
00
K
7
0
0
2.8143
E
+
00
2.9320
E
+
00
f
7
0
0
5.0000
E
+
00
5.0000
E
+
00
A
r
7
0
0
3.5972
E
+
00
2.7190
E
+
00
φ
7
0
0
-
3.1416
E
+
00
-
3.1416
E
+
00
Please refer to Table1 for the meanings of the terms.This simple operating equation contains 40 parameters and includesn1=6 smooth waveform components, n2=6 paired oscillation waveform components, and n3-n2=1 unpaired oscillation waveform components. Since the minimum value of the samples of the multiples of tap value current of the characteristic curve is 1.3 and the sampling step ΔM=0.1, M-=1.3-0.1=1.2.Table5 compares the 17 relay operating times and the corresponding SVD fitted values for each TDS. The average of absolute values of errors for each TDS ranges from 0.37 to 0.95 ms.Table 5
The curve-fitting errors of 17 actual operating times for each TDS of the CO-8 relay in milliseconds.
TDS
0.5
3
6
9
M
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
Top
SVD
Err
1.5
1682.2
1682.1
0.11
15254.0
15254.0
0.10
31907.0
31907.0
0.26
50681.0
50681.0
0.16
2
760.5
764.2
3.69
6182.1
6182.6
0.49
13232.0
13232.0
0.21
21256.0
21257.0
1.11
3
343.9
343.6
0.25
2577.6
2581.6
3.98
5334.0
5337.1
3.10
8449.4
8450.2
0.78
4
213.7
214.3
0.60
1578.5
1578.0
0.44
3341.9
3339.4
2.51
5222.6
5218.7
3.92
6
128.1
127.6
0.55
1004.0
1002.7
1.24
2062.6
2060.3
2.32
3221.3
3218.3
3.02
8
99.7
99.7
0.01
821.3
820.2
1.12
1618.9
1618.4
0.44
2542.7
2540.0
2.66
10
86.1
86.3
0.19
711.7
712.9
1.13
1407.3
1409.4
2.16
2210.1
2211.9
1.79
13
75.8
75.7
0.04
618.3
617.9
0.42
1244.0
1243.5
0.56
1944.8
1944.7
0.10
16
70.1
70.0
0.10
563.1
563.1
0.01
1145.0
1145.3
0.30
1783.3
1783.9
0.62
19
66.5
66.6
0.13
526.5
527.1
0.56
1078.3
1078.6
0.28
1672.8
1672.7
0.06
22
64.7
64.6
0.07
501.4
501.2
0.12
1030.0
1029.9
0.06
1589.8
1589.9
0.13
26
62.9
62.9
0.01
476.0
476.2
0.23
982.0
982.0
0.02
1506.4
1506.5
0.02
30
61.4
61.4
0.00
458.5
458.1
0.47
946.5
946.1
0.36
1442.9
1442.6
0.33
34
60.1
60.4
0.35
443.7
444.2
0.45
917.4
917.8
0.38
1390.5
1391.4
0.93
39
59.6
59.6
0.02
430.7
430.6
0.19
889.1
889.1
0.02
1339.7
1339.5
0.14
44
59.0
58.9
0.09
419.3
419.5
0.20
865.5
865.5
0.03
1297.0
1296.7
0.27
49
58.6
58.7
0.12
409.8
409.9
0.11
845.0
845.1
0.14
1260.1
1259.9
0.15
AV
=
0.37
AV
=
0.66
AV
=
0.77
AV
=
0.95
Please refer to Table2 for the meanings of the terms.Table6 summarizes the complete fitting results. The maximum absolute values of errors range from 3.69 to 8.22 ms, and all occur at smaller M values (i.e., 2.0, 3.9, 2.5, and 4.4). The maximum absolute values of percentage errors ranges from 0.18 to 0.73, and also occur at smaller M values. The average of absolute values of errors ranges from 0.17 to 0.58 ms, while the averages of absolute values of percentage errors ranges from 0.03 to 0.16.Table 6
Summary of the SVD results for the CO-8 relay.
TDS
Max_Err/M
Max_Err%/M
AV
AV%
0.5
3.69/2.0
0.73/3.9
0.17
0.16
3
8.22/3.9
0.50/3.9
0.42
0.05
6
7.72/2.5
0.29/5.2
0.53
0.03
9
6.08/4.4
0.18/6.9
0.58
0.03
Please refer to Table3 for the meanings of the terms.The absolute values of error of 488 sample points with TDS 3 CO8 relay are shown in Figure4. It can be seen that the absolute values of errors are within preset range, and the absolute values of errors at smaller values of M are larger, in accordance with what was mentioned in Section 1 that, for EM OC relays, curve fitting is harder at small values of M where the relay operating time changes nonlinearly and drastically.Figure 4
Absolute values of error forTDS=3 CO-8 relay.Finally, Figure5 shows the samples of the actual relay operating times and the corresponding values from the curve-fitting equations obtained by SVD. This figure shows that these two sets of values are also closely matched.Figure 5
The actual characteristic curves and the curves fitted by (4.2) for the CO-8 relay.


## 4.3. Result Analysis

The following paragraphs analyze and compare the SVD fitting results for the IAC51 and CO-8 EM OC relays.The maximum absolute values of errors of the hundreds of sample points calculated by SVD are within the millisecond range, which is accurate enough for practical purposes. However, more fitted waveform components can be selected if greater precision is desired. Although the maximum absolute values of errors and maximum absolute values of percentage errors primarily occur at smallM values where the relay operating time is highly nonlinear and changes rapidly, the curves fitted by SVD method still match the relay characteristic very well in the case study.The SVD method can be applied to different types of inverse-time EM OC relay. As shown in Table1 for the IAC51 relay, the numbers of fitted waveform components for the four characteristic curves with TDS 1, 4, 7, and 10 are 4, 9, 12, and 13, respectively. Table 4 shows the numbers of fitted waveform components for the four characteristic curves with TDS 0.5, 3, 6, and 9 are 9, 12, 15, and 19, respectively for the CO-8 relay. Relays with more inverse-time characteristics require more fitted waveform components to achieve the same range of maximum absolute values of errors. For a particular inverse-time EM OC relay, a larger number of fitted waveform components is also required to achieve the same range of maximum absolute values of errors for higher TDS. In any case, the fitted results for both types of relays are accurate enough for practical applications.SVD can be applied to achieve excellent protection coordination between widely used digital relays and conventional EM OC relays. The procedures may be outlined below using CO-8 relay as an example, with a coordination time interval of 0.3 second.(a) Use the CO-8 Relay as the Primary Relay
SetTDS=5 so that (t,M) = (1.098, 18). The setting of the backup digital relay should be (t,M)=(1.398,18) and can be easily accomplished by setting TDS=6.49. This makes it easy to set the digital relay using the curve-fitted equation of the EM OC relay’s characteristic curves with TDS=6 and TDS=7 as a base reference.(b) Use the Digital Relay to Achieve Adaptive Relay [1, 8] Protection Coordination
First, apply SVD to determine the settings of the CO-8 relay and then calculate the operating time for the digital relay by adding 0.3 s to, or subtracting 0.3 s from, that of the CO-8 relay if the digital relay is to function as a backup or primary relay, respectively.


## 5. Comparison between SVD and ANFIS

In 2008, Geethanjali and Slochanal applied fuzzy logic and artificial neural network on six adaptive network and fuzzy inference system (ANFIS) algorithms (tri-mf5, gauss-mf5, gbell-mf5, tri-mf7, gauss-mf7, and gbell-mf7) for the characteristic curve fitting of one EM OC relay: the very inverse-time CRP9 relay withTDS=4, 7, and 10 [9]. The current study applies the SVD method to the relay for the same TDS characteristic curves for comparison purposes, with sampling step 1.0, and 18 relay operating times [9]. Table 7 lists the calculated mathematic parameters. Equation (5.1) represents the characteristic curves as
(5.1)t(M)=∑i=13Cie-αi(M-M-)+2∑i=13KiAri-fi(M-M-)cos(2πfi(M-M-)+φi)+∑i=45KiAri-fi(M-M-)cos(2πfi(M-M-)+φi).Table 7
The 26 parameters in (5.1) for the CRP9 relay.
TDS
4
7
10
n
9
9
9
C
1
7.7404
E
+
00
1.0543
E
+
03
8.2698
E
+
02
α
1
-
1.6459
E
-
01
2.6177
E
-
02
-
5.0138
E
-
03
C
2
9.5174
E
+
02
2.2089
E
+
03
2.3671
E
+
03
α
2
8.3416
E
-
02
3.5173
E
-
01
1.6485
E
-
01
C
3
5.0011
E
+
03
7.2570
E
+
03
1.0160
E
+
04
α
3
1.1198
E
+
00
1.4926
E
+
00
1.0356
E
+
00
K
1
6.9836
E
+
00
1.1067
E
+
01
5.1534
E
+
01
f
1
1.7548
E
-
01
1.5364
E
-
01
1.8441
E
-
01
A
r
1
1.7617
E
+
00
1.2087
E
+
00
2.6035
E
+
00
φ
1
-
3.9473
E
+
00
-
2.8492
E
+
00
-
2.2634
E
+
00
K
2
3.3350
E
+
01
2.4708
E
-
03
3.9790
E
+
01
f
2
2.9238
E
-
01
3.3275
E
-
01
2.9308
E
-
01
A
r
2
3.4628
E
+
00
2.8111
E
-
01
2.2240
E
+
00
φ
2
9.8357
E
-
01
2.2210
E
-
01
-
3.2530
E
+
00
K
3
0
0
3.8059
E
-
01
f
3
0
0
4.5246
E
-
01
A
r
3
0
0
7.3894
E
-
01
φ
3
0
0
-
3.3901
E
+
00
K
4
1.4208
E
-
03
5.0216
E
+
00
0
f
4
5.0000
E
-
01
5.0000
E
-
01
0
A
r
4
4.6291
E
-
01
1.0995
E
+
00
0
φ
4
-
3.1416
E
+
00
0.0000
E
+
00
0
K
5
1.3780
E
+
01
5.2586
E
+
01
0
f
5
5.0000
E
-
01
5.0000
E
-
01
0
A
r
5
2.0196
E
+
00
2.8355
E
+
00
0
φ
5
-
6.2832
E
+
00
-
3.1416
E
+
00
0
Please refer to Table1 for the meanings of the terms.This simple operating equation includesn1=3 smooth waveform components, n2=3 paired oscillation waveform components, and n3-n2=2 unpaired oscillation waveform components. The minimum value of the samples of the multiples of tap value current of the characteristic curve is 3 for the CRP9 EM OC relay. Thus, the value of M- is 2.For each TDS, Table8 compares the best results of each of the six ANFIS and those of SVD. For the CRP9 relay, the table shows that the ranges of maximum absolute values of errors are 0.00 ms and 16.00–28.00 ms, the ranges of maximum absolute values of percentage errors are 0.00 and 0.91–1.82, the ranges of averages of absolute values of errors are 0.00 ms and 5.20–11.80 ms and the ranges of averages of absolute values of percentage errors are 0.00 and 0.61–1.59, for the SVD and ANFIS, respectively.Table 8
Comparison of results by ANFIS and SVD for the CRP9 relay.
M
Top
ANFIS
Err1
Err1%
SVD
Err
Err%
TDS
=
4
4
1352
1368
16
1.18
1352.00
0.00
0.00
8
605
609
4
0.66
605.00
0.00
0.00
12
459
456
3
0.65
459.00
0.00
0.00
16
375
374
1
0.27
375.00
0.00
0.00
20
358
356
2
0.56
358.00
0.00
0.00
AV = 5.2
AV% = 1.59★
AV = 0.00
AV% = 0.00
TDS
=
7
4
2459
2477
18
0.73
2459.00
0.00
0.00
8
1153
1132
21
1.82
1153.00
0.00
0.00
12
894
883
11
1.23
894.00
0.00
0.00
16
744
739
5
0.67
744.00
0.00
0.00
20
669
665
4
0.60
669.00
0.00
0.00
AV = 11.8
AV% = 0.61★
AV = 0.00
AV% = 0.00
TDS
=
10
4
3935
3963
28
0.71
3935.00
0.00
0.00
8
1754
1738
16
0.91
1754.00
0.00
0.00
12
1304
1305
1
0.08
1304.00
0.00
0.00
16
1123
1124
1
0.09
1123.00
0.00
0.00
20
1024
1027
3
0.29
1024.00
0.00
0.00
AV = 9.8
AV% = 0.64★
AV = 0.00
AV% = 0.00
M: multiples of tap value current; TDS: time dial setting; top: actual operating time in milliseconds; ANFIS: fitted value by adaptive network and fuzzy inference system; Err1: absolute value of difference between top and ANFIS value; Err1%: absolute value of difference between top and ANFIS in percentage; SVD: fitted value by singular value decomposition; Err: absolute value of difference between top and SVD value; Err%: absolute value of difference between top and SVD value in percentage; AV: average of absolute values of errors in milliseconds; AV%: average of absolute values of percentage errors; ★refer to [9].Nine methods, which include the coefficient analytical model and fuzzy model in [7], SVD method and six ANFIS methods [9], are used to fit the same three TDS characteristic curves for CRP9 relay. The average values of absolute percentage errors are shown in Figure 6. SVD is clearly the best.Figure 6
Average of absolute values of percentage errors of coefficient analytical model, fuzzy model, SVD method, and six ANFIS methods for CRP9 relay.In summary, SVD outperforms ANFIS in curve fitting the operating time characteristics of the EM OC relay CRP9. Once again, this confirms the accuracy and identification robustness of the SVD method.


## 6. Conclusion

This study proposes an algorithm based on singular value decomposition (SVD) and fits the characteristic curves corresponding to eight TDS of two different types of inverse-time EM OC relays. The method decomposes the waveforms of the curves, according to their eigenvalues and corresponding eigenvectors, into various smooth and oscillating components without converting them to the frequency domain. Results show that the SVD performs exceedingly well in every respect regarding four different types of absolute values of errors. Although the maximum absolute values of errors and maximum absolute values of percentage errors primarily occur at smallM values where the relay operating time is highly nonlinear and changes rapidly, the curves fitted by the SVD method can still match the relay characteristic very well. The SVD is also superior to ANFIS.The ability to accurately represent the characteristic curves of the EM OC relays by a simple equation can provide good protection coordination between conventional EM OC relays and digital relays for subtransmission systems and distribution systems. The formula is a unified equation by SVD algorithm, and it can fit the characteristics of any piecewise nonlinear continuous smooth descending curves by simply changing the values of its parameters. Such convenience or advantage can be exploited not only for OC protection equipments such as EM OC relays or power fuses in power systems for practical applications or future studies, but also for any subject in any field with such characteristic.



---

*Source: 104952-2012-12-17.xml*