Kalman Filter Estimation of Impedance Parameters for Medium Transmission Line

TELKOMNIKA, Vol.16, No.2, April 2018, pp.900~908
ISSN: 1693-6930, accredited A by DIKTI, Decree No: 58/DIKTI/Kep/2013
DOI: 10.12928/TELKOMNIKA.v16i2.9026  900
Received November 16, 2017; Revised February 2, 2018; Accepted February 22, 2018
Kalman Filter Estimation of Impedance Parameters for
Medium Transmission Line
Siti Nur Aishah Mohd Amin*, Hamzah Ahmad, Mohd Rusllim Mohamed, Mohd Mawardi
Saari, Omar Aliman
Faculty of Electrical & Electronics Engineering,Universiti Malaysia Pahang,Pekan Campus,26600 Pekan,
Pahang Darul Makmur, 09-4246024, Malaysia
*Corresponding author, e-mail: aishahamin_92@yahoo.com
Abstract
Accurate knowledge of impedance parameters in transmission line helps to improve the system
efficiency and performance.Nowadays,the estimation of impedance parameters in transmission line has
become possible with the availability of computational method. This paper aims to develop Kalman filter
model by using Matlab simulink to estimate accurate values of resistance (R), reactance (X), and
susceptance (B) in medium transmission line. The accuracy of the parameters can be improved by
reducing the unknown errors in the system. To demonstrate the effectiveness of the Kalman filter method,
a case study of simulated medium transmission line is presented and comparison between Kalman Filter
(KF) and Linear Least Square (LLS) method is also considered to evaluate their performances.
Keywords:RXB, kalman filter,matlab simulink,accuracy,transmission line
Copyright © 2018 Universitas Ahmad Dahlan. All rights reserved.
1. Introduction
An electric power system consists of three principal divisions which are the generating
station, transmission station and distribution station. Generally these three stations is a network
of conductors between the power stations and the consumer in order to utilize the power
generated. The power transmission is one of the major components in the electrical power
system [1]. All transmission lines in a power system exhibit the electrical properties of resistance
(R), reactance (X) and susceptance (B). It is important to consider all the R, X and B parameters
in the design and operation of transmission line which these parameters will indicate the power
losses and efficiency of transmission line [2]-[3]. By reducing losses and improving the
efficiency, values of voltage profile and power factor in a system will be corrected, hence
following the standard of IEEE an IEC system. The important consideration of R, X and B in the
design and operation in transmission line are the determination of line losses, voltage drop and
efficiency of transmission line. According to the previous paper, R, X and B are greatly
influenced the transmission line performance in a system [5].
The frequency involve in this three constant parameter indicates the performance of a
system. When considering high frequency, the magnetic field generated drive the currents to the
outer edge of the conductor that carries them. As a result, the higher the frequency, the thinner
the layer of metal available to carry the current and the higher the effective resistance of the
cable [6]. In fact, the voltage drop in the line depends on the values of above three lines
constant [1]. The resistance involved in the transmission line is the most important cause of
power losses in the line and it determine the transmission line efficiency. Since there are many
research has been done on R, L and C parameters value [7–9] in medium and long
transmission line, so this paper will focuses on R, X an B parameters value medium
transmission line. The difference of this paper to other papers is that the analysis includes the
different noises conditions and it is applied to the real system which is in East Midland to
analyse the performance of the system.
Power system state estimation is a main tool of energy management system. State
estimation was first applied to power systems by Schweppe and Wildes in the late 1960s
[4],[10]. Since then, several applications and techniques of KF method have been one of the
famous methods in power system estimation. It has been used for damping and frequency [1],
modal changes [11], state estimation [12], power factor detection and etc. However all of these

