Open Collections

UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

High-frequency transformer model for switching transient studies Chimklai, Suthep 1995

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata

Download

Media
831-ubc_1995-982683.pdf [ 4.27MB ]
Metadata
JSON: 831-1.0065256.json
JSON-LD: 831-1.0065256-ld.json
RDF/XML (Pretty): 831-1.0065256-rdf.xml
RDF/JSON: 831-1.0065256-rdf.json
Turtle: 831-1.0065256-turtle.txt
N-Triples: 831-1.0065256-rdf-ntriples.txt
Original Record: 831-1.0065256-source.json
Full Text
831-1.0065256-fulltext.txt
Citation
831-1.0065256.ris

Full Text

HIGH-FREQUENCY TRANSFORMER MODEL FOR SWITCHING TRANSIENT STUDIES by SUTHEP CHIMKLAI B.Eng., Chulalongkorn University, 1977 M.Sc. inEE., Carnegie Mellon University, 1984 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES Department of Electrical Engineering We accept this thesis as conforming tQihe required standard  THE UNIVERSITY OF BRITISH COLUMBIA April 11, 1995 © Suthep Chimklai, 1995  In presenting this  thesis  in partial  degree at the University of  fulfilment  of  the  requirements  for an advanced  British Columbia, I agree that the Library shall make  it  freely available for reference and study. I further agree that permission for extensive copying of  this thesis for scholarly  department  or  by  his  or  her  purposes may be granted by the head of  representatives.  It  is  understood  that  copying  my or  publication of this thesis for financial gain shall not be allowed without my written permission.  (Signature)  Department of  £ \ec W  CAL  & A & I r\ € e n Tt <f  The University of British Columbia Vancouver, Canada  Date  DE-6 (2/88)  fepol  Gil),  iWS'  ABSTRACT  The objective of this thesis is to develop a simplified high-frequency model for three-phase, two- and three-winding transformers. The model is an extension of the classical 60 Hz model which includes two important factors prevailing in transformers under transient conditions: stray capacitances which cause transformers to resonate and frequency dependent characteristics of the leakage flux and winding resistances due to skin effects. The model is not aimed to represent internal details of the transformer and only lumped circuit parameters are used in order to simulate terminal behaviours of the transformer. However, it is different from other terminal models in that it is not just an impedance or admittance black box derived from measured transfer functions. Only the meaningful parameters which correspond to the physical components in the real transformer are included in the model. The short-circuit impedances T-form of the classical model is retained which makes it possible to separate thefrequency-dependentseries branch form the constantvalued capacitances. In addition, it enables the model to be built at the coil level which is independent of winding connections. The model stray capacitances are placed at the corresponding coils terminals. If they link two coils they will be split into two halves with one half connected at the upper ends and the other half at the lower ends. The frequency ii  dependent series branch is divided into sections corresponding to various sections in the transformer coil which can be assumed uniform. An RL equivalent network is used to synthesise the frequency dependent behaviour of each section. The values of R's and L's are calculated from minimum-phase-shift approximations which guarantees numerical stability of the resulting network. With the use of symmetrical components, mathematical complications of fitting mutual impedance functions are avoided and also the number of impedance functions to befittedby rational functions is reduced. A number of short-circuit tests on the actual power transformers installed in the Thailand 's power system were performed to determine the parameters of the model. The frequency responses calculated from the model are compared with the tests. Also, a timedomain test was conducted and the result was used for comparison with the simulation from the model.  iii  TABLE OF CONTENTS  ABSTRACT  ii  TABLE OF CONTENTS  iv  LIST OF ILLUSTRATIONS  vi  LIST OF TABLES  x  ACKNOWLEDGMENT Chapter One  Chapter Two  Chapter Three  xii  INTRODUCTION  1  1. Previous Work  3  2. Proposed Model  10  NON-FREQUENCY DEPENDENT PART  14  1. Stray capacitances in Three-Phase Transformer  14  2. Types of Stray Capacitances  16  3. Discrete-Time Model  20  4. Capacitance Model for Multi-Resonance Transformer  23  5. Implementaiton of Interwinding Capacitance Modification  34  MODELLING OF FREQUENCY DEPENDENT COMPONENTS  36  1. Short-Circuit Responses of Transformers  36  2. Equivalent Network Representation for Z,,^^  39  3. Discretization of Equivalent Network for Z , ^ ^  44  iv  4. Updating History Terms Chapter Four  Chapter Five  THE COMPLETE MODEL  50 53  1. Discrete-time Model in Sequence Domain  54  2. Discrete-time Model in Phase Co-ordinates  59  3. Nodal Admittance Matrix  64  MEASUREMENTS AND NUMERICAL RESULTS  68  1. Measurement of Frequency Dependent Leakage Impedance 2. Measurement of Stray Capacitances Chapter Six  NUMERICAL EXAMPLES  68 84 90  1. Two-winding Transformer  90  2. Three-winding Transformer  112  3. Simulation of Transient Recovery Voltage  137  Chapter Seven  PROPOSED FUTURE WORK  141  Chapter Eight  CONCLUSIONS  143  REFERENCE LIST  145  Appendix 1  148  FITTING OF THE IMPEDANCE FUNCTION  v  LIST OF ILLUSTRATIONS Figure 1.1  Page Amplitudefrequencyresponse of a 50 MVA 115/23 kV transformer versus calculated impedancefromits 60 Hz model  2  1.2  Transformer model developed by [2]  6  1.3  P. T. M. Vaessen's highfrequencytransformer model  7  1.4  The proposed transformer model which constists of the non-frequency dependent stray capacitances and thefrequencydependent part  10  2.1  A cutaway view of a single phase, two-winding transformer  14  2.2  Stray capacitances in a three-phase, three-winding transformer  15  2.3  Stray capacitances in any phase of a three-phase, three-winding transformer  17  2.4  Impedance of stray capacitances in the 115/23 kV, 50 MVA transformer  18  2.5  Phase-to-phase capacitances in a three-phase transformer  19  2.6  Discrete time representation of a capacitance  20  2.7  Stray capacitances in a single phase, two-winding transformer  21  2.8  Representation of interwinding capacitances for a two-winding transformer 2.9 Transformation of a three-winding transformer into a delta representation 2.10 Equivalent circuit after interwinding capacitances have been moved to the same side as the leakage impedances vi  24 27 31  Figure  Page  3.1  Typical short-circuit frequency response of a transformer  3.2  Frequency dependent characteristics of the real part of the series impedance Z ^ ^  37  39  3.3  RLC synthesis network to approximate Z ^ ^  44  3.4  Discretization of a section of the equivalent network for Z ^ ^  46  4.1  Discrete-time equivalent circuit for a sequence model of a transformer  54  4.2  Transformer node labelling for external network connection  65  5.1  Measurement of positive-sequence impedance in a Yd transformer  75  5.2  Application of the Capacitance & D. F. Bridge to measure stray capacitances in a transformer Measurement of the sum of interwinding and winding-to-ground capacitances in a three-winding, three-phase transformer with CB100  5.3  85  86  6.1  Transformer capacitance impedance measured with an Impedance Analyser 92  6.2  Schematic of short-circuit frequency response tests  94  6.3  Real part of short-circuit impedance (zero-sequence) for the two-winding transformer Short-circuit impedance (zero-sequence) for the two-winding transformer  95 96  Short-circuit impedances (positive-sequence) of the two-winding transformer measured on different windings  99  6.4  6.5  6.6  6.7  Real part of the short-circuit impedance (positive-sequence) for the two-winding transformer  100  Transformer short-circuit impedance (positive-sequence) for the two-winding transformer  101  vii  re 6.8  Fitting of the real part of the Zv/iIl^Il^ impedance for the zero-sequence circuit of the two-winding transformer  104  Fitting of the magnitude and phase angle of the frequency dependent Z , , ^ ^ impedance for the zero-sequence circuit of the two-winding transformer  105  Fitting of the real part of the Zwinding impedance for the positive-sequence circuit of the two-winding transformer  108  Fitting of the magnitude and phase angle of the frequency dependent Z ^ ^ impedance for the positive-sequence circuit of the two-winding transformer  109  6.12  Impulse response test  111  6.13  Effective stray capacitances of a three-winding transformer during the short-circuit tests  114  Fitting of the short-circuit impedance measured on the low-voltage side with the high and tertiary-voltage windings short-circuited (zero-sequence)  115  Fitting of the short-circuit impedance measured on the high-voltage side with the low and tertiary-voltage windings short-circuited (zero-sequence)  116  Fitting of the short-circuit impedance measured on the low-voltage side with only the tertiary-voltage winding short-circuited (zero-sequence)  117  Fitting of the short-circuit impedance measured on the tertiary-voltage side with the high and low-voltage windings short-circuited (positive-sequence)  118  Fitting of the short-circuit impedance measured on the tertiary-voltage side with only the high-voltage winding short-circuited (positive-sequence)  119  6.9  6.10  6.11  6.14  6.15  6.16  6.17  6.18  viii  Figure 6.19  Page Fitting of the short-circuit impedance measured on the tertiary-voltage side with only the low-voltage winding short-circuited (positive-sequence)  120  Fitting of the sum of the high and low-voltage winding impedances (zero-sequence)  124  Fitting of the sum of the high and tertiary-voltage winding impedances (zero-sequence)  125  Fitting of the sum of the low and tertiary-voltage winding impedances (zero-sequence)  126  Fitting of the sum of the high and low-voltage winding impedances (positive-sequence)  127  Fitting of the sum of the high and tertiary-voltage winding impedances (positive-sequence)  128  Fitting of the sum of the low and tertiary-voltage winding impedances (positive-sequence)  129  Fitting of the winding impedance of the high-voltage winding (positive-sequence)  136  Winding impedance of the low-voltage winding (positive-sequence)  138  Winding impedance of the tertiary-voltage winding (positive-sequence)  139  6.29  Transient rocovery voltage in a circuit breaker after clearing a fault  140  7.1  Capacitance model with dielectric losses  141  Equivalent network representation of the frequency dependent branch of the sequence model  148  6.20  6.21  6.22  6.23  6.24  6.25  6.26  6.27  6.28  Al. 1  ix  LIST OF TABLES Table 1.1  1.2 3.1  Page Comparison of the no. impedance functions required by the proposed model and the Ontario Hydro's model for a three-phase transformer  4  No. of history terms to be updated in the proposed model and in the Ontario Hydro's model  5  Summary of discrete-time models for inductances and capacitances using trapezoidal and backward Euler integration rules  45  6.1  Measured values of capacitances for the tested transformer  91  6.2  Parameters of the zero-sequence short-circuit equivalent network for the tested two-winding transformer Parameters of the positive-sequence short-circuit equivalent network for the tested two-winding transformer  97 102  Parameters of the zero-sequence Z,,^^ impedance for the tested two-winding transformer  106  6.3 6.4 6.5  Parameters of the positive-sequence Zvlindjng impedance for the tested two-winding transformer  110  6.6  Measured values of capacitances for the three-winding transformer  113  6.7  Parameters of the zero-sequence short-circuit equivalent network for the tested three-winding transformer Parameters of the positive-sequence short-circuit equivalent network for the tested three-winding transformer  121  6.8  x  122  6.9  Parameters of the zero-sequence Z ^ ^ impedance for the tested three-winding transformer  6.10 Parameters of the positive-sequence Z ^ ^ impedance for the tested three-winding transformer  xi  130 131  ACKNOWLEDGMENT  I would like to thank the Electricity Generating Authority of Thailand (EGAT) which granted me the opportunity to come to Canada on an education leave. I would like also to thank Mr. Siridat Glankhamdee, Director of Systems Planning Department, for his very much appreciated guidance regarding my research work. A very special thank you to the Canadian International Development Agency (CIDA) for providing me the scholarship for my graduate study at UBC. My deepest gratitude to Dr. Jose R. Marti, my advisor, for initiating the research idea and for his supervision afterwards. During the entire course of my study, I received invaluable suggestions, support and time from him. With his constant encouragement, my achievements are even more special to me. To all the professors whose courses I attended. I learned a lot from them during the time I did my coursework. I found their lectures extremely worthwhile for my research work. To my special friend Supatra for her excellent work on transformer data which was a key to the success of my project. I will always hold this exceptionalfriendshipin high regard. To all myfriendsin the Power Group who in one way or another contributed to the accomplishment of my study, I sincerely thank all of them. Last but not least, I want to thank my wife Porntip for her understanding and patience. Without these, I would have not been able to stay in Canada until the end of my study program.  xii  Chapter One INTRODUCTION  Switching transients have been a concern in most power utilities. They could cause damages to various equipment in the power system. The problem may be prevented or avoided if simulations on the power system under various transient conditions can be conducted with sufficient accuracy so that the effects of the transient can be predicted in advance. However, simulation studies need good models of the electrical equipment including power transformers which are very common in the power system. Unfortunately, a mathematical model of the power transformer is not common, and a general model is being under research by many interested people. When the power transformer is in operation, it will be subject to both steady-state and transient conditions. At steady state, the coupling in the transformer is primarily due to the inductive interaction of various coils. Under this circumstance, the transformer may be accurately modelled with inductances and resistances without the need to include any capacitances. The low frequency model is good for an electromagnetic transients study if the frequency components under consideration are below a certain frequency limit. For example, the low frequency model of the 50 MVA, 115/23 kV power transformer whose data are used in this thesis will be good to about 5 kHz, as shown in fig. 1.1. At a higher 1  20 k -  IT o  •§  II  -  /^Constant 60-Hz -  10 k Measured Response  ^  /  >  _ VV_  . v ^  0 100  Ik  10 k  100 k  1M  Frequency (Hz)  Fig. 1.1 Ampitude frequency response of a 50 MVA 115/23 KV transformer versus the calculated impedancefromits 60 Hz model.  frequency range, such as switchingfrequencies,the lowfrequencymodel should be modified so as to improve its accuracy. In recent years, many authors have presented a number of good transformer models. Nonetheless, some of these models are limited to single phase transformer while some others are "black box" representations. The problems which are mostly encountered in transformer modelling are numerical stability, difficulty in data acquisition to realise parameters for the model, size of the matrix to describe the model, and complexity of the model itself which requires sophisticated knowledge of electromagnetic theory. This implies that further improvements are still required before satisfactory models finally 2  become available. The proposed model in this thesis is designed to be simple, easy to implement and requires minimum understanding of electromagnetic theory, yet it is sufficiently accurate. The model can be used in various analyses associated with switching transients such as surge arresters selection, insulation level evaluation, etc. In addition, the model may be used to carry out simulations under different transient conditions to investigate causes of problems which already occurred or are expected to take place without the need to conduct the actual field tests. Before further details of the proposed model are discussed, it is helpful to look back at some of the previous works in transformer modelling.  1. Previous Work It may be accurate to say that there are hundreds of transformer models which have been proposed for transient simulation studies. However, the objective of this thesis is on improving EMTP models. Therefore, transformer models will be grouped into EMTP models and models belonging to other sources.  1.1 EMTP A high frequency model for EMTP has been recently developed by Ontario Hydro, Canada [1] . This model was derivedfromthe nodal admittance matrix equation which relates transformer currents and voltages at its external terminals. For example, a wye-deltawye transformer in which the delta winding has no accessible terminals is modelled by a  3  [6x6] matrix equation. This representation is in contrast with the proposed model in this thesis in which the modelling of transformers begins at the coil level. In Ontario Hydro's model, the coupling between different windings of the same phase is averaged so as to make use of the symmetrical components transformation. The distinct elements of the admittance matrix after transformation which are needed to be fitted by rational functions are reduced to m(m+l) from 3m(3m+l)/2, where m is the no. of windings. Comparison of the no. of impedance functions required by the proposed model and Ontario Hydro's model is given in table 1.1 below.  Table 1.1—Comparison of the no. of impedance functions required by the proposed model and Ontario Hydro's model for a three-phase transformer No. of windings  Proposed model  Ontario Hydro's  Times  2  2  6  3  3  6  12  2  Even though it is required to fit only m(m+l) elements, the total number of measurements required is still 3m(3m+l)/2. Therefore, the transformation only saves fitting time not the time spent for measurements. In terms of time-domain solution costs, the differences are even greater. For example, for a three-phase two-winding transformer, only 3 history terms are needed to be  4  updated in the proposed model. Whereas in Ontario's Hydro model, 12 history terms must be updated due to the use of a pi-equivalent circuit which requires 12 impedances for the complete three-phase two winding transformer. The savings of simulation time in the proposed model is clearly demonstrated in table 1.2.  Table 1.2~No. of history terms to be updated in the proposed model and in Ontario Hydro's model No. of windings  Proposed model  Ontario Hydro's  Times  2  3  12  4  3  9  27  3  1.2 Other Models There are a number of models developed by other sources. However, they may be classified into 2 major groups: the terminal models and the detailed models.  1.2.1 Terminal Models Most of the previous works in transformer modelling fall into this type. Some of these works are: (a) M. D'Amore, et al., [2] modelled a transformer as a two port network consisting of a series impedance Z(s) and a parallel capacitance connected at the output port, as 5  Z(s) -o  V.in  V  out  Fig 1.2 Transformer model developed by [2].  shown in fig. 1.2. A number of networks for Z(s), made up of R, L and C's were used to reproduce the voltage gain functions of single and double resonance transformers. There was no restriction in the sign of the R, L, C parameters used. Consequently, there is no guarantee of numerical stability in the model. Furthermore, the choice of Z(s) which would produce more than two resonances was not shown in the work. (b) P. T. M. Vaessen [3] also constructed a transformer model from the voltage function. He used a two-port network consisting of two ideal transformers to duplicate each peak in the voltage gain function. Many such two-port networks were connected in parallel if multiple peaks were to be matched. The diagram of Vaessen's model is shown in fig. 1.3. (c) R. C. Degeneff [4] proposed a transformer model which was made of RL parallel components as the basic building block. The model was constructed by joining a pair of transformer terminals with one RL parallel block and also from each terminal to  R o-  -o  a  a  R,  Vi  X  XX  X  l:m  l:n  X  XX  XX  Fig. 1.3 P. T. M. Vaessen's high frequency transformer model.  ground. In his paper, he demonstrated the model for a 250 MVA single phase auto-transformer with tertiary winding. For this single phase transformer where the common winding was grounded, there were a total of five terminals including ground to be considered. The total number of RL parallel branches required to build the model was ten, or 0.5(n-l)n where n is the number of terminals in the model. This type of model would be almost impossible to realise for the three-phase transformer because the number of branches required would be too many. Furthermore, losses were not included in the model. (d) T. Adielson et al. [5] developed a transformer model from a matrix equation which describes a multi-winding transformer as n magnetically coupled coils. The model was designed for transient studies within the frequency range of tens of kHz. Adielson  proposed the separate modelling of stray capacitances from the frequency dependent part of the transformer model. However, he used an RLC equivalent circuit to represent the frequency dependent part without consideration on numerical stability. He stated in his paper that he encountered many cases of numerical problems in the time-domain simulations. (e) CIGRE WG 33.02 Report [6] presented several models for transient studies. Each model was designed to serve the study within a particular frequency range. The effects of stray capacitances, frequency dependency of the windings and leakage flux were introduced. These models, therefore, were more accurate for transient studies than the 60 Hz model. However, the frequency dependent effect of the winding was represented by only RL parallel branches which is not sufficient to represent a transformer with multi-resonance characteristics. The models described in (a) to (e) above were developed for the single phase transformer. Most of them, except in (e), were based on the concept of the black box, without paying attention to the internal conditions of the transformer under transients. We believe the transformer model should be more representative of the real transformer. It should be general and every parameter constituting the model should be included only if it is related to what it would represent in the real transformer. This is the core idea of the model proposed in this thesis.  8  1.2.2 Detailed Model Not many works are devoted to modelling the transformer in detail. Two of these works are reviewed here: (a) R. C. DegenefF [7] has worked at modelling a transformer at the turn level. Each turn of the coil was represented with a loss resistance, a turn-to-turn capacitance and an inductance with loss. Also the capacitive coupling of the turn to ground was taken into account and was represented with a capacitance with loss. The parameters for the model were obtained from measurements. (b) Francisco de Leon and Adam Semlyen [8, 18], instead of finding the transformer parameters from measurements, applied electromagnetic field theory to calculate the winding leakage inductance and capacitance on a turn basisfromthe knowledge of the physical layout of the transformer. Not only was the leakage impedance modelled in detail, the iron core was also mathematically represented and included in the overall model. Models of this type are good for the study of internal winding stresses. The difficulty which may be encountered is that detailed transformer data are not normally made available by the manufacturers. Also, the size of the matrices involved is too large to make it practical to interface the transformer model with the EMTP.  9  CHL  "HL  Fig. 1.4 The proposed transformer model which consists of the non-frequency dependent stray capacitances (outside the dotted frame ) and thefrequencydependent part (inside the dotted frame).  2. Proposed Model The proposed transformer model could be classified as belonging to the group of "terminal models" but the principles upon which the model is developed are totally different. The model is not a direct translation from the admittance matrix equation nor a black box. On the contrary, the admittance matrix equation could be derived from the model. The model is made to be representative of the real transformer in a simplified way. That is, the real physical components in a transformer under transient conditions are incorporated into the model with the use of lumped R, L and C components. The two main factors influencing the behaviour of the transformer at high frequencies which are - 1) the stray 10  capacitances and 2) the frequency dependence of the leakage flux and losses due to skin effect are taken into consideration. These two components are separately treated in the proposed model which can be seen in fig. 1.4. The magnetising effects such as saturation and hysteresis in the core are not considered. Under switching frequencies, these effects have a lesser influence and can be disregarded [5]. Stray capacitances are assumed to be constant and are represented in the model by non-frequency dependent capacitances connected between the outside terminals of the model. Unlike other previous models in which "mathematically equivalent" capacitances are connected between any two terminals possible, only the physical capacitances are considered in the proposed model. No fictitious capacitances are then present to complicate the model. Leakage fluxes and losses are represented in the same way as it is found in the traditional transformer model. They constitute a series branch of resistances and inductances but both elements are now frequency dependent. In the model, this frequency dependent branch is subdivided into a number of sections, each associated with a resonance region in the transformer response. Each frequency dependent section is approximated with a rational function of the minimum-phase-shift type. With the use of symmetrical components to decouple the model into three sequence-models, direct fitting of the mutual impedances is avoided. The mutual impedances are non-minimum-phase-shift functions and very strict conditions regarding the location of poles and zeroes of  11  approximating rational functions would have to be satisfied to prevent numerical stability problems [13]. RLC equivalent networks, each consisting of a capacitance connected in parallel with a branch of a resistance connected in series with an RL parallel circuit, are used to synthesise the rational functions. The resistance of the series branch of the RLC equivalent network corresponds to the dc resistance of each coil section the network is representing. For the frequency range up to the first resonance peak, several RL parallel circuits are used in the section so as to reproduce the transition from low frequency to the first resonance peak more effectively. In the real transformer, there are capacitances between each turn and between turns on one coil to the others. The former are lumped at the coil terminals and the latter are lumped across each section in the model. Detailed discussion on this issue will be presented later in chapter 2. Modelling the transformer this way has the following advantages: 1) The model is simple yet representative of the real transformer. 2) The model is numerically stable because only minimum-phase-shift functions are considered. 3) The number of impedances to be synthesised is minimal. Only two impedances are required to be fit in a three-phase, two-winding transformer and only six for a three-winding transformer. This can be accomplished with the use of symmetrical components.  12  4) The representation is independent of the particular external connection among windings (wye, delta, etc.) The connection is taken care of by node labelling at the EMTP level. In the proposed model, short-circuit tests and capacitance measurements are needed for calculations of the model parameters. For this thesis, the experiments were carried out at the high-voltage laboratory of the Electricity Generating Authority of Thailand. The type of experiments needed for the model were proposed by the author.  13  Chapter Two NON-FREQUENCY DEPENDENT PART  1. Stray Capacitances in Three-Phase Transformer The non-frequency dependent parts of the proposed transformer model are the stray capacitances which exist in every transformer. Fig. 2.1 shows a simplified cutaway view of a single-phase, two-winding transformer. The high voltage and the low voltage windings are shown as two concentric cylinders. In this arrangement, there is capacitive coupling between the high voltage and the low voltage coils as there would be between  Low Voltage High Voltage  Fig. 2.1 A cutaway view of a single phase, two-winding transformer. 14  any two concentric metallic cylinders. Also, the cylinders and the ground (core and tank) act as capacitors and there will be capacitancesfromtransformer coils to ground. There is another type of stray capacitance which is not shown in fig. 2.1. This capacitance is the turn-to-turn capacitance which is distributed throughout each transformer coil. In the case of three phase transformers, there are also capacitances between outer windings on different phases. Ideally, the stray capacitances for a three-phase, three-winding transformer will be as shown in fig. 2.2 which is drawn to represent a transformer with three-limb core. In each phase, there exist capacitances between the inside winding to the core, the  Phase 1  "^(Tank)  =F  Phase 2  =J=  Phase 3  (Core)  Fig. 2.2 Stray capacitances in a three-phase, three-winding transformer. 15  middle winding to the inside and the outside windings, and between the outside winding and the transformer tank, as well as between the outside coil on the middle limb to the corresponding coils on the other two outer limbs. In addition, there exists additional capacitances which are not shown in fig. 2.2. For example, there are some small capacitances from the outside coil to the inside coil of the same phase and also between the outside coils of phase 1 and phase 3 as well.  2. Types of Stray Capacitances Within the scope of this thesis, the capacitances will be lumped into the following four groups: (1) Capacitances between windings on the same limb (phase) or "interwinding capacitances" ( numbered 2 in fig. 2.2), (2) Capacitances between windings and the tank and between windings and the core or "winding-to-ground capacitances" (numbered 1 in fig. 2.2), (3) Capacitances between turns of the same coil or "turn-to-turn capacitances" (numbered 3 in fig. 2.2), (4) Capacitances among the outer coils on different  limbs (phases) to the  windings on the other two limbs or "phase-to-phase capacitances" (numbered 4 in fig. 2.2). The turn-to-turn capacitances are modelled as a single lumped capacitance connected between the two ends of the coil. The rest of the capacitances are split into two  16  equal parts. As shown in fig. 2.7, the winding-to-winding capacitances are connected between each pair of coils, while the winding-to-ground capacitances are connected from each end of the coil to ground. An alternative way of splitting the winding-to-ground capacitances has been worked out by Allan Greenwood [9]. He has shown that if one end of the coils is grounded, the effective winding-to-ground capacitance will be reduced to only one-third of the total value instead of one half. This idea was implemented in the software developed for the model. If, at any instant in time, a transformer coil is grounded, the  Fig. 2.3 Stray capacitances in any phase of a three-phase, three-winding transformer (turn-to-turn capacitances not shown).  winding-to-ground capacitances of that coil, which were previously composed of two equal parts of value one-half each, will both be changed to one-third.  17  If the three windings of a three-phase transformer are labelled as the high-voltage, the tertiary-voltage and the low-voltage windings, denoted by H, T and L respectively, the averaged stray capacitances which are formed between H, T, L and ground in any of the three phases will be as shown in fig. 2.3. Even though they belong to different phases, these capacitances are assumed to be equal (to their averaged values) so as to achieve a  I  —  I  I  "•*—  I  100 k 90 10 k  45  o  <u  1 3  Ik  a a>  "3>  -45  5<u  -90  cd J3 PL,  Phase Angle  CO  100  10  100  Ik  10 k  100 k  1M  Frequency (Hz)  Fig. 2.4 Impedance of stray capacitances in a 115/23 kV, 50 MVA transformer.  balanced condition. For the same reason, the phase-to-phase capacitances, shown in fig. 2.5, are assumed also to be equal among each other.  18  (a)  (b)  (c)  Fig. 2.5 Phase-to-phase capacitances in a three-phase transformer, (a) unconnected windings, (b) delta-connected windings and (c) wye-connected windings.  In the proposed model, all stray capacitances are assumed to be constant within the frequency range of interest (which is sufficiently high to cover switching transients). Although capacitance isfrequencydependent to some extent, the idea of modelling stray capacitances as constant is not overly simplified. An analysis of the measured frequency response of stray capacitances shows that they are basically constant over a wide frequency range. Figure 2.4 shows an experimental measurement of the combined winding-to-ground and the winding-to-winding capacitances of a 115/23 kV, 50 MVA transformer. As can be seen in this graph, the capacitance is constant to about 100 kHz. The irregularities beyond 100 kHz may probably be due to inaccuracies in the measurements or to local resonances in the winding at high frequencies.  19  Mi-C  ir(t)  (t-At)  & ^> -vc(t)  o  W\A-  R,eq-C  vr(^  Fig. 2.6 Discrete-time representation of a capacitance.  Stray capacitances are modelled directly in the phase domain without the use of symmetrical components. Although, there is mutual coupling of stray capacitances between the three phases due to the phase-to-phase capacitances, it is not necessary to decouple them. Each of these capacitances can be represented as two equal capacitances and connected as shown in fig. 2.5 for wye or delta configuration of the outermost windings.  3. Discrete-Time Model In the electromagnetic transient program, such as the EMTP, capacitances must be represented in discrete time form for a numerical step-by-step solution of the power system network with digital computers. The discrete-time model of a capacitance, shown in fig. 2.6, consists of a resistance and a current source connected in parallel. The values of the resistance and the current source depend on the integration rule used to discretize  20  the capacitance. If the trapezoidal rule is used, the resistance and the current source will be given by  Req-C  A*. 2C  -  (2.1)  -HL  (a)  -LG  N,  N,  © =r  cL  -LG -HL  (b)  -HL  Fig. 2.7 Stray capacitances in a single phase, two-winding transformer, (a) Continuous time model, and (b) Discrete time model. 21  and,  In-cit-At)  =  -2£.v c (/-A/)-»c(f-A/)  (2.2)  where, vc(t) and icO)  = voltages across the capacitance and current flowing into the capacitance, respectively,  Req-c  =  Ih-cif - At)  = known current source.  equivalent internal resistance of the current source,  If the backward Euler is used, the parameters for the discretized capacitance will be,  Req-C = -g  (2.3)  and, 7*-c(/-A/)  =  _£.vc(/-A*)  (2.4)  An admittance matrix equation for the entire stray capacitances network can be formulated in the regular manner. Each diagonal element of the admittance matrix comes from the sum of all the admittances connected to that node and the off-diagonal elements are the negative of the admittances joining a pair of nodes. The external current entering a node will be added to the current vector on the right-hand side of the equation. Figure 2.7 illustrates the discrete-time model for the stray capacitances of a single phase, two winding transformer. C^is half the total winding-to-winding capacitance, CHGis half the  22  total winding-to-ground capacitance for the high-voltage winding and so is C ^ for the low-voltage winding. CH and CL are the total turn-to-turn capacitances for the high and low-voltage winding, respectively. The model of fig. 2.7 is used to show how a nodal equation for stray capacitances can be formulated. An equation for the more complicated stray capacitances configuration of a three-phase transformer can be worked out in a similar fashion. The admittance matrix for the network of fig. 2.7, using the trapezoidal rule of integration, is as follows: CHG  + CHL + CH  2_ At  -CHL -CH 0  -CHL CHL  + CLG + CL 0  0  -CH  0  -CL  CHL + CHG + CH  -CHL  -CHL  CHL+CLG + CL  -CL  -Ihi-oit - At) - Ihx-^it - At) - Ihi-i.it - At) Ih\-i(t - At) - Ih2-4(t - At) - Ih2-G(t - At) Ih\-zif-At) - Ihi-G(t-At) - Ifa-4(t-At) Ihi-A^t - At) + Ih3-4(t - At) - Ih^-oit - At)  vifl) v2(t) V3© V40  (2.5)  where, VI©  = voltage at terminal i of the transformer at the present time step,  IhH(f-At)  =  current source at the previous time step resulting form discretizing the capacitance joining nodes / and j ,  At  =  duration of the time step.  4. Capacitance Model for Multi-Resonance Transformer It should be noted that the capacitance model shown in fig. 2.7 produces only a single resonance during the short-circuit tests. In the short-circuit test, the transformer  23  -HL  (a) Hx  -'leak  NI:N2  n=  -» N,  'HL ^-HL  0>)  5  Hx O  C H L  l  H P2P ^ 'leak  N. Jt!)  'HL  BL  -HL  Fig. 2.8 Representation of interwinding capacitances for a two-winding transformer, (a) Original model and (b) The proposed modification.  model is reduced to an equivalent network consisting of a capacitance in parallel with the leakage impedance, R(a))+j©L(w), which has only one resonance. It is more likely, however, for a transformer to have multiple resonances. These are due to non-uniform winding parts produced, for instance, by the presence of auxiliary tap changing windings [19]  24  and by voltage grading schemes [20], therefore, modifications are required to take this effect into consideration. In the two-winding transformer, the capacitances between the high voltage and low voltage winding can be moved to the side where Zleak is located without altering the terminal characteristics of the model. The transferred capacitance can then be combined with the leakage impedance and treated as a single unit. This process is illustrated in fig. 2.8b.  4.1 Two-Winding Transformer A matrix equation to describe the transformer of fig. 2.8a, in which the interwinding capacitance between the two coils is represented as two equal parts - each with half of the total capacitance value, is 1 Z  hak  "/*,"  hx  =  IH2 .***  +JG>CHL  Z  c + c 5r> ^ i > ^ -i n  Zleak  .  -1  J*--J®CHL  Z  hak  Uak  n  zh. ltak  Z  Uak  Zleak  ^h+JnCHL Zleak  ztl leak  TZ-J^HL Zltak  Z  Z  n  n u<* -rP-  z  (2.6) TTZ-JVCHL Zuak  VH2  -jh+jaCHL Zleak  where, =  currents entering respectively,  VHx, VH2, VLi and VL2  =  voltages at nodes Hi, H 2 , Li and L 2 , respectively,  CHL  =  half of the interwinding capacitance between the high and low-voltage windings.  I Hi, I Hi, hi  and/ij  Similarly, from fig. 2.8b,  25  nodes  Hi, H 2 , Li and L2,  ^+-/(°—  "/*.~  2^-^  C  ^  Z^-^—  ZZ+JOC**  hx  -zh-J®CHL £:+J(onC«L zz+JaCu.  IH2  £-•** i + > c - i + >^ e->c-  I Li  . ^  ^-J^CHL  0  »CHL  0  C//£ - TICHL  -CHL  T;CHL  -CHL  -CHL  nCHL  1  -w Z  £-J®"CHL  > < ^  CHL-J;CHL  + J®  '  +  leak  -n  +j®CHL  Z  -JG>CHL  Z  leak  n2 kak  0  Z  Z  fca*  Z  kak  n leak  VHX  »CHL  VLX  0  VH2  CHL-HCHL  ^  2  VL2  VL2  n 2leak  -n2 Z  n  fco*  -CHL  -1  +j®CHL  Z  £+M>CHL  CHL-^CHL  -JG>CHL  -1  £-JmCHL  1 Z  fca* -W  Z  kak  fra*  +j®CHL -J®CHL  -n Z  ltak  -J(OCHL  JL. +JG>CHL Z leak  VLX VH2  (2.7)  VL2  Equation 2.7 is exactly the same as equation 2.6. This indicates that the transformer model offig.2.8b is identical to that offig.2.8a with parameters as shown. Next, C^/n and Zleak, which are now connected in parallel, will be combined into a new frequency dependent impedance which consists of R, L and C and will be, thereafter, referred to as 117  ii  winding •  4.2 Three-Winding Transformer The same idea can be applied to a three-winding transformer but it takes more steps to get to the final result. The transformer without stray capacitances, which is  26  (a)  (b)  (c)  I  Fig. 2.9 Transformation of a three-winding transformer into a delta representation, (a) Transformer is modelled as three coupled coils, (b) Leakage impedances are referred to the VH side, (c) Leakage impedances are transformed into a delta representation consisting  27  initially modelled as a "T-circuit", must be transformed into a delta representation. The transformation process is shown in fig. 2.9. The model in fig. 2.9c can be described by a matrix equation, where all voltages and currents are associated with the winding, as follows:  1  ZffT  IH  1_  l  +  ZHL  1  •IT  ZHT  1 ZHL  ZHT 1 L + ZHT ZTL  1  1  ZTL  •<HL  VH  1  %'VL  -J- + J ZHL  (2.8)  %-VT  ZTL ZTL  or.  J_ IH IT II  ZHT  +  J_ ZHL  »»_ 1 T  "  na  1 '" ' ZHL L  ZHT  4 "T ZHT "T \-ZHT ZTL-'  I  "H  ~"L  7 HL  nTnL  T L  ~" " 'ZTL  (i  1 •  VH  ^  VT  .+ n  VL  (2.9)  4'^ZHL ZTLJ  where, IH, IT and II  = currents in high, tertiary and low-voltage windings, respectively,  VH, VT and VL  = terminal voltages of high, tertiary and low-voltage windings, respectively,  ZHT, ZTL and ZHL  = leakage impedance of high, tertiary and low-voltage windings, as represented in the delta configuration (fig. 2.9c),  rtH, nT and nL  =  no. of turns of high, tertiary and low-voltage windings, respectively.  28  Equation 2.9 is not the nodal admittance matrix equation, therefore, it cannot be combined with the nodal equation indicated previously to describe the interwinding capacitances. However, a change in variables to express all the voltages and currents as node quantities takes only a small additional effort.  Based on equation 2.9, the nodal  admittance matrix equation for a three-winding transformer can be written as:  VHX  /*, ITX IH2  1  [YD]  -[YD]  =  h2  VTl  vHl  (2.10)  [YD]  -[YD]  l^ J  VL2  where, [YD]  1 and 2  =  admittance matrix of size [3 x 3 ] in equation 2.9, designate variables associated with node 1 and node 2 of each winding. For example, / # , and IH2 are the currents entering nodes H i and H 2 , where / # , = IH, IH2 = -IH- Note that after the interwinding capacitances have been included, IHX is no longer equal to the negative of IM2 .  The interwinding capacitances can be added to the transformer in fig. 2.9c in a similar fashion. This is shown in fig. 2.8a_ where half of the capacitance is connected from "node 1" of one winding to "node 1" of the other winding, and the other half is connected between their corresponding "node 2". After the interwinding capacitances are included, equation 2.10 becomes  29  /*, IT,  [YD]  -[YD]  -[YD]  [YD]  vTl  K VH2  IH2  h2  VT2  J* CHT  + CHL  -CHT —CHL  0 0 0  VL2  0 0 0 0 0 0 CHT + CTL -CTL 0 0 0 -CTL CHL + CJL 0 0 CHT + CHL —CHT —CHL 0 0 -CHT CHT + CTL -CTL 0 0 —CHL —CTL CHL + CTL -CHT  —CHL  vTi vLl  (2.11)  VH2  vTl vrj.  vhere, CHT,  CTL  one half of the high -to-tertiar  and CHL  high-to-low interwinding capacitances, respectively. We consider now taking a portion of the interwinding capacitances and lumping them with the delta-equivalent leakage impedance, as shown in fig. 2.10a. In so doing, some extra capacitances are generated as by-products which must also be included in the equivalent circuit in order to retain similarity between the modified and original network (as shown in fig. 2.10b). All these fictitious capacitances can be chosen properly to be as follows: C'  "T  r  »L  r  ^HT  r' ^HL  r' ^TL  nTnL 2  nH  and, 30  r, (-TL  According to this choice of capacitances, a matrix equation describing the transformer of fig. 2.10a can be expressed as:  (a)  (b) H,  H,  Fig. 2.10 Equivalent circuit after the interwinding capacitances have been moved to the same side as the leakage impedances, (a) Capacitances which will be lumped with the leakage impedances and (b) Additional capacitances resultingfromthe modification. 31  1HX IT.  [YD]  VT,  -[YD]  IL,  VLX  IH2  VH2 -[YD]  IT2  VT2 VL2  [YD]  J* . -[Yo]  [YD]  +  VTX  j®  VH2 -[Y'D]  VT2  [YD]  VL2  [YCl]  VTl  [Yc2]  + >  (2.12)  VHl [Yc2]  [YCl]  VT2 VL2  where, CHT  + CHL HT  [YD]  —fiJ^HT n2 ~T(CHT + CTL)  HL  TL  -HZ^HL nTnL *- TL  „2 (PHL "L  -CH  -CHT  -CHL  -CHT  -CT  -CTL  -CHT  -CTL  -CL  32  +  CTL)  c'H+cHL+cHT [YCl]  o  0  0  CT + CTL + CHT  0  0  0  C'L+CTL + CHL  and -CH  [Yc2]  -CHT —CHL  -CHT  CT  -CTL  -CHL  -CTL  -CL  The first term on the right-hand side of equation 2.12 comes from the leakage impedances Zjjp Z ^ and Z^ and is identical to the first term on the right-hand side of equation 2.11. The second and third terms on the right-hand side of equation 2.12 come from the interwinding capacitances. Equation 2.12 will be identical to equation 2.11 only if the sum of these two terms is the same as the second term on the right-hand side of equation 2.11, which also comes from the interwinding capacitances. The sum of the second and third terms on the right-hand side of equation 2.12 is  [Yc]  =  WDIHYCA -[Y'D] + [YC2]  -[Y'DIHYCA  (2.13)  [^1 + ^C,]  There are only two distinct elements in the above matrix, i.e. the diagonal and off-diagonal blocks. By substitution, each of the matrices in the diagonal position of equation 2.13 becomes  33  -CH  -CHT  -CHL  CH + CHL+CHT  0  0  -CHT  -CT  -CTL  0 0  CT + CTL + CHT  0  0  C'L + CTL + CHL  -CHT -CTL  -CL  CHL  + CHT  -CHT  -CHT  -CHL  CHT + CTL  -CTL  —CHT  -CTL  CHL + CTL  and each of the off-diagonal matrices is -CH  -CHT  -CHL  -CHT -CT —CTL -CHT -CTL -CL  -CH  + -CHT -CHL  -CHT -CHL CT -CTL -CTL  -C'L  =  .  0 0 0 0 0 0 0 0 0  Substituting the diagonal matrices and the off-diagonal matrices into [Yc] makes it exactly the same as the matrix in the right-hand side of equation 2.11. Therefore, it can be concluded that the representation in fig. 2.10 is also the valid model for a transformer when interwinding capacitances are included. The impedances in delta of fig. 2.10a are transformed back to the T-circuit and then back to the original circuit of fig. 2.9a. The only difference now is that the leakage impedance of each winding has a part of the interwinding capacitance in it. This makes it consistent to synthesise each of them with a number of RLC blocks, each block for each resonance to be matched.  5. Implementation of Interwinding Capacitance Modification In the model, modification of interwinding capacitance connection is done in the nodal decoupled domain. The capacitance value in the decoupled mode is the same as that in the coupled mode because the interwinding capacitance exists only between 34  windings of the same phase. Only the part of the interwinding capacitance which is combined with the decoupled leakage impedance is modelled in the decoupled mode. All other capacitances resulting form the modification process will be additionally entered into equation 2.5 in the phase domain as other stray capacitances are treated. Therefore, there will be either four (two-winding) or nine (three-winding) more capacitances to be included in each phase. Eventually, equation 2.5 can be written symbolically as:  ^[Cstray\[Vnode(i)]  =  [Ihist-c(t - M)}  (2.14)  Equation 2.14 describes the model for stray capacitances in phase co-ordinates and, along with the model built from the frequency dependent part, (to be explained next in chapter 4) constitute the completefrequencydependent model for the transformer.  35  Chapter Three MODELLING OF FREQUENCY DEPENDENT COMPONENTS It was mentioned briefly in the previous chapter that the frequency dependent components in the proposed transformer model consist of the leakage impedance in combination with part of the interwinding capacitances. Also included in the measured parameters are the hidden effects arising from the simplified modelling of the actual stray capacitances with constant capacitances. In this chapter, the detailed modelling of the frequency dependent part will be discussed at length. The objective in the modelling of the frequency dependent part is to find a representative network consisting of a combination of constant parameter components, which are basically resistances, inductances and capacitances. The parameters will be chosen such that the network producesfrequencyresponses which closely match the actual data obtained from laboratory measurements.  1. Short-Circuit Responses of Transformers The measurements required to determine the characteristics of thefrequencydependent branch of the transformer model are short-circuit tests. These tests need to be performed with a device which is capable of generating a signal of variablefrequencyso that the measurements can be taken in a broadfrequencyrange, extending from a few Hz  36  up to a few MHz. Most transformers will exhibit their frequency dependent behaviour clearly in this frequency range. Figure 3.1 shows a frequency response of a distribution transformer under a short-circuit test. This type of response is typical for most power transformers [16, 17]. This measurement was taken at the high voltage laboratory of the Electricity Generating Authority of Thailand (EGAT). The measuring device available at EGAT was a network analyser model HP4192A which has a working frequency range between 5 Hz and 13 MHz. The four-terminal pair measuring technique (known as "Kelvin connection") was used in order to achieve a wide range of impedance measurements (1 mQ to 10 MQ) and to minimise measuring errors due to parasitic coupling with the test  i  i  Magnitude  (U  100 50  -50 •100 100  Ik  10k  100k  2 &o <o "O  I CO  m  JS  OH  1M  Frequency (Hz)  Fig. source  3.1  Typical short-circuit frequency responses of a transformer.  source"  37  Zmeasured  leads. The network analyser was connected to a personal computer (PC). Variation of frequency in discrete steps, either on a linear scale or a log scale, can be accomplished via a software installed in the personal computer. The network analyser measures the magnitude and phase angle of the device under test (DUT) at each discrete frequency and transmits the measured data to the PC to be stored. The data can then be retrieved for plotting or can be transferred to a diskette for later processing. The short-circuit test data so obtained are not yet in a form useful for transformer modelling. In order to obtain the data for the frequency dependent branch of the model, the value of stray capacitances (known from prior measurements and/or calculations) must be deducted from the raw data. Under short-circuit condition, the measured impedance is  where, Zshorti®)  =  Cext  = Sum of all capacitances except part of the interwinding capacitance which becomes combined with the leakage impedance. = Impedance of the frequency dependent series branch,  Zwindingio)  Measured short-circuit impedance,  and co Then, Z w w,ng(co)  = 2%i can b e c o m p u t e d a s  38  £>winthng\f&)  —  (3.2)  1  j(0Ce  Zshorti®)  The stray capacitances ( Q J that is deducted must not include the part of the interwinding capacitance (Cm/n in fig. 2.8b) which has been combined with the leakage impedance to make 2YtMs^  2. Equivalent Network Representation for ZvjMStlg The magnitude of ZwMtn^ is similar in shape to the magnitude of the measured short-circuit response. It possesses many portions with sharp peaks, these portions are not symmetrical around their peak and there are many peaks clustering together,  100k  11 i i I T 1 " " — i — r - f - m T  -r—r'TTiTin  10k Ik T3 3  100  •4-J  10 1 0.1  —>—I.I m i l  10  1—i—••*•••'  100  Ik  •  •  10k  i 4 ii m i  100k  <  •  i  1M  • i—i  10M  Frequency (Hz)  Fig. 3.2 Frequency dependent characteristics of the real part of the series impedance winding'  39  especially in the highfrequencyregion. Due to this nature, it is somewhat difficult to directly apply the pole-zero approximation method [10] to find a rational function representing Z^^g. However, it has been found that the real part of 2 ^ ^  reveals more  recognisable information. It displays more clearly the resonance characteristics which consist of several peaks superimposed on each other, similar to those shown in fig. 3.2. The peaks of the real parts are almost symmetrical around their resonance frequencies which are similar to the response of the real part of a parallel RLC network. This suggests a simple solution to the problem which otherwise could have been significantly hard to solve. As a first approximation, those peaks in the real part of Zw/ldnir to be duplicated are each fitted with a parallel RLC network. The resistance of this network can be viewed as representing ohmic losses in the windings and parasitic losses in the metallic parts caused by leakage fluxes. The inductance represents the leakage fluxes and the capacitance represents part of the interwinding capacitance. Therefore, all the circuit elements have associated physical meaning. The resonancefrequencyof the block is can be read off from the plot of ZwbuSng at the peak it duplicates. For an RLC parallel block, its real part is  Re{Z(a>)} =  ^ ^ R2(l-®2LC)2+a>2L2  (3.3)  where, Re{Z(&)  = Real part of the impedance of an R L C parallel block,  co  = angular frequency in radian/second,  40  L, R and C  = inductance, resistance and capacitance respectively of the RLC parallel block.  The product LC in equation 3.3 is known, it is the inverse of the square of the resonance frequency (angular frequency). Therefore, it is only necessary to estimate two more variables: R and L . If the real part of each resonance peak is expressed as in equation 3.3, the sum of all these functions will, approximately, represent the real part of the entire series branch. If there are n peaks, the expression for the desired function will be  <Q2L2Rt  Re{Zynn<ti„g{Gi)}  \  ;=1  ti  1-^r 2 CI )  (3.4)  2  + <o2L2  where, Re{Zwi„dmg(G))}  = Function representing the real part of the frequency dependent impedance (Zwmfi ) of the transformer model,  n  = no. of peaks to be matched,  Rt and Lt  = R and L of the block representing peak i, = angular frequency corresponding to the resonant frequency of peak /'.  In equation 3.4, there are 2 x n unknown variables, namely, Ri and Lt. The variable Q, is known, and the capacitance for each block can be calculated from Q, after the inductance Li has been found. The variables Rt and Lt  are estimated by fitting the function  ReiZwindingi®)} t 0  Re z  the real  P^  of  Ending • s i n c e  41  i winding^))  ls a  non-linear function, a  non-linear curvefittingmethod must be used tofindthe unknown variables. In this thesis, a non-linear least-squares fit with the Levenberg-Marquardt method [11] has been applied. Even though the effort is spent entirely onfittingthe real part of the function, the imaginary part of the function will automatically be matched. As explained in [1], the imaginary part of any analytical function is uniquely determined by its real part. The result from the first approximation technique described above will match the response of Z^^  moderately well except in the neighbourhood of the "major peak". The  major peak is the one with maximum amplitude, usually occurring as the first peak in the frequency response, and contains most of the information available from the short-circuit tests. Because only one RL block is used for the major peak in the first approximation, the fitting will not match this region properly and will need to befinetuned. The asymptotic pole-zero approximation method of [1,10] is used to find more RL blocks (series Foster) for the major peak region. Capacitance of this section (known from the first approximation) is first removed from thefrequencyresponse (as in equation 3.2) leaving a non-resonant response which can be fitted with simple poles and zeroes. The rational function used to approximate the impedance after extracting the C is  (s +p0)(s +Pl)(s +P3).(S  where, Ffirsti®)  -  +/?„)  rational function to approximate the remaining magnitude of the first peak after its capacitance has been removed, 42  K  = a constant,  s  = j®,  z. and p.  = zeroes and poles of the function,  m and n  = no. of zeroes and poles, respectively, m < n.  The poles and zeroes in equation 3.5 can be taken from the plot of the magnitude function  (/R2(G)) + © 2 Z 2 ( © ) )  on a log-log scale. Applying partial fractions to equation 3.5  yields  fi tK  "  (S+Pl)  (S+P2)  (S+P3)  or,  [>+t) ("^ L\iJ  (s+a \  L\i  Equation 3.7 gives all the desired parameters for the first peak. Using these RL blocks instead of a single RL block will improve the approximation a great deal especially in the region from dc to the first peak of the frequency response. However, due to the fact that parameters for the first peak and the remaining peaks are obtained at different stages in the optimisation process, the resulting combination will not be globally optimised. Further improvement is required to increase the quality of the fitting. The previous results from the first optimisation and the fine-tuning will be used as initial condition for another optimisation. This time, the model function is based on the real part of the circuit shown in fig 43  3.3. Resistances Rw, R^...., representing dc resistances for each peak are inserted because in reality each coil section will have a small dc resistance. Insertion of those resistances makes it possible to satisfy the dc response of the network. Examples of how  C, 'ii  A/W  C)W 11  R2o  R12  L31  '21  '12  ^TPh  R3O  R21  /m R31  Fig. 3.3 RLC synthesis network to approximate Z,'winding'  % winding1S realised from the test data will be given in chapter 6. Next, the parameters of the equivalent network for Z^^  must be discretized. This can be done in the normal,  straightforward manner.  3. Discretization of Equivalent Network for Zwjn£ng The equivalent network for ZyitnA consists of a number of constant R, L and C's. In spite of the fact that all the R, L and C's are constant, the overall characteristics of the network is frequency dependent and follows that of the impedance Z^.^  it represents. In  the EMTP, as well as in other electromagnetic transient programs based on time domain simulation, all circuit elements must be represented in a discrete time form. The discretetime representation of L and C depends on the integration rule used in performing the  44  discretiztion [10]. Table 3.1 summarises the discrete-time representations for a capacitance and an inductance with respect to trapezoidal and backward Euler rules.  Table 3.1—Summary of discrete-time models for inductances and capacitances using the trapezoidal and backward Euler integration rules. Circuit Element  Integration  Difference Equation  Rule Trapezoidal  Kt)  =  ^••v{t)+^-vit-At)  Keq  ~  At  IUt-At)  w«- At >  -e-  Kt) = f-v(t) + Kt-At)  Euler  R eq  Trapezoidal -v(t;  L_ At =  Kt)  =  p  -  M.  ~  2C  eq  •e "eq  Backward  Kt)  Euler  R  =  + Kt-At)  i(t-At) ?£•*)-?£•*-At)-i(t-*t)  Ihistit-At) W1-^  At ^-vit-AD  Backward  IuM-At) Kt)  =  + Kt-At)  =  2C v(t-Af)-i(t-At) At  £-[v(r)-v(f-A/)] AT  e  9  C  htstit-At)  =  -±.v(t-At)  Figure 3.4 shows one of the blocks used to model the frequency dependent branch. This block has to be transformed into its corresponding discrete time  45  c, (a)  rnfT\  /YT\  r^TPn  ^VW ^WV R,;  WW  R,  ct  U*) (b)  Ht)  ut) W W  WW  Rii  WW  R«  R«  vt(t)  Ut) (c)  Kt)  r ^ ^ , ^  ut>WW 4 W RLR-II  RLR-U  WW RLR-13  V4(t)  Ihist-l  (e)  (d) i(t)  Kt) -o  o-  eq-1  woLAAAd  v i(t )  KLR-I  Fig. 3.4 Discretization of a section of the equivalent network for Z^^ represents a general network used to duplicate a resonant peak. 46  . This section  representation. Each block of the approximating network has a single C regardless of how many RL parallel blocks the section contains. Parameters of the network in fig. 3.4a are found from non-linear fitting and/or from the pole-zero approximation described earlier in this chapter. This general network can also include the case of a section containing a dc resistance: an R can be regarded as an RL block without an L. The discretization process can be carried out as outlined in fig. 3.4. Initially, each inductance is replaced by its discrete time equivalent (fig. 3.4b). After that the equivalent resistance of the discretetime inductance is combined with the physical resistance of the RL block (fig. 3.4c), resulting in a combined circuit which has the same form as that of a single C or L (table 3.1). Subsequently, all these circuits can be added mathematically leading to afinalequivalent current source connected in parallel with its source resistance. At this point, the capacitance C, is replaced with its discretized form to produce the circuit of fig. 3.4d. After another combination, the network becomes again one dc current source with its internal resistance as shown infig.3.4e. The sequence described above is repeated for all blocks of the approximating network. After the process has been completed, all the equivalent circuits (which are connected in series) are again combined into a single equivalent similar to the one shown in fig. 3.4^ If the equivalent resistance of L, • is R^,. and its current source (or history term) is I, j(t-At), j = 1, n, where n is the number of RL parallel blocks, the following equations can be written (with reference to figs. 3.4b and 3.4c),  47  RLR-U  !t1J'Ri~lJ  =  , /JV*0,ZV*0,  (3.8)  Aly + ./U-l/  Rv  , Ly = 0,  /fc. v  (3.9)  , Rv = 0,  (3.10)  and, /L-V  =  h  , Lv*0,  (3.11)  0 , Lu = 0  (3.12)  and with reference to fig. 3.4-c,  vi(0  =  (iB-i(i) -IL-U)  • RLR-U + (iB-i(t)-IL-n)  • /?ZR-I2  + (Jw(d-/wj) •^-13 + • • • =  *B-I(0 • (RLR-U + RLR-U +RLR-U + •••) - (JL-U • RLR-U +IL-U • RLR-U +IL-I3 • RLR-U  =  h-\ (t) • S RLR-U -  H—)  S IL-U • RLR-U  or, n  /A  iB-iKt) -  -s  Vl(/)  RLR-  £ /t-i;' RLR-U + s  +/«-!  (3.13)  with reference to fig. 3.4d,  48  i(t)  = ic-i(t) + iB-i(i)  -T 4. V l ^ ) 4. A 4. ^ -•/C-l + — + lLR-i + — Rc-i  KLR-I  =vi(0 • ( ^ - + - ^ — J + /c-i + /L/M VKC-1  KLR-l S  = £&+/WlM(f_A/)  (3.14)  Keq-l  In deriving equation 3.13 and equation 3.14, the term (t-Af) was omitted from the history terms for simplicity. Equations 3.8 to 3.14 can be used to compute the equivalent resistance and the history term for peaks no. 2, 3, and so on. If Req and I^Jt-M) denote the equivalent resistance and history term of the discrete time representation of the entirefrequencydependent branch, the computations needed for Req and IMst(t-At) will be the same as those needed for the corresponding terms in equation 3.14, which are as follows:  Req  =  m TsReq-i  ,  /' = 1, JW  (315)  1=1  m 2-i ihist-i ' K-eq-i  hist(t-At)  =  &—-  ,  ts-ea  where, m  = no. of peaks to be fitted.  49  /=1,m  (3.16)  and, R,eq-i  Rc-i  Rm-i  The process described must be implemented three times - once for the zero-sequence network and repeated for the positive- and negative-sequence networks. Since the parameters of the positive-sequence network are identical to those of the negative-sequence network, R^ are the same for both networks. The history terms, though, must be calculated independently. The next and the last action to be taken in modelling Zvllndllv is to update the history terms of its discrete time representation so that the simulation can proceed to the next time step (t+At).  4. Updating History Terms To update the history terms for a capacitance or an inductance, the voltage across its terminals and the currents flowing through it must be known before any calculations can be performed. From the power network solution, the voltage at every node in the system in phase co-ordinates is known at each particular time step. Terminal voltages of the transformer are, therefore, known at the current time step in phase co-ordinates. These values have to be transformed into sequence co-ordinates in order to update the history currents of the model discrete-time circuits. The steps to update the history terms are as follows:  50  (1) Calculate vi(0 (fig. 3.4e), Vl(0  =  VW-hisMt-AW-Req-!  (3.17)  (2) Compute iC-i(i) and/ 5 -i(0 (fig. 3.4d), ic-xif) iB-i(t)  = Ic-i(t-At) + %@= lLR-i(t-At) + £@-  (3.18) (3.19)  (3) Evaluate the voltage of each RL block, v^iy(f) (fig. 3.4c), VRL-u(i) = [iB-iW-h-yit-AfyRLR-u  (3.20)  4) Find the current of the element Ly (fig. 3.4-b), ii-yit)  = iB-i(t)-^Q  (3.21)  Ky  (5) Determine the history term of each individual L and C at the current time step, h-y(t) and Ic-i (0, using the formulas given in table 3.1. (6) Combine the history terms of all inductances, n  /«-I©  2 h-y(t) • RLR-IJ = ^—p  (322)  where, n  = no. of RL parallel blocks in the section for peak no. 1. (7) Finally, add Ic-iO) and ILR-IQ) together to obtain history terms associated  with peak no. 1.  51  The history terms for the blocks associated with the other peaks are calculated in the same manner. In the end, all the history terms can be combined together to get the total equivalent history term which will be used along with R to formulate the nodal admittance matrix at the next time step t+At, m  £ 'hist-i\f) " Req-i  hiM  =  ^  Reg  (3.23)  where, m  = no. of peaks to be fitted, m la K-eq—i  ty-eq  The history term in equation 3.23 need to be computed three times for the three symmetrical component networks. Although the equivalent networks for ZYlitMni in the positive and negative sequences are identical, the transformer voltages and currents are different in all the three sequences. The discrete-time model for the frequency dependent part of the transformer described in this chapter will be used subsequently to construct the full transformer model. Details will be found in the next chapter.  52  Chapter Four THE COMPLETE MODEL All the component parts of the transformer model were discussed in the previous chapters. The equivalent circuits for these components were also developed for both the frequency domain and the discrete-time domain. In this chapter, all the pieces are put together to form the full transformer model. Equivalent networks for the leakage impedances and lumped terminal capacitances, together with ideal transformers, are the basic building blocks for the full transformer model. In the derivation of the full three-phase model, it will be assumed that the transformer is physically symmetrical so that symmetrical components can be applied. This eliminates the need to model directly the mutual coupling between phases, which may cause numerical stability problems [13]. The proposed approach not only does make the modelling simple but also reduces the number of elements to be mathematically synthesised. Although there are no transformers in existence which are truly symmetrical, tests performed on actual power transformers show that there are no significant errors in making such an assumption [1].  53  1. Discrete-Time Model in Sequence Domain Following the discretization procedure described in section 3 of chapter 3, the discrete-time model for a three-winding three-phase transformer in either of the three se-  Fig. 4.1 Discrete-time equivalent circuit for a sequence model of a transformer.  quences [14] (positive, negative, or zero sequence) is depicted in fig. 4.1. The transformer is visualised as three coupled windings [15]. This representation is also know as a T-circuit [21]. The magnetising branch is not shown in the model of fig. 4. As an approximation, this branch can be added externally to the model. The exact placement of the magnetising branch is not critical in switching transient studies in which the frequencies 54  are beyond the kHz. At these frequencies, the transformer core behaves close to ideal since the core flux is inversely proportional to the frequency and the incremental values for these high frequency harmonics will be very small. However, if the magnetising branch model at the power frequency is available, it can be incorporated to the proposed model in the same way as it is implemented in the conventional model. (The particular modelling of the magnetising branch is beyond the scope of this thesis). The models for the positive and negative sequence circuits will be exactly the same in the frequency domain but they will be different in the discrete time representation due to the history terms. Terminal capacitances are not shown in the network of fig. 4.1 because it is more convenient to model the capacitances directly in phase co-ordinates. Once the circuits of fig. 4.1 are transferred to the phase domain, the capacitance can be added to obtain the complete model. With reference to fig. 4.1, the following equations can be written:  ii(i)-Ri+ei(i)  =  h(t)R2+e2(t)  = v2(f) + ih2(t-Af)R2  (4.2)  h(t)R3+e3(t)  =  v3(t) + ih3(t-At)R3  (4.3)  + n3h(t)  (4.4)  niii(t)+n2i2(t)  vi(J) + /Ai(*-A/)-tfi  »2-«i(/)-«i-ea(0  55  =  =  0  0  (4.1)  (4.5)  m-ei(t)-m-e3(f)  (4.6)  = 0  where,  hit)  = current in winding /,  ihiit-At)  = history terms of frequency dependent part for winding  R,  = equivalent resistance of frequency dependent part for winding /',  nt  = no. of turns of winding /',  etit)  = voltage at the ideal transformer terminal on the side of winding /',  v,(0  = terminal voltage of winding i.  The terminal voltage vft) in equations 4.1 to 4.3 becomes known after the system of equations has been solved at the time step t. Therefore, there are only two unknowns, /\(f) and et(t), in equations 4.1 to 4.6. The unknown variables can be found systematically if equations 4.1 to 4.6 are rewritten as matrix equation as follows:  Ri 0 0 «1 0 0  0 R2 0 «2 0 0  0 0 R3 «3 0 0  1 0 0 1 0 0 0 0 n2 -Tii «3 o  0 0 1 0 0 -Ml  hit) hit) hit) exit) e2it) e>it)  1 0 0 0 0 0  0 1 0 0 0 0  or, in a more compact form,  56  0 0 1 0 0 0  Ri 0 0 0 R2 0 0 0 R3 0 0 0 0 0 0 0 0 0  vi(0 v2(f)  v3(0 ihiit-At) ih2it-At) ih3it-At)  (4.7)  POM  [R] U\ Wi] [N2]  U] [R] [*] [*]  Ht)]  MO] [ih(t-At))  (4.8)  where, Ri 0 0 0 i?2 0 0 0 R3  [R]  « i n 2 H3  0 0 0 0 0 0  Wi\  0  [N2]  0  H2 -Hi  m  0  0 0  -HI  [Km  = [ hit) hit) hit) l^,  MM  = [ *i(0 e2(0 e3(0 J^,  MO]  = [ V!0 V2(0 V3(0 J ,  Mt-At))  = [/Hi(r-AO ih2(t-At) ihiit-At) ] ,  [7] and [<D]  =  and unity matrix and zero matrix, respectively.  Equation 4.8 is also valid for the two winding transformer. In this case, submatrices [R], [NJ and [N2] will be of order 2 and their elements will not contain the third row and third column of the submatrices for the three winding transformer. Also, [i(t)], M0L 1X0]  an  ^ Vhit-At)] wi^ contain only the first and second elements. Expanding  equation 4.8 gives the following two equations: 57  [R][Kt)} + [e(i)] =  Ht)] + [R][ih(t-At)]  (4.9)  and, [Nx][Kt)] + [N2][e(t)] = [0]  (4.10)  Multiplying equation 4.9 by [N2] and subtracting equation 4.10 yields  [N2][R][Kt)]-[NrUKt)]  = [N2][v(t)] + [N2][R][ih(t-At)]  (4.11)  or, BOM =  m MO] + [h(t- At)]  (4.12)  where, W2][R]-m]\-i[N2],  m  =  [h(t-Aty\  = [Y\[R][W-Aty\.  Equation 4.12 relates the currents of each transformer winding to the terminal voltages and the history terms of all three windings for any of the three sequences. Matrix [Y] is the admittance matrix for the corresponding sequence and its order is equal to the number of transformer windings. Three equations similar to equation 4.12 will be required to construct the complete model in the sequence domain, for the zero, positive and negative sequences. The full model in symmetrical components representation is  58  [izero(t)] [iposit)] {inegii)]  [Yzer0] [0] [0] [0] [Ypos] [0] [0] [0] [Y„eg]  [Vzero(t)] [VposV)]  [v neg (0]  [ihzero{t-At)] [ihpos(t-ti)\ [ih„egif-kt)]  (4.13)  where, zero, pos and neg  signify the variables for zero, positive and negative sequences, respectively, and = zero matrix (all elements are zeroes).  [0]  If ['«,(0L [Yseq], [ v „,(0] and {ih^t-At)] denote the currents, the admittance matrix, the terminal voltages and the history terms in equation 4.13, then the transformer model in the sequence domain can be expressed as  [iseq(t)]  =  [Yseq] [vseq(t)] + [ihseq(t - At)]  (4.14)  2. Discrete-Time Model in Phase Co-ordinates Equation 4.14 is not yet ready to be used in transient simulation programs, such as the EMTP, in which all the voltages and currents must be expressed in phase co-ordinates. Equation 4.14 can be transformed into phase co-ordinates by means of a transformation matrix. For a balanced systems, the transformation matrix is not unique. Due to the fact that the symmetrical components transformation involves complex numbers, other transformations with real matrix elements, such as the one given by Edith Clarke [12], will be used for the numerical transformation procedure. Clarke's transformation matrix for a three-phase system is  59  1  1  1  73 72 76 [r 3 x 3 ]  1  =  1  -1  (4.15)  73 72 76 1  o  . V3  -2  76 _  and,  i-i  [^3x3]"  =  (4.16)  [7^3x3]  Matrix [r 3x3 ] given in equation 4.15 can be modified to fit the transformer problem by replacing each of its scalar elements with a diagonal submatrix. The order of the submatrix is determined by the number of windings of the transformer, i.e., two for a two-winding transformer and three for a three-winding transformer. All the diagonal elements in these submatrices are identical and equal to the corresponding scalar value of the original transformation matrix. For example, the transformation matrix for a three-phase, three-winding transformer is the following [9x9] matrix  m =  -±r 0 0 -± 0 „ „0 j . 0 0 7? 72 J6 0 - ^ 0 0 - ^ 0 0 -pv 0 71 72 76 JJ J2 J6 -^ 0 0 - ^ 0 0 - ^ 0 0 JJ Jf -1 J6 0 ~=r 0 0 0 0 1 0 72 76 0 0 1 0 0 -1 0 0 1 73" 72 76 - ^ 0 0 0 0 0 ^ 0 0 7376 1 0 -7=- 0 0 0 0 0 4 0 0  0 -=W 71 0  0 60  0  0  76 0 ^pr 7?  (4.17)  [7] -1  = [T\T  (4.18)  The voltages and currents of each transformer winding in phase co-ordinates are related to the corresponding sequence voltages and currents by [14]  [vpka,.®] = [7][v^(0]  (4.19)  and, [iphaseii)]  = [T\[iseq(t)]  (4.20)  where, [v/*fl«(0] ^ d [iphase(t)]  = voltages and currents in phase co-ordinates,  [v (f)] and [^eg(01  =  voltages and currents in sequence co-ordinates.  Applying equations 4.19 and 4.20 to equation 4.14 yields  [iphaseW] = [T\[Yseq]Urx{vphase{t)\ + [T\[ihseq(t-Ai)]  (4.21)  or, [iphaseif)] = [Yphase][Vphase(t)]  + [ikphaseif - At)]  where,  [YPhase]  = mt^im-1,  [ihphase(t-At)]  = [T][ihseq(t-At)].  and  In full matrix form, equation 4.22 for a three-winding transformer is  61  (4.22)  Vll(f)  /11(f) [Ym]  [iy  /12(f)  [Ym]  vi 2 (f)  in(t)  v i3 (f) V2l(f)  /21(f) 122(f)  =  [Ys]  [Ym]  [Ym]  V22(f)  v23(f) V3l(0  /'23(f) '31(f) /32(f) . /'33(f) .  [Ym]  [Ym]  [Ys]  -  -  v32(f) . v33(f) _  ///n(f-Af) Ihn(t-M) Ihn(t-M) Ih2\(t-At) + Ih22(t-At) Ihnit-At)  (4.23)  7/j3i(f-Af)  Ihn(t-At) Ihn(t-At)  where, i«(t)  current of phase /', windings,  vyit)  terminal voltage of phase /', winding j ,  [Ys] and [Y„.]  =  [3x3] submatrices (or [2x2] for a two-windii  former) for the self and mutual terms of the admittance matrix [Yphase] in equation 4.22, element of history vector associated with phase i, winding/  Ihijit-At)  In the actual implementation of the model, there are only two transformations involved. Phase voltages need to be transformed into sequence voltages in order to calculate and update the history terms. After that, the history terms are transformed from sequence co-ordinates back into phase co-ordinates and are entered into the right hand side of the system equations. There is no need to calculate [Yphase ] in equation 4.23 directly from [ritl^Jt?]" 1 because it is already known that [Ys]  =  J • ([Yzero] + 2[Ypos])  and, [Ym]  Also, from equation 4.19, [vMg(f)]  =  =  \-([Yzero]  -[Ypos])  U\ 1[vPhase(t)], or in a full matrix form 62  1  I  0  0  1  0  0  1  0  I  0  0  0  0  V2l(0  0  0  0  0  V22W  -1  0  0  0  V23W  -2  0  0  V3l(0  -2  0  V32(0  -2  V33W  0  0  0  1  0  0  1  0  0  i  0  0  I  0  0  1  0  =  Vpos-sit) V„eg-l(t)  72  0  0  1  0  0  1  0  0  1  0  76  V„eg-2(t)  0  v weg - 3 (0  0  76  0  J2  76  76  Vl3(0  0  I  Vpos-2(t)  1  0  0  -l  V12C)  1  0  Vpos-l(t)  0  0  V  73-  1  0  0  zero-l{j)  Vll(0  1  Vzero-lif)  73"  0  0  75-  73"  0  0  Vzero-\if)  73-  0 0 -l  72"  76  0  73"  72  76  7T  76  73"  76  0  73"  76  (4.24)  If the elements of the sequence voltage vector [vseq(t)] in equation 4.24 are arranged in the order of windings, instead of sequences, equation 4.24 becomes  Vzero-\(t)  Vpos-lit)  1  1  1  71  73  73"  1  -1  J2  0  0  0  0  0  0 " " Vn(f)  0  0  0  0  0  0  0  V2l(0  0  0  0  0  0  0  V31W  0  0  0  vnif)  vneg-i(t)  1  1  -2  76  76  76  Vzero-lif)  0  0  0  0  0  0  Vneg-2(t)  0  0  0  Vzero-lif)  0  0  Vpos-zij)  0  Vneg-lif)  0  Vpos-lit)  =  1  1  1  73-  73-  1  -1  72-  72  0  0  0  0  V22W  0  0  0  V320  1  1  -2  76  76  76  0  0  0  0  0  0  0  0  0  0  0  0  0  0  63  1  1  1  73  J3  75-  Vl 3 (0  1  -1  72  72  0  V23W  -2  v 3 3 (0  1  1  76  76  76" _  (4.25)  where, Vzero-i(f), Vpos-tii), V„ eg _,(0 :  Zero-, positive- and negative-sequence voltages for winding /' respectively.  Equation 4.25 can be further simplified, if the voltages are rearranged into matrices instead of vectors. Now, the sequence voltages can be computed with a single matrix multiplication, which is more efficient, as follows:  1  ( 0 V'zero-2{t)  1  1  Vzero-3(t)  yr yy yj  Vpos-3(t)  -J- ^ -  0  _J  -2  Vpos-l(t)  VpoS-2(fy  rneg-1 ( 0  Vneg-2(t) v neg _ 3 (0  JZ  1  J6 J6  vn(/) vi 2 (0 vn(t) V2l(t) v22(t) v 23 (0 v3i(f) v 32 (0 v 33 (0  (4.26)  The sequence voltages in equation 4.26 can be reshuffled into a vector as in equation 4.24, which is in the desired order in sequence co-ordinates. The sequence currents can be calculated likewise. For the history terms, transformation from sequence variables to phase variables are also required. A similar implementation to that used to derive equation 4.26 can be applied to obtain the history terms in phase co-ordinates. But the transformation matrix [r3x3] is used in this case.  3. Nodal Admittance Matrix Equation 4.22, however, is not yet of the right form for an electromagnetic transient program solution because the variables are associated with voltages and currents in the transformer windings. It must be further modified so that voltages and currents are  64  Winding 2 Winding 1  Winding 3 o 3  1 o-  Phase 1 10 o-  -o 12  4 o-  -o 6 Phase 2  13 o-  -° 15  7 »-  -o 9 Phase 3  16 o-  -o 18  Fig. 4.2 Transformer node labelling for external network connection.  the node voltages and the currents flowing into the node of the transformer. This involves labelling of the nodes in a certain way, as shown in fig. 4.2. Since the terminal voltage is the difference of the voltages of the two nodes at both ends of the winding, it follows that  Vll(/)  =  Vl(/)-Vio(fl  (4.27)  V12O)  -  v 2 (f)-vn(f)  (4.28)  Vl3(0  =  V3(*)-V,2(*)  (4.29)  65  V2l(t)  =  V 4 (/)-Vi 3 (0  (4.30)  V220  =  V 5 (/)-V 1 4 (0  (4.31)  V 23 (0  =  V 6 (/)-Vi 3 (d  (4.32)  V3l(0  =  V 7 (0~Vl6(0  (4.33)  V32O)  =  v 8 (r) -  vn(i)  (4.34)  V33<0  =  V 9 (0 - Vig(f)  (4.35)  where, = node voltage with respect to fig. 4.2.  vi(0  In addition, it is also true that /10(0, 'u(0> i\$)>—> 'i$)> a r e  ec ua to  l l  the negative of i^t),  /2(f), /3(0,--, '9(0> respectively. Arranging the node voltages and currents in this manner makes it possible to construct the final equation in phase co-ordinates without much work. It is just an extension of equation 4.22, which is as follows,  [tnodeO)]  I phase ~ I phase  J- phase  [v„ode(t)] +  *phase  ihphaseit-At) -ihphaseif - At)  (4.36)  where,  ['«*(01  vector of currents flowing into node 1, node 2, node 3,..., node 18,  66  [vno(fe(f)]  = vector of voltages of node 1, node 2, node3,..., node 18,  [^*aJi  =  [ihphase(t - At)]  = history terms in equation 4.22.  admittance matrix in equation 4.22,  Equation 4.36 is the nodal equation without terminal capacitances. The complete frequency dependent model for the three-phase transformer comprises equations 4.36 and 2.14 combined. Elements of the admittance matrices from both models can be added algebraically to the corresponding elements of the system admittance matrix. The history terms in equation 4.36, though, must be subtracted from the right hand side of the system nodal equation because the history terms are considered known and must, therefore, be on the opposite side of the equation.  67  Chapter Five MEASUREMENT OF THE MODEL PARAMETERS  This chapter describes the techniques used to measure the transformer parameters and to claculate the circiut parameters from these measurements. As explained earlier, there are two types of parameters in the model - stray capacitances, represented by constant capacitances, and the leakage impedances represented by branches of constant RLC sections. Numerical examples on the network synthesis realization of the model parameters will be presented in chapter 6.  1. Measurement of the Frequency Dependent Leakage Impedance The measurements required to obtain the leakage impedances are the short-circuit frequency response tests. Measurements must be conducted such that both the zerosequence and the positive-sequence impedances can be determined. The measuring device used in the present research work for impedance measurements was an "Impedance Analyzer", which produces a low voltage signal with variablefrequencies.The signal voltage energizes the transformer under test. At the same time the Analyzer senses the current that flows into the transformer. Voltages and currents are then used by the device to calculate the magnitude and phase angle of the transformer impedance within a band of 68  frequencies specified by the users. The results are both displayed and stored in a diskette. The "single-phase" Impedance Analyzer can be used to measure both the zero-sequence and the positive sequence impedances. Details of the measurement technique are explained next, with special note on transformers with delta-connected windings.  1.1 Transformer without Delta Windings The zero-sequence tests can be performed as usual. However, since the Impedance Analyzer is a single phase source, the positive-sequence tests cannot be conducted directly. In this case, the positive-sequence information will be calculated from the zero-sequence data and some additional data obtained from separated measurements. A three-phase, three-winding transformer (without terminal capacitances) can be described by the following matrix equation [14]:  IA-I IA-2  [Ys]  [Ym]  [Ym]  IA-3 IB-X IB-2  [Ym]  [Ys]  [Ym]  VA-\ VA-2 VA-s VB-i VB-2  IBS  VB-3  IC-I  Vc-x Vc-2 Vc-3  Ic-2 Ic-3  [Ym]  [Ym]  [Ys]  (5.1)  where. IA-J, IB-J  and IC-j  transformer currents in winding j of phases A, B and C respectively, at a particular frequency.  69  VA-J, VB-J and VQ-J  = terminal voltages of winding j of phases A, B and C respectively, at a particular frequency.  [Ys] and [Ym]  = [3x3] self and mutual submatrices of the transformer admittance matrix. Both [Ys] and [Ym] are symmetric.  From the zero-sequence short-circuit test, a zero-sequence admittance matrix is obtained. This matrix is related to [FJ and (YJ by  [Yzero] =  [Ys]+2[Ym]  (5.2)  If either [Ys ] or [Ym] is available, it can also be used to compute the positive-sequence admittance matrix. Suppose that [Ys ] is available from the measurement. The positivesequence admittance matrix can then be calculated form the zero-sequence admittance as  [Ypos] = ±0[Ys]-[Yzero])  (5.3)  In fact, it is not necessary to know all the elements of the matrix [Ys]. Only the diagonal elements are sufficient for the calculation of [Y^]. The diagonal elements of [Ys ] can be measured easily with the impedance analyzer. The readings will be taken at one coil while all the other coils are short-circuited at their terminals. The measured admittance can be expressed, for example for phase A, as:  jV*  =  !±L  f  vA_iMj = 0, VB-I = 0, Fc_, = 0; / = 1,3.  *A-j  where, 70  (5.4)  IA-j and VA-j  = current and voltage measured at the terminals of coil j , 7=1,3.  Measurement of ys_M should be taken for phases A, B and C and the averaged values of the three measurements should be used. The diagonal elements of [Y ] can be computed now with equation 5.3. If the actual measurements could be performed, these values would correspond to the positive-sequence impedances obtained from the short-circuit tests measured on one winding while the other two windings are shorted. If all stray capacitances are not taken into account ( modelled seperately), the positive-sequence short circuit admittances can be related to the winding impedances as follows (circuit representation similar to fig. 2.9b) :  1 YH  TA*>  T L  = z +  '^  ZT+ZL  - *T + z77FL - ¥T  (5 6)  '  where, YH, YT and YL  = diagonal elements of the positive-sequence admittance matrix when measurements are taken on the high, low and tertiary voltage windings, respectively. The windings where measurements are not taken are all shorted,  Z'TmdZ'L  =  (^)2Zrand(V)2ZL  71  ZH  -  positive-sequence frequency dependent impedance for high-, tertiary- and low-voltage windings, respectively,  nH, w^and nL  = no. of turns for high-, tertiary- and low-voltage windings, respectively.  Since Ym YT and YL can be taken from the corresponding diagonal elements of [Y^], ZH, Z\ and Z'L , or, alternatively, ZH+Z'T , Z'T+Z'L and ZH+ Z'L can be found by solving equations 5.5, 5.6 and 5.7. Rearranging terms in equation 5.5 gives  ZT + ZL  =  YH  • (ZH • ZL + ZH • ZT + ZL • ZT)  (5.8)  Similarly, from equation 5.6 and equation 5.7,  ZH + Z'L  =  Y'T-<ZH-Z'L+ZH-Z,T  Z'T + ZH  =  YiiZH-Zi  + Z,L-Z'T)  + ZH-Z'r + Zi.Z'r)  (5.9)  (5.10)  If ZJJZ'T +ZHZ'L + Z'JZ'J is denoted by Zn, then Equations 5.8, 5.9 and 5.10 can be solved for ZH, Z'T and Z'L, which are given by  ZH  = \(-Y„+rT+rL)-Zn  (5.11)  Z'T = \(YH-fT+rL)-Zn  (5.12)  Z'L = \{YH+Y,r-Y'L)-ZTl  (5.13)  72  However, Zn is still unknown.  From direct multiplication, the following equations  follows:  ZHZfT = \Zl\-Y2H-Y2T+Y2L+2YH-YT)  Z'TZ[  ZHZ!L  (5.14)  = \ZI.{Y2H-Y2T-Y2L+2YL.YT)  (5.15)  = \ZI-(-Y2H  (5.16)  + Y2T-Y2L+2YH-YL)  Adding equations 5.14, 5.15 and 5.16 together yields Z% = \Zl-[-Y2H-Y2T-Yl+2(YH-YT  + YT-YL + YH-YL)~\  (5.17)  Dividing equation 5.17 by Zn * 0 and rearranging terms results in  Z„  =  —5 4 2 2 -Y\ -Y T-Y L + HJH-YT + YT• YL + YHYL)  (5.18)  In conclusion, the positive-sequence impedances for the transformer without delta windings can be obtained in the following manner: 1) Find the zero-sequence impedance in the usual way from the short-circuit tests. 2) Measure the diagonal elements of the admittance matrix in phase coordinates for the three phases.  73  3) Average the corresponding measurements from the three phases to obtain the diagonal element of the self-admittance sub-matrix. 4) Calculate the diagonal elements of the positive-sequence admittance matrix with equation 5.3. 5) Compute Zw with equation 5.18. 6) Substitute Zn back into Equations 5.8, 5.9 and 5.10 to obtain the short-circuit impedances for the positive sequence. In this case, the short-circuit impedances are those which would be obtained if the short-circuit tests were conducted between pairs of windings at a time. Equations 5.11,5.12 and 5.13 can also be used to calculate the impedances of the individual windings of the "T" circuit. If this was done, however, one or more of the individual impedances might not be minimum-phase-shift. This problem can be avoided if the combined impedance of every two windings isfittedinstead. By applying the above outlined steps, it is only necessary to measure nine functions for a three-winding transformer and six functions for a two-winding transformer in order to get the positive-sequence data for a transformer without delta windings. The total number of measurements to obtain both the positive and zero-sequence data for the three and two-winding transformers comes, therefore, to a total of twelve and eight, respectively. This indicates a dramatic reduction in the number of tests which would normally amount to ninety-one for a three-winding transformer and to thirty-six for a two-winding transformer if all the elements of the admittance matrix were to be  74  Winding 1  Winding 2  Fig. 5.1 Measurement of positive-sequence impedance in a Yd transformer.  measured. This latter case is, for instance, where the transformer is modeled directly from the [Y] matrix.  1.2 Transformer with Delta Winding If one of the windings of a three-phase transformer is connected in delta and all the delta terminals are accessible from the outside, the positive-sequence test can be conducted with a short-cut. It is very common for a power transformer to have one of the windings connected in delta, but the delta winding might have just one accessible terminal to be grounded. If all the delta winding terminals are accessible, measurements can be taken on one winding while the other windings are either shorted or open.  75  1.2.1 Two-Winding Transformer A two-winding transformer (without terminal capacitances) can be described by the following matrix equation:  " IA-X '  )>sXX  )>sX2  ymXX  )>mX2  )>mXX  )>mX2  VA-x  IA-2  )>s2X )>s22 ym2\  J V 2 2 ^ J 2 1 ,V»I22  VA-2  h-l  ymu ymxi y*w ysYl ym\\ ym\2  VB-X  IB-2  y&x ym22 y m y$22 yax yntn  VB-2  ymXX ymX2  Vc-x Vc-2  Ic-x . !c-x .  ymlX  ymX2  ySXX  ysX2  y m ym22 yS2x yna% yS2\ y*22  (5.19)  where, IA-i, Is-i and Ic-  Currents entering terminals i, i = 1,2 of phase A, phase B and phase C, respectively,  VA-U VB-I and Vc-t  = Terminal voltages of terminals i, i = 1,2 of phase A, phase B and phase C, respectively,  Assume a short-circuit test is performed on this transformer and that the measurement is taken on winding 1, which is delta-connected, while winding 2 is shorted. Under this condition, VA.2 = VB_2 = Vc_2 = 0. Consequently, the measured currents in winding 1 are  IA-X  ysu ymu ymu  VA-X  IB-I  ymXX ysll  ymXX  VB-X  Ic-x  ymXX ymXX  ysXX  Vc-x  76  (5.20)  If VA_, , Vg^ and Vc., are in positive sequence, i.e., having equal magnitude and 120 degrees apart in phase angles, the ratio of IA_, to VA_, would give the positive-sequence admittance of the transformer. This could be accomplished if a variablefrequency,threephase source were available. However, a single-phase source with variable frequency such as an impedance analyzer can be used also to obtain the same result. Suppose phase C of winding 1 is shorted to ground, as shown in the arrangement offig.5.1. Because phase C is shorted, VCA is zero. VBA is connected to the source in its nornal polarity while VAA is in the reverse polarity. We then have that VAA = -Vs, VBA = +VS, IAA = (-ysU+ymU) and IBA = (+ysU-ymu)- The measured current Isflowingfromthe source is  h  = IB-I — IA-I =  2-(ysn-ymn)-Vs  or, |r  = 2-(ym-ymii)  (5.21)  Equation 5.21 is, in fact, an admittance which is twice as much the short-circuit admittance of a two-winding transformer in the positive sequence mode. The model's series impedance is then twice the impedance measured in the way described. The measured impedance contains both the leakage impedance and stray capacitances. The stray capacitances can be obtained from short-circuit measurements in the high frequency spectrum where they become the dominant factor. What remains after removing these capacitances  77  can be combined with part of the interwinding capacitance to become the frequency dependent series impedance of the model, as explained in chapter 4.  1.2.2 Three-Winding Transformer The measurement technique explained in the previous section for the two-winding transformer can be applied to the three-winding transformer as well, but the measurements will involve two or three winding at a time. For example, if a transformer is connected in YYd (high-voltage, low-voltage and tertiary-windings are connected in wye, wye, and delta, respectively), the zero-sequence measurements can be performed on the wye-connected windings while the other two windings are either both shorted or only one of them is shorted. Two tests can be conducted to measure the short-circuit impedances between each of the two wye-connected windings with the delta winding shorted while the other wye-connected winding is left open. The third measurement will be taken on either of the wye-connected winding while the two other windings are shorted. A positive-sequence measurement can be taken on the delta winding by first shorting one of the wye-connected windings. A second measurement may then be conducted with the previously open winding short-circuited and the previously shortcircuited winding open. The third measurement should be taken on the delta winding while the other two wye-connected windings are short-circuited.  78  From these measurements, the leakage impedances for each individual winding or the sum of two of them can be computed as follows (assuming that the tertiary winding is connected in delta and the equivalent circuit is similar to that of fig. 2.9b):  1.2.2.1 Zero-Sequence Tests There are a few possible combinations of ways to perform short-circuit tests for the three-winding transformer. The following combination involves two complete shortcircuit tests and one partial short-circuit test: (a) Perform a short-circuit test by shorting the high-voltage and low-voltage windings to ground. The impedance ZA seen from the low-voltage side (after the terminal capacitances have been removed) is:  or,  ©* = ^frf  522  < >  (b) Perfor short-circuit test with high-voltage winding open. The measured impedance ZB on the low-voltage side (after removing the terminal capacitances) is  * = fe)Vr + 2t) or, 79  {jjffzB  = Z'T + Z'L  (5.23)  (c) Perform a short-circuit test with both tertiary and low-voltage windings shorted. The measured impedance Zc on the high-voltage side (after removing terminal capacitances) is  Zc = ^ + 5 - 4  (5-24)  All the symbols used in equations 5.22, 5.23 and 5.24 have the same meanings as those in equations 5.5, 5.6 and 5.7. Equations 5.22, 5.23 and 5.24 can be solved for Z'T as follows: (1) Multiply equation 5.22 and 5.24 by (ZH+Z'T) and (Z'L+Z'T) respectively, ZTZH  + ZTZL + ZHZL  = (ZH + ZT) ZA • {jfrj  and ZJZH  + ZTZi + ZHZL  — (Zf+Zi) Zc  which, together with equation 5.23, results in  (ZH + Z'T)ZA-(%f)2  = (Z'T + Z'L)ZC =  {j^)2zBZc  or ZH+Z'T  =  % ^ ZA  80  (5.25)  (2) Subtract equation 5.22fromequation 5.23,  &)'<*-*> = z - S -•&  <»*  (3) Solve equation 5.25 with equation 5.26,  Z  T = ±{w)rZc(z7B~ZA) Z  ( 5 - 2? >  A  Then, ZHand Z'L can be solved by substituting Z'T into equations 5.23 and 5.25 which are,  Z  £  =  intY2*-2*  (5 28)  -  and, ZH  = ^—Zff  (5.29)  1.2.2.2 Positive-Sequence Tests For the positive-sequence tests, all measurements are taken on the delta winding. The tests may be carried out as follows: (a) Perform a the short-circuit test with the low-voltage winding shorted. The measured impedance (after the terminal capacitances have been removed) is:  = 81  ZT + Z'L  (5.30)  (b) Perform a short-circuit test with both high-voltage and low-voltage windings shorted. The measured impedance (after removing the terminal capacitances) is  \nT)  \JIH)  = ZT +  ZT +  zH+(%)2zL  * • & ) ' & * + ( & ) * * •  ZT +  ZT +  7' .7'  (5.31)  Z'H + Z'L  (c) Perform a short-circuit test with only the high-voltage windings shorted. The measured impedance (after removing the terminal capacitances) is  * - &)T(&r*+z„  = zT + z'H  (5.32)  Equations 5.30, 5.31 and 5.32 can be solved simultaneously for ZT . By subtracting equation 5.30 from equation 5.3, and equation 5.32 from equation 5.31, the following two equations are obtained:  82  Zb - Za  -(zj)2 z'„ + z'L  =  (5.33)  and  Zh — Zc —  (5.34) Z'H+Z'L  Multiplying equation 5.33 with equation 5.34 one obtains:  {Zb - Za)(Zb - Zc)  7'7'  -  "12  z'H + z'L\  (5.35)  From equation 5.31, it follows that  7'7' Z'H + Z'L]  (zb - zTy  (5.36)  or  (Zb-ZT)  =  ± J(Zb - Za)(Zb - Zc)  (5.37)  which gives  ZT  -  Zb ± JiZb - Za)(Zb - Zc)  (5.38)  The leakage impedances for both zero and positive-sequences for the threewinding transformer are next converted into delta, as shown in fig 2.9c. The interwinding  83  capacitances are then added (fig. 2.10a) so that in the network synthesis procedure they from part of the RLC equivalent networks.  2. Measurement of Stray Capacitances Some of the stray capacitances of the model can be directly measured but some of them must be calculated from the short-circuit test data. The method on how these stray capacitances are obtained is described next.  2.1 Direct Measurement Stray capacitances which can be directly measured with a capacitance measuring device are the capacitance-to-ground and the interwinding capacitances. The device used at the high voltage lab of EGAT is the "Capacitance and D.F Bridge", model CB100. The device has three connection leads marked CL, CH and G which should be connected to the high-voltage, low-voltage windings of the transformer and ground respectively to ensure correct measurements. When the device is adjusted until a balanced condition is achieved, the values of capacitances can be read off the meter dial. Fig. 5.2 depicts a connection of the device to the transformer under test for a standard measurement of capacitances between the high-voltage and low voltage windings (C^) and capacitances between the high-voltage winding and ground (CHG ) of a two-winding, three-phase transformer. In the measurements, all six terminals of the high voltage winding are connected together (terminal H in fig. 5.2) and so are all terminals of the low-voltage winding (terminal L in  84  APPLICATION Transformer Under Test L C Standard  Measurement taken is CH-L + CH-G. CB100 internally grounds CL in mis configuration, thus effectively shorting out CL-G.  Capacitance & D.F. Bridge (Model CB100)  Fig. 5.2 Application of the Capacitance & D.F. Bridge to measure stray capacitances in a transformer. The configuration shown is for the measurement of Cm plus CHG.  fig. 5.2). Two more standard measurements are required to get C^ and CHG . From these measurements, the capacitances to ground and the interwinding capacitances of the transformer can be found. For a three winding transformer, the interwinding and windingto-ground capacitances present in the transformer are as shown in fig 2.3, chapter 2. In this case, the measureing device CB100 is connected to a pair of the windings at a time while the third winding is grounded. An illustrative arrangement for measuring stray capacitances associated with the high-voltage and low-voltage windings of a three-winding transformer is shown infig.5.3.  85  C Standard  JZJLX  -r-  / s/ /  s /s /  s N N  f- -t  \  / s s / ss / ssss / f- H*s f"-I  Transformer under Test Capacitance & D.F. Bridge (Model CB100)  Fig. 5.3 Measurment of the sum of interwinding and winding-to-ground capacitances in a three-winding, three-phase transformer with CB100.  2.2 Stray Capacitances Calculated from Short-Circuit Tests Stray capacitances for the model which are not available from direct measurements can be calculated from the short-circuit frequency responses. At the high end of the frequency range, the behavior of the short-circuit impedance becomes close to that of a capacitance. In this frequency range, all the inductive elements  have very high im-  pedances and can be considered open-circuited. The stray capacitances that can be obtained from the short-circuit tests data are the turn-to-turn and the phase-to-phase capacitances. The turn-to-turn capacitance are always present in the short-circuit tests but the phase-to-phase capacitance will be present only when the measurements are taken from the outermost winding on each transformer limb.  86  2.2.1 Turn-To-Turn Capacitances In the proposed transformer model, the turn-to-turn capacitance is lumped as a single element and is connected across the two terminals of each transformer coil. Firstly, suppose that the winding where the short-circuit test is conducted is not the outermost winding so that the phase-to-phase capacitance is not present (the case in which the measurement is made on the outermost winding will be addressed next in the Phase-to-Phase Capacitance section). In practice, the zero-sequence measurement is taken on the wye-connected winding because it is not possible to conduct the test on the delta winding with a single power source. In making the measurement, all three terminals of the wye are joined together. As a result, the measured turn-to-turn capacitance is three times that of a single coil. At the high end of the frequency response, a value of capacitance can be calculated from the magnitude of the impedance. By subtracting the known values of winding-to-ground capacitances and interwinding capacitances, which were previously obtained from the direct measurements, the value of the turn-to-turn capacitance of the winding under test can be determined. The turn-to-turn capacitance of a delta winding, on the other hand, can be obtained from a positive-sequence test conducted on the delta winding, as explained earlier. The turn-to-turn capacitance obtained in the measurement is two times that of a single winding because two coils are connected in parallel while the third coil is shorted. With  87  simple subtraction, as in the case of the wye winding, the turn-to-turn capacitance for the model can be determined. For a three winding transformer, the measurements should be taken while the other two windings are simultaneously short-circuited. This simplifies the calculation of the turn-to-turn capacitance because all apacitances in the model that are not shorted become connected in parallel.  2.2.2 Phase-To-Phase Capacitance The phase-to-phase capacitance which exists between the outermost windings of the transformer legs are schematically shown in fig. 2.5. Figure 2.5a shows the case in which the windings are not connected, fig. 2.5b corresponds to the case in which the windings are connected in delta, and fig. 2.5c corresponds to the case in which the windings are connected in wye. Phase-to-phase capacitances of full value are present across each coil in the delta configuration, while only halves are effectively present between the wye terminals. In the positive-sequence tests (which are conducted on the delta windings), the phase-to-phase capacitances are combined with the turn-to-turn capacitances. For the delta connection, the phase-to-phase capacitance has to be modeled together with the turn-to-turn capacitance with both connected across the winding terminals, as shown in fig. 2.5b, because there is not enough information from the tests to separate them. For the wye connection, on the contrary, the phase-to-phase capacitance will not be present in the zero-sequence test even though the winding is the outermost one. This has the  88  advantage that it makes it possible to isolate the turn-to-turn capacitances form the test results. With reference to fig. 2.5c, if a short-circuit test is carried out with all the other coils short-circuited, the phase-to-phase capacitance will be the only unknown capacitance present in the test. The phase-to-phase capacitance can then be isolated and modeled as in fig. 2.5a.  89  Chapter Six NUMERICAL EXAMPLES  The proposed method to calculate all the parameters required for the proposed frequency dependent transformer model was presented in chapter 5. This chapter is aimed at putting the method into practice. It covers some examples on how the model parameters are realised from the information gathered from the various short-circuit tests. The frequency responses calculated from the model are compared to those obtained from the tests. In addition, time-domain simulations with the model using the estimated parameters are verified by comparison with the results recorded in the laboratory. Examples are given for both two- and three-winding transformers.  1. Two-Winding Transformer The numerical data for a two-winding transformer is based on measurements taken on a power transformer rated 30/40/50 MVA. The transformer is a two-winding transformer with the high-voltage winding connected in delta and the low-voltage winding connected in wye. The voltage rating is 115A/23Y kV. The transformer was manufactured by Asia Brown Boveri (ABB).  90  1.1 Capacitances The capacitances which can be directly measured with the Capacitance & D. F. Bridge for the tested transformer are listed in table 6.1.  Table 6.1—Measured values of capacitances for the tested transformer Type of capacitance  Capacitance for 3 phases (pF)  Winding-to-ground -High-voltage (CHG)  3,418  -Low-voltage (CLG)  12,395  Winding-to-winding 6,441  -High-voltage to low-voltage (C^)  To reassure that stray capacitances can be assumed to be constant within a specific range of frequencies, a frequency response test was carried out by shorting all the high-voltage terminals (delta) to ground and connecting the three terminals of the low-voltage terminals (wye) to the wye neutral. The measurement was taken between the low-voltage terminals and ground from which the sum of C^ and CLG was obtained. The result of this test is shown in fig. 6.1 together with the equation Z = 8.10x106 f0994 which resulted from fitting the straight line portion of the capacitance curve to a function Z = afb. As expected, the power of the frequency (f) is close to -1.0 because the 91  Frequency Response of C-HL + C-LG 1,000,000  1,000,000 6 . 0 994 Z = 8.10xl0f  E  o  100,000  100,000  10,000  10,000  1,000  1,000  T3 3  to  100  100  i N 11  10 10  100  1,000  10,000  100,000  I  1,000,000  10 10,000,000  I I I I I 11  Frequency (Hz)  Fig. 6.1  Transformer capacitance impedance measured with the impedance analyser  magnitude of the impedance of a capacitance is inversely proportional to the frequency. The coefficient of f gives an estimate of (27CC)'1 and C is found to be 19.644e-09 farad. The sum of C^ and CLG taken from table 6.1 is 18.836e-09 farad which is relatively close to C. Therefore, the values of the stray capacitances taken from the two different measurements agree with each other. Other capacitances which must be calculated from the short-circuitfrequencyresponses will be given later.  1.1.1 Capacitances Determined from Zero-Sequence Test A zero-sequence short-circuit test was carried out with the circuit configuration shown in fig. 6.2a. Only one test is sufficient to get the zero-sequence data of a two-winding transformer. The test result was fitted with the network of fig. 3.3. The graphs in fig. 6.4 show the measured and fitted responses. The circuit parameters for the equivalent network are given in table 6.2. Numbers in table 6.2, and in other similar tables, are shown to five significance digits since sufficient accuracy in these values is needed for the correct placement of poles and zeroes of thefittingfunction. For the tested transformer, the capacitances which are present in the zero-sequence test are the turn-toturn capacitance of the low-voltage winding, the interwinding capacitance, and the capacitance of the low-voltage winding to ground. The purpose offittingthe test results is to find out how much capacitance is included in the measured impedance. The measured short-circuit impedance was fitted up to thefrequencyat which the data are believed to  93  x  X  (b)  2  l  H  N  0  / X  !  c»—-&— ^H3"  3J  _nnnru\  H  2  tank  Fig. 6.2 Schematic of short-circuitfrequencyresponse tests: (a) Zero-sequence test; (b) Positive-sequence test. 94  10,000  10,000  Estimated 1,000  1,000  Measured 100  100  10  10  0.1  0.1  0.01  0.01  0.001  0.001  J  0.0001  10  I  I  M I N I  100  1,000  I I II I  10,000  i  i i n  100,000  i  i  i  i i i Mi  1,000,000  i  0.0001 10,000,000  i i i  Frequency (Hz)  Fig. 6.3 Real part of short-circuit impedance (zero-sequence) for the two-winding transformer.  180 Estimated 135 Measured Phase angl& 90  -/V~~  45  in <0 <D  feb 0> w CD  "bb  c  -45  -90  -135  J  10  I  M I N I  100  I  1  M  1,000  I  I  I  J  I I  I  10,000  I Mill  100,000  1  L  1,000,000  Frequency (Hz) Fig. 6.4 Short-circuit impedance (zero-sequence) for the two-winding transformer.  -180 10,000,000  a  Table 6.2— Parameters of the zero-sequence short-circuit equivalent network for the tested two-winding transformer Peak No.  1  2  Capacitance  RL parallel block  (farad)  Resistance (ohms)  2.96201e-08  4.69634e-05  9.54346e-06  1.02500e+04  6.52899e-05  2.21436e-01  6.46047e-06  1.66978e+00  6.36733e-05  1.61492e+00  3.80926e-05  7.79935e+02  1.63820e-04  1.18106e-02 7.54854e+00  3  1.59905e-07  3.64525e-08  3.90087e-08  7.71330e-08  3.52545e-08  4.72049e-06  3.86791e-05 2.47888e+01  7  1.65258e-05  9.79614e-05 7.46147e+01  6  8.13311e-06  2.00860e-05 2.22654e+02  5  7.97444e-06  2.10021e-07 1.94427e+01  4  Inductance (henry)  1.38830e-06  4.31114e-05 5.09657e+01  1.97948e-06  High-frequency equivalent capacitance = 7.43398e-09 farad  97  be reliable, which is about 1 MHz. This capacitance is subtracted from the short-circuit impedance to obtain the leakage impedance of the transformer. This impedance is later combined with part of the interwinding capacitance (^- CHL) to form thefrequencyde2n pendent series branch of the transformer model (Z^^ ). The combined impedance is then fitted once more to obtain the RLC network offig.3.3. From the previous procedure, the turn-to-turn capacitance of the low-voltage winding can be calculated as follows: Total capacitancefromtable 6.2 iCL-o  = = =  7.434e-09 3.220e-09 4.132e-09  farad farad farad  Turn-to-turn capacitance  =  0.082e-09  farad  ;*-CH-L  The factor 1/3 is used instead of 1/2 for C ^ is because the low-voltage winding was grounded while the measurement was carried out. When the winding is grounded, the effective winding-to-ground capacitance of the winding is reduced to 1/3, as suggested by Allan Greenwood [9].  1.1.2 Capacitance Determined from Positive-Sequence Tests Three tests were conducted to determine the positive-sequence short-circuit impedance of the transformer. Each test is performed with one of the coils of the delta winding shorted, as shown in fig. 2b. Results of the tests are shown in fig. 6.5 in which all the three tests are superimposed. It can be seen in fig. 6.5 that there is no significant difference between results of the three tests. The average of the three tests is used as input data 98  100,000 ;  180 HI -H2 shorted =  — 135  HI -H3 shorted 10,000 -  \— Phase angle  i  H2-H3 shorted  \^jiM \ 1,000 -  fa \ I V\  /  90  I  *  45  GO  <U <D s-.  tx>  D  -a c  0>  100 jr JT  -  :  \ 1  x  '/ ' /  \ \ I s  \  S  Magnitude  /  k, j  IA  I  !  -45  * W c¥ I I I ell  \ i IV <v  /  \\\j  \ -90  i  10 -135  i  10  i  i  i i i i i 1  100  i  i  i  i i i i i 1  1,000  i  i  i  i i i i i 1  i  i  10,000 Frequency (Hz)  i  i i i i i 1  100,000  i  i  i  i i i i  11  1,000,000  Fig. 6.5 Short-circuit impedances (positive-sequence) of the two-winding transformer measured on different windings.  111  -180 10,000,000  </> c« PH  100,000  100,000 Estimated  10,000  10,000  Measured  1,000  1,000  100  100  10  10  0.1  -  0.01  i  10  100  i  i  i  1 1 1 1  1,000  J  iii  10,000  100,000  1,000,000  0.1  0.01 10,000,000  Frequency (Hz) Fig. 6.6 Real part of the short-circuit impedance (positive-sequence) for the two-winding transformer.  100,000  180 Estimated 135 Measured  10,000 90  45  1,000  0 100  -45  -90 10 -135  i  10  i  i  M  100  i  i  i  i  1111  1,000  i  i  i  i  1111  j  10,000  i  i  i  1111  100,000  1,000,000  -180 10,000,000  Frequency (Hz)  Fig. 6.7 Transformer short-circuit impedance (positive-sequence) for the two-winding transformer.  Table 6.3~Parameters of the positive-sequence short-circuit equivalent network for the tested two-winding transformer Peak No.  1  RL parallel block  Capacitance (farad)  Resistance (ohms)  2.34919e-08  5.21162e-01 1.74892e+03  2  1.36025e-08  6.27161e-09  6  2.73908e-03  0.0 7.71918e+02  5  2.94064e-02  0.0  1.91069e-08  3.63068e+03 4  2.70085e-02  0.0 2.17093e+04  3  Inductance (henry)  1.48156e-07  1.06217e-04  0.0 5.72532e+02  1.83790e-04  0.0  -  7.35160e+01  1.17626e-05  1.27424e-07  High frequency equivalent capacitance = 2.91999e-09 farad  to find the circuit parameters. The same process implemented for the zero-sequence test described earlier is applied to the positive-sequence test to obtain the capacitances to be subtracted from the short-circuit branch. In this case these capacitances are the phase-tophase capacitances and the turn-to-turn capacitance of the high-voltage winding. The  102  result of the RLC synthesis procedure are shown in fig. 6.6 and fig. 6.7, while the circuit parameters are listed in table 6.3. Since, for the positive-sequence test, only two phases are involved, the capacitance values shown in table 6.1 must be multiplied by a factor of two thirds. The calculation of the unknown capacitances is as follows: Total capacitance from table 6.3  =  2.920e-09  farad  ;kfcH.L)  =  2.147e-09  farad  J(|CH-G)  =  0.760e-09  farad  Turn-to-turn plus phase-to-phase capacitances  =  0.013e-09  farad  0.020e-09 (3 phases)  1.2 The Series Branch Impedance  (Z^^  Each short-circuit impedance is composed of two components of the transformer model, i.e., the leakage impedance and the short-circuit capacitance. The short-circuit capacitance can be obtained from the test result with the procedure described earlier in sections 1.1.1 and 1.1.2. The short-circuit capacitance will be subtracted from the measured impedance to obtain the leakage impedance. The leakage impedance is combined with part of the interwinding capacitance which is transferred to the side of the leakage impedance to become the series branch impedance used in the model. For the zero-sequence, this capacitance is  Ul  xo.5  x6441  x 10 9  ~  = 5.367xl0-9 farad/phase.  103  1,000,000  1,000,000 Estimated  100,000  100,000 Measured  10,000  10,000  1,000  1,000  100  100  10  10  0.1  0.1  J  0.01 10  I 1 I I IIII  100  I  I  III  1,000  J  I I I I II  i  10,000  i  i  i 1111  100,000  i  i l l  j  1,000,000  i  i  0.01 10,000,000 i i 111  Frequency (Hz) Fig. 6.8 Fitting of the real part of the Zv/ind- impedance for the zero-sequence circuit of the two-winding transformer.  1,000,000  180 Estimated 135  100,000  Measured Phase angle 90  10,000 45  CO  <U <D  & CD  -a  1,000  0  c«J <u  -45 100 -90 10 -135  i  10  i  i i i 111  100  i  i  i  i i i iii  1,000  i  i  i  i i 1111  i  i  10,000 Frequency (Hz)  i  i M 111  100,000  i  M I  1,000,000  -180 10,000,000  Fig. 6.9 Fitting of the magnitude and phase angle of the frequency dependent Z winding impedance for the zero-sequence circuit of the two-winding transformer.  a) aJ J3 PL,  Table 6.4~Parameters of the zero-sequence ZytinAm impedance for the tested two-winding transformer (per phase values referred to high-voltage side) Peak No.  1  2  RL parallel block  Capacitance (farad)  Resistance (ohms)  4.76059e-10  9.02003e-03  1.27529e-07  2.81494e+04  8.64570e-03  3.44792e+01  1.91425e-03  2.3675 le+02  6.28649e-03  5.59907e+05  1.00467e-02  7.78905e-01 6.12763e+02  3  2.92333e-09  1.00064e-09  1.34472e-09  3.60143e-09  9.50437e-09  2.67993e-09  3.18966e-05  2.80912e-03 2.57767e+02  8  1.57439e-04  7.55318e-04 6.08895e+02  7  7.35136e-04  1.50116e-04 2.56361e+03  6  4.54154e-04  1.04837e-03 1.00262e+04  5  5.96037e-04  7.50465e-02 1.30156e+03  4  Inductance (henry)  8.18850e-06  3.21350e-03 5.25896e+02  2.98555e-05  High-frequency equivalent capacitance = 2.02039e-10 farad 106  And for the positive-sequence, the capacitance is  After the calculated amount of capacitances have been added to the leakage impedances, they are fitted with the synthesis RLC network. Figures 6.8 and 6.9 show the fitting results for the zero-sequence winding impedance. It is seen from these results that the calculated winding impedance of the tested transformer can be matched very well to 1MHz. The circuit parameters that resulted from the fitting are given in table 6.4. Five sections are used to simulate the response for the first peak and a total of eight peaks are matched all together. The parameters given in table 6.4 were converted to per phase values and transferred to the high-voltage side so that they can be used in the transients program directly. The positive-sequence winding impedance is treated the same way as the zero-sequence winding impedance. The fitting results for the real part are shown in fig. 6.10. Five sections were used to fit thefivedominant peaks of the real part and two RL parallel blocks were used to refine thefittingfor the principal peak which occurred in the vicinity of 10 kHz. The calculated results for the magnitude and phase angle of the positive-sequence Z ^ ^ are plotted in fig. 6.11. The fitted and measured results agree very well from 50 Hz to about 0.5 MHz. Even though the fitting results over 0.5 MHz are not as good as in the case of the zero-sequence impedance, the overall outcome can be considered sufficiently good. The circuit parameters for the network producing the response  107  100,000  100,000 Estimated  10,000  10,000  Measured  1,000  1,000  100  100  10  10  0.1  0.1  J  0.01 10  I  100  L  I  1,000  I  I  I I I  I  I  10,000  I  I  I  I  I I I I  100,000  1,000,000  I  0.01 10,000,000 I  I  I I I  Frequency (Hz) Fig. 6.10 Fitting of the real part of the Z^ding impedance for the positive sequence circuitsequence circuit of the two winding transformer.  100,000  180 Estimated 135 Measured Phase angle  10,000  90  45  CO  S-H  1,000 C D  -45  100  -90  -  10  j  10  i  i  i  i i  100  i  i  1 1 1 1  1,000  i  i  i  i  i i  i  i  10,000 Frequency (Hz)  i  i  i  1 1 1 1  100,000  i  i  i  i  1 1 1 1  1,000,000  -135  -180 10,000,000  Fig. 6.11 Fitting of the magnitude and phase angle of the frequency dependent Z winding impedance for the positive-sequence circuit of the two-winding transformer.  a  Table 6.5~Parameters of the positive-sequence Z^^ impedance for the tested two-winding transformer (per phase values referred to high-voltage side) Peak No.  1  2  RL parallel block  Capacitance (farad)  Resistance (ohms)  3.83689e-09  9.66480e-01 2.14788e+02  1.58297e-02  7.79506e+04  7.16340e-02  6.63480e-08  0.0 2.08566e+03  3  4.24785e-09  1.38289e-08  1.03526e-02  0.0 4.09054e+03  5  1.15472e-02  0.0 1.50852e+04  4  Inductance (henry)  8.29235e-10  1.90603e-03  0.0 1.01916e+03  4.51872e-04  High-frequency equivalent capacitance = 5.58860e-10 farad  shown in fig. 6.11 are given in table 6.5. The numbers in the table are the per phase values which are ready to be used in the transformer model software.  1.3 Time-Domain Simulation An impulse response test was also performed on the 50 MVA, 115/23 kV transformer from which the two-winding model data are derived. A standard 1.2/50 us voltage  110  (a) H  (b) 1,000  1,000 -1  -  i  1  500  -__.- 500 1A  -  e  I'C  k-l  u  /  Estimated -- -500  -500 Measured  -  -1,000  i  20  1  ,  1  40 60 Time (micro sec.)  ,  i  80  -1,000 100  Fig. 6.12 Impulse response test, (a) Circuit diagram; (b) Simulation result. Ill  impulse with a magnitude of 550 kV was applied to the delta winding while the wye winding was shorted. The schematic of the impulse test is shown in fig. 6.12a. The current in the 1 Q-resistor was recorded. Figure 6.12b shows the measured current response of the tested transformer against the model response using the parameters given in tables 6.4 and 6.5. The agreement between measured and computed results is reasonably good.  2. Three-Winding Transformer The three-winding transformer which was available for short-circuit tests was rated 25 MVA and 115/22/11 kV. The tertiary winding was connected in delta and the high and low-voltage windings were connected in wye. The transformer was manufactured by "Volta-Werke". The high-voltage winding is the outermost winding and the tertiary winding is located next to the core.  Zero and positive-sequence tests were  conducted on this transformer as explained in chapter 5.  2.1 Stray Capacitances Similar to the two-winding transformer case, the stray capacitances that can be directly measured with the capacitance meter are the winding-to-winding capacitances and the capacitances from windings to ground. The terminals of the windings belonging to the same phase are connected together, therefore, the capacitance values obtained are the sum of the three phases. The measured capacitances are as shown in table 6.6 below.  112  Table 6.6~Measured values of capacitances for the three-winding transformer Type of capacitance  Capacitance for 3 phases (pF)  Winding-to-ground -High-voltage (CHG)  531.5  -Low-voltage (CLG)  394.4 3.7  -Tertiary-voltage (CTG) Winding-to-winding -High-voltage to low-voltage (C^)  5,400.0  -High-voltage to tertiary-voltage (Cm)  87.5  -Low-voltage to tertiary-voltage (CLT)  7,900.0  In performing the short-circuit tests, one of the terminals of each winding is grounded. Therefore, half of the winding-to-winding capacitance and one-third of the stray capacitance to ground will appear in the tests. The capacitances that correspond to the short-circuit condition are shown in fig. 6.13. The turn-to-turn capacitances included infig.6.13 can be calculated from the measured short-circuit impedances. Since the tertiary-winding, which is connected in delta, is not the outermost winding, the phase-to-phase capacitances will not be present in the positive-sequence tests. Therefore, the phase-to-phase capacitances cannot be calculated from the positive-sequence results and the identification of those capacitances requires some additional tests. Due to the fact that these tests were not carried out for the tested transformer, calculation  113  —c + c  c  + C  3 ^ LG *" L  3 ^ HG ^ H  Fig. 6.13 Effective stray capacitances of a three-winding transformer during the short-circuit tests.  of the phase-to-phase capacitance will not be demonstrated here (an explanation on how to obtain the phase-to-phase capacitance is given in section 2.2.2 of chapter 5). The turn-to-turn capacitances CH, CL and CT in fig. 6.13 for the high-voltage, low-voltage and tertiary-voltage windings of the three-winding transformer can be calculated from the short-circuit tests. Six tests were conducted as follows: Zero-sequence ZoL HT ZoH LT: ZoL T:  impedance measured on the low-voltage winding with both high-voltage and tertiary-voltage windings short-circuited. impedance measured on the high-voltage winding with both low-voltage and tertiary-voltage windings short-circuited. impedance measured on the low-voltage winding with only the tertiary-voltage winding short-circuited. 114  10,000  135 Estimated Phase angle  Measured  90  1,000  45 <D <D  100  & <D  CD  Ho C  a  10  CD  -45  -90  0.1  I  I I I I I I  10  L  J  100  i  i  i 111ii  1,000  i  10,000  100,000  i  i  i  i 1111  1,000,000  i  -135 10,000,000 i 1111  Frequency (Hz) Fig. 6.14 Fitting of the short-circuit impedance measured on the low-voltage side with the high and tertiary-voltage windings short-circuited (zero-sequence).  a*  1,000,000  135 Estimated Phase angle  100,000  Measured  90  10,000  45 in CD CD CD w  1,000  CD  '60  c  100  -45  10  -90  J  i  i i 11 i n  10  m  100  i  i i i Mii  i  i i M i  1,000 10,000 Frequency (Hz)  i  i  i  i 11111  100,000  i  i  . i n  i  1,000,000  i  -135 10,000,000  Fig. 6.15 Fitting of the short-circuit impedance measured on the high-voltage side with the low and tertiary-voltage windings short-circuited (zero-sequence).  10,000  135 Estimated Measured  90  1,000 Phase angle  45 <L> <D  <D w 0)  100  c <u a) cS PL,  -45 10 -90  Mil  10  I  1 1 I I I I ll  100  I  1,000  [  I  10,000  I  I  I I I I I II  100,000  I  I  1,000,000  I  -135 10,000,000  L  Frequency (Hz) Fig. 6.16 Fitting of the short-circuit impedance measured on the low-voltage side with only the tertiary-voltage winding short-circuited (zero-sequence).  10,000  135 Estimated Phase angle  1,000  Measured  90  100  45 in  <U <D >-C  0)  T3 0>  10  c as 0) 1/3  aS  -45  -  0.1  iL  i  0.01 10  100  1,000  10,000  100,000  1,000,000  -90  -135 10,000,000  Frequency (Hz) Fig. 6.17 Fitting of the short-circuit impedance measured on the tertiary-voltage side with the high and low-voltage windings short-circuited (positive-sequence).  PL,  10,000  135 Estimated Phase angle  1,000  Measured  90  100  45 <U CU  10  c as CO  si  -45  0.1  0.01  -90  j  i  i i i  i  10  100  i  I  i i 111ii  1,000  10,000  I  III  100,000  L  I I I  I  1,000,000  I  I  -135 10,000,000 I I I I I  Frequency (Hz) Fig. 6.18 Fitting of the short-circuit impedance measured on the tertiary-voltage side with only the high-voltage winding short-circuited (positive-sequence).  1,000  135  Phase angle  90  100  45 10  <U  feb <U  T3  <D  "Hb  c a  <u CO  a  -45  0.1 -90 Estimated 0.01  J  i  i i 11111  10  100  J_ 1,000  Measured I  I I M I I  10,000  I  I  I  I II  I I I  100,000  I  I  I  I I I I I I  L  1,000,000  -135 10,000,000  Frequency (Hz) Fig. 6.19 Fitting of the short-circuit impedance measured on the tetiary-voltage side with only the low-voltage winding short-circuited (positive-sequence).  PL,  Table 6.7 — Parameters of the zero-sequence short-circuit equivalent network for the tested three-winding transformer (per phase values) Case: ZoL_HT  Case: ZoH_LT C (farad)  Case: ZoLT  Peak no.  C (farad)  K.(Q)  R,(Q)  L,(G)  Peak no.  1  7.0913e-9  6.7400e-5  9.0443e-l  2.4790e-3  1  2.8359e+l  5.8812e-4  1.2652e+3 4.4174e-2  7.1257e+0 5.1101e-4  4.2490e+4  5.5459e-4  2.4716e+5  1.6921e-l  2.1473e+4 7.6706e-4  **(Q)  R,(G)  2.2592e-9 2.1720e+0 1.7950e+l  L, (henry)  Peak no.  C (farad)  R*(n)  R,(Q)  L, (henry)  8.7832e-2  1  3.5442e-9  2.9224e-l  3.8500e+l  3.8599e-4  2  8.4926e-8  1.3492e-l  3.7518e+2  2.9200e-4  2  9.6098e-8  1.4946e+0 4.6004e+2 8.1353e-4  2  9.7392e-7  3  2.6750e-8  1.8318e-4 8.6748e+2  2.685 le-4  3  9.5892e-8 0.0000e+0 3.1787e+2 1.9087 e-4  3  1.0361e-7 2.1972e+0 4.7329e+2 2.2034 e-4  4  1.3775e-8  3.5086e-2  2.1U6e-4  4  2.2672e-8  1.8869 e-2 7.5099 e+2 2.4447 e-4  4  2.8894 e-8 5.4405e-l  5  8.3690e-9  4.2042e-2 4.7543e+2 7.5667e-5  5  7.6848 e-9 1.2863e+0 1.0750 e+3 2.6704 e-4  5  1.7176e-8  1.9244 e-3 2.2344e +3 2.0885e-4  6  1.0580e-8  2.4801e-4 5.8732e+2 9.0962e-5  8.0248e+2  Higji-frequency equivalent capacitance = 2.61602e-9 farad  High-frequency equivalent capacitance = 1.56810e-9 farad  5.8695e-l  2.0642e+2 9.6540e-4  5.3471e+2 1.8971 e-4  High-frequency equivalent capacitance = 2.08256e-9 farad  Table 6.8 — Parameters of the positive-sequence short-circuit equivalent network for the tested three-winding transformer (per phase values) Case: ZpT_HL  Case: ZpTH  Peak no.  C (farad)  R*(")  R,(Q)  L,(Q)  Peak no.  C (farad)  R*(")  1  5.8024e-9  2.3667e-3  3.4367e-2  3.7188e-5  1  3.2324e-8  4.8274e-3  4.4709e+0  5.7350e-5  2.7302e+0  8.4852e-5  Case: ZpT_L Peak no.  C (farad)  »*(Q>  R,(Q)  L, (henry)  1  6.7837e-8  1.0602e-2  1.0058e-l  6.8052e-6  1.2728e+4 5.1797e-4  4.2487e-K)  1.0914e-4  8.7523e-2  6.3950e-5  2.9290e-2  3.5993e-5  3.3721e-H) 4.6749e-6  3.2372e+l  8.0960e-4  1.0878e+l  2.7954e-4  1.2007e+3  1.0783e+3 6.5569e-4  6.6783e+3  5.5477e-4  1.5906e-4  R,(G)  L, (henry)  1.4251e+4 2.2247e-13  2  1.2704e-7  4.0226e-3  1.4424e+2 2.3170e-5  2  2.7388e-8  1.0324e-2  1.2219e+3 5.0313e-4  2  6.1203e-8  3  8.0675e-9  5.4044e-3  4.1702e+2  2.0476e-5  3  7.1427e-8  2.2936e-5 1.773 le+2 8.2175e-5  3  1.2609e -7 3.0937e-3  4  1.2732e-8  5.4202e-3  1.8349e+2  8.4803e-6  4  6.3884e-9  7.4689e-3  1.3138e-4  4  7.9259e-9  3.4538e-3 4.2665e+2 1.0958 e-4  5  7.1573e-8  3.1479e-3  5.7049e+l  8.7374e-7  5  5.0629 e-9  1.1541e-2 5.9023e+2 3.1310e-5  5  5.7841e-9  3.2588 e-3 2.6653e +2 2.4716e-5  6  4.8523e-9  5.3440e-3  2.0570e+2  8.0285e-6  6  1.0616e-8  1.4553e-2 3.0960e+2  1.0052e-5  6  1.651 le-8  2.9264e-3  1.2482e+2 3.1629e-6  7  1.2366e-8  8.2188e-6  2.8295e+2 4.3123e-6  7  5.3568e-9  3.2541e-3  1.7151e+2 6.7650e-6  8  3.0125e-9  8.8793e-7  2.4023e+2  High-frequency equivalent capacitance = 1.65901e-9 farad  5.6647e+2  3.0982e-3 5.5532e+2 2.2484e-4 8.0822e+l  4.5678e-5  U952e-5  High-frequency equivalent capacitance = 1.06100e-9 farad  High-frequency equivalent capacitance = 1.70849e-9 farad  Positive-sequence ZpTHL:  impedance measured on the tertiary-voltage winding with both high-voltage and low-voltage windings short-circuited. impedance measured on the tertiary-voltage winding with only high-voltage winding short-circuited. impedance measured on the tertiary-voltage winding with only low-voltage winding short-circuited.  ZpTH: ZpT_L:  To determine the capacitances included in the measured short-circuit impedances, the impedances are fitted by an RLC synthesis network. The results of this fitting for the six measured cases are shown in figs. 6.14 to 6.19. The parameters of the RLC networks are indicated in tables 6.7 and 6.8. All quantities in fig. 6.14 through fig. 6.19 and in tables 6.7 and 6.8 has been converted to per phase values. To determine the turn-to-turn capacitances from the short-circuit tests, it preferable to use those tests in which the remaining tow windings are both shorted. These correspond to cases :ZoL_HT, Z o H L T and ZpTHL. Using the total capacitances given in figs. 6.20, 6.21 and 6.24, together with the measured capacitances from table 6.6 and the capacitance diagram of fig. 6.13, CH, CL and CT are found to be: CH (from case:ZoH_LT) =  1.56810e-09 -  J(JCHL  + JCHT + }C HG )  = 5.94461 e-10 farad (per phase) CL(fromcase:ZoL_HT)  =  2.61602e-09-±(jC  = 3.5553le-10 farad (per phase). CT (from case:ZpT_HL) =  1.65901e-09 - j ( j C H T + jC L T + }C TG  = 3.27349e-10 farad (per phase).  123  300,000  135 Estimated  100,000  Measured  90  30,000  45 Phase angle  10,000  3,000  -45  1,000  -90  300  ii  10  100  i  i  i  i 111M  1,000  i  i 11111  10,000  i  i  i  i i M i  100,000  I I I I II  L  1,000,000  -135 10,000,000 _LLL  Frequency (Hz) Fig. 6.20  Fitting of the sum of the high and low-voltage winding impedances (zero-sequence).  135 100,000 Estimated Measured  90  10,000 Phase angle  45  1,000  -45 Magnitude  100  -90  I MINI  10 10  100  I  I  I I II  1,000  10,000  I I III  100,000  i  i i 111 u  i  1,000,000  i  -135 10,000,000  i i 1111  Frequency (Hz) Fig. 6.21 Fitting of the sum of the high and tertiary-voltage winding impedances (zero-sequence).  100,000  135 Estimated  30,000  Measured  90  10,000 Phase angle  45  3,000  1,000 -45 300 Magnitude  -90  100  30  LU  10  I  I  I  I I I I 11  100  I  I  I  I I 111  1,000  l_  i  10,000  100,000  i  i  i  1,000,000  j  -135 10,000,000  Frequency (Hz) Fig. 6.22 Fitting of the sum of the low and tertiary-voltage winding impedances (zero-sequence).  1,000,000  135 Estimated Phase angle  Measured  90  100,000  45 10,000  1,000 -45  100 -90  I I I I III  10  10  I  I  I  M i l l  100  1,000  l  I I III  10,000  I  I  I I I I II  100,000  1,000,000  -135 10,000,000  Frequency (Hz)  Fig. 6.23  Fitting of the sum of the high and low-voltage winding impedances (positive-sequence).  1,000,000  135 Estimated Phase angle  100,000  Measured  90  10,000  45  1,000  0  100  -45  10  -90  I  I  I  I  I I 111  10  i  11111  100  i  i  i  i  i i n  1,000  i  i  i  i  i  1111  10,000  100,000  i  i  1111 [  i  1,000,000  i  i  -135 10,000,000 i  1111  Frequency (Hz)  Fig. 6.24  Fitting of the sum of the high and tertiary-voltage winding impedances (positive-sequence).  1,000,000  135 Estimated Phase angle  100,000  Measured  10,000  90  45  1,000  100  -45  10  -90  i  i  i 11111  10  i  i  i  i  i 11  100  i  i  i 111 n  1,000  i i  i  10,000  i  i  i  i 111 n  100,000  i  i  i  i i  1,000,000  -135 10,000,000  Frequency (Hz)  Fig. 6.25 Fitting of the sum of the low and tertiary-voltage winding impedances (positive-sequence).  Table 6.9 ~ Parameters of the zero-sequence Z , ^ ^ impedance for the tested three-winding transformer (per phase values referred to high-voltage side) Case: Zo_HL Peak no. 1  C (farad)  R*W  5.5878e-10 5.2597e+2  Case: Zo_HT R,(Q)  L,(")  1.7928e+2 2.2965e-l  Peak no. 1  C (farad)  R*(n)  Case: Zo_LT R,(Q)  7.3880e-lO 1.0235e+2 6.8885e+l  L, (henry)  Peak no.  C (farad)  8.6127e-l  1  3.1390e-8  R*(")  R,(«)  L,(henry)  2.1268e+0 9.3466e+l  5.2827e-3 6.1356e-3  3.8902e+3  9.4942e-2  1.9761e+2 3.7134e-2  1.1302e+l  3.0746e+7  2.1162e-l  4.3893e+3  5.3395e-2  8.005 le+1 5.4825e-3  -  -  1.4622e+7  1.4181e-l  8.4193e+4 2.8075e-2  2.5810e+4  1.4326e-l  2  9.7171e-9  5.1643e-2 2.1409e+3 3.1125e-4  3  1.4427e-9 6.6190e+l  2  2.6048e-9  -  5.8748e+4  3.5960e-l  2  6.8965e-9  3  1.1945e-9  -  2.0818e+4 2.0708e-2  3  2.5250e-8  4  3.4651e-9  -  1.3832e+4 3.1728e-3  4  3.6330e-9  3.0161e+0 7.6270e+3  1.1210e-3  4  6.9580e-9  -  2.3692e+4  5  4.8565e-10  -  2.9141e+4  1.5505e-2  5  3.0130e-9  4.4105e+l  4.8946e+3  4.3305e-4  5  3.9650e-10  -  4.5959e+4 1.8990e-2  6  9.1151e-10  -  2.0995e+4  4.5676e-3  -  -  -  -  -  6  1.4772e-9  -  5.1015e+4  2.8185e-3  7  4.7763e-10  -  4.9790e+3  2.3018e-4  -  -  -  -  -  7  8.9075e-10  -  4.2739e+3  1.2343e-4  High-frequency equivalent capacitance = 1.16971e-10 farad  2.4032e-2  High-frequency equivalent capacitance = 4.66133e-10 farad  2.8105e+l  1.4079e+3  1.2851e-2  1.7237e+4 1.7146e-2 1.580U-3  High-frequency equivalent capacitance = 1.88918e-10 farad  Table 6.10 — Parameters of the positive-sequence Z ^ ^ impedance for the tested three-winding transformer (per phase values referred to high-voltage side) Case : Zp_HL Peak no. 1  C (farad)  R*(S)  Case: Zp_HT R,(Q)  2.3803e-10 7.3622e+0 5.5878e+l  Ii(")  Peak no.  1.6394e-2  1  C (farad)  R*(n)  Case: Zp_LT R,(Q)  L, (henry)  7.1241e-10 3.5828e+0 7.0207e+0 6.4685e-3  Peak no. 1  C (farad)  RdcW  R,(Q)  L, (henry)  1.0879e-9 2.3638e+0 5.5462e+0 5.5567e-3  1.5984e+4 2.9495e-l  1.9539e+3 6.3287e-2  1.3610e+3 4.4632e-2  4.8474e+8  1.1581e+5 2.0989e-l  7.9797e+4  5.2870e-l  1.3632e-l  2  3.0289e-9  2.3280e-l  1.201 le+4 7.0971e-3  2  5.9693e-9  5.1206e-3  3.5707e-3  2  6.0227e-9  3  1.8657e-8  2.4717e-l  1.2502e+3  3.5098e-4  3  6.7571e-10  1.3907e+0 1.6622e+4 4.3464e-3  3  7.3699e-10 1.8220e+0 1.6196e+4 3.9964e-3  4  2.3904e-9  -  1.1097e+3  1.1774e-4  4  5.2252e-ll  2.0532e-l  1.2672e+5  2.3163e-2  4  5.6698e-ll  1.3792e-l  1.2584e+5 2.1392e-2  5  2.1355e-8  -  5.3334e+2  3.2948e-6  5  7.4005e-ll  4.7657e-2  5.0967e+4  2.5542e-3  5  8.0709e-ll  2.2096e-2  4.8214e+4 2.3490e-3  6  7.1496e-9  -  9.5134e+2  5.8233e-6  6  2.9646e-10 7.9113e-2  1.1106e+4 3.8644e-4  6  3.2834e-10 4.8720e-2  9.6409e+3 3.4997e-4  7  4.4640e-ll  1.5070e+4  7  4.7813e-ll  1.2842e+4 9.9546e-4  High-frequency equivalent capacitance = 1.92682e-10 farad  2.1277e-4  3.0797e+3  1.0495e-3  High-frequency equivalent capacitance = 1.62664e-l 1 farad  7.766 le-3  2.1277e-4  1.7928e+3 3.5436e-3  High-frequency equivalent capacitance = 1.77211e-l 1 farad  2.2 Series Branch Impedance (Z^^) Before the measured impedances are used in the computation of the leakage impedances, the total capacitances derived from thefittinggiven in tables 6.7 and 6.8 have to be deducted from each of them. Equations 5.27, 5.28, 5.29 and 5.38 can then be applied to calculate the individual leakage impedance for each winding. After this, the three leakage impedances for the transformer model can be transformed into delta connection, as shown in fig. 2.9c. Part of the winding-to-winding capacitances can then be combined with the leakage impedances (as shown in fig. 2.10) to obtain the 2yihlSng impedance to be synthesized with multiple RLC equivalent blocks. When this is done, the delta-connected winding impedances can then be converted back to their original configuration ( fig. 9.2a). The manipulations explained above, however, do not always result in winding impedances which are minimum-phase-shift functions. An alternative way was sought so that only minimum-phase-shift functions are dealt with in thefittingof impedances. It was found that the sum of two individual leakage impedances produces a minimum-phase-shift function. Physically, this may correspond to the fact that it is always possible to measure the sum of two leakage impedances of a transformer (if the windings can be disconnected) while it is not possible to directly measure the individual branches of the T circuit. The sum of two winding impedances (for both the zero and positive sequence) were fitted with RLC equivalent networks. The corresponding results are shown in figs. 6.20 to 6.25  132  and tables 6.9 and 6.10. For the reference purposes, the fitted results are named as follows: Zero-sequence Zo_HL: Zo_HT: Zo LT:  Sum of the high-voltage and low-voltage winding impedances. Sum of the high-voltage and tertiary-voltage winding impedances. Sum of the low-voltage and tertiary-voltage winding impedances.  Positive-sequence Zp_HL: Zp_HT: Zp_LT:  Sum of the high-voltage and low-voltage winding impedances. Sum of the high-voltage and tertiary-voltage winding impedances. Sum of the low-voltage and tertiary-voltage winding impedances.  Impedances for both the zero-sequence and the positive-sequence are referred to the high-voltage winding side. The sum of the Z^^  of two windings cannot be used right away since the model  requires a seperate series branch impedance for each winding. However, further fitting is not necessary because Zy/in£ for each winding is related directly to those fitted results explained earlier. In terms of equations, the relationship can be written as follows: Zo HL = Zo H+Zo L  (6.1)  Zo HT = Zo H+Zo T  (6.2)  133  ZoLT  = Zo_L+Zo_T  (6.3)  where, Zo_H, Zo T and Zo_L  = the zero-sequence Z^^^ impedances for the high, tertiary and low-voltage windings, respectively. The impedances are referred to the high-voltage winding side.  Then, it can be shown that  Zo_H  = \(Zo_HL + Zo_HT+ Zo_LT) - Zo_LT  = \Zo_HL + \Zo_HT-\Zo_LT  (6.4)  Similarly, ZoJ  = \Zo_HT+\Zo_LT-\Zo_HL  (6.5)  Zo_L  = \Zo_HL + \Zo_LT-\Zo_HT  (6.6)  and  In terms of circuit representation, one-half impedance can be realised from the original impedance by simply reducing the constituent resistances and inductances to half their values and increasing the capacitances to double their values. Likewise, a negative impedance can be accomplished by replacing each and every parameter of the fitted result with a negative value of an equal magnitude. The winding impedances for the positive-sequence model can be constructed exactly in the same manner.  134  However, if it happens that one of the winding impedances is a minimum-phaseshift function, it can be fitted and used to obtain the other two winding impedances by simple subtraction. For instance, it was found that the winding impedance for the high-voltage winding in the positive sequence of the tested three-winding transformer was a minimum-phase-shift function, as shown in fig. 6.26. Zp_T and Zp L can be found as follows: Zp_T = Zp_HT-Zp_H  (6.7)  Zp_L = Zp_HL-Zp_H  (6.8)  where, ZpTand Zp_L  = the winding impedances of the tertiary and low-voltage windings, respectively.  Zp_H  = the winding impedance of the high-voltage winding which was found to be minimum-phase-shift and was previously fitted.  It is obvious that equations 6.7 and 6.8 are simpler than equations 6.5 and 6.6. Only two impedances are required to be connected in series to make up the winding impedances. Therefore, it is advantageous to use equations 6.7 and 6.8 whenever it is possible. Figures 6.27 and 6.28 show the fitting results for the low-voltage and the tertiary-voltage Z ^ ^ of the positive sequence model computed with equations 6.7 and 6.8, respectively.  135  1,000,000  135 Estimated Phase angle  100,000  Measured  90  10,000  45  1,000  100  -45  10  -90  i  i 111 ii i  10  i  i i 11111  100  i  i  i  i M 11  1,000  10,000  100,000  i  1,000,000  i  i  -135 10,000,000 " i  Frequency (Hz) Fig. 6.26 Fitting of the winding impedance of the high-voltage winding (positive-sequence).  From the numerical examples for both the two-winding and the three-winding transformers presented in this chapter, it can be clearly seen that the proposed technique of impedance derivation to process the experimental data leads to very accurate results. Furthermore, it can be concluded that the developed transformer model is capable of producing satisfactory simulation results both in thefrequency-domainand in the timedomain (experimental data was not available for the three-winding transformer in the time-domain).  3. Simulation of Transient Recovery Voltage To compare the proposedfrequencydependent transformer model with the constant-parameter model, a fault interruption case (Fig. 6.29) was simulated. The plots in Fig. 6.29 compare the voltage across contancts of the circuit breaker obtained with the proposed frequency dependent model (with parameters derived from the tested 50 MVA 115/23 kV transformer) and with a constant parameters model in which the short-circuit impedance i represented with the 60-Hz resistance and inductance. The same external capacitances network is used for both models. The difference in the results illustrate the importance of more accurate transformer modelling in fast switching transients.  137  1,000,000  135 Phase angle  100,000  Estimated  90  Measured  45 0  10,000 -45 -90 1,000  CD CD CD w CD  -135  "Sh  c c« 03 CD  -180 100 -225 -270  10 j  i  i  10  i  i i 11111  100  i  i  i i 111ii  1,000  i  i i 11111  10,000  i  i  i i in  100,000  i  i  i  1,000,000  Frequency (Hz)  Fig. 6.27 Winding impedance of the low-voltage winding (positive-sequence).  i  -315  -360 10,000,000  i i 1111  £  1,000,000  135  100,000  Estimated  90  Measured  45 0  10,000 -45  CO  0>  Phase angle  feb -90  1,000 -135  -180 100 -225 -270  10  -315 -1  |  1  l  l  l I ii I I  10  l  l  I I I ; ;  100  I  I  l  1,000  10,000  100,000  I  I  I  I  1,000,000  Frequency (Hz)  Fig. 6.28  Winding impedance of the tertiary-voltage winding (positive-sequence).  I  -360 10,000,000 I  I I  •5b  a as  £  (a) system voltage  300 200  100 0  •A/^*  transient recovery voltage  300  (b)  200  system voltage  100 0 transient recovery voltage 0  0.2  0.4 0.6 Time (ms)  0.8  Fig. 6.29 Transient recovery voltage in a circuit breaker after clearing a fault, (a) Frequency-dependent model, (b) Constant-parameter model.  140  Chapter Seven PROPOSED FUTURE WORK  Although much effort has been put into the research work for this thesis to make it as complete as possible, there are still a number of issues missing which might be considered worthwhile for continued future research. Three possible areas of further research are envisioned: (1) Modeling of stray capacitances. Stray capacitances in the present model are assumed to be non-frequency dependent and lossless. It might be possible to model stray capacitances in more detailed by introducing dielectric losses into the model, as shown in fig. 7.1. In addition, some of the aspects in the measurement of stray capacitance may need to be improved. Since capacitances are modeled as separate entitiesfromthe rest of the model, balanced conditions for stray capacitances may not need to be assumed. Rn  AA/V R. -o  Fig. 7.1 Capacitance model with dielectric losses. 141  (2) Calculation of model parameters. As an alternative to measurements, it might be possible to determine the model parameters from the transformer physical dimension and construction. These data would have to be obtainedfromthe manufacturer. Although most of the transformer data may be confidential, some manufacturers might still be willing to provide some useful information for educational purposes. (3) Extend the capability of the model to cover any number of windings. It may not require much extra work to make the model valid for transformers having more than three windings. In this more general case, the concept of a T-equivalent circuit to model the coupling among N-coils (N > 3) may have to be abandoned in favor of "delta" branches connected among all nodes (similar to the connection infig.2.9c).  142  Chapter Eight CONCLUSIONS  A wide-band general-purpose model has been developed for three-phase power transformer in the two and three windings per phase. Instead of using a [Y(G>)] matrix formulation for the transformer as seenfromits external terminals, the model uses the classical 60-Hz T-circuit to represent the electric and magnetic interaction among coils belonging to the same phase. For three-phase common-core units, the mutual interaction among different phases is decoupled through a modal transformation matrix. Even though the concept is general, a balanced-system transformation matrix is assumed in order to simplify the modelling and the test data requirements. The decoupled sequence networks consist of a frequency dependent short-circuit branch and constant-valued terminal capacitances. These parameters were measured experimentally on a two-winding and a three-winding, three-phase core-type transformers. The two-winding transformer is rated at 50 MVA, 115/23 kV and the three-winding transformer is rated at 25 MVA , 115/22/11 kv. A technique was suggested that allows positive sequence impedance measurements to be made with a single-phase impedance analyser. The model has been shown to give good results for both two-winding and three-winding transformers. Although the calculation of the series impedances from the 143  test data for the three-winding transformer is more complicated than for the two-winding case, the resulting impedances can still be synthesised with RLC networks with a high degree of accuracy. As a result of the simplified topology, thefrequencydependence modelling problem is reduced to the fitting of simple minimum-phase-shift impedances (the short-circuit impedances). Therefore, the possible numerical stability problems associated with the synthesis of the mutual terms in the [Y(©)] formulation have been eliminated. Also, the model has fewer and simplerfrequencydependent functions to synthesise, making it much faster than other models in time domain simulations. Although the model presented in this thesis depends upon data obtained from measurements, the underlying concepts of the model are not limited to how the data is obtained. Once the data is available, regardless of its source, it can be processed in the same way explained in this thesis to get the parameters of the proposed model.  144  REFERENCE LIST  [1]  Morched, A., Marti, L., Ottevangers, J.: "A High Frequency Transformer Model for the EMTP", IEEE Transactions on Power Delivery, vol. 7, No. 1, July 1992. pp. 1615-26.  [2]  D'Amore, M., and Salerno, M.: "Simplified Models for Simulating Transformer Windings Subject to Impulse Voltage", Paper A 79431-8, presented at IEEE PES Summer Meeting, Vancouver, British Columbia, Canada, July 1979.  [3]  P. T. M. Vaessen.: "Transformer Model for High Frequencies", IEEE Transactions on Power Delivery, vol. 3, No. 4, October 1988, pp. 1761-68.  [4]  R. C. Degeneff.: "A Method for Constructing Terminal Models for Single-Phase n-Winding Transformers", Paper A 78 539-9 presented at IEEE Pes Summer Meeting, Los Angeles California, July 1987.  [5]  T. Adielson, A. Carlson, H. B. Margolis, and J. A. Halladay.: "Resonant Overvoltages in EHV Transfomers - Modelling and Application ", IEEE Transactions on Power Delivery, vol. 1992 PAS-100, No. 7, July 1992 , pp. 3563-72.  [6]  CIGRE WG 33.02 (Internal Voltage).: "Guidelines for Representation of Network Elements When Calculating Transients", CIGRE, 1989.  [7]  R. C. Degeneff.:. "A General Method for Determining Resoances in Transformer Windings", IEEE Transactions on Power Apparatus and Systems, vol. PAS-96, No. 2, March/April 1977, pp. 423-30.  [8]  F. deLeon and A. Semlyen.: "Complete Transformer Model for Electromagnetic Transients", IEEE Transactions on Power Delivery, vol. 9, No. 1, January 1993, pp. 231-9.  [9]  Allan Greenwood.: " Electrical Transients in Power Systems", (Book) 2d ed., New York: John Wiley & Sons, Inc., 1991.  145  [10]  J. R. Marti. 1982.: "Accurate Modelling of Frequency-Dependent Transmission Lines in Electromagnetic Transient Simulations", IEEE Transactions on Power Apparatus and Systems, vol. PAS-101, No. 1, January 1982, pp. 147-55.  [11]  William H. Press, Saul A. Teukolsky, William T. Vetterling and Brian P. Flannery.: "Numerical Recipes in C : The Art of Scientific Computing", 2d ed. Cambridge: Cambridge University Press, 1992.  [12]  Clarke, Edith.: "Circuit Analysis of A-C Power Systems, Volume I: Symmetrical and Related Components", (Book) New York: John Wiley & Sons, Inc., 1943.  [13]  H. Baher.: "Synthesis of Electric Networks", (Book) New York: John Wiley & Sons, Inc., 1984  [14]  Hermann W. Dommel.: "Electromagnetic Transients Program Reference Manual, Vancouver, B.C., By the author, 1986.  [15]  J. R. Marti.: "Modelling of Power Transformers (in Spanish)", Central University of Venezuela, Lecture note (1975).  [16]  Ross Caldecot, Yilu Liu and Selwyn E. Wright. 1990.: "Measurement of the Frequency Dependent Impedance of Major Station Equipment", IEEE Transactions on Power Delivery vol. 5, No. 1, January 1990, pp. 474-80.  [17]  Yilu Liu, Stephen A. Sebo and Selwyn E. Wright.: "Power Transformer Resonance - Measurements and Prediction", IEEE Transactions on Power Delivery, vol. 7, No. 1, January, 1992, pp. 245-53.  [18]  Francisco de Leon and Adam Semlyen.: "Reduced Order Model for Transformer Transients", IEEE Transactions on Power Delivery, vol. 7, No. 1, January 1992, pp. 361-9.  [19]  P. T. M. Vaessen andE. Hanique.: "A New Frequency Response Analysis Method for Power Transformer", IEEE Transactions on Power Delivery, vol. 7, No. 1, January 1992, pp. 384-91.  [20]  Richard L. Bean, Nicholas Chackan, Jr., Harold R. Moore, Edward C. Wentz. 1959.: "Transformers for the Electric Power Industry", (Book) McGraw-Hill, 1959.  146  [21]  L. F. Blume, A. Bayajian, G. Camilli, T. C. Lennox, S. Minneci and V. M. Montsinger.: "Transformers Engineering", (Book) 2d ed. New York: John Wiley & Sons, Inc., 1951  147  Appendix 1 FITTING OF THE IMPEDANCE FUNCTION  The frequency dependent part of the proposed transformer model is approximated with a series connection of several RLC network sections. One such section is shown in fig. 1. Each of the sections is used to duplicate one peak of the frequency response. The RL blocks can be viewed as representing the leakage inductance and associated losses, (a)  R„  AA/V  ,-onhs r^rirs ^vA/V  ^\A/V  R.  R2  (b)  L(oo)  R(G>)  AA/V  Fig. Al. 1 Equivalent network representation of the frequency dependent branch of the sequence model, (a) A section representing any peak, (b) The frequency dependent resistance and inductance equivalent of (§). 148  which are both frequency dependent in nature due to skin effects in the transformer windings and associated physical enclosure. The capacitances of the synthesis blocks represent the interwinding capacitances which become reflected into the leakage impedance for modeling purposes. The combined impedance of the network of fig. 1 can be expressed as:  =  Z(Q)  [R{&) + jd)L(G>)]/JG)C R(G>) +JG)L(6)) + - ^ i?((o) +j(oL((o) [1-G>2CZ(0))] +JG>R(G>)C  (1)  where,  R(.®)  = Ro+t^-t 1=1 ffl2 +pf  Z((D)  = i  RiPi  /=1(02+p2 '  and, EL  Pi  n  L,' >  1 = no. of RL parallel blocks for the principal peak,  =  1  for remaining peaks.  The overall impedance of the frequency dependent branch is the sum of each block's impedance, as expressed in equation 1. To realize the parameters Rff Rf P. and C of each block, it is equally valid to fit either the magnitude or only the real part of the function 149  [1]. Since the impedance is a complex number, finding its magnitude involves operations of both squares and square roots of the real and imaginary parts. Furthermore, the magnitudes of two impedance functions cannot be added directly in order to get the total magnitude, while real parts of two functions can be combined by simple addition. Therefore, it is less complicated to fit the real part of the impedance function than it is to fit its magnitude. The real part of the combined frequency dependent impedance branch can be expressed as the sum of the real part of each impedance block. Multiplying both the denominator and numerator of equation 1 by the complex conjugate of the denominator, the real part of Z(co) is found to be  XeiZfa)}  =  ^  ;  [1 - (Q2CL((0)f + [(Di?(G))Cr (2)  fiCo) where, Q(G>)  =  [1 - ©2CZ(©)]2 +  [G>R(G>)Q2.  The Levenberg-Marquardt method (also called Marquardt method) for nonlinear functions [11] is used in the fitting of the real part of the impedance function. This method has been put forth by Marquardt based on an earlier suggestion by Levenberg. This method works very well in practice and has become the standard in nonlinear leastsquares routines. The algorithm requires, at each iteration, an evaluation of partial derivatives of the model function with respect to each and every parameter on which the model 150  is dependent. The model function for the present problem is simply the sum of the real parts of all the blocks, which may be expressed as  F(<o)  =  £/te{Z,(o>)  (3)  /=i  where, m  = no. of RLC impedance blocks or no. of peaks to be fit.  Although F((Q) contains variables which come from all the impedance blocks, the partial derivatives of F(<o) with respect to the variables of one block will not include any variables from other blocks because all the blocks are independent of each other. The following are the partial derivatives of F((o) with respect to R9 Rp Pt and C which are applicable to any block (from block # 1 to block # m) of the frequency dependent impedance of the transformer model: 1 '  2)  dRo  dFM BRi  32(co Q (a>) dRo  3R(G>)  R(G>) 2  0(G>)  dRo  1 0(G>)  2R2(®)o>2C2 02(o>) 7?(to) dQ(e>) e 2 ((o) dRt  dR((o) 0(co) dRi  1  2  Qfa)  or ®2+P2  ••jgpr- {2[1 - ( 0 2 C Z ( o ) ] r_ f l ) 2 C _l^ ] Q2(G>)  <a2 + P2  151  + 2<B2C2/?(CO)-^-I}  CO2 3)  dPt 3F((o)  _  Q((o) dPt _ j _ dRja)  +P*  g » R(&)  Q(a)  (CO 2 +P 2 ) 2  1  -ItfRtPi  dPt dQ((o)  .^.{2Il-^Ct(a,)][HD»C^] -2(D2fl,P, +2(D C i?(Q)-~, f;}, (co2 +P 2 ) 2 2  2  where,  az(co)  =  co2+/>2  aP/  4v  aF(co) ' 5C  /;,  =  2J?,P2 (co 2 +P 2 ) 2  _ J _ 8R(G>) _ R(<a) a(9(cp) (9(co)' dC (92(co)' 5C  ^ 2• {2©2I(co)[l -co 2 CI(©)] + 2Cco2/?2(©)} Q (o)  In implementing the algorithm, all the peaks are initially approximated with simple blocks of RLC elements. The dc resistance R0 is not initially included. The purpose of this simplification is to find a set of circuit parameters which will be used as an initial guess for further detailed fitting. For these simple blocks, there are only two unknown parameters to be fitted for each block: although there are three variables, i.e., R, L and C, the variables L and C can be assumed to be dependent. They are related to the resonant 152  frequency of the block, ©0, by the well known relationship co02 = (LC)'1. The resonant frequency of each block is readily obtainable from the test data. Therefore, the real part of each block can be reduced from the full expression  Re{Z((o)}  =  ^ ^  (4) +<o2L2  R2(1-<Q2LQ2  down to Re{Z(<o)} =  f  2^ 2  (5) 2  <o RL+ <Q L  R  2 2  which contains only two parameters to be fitted, R and L. The initial value of R is assumed to be the peak value of the real part at the frequency where the resonance of the block occurs. The initial value for L is chosen from the Z(s) function as follow. Let  Z(s)  S  =  -^—~. s +— + — RC LC 2  (6)  Then, in order for Z(co) to have a resonant peak, the roots of the denominator of Z(s) must be complex, or equivalently  2  4  iw) -LC- < ° <.RC;  153  which is the same as  1  RC  (7)  < -h= JIc  At resonance, it is already known that  ©o  =  .  and,  JLC and 1  — C  _  =  „2,  co0^  Substituting into equation 7, one obtains  -£-  < 2©o  which implies that  L < Tlfo 4  (8)  Based on equation 8, an initial value of L is taken as  O.li? ••initial  Tlfo  (9)  After the initial fit, the low frequency region before the principal peak is fiirther refined by adding several RL parallel blocks. After removing the capacitance of the first 154  peak (obtained from the initial global fitting procedure), the asymptotic approximation technique of [1, 10] is used to fit the magnitude rather than the real part of the impedance function. The RL parameter of the initialfrequencyregion, together with the parameter for the peaks from the initial global fitting procedure are combined together as an initial guess for a re-optimization. This second optimization will now include dc resistances for all the impedance blocks. The initial values for these R0 's can be any small numbers and are arbitrarily chosen to be 0.1 ohm each. The correct values will be given by the optimization result. Since numerical stability is an important issue for the electromagnetic transients simulation program, the least-square fits must be constrained so that the result is numerically stable. At the same time, the result must be compatible with the test data. One of the advantages of parameters realization byfittingthe real part of the impedance function is that the real part is always positive for a minimum-phase-shift function. Referring to equation 2, it is seen that Q(G>) is always positive due to the square operations. If the whole function is to be positive, then R{&) must be positive too. One way that this condition can be satisfied, as indicated by equation 1, is to select all i?, as positive numbers. Since the synthesis impedance block consists of two branches, namely, the capacitance branch and the RL branch, the total network will be numerically stable if both branches are numerically stable themselves. The RL parallel network will be numerically stable if its pole (P) is positive. Keeping P. positive and using only positive C 's makes the poles of  155  the rational function approximation positive, therefore guaranteeing that the overall function is numerically stable. The result of keeping P. and Rt positive automatically leads to positive L(.  Having C and Lt positive makes it certain that Q((o) has a minimum  (equation 2) which will in turn cause the function Re{Z((o)} to possess a maximum which is the objective of performing the fitting.  156  

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0065256/manifest

Comment

Related Items