張麗，李磊，蘭紅平

（1.深圳市國(guó)家氣候觀(guān)象臺(tái)，廣東深圳 518040；2.深圳南方強(qiáng)天氣研究重點(diǎn)實(shí)驗(yàn)室，廣東深圳 518040；3.中山大學(xué)大氣科學(xué)學(xué)院，廣東珠海 519082；4.深圳市氣象局，廣東深圳 518040）

自20世紀(jì)60年代開(kāi)始，氣象條件對(duì)空氣污染作用的問(wèn)題就一直備受關(guān)注。污染物濃度的變化受到多種要素的影響，如污染排放源、治理措施、氣象條件等，這些要素錯(cuò)綜復(fù)雜，在不同尺度的時(shí)間周期所起的作用也不相同。目前大多數(shù)研究氣象條件對(duì)污染的影響都是基于同一種假設(shè)，在短期內(nèi)（月尺度以下）不考慮污染源的變化，氣象條件是重要的影響因子。很多分析表明，引起污染質(zhì)量濃度日均值波動(dòng)因素的70%是氣象條件變化；污染物質(zhì)量濃度與風(fēng)速、混合層高度和降水呈負(fù)相關(guān)［1-4］。通過(guò)分析空氣污染過(guò)程對(duì)應(yīng)的天氣形勢(shì)，歸納總結(jié)出容易導(dǎo)致空氣污染的局地天氣類(lèi)型［5-7］。例如，鄒旭東等［8］提出我國(guó)北方地區(qū)大氣中污染物濃度過(guò)高多出現(xiàn)在地面高壓系統(tǒng)的控制下；廖碧婷等［9］計(jì)算了垂直交換系數(shù)，對(duì)空氣污染物的垂直輸送能力進(jìn)行評(píng)估，并嘗試了對(duì)灰霾天氣和無(wú)視程障礙天氣進(jìn)行預(yù)報(bào)；張恒德等［10］挑選發(fā)生大氣污染時(shí)的氣象要素及其閾值條件，通過(guò)權(quán)重求和初步構(gòu)建了靜穩(wěn)天氣綜合指數(shù)（SWI），但將污染狀況放到年以上的更長(zhǎng)的時(shí)間考慮時(shí)，必須考慮污染源排放的變化情況；張麗等［11］對(duì)深圳市近35年灰霾日變化的成因分析中提出，氣象條件不是深圳灰霾日年際變化的主要原因，產(chǎn)業(yè)政策引導(dǎo)下產(chǎn)業(yè)結(jié)構(gòu)升級(jí)才是空氣質(zhì)量轉(zhuǎn)變的根本原因。因此在評(píng)估氣象條件對(duì)年以上尺度的大氣污的貢獻(xiàn)時(shí)，需要剔除因大氣治理及排放源變化等非氣象要素對(duì)污染的影響。現(xiàn)實(shí)中二者是相互交織相互作用，難以定量分離的。本研究通過(guò)回歸去傾斜法計(jì)算了一個(gè)污染氣象年景指數(shù)（PMI），并進(jìn)行了相關(guān)性驗(yàn)證，可以定量分析氣象條件對(duì)污染物年均值質(zhì)量濃度變化的影響程度，將污染氣象條件給出具有物理意義的、量化的表達(dá)，可定量評(píng)估氣象條件對(duì)大氣污染的影響，同時(shí)也便于開(kāi)展污染氣象條件的年度預(yù)報(bào)。

1 污染氣象年景指數(shù)（PMI）構(gòu)建方法

氣象數(shù)據(jù)選用深圳市國(guó)家氣候觀(guān)象臺(tái)2011—2018年全市海拔高度70 m以上的自動(dòng)氣象站點(diǎn)風(fēng)速、降雨數(shù)據(jù)。污染物數(shù)據(jù)選用深圳市環(huán)境監(jiān)測(cè)站2011—2018年平均PM2.5質(zhì)量濃度值。

污染氣象年景指數(shù)（PMI）定義為氣象條件變化所導(dǎo)致的細(xì)顆粒物PM2.5年值變化的指標(biāo)。

其中，x1、x2、xn是敏感氣象因子；α、β、λ、μ是因子權(quán)重；ρ（PM2.5氣象）是指受氣象條件影響的部分質(zhì)量濃度（μg/m3）；ρ（PM2.5非氣象）是指不受氣象條件影響的剩余部分質(zhì)量濃度（μg/m3）。

回歸去傾法簡(jiǎn)介：

（1）將2011—2018年 ρ（PM2.5）年值進(jìn)行回歸分析，得到最優(yōu)回歸方程并計(jì)算逐年回歸值即ρ（PM2.5非氣象）；

（2）將逐年 ρ（PM2.5）年值減去回歸值ρ（PM2.5非氣象），得到，負(fù)值代表氣象條件有利于PM2.5的減少，正值代表氣象條件有利于PM2.5的增加；

（3）將 PM2.5氣象進(jìn)行量綱為一化，得到標(biāo)準(zhǔn)化ρ（PM2.5氣象）（圖1）。無(wú)量綱一化方法：標(biāo)準(zhǔn)化ρ（PM2.5氣象）=ρ（PM2.5氣象）/ρ（PM2.5氣象max）×100。

2 回歸去傾法計(jì)算ρ（PM2.5氣象）

2011年以來(lái)全市ρ（PM2.5）年值受排放源、減排措施和氣象因素的綜合影響，呈現(xiàn)明顯的波動(dòng)下降趨勢(shì)，最優(yōu)回歸方程y=-2.440 5x+44.107（R2=0.91）。經(jīng)過(guò)回歸去傾法計(jì)算2011—2018年 ρ（PM2.5氣象）占 ρ（PM2.5）年均值-9% ～8%（圖 1），說(shuō)明氣象條件對(duì) ρ（PM2.5）年值的影響最大接近10%。

圖1 ρ（PM2.5）年值及 ρ（PM2.5氣象）年變化

3 ρ（PM2.5氣象）與氣象因子的相關(guān)性驗(yàn)證

3.1 無(wú)降雨弱風(fēng)日數(shù)與ρ（PM2.5氣象）的相關(guān)性分析

無(wú)明顯降雨弱風(fēng)日是指當(dāng)日降雨量5 mm以下，且24 h 10 min平均風(fēng)速≤2 m/s。通過(guò)對(duì)全市70 m高度以上27個(gè)站點(diǎn)的相關(guān)性篩選，選出南山站（海拔104.4 m）和大浪站（海拔81.5 m）的無(wú)降雨弱風(fēng)日數(shù)與 ρ（PM2.5氣象）相關(guān)性最高，相關(guān)系數(shù)分別為0.64和0.60（圖2），并通過(guò)顯著性檢驗(yàn)，因此選取無(wú)降雨弱風(fēng)日作為ρ（PM2.5氣象）影響因子。

圖2 全市海拔高度70 m以上站點(diǎn)無(wú)降雨弱風(fēng)日數(shù)與ρ（PM2.5氣象）相關(guān)系數(shù)分布

3.2 降雨日數(shù)與ρ（PM2.5氣象）的相關(guān)性分析

降雨日數(shù)分日雨量5、10及15 mm 3種情況進(jìn)行對(duì)比分析，結(jié)果表明，5 mm以上降雨日數(shù)與全市ρ（PM2.5氣象）的相關(guān)性最好，最高相關(guān)系數(shù)達(dá)-0.75以上，并通過(guò)顯著性檢驗(yàn)，表明5 mm以上降雨日數(shù)越多，對(duì)ρ（PM2.5）年值的降低越有利。從站點(diǎn)選擇來(lái)看，福田、羅湖、坪山站的相關(guān)性最好，相關(guān)系數(shù)＞-0.7，大浪站（-0.62）、南山站（-0.64）的相關(guān)性次之。結(jié)合弱風(fēng)的相關(guān)性綜合來(lái)看，選擇南山冬半年5 mm以上降雨日數(shù)作為全市 ρ（PM2.5氣象）的影響因子（圖3）。

圖 3 2011—2018年雨量達(dá) 5（a）、10（b）、15 mm（c）的年降雨日數(shù)與 ρ（PM2.5氣象）的相關(guān)系數(shù)

通過(guò)相關(guān)性分析，選取了無(wú)降雨弱風(fēng)日數(shù)（d0）和 5 mm以上降雨日數(shù) （d5）作為ρ（PM2.5氣象）變化因子。無(wú)降雨弱風(fēng)日數(shù)越多對(duì)污染物擴(kuò)散越不利，為正相關(guān)；5 mm以上降雨日數(shù)越多，污染物清除作用越明顯，為負(fù)相關(guān)。

4 污染氣象年景指數(shù)（PMI）構(gòu)建

南山站無(wú)降雨弱風(fēng)日數(shù)與5 mm以上降雨日數(shù)2個(gè)氣象因子各占比重分別為80%和20%時(shí)，綜合相關(guān)性最高，相關(guān)系數(shù)可達(dá)0.75并通過(guò)顯著性檢驗(yàn)。因此建立污染氣象年景指數(shù)（PMI）：

其中，d0為無(wú)降雨弱風(fēng)日數(shù)的量綱為一的值；d5為5 mm以上降雨日數(shù)的量綱為一的值

2011—2018年污染氣象年景指數(shù)變化表明（圖4），2012、2016和 2017年污染氣象年景較好。2013、2015、2018年污染氣象年景較差。

圖4 2011—2018年污染氣象年景指數(shù)年變化

5 業(yè)務(wù)運(yùn)用

2019年初受厄爾尼諾影響，赤道中東太平洋海溫偏高，副高持續(xù)偏強(qiáng)偏西，同時(shí)西南季風(fēng)偏強(qiáng)，整個(gè)華南水汽輸送偏多，廣東省入汛偏早，1—3月大氣擴(kuò)散條件較有利，預(yù)計(jì)厄爾尼諾事件將持續(xù)到夏季。近10年深圳開(kāi)汛期偏早相似年如2009年（3月6日）、2013年（3月26日）和2014年（3月30日）。厄爾尼諾相似年有2005、2010、2016年。從這些年份中選取1—3月的500 hPa環(huán)流形勢(shì)、水汽輸送和2019年最為相似的年份2010年（圖5）。

圖5 2019年（a）和2010年（b）1—3月500 hPa高度場(chǎng)（等值線(xiàn)，單位：gpm）和距平（陰影，單位：gpm）

綜上所述，入汛偏早的3年環(huán)境氣象綜合指數(shù)分別為62、92和79（均值78）；大氣環(huán)流最相似的2010年環(huán)境氣象綜合指數(shù)70。綜合給出2019年污染氣象條件指數(shù)70～78，好于2018年（93）。可使 2019年 ρ（PM2.5）年值在 ρ（PM2.5非氣象）上增加3% ～5%，約23～24μg/m3，同比仍減少10%。其中夏半年降雨偏多，大氣環(huán)境同比較好；冬半年降雨偏少，無(wú)降雨弱風(fēng)日數(shù)偏多的概率大，大氣環(huán)境較去年同期略差（圖6）。

圖6 2019年污染氣象年景指數(shù)年變化及預(yù)測(cè)

6 結(jié)論

利用2011—2018年深圳市PM2.5質(zhì)量濃度及全市60個(gè)氣象站點(diǎn)氣象數(shù)據(jù)，開(kāi)展氣象條件對(duì)ρ（PM2.5）年值的影響分析，建立了深圳市污染氣象年景指數(shù)，并開(kāi)展了業(yè)務(wù)運(yùn)用。

1）基于回歸去傾法構(gòu)建了污染氣象年景指數(shù)，該指數(shù)計(jì)算結(jié)果表明近8年來(lái)環(huán)境氣象條件對(duì)全市平均PM2.5質(zhì)量濃度的影響接近10%，即在排放源不變的情況下，污染氣象年景指數(shù)最低的年份，大氣擴(kuò)散條件最好，年平均PM2.5質(zhì)量濃度可下降10%。

2）2017—2018年污染氣象條件指數(shù)由71增加至93，大氣環(huán)境氣象條件轉(zhuǎn)差，氣象條件使得全市平均PM2.5在非氣象值得基礎(chǔ)上增了5.4%。

3）預(yù)計(jì)在厄爾尼諾影響和入汛偏早的情況下，2019年大氣環(huán)境綜合指數(shù)70～78，大氣擴(kuò)散條件介于2017和2018年之間，氣象條件可使全市平均 PM2.5質(zhì)量濃度在 ρ（PM2.5非氣象）基礎(chǔ)上增加3%～5%，約23～24μg/m3。