A New Method for Detection and Evaluation of Winding Mechanical

0 downloads 0 Views 362KB Size Report
Jun 1, 2011 - fault diagnosis in actual power transformer windings. Index Terms—Transformer ..... short circuits and their spatial arrangement on the TF. For.
[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering

Volume 11, Number 2, 2011

A New Method for Detection and Evaluation of Winding Mechanical Faults in Transformer through Transfer Function Measurements Mehdi BIGDELI1*, Mehdi VAKILIAN2, and Ebrahim RAHIMPOUR3 1 Department of Electrical Engineering, Science and Research Branch, Islamic Azad University, Tehran, Iran 2 Department of Electrical Engineering, Sharif University of Technology, Tehran, Iran 3 ABB AG, Power Products Division, Transformers, Bad Honnef, Germany * Email: [email protected] 1 Abstract—Transfer function (TF) is an acknowledged method for power transformer mechanical faults detection. However the past published works mostly discovered how to specify the faults levels and paid less attention to detection of the type of faults using comparison of TFs. whereas, it seems important for most of the applications to specify the type of fault without opening the unit. This paper presents a new method based on vector fitting (VF) to compare the TFs and specify the type, level and location of the fault. For development of the method, and its verification the required measurements are carried out on four model transformers; under intact condition, and under different fault conditions (axial displacement, radial deformation, disc space variation and short circuit of winding) and the TFs are determined. Employing VF, the coefficients of TFs are determined with the required accuracy. Using those coefficients, a new index is introduced to specify the type, level and location of the fault in the winding. Convincingly good results were obtained. Therefore it is believed that this finding could be helpful in fault diagnosis in actual power transformer windings.

Index Terms—Transformer, Fault Diagnosis, Measurement, Transfer Function, Vector Fitting

I. INTRODUCTION Power transformers are among the most important equipments in the electrical power transmission and distribution systems and occurrence of any fault in these transformers will reduces the power system reliability, in other words increases the possibility of power supply interruption. In addition, it will impose extra costs required for their maintenance and for their transportation to a repair factory. Due to the existence of a strong competition in the electrical power supply industry, the importance of application of the transformer monitoring systems has increased over the time [1]. Different fault modes in transformers are classified into two main groups; the mechanical types, and the electrical types. Different methods have been offered for detection of each of these faults, such as: a – Oil analysis (Dissolved Gas Analysis, Furfural) b – Partial discharge measurements c – Dielectric response analysis d – Transfer function analysis. Among these methods, TF method is a well-known

method for detection of the winding mechanical faults [2][4]. TF method is a comparative method; in this method a new measurement is put besides a referential measurement. If the deviations were remarkable, the direction and magnitude of deviations should be studied and analyzed. Mechanical faults such as; axial displacement (AD), radial deformation (RD), disc space variation (DSV) and short circuit (SC) of the winding affect the TFs and cause dislocation of resonance frequencies and moreover will decrease or increase the magnitude of resonance frequencies in the TF. These changes will cause a relevant shift in the poles and zeros of TF. Therefore, by evaluation and comparison of the TFs, any variation in the structure of the winding due to a mechanical damage can be detected. TF evaluation and its comparison can be done in time domain or in frequency domain. The identification process is normally carried out in frequency domain through inspection of the rational functions. Some studies [5]-[6] have shown that VF method is an effective and accurate method for estimating the transformer TFs. In [7]-[8] have proposed an index (based on TF estimation using VF) for determination of AD and RD extent. That index played its role only to specify the level of AD or RD and does not specify the type of fault. In addition, their proposed index only perceives the variations in the poles of the rational function where as, during incidence of winding mechanical faults both the poles and the zeros change. Therefore continuation of the researches on VF application for fault detection is necessary. In this paper employing VF, the TFs of the transformer are estimated both during faulty condition and during normal condition. Introducing a new algorithm the TF is identified and compared with the reference one. Not only this algorithm is able to recognize the type of fault, but also it can specify the level and location of fault. Comparing this new method with one of the past methods, it is shown that the proposed index is more efficient, even in determination of the fault level. II. LITERATURE REVIEW AND PROBLEM DEFINITION TF analysis is a reliable method for winding failures detection in power transformers. Different methods are introduced for evaluation of the TFs. Some of these methods are verified against the practical samples of mechanical

23

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering faults [9]. Also, some analytical investigations have been performed using the equivalent circuit models [10]. On the other hand, transformer modeling using its internal geometry and the material properties provided the possibility of evaluating the sensitivity of the TF parameters against the type of faults [11]-[12]. In this paper, the important winding faults, which are most likely to be detected using the TF analysis, have been studied. These faults can be classified as follows: a- Axial Displacement [11]-[14], b- Radial Deformation [8], [12] c- Disc Space Variation [15] and d- Short circuit [14], [16] It should be noted that SC is not mechanical faults per se, but usually it do stem from them, because it can be detected using TF method. The above faults have been studied analytically and experimentally in [7]-[17]. These researches have studied the sensitivity of TFs parameters against those faults separately and have not discussed that how a specific mechanical fault can be found and can be classified as to be in one of the AD, RD, DSV or SC categories. In this analysis, the TF method is analyzed by a comparative method. However, in general the comparative algorithms introduced in the literature, can be classified into four major categories as follows. A. Algorithms Based on Artificial Intelligence Methods The methods such as neural networks, genetic algorithm, fuzzy logic, and so on can be included in this group [18][20]. Conclusions in these works emphasize the necessity for further investigations to classify the type and location of the fault. Moreover a limitation exists in these methods which made them applicable only for transformers like the one used as test transformer in the development of method. B. Algorithms Based on Mathematical and Statistical Methods Some of the methods of this category have been used for fault diagnosis in [21]-[24]. To establish the necessary measures for detection of the faults some parameters have been defined, namely: correlation coefficient (CC), standard deviation (SD). Disadvantages of these methods discussed in [21], [25]. Meanwhile all of the algorithms that employ specific equations with single-valued solutions for the evaluation of the measured TFs belong to this category. The investigations carried out in [26]-[30] introduce some additional parameters to evaluate the fault level. However, those works do not determine the type and location of the fault. Since in this class of algorithms the phase response is not evaluated, therefore some of them contain unreliable comparison parameter [25]. C. Algorithms Based on Electric Circuit Models Some investigations have been done to get information about the faults from the measured TFs by describing the TFs using an equivalent electric circuit model [31]. The parameters of the model representing the TFs can be calculated by the means of the details of the transformer design [12] or can be estimated using terminal measurements [32]. Determining the parameters of such models, which 24

Volume 11, Number 2, 2011 usually includes errors [12] and is more difficult for real windings and transformers, is ruled out by avoiding electric circuit models. D. Algorithms Based on Estimation Methods TFs of transformers can be approximated by rational functions consisting of two polynomials with real coefficients [5]. The estimation of such a rational function approximation to describe the measured terminal characteristic is employed for comparing TFs [7]-[8]. However, in this works only the AD and RD are considered; thus, specifying the type of fault by comparing different types of faults with each other is left as a subject to be studied. VF is a widely applicable and general methodology for estimation of TF parameters from frequency response measurement results [5]-[6]. So, in this work VF is used for estimation of TFs. The algorithms based on VF methods to analyze the process of fault detection are preferred in this work, because: a- Applying these algorithms is very simple and the results have an excellent match, b- The estimated functions are comparable in an easier way, c- As it is discussed in [25], to have an efficient diagnostic tool, the magnitude and the phase responses must be taken into account for the analysis, and this is possible with VF method, and d- Since the employed rational functions are not limited to have any specified order; it allows finding a function having a perfect fit with the measurements. A comprehensive windings (mechanical) fault detection procedure, which entails the simultaneous identification of all the three faults characteristics (type, level and location), can be carried out by the algorithms based on application of VF. The important windings (mechanical) faults, classified as: AD, RD, DSV and SC are simulated by application of those faults on four different model transformers of a similar size. In the first step, introducing an algorithm based on VF, the type of fault among the four classified mechanical faults, (using a new presented index) is detected. Next, the level and location of the mechanical fault are determined. III. TF IDENTIFICATION BY VECTOR FITTING In principle, a TF approximation (of a given order) can be found by fitting it with a function which is made of the ratio of the two polynomials [5] as shown in equation (1): f s  

am s m      a1s  a0 bn s n      b1s  b0

(1)

Where m and n are the number of numerator and denominator coefficients, ai and bi are the values of ith coefficient for the numerator and denominator, respectively. However, generally the TF of a passive system can be approximated by a rational function f(s) in the form [5] shown by equation (2):

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering

f s  

n

 i 1

ri  d .s  e s  pi

Volume 11, Number 2, 2011 n 1

(2)

The residues r i and poles p i are either real quantities or appear in complex conjugate pairs, whereas d and e are real. The problem at hand is to estimate all of the coefficients (r i , p i , d and e) in equation (2), so that a least squares approximation of f(s) is obtained over a given frequency interval. It should be noted that equation (2), represents a nonlinear problem in terms of the unknowns, because the unknowns p i appear in the denominator. VF solves the problem of equation 2, sequentially as a linear problem in six stages, as follows [5]: Stage1- Specifying the starting poles pi in (2) The starting poles should be complex conjugate with imaginary parts linearly distributed over the frequency range of interest. Each pair is chosen as follows: p=-α±jβ The real part of the poles should be small enough to avoid any ill conditioning problem, i.e., α=0.01β can be a good choice. Stage2- Multiplying f(s) by an unknown function of σ(s) having poles similar with f(s), as shown in the following:  n ri   d .s  e   s  f s   i 1 s  pi    s     n  R   i 1    i 1 s  pi 





(3)

Multiplying the second row in (3) with f(s) and using first row yields the following relation:   n Ri ri  d s e .     1 f ( s )  s  pi   s p i i 1   i 1

f  fit s   d

 s  zk  k 1 n

 s  pk 

n

,  fit s  

k 1

 s  ~zk  k 1 n

 s  pk 

(7)

k 1

Then using equations (5), & (7) the equation (8), is derived as follows: n 1

f s  

f  fit s  d  fit s 

 s  zk  k 1 n

 s  ~zk 

(8)

k 1

Equation (8) is employed for calculating the poles of f(s), it is sufficient to calculate the zeros of σ fit (s). Stage 5- Resolving equation (6) with new poles and repeating the operations until the best approximation is obtained. Stage 6- After achieving a good approximation of f(s), coefficients of ai and bi in equation (1) can be calculated using r i , p i , d and e. IV. TEST OBJECTS (THE FAULTS MODEL TRANSFORMERS) AND MEASUREMENTS In the present investigation, all measurements were executed in the time domain then using FFT acquired the frequency response of TF which is illustrated in [12]. Four model transformers were considered in this study and some tests were performed on them to development of the proposed method for finding type, level and location of the faults, as follow:

n

(4)

Or

f  fit s    fit s  f s 

(5)

Equation (4) is linear in its unknowns r i , d, e and R i . Deriving (4) for several single frequencies, leads to a linear relation of equation (6):

Ak X  Bk

(6)

Where:  1  f sk   f sk   1 Ak   sk 1    sk  pn sk  p1 sk  pn   sk  p1 X  r1    rn d e R1    Rn 

Bk  f sk 

Stage 3- Solving (6) by least square method and identifying the coefficients r i , d and e of (2) Stage 4- Calculation of new poles If the functions in equation (5) are reformed as shown in equation (7):

A. Study of AD As a model transformer for the study of AD a high voltage winding with 31 double inverted discs, where 6 turns are present in each disc, and a four layer concentric low voltage winding, where 99 turns are present in each layer, were used. These particular windings were manufactured for the special experimental purposes, and they correspond to windings of a transformer with a rated voltage of about 10 kV and a rated output of 1.3 MVA. Its specific construction permits a gradual axial movement of the internal layer winding with respect to the outer winding. The test object has 82.7 cm height and therefore a 1cm axial displacement is equivalent of 1.2% displacement. B. Study of RD As a model transformer for the study of RD a high voltage winding with 30 double inverted discs, where 11 turns exists in each disc, and a one layer low voltage winding, having 23 turns are used. The double disc winding has a rated voltage of 10 kV and a rated output of 1.2 MVA. The deformation has occurred on the double disc winding in four degrees, as follow:  Degree 1: The sixth up to the 54-th discs were all radially deformed on one side. Deformation was around 7% of the disc radius  Degree 2: The sixth up to the 54-th discs were all radially deformed on two opposite sides. Deformation was around 7% of the disc radius

25

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering 



Volume 11, Number 2, 2011

Degree 3: The sixth up to the 54-th discs were all radially deformed on three sides with 90º with respect to each another. Deformation was around 7% of the disc radius Degree 4: The sixth up to the 54-th discs were all radially deformed on four sides with 90º with respect to each another. Deformation was around 7% of the disc radius

C. Study of DSV For studying DSV the same model winding that is used in section 4-2 experimented. For this purpose another intact winding that space between its discs is 5mm, is selected. For research about the effect of DSV on TFs, the space of discs is changing to 7.5, 10, 15, 20 and 25mm, then TF of them is measured. To better studying the following states are experimented: State 1 of DSV: space variation between disc 2 and disc 3 State 2 of DSV: space variation between disc 4 and disc 5 State 3 of DSV: space variation between disc 8 and disc 9 State 4 of DSV: space variation between disc 12 and disc 13 State 5 of DSV: space variation between disc 16 and disc 17 Although some of the DSV faults given above seem to be exaggerated and impossible to occur on real transformers, they have been executed in this work in order to accumulate more data about DSV and find out precise information concerning TF variations. Another important point is that the higher distances between discs implemented here are compensable with less circumference of winding covered with the DSV, so that the DSVs carried out in this work could be more similar to real faults in transformers.

H1

L1

L1

H1

L2

H2

VLV

L2

H2

c) TF3 

