, 李玉铭中国民航大学航空工程学院, 天津 300300
收稿日期: 2020-05-26; 修回日期: 2020-07-16; 录用日期: 2020-07-16
基金项目: 中国民航大学开放基金（No.000031020102）
作者简介: 曹惠玲(1962-), 女, 教授, E-mail:hlcao@cauc.edu.cn
通讯作者（责任作者）: 曹惠玲
摘要：民航飞机在巡航阶段排放的黑碳（BC）可产生区域增温效应.为提高飞机巡航阶段黑碳排放的估算精度，采用一种基于国际民航组织一阶近似方法（FOA）的修正计算方法（FOA-MOD），使用B777-200机载飞行数据（QAR）进行实例分析，并结合当量比、燃烧室入口压力、主燃区火焰温度等热力学参数，计算了某次飞行巡航阶段的BC排放总量.同时，在排放指数和燃油流量的基础上引入排放强度的概念，以更全面地表征污染物排放特性.最后分析了巡航高度、巡航马赫数和巡航重量等飞机性能参数对排放特性的影响.结果表明：修正方法的计算结果更加贴近实际，相对一阶近似方法的计算结果降低19.28%；通过增加巡航高度、减小巡航马赫数和巡航重量，可以不同程度地降低燃油消耗和污染物排放.本研究所做工作能够为评估飞机污染物排放特性，提高飞机运行的经济性和环保性提供参考.
关键词：民航飞机巡航阶段修正计算方法排放特性飞机性能
Research on BC emission characteristics and influencing factors of aircraft during cruise phase based on FOA-MOD
CAO Huiling
, LI YumingCollege of Aeronautical Engineering, Civil Aviation University of China, Tianjin 300300
Received 26 May 2020; received in revised from 16 July 2020; accepted 16 July 2020
Abstract: The black carbon (BC) emitted by civil aviation aircraft during the cruise phase could produce a regional warming effect. In order to improve the estimation accuracy of the BC emissions during the cruise phase of the aircraft, a modified calculation method (FOA-MOD) based on the ICAO first order approximation method (FOA) is used; the B777-200 airborne flight data (QAR) is used in a case analysis, where the total BC emissions of a certain cruise phase are calculated by combining thermodynamic parameters such as equivalent ratio, combustor inlet pressure and flame temperature in the main combustion zone; the concept of emission intensity is introduced on the basis of emission index and fuel flow to characterize pollutant emission characteristics more comprehensively; the effects of aircraft performance parameters such as cruise altitude, cruise Mach number, and cruise weight on emission characteristics are analyzed. The result shows that the calculation result using FOA-MOD method is closer to reality, with a 19.28% reduction relatively comparing to FOA method. By increasing the cruise altitude, reducing the cruise Mach number and cruise weight, the fuel consumption and pollutant emissions can be reduced to varying degrees. The work can provide a reference for evaluating aircraft pollutant emission characteristics and improving the economic and environmental protection of aircraft operations.
Keywords: civil aircraftcruise phasemodified calculation methodemission characteristicsaircraft performance
1 引言(Introduction)民航飞机排放的氮氧化物(NOx)、一氧化碳(CO)、未燃碳氢化合物(UHC)、硫氧化物(SOx)与颗粒物(PM)等污染物(ICAO, 2011), 能够对全球气候变化、局部空气质量产生影响.民航飞机在“安全性、经济性、环保性、舒适性”(简称“四性”)等方面有特殊要求, 就“环保性”而言, 必须满足适航的强制性规定(张驰等, 2014).
作为民航发动机排放的适航标准, CCAR-34部(涡轮发动机飞机燃油排泄和排气排出物规定)规定了各种污染物的适航标准.目前我国正在进行国产商用飞机的设计研发工作, 为了提高飞机在国际市场的竞争力, 有必要在设计和适航验证阶段对飞机排放的污染物做出准确评估.由于专用的测试程序复杂, 且对测试设备的要求较高, 因此, 寻求一种利用飞行数据针对在役或适航验证阶段的航空器排放适航符合性验证的等效替代方法值得探讨, 其结果对设计阶段的航空器也具有较高的参考价值.
目前, 国内外****对飞机污染物排放方面的研究主要集中在排放量的估算上:针对起飞着陆循环(Landing and Take-off, LTO)排放的污染物对机场周边环境影响的研究较多(王瑞宁等, 2018；韩博等, 2020), 针对高空巡航阶段排放的污染物对气候变化影响的研究较少.黑碳(Black Carbon, BC)气溶胶是飞机排放颗粒物的重要组成成分, 由于飞机需要经历长时间的巡航飞行, 该阶段排放的黑碳在大气层形成的辐射强迫(Radiative Forcing, RF)能够极大地减弱对地球的冷却效果, 产生温室效应(Ramanthan et al., 2008), 且飞行中绝大多数的燃油消耗也发生在该阶段.因此, 评估飞机巡航阶段的黑碳排放, 提高飞机运行的经济性和环保性具有重要的意义.
在飞机黑碳排放的研究方面, Stettler等(2013)开发了基于黑碳形成和氧化过程的形成氧化法, 使飞机的黑碳排放估算精度相比国际民航组织一阶近似方法有所提高；Petzold等(1999)开发了估算飞机发动机在给定巡航高度和速度下黑碳排放指数(每燃烧1 kg燃油的黑碳排放量)的算法；Dopelheuer等(1998)分析了发动机性能参数、大气条件、飞行剖面等因素对飞机黑碳排放的影响.我国在20世纪80年代后也开展了相关研究, 目前主要集中于城乡大气方面, 关注较多的领域是黑碳的来源分析(井安康等, 2019)、浓度分布(魏夏潞等, 2019)和辐射强迫(Miyazaway et al., 2013).
飞机排放的污染物主要是排放指数(Emission Index, EI)、运行时间(Time In Mode, TIM)和燃油流量(Fuel Flow, FF)3个参数的函数.飞机运行过程中, 排放指数会随发动机实际运行工况发生变化, 运行时间和燃油流量与国际民航组织推荐的发动机排放数据库(Engine Emission Data Bank, EEDB)标准值也存在较大差异(曹惠玲等, 2018).
快速存取记录器(Quick Access Recorder, QAR)中既包含飞机的运行时间, 还包含表征操纵飞机和发动机的诸多参数, 如燃油流量、推力值等, 能够为待评估发动机排放量的准确评估提供数据支撑.由于其记录的参数种类齐全, 现阶段也广泛用于航空发动机的状态监控(孙见忠等, 2013)、故障诊断(钟诗胜等, 2019)与趋势预测(李艳军等, 2017).
本研究为更加全面地考虑发动机性能差异和飞行操纵差异对排放的影响, 结合当量比(Equivalent Ratio)、燃烧室入口压力、主燃区火焰温度等影响排放的参数, 采用基于一阶近似方法的修正计算方法, 对B777-200飞机某次飞行巡航阶段的黑碳排放总量进行计算.在排放指数和燃油流量的基础上, 引入排放强度(Emission Intensity, EI′)的概念来综合表征排放特性, 在此基础上对巡航高度、巡航马赫数和巡航重量等影响排放的性能参数进行讨论.以期为评估飞机污染物排放, 开展飞机的节能减排工作提供理论支撑.
2 排放特性计算方法(Calculation method of emission characteristics)民航飞机排放的颗粒物主要包括非挥发性颗粒物(PMnvol, 主要成分是黑碳)、挥发性硫酸盐颗粒物(PMvol-FSC, 由燃料中的硫成分导致)和挥发性有机颗粒物(PMvol-FuelOrganics, 由燃料的不完全燃烧导致)3种.本研究主要针对PMnvol展开计算, 其粒径约为0.02~0.06 μm, 属于PM2.5范畴.
2.1 排放指数2.1.1 一阶近似方法一阶近似方法(First Order Approximation, FOA)是航空环境保护委员会制定并批准用于估算飞机起飞着陆循环颗粒物排放的方法.它基于国际民航组织公布的发动机在地面试车台不同推力等级下的污染物排放数据, 可以作为下文中地面基准状态黑碳浓度的估算依据.首先通过实验关联式将烟度值(Smoke Number, SN)转换成黑碳浓度CBC(mg·m-3), 然后使用空气燃料比(Air Fuel Ratio, AFR)和涵道比(Bypass Ratio, BPR)得到每千克燃油的废气容积流量QMixed(m3·kg-1), 最后由黑碳浓度和废气容积流量得出黑碳排放指数EIBC(mg·kg-1).
 | (1) |
 | (2) |
 | (3) |
本研究为有效提高黑碳排放的估算精度, 首先利用发动机排放数据库4个推力级别下的烟度值、空燃比和燃油流量分别建立零推力到全推力范围内SN~F00、AFR~F00和FF~F00的函数关系, 然后通过统计飞机发动机的运行数据得到巡航阶段推力值的频数分布, 将频数分布最大的推力值(用Fa表示)所对应的参数值作为估算依据.下文以B777-200飞机的PW4077D型发动机排放计算为例进行研究.
2.1.2 修正计算方法一阶近似方法可利用国际民航组织公布的不同推力等级下的污染物排放数据, 对各飞行阶段的污染物进行评估, 但因其未考虑发动机的个体差异及操纵品质, 因此只用于统计意义上的快速评估.为提高个体发动机黑碳排放的估算精度, 且无需使用专用的取样及测试设备和程序, 而是借助于实际飞行数据, 采用一种修正方法(Modified Calculation Method, FOA-MOD)可提高飞机巡航阶段黑碳浓度计算精度.它在一阶近似方法的基础上, 考虑了发动机性能及外界环境条件对排放特性的影响, 同时将燃烧室入口压力P3、当量比Φ、主燃区火焰温度Tfl等与燃烧品质直接相关的参数考虑在内, 可以在已知地面基准状态下各参数值的基础上, 对不同寿命周期(新服役或性能退化)的飞机发动机黑碳浓度排放值进行修正, 使计算值更加接近排放的真实值(Stettler et al., 2013):
 | (4) |
如2.1.1节所述, 选取推力值为Fa的地面车台实验状态点作为巡航基准状态, 使用式(1)计算得到CBC, ref.由于实验程序复杂, 且对测试设备的要求较高, 因此为获取P3, ref值, 基于大量的历史飞行数据, 建立了零推力到全推力范围内的燃烧室入口压力P3与推力值F/F00之间的函数关系:在零推力时, P3近似等于大气压力；在全推力时, P3主要由压气机增压比π00决定.而推力值为Fa对应的P3即为所需的P3, ref:
 | (5) |
为获取地面基准状态下燃烧室主燃区火焰温度Tfl, ref, 首先需要计算燃烧室入口燃气温度T3, ref.在得到P3, ref后, 根据热力过程方程, 可计算T3, ref:
 | (6) |
下面可以计算主燃区火焰温度Tfl.由于其与燃烧室入口温度T3直接相关, 从航空燃料(C12H22)的燃烧反应方程式出发, 假定完全燃烧, 不发生离解现象, 结合燃烧反应物和产物(O2、N2、CO2、H2O)在不同温度下的绝对焓值, 可得到主燃区火焰温度Tfl与燃烧室入口温度T3的线性关系(R2=1):
 | (7) |
 | (8) |
 | (9) |
2.2 排放量黑碳排放量由发动机运行时间、燃油流量和黑碳排放指数三者的乘积决定.为确定巡航阶段的实际运行时间, 结合飞行数据中飞行阶段指示(Flight Phase)、低压转子转速(N1)、飞行高度(Altitude)、飞行马赫数(Mach)等参数对巡航阶段进行准确划分.
采用一阶近似方法计算黑碳排放量时, 使用频数分布最大的推力值(Fa)所对应的燃油流量(FF)与排放指数(EIBC)及划分后的巡航阶段总工作时间(T)做乘积, 得到巡航阶段的黑碳排放量Em:
 | (10) |
 | (11) |
2.3 排放强度为更加全面地表征黑碳排放量随时间的变化特性, 以排放指数和燃油流量组成的综合指标排放强度来评估黑碳排放特性.排放强度定义为每飞行单位距离的黑碳排放量, 计算时可通过排放指数和燃油流量的乘积除以真实空速(True Air Speed, TAS)得出, 具体如式(12)所示.
 | (12) |
2.4 对比说明黑碳浓度值会直接影响黑碳排放指数的计算值, 进而影响黑碳排放量的计算结果.因此, CBC在计算过程中尤为关键.对实际飞行的各飞行阶段而言, 并不像排放数据库推荐的那样始终保持推力固定在某一数值, 实际上不同时刻的推力值是上下波动的.因此, 反映黑碳排放特性的各个参数都会随发动机实际运行工况发生变化.采用一阶近似方法计算黑碳浓度值时, 将烟度值作为独立变量, 通过经验公式将烟度值转换成黑碳浓度, 未考虑发动机的个体差异和操纵品质；而修正方法用一阶近似方法首先计算地面基准状态黑碳浓度, 在此基础上, 结合实际飞行数据中与发动机性能、燃烧品质、环境条件相关的热力学参数, 对前面得到的黑碳浓度值进行修正, 从而提高黑碳排放的估算精度, 计算流程见图 1.
图 1(Fig. 1)
 |
| 图 1 黑碳排放量计算流程 Fig. 1Black carbon emissions calculation flowchart |
3 案例分析(Case Analysis)为说明黑碳排放量计算及修正方法, 以PW4077D型商用航空涡轮风扇发动机(Turbo Fan, TF)为例计算某次飞行巡航阶段的黑碳排放量.
3.1 一阶近似方法PW4077D型航空发动机的设计参数及LTO循环各阶段的推荐烟度值、空燃比、燃油流量等数据见表 1, 这些数据取自发动机排放数据库和机场空气质量手册, 可为一阶近似方法的计算提供数据支撑.
表 1(Table 1)
| 表 1 PW4077D发动机参数 Table 1 PW4077D engine parameters | |||||||||||||||||||||||||||||||
表 1 PW4077D发动机参数 Table 1 PW4077D engine parameters
| |||||||||||||||||||||||||||||||
 |
因排放数据库不包含巡航阶段的标准工作时间, 为获取准确有效的巡航数据, 采用严格的巡航数据筛选标准, 结合实际飞行数据中飞行阶段指示、飞行高度、低压转子转速、飞行马赫数等参数描述划分后的巡航阶段特征, 具体如表 2所示.本研究计算黑碳排放量时均采用此工作时间.由表 2可见, 巡航阶段的实际运行参数并非稳定在某一数值, 而是在一定范围内波动.
表 2(Table 2)
| 表 2 巡航阶段特征参数 Table 2 Characteristic parameters of cruise phase | ||||||||||||
表 2 巡航阶段特征参数 Table 2 Characteristic parameters of cruise phase
| ||||||||||||
|
如2.1.1节所述, 使用式(1)~(3), 结合表 1中的相关参数, 分别使用爬升阶段及频数分布最大的推力值(经统计, Fa=79.50%F00, F00为发动机的全推力值)所对应的烟度值和空燃比, 计算巡航阶段的黑碳排放指数.为获取Fa所对应的SN和AFR, 首先利用表 1的推力设置、烟度值和空燃比分别建立零推力到全推力范围内SN~F00、AFR~F00的函数关系, 然后将Fa所对应的SN和AFR作为估算依据.由于计算时使用的烟度值和空燃比均为固定值, 因此, 巡航阶段的排放指数也为恒定值.
为得到巡航阶段频数分布最大的推力值所对应的燃油流量值, 同样需要使用表 1中的推力设置和燃油流量建立零推力到全推力范围内FF~F00的函数关系, 然后将Fa所对应的FF作为估算依据.将两种推力设置下(85%F00、79.50%F00)的相关参数及计算结果进行对比, 结果如表 3所示.
表 3(Table 3)
| 表 3 不同推力设置下的计算结果对比 Table 3 Comparison of calculation results under different thrust settings | |||||||||||||||||||||
表 3 不同推力设置下的计算结果对比 Table 3 Comparison of calculation results under different thrust settings
| |||||||||||||||||||||
|
85%F00为很多文章推荐采用的巡航推力, 为使巡航状态下排放的计算结果更加符合实际, 本文推荐采用每次飞行实际使用最多的巡航推力.由表 3可见, 黑碳排放量随着推力设置的下降而下降, 两种推力设置下黑碳排放量的相对偏差为19.02%.上述计算过程中均假定巡航阶段推力稳定在某一固定数值, 而实际飞行过程中不同时刻的推力值是上下波动的, 反映排放特性的各个参数都会随发动机实际运行工况发生变化, 加之一阶近似方法没有考虑发动机性能和环境状态的影响, 因此仅可用于快速评估.为更精确地计算飞机巡航阶段的黑碳排放指数, 需要采用下面的修正计算方法.
3.2 修正计算方法如2.1节所述, 使用式(1)~(9), 结合当量比、燃烧室入口压力、主燃区火焰温度等热力学参数, 计算飞机巡航阶段的黑碳排放指数, 修正方法同一阶近似方法的计算结果对比见图 2a.由图 2a可见, 修正方法与一阶近似方法计算的黑碳排放指数存在较大差异.首先对比算例1和算例2, 二者的本质区别在于推力设置不同, 可以看出排放指数随推力设置的下降而下降, 说明黑碳排放指数主要受推力设置的影响.然后对比算例2和算例3, 可以发现修正方法计算的黑碳排放指数在一阶近似方法计算值的附近波动, 且波动幅度较大, 表明修正方法能够更好地反映黑碳排放指数随飞行实时工况的变化, 同时将发动机性能差异及外界环境条件对排放特性的影响考虑在内, 最终使计算结果更加接近实际的排放值.
图 2(Fig. 2)
 |
