发明名称 一种新型电流互感器饱和识别方法
摘要 一种新型电流互感器饱和识别方法,包括如下步骤:1、采集电流互感器保护安装处的三相电压、电流瞬时值数据,将采样结果记为u<sub>a</sub>(k),u<sub>b</sub>(k),u<sub>c</sub>(k),i<sub>a</sub>(k),i<sub>b</sub>(k),i<sub>c</sub>(k);2、计算三相电压、电流故障分量Δu<sub>a</sub>(k),Δu<sub>b</sub>(k),Δu<sub>c</sub>(k),Δi<sub>a</sub>(k),Δi<sub>b</sub>(k),Δi<sub>c</sub>(k);3、计算等效电感和电阻参数L<sub>eq.j</sub>,R<sub>eq.j</sub>(j=a,b,c);4、计算等效电感参数的离散度<img file="DDA00003242389400011.GIF" wi="90" he="63" />(j=a,b,c);5、判据识别各相电流互感器是否饱和:若<img file="DDA00003242389400012.GIF" wi="188" he="55" />则判断a相电流互感器发生饱和;若<img file="DDA00003242389400013.GIF" wi="186" he="55" />则判断b相电流互感器发生饱和;若<img file="DDA00003242389400014.GIF" wi="180" he="55" />则判断c相电流互感器发生饱和;本发明方法具有简单可靠、仅使用短数据窗、不受系统运行方式和故障初相角影响的的特点。
申请公布号 CN103323698B 申请公布日期 2016.03.02
申请号 CN201310196547.9 申请日期 2013.05.24
申请人 西安交通大学 发明人 康小宁;姚旭;焦在滨;宋国兵;马超
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 何会侠
主权项 一种新型电流互感器饱和识别方法,其特征在于,包括如下步骤:步骤一、采集电流互感器保护安装处的三相电压、电流瞬时值数据,将采样结果记为u<sub>a</sub>(k),u<sub>b</sub>(k),u<sub>c</sub>(k),i<sub>a</sub>(k),i<sub>b</sub>(k),i<sub>c</sub>(k);步骤二、计算三相电压、电流故障分量Δu<sub>a</sub>(k),Δu<sub>b</sub>(k),Δu<sub>c</sub>(k),Δi<sub>a</sub>(k),Δi<sub>b</sub>(k),Δi<sub>c</sub>(k),计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close='' separators=''><mtable><mtr><mtd><msub><mi>&Delta;u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>pN</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>pN</mi><mo>)</mo></mrow></mtd></mtr></mtable><mo>,</mo></mfenced><mi>j</mi><mo>=</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>c</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003242389100011.GIF" wi="1406" he="156" /></maths>式中:p取正整数,数值大小由保护装置决定;N为每工频周期采样点数;步骤三、计算等效电感和电阻参数L<sub>eq.j</sub>,R<sub>eq.j</sub>(j=a,b,c),计算方法如下:取出一定数据窗长的某相电压电流故障分量,构成以下超定方程组:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mfrac><mrow><mi>&Delta;</mi><msub><mi>i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msub><mi>T</mi><mi>s</mi></msub></mrow></mfrac></mtd><mtd><mi>&Delta;</mi><msub><mi>i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msub><mi>T</mi><mi>s</mi></msub></mrow></mfrac></mtd><mtd><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msub><mi>T</mi><mi>s</mi></msub></mrow></mfrac></mtd><mtd><msub><mi>&Delta;i</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>L</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&Delta;u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&Delta;u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msub><mi>&Delta;u</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003242389100012.GIF" wi="1558" he="465" /></maths>式中:T<sub>s</sub>为采样时间间隔;将式(2)简写为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>[</mo><mi>i</mi><mo>]</mo><mo>&CenterDot;</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>L</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mo>[</mo><mi>u</mi><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003242389100013.GIF" wi="1260" he="160" /></maths>采用最小二乘算法计算参数L<sub>eq.j</sub>,R<sub>eq.j</sub>,计算公式如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>L</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mrow><mo>(</mo><msup><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msup><mo>&CenterDot;</mo><mo>[</mo><mi>i</mi><mo>]</mo><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&CenterDot;</mo><msup><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow><mi>T</mi></msup><mo>&CenterDot;</mo><mo>[</mo><mi>u</mi><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003242389100014.GIF" wi="1285" he="150" /></maths>步骤四、计算等效电感参数的离散度<img file="FDA00003242389100017.GIF" wi="373" he="78" />计算公式如下:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><msub><mi>D</mi><mi>L</mi></msub><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>|</mo><msub><mi>L</mi><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mover><mi>L</mi><mo>&OverBar;</mo></mover><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mover><mi>L</mi><mo>&OverBar;</mo></mover><mrow><mi>eq</mi><mo>.</mo><mi>j</mi></mrow></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00003242389100015.GIF" wi="1270" he="205" /></maths>式中:L<sub>eq.j</sub>(k)(k=1,2,…,M)为计算出的一定数据窗内的M个电感参数;<img file="FDA00003242389100016.GIF" wi="91" he="84" />为相同数据窗内M个电感参数的平均值;步骤五、按照以下判据识别各相电流互感器是否饱和:若<img file="FDA00003242389100021.GIF" wi="204" he="78" />则判断a相电流互感器发生饱和;若<img file="FDA00003242389100022.GIF" wi="196" he="75" />则判断b相电流互感器发生饱和;若<img file="FDA00003242389100023.GIF" wi="203" he="78" />则判断c相电流互感器发生饱和;其中:ε为整定门槛,取为0.2。
地址 710049 陕西省西安市咸宁路28号