上海市环境科学研究院, 国家环境保护城市大气复合污染成因与防治重点实验室, 上海 200233
收稿日期: 2018-10-23; 修回日期: 2018-12-29; 录用日期: 2018-12-29
基金项目: 国家重点研发计划项目（No.2017YFC0212100）；上海市科委项目（No.16DZ1206704，18DZ1203100）
作者简介: 胡磬遥(1986-), 男, E-mail:huqy@saes.sh.cn
通讯作者（责任作者）: 胡磬遥
摘要: 为了能够低成本地获得包括工程机械在内的非道路移动机械的NOx排放因子，基于目前广泛使用的车用氮氧化物传感器开展了NOx排放因子的测量系统搭建和计算模型研究.基于柴油机尾气O2浓度、NOx浓度、环境温湿度和燃料特性等参数建立了非道路移动机械单位油耗NOx排放因子计算模型.为了验证测量系统与计算模型，开展了4台工程机械的实际道路验证试验，并选取其中一台开展了基于车载排放测试系统（PEMS）的NOx排放比对试验.比对和实测试验结果表明建立的NOx排放因子测量方法具有较低的相对误差和稳定性，与车载排放测量系统的测量结果相比具有很好的相关性.被测2台挖掘机和1台推土机、1台压路机的平均单位油耗NOx排放因子分别为（22.8±4.0）、（35.8±16.0）和（55.0±18.0）、（68.8±13.0）g·kg-1.基于车用氮氧化物传感器的NOx排放因子测量结果与目前国内外已有研究的NOx排放因子测量结果接近.
关键词:车用氮氧化物传感器工程机械NOx排放因子测量系统
Research on construction machinery NOx emission factor measurement method based on vehicular NOx sensor
HU Qingyao
Shanghai Academy of Environmental Sciences, State Environment Protection Key Laboratory of Formation and Prevention of the Urban Air Complex, Shanghai 200233
Received 23 October 2018; received in revised from 29 December 2018; accepted 29 December 2018
Abstract: In order to obtain the NOx emission factors of non-road mobile machinery which including construction machinery at low cost, the establishment and algorithm of the measurement system of transient fuel-specific NOx emission factor is carried out based on the currently widely used vehicular NOx sensor. The formula of instantaneous fuel specific NOx emission factor was established by the oxygen concentration, nitrogen oxide concentration and ambient temperature, humidity. In order to verify the measurement system and the algorithm, the real operation condition test of four construction machineries were carried out, and one of them was selected to carry out the emission comparison experiment based on the portable emission measurement system (PEMS). The comparison and the verification experiments results show that the established NOx emission factor measurement method has low relative error and excellent stability. The established measurement method has a good correlation with the measurement results of the PEMS. The average fuel-specific NOx emission factors of the two excavators, one excavator and one roller measured as shown in the verification test results are (32.4±4.0), (35.8±16.0) and (55.0±18.0), (68.8±13.0) g·kg-1 respectively. NOx emission factors obtained by the measurement method based on the vehicle NOx sensor are similar to the existing worldwide research results.
Keywords: vehicular NOx sensorconstruction machineryNOx emission factorsmeasurement system
1 引言(Introduction)挖掘机、推土机、压路机等以柴油机为动力的工程机械, 在基建领域有着不可替代的作用(Zhao et al., 2013).由于柴油机的运行特点是缸内直喷、稀薄半预混非均质燃烧, 导致其燃烧过程高温、富氧且燃烧反应时间充分, 更有利于NOx的产生.包括工程机械在内的各类非道路移动柴油机械是我国城市NOx的重要来源之一, 已成为制约我国大城市空气质量持续改善的重要因素, 对城市大气环境和居民健康造成了重大的威胁(张楚莹等, 2008; Anenberg et al., 2017).如何准确、高效地获得各类工程机械的NOx排放水平, 是开展工程机械等非道路移动机械NOx排放治理的基础工作(张礼俊等, 2010; 李东玲等, 2012; Wang et al., 2016; 鲁君等, 2017; 黄成等, 2018).
目前大多数非道路移动机械排放因子实测研究都是基于车载排放测量系统(PEMS)进行的, 其特点在于NOx测量设备是基于测量精度较高的非分散紫外线法(NDUV)、非分散红外线法(NDIR)或化学发光法(CLD)等原理开发的, 但测量系统复杂且设备成本较高.国外研究机构和环保主管部门已针对非道路移动机械开展了大量的实测, CAO等(2016; 2018)、Frey等(2008; 2010)和Zavala等(2017)基于PEMS方法研究了不同典型工程机械的排放因子.国内, 付明亮等(2012; 2014; 2013)也开展了工程机械、农业机械的排放因子研究, 叶潇等(2018)则开展了场内叉车的排放因子实测.近年来随着车用氮氧化物传感器的出现, 使得实时尾气NOx监控成为可能.Yang等(2016)的研究, 基于国四排放标准混合动力公交车的OBD通讯接口获得了累计87 h的逐秒NOx排放特性.杨强等(2018)利用车用NOx传感器和车载诊断(OBD)接口获得的车辆信息, 对10辆国4排放标准公交车NOx排放情况进行了为期1年的研究.另外, 胡磬遥等(2015)还利用车用NOx传感器对柴油机燃用生物柴油进行了研究.
我国目前非道路移动机械用发动机排放标准普遍在国0~国2水平, 安装SCR系统的机械数量占比极低, 均不具备通过OBD接口获取NOx排放水平的条件.本研究通过加装在工程机械排气管上的NOx排放测量系统, 建立了非道路移动机械单位油耗NOx排放因子测量方法.测量方法由基于车用NOx传感器的测量系统和NOx单位油耗排放因子算法组成.为了验证NOx排放测量方法, 开展了4台工程机械的实测试验, 并选取其中一台开展了基于车载排放测量系统(PEMS)的排放比对试验.
2 材料与方法(Materials and methods)2.1 测量系统测试及比对系统如图 1所示.NOx排放因子测量系统由车用NOx传感器、环境温湿度传感器和控制器局域网(CAN)通讯解析器、模拟数字量转换模块组成.NOx传感器安装在发动机排气管上, 温度、湿度传感器被安装在被测机械上不会受到尾气干扰的地方.上位机通过CAN通讯解析器获得实时的O2和NOx体积浓度.研究选用的NOx传感器为德国大陆公司生产的5WK9型车用NOx传感器, 传感器同时具有测定排气O2浓度的功能.
图 1(Fig. 1)
 |
| 图 1 测试及比对系统配置 Fig. 1Test and comparison system configuration |
为了验证本研究搭建的测量系统和计算方法, 结合车载排放设备和车用NOx传感器搭建了排放比对测量系统.被测机械发动机产生的尾气经过颗粒物过滤器和保温采样管路进入车载排放测量系统的气体分析单元.尾气测试设备被安装在可移动拖车上, 拖车与被测工程机械刚性连接, 进行比对试验时小车跟随被测机械移动.车载排放系统中的NOx、O2分析模块分别为美国Sensors公司基于非分散红外原理(NDIR)测和电化学原理开发的.NOx传感器和车载排放测量系统的相关测试模块参数如表 1所示.进行试验前车载排放测量系统基于高纯氮气和标准气体进行了3次严格的标零和标定.
表 1(Table 1)
| 表 1 测试及比对系统技术指标 Table 1 Test and comparison system specifications | ||||||||||||||||||||||||
表 1 测试及比对系统技术指标 Table 1 Test and comparison system specifications
| ||||||||||||||||||||||||
 |
2.2 测试机械和工况表 2为测试的4台工程机械的基本信息.测量系统验证试验均按照测试实际工作情形开展.挖掘机的实际工况有怠速工况、工作工况, 而推土机和压路机的实际工况有怠速工况、行驶工况和工作工况(李东玲等, 2012; 付明亮等, 2014).怠速工况:机械启动后和工作循环间歇发动机以700~800 r·min-1的转速维持工作但不输出功.行驶工况:机械依靠发动机提供动力从一个工作地点转移到另一个工作地点, 期间工作机构(铲斗、抓斗等)不工作.工作工况：机械按照实际作业情形在试验场地内反复进行挖掘、推土和平地等作业.
表 2(Table 2)
| 表 2 测试工程机械基本信息 Table 2 Basic information of tested construction machinery | |||||||||||||||||||||||||||||||||||||||||||||
表 2 测试工程机械基本信息 Table 2 Basic information of tested construction machinery
| |||||||||||||||||||||||||||||||||||||||||||||
|
3 排放因子计算模型(Emission factor calculation model)3.1 模型简介模型建立的基本思路是基于进排气氧元素、氢元素平衡, 通过发动机进气状态、燃料碳氢比和排气NOx浓度ΦNOxout(10-6)、O2浓度ΦO2out(%)获得空燃比AFR和单位油耗NOx排放因子.模型推导有3个简化条件：模型忽略CO、THC、NOx和颗粒物排放对排气氧含量的影响；发动机进气中的H2O、CO2、N2和Ar均为理想气体且不参与反应；本模型中燃料的品质是稳定的, 燃料中的碳、氢比例在一次试验中是稳定的.
3.2 发动机进气主要组分空燃比AFR, 式(1)通过空燃比AFR将发动机进气流量mair(kg·h-1)和发动机小时油耗mfuel(kg·h-1)联系起来.为了使得模型的精度更高, 模型包含了H2O、CO2在内的发动机进气中的各种主要气体.如式(2)所示, 通过由环境温、湿度传感器获得的环境温度T0(℃)和环境相对湿度φ可以获得进入发动机的单位质量干空气所携带的水, 即含湿量d(g·kg-1)(盛裴轩, 2013).如式(3)所示, d是以饱和水蒸气分压力ps(hPa)和φ构成的二元函数(盛裴轩, 2013).其中, ps是以环境温度T0为变量的方程(胡磬遥等, 2015).如式(4)、(5)所示, 由含湿量定义可以分别获得发动机进气的干空气和水蒸气的物质量流量.发动机进气中的N2、O2和Ar、CO2物质量流量可由式(6)获得.进气干空气中的N2、O2和Ar、CO2浓度分别为78.084%、20.946%和0.9340%和0.036%.
 | (1) |
 | (2) |
 | (3) |
 | (4) |
 | (5) |
 | (6) |
3.3 发动机排气主要组分根据2.1节中的简化条件, 如式(7)所示, 尾气中的N2和Ar不参与反应, 排气与进气不发生变化.如式(8)所示, 排气中CO2由进气中的CO2加上燃烧生成的CO2组成.燃料中1 kmol·h-1的C元素可以生成1 kmol·h-1的CO2, 将燃料含碳量、AFR和发动机排气中的CO2建立了联系.如式(9)所示, 排气中H2O由进气中的H2O加上燃烧生成的H2O组成.燃料中1 kmol·h-1的H元素可以生成0.5 kmol·h-1的H2O, 如式(9)将燃料含氢量、AFR和发动机排气中的H2O建立了联系.式(10)所示为排气中O2的物质量, 由进气和燃料消耗量之和减去排气中的N2、Ar和CO2、H2O获得.如式(11)所示, 将上述发动机尾气中的各种气体累加, 即可以得到发动机尾气的总物质量流量.
 | (7) |
 | (8) |
 | (9) |
 | (10) |
 | (11) |
3.4 空燃比和单位油耗NOx排放因子由于本文将排气中的各成分近似为理想气体, 所以ΦO2out可近似等于排气中nO2out与nexhaust之比.式(12)是包含了未知ARF和已知的φ、d、ωC、ωH的待解等式.将式(10)~式(11)代入式(12)中对应的位置.而式(12)的右边内含有求解目标ARF, 对式(12)进行求解可以得到式(13)——由ΦO2out和d、ωC、ωH组成的空燃比ARF计算模型.式(14)为瞬态单位燃料NOx排放因子计算方法, 式中的各个量可根据式(11)、(13)求得.最终经过计算EFNOxfuel是由自己机械发动机进气温度T0、湿度φ和排气NOx浓度、O2浓度的函数, 与瞬时发动机进气量无关.需要说明的是, 由于柴油机NOx中NO的比例在90%~95%范围, 所以式中氮氧化物的分子量取31.5 g·mol-1.
当环境温度T0为25 ℃, 相对湿度φ为50%, 燃料的C、H和O含量分别为86%、13.55%、0.45%时.ARF和EFNOxfuel如式(15)和式(16)所示.可以看到空燃比AFR随着排气氧浓度的上升而上升, 单位油耗NOx排放因子EFNOxfuel随排气NOx浓度和排气氧浓度的上升而上升.
 | (12) |
 | (13) |
 | (14) |
 | (15) |
 | (16) |
4 结果讨论(Results and discussion)4.1 排放测试设备比对结果基于挖掘机1开展了怠速和工作工况下, NOx排放因子测量系统和车载排放测量系统的发动机尾气O2浓度及NOx浓度比对试验, 图 2所示为比对测试结果.NOx传感器的响应较快、相对偏差较低, 与车载排放测量系统的测量结果相比具有很好的相关性.
图 2(Fig. 2)
 |
| 图 2 NOx排放因子测量系统与车载排放测量系统比对结果 Fig. 2Test result comparison between NOx emission factor measurement system and portable emission measurement system |
怠速、工作工况下NOx传感器测量得到的O2浓度较车载排放测量系统获得的O2浓度的平均相对偏差分别为-3%和-2%, 此同时NOx浓度的平均相对偏差分别为-6%和2%.两种气体怠速工况的平均相对偏差略大于工作工况.车载排放测量系统与NOx传感器获得的发动机尾气O2浓度及NOx浓度具有很高的相关性, 怠速和工作工况下O2浓度的相关系数r分别为0.94和0.95, 与此同时NOx浓度的相关系数r分别为0.92和0.98.怠速工况的O2、NOx浓度相关系数略低于工作工况.
4.2 排放因子测量方法验证为了验证排放因子测量方法, 开展了4台工程机械的实测研究.图 3所示为测量系统获得的NOx浓度和O2浓度关系.可以看到4台工程机械的排气NOx浓度范围为100×10-6~1300×10-6, 排气O2浓度范围为9%~20%.测试结果表明搭建的NOx排放因子测量系统能够适用于不同排气O2、NOx浓度范围的各类型非道路移动机械.测试结果也反映了被测试机械用发动机的燃烧过程均属于典型的稀薄燃烧, NOx排放水平较高.测量系统能够比较好地反映不同机械尾气中的NOx排放特征, 测试结果显示4台被测工程机械的NOx排放浓度均随着O2浓度的上升而下降, 但不同机械排气NOx浓度随着O2浓度的变化幅度不同.尾气中O2浓度越高表明燃烧反应产物CO2浓度越低, 反映发动机负荷低, 燃烧剧烈程度低, 最终可以看到表征发动机燃烧剧烈程度的NOx浓度下降.
图 3(Fig. 3)
 |
| 图 3 发动机排气氧浓度和氮氧化物浓度关系 Fig. 3Relationship between oxygen concentration and NOx concentration of engine exhaust |
图 4所示为挖掘机1、2和推土机、压路机的瞬态排气O2、NOx浓度和逐秒NOx单位油耗排放因子情况.实测结果表明搭建的NOx排放因子测量方法具有较好的稳定性, 可以用来获得不同类型的非道路移动机械不同工况下的瞬态NOx排放特性和排放因子怠速工况可以看到基于车用NOx传感器的测量方法具有较好的稳定性.怠速工况的300 s内4台工程机械在怠速工况下的NOx、O2浓度的相对误差分别为1.4%±1%、2.2%±1%, 而由NOx、O2浓度获得的逐秒单位油耗NOx排放因子相对误差为2.3%±2%.虽然由于NOx、O2浓度的波动导致NOx排放因子的相对误差高于NOx、O2浓度的相对误差, 但稳定工况下2%左右的相对误差能够满足实测需求.
图 4(Fig. 4)
 |
| 图 4 测试工程机械尾气氧浓度、氮氧化物浓度和瞬态单位油耗NOx排放因子 Fig. 4Exhaust O2 concentration, NOx concentration and transient fuel consumption specific NOx emission factor of tested construction machinery |
由于测量方法中采用的NOx传感器具有较高的时间分辨率, 通过测量方法可以获得不同类型工程机械每个工作周期的瞬态工作和排放特征.挖掘机1、2工作工况、压路机的行走、工作工况下NOx排放因子呈明显的周期性.以工作工况为例, 挖掘机1每个工作循环的周期分别为8.5 s, NOx排放因子在19 g·kg-1和25 g·kg-1之间周期变化, 压路机的周期为27 s, NOx排放因子在5 g·kg-1和120 g·kg-1之间周期变化.而推土机由于工作场地高低不平且土地的硬度没有规律, 所以其行走和工作工况的瞬态排放特性剧烈变化但没有明显的周期性.
由测量方法获得的怠速、行走、工作工况下的瞬态NOx排放因子可以看到：随着柴油机负荷不断上升, 循环喷油量不断上升的同时空燃比下降, 使得发动机缸内燃烧温度、压力均上升, 最终导致热力型NO大量生成, 排气NOx浓度和NOx排放因子不断上升.瞬态NOx排放因子的变化幅度也随着工况负荷的上升而加剧.虽然挖掘机1工作工况下的NOx浓度高于怠速工况, 但由于挖掘机1怠速工况、工作工况下NOx浓度差异小于O2浓度差异(即空燃比差异), 导致挖掘机1怠速工况的NOx单位油耗排放因子高于工作工况.
4.3 NOx排放因子与已有研究的比较为了将单位油耗NOx排放因子测量系统获得的结果与已有研究进行比较, 表 4汇总了被测试工程机械不同工况下平均空燃比AFR、排气O2浓度和NOx浓度和平均单位油耗NOx排放因子的情况.和上述瞬态工况的分析结果一致, 随着怠速、行走和工作工况的负荷不断上升, 4台机械的AFR、排气O2浓度不断下降而NOx浓度上升.与此同时除了挖掘机1以外NOx单位油耗因子也依次上升.挖掘机1、挖掘机2、压路机、推土机的平均排气NOx浓度为(494.8±100)×10-6、(625.0±211)×10-6、(513.6±308)×10-6、(720.1±253)×10-6.单位油耗NOx排放因子分别为(22.8±16)、(32.4±4)、(50.7±13)、(68.8±13) g·kg-1.
表 4(Table 4)
| 表 4 被测工程机械尾气基本情况 Table 4 Exhaust emissions basic characteristic of the measured construction machinery | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
表 4 被测工程机械尾气基本情况 Table 4 Exhaust emissions basic characteristic of the measured construction machinery
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
在广泛调研大量文献的基础上, 本文总结了国内外****得到的各类工程机械在实际运行工况下获得的单位燃料NOx排放因子, 文献调研得到的排放因子结果如表 5所示.可以看到本研究建立测量方法获得的NOx的排放因子与文献调研结果较为接近.挖掘机和国内外实测结果比较一致, 略低于非道路污染物排放清单指南.压路机国内外实测较少, 与Wang等(2016)的研究较为一致, 但远高于清单指南.本研究实测推土机NOx排放因子略高于现行清单指南和国外实测结果.
表 5(Table 5)
| 表 5 国内外研究结果 Table 5 Comparison of domestic and foreign research results | |||||||||||||||||||||||||||||||||||||||||||
表 5 国内外研究结果 Table 5 Comparison of domestic and foreign research results
| |||||||||||||||||||||||||||||||||||||||||||
| |||||||||||||||||||||||||||||||||||||||||||
5 结论(Conclusions)1) 基于车用NOx传感器建立的NOx排放因子测量方法可以高效、低成本地获得非道路移动机械的单位油耗NOx排放因子.
2) 基于环境温度、湿度、燃料特效和排气O2浓度、NOx浓度, 通过氧元素守恒获得了带有环境温、湿度修正的空燃比AFR计算方法和瞬态单位油耗NOx排放因子计算方法.
3) 比对和实测试验结果表明建立的NOx排放因子测量方法具有较低的相对误差和稳定性, 与车载排放测量系统的测量结果相比具有很好的相关性.
4) 包括两台挖掘机、一台推土机和一台压路机在内的4台工程机械的平均瞬态单位油耗NOx排放因子分别为(22.8±16.0)、(32.4±4.0)、(50.7±13.0)和(68.8±13.0) g·kg-1.基于车用氮氧化物传感器NOx排放因子测量方法获得的工程机械排放因子与目前国内外已有实测研究结果接近.
参考文献
| Abolhasani S, Frey H, Kim K, et al. 2008. Real-world in-use activity, fuel use, and emissions for nonroad construction vehicles:a case study for excavators[J]. Journal of the Air & Waste Management Association, 58: 1033–1046. |
| Anenberg S C, Miller J, Injares R M, et al. 2017. Impacts and mitigation of excess diesel-related NOx emissions in 11 major vehicle markets[J]. Nature, 545(7655): 467.DOI:10.1038/nature22086 |
| Cao T, Durbin T D, Russell R L, et al. 2016. Evaluations of in-use emission factors from off-road construction equipment[J]. Atmospheric Environment, 147: 234–245.DOI:10.1016/j.atmosenv.2016.09.042 |
| Cao T, Russell R L, Durbin T D, et al. 2018. Characterization of the emissions impacts of hybrid excavators with a portable emissions measurement system (PEMS)-based methodology[J]. The Science of the Total Environment, 635: 112–119.DOI:10.1016/j.scitotenv.2018.04.011 |
| Frey H C, Kim K, Pang S H, et al. 2008. Characterization of real-world activity, fuel use, and emissions for selected motor graders fueled with petroleum diesel and B20 biodiesel[J]. Journal of the Air & Waste Management Association, 58(10): 1274–1287. |
| Frey H C, Rasdorf W, Lewis P. 2010. Comprehensive field study of fuel use and emissions of nonroad diesel construction equipment[J]. Transportation Research Record, 2158(1): 69–76.DOI:10.3141/2158-09 |
| Fu M, Ge Y, Tan J, et al. 2012. Characteristics of typical non-road machinery emissions in China by using portable emission measurement system[J]. Science of the Total Environment, 437: 255–261.DOI:10.1016/j.scitotenv.2012.07.095 |
| 付明亮, 丁焰, 尹航, 等. 2013. 实际作业工况下农用拖拉机的排放特性[J]. 农业工程学报, 2013, 29(6): 42–48. |
| 付明亮, 丁焰, 尹航, 等. 2014. 基于燃油消耗的工程机械排放因子研究[J]. 北京理工大学学报, 2014, 34(2): 138–142. |
| 胡磬遥, 楼狄明, 胡志远, 等. 2015. 基于排气氧浓度的车用生物柴油混合比例判断策略[J]. 内燃机工程, 2015, 36(6): 14–19. |
| 黄成, 安静宇, 鲁君. 2018. 长三角区域非道路移动机械排放清单及预测[J]. 环境科学, 2018, 39(9): 1–13.DOI:10.3969/j.issn.1673-1212.2018.09.001 |
| Heidari B, Marr L C. 2015. Real-time emissions from construction equipment compared with model predictions[J]. Journal of the Air & Waste Management Association, 65(2): 115–125. |
| 李东玲, 吴烨, 周昱, 等. 2012. 我国典型工程机械燃油消耗量及排放清单研究[J]. 环境科学, 2012, 33(2): 518–524. |
| 鲁君, 黄成, 胡磬遥, 等. 2017. 长三角地区典型城市非道路移动机械大气污染物排放清单[J]. 环境科学, 2017, 38(7): 2738–2746. |
| 盛裴轩, 毛节泰, 李建国. 2013. 大气物理学[M]. 北京: 北京大学出版社. |
| Wang F, Li Z, Zhang K, et al. 2016. An overview of non-road equipment emissions in China[J]. Atmospheric Environment, 132: 283–289.DOI:10.1016/j.atmosenv.2016.02.046 |
| Yang L, Zhang S, Wu Y, et al. 2016. Evaluating real-world CO2 and NOx emissions for public transit buses using a remote wireless on-board diagnostic (OBD) approach[J]. Environmental Pollution, 218: 453–462.DOI:10.1016/j.envpol.2016.07.025 |
| 杨强, 胡磬遥, 黄成, 等. 2018. 基于远程通讯技术的混动公交车SCR系统运行及NOx排放特征[J]. 环境科学, 2018, 39(10): 1–15.DOI:10.3969/j.issn.1673-1212.2018.10.001 |
| 叶潇, 李成, 夏泽群, 等. 2018. 内燃叉车尾气气态及细颗粒污染物排放特征研究[J]. 环境科学学报, 2018, 38(4): 1392–1403. |
| Zavala M, Huertas J I, Prato D, et al. 2017. Real-world emissions of in-use off-road vehicles in Mexico[J]. Journal of the Air & Waste Management Association, 67(9): 958–972. |
| Zhao B, Wang S X, Liu H, et al. 2013. NOx emissions in China:historical trends and future perspectives[J]. Atmospheric Chemistry and Physics, 13(19): 9869–9897.DOI:10.5194/acp-13-9869-2013 |
| 张楚莹, 王书肖, 邢佳, 等. 2008. 中国能源相关的氮氧化物排放现状与发展趋势分析[J]. 环境科学学报, 2008, 28(12): 2470–2479.DOI:10.3321/j.issn:0253-2468.2008.12.011 |
| 张礼俊, 郑君瑜, 尹沙沙, 等. 2010. 珠江三角洲非道路移动源排放清单开发[J]. 环境科学, 2010, 31(4): 886–891. |
| 中华人民共和国环境保护部.2015.非道路移动源大气污染物排放清单编制技术指南(试行)[S].北京 |
