环境科学  2022, Vol. 43 Issue (1): 113-122   PDF    
拉萨市挥发性有机物的组成特征、季节变化和来源解析
余家燕1,2,3, 韩燕4, 陈木兰4, 张惠芳5, 陈阳4, 刘建国1     
1. 中国科学院安徽光学精密机械研究所环境光学与技术重点实验室, 合肥 230031;
2. 中国科学技术大学, 合肥 230026;
3. 重庆市生态环境监测中心, 重庆 401147;
4. 中国科学院重庆绿色智能技术研究院, 重庆 400714;
5. 西藏自治区生态环境监测中心, 拉萨 850031
摘要: 目前在世界范围内,由于高原城市海拔高、紫外线照射强烈,造成大气中挥发性有机物的源汇归趋呈现较独特的特点.本研究对2019年拉萨市2个城市站点和1个背景站进行大气挥发性有机物(VOCs)的离线罐采样及实验室组分分析,探究了拉萨地区VOCs体积分数水平、组成特征和来源贡献.拉萨市平均φ(VOCs)为49.83×10-9,其中烷烃类最高,占比为61%,其次为含氧VOCs(OVOCs,12%)、卤代烃(9%)、烯烃(9%)、芳香烃(5%)和炔烃(4%).拉萨市市区站点中,八廓街站和区辐射站的VOCs源贡献顺序为:燃烧源(分别为64%和48%,下同)>交通排放(17%和31%)>工业排放(14%和14%)>溶剂涂料(3%和3%)≈植物+背景(2%和4%).八廓街站燃烧源贡献大,与当地贡香焚烧和供暖排放有关;区辐射站的交通排放贡献较高,与距离交通枢纽近且上风向物流中心有关;工业排放对VOCs具有区域性影响.拉萨地区的VOCs体积分数、组分特征、源类贡献受气象和排放协同影响,呈现明显季节变化和地区差异.本研究可以提供我国典型高原城市的VOCs组成特征,可为高源城市臭氧形成和驱动因素的研究积累基础数据.
关键词: 挥发性有机污染物(VOCs)      正定矩阵因子分解法(PMF)      拉萨      组成特征      来源解析     
Characteristics and Source Apportionment of Ambient VOCs in Lhasa
YU Jia-yan1,2,3 , HAN Yan4 , CHEN Mu-lan4 , ZHANG Hui-fang5 , CHEN Yang4 , LIU Jian-guo1     
1. Key Laboratory of Environmental Optics and Technology, Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Hefei 230031, China;
2. University of Science and Technology of China, Hefei 230026, China;
3. Chongqing Eco-Environmental Monitoring Center, Chongqing 401147, China;
4. Chongqing Institute of Green and Intelligent Technology, Chinese Academy of Sciences, Chongqing 400714, China;
5. Tibet Eco-Environmental Monitoring Center, Lhasa 850031, China
Abstract: Due to the high altitude of plateau cities and strong ultraviolet radiation, the sources and fates of volatile organic compounds show unique characteristics. In this study, the atmospheric volatile organic compound (VOCs) samples were collected at two urban sites and one background site using tank sampling in Lhasa in 2019, and then the composition, concentration, and sources were characterized. The results showed that the average φ(VOCs) in Lhasa was 49.83×10-9, of which the proportion of alkanes was the highest (61%), followed by OVOCs (12%), halogenated hydrocarbons (9%), olefin (9%), aromatic hydrocarbons (5%), and alkynes (4%). The respective contributions of VOCs sources at urban sites, such as Barkhor Street and Radiation Station in Lhasa, were as follows: combustion (64% and 48%) > traffic emission (17% and 31%) > industrial emission (14% and 14%) > solvents and coatings (3% and 3%) ≈plant+background (2% and 4%). The contribution of combustion was large mostly due to local incense burning (especially at Barkhor Street) and heating emissions. Traffic emissions contributed about one third to the VOCs at Radiation Station, which is related to its proximity to the transportation hub and the storage and logistics center upwind. Industrial emissions have a regional impact on ambient VOCs. Under the synergistic influence of meteorology and emissions, VOCs concentration, composition characteristics, and source contribution showed obvious seasonal variations and site differences in the Lhasa area.
Key words: volatile organic compounds(VOCs)      positive matrix factorization(PMF)      Lhasa      component and characteristics      source apportionment     

挥发性有机物(volatile organic compounds, VOCs)是大气中普遍存在的痕量活性物质, 是生成近地面臭氧和二次有机气溶胶(secondary organic aerosol, SOA)的重要前体物[1, 2].VOCs的来源有天然源和人为源, 天然源主要来自植被排放和林草火灾等的排放[3, 4]; 人为源主要包括化石燃料燃烧、生物质燃烧、油料储存运输、溶剂与涂料的使用、石油精炼、化工生产和餐饮油烟等生产生活活动[5].研究VOCs的组成特征、时空变化和来源是探究臭氧形成机制的重要依据[6].

VOCs的化学组成十分复杂, 其组成受源排放、气象和光化学过程的共同影响.VOCs源排放具有显著的地域差异, 受人为活动和自然排放的共同影响, 且人为排放随着法规政策和控制措施的改变而呈动态变化[7~13].不同VOCs物种在大气中生成O3和SOA的能力存在显著差异[14, 15].因此, 准确理解一个地区VOCs的化学组成特征, 是评估其大气污染形成的潜力和健康效应的基础.这些特点都给VOCs的来源解析工作带来挑战[16].

