An Equivalent Pi Network Model For PDF
An Equivalent Pi Network Model For PDF
equivalent value itself. The model makes true sense when it where R is the measurement error covariance matrix, H is the
comes to transformer connected buses. measurement Jacobian and G is the state estimation gain
It is not always necessary to estimate all the six parameters matrix.
of the equivalent shown in (6), (7) and (8). A transformer is The measurement residuals detect any error only when
plainly used to shift the voltage level and is generally not a there is enough local redundancy. If the measurement is
part of a line having line impedance and line charging critical in nature, the error in it goes unidentified.
admittances. Such lines mostly have only the transformer
B. Identification by Branch Index
reactance in place. If this is the case, all the six terms
mentioned need not be estimated. If this is the case, only the The high normalized residuals need to be linked to the
three reactive components would have the transformer effect, correct branch to ensure correct estimation. The work
and would have to be estimated. On the other hand, if the line pertaining to this paper identifies branches by ranking them in
has no transformer, the π model would not have different order of their respective indices. For a given branch,
leakage admittances on either side thus reducing the number normalized residuals would arise from its branch flows or the
of parameters to be estimated, to four. power injections at the end nodes. The approach used here
The estimation of the equivalent parameters would yield in takes into account, all normalized residuals greater than 3.0,
the same state estimate as the estimation of the individual pertaining to a branch. The maximum among these is the
parameters. However, if the estimation of the individual index of that branch, as shown in (12).
parameter is also required, they can always be obtained from Indexi = max ( rjN )
the equivalent anytime, as shown in the further section on (12)
normal equations method of parameter estimation. i ∈ [1, m ] , j ∈ {rjN > 3.0}
where m is the number of branches and j is corresponding to a
III. PARAMETER ESTIMATION
measurement associated with the branch i. The branch with
A. Detection of Parameter Error the maximum index is chosen to have its parameters
estimated.
The method employs the normalized measurement
Critical Tuples
residuals [2] to detect any parameter error. A parameter error
In case, two branches end up with identical values for their
is analogous to the correlated errors in any of the
indices, the procedure has not been able to clearly identify any
measurements adjacent to the erroneous branch. These
one branch. The parameter error could be in either branch.
measurements are the power flow on the erroneous branch and
This generally happens in case of critical tuples which has
the power injections on the end nodes. Thus, the measurement
been discussed in [8]. The identical value may also be
residuals provide the window to lookout for parameter errors
attributed to the normalized residual of a power injection
as well. The normalized measurement residual, when
common to both the branches. In such cases, counter-checking
calculated to be above a threshold value, generally 3.0 [2],
the normalized residuals of the branch flows so as to identify
indicates an indirect error in the calculation of that
which branch has a higher value may help provided they too
measurement due to the actual error in one or more branch
are not similar.
parameters related to that measurement. This only detects the
error but does not identify the branch. The erroneous branch is C. Estimation by Normal Equations Method
identified by the branch index described in the immediately The normal equations method of augmenting the state
next section. vector clubs the conventional state variables with the
The residual for a measurement zi is given by (9). The suspicious parameters to be estimated together using the
normalized value is given by (10). Newton-Raphson technique, as shown in (13).
∧
ri = zi − hi ( x) (9)
∧
x = ⎡⎣ xV , xθ | x p ⎤⎦ (13)
where x is the estimated state of the prior state estimation, zi is where V, θ and p represent voltage magnitude, its angle and
the vector of given measurement set and hi is the vector of parameter respectively.
calculated measurements. The residuals can be determined One area of concern with this approach is regarding the flat
only after the state estimation process is done. start of the variables generally done before the start of the first
iteration. The flat start leads to Jacobian singularity on
ri account of last column of the Jacobian turning null very often.
ri N = (10) To counter this, the parameters are augmented to the state
Ωii vector after the first iteration.
Ω is the covariance matrix of the measurement residuals,
D. Equivalent and Individual Estimates
given in (11).
The approach in this paper directly gives the final estimates
Ω = R − HG −1 H T (11) eq eq eq
of the equivalent network, given by y , yoi , yoj . This
4
would suffice for the purpose of state estimation to yield Table II gives the ranking of the branch indices for each of
correct states because at each iterative step, the network the three cases. The indices identify the branches distinctly
admittance matrix is updated with the equivalent element and correctly in each of the three cases.
estimates. However, the method is not incapable of updating
the individual parameters. The equivalent parameters are TABLE II
BRANCH INDEX RANKING
nothing but functions of these. Thus, if the equivalent network
estimates are given by k1, k2 and k3, the solution of (14) in 3
variables would yield the individual parameter estimates. Any Branch
Rank of Highest Corresponding
simple non-linear solution technique would yield the required Subjected
Normalized Residual Branch
results. To Error
k1 = ayij
12.5758 13-14
k 2 = a yo + a ( a − 1) yij
2
(14) 13-14 10.7922 9-14
This set of equations can be written twice to solve for g and b 6.63021 7-9
separately if required. 7-9 3.75991 4-7
Transformer Tap 3.75991 7-8
The variable estimated in (14) is a, which may not
necessarily be the tap ratio. If the original ratio was inverted to 6.46194 3-4
convert to 1:a form for developing this equivalent model, the 3-4 4.8282 2-3
true estimate would be 1 / a. 3.78833 2-4
13-14 r 0.06701 1.06701 Table IV gives the ranking of the branch indices for each
of the three cases. The indices identify the branches distinctly
7-9 x 0.11001 1.11001 and correctly in each of the three cases.
3-4 bs 0 0.7
5
x 0.25202 0.55202
5-6
TABLE IV tap 0.932 0.862
BRANCH INDEX RANKING
TABLE VII
STATE ESTIMATION RESULTS
8.69994 29-30
Table VII shows clearly that the state estimation results
29-30 8.69994 27-30
after the estimation of equivalent parameters are the same
as without the error. In addition, when (14) is solved for 8.47777 27-29
these numerous cases given above, the individual
parameters match exactly with the original values 2) Case II: Multiple Errors on Single Branch
mentioned in the tables. This case deals with introduction of error in more than
B. IEEE 30 bus system one component except transformer tap, on a single branch.
Table X shows three different branches having errors on r,
1) Case I: Single Error on Single Branch
x and bs at the same time. But only one branch is
Table VIII shows particulars of three different branches
considered to have the errors.
having a single error, but not simultaneously. They are
three separate cases and are handled individually.
TABLE VIII TABLE X
SINGLE ERROR ON SINGLE BRANCH MULTIPLE ERRORS ON SINGLE BRANCH
Branch
Branch Component Original Error
Component Original Error Subjected
Subjected in Error Value Value
in Error Value Value To Error
To Error
r 0.0472 1.0472
Table IX gives the ranking of the branch indices for bs 0.0170 0.2170
each of the three cases. The indices do not identify the r 0.0169 1.0169
branches distinctly in the case of error in the leakage 6-28 x 0.0599 1.0599
susceptance of branch 29-30. bs 0.0130 0.2130
Identical Normalized Residuals:
The maximum normalized residuals of both branches
Table XI gives the ranking of the branch indices for
27-30 and 29-30 are identical. This is due to the maximum
each of the three cases. The indices identify the branches
residual arising from the common power injection at bus
distinctly and correctly in each of the three cases.
30. Bus 30 is connected to only two buses 27 and 29. So,
the power injection at bus 30 is a measurement adjacent to TABLE XI
both buses 27 and 29. However, on comparing the BRANCH INDEX RANKING
normalized residuals of the branch flows on the two
branches, the branch flow on 29-30 has a normalized Branch
residual greater then 8 whereas that on branch 27-30 has Rank of Highest Corresponding
Subjected
just above 0.5. It is clear from this that the erroneous Normalized Residual Branch
To Error
branch is 29-30.
21.6674 2-5
TABLE IX
BRANCH INDEX RANKING 2-5 16.5677 5-7
12.786 2-4
Branch
Rank of Highest Corresponding 10.6371 6-7
Subjected
Normalized Residual Branch 6-7 7.97217 5-7
To Error
5.53609 2-6
11.9403 2-4 5.42177 6-28
2-4 7.78801 3-4 6-28 3.44432 8-28
7.78801 4-6 3.44432 28-27
20.2966 4-6
4-6 11.7725 3-4 3) Case III: Transformer Tap Error on Single Branch
11.7725 4-12 Table XII shows two different branches having the tap
7
error. This case too pertains to one branch at a time. XIV compares the state estimation results obtained when
there is no error with those obtained with error and
TABLE XII
corrected equivalent parameters.
TRANSFORMER TAP ERROR ON SINGLE BRANCH
TABLE XIII
BRANCH INDEX RANKING
Branch
Component Original Error
Subjected
in Error Value Value Branch Highest
To Error Corresponding
Subjected Normalized
Branch
x 0.2560 0.5560 To Error Residual
4-12
tap 0.932 0.832
5.0062 4-12
x 0.5560 0.2560
6-10 4-12 2.85563 2-4
tap 0.969 0.869
2.85563 3-4
Table XIII gives the ranking of the branch indices for 10.0644 6-10
both the cases. The indices identify the correct branch 6-10 5.64865 6-9
clearly on both occasions. 4.38077 9-10
The branches indicated by the branch index based on
highest normalized residuals are estimated for their
equivalent parameters. Due to space constraints, the state
estimation results for only the case of transformer tap error
in branch 4-12 from Table XII have been shown. Table
TABLE XIV
STATE ESTIMATION RESULTS