魏 玉， 胡二琴， 穆新宇

(湖北工業(yè)大學(xué)理學(xué)院， 湖北 武漢 430068)

大氣中高濃度的NO2與氨、水分等在日照等條件下會(huì)產(chǎn)生化學(xué)反應(yīng)生成臭氧等二次污染物，不僅嚴(yán)重影響空氣質(zhì)量，而且危害環(huán)境效益[1-2]，是我國(guó)及歐美國(guó)家環(huán)保部門監(jiān)測(cè)的主要大氣污染物。研究發(fā)現(xiàn)，長(zhǎng)時(shí)間暴露在富含高濃度NO2或O3大氣環(huán)境中，極易導(dǎo)致人群尤其是兒童和老人產(chǎn)生肺部及呼吸系統(tǒng)疾病[3-4]。隨著經(jīng)濟(jì)的快速發(fā)展，我國(guó)已經(jīng)成為全球氮氧化物污染最為嚴(yán)重的地區(qū)之一[5-6]。

2017年國(guó)家環(huán)保部城市質(zhì)量報(bào)告顯示，74個(gè)實(shí)施新標(biāo)準(zhǔn)第一階段監(jiān)測(cè)的城市中，湖北省武漢市空氣質(zhì)量綜合指數(shù)排名48，在74個(gè)城市中處于下游。湖北省環(huán)境質(zhì)量公告顯示，按照《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB 3095—2012)評(píng)價(jià)，除神農(nóng)架以外，2016年至2018年全省其余城市空氣質(zhì)量均未達(dá)到二級(jí)標(biāo)準(zhǔn)，其中2017年NO2濃度較2016年上升7.7%，O3濃度與2016年持平；2018年，NO2濃度與2017年持平，O3濃度值較2017年上升10.8%。由此可以看出，湖北省治理大氣污染問(wèn)題已刻不容緩。

不少學(xué)者已對(duì)我國(guó)空氣質(zhì)量進(jìn)行了探究[5,7-14]，但大多數(shù)是基于傳統(tǒng)的統(tǒng)計(jì)方法，少有從函數(shù)型數(shù)據(jù)尤其是多元函數(shù)型主成分的角度對(duì)湖北省的空氣質(zhì)量進(jìn)行分析。由于空氣質(zhì)量及其相關(guān)指標(biāo)具有函數(shù)特征，所以本文從函數(shù)型高維數(shù)據(jù)的角度，采用函數(shù)型二元主成分挖掘出相互影響的兩種污染氣體之間的聯(lián)合變異方式，最大程度地利用變量之間的關(guān)系以及數(shù)據(jù)提供的信息，這樣，根據(jù)分析結(jié)果去制定相應(yīng)治理方案會(huì)更加準(zhǔn)確有效。

1 數(shù)據(jù)處理

1.1 數(shù)據(jù)情況

原始數(shù)據(jù)采自中國(guó)環(huán)境監(jiān)測(cè)總站。選取湖北省12個(gè)省轄市和1個(gè)自治州共13個(gè)行政區(qū)的大氣NO2與O3日均濃度數(shù)據(jù)，時(shí)間跨度為2018年10月21日至2019年10月20日，共365天。數(shù)據(jù)狀態(tài)良好，每個(gè)城市365天的數(shù)據(jù)缺失值均少于8個(gè)，采用該缺失值的相鄰前后共4天數(shù)據(jù)的平均值進(jìn)行填補(bǔ)。

1.2 曲線擬合

(a)NO2

(b) O3圖 1 湖北省13個(gè)城市大氣中NO2和O3變化趨勢(shì)

湖北省13個(gè)城市大氣中NO2和O3變化趨勢(shì)如圖1所示。采用傅里葉基函數(shù)對(duì)原始數(shù)據(jù)進(jìn)行懲罰擬合，通過(guò)廣義交叉驗(yàn)證法(GCV)確定基函數(shù)個(gè)數(shù)為21個(gè)。從圖2可見(jiàn)，武漢空氣中NO2的含量全年都明顯高于其他城市，而咸寧市和孝感市的NO2含量較低?？偟膩?lái)看，湖北省13個(gè)城市NO2含量在冬天要更高一些，并且在2月份的時(shí)候濃度有一個(gè)明顯的低谷區(qū)域。O3情況恰好與NO2相反，夏季空氣中的含量要明顯高于冬季，這是因?yàn)橄募居绕渲形绻饩€較強(qiáng)的時(shí)候，光化學(xué)反應(yīng)會(huì)更加劇烈空氣中O3含量會(huì)更高。13個(gè)城市中大氣中O3含量較低的城市是恩施市。春夏秋三季的O3濃度要高于冬季，但是在7、8兩個(gè)月O3濃度較低，形成了一個(gè)明顯的波谷。

