Very Fast Transient Overvoltages (VFTO) in SF6 Gas-Insulated Substations

Phenomenology of Very Fast Front Transients (VFT) in SF6 Gas-Insulated Substations (GIS)

In sulfur hexafluoride (SF6) gas-insulated substations (GIS), disconnector switch (DS) operations represent one of the most severe and complex sources of electromagnetic disturbance within the power system. Unlike high-voltage circuit breakers (CB), disconnector switches are designed to operate at relatively slow speeds, with contact translation velocities typically ranging between 0.1 m/s and 1.0 m/s. This low mechanical speed, combined with the progressively decreasing inter-electrode dielectric distance during opening or closing sequences, inevitably leads to the occurrence of multiple pre-strikes and re-strikes across the SF6 gas gap.

When the disconnector switch closes, the dielectric strength of the gas gap between the moving and stationary contacts decreases continuously until the local electric field exceeds the critical breakdown strength of the SF6 gas. At this precise instant, an abrupt dielectric breakdown (voltage collapse) occurs within an extremely short time interval, typically ranging from 3 ns to 20 ns. This rapid voltage collapse acts as a near-vertical step-voltage wave that propagates in both directions from the electric arc point along the coaxial busbars of the GIS.

The resulting traveling voltage wave is characterized by ultra-high-frequency spectral components, spanning from a few hundred kilohertz (kHz) to well over 100 megahertz (MHz). These electromagnetic phenomena are defined as Very Fast Front Transients (VFT). When analyzed specifically from the perspective of the resulting overvoltage amplitude, they are designated as Very Fast Front Overvoltages (VFTO).

The propagation of these traveling waves within the GIS is governed by coaxial transmission line theory. Because the GIS is a highly compact system with extremely low propagation losses in the high-frequency spectrum, the traveling waves undergo multiple reflections and refractions at every physical and impedance discontinuity within the substation. These discontinuities include elbows, epoxy resin insulator spacers, T-junctions, instrument transformers (CTs/VTs), cable terminations, and SF6-to-air outdoor bushings. The constructive superposition of these reflected waves gives rise to highly complex waveforms with local overvoltages that can significantly exceed the nominal system voltage, reaching peak values of up to 2.5 p.u. (per unit) of the maximum system crest voltage. These overvoltages exhibit extremely high voltage gradients ( $dV/dt$ ), posing a severe threat to the dielectric integrity of adjacent primary equipment insulation.

SF6 Electric Arc Dynamics during Disconnector Operation

The dielectric breakdown process in SF6 during a disconnector operation is not a single, isolated event, but rather a quasi-periodic sequence of electrical discharges. During a closing operation, as the physical distance between the contacts decreases, the voltage withstand capability of the gas gap degrades. When the instantaneous power-frequency voltage across the contacts exceeds the breakdown voltage of the gas at that specific gap distance, a highly conductive plasma channel is established.

The voltage collapse time during the formation of the spark channel ( $\tau_f$ ) can be modeled using semi-empirical approximations based on plasma physics, where the SF6 gas pressure and the local electric field intensity are the dominant variables. The rapid voltage drop across the spark gap generates an electromagnetic wave that propagates through the coaxial GIS enclosure. Following the first pre-strike, the parasitic capacitances of the GIS busbars charge or discharge through the arc channel until the high-frequency transient current is extinguished—either due to potential equalization across the contacts or when the high-frequency current passes through zero.

However, because the mechanical contacts continue their slow physical movement, the power-frequency voltage (50/60 Hz) on the source side continues to oscillate, while the disconnected busbar section retains a trapped charge (resulting in a "trapped charge voltage"). This creates a renewed potential difference across the contacts. A second pre-strike is triggered as soon as the local electric field again exceeds the dielectric strength of the gas, which now occurs at a shorter inter-electrode distance. This cycle repeats sequentially—potentially dozens of times during a single mechanical operation—until the physical contacts make solid, low-resistance metallic contact.

Mathematical Modeling and Wave Propagation Physics in Coaxial Media

To analyze the propagation of VFT within a gas-insulated substation with scientific rigor, the GIS must be treated as a multi-phase system of coupled coaxial transmission lines. The inner conductor (typically made of aluminum or copper) and the outer metal enclosure (aluminum or stainless steel) form a coaxial waveguide where the dielectric medium is the SF6 gas ( $\epsilon_r \approx 1.002$ at standard operating pressures).

The characteristic impedance ( $Z_c$ ) of a single-phase coaxial GIS busbar is defined by the geometry of its conductors according to the following analytical expression:

Z_c = \frac{\eta_0}{2\pi \sqrt{\epsilon_r}} \ln\left(\frac{b}{a}\right) \approx \frac{60}{\sqrt{\epsilon_r}} \ln\left(\frac{b}{a}\right)

Where:

$\eta_0 \approx 376.73\ \Omega$ is the intrinsic impedance of free space.
$\epsilon_r$ is the relative permittivity of the SF6 gas at the operating pressure.
$b$ is the internal radius of the outer GIS enclosure or shield.
$a$ is the external radius of the inner cylindrical conductor.

Typically, the characteristic impedances of GIS busbars range between $60\ \Omega$ and $80\ \Omega$ . These values are significantly lower than the characteristic impedances of high-voltage overhead lines (which range from $300\ \Omega$ to $450\ \Omega$ ) and higher than those of XLPE-insulated underground power cables (which typically range between $15\ \Omega$ and $30\ \Omega$ ).

Telegrapher's Equations and Dispersive Frequency Effects

The propagation of the transient voltage wave $v(x,t)$ and current wave $i(x,t)$ along the coaxial structure is modeled using the time-domain partial differential Telegrapher's equations:

-\frac{\partial v(x,t)}{\partial x} = R' \cdot i(x,t) + L' \frac{\partial i(x,t)}{\partial t}

-\frac{\partial i(x,t)}{\partial x} = G' \cdot v(x,t) + C' \frac{\partial v(x,t)}{\partial t}

Where $R'$ , $L'$ , $G'$ , and $C'$ represent the distributed parameters per unit length of the coaxial GIS busbar:

$L' = \frac{\mu_0 \mu_r}{2\pi} \ln\left(\frac{b}{a}\right)$ is the inductance per unit length.
$C' = \frac{2\pi \epsilon_0 \epsilon_r}{\ln\left(\frac{b}{a}\right)}$ is the capacitance per unit length.
$R'$ is the resistance per unit length, which is highly frequency-dependent due to the skin effect in both the inner conductor and the outer enclosure.
$G'$ is the conductance per unit length, representing the dielectric losses in the epoxy resin spacers (generally negligible in the SF6 gas itself, but critical within the solid supports at MHz frequencies).

At the ultra-high frequencies characteristic of VFT (up to 100 MHz), the skin effect restricts current flow to a very thin surface layer of the conductors. The current penetration depth or skin depth ( $\delta$ ) is calculated as:

\delta = \sqrt{\frac{1}{\pi f \mu \sigma}}

Where $f$ is the frequency of the transient component, $\mu$ is the magnetic permeability of the conductor material, and $\sigma$ is its electrical conductivity. For aluminum at 100 MHz, $\delta$ is on the order of a few micrometers. This restriction drastically increases the distributed resistance parameter $R'$ , causing frequency-dependent attenuation of the ultra-high-frequency components as the wave propagates. However, this attenuation is generally insufficient to suppress the initial, highly damaging peak voltage crests in the immediate vicinity of the disconnector switch.

Modeling of Disconnector Arc Resistance

The spark channel established between the disconnector contacts exhibits a time-varying resistance, $R_g(t)$ , whose dynamic behavior is a critical factor in determining the rise time of the VFT wave. Various mathematical models exist to represent this transition from an insulating state to a highly conductive state. One of the most widely accepted physical models is Toepler's Spark Law:

R_g(t) = \frac{k_T \cdot d}{\int_0^t i(\tau) d\tau}

Where:

$k_T$ is Toepler's constant for SF6 gas (typically in the range of $0.5 \times 10^{-4}$ to $1 \times 10^{-4}\ V \cdot s / m$ ).
$d$ is the instantaneous spark gap distance between the contacts.
$i(\tau)$ is the transient current flowing through the arc channel.

