环境科学  2022, Vol. 43 Issue (4): 1777-1787   PDF    
基于隧道测试的机动车VOCs排放特征及源解析
刘鑫会1, 朱仁成1, 金博强1, 梅慧1, 祖雷2, 尹沙沙1,3, 张瑞芹1,3, 胡京南2     
1. 郑州大学生态与环境学院, 郑州 450001;
2. 中国环境科学研究院, 北京 100012;
3. 郑州大学环境科学研究院, 郑州 450001
摘要: 为探究以乙醇汽油(E10)为主要燃料的机动车尾气源和蒸发源挥发性有机物(VOCs)排放特征, 于2019年12月在郑州市北三环隧道内展开了连续两周的VOCs采样, 并对隧道内车流特征和环境参数等进行在线监测.首先, 利用气相色谱/质谱(GC/MS)法定量出106种VOCs组分, 然后采用正交矩阵因子分析(PMF5.0)-化学质量平衡(CMB8.2)复合模型对机动车尾气源和蒸发源VOCs排放的贡献率进行定量解析, 最后采用最大增量反应活性(MIR)和气溶胶生成系数(FAC)分别测算了尾气源和蒸发源的臭氧生成潜势(OFP)和二次有机气溶胶生成潜势(SOAFP).结果表明, 采样期间隧道环境空气中ρ(VOCs)为(2794.5±147.4) μg ·m-3, 其中卤代烃类的质量分数最高[(32.4±2.0)%], 其次为芳烃类[(27.5±0.6)%]和烷烃类[(23.3±0.8)%]; 基于PMF5.0-CMB8.2复合模型的机动车源VOCs解析结果为:尾气排放(62.5%)>蒸发排放(37.5%); 机动车源VOCs排放的OFP贡献率为:尾气排放(71.9%)>蒸发排放(28.1%), SOAFP贡献率为:尾气排放(75.8%)>蒸发排放(24.2%); 蒸发源OFP的优势组分有间/对-二乙苯、异戊二烯和反-2-戊烯等, 蒸发源SOAFP的优势组分有间/对-二乙苯、间/对-二甲苯和1, 2, 3-三甲基苯等, 尾气源OFP的优势组分有间/对-二甲苯、1, 2, 4-三甲基苯和1, 3, 5-三甲基苯等, 尾气源SOAFP的优势组分有间/对-二甲苯、间/对-二乙苯和1, 3, 5-三甲基苯等.建议E10使用区域, 在加强机动车尾气排放控制的同时, 也应重视蒸发VOCs排放, 尤其是其中芳烃类和烯烃类等高活性组分.
关键词: 隧道测试      挥发性有机物(VOCs)      源解析      蒸发排放      尾气排放      臭氧生成潜势(OFP)      二次有机气溶胶生成潜势(SOAFP)     
Characteristics and Source Apportionment of Vehicular VOCs Emissions in a Tunnel Study
LIU Xin-hui1 , ZHU Ren-cheng1 , JIN Bo-qiang1 , MEI Hui1 , ZU Lei2 , YIN Sha-sha1,3 , ZHANG Rui-qin1,3 , HU Jing-nan2     
1. School of Ecology and Environment, Zhengzhou University, Zhengzhou 450001, China;
2. Chinese Research Academy of Environmental Sciences, Beijing 100012, China;
3. Institute of Environmental Science, Zhengzhou University, Zhengzhou 450001, China
Abstract: To explore the emission characteristics of volatile organic compounds (VOCs) from vehicular exhaust sources and evaporative sources with ethanol gasoline (E10) as the main fuel, VOCs sampling campaigns were carried out in the north third ring tunnel of Zhengzhou city for two consecutive weeks in December 2019. In addition, the characteristics of traffic flow and environmental information were also monitored in the tunnel. Firstly, 106 VOCs were quantified using gas chromatography/mass spectrometry (GC/MS), and then source apportionment of VOCs in the tunnel was carried out using a positive matrix factorization (PMF5.0)-chemical mass balance (CMB8.2) composite model. Finally, the ozone formation potential (OFP) and secondary organic aerosol formation potential (SOAFP) of vehicle exhaust sources and evaporative sources were analyzed using the maximum incremental reactivity (MIR) and fractional aerosol coefficient (FAC). The results showed that ρ(VOCs) in the tunnel was (2794.5±147.4) μg ·m-3 during the experiment, among which halogenated hydrocarbons [(32.4±2.0)%] accounted for the highest proportion, followed by aromatic hydrocarbons [(27.5±0.6)%] and alkanes [(23.3±0.8)%]. Source apportionment of vehicular VOCs showed that exhaust emissions (62.5%)>evaporative emissions (37.5%), whereas the contribution of OFP was that exhaust emissions (71.9%)>evaporative emissions (28.1%), and the contribution of SOAFP was that exhaust emissions (75.8%)>evaporative emissions (24.2%). The dominant components of OFP in evaporative sources were m, p-diethylbenzene, isoprene, and trans-2-pentene, whereas m, p-diethylbenzene, m, p-xylene, and 1, 2, 3-trimethylbenzene were the dominant components of SOAFP. The major components of OFP in exhaust sources were m, p-xylene, 1, 2, 4-trimethylbenzene, and 1, 3, 5-trimethylbenzene, whereas m, p-xylene, m, p-diethylbenzene, and 1, 3, 5-trimethylbenzene were the dominant components of SOAFP. In regions where ethanol gasoline is used, special attention should be paid not only to the exhaust emissions control but also to strengthening the emissions reduction of VOCs from vehicle evaporative sources, especially the high active components such as aromatic hydrocarbons and alkenes.
Key words: tunnel test      volatile organic compounds (VOCs)      source apportionment      evaporative emissions      exhaust emissions      ozone formation potential (OFP)      secondary organic aerosol formation potential (SOAFP)     

近年来, 臭氧(ozone, O3)与细颗粒物(fine particulate matter, PM2.5)协同控制成为提升我国环境空气质量的关键所在[1].作为大气O3和二次有机气溶胶(secondary organic aerosol, SOA)共同的关键前体物[2, 3], 挥发性有机物(volatile organic compounds, VOCs)的污染源排放特征、来源解析和二次污染物生成潜势等逐渐成为国内外研究热点.北京[4]、天津[5]、河北[6]、上海[7]、厦门[8]和西安[9]等地区2014~2019年间开展的大气VOCs源解析结果表明, 机动车源的贡献率高达23.1% ~48.4%, 是城市大气VOCs的主要来源.因此, 加强机动车VOCs排放控制是目前我国大气污染防治的重要组成部分.

作为我国重要的综合交通枢纽和商贸物流中心, 河南省省会郑州市目前机动车保有量已超过400万辆, 在全国城市排名中位居第六[10].为降低机动车污染排放, 郑州市作为首批试点城市于2004年12月1日已实现了乙醇汽油(E10)的全面封闭式推广使用[11].在减污降碳协同背景下, 国家能源局已发布实施方案要求积极做好车用乙醇汽油推广的重点工作[12].但由于E10汽油的含氧量和饱和蒸气压与普通汽油存在差异, 其降污减排效果在实际应用过程中仍存在一定的不确定性[13~16].例如, Wu等[13]的研究结果显示与普通汽油相比, E10能降低尾气颗粒物排放, 但对THC并无明显减排效果; Zhu等[14]的研究也发现常温下E10能降低大多数污染物排放, 但是低温下E10的节能减排效果会显著下降; 而Koupal等[15]基于MOVES-Mexico模型研究则认为推广E10后城市区域VOCs将增加13% ~18%, 此外还会增加有毒污染物排放.

目前, 用于机动车VOCs排放特征研究的测试方法主要有台架试验、车载实验、路边采样、遥感技术和隧道测试等.与其他方法相比, 隧道测试可以反映车队实际运行状态下蒸发源与尾气源的综合排放水平, 且具有光化学反应弱、受外源传输影响小等优点.我国上海[17]、天津[18]、广州[19]、合肥[20]、南京[21]和成都[22]等城市典型隧道已开展了相关测试研究, 获得了本地化的机动车VOCs排放特征, 但缺少对尾气源和蒸发源贡献率的进一步探究.此外, 上述研究区域的车用燃油多为国五及以下标准的传统燃油, 不能反映机动车使用我国目前大力推广的国六阶段乙醇汽油的VOCs排放特征.

为获得以E10为主要燃料的机动车VOCs排放特征, 本研究于2019年12月在郑州市北三环隧道开展了连续两周的VOCs采样, 利用气相色谱-火焰离子化检测器/质谱(gas chromatography-flame ionization detector/mass spectrometry, GC-FID/MS)方法对隧道内106种VOCs组分进行定量分析, 结合正交矩阵因子分析(positive matrix factorization, PMF5.0)-化学质量平衡(chemical mass balance, CMB8.2)复合模型探究机动车尾气源与蒸发源VOCs排放贡献率, 并采用最大增量反应活性(maximum incremental reactivity, MIR)和气溶胶生成系数(fractional aerosol coefficient, FAC)测算尾气源与蒸发源的臭氧生成潜势(ozone formation potential, OFP)和二次有机气溶胶生成潜势(secondary organic aerosols formation potential, SOAFP).本研究充分利用了隧道环境受外源传输影响小、光化学反应弱和机动车排放贡献大等优势, 将隧道测试与受体模型源解析相结合, 提高解析结果的准确性和可靠性, 力图探明E10乙醇汽油车在实际行驶状态下尾气源和蒸发源VOCs排放的贡献率, 以期为制定科学有效的机动车VOCs管控政策提供依据.

1 材料与方法 1.1 样本采集

本研究基于郑州市北三环隧道(113°43′27″ E, 34°48′50″ N)开展连续两周的隧道环境空气VOCs采样.该隧道采用独立双孔箱涵形式, 隧道单孔宽13.5 m, 单孔横截面积83.3 m2, 全长4 520 m, 是目前河南省最长的城市隧道.采样点布设在隧道南侧孔道入口870 m处, 采样高度为1.5 m, 北三环隧道地理位置及采样点布设如图 1所示.VOCs样本采集时间为2019年12月5~18日早中晚3个时段(08:00~09:00、11:00~12:00和17:00~18:00), 包含工作日和非工作日.采样期间隧道顶部设置的全射流式轴流风机处于关闭状态, 除隧道进出口外无其他通风口, 隧道附近无特别工业源, 因此采样点VOCs主要来自机动车排放累积, 受外源传输影响较小.

图 1 北三环隧道地理位置及采样点布设示意 Fig. 1 Location of the north third ring tunnel and sampling site

1.2 仪器设备

采用数字高清双雷达微波车辆检测器(Wavetronix, USA)结合摄像机录像获取采样期间车流量和车型占比等信息, 采用便携式手持气象仪(YGY-QXY, China)对采样点处风速进行在线监测.VOCs样本采用3.2 L内壁硅烷化、抽真空的Summa罐(Entech, USA)以0.04 L·min-1恒定流量采集, 采样时间为60 min. VOCs样本浓缩除杂后进入GC-FID/MS (Agilent, USA)进行检测分析.根据色谱保留时间和质谱图进行定性分析, 以EPA TO-15和PAMS作为外标气体, 以溴氯甲烷、1, 2-二氯苯、氯苯-d5和4-溴氟苯作为内标气体进行定量分析, 共检测出106种VOCs组分.为降低高反应活性组分的检测误差, 所有样本于15 d内完成分析[23].

1.3 分析方法 1.3.1 PMF5.0模型

PMF5.0模型的基本原理是将样本组分浓度矩阵分解为源贡献矩阵、源成分谱矩阵和残差矩阵, 并使用多次迭代算法使目标函数Q最小化, 进而确定最优解析结果, 计算公式如式(1)所示:

(1)

式中, xij为样本i中组分j的浓度; p为污染源个数; gik为污染源k对样本i的贡献; fkj为污染源k中组分j的含量; eij为残差; Q为累计残差; uijxij的不确定度.

PMF5.0模型需要受体点检测浓度和不确定度数据, 当检测浓度低于方法检出限(method detection limit, MDL)时, 浓度取(1/2)MDL, 不确定度取(5/6)MDL; 当检测浓度高于MDL时, 不确定度计算公式如式(2)所示[24]

(2)

式中, EF为误差分数, 以相对标准偏差(relative standard deviation, RSD)表示; c为检测浓度, μg·m-3; MDL为方法检出限, μg·m-3.

1.3.2 CMB8.2模型

CMB8.2模型的基本原理是假设各污染源排放组分与受体点检测组分遵循质量守恒原则, 将受体点处各组分检测浓度分配至不同污染源中, 并依据有效方差最小二乘法求解出模型拟合的最优解[25], 计算公式如式(3)所示:

(3)

式中, cik为采样点处第k个样本中组分i的检测浓度, μg·m-3; aij为第j个污染源成分谱中组分i的质量分数; Sjk为第j个污染源对样本k的贡献.

CMB8.2模型需要污染源成分谱和不确定度数据、受体点VOCs检测浓度和不确定度数据.成分谱数据的不确定度设为自身质量分数的15%, 受体点浓度数据的不确定度计算公式如式(4)所示:

(4)

式中, σcik为第k个样本中组分i的不确定度, μg·m-3; MDLi为组分i的方法检出限, μg·m-3; CVi为组分i的相对标准偏差; cik为第k个样本中组分i的检测浓度, μg·m-3.

1.3.3 OFP分析

采用MIR系数对机动车尾气源和蒸发源VOCs排放的OFP进行测算[26], 计算公式如式(5)所示:

(5)

式中, OFPi为组分i的臭氧生成潜势, μg·m-3; ci为组分i的质量浓度, μg·m-3; MIRi为组分i的最大增量反应活性系数, μg·μg-1.

1.3.4 SOAFP分析

采用FAC方法对机动车尾气源和蒸发源VOCs排放的SOAFP进行测算[27, 28], 计算公式如式(6)所示:

(6)

式中, SOAFPi为组分i的二次有机气溶胶生成潜势, μg·m-3; ci为组分i的检测浓度, μg·m-3; FACi为组分i的气溶胶生成系数, μg·μg-1; Fri为组分i参与反应的质量分数, %.

2 结果与讨论 2.1 测试参数日变化特征

对实验期间(包含工作日和非工作日各两天)的车流参数及风速日变化特征进行统计分析, 结果如图 2所示.隧道内日均车流量达(43 228±2 643) veh·d-1, 并在上下班高峰期呈明显双峰分布特征.隧道内日均风速为(3.6±0.1) m·s-1, 日变化特征与车流量特征相近且白天风速显著高于夜晚, 这是由于白天车辆行驶促进了隧道内空气对流, 而夜晚仅受气象条件影响, 隧道内空气流通性较弱.实验期间日均车型占比为:小型车[(85.6±1.2)%]>中型车[(12.9±1.0)%]>大型车[(1.5±0.3)%].凌晨03:00附近大中型车占比较高, 主要是因为白天北三环隧道对大型货车限行, 而夜晚有少量大型货车与清扫车通过.

图 2 隧道内车流参数及风速日变化特征 Fig. 2 Daily change characteristics of traffic flow parameters and wind speed in the tunnel

本实验期间3个VOCs采样时段的车流参数、风速和ρ(VOCs)平均变化特征汇总如表 1所示.下午风速略低于上午和中午, 主要是由于下午隧道内外温差略低于上午和中午.中午小型车占比略高于上午和下午, 但车流量略低于上午和下午.由于白天隧道内对大型货车限行且3个采样时段均为车流高峰期, 因此大型车占比均处于较低水平. ρ(VOCs)呈早中晚逐渐升高趋势, 主要受隧道内机动车VOCs排放累积所致.

表 1 不同采样时段的车流参数、风速及VOCs浓度变化特征 Table 1 Characteristics of traffic flow parameters, wind speed and VOCs concentration at different sampling periods

2.2 隧道环境空气VOCs成分谱特征

对测试期间采集的42个VOCs样品平均, 所得隧道环境空气VOCs成分谱如表 2所示.隧道内ρ(VOCs)的平均值为(2 794.5±147.4) μg·m-3, 远高于天津五经路隧道(257.4 μg·m-3)和南京富贵山隧道(225.6 μg·m-3)的测试结果, 与2005年广州珠江隧道测试结果相近(2 216.2 μg·m-3)[18], 这主要与隧道内车流量大且分析的VOCs组分多等因素有关.在所检测的106种VOCs组分中, 卤代烃类的质量分数最高[(32.4±2.0)%], 其后依次为芳烃类[(27.5±0.6)%]、烷烃类[(23.3±0.8)%]、含氧挥发性有机物(oxygenated volatile organic compounds, OVOCs) [(9.8±0.8)%]、烯烃类[(6.2±0.3)%]、二硫化碳[(0.6±0.1)%]和乙炔[(0.2±0.1)%].

