MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2017/2095896 2095896 Research Article A Method for Identification of Transformer Inrush Current Based on Box Dimension Gong Maofa 1 http://orcid.org/0000-0002-4677-6424 Zheng Ran 1 Hou Linyuan 1 2 Wei Jingyu 1 Wu Na 1 Gil-Lafuente Anna M. 1 College of Electrical Engineering and Automation Shandong University of Science and Technology Qingdao Shandong 266590 China sdust.edu.cn 2 State Grid Shandong Electric Power Maintenance Company Jinan Shandong 250021 China 2017 19112017 2017 10 07 2017 23 10 2017 19112017 2017 Copyright © 2017 Maofa Gong 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.

Magnetizing inrush current can lead to the maloperation of transformer differential protection. To overcome such an issue, a method is proposed to distinguish inrush current from inner fault current based on box dimension. According to the fundamental difference in waveform between the two, the algorithm can extract the three-phase current and calculate its box dimensions. If the box dimension value is smaller than the setting value, it is the inrush current; otherwise, it is inner fault current. Using PSACD and MATLAB, the simulation has been performed to prove the efficiency reliability of the presented algorithm in distinguishing inrush current and fault current.

Natural Science Foundation of Shandong Province ZR2016EEP10 Shandong University of Science and Technology 2015RCJJ073 Shandong Provincial Planning and Project, China J17KA074
1. Introduction

Longitudinal differential protection, as a part of main protection, can identify the inner fault and outer fault effectively in transformer. When an external fault of transformer is cleared or a no-load transformer is input, the high amplitude of inrush current will be generated , which can lead to maloperation of differential protection. In recent years, the accelerated construction of global energy interconnection, the increased capacity of transformers, and improved level of voltage have put forward higher requirements for transformer protection. Therefore, the accurate identification of inrush current is particularly critical .

Nowadays, there are three algorithms that are widely used in identification of inrush current: second harmonic discriminance, dead angle discriminance, and waveform symmetry discriminance. Due to the increase of the transformer capacity and the change of the core material, the second harmonic content in the inrush current will reduce. And long transmission line inductance and distributed capacitance may increase the second harmonic content. So, the change of second harmonic in power system reduces the accuracy of the second harmonic discriminance. The dead angle discriminance and the waveform symmetry discriminance have high requirements on the sampling frequency and are subject to the saturation of the current transformer. In view of shortages of above algorithms, many scholars have proposed improved algorithms or new algorithms. In literature , after the wavelet transform of signal, modular maximum can be extracted to identify inrush current. Although the algorithm has a high accuracy, the wavelet transform is of heavy computation and difficult when selecting basis function. Literature  proposed that Hilbert marginal spectrum can identify the inrush current, but there is a problem of oversampling. Literature  proposed grid fractal to obtain the changing curve of grid. The frequency characteristic of the changing curve can identify the inrush current, but the algorithm is short of fast motion. In literature , the proposed algorithm of combining wavelet transform with neural network has a high accuracy. However the neural network needs a lot of samples to train classifier. Literature  proposed that the morphological gradient of the sampled signal can be extracted by mathematical morphology. And its waveform area can be used to identify the inrush current.

Fractal is proposed by Mandelbrot in 1982. Fractal is widely used in complex and irregular graphics description . And it can process digital signal simply and quickly. Fractal also can quantify unordered behavior and does not require the same pace between signal frequency and sampling frequency and does not need time-frequency transformation. So, it is suitable for real-time signal processing . In this paper, the box dimension of fractal is applied to the identification of inrush current.

2. Theoretical Analysis of Inrush Current

In this paper, the single-phase transformer is used to analyze the producing mechanism of the inrush current. Its equivalent circuit is shown in Figure 1.

Transformer equivalent circuit after converting to primary side.

Denote winding voltage as u1(t)=Umsin(ωt+α) in steady-state operation. According to Faraday’s law of induction, as transformer is no-load and closing, using differential equation of u1t=dΦ/dt, the core flux Φ can be calculated as(1)Φ=u1tdt=-Φmcosωt+α+C,where ω=100π and Φm=Um/ω. C is integral constant, and its value is connected with voltage and magnetic flux in no-load closing time (t=0). Let the core remnant be Φr; then(2)C=Φmcosα+Φr.

When the no-load transformer is closed, the magnetic flux in the transformer core is(3)Φ=-Φmcosωt+α+Φmcosα+Φr,where Φ=-Φmcos(ωt+α) is steady-state flux. It is a periodic sinusoidal quantity. Φmcosα+Φr is a transient flux. Considering the transformer loss, it is a aperiodic component and decays with time. If Φmcosα is in phase with Φr, then the core flux reaches the maximum Φp=2Φmcosα+Φr at t=π. Obviously, the maximum flux generated during no-load closing may be much larger than the saturated flux Φs of the transformer, which makes the iron core saturation.

In Figure 2, when Φ<Φs, the core does not reach saturation, so the impedance of magnetizing branch is large. And the excitation current ie is far less than is. In this case, ie can be neglected. When Φ>Φs, the core enters into the saturated state. Then the impedance of excitation branch decreases. The excitation current increases sharply and forms magnetizing inrush current where amplitude ip may be up to 6–8 times the rated current.

Generating mechanism of inrush current.

Based on the above theoretical analysis, the magnitude of inrush current is related to the initial phase angle α of voltage when closing breaker, the voltage amplitude Um, the saturation flux Фs, and the size and direction of the remanence Фr in the core. In addition, inrush current of three-phase transformer is also related to its wiring methods .

3. Introduction of Box Dimension Algorithm

The basic property of fractal is the local self-similarity of the signal at different scales. The signal of inrush current is a nonstationary signal with singularity. And the singular signal has statistical self-similarity. This kind of self-similarity provides the basis for fractal analysis.

The singularity of the signal is usually described by the Lipschitz (Lip) exponent a. For the function f(t), and 0a1, given constant A and h0>0, when 0<h<h0, if(4)ft0+h-ft0Ahathen function f(t) has the Lip exponent a at t0. If a is a positive number, assume na<n+1, n is an integer, the function f(t) will be differentiable of n order at t0, and the nth derivative is point singularity at t0. If a=1, the function f(t) also will be differentiable at t0, and the derivative is bounded and discontinuous. If 0<a<1, the function f(t) will be singularity at t0. If a=0, the function f(t) is discontinuous and bounded at t0. If a<0, the function f(t) is singularity at t0.

Lip exponent can reflect the degree of signal similarity. The greater the value of a, the higher the degree of similarity about the self-similar process. Self-similar process has fractal features. Its fractal dimension D is associated with Lip index a. There are many algorithms to calculate the fractal dimension, and the calculation and experience estimation of box dimension are simpler than other algorithms.

This paper defines a simplified method for calculating the box dimension D of a one-dimensional discrete signal. Before calculating the box dimension in the plan set F, we can first make a δ coordinate grid intersecting set F on the plane, as shown in Figure 3 .

Definition of box number.

For a discrete digital signal, the size of the square grid δ is chosen as the reciprocal of sampling frequency f.

The formula of Nδ(F) is as follows:(5)NδF=1δi=1nxi-xi+1,where xi, xi+1 are the sampling points and n is the number of sampling points in a cycle and is even.

The grid size is increased to 2δ; the number N2δ(F) intersecting grids can be calculated as follows:(6)N2δF=12δi=1n/2maxx2i-1,x2i,x2i+1-minx2i-1,x2i,x2i+1.

The box dimension D is defined as follows:(7)D=lgNδF-lgN2δFlgδ-1-lg2δ-1.

Each sampled signal has a unique value about box dimension D. The values about different signals are generally not the same. Thus, the box dimension D can be used as an eigenvalue for signal identification.

4. Simulation Study on Inrush Current and Inner Fault Current

Through the PSCAD, the model of inrush current about three-phase transformer  was established as shown in Figure 4.

Simulation model of three-phase transformer inrush current.

The transformer of SFPZ7-120000/220 is selected as power transformer in simulation. Its wiring way is YNd11 and variable ratio is 220/35 kV. Modules (BRK1 and BRK2) of circuit breaker are used to control the breaking of circuit . And the different transformer parameters are set to simulate the inrush current in different circumstances. At the beginning of the simulation, the circuit breakers are disconnected. And the circuit breaker BRK1 is closed at 0.1 s. The current waveform of A phase is shown in Figure 5.

Waveform of inrush current.

From Figure 5, we can see that the inrush current waveform has serious distortion. Its amplitude is large and dead angle exists.

Through setting four fault modules on both sides of transformer, different types of inner fault and outer fault in Δ side and Y can be simulated. The inner faults of transformer differential protection (including interturn short circuit, ground fault, and phase short circuit) can be divided into the inner faults in normal operation and the inner faults in transformer no-load closing. Then the simulation results of the above two cases are analyzed to identify the difference of differential currents.

At 0.1 s, the Δ side in normal operation is set to have three-phase ground fault. Differential current of A phase is shown in Figure 6.

Waveform of inner fault current (fault occurs in normal operation).

A phase differential current, when the Δ side in no-load closing has the three-phase ground fault, is shown in Figure 7.

Waveform of transformer inner fault current (fault occurs in no-load closing).

The simulation model shown in Figure 4 cannot simulate interturn fault of three-phase transformer. Thus, we use three single-phase three-winding transformers to form a three-phase transformer. The secondary windings of two transformers are connected together as a new secondary winding . The fault point is set in the wiring of secondary winding. So changing the parameters of resistance and inductance in secondary winding can simulate turn-to-turn faults in different positions. The simulation model is shown in Figure 8.

Simulation model of three-phase transformer inner fault.

The three-phase differential current in the 10% of turn-to-turn fault is shown in Figure 9.

Waveform of turn-to-turn fault current.

According to the simulation results of the inner fault in Figures 6, 7, and 9, it is found that the amplitude of the differential current increases. But the waveform is basically maintained as sine wave and the amplitude of that does not decay after the fault occurs.

Compared with simulation results of Figures 5, 6, 7, and 9, the waveforms of inrush current and the inner fault current all have large amplitude. However there is a big difference in the waveforms shape. The difference in waveforms can be detected by their box dimension D which is a basis for identifying the inrush current .

5. Identification Criterion of Inrush Current Based on the Box Dimension

It is analyzed that the box dimension D can be used as the eigenvalue to distinguish the inrush current from the inner fault current. Box dimension D can be calculated according to (7). And it is compared with setting value DZD to determine whether the sampling signal is inrush current or inner fault current. The identification of the inrush current is(8)DDZD.

According to the various operations of the transformer, DZD=1.38 is determined ultimately.

D A , DB, DC are box dimension of A phase, B phase, and C phase, respectively. When at least two of them satisfy formula (8), it is determined that the current signal is inrush current, and then the differential protection device should be blocked. The algorithm flow of identifying the inrush current is shown in Figure 10.

The flow of inrush current identification.

6. Analysis of Simulation Results

The results of box dimension about fault currents are shown in Tables 1 and 2.

Box dimensions of inner fault current (faults occur in normal operation).

Fault type D A D B D C
Single-phase ground 1.4178 1.3780 1.4177
Two-phase short circuit 1.4174 1.4174 1.4174
Two-phase ground 1.3839 1.4174 1.3929
Three-phase short circuit 1.3804 1.4174 1.3863
Three-phase ground 1.3844 1.4174 1.3903
Turn-to-turn fault (10%) 1.4063 1.4176 1.4153
Turn-to-turn fault (20%) 1.3822 1.3955 1.4088
Turn-to-turn fault (50%) 1.3799 1.3895 1.4131

Box dimensions of inner fault current (faults occur in no-load closing).

Fault type D A D B D C
Single-phase ground 1.3895 1.3913 1.3603
Two-phase short circuit 1.4117 1.4116 1.4063
Two-phase ground 1.3812 1.4116 1.3768
Three-phase short circuit 1.3839 1.4116 1.3755
Three-phase ground 1.3840 1.4116 1.3755
Turn-to-turn fault (10%) 1.3918 1.4079 1.411
Turn-to-turn fault (20%) 1.3757 1.3948 1.4087
Turn-to-turn fault (50%) 1.3844 1.3895 1.4153

From Tables 1 and 2, it can be seen that at least two phase currents’ box dimensions do not satisfy formula (8), so the signal is inner fault current.

Since the inrush current will be affected by various factors, the parameters of the simulation model are changed to obtain inrush current and box dimension D in different cases. Three-phase box dimension value D, when the transformer remanence is 0 and the closing phase varies, is shown in Table 3.

Box dimensions of inrush current (remanence is 0).

No-load closing angle D A D B D C
1.3565 1.3026 1.3559
30° 1.3772 1.3464 1.3712
60° 1.3011 1.3563 1.3561
90° 1.3714 1.3740 1.3766
120° 1.3558 1.3564 1.3017
150° 1.3172 1.3771 1.3464
180° 1.3561 1.3039 1.3562

Three-phase box dimensions D, when the no-load closing angle is set to 0° and the transformer remanence varies, are shown in Table 4.

Box dimensions of inrush current (closing angle is 0°).

Remanence (T) D A D B D C
0.2 1.3224 1.3060 1.3229
0.4 1.3523 1.2938 1.3528
0.6 1.3418 1.3128 1.3447
0.8 1.3691 1.3351 1.3722

From Tables 3 and 4, it can be seen that the box dimensions DA, DB, and DC satisfy formula (8). Therefore, the criterion can accurately identify the inrush current. The differential protection device can be blocked correctly.

7. Conclusion

In this paper, after establishing the simulation model, the difference between the waveform of the inrush current and the fault current is analyzed. The proposed algorithm is used to calculate the box dimension of the different signals. Then the box dimension is compared with the setting value DZD. The algorithm has been tested by MATLAB; the results show that the algorithm is simple to implement and calculate and also has high accuracy.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors acknowledge the financial support from Shandong Provincial Natural Science Foundation, China (ZR2016EEP10), Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (2015RCJJ073), and Shandong Provincial Planning and Project, China (J17KA074).

Chandrasena W. McLaren P. G. Annakkage U. D. Jayasinghe R. P. An improved low-frequency transformer model for use in GIC studies IEEE Transactions on Power Delivery 2004 19 2 643 651 2-s2.0-1942504320 10.1109/TPWRD.2004.824429 Zhang L. Lv P. Shen H. Wang D. Zhou Z. Li Y. Liu Y. Wu C. Shen X. Yang G. Wang W. Analysis on protective relaying and its operation conditions of SGCC in 220 kV and above voltage AC systems in 2013 Power System Technology 2015 39 4 1153 1159 10.13335/j.1000-3673.pst.2015.04.042 2-s2.0-84928237018 Maofa G. Linyuan H. Danmiao S. Research of inrush current identification based on wavelet and modulus maximum Electrical Measurement & Instrument 2016 53 17 29 33 Maofa G. Wenhua X. Guoliang L. Identification of transformer inrush current based on hilbert spectrum Electrical Measurement & Instrument 2012 49 556 47 50 Ma J. Wang Z. Wang Z. Identify inrush current based on frequency characteristics of grille curve Electric Power Automation Equipment 2008 28 7 29 36 2-s2.0-49049119145 Maofa G. Meirong L. Fanjiao Y. Identification of transformer inrush current based on hilbert spectrum Electrical Measurement & Instrument 2015 52 6 124 128 Jiadong H. Weiqiang L. New algorithm to identify inrush current based on improved mathematical morphology Proceedings of the CSEE 2009 29 7 98 105 Quanmin X. Yuan L. Minshou Z. Blasting vibration signal analysis with wavelet and fractal portfolio analysis technique Journal of Vibration and Shock 2011 30 12 120 124 Jayaprakash A. Reddy G. R. Prasad N. Small target detection within sea clutter based on fractal analysis Procedia Technology 2016 24 988 995 10.1016/j.protcy.2016.05.217 Qun Y. Na C. Modeling and Simulation of MATLAB/Simulink Power System 2013 Beijing, China China Machine Press Bouda M. Caplan J. S. Saiers J. E. Box-counting dimension revisited: presenting an efficient method of minimizing quantization error and an assessment of the self-similarity of structural root systems Frontiers in Plant Science 2016 7 402, article 149 10.3389/fpls.2016.00149 2-s2.0-84960082649 Qi X. Yin X. Zhang Z. Sympathetic inrush current in a transformer and a method for its identification IEEJ Transactions on Electrical and Electronic Engineering 2016 11 4 442 450 2-s2.0-84973594889 10.1002/tee.22261 Wu Q. H. Zhang J. F. Zhang D. J. Ultra-high-speed directional protection of transmission lines using mathematical morphology IEEE Transactions on Power Delivery 2003 18 4 1127 1133 10.1109/tpwrd.2003.817513 2-s2.0-0141918387 Kim S. T. Jin B. G. Lee S. J. ν-i trend-based protective relaying algorithm for 3-phase power transformer Proceedings of the Power Engineering Society Summer Meeting 2001 IEEE 605 610 Kang Y.-C. Ok S.-H. Kang S.-H. A CT saturation detection algorithm IEEE Transactions on Power Delivery 2004 19 1 78 85 10.1109/TPWRD.2003.820200 2-s2.0-0347968154