1、第 45 卷（2023）第 3 期Vol.45（2023）No.3柴油机Diesel EngineDOI：10.12374/j.issn.1001-4357.2023.03.003船用中速柴油机预测模型的自动标定与优化刘岱1，李玉超1，刘龙1，马修真1，桂勇2，夏倩2（1.哈尔滨工程大学动力与能源工程学院，哈尔滨 150001；2.中船动力研究院有限公司，上海 200120）摘要：为解决柴油机预测模型参数过多导致的标定困难问题，采用遗传基因算法对搭建的船用中速机DIjet预测燃烧模型进行智能标定。借助最大缸内压力以及对应的曲轴转角快速完成缸压曲线标定，相较于传统手动标定，标定时间可缩短到几十分
2、之一。在进行柴油机整体功率验证后，对NOx排放进行标定，形成一套简便可行的智能标定流程。在完成标定的模型上进行NOx排放及有效燃油消耗率（brake specific fuel consumption，BSFC）优化，在考虑两者不同权重系数的情况下得到更优的帕累托前沿面。关键词：中速柴油机；预测模型；遗传基因算法；自动标定；智能优化中图分类号：TK422文献标志码：A文章编号：1001-4357（2023）03-0014-06Automatic Calibration and Optimization of Predictive Model for Marine Medium Speed Di
3、esel EngineLIU Dai1，LI Yuchao1，LIU Long1，MA Xiuzhen1，GUI Yong2，XIA Qian2（1.College of Power and Energy Engineering，Harbin Engineering University，Harbin 150001，China；2.China Shipbuilding Industry Power Engineering Institute Co.，Ltd.，Shanghai 200120，China）Abstract：In order to solve the calibration dif
4、ficulties caused by excessive parameters in the diesel engine prediction model，a genetic algorithm was used to intelligently calibrate the combustion prediction model DIjet for a marine medium speed engine.By utilizing the maximum cylinder pressure and corresponding crankshaft angle，the cylinder pre
5、ssure curve calibration was quickly completed，which reduced the calibration time to the several tenth of its original time compared to traditional manual calibration.After verifying the overall power of the diesel engine，the NOx emissions was calibrated to form a simple and feasible intelligent cali
6、bration process.The NOx and brake specific fuel consumption（BSFC）was optimized on the calibrated model，and a better Pareto front surface considering the different weight coefficients of the two was obtained.Key words：medium speed diesel engine；prediction model；genetic algorithm；automatic calibration
7、；intelligent optimization0引言柴油机因输出转矩大、热效率高、经济性好、工作可靠等特点，广泛用于船舶动力系统。国际海 事 组 织（International Maritime Organization，IMO）于2016年制定严格降低NOx等污染物的Tier 排放指标1。Tier 标准将NOx排放指标限制在3.4 g/（kWh）以内，与原Tier 标准相比，船舶柴油机的NOx排放限值降低约79%。2020年9月，习近平主席在第七十五届联合国大会上郑重宣布：2022年内燃动力碳中和与排放控制学术年会绿色船舶动力系统分会场专栏收稿日期：2022-09-23基金项目：科技部国
8、家重点研发计划“船用柴油机可控燃烧技术研究”（2017YFE0116400）2023 年 5 月 15 刘岱等：船用中速柴油机预测模型的自动标定与优化中国将提高国家自主贡献力度，二氧化碳排放力争于2030年前达到峰值，争取2060年前实现碳中和。严格的排放法规与国家发展战略推动船舶柴油机向着高效率、低排放的方向发展2。在这种形势下，柴油机研发人员开展基于仿真模型的智能化优化研究3-5。对柴油机进行仿真时，使用较多的燃烧模型是Wiebe模型6，该模型能够较好地拟合燃烧曲线，但是需要收集较多的数据。船用柴油机的试验成本较高、试验数据较少，所以在对船用柴油机进行仿真时，多使用 DIjet 预测模型7
9、-9，通过现象学建模方法对燃烧的各个过程（喷雾、破碎、混合、着火滞燃期、燃烧速度、排放计算等）建模与拟合，能在数据较少的情况下得到更好的效果。该模型须对较多参数进行标定，即使是经验丰富的专家也须花费大量的时间。本文使用一种遗传算法对模型进行燃烧与NOx排放标定，然后在标定好的模型上借助遗传算法对有效燃油消耗率和 NOx排放进行优化，形成简便、快捷的智能标定与优化流程。通常采用手动法仅对燃烧过程进行标定就须耗时1个月以上，而采用本文的方法，燃烧与NOx排放的整体标定时间可缩短至2天。因此，该方法具有较大的优势。1模型介绍本文的研究对象是320中速柴油机，表1为发动机的主要参数。在 GT-POWE
10、R 软件中建立发动机模型，考虑到标定与优化须进行大量的重复运算，为减少模型的复杂程度，缩短运算时间，在标定时仅建立单气缸模型，见图1。该模型采用软件自带的DIjet燃烧模型对柴油燃烧过程进行仿真，通过输入喷油压力曲线计算喷射速率、平均液滴直径等，模拟液滴蒸发、着火延迟等过程，从而得到燃烧放热曲线，用于预测缸压与NOx排放。具体过程如下。（1）贯穿距贯穿距是喷油器的一个十分重要的参数，该参数对喷入气缸中燃料的分布区域有很大的影响。通常为了避免类似“火包油”现象的发生，对该参数有一定的要求：Sp=Uitbr(tbrt)0.5-0.06（1）式中：Sp为贯穿距；tbr为破碎时间；Ui为初始燃油喷射速
11、度；t为时间。涡流会影响贯穿距离，影响关系为：Ss=S()1+RsnSUi（2）式中：n为发动机转速；S为涡流前贯穿距离；Ss为实际贯穿距离；Rs为涡流比。（2）破碎在破碎之前，燃料以喷射柱的形式存在，随后燃料分解成许多小液滴。燃油的破碎/雾化对燃烧放热规律和着火点位置等参数均有重要影响。tbr=0.29(3+LnDn)(fuelair)(DnUi)Kbr（3）式中：tbr为燃料的破碎时间；Dn为喷油器的喷孔直径；Ln为喷油器的喷孔长度；Kbr为破碎时间乘数；x为流体密度。（3）空气卷吸随着射流的运动，气缸中的每一个区域都会卷吸气缸新吸入的空气以及上个循环未排除的残余气体。ma=mf UidS
12、sdt-1 f(EGR，qm)（4）式中：ma和mf分别为卷吸量和区域中燃油的质量；dSsdt为贯穿速率；EGR为燃烧废气的体积分数；qm为喷油的质量流量。每个子区域的空气卷吸率为：dma，xdt=Cxdmadt（5）根据情况的不同，空气卷吸率dma，xdt 也会随之变化。在子区域的点火点之前，Cx表示燃烧前的表1发动机主要参数项目缸径/mm冲程/mm连杆长度/mm排量/L压缩比进气阀直径/mm排气阀直径/mm参数32042095033.81393.389.1图1GT-POWER发动机模型第 45 卷第 3 期 16 柴 油 机卷吸率；在子区域开始燃烧后，Cx表示开始燃烧后的卷吸率；在燃油喷射
13、射流撞击气缸顶部后；Cx表示撞击后的卷吸率。（4）点火延迟点火延迟对整个燃烧过程乃至柴油机的功率和指示燃油消耗率都有着十分重要的影响，该值的大小会影响可燃混合气的形成数量和质量。每个未燃烧子区域的点火延迟是单独计算的。I=1dt（6）=C1()3-2pC2eC4T1.210-6f()f，C8，3.0，3.0（7）式中：C1为点火延迟乘数；C2为压力延迟系数；C4为活化温度延迟；C8为稀释效应乘数；f为所有不是O2、N2或燃料工质的质量分数；I为点火积分；为等效比；为点火时间常数；p为压力。（5）燃烧燃烧是影响柴油机动力性、经济性和排放特性的关键因素。该燃烧模型可以处理混合燃烧以及有动力学限制的
14、燃烧，该模型燃用未燃烧放热分区中的变成蒸气的燃料或动力学限制量中的较小量的燃料。在任何情况下，燃烧区的燃料和适量的空气都是从未燃烧区转移过去的。在燃烧结束后，由于子区域的物质被氧化、区域中卷吸的空气被稀释，燃料的能量可能会继续释放。dmkdt=()3-2p2.5e-4 000T 3.0 （8）式中：dmkdt为有限燃烧动力学燃烧速率；为燃烧率乘数；p为压力；T为子区域的温度。2方法介绍使用遗传算法对发动机模型进行自动标定和优化，标定的数据来源于试验得到的缸压等性能与排放数据。借助GT-POWER软件中自带的优化工具能够方便地使用遗传算法以及其他的一些优化方法进行优化。2.1遗传算法遗传算法10
15、是模仿自然界生物进化机制发展起来的随机全局搜索和优化方法，借鉴达尔文的进化论和孟德尔的遗传学说，其本质是一种高效、并行、全局搜索的方法，能在搜索过程中自动获取和积累有关搜索空间的知识。图2是遗传算法流程图。在遗传算法运算的过程中，首先生成初始化种群，种群的规模由前面设置的N来确定。生成种群后，将对所有个体进行适应度计算，然后依概率进行选择，留下适应度较高的个体进行交叉、变异，产生下一代种群，同时保留适应度最高的个体，之后重新对新产生的种群进行适应度的计算与比较，直到算法满足终止条件。对于发动机的标定，适应度函数是模型运算值与试验测量值的差值。在使用遗传算法进行求解时，根据所需求解自变量数目的不
16、同，推荐设置的种群规模也不同，表2是根据自变量数目推荐的种群规模。2.2适应度函数直接借助 GT-POWER 软件中的遗传算法模块将模型输出的参数与试验得到的参数组合成适应度函数，并进行迭代计算。在对缸压曲线标定时使用的适应度函数如下：F(x)=|Pmaxm-Pmaxe|+|Cm-Ce|（9）式中：F（x）为适应度；Pmaxm为模型运算得到的缸内最大压力；Pmaxe为试验数据的缸内最大压力，Cm为最大压力对应的曲轴转角；Ce为试验得到的最大压力对应的曲轴转角。在对缸压曲线进行标定时，准确的边界条件通常可以保证压缩曲线的准确性。由于对燃烧的情况进行标定分析比较困难，使用最大缸压以及最大缸压对应的曲轴转角作为评价燃烧曲线的标准进行自动标定。这两个图2遗传算法流程示意图表2不同变量的推荐种群规模项目变量数种群规模参数3104165206267308409502023 年 5 月 17 刘岱等：船用中速柴油机预测模型的自动标定与优化参数在实际的缸压曲线以及模型运算的缸压曲线中均较易获取。完成缸压曲线的标定后，对 NOx排放进行标定，使用的适应度函数为：F(x)=|GNOxe-GNOxm|（10
