发明名称 自适应故障电流检测方法
摘要 一种自适应故障电流检测方法,其特征在于,所述的故障电流检测方法采用自适应滤波器对电流信号进行预测,将预测值作为参考依据,并设定门限值。将电流互感器检测到的电流信号值与自适应滤波器输出的电流信号的预测值比较,如果电流互感器检测到的实际的电流信号值与预测值之间的偏差超过预先设定的门限值,则认为产生故障。本发明适用于电力系统与电力设备等装置的电流故障诊断。
申请公布号 CN104198893B 申请公布日期 2017.03.15
申请号 CN201410496062.6 申请日期 2014.09.24
申请人 中国科学院电工研究所;江苏中天科技股份有限公司 发明人 马韬;张京业;戴少涛;赵连岐;薛弛;滕玉平;朱志芹;许熙
分类号 G01R31/08(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 关玲
主权项 一种自适应故障电流检测方法,其特征在于,所述的故障电流检测方法采用自适应滤波器对电流信号进行预测,将预测值作为参考依据,并设定门限值;将电流互感器检测到的电流信号值与自适应滤波器输出的电流信号的预测值比较,如果电流互感器检测到实际的电流信号值与所述预测值之间的偏差超过预先设定的门限值,则认为产生故障;具体步骤如下:(1)将当前t时刻采集到的电流互感器输出值x<sub>t</sub>与自适应滤波器的输出值<img file="FDA00011824392800000111.GIF" wi="41" he="59" />进行比较:<img file="FDA0001182439280000011.GIF" wi="797" he="226" />其中,t为数据采集时刻,x<sub>t</sub>为t时刻采集的电流互感器输出值,<img file="FDA0001182439280000012.GIF" wi="40" he="59" />为t时刻自适应滤波器的输出值;δ<sub>l</sub>、δ<sub>h</sub>分别为预先设定的常数,且满足δ<sub>l</sub>&lt;δ<sub>h</sub>;符号|·|表示取绝对值;上式表明三种情况:a、t时刻电流互感器输出值x<sub>t</sub>与自适应滤波器输出值<img file="FDA0001182439280000013.GIF" wi="41" he="62" />之间的差别小于设定值δ<sub>l</sub>时,认为电网电流正常;b、t时刻电流互感器输出值x<sub>t</sub>与自适应滤波器输出值<img file="FDA0001182439280000014.GIF" wi="40" he="61" />之间的差别不小于设定值δ<sub>l</sub>而小于设定值δ<sub>h</sub>,认为电网发生故障,需要立即使断路器动作;c、t时刻电流互感器输出值x<sub>t</sub>与自适应滤波器输出值<img file="FDA0001182439280000015.GIF" wi="43" he="62" />之间的差别大于等于设定值δ<sub>h</sub>,认为电流互感器输出受到外部干扰,该信号为虚假信号,可以忽略;(2)若判断属于情况b,则由微机继电保护装置立即给出故障发生信号至继电保护系统的继电器;若判断属于情况a或c,则认为未发生故障,根据采集到的电流互感器输出值计算此时自适应滤波器输出值:<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>W</mi><mi>t</mi><mi>T</mi></msubsup><msub><mi>X</mi><mi>t</mi></msub></mrow>]]></math><img file="FDA0001182439280000016.GIF" wi="242" he="63" /></maths>其中,<maths num="0002"><math><![CDATA[<mrow><msub><mi>W</mi><mi>t</mi></msub><mo>=</mo><msub><mi>W</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msubsup><mi>&Phi;</mi><mi>t</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>x</mi><mi>t</mi></msub><mo>-</mo><msubsup><mi>W</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><msub><mi>X</mi><mi>t</mi></msub><mo>)</mo></mrow><msub><mi>X</mi><mi>t</mi></msub><mo>-</mo><mi>C</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001182439280000017.GIF" wi="738" he="87" /></maths><maths num="0003"><math><![CDATA[<mrow><msubsup><mi>&Phi;</mi><mi>t</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msup><mi>&lambda;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>K</mi><mi>t</mi></msub><msubsup><mi>X</mi><mi>t</mi><mi>T</mi></msubsup><mo>)</mo></mrow><msubsup><mi>&Phi;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA0001182439280000018.GIF" wi="493" he="70" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mi>K</mi><mi>t</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Phi;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>X</mi><mi>t</mi></msub></mrow><mrow><mi>&lambda;</mi><mo>+</mo><msubsup><mi>X</mi><mi>t</mi><mi>T</mi></msubsup><msubsup><mi>&Phi;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>X</mi><mi>t</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0001182439280000019.GIF" wi="389" he="139" /></maths><maths num="0005"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><msup><mfenced open = "[" close = "]"><mtable><mtr><mtd><msubsup><mi>&sigma;</mi><mi>b</mi><mn>2</mn></msubsup></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA00011824392800000110.GIF" wi="443" he="93" /></maths>式中,符号<img file="FDA0001182439280000021.GIF" wi="75" he="63" />表示求取矩阵或向量的转置,符号<img file="FDA0001182439280000022.GIF" wi="75" he="63" />表示求取矩阵或向量的逆,<img file="FDA0001182439280000023.GIF" wi="142" he="62" /><img file="FDA0001182439280000024.GIF" wi="355" he="63" />分别表示“t”、“t‑1”、“t+1”时刻该符号变量对应的值;λ为接近1但又小于1的正数,取λ=0.99;<img file="FDA0001182439280000025.GIF" wi="54" he="63" />为传感器输出信号的噪声方差,由传感器的特性参数可以获得;I为d维单位方阵;<img file="FDA0001182439280000026.GIF" wi="74" he="63" />初值设置为单位方阵;W<sub>t</sub>=[w<sub>t</sub> w<sub>t‑1</sub> … w<sub>t‑d‑1</sub>]<sup>T</sup>为自适应滤波器的d维权值向量,其初值选为零向量,d为滤波器阶数;<img file="FDA0001182439280000027.GIF" wi="552" he="78" />为自适应滤波器d维输入向量,当采样时刻t&lt;d时其值选为0,且有:<img file="FDA0001182439280000028.GIF" wi="677" he="147" />(3)采集t+1时刻的电流互感器输出值x<sub>t+1</sub>,将电流互感器输出值x<sub>t+1</sub>与步骤(2)中自适应滤波器的输出值<img file="FDA0001182439280000029.GIF" wi="65" he="62" />相减,并根据步骤(1)中的公式进行判别;由此反复迭代计算、判断,直至检测到故障发生。
地址 100190 北京市海淀区中关村北二条6号
您可能感兴趣的专利