CN101825677B - 基于环流量电流相角差的同塔四回输电线路故障选相方法 - Google Patents

Abstract

本发明公开了一种基于环流量电流相角差的同塔四回输电线路故障选相方法，包括以下步骤：(1)数据处理：计算线路测量端的各回线各相的故障电流工频相量；(2)解耦处理：求出同塔四回输电线路的四回线解耦变换矩阵，对各回线各相的故障电流工频相量进行解耦变换，求出同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量；(3)计算环流量电流的相角差；(4)故障选线：根据不同回线故障时环流量故障电流的相角特征，建立不同回线故障的判据，判断故障回线；(5)计算序电流的相角差；(6)故障选相：建立单回线不同相故障的判据，判别故障相。本发明方法不受系统阻抗的影响，具有准确性高的优点。

Description

基于环流量电流相角差的同塔四回输电线路故障选相方法

技术领域

本发明涉及同塔四回输电线路故障选相方法，具体涉及一种基于环流量电流相角差的同塔四回输电线路故障选相方法。 

背景技术

随着我国经济协调快速发展，作为国民经济命脉的电力工业正处于大发展阶段，电网规模得到快速扩张，电网结构、输电方式等方面都发生了巨大的变化。其中，由于我国特别是华东、华北、华南等沿海地区人口密集，民房建筑多，可用耕地日益减少，拆迁赔偿费越来越高，输电线路走廊通道资源受到制约的矛盾越发突出。因此，为了提高线路单位走廊的输电容量和土地利用率，降低电力建设成本，同塔多回输电方式已成为我国电网建设的必然趋势。 

输电线路担负着传送电能的重要任务，是电力系统的命脉，国内外都发生过由于输电线路故障而诱发的电力系统瓦解的灾难性大事故。可见，输电线路故障直接威胁到电力系统的安全可靠运行，关系到国民经济能否稳定快速发展。所以，如果输电线路发生故障，必须快速准确地实现故障的识别与切除。而故障选相则是高压输电线路继电保护的一个关键环节，其快速准确地选相对于确保输电线路安全和电力系统稳定有着积极意义。此外，故障选相在常规保护装置中为综合重合闸服务。但是，由于同塔四回线路导线间存在非常复杂的耦合互感，现有的基于12序分量的故障选相原理，受系统阻抗的影响，当系统发生振荡，尤其发生快速振荡时，选相元件可能误判。因此亟需不受系统阻抗影响的同塔四回输电线路故障选相方法。 

发明内容

本发明的目的在于克服现有技术的缺点和不足，提供了一种基于环流量电流相角差的同塔四回输电线路故障选相方法，该方法不受系统阻抗的影响，具有准确性高的优点。 

本发明的目的通过下述技术方案实现：一种基于环流量电流相角差的同塔四回输电线路故障选相方法，包括以下步骤： 

(1)数据处理：根据已知的同塔四回输电线路测量端的各回线各相电流瞬时值，采用傅氏算法计算线路测量端的各回线各相的故障电流工频相量；其中四回线分别用I、II、III、IV表示，各回线三相分别用A、B、C表示； 

(2)解耦处理：求出同塔四回输电线路的四回线解耦变换矩阵，并对各回线各相的故障电流工频相量进行解耦变换，求出同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量，其中E表示同向分量，环流分量分别用F、G、H表示，1表示正序分量，2表示负序分量，0表示零序分量；提取解耦变换后的正、负序环流分量； 

(3)计算环流量电流的相角差：根据步骤(2)解耦变换后的正、负序环流分量，分别计算F模正序与G模正序故障电流之间的相角差、F模负序与G模负序故障电流之间的相角差、G模正序与H模正序故障电流之间的相角差、G模负序与H模负序故障电流之间的相角差； 

(4)故障选线：由步骤(3)计算得到的环流量电流的相角差，在步骤(2)解耦处理的基础上，根据不同回线故障时环流量故障电流的相角特征，建立不同回线故障的判据，判断故障回线； 

(5)计算序电流的相角差：由步骤(2)得到的正、负序环流分量，分别计算F模正序与F模负序故障电流之间的相角差、G模正序与G模负序故障电流之间的相角差、H模正序与H模负序故障电流之间的相角差； 

(6)故障选相：由步骤(4)判别的故障回线，和步骤(5)计算得到的序电流相角差，基于对称分量法和单回线故障的正、负序电流相角特征，建立单回线不同相故障的判据，判别故障相，从而确定故障回线与故障相。 

为更好的实现本发明，所述步骤(1)数据处理，具体包括以下步骤： 

S1.1采集同塔四回输电线路测量端的各回线三相故障电流瞬时值，i_mn(t)表示t时刻m回线n相的电流瞬时值，其中m∈(I，II，III，IV)，n∈(A，B，C)； 

S1.2利用傅氏算法分别计算步骤S1.1中各回线三相的故障电流工频相量， 

表示m回线n相的故障电流工频相量，其中m∈(I，II，III，IV)，n∈(A，B，C)。 

所述步骤S1.2利用傅氏算法分别计算步骤S1.1中各回线三相的故障电流工频相量，具体是指： 

设故障时刻为0时刻，取k＝1，2，3，...，N_C，其中N_C为一工频周期采样点数，由采样频率f_s决定，即N_C＝T_Cf_S，T_C为工频周期，采样间隔ΔT＝1/f_S，对任意回线任意相电流的瞬时值都有： 

${\stackrel{\·}{I}}_{\mathrm{mn}}=\frac{2}{N_C}\underset{k=1}{\overset{N_C}{\mathrm{\Σ}}}[i_{\mathrm{mn}}\left(\mathrm{k\ΔT}\right)\·e^{-\mathrm{j\ωk\ΔT}}]$

