基于瞬时电流特征的电流互感器饱和识别改进方法
【专利摘要】本发明公开一种基于瞬时电流特征的电流互感器(CT)饱和改进识别方法，在传统时差法配合虚拟制动电流法识别电流互感器饱和的基础上，根据饱和前的故障电流瞬时采样值，利用离散滑窗迭代傅立叶算法拟合出CT饱和前故障电流的交直流分量情况,然后利用遗传算法(GA)训练后的BP神经网络拟合故障电流与虚拟制动电流门槛值D的隐式关系，最后将实时测量的故障电流瞬时采样值输入算法程序，自适应的选取虚拟制动电流门槛值D。本发明能够有效解决传统虚拟制动电流法选取固定门槛值在严重饱和时闭锁差动保护时间过长问题，并且能够在发生区外转区内故障时加快开放差动保护。
【专利说明】
基于瞬时电流特征的电流互感器饱和识别改进方法
技术领域
[0001] 本发明专利属于电力系统继电保护领域，特别设及一种基于瞬时电流特征的电流 互感器保护识别改进方法。 技术背景
[0002] 电流差动保护原理简单，具有很高的灵敏度及速动性，已被广泛用于线路及电气 设备的主保护，其能否可靠识别出区内故障，主要取决于对应的电流互感器能否正确的传 变故障和非故障电流。当电流互感器发生铁屯、饱和时就会使二次电流发生崎变而影响差动 保护的正确动作，准确识别电流互感器饱和与非饱和状态并实际闭锁和开放差动保护，对 保证差动保护的可靠性与速动性有重要意义。
[0003] 目前饱和识别方法主要有差分法、时差法、小波变换法、谐波比法等，其中同步识 别法（即时差法)应用比较普遍。它是根据区外故障引起互感器饱和时故障发生时刻与饱和 发生时刻（即差流出现时刻)是否存在时间差，来识别是否为区外故障，若识别为区外故障 则闭锁差动保护一段时间，但如果在闭锁期间发生区外转区内故障时，该方法易产生延时 动作。为此有学者根据差动电流自身波形特征提出了一种基于虚拟制动电流采样点差动的 CT饱和识别方法。该方法能够可靠识别出互感器退出饱和后即开放差动保护，在发生转换 性故障时也能可靠开放差动保护动作，但是其开放保护口槛值选取固定，引起差动保护灵 敏度降低，在转换性故障时反应不够迅速。
【发明内容】
[0004] 本发明的目的解决传统虚拟制动电流互感器饱和识别方法口槛值设置固定导致 由区外转区内故障时差动保护开放延时，灵敏度低的问题，提出一种根据电流互感器饱和 程度自适应调整口槛值的方法。
[0005] 本发明采用的具体技术方案如下：
[0006] -种基于瞬时电流特征的电流互感器饱和识别改进方法，包括如下步骤：
[0007] 第一步，利用时差法，判定故障是否为区外故障，若为区外故障则闭锁差动保护， 否则动作与断路器跳闽，本步骤分为两个小部分。
[000引第1步，检测线路是否出现故障，记录故障时刻Ti。
[0009] 采用离散滑窗迭代傅立叶（DFT)检测法，检测故障时刻Tl, W检测电流互感器A相 电流为例，基波可W表达为式(1)所示：
[0010] Jia化〇=Aicos( 〇k〇+Bisin( 〇k〇 (I)
[ocm]当采样数据大于一个采样周期时，Ai、B河W有式(2)表示如下：
[0012]
[001引式中：i^laao表示采样数据，Nnew表示新的采样点。
[0014]通过式（1)、（2)计算可W得到瞬时基波数据，用采集的瞬时数据减去瞬时基波数 据就得到瞬时谐波数据如式(3 )所示：
[001 引 i"la(k0 = i'la(k0-Jla(k0 (3)
[0016] 正常状态，谐波含量少，在不正常状态时，谐波含量大幅度增加，通过运个现象判 断不正常状态时刻Tl,即如式(4)所示：
[0017] i"ia(k〇>io (4)
[001引式中：i质示谐波和的口槛值。
[0019] 第2步，检测差动保护继电器是否出现差流，如果有差流出现，记录差流越限时刻 T2,若T2-Ti>2ms，则闭锁差动保护，否则在差流越限的时候发出跳闽指令。
[0020] 差流计算公式如公式(5)所示：
[0021]
[0022] 差流越限的判断条件如式(6)所示：
[0023] Id>KrIr (6)
[0024] 其中：与-4|为制动电流，Kr=IdZlr为制动系数。
[0025] 第二步，提取故障电流未饱和区区数据，拟合故障电流交流分量A、直流分量B、相 角巧和常数A四个参数。
[0026] 故障电流由周期的交流分量和衰减的直流分量组成，未饱和电流采样数据Ik可用 式(7)表示如下：
[0027]
[002引其中:fsamp为电流采样频率、f为系统频率。
[0029] 因此只要计算出A、B、气A四个未知参数就可W拟合出饱和前故障电流波形。公式 (7)经过CO-S ine变换可W得到公式(8)如下所示：
[0030]
姆）
[00川用Cl代替(如OW)、C2代替(-^ sin切并且用角频率O代替如f/fsamp，将公式（8)改 写为如下形式如公式(9)所示：
[0032] Ik = Cicos?k+C2sin?k+I3+Ak (9)
[0033] 提取两个Tl到T2之间电流的连续部分，第一部分用标号ki到ki+mi表示，第二部分用 标号k剧k伸康示，由式(9)可W写出如式（10)所示矩阵方程：
[0034] MP = I (10)
[0035] 其中；
[0036] 由于矩阵M与矩阵I中的量均为已知量，由式（10)可W解出A、B、疗、A四个未知参数 如式(11)所示：
[0037] P= (1?)" Vl (11)
[0038] 第S步，先采用遗传算法优化BP神经网络参数，然后用数据样本训练BP神经网络 模型，拟合A、B、0、A四个参数与虚拟制动电流口槛值D对应的隐式关系，本步骤具体分为两 小步，在算法初始化时进行一次。
[0039] 第1步，采用遍历实验的方法，基于虚拟制动电流保护法的原理，针对所有对称故 障与不对称故障(共10种故障类型)做仿真实验，得到一系列故障电流采样波形与实际应该 选择虚拟制动电流口槛值的数据样本。例如:针对一种故障类型，在线路长度Xl到X2间隔A X 做一次仿真实验，每种电路长度一周期内间隔Ims做一次仿真实验，一共可W得到 、一V 组故障电流采样数据与实际应该选择虚拟制动电流口槛值D的数据样本，所有故 'tSX 障类型一共可^得到^^^^^^^^^组数据样本，利用步骤二所述的方法得到每组未饱和故障 Ax 电流A、B、g、A四个参数与虚拟制动电流口槛值D-一对应的数据样本。
[0040] 第2步，基于遗传算法(GA)训练BP神经元网络，拟合A、B、(6^A四个参数与虚拟制动 电流口槛值D之间的隐式关系。将^^^^^^^^^组A、B、(6^A四个参数与虚拟制动电流口槛值D Ar 对应的数据样本输入遗传算法得到最佳的初始权值和阀值，用得到的最佳初始权值和阀值 O OOfv _ V 1 来构建BP神经元网络，再用-组数据样本训练BP神经网络。 Ax
[0041] 第四步，实时提取故障电流采样数据，自适应的调整虚拟制动电流口槛值D。
[0042] 由第二步可W得到实时故障电流采样数据的A、B、巧、A四个参数，将A、B、^、A四个 参数输入到训练好的8?神经元网络，根据8?神经元网络拟合好的4、8、^\^与〇的隐式关系， 自适应地输出相对应的虚拟制动电流口槛值D。
[0043] 第五步，计算一个周期内差流值大于虚拟制动电流值的电流采样点数N，决定差动 保护开放与闭锁时刻。
[0044] 基于虚拟制动电流采样差动的CT饱和识别方法，虚拟制动电流值选取为一个周波 内第一个电流采样点的0.8倍，计算公式如式(12)所示：
[0045] Ixr = O. 8ImaxSin(2 地)(12)
[0046] 其中：0.8为可靠系数;k为一个周波的采样点数;Imax为正弦波电流的峰值。
[0047] 计算一个周期内差流值大于虚拟制动电流值的采样点数N，并与步骤四求出的虚 拟制动电流口槛值D进行比较，当N>D时开放差动保护，当N《D时闭锁差动保护。
【附图说明】
[0048] 图1是基于瞬时电流特征的CT饱和识别方法流程图
[0049] 图2是虚拟制动电流法采样点数图
[0050] 图3是线路差动电流波形图
[0051 ] 图4是GA训练BP神经元网络流程图
[0052] 图5是BP神经元输出预测结果W及误差图
[0053] 图6是差动保护状态图
【具体实施方式】
[0054] 本发明提出一种基于瞬时电流特征的电流互感器饱和识别改进方法，该方法在基 于虚拟制动电流电流互感器饱和识别方法的基础上，通过遗传算法优化BP神经元网络，拟 合故障电流A、B、^、A四个参数与虚拟制动电流口槛值D之间的隐式关系，在实际运行中自 适应的调整虚拟制动电流口槛值D，使得区外故障造成电流互感器饱和后，故障又由区外转 区内时，能更快的开放差动保护。
[0055] W500kV，30km输电线路在0.304s发生A相瞬时区外故障，在0.45s转为区内故障来 具体说明本方法，实施步骤如下：
[0056] 第一步，利用时差法，判定故障是否为区外故障，若为区外故障则闭锁差动保护， 否则动作与断路器跳闽，本步骤分为两个小部分。
[0057] 第1步，检测线路是否出现故障，记录故障时刻Ti。
[0058] 提取电流互感器一次侧A相电流采样数据，采用离散滑窗迭代傅立叶(DFT)检测法 检测故障发生时刻，A相瞬时基波可W表达为式(1)所示：
[0059] Jia化〇=Aicos( 〇k〇+Bisin( 〇k〇 (I)
[0060] 设置一个周期电流采样24个点，当采样数据大于一个周期采样点数24时，Ai、Bi可 W有式(2)表示如下：
[0061]
[0062] 式中：i/ la( i O表示采样数据，Nnew表示新的采样点。将新加入数据进入和值，最前 一个数据退出和值，Al、Bi可式(3)计算所得：
[0063]
[0064]通过上述计算可W得到瞬时基波数据，用采集的瞬时数据减去瞬时基波数据就得 到瞬时谐波数据如式(4)所示：
[00化]i"la 化 0 = i' Ia 化 〇-Jla(k〇 (4)
[0066] 设置io为正常时谐波和的10倍，记录满足式巧)的时刻即为Tl:
[0067] i" Ia 化 O^io (5)
[006引式中：io为谐波和的口槛值。
[0069] 通过步骤一判断故障发生时刻如附图2所示，在0.304s时检测出线路发生故障。
[0070] 第2步，计算差流越限时刻T2,并判断是否为区外故障。
[0071] 根据差动保护原理，计算流过差动保护继电器的差动电流Id和制动电流Ir,记录满 足式(6)的时刻T2:
[0072] Id>KrIr (6)
[0073] 其中:Kr 取 /, =|4-4|。
[0074] 当满足T2-Ti>2ms时，判定为区外故障，闭锁差动保护，否则作用与断路器跳闽。
[0075] 差流的波形如附图3所示，可W检测出在0.3104s出现差流越限，所WTs-Ti = 6.4ms > 2ms，闭锁差动保护。
[0076] 第二步:提取故障发生时刻Tl到差流越限时刻T2之间第一个周期和第S个周期的 电流采样数据，近似得到故障电流交流分量A、直流分量B、相角庐和常数入。
[0077] 故障电流由周期的交流分量和衰减的直流分量组成，未饱和电流数据采样点Ik可 用式(7)表示如下：
[007引
[0079] 其中：fsamp为电流采样频率，取为系统频率，取50化，因此只要计算出A、 B、0、A四个未知参数就可W估算出饱和前故障电流波形。公式(7)经过CO-Sine变换可W得 到公式(8)如下所示：
[0080]
位）
[00川用Cl代替COS树、C2代替(-人妃的并且用角频率O代替2村化amp,将公式(8)改 写为如下形式如公式(9)所示：
[0082] Ik = Cicos ? k+C2sin ? k+E5+入k (9)
[0083] 将提取的两部分电流采样数据第一部分用标号ki到ki+23表示，第二部分用标号kj 到kw3表示，由式(9)可W写出如式(10)所示矩阵方程：
[0084] MP = I (10)
[0085] 其中：
[0086] 由于矩阵M与矩阵I中的量均为已知量，由式（10)可W解出A、B、^^、A四个未知参数 如式(11)所示：
[0087] P=(MTm)-Vi (11)
[0088] 第S步，采用遗传算法优化BP神经网络参数，训练BP神经网络模型，拟合A、B、0、入 四个参数与虚拟制动电流口槛值D对应的隐式关系，本步骤分为两小步，具体流程如附图4 所示。
[0089] 第1步，WA相接地故障实验为例，通过调节线路长度达到改变线路参数的目的。基 于虚拟制动电流电流互感器饱和识别法，线路长度30km到IlOkm间隔5km做一次仿真实验， 每种线路长度在一个周期内间隔Ims做一次实验，将每次仿真的实验数据经过步骤二可W 得到每次实验对应的故障电流A、B、巧、A四个参数，通过观测每次仿真实验得到的虚拟制动 电流法的采样点数波形，可W得到实际应该取得的虚拟制动电流口槛值D，运样通过遍历仿 真可W得到320组A、B、口、A与D-一对应的数据样本。
[0090] 第2步，选取其中200的数据样本作为训练样本，剩余120组作为测试样本。首先将 320组样本数据输入遗传算法得出最佳初始权值和阀值，用得到的最佳初始权值和阀值构 造BP神经元网络，再将200组训练数据样本输入优化后BP神经元网络，W此训练BP神经元网 络，最后用120组测试数据样本测试训练好的BP神经元网络，分析预测数据与期望数据之间 的误差，预测数据与期望数据之间的误差如附图5所示。由于口槛值是选取只要满足比固定 口槛值小，比实际的虚拟制动采样值的口槛值大就可W满足改进要求，所W计算准确率时 只要误差不出现负值就可W认为预测准确，上述训练结果的准确率可W达到95.93%，若增 加 1-2个点的裕度，测试结果接近100%的准确率。
[0091] 第四步:提取实时故障电流采样数据，自适应的调整虚拟制动电流口槛值D。
[0092] 由第二步可W得到实时故障电流采样数据的A、B、0、A四个参数，将A、B、巧、A四个 参数输入到由第S步训练好的BP神经网络，根据BP神经网络事先拟合好的得A、B、f、A四个 参数与虚拟制动电流口槛值D的关系，可W得到此时的口槛值0 = 7。
[0093] 第五步:计算一个周期内差动电流值大于虚拟制动电流值的电流采样点数N，决定 差动保护开放与闭锁时刻。
[0094] 取一个周期电流采样点数k = 24,由公式(12)计算出虚拟制动电流值：
[00 巧]Ixr = O. 8Imaxsin(2 地)(12)
[0096] 其中：0.8为可靠系数;Imax为正弦波的峰值。
[0097] 计算出一个周期内采样电流值大于虚拟制动电流值的采样点数，虚拟制动电流采 样点数如附图2所示。
[0〇9引由附图2可见在O . 304s时刻发生区外故障，经过6.4ms的线性传变区CT开始饱和， 时差法识别为区外故障，在0.3104s时刻闭锁差动保护，在0.45s时刻转区内故障，由于传统 虚拟制动电流法口槛值固定选取16,在出现转换性故障时，存在一定的延时才能开放差动 保护，在0.4648s时刻才检测到采样电流值大于虚拟制动电流值的采样点数N超过虚拟制动 电流口槛值D，此时才能开放差动保护动作;而本发明根据故障电流情况得到此时的最低口 槛值为7，在0.4568s时刻就能测到采样电流值大于虚拟制动电流值的采样点数N超过虚拟 制动电流口槛值D，此时开放差动保护并动作，相比较固定口槛值情况本发明所提的方法提 前了 8ms开放保护，差动保护状态如附图6所示，0表示保护闭锁，1表示保护开放。
【主权项】
1. 一种基于瞬时电流特征的电流互感器饱和识别方法，其特征在于，包括W下步骤： 1) 利用时差法确定CT未饱和区，判定故障是否为区外故障，若为区外故障则闭锁差动 保护，否则动作与断路器跳闽； 2) 提取故障电流非饱和区数据，拟合故障电流交流分量A、直流分量B、相角口和常数λ四 个参数； 3) 采用遗传算法优化ΒΡ神经网络参数，然后用数据样本训练ΒΡ神经网络模型，拟合Α、 Β、0、λ四个参数与虚拟制动电流口槛值D对应的隐式关系； 4) 实时提取故障电流采样数据，自适应的调整虚拟制动电流口槛值D; 5) 计算一个周期内差流值大于虚拟制动电流值的电流采样点数Ν，决定差动保护开放 与闭锁时刻。2. 根据权利要求1所述的基于瞬时电流特征的电流互感器饱和识别方法，其特征在于， 所述的步骤1)中利用时差法确定CT未饱和区，判定故障是否为区外故障，具体为： 步骤1-1，利用采用离散滑窗迭代傅立叶检测法，检测故障时刻Ti，由式（1)、（2)计算瞬 时基波数据：其中，la(i〇表示采样数据，Nnew表示新的采样点； 由式(3)计算计算瞬时谐波数据，根据式(4)计算瞬时谐波发生突变的时刻，即为故障 发生时刻Τι; i"la 化 〇 = i'la(k〇-Jla 化 0 (3) i"ia 化 〇>io (4) 其中，io表示谐波和的口槛值； 步骤1-2,检测差动保护继电器是否出现差流，如果有差流出现，记录差流越限时刻T2， T2-T1为CT未饱和区，若T2-Ti〉2ms，则闭锁差动保护，差流计算公式W及差流越限判别式分 另帖曰式(5)、（6)所示： /.=14-41 傾 Id>KrIr (6) 其中，/,.二-人-|为制动电流，Kr=ld/Ir为制动系数。3. 根据权利要求1所述的基于瞬时电流特征的电流互感器饱和识别方法，其特征在于， 所述的步骤2)中根据步骤1)确定的CT非饱和区，提取故障电流非饱和区数据，拟合故障电 流交流分量A、直流分量B、相角巧和常数λ四个参数，具体为： 将式(7)表示的未饱和电流采样数据Ik表示为式(8)的形式：/左二（/? cos《>)cos w/(小（一/? sin 科）sin w/(小公小 λ/( @) 其中，fsamp为电流采样频率、f为系统频率；C-'i二（.4cos批)、Q二(-/Isin知、ω二2村/ fsamp ； 提取两个Τι到Τ2之间电流的连续部分，第一部分用标号ki到ki+mi表示，第二部分用标号 kj到k片表示，由式(8)写出如式(9)所示矩阵方程，由式（10)可W解出故障电流交流分量 A、直流分量B、相角和常数λ四个参数：4. 根据权利要求1所述的基于瞬时电流特征的电流互感器饱和识别方法，其特征在于， 所述的步骤3)中，通过遍历仿真得到样本数据，将样本数据输入遗传算法，求出最佳初始权 值和阀值，用求出的初始权值和阀值来构造 ΒΡ神经元网络，最后将样本数据输入优化好的 ΒΡ神经元网络。训练ΒΡ神经元网络，具体方法为： 基于虚拟制动电流法，对所有对称故障与不对称故障，每种故障类型在线路长度XI到Χ2 间隔Αχ做一次仿真实验，每种电路长度一周期内间隔1ms做一次仿真实验，一共可W得到 ^2^^>组故障电流采样数据与实际应该选择虚拟制动电流口槛值D一一对应的数据样 Δ,ν 本，用组数据样本训练GA优化参数后的ΒΡ神经元网络。 Δλ'5. 根据权利要求1所述的基于瞬时电流特征的电流互感器饱和识别方法，其特征在于， 所述的步骤4)中，将实时采样的故障电流先经过步骤2)计算出故障电流交流分量Α、直流分 量Β、相角巧和常数λ四个参数，将Α、Β、料、λ四个参数输入有步骤3)训练好的ΒΡ神经元网络， 由ΒΡ神经元网络输出对应的虚拟制动电流口槛值D。6. 根据权利要求1所述的基于瞬时电流特征的电流互感器饱和识别方法，其特征在于， 所述的步骤5)中，通过比较电流采样值大于虚拟制动电流值的采样点数Ν与虚拟制动电流 口槛值D的关系，决定差动保护开放还是闭锁，具体为：
【文档编号】H02H3/26GK106099850SQ201610399721
【公开日】2016年11月9日
【申请日】2016年6月6日
【发明人】卜京, 马迪, 刘馨月, 孙震宇, 杨阳, 周斌, 殷明慧, 谢云云, 王俊, 张梦月
【申请人】南京理工大学