删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

内燃叉车尾气气态及细颗粒污染物排放特征研究

本站小编 Free考研考试/2021-12-31

叶潇1, 李成2, 夏泽群1, 王延龙1, 卞雅慧1, 肖笑1, 廖松地1, 郑君瑜1
1. 华南理工大学环境与能源学院, 广州 510006;
2. 暨南大学环境与气候研究院, 广州 511486
收稿日期: 2017-10-24; 修回日期: 2017-12-15; 录用日期: 2017-12-15
基金项目: 国家自然科学基金项目(No.41627809);国家****科学基金项目(No.41325020);广东省省级科技计划项目(No.2016B050502005)
作者简介: 叶潇(1992-), 女, E-mail:yyexiao@foxmail.com
通讯作者(责任作者): 郑君瑜, E-mail:zhengjunyu_work@hotmail.com

摘要: 内燃叉车在我国保有量大、分布广,其尾气排放严重危害空气质量和人体健康.本文以7台典型内燃叉车为研究对象,搭建工程机械尾气测试系统,识别了实时的"燃油-机械特征-污染排放"关系,评估了不同发动机功率、排放标准、工况对排放因子的影响,利用自展模拟的方法量化排放因子的不确定性.结果表明:内燃叉车NO和NO2的单位燃油的排放速率变化趋势相同,NO与NO2的峰值出现时间略有差距.CO、NO、NO2、NOx、PM2.5和SO2的综合排放因子与不确定性范围分别是24.1 g·kg-1(-25.4%,26.1%)、26.5 g·kg-1(-24.2%,24.4%)、10.2 g·kg-1(-22.3%,21.9%)、36.7 g·kg-1(-21.5%,22.6%)、2.8 g·kg-1(-25.5%,24.7%)和0.9 g·kg-1(-55.7%,80.6%).同时对作业工况下PM2.5的组分分析发现,TC为主要成分,质量占比达77%~92%;OC/EC的平均比值为4.4;水溶性离子以NO3-和PO43-为主;金属组分中Ca、Na、Mg和Fe占比较多.
关键词:内燃叉车排放特征气态污染物排放因子细颗粒物
Characteristics of gaseous pollutants and fine particulates from diesel forklifts
YE Xiao1, LI Cheng2, XIA Zequn1, WANG Yanlong1, BIAN Yahui1, XIAO Xiao1, LIAO Songdi1, ZHENG Junyu1
1. School of Environment and Energy, South China University of Technology, Guangzhou 510006;
2. Institute for Environmental and Climate Research, Jinan University, Guangzhou 511486
Received 24 October 2017; received in revised from 15 December 2017; accepted 15 December 2017
Supported by the National Natural Science Foundation of China(No. 41627809), the National Science Fund for Distinguished Young Scholars(No. 41325020) and the Key Technology Research and Development Program of Guangdong Province, China(No. 2016B050502005)
Biography: YE Xiao (1992—), female, E-mail: yyexiao@foxmail.com
*Corresponding author: ZHENG Junyu, E-mail:zhengjunyu_work@hotmail.com
Abstract: Due to the large number, emissions from diesel forklifts in China are high, posing an important impact on human health and air quality. In this study, seven typical forklifts were selected to characterize the emission pattern of gaseous pollutants and fine particulates by using a set of on-board test system. This study not only identified the real-time "fuels-machines-emissions" relationship, but also evaluated the influences of different engine power, emission standard and operation models on emission factors. Monte Carlo Method was used to quantify uncertainties in emission factors. Results showed that nitric oxide (NO) and nitrogen dioxide (NO2) have a similar emission trend, but their peak periods were not identical. The average emission factors of carbon monoxide (CO), NO, NO2, nitrogen oxides (NOx), particulate matter (PM2.5) and sulfur dioxide (SO2) were 24.1 g·kg-1 (-25.4%, 26.1%), 26.5 g·kg-1 (-24.2%, 24.4%), 10.2 g·kg-1 (-22.3%, 21.9%), 36.7 g·kg-1 (-21.5%, 22.6%), 2.8 g·kg-1 (-25.5%, 24.7%) and 0.9 g·kg-1 (-55.7%, 80.6%) respectively. Additionally, this study obtained PM2.5 compositions in the working conditions. The total carbon (TC) was the dominant component in the exhaust PM2. 5 emitted from forklifts, accounting for 77%~93% of the total PM2.5 emissions. The ratio of organic carbon (OC) to elemental carbon (EC) was 4.4, and NO3- and PO43- were the most important contributor to water soluble ion. Moreover, Ca、Mg、Na and Fe were the main contributors in inorganic elements.
Key words: diesel-fueled forkliftemission characteristicsgaseous pollutantsemission factorPM2.5
1 引言(Introduction)近年来, 我国大气正面临严峻的区域性复合污染问题, 工程机械已成为国内诸多区域NOx和PM的重要污染排放贡献源(张礼俊等, 2010Yan et al., 2014鲁君等, 2017).据《中国工程机械工业统计年鉴》统计, 截止到2015年, 我国内燃叉车保有量多达120万台, 占工程机械比例15%以上, 且近5年年均增长率超过10% (中国机械工业联合会, 2010—2015).内燃叉车在工业生产、物流搬运、工程作业等过程中排放的污染物不但对局部空气质量影响较大, 也严重危害人体健康.
污染源排放特征识别是精细化环境管理的基础, 开展排放因子和排放源成分谱特征研究对于识别污染来源、提升空气质量管理有重要作用.近些年来, 国内外****基于台架实验(Lindgren et al., 2011)、模型模拟(Abolhasani et al., 2008)、现场测试(Abolhasani et al., 2008Frey et al., 2008Fu et al., 2012)等方法对主要类型的工程机械尾气排放特性开展了一系列工作.台架实验方便操作, 能够识别影响柴油机排放因子的主要因素, Lindgren等利用台架实验研究了不同燃料在不同工况条件下, 常规气态污染物(Lindgren et al., 2010)和VOCs组分(Lindgren et al., 2011)的排放特征.然而, 台架试验难以反映机械的实际排放特征, 现场测试方法更能准确探究真实作业环境下机械、燃油、工况等条件对污染物排放特征的影响(Abolhasani et al., 2008Lewis, 2009).随着源排放测试研究的深入, 现场测试成为工程机械污染物排放特征研究的重要方法.例如, 北卡罗来纳州立大学研究团队基于工程机械排放源现场测试的方法, 测试了多台平地机、反铲机、挖掘机、推土机、装载机等机械, 识别了发动机参数、实际工况、进气压变化等因素对燃油消耗率和尾气中NOx、HC、CO等污染物的排放特征的影响关系(Abolhasani et al., 2008Frey et al., 2008, 2010Lewis, 2009;Phil et al., 2011).崔敏(2016)在对挖掘机PM2.5排放因子进行现场测试的基础上, 进一步分析了组成特征.但是现场测试容易受到测试环境的限制, 仪器难以安装和稳定, 测试的不确定性较大.不仅如此, 现阶段已建立的工程机械排放源研究多集中在挖掘机、装载机等类型的机械中, 内燃叉车作为保有量最大的一类工程机械, 基于实测的排放特征研究仍鲜有报道(Durbin et al., 2007朱倩茹等, 2015Wang et al., 2016).同其他工程机械相比, 内燃叉车的自重较小, 设计构造和功能也有显著差异, 基于此, 亟需开展基于实测的内燃叉车排放特征研究.
本文选取典型内燃叉车作为研究对象, 搭建现场尾气测试系统, 对不同工况下的尾气气态污染物的实时排放量进行测试研究, 分析了影响排放因子的关键因素, 建立CO、NO、NO2、NOx、PM2.5、SO2的排放因子, 采用自展模拟的方法定量评估排放因子的不确定性.同时, 深入分析尾气样品中的PM2.5的组成特征, 建立内燃叉车的PM2.5成分特征谱.本研究旨在全面了解国内内燃叉车的污染物排放特征, 为非道路移动源清单开发及污染源来源解析提供基础数据, 为非道路移动源的控制策略提供重要依据和参考.
2 材料与方法(Materials and methods)2.1 测试系统搭建与颗粒物样品分析现场测试搭建的烟气采集系统主要由气体采集与分析模块、流量监测模块、PM2.5采样模块3个部分组成.
气体采集与分析模块采用奥地利Madur品牌Photon便携式烟气分析仪, 测试实时的尾气气态污染物浓度.烟气分析仪单独配置烟气预处理系PGD-100, 用于干燥采样气体并移除采样气体中的颗粒物, 加热采样管线能够保持采集气体的温度恒定在120 ℃.烟气分析仪不仅可以自动测量尾气管道中CO、CO2、NO、NO2、SO2的浓度, 还能同时监测环境温度、烟气温度、大气压力等参数.烟气分析仪在线测试组分的具体参数如表 1所示.
表 1(Table 1)
表 1 烟气分析仪在线测试组分参数表 Table 1 Parameters in online component of handy flue gas analyzer
表 1 烟气分析仪在线测试组分参数表 Table 1 Parameters in online component of handy flue gas analyzer
测试指标 测试原理 检测限 分辨率 精度
CO 非分光红外(NDIR) 1~10000 ppm 1 ppm 读数的2%
CO2 非分光红外(NDIR) 0.01%~25% 0.01% 读数的2%
NO 非分光红外(NDIR) 1~5000 ppm 0.2 ppm 读数的2%
NO2 非分光红外(NDIR) 1~1000 ppm 0.2 ppm 读数的2%
O2 顺磁法 0.1%~25% 0.01% 读数的2%
SO2 非分光红外(NDIR) 1~1000 ppm 1 ppm 读数的2%


