WO2008034936A1

WO2008034936A1 - Method and apparatus for determining circular characteristic

Info

Publication number: WO2008034936A1
Application number: PCT/FI2006/050401
Authority: WO
Inventors: Ari Wahlroos; Janne Altonen
Original assignee: Abb Technology Ag
Priority date: 2006-09-19
Filing date: 2006-09-19
Publication date: 2008-03-27
Also published as: WO2008034936A8; US8022709B2; EP2070170A4; CN101529683B; EP2070170B1; CN101529683A; US20090309612A1; EP2070170A1

Abstract

A method and a system for determining a circular characteristic for distance protection of a three-phase electric line (30), the system comprising means for detecting a fault (F) on the electric line (30), means for identifying a faulted phase or phases of the electric line, means for determining, at a measuring point, a first fault loop impedance by using voltage(s) of the faulted phase(s), a second fault loop impedance by using a polarization voltage and a third fault loop impedance by using predetermined line parameters, and means for determining a radius and midpoint of the circular characteristic.

Description

METHOD AND APPARATUS FOR DETERMINING CIRCULAR CHARACTERISTIC

FIELD OF THE INVENTION

The invention relates to determining a circular characteristic, and more particularly to determining a circular characteristic for distance protection of a three-phase electric iine.

BACKGROUND OF THE INVENTION

A distance relay is a protection relay, which is used to protect eiectric systems and parts thereof, such as electric ϋnes, against faults. A traditiona! method of generating a distance protection function with circular characteristics is to compare the angle between two voltage phasors: an operating voltage phasor S₁ and a polarizing voltage phasor S₂. Phasor S₁ is also commonly known as "the line drop compensated voltage". They take the following form:

S₁ = U_Lx + dir · (I_Lx+I_N ·K_N) ^· Z_1set (eq 2.1 )

S₂ = U_Lxpola (eq 2.2)

where

U_Lx = voltage phasor of faulted phase(s) x. For phase-to- earth elements the voltage is U_L1, U_L2 or U_L3. For phase-to-phase elements the voltage is U_L12> U_L23 or

U_L31. dir = -1 , if operation direction is forward, +1 , if operation direction is reverse. ~ I_Lx = current phasor of faulted phase(s). For phase-to- earth elements the current is I_L1, I_L2 or I_L3. For phase- to-phase elements the current is I_L12, I_L23 or I_L31.

I_N = residual current (I_L1 + I_L2 + I_L3) phasor. This term is zero for phase-to-phase elements. K_N = residual compensation factor = (Z_0set - Z_1set)/(3 ·

Z_1set). This term is zero for phase-to-phase elements.

Z_{1 set} = positive sequence line replica impedance.

Z_0set = zero sequence line replica impedance.

U_Lxpola = polarization voltage. Figure 1 illustrates examples of forward directional self-polarized circular characteristics in two cases: A) a fault inside the protection zone and B) a fault outside the protection zone. Note that since S₁ and S₂ are voltage phasors, the characteristic is drawn in the voltage plane. From Figure 1 it can be seen that the angle a between Sj and S₂ can be used to detect whether the measured impedance lies inside the circle. If the angle α becomes greater than 90 degrees, the measured impedance lies inside the circle and a trip signal should be generated. At the circumference the angle α is 90 degrees. Typically the angle comparison is implemented based on a torque- like relay algorithm utilizing either a cosine or sine phase comparator. These phase comparators emulate the behavior of an induction cup element, the amplitude representing the rotating cup torque and sign rotation direction.

With the cosine-comparator the torque-like equation is:

T_mho = re(S₁) · re(S₂) + im(S₁) - im(S₂) (eq 2,3)

If T_mho < 0 then impedance is inside the zone.

