留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于动态贝叶斯网络的电源系统可靠性分析与故障诊断

李享 黄洪钟 黄鹏 李彦锋

李享, 黄洪钟, 黄鹏, 李彦锋. 基于动态贝叶斯网络的电源系统可靠性分析与故障诊断[J]. 电子科技大学学报, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
引用本文: 李享, 黄洪钟, 黄鹏, 李彦锋. 基于动态贝叶斯网络的电源系统可靠性分析与故障诊断[J]. 电子科技大学学报, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
LI Xiang, HUANG Hong-zhong, HUANG Peng, LI Yan-feng. Reliability Analysis and Fault Diagnosis for Power System via Dynamic Bayesian Network[J]. Journal of University of Electronic Science and Technology of China, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
Citation: LI Xiang, HUANG Hong-zhong, HUANG Peng, LI Yan-feng. Reliability Analysis and Fault Diagnosis for Power System via Dynamic Bayesian Network[J]. Journal of University of Electronic Science and Technology of China, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416

基于动态贝叶斯网络的电源系统可靠性分析与故障诊断

doi: 10.12178/1001-0548.2020416
基金项目: 国家自然科学基金(51875089)
详细信息
    作者简介:

    李享(1989-),男,博士生,主要从事基于动态贝叶斯网络的系统可靠性分析方面的研究

    通讯作者: 黄洪钟,E-mail:hzhuang@uestc.edu.cn
  • 中图分类号: TM417

