李冰月 周新宇 王楠 梁冰

摘要：為了實現(xiàn)對白塔寺歷史文化區(qū)的精細化管理，在該區(qū)域內(nèi)部署了60套CityGrid -體化傳感器進行實時監(jiān)測，選取2017年11月份二氧化硫、二氧化氮、一氧化碳、PM2.5、PMio和臭氧的監(jiān)測數(shù)據(jù)作為研究對象，通過與北京監(jiān)測站——西城官園的監(jiān)測數(shù)據(jù)對比發(fā)現(xiàn)，兩地的各個指標變化趨勢大體一致，說明傳感器滿足監(jiān)測的精度要求。通過線性插值法和主成分分析分別構(gòu)造空氣質(zhì)量指數(shù)和空氣質(zhì)量綜合指數(shù)，發(fā)現(xiàn)兩者的相關(guān)系數(shù)達到0.831，在一定程度上都能反映白塔寺街區(qū)的空氣質(zhì)量狀況。與官園監(jiān)測站的空氣質(zhì)量指數(shù)對比發(fā)現(xiàn)空氣質(zhì)量綜合指數(shù)和官園空氣質(zhì)量指數(shù)的相關(guān)性更是高達0.883，特別是在11月2日、7日和22日指數(shù)大幅度下降時，空氣質(zhì)量綜合指數(shù)和官園空氣質(zhì)量指數(shù)更加吻合。綜上可知主成分分析比線性插值法更適合于白塔寺街區(qū)的空氣質(zhì)量研究，并建立新的評價等級，以切實反映白塔寺街區(qū)的空氣狀況。

關(guān)鍵詞：空氣質(zhì)量;歷史文化保護區(qū);線性插值法;主成分分析;評價等級

中圖分類號：X8 20.2

文獻標識碼：A

文章編號：1674-9944（2018）4-0066-05

1 引言

歷史文化街區(qū)是城市中寶貴的歷史文化遺產(chǎn)，它承載了這個城市的傳統(tǒng)意蘊和精神內(nèi)涵。北京西城區(qū)內(nèi)的阜內(nèi)白塔寺歷史文化街區(qū)（簡稱白塔寺街區(qū)）是北京舊城的典型樣本，坐落其中的妙應白塔寺是我國現(xiàn)存最早、最大的喇嘛塔，周邊的宮門口頭條、白塔寺東夾道、蘇蘿卜胡同等構(gòu)成了富有特色的“胡同風情”。但近幾年來，由于市民機動車保有量的增加，機動車與非機動車混行、胡同亂停車現(xiàn)象普遍，環(huán)境衛(wèi)生令人堪憂，隨著“白塔寺再生計劃”的實施，關(guān)注其發(fā)展建設的人越來越多。

目前我國大氣污染正由以SOz和PMio為特征的煤煙型污染發(fā)展成多種污染物復合型污染[1-4]，單一污染物的研究[5-7]已不能完全解釋我國城市大氣污染狀況，空氣質(zhì)量的問題變得十分嚴峻。國外學者對空氣質(zhì)量的研究趨勢是采用多種方法進行綜合分析，如Powell和Lee[8]提出一種通用的貝葉斯層次框架，Bruno和Cocchic9]通過連續(xù)地應用一些聚合函數(shù)來合成空氣質(zhì)量指數(shù)。國內(nèi)傳統(tǒng)的計算方法是《環(huán)境空氣質(zhì)量標準》（ GB 3095 -2012）[10]里的線性插值法，如對保定市、秦皇島市和北京市的空氣質(zhì)量的研究[11—13]。潘本鋒籌[14]發(fā)現(xiàn)此方法在應用中存在一些問題，并提出相應的改進建議。目前常用的空氣質(zhì)量綜合評價方法包括綜合指數(shù)法、灰色聚類法、模糊綜合評價法、主成分分析法等，如姜新華等[15]基于主成分分析法對呼和浩特市的空氣質(zhì)量進行研究，得出主要影響因素為PM₁₀、CO和S02。張茹[16]采用層次分析和主成分方法對徐州市大氣監(jiān)測數(shù)據(jù)進行研究，發(fā)現(xiàn)兩種方法均可作為質(zhì)量評價的方法。

由于國家監(jiān)測站點有限且間距較大，很難做到監(jiān)測指標的全區(qū)域覆蓋。鑒于國測站造價高且分布不夠精細等缺陷，北京市科委借示范項目系統(tǒng)建設的機會，在白塔寺街區(qū)內(nèi)部署了60套CityGrid -體化傳感器，實時采集多項大氣環(huán)境數(shù)據(jù)，以便觀察白塔寺街區(qū)的實時運行狀態(tài)。作為階段性成果，本文選取北京市監(jiān)測站中距白塔寺街區(qū)最近的官園監(jiān)測點作為對照，分別采用線性插值法和主成分分析進行統(tǒng)計建模，經(jīng)分析發(fā)現(xiàn)主成分分析一方面能體現(xiàn)白塔寺街區(qū)的區(qū)域特點，另一方面和官園監(jiān)測站更加吻合，說明主成分分析比線性插值法更適合白塔寺街區(qū)的空氣質(zhì)量建模。

2 主成分分析

很多情況下，在所研究課題的眾多變量之間存在一定的相關(guān)關(guān)系，導致部分信息的重疊，同時變量個數(shù)過多也會增加研究的難度。主成分分析恰好能解決此類問題，其主要目的便是為了降維，基本思想是通過構(gòu)造原變量的線性組合，以產(chǎn)生一系列互不相關(guān)的新變量，從中選出為數(shù)較少的新變量并使它們含有原變量所提供的大部分信息，從而使得用這幾個新變量來分析問題成為可能。分析步驟分5步：首先，對原始數(shù)據(jù)進行標準化，以消除量綱的不同帶來的影響;接著，計算樣本協(xié)方差陣或相關(guān)系數(shù)陣;然后，求樣本協(xié)方差陣或相關(guān)系數(shù)陣的特征值和特征向量，從而計算各變量的貢獻率和累計貢獻率;再后，確定主成分的個數(shù)（一般要求累計貢獻率大于80%即可）;最后，對所選擇的主成分給出合理解釋，計算主成分得分和綜合得分。

3 實例應用與分析

3.1 數(shù)據(jù)來源與預處理

白塔寺街區(qū)的監(jiān)測數(shù)據(jù)來自該街區(qū)范圍內(nèi)部署的60套CityGrid -體化傳感器，該設備基于定電位電解傳感器原理及激光檢測技術(shù)監(jiān)測污染氣體、光散射原理檢測粉塵的技術(shù)，可監(jiān)測19項參數(shù)，包括：光照強度、紫外線指數(shù)、溫度、濕度、風向、風速、PM2.5、PM10、大氣壓強、噪聲、S02、NO2、CO、C02、03、VOC濃度、人流量、車流量、車速分析。

另一方面，了解到當前北京市西城區(qū)共建有2個用于采集空氣質(zhì)量指標數(shù)據(jù)的監(jiān)測站，選擇距白塔寺街區(qū)最近的官園監(jiān)測站，從北京市環(huán)境保護監(jiān)測中心網(wǎng)站上可收集北京市各監(jiān)測站的空氣質(zhì)量時報數(shù)據(jù)（S02、NO2、CO、PM2.5、PM10、03、AQI）.

白塔寺街區(qū)和官園監(jiān)測站在地圖上的位置如圖1所示。

選取影響空氣質(zhì)量的六項主要污染物指標作為研究對象，首先對數(shù)據(jù)進行預處理，將60套監(jiān)測設備的數(shù)據(jù)匯總，然后以小時為單位取平均值，得出白塔寺街區(qū)2017年11月的整體污染物監(jiān)測情況，共收集720條數(shù)據(jù)，部分數(shù)據(jù)如表1。

3.2 白塔寺街區(qū)污染物指標的特征分析

首先對白塔寺街區(qū)污染物指標的監(jiān)測數(shù)據(jù)進行逐一分析，并通過與官園的數(shù)據(jù)對比，得出白塔寺街區(qū)污染物指標的區(qū)域特征。通過R語言繪圖，觀察兩地污染物濃度的變化規(guī)律，見圖2。

圖2為PM2.5、PMio、S02、NO2、CO、03的時間序列圖，水平直線代表空氣質(zhì)量指數(shù)為100時對應的各污染物濃度限值。從圖中可以看出，白塔寺街區(qū)在11月期間主要污染物為PM2.5，而PM10、S02和03只出現(xiàn)短暫的超標，NO2和CO濃度一直未超標。另一方面PM2.5和PMio波動較大，且具有同升同降的現(xiàn)象，而S02、NO2、CO、03的濃度變化相對穩(wěn)定，且CO濃度有逐漸降低的趨勢，整體上看，白塔寺街區(qū)的空氣質(zhì)量正在逐漸轉(zhuǎn)好。經(jīng)對比可知，白塔寺街區(qū)和官園監(jiān)測點的為PM2.5、PM10、S02、CO的變化趨勢保持高度一致，考慮到官園距白塔寺街區(qū)僅相隔1.2 km，此現(xiàn)象符合實際情況，說明傳感器在精度上滿足監(jiān)測要求。

3.3 影響空氣質(zhì)量因素的主成分分析法

對污染物指標SO2 （Xl）、NO2 （Xz）、CO（ X3）、PM2.5（X4）、PM10 （X5）和03（X6）進行主成分分析，首先計算各變量之間的Pearson相關(guān)系數(shù)，通過R語言編程實現(xiàn)，計算結(jié)果見表2。

從表2可以看出CO與N02、PM25、PMio的相關(guān)系數(shù)均在0.5以上，存在較強的相關(guān)性，另外PMz.s和PMio的相關(guān)系數(shù)高達0.996，這是因為PMio濃度中包含了PM2.5，兩者的強相關(guān)與實際相符。接下來計算相關(guān)系數(shù)陣的特征值和特征向量，從而得到各主成分的貢獻率和累計貢獻率見表3。

依據(jù)表3的結(jié)果及累計貢獻率大于80%的原則，選取前3個主成分進行分析，表達式如下：

Yl=- 0.0614Xl+0.3688X2 - 0.4920X3 - 0.5492X4 -0.5490X 5 +0.1231X 6

Y2=0.8558X1+ 0.4416X2- 0.0 918X 3+0.1216X 4+0.1196X5 -0.1872X 6

Y3=-0. 2169X1+0.0637X 2+0.0267X3 - 0.094]X 4 -0.0794X 5 -0.9659X 6

這里總共選取了三個主成分，第一主成分包含了47.1%的信息量，主要代表因素是PM2.5和PMio，反映了空氣中微小顆粒物的濃度;第二主成分包含17.8%的信息量，代表因素為二氧化硫;第三主成分包含了16.3%的信息量，代表因素為臭氧。由此可見，空氣質(zhì)量的主要影響因子包括PM25、二氧化硫和臭氧。

以各主成分相應的貢獻率作為權(quán)重求和，即Z=0.4710 Y1+0.1778Y2+0.1629Y3，將Yl、Y2、Y3代入可得

Z=0.0879X1+0.2626X 2 - 0.2437X 3 - 0.2524X 4-0.2502X 5 -0.1326X 6

將上式的系數(shù)進行標準化得出污染物指標的權(quán)重，可得空氣質(zhì)量綜合指數(shù)AQCI模型為：

F=0.0276X1+0.2463X 3+0.2121X 3+0.2275X 4-l-0.2237X5+0.0629X6

3.4 影響空氣質(zhì)量因素的線性捶值法

本節(jié)針對白塔寺街區(qū)的污染物指標監(jiān)測數(shù)據(jù)，根據(jù)《環(huán)境空氣質(zhì)量指數(shù)（AQI）技術(shù)規(guī)定（試行）》（HJ633-2012）的空氣質(zhì)量分指數(shù)IAQI及對應的污染物項目濃度限值，通過構(gòu)造S02、NO2、CO、PM2.5、PM10、03的1h平均，PM2.5和PM10的24 h滑動平均和03的8h滑動平均這9個指標進行建模，根據(jù)線性插值法可知，污染物指標P的空氣質(zhì)量分指數(shù)IAQIP計算公式為：

其中IAQI表示空氣質(zhì)量分指數(shù)，n表示污染物項目的數(shù)量。

將監(jiān)測數(shù)據(jù)帶入式（1），用R語言統(tǒng)計軟件編程實現(xiàn)，可得到空氣質(zhì)量指數(shù)AQI和首要污染物，其中11月份白塔寺街區(qū)首要污染物的頻率見表4。

從表4可以看出，白塔寺街區(qū)首要污染物為PM2.5的頻率為0.947，其次是03，頻率為0.031，說明PM2.5是白塔寺街區(qū)空氣污染的主要來源。

3.5 結(jié)果分析

首先分析白塔寺街區(qū)AQCI、白塔寺街區(qū)AQI與官園AQI的相關(guān)性，計算三者之間的Pearson相關(guān)系數(shù)見表5。

從表5可以看出白塔寺街區(qū)AQCI、白塔寺街區(qū)AQI的相關(guān)系數(shù)為0.831，說明兩種方法的相關(guān)性很強，在一定程度上都能用來作為白塔寺街區(qū)的空氣質(zhì)量分析模型。而當與官園AQI進行對比時，發(fā)現(xiàn)AQCI與官園AQI的相關(guān)性更強，為0.883，鑒于白塔寺街區(qū)距官園僅1.2 km，而國測站要求監(jiān)測能反映周圍5 km的空氣狀況，因此主成分分析更適用于白塔寺街區(qū)的空氣質(zhì)量建模，下面畫出三者的折線圖，結(jié)果見圖3。

從圖3可以看出與線性插值法相比，主成分分析得出的AQCI與官園監(jiān)測站AQI更加貼合，特別是在2、7和22日空氣質(zhì)量指數(shù)大幅度下降的時候，AQCI與官園的下降斜率基本相同，而插值法得到的AQI下降的相對緩慢，鑒于線性插值法僅以最大污染物的空氣質(zhì)量分指數(shù)為準，而實際生活中，次大污染物對空氣質(zhì)量的影響也不容小視，而主成分分析方法能綜合各個污染物指標分析計算，更適合白塔寺街區(qū)的空氣質(zhì)量建模。另一方面，無論是線性插值還是主成分分析，得出的白塔寺街區(qū)空氣質(zhì)量數(shù)據(jù)在19～22日都比官園的低，這是因為據(jù)中國新聞網(wǎng)11月20日報道：根據(jù)環(huán)保局空氣質(zhì)量預報和實時監(jiān)測情況，目前，京津冀地區(qū)污染擴散條件不利，北京市受區(qū)域性逆溫、西南風傳輸及疊加本地排放影響，預計今日夜間至明日傍晚空氣質(zhì)量將達到重度污染水平，故官園污染較嚴重。由于主成分分析得出的AQCI相對較小，將其擴大兩倍，見圖4。

從圖4發(fā)現(xiàn)，擴大兩倍的AQCI與官園的相似性更加明顯，且基本量綱相同，據(jù)此根據(jù)《環(huán)境空氣質(zhì)量指數(shù)（AQI）技術(shù)規(guī)定（試行）》（HJ633 - 2012）的分級方法，將其標準縮小一半，制定為白塔寺街區(qū)的評價等級，見表6。

4 結(jié)論

通過分析發(fā)現(xiàn)，用低成本的傳感器設備監(jiān)測的白塔寺街區(qū)污染物指標數(shù)據(jù)和北京官園的數(shù)據(jù)變化趨勢基本保持一致，在精度上達到監(jiān)測的要求，同時白塔寺街區(qū)具有自己的區(qū)域特征，如PM2.5和PM10濃度值相對較低，其它污染物濃度變化較平穩(wěn)。分別采用線性插值法和主成分分析法得出空氣質(zhì)量指數(shù)和空氣質(zhì)量綜合指數(shù)，發(fā)現(xiàn)兩者高度相關(guān)，相關(guān)系數(shù)達到0.831。但主成分分析方法與官園空氣質(zhì)量指數(shù)的變化幅度更加貼合，相關(guān)系數(shù)也比較高，為0.883。且受區(qū)域性逆溫、西南風傳輸及疊加本地排放影響，官園的污染較重，而無論是線性插值法還是主成分分析法，都能準確反映這一變化，得出的白塔寺街區(qū)空氣質(zhì)量數(shù)據(jù)都比官園的低。整體上看，主成分分析法一方面能體現(xiàn)白塔寺街區(qū)的區(qū)域特點，另一方面和官園監(jiān)測站更加吻合，因此主成分分析法比線性插值法更適合白塔寺街區(qū)的空氣質(zhì)量建模，最后為空氣質(zhì)量綜合指數(shù)制定新的評價等級，以切實反映白塔寺的空氣狀況。

參考文獻：

[l]Chen C K， Yao X H.Air pollution in mega cities in China [J]. Atmospheric Environment， 2008， 42（1）：1～42.

[2]Deng X J，Tie X X，Wu D，et al.Long-term trend of visibilityand its characterizations in the Pearl River Delta （PRD） region，China [J]. Atmospheric Environment， 2008， 42（7）：1424～1435.

[3]Gao L N， Jia G S，Zhang R J，et al.Visual range trends in theYangtze River Delta region of China during 1981- 2005 [J]. Journal of the Air&Waste Management Association， 2011， 61（8）：843～849.

[4]許亞宣，李小敏，于華通，等，邯鄲市大氣復合污染特征的監(jiān)測研究[J].環(huán)境科學學報，2015，35（9）：2710～2722.

[5]奇奕軒，胡 君，張鶴豐，等，北京市郊區(qū)夏季臭氧重污染特征及生成效率[J].環(huán)境科學研究，2017，30（5）：663～671.

[6]夏澤群，范小莉，黃志炯，等，國內(nèi)外PM2.5源譜對比及其對空氣質(zhì)量模擬效果的影響[J].環(huán)境科學研究，2017，30（3）：359～367.

[7]Chen P F， Bi X H， Zhang J Q， et al_Assessment of heavy metalpollution characteristics and human health risk of exposure to ambient PMz.s in Tianjin， China [J]. Particuology， 2015， 20（3）：104～l09.

[8]Powell H， Lee D.Modelling spatial variability in concentrations ofsingle pollutants and composite air quality indicators in healtheffects studies [J]. Journal of the Royal Statistical Society， 2014，177（3）：607～623.

[9]Bruno F， Cocchi D.A unified strategy for building simple air quality indices [J]. Environmetrics， 2010， 13（3）：243～261.

[10]中華人民共和國環(huán)境保護部，中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局.環(huán)境空氣質(zhì)量標準：GB 3095-2012 [S].北京：中國環(huán)境科學出版社，2012.

[11]趙靜，齊芳，張輝，等.保定市空氣質(zhì)量現(xiàn)狀及綜合治理措施[J].河北建筑工程學院學報，2014，32（1）：55～59.

[12]馮建社，張婷，張明月，等.秦皇島市空氣質(zhì)量指數(shù)AQI現(xiàn)狀分析[J].中國環(huán)境管理干部學院學報，2014，24（5）：65～69.

[13]劉杰，楊鵬，呂文生，等.基于北京市6類污染物的環(huán)境空氣質(zhì)量評價方法[J].安全與環(huán)境學報，2015，15 （1）：310～315.

[14]潘本鋒，宮正宇，王帥，等.環(huán)境空氣質(zhì)量指數(shù)在應用中存在的問題及建議[J].中國環(huán)境監(jiān)測，2015，31（1）：64～67.

[15]姜新華，薛河儒，張存厚，等.基于主成分分析的呼和浩特市空氣質(zhì)量影響因素研究[J].安全與環(huán)境工程，2016，23 （1）：75 ～79.

[16]張茹，張學楊，陸洪光，等.基于層次分析和主成分分析的城市空氣質(zhì)量評價——以徐州市為例[J].安全與環(huán)境工程，2017，24（3）：103～107.