其中e为自然常数，j为虚数符号，ω为工频角频率，ω＝2πf_S，i_mn(kΔT)表示m回线n相电流在kΔT时刻的瞬时值， 

表示m回线n相电流的工频相量，m∈(I，II，III，IV)，n∈(A，B，C)。 

所述步骤(2)解耦处理，具体包括以下步骤： 

S2.1求出同塔四回输电线路的四回线解耦变换矩阵： 

根据任意多个不平衡的相量系统分解为平衡相量系统的变换形式，将同塔四回输电线路的各回线三相看成是一组导线，则同塔四回输电线路等效为四个不平衡的相量系统，有如下矩阵式表示： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_I\\ {\stackrel{\·}{I}}_{\mathrm{II}}\\ {\stackrel{\·}{I}}_{\mathrm{III}}\\ {\stackrel{\·}{I}}_{\mathrm{IV}}\end{array}\right]=\left[\begin{array}{cccc}1& 1& 1& 1\\ 1& -j& -1& j\\ 1& -1& 1& -1\\ 1& j& -1& -j\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_E\\ {\stackrel{\·}{I}}_F\\ {\stackrel{\·}{I}}_G\\ {\stackrel{\·}{I}}_H\end{array}\right]$

即： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_E\\ {\stackrel{\·}{I}}_F\\ {\stackrel{\·}{I}}_G\\ {\stackrel{\·}{I}}_H\end{array}\right]=\frac{1}{4}\left[\begin{array}{cccc}1& 1& 1& 1\\ 1& j& -1& -j\\ 1& -1& 1& -1\\ 1& -j& -1& j\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_I\\ {\stackrel{\·}{I}}_{\mathrm{II}}\\ {\stackrel{\·}{I}}_{\mathrm{III}}\\ {\stackrel{\·}{I}}_{\mathrm{IV}}\end{array}\right]$

式中， 

分别对应表示I、II、III、IV回线的电流相量， 

表示同向模量， 分别对应表示为F、G、H环流模量，由于I、II、III、IV回线的电流相量为三相系统，同向模量和F、G、H环流模量也为三相系统，则上式应写成： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EP}}\\ {\stackrel{\·}{I}}_{\mathrm{FP}}\\ {\stackrel{\·}{I}}_{\mathrm{GP}}\\ {\stackrel{\·}{I}}_{\mathrm{HP}}\end{array}\right]=\frac{1}{4}\left[\begin{array}{cccc}{E}_{3\×3}& E_{3\×3}& E_{3\×3}& E_{3\×3}\\ E_{3\×3}& j{E}_{3\×3}& -E_{3\×3}& -j{E}_{3\×3}\\ E_{3\×3}& -E_{3\×3}& E_{3\×3}& -E_{3\×3}\\ E_{3\×3}& -j{E}_{3\×3}& -E_{3\×3}& j{E}_{3\×3}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{IP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IVP}}\end{array}\right]$

式中，E_3×3为3×3的单位矩阵， 

分别对应表示为I、II、III、IV回线的三相电流相量，即 

其中 

表示I回线A相的电流工频相量， 

表示I回线B相的电流工频相量， 

表示I回线C相的电流工频相量， 

表示II回线A相的电流工频相量， 

表示II回线B相的电流工频相量， 

表示II回线C相的电流工频相量， 

表示III回线A相的电流工频相量， 

表示III回线B相的电流工频相量， 

表示III回线C 相的电流工频相量， 表示IV回线A相的电流工频相量， 

表示IV回线B相的电流工频相量， 表示IV回线C相的电流工频相量； 

表示为同流量的三相电流相量，即 

分别对应表示为F、G、H环流量的三相电流相量，即 

其中 

表示同向E模量A相的电流工频相量， 

表示同向E模量B相的电流工频相量， 

表示同向E模量C相的电流工频相量， 

表示环流F模量A相的电流工频相量， 

表示环流F模量B相的电流工频相量， 

表示环流F模量C相的电流工频相量， 

表示环流G模量A相的电流工频相量， 

表示环流G模量B相的电流工频相量， 

表示环流G模量C相的电流工频相量， 

表示环流H模量A相的电流工频相量， 

表示环流H模量B相的电流工频相量， 

表示环流H模量C相的电流工频相量； 

对同向模量和环流模量的三相系统进行解耦，采用对称分量变换，则有： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EM}}\\ {\stackrel{\·}{I}}_{\mathrm{FM}}\\ {\stackrel{\·}{I}}_{\mathrm{GM}}\\ {\stackrel{\·}{I}}_{\mathrm{HM}}\end{array}\right]=\begin{array}{c}\left[\begin{array}{c}A{\stackrel{\·}{I}}_{\mathrm{EP}}\\ A{\stackrel{\·}{I}}_{\mathrm{FP}}\\ A{\stackrel{\·}{I}}_{\mathrm{GP}}\\ A{\stackrel{\·}{I}}_{\mathrm{HP}}\end{array}\right]\end{array}=\frac{1}{4}\left[\begin{array}{cccc}A& A& A& A\\ A& \mathrm{jA}& -A& -\mathrm{jA}\\ A& -A& A& -A\\ A& -\mathrm{jA}& -A& \mathrm{jA}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{IP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IVP}}\end{array}\right],\left(A=\frac{1}{3}\left[\begin{array}{ccc}1& a& a^2\\ 1& a^2& a\\ 1& 1& 1\end{array}\right]\right)$

式中，A为对称分量变换矩阵， 

其中e为自然常数，j为虚数符号， 

表示为同向E模量的序分量电流，即 

分别对应表示为环流F、G、H模量的序分量电流相量，即 

其中 

表示同向E模量正序电流相量， 

表示同向E模量负序电流相量， 

表示同向E模量零序电流相量， 表示环流F模量正序电流， 

表示环流F模量负序电流， 

表示环流F模量零序电流， 

表示环流G模量正序电流， 

表示环流G模量负序电流， 

示环流G模量零序电流， 

表示环流H模量正序电流， 

表示环流H模量负序电流， 

示环流H模量零序电流； 

则有同塔四回输电线路的解耦变换式： 

设M为同塔四回输电线路的解耦变换矩阵，则： 

S2.2对线路测量端的各回线各相故障电流工频相量进行解耦变换： 

${\stackrel{\·}{I}}_M=M{\stackrel{\·}{I}}_P$

其中， 

为同塔四回输电线路测量端流入线路的四回线各相故障电流工频相量， 

为同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量， 

${\stackrel{\·}{I}}_M={[{\stackrel{\·}{I}}_{E1},{\stackrel{\·}{I}}_{E2},{\stackrel{\·}{I}}_{E0},{\stackrel{\·}{I}}_{F1},{\stackrel{\·}{I}}_{F2},{\stackrel{\·}{I}}_{F0},{\stackrel{\·}{I}}_{G1},{\stackrel{\·}{I}}_{G2},{\stackrel{\·}{I}}_{G0},{\stackrel{\·}{I}}_{H1},{\stackrel{\·}{I}}_{H2},{\stackrel{\·}{I}}_{H0}]}^T;$

S2.3提取解耦变换后的正、负序环流分量，即提取 

和 

所述步骤(3)计算环流量电流的相角差，具体是指： 

根据步骤S2.3提取的解耦变换后的正、负序环流分量 

和 

分别计算 

与 相角差 

与 

的相角差 

与 

的相角差 与 

的相角差 

其中，Angle()表示相量的相角函数。 

所述步骤(4)故障选线，具体是指： 

S4.1分析不同回线故障的故障特征： 

将解耦变换矩阵M写为： 

$M=\left[\begin{array}{cccc}{A}^{\′}& A^{\′}& A^{\′}& A^{\′}\\ A^{\′}& {\mathrm{jA}}^{\′}& {-A}^{\′}& {-\mathrm{jA}}^{\′}\\ A^{\′}& -A^{\′}& A^{\′}& {-A}^{\′}\\ A^{\′}& {-\mathrm{jA}}^{\′}& {-A}^{\′}& {\mathrm{jA}}^{\′}\end{array}\right],$ 其中 $A^{\′}=\frac{1}{12}\left[\begin{array}{ccc}1& a& a^2\\ 1& a^2& a\\ 1& 1& 1\end{array}\right]$

代入S2.2中对线路测量端的各回线各相故障电流工频相量进行解耦变换公式，得解耦变换公式2： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EM}}\\ {\stackrel{\·}{I}}_{\mathrm{FM}}\\ {\stackrel{\·}{I}}_{\mathrm{GM}}\\ {\stackrel{\·}{I}}_{\mathrm{HM}}\end{array}\right]=\left[\begin{array}{cccc}{A}^{\′}& A^{\′}& A^{\′}& A^{\′}\\ A^{\′}& j{A}^{\′}& -A^{\′}& -j{A}^{\′}\\ A^{\′}& -A^{\′}& A^{\′}& -A^{\′}\\ A^{\′}& -j{A}^{\′}& -A^{\′}& j{A}^{\′}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{IP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IVP}}\end{array}\right]$

根据同塔四回输电线路发生故障时，不同回线故障的不同故障边界条件，建立不同回线故障的判据： 

(a)若I回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(b)若II回线发生任意类型故障，则有： 

代入解耦变换公式2得： ${\stackrel{\·}{I}}_{\mathrm{EM}}=A^{\′}{\stackrel{\·}{I}}_{\mathrm{IIP}},{\stackrel{\·}{I}}_{\mathrm{FM}}=j{A}^{\′}{\stackrel{\·}{I}}_{\mathrm{IIP}},{\stackrel{\·}{I}}_{\mathrm{GM}}=-A^{\′}{\stackrel{\·}{I}}_{\mathrm{IIP}},{\stackrel{\·}{I}}_{\mathrm{HM}}=-j{A}^{\′}{\stackrel{\·}{I}}_{\mathrm{IIP}},$

则有： 

(C)若III回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(D)若IV回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

${\stackrel{\·}{I}}_{\mathrm{EM}}=A^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},$ ${\stackrel{\·}{I}}_{\mathrm{FM}}=-j{A}^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},$ ${\stackrel{\·}{I}}_{\mathrm{HM}}=-A^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},$ ${\stackrel{\·}{I}}_{\mathrm{GM}}=j{A}^{\′}{\stackrel{\·}{i}}_{\mathrm{IVP}},$ 则有： 

S4.2根据同塔四回输电线路环流量故障特征，构建故障回线判据，设不同回线故障判据的相位裕度为 

则得下列故障选线的判据： 

A)当 

则为I回线故障； 

B)当 

则为II回线故障； 

C)当 

则为III回线故障； 

D)当 

则为IV回线故障。 

所述步骤(5)计算序电流的相角差，具体为： 

在判断故障回线后，对故障回线进行故障相判别，分别计算环流量电流 

与 

的相角差 

与 

的相角差 

与 

的相角差 

其中，Angle()表示相量的相角函数。 

所述步骤(6)故障选相，具体为： 

基于现有的对称分量法分析单回线各种故障类型的正、负序分量的故障特征，建立单回线不同相故障的判据，设单回线不同相故障判据的相位裕度为 

则有故障选相的判据为： 

A)当 

则为A相故障； 

B)当 

则为B相故障； 

C)当 则为C相故障； 

D)当 则为AB相故障； 

E)当 

则为BC相接地故障； 

F)当 

则为AC相接地故障； 

根据步骤(5)计算出的环流量电流 与 

的相角差 

与 的相角差 

与 

的相角差 

及上述单回线不同相故障的判据，判别故障相。 

优选的，所述相位裕度 

的取值范围从0°至45°。 

优选的，所述单回线不同相故障判据的相位裕度 

取值范围从0°至30°。 

本发明的工作原理：本发明一种基于环流量电流相角差的同塔四回输电线路故障选相方法，是利用同塔四回输电线路解耦变换和同塔四回输电线路环流量实现故障回线判别的方法。由于同塔四回输电线路存在复杂的线间互感和相间互感，需要对其进行解耦，通过建立同塔四回输电线路解耦变换矩阵，对同 塔四回输电线路故障电流进行解耦变换，将12个相分量变换为12个序分量，当输电线路处于故障状态下，各序分量将包含不同的故障特性，分析不同回线故障时故障电流环流量的相角特征，构建同塔四回输电线路故障选线判据，判断故障回线，然后利用对称分量法和单回线不同相故障的正、负序电流的相角特征，构建同塔四回输电线路故障选相判据，从而得到一种基于环流量电流相角差的同塔四回输电线路故障选相方法。 

与现有技术相比，本发明具有以下有益效果： 

第一、可靠性高：现有的同塔四回输电线路故障选相方法，在系统发生振荡尤其是快速振荡时，可能发生误选，可靠性较低；由于本发明中线路两端电压的环流分量为零，环流量故障等效序网不含线路外系统，即与线路外系统无关，因此本发明方法不受系统阻抗的影响，在系统发生振荡避免了选相元件发生误选的可能，提高了系统的稳定性及故障选相的可靠性； 

第二、运算量少：本发明方法基于所提的同塔四回输电线路解耦矩阵，仅需提取故障环流量即可实现利用电流相角差的故障选相，运算量少，易于实现快速的故障选相元件； 

第三、满足自动重合闸的需要：本发明方法利用环流模量和正、负序分量相互间的相角关系实现的电流选相原理，完全能满足自动重合闸的需要； 

第四、实用性强：本发明方法原理简单，确定故障类型和故障相别时裕度角较大，可以根据具体的设备情况及需要的精度，确定故障类型和故障相别时裕度角，有很高的选相精度和选相可靠性。 

附图说明

图1是本发明一种同塔四回输电系统的示意图； 

图2是本发明本发明基于环流量电流相角差的同塔四回输电线路故障选相方法的流程图。 

具体实施方式

下面结合实施例及附图，对本发明作进一步地详细说明，但本发明的实施方式不限于此。 

实施例 

如图1所示，一种同塔四回输电系统，设输电线路左端为测量端，以测量端流入输电线路的方向为各回线电流的正方向，四回线分别用I、II、III、IV表 示，各回线三相分别用A、B、C表示，如：IA表示I回线A相、IIB表示II回线B相。该系统应用本发明基于环流量电流相角差的同塔四回输电线路故障选相方法，如图2所示，包括以下步骤： 

S1数据处理：根据已知的同塔四回输电线路测量端的各回线各相电流瞬时值，采用傅氏算法计算线路测量端的各回线各相的故障电流工频相量； 

步骤S1具体为： 

S1.1采集同塔四回输电线路测量端的各回线三相故障电流瞬时值，i_mn(t)表示t时刻m回线n相的电流瞬时值，其中m∈(I，II，III，IV)，n∈(A，B，C)；即i_IA(t)，i_AB(t)，i_IC(t)，i_IIA(t)，i_IIB(t)，i_IIC(t)，i_IIIA(t)，i_IIIB(t)，i_IIIC(t)，i_IVA(t)，i_IVB(t)，i_IVC(t)，其中i_IA(t)表示t时刻I回线A相的电流瞬时值，i_IB(t)表示t时刻I回线B相的电流瞬时值，i_IC(t)表示t时刻I回线C相的电流瞬时值，i_IIA(t)表示t时刻II回线A相的电流瞬时值，i_IIB(t)表示t时刻II回线B相的电流瞬时值，i_IIC(t)表示t时刻II回线C相的电流瞬时值，i_IIIA(t)表示t时刻III回线A相的电流瞬时值，i_IIIB(t)表示t时刻III回线B相的电流瞬时值，i_IIIC(t)表示t时刻III回线C相的电流瞬时值，i_IVA(t)表示t时刻IV回线A相的电流瞬时值，i_IVB(t)表示t时刻IV回线B相的电流瞬时值，i_IVC(t)表示t时刻IV回线C相的电流瞬时值； 

S1.2利用傅氏算法分别计算步骤S1.1中各回线三相的故障电流工频相量，设故障时刻为0时刻，取k＝1，2，3，...，N_C，其中N_C为一工频周期采样点数，由采样频率f_s决定，即N_C＝T_Cf_S，T_C为工频周期20ms，ΔT＝1/.f_S，对任意回线任意相电流的瞬时值都有： 

${\stackrel{\·}{I}}_{\mathrm{mn}}=\frac{2}{N_C}\underset{k=1}{\overset{N_C}{\mathrm{\Σ}}}[i_{\mathrm{mn}}\left(\mathrm{k\ΔT}\right)\·e^{-\mathrm{j\ωk\ΔT}}]$

其中e为自然常数，j为虚数符号，ω为工频角频率，ω＝2πf_S，i_mn(kΔT)表示m回线n相电流在kΔT时刻的瞬时值， 

表示m回线n相电流的工频相量，m∈(I，II，III，IV)，n∈(A，B，C)，其中 

表示I回线A相的故障电流工频相量， 

表示I回线B相的故障电流工频相量， 

表示I回线C相的故障电流工频相量， 

表示II回线A相的故障电流工频相量， 

表示II回线B相的故障电流工频相量， 

表示II回线C相的故障电流工频相量， 

表示III回线A相的故障电流工频相量， 

表示III回线B相的故障电流工频相量， 表示III回线C相的故障电流工频相量， 

表示IV回线A相的故障电流工频相量， 

