基于MCMC方法的船舶柴油机磨损故障分布研究
时间：2020-12-15 04:01:47 来源：雅意学习网 本文已影响 人 
摘要：该文借助WINBUGS软件运用MCMC（MarkovChainMonteCarlo，马尔科夫链蒙特卡洛）估计方法，采用Gibbs抽样原理，在贝叶斯理论的框架下，结合某船舶柴油机气缸盖历史的和更新的磨损故障时间数据，有效提高了磨损故障概率分布的精度。首先介绍了柴油機对船舶航行安全的重要性及其常见故障类型，确定了采用贝叶斯框架下的MCMC抽样方法对磨损故障分布进行仿真计算；其次阐述了通过最大似然估计法获得先验分布参数的过程，并说明了MCMC方法的原理及其常用的抽样方法；最后，引
WINBUGS软件并采用此软件进行建模分析，进行基于Gibbs抽样的MCMC方法的柴油机磨损故障分布计算，以实际算例证明此方法提高了柴油机磨损故障概率分布的计算精度。
关键词：柴油机；磨损故障；贝叶斯网络；MCMC；WinBUGS
中图分类号：TP311
文献标识码：A
文章编号：1009-3044（2017）10-0190-05
1.概述
柴油机是船舶动力装置的关键设备，一旦其发生故障，将会严重影响船舶的正常运行，并有可能直接或间接地造成无法估量的经济损失，甚至可能损坏其他关键设备，严重时还会危及到船员和旅客的人身安全。船舶柴油机的工作环境和条件往往多变而恶劣，由于不稳定的环境与人为操作，船舶柴油机常在负荷与转速多变的工况下运转。柴油机的结构复杂性决定了柴油机故障类型的复杂性，负荷与转速多变又会直接导致柴油机的许多零部件常处于高温、高压、高速的环境下工作，加之人为因素导致的润滑不足、操作不规范等，柴油机零部件故障多发、复杂，具有不确定性和多样性。
船舶柴油机的故障按照子系统来分，主要分为喷油系统故障、燃料系统故障、涡轮增压系统故障、冷却系统故障以及润滑系统故障等。如果按照导致柴油机故障的原因来分类，那么主要有磨损、冲击振动、腐蚀、老化等。其中，磨损故障时柴油机的主要故障类型，约占船舶柴油机故障比率的37.5%t31。船舶柴油机磨损故障一旦发生，若没有及时、准确地排故和维修，就无法保障船舶动力装置的正常运行，甚至会导致严重的安全事故和经济损失。但是，过于频繁地对船舶柴油机进行检修，不仅会带来不必要的人为差错，还会导致维护费用急剧上升，不符合船舶管理的经济性要求。因此，掌握柴油机磨损故障的发生规律、合理地安排检修工作至关重要。
国内外在这一领域有不少相关研究，但大多集中在柴油机故障诊断技术可靠性分析方面，对磨损故障的具体分布研究不够细化，本文在这方面进行了一定的研究，结合MC-MC方法有效提高了船舶柴油机磨损故障分布的计算精度。
柴油机磨损故障的失效前平均时间（Mean Time To Fail-tire，MTI"F）的概率分布可以通过历史的统计数据获得，但由于船舶柴油机工作环境复杂多变，根据历史数据计算得出的概率分布仅能代表综合环境下的情况，对于近期具体的工作环境下的故障分布来说不够准确。如果能结合近期的MTTF数据，对由历史数据获得的柴油机磨损故障MTTF分布进行更新，就可以得到更加符合现有实际情况的概率分布，有效提高磨损故障发生的预测精度。
在信息更新方面，贝叶斯网络（Bayesian Network，BN）具有独特的优势。贝叶斯网络是表示随机变量间依赖和独立关系的网络模型，也称为贝叶斯网、因果概率网络或信念网络，它是一个有向非循环的图形结构，其中节点代表问题域的随机变量，节点间的有向边代表变量之间的直接依赖关系，每个节点都附有一个概率分布。贝叶斯定理可以结合先验信息和样本信息，获得更新的后验分布，提高预测精度，如图1所示。但是，贝叶斯统计的困难在于，后验分布是复杂的、高维的分布，难以进行直接计算，所以随着计算机技术与MCMC方法的发展，后验分布的计算被大大简化。
常见的船舶柴油机零部件故障分布规律有指数分布、正态分布、指数正态分布和威布尔分布。超载下工作或过热导致的柴油机故障、人为失误以及偶然操作不当引起的突发故障等均适用于指数分布的故障分布规律；由几种相对独立、作用均匀的微小差异量因素造成的零部件故障，如燃烧室组件、气缸和活塞等的磨损、喷油系统沉淀故障等，其故障概率分布函数大多为正态分布；油水系统和齿轮传动系统故障、滚珠轴承的磨损故障等服从威布尔分布规律。由于不同故障分布问题的在本文的研究方法下类似，所以，选取正态分布进行后续计算说明。
3.后验分布参数获取
在贝叶斯估计中，后验均值即为参数的估计。但在实际计算时，参数后验分布往往难以计算，解析法求解几乎不可能，计算量庞大而过程复杂。近十几年来，计算机技术高速发展，加之贝叶斯方法的不断改进完善，特别是MCMC方法以及Win-BUGS（Bayesian inference Using Gibbs Sampling）软件的发展和应用，原本复杂的高维计算问题可以通过从后验密度抽样取平均的方法方便快速地得到解决。原本传统蒙特卡洛的基本思想是，当所要求解的问题是某种事件出现的概率，或者是某个随机变量的期望值时，它们可以通过某种“试验”的方法，得到这种事件出现的频率，或者这个随机变数的平均值，并用它们作为问题的解。在处理高维、复杂的问题时，MCMC方法相对于经典的蒙特卡洛方法有明显的优越性，可以在非标准形式、各变量之间互相不独立的情况下进行分布的模拟，在近十几年内被广泛应用于贝叶斯统计、显著性检验等方面。
3.1 MCMC方法
MCMC方法的核心思想是反复对目标分布进行抽样，得到目标分布的一组样本数据，然后再利用抽样得到的样本对目标分布的性质进行分析和讨论。也就是说，通过建立一个平稳分布为π（x）的马尔科夫链来得到π（x）的样本，基于这些样本对目标分布进行各种统计推断。
在MCMC方法中，可以根据转移核的不同构造方法，将主要的MCMC方法分为三种：
