Advanced
Transformer Winding Modeling based on Multi-Conductor Transmission Line Model for Partial Discharge Study
Transformer Winding Modeling based on Multi-Conductor Transmission Line Model for Partial Discharge Study
Journal of Electrical Engineering and Technology. 2014. Jan, 9(1): 154-161
Copyright © 2014, The Korean Institute of Electrical Engineers
  • Received : February 25, 2013
  • Accepted : July 19, 2013
  • Published : January 01, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Seyed Mohammad Hassan Hosseini
Corresponding Author: Dept. of Electrical Engineering, Islamic Azad University, South Tehran Branch, Iran. (Smhh110@azad.ac.ir)
Peyman Rezaei Baravati
Dept. of Electrical Engineering, Islamic Azad University, South Tehran Branch, Iran. (Peyman.rezaei@gmail.com)

Abstract
To study and locate partial discharge(PD) and analyze the transient state of power transformer, there is a need for a high frequency model of transformer winding and calculation of its parameters. Due to the high frequency nature of partial discharge phenomenon, there is a need for an accurate model for this frequency range. To attain this goal, a Multi-Conductor Transmission Line (MTL) model is used in this paper for modeling this transformer winding. In order that the MTL model can properly simulate the transformer behavior within a frequency range it is required that its parameters be accurately calculated. In this paper, all the basic parameters of this model are calculated by the use of Finite Element Method (FEM) for a 20kV winding of a distribution transformer. The comparison of the results obtained from this model with the obtained shape of the waves by the application of PD pulse to the winding in laboratory environment shows the validity and accuracy of this model.
Keywords
1. Introduction
With regard to the importance of transformers to be in circuit in power systems, their maintenance is very important. On the other hand, the insulation breakage of power transformers is of their major problems resulting from partial discharge gradually destroying the insulation. By early detection of partial discharge, a bulk of the repair costs is reduced [1] . The major methods for partial discharge locating can be divided into electric and non-electric methods. Most of the studies have been done on electric methods. Detection and locating of partial discharge in this method is based on obtaining the intended signals at the transformer terminals. By analyzing the recorded signals, it is possible to find out the exact place of partial discharge (PD) in the winding. For locating the partial discharge by electric method there is a need for a high frequency model of transformer winding. Many research activities have been performed on transformer winding modeling for the study and analysis of the propagation of PD signals for their location. With regard to these studies, two major methods which have been used can be proposed for the winding modeling: winding modeling by MTL model [2] , and detailed winding model [3] .
In the detailed model, called internal model or ladder RLC, the smallest space in a winding of a turn or a group of turns such as a disc or a pair of discs is modeled by an RLC circuit and the collection of the winding is in the form of a ladder network. Each element in the network has a capacitor relative to the ground and the inductances have also mutual induction. It is difficult in this method to estimate the circuit parameters with regard to the data related to the geometrical dimensions of transformer and the materials employed [4] .
Contrary to the above method, the distributed parameters of transmission line are used for winding model in the MTL modeling method. Due to the presence of inductances and capacitors in transient states, traveling waves are created which are propagated at different speeds inside the winding and based on the total effects of these transient waves, a state of resonance is created. Based on the theory of multi-conductor lines, the winding modeling is in the form of transmission lines which are parallel in positioning, but are in series electrically [5 , 6] . In this modeling, each winding unit is considered as one transmission line with distributed parameters and their collection forms a multi-conductor transmission system. This method is also applicable for electric machines, transformers, and cables [7] .
In this paper, the model of MTL is used for the simulation of transformer winding in the study of high frequency phenomenon effects. A new solution approach to obtain the MTL model parameters by the use of finite element method (FEM) is also presented.
This paper is organized as follows: MTL model is introduced in the second part and the manner of obtaining the electrical parameters is explained. In the third part, a new method is explained for the manner of obtaining the electrical parameters of the model through finite element method (FEM). Then, a comparison is made in the fourth part between the model simulation results and the laboratory results to assess and validate the model, and at the end the final conclusion is presented.
2. Multi-Conductor Transmission Line Model (MTL)
In the study of the partial discharge along the winding, the transformer winding is considered as a single input multiple output system (SIMO). Its input is the PD signal and its outputs are the measurement current signals of the terminals. Fig. 1 shows the transformer winding with n turns modeled by MTL method. Each turn of the transformer winding is modeled as one single-conductor transmission line in this model [8] .
PPT Slide
Lager Image
Multi-conductor transmission line model for transformer winding
In Fig. 1 when each turn of the winding is expressed by one transmission line, the propagation phenomenon in the transformer winding can be expressed through telegraphic equations in time domain:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where u(x,t) and i(x,t) are the voltage and current vectors at position x along the line; R,L,G and C are the per unit length resistance, inductance, conductance and capacitance matrix, respectively.
Due to the symmetry of electrical parameters, the telegraphic equations in frequency domain are as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where
PPT Slide
Lager Image
( x ),
PPT Slide
Lager Image
( x ) are phasor voltage and currentvectors; [P] 2 = [Z][Y], ([P t ] 2 ) = [Y][Z],[Z] = R + j. 2πf.[L] and [Y] = [G] +j. 2πf. [C] .[Z] and [Y] are impedance and admittance matrixes of the model, f is the frequency.
By the use of this model based on the theory of traveling waves, the calculation of the voltage and transient current is possible through all the points along the winding.
Voltage vectors and current at x distance from the beginning of the line can be obtained by Eqs. (3) and (4) as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where,
PPT Slide
Lager Image
i and
PPT Slide
Lager Image
r are the voltage vectors with regard to the boundary conditions, [ Y 0 ] = [ Y ][ P ] –1 of model specification admittance matrixand [ P ] Product of matrixes Z and Y .
By the use of boundary conditions which are as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
Where uh ( i ) and Ih ( i ) are sending end voltage and current of i th transmission line and ut ( i ) and It ( i ) are its receiving end voltage and current. If PD pulse is applied to K th turn of the winding, Eq. (7) for i = k -1 is corrected as follows:
PPT Slide
Lager Image
Where I PD is the PD pulse current. By the application of the boundary conditions, the currents of the sending end and receiving end can be expressed in terms of voltage vectors:
PPT Slide
Lager Image
By inversing matrix [Y], the transfer function matrix [T] of transformer will be as follows:
PPT Slide
Lager Image
The transformer bushing can be included in the equations by C B capacitor connected to the beginning of the line.
After solving the stated equations, it is possible to calculate the current equation produced due to the occurrence of PD at the head and at the end of the winding. By considering the above conditions, the transfer functions for the location of PD occurrence toward the winding phase ( TFL ) and up to the neutral end of the winding ( TFN ) is as follows [9] :
PPT Slide
Lager Image
PPT Slide
Lager Image
Where, in Eqs. (12) and (13) N = n +1 and n is the number of the conductors in the model.
- 2.1 Calculation of the model parameters
- A. Capacity matrix
Calculation of the capacitor capacity plays a significant role in the accuracy of this model. The capacitor capacity between each turn of the winding and the next turn is obtained as follows [10] :
PPT Slide
Lager Image
Where, ε 0 is vacuum dielectric constant and ε p is the relative dielectric coefficient of the paper and α p is the thickness of the paper at the two sides of the conductor while, d a is the average diameter of the winding and h is its height of conductor.
- B. Inductance matrix
Self-inductance for each turn of the winding and its mutual inductances form the inductance matrix with other turns of the winding. Magnetic flux penetration inside the iron core can be ignored in high frequencies assuming the winding with an air core. In mutual inductance calculation we have:
PPT Slide
Lager Image
Where, L n is the mutual matrix; C o is the capacity matrix in vacuum medium;
PPT Slide
Lager Image
is the propagation speed of electromagnetic waves in vacuum; ε 0 is vacuum dielectric constant; μ 0 is permeability of the vacuum and E is the unit matrix.
With regard to the skin effect, the self-inductance at high frequencies is calculated as follows:
PPT Slide
Lager Image
Where, in Eq. (16) Rs is the resistance due to the skin effect, and f is the frequency. The whole inductance is calculated by Eq. (17):
PPT Slide
Lager Image
- C. Resistance calculation
The resistance equation of the conductor length unit considering the skin effect can be represented as follows:
PPT Slide
Lager Image
Where, d 1 and d 2 are the width dimensions of a rectangular conductor; μ is the conductor magnetic permeability coefficient; σ is the electric conductivity capability; and f is the frequency.
3. Calculation of Parameters by Finite Elements Method (FEM)
The Finite Elements Method (FEM) is a numerical method used for approximate solution of partial differential equations and also for the solution of integrals. The finite elements method as an accurate and valid method has great applications in simulating the behavior of transformers in different conditions. Of these applications, it can be referred to the calculation of leaking fluxes and electromagnetic forces. By the use of FEM method, the circuit parameters equal to winding are calculable by considering the details of the transformer structure. In this paper, ANSOFT MAXWELL software, which solves problems on the basis of the analysis of the elements, is employed to obtain the winding parameters.
- 3.1 Calculation of capacity matrix
Based on electromagnetic laws, the total stored electric energy W in a space with a volume of V can be obtained by the following Eq. [11].
PPT Slide
Lager Image
Where, Ui and Uj are the voltage relative to conducting ground i and j; Cii the self-capacitor, and Cij is the mutual capacitor capacity, respectively. The energy method is a method by which the capacity capacitor is calculated. To calculate the coefficients in this method, Cii is applied to the desired constant voltage i th conductor and the remaining voltage of conductors is equal to zero and the energy is calculated in each of the three states. Then, the capacitor capacities are simultaneously obtained. In this state, as it is shown in Fig. 3 the total energy W is calculated as it was mentioned above. The diametric elements of the matrix C can be obtained by the following equation:
PPT Slide
Lager Image
PPT Slide
Lager Image
Energy around the 32th turn during calculation of capacitor capacity
Where, D i is the density of the generated electric flux in conductor i when a 1-volt voltage is applied to the conductor; and E j is the intensity of the generated electric field in conductor j when conductor j voltage is 1. To calculate the C ij coefficients, constant voltages are applied to conductors i and j and the remaining conductors are kept at zero voltage. Like the calculation of C ii , by calculating the stored electric energy in space, the non-diametric elements of C matrix are obtained. The major non-diameter elements of C matrix can be obtained as follows:
PPT Slide
Lager Image
- 3.2 Calculation of inductance matrix
The stored energy in the magnetic field which couples two conductors is obtained as follows:
PPT Slide
Lager Image
Where, W ij is the stored energy in the magnetic field of i and j couplers; I is the i conductor current; B i is the density of the magnetic field resulting from the passage of 1-ampere current through i conductor; and H j is the magnetic field resulting from the passage of 1-ampere current through j conductor. To calculate the inductance matrix at each stage, by applying the current to the conductors and considering other conductors as zero, energy (W) is calculated and diametric and non-diametric elements of L matrix are obtained. The calculated values of inductance for different turns of winding are shown in Figs. 6 to 8 .
PPT Slide
Lager Image
Inductance diagram in terms of turns for the first turn at frequency 2 and 5 MHz, L1j (j= 1~38)
PPT Slide
Lager Image
Inductance diagram in terms of turns for the 20th turn at frequency 2 and 5 MHz, L20j (j= l ~ 38)
PPT Slide
Lager Image
3-D view of inductance for 38 turns of winding
- 3.3 Calculation of resistance matrix
To calculate the resistance, we can first calculate ohms losses in a conductor. In fact, by the use of FEM method and solving the Maxwell equations, it is possible to obtain the field around every conductor and after calculating the current density by the use of equation J =∇ × H , the losses generated in the conductor can be calculated. The generated losses in the conductor are equal to:
PPT Slide
Lager Image
Where, is the electric conductive of the conductor and J is the current density. Then, by using the following equation we can obtain the resistance [12] :
PPT Slide
Lager Image
Simulator considers the peak of flowing current in each conductor equal to one ampere, so the resistance is equal to 2P. To calculate inductance and resistance in MHz range, the skin effect should be considered in calculations.
Fig. 9 shows the manner of resistance size change in terms of turns for the 20 th turn.
PPT Slide
Lager Image
Resistance diagram in terms of turns for the winding 20th turn at frequency 2 and 5 MHz, R20j (j= 1~38)
The specifications of the tested transformer winding depicted in Fig. 2 is presented in Table 1 .
PPT Slide
Lager Image
The tested distribution 20kV transformer winding
Specifications of tested winding
PPT Slide
Lager Image
Specifications of tested winding
- 3.4 Modeling results
Magnetic results obtained from 20 kV transformer winding modeling are shown in Figs 3 to 5 .
PPT Slide
Lager Image
2-D magnetic flux lines in the model
PPT Slide
Lager Image
A view of the sample transformer meshing (winding 2-D meshing) at 2 MHz frequency
In order to study the amount of accuracy of presented parameters, the measured and calculated results are compared with one another. The measured and calculated magnitude of transfer function of turn 20 compared to the beginning of the winding has been shown in Fig. 10 . The blue diagram shows the transfer function magnitude of the transformer winding obtained in laboratory, and the red diagram illustrates the results of simulation of transfer function obtained from parameters.
PPT Slide
Lager Image
The measured and calculated magnitude of transfer function of turn 20 compared to the beginning of the winding
Some of the calculated parameters for the winding in the experiment using analytical relations and finite elements method (FEM) have been presented in Table 2 .
The calculation and FEM results
PPT Slide
Lager Image
The calculation and FEM results
4. Study of the Model Validity
The distribution 20kV transformer winding was modeled and simulated by using the theory of MTL model and with the help of MATLAB software, and the parameters of this model were calculated by the Finite Elements Method in this paper. Having the quantities of transfer functions and pulse supposition for the partial discharge current using relations (12) and (13) we can obtain the current amounts in the beginning and in the end of the coil for various quantities of injection locations of partial discharge pulse where there are selected points for the injection of partial discharge pulse. To study the model, the partial discharge pulse similar to the data obtained from lab measurements was injected into the turns of 6 and 24 (K=6,24) of the simulated winding and the signal of the winding head which was a response to the PD pulse was measured. Fig. 11 shows the partial discharge pulse and Fig.12 shows the frequency spectrum of this signal.
PPT Slide
Lager Image
Diagram of partial discharge pulse (Rise time 0.2 μs , pulse width 8.5μs)
PPT Slide
Lager Image
Frequency spectrum of partial discharge calibrate signal
Also, in laboratory the PD pulse with the specifications of Fig. 11 is applied to 6 and 24 turns and is measured from the head of the winding. The partial discharge measurement circuit is as shown in Fig. 13 .
PPT Slide
Lager Image
Employing PD pulse to various points in the coil and measuring them at the two ends of the coil
These measured signals are in time domain according to Fig.14 , and are transferred to frequency domain to be compared with the simulation results.
PPT Slide
Lager Image
The measured signals at the first of the 20 kV winding resulting from the application of the PD pulse to turns 6 and 24
Fig. 15 shows the measured signals transferred to frequency domain.
PPT Slide
Lager Image
The measured signals transferred to frequency domain
Fig. 16 shows the output results obtained from simulation and laboratory. By comparing the obtained results in the figure, it is observed that the result of simulation is very close to the results obtained from laboratory. The small difference between the results shows the accuracy and good validity of this model.
PPT Slide
Lager Image
Results of simulation and measurement of the shape of the waves of the first of the winding resulting from the application of the PD pulse. A: turn 6; B: turn 24.
To perform more studies on the results, the correlation coefficient was applied by the use of formula (25) for all the turns resulting from simulation and measurement depicted in Table 3 .
PPT Slide
Lager Image
Correlation coefficient of 20kV transformer winding
PPT Slide
Lager Image
Correlation coefficient of 20kV transformer winding
The results of Table 3 show that the diametric coefficients (the correlation coefficients in the simulated turns with the measured turns of the same location) show the most dependence on each other. This shows that the simulation results enjoy good accuracy.
In general, these results show that the proposed model can properly study the partial discharge pulse propagation along the length of the winding.
5. Conclusion
The transformer winding was simulated by the use of multi-conductor transmission line in this paper. The efficiency of the multi-conductor transmission line model is greatly dependent on the MTL model parameters in the studies of transformer internal phenomena. With regard to inevitable approximations present in the calculation of the MTL model formula and also transformer generated tolerances along with the model intrinsic limitations, it is not practically possible to calculate the MTL model parameters with high accuracy in most cases. It was shown in this paper that it was possible to estimate the MTL parameters with high accuracy by the use of finite elements model (FEM) through An soft Maxwell software. The comparison of the results obtained through a test on a 20 kV winding and the results obtained from simulation of this winding is indicative of the accuracy and validity of this model.
BIO
Seyed Mohammad Hassan Hosseini was born in Tehran, in 1969. He received his M.Sc. and Ph.D. degrees in electrical power engineering in 2000 and 2005 from Azad University South-Tehran Branch and Science & Research Branch, Tehran, Iran, respectively. His research interest is transient modeling of transformer, Partial Discharge, High Voltage Engineering, Electrical Insulation and Substation.
Peyman Rezaei Baravati was born in Tehran, in 1985. He received his M.Sc. degree in Electrical power Engineering from the Azad University South-Tehran Branch in 2013, His research interests include partial discharge measurement and localization in transformers, high frequency modeling of transformer windings and high frequency measurements on power transformers.
References
Morched A.S. , Marti L. , Brierley R.H. , Lackey J.G 1996 “Analysis of internal winding stresses in EHV generator step-up transformer failures” IEEE Transactions on Power Delivery 11 (2) 888 - 894    DOI : 10.1109/61.489348
Lv Fangcheng , Liu Yunpeng , Liu Lei , Li Chengrong 2005 “Pulse Propagation Model of Partial Discharge in Transformer Winding” Electrical Insulating Materials. (ISEIM 2005). Proceedings of 2005 International Symposium 3 729 - 732
Wilcox D.J. , Hurley W.G. , Conlon M 1989 “Calculation of self and mutual impedances between sections of transformer windings” IEE Proceedings 136 (5) 308 - 314
Hosseini S.M.H. , Vakilian M. , Gharehpetian G.B 2008 “Comparison of Transformer Detailed Models for Fast and Very Fast Transient Studies” IEEE Transactions on Power Delivery 23 (2) 733 - 741    DOI : 10.1109/TPWRD.2008.915795
Hosseini S.M.H. , Vakilian M. , Gharehpetian G.B 2009 “Partial discharge location in transformers through application of MTL model” international conference on power systems transients (IPST2009) Kyoto, Japan
Hosseini S. M. H , Vakilian M , Ghaffarian M 2012 “Modeling of Partial Discharge in Transformer Winding by Multi-Conductor Transmission Line Model” Praise Worthy Prize, (I.R.E.E)
Guardado J.L. , Cornick K.J 1989 “A computer model for calculating steep-fronted surge distribution in machine windings” IEEE Transactions on Power Delivery 4 (1) 95 - 101    DOI : 10.1109/61.19193
Du Zhi-Ye , Zhao Chun , Ruan Jiang-Jun , Yu Shi-Feng , Yun Wen-Bing 2009 “Calculation of distribution parameters for Research on Propagation characteristic of PD in Transformer Winding” Power and Energy Engineering Conference 1 - 5
Hettiwatte S.N. , Crossley P.A. , Wang Z.D. , Darwin A. , Edwards G 2002 “Simulation of a transformer winding for partial discharge propagation studies” IEEE 2 1394 - 1399
Zhao Chun , Ruan Jiangjun , Du Zhiye , Liu Shoubao , Yu Yefeng , Zhang Yuanhang 2008 “Calculation of Parameters in Transformer Winding Based on The Model of Multi-Conductor Transmission Line” Electrical Machines and Systems 463 - 467
Li Yan , Du Jianping , Li Xiaohui , Li Dongxue 2011 “Calculation of capacitance and inductance parameters based on FEM in high-voltage transformer winding” Electrical Machines and Systems (ICEMS) 1 - 4
Venegas V. , Escarela R. , Mota R. , Melgoza E. , Guardado J.L 2003 “Calculation of electrical parameters for transient overvoltage studies on electrical machines” Electric Machines and Drives Conference 3 1978 - 1982