马丽娜，田 勇，王 倩，徐 灿

（南京航空航天大学民航学院，江苏 南京 211106）

近年来，全球航空运输业迅猛发展，航班数量显著增加，然而，航空器运行带来的环境问题却日益凸显。据联合国政府间气候变化专门委员会（Intergovernmental Panel on Climate Change, IPCC）统计，目前全球人为排放二氧化碳（CO2）总量的2%来源于航空运输排放，且该比例保持较快增速[1]。2019年9月，中国向国际民航组织（International Civil Aviation Organization, ICAO）正式提交第三版航空减排行动计划书，宣示了在绿色民航发展方面的决心与信心。航空器巡航阶段时间占比高，气体排放产生的温室效应具有放大效果[2]，对此，有必要研究一种有效的航空器巡航温室效应评估方法。温室效应表征和燃油消耗计算是温室效应评估方法研究的主要议题。

国内外学者针对航空器运行温室效应的研究多集中在气体排放计算方面。曹惠玲等[3]使用P3-T3方法对ICAO基准排放模型中的气体排放指数进行修正，在此基础上实现了标准起飞着陆循环（land and take off, LTO）过程氮氧化物（NOx）排放量的获取，随后对黑碳排放量进行了计算和分析[4]；OWEN et al[5]使用飞机性能和污染物排放模型对未来航空运输中CO2和NOx排放情况做出预测；罗雯莉等[6]使用线源扩散模型，比较了不同减排措施下污染物排放减少作用。针对航空器运行对机场周边区域造成的环境影响，学者们也开展了诸多研究，高垒等[7]对南京禄口国际机场2017年NOx、碳氢化合物（HC）和一氧化碳（CO）排放量进行了有效计算；韩博等[8]从系统水平出发，获取了京津冀机场群2018～2019年大气污染物排放清单。上述研究实现了对航空器运行污染物排放量的计算，但对时间占比较长的巡航阶段关注较少，且尚未针对污染物排放产生的综合温变影响进行进一步研究。

航空排放中的温室气体主要来源于燃油消耗，在燃油消耗获取方面，学者们大多通过数学建模的方法展开研究。TURGUT et al[9]针对航空器下降过程构建了燃油消耗与高度之间的指数模型；LAWRANCE et al[10]基于飞行数据构建了新型的动力学参数模型、高斯回归模型及K近邻回归模型；张军等[11]使用最小二乘支持向量机（support vector machines,SVM）的方法，提出了基于横向与纵向二维驱动的动态预测方法。随着研究的深入，气象因素的影响也引起了学者的关注，并由此形成了飞机油耗修正模型[12]。但，这些油耗估算方法多依赖于大量的飞行参数，其中绝大多数数据公开程度较低、获取难度较大，使得燃油消耗计算成为航空器运行环境影响和节能减排效果评估过程中的难点。

本研究将使用机器学习的方法，针对航空器巡航温室效应评估方法展开研究：根据温室气体排放和表征模型，建立航空器巡航综合温室效应评估模型；构建用于获取巡航燃油流量的三层神经网络（neural network, NN）模型，引入遗传算法优化网络参数；从预测精度和优化质量2方面评价改进后的网络模型，并选取深圳航空公司ZH9164航班进行实例应用，验证评估方法的可行性和可靠性。

1 综合温室效应评估模型

CO2和NOx是航空器巡航阶段气体排放的重要组成，两者具有不同的增温效应。因此，本文首先从温变影响的角度构建综合温室效应评估模型，以对航空排放中不同污染物的环境影响进行统一度量。

1.1 温室气体排放模型

1.1.1 CO2排放模型 根据美国联邦航空局（Federal Aviation Administration, FAA）气体排放扩散模型（emissions and dispersion modeling system, EDMS），CO2排放量与污染物排放指数和燃油消耗量相关，见式（1）。

式（1）中：ECO2为CO2排放量，kg；Tc为巡航总时长，h；FCR(t)为t时刻的燃油流量，kg/h；EICO2为CO2排放指数，kg/kg，表示消耗 1 kg航空燃油产生的CO2量，在发动机类型与航空燃油种类不变的条件下，EICO2保持恒定。表1可知，B737-800飞机选装的CFM56-7B26发动机CO2排放指数为3.155。