Figure 1 is drawn assuming that the polarizing phasor S₂ was the voltage of the faulty phase. Such polarization method is called "self polarization". The drawback of self polarization is that in case of a close-in fault, the measured voltage may become too smail and trip decision may become uncertain or delayed. Therefore, the polarizing voltage is typically chosen such as to have influence from the healthy phases. The most common types of polarization methods are cross- (or quadrature) polarization and positive sequence polarization. The polarization voltage may additionally include a memorized voltage part prior to fault inception to cope with close-in three phase faults. Table 1 shows cross- (or quadrature) polarization voltages for different fault loops. Table 1

Figure 2 iilustrates voltage triangles of a symmetrical three-phase system. The positive sequence polarization voltage U-j can be calculated based on either phase-to-earth or phase-to-phase voltages as indicated in Tabie 2, which shows positive sequence polarization voltages for different fault ioops.

Tabie 2

An often undocumented feature of circular impedance characteristics is that the selection of polarization voltage affects the shape of the characteristics. In case of short-circuit faults the circle expands as a function of source impedance magnitude, In case of an earth fault, the circle expands as a function of source impedance and earthing impedance magnitude. The magnitude of expansion depends on the poiarization method and it is different for cross-polarization and for positive sequence polarization. An example of circular characteristic expansion is illustrated in Figure 3.

From Figure 3 it can be noticed that the ohmic reach at angle θ is not affected by the expansion. The expansion of the characteristic is in fact a desired feature as it offers a better resistive coverage for close-in faults. Also the directionality is not affected by the expansion.

The circular characteristic is commonly used in distance protection relays. The application is typically short-circuit protection.

For earth-fault protection the desired shape of the protection zone is typically quadrilateral (polygonal), which enables detection of earth-faults with high fault resistance for the total line reach. The distance protection function with quadrilateral characteristics is typically implemented using an impedance mapping approach. This means that an impedance estimate is first calculated and the result is then compared in the impedance plane with operation zone boundaries. A typical quadrilateral characteristic has four or five boundary lines.

Typically distance protection terminals or relays include both circular and quadrilateral characteristics such that circular characteristics are implemented utilizing torque-like algorithms, while quadrilateral characteristics are implemented using the impedance mapping approach. A problem with this kind of arrangement is that the distance protection terminal or relay becomes complex.

BRIEF DESCRIPTION OF THE INVENTION

An object of the present invention is thus to provide a method and an apparatus for implementing the method so as to overcome the above problem or at least to alleviate it. The object of the invention is achieved by a method, a computer program product, a system and an apparatus which are characterized by what is stated in independent claims 1 , 6, 7 and 12. Preferred embodiments of the invention are disclosed in the dependent claims. The invention is based on the idea that fault !oop impedance can be calculated independently from a zone shape calculation also in case of circular characteristics.

An advantage of the invention is that it reduces the required computational burden as fault loop impedance(s) can be calculated centrally. Furthermore, they can be utilized e.g. for both quadrilateral and circular characteristics whereby the circular and quadrilateral characteristics can both be implemented similarly. This is advantageous in case possible short-circuit and earth-fault measuring elements have different trip zone shapes. The invention further enables the circular characteristic expansion to be included in the impedance mapping approach. The invention also allows any kind of modeling for fault impedance calculation. As a result, operation analysis becomes simpler and more straightforward as different characteristics can all be analyzed similarly in the impedance plane.

BRIEF DESCRIPTION OF THE DRAWINGS

In the following the invention will be described in greater detail by means of preferred embodiments with reference to the attached drawings, in which

Figure 1 shows examples of forward directional self-polarized circular characteristics in two cases;

Figure 2 shows voltage triangles of a symmetrical three-phase system;

Figure 3 shows an example of circular characteristic expansion;

Figure 4 shows examples of circular characteristics; Figure 5 illustrates an example of the adjustment of a circular characteristic according to an embodiment of the invention; and

Figure 6 is a diagram illustrating an electric network in which the invention can be used.

DETAILED DESCRIPTION OF THE INVENTION, The use of the method and system of the invention is not limited to any specific system, but they can be used in connection with various three- phase electric systems. The electric line can be a feeder, for example, and may be an overhead-line or a cable or a combination of both. The electric power system in which the invention is implemented can be an electric transmission or distribution network or a component thereof, for example, and may comprise several feeders. Moreover, the use of the invention is not limited to systems employing 50-Hz or 60-Hz fundamental frequencies or to any specific voltage !evef.

