中國(guó)電力科學(xué)研究院有限公司 天津大學(xué) 成 嶺中國(guó)電力科學(xué)研究院有限公司 郭炳慶 鐘 鳴
關(guān)于農(nóng)村清潔供暖及其影響,眾多學(xué)者進(jìn)行了研究。周淑慧等人對(duì)北方農(nóng)村地區(qū)主要的清潔供暖方式進(jìn)行了對(duì)比分析并提出了建議[1];孫賦敏等人以山東農(nóng)村地區(qū)為例計(jì)算了清潔供暖的環(huán)境效益[2]。本文通過(guò)農(nóng)村清潔供暖實(shí)施情況,建立了供暖排放對(duì)大氣環(huán)境質(zhì)量影響的因子模型(A-H模型),分析清潔供暖對(duì)京津冀及周邊地區(qū)空氣質(zhì)量的影響程度,并綜合考慮人口、建筑面積、清潔能源等發(fā)展趨勢(shì),預(yù)測(cè)清潔供暖技術(shù)推廣對(duì)空氣質(zhì)量的改善。
本研究建立了A-H模型,以分析散煤替代對(duì)空氣質(zhì)量的改善程度,模型中考慮的因素包括氣象條件、城市特性累積污染物的二次生成和城市的自凈能力[3]。
因城市自凈能力而減少的細(xì)顆粒物(PM2.5)的濃度ρ1可通過(guò)下式計(jì)算:
ρ1=me-nt
(1)
式中m為城市最大凈化能力,μg/m3;n為凈化能力衰減系數(shù);t為時(shí)間,s。
供暖產(chǎn)生的PM2.5的濃度ρh為
ρh=AH(t)
(2)
式中A為直接排放因子;H(t)為燃煤消耗量,萬(wàn)t。
PM2.5的二次生成濃度ρ2為
ρ2=ξ(AH(t)-me-nt)
(3)
式中ξ為城市再生系數(shù)。
供暖燃煤減少與去除其他因素的PM2.5濃度之間的關(guān)系為
ρ′h=ρh-ρ1+ρ2
(4)
式中ρ′h為供暖燃煤減少而降低的PM2.5濃度,μg/m3。
排放空氣污染物的行業(yè),例如交通行業(yè)和工業(yè)等對(duì)PM2.5濃度有直接影響。但是,這類行業(yè)空氣污染物的排放基本不隨季節(jié)發(fā)生明顯變化[3]。因此,本研究以過(guò)渡季節(jié)的PM2.5濃度ρb為基準(zhǔn),并假定PM2.5濃度在供暖季的上升是由建筑供暖引起的。因此,去除城市特性干擾的PM2.5濃度根據(jù)下式計(jì)算:
ρh=ρs-ρb
(5)
式中ρs為供暖季實(shí)際PM2.5濃度,μg/m3。
與城市特性的影響不同,氣象因素變量彼此相關(guān)。初始數(shù)據(jù)包括風(fēng)向、風(fēng)速、壓力、降水和相對(duì)濕度。因此,本研究引入主成分分析方法,以減少氣象參數(shù)變量。以天津市為例,經(jīng)計(jì)算,5個(gè)氣象研究變量可以轉(zhuǎn)化為3個(gè)主成分因子,其相關(guān)的因子得分矩陣如表1所示。
表1 氣象因素的3個(gè)主成分及其系數(shù)
為了消除主成分對(duì)PM2.5濃度的影響,使用灰色關(guān)聯(lián)投影法對(duì)其進(jìn)行分解。每個(gè)主成分都可以分解為2個(gè)正交向量f和s。向量f與PM2.5濃度有關(guān),向量s與PM2.5濃度無(wú)關(guān),因此3個(gè)主成分定義式為
(6)
式中 cosα、cosβ、cosγ分別為R、P、F與ρh之間的相關(guān)系數(shù)。本研究引入綜合影響因子CIF表征3個(gè)主成分對(duì)PM2.5濃度的總影響,其表達(dá)式為
CIF=fR+fP+fF=Rcosα+Pcosβ+Fcosγ
(7)
cosα、cosβ和cosγ用皮爾遜相關(guān)系數(shù)進(jìn)行計(jì)算,通過(guò)計(jì)算,fR、fP、fF與ρh之間的皮爾遜相關(guān)系數(shù)分別為-0.262、-0.34和-0.2。因此CIF可以通過(guò)下式計(jì)算:
CIF=0.262R+0.34P+0.2F
(8)
為了消除氣象條件的影響,可計(jì)算ρh在氣象條件最差時(shí)的值ρ′h,即CIF最小的情況。
(9)
式中CIF,max和CIF,min分別為綜合影響因子的最大值和最小值;ρh,max和ρh,min分別為PM2.5濃度的最大值和最小值,μg/m3。
利用1stopt軟件通用全局優(yōu)化(UGO)計(jì)算直接排放因子A、城市最大凈化能力m、凈化能力衰減系數(shù)n和城市再生系數(shù)ξ。本文利用2006—2017年的逐月數(shù)據(jù)計(jì)算上述值,計(jì)算結(jié)果見(jiàn)表2;利用2018年的監(jiān)測(cè)數(shù)據(jù)驗(yàn)證模型。
表2 A、m、n和ξ計(jì)算結(jié)果
直接排放因子A的擬合公式為
A=1.445 42e-0.128 28t(R2=0.979)
(10)
城市最大凈化能力m的擬合公式為
m=7.880 8e0.001 38x(R2=0.93)
(11)
式中x為每年的供暖消耗煤量,萬(wàn)t。
凈化能力衰減系數(shù)n的擬合公式為
n=-0.000 7t2+0.004 1t+0.093 78
(R2=0.84)
(12)
城市再生系數(shù)ξ的擬合公式為
ξ=0.008 22e0.002 053x(R2=0.913)
(13)
利用2018年發(fā)布的PM2.5濃度值驗(yàn)證A-H模型的準(zhǔn)確性。實(shí)際值與模擬值的比較如圖1所示。引入平均絕對(duì)誤差EMA評(píng)估模型的準(zhǔn)確性,其定義為
圖1 模型計(jì)算結(jié)果的驗(yàn)證
(14)
式中yi為觀測(cè)值;y′i為計(jì)算值;p′為測(cè)試數(shù)據(jù)集中的數(shù)據(jù)數(shù)。
EMA為3.78 μg/m3,約為測(cè)試值的6%,說(shuō)明A-H模型的整體性能較好,可以較準(zhǔn)確評(píng)價(jià)供暖排放量對(duì)環(huán)境的影響。
圖2顯示了截至2018年京津冀大氣污染傳輸通道“2+26”城市清潔供暖改造情況。6個(gè)省/直轄市總計(jì)1 130萬(wàn)戶完成了清潔供暖改造,散煤消耗量減少了2 024萬(wàn)t。其中,天然氣、熱泵、直接電加熱和集中供暖改造占比分別為52.6%、14.1%、17.4%和15.9%。
圖2 截至2018年“2+26”城市的散煤替代改造情況
在2016年和2017年,分別有100萬(wàn)戶和530萬(wàn)戶的家庭從傳統(tǒng)的散煤燃燒供暖改造為清潔供暖,減少PM2.5排放超過(guò)20萬(wàn)t,從而顯著改善了京津冀及其周邊地區(qū)的空氣質(zhì)量。圖3顯示了2015年和2018年供暖季6個(gè)省/直轄市的PM2.5濃度分布。在大多數(shù)省份,PM2.5濃度得到明顯改善。
本研究利用犃A-H犎模型分析了到2018年由供暖燃煤消耗量減少引起PM2.5濃度下降的情況。如圖4所示,北京、石家莊、保定等城市的PM2.5濃度降低量較大,其余大部分城市的PM2.5濃度減少量也都約為20ug/m3。由供暖造成的空氣污染逐漸減少,清潔供暖對(duì)改善空氣環(huán)境具有積極的作用。
圖4 2018年“2+26”城市的PM2.5濃度減少量
對(duì)于經(jīng)濟(jì)發(fā)達(dá)、人口密集的城市群體,其空氣污染不再局限于一個(gè)城市,城市之間的空氣污染變化具有明顯的同步性和區(qū)域污染特征[4]。因此,為了研究“2+26”城市區(qū)域性的空氣污染相關(guān)性,基于地理相關(guān)性方法對(duì)城市群進(jìn)行了劃分,這些城市群的詳細(xì)信息及其2018年的散煤替代情況如表3所示。
表3 城市群的詳細(xì)信息
圖5給出了6個(gè)城市群2018年的平均散煤替代量、平均PM2.5濃度減少量和實(shí)際PM2.5濃度。其中城市群A和B的散煤替代量最多,平均替代量都超過(guò)了450萬(wàn)t,其他城市群的替代量平均為150萬(wàn)t左右。這是因?yàn)锳和B城市群中的大多數(shù)城市屬于京津冀地區(qū),這些地區(qū)的清潔供暖改造要求最為嚴(yán)格[5]。6個(gè)城市群的PM2.5濃度減少量與各自的散煤替代量密切相關(guān),城市群A和B的PM2.5濃度減少量最大。
從圖5可以看出,6個(gè)城市群體的實(shí)際PM2.5濃度差異很大。以散煤替代量相近的城市群A和B為例進(jìn)行分析,城市群A中北京和天津是最早進(jìn)行空氣污染治理的城市,并在2018年之前就取得了很大進(jìn)展,因此,城市群A在2018年供暖季的平均PM2.5濃度最小。城市群B包括石家莊、保定等城市,大規(guī)模的散煤替代已在一定程度上改善了空氣質(zhì)量,但受工業(yè)等因素影響,PM2.5濃度仍處于較高水平。
圖5 2018年6個(gè)城市群空氣質(zhì)量的改善情況
為了研究京津冀及其周邊地區(qū)清潔供暖政策對(duì)空氣質(zhì)量的深遠(yuǎn)影響,在考慮人口變化情況和可再生能源電力發(fā)展的情況下,本文利用A-H模型對(duì)空氣質(zhì)量進(jìn)行預(yù)測(cè)分析,并參考了清潔供暖政策設(shè)定的到2021年城市城區(qū)全部實(shí)現(xiàn)清潔供暖、縣城和城鄉(xiāng)結(jié)合部清潔供暖率達(dá)到80%以上、農(nóng)村地區(qū)清潔供暖率達(dá)到60%以上的目標(biāo)。在改造方式上,煤改電所需的電力終端燃煤量計(jì)入,而天然氣的燃燒被認(rèn)為只會(huì)對(duì)CO2的排放量產(chǎn)生影響,對(duì)于PM2.5的濃度影響不大。
人口變化趨勢(shì)會(huì)影響建筑供暖面積和供暖能源消耗,根據(jù)文獻(xiàn)[6],本文利用人口-發(fā)展-環(huán)境(PDE)模型預(yù)測(cè)未來(lái)人口,并假定城鄉(xiāng)人口比率(本文取1.33)不變。到2030年能源清潔化率達(dá)到40%[7],因此在本研究中,設(shè)定非化石能源的平均年增長(zhǎng)率為8.2%。
清潔供暖可以有效減少散煤的消耗,提高煤炭利用效率,減少污染物排放,但與之相對(duì)應(yīng)的用電需求也大大增加。根據(jù)計(jì)算,由于2018年的清潔供暖改造,“2+26”城市的平均電力需求增加了約20億kW·h。根據(jù)電力增長(zhǎng)情況分析,本文分別從燃煤電力(場(chǎng)景1)和清潔電力(場(chǎng)景2)2個(gè)場(chǎng)景進(jìn)行對(duì)比分析。
圖6顯示了2個(gè)場(chǎng)景PM2.5濃度與供暖消耗有關(guān)的趨勢(shì),以開(kāi)始實(shí)施清潔供暖改造的2016年作為基準(zhǔn)年,2016年供暖季北京、天津、河北及其周邊地區(qū)的PM2.5濃度分別為96.6、100.6、128.9和109.7 μg/m3。在場(chǎng)景1中,PM2.5濃度呈現(xiàn)先上升后下降的趨勢(shì)。這一趨勢(shì)與這些地區(qū)的人口變化趨勢(shì)相吻合,由于沒(méi)有考慮清潔電力的發(fā)展,PM2.5濃度將主要受供暖面積和人口的影響。相比2016年供暖季,預(yù)計(jì)到2030年,4個(gè)地區(qū)PM2.5濃度將分別下降48%、35%、29%和23%,主要原因是后期人口的下降及清潔供暖。在場(chǎng)景2中,考慮清潔電力的發(fā)展,4個(gè)地區(qū)的PM2.5濃度比場(chǎng)景1平均低3~4 μg/m3,而且這一差距在后期還會(huì)繼續(xù)增大;場(chǎng)景2的PM2.5濃度幾乎不存在場(chǎng)景1先上升的情況,也就是說(shuō)清潔電力的推廣基本能抵消人口增長(zhǎng)的壓力。到2030年,4個(gè)地區(qū)供暖季PM2.5濃度相比基準(zhǔn)年將分別下降55%、40%、32%和27%??梢?jiàn)場(chǎng)景2治理空氣污染的效果有明顯的提升。
本研究只對(duì)供暖相關(guān)的因素進(jìn)行了分析,因此所得的PM2.5濃度只考慮了供暖因素造成的污染排放降低,若2030年實(shí)現(xiàn)年平均PM2.5質(zhì)量濃度35 μg/m3的目標(biāo)[8],在考慮其他行業(yè)如交通運(yùn)輸業(yè)和重工業(yè)等的空氣污染防治計(jì)劃之后,實(shí)際的PM2.5濃度將會(huì)得到更加明顯的降低,可能完成這一目標(biāo)。
1) 通過(guò)推行散煤替代政策,京津冀及其周邊地區(qū)的散煤消耗已大大減少。到2018年,6個(gè)省/直轄市的1 130萬(wàn)戶家庭已經(jīng)實(shí)施了清潔供暖改造,通過(guò)用天然氣、電力和集中供暖取代散煤燃燒,使煤炭消耗減少了2 024萬(wàn)t,減少PM2.5排放超過(guò)20萬(wàn)t。
2) 清潔供暖改造可以明顯改善空氣質(zhì)量。通過(guò)A-H模型預(yù)測(cè),京津冀及其周邊地區(qū)的大多數(shù)城市的供暖季平均PM2.5質(zhì)量濃度減少量達(dá)到20 μg/m3。
3) 對(duì)PM2.5濃度的預(yù)測(cè)表明,考慮到人口的變化,充分發(fā)展利用清潔能源電力有利于更好地控制空氣污染,提升空氣質(zhì)量??紤]清潔能源電力比傳統(tǒng)燃煤電力PM2.5質(zhì)量濃度低3~4 μg/m3,北京、天津、河北及其周邊地區(qū)在采用傳統(tǒng)燃煤電力時(shí)PM2.5濃度將分別下降48%、35%、29%和23%。