基于氣象分型的大氣環(huán)境允許排放量測算方法

高 爽,李時蓓*,馬 巖,李 剛,伯 鑫,屈加豹,雷團(tuán)團(tuán),毛 娜,路瑞娟,任麗紅 (1.生態(tài)環(huán)境部環(huán)境工程評估中心,北京 100012；2.生態(tài)環(huán)境部國家環(huán)境保護(hù)環(huán)境影響評價數(shù)值模擬重點(diǎn)實(shí)驗(yàn)室,北京 100012；.中國環(huán)境科學(xué)研究院,北京 100012；.滄州市生態(tài)環(huán)境保護(hù)科學(xué)研究院,河北 滄州 061000)

以環(huán)境質(zhì)量改善為核心的大氣污染防治模式推動了大氣環(huán)境允許排放量研究的快速發(fā)展,近年來,探索基于PM2.5濃度達(dá)標(biāo)約束條件下的多污染物允許排放量計(jì)算方法成為研究熱點(diǎn).隨著對環(huán)境理論認(rèn)知水平的不斷提升,主要形成了A 值法、線性優(yōu)化法和模型模擬法等大氣環(huán)境允許排放量計(jì)算方法.但A 值法更適合計(jì)算理想狀態(tài)下的大氣容量,誤差很大[1-10];而線性優(yōu)化法雖能夠反映“排放-受體”的響應(yīng)關(guān)系,但該方法缺少邊界外區(qū)域的制約,也不能處理二次轉(zhuǎn)化污染問題[11-13];之后多位學(xué)者建立起來的模型模擬法可以綜合考慮地形、氣象、污染源排放、二次轉(zhuǎn)化等因素,反應(yīng)復(fù)雜的大氣物理化學(xué)過程,被廣泛應(yīng)用于允許排放量的計(jì)算研究中[14-21].例如李莉等[22]和盧燕宇等[23]從應(yīng)用的角度利用模型模擬了城市主要污染物的大氣環(huán)境容量;薛文博等[24]和李敏等[25]從算法的角度進(jìn)行優(yōu)化,基于CAMx 模型,通過優(yōu)化其內(nèi)部算法進(jìn)而得到主要污染物的環(huán)境容量.然而,模型模擬法技術(shù)復(fù)雜、計(jì)算量大、耗時較長,有一定局限性,不能對特定時段的大氣環(huán)境允許排放量進(jìn)行快速地估算.此外,目前允許排放量核算方法大多針對年際環(huán)境容量,缺少對于重污染天氣下的允許排放量核算研究.不利氣象條件是發(fā)生重污染過程的重要原因之一.研究不同氣象條件下的污染程度和特征比值變化對測算重污染天氣下的大氣環(huán)境容量、降低和消除重污染天氣有重要意義[26-30].

基于此,本文提出了一種快捷的大氣環(huán)境允許排放量測算方法,將基準(zhǔn)年氣象條件分型,研究分污染天氣類型的PM2.5組分濃度、二次轉(zhuǎn)化比例和遷散系數(shù),基于擴(kuò)散原理,結(jié)合污染物目標(biāo)控制濃度,快速測算特定氣象條件下達(dá)空氣質(zhì)量標(biāo)準(zhǔn)的大氣環(huán)境允許排放量.滄州市位于京津冀重點(diǎn)區(qū)域,屬于“2+26”城市之一,且PM2.5年均環(huán)境空氣質(zhì)量不達(dá)標(biāo), 秋冬季重污染頻發(fā),是典型的秋冬季高污染城市, 通過將本文提出的測算方法應(yīng)用于滄州市, 測算空氣質(zhì)量由優(yōu)～嚴(yán)重污染程度的7 種天氣類型的允許排放量,探索測算特定氣象條件下區(qū)域允許排放量的可行性,為制定城市重污染天氣應(yīng)急預(yù)案和分污染類型應(yīng)急響應(yīng)措施提供科學(xué)依據(jù).

1 材料與方法

本研究提出一種基于氣象條件分型的大氣污染物允許排放量測算方法,該方法基于排放清單和來源解析結(jié)果,通過分污染天氣類型的大氣擴(kuò)散規(guī)律和二次轉(zhuǎn)化規(guī)律研究,測算n 種天氣類型下的二次轉(zhuǎn)化比例和遷散系數(shù),并通過建立排放強(qiáng)度和空氣質(zhì)量濃度之間的關(guān)系快速測算大氣污染物允許排放量.算法過程如下:首先,綜合利用區(qū)域網(wǎng)格模型、受體模型,融合環(huán)境空氣質(zhì)量監(jiān)測數(shù)據(jù),快速提取各類源對PM2.5的貢獻(xiàn)濃度;第二步,通過統(tǒng)計(jì)研究區(qū)域基準(zhǔn)年的氣象條件,按污染程度劃分n 種天氣類型,結(jié)合基準(zhǔn)年污染源排放強(qiáng)度,計(jì)算n 種天氣類型m 類污染源的遷移擴(kuò)散系數(shù)F(m,n);第三步,基于區(qū)域網(wǎng)格模型模擬結(jié)果、多年觀測的6 因子質(zhì)量監(jiān)測數(shù)據(jù)和PM2.5組分?jǐn)?shù)據(jù),獲得n種天氣類型達(dá)PM2.5規(guī)劃目標(biāo)下的PM10、SO42-和NO3-的控制濃度;進(jìn)一步利用遷散系數(shù)計(jì)算第j 種天氣類型下達(dá)控制濃度目標(biāo)的允許排放量.運(yùn)用該創(chuàng)新方法,可以計(jì)算特定氣象條件下PM2.5的允許排放量,尤其是重污染條件下的允許排放量,可為制定重污染應(yīng)急管控措施提供依據(jù).技術(shù)路線如圖1 所示.

圖1 以濃度為控制目標(biāo)的大氣環(huán)境允許排放量核算方法技術(shù)路線Fig.1 Technical route of atmospheric environmental capacity accounting method aimed at concentration control target

1.1 多尺度模型融合監(jiān)測數(shù)據(jù)的高分辨解析方法

1.1.1 區(qū)域網(wǎng)格模型模擬 采用WRF-CAMx模式雙層嵌套的方法,外層模擬區(qū)域?yàn)?7°～161°E, 1°～59°N,網(wǎng)格水平分辨率為 36km×36km,共劃分為200×160個網(wǎng)格;內(nèi)層模擬網(wǎng)格水平分辨率為12km×12km,共劃分為119×101 個網(wǎng)格,涵蓋京津冀及周邊城市.垂直方向分為20 層.污染源清單選取課題組建立的滄州市2018 年污染源清單.通過區(qū)域網(wǎng)格模型模擬,得到本地源貢獻(xiàn)P區(qū)域模型本地.

1.1.2 CMB(Chemical Mass Balance,化學(xué)質(zhì)量平衡)模型來源解析 本研究選擇滄州市監(jiān)測站、青縣海子公園、黃驊財(cái)政局培訓(xùn)中心和泊頭市政府4 個環(huán)境受體點(diǎn)開展PM2.5樣品采集,采樣時間為2017 年10 月15 日～2018 年1 月31 日及2018 年10 月15日～2019 年1 月31 日,共采集有效樣品1948 個,對各點(diǎn)位PM2.5樣品中的水溶性離子、無機(jī)元素和碳組分進(jìn)行化學(xué)分析和討論.基于顆粒物采樣和環(huán)境受體PM2.5的化學(xué)分析,綜合利用CMB 受體模型解析得到滄州市工業(yè)源、揚(yáng)塵源、交通源和民用燃燒源的分擔(dān)率.

1.1.3 融合解析方法 首先利用區(qū)域網(wǎng)格模型獲得本地貢獻(xiàn)占比P區(qū)域模型本地,然后基于受體模型獲得第i 類源貢獻(xiàn)占比P受體解析i占比,再融合環(huán)境空氣質(zhì)量監(jiān)測數(shù)據(jù),提取第i 類源在第j 種天氣類型下的貢獻(xiàn)濃度,見公式(1).

式中: C(i,j)為第j 種天氣類型下第i 類源對PM2.5的貢獻(xiàn)濃度,μg/m3;Cj監(jiān)測為第j 種天氣類型下PM2.5環(huán)境質(zhì)量監(jiān)測濃度,μg/m3;P區(qū)域模型本地為本地源對PM2.5濃度的貢獻(xiàn)占比,%;P受體解析i占比為受體模型解析第i類源對PM2.5濃度的貢獻(xiàn)占比,%.

1.2 基于氣象條件分型的污染特征和轉(zhuǎn)化規(guī)律研究方法

