发明名称 一种电力异常故障数据分析装置及诊断方法
摘要 本发明涉及一种电力异常故障数据分析装置及诊断方法,属于设备故障检测领域,分析装置包括信号采集单元、旋转电容滤波器电路、信号处理器、电源模块、通讯模块和上位机。采用该装置对异常故障数据诊断方法为:1、利用相空间重构的方法,求出嵌入延迟τ和最佳嵌入维d;2、对采集数据进行建模,得到神经网络预测原始模型;3、预测的下一步输出xi+1-;4、更新下一步预测的BP网络的权值;5、判断误差。本发明地优点为把BP神经网络和相空间重构的方法进行了巧妙的结合,适合多种采集信号类型、处理速度准确快速、适用于各种复杂储运工况,能够有效地对信号故障特别是细微的故障数据进行分析处理。旋转滤波电容器可以对噪音有很强的抑制能力。
申请公布号 CN101533058B 申请公布日期 2012.02.01
申请号 CN200910011275.4 申请日期 2009.04.24
申请人 东北大学 发明人 杨东升;张化光;李爱平;孙秋野;李营;刘博;邢颖;王迎春
分类号 G01R31/08(2006.01)I;G06N3/02(2006.01)I 主分类号 G01R31/08(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 朱光林
主权项 1.一种采用电力异常故障数据分析装置的异常故障数据诊断方法,其特征在于:所述电力异常故障数据分析装置包括信号采集单元、旋转电容滤波器电路、信号处理器、电源模块、通讯模块和上位机,信号采集单元包括第一级运算放大器、第二级反相比例运算放大器、第三级运算放大器、电阻、电容和二极管;旋转电容滤波器电路包括电阻、第四级运算放大器、第五级运算放大器、RC电路、第二双刀双掷电子开关、第三双刀双掷电子开关、第四双刀双掷电子开关和第五双刀双掷电子开关;电路连接是信号采集单元的输出端连接旋转电容滤波器电路的输入端;旋转电容滤波器电路的输出端连接信号处理器的输入端,信号处理器的输出端通过串行通讯连接上位机的输入端,电源模块连接信号处理器;其中信号采集单元由两个电阻组成的分压电路连接第一级运算放大器的同相输入端,第一级运算放大器的输出端通过电阻连接第二级反相比例运算放大器反相输入端,第二级反相比例运算放大器输出端通过电阻连接旋转电容滤波器电路的输入端,旋转电容滤波器电路的输出端连接第三级运算放大器的反相输入端,第三级运算放大器的输出端连接由电阻和电容组成的无源滤波器,无源滤波器的输出端连接二极管组成的限幅电路;所述的电力异常故障数据分析装置的异常故障数据诊断方法,按如下步骤进行:步骤一、采集数据,包括电压,利用相空间重构的方法,求出嵌入延迟τ和最佳嵌入维d,利用原始时间序列V重构相空间;步骤二、局部对电压的时间序列利用神经网络建模,对采集数据进行建模,得到神经网络预测原始模型;步骤三、采集数据x<sub>i</sub>联合原始数据组成输入向量,作为当前神经网络的输入,通过建模,得到预测的下一步输出<img file="FSB00000528825200011.GIF" wi="92" he="71" />步骤四、采集数据x<sub>i</sub>在内共k个向量x<sub>i</sub>,x<sub>i-1</sub>,…x<sub>i-k+1</sub>实时训练BP网络,得到的权值用来实时更新下一步预测的BP网络的权值;步骤五、x<sub>i</sub>和其预测值<img file="FSB00000528825200012.GIF" wi="40" he="65" />的预测误差值err作为故障判断的依据,如果err>RT,则表示当前数据有异常,否则继续判断,其中RT为预先设定或者通过自动识别得到的阈值;其中k表示需要确定在线训练的原始数据的个数,k=0,...,N,N为自然数,该参数k如果太大,计算机将无法在规定的时间内完成网络参数的训练,如果太小则无法体现要训练数据的特点,导致预测出现较大误差甚至错误;另外需要确定的一个参数是阈值RT,RT为预先设定或者通过自动识别得到的阈值,该参数RT作为判断数据是否发生异常的触发器,要根据神经网络的预测误差err来决定;如果err较小则适当减小RT,如果err变大,则要相应增大RT,事实上,k、err和RT是有关系的,在线训练的数据个数k越大,那么预测的误差err越小,那么RT也就能够适当的减小,并且这样能够发现更小的故障; 其中,步骤一中所述相空间重构的方法通过测量可得到采样信号的时间序列矩阵,采集数据{X<sub>k</sub>}={X(kT)},其中k=0,...,N,T为采样周期,每个X(kT)为kT时刻实测的m维向量,针对X(kT)的第i个分量构造形式上的状态向量<img file="FSB00000528825200021.GIF" wi="1347" he="82" />其中τ=[τ<sub>1</sub> τ<sub>2</sub>,...,τ<sub>m</sub>]<sup>T</sup>是未知的延迟时间向量,τ<sub>i</sub>是τ的第i个分量,i=1,...,m;n=[n<sub>1</sub> n<sub>2</sub>,...,n<sub>m</sub>]<sup>T</sup>,n<sub>i</sub>是n的第i个分量,是未知的嵌入维数,i=1,...,m;为了确定τ<sub>i</sub>和n<sub>i</sub>值,从而确定τ和n,首先利用时间序列计算<img file="FSB00000528825200022.GIF" wi="1314" he="78" />和<img file="FSB00000528825200023.GIF" wi="1470" he="142" />其中k′<k,k′=0,1,...,E(·)是期望值,<img file="FSB00000528825200024.GIF" wi="59" he="62" />是平均值,<img file="FSB00000528825200025.GIF" wi="156" he="62" />是线性相关的函数,<img file="FSB00000528825200026.GIF" wi="188" he="68" />是非线性相关的函数,设τ<sub>iX</sub>,<img file="FSB00000528825200027.GIF" wi="64" he="45" />分别是<img file="FSB00000528825200028.GIF" wi="153" he="59" />和<img file="FSB00000528825200029.GIF" wi="189" he="66" />达到第一个最小值的时间,则得到如下时间单元:<img file="FSB000005288252000210.GIF" wi="1055" he="85" />延迟时间常数τ<sub>i</sub>选择为<img file="FSB000005288252000211.GIF" wi="1188" he="94" />其中Int[·]表示实数的整数部分,通过采用相关函数来计算吸引子的相关维数:<img file="DEST_PATH_FSB000006049415000112.GIF" wi="820" he="158" />i=1,...,m    (6)其中v是一个常参数;θ为Heaviside函数,计算<img file="DEST_PATH_FSB000006049415000113.GIF" wi="459" he="99" />可得到吸引子的估计维数<img file="FSB000005288252000214.GIF" wi="57" he="50" />则第i个分量的嵌入维数为<img file="FSB000005288252000215.GIF" wi="228" he="51" />步骤四中所述BP网络的训练过程:1网络初始化,用一组随机数对网络赋初始权值,设置学习因子η、允许误差ε、网络结构,即网络层数L和每层节点数n<sub>1</sub>,令样本模式计数器P=1,训练次数计数器q=1;2为网络输入一组学习样本; 3对每个学习样本p循环;(1)逐层正向计算网络d个输入节点的输入X(t)和输出;(2)计算第p个样本的输出的误差E<sub>P</sub>和网络的总误差E,<img file="FSB00000528825200031.GIF" wi="1068" he="151" /><img file="FSB00000528825200032.GIF" wi="864" he="127" />式中:T<sub>jdp</sub>和<img file="FSB00000528825200033.GIF" wi="52" he="80" />分别为输出层第j个节点的期望输出和实际输出,E<sub>P</sub>为输出的误差,E为网络的总误差;(3)当E小于允许误差ε或者达到指定的迭代次数时,学习过程结束,否则,进行误差反向传播;(4)反向逐层计算网络各节点误差<img file="FSB00000528825200034.GIF" wi="86" he="64" />即对于输出层有:<img file="FSB00000528825200035.GIF" wi="1005" he="99" />对于隐含层有:<img file="DEST_PATH_FSB00000604941500026.GIF" wi="1160" he="114" />(5)修正网络连接权值:<img file="DEST_PATH_FSB00000604941500027.GIF" wi="1170" he="109" />式(9)、(10)、(11)中,z为学习次数,η为学习因子;<img file="DEST_PATH_FSB00000604941500028.GIF" wi="69" he="65" />表示第l层第j个节点p次学习的输出值,l取0,1,2,3时分别代表输入层、第一、第二隐含层和输出层;y<sub>jdp</sub>表示网络期望输出值;<img file="DEST_PATH_FSB00000604941500029.GIF" wi="107" he="86" />表示为第l+1层第s个节点与上一层第j个节点学习时的连接权重;<img file="FSB000005288252000310.GIF" wi="92" he="67" />表示第l-1层第j个节点p次学习的输出值;W<sub>ij</sub>(z)表示第z次修正网络连接权值,W<sub>ij</sub>(z+1)表示第z+1次修正网络连接权值;所述诊断方法的诊断过程最终由信号处理器DSPIC30F6010A处理器的控制实现,按以下控制步骤执行: 步骤一:开始;步骤二:定义程序出口地址并初始化I/O设备、同时初始化看门狗、中断向量、装置时钟;步骤四:设置异常响应,并设置中断响应地址和开中断;步骤五:初始化存储系统,为采样数据分配内存;步骤六:设置采样间隔及使能各采样通道;步骤七:DSPIC30F6010A微处理器对常规数据进行信号采集和利用BP神经网络进行分析处理;步骤八:对数据结果打包汇总,上传到上位机;步骤九:结束;其中步骤七中所述BP神经网络的分析过程包括以下步骤:步骤一:网络初始化,用一组随机数对网络赋初始权值;步骤二:为网络提供一组学习样本;步骤三:逐层正向计算网络各节点的输入和输出;步骤四:计算第p个样本的输出的误差E<sub>P</sub>和网络的总误差E;步骤五:修正网络连接权值;所述的旋转电容滤波器电路中第四运算放大器与第五级运算放大器组成并联电路,在第四级运算放大器的反相输入端连接第二双刀双掷电子开关,第二双刀双掷电子开关通过RC电路连接第三双刀双掷电子开关;在第五级运算放大器的反相输入端连接第四双刀双掷电子开关,第四双刀双掷电子开关通过RC电路连接第五双刀双掷电子开关,旋转电容滤波器电路中第二双刀双掷电子开关、第三双刀双掷电子开关、第四双刀双掷电子开关和第五双刀双掷电子开关在电路运行时相继循环连接。 
地址 110004 辽宁省沈阳市和平区文化路3号巷11号