Elektrotehniški vestnik 83(5): 259-265, 2016 Original scientific paper Ferroresonance in 35 kV isolated networks: causes and mitigations Amir Tokic1, Mensur Kasumovic1, Damir Demirovic1, Irfan Turkovic2 Affiliation: 1University of Tuzla, Faculty of electrical engineering, Bosnia and Herzegovina, 2Faculty of electrical engineering Sarajevo,Bosnia and Herzegovina E-mail: amir.tokic@untz.ba Abstract. The ferroresonant states associated with inductive voltage transformers may be generated by elimination of single-phase short-circuits in isolated distribution networks. The paper presents a qualitative analysis of causes of ferroresonance initiation, and significant parameters impacting ferroresonance, like zero-sequence system capacitance and the number of connected voltage transformers. Using simulations, the causes of short-circuit occurrence for different system scenarios are investigated and damping devices for ferroresonance mitigations in a 35 kV isolated network are presented. Keywords: ferroresonance, isolated distribution networks, inductive voltage transformer Feroresonanca v 35 kV izoliranem omrežju: vzroki in zmanjševanje V izoliranih distribucijskih omrežjih lahko pojavi feroresonanca pri izločitvi enofaznega kratkega stika pri induktivnih napetostnih transformatorjih. V članku analiziramo vzroke za začetek feroresonance in bistvene parametre, ki vplivajo na ta pojav: zmogljivost sistema ničtega reda in število priključenih napetostnih transformatorjev. S simulacijami kratkega stika v omrežju smo analizirali vzroke za nastanek feroresonance pri različnih scenarijih. Predstavili smo tudi gradnike za zmanjšanje feroresonance v 35 kV izoliranem distribucijskem omrežju. 1 Introduction Ferroresonance is a nonlinear dynamic phenomenon which belongs to the group of low-frequency electromagnetic transients [1]-[3]. Basic preconditions for the ferroresonance occurrence are a sine-wave source, capacitance, nonlinear inductance and relatively small loss in the system. Ferroresonance appears in the following modes: a periodical, pseudo-periodical, subharmonic and chaotic mode [4]-[5]. Results of ferroresonance are permanent overvoltages and overcurrents in the system, which can affect the normal operation mode of the equipment, equipment malfunction, or its total damage [6]-[7]. The most common examples of ferroresonances are serial feroresonance due to compensation, ferroresonance with irregular operations in three-phase breakers, and ferroresonance with an inductive voltage transformer [1]-[2]. With a high voltage inductive transformer, single-phase ferroresonance can be initiated between a grading and bus capacitance together with a nonlinear inductance voltage transformer. In the middle voltage isolated distributed networks, initiation of ferroresonance occurs between the zero-sequence system capacitance and non-linear inductance of the voltage transformer. In a case of isolated distribution networks, initiation of ferroresonance is usually due to elimination of a single-phase short-circuit from the system. This particular type of ferroresonance is the topic of this research. A special attention is devoted to the analysis of significant parameters on this type of ferroresonance. The remainder of this paper is organized as follows. Section II presents conditions of initiation of ferroresonance and analysis of important parameters in the electrical circuit for the ferroresonance occurrence. Section III describes ferroresonance in three-phase isolated networks due to single-phase short-circuits. Section IV shows the simulated results for the different scenarios in the electrical distribution network. Section V concludes this paper. 2 Parallel Ferroresonance For a simplified description of the ferroresonance phenomena, in most literature, a serial ferroresonance circuit [1]-[2], [5] is analyzed. However, in isolated three-phase networks, ferroresonance can be classified into parallel ferroresonance due to the parallel connection zero-sequence system capacitance (line or cable to the ground capacitance) and nonlinear inductance of the inductive voltage transformer. In order to qualitatively analyze conditions for initiation of ferroresonance in the inductive voltage transformer, a simplified nonlinear circuit shown in Fig. 1 is analyzed. Received 12 July 2016 Accepted 16 September 2016 260 TOKIC, KASUMOVIC, DEMIROVIC, TURKOVIC Figure 1. Simplified parallel ferroresonant circuit. In Fig. 1, u is the input terminal voltage, C is the system capacitance of lines or cables, R is the resistance due to system losses, and L is the nonlinear inductance of the voltage transformer iron core. For the electrical system in Fig. 1, in the time domain, we have the following: is(t) = iR(t) + [¿c(0 + ¿l(0] = *a(0 + i(t) (1) Let us suppose that the input current is a sine wave with effective value h. After transformation into a phasor, and in the domain of effective values, we get the following system of nonlinear algebraic equations: I = | f (U) -uc -u | (2) (3) corresponds to relatively small currents (voltages) in the electric circuit. When the current rises above critical value Ic, the operating point moves in the position of point C where the values of currents (voltages) in the circuit are large. When the input currents decrease to the starting values, large currents (voltages) remain in the circuit. Therefore, the circuit is in ferroresonance, which, by the theory of settings, remains indefinitely. To make it simple, to initiate ferroresonance, a change is sufficient, such as.an increased input current of the minimal value: ^Anin Ac Is (4) Whenever the input current changes above the critical value of A/ > A/min, ferroresonance is initiated. 2.2 Impact of the system capacitance Let us now consider the case of increasing the capacitance in a system from C to Ci > C, as shown in Fig. 3. The change in the system capacitance is due to different places of elimination of the ground fault, as well as to a change in the number of connected lines and cables in a system. In relations (2) and (3), the capacitance and resistor are shown with linear lines Ic = rnC-U and U = R-I. We can suppose that the inductance curve is defined with nonlinear relation Il = f(U). 2.1 Normal ferroresonant states At any time, the state of an electrical system can be obtained in the intersection on non-linear curves (2) and (3) in the voltage-current (U-I) the coordinate system as shown in Fig. 2. Figure 3. Nonlinear curves of a parallel ferroresonant circuit: a case of an increased system capacitance. The system capacitance changes the mutual position of nonlinear curves (2) and (3). From Fig. 3 it is clearly seen that increasing the system capacitance leads to a negligible change of the operating point in a normal operating regime (point A1). Increasing the input current for a minimal value leads to ferroresonance: ^Anin! ^ Ci Is Figure 2. Nonlinear curves of a parallel ferroresonant circuit: basic case. In a normal state, when the input current of the circuit is /s, intersection of curves (2) and (3) are point A which According to Figs. 2 and 3, the following is applies: ^Anin! — ^Anin (5) (6) Therefore, it can be concluded that during an increase in the system capacitance, a bigger change in the value of FERRORESONANCE IN 35 KV ISOLATED NETWORKS: CAUSES AND MITIGATIONS 261 the input current to initiate ferroresonance in the circuit is needed. In other words, while the system capacitance increases, the probability of ferrorresonance initiation in a parallel non-linear ferroresonant circuit decreases. This statement is true for any system capacitance for ferroresonance initiation, i.e. curves 1l and Ic intersect. For a very small parameter C, for example, for which Ic is tangential on 1l, there is only one point of intersection of curves (2) and (3). So, in these scenarios, ferroresonance cannot be initiated at all. 2.3 Impact of the nonlinear inductance curve Now consider the case of changing ("lowering") the curve of a nonlinear inductance (Fig. 4). This scenario happens when two or more inductive voltage transformers are connected on the same bus-bar, which makes the inductance at each point on a new curve smaller than the original inductance, i.e. on the corresponding points of the nonlinear curves, the relations are Li < L (Fig. 4). Figure 4. Nonlinear curves of a parallel ferroresonant circuit: the case of changing (decreasing) the inductance nonlinear curve. Also in this case, with regard to a normal case, the mutual position of nonlinear curves (2) and (3) changes. As clearly seen from Fig. 4, by decreasing the inductance of the nonlinear curve, there is again, a small change in the location of the operating point in the normal operating regime (point A1). In this case, ferroresonance will be initiated by an increased input current for a minimal value of the current: ^mi„2 = Ic2 - 's (7) Based on Figs. 4 and 2, the following is true: 4/mi„2 < ^/mi„ (8) One can conclude that for decreasing the non-linear inductance, there is a smaller change in the input current value needed to initiate ferroresonance in the circuit. In other words, decreasing the inductance at any point of the non-linear curve ("lowering" the nonlinear magnetizing curve) increases the probability of initiation of ferroresonance in a parallel nonlinear ferroresonant circuit. It should be noted that a change in the resistance due to the system loss lead to an eventual damping or elimination of ferroresonance in the system. Namely, as shown in Fig. 2 due to a significant decrease in resistance R, curves (2) and (3) can have only one interception point in the normal operating regime. Therefore, a significant decrease in resistance R can assist in ferroresonance damping in the parallel circuit according to Fig. 1. This is achieved by connecting a resistor in an open delta on the secondary winding of an inductive voltage transformer. 3 Ferroresonance in A Three-phase Isolated Network Due to Short-Circuit Elimination In this part we analyze a concrete example of a ferroresonance occurrence in a real part of an isolated electrical distribution network (Fig. 5). Ferroresonance is initiated by eliminating a single-phase short-circuit from the system, thus establishing a ferroresonance circuit over the neutral point of the system. In the ferroresonant circuit, the ground to phase capacitance of the cable or lines participates together with the associated nonlinear inductance of the inductive voltage transformers. Depending on parameters of the electrical system, and as well as on the time of occurrence and elimination of a short-circuit in this three-phase system, ferroresonance can be initiated. Dominant parameters for ferroresonance initiation are the zero-sequence system capacitance and nonlinear curve of the voltage transformer. Other important factors are the residual magnetic flux and loss in the system [1]-[2]. Ferroresonance occurs upon elimination of a singlephase short-circuit from the system. Then, at each system phase, ferroresonant overvoltages appear. On the other side, the ferroresonance initiation is verified by the voltage of the neutral point of the system, which after elimination of the short-circuit turns back to the zero value. However, this voltage is held during the short-circuit (as expected) and after its elimination (which is an abnormal i.e. ferroresonant state). For this reason, it is practical to record the phase voltages, as well as the neutral point voltage, to determine the character of the electromagnetic transients. The parameters of the real part of the analyzed 35 kV distribution network are: Parameters of the system, (j = 1,2,3): 262 TOKIC, KASUMOVIC, DEMIROVIC, TURKOVIC • source voltage: e7-(t) = V2 ^sin [wt + (j - 1) • 2^], • voltage and frequency: E = 35 kV, f = 50 Hz, • network resistance: Rs = 0.055 Q, • network inductance: Ls = 0.5 mH, • zero sequence system capacitance: C = 150 nF. Parameters of the inductive voltage transformer: • rated power: S = 50 VA, • winding resistance: Rp = 13.7 kD, • leakage inductance: Lp = 23.5 H, • iron core loss equivalent resistance: Rm = 65.95 MD, The measured effective nonlinear curve of the current-voltage (I-U) of the inductive voltage transformer is converted to current-flux (/-9) magnetizing curve following the procedure described in paper [8]. By using the fitting procedure, the nonlinear magnetizing curve of the voltage transformer is described using the following polynomial function: i = aq + by5 + cq1 (9) Constants in Eq. (9) are: a = 1.91 • 10-5, b = 2.47 • 10-15, c = 4.07 • 10-27. The single-phase short circuit in phase 1 is modeled as an ideal breaker with switched on at time T0 = 40 ms and switched off at time T1 = 295 ms. 4 Ferroresonance Simulation in aThree-phase Isolated Network In this part of the paper, results of a ferroresonance simulation made with a model of a three-phase distributed network given in Fig. 5 are described. The model is developed in the software package SimPowerSystems [9]. The model of a nonlinear inductance voltage transformer is proposed according to Eq. (10). When simulating a power transformer and inductive voltage transformers a special attention must be paid to solvers and integration steps as well as to the relative tolerance for the simulation configuration parameters. Namely, the form of the state space of a dynamic system with transformers exhibits an extremely stiff differential equation system [10]-[11]. So for this type of dynamic system simulations, the A and L stable numerical methods must be used as they are not prone to numerical oscillations and numerical instabilities. For this reason, for the ferroresonance simulation in a three-phase, a BDF-based solver ode 23tb is used, to implement a combination of the trapezoidal and BDF2 numerical method. Three different scenarios of the electrical distribution network are analyzed: (a) basic case, (b) impact of the nonlinear magnetizing curve of the inductive voltage transformer and (c) impact of the system capacitance. In each scenario, the phase voltages of the system and voltage of neutral point of the network are simulated. The total simulation time is set at Ttot =2 s. The maximum value of the integration step is Atmax = 1 l^s, and the relative tolerance is e = 10-6. 4.1 Basic case of the isolated networks In the basic case, an inductive voltage transformer is connected per phase and the basic zero-sequence system capacitance is C = 150 nF. All other parameters are the same as above in this paper. The simulated phase and voltage of the neutral point are shown in Fig. 6. As seen from Fig. 6, there is no ferroresonance after elimination of a short-circuit from the system. Figure 5. Equivalent model of ferroresonance in a three-phase isolated network: the case of initiation and elimination of singlephase short circuits. FERRORESONANCE IN 35 KV ISOLATED NETWORKS: CAUSES AND MITIGATIONS 263 60 40 r-i 2U g tu n 60 O á >-20 -40 -60 neutral point u o 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 Time [s] Figure 6. Phase voltages and the voltage of the neutral point, one inductive voltage transformer, C = 150 nF. Figure 7. Phase voltages and voltage of the neutral point, two inductive voltage transformers, C = 150 nF. 264 TOKIC, KASUMOVIC, DEMIROVIC, TURKOVIC 4.2 Impact of nonlinear curve magnetizing of an inductive voltage transformer This case involves two inductive voltage transformers per phase; the basic zero-sequence system capacitance is C = 150 nF. Again, all other parameters are the same as above in this paper. The simulated phase voltage and voltage of the neutral point are shown in Fig 7. Fig. 7 clearly shows that there is ferroresonance after elimination of a short-circuit from the system. The enlarged part of the voltage phase 1 (time interval from 1.8 to 2 s) shows the dominant seventh subharmonic component. The main reason for the ferroresonance occurrence is lowering the equivalent nonlinear magnetizing curve of two parallel inductive voltage transformers. Therefore, in each point of the equivalent nonlinear curve of the two transformers, the iron core equivalent inductance of the voltage transformer is decreased. 4.3 Impact of the system capacitance This case involves two per phase connected inductive voltage transformers and an increased zero-sequence system capacitance at the point of a single-phase short-circuit of value C = 575 nF. The results of the simulated phase voltages and of voltage of the neutral point are shown in Fig. 8. From Fig. 8 it is shown that in this case there is no ferroresonance after elimination of a short-circuit from the system. The main reason for the ferroresonance disappearance is the increasing zero- sequence system capacitance which significantly decreases the probability of initiation of the ferroresonance state. 4.4 Ferroresonance mitigation There are three basic types of the ferroresonance damping devices, usually connected in an open delta auxiliary winding of the voltage transformers: (a) damping resistor whose resistance can be calculated based on a permisible thermal load of a voltage transformer secondary winding; (b) electronically controlled resistor with a lower resistance than that of the standard damping resistor. The device protects the equipment against the ferroresonant phenomena, after ferroresonance detection, by a proper and prompt damping action; (c) serial connection of a very low-value resistor and saturable reactor, which must be separately dimensioned for each voltage transformer. Under normal operating conditions (no fault), the damping device is disconnected from the system. 6o 40 ?20 A & 0 -M >"20 -40 -60 —i-1— —i-rr 6o 40 r—i au g V n M ° S p ¡>-20 -40 -60 0.2 0.4 0.6 0.8 1 1.2 Time [s] 1.4 1.6 1.8 2 0.8 1 1.2 Time [s] 1.8 2 —I-IT- - phase 3 0.2 0.4 0.6 0.8 1 1.2 Time [s] 1.4 1.6 1.8 2 0.8 1 1.2 Time [s] 1.8 2 Figure 8. Phase voltages and voltage of the neutral point, two inductive voltage transformers, C = 575 nF. FERRORESONANCE IN 35 KV ISOLATED NETWORKS: CAUSES AND MITIGATIONS 265 Fig. 9 shows a simulation result for the ferroresonance case presented in section 4.2 (two voltage transformers and system capacitance of C =150 nF) with a connected damping resistor Rd = 25 Q in the auxiliary winding of a voltage transformer. Figure 9. Neutral point voltage, two inductive voltage transformers, C = 150 nF, the damping resistor connected. Grounded-Neutral Voltage Transformers", IEEE Trans. on Power Delivery, 11 (3), pp. 1546-1553, Jul. 1996. [5] M. Pejic, A. Tokic, "Impact of the system parameters on the ferroresonant modes", Elektrotehniški vestnik, 80 (1-2), pp. 8-12, 2013. [6] D. A. Jacobson, D. R. Swatek and R. W. Mazur, "Mitigating Potential Transformer Ferroresonance in a 230 kV Converter Station", Computer Analysis of Electric Power System Transients: Selected Readings, IEEE Press, pp. 359-365, 1997. [7] W. Piasecki, M. Florkowski, M. Fulczyk, P. Mahonen, W. Nowak, "Mitigating Ferroresonance in Voltage Transformers in Ungrounded MV Networks", IEEE Trans. on Power Delivery, 22 (4), pp. 2362-2369, Oct. 2007. [8] A. Tokic, I. Uglešic, V. Milardic, G. Štumberger, "Conversion of RMS to instantaneous saturation curve: inrush current and ferroresonance cases", The 14th Int. IGTE Symposium on Numerical Field Calculation in El. Engineering, Graz, 2010. [9] "SimPowerSystems", User's Guide, Mathworks, 2013. [10] A. Tokic, J. Smajic, "Modeling and Simulations of Ferroresonance by Using BDF/NDF Numerical Methods", IEEE Trans. on Power Delivery, 30 (1), pp. 342-350, Jan. 2015. [11]A. Tokic, V. Milardic, I. Uglešic, A. Jukan, "Simulation of three-phase transformer inrush currents by using backward and numerical differentiation formulae", Electric Power Systems Research, 127, pp. 177-185, Oct. 2015. As seen ferroresonance is eliminated after switching on the damping resistor. 5 Conclusion The ferroresonance state in the inductive voltage transformers operating in neutral electrical distribution networks can be initiated as a result of clearing operations of single-phase short-circuits. In this paper, some significant parameters for initiation of ferroresonance after elimination of a short-circuit are researched. It is shown that increasing the zero-sequence system capacitance decreases the probability of initiation of ferroresonance. Also, increasing the number of the connected inductive voltage transformers increases the probability of the ferroresonance entering the neutral networks. Damping devices used to mitigate ferroresonance in neutral networks are presented. Our future work will include measurements of ferroresonance in three- phase neutral networks occuring to the elimination of short-circuits. References [1] A. Tokic, V. Milardic, "Power Quality", Printcom, Graficki inzenjering, Tuzla, 2015. [2] WG on Modelling and Analysis of System Transients Using Digital Programs. "Modelling and Analysis Guidelines for Slow Transients - Part III: The Study of Ferroresonance", IEEE Trans. on Power Delivery, 15 (1), pp. 255-265, Jan. 2000. [3] WG C-5, "Mathematical Models for Current, Voltage, and Coupling Capacitor Voltage Transformers", IEEE Trans. on Power Delivery, 15 (1), pp. 62-72, Jan. 2000. [4] N. Janssens, Th. Van Craenenbroeck, D. Van Dommelen F., Van De Meulebroeke, "Direct Calculation of the Stability Domains of Three-phase Ferroresonance in Isolated Neutral Networks with Amir Tokic received his M.Sc. and Ph.D. degrees in electrical engineering and computing from the University of Zagreb, Zagreb, Croatia, in 2001 and 2004, respectively. Currently, he is a Professor at the University of Tuzla, Tuzla, Bosnia and Herzegovina. His areas of interest include power system transients, power quality, and applied numerical and optimization methods. Mensur Kasumovic received his M.Sc. and Ph.D. degrees in electrical engineering from the University of Tuzla, Tuzla, Bosnia and Herzegovina, in 2006 and 2012, respectively. Currently he is an assistant Professor at the University of Tuzla, Tuzla, Bosnia and Herzegovina. His areas of interest include electric motor drives and power electronics. Damir Demirovic received his M.Sc. and Ph.D. degrees in electrical engineering and computing from the University of Tuzla, Tuzla, Bosnia and Herzegovina, in 2006 and 2011, respectively. Currently he is an assistant Professor at the University of Tuzla, Tuzla, Bosnia and Herzegovina. His areas of interest include modeling and simulations, pattern recognition and signal processing. Irfan Turkovic received his M.Sc. and Ph.D. degrees in electrical engineering from the University of Sarajevo, Bosnia and Herzegovina, in 2003 and 2010, respectively. Currently, he is an associate Professor at the University of Sarajevo -Faculty of Electrical Engineering. His topics of interest include cathodic protection, electrical measurements, safety in low voltage network and EMC in power systems.