流量监测模块采用耐高温高湿实时便携式2445型Kurz流量计及自制不锈钢烟道管, 以满足该流量计测试对烟气稳流状态的要求, 流速测量的时间间隔为2 s.
PM2.5采样模块采用多功能便携式稀释通道PM2.5采样器(MFD25)采集叉车尾气的PM2.5样品. MFD25采样器由采样管、稀释系统、中控系统、供气系统和抽气系统组成.稀释通道主要由对流折叠的停留室构成, 便于尾气和洁净空气形成湍流状态, 易于混合, 再经PM2.5切割器切割并分流, 由石英滤膜和特氟龙膜采集得到颗粒物样品.把采集到的PM2.5滤膜先放入恒温恒湿箱内以温度25 ℃、湿度41%的条件恒温恒湿24 h, 再称量采样前后滤膜的质量确定样品质量.针对PM2.5的碳质组分、水溶性离子、金属组分分别利用热光反射法\热光透射法、离子色谱法分析、电感耦合等离子体质谱法测试分析.
2.2 对象选取及油品测试为全面了解内燃叉车的尾气排放特征, 采用正交的方法选取了涵盖不同国标、净功率、载重量的典型内燃叉车7台进行现场测试, 测试对象基本信息如表 2所示.每辆叉车测试典型工况包括怠速、行走(空载)、作业(空载) 3个典型工况, 每个工况测试时间不少于15 min.
表 2(Table 2)
表 2 测试对象基本信息 Table 2 Details of tested forklifts
表 2 测试对象基本信息 Table 2 Details of tested forklifts
编号 排放标准 厂商 型号 出厂年份 吨位/kg 净功率/kW 载重量/kg 累计工作时间/h
1 国一 杭叉 CPC30 2010.08 4300 36.8 3000 3772
2 国一 杭叉 CPCD40 2010.05 6300 57 4000 1473.5
3 国二 杭叉 CPCD100 2013 12895 81 10000 2400+
4 国二 杭叉 CPCD70 2012.07 9300 81 7000 2356
5 国二 杭叉 CPCD50 2012.03 8000 57 5000 2582.7
6 国二 合力 CPCD70 2013.11 9350 85 7000 794.5
7 国二 合力 CPC35 2015.12 4700 40 3500 469+


叉车的燃油均为广州地区中国石油公司标准的92号柴油, 为了检验实际使用的柴油品质, 保障碳平衡法计算排放因子时的数据可靠性, 分别对每次测试的内燃叉车的柴油品质进行检测, 柴油测试方法及测试结果分别如表 3表 4所示. 表 4中A代表编号为3号的叉车燃油, B代表 1、4、6、7号叉车的燃油, C代表 2、5号叉车的柴油.相比于《非道路移动机械用柴油机排气污染物排放限值及测量方法》 (GB252—2005)规定的油品质量标准, 有2台叉车的柴油含硫率偏高, 不符合国三油质标准, 其余5台均符合国三油质标准.
表 3(Table 3)
表 3 油品检测方法 Table 3 Test methods of fuels
表 3 油品检测方法 Table 3 Test methods of fuels
测试指标 测试原理/方法 仪器
含硫率 GB/T17040—2008 X-荧光测硫仪Lab-X3500
灰分 GB/T508—1985 (2004) 马弗炉
水分 GB/T260—1977 (2004) 831KF 832KF瑞士万通卡氏水分测定仪
运动粘度 GB/T265—1988 (2004) PSL TCB-7Mkll恒温浴
碳质含量 GB/T384计算 5E-1AC/ML自动量热仪
金属组分 ASTM D5185-13e1 SPECTRO CIROS VISION电感耦合等离子体发射光谱仪



表 4(Table 4)
表 4 燃油样品分析结果 Table 4 Specifications of fuels used in emission testing
表 4 燃油样品分析结果 Table 4 Specifications of fuels used in emission testing
编号 灰分 水分 运动粘度@20 ℃/
(mm2·s-1)
C S 质量分数/(mg·kg-1)
A 0.006% 痕迹 5.501 86.00% <0.015% <1 43 <1 <1 <1 <1 <1 <1 <1
B 0.004% 痕迹 5.132 86.04% <0.015% <1 13 <1 <1 <1 <1 <1 <1 <1
C 0.005% 痕迹 5.594 87.33% 0.108% <1 <1 <1 <1 <1 <1 <1 <1 <1


2.3 数据分析方法2.3.1 排放因子本研究中, 叉车的气态尾气及PM2.5的排放因子利用碳平衡法计算.假设燃料燃烧后, 所有的碳以OC、EC、CO和CO2的形式排出.公式如下:
(1)
式中, EFi为污染物i的排放因子(g·kg-1);Ciw为污染物i的质量浓度(g·m-3);RFG为单位质量燃油气体排放量(m3·kg-1).
(2)
式中, cF为单位质量燃料中含碳量(g·kg-1) (实验测量值); cCOcCO2cECcOC代表尾气排放的CO、CO2、EC、OC中碳的质量浓度(g·m-3).
综合排放因子计算公式如式(3)所示:
(3)
其中, EFj, i表示第j台叉车污染物i的综合排放因子(g·kg-1);EFj, i, k表示第j台叉车在k工况下污染物i的排放因子;Pi, k表示为一般情况下, 第j台叉车k工况的时间权重.时间权重以30台内燃叉车的现场调研及问卷调查为基础, 获取每台内燃叉车在实际运行条件下不同工况的时长, 计算得出平均状态下3种工况的时间权重分别为怠速工况5%, 行走工况55%, 作业工况40%.
2.3.2 不确定性分析方法常规研究中均默认测试结果应该满足正态分布, 所以偏差被用来作为测试结果不确定性的表达方式, 由于测试的样品数量有限, 数值模拟是一种比概率统计方法更灵活的不确定性分析方法(Zheng, 2012;Zheng and Frey, 2004).采用自展模拟的方式, 能够在样本量较少的情况下, 拟合出最佳的分布模型, 通过自展抽样与统计分析, 得到可变性和不确定性的大小.自展模拟的方法便于准确表达测试的不确定性范围, 有效评估所选影响因素对排放因子测试的合理性, 为测试的质量保证和质量控制提供评价指标(郑君瑜, 2014).
利用不确定性分析软件Analysis of Variability and Uncertainty Tool Progress (AuvToolPro)定量分析排放因子的不确定性, AuvToolPro内设计了7种概率分布模型、2种参数估值方法和2种检验方法.概率分布类型包括正态分布(Normal)、对数正态分布(Lognormal)、贝塔分布(Beta)、伽马分布(Gamma)、韦伯分布(Weibul)、均匀分布(Uniform)和对称三角分布(Symmetric triangle); 参数估计方法分别是矩匹配法(Method Of Moment)和最大似然估算法(Max Likelihood Method); 检验方法包括Kolmogorov-Smirnov (K-S)检验和Anderson-Darling (A-D)检验, 对于输入的变量, AuvToolPro内设计了自展模拟的模块, 能够自行抽样和模拟运算(郑君瑜, 2014).自展模拟分析不确定性的方法如下:①模型拟合排放因子的概率分布以Kolmogorov-Smirnov检验(K-S检验)为主要判据, K-S值越小表示选择的概率分布类型越能更好的反映总体的分布特征, 以K-S检验值最小的分布类型作为选定的概率分布类型;②根据选取的分布类型, 用自展模拟的方式模拟数据的随机抽样过程, 抽样次数不少于1000次, 量化测试结果的不确定性, 以样本容量的第2.5%和97.5%位数作为最低值和最高值;③分析自展模拟的结果, 描述排放因子的不确定性.不确定性的范围上、下限的计算如公式(4)、(5)所示:
(4)
(5)
其中, 最低值、样本均值、最高值是统计量95%置信区间的最低值、样本均值和最高值.
3 结果与讨论(Results and discussion)3.1 气态污染物瞬时排放特征7台内燃叉车在3种工况下的气态污染物瞬时排放特征, 均呈现出怠速时较为平稳, 行走和作业工况下波动较大的特点.从整体分布来看怠速工况下, 7台叉车CO2的排放速率波动最小, 标准偏差均低于0.020 g·s-1, 行走和作业工况下CO2的排放速率的波动较大, 标准偏差在0.12~0.50 g·s-1之间. 图 1是以1号叉车为例的内燃叉车在不同工况下的气态污染物瞬时排放特征及燃油的瞬时消耗速率图, 包括CO2、CO、NO、NO2的实时排放速率、燃油燃烧的实时消耗量、每单位燃油燃烧产生CO2、CO、NO、NO2的实时排放量, 即单位燃油的排放速率.由图可知, CO、NO、NO2的单位燃油的排放速率在怠速时波动最小, 例如CO在怠速时的单位燃油的排放速率范围是39.6~43.4 g·kg-1·s-1, 作业工况下波动达24.2~37.8 g·kg-1·s-1.由于叉车在怠速时发动机无负荷运转, 发动机转速稳定, 因此油耗率低且各类污染物排放平稳;在作业工况时需上下移动前叉, 导致叉车负荷不断变化, 发动机转速和扭矩发生瞬间变化(Lindgren and Hansson, 2004), 消耗的燃油速率也随之增加且剧烈波动, 引起各类污染物瞬时排放量出现较多峰值.在挖掘机、装载机的瞬时排放特征研究中, 也发现各类污染物排放速率在怠速时波动幅度最小, 作业时波动幅度最大, 同本研究的结论类似(Abolhasani et al., 2008;付明亮等, 2014;Fu et al., 2012).
图 1(Fig. 1)
图 1 3种工况下尾气排放及燃油消耗率的瞬时变化图 Fig. 1Temporal variation of fuel consumption and pollutant emission characteristics in three operating modes