Figure 6 is a simplified diagram illustrating an electric network in which the invention can be applied. The figure shows only the components necessary for understanding the invention. The exemplary network can be a network fed through a substation comprising a transformer 10 and a busbar 20. The illustrated network also comprises electric line outlets, i.e. feeders, of which one 30 is shown separately. The figure aiso shows a protective relay unit 40 at the beginning of line 30, and a point of fault F. it should be noted that there may be any number of feeders or other network elements in the network.

System current and voltage values used in the following are preferably obtained by a suitable measuring arrangement including e.g. current and voltage transducers (not shown in the figures) in connection with the electricity system. In most of the existing protection systems such values are readily available and thus the implementation of the invention does not necessarily require any separate measuring arrangements. How these values are obtained is not relevant to the basic idea of the invention and depends on the particular electricity system in which the invention is implemented. According to an embodiment of the invention, the calculation procedure for a circular characteristic for distance protection of a three-phase electric line goes as follows:

First a fault F is detected on the electric line 30 and the faulted phase or phases of the electric line is determined based on information e.g. from a phase selector (not shown). The particular way in which the fault is detected and the faulted phase or phases are determined is not relevant to the basic idea of the invention. Next, at a measuring point which can be at the beginning of the eiectric line or at some other location, the following impedances are determined:

Z_loop = Fault loop impedance calculated by using voltage(s) of faulted phase(s), Z_{loop_polarization} = Fault loop impedance calculated by using polarization voltage, and Z_loopset = Fault loop impedance calculated by using predetermined line parameters The analysis here assumes that the impedance of short-circuit loops includes iine impedance (Z₁) and fault resistance (RF) and the impedance of earth fault loops includes line impedance (Z₁), fault resistance (R_F) and earth return path impedance (Z_N) apart from Z_loopset, which does not include fault resistance. The particular way in which the impedances are determined or calculated is not relevant to the basic idea of the invention. However, one exemplary method for calculating the fault loop impedances (Z_loop) in case of a short-circuit is as follows (d is the per unit fault location):

Z_L12 = U_L12/I_L12 = d·Z₁ + R_FPP (Fault loop L1 to L2) Z_L23 = U_L23/I_L23 = d·Z ₁ + R_FPP (Fault loop L2 to L3) Z_L31 = U_L31/I_L31 = d·Z ₁ + R_FPP (Fault loop L3 to L1 ) Z_L123 = U_L12/I_L12 = d·Z ₁ + R_FPP (Fault loop L1 to L2 to L3)

Where in case of two-phase short-circuit, R_FPP is half of the total physical fault resistance between phases. In case of three-phase short-circuit R_FPP is the physical fauit resistance per phase. Three-phase short-circuit due to its symmetrical nature could be calculated with any fault loop (here it is calculated similarly as the fault loop L1 to L2).

Similarly the fault ioop impedance (Z_loop) in case of an earth fault can be calculated as follows:

Z_{L1 E} = U_L1/I_L1 = d·(Z ₁+Z _N) + R_FPE (Fault loop L1 to earth) Z_L2E = U_L2/I_L2 = d·(Z ₁+Z _N) + R_FPE (Fault loop L2 to earth) Z_L3E = U_L3/I_L3 = d·(Z ₁+Z _N) + R_FPE (Fault loop L3 to earth)

Where R_FPE is the physical fault resistance between phase and earth inciudiπg arc and earthing resistances.

Z_loop can then be obtained from the respective equation of fault loop impedance in question. For example, in case of fault loop from phase L1 to phase L2, Z_loop = Z_L12 = U_L12/ILI2. Z_{loop_polarization} can be obtained from the respective equation of fault loop impedance in question by replacing the faulted phase voltage in the equation with polarization voitage selected from table 3. For example, in case of positive sequence polarization, the fault loop impedance from phase L1 to phase L2 is: Z_{loop_polarization} = U₁·√3.1<30º/I_L12.

