MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2014/384790 384790 Research Article Vacuum Large Current Parallel Transfer Numerical Analysis http://orcid.org/0000-0003-1810-4644 Dong Enyuan 1 http://orcid.org/0000-0001-9551-7579 Hou Shengkai 1 Zheng Xiang 2 Qin Taotao 1 Liu Guixin 1 Zhao Bing 1 Shi Peng 1 College of Electrical Engineering Dalian University of Technology Dalian 116024 China dlut.edu.cn 2 School of Electronics and Information Engineering Dalian Jiaotong University Dalian 116028 China djtu.edu.cn 2014 382014 2014 17 04 2014 19 05 2014 19 05 2014 3 8 2014 2014 Copyright © 2014 Enyuan Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The stable operation and reliable breaking of large generator current are a difficult problem in power system. It can be solved successfully by the parallel interrupters and proper timing sequence with phase-control technology, in which the strategy of breaker’s control is decided by the time of both the first-opening phase and second-opening phase. The precise transfer current’s model can provide the proper timing sequence to break the generator circuit breaker. By analysis of the transfer current’s experiments and data, the real vacuum arc resistance and precise correctional model in the large transfer current’s process are obtained in this paper. The transfer time calculated by the correctional model of transfer current is very close to the actual transfer time. It can provide guidance for planning proper timing sequence and breaking the vacuum generator circuit breaker with the parallel interrupters.

1. Introduction

With the rapid development of the power system and the increase of the transmission capacity, it requires safer and more stable environment. The fault current’s breaking capacity and the longevity of high voltage circuit breaker that controls and protects power system are essential to the reliable operation of electric power systems. It can effectively reduce the average arcing time and peak arc current to use synchronous technology  to control the separate moment of breaker’s contacts through the zero point and it is useful to improve the breaking capacity of circuit breaker and reduce contact’s wear .

The rated current of single vacuum circuit breaker is less than 5 ka. And the single breaker is unable to burden the high rated continuous current and break the large short-circuit current . In order to improve the vacuum circuit breaker’s rated current and the capacity of the breaker, two parallel interrupters are used to share a large current especially for the generator circuit breaker and the occasion that needs to break the large current [7, 8]. However two pairs of contacts are impossible to absolutely separate at the same time due to the manufacture of breakers, the dispersion of actuators, and the arrangement which leads to a failure to break the whole fault current .

The parallel vacuum interrupters are used to share the rated large current. When the fault current occurs, main vacuum interrupter (MVI) is opened firstly. The fault current gradually is transferred from MVI to AVI (auxiliary vacuum interrupter). At the same time, the current constantly decays.

At the end of the current-transfer process, the fault current through AVI is removed with phase-control technology . But in the transfer process, the multiplication of AVI’s inductance and the instantaneous current is greater than MVI’s arc voltage, which leads to uncontrollable transfer process and the current cannot transfer from MVI to AVI. The successful breaking of the AVI depends much on the precise calculated transfer time. So an accurate model of the transfer current is needed.

In the paper , the mathematical model of transfer current is provided through the relevant theoretical derivation. It is helpful to analyze the transfer process. The arc resistance’s value in the model is constant, which leads to a large deviation for transfer time. So this model cannot provide accurate timing sequence for using the phase-control technology . The success rate of transferring is related to transfer current’s value, the time of breaking MVI, and so forth. This transfer current’s model needs to be corrected. In this paper, the vacuum arc resistance’s formula is obtained by Matlab through the analysis of experimental data of transfer current. The arc resistance changes with time and transfer current. Finally, the correctional model of transfer current is derived with the real arc resistance.

2. The Model of Transfer Current

At the beginning, the large rated current is shared by MVI and AVI as shown in Figure 1.

MVI and AVI share the large current.

When the short-circuit fault occurs, MVI is opened firstly. The large current gradually is transferred from MVI to AVI. The equivalent circuit is shown in Figure 2.

Equivalent circuit when MVI’s contacts are opened.

Short-circuit current of the system : (1)Is(t)=2Iesin(ωt+α-φd)+[Iesin(α-φ)-2Iesin(α-φd)]e-t/τsys.

In (1), Ie is the effective value before short circuit occurring. τsys is the decay constant of the system. ω is the current’s angular frequency. α is the initial phase angle of power source’s voltage when short circuit occurs instantly. φd is the impedance angle after short circuit occurring. And φ is the impedance angle before short circuit occurring.

The impedance of power source can be ignored when it is compared with the interrupter’s resistance and the system capacity is much larger. Considering the constant arc voltage Ua when the MVI is opened, the transfer current in Figure 2 can be expressed with  (2)i1+i2=Is(t),L1di1dt+R1i1+Ua=L2di2dt+R2i2ddddddddddddddd(T1tT2).

T 1 is the moment of breaking MVI. T2 is the moment of breaking AVI. R1 does not contain contact resistance between a pair of contacts when MVI is separated in (2) and it mainly consists of arc resistance and the resistance of internal conductive rod in MVI. So the differential equation (3) is derived from (2): (3)di1dt+(R1+R2)(L1+L2)i1=L2(L1+L2)dIs(t)dt+R2(L1+L2)Is(t)-Ua(L1+L2).

i 1 is the instantaneous value of the transfer current  in the following: (4)i1=Ce-((R1+R2)/(L1+L2))t-UaR1+R2+2IeL2-R2τsys(L1+L2)-(R1+R2)τsyse-t/τsys-2IeωL1R2-ωL2R1(R1+R2)2+ω2(L1+L2)2sin(ωt)-2Ie(R1+R2)R2+ω2L2(L1+L2)(R1+R2)2+ω2(L1+L2)2cos(ωt).

In (4), φd is π/2. R1 is the equivalent resistance of MVI’s branch. L1 is the equivalent inductance of MVI’s branch. R2 is the equivalent resistance of AVI’s branch. L2 is the equivalent inductance of AVI’s branch. Tk is the moment at which MVI is opened. The differential constant C is shown in the following: (5)C=e((R1+R2)/(L1+L2))Tk×{UaR1+R2-2IeL2-R2τsysL1+L2-(R1+R2)τsyse-Tk/τsysddddd+2IeωL1R2-ωL2R1(R1+R2)2+ω2(L1+L2)2sin(ωTk)ddddd+2Ie(R1+R2)R2+ω2L2(L1+L2)(R1+R2)2+ω2(L1+L2)2cos(ωTk)ddddd+2Iesin(ωTk+α-π2)ddddd+[Iesin(α-φd)-2Iesin(α-π2)]e-Tk/τsysUaR1+R2}.

The emulational value is deviated from the actual value, because arc resistance is assumed to be constant and unchanged with time in the model of transfer current.

3. Mathematical Model of Vacuum Arc Resistance

In the experiment, two parallel breakers are used to share the large current as shown in Figure 3. The interrupter’s model is BD-12/1250-20A produced by the Kaisaier. The large current is generated by the Synthesis loop. Its half cycle is 12.5 ms and the current’s value is adjustable.

The two parallel breakers.

The two breakers controlled by the controller are opened according to the timing sequence. The current signal of the two breakers is collected by Rogowski coil and the arc voltage is collected by the resistance-type voltage divider that is parallel with the MVI. In Figures 4 and 5, the current value of 1200 A corresponds to 1 V, which is the output of Rogowski coil. The curve 2 which is collected by resistance-type voltage divider amplified 22.6 times is the voltage between MVI. In this paper, the times of breaking MVI are 1.4 ms, 4.4 ms, 5.0 ms, 6.2 ms, 7.0 ms, and 8.3 ms. The transfer currents are 3.7 ka, 5.2 ka, 6.4 ka, 7.7 ka, and 9.6 ka. All current values are the peak of transfer currents.

Transfer current of 5.2 ka and voltage.

Transfer current of 9.6 ka and voltage.

The current that transfers from MVI to AVI is prevented by AVI’s reactance. So the moment of breaking MVI needs to avoid the current’s peak and it is better to choose the time in which current’s value is relatively low. In the moment that is close to current’s zero crossing, the success rate of transferring is higher. The enough computation-time of the controller is required to break AVI with the phase-control technology and the result in the moment is more obviously and more representatively compared with all experimental data. So the time of breaking MVI is selected at 8.3 ms in this paper. The waveforms of the oscilloscope are only listed in 5.2 ka and 9.6 ka due to limited space. MVI’s current (curve 1) and the voltage (curve 2) across the interrupter are shown in Figures 4 and 5.

In Figures 4 and 5, curve 2 is the voltage between MVI. At the beginning, the voltage is zero. When large current flows through MVI, its voltage has a slight uplift and maintains a stable value. At the moment at which MVI is opened, MVI’s current drops rapidly and the arc is generated and the voltage pulse  increases immediately. The arc voltage is higher with the increase of transfer current. But the voltage’s growth rate is low with the increase of transfer current. At the end of the transfer process, the arc is extinguished. The sinusoidal voltage is not the arc voltage but the voltage of AVI’s branch in Figures 4 and 5. The voltage declines until AVI is opened. At the same time, the voltage pulse is the superposition of AVI’s voltage and AVI’s arc voltage in the negative half-cycle. Finally, the voltage is raised to the voltage of current source’s capacitor and remains unchanged as shown in Figures 4 and 5.

It is shown how the MVI’s and AVI’s currents change in the first half-cycle when MVI is opened in Figures 6 and 7. MVI’s current is 1.8 times larger than AVI’s current before 8.3 ms. The large current of the MVI is transferred from MVI to AVI at 8.3 ms. In the transfer process, the MVI’s arc voltage in 9.6 ka is shown in Figure 8. Figure 8 is the amplification of arc voltage in Figure 5.

MVI’s and AVI’s currents in transfer current of 5.2 ka.

MVI’s and AVI’s currents in transfer current of 9.6 ka.

The arc voltage in the transfer process.

The real vacuum arc resistance is calculated with transfer current and arc voltage as shown in Figures 9 and 10.

The vacuum arc resistance in transfer current of 9.6 ka.

The vacuum arc resistance in all transfer currents.

The vacuum arc resistance is close to exponential growth. The arcing time becomes longer with the increase of transfer current. At the beginning of the transfer process, arc resistance grows very slowly at a low level. On the contrary, the arc resistance increases rapidly to a high value when the arc is extinguished.

The mathematical formula of the arc resistance is expressed with (6) through the analysis of experimental data: (6)Rarc(t)=aebt+600μΩ,a=2.794×106×e[-3.024×10-3(2.892×10-4I2MVI-0.8731IMVI+2729)]b=2.124×IMVI+2376.

In (6), t is the arcing time. IMVI is the transfer current.

The arc resistance is related with transfer current and arcing time. The arcing time becomes longer with the increase of the transfer current and the accumulation of heat becomes more and more. In (6), a becomes smaller and b becomes greater with the increase of the transfer current. In fact, arc resistance is related to contact material, the structure, the spacing of the contact, and so forth. The constants in (6) need to adjust in quite different experimental condition.

The red points are the experimental data and the blue curve is the calculated value by the mathematical model of vacuum arc resistance. The calculated arc resistance is very close to the actual value with the model as shown in Figures 11, 12, 13, and 14.

Arc resistance fitting in transfer current of 3.7 ka.

Arc resistance fitting in transfer current of 6.4 ka.

Arc resistance fitting in transfer current of 7.7 ka.

Arc resistance fitting in transfer current of 9.6 ka.

4. The Correctional Model of Transfer Current

The I(T1) is the part of particular solution of differential equation (3). At this moment, the arc is not generated and the R1  remains constant. The constant R1 in other items of (4) is replaced with R1+Rarc(t) and the correctional model of transfer current is deduced as shown in the following: (7)I(T1)=Ce-((R1+R2)/(L1+L2))t-UaR1+R2(t=T1),i1(t)=I(T1)+2IeL2-R2τsys(L1+L2)-(R1+Rarc(t)+R2)τsyse-t/τsys-2IeωL1R2-ωL2(R1+Rarc(t))(R1+Rarc(t)+R2)2+ω2(L1+L2)2sin(ωt)-2Ie(R1+Rarc(t)+R2)R2+ω2L2(L1+L2)(R1+Rarc(t)+R2)2+ω2(L1+L2)2×cos(ωt)T1<tT2.

The correctional transfer time is calculated with the correctional model of transfer current. The primary transfer time is calculated with the model of transfer current ((4)-(5)) in Figures 15 and 16. The actual transfer time is the experimental data.

Three kinds of transfer time’s comparison in different currents.

Three kinds of transfer time’s comparison at different time of breaking MVI.

The comparison of the three kinds of transfer time in different current is shown in Figure 15 and the trend of transfer time at different times of breaking MVI is shown in Figure 16.

The correctional model of transfer current by analyzing experimental data with Matlab can accurately reflect the changing of the real vacuum arc resistance. The correctional model is more accurate and the correctional transfer time is closer to the actual transfer time. It can provide a proper timing sequence for breaking AVI with the phase-control technology.

So the correctional model can provide more accurate transfer time than the primary model for using the phase-control technology. The main contribution of the models can provide guidance for planning proper timing sequence and breaking a vacuum generator circuit breaker with the parallel interrupters.

Because the vacuum arc is affected by many factors, there is much work to research the more accurate vacuum arc’s model.

5. Conclusion

The mathematical model of vacuum arc resistance and the correctional model of transfer current are established in this paper. By the analysis of experimental data and the results of simulation, the changing of real arc resistance is described and the real vacuum arc resistance is close to exponential growth. At the beginning of the transfer process, the arc resistance’s growth rate is very low. But when the arc is extinguished, the real arc resistance increases rapidly to a high value. The arcing time becomes longer and the arc voltage is higher with the increase of the transfer current. At the same time, the duration of the low arc resistance’s state maintains more time.

Although the mathematical model of vacuum arc resistance and the correctional model of transfer current are obtained under specific experimental condition, each of the important parameters in the models is adjustable according to specific circumstances. The real vacuum arc resistance’s model and the correctional model of transfer current are widely applicable.

Because the transfer time calculated by the correctional model of transfer current is very close to the actual transfer time and the deviation is smaller, the models can provide proper guidance for breaking a large generator current with the phase-control technology.

As future work, there is much work to integrate the robust control algorithm  into the breakers and regulate the speed of breakers’ movement. The authors plan to improve the scalability and fault tolerance  of the models at a satisfactory quality level and add the reciprocal determination process  where the models can take into account the speed of breakers and the type of fault current.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (no. 51277019).

CIGRE Working Group 13.07 Controlled Switching of HV AC Circuit Breaker Guidance For Further Application Including UnLoaded transFormer Switching, Load and Fault Interruption and Circuit Breaker Uprating 1999 Chengdu, China CIGRE Poltl A. Frohlich K. A new algorithm enabling controlled short circuit interruption IEEE Transactions on Power Delivery 2003 18 3 802 808 2-s2.0-0037480052 10.1109/TPWRD.2003.813635 Chen X.-S. Zou J. Generator protection vacuum circuit breaker High Voltage Apparatus 1998 36 1 34 37 Homma M. Sakaki M. Kaneko E. Yanabu S. History of vacuum circuit breaker and recent development in Japan Proceedings of the 21th International Symposium on Discharges and Electrical Insulation in Vacuum (ISDEIV '04) October 2004 Yalta, Ukraine 378 383 2-s2.0-16244366455 10.1109/DEIV.2004.1422626 Zhou P.-H. Dai M. Lou Y. Xiu M.-H. DC component time constant and zero off set phenomena of breaking current for 1000 kV circuit breaker High Voltage Engineering 2009 35 4 722 730 2-s2.0-65749098389 Cao L.-J. Li J. Exploration on high capacity generate or outlet protection vacuum circuit breaker Electrical Device 2004 5 9 84 85 Chen X.-S. Zou J.-Y. - He J.-J. Parallel high-current vacuum arc breaking High Voltage Apparatus 1997 4 29 32 Matsukawa M. Miura Y. M. Terakado T. Development of a vacuum switch carrying a continuous current of 36 kA DC 2 Proceedings of the 19th International Symposium on Discharges and Electrical Insulation in Vacuum 2000 415 418 10.1109/DEIV.2000.879015 Chen X.-S. Wang Z.-W. Zou J.-Y. Experiment al research on paralleled breaking process of vacuum circuit breaker High Voltage Apparatus 1998 33 5 33 37 Duan X.-Y. Liao M.-F. Ding F.-H. Zou J.-Y. Application and key technology analysis of controlled switching High Voltage Apparatus 2007 43 2 113 117 2-s2.0-34248216365 Duan X. Liao M. Ding F. Liu B. Zou J. Controlled vacuum circuit breaker for shunt capacitor banks Transactions of China Electrotechnical Society 2007 22 10 78 84 2-s2.0-39749157004 Ding F.-H. Zou J.-Y. Duan X.-Y. Research on optimum switching phase of phasing switches High Voltage Apparatus 2005 41 6 408 411 2-s2.0-30044437635 Zheng X. Zou J.-Y. Cheng Z. Current transfer process of paralleled vacuum generator circuit breaker High Voltage Engineering 2009 35 11 2709 2715 2-s2.0-77950578262 Liu Z.-Y. Wang J.-M. Wang Z.-Y. Yuan S. Current transfer process in open operation between parallel breaks of a vacuum generator circuit breaker Proceedings of the Chinese Society of Electrical Engineering 2000 20 12 20 23 2-s2.0-0034592247 Veselly V. Rosinova D. Robust MPC controller design for switched systems using multi- parameter dependent lyapunov function International Journal of Innovative Computing, Information and Control 2014 10 1 269 280 Dobre C. A cluster-enhanced fault tolerant peer-to-peer systems International Journal of Innovative Computing, Information and Control 2014 10 2 417 436 Harada K. Ishida Y. Complex spatial-temporal patterns developed from an interaction of functions and data International Journal of Innovative Computing, Information and Control 2014 10 1 405 416