城市是人类活动的密集区域, 其VOCs体积分数一般高于农村和背景点.整体来看, 城市环境下含氧VOCs(OVOCs)物种的体积分数水平远高于BVOCs和芳香烃, 含量最丰富的OVOCs物种为甲醛、乙醛和丙酮, 可占到美国城市总醛酮类化合物的48% ~64%[17].一般而言, 由于农村、背景点和高山等区域受天然源影响大.城市环境中, 我国北京、广州、成都、重庆和深圳等城市均呈现VOCs体积分数较高的特征, 主要优势种为OVOCs中的甲醛、乙醛和丙酮等物种[18~21].世界范围内, 高原城市地区的VOCs组成和演化特征目前所知尚有限.新西兰有研究指出, 其中央高原的VOCs以癸醛、苯甲醛和苯乙醛为主[22].

拉萨市位于西藏高原中部, 是典型的内陆高原城市, 其海拔3 650 m, 城区为谷地, 群山环绕, 气团流动易受山脉阻隔影响.拉萨全年日照时间在2 400~3 150 h, 年总辐射7 624 MJ ·m-2[23].产业结构上, 拉萨地区主要以第三产业为主, 工业产品主要以水泥、藏药等中成药、电、酒和自来水瓶装饮用水为主; 生活源排放方面, 存在秋冬季节取暖和宗教焚烧等习惯, 使得其VOCs排放来源与内陆城市有较大差异.

本研究的目的为探究典型高原城市拉萨地区的VOCs的全年组成特征, 分析当地自然排放、产业结构和宗教活动对排放的影响并识别其关键来源.本研究对VOCs进行全年观测和典型季节加密观测相结合, 结合正定矩阵因子分解法(positive matrix factorization, PMF)分析拉萨市VOCs的污染特征、识别其主要来源和对臭氧生成的影响.本研究可以增进对高原城市VOCs排放特征的认识, 探索高海拔和强太阳辐射对高原臭氧形成规律的认识.

1 材料与方法 1.1 观测概况

拉萨市位于中国西藏自治区东南部, 主导风向为偏西风.拉萨城区平原地带与周围山脉高差较大(最大可相差3 500 m), 不利于污染物扩散.根据拉萨市主导风向及城区各监测点位分布情况(图 1), 选择城市上风向点1个和城市混合区点位1个, 分别是区辐射站和八廓街.区辐射站(E 90.981°, N 29.656°)位于拉萨城区西部, 位于城市上风向和谷地边缘, 周围以山地和植被为主, 可代表城市周边的环境, 在一定程度上反映了城市生产生活排放VOCs的总体情况; 八廓街站(E 91.138°, N 29.651°)是城市O3高体积分数点, 邻近大昭寺等藏式建筑, 人口密集, 可综合反映城市交通、宗教文化和生产生活的排放VOCs污染水平, 以上两个点位于2019年进行逐月观测.同时, 为获取青藏高原背景VOCs组成信息, 在藏北高原的纳木错湖旁点位(E 90.961°, N 30.772°)进行两次VOCs采样(2019年6月19日的12:00和16:00), 作为背景点数据.

图 1 拉萨市环境空气质量监测站点分布示意 Fig. 1 Distribution of national air quality monitoring stations in Lhasa

1.2 采样方法和频次

为掌握拉萨市区VOCs在全年的组成特点, 在区辐射站和八廓街点位, 于2019年的1月至12月, 每月采集VOCs样品2个, 采样时间为09:00~10:00和15:00~16:00时, 采样方式为1 h累计采样.为了掌握臭氧体积分数水平较高时期的VOCs组成, 2019年5月25日至6月1日在八廓街站点开展了为期8 d的加密离线监测, 采样时间为09:00~19:00, 每隔2 h采1次样, 每次为1 h累积采样.在纳木错站点进行两次VOCs监测(2019年6月19日的12:00和16:00), 作为背景点数据.

采样方法按照文献[24]中规定的方法进行VOCs离线样品采集, 所用采样罐由美国Entech公司生产, 体积为3.2 L, 最大承受压力344 737 Pa, 采样流量设置为26.22 mL ·min-1.

使用前对采样罐反复清洗并抽真空, 降低其内表面的活性, 采样前开展真空度检验, 防止真空度泄漏导致采样过程样品受到污染, 每批清洗过的采样罐随机选择1~2个进行清洗空白检验.采样期间要求采样人员禁止使用香水、护手霜等护肤品/化妆品, 避开局地污染源, 采样完成后应按照采样记录表进行样品信息的标识和记录; 采样罐存储和运输均置于实验室定制的高密度海绵运输包, 空运回实验室进行分析.

1.3 VOCs测定方法

VOCs组分分析采用气相色谱质谱法, 参考文献[24, 25]的方法, 利用低温预浓缩技术对大气中的VOCs进行富集, 加热解析后通过GC-FID/MS进行分析测量[26].分析设备GC-MS/FID型号为安捷伦GC(7820)和MS(5977), 由气路系统、样品采集系统、冷冻捕集与热解析系统、GC-MS-FID分析系统和数据处理系统构成. 环境样品在采样泵的动力下被抽入超低温预浓缩系统, 在低温下被冷冻富集(-160℃)在捕集柱上; 两个通道的捕集柱分别被加热至120℃, 热解析后目标化合物在载气的吹扫下分别进入两路的分析柱和检测器, 气路I通过PLOT(Al2 O3/KCl)柱进行分离, 利用FID进行检测, 其检测的组分主要是C2~C5的碳氢化合物; 气路Ⅱ通过DB-624柱进行分离, 利用MS进行检测, 其检测组分主要是C5~C10的碳氢化合物、卤代烃、含氧有机物等.采样流量60 mL ·min-1, 色谱柱流速MS: 1.3 mL ·min-1; FID: 1.5 mL ·min-1, 共分析107种VOC化合物, 具体组分名称如表 1所示.

表 1 本研究的目标VOCs化合物1) Table 1 Target VOCs compounds for this study

1.4 质量保证/质量控制(QA/QC)

本研究中的QA/QC主要包括两部分: VOCs采样和组分分析.采样期间的QA/QC主要有: 采样罐清洁度检验、采样罐真空度检验、采样过程避免突发局地源; 采样后样品保存等方面.组分分析的QA/QC包括分析前、分析中和分析后三大环节, 分析前对检测系统进行MS调谐、系统验漏、流量校正、压力校正和标气测试等, 利用配气仪配制不同体积分数梯度(0.5×10-9、1×10-9, 2×10-9、4×10-9和8×10-9)的标准使用气体和4×10-9的内标使用气, 标准气体(PAMS、TO15和ISTD)为德国Linde公司, 可溯源至NIST; 测量前对同一体积分数标气开展重复测定, 计算其相对标准偏差获得精密度, 本研究各组分精密度在0.36% ~4.92%之间.

1.5 PMF模型

PMF是一种识别和解析受体点大气污染物不同来源及其贡献的数学模型.本研究基于2019年VOCs监测数据, 采用PMF模型对拉萨市两个城市站点开展来源解析[27].根据观测物种信噪比(S/N)、物种体积分数、示踪意义和活性大小筛选了33个拟合物种用于PMF分析, 包括11种烷烃(乙烷、丙烷)、5种烯烃、乙炔、8种芳香烃、4种卤代烃和4种OVOCs(表 1).

2 结果与讨论 2.1 拉萨市VOCs污染特征分析 2.1.1 常规监测期间VOCs体积分数及化学组成

表 2统计了观测期间(2019年全年)拉萨市VOCs总体积分数及各类组分平均体积分数, 与国内其他典型地区VOCs体积分数水平比较.拉萨市φ(VOCs)(49.83 ×10-9)高于北方城市郑州(30.66×10-9)和西南的成都(53.39 ×10-9), 与东南的杭州(56.72×10-9)相差不大, 但φ(烷烃)(30.46×10-9)高于这3个城市, φ(烯烃)(4.45 ×10-9)与杭州(4.45×10-9)接近, 低于成都(7.15×10-9)和南京(7.46 ×10-9), φ(芳香烃)(2.47×10-9)低于列出的5个城市, 卤代烃和OVOCs低于杭州.与同为高原站点的青海门源站相比, 拉萨市的φ(VOCs)是高原背景站(11.03×10-9)的4.5倍, 尤其是烷烃和卤代烃, 明显高于门源站(分别为6.46×10-9和0.19 ×10-9).

表 2 拉萨市与国内其他典型地区VOCs体积分数水平比较1)×10-9 Table 2 Comparison of VOCs concentration levels between Lhasa and other typical areas in China×10-9

拉萨市2019年3个站点的VOCs体积分数水平及化学组成见图 2, 八廓街φ(VOCs)(58.90 ×10-9)最高, 其次为区辐射站(40.80 ×10-9), 最后为纳木错背景点(10.10 ×10-9), 体积分数最低.总体上, 两个城市站点的化学组成相似度较高, 说明受到相似的VOCs排放源的影响.八廓街和区辐射站ω(烷烃)均最高, 分别为63.5%和57.8%, 其次为OVOCs, 质量分数相差不大, 在11% ~12%范围.两个城市站点烯烃和芳香烃体积分数相差不大, 区辐射站的卤代烃和炔烃较八廓街稍低.纳木错的ω(OVOCs)最高, 其次ω(卤代烃)较高, 分别为45.9%和24.9%, 由于部分卤代烃在环境空气中比较稳定, 所以纳木错站的OVOCs和卤代烃体积分数更接近城市站点.八廓街站体积分数排名前5的组分依次是丙烷、正丁烷、丙酮、异丁烷和乙烯, 其中丙烷体积分数远高于其他组分; 而区辐射站VOCs排名前5的组分与八廓街稍有不同, 虽然也是丙烷体积分数最高, 但高值并不突出, 其次是异丁烷、异戊烷、丙酮和正丁烷.鉴于以上组分通常来自机动车尾气和燃料过程的排放, 表明拉萨市典型区域的VOCs组成特征受人为排放源影响较大[33].具体来源可参见2.2节的来源解析.

图 2 2019年拉萨市3个站点VOCs体积分数水平及化学组成 Fig. 2 VOCs concentration and chemical composition at three sites in Lhasa during 2019