IH2

IH 2 VH 1

d) TF4 

VLV VH 1

Vinput

Ioutput

I output

e) TF5 

Vinput

Figure 1. Different terminals connections, for sensitivity analysis of transfer functions measurements when examining: a) AD, b) AD, c) RD and DSV, d) RD, and e) SC

The measured results show that the mechanical faults (AD, RD, DSV and SC) affect the TFs and modify them. For example, the measured TFs (TF 1 , TF 3 and TF 4 ) are showed in Figures 2 to 4. 4 0cm Displacement 1cm Displacement 2cm Displacement 3cm Displacement 4cm Displacement 5cm Displacement 6cm Displacement 7cm Displacement 8cm Displacement

3.5 1.5

3

L1

H1

H1

2.5

0.5 0.02

0.04

0.06

0.08

0.1

0.4

2

0.2 0.45

0.5

0.55

0.5

0.6

1.5

1

0.5

0.01

0.1

0.2

0.4

0.3

0.7

0.8

0.9

1

Frequency [MHz]

Figure 2. The measured TF of transformer winding when AD occurred (setup TF 1 ) 1.4 Without Deformation

1.2

Degree1 of Deformation Degree2 of Deformation

1

1.2

0.8

Degree3 of Deformation

0.6

Degree4 of Deformation