TELKOMNIKA ISSN: 1693-6930 
Kalman Filter Estimation of Impedance Parameters for Medium… (Siti Nur Aishah Mohd Amin)
901
researches did not highlighted important indicators that give great effects in the system which
are R, X and B parameters. Theoretically, values of R, X and B in long transmission line are
supposedly higher than short transmission line. But due to some condition, the values of R, X
and B can be opposed the actual properties such as type of cable, size of cable, length of cable,
number of winding and etc. All of these properties will indicate the values of R, X and B in the
transmission line.
Traditionally, all the impedance parameters involved were calculated manually by using
handbook formulas based or through fault record analysis. Accounting, for any system, can be a
complex undertaking. A manual accounting system requires an understanding in the accounting
process in a way that may be unnecessary with a computerized accounting system. This can be
an advantage or a disadvantage, depending on the person doing the bookkeeping, usually a
professionally trained is needed to ensure that accounting in the system is done properly.
Therefore, with the availability of computational method, it is possible to calculate and measure
the impedance parameters by using computational method such as KF and LLS method.
Previous paper has shown the estimation of R, X and B by using LLS method [5]. Thus, in this
paper, with the same information used the estimation of R, X and B value is carried out by using
KF method. Both of these methods are used for parameters estimation purposes in a system
with the different in accuracy of estimation.
KF works on prediction-correction model applied for linear and time-variant or time-
invariant systems [13]. The prediction model involves the actual system which is predictable
with the process noise in the system [13]. While LLS method is based on weighted least
squares in which past values taken in account for determining the future value. The weights are
updated recursively based on memory. LLS method will only use the measurements to update
the position estimate every second but, KF will provide additional information and yield a more
accurate estimation. Kalman filter is a powerful tool especially when it comes to control noise in
a system. Thus, in this paper the development of KF model by using Matlab simulink has been
done. The mathematical model and its standard parameters are also used. The formula for the
KF method estimation will be examined in the next section.
The rest of this paper is structured as follows. In section 2, the transmission line model
is presented. In section 3, KF method based impedance estimation is discussed while LLS
method based impedance estimation in section 4. Simulink model of KF method is shown in
section 5. System model analysis and discussion are made in section 6. Section 7 concludes
this paper.
2. Transmission Line Model
A transmission line can be defined as a medium length line if its length is more than 50
km and the voltage is comparatively higher than 20kV. In this case, the resistance, inductive
reactance and shunt capacitor are taken into account. Balanced conditions of the line can be
represented by the equivalent circuit of a single phase with resistance (R), and inductance (L) in
series (series impedance) as well as shunt capacitor [2]. The nominal pi circuit of standard
model for medium transmission line is shown in Figure 1 [8]. Based on this figure, phasor values
and represents the sending and receiving end voltages and current respectively.
Figure 1. Nominal pi circuit diagram for medium length of transmission line [8]

 ISSN: 1693-6930
TELKOMNIKA Vol. 16, No. 2, April 2018 : 900 – 908
902
For this circuit, the equations are derived based on Kirchhoff’s voltage and current law [14]. The
following equations are:
a. Kirchhoff’s voltage law
( ) (1)
b. Kirchhoff’s current law
( ) (2)
Where , , and . X is the inductive reactance and
B is the capacitive susceptance. Formula for impedance, Z and admittance, Y can be written as
below by substituting equation (2) into (1)
(3)
(4)
The values of Z and Y can be calculated. R, X, G and B are obtained from real and imaginary
parts of Z and Y respectively.
3. Kalman filter based impedance estimation
a. R, X and B model
For a linear system, the process model of Kalman filter is described as
(5)
(6)
Where is the estimation value of X, is the estimate value on the previous state, is the
control input, and is noise covariance matrix, A, B and C will be numerical constant. The
estimation value of is represented by R, X and B as
[ ] (7)
Where R is resistance, X is reactance and B is susceptance. In this study, current, A is used as
a control input, .
[ ] (8)
The estimated value on the previous state, is represented by
[ ] (9)
b. Kalman filter algorithm
Kalman filter is used to provide an estimation of R, X and B values in the system.
Kalman filter recursively computes estimation for a state according to the process and
observation model. The stages of Kalman filter algorithm are as follows:
1) Prediction (time update) to estimate priori estimation of state and its error covariance matrix:

903
(10)
(11)
Where is the prior error covariance. This prior value is used in the measurement update
equations.
2) Correction (measurement update) to provide correction based on the measurement to
yield the posteriori state estimate and its error covariance:
( ) (12)
( ) (13)
( ) (14)
where is Kalman gain.
4. Linear Least Square Method Based Impedance Estimation
The proposed method was applied to improve the accuracy of the calculated Z and Y
values. For comparison, an existing method which is linear squares method was applied to
obtain the parameter estimation [8]. The transmission line was modelled as a general two port
network, in which voltages and currents are related to.
(15)
(16)
Where A, B, C and D are constant. The real and imaginary parts of A and B were computed
through unbiased linear LS estimation [5]. By assuming a pi circuit in Figure 1, the constant A,
B, C and D inside the two port network can be expressed in terms of impedance, Z and
admittance, Y as shown below.
(17)
(18)
( ) (19)
(20)
Z and Y are calculated from LS estimation of A and B by using (3) and (4).
(21)
( )
(22)
5. Simulink Model of Kalman Filter
Matlab simulink is a practical tool for modelling the real scenario in computational
method. Figure 2 shows the general development of Kalman filter model by using Matlab
simulink.

 ISSN: 1693-6930