图 1中作业工况下的局部放大得到图 2.从气态污染物和燃油消耗速率的瞬时变化一致性上观察, 可以发现在任一时刻, 燃油消耗速率的波动同CO、NO、NO2的排放速率变化基本同步, 由于CO来源于燃料中的含碳部分的燃烧, 而NO和NO2则来源于燃烧时空气中的氧气与氮气的反应, 因而瞬时的燃油消耗速率直接影响尾气的排放速率.
图 2(Fig. 2)
图 2 作业工况下尾气排放及燃油消耗率的局部放大(a.单位燃油的排放速率;b.排放速率及燃油消耗率) Fig. 2Temporal variation of fuel consumption and pollutant emission characteristics in working conditions (a. fuel-based pollutant emission rates; b. pollutant emission rates and fuel consumption rates)

对比CO2、CO、NO、NO2的单位燃油的排放速率波动幅度可以看出, CO、NO和NO2的瞬时变化比较明显, CO2的每单位燃油燃烧产生量的波动最小, 因为燃料中99.8%的碳转化为了CO2 (Abolhasani et al., 2008), 所以CO2的整体排放情况比较稳定, 而CO、NO、NO2、NOx的排放则更容易受到瞬时工况变化、温度变化、油气混合程度等因素的影响. CO2与CO的单位燃油的排放速率的变化曲线刚好相反, CO2的峰值对应于CO的谷值.当某时刻CO2的单位燃油产生量高时, 表明燃料充分燃烧, 由燃烧过程中碳元素守恒(PM中的碳占比较少对整体的影响较小), 导致CO的单位燃油产生量低. NO和NO2的单位燃油的排放速率的变化趋势基本相同, 但是峰值出现时间不完全相同, 说明在作业工况下燃烧室中的温度较高, 存在部分NO2与NO的相互转化.
3.2 气态污染物综合排放特征3.2.1 单位时间排放特征为消除工况的影响, 研究不同净功率下内燃叉车尾气气态的小时排放量, 获取CO、NO、NO2和NOx 3种工况下的小时排放量, 再根据调研获得的3种工况的时间权重分别计算整体情况下单位小时中内燃叉车尾气排放的CO、NO、NO2和NOx的排放量如图 3所示.总体而言, 当发动机的净功率增大时, CO、NO、NO2和NOx的小时排放量增加.当37 kW<P<75 kW时, CO、NO、NO2、NOx的小时排放量分别是P≤37 kW时的1.2倍、5.0倍、2.5倍、4.2倍.更大的发动机需要更多的空气和燃料燃烧, 来满足更大的负荷要求, 因此燃油消耗率更大, 相应地, CO、NO、NO2、NOx的小时排放量更高.该结论在其他****研究中也有发现(Fu et al., 2012).当P≥75 kW时, NO和NOx的小时排放量较37 kW<P<75 kW时有所下降, 在发动机功率大于75 kW时, 会采用高压共轨技术(刘旭刚等, 2007)等方式减少NO和NOx的排放, 以满足国家标准的污染物限值.
图 3(Fig. 3)
图 3 不同功率下叉车气态污染物小时排放量 Fig. 3Hourly emissions of gaseous pollutants from forklifts with different engine power

3.2.2 单位燃油消耗的排放特征为了消除发动机的功率对排放因子的影响, 用基于燃料的排放因子表征各污染物的排放水平(Frey et al., 2010). 图 4表示不同工况和不同排放标准对各排放因子的影响.
图 4(Fig. 4)
图 4 不同工况(a)与不同排放标准(b)下的排放因子 Fig. 4Comparison of fuel-based emission factors on different operating models (a) and emission standards (b)

图 4为不同工况下的排放因子对比图, CO、NO、NOx的排放因子均呈现怠速>行走>作业的特征. CO、NO、NOx在怠速时的排放因子分别是作业时的1.1、1.2、1.2倍, SO2怠速时的排放因子为行走时的2.6倍.这是由于怠速时内燃叉车燃烧室内的燃油与空气混合不均匀, 导致柴油燃烧不充分, CO的排放因子较高.与此同时, 燃烧室内会产生较多的低温区, 使得NO2在怠速时排放量高于其他工况. NO的排放因子在作业工况下最低, 这是由于作业工况下内燃叉车的负荷较大, 进气量较多, 燃油燃烧更加充分, 因此基于燃油的排放因子在作业工况下最低.而尾气排放的NOx主要由NO和NO2组成, 因而NOx的排放因子也在怠速时最高.
本研究测试的内燃叉车中有两台符合国一标准(简称国一叉车); 五台叉车符合国二标准(简称国二叉车).由图 4b可以看出, 随着非道路机械排放标准的提高, CO、NO、NOx、SO2的排放因子均有不同程度的降低.国一叉车CO、NO、NOx、SO2的平均排放因子分别是25.6 g·kg-1、44.0 g·kg-1、57.4 g·kg-1、1.6 g·kg-1, 国二叉车的排放因子则分别下降到了23.6 g·kg-1、27.6 g·kg-1、37.1 g·kg-1、0.7 g·kg-1.从国一排放标准到国二排放标准, CO、NO、NOx、SO2的排放因子分别降低了7.7%、37.2%、35.4%和53.4%.新排放标准的实施促进了燃烧技术的推广, 为了满足更高的排放标准, 国二阶段的发动机采用增压中冷技术, 通过增加进气量、提高燃烧效率来减少CO的排放;在增压的同时, 利用中冷器降低进气温度以减少NOx的产生(魏名山等, 2002潘跃和吴子文, 2008).
3.2.3 综合排放因子及不确定范围由公式(3)分别得到7台内燃叉车CO、NO、NO2、NOx、SO2的综合排放因子, 再用自展模拟对CO、NO、NO2、NOx、SO2的排放因子进行不确定性分析. 7台叉车的综合排放因子及不确定性范围如表 6所示, 其中, SO2的排放因子不确定性范围最大, 原因是6、7号叉车的柴油中含硫率超过国家标准, 燃油中的硫随着燃烧会以SO2的形式排放, 因而这两台叉车SO2的排放因子偏高, 导致整体SO2排放因子不确定性稍高. CO、NO、NOx的排放因子不确定性的上下波动范围均在20%~30%之间, 这部分的不确定性可能来源于发动机功率、吨位、车龄、测试环境等因素.
表 5(Table 5)
表 5 各叉车的平均排放因子分布 Table 5 Emission factors and uncertainty of forklifts
表 5 各叉车的平均排放因子分布 Table 5 Emission factors and uncertainty of forklifts
g·kg-1
编号 CO NO NO2 NOx SO2
1 22.67 16.28 4.85 21.13 0.12
2 30.16 30.17 9.97 40.14 0.26
3 37.17 15.31 11.56 26.88 0.98
4 28.64 36.27 11.60 47.87 0.33
5 14.29 27.62 11.17 38.79 0.71
6 13.54 37.13 14.40 51.53 1.97
7 22.42 22.87 8.01 30.88 1.71



表 6(Table 6)
表 6 综合排放因子及不确定性范围 Table 6 Emission factors and uncertainty of forklifts
表 6 综合排放因子及不确定性范围 Table 6 Emission factors and uncertainty of forklifts
污染物 平均值/
(g·kg-1)
分布类型 95%置信区间/
(g·kg-1)
不确定性范围
CO 24.1 Weibull (18.1, 30.7) (-25.4%, 26.1%)
NO 26.7 Weibull (20.2, 33.2) (-24.2%, 24.4%)
NO2 10.2 Normal (8.0, 12.5) (-22.3%, 21.9%)
NOx 37.0 Weibull (28.8, 45.0) (-21.5%, 22.6%)
SO2 0.9 Weibull (0.4, 1.6) (-55.7%, 80.6%)


叉车综合排放因子的比较如表 7所示, 由于机械类型众多, 且缺乏统一标准的现场测试, 各研究的工程机械、港口机械的平均排放因子差异明显, 具有较大的不确定性.本研究中内燃叉车CO、NOx和SO2的排放因子均处于工程机械类综合排放因子的低值区, 可能是因为相比于挖土机、推土机等机械, 叉车的自重较小, 排放因子偏低.同NONROAD模拟出的叉车排放因子(朱倩茹等, 2015)相比, 本研究的SO2的排放因子与模拟结果较为接近, CO和NOx的实测排放因子高于模型模拟的结果.
表 7(Table 7)
表 7 叉车尾气污染物排放因子比较 Table 7 Fuel-used emission factors for forklift
表 7 叉车尾气污染物排放因子比较 Table 7 Fuel-used emission factors for forklift
排放因子/(g·kg-1) 机械类型 来源
CO NOx NO SO2
24.10 37.00 26.7 0.9 叉车 本研究
17.70 29.70 N/A 0.8 叉车 朱倩茹等, 2015
10.72 32.79 N/A N/A 工程机械 鲍晓峰等, 2014
31.38 49.76 N/A N/A 工程机械 殷小鸽, 2010
41.56 66.10 N/A N/A 工程机械 张礼俊等, 2010
N/A 52.60 N/A N/A 工程机械 李东玲等, 2012
42.65 51.62 N/A N/A 港口机械 贾旭等, 2014
14.50 50.30 N/A N/A 港口机械 谭华等, 2013


3.3 细颗粒物排放因子与组分特征3.3.1 排放因子不同工况及不同国标下PM2.5的分布特征如图 5所示, 国二叉车PM2.5的排放因子较国一叉车明显降低, 在行走和作业工况下降幅度达19%和22%.由于国二阶段采用的高压燃油喷射等技术能够有效提高燃油雾化质量, 降低PM2.5的排放.对比不同工况下的排放特征, 内燃叉车PM2.5在行走和作业工况下的排放因子高于怠速.这是因为PM2.5的生成同机械的负荷相关(Liu et al., 2005), 内燃叉车行走和工作时负荷较大, 加速了PM2.5的生成.
图 5(Fig. 5)
图 5 不同国标和不同工况下PM2.5的排放因子 Fig. 5Comparison of fuel-based emission factors of PM2.5 on different operating models and emission standards

选取6台叉车在作业工况下采集的PM2.5样品, 分析其中的碳质组分、水溶性离子、金属离子的含量. 表 8表示每台叉车中的OC、EC、水溶性离子、金属组分及PM2.5的排放因子, 图 6给出了每台叉车的OC与EC的质量百分比及OC/EC的比值.将测得的组分浓度归一化处理后, 得到各组分的质量百分比, 建立叉车PM2.5源排放特征谱如表 9所示.本研究对PM2.5的研究发现, PM2.5排放因子的范围在1.71~3.76 g·kg-1, 平均值为2.8 g·kg-1, 蒙特卡罗分析得到不确定性范围在(-25.5%, 24.7%);对PM2.5的组分排放因子分析得出OC、EC、水溶性离子、金属组分的排放因子范围为别为2319~3983 mg·kg-1、465~3639 mg·kg-1、49~451 mg·kg-1和17~256 mg·kg-1.
表 8(Table 8)
表 8 PM2.5中各组分的排放因子 Table 8 Emission inventories of chemical composition in PM2.5
表 8 PM2.5中各组分的排放因子 Table 8 Emission inventories of chemical composition in PM2.5
组分 排放因子/(mg·kg-1)
1号 2号 3号 5号 6号 7号
碳质组分 OC 3748.04 3235.18 3509.95 2318.94 4739.70 3982.51
EC 638.10 799.45 465.07 3638.81 520.91 1287.75
F- nd 1.07 1.74 4.23 1.71 nd
Cl- 2.75 nd 1.31 nd 2.78 6.46
NO2- nd 8.58 3.78 28.62 11.31 nd
NO3- 15.68 89.86 22.82 7.15 56.14 176.78
水溶性离子 PO43- 2.99 56.87 16.86 48.78 45.47 53.63
SO42- 1.56 59.28 1.45 31.87 5.55 18.87
Na+ 15.44 18.78 nd 6.18 5.98 128.61
NH4+ nd 3.22 nd nd 1.07 nd
Mg2+ 0.48 1.61 nd nd 0.64 5.46
Ca2+ 10.53 16.90 8.72 11.38 21.77 61.08
Na 3.47 11.95 5.78 8.61 10.29 31.93
Mg 3.28 10.05 5.25 8.26 9.29 33.95
K nd nd nd 0.12 nd 5.91
Ca 3.90 3.91 3.71 5.49 3.19 130.69
V 0.00 0.01 0.00 0.01 0.00 0.08
Cr 0.20 0.40 0.34 0.55 0.20 2.73
Mn 0.17 0.01 0.02 0.03 0.03 0.82
Fe 0.62 1.04 1.12 1.55 0.43 41.00
Co 0.00 nd nd 0.00 nd 0.02
金属组分 Ni 1.76 0.00 0.03 0.04 0.00 0.84
Cu 12.75 0.04 0.03 0.10 0.03 0.73
Zn 15.66 0.13 0.73 1.28 0.09 5.59
Ga 0.00 0.01 0.00 0.00 0.00 0.08
As 0.00 0.00 0.00 0.01 0.00 0.01
Rb 0.00 0.00 0.00 0.00 0.00 0.04
Sr 0.01 0.01 0.01 0.01 0.01 1.11
Ag 0.00 0.00 0.00 0.00 0.00 0.00
Cd 0.01 0.00 0.00 0.00 nd 0.00
Cs 0.00 0.00 0.00 0.00 0.00 0.00
Ba 0.01 0.06 0.03 0.02 0.01 0.72
Pb 0.04 0.01 0.00 0.02 0.00 0.12
颗粒物 PM2.5 3730 2840 3760 1710 2700 3350
注:nd表示低于检出限.



