0 引　言

近年来，我国经济的快速发展促进了船舶事业的繁荣，机舱作为船舶的动力场所，存放的各种油料和电器设施，成为火灾的高发区^{[1, 2]}。因此，研究船舶机舱火灾对于消防训练和火灾逃生有重要意义。

测量火灾实验过程中空间内的温度分布，对研究火灾行为规律有重要意义^{[3, 4]}。就固定空间而言，由美国NIST研究所开发的基于火灾动力学的大涡模拟软件FDS（Fire Dynamics Simulator）对火灾温度流动模拟的准确性已得到广泛认可^[5]。但若改变机舱火源功率、火源面积、补风口面积、风机流速等条件，FDS就不能满足要求，所以有必要建立火灾温度预测模型。20世纪90年代Vapnik等提出的支持向量机（Support Vector Machines，SVM），是一种新的机器学习算法，并迅速在多个领域得到快速发展^[6]。相比其他传统学习方法，具有良好的泛化能力，能够很好地克服维数灾难、小样本、局部极小点以及过拟合等问题^[7]。因此，运用支持向量机技术求解实际工程中有关火场温度预测的问题非常适宜。

本文以船舶机舱为研究对象，采用大涡模拟方法，以火灾动力学分析软件FDS为平台，计算不同条件下的船舶机舱火场温度变化规律。以计算结果为训练样本建立船舶机舱火灾温度的支持向量机预测模型，最后将遗传算法优化支持向量机（GA-SVM）理论运用于船舶机舱火灾温度的预测。为了验证本文模型的准确性，对比了SVM模型和BP神经网络方法的船舶机舱火灾温度预测结果。预测结果表明，利用遗传算法优化支持向量机的方法的预测效果较好且准确度较高，可以作为快速预测火场温度的工程计算方法。

1 船舶机舱火灾数值模拟

本文运用FDS软件对实习公司的船舶机舱实验室进行舱底油水火灾数值模拟。

1.1 船舶机舱火灾场景设计

以实习公司船舶火灾实验室为模拟对象，该实验室内部尺寸15 m（长）×7 m（宽）×5 m（高），机舱内部尺寸是6 m（长）×7 m（宽）×5 m（高），舱内有2台柴油机和1台柴油发电机，以及其他一些常用机舱设备。建立的模型如图1所示。采用机械通风和自然通风，计算场景中在四面墙壁设有多个格栅作为补风口，每个格栅是0.8 m（长）×0.6 m（宽），顶部设有机械抽风口，抽风口是0.5 m（长）×0.5 m（宽）。

本文将火源置于实验舱地面左边，燃料为柴油，火源为超快速t²火，最大火源功率为1.8MW。计算场景的环境温度为20℃。墙面、地面和天花板都为304钢板表面，墙面中间填充耐热材料。通过对该场景的网格独立性分析^[8]，本文使用了72×36×25非均一网格系统对所有涉及样本进行计算。

1.2 船舶机舱火灾样本计算

在整个船舶机舱火灾场景中，火场温度变化的过程是一个动态过程，受多种因素的影响，本文主要考虑火源功率、补风口面积、风机流量、火灾面积这4个因素对温度场的影响，构建了表1中15组火灾场景，每组只改变一个影响因素，采用FDS分别对这15个火灾场景下的船舶机舱火灾温度变化过程进行模拟，计算其不同时刻的温度分布。在主频为2.80 GHz的Intel CPU计算速度下，表1中每组火灾场景计算到400 s需要使用CPU的时间不少于4 h。15组数据计算结果如表2所示。

表1中各符号的含义如下：Q为火灾功率，kW；T为风机启动时间，s；V为风机流速，m/s；S为补风口面积，m²。

2 船舶机舱火灾温度预测模型的建立 2.1 基于SVM模型

为检验算法训练的精度，需要将第1节中15组FDS计算数据划分为13个训练样本和2个测试样本，将表1中的Q，T，V，S四个主要影响因素作为支持向量机模型的输入量 ${x_i} = {\left[ {Q,T,V,S} \right]_i}$ ，表2中FDS计算结果作为模型的输出 ${y_i} = {\left[ {{t_1},{t_2},{t_3}, \ldots {t_{21}}} \right]_i}$ ，因此15组训练样本集为 $\left\{ {\left( {{x_i},{y_i}} \right),i = 1,2, \cdots ,15} \right\}$ 。支持向量机将低维复杂的输入通过非线性映射 $\phi \left( x \right)$ ，投影至高维空间中，并在高维空间中利用结构风险最小化原则建立线性回归函数

