近几年，臭氧(O₃)已成为影响我国环境空气质量的重要因素，其中京津冀及周边地区、长三角地区以O₃为首要污染物的超标天数占比已经超过PM_2.5^[1-3]。研究表明，挥发性有机物（VOCs）可在紫外线照射下与氮氧化物（NO_x）发生光化学反应，产生光化学烟雾，光化学烟雾的主要成分为O₃^[4-5]。因此，作为O₃重要前体物的VOCs到研究学者的广泛关注。

VOCs种类繁多，不同种类的VOCs化学反应活性也不相同，研究VOCs的组成和来源特征对控制O₃污染和揭示复合型大气污染的形成都具有重要意义^[6]。目前国内关于VOCs的监测和研究主要集中在长江三角洲^[7-8]、珠江三角洲^[9-10]和京津冀^[11-12]等地区。山东半岛相关的研究较少，刘泽常等^[13]研究表明，济南市区VOCs的优势组分为C3—C5的烷烃、丙烯、顺-2-丁烯、间/对二甲苯和甲苯等，主要来源为汽车尾气、工业源和燃烧源。薛莲等^[14]发现青岛市大气VOCs中烯烃对臭氧的生成贡献远高于烷烃和芳香烃。张桢超^[15]发现威海市大气中，C2—C4烯烃类、烷烃类和苯系物对臭氧的生成贡献率较高，VOCs主要来源于机动车排放、工艺过程和溶剂使用。

泰安市地处山东省中部的泰山南麓，三面环山，属于内陆中小型城市。2016—2017年，泰安市O₃最大8 h平均浓度分别为197 μg·m⁻³和210 μg·m⁻³，在全省分别排名第二位和第一位。臭氧已成为泰安市夏、秋季节环境空气的首要污染物^[16]。了解臭氧前体物VOCs的污染现状及来源对泰安市采取适当措施改善空气质量具有重要意义。

本研究在泰安市城区建立一个观测站点，采用在线观测法，连续对站点大气中的VOCs进行监测，分析其浓度特征，并利用特征比值和模型分析对VOCs进行来源解析，同时评估其臭氧生成潜势，以期为泰安市大气环境VOCs和O₃污染管控提供科学支撑。

1.   实验与方法（Experiments and methods）

1.1   样品采集与分析

本次观测时间为2018年6月1日—7月11日，可以反应泰安市夏季大气中VOCs的污染特点。监测地点位于泰安市泰山区的山东电力高等专科学校校院内(36.18°N，117.11°E)，该观测点是泰安市的国控监测点，周边紧邻交通干线，同时分布着农贸市场、工业区、商业区和居民区，是典型的城市中心站点，观测点位置如图1所示。

图 1  监测点位图

Figure 1.  Location of the monitoring site

下载: 全尺寸图片 幻灯片

NMHCs的观测采用由中国科学院生态环境研究中心自主研发的GC-FID-VOCs在线监测仪24 h连续监测，采样时间分辨率为1 h，毛细管色谱柱型号为OV-1(30 m 柱长× 0.32 mm直径 × 1.0 μm厚度)；采样时，通过采样泵以50 mL·min⁻¹的流速将环境气体浓缩至温度为−80 ℃的吸附管中，然后升温加热至100 ℃进行热脱附，保持6 min；同时以5 mL·min⁻¹的N₂流速将解吸的样品吹入GC毛细管色谱柱中进行分离，此时将吸附管的温度升高至220 ℃，以60 mL·min⁻¹反吹10 min以清除残留；色谱柱的程序升温如下：初始温度为−60 ℃，保持3 min；以12 ℃·min⁻¹升温至−20 ℃；以6 ℃·min⁻¹升温至30 ℃；以10 ℃·min⁻¹升温至170 ℃，保持2 min；FID检测器的温度为250 ℃，仪器共检测到51种物质，其中包含27种烷烃、9种烯烃和15种芳香烃^[16]。醛酮类化合物(OVOCs)的观测采用涂有2，4-二硝基苯肼(DNPH)衍生化试剂的硅胶小柱采集，每2 h采集1个样品，并采用高效液相色谱(HPLC)方法检测，共检测出15种OVOCs。

为了保证观测数据的有效性和可靠性，GC-FID-VOCs在线监测仪每2 d采用美国环保署认可的Linde SPECTRA Environmental Gases标准气体进行5点校准，校准时相关系数均在0.992—0.995；同时，为了避免一些高反应性的VOCs物种的氧化损失，在采样吸附管的前端连接填充亚硫酸钠的捕集器，用于去除空气中的氧化剂，每2 d更换1次亚硫酸钠捕集器；高效液相色谱每2 d进行曲线校准，每20个样品分析一次校准曲线中间浓度点，每个目标化合物的测定结果与初始浓度值相对偏差≤30%^[16-17]。除VOCs的观测外，同时观测环境空气中的CO、SO₂和NO_x等参数，监测仪器均采用赛默飞世尔科技公司i系列的自动连续检测仪。

1.2   VOCs化学活性评价

不同城市中大气VOCs的来源各异，VOCs中各组分的浓度水平和化学活性也不同，对大气O₃生成的贡献也有差异。臭氧生成潜势(OFP)是用最大增量反应活性方法评估挥发性有机化合物的光化学反应性，并估算臭氧形成过程中单个有机化合物的贡献率^[18]，计算公式为：

式中，OFP_i为第i个VOCs物种的臭氧生成潜势，μg·m⁻³；[VOCs_i]表示物种i的环境质量浓度，μg·m⁻³；MIR_i为VOCs第i个物种最大增量反应中臭氧生成系数，可在文献[19]中查出。

1.3   正交矩阵因子模型(PMF)分析

正交矩阵因子分析模型(positive matrix factorization，PMF)作为受体模型，根据长时间序列的受体化学组分数据集进行VOCs来源解析^[20]。PMF计算过程中的基本公式为：

式中，X_ij为样本i中污染物j的浓度，×10⁻⁹；p表示污染源的数量；g_ik为第k个来源对第i个因子的贡献量，%；f_kj为第k个源中第j个组分的分布占比，%；e_ij为样本残差。PMF模型主要是将目标函数Q最小化^[21-22]，目标函数Q定义为：

式中，n为样本个数，m为物种个数；u_ij表示样本中物种的不确定性。根据PMF5.0指导方法要求，不确定度的计算公式为：

式中，unc. 表示样本中物种的不确定度；C表示样本中物种的浓度；RSD表示相对标准偏差；MDL表示检出限。

2.   结果与讨论（Results and discussion）

2.1   大气VOCs的浓度水平及组成特征

观测期间采样频率为1 h，VOCs浓度平均值为(16.57±7.99)×10⁻⁹（体积分数）。由表1可知，观测期间VOCs浓度水平最高的物种是甲醛(3.18±2.09)×10⁻⁹和丙酮(2.02±1.27)×10⁻⁹，其次为丙烷(1.71±1.41)×10⁻⁹、乙醛(1.39±0.61)×10⁻⁹和丁烷(0.92±0.88)×10⁻⁹。由图2可以看出，整个监测期间VOCs四大组分浓度顺序依次为：OVOCs(41.9%)> 烷烃(30.8%)> 芳香烃(19.5%) > 烯烃(7.8%)。

表 1  观测期间主要VOCs物种的浓度和OFP值

Table 1.  Concentrations and OFP of main VOCs species during the observation period

序号Serial number	平均浓度Average concentration			臭氧生成潜势OFP
	组分Species	数值(×10−9)Value		组分Species	数值/(μg·m−3)Value
1	甲醛	3.18±2.09		甲醛	40.34±26.52
2	丙酮	2.02±1.27		间/对-二甲苯	20.51±14.96
3	丙烷	1.71±1.41		乙醛	17.87±7.90
4	乙醛	1.39±0.61		异戊二烯	15.44±14.10
5	丁烷	0.92±0.88		间二乙基苯	14.59±11.71
6	甲苯	0.88±0.67		甲苯	14.52±10.96
7	异戊烷	0.81±0.43		戊醛	12.10±8.01
8	苯乙烯	0.72±1.13		正戊烯	11.90±7.36
9	戊醛	0.62±0.41		丁醛	10.74±2.46
10	苯	0.60±0.39		间-甲基苯甲醛	10.01±7.20
11	丁醛	0.56±0.13		邻-二甲苯	8.54±8.36
12	正戊烯	0.53±0.33		1,2,4-三甲基苯	7.86±4.77
13	异丁烷	0.50±0.44		1,3,5-三甲基苯	7.77±6.62
14	异戊二烯	0.48±0.44		对二乙基苯	7.76±5.27
15	丁烯醛	0.47±0.33		反-2-戊烯	6.15±3.47
16	乙苯	0.33±0.30		丙烯	5.83±5.59
17	十二烷	0.32±0.13		苯乙烯	5.77±9.12
18	戊烷	0.31±0.27		顺-2-丁烯	5.62±3.80
19	间/对-二甲苯	0.28±0.20		苯甲醛	5.28±2.75
20	丙醛	0.27±0.14		丙醛	5.02±2.57

 | Show Table

DownLoad: CSV

图 2  VOCs各组分浓度占比

Figure 2.  Concentration ratio of VOCs components

下载: 全尺寸图片 幻灯片

由于VOCs来源的不同和化学活性的差异，导致VOCs组分浓度的日变化特征也不同，分析VOCs浓度日变化特征是探讨其来源的重要手段之一。图3给出了观测站点大气中烷烃、烯烃、芳香烃和OVOCs的日变化趋势。烷烃、芳香烃和OVOCs日变化趋势较为一致，整体呈现夜间高白天低的变化特征；因为早晚存在较强的人类活动，如城市地区机动车尾气排放等，同时早晚大气较稳定，不利于VOCs的扩散；中午及下午对流强，边界层抬升，有利于污染物的扩散，同时中午和下午太阳辐射强，大气光化学反应活性剧烈，也会消耗一定量的VOCs，造成大气中VOCs的浓度下降。然而，烯烃在白天出现了显著抬升的变化特征，这主要是来自于植物活动排放异戊二烯的重要贡献^[23-24]，在白天随着太阳辐射的增强，异戊二烯排放量明显增加，在午后达到峰值，16:00以后由于植物活性等影响，其浓度快速下降。

图 3  VOCs各组分日变化特征

Figure 3.  Daily variation characteristics of VOCs components

下载: 全尺寸图片 幻灯片

对于烷烃，由于观测站点靠近蔬菜批发市场，运输车辆工作较早，受交通早高峰影响，在06:00出现峰值，随后太阳辐射增强，光化学反应消耗增加，浓度逐渐降低，在17:00以后，光化学反应消耗停止并随着市内交通晚高峰到来，污染物的浓度逐渐积累。芳香烃和OVOCs相对于烷烃和烯烃峰型规律没那么明显，芳香烃在中午出现了浓度升高的趋势，说明芳香烃类除来源于机动车排放外，还受溶剂挥发和化工企业排放的影响；OVOCs夜晚浓度较高，日出后出现微弱下降，随后缓慢抬升，在正午12:00左右发生快速下降达到最低，而后又快速抬升一直维持到深夜，说明除来自一次排放外，二次贡献生成和区域气象因素对其也有重要影响。

2.2   大气VOCs臭氧生成潜势特征

环境空气中VOCs组分的化学反应活性不同，对O₃的形成影响也不同，识别对大气O₃贡献较大VOCs物种，对于制定有效的减排控制措施意义重大^[25]。由表1可知，观测点观测期间甲醛(40.34±26.52) μg·m⁻³、间/对-二甲苯(20.51±14.96) μg·m⁻³、乙醛(17.87±7.90) μg·m⁻³、异戊二烯(15.44±14.10) μg·m⁻³和间二乙基苯(14.59±11.71) μg·m⁻³是OFP值最高的5种VOCs物种。由图4可知，观测期间大气VOCs四大类别对OFP贡献率顺序为：芳香烃(35.6%) > OVOCs(35.5%) > 烯烃(18.5%) > 烷烃(10.5%)。烷烃化合物VOCs浓度占比最高，但化学反应活性低，故对OFP的贡献较小^[26-27]；烯烃中异戊二烯的浓度水平较高，且所含的碳碳双键化学反应活性强，对OFP的贡献较大，异戊二烯主要来源于植物活动排放，说明植物源对大气中O₃的生成有重要贡献^[20]；芳香烃和OVOCs浓度水平和臭氧生成潜势都较高，因此，控制芳香烃和OVOCs的排放是未来控制泰安市臭氧污染的关键。

图 4  VOCs各组分对OFP的贡献

Figure 4.  Contribution of VOCs components to OFP

下载: 全尺寸图片 幻灯片

2.3   泰安市VOCs来源解析

2.3.1   典型物种特征比值分析

VOCs组分中，与·OH具有相似反应速率的的特征污染物之间的比值，可以反映其来源特征。戊烷在环境中主要来源于机动车尾气、天然气排放、燃料和液体汽油挥发，正戊烷(n-pentane)和异戊烷(i-pentane)具有相似的物理和化学性质，异戊烷/正戊烷的值可以初步判断其来源，当比值范围为0.82—0.89时，来源为天然气排放，比值范围为2.20—3.80时，来源为机动车尾气，比值范围在1.50—3.00时为液态汽油排放，比值范围在1.80—4.60时为燃料挥发，比值范围大于4.60时为其他源^[28]。图5(a)分别给出了异戊烷/正戊烷不同的比值线，可知观测点异戊烷和正戊烷的比值大部分分布在1.80—4.60之间，表明机动车尾气、液态汽油排放和燃料挥发都对其有贡献作用，而天然气排放对戊烷的来源贡献较少，其它源对其来源也有贡献。利用特征比值法分析可以看出异戊烷和正戊烷的比值较为分散，异戊烷和正戊烷的来源复杂，并不是单一来源，特征比值法无法很好的解析其来源。

图 5  典型物种特征比值

Figure 5.  Characteristic ratios of typical species

下载: 全尺寸图片 幻灯片

芳香烃是泰安市大气VOCs中对O₃生成贡献最大的物种，苯(benzene)与甲苯(toluene)的比值(B/T)常用来判断芳香烃在环境中的来源，当苯/甲苯值不大于0.20时，判断其来源为工业溶剂，机动车尾气源为0.50—0.60，燃煤源为1.50—2.20，生物质燃烧源比值约为2.50^[29]。由图5(b)分别给出了苯/甲苯的不同比值线，可以看出苯和甲苯比值分布在0.20—2.20之间，表明机动车尾气和燃煤源对其有贡献作用，经线性拟合的比值为0.56，与机动车尾气排放比值相近，进一步表明机动车尾气排放对观测点大气VOCs和O₃有重要贡献。B/T的值不只受到其来源的影响，还受到大气氧化性的影响，甲苯的光化学反应活性要大于苯，B/T的值相较于异戊烷/正戊烷也更为分散；另一方面B/T的值受风向影响较大，西南风向上有较高的甲苯浓度，表明污染物传输对其比值有影响。

综上，观测点位于城市中心点，VOCs的来源更为复杂，特征污染物的比值范围较为分散，简单的特征比值法无法准确的解析出污染物的来源，需要与PMF模型解析的结果进行相互印证。

2.3.2   PMF模型解析

针对观测点的VOCs数据，选取了来源指示性强和监测数据相对完整VOCs物种输入到PMF模型中，同时将观测的CO、NO_x和SO₂数据纳入计算，帮助识别排放源，最终共识别出6类因子，各类因子的源成分谱特征如图6所示。

图 6  观测期间大气VOCs源成分谱图

Figure 6.  Source composition spectrum of VOCs species during observation period

下载: 全尺寸图片 幻灯片

第一类因子中甲醛、乙醛和丙酮等OVOCs所占的百分比较高，因此将此类因子归为OVOCs源。第二类因子中，C2—C4烯烃和烷烃为优势组分，C2—C4烯烃和烷烃是LPG和溶剂挥发的关键物种^[12,30]，因此，将此类因子归为LPG和溶剂挥发源。第三类因子具有高组分的CO和苯，工业燃烧过程中可排放大量的CO，苯也是工业燃烧过程中重要的特征指示物^[31]，故将此类因子归为工业排放。第四类因子具有高组分的SO₂，电厂可排放大量的SO₂，因此，将此类因子归为电厂排放。第五类因子的优势组分为苯系物和烷烃，根据之前特征比值分析出异戊烷、正戊烷、苯和甲苯的主要来源都有机动车尾气排放，且C8—C10烷烃是柴油发动机排放尾气的标志^[32]，因此，将此类因子归为包括汽油车与柴油车的机动车尾气排放。第六类因子的优势组分为异戊二烯，城市中的异戊二烯大部分来源于植物活动排放，另外机动车尾气也会排放一部分排放，但该因子中与机动车尾气相关的其它VOCs物种贡献率都不高，因此将此类因子定义为植物排放源。

图7为PMF解析出6类VOCs排放源在观测期间对VOCs的相对贡献结果，可以看出，观测期间VOCs最大的排放源为LPG和溶剂挥发源(40.2%)，其次分别为OVOCs源(17.8%)、机动车排放(17.4%)、工业排放(11.8%)、植物源(10.5%)和电厂排放(2.3%)。因此，控制燃烧源和工业溶剂排放是控制泰安市夏季VOCs污染的重要途径。

图 7  观测期间不同源排放对VOCs的贡献率

Figure 7.  Contribution rate of emissions from different sources to VOCs during the observation period

下载: 全尺寸图片 幻灯片

3.   结论（Conclusion）

（1）观测期间VOCs浓度平均值为(16.57±7.99)×10⁻⁹，其中OVOCs占比最高为41.9%，烷烃占比为30.8%，芳香烃为19.5%，烯烃为7.8%。烷烃、芳香烃和OVOCs日变化趋势较为一致，整体呈现夜间高白天低的变化特征，而烯烃受到植物源排放异戊二烯的影响，在白天出现了显著抬升的变化特征。VOCs物种中平均浓度水平最高的前5种是甲醛、丙酮、丙烷、乙醛和丁烷。观测期间大气VOCs四大类别对OFP贡献率顺序为：芳香烃 > OVOCs > 烯烃 > 烷烃，其中甲醛、间/对-二甲苯、乙醛、异戊二烯和间二乙基苯是OFP水平最高的5种VOCs物种。

（2）观测点VOCs来源解析结果显示，观测期间泰安市VOCs最大的排放源为LPG和溶剂挥发(40.2%)，其次分别为OVOCs源(17.8%)、机动车排放(17.4%)、工业排放(11.8%)、植物源(10.5%)和电厂排放(2.3%)。因此，控制LPG和溶剂挥发是控制泰安市夏季VOCs和O₃污染的重要途径。