表1 CFM56-7B26发动机基准排放数据

1.1.2 NOx排放模型 对于 NOx排放模型，与 CO2排放模型不同，其排放指数随巡航条件的不同而实时变化。此外，基准排放模型中的NOx排放数据，只包括4种典型工作状态下的排放指数，且仅适用于基准条件（ISA、0 m），因此，针对航空器巡航NOx排放量的计算，需要根据具体燃油流量和巡航条件，对其排放指数进行插值和修正，见图1。

图1 NOx排放指数修正方法

①将实际燃油流量按式（2）修正到基准条件下的相应值。

式（2）中，FCR’为修正后的基准燃油流量；δ为外界大气压强与标准海平面大气压强（1 013.2 hPa）的比值；θ为外界大气温度与标准海平面大气温度（288 K）的比值；e为自然常数；M为巡航马赫数。

②根据发动机基本排放数据，对基准燃油流量和NOx排放指数进行拟合，以便能够通过拟合关系，获取基准燃油流量对应的NOx基准排放指数表示在基准条件下消耗1 kg航空燃油产生的NOx量。

③将NOx基准排放指数按照公式（3）转换回实际排放指数。

式（3）中：EINOx为NOx的实际排放指数，kg/kg，表示在实际运行条件下消耗1 kg航空燃油产生的NOx量；φ为外界大气相对湿度，%；pv为饱和蒸汽压，Pa，可采用公式（4）根据外界大气温度计算获得。

式（4）中：T为外界大气温度，K。综合EDMS排放模型和NOx排放指数修正模型，可得出NOx排放量表达式（5）。

式（5）中：ENOx为NOx排放量，kg；EINOx(t)为t时刻的NOx排放指数。

1.2 温室效应表征模型

在温室气体排放定量计算的基础上，使用全球绝对温变潜势（absolute global temperature potential,AGTP）这一表征参数进一步探究气体排放造成的大气温变影响。根据推导和验证，见式（6）。

式（6）中， A GTPG(H)为温室气体G（CO2或NOx）在决策时间水平H下的全球绝对温变潜势，K/kg，是其造成全球平均地表温度的变化；AG为单位质量G的 辐射强迫,Wm2/kg；a0、ai、αi、cj和dj为给定参数[13]。

由此，航空器巡航阶段CO2和NOx排放造成的综合温变影响，见式（7）。

式（7）中， ∆T(H)为全球总温变潜势，K。

2 巡航燃油流量神经网络模型

获取各时刻的燃油流量序列是准确评估巡航温室效应的关键步骤。本研究打破基于发动机性能模型对燃油流量进行建模的一贯研究，利用神经网络对燃油流量进行建模，以期减少计算过程中对公开程度较低的飞行参数的依赖。

2.1 基于神经网络的模型建立

人工神经网络是由大量称为神经元的节点和权重边组成的信息或信号处理系统，可通过模拟生物大脑的结构和功能，共同完成分布式处理和运算任务，具有学习能力强、可操作程度高的优点。本研究根据影响飞机巡航燃油流量的性能和气象因素，构建基于三层神经网络的巡航燃油流量估算模型FCR=N(h、v、d、ws、wd)，网络的拓扑结构，见图2。

图2 三层神经网络拓扑图

网络由输入层X、隐含层Y和输出层Z组成。其中：X包含5个神经元h、v、d、ws、wd，分别表示气压高度、地速、航向、风速和风向；Y层神经元用以传递X与Y间的关系；Z含1个神经元FCR，为燃油流量；wxy、b1和f1分别表示X到Y的连接权值、阈值和激活函数；wyz、b2和f2分别表示Y到Z的连接权值、阈值和激活函数。

2.2 基于遗传优化的算法设计