图 6(Fig. 6)
图 6 作业工况下OC与EC分布图 Fig. 6The distribution of OC and EC in PM under working condition


表 9(Table 9)
表 9 金属元素的相关性分析 Table 9 The elements correlation analysis
表 9 金属元素的相关性分析 Table 9 The elements correlation analysis
Na Mg K Ca V Cr Mn Fe Co Ni Cu Zn Ga As Rb Sr Ag Cd Cs Ba Pb
Na 1
Mg 0.776 1
K 0.388 0.682 1
Ca 0.325 0.606 0.994 1
V 0.346 0.616 0.979 0.976 1
Cr 0.236 0.474 0.863 0.872 0.921 1
Mn -0.361 -0.175 0.559 0.629 0.570 0.464 1
Fe 0.373 0.642 0.996 0.997 0.979 0.890 0.575 1
Co -0.377 0.235 0.160 0.132 0.102 0.061 0.080 0.115 1
Ni -0.699 -0.732 -0.099 -0.012 -0.063 -0.103 0.767 -0.076 -0.065 1
Cu -0.722 -0.771 -0.172 -0.086 -0.136 -0.171 0.718 -0.150 -0.064 0.997 1
Zn -0.719 -0.760 -0.133 -0.045 -0.090 -0.117 0.745 -0.109 -0.073 0.999 0.998 1
Ga 0.436 0.665 0.986 0.986 0.956 0.860 0.550 0.991 0.091 -0.094 -0.166 -0.128 1
As 0.362 0.473 0.509 0.467 0.641 0.673 0.064 0.494 0.019 -0.275 -0.309 -0.275 0.467 1
Rb 0.317 0.706 0.955 0.933 0.967 0.876 0.456 0.939 0.292 -0.192 -0.259 -0.218 0.902 0.652 1
Sr 0.374 0.665 0.998 0.996 0.967 0.843 0.584 0.995 0.161 -0.070 -0.143 -0.105 0.988 0.456 0.937 1
Ag 0.288 0.106 -0.168 -0.184 -0.172 -0.032 -0.347 -0.149 0.009 -0.245 -0.228 -0.243 -0.044 0.165 -0.228 -0.171 1
Cd -0.696 -0.746 -0.136 -0.050 -0.098 -0.136 0.740 -0.114 -0.067 0.998 0.998 0.998 -0.126 -0.275 -0.229 -0.107 -0.198 1
Cs -0.186 -0.287 0.063 0.130 0.067 0.332 0.148 0.147 -0.314 0.102 0.083 0.113 0.128 -0.285 -0.037 0.079 -0.085 0.073 1
Ba 0.449 0.669 0.977 0.977 0.945 0.861 0.529 0.985 0.089 -0.111 -0.183 -0.146 0.999 0.462 0.888 0.980 0.002 -0.142 0.146 1
Pb -0.669 -0.325 0.248 0.305 0.261 0.195 0.822 0.244 0.496 0.787 0.767 0.774 0.221 0.000 0.235 0.269 -0.140 0.780 -0.100 0.206 1