拉萨市VOCs体积分数及组成的日变化如图 3所示, 八廓街下午的VOCs体积分数较上午略有上升, 区辐射站则下降了23%, 从化学组成看, 八廓街的烷烃占比始终高于区辐射站, 烯烃质量分数始终低于区辐射站; 八廓街站点下午的烷烃较上午有上升, 主要升高物种为异丁烷、环戊烷、3-甲基己烷、甲苯和苯乙烯, 初步表明来自燃烧和溶剂等一次排放源在下午有增加, 烯烃则均有下降, 原因是烯烃光化学活性更高, 白天易被消耗所导致[12], 芳香烃和OVOCs质量分数均是下午高于上午.区辐射站趋势与八廓街正好相反, 说明同为城市点的两站点排放、扩散和光化学反应模式上都存在差异.

图 3 2019年拉萨市城市站点上午、下午VOCs体积分数水平及化学组成 Fig. 3 Concentration and chemical composition of VOCs in the morning and afternoon at urban sites in Lhasa during 2019

拉萨市两个站点VOCs逐月体积分数水平和化学组成存在明显的变化(图 4).八廓街φ(VOCs)曲线呈波动性变化, 介于(30.7~105.8)×10-9. 11月VOCs体积分数最高, 夏季5、6和7月体积分数较低, 这是因为夏季扩散条件好于冬季, 且无采暖排放影响.烷烃是最主要组分, 质量分数介于44.1% ~77.2%; 烯烃在6月质量分数最低(3.3%), 12月最高(12.4%); 芳香烃、卤代烃和OVOCs对体积分数的贡献分别介于1.5% ~10.2%、3.7% ~22.3%和3.2% ~22.9%.区辐射站VOCs体积分数在3月最高, 次高为4月, 6月最低.从化学组成来看, ω(烷烃)介于38.1% ~74.4%, 同样是各月质量分数最高的组分; ω(烯烃)最低为2.5%, 最高达到20.4%; OVOCs在5月质量分数最大(35.3%), 最低是在10月(6.3%); ω(芳香烃)和ω(卤代烃)分别介于2.2% ~10.6%和3.9% ~19.0%.除个别月份外, 偏冷季节(如1、2、11和12月)的VOCs中烯烃和炔烃的体积分数和质量分数要高于温暖季节(5~8月), 卤代烃和OVOCs的整体趋势则相反.

(a)八廓街站, (b)区辐射站 图 4 八廓街站和区辐射站VOCs逐月体积分数水平及化学组成 Fig. 4 Monthly VOCs concentration and chemical composition at Barkhor Street Station and Radiation Station

2.1.2 加密监测期间VOCs体积分数及化学组成

5月25日至6月1日加密监测期间, 小时φ(VOCs)范围在(12.90~237.90)×10-9之间, 波动较大, 平均值为53.2 ×10-9(如图 5).VOCs体积分数高值(>平均值)主要集中在09:00、11:00、17:00和19:00(各覆盖5个时次), 而13:00和15:00体积分数相对较低.VOCs呈现“早晚高, 白天低”的日变化规律.由于本研究缺乏探空资料, 故用已有报道的拉萨边界层高度进行替代分析[34], 夏季07:00边界层高度约为1.2~1.6 km, 10:00为1.2~1.6 km, 19:00为1.4~1.8 km.由此可知, 高原城市夏季(09:00~19:00)边界层高度变化不显著, 近地面VOCs体积分数主要受排放情况影响.由于前1d夜间排放的污染物得以积累, 而早晨因为上班出行高峰交通排放量大, 使得VOCs体积分数迅速增加.白天随着紫外辐射和湍流加强, 大气光化学反应程度加剧, 混合层高度也小辐攀升, 污染物水平和垂直扩散加强, VOCs体积分数逐渐下降; 傍晚光照减弱温度降低, 光化学反应逐步减弱, 并伴随交通晚高峰的影响, VOCs体积分数开始回升.

图 5 八廓街站2019年5月25日至6月1日的VOCs体积分数变化 Fig. 5 Variation of VOCs concentration at Barkhor Street Station from May 25 to June 1, 2019

加密观测期间各时次的VOCs化学组分占比如图 6(a)所示, 与每月常规监测比例相近, 烷烃是最主要的物种, 其质量分数介于56.9% ~66.7%之间, 平均质量分数为62.6%; 其次是OVOCs, 平均质量分数为17.4%, 09:00质量分数最低(14.2%), 19:00最高(24.0%), 最主要的OVOCs是丙酮; ω(卤代烃)、ω(芳香烃)、ω(烯烃)和ω(炔烃)的平均值分别为7.3%、6.0%、4.4%和1.9%, 其他组分比如二硫化碳仅占0.4%.排名前10的VOCs物种有7种都是烷烃[图 6(b)], 丙烷、异戊烷、异丁烷、正戊烷和正丁烷是重点关注组分, 其次丙酮体积分数也较高.具体来源分析详见2.2节.

(a)化学组分占比, (b)排名前物种 图 6 八廓街站2019年5月25日至6月1日的VOCs化学组成占比和主要组分 Fig. 6 VOCs chemical composition ratio and main components in Barkhor Street Station from May 25 to June 1, 2019

2.2 拉萨市2019年VOCs来源解析 2.2.1 PMF因子识别

本研究使用PMF模型对输入物种反复调试, 使解析得到Qrobust收敛, 且Qtrue/Qrobust小于1.5, 物种残差介于-3.0~3.0之间, 最终确定5个因子为解析结果, 分别为工业排放、溶剂涂料、交通排放、燃烧源和植物源(图 7).