904
Figure 2. Kalman filter model
6. System Model Analysis and Discussion
This subsection investigates the results of having different noise conditions to the
estimation of medium transmission line. This paper takes into account system noises of 0.1 and
50 to simulate the unknown parameters which can effects the estimation results such as size,
length and type of cables. The results obtained are described as follows.
A single phase of the 400 kV, 102 km long transmission line data located between
substations Grendon and Staythorpe, East Midlands, England was used and simulated by using
Matlab simulink [5],[15]. This section analyse how Kalman Filter works to estimate those value
efficiently.
Simulation was conducted for the equivalent circuit of medium transmission line and
tuning process also take place in order to gain converge graph of each parameters of the circuit.
Figure 3, Figure 4 and Figure 5 below shows the real and estimation values of R, X and B over
a period. The real values of R, X and B are 2.96 Ω, 32.40 Ω and 3.69 x S while the values
of estimated are 2.94 Ω, 30.92 Ω and 3.50 x respectively with noise 0.1. There are small
differences between the real and new estimation values of R, X and B. Both of the real and
estimated parameters are seen overlapping with each other due to the approximate parameter
values. The result shows the higher consistency involved in R, X and B value due to the
approximate results between real and estimated values.
Figure 3. Graph of real and estimated value of resistance, R

905
Figure 4. Graph of real and estimated value of reactance, X
Figure 5. Graph of real and estimated value of susceptance, B
Table 1 defines the system parameters of the simulink model with the small value of
noise which is 0.1. It shows the accuracy of R, X and B in the system more precise with the
small value of noise which are 99.32 %, 95.43 % and 94.85 % respectively. Thus, it is shown
that higher estimation accuracy will be achieved with the decreasing values of the noise.
Table 1. Parameter values
Real
values
Estimation
values
Accuracy
Resistance, R 2.96 Ω 2.94 Ω 99.32
Reactance, X 32.4 Ω 30.92 Ω 95.43
Susceptance,
B
3.69 x
S
3.50 x S 94.85
Figure 6, Figure 7 and Figure 8 shows the real and estimated values of R, X and B in
the system. The real values of R, X and B are 2.96 Ω, 32.40 Ω and 3.69 x S while the
estimated values R, X and B are 2.66 Ω, 33.02 Ω and 2.84 x S respectively with higher
system noise of 50. Both of the calculated and estimated parameters are seen less consistency
and the results are slightly different due to the greater noise in the system.

 ISSN: 1693-6930
906
Figure 6. Graph of real and estimated value of resistance, R
Figure 7. Graph of real and estimated value of reactance, X
Figure 8. Graph of real and estimated value of susceptance, B
Table 2 defines the system parameters of the simulink model with the greater value of
noise which is 50. It shows the accuracy of R, X and B in the system are inaccurate with the