3.3.2 组分特征PM2.5中的OC与EC的质量分布及OC/EC的比值是污染物来源解析的重要参数, 可以用来评判大气中各个污染物的相对贡献.不同国家和地区对柴油机的PM2.5组分研究结果基本一致, 即碳质组分为排放的PM2.5中最主要组分(冯小琼, 2013崔敏, 2016Vega et al., 2001).本研究中, 总含碳量(TC)占PM2.5总质量的77%~92%, 同其他非道路机械(Cui et al., 2016Liu et al., 2005)的碳质组分含量(70%~91%)接近.作业工况下OC、EC在PM2.5中的占比范围变化较大, OC/EC的比值的平均值为4.43, 用自展模拟的方法分析得到OC/EC的不确定范围是(-47%, 51%). OC/EC的比值受含硫率、空燃比、负荷等因素的影响较大(Liu et al., 2005), 在作业工况下受驾驶员的习惯和叉车负载变化的影响, 导致OC/EC的比值波动范围较大.崔敏(2016)利用车载试验得到挖掘机在作业工况下OC/EC的比值为0.6, 比本研究结果低, 但曲量等(2015)研究发现挖掘机OC/EC的比值的范围在2.5~3.0之间, 与本研究相近.由于OC/EC的比值随测试机械负荷的增加而降低(Liu et al., 2005), 本研究中测试的叉车在作业工况时是空载作业, 因此本研究的OC/EC的比值相对较高.
本研究检测出各类水溶性离子共计10种, 各种水溶性离子之和占PM2.5总质量的0.9%~6.7%.在水溶性离子中, NO3-和PO43-为主要组分, 分别占水溶性离子的33%和21%.其中, NO3-主要来源于燃油燃烧, 是尾气燃烧的二次转化产物, PO43-主要来自润滑油中的含磷添加剂(敖广, 2015).在CPCD40和CPCD50型号的叉车尾气中检测到了较高的SO42-含量, 由表 4的油品分析可知, 这两台叉车的柴油中含硫率显著高于其他叉车, 而尾气中的SO42-主要来源于燃料中的硫, 分析这两台叉车尾气的SO42-含量约为其他叉车平均值的13倍.同挖掘机的研究相比, 工作状态下挖掘机的各种水溶性离子之和占PM2.5总质量的0.5%~1.2% (崔敏, 2016), 小于叉车中水溶性离子的比例.计算每台叉车的SO42-/ NO3-的比值, 除CPCD50的叉车外, 其余均在0.06~0.66之间, 符合源解析中移动源SO42-/ NO3-<1的结果(Liu et al., 2014Yao et al., 2002).
PM2.5中共检测出21种金属组分, 其质量之和占总质量的0.33%~3.78%, 比例最高的组分分别是Ca、Na、Mg和Fe, 这4种组分的总质量占总体金属组分的85%.其中, Na为地壳元素(宋少洁等, 2011); Fe一方面来源于柴油中添加的助燃剂和消烟剂, 一方面来源于铁片的磨损(Srimuruganandam et al., 2012); Ca和Mg主要来源于钙基润滑油和清洁剂(崔敏, 2016).对于污染元素来说, Zn、Cu、Cr和Ni的比例较高, 对金属元素做相关性分析发现(表 9), Pb和Ni (r=0. 787), Pb和Cu (r=0.767), Pb和Zn (r=0.774), 而Pb被认为是来源于柴油车的排放(Manalis et al., 2005), 因此Zn、Cu和Ni这几种人为元素可能都来源于燃料油.在挖掘机的PM2.5组分中, 也检测出了较高比例的Cu、Zn、Ni和Cr, 与叉车尾气的PM2.5测试结果相近(崔敏, 2016).空气中的Cu、Cr等元素可以导致呼吸道疾病, 危害人体健康, 此类金属组分的排放应该受到重视.本研究得到V/Ni的平均值为0.59, 同PM2.5的源解析中V/Ni小于2相一致(Qin et al., 1997;Watson et al., 2001).
4 结论(Conclutions)1) 研究建立了实时燃油消耗速率-尾气气态污染物瞬时排放速率-单位燃油排放速率的关系.结果发现所有尾气的气态污染物变化均在怠速时波动最小, 作业时波动最大.燃油消耗速率与CO2、CO、NO和NOx的瞬时排放速率的波动变化完全一致;CO的单位燃油排放速率波动变化同CO2相反, NO和NO2的单位燃油排放速率的变化趋势基本相同, 由于部分NO2与NO的相互转化, NO2和NO的峰值发生时间不完全相同.
2) 基于现场测试的方法, 研究了功率、国标、工况对内燃叉车尾气污染物的排放影响.随着净功率的增加, 小时油耗量及CO、NO、NO2、NOx的小时排放量增大.国二叉车CO、NOx、PM2.5、SO2的排放因子较国一叉车分别降低了7.7%、16.0%、3.2%和53.6%.在3种工况下, CO、NO、NOx的排放因子都呈现出怠速>行走>作业的特征;而PM2.5的排放因子在行走和工作时较高.
3) 利用自展模拟的方法定量分析每种污染物的不确定性范围.结果表明, CO、NO、NO2、NOx、PM2.5、SO2的综合排放因子与不确定性范围分别是24.1 g·kg-1 (-25.4%, 26.1%)、26.5 g·kg-1 (-24.2%, 24.4%)、10.2 g·kg-1 (-22.3%, 21.9%)、36.7 g·kg-1 (-21.5%, 22.6%)、2.8 g·kg-1 (-25.5%, 24.7%)和0.9 g·kg-1 (-55.7%, 80.6%).
4) 内燃叉车PM2.5组分中, OC、EC、水溶性离子、金属组分的排放因子范围为别为2319~3983 mg·kg-1、465~3639 mg·kg-1、39~390 mg·kg-1和17~256 mg·kg-1. OC/EC的比值平均值为4.4, 不确定性范围是(-47%, 51%).水溶性离子中NO3-和PO43-最多, SO42-含量与燃油的含硫率相关.金属组分中Ca、Na、Mg和Fe占比最多, 其中Na为地壳元素, Ca、Mg和Fe则来源于燃料油, V/Ni的平均值为0.59. Cu、Zn、Ni和Cr等人为元素的排放值得关注.

