发明名称 采用HS变换和单类支持向量机的高压断路器机械状态监测方法
摘要 本发明是一种采用HS变换和单类支持向量机的高压断路器机械状态监测方法,其特点是,包括高压断路器振动信号采集、将振动信号利用HS变换进行处理、振动信号特征提取和单类支持向量机判断断路器机械状态是否正常等步骤。与以往的高压断路器机械状态监测方法相比,本方法中使用的HS变换方法具有较高的抗噪性和鲁棒性,而且相对容易在实际工程环境下实现。采用单类支持向量机方法判断高压断路器机械状态是否正常,不依赖各种故障样本信号,只用一类正常振动信号样本即可完成分类器的训练,具有科学合理,实用性强,推广应用价值较高等优点。
申请公布号 CN105241643A 申请公布日期 2016.01.13
申请号 CN201510599527.5 申请日期 2015.09.19
申请人 东北电力大学 发明人 黄南天;张书鑫;蔡国伟;陈怀金;方立华
分类号 G01M13/00(2006.01)I;G01H11/08(2006.01)I 主分类号 G01M13/00(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种采用HS变换和单类支持向量机的高压断路器机械状态监测方法,其特征是,它包括以下步骤:1)高压断路器振动信号采集通过压电式加速度传感器和数据采集系统采集高压断路器振动信号;2)将振动信号利用Hyperbolic S变换进行处理,亦称HS变换进行处理Hyperbolic S变换是S变换的扩展,其关键在于采用Hyperbolic窗函数代替传统S变换使用的高斯窗函数,具体计算方法如下:离散HS变换公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mo>&lsqb;</mo><mi>n</mi><mo>,</mo><mi>j</mi><mo>&rsqb;</mo><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>H</mi><mo>&lsqb;</mo><mi>m</mi><mo>+</mo><mi>n</mi><mo>&rsqb;</mo><mi>G</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>i</mi><mn>2</mn><mi>&pi;</mi><mi>m</mi><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806612780000011.GIF" wi="1404" he="127" /></maths>其中N为总采样点数;n=0,1,...,N‑1;m=0,1,...,N‑1;j=0,1,...,N‑1;G(m,n)为变换过程中所使用的Hyperbolic窗函数,可记做<img file="FDA0000806612780000012.GIF" wi="748" he="150" />其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>&alpha;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>&beta;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn><msub><mi>&alpha;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><msub><mi>&beta;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub></mrow></mfrac><mi>t</mi><mo>+</mo><mfrac><mrow><mo>(</mo><msub><mi>&alpha;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><mo>-</mo><msub><mi>&beta;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn><msub><mi>&alpha;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub><msub><mi>&beta;</mi><mrow><mi>h</mi><mi>y</mi></mrow></msub></mrow></mfrac><msqrt><mrow><msup><mi>t</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>&lambda;</mi><mrow><mi>h</mi><mi>y</mi></mrow><mn>2</mn></msubsup></mrow></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000806612780000013.GIF" wi="766" he="135" /></maths>3)振动信号特征提取通过统计实验分析发现,频率高于10kHz的振动信号成分振幅很小,因此,主要分析0‑10kHz频段,将0‑10kHz、0‑150ms的时‑频平面划分为300个频宽为1kHz、时宽为5ms的等大小时频块,通过分析振动信号在该部分的分布情况,分析断路器工作状态;通过时‑频能量熵描述断路器机械振动信号能量时‑频分布特征;设整个时‑频平面的总能量为E,时‑频块S<sub>ij</sub>的能量为E<sub>ij</sub>(i=1,…,10;j=1,…,30),首先对E<sub>ij</sub>做归一化处理如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>E</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mi>E</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806612780000014.GIF" wi="1125" he="126" /></maths>那么对于10个频带中的每个频带,计算该频带中时‑频块能量相对于时间的信息熵公式为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>H</mi><mi>i</mi></msub><mo>=</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>30</mn></munderover><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>lg</mi><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806612780000015.GIF" wi="1300" he="142" /></maths>其中H<sub>i</sub>反映第i个频带上的能量随时间分布的均匀程度,此外,由实验发现振动信号能量集中在收到触发信号后的前50ms内;因此,为了提高延时类故障的识别能力,计算前10个时宽对应的时‑频块相对于频率的信息熵,其公式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>L</mi><mi>j</mi></msub><mo>=</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>10</mn></munderover><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&CenterDot;</mo><mi>lg</mi><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>...</mo><mo>,</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000806612780000016.GIF" wi="1302" he="142" /></maths>由以上20个时‑频能量熵构成信号的时‑频能量熵向量X=[H L],其中H=[H<sub>1</sub>,…,H<sub>10</sub>],L=[L<sub>1</sub>,…,L<sub>10</sub>],向量X作为高压断路器机械状态监测的特征向量;4)单类支持向量机判断断路器运行状态每个训练样本通过3)中所介绍的特征提取方法可得到20维输入特征向量样本集,将该训练样本集作为单类支持向量机的输入变量,计算得到特定的决策超平面,完成单类支持向量机的训练过程;当需要判断断路器运行状态时,将采集到的需诊断的样本按3)中介绍的特征提取方法得到输入向量,计算该样本与决策超平面的关系,判断高压断路器机械状态是否正常。
地址 132012 吉林省吉林市船营区长春路169号