Open access

Numerical Methods for Analyzing the Transients in Medium Voltage Networks

Written By

Dumitru Toader, Stefan Haragus and Constantin Blaj

Submitted: 12 May 2012 Published: 19 September 2012

DOI: 10.5772/49953

From the Edited Volume

Numerical Simulation - From Theory to Industry

Edited by Mykhaylo Andriychuk

Chapter metrics overview

3,413 Chapter Downloads

View Full Metrics

1. Introduction

The analysis of transient regimes in electric networks is a complex problem due to great number of elements, some nonlinear, as well as, due to the nonsinusoidal variation of currents and voltages. Numerical simulation can solve properly such a problem

There are available a lot of simulation programs for transients. One of the first of them, dedicated to electro energetic systems, was initiated by H.W.Dommel. Electromagnetic Transient Program (EMTP) was bought by Bonneville Power Administration (USA) (Dommel, 1995). Later on, from this initial program were developed several versions, such as MicroTrans, EMTP-RV, ATP (Alternative Transients Program), PSCAD (Chuco, 2005; DeCarlo & Lin, 2001).

Comparing some of the main simulation programs is not very conclusive, because each one of them are used in some versions trying to solve a large group of problems (Iordache & Mandache, 2004; Istrate et al., 2009).

Interface becomes more and more friendly, the library becomes larger and the facilities for creating own models are simpler to use. The speed of calculations and the storage capacity do no more represent a problem. It is possible now to implement some complicated algorithms able to solve problems with high speed in the variation of variables, able to solve all kinds of nonlinearities or to handle eigenvalues highly distanced one from another (large stiffness ratio dynamical systems).

An important criterion for the selection is represented by the integration of complex dynamical systems like rotating electric machines, relays, power electronic devices, FACTS, controllers etc. in the simulation. Some of the programs, depending on the way that were conceived, allow a better integration of these subsystems.

Choosing one or other of the simulation programs depends on the previous experience of the user, as well as, on the manner the program responds to the specific demands (Chuco, 2005; Danyek et al., 2002; Foltin et al., 2006; Karlsson, 2005; Rashid & Rashid, 2006).

No mater what program is used, the numerical analysis of transients eventually must be able to solve a set of differential – algebraic equations that models the physical system. By accounting the stage the numerical integration takes place at, there are two main classes of simulation programms.

One class consists of those programs that make the integration at element level, meaning that for each step of time discretization the differential equations associated to dynamic elements are transformed in finite difference relations. Several ways of approximation can be used, more frequently the trapeze rule being used. All these algebraic relations are then assembled and as a result an algebraic system of equations is obtained and consequently solved using specific algorithms. This procedure is repeated at each iteration step, the parameters involved in the equations being modified by the results obtained at the previous step. Usually the integration step is fixed, but if this is required by a lack of convergence, the integration step can be split in half and the computation process is resumed. This method is known under the name of the implicit integration method and is used in programs such as SPICE and EMTP (Blume, 1986).

A second class of programs uses a two step procedure. In the first step the mathematic model is expressed as a system of first order differential equations, known as the system’s state equations. During the next step this system of equations is integrated using algorithms with fixed step or with variable steps, depending on the systems particularities. The advantage of this method, called the state variables method, is to treat in the same way, in a unitary manner, electric networks, electric machines, drives, control devices or any other device that allows state equations. A representative of this method is SimPowerSystemsTM which is an extension of Simulink® with tools for modeling and simulating of electrical power systems. It provides models of many components used in these systems, including three-phase machines, electric drives, and libraries of application-specific models such as Flexible AC Transmission Systems (FACTS) and wind-power generation. Harmonic analysis, calculation of Total Harmonic Distortion (THD), load flow, and other key power system analyses are automated (Mathworks®).

The precision of the results obtained by integration of the equations describing the simulated system is quite remarkably high. It is obvious that none of the simulation programs, no mater which of the methods is using, cannot be more precise than the mathematic models used for the simulated components. The accuracy of the values of the parameters of these models is of high importance.

Advertisement

2. The state-variable method for electrical circuits

Medium voltage electric networks are in fact complex circuits mainly made of power sources, resistors, inductors, capacitors and, sometimes, other electric components. Analyzing transient regimes in such networks is the same as the analysis of any other complex electric circuit (Dessaint, et al. 1999; Mandache & Topan, 2009; SIMULINK, 1997).

Generally, an electric circuit, no mater how complex it is, is described by an algebraic linear system of equations, obtained by applying the Kirchhoff's laws. This system of equations reflects the circuit's topology. To complete the model, the voltage-current equation at the terminals of each element of circuit (called also constitutive relations) must be added.

The inductors and the capacitors are described by constitutive relations in which are involved the derivatives of currents, respectively voltages. As a result, the mathematic model of an electric circuit consists of a system of algebraic and differential equations (DAE).

In the electric circuit theory it is demonstrated that each current or voltage can be expressed as function of inductors currents and capacitors voltages. By other words if the inductors currents and the capacitors voltages are known, all remaining currents and voltages are uniquely determined and because of this they are called state-variables.

The state-variable equations are obtained by elimination of the algebraic equations in the initial DAE.

The modelization of an electric circuit made using state-equations has several important advantages.

First of all, electric circuits can be integrated with other dynamic systems, of completely different physical nature, as long as the last ones are described also by state-equations.

By the other hand, the procedures for integrating state-equations have a well-established theoretic support, being the domain of interest for quite a long period of time for famous mathematicians.