人工神经网络具有较强泛化能力和可操作性，但也存在明显缺陷。在网络训练过程中，初始连接权值的和阈值对训练结果影响很大，但由于随机选取的方法，传统人工神经网络存在收敛速度慢、易陷入局部最优的缺点[14]；遗传算法是模拟自然选择和遗传学机理的最优解搜索技术，具有较强的全局搜索能力。为了有效克服传统人工神经网络的缺点，本研究将进化计算的自适应机制与神经网络的学习机制有机结合，使用遗传算法对神经网络参数进行优化：在NN基础上，以网络实际输出与期望输出间的均方误差构造适应度函数，通过遗传算法的进化操作和全局搜索能力，获取表现最优进化初始权值和阈值；利用神经网络的局部寻优能力，训练出性能更优的燃油流量估算模型。改进的神经网络（improved neural network, INN）算法流程，见图 3。

图3 遗传算法优化的改进神经网络算法流程图

2.3 基于运行数据的网络训练

使用实际运行数据进行网络训练。针对收集到的B737-800机型的实际飞行记录数据进行清洗和处理，原始数据包括飞机高度、马赫数、地速、时间、左发燃油流量和右发燃油流量等信息。根据发动机转速、爬升率等信息，提取出所需的飞机巡航阶段相关数据（共3 025条），将其随机分为训练集（70%）、测试集（20%）和验证集（10%）。

2.3.1 隐含层神经元数目和学习算法的确定 隐含层神经元数目和学习算法是决定网络性能的重要参数。为尽可能提升网络性能，本研究在学习率为 0.01、训练目标为 1×10−3、最大训练次数为 1 000、f1使用Tansig函数、f2使用Purelin函数的条件下，针对隐含层神经元数目为6～20间的不同取值，学习算法分别使用梯度下降（Traingd）、贝叶斯正则化（Trainbr）和 Levenberg-Marquardt（Trainlm）的方法对NN模型进行试验。图4为各自50次训练中测试集实际输出与期望输出间均方误差（mean square error, MSE）的平均结果，当隐含层神经元数目取18，使用贝叶斯正则化学习算法时网络精度较高。

图4 不同隐含层神经元数目和学习算法组合下的神经网络性能比较

2.3.2 进化初始权值和阈值的获取 在优选的隐含层神经元数目和学习算法的基础上，通过2.2中的算法流程，对网络初始权值和阈值进行优化。遗传算法部分使用二进制进行编码；以网络实际输出和期望输出的均方误差为适应度函数；以轮盘选择、两点交叉、离散变异和精英保留的策略进行进化操作；经过多次试验，综合考虑设置遗传算法参数：种群大小定义为100，最大遗传代数为100，变量的二进制位数为16，交叉、变异概率分别为0.85和0.01，代沟为0.9。通过遗传算法优化的最佳适应度曲线，见图5，当进化到第42代时，测试集的最小MSE开始收敛，此时进化初始权值和阈值。

图5 遗传算法适应度曲线

使用进化初始权值和阈值对网络进行训练，可得出精度更高的巡航燃油流量估算INN模型。使用INN时测试集各个数据点实际输出与期望输出间绝对值相对百分比误差σ的统计情况，见图6。

图6 测试集实际输出与期望输出对比

图6可知，74%的数据点绝对值相对百分比误差在5%以内，初步说明模型取得了较好的训练效果。

3 网络性能评价与案例分析

3.1 模型评价

从验证集数据中，随机选取100个样本进行实验，为确保结果可靠，样本选取分5次进行，形成的G1、G2、G3、G4和G5共5组数据的网络输出情况，见图7。总体来看，5组数据均体现出优异的拟合效果，说明INN具有良好的泛化能力。

图7 验证数据实际输出与期望输出对比

为评估遗传算法的优化效果，针对上述5组数据，分别使用NN和INN进行实验，通过拟合优度（R2）、均方误差（MSE）、相对百分比误差的最大值（max-σ）和平均值（mean-σ）4 项指标比较其预测性能，对比结果见表2。

表2 NN与INN性能指标对比

