# Analyzing Ferroresonance Phenomena in Power Transformers Including Zinc Oxide Arrester and Neutral Resistance Effect

Authors: Hamid Radmanesh; Fathi Seyed Hamid
Journal: Applied Computational Intelligence and Soft Computing
 (2012)

---

## Abstract

This paper studies the effect of zinc oxide arrester (ZnO) and neutral earth resistance on controlling nonconventional oscillations of the unloaded power transformer. At first, ferroresonance overvoltage in the power system including ZnO is investigated. It is shown this nonlinear resistance can limit the ferroresonance oscillations but it cannot successfully control these phenomena. Because of the temperature dissipation of ZnO, it can withstand against overvoltage in a short period and after that ferroresonance causes ZnO failure. By applying neutral earth resistance to the system configuration, mitigating ferroresonance has been increased and chaotic overvoltage has been changed to the smoother behavior such as fundamental resonance and periodic oscillation. The simulation results show that connecting the neutral resistance exhibits a great mitigating effect on nonlinear overvoltage.


---

## Body

## 1. Introduction

Ferroresonance is a complex electromagnetic phenomenon which may be neglected in power system studies which is carried out for routine designs, planning, and operations [1]. A stability domain calculation of period-1 ferroresonance in nonlinear resonant circuit power system elements is given in [2]. In this case, quasistatic analytical approaches can be used to give a quick indication of the locations of domains of different ferroresonant states as a function of a set of parameters. Fast ferroresonance suppression of coupling capacitor voltage transformers (CCVT) is studied in [3]. This paper describes a procedure for fast suppression of the phenomenon of ferroresonance in CCVT without major change in the design. The design of a hall effect current transformer and examination of the linearity with real time parameter estimation is given in [4]. The aim of “blind source separation” (BSS) is to recover mutually independent unknown source signals only from observations obtained through an unknown linear mixture system [5]. Sensitivity studies on power transformer ferroresonance of a 400 kV double circuit are given in [6]. Novel analytical solution to fundamental ferroresonance in [7] investigated a major problem with the traditional excitation characteristic of nonlinear inductors. Application of wavelet transform and MLP neural network for ferroresonance identification was done in [8]. Impacts of transformer core hysteresis formation on stability domain of ferroresonance modes were done in [9]. The principle of AC current transformers (CT) based on the magnetic coupling principle is given in [10]. Current paper studies the effect of neutral resistance on the global behavior of a ferroresonance circuit including ZnO in parallel to the transformer with linear core losses. Then, neutral earth resistance effect is discussed.


## 2. Power System Modeling Connecting ZnO

Transformer is connected to the power system while one of the three switches are open and only two phases of it are energized, which induced voltage in the open phase. This voltage back feeds the distribution line. Ferroresonance phenomenon is occurring if the distribution line is capacitive.Base system model is adopted from [1] while ZnO is added to the initial ferroresonance circuit. Linear approximation of the peak current of the magnetization reactance can be presented by the following:
(1)iL=aλ,
where λ is flux of the transformer coil and qis transformer nonlinear curve index. However, for very high currents λ-i characteristic of the transformer can be demonstrated by the polynomial in the following:
(2)iL=aλ+bλq.
The differential equation for the circuit in Figure 1 can be derived as follows. Polynomial coefficient of (2) is tabulated in Table 1:
(3)vl=dλdt,dvldt=dEdt-1RCdλdt-1C(aλ+bλq)-1C(vlk)α,dvldt=d2λdt2,
where vl=dλ/dt, ω represents the power frequency and E is the peak value of the voltage source. Also, α and k are ZnO parameters as shown in Figure 1.Table 1
Typical value ofq and its coefficients.
q
(a)
(b)
11
0.0028
0.0072
7
3.14
0.41Figure 1
Final reduced model of power system for ferroresonance studies.


## 3. Simulation Results