表示IV回线B相的故障电流工频相量， 表示IV回线C相的故障电流工频相量。 

S2解耦处理：求出同塔四回输电线路的四回线解耦变换矩阵，并对各回线各相的故障电流工频相量进行解耦变换，求出同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量，其中E表示同向分量，环流分量分别用F、G、H表示，1表示正序分量，2表示负序分量，0表示零序分量；提取解耦变换后的正、负序环流分量； 

步骤S2具体为： 

S2.1求出同塔四回输电线路的四回线解耦变换矩阵： 

根据C.L.Forstescue在1918年提出的任意多个不平衡相量系统分解为平衡相量系统的变换形式，在此基础上，对同塔四回输电线路的各回线三相看成是一组导线(线路数取4)，则同塔四回输电线路可等效为四个不平衡的相量系统，可有如下矩阵式表示： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_I\\ {\stackrel{\·}{I}}_{\mathrm{II}}\\ {\stackrel{\·}{I}}_{\mathrm{III}}\\ {\stackrel{\·}{I}}_{\mathrm{IV}}\end{array}\right]=\left[\begin{array}{cccc}1& 1& 1& 1\\ 1& -j& -1& j\\ 1& -1& 1& -1\\ 1& j& -1& -j\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_E\\ {\stackrel{\·}{I}}_F\\ {\stackrel{\·}{I}}_G\\ {\stackrel{\·}{I}}_H\end{array}\right]$

可得： $\left[\begin{array}{c}{\stackrel{\·}{I}}_E\\ {\stackrel{\·}{I}}_F\\ {\stackrel{\·}{I}}_G\\ {\stackrel{\·}{I}}_H\end{array}\right]=\frac{1}{4}\left[\begin{array}{cccc}1& 1& 1& 1\\ 1& j& -1& -j\\ 1& -1& 1& -1\\ 1& -j& -1& j\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_I\\ {\stackrel{\·}{I}}_{\mathrm{II}}\\ {\stackrel{\·}{I}}_{\mathrm{III}}\\ {\stackrel{\·}{I}}_{\mathrm{IV}}\end{array}\right]$

式中， 分别对应表示I、II、III、IV回线的电流相量， 

表示同向模量， 

分别对应表示为F、G、H环流模量，由于I、II、III、IV回线的电流相量为三相系统，同向模量和F、G、H环流模量也为三相系统，则上式应写成： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EP}}\\ {\stackrel{\·}{I}}_{\mathrm{FP}}\\ {\stackrel{\·}{I}}_{\mathrm{GP}}\\ {\stackrel{\·}{I}}_{\mathrm{HP}}\end{array}\right]=\frac{1}{4}\left[\begin{array}{cccc}{E}_{3\×3}& E_{3\×3}& E_{3\×3}& E_{3\×3}\\ E_{3\×3}& {\mathrm{jE}}_{3\×3}& -E_{3\×3}& -j{E}_{3\×3}\\ E_{3\×3}& -E_{3\×3}& E_{3\×3}& {-E}_{3\×3}\\ E_{3\×3}& -j{E}_{3\×3}& -E_{3\×3}& j{E}_{3\×3}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{IP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IVP}}\end{array}\right]$

式中，E_3×3为3×3的单位矩阵， 

分别对应表示为I、II、 III、IV回线的三相电流相量，即 

${\stackrel{\·}{I}}_{\mathrm{IIIP}}={[{\stackrel{\·}{I}}_{\mathrm{IIIA}},{\stackrel{\·}{I}}_{\mathrm{IIIB}},{\stackrel{\·}{I}}_{\mathrm{IIIC}}]}^T,$ ${\stackrel{\·}{I}}_{\mathrm{IVP}}={[{\stackrel{\·}{I}}_{\mathrm{IVA}},{\stackrel{\·}{I}}_{\mathrm{IVB}},{\stackrel{\·}{I}}_{\mathrm{IVC}}]}^T;$

表示为同流量的三相电流相量，即 

分别对应表示为F、G、H环流量的三相电流相量，即 

其中 

表示同向E模量A相的电流工频相量， 

表示同向E模量B相的电流工频相量， 

表示同向E模量C相的电流工频相量， 

表示环流F模量A相的电流工频相量， 

表示环流F模量B相的电流工频相量， 

表示环流F模量C相的电流工频相量， 

表示环流G模量A相的电流工频相量， 

表示环流G模量B相的电流工频相量， 

表示环流G模量C相的电流工频相量， 

表示环流H模量A相的电流工频相量， 表示环流H模量B相的电流工频相量， 

表示环流H模量C相的电流工频相量； 

对同向模量和环流模量的三相系统进行解耦，采用对称分量变换，则有： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EM}}\\ {\stackrel{\·}{I}}_{\mathrm{FM}}\\ {\stackrel{\·}{I}}_{\mathrm{GM}}\\ {\stackrel{\·}{I}}_{\mathrm{HM}}\end{array}\right]=\begin{array}{c}\left[\begin{array}{c}A{\stackrel{\·}{I}}_{\mathrm{EP}}\\ A{\stackrel{\·}{I}}_{\mathrm{FP}}\\ A{\stackrel{\·}{I}}_{\mathrm{GP}}\\ A{\stackrel{\·}{I}}_{\mathrm{HP}}\end{array}\right]\end{array}=\frac{1}{4}\left[\begin{array}{cccc}A& A& A& A\\ A& \mathrm{jA}& -A& -\mathrm{jA}\\ A& -A& A& -A\\ A& -\mathrm{jA}& -A& \mathrm{jA}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{Y}}_{\mathrm{IP}}\\ {\stackrel{\·}{Y}}_{\mathrm{IIP}}\\ {\stackrel{\·}{Y}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{Y}}_{\mathrm{IVP}}\end{array}\right],$ $(A=\frac{1}{3}\left[\begin{array}{ccc}1& a& a^2\\ 1& a^2& a\\ 1& 1& 1\end{array}\right])$