1.2.1 氣象條件分型方法 統(tǒng)計(jì)研究區(qū)域基準(zhǔn)年的風(fēng)速、風(fēng)向、穩(wěn)定度、相對濕度等氣象要素,根據(jù)PM2.5濃度從優(yōu)～嚴(yán)重污染程度,確定對應(yīng)的氣象要素,劃分n種天氣類型(剔除雨雪沙塵天).氣象要素統(tǒng)計(jì)原則以該種天氣類型下風(fēng)向頻次超過80%的風(fēng)向?yàn)橹鲗?dǎo)風(fēng)向、以該種天氣類型下2min 風(fēng)速平均值為平均風(fēng)速、以該種天氣類型下相對濕度的最小值和最大值為相對濕度范圍.

1.2.2 分天氣類型的污染特征和轉(zhuǎn)化規(guī)律研究方法 對于每種天氣類型,計(jì)算該類型下的天數(shù)(或小時數(shù))占全部統(tǒng)計(jì)時間段的比例、常規(guī)污染物平均濃度、水溶性離子濃度及其與前體物比值(二次轉(zhuǎn)化系數(shù))等,據(jù)此研究分污染天氣類型的二次轉(zhuǎn)化規(guī)律.

1.3 不同天氣類型下的遷散系數(shù)確定方法

根據(jù)擴(kuò)散原理,污染濃度不僅與排放源強(qiáng)成正比關(guān)系,同時還與污染源排放布局和排放方式、氣象條件、污染物的遷移轉(zhuǎn)化以及地形地表等條件密切相關(guān),在特定區(qū)域排放布局和排放方式基本不變的假設(shè)前提下,污染濃度的大小取決于排放強(qiáng)度、氣象條件和遷移、轉(zhuǎn)化等因素[31-32].

對于PM2.5,基于融合解析方法和氣象分型研究,按式(2)計(jì)算得到第j 種天氣類型第i 類源的遷散系數(shù).

式中: F(i,j)為第i 類污染源在第j 種天氣類型的遷散系數(shù),(μg?d)/(m3?t); C(i,j)基準(zhǔn)為基準(zhǔn)年第i 類污染源在第j 種天氣類型下的質(zhì)量濃度,μg/m3; Q(i,j)基準(zhǔn)為基準(zhǔn)年第i 類污染源在第j 種天氣類型下的排放強(qiáng)度,t/d.

1.4 特定天氣類型下的允許排放量測算方法

1.4.1 多種污染物控制濃度確定方法 基于區(qū)域網(wǎng)格模型模擬的P區(qū)域模型本地,結(jié)合6 因子質(zhì)量監(jiān)測數(shù)據(jù)和PM2.5組分?jǐn)?shù)據(jù),分析n 種天氣類型下二次污染物占比和二次轉(zhuǎn)化系數(shù).基于PM2.5規(guī)劃目標(biāo),綜合上述分析方法,獲得n 種天氣類型下SO42-和NO3-的控制濃度.

1.5 數(shù)據(jù)來源

本研究區(qū)域網(wǎng)格模型和CMB 受體模型解析結(jié)果引用《大氣重污染成因與治理攻關(guān)項(xiàng)目——滄州市大氣污染防治綜合解決方案研究》(DQGG-05-滄州07)項(xiàng)目成果,清單數(shù)據(jù)采用課題組編制的滄州市2018 年源清單,滄州市2018 年全年環(huán)境空氣質(zhì)量監(jiān)測數(shù)據(jù)來自大氣環(huán)境科學(xué)綜合數(shù)據(jù)采集與共享平臺,氣象條件數(shù)據(jù)由滄州市氣象局提供.

2 滄州市應(yīng)用實(shí)例與結(jié)果討論

2.1 融合解析方法確定滄州市主要源類貢獻(xiàn)濃度

表1 滄州市主要污染源類對PM2.5 貢獻(xiàn)濃度Table 1 Contributions of main sources to PM2.5

根據(jù)1.1 節(jié)方法,采用區(qū)域網(wǎng)格模型模擬得到滄州市本地源貢獻(xiàn)75%;根據(jù)CMB 受體模式解析結(jié)果,滄州市工業(yè)源貢獻(xiàn)64%,揚(yáng)塵源貢獻(xiàn)20%,交通源和民用燃燒源均貢獻(xiàn)8%.2018 年滄州市PM2.5年均濃度為59μg/m3,運(yùn)用公式(1),計(jì)算得到工業(yè)源對PM2.5貢獻(xiàn)濃度為28.32μg/m3,同理得到揚(yáng)塵源、交通源和民用燃燒源貢獻(xiàn)濃度,結(jié)果見表1.

2.2 氣象條件分型結(jié)果