| 图 2 排放指数(a)和燃油流量(b)变化趋势 Fig. 2Emission index(a) and fuel flow(b)trend chart |
为了进一步讨论排放总量, 还需要知道燃油流量的变化.图 2b所示为燃油流量随时间的变化趋势图.算例1~3分别对应一阶近似方法的燃油流量推荐值、频数分布最大的推力值所对应的燃油流量值和飞机实际运行数据中的燃油流量值.算例2和算例1的区别在于推力设置的不同, 说明燃油流量也主要受推力设置的影响, 且随推力设置的下降而下降.由算例3可见, 巡航阶段的燃油流量在0.60~0.85 kg·s-1的范围内波动.由于不同时刻的低压转子转速是波动的, 因此, 燃油流量也会随之波动, 波动情况能够充分反映飞行的实时工况.不同算例下的燃油流量差异及形成原因与排放指数相对应, 均由实际运行工况与标准工况下的推力设置导致.
3.3 黑碳排放量在已知排放指数、燃油流量和工作时间的基础上, 使用式(10)~(11)计算巡航阶段的黑碳排放总量, 计算结果对比如表 4所示.由表 4可见, 修正方法和两种推力设置下的一阶近似方法计算结果均存在较大差异.由于3个算例使用的工作时间相同, 因此, 排放量的差异是排放指数和燃油流量二者共同作用的结果.
表 4(Table 4)
| 表 4 巡航阶段黑碳排放计算结果 Table 4 Calculation results of black carbon emissions during cruise phase | ||||||||||||||||||||||||||||||||||||||||||||||||
表 4 巡航阶段黑碳排放计算结果 Table 4 Calculation results of black carbon emissions during cruise phase
| ||||||||||||||||||||||||||||||||||||||||||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||
算例2和算例1的本质区别在于推力设置的不同.为了使计算结果更加接近实际, 应选用每次飞行实际采用最多的巡航推力设置.算例3的平均排放指数、平均燃油流量和黑碳排放量均低于算例2, 相对偏差分别为-3.88%、-4.97%和-19.28%.一阶近似方法选用频数分布最大的推力值作为估算依据, 能够快速建立巡航阶段黑碳排放的统计性计算；而修正方法结合与发动机性能、燃烧品质、环境条件相关的热力学参数, 能够更好地反映出黑碳排放特性随实时工况的变化, 进一步提高了黑碳排放的估算精度, 因此, 计算结果偏差反映的是实际运行工况与选用标准工况的差异.
为了进一步探讨黑碳排放的变化规律及影响因素, 前面引入了排放强度的概念, 下面就此问题展开进一步讨论.
4 影响因素分析(Analysis of influencing factors)巡航是飞机的主要飞行阶段之一, 通过优化巡航性能参数可以节省燃油消耗和降低污染物排放, 提高飞机运行的经济型和环保性, 具有重要的意义.飞机发动机型号、发动机燃烧室类型一定的情况下, 飞机巡航高度、巡航马赫数和巡航重量等运行参数会对排放特性产生影响(Turgut et al., 2017), 同时这些参数与飞行操纵品质直接相关, 因此, 研究其对黑碳排放特性的影响规律, 可以通过优化飞行操纵程序达到节能减排的效果.
选取B777-200飞机的历史飞行数据作为数据源, 运行期间发动机没有因性能衰退、时寿件到寿、部件损伤等原因拆下送修, 因此, 不考虑发动机型号变更等因素对黑碳排放特性的影响.统计数据中, 飞机运行期间的巡航高度所在区间为24000~34000 ft, 巡航马赫数所在区间为0.75~0.90, 巡航重量所在区间为175~195 t.选取2000 ft作为巡航高度的组距, 0.03作为巡航马赫数的组距, 5 t作为巡航重量的组距, 分析3个参数的变化对黑碳排放特性的影响.
由于排放强度综合了排放指数和燃油流量的影响, 下面在保证其他运行参数不变的基础上, 计算不同巡航高度、巡航重量和巡航马赫数下的排放强度, 并分析其变化规律.图 3为不同巡航重量和不同巡航马赫数下排放强度随飞行高度的变化规律.由图 3a可见, 其他条件不变, 随着巡航高度的增加, 排放强度呈下降趋势.巡航高度每增加2000 ft, 排放强度下降约2.2%.巡航高度越高, 空气越稀薄, 空气密度越小, 飞行阻力越小, 在同等工况下推力设置越小, 燃油流量和排放指数下降, 排放强度随之下降.
图 3(Fig. 3)
 |