0.4 0.2

1

0.02

0.04

0.06

0.08

0.1 0.6

Admittance [mA/V]

L1

Admittance [mA/V]

1

D. Study of SC For the study of SC a high voltage winding with 30 double inverted discs, where 9 turns exists in each disc is used. The double disc winding has a rated voltage of 10 kV and a rated output of 1.2 MVA. All discs of the model winding are provided with a tap, whereby the recording of the impulse voltage distribution is enabled along the winding. The input terminal of the winding is subjected with the impulse voltage, as response signals is the earth current at a selected tap along the winding recorded. A short circuit between two discs is applied to different places along the winding for the experimental determination of the effect of short circuits and their spatial arrangement on the TF. For each individual position of the SC the TF is measured. To investigate the sensitivity of transfer TF measurements, for AD, RD, DSV and SC, different terminal conditions have been studied, as shown in Figure 1.

0.5

0.8

0.4 0.3

0.6

0.2 0.1 0.42 0.44 0.46

0.48

0.5

0.52

0.4

0.2

L2

H2

IH2

I a) TF1  H 2 VH 1

26

L2

H2

V b) TF2  L1 VH 1

0 0.01

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

Frequency [MHz]

Figure 3. The measured TF of transformer winding when RD occurred (setup TF 3 )

1

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering

Volume 11, Number 2, 2011

0.8

  FI    

Without SC SC between units 1, 2

0.7

SC between units 2, 3

Admittance [mA/V]

0.6

m

n

(10)

Where m and n are the number of numerator and denominator coefficients (in equation (1)), a fi and b fi are

0.5

the values of ith coefficient for the faulted conditions, ani

0.4

0.3

0.2

0.1

0

i 1 a fi    i 1 b fi    n  m i 1 ani   i 1 bni 

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

Frequency [MHz]

Figure 4. The measured TF of transformer winding when SC occurred (setup TF 5 )

and bni are the values of ith coefficient for the normal conditions, respectively. This new index is the ratio of sum of the numerator coefficients in faulted condition to sum of the numerator coefficients in normal condition multiply by the ratio of sum of the denominator coefficients in faulted condition to sum of the denominator coefficients in normal condition. The proposed index is used for detection of mechanical faults as follow: 15

In this paper, fault diagnosis in transformer will be done with underneath assumptions: Measurements is carried out on a frequency range of 101000 kHz Measurements is limited to layer- disc type winding of transformer Impact of other two phases on the TF of inspected phase is considered negligible Using the three assumptions introduced in this section, there can be some misleading results only if those assumptions were not applicable in special cases. However to evaluate and identity those cases we plan to extend our verification work to cover some special type transformers in future. Reminding that VF method is one of the best methods for estimating TF, in this paper using the measured data, the TFs are estimated by VF method. A 40-th order approximation has been fitted to these functions. This is employed to develop TFs for possible form of faults and to compare them with the TFs of the related normal conditions. Besides detecting the poles and residues of TFs, the coefficients of TFs’ polynomials in numerator and denominator are determined. A maximum error of less than 10-6 in all cases is the proof of a good fit. References [7-8], introduce one index for comparison of TF using VF, as follows:

10

SDP 

i 1 n

Pdi  Pwi Pdi

(9)

Where n is the number of poles, P di and P wi are the values of ith pole for displaced (deformed) and intact cases, respectively. SDP is the sum of absolute differences of poles and doesn’t perceive the effect of zeros. However, Figure 5 shows that the mechanical faults cause dislocation of both poles and zeros. Therefore to compare the TFs, the effect of both of them (pole and zero) is considered in this method. For this purpose, a new index (FI: Faulted-Intact relation) is introduced as follows:

Imaginary Axis

V. PROPOSED METHOD FOR FAULT DIAGNOSIS

Pole-Zero Map Poles of intact case Zeros of intact case Poles of degree1 of deformation Zeros of degree1 of deformation Poles of degree2 of deformation Zeros of degree2 of deformation Poles of degree3 of deformation Zeros of degree3 of deformation Poles of degree4 of deformation Zeros of degree4 of deformation

5

0

-5

-10

-15 -14

-12

-10

-8

Real Axis

-6

-4

-2

0

Figure 5. Poles and zeros of the estimated TF 4 for study of RD

A. Development of Winding Model for a Layer-Disc Type Winding The test objects that is introduced in section 4, is employed for development of this winding model. Here it will be shown that all of the four fault types (AD, RD, DSV and SC) can be detected by application of the new method. Based on our experiments as it will be shown in next section, for each of the four fault types the FI will preserve a value that is distinguishable from the other type easily. Although this finding can be used in fault detection of similar windings in real transformers, however this method will be extending to cover all types of windings in future. A.1. Determining the Type of Fault The related results demonstrate that the value of FI if AD or RD occurs is less than 1 however, if DSV or SC occurs is more than 1. A careful inspection of Figures 2, 3 implies that the changes in TFs do not spread equally over the entire frequency range, and those are often bounded to several frequency ranges. On the other hand it is clear that by subdividing the entire frequency range into smaller ranges, it can increase the sensitivity of the comparison index, such as FI. In other words, existence of any similarity between the faulted TF and the intact TF depends on the frequency range of interest and the type of fault. Hence, evaluation of the TFs in a set of smaller frequency ranges can be more helpful in determining the type of fault. Examining the simulation results, it has been shown in [12] that the variation of TF characteristics due to RD in the

27

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering

Volume 11, Number 2, 2011

TABLE I. FI LF WITHIN THE RANGE OF LOW FREQUENCY (10-100 KHZ) FOR AD FI lf Displacement (cm) TF 1 TF 2 0 1 1 1 0.996 0.998 2 0.993 0.987 3 0.987 0.985 4 0.983 0.983 5 0.979 0.981 6 0.974 0.976 7 0.970 0.973 8 0.968 0.971 TABLE II. FI LF WITHIN THE RANGE OF LOW FREQUENCY (10-100 KHZ) FOR RD FI lf Degree of Deformation TF 3 TF 4 Without RD 1 1 Degree1 0.908 0.935 Degree2 0.885 0.851 Degree3 0.831 0.793 Degree4 0.753 0.771 TABLE III. CALCULATED FI USING TF 3 FOR DIFFERENT DISC SPACES IN FIVE STATES

State of DSV State 1 State 2 State 3 State 4 State 5

28

Disc Space (mm) 5 1 1 1 1 1

7.5 1.0191 1.0212 1.0252 1.0292 1.0334

10 1.0404 1.0454 1.0488 1.0555 1.0636

15 1.0743 1.0777 1.0803 1.0839 1.0881

20 1.0911 1.1054 1.1091 1.1144 1.1199

25 1.1283 1.1316 1.1397 1.1478 1.1553

5.5 5

Amount of FI

winding occurs in the entire frequency range; whereas the variations due to occurrence of AD are located in the frequency ranges (approximately) above 100 KHz. Consequently, to distinguish the occurrence of an AD fault from an RD fault, FI index should be determined within the low frequency range of (10-100 KHz), which is named FI lf . In tables I and II, FI lf results are determined under different possible faults types of AD and RD, and under different terminal connections. Studying the tables shows that the maximum variation due to occurrence of AD is 3.2% which is less than the minimum variation under occurrence of RD which is 6.5%; moreover these tables show that if the measured TF be admittance function, then the expected minimum variation of RD is 9.2%. Therefore it would be simple to distinguish between AD and RD. Meanwhile as it is shown in Table III and Figure 6, due to occurrence of some other faults, such as; DSV and SC, the value of FI index will be more than 1. Reminding that SC causes a substantial variation, while DSV makes a little variation in the measured TF, it is easy to distinguish between them. Also they can be distinguished by examination of FI index, since the minimum value of FI during SC is 1.5501, whereas the maximum of FI during DSV is 1.1553. It means that during a SC, the minimum variation in FI is 55%, whereas the variation in FI during a DSV is 15.5% at most. Therefore, to specify the type of a fault, the FI index is determined, if this index is more than 1, that fault is SC or DSV and otherwise the fault is AD or RD. Additionally, if FI is less than 1, the FI index in low frequency range should be determined (named FI lf ) to specify the type of fault. However, if it’s more than 1, by determining the value of FI the type of fault can be specified.

SC between discs 2 and 3 SC between discs 4 and 5

4.5 4 3.5 3 2.5

SC between discs 1 and 2

SC between 2 discs which 2 turns short circuited SC between 2 discs which 18 turns short circuited

SC between discs 3 and 4

