ACTA Scientiarum Naturalium Universitatis Pekinensis
天津市夏季臭氧污染特征及来源的模拟分析
洪礼楠 黄争超 秦墨梅 王雪松†
环境模拟与污染控制国家重点联合实验室, 北京大学环境科学与工程学院, 北京 100871; † 通信作者, E-mail: xswang@pku.edu.cn
利用MODELS-3/CMAQ模式系统和高阶去耦合直接技术(HDDM-3D), 对天津市2014年7月的臭氧(O3)污染进行模拟, 分析臭氧生成的前体物控制区分布规律, 量化天津市及周边地区排放的影响。研究结果表明,天津市臭氧浓度分布从中心城区向外部郊区逐渐增高, 东南部临海的滨海新区浓度水平最高。天津全市7月的臭氧生成以VOCS控制为主, 中部地区(中心城区、北辰区、东丽区、滨海新区) 95%以上天数的臭氧处于VOCS控制区, 北部远郊(蓟县、宝坻区、宁河区、武清区)以及中部津南区有2/3以上天数受VOCS控制, 南部郊区(西青区、静海区、大港区)的VOCS控制区与共同控制区出现的比例相近。山东省排放是天津市夏季臭氧的主要来源, 平均贡献占比约为1/4, 对天津东南沿海地区的影响尤为突出; 河北省排放平均贡献占比约为1/6, 主要影响天津西部区县的臭氧水平; 天津排放贡献则主要分布在北部的宝坻和蓟县。关键词 天津; 臭氧; 前体物控制区; 来源; HDDM中图分类号 X513
臭氧(O3)是对流层大气的重要组成部分, 在地球大气化学中起着非常重要的作用[1]。近年来, 随着城市化和工业化的不断发展, 在中国的一些大型城市群, 光化学污染事件屡屡发生, 臭氧污染成为严重的空气污染问题[2–3]。高浓度的臭氧会对人体健康和生态系统造成诸多不利影响[4]。对臭氧污染的控制需要掌握其来源及其生成对前体物排放变化的响应规律, 从而通过对污染源的合理规划与控制, 达到有效降低臭氧浓度水平的目的。空气质量模式常用于研究对流层臭氧的生成机制与来源, 主要的数值分析方法包括强力扰动法、去耦合直接法[5]和标记示踪技术[6]。其中, 高阶去耦合直接法(High-order Decoupled Direct Method, HDDM)能够计算目标污染物对不同污染源排放的一阶和二阶敏感性系数, 既可用于分析臭氧与前体物之间的敏感性响应, 也可带入泰勒展开公式来计算不同源排放对臭氧生产的贡献[7]。与传统的强力扰动法相比, 该方法能够有效地避免数值计算误差的影响, 同时具有较高的计算效率, 近年来越来越多地应用于东亚及北美的相关研究中[8–14]。
天津市位于华北平原东部, 东临渤海, 北倚燕山, 是环渤海经济中心, 拥有京津冀地区最多的工
[15]业企业 。臭氧污染是影响天津市夏季空气质量 的一个重要问题。2014 年, 天津市全年有 52 天首要污染物为臭氧, 占比为 15.2%, 且均出现在夏季(6—8 月)。目前, 对天津市臭氧浓度的时空分布以及与气象要素的关系等研究显示, 天津市臭氧受外来源传输影响较大[16–20]。但是, 对周边地区排放如何影响天津臭氧浓度以及天津地区臭氧形成机制的研究存在不足。本研究选择 2014 年 7 月天津地区的臭氧污染过程, 运用 Models-3/cmaq 模式系统和 HDDM-3D 技术, 对天津及周边地区的空气质量进行模拟, 分析天津市臭氧的生成对前体物排放变化的响应规律, 量化不同地区的排放影响, 以期对天津市夏季的臭氧污染控制提供科学参考。
1 资料与方法 1.1模型系统介绍
本文采用中尺度气象模式 WRF[21]、源排放处理模式 SMOKE[22]和化学传输模式 CMAQV5.0.2[23],建立适用于京津冀及周边地区的大气污染模式系统。模拟区域以(34.53°N, 108.92°E)为中心, 采用二重嵌套网格, 网格分辨率分别为 36 km 和 12 km,第二重网格覆盖京津冀及周边地区(图 1), 模式垂直方向分为 13 层。采用 SAPRC07 气相化学机理和 AERO6 气溶胶机理。在模拟区域内, 我国的源
排放清单来自清华大学建立的中国多尺度排放清单(MEIC), 数据年份为 2012 年, 数据网格分辨率为0.25 (www.meicmodel.org); 中国以外地区的排放数据使用 INTEX-B 项目建立的东亚地区 0.5°×0.5°分辨率的排放清单[24]; 此外, 运用 BEIS3 (Biogenic Emissions Inventory System)模式, 结合模拟区域内的植被与土地利用分布数据以及模拟期间的气温、辐射等气象场模拟结果, 对天然源的污染物排放量进行估算[25]。图 2 给出第二重模拟范围内 NOX 和VOCS 人为源排放的分布情况, 其中, 天津市的主要排放集中在天津市中心城区和滨海新区。
1.2 排放源分类和地理分区
为研究不同地区、不同类型污染源对天津市臭氧浓度的影响, 将源排放划分为电厂源、工业源、移动源、居民源、农业源和天然源六类, 同时将第二重模拟区域内的源排放分为 12 个地区: 天津、北京、河北省北部(包括张家口和承德)、河北省东部(包括唐山和秦皇岛)、河北省中部(包括廊坊、保定、石家庄、沧州、衡水)、河北省南部(包括邢台和邯郸)、内蒙古、辽宁、山西、山东、河南以及其他地区。
1.3 臭氧来源的计算方法
本文使用建立于 CMAQ 模式的 HDDM-3D 技术, 计算臭氧浓度对各地区不同类型源排放的敏感性系数。臭氧浓度 c 对于某类源排放 Eij (第 i 个地区第 j 类源排放)的一阶敏感性系数sij 和二阶敏感(1) (2)性系数 sij定义如下:
式中, Eij0 表示第 i 个地区第 j 类源的基础排放量。当该类源排放变化后(相对于基础排放的变化比例为 ∆ε 后), 将臭氧浓度 c 进行泰勒级数展开, 可表示为
式中, ccase代表源排放发生变化的情景浓度, c0代表基准情景浓度, H.O.T 代表高阶项。那么, 发生变化的这部分源排放对污染物浓度的影响大小(或称贡献)即为
需要说明的是, 本文计算的源贡献中包含某个地区某类源排放的全部污染物对臭氧浓度水平的总影响, 并不局限于臭氧前体物(NOX 和 VOCS)的作用。Itahashi 等[12]在对东亚地区的研究中比较了强
力扰动法和 HDDM 方法对源贡献计算的结果差别,认为对于中国的高排放地区, 当源排放变化比例过大时, 式(3)中的高阶项(即 H.O.T)对浓度的影响增大(即高阶的非线性作用对臭氧浓度的影响不可忽略), 导致利用式(4)计算得到的源贡献与强力扰动法的结果偏差较大, 但源排放改变程度在 50%以内时, 两种方法的结果较为一致。鉴于此, 同时考虑到实际工作中源排放控制的可行性, 本文将源排放变化系数(∆ε)设置为−50% (即计算各类源在削减排放量 50%的条件下对臭氧浓度的贡献)。某个地区某类源对天津市臭氧浓度贡献率可表示为
式中, m 为排放地区数, n 为源类别数, Xij 表示第 i个地区第 j 类源对天津市臭氧污染的贡献比例[26]。
1.4 臭氧前体物控制区的划分
使用 HDDM-3D 方法计算出臭氧对第二重模拟区域内全部 NOX 和 VOCS 排放的敏感性系数, 利用式(4)计算 NOX 和 VOCS 排放分别削减 50%条件下的臭氧浓度改变量, 根据浓度下降程度判断臭氧生成的控制区类型[27–28]。1) NOX 控制区: NOX 减排使得 O3下降大于 5 μg/m3, 且大于 VOCS 减排引起 O3下降量两倍以上的地区; 2) VOCS 控制区: VOCS 减排使得 O3下降大于 5 μg/m3, 且大于 NOX 减排引起O3下降量两倍以上(或 NOX 减排引起 O3上升)的地区; 3) 共同控制区: NOX 和 VOCS 减排均能使 O3下降 5 μg/m3 以上, 且二者差距在两倍以内的地区; 4) NOX 滴定区: NOX 减排使得 O3 升高 5 μg/m3 以上,且 VOCS 减排使得 O3降低不足 5 μg/m3的地区; 5)非控制区: NOX 和 VOCS 减排使得 O3的变化均在 5 μg/m3以内的地区。
2 结果与讨论
2.1模拟效果评估
臭氧浓度模拟效果评估的观测数据来源于中国环境监测总站在全国城市空气质量实时发布平台(http://113.108.142.147:20035/emcpublish/)发布的污染物小时浓度数据。使用标准化平均偏差(NMB)、标准化平均误差(NME)以及相关系数(COR), 评估模拟结果与观测数据在浓度水平和变化趋势上的吻合程度。
表 1 给出天津市 15 个国控站点臭氧小时浓度模拟值与观测值对比的平均统计指标。虽然不同站点的模拟效果存在一定差异, 但就总体而言, 本文与国内外一些相关研究报道的模拟效果相近, 在可接受范围内。图 3 是天津市国控站点的平均臭氧模拟浓度与观测值对比的时间序列, 可以看到, 模拟得到的 O3浓度和变化趋势与观测数值基本上一致,峰值出现时间也较为吻合。
2.2 臭氧污染特征
根据天津市 2014 年 7 月臭氧最大 8 小时观测浓度, 计算该月超过国家空气质量二级标准的污染日中 10:00—17:00 臭氧高值时段模拟浓度空间分布的平均结果, 如图 4 所示。在 7 月的污染日中, 天津市中心城区臭氧平均浓度为 150 μg/m3 左右; 天津近郊区中, 北部(北辰、东丽)臭氧平均浓度(约167 μg/m3)低于南部(西青、津南); 在远郊区中, 宝坻、武清、宁河、静海、大港臭氧平均浓度均高于近郊地区, 西南部的静海浓度水平相对略低, 而宝坻浓度较高, 北部的蓟县浓度与天津市近郊相当;东南沿海的滨海新区平均浓度>200 μg/m3, 为天津市最高值。总体而言, 天津市臭氧分布呈现中心城区最低, 沿市区向外逐渐增高的环形分布, 与刘彩
[20]霞等 在天津市夏季臭氧的监测实验研究结果一致。
图 3 显示, 天津市 2014 年 7月经历了上旬、中旬、下旬 3 次光化学烟雾污染过程, 臭氧浓度为中度到重度污染水平。图 5 给出 7 月天津中心城区气象站点的风玫瑰图(站点位置如图 1(b)所示, 数据来自 http://rda.ucar.edu/datasets/ds461.0), 其主导风向以偏南风为主, 4 m/s以下的低风速占很大比例。图6 给出 7 月部分污染日 14:00 近地面臭氧浓度及风 场的模拟结果。7 日京津冀区域整体大范围南风,天津市北部地区臭氧浓度较高, 南部郊区和滨海新区浓度相对较低; 17 日天津市渤海海面东南风, 河北省中部地区静风, 市内东南风为主, 臭氧污染高值集中在天津市北部的宝坻、北辰等地, 南部静海浓度较低; 26 日天津市以西南风为主, 天津市区及南部的静海、津南等地臭氧浓度较低, 中部和北部地区臭氧浓度较高; 29 日天津市整体南风, 渤海海面东南风, 臭氧高值出现在宝坻、宁河和滨海新区。总体而言, 当天津市出现南、西南、东南风向时易发生臭氧污染, 这与 Han 等[34]对天津夏季的臭氧观测的研究结果基本上一致。
2.3 臭氧生成控制区分布
按照 1.4 节的方法, 对 2014 年 7月天津污染日10:00—17:00 的臭氧平均浓度进行前体物敏感性分析。图 7为分别削减第二重模拟区域内 50%的NOX 和 VOCS 排放情景下臭氧浓度下降情况。当NOX 削减 50%后, 市中心城区和滨海新区臭氧平均浓度出现不同程度的上升。这是因为两地工业企业很多, 导致NOX排放较大, 大量的NO还原臭氧,同时高浓度NO2中止OH自由基在大气中的氧化循环(OH+NO2→HNO3), 因此, 降低 NOX 排放会带来O3 浓度的上升(即NOX对O3的滴定作用)。在其他区域, NOX排放下降50%, 臭氧浓度变化在
20 μg/m3 以内。当 VOCS 下降 50%, 滨海新区的臭氧会相应下降 50~80 μg/m3, 在天津市区及周边郊区也会下降 20~30 μg/m3。根据 1.4 节臭氧前体物控制区分区定义, 计算天津市各区县 2014 年 7 月每日的前体物控制区类别, 表 2 统计了模拟期间各类臭氧前体物控制区出现的概率。天津市中心城区、北部近郊区和滨海新区均以 VOCS 控制区为主, 南部近郊区和远郊以 VOCS 控制区和共同控制区为主, 蓟县臭氧则存在一定比例(10%)的 Nox 控制情况, 其余日期为 VOCS 控制或两类前体物共同控制。
2.4 各地区排放的影响
选取天津市环境监测中心(中心城区站点)、蓟县环保局(北部郊区站点)、团泊洼(南部郊区站点)和中新生态城(东部沿海站点)4 个代表性站点(如图1(b)所示), 分析天津及周边各省市排放对上述站点O3浓度的影响。根据 1.3 节的计算方法, 表 3 给出4 个站点 O3浓度受到不同地区排放的贡献情况。市环境监测中心臭氧主要来源于山东省(26%)、河北省(24%)和河南省(7%)的排放影响, 而天津源排放对臭氧水平具有抑制作用, 这与图 8 所示的在中心城区 NOX 排放对臭氧生成具有滴定作用是一致
的。蓟县环保局站点受主导风向南风影响, 天津市本地排放对该站点的贡献占 8%, 山东省贡献 15%,河北省贡献 13%。团泊洼站点主要受山东(17%)、河北(18%)和河南(5%)排放影响。中心生态城临近渤海湾, 在海洋东南风传输作用下, 受山东省源排放影响相对较高(35%), 相比其他站点, 河北省贡献 比例(12%)有所下降。总体而言, 虽然天津不同地区的臭氧来源存在差别, 但河北和山东两省的排放在污染日高值时段影响最为突出, 天津本地排放对除中心城区以外大部分地区的贡献也占有重要份额, 北京的排放对天津各站点影响不大, 平均接近3%。此外, 跨区域背景的贡献约占一半的比例, 反
映出 O3环境背景和跨区域传输的重要作用。
图 8 给出天津市、河北省中部、山东省 50%源排放扰动对天津市臭氧浓度影响的分布情况。天津市源排放对于市区及周边市郊的小部分区域以及滨海新区部分地区的 O3水平具有抑制作用。在 7月份南风的主导风向影响下, 天津市排放对其北部的宝坻、蓟县贡献相对较高。在南风作用下,天津市北部的北京、河北东部、河北北部的影响相对较小, 河北中部排放主要影响到天津市武清区、静海区、西青区等西部郊县地区(平均 7~12 μg/m3 的贡献)。山东省排放主要影响天津东部地区, 对滨海新区有约 20 μg/m3的臭氧贡献, 对天津市区和城市近郊地区影响在 10~15 μg/m3 范围内,对天津北部(蓟县、宝坻)和西部(武清、静海)的臭氧影响有所减小。
山东省是对天津市臭氧浓度影响最大的地区,图 9 给出山东省各类源排放对天津市环境监测中心站点的贡献比例(其余站点结果相似)。其中, 来自工业源的影响最大, 占 60%, 移动源、居民源和天然源的贡献相近(约 10%), 电厂源的影响则相对较小, 农业源对臭氧影响可忽略不计。
3 结论
本文运用 MODEL-3/CMAQ 模式和高阶去耦合直接技术(HDDM-3D), 模拟研究天津市 2014 年 7月的臭氧污染特征以及前体物控制区分布, 量化了 天津及周边地区源排放的影响。天津市的臭氧污染多在偏南风传输条件下发生。在污染日的臭氧高值时段, 中心城区的臭氧浓度最低, 从近郊区到远郊臭氧浓度逐渐升高, 近郊区中南部郊区高于北部郊区, 远郊区中临近渤海湾的滨海新区浓度最高, 北部蓟县臭氧浓度与近郊区水平相当。
除天津西南静海的臭氧生成主要受到 NOX 和VOCS 排放的共同控制外, 天津市大部分地区主要处于 VOCS 控制区, 其中中心城区、北辰区和东丽区以及滨海新区基本上均为 VOCS 控制区, 北部远
郊区及津南区的臭氧约 2/3 的天数受 VOCS 控制、1/3 的天数受 NOX 和 VOCS 共同控制, 南部郊区VOCS控制区和共同控制区出现比例相近。天津市的夏季臭氧污染主要受山东省、河北省以及天津本地的排放影响。山东省排放的平均贡献约占 1/4, 尤其对天津东南沿海地区影响最为突出, 山东省贡献中有近 2/3 来自工业活动的排放; 河北省排放对天津臭氧的贡献约占 1/6, 主要影响天津西部的区县; 天津排放对中心城区和滨海新区部分地区的臭氧因 NOX 滴定效应而体现出抑制作用, 其对臭氧的生成贡献主要分布于天津市北部郊区(宝坻、蓟县)。 致谢 感谢清华大学 MEIC 团队提供全国污染源排放清单。