Industrial Electrical Engineering and Automation

Similar documents
POLYPHASE CIRCUITS. Introduction:

Network Analysis and Synthesis. Chapter 5 Two port networks

ELE B7 Power Systems Engineering. Power System Components Modeling

ELE B7 Power System Engineering. Unbalanced Fault Analysis

u( t) + K 2 ( ) = 1 t > 0 Analyzing Damped Oscillations Problem (Meador, example 2-18, pp 44-48): Determine the equation of the following graph.

2.4 Linear Inequalities and Interval Notation

1B40 Practical Skills

332:221 Principles of Electrical Engineering I Fall Hourly Exam 2 November 6, 2006

Math 259 Winter Solutions to Homework #9

Review of Gaussian Quadrature method

Designing Information Devices and Systems I Anant Sahai, Ali Niknejad. This homework is due October 19, 2015, at Noon.

Designing Information Devices and Systems I Discussion 8B

Designing Information Devices and Systems I Spring 2018 Homework 7

Special Relativity solved examples using an Electrical Analog Circuit

set is not closed under matrix [ multiplication, ] and does not form a group.

p-adic Egyptian Fractions

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

Lec 3: Power System Components

5.4 The Quarter-Wave Transformer

1. For each of the following theorems, give a two or three sentence sketch of how the proof goes or why it is not true.

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Things to Memorize: A Partial List. January 27, 2017

Chapter E - Problems

Lecture 3: Equivalence Relations

Chapter 9 Definite Integrals

EECS 141 Due 04/19/02, 5pm, in 558 Cory

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17

Quadratic Forms. Quadratic Forms

Designing Information Devices and Systems I Fall 2016 Babak Ayazifar, Vladimir Stojanovic Homework 6. This homework is due October 11, 2016, at Noon.

Analytically, vectors will be represented by lowercase bold-face Latin letters, e.g. a, r, q.

Interpreting Integrals and the Fundamental Theorem

4.1. Probability Density Functions

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17

Physics 1402: Lecture 7 Today s Agenda

1 Nondeterministic Finite Automata

Vectors , (0,0). 5. A vector is commonly denoted by putting an arrow above its symbol, as in the picture above. Here are some 3-dimensional vectors:

Solution of Tutorial 2 Converter driven DC motor drive

September 13 Homework Solutions

Section 6: Area, Volume, and Average Value

10. AREAS BETWEEN CURVES

Chapter 10: Symmetrical Components and Unbalanced Faults, Part II

9.4. The Vector Product. Introduction. Prerequisites. Learning Outcomes

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Parse trees, ambiguity, and Chomsky normal form

Continuous Random Variables Class 5, Jeremy Orloff and Jonathan Bloom

New data structures to reduce data size and search time

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Reading from Young & Freedman: For this topic, read the introduction to chapter 24 and sections 24.1 to 24.5.

Matching patterns of line segments by eigenvector decomposition

Describe in words how you interpret this quantity. Precisely what information do you get from x?

A Thermal Model of a Forced-Cooled Heat Sink for Transient Temperature Calculations Employing a Circuit Simulator

8 THREE PHASE A.C. CIRCUITS

Chapter 6 Techniques of Integration

2. VECTORS AND MATRICES IN 3 DIMENSIONS

5 Equivalent Circuits and Parameters of Power System Plant

MA123, Chapter 10: Formulas for integrals: integrals, antiderivatives, and the Fundamental Theorem of Calculus (pp.

Improper Integrals. The First Fundamental Theorem of Calculus, as we ve discussed in class, goes as follows:

Intermediate Math Circles Wednesday, November 14, 2018 Finite Automata II. Nickolas Rollick a b b. a b 4

5: The Definite Integral

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Designing Information Devices and Systems I Spring 2018 Homework 8

Farey Fractions. Rickard Fernström. U.U.D.M. Project Report 2017:24. Department of Mathematics Uppsala University

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Section 4: Integration ECO4112F 2011

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Chapter 4: Techniques of Circuit Analysis. Chapter 4: Techniques of Circuit Analysis

Coalgebra, Lecture 15: Equations for Deterministic Automata

Lecture 2e Orthogonal Complement (pages )

Dorf, R.C., Wan, Z. T- Equivalent Networks The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000

Polynomials and Division Theory

Math& 152 Section Integration by Parts

13: Diffusion in 2 Energy Groups

Chapters Five Notes SN AA U1C5

This chapter will show you. What you should already know. 1 Write down the value of each of the following. a 5 2

AN020. a a a. cos. cos. cos. Orientations and Rotations. Introduction. Orientations

MTH 505: Number Theory Spring 2017

The Periodically Forced Harmonic Oscillator

1.2 What is a vector? (Section 2.2) Two properties (attributes) of a vector are and.

Linear Inequalities. Work Sheet 1

Note 12. Introduction to Digital Control Systems

8. Complex Numbers. We can combine the real numbers with this new imaginary number to form the complex numbers.

Torsion in Groups of Integral Triangles

COSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature

Designing finite automata II

Lecture 3. Introduction digital logic. Notes. Notes. Notes. Representations. February Bern University of Applied Sciences.

2 b. , a. area is S= 2π xds. Again, understand where these formulas came from (pages ).

The Properties of Stars

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Lecture 6. Notes. Notes. Notes. Representations Z A B and A B R. BTE Electronics Fundamentals August Bern University of Applied Sciences

Homework Solution - Set 5 Due: Friday 10/03/08

1 1D heat and wave equations on a finite interval

STRAND J: TRANSFORMATIONS, VECTORS and MATRICES

7.2 The Definite Integral

Assignment 1 Automata, Languages, and Computability. 1 Finite State Automata and Regular Languages

SOME INTEGRAL INEQUALITIES OF GRÜSS TYPE

PH12b 2010 Solutions HW#3

Transcription:

CODEN:LUTEDX/(TEIE-719)/1-7/(7) Industril Electricl Engineering nd Automtion Estimtion of the Zero Sequence oltge on the D- side of Dy Trnsformer y Using One oltge Trnsformer on the D-side Frncesco Sull Dept. of Industril Electricl Engineering nd Automtion Lund University

Tle of Contents 1 INTRODUCTION...3 THEORY...3 3 ERIFICATION OF THE METHOD THROUGH SIMULATIONS IN MATLAB SIMPOWER...5 4 LABORATORY TESTS AND RESULTS...6

1 Introduction This report dels with the estimtion of the zero sequence voltge on the D-side of Dy trnsformer y the use of only one voltge trnsformer (T) connected etween one phse nd ground on the D-side. Zero sequence voltge cn e used for phse to ground fult detection in non-directly erthed systems. Normlly, on the D-side of trnsformer the zero sequence voltge is mesured y the use of three Ts, ech of them connected etween ech phse nd ground on the primry side nd open delt connected on the secondry side. The method descried in this report llows D-side zero sequence voltge estimtion y the use of one T on the D-side nd three Ts on the y-side. If the y-side is t lower nominl voltge, the descried method could e ttrctive on n economic point of view. Though the descried method is generl nd works for every Dy trnsformer (nd with some djustments for ny trnsformer connection), the cses where it is economiclly dvntgeous re few. An exmple is the cse where distriuted genertor is connected to medium voltge non-directly erthed network through Dy step-up trnsformer nd zero sequence voltge on the D-side is necessry for protection purposes. 3 Theory The min ide ehind the method descried in this section for estimting the zero sequence voltge on the D-side is to use only one T on the D-side nd to use the informtion from three Ts connected on the y-side of the trnsformer. The estimtion of the zero sequence voltge is performed s descried elow. Two premises re necessry (Zorn Gcjc s PhD thesis): - Dy trnsformer trnsforms sequence quntities independently - positive nd negtive sequence re trnsformed ccordingly to: N j j = ϑ ϑ 1 y1e τy1e N = y (1) N = e = τ e jϑ y y N y jϑ where - ϑ is the ngle y which the primry-secondry quntities re shifted in the trnsformtion nd it is determined y the trnsformer group numer - τ is the trnsformer rtio - N nd N y re the numer of windings on the D- nd y-side respectively

4 Zero sequence quntities re not trnsformed in cse of Dyn trnsformer. Eqution (1) does not tke into ccount the voltge drop inside the trnsformer. This cn e esily ccounted for y modifying Eqution (1) s follows: τ e τz I e τe Z I = jϑ jϑ = jϑ () 1 y1 tr y1 y1 tr y1 = τ e τz I e = τe Z I jϑ jϑ jϑ y tr y y tr y where - Z tr is the trnsformer impednce s seen from the low voltge y-side If we express the phse voltges on the D-side in terms of their sequence components, we get: 1 1 1 1 1 = c 1 (3) In this eqution the unknown is the zero sequence voltge. 1 nd re otined from Eqution () through the informtion collected on the y-side. By inserting Eqution () into Eqution (3) we get: ( y1 tr y1) jϑ 1 1 1 1 1 1 1 τe Z I j 1 1 ϑ τ e y ZtrIy = = c 1 1 The previous is system of three equtions. Solving the first of these equtions for we find: jϑ τ = τe Z I e Z I (5) jϑ y1 tr y1 y tr y This is n estimtion of the zero sequence voltge on the D-side nd voltge is mesured. Any of the three equtions in (4) cn e used for the estimtion, ech requiring different phse voltge to ground to e mesured on the D-side. The zero sequence voltge on the D-side cn e ccurtely estimted if: - positive nd negtive sequence currents nd voltges on the low voltge y-side re known - step-up trnsformer impednce, rtio nd ngle group re known - one phse to ground voltge on high voltge D-side is known The dvntge of this method is the use of only one T on the high voltge D-side to estimte the zero sequence voltge rising fter ground fult in non-directly erthed system. The mjority of the distriuted genertors re lredy provided with current nd voltge trnsformers on the low voltge y-side, i.e. there re no dditionl costs to get the positive nd negtive sequence quntities used in (5). However, to dte, distriuted genertors re not commonly operted in islnd nd therefore, sensing the zero sequence voltge on the D-side of the step-up trnsform is not lwys needed. Insted, if islnd opertion is desired, sensing the zero sequence (4)

5 voltge on the D-side cn e the most prcticl wy to detect ground fults on tht side. The word estimtion hs een used through ll the report insted of clcultion. The zero sequence voltge cn only e estimted ecuse of the uncertinties in trnsformer prmeters nd mesurement errors from CTs nd Ts involved in Eqution (5). Another source of error in the estimtion derives from neglecting the trnsformer mgnetizing current in Eqution (5). However, s it is shown in the following, the estimtion cn e very ccurte if the uncertinties re low. Also, the mgnetizing current of trnsformer is usully very smll. The physicl mening of the method is simple. On the y-side, ny informtion is lost on the neutrl point s (which is not physicl point) shift on the D-side. By mesuring one phse-to-ground voltge on the D-side we regin this informtion. All informtion out positive nd negtive sequence voltges is identicl on oth sides of the step-up trnsformer. Therefore the zero sequence voltge on the D-side cn e estimted. 3 erifiction of the method through simultions in MATLAB SimPower The method hs een implemented in MATLAB SimPower nd it proved to e ccurte. The estimted zero sequence voltge hs een compred with the mesured one on the D-side of the trnsformer, during phse to ground fult on tht side t k voltge level. Results from simultions re shown elow. x 14 1.5 [] 1.5 5 error [] Figure 3.1-5 4 6 8 1 1 Time [s] Simultion of phse to ground fult on the D-side. ) mesured (dshed) nd estimted (dotted) zero sequence voltge. ) error etween the estimted nd the mesured zero sequence voltge