Finally, but not last, the method can be applied to linear circuits, as well as, to nonlinear circuits.

Of course there is a price to be paid for all these advantages: a certain difficulty in the elimination of the algebraic equations from the system.

In order to give an example for the basics of the state-variable method, we will analyze the simple RLC series connection circuit supplied by an ideal voltage source u(t). The Kirchhoff equations and the constitutive relations for such a circuit are, respectivelly:

iR(t)=iL(t)=iC(t),u(t)=uR(t)+uL(t)+uC(t),uR(t)=Ri(t),uL(t)=LdiL(t)dt,iC(t)=CduC(t)dt.E1

It is customary to write the last two equation in the form

diL(t)dt=1LuL(t),duC(t)dt=1CiC(t).E2

With the remaining equations uL(t) and iC(t) are obtained as functions of iL(t), uC(t):

uL(t)=u(t)RiL(t)uC(t),iC(t)=iL(t).E3
,

By the substitution in (2) it follows

diL(t)dt=1Lu(t)RLiL(t)1LuC(t),duC(t)dt=1CiL(t),E4

or, in matrix form,

ddt[iL(t)uC(t)]=[RL1L1C0][iL(t)uC(t)]+[1L0][u(t)0].E5

This is a first order differential equations system, with the variables iL(t) and uC(t).

If the initial conditions iL(0) and uC(0) are known, then, according to the uniqueness theorem for the solution of a first order differential equation, the solution iL(t) and uC(t), t0, exists and is unique.

Once the solutions for the variables iL(t) and uC(t) are available, immediately all the other unknowns are solved. Therefore, iL(t) and uC(t) are the state-variables of the circuit, and (5) is the matrix state equation of the circuit.

Generally, for a linear electric circuit, the matrix state equation is of the type (6):

dx(t)dt=Ax(t)+Bu(t),E6

where x(t) is the column matrix of the inductor currents and of the capacitor voltages (state-vector),

u(t) is the column matrix of independent supply sources (sources vector), and, A and B are matrix depending on the circuits topology and of its parameters.

The total number of the state variables is equal with the sum of the number of the inductors and capacitors from the circuit, representing the order of the circuit.

All the other unknown variables are represented by the column matrix y(t) (also called output vector) represented by a linear combination of the state-vector and the sources vector):

y(t)=Cx(t)+Du(t),E7

C and D being matrix of the same origin as A and B.

The equation given by (6) describes the model, in state-space, of a linear electric circuit.

It is possible, as the result of using ideal elements in the modelization, to occur capacitive loops or inductive sections.

Figure 1.

Capacitive loop and, respectively, inductive section

A capacitive loop (see Fig. 1a) is a loop consisting only from capacitors and, possibly, independent voltage sources.

In this case the voltages at the terminals of capacitors from the loop are no more linearly independent variables, because the second Kirchhoff’s law gives:

uS+uC1+uC2+uC3=0E8
.

One, arbitrarily selected, of the capacitor voltages can be expressed as a linear combination of the other voltages.

Such a situation is present also for inductive sections.

An inductive section is represented in Fig.1b, being composed by the convergence in a node of ideal inductors and, possibly, independent current sources. In this case the first Kirchhoff’s law states that

iS+iL1+iL2+iL3=0E9
.

Again anyone of the state variables can be expressed as a function of the others.

The order of an RLC circuit, meaning the number of linearly independent variables (and by consequence the number of state-equations) is equal to nLC - nbC - nsL:

nLC is the total number (sum) of inductors and capacitors

nbC is the number of the capacitive loops

nsL is the number of inductive sections.

The programs performing the analysis make the detection of capacitive loops and of the inductive sections presents in the electric circuit. Usually, it is recommended to introduce a very small value resistor in series connection with one of the elements of the capacitive loop, or to connect a high value resistor in parallel to one of the elements of the inductive section. By this the circuits differs from the original, but the difference in currents values and voltages values are insignificant

The method of state variables can be naturally extended for nonlinear circuits. In this case the circuit equations are of the form

fR(uR,iR)=0,uC=dqdt;fC(uC,q)=0,iL=dΦdt;fL(iL,Φ)=0,E10

with fR,fC,fL being the characteristics of nonlinear elements and, q and Φ, the electric charge of the capacitor, respectively the magnetic flux trough the inductor.

The procedure for obtaining the state equations is similar as for linear circuits, with the difference that the state variables are q and Φ. These variables are present in the dynamic constitutive equations.

The variables uC and iL are obtained as functions of q and Φ and afterwards :

dxdt=F(x,t)E11
,

In (11) x is the column matrix of state variables (state vector), and F is a matrix depending on the circuits topology and on the nonlinear characteristics.

If for linear circuits, the circuits equations can be always reduced to state equations, for certain nonlinear circuits it is not possible to allow state equations (Hasler & Neirynck, 1985).

Anyhow, for a modelization not excessively idealized and a reasonable choice of state variable it is possible to obtain always state equations. For nonlinear state equations even the problem of existing or not of a solution might occur.

State equations for linear circuits have exact analytically determined solution, available in any book of college mathematics. If the order of the circuit is higher than three, analytic solution implies calculations so complicated that this method becomes inefficient. For this kind of situations a large offer of numeric integrations methods is available embedded in the commercial programs of numeric simulation.

Advertisement

3. Simulation of electric circuits using PSPICE

The levels of the analysis of electric circuits performed by PSPICE are, mainly, the following (Radoi, 1994; Vladimirescu, 1999):

  1. Solving a linear electric circuit, invariable with time,

  2. Solving a nonlinear electric circuit,

  3. Solving a time variable electric circuit.

For solving the electric circuit problems PSPICE uses the branch current method. Each element of the circuit is considered being placed between two nodes. The nodes are numbered from “0”, the reference node, to N and the position of an element is given by the numbers of its terminals. A very simple circuit example is given in Fig. 2.

Figure 2.

Simple circuit for SPICE application

{V0V3R3=Igfornode0,V1V2R1=Igfornode1,V1V2R1=V2V3R2fornode2,V2V3R2=V3V0Rfornode3.E12

These equations can be written in matrix form also as (13), where [ I ] is the column matrix of the sources for all the nodes, in the same order as the unknown nodes potentials in matrix [ V ]. [ G ] is the square conductance matrix.

The potential of the reference node, “0”, is taken with “zero” value so all the other N-1 nodes potentials remain as unknown values in the equation system.

With V0 = 0, system (13) is transformed into a system (14) with N-1 equations.

In this very simple example remain only 3 unknown values of nodes potentials (14).

[Ig+Ig00]=[1R3001R301R11R1001R11R1+1R21R21R301R21R2+1R3][V0V1V2V3]or,[I]=[G][V]E13
[+Ig00]=[1R11R101R11R1+1R21R201R21R2+1R3][V1V2V3]E14
.

The system (14) can be solved using the Gaussian elimination method (Ross, 2004).

An equivalent method that can be used for solving the system is the LU factorization, by decomposing the [G] matrix.

Matrix [G]=[L]∙[U], with [L] being a matrix with “1” on the principal diagonal, and [U] is a matrix having only “0” below the principal diagonal. This method is more advantageous than the “classic” Gauss elimination method, mainly when the equation [I] = [G]∙[V] has to be solved several times: several sources [I] for the same conductance matrix. This is the situation for solving our problem using PSPICE, so that the LU factorization method shall be used.

In the analysis of dynamical circuits either the sources, voltage or current sources, have time variation, either the topology of the circuit is changing. If the passive elements are resistors, linear or nonlinear, using a time step ∆t for the discretization of the time variation (sinusoidal, exponential, linear) of the sources. For each time moment a set of values [V(t)] is calculated and memorized. These values can be used in further graph representations or printed.

The presence of inductors and/or capacitors implies replacing them with equivalent elements.

Figure 3.

Thevenin and Norton generators for a capacitor, when trapeze method is used

For an iterative calculation process the capacitor can be replaced by a Thevenin or a Norton generator. The Norton generator presented in the Fig. 3 has the parameters calculated according to the trapeze method of integration. The method using the Euler’s regressive algorithm and the Gear’s second order algorithm can be used also, but the equivalences (parameters of the Norton or Thevenin generators) are different from those presented in Fig. 3.

Figure 4.

Norton equivalent generator representation of a inductor, when trapeze method is used

In a similar manner for an inductor can be used a Norton equivalent generator or a Thevenin equivalent generator, with parameters depending on the method used for solving the problem. Taking the nodes potentials as unknown values it is more suitable to use the Norton generator equivalent to the inductor in a dynamic circuit.

The transient regime analysis is the most important analysis that can be performed using PSPICE.

When the transient regime begins with non zero initial conditions this fact can be introduced in the program by the option UIC (use initial conditions).

The situation of the modelization of a complex network is shown in the Fig. 5.

The transient regime is triggered by closing the switch, meaning that phase one is grounded, as if a fault phase to ground is produced.

Resistor R17 has a value that is corresponding to the conditions of the fault.

Figure 5.

MVN represented as a complex circuit in SPICE drawing

Currents and voltages can be represented, as on a “soft” oscilloscope, after solving the transient regime, using the elements of the V[t] matrix, calculated for each time step, and the resulting, calculated, currents values.

Figure 6.

Current and voltages representation as SPICE output oscillogrames

Advertisement

4. Numeric simulation of transients triggered by faults in medium voltage networks

Faults in medium voltage networks (MVN) generate transient regimes with duration depending on the networks parameters, on the contact type at the fault place and on the default type itself. The numeric simulation performed in this chapter is based on PSPICE and the transient regimes are caused by simple and double grounding faults, as well as, by the fault produced by a broken conductor, grounded towards the consumer. A fault of the type broken conductor grounded towards the source can be treated similar as the simple grounding fault, with an insignificant error, so shall be no more analyzed in the followings (Sybille et al., 2000).

The elements that are important in the analysis are:

  • the initial phase α of the faulty phase,

  • the electric resistance at the fault, grounding point,

  • the method used for grounding the neutral of the network,

  • the functioning regime of the MVN.

4.1. The structure of MVN

The single line diagram of the MVN whose behavior is studied in case of transient caused by faults is presented in Fig. 7.

Figure 7.

Single Line Diagram of the MVN

The MVN in Fig. 7 contains the following elements:

  • The source, the 110kV line supplying the transformer 110kV/20kV, considered as being of infinite power supply and internal impedance zero,

  • The transformer (Tr) 110kV/20kV star connection with neutral point on 110kV side and delta connection on the 20kV side,

  • The lines of 20kV (L1... Ln),

  • The transformers TR.20kV/0.4kV supplying low voltage consumers, transformers having star or delta connection on the 20kV side and star with null or zigzag with null on the low voltage side,

  • Low voltage consumers (C),

  • The transformer for internal services (TIS) ; this transformer is used also for creating the artificial null of the 20kV network (has zigzag connection with null on the 20kV side and star with null on the low voltage side),

  • The grounding inductor (GC), having zigzag with null and creates the artificial null of the 20kV network,

  • The impedance for treatment of the null (Zn), impedance that might be a compensation inductor, a resistor or might be missing in the case of isolated null,

  • The resistor for the null’s treatment (Zn1),

  • The switch for the connection to the ground of the resistor Zn1.

The model is designed with the possibility of changing the values of parameters.

4.2. Numerical simulation of some fault types

The numeric simulator of the three phased circuit corresponding to a MVN is shown in Fig.8. The PSPICE program is used to perform the modelization. Several types of faults are simulated using the model:

  • simple grounding fault,

  • double grounding fault,

  • broken conductor and grounded on the consumer side.

The MVN where the faults are modelized has the following characteristics:

  • nominal line voltage 20kV,

  • S = 25MVA, the apparent power of the transformer 110kv/MV,

  • the capacitive current of the MVN is 100A,

  • the current of the two faulty lines, 8A each,

  • the apparent power of the consumers is 80% of the nominal one and is uniformly distributed on the MV lines.

The simulator contains three lines, denoted with “a”, “b” and “c”, so that to each index of parameter such as resistance or inductance is added a, b or c. The line “c” is considered as the healthy line and is the equivalent of all the 20kV lines without fault (Bucatariu et al., 2009; Drapela, 2009)

The line “a” is provided with the possibility of the simulation of a simple grounding fault, by connecting Vsa. The same line “a” might simulate a fault of the type broken conductor (by opening the switch VSGOL). Even more, by opening the switch VSGOL and closing, in the same time, the switch VSa it is possible to simulate a fault of the type broken conductor and grounded on the consumer’s side.

The line “b” can simulate a simple grounding fault by closing the switch VSb.

The fault of the type double grounding is simulated by closing the switch VSa connected to the line “a” and closing the switch VSb connected to the line “b”, as well.

If a fault of the type double grounding on the same line is wanted to be simulated, the switch VSb should be moved from line “b” to line “a”.

The parameters of the Earth (resistance, inductance) are encompassed in the parameters of the electric line.

Figure 8.

Three phased circuit simulating a MVN

4.2.1. The single phase-to-ground fault

Using the numerical simulator a multitude of situation could be analyzed.

One of those situations is the one when the MVN has the null grounded trough compensation inductor and the MVN functions under resonance conditions, respectively overcompensated 10%.

In these circumstances the evolution of voltages and currents during the transient is analyzed according to:

  • the initial phase of the faulty phase (α) in the very moment when the fault is produced

  • the electric resistance at the fault place (Rt)

  • the way the null is treated

  • the functioning regime (resonance, under or over compensated)

The first oscillogram, left side, is corresponding to the currents and the second corresponding to the voltages.

In these oscillograms I(S) represents the time variation of the current at the fault place, the zero sequence current of the faulty line is denoted as (I(R1a)+I(R2a)+I(R3a))/3, the zero sequence current is (I(R1b)+I(R2b)+I(R3b))/3 on the healty line, and the current trough the element for grounding the neutral point is I(L7).

In the oscillograms of the voltages, V(7), V(8) and V(9) represent, respectively, the voltages of the three phases and the zero sequence voltage of the MVN bars is represented by (V(7) + V(8) + V(9))/3.

In figures 9….20 the results correspond to the 10% overcompensated MVN with the neutral point grounded with compensation inductor. Only two values, 0 and 90°, of the initial phase of the voltage of the faulty line are taken into account. An intermediate value of 45° for this angle does not modify significantly the transient of the simple grounding fault.

From the oscillograms it can be observed that the initial phase of the faulty phase has great importance in the evolution of voltages and currents during the transient regime caused by the simple gounding fault.

Knowing, as accurate as possible, the evolution in time of the voltages gives the possibility of a good design for the MVN, strong enough to resist to the maximal applied voltage.

The maximal values of the currents are important for the correct evaluation of the mechanical forces acting on lines conductors, on the insulators and on the bars of MVN from the transformer station 110kV/MV.

The way that the electric contact is established, Rt, at the fault place, is very important in the transient evolution of the voltages and currents. If the resistance is over 100Ω the transient vanishes rapidly, in less than 5ms, and the maximal values of voltages and currents are not significant. This conclusion results by comparing oscillograms from Fig. 9 to Fig. 19, and, respectively, from Fig. 10 to Fig. 20.

Figure 9.

Time variation of the currents when MVN is 10% overcompensated and Rt = 1Ω, α = 0°

Figure 10.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 1Ω, α = 0°

Figure 11.

Time variation of the currents when MVN is 10% overcompensated and Rt = 1Ω, α = 90°

Figure 12.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 1Ω, α = 90°

Figure 13.

Time variation of the currents when MVN is 10% overcompensated Rt = 10Ω, α = 0°

Figure 14.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 10Ω, α = 0°

Figure 15.

Time variation of the currents when MVN is 10% overcompensated Rt = 10Ω, α = 90°

Figure 16.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 10Ω, α = 90°

Figure 17.

Time variation of the currents when MVN is 10% overcompensated Rt = 100Ω, α = 0°

Figure 18.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 100Ω, α = 0°

Figure 19.

Time variation of the currents when MVN is 10% overcompensated Rt = 100Ω, α = 90°

Figure 20.

Time variation of the voltages when MVN is 10% overcompensated and Rt = 100Ω, α = 90°

If the MVN is functioning at resonance, the same simulation results are presented and, in addition, also the 45° value for the initial phase of the faulty phase voltage is presented.

The oscillograms show that the differences between the 10% overcompensated regime and the resonant regime are not significant.

If the functioning regime is far from being at resonance the transient regim produced by simple grounding fault is very short and the variation of voltages and currents is less important.

The maximal value of the current at fault place is obtained for 90° (comparing α values 0°, 45°, respectively 90°).

For Rt = 1Ω at α = 45° the current is twice its value at 0° and at 90° the current at the fault place is four times greater than the same current if α = 0°.

If Rt becomes 10Ω the conclusions regarding the dependence of the current at the fault place on the phase α remain the same (higher α, higher value of the current) but with a decrease of about 25% of the maximal value of the current for a 10 times increase of the resistance at the fault place.

Figure 21.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 0°

Figure 22.

Time variation of the voltages when MVN is at resonance, Rt = 1Ω, α = 0°

Figure 23.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 45°

Figure 24.

Time variation of the voltages when MVN is at resonance, Rt = 1Ω, α = 45°

Figure 25.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 90°

Figure 26.

Time variation of the voltages when MVN is at resonance, Rt = 1Ω, α = 90°

Figure 27.

Time variation of the currents when MVN is at resonance, Rt = 10Ω, α = 0°

Figure 28.

Time variation of the voltages when MVN is at resonance, Rt = 10Ω, α = 0°

Figure 29.

Time variation of the currents when MVN is at resonance, Rt = 10Ω, α = 45°

Figure 30.

Time variation of the voltages when MVN is at resonance, Rt = 10Ω, α = 45°

Figure 31.

Time variation of the currents when MVN is at resonance, Rt = 10Ω, α = 90°

Figure 32.

Time variation of the voltages when MVN is at resonance, Rt = 10Ω, α = 90°

Figure 33.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 0°

Figure 34.

Time variation of the voltages when MVN is at resonance, Rt = 100Ω, α = 0°

Figure 35.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 45°

Figure 36.

Time variation of the voltages when MVN is at resonance, Rt = 100Ω, α = 45°

Figure 37.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 90°

Figure 38.

Time variation of the voltages when MVN is at resonance, Rt = 100Ω, α = 90°

4.2.2. The double phase-to-ground fault

The same simulation model from Fig. 8 is used with the switch VSa closed, and after 20ms the switch VSb is connected also. VSGOL remains in „closed” position.

The double phase to ground fault was simulated for neutral point grounded by compensation inductor and, respectively, by resistor.

In both cases Rt was considered consecutively with 1Ω, 100Ω, 1000Ω trough VSa switch and with 1Ω trough VSb.

In both cases α = 90° was considered, but also α = 0° was imposed for the situation with compensation inductor.

Also the situation with isolated neutral point was simulated, with same Rt values, but only with α = 90°. (Curcanu et al., 2006).

The same symbols are used in oscillograms for the currents, just in addition the index “a” is used for the first faulty line (simulated by switch VSa) and the index “b” for the second faulty line (VSb). For the voltages the symbols are quite the same as those used previously.

Figure 39.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 0°

Figure 40.

Time variation of the voltages when MVN is at resonance, Rt = 1Ω, α = 0°

Figure 41.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 0°

Figure 42.

Time variation of the voltages when MVN is at resonance, Rt = 100Ω, α = 0°

Figure 43.

Time variation of the currents when MVN is at resonance, Rt = 1000Ω, α = 0°

Figure 44.

Time variation of the voltages when MVN is at resonance, Rt = 1000Ω, α = 0°

Figure 45.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 90°

Figure 46.

Time variation of the voltages when MVN is at resonance, Rt = 1Ω, α = 90°

Figure 47.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 90°

Figure 48.

Time variation of the voltages when MVN is at resonance, Rt = 100Ω, α = 90°

Figure 49.

Time variation of the currents when MVN is at resonance, Rt = 1000Ω, α = 90°

Figure 50.

Time variation of the voltages when MVN is at resonance, Rt = 1000Ω, α = 90°

The double phase to ground fault has less important transient effects than the single phase to ground fault.

When the simple phase to ground fault occurs than the value of initial phase α is importans, but its variation is totaly insignifiant in the double phase to ground fault. The value of the resistance at the fault place is important in both fault types, the maximal value of the fault currents and the zero sequence currents depending strongly on the electric resistance at the fault place.

The oscillograms for the MVN grounded trough resistor are presented in Fig. 51... Fig. 54. The fault is simulated with the switch VSa and for α = 90° the transient is most important.

The current flowing trough the grounding resistor is calculated also and this is important because the value of this current is taken into account for the protection.

The resistance at the fault place is simulated with the switch VSa.

Figure 51.

Time variation of the currents when MVN is grounded trough resistor, Rt = 1Ω, α = 90°

If the resistance at the fault place is increasing from 1Ω to 100Ω the current at the fault place is strongly decreasing from 1.8kA to 500A.

Oscillograms from Fig. 55Fig. 60, for isolated neutral point of the MVN, show that the double phase to ground fault gives a heavier transient for voltages than for currents.

There is not an important difference between the transient produced by the double phase to ground fault in MVN with isolated neutral point compared with the MVN grounded trough compensation inductor, but both cases are much heavier than the transient for MVN grounded trough resistor.

Figure 52.

Time variation of the voltages when MVN is grounded trough resistor Rt = 1Ω, α = 90°

Figure 53.

Time variation of the currents when MVN is grounded trough resistor, Rt = 100Ω, α = 90°

Figure 54.

Time variation of the voltages when MVN is grounded trough resistor Rt = 100Ω, α = 90°

Figure 55.

Time variation of the currents when the neutral point of MVN is isolated, Rt = 1Ω, α = 90°

Figure 56.

Time variation of the voltages when the neutral point of MVN is isolated, Rt = 1Ω, α = 90°

Figure 57.

Time variation of the currents when the neutral point of MVN is isolated, Rt = 100Ω, α = 90°

Figure 58.

Time variation of the voltages when the neutral point of MVN is isolated, Rt = 100Ω, α = 90°

Figure 59.

Time variation of the currents when the neutral point of MVN is isolated, Rt = 1000Ω, α = 90°

Figure 60.

Time variation of the voltages when the neutral point of MVN is isolated, Rt = 1000Ω, α = 90°

The greater values of the currents at the two fault places, when the double phase to ground fault is produced, give high thermic solicitation to the instalation, as well as, higher values of the step voltages. This situatiom is dangerous for human being and animals, especially if this events are situated near the earth plate.

The time variation at the two fault places is far from being the same, no matter what grounding methos is used for the neutral point of the MVN.

4.2.3. The interrupted conductor with ground contact towards the customer fault

For this simulation VSGOL is opened, as if the phase conductor would be interrupted, and with a 20ms delay VSa is connected (simulating the ground contact of the broken conductor). The switch VSb remains opened.

The MVN functioning at resonance, with compensation inductor, was simulated with three values of Rt, namely 1Ω, 10Ω, 100Ω. Two values for α were taken into account, 0°, respectively 90°.

The equivalent capacitance of the phase conductor from the fault to the consumers is simulated by :

  • C101a is simulating the capacitance of phase 1 from the fault place to the bars of the MV transformer station

  • C102a represents the model for the capacitance of the phase conductor from the fault place to the consumers supplied by the faulty line.

Two possible situations are considered: C101a = C102a, and C102a = 0,1C101a.

Figure 61.

Time variation of the currents, neutral point of MVN is grounded by resistor and Rt = 1Ω, α = 90°, C101a = 9C102a

Figure 62.

Time variation of the voltages, neutral point of MVN is grounded by resistor and Rt = 1Ω, α = 90°, C101a = 9C102a

Figure 63.

Time variation of the currents, neutral point of MVN is grounded by resistor and Rt = 100Ω, α = 90°, C101a = 9C102a

Figure 64.

Time variation of the voltages, neutral point of MVN is grounded by resistor and Rt = 100Ω, α = 90°, C101a = 9C102a

In the oscillograms of the currents I(S) is the current to ground at the place of the fault, I(R7) is the current flowing trough the grounding resistor and the zero sequence current of the faulty line is (I(R1a) + I(R2a) + I(R3a))/3. In the oscillogramms (I(R1b + I(R2b) + I(R3b))/3 is the zero sequence current of the healty “b” line.

The phase voltages are V(7), V(8), V(9), the zero sequence voltage on the bars of the MV transform station is (V(7) + V(8) + V(9))/3 and V(71) is the voltage of the broken conductor behind the fault place.

In each situation, either the grounding resistance is 1Ω, or 100Ω, the transient regime is stronger in the very moment of broking the phase conductor (t = 40 ms) than in the moment when the broken conductor contacts the ground (t = 60 ms). The voltage on the broken conductor might reach very high, dangerous values, jeopardizing the insulation of the line as well as the insulation of the consumer.

In Figs. 65 to 78 the oscillograms correspond to the MVN at resonance, with compensation inductor for grounding the neutral point. In this cases the current trough the compensation inductor is not represented because the protections used in MVN do not survey this current.

By the point of view of the overvoltages that might appear when the conductor is broken are more dangerous than the overvoltages in the case when the broken conductor is in contact with the ground.

Figure 65.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 0°, C101a = C102a

Figure 66.

Time variation of the voltages, when MVN is at resonance, Rt = 1Ω, α = 0°, C101a = C102a

Figure 67.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 90°, C101a = C102a

Figure 68.

Time variation of the voltages, when MVN is at resonance, Rt = 1Ω, α = 90°, C101a = C102a

Figure 69.

Time variation of the currents when MVN is at resonance, Rt = 1Ω, α = 90°, C101a = 9C102a

Figure 70.

Time variation of the voltages, when MVN is at resonance, Rt = 1Ω, α = 90°, C101a = 9C102a

Figure 71.

Time variation of the currents when MVN is at resonance, Rt = 1MΩ, α = 90°, C101a = 9C102a

Figure 72.

Time variation of the voltages. when MVN is at resonance, Rt = 1MΩ, α = 90°, C101a = 9C102a

Figure 73.

Time variation of the currents when MVN is at resonance, Rt = 100Ω, α = 90°, C101a = 9C102a

Figure 74.

Time variation of the voltages, when MVN is at resonance, Rt = 100Ω, α = 90°, C101a = 9C102a

When Rt = 1MΩ the fault becomes of the type broken conductor and at t=60ms does not occur a transient regime.

When the fault is of the type interrupted conductor with ground contact towards the customer the heavyest situation is shown to be when the phase is near 90°.

For higher values of Rt the transient regime is less dangerous.

Advertisement

5. Measurement on real MVN

The results obtained during the measurements in the experiment area, while faults were produced on purpose, are shown in the followings.

The fault produced was of type single phase grounding.

The experimental aria contains 7 electric lines of 6kV, with total value of the capacitive current of 27A From the 7 electric linea 6 are connected at first system of bars and at the second bar system was connected the line on which the faults were produced. The capacitive current of this line is 2.7A.

The transformer station contains a transversal switch in connected position.

For recording the currents and the voltages during the experiment was used a CDR oscilloperturbograph.

Figures 75 to 83 show the following:

  • the zero sequence current of the faulty line I0PT Turn,

  • the current trough the compensation inductor IB,

  • the star voltages, Va, Vb, Vc; - the zerosequence voltage of MV bars from the transformer station U0.

The experiment was made in order to verify the concordance with the numeric simulation and to validate the accuracy of the simulator.

The experiment consisted on producing on purpose faults of the type single-phase grounding, namely:

  • metallic grounding of a phase Rt = 0,

  • grounding through the passing resistance Rt= 250 Ω,

  • grounding through the passing resistance Rt= 500 Ω.

The functioning regimes of the MV network were:

  • MV network functioning at resonance; -MV network functions at 10% overcompensation,

  • MV network functions at 10% undercompensation.

The three recordings, corresponding to the three functioning regimes for metallic grounding are shown in Figs. 75, 76 and 77.

Figure 75.

Single phase metallic grounding, Rt= 5 Ω, network at resonance

Figure 76.

Single phase metallic grounding, Rt= 5 Ω, network at 10% over-compensation network

Figure 77.

Single phase grounding with Rt=5 Ω, 10% under-compensated network

By comparing the shape of the current of the faulty line during the transient regime (the zero sequence current of the faulty line) obtained by numeric simulation with the one obtained by recording the same current during experiment it can be observed that the simulated current has an oscillating component greater than in recorded experiment.

This difference is due to a greater equivalent resistance of the experimental equivalent circuit than the simulated one.

In the same maner, as for the metallic grounding, but with a 250Ω grounding resistance for the three functioning regimes of the network the, results are presented in Figs. 78, 79 and 80.

Figure 78.

Single phase grounding with Rt=250 Ω, network at resonance

Figure 79.

Single phase grounding with Rt=250 Ω, 10% over-compensated network

Figure 80.

Single phase grounding with Rt=250 Ω, 10% under-compensated network

The third value for the grounding resistance used in the experiment is 500Ω and the results for the resonant regime, for 10% overcompensated and for 10% undercompensated are presented, respectively, in the figures 81, 82 and 83.

Figure 81.

Single phase grounding resistance, Rt= 500 Ω, network at resonance

Figure 82.

Single phase grounding resistance Rt= 500 Ω, network at 10% overcompensation

Figure 83.

Single phase grounding with Rt=500 Ω network at 10% undercompensation

Table 1, campare some of the values measured during the experiment with the corresponding values obtained by numeric simulation.

The significance of the symbols from the table is as follows:

  • I0max – maximal value of the zero sequence current on the faulty line,

  • IB – current trough the compensation inductor,

  • I0stab – the zero sequence current of the faulty line, after the dumping of the dead-beat component,

  • tam – dumping time of the dead-beat component,

  • U0 – the zero sequence voltage at the MV bars of the transform station,

  • Uf – voltage of the healty lines.

From Table 1. results that the differences between the values obtained by the numerical simulation and the values measured during the experiment are reasonably small.

Regime of the networkSymbol of
quantities
Rt[Ω]
5250500
Values obtained byError
[%]
Values obtained byError
[%]
Values obtained byError
[%]
Experi-
ment
Simu-
lation
Experi-
ment
Simu-
lation
Experi-
ment
Simu-
lation
ResonanceI0max[A]21,724,613,44,84,56,32,82,67,1
IB [A]26,1624,37,117,8116,29,114,3712,8410,6
I0stab[A]2,282,426,11,941,759,81,641,499,1
tam [ms]1261303,2383410,557545,3
U0 [kV]3,63,551,42,622,437,32,312,176,1
Uf [kV]6,226,033,15,65,462,55,34,789,8
10%
Over-compensated
I0max[A]17,8582264,24,57,13,22,99,4
IB [A]26,6924,827,019,618,36,615,1613,868,6
I0stab[A]2,22,066,42,642,428,32,532,2710,3
tam ms]80867,535328,641379,8
U0 [kV]3,633,552,22,732,479,51,911,79,9
Uf [kV]6,245,964,55,695,218,44,384,768,7
10%
Under - compensated
I0max[A]20,99431054,34,54,73,12,89,7
IB [A]22,8221,286,717,7716,57,19,449,222,3
I0stab[A]2,52,327,22,051,916,82,532,318,7
tam ms]1982106,131288,623218,7
U0 [kV]3,593,560,842,82,539,61,841,679,1
Uf [kV]6,25,973,75,715,395,65,354,8110,1