表2可知，5组数据输出结果中，INN的R2均在 0.95以上，平均R2高达0.973 1，大于 NN的0.954 6；MSE 平均值由 0.002 7 降低至 0.001 4，说明了遗传算法的优化有效提升了网络的精度，使得燃油流量的获取更加可靠；max-σ和mean-σ的平均值分别降低了43.45%和29.30%，表明优化后的INN具备更加稳定的预测性能。

3.2 算例分析

根据综合温室效应计算模型和巡航燃油流量神经网络模型，通过公开的航迹和气象信息，可实现航空器巡航温室效应的定量计算。本研究以B737-800机型为例，选取2019年10月7日执飞成都→北京航段的ZH9164航班，航空器巡航温室效应评估框架进行算例计算和分析，见图8。

图8 航空器巡航温室效应评估框架

已知执飞该航班的航空器机型为B737-800，根据该机型选装发动机的NOx基准排放指数（表1），利用1stOpt软件工具，使用通用全局优化算法，获取该机型基准燃油流量在[300, 5 000]取值范围，其与NOx基准排放指数间的拟合关系，见式（8）。

式（8）中，x为NOx的基准排放指数；p1、p2、p3和p4为系数，系数取值与回归统计，见表3。残差平方和与均方误差均处于较低水平，拟合效果较好，决定系数达到 0.999 9，可见式（8）的拟合关系较为可靠。由此，可以通过修正后的基准燃油流量，获取NO基准排放指数

表3 拟合系数取值及回归统计

从VariFlight飞行轨迹实时跟踪雷达记录中，获取该航班巡航阶段的ADS-B数据（包括航班号、应答机编码、时刻、航向、气压高度、地速、经度和纬度）。根据 Real World Flightplan Database 数据，该日成都→北京航段的巡航飞行中，飞机途经的重要航路点位置信息，见表4。

表4 成都至北京航路上重要航路点位置信息

针对气象数据，本研究以每段航路中点的气象数据作为该段航路的气象条件标准。根据航路点的位置信息，从美国怀俄明州立大学发布的高空探空气象数据（包括大气压强、温度、相对湿度、风速和风向等）中，获取航路周围探测站点的气象数据并以此为样本点，使用反距离加权（inverse distance weighting, IDW）的空间插值方法，获得航迹中各点的气象信息。

将航迹数据与气象数据进行归一化处理，输入到巡航燃油流量估算INN模型中，获取与雷达记录时刻点对应的燃油流量序列；进一步通过式（7）可实现航空器巡航温室效应定量评价，见表5。

表5 算例计算结果及对比

表5列出了现有文献[15]中利用基于BADA性能模型的仿真软件，以B737-800机型为例，在初始巡航质量 67 t、巡航高度 9 500 m、航程 1 500 km条件下仿真计算出的巡航油耗及温室气体排放量数据。可以看出，在初始条件近似的情况下，通过2种方法获得的油耗和气体排放结果较为接近，说明了本研究提出的计算方法在减少对飞行参数依赖的同时具备可靠性，且能进一步给出气体排放产生的具体温变影响。

4 讨论与结论

本研究通过建立综合温室效应评估模型和巡航燃油流量神经网络模型，实现了对航空器巡航阶段污染物排放温变影响的定量计算，结论如下。

（1）基于航空器实际运行数据，使用机器学习的方法构建三层神经网络用于巡航燃油流量获取，与现有文献中常用的通过BADA性能模型计算的方法相比，本研究提出的方法具备操作简便、符合实际运行特点的优势。

（2）引入遗传算法优化经典神经网络中的初始权值和阈值，改进后的神经网络在精确性和稳定性方面均有明显提升，由此获取的油耗结果更加可靠，为进一步的航空器巡航污染物排放量计算和温室效应评估奠定基础。

（3）使用本研究提出的方法，能够将航空器巡航阶段不同污染物排放的影响进行统一度量，相比于污染物排放清单，给出的具体温变结果在环境影响评估深度方面有所突破。

本研究成果可为实际应用中的燃油消耗计算和环境影响评价提供新思路，在未来的研究中，可将凝结尾等更多的温室效应影响因素纳入考虑范围，以进一步分析现代航空运输的环境影响和可持续发展能力。