发明名称 一种基于ESPRIT与PSA的笼型异步电动机转子断条故障检测方法
摘要 一种基于ESPRIT与PSA的笼型异步电动机转子断条故障检测方法,它首先对按一定频率采集的定子电流瞬时信号应用旋转不变技术,得到其基波分量与边频分量的准确的频率值以及粗糙的幅值、粗糙的初相角;再应用模式搜索算法估计定子电流瞬时信号基波分量与边频分量的准确的幅值、初相角;进而求得当前边频分量与基波分量幅值的比值并把它作为故障特征;然后求出该比值与检测阈值之比确定故障指数;最后依故障指数判断是否存在转子断条故障。本发明利用很少的定子电流信号采样点数即可高灵敏度、高可靠性地在线检测异步电动机转子断条故障,有效克服了负荷波动、噪声等不利因素的影响,非常适用于异步电动机低转差率运行情况。
申请公布号 CN102279341A 申请公布日期 2011.12.14
申请号 CN201110207098.4 申请日期 2011.07.23
申请人 华北电力大学(保定) 发明人 许伯强;孙丽玲
分类号 G01R31/02(2006.01)I;G01R31/34(2006.01)I 主分类号 G01R31/02(2006.01)I
代理机构 石家庄冀科专利商标事务所有限公司 13108 代理人 李羡民;高锡明
主权项 1.一种基于ESPRIT与PSA的笼型异步电动机转子断条故障检测方法,其特征是,它首先对按一定频率采集的定子电流瞬时信号应用旋转不变技术,得到其基波分量与边频分量的准确的频率值以及粗糙的幅值、粗糙的初相角;再应用模式搜索算法估计定子电流瞬时信号基波分量与边频分量的准确的幅值、初相角;进而求得当前边频分量与基波分量幅值的比值并把它作为故障特征;然后求出该比值与检测阈值之比确定故障指数;最后依故障指数判断是否存在转子断条故障;具体步骤如下:a. 测取一相定子电流瞬时信号<img file="170007DEST_PATH_IMAGE001.GIF" wi="14" he="25" />:对于高压电机,采用一只电流钳在电流互感器CT二次侧测取一相定子电流瞬时信号;对于低压电机,采用一只电流钳直接在电机接线端子处测取一相定子电流瞬时信号;b. 采用单工频周期滑动窗方法计算定子电流瞬时信号<img file="113693DEST_PATH_IMAGE002.GIF" wi="14" he="25" />的有效值<img file="141691DEST_PATH_IMAGE003.GIF" wi="17" he="25" />,通过分析有效值的变化趋势,提取其最平稳亦即波动最小的一段数据,记为<img file="605165DEST_PATH_IMAGE004.GIF" wi="16" he="26" />;c. 对提取的定子电流瞬时信号<img file="325996DEST_PATH_IMAGE004.GIF" wi="16" he="26" />应用旋转不变技术,确定其基波分量与边频分量的准确的频率值以及粗糙的幅值、初相角:提取的定子电流瞬时信号<img file="225819DEST_PATH_IMAGE004.GIF" wi="16" he="26" />可以表示为一系列余弦谐波分量之组合:<img file="323088DEST_PATH_IMAGE005.GIF" wi="282" he="46" />,其中,<img file="155915DEST_PATH_IMAGE006.GIF" wi="18" he="25" />表示采样周期;<img file="516489DEST_PATH_IMAGE007.GIF" wi="20" he="20" />表示采样点数;<img file="684034DEST_PATH_IMAGE008.GIF" wi="17" he="18" />表示谐波个数;<img file="584994DEST_PATH_IMAGE009.GIF" wi="18" he="25" />、<img file="475590DEST_PATH_IMAGE010.GIF" wi="17" he="25" />、<img file="803803DEST_PATH_IMAGE011.GIF" wi="16" he="25" />分别表示第<img file="943797DEST_PATH_IMAGE012.GIF" wi="9" he="18" />个谐波的幅值、频率、初相角,则基波分量与边频分量的准确的频率值以及粗糙的幅值、初相角按下列步骤求得:(a) 、定义<img file="399180DEST_PATH_IMAGE013.GIF" wi="99" he="22" />,引入以下<img file="409862DEST_PATH_IMAGE014.GIF" wi="44" he="20" />阶矩阵(保证<img file="908976DEST_PATH_IMAGE015.GIF" wi="49" he="20" />):<img file="536267DEST_PATH_IMAGE016.GIF" wi="276" he="31" />;<img file="247871DEST_PATH_IMAGE017.GIF" wi="280" he="31" />;构造<img file="644217DEST_PATH_IMAGE018.GIF" wi="36" he="22" />的自相关矩阵为:<img file="563500DEST_PATH_IMAGE019.GIF" wi="145" he="27" />;而<img file="146928DEST_PATH_IMAGE018.GIF" wi="36" he="22" />和<img file="927803DEST_PATH_IMAGE020.GIF" wi="36" he="22" />的互相关矩阵为:<img file="444235DEST_PATH_IMAGE021.GIF" wi="145" he="27" />;其中,<img file="753993DEST_PATH_IMAGE022.GIF" wi="16" he="18" />表示数学期望,<img file="355876DEST_PATH_IMAGE023.GIF" wi="18" he="18" />表示共轭转置;(b) 、对<img file="691173DEST_PATH_IMAGE024.GIF" wi="33" he="25" />进行特征值分解,确定其最小特征值<img file="530953DEST_PATH_IMAGE025.GIF" wi="22" he="22" />;(c) 、计算<img file="542772DEST_PATH_IMAGE026.GIF" wi="105" he="26" />,<img file="366371DEST_PATH_IMAGE027.GIF" wi="13" he="18" />表示<img file="223469DEST_PATH_IMAGE028.GIF" wi="18" he="16" />阶单位阵;(d) 、计算<img file="448914DEST_PATH_IMAGE029.GIF" wi="110" he="26" />,<img file="880901DEST_PATH_IMAGE030.GIF" wi="17" he="18" />为一<img file="926217DEST_PATH_IMAGE028.GIF" wi="18" he="16" />阶阵,<img file="852585DEST_PATH_IMAGE031.GIF" wi="78" he="49" />(此处,<img file="932537DEST_PATH_IMAGE027.GIF" wi="13" he="18" />表示<img file="20578DEST_PATH_IMAGE032.GIF" wi="37" he="20" />阶单位阵);(e) 、对<img file="818770DEST_PATH_IMAGE033.GIF" wi="22" he="25" />进行奇异值分解<img file="299561DEST_PATH_IMAGE034.GIF" wi="81" he="26" />,此处<img file="702861DEST_PATH_IMAGE035.GIF" wi="94" he="28" />,<img file="492962DEST_PATH_IMAGE036.GIF" wi="96" he="52" />(<img file="778450DEST_PATH_IMAGE037.GIF" wi="21" he="25" />由<img file="515462DEST_PATH_IMAGE008.GIF" wi="17" he="18" />个主奇异值组成),<img file="570005DEST_PATH_IMAGE038.GIF" wi="80" he="52" />;(f) 、计算矩阵<img file="791994DEST_PATH_IMAGE039.GIF" wi="58" he="26" />;(g) 、对<img file="768041DEST_PATH_IMAGE040.GIF" wi="95" he="27" />进行广义特征值分解,确定<img file="839902DEST_PATH_IMAGE008.GIF" wi="17" he="18" />个广义特征值<img file="14531DEST_PATH_IMAGE041.GIF" wi="104" he="25" />(其余<img file="349698DEST_PATH_IMAGE042.GIF" wi="42" he="18" />个广义特征值恒等于0);(h) 、根据广义特征值,确定采样信号各个分量的频率<img file="609778DEST_PATH_IMAGE043.GIF" wi="146" he="54" />,<img file="970483DEST_PATH_IMAGE044.GIF" wi="50" he="28" />、<img file="468461DEST_PATH_IMAGE045.GIF" wi="50" he="28" />分别表示特征值<img file="505687DEST_PATH_IMAGE046.GIF" wi="17" he="25" />的虚部、实部;(i) 、计算矩阵<img file="253063DEST_PATH_IMAGE047.GIF" wi="189" he="100" />;(j) 、计算矩阵<img file="135568DEST_PATH_IMAGE048.GIF" wi="112" he="33" />,此处<img file="19211DEST_PATH_IMAGE049.GIF" wi="13" he="16" />为一列向量<img file="742185DEST_PATH_IMAGE050.GIF" wi="148" he="33" />,而<img file="180119DEST_PATH_IMAGE051.GIF" wi="18" he="18" />为列向量<img file="131895DEST_PATH_IMAGE052.GIF" wi="166" he="30" />;(k) 确定采样信号各个分量的幅值与初相角<img file="135623DEST_PATH_IMAGE053.GIF" wi="62" he="28" />、<img file="983493DEST_PATH_IMAGE054.GIF" wi="126" he="54" />,<img file="439882DEST_PATH_IMAGE055.GIF" wi="77" he="22" />;d. 应用模式搜索算法估计提取的定子电流瞬时信号<img file="946081DEST_PATH_IMAGE004.GIF" wi="16" he="26" />基波分量与边频分量的准确的幅值、初相角:首先构造目标函数:提取的定子电流瞬时信号<img file="273157DEST_PATH_IMAGE004.GIF" wi="16" he="26" />可以表示为:<img file="88667DEST_PATH_IMAGE056.GIF" wi="464" he="46" />,生成<img file="32352DEST_PATH_IMAGE057.GIF" wi="42" he="22" />矩阵<img file="794771DEST_PATH_IMAGE058.GIF" wi="40" he="25" />、<img file="507513DEST_PATH_IMAGE059.GIF" wi="41" he="25" />,具体如下:<img file="743191DEST_PATH_IMAGE060.GIF" wi="494" he="36" />,<img file="643014DEST_PATH_IMAGE061.GIF" wi="489" he="36" />,令状态<img file="740283DEST_PATH_IMAGE062.GIF" wi="86" he="28" />,其中<img file="307530DEST_PATH_IMAGE063.GIF" wi="18" he="25" />、<img file="933684DEST_PATH_IMAGE064.GIF" wi="20" he="25" />分别为:<img file="851961DEST_PATH_IMAGE065.GIF" wi="273" he="30" />;<img file="238074DEST_PATH_IMAGE066.GIF" wi="268" he="30" />;构造如下目标函数:<img file="394249DEST_PATH_IMAGE067.GIF" wi="208" he="30" />;此处,<img file="722462DEST_PATH_IMAGE051.GIF" wi="18" he="18" />为列向量<img file="862456DEST_PATH_IMAGE052.GIF" wi="166" he="30" />,<img file="35949DEST_PATH_IMAGE062.GIF" wi="86" he="28" />为使<img file="577789DEST_PATH_IMAGE068.GIF" wi="37" he="22" />取最小值的待定状态;之后,按如下步骤估计提取的定子电流瞬时信号<img file="326171DEST_PATH_IMAGE004.GIF" wi="16" he="26" />基波分量与边频分量的准确的幅值、初相角:(a’)、 根据旋转不变技术的计算结果设定初始状态<img file="422303DEST_PATH_IMAGE069.GIF" wi="20" he="22" />,并给定轴向方向<img file="399486DEST_PATH_IMAGE070.GIF" wi="16" he="25" />,<img file="61411DEST_PATH_IMAGE070.GIF" wi="16" he="25" />,…<img file="200269DEST_PATH_IMAGE071.GIF" wi="17" he="25" />,步长<img file="314855DEST_PATH_IMAGE072.GIF" wi="16" he="20" />,减缩率<img file="846462DEST_PATH_IMAGE073.GIF" wi="65" he="28" />,终止参数<img file="831735DEST_PATH_IMAGE074.GIF" wi="14" he="16" />,令<img file="672652DEST_PATH_IMAGE075.GIF" wi="52" he="25" />;(b’)、探索移动, 对<img file="274535DEST_PATH_IMAGE076.GIF" wi="22" he="26" />(<img file="327942DEST_PATH_IMAGE077.GIF" wi="80" he="22" />)依次进行如下轴向搜索:令<img file="698880DEST_PATH_IMAGE078.GIF" wi="92" he="26" />,若<img file="959966DEST_PATH_IMAGE079.GIF" wi="92" he="25" />,则令<img file="252407DEST_PATH_IMAGE080.GIF" wi="97" he="26" />;否则,令<img file="640663DEST_PATH_IMAGE081.GIF" wi="90" he="26" />,若<img file="866108DEST_PATH_IMAGE079.GIF" wi="92" he="25" />,则令<img file="517669DEST_PATH_IMAGE082.GIF" wi="97" he="26" />;(c’)、模式移动,若<img file="94144DEST_PATH_IMAGE083.GIF" wi="98" he="25" />,则令<img file="771244DEST_PATH_IMAGE084.GIF" wi="122" he="25" />,以<img file="320037DEST_PATH_IMAGE085.GIF" wi="18" he="22" />为新的初始状态,转(b’),得到新的迭代点<img file="939237DEST_PATH_IMAGE086.GIF" wi="18" he="25" />——若<img file="737429DEST_PATH_IMAGE087.GIF" wi="94" he="25" />,则令<img file="670750DEST_PATH_IMAGE088.GIF" wi="49" he="25" />;否则,令<img file="870787DEST_PATH_IMAGE089.GIF" wi="64" he="22" />;(d’)、 若<img file="910156DEST_PATH_IMAGE090.GIF" wi="40" he="20" />,则停止;否则,转(b’);最终获得的<img file="398907DEST_PATH_IMAGE062.GIF" wi="86" he="28" />可使<img file="932656DEST_PATH_IMAGE068.GIF" wi="37" he="22" />取最小值,实际上就是确定了采样信号各个分量准确的幅值与初相角<img file="987200DEST_PATH_IMAGE009.GIF" wi="18" he="25" />、<img file="151465DEST_PATH_IMAGE011.GIF" wi="16" he="25" />,<img file="924249DEST_PATH_IMAGE055.GIF" wi="77" he="22" />;e. 确定提取的定子电流瞬时信号<img file="746842DEST_PATH_IMAGE004.GIF" wi="16" he="26" />的<img file="124734DEST_PATH_IMAGE091.GIF" wi="65" he="24" />边频分量与<img file="256638DEST_PATH_IMAGE092.GIF" wi="18" he="24" />基频分量幅值之比<img file="251139DEST_PATH_IMAGE093.GIF" wi="72" he="26" />+<img file="595533DEST_PATH_IMAGE094.GIF" wi="72" he="26" />,其中,<img file="624669DEST_PATH_IMAGE093.GIF" wi="72" he="26" />是<img file="911162DEST_PATH_IMAGE095.GIF" wi="65" he="24" />边频分量与<img file="127380DEST_PATH_IMAGE092.GIF" wi="18" he="24" />基频分量幅值之比,<img file="541044DEST_PATH_IMAGE094.GIF" wi="72" he="26" />是<img file="424686DEST_PATH_IMAGE096.GIF" wi="65" he="24" />边频分量与<img file="101655DEST_PATH_IMAGE092.GIF" wi="18" he="24" />基频分量幅值之比;f. 确定故障指数:根据常规经验设置检测阈值2% ,<img file="70748DEST_PATH_IMAGE093.GIF" wi="72" he="26" />+<img file="38835DEST_PATH_IMAGE094.GIF" wi="72" he="26" />与检测阈值的比值即为故障指数;g. 根据故障指数判断故障存在与否:故障指数数值&lt;1,表示电机处于健康状态,且其数值愈小,健康状态愈明确;故障指数数值&gt;1,表示电机处于故障状态,且其数值愈大,故障状态愈严重。
地址 071003 河北省保定市北市区永华北大街619号