葛禧,戴馨,黃嬌嬌,吳嘉儀
(1.宿遷學院文理學院,江蘇宿遷 223800;2.浙江海洋大學,浙江舟山 316021)
在全國城鎮化的進程中,環境和生態不可避免地會遭到不同程度的破壞。以煤炭為主的能源消耗大幅攀升,機動車、汽車尾氣等急劇增加,二氧化硫、氮氧化物、臭氧、大氣顆粒物等污染物嚴重影響了空氣質量[1],危害了人們的身體健康。PM2.5通過呼吸進入人體后,容易刺激和腐蝕肺泡壁,損害肺部功能[2]。此外,PM2.5的危害不僅局限于哮喘、呼吸系統驗證、心腦血管疾病,而且還與大腦老化、記憶衰退等多種神經退化性疾病有關[3]。黨的十九大報告提出,“要像對待生命一樣對待生態環境,保護好生態環境就是保護未來”。雖然宿遷市近幾年的空氣質量比北方其他城市更好,但也存在部分月份PM2.5含量較高的現象,所以我們有必要對其進行研究。筆者采用GM(1,1)灰色預測模型對宿遷的PM2.5水平進行預測。
宿遷市逐日空氣PM2.5濃度數據來源于真氣網,該網站是全國公眾空氣質量監測平臺,專注于環境大數據的互聯網化運用,包含PM2.5在內的多種有關數據。真氣網能真實地評估空氣污染水平和空氣質量水平,可靠性較高。筆者截取該網站2013—2022年的數據,通過計算得其月平均值,具體數據如表1所示。

表1 宿遷市2013—2022 年PM2.5 濃度月均值
從表1中我們可以看出,宿遷市的PM2.5平均值較低,并且有著逐年降低的趨勢。對比北京等地區,宿遷市的空氣質量較好。在2013年,北京市有重度污染58天。在2016年,北京市有重度污染39天。而從2013—2021年,北京的PM2.5濃度以每年7μm/m3左右的速度下降。宿遷市近幾年的空氣質量都保持在優良以上,特別是建設全國文明城市以來,宿遷市的空氣質量有著大幅提高。
進行一次累加得到生成時間序列:
式子中X(1)的緊鄰均值生成序列:
建立方程:
白化微分方程為:
最后得出時間響應式為:
1.級比檢驗
建立PM2.5濃度時間序列:
根據公式
求級比:
2.GM(1,1)模型的建立
筆者以表1所示宿遷市2014—2021年的PM2.5濃度年平均值數據為原始數據列,利用灰色預測GM(1,1)對PM2.5濃度進行預測,統計監測點的數據見表1。筆者選用要預測時間之前的1~4組數據來分別建立GM(1,1)模型,并對所得的預測值進行加權平均融合[4]。筆者運用Matlab 2020b編寫相應的代碼,得到預測結果見表2。

表2 PM2.5 的預測結果
由Matlab編寫的代碼可得:
X(1)=(67.916 7,128.750 0,184.500 0,238.416 7,288,335.750 0,381.416 7,417.916 7)
k=0,1,2,…,n-1,…
筆者多次調整程序,使得擬合程度最好,以便預測得出下一年的PM2.5濃度。隨后,筆者將預測出來的年份作為原始數據,重新構建GM(1,1)模型,以此類推,得出宿遷市未來六年的PM2.5濃度。
筆者運用相對誤差檢驗法來將預測值與實際數值進行比較,以此來判斷預測值是否準確。
1.殘差檢驗
2.相對誤差檢驗
3.級比偏差檢驗
筆者選取前四年(2015—2018年)的數據來檢驗模型的誤差。
從表3可以看出,利用該模型所得到的PM2.5濃度預測值與真實值的誤差控制在[-0.02%,3.29%],平均相對誤差為0.017 1。模型擬合精度Q(0)=98.288 4%>90%,說明擬合精度較高。級比偏差P(k)<0.1,也達到了較高要求。

表3 PM2.5 預測的殘差和相對誤差
預測圖如圖1所示(橫軸表示年份,縱軸表示PM2.5濃度)。

圖1 對2022 年的PM2.5 的預測
筆者利用GM(1,1)模型得到未來六年宿遷市PM2.5的濃度情況。由圖1、圖2、圖3、圖4、圖5、圖6可知,PM2.5的濃度總體呈遞減趨勢,預測值與實際值的誤差也較小。由此可見,GM(1,1)模型的預測精度良好。借助GM(1,1)模型,我們可以在生活中有針對性地提出預防措施,從而提高空氣質量。我們借助GM(1,1)模型預測PM2.5的濃度可以減少由空氣污染所引發的疾病。PM2.5的數值變化也可以增強人們的環保意識,如此循環往復,促進環境效益與社會效益的統一。

圖2 對2023 年的PM2.5 的預測

圖3 對2024 年的PM2.5 的預測

圖4 對2025 年的PM2.5 的預測

圖6 對2027 年的PM2.5 的預測
筆者收集了2013—2022年PM2.5的數據,利用GM(1,1)模型和Matlab等數學軟件對宿遷市的PM2.5濃度進行了有效分析,并且預測了宿遷市未來六年PM2.5的濃度情況。