首页 / 专利库 / 电脑编程 / 算法 / 期望最大化算法 / 基于时间序列分析的阀门故障实时诊断系统及诊断方法

基于时间序列分析的故障实时诊断系统及诊断方法

阅读:309发布:2020-05-19

专利汇可以提供基于时间序列分析的故障实时诊断系统及诊断方法专利检索,专利查询,专利分析的服务。并且本 发明 属于 阀 门 故障检测与诊断技术领域,具体涉及基于时间序列分析的阀门故障实时诊断系统及诊断方法。具体包括:通过振动 加速 度 传感器 和声 压传感器实时监测阀门;利用 信号 采集芯片将监测到的振动加速度信号和声压信号采集到诊断处理器中;对数据进行五点三次平滑降噪处理;生成不同故障诊断类型的诊断数据训练集;建立故障诊断时序模型;使用前向 算法 依次计算各个HMM模型的似然概率值,得出故障诊断结果;通过蓝牙模 块 将故障诊断结果传输到工业显示屏上进行显示。本发明通过时间序列分析对阀门进行故障诊断,在能在故障产生初期识别故障信息并得出故障类别,突破仅能在单一时刻观测得出故障的限制,同时降低故障误检率,显著提高阀门的使用性能。,下面是基于时间序列分析的故障实时诊断系统及诊断方法专利的具体信息内容。

1.一种基于时间序列分析的故障实时诊断系统,其特征在于:包括振动加速传感器和声压传感器,所述振动加速度传感器和声压传感器安装在待诊断的阀门上,用于从阀门上监测振动加速度信号和声压信号;信号采集芯片,通过电缆与所述振动加速度传感器和声压传感器连接,用于将振动加速度传感器和声压传感器监测到时的信号采集到诊断处理器中;诊断处理器,与所述信号采集芯片电连接,用于接收所述信号采集芯片的信号进行处理得出故障诊断类型,并通过蓝牙模传输到工业显示屏上进行显示。
2.根据权利要求1所述的基于时间序列分析的阀门故障实时诊断系统,其特征在于:所述故障诊断类型具体为:正常、早期故障、磨损故障和卡阀故障。
3.一种基于时间序列分析的阀门故障实时诊断方法,其特征在于:具体包括以下步骤:
A.通过振动加速度传感器和声压传感器实时监测振动加速度信号和声压信号;
B.利用信号采集芯片将监测到的振动加速度信号和声压信号采集到诊断处理器中;
C.利用诊断处理器对步骤B中获得的振动加速度信号和声压信号进行五点三次平滑降噪处理;
D.利用步骤C中获得的数据生成不同故障诊断类型的诊断数据训练集;
E.建立故障诊断时序模型;
F.使用前向算法依次计算各个HMM模型的似然概率值,得出阀门的故障诊断结果;
G.通过蓝牙模块将故障诊断结果传输到工业显示屏上进行显示。
4.根据权利要求3所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述步骤B中的采集步长为0.2s。
5.根据权利要求3所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述步骤C中五点三次平滑降噪处理为使用Savizkg-Golag五点三次滤波器对采集到的振动加速度a与声压τ进行平滑与降噪处理,具体包括以下步骤:
C.1、建立振动加速度a与采集时间t的三次项关系式
y=a0+a1x+a2x2+a3x3
其中,a0、a1、a2、a3为多项式各系数,y对应振动加速度a,x对应采集时间t;
C.2、设定动态时窗为1s,在每个动态时窗内采集5个点数据,分别为(x-2,y-2),(x-1,y-1),(x0,y0),(x1,y1),(x2,y2),将五个点的坐标一一代入,即有方程组
C.3、基于最小二乘法,方程组可转为
上列方程组可用矩阵表示为Y5×1=X5×4·A4×1+E5×1;
C.4、解得A的最小二乘解 则滤波后的值
从而对振动加速度a进行平滑与降噪处理,同理,对声压
τ进行平滑与降噪处理。
6.根据权利要求3所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述步骤D具体为:
D.1、对振动加速度a进行均匀量化编码,即对[amin,amax]分8个等分的区间,依次从0~7编码;将声压τ分为[0,0.1]、(0.1,+∞)两个区域,依次编码为I,II;
D.2、将振动加速度a、声压τ特征联合编码:将振动加速度特征a和声压τ特征结合,形成在t时刻同时包含两个特征的特征值,构成新的码字:I0,II0,I1,…,II6,I7,II7;
D.3、将步骤D.2处理后的数据根据故障诊断类型分为正常、早期故障、磨损故障和卡阀故障4类数据,每条数据由20连续步长的序列和其对应故障诊断类型标签组成,序列矩阵为
7.根据权利要求3所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述步骤E具体为:
E.1、针对每类故障诊断类型分别建立HMM(隐尔可夫模型,Hidden Markov Model,简称HMM)模型,每个HMM模型由一个五元组μ=(Q,V,A,B,π)组成,其中,隐藏状态Q={Q1,Q2,…,QN},N为隐藏状态的数,其中N为4;可观察状态V={V1,V2,…,VM}为振动加速度a与声压τ,M为观察状态的数目,M为16;隐藏状态转移概率矩阵A=[aij]N×N的元素表示HMM模型中各个隐藏状态之间的转移概率,aij是在t时刻隐藏状态为Qi、在t+1时刻隐藏状态为Qj的概率,aij=P(It+1=Qj|It=Qi),i=1,2…,N;j=1,2…,N,I是长度为T的状态序列,且I={I1,I2,…,IT};混淆矩阵B=[bj(k)]N×M的元素表示HMM模型中各个隐藏状态和观察状态之间的转移概率,bj(k)表示在t时刻,隐藏状态为Qj、观察状态为Ot的概率,bj(k)=P(Ot=Vk|It=Qj),k=1,2…,M;j=1,2…,N,O是对应的观测序列;初始状态概率矩阵π=(πi),其中πi=P(I1=Qi),i=1,2,…,N,表示初始时刻t=1各个隐藏状态Qi的概率;
E.2、对每类HMM模型初始化:随机的给参数πi,aij,bj(k)赋值,使其满足约束:
由此得到模型μ0;
E.3、将同一类别故障诊断类型的观测序列 作为对应的HMM分类模型的输入,依据模型初始化后的参数,采用EM期望最大化算法(Expectation Maximum)调整模型μ的参数,使概率函数 最大化即 逐步更新模型参数,最终获取各个故障诊
断类型对应的最优HMM模型。
8.根据权利要求7所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述HMM模型具体包括阀门正常HMM模型、阀门早期故障HMM模型、阀门磨损故障HMM模型和阀门卡钳故障HMM模型。
9.根据权利要求3所述的基于时间序列分析的阀门故障实时诊断方法,其特征在于:所述步骤F具体为:
α1(i)=πibi(O1),1≤i≤N
其中αt(i)为前向中间变量,表示在时间t、HMM模型输出了序列O1O2ΛOt、并且位于状态si的概率,最终取似然概率值最大者为阀门的故障诊断结果。

说明书全文

基于时间序列分析的故障实时诊断系统及诊断方法

技术领域

[0001] 本发明属于阀门故障检测与诊断技术领域,具体涉及基于时间序列分析的阀门故障实时诊断系统及诊断方法。

背景技术

[0002] 在大多数工业领域,阀门是必不可少的元件之一,若阀门发生严重故障,则整个运作系统将无法正常工作,易引发严重的工作事故。时间序列是将某种统计指标的数值,按时间先后顺序排序所形成的数列。时间序列的识别与预测就是通过分析时间序列,根据时间序列所反映出来的发展过程、方向和趋势,进行类推或延伸,识别或预测下一段时间或以后若干年内可能达到的平。时间序列的异常检测就是通过历史的数据分析,查看当前的数据是否发生了明显偏离了正常的情况。
[0003] 传统的阀门故障诊断机理,是以检修工人以一定周期地去巡检与查验。随着工业技术的不断发展,经常拆卸阀门的诊断方法已经远不能适应要求,并且增加了维修成本和检修周期。而且,传统的常规手段在阀门故障初期难以通过直接监测,而检测出早期的故障信号对提升整个系统的可靠性才有意义。阀门在实际运作中,故障可以通过内部工况及外界因素去加以观测。而随着现代信号分析与处理技术的发展,这些观测数据可以准确地获取,这为阀门故障的时序分析奠定了基础。同时现有的方法,仅仅能够检测阀门是否有故障,并不能得出故障的类别,存在较大的局限性,不方便后期的维修。

发明内容

[0004] 本发明要解决的技术问题在于克服现有技术的不足,提供一种基于时间序列分析的阀门故障实时诊断系统及诊断方法,通过时间序列分析对阀门进行故障诊断,在能在故障产生初期识别故障信息并得出故障类别,突破仅能在单一时刻观测得出故障的限制,同时降低故障误检率,显著提高阀门的使用性能。
[0005] 本发明解决上述技术问题的技术方案如下:基于时间序列分析的阀门故障实时诊断系统,包括振动加速传感器和声压传感器,所述振动加速度传感器和声压传感器安装在电路板上,用于从阀门上监测振动加速度信号和声压信号;信号采集芯片,通过电缆与所述电路板连接,用于将振动加速度传感器和声压传感器监测到时的信号采集到诊断处理器中;诊断处理器,与所述信号采集芯片电连接,用于接收所述信号采集芯片的信号进行处理得出故障诊断类型,并通过蓝牙模传输到工业显示屏上进行显示。
[0006] 进一步地,所述故障诊断类型具体为:正常、早期故障、磨损故障和卡阀故障。
[0007] 本发明还提供了一种基于时间序列分析的阀门故障实时诊断方法,具体包括以下步骤:
[0008] H.通过振动加速度传感器和声压传感器实时监测振动加速度信号和声压信号;
[0009] I.利用信号采集芯片将监测到的振动加速度信号和声压信号采集到诊断处理器中;
[0010] J.利用诊断处理器对步骤B中获得的振动加速度信号和声压信号进行五点三次平滑降噪处理;
[0011] K.利用步骤C中获得的数据生成不同故障诊断类型的诊断数据训练集;
[0012] L.建立故障诊断时序模型;
[0013] M.使用前向算法依次计算各个HMM模型的似然概率值,得出阀门的故障诊断结果;
[0014] N.通过蓝牙模块将故障诊断结果传输到工业显示屏上进行显示。
[0015] 进一步地,所述步骤B中的采集步长为0.2s。
[0016] 进一步地,所述步骤C中五点三次平滑降噪处理为使用Savizkg-Golag五点三次滤波器对采集到的振动加速度a与声压τ进行平滑与降噪处理,具体包括以下步骤:
[0017] C.1、建立振动加速度a与采集时间t的三次项关系式
[0018] y=a0+a1x+a2x2+a3x3
[0019] 其中,a0、a1、a2、a3为多项式各系数,y对应振动加速度a,x对应采集时间t;
[0020] C.2、设定动态时窗为1s,在每个动态时窗内采集5个点数据,分别为(x-2,y-2),(x-1,y-1),(x0,y0),(x1,y1),(x2,y2),将五个点的坐标一一代入,即有方程组
[0021]
[0022] C.3、基于最小二乘法,方程组可转为
[0023]
[0024] 上列方程组可用矩阵表示为Y5×1=X5×4·A4×1+E5×1;
[0025] C.4、解得A的最小二乘解 则滤波后的值
[0026] 从而对振动加速度a进行平滑与降噪处理,同理,对声压τ进行平滑与降噪处理。
[0027] 进一步地,所述步骤D具体为:
[0028] D.1、对振动加速度a进行均匀量化编码,即对[amin,amax]分8个等分的区间,依次从0~7编码;将声压τ分为[0,0.1]、(0.1,+∞)两个区域,依次编码为I,II;
[0029] D.2、将振动加速度a、声压τ特征联合编码:将振动加速度特征a和声压τ特征结合,形成在t时刻同时包含两个特征的特征值,构成新的码字:I0,II0,I1,…,II6,I7,II7;
[0030] D.3、将步骤D.2处理后的数据根据故障诊断类型分为正常、早期故障、磨损故障和卡阀故障4类数据,每条数据由20连续步长的序列和其对应故障诊断类型标签组成,序列矩阵为
[0031] 进一步地,所述步骤E具体为:
[0032] E.1、针对每类故障诊断类型分别建立HMM(隐尔可夫模型,Hidden Markov Model,简称HMM)模型,每个HMM模型由一个五元组μ=(Q,V,A,B,π)组成,其中,隐藏状态Q={Q1,Q2,…,QN},N为隐藏状态的数,其中N为4;可观察状态V={V1,V2,…,VM}为振动加速度a与声压τ,M为观察状态的数目,M为16;隐藏状态转移概率矩阵A=[aij]N×N的元素表示HMM模型中各个隐藏状态之间的转移概率,aij是在t时刻隐藏状态为Qi、在t+1时刻隐藏状态为Qj的概率,aij=P(It+1=Qj|It=Qi),i=1,2…,N;j=1,2…,N,I是长度为T的状态序列,且I={I1,I2,…,IT};混淆矩阵B=[bj(k)]N×M的元素表示HMM模型中各个隐藏状态和观察状态之间的转移概率,bj(k)表示在t时刻,隐藏状态为Qj、观察状态为Ot的概率,bj(k)=P(Ot=Vk|It=Qj),k=1,2…,M;j=1,2…,N,O是对应的观测序列;初始状态概率矩阵π=(πi),其中πi=P(I1=Qi),i=1,2,…,N,表示初始时刻t=1各个隐藏状态Qi的概率;
[0033] E.2、对每类HMM模型初始化:随机的给参数πi,aij,bj(k)赋值,使其满足约束:由此得到模型μ0;
[0034] E.3、将同一类别故障诊断类型的观测序列 作为对应的HMM分类模型的输入,依据模型初始化后的参数,采用EM期望最大化算法(Expectation Maximum)调整模型μ的参数,使概率函数 最大化即 逐步更新模型参数,最终获取各个故障诊断类型对应的最优HMM模型。
[0035] 进一步地,所述HMM模型具体包括阀门正常HMM模型、阀门早期故障HMM模型、阀门磨损故障HMM模型和阀门卡钳故障HMM模型。
[0036] 进一步地,所述步骤F具体为:
[0037] α1(i)=πibi(O1),1≤i≤N
[0038]
[0039]
[0040] 其中αt(i)为前向中间变量,表示在时间t、HMM模型输出了序列O1O2ΛOt、并且位于状态si的概率,最终取似然概率值最大者为阀门的故障诊断结果。
[0041] 本发明的有益效果是:
[0042] 1、本发明通过建立时序分析模型,对阀门有无故障进行检测,实时性与智能性好,可适应阀门在实际使用中的稳定性与可靠性的要求,减免了人工检测的大量工作;
[0043] 2、本方法不仅能检测出阀门是否有故障,还能得出在存在故障时给出故障类型,方便了后期的维修与整改,实用性高;
[0044] 3、本发明采用了振动加速度与声压双特征联立编码,提高了阀门故障检测的准确性,同时简化了特征值,减免了时序模型的计算复杂度,大大省去了模型训练时间;
[0045] 4、本发明突破了传统的常规手段难以在阀门故障早期获取故障信息的限制,提升了整个阀门系统的可靠性,降低了维修难度。附图说明
[0046] 为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
[0047] 图1为本发明具体实施例所述的基于时间序列分析的阀门故障实时诊断系统的结构框图
[0048] 图2为本发明具体实施例所述的基于时间序列分析的阀门故障实时诊断方法的流程示意图;

具体实施方式

[0049] 下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只作为示例,而不能以此来限制本发明的保护范围。
[0050] 在本申请中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
[0051] 实施例
[0052] 如图1所示,本发明所提供的基于时间序列分析的阀门故障实时诊断系统,基于时间序列分析的阀门故障实时诊断系统,包括振动加速度传感器和声压传感器,所述振动加速度传感器和声压传感器集装在电路板上,电路板安装在待检测的阀门上,用于从阀门上监测振动加速度信号和声压信号;信号采集芯片,通过电缆与所述电路板连接,用于将振动加速度传感器和声压传感器监测到时的信号采集到诊断处理器中;诊断处理器,与所述信号采集芯片电连接,用于接收所述信号采集芯片的信号进行处理得出故障诊断类型,并通过蓝牙模块传输到工业显示屏上进行显示。
[0053] 进一步地,所述故障诊断类型具体为:正常、早期故障、磨损故障和卡阀故障。
[0054] 本发明还提供了一种基于时间序列分析的阀门故障实时诊断方法,如图2所示,具体包括以下步骤:
[0055] 步骤1:通过振动加速度传感器和声压传感器实时监测振动加速度a和声压τ,并通过信号采集芯片采集到故障诊断处理器中,其中采集步长为0.2s。
[0056] 步骤2:数据处理
[0057] 传感器收取的数据存在较大的噪声,使用Savizkg-Golag五点三次滤波器对采集到的阀门的振动加速度a与声压τ进行平滑与降噪处理,具体包括以下步骤:
[0058] 将采集到的振动加速度a与采集时间t一一对应,有关系式y=f(x),其中,y对应振动加速度a,x对应采集时间t,设定动态时窗为1s,每个动态时窗包含5个点数据,分别为(x-2,y-2),(x-1,y-1),(x0,y0),(x1,y1),(x2,y2),f(x)采用三次多项式对振动加速度a拟合,即有
[0059] y=a0+a1x+a2x2+a3x3
[0060] 其中,a0、a1、a2、a3为多项式各系数,将五个点的坐标一一代入,即有方程组[0061]
[0062] 基于最小二乘法,方程组可转为
[0063]
[0064] 上列方程组可用矩阵表示为Y5×1=X5×4·A4×1+E5×1;
[0065] 解得A的最小二乘解 则滤波后的值
[0066] 从而对振动加速度a进行平滑与降噪处理,同理,对声压τ进行平滑与降噪处理。
[0067] 步骤3:生成诊断数据训练集
[0068] 对振动加速度a进行均匀量化编码,即对[amin,amax]分8个等分的区间,依次从0~7编码;将声压τ分为[0,0.1]、(0.1,+∞)两个区域,单位Pa,依次编码为I,II;
[0069] 将振动加速度a、声压τ特征联合编码:将振动加速度特征a和声压τ特征结合,形成在t时刻同时包含两个特征的特征值,构成新的码字:I0,II0,I1,…,II6,I7,II7;
[0070] 根据阀门的工作特点,将处理后的数据根据故障诊断类型分为正常、早期故障、磨损故障和卡阀故障4类数据,每条数据由20连续步长的序列和其对应故障诊断类型标签组成,序列矩阵为
[0071] 步骤4:建立故障诊断时序模型
[0072] 本发明是基于隐马尔可夫模型故障诊断时序模型,隐马尔科夫模型(Hidden Markov Model,HMM)是关于时序的概率模型,通过可观测到的变量去研究不可观测的变量。针对每类故障诊断类型分别建立HMM模型,分别为阀门正常HMM模型、阀门早期故障HMM模型、阀门磨损故障HMM模型、阀门卡钳故障HMM模型,每个HMM模型由一个五元组μ=(Q,V,A,B,π)组成,其中,隐藏状态Q={Q1,Q2,…,QN},N为隐藏状态的数,其中N为4;可观察状态V={V1,V2,…,VM}为振动加速度a与声压τ,M为观察状态的数目,M为16;隐藏状态转移概率矩阵A=[aij]N×N的元素表示HMM模型中各个隐藏状态之间的转移概率,aij是在t时刻隐藏状态为Qi、在t+1时刻隐藏状态为Qi的概率,aij=P(It+1=Qj|It=Qi),i=1,2…,N;j=1,2…,N,I是长度为T的状态序列,且I={I1,I2,…,IT};混淆矩阵B=[bj(k)]N×M的元素表示HMM模型中各个隐藏状态和观察状态之间的转移概率,bj(k)表示在t时刻,隐藏状态为Qj、观察状态为Ot的概率,bj(k)=P(Ot=Vk|It=Qj),k=1,2…,M;j=1,2…,N,O是对应的观测序列;初始状态概率矩阵π=(πi),其中πi=P(I1=Qi),i=1,2,…,N,表示初始时刻t=1各个隐藏状态Qi的概率;
[0073] 对每类HMM模型初始化:随机的给参数πi,aij,bj(k)赋值,使其满足约束:由此得到模型μ0;
[0074] 将同一类别故障诊断类型的观测序列 作为对应的HMM分类模型的输入,依据模型初始化后的参数,采用EM期望最大化算法(Expectation Maximum)调整模型μ的参数,使概率函数 最大化即 逐步更新模型参数,最终获取各个故障诊断类型对应的最优HMM模型。
[0075] 步骤5:阀门故障实时诊断
[0076] 实时采集振动加速度与声压传感器的数据,经过数据处理后,使用前向算法依次计算各个HMM模型的似然概率值:
[0077] α1(i)=πibi(O1),1≤i≤N
[0078]
[0079]
[0080] 其中αt(i)为前向中间变量,表示在时间t,HMM输出了序列O1O2ΛOt,并且位于状态si的概率,最终利用贝叶斯判别法对每个HMM模型的似然概率值进行判断,取似然概率值最大者为阀门的故障诊断结果。
[0081] 步骤6:通过蓝牙模块将故障诊断结果传输到工业显示屏上进行显示。
[0082] 最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