参考文献
Abolhasani S, Frey H C, 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(8): 1033–1046.
敖广. 2015. 润滑油常用添加剂种类及作用[J]. 设备管理与维修, 2015(11): 30–33.
鲍晓峰, 丁焰, 尹航等. 2014. 非道路移动污染源排放清单编制技术指南(试行)[OL]. 北京: 环境保护部. 2014-12-31. http://www.qhepb.gov.cn/xwzx/zqyj/201407/W020140711608259681099
Cui M, Chen Y, Feng Y, et al. 2016. Measurement of PM and its chemical composition in real-world emissions from non-road and on-road diesel vehicles[J]. Atmospheric Chemistry & Physics, 17(11): 1–41.
崔敏. 2016. 基于实测的道路与非道路柴油车颗粒物排放特征初步研究[D]. 北京: 中国科学院大学http://www.irgrid.ac.cn/handle/1471x/1066412
Durbin T D, Cocker D R I, Sawant A A, et al. 2007. Regulated emissions from biodiesel fuels from on/off-road applications[J]. Atmospheric Environment, 41(27): 5647–5658.DOI:10.1016/j.atmosenv.2007.02.048
冯小琼. 2013. 珠江三角洲机动车VOCs和PM2. 5排放特征谱研究[D]. 广州: 华南理工大学http://cdmd.cnki.com.cn/Article/CDMD-10561-1013317566.htm
Frey H C, Rasdorf W, Kim K, et al. 2008. Comparison of Real-World Emissions of B20 Biodiesel Versus Petroleum Diesel for Selected Nonroad Vehicles and Engine Tiers[J]. Transportation Research Record Journal of the Transportation Research Board, 2058(2058): 33–42.
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 Journal of the Transportation Research Board, 2158(2010): 69–76.
Frey H C, Zheng J. 2002. Quantification of variability and uncertainty in air pollutant emission inventories:method and case study for utility NOx emissions[J]. Air Repair, 52(9): 1083–1095.
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(20): 255–261.
贾旭, 封学军, 蒋柳鹏. 2014. 港口作业机械大气污染物排放研究[J]. 华东交通大学学报, 2014(3): 12–17.
Lewis P. 2009. Estimating fuel use and emission rates of nonroad diesel construction equipment performing representative duty cycles[D]. Dissertations & Theses-Gradworkshttps://search.proquest.com/docview/304964388
Lewis P, Asce S M, Rasdorf W, et al. 2009. Requirements and incentives for reducing construction vehicle emissions and comparison of nonroad diesel engine emissions data sources[J]. Journal of Construction Engineering & Management, 135(5): 341–351.
Lewis P, Leming M, Rasdorf W. 2011. Impact of engine idling on fuel use and CO2 emissions of nonroad diesel construction equipment[J]. Journal of Management in Engineering, 28(1): 31–38.
李东玲, 吴烨, 周昱, 等. 2012. 我国典型工程机械燃油消耗量及排放清单研究[J]. 环境科学, 2012, 33(2): 518–524.
Lindgren M, Arrhenius K, Larsson G, et al. 2011. Analysis of unregulated emissions from an off-road diesel engine during realistic work operations[J]. Atmospheric Environment, 45(30): 5394–5398.DOI:10.1016/j.atmosenv.2011.06.046
Lindgren M, Hansson P A. 2004. Effects of transient conditions on exhaust emissions from two non-road diesel engines[J]. Biosystems Engineering, 87(1): 57–66.DOI:10.1016/j.biosystemseng.2003.10.001
Lindgren M, Larsson G, Hansson P A. 2010. Evaluation of factors influencing emissions from tractors and construction equipment during realistic work operations using diesel fuel and bio-fuels as substitute[J]. Biosystems Engineering, 107(2): 123–130.DOI:10.1016/j.biosystemseng.2010.07.010
Liu J, Li J, Zhang Y, et al. 2014. Source apportionment using radiocarbon and organic tracers for PM2.5 carbonaceous aerosols in Guangzhou, South China:Contrasting local-and regional-scale haze events[J]. Environmental Science & Technology, 48(20): 12002–12011.
Liu Z F, Lu M, Birch M E, et al. 2005. Variations of the particulate carbon distribution from a nonroad diesel generator[J]. Environmental Science & Technology, 39(20): 7840–7844.
刘旭刚, 张国彬, 蔚建玺, 等. 2007. 柴油发动机高压共轨技术的产生与发展[J]. 汽车运用, 2007(11): 30–31.
鲁君, 黄成, 胡磬遥, 等. 2017. 长三角地区典型城市非道路移动机械大气污染物排放清单[J]. 环境科学, 2017, 38(7): 1–14.
Manalis N, Grivas G, Protonotarios V, et al. 2005. Toxic metal content of particulate matter (PM10), within the Greater Area of Athens[J]. Chemosphere, 60(4): 557–566.DOI:10.1016/j.chemosphere.2005.01.003
潘跃, 吴子文. 2008. 柴油机涡轮增压中冷技术及对排放控制分析[J]. 农业装备与车辆工程, 2008(1): 48–49.
Qin Y, Chan C K, Chan L Y. 1997. Characteristics of chemical compositions of atmospheric aerosols in Hong Kong:spatial and seasonal distributions[J]. Science of the Total Environment, 206(1): 25–37.DOI:10.1016/S0048-9697(97)00214-3
曲亮, 何立强, 胡京南, 等. 2015. 工程机械不同工况下PM2.5排放及其碳质组分特征[J]. 环境科学研究, 2015, 28(7): 1047–1052.
宋少洁, 吴烨, 蒋靖坤, 等. 2012. 北京市典型道路交通环境细颗粒物元素组成及分布特征[J]. 环境科学学报, 2012, 32(1): 66–73.
Srimuruganandam B, Nagendra S M. 2012. Source characterization of PM10, and PM2.5, mass using a chemical mass balance model at urban roadside[J]. Science of the Total Environment, 433(7): 8–19.
谭华, 刘娟, 沈寅, 等. 2013. 码头港作机械大气污染物排放清单研究[J]. 环境科学与管理, 2013, 38(6): 82–87.
Vega E, Mugica V, Reyes E, et al. 2001. Chemical composition of fugitive dust emitters in Mexico City[J]. Atmospheric Environment, 35(23): 4033–4039.DOI:10.1016/S1352-2310(01)00164-9
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
魏名山, 马朝臣, 张富红, 等. 2002. 云内4100柴油机增压中冷机型的试验研究[J]. 内燃机工程, 2002, 23(4): 46–47.
Yan F, Winijkul E, Streets D G, et al. 2014. Global emission projections for the transportation sector using dynamic technology modeling[J]. Atmospheric Chemistry & Physics, 13(9): 23373–23419.
Yao X, Chan C K, Fang M, et al. 2002. The water-soluble ionic composition of PM2.5 in Shanghai and Beijing, China[J]. Atmospheric Environment, 36(26): 4223–4234.DOI:10.1016/S1352-2310(02)00342-4
殷小鸽. 2010. 非道路移动源燃油品质改善对污染排放影响的成本效益分析[D]. 天津: 南开大学http://www.doc88.com/p-0049218138826.html
张礼俊, 郑君瑜, 尹沙沙, 等. 2010. 珠江三角洲非道路移动源排放清单开发[J]. 环境科学, 2010, 31(4): 886–891.
Zheng J, Frey H C. 2004. Quantification of variability and uncertainty using mixture distributions:evaluation of sample size, mixing weights, and separation between components[J]. Risk Analysis, 24(3): 553–571.DOI:10.1111/risk.2004.24.issue-3
郑君瑜. 2014. 区域高分辨率大气排放源清单建立的技术方法与应用[M]. 北京: 科学出版社.
中国机械工业联合会. 2010-2015. 中国工程机械工业年鉴[M]. 北京: 机械工业出版社
朱倩茹, 廖程浩, 刘剑筠, 等. 2015. 广东省典型内河港区主要大气污染物排放特征研究[J]. 安全与环境学报, 2015, 15(3): 191–195.




相关话题/测试 质量 机械 金属 综合