${y'} = w \cdot \phi (x) + b\text{。}$

(1)

式中：w为权值向量；b为偏差。参数w和b的值由SVM采用最小化结构风险 ${R_{STR}}$ 来确定， ${y'}$ 为预测值，y为实际值。

2.2 遗传算法参数寻优

为了使SVM有最优的性能，本文采用遗传算法对SVM模型的惩罚因子C和核参数g进行参数寻优，避免认为选择参数的盲目性。遗传算法是一种借鉴生物界自然选择和进化机制发展起来的高度并行、自适应的全局优化概率搜索算法，已被广泛应用于组合优化、机器学习、信号处理等领域^[10]。

3 结果对比及分析 3.1 参数优化对比

本文采用网格搜索寻优和GA参数寻优对支持向量机的惩罚系数C和核参数g进行优化分析和选择，见表3。利用网格搜索寻优法时，先粗略寻找最佳参数，再进行精细的参数选择。

从表3可以看出，以上2种寻优方法对模型参数C和g寻优后，GA的适应度值0.046 691最小，模型的适应能力最好，对应的参数[C，g]=[81.192 6，0.300 69]组合最佳。由GA寻找最佳参数的适应度变化曲线图2 可知，进化代数为150次时终止进化，此时适应度值0.046 691就接近理想最优值。

3.2 预测结果分析

为比较GA-SVM模型的船舶机舱火灾温度预测的性能，分别采用SVM模型和前馈型BP神经网络建模对相同的样本集进行训练和预测。GA-SVM模型选用RBF核函数，经反复参数寻优，得到满足要求的C和g分别为 81.192 6和0.300 69。利用神经网络预测时，经过反复试验，网络的输入层、隐含层、输出层的节点数分别为 4，30，21，隐含层传递函数选用tansig函数，输出层传递函数选用purelin函数，训练函数选用Levenberg-Marquardt算法，学习速率Lr=0.2，目标误差Eg=0.000 01。3种模型下船舶火灾温度预测值与FDS计算结果比较见表4，可以得出，GA-SVM模型有较强的泛化能力，对训练样本学习后，对预测输入样本进行预测，得出相应预测输出值与FDS计算结果吻合很好。

图3列举了采用GA-SVM模型、SVM模型和BP神经网络方法进行火灾温度预测的效果图。

图4~图6分别列举了GA-SVM模型、SVM模型和BP神经网络预测值与FDS计算结果的相对误差曲线，由图4~图6可知，3种模型的预测值与FDS计算结果的最大相对误差分别为0.005 63，0.083 81和0.200 35。最小相对误差分别为0，0和0.008 63，GA-SVM模型和SVM模型均可以达到消防工程的精度要求，且GA-SVM模型预测精度明显高于SVM模型。

3种模型预测性能结果对比见表5，用MSE，e_MAPE和R三个指标进行预测性能评价。

从表5 可看出，GA-SVM模型的均方误差MSE和平均相对误差e_MAPE值分别为 0.046 691，0.001 6，远小于SVM模型和BP网络预测的相应指标值，且GA-SVM模型的相关系数R为 0.971 2，比SVM模型和BP网络模型的相关系数更接近1，其预测结果精度更高，适应性更好，表明本文模型预测性能更优。

4 结　语

本文利用FDS软件对不同条件下的船舶机舱火灾场景进行模拟，以火灾功率、风机启动时间、风机流速和补风口面积这4个影响因素为输入参数，不同时刻的火灾温度为输出构造了船舶机舱火灾温度的GA-SVM预测模型。在建模过程中，采用GA对SVM回归模型的惩罚因子C和核参数g自动寻优选择，避免了人为主观经验地确定参数。通过输入火灾功率、风机启动时间、风机流速和补风口面积到GA-SVM模型中，其对火灾温度预测的效果较好。与SVM模型和BP神经网络模型的相对误差对比分析表明，GA-SVM模型具有较好的泛化能力，为船舶机舱火灾温度间接计算及预测提出了一种新方法。

本文建立的船舶机舱火灾温度的支持向量机计算模型包括4种条件参数对火灾温度的影响。在今后的工作中，可以设计更多的火灾场景，例如考虑火源位置、补风口位置、火灾空间高度、面积等因素对船舶机舱火灾温度的影响，通过FDS数值模拟构造更加全面系统的训练样本，从而完善模型功能并进一步提高预测精度。