①乙烷, ②丙烷, ③异丁烷, ④正丁烷, ⑤异戊烷, ⑥正戊烷, ⑦2, 2-二甲基丁烷, ⑧2, 3-二甲基丁烷, ⑨2-甲基戊烷, ⑩3-甲基戊烷, B11正己烷, B12乙烯, B13丙烯, B14 1-丁烯, B15 1, 3-丁二烯, B16异戊二烯, B17乙炔, B18苯, B19甲苯, B20乙苯, B21间/对-二甲苯, B22邻-二甲苯, B23苯乙烯, B24 1, 3, 5-三甲基苯, B25 1, 2, 4-三甲基苯, B26 1, 2-二氯丙烷, B27 1, 2-二氯乙烷, B28二氯甲烷, B29氯甲烷, B30丙酮, B31 2-丁酮, B32甲基叔丁基醚, B33乙酸乙酯 图 7 各因子中VOCs物种浓度及对各因子的贡献 Fig. 7 Concentrations and contributions of VOCs species for each other

因子1主要组分包含C2~C4烷烃、1-3丁二烯和二氯甲烷等卤代烃和丙酮, 其中与卤代烃相关的工业包含农药、医药、化学纤维、塑料、橡胶和日用化学品等[35, 36], 丙酮也是重要的有机合成原料, 广泛应用于塑料、橡胶、纤维和喷漆等行业[37], 因此将因子1识别为工业排放源.因子2主要包含乙苯、间/对-二甲苯和邻-二甲苯等芳香烃和乙酸乙酯, 涂料和溶剂使用过程中会排放大量的芳香烃, 而乙酸乙酯是许多类树脂的高效溶剂[38, 39], 广泛应用于油墨、人造革生产中, 将因子2识别为溶剂涂料相关的排放源.因子3主要是C4~C5烷烃, 这些组分多来自于机动车尾气, 比如异戊烷是机动车尾气和汽油挥发的主要成分之一[40, 41], 且因子3对甲基叔丁基醚(MTBE)的解释率很高, MTBE是一种汽油添加剂, 很多研究中都把MTBE作为汽油挥发源的示踪物[42, 43], 因子3可识别为交通排放源.因子4中低碳烷烃、氯甲烷、乙烯、丙稀、乙炔和苯等组分解释率较高, 具有燃烧源特征, 由于拉萨地区贡香焚烧现象较为突出, 本研究测定了当地寺庙中使用较多的4类贡香燃烧产生的主要VOCs组分的质量分数, 其主要成分分别为乙烷、丙烷、乙烯、丙烯、甲苯、丙酮、1-丁烯和丙烯醛等物种, 与因子4中物质有较好吻合, 因此因子4识别为燃烧源.因子5中主要物种为异戊二烯, 其主要来自植被排放[44], 识别为植物源.

2.2.2 八廓街VOCs来源及反应活性

拉萨市八廓街站2019年VOCs的PMF解析结果如图 8(a)所示, 燃烧源是八廓街站最主要的排放源(64%), 八廓街站方圆500 m内有近10家寺庙, 且距离拉萨著名景点大昭寺不到300 m, 受到贡香焚烧影响比较突出; 其次为交通排放和工业排放, 其贡献比例分别为17%和14%, 该点位地处繁华商业街, 交通排放量比较大, 且周边有多家制药厂, 北边有多家家具制造企业、印刷厂等, 因此这两类污染源的贡献也比较显著; 溶剂涂料和植物+背景的贡献较低, 仅占3%和2%.由于各排放源排放的VOCs物种不同, 其对臭氧生成潜势(ozone formation potential, OFP)的贡献也不一样.

图 8 八廓街站2019年VOCs来源及其对OFP的贡献 Fig. 8 VOCs sources at Barkhor Street Station and their contribution to OFP in 2019

基于最大增量反应活性方法(maximum incremental reactivity)[45, 46]分析了拉萨市八廓街站2019年各种来源的VOCs对OFP的贡献, 如图 8(b)所示.燃烧源对OFP的贡献依旧最大, 其次是交通排放, 工业排放对OFP贡献低于其对VOCs体积分数的贡献, 这与工业排放的VOCs组分多为低活性的卤代烃有关[47]; 溶剂涂料和植物+背景则相反, 这是由于溶剂涂料排放的VOCs多为芳香烃, 植物源主要排放异戊二烯, 均为高活性组分.

2.2.3 区辐射站VOCs来源及反应活性

区辐射站2019年VOCs的PMF解析结果如图 9(a)所示, 在该点位燃烧源仍然是VOCs的最大贡献源, 但明显低于在八廓街的贡献; 其次为交通排放(31%), 明显高于在八廓街的贡献, 区辐射站距离青藏公路和拉萨西站较近, 且上风向是当地的仓储物流中心, 交通运输繁忙, 尤其有较多货车运行, 因而贡献较大; 工业排放的贡献与对八廓街站的贡献相当, 该点位的上风向是当地的工业集中区, 有较多的建材、食品和制药企业, 对下风向的VOCs体积分数有一定贡献; 植物+背景和溶剂涂料分担率较低.从区辐射站2019年各类来源的VOCs对OFP的贡献见图 9(b), 燃烧源贡献了50%的OFP; 交通排放对OFP贡献略高于25%, 明显高于八廓街站; 工业排放、植物+背景和溶剂涂料对区辐射站的OFP贡献较为接近.