式中，A为对称分量变换矩阵， 其中e为自然常数，j为虚数符号， 

表示为同向E模量的序分量电流，即 

分别对应表示为环流F、G、H模量的序分量电流相量，即 

其中 

表示同向E模量正序电流相量， 

表示同向E模量负序电流相量， 

表示同向E模量零序电流相量， 

表示环流F模量正序电流， 

表示环流F模量负序电流， 

表示环流F模量零序电流， 

表示环流G模量正序电流， 

表示环流G模量负序电流， 示环流G模量零序电流， 

表示环流H模量正序电流， 

表示环流H模量负序电流， 

示环流H模量零序电流； 

则有同塔四回输电线路的解耦变换式： 

设M为同塔四回输电线路的解耦变换矩阵，则： 

其中 其中e为自然常数，j为虚数符号；经过严谨的推导，矩阵M的系数应该为1/12；但实际计算中，一般不需要乘上系数1/12也具有比较高的精度。 

S2.2对线路测量端的各回线各相故障电流工频相量进行解耦变换： 

${\stackrel{\·}{I}}_M=M{\stackrel{\·}{I}}_P$

其中， 

为同塔四回输电线路测量端流入线路的四回线各相故障电流工频相量， ${\stackrel{\·}{I}}_P={[{\stackrel{\·}{I}}_{\mathrm{IA}},{\stackrel{\·}{I}}_{\mathrm{IB}},{\stackrel{\·}{I}}_{\mathrm{IC}},{\stackrel{\·}{I}}_{\mathrm{IIA}},{\stackrel{\·}{I}}_{\mathrm{IIB}},{\stackrel{\·}{I}}_{\mathrm{IIC}},{\stackrel{\·}{I}}_{\mathrm{IIIA}},{\stackrel{\·}{I}}_{\mathrm{IIIB}},{\stackrel{\·}{I}}_{\mathrm{IIIC}},{\stackrel{\·}{I}}_{\mathrm{IVA}},{\stackrel{\·}{I}}_{\mathrm{IVB}},{\stackrel{\·}{I}}_{\mathrm{IVC}}]}^T;$

为同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量， ${\stackrel{\·}{I}}_M={[{\stackrel{\·}{I}}_{E1},{\stackrel{\·}{I}}_{E2},{\stackrel{\·}{I}}_{E0},{\stackrel{\·}{I}}_{F1},{\stackrel{\·}{I}}_{F2},{\stackrel{\·}{I}}_{F0},{\stackrel{\·}{I}}_{G1},{\stackrel{\·}{I}}_{G2},{\stackrel{\·}{I}}_{G0},{\stackrel{\·}{I}}_{H1},{\stackrel{\·}{I}}_{H2},{\stackrel{\·}{I}}_{H0}]}^T;$

S2.3提取解耦变换后的正、负序环流分量，即提取 

和 

S3计算环流量电流的相角差：根据步骤(2)解耦变换后的正、负序环流量故障电流，分别计算F模正序与G模正序故障电流之间的相角差、F模负序与G模负序故障电流之间的相角差、G模正序与H模正序故障电流之间的相角差、G模负序与H模负序故障电流之间的相角差；具体是指： 

根据步骤S2.3提取的解耦变换后的正、负序环流分量 

和 

分别计算 与 

相角差 

与 

的相角差 

与 

的相角差 

与 

的相角差 

其中，Angle()表示相量的相角函数。 

S4故障选线：由步骤(3)计算得到的环流量电流的相角差，在步骤(2)解耦处理的基础上，根据不同回线故障时环流量故障电流的相角特征，建立不同回线故障的判据，判断故障回线；步骤S4具体包括： 

S4.1分析不同回线故障的故障特征： 

将解耦变换矩阵M写为： 

$M=\left[\begin{array}{cccc}{A}^{\′}& A^{\′}& A^{\′}& A^{\′}\\ A^{\′}& {\mathrm{jA}}^{\′}& {-A}^{\′}& {-\mathrm{jA}}^{\′}\\ A^{\′}& -A^{\′}& A^{\′}& {-A}^{\′}\\ A^{\′}& {-\mathrm{jA}}^{\′}& {-A}^{\′}& {\mathrm{jA}}^{\′}\end{array}\right],$ 其中 $A^{\′}=\frac{1}{12}\left[\begin{array}{ccc}1& a& a^2\\ 1& a^2& a\\ 1& 1& 1\end{array}\right]$

代入S2.2中对线路测量端的各回线各相故障电流工频相量进行解耦变换公式， 得解耦变换公式2： 

$\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{EM}}\\ {\stackrel{\·}{I}}_{\mathrm{FM}}\\ {\stackrel{\·}{I}}_{\mathrm{GM}}\\ {\stackrel{\·}{I}}_{\mathrm{HM}}\end{array}\right]=\left[\begin{array}{cccc}{A}^{\′}& A^{\′}& A^{\′}& A^{\′}\\ A^{\′}& j{A}^{\′}& -A^{\′}& -j{A}^{\′}\\ A^{\′}& -A^{\′}& A^{\′}& -A^{\′}\\ A^{\′}& -j{A}^{\′}& -A^{\′}& j{A}^{\′}\end{array}\right]\left[\begin{array}{c}{\stackrel{\·}{I}}_{\mathrm{IP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IIIP}}\\ {\stackrel{\·}{I}}_{\mathrm{IVP}}\end{array}\right]$

根据同塔四回输电线路发生故障时，不同回线故障的不同故障边界条件，建立不同回线故障的判据： 

(a)若I回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(b)若II回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(C)若III回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(D)若IV回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

