Next Article in Journal / Special Issue
On the Efficiency of a Two-Power-Level Flywheel-Based All-Electric Driveline
Previous Article in Journal
Deacidification of Pistacia chinensis Oil as a Promising Non-Edible Feedstock for Biodiesel Production in China
Previous Article in Special Issue
The Impact of Power Switching Devices on the Thermal Performance of a 10 MW Wind Power NPC Converter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Steady-State and Faulty Regimes of Overhead Lines by Means of Multiconductor Cell Analysis (MCA)

by
Roberto Benato
*,
Sebastian Dambone Sessa
and
Fabio Guglielmi
Department of Industrial Engineering, University of Padova, Via Gradenigo, 6/A, Padova 35131, Italy
*
Author to whom correspondence should be addressed.
Energies 2012, 5(8), 2771-2793; https://doi.org/10.3390/en5082771
Submission received: 15 June 2012 / Revised: 24 July 2012 / Accepted: 26 July 2012 / Published: 31 July 2012
(This article belongs to the Special Issue Electrical Power and Energy Systems)

Abstract

:
Single-phase positive sequence modelling is often used in power systems when power flows and short circuit analysis are assessed. Of course, the use of single-phase positive sequence modelling assumes purely three-phase configurations and perfectly symmetrical ones so that single-phase modelling considers that all the phase conductors behave in the same way. When considering the physical reality of power networks, this assumption can be questionable and the behaviors of all the system conductors including the passive ones (earth wires for overhead lines, metallic screens and armours for cables and enclosures for gas insulated lines) is completely unknown. Therefore, the present multiconductor cell analysis (MCA) becomes necessary, since it allows one to achieve great precision results on the regimes of both phase conductors and passive conductors. MCA offers a powerful tool in order to validate (or less) approximated and simplified computation methods. In particular, for single and double circuit overhead lines (OHLs), the current phasors induced in the earth wires and the ground return current alongside the line can be directly computed by MCA in steady state and faulty regimes. It is worth noting that, for faulty regimes, MCA allows also evaluating the approximation degree and validity field of screening factors k.

1. Introduction

Multiconductor cell analysis (MCA) takes advantage of the extensive use of the admittance matrices in order to consider any asymmetrical system with all the active and passive (earth wires for OHLs, screens and armors for underground insulated cables (UGCs), enclosures for gas insulated lines (GILs), keraunic wires, buried conductors, rails, feeders, overhead contact line with catenary suspension for AC high speed railway systems) conductors. It allows knowing directly the electrical behaviors of all the electrical quantities alongside an OHL. In the following, the steady state regime and the faulty regime of OHLs are considered. Furthermore a comparison between MCA and the commercial software EMTP-RV results is presented, which has given a very good agreement.

2. Application of MCA to Overhead Lines

