CN103248025B - 一种有损耗输电线路电流行波差动保护方法 - Google Patents

Abstract

本发明公开了一种有损耗输电线路电流行波差动保护方法。该方法首先采集输电线路两端保护安装处各个采样时刻的电压、电流行波分量，计及线路损耗的影响，通过预测算法计算输电线路一端的电流行波分量，根据得到的电流行波分量与采样得到的该端线路电流行波分量构成电流行波差动保护，若任一相电流行波差动保护动作，则向该相线路两端的断路器发跳闸信号。本发明方法无需通过插值运算求取各采样时刻的电流行波分量，减少了电流行波保护运算量，提高了电流行波差动保护动作速度。本发明可在1/4周期的短数据窗内完成输电线路电流行波保护功能，保护动作速度快，且计及线路损耗影响，因此适用于各种电压等级有损耗输电线路的整个故障暂态过程的行波保护。

Description

一种有损耗输电线路电流行波差动保护方法

技术领域

本发明涉及电力系统继电保护技术领域，具体地说是涉及一种有损耗输电线路电流行波差动保护方法。

背景技术

由于不受系统运行方式和电网结构影响且具有天然的选相功能，电流差动保护一直是各种电压等级输电线路的主保护。在220kV及以下电压等级输电线路中，由于输电线路沿线分布电容电流很小，分布电容对电流差动保护动作性能影响很小。然而，超/特高压交流输电线路的电压、电流传输具有明显的波过程，沿线分布电容电流很大，利用两端基频稳态电流分量的矢量和的幅值作为动作量的传统电流差动保护面临着电流差动保护启动电流大，而为了防止保护误动，提高启动设定值又会导致保护灵敏度不足，制约着传统电流差动保护在超/特高压交流输电线路上的应用。

行波差动保护在保护算法数学模型中考虑了分布电容的影响，行波差动保护原理上不受输电线路分布电容的影响，具有很高的动作性能。申请号200910034669.1发明专利“适用于串联电容补偿线路的行波差动保护方法”解决了分布电容对差动保护动作性能的影响，但对于行波传播延时为非整数倍采样间隔的情况，需要通过插值运算得到各个时间点上的电气量，对保护装置采样频率的要求非常高，因此对保护装置硬件要求很高，且每个采样时间都要进行插值运算，保护算法本身所需运算量大，无法满足保护速动性的要求。苏斌、董新洲和孙元章发表的《特高压带并联电抗器线路的行波差动保护》以及苏斌、董新洲和孙元章发表的《基于小波变换的行波差动保护》和申请号200410079501.X发明专利“行波保护中电压过零附近故障的检测方法”对于行波传播延时为非整数倍采样间隔的情况也需要通过插值运算得到各个时间点上的电气量，同样存在着运算量大的问题；需要进行小波变换，所需数据窗大，保护检测到故障发生所需时间长。

目前，许多学者已提出的输电线路行波差动保护方法对行波传播延时为非整数倍采样间隔的情况都需要进行插值运算求其各个时间点上的电气量，保护算法本身运算量大，对保护装置采样硬件要求高。其中部分输电线路行波差动保护方法甚至需要进行小波变换，所需数据窗大，延长了保护检测到故障发生的时间，无法满足继电保护对速动性的要求。

发明内容

本发明的目的在于克服已有技术存在的不足，提供一种有损耗输电线路电流行波差动保护方法。

为实现上述目的，本发明方法采用如下技术方案：

一种有损耗输电线路电流行波差动保护方法，其要点在于，包括如下依序步骤：

（1）提供一种保护装置，其利用位于输电线路两端的m变电站保护安装处和n变电站保护安装处的各采样时刻的行波电气量计算t采样时刻的m变电站保护安装处的0、α、β模电流行波分量

$i_{m0}^{\′}\left(t\right)=\frac{u_{m0}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_0\right))-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\ω\τ}}_0\right))}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}$

$-i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}+\frac{u_{n0}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_0\right))-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\ω\τ}}_0\right))}$

$i_{\mathrm{m\α}}^{\′}\left(t\right)=\frac{u_{\mathrm{m\α}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))-u_{\mathrm{m\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{\mathrm{c\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))}-\frac{i_{\mathrm{m\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}$

$-i_{\mathrm{n\α}}\left(t\right)-\frac{i_{\mathrm{n\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}+\frac{u_{\mathrm{n\α}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))-u_{\mathrm{n\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{\mathrm{c\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))}$

$i_{\mathrm{m\β}}^{\′}\left(t\right)=\frac{u_{\mathrm{m\β}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))-u_{\mathrm{m\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{\mathrm{c\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))}-\frac{i_{\mathrm{m\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}$

$-i_{\mathrm{n\β}}\left(t\right)-\frac{i_{\mathrm{n\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}+\frac{u_{\mathrm{n\β}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))-u_{\mathrm{n\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{\mathrm{c\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))}$

其中，t为采样时间；l为连接m变电站和n变电站的输电线路长度；T为基频分量的周期时间；Z_c0、Z_cα、Z_cβ分别为输电线路0、α、β模行波分量的特性阻抗；ν_c0、ν_cα、ν_cβ分别为输电线路0、α、β模行波分量的传播速度；ω为电力系统角频率；R₀、R_α、R_β分别为输电线路0、α、β模行波分量的电阻；u_m0(t)、u_mα(t)、u_mβ(t)分别为m变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；u_n0(t)、u_nα(t)、u_nβ(t)分别为n变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；i_n0(t)、i_nα(t)、i_nβ(t)分别为n变电站保护安装处的t采样时刻的0、α、β模的电流行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电流行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电流行波分量；

（2）保护装置将进行相模反变换得到t采样时刻的m变电站保护安装处的三相电流行波分量对采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相基频电流分量对t采样时刻的m变电站保护安装处的三相实测的电流行波分量i_mA(t)、i_mB(t)、i_mC(t)采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相实测的基频电流分量

（3）若成立，则保护向A相输电线路两端的断路器发送跳闸信号，跳开A相输电线路两端的断路器；其中，I_set为电流整定值；

（4）若成立，则保护向B相输电线路两端的断路器发送跳闸信号，跳开B相输电线路两端的断路器；其中，I_set为电流整定值；

（5）若成立，则保护向C相输电线路两端的断路器发送跳闸信号，跳开C相输电线路两端的断路器；其中，I_set为电流整定值。

本发明与现有技术相比较，具有以下积极成果：

本发明方法无需通过插值运算求取各采样时刻的电流行波分量，减少了电流行波保护运算量，提高了电流行波差动保护动作速度。本发明方法可在1/4周期的短数据窗内完成故障行波保护功能，保护动作速度快，且计及线路损耗影响，因此适用于各种电压等级输电线路的整个暂态故障过程的行波保护。

附图说明

图1为应用本发明方法的一种有损耗输电线路电流行波差动保护方法流程图。

具体实施方式

下面根据说明书附图对本发明的技术方案做进一步详细表述。

图1为应用本发明方法的一种有损耗输电线路电流行波差动保护方法流程图。本实施例首先采集输电线路在m变电站保护安装处的各t采样时刻的三相电压行波分量和三相电流行波分量；采集输电线路在n变电站保护安装处的各t采样时刻的三相电压行波分量和三相电流行波分量。

采用相模变换将m变电站保护安装处的各t采样时刻的三相电压行波分量、三相电流行波分量转换成各t采样时刻的0、α、β模电压行波分量u_m0(t)、u_mα(t)、u_mβ(t)和各t采样时刻的0、α、β模电流行波分量i_m0(t)、i_mα(t)、i_mβ(t)。

采用相模变换将n变电站保护安装处的各t采样时刻的三相电压行波分量、三相电流行波分量转换成各t采样时刻的0、α、β模电压行波分量u_n0(t)、u_nα(t)、u_nβ(t)和各t采样时刻的0、α、β模电流行波分量i_n0(t)、i_nα(t)、i_nβ(t)。

利用位于输电线路两端的m变电站保护安装处和n变电站保护安装处的各采样时刻的行波电气量计算t采样时刻的m变电站保护安装处的0、α、β模电流行波分量

$i_{m0}^{\′}\left(t\right)=\frac{u_{m0}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_0\right))-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\ω\τ}}_0\right))}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}$

$-i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}+\frac{u_{n0}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_0\right))-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos \left({\mathrm{\ω\τ}}_0\right))}$

$i_{\mathrm{m\α}}^{\′}\left(t\right)=\frac{u_{\mathrm{m\α}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))-u_{\mathrm{m\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{\mathrm{c\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))}-\frac{i_{\mathrm{m\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}$

$-i_{\mathrm{n\α}}\left(t\right)-\frac{i_{\mathrm{n\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}+\frac{u_{\mathrm{n\α}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))-u_{\mathrm{n\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{\mathrm{c\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right))}$

$i_{\mathrm{m\β}}^{\′}\left(t\right)=\frac{u_{\mathrm{m\β}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))-u_{\mathrm{m\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{\mathrm{c\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))}-\frac{i_{\mathrm{m\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}$

$i_{\mathrm{n\β}}^{\′}\left(t\right)-\frac{i_{\mathrm{n\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}+\frac{u_{\mathrm{n\β}}\left(t\right)(1-\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))-u_{\mathrm{n\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{\mathrm{c\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right))}$

其中，t为采样时间；l为连接m变电站和n变电站的输电线路长度；T为基频分量的周期时间；Z_c0、Z_cα、Z_cβ分别为输电线路0、α、β模行波分量的特性阻抗；ν_c0、ν_cα、ν_cβ分别为输电线路0、α、β模行波分量的传播速度；ω为电力系统角频率；R₀、R_α、R_β分别为输电线路0、α、β模行波分量的电阻；u_m0(t)、u_mα(t)、u_mβ(t)分别为m变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；u_n0(t)、u_nα(t)、u_nβ(t)分别为n变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；分别为n变电站保护安装处的t采样时刻的0、α、β模的电流行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电流行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电流行波分量。

将进行相模反变换得到t采样时刻的m变电站保护安装处的三相电流行波分量

对采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相基频电流分量

对t采样时刻的m变电站保护安装处的三相实测的电流行波分量采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相实测的基频电流分量

若成立，则保护向A相输电线路两端的断路器发送跳闸信号，跳开A相输电线路两端的断路器；其中，I_set为电流整定值。

若成立，则保护向B相输电线路两端的断路器发送跳闸信号，跳开B相输电线路两端的断路器；其中，I_set为电流整定值。

若成立，则保护向C相输电线路两端的断路器发送跳闸信号，跳开C相输电线路两端的断路器；其中，I_set为电流整定值。

本发明方法无需通过插值运算求取各采样时刻的电流行波分量，减少了电流行波保护运算量，提高了电流行波差动保护动作速度。本发明方法可在1/4周期的短数据窗内完成故障行波保护功能，保护动作速度快，且计及线路损耗影响，因此适用于各种电压等级输电线路的整个暂态故障过程的行波保护。

以上所述仅为本发明的较佳具体实施例，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种有损耗输电线路电流行波差动保护方法，其特性在于，包括如下依序步骤：

(1)提供一种保护装置，其利用位于输电线路两端的m变电站保护安装处和n变电站保护安装处的各采样时刻的行波电气量计算t采样时刻的m变电站保护安装处的0、α、β模电流行波分量

$\begin{array}{c}{i}_{m0}^{\′}\left(t\right)\frac{u_{m0}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_0\left)\right)-u_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos ({\mathrm{\ω\τ}}_0\left)\right)}-\frac{i_{m0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}\\ -i_{n0}\left(t\right)-\frac{i_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{1+\cos \left({\mathrm{\ω\τ}}_0\right)}+\frac{u_{n0}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_0\left)\right)-u_{n0}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_0\right)}{(Z_{c0}+\frac{R_0}{4})(1+\cos ({\mathrm{\ω\τ}}_0\left)\right)}\end{array}$

$\begin{array}{c}{i}_{m\mathrm{\α}}^{\′}\left(t\right)=\frac{u_{m\mathrm{\α}}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_{\mathrm{\α}}\left)\right)-u_{m\mathrm{\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{c\mathrm{\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos ({\mathrm{\ω\τ}}_{\mathrm{\α}}\left)\right)}-\frac{i_{m\mathrm{\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}\\ -i_{n\mathrm{\α}}\left(t\right)-\frac{i_{n\mathrm{\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}+\frac{u_{n\mathrm{\α}}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_{\mathrm{\α}}\left)\right)-u_{n\mathrm{\α}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\α}}\right)}{(Z_{c\mathrm{\α}}+\frac{R_{\mathrm{\α}}}{4})(1+\cos ({\mathrm{\ω\τ}}_{\mathrm{\α}}\left)\right)}\end{array}$

$\begin{array}{c}{i}_{m\mathrm{\β}}^{\′}\left(t\right)=\frac{u_{m\mathrm{\β}}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_{\mathrm{\β}}\left)\right)-u_{m\mathrm{\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{c\mathrm{\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos ({\mathrm{\ω\τ}}_{\mathrm{\β}}\left)\right)}-\frac{i_{m\mathrm{\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}\\ -i_{n\mathrm{\β}}\left(t\right)-\frac{i_{n\mathrm{\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{1+\cos \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}+\frac{u_{n\mathrm{\β}}\left(t\right)(1-\cos ({\mathrm{\ω\τ}}_{\mathrm{\β}}\left)\right)-u_{n\mathrm{\β}}(t-\frac{T}{4})\sin \left({\mathrm{\ω\τ}}_{\mathrm{\β}}\right)}{(Z_{c\mathrm{\β}}+\frac{R_{\mathrm{\β}}}{4})(1+\cos ({\mathrm{\ω\τ}}_{\mathrm{\β}}\left)\right)}\end{array}$

其中，t为采样时刻；l为连接m变电站和n变电站的输电线路长度；T为基频分量的周期时间；Z_c0、Z_cα、Z_cβ分别为输电线路0、α、β模行波分量的特性阻抗；ν_c0、ν_cα、ν_cβ分别为输电线路0、α、β模行波分量的传播速度；ω为电力系统角频率；R₀、R_α、R_β分别为输电线路0、α、β模行波分量的电阻；u_m0(t)、u_mα(t)、u_mβ(t)分别为m变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；u_n0(t)、u_nα(t)、u_nβ(t)分别为n变电站保护安装处的t采样时刻的0、α、β模的电压行波分量；i_n0(t)、i_nα(t)、i_nβ(t)分别为n变电站保护安装处的t采样时刻的0、α、β模的电流行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电压行波分量；分别为m变电站保护安装处的采样时刻的0、α、β模的电流行波分量；分别为n变电站保护安装处的采样时刻的0、α、β模的电流行波分量；

(2)保护装置将i′_m0(t)、i′_mα(t)、i′_mβ(t)进行相模反变换得到t采样时刻的m变电站保护安装处的三相电流行波分量i′_mA(t)、i′_mB(t)、i′_mC(t)；对i′_mA(t)、i′_mB(t)、i′_mC(t)采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相基频电流分量对t采样时刻的m变电站保护安装处的三相实测的电流行波分量i_mA(t)、i_mB(t)、i_mC(t)采用傅里叶算法计算m变电站保护安装处的t采样时刻的三相实测的基频电流分量

(3)若成立，则保护向A相输电线路两端的断路器发送跳闸信号，跳开A相输电线路两端的断路器；其中，I_set为电流整定值；

(4)若成立，则保护向B相输电线路两端的断路器发送跳闸信号，跳开B相输电线路两端的断路器；其中，I_set为电流整定值；

(5)若成立，则保护向C相输电线路两端的断路器发送跳闸信号，跳开C相输电线路两端的断路器；其中，I_set为电流整定值。