表 2 隧道环境空气VOCs成分谱特征 Table 2 Characteristics of VOCs composition profile in the north third ring tunnel

卤代烃类的优势组分主要有1, 2-二氯丙烷、1, 2, 4-三氯苯和1, 1, 2, 3, 4, 4-六氯-1, 3-丁二烯等, 这些组分化学活性较低, 可在外界环境中实现远距离输送[29, 30], 主要来自外界氯化工业的气团传输和机动车排放[31, 32], 其中氯甲烷等卤代烷烃可作为生物质燃烧示踪组分[33], 卤代烯烃主要来自石油化工排放[34].芳烃类的优势组分有间/对-二甲苯、间/对-二乙苯、1, 2, 4-三甲基苯和苯乙烯等, 其中BTEX(苯、甲苯、乙苯和二甲苯的统称)主要来自溶剂使用和机动车排放[6, 33].烷烃类的优势组分有十二烷、十一烷、环己烷、正癸烷和正己烷等, 其中长链烷烃是重型柴油车排放的特征组分, 而C4~C7烷烃类主要来自机动车蒸发排放[9, 35]. OVOCs的优势组分有2-己酮、乙酸乙烯酯、甲基丙烯酸甲酯和4-甲基-2-戊酮等, 其中乙醇和甲基叔丁基醚(methyl tert-butyl ether, MTBE)通常作为提高汽油辛烷值的添加剂, 是蒸发排放的示踪性组分[30, 36], 醛酮化物与E10汽油不完全燃烧排放有关[13].烯烃类化学反应活性较高, 在外界环境中容易被消耗, 因此受气团传输影响小, 主要来自机动车排放累积, 优势组分有1-己烯、异戊二烯和1-戊烯等, 其中异戊二烯受外界植物排放和机动车排放的共同影响[36, 37], 1, 3-丁二烯可作为机动车尾气排放示踪组分[34], C3~C4烯烃类主要来自生物质燃烧[35].乙炔通常作为尾气排放示踪组分[30, 37].综上所述, 隧道环境空气VOCs主要受外源传输(背景源)、机动车尾气排放和蒸发排放这3个方面的贡献.

2.3 隧道环境空气VOCs源解析 2.3.1 PMF5.0模型解析

为开展隧道内VOCs源解析, 本研究筛选了37种VOCs组分并采用PMF5.0模型多次调试, 解析结果如图 3所示.因子数为6时, Q值发生突减, 此后继续增加因子数Q值变化不大, Qtrue=Qrobust, 且解析结果稳定性较好, 因此适宜因子个数为6.运行结果显示, 所有组分模型的测算数据与实测数据的拟合优度(R2)均超过0.85.

①氯甲烷, ②二氯甲烷, ③ 1, 1-二氯乙烷, ④三氯甲烷, ⑤异丁烷, ⑥正丁烷, ⑦异戊烷, ⑧丙烷, ⑨ 2, 4-二甲基戊烷, ⑩ 2, 3-二甲基戊烷, ⑪ 甲基环己烷, ⑫ 1-丁烯, ⑬ 1, 3-丁二烯, ⑭ 反-2-丁烯, ⑮ 顺-2-丁烯, ⑯ 1-戊烯, ⑰ 异戊二烯, ⑱ 反-2-戊烯, ⑲ 顺-2-戊烯, ⑳ 苯, ㉑ 甲苯, ㉒ 乙苯, ㉓ 间/对-二甲苯, ㉔ 苯乙烯, ㉕ 邻-二甲苯, ㉖ 异丙苯, ㉗ 正丙苯, ㉘ 间-乙基甲苯, ㉙ 对-乙基甲苯, ㉚ 1, 3, 5-三甲基苯, ㉛ 邻-乙基甲苯, ㉜ 1, 2, 4-三甲基苯, ㉝ 1, 2, 3-三甲基苯, ㉞ 间/对-二乙苯, ㉟ 2-丙烯醛, ㊱ 丙酮, ㊲ 甲基叔丁基醚 图 3 PMF5.0模型解析结果 Fig. 3 Source apportionment results of PMF5.0 model

因子1中优势组分有三氯甲烷、苯、间/对-二甲苯、邻-二甲苯和间/对-二乙苯等, 以上组分主要来自橡胶工业和有机合成等工业源[38, 39], 将其识别为“工业排放”.因子2中氯甲烷、二氯甲烷、正丁烷、异丁烷和丙烷的贡献率较高, 其中氯甲烷和二氯甲烷是生物质燃烧示踪组分[33], 丙烷、正丁烷和异丁烷主要来自居民日常生活中液化石油气(LPG)/天然气(NG)等燃料的使用[40], 将其识别为“燃烧+燃气使用”.因子3中丙酮和二氯甲烷的贡献率较高, 这两种组分的化学反应活性较低, 主要来自外界老化气团传输[38], 将其识别为“老化气团”.因子4中烯烃类和芳烃类贡献率较高, 其中BTEX和1, 3-丁二烯为机动车尾气排放特征组分[5, 33], 甲基叔丁基醚为蒸发排放示踪组分[36], 将其识别为“机动车排放1”.因子5的优势组分主要是间/对-二乙苯、苯乙烯、1, 3, 5-三甲基苯和1, 2, 4-三甲基苯等芳烃类, 主要来自油漆和涂料等溶剂使用[34, 39], 将其识别为“溶剂使用”.因子6中低碳烷烃、异戊烷、苯、甲苯和2-丙烯醛等贡献率较高, 为机动车排放的特征组分[37, 40, 41], 将其识别为“机动车排放2”.各因子累计贡献率分别为:机动车排放(62.4%)>燃烧+燃气使用(15.8%)>溶剂使用(9.5%)>工业排放(8.4%)>老化气团(3.9%).

3个采样时段的PMF5.0模型解析结果如图 4所示.“机动车排放”的贡献率特征为:下午>上午>中午, 与车流量变化特征一致.“老化气团”的贡献率与“机动车排放”特征相反, 这主要是受隧道内光化学反应弱且机动车VOCs排放逐渐累积等因素影响.“燃烧+燃气使用”的贡献率在中午最高, 这主要受周围餐饮业燃气使用等因素影响.“工业排放”和“溶剂使用”的贡献率相近.与郑州市环境空气VOCs来源解析结果相比[38, 39], 本研究中“机动车排放”的贡献率显著上升, 这主要与隧道内车流量大且空间相对密闭有关, 而其余源(如工业排放和溶剂使用等)的相对贡献率较为相符, 说明解析结果基本可信.

图 4 不同采样时段PMF5.0模型解析结果 Fig. 4 Source apportionment results of PMF5.0 model at different sampling periods

2.3.2 CMB8.2模型解析

由于机动车尾气源与蒸发源的成分谱相近, 本研究中PMF5.0模型无法将二者有效分离.Shi等[42]的研究采用PMF-CMB复合模型对郑州市PM2.5进行源解析, 发现该法可较好解决共线性问题.因此, 本研究将PMF5.0模型解析所得机动车源成分谱纳入CMB8.2模型, 对尾气源与蒸发源VOCs排放的贡献率进行深入解析.

结合成分谱数据库[41, 43], 采用CMB8.2模型进行多次调试并筛除共线性较强的源成分谱, 解析结果如图 5所示.蒸发排放与尾气排放的贡献值均为正, 浓度平衡值(%MASS)达84.5%, 模型拟合效果较优.蒸发源和尾气源VOCs排放的贡献率分别为37.5%和62.5%, 其中蒸发源贡献率略高于Liu等[44]基于排放清单的测算结果(31.8%), 明显高于高丹丹等[45]测算的河南省2016年机动车排放清单中汽油车蒸发源VOCs的排放占比(24.6%), 这主要与2018年7月以后河南省实现国六车用乙醇汽油的全面供应等因素有关.考虑到本研究采样时间为低温冬季, 年均蒸发排放水平较本研究结果更高, 因此未来应进一步加强对E10乙醇汽油车蒸发源VOCs排放的管控力度.

横坐标符号含义同图 3 图 5 CMB8.2模型解析结果 Fig. 5 Source apportionment results of CMB8.2 model

在早、中、晚这3个采样时段基于CMB8.2模型的机动车源VOCs解析结果如图 6所示. 08:00~09:00时段蒸发排放贡献率略高, 可能与早高峰期间车流量大, 交通拥堵, 导致车辆怠速行驶等因素有关. 11:00~12:00时段蒸发排放贡献率高于17:00~18:00时段, 主要是由于中午温度升高使机动车蒸发排放强度略有增加.综合来看, 蒸发源与尾气源VOCs排放的贡献率变化波动较小, 主要与3个采样时段的车型占比相对稳定有关.需要说明的是, 本研究CMB8.2模型采用了国外的尾气源谱数据, 由于国内外机动车排放标准与发动机工况存在差异, 会给本研究的解析结果带来一定的不确定性, 未来可进一步加强基于国内汽车燃用E10汽油的VOCs排放源谱的完善和更新, 为源解析研究提供支持.

图 6 不同采样时段CMB8.2模型解析结果 Fig. 6 Source apportionments of CMB8.2 model at different sampling periods

2.4 尾气源与蒸发源二次污染特征分析

为判断机动车尾气源和蒸发源VOCs排放对O3和SOA二次生成的优势组分和贡献率, 根据PMF5.0-CMB8.2复合模型的解析结果, 采用MIR和FAC方法对两源中各组分和各类别的OFP和SOAFP的贡献特征展开分析, 如图 7图 8所示.

横坐标符号含义同图 3 图 7 尾气排放与蒸发排放各组分的OFP和SOAFP Fig. 7 OFP and SOAFP of various VOCs species in exhaust emissions and evaporative emissions

贡献率低于1%的类别数据已隐藏标签 图 8 尾气排放与蒸发排放各类别的OFP和SOAFP Fig. 8 OFP and SOAFP of different groups in exhaust emissions and evaporative emissions

所解析37种组分的总OFP为3 707.2 μg·m-3, 其中蒸发源和尾气源的贡献率分别为28.1%和71.9%.蒸发源OFP贡献较高的前5种组分及其贡献率为:间/对-二乙苯(12.0%)、异戊二烯(11.0%)、反-2-戊烯(8.8%)、2-丙烯醛(8.4%)和间/对-二甲苯(7.0%), 各类别贡献率为:芳烃类(52.2%)>烯烃类(33.8%)>OVOCs(9.7%)>烷烃类(4.1%)>卤代烃类(0.2%).尾气源OFP贡献较高的前5种组分及其贡献率为:间/对-二甲苯(14.2%)、1, 2, 4-三甲基苯(12.4%)、1, 3, 5-三甲基苯(10.8%)、邻-二甲苯(6.5%)和1, 2, 3-三甲基苯(6.2%), 各类别贡献率为:芳烃类(77.7%)>烯烃类(20.2%)>烷烃类(1.7%)>OVOCs (0.4%).

所解析37种组分的SOAFP共计28.0 μg·m-3, 其中蒸发源与尾气源的贡献率分别为24.2%和75.8%.蒸发源SOAFP贡献较高的前5种组分及其贡献率为:间/对-二乙苯(37.9%)、间/对-二甲苯(9.8%)、1, 2, 3-三甲基苯(5.8%)、甲基环己烷(5.2%)和对-乙基甲苯(5.1%), 各类别贡献率为:芳烃类(91.1%)>烷烃类(5.2%)>烯烃类(3.7%).尾气源SOAFP贡献较高的前5种组分及其贡献率为:间/对-二甲苯(16.3%)、间/对-二乙苯(14.2%)、1, 3, 5-三甲基苯(12.9%)、1, 2, 4-三甲基苯(8.4%)和间-乙基甲苯(7.8%), 各类别贡献率为:芳烃类(99.7%)>烷烃类(0.2%)>烯烃类(0.1%).

与蒸发源和尾气源VOCs排放的贡献率对比可知, 蒸发源VOCs活性低于尾气源, 这是由于高反应活性的C9芳烃类和C4烯烃类主要来自尾气排放, 而低反应活性的烷烃类及卤代烃类在蒸发源中占比较高.综上可知, 着重加强尾气源中芳烃类的减排力度是治理郑州市冬季雾-霾污染与夏季臭氧污染的重中之重.此外, 考虑到夏季高温会大大增加蒸发源的排放强度, 因此蒸发源中芳烃类和烯烃类等高活性组分也应引起高度重视.

3 结论

(1) 测试期间隧道环境空气的ρ(VOCs)为(2 794.5±147.4) μg·m-3, 其中卤代烃类和芳烃类的贡献率较高, 成分谱的优势组分有1, 2-二氯甲烷、间/对-二乙苯、间/对-二甲苯、1, 2, 4-三甲基苯和十二烷等.

(2) PMF5.0-CMB8.2复合受体模型的解析结果为尾气源(62.5%)>蒸发源(37.5%), 由于采样时间为低温冬季, 蒸发排放强度低于年均水平, 因此未来在加严机动车尾气排放标准的同时, 蒸发排放管控同样不容忽视.

(3) 尾气源和蒸发源VOCs排放的OFP贡献率分别为71.9%和28.1%, SOAFP贡献率分别为75.8%和24.2%, 其中芳烃类在尾气源和蒸发源中对OFP和SOAFP的贡献率均居于首位, 是汽油车VOCs排放的优先控制组分.

(4) 由于采样点选取在车流量大、空间相对密闭且光化学反应弱的隧道内, 基于PMF5.0模型解析的机动车源贡献率显著高于郑州市环境空气中VOCs的解析结果, 而工业排放和溶剂使用等来源的相对贡献率与前人研究结论相符.

致谢: 感谢郑州市生态环境局为本研究顺利开展工作提供的大力支持和帮助; 感谢郑州市郑东新区龙湖区隧道管理中心为隧道内大气VOCs采样提供的实验条件和安全保障措施.

参考文献
[1] 中华人民共和国生态环境部. 2020中国生态环境状况公报[EB/OL]. http://www.mee.gov.cn/hjzl/sthjzk/zghjzkgb/, 2021-05-26.
[2] 陈文泰, 邵敏, 袁斌, 等. 大气中挥发性有机物(VOCs)对二次有机气溶胶(SOA)生成贡献的参数化估算[J]. 环境科学学报, 2013, 33(1): 163-172.
Chen W T, Shao M, Yuan B, et al. Parameterization of contribution to secondary organic aerosol (SOA) formation from ambient volatile organic compounds (VOCs)[J]. Acta Scientiae Circumstantiae, 2013, 33(1): 163-172.
[3] Rohrer F, Lu K D, Hofzumahaus A, et al. Maximum efficiency in the hydroxyl-radical-based self-cleansing of the troposphere[J]. Nature Geoscience, 2014, 7(8): 559-563. DOI:10.1038/ngeo2199
[4] Wang G, Cheng S Y, Wei W, et al. Characteristics and source apportionment of VOCs in the suburban area of Beijing, China[J]. Atmospheric Pollution Research, 2016, 7(4): 711-724. DOI:10.1016/j.apr.2016.03.006
[5] Liu B S, Liang D N, Yang J M, et al. Characterization and source apportionment of volatile organic compounds based on 1-year of observational data in Tianjin, China[J]. Environmental Pollution, 2016, 218: 757-769. DOI:10.1016/j.envpol.2016.07.072
[6] Song C B, Liu B S, Dai Q L, et al. Temperature dependence and source apportionment of volatile organic compounds (VOCs) at an urban site on the north China plain[J]. Atmospheric Environment, 2019, 207: 167-181. DOI:10.1016/j.atmosenv.2019.03.030
[7] Liu Y H, Wang H L, Jing S G, et al. Characteristics and sources of volatile organic compounds (VOCs) in Shanghai during summer: implications of regional transport[J]. Atmospheric Environment, 2019, 215. DOI:10.1016/j.atmosenv.2019.116902
[8] Hu B Y, Xu H, Deng J J, et al. Characteristics and source apportionment of volatile organic compounds for different functional zones in a coastal city of southeast China[J]. Aerosol and Air Quality Research, 2018, 18(11): 2840-2852. DOI:10.4209/aaqr.2018.04.0122
[9] Sun J, Shen Z X, Zhang Y, et al. Profiles and source apportionment of nonmethane volatile organic compounds in winter and summer in Xi'an, China, based on the hybrid environmental receptor model[J]. Advances in Atmospheric Sciences, 2021, 38(1): 116-131. DOI:10.1007/s00376-020-0153-0
[10] 中华人民共和国公安部. 2021年一季度新注册登记机动车966万辆[EB/OL]. https://www.mps.gov.cn/n7944517/n7944597/n7945888/c7829360/content.html, 2021-04-06.
[11] 河南省人民政府. 河南省车用乙醇汽油管理办法[EB/OL]. http://www.henan.gov.cn/zwgk/system/2007/01/18/010020281.shtml, 2004-11-17.
[12] 国家能源局. 关于扩大生物燃料乙醇生产和推广使用车用乙醇汽油的实施方案[EB/OL]. http://www.nea.gov.cn/2017-09/13/c_136606035.htm, 2017-09-13.
[13] Wu X, Zhang S J, Guo X, et al. Assessment of ethanol blended fuels for gasoline vehicles in China: fuel economy, regulated gaseous pollutants and particulate matter[J]. Environmental Pollution, 2019, 253: 731-740. DOI:10.1016/j.envpol.2019.07.045
[14] Zhu R C, Hu J N, He L Q, et al. Effects of ambient temperature on regulated gaseous and particulate emissions from gasoline-, E10- and M15-fueled vehicles[J]. Frontiers of Environmental Science & Engineering, 2021, 15(1). DOI:10.1007/s11783-020-1306-1
[15] Koupal J, Palacios C. Impact of new fuel specifications on vehicle emissions in Mexico[J]. Atmospheric Environment, 2019, 201: 41-49. DOI:10.1016/j.atmosenv.2018.12.028
[16] Zhu R C, Hu J N, Bao X F, et al. Investigation of tailpipe and evaporative emissions from China Ⅳ and Tier 2 passenger vehicles with different gasolines[J]. Transportation Research Part D: Transport and Environment, 2017, 50: 305-315. DOI:10.1016/j.trd.2016.10.027
[17] Tang J Y, Li Y J, Li X L, et al. Intermediate volatile organic compounds emissions from vehicles under real world conditions[J]. Science of the Total Environment, 2021, 788. DOI:10.1016/j.scitotenv.2021.147795
[18] 孙露娜, 刘妍, 赵静波, 等. 天津隧道机动车VOCs污染特征与排放因子[J]. 环境科学, 2019, 40(1): 104-113.
Sun L N, Liu Y, Zhao J B, et al. Pollution characteristics and emission factors of VOCs from vehicle emissions in the Tianjin tunnel[J]. Environmental Science, 2019, 40(1): 104-113.
[19] Zhang Y L, Yang W Q, Simpson I, et al. Decadal changes in emissions of volatile organic compounds (VOCs) from on-road vehicles with intensified automobile pollution control: case study in a busy urban tunnel in South China[J]. Environmental Pollution, 2018, 233: 806-819. DOI:10.1016/j.envpol.2017.10.133
[20] Deng C X, Jin Y J, Zhang M, et al. Emission characteristics of VOCs from on-road vehicles in an urban tunnel in eastern China and predictions for 2017-2026[J]. Aerosol and Air Quality Research, 2018, 18(12): 3025-3034. DOI:10.4209/aaqr.2018.07.0248
[21] Zhang Q J, Wu L, Fang X Z, et al. Emission factors of volatile organic compounds (VOCs) based on the detailed vehicle classification in a tunnel study[J]. Science of the Total Environment, 2018, 624: 878-886. DOI:10.1016/j.scitotenv.2017.12.171
[22] Feng M, Hu X, Zhou L, et al. Real-world vehicle volatile organic compound emissions and their source profile in Chengdu based on a roadside and tunnel study[J]. Atmosphere, 2021, 12(7). DOI:10.3390/atmos12070861
[23] Jin B Q, Zhu R C, Mei H, et al. Volatile organic compounds from a mixed fleet with numerous E10-fuelled vehicles in a tunnel study in China: emission characteristics, ozone formation and secondary organic aerosol formation[J]. Environmental Research, 2021, 200. DOI:10.1016/j.envres.2021.111463
[24] United States Environmental Protection Agency. EPA positive matrix factorization (PMF) 5.0 fundamentals and user guide[EB/OL]. https://www.epa.gov/sites/default/files/2015-02/documents/pmf_5.0_user_guide.pdf, 2014-04-01.
[25] Coulter C T. EPA-CMB8.2 Users Manual[EB/OL]. https://www.epa.gov/sites/default/files/2020-10/documents/epa-cmb82manual.pdf, 2004-12-01.
[26] Carter W P L. Updated maximum incremental reactivity scale and hydrocarbon bin reactivities for regulatory applications[EB/OL]. http://ww2.arb.ca.gov/sites/default/files/barcu/regact/2009/mir2009/mir10.pdf, 2010-01-28.
[27] Grosjean D, Seinfeld J H. Parameterization of the formation potential of secondary organic aerosols[J]. Atmospheric Environment (1967), 1989, 23(8): 1733-1747. DOI:10.1016/0004-6981(89)90058-9
[28] Grosjean D. In situ organic aerosol formation during a smog episode: estimated production and chemical functionality[J]. Atmospheric Environment. Part A. General Topics, 1992, 26(6): 953-963. DOI:10.1016/0960-1686(92)90027-I
[29] 赵金帅, 于世杰, 王楠, 等. 郑州市少数民族运动会期间O3及VOCs污染特征的演变和评估[J]. 环境科学, 2020, 41(10): 4436-4445.
Zhao J S, Yu S J, Wang N, et al. Evolution and evaluation of O3 and VOCs in Zhengzhou during the national traditional games of ethnic minorities period[J]. Environmental Science, 2020, 41(10): 4436-4445.
[30] Han L X, Chen L H, Li K W, et al. Source apportionment of volatile organic compounds (VOCs) during ozone polluted days in Hangzhou, China[J]. Atmosphere, 2019, 10(12). DOI:10.3390/atmos10120780
[31] 邵敏, 袁斌, 王鸣, 等. 挥发性有机物(VOCs)来源及其大气化学作用[M]. 北京: 科学出版社, 2020: 158-159.
[32] Wang M L, Li S Y, Zhu R C, et al. On-road tailpipe emission characteristics and ozone formation potentials of VOCs from gasoline, diesel and liquefied petroleum gas fueled vehicles[J]. Atmospheric Environment, 2020, 223. DOI:10.1016/j.atmosenv.2020.117294
[33] Guan Y N, Wang L, Wang S J, et al. Temporal variations and source apportionment of volatile organic compounds at an urban site in Shijiazhuang, China[J]. Journal of Environmental Sciences, 2020, 97: 25-34. DOI:10.1016/j.jes.2020.04.022
[34] Liu Y, Shao M, Fu L L, et al. Source profiles of volatile organic compounds (VOCs) measured in China: part Ⅰ[J]. Atmospheric Environment, 2008, 42(25): 6247-6260. DOI:10.1016/j.atmosenv.2008.01.070
[35] Gu Y, Liu B S, Li Y F, et al. Multi-scale volatile organic compound (VOC) source apportionment in Tianjin, China, using a receptor model coupled with 1-hr resolution data[J]. Environmental Pollution, 2020, 265. DOI:10.1016/j.envpol.2020.115023
[36] 徐晨曦, 陈军辉, 姜涛, 等. 成都市区夏季大气挥发性有机物污染特征及来源解析[J]. 环境科学, 2020, 41(12): 5316-5324.
Xu C X, Chen J H, Jiang T, et al. Characteristics and sources of atmospheric volatile organic compounds pollution in summer in Chengdu[J]. Environmental Science, 2020, 41(12): 5316-5324.
[37] 王成辉, 陈军辉, 韩丽, 等. 成都市城区大气VOCs季节污染特征及来源解析[J]. 环境科学, 2020, 41(9): 3951-3960.
Wang C H, Chen J H, Han L, et al. Seasonal pollution characteristics and analysis of the sources of atmospheric VOCs in Chengdu urban area[J]. Environmental Science, 2020, 41(9): 3951-3960.
[38] 张翼翔, 尹沙沙, 袁明浩, 等. 郑州市春季大气挥发性有机物污染特征及源解析[J]. 环境科学, 2019, 40(10): 4372-4381.
Zhang Y X, Yin S S, Yuan M H, et al. Characteristics and source apportionment of ambient VOCs in spring in Zhengzhou[J]. Environmental Science, 2019, 40(10): 4372-4381.
[39] 李一丹, 尹沙沙, 张瑞芹, 等. 郑州市某城区冬季不同污染水平大气VOCs特征及源解析[J]. 环境科学, 2020, 41(8): 3500-3510.
Li Y D, Yin S S, Zhang R Q, et al. Characteristics and source apportionment of VOCs at different pollution levels during the winter in an urban area in Zhengzhou[J]. Environmental Science, 2020, 41(8): 3500-3510.
[40] 叶露, 邰菁菁, 俞华明. 汽车工业区大气挥发性有机物(VOCs)变化特征及来源解析[J]. 环境科学, 2021, 42(2): 624-633.
Ye L, Tai Q Q, Yu H M. Characteristics and source apportionment of volatile organic compounds (VOCs) in the automobile industrial park of Shanghai[J]. Environmental Science, 2021, 42(2): 624-633.
[41] George I J, Hays M D, Herrington J S, et al. Effects of cold temperature and ethanol content on VOC emissions from light-duty gasoline vehicles[J]. Environmental Science & Technology, 2015, 49(21): 13067-13074.
[42] Shi G L, Li X, Feng Y C, et al. Combined source apportionment, using positive matrix factorization-chemical mass balance and principal component analysis/multiple linear regression-chemical mass balance models[J]. Atmospheric Environment, 2009, 43(18): 2929-2937. DOI:10.1016/j.atmosenv.2009.02.054
[43] Yue T T, Yue X, Chai F H, et al. Characteristics of volatile organic compounds (VOCs) from the evaporative emissions of modern passenger cars[J]. Atmospheric Environment, 2017, 151: 62-69. DOI:10.1016/j.atmosenv.2016.12.008
[44] Liu H, Man H Y, Cui H Y, et al. An updated emission inventory of vehicular VOCs and IVOCs in China[J]. Atmospheric Chemistry and Physics, 2017, 17(20): 12709-12724. DOI:10.5194/acp-17-12709-2017
[45] 高丹丹, 尹沙沙, 谷幸珂, 等. 河南省2016~2019年机动车大气污染物排放清单及特征[J]. 环境科学, 2021, 42(8): 3663-3675.
Gao D D, Yin S S, Gu X K, et al. Vehicle air pollutant emission inventory and characterization in Henan province from 2016 to 2019[J]. Environmental Science, 2021, 42(8): 3663-3675.