发明名称 汽轮发电机组低频振动实时预警方法
摘要 本发明公开了旋转机械振动状态监测与故障诊断领域的一种汽轮发电机组低频振动实时预警方法。通过采集计算存储汽轮发电机组转子轴相对振动中的低频振动数据。在此基础上,结合FFT频谱分析方法及ANFIS自适应神经模糊推理系统计算方法,通过计算分析低频振动统计性内禀参数变化情况,对汽轮发电机组转子的低频振动异常状态进行预警。本发明提供的汽轮发电机组低频振动实时预警方法,利用机组运行中转子振动数据及ANFIS方法,对机组运行中转子振动数据进行实时自动在线监测、分析及判别,提高大型汽轮发电机组低频振动实时监测预测分析工作的效率和准确度。
申请公布号 CN103335708A 申请公布日期 2013.10.02
申请号 CN201310233932.6 申请日期 2013.06.13
申请人 华北电力大学 发明人 宋光雄
分类号 G01H17/00(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 陈波
主权项 1.一种汽轮发电机组低频振动实时预警方法,其特征在于,所述方法具体包括以下步骤:步骤1:采用高速振动数据采集卡实时采集汽轮发电机组转子一侧支持轴承的轴相对振动数据、转子的转速信号和键相信号;步骤2:利用FFT频谱分析方法对当前时刻支持轴承的轴相对振动数据进行频谱分析,得到从低频到高频的不同振动频率所对应的振动幅值数据序列;步骤3:从步骤2得到的振动幅值数据序列中,截取当前时刻所有小于机组工作转速对应频率的低频频率,并计算所有低频频率对应的振动幅值之和A;存储低频振动幅值A,每隔Δ<sub>s</sub>秒存储一次;步骤4:判断低频振动幅值存储时间是否大于预设时间段长度P<sub>MN</sub>,如果大于预设时间段长度,那么记录当前时刻T<sub>N</sub>前的低频振动幅值数据A,进入步骤5;否则,继续存储数据;步骤5:按照数据存储时间先后顺序排序,不同时刻下低频振动幅值数据以下标i表示数据存储时间先后顺序,i=1,2,3,…,m;从T<sub>N</sub>时刻向前截取至T<sub>M</sub>时刻的低频振动幅值数据A,将T<sub>M</sub>时刻至T<sub>N</sub>时刻的低频振动幅值数据表示为A<sub>i</sub>(i=1,2,3,…,m);|T<sub>N</sub>-T<sub>M</sub>|=p<sub>MN</sub>,p<sub>MN</sub>为预设时间段长度,<maths num="0001"><![CDATA[<math><mrow><mi>m</mi><mo>=</mo><mfrac><msub><mi>p</mi><mi>MN</mi></msub><msub><mi>&Delta;</mi><mi>s</mi></msub></mfrac><mo>;</mo></mrow></math>]]></maths>步骤6:根据低频振动幅值数据A<sub>i</sub>(i=1,2,3,…,m),建立ANFIS预测学习训练数列<maths num="0002"><![CDATA[<math><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>A</mi><msub><mi>s</mi><mn>1</mn></msub></msub></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mn>2</mn></msub></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub></msub></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mi>j</mi></msub></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>p</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>p</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>并进行训练,得到模糊推理系统FIS的规则组R<sub>s</sub>;步骤7:将步骤6中的学习训练数列中的前j列构成新的数列,<maths num="0003"><![CDATA[<math><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>A</mi><msub><mi>s</mi><mn>1</mn></msub></msub></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mn>2</mn></msub></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub></msub></mtd><mtd><msub><mi>A</mi><msub><mi>s</mi><mi>j</mi></msub></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>1</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mn>2</mn></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mrow><mi>j</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced></math>]]></maths>利用规则组R<sub>s</sub>,进行基于ANFIS方法的模糊推理计算,得到训练学习输出项序列<img file="FDA00003340086500023.GIF" wi="527" he="88" />步骤8:计算<img file="FDA00003340086500024.GIF" wi="515" he="79" />与训练预测项序列<maths num="0004"><![CDATA[<math><mrow><msub><mi>A</mi><mrow><msub><mi>s</mi><mi>j</mi></msub><mo>+</mo><mi>p</mi><mo>+</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>r</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mrow></math>]]></maths>的偏差<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>&Delta;</mi><mi>k</mi><mi>T</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>r</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>步骤9:计算<img file="FDA00003340086500027.GIF" wi="493" he="82" />的标准偏差σ<sup>T</sup>及最大值<img file="FDA00003340086500028.GIF" wi="132" he="80" />步骤10:根据任意当前时刻低频振幅数据<img file="FDA00003340086500029.GIF" wi="80" he="77" />及其前的低频振动幅值数据,得到当前时刻的标准偏差σ<sup>N</sup>及最大值<img file="FDA000033400865000210.GIF" wi="125" he="86" />步骤11:如果满足σ<sup>N</sup>&gt;w<sub>1</sub>×σ<sup>T</sup>或<img file="FDA000033400865000211.GIF" wi="357" he="87" />其中w<sub>1</sub>、w<sub>2</sub>为预警阈值系数,判定所测得的机组轴系转子一侧的当前时刻低频振动出现异常状态。
地址 102206 北京市昌平区朱辛庄北农路2号