CN103824137B - 一种复杂机械设备多工况故障预测方法 - Google Patents
一种复杂机械设备多工况故障预测方法 Download PDFInfo
- Publication number
- CN103824137B CN103824137B CN201410103676.3A CN201410103676A CN103824137B CN 103824137 B CN103824137 B CN 103824137B CN 201410103676 A CN201410103676 A CN 201410103676A CN 103824137 B CN103824137 B CN 103824137B
- Authority
- CN
- China
- Prior art keywords
- fault
- spe
- value
- amplitude
- representing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 230000008569 process Effects 0.000 claims abstract description 45
- 238000001514 detection method Methods 0.000 claims abstract description 29
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 42
- 230000007704 transition Effects 0.000 claims description 33
- 238000012706 support-vector machine Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 abstract 3
- 238000000513 principal component analysis Methods 0.000 description 21
- 238000004088 simulation Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013499 data model Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 101000620653 Homo sapiens Serine/threonine-protein phosphatase 5 Proteins 0.000 description 1
- 102100022346 Serine/threonine-protein phosphatase 5 Human genes 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明涉及一种复杂机械设备多工况故障预测方法,其步骤:1)针对多工况过程建立多PCA模型,对每个PCA模型计算出相应的检测指标T2统计量和SPE;2)对多PCA模型中的各个T2统计量和SPE两检测指标进行优化,并对机械设备进行故障检测,检测得到过渡过程机械设备的故障数据;3)对优化后的两检测指标T2统计量和SPE检测到的过渡过程机械设备故障数据进行故障重构,得到使重构后SPE最小化的幅值估计值;4)对过渡过程中不同工况下的同一故障重构后的幅值估计值进行一致性幅值估计;5)根据进行一致性幅值估计后的幅值估计值
Description
技术领域
本发明涉及一种机械设备故障预测方法,特别是关于一种复杂机械设备多工况故障预测方法。
背景技术
随着科学技术和工业发展,机械设备向着大型化、高速化、复杂化发展。因此,现在企业生产对设备及***的可靠性、连续性、经济性等要求日益提高,在以往对设备及***故障进行有效诊断和提出解决方案的基础上,进一步要求,在故障只出现微小异常征兆时即可实现对故障进行预报并提出相应紧急处理措施。故障预测的方法有多种多样,其中统计过程监控技术已经发展20余年,并且广泛应用于工业过程的故障检测、诊断和估计等。在最近的研究中,基于主元分析(Principal Component Analysis,PCA)的故障估计方法被成功用于故障预测中,但对于多工况过程的数据,基于主元分析的故障估计方法却不能够非常准确地进行故障预测,当复杂***运行在多个工况时,变量之间的关系会根据***当前的工作模式发生相应的变化,如果按照单工况的数据模型描述***,会造成大量的误报和漏报。因此需要针对多工况过程的数据提出新的故障预测方案。
发明内容
针对上述问题,本发明的目的是提供一种复杂机械设备多工况故障预测方法,该方法的预测精度较高,能有效避免故障误报和漏报情况,避免了不同工况下同一种故障的估计误差。
为实现上述目的,本发明采取以下技术方案:一种复杂机械设备多工况故障预测方法,其包括以下步骤:1)针对多工况过程建立多PCA模型,并针对每个PCA模型计算出其相应的检测指标Hotelling’s T2统计量和SPE,其中Hotelling’s T2统计量简称T2统计量;(1)假设x∈Rm表示具有m个测量变量的样本向量,正常运行时的样本有n个;数据矩阵X∈Rn×m由n个样本组成,其中每一行代表了一个样本,每一列代表了一个测量变量共有n次采样;将数据矩阵X各列经过标准化处理成零均值和单位方差的变量,得到进行标准化后的样本x的协方差矩阵S,并对该协方差矩阵特征值分解并按大小降序排列;协方差矩阵S为:
(2)PCA模型将样本矩阵X∈Rn×m分解成建模部分和残差部分E两个部分:,其中,表示被建模部分;E表示残差部分;P∈Rm×A为负载矩阵,是由S的前A个特征向量组成的,A表示主元的个数;T∈Rn×A为得分矩阵,T=XP;(3)根据已有的测量变量在每个稳态工况下的历史数据建立相对应的主元模型,组成包含所有工况的多主元模型:,其中,q表示稳定工况个数;(4)在多PCA模型中,针对每个PCA模型需要计算出其相应的检测指标T2统计量和SPE,用SPE指标衡量样本向量在残差空间投影的变化,用T2统计量衡量测量变量在主元空间中的变化;2)采用随时间加权算法对多PCA模型中的各个T2统计量和SPE两检测指标进行优化,并根据优化后的检测指标T2统计量和SPE对机械设备进行故障检测,检测得到过渡过程机械设备的故障数据:(1)将多工况过渡过程中的上一工况的权重w1分和下一工况的权重w2分别取为:
式中,t1表示过渡过程开始的时刻,t2表示过渡过程结束的时刻,Δt=t2-t1;w1+w2=1;(2)假设t=k时刻,分别表示上一工况和下一工况的测量变量均值,σ1和σ2分别表示上一工况和下一工况的测量变量标准差,则此时过渡过程的均值和标准差σt=k分别为:
σt=k=w1×σ1+w2×σ2;
(3)样本向量通过加权后的均值和标准差σt=k进行预处理,求出样本向量的协方差矩阵St=k和过渡过程的负载矩阵Pt=k,进而得到优化后的SPE指标和T2统计量:
3)对步骤2)中优化后的两检测指标T2统计量和SPE检测到的过渡过程机械设备故障数据进行故障重构,得到使重构后SPE最小化的幅值估计值;4)对过渡过程中不同工况下的同一故障重构后的幅值估计值进行一致性幅值估计:(1)假设,式中,表示同一故障在上一工况和下一工况下的幅值估计值,Ξi_1、Ξi_2表示各自故障在残差空间的投影矩阵,假设维数为a;(2)根据步骤(1)得到同一故障在不同工况下的幅值估计值为:求解使最小化的d和e,得到同一故障在不同工况下的一致估计5)根据进行一致性幅值估计后的幅值估计值,利用支持向量机预测模型对故障幅值进行趋势预测。
所述步骤3)中,所述故障重构的过程如下:(1)假设检测得到的故障数据用x表示,正常数据用x*表示,x*,x∈[x1,x2,...,xm],m为每一个测量样本的维数,故障幅值的大小用f表示,则:x=x*+Ξf,式中,Ξ表示故障子空间;(2)重构正常数据以消除故障的影响,得到正常数据的估计值为:xx *=x-Ξifi,式中,fi表示对故障幅值大小f的估计;(3)根据步骤(2)重构后的正常数据的估计值得到使重构后SPE最小化的幅值估计值,式中,表示故障子空间Ξ到残差子空间的投影,C表示残差子空间。
所述步骤5)中,所述支持向量机预测方法如下:(1)对输入的测量变量采用降噪平滑处理进行数据的预处理,根据FPE准则计算嵌入维数,然后分别形成训练和测试样本;(2)对支持向量机模型进行初始化,将拉格朗日乘子αi和以及阈值b赋以随机初始值;(3)利用训练样本建立目标函数,得到拉格朗日乘子αi和以及阈值b,最后选取最优的惩罚因子C*和核函数参数g;(4)将步骤(3)得到的参数值代入预测模型,用测试样本计算在未来某时刻的预测值;(5)计算误差函数,当误差的绝对值小于预先设定值时,结束学习过程,否则返回步骤(1)继续学习。
所述步骤5)中,所述支持向量机预测模型对故障幅值进行趋势预测精度采用均方误差和平均相对预测误差来评价。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于采用基于多PCA故障重构的方法对多工况过程的故障进行预测,可以很好地解决过程数据的多工况问题,从隐含故障的数据中挖掘出故障方向并估计出故障幅值,也考虑到了故障的多维特性,可以获得更加精确的故障预测结果。2、由于复杂机械设备运行在多个工况时,测量变量之间的关系会根据设备当前的工作模式发生相应的变化。如果按照单工况的数据模型描述设备状态,会造成大量的误报和漏报。因此,本发明从过程监控的角度采用多主元分析模型,利用过渡过程的加权检测方法,提出了基于多主元分析模型的重构、故障估计和故障一致性估计,避免了不同工况下同一种故障的估计误差。3、本发明采用基于重构的多工况故障预测方法,当故障幅值被估计出来以后,采用支持向量机对故障幅值进行趋势预测。由于故障重构的目的是沿故障的方向用故障数据构造正常数据,消除故障对正常数据的影响,因此故障重构后进行预测的结果更为准确。4、本发明采用支持向量机对故障幅值进行直接预测,解决了传统预测方法无法利用多维有效数据的问题,在长期预测中,支持向量机模型具有很好的预测精度,具有很好的泛化性能。本发明可以广泛在机电设备故障预测中应用。
附图说明
图1是本发明在多PCA模型下的SPE指标检测仿真曲线示意图;
图2是本发明在多PCA模型下的T2统计量指标检测仿真曲线示意图;
图3是本发明的两个工况下的分别故障估计曲线示意图;
图4是本发明的两个工况下的分别故障一致估计曲线示意图;
图5是本发明基于重构后的故障预测曲线示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供一种机械设备多工况故障预测方法,其包括以下步骤:
1)针对多工况过程建立多PCA模型(主元分析模型),并针对每个PCA模型计算出其相应的检测指标Hotelling’s T2统计量(以下简称T2统计量)和SPE(平方预测误差,也称Q统计量);
(1)假设x∈Rm表示具有m个测量变量的样本向量(即m为样本x的维数),正常运行时的样本有n个。数据矩阵X∈Rn×m由n个样本组成,其中每一行代表了一个样本,每一列代表了一个测量变量共有n次采样。将数据矩阵X各列经过标准化处理成零均值和单位方差的变量,可以得到进行标准化后的样本x的协方差矩阵S,并对该协方差矩阵特征值分解并按大小降序排列。协方差矩阵S为:
其中,将数据矩阵X各列经过标准化处理成零均值和单位方差的变量的方法是将数据矩阵X的每一列减去相应的变量均值并且除以相应的变量标准差。
(2)根据PCA模型将测量变量空间分成主元子空间和残差子空间两个正交且互补的子空间,任意一个样本向量均可分解成为在主元子空间和残差子空间上的投影,即PCA模型将样本矩阵X∈Rn×m分解成建模部分和残差部分E两个部分:
其中,表示被建模部分;E表示残差部分;P∈Rm×A为负载矩阵,是由S的前A个特征向量组成的,A表示主元的个数;T∈Rn×A为得分矩阵,T=XP。
(3)根据已有的测量变量在每个稳态工况下的历史数据建立相对应的主元模型,组成包含所有工况的多主元模型:
其中,q表示稳定工况个数。
(4)在多PCA模型中,针对每个PCA模型需要计算出其相应的检测指标T2统计量和SPE,即T2(q)和SPE(q),用SPE指标衡量样本向量在残差空间投影的变化,用T2统计量衡量测量变量在主元空间中的变化:
其中,SPE指标表达式为:
式中,I为单位矩阵;表示置信水平为α时SPE的控制限。当SPE在控制限内时,认为当前运转过程处于正常状态。当SPE值超出了统计控制限时,代表当前运转过程发生了故障,SPE值的变化代表着数据间相关性的变化。该控制限的计算公式为:
式中,,λj为样本矩阵X的协方差矩阵∑的特征值,cα为标准正态分布在置信水平α下的阈值,m是样本x的维数。
T2统计量表达式为:
其中,Λ=diag{λ1,...,λA},表示置信度为α的T2统计限。当T2位于控制限内时,认为当前运转过程处于正常工作状态。
由此可知,主元分析是根据数据变化的方差大小来确定变化方向的主次地位,按主次顺序得到各个主元,这些主元彼此是相互独立的。本发明采用SPE和T2统计量两种指标来检测机械设备运转过程是否发生异常,T2统计量指标衡量的是样本距离主元子空间原点的距离,而SPE指标主要体现了正常过程变量之间的相关性被改变的程度,显示了异常的过程状态。主元子空间主要反映正常数据变化的测试,残差子空间主要反映噪声情况,因此T2统计量定义的控制限比SPE定义的要大得多,较小的故障可能不会超过T2统计量控制限但却可以超过SPE控制限。假设在正常的状态下,各个过程变量之间的关系是不变的,***的动态变化主要表现在主元子空间中。当有故障发生时,残差子空间将会发生变化。因此采用SPE指标来检测残差空间中是否发生异常,T2统计量衡量测量变量在主元空间中的变化。
2)采用随时间加权算法对多PCA模型中的各个T2统计量和SPE两检测指标进行优化,并根据优化后的检测指标T2统计量和SPE对机械设备进行故障检测,检测得到过渡过程机械设备的故障数据,通过优化后的检测指标进行故障检测可以有效避免在工况过渡过程中出现故障误报;
其中,对T2统计量和SPE两检测指标进行优化的方法包括以下步骤:
(1)由于在多工况的过渡过程中每一个样本向量点的权重值都是不同的,都会随时间逐渐增加,上一工况的权重值从1逐渐到0,下一工况的权重值从0逐渐到1。因此,将多工况过渡过程中的上一工况的权重w1分和下一工况的权重w2分别取为:
式中,t1表示过渡过程开始的时刻,t2表示过渡过程结束的时刻,Δt=t2-t1;w1+w2=1。
当一个稳定运行的工况刚刚向下一个工况进行过渡时,过渡过程的特性受该工况的影响较大,随着时间的推移,越靠近新的工况的时间点上,过渡过程的特性越贴近下一个工况,直到进入下一工况,因此过渡过程中的每一个点的权重值都是不同的。
(2)假设t=k时刻,分别表示上一工况和下一工况的测量变量均值,σ1和σ2分别表示上一工况和下一工况的测量变量标准差,则此时过渡过程的均值和标准差σt=k分别为:
σt=k=w1×σ1+w2×σ2。 (10)
(3)样本向量通过步骤(2)中求出的加权后的均值和标准差σt=k进行预处理,求出样本向量的协方差矩阵St=k和过渡过程的负载矩阵Pt=k,进而得到优化后的SPE指标和T2统计量:
如图1所示,SPE指标的检测仿真曲线,横坐标表示样本数,纵坐标表示SPE值。虚线代表各个统计量99%的控制限,若超过控制限则表示***出现故障。由图可知,基于多PCA方法的SPE在第166个样本点超过控制限。
如图2所示,T2统计量指标的检测仿真曲线,横坐标表示样本数,纵坐标表示T2值。虚线代表各个统计量99%的控制限,若超过控制限则表示***出现故障。由图可知,基于多PCA方法的T2在第164个样本点超过控制限。
3)对步骤2)中优化后的两检测指标T2统计量和SPE检测到的过渡过程机械设备故障数据进行故障重构,得到使重构后SPE最小化的幅值估计值,以使故障因素对故障数据内正常部分的影响降至最小,进而构造出正常状态下的过程数据。
故障重构的具体过程如下:
(1)假设检测得到的故障数据用x表示(即x为样本),正常数据用x*表示,x*,x∈[x1,x2,...,xm],m代表传感的个数,即每一个测量样本的维数,故障幅值的大小用f表示,则:
x=x*+Ξf, (13)
式中,Ξ表示故障子空间(也称为故障方向矩阵)。
(2)重构正常数据以消除故障的影响,得到正常数据的估计值为:
xi *=x-Ξifi (14)
式中,fi表示对故障幅值大小f的估计,从几何意义上讲,是将样本x沿着故障所在的子空间拉回到主元子空间。
(3)根据步骤(2)重构后的正常数据的估计值,得到使重构后SPE最小化的幅值估计值
式中,,表示故障子空间Ξ到残差子空间的投影,C表示残差子空间。
其中,重构后的SPE指标为:
重构就是根据公式(16)搜索,使得:
式(17)的最优解为:
其中,表示的Moore-Penrose伪逆,表示使重构后SPE最小化的故障幅值的估计。
4)对过渡过程中不同工况下的同一故障重构后的幅值估计值进行一致性幅值估计:当故障出现在工况一的过程中,而且延续到了工况二的过程中,这时的故障方向矩阵Ξ发生了变化,但是故障的幅值没有发生变化。由于是同一个故障,即使在不同的工况下被重构估计,理论上也应该会得到一致估计,可是在实际工业过程中,数据的噪声和不同工况下机器的干扰等因素都会使结果产生偏差,因此需要对不同工况下的同一故障进行一致性幅值估计。
(1)假设
式中,表示同一故障在上一工况和下一工况下的幅值估计值(即不同工况下的幅值估计值),Ξi_1、Ξi_2表示各自故障在残差空间的投影矩阵,假设维数为a。
(2)根据公式(18)、(19)可以得到同一故障在不同工况下的幅值估计值为:
求解使最小化的d和e,其求解方法如下:
令矩阵G=[A,B],G∈Rn×2k对GT进行奇异值分解(SVD):
UHGTV=Σ(22)
其中,U和V都是酉矩阵维数分别为2k和n,Σ是对角阵,对角线上元素为奇异值。按照由大到小对奇异值进行排序,取左奇异向量U的后a列特征向量,此时阶数是2k×k,前k行即为d的估计,后k行即为e的估计,得到d和e后代入公式(20)和公式(21)中,得到的即为同一故障在不同工况下的一致估计。
如图3所示,假设工况一先发生故障,然后转换工况,故障进入工况二。如图所示,样本点从0到900为工况一,故障是在样本点160处发生,样本点900到1100是过渡过程,从样本点1100开始进入工况二。故障估计的曲线理论上应该是在工况一和工况二时误差较小,但如图3所示,工况一的故障估计和工况二的故障估计误差较大。
如图4所示,假设工况一先发生故障,然后转换工况,故障进入工况二。如图所示,样本点从0到900为工况一,故障是在样本点160处发生,样本点900到1100是过渡过程,从样本点1100开始进入工况二。经过一致幅值估计的算法后,得到的故障估计如图4所示,这时的故障估计达到一致,经过过渡过程后曲线更平滑。
5)根据进行一致性幅值估计后的幅值估计值,利用支持向量机(SVM)预测模型对故障幅值进行趋势预测。其预测方法如下:
(1)对输入的测量变量采用降噪平滑处理进行数据的预处理,根据FPE准则计算嵌入维数,然后分别形成训练和测试样本。
(2)对支持向量机模型进行初始化,将拉格朗日乘子αi和以及阈值b赋以随机初始值。
(3)利用训练样本建立目标函数,得到拉格朗日乘子αi和以及阈值b,最后选取最优的惩罚因子C*和核函数参数g。
(4)将步骤(3)得到的参数值代入预测模型,用测试样本计算在未来某时刻的预测值。
(5)计算误差函数,当误差的绝对值小于预先设定值时,结束学习过程,否则返回步骤(1)继续学习。
例如,对于给定的时间序列{x1,x2,...,xn},i=1,2,...,n,xn是预测的目标值,将先前的{xn-1,xn-2,...,xn-m}作为输入值,建立自相关输入与输出yn={xn}之间的映射关系f:,m为嵌入维数,经过变换之后,得到用于支持向量机的学习样本:
对向量机进行训练的回归函数为:
其中,αi和是拉格朗日乘子,b是阈值,K(xi,xt)是核函数,本发明采用径向基核函数(Radial Basis Function):
则一步预测模型为:
由此得到一个样本,表示第n+1个数据的预测值。
进而有
式(27)表示第p步的预测模型,得xn表示第n个数据的实际值,表示第n个数据的预测值。
上述步骤5)中,SVM预测模型对故障幅值进行趋势预测精度可以采用均方误差(MSE)和平均相对预测误差(ARE)来评价,MSE和ARE越小预测精度越精确,本发明优选将MSE和ARE控制在5%以内。其中,均方误差和平均相对预测误差分别为:
均方误差为:
平均相对预测误差为:
例如,如图5所示,选取2000个样本数据用于故障检测,故障在第166个样本点被SPE指标检测出来。当故障被检测出来以后,利用故障重构对故障幅值进行估计,得到相应的故障估计。采用重构后得到的故障估计进行故障预测,利用支持向量机预测模型对未来一步和十步的故障估计进行预测。图5中实线表示实际曲线,点划线表示一步预测曲线,星型线表示十步预测曲线,则一步预测和十步预测的相对预测误差如表1所示。
表1不同预测步数下的相对预测误差
上述各实施例仅用于说明本发明,各部件的连接和结构都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别部件的连接和结构进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (4)
1.一种复杂机械设备多工况故障预测方法,其包括以下步骤:
1)针对多工况过程建立多PCA模型,并针对每个PCA模型计算出其相应的检测指标Hotelling’s T2统计量和SPE,其中Hotelling’s T2统计量简称T2统计量;
(1)假设x∈Rm表示具有m个测量变量的样本向量,正常运行时的样本有n个;数据矩阵X∈Rn×m由n个样本组成,其中每一行代表了一个样本,每一列代表了一个测量变量共有n次采样;将数据矩阵X各列经过标准化处理成零均值和单位方差的变量,得到进行标准化后的样本x的协方差矩阵S,并对该协方差矩阵特征值分解并按大小降序排列;协方差矩阵S为:
(2)PCA模型将样本矩阵X∈Rn×m分解成建模部分和残差部分E两个部分:
其中,表示被建模部分;E表示残差部分;P∈Rm×A为负载矩阵,是由S的前A个特征向量组成的,A表示主元的个数;T∈Rn×A为得分矩阵,T=XP;
(3)根据已有的测量变量在每个稳态工况下的历史数据建立相对应的主元模型,组成包含所有工况的多主元模型:
其中,q表示稳定工况个数;
(4)在多PCA模型中,针对每个PCA模型需要计算出其相应的检测指标T2统计量和SPE,用SPE指标衡量样本向量在残差空间投影的变化,用T2统计量衡量测量变量在主元空间中的变化;
2)采用随时间加权算法对多PCA模型中的各个T2统计量和SPE两检测指标进行优化,并根据优化后的检测指标T2统计量和SPE对机械设备进行故障检测,检测得到过渡过程机械设备的故障数据:
(1)将多工况过渡过程中的上一工况的权重w1分和下一工况的权重w2分别取为:
式中,t1表示过渡过程开始的时刻,t2表示过渡过程结束的时刻,Δt=t2-t1;w1+w2=1;
(2)假设t=k时刻,分别表示上一工况和下一工况的测量变量均值,σ1和σ2分别表示上一工况和下一工况的测量变量标准差,则此时过渡过程的均值和标准差σt=k分别为:
σt=k=w1×σ1+w2×σ2;
(3)样本向量通过加权后的均值和标准差σt=k进行预处理,求出样本向量的协方差矩阵St=k和过渡过程的负载矩阵Pt=k,进而得到优化后的SPE指标和T2统计量:
3)对步骤2)中优化后的两检测指标T2统计量和SPE检测到的过渡过程机械设备故障数据进行故障重构,得到使重构后SPE最小化的幅值估计值
4)对过渡过程中不同工况下的同一故障重构后的幅值估计值进行一致性幅值估计:
(1)假设
式中,表示同一故障在上一工况和下一工况下的幅值估计值,Ξi_1、Ξi_2表示各自故障在残差空间的投影矩阵,假设维数为a;
(2)根据步骤(1)得到同一故障在不同工况下的幅值估计值为:
求解使最小化的d和e,得到同一故障在不同工况下的一致估计f1 *、
5)根据进行一致性幅值估计后的幅值估计值利用支持向量机预测模型对故障幅值进行趋势预测。
2.如权利要求1所述的一种复杂机械设备多工况故障预测方法,其特征在于:所述步骤3)中,所述故障重构的过程如下:
(1)假设检测得到的故障数据用x表示,正常数据用x*表示,x*,x∈[x1,x2,...,xm],m为每一个测量样本的维数,故障幅值的大小用f表示,则:
x=x*+Ξf,
式中,Ξ表示故障子空间;
(2)重构正常数据以消除故障的影响,得到正常数据的估计值为:
xi *=x-Ξifi,
式中,fi表示对故障幅值大小f的估计;
(3)根据步骤(2)重构后的正常数据的估计值得到使重构后SPE最小化的幅值估计值
式中,表示故障子空间Ξ到残差子空间的投影,C表示残差子空间。
3.如权利要求1或2所述的一种复杂机械设备多工况故障预测方法,其特征在于:所述步骤5)中,所述支持向量机预测方法如下:
(1)对输入的测量变量采用降噪平滑处理进行数据的预处理,根据FPE准则计算嵌入维数,然后分别形成训练和测试样本;
(2)对支持向量机模型进行初始化,将拉格朗日乘子αi和以及阈值b赋以随机初始值;
(3)利用训练样本建立目标函数,得到拉格朗日乘子αi和以及阈值b,最后选取最优的惩罚因子C*和核函数参数g;
(4)将步骤5)中的步骤(3)得到的参数值代入预测模型,用测试样本计算在未来某时刻的预测值;
(5)计算误差函数,当误差的绝对值小于预先设定值时,结束学习过程,否则返回步骤(1)继续学习。
4.如权利要求1或2所述的一种复杂机械设备多工况故障预测方法,其特征在于:所述步骤5)中,所述支持向量机预测模型对故障幅值进行趋势预测精度采用均方误差和平均相对预测误差来评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410103676.3A CN103824137B (zh) | 2014-03-20 | 2014-03-20 | 一种复杂机械设备多工况故障预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410103676.3A CN103824137B (zh) | 2014-03-20 | 2014-03-20 | 一种复杂机械设备多工况故障预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103824137A CN103824137A (zh) | 2014-05-28 |
CN103824137B true CN103824137B (zh) | 2016-08-17 |
Family
ID=50759188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410103676.3A Expired - Fee Related CN103824137B (zh) | 2014-03-20 | 2014-03-20 | 一种复杂机械设备多工况故障预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103824137B (zh) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101561026B1 (ko) * | 2014-11-27 | 2015-10-16 | 주식회사 베가스 | 설비의 예측 정비 시스템 및 방법 |
CN104898646A (zh) * | 2015-04-30 | 2015-09-09 | 东北大学 | 一种基于kpca进行故障分离与重构的电熔镁炉故障诊断方法 |
CN104950675B (zh) * | 2015-06-12 | 2018-01-23 | 华北电力大学 | 多工况电力***自适应控制方法及装置 |
CN105403705B (zh) * | 2015-12-05 | 2017-04-12 | 浙江大学 | 一种包含血糖分类功能故障检测模块的连续血糖监测设备 |
CN105700518B (zh) * | 2016-03-10 | 2018-06-12 | 华中科技大学 | 一种工业过程故障诊断方法 |
CN105893700B (zh) * | 2016-04-26 | 2019-05-31 | 陆新建 | 基于物理-大数据混合模型的化工生产在线故障检测与诊断技术 |
CN106200624B (zh) * | 2016-08-26 | 2019-01-04 | 大连海事大学 | 基于交叉分段pca的工业锅炉故障诊断方法 |
KR101863781B1 (ko) * | 2016-09-08 | 2018-06-01 | 두산중공업 주식회사 | 로터 진동 이상 감지 장치 및 방법 |
CN106952028A (zh) * | 2017-03-13 | 2017-07-14 | 杭州安脉盛智能技术有限公司 | 机电装备故障预诊与健康管理方法及*** |
JP6800074B2 (ja) * | 2017-03-31 | 2020-12-16 | アズビル株式会社 | 劣化診断装置 |
CN108416386A (zh) * | 2018-03-08 | 2018-08-17 | 云南电网有限责任公司电力科学研究院 | 一种用于水电机组轴承温度异常判断的方法及*** |
CN108490923B (zh) * | 2018-04-28 | 2020-09-15 | 南京航空航天大学 | 用于电气牵引***的微小故障检测与定位的***设计方法 |
JP7143148B2 (ja) * | 2018-08-23 | 2022-09-28 | 三菱重工業株式会社 | 予測装置、予測方法、及びプログラム |
CN110009126B (zh) * | 2019-01-23 | 2021-10-01 | 北京化工大学 | 基于pls模型与pca贡献度融合的在线报警分析方法 |
CN110187206B (zh) * | 2019-05-22 | 2021-11-02 | 中国人民解放军国防科技大学 | 一种复杂工况下非高斯过程的悬浮***的故障检测方法 |
CN110045617B (zh) * | 2019-05-22 | 2021-12-07 | 杭州电子科技大学 | 一种工业过程约束预测先进控制方法 |
CN110543153B (zh) * | 2019-08-20 | 2020-11-17 | 华中科技大学 | 一种多工况工业过程的过渡工况故障检测方法 |
CN110794814B (zh) * | 2019-11-27 | 2022-06-28 | 中国人民解放***箭军工程大学 | 一种基于广义主成分的故障确定方法及*** |
CN111061252A (zh) * | 2019-12-24 | 2020-04-24 | 浙江大学 | 一种基于数据驱动的精馏塔智能看板 |
CN112836431A (zh) * | 2021-02-05 | 2021-05-25 | 河北工业大学 | 基于pso-lstm的青霉素发酵过程故障预测方法 |
CN115144182B (zh) * | 2022-09-01 | 2023-01-17 | 杭州景业智能科技股份有限公司 | 轴承健康状态监测方法、装置、计算机设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201295A (zh) * | 2006-12-13 | 2008-06-18 | 上海海事大学 | 旋转机械小波灰色故障预测方法及装置 |
CN101339434A (zh) * | 2008-08-22 | 2009-01-07 | 唐山学院 | 机械设备故障在线诊断装置及方法 |
CN101799369A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障监测预报中的灰色建模方法 |
CN101799359A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障监测预报方法及其*** |
CN101798957A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障诊断方法 |
CN101799366A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种机械故障预测的特征提取方法 |
-
2014
- 2014-03-20 CN CN201410103676.3A patent/CN103824137B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201295A (zh) * | 2006-12-13 | 2008-06-18 | 上海海事大学 | 旋转机械小波灰色故障预测方法及装置 |
CN101339434A (zh) * | 2008-08-22 | 2009-01-07 | 唐山学院 | 机械设备故障在线诊断装置及方法 |
CN101799369A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障监测预报中的灰色建模方法 |
CN101799359A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障监测预报方法及其*** |
CN101798957A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种动力设备故障诊断方法 |
CN101799366A (zh) * | 2010-01-27 | 2010-08-11 | 北京信息科技大学 | 一种机械故障预测的特征提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103824137A (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103824137B (zh) | 一种复杂机械设备多工况故障预测方法 | |
CN106951695B (zh) | 多工况下的机械设备剩余使用寿命计算方法及*** | |
CN108062565B (zh) | 基于化工te过程的双主元-动态核主元分析故障诊断方法 | |
Harmouche et al. | Incipient fault detection and diagnosis based on Kullback–Leibler divergence using principal component analysis: Part I | |
EP3422518B1 (en) | A method for recognizing contingencies in a power supply network | |
CN105259895B (zh) | 一种工业过程微小故障的检测和分离方法及其监测*** | |
CN104914847B (zh) | 基于方向核偏最小二乘的工业过程故障诊断方法 | |
CN104714537B (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
Jiang et al. | Weighted kernel principal component analysis based on probability density estimation and moving window and its application in nonlinear chemical process monitoring | |
EP3422517B1 (en) | A method for recognizing contingencies in a power supply network | |
CN109522948A (zh) | 一种基于正交局部保持投影的故障检测方法 | |
CN108490923A (zh) | 用于电气牵引***的微小故障检测与定位的***设计方法 | |
CN103197663B (zh) | 一种故障预测方法及*** | |
CN111680725A (zh) | 基于重构贡献的气体传感器阵列多故障隔离算法 | |
Shi et al. | Online event detection in synchrophasor data with graph signal processing | |
CN109298633A (zh) | 基于自适应分块非负矩阵分解的化工生产过程故障监测方法 | |
Wang et al. | Fault detection and identification using a Kullback-Leibler divergence based multi-block principal component analysis and Bayesian inference | |
CN109799808A (zh) | 一种基于重构技术的动态过程故障预测方法 | |
Bo et al. | An integrated method of independent component analysis and support vector machines for industry distillation process monitoring | |
CN103631145A (zh) | 基于监控指标切换的多工况过程监控方法和*** | |
CN105137354A (zh) | 一种基于神经网络电机故障检测方法 | |
CN115081673A (zh) | 一种油气管道的异常预测方法、装置、电子设备及介质 | |
Zhang et al. | A local and global statistics pattern analysis method and its application to process fault identification | |
CN109213120A (zh) | 基于高维特征表示下多级主元空间的故障诊断方法 | |
CN111474476A (zh) | 一种电机故障预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |