Chinese Journal of Ship Research
Evaluation of main engine energy efficiency based on BP neural network
陈伟南，黄连忠，张勇，路通 116026大连海事大学 轮机工程学院，辽宁 大连
摘 要：［目的］在船舶航行期间，需要通过分析船舶和主机的运行参数来客观判断主机当前的工作情况，从而准确评估主机的能效状态。［方法］以状态良好的船舶运行记录为样本，结合主成分分析法和BP神经网络算法，构建船舶的航行状态识别模型和主机油耗模型，并在船舶航行期间对船舶实时运行参数进行分析，得出船舶主机在当前工况下的油耗量正常值。以某30万吨级远洋散货船为例开展模型计算验证，将正常油耗值与实际油耗值进行对比，以二者的残差值为依据，进而评估当前的主机能效状态。［结果］计算结果显示，航行状态识别模98.05%，油耗模型的平均误差为3.47%，2型的正确率为 种模型的可靠性均较高。［结论］研究成果可为智能船舶的能效管理提供一定的参考。关键词：BP神经网络；主成分分析法；油耗模型；船舶能效管理中图分类号：U676.4+1 文献标志码：A DOI：10.19693/j.issn.1673-3185. 01219
Evaluation of main engine energy efficiency based on BP neural network CHEN Weinan，HUANG Lianzhong，ZHANG Yong，LU Tong College of Marine Engineering，Dalian Maritime University，Dalian 116026，China
Abstract：［Objectives］During the voyage of a ship，the operating parameters of the ship and its main engine need to be analyzed so as to objectively judge the current working state of the main engine and accurately evaluate its energy efficiency.［Methods］Taking good ship operation records as samples，and combined with principal component analysis and the BP neural network algorithm，a ship navigation state identification model and main engine fuel consumption model are built. During the voyage of a ship，these two models are used to analyze the real-time operating parameters of the ship in order to obtain the normal value of the main engine fuel consumption under current working conditions. The model of a 300 000-ton ocean bulk carrier is used to calculate and verify the models，the normal value is compared with the actual fuel consumption value of the main engine，which is based on the residual value of the two，and the current energy efficiency state of the diesel engine is then evaluated.［Results］The validation results show that the correct rate of the navigation state identification model is 98.05% and the average error margin of the fuel consumption model is 3.47% ，thus proving the two models to be more reliable.［Conclusions］ The results of this research can provide references for the energy efficiency management of intelligent ships. Key words：back-propagation neural network；principal component analysis；oil consumption model； ship energy efficiency management
0引言
对船舶主机而言，其能效状态由主机的设定转速、外界的海洋气象状况、船舶载重、船舶姿态等多个因素共同决定。在传统的船舶管理模式中，船舶操纵人员大多依靠自身经验对主机状态进行评估，并依据评估结果进行操作。这种评估具有极大的不确定性，并降低了船舶的能源利用效率［1］。刘伯运等［2］和孙宜权等［3］针对车用柴油机的状态评估进行了分析。在船用柴油机方面，周根明等［4］运用层次分析法对船舶主机的整体健康状态建立了评估体系；孙峰等［5］基于数据挖掘技术提出了一种评估船用柴油机性能的方法。这些方法虽然对船舶主机的能效状态评估进行了研究，但并不能对航行中的船舶主机状态进行实时、迅等［6］利速的评估。叶睿 用神经网络构建了客滚船的主机油耗模型，但其输入参数较为复杂，没有进行简化处理，同时其输入参数仅为船舶航行时的船体姿态数据，并未包含主机的运行参数，这也影响了模型的可靠度与准确性。基于此，本文将以状态良好的船舶运行记录作为样本，结合主成分分析法和反向传播（Back-Propagation，BP ）神经网络智能算法，构建船舶航行状态识别模型和主机油耗模型。在构建油耗模型时，将利用主成分分析法对模型输入参数进行简化，从而降低数据噪声和网络复杂性。此外，本文油耗模型的输入参数还将包含主机的运行参数和船舶的航行姿态参数，用以更准确地反映主机的能效状态。本文构建的航行状态识别模型可用于识别船舶当前的航行状态，而油耗模型可用于判断船舶主机当前的能耗状态是否正常，从而对主机能效进行智能、准确的评估，可为智能船舶的能效评估研究提供一定的参考。
1 构建模型
采用船舶的航行状态模型和主机油耗模型对1主机状态进行评估的流程如图 所示。首先，采集船舶的实时运行数据，将其输入到基于神经网络构建的主机航行状态识别模型中，运用该模型识别出主机不同运行数据对应的航行状态；然后，选取正常航行状态数据，运用主成分分析法对数据进行预处理，用以降低后续计算的复杂度，加快神经网络的收敛速度并减小误差；接着，将处理后的航行状态数据输入到基于神经网络构建的主机油耗模型中，得出主机的功率、转速、油耗等参数； 最后，将主机的实际油耗值 x2 与理论油耗值 x1 进行实时对比，得到油耗残差值 δ ，即可判断主机当前的运行状态是否正常。
1.1 目标船及数据
本文的船舶实时运行数据均来源于某目标船——30 327 m，万吨级远洋矿砂运输船。该船总长55 m 29 m 21 m型宽 ，型深 ，设计吃水 ，航速14.5 kn 298 000 t MAN B&W ，载重 ，主机型号为6S80MC-C 22 360 kW ，额定功率为 ，额定转速为73 r/min。通过船载传感器进行数据采集工作，目标船1的状态参数如表 所示。其中，风速分解为沿船舶航行方向（X方向）和垂直于航行方向（Y方向） 2个风速，并分别以船艏来风和右舷来风为正。 30s传感器每隔 采集一次数据，为保证模型5 min的稳定性，以 时长为间隔来处理所采集的数5 min据。其中，除主机油耗量数据为 内的累计值5 min之外，其他数据均为 内的均值。X方向和Y方向的风速分别为N X ˉ = 1 åV cos θc w c N （1） c =1 N Yˉ= 1 åV sin θc w c N c =1
5 min式中：N为 内传感器采集的总次数；V 为风w c c=1，2，…，N；速仪第c次采集的相对风速，其中θc为船舶航向与风向之间的夹角。主机缸套水的进口、出口温差平均值ΔT为6 ΔT å 2 = 1 (T - T in) （ ） outg 6 g =1 1～6式中：T 为主机缸套水的进口温度；T 为in outg =1，2，…，6。号气缸缸套水的出口温度，其中 g 1～6号气缸排气温度的平均值Tˉ 为6 Tˉ = 1 åTg 3 （ ） 6 g =1
1～6式中，Tg 为 号气缸的排气温度。
1.2 BP神经网络
对于输入参数繁多的船舶航行状态识别模型和油耗模型而言，若构建一般的函数，则难以准确
描述各物理量之间的因果关系，且普通函数对各6］。 BP型船舶的通用性较差［ 本文将采用 神经网BP络来构建模型。 神经网络是一种非线性智能算法，具有运算速度快、容错能力高、自学能力强
等特点，是目前应用最广泛的神经网络模型之BP一。 神经网络可以学习和存储大量的输入—输出模式映射关系，而无需事前揭示描述这种映射关系的数学方程［7-9］。BP神经网络是一种多层前馈型神经网络，其2拓扑结构如图 所示。图中：P1 ， P ，…，Pn 为2 BP神经网络的输入参数；b1 ，b ，…，bm 和 θ ， 2 1 θ ，…，θk 为网络的偏置；a1 ，a ，…，ak 为网络2 2 BP的预测值； ωij 和 ω 为 神经网络的权值。其jh i=1，2，…，n；j=1，2，…，m；h=1，2，…，k。中 2图 中的网络包括n个输入节点和k个输出节点，表达了从n个自变量到k个因变量的函数映射关系。网络中隐含层第j个节点的输出结果 yj和输出层第h个节点的输出结果 ah 分别为 n 4 yj = f1 åωij Pi + bj （ ） i =1
k åω 5
ah = f2 y +θ （ ） jh j h h =1
式中，f1和 f2分别为输入关系和输出关系的传递函数。
Matlab本文将采用 软件来构建神经网络，选tansig，f2取 f1为双曲线正切函数 为线性变换函purelin。BP 2数 神经网络的学习分为 个阶段：第1 1个阶段是根据输入的已知样本，从第 层向后计2算各神经元的输出；第 个阶段是从最后一层向前计算各权值和阈值对总误差的影响，进而修改2各权值和阈值。这 个阶段反复交替进行，直至函数收敛。为了提高网络的收敛速度，本文选用特（Levenberg-Marquardt，LM）列文伯格—马夸尔2算法作为第 个阶段的计算方法，并以均方误差函数作为性能函数。
1.3 主成分分析法
5航行状态识别模型的输出结果分为 种，即停泊、机动航行、正常航行、冲车航行和大风浪航行，其设计输入为相对风速 Xˉ 和 Yˉ 、船舶航速Vs 、主机转速 ne 、主机曲轴扭矩M、主机曲轴功率BP Pe 及主机燃油消耗量F，运用 神经网络算法可以实现较高准确度的航行状态识别。
然而，对于主机油耗模型而言，其输出结果是10 1中具体的燃油消耗量，其输入参数多达项（表10前项）。在此情况下，变量之间有一定的关联BP性，故变量提供给神经网络的信息会存在重叠，从而增加了计算的复杂性，减缓了网络的收敛BP速度，降低了模型的准确性并提高了 神经网络性［10-11］。如果神经网络过于复杂，过拟合的可能会使信号和噪声同时进行拟合，从而导致神经网
络的预测值与实际值严重不符，而采用大量的数
据对网络进行训练并降低训练数据的噪声则是避免出现过拟合问题的有效方法［12］。因此，本文除了采用大量数据对模型进行训练之外，还将应用主成分分析法（Principal Component Analysis）对油耗模型的输入数据进行预处理以降低数据噪声。
主成分分析法的实质是在不改变样本数据结构的情况下，通过旋转多维空间的坐标将原变量
转换成两两不相关的主成分，同时尽可能反映原变量所包含的信息，从而简化计算［13-14］。本文将BP采用主成分分析法与 神经网络算法相结合的方法，假设样本矩阵 X 为
式中：样本矩阵 X 中第 v 个参数类型中的第 u 个=1，2，…，n =1，样本数据为 xuv ，其中 u 且 v 2，…，p；x =［xu1 xu2 xup］ T ，为p维列向量。u样本矩阵的均值 xˉ 为å n 7 xˉ = 1 x = (xˉ1 xˉ2 xˉ p)T （ ） n u u =1 n å
式中：x ˉ1 x ˉ2 xp ˉ 中的任一元素 xˉ = 1 xuv 。v n u =1协方差矩阵 S 的表达式为n 1 å( )( )T （8）
S= xu - xˉ x - xˉ n - 1 u u =1对于n行 p列的样本矩阵 X 而言，经过矩阵S ，其为pp变换后即得到样本协方差矩阵 行p列的8 =1，2，…， =1，方阵，故式（ ）之后的 u 且 v 2，…，p。协方差矩阵 S 中的任一元素 Suv 为å( n 1 （9） Suv = xu′ - xˉu )( xu′ - xˉ n - 1 u v v) u′ =1式中： xu′ 和xu′ 分别为样本矩阵 X 中的某个元u v p =1，2，…，p；xˉu å素，其中 u′ = 1 xuv 。p v =1样本相关矩阵R为=( ruv )p 10 R （ ） ´p
其中
2 航行状态识别模型
对正常航行的船舶而言，在实时评估其主机能耗状态之前，应先识别当前船舶的航行状态。在停泊状态下，船舶航速、主机转速、耗油量基本0；在机动航行状态下，各项参数基本在停泊状为态和正常航行状态的数值特征范围内波动；在正常航行状态下，各项参数基本稳定在一定范围内；在冲车航行状态下，各项参数将出现波动且有所 增加；在大风浪航行状态下，相对风速将明显增加，而船舶航速将有所降低。根据不同航行状态的参数特征，将现有的历史数据进行统计分类，即可得到各个航行状态的样本数据。为了实现较高准确度的航行状态识别，本文将以船舶航速、相对风速、主机转速、曲轴扭矩、曲轴功率和主机油耗量作为输入，船舶航行状态作为输出。同时，还需要设定神经网络计算矩阵的10000 010目标值，即停泊为［ ］、机动航行为［ T 00 00100 000 ］、正常航行为［ T ］、冲车航行为［ T 10 00001 ］、大风浪航行为［ T ］，然后进行训练学T Matlab BP习。在 软件中建立并训练 神经网络时， trainlm、学习函数为 learngdm、性设定训练函数为mse、隐含层和输出层的传递函数分别能函数为tansig purelin，其中 2为 和 层隐含层的神经元数量10。训练函数trainlm LM均为 采用 算法，该方法结合了梯度下降法和牛顿法的优点，可以有效避免陷入局部最优解并保证收敛速度［15］。本文将选取目标船在新加坡—圣路易斯航段（2015 03 20 2015 04 06年 月 日至 年 月 日）作为训练样本，即该时间段内船舶和主机均处于良好状22态下的航行数据，其中停泊记录 条、机动航行33 153记录 条、正常航行记录 条、冲车航行记录14 114条、大风浪航行记录 条。4为了验证模型的准确性，本文选择了 个月后的航行记录模型作为检验样本，即目标船在马段（2015 07 19 2015迹山—新加坡航 年 月 日至 年07 22 77h月 日）内连续 的样本数据，仿真结果如3图 所示。图中，各分图的横坐标均为传感器数据采集点的序号。3（g）所示的评估航行状态图即模型识别所图 3（h）所示的实际航行状态得的船舶航行状态，图图即通过查询该目标船轮机日志和航海日志后统1，2，计得出的船舶实际航行状态，其纵坐标数值3，4，5分别代表停泊、机动航行、正常航行、冲车5 3航行、大风浪航行这 个状态。由图 可知，在第0～120采集点之间，船舶航速、主机转速、扭矩、0，识别模型判断为停泊状态，功率、油耗等均为120～135判断准确；在第 采集点之间，船舶和主
机的各项参数均出现了波动，模型将其识别为机128动航行状态，判断准确；在第 点左右，由于船舶的短暂提速，模型将其识别为正常航行状态，判135～486断错误；在第 点之间，模型判断船舶为正常航行状态，经查询航海日志，该目标船确实为486～719正常航行状态，判断准确；在第 点之间，风速提高、船速降低，模型判断为大风浪航行状885～892态，经查询航海日志，判断准确；在第 点之间，各项参数均出现了波动且有所增加，模型判断为冲车航行状态，经查询轮机日志，该目标船在此时间段确实进行了冲车操作，故模型判断准确。由此可知，模型判断整个样本的识别正确率98.05%，从而为下一步建立主机油耗模型提供为了基础。
3 主机油耗模型 3.1 主成分分析处理
BP应用 神经网络构建主机油耗模型时，需要先对训练样本进行主成分分析。本文选取了目标2015 3～5船在 年 月的船舶和主机运行数据作为训练样本，该时间段内目标船和主机均处于良好1 10运行状态。选取表 的前 个参数作为主成分分析的原始变量，设定主成分Z1～Z10，计算所得的成2分矩阵如表 所示，各主成分对应的特征值、方差3贡献率及累计方差贡献率如表 所示。3可知，Z1～Z5 93.87%，由表 的累计贡献率为而 Z 6 ～Z10的贡献率则很小。因此，本文将选择Z1～Z5作为表述样本内容的主要数据来计算各项4主成分的系数，结果如表 所示。根据各主成分10 5的系数，将原训练样本的 种原始变量转换为BP项主成分，即可作为 神经网络的输入数据。
3.2 建立油耗模型
matlab BP在 中设定 网络油耗模型的训练函trainlm learngdm数为 、学习函数为 、性能函数为mse、隐 tansig含层和输出层的传递函数分别为 和purelin。对于隐含层的层数和神经元节点数目，目前缺乏设定最优数目的明确指导方法，本文将0～20，然设定神经元节点数目的取值范围为 后2采用循环程序遍历 层隐含层中所有神经元数目组合的模型训练误差，最后根据计算所得的训练4所误差值确定最优的神经元数目设定。如图3.1 BP示，将 节中已处理的主成分数据信息作为网络的设计输入，通过循环程序即可得到不同神经元数目组合的训练误差结果。
1 16经计算，当第 层隐含层神经元数目为 且2 11第 层为 时，数据整体的均方误差最小，其值为6.325 9×10-5。因此，本文设定BP 2层神经网络中16 11，然后，即可隐含层的神经元数目分别为 和建立并训练网络，最终得出正常航行状态下的主 机油耗模型。需要指出的是，由于不同的船舶和主机的匹配关系与性能均有所差别，所以针对其他船舶构建航行状态识别模型和主机油耗模型时，还需要重新通过计算机程序进行神经网络训练，所需的时间成本也应在可接受范围内。本文已经确定了神经网络训练的模式和函数，故可以推广应用至其他船舶和主机。
3.3 模型验证与分析
选取目标船在圣路易斯—马迹山—新加坡航段（2015 05 19 2015 05 24年 月 日至 年 月 日）的部分航程作为验证样本，经航行状态识别模型提取正常航行状态数据，并根据各主成分系数计算得5 BP出 项主成分参数，输入 神经网络即可得出对5应的理论油耗值，结果如图 所示。
由于主机燃油流量计数器的输出数据精度为0.01 m3，故图 5中的实际油耗值曲线存在幅值为0.01左右的波动，但总体趋势与计算油耗值曲150线一致。在第 个数据点附近，计算油耗值出现了小幅降低，船舶的曲轴扭矩和功率分别1 620 kN·m 9 600 kW在 和 左右波动，略低于整体1 673 kN·m 10 022 kW，而实际油耗平均水平 和值也出现了降低现象，可见计算结果与实际运行500情况相符。在第 个数据点附近，实际油耗值出现了异常高峰点，经查为传感器的暂时故障所5可致，而其他数据点的传感器工作正常。由图3.47%，知，主机油耗模型的总体平均相对误差为且未产生明显的过拟合现象。将实际油耗值和计算油耗值进行对比，以残6差形式体现主机的能效状态，如图 所示。 6由图 可知，在该段航程内的各项工况下，主机的性能特征与良好能效状态相吻合，即主机能效状态正常。同时，由于主机燃油流量计数器的0.01 m3，故大部分样本数据点的输出数据精度为耗油量残差都在-0.01～0.01范围内波动。
4结语
BP本文结合实船航行数据，运用 神经网络算法建立了船舶航行状态识别模型，可以有效识别5 98.05%。同时，船舶的 种航行状态，其准确率达10依据船舶和主机的 个运行参数与主机燃油消耗量之间的关系，采用主成分分析法降低了各个BP变量之间的关联性和复杂性，并运用 神经网络算法构建了船舶主机正常值油耗模型，其总体平3.47%。通过综合运用这2均相对误差为 种模型，可以对正常航行状态下的主机实时状态进行评估，可为智能船舶的能效管理提供参考。
参考文献：
［1］ 董建华，吴雨华，郑士君.船舶状态评估及预算系统［J］. 中国航海，2009，32（2）：26-29. DONG J H， WU Y H ，ZHENG S J. A system for ship condition assessment and budgeting［J］. Navigation of China，2009，32（2）：26-29（in Chinese）. 2 刘伯运，欧阳光耀，常汉宝. D-S ［］基于神经网络和证据理论的柴油机状态评估［J］. 车用发动机，2005 （5）：14-17. LIU B Y ，OUYANG G Y，CHANG H B. Diesel engine condition assessment based on neural network and D-S evidence theory ［J］. Vehicle Engine， 2005 （5 ）： 14-17（in Chinese）.
［ 3］ 孙宜权，张英堂，李志宁，等. 基于核主成分分析的
柴油机技术状态评估［ J］. 车用发动机， 2012（2）： 89-92. SUN Y Q ，ZHANG Y T， LI Z N ，et al. Evaluation of diesel engine technical state based on KPCA［J］. Vehi⁃ cle Engine，2012（2）：89-92（in Chinese）. ［ 4］ 周根明，唐晓霞，郭霆，等. 基于AHP方法的船舶柴油机健康状态评估［ J］. 柴油机， 2013，35（6）： 9-13. ZHOU G M，TANG X X，GUO T，et al. Evaluation of the health state of marine diesel engine based on AHP method［J］. Diesel Engine，2013，35（6）：9-13（in Chinese）. ［ 5］ 孙峰，黄连忠，刘伊凡，等. 一种应用数据挖掘技术评估柴油机性能的方法［ J］. 大连海事大学学报， 2017，43（3）：83-88. SUN F，HUANG L Z， LIU Y F ，et al. A method of evaluating diesel engine performance by using data mining technology［J］. Journal of Dalian Maritime Uni⁃ versity，2017，43（3）：83-88（in Chinese）. 6］ 叶睿，许劲松. ［ 基于人工神经网络的船舶油耗模型［J］. 船舶工程，2016，38（3）：85-88. YER ，XU J S. Vessel fuel consumption model based on neural network［J］. Ship Engineering，2016，38 （3）：85-88（in Chinese）. 7］温正. MATLAB 法［M］. ［ 精通 智能算 北京：清华大学出版社，2015. WEN Z. Proficient in MATLAB intelligent algorithm ［M］. Beijing：Tsinghua University Press，2015 （in Chinese）. 8］ 闻新，李新，张兴旺，等. MATLAB ［ 应用 实现神经网络［M］. 北京：国防工业出版社，2015. WEN X， LI X ，ZHANG X W，et al. Application of MATLAB neural network［M］. Beijing：National De⁃ fense Industry Press，2015（in Chinese）. 9］ 董建华，郑士君，王伟彬，等. ［ 船舶状态检测技术与
评估方法探讨［J］. 机电设备，2004，21（3）：1-4. DONG J H，ZHENG S J，WANG W B，et al. Discus⁃ sion on the test technology and assessment method of ship status［J］. Mechanical and Electrical Equipment，
160 （下转第 页）