Reliability Analysis and Fault Diagnosis for Power System via Dynamic Bayesian Network

  • 摘要: 动态系统的可靠性分析与故障诊断一直是可靠性领域的热点及难点问题,作为该领域热门的分析工具之一,动态贝叶斯网络(DBN)得到了充分的应用与开发。但是,现有的DBN算法受限于系统的失效分布类型,且建模难度也随着系统复杂度的增加而呈指数增长。针对以上问题,该文提出一种改进的动态贝叶斯网络概率表建模方法,在连续任务时间的条件下,实现动态系统的可靠性分析。然后,结合DBN双向推理算法,求解系统失效时部件失效的后验概率,并将计算结果应用于系统故障诊断及薄弱部件定位。最后,结合某电源系统的可靠性分析与故障诊断,验证了该方法的实用性。
  • 图  1  CTBN中“与”门的CPT&MPT

    图  2  “与”门的CPT&MPT

    图  3  温备份门的CPT&MPT

    图  4  结果对比分析

    图  5  电源系统结构框图

    图  6  电源系统的DBN

    图  7  电源系统可靠度曲线

    图  8  电源系统故障诊断

    表  1  温备份节点B的CPT

    A 12
    B1$1 - \dfrac{{{{\rm{e}}^{ - \lambda t}}\left( {1 - {{\rm{e}}^{ - \alpha \lambda t}}} \right)}}{{\alpha \left( {1 - {{\rm{e}}^{ - \lambda t}}} \right)}}$$1 - {{\rm{e}}^{ - \alpha \lambda t}}$
    2$\dfrac{{{{\rm{e}}^{ - \lambda t}}\left( {1 - {{\rm{e}}^{ - \alpha \lambda t}}} \right)}}{{\alpha \left( {1 - {{\rm{e}}^{ - \lambda t}}} \right)}}$${{\rm{e}}^{\left( { - \alpha \lambda t} \right)}}$
    下载: 导出CSV

    表  2  事件节点编号及失效率

    节点
    编号
    名称节点
    编号
    名称失效率(${\lambda _i}$)/(10−6/h)
    TM电源系统失效AM配电单元2.67900
    BM1充电失效EM调波单元0.01120
    BM2充电模块失效B1低压供电模块0.29400
    CM1储能失效B2主充电模块2.16040
    CM2储能电容失效B3预充电模块1.08020
    DM开关失效B4控制采集模块7.62380
    C1主储能电容模块0.00580
    C2预储能电容模块0.00058
    C3泄放单元模块0.65470
    C4充电保护单元模块0.00780
    D1主开关单元模块7.53300
    D2预开关单元模块5.27310
    下载: 导出CSV
  • [1] 肖建平. 高压开关电源的拓扑研究[J]. 电子科技大学学报, 2007, 36(4): 726-729. doi:  10.3969/j.issn.1001-0548.2007.04.021

    XIAO Jian-ping. Research on the topology of high-voltage switching power supply[J]. Journal of University of Electronic Science and Technology of China, 2007, 36(4): 726-729. doi:  10.3969/j.issn.1001-0548.2007.04.021
    [2] 王京梅, 兰中文, 余忠, 等. 高频开关电源变压器的优化设计[J]. 电子科技大学学报, 2002, 31(4): 362-365. doi:  10.3969/j.issn.1001-0548.2002.04.007

    WANG Jing-mei, LAN Zhong-wen, YU Zhong, et al. Optimum design of transformer in high frequency switch mode power[J]. Journal of University of Electronic Science and Technology of China, 2002, 31(4): 362-365. doi:  10.3969/j.issn.1001-0548.2002.04.007
    [3] HUANG C G, HUANG H Z, LI Y F. A bidirectional LSTM prognostics method under multiple operational conditions[J]. IEEE Transactions on Industrial Electronics, 2019, 66(11): 8792-8802. doi:  10.1109/TIE.2019.2891463
    [4] PENG W, LI Y F, YANG Y J, et al. Bayesian degradation analysis with inverse Gaussian process models under time-varying degradation rates[J]. IEEE Transactions on Reliability, 2017, 66(1): 84-96. doi:  10.1109/TR.2016.2635149
    [5] JIANG T, LIU Y. Parameter inference for non-repairable multi-state system reliability models by multi-level observation sequences[J]. Reliability Engineering & System Safety, 2017, 166: 3-15.
    [6] MI J, LI Y F, YANG Y J, et al. Reliability assessment of complex electromechanical systems under epistemic uncertainty[J]. Reliability Engineering & System Safety, 2016, 152: 1-15.
    [7] ZHANG Y, ZHOU G, JIN J, et al. Sparse Bayesian multiway canonical correlation analysis for EEG pattern recognition[J]. Neurocomputing, 2017, 225: 103-110. doi:  10.1016/j.neucom.2016.11.008
    [8] PINTO H, GATES I, WANG X. Bayesian biclustering by dynamics: A clustering algorithm for SAGD time series data[J]. Computers & Geosciences, 2019, 133: 104304.
    [9] 米金华, 李彦锋, 彭卫文, 等. 复杂多态系统的区间值模糊贝叶斯网络建模与分析[J]. 中国科学: 物理学 力学 天文学, 2018, 48(1): 54-66.

    MI Jin-hua, LI Yan-feng, PENG Wei-wen, et al. Reliability modeling and analysis of complex multi-state system based on interval fuzzy Bayesian network[J]. Scientia Sinica Physica, Mechanica & Astronomica, 2018, 48(1): 54-66.
    [10] LI H, SOARES C G, HUANG H Z. Reliability analysis of a floating offshore wind turbine using Bayesian Networks[J]. Ocean Engineering, 2020, 217: 107827. doi:  10.1016/j.oceaneng.2020.107827
    [11] REBELLO S, YU H, MA L. An integrated approach for system functional reliability assessment using dynamic Bayesian network and hidden markov model[J]. Reliability Engineering & System Safety, 2018, 180: 124-135.
    [12] 李彦锋, 黄洪钟, 刘宇, 等. 基于贝叶斯网络的卫星太阳翼驱动机构系统可靠性建模与评估[J]. 中国科技论文, 2012, 7(8): 583-588.

    LI Yan-feng, HUANG Hong-zhong, LIU Yu, et al. System reliability modeling and assessment of satellite solar array drive assembly based on a Bayesian network[J]. China Sciencepaper, 2012, 7(8): 583-588.
    [13] 王晓明, 李彦锋, 李爱峰, 等. 模糊数据下基于连续时间贝叶斯网络的整流回馈系统可靠性建模与评估[J]. 机械工程学报, 2015, 51(14): 167-174. doi:  10.3901/JME.2015.14.167

    WANG Xiao-ming, LI Yan-feng, LI Ai-feng, et al. Reliability modeling and evaluation for rectifier feedback system based on continuous time Bayesian networks under fuzzy numbers[J]. Journal of Mechanical Engineering, 2015, 51(14): 167-174. doi:  10.3901/JME.2015.14.167
    [14] FAN D, WANG Z, LIU L, et al. A modified GO-FLOW methodology with common cause failure based on discrete time Bayesian network[J]. Nuclear Engineering and Design, 2016, 305: 476-488. doi:  10.1016/j.nucengdes.2016.06.010
    [15] SONG Y, MI J, CHENG Y, et al. Application of discrete-time Bayesian network on reliability analysis of uncertain system with common cause failure[J]. Quality and Reliability Engineering International, 2019, 35(4): 1025-1045. doi:  10.1002/qre.2443
    [16] LI X, LI YF, LI H, et al. An algorithm of discrete-time Bayesian network for reliability analysis of multilevel system with warm spare gate[J]. Quality and Reliability Engineering International, 2020, 37(3): 1-19.
    [17] STURLAUGSON L, PERREAULT L, SHEPPARD J W. Factored performance functions and decision making in continuous time Bayesian networks[J]. Journal of Applied Logic, 2017, 22: 28-45. doi:  10.1016/j.jal.2016.11.030
    [18] LI Y F, MI J, LIU Y U, et al. Dynamic fault tree analysis based on continuous-time Bayesian networks under fuzzy numbers[J]. Journal of Risk and Reliability, 2015, 229(6): 530-541. doi:  10.1177/1748006X15588446
    [19] BOUDALI H, DUGAN J B. A continuous-time Bayesian network reliability modeling, and analysis framework[J]. IEEE Transactions on Reliability, 2006, 55(1): 86-97. doi:  10.1109/TR.2005.859228
  • [1] 巩彬, 安爱民, 石耀科, 杜先君.  一种基于IMODA自适应深度信念网络的复杂模拟电路故障诊断方法 . 电子科技大学学报, 2024, 53(): 1-18. doi: 10.12178/1001-0548.2023047
    [2] 黄洪钟, 童辉, 邓智铭, 钱华明, 李彦锋, 郭超.  基于ALK-MCS算法的航空发动机轴承可靠性分析 . 电子科技大学学报, 2022, 51(3): 438-443. doi: 10.12178/1001-0548.2022021
    [3] 杜娟, 刘志刚, 宋考平, 杨二龙.  基于卷积神经网络的抽油机故障诊断 . 电子科技大学学报, 2020, 49(5): 751-757. doi: 10.12178/1001-0548.2019205
    [4] 肖宁聪, 袁凯, 王永山.  基于序列代理模型的结构可靠性分析方法 . 电子科技大学学报, 2019, 48(1): 156-160. doi: 10.3969/j.issn.1001-0548.2019.01.023
    [5] 邵继业, 谢昭灵, 杨瑞.  基于GA-PSO优化BP神经网络的压缩机气阀故障诊断 . 电子科技大学学报, 2018, 47(5): 781-787. doi: 10.3969/j.issn.1001-0548.2018.05.023
    [6] 周启忠, 谢永乐, 毕东杰, 李西峰.  种模拟电路参数型故障诊断新方法 . 电子科技大学学报, 2017, 46(2): 386-391. doi: 10.3969/j.issn.1001-0548.2017.02.012
    [7] 姜书艳, 罗刚, 夏登明, 李琦, 宋国明.  片上网络互连线延迟故障测试方法研究 . 电子科技大学学报, 2016, 45(4): 557-563. doi: 10.3969/j.issn.1001-0548.2016.04.008
    [8] 闫涛, 赵文俊, 胡秀洁, 宋家友.  基于信息融合技术的航空电子设备故障诊断研究 . 电子科技大学学报, 2015, 44(3): 392-395. doi: 10.3969/j.issn.1001-0548.2015.03.013
    [9] 李从善, 刘天琪, 李兴源, 曹喜民, 刘利兵.  基于排列熵算法的电力系统故障信号分析 . 电子科技大学学报, 2015, 44(2): 233-238. doi: 10.3969/j.issn.1001-0548.2015.02.013
    [10] 李旻, 咸卫明, 龙兵, 王厚军.  基于特征优选模拟电路故障诊断方法 . 电子科技大学学报, 2014, 43(4): 557-561. doi: 10.3969/j.issn.1001-0548.2014.04.015
    [11] 张旭东, 黄洪钟, XIE Min.  单环法求解基于可靠性的设计优化问题 . 电子科技大学学报, 2011, 40(1): 157-160. doi: 10.3969/j.issn.1001-0548.2011.01.030
    [12] 杜丽, 肖宁聪, 黄洪钟, 何俐萍, 李海庆.  认知不确定性的谐波齿轮减速器可靠性分析研究 . 电子科技大学学报, 2011, 40(3): 470-475. doi: 10.3969/j.issn.1001-0548.2011.03.028
    [13] 张健, 徐红兵, 王情.  多重分形谱在集成电路动态电流故障诊断中的应用 . 电子科技大学学报, 2011, 40(4): 544-548.
    [14] 孙永奎, 陈光礻禹, 李辉.  灵敏度分析和SVM诊断模拟电路故障的方法 . 电子科技大学学报, 2009, 38(6): 971-974. doi: 10.3969/j.issn.1001-0548.2009.06.016
    [15] 雷霖, 代传龙, 王厚军, 赵旭.  粗糙集-神经网络集成的WSN节点故障诊断 . 电子科技大学学报, 2008, 37(4): 565-568.
    [16] 杨俊华, 尚志恩, 吕锋.  基于布尔差分的数字逻辑电路故障诊断 . 电子科技大学学报, 2005, 34(4): 517-520.
    [17] 冯兰君, 李建, 吕炳朝.  燃气轮机故障诊断专家系统的设计与实现 . 电子科技大学学报, 2000, 29(3): 278-281.
    [18] 骆德汉, 郎文辉, 黄胜.  基于人工神经网络的在线设备状态监测系统的研究 . 电子科技大学学报, 1999, 28(5): 527-532.
    [19] 王厚军, 陈相富, 陈光.  基于总线观测的微机故障诊断专家系统研究 . 电子科技大学学报, 1998, 27(2): 161-165.
    [20] 潘中良, 陈光.  模拟电路故障诊断的神经网络专家系统方法 . 电子科技大学学报, 1997, 26(4): 405-408.
  • 加载中
图(8) / 表(2)
计量
  • 文章访问数:  5128
  • HTML全文浏览量:  1524
  • PDF下载量:  50
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-11-23
  • 修回日期:  2021-03-03
  • 网络出版日期:  2021-07-23
  • 刊出日期:  2021-06-28

基于动态贝叶斯网络的电源系统可靠性分析与故障诊断

doi: 10.12178/1001-0548.2020416
    基金项目:  国家自然科学基金(51875089)
    作者简介:

    李享(1989-),男,博士生,主要从事基于动态贝叶斯网络的系统可靠性分析方面的研究

    通讯作者: 黄洪钟,E-mail:hzhuang@uestc.edu.cn
  • 中图分类号: TM417

摘要: 动态系统的可靠性分析与故障诊断一直是可靠性领域的热点及难点问题,作为该领域热门的分析工具之一,动态贝叶斯网络(DBN)得到了充分的应用与开发。但是,现有的DBN算法受限于系统的失效分布类型,且建模难度也随着系统复杂度的增加而呈指数增长。针对以上问题,该文提出一种改进的动态贝叶斯网络概率表建模方法,在连续任务时间的条件下,实现动态系统的可靠性分析。然后,结合DBN双向推理算法,求解系统失效时部件失效的后验概率,并将计算结果应用于系统故障诊断及薄弱部件定位。最后,结合某电源系统的可靠性分析与故障诊断,验证了该方法的实用性。

