基于聚類(lèi)分析對(duì)我國(guó)31個(gè)主要城市空氣質(zhì)量的研究

一、緒論（一）選題背景和意義近些年，隨著經(jīng)濟(jì)的發(fā)展和社會(huì)的進(jìn)步，人們的物質(zhì)生活水平越來(lái)越高，但在生態(tài)環(huán)境和空氣質(zhì)量方面也付出了相當(dāng)大的代價(jià)。在注重生態(tài)文明建設(shè)強(qiáng)調(diào)綠色低碳發(fā)展的大背景下，研究我國(guó)城市空氣污染情況，空氣質(zhì)量檢測(cè)數(shù)據(jù)的分類(lèi)方法、類(lèi)別特征、優(yōu)劣排序等，挖掘空氣質(zhì)量數(shù)據(jù)隱含的重要信息，尋找改善空氣質(zhì)量的對(duì)策，具有重要的意義。（二）數(shù)據(jù)選取按照《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)（GB3095-2012）》，空氣質(zhì)量的六項(xiàng)主要監(jiān)測(cè)污染物分別為二氧化硫（SO2）、二氧化氮（NO2）、可吸入顆粒物（PM10）、一氧化碳（CO）、臭氧（O3）和細(xì)顆粒物（PM2.5）。因此，本文考慮分析31個(gè)環(huán)保重點(diǎn)城市這7個(gè)具體監(jiān)測(cè)指標(biāo)，分別為SO2年平均濃度x1（μg/m3）、NO2年平均濃度x2（μg/m3）、PM10年平均濃度x3（μg/m3）、CO日均值第95百分位濃度x4（mg/m3）、O3日最大8小時(shí)第90百分位濃度x5（μg/m3）、PM2.5年平均濃度x6（μg/m3），再加上一年內(nèi)空氣質(zhì)量達(dá)到及好于二級(jí)的天數(shù)x7（天），記作G-day。二、聚類(lèi)分析（一）k-means聚類(lèi)1．k-means聚類(lèi)簡(jiǎn)介k-means聚類(lèi)，又稱(chēng)快速聚類(lèi)法，是一種按樣本數(shù)據(jù)的相似程度將樣本分類(lèi)的方法.需要提前設(shè)定聚類(lèi)個(gè)數(shù)。相似程度的度量一般按樣本間的距離劃分，兩個(gè)樣本間距離越小，樣本的相似程度越高。2．k-means聚類(lèi)結(jié)果仍然以上述2019年的空氣質(zhì)量數(shù)據(jù)為例，讀入數(shù)據(jù)，并對(duì)數(shù)據(jù)畫(huà)出散點(diǎn)圖，從樣本點(diǎn)的分布的集中程度以看出兩個(gè)變量之間的線性相關(guān)性的強(qiáng)弱。從上圖1可以看出，PM2.5與PM10的濃度呈正相關(guān)，而PM2.5和PM10的濃度均與一年內(nèi)空氣質(zhì)量達(dá)到及好于二級(jí)天數(shù)（圖中記為day）呈較強(qiáng)的負(fù)相關(guān)，PM2.5和PM10的濃度越高，好天氣數(shù)越少，空氣質(zhì)量越差。在實(shí)際生活中，表現(xiàn)為霧霾嚴(yán)重，空氣可見(jiàn)度低，容易引發(fā)呼吸系統(tǒng)疾病，危害人們的身體健康。同時(shí)可見(jiàn)，G-day與二氧化氮年平均濃度呈負(fù)相關(guān)。NO2的排放源于汽車(chē)尾氣，日常生活中，汽車(chē)尾氣排放的越多，空氣質(zhì)量越差。以下利用k-means函數(shù)對(duì)31個(gè)城市2019年的空氣質(zhì)量數(shù)據(jù)進(jìn)行k-means聚類(lèi)，令聚類(lèi)數(shù)k由小到大逐漸增加，同時(shí)計(jì)算類(lèi)間平方和與總平方和的比值。該比值越大，說(shuō)明兩個(gè)類(lèi)間的差異越大，則類(lèi)內(nèi)部的差異越小，聚類(lèi)的效果越理想。利用R程序可計(jì)算并構(gòu)造類(lèi)間平方和與總平方和的比值序列，當(dāng)k=3之后，組間平方和占比的增幅就很小了，變化幅度不超過(guò)0.01。組間平方和占比趨近于1。作聚類(lèi)分析時(shí)既要使組間平方和占比盡量大，又要讓聚類(lèi)數(shù)不能太大（否則太雜亂），權(quán)衡考慮后取k=3為最佳聚類(lèi)數(shù)。利用聚類(lèi)所得的3類(lèi)城市的均值點(diǎn)的最后兩列數(shù)據(jù)，將空氣質(zhì)量達(dá)到及好于二級(jí)的天數(shù)（即最后一列day）對(duì)三個(gè)聚類(lèi)結(jié)果進(jìn)行排序，如圖2所示。（二）PAM聚類(lèi)1．PAM聚類(lèi)簡(jiǎn)介PAM聚類(lèi)針對(duì)k-means聚類(lèi)法的不足，用類(lèi)中心點(diǎn)代替質(zhì)心。PAM算法鑒于k均值算法的這一不足，選取一個(gè)類(lèi)中位置最接近質(zhì)心的數(shù)據(jù)點(diǎn)作為類(lèi)中心點(diǎn)，來(lái)代替質(zhì)心，其他數(shù)據(jù)點(diǎn)就根據(jù)距離被分到離數(shù)據(jù)中心點(diǎn)中最近的那一類(lèi)。反復(fù)根據(jù)一類(lèi)的數(shù)據(jù)點(diǎn)坐標(biāo)來(lái)調(diào)整類(lèi)中心點(diǎn)，解決了離群點(diǎn)對(duì)結(jié)果的影響。2．PAM聚類(lèi)結(jié)果本文運(yùn)用R軟件對(duì)數(shù)據(jù)做PAM聚類(lèi)，運(yùn)用cluster程序包中的pam函數(shù)進(jìn)行聚類(lèi)。從各類(lèi)中選取的距離質(zhì)心位置最近的中心點(diǎn)列表如下表1。表1距離質(zhì)心位置最近的中心點(diǎn)情況上面所得結(jié)果是各類(lèi)的中心點(diǎn)坐標(biāo)，分別為成都、濟(jì)南、昆明。聚類(lèi)中心點(diǎn)是各類(lèi)城市的代表，從中心點(diǎn)就可以大致看出各類(lèi)的空氣質(zhì)量。將城市空氣質(zhì)量由好到差排序，結(jié)果如下。優(yōu)（第三類(lèi)）：福州、南寧、?？?、貴陽(yáng)、昆明、拉薩、西寧。良（第一類(lèi)）：北京、呼和浩特、沈陽(yáng)、長(zhǎng)春、哈爾濱、上海、南京、杭州、合肥、南昌、武漢、長(zhǎng)沙、廣州、重慶、成都、蘭州、銀川、烏魯木齊。差（第二類(lèi)）：天津、石家莊、太原、濟(jì)南、鄭州、西安。（三）兩種聚類(lèi)分析結(jié)果比較對(duì)幾種聚類(lèi)結(jié)果觀察可以發(fā)現(xiàn)，按照城市的空氣質(zhì)量好壞分類(lèi)，幾種聚類(lèi)結(jié)果具有一致性，大致客觀反映了全國(guó)31個(gè)環(huán)保重點(diǎn)城市空氣質(zhì)量狀況的分類(lèi)、空間分布和空氣檢測(cè)指標(biāo)的實(shí)際情況。k-means聚類(lèi)就不變作圖，但k-means聚類(lèi)能快速聚類(lèi)，節(jié)省時(shí)間，給出每一類(lèi)中的具體成員。PAM聚類(lèi)是一種基于質(zhì)心的劃分型聚類(lèi)算法，與K-Means聚類(lèi)的主要不同在于：第一，距離測(cè)度采用絕對(duì)聚類(lèi)。聚類(lèi)目標(biāo)是找到類(lèi)內(nèi)絕對(duì)距離之和最小下的類(lèi)；第二，增加判斷本次迭代類(lèi)質(zhì)心合理性的步驟。三、結(jié)束語(yǔ)本文對(duì)全國(guó)31個(gè)環(huán)保重點(diǎn)城市的空氣質(zhì)量數(shù)據(jù)通過(guò)R軟件做了聚類(lèi)分析。所做的兩種聚類(lèi)分析的結(jié)果大體一致。6種主要的空氣污染物中，SO2與CO之間相關(guān)性相對(duì)較強(qiáng)，大致主要來(lái)源于燃煤的燃燒，NO2、PM10、O3與PM2.5之間的相關(guān)性較強(qiáng)，大致來(lái)源于汽車(chē)尾氣及揚(yáng)塵污染。從31個(gè)環(huán)保重點(diǎn)城市空氣質(zhì)量數(shù)據(jù)的雙坐標(biāo)圖中能夠看出每個(gè)城市