MCA has been already applied to a great number of asymmetrical systems: cross-bonded underground insulated cables [1,2], shunt compensated ones, distribution line carrier, solid bonded or multiple-point bonded gas insulated lines [3], Milliken type EHV insulated cables. In this paper a thorough exposition of the application of MCA to overhead lines is presented. It allows considering:
(1)
Single or double circuit in electrical parallel with any type of phasing (low reactance or super bundle ones);
(2)
Any tower type, with different earth resistivities along the line and different tower resistances;
(3)
The steady-state regime of an OHL supplying a balanced (or not) three-phase load;
(4)
Faulty regimes (line-to-earth, line-to-line, line-to-line-to-earth, three phase short-circuit), with two end substations each characterized by its own fault levels.
For the sake of a clear explanation, let us consider a single-circuit OHL composed of five conductors (n = 5; see Figure 1) parallel to themselves and to the ground: 1, 2, 3 are the phase conductors (in Figure 1, bundle ones) and 4, 5 are the earth wires.
Figure 1. Multiconductor representation of OHL supplying a static balanced load; SEND—sending-end; REC—receiving-end.
Figure 1. Multiconductor representation of OHL supplying a static balanced load; SEND—sending-end; REC—receiving-end.
Energies 05 02771 g001
This multiconductor line may be represented as a cascade of N elementary cells of length Δ (e.g., Δ = span length), modelled by a lumped π circuit (see Figure 2).
Figure 2. Cascade of N elementary cells for OHL modelling.
Figure 2. Cascade of N elementary cells for OHL modelling.
Energies 05 02771 g002
Being Δ sufficiently small, it is possible to lump distributed shunt admittances at the ends of the cell (YTS and YTR are the transverse blocks), thus allowing the separate study of longitudinal elements. In order to form ZL, it is necessary to compute the self and mutual longitudinal impedances, which account for the earth return currents. They can be obtained by applying the simplified Carson-Clem formulae [4] or the complete Carson formulae [5] with Carson’s corrections in terms of series [6]. It has been widely demonstrated by one of the authors that at power frequency (50 or 60 Hz) the Carson-Clem formulae are very precise and consequently they are a powerful tool so giving the same results of complete Carson ones (Section 5).
In order to give the above mentioned sentence a theoretical justification, the simplified Carson-Clem formulae can be licitly applied if the spacing between conductors di,j ≤ 0.135De when D e = 658 ρ / f [m] “Carson Depth”; ρ = electrical resistivity of the soil [Ωm]; f = power frequency (50–60 Hz). By supposing f = 50 Hz and a low earth resistivity e.g., 10 Ωm, it yields:
D e = 658 ρ f = 658 10 50 = 294 [ m ]
The spacings between phases in OHL are never greater than 40 m, therefore it is rather unlikely to have the condition di,j > 0.135De 40 m even with very low resistive earth; hence the Carson-Clem formulae can be always used. With great mutual spacings between phases belonging to different circuits in the same corridor, the complete Carson formulae must be used. The self and mutual impedances by means of Carson-Clem formulae are given by:
z ¯ i , i = r i + π 2 10 4 f + j 4 π 10 4 f ln ( 2 D e d i ) [ Ω/km ]
z ¯ i , j = π 2 10 4 f + j 4 π 10 4 f ln ( D e d i , j ) [ Ω/km ]
where ri = kilometric resistance of conductor i [Ω/km] (if bundle conductor with n subconductors per phase r i = r i _ s u b / n ); di = 2·GMR = k″d [m] (GMR = Geometrical Mean Radius) where k′′ are listed in tables depending upon the composition of the conductor and d [m] is the diameter of the conductor (if bundle conductor:
d i ´ = 2 n k d s u b 2 R c i r c l e n 1 n
where Rcircle [m] is the radius of the circle on which the subconductors lay and dsub [m] is the diameter of the subconductors); di,j = distance between conductor i and j [m]. Once zi,i and zi,j have been computed, it is possible to form the matrix ZL (5 × 5) and to characterize, by means of Equation (4) the steady state regime of longitudinal block L of Figure 2 (where all voltage and current vectors are shown):
u ¯ S u ¯ R = Z ¯ L i ¯ S L
and by considering the evident relation Equation (5):
i ¯ S L = i ¯ R L
it yields:
Z ¯ L 1 u ¯ S Z ¯ L 1 u ¯ R = i ¯ S L
Z ¯ L 1 u ¯ S + Z ¯ L 1 u ¯ R = i ¯ R L
The Equations (6,7) can be written in an elegant matrix form as:
Energies 05 02771 i001
The matrix YLΔ (10 × 10) is the admittance matrix regarding block L circuit formed by five longitudinal links.
The shunt capacitive matrix computation can be obtained from Maxwell potential coefficient matrix. Its elements are easy to compute from the tower geometry and from the conductor radii [7]. It is worth remembering that the average heights of conductors can be used by means of Equation (9):
h a v e r a g e = h t o w e r 2 3 s a g [ m ]
where htower is the conductor height at tower.
Then, as already mentioned, the shunt capacitive and conductive links are divided into two halves and lumped at the ends of the elementary cell. The capacitive susceptances (of half cell) between each conductor are considered in the (5 × 5) matrix YTS = YTR. The vectors of the shunt currents at both ends of the cell are (see Figure 2):
Energies 05 02771 i002
Hence the steady state regime of the elementary cell is represented by YΔ which is the superimposition of YLΔ and YTΔ:
Energies 05 02771 i003

2.1. Steady State Regime of Single Circuit EHV OHL

In order to achieve a complete analysis of an overhead line it is necessary to obtain the admittance matrix Y of the whole multiconductor system shown in Figure 3 (where Sg, S0,…. SL indicate the sections or ports). This involves all the elements of the line which can be represented by their admittance matrix: (by starting from the left side) the equivalent supply network matrix Ypha (3 × 3) (see Appendix I), the load matrix YL (3 × 3) (see Appendix II), the tower matrix YTW (2 × 2) (see Appendix III), and, of course, the already mentioned elementary cell matrices YΔi. In short-circuit occurrence (Section 4), the equivalent substation impedance YSUB (3 × 3) and fault matrix Yfault (see Appendix IV) must be considered.
On the basis of the equation which leads to the general formulation of the current vectors injected in the different ports of the system (Sg, S0,…. SL), it can be ascertained that the admittance matrix Y is formed by partial overlapping ("tile style") of the above mentioned matrices as shown in Figure 3b. The resulting admittance matrix Y can be of large dimensions (depending upon the ratio OHL length/Δl) but it can be easily managed because of it is structurally sparse. Anyway, the use of the admittance frame of reference allows introducing easily local modifications of the system structure without reconstructing its entire matrix.
Figure 3. (a) Single circuit overhead line; (b) Structure of the matrix Y.
Figure 3. (a) Single circuit overhead line; (b) Structure of the matrix Y.
Energies 05 02771 g003aEnergies 05 02771 g003b
If the modification is between two sections, it is sufficient to substitute the matrix representing the modification for that existing before. If the modification is derived from a section, the matrix representing the modification must be superimposed in that section (e.g. the tower presence anywhere along the OHL). In fact, they can be represented by their admittance matrices which must be superimposed in the right location of matrix Y so to give a new total matrix YTOT.

2.2. Steady State Matrix Solution of Multiconductor System

For the system of Figure 3b, the general equation i = YTOT u can be partitioned as in Figure 4. It is worth noting that only ipha is a non zero current vector whereas ix ≡ 0 so that it yields:
i ¯ x 0 = Y ¯ 3 u ¯ p h a + Y ¯ 4 u ¯ x
and thus:
u ¯ x = Y ¯ 4 1 Y ¯ 3 u ¯ p h a
where Y4 is not singular and upha is fixed. So by knowing the subvector ux the steady state regime of any cell, and hence of the whole multiconductor system, is completely available without iterations.
Figure 4. Partitioned matrix of i = YTOTu.
Figure 4. Partitioned matrix of i = YTOTu.
Energies 05 02771 g004
The inversion of Y4 in Equation (13) could give computational drawbacks, but even if Y4 is of great dimensions (i.e., with great number of N cells), its inversion can benefit from the partition in submatrices [8] and the use of parallel computing.

3. Application to HV OHL on Steady State Regime

The proposed multiconductor algorithms are applied to a 275 kV 2 × 400 mm2 single circuit aluminum conductor steel reinforced (ACSR) OHL (Figure 5), whose geometrical and electrical characteristics are reported in Table 1.
Figure 5. Single circuit OHL tower arrangement.
Figure 5. Single circuit OHL tower arrangement.
Energies 05 02771 g005
Table 1. OHL geometrical and electrical characteristics.
Table 1. OHL geometrical and electrical characteristics.
Overhead lineSymbolUnitValue
Voltage level-kV275
Line length-km55.632
Span length-km0.366
Earth resistivity (assumed equal along the line)-Ω m100
Symmetrical loading at receiving end-MW + j Mvar508.8 + j 103.32
Substation earthing resistance-Ω0.1
Tower earthing resistance (assumed equal for all the tower)-Ω10
Conductor diameter (bundle conductor)-mm2 sub-cond. ACSR Lynx Ф = 19.53 mm
Per unit length resistance at 75° operation (50 Hz)RmΩ/km96.10
Per unit length series inductancelmH/km1.1051
Per unit length shunt leakance (50 Hz)gavnS/km44.5
Per unit length capacitancecμF/km0.010215
Ampacity referred to winter rating [9]IaA545 × 2
Winter rating [9]- MVA519.2
By means of the procedures of Section 2.2 it is possible to completely know the steady state behaviors of the electrical quantities along the line. Figure 6 shows the phase voltage behaviors along the line. There is a voltage drop which is different from phase to phase due to the asymmetrical nature of the un-transposed OHL.
Figure 6. Line-to-earth voltage magnitudes along the OHL.
Figure 6. Line-to-earth voltage magnitudes along the OHL.
Energies 05 02771 g006
The phase current magnitudes are almost constant along the line ( 1090 A) (see Figure 7) due to the very low capacitive susceptance of OHL. The earth wire voltage along the line is almost null (see Figure 8). This is a confirmation that it is licit to use the elimination technique of grounded conductors.
It is meaningful to observe that, even in steady state regime, both the induced currents in the earth wires and the ground return current are not negligible. The earth wire induced current magnitudes along the line have an almost constant value of 96 A (see Figure 9).
Figure 7. Phase current magnitudes along the OHL.
Figure 7. Phase current magnitudes along the OHL.
Energies 05 02771 g007
Figure 8. Earth wire voltage magnitudes along the OHL.
Figure 8. Earth wire voltage magnitudes along the OHL.
Energies 05 02771 g008
Figure 9. Earth wire current magnitudes (left axis) and angles (right axis) along the OHL.
Figure 9. Earth wire current magnitudes (left axis) and angles (right axis) along the OHL.
Energies 05 02771 g009
One of the most important features of MCA is that it allows directly achieving the ground return current. Since MCA can compute all the conductor current phasors in each section, IGR is given by the 1st Kirchhoff’s Law:
I ¯ G R + x = 1 n I ¯ x 0 I ¯ G R = x = 1 n I ¯ x
With the constant loading of Table 1 (S = 508.8 MW + j103.32 MVar), |IGR |= 25 A (see Figure 10). This current could be responsible for ac metallic corrosion [10] and induced voltages on parallel metallic pipes buried in the neighborhood of OHL.
Figure 10. Ground return current magnitude and angle along the OHL.
Figure 10. Ground return current magnitude and angle along the OHL.
Energies 05 02771 g010
Another important result from MCA is the possibility to compute the whole power losses which are 19.269 MW. In order to roughly validate this value, it is possible to consider a constant current in the phase conductor equal to 1090 A so that the three-phase power losses result:
P J _ P H = 3 r L I 2 = 3 0.0961 55.632 1090 2 = 19.056 MW
To these phase conductor power losses, the insulator and corona losses must be added so that it has:
P i n s + c o r = g L U 2 = 0.18722 MW
The power losses in the earth wires are not negligible and can be computed (by considering that the induced currents are 96 A constant along the line as in Figure 9) by means of:
P E W = 2 r E W L I E W 2 = 2 0.1575 55.632 96 2 = 0.1615 MW
The whole power losses are given by the sum of the three abovementioned components i.e.:
P = P J _ P H + P i n s + c o r + P E W = 19.4047 MW
which is in very good agreement with the power losses computed by MCA.