907
greater value of noise which are 89.86 %, 98.08 % and 76.96 % respectively. Thus, it is proved
that the accuracy in the system will decrease with the increasing noise.
Table 2. Parameter values
Real
values
Estimation
values
Accuracy
Resistance, R 2.96 Ω 2.66 Ω 89.86
Reactance, X 32.4 Ω 33.02 Ω 98.08
Susceptance,
B
3.69 x
S
2.84 x S 76.96
According to these table and graphs presented previously, the value of noise affected
the accuracy of R, X and B estimation in the system. The higher the value of noises or
disturbances in the system, the performance will degrades. These results of R, X and B
estimation indicate the possibility for parameters estimation even in different places. Noted that,
it is important to consider that the values of R, X and B affected the determination of system
efficiency in the transmission line.
Comparison of parameter errors between Kalman filter and Least Square method also
has been done [5]. The line parameter errors of R, X and B by using Kalman filter are 0.02 Ω,
1.48 Ω and 0.16 S while the line parameter errors of R, X and B by using Linear Square method
are 11.9 Ω, 0.212 Ω and 38.7 S. Generally, from the results presented in Table 3 below, Kalman
filter give a higher accuracy compared to the Least Square method. From these two results of
accuracy, it can be proven that Kalman filter is a better method that can be used for estimation.
Table 3. Parameter errors
R X B
KF 0.02 1.48 0.16
LS 11.9 0.212 38.7
KF: Kalma filter method
LS: Least Squares Method
Remark that, it is important to consider of R, X and B as they affected the determination
of power transfer, voltage drop and efficiency in the transmission line. The voltage drop in the
line depends upon the values of above three lines constant. Similarly, the resistance of
transmission line is the most important cause of power losses in the line and determine the
transmission efficiency.
7. Conclusion
This paper presented the analysis of R, X and B estimation in a system. The analysis
has proven that Kalman filter is one of the optimal methods that can be used for parameter
estimation in transmission line. From this research, it has been revealed that noise in the
system indicates the performance of medium transmission line.
Acknowledgement
Author would like to thank University Malaysia Pahang and Ministry of Higher Education
for the supports in conducting the research under RDU160379 and RDU160145.
References
[1] P Korba, M Larsson, C Rehtanz. Detection of Oscillations in Power Systems using Kalman Filtering
Techniques. Control Applictions 2003. Istanbul, Turkey. 2003: 183-188.
[2] S Gajare, AK Pradhan, V Terzija. A Method for Accurate Parameter Estimation of Series
Compensated Transmission Lines Using Synchronized Data. IEEE Transactions on Power Systems.
2017; 32(6): 4843-4850.

 ISSN: 1693-6930
908
[3] MV Venkatasubramanian. Coordination ofTransmission Line Transfer Capabilities. Washington State
University. 2002.
[4] D Alazard. Introduction to Kalman Filtering. France. ISAE-SUPAERO (FRANCE). 2005: 1-74.
[5] D Ritzmann, PS Wright, W Holderbaum, B Potter. A Method for Accurate Transmission Line
Impedance Parameter Estimation. IEEE Transactions on Instrumentation Measurement.2016;65(10):
2204-2213.
[6] B Curran, I Ndip, K Lang. A Comparison of Typical Surface Finishes on the High Frequency
Performances of Transmission Lines in PCBs. Signal and Power Integrity. Baveno, Italy. 2017: 1-3.
[7] S Poornima,R Sharma.Estimation ofLine Parameters ofan IEEE 14 Bus System.RSIS International.
2014; 1: 14-16.
[8] D Shi, DJ Tylavsky, N Logic. Identification of short transmission-line parameters from synchrophasor
measurements. Power Symposium, 2008. Calgary, Canada. 2008: 1-8, 2008.
[9] C Li, Y Zhang, H Zhang, Q Wu. Measurement-Based Transmission Line Parameter Estimation with
Adaptive Data Selection Scheme. IEEE Transactions on Smart Grid. 2017: (99): 1-1.
[10] A Saikia, RK Mehta. Power system static state estimation using Kalman filter algorithm . EDP
Sciences. 2016; 7: 1-7.
[11] RA Wiltshire, G Ledwich, S Member, PO Shea. A Kalman Filtering Approach to Rapidly Detecting
Modal Changes in Power Systems. IEEE Transactions on Power System. 2007; 22(4): 1698-1706.
[12] K Vadirajacharya, OS Mankar. Kalman Filter Analysis in Dynamic State of Power System.
International Research Journal of Engineering and Technology. 2015; 2(6): 1-4.
[13] G Bishop, G Welch. An Introduction to the Kalman Filter. Los Angeles: ACM. 2001: 1-81.
[14] G Sivanagaraju, S Chakrabarti, SC Srivastava. Uncertainty in transmission line parameters:
Estimation and impact on line current differential protection. IEEE Transactions on Instrumentation
Measurement. 2014; 63(6): 1496–1504.
[15] R Smith. Electricity Ten Year Statement 2015. United Kingdom. National Grid. 2015: 2-13.

Kalman Filter Estimation of Impedance Parameters for Medium Transmission Line

More Related Content

Similar to Kalman Filter Estimation of Impedance Parameters for Medium Transmission Line (20)

More from TELKOMNIKA JOURNAL (20)

Recently uploaded (20)

Kalman Filter Estimation of Impedance Parameters for Medium Transmission Line