Next Article in Journal
Computational Fluid Dynamic Analysis of Co-Firing of Palm Kernel Shell and Coal
Next Article in Special Issue
Data Distribution Service-Based Interoperability Framework for Smart Grid Testbed Infrastructure
Previous Article in Journal
Closed Solar House with Radiation Filtering Roof for Transplant Production in Arid Regions: Energy Consumption
Previous Article in Special Issue
Shielding Effectiveness Simulation of Small Perforated Shielding Enclosures Using FEM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate

Department of Astronautics, Electric and Energy Engineering, Sapienza University of Rome, Rome, 00184, Italy
*
Author to whom correspondence should be addressed.
Energies 2016, 9(3), 139; https://doi.org/10.3390/en9030139
Submission received: 12 January 2016 / Revised: 15 February 2016 / Accepted: 19 February 2016 / Published: 26 February 2016

Abstract

:
The paper deals with the impact of subsequent strokes on the backflashover rate (BFR) of HV overhead transmission lines (OHLs), assessed by means of an ATP-EMTP Monte Carlo procedure. The application to a typical 150 kV Italian OHL is discussed, simulating several tower grounding system arrangements. Subsequent strokes parameters are added to the statistical simulation variables: peak current, front time, time-to-half value, lightning polarity, line insulation withstand, lightning location and phase angle of the power frequency voltage. The input data are fed to an ATP-EMTP complete circuit model of the OHL, including line insulation, lightning representation and tower grounding system, the latter simulated by a pi-circuit model able to simulate the effects due to propagation and soil ionization, at modest computational costs. Numerical results evidence a non-negligible BFR increase (in relative terms) due to subsequent strokes: for spatially concentrated grounding systems the BFR increase approximatively vary in inverse proportion with the low frequency grounding resistance, whereas for spatially extended grounding systems the BFR increase depends on the grounding system behavior at high frequencies.

1. Introduction

Most faults affecting overhead transmission lines (OHLs) at HV and EHV level are caused by lightning strokes to the line itself (i.e., direct lightning). Due to the presence of shield wires, strokes to tower peaks and the shield wires largely prevail on those to phase conductors; as a consequence, faults associated with tower/shield wire strokes (backflashovers) are also more numerous than those caused by strokes to phase conductors (shielding failures). Due to the random nature of lightning and line insulation behavior, quantitative analyses aimed at assessing backflashovers require, in principle, a multivariable, statistical approach, ideally suited to Monte Carlo methods.
The main drawback of every Monte Carlo-based approach is that a large number of simulations is required in order to obtain convergence. Since the tower grounding system impulse response plays a key role in the lightning performance of OHLs, especially if spatially extended grounding system arrangements are employed (e.g., in medium-to-high resistivity soils), backflashover analyses become too computationally expensive when accurate grounding system models [1,2,3,4] are implemented. Therefore, on one hand simplified methods not based on a Monte Carlo approach were developed, notably the CIGRE [5] and IEEE [6,7] procedures. On the other hand, several procedures based on a Monte Carlo approach have been proposed [8,9,10,11,12,13], but adopting very simple grounding system models. In [8], a user-defined tower footing resistance distribution is considered. In [10], impulse current resistances, derived from laboratory test data, simulate concentrated grounds. In [11], a constant resistance is used, whereas in [12] grounding system model consists of a linear lumped resistance R and an inductance L. Lastly, in [13] the tower footing “impulse resistance” is approximated by the non-linear resistance model proposed in [1] and a log-normal distribution of the low-frequency resistance is given. In previous papers, the authors proposed a simplified (yet accurate) grounding system model [14,15,16,17], able to reproduce the lightning responses, including soil ionization phenomena, at modest computational costs; the model has been subsequently employed in an ATP-EMTP-based Monte Carlo procedure in order to assess the backflashover rate (BFR) of OHLs. All these studies, however, only focus on the first stroke, since it almost invariably has the highest peak current in a given lightning flash. In recent years, the influence of the smaller (but steeper) subsequent strokes on BFR was studied in several non-Monte Carlo-based papers [18,19,20], showing a non-negligible BFR increase. In [18,19,20], however, grounding systems are simulated by linear lumped resistances, thus not taking into account both soil ionization (relevant especially in low-resistivity soils) and the high-frequency response to the steeper wavefronts associated to subsequent strokes (relevant especially in high-resistivity soils, where spatially extended grounding systems are employed).
The paper presents a comprehensive ATP-EMTP-based backflashover performance evaluation tool, building up subsequent strokes into the authors’ Monte Carlo procedure [21,22,23]. The latter previously only took into account first strokes, whereas the presented procedure assumes a number of subsequent strokes statistically variable in accordance with data reported in [24]. Results pertaining to typical 150 kV-50 Hz Italian OHLs equipped with a wide variety of grounding systems, both spatially concentrated and extended (i.e., counterpoises and crowfoot arrangements), are presented, showing the effect of subsequent strokes on BFR estimation. The paper is structured as follows: Section 2 recalls the BFR calculation procedure and Section 3 details ATP-EMTP simulation models, whereas Section 4 reports the statistical variables and parameters used in the procedure; lastly, simulation results are shown and analyzed in Section 5.