图 9 区辐射站2019年VOCs来源及其对OFP的贡献 Fig. 9 VOCs sources at Radiation Station and their contribution to OFP in 2019

3 结论

(1) 本文研究了2019年拉萨市典型区域的VOCs组成和时空分布特征.在城区, 烷烃质量分数均最高, 其次为OVOC; 而背景区域以OVOCs和卤代烃为主要化学物种.拉萨市区VOCs体积分数水平和化学组成存在明显的季节变化特征, 冬季最高, 主要受周边贡香焚烧、采暖燃烧的影响, 其次为春季和秋季, 夏季最低.在上风向地区, VOCs体积分数在春季最高, 夏季体积分数最低.夏季拉萨市VOCs日变化呈“早晚高, 白天低”的日变化规律, 主要受气象、光化学反应和排放的影响.

(2) 来源解析结果表明, 拉萨市大气VOCs主要有5大类源, 分别为燃烧源、交通排放、工业排放、溶剂涂料和植物源.其中燃烧源的贡献最大, 交通排放次之, 其三为工业排放, 各源类贡献季节变化特征明显.地域分布上, 市区主要受周边居民生活、商业活动、贡香焚烧的影响, 而上风向点受仓储物流中心的影响较为显著.

(3) 本研究可增进对我国西藏自治区典型高原城市VOCs组成和活性的科学认识, 为深入理解高原城市的VOCs排放现状和臭氧生成活性提供依据, 也可为管控臭氧污染提供依据.

致谢: 感谢中国科学院重庆绿色智能技术研究院三峡所大气研究中心团队对本研究的支持.感谢西藏自治区生态环境监测中心的王彩红和贡桑多吉为本研究开展的严谨和负责的采样工作, 在此一并致谢!

