陳謙

摘 ?要：2020年初，中國面臨新型冠狀病毒（COVID-19）的威脅，為了對抗突然爆發(fā)的疫情，中國大陸地區(qū)采取了積極的應(yīng)對措施，其中包括持續(xù)數(shù)十日的個(gè)人隔離與城市隔離【1】，隔離期間，工業(yè)、交通等設(shè)施停止或以低功率運(yùn)行，企業(yè)暫緩復(fù)工，大范圍的交通出行降低，對于大氣質(zhì)量有一定的影響，本文以鹽城市近年的前三月的相關(guān)空氣質(zhì)量評價(jià)指標(biāo)的數(shù)據(jù)為基礎(chǔ)，運(yùn)用R語言進(jìn)行處理，從而獲得疫情對于大氣污染的直觀影響。

關(guān)鍵詞：大氣污染;R語言

1：介紹

對于大氣質(zhì)量的評價(jià)常常有PM2.5、PM10、SO2等指標(biāo)，其數(shù)據(jù)量較大處理過程較為繁瑣，并且AQI指數(shù)只適合用來表征短期內(nèi)的空氣質(zhì)量【2】，且考慮到AQI指數(shù)在處理顆粒物時(shí)所暴露的問題，所以對于分析擁有較大的時(shí)間跨度的大氣質(zhì)量變化，R語言是一個(gè)很好的工具。R語言【3】（The R Programming）誕生與1980年左右，其前身為貝爾實(shí)驗(yàn)室開發(fā)的S語言，目前R語言（以下簡稱R）廣泛的用于統(tǒng)計(jì)計(jì)算、統(tǒng)計(jì)畫圖等領(lǐng)域，由于R開源且具有豐富的程序包，在處理大氣污染相關(guān)的數(shù)據(jù)時(shí)具有速度快、操作簡單、分析準(zhǔn)確等其他軟件不及的優(yōu)勢。

本文將以鹽城市近五年前三個(gè)月的PM2.5、PM10、SO2、CO、NO2以及O3的濃度作為基礎(chǔ)數(shù)據(jù)，并運(yùn)用到R中包含的“ggplot2”、“OpenAir”、“stats”、“randomForest”程序包進(jìn)行分析，同時(shí)本文將具體的探究R語言在執(zhí)行環(huán)境數(shù)據(jù)時(shí)的使用，并不涉及基礎(chǔ)的R的使用知識。

2：方法過程

大氣污染的數(shù)據(jù)來自https：//www.aqistudy.cn/，導(dǎo)入數(shù)據(jù)的方法有很多，可以將網(wǎng)頁中的表格復(fù)制到Excel中，再將其導(dǎo)入到R或者可以使用“rvest”包進(jìn)行網(wǎng)絡(luò)爬蟲直接讀取數(shù)據(jù)，由于爬蟲的代碼相比于前者較為復(fù)雜所以我們采取第一種方法即讀取Excel文件，R自帶的函數(shù)只可以讀取以“.csv”為后綴的文件所以我們此時(shí)需要加載“openxlsx”包才可以讀取以“.xlsx”為后綴的文件，代碼如下圖【圖1】。

將數(shù)據(jù)導(dǎo)入后可以先對原始數(shù)據(jù)做初步處理，我們先以時(shí)間為依據(jù)將數(shù)據(jù)分組進(jìn)行處理，得到各個(gè)污染物數(shù)值的分布，分布圖與代碼見下圖【圖2】。

根據(jù)上述圖像我們可以看出，2020年一、二、三月份的污染物濃度同前四年相比散點(diǎn)的位置較接近X軸，說明污染物的濃度低于前四年所測的數(shù)據(jù)，從直觀上觀測我們可以得出污染物濃度降低，大氣質(zhì)量改善的結(jié)論，但是從科學(xué)的角度來說這種初步的、定性的結(jié)論是粗糙的，下面我們將進(jìn)一步的從定量的角度來說明這一變化。

我們采用t-檢驗(yàn)來比較2020年的數(shù)據(jù)與其他四個(gè)月份的區(qū)別 ，一般的t-檢驗(yàn)需要考慮兩樣本的方差是否相等，在這里我們采用Welch兩樣本t-檢驗(yàn)來檢驗(yàn)不同方差的樣本來節(jié)約分析時(shí)的工作量【4】，由于t-檢驗(yàn)要求只有兩類樣本所以我們依據(jù)年份將數(shù)據(jù)分為a組（年份為2020）和b組（年份為2020年以外的年份）來運(yùn)行t-檢驗(yàn)，得到的結(jié)果與相關(guān)的代碼如下圖所見【圖3】。

t-檢驗(yàn)顯示各污染物的p值全部超過0.05，在統(tǒng)計(jì)學(xué)上及認(rèn)為此次假設(shè)不具備顯著性【5】，表明2020年污染物數(shù)據(jù)與其他年份存在顯著性差異即2020年大氣質(zhì)量的平均水平與往年同時(shí)段的平均水平存在差異，同時(shí)從t檢驗(yàn)所示的各污染物平均數(shù)據(jù)可以看出2020的污染物平均水平低于其余年份，所以從我們從統(tǒng)計(jì)的角度可以初步確定2020年大氣質(zhì)量優(yōu)于往年，但是由于數(shù)據(jù)選取的不足，我們難以排除偶然因素的影響，所以對于數(shù)據(jù)我們將做進(jìn)一步的處理，這里我們要介紹一種統(tǒng)計(jì)學(xué)和機(jī)器學(xué)習(xí)中的新的模型——隨機(jī)森林【6】（randomForest），隨機(jī)森林是指利用多棵決策樹對樣本進(jìn)行訓(xùn)練并預(yù)測的一種分類器，在機(jī)器學(xué)習(xí)中，隨機(jī)森林是一個(gè)包含多個(gè)決策樹的分類器，并且其輸出的類別是由個(gè)別樹輸出的類別的眾數(shù)而定。Leo Breiman和Adele Cutler發(fā)展出推論出隨機(jī)森林的算法。而 "Random Forests" 是他們的商標(biāo)。這個(gè)術(shù)語是1995年由貝爾實(shí)驗(yàn)室的Tin Kam Ho所提出的隨機(jī)決策森林（random decision forests）而來的。這個(gè)方法則是結(jié)合 Breimans 的 "Bootstrap aggregating" 想法和 Ho 的"random subspace method"以建造決策樹的集合。關(guān)于隨機(jī)森林的原理由于和本文無關(guān)就不做過多的說明了，我們需要了解的是隨機(jī)森林在處理大量的變數(shù)以及在內(nèi)部對于一般的化后的誤差產(chǎn)生不偏差的估計(jì)時(shí)具有很大的優(yōu)勢，所以在本文我們選取隨機(jī)森林為進(jìn)一步處理的模型，首先我們根據(jù)六類污染物的濃度做加和運(yùn)算，用得到的數(shù)據(jù)來表示大氣質(zhì)量的優(yōu)劣程度，在選擇判斷數(shù)時(shí)我們參考AQI的計(jì)算方法得到AQI為優(yōu)情況下各污染物濃度并考慮到鹽城市本地全年的平均大氣質(zhì)量，我們將判斷數(shù)設(shè)為217，當(dāng)加和得到的數(shù)小于217時(shí)我們認(rèn)為大氣質(zhì)量為優(yōu)（用T表示），反之則為劣（用F表示），這樣我們的原始數(shù)據(jù)就完成了二分類的改造，將其用于隨機(jī)森林的模型之中，需要注意的是由于我們改變了原始數(shù)據(jù)所以我們需要重新導(dǎo)入數(shù)據(jù)，同時(shí)由于隨機(jī)森林在標(biāo)簽變量為numeric與factor時(shí)執(zhí)行的分別為回歸與分類，所以我們還需將標(biāo)簽數(shù)據(jù)轉(zhuǎn)變?yōu)閒actor類型，確定執(zhí)行的為分類時(shí)我們還需將原始數(shù)據(jù)分為訓(xùn)練數(shù)據(jù)與測試數(shù)據(jù)，我們先用訓(xùn)練數(shù)據(jù)運(yùn)行隨機(jī)森林接著用測試數(shù)據(jù)對模型進(jìn)行預(yù)測從而判斷分類的準(zhǔn)確性，當(dāng)這些前期工作都就緒時(shí)我們就可以執(zhí)行隨機(jī)森林分類了，運(yùn)行后的結(jié)果與運(yùn)行代碼如下圖所見【圖4】【圖5】【圖6】。

圖4表明在大氣質(zhì)量為劣的分組中我們認(rèn)為進(jìn)行的分組與隨機(jī)森林進(jìn)行的分類相比誤差在0.003349282，而在大氣質(zhì)量為優(yōu)的分組中錯(cuò)誤率為0.03703704，我們進(jìn)行的分類與機(jī)器進(jìn)行的分類誤差在5%以內(nèi)，這表明我們進(jìn)行的分類具有較強(qiáng)的準(zhǔn)確性，與測試數(shù)據(jù)的對比也證明了這一點(diǎn)，此時(shí)通過進(jìn)一步的數(shù)據(jù)處理我們將數(shù)據(jù)誤差控制在了5%以內(nèi)，基本排除了由于數(shù)據(jù)偶然性造成的誤判，剩余的工作就非常的簡單了，我們只需統(tǒng)計(jì)出2020年份中T的占比并將其與其他年份T的占比進(jìn)行比較，這一需求在其他軟件中也可以簡單的實(shí)現(xiàn)，在這里我們還是使用R進(jìn)行處理，結(jié)果及代碼如下圖所示【圖7】。

我們可以看出2016——2019四個(gè)年份中質(zhì)量為T的天數(shù)為98約占總數(shù)的0.2715，2020年質(zhì)量為T的天數(shù)為57約占總數(shù)的0.6264，后者約為前者的三倍，至此我們可以從定量的角度得到結(jié)論，疫情期間大氣質(zhì)量有較好的改善。

3：解釋與思考

大氣中污染物的類型大致有兩種，即煤煙型與交通型，前者的主要污染物為二氧化硫、氮氧化物、懸浮物、PM0、一氧化碳，后者的污染物主要有氮氧化物、有機(jī)物，這兩種不同的污染類型所發(fā)生時(shí)間也不同，前者多發(fā)生在冬季，后者發(fā)生于夏季，本次分析選取的數(shù)據(jù)除臭氧外都為煤煙型污染物且因?yàn)橐咔榘l(fā)生在冬季所以同其他年份同時(shí)段的數(shù)據(jù)相比有較為明顯的改善，而由于低空臭氧的產(chǎn)生往往都是由交通型污染產(chǎn)生的光化學(xué)煙霧二次反應(yīng)生成的，所以在交通型污染不占主導(dǎo)的情況下，臭氧的數(shù)據(jù)與往年同時(shí)段相比，差異較之其他污染物略微不明顯。

與此同時(shí)在本實(shí)驗(yàn)中沒有采用比較AQI的數(shù)值來判斷的原因有如下幾點(diǎn)：首先前文中已經(jīng)提到在比較長時(shí)段的數(shù)據(jù)時(shí)缺少精確度，其次如上文所講因?yàn)槿鄙傧募镜臄?shù)據(jù)所以臭氧的數(shù)據(jù)缺少顯著的差異，所以在進(jìn)行AQI的計(jì)算中時(shí)某些情況下AQI數(shù)值也缺少顯著的差異。

在某些方面該分析存在一定的問題，首先受限于人力與物力本次分析采用的數(shù)據(jù)雖然在時(shí)間上有一定的跨度但是在空間少缺少相應(yīng)的跨度，僅僅考慮了鹽城市的大氣質(zhì)量，本文的結(jié)論在全國范圍內(nèi)能否成立還需要進(jìn)一步的驗(yàn)證。同時(shí)在做最后大氣質(zhì)量的優(yōu)劣時(shí)僅僅對各項(xiàng)數(shù)據(jù)進(jìn)行加和處理，而沒有考慮到具體污染物的比重問題，這在經(jīng)一步的探究時(shí)需要改進(jìn)。

本次疫情對于全國乃至全球都是一個(gè)嚴(yán)峻的挑戰(zhàn)，我們?yōu)榇舜螒?zhàn)疫付出了許多，但是作為一名環(huán)境科學(xué)的研究者來說，我們很高興的看到大氣質(zhì)量得到了改善，雖然這次改善犧牲了部分經(jīng)濟(jì)水平，但是我們發(fā)現(xiàn)了以前困擾我們許久的大氣問題并不再是一味的變壞或者是維持原狀，而是有了提高的可能性，下一步我們要做的就是在經(jīng)濟(jì)和環(huán)境之中找到一個(gè)更好的平衡點(diǎn)。

4：總結(jié)

在疫情期間，由于工廠的停工以及交通運(yùn)輸?shù)牡凸β蔬\(yùn)行，在冬季占主導(dǎo)作用的煤煙型污染與不占主導(dǎo)的交通型污染收到抑制，使大氣質(zhì)量得到有效的提高，同時(shí)對于進(jìn)一步探究人類活動(dòng)（主要為經(jīng)濟(jì)活動(dòng)）與環(huán)境關(guān)系提出了思考。

參考文獻(xiàn)

[1] ?人民日報(bào)http：//society.people.com.cn/n1/2020/0123/c1008-31561292.html

[2] ?百度百科https：//baike.baidu.com/item/AQI/463245？fr=aladdin

[3] ?百度百科https：//baike.baidu.com/item/R語言

[4] ?賈里德·P·蘭德.R語言實(shí)用數(shù)據(jù)分析和可視化技術(shù)[M]208-209

[5] ?浙江大學(xué)，盛驟，謝式千等.概率論與數(shù)理統(tǒng)計(jì)[M]178-182

[6] ?百度百科https：//baike.baidu.com/item/隨機(jī)森林