In case of a close-in three-phase fault directly in front of or behind the relay location, or other measuring point, the polarization voltage may become too small for secure determination of fault direction. Therefore, the polarization voltage phasor S₂ for three-phase fauit loop should preferably contain a memorized voltage part prior to fault. For example in case of positive sequence polarization, the poiarization voltage could include 15% of memory voltage: S₂ = 0.85-U₁·√3·1 <30º + 0.15-U_1memory·√3·1<30º. Memorized voltage part here refers to voltage prior to fault occurrence, which is stored in a voltage memory (cycϋc buffer).

Z_loopset is preferably obtained from the following equations by using predetermined line parameters:

Z_loopset = Ri+ j-Xi (for short circuit loops) = (R₁+ R_N) + j-(X₁ + X_N) (for earth fauit loops)

where

Ri = positive sequence line resistance

RN = earth return path resistance = (R₀- R₁}/3 Ro = zero sequence line resistance

Xi = positive sequence tine reactance

XN = earth return path reactance = (X₀ - X₁)/3 X₀ = zero sequence line reactance.

Finally, a radius (r) and a midpoint (m) of the circular characteristic are calculated using the following equations:

r = |(K₂ - K₁)/2 | (eq 5.1 ) m = dir·(Z_loopset + (K₁ - K₂)H) (eq 5.2)

where K₂ = polarization impedance phasor = Z_{|oop_polarization} K₁ = operating impedance phasor = Z_loop - Z_loopset dir = +1 , if operation direction is forward; -1 , if operation direction is reverse.

Figure 4 shows examples of circular characteristics; a general case on the left and a self-polarized characteristic on the right. In case of self polarization, Z_{loop_po|arization} = Z_loop, and thus the equations for radius (r) and midpoint (rn) above can be reduced as follows:

r = |0.5· Z_loopset| m = dir·0.5·Z_loopset

In other words, the radius and midpoint can be determined based entirely on predetermined line parameters of the protected electric line. Thus the circle is static in the impedance plane.

