new principle for transmission line protection using phase portrait plane

Upload: batchuraja

Post on 10-Apr-2018

213 views

Category:

Documents


0 download

TRANSCRIPT

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    1/8

    Published in IET Generation, Transmission & Distribution

    Received on 11th October 2007

    Revised on 10th September 2008

    doi: 10.1049/iet-gtd:20080456

    ISSN 1751-8687

    New principle for transmission line protection

    using phase portrait planeM.M. Eissa*Department of Electrical Engineering, Helwan University, Helwan, Cairo, Egypt

    *Projects Department, King Abdulaziz University, Jeddah, P.O. Box 80200, The Kingdom of Saudi ArabiaE-mail: [email protected]

    Abstract: Phase portraits are a powerful mathematical model for describing oriented textures. An isotangent-

    based approach in a phase portrait is introduced to discriminate between internal and external faults. The

    scheme described is new and is a different approach to the problem of relaying ground faults on transmission

    lines. The geometrical theory of differential equations is used to drive a symbol set on the basis of the visual

    appearance of phase portraits. The phase portrait of the instantaneous rate of change of current against

    voltage gives valuable information on the transient and stability characteristics of system configuration with

    high fault resistance and shunt susceptance of the transmission lines. Hundred percent of the protected zone

    is successfully discriminated rather than the external zone. An additional parameter that is considered, which

    is sometimes neglected in protection studies, is the shunt susceptance of the transmission lines.

    1 Introduction

    Transmission line protection is generally regarded as themost complex of all relaying applications. This is because ofthe distances involved and the variety of possibilities thatexist in system configuration and opportunities for bothcorrect and incorrect operations. A distance relay isdesigned to only operate for faults occurring between therelay location and the selected reach point, and remainsstable for all faults outside this region or zone [1].

    In developing distance relay equations, the fault underconsideration is assumed to be ideal (i.e. zero resistance)[210]. In reality, the fault resistance will be between twohigh-voltage conductors, whereas for ground faults, the faultpath may consist of an electrical arc between the high-voltageconductor and a grounded object. The fault resistanceintroduces an error in the fault distance estimate and, hence,may create an unreliable operation of a distance relay [11].

    The impedance seen by the relay is not proportional to thedistance between the relay and the fault, in general, becauseof the presence of resistance at the fault location.

    Some techniques for arcing faults detection and faultdistance estimation are introduced in [12, 13]. The

    techniques are based on the voltage and current at oneterminal in the time domain. The overhead line parametersand arc voltage amplitude during fault are given. Thetechniques have optimal application in the medium-voltagenetworks and symmetrical faults. Some techniques [1419]are suggested for enhancing the high fault resistanceproblem. These techniques accommodate this problem byshaping the trip zone of the distance relay, to ensure thatthe apparent impedance is included inside the trip zone.New techniques using high-frequency components of thefault-generated transient signals are also studied in [2022].

    Distance-based techniques are applied using the currentand voltage at the relaying point by differential equations.

    The fundamental assumption used in this approach is thatthe transmission line can be represented by either a lumped-series impedance or a single PI-section. An additionalparameter that is sometimes neglected in protection studiesis the shunt susceptance of the transmission lines. Forthe EHV lines, the shunt susceptance is rather large. If thesusceptance is omitted from the study, the results areconsidered to be errors.

    In this paper, a new isotangent-based approach in a phaseportrait to discriminate between an internal and external faultis introduced. Phase portrait analysis is done to obtain full

    IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 49 56 49

    doi: 10.1049/iet-gtd:20080456 & The Institution of Engineering and Technology 2008

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    2/8

    insight into the fault trajectory during the steady state andtransients concerning the effect of high fault resistance andshunt susceptance. The geometrical theory of differentialequations to describe the transmission line during faultis obtained. The instantaneous rate of change of currentagainst voltage in a phase portrait is given. Hundred percent

    of the protected zone is successfully discriminated.

    2 Phase portraits

    Phase portraits provide an analytical tool to study systems offirst-order differential equations [23]. The method hasproved to be useful in characterising oriented texture [24,25]. Let _x and _y denote two differentiable functions oftime t, related by a system of first-order differentialequations as

    _x F(x, y) cx dyf (1)

    _y G(x, y) ax by e (2)

    to describe a set of oriented textures comprising saddles,star-nodes, nodes, improper nodes, centres and spirals(Fig. 1). Equations (1) and (2) can be represented in matrixnotation as

    _~X ~A~X ~B (3)

    where

    ~X xy ,

    ~Aca

    db and

    ~Bfe

    ~A is the characteristic matrix of the system.

    A point at which both _xand _yare zero is said to be a criticalpoint (xo, yo). It means that there is no flow orientation at thispoint. If there is only one critical point, it carries the position

    information about the flow pattern. The critical points of asystem with differential equations (1) and (2) are theintersection points of the following equations

    F(x, y) cx dyf 0 (4)

    G(x, y) ax by e 0 (5)

    The purpose of the phase portrait analysis is to obtain fullinsight for the fault trajectory during the steady state andtransients concerning the effect of high fault resistance andshunt susceptance.

    3 Phase portrait classification

    The elements of the characteristic matrix are used todetermine six flow patterns if the characteristic matrix isnon-singular. The two classifications may be added forsingular characteristic matrices. The eight classifications

    derived from the characteristic matrix are no-flow, constantflow, star-node, node, improper node, saddle, centre andspiral. Some of these shapes are shown in Fig. 1 and canbe obtained as follows:

    Node: the components x(t) and y(t) are exponentials thateither simultaneously converge to, or diverge from, thefixed point coordinates x0 and y0.

    Saddle: the components x(t) and y(t) are exponentials, whereas one of the components [either x(t) or y(t)]converges to the fixed point, and the other diverges from it.

    Spiral: the components x(t) and y(t) are exponentiallymodulated sinusoidal functions the resulting streamlineforms a spiral curve.

    Throughout the paper, a set of differential equations, whichdescribe the transmission line during the fault, is given. Thedifferent events of fault condition can be described bythe phase trajectory, which passes through the phaseportraits. The transition between different phase portraitsis determined since phase trajectories are continuous lines.

    The isotangent-based approach in the phase portrait willsuccessfully discriminate between the internal and external

    faults, as will be explained in Section 7.

    4 Conventional problems

    Fig. 2 shows the tripping quadrilateral characteristic in thecase of high fault resistance. It will be noted that theresistance of the fault arc takes the fault impedance outsidethe relays tripping characteristic and, hence, it does notdetect this condition. Alternatively, it is only picked upeither by zone 2 or zone 3 in which case tripping will beunacceptably delayed.

    The infeed effect can be seen from Fig. 3, where there areother lines and sources feeding current to a fault atFfrom busR. The relays at bus Sare set beyond this fault point to F0. As

    Figure 1 Phase portraits for a system of linear first-order

    differential equations

    50 IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 4956

    & The Institution of Engineering and Technology 2008 doi: 10.1049/iet-gtd:20080456

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    3/8

    theinfeed increases in proportionaltoIS, theapparentimpedanceincreases. Since this impedance, as measured by the distancerelay, is larger than the actual, the reach of the relay decreases;that is, the relay protects less of the line as infeed increases.

    The following equation proves that the voltage at bus Sis related to the current at bus Sby the equation ES ZLISZR(ISIR) and the apparent impedance seen by the relay at Sis Zapp ES=IS ZL ZR(ISIR)=IS the current IR, thecontribution to the fault from the tap, is known as the infeedcurrent and causes error in measurements.

    The paper introduces a new approach on the basis of the

    instantaneous rate of change of current against voltage usingthe phase portrait. The relay uses one end transmission linedata taking in consideration all fault events such as highfault resistance and remote-end infeed. The new approachdiscriminates precisely between internal and external faults onthe basis of the isotangent lines.

    5 Transmission line model

    The technique takes a more fundamental approach withregard to the inclusion of high-frequency oscillations thatcan occur during faults. This can be achieved by including

    the effect of capacitance of the transmission line using asingle PI-section transmission model.

    Fig. 4 shows the protection equivalent for an EHVline. The length of the transmission and the effect of theline susceptance will be variables in this case. Thetransmission line is divided into two sections, separated at

    the point of fault application. Both segments of thetransmission line are represented as a long line equivalentPI-section.

    6 Development of the algorithm

    As explained above, a set of differential equations to describethe transmission line with the susceptance effect during thefault will be given. Fig. 4 shows the transmission lineconfiguration with the shunt capacitance and the faultresistance. It is assumed that the magnitude of any arc

    resistance Rf is small and that the effect of the branchedcapacitance can be neglected. Therefore the circuit equationcan be described as follows

    v nR(is ic) nLd(is ic)

    dt Rf(is ic ir) (6)

    Rearranging (6)

    d(is ic)

    dt

    v

    nL

    R

    L

    RfnL

    (is ic)

    RfnL

    (ir) (7)

    The current ic is the capacitive current flowing through thelumped equivalent capacitance at relay location (S), and itis related simply to the voltage at the relaying point bynoting that

    ic C

    n

    dv

    dt(8)

    Rearranging (8)

    dv

    dtn

    C(is ic)

    n

    C(is) (9)

    Figure 3 Effect of the infeed on impedance measured by

    distance relays

    Figure 2 Under reach of the distance relays (RF , fault

    impedance; RS , Z1, quadrilateral characteristics and ZL line

    impedance)

    Figure 4 Protection equivalent for an EHV line (n fault

    location on protected line as portion of the total)

    IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 49 56 51

    doi: 10.1049/iet-gtd:20080456 & The Institution of Engineering and Technology 2008

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    4/8

    From (7) and (9) in the matrix form

    d(is ic)dt

    dv

    dt

    0

    B@

    1

    CA

    R

    L

    RfnL

    1

    nLn

    C 0

    0

    B@

    1

    CA

    (is ic)

    v

    0

    B@

    1

    CA

    RfnL

    ir

    n

    Cis

    0BB@

    1CCA

    (10)

    Equation (10) describes the characteristic of the transmissionline. To study the dynamics of the transmission line, Rf, L, Rand n are changed. For every given Rf, the correspondingphase portrait is ic v is shown in Fig. 5. Also, for everygiven n, the corresponding phase portrait is ic v is

    shown in Fig. 6. The phase portrait of the external andinternal faults is seen in Fig. 7.

    As shown in Figs. 5 7, the instantaneous current against voltage of (10) gives valuable information on the transientand stability characteristics at different fault conditions. Afocus on the transient part is circled to show that valuableinformation and features can be obtained from the phaseportrait.

    From the recognisable shapes in the portrait, some fault

    events can be easily identified. However, other fault eventshave valuable information, but the phase portrait is not easilyidentified. For this reason, another phase portrait producedfrom is ic v is used to assist the discrimination betweendifferent fault events, specifically the internal and externalfaults. The isotangent-based approach in a phase portraitto discriminate between the internal and external faults isused.

    7 Isotangent-based approach

    An isotangent curve is the locus of points in a flow field,which share the same flow orientation [26, 27]. All curvesare always straight lines for non-singular, two-dimensionalfirst-order phase portraits. Also, the relationship betweenthe flow orientation and the slope of the isotangent lines isunique for every fault condition along the transmission

    lines. The following section can prove these facts.

    Since the systems of (9) and (7) corresponding to (2) and(1), respectively, are linear in the parameter space, parameterasset is unique for any given flow pattern. For a

    Figure 5 Phase portrait of is ic v for different R (R 0

    and R 100) and SLG fault on phase-A

    Figure 6 Phase portrait of is ic v for different fault

    location and SLG fault on phase-A

    Figure 7 Phase portrait illustration between internal and

    external faults and SLG fault on phase-A

    52 IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 4956

    & The Institution of Engineering and Technology 2008 doi: 10.1049/iet-gtd:20080456

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    5/8

    two-dimensional first-order system, the tangent is given by

    tan(w(x, y)) tan(w(is ic, v)) _y

    _x

    _v_is _ic

    G(x, y)

    F(x, y)

    ax by e

    cx

    dy

    f

    (11)

    if tan w b (11) becomes

    (a cb)x (b db)y (efb) 0 (12)

    which is a straight line for any given w. Since the lineequation is a function ofb, every straight line is unique fora given point set in a given flow pattern. For the matrix ~Agiven in (10), there is a unique straight line with a uniqueisotangent, which depends on the system parameters.

    The paper presents a phase portrait on the basis of the

    instantaneous rate of change of current against voltage, thatis, (d(is ic))=dt dv=dt rather than the instantaneouscurrent and voltage (is ic v). The phase portrait givesrecognisable different isotangents for different cases ofinternal faults.

    The X/R ratio of the transmission line is fixed. The(d(is ic))=dt dv=dt portrait can be considered as a locusof the instantaneous _Z. The portrait gives more visualisationof the instantaneous part of the _Z. The portrait in such acase is concerned with the rate of change of the impedancein the transient part, which is highly affected by the

    PX

    and PR of the transmission line configuration at a point(i.e. equivalent of X/equivalent of R at the point of fault).

    This ratio can be calculated from the input impedance at thefault point as will be explained later.

    In case of external faults, the phase portrait gives an almostfixed isotangent. Thus, the process of discrimination betweeninternal and external faults can be easily identified.

    The isotangent given in (11) depends on theP

    X=P

    Rratio (i.e. equivalent of X/equivalent of R at the point offault) of system configuration. This can be easily calculatedas the input impedance at the fault point. The isotangent

    or theP

    X=P

    R ratio is changed during internal faultsbased on the fault location. During external faults, theisotangent or the

    PX=

    PR ratio seen at the fault point is

    almost fixed. Consequently, the isotangent in the phaseportrait during internal faults is different from theisotangent in case of external faults.

    To prove this fact, a simplified configuration given inFig. 8 is considered. For a three-phase transmission lineoperating at a certain voltage and having an impedanceconnected to the generating station busbar, for a certainfault at Fi, the input impedance is calculated using

    Thevenins as

    ZinFi (n(RjX)jXg1)==((1 n)(RjX)jXg2) (13)

    Equation (13) has theP

    X=P

    R ratio determined as

    Xt(n(1 n)R2X0X00)Rt(2n(1 n)RXX

    000)

    Rt(n(1 n)R2 X0X00)Xt(2n(1 n)RXX

    000)(14)

    where Rt nR (1 n)R, Xt nXXg1 (1 n)XXg2, X

    0 (nXXg1), X

    00

    ((1 n)XXg2) and X000

    2n(1 n)RX (1 n)RXg1 (1 n)RXg2. TheP

    X=P

    R

    ratio given above becomes totally different if the faultlocation is changed along the transmission line, and thatcan be easily recognised using the phase portrait based onthe rate of change of the current against voltage.

    Similarly, the input impedance Zin for an external fault atFe can be calculated as follows

    ZinFe (RjX) jXg1)==( jXg2) (15)

    Equation (15) has theP

    X=P

    R ratio determined as

    R2

    (XXg1)(XXg1 Xg2)RXg2

    (16)

    TheP

    X=P

    R ratio given in (16) is different from thePX=

    PR ratio given in (14). For extra and ultra high-

    voltage transmission lines, the installed source capacitieslocated at the sending and receiving ends of the TL arealmost constant. For an internal fault along the transmissionline, the value of the (14) changes according to the faultlocation and, consequently, the isotangent is also changed.

    However, for an external fault along the transmission line,

    the source capacity at the far end (Fig. 8) will have very highsource capacity compared with the transmission lineparameters. The rate of change of the impedance in thephase portrait during the transient part depends on theP

    X=P

    R ratio given in (16), which is totally differentfrom that in (14) and almost with a fixed value.Consequently, the isotangent in the phase portrait is alsodifferent and fixed for external faults.

    Figs. 911 show the instantaneous rate of change ofcurrent against voltage using a phase portrait for the samefault conditions as shown in Figs. 5 7. As shown in these

    figures, the isotangent for different cases of fault events isdifferent, and the process of discrimination is easy to beidentified.

    Figure 8 Input impedance seen at fault point for internaland external faults

    IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 49 56 53

    doi: 10.1049/iet-gtd:20080456 & The Institution of Engineering and Technology 2008

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    6/8

    The criterion for the protection relay to initiate a trip signalis set. For the studied configuration system given in Fig.8, (16)is used to set up the relay. Thus, the isotangents in the case ofinternal faults respond up or down away from the thresholdboundary THR. In this respect, the isotangent for externalfaults has a unique value equal to THR. This criterionsuccessfully recognises the internal and external faults, and100% of the transmission line is protected. THR can be

    easily calculated as given in (11) from the slop of the lineafter fault inception.

    Fig. 12 shows the difference between the conventionaldistance relay with three zones and the new techniquebased on the phase portrait principle.

    8 Power system configuration

    The power system used for testing the proposed new methodis a part of 500 kV power systems shown in Fig. 13. Thesystem includes two generating stations. A relay is locatedat Bus S as shown in Fig. 13. The voltage and current

    Figure 9 Phase portrait based on rate of change fordifferent fault resistances (R 0 and R 100)

    Figure 11 Phase portrait illustration between internal and

    external faults in case of rate of change-based portrait

    Figure 10 Phase portrait based on rate of change for

    different fault locations

    Figure 12 Zone of operation between the distance relay

    and rate of change using phase portrait

    a Distance relayb The rate of change using phase portrait

    Figure 13 Studied configuration system and different

    events of faults

    54 IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 4956

    & The Institution of Engineering and Technology 2008 doi: 10.1049/iet-gtd:20080456

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    7/8

    signals are the inputs to the relays, and 250 km is the linelength. The results are described on a phase portrait using

    the (d(is ic))=dt dv=dt diagram. The relay is set toprotect 100% of the line (250 km); it forms the first zoneof the relay. The power system is modelled and differentsymmetrical and unsymmetrical faults with solid and faultresistance are simulated using the ATP-EMTP.

    The performance of the proposed technique was evaluatedfor different types of internal and external faults, faultlocation and fault resistance. Results showed faults aretaken with a fault resistance ranged from 0 to 300 V.

    Table 1 shows the different cases of the fault conditions. The THR is adjusted using (16). The parameters of the

    studied configuration system is given in the Appendix.

    Fig. 14 shows the phase portrait of the rate of change of thecurrent against voltage for different cases of fault conditionsmentioned in ( Table 1). The isotangent in the case of anexternal fault is almost fixed. The suggested technique givesthe solutions for the symmetrical and unsymmetrical faults

    with solid and fault resistances. The threshold boundary

    THR for 100% of the transmission line is set. As shown inthe figure, all internal faults rather than the external faultshave different isotangents (i.e. the isotangent in the case ofinternal faults lies up and down away from the THRboundary). This indicates that the relay is able to distinguishsuccessfully between internal and external faults.

    9 Conclusions

    The paper presented a technique for discrimination betweeninternal and external faults using a phase portrait. The paperused the geometrical theory of differential equations to drive a

    symbol set on the basis of the visual appearance of the phaseportraits. The phase portrait of the instantaneous rate ofchange of current against voltage gives valuable informationon the transient and stability characteristics of the systemconfiguration with high fault resistance. The paper showedthat for non-singular, two-dimensional first-order phaseportraits, every isotangent curve is a unique straight line. Theresults showed the flow patterns for different internal faults

    with different isotangent lines. In case of external faults, aunique isotangent line is given. The discrimination betweeninternal and external faults can be identified. Furthermore, theimpact of high fault resistance and remote-end feed isassessed. The results indicated different events with types offault condition, fault resistance and fault location using thephase trajectory which passes through the phase portrait.Hundred percent of the protected zone is successfullydiscriminated rather than the external zone.

    10 References

    [1] HOROWITZ S.H., PHADKE A.G.: Power system relaying

    (Research Studies Press, Taunton, Somerset, England, 1992)

    [2] SIDHU T.S., GHOTRA D.S., SACHDEV M.S.: An adaptive distance

    relay and its performance comparison with a fixed data

    window distance relay, IEEE Trans. Power Deliv., 2002,17, (3), pp. 691697

    [3] SIDHU T.S., GHOTRA D.S., SACHDEV M.S.: A fast distance relay

    using adaptive data window filters. Power Engineering

    Society Summer Meeting, 2000, 1620 July 2000, vol. 3,

    pp. 14071412

    [4] KUDO H., SASAKI H., SEO K., TAKAHASHI M., YOSHIDA K., MAEDA T.:

    Implementation of a digital distance relay using an

    interpolated integral solution of a differential equation,

    IEEE Trans. Power Deliv., 1988, 3, (4), pp. 14751484

    [5] SACHDEV M.S., SIDHU T.S., GHOTRA D.S.: Implementation of an

    adaptive data window technique in a distance relay.

    Figure 14 Phase portrait discrimination between internal

    and external faults at different fault events given in Table 1

    Table 1 Different cases of fault conditions

    Case Fault resistance,V Fault location, km

    F1 (SLG) 200 50

    F2 (SLG) 0 150

    F3 (SLG) 100 150

    F4 (SLG) 0 240

    F5 (SLG) 0 10 (external)

    F6 (SLG) 50 10 (external)

    F7 (2LG) 0 50 (external)

    F8 (SLG) 200 50 (external)

    IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 49 56 55

    doi: 10.1049/iet-gtd:20080456 & The Institution of Engineering and Technology 2008

    www.ietdl.org

  • 8/8/2019 New Principle for Transmission Line Protection Using Phase Portrait Plane

    8/8

    Developments in Power System Protection, 2001, 7th Int.

    Conf. (IEE), 912 April 2001, pp. 161164

    [6] ALTUVE FERRER H.J., DIAZ VERDUZCO I., VAZQUEZ MARTINEZ E.:

    Fourier and Walsh digital filtering algorithms for distance

    protection, IEEE Trans. Power Syst., 1996, 11, ( 1) ,

    pp. 457462

    [7] WAIKA R D.L ., L IE W A .C ., EL ANGOVA N S.: Design,

    implementation and performance evaluation of a new

    digital distance relaying algorithm, IEEE Trans. Power Syst.,

    1996, 11, (1), pp. 448 456

    [8] EISSA M.M., MASOUD M.: A novel digital distance relaying

    technique for transmission line protection, IEEE Trans.

    Power Deliv., 2001, 16, (3), pp. 380384

    [9] ZHNAG N., MLADEN K.: Transmission line boundary

    protection using wavelet transform and neurasl network,IEEE Trans. Power Deliv., 2007, 22, (2), pp. 859869

    [10] VAZQUEZ E., CASTRUITA J., CHACON O., CONDE A.: A new

    approach traveling wave distance protection-part I:

    algorithm, IEEE Trans. Power Deliv., 2007, 22, ( 2) ,

    pp. 795800

    [11] WAIKAR D.L., ELANGOVAN S., LIEW A.C.: Fault impedance

    estimation algorithm for digital distance relaying, IEEE

    Trans. Power Deliv., 1994, 9, (3), pp. 13751383

    [12] DJURI M., RADOJEVI Z., TERZIJA V.: Time domain solution of

    fault distance estimation and arcing faults detection on

    overhead lines, IEEE Trans. Power Deliv., 1999, 14, (1),

    pp. 6067

    [13] TERZIJA V., RADOJEVIC Z.: Numerical algorithm for adaptive

    autoreclosure and protection of medium voltage overhead

    lines, IEEE Trans. Power Deliv., 2004, 19, (2), pp. 554559

    [14] XIA Y.O., LI K.K., DAVID A.K.: Adaptive relay setting for stand-

    alone digital distance protection, IEEE Trans. Power Deliv.,

    1994, 9, (1), pp. 480491

    [15] EREZZAGHI M.E., CROSSLEY P.A.: The effect of high resistancefaults on a distance relay. Power Engineering Society

    General Meeting, 2003, 1317 July 2003, vol. 4

    [16] LI Y., DE-SHU C., XIANG-GEN Y., ZHANG Z.: The uniform formula

    of single phase earth fault distance relay with

    compensation. Transmission and Distribution Conference

    and Exposition, 2003, IEEE PES, 712 September 2003,

    vol. 1, pp. 245249

    [17] SIDHU T.D., BALTAZAR D.S., PALOMINO R.M., SACHDEV M.S.: A new

    approach for calculating zone-2 setting of distance relays

    and its use in an adaptive protection system, IEEE Trans.Power Deliv., 2004, 19, (1), pp. 7077

    [18] LEE J.B., JUNG C.H., KIM I.D., BAEK Y.K.: Protective relay testing

    and characteristic analysis for high impedance faults in

    transmission lines. Power Engineering Society Summer

    Meeting, 1999, 1822 July 1999, vol. 2, pp. 10761081

    [19] CHUL-HWAN K., JEONG-YONG H., AGGARWAL R.K.: An enhanced

    zone 3 algorithm of a distance relay using transientcomponents and state diagram, IEEE Trans. Power Deliv.,

    2005, 20, (1), pp. 3946

    [20] BO Z.Q., AGGARWAL R.K., JOHNS A.T.: A new approach to

    transmission protection using fault generated high

    frequency current signals. Proc. 12th Power Systems

    Computation Conf., Dresden, Germany, 19 23 August 1996

    [21] JOHNS A.T., AGARAWAL A.P.: New approach to power line

    protection based upon the detection of fault induced

    high frequency signals, Proc. Inst. Elect. Eng. Gen. Trans.

    Distrib., 1990, 137, (4), pp. 307313

    [22] PATHIRANA V., MCLAREN P.G.: A hybrid algorithm for high

    speed transmission line protection, IEEE Trans. Power

    Deliv., 2005, 20, (4), pp. 24222428

    [23] W YL IE C .R ., B AR RET L .C.: Advanced engineering

    mathematics (McGraw-Hill, 1995, 6th edn.)

    [24] RAO A.R., JAIN R.C.: A taxonomy for texture description

    and identification (Springer-Verlag, 1990)

    [25] CHO K.H., KIM J.R., BAKE S., CHOI H.S., CHOO S.M.: Inferring

    bimolecular regulatory networks from phase portraits of time

    series expression profiles, FEBS Lett., 2006, 580, pp. 35113518

    [26] ROUSSEUW P.J.: Leat median of squares regression,

    J. Am. Stat. Assoc., 1984, 79, (388), pp. 871880

    [27] RAO R., SCHUNK B.: Computing oriented texture fields.

    Proc. Conf., Computer Vision and Pattern Recognition, 1989

    11 Appendix

    The parameters of the three-phase 500 kV transmission lineR0 0:1896V=km

    R1 0:018V=km

    L0 3:456mH=km

    L1 0:936mH=km

    C0 0:00828mF=km

    C1 0:01134mF=km

    The source capacity rating is

    GVA 10 (sending end)

    GVA 6 (receiving end)

    56 IET Gener. Transm. Distrib., 2009, Vol. 3, No. 1, pp. 4956

    & The Institution of Engineering and Technology 2008 doi: 10.1049/iet-gtd:20080456

    www.ietdl.org