CN103259250A - 基于智能变电站站域信息识别变压器和应涌流的方法 - Google Patents

Abstract

基于智能变电站站域信息识别变压器和应涌流的方法：（1）判断有无变压器处于空载停运状态；（2）判断有无空载变压器合闸；（3）获取空载合闸变压器合闸侧和其他运行变压器原副两侧的线电流，计算运行变压器两侧差流；（4）持续计算运行变压器二次谐波含量并从空载合闸时刻开始计时，当二次谐波含量大于m%且基波超过额定电流的n%时停止计时；（5）比较合闸变压器和运行变压器二次谐波相位；（6）将运行变压器差动保护启动电流定值提高到原启动值的k₁倍并计时，当满足三个条件之一时，返回为原定值，并转回（1）。本发明在和应涌流期间，发生变压器内部故障保护仍能有效出口；和应涌流结束，差动保护动作特性能迅速恢复至原有动作状态。

Description

基于智能变电站站域信息识别变压器和应涌流的方法

技术领域

本发明属于电力系统继电保护技术领域，具体涉及一种利用智能变电站站域信息识别变压器和应涌流的方法，可用于防止智能变电站中变压器差动保护因和应涌流而发生误动。

背景技术

变压器是电力系统中重要的电气设备，其安全运行影响到电力系统的安全稳定。差动保护作为变压器保护的主保护被长期使用，其优点是动作速度快、选择性好，但受到外界因素和保护原理等方面的影响，其正确动作率并不高。变压器空载合闸，引起与其相邻的运行变压器产生和应涌流，导致运行变压器差动保护误动的事故已经被多次报道。由于和应涌流产生的隐蔽性、形式的多样性和特征的复杂性，加之实际运行中电流互感器饱和的影响，传统的涌流识别判据可能会失效而导致差动保护误动。

目前，关于和应涌流的研究已经变成变压器保护研究中的热点问题之一，也有一些防止和应涌流造成差动保护误动的建议与措施被提出：将运行变压器中性点不接地（见上官帖,谌争鸣,郭军燕.和应涌流对变压器差动保护的影响及对策.华中电力,2004,17(5):50-52.）；选择性能更好的电流互感器以防止其在和应涌流期间暂态饱和，利用零序二次谐波分量构成制动判据防止差动保护误动（见束洪春,贺勋,李立新.变压器和应涌流分析.电力自动化设备,2006,26(10):7-12.）；加强二次谐波判据，并针对TA饱和问题改善比例制动特性（见郑涛,赵萍.和应涌流对差动保护的影响因素分析及防范措施.电力系统自动化,2009,33(3):74-77）；针对Y/△接线变压器，采用时差法鉴别和应涌流（见谷君,郑涛,肖仕武等.基于时差法的Y/△接线变压器和应涌流鉴别新方法[J].中国电机工程学报,2007,27(13):6-11）。上述方法对于和应涌流的识别有一定指导意义，但实际应用中仍存在问题。断开运行变中性点并不能避免和应涌流的产生；选择性能更好的电流互感器增加了运行成本；零序二次谐波分量制动判据只适用于运行变中性点接地的情况，其有效性也有待验证；加强二次谐波闭锁判据可能会造成故障情况下保护误闭锁；时差法的有效性仍需经实际运行的检验。

发明内容

本发明所要解决的技术问题，就在于提供一种基于智能变电站站域信息识别变压器和应涌流的方法，防止变压器差动保护在发生和应涌流期间误动，可以在智能变电站使用。

解决上述技术问题，本发明采用如下的技术方案：

一种基于智能变电站站域信息识别变压器和应涌流的方法，包括以下步骤：

步骤（1）利用站域信息判断站内有无变压器处于空载停运状态

检查各变压器各相电流，若一段时间t_set1内某变压器三相线电流小于I_set1，并且开关量处于断开状态，即（I_a<I_set1）∩（I_b<I_set1）∩（I_c<I_set1）∩（开关处于断开状态）∩（t<t_set1），则判定该变压器处于空载停运状态，转到步骤（2），若不满足则继续步骤（1）；

步骤（2）利用站域信息判断站内有无空载停运变压器合闸

若空载变压器某相线电流大于I_set2，即（I_a>I_set2）∪（I_b>I_set2）∪（I_c>I_set2），则判断为空载合闸，转到步骤（3），否则返回步骤（1）；

步骤（3）获取该空载合闸变压器合闸侧和运行变压器两侧电流，计算出运行变两侧差流；

步骤（4）持续计算运行变压器二次谐波含量并从空载合闸时刻开始计时，当二次谐波含量大于m%且基波超过额定电流的n%时停止计时：

若计时时间小于t_set2则转到步骤（5），计时时间达到t_set2则转到步骤（1）；

步骤（5）比较合闸变压器和运行变压器二次谐波相位

从步骤（4）结束时开始计时，若运行变压器二次谐波与空载合闸变压器二次谐波相角差□θ在θ_set1°～θ_set2°范围内且计时时间小于t_set3，即满足（θ_set1°<□θ<θ_set2°）∩（t<t_set3），则转入步骤（6）并停止计时；若运行变压器二次谐波与空载合闸变压器二次谐波相角差□θ不在θ_set1°～θ_set2°范围内则继续计时，计时时间到t_set3则转入步骤（1）；

步骤（6）将运行变压器差动保护启动电流定值提高到原启动值的k₁倍并计时，当满足以下三个条件之一时，返回为原定值，并转到步骤（1）；

返回条件为：

A．运行变压器差流大于k₂I_n则返回；

B．若计算到运行变压器差流落在原有差动保护动作区外，则继续计时，计时达到t_set4则返回；

C．若运行变压器差流先进入原有差动保护动作区内，则在出动作区时返回。

所述的步骤（1）中的I_set1取2%I_n，t_set1取5s，三相线电流计算公式为：

$I_{R1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(26\right)$

$I_{I1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(27\right)$

$I_{1,\mathrm{\&phi;}}=\sqrt{{I_{R1,\mathrm{\&phi;}}}^2+{I_{I1,\mathrm{\&phi;}}}^2}---\left(28\right)$

其中：φ可以为A，B和C；i_φ代表三相线电流的采样值；I_1,φ代表三相线电流的基波幅值；N为一个工频周期内电流采样点数。

所述的步骤（2）的I_set2取5%I_n。

所述的步骤（3）中，Y/△接法的差流计算公式为：

i_CA=i_C-i_A  (29)

i_da=i_CA-i_a  (30)

i_AB=i_A-i_B  (31)

i_db=i_AB-i_b  (32)

i_BC=i_B-i_C  (33)

i_dc=i_BC-i_c  (34)

其中：i_A、i_B和i_C为变压器Y侧各相电流采样值；i_a、i_b和i_c为变压器△侧各相电流采样值；

空载合闸变压器绕组电流Y/△折算如式（10）-（11）所示，

i_tA=i_C-i_A  (35)

i_tB=i_A-i_B  (36)

i_tC=i_B-i_C  (37)

其中：φ可以为A，B和C；i_φ代表空载合闸变压器三相绕组线电流。

所述的步骤（4）中的t_set2取3s，m%取6%，n%取20%，差流基波电流与额定电流比例n%算法如下：

$I_{R1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(38\right)$

$I_{I1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(39\right)$

$I_{1,\mathrm{d\&phi;}}=\sqrt{{I_{R1,\mathrm{d\&phi;}}}^2+{I_{I1,\mathrm{d\&phi;}}}^2}---\left(40\right)$

$n\%=\frac{I_{1,\mathrm{\&phi;}}}{I_n}---\left(41\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_1,dφ代表三相差流的基波幅值；I_n为额定电流；N为一个工频周期内电流采样点数；

运行变差流的二次谐波含量m%算法如下：

$I_{R2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(42\right)$

$I_{I2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(43\right)$

$I_{2,\mathrm{d\&phi;}}=\sqrt{{I_{R2,\mathrm{d\&phi;}}}^2+{I_{I2,\mathrm{d\&phi;}}}^2}---\left(44\right)$

$m\%=\frac{I_{2,\mathrm{d\&phi;}}}{I_{1,\mathrm{d\&phi;}}}---\left(45\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_2,dφ代表三相差流的二次谐波幅值；I_1,dφ代表根据式(13)～(16)计算得到的基波电流幅值；N为一个工频周期内电流采样点数。

所述的步骤（5）中的t_set3取15s，θ_set1°取120°，θ_set2°取240°；

运行变差流二次谐波相位θ₁为：

${\mathrm{\&theta;}}_1=\mathrm{arctg}\frac{I_{I2,\mathrm{d\&phi;}}}{I_{R2,\mathrm{d\&phi;}}}---\left(46\right)$

其中：I_I2,dφ，I_R2,dφ由式（17）（18）可得，为运行变差流的二次谐波的实部和虚部；

空载合闸变压器合闸侧电流二次谐波相位θ₂算法：

$I_{R2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(47\right)$

$I_{I2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(48\right)$

${\mathrm{\&theta;}}_2=\mathrm{arctg}\frac{I_{I2,\mathrm{t\&phi;}}}{I_{R2,\mathrm{t\&phi;}}}---\left(49\right)$

其中：φ可以为A，B和C；i_tφ代表空载合闸变压器绕组电流经后Y/△折算的采样值，如式（10）-（12）所示；N为一个工频周期内电流采样点数；

Δθ=θ₁-θ₂  (50)。

所述的步骤（6）中的k₁取2，t_set4取为15s。

有益效果：本发明利用智能变电站的站域信息判断变压器空载合闸，并启动和应涌流识别判据，通过比较空载合闸变压器和运行变压器励磁电流的二次谐波相位差来识别和应涌流。本发明具有如下技术特征：能有效识别和应涌流，所需的时间较短；变压器在和应涌流期间，如果发生内部故障，保护仍能按步骤（7）有效出口；和应涌流结束后，差动保护动作特性能迅速恢复至原有灵敏动作状态。

附图说明

图1为本发明实施例中变电站的连接示意图；

图2为调整后的差动保护制动特性，其中---------表示抬高后的启动值；

图3为本发明实施例的程序流程图。

具体实施方式

下面结合附图详细解释本发明的原理与具体实施方式。

图1为本发明实施例中变电站的连接图。某时刻6#变压器空载合闸产生励磁涌流，在2#运行变压器中将产生和应涌流现象。

本发明的基于智能变电站站域信息识别变压器和应涌流的方法实施例，包括以下步骤：

步骤（1）首先利用站域信息判断站内有无变压器处于空载停运状态

判断站内有无变压器空载停运判据为：检查各变压器三相线电流，一段时间t_set1内某变压器满足电流小于I_set1并且开关量处于断开状态，其中I_set1可取2%I_n，即满足条件（I_a<I_set1）∩（I_b<I_set1）∩（I_c<I_set1）∩（开关处于断开状态）∩（t<t_set1），则判断该变压器处于空载停运状态，时间t_set1可取5s。

三相线电流计算公式为：

$I_{R1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(51\right)$

$I_{I1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(52\right)$

$I_{1,\mathrm{\&phi;}}=\sqrt{{I_{R1,\mathrm{\&phi;}}}^2+{I_{I1,\mathrm{\&phi;}}}^2}---\left(53\right)$

其中：φ可以为A，B和C；i_φ代表三相线电流的采样值；I_1,φ代表三相线电流的基波幅值；N为一个工频周期内电流采样点数。

若满足判据则转到步骤（2），否则继续步骤（1）。

步骤（2）利用站域信息判断站内有无空载停运变压器合闸

判断变压器空载合闸判据为：空载变压器满足某相线电流大于I_set2，I_set2可取5%I_n，即（I_a>I_set2）∪（I_b>I_set2）∪（I_c>I_set2）。若满足判据则转到步骤（3），否则则转到步骤（1）。

步骤（3）获取该空载合闸变压器合闸侧和运行变压器两侧电流，计算出运行变两侧差流；

作为一个优选的实施例，当6#变压器空载合闸时，将转入步骤（3）获取该空载合闸变压器合闸侧和其他运行变压器原副两侧的相电流，并计算运行变差流。本实施例的2#变压器为Y/△接法，其差流计算公式为：

i_CA=i_C-i_A  (54)

i_da=i_CA-i_a  (55)

i_AB=i_A-i_B  (56)

i_db=i_AB-i_b  (57)

i_BC=i_B-i_C  (58)

i_dc=i_BC-i_c  (59)

其中：i_A、i_B和i_C为2#变压器Y侧各相线电流采样值；i_a、i_b和i_c为2#变压器△侧各相线电流采样值。

6#空载合闸变压器绕组电流Y/△折算如式（10）-（11）所示，

i_tA=i_C-i_A  (60)

i_tB=i_A-i_B  (61)

i_tC=i_B-i_C  (62)

其中：φ可以为A，B和C；i_φ代表6#空载合闸变压器三相绕组线电流。

转入步骤（4）。

步骤（4）持续计算运行变压器二次谐波含量并从空载合闸时刻开始计时，当二次谐波含量大于m%且基波超过额定电流的n%时停止计时

如果计时时间大于t_set2，说明没有发生和应涌流转回步骤（1）。t_set2可取3s，m%可取6%，n%可取20%。作为一个优选的实施例，2#变压器差流二次谐波含量大于6%，且基波超过额定电流的20%时转到步骤（5）并停止计时。差流基波电流与额定电流比例n%算法如下：

$I_{R1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(63\right)$

$I_{I1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(64\right)$

$I_{1,\mathrm{d\&phi;}}=\sqrt{{I_{R1,\mathrm{d\&phi;}}}^2+{I_{I1,\mathrm{d\&phi;}}}^2}---\left(65\right)$

$n\%=\frac{I_{1,\mathrm{\&phi;}}}{I_n}---\left(66\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_1,dφ代表三相差流的基波幅值；I_n为额定电流；N为一个工频周期内电流采样点数。

2#运行变差流的二次谐波含量m%算法如下：

$I_{R2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(67\right)$

$I_{I2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(68\right)$

$I_{2,\mathrm{d\&phi;}}=\sqrt{{I_{R2,\mathrm{d\&phi;}}}^2+{I_{I2,\mathrm{d\&phi;}}}^2}---\left(69\right)$

$m\%=\frac{I_{2,\mathrm{d\&phi;}}}{I_{1,\mathrm{d\&phi;}}}---\left(70\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_2,dφ代表三相差流的二次谐波幅值；I_1,dφ代表根据式(13)～(16)计算得到的基波电流幅值；N为一个工频周期内电流采样点数。

步骤（5）比较6#空载合闸变压器合闸侧电流二次谐波和2#运行变压器差流二次谐波相位

从步骤（4）结束时开始计时，若t_set3之内运行变压器二次谐波与空载合闸变压器二次谐波相角差不满足在θ_set1°～θ_set2°范围内，则转入步骤（1），其中t_set3可取15s，θ_set1°可取120°，θ_set2°可取240°。作为一个优选的实施例，计算6#空载合闸变压器合闸侧电流二次谐波与2#运行变压器差流二次谐波相位差，满足运行变压器二次谐波与空载合闸变压器二次谐波相角差□θ在120°～240°范围内且计时时间t小于15s，即（120°<□θ<240°）∩（t<15s），转入步骤（6）并停止计时。6#空载合闸变压器与2#运行变压器差流的二次谐波相位差Δθ算法。

2#运行变差流二次谐波相位θ₁为：

${\mathrm{\&theta;}}_1=\mathrm{arctg}\frac{I_{I2,\mathrm{d\&phi;}}}{I_{R2,\mathrm{d\&phi;}}}---\left(71\right)$

其中：I_I2,dφ，I_R2,dφ由式（17）（18）可得，为2#运行变差流的二次谐波的实部和虚部。

6#空载合闸变压器差流二次谐波相位θ₂算法：

$I_{R2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(72\right)$

$I_{I2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(73\right)$

${\mathrm{\&theta;}}_2=\mathrm{arctg}\frac{I_{I2,\mathrm{t\&phi;}}}{I_{R2,\mathrm{t\&phi;}}}---\left(74\right)$

其中：φ可以为A，B和C；i_tφ代表6#空载合闸变压器三相绕组线电流经后Y/△折算的采样值，如式（10）-（12）所示；N为一个工频周期内电流采样点数。

Δθ=θ₁-θ₂  (75)

步骤（6）提高运行变压器差动保护启动电流定值到原启动值的k₁倍并计时，其中k₁可取2，调整后差动保护制动特性如图2所示。当满足以下三个条件之一，则定值返回为原定值，并转到步骤（1）。返回条件为：

A．运行变压器差流大于k₂I_n则返回。

B．若计算到运行变压器差流落在原有差动保护动作区外，则继续计时，计时达到t_set4则返回，t_set4建议取为15s。

C．若运行变压器差流先进入原有差动保护动作区内，则在出动作区时返回。

其中条件A表示在此过程中6#变压器又发生了区内故障，则该变压器差动保护制动特性返回为原制动特性；条件B表示和应涌流未引起6#变压器原有差动保护误动，15s后返回原有动作特性；条件C表示和应涌流发展到一定阶段或结束后，6#变压器差流低于原本的差动保护制动特性，则6#变压器差动保护制动特性返回为原制动特性，返回步骤（1）继续进行下一次判别。

Claims (7)

1.一种基于智能变电站站域信息识别变压器和应涌流的方法，其特征是包括以下步骤：

步骤（1）利用站域信息判断站内有无变压器处于空载停运状态

检查各变压器各相电流，若一段时间t_set1内某变压器三相线电流小于I_set1，并且开关量处于断开状态，即（I_a<I_set1）∩（I_b<I_set1）∩（I_c<I_set1）∩（开关处于断开状态）∩（t<t_set1），则判定该变压器处于空载停运状态，转到步骤（2），若不满足则继续步骤（1）；

步骤（2）利用站域信息判断站内有无空载停运变压器合闸

若空载变压器某相线电流大于I_set2，即（I_a>I_set2）∪（I_b>I_set2）∪（I_c>I_set2），则判断为空载合闸，转到步骤（3），否则返回步骤（1）；

步骤（3）获取该空载合闸变压器合闸侧和运行变压器两侧电流，计算出运行变两侧差流；

步骤（4）持续计算运行变压器二次谐波含量并从空载合闸时刻开始计时，当二次谐波含量大于m%且基波超过额定电流的n%时停止计时：

若计时时间小于t_set2则转到步骤（5），计时时间达到t_set2则转到步骤（1）；

步骤（5）比较合闸变压器和运行变压器二次谐波相位

从步骤（4）结束时开始计时，若运行变压器二次谐波与空载合闸变压器二次谐波相角差□θ在θ_set1°～θ_set2°范围内且计时时间小于t_set3，即满足（θ_set1°<□θ<θ_set2°）∩（t<t_set3），则转入步骤（6）并停止计时；若运行变压器二次谐波与空载合闸变压器二次谐波相角差□θ不在θ_set1°～θ_set2°范围内则继续计时，计时时间到t_set3则转入步骤（1）；

步骤（6）将运行变压器差动保护启动电流定值提高到原启动值的k₁倍并计时，当满足以下三个条件之一时，返回为原定值，并转到步骤（1）；

返回条件为：

A．运行变压器差流大于k₂I_n则返回；

B．若计算到运行变压器差流落在原有差动保护动作区外，则继续计时，计时达到t_set4则返回；

C．若运行变压器差流先进入原有差动保护动作区内，则在出动作区时返回。

2.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：

所述的步骤（1）中的I_set1取2%I_n，t_set1取5s，三相线电流计算公式为：

$I_{R1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(1\right)$

$I_{I1,\mathrm{\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(2\right)$

$I_{1,\mathrm{\&phi;}}=\sqrt{{I_{R1,\mathrm{\&phi;}}}^2+{I_{I1,\mathrm{\&phi;}}}^2}---\left(3\right)$

其中：φ可以为A，B和C；i_φ代表三相线电流的采样值；I_1,φ代表三相线电流的基波幅值；N为一个工频周期内电流采样点数。

3.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：所述的步骤（2）的I_set2取5%I_n。

4.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：所述的步骤（3）中，Y/Δ接法的差流计算公式为：

i_CA=i_C-i_A  (4)

i_da=i_CA-i_a  (5)

i_AB=i_A-i_B  (6)

i_db=i_AB-i_b  (7)

i_BC=i_B-i_C  (8)

i_dc=i_BC-i_c  (9)

其中：i_A、i_B和i_C为变压器Y侧各相电流采样值；i_a、i_b和i_c为变压器△侧各相电流采样值；

空载合闸变压器绕组电流Y/Δ折算如式（10）-（11）所示，

i_tA=i_C-i_A  (10)

i_tB=i_A-i_B  (11)

i_tC=i_B-i_C  (12)

其中：φ可以为A，B和C；i_φ代表空载合闸变压器三相绕组线电流。

5.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：所述的步骤（4）中的t_set2取3s，m%取6%，n%取20%，差流基波电流与额定电流比例n%算法如下：

$I_{R1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(13\right)$

$I_{I1,\mathrm{d\&phi;}}=\frac{2}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{2\mathrm{\&pi;}}{N}\right)---\left(14\right)$

$I_{1,\mathrm{d\&phi;}}=\sqrt{{I_{R1,\mathrm{d\&phi;}}}^2+{I_{I1,\mathrm{d\&phi;}}}^2}---\left(15\right)$

$n\%=\frac{I_{1,\mathrm{\&phi;}}}{I_n}---\left(16\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_1,dφ代表三相差流的基波幅值；I_n为额定电流；N为一个工频周期内电流采样点数；

运行变差流的二次谐波含量m%算法如下：

$I_{R2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(17\right)$

$I_{I2,\mathrm{d\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{d\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(18\right)$

$I_{2,\mathrm{d\&phi;}}=\sqrt{{I_{R2,\mathrm{d\&phi;}}}^2+{I_{I2,\mathrm{d\&phi;}}}^2}---\left(19\right)$

$m\%=\frac{I_{2,\mathrm{d\&phi;}}}{I_{1,\mathrm{d\&phi;}}}---\left(20\right)$

其中：φ可以为A，B和C；i_dφ代表三相差流的采样值；I_2,dφ代表三相差流的二次谐波幅值；I_1,dφ代表根据式(13)～(16)计算得到的基波电流幅值；N为一个工频周期内电流采样点数。

6.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：所述的步骤（5）中的t_set3取15s，θ_set1°取120°，θ_set2°取240°；

运行变差流二次谐波相位θ₁为：

${\mathrm{\&theta;}}_1=\mathrm{arctg}\frac{I_{I2,\mathrm{d\&phi;}}}{I_{R2,\mathrm{d\&phi;}}}---\left(12\right)$

其中：I_I2,dφ，I_R2,dφ由式（17）（18）可得，为运行变差流的二次谐波的实部和虚部；

空载合闸变压器合闸侧电流二次谐波相位θ₂算法：

$I_{R2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\cos \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(22\right)$

$I_{I2,\mathrm{t\&phi;}}=\frac{4}{N}\underset{k=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}_{\mathrm{t\&phi;}}\left(k\right)\sin \left(k\frac{4\mathrm{\&pi;}}{N}\right)---\left(23\right)$

${\mathrm{\&theta;}}_2=\mathrm{arctg}\frac{I_{I2,\mathrm{t\&phi;}}}{I_{R2,\mathrm{t\&phi;}}}---\left(24\right)$

其中：φ可以为A，B和C；i_tφ代表空载合闸变压器绕组电流经后Y/Δ折算的采样值，如式（10）-（12）所示；N为一个工频周期内电流采样点数；

Δθ=θ₁-θ₂  (25)。

7.根据权利要求1所述的基于智能变电站站域信息识别变压器和应涌流的方法，其特征是：所述的步骤（6）中的k₁取2，t_set4取为15s。

Images