Comparison between the MCA and EMTP-RV

In this section some of the above presented MCA results are compared with the EMTP-RV results for the OHL of Section 3. The EMTP-RV electrical layout (the first and the last of four pages) is shown in Figure 11.
Figure 11. EMTP-RV single circuit OHL electrical layout (pages 1/4 and 4/4).
Figure 11. EMTP-RV single circuit OHL electrical layout (pages 1/4 and 4/4).
Energies 05 02771 g011
Figure 12 shows the phase voltage magnitude behavior along the single circuit OHL by means of MCA and EMTP-RV. The per cent voltage drops are reported in Table 2 where it infers that the maximum difference between the MCA and EMTP-RV is equal to 5.28%. The cpu-time in the simulation with MCA is 0.2087 s and with EMTP-RV is 0.3125 s (INTEL CORE 2 QUAD [email protected] GHz, 3.25 GB RAM).
Figure 12. Voltage magnitude comparison between MCA and EMTP-RV along the OHL.
Figure 12. Voltage magnitude comparison between MCA and EMTP-RV along the OHL.
Energies 05 02771 g012
Table 2. Single circuit OHL voltage drop comparison.
Table 2. Single circuit OHL voltage drop comparison.
EliminateVoltage drop [%]
Phase 1-RPhase 2-SPhase 3-T
MCA7.4%6.1%5.49%
EMTP-RV7.4%6.26%5.2%
Percentage difference0.0%2.62%5.28%
In Table 2 “percentage difference” means:
100 M C A E M T P R V M C A
In Figure 13 the steady state |IGR| for OHL of Section 3 is shown by using the MCA and EMTP-RV results. The maximum current difference is 1 A (4%).
Figure 13. Steady state ground return current magnitude comparison between the MCA and EMTP-RV.
Figure 13. Steady state ground return current magnitude comparison between the MCA and EMTP-RV.
Energies 05 02771 g013

4. Short-Circuit Analysis of Double-Circuit Overhead Line 400kV

