IPST 2015 CAVTAT CROATIA

Investigation into Transient SFO, FFO, VFTO

Overvoltage Characteristics for Typical
Gas Insulated Substations
L. Czumbil, J. Kim, H. Nouri

Abstract--Overvoltage characteristics of typical single bus, conventional substations [3] [4].

double bus and one and a half bus GIS configurations under The paper aims to identify transient overvoltage
transient SFO, FFO and VFTO conditions are studied. The characteristics of typical Gas Insulated Substations under
transient conditions are simulated through load rejections,
various transient conditions, namely slow front overvoltages
lightning and the opening and closing of circuit breakers or
disconnectors. Surge impedance and travel time are used for (SFO), fast front overvoltages (FFO) and respectively very fast
defining the distributed parameter models of the GIS. The results transient overvoltages (VFTO).
suggest that for VFTOs the magnitude of the generated
overvoltage solely depends on the switching sequences while for II. GIS MODELING
SFOs and FFOs the overvoltage depends on the degree of the
In order to investigate SFO, FFO and VFTO overvoltage
inductive load and the type and length of transmission lines.
characteristics of typical GIS configurations a 380 kV single
Keywords: Gas insulated substation (GIS), transient power and single outlet GIS is considered (see Fig. 1). On the
overvoltages, SFO, FFO, VFTO, modeling. power side the GIS substation is connected to the grid through
a 100 km long single circuit 380 kV/50 Hz overhead line
I. INTRODUCTION (OHL) to the grid, while on the outlet side the GIS is
connected to a medium voltage distribution station with a total
T HE rapid advancement in industrial technologies and
population growth during the last decades has
considerably increased the energy demand. Due to
load of 100 MVA through a 120 MVA, 380/33 kV star-star
configuration transformer and respectively a 10 km long 33 kV
environmental constraints and reliability of the system Gas transmission line.
Insulated Substations (GIS) have become a major component
in today’s power networks [1].
The GIS, which is filled with pressurized SF6 gas for
electrical insulation and rapid arc extinction presents a series
of advantages against classical air insulated substations, like: Fig. 1. Investigated GIS substation.
small ground space requirements, reduced maintenance, high
reliability and protection from pollution. Therefore, these are Due to the traveling nature of transients, the investigated
particularly used in large cities, industrial townships, deserts GIS is modeled by an equivalent electrical circuit composed
and arctic areas [2]. by distributed parameter elements. Surge impedance and travel
Despite these merits, the GIS has its own unique problems. time have been used for defining the distributed parameter
These problems include an increase of overvoltages caused by model. The inner system, which consists of the high-voltage
transient waves reflected from different connections, low surge bus duct and the inner surface of the encapsulation, has been
impedance, and the decreased length of the conductors in these thoroughly represented by line sections, modeled as the
substations. Considering these issues, the propagation of these distributed parameter transmission lines [1]:
waves along the conductor will increase compared with
L l 1
ZS  ,   ,  , (1)
This work was supported by Hyundai Heavy Industries Co., Ltd. Dong-gu
C  L C
Ulsan, South Korea.   lnb a  2   0  R
L 0 , C
lnb a 
(2)
L. Czumbil and H. Nouri are with the Power Systems, Electronics and 2 
Control research Laboratory at University of the West of England (UWE), where: Z S is equipment surge impedance,  is travel time, 
Bristol, U.K.
J. Kim is with Electro Electric Systems Research Institute, Hyundai Heavy is waveform propagation speed, a is bus duct conductor outer
Industries Co., Ltd. Dong-gu Ulsan, South Korea and is presently a visiting radius, b is enclosure inner radius, l is equipment length, L and
research scholar at the Power Systems, Electronics and Control research C are the inductance and capacitance of the equipment,
Laboratory, University of the West of England (UWE), Bristol, U.K.
respectively.
Paper submitted to the International Conference on Power Systems
Transients (IPST2015) in Cavtat, Croatia June 15-18, 2015
IPST 2015 CAVTAT CROATIA
TABLE I loads is studied. Therefore, an 80 MVA load with various
EQUIVALENT ELECTRICAL PARAMETERS OF GIS COMPONENTS [4] inductive and capacitive power factors between 0.97 and 0.4 is
Component Notes
rejected while a mostly resistive 20 MVA load (with a power
GIS Bus Duct ZS = 95 Ω
factor of 0.95) is maintained continuously connected to the
in closed position ZS = 42
medium voltage distribution station. The overvoltages at both
Disconnectors (DS) Ω
in open position C = 4 pF
medium and high voltage sides of the 380/33 kV, 120 MVA
in closed position ZS = 66 transformer placed at the outlet side of the GIS are recorded.
Circuit Breakers (CB) Ω Two different scenarios are considered when the medium
in open position C = 4 pF voltage distribution station is connected to the GIS through an
Potential Transformers (PT) ZS = 25 Ω, C = 10 pF overhead 33kV power line and three single core 33kV
Current Transformers (CT) ZS = 42 Ω underground power cables, respectively:
Spacers, Elbows C = 10 pF

For the transient overvoltage studies carried out on the

investigated 380 kV gas insulated substation three different
GIS configurations were analyzed: a single bus (Fig. 2a),
double bus (Fig. 2b) and a one and a half bus (Fig. 2c)
configuration, respectively.

Fig. 3. SFO variation with the power factor of the rejected load.

Fig. 3 presents the variation of the recorded overvoltage at

the medium voltage side of the transformer according to the
power factor of the rejected inductive and capacitive loads. It
can be observed that the rejection of an inductive load
produced a 1.25 p.u. overvoltage with a slight decrease with
the increase of the inductive part of the rejected load, while the
overvoltage produced by the rejection of the capacitive load
presents a ‘V’ curve with a minimum overvoltage at a 0.65
power factor. Higher overvoltage levels could be recorded for
the capacitive load rejection when the distribution station is
connected to the GIS through the 33kV OHL. Due to the
galvanic separation between the two winding of the 380/33 kV
transformer the recorded overvoltage at the high voltage side
(GIS outlet) is less than 1.05 p.u.
The influence of the transmission line length connecting the
medium voltage distribution station to the GIS, on the
overvoltage produced by load rejection has also been
investigated. The obtained overvoltage values at the medium
voltage side for different underground cable and OHL lengths
are presented in Fig. 4 and relate to the 100 MVA load of the
medium voltage distribution station being rejected:

Fig. 2. Investigated GIS configurations: a) single bus GIS; b) double bus GIS;
c) one and half bus GIS;

III. SLOW FRONT OVERVOLTAGE STUDY

In high voltage power systems slow front transient
overvoltage could be produced by load rejection or phase to
ground faults [5].
In the following, the overvoltage seen at the outlet side of
the GIS, due to load rejection at the medium voltage
distribution station connected to the investigated 380 kV GIS
is analyzed. The rejection of different inductive and capacitive
IPST 2015 CAVTAT CROATIA

Fig. 4. SFO variation with transmission line length.

Fig. 6. Delta shape 380 kV tower model
IV. FAST FRONT OVERVOLTAGE STUDY
One of the principal causes of fast front overvoltages in R0 1 E0  
RT  ; IC  
power systems is lightning stroke to transmission lines. I 2 R0 2 (3)
Lightning overvoltages have a wave head of several micro 1
IC
seconds and are one of the important factors to determine the
insulation design of substation equipment, especially in case of where: RT is the current dependent tower footing resistance,
GIS. Therefore, in the following the overvoltage seen by the R 0 is the low current and low frequency footing resistance
GIS substation at its entrance due to lightning strokes to the (10 Ω in our case), I C is the limiting current to initiate soil
380 kV OHL is investigated. To evaluate lightning
ionization, I is the stroke current,  is the soil resistivity
overvoltages both direct lightning and back flashover
situations are analyzed. (100 Ω·m) and E 0 is the soil ionization gradient (300 kV/m).
For insulation coordination purposes both situations are
investigated when the 380 kV OHL is connected to the GIS
substation with and without surge arresters mounted at GIS
gantry tower are investigated and obtained overvoltage values
are compared to GIS basic insulation level (BIL = 3.5 p.u.
considering a 1.1 p.u. safety margin). To model the 420 kV
surge arresters the IEEE Std. C62.22 [8] proposed frequency-
dependent model (see Fig. 7) has been implemented. Model
parameters have been computed based on surge arrester
Fig. 5. First five towers of the 380kV OHL modeled for FFO analysis.
datasheet provided by the manufacturer.
TABLE II
For the fast front overvoltage study the first five SURGE ARRESTER CHARACTERISTICS
transmission line towers from the GIS entrance were modeled System Rated Residual voltage
in detail. The two scenarios when lightning stroke hits the first voltage voltage 10 kA, 8/20 µs 10 kA, 1/5 µs
and respectively the third tower from GIS are taken into 420 kV 360 kV 783 kV 824 kV
consideration (see Fig. 5). A combined Hara [9] / Ametani
[10] multi-story model has been implemented for the 380 kV
delta shape towers to take into consideration both the bracings
and the damping effect of each tower section (see Fig. 6).
Tower footing is represented by a current dependent variable
resistance driven by equation (3) proposed by CIGRE [12] and
IEEE [4] to take into account soil ionization phenomena [13].
Fig. 7. Implemented IEEE model for 420 kV surge arrester.

