一种基于分块的船舶柴油机故障监测方法，属于船舶柴油机故障监测领域，包括：在船舶柴油机各子系统安装传感器监测运行状态，同时根据结构和工作原理对监测变量进行分块；分别构建每一分块对应的主成分分析监控空间和支持向量机分类模型；使用T2和SPE统计量实时监测每一分块中的运行数据，将超出对应控制限的监测变量导入构建的对应支持向量机分类模型中；使用支持向量机确定滑动窗口内出现故障的概率，根据故障阈值确定故障是否发生；最后根据故障贡献率确定故障发生分块。本发明使用分块方法可以提升故障监测过程的计算效率，根据主成分分析和支持向量机实现二次故障检测降低故障误判，最后利用故障贡献率确定故障发生位置，降低船舶柴油机维护难度。
1.一种基于分块的船舶柴油机故障监测方法，其特征在于，该方法根据船舶柴油机的结构和工作原理对其监测变量进行分块，分别使用主成分分析方法和支持向量机分类模型实现每一分块监测数据的二次故障检测，随后根据滑动窗口内出现故障的概率和人为设定的故障阈值确定故障发生的可能性以减小故障误判，最后利用故障贡献率确定故障发生的具体分块；具体包括以下步骤：步骤S1：在船舶柴油机重要节点安装状态参数监测传感器，并对监测变量进行分块；
对柴油机重要节点的状态参数进行采集和处理，得到监测变量；所述的状态参数包括温度、压力、各种介质的质量流量及各子系统的工作效率；所述的重要节点包括不同介质流经各子系统前后位置节点和各个工作步骤的发生位置节点，其中所述介质包括空气、冷却水及滑油，工作步骤包括燃烧过程、冷却过程及空气增压过程，子系统包括涡轮器、增压器、机械增压器、空冷器、滤清器、滑油冷却器、扫气箱及气缸；
基于对监测变量的分块实现故障监测，分块方法是根据船舶柴油机的结构和工作原理实现，将每一个单独的工作步骤获得的监测变量分为一块；
步骤S2：在每一分块中使用主成分分析方法将监测变量从原始监测空间映射到主成分分析监测空间，并同时训练获得支持向量机故障分类模型；
所述的主成分分析监测空间可以表示为不同分块中主元空间和残差空间的负载矩阵；
所述的支持向量机故障分类模型是根据历史正常和故障数据确定的可以实现正常监测变量和故障监测变量分类的函数；
步骤S3：分别在每一分块的主元空间和残差空间中构建两个监控统计量及其对应的控制限，实现监测变量的一次故障检测，随后将检测出可能存在故障的监测变量传递给对应分块的支持向量机分类模型；
其中，T2用于表征当前主元空间中监测数据与正常数据的距离，SPE用于表征残差空间中监测变量之间关系的变化程度；X表示分块中的监测数据矩阵，Pt表示主成分分析方法中主元空间的负载矩阵，Pe表示主成分分析方法中残差空间的负载矩阵；
分析获得两种监控统计量的正常情况下的分布规律如下式所示：
其中，F为统计学中的F分布，β是置信度系数，两个统计量可以根据其取值确定对应的控制限；χ2为统计学中的卡方分布，g＝μ/δ，h＝2δ2/μ，其中μ为正常情况下监控统计量的均值，δ为正常情况下监控统计量的标准差；n表示监测数据矩阵X的行数，Jt表示负载矩阵Pt的列数；当实时监测数据的监控统计量均超出控制限时说明当前数据中可能存在故障因素，完成一次故障检测；将检测出的可能存在故障的监测数据传递给对应分块的支持向量机分类模型，继续监测后续数据；
步骤S4：当连续t个采样点的监控统计量均超出控制限时则使用对应分块的支持向量机分类模型对步骤S3中传递过来的监测数据进行二次故障检测，使用设定的长度为N的滑动窗口中出现故障因素的概率和人为设定的故障阈值最终确定故障是否发生；
所述的连续监控统计量超出控制限的采样点个数t必须大于或等于滑动窗口长度N，因此认为只有当t≥N时才开始二次故障检测；所述的二次故障检测是指当连续n个采样点的监测变量均可以被主成分分析方法检测出故障因素时，将这些数据输入到对应分块的支持向量机分类模型中再次进行故障检测，判断是否可能存在故障因素；随后建立一个长度为N的滑动窗口，当滑动窗口内出现故障的概率超出人为设定的阈值时即可确定系统此时存在故障；所述的出现故障的概率是指滑动窗口内检测出故障的采样点个数与滑动窗口长度的比值；
步骤S5：构建故障贡献率参数，分别计算出现故障时每一分块的故障贡献率，确定故障发生的具体分块；
其中，下标b是指具体的分块编号， 为第b个分块的T2统计量； 为对应的第b个分块的T2统计量控制限；SPEb为第b个分块的SPE统计量；SPEb,lim为对应的第b个分块的SPE统计量控制限；
使用T2和SPE统计量超出控制限的倍数作为此时对应分块的故障贡献率，确定其中故障贡献率最高的分块，选定故障贡献率最大的分块即为故障发生分块。
2.根据权利要求1所述的一种基于分块的船舶柴油机故障监测方法，其特征在于，所述的步骤S2中：所述的主成分分析方法对分块中的监测数据矩阵X的处理方式如下式所示，其中监测数据矩阵X的维度为n×J：T＝XPt
其中，T为得分矩阵，Pt是主元空间的负载矩阵，维度为J×Jt，即矩阵行数为J，列数为Jt，E为残差空间，Pe是残差空间中的负载矩阵，维度为J×Je，即矩阵行数为J，列数为Je；根据累计方差贡献率准则选定主元空间中负载矩阵，选定的累计方差贡献率为85％时可以得到较好的结果；由于这些负载矩阵是主成分分析监测空间的方向向量，因此所述的主成分分析监测空间可以表述为根据上述说明确定的每一分块中主元空间和残差空间的负载矩阵Pt和Pe；
其中，αk＝[α1，α2，…，αn]是支持向量机分类函数的系数矩阵，c是对应的偏差量，K(xi,xj)＝exp(-||xi-xj||2/2σ2)是选定的高斯核函数，σ为xi和xj的协方差矩阵；这里是利用每一分块中历史正常数据和故障数据确定的系数αk和c，获得故障分类模型即为根据历史正常和故障数据确定的可以实现正常监测变量和故障监测变量分类的函数。
一种基于分块的船舶柴油机故障监测方法
技术领域
[0001] 本发明属于船舶柴油机故障监测领域，涉及一种分块实现船舶柴油机故障监测的方法。
背景技术
[0002] 船舶柴油机是当前最主要的船舶动力装置，在完整的船舶运输中占据十分重要的地位。其故障占完整船舶故障的60％～80％，而当船舶发生故障时会造成巨大的经济损失，严重的甚至会影响操作人员的人身安全。此外，随着“工业4.0”背景下的“E-航海”战略的实施，人们对船舶长期稳定运行提出了更高的要求，因此十分有必要实时监控船舶柴油机健康状态以保证其正常的工作。基于数据驱动的故障监测方法主要是发现和监测由故障引起的监测变量波动，进而实时监测其健康状态，由于这种方法不需要过多的船舶柴油机机理知识和专家经验，已经受到了人们的广泛关注。
[0003] 船舶柴油机是典型的多系统、多层次的复杂系统，各种监测变量携带不同的系统健康状态信息，为尽量的扩大故障监控范围并提升故障检测精度，通常需要同时监测多个变量，因此可以认为多种不同的监测变量是实现基于数据驱动的故障监测技术的基础。然而并不是每一个监测变量均与所有的故障相关，大量监测变量的存在也必然会引起计算复杂度的增加。同时某些监测变量之间关联关系较小，没有必要同时处理，因此发明一种基于分块的船舶柴油机故障监测方法，用以高效精确的确定故障发生位置具有十分重要的理论和实际意义。
[0004] 经对现有技术的文献检索发现，发明名称为“基于离群分析的柴油机故障诊断方法”(公开号为CN103969052A，公开日为2014.08.06)的专利中提出一种使用群体智能实现故障诊断的方法。该专利申请自述为：“本发明包括如下步骤：1)实时采集柴油机运行数据，包括各气缸排气温度、各气缸的爆发压力及各缸发货时刻对应的转速；2)各气缸参数标准化处理；3)对气缸运行状态进行聚类分析；4)通过计算异常状态气缸的离群因子定量的描述其故障程度”。其不足之处是：选用的监测变量较少，易受噪声因素的影响，此外该故障诊断方法不利于大型柴油机故障因素的确定，因此难以实现高精度的复杂船舶柴油机系统的故障监测。
发明内容
[0005] 本发明致力于解决现有船舶柴油机存在大量监测变量引起的故障监测过程计算效率较低的问题，提出一种基于分块的船舶柴油机故障监测方法，可以在保证故障检测精度的同时提升计算效率。
[0006] 为了达到上述目的，本发明采用的技术方案是：
[0007] 一种基于分块的船舶柴油机故障监测方法，该方法根据船舶柴油机的结构和工作原理对其监测变量进行分块，分别使用主成分分析方法和支持向量机分类模型实现每一分块监测数据的二次故障检测，随后根据滑动窗口内出现故障的概率和人为设定的故障阈值确定故障发生的可能性以减小故障误判，最后利用故障贡献率确定故障发生的具体分块。上述人为设定的故障阈值是根据具体的故障情况和故障危害程度一起确定的，当故障危害严重或故障发生频繁时可以选取较小的故障阈值以及时发现故障，当要求故障误判率低时则可以选择较大的故障阈值。具体包括以下步骤：
[0008] 步骤S1：在船舶柴油机重要节点安装状态参数监测传感器，并按照结构和工作原理对监测变量进行分块。
[0009] 为扩大本发明所提方法实现的故障监控范围，需要对柴油机重要节点的状态参数进行采集和处理。所述的状态参数主要包括温度、压力、各种介质的质量流量及各子系统的工作效率等。所述的重要节点主要包括不同介质流经各子系统前后位置节点和各个工作步骤的发生位置节点，其中所述的介质主要包括空气、冷却水及滑油等，工作步骤包括燃烧过程、冷却过程及空气增压过程等，所述的子系统包括涡轮器、增压器、机械增压器、空冷器、滤清器、滑油冷却器、扫气箱及气缸等。对上述重要节点的状态参数进行采集和处理之后即可获得监测变量。
[0010] 本发明主要是基于对监测变量的分块实现故障监测，所述的分块方法主要是根据船舶柴油机的结构和工作原理实现的。由于每一单独的工作步骤均会引起对应节点的状态参数变化而对其他节点的状态参数影响较小，因此可以将每一个单独的工作步骤获得的监测变量分为一块，即最简单的分块方案就是分别将每一子系统的监测变量分为一块。
[0011] 步骤S2：在每一分块中使用主成分分析方法将监测变量从原始监测空间映射到主成分分析监测空间，并同时训练获得支持向量机故障分类模型。
[0012] 所述的主成分分析方法对分块中的监测数据矩阵X(维度为n×J，即矩阵行数为n，列数为J)的具体处理方式如下式所示：
[0013] T＝XPt
[0014]
[0015] 其中，T为得分矩阵，Pt是主元空间的负载矩阵，维度为J×Jt(即矩阵行数为J，列数为Jt)，E为残差空间，Pe是残差空间中的负载矩阵，维度为J×Je(即矩阵行数为J，列数为Je)；根据累计方差贡献率准则选定主元空间中负载矩阵，一般认为选定的累计方差贡献率为85％时可以得到较好的结果。由于这些负载矩阵是主成分分析监测空间的方向向量，因此所述的主成分分析监测空间可以表述为根据上述说明确定的每一分块中主元空间和残差空间的负载矩阵Pt和Pe。
[0016] 所述的支持向量机分类函数表达式如下所示：
[0017]
[0018] 其中，αk＝[α1,α2,…,αn]是支持向量机分类函数的系数矩阵，c是对应的偏差量，K(xi,xj)＝exp(-||xi-xj||2/2σ2)是选定的高斯核函数，σ为xi和xj的协方差矩阵。这里主要是利用每一分块中历史正常数据和故障数据确定的系数α和c，获得故障分类模型即为根据历史正常和故障数据确定的可以实现正常监测变量和故障监测变量分类的函数。
[0019] 步骤S3：分别在每一分块的主元空间和残差空间中构建两个监控统计量及其对应的控制限，实现监测变量的一次故障检测，随后将检测出可能存在故障的监测变量传递给对应分块的支持向量机分类模型。
[0020] 所述的监控统计量构建公式如下式所示：
[0021] T2＝(XPt)Λ-1(XPt)T
[0022]
[0023] 其中，T2用于表征当前主元空间中监测数据与正常数据的距离，SPE用于表征残差空间中监测变量之间关系的变化程度。分析获得两种监控统计量的正常情况下的分布规律如下式所示：
[0024]
[0025]
[0026] 其中，F为统计学中的F分布，β是置信度系数，两个统计量可以根据其取值确定对应的控制限。χ2为统计学中的卡方分布，g＝μ/δ，h＝2δ2/μ，其中μ为正常情况下监控统计量的均值，δ为正常情况下监控统计量的标准差。当实时监测数据的监控统计量均超出控制限时说明当前数据中可能存在故障因素，完成一次故障检测。将检测出的可能存在故障的监测数据传递给对应分块的支持向量机分类模型，继续监测后续数据。
[0027] 步骤S4：当连续t个采样点的监控统计量均超出控制限时则使用对应分块的支持向量机分类模型对步骤S3中传递过来的监测数据进行二次故障检测，使用设定的长度为N的滑动窗口中出现故障因素的概率和人为设定的故障阈值最终确定故障是否发生。
[0028] 所述的连续监控统计量超出控制限的采样点个数t必须大于或等于滑动窗口长度N，因此可以直接认为只有当t≥N时才会开始二次故障检测。所述的二次故障检测是指当连续n个采样点的监测变量均可以被主成分分析方法检测出故障因素时，将这些数据输入到对应分块的支持向量机分类模型中再次进行故障检测，判断是否可能存在故障因素。随后建立一个长度为N的滑动窗口，当滑动窗口内出现故障的概率超出人为设定的阈值时即可确定系统此时存在故障。所述的出现故障的概率是指滑动窗口内检测出故障的采样点个数与滑动窗口长度的比值。当此时的故障发生概率超出人为设定的故障阈值时即可确定故障已经发生。此外操作人员可以根据实际故障情况与故障监测要求调整滑动窗口长度与对应故障阈值获得更合理的故障检测误判率和漏判率。
[0029] 步骤S5：构建故障贡献率参数，分别计算出现故障时每一分块的故障贡献率，确定故障发生的具体分块。
[0030] 所述的故障贡献率主要用于表现不同分块中监测变量内的故障严重程度，具体的可以通过监控统计量超出控制限的倍数确定，因此具体的第b个分块的故障贡献率计算方法如下式所示：
[0031]
[0032] 其中，下标b是指具体的分块编号， 为第b个分块的T2统计量； 为对应的第b个分块的T2统计量控制限；SPEb为第b个分块的SPE统计量；SPEb,lim为对应的第b个分块的SPE统计量控制限。
[0033] 使用T2和SPE统计量超出控制限的倍数作为此时对应分块的故障贡献率，确定其中故障贡献率最高的分块，选定故障贡献率最大的分块即为故障发生分块。
[0034] 与现有技术相比，本发明的有益效果为：将分块理论应用于船舶柴油机故障监测领域，提升故障监测过程的计算效率；使用主成分分析和支持向量机两种方法实现故障的二次检测，降低故障误判率；使用滑动窗口方法评估出现故障的可能性，可以根据人为设定的故障阈值和实际要求调整故障检测结果的误判率和漏判率；使用贡献率方法确定船舶柴油机故障发生分块，降低故障维护难度。
附图说明
[0035] 图1是基于分块的船舶柴油机故障监测方法流程图；
[0036] 图2是六缸二冲程船舶柴油机分块情况说明；
[0037] 图3是涡轮增压器分块情况示意图；
[0038] 图4是多缸柴油机气缸分块情况示意图；
[0039] 图5是空冷器冷却不足故障不同分块故障贡献率示意图。
[0040] 图中：1涡轮增压器；2空冷器分块；3气缸；1-1增压器分块；1-2涡轮器分块；3-1气缸分块；3-2气缸分块；3-3气缸分块；3-4气缸分块；3-5气缸分块；3-6气缸分块。
具体实施方式
[0041] 以下结合空冷器冷却不足故障具体实施例对本发明做进一步说明，但本发明的保护范围不限于下述的实施例。完整的故障监测过程包括如下步骤：
[0042] 步骤S1：分别在船舶柴油机的不同子系统的不同位置装入传感器。以最简单的船舶柴油机为例，如图2图3和图4所示，对所述的船舶柴油机监测变量具体分块进行说明，当前柴油机主要由涡轮增压器1、空冷器2、气缸3组成。所述涡轮增压器1包括增压器1-1和涡轮器1-2两个分块，其中增压器1-1分块的监测变量包括流经增压器的空气质量流量、增压器前后空气温差、增压器前后空气压差和增压器转速；涡轮器1-2分块的监测变量包括涡轮器工作效率、涡轮器前后空气温差、涡轮器前后空气压差、流经涡轮器的空气质量流量。所述的空冷器2分块监测变量包括空冷器前后温差、空冷器前后空气压差、流经空冷器的空气质量流量和空冷器冷却液质量流量等。所述的气缸3是柴油机中燃烧过程发生的场所，可以被看作是整个柴油机的核心，因此将每个气缸作为单独的分块，这里以六缸船舶柴油机为例，气缸包括1号气缸3-1、2号气缸3-2、3号气缸3-3、4号气缸3-4、5号气缸3-5和6号气缸3-6六个分块，每个气缸分块中提取的监测变量包括扫气温度、扫气压力、最大爆发压力、排气温度、排气压力、滑油流量和滑油温度。
[0043] 步骤S2：构建不同分块对应的主成分分析监测空间，本实施例中所述的增压器1-1分块中主元空间负载矩阵为 残差空间负载矩阵为所述的涡轮器1-2分块中主元空间负载矩阵为
残差空间负载矩阵为
所述的空冷器2分块中主元空间负载矩阵为[0.6650 0.7407 0.0832 0.0469]T，残差空间负载矩阵为 所述的1号气缸3-1分块选用主元空间负载矩阵
为 残差空间负载矩阵为
所述的2号气缸3-2分块中主元空
间负载矩阵为 残差空间负
载矩阵为 3号气缸3-3分
块中主元空间负载矩阵为
残差空间负载矩阵为 4号气缸3-4分
块中主元空间负载矩阵为 残
差空间负载矩阵为 5号气缸3-5分
块主元空间负载矩阵为 残差
空间负载矩阵为 6号气缸3-6分
块中主元空间负载矩阵为 残
差空间负载矩阵为
[0044] 步骤S3：在每个不同的分块中构建监控统计量，设定β＝95％以确定此时的控制限，分别监控不同分块中的监控统计量超出控制限的情况，认为分块中两个监控统计量均超出控制限时可能存在故障因素，将此时的监测变量传递给对应分块的支持向量机分类模型。
[0045] 步骤S4：设定滑动窗口长度N为10，阈值为90％，当连续10个采样点的监控统计量均超出控制限时使用支持向量机对上述对应的监测变量进行二次故障检测。当滑动窗口内检测出故障的采样点个数与滑动窗口长度的比值超出设定的故障阈值90％时则确定出现故障，计算此时的故障检测精度，可以发现获得的误判率为0，漏判率均为0.002.[0046] 步骤S5：分别计算并比较每一分块的故障贡献率，确定具有最大故障贡献率的分块即为故障发生分块。具体的如图5所示，可以发现对于空冷器冷却不足故障影响最大的是空冷器2分块。此外由于空冷器冷却不足造成进入气缸空气温度上升，同样会影响气缸工作，可见对于不同的气缸分块，故障的影响相似。随后进入涡轮器，影响开始变小，而对于增压器的影响最小。
[0047] 以上所述实施例仅表达本发明的实施方式，但并不能因此而理解为对本发明专利的范围的限制，应当指出，对于本领域的技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些均属于本发明的保护范围。