Alternatively, the Rompe-Weizel model describes the conductivity of the plasma channel based on the thermal energy balance of the arc column:

R_g(t) = \frac{a_R \cdot d}{\left[ \int_0^t i^2(\tau) d\tau \right]^{2/3}}

Where $a_R$ is an empirical coefficient characteristic of SF6 gas. These models demonstrate that the rate of decay of the arc resistance directly dictates the front time of the initial electromagnetic disturbance; higher SF6 gas pressures and smaller dielectric gap distances yield shorter front times ( $t_f$ ), thereby increasing the high-frequency harmonic content of the electromagnetic transient.

Wave Propagation and Reflection Phenomena at Discontinuities

When a VFT traveling wave of amplitude $V_i$ propagates along the coaxial GIS busbar and encounters an impedance discontinuity (such as a junction with another power system component), reflected ( $V_r$ ) and transmitted ( $V_t$ ) waves are generated.

The amplitudes of these waves are calculated using the reflection ( $\Gamma$ ) and transmission ( $T$ ) coefficients at the junction node between a zone with characteristic impedance $Z_1$ and a zone with $Z_2$ :

\Gamma = \frac{Z_2 - Z_1}{Z_2 + Z_1}

T = \frac{2 Z_2}{Z_2 + Z_1} = 1 + \Gamma

Let us analyze the critical impedance transition scenarios that occur within a GIS substation:

Transition from GIS Busbar to Power Transformer

The high-frequency input impedance of a power transformer at MHz frequencies behaves primarily as a complex capacitive network. This network is dominated by the input capacitance of the high-voltage bushings and the series-parallel capacitance of the windings. Consequently, the transformer presents a very high equivalent impedance, which is modeled in high-frequency transient studies as an equivalent open circuit ( $Z_2 \gg Z_1$ ).

Approximating $Z_2 \to \infty$ :

\Gamma \approx 1, \quad T \approx 2

This phenomenon results in a theoretical doubling of the incident transient voltage wave amplitude at the transformer terminals. The reflected wave superimposes constructively with the incoming wave, generating extreme dielectric stress on the bushing insulation and the initial turns of the transformer winding.

Transition from GIS Busbar to High-Voltage Underground Cable

When the GIS is connected to an XLPE-insulated underground power cable, the transition occurs from a higher impedance $Z_1 \approx 70\ \Omega$ to a much lower cable impedance, $Z_2 \approx 20\ \Omega$ . In this scenario:

\Gamma = \frac{20 - 70}{20 + 70} = -0.55

T = \frac{2 \times 20}{20 + 70} = 0.44

The voltage wave transmitted into the cable is reduced to 44% of the incident wave amplitude, and a reflected wave of inverted polarity is generated. While this transition mitigates the transient amplitude within the cable itself, the negative reflected wave travels back into the GIS. There, it can superimpose with other positive reflections from open disconnectors (which act as open circuits with $\Gamma = 1$ ), leading to complex resonance phenomena and localized voltage amplification within short GIS busbar sections.

Forensic Analysis of Engineering Failures Induced by VFTO

The severity of VFTO is not solely a function of its peak amplitude, but is heavily driven by its extremely high rate of voltage rise ( $dV/dt$ ), which can exceed $1000\ kV /\mu s$ (or $1\ kV / ns$ ). This section details the physical failure mechanisms and dielectric degradation processes in critical substation components.

Power and Instrument Transformers

Power transformers connected directly to GIS substations are subjected to severe dielectric stresses for which conventional insulation systems are not designed. Under power-frequency voltages or standard lightning impulses (with a front time of $1.2\ \mu s$ ), the voltage distribution along the transformer winding is predominantly linear, governed by the inductive coupling of the coils.

However, when subjected to a very fast front transient (with a front time of $10\ ns$ ), the winding behavior is governed almost exclusively by its parasitic capacitance network, consisting of inter-turn capacitances ( $C_s$ ) and capacitances to ground ( $C_g$ ). The initial voltage distribution along the winding is determined by the capacitive distribution factor $\alpha$ :

\alpha = \sqrt{\frac{C_g}{C_s}}