Case 1 (Power System Behavior Connecting ZnO).
In this section, effect of ZnO on mitigating ferroresonance phenomenon is investigated. Power system parameters values are tabulated in Table2.
It was shown that for proper representation of the saturation characteristics of a power transformer the values ofq are 5, 7, and 11 [1].Table 2
Power system parameters considered for simulation.
Parameters
Actual value
Per unit value
S
base
50 MVA,
—
V
base
635.1 kv
—
I
base
78.2 A
—
R
base
8121.4 Ω
—
C
4.9μF
0.07955 pu
R
core
4.5 MΩ
554.09 pu
ω
377 (rad/sec)
1 pu
K
—
2.5101 pu
α
—
25Behavior of the system is analyzed by consideringq degrees 7 and 11.Figures2 and 3 show time domain and phase plan simulation of the system behavior considering ZnO effect. By connecting ZnO chaotic ferroresonance is changed to fundamental resonance for q=7 and subharmonic behavior for q=11. ZnO clamps this overvoltage successfully. ZnO effect is shown in Figures 2 and 3 for q=7 and in Figures 4 and 5 for q=11.Figure 2
Time domain simulation forq=7 (ZnO effect).Figure 3
Phase plan diagram forq=7 (ZnO effect).Figure 4
Time domain simulation forq=11 (ZnO effect).Figure 5
Phase plan diagram forq=11 (ZnO effect).Changing subharmonic resonance and ferroresonance overvoltage in the purpose power system by connecting ZnO is clearly shows by using phase plan diagrams in Figures4 and 5. Effect of ZnO surge arrester is clearly obvious, ZnO limit the ferroresonance overvoltage.By increasing degree ofq, amplitude of overvoltage remains in 2 p.u. Some sudden changes in nonlinear systems parameters may cause to the chaotic behavior. The noun “chaos” is used to describe the time behavior of a system when the behavior is nonperiodic and it never exactly repeats. In fact, most of the systems that will be studying are completely deterministic. In general we need these three factors to determine the behavior of the proposed power system. The time equations which are derived from ferroresonance equivalent circuit are given in (4) and (5). The value of parameters describing the system which all parameters values are tabulated in Table 2, then the initial conditions are most important factor for initiating ferroresonance and given in (6).For studying the nonlinear dynamical systems as like power system, we need some tools such as time domain simulation, phase plan, and bifurcation diagrams. So, one of the most useful nonlinear dynamical tools is the bifurcation diagram [11]. To produce this diagram, we record the value of the peak voltage of the transformer as a function of control parameter being varied [11]. Here, control parameter is input voltage of the power system. In practice, a computer records the sampled values and then plotted as a function of the value of the control parameter. Two of these bifurcation diagrams are shown in Figures 6 and 7, where we have recorded sampled value of transformer voltage as a function of input voltage applied to the ferroresonance circuit with the frequency fixed.Figure 6
Bifurcation diagram withq=7 (ZnO effect).Figure 7
Bifurcation diagram withq=7 (ZnO effect).According to the bifurcation diagram, it shows the ferroresonance overvoltage has been begun after 3 p.u. Tendency to chaos exhibited by the system voltage increases whileq increases too.


## 4. Power System Modeling Considering Earth Resistance Effect

The main purpose of connecting resistance to the star point of a transformer is to limit ferroresonance. Low impedance is conventionally defined as impedance that limits the prospective ferroresonance to the no load of the transformer. The value of impedance required is easily calculated to a reasonable approximation by dividing the rated phase voltage by the rated phase current of the transformer. Neutral earthing resistance is achieved using resistors so as to limit the tendency for the ferroresonance to persist due to inductive energy storage. These resistors will limit heat when ferroresonance current flows and are usually only short term rated (typically 30 secs) so as to achieve an economic design.So the typical values for various system parameters which have been considered for simulation were kept the same by Case1, while earth resistance is added to the system and its value is given below and new system configuration is shown in Figure 8:
(4)Rnatural=40kΩ.Figure 8
Equivalent circuit of the power system connecting ZnO and earth resistance.The differential equation for the circuit in Figure8 is given by the following:
(5)dvldt=dEdt-1RCdλdt-1C(aλ+bλq)-RnRd2λdt2=-Rn⋅adλdt-Rnqbλq-1dλdt-α(1k)α(d2λdt2)α-1.
For the initial conditions, we have:
(6)λ(0)=0.0;vl=dλdt(0)=2,
where Rn is the neutral earth resistance.


## 5. Simulation Results

Case 2 (Power System Behavior Considering Neutral Earth Resistance Effect).
Figures9 and 10 show the corresponding time domain simulation which shows the effect of earth resistance. Also Figures 11 and 12 show the phase plan diagrams for corresponding system including earth resistance effect. It is shown that subharmonic regions mitigates by applying earth resistance. It is shown that tendency to chaos reduce in this case. In the big degree of q, system has been simulated when neutral resistance is considered. It shows the ferroresonance carried out and there are no abnormal phenomena and chaotic region mitigates by applying neutral resistance.
By increasing in the degree ofq, there is no more change in the system behavior. By comparing this case of bifurcation diagram with the previous case, it is clearly obvious that neutral resistance successfully clamps the ferroresonance overvoltage. This effect is shown in Figures 13 and 14.Figure 9
Time domain simulation forq=7 (neutral resistance effect).Figure 10
Time domain simulation forq=11 (neutral resistance effect).Figure 11
Phase plan diagram withq=7 (neutral resistance effect).Figure 12
Phase plan diagram withq=11 (neutral resistance effect).Figure 13
Bifurcation diagram withq=7 (neutral resistance effect).Figure 14
Bifurcation diagram withq=11 (neutral resistance effect).


## 6. Conclusions

The nonlinear behavior of ferroresonance falls into two main categories. In the first, the response is a distorted periodic waveform, containing the fundamental harmonics of the fundamental frequency. The second type is a chaotic response. In both cases response’s phase plan diagram contains fundamental and odd harmonic frequency components. In the chaotic response, there is also distributed subharmonics frequency. The simulation results confirm that system goes to chaos and bifurcation occurs in the proposed power system. The presence of the neutral resistance causes to clamp the ferroresonance overvoltage. The neutral resistance successfully, reduces the chaotic region for higher exponents ofq. Simulation of system consists of two cases, at first, system modeling of power transformer including ZnO arrester and second, system contains ZnO and neutral earth resistance. Finally we compare the result of these two cases.



---

*Source: 525494-2012-07-30.xml*