参考文献
[1] 陆思华, 白郁华, 张广山, 等. 大气中挥发性有机化合物(VOCs)的人为来源研究[J]. 环境科学学报, 2006, 26(5): 757-763.
Lu S H, Bai Y H, Zhang G S, et al. Source apportionment of anthropogenic emissions of volatile organic compounds[J]. Acta Scientiae Circumstantiae, 2006, 26(5): 757-763. DOI:10.3321/j.issn:0253-2468.2006.05.010
[2] Li J L, Zhang M G, Wu F K, et al. Assessment of the impacts of aromatic VOC emissions and yields of SOA on SOA concentrations with the air quality model RAMS-CMAQ[J]. Atmospheric Environment, 2017, 158: 105-115. DOI:10.1016/j.atmosenv.2017.03.035
[3] Wu F K, Yu Y, Sun J, et al. Characteristics, source apportionment and reactivity of ambient volatile organic compounds at Dinghu Mountain in Guangdong Province, China[J]. Science of the Total Environment, 2016, 548-549: 347-359. DOI:10.1016/j.scitotenv.2015.11.069
[4] 杨夏捷, 马远帆, 彭徐剑, 等. 南方林区2000~2016年林火释放污染物动态变化研究[J]. 中国环境科学, 2018, 38(12): 4687-4696.
Yang X J, Ma Y F, Peng X J, et al. Dynamic changes of pollutants released from forest fire in southern forested region during 2000-2016[J]. China Environmental Science, 2018, 38(12): 4687-4696. DOI:10.3969/j.issn.1000-6923.2018.12.037
[5] 代伶文, 孟晶, 李倩倩, 等. 长江经济带湖北省人为源VOCs排放清单及变化特征[J]. 环境科学, 2021, 42(3): 1039-1052.
Dai L W, Meng J, Li Q Q, et al. VOCs Emission inventory and variation characteristics of artificial sources in Hubei Province in the Yangtze River Economic Belt[J]. Environmental Science, 2021, 42(3): 1039-1052.
[6] Zheng J Y. Zhang L J, Che W W, et al. A highly resolved temporal and spatial air pollutant emission inventory for the Pearl River Delta region, China and its uncertainty assessment[J]. Atmospheric Environment, 2009, 43(32): 5112-5122. DOI:10.1016/j.atmosenv.2009.04.060
[7] 王俏丽. 浙江省沿海地区PM2.5和VOCs源成分谱构建及大气复合污染来源解析[D]. 杭州: 浙江大学, 2019.6-7.
[8] 苏榕, 陆克定, 余家燕, 等. 基于观测模型的重庆大气臭氧污染成因与来源解析[J]. 中国科学: 地球科学, 2018, 61(1): 102-112.
Su R, Lu K D, Yu J Y, et al. Exploration of the formation mechanism and source attribution of ambient ozone in Chongqing with an observation-based model[J]. Science China Earth Sciences, 2018, 61(1): 23-32.
[9] 王伶瑞, 李海燕, 陈程, 等. 长三角北部沿海城市2018年大气VOCs分布特征[J]. 环境科学学报, 2020, 40(4): 1385-1400.
Wang L R, Li H Y, Chen C, et al. Distributions of VOCs in a coastal city in the northern Yangtze River delta during 2018[J]. Acta Scientiae Circumstantiae, 2020, 40(4): 1385-1400.
[10] 刘郁葱. 大气排放源化学成分谱的对比分析及其对空气质量模拟的影响研究[D]. 广州: 华南理工大学, 2016.
[11] Wang B G, Shao M, Zhang Y H, et al. A study of volatile organic compounds and its emission factors from city vehicles[J]. Research of Environmental Sciences, 2006, 19(6): 75-80.
[12] 区家敏. 珠江三角洲VOCs排放来源识别、验证与基于反应活性的控制对策研究[D]. 广州: 华南理工大学, 2014.
[13] 刘金荣. 长三角城市群大气VOCs污染特征及来源解析[D]. 南京: 南京信息工程大学, 2018.
[14] Li J, Xie S D, Zeng L M, et al. Characterization of ambient volatile organic compounds and their sources in Beijing, before, during, and after Asia-Pacific Economic Cooperation China 2014[J]. Atmospheric Chemistry and Physics, 2015, 15(14): 7945-7959. DOI:10.5194/acp-15-7945-2015
[15] 胡崑. 长三角典型城市大气羰基化合物污染特征和来源解析[D]. 南京: 南京信息工程大学, 2018.
[16] 莫梓伟, 邵敏, 陆思华. 中国挥发性有机物(VOCs)排放源成分谱研究进展[J]. 环境科学学报, 2014, 34(9): 2179-2189.
Mo Z W, Shao M, Lu S H. Review on volatile organic compounds (VOCs) source profiles measured in China[J]. Acta Scientiae Circumstantiae, 2014, 34(9): 2179-2189.
[17] Grosjean E, Grosjean D, Fraser M P, et al. Air quality model evaluation data for organics. 2. C1-C14 carbonyls in Los Angeles air[J]. Environmental Science & Technology, 1996, 30(9): 2687-2703.
[18] 刘芮伶, 翟崇治, 李礼, 等. 重庆主城区夏秋季挥发性有机物(VOCs)浓度特征及来源研究[J]. 环境科学学报, 2017, 37(4): 1260-1267.
Liu R L, Zhai C Z, Li L, et al. Concentration characteristics and source analysis of ambient VOCs in summer and autumn in the urban area of Chongqing[J]. Acta Scientiae Circumstantiae, 2017, 37(4): 1260-1267.
[19] 刘芮伶, 黄晓锋, 何凌燕, 等. 质子转移反应质谱在线测量大气挥发性有机物及来源研究——以深圳夏季为例[J]. 环境科学学报, 2012, 32(10): 2540-2547.
Liu R L, Huang X F, He L Y, et al. On-line measurements and source analysis of VOCs in the atmosphere using proton transfer reaction-mass spectrometry: case study in summer of Shenzhen[J]. Acta Scientiae Circumstantiae, 2012, 32(10): 2540-2547.
[20] Liu Y, Yuan B, Li X, et al. Impact of pollution controls in Beijing on atmospheric oxygenated volatile organic compounds (OVOCs) during the 2008 Olympic Games: observation and modeling implications[J]. Atmospheric Chemistry and Physics, 2015, 15(6): 3045-3062. DOI:10.5194/acp-15-3045-2015
[21] Tan Q W, Liu H F, Xie S D, et al. Temporal and spatial distribution characteristics and source origins of volatile organic compounds in a megacity of Sichuan Basin, China[J]. Environmental Research, 2020, 185. DOI:10.1016/j.envres.2020.109478
[22] Effah E, Barrett D P, Peterson P G, et al. Herbivory and attenuated UV radiation affect volatile emissions of the invasive weed Calluna vulgaris[J]. Molecules, 2020, 25(14): 3200. DOI:10.3390/molecules25143200
[23] 赵地, 诺桑, 措加旺姆, 等. 西藏拉萨地区太阳紫外辐射观测[J]. 大气与环境光学学报, 2018, 13(2): 81-87.
Zhao D, Norsang G, Tsoja W, et al. Measurements of solar UV radiation in Lhasa, Tibet[J]. Journal of Atmospheric and Environmental Optics, 2018, 13(2): 81-87.
[24] HJ 759-2015, 环境空气挥发性有机物的测定罐采样/气相色谱-质谱法[S].
[25] US EPA. Analysis of volatile organic compounds in air contained in canisters by method TO-15[EB/OL]. https://www.epa.gov/sites/production/files/2018-01/documents/sop_hw_31_revision_6.pdf, 2014-06-05.
[26] 梁柳玲, 李青倩, 章洁文, 等. 苏玛罐采样-GC-FID/MS同时测定环境空气中多种VOCs[J]. 环境监测管理与技术, 2021, 33(1): 38-42.
Liang L L, Li Q Q, Zhang J W, et al. Simultaneous determination of VOCs in Air by SUMMA canister sampling-GC-FID/MS[J]. The Administration and Technique of Environmental Monitoring, 2021, 33(1): 38-42.
[27] U.S. Environmental Protection Agency Office of Research and Development. EPA Positive Matrix Factorization (PMF) 5.0 fundamentals and user guide[EB/OL]. https://www.epa.gov/sites/production/files/2015-02/documents/pmf_5.0_user_guide.pdf, 2014-08-01.
[28] 何丽. 成都城区大气挥发性有机物污染特征及来源研究[D]. 成都: 西南交通大学, 2018.
[29] 李晓茜. 重庆市北碚城区大气中VOCs的组成特征与长期变化规律研究[D]. 重庆: 西南大学, 2020.
[30] 杨笑笑, 汤莉莉, 张运江, 等. 南京夏季市区VOCs特征及O3生成潜势的相关性分析[J]. 环境科学, 2016, 37(2): 443-451.
Yang X X, Tang L L, Zhang Y J, et al. Correlation analysis between characteristics of VOCs and ozone formation potential in Summer in Nanjing urban district[J]. Environmental Science, 2016, 37(2): 443-451. DOI:10.3969/j.issn.1000-6923.2016.02.021
[31] 张翼翔, 尹沙沙, 袁明浩, 等. 郑州市春季大气挥发性有机物污染特征及源解析[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.
[32] 白阳. 青藏高原背景站VOCs特征及来源分析[D]. 太原: 山西大学, 2013.
[33] Zhang X F, Yin Y Y, Wen J H, et al. Characteristics, reactivity and source apportionment of ambient volatile organic compounds (VOCs) in a typical tourist city[J]. Atmospheric Environment, 2019, 215. DOI:10.1016/j.atmosenv.2019.116898
[34] 刘红燕, 苗曼倩. 青藏高原大气边界层特征初步分析[J]. 南京大学学报(自然科学), 2001, 37(3): 348-357.
Liu H Y, Miao M Q. Preliminary analysis on character isiics of boundary layer in Qinghai-Tibet Plaeaut[J]. Journal of Nanjing University (Natural Sciences), 2001, 37(3): 348-357.
[35] 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
[36] Mo Z W, Shao M, Lu S H. Compilation of a source profile database for hydrocarbon and OVOC emissions in China[J]. Atmospheric Environment, 2016, 143: 209-217. DOI:10.1016/j.atmosenv.2016.08.025
[37] 周子航, 邓也, 吴柯颖, 等. 成都市典型工艺过程源挥发性有机物源成分谱[J]. 环境科学, 2019, 40(9): 3949-3961.
Zhou Z H, Deng Y, Wu K Y, et al. Source profiles of VOCs associated with typical industrial processes in Chengdu[J]. Environmental Science, 2019, 40(9): 3949-3961.
[38] 王学臣, 王帅, 刘大喜, 等. 典型工业源挥发性有机物排放特征及臭氧生成潜势分析[J]. 环境污染与防治, 2020, 42(11): 1387-1391.
Wang X C, Wang S, Liu D X, et al. VOCs emission characteristics and ozone formation potential analysis of typical industrial sources[J]. Environmental Pollution and Control, 2020, 42(11): 1387-1391.
[39] 潘昕, 张巍, 黄银芝, 等. 典型涂料制造企业VOCs排放量核算与排放特征分析[J]. 环境工程学报, 2021, 15(3): 1049-1059.
Pan X, Zhang W, Huang Y Z, et al. Emissions accounting and characteristics of volatile organic compounds in typical paint manufacturing enterprises[J]. Chinese Journal of Environmental Engineering, 2021, 15(3): 1049-1059.
[40] 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
[41] Wang P, Ying Q, Zhang H L, et al. Source apportionment of secondary organic aerosol in China using a regional source-oriented chemical transport model and two emission inventories[J]. Environmental Pollution, 2018, 237: 756-766. DOI:10.1016/j.envpol.2017.10.122
[42] 韩丽, 陈军辉, 姜涛, 等. 成都市冬季挥发性有机物污染特征及来源研究[J]. 四川环境, 2019, 38(3): 120-125.
Han L, Chen J H, Jiang T, et al. Pollution characteristics and source apportionments of volatile organic compounds of Chengdu in winter[J]. Sichuan Environment, 2019, 38(3): 120-125.
[43] Hsieh L T, Yang H H, Chen H W. Ambient BTEX and MTBE in the neighborhoods of different industrial parks in southern Taiwan[J]. Journal of Hazardous Materials, 2006, 128(2-3): 106-115. DOI:10.1016/j.jhazmat.2005.08.001
[44] 高超, 张学磊, 修艾军, 等. 中国生物源挥发性有机物(BVOCs)时空排放特征研究[J]. 环境科学学报, 2019, 39(12): 4140-4151.
Gao C, Zhang X L, Xiu A J, et al. Spatiotemporal distribution of biogenic volatile organic compounds emissions in China[J]. Acta Scientiae Circumstantiae, 2019, 39(12): 4140-4151.
[45] An J L, Zhu B, Wang H L, et al. Characteristics and source apportionment of VOCs measured in an industrial area of Nanjing, Yangtze River Delta, China[J]. Atmospheric Environment, 2014, 97: 206-214. DOI:10.1016/j.atmosenv.2014.08.021
[46] 张俊刚, 王跃思, 王珊, 等. 京津地区大气中非甲烷烃(NMHCs)质量浓度水平和反应活性研究[J]. 环境科学研究, 2008, 21(5): 158-162.
Zhang J G, Wang Y S, Wang S, et al. Ambient mass concentration and reactivity of non-methane hydrocarbons (NMHCs) in Beijing and Tianjin cities[J]. Research of Environmental Sciences, 2008, 21(5): 158-162.
[47] 孙杰, 王跃思, 吴方堃, 等. 唐山市和北京市夏秋季节大气VOCs组成及浓度变化[J]. 环境科学, 2010, 31(7): 1438-1443.
Sun J, Wang Y S, Wu F K, et al. Concentration and change of VOCs in summer and autumn in Tangshan[J]. Environmental Science, 2010, 31(7): 1438-1443.