基于滄州市近2a 秋冬季氣象條件數(shù)據(jù),根據(jù)1.2.1 節(jié)所述方法,按優(yōu)、良、輕度污染、中度污染、重度污染和嚴(yán)重污染分別統(tǒng)計(jì)風(fēng)向、風(fēng)速、相對濕度、靜小風(fēng)持續(xù)時間4 項(xiàng)氣象要素,由于污染水平為良時可以歸納為2 種不同的氣象條件,因此最終劃分為7 種天氣類型(已剔除雨雪沙塵天).氣象條件分型結(jié)果見表2.

由表2 可知,天氣類型1 以東北風(fēng)主導(dǎo),平均風(fēng)速最大,相對濕度最低,擴(kuò)散條件最有利;天氣類型7靜小風(fēng)時間大于72h,相對濕度在80%～100%,擴(kuò)散條件最不利.從天氣類型1～7,靜小風(fēng)持續(xù)時間不斷升高,平均風(fēng)速越來越小、相對濕度逐漸增大,污染程度呈加重趨勢.污染發(fā)生時段多以偏南風(fēng)和本地靜小風(fēng)主導(dǎo).

表2 滄州市氣象條件分型結(jié)果Table 2 7 typical meteorological conditions of the city of Cangzhou

2.3 基于不同氣象條件分型的遷移擴(kuò)散和轉(zhuǎn)化規(guī)律研究

2.3.1 不同天氣類型下的污染規(guī)律研究 由表3可知,從天氣類型1～7 污染呈加重趨勢,隨著污染的加重,PM2.5/PM10比值不斷增大,細(xì)顆粒物占比增高說明二次轉(zhuǎn)化加劇;硝酸鹽/NO2和硫酸鹽/SO2的濃度比總體呈上升趨勢,硫酸鹽/SO2平均比值要高于硝酸鹽/NO2,說明本地高濕度環(huán)境下,更有利于SO2的二次轉(zhuǎn)化,這與相關(guān)研究結(jié)論有一致性[33].

表3 滄州市7 種天氣類型下的污染特征Table 3 Pollution characteristics under 7 typical meteorological conditions in the city of Cangzhou

表4 滄州市7 種天氣類型下主要污染物的遷散系數(shù)F[(μg?d)/(m3?t)]Table 4 Migration and diffusion coefficients under 7typical meteorological conditions

2.3.2 不同天氣類型下遷散系數(shù)確定 根據(jù)排放清單,滄州市2018 年一次PM2.5排放量為156.71t/d、PM10為328.22t/d、SO2為67.40t/d、NOx為232.05t/d.按公式(2)計(jì)算得到滄州市7種天氣類型下四種污染物的遷散系數(shù)F,如表4 所示.可以看到,從天氣類型1～7,4 種污染物的遷散系數(shù)總體不斷增大,表明隨著污染程度的加重,單位排放量的污染物對環(huán)境空氣質(zhì)量濃度的貢獻(xiàn)逐漸增大.

2.4 滄州市不同天氣類型下允許排放量測算

2.4.1 PM2.5、PM10、NO2與SO2的控制濃度確定《環(huán)境質(zhì)量標(biāo)準(zhǔn)》規(guī)定的 PM2.5達(dá)標(biāo)濃度為35μg/m3,PM2.5主要由一次排放貢獻(xiàn)和二次轉(zhuǎn)化構(gòu)成,由滄州市秋冬季PM2.5組分?jǐn)?shù)據(jù)及氣象條件分型結(jié)果,獲得7 種天氣類型下NO3-和SO42-的占比,天氣類型1～7,NO3-和SO42-濃度的占比在39.3%～67.3%之間.根據(jù)式(3)、(4),獲得7 種天氣類型下NO2、SO2和PM10控制濃度見表5.

表5 滄州市主要污染物目標(biāo)控制濃度(μg/m3)Table 5 Target concentrations of main pollutants(μg/m3)

2.4.2 天氣類型1～7 達(dá)控制濃度目標(biāo)下的允許排放量測算 由表6可知,從天氣類型1～7,各項(xiàng)主要污染物的允許排放量遞減.天氣類型1 下,平均風(fēng)速約為2.5m/s,平均相對濕度50%以下,擴(kuò)散條件最好,滄州市一次PM2.5、PM10、NOx和SO42-的允許排放量分別為273.36, 486.25, 237.05 和73.27t/d;天氣類型7下,靜小風(fēng)持續(xù)時間分別在72h 以上,平均相對濕度為80%～100%,擴(kuò)散條件最不利,滄州市一次PM2.5、PM10、NOx和SO2的允許排放量急劇下降,分別為27.23, 61.82, 57.89 和17.16t/d.