2 1.5

5

10

15

20

25

30

35

40

45

50

55

Number of Disc

Figure 6. Calculated FI for different locations of SC between discs

A.2. Detection of Fault Level To specifying the fault level, the suggested coefficients are calculated in both faulted and intact conditions within the whole frequency range. Using the value of FI index, the amount of AD, RD, DSV or SC occurred in a winding can be determined. As it is shown in the tables IV, and V any decrease in FI index, corresponds to an increase in the degree of deformation, or axial displacement of the winding. Whereas, Table III shows an increase in FI index is related to an increase in DSV. Meanwhile, Table IV also shows that the values of SDP, for larger displacements (such as 6, 7 and 8 cm), are very closer to each other, and it is very hard to distinguish between these displacements. While, the value of FI index in these conditions are more desirable to be employed. Not only the FI index can specify the type of fault, but also it is a better measure in determination of the fault level than the other index. Similarly, the level of SC can be detected using FI index. As shown in Figure 6, occurrence of SC between 2 discs changes the proposed index. However if SC occurs between 2 discs where 2 turns are short circuited, the FI index changes less than when 18 turns are short circuited. Moreover for all states of SC, the FI index has a unique and specific value for any SC between discs. Therefore if the value of FI is between 1.5501 and 1.9543, the SC between 2 discs is occurred across two turns and if the value of FI is between 2.9688 and 5.3011, 18 turns are short circuited. Therefore, the level of the SC can be determined with the help of the FI as shown in Figure 6. TABLE IV. FI AND SDP WITHIN THE WHOLE MEASURED FREQUENCY RANGE (10-1000 KHZ) FOR AD Displacement (cm) 0 1 2 3 4 5 6 7 8

TF 2

TF 1 SDP 0 0.0042 0.0092 0.0235 0.0388 0.0533 0.0696 0.0722 0.0790

FI 1 0.959 0.836 0.692 0.575 0.493 0.380 0.289 0.197

SDP 0 0.0085 0.0159 0.0314 0.0566 0.1052 0.5048 0.5089 0.5175

FI 1 0.938 0.857 0.772 0.652 0.412 0.274 0.183 0.103

TABLE V. FI AND SDP WITHIN THE WHOLE MEASURED FREQUENCY RANGE (10-1000 KHZ) FOR RD TF 3 TF 4 Degree of Deformation SDP FI SDP FI

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering Without RD Degree1 Degree2 Degree3 Degree4

0 0.1244 0.3083 0.7544 0.8891

1 0.925 0.749 0.525 0.315

0 0.7740 1.0991 2.7110 3.7917

Volume 11, Number 2, 2011 1 0.850 0.732 0.422 0.203

A.3. Fault Localization To demonstrate the ability of FI index (to recognize location of fault) in Table III, its value is given for different locations and levels of DSV. Table III demonstrates FI has a unique and specific value for any level and location of the fault. If the value of FI is (for example) between 1.0743 and 1.0881, the distance between two successive discs of the winding is 15 mm. If the fault level is known, the location of the fault can be determined with the help of the FI index and one column of the Table III for the known fault levels. It means the proposed index, in addition to detection of the type and the level of the fault, can be employed to detect the location of fault. B. Verification of Winding Model In order to evaluate and validate the proposed method, another test object of similar type (layer-disc winding) with a different rating is examined. Specifications of this test object are: - high voltage winding with 40 double inverted discs, where 12 turns are present in each disc - low voltage winding with four layer concentric, where 88 turns are present in each layer A short circuit is applied between two discs of the winding, in two different locations along the winding separately, and the TFs of them (according to fig. 1 (e)) is acquired. Then for fault diagnosis, the proposed algorithm is used. The proposed method determines the values of FI for these places as 1.6602 and 5.1715, respectively. For these values of FI, the proposed algorithm determines the fault to be Short Circuit, correctly. Moreover these result show precisely that the first case (FI=1.6602) is related to a SC occurred between 2 discs where 2 turns are short circuited, and the second case (FI=5.1715) is related to 24 turns being short circuited. It means the proposed algorithm can be employed to detect the type and level of fault in different transformer windings. C. Proposed Algorithm In this paper to specify the type, level and location of the mechanical fault in transformer windings, the following algorithm has been proposed: 1- Obtaining TF of transformer in both intact and faulted conditions using the results of measurements 2- Estimating coefficients of TF in both intact and faulted conditions using VF method 3- Calculation of FI 4- if FI>1 then GOTO step 5 5- if percentage of changes in FI is more than 30% then Fault is SC 6- else fault is DSV 7- end if 8- else Estimating coefficients of TF in both intact and faulted conditions using VF method within the range of low frequency (10-100 KHz) 9- Calculation of FI lf

10- if percentage of changes in FI lf is more than 5% then fault is RD 11- else fault is AD 12- end if 13- end if 14- Specifying the level of AD, RD, DSV and SC using FI index 15- Localization of DSV using FI index VI. CONCLUSION The past methods suggested by other researcher only focus on detection of the level of mechanical fault in transformer and there is no well defined method to specify the type of fault, by comparison of the TFs. Therefore a new method is presented in this work to specify the type, level and location of mechanical fault in windings of transformers with a relatively well defined and accurate procedure. For this purpose, the required measurements have been carried out on four proper Models in both faulted and intact conditions with different degrees of AD, RD, DSV and SC of winding. The results of measurements clearly demonstrate how the mechanical faults modify the TFs. Therefore, using vector fitting method, appropriate coefficients of TFs has been estimated in intact and faulted conditions. Subsequently by introducing a new index based on the calculated coefficients named FI, the TFs are compared. The obtained results show that: a) The value of FI is less than 1 during AD or RD faults, and more than 1 during DSV or SC faults, b) The difference between the FI indexes, for different degrees of faults, during DSV or SC faults is large and can easily distinguish between them, c) Variation in TF in the low frequency range (10-100 KHz) during AD is less than during RD. Thus, to distinguish AD from RD, the FI index is calculated in the low frequency range. The measured results showed that it is easy to distinguish easily between AD and RD faults, d) FI index has shown its effectiveness, even in the determination of the fault level, and e) Results of the application of the proposed index for DSV fault show that the FI index can be used for detection of fault location. In order to evaluate and validate the proposed method, another test object of similar type (layer-disc winding) with a different rating is examined. The related result showed the proposed method is valid for windings which have the same design. Applicability of the proposed method on fault detection of different type of transformers (including several windings of different types) will be studied in the future works, too. ACKNOWLEDGMENT Authors would like to thank Mr. Kharazi and Mr. Firoozi for provision of the second sample winding, permission and assistance in performing the required measurements on the windings. REFERENCES [1]

C. Bengtsson, “Status and trends in transformer monitoring,” IEEE Trans. Power Delivery, vol. 11, no.3, pp. 1379-1384, 1996.

29

[Downloaded from www.aece.ro on Wednesday, June 01, 2011 at 16:04:09 (UTC) by 217.218.226.142. Redistribution subject to AECE license or copyright. Online distribution is expressly prohibited.]

Advances in Electrical and Computer Engineering [2]

[3] [4]

[5] [6] [7]

[8]

[9] [10] [11]

[12]

[13] [14]

[15]

[16]

[17]

30

J. Christian, and K. Feser, “Procedures for detecting winding displacements in power transformers by the transfer function methods,” IEEE Trans. Power Delivery, vol. 19, no.1, pp. 214-220, 2004. T. Leibfried, and K. Feser, “Monitoring of power transformers using the transfer function method,” IEEE Trans. Power Delivery, vol. 14, pp. 1333–1341, 1999. K. Feser, J. Christian, C. Neumann, U. Sundermann, T. Leibfried, A. Kachler, and M. Loppacher, “The transfer function method for detection of winding displacements on power transformers after transport,” in Proc. Short Circuit or 30 Years of Service, CIGRE 12/33-04, 2000. B. Gustavsen, and A. Semlyen, “Rational approximation of frequency domain responses by vector fitting,” IEEE Trans. Power Delivery, vol. 14, pp. 1052–1061,1999. B. Gustavsen, “Wide band modeling of power transformers,” IEEE Trans. Power Delivery, vol. 19, pp. 414–422, 2004. P. Karimifard, G. B. Gharehpetian, and S. Tenbohlen, “Determination of axial displacement extent based on transformer winding transfer function estimation using vector-fitting method,” European Trans. Electrical Power (ETEP), vol. 18, pp. 423-436, 2008. P. Karimifard, G. B. Gharehpetian, and S. Tenbohlen, “Localization of winding radial deformation and determination of deformation extent using vector fitting-based estimated transfer function,” European Trans. Electrical Power (ETEP), vol. 19, pp. 749-762, 2009. Cigre Working Group A2.26, “Mechanical-Condition assessment of transformer windings using frequency response analysis (FRA),” 2007 L. Satish, and S. K. Sahoo, “An effort to understand what factors affect the transfer function of a two-winding transformer,” IEEE Trans. Power Delivery, vol. 20, no.2, pp. 1430-1440, 2005. A. Singh, F. Castellanos, J. R. Marti, and K. D. Srivastava, “A comparison of trans-admittance and characteristic impedance as metrics for detection of winding displacements in power transformers,” Electric Power Systems Research, vol. 79, no.1, pp. 871-877, 2009. E. Rahimpour, J. Christian, K. Feser, and H. Mohseni, “Transfer function method to diagnose axial displacement and radial deformation of transformer winding,” IEEE Trans. Power Delivery, vol. 18, no.2, pp. 493–505, 2003. M. Florkowski, and J. Furgał, “Modelling of winding failures identification using the frequency response analysis (FRA) method,” Electric Power Systems Research, vol. 79, no.1, pp. 1069-1075, 2009. J. Jayasinghe, Z. D. Wang, P. N. Jar, and A. W. Darwin, “Investigations on sensitivity of FRA technique in diagnosis of transformer winding deformations,” in Proc. Int. Symp. Electrical Insulation, Indianapolis, USA, September 2004. E. Rahimpour, and S. Tenbohlen, “A mathematical model to investigate disc space variation in power transformer using transfer function analysis,” in Proc. Int. Symp. High Voltage Engineering, Ljubljana, Slovenia, August 2007. S. M. Islam, and G. Ledwich, “Locating transformer faults through sensitivity analysis of high frequency modeling using transfer function approach,” in Proc. Int. Symp. Electrical Insulation, Montreal, Canada, June 1996. A. Shintemirov, W. H. Tang, and Q. H. Wu, “Transformer winding condition assessment using frequency response analysis and

Volume 11, Number 2, 2011

[18]

[19]

[20]

[21]

[22] [23] [24]

[25]

[26]

[27]

[28] [29]

[30] [31] [32]

evidential reasoning,” IET Elec. Pow. Applications, vol. 4, no.3, pp. 198 – 212, 2010. D. K. Xu, C. Z. Fu, and Y. M. Li, “Application of artificial neural network to the detection of the transformer winding deformation,” in Proc. IET Conf. Int. Symp. High Voltage Engineering, August 1999. pp. 220–223. J. Zhijian, L. Jingtao, and Z. Zishu, “Diagnosis of transformer winding deformation on the basis of artificial neural network,” in Proc. Int. Conf. Properties and Applications of Dielectric Materials, June 2000, pp. 173–176. S. Birlasekaran, Y. Xingzhou, F. Fetherstone, R. Abell, and R. Middleton, “Diagnosis and identification of transformer faults from frequency response data,” in Proc. IEEE Conf. Power Engineering Society Winter Meeting, Jan. 2000, pp. 2251-2256. K. Jong-Wook, P. ByungKoo, J. Seung, K. S. Woo, and P. PooGyeon, “Fault diagnosis of a power transformer using an improved frequency-response analysis,” IEEE Trans. Power Delivery, vol. 20, no.1, pp. 169–178, 2005. S. A. Ryder, “Transformer diagnosis using frequency response analysis: results from fault simulations,” in Proc. IEEE Conf. Power Engineering Society Summer Meeting, July 2002, pp. 399–404. S. A. Ryder, “Diagnosing transformer faults using frequency response analysis,” IEEE Electrical Insulation Magazine, vol. 19, no.2, pp. 1622, 2003. P. M. Nirgude, D. Ashokraju, A. D. Rajkumar, and B. P. Singh, “Application of numerical evaluation techniques for interpreting frequency response measurements in power transformers,” IET Sci. Meas. Technology, vol. 2, no.2, pp. 275–285, 2008. J. R. Secue, and E. Mombello, “Sweep frequency response analysis (SFRA) for the assessment of winding displacements and deformation in power transformers,” Electric Power Systems Research, vol. 78, no.6, pp. 1119-1128, 2008. J. Gui, W. Gao, K. Tan, and S. Gao, “Deformation analysis of transformer winding by structure parameter,” in Proc. Int. Conf. Properties and Applications of Dielectric Materials, paper no. P2-47, 2003. P. Karimifard, and G. B. Gharehpetian, “A new algorithm for localization of radial deformation and determination of deformation extent in transformer windings,” Electric Power Systems Research, vol. 78, no.10, pp. 1701-1711, 2008. R. Wimmer, S. Tenbohlen, M. Heindl, A. Kraetge, M. Krüger, and J. Christian, “Development of algorithms to assess the FRA,” in Proc. Int. Symp. High Voltage Engineering, paper no. T7-523, 2007. E. Rahimpour, and D. Gorzin, “A new method for comparing the transfer function of transformers in order to detect the location and amount of winding faults,” Electrical Engineering, vol. 88, no.5, pp. 411-416, 2005. M. Wang, A. J. Vandermaar, and K. D. Srivastara, “Evaluation of frequency response analysis data,” in Proc. Int. Symp. High Voltage Engineering, 2001, pp 904–907. S. K. Sahoo, and L. Satish, “Discriminating changes introduced in the model for the winding of a transformer based on measurements,” Electric Power Systems Research, vol. 77, no.7, pp. 851-858, 2007. K. Ragavan, and L. Satish, “Localization of changes in a model winding based on terminal measurements: experimental study,” IEEE Trans. Power Delivery, vol. 22, no.3, pp. 1557-1565, 2007.