${\stackrel{\·}{I}}_{\mathrm{EM}}=A^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},{\stackrel{\·}{I}}_{\mathrm{FM}}=-j{A}^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},{\stackrel{\·}{I}}_{\mathrm{HM}}=-A^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},{\stackrel{\·}{I}}_{\mathrm{GM}}=j{A}^{\′}{\stackrel{\·}{I}}_{\mathrm{IVP}},$ 则有： 

表1为利用本发明提出的解耦变换矩阵，对同塔四回线各回线故障进行分析，得出同塔四回输电线路环流量故障特征，并构建故障回线判据。 

表1各回线故障的环流量特征 

S4.2根据同塔四回输电线路环流量故障特征，构建故障回线判据，设不同回线故障判据的相位裕度为 

所述相位裕度 

的取值范围从0°至45°，本实施例中 

取值为45°，则得下列故障选线的判据： 

A)当 

则为I回线故障； 

B)当 

则为II回线故障； 

C)当 则为III回线故障； 

D)当 

则为IV回线故障。 

S5计算序电流的相角差：由步骤(2)得到的环流量故障电流，分别计算F模正序与F模负序故障电流之间的相角差、G模正序与G模负序故障电流之间的相角差、H模正序与H模负序故障电流之间的相角差；步骤S5具体为： 

在判断故障回线后，对故障回线进行故障相判别，分别计算环流量电流 

与 

的相角差 

与 的相角差 

与 的相角差 

其中，Angle()表示相量的相角函数。 

S6故障选相：由步骤(4)判别的故障回线，和步骤(5)计算得到的序电流相角差，基于对称分量法和单回线故障的正、负序电流相角特征，建立单回线不同相故障的判据，判别故障相，从而确定故障回线与故障相。 

步骤S6具体为： 

基于现有的对称分量法分析单回线各种故障类型的正、负序分量的故障特征，得表2所示的单回线各相故障时的故障特征： 

表2各相故障的环流量特征 

设单回线不同相故障判据的相位裕度为 

所述单回线不同相故障判据的相位裕度 

取值范围从0°至30°，本实施例中 

取值为30°，则有故障选相的判据为： 

A)当 则为A相故障； 

B)当 

则为B相故障； 

C)当 则为C相故障； 

D)当 

则为AB相故障； 

E)当 

则为BC相接地故障； 

F)当 

则为AC相接地故障； 

根据步骤(5)计算出的环流量电流 

与 

的相角差 

与 的相角差 

与 

的相角差 

及上述单回线不同相故障的判据，判别故障相。 

采用电磁暂态仿真程序ATP/EMTP，构建如图1所示的同塔四回输电系统仿真模型，通过全面的故障仿真分析计算，验证本发明所提的故障选相方法。表3列出了不同故障情况下，环流量电流相角差值和故障选相结果，结果表明，本发明所提的故障选相方法是准确的。 

表3环流量电流相角差值和故障选相结果(度) 

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受所述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。 

Claims (10)

1.一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，包括以下步骤：

(1)数据处理：根据已知的同塔四回输电线路测量端的各回线各相电流瞬时值，采用傅氏算法计算线路测量端的各回线各相的故障电流工频相量；其中四回线分别用I、II、III、IV表示，各回线三相分别用A、B、C表示；

(2)解耦处理：求出同塔四回输电线路的四回线解耦变换矩阵，并对各回线各相的故障电流工频相量进行解耦变换，求出同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量，其中E表示同向分量，环流分量分别用F、G、H表示，1表示正序分量，2表示负序分量，0表示零序分量；提取解耦变换后的正、负序环流分量；

(3)计算环流量电流的相角差：根据步骤(2)解耦变换后的正、负序环流分量，分别计算F模正序与G模正序故障电流之间的相角差、F模负序与G模负序故障电流之间的相角差、G模正序与H模正序故障电流之间的相角差、G模负序与H模负序故障电流之间的相角差；

(4)故障选线：由步骤(3)计算得到的环流量电流的相角差，在步骤(2)解耦处理的基础上，根据不同回线故障时环流量故障电流的相角特征，建立不同回线故障的判据，判断故障回线；

(5)计算序电流的相角差：由步骤(2)得到的正、负序环流分量，分别计算F模正序与F模负序故障电流之间的相角差、G模正序与G模负序故障电流之间的相角差、H模正序与H模负序故障电流之间的相角差；

(6)故障选相：由步骤(4)判别的故障回线，和步骤(5)计算得到的序电流相角差，基于对称分量法和单回线故障的正、负序电流相角特征，建立单回线不同相故障的判据，判别故障相，从而确定故障回线与故障相。

2.根据权利要求1所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述步骤(1)数据处理，具体包括以下步骤：

S1.1采集同塔四回输电线路测量端的各回线三相故障电流瞬时值，i_mn(t)表示t时刻m回线n相的电流瞬时值，其中m∈(I，II，III，IV)，n∈(A，B，C)；

S1.2利用傅氏算法分别计算步骤S1.1中各回线三相的故障电流工频相量， 

表示m回线n相的故障电流工频相量，其中m∈(I，II，III，IV)，n∈(A，B，C)。

3.根据权利要求2所述一种基于环流量电流相角差的同塔四回输电线路故 障选相方法，其特征在于，所述步骤S1.2利用傅氏算法分别计算步骤S1.1中各回线三相的故障电流工频相量，具体是指：

设故障时刻为0时刻，取k＝1，2，3，...，N_C，其中N_C为一工频周期采样点数，由采样频率f_s决定，即N_C＝T_Cf_s，T_C为工频周期，采样间隔ΔT＝1/f_S对任意回线任意相电流的瞬时值都有：

其中e为自然常数，j为虚数符号，ω为工频角频率，ω＝2πf_S，i_mn(kΔT)表示m回线n相电流在kΔT时刻的瞬时值， 

表示m回线n相电流的工频相量，m∈(I，II，III，IV )，n∈(A，B，C)。

