发明名称 基于CWD谱峭度的暂态电能质量扰动分类识别方法
摘要 本发明提出了一种基于CWD谱峭度的暂态电能质量扰动分类识别方法,并将其与有效值相结合应用于五种单一和复合的暂态电能质量扰动的识别。该识别方法首先将电压暂升、暂降和中断这三种幅值类扰动看作一类,通过CWD谱峭度值将电力局数据采集系统采集到的数据分为暂态脉冲、暂态振荡和幅值类扰动三类,再通过计算幅值类扰动的有效值将数据分为电压暂升、电压暂降和电压中断,并向后续处理设备输出所求分类结果。本发明直接利用数值的大小区分不同的扰动类型,不需要使用任何分类器,极大的简化了识别的流程和识别的时间。本发明方法能够不仅能准确区分暂态脉冲、暂态振荡及暂升、暂降和中断这五种单一暂态扰动及其组合而成的复合扰动,还具有良好的抗噪性。
申请公布号 CN103018537A 申请公布日期 2013.04.03
申请号 CN201210495025.4 申请日期 2012.11.28
申请人 西南交通大学 发明人 刘志刚;朱玲;胡巧琳;张巧革
分类号 G01R19/02(2006.01)I;G01R23/16(2006.01)I 主分类号 G01R19/02(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 1.一种基于CWD谱峭度的暂态电能质量扰动分类识别方法,在电力系统中对暂态扰动进行识别,其主要步骤是:A、提取扰动特征信号输入含有扰动信号的电压信号u(n),n=1,2,...,N,N为数据长度,将u(n)进行小波变换,提取扰动特征信号u<sub>r</sub>(n);B、计算扰动特征信号谱峭度计算由A获得的扰动特征信号u<sub>r</sub>(n)的Choi-Williams分布Cu<sub>r</sub>(t,f),再根据Cu<sub>r</sub>(t,f)的2n阶瞬时谱矩得出2阶和4阶瞬时谱距,最后根据谱峭度的定义,求得u<sub>r</sub>(n)的谱峭度<img file="FDA00002484106000011.GIF" wi="164" he="68" />具体作法为:B1、扰动特征信号u<sub>r</sub>(n),计算其Choi-Williams分布,结果为Cu<sub>r</sub>(t,f),<maths num="0001"><![CDATA[<math><mrow><msub><mi>Cu</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></munderover><msqrt><mfrac><mi>&sigma;</mi><mrow><mn>4</mn><msup><mi>&pi;&tau;</mi><mn>2</mn></msup></mrow></mfrac></msqrt><msup><mi>e</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>&sigma;t</mi><mn>2</mn></msup><mrow><mn>4</mn><msup><mi>&tau;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></msup><msub><mi>u</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>f</mi><mo>+</mo><mfrac><mi>&tau;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><msup><msub><mi>u</mi><mi>r</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mi>f</mi><mo>-</mo><mfrac><mi>&tau;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;f&tau;</mi></mrow></msup><mi>d&tau;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>B2、根据Choi-Williams分布Cu<sub>r</sub>(t,f)求取的2阶瞬时谱距<img file="FDA00002484106000013.GIF" wi="116" he="63" />和4阶瞬时谱距<img file="FDA00002484106000014.GIF" wi="139" he="63" /><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>S</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>&lt;</mo><msup><mrow><mo>|</mo><msub><mi>Cu</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>></mo></mrow><mi>k</mi></msub></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>S</mi><mo>^</mo></mover><mn>4</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>&lt;</mo><msup><mrow><mo>|</mo><msub><mi>Cu</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>4</mn></msup><mo>></mo></mrow><mi>k</mi></msub></mrow></math>]]></maths>其中,&lt;·&gt;<sub>k</sub>表示的是k阶时平均;B3、根据谱峭度定义,求得基于CWD的谱峭度,<maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>K</mi><mo>^</mo></mover><mi>x</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mover><mi>S</mi><mo>^</mo></mover><mrow><mn>4</mn><mi>x</mi></mrow></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>S</mi><mo>^</mo></mover><mrow><mn>2</mn><mi>x</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mn>2</mn><mo>,</mo></mrow></math>]]></maths>f≠0    (3);C、根据阈值分类识别由B算出的谱峭度的最大值K<sub>max</sub>和谱峭度的最小值K<sub>min</sub>的大小,设置合适的阈值,判定输入信号属于幅值类扰动、暂态脉冲或者是暂态振荡;D、有效值方法区分幅值类扰动将未区分出来的幅值类扰动,计算其有效值曲线,通过幅值的大小范围来判断是电压暂升、电压暂降或者是电压中断,得出并向后续处理设备输出所求分类结果。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处