Calculation of Transmission Line Parameters From Synchronized Measurements
Calculation of Transmission Line Parameters From Synchronized Measurements
Calculation of Transmission
Line Parameters From
Synchronized Measurements
Robert E. Wilson, Gary A. Zevenbergen, Daniel
L. Mah, A. Jay Murphy
Published online: 30 Nov 2010.
To cite this article: Robert E. Wilson, Gary A. Zevenbergen, Daniel L. Mah, A. Jay
Murphy (1999) Calculation of Transmission Line Parameters From Synchronized
Measurements, Electric Machines & Power Systems, 27:12, 1269-1278, DOI:
10.1080/073135699268560
To link to this article: http://dx.doi.org/10.1080/073135699268560
A. JAY MURPHY
Macrodyne, Inc.
Clifton Park, NY, USA
Accurate knowledge of transmission system parameters, such as series impedance, optimizes distance relay settings and impedance-based fault location.
A new method is developed to measure transmission line impedances and admittances from synchronized phasor measurements. Power system voltage and
current phasors were recorded during commissioning tests of a 525 kV transmission system containing shunt and series compensation. Steady-state phase
angles ( d ) and magnitudes of line-side voltages and currents were measured
relative to a satellite-based global time standard at two substations. Pairs of
synchronized measurements were used to directly calculate real-time steadystate ABCD parameters. For comparison, two-port ABCD transmission system
parameters were derived by classical methods and from an ElectroMagnetics
Transients Program (EMTP) transmission system design model.
To verify the accuracy of the measured parameters, a new method comparing eld measured with calculated power ows was developed. Sending-end
power ow was measured from accurate recordings and was compared with
power ow calculated using the three estimates of ABCD parameters, sendingend voltages, and receiving-end voltages. Close agreement was found between
eld-measured power ows and calculations using the measured ABCD parameters.
Introduction
1269
1270
R. E. Wilson et al.
1271
Experimental M ethods
V (l n)
Send (kV)
52 .76 + j
311 .4 + j
187 .6 + j
129 .8 j
294 .6 + j
269 .9 j
68 .43 + j
287 .9 j
305 .9
1 .80
244 .4
280 .6
309 .9
158 .3
302 .2
111 .2
Send (A)
84 .29 + j
471 .84 j
354 .6 + j
351 .6 j
73 .5 + j
399 .4 + j
133 .6 + j
666 .6 j
413 .2
22 .10
499 .7
612 .7
316 .6
85 .3
480 .6
453 .8
V (l n)
Recv. (kV)
83 .61 + j
307 .7 j
219 .8 + j
174 .6 j
5 .20 + j
259 .4 j
104 .3 + j
302 .3 j
297 .4
33 .3
215 .9
251 .4
313 .9
178 .7
289 .7
53 .8
1272
R. E. Wilson et al.
All CCVTs were of identical manufacturer and model number. The authors
estimated that the signal delay through the CCVTs would be roughly equal. The
stated accuracy of the CCVTs is 1.2% for voltage magnitudes and 1 for phase
angles. (Burnett et al., 1994) estimated that the phase angle errors introduced by
the instrument transformers would be small at load levels.
The phasor measurement units use GPS to synchronize within one microsecond
digital sampling of currents and voltages within one phasor measurement unit at
dierent substations nearly simultaneously, providing magnitude and phase angles
(phasors) of power system quantities (Slavinsky et al., 1992). An integral GPS receiver recovers Coordinated Universal Time (UTC) and produces a timing pulse
that phase locks the sampling hardware of the instrument. The delay through each
channel of the measuring unit is very stable and well known. An internal algorithm
estimates phasor quantities. Phase angles of input signals are measured to better
than 0 .02 relative to UTC at each substation (60 Hz systems). Because measurements at dierent substations are synchronized to UTC, the phase angle across
a transmission system is found by subtracting receiving values from sending values. Thus, wide-area phase angle measurements can be made to 0 .04 . The phasor
magnitude is accurate to 0 .5% of full scale amplitudes.
Only prefault periods, when the power system is in steady-state and usual
phasor methods can be used, were analyzed. All work in this paper uses positive sequence values. Future work will analyze zero-sequence values. Out of 11
staged faults, usable simultaneous prefault data were recorded for eight events. The
tested transmission system was in two states. For six events, the 190 MVAr Olinda
shunt reactor was in service; for two events the Olinda shunt reactor was out of
service.
During the eld testing, no absolute measurements using a traceable standard were performed because laboratory-grade instrument transformers were not
installed nor were the substations instrument transformers calibrated.
In the rst subsection, the results from a standard transmission line parameter
calculation program are present. The second subsection shows how ABCD parameters were calculated from a carefully constructed EMTP design model. The third
subsection develops a new method using synchronized data to measure the ABCD
parameters. The zone of measurement was identical to the zone of protection of
the transmission line relays. The zone included the series capacitors at Maxwell
Substation, the shunt reactors at Tracy, and for all but Events E9 and E11, the
shunt reactors at Olinda.
During steady-state operation any system element, such as a transmission line,
reactor, or series capacitor, may be characterized by complex-valued two-port parameters, such as the ABCD transmission parameters. As given by (Stevenson,
1982), sending- and receiving-end quantities are related by:
V s = A V r + B I r volts,
I s = C V r + D I r amps.
(1)
1273
To aid in the design of the COTP, a complete model of the transmission system
was developed within the EMTP. From equation (1):
A = V s =V r for I r = 0,
(2)
B = V s =I r for V r = 0 ohm,
C = I s =V r for I r = 0 siemens,
D = I s =I r for V r = 0 .
To evaluate the parameters, the sending and receiving ends of the complete model
were driven by current or voltage sources at one end; the opposite end is open or
short circuited as required by equation (2). The results are shown in column 3 of
Table 2.
Table 2
Transmission parameters
With
reactors
A
B (ohms)
C (S)
D
No
reactors
A
B (ohms)
C (S)
D
Classical
0.9744 @
76.00 @
4.23E-4 @
0.9933 @
0.0996
86.73
89.92
0.012
Classical
0.9744 @
76.00 @
0.0012 @
0.9364 @
0.0996
86.73
90.03
0.21
EMTP-Derived
0.9746 @
75.18 @
4.25E-4 @
0.9932 @
0.102
86.57
89.97
0.01
EMTP-Derived
0.9747 @
75.18 @
0.00116 @
0.9367 @
0.10
86.6
90.11
0.21
Measured
0.967 @
75.86 @
4.52E-4 @
1.004 @
0.05
87.30
91.82
0.11
Measured
0.968 @
74.52 @
0.001 @
0.9400 @
0.12
87.18
90.63
0.30
1274
3.3
R. E. Wilson et al.
M easurement of A BCD Parameters
I r 2 * V s 1 ) = det,
B = ( Vr 2 * V s1
V r 1 * V s 2 ) = det ohm,
(3)
C = ( I r 1 * I s2
I r 2 * I s 1 ) = det S ,
D = ( I s 1 * Vr 2
I s 2 * V r 1 ) = det,
where det = I r 1 * V r 2 I r 2 * V r 1 .
Measured complex sending-end and receiving-end voltages from two events were
used in Cramers Rule to solve for the ABCD parameters. The results are shown
in column 4 of Table 2. Events E6 and E8 were used to calculate the parameters
because these events diered the most in power ow and would be the most linearly independent. The determinant in nonzero (det = 1.012E8 + j 7.025E7 watts),
showing that the events are not linearly dependent. The entire Western Systems
Coordinating Council transmission system itself forced dierent watt and VAr ows
from hour to hour on the COTP system over the 2-day commissioning period and
produced the linear independence needed for solution.
Verication of Results
The question remains of how to determine which set of ABCD parameters is most
accurate. To answer this question, the following method of comparing recorded
power ows with power ows calculated from ABCD parameters was developed.
Accurate values of sending-end voltages and current were available from the recordings made on the phasor measuring unit (in eect, the unit was used as a 16-bit watt
and VAr transducer). Because all channels are sampled within a microsecond of one
another, angle information, as well as magnitude information, can be recovered from
the recordings of the real and imaginary parts of the positive sequence quantities.
Measured sending-end power was calculated by simply multiplying the sendingend (Olinda) voltage times, the sending-end current times, the angle between the
current and voltage times appropriate instrument, and three-phase scaling factors.
The following shows the new technique for the calculation of sending-end power
S s using ABCD two-port parameters and wide-area-synchronized measurements of
the magnitude and angle of the sending- and receiving-end voltages. By denition,
*
S s = P s + j Q s = V s I s MVA.
(4)
*
2
*
A D = B ) + | V s | ( D =B )
*
= V s V r ( 1= B ) + | V s | ( D =B ) .
(5)
1275
The three values of ABCD parameters were then used in equation (5) with measured
sending-end voltage V s and receiving-end voltage V r , to calculate sending-end total
power.
Comparisons between recorded power ows and power ows calculated from
ABCD parameters developed by classical methods are listed in Table 3. The row
labeled Avg lists the average dierence between recorded and calculated power
ow values, while the row labeled Var lists the variance of data. Events E9 and
E11 are separated in Tables 3, 4, and 5 because the 190 MVAr Olinda shunt reactor was out of service during these events. With the reactors in service, the average
dierence was 5 .40 + j 17 .11 MVA, but the variance of the dierences was 6 .75 + j
271 .0 (MVA) 2 . The large MVAr variance was due to heavy megawatt ow (E10)
when the classical parameters overestimated the megavar ow. Without the reactors, the average was 7 .35 j 20 .1 MVA, with variances of 19 .80 + j 3 .24 (MVA) 2 .
Using the EMTP-derived parameters, columns ABCD/ EMTP and RecordedABCD/ EMTP of Table 4 show the results of the comparison. The average dierence was 4.53 + j 21 .85 MVA, with variances of 3 .84 + j 4 .28 (MVA) 2 . The authors
believe the improved megavar estimates were obtained because the EMTP-ABCDs
were a time-domain solution of the COTP system with shunt eects uniformly
distributed along the transmission line. The classical technique uses pi-equivalent
models of the transmission line with shunt-connected eects modeled as a lumped
element at each end of the line.
Columns ABCD/ Meas. and Recorded (ABCD/ Meas.) of Table 5 show the
results. The results from Events E6 and E8 are not signicant because these events
were used to calculate the parameters. Inclusion of these events would have produced a circular argument. The row labeled Avg was calculated using four independent events: E1, E2, E5, and E10. When the calculated ABCD parameters
were compared with four independent measurements, the measured and calculated
Table 3
First comparison of power ows
Recorded
MVA
P1
Classical
ABCD
Q1
P2
Q2
Recorded
classical
P1-P2
Q1-Q2
E1
E2
E5
E6
E8
E10
Avg
Var
395.1
439.7
530.09
654.1
301.3
464.0
12.0
24.0
7.0
57.9
40.3
21.6
391.1
432.6
522.3
649.5
294.0
463.2
9.3
2.2
18.3
31.7
63.4
41.1
4.00
7.10
8.60
4.60
7.30
0.80
5.40
0.00
21.30
26.20
25.30
26.20
23.10
19.50
17.10
271.01
E9
E11
Avg
Var
283.9
730.2
259.5
170.3
281.0
718.4
277.8
192.2
2.90
11.80
7.35
19.80
18.30
21.90
20.10
3.24
1276
R. E. Wilson et al.
Table 4
Second comparison of power ows
Recorded
MVA
ABCD/ EMTP
Event
P1
Q1
P3
E1
E2
E5
E6
E8
E10
Avg
Var
395.1
439.7
530.9
654.1
301.3
464.0
12.0
24.0
7.0
57.9
40.3
21.6
398.9
443.6
532.6
661.1
304.9
471.2
E9
E11
283.9
730.2
259.5
170.3
284.7
726.9
Q3
Recorded
(ABCD/ EMTP)
P1-P3
Q1-Q3
6.0
0.6
16.1
33.8
62.3
1.1
3.80
3.90
1.70
7.00
3.60
7.20
4.53
3.84
18.00
23.40
23.10
24.10
22.00
20.50
21.85
4.28
281.2
194.9
0.80
3.30
21.70
24.60
Table 5
Third comparison of power ows
Recorded
MVA
ABCD/ Meas.
Event
P1
Q1
P4
Q4
E1
E2
E5
E6
E8
E10
Avg
Var
395.1
439.7
530.9
654.1
301.3
464.0
12.0
24.0
7.0
57.9
40.3
21.6
394.2
438.5
527.1
654.1
301.4
465.9
16.4
23.5
7.4
58.3
40.3
24.1
E9
E11
283.9
730.2
259.5
170.3
See
See
text
text
Recorded
(ABCD/ Meas.)
P1-P4
Q1-Q4
0.90
1.20
3.80
0.00
0.10
1.90
1.00
4.07
4.40
0.50
0.40
0.40
0.00
2.50
1.70
3.62
Conclusions
This study has shown that synchronized measurements at opposite ends of a transmission system can be used to actually measure transmission system parameters.
The measured parameters produced improvements in calculations of watt ow and
1277
A cknowledgment
When the rst author was with the University of Wyoming, the work was supported
by Western Area Power Administration, U.S. DOE, under contract numbers AAPO10841-16814 and AA-PO-10593-19197.
References
Burnett, R. O., Butts, M., Cease, T. W., Centeno, V., Michel, G., Murphy, R. J., Phadke,
A., 1994, Synchronized Phasor Measurements of a Power System Event, IEEE Winter
Power Meeting, New York, NY, Paper No. 94-WM-178-4.
Chan, S. M., 1993, Computing Overhead Line Parameters, Computer Applications in
Power, Vol. 6, No. 1, pp. 4345.
Denys, P. H., Counan, C., Hossenlopp, L., Holweck, C., 1992, Measurement of Voltage
Phase for the French Future Defense Plan Against Loss of Synchronism, IEEE Trans.
on Power Delivery, Vol. 7, No. 1, pp. 6269.
Dommel, H., 1985, Overhead Line Parameters From Handbook Formulas and Computer
Programs, IEEE Trans. on PAS, Vol. PAS-104, No. 4, pp. 366370.
Koglin, H., 1994, Estimation of Transmission Line Parameters by Evaluating Fault Data
Records, Discussion of CIGRE Group 34, Paris, pp. 2930.
Liu, H.-W., Lim, S.-L., 1995, Parameter Error Identication and Estimation in Power
System State Estimation, IEEE Trans. on Power Systems, Vol. 10, No. 1, pp. 200
207.
1278
R. E. Wilson et al.
Meliopoulos, A. P. S., Fardanesh, B., 1995, Phasor Measurement ApplicationsTransmission Line Loss Measurement and Parameter Identication; Monitoring for Transformer Unbalance Operation, Fault and Disturbance Analysis & Precise Measurements in Power Systems Conference, Arlington, VA.
Muler, C., MATLAB Users Guide, pp. 3234.
Overholt, P., Rizy, D. T., Mittelstadt, W. A., Litzenberger, W., Krause, P. E., W ilson,
R. E., 1995, The DOE Wide-Area Measurement System, A Demonstration Project
for Monitoring of Power System Dynamics, Fault and Disturbance Analysis & Precise
Measurements in Power Systems Conference, Arlington, VA.
Phadke, A., Thorp, J. S., Karimi, N. J., 1985, Real Time Voltage Measurements for Static
State Estimation, IEEE Trans. on PAS, Vol. PAS-104, No. 11, pp. 30983104.
Phillippot, L., Maun, J. C., 1995, An Application of Synchronized Phasor Measurement
to the Estimation of the Parameters of an Overhead Transmission Line, Fault and
Disturbance Analysis & Precise Measurements in Power Systems Conference, Arlington, VA.
Slavinsky, C., Stanton, S., Esztergalyos, J., Nordstrom, J., Martin, K., Sterlina, P. S.,
Lemme, P. K., Centano, V. A., 1992, Phase Angle Measurement Applications, A New
Transducer that Measures Sate Variables in Real Time, Western Protective Relay
Conference, Spokane, WA.
Slutsker, I. W., Mokhtari, S., Clements, K. A., 1996, Real Time Recursive Estimation
in Energy Management Systems, IEEE Trans. on Power Systems, Vol. 11, No. 3,
pp. 13931399.
Slutsker, I. W ., Mokhtari, S., 1995, Comprehensive Estimation in Power Systems: State,
Topology, and Parameter Estimation, 1995 American Power Conference, Chicago, IL.
Stanton, S. E., Slivinski, C., Martin, K., Nordstrom, J., 1995, Application of Phasor
Measurements and Partial Energy Analysis in Stabilizing Large Disturbances, IEEE
Transactions on Power Systems, Vol. 10, No. 1, pp. 297302.
Stevenson, W . D., 1982, Elements of Power System Analysis, 4th ed., McGraw Hill, pp. 94
107.
Wilson, R. E., Sterlina, P. S., Griess, B. W., 1993, GPS Synchronized Power System
Phase Angle Measurements Recorded During 500 kV Staged Fault Testing, Precise
Measurements in Power Systems Conference, Washington, DC.