In this section a double circuit OHL 400 kV 2 × 400 mm2 ACSR, whose electrical and geometrical characteristics are reported in Table 3, has been studied on faulty regime by means of MCA (see Appendix IV). The tower arrangement is reported in Figure 14.
Table 3. Data assumed in the MCA for double-circuit OHL.
Table 3. Data assumed in the MCA for double-circuit OHL.
Overhead lineSymbolUnitValue
Voltage level-kV400
Line length-km135.42
Span length-km0.366
Earth resistivity (assumed equal along the line)-Ω m100
Symmetrical loading at receiving end-MW + j Mvar2 × (2172.7 + j 441.2)
Substation earthing resistance-Ω0.1
Tower earthing resistance(equal for all the towers)-Ω10
Conductor diameter (bundle conductor)-mm24 sub-cond. ACSR Zebra Ф = 28.62 mm
Per unit length resistance at 25° (50 Hz)rmΩ/km8.25
Per unit length series inductancelmH/km0.4218
Per unit length shunt leakance (50 Hz)gavnS/km33.56
Per unit length capacitancecμF/km0.0266521
Ampacity referred to winter rating [11]IaA800 × 4
Winter rating [11] of double-circuitMVA-2 × 2217
Figure 14. Tower characteristics of 400 kV double-circuit OHL in electric parallel with low reactance phasing.
Figure 14. Tower characteristics of 400 kV double-circuit OHL in electric parallel with low reactance phasing.
Energies 05 02771 g014
Let us suppose the fault occurring at midline (i.e., at 67.71 km from both ends) between one phase (e.g., 1R phase) and earth wire (line-to-earth short-circuit). It is worth remembering that in order to characterize the end substations it is necessary to know the fault levels (three-phase x1, x2 and single-phase y1, y2 sub-transient short circuit currents usually given by TSO) as in Figure 15.
Figure 15. Multiconductor line in short-circuit condition supplied by end substations.
Figure 15. Multiconductor line in short-circuit condition supplied by end substations.
Energies 05 02771 g015
Once the fault levels are known, the equivalent sequence impedances and phase impedances of each substation can be computed by means of the procedure in Appendix I, in order to create the equivalent supply network matrix. The short-circuit between a phase and the earth wire can be considered by means of the matrix Yfault reported in Appendix IV. In accordance with IEC 60909-0 [12], a voltage factor must be applied to the nominal line-to-earth voltage. Since maximum currents in HV and EHV installations are considered, the voltage factor is 1.1. With regard to the resistances of the conductors, IEC-60909-0 states that, when computing the maximum currents, they have to be computed at 20 °C (cold conductors). Let us suppose the fault levels (for computation of equivalent sequence and phase impedances see Appendix IV) of Table 4. In the same table the equivalent sequence impedances have been reported.
Table 4. Data assumed in the MCA for fault levels.
Table 4. Data assumed in the MCA for fault levels.
SubstationFault levels and equivalent sequence impedances
SEND. | I ¯ 3 P | = 42.62    k A ;       | I ¯ 1 P | = 43.42    k A     Z ¯ 1 = Z ¯ 2 = j 5.9604    Ω ;     Z ¯ 0 = j 5.6310    Ω
REC. | I ¯ 3 P | = 32.27    k A ;       | I ¯ 1 P | = 28.27    k A     Z ¯ 1 = Z ¯ 2 = j 7.8721    Ω ;     Z ¯ 0 = j 11.2137    Ω
Figure 16 shows the phase voltage magnitudes along the OHL. At midline the phase 1-R voltage magnitude is equal to the tower voltages (which is the same voltage of the earth wire). Also the phase 6-R of healthy circuit has a voltage drop whereas all the other phase voltage magnitudes have an overvoltage. These ones are the so called “temporary overvoltages due to earth faults” or power frequency overvoltages of relatively long duration.
Figure 17 shows the current sharing between the six phases: 6.85 kA coming from substation SEND. and 6.45 kA from substation REC. due to the different fault levels. The fault current |I1P-1R| is equal to 13.304 kA whereas healthy phases (or un-faulted phases) have almost zero current. Beyond the faulted phase 1-R, the only not-negligible short circuit current flows in the phase which is in parallel with 1-R i.e., 6-R (since “low reactance” phasing has been applied). The short circuit currents in all the other phases are almost zero since, according to IEC 60909-0, the shunt load admittance can be neglected so that the faulted network is at no-load.
Figure 16. Line-to-earth voltage magnitudes along the line.
Figure 16. Line-to-earth voltage magnitudes along the line.
Energies 05 02771 g016
Figure 17. Phase current magnitudes along the line.
Figure 17. Phase current magnitudes along the line.
Energies 05 02771 g017
The earth wire current magnitude |IEW| (see Figure 18) decreases from the fault location up to a given number of spans injecting current into the earth. This number of spans, after which there is not current injection into the earth anymore, can be theoretically evaluated [7,13]. By observing Figure 18 it can be ascertained that after 8 km from the fault location the earth wire stops injecting current into the earth. Figure 19 shows the ground return current magnitudes and the angles for a midline short-circuit along the line.
It is worth noting that the exact knowledge of the ground return current behaviour results fundamental when evaluations of EMI must be assessed in a given point along the line so that the only knowledge of the sending and receiving end values |IGR| are not enough.
Figure 18. Earth wire current magnitude |IEW| along the line in fault condition.
Figure 18. Earth wire current magnitude |IEW| along the line in fault condition.
Energies 05 02771 g018
Figure 19. Ground return current magnitude |IGR| and angle along the line in fault condition.
Figure 19. Ground return current magnitude |IGR| and angle along the line in fault condition.
Energies 05 02771 g019

Comparison with EMTP-RV

Some of the above presented MCA results are compared with the EMTP-RV ones for the OHL of previous section whose EMTP-RV electrical layout is shown in Figure 20.
Figure 20. EMTP-RV double-circuit OHL electrical layout (pages 1/6 and 6/6).
Figure 20. EMTP-RV double-circuit OHL electrical layout (pages 1/6 and 6/6).
Energies 05 02771 g020
Figure 21 shows the phase voltage magnitude behaviors along the double circuit OHL by means of MCA and EMTP-RV. The maximum percentage voltages for each phase are reported in Table 5 where the maximum value is equal to 0.95%. The cpu-time in the simulation with MCA is 3.41 s and with EMTP-RV is 3.5 s (INTEL CORE 2 QUAD [email protected] GHz, 3.25 GB RAM). |IGR| comparison is shown in Figure 22. The maximum difference between the MCA and EMTP-RV ground current magnitude is 77 A (−1.9%). In conclusion, the agreement between MCA and EMTP-RV is extremely good.
Figure 21. Line-to-earth voltage magnitude comparison between the MCA and EMTP-RV under phase-to-earth short circuit in a double-circuit OHL.
Figure 21. Line-to-earth voltage magnitude comparison between the MCA and EMTP-RV under phase-to-earth short circuit in a double-circuit OHL.
Energies 05 02771 g021
Table 5. Double-circuit line maximum voltage magnitude comparison.
Table 5. Double-circuit line maximum voltage magnitude comparison.
EliminateMax. faulted circuit voltage magnitudes [kV]Max. healthy circuit voltage magnitudes [kV]
Phase 1-R (faulted phase)Phase 2-SPhase 3-TPhase 6-RPhase 5-SPhase 4-T
MCA6,68307,7304,6140,4294,8303,8
EMTP-RV6,74308,8301,7141,4295,3301,9
Percentage difference−0.9%−0.36%0.95%−0.71%0.16%0.62%
Figure 22. Ground return current comparison between the MCA and EMTP-RV under line-to-earth short circuit in a double-circuit OHL.
Figure 22. Ground return current comparison between the MCA and EMTP-RV under line-to-earth short circuit in a double-circuit OHL.
Energies 05 02771 g022

5. Conclusions

MCA has already had a great number of successful applications to multi-phase asymmetrical systems [1,2,3]. In this paper, it has been applied to overhead lines with bundled phase conductors and one or more earth wires. The steady state regime of a typical un-transposed HV OHL has highlighted that there is a not negligible ground return current which may lead to electromagnetic interferences and ac corrosion issues. In faulty analysis, the behaviour of ground return current can be directly achieved so to recognize the number of spans after which the injection of currents in the ground terminates and the ground return current becomes constant. MCA is also extremely useful to evaluate the rise of earth potential at substation sites. Further researches are needed in order to validate the screening factor approach by means of MCA tool. Comparisons with EMTP-RV commercial software has demonstrated a very good agreement.

Appendix I. Equivalent Supply Network Matrix

In the study of multiconductor steady state regimes the following positive sequence voltage vector must be directly applied:
u ¯ p h = [ 1 α 2 α ] U N 3           α = e j 2 3 π
(e.g., for UN = 400 kV → U N / 3 230  kV ) and for the study of the multiconductor no-load energization, also a phase model of the supply network must be introduced by means of the matrix Zpha.
By assuming a perfectly symmetric structure of the three phase supply network it can be fully characterized by the corresponding zero, positive and negative longitudinal impedances derived from the sub-transient short-circuit currents (single-phase and three-phase ones). The corresponding phase matrix Zpha can be obtained by means of the Fortescue’s transformations:
u ¯ S = F ¯ u ¯ p h a    i ¯ S = F ¯ i ¯ p h a
u ¯ p h a = F ¯ 1 u ¯ S    i ¯ p h a = F ¯ 1 i ¯ S
If a passive symmetric tripole is considered the following diagonal matrix relations can be formed
Energies 05 02771 i004
By using Equations (A.2,A.3) and by arranging Equation (A.4), it is possible to compute the phase impedance matrix Zpha from the sequence impedance matrix ZS, as in Equation (A.8):
F ¯ u ¯ p h a = Z ¯ S i ¯ S
u ¯ p h a = F ¯ 1 Z ¯ S i ¯ S
u ¯ p h a = F ¯ 1 Z ¯ S F ¯ i ¯ p h a = Z ¯ p h a i ¯ p h a
Z ¯ p h a = F ¯ 1 Z ¯ S F ¯
By following the same criteria used for the longitudinal block L of Section 2 it is possible (see Figure A1) to achieve for the block A the matrix Ypha which characterizes the power supply in the sub-transient regimes.
Figure A1. (a) Supply network modelling (sequence); (b) Supply network modelling (Multiconductor).
Figure A1. (a) Supply network modelling (sequence); (b) Supply network modelling (Multiconductor).
Energies 05 02771 g023
If faulty conditions are considered, this model must be applied at both ends in order to take in account the two end substations (with their sub-transient fault levels).
Figure A2. Model of complex load at receiving end.
Figure A2. Model of complex load at receiving end.
Energies 05 02771 g024

Appendix II. Equivalent Load Matrix

In order to model a complex power loading at receiving-end, a three-phase balanced static load can be considered as in Figure A2. The value of yL is calculated so that the three admittances absorb, once applied the rated positive sequence voltage u, the rated receiving-end three-phase complex power SL:
y ¯ L = S ¯ L * 3 u 2

Appendix III. Earthing Technique

The earthing of own two earth wires (see Figure A3) in correspondence of the tower must be considered with its admittance matrix YTW (2 × 2). This matrix presents different composition depending on the earth wire earthing technique.
Figure A3. Two earth wires earthing at tower.
Figure A3. Two earth wires earthing at tower.
Energies 05 02771 g025
The formation of YTW is rather simple and involves the inversion of nodal impedance matrix i.e.,:
Energies 05 02771 i005
The “contact resistance” rtw accounts for the resistance due to the tower presence and can be set to a very small value i.e., rtw = 10−3 Ω. Differently Rt can range between about 5 and 20 Ω.
Figure A4. EW4 insulation from tower.
Figure A4. EW4 insulation from tower.
Energies 05 02771 g026
In some transmission grids, the EWs are insulated from the tower. This can be made in order to supply discharge lamps [14] which warn the presence of tower itself or to zero the EW power losses (usually for very long OHLs as in Brazil). For example, if EW4 is insulated from the tower (see Figure A4), the suitable admittance must be considered.

Appendix IV. Short-Circuit Analysis

The short-circuit, for OHL in Section 4, between phase 1-R and earth wire can be easily considered by constructing the following matrix, where numbers 1, 2, 3, … 6 are the phase conductors and number 7 is the earth wire (Figure A5):
Figure A5. Admittance matrix representing a line-to-earth short circuit.
Figure A5. Admittance matrix representing a line-to-earth short circuit.
Energies 05 02771 g027
where g must account for the short-circuit conductance (e.g., 1 kS). The matrix Yfault must be superimposed in YTOT in the location corresponding to the faulty section. By means of this matrix technique other types of fault can be easily assessed e.g. phase interruptions, multiple dead and not short-circuits. By considering the fault levels, the sequence impedances for each substation can be computed by the well-known formulae:
I ¯ 3 P = 1.1 E ¯ P V Z ¯ 1
I ¯ 1 P = 1.1 3 E ¯ P V Z ¯ 1 + Z ¯ 2 + Z ¯ 0
By solving the Equations (A11,A12) it yields:
Z ¯ 1 = Z ¯ 2 = 1.1 E ¯ P V I ¯ 3 P
Z ¯ 0 = 1.1 E ¯ P V ( 3 I ¯ 1 P 2 I ¯ 3 P )
.
Once the diagonal sequence matrix ZS has been computed, it is easy to compute the phase matrix Zpha by means of Fortescue matrix F (as clearly shown in Appendix I).
Once the YTOT has been computed the faulty regime can be evaluated with a suitable partition of it as in Figure A6.
Figure A6. Partitioned matrix of i = YTOT·u in short-circuit occurrence.
Figure A6. Partitioned matrix of i = YTOT·u in short-circuit occurrence.
Energies 05 02771 g028
The study of short-circuit leads to the following equations:
i ¯ p h a S = Y ¯ 1 u ¯ p h a S + Y ¯ 2 u ¯ x + Y ¯ 3 u ¯ p h a R
0 = Y ¯ 4 u ¯ p h a S + Y ¯ 5 u ¯ x + Y ¯ 6 u ¯ p h a R
i ¯ p h a R = Y ¯ 7 u ¯ p h a S + Y ¯ 8 u ¯ x + Y ¯ 9 u ¯ p h a R
where the vector ix is null (see Figure A6). From Equation (A16), it is rather immediate to obtain the unknown vector ux:
u ¯ x = Y ¯ 5 1 ( Y ¯ 4 u ¯ p h a S + Y ¯ 6 u ¯ p h a R )
Once again, if Y5 is of great dimensions, its inversion can benefit from the partition techniques and parallel computing.

Acknowledgments

This research has been financed by National Grid (contract reference TAO 22000, Title: Fault Current Distribution in Overhead Line and new Types of EHV Cables).

References

  1. Benato, R.; Paolucci, A. EHV AC Underground Electrical Power. Performance and Planning; Springer: Berlin, Germany, 2010. [Google Scholar]
  2. Benato, R. Multiconductor analysis of underground power transmission systems: EHV AC Cables. Electr. Power Syst. Res. 2009, 79, 27–38. [Google Scholar] [CrossRef]
  3. Benato, R.; Brunello, P.; Carlini, E.M.; di Mario, C.; Fellin, L.; Knollseisen, G.; Laußegger, M.; Muhr, M.; Paolucci, A.; Stroppa, W.; Wörle, H.; Woschitz, R. Italy-Austria GIL in the new planned railway galleries Fortezza-Innsbruck under Brenner Pass. Elektrotech. Informationstechnik 2006, 123, 551–558. [Google Scholar] [CrossRef]
  4. International Telecommunication Union. Directives Concerning the Protection of Telecommunication Lines against Harmful Effects from Electric Power and Electrified Railway Lines; ITU: Geneva, Switzerland, 1989. [Google Scholar]
  5. Carson, J.R. Wave propagation in overhead wires with ground return. Bell Syst. Tech. J. 1926, 5, 539–554. [Google Scholar] [CrossRef]
  6. Dommel, H.W. Electromagnetic Transients Program Reference Manual (EMTP Theory Book); Bonneville Power Administration: Portland, OR, USA, 1986. [Google Scholar]
  7. Nasser, D.; Tleis, A. Power Systems Modelling and Fault Analysis: Theory and Practice; Newnes: London, UK, 2008. [Google Scholar]
  8. Braae, R. Matrix Algebra for Electrical Engineers; Pitman & Sons: London, UK, 1963. [Google Scholar]
  9. Booth, E.S.; Clark, D.; Egginton, J.L.; Forrest, J.S. The 400 kV grid system for England and Wales. Proc. IEE Part A Power Eng. 1962, 109, 493–511. [Google Scholar] [CrossRef]
  10. AC Corrosion on Metallic Pipelines due to Interference from AC Power Lines—Phenomenon, Modelling and Countermeasures; Cigré Technical Brochure # 290; Cigré: Paris, France, 2006.
  11. British Electricity International. Modern Power Station Practice, EHV Transmission; Permagon: Oxford, UK, 1991. [Google Scholar]
  12. International Electrotechnical Commission. Short-Circuit in Three-Phase A.C. Systems. Part 1: Calculations of Current, 1st ed.; IEC 60909-0; International Electrotechnical Commission: Geneva, Switzerland, 2001. [Google Scholar]
  13. Paolucci, A. Lezioni di Trasmissione Dell’energia Elettrica; CLEUP: Padova, Italy, 1998. [Google Scholar]
  14. Benato, R.; Paolucci, A.; Turri, R. Insulated Ground wire capacitive currents for tower discharge warning lamp supplying. Electr. Power Syst. Res. 2004, 71, 211–221. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Benato, R.; Dambone Sessa, S.; Guglielmi, F. Determination of Steady-State and Faulty Regimes of Overhead Lines by Means of Multiconductor Cell Analysis (MCA). Energies 2012, 5, 2771-2793. https://doi.org/10.3390/en5082771

AMA Style

Benato R, Dambone Sessa S, Guglielmi F. Determination of Steady-State and Faulty Regimes of Overhead Lines by Means of Multiconductor Cell Analysis (MCA). Energies. 2012; 5(8):2771-2793. https://doi.org/10.3390/en5082771

Chicago/Turabian Style

Benato, Roberto, Sebastian Dambone Sessa, and Fabio Guglielmi. 2012. "Determination of Steady-State and Faulty Regimes of Overhead Lines by Means of Multiconductor Cell Analysis (MCA)" Energies 5, no. 8: 2771-2793. https://doi.org/10.3390/en5082771

APA Style

Benato, R., Dambone Sessa, S., & Guglielmi, F. (2012). Determination of Steady-State and Faulty Regimes of Overhead Lines by Means of Multiconductor Cell Analysis (MCA). Energies, 5(8), 2771-2793. https://doi.org/10.3390/en5082771

Article Metrics

Back to TopTop