In case of cross- or positive sequence polarization, the size and location of the circle are determined by the prevailing system conditions (see Figure 4). In case of short-circuit faults the circle expands as a function of source impedance magnitude (cross-polarization: Z_sourceloop = Z_1source, positive sequence polarization: Z_sourceloop = Z_1source/2). In case of earth fault, the circle expands as a function of source impedance and earthing impedance magnitude (cross-polarization: Z_sourceloop = (2·Z_1source+ Z_0source)/3 + Z_E, positive sequence polarization: Z_sourceloop = (Z_1source+ Z_0source)/3 + ZE₎.

According to an embodiment of the invention the calculated fault loop impedance Z_loop is compared against the geometric shape of the determined circular characteristic; and it is determined whether the calculated fault loop impedance Z_loop is inside the drcfe of the determined circular characteristic. According to an embodiment of the invention, the calculated fault !oop impedance Z_loop - R_loop + j·X_loop is determined to be inside the circie of the determined circular characteristic, if the following equation is true:

(R_loop - Re(m))² + (X_loop - Im(m))² - r < 0.

After the circular characteristic has been determined i.e. the radius (r) and midpoint (m) of the expanded circle have been solved, the characteristic could be freely adjusted in the impedance plane. Figure 5 illustrates an example of an embodiment of the invention in which the circie is moved from the balance point 51 to a new location 52 in order to provide increased coverage for fault resistance. The characteristic can aiso be adjusted in some other way. A system for determining a circuiar characteristic for distance protection of a three-phase electric line according to an embodiment of the invention may be implemented such that it comprises a detection unit which detects a fault on the electric line, an identification unit, which identifies a faulted phase or phases of the electric line, and a calculation unit which determines the impedances and the radius and midpoint of the circular characteristic as described above. Here the term 'unit' refers generally to a physical or logical entity, such as a physical device or a part thereof or a software routine. The other embodiments of the invention described above may be implemented e.g. with the calculation unit or one or more additional units. The above-mentioned detection, identification and calculation units and possible additional units may be physically separate units or implemented as one entity. One or more of these units may reside in the protective relay unit 40 of Figure 6, for example.

When the fault detection and identification of the faulted phase are performed by a separate unit or units which are possibly already present in the system to be protected, the invention can be implemented with an apparatus which receives the fault detection information and indication of the faulted phase from such units. According to an embodiment of the invention, such an apparatus is configured to determine, in response to receiving an indication of a fault on the electric line and of a faulted phase or phases of the electric line, the impedances and the radius and the midpoint of the circular characteristic as described above. Such an apparatus may also be configured to implement the other embodiments of the invention described above.

The system and apparatus according to embodiments of the invention can be implemented by means of a computer or corresponding digital signal processing equipment with suitable software therein, for example. Such a computer or digital signal processing equipment preferably comprises at least a memory providing storage area used for arithmetical operations and a processor, such as a general-purpose digital signal processor (DSP), for executing the arithmetical operations. It is also possible to use a specific integrated circuit or circuits, or corresponding components and devices for implementing the functionaiity according to various embodiments of the invention.

The invention can be implemented in existing system elements, such as various protective relays or relay arrangements, or by using separate dedicated elements or devices in a centralized or distributed manner. Present protective devices for electric systems, such as protective relays, typically comprise processors and memory that can be utilized in the functions according to embodiments of the invention. Thus, all modifications and configurations required for implementing an embodiment of the invention e.g. in existing protective devices may be performed as software routines, which may be implemented as added or updated software routines. If the functionality of the invention is implemented by software, such software can be provided as a computer program product comprising computer program code which, when run on a computer, causes the computer or corresponding arrangement to perform the functionality according to the invention as described above. Such a computer program code can be stored on a computer readable medium, such as suitable memory means, e.g. a flash memory or a disc memory from which it is loadable to the unit or units executing the program code. In addition, such a computer program code implementing the invention can be loaded to the unit or units executing the computer program code via a suitable data network, for example, and it can replace or update a possibly existing program code.

It will be obvious to a person skilled in the art that, as the technology advances, the inventive concept can be implemented in various ways. The invention and its embodiments are not limited to the examples described above but may vary within the scope of the claims.

Claims

1. A method for determining a circular characteristic for distance protection of a three-phase electric line, the method comprising: detecting a fault on the electric line; identifying a faulted phase or phases of the electric line; and determining, at a measuring point, a first fault loop impedance by using voltage(s) of the faulted phase(s), a second fault loop impedance by using a polarization voltage and a third fault loop impedance by using predetermined line parameters, characterized in that the method further comprises: determining a radius r and a midpoint rn of the circular characteristic according to the following equations:

r = I(K₂- K₁)/2| m = dir·(Z_loopset + (K₁ - K₂}/2}, where

K2 ⁼ polarization phasor = Z_loop__polarization Z_loop_poiarization - the second fault loop impedance determined by using polarization voltage K-I = operating phasor = Z_loop - Z_loopset Z_loop = the first fault loop impedance determined by using voltage(s) of the faulted phase(s)

Z_loopset = the third fault loop impedance determined by using predetermined line parameters dir - +1, if operation direction is forward; -1, if operation direction is reverse.

2. A method as claimed in claim 1, characterized in that the method further comprises comparing the determined first fault loop impedance against the geometric shape of the determined circular characteristic; and determining whether the determined first fault loop impedance is inside the circle of the determined circular characteristic.

3. A method as claimed in claim 2, characterized by determining that the determined first fault loop impedance Z_loop = R_loop + j·X_loop is inside the circle of the determined circular characteristic if the following equation is true:

(R_loop - Re(m))²+ (X_loop- lm(m))² - r < 0.

4. A method as claimed in ciaim 1,2 or 3, characterized by adjusting the determined circuiar characteristic.

5. A method as claimed in claim 4, characterized in that the determined circular characteristic is adjusted by adjusting the location of its midpoint in the impedance plane.

6. A computer program product comprising computer program code, wherein the execution of the program code in a computer causes the computer to carry out the steps of the method according to any one of claims 1 to 5.

7. A system for determining a circular characteristic for distance protection of a three-phase electric line, the system comprising: means for detecting a fault (F) on the electric line (30); means for identifying a faulted phase or phases of the electric line (30); and means for determining, at a measuring point, a first fault loop impedance by using voltage(s) of the faulted phase(s), a second fault loop impedance by using a polarization voltage and a third fault loop impedance by using predetermined line parameters, characterized in that the system further comprises: means for determining a radius r and a midpoint m_ of the circular characteristic according to the following equations:

r = l(K₂-K₁)/2| m = dir-(Z_loopset + (K₁ - K2)/2), where

K2 = polarization phasor = Z_loop__polarization Z_loop__polarization = the second fault loop impedance determined by using polarization voltage K₁ = operating phasor = Z_loop - Z_loopset Z_loop = the first fault loop impedance determined by using voltage(s) of faulted phase(s)

8. A system as claimed in claim 7, characterized in that the system comprises means for comparing the determined first fauit loop impedance against the geometric shape of the determined circular characteristic; and means for determining whether the determined first fauit loop impedance is inside the circle of the determined circular characteristic.

9. A system as claimed in claim 8, characterized in that the system comprises means for determining that the determined first fault ioop impedance Z_loop - R_loop + j·X_loop is inside the circle of the determined circular characteristic if the following equation is true:

(R_loop - Re(m))² + (X_loop - im(m))² - r < 0.

10. A system as claimed in claim 7, 8 or 9, characterized in that the system comprises means for adjusting the determined circular characteristic.

11. A system as claimed in claim 10, characterized in that in that the system comprises means for adjusting the determined circular characteristic by adjusting the location of its midpoint in the impedance plane,

12. An apparatus for determining a circular characteristic for distance protection of a three-phase electric fine (30), the apparatus (40) being configured to determine, in response to receiving an indication of a fault (F) on the electric line and_ of a faulted phase or phases of the electric tine, at a measuring point a first fault loop impedance by using voftage(s) of the faulted phase(s), a second fauit loop impedance by using a polarization voltage and a third fault loop impedance by using predetermined line parameters, characterized in that the apparatus (40) is further configured to determine a radius r and a midpoint rn of the circular characteristic according to the following equations:

r= I(K₂-K₁)/2I m = dir-(Z_loopset + (K₁ - K₂)/2), where

K₂ = polarization phasor = Z_{looppolarization} Z_{loop_polarization} = the second fault loop impedance determined by using polarization voltage K₁ = operating phasor = Z_loop - Z_loopset Z_loop = the first fault loop impedance determined by using voltage(s) of faulted phase(s) Z_loopset = the third fault !oop impedance determined by using predetermined line parameters dir = +1, if operation direction is forward; -1, if operation direction is reverse.

13. An apparatus as claimed in claim 12, characterized in that the apparatus is further configured to compare the determined first fault loop impedance against the geometric shape of the determined circular characteristic and to determine whether the determined first fault loop impedance is inside the circle of the determined circuiar characteristic.

14. An apparatus as claimed in claim 13, characterized in that the apparatus is further configured to determine that the determined first fault loop impedance Z_loop = R_loop + j -X_loop is inside the circle of the determined circular characteristic if the following equation is true: (R_loop - Re(m))² + (X_loop - lm(m))² - r < 0.

15. An apparatus as claimed in claim 12, 13 or 14, characterized in that apparatus is further configured to adjust the determined circular characteristic,

16. An apparatus as claimed in claim 15, characterized in that in that the apparatus is further configured to adjust the determined circular characteristic by adjusting the location of its midpoint in the impedance plane.

17. An apparatus as claimed in any one of claims 12-16, characterized in that the apparatus is a protective relay.