表6 滄州市不同氣象條件下4 種污染物的減排比例Table 6 Reduction percentages of main pollutants under 7 typical meteorological conditions

從減排比例結(jié)果來看,要達(dá)到PM2.5目標(biāo)控制濃度,從天氣類型1～7,各項(xiàng)主要污染物的減排比例遞增.天氣類型1 下各污染物已達(dá)到控制目標(biāo),不需要再做削減;天氣類型7 可認(rèn)為發(fā)生嚴(yán)重污染過程,氣象條件最為不利,一次PM2.5、PM10、NOx和SO2減排比例需分別在82.62%、81.17%、75.05%和74.54%以上,才可達(dá)到環(huán)境質(zhì)量控制目標(biāo).可以看出最不利氣象條件下大氣環(huán)境允許排放量很小,允許排放量約是有利氣象條件的1/6～1/10.因此,在該氣象條件下,需要更大力度削減污染物排放,才能避免發(fā)生重污染過程.

通過該方法獲得了滄州市基于氣象條件分型的大氣環(huán)境允許排放量.該方法是以源清單、源解析研究結(jié)果為基礎(chǔ)的,由于京津冀、汾渭平原、長三角等重點(diǎn)區(qū)域大部分城市均已開展源清單和源解析研究,因此該方法對這些城市具有可行性,可以快速測算特定氣象條件下的允許排放量.

2.5 研究方法的不確定性分析

由于污染濃度與排放源強(qiáng)、污染源排放布局、排放方式、氣象條件、污染物的遷移沉降、地形地表和物理化學(xué)反應(yīng)等條件存在著密切的關(guān)系,大氣污染物允許排放量的計(jì)算十分復(fù)雜,本研究旨在排放布局和排放方式不變的前提下,提出一種簡便的快速測算方法.本研究的不確定性主要體現(xiàn)在以下幾方面.

2.5.1 基準(zhǔn)年選取的不確定性 本研究是統(tǒng)計(jì)區(qū)域內(nèi)基準(zhǔn)年的氣象條件、環(huán)境質(zhì)量監(jiān)測數(shù)據(jù)、顆粒物組分?jǐn)?shù)據(jù)、模型模擬和解析結(jié)果等,但是氣象條件分型研究可能受到基準(zhǔn)年選取的影響,同時,污染源排放布局和排放方式、PM2.5來源解析結(jié)果也與基準(zhǔn)年選取有關(guān),因此需要考慮基準(zhǔn)年選取對測算結(jié)果帶來的不確定性影響.

2.5.2 數(shù)據(jù)樣本代表性影響 本研究案例中PM2.5組分?jǐn)?shù)據(jù)是采用近2a 秋冬季的膜采樣分析數(shù)據(jù),采樣時間為秋冬季4 個月(當(dāng)年10 月份～次年1 月份),對于特定天氣類型下(例如嚴(yán)重污染日)的樣本數(shù)有限,后續(xù)可進(jìn)一步增加采樣分析數(shù)量,提高數(shù)據(jù)代表性.

2.5.3 模型計(jì)算的不確定性 本研究方法基于區(qū)域網(wǎng)格模型和受體模型解析結(jié)果.由于模型模擬本身在參數(shù)選取、模型計(jì)算等方面具有一定不確定性,模型計(jì)算結(jié)果疊加也會對測算結(jié)果帶來影響.

3 結(jié)論

3.1 本研究提出一種基于氣象條件分型的大氣污染物允許排放量測算方法,該方法能快速測算特定污染天氣類型下的允許排放量,對以改善環(huán)境空氣質(zhì)量為目標(biāo)核算最大允許排放量探索了一條可行的技術(shù)路線.

3.2 本研究以滄州市為例,得到滄州市秋冬季7 種天氣類型,天氣類型1 的氣象條件最有利于擴(kuò)散,天氣類型7最為不利,從天氣類型1～7,污染程度不斷加重,靜小風(fēng)持續(xù)時間不斷升高,平均風(fēng)速不斷減小、相對濕度增大,硝酸鹽/NO2和硫酸鹽/SO2的濃度比總體上升.

3.3 以滄州市PM2.5控制目標(biāo)35μg/m3為約束,在最不利氣象條件下,滄州市一次PM2.5、PM10、NOx和SO42-排放量削減率應(yīng)分別在82.62%、81.17%、75.05%和74.54%以上,才能達(dá)到空氣質(zhì)量標(biāo)準(zhǔn).不利氣象條件下大氣環(huán)境允許排放量很小,需要更大力度減少污染物排放,才能避免發(fā)生重污染天氣.