A. Direct Lightning
Direct lightning strokes to phase wires in the case of
shielded transmission lines could appear in the case of
shielding failure. When a relatively low magnitude lightning
stroke bypasses the overhead ground (shield) wire and attaches
to one of the transmission lines phase conductors, shielding
failure occurs [13].
Direct lightning due to shielding failure could occur on the
IPST 2015 CAVTAT CROATIA
upper or the most outward phase wire. The maximum intensity GIS entrance for 40 kA and 20 kA lightning stroke currents
of a lightning stroke that could produce shielding failure can with and without surge arresters mounted at the GIS entrance:
be evaluated based on tower geometry, striking distance and
the implemented lightning attachment model [15]. For the
380 kV transmission line that connects the investigated GIS
substation to the grid, the Eriksson model [16] has been used
to evaluate the maximum shielding failure lightning current:

 
1 0.74
 R  R 2   2   2  1 
I MSF   
 
(4)
 0.67  H PhW 0.6   2  1 
 
where: I MSF is the maximum shielding failure current in (kA),
  H GrdW H PhW  ,  2  H GrdW  H PhW   R 2 , R is
0.6 2

the horizontal distance between phase and grounding wire, Fig. 9. Phase voltage at GIS entrance in case of direct lightning
H GrdW and H PhW are the height of the grounding and phase
wire respectively. B. Back Flashover
In order to determine the overvoltage at the GIS entrance
In order to evaluate the overvoltage seen at the GIS
due to direct lightning strokes to the 380 kV OHL phase
entrance due to back flashover on the 380 kV overhead
conductors, based on the evaluated maximum shielding failure
transmission line connected to the GIS, a 200 kA lightning
current values (25.6 kA for the upper phase and 20.4 kA for
current with a 8/20 µs waveform is considered. The black
the lower phase wires, respectively), lightning currents with
flashover mechanism has been modeled by an open switch in
1/5 µs waveform and an amplitude between 10 kA and 40 kA,
parallel to the phase wire insulator string capacitance. Once
are considered. Two different situations were analyzed: when
the voltage across the insulator string reaches the withstanding
the lightning stroke hits the upper phase wire (phase A see Fig.
capability of the insulator, back flashover occurs and the
6) in the vicinity of the first tower and when it hits nearby the
switch is closed.
third tower from GIS, respectively (see Fig. 5). For both
Due the fact that the insulator string may withstand a high
situations the worst case scenario is applied: lightning hits the
transient voltage for a short duration, but it could fail to
wire when the phase voltage reaches its positive peak and thus
withstand a lower transient voltage with a longer duration, the
the maximum overvoltage is produced in the 380 kV OHL.
volt-time characteristic proposed by CIGRE [12] has been
implemented for the back flashover simulations:
K
VWIns  K 1  0.752 (5)
t
where: VWIns is the insulator string withstand voltage in (kV),
K1  400  L , K 2  710  L , L is the length of the insulator
string in (m), and t is the elapsed time in (µs) from the lighting
stroke occurring.
Two different situations were investigated: when the
lightning stroke hits the ground wire at the first tower from the
GIS substation (Tower 1 see Fig. 5) and respectively when the
lightning hits the ground wire at the third tower (Tower 3). For
Fig. 8. Maximum overvoltage according to direct lightning current amplitude
both situations the phase overvoltage at the GIS entrance has
been evaluated with and without surge arrestors mounted at the
Fig. 8 presents the evaluated maximum overvoltage levels
380 kV OHL gantry tower. In order to obtain the highest
at the GIS entrance with and without surge arresters mounted
overvoltage values that could occur due to back flashover the
on the 380 kV OHL at GIS entrance. In the worst case
worst case scenario is applied: the lightning stroke hits the
scenario a 3.14 p.u. overvoltage is recorded for a 40 kA, 1/5
tower when the phase voltage on the upper conductor (phase
µs lightning current (which is almost double that the evaluated
A) reaches its negative peak. In this situation voltage across
maximum shielding failure current), and this could be reduced
the insulator string reaches its withstand capability earlier.
to 2.58 p.u. by placing surge arresters at the GIS entrance. For
Computed overvoltage waveforms are presented in Fig. 10:
less than 20 kA, respectively 25 kA direct lightning strokes
near the first tower and respectively the third tower or further
away the produced overvoltage levels at the GIS entrance will
be even lower than the discharge voltage of the surge arrester.
Fig. 9 represents the overvoltage waveform observed at the
IPST 2015 CAVTAT CROATIA
waveform recorded in the middle of the main bus during the
switching operation:

Fig. 11. GIS overvoltage produced by switching 380 kV TLine side CB.

In order to energize the outlet of the investigated GIS from

Fig. 10. Phase voltage at GIS entrance in the case of back flashover.
the 380kV transmission line the following sequence of closing
switching operations has been identified for the single bus GIS
It can be observed that the highest 4.97 p.u. overvoltage
configuration:
(1546 kV) occurs when the lightning stroke hits the first tower
(no surge arrester). If the lightning hits the third tower the OP01. Closing Load side DS1, DS2 and CB - Closing 380
overvoltage impulse needs a 2 µs travel time to reach the GIS kV TLine side DS1, DS2 and CB;
entrance and it is reduced to 3.44 p.u. (1069 kV). By mounting OP02. Closing 380 kV TLine side DS1, DS2 and CB -
surge arresters at the GIS entrance these overvoltage values Closing Load side DS1, DS2 and CB;
will be reduced to 2.9 p.u. and 2.3 p.u. respectively, lower than
the basic insulation level of the GIS (3 p.u.).

V. VERY FAST TRANSIENT OVERVOLTAGE STUDY

Very fast transient overvoltage could be generated in GIS
substations during the opening and closing of circuit breakers
or disconnectors. The disconnector restriking surge is an Fig. 12. Maximum overvoltage / number of energized switching for the one
oscillation surge with a very high frequency of several MHz. bus GIS configuration.
The frequency of the disconnector restriking surge is much Fig. 12 presents the produced maximum overvoltage levels
higher than that of lightning surges because every disconnector inside the one bus GIS and the number of energized switching
operation potentially generates the overvoltage, and the surge need by each of the above presented closing operation
could impose negative effects not only on the main circuit sequences.
insulation but also on a secondary system such as EMC [5]. B. Double Bus Configuration
Therefore, all the possible disconnector and circuit breaker
If Bus B of the double bus GIS configuration is totally
closing switching operations have been analysed for the three
disconnected from the 380 kV TLine feeder, Load feeder and
investigated GIS configurations (single bus Fig. 2a, double bus
respect Bus A (see Fig. 2b) than the double bus GIS will work
Fig. 2b and one and a half bus Fig. 2c). It is considered that the
in a one bus configuration similar to that presented above. In
switching order for a feeder connected to the main bus is
this case the overvoltage introduced by switching the 380 kV
closing the disconnectors and then the circuit breaker.
TLine feeder CB is reduced to 2.31 p.u. Fig. 13 presents the
During the closing operation of a disconnector (DS) or
overvoltage waveform recorded at Bus A, measurement point
circuit breaker (CB), the sparks are modeled by a fixed
2 (in the middle of Bus Duct A see Fig. 2b):
resistance in series with an exponentially decreasing one:
 t
Rt   R0  R1  exp   (6)
 
where: R0  0.5 , R1  10 6  and   1ns .
Obtained overvoltage values due to DS and CB switching
were analyzed at several measurement points along the
investigated GIS configurations (see Fig. 2).
A. Single Bus Configuration
Fig. 13. GIS overvoltage produced by switching 380 kV TLine side CB.
Analyzing all the possible disconnector and circuit breaker
switching inside the single bus configuration GIS revealed that
In order to energize the outlet of the investigated GIS from
the highest VFTO value is produced by switching (closing) the
the 380kV transmission line the following sequence of closing
380 kV TLine feeder circuit breaker with the load side feeder
switching operations has been identified for the double bus
connected to the bus duct. Fig. 11 presents the overvoltage
IPST 2015 CAVTAT CROATIA
GIS configuration (see Fig. 2b): OP04. Closing 380 kV TLine DS - Closing 380 kV TLine to
OP01. Closing Load side DS1, DS2 and CB - Closing 380 Load DS and CB - Closing Load side DS with both bus bars
kV TLine side DS1, DS2 and CB (energizing the load side disconnected;
through Bus A with Bus B disconnected); OP05. Closing Load side DS - Closing 380 kV TLine to
OP02. Closing 380 kV TLine side DS1, DS2 and CB - Load DS and CB - Closing 380 kV TLine DS with both bus
Closing Load side DS1, DS2 and CB (energizing the load side bars already connected;
through Bus A with Bus B disconnected); OP06. Closing 380 kV TLine to Load DS and CB - Closing
380 kV TLine DS - Closing Load side DS with bus bars
OP03. Closing Load side DS1, DS3 and CB - Closing 380
already connected;
kV TLine side DS1, DS3 and CB (energizing the load side
OP07. Closing Load side DS - Closing 380 kV TLine DS -
through Bus B with Bus A disconnected);
Closing 380 kV TLine to Load DS and CB with both bus bars
OP04. Closing 380 kV TLine side DS1, DS3 and CB -
already connected;
Closing Load side DS1, DS3 and CB (energizing the load side
OP08. Closing 380 kV TLine DS - Closing 380 kV TLine to
through Bus B with Bus A disconnected);
Load DS and CB - Closing Load side DS with both bus bars
OP05. Closing Load side DS1, DS2 and CB - Closing 380 already connected;
kV TLine side DS1, DS2 and CB (energizing the load side
through Bus A with Bus B already connected to Bus A);
OP06. Closing 380 kV TLine side DS1, DS2 and CB -
Closing Load side DS1, DS2 and CB (energizing the load side
through Bus A with Bus B already connected to Bus A);
OP07. Closing 380 kV TLine side DS1, DS3 and CB -
Closing Load side DS1, DS3 and CB (energizing the load side
through Bus A with Bus B already connected to Bus B);
OP08. Closing Load side DS1, DS3 and CB - Closing 380
Fig. 15. Maximum overvoltage / number of energized switching for the one
kV TLine side DS1, DS3 and CB (energizing the load side and a half bus GIS configuration.
through Bus B with Bus A already connected to Bus B); From Fig. 15 it can be observed that closing sequence
OP04 produces the lowest overvoltage level in the one and
half bus GIS configuration (1.85 p.u.) but it needs 4 energized
switching, while sequence OP1 needs only one energized
switching with a maximum overvoltage of 2.23 p.u.

VI. CONCLUSIONS
Overvoltages produced by the rejection of inductive loads
decrease slightly with the increase of the inductive part of the
Fig. 14. Maximum overvoltage / number of energized switching for the
double bus GIS configuration.
rejected loads, while overvoltages produced by the rejection of
the capacitive load exhibits a ‘V’ curve.
From Fig. 14 it can be observed that closing sequence Different line length cables and/or overheads connecting
OP06 produces the lowest overvoltage level in the double bus the medium voltage distribution station to the GIS have
GIS configuration (1.76 p.u.) but it needs 4 energized revealed that the overvoltage produced by the overhead line is
switchings, while sequence OP1 needs only two energized greater.
switching with a maximum overvoltage of 2.31 p.u. The combined Hara [9] and Ametani [10] multi-story
tower model considers both the bracings and the damping
C. One and a Half Bus Configuration effect and used in conjunction with an adaption of CIGRE
In the case of the one and a half bus GIS configuration (Fig. tower footing variable resistance yields realistic results.
2c) the following sequence of closing switching operations has Results from all possible disconnectors and circuit breakers
been identified: closing switching operations for the three GIS configurations
OP01. Closing Load side DS - Closing 380 kV TLine to show that the switching sequence plays a significant role in the
Load DS and CB - Closing 380 kV TLine DS with both bus magnitude of the generated overvoltage. For example, for
bars disconnected; single bus, double bus, one and a half bus GIS configurations,
OP02. Closing 380 kV TLine to Load DS and CB - Closing the switching operations OP02, OP06 and OP04, respectively,
380 kV TLine DS - Closing Load side DS with both bus bars generate the minimum overvoltage.
disconnected;
OP03. Closing Load side DS - Closing 380 kV TLine DS - VII. REFERENCES
Closing 380 kV TLine to Load DS and CB with both bus bars [1] A. Tavakoli, A. Gholami, H. Nouri, M. Negnevitsky, “Comparison
disconnected; Between Suppressing Approaches of Very Fast Transients in Gas-
IPST 2015 CAVTAT CROATIA
Insulated Substations (GIS)”, IEEE Trans. on Power Delivery, vol. 28,
pp. 303-310, January 2013.
[2] G. H. C. Oliveira, S. D. Mitchell, “Comparison of Black-Box Modeling
Approaches for Transient Analysis: A GIS Substation Case Study”, in
Proc. International Conference on Power System Transients, (IPST),
Vancouver, Canada, 18-20 July, 2013
[3] M. M. Rao, M. J. Thomas, B. P. Singh, “Frequency Characteristics of
Very Fast Transient Currents in a 245-kV GIS”, IEEE Trans. on Power
Delivery, vol. 20, pp. 2450-2457, October 2005
[4] IEEE Fast Front Transients Task Force, “Modeling Guidelines for Fast
Front Transients”, IEEE Trans. on Power Delivery, vol. 11, pp.
493-506, January 1996.
[5] IEEE Std. 1313.2, “IEEE Guide for the Application of Insulation
Coordination”, IEEE-SA Standards Board, 1999.
[6] IEC 60071-2, “Insulation Coordination - Part 2: Appication Guide”,
International Electrotechnical Commission, (IEC), Geneva, Swizerland,
1996.
[7] H. C. Seo, W. H. Jang, C. H. Kim, T. Funabashi, T. Senju, “Analysis of
rate-of-rise of VFTO according to Switching Conditions in GIS”, in
Proc. International Conference on Power System Transients, (IPST),
Delft, The Netherlands, 14-17 June, 2011
[8] IEEE Std. C62.22, “IEEE Guide for the Application of Metal-Oxide
Surge Arresters for Alternating-Current Systems”, IEEE-SA Standards
Board, New York, USA, 2009.
[9] T. Hara, O. Yamamoto, “Modelling of a Transmission Tower for
Lightning Surge Analysis”, IEE Proc. - Generation, Transmission and
Distribution, vol. 143, pp. 283-289, May 1996.
[10] A. Ametani, T. Kawamura, “A Method of a Lightning Surge Analysis
Recommended in Japan Using EMTP”, IEEE Trans. on Power
Delivery, vol. 20, pp. 867-895, April 2005.
[11] A. Ceclan, V. Topa, D.D. Micu, A. Andreotti, “Lightning-Inverse
Reconstruction by Remote Sensing and Numerical-Field Synthesis”,
IEEE Trans. on Magnetics, vol. 49 , no. 5, pp. 1657-1660, May 2013.
[12] CIGRE WG33.01, “Guide to Procedures for Estimating the Lightning
Performance of Transmission Lines”, CIGRE Technical Report 063,
October 1991.
[13] Z. G. Datsios, P. N. Mikropoulos, T. E. Tsovilis, “Impulse Resistance of
Concentrated Tower Grounding Systems Simulated by an ATPDraw
Object”, in Proc. International Conference on Power System
Transients, (IPST), Delft, Netherlands, 2011.
[14] P. N. Mikropoulos, T. E. Tsovilis, “Lightning Attachment Models and
Maximum Shielding Failure Current of Overhead Transmission Lines:
Implications in Insulation Coordination of Substations”, IET
Generation, Transmission & Distribution, vol. 4, pp. 1299-1313,
December, 2010.
[15] IEEE Std. 1243, “Guide for Improving the Lightning Performance of
Transmission Lines”, June, 1997.
[16] A. J. Eriksson, “An improved electrogeometric model for transmission
line shielding analysis”, IEEE Trans. Power Delivery, vol. 2, no. 3, pp.
871–886, 1987.
[17] J. Mahseredjiana, S. Dennetièreb, L. Dubéc, B. Khodabakhchiand, L.
Gérin-Lajoiee, “On a New Approach for the Simulation of Transients in
Power Systems”, Electric Power Systems Research, vol. 77, pp. 1514 –
1520, September 2007.
[18] L. Zhang, Q. Zhang, S. Liu, F. Liu, “Insulation Characteristics of 1100
kV GIS under Very Fast Transient Overvoltage and Lightning Impulse”,
IEEE Transactions on Dielectrics and Electrical Insulation, vol. 19, no.
3, pp. 1029-1036, June 2012.