2. Basic Considerations

The statistical ATP-EMTP-based Monte Carlo procedure generates Ntot lightnings, which are supposed to fall within a 1 km-wide zone having the OHL placed in the center. The NL lightnings impinging on the line are found via the Eriksson’s electrogeometric model [25], assuming that all strokes terminate on towers. Such strokes are simulated with ATP-EMTP, in order to check the occurrence of backflashover. The BFR of the line (faults/100 km·year) is assessed as:
B F R = 0.6 N B F O N t o t N g 100
where NBFO is the total number of backflashovers resulting from simulations, Ng is the ground flash density (flashes/km2·year) and 0.6 is a numerical coefficient which multiplies the BFR calculated considering only strokes to tower, in order to account for strokes within the span, as in [26]. Figure 1 shows a flow-chart of the procedure [21,22,23] (in the flow-chart N is the number of the current simulation, TgN is the NBFO/N ratio calculated by the procedure and εN is the convergence error).

3. System Modeling

The simulated system is a 10 km long stretch of a typical Italian 150 kV-50 Hz OHL. Seven different tower grounding arrangements were considered: three concentrated ones, shown in Figure 2, and four spatially extended ones, shown in Figure 3; main features of the simulated grounding systems are reported in Table 1.

3.1. OHL Model

The ATP-EMTP OHL model consists of 25 line spans, each 400 m long and simulated by the “JMarti” frequency-dependent model. The single circuit OHL is equipped with single 31.5 mm ACSR conductors (phases) and a single 11.5 mm steel conductor (ground wire): Figure 4 shows the outline of the tower. At both ends of the simulated line stretch, the OHL is connected to the line surge impedance matrix: phase conductors are then terminated on a three-phase 150 kV-50 Hz voltage system, whereas the shield wire is solidly grounded [15,16,27]. Table 2 reports the phase and shield wire conductors coordinates, whereas the phase conductor and the shield wire sags are 11.4 m and 9.7 m, respectively.

3.2. Line Insulation Model

Line insulation has been represented by means of the well-established CIGRE Leader Progression Model (LPM), programmed into ATP-EMTP by using “MODELS” language [28]:
d l d t = k u ( t ) [ u ( t ) d G l ( t ) E 0 ]
being dG (m) the gap length (1.46 m in the simulations reported in the paper), l (t) (m) the leader length, u(t) (kV) the voltage across tower insulation, E0 (kV/m) the critical electric field, whereas the speed parameter k (m2·kV−2·s−1) is 1.2·10−6 and 1.3·10−6 for positive and negative flashes, respectively.

3.3. Lightning Model