4.根据权利要求2所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述步骤(2)解耦处理，具体包括以下步骤：

S2.1求出同塔四回输电线路的四回线解耦变换矩阵：

根据任意多个不平衡的相量系统分解为平衡相量系统的变换形式，将同塔四回输电线路的各回线三相看成是一组导线，则同塔四回输电线路等效为四个不平衡的相量系统，有如下矩阵式表示：

即：

式中， 

分别对应表示I、II、III、IV回线的电流相量， 表示同向模量， 

分别对应表示为F、G、H环流模量，由于I、II、III、IV回线的电流相量为三相系统，同向模量和F、G、H环流模量也为三相系统，则上式应写成： 

式中，E_3×3为3×3的单位矩阵， 

分别对应表示为I、II、III、IV回线的三相电流相量，即 

其中 

表示I回线A相的电流工频相量， 

表示I回线B相的电流工频相量， 表示I回线C相的电流工频相量， 

表示II回线A相的电流工频相量， 

表示II回线B相的电流工频相量， 

表示II回线C相的电流工频相量， 

表示III回线A相的电流工频相量， 

表示III回线B相的电流工频相量， 表示III回线C相的电流工频相量， 

表示IV回线A相的电流工频相量， 表示IV回线B相的电流工频相量， 

表示IV回线C相的电流工频相量；

表示为同流量的三相电流相量，即 

分别对应表示为F、G、H环流量的三相电流相量，即 

其中 

表示同向E模量A相的电流工频相量， 

表示同向E模量B相的电流工频相量， 

表示同向E模量C相的电流工频相量， 

表示环流F模量A相的电流工频相量， 

表示环流F模量B相的电流工频相量， 

表示环流F模量C相的电流工频相量， 

表示环流G模量A相的电流工频相量， 

表示环流G模量B相的电流工频相量， 

表示环流G模量C相的电流工频相量， 

表示环流H模量A相的电流工频相量， 

表示环流H模量B相的电流工频相量， 

表示环流H模量C相的电流工频相量；

对同向模量和环流模量的三相系统进行解耦，采用对称分量变换，则有：

式中，A为对称分量变换矩阵，a＝e^j120°，其中e为自然常数，j为虚数符号， 

表示为同向E模量的序分量电流，即 

分 别对应表示为环流F、G、H模量的序分量电流相量，即 

其中 

表示同向E模量正序电流相量， 

表示同向E模量负序电流相量， 

表示同向E模量零序电流相量， 表示环流F模量正序电流， 

表示环流F模量负序电流， 

表示环流F模量零序电流， 

表示环流G模量正序电流， 

表示环流G模量负序电流， 

示环流G模量零序电流， 

表示环流H模量正序电流， 

表示环流H模量负序电流， 

示环流H模量零序电流；

则有同塔四回输电线路的解耦变换式：

设M为同塔四回输电线路的解耦变换矩阵，则： 

S2.2对线路测量端的各回线各相故障电流工频相量进行解耦变换：

其中， 为同塔四回输电线路测量端流入线路的四回线各相故障电流工频相量， 

为同塔四回输电线路测量端流入线路的四模量各序故障电流工频相量， 

S2.3提取解耦变换后的正、负序环流分量，即提取 

和 

5.根据权利要求4所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述步骤(3)计算环流量电流的相角差，具体是指：根据步骤S2.3提取的解耦变换后的正、负序环流分量 

和 

分别计算 

与 

相角差 

与 

的相角差 

与 

的相角差 

与 的相角差 

其中，Angle()表示相量的相角函数。

6.根据权利要求5所述一种基于环流量电流相角差的同塔四回输电线路故 障选相方法，其特征在于，所述步骤(4)故障选线，具体是指：

S4.1分析不同回线故障的故障特征：

将解耦变换矩阵M写为：

其中

代入S2.2中对线路测量端的各回线各相故障电流工频相量进行解耦变换公式，得解耦变换公式2：

根据同塔四回输电线路发生故障时，不同回线故障的不同故障边界条件，建立不同回线故障的判据：

(a)若I回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有：

(b)若II回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有： 

(C)若III回线发生任意类型故障，则有： 

代入解耦变换公式2得： 

则有：

(D)若IV回线发生任意类型故障，则有： 

代入解耦变换公式2得：

则有：

S4.2根据同塔四回输电线路环流量故障特征，构建故障回线判据，设不同回线故障判据的相位裕度为 

则得下列故障选线的判据： 

A)当 

则为I回线故障；

B)当 

则为II回线故障；

C)当 则为III回线故障；

D)当 

则为IV回线故障。

7.根据权利要求6所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述相位裕度 的取值范围从0°至45°。

8.根据权利要求6所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述步骤(5)计算序电流的相角差，具体为：

在判断故障回线后，对故障回线进行故障相判别，分别计算环流量电流 

与 

的相角差 

与 

的相角差 

与 的相角差 

其中，Angle()表示相量的相角函数。

9.根据权利要求8所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述步骤(6)故障选相，具体为：

基于现有的对称分量法分析单回线各种故障类型的正、负序分量的故障特征，建立单回线不同相故障的判据，设单回线不同相故障判据的相位裕度为 

则有故障选相的判据为： 

A)当 

则为A相故障；

B)当 

则为B相故障；

C)当 则为C相故障；

D)当 

则为AB相故障；

E)当 

则为BC相接地故障；

F)当 

则为AC相接地故障；

根据步骤(5)计算出的环流量电流 

与 的相角差 

与 

的相角差 

与 的相角差 

及上述单回线不同相故障的判据，判别故障相。

10.根据权利要求9所述一种基于环流量电流相角差的同塔四回输电线路故障选相方法，其特征在于，所述单回线不同相故障判据的相位裕度 取值范围从0°至30°。 

Images