2020, Vol. 12 
2. 河北省生态环境应急与重污染天气预警中心, 河北 石家庄 050051;
3. 中国人民解放军第3302工厂，河北 石家庄 050051
2. Hebei Environmental Emergency and Heavy Pollution Weather Warning Center, Shijiazhuang, Hebei 050051, China;
3. CPLA 3302 Factory, Shijiazhuang, Hebei 050051, China
空气质量模型是评估特定污染源对空气质量影响的重要工具。当前国际上典型的空气质量模式主要包括：哥大德地球观测系统-大气化学传输模型(GEOS-CHEM)这样全球尺度空气质量模型，嵌套网格空气质量预报模式(NAQPMS)、三维欧拉型区域空气质量模式(CAMx)、完全耦合气象模式和化学模式的区域空气质量模式(WRF-CHEM)、第三代空气质量预报和评估系统(CMAQ)等综合型区域尺度模型以及大气污染物扩散模型(ISC3)、环安大气影响评价系统(AERMOD)、非定常三维拉格朗日烟团输送模式(CALPUFF)等中小尺度模型[1]。薛文博等[2]应用CMAQ模型, 基于2008年排放清单分析了电力行业二氧化硫(SO2)、二氧化氮(NO2)、烟粉尘协同控制情境下硫沉降、氮沉降及细颗粒物(PM2.5)污染的改善效果。伯鑫等[3]利用CAMx模型分析了重点煤电基地大气污染物扩散对京津冀的影响。但CMAQ以及CAMx这类区域尺度模型比较适合模拟区域尺度的大气问题。针对中小尺度的大气问题则适宜采用中小尺度模型来模拟。国内外学者先后采用ISC3、CALPUFF、致酸污染物长距离传输模型(ATMOS)等空气质量模型，模拟了电厂对局部地区空气中SO2、NO2、烟尘、汞、酸沉降等污染指标的影响[4-12]。
河北省是目前我国空气污染最严重的地区之一[13]。据生态环境部公布的城市排名，石家庄、保定、邢台、邯郸、唐山等城市均位于全国空气质量较差的城市之列[14]。河北省偏煤炭的能源结构、偏重化工的产业结构和偏公路运输的交通运输结构造成区域内大气污染物排放超过环境容量，是该省空气质量较差的根本原因，其中电力行业，特别是煤电对污染物的贡献依然较大[3, 15]。现选取河北省南部某热电厂作为研究对象，利用热电厂在线排放数据，基于高斯烟羽模型的输入频率分布模型(Immision Frequency Distribution Model，IFDM)系统，对该电厂在2019年4—6月对周边地区及东北方向6 km的2个省控监测站点(S1和S2)空气质量影响进行模拟评估，以期为该地区空气质量的精细化管控提供技术支持。
1 IFDM系统介绍 1.1 方法原理IFDM是一款基于高斯烟羽扩散模型的精细化管控系统，设计用于在局地尺度上模拟反应活性较弱的污染物的扩散。根据布特林克(Bultynck)和马莱特(Malet)公式，高斯弥散参数取决于大气的稳定性和风速，基于体积理查森(Richardson)参数[16]。该模型不仅可以模拟点源对环境空气质量的影响，还可以模拟道路交通等线源及居民生活、工业区等面源对环境空气质量的影响[17]。线源在文卡特拉姆(Venkatram)和霍斯特(Horst)(2006年)中进行了处理[18]，通过将线源划分为一系列点源，并针对这一系列点源的数值积分。而面源则被视为一组垂直于风的等效平行线源，使得其不仅能够模拟点源对空气质量的影响，同时还可以模拟线源以及面源对空气质量的影响，使人们能够直接了解到不同污染源到空气质量的贡献情况。
模型基于Bultynck等[16]基本方程(1972年)。该方程直接适用于点源，但需要更改参考系统。新参考系统的原点位于污染源处，x轴与水平风矢量平行，y轴在垂直于x轴的水平面上。z轴向上。其模型见式(1)。
| $ C = \frac{Q}{{2\pi {\sigma _y}{\sigma _z}u}}{e^{\frac{{ - y2}}{{2\sigma _y^2}}}}{e^{\frac{{ - {{\left( {z - h} \right)}^2}}}{{2\sigma _z^2}}}} $ | (1) |
式中：C——质量浓度，g/m3；Q——排放量，g/s；h——有效烟羽抬升高度，m；u——有效烟羽抬升高度的风速，m/s；σy，σz——水平和垂直扩散参数，m；y——距烟羽中心线的距离，m；z——空间点高度，m。
1.2 线源处理IFDM系统也可以模拟线源的扩散，主要是将线源作为点源的一个整体来处理。具体有2种不同的情况，一是平行或几乎平行于风向的线源，在这种情况下，不可能用解析方法求解积分，可以用辛普森规则来确定；二是线源与风矢量之间的角度足够大，可以用解析方法求解积分，该方法比辛普森规则要快得多。
1.3 面源处理IFDM系统也可以模拟面源的扩散，对于面源的处理也有2种方法：一是将面源作为一个大的点源来处理，如果只需要确定距离污染源较远位置的浓度值，在这种理想情况下可以将面源作为一个大的点源来处理。这种方法几乎不增加计算时间，速度非常快；二是将面源作为线源的集合来处理，这种情况非常适合更详细地描述面源，然后将面源视为线源的集合。针对一个面源将创建多个线源，并且所有线源都垂直于风向，面源的排放分布在这些线源上，受体点的浓度则由单个线源影响的总和确定。由于线源都垂直于风向，可按上述线源的第二种方法进行处理，得出可接受的计算时间。
1.4 与区域尺度模型的耦合IFDM系统可以与区域尺度模型进行耦合，区域尺度模型可以为IFDM系统的输出提供背景浓度，具体方法见图 1。图 1右下角G中黑色网格为区域尺度模型的网格，红色网格为IFDM系统的网格，在该示例中IFDM系统的网格为规则网格。绿色的线为想要研究的道路源，在区域尺度模型和IFDM系统中都有该道路源的排放量。
|
图 1 区域尺度模型与IFDM系统的耦合方法 |
区域尺度模型的模拟结果见图 1 A所示，没有修正过的IFDM模拟结果见图 1 B所示。耦合方法：首先将IFDM模拟结果插值到区域尺度模型的网格当中，见图 1 C所示。然后用来修正区域尺度模型的模拟结果D=A-C，此时的区域尺度模拟结果为不含该道路排放源的结果。然后再将该结果插值到IFDM系统的网格当中，见图 1 E所示。最后，将IFDM的模拟结果与上述E结果相加，即为IFDM的最终模拟结果(见图 1 F)。通过这种方法可以避免重复计算，同时模拟的结果可以包含本地污染源贡献值以及区域背景值[19]。
1.5 基于IFDM系统的精细化管理技术通过上述针对线源和面源处理的IFDM系统可以模拟局地尺度所有污染源的扩散情况，同时还可以耦合区域尺度模型。基于该模型开发了地区精细化管理技术。该技术可对某个区域未来时段的污染来源进行解析，包括区域背景的贡献以及当地重点企业、重点道路和生活源对任意受体点的贡献浓度。当地管理部门可以通过该技术掌握所管辖地区不同污染源的贡献值和贡献比例，特别是重污染天气时可以做到有的放矢，精准治理。同时该技术还可以模拟不同减排情景的减排效果，为地方政府制定限产减排措施，提供科学的决策支持工具。
1.6 IFDM系统的验证文献[19-21]对进行以上技术处理的IFDM系统做过验证工作，表明IFDM系统具有良好的模拟性能。现获取河北某城市的污染源排放清单以及气象数据，输入该IFDM系统中模拟2015年11月15日—2016年3月15日采暖季的ρ(PM2.5)，验证表明，IFDM系统具有良好的模拟性能。图 2为IFDM模拟值与实际监测值的时间序列对比；图 3为IFDM模拟值与实际监测值对比的散点图，对模拟值与实际监测值进行了线性回归拟合，发现IFDM模拟值与实际监测值的相关系数为0.76。
|
图 2 河北某城市某站点2015年采暖季IFDM系统对PM2.5的模拟值与监测值对比 |
|
图 3 河北某城市2015年采暖季IFDM系统对PM2.5模拟值与监测值对比散点图及拟合度检验 |
河北省南部某热电厂所处地区第二季度的主导风向为西南风[22]，利用IFDM系统模拟了2019年4月1日—6月30日该热电厂对所在地区空气质量的影响，模拟区域覆盖所在地区整个区域，其中设置了包含2个省控监测站点在内的4 539个受体点。
S1和S2站点的实时污染物监测数据来自于河北省空气质量实时发布平台，风向、风速及温度等气象数据来自于该地区的气象地面观测站；该电厂分为2个机组，热电厂的烟囱信息及排放量数据均来自于该企业的烟气在线监测数据，见表 1。
|
表 1 热电厂烟囱信息及排放量数据 |
基于模拟分析结果计算了该热电厂对所在地区空气质量的月均贡献情况，提取了该热电厂对S1、S2 2个省控监测站点各项污染物的月均贡献值。图 4(a)(b)(c)为该热电厂2019年4—6月对所在地区SO2的影响均值。
|
图 4 2019年4—6月该热电厂对所在地SO2的影响均值 |
因4月份该地区主导风向为南风和东北风，该热电厂对其北部地区和西南部地区空气质量影响较大；5月份该地区的主导风向为南风，该热电厂对其北部地区空气质量影响较大；6月份的主导风向为西南风和东北风，该热电厂对其东北部地区和西南部地区空气质量的影响较大。由于双高斯烟羽模型模拟的是烟气的扩散过程主要受风向、温度等气象条件的影响，因此模拟的其他污染物的扩散形态及影响范围与SO2一致。热电厂对2个监测站点SO2、NO2、PM2.5、可吸入颗粒物(PM10)的月均贡献值见表 2。
|
|
表 2 热电厂对2个监测站点污染物的月均贡献 |
S1和S2省控监测站点距离该热电厂分别为6.4和6.6 km。在二季度多发西南风的情况下，易处于下风向的位置，且距离正处在下风向6~10 km热电厂落地值较大的区域内[11]。通过模拟发现2019年5月24日04：00，该热电厂对上述2个站点各项污染物的贡献值最大，见表 3。这段时间的气象条件为风速1.3 m/s，风向为229°的西南风。
|
|
表 3 模拟期间热电厂对S1和S2站点的最大贡献值和贡献率 |
利用热电厂所在地区2019年4月1日—6月30日S1和S2站点SO2、NO2、PM10及PM2.5质量浓度与气象数据，分析了其各个风向时污染物均值，见表 4。由表 4可见，当风向为南西南、西南及西西南时，各项污染物均值均处于相对较高水平，这与上述热电厂在西南风时对所在地区各项污染物的贡献值较大一致。
|
|
表 4 S1和S2站点的各个风向下污染物均值μg/m3 |
(1) IFDM系统在双高斯烟羽模型的基础上，对线源和面源进行相应的处理，并可以耦合区域尺度模型，这样集成的模拟系统不仅可以模拟局地尺度所有污染源的扩散情况，同时还可以对局地的污染来源进行解析。
(2) 通过实际运用模型发现，该热电厂4—6月对S1、S2站点NO2的平均贡献值分别为2.0，3.2和4.4 μg/m3；对SO2的平均贡献值分别为1.9，3.1和3.5 μg/m3；对PM10的平均贡献值分别为0.45，0.7和0.85 μg/m3；对PM2.5的平均贡献值分别为0.25，0.45和0.55 μg/m3。
(3) 当风速为1.3 m/s, 风向为229°的西南风向时，该热电厂对S1和S2站点的污染物贡献值最大。
(4) 所在地区二季度的主导风向为西南风，因此建议在这段时间应注意该热电厂对所在地区空气质量的影响。在气象条件不利时，尤其是风向为西南风时，建议对该热电厂实行管控措施。
| [1] |
薛文博, 王金南, 杨金田, 等. 国内外空气质量模型研究进展[J]. 环境与可持续发展, 2013(3): 14-21. DOI:10.3969/j.issn.1673-288X.2013.03.004 |
| [2] |
薛文博, 许艳玲, 王金南, 等. 全国火电行业大气污染物排放对空气质量的影响[J]. 中国环境科学, 2016, 36(5): 1281-1288. DOI:10.3969/j.issn.1000-6923.2016.05.001 |
| [3] |
伯鑫, 田飞, 唐伟, 等. 重点煤电基地大气污染物扩散对京津冀的影响[J]. 中国环境科学, 2019(2): 514-522. DOI:10.3969/j.issn.1000-6923.2019.02.009 |
| [4] |
黄青, 程水源, 陈东升, 等. 近周边电厂源对北京市采暖期间SO2的贡献分析[J]. 环境科学研究, 2009, 22(5): 567-573. |
| [5] |
许艳玲, 蒋春来, 杨金田, 等. 我国火电行业NO2减排进展及控制对策研究[J]. 环境与可持续发展, 2014(5): 17-20. DOI:10.3969/j.issn.1673-288X.2014.05.004 |
| [6] |
王志轩, 潘荔, 赵鹏, 等. 现有燃煤电厂二氧化硫治理"十一五"规划研究[J]. 环境科学研究, 2007, 20(3): 142-147. DOI:10.3321/j.issn:1001-6929.2007.03.025 |
| [7] |
高兰兰, 傅成诚, 冯新斌, 等. 贵州东部某燃煤电厂汞排放对周边环境空气及土壤的影响[J]. 中国环境监测, 2018, 34(3): 51-58. |
| [8] |
薛志钢, 柴发合, 段宁, 等. 运用ISC3模型模拟电厂脱硫后的大气环境影响[J]. 环境科学研究, 2003, 16(5): 62-64. DOI:10.3321/j.issn:1001-6929.2003.05.017 |
| [9] |
黄青, 程水源, 陈东升, 等. 华北电厂源对北京及周边地区SO2影响[J]. 北京工业大学学报, 2009, 35(10): 1389-1395. DOI:10.11936/bjutxb2009101389 |
| [10] |
马英, 陶申鑫. 火电厂烟尘排放对重点城市的环境影响分析[J]. 中国环境监测, 1990, 6(5): 38-42. |
| [11] |
张秀义, 周爱国, 闫来洪, 等. 胜利电厂一期工程环境影响回顾评价——污染源和污染治理措施回顾评价[J]. 中国环境监测, 1996, 12(6): 30-34. |
| [12] |
王志勤, 颜文洪. 燃煤电力工程对城区环境空气质量影响分析[J]. 湖南理工学院学报, 2008, 21(2): 84-87. DOI:10.3969/j.issn.1672-5298.2008.02.027 |
| [13] |
张明华, 谢文勇, 唐晓青, 等. 河北省大气污染状况分析及对策研究[J]. 中国环境管理干部学院学报, 2015(1): 47-50. |
| [14] |
生态环境部. 2019年7月和1—7月全国空气质量状况[EB/OL].(2019-08-27)[2019-11-12].http://www.mee.gov.cn/.
|
| [15] |
韩立红. 河北省能源结构分析[J]. 河北经贸大学学报, 2010, 31(3): 64-68. |
| [16] |
BULTYNCK H, MALET L M, Van der Parren. Evaluation of atmospheric dilution factors for effluents diffused from an elevated continuous point source[J]. Tellus, 1972(5): 465-467. |
| [17] |
LEFEBVRE W, FIERENS F, TRIMPENEERS E, et al. Modeling the effects of a speed limit reduction on traffic-related elemental carbon (EC) concentrations and population exposure to EC[J]. Atmospheric Environment, 2011, 45: 197-207. DOI:10.1016/j.atmosenv.2010.09.026 |
| [18] |
VENKATRAM A, HORST T W. Approximating dispersion from a finite line source[J]. Atmospheric Environment, 2006, 40: 2401-2408. DOI:10.1016/j.atmosenv.2005.12.014 |
| [19] |
LEFEBVRE W, VERCAUTEREN J, SCHROOTEN L. Validation of the MIMOSA-AURORA-IFDM model chain for policy support: Modeling concentrations of elemental carbon in Flanders[J]. Atmospheric Environment, 2011, 45: 6705-6713. DOI:10.1016/j.atmosenv.2011.08.033 |
| [20] |
COSEMANS G, KRETZSCHMAR J G. The 1998 IFDM evaluation with the model validation kit[J]. International Journal of Environment and Pollution, 2000, 14: 122-130. DOI:10.1504/IJEP.2000.000533 |
| [21] |
LEFEBVRE W, VAN POPPEL M, MAIHEN B. Evaluation of the RIO-IFDM-street canyon model chain[J]. Atmospheric Environment, 2013, 77: 325-337. DOI:10.1016/j.atmosenv.2013.05.026 |
| [22] |
所在地区地方志编纂委员会. 所在地区志(1979—2009)[M]. 石家庄: 河北人民出版社, 2015: 75-78.
|