At ultra-high frequencies, $\alpha$ is very large, resulting in an extremely non-linear initial voltage distribution. Virtually the entire transient voltage drop is concentrated across the first 5% to 10% of the entry winding turns adjacent to the bushing. This creates an unsustainable inter-turn potential gradient, initiating localized surface partial discharges and, eventually, puncturing the paper-oil or solid resin insulation, triggering a catastrophic inter-turn fault.

Furthermore, if the fundamental frequency of the VFTO spectrum coincides with one of the natural resonant frequencies of the transformer winding (internal resonance), internal voltage amplification occurs. This can destroy the insulation in deep or intermediate sections of the winding—a failure mode that is exceptionally difficult to detect via routine field maintenance testing.

Power Cable Terminations and Joints

Power cables interfaced with GIS experience severe dielectric stress at their transition terminations (cable bushings). The geometric and permittivity discontinuities at the interface between the cable insulation (XLPE), the epoxy resin of the bushing, and the SF6 gas distort the local electric field distribution.

The high $dV/dt$ rate of the VFTO dramatically increases the displacement current density through these dielectric interfaces:

J_d = \epsilon \frac{\partial E}{\partial t}

This displacement current leads to localized dielectric loss heating and concentrates electrical stress at the edge of the stress cone. If the stress relief system is not optimized for high-frequency performance, surface partial discharges are initiated along the XLPE/epoxy interface, progressively degrading the material until a solid phase-to-ground flashover occurs.

GIS Internal Insulation: Spacer Degradation and Metallic Particle Dynamics

Epoxy resin insulator spacers are the primary mechanical components responsible for centering and supporting the inner conductor relative to the outer enclosure. Although SF6 gas possesses excellent dielectric strength, the presence of microscopic metallic contaminants (introduced during site assembly or generated by mechanical wear from disconnector operations) severely distorts the local electric field.

Under the high-frequency electric fields associated with VFTO, these metallic particles experience rapid electromagnetic levitation and oscillation. The VFTO can trigger gas ionization at the sharp tips of these particles, initiating pre-discharges that propagate along the surface of the epoxy spacer. This phenomenon, known as spacer surface flashover, is further exacerbated by the accumulation of static charges on the epoxy surface, which can reduce the dielectric strength of the spacer support by up to 50% during disconnector switching.

Transient Ground Potential Rise (TGPR) and Control Systems

One of the most critical secondary effects of VFTO is Transient Ground Potential Rise (TGPR). When the VFT wave reaches the transition from the shielded GIS enclosure to the external grounding system (for example, at SF6-to-air bushings), the high-frequency transient current must return to ground.

Because conventional grounding conductors possess non-negligible inductance (typically $1\ \mu H / m$ ), the impedance of a standard 2-meter grounding strap at 50 MHz is extremely high:

X_L = 2\pi f L = 2\pi \times (50 \times 10^6\ Hz ) \times (2 \times 10^{-6}\ H ) \approx 628\ \Omega

This high impedance prevents the transient current from flowing freely to the deep grounding grid, forcing it to seek alternative low-impedance paths. Consequently, the external metallic enclosure of the GIS experiences a transient potential rise, with amplitudes reaching tens of kilovolts (kV) relative to remote ground.

This TGPR phenomenon induces severe inductive and capacitive coupling onto control, instrumentation, and auxiliary power cables running parallel to the GIS. The induced transient currents penetrate Local Control Panels (LCP), resulting in:

Physical destruction of RS-485, Ethernet, or fiber-optic transceivers with metallic connectors.
Erratic tripping of digital protection relays due to electromagnetic noise injection into analog current and voltage inputs.
Insulation breakdown in the secondary windings of instrument transformers (CTs and VTs).

Technical Parameter	Typical Range / Standard Limit	Reference Standard	Consequence of Exceedance
Front Time ( $t_f$ )	3 ns to 100 ns	IEC 60071-1 / IEC 62271-102	Highly non-linear voltage distribution across transformer windings.
Maximum VFTO Amplitude	1.5 p.u. to 2.5 p.u.	IEEE C37.122 / IEC 62271-203	Dielectric breakdown of SF6 and surface flashover of epoxy spacers.
Rate of Rise ( $dV/dt$ )	Up to 2000 kV/µs	CIGRE WG A3.22	High displacement currents, degradation of cable termination interfaces.
Enclosure TGPR Amplitude	5 kV to 30 kV	IEEE 80 / IEC 61000-4-5	Permanent damage to secondary control equipment and electric shock hazards.
Frequency Spectrum	100 kHz to 100 MHz	IEC 61000-4-25 / Cigre Brochure 519	Internal electromagnetic resonance and severe radiated EMI.

Mitigation Strategies and Electrical Design Criteria

Mitigating the effects of VFTO requires a coordinated implementation of design countermeasures across both the primary high-voltage circuit and the secondary grounding and shielding systems.

Pre-Insertion Resistors (PIR) in Disconnector Switches

The most effective mitigation technique at the source is to equip disconnector switches with Pre-Insertion Resistors (PIR). The operating principle involves temporarily inserting a damping resistor in series with the arc circuit before the main contacts establish solid metallic contact.

When the disconnector closes, the auxiliary contact closes first, forcing the initial pre-strikes to occur through the resistor $R_{ PIR }$ . The optimal value of this resistor is selected to approximate the characteristic impedance of the GIS busbar:

R_{ PIR } \approx Z_c \approx 60\ \Omega\ to \ 100\ \Omega

By matching the circuit resistance to the characteristic impedance, the reflection coefficient at the arc point is drastically reduced, damping the traveling wave at its source and reducing the rate of rise ( $dV/dt$ ) by more than 70%. However, the inclusion of PIR significantly increases the mechanical complexity, cost, and maintenance requirements of the disconnector switch, which typically limits its application to Ultra-High Voltage (UHV) substations rated at 500 kV and above.

Ferrite Rings (Magnetic Rings)

A highly reliable, non-invasive passive alternative is the installation of nanocrystalline ferrite rings around the inner conductor of the GIS at strategic locations (typically near the disconnector terminals).

The ferrite acts as a frequency-dependent impedance. At power frequency (50/60 Hz), its relative magnetic permeability is low to avoid hysteresis and induction losses, behaving essentially as a lossless conductor. However, at MHz frequencies (the range of VFTO), the ferrite exhibits a complex, high-loss permeability:

\mu(f) = \mu'(f) - j\mu''(f)

The imaginary component $\mu''(f)$ represents magnetic losses within the core. At high frequencies, this introduces a high equivalent series resistance that absorbs the electromagnetic energy of the high-frequency transient components, attenuating the wave amplitude and smoothing the front rise time ( $dV/dt$ ) without affecting normal power-frequency operation.

RC Snubbers and Zinc Oxide (ZnO) Surge Arresters

Conventional zinc oxide (ZnO) surge arresters installed in the GIS offer limited protection against VFTO. This limitation arises because the intrinsic response time of the ZnO ceramic material and the parasitic inductance of the internal connection leads restrict their effectiveness for wave fronts faster than 50 ns.

To address this, high-frequency ZnO surge arresters with direct capacitive coupling or ultra-low parasitic inductance (less than 10 nH) are utilized. Additionally, installing pure capacitive filters (high-frequency SF6-insulated damping capacitors rated at $1\ nF$ to $5\ nF$ ) near power transformers diverts high-frequency spectral components to ground, increasing the front time of the traveling wave and protecting the windings.

High-Frequency Grounding System Optimization

To mitigate the TGPR phenomenon and protect secondary control systems, high-frequency grounding design principles must be applied:

Low-Inductance Grounding Connections: Replace conventional round copper conductors with wide, flat copper straps. Flat straps exhibit significantly lower internal and external inductance at high frequencies due to the increased effective perimeter available for skin effect current flow.
Multi-Point Enclosure Interconnection: The metallic enclosure of the GIS must be connected to the structural ground grid at multiple, closely spaced points, particularly at physical discontinuities (bushings, current transformers, cable terminations). This creates parallel paths, reducing the equivalent grounding inductance:
$L_{ eq } = \frac{1}{\sum \frac{1}{L_i}}$
Coaxial Enclosure Grounding Rings: At the insulating flanges of cable terminations, fast-acting gas discharge tubes or low-voltage varistors should be installed in parallel with the flange. This transiently short-circuits the discontinuity during high-frequency events, maintaining a continuous shield for the traveling wave.
Control Cable Shielding: All control and signal cables routing from the GIS to the LCP must feature double metallic shielding (high-coverage copper braid). The shields must be grounded at both ends using 360° concentric EMC cable glands to ensure the continuity of the Faraday cage.