Table 1.

Comparison of measured and simulated values

Advertisement

6. Conclusions

The numerical simulator designed by us allows the analysis of transients caused by diferent types of faults, such as simple groundings, double groundings, or broken conductor grounded towards the consumer.

We have compared results obtained by simulation with measured values only for simple grounding faults.

In this type of fault the simulator is validated by the measurements, no matter how the neutral point of the medium voltage network is grounded and whatever are the functioning conditions of the electrical network.

The results from Table 1. show that the simplificatory conditions, taken into account for developing the simulator, are correct.

The numerical simulation of the transient regimes caused by simple grounding faults produced in medium voltage networks shows to be an efficient method for analyzing such faults.

The most dangerous transient regimes occur when the initial phase of the voltage of the faulty line is near 90°.

The initial phase of the voltage, that happen to be in the very moment of producing the fault, was taken with the same value in simulation. The differences between the maximal values of zero sequence component, obtained by the two methods, are quite acceptable from a technical point of view and this concordance validates the model.

The resistance of the broken conductor to ground is extremely important, small values of Rt implying long damping periods. If this undesirable condition happens the currents might have high values, as well as important values of the over voltages leading to important supplementary mechanical stress and damages of the insulating devices and, eventually a simple fault can turn to a multiple fault situation.

For small values of the grounding resistance at the fault location (Rt < 10Ω) difference between measured and simulated values is a litle higher than in the case of greater values (Rt > 100Ω).