(a)NO2

(b)O3圖 2 各城市NO2與O3的擬合數(shù)據(jù)

2 二元函數(shù)型主成分分析

2.1 二元函數(shù)型主成分

設(shè)有兩個(gè)相互聯(lián)系的函數(shù)型變量(隨機(jī)函數(shù))x(s)和y(s)，s∈T，假設(shè)x(s)和y(s)的N次觀察分別為xi(s)和yi(s)，i=1,2,…,N，s∈T，且二維向量函數(shù)mi(s)=(xi(s),yi(s))′。定義x(s)和y(s)主成分權(quán)函數(shù)為向量函數(shù)ξ=(ξ(x)(s),ξ(y)(s))′，其中ξ(x)(s)和ξ(y)(s)分別表示x(s)與y(s)各自的主成分權(quán)函數(shù)，本文稱之為邊緣權(quán)重函數(shù)。定義兩個(gè)函數(shù)型變量的第j個(gè)主成分得分

(1)

x(s)和y(s)各自的方差函數(shù)記為vxx和vyy，x(s)和y(s)的交叉協(xié)方差函數(shù)記為vxy(s,t)=cov(x(s),y(t))(t,s∈T)。

與一元隨機(jī)函數(shù)的主成分分析過(guò)程類似，二元函數(shù)型主成分分析也可以轉(zhuǎn)化為特征方程的求解問(wèn)題，具體表示為[15]

(2)

(3)

(4)

(5)

2.2 二元主成分權(quán)函數(shù)的求解

(6)

(7)

進(jìn)一步，令b(1)*=W1/2b(1)，b(2)*=W1/2b(2)，γ*=[(b(1)*)T,(b(2)*)T]，并帶入式(6)可解出γ*，之后就可求出b(1)、b(2)，將其帶入ξ(x)(s)和ξ(y)(s)的基函數(shù)展開(kāi)式即可得到x(s)和y(s)的權(quán)函數(shù)ξ(x)(s)和ξ(y)(s)。

3 模型應(yīng)用

3.1 聯(lián)合主成分分析

參與O3光化學(xué)反應(yīng)的氮氧化物以NO2為典型代表，所以在實(shí)際的空氣監(jiān)測(cè)中，O3和NO2大部分時(shí)刻是呈現(xiàn)相互轉(zhuǎn)化的情況。因此相對(duì)于一元函數(shù)主成分分析來(lái)說(shuō)，將這兩個(gè)變量作為一個(gè)整體來(lái)進(jìn)行二元函數(shù)型主成分分析會(huì)更加科學(xué)可靠。

(a)NO2

(b)O3圖 3 NO2和O3的特征函數(shù)

圖3展示了對(duì)數(shù)據(jù)進(jìn)行聯(lián)合主成分后提取的NO2和O3前三個(gè)特征函數(shù)即權(quán)函數(shù)。黑色實(shí)線代表第一主成分權(quán)函數(shù),記為PC1；紅色虛線和綠色虛線分別代表第二、三主成分權(quán)函數(shù)，記為PC2和PC3。NO2的3個(gè)主成分依次解釋了秋季、春夏、秋冬時(shí)期對(duì)大氣中NO2波動(dòng)特征的貢獻(xiàn)；O3的3個(gè)主成分依次解釋了夏秋、春季、秋季時(shí)期對(duì)大氣中臭氧波動(dòng)的貢獻(xiàn)。為了更加清楚地展示權(quán)函數(shù)對(duì)于均值函數(shù)的影響，繪制聯(lián)合主成分權(quán)函數(shù)偏離均值函數(shù)的效果圖(圖4)，即在同一圖中繪制三條曲線，中間由黑色實(shí)心點(diǎn)繪制的曲線為13條擬合曲線的均值函數(shù)曲線，由“+”構(gòu)成的曲線是在均值曲線的基礎(chǔ)上加上權(quán)函數(shù)的適當(dāng)倍數(shù)繪制而成，相對(duì)地由“-”構(gòu)成的曲線是在均值曲線的基礎(chǔ)上減去權(quán)函數(shù)的適當(dāng)倍數(shù)繪制而成，這里倍數(shù)數(shù)值等于權(quán)函數(shù)對(duì)應(yīng)的特征值開(kāi)平方后的數(shù)值。“+” 、“-”形成的曲線可以理解為湖北省13個(gè)城市大氣中NO2和O3全年含量的變化特征，前3個(gè)聯(lián)合主成分的累積貢獻(xiàn)率達(dá)到90.7%。其中：主成分1刻畫(huà)了NO2秋季和O3夏秋季的變動(dòng)趨勢(shì)，解釋了變動(dòng)特征的55.8%，是這兩種污染氣體隨時(shí)間變化的主要方式；主成分2刻畫(huà)了NO2春夏和O3春季的變動(dòng)趨勢(shì)，解釋了變動(dòng)特征的24.9%，是兩種污染氣體隨時(shí)間變化的第二方式；主成分3刻畫(huà)了NO2秋冬和O3秋季的變動(dòng)趨勢(shì)，其解釋能力為10%。

(a)主成分1(1)

(c)主成分2(1)

(d)主成分2(2)

(e)主成分3(1)

(f)主成分3(2)圖 4 NO2與O3的聯(lián)合主成分

3.2 函數(shù)型主成分聯(lián)合變異分析

在二元情況下顯示聯(lián)合主成分的有效方法是構(gòu)造一個(gè)變量與另一個(gè)變量的關(guān)系圖(圖5)。圖中，以NO2和O3兩個(gè)變量的均值作為橫縱坐標(biāo)即可得到聯(lián)合均值曲線，在均值曲線的基礎(chǔ)上加上聯(lián)合主成分權(quán)函數(shù)倍數(shù)，得到聯(lián)合主成分權(quán)函數(shù)對(duì)平均水平的影響情況，線段的角度反映在不同時(shí)間點(diǎn)兩個(gè)變量對(duì)均值曲線影響的強(qiáng)弱區(qū)別。通過(guò)計(jì)算NO2和O3的特征函數(shù)占兩者總特征函數(shù)的比例，可得NO2和O3分別對(duì)于第一聯(lián)合主成分變異程度的貢獻(xiàn)。

第一步觀察紅色數(shù)字組成的曲線。它表示兩變量在區(qū)間內(nèi)的平均水平(即為聯(lián)合均值曲線)，數(shù)字大小表示均值曲線在一年中的走勢(shì)情況。一月初NO2濃度達(dá)到峰值，可以看到NO2濃度達(dá)到最大的時(shí)候反而是臭氧濃度的最小值點(diǎn)，接著NO2濃度降低并在8月達(dá)到最小值，O3濃度在9月取得最大值，但是產(chǎn)生的O3累計(jì)到秋季初才達(dá)到峰值。符號(hào)“+”構(gòu)成的曲線是在均值函數(shù)的基礎(chǔ)上加上聯(lián)合主成分權(quán)函數(shù)的若干倍后得到，并用藍(lán)色線段將其連接。

從2018年10月開(kāi)始，箭頭指向右上方，第一聯(lián)合主成分對(duì)均值函數(shù)有正向影響，6月箭頭方向由右偏到左邊，聯(lián)合主成分對(duì)于NO2的作用由正向變?yōu)樨?fù)向；9月之后第一聯(lián)合主成分對(duì)于NO2的影響又轉(zhuǎn)為正向，并且保持該狀態(tài)直至結(jié)束，在這期間效應(yīng)是逐漸凸顯的，反映在藍(lán)色線段慢慢加長(zhǎng)。通過(guò)計(jì)算可以得到，NO2和O3兩個(gè)變量的變化對(duì)于第一聯(lián)合主成分變異的貢獻(xiàn)率分別為2.36%和97.64%。因此第一聯(lián)合主成分的變異絕大部分來(lái)自與臭氧的變異貢獻(xiàn)。

第二聯(lián)合主成分變異主要來(lái)自NO2，箭頭方向基本指向右方，計(jì)算得NO2和O3對(duì)于第二聯(lián)合主成分的貢獻(xiàn)率分別為70.4%和29.6%。從數(shù)字“1”開(kāi)始(即2018年10月份開(kāi)始到2019年6月)箭頭一直指向右方偏下的位置，這個(gè)過(guò)程中第二聯(lián)合主成分對(duì)臭氧存在微弱的負(fù)向作用；7月份主成分對(duì)O3的負(fù)向影響變?yōu)檎颍Ч廊晃⑷?；該影響?0月轉(zhuǎn)回負(fù)向，并保持直至結(jié)束。整個(gè)過(guò)程中第二聯(lián)合主成分都對(duì)NO2保持較為強(qiáng)烈正向作用。

(a)第一主成分

(b)第二主成分

(c)第三主成分圖 5 前3個(gè)主成分的聯(lián)合變異效果

第三聯(lián)合主成分的變化形式相對(duì)于前兩個(gè)聯(lián)合主成分來(lái)說(shuō)較為復(fù)雜，NO2和O3對(duì)于第三聯(lián)合主成分的貢獻(xiàn)率分別為22.9%和77.1%。曲線剛開(kāi)始有微弱負(fù)向影響作用于O3，但很快轉(zhuǎn)為正向，并且隨著春季的到來(lái)，該影響逐漸加強(qiáng)，NO2在上述過(guò)程中一直是受到聯(lián)合主成分的正向作用；聯(lián)合主成分對(duì)NO2的作用在5月完成由正到負(fù)的轉(zhuǎn)換，而且對(duì)O3的正向影響也變?nèi)?，但是該狀態(tài)維持時(shí)間很短，6月初時(shí)箭頭又指回右上方；8月份及以后，聯(lián)合主成分對(duì)于O3保持負(fù)向影響，對(duì)NO2施加正向作用直至結(jié)束。

通過(guò)對(duì)NO2和O3進(jìn)行二元聯(lián)合主成分分析發(fā)現(xiàn)，春夏兩季NO2的變動(dòng)為其全年變異的主要方式，夏秋兩季則為臭氧變動(dòng)的主要方式。因此根據(jù)全年變異程度的不同，針對(duì)性地選擇時(shí)間段進(jìn)行大氣治理會(huì)得到較好的效果。

3.3 基于函數(shù)型主成分的K-均值聚類

在利用二元函數(shù)型主成分分析對(duì)13個(gè)城市空氣質(zhì)量變化曲線族進(jìn)行分析之后，得到NO2和O3的主要變異方式，但是由于每個(gè)城市大氣中這兩種氣體的污染程度是不同的，因此采用K-Means聚類對(duì)主成分得分對(duì)城市進(jìn)行分類，聚類結(jié)果繪制如圖6所示。

圖 6 K-均值聚類結(jié)果

聚類結(jié)果顯示NO2的污染主要集中在個(gè)別受工業(yè)污染與城市汽車尾氣排放影響比較突出的大中型城市，如武漢、宜昌、鄂州、黃石等地區(qū)。第一類的恩施市和十堰市是空氣質(zhì)量較優(yōu)的，相較于其他11個(gè)城市，這兩個(gè)城市大氣中的NO2和O3都是低于平均值的；第二類屬于居中水平，被聚到該類的5個(gè)城市具有一個(gè)普遍的特點(diǎn)就是大氣中的O3濃度高、NO2濃度低或是居中水平；最后一類的6個(gè)城市，是三類中空氣質(zhì)量最不好的，NO2濃度高、臭氧濃度高或是處于居中水平，治理任務(wù)重。

分類后總體可以概括如下：NO2和O3低；O3高但NO2中等或低；NO2高且O3高或中等。對(duì)于污染最嚴(yán)重的第三類城市來(lái)說(shuō)，需要專注于春夏兩季NO2和夏秋兩季O3的治理，污染程度中等的第二類城市則應(yīng)該先投入到夏秋季O3的治理，由于第一類城市空氣質(zhì)量最優(yōu)，其余城市可參考第一類城市的大氣治理方案，比如恩施地區(qū)在構(gòu)建鄂西生態(tài)旅游圈過(guò)程中注重旅游業(yè)發(fā)展，并且煙氣脫硝使得減排改善效果明顯。

3 結(jié)束語(yǔ)

對(duì)近一年湖北省大氣中NO2和O3濃度數(shù)據(jù)做二元函數(shù)型主成分分析，不僅能夠反映各個(gè)城市各個(gè)時(shí)期的空氣質(zhì)量，而且能夠反映相互聯(lián)系的變量之間的聯(lián)合變異情況。對(duì)聯(lián)合主成分圖和聯(lián)合變異圖進(jìn)行分析后，得出全年中NO2和O3的主要變異時(shí)期，同時(shí)對(duì)主成分得分進(jìn)行聚類發(fā)現(xiàn)，湖北省空氣污染在高原和山地較輕，污染主要集中于平原地區(qū)。根據(jù)分類結(jié)果，再結(jié)合兩個(gè)污染氣體的主要變異時(shí)期，即可對(duì)不同地區(qū)量身定制治理計(jì)劃。