6 The mesured nd estimted zero sequence voltges re prcticlly undistinguishle. The mximum error is round 5, i.e. round,35% of system nominl phse voltge nd it cn e totlly neglected. The error depends on neglecting the mgnetizing current in the trnsformer. In the simultions it hs een ssumed mgnetizing current equl to,1 pu, i.e. mgnetizing impednce of 1 p.u. Angle errors, not reported, re lso very smll. 4 Lortory tests nd results Some tests hve een performed in the l to verify the method descried ove for estimting the zero-sequence voltge on the d-side of Ynd trnsformer. The fcts tht the d-side is the low voltge side nd tht the Y-side is grounded do not ffect the conclusions. The lortory experiment is descried in the figure elow. Figure 4.1 Set-up of the lortory experiment An Ynd trnsformer is fed through tp-chnger trnsformer from the net. The voltge on the Yn side is regulted so to get 11 phse-to-phse voltge on the d- side. On the d-side resistive lod is fed nd cpcitnces re connected so to simulte long line. On the Yn side ll phse currents nd voltges re mesured. On the d-side phse to ground voltge is mesured. Also the zero-sequence voltge is mesured on the d-side through three open delt connected voltge trnsformers, not shown in the figure for ske of clrity. All mesurements re fed in DSpce nd phsors clcultions re performed in DSpce through one-cycle Fourier lgorithm. Finlly, the new lgorithm for zero sequence estimtion is lso implemented in DSpce. A phse to ground fult is simulted on the d-side y connecting one phse directly to the cpcitnces neutrl point. Results elow show the mesured zero-sequence voltge, the zero sequence voltge estimted through the new lgorithm nd the fulted phse voltge. It is noted tht efore fult occurrence, the ngles of the estimted nd mesured zero sequence voltge hve ctully no mening, eing the voltge idelly equl to zero.

7 It is seen tht the zero sequence voltge estimted through the lgorithm reproduces firly the mesured zero sequence voltge, lredy fter less thn cycle. After the trnsient, the two voltges re prcticlly equl. The smll differences re minly due to uncertinties in the converting fctors for DSpce input signls. 8 6 [] 4 [degrees] -5-1 -15..4.6.8.1.1 time (s) Figure 4. Phse to ground fult on the d-side. ) Phse voltge (solid), mesured (dshed) nd estimted (dotted) zero sequence voltge. ) mesured (solid) nd estimted (dshed) zero sequence voltge ngle 15 1 [] 5-5 [degrees] -1-15 -..4.6.8.1.1 time (s) Figure 4.3. Phse to ground fult on the d-side. ) Phse voltge (solid), mesured (dshed) nd estimted (dotted) zero sequence voltge. ) mesured (solid) nd estimted (dshed) zero sequence voltge ngle