Advanced Modeling and Analysis using Vexten Suite

The analysis of very fast front electromagnetic transients (VFT) requires high-precision digital simulation tools capable of modeling the complex physics of multi-phase transmission lines, electric arc dynamics, and the high-frequency behavior of substation equipment. The Vexten Suite platform integrates specialized modules that enable design engineers to perform these studies with rigorous accuracy and automation.

Network Initialization via Short-Circuit Analysis (IEC 60909 / IEEE 141)

Electromagnetic transient analysis cannot be performed in isolation from the short-circuit parameters of the surrounding power grid. The Vexten Short-Circuit module calculates the sequence impedances of the network at the Point of Common Coupling (PCC) of the GIS substation under the IEC 60909 and IEEE 141 standards.

These calculated impedances establish the initial conditions for the transient simulation, including the initial symmetrical short-circuit current ( $I''_k$ ), the reactance-to-resistance ratio ( $X/R$ ), and the short-circuit capacity of the source ( $S''_k$ ). These parameters define the dynamic strength of the grid and are used to initialize the AC source voltage prior to disconnector switching, ensuring that the transient model represents the worst-case operating scenario (e.g., switching under peak load or grid contingency conditions).

Power Cable Modeling and Skin Effect Analysis (IEC 60287 / NEC 310)

When VFTO waves propagate into high-voltage underground cables, the thermal and impedance behavior of the cable at extreme frequencies must be modeled with high fidelity. The Vexten Thermal & Cable module implements the calculation algorithms of the IEC 60287 standard to determine cable ampacity and the distribution of losses within metallic shields and armors.

At MHz frequencies, the skin and proximity effects drastically alter the AC resistance of the cable conductor ( $R_{ AC }$ ). The software calculates the high-frequency thermal derating factor by numerically solving the current diffusion equations:

R_{ AC }(f) = R_{ DC } \left( 1 + y_s(f) + y_p(f) \right)

Where $y_s(f)$ and $y_p(f)$ are the skin and proximity effect factors, respectively, which are dynamically recalculated by Vexten Suite across the entire harmonic spectrum of the VFTO. This enables precise sizing of the insulation thickness and shield design under the NEC 310 standard, preventing localized hot spots caused by repetitive switching transients or high-frequency resonances.

Resonance Mitigation and Filter Optimization (Vexten Resonance Mitigation)

The resonance mitigation module within Vexten Suite performs high-resolution frequency scans of the GIS busbar system and connected primary equipment. The software calculates the nodal impedance matrix $Z(f)$ as seen from the disconnector switch contacts:

[V(f)] = [Z(f)] \cdot [I(f)]

If the impedance profile $Z(f)$ exhibits sharp peaks (resonant poles) that coincide with the dominant frequencies of the VFTO calculated by the transient engine, the software alerts the design engineer to the risk of internal resonance.

Using its built-in optimization algorithms, Vexten Suite automatically generates the required technical specifications for mitigation equipment:

Ferrite Rings: Determines the required volume of ferrite, the optimal complex magnetic permeability $\mu'(f) - j\mu''(f)$ , and the ideal physical location along the conductor to maximize wave-front damping.
RC Snubbers: Calculates the optimal capacitance ( $C_{ snub }$ ) and resistance ( $R_{ snub }$ ) values to tune the damping filter, reducing resonant peak amplitudes below the dielectric withstand limits of the transformer.
Pre-Insertion Resistors (PIR): Evaluates the impact of various $R_{ PIR }$ values on reducing the voltage gradient ( $dV/dt$ ) and peak overvoltage at the power transformer terminals.

This integrated engineering approach ensures that the design of the gas-insulated substation meets the highest standards of reliability, safety, and electromagnetic compatibility, minimizing the risk of catastrophic failures induced by routine operational switching.