Both first and subsequent strokes have been simulated by the well known “Heidler” current surge function [29] implemented in ATP-EMTP:
i ( t ) = I P η k s n 1 + k s n e t τ 2
being IP the peak value of the lightning stroke, ks the t1 ratio, τ1 and τ2 time constants establishing the front time and the time-to-half value, respectively, η the adjustment constant and n a factor influencing the rate of rise and the time instant of the maximum rate of rise.

3.4. Grounding System Model

Tower grounding systems were simulated by means of the simplified pi-circuit model described in [14,15,17] and shown in Figure 5. The pi-circuit model is obtained by synthesizing the full circuit model described in [1] (also accounting for soil ionization), by means of a procedure based on the micro-Genetic Algorithm. Linear circuit elements in Figure 5 (R1, R2, C1, C2, R, L) reproduce the frequency response of the non-ionized grounding system in the frequency range 1 Hz–1 MHz.
The ideal voltage-controlled current sources G1 and G2 simulate the non-linear behavior of the grounding system due to soil ionization and are defined as (i = 1, 2):
G i ( t ) = V R i ( t ) F i ( t ) V R i ( t ) R i
where:
F i ( t ) = R i α i log ( 10 4 + β i V R i ( t ) R i ) ,    F i [ 10 4 ; R i ]
and V R i ( t ) is the voltage across the resistor Ri. The pi-circuit model has been also validated by comparison with two different models, one based on an electromagnetic full-wave approach, the other adopting hybrid circuit-field approach [30], yielding very good results even when subsequent strokes are injected. Numerical values of pi-circuit parameters of the simulated grounding systems are reported in Table 3 and Table 4, whereas the impedances in frequency domain (in the frequency range 1 Hz–1 MHz) of the concentrated and spatially extended grounding systems are plotted in Figure 6. Figure 6a shows that concentrated grounding systems GS 1, GS 2 and GS 3 have a predominant resistive behavior in the whole frequency range of interest, whereas spatially extended ones have a predominant inductive (GS 4, GS 6 and GS 7) and capacitive (GS 5) behavior at high frequencies.

4. Statistical Inputs of the Procedure

4.1. Lightning Polarity

Since lightning polarity influences the generation of the other statistical inputs, it is sampled first, as a random variable uniformly distributed between 0 and 1, assuming that 10% of flashes to ground are positive.

4.2. Lightning Parameters

Log-normal distributions were assumed for lightning stroke parameters (peak current IP, front time tF and time-to-half-value tH), with the values reported in Table 5 [24] for first strokes. Parameters for subsequent negative return-strokes are summarized in Table 6 [24,31], whereas the statistical distribution of multistroke negative lightning flashes is reported in Table 7 [24,31]. According to [24], no correlation between the first and the subsequent stroke peak amplitudes has been assumed.

4.3. Line Insulation Parameters

Critical field E0 in the LPM is also treated as a log-normal distributed variable. Median values E0m and standard deviations [5] are reported in Table 8.

4.4. Lightning Location and Attachment

To check the actual occurrence of a lightning stroke to the OHL, the position of the lightning in a 1 km wide strip, centered on the OHL (i.e., its initial distance from the line, assuming a vertical channel), is generated as a random, uniformly distributed variable in the −500 m, +500 m interval. To check lightning attachment to the line, this distance is compared to the attractive radius Ra of the OHL, calculated by using the Eriksson’s electrogeometric model [25] as:
R a = 0.67 H 0.6 I P 0.74
where H is the tower height (m).

4.5. Phase Angle of the Power Frequency Voltage

A variable uniformly distributed in the range 0°–360° simulates the phase angle of the impressed power frequency voltages (three-phase positive-sequence system).

5. Results