| 图 3 排放强度变化趋势 Fig. 3Emission intensity trend chart |
在其他条件不变的前提下, 排放强度随巡航重量的增加而增大, 巡航重量每增加5 t, 排放强度上升约3.7%.巡航重量越大, 在同等工况下推力设置越大, 燃油流量和排放指数增大, 排放强度随之增大.
由图 3b可见, 其他条件不变, 随着巡航马赫数的增加, 排放强度呈上升趋势, 巡航马赫数每增加0.03, 排放强度上升约3.1%.巡航马赫数越大, 相应时刻的飞行速度越大, 在同等工况下推力设置越大, 燃油流量和排放指数增大, 排放强度随之增大.研究结果可为飞机巡航性能参数优化, 提高航班运行的经济环保性提供参考.
5 结论(Conclusions)1) 结合飞行数据的修正计算方法能够充分反映飞行的实时工况.采用修正方法计算的黑碳排放量低于一阶近似方法, 相对偏差为-19.28%.黑碳排放量是排放指数、燃油流量和运行时间3个参数的函数.由于计算过程中使用的工作时间相同, 因此, 结果偏差是排放指数和燃油流量二者作用的结果.对各飞行阶段而言, 并不像排放数据库推荐的那样始终保持推力固定在某一数值, 实际上不同时刻的推力值是上下波动的, 因此, 计算结果偏差反映的是实际运行工况与标准工况的差异.
2) 对于确定的发动机型号而言, 燃油流量和排放指数主要取决于发动机的推力状态.基于实际飞行数据, 使用频数分布最大的推力值作为估算依据, 能够使影响排放特性的推力设置取值更加合理, 从而提高黑碳排放快速估算的精度.
3) 飞机巡航高度、巡航马赫数和巡航重量等运行参数对黑碳排放特性影响显著.排放强度随巡航重量和巡航马赫数的增大而增大, 随巡航高度的增大而减小.实际运行过程中, 通过优化飞行操纵程序, 适当增加巡航高度、减小巡航马赫数和巡航重量, 可以不同程度地节省燃油消耗和降低污染物排放, 提高飞机运行的经济性和环保性.
4) 一阶近似方法基于大量飞行数据, 能够快速建立飞机颗粒物排放的统计性计算与分析, 可用于飞机排放的整体快速评估.修正方法考虑了飞机发动机型号、发动机性能差异、外界环境条件、飞行操纵品质、燃烧品质等诸多因素, 能够真实准确地评估单架飞机巡航阶段的黑碳排放量.这样不仅可用于发动机排放的精确评估, 也可为今后进一步进行发动机排放适航符合性验证提供一种等效替代方案.
参考文献
| 曹惠玲, 汤鑫豪, 苗佳禾. 2018. 基于QAR数据的发动机LTO阶段氮氧化物排放量计算与分析[J]. 环境科学学报, 38(10): 3900-3904. |
| Dopelheuer A, Lecht M.1998.Influence of engine performance on emission characteristics[C].RTO AVT Symposium on"Gas Turbine Engine Combustion, Emissions and Alternative Fuels"Proceedings, RTO-MP-14.Lisbon |
| 韩博, 孔魏凯, 姚婷玮, 等. 2020. 京津冀机场群飞机LTO大气污染物排放清单[J]. 环境科学, 41(3): 1143-1150. |
| ICAO.2011.Airport air quality manual[R].Montreal: International Civil Aviation Organization |
| ICAO.2013.Engine Emission Data Bank[R]. Doc 9646-AN/943. Montreal: International Civil Aviation Organization. |
| 井安康, 朱彬, 丁德平, 等. 2019. 中国长江三角洲地区黑碳特征和来源分析[J]. 中国环境科学, 39(9): 3585-3594. DOI:10.3969/j.issn.1000-6923.2019.09.001 |
| 李艳军, 张建, 曹愈远, 等. 2017. 基于模糊信息粒化和优化SVM的航空发动机性能趋势预测[J]. 航空动力学报, 32(12): 3022-3030. |
| Miyazawa Y, Harada A, Wickramasinghe N K, et al. 2013. Effect of jet passenger aircraft performance model on the optimal periodic cruise maneuver[J]. Space Technology Japan the Japan Society for Aeronautical & Spaceences, 12: 99-105. |
| Petzold A. 1999. In situ observations and model calculations of black carbon emission by aircraft at cruise altitude[J]. Journal of Geophysical Research:Atmospheres, 104(D18): 22171-22181. DOI:10.1029/1999JD900460 |
| Ramanthan V, Carmichael G. 2008. Global and regional climate changes due to black carbon[J]. Nature, 1: 221-227. |
| Stettler M E J, Boies A M, Petzold A, et al. 2013. Global civil aviation black carbon emissions[J]. Environmental Science & Technology, 47(18): 10397-10404. |
| 孙见忠, 左洪福, 刘鹏鹏, 等. 2013. 基于动态线性模型的民航发动机性能状态监控与在翼寿命预测[J]. 系统工程理论与实践, 33(12): 3243-3250. DOI:10.12011/1000-6788(2013)12-3243 |
| 孙见忠, 左洪福, 刘鹏鹏, 等. 2012. 航空发动机污染物排放量估算方法[J]. 交通运输工程学报, 12(2): 53-61. |
| Turgut E T, Usanmaz O. 2017. An assessment of cruise Nox emissions of short-haul commercial flights[J]. Atmospheric Environment, 171: 191-204. DOI:10.1016/j.atmosenv.2017.10.013 |
| 王瑞宁, 黄成, 任洪娟, 等. 2018. 长三角地区民航飞机起飞着陆(LTO)循环大气污染物排放清单[J]. 环境科学学报, 38(11): 4472-4479. |
| 魏夏潞, 王成刚, 凌新峰, 等. 2019. 安徽寿县黑碳气溶胶浓度观测分析研究[J]. 环境科学学报, 39(11): 3630-3638. |
| 张弛, 林宇震, 徐华胜, 等. 2014. 民用航空发动机低排放燃烧室技术发展现状及水平[J]. 航空学报, 35(2): 332-350. |
| 钟诗胜, 李旭, 张永健. 2019. 基于DBN的不均衡样本驱动民航发动机故障诊断[J]. 航空动力学报, 34(3): 708-716. |