The model is very useful by giving the values of the voltages on the healthy lines (over voltages that might jeopardize the insulation) and zero sequence currents of the faulty lines (fixing the condition of the fault detection by protection devices).

Using the conclusions of the simulation of a MVN for the values of the currents and voltages during different types of faults it is possible to adjust the prescribed values of the protection.

Simulation of the broken conductor fault, either connected to the ground or not, shows that resonance might occur, over voltages due to this resonances being dangerous to the equipments.

When double phase to ground fault was simulated the values of the voltages showed dangerous voltages for step or touch voltage.

The numerical simulator has the advantage of analyzing rapidly several variants as models of different possible situations.

The simulator is flexible and all kinds of faults can be simulated for different MVN, just by making the proper modification in the parameters of the simulator. These modifications are slightly easier to be performed when PSPICE simulating model is used.

Either using PSPICE or Mathlab-Simulink the accuracy of the analysis is similarly good.

The precision of the results obtained by simulation depends essentially on the accuracy of the MVN parameters.

References

  1. 1. BlumeW.1986Computer Circuit Simulation, BYTE, July 1986
  2. 2. BucatariuI.SurianuF.2009Optimal Placement of Fixed Series Capacitor in Distribution Networks, Proceedings of the 9thWSEA/IASMES International Conference on Electric Power Systems High Voltages, Electric Machines Control & Signal Processing, Genova, Italy, October 17-19, 2009, 6570
  3. 3. ChucoB.2005Electrical Software Tools Overview in SINATEC-IEEE
  4. 4. CurcanuG.ToaderD.PandiaT.2006Determination of overvoltages in high voltage networks at single phase faults by numarical simulation and experiments, Proceedings 12th International IGTE Symposium, Gratz
  5. 5. De CarloR.LinP.2001Linear Circuit Analysis, Oxford University Press
  6. 6. DanyekM.HandlP.RaiszD.2002Comparison of Simulation Tools ATP-EMTP and MATLAB-Simulink for Time Domain Power System Transient Studies, Proceedings of the European EMTP-ATP Conference, Sopron, Hungary
  7. 7. DessaintL.A.Al-HaddadK.Le -HuyH.SybilleG.BrunelleP.1999A power system simulation tool based on Simulink, IEEE Trans. on Ind. Electronics, 46612521254Dec. 1999
  8. 8. DommelH. W.1995ElectroMagnetic Transient Program. Theory Book, Bonneville Power Administration, Portland
  9. 9. DrapelaJ.2009Performance of a Voltage Peak Detection-Based Flichermeter, Proceedings of the 8thWSEAS International Conference on Circuits, Systems, Electronics, Control & Signal Processing, Teneriffe, Spain, December 14-16, 2009, 296301
  10. 10. FoltinM.ErnekM.HnatJ.2006SimPowerSystems in Education, Mezinarodni Konference Technical Computing, Prague
  11. 11. HaslerM.NeirynckJ.1985Circuits non lineares, Presses Polytechniques Romandes, Lausannes.
  12. 12. IordacheM.MandacheA.2004Computer Aided Analyse of Nonlinear Circuits (in Romanian), Politehnica Press Bucharest
  13. 13. IstrateM.GavrilasM.IstrateC.UrsuleanuR.2009Algorithms in Transmission Grids Using ATP Simulations, Proceedings of the 9thWSEA/IASMES International Conference on Electric Power Systems High Voltages, Electric Machines Control & Signal Processing, Genova, Italy, October 17-19, 2009, 109114
  14. 14. KarlssonA.2005Evaluation of Simulink/SimPowerSystems and other Commercial Simulation Tools for the Simulation of Machine System Transient, Swedish Royal Institute of Technology, Stockholm
  15. 15. MandacheL.TopanD.2009Algorithms for Electric Circuit Simulation, Universitaria Press Craiova
  16. 16. RadoiC.1994SPICE Simulation and Analysis of Electronic Circuits, Ameo Press Bucureşti.
  17. 17. RashidM. H.RashidH. M.2006SPICE for Power Electronics and Electric Power, Taylor & Francis Group, Boca Raton
  18. 18. RossC. C.2004Differential Equations, Springer Science, 5156
  19. 19. SybilleG.BrunelleP.Hoang-Huy LeDessaint. L. A.Al-HaddadK.2000Theory and applications of power system blockset, a MATLAB/Simulink-based simulation tool for power systems, Power Engineering Society Winter Meeting 2000, IEEE, 1774779
  20. 20. VladimirescuA.1999SPICE, Ed.Tehnica, Bucuresti
  21. 21. VladimirescuA.1997SIMULINK ; Dynamic System Simulation for MATLAB, User’s Guide, version 2.1, MathWorks Inc.

Written By

Dumitru Toader, Stefan Haragus and Constantin Blaj

Submitted: 12 May 2012 Published: 19 September 2012