English Abstract

李享, 黄洪钟, 黄鹏, 李彦锋. 基于动态贝叶斯网络的电源系统可靠性分析与故障诊断[J]. 电子科技大学学报, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
引用本文: 李享, 黄洪钟, 黄鹏, 李彦锋. 基于动态贝叶斯网络的电源系统可靠性分析与故障诊断[J]. 电子科技大学学报, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
LI Xiang, HUANG Hong-zhong, HUANG Peng, LI Yan-feng. Reliability Analysis and Fault Diagnosis for Power System via Dynamic Bayesian Network[J]. Journal of University of Electronic Science and Technology of China, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
Citation: LI Xiang, HUANG Hong-zhong, HUANG Peng, LI Yan-feng. Reliability Analysis and Fault Diagnosis for Power System via Dynamic Bayesian Network[J]. Journal of University of Electronic Science and Technology of China, 2021, 50(4): 603-608. doi: 10.12178/1001-0548.2020416
  • 电源系统是由配电设备、储能设备、调波设备、开关设备、充电控制设备以及相关电线电路组成的总体[1],为各种电机提供高/低频、交/直流电源,其可靠性直接影响整个电机系统的平稳运行[2]。然而,电源系统结构复杂、设备繁多,如何理清其系统结构,识别主要脆弱点和关键薄弱点,影响并制约着供电任务的平稳实施。传统的静态可靠性分析方法只能应用于固定时刻的系统可靠性分析,在所有电源设备集成为一个动态系统时,其复杂结构、高危高压、动态失效的运行特点,使得该系统可靠性问题需采用动态的分析方式。

    动态系统可靠性分析方法是近年的研究热点及难点,DBN从这些方法中脱颖而出。机器学习[3]、数据挖掘[4]等技术的兴起也为DBN的发展和应用提供了更为广阔的空间。目前,DBN的理论成果得到了进一步的发展完善,包括因果推断[5]、不确定性知识表达[6]、模式识别[7]和聚类分析[8]等。DBN也因其在双向推理及故障诊断方面的优势,被广泛应用于复杂系统的动态可靠性分析。文献[9]通过引入$\,\beta $因子,结合模糊DBN,提出了一种针对数据缺失的多态系统可靠性分析方法。文献[10]基于贝叶斯网络对海上浮式风机进行了可靠性分析,大幅提升了故障诊断的准确率。文献[11]利用隐马尔科夫模型度量证据节点的先验概率,结合DBN,实现了化工设备的任务可靠性预计与评估。

    现有的DBN算法主要分为:离散时间贝叶斯网络算法(discrete-time Bayesian network, DTBN)[12]与连续时间贝叶斯网络算法(continues-time Bayesian network, CTBN)[13]。其中,DTBN以离散任务时间的方式,通过定义节点的条件概率表(conditional probability table, CPT)及边缘概率表(marginal probability table, MPT),求解出系统在不同任务时间片的后验概率。其中,文献[14]基于DTBN与改进的GO-FLOW方法,研究分析了多阶段多状态系统的共因失效问题。文献[15]通过引入区间分析理论,对无人机中表决系统的不确定性问题进行了量化处理,并基于DTBN算法对该系统进行了可靠性建模分析。然而,随着系统节点数与时间片数的增加,CPT&MPT的维度呈指数增加[16],导致求取系统可靠度所需的计算量也指数增长。

    CTBN可在连续任务时间条件下,建立各个节点的概率密度函数(probability density function, PDF)的解析解,并借此得出任意时刻节点的后验概率,极大地减少了计算时间。文献[17]基于CTBN构建车辆系统的性能函数,并将其应用于该系统的可靠性优化设计。文献[18]通过模糊函数量化了失效数据的参数不确定性,并针对DBN中的动态逻辑门,构建了基于CTBN算法的模糊函数解析解模型。然而,现有的CTBN算法,需要针对不同的分析对象建立特定的分析模型,缺少通用的建模方法[19]

    因此,本文提出一种改进的DBN概率表建模方法,在无需离散任务时间的情况下,实现电源系统的动态可靠性分析与故障诊断。

    • 在DTBN算法中,任务时间被离散为多个时间片段,且节点的条件概率与边缘概率以CPT和MPT的形式给出。以2时间片段的DTBN为例,“与”门的CPT&MPT如图1所示。对于节点AB${P_1}$${P_2}$分别代表节点在对应时间片中的失效概率,而状态3的概率$1 - {P_1} - {P_2}$表示该节点在任务时间内未发生失效的概率。而在“与”门的CPT中,$P\left( {T{\rm{ = }}1\left| {A = 1,B = 1} \right.} \right) = 1$表示节点AB在第一个时间片中失效时,该“与”门在此时间片内失效的条件概率为1,其他的概率也类似定义。基于图1的CPT,可通过式(1)~式(3)求解出节点T在各个状态的边缘概率。

      图  1  CTBN中“与”门的CPT&MPT

      $$\begin{split} P\left( {T = 1} \right) &= \sum\limits_{i = 1}^3 {\sum\limits_{j = 1}^3 {P\left( {\left. {T = 1} \right|A = i,B = j} \right)} } \times \\ & P\left( {A = i} \right)P\left( {B = j} \right) = {P_1}^2 \\ \end{split} $$ (1)
      $$\begin{split} P&\left( {T = 2} \right) = \sum\limits_{i = 1}^3 {\sum\limits_{j = 1}^3 {P\left( {\left. {T = 2} \right|A = i,B = j} \right)} } \times \\ &P\left( {A = i} \right)P\left( {B = j} \right) = {P_1}{P_2}{\rm{ + }}{P_2}{P_1}{\rm{ + }}{P_2}^2 \end{split} $$ (2)
      $$\begin{split} P \left( {T = 3} \right) = \sum\limits_{i = 1}^3 {\sum\limits_{j = 1}^3 {P\left( {\left. {T = 3} \right|A = i,B = j} \right)} } \times \quad \\ P\left( {A = i} \right)P\left( {B = j} \right) = \left( {1 - {P_1} - {P_2}} \right)\left( {1 + {P_1} + {P_2}} \right) = \\ 1 - {\left( {{P_1} + {P_2}} \right)^2} = 1 - P\left( {T = 1} \right) - P\left( {T = 2} \right) \quad \end{split}$$ (3)
    • 与DTBN不同的是,在CTBN中,引入了单位阶跃函数和脉冲函数[13],所有节点的概率分布以PDF的形式给出。以图1中的“与”门为例,T${f_{\left. T \right|A,B}}\left( {\left. T \right|A,B} \right)$为:

      $${f_{\left. T \right|A,B}}\left( {\left. T \right|A,B} \right) = u\left( {b - a} \right)\delta \left( {t - b} \right) + u\left( {a - b} \right)\delta \left( {t - a} \right)$$ (4)

      式中,$u\left( x \right) = \left\{ {\begin{array}{*{20}{c}}{0\quad\;\;\;\; {\rm{if}}\;x < 0}\\{1/2\quad {\rm{if}}\;x = 0}\\{1\quad\;\;\;\; {\rm{if}}\;x > 0}\end{array}} \right.$

      $x \ne {\rm{0}}$$\delta \left( x \right) \!=\! 0$,且$\displaystyle\int\nolimits_0^\infty {\delta \left( t \right){\rm{d}}t} \!=\! 1$$\displaystyle\int\nolimits_0^\infty {f\left( t \right)\delta \left( {t - \tau } \right)} $$ {{\rm{d}}x} \!=\! f\left( \tau \right)$

      则求解出节点T的PDF为:

      $$\begin{split} {f_T}\left( t \right) = \int_0^\infty {\int_0^\infty {{f_{\left. T \right|A,B}}\left( {\left. T \right|A,B} \right){f_B}\left( b \right){f_A}\left( a \right){\rm{d}}a{\rm{d}}b} } = \\ \int_0^\infty {\int_0^\infty {\left[ {u\left( {b - a} \right)\delta \left( {t - b} \right) + u\left( {a - b} \right)\delta \left( {t - a} \right)} \right]}}\times\;\;\\ {{{f_B}\left( b \right){f_A}\left( a \right){\rm{d}}a{\rm{d}}b} } = {f_B}\left( t \right){F_A}\left( t \right) + {f_A}\left( t \right){F_B}\left( t \right) =\;\;\;\\ [{F_A}\left( t \right){F_B}\left( t \right)]' \qquad\qquad\qquad \end{split} $$ (5)

      但是,正如引言中介绍的,现有的CTBN算法,需要针对不同的分析对象建立特定的解析解模型。并且随着动态逻辑门和节点数的增加,该方法的建模难度也会随之增加。因此,本文提出一种改进的CPT&MPT建模方法。

    • 以“与”门的CPT&MPT为例,如图2所示,$P\left( {A = 1} \right) = {F_A}\left( t \right)$表示节点A在该时间点的失效概率,而$P\left( {A = 2} \right) = {R_A}\left( t \right) = 1 - {F_A}\left( t \right)$表示节点A在该时间点的可靠度。通过图2中节点T的CPT,则可由式(6)得出该节点的边缘概率分布。

      $$\begin{split} {F_T}\left( t \right) = &P\left( {T = 1} \right) = \sum\limits_{i = 1}^2 {\sum\limits_{j = 1}^2 } {{P\left( {\left. {T = 1} \right|A = i,B = j} \right)} } \times \\ &P\left( {A = i} \right)P\left( {B = j} \right) = {F_A}\left( t \right){F_B}\left( t \right) \\ \end{split}$$ (6)

      图  2  “与”门的CPT&MPT

      可以看出,该CPT与传统CTBN算法中的单位阶跃函数和脉冲函数有相同的功能。且通过算法1,可以计算出任意数量节点的“与”门CPT,且该“与”门的父节点可以服从任意的失效分布类型。

      算法 1 由$i$个节点组成的“与”门的$\rm{CPT}$计算

      Input: $i$← the number of parent nodes

      $\rm{CPT}$← zeros($\left[ {2,{2^i}} \right]$)

      for ${P_1}$ ← 1 to 2

      for ${P_2}$$\left( {\left( {{P_1} - 1} \right) \times 2 + 1} \right):{P_1} \times 2$

      for ${P_i}$$\left( {\left( {{P_{i - 1}} - 1} \right) \times 2 + 1} \right):{P_{i - 1}} \times 2$

      $h = \max \left[ {{P_1},{P_2} - \left( {{P_1} - 1} \right) \times 2,\cdots,{P_i} - \left( {{P_{i - 1}} - 1} \right) \times 2} \right]$

      ${\rm{CPT}}\left( {\left[ h,{P_i} \right]} \right) = 1$

      end for

      end for

      end for

      Output $\rm{CPT}$

      而针对动态贝叶斯网络中的动态逻辑门,以温备份门为例(warm spare, WSP),如图3所示,节点B为节点A的温备份节点。当节点A正常运行时,节点B处于温备份状态,其失效率${\lambda _B}$下降为$\alpha {\lambda _B}$$\alpha $为备份因子,且$0 \leqslant \alpha \leqslant 1$。当$\alpha = 0$时,节点B转变为冷备份节点;$\alpha = 1$时,节点B转变为热备份节点。当A失效时,B转变为工作状态,且失效率转变为${\lambda _B}$。则节点B的条件概率可由式(7)~式(10)得出。

      图  3  温备份门的CPT&MPT

      $$P\left( {B = 1\left| {A = 1} \right.} \right) = \frac{{\displaystyle\int_0^t {{f_A}\left( {{t_A}} \right)\int_0^{{t_A}} {{f_{\alpha B}}\left( {{t_B}} \right){\rm{d}}{t_B}{\rm{d}}{t_A} + \int_0^t {{f_A}\left( {{t_A}} \right)\left[ {\left( {1 - \int_0^{{t_A}} {{f_{\alpha B}}\left( {{t_B}} \right){\rm{d}}{t_B}} } \right)\int_{{t_A}}^t {{f_B}\left( {{t_B} - {t_A}} \right)} {\rm{d}}{t_B}} \right]{\rm{d}}{t_A}} } } }}{{\displaystyle\int_0^t {{f_A}\left( {{t_A}} \right){\rm{d}}{t_A}} }}$$ (7)
      $$P\left( {B = 2\left| {A = 1} \right.} \right) = \frac{{\displaystyle\int_0^t {{f_A}\left( {{t_A}} \right)\left[ {\left( {1 - \int_0^{{t_A}} {{f_{\alpha B}}\left( {{t_B}} \right){\rm{d}}{t_B}} } \right)\int_t^{ + \infty } {{f_B}\left( {{t_B} - {t_A}} \right)} {\rm{d}}{t_B}} \right]{\rm{d}}{t_A}} }}{{\displaystyle\int_0^t {{f_A}\left( {{t_A}} \right){\rm{d}}{t_A}} }}$$ (8)
      $$P\left( {B = 1\left| {A = 2} \right.} \right) = \displaystyle\int_0^t {{f_{\alpha B}}\left( {{t_B}} \right){\rm{d}}{t_B}} $$ (9)
      $$P\left( {B = 2\left| {A = 2} \right.} \right) = \displaystyle\int_t^\infty {{f_{\alpha B}}\left( {{t_B}} \right){\rm{d}}{t_B}} $$ (10)

      式中,${f_A}\left( t \right)$${f_B}\left( t \right)$分别为节点AB工作时失效概率的PDF;${f_{\alpha B}}\left( t \right)$为节点B在温备份状态下的PDF;${t_A}$${t_B}$分别为节点AB的失效时间。针对电子系统,当节点AB的失效类型服从指数分布时,即${f_A}\left( t \right) = {f_B}\left( t \right) = \lambda {{\rm{e}}^{ - \lambda t}}$${f_{\alpha B}}\left( t \right) = \alpha \lambda {{\rm{e}}^{ - \alpha \lambda t}}$,则节点B的CPT如表1所示。

      表 1  温备份节点B的CPT

      A 12
      B1$1 - \dfrac{{{{\rm{e}}^{ - \lambda t}}\left( {1 - {{\rm{e}}^{ - \alpha \lambda t}}} \right)}}{{\alpha \left( {1 - {{\rm{e}}^{ - \lambda t}}} \right)}}$$1 - {{\rm{e}}^{ - \alpha \lambda t}}$
      2$\dfrac{{{{\rm{e}}^{ - \lambda t}}\left( {1 - {{\rm{e}}^{ - \alpha \lambda t}}} \right)}}{{\alpha \left( {1 - {{\rm{e}}^{ - \lambda t}}} \right)}}$${{\rm{e}}^{\left( { - \alpha \lambda t} \right)}}$

      为验证所提方法的准确性,以图3中的温备份门(节点T)为例,将各节点的CPT&MPT代入算法2中,当${\lambda _A} = {\lambda _B} = 0.000\;01$$\alpha = 0.5$时,可计算出T节点的可靠度曲线,分别与蒙特卡罗仿真(Monte Carlo simulation, MCS)[18]以及DTBN的计算结果进行了对比,如图4所示。

      图  4  结果对比分析

      算法 2:DBN后验概率求解算法

      Input: $N$ ← number of nodes

      dag($N$,$N$) ← adjacent matrix for DBN

      state number of each node ← $2\times {\rm{ones}}\left( {1,N} \right)$

      for $i$ ← 1 to 50 do

      $t$$i \times 10\;000$

      bnet.CPD{$N$} ←CPT&MPT;

      define inference engine;

      define evidence;

      end for

      Output posterior probabilities of each node

      图4结果可知,所提方法可在未离散任务时间的情况下,获取与DTBN相同的系统可靠度计算结果。而且,本方法B节点的CPT维度仅为$\left( {2 \times 2} \right)$,即便是与最简单的2时间片段DTBN相比,其建模难度也远小于DTBN的CPT[16]。与MCS的计算结果的对比分析结果也验证了所提方法的准确性。

      表1的CPT可求解出B节点${F_B}\left( t \right)$的解析解,如式(11)所示,其结果与CTBN[19]一致。并且,通过将构建的CPT&MPT与算法2相结合,即可计算得到观测节点的后验概率,而无需针对不同分析对象建立特定的解析解模型。

      $$\begin{split} {F_B}\left( t \right) = &P\left( {B = 1} \right) = \sum\limits_{i = 1}^2 {P\left( {\left. {B = 1} \right|A = i} \right)} P\left( {A = i} \right) = \\ &1 - \frac{1}{\alpha }{{\rm{e}}^{ - \lambda t}} + \frac{{1 - \alpha }}{\alpha }{{\rm{e}}^{ - \lambda \left( {1 + \alpha } \right)t}} \end{split} $$ (11)
    • 本文以某电源系统为研究对象,该系统组成结构如图5所示。

      图  5  电源系统结构框图

      该系统通过将民用电转换成高压直流电,给规定容量的主预储能电容器充电。在规定的时间范围内达到指定值,在规定的时序触发下,主、预储能电容通过晶闸管固态开关、平波电感给负载放电,为负载提供合适的能量脉冲,使它们的形状、幅度和时序保持一致。根据该电源系统的结构框图,可建立其DBN如图6所示,事件编号及失效率如表2所示。

      图  6  电源系统的DBN

      表 2  事件节点编号及失效率

      节点
      编号
      名称节点
      编号
      名称失效率(${\lambda _i}$)/(10−6/h)
      TM电源系统失效AM配电单元2.67900
      BM1充电失效EM调波单元0.01120
      BM2充电模块失效B1低压供电模块0.29400
      CM1储能失效B2主充电模块2.16040
      CM2储能电容失效B3预充电模块1.08020
      DM开关失效B4控制采集模块7.62380
      C1主储能电容模块0.00580
      C2预储能电容模块0.00058
      C3泄放单元模块0.65470
      C4充电保护单元模块0.00780
      D1主开关单元模块7.53300
      D2预开关单元模块5.27310
    • 该电源系统的DBN模型中,包括3个WPS门(BM2、CM2、DM),其组成节点的CPT&MPT如图3所示,节点B3、C2、D2的CPT可由表1计算。将各节点的CPT&MPT代入算法2中,即可求解出各个节点的后验概率。各节点状态为2时的边缘概率即为该节点对应事件的可靠度值。通过算法2,可以求解出电源系统顶事件及中间事件在任务时间内(5×105 h)的可靠度变化曲线,如图7所示。

      图  7  电源系统可靠度曲线

      其中,节点AM为底事件节点,且服从指数分布,所以其可靠度曲线是典型的指数分布曲线。根据表2中的失效率数据,因为与节点EM、CM1、CM2相关联底事件的失效率太低,所以在任务时间内上述3个节点的可靠度没有明显变化。其他的节点(TM、BM1、BM2、DM)则由各种逻辑门组合而成,因可靠度不服从特定分布,所以曲线呈现出图7所示的变化趋势。

      基于算法2,可以求解出顶事件TM状态为1时,其他节点i状态为1的后验概率$P\left( {i = 1\left| {\rm{TM} = 1} \right.} \right)$(i=AM、BM1、BM2、CM1、CM2、DM、EM)。即电源系统失效时,其他节点失效的条件概率,该条件概率可用作系统失效下部件故障的诊断分析。诊断结果如图8所示,在0.875×105 h前,节点AM的后验失效概率最大,此时段若系统失效,则配电单元失效的概率最大。而在(0.875~2.9)×105 h之间,则是DM节点(开关失效)的后验失效概率最大。在1.925×105 h时刻,BM2(充电模块失效)的后验失效概率超过AM节点,且在2.9×105 h之后成为后验失效概率最大的节点。

      图  8  电源系统故障诊断

    • 本文针对DBN中的“与”、“或”、WSP逻辑门,提出了一种改进的DBN概率表建模方法,并应用于电源系统可靠性分析与故障诊断,由前面分析可得出以下结论:

      1) 与DTBN相比,本文所提出的概率表建模方法可在无需离散任务时间的条件下建立节点的CPT&MPT,降低建模难度。而且,所构建的CPT&MPT可单独赋值于各节点,减少针对特定系统建模的工作量。

      2) 结合文中所提概率表建模方法及边缘概率求解算法,可实现复杂系统任务时间内各个事件的动态可靠性分析,并同时求解出所有观察节点的可靠度曲线。

      3) 通过贝叶斯双向推理,可求解出顶事件发生时,各子系统失效的后验概率。基于所求后验概率,可实现系统任务时间内的动态故障诊断,评估各时段内子系统失效的后验概率值,定位出最可能失效的子系统及部件。

参考文献 (19)

目录

    /

    返回文章
    返回