CN114063456B - 利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 - Google Patents
利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 Download PDFInfo
- Publication number
- CN114063456B CN114063456B CN202111350431.7A CN202111350431A CN114063456B CN 114063456 B CN114063456 B CN 114063456B CN 202111350431 A CN202111350431 A CN 202111350431A CN 114063456 B CN114063456 B CN 114063456B
- Authority
- CN
- China
- Prior art keywords
- fault
- vector
- calculating
- moment
- state
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了一种利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,所述方法包括如下步骤:一、建立故障自回归退化模型的状态方程与测量方程;二、计算状态向量预测值三、计算***状态协方差矩阵的预测值P(k|k‑1);四、计算滤波增益系数K(k);五、计算***状态协方差矩阵的估计值P(k);六、计算状态向量估计值七、计算测量向量的预测值到八、计算测量向量预测值的上边界到和下边界y(k)到y(k+m);九、确定故障预警结果;十、每个k时刻重复二~九,迭代实现故障预测和预警。该方法能够在控制***发生故障之前给出可靠的预警信号,有效保障***的安全性及稳健性,其原理清晰,算法简单,易于实际工程实现。
Description
技术领域
本发明涉及一种控制***的故障预警方法,具体涉及一种利用自回归模型和卡尔曼滤波算法进行控制***故障预测和预警的方法。
背景技术
现代控制***变得越来越复杂,相应地,控制***的安全性和可靠性就显得尤为重要。但是,实际***运行时发生的各种故障会降低其可靠性甚至破坏***的稳定性从而造成严重的安全事故。因此,为了提高控制***的安全性和可靠性,必须及时地诊断出故障并采取应对措施,让故障对***造成的损害最小。通常来说,对控制***安全性影响更多的是存在其内部的各种缓变故障,这些缓变故障是由***部件的性能退化所引起的。缓变故障在初期,幅值较小,不会对控制***的性能造成严重影响,但随着时间维度的不断积累,其故障幅值会逐渐增大,对***安全性和可靠性的威胁也就越高。特别地,当缓变故障演变至***部件性能失效阈值时,则会严重威胁整个控制***的安全性与稳定性。因此,针对控制***缓变故障进行退化过程建模,预测其未来的故障发展趋势并在性能失效前及时做出故障预警,对于提升控制***的可靠性和安全性具有重要意义。
根据国内外已发表的相关文献可知,目前的故障预测和预警算法大都基于模糊网络或深度学习技术来进行故障退化建模和预测,对故障的数据量和模型训练的计算量要求较高,适用范围有限。同时,现有方法大都需要线下预先训练故障模型,无法在线上实时更新和调整模型参数,难以提供准确的故障预测结果,这会在一定程度上影响故障预警的可靠性。
由于使用模糊网络或深度学习技术来进行故障退化建模和预测需要高额的计算量,因此这类算法仅适用于组成结构简单、部件数量较少的控制***。面对现代控制***结构日益复杂,组成部件的数目不断增加的发展趋势,这些故障预测和预警算法已不能满足高效准确故障预警性能的需求,现代控制***迫切需要适用性好,计算量小且性能可靠的故障预警算法。
发明内容
本发明的目的是提供一种利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,该方法能够在控制***发生故障之前给出可靠的预警信号,有效保障***的安全性及稳健性,其原理清晰,算法简单,易于实际工程实现。
本发明的目的是通过以下技术方案实现的:
一种利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,包括如下步骤:
步骤一、建立故障自回归退化模型的状态方程与测量方程:
式中,x(k)、x(k-1)分别为k、k-1时刻***的状态向量;A(k-1)为退化模型状态的一步转移阵;w(k-1)为***的扰动向量;C为测量矩阵;v(k)为***的测量噪声向量;y(k)为测量向量;
步骤三、计算***状态协方差矩阵的预测值P(k|k-1):
P(k|k-1)=A(k-1)P(k-1)A(k-1)T+Q;
式中,P(k-1)为k-1时刻***误差协方差阵的估计值;Q为***噪声协方差阵;
步骤四、计算滤波增益系数K(k):
K(k)=P(k|k-1)CT/(CP(k|k-1)CT+R);
式中,R为测量噪声向量的协方差矩阵;
步骤五、计算***状态协方差矩阵的估计值P(k):
P(k)=P(k|k-1)-K(k)CP(k|k-1);
式中,P(k)为k时刻状态协方差矩阵的估计值;
步骤九、确定故障预警结果:
式中,yth为设定的故障检测阈值;
若对j从0取到m,所有的j都使得上式成立,则说明在k时刻未预测到故障;
若对j从0取到m,存在一个j使得上式不成立,则说明k时刻预测到了故障,预测的未来发生故障的时刻为k+j*;
步骤十、每个k时刻重复步骤二~步骤九,迭代实现故障预测和预警。
本发明提供了一种利用自回归模型和卡尔曼滤波算法的故障预警方法,该方法能依据故障的退化特性来预测故障未来的变化值,并在此基础上进行有效的故障早期预警,其优点和有益效果总结如下:
(1)采用自回归模型描述故障的退化过程,模型形式简单且描述范围广泛,故障预警方法适用性好;
(2)利用卡尔曼滤波算法实时在线更新模型参数,参数估计精度高,线上运行计算量小,易于实际硬件平台实现;
(3)故障预警信号通过故障预测值及其上下界给出,保证了预警结果的可靠性,进而提高了控制***的安全性。
附图说明
图1为本发明利用自回归模型和卡尔曼滤波算法的故障预测与预警方法的流程图。
图2为直流电机缓变故障数据结果。
图3为k=180时电机缓变故障预测与预警结果。
图4为k=185时电机缓变故障预测与预警结果。
图5为k=190时电机缓变故障预测与预警结果。
图6为k=195时电机缓变故障预测与预警结果。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围中。
本发明提供了一种利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,所述方法包括如下步骤:
步骤一、建立故障自回归退化模型的状态方程与测量方程,其公式如下:
式中,x(k)、x(k-1)分别为k、k-1时刻***的状态向量;A(k-1)为退化模型状态的一步转移阵;w(k-1)为***的扰动向量;C为测量矩阵;v(k)为***的测量噪声向量;y(k)为测量向量。
步骤三、计算***状态协方差矩阵的预测值P(k|k-1),其公式如下:
P(k|k-1)=A(k-1)P(k-1)A(k-1)T+Q (3);
式中,P(k-1)为k-1时刻***误差协方差阵的估计值;Q为***噪声协方差阵。
步骤四、计算滤波增益系数K(k),其公式如下:
K(k)=P(k|k-1)CT/(CP(k|k-1)CT+R) (4);
式中,R为测量噪声向量的协方差矩阵。
步骤五、计算***状态协方差矩阵的估计值P(k),其公式如下:
P(k)=P(k|k-1)-K(k)CP(k|k-1) (5);
式中,P(k)为k时刻状态协方差矩阵的估计值。
步骤九、确定故障预警结果,其中,故障预警策略如下:
若对j从0取到m,所有的j都使得下式成立:
则说明在k时刻未预测到故障。
若对j从0取到m,存在一个j使得式(9)不成立,则说明k时刻预测到了故障,预测的未来发生故障的时刻为k+j*,其中,yth为设定的故障检测阈值。
步骤十、每个k时刻重复步骤二~步骤九,迭代实现故障预测和预警。
优选地,在步骤一中,***状态向量由公式(10)给出:
x=[f c1 c2 …cn]T (10);
其中,f为工程***的故障值,由测量得到,c1、c1、…、cn为自回归模型的系数;测量值y=f;矩阵A(k-1)和C分别由公式(11)和公式(12)给出:
C=[1 0 0 … 0] (12)。
优选地,在步骤二中,从k=n+1开始迭代计算,前n+1个时刻的状态向量估计值分别由公式(13)和公式(14)给出:
优选地,在步骤九中,故障预警结果通过测量向量预测值及其上下边界与故障检测阈值的大小关系来共同确定。
综上所述,本发明首先利用自回归模型来描述故障的动态退化过程;再利用卡尔曼滤波算法来在线实时更新模型参数,进而实现对未来时刻故障值及其上下界的预测;最后,结合故障预测结果和设置的故障阈值即可实现对故障的有效预警。本发明利用自回归模型描述故障退化动态,模型形式简单且适用好,拓宽了算法的应用范围;该方法利用卡尔曼滤波实时在线更新模型参数,参数估计精度高,算法简单,原理清晰,便于实际硬件实现;该算法利用故障预测值及其上下界给出故障预警信号,保证了预警结果的可靠性。
应用实例
下面以一个直流电机的实例来说明利用自回归模型和卡尔曼滤波算法的故障预测与预警方法的应用过程。考虑直流电机随运行时间的增加,其性能会逐渐退化从而导致正常运行条件下的电机温度会不断增高,因此可以将正常运行条件下的电机温度和标称温度的差值作为直流电机缓变故障的指示信号,通过采集的正常运行条件下的电机温度数据则可绘制如图2所示的的缓变故障数据结果。
故障预测算法仿真中,设置故障从当前时刻往后进行预测的步长为m=20,故障自回归退化模型的阶数为n=20,退化模型的输出矩阵C=[1 0 0 … 0]。故障检测阈值选取为yth=3.8,卡尔曼滤波算法的初值选为/>状态初值/>的协方差矩阵选为P(k)=10-4In+1,In+1为n+1维的单位矩阵;扰动w(k)和噪声v(k)的协方差矩阵为Q=10-4In+1和R=10-4。基于给定的仿真参数和测量的缓变故障数据可以得到如图3-图6所示的故障预测仿真结果。
从图3-图6中的仿真结果可以看出,在k=180和k=185时,直流电机的实际缓变故障数据未超出给定的故障检测阈值,此时所设计的方法同样没有给出故障预警信号。在k=190和k=195时,直流电机的实际缓变故障数据超出了给定的故障检测阈值,此时所提出的方法均可以快速准确的预测故障并给出预警信号,验证了所提出方法的有效性。
Claims (4)
1.一种利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,其特征在于所述方法包括如下步骤:
步骤一、建立故障自回归退化模型的状态方程与测量方程:
式中,x(k)、x(k-1)分别为k、k-1时刻***的状态向量;A(k-1)为退化模型状态的一步转移阵;w(k-1)为***的扰动向量;C为测量矩阵;v(k)为***的测量噪声向量;y(k)为测量向量;
***的状态向量由如下公式给出:
x=[f c1 c2…cn]T;
其中,f为工程***的故障值,c1、c1、…、cn为自回归模型的系数;
矩阵A(k-1)由如下公式给出:
其中,f为工程***的故障值;
步骤三、计算***状态协方差矩阵的预测值P(k|k-1):
P(k|k-1)=A(k-1)P(k-1)A(k-1)T+Q;
式中,P(k-1)为k-1时刻***误差协方差阵的估计值;Q为***噪声协方差阵;
步骤四、计算滤波增益系数K(k):
K(k)=P(k|k-1)CT/(CP(k|k-1)CT+R);
式中,R为测量噪声向量的协方差矩阵;
步骤五、计算***状态协方差矩阵的估计值P(k):
P(k)=P(k|k-1)-K(k)CP(k|k-1);
式中,P(k)为k时刻状态协方差矩阵的估计值;
式中,到/>为k到k+m时刻测量向量预测值的上边界;y(k)到y(k+m)为k到k+m时刻测量向量预测值的下边界;测量向量预测值的上边界/>到/>和下边界y(k)到y(k+m)通过卡尔曼滤波的3σ区间确定,/>为/>的标准差;
步骤九、确定故障预警结果:
式中,yth为设定的故障检测阈值;
若对j从0取到m,所有的j都使得上式成立,则说明在k时刻未预测到故障;
若对j从0取到m,存在一个j使得上式不成立,则说明k时刻预测到了故障,预测的未来发生故障的时刻为k+j*;
步骤十、每个k时刻重复步骤二~步骤九,迭代实现故障预测和预警。
2.根据权利要求1所述的利用自回归模型和卡尔曼滤波算法的故障预测与预警方法,其特征在于所述步骤一中,C分别由如下公式给出:
C=[1 0 0…0]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111350431.7A CN114063456B (zh) | 2021-11-15 | 2021-11-15 | 利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111350431.7A CN114063456B (zh) | 2021-11-15 | 2021-11-15 | 利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114063456A CN114063456A (zh) | 2022-02-18 |
CN114063456B true CN114063456B (zh) | 2023-06-02 |
Family
ID=80272295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111350431.7A Active CN114063456B (zh) | 2021-11-15 | 2021-11-15 | 利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114063456B (zh) |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4812995A (en) * | 1987-05-19 | 1989-03-14 | Girgis Adly A | Adaptive Kalman Filtering in fault classification |
US5210704A (en) * | 1990-10-02 | 1993-05-11 | Technology International Incorporated | System for prognosis and diagnostics of failure and wearout monitoring and for prediction of life expectancy of helicopter gearboxes and other rotating equipment |
TR199600527A2 (xx) * | 1996-06-24 | 1998-01-21 | Ar�El�K A.�. | Elektrik motorlar� i�in model bazl� hata tespit ve te�his sistemi. |
US20050114023A1 (en) * | 2003-11-26 | 2005-05-26 | Williamson Walton R. | Fault-tolerant system, apparatus and method |
CN101859146B (zh) * | 2010-07-16 | 2011-11-30 | 哈尔滨工业大学 | 一种基于预测滤波和经验模态分解的卫星故障预测方法 |
CN108845495B (zh) * | 2018-04-03 | 2021-08-03 | 南通大学 | 基于双层Kalman滤波器的间歇故障诊断与主动容错控制方法 |
CN110647721A (zh) * | 2018-06-27 | 2020-01-03 | 北京自动化控制设备研究所 | 一种多元信息融合故障检测优化方法 |
CN109508475A (zh) * | 2018-10-17 | 2019-03-22 | 北京航空航天大学 | 一种基于多维卡尔曼滤波修正的故障主动预测方法 |
CN112949026B (zh) * | 2021-01-19 | 2023-05-23 | 中国人民解放***箭军工程大学 | 一种考虑年龄和状态依赖的退化设备剩余寿命预测方法 |
-
2021
- 2021-11-15 CN CN202111350431.7A patent/CN114063456B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114063456A (zh) | 2022-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110298455B (zh) | 一种基于多变量估计预测的机械设备故障智能预警方法 | |
CN107145645B (zh) | 带不确定冲击的非平稳退化过程剩余寿命预测方法 | |
Nozari et al. | Model-based robust fault detection and isolation of an industrial gas turbine prototype using soft computing techniques | |
WO2019080367A1 (zh) | 一种机械设备健康状态评估方法 | |
CN112305418B (zh) | 一种基于混合噪声双重滤波的电机***故障诊断方法 | |
CN102829967A (zh) | 一种基于回归模型系数变化的时域故障识别方法 | |
CN114239377A (zh) | 一种城轨机电设备健康状态评估方法、***及存储介质 | |
CN103197663A (zh) | 一种故障预测方法及*** | |
CN116165885B (zh) | 一种高速列车的无模型自适应鲁棒控制方法及*** | |
CN110187206B (zh) | 一种复杂工况下非高斯过程的悬浮***的故障检测方法 | |
CN113791351B (zh) | 基于迁移学习和差值概率分布的锂电池寿命预测方法 | |
CN109324591B (zh) | 基于鲁棒估计器的随机混杂***的故障诊断装置及方法 | |
CN114063456B (zh) | 利用自回归模型和卡尔曼滤波算法的故障预测与预警方法 | |
CN112016208B (zh) | 一种考虑扰动的隐含故障诊断方法及*** | |
Le Pocher et al. | Sensor fault detection of a real undershot/overshot gate based on physical and nonlinear black-box models | |
CN112488295A (zh) | 交叉验证算法优化lstm网络继电器贮存寿命预测方法 | |
CN111563626A (zh) | 一种电力***预测辅助状态估计方法及*** | |
CN108282360B (zh) | 一种长短期预测融合的故障检测方法 | |
CN115407662A (zh) | 一种基于扩张故障检测观测器和集员估计的鲁棒故障检测方法 | |
CN114135580B (zh) | 磁轴承转子的位置评估方法及装置 | |
CN113221248B (zh) | 一种基于pf-garch模型的船舶***设备状态参数预测方法 | |
CN113236506B (zh) | 一种基于滤波的工业时滞***故障检测方法 | |
CN114254555A (zh) | 基于深度学习的反作用轮故障检测和健康评估***及方法 | |
CN114676384B (zh) | 一种基于粒子滤波的性能状态百分位值估计方法 | |
CN109978349B (zh) | 一种武器装备指数型和联***的可靠度置信下限计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |