淮南市PM2.5濃度變化分析及預測

來鵬飛

摘 要：本研究對淮南市2015—2021年PM2.5濃度進行了時間序列分析和算法預測。結果表明：PM2.5濃度最高值出現在每年12月到次年1月，最低值出現在每年7、8月。濃度值在2020年之前處于上升趨勢，2020年開始有所下降。為驗證預測結果的準確性，使用了CNN-GRU預測模型進行預測，預測結果的MSE、RMSE、MAE、MAPE和SMAPE分別為122.07、11.05、9.67、24.64%和21.18%。預測結果表明，PM2.5濃度的實際值變化趨勢有所下降，可為日后環境治理提供參考。

關鍵詞：淮南市;PM2.5預測;CNN-GRU;時間序列

中圖分類號：X513 ? ? 文獻標志碼：A ? ? 文章編號：1003-5168（2022）12-0128-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2022.12.027

Analysis and Prediction of PM2.5 Concentration Change in Huainan City

LAI Pengfei

（School of Electrical and Information Engineering， Anhui University of Technology， Huainan 232001，China）

Abstract： The time series analysis and algorithm prediction of PM2.5 concentration in Huainan city from 2015 to 2021 were carried out. The results show that the highest value of PM2.5 concentration occurs from December to January， and the lowest value occurs in July and August. The concentration was in an upward trend before 2020， and began to decline in 2020. In order to verify the accuracy of prediction results， CNN-GRU prediction model is used for prediction. The predicted MSE， RMSE， MAE， MAPE and SMAPE were 122.07， 11.05， 9.67， 24.64% and 21.18%， respectively. The prediction results show that the actual value of PM2.5 concentration decreases， which can provide reference for future environmental governance.

Keywords： Huainan City; PM2.5 prediction; CNN-GRU; time series

0 引言

霧霾天氣在我國多個地區頻繁出現，導致環境問題日益加重。淮南市作為我國的主要煤生產基地之一，在煤礦開采的過程中，由于運輸燃燒等方式會導致許多煙塵釋放到空氣中，其中包含了大量的細顆粒物（PM2.5）[1]。本研究采用2015—2021年的數據研究PM2.5的特征及變化，了解PM2.5的污染狀況，為淮南市大氣污染防控提供數據參考。

PM2.5是指環境空氣中空氣動力學直徑小于等于2.5 μm、大于0.1 μm的顆粒物。PM2.5對能見度以及空氣質量有著很大的影響，并且PM2.5中含有的重金屬、微生物等有害物質，也對人類的健康有著嚴重的危害，可能會引發呼吸疾病[2-4]，同時對氣候的變化也有著很大的影響[5-7]。因此，對PM2.5濃度變化進行分析以及準確預測對環境的治理和人類的健康有著重要的意義。

以往研究多是基于多元時間序列模型或多元線性回歸的方式實現對PM2.5濃度的預測。隨著人工神經網絡的發展，出現了如卷積神經網絡（CNN）[8-10]和循環神經網絡（RNN）[10-11]這一類的深度神經網絡。CNN的優勢主要體現在圖像處理中，在時間序列上的預測效果并不明顯，RNN本身是處理時間序列的模型，可以將多種因素和時間序列考慮在內，但只能對短序列進行建模。本研究采用CNN-GRU[12]網絡對PM2.5濃度進行預測和趨勢分析，并與GRU[13-14]、LSTM[15-16]和ARIMA[17]的結果進行比較分析，以期為PM2.5污染治理提供依據。

1 數據來源與方法

1.1 數據來源

本研究使用2016年1月至2021年12月的每日數據為研究基礎。數據源自美國國家氣候數據中心（NCDC），隸屬于美國國家海洋及大氣管理局（NOAA）。提供自1942年以來的中國地面氣象數據下載。下載的數據格式為ISD-Lite，是一種簡化了的ISD（Integrated Surface Data）數據。每個文件按天保存，內容包含時間間隔為1 h的單日AQI、PM2.5、PM10、SO2、NO2、O3以及CO的值。

1.2. 數據處理方法

首先通過python代碼對按日排列的數據進行整合，使其從按日存儲的表格變為按年存儲的表格，然后通過python與MySQL數據庫建立連接，將其存儲到MySQL數據庫中，按照需求將淮南市每日12 h的月均值PM2.5數據查詢出來進行數據處理分析。根據最終結果繪制月均濃度趨勢圖并進行預測分析，以便研究淮南市PM2.5的變化。

2 結果與分析

2.1 PM2.5濃度的變化

如圖1所示，近年來PM2.5濃度呈周期性變化，每年PM2.5濃度呈現出“V”字形變化。最高值出現在每年12月到次年1月，濃度最高出現在2018年1月，為108.65 μg/m3。最低值出現在每年7、8月，濃度最低出現在2021年7月，為16.29 μg/m3。淮南市PM2.5的濃度主要分布在35～75 μg/m3這個范圍內。從圖1中可以看出，2020年之前PM2.5月均值濃度一直處于上升的趨勢，2020年開始有明顯的下降趨勢。

如表1所示，將一天24 h里檢測到的PM2.5濃度的值求平均值，按照空氣質量標準，將空氣質量等級共分為6個等級，分別為優、良、輕度污染、中度污染、重度污染以及嚴重污染。從表1中可以看出，近年來淮南市空氣質量等級主要處于優、良以及輕度污染三個水平，并集中在良這個等級。2020—2021年出現輕度污染的月份每年僅有1個月，而空氣質量為優的月份增加到了每年5個月，這說明了淮南市的空氣質量在近年有所改善。

2.2 時間序列預測分析

2.2.1 卷積神經網絡。卷積神經網絡（Convolutional Neural Networks，CNN）是深度學習神經網絡模型中的一種。在圖像處理、視頻處理等方面有著廣泛的應用，并表現出優異的性能，實現對輸入特征的提取。CNN一般有卷積層、池化層和全連接層組成。卷積層是CNN網絡的核心，通過對輸入的數據進行卷積運算，降低噪聲，并增強輸入數據的特征;池化層的功能是將卷積層的輸出作為輸入，通過最大值、平均值池化等操作，保留特征的同時，減少參數實現對數據的降維;全連接層是將經過卷積層、池化層處理過的數據進行整合，得到分類或者回歸的結果

2.2.2 GRU網絡。GRU網絡是在LSTM網絡的基礎上進行改進的，組成有更新門以及重置門。更新門的作用是判斷目前時間的輸入與前面的隱藏狀態是否更新到候選隱藏狀態。重置門是將目前時間的輸入與前面的隱藏狀態相結合，來判斷是否將候選隱藏狀態中的信息遺忘。更新門更適用于中長期序列中，而重置門更適用于短期序列中。更新門和重置門的計算公式如式（1）、式（2），候選隱藏狀態的計算公式如公式（3）所示。

式中：zt、rt、xt、ht分別是更新門的輸入、重置門的輸入、當前時間狀態的輸入以及上一時間狀態的隱藏狀態;xt是t時刻的輸入;ht-1是t-1時刻的隱藏狀態;Wz、Uz、Wr、Ur是權重矩陣;br、bz是偏差矩陣;σ（·）表示sigmoid函數，它將得到的激活結果變換到0與1之間。

最后利用網絡計算出當前時刻的最終狀態ht，公式如式（4）。

2.2.3 CNN-GRU。本研究使用的是一種融合CNN和GRU的神經網絡混合模型，可以更好地學習數據中的特征和規律，準確地預測出結果。模型結構如圖2所示。輸入數據，經過卷積層的計算來得出輸入數據中的特征信息;池化層中采用最大池化的方法進行池化，保留最主要的特征;Dropout層隨機丟棄一些神經元來防止出現過擬合現象;在GRU層中，通過不斷地訓練調整自身的參數，不斷優化并進行預測，最終在全連接層中對輸出進行處理得到預測的結果。

在進行預測時，首先對原始數據進行處理，對淮南PM2.5時間序列數據的數據特征進行歸一化處理，并將數據映射到區間[0，1];然后將數據集按照6∶1分為兩個部分：訓練集和測試集，當訓練結束后對預測出的結果進行反歸一化;最后將2021年數據作為測試集與預測結果進行對比。

為驗證預測結果的可靠性，另外使用了GRU、LSTM、ARIMA三種算法進行預測對比。圖3展示了GRU、LSTM和ARIMA算法對未來1年時間的預測結果與原始數據的對比情況。結果表明，三種算法預測的趨勢基本一致，都呈現出“V”字形的變化。但GRU算法預測結果只有3月和6月的濃度比真實值低，其他月份都略高于真實值;LSTM預測結果全部高于真實值;ARIMA預測3月、8月、10月、11月、12月結果比真實值低，其余月份高于真實值。這說明了2021年PM2.5濃度值相較于之前的變化有所降低。

表2為CNN-GRU、GRU、LSTM和ARIMA預測結果的評價指標比較。分別選取了均方誤差（MSE）、均方根誤差（RMSE）、平均絕對誤差（MAE）、平均絕對百分比誤差（MAPE）和對稱平均絕對百分比誤差（SMAPE）作為評判標準。從表1中可以清晰地看出，四種算法預測的結果中，CNN-GRU預測的結果最為準確，分別是119.53、10.93、8.77、24.42%和20.46%。

3 結語

本研究分析了2015—2021年淮南市PM2.5的變化趨勢并使用了相應的算法對數據進行預測分析。結果表明，2020年之前淮南市PM2.5濃度具有緩慢上升的趨勢。其中最高值出現在每年12月與次年的1月，最低值出現在每年7、8月;空氣質量等級主要處于優、良以及輕度污染三個水平，主要集中在良的水平。2020年開始，淮南市PM2.5濃度開始下降，出現輕度污染的月份有所降低。使用CNN-GRU算法與其他算法預測結果比較得出，CNN-GRU算法相較于其他算法預測準確性更高。預測結果顯示，預測值略高于真實值，表明淮南市近1年PM2.5濃度有所降低。

參考文獻：

[1] 胡煜，鄭劉根，程樺，等.淮南市PM2.5中PAHs污染特征及來源分析[J].環境監測管理與技術，2016，28（6）：33-37.

[2] 王瑋，湯大鋼，劉紅杰，等.中國PM2.5污染狀況和污染特征的研究[J].環境科學研究，2000（1）：1-5.

[3] BELL M L，FRANCESCA D，KEITA E，et al. Spatial and temporal variation in PM2.5 chemical composition in the United States for health effects studies [J]. Environmental health perspectives， 2007，115（7）：989-995.

[4] 王庚辰，王普才.中國PM2.5污染現狀及其對人體健康的危害[J].科技導報，2014，32（26）：72-78.

[5] TAI A P K，MICKLEY L J，JACOB D J. Correlations between fine particulate matter （PM 2.5 ） and meteorological variables in the United States： Implications for the sensitivity of PM 2.5 to climate change [J]. Atmospheric Environment， 2010，44（32）：

[6] 王薇，陳明.城市綠地空氣負離子和PM_（2.5）濃度分布特征及其與微氣候關系：以合肥天鵝湖為例[J].生態環境學報，2016，25（9）：1499-1507.

[7] 陳波，魯紹偉，李少寧.北京城市森林不同天氣狀況下PM2.5濃度變化[J].生態學報，2016，36（5）：1391-1399.

[8] 陸繼翔，張琪培，楊志宏，等.基于CNN-LSTM混合神經網絡模型的短期負荷預測方法[J].電力系統自動化，2019，43（8）：131-137.

[9] 羅文慧，董寶田，王澤勝.基于CNN-SVR混合深度學習模型的短時交通流預測[J].交通運輸系統工程與信息，2017，17（5）：68-74.

[10] 范竣翔，李琦，朱亞杰，等.基于RNN的空氣污染時空預報模型研究[J].測繪科學，2017，42（7）：76-83，120.

[11] 李潔，林永峰.基于多時間尺度RNN的時序數據預測[J].計算機應用與軟件，2018，35（7）：33-37，62.

[12] 姚程文，楊蘋，劉澤健.基于CNN-GRU混合神經網絡的負荷預測方法[J].電網技術，2020，44（9）：3416-3424.

[13] 王增平，趙兵，紀維佳，等.基于GRU-NN模型的短期負荷預測方法[J].電力系統自動化，2019，43（5）：53-58.

[14] 趙兵，王增平，紀維佳，等.基于注意力機制的CNN-GRU短期電力負荷預測方法[J].電網技術，2019，43（12）：4370-4376.

[15] 王鑫，吳際，劉超，等.基于LSTM循環神經網絡的故障時間序列預測[J].北京航空航天大學學報，2018，44（4）：772-784.

[16] 白盛楠，申曉留.基于LSTM循環神經網絡的PM2.5預測[J].計算機應用與軟件，2019，36（1）：67-70，104.

[17] 吳玉霞，溫欣.基于ARIMA模型的短期股票價格預測[J].統計與決策，2016（23）：83-86.