For the case under study, the total number of generated lightnings is Ntot = 720,345, corresponding to 2,153,833 lightning strokes, of which NL = 200,000 lightning strokes impinge on the line. Figure 7 plots the peak current distribution of the NL lightning strokes. In order to show the self-consistency of the Monte Carlo procedure, i.e., the actual convergence within the given number Ntot of lightning strokes, Figure 8 shows the 100·(NBFO/N) ratio vs. N yielded by the procedure for GS 6 grounding arrangement, for both cases of subsequent strokes disregarded or accounted for (in [32] 100·(NBFO/N) ratio vs. N trends for GS 1, GS 2 and GS 3 grounding systems are reported). Table 9 summarizes the smallest lightning currents (also the corresponding front times are indicated) causing BFR, both when single flashes and multiple flashes are simulated. Table 10 reports the calculated 100·(NBFO/N) ratios obtained for the seven different grounding systems with or without subsequent strokes, as well as the corresponding increase of the outage rate due to subsequent strokes.
The “overall” peak current distribution for strokes to the OHL in Figure 7 includes both polarities of first strokes, as well as subsequent negative strokes. It still resembles a log-normal probability distribution function, although three different log-normal distributions are superimposed and, moreover, original values are filtered by the electrogeometric model assessment of impingement on the line.
Results in Table 9 show some points:
  • the minimum current causing backflashover strictly depends on R@50Hz for concentrated grounding systems, decreasing as R@50Hz increases;
  • for spatially extended grounding systems, the minimum current causing backflashover also depends on the their frequency behavior: for instance, GS 5 has a larger R@50Hz than GS 7, but even a larger minimum current causing backflashover, due to its capacitive behavior at high frequency with respect to the inductive behavior of GS 7.
BFR results in Table 10 show some important points:
  • for concentrated, predominantly resistive, grounding systems GS 1, GS 2 and GS 3, the incorporation of subsequent strokes into the Monte Carlo procedure causes a non-negligible increase of BFR values. A qualitative trend may be evidenced: the impact is noticeably higher (in relative terms) for lower values of low-frequency grounding resistance, as explained in [18];
  • for spatially extended grounding systems, the impact of subsequent strokes on BFR depends on the grounding system behavior at high frequencies and is not related to its R@50Hz value. Moreover, the qualitative trend found for concentrated grounding systems does not appear: the impact of subsequent strokes is generally smaller (in relative terms) than for concentrated ones, but not negligible. The largest BFR increase is found in conjunction with the predominantly inductive "counterpoise" grounding system GS 7, caused by the steeper subsequent stroke lightning waves with respect to first strokes, reflecting on a larger L d ( i ) d t contribution to the ground potential rise.
Bar charts in Figure 9 and curves in Figure 10 show the distinctive effect of subsequent strokes on BFR for concentrated and spatially extended grounding arrangements, respectively, i.e., the appearance of a narrow but not negligible band of “critical” peak current values to the left of the plot.
A comparison with prior work, notably the wide-range study in [20], shows that the qualitative trends (BFR vs. low-frequency grounding resistance) are fully consistent, whereas the BFR increases calculated by the authors’ procedure are smaller than those in [20]. This quantitative difference is probably due to the different modeling approaches, since [20] simulates grounding systems by a linear lumped resistance, thus disregarding soil ionization, which might be quite significant for spatially concentrated grounding systems in low-resistivity soils, as well as the inductive/capacitive behavior at high frequencies, which is relevant for spatially extended grounding systems.

6. Conclusions

The paper presented a comprehensive ATP-EMTP-based Monte Carlo procedure for backflashover assessment, enhancing previous studies with the consideration of the effect of subsequent strokes. Computationally cheap reproduction of the high-frequency behavior of tower grounding systems subjected to the steep lightning current wavefronts of subsequent strokes is allowed by the authors’ simple, non-linear, wideband grounding system circuit model. A preliminary study was carried out for a typical Italian 150 kV OHL, considering seven different tower grounding arrangements, three spatially concentrated and four extended. For all cases, a non-negligible impact of subsequent strokes on the BFR estimation was evidenced, with calculated increases ranging between 6% and 25%. For concentrated, predominantly resistive, grounding systems the incorporation of subsequent strokes into the Monte Carlo procedure causes a non-negligible increase of BFR values. A qualitative trend may be evidenced: the effect of subsequent strokes is noticeably higher (in relative terms) as the "grounding resistance" decreases, as found and explained in previous papers available in literature. From a quantitative point of view, the authors’ procedure assesses a smaller increase than those reported in literature, probably due to different modeling approaches: previous papers model grounding systems as linear lumped resistances, thus disregarding soil ionization, which is quite significant for spatially concentrated grounding systems in low-resistivity soils. The response of the different studied spatially extended grounding systems depends on their high frequency behavior, which can be either inductive or capacitive. The largest BFR increase is found in conjunction with the predominantly inductive “counterpoise” grounding systems, suggesting that in these cases subsequent strokes should not be neglected. On the other hand, the mainly capacitive "crowfoot" grounding system is less affected, in relative terms, by the inclusion of subsequent strokes: the behavior of such capacitive grounding systems is quite similar to that of resistive grounding systems.

Author Contributions

The authors contributed equally to the work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Geri, A. Behaviour of grounding systems excited by high impulse currents: The model and its validation. IEEE Trans. Power Deliv. 1999, 14, 1008–1017. [Google Scholar] [CrossRef]
  2. Mentre, F.E.; Grcev, L. EMTP-based model for grounding system analysis. IEEE Trans. Power Deliv. 2005, 20, 2200–2210. [Google Scholar] [CrossRef]
  3. Cidras, J.; Otero, A.F.; Garrido, C. Nodal frequency analysis of grounding systems considering the soil ionization effect. IEEE Trans. Power Deliv. 2000, 15, 103–107. [Google Scholar] [CrossRef]
  4. Visacro, S.; Soares, A. HEM: A model for simulation of lightning-related engineering problems. IEEE Trans. Power Deliv. 2005, 20, 1206–1208. [Google Scholar] [CrossRef]
  5. CIGRE Working Group 01 of SC 33. Guide to Procedures for Estimating the Lightning Performance of Transmission Lines; International Council on Large Electric Systems: Paris, France, 1991. [Google Scholar]
  6. IEEE Working Group. A simplified method for estimating lightning performance of transmission lines. IEEE Trans. Power App. Syst. 1985, 104, 919–932. [Google Scholar]
  7. IEEE Working Group. Estimating lightning performance of transmission lines II—Updates to analytical models. IEEE Trans. Power Deliv. 1993, 8, 1254–1267. [Google Scholar]
  8. Anderson, J.G. Monte Carlo computer calculation of transmission line lightning performance. IEEE Trans. Power App. Syst. 1961, 80, 414–419. [Google Scholar] [CrossRef]
  9. Sargent, M.A.; Darveniza, M. The calculation of double circuit outage rate of transmission lines. IEEE Trans. Power App. Syst. 1967, 86, 665–678. [Google Scholar] [CrossRef]
  10. Darveniza, M.; Sargent, M.A.; Limbourn, G.J.; Liew, A.C.; Caldwell, R.O.; Currie, J.R.; Holcombe, B.C.; Stillman, R.H.; Frowd, R. Modeling for lightning performance calculations. IEEE Trans. Power App. Syst. 1979, 98, 1900–1908. [Google Scholar] [CrossRef]
  11. Ekonomou, L.; Gonos, I.F.; Stathopulos, I.A. An improved backflashover model for estimating the lightning performance of transmission lines. In Proceedings of the 38th International Universities Power Engineering Conference (UPEC), Thessaloniki, Greece, 1–3 September 2003; pp. 109–112.
  12. Ekonomou, L.; Gonos, I.F.; Stathopulos, I.A.; Dialynas, E.N. A backflashover model for calculating the transmission lines lightning performance. In Proceedings of the 12th International Symposium on High-Voltage Engineering (ISH), Bangalore, India, 20–24 August 2001.
  13. Martinez, J.A.; Castro-Aranda, F. Lightning performance analysis of overhead transmission lines using the EMTP. IEEE Trans. Power Deliv. 1994, 9, 1838–1849. [Google Scholar] [CrossRef]
  14. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Equivalent Lumped Parameter Π-Network of Typical Grounding Systems for Linear and Non-Linear Transient Analysis. In Proceedings of the 2009 IEEE Bucharest Powertech, Bucharest, Romania, 28 June–2 July 2009; pp. 1–6.
  15. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Simplified HV tower grounding system model for backflashover simulation. In Proceedings of the 30th International Conference on Lightning Protection (ICLP), Cagliari, Italy, 13–17 September 2010; Volume 85, pp. 16–23.
  16. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Simplified HV tower grounding system model for backflashover simulation. Electr. Power Syst. Res. 2012, 85, 16–23. [Google Scholar] [CrossRef]
  17. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Generalized pi-circuit tower grounding model for direct lightning response simulation. Electr. Power Syst. Res. 2014, 116, 330–337. [Google Scholar] [CrossRef]
  18. Silveira, F.H.; Visacro, S.; De Conti, A.; de Mesquita, C.R. Backflashovers of Transmission Lines Due to Subsequent Lightning Strokes. IEEE Trans. Electromagn. Compat. 2012, 54, 316–322. [Google Scholar] [CrossRef]
  19. Silveira, F.H.; Visacro, S.; de Souza, R.E.; De Conti, A. Assessing the Impact of Subsequent Strokes on the Lightning Performance of Transmission Lines of 138 kV. In Proceedings of the 2012 International Conference on Lightning Protection (ICLP), Vienna, Austria, 2–7 September 2012.
  20. Silveira, F.H.; Visacro, S.; De Conti, A. Lightning Performance of 138-kV Transmission Lines: The Relevance of Subsequent Strokes. IEEE Trans. Electromagn. Compat. 2013, 55, 1195–1200. [Google Scholar] [CrossRef]
  21. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M.; Santarpia, A. An ATP-EMTP Monte Carlo procedure for backflashover rate evaluation. In Proceedings of the 2012 International Conference on Lightning Protection (ICLP), Vienna, Austria, 2–7 September 2012.
  22. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M.; Santarpia, A. An ATP-EMTP Monte Carlo procedure for backflashover rate evaluation: A comparison with the CIGRE method. Electr. Power Syst. Res. 2014, 113, 134–140. [Google Scholar] [CrossRef]
  23. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M.; Palone, F. Tower grounding improvement vs. line surge arresters: comparison of remedial measures for high-BFOR subtransmission lines. IEEE Trans. Ind. Appl. 2015, 51, 4952–4960. [Google Scholar] [CrossRef]
  24. Chowdhuri, P.; Anderson, J.G.; Chisholm, W.A.; Field, T.E.; Ishii, M.; Martinez, J.A.; Marz, M.B.; McDaniel, J.; McDermott, T.R.; Mousa, A.M.; et al. Parameters of lightning strokes: A review. IEEE Trans. Power Deliv. 2005, 20, 346–358. [Google Scholar]
  25. Eriksson, A.J. An Improved Electrogeometric Model for Transmission-Line Shielding Analysis. IEEE Trans. Power Deliv. 1987, 2, 871–886. [Google Scholar] [CrossRef]
  26. Hileman, A.R. Insulation Coordination for Power Systems, 1st ed.; Marcel Dekker, Inc.: New York, NY, USA, 1999. [Google Scholar]
  27. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Backflashover simulation of HV transmission lines with enhanced counterpoise groundings. Electr. Power Syst. Res. 2009, 79, 1076–1084. [Google Scholar]
  28. Alternative Transients Program (ATP) Rule Book. Canadian/American EMTP User Group, Oregon, USA, 1995.
  29. Heidler, F.; Cvetić, J.M.; Stanić, B.V. Calculation of lightning current parameters. IEEE Trans. Power Deliv. 1999, 14, 399–404. [Google Scholar]
  30. Araneo, R.; Maccioni, M.; Lauria, S.; Geri, A.; Gatta, F.M.; Celozzi, S. Hybrid and Pi-Circuit Approaches for Grounding System Lightning Response. In Proceedings of the 2015 IEEE Eindhoven Powertech, Eindhoven, Netherlands, 29 June–2 July 2015; pp. 1–6.
  31. Anderson, R.B.; Eriksson, A.J. A Summary of Lightning Parameters for Engineering Applications. In Proceedings of the International Council on Large Electric Systems, Paris, France, August 1980.
  32. Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate. In Proceedings of the 2015 IEEE International Conference on Environment and Electrical Engineering (EEEIC), Rome, Italy, 10–13 June 2015; pp. 1210–1215.
Figure 1. Flow-chart of the procedure.
Figure 1. Flow-chart of the procedure.
Energies 09 00139 g001
Figure 2. Sketches of the simulated concentrated grounding systems: (a) GS 1; (b) GS 2; (c) GS 3.
Figure 2. Sketches of the simulated concentrated grounding systems: (a) GS 1; (b) GS 2; (c) GS 3.
Energies 09 00139 g002
Figure 3. Sketches of the simulated spatially extended grounding systems: (a) GS 4; (b) GS 5; (c) GS 6; (d) GS 7.
Figure 3. Sketches of the simulated spatially extended grounding systems: (a) GS 4; (b) GS 5; (c) GS 6; (d) GS 7.
Energies 09 00139 g003
Figure 4. OHL tower head (dimensions in m; the tower height is 31.1 m).
Figure 4. OHL tower head (dimensions in m; the tower height is 31.1 m).
Energies 09 00139 g004
Figure 5. Pi-circuit model.
Figure 5. Pi-circuit model.
Energies 09 00139 g005
Figure 6. Grounding system impedance vs. frequency (no soil ionization): (a) Concentrated grounding systems; (b) Spatially extended grounding systems.
Figure 6. Grounding system impedance vs. frequency (no soil ionization): (a) Concentrated grounding systems; (b) Spatially extended grounding systems.
Energies 09 00139 g006
Figure 7. Peak current distribution of lightning strokes impinging on the line.
Figure 7. Peak current distribution of lightning strokes impinging on the line.
Energies 09 00139 g007
Figure 8. 100·(NBFO/N) ratio vs. N calculated by the Monte Carlo procedure for GS 6 grounding system.
Figure 8. 100·(NBFO/N) ratio vs. N calculated by the Monte Carlo procedure for GS 6 grounding system.
Energies 09 00139 g008
Figure 9. Probability distribution of peak currents resulting in a backflashover. (a) GS 1; (b) GS 2; (c) GS 3 grounding systems.
Figure 9. Probability distribution of peak currents resulting in a backflashover. (a) GS 1; (b) GS 2; (c) GS 3 grounding systems.
Energies 09 00139 g009
Figure 10. Probability distribution of peak currents resulting in a backflashover. (a) GS 4; (b) GS 5; (c) GS 6; (d) GS 7 grounding systems.
Figure 10. Probability distribution of peak currents resulting in a backflashover. (a) GS 4; (b) GS 5; (c) GS 6; (d) GS 7 grounding systems.
Energies 09 00139 g010
Table 1. Main features of the simulated grounding systems.
Table 1. Main features of the simulated grounding systems.
Grounding system configurationSoil resistivity ρg (Ω·m)Soil relative permittivity εrBurying depth (m)R@50Hz (Ω)
GS 1100500.86.3
GS 2100500.811.42
GS 3500430.814.95
GS 41000350.819.1
GS 52000190.830.4
GS 61000350.824.05
GS 72000190.824.5
Table 2. Phase and shield wire conductor coordinates.
Table 2. Phase and shield wire conductor coordinates.
CoordinateABCSW
x (m)−2.93−3.50
y (m)25.423.421.431.1
Table 3. Linear parameter values for pi-type simplified equivalent circuits.
Table 3. Linear parameter values for pi-type simplified equivalent circuits.
Grounding SystemR1 (Ω)R2 (Ω)R (Ω)L (μH)C1 (nF)C2 (nF)
GS 127.71.656.550.320.6625
GS 214.410.944.095.470.0115.7
GS 332.125.92.224.853.188.66
GS 434.1243.27012.53.329
GS 542.4969.3937.5410.26.7812.5
GS 628.75137.359.5412.51.085.67
GS 740.9660.990.0728.22.8910.5
Table 4. Non-linear coefficient values for pi-type simplified equivalent circuits.
Table 4. Non-linear coefficient values for pi-type simplified equivalent circuits.
Grounding Systemα1 (Ω)α2 (Ω)β1 (A1)β2 (A1)
GS 11.2040.045326.6163.1
GS 20.70.70.633850.5
GS 30.3220.630.633850.5
GS 40.4060.2851654.9931.0
GS 50.0053.1421090.71847.6
GS 60.4194.9153.685172.9
GS 702.832200.10.074
Table 5. Statistical parameters of first negative and positive return-strokes current.
Table 5. Statistical parameters of first negative and positive return-strokes current.
ParameterMedian ValueStandard Deviation
++
IP31.1 kA35 kA0.481.21
tF3.83 μs22 μs0.551.23
tH77.5 μs230 μs0.581.33
Table 6. Statistical parameters of subsequent negative return-strokes current.
Table 6. Statistical parameters of subsequent negative return-strokes current.
ParameterMedian ValueStandard deviation
IP12.3 kA0.5296
tF0.67 μs1.0131
tH30.2 μs0.93
Table 7. Statistical distribution of multistroke negative lightning flashes.
Table 7. Statistical distribution of multistroke negative lightning flashes.
Number of Strokes per FlashFrequency of Occurrence (%)Number of Strokes per FlashFrequency of Occurrence (%)
14564
21473
3983
4892
5810 or more4
Table 8. Parameters of line insulation.
Table 8. Parameters of line insulation.
PolarityMedian Value (kV/m)Standard Deviation (kV/m)
60518.15
+56016.80
Table 9. Minimum currents causing backflashover.
Table 9. Minimum currents causing backflashover.
Grounding System ArrangementR@50Hz (Ω)Minimum Current Causing Backflashover Single FlashMinimum Current Causing Backflashover Multiple Flash
IP (kA)tF (μs)IP (kA)tF (μs)
GS 16.348.830.5134.730.29
GS 211.4245.110.5934.730.29
GS 314.9543.710.5334.020.23
GS 419.0743.710.5332.280.24
GS 624.0541.910.6132.110.23
GS 724.535.100.4129.030.21
GS 530.437.230.6329.030.21
Table 10. Calculated 100·(NBFO/Ntot) ratios and corresponding increases due to the contribution of subsequent strokes.
Table 10. Calculated 100·(NBFO/Ntot) ratios and corresponding increases due to the contribution of subsequent strokes.
Grounding System Arrangement100·(NBFO/Ntot) Single Flash100·(NBFO/Ntot) Multiple FlashΔ (%)
GS 10.160.199+24.375
GS 20.6670.765+14.69
GS 30.7810.844+8.06
GS 41.4541.545+6.26
GS 61.9712.101+6.60
GS 72.6782.913+8.78
GS 53.3733.594+6.55

Share and Cite

MDPI and ACS Style

Gatta, F.M.; Geri, A.; Lauria, S.; Maccioni, M. Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate. Energies 2016, 9, 139. https://doi.org/10.3390/en9030139

AMA Style

Gatta FM, Geri A, Lauria S, Maccioni M. Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate. Energies. 2016; 9(3):139. https://doi.org/10.3390/en9030139

Chicago/Turabian Style

Gatta, Fabio Massimo, Alberto Geri, Stefano Lauria, and Marco Maccioni. 2016. "Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate" Energies 9, no. 3: 139. https://doi.org/10.3390/en9030139

APA Style

Gatta, F. M., Geri, A., Lauria, S., & Maccioni, M. (2016). Monte Carlo Evaluation of the Impact of Subsequent Strokes on Backflashover Rate. Energies, 9(3), 139. https://doi.org/10.3390/en9030139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop