空氣質(zhì)量指數(shù)預(yù)測方法綜述

俞婧婧 唐立力 吳 浩 王志舜 鐘長華

(1.四川輕化工大學(xué)自動化與信息工程學(xué)院,四川 宜賓 644000;2.人工智能四川省重點實驗室,四川 宜賓 644000)

空氣質(zhì)量指數(shù)(AQI)是描述空氣質(zhì)量的一個綜合指數(shù),是衡量空氣質(zhì)量的重要指標(biāo)。該指標(biāo)的計算通常根據(jù)PM2.5、PM10、臭氧等多種污染物濃度進行計算。AQI的取值范圍在0～500,數(shù)值越高,空氣質(zhì)量越差,對人體的危害也越大[1]。AQI是為了讓公眾易于評估空氣質(zhì)量的好壞,并以數(shù)據(jù)的形式協(xié)助環(huán)境管理做出一些重要的決策。1997年6月,我國頒布了國標(biāo)文件《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB3095-1996)采用空氣污染指數(shù)(Air Pollution Index,API)來評價空氣質(zhì)量,API評價體系中僅包含SO2、NO2、PM10這三大污染物。但在實際運用中,我國人民在日常生活中感受到的空氣狀況與API顯示出的空氣狀況存在非常大的反差,為了解決這一問題,我國于2012年頒布了新的國標(biāo)文件《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB3095-2012),用AQI替代API。AQI評價體系在API評價體系的基礎(chǔ)上增加了PM2.5、O3、CO這三種污染物。

隨著全球工業(yè)化與城市化的不斷發(fā)展,空氣污染已成為全球范圍內(nèi)人們廣泛關(guān)注的問題。近年來,在政府的積極推動下,人們的環(huán)保意識有所提高,對降低空氣污染付諸行動,但效果甚微。根據(jù)相關(guān)研究,空氣質(zhì)量的好壞會影響公眾的身體健康,暴露在劣質(zhì)空氣中可能引起過敏反應(yīng),增加公眾患呼吸系統(tǒng)疾病、心血管疾病、神經(jīng)系統(tǒng)等疾病的概率[2-3]。因此,AQI的預(yù)測可以為公眾了解當(dāng)前的空氣質(zhì)量情況,避免長時間暴露在有害氣體中,保護公眾健康。對于保障人類健康和城市空氣污染防治起著至關(guān)重要的作用。

Xu 等[4]選擇了58個氣象因子與多個社會經(jīng)濟因素研究發(fā)現(xiàn)氣象和社會經(jīng)濟因素都與空氣質(zhì)量密切相關(guān),Ye等[5]針對傳統(tǒng)空氣污染研究只關(guān)注單一污染物問題,研究發(fā)現(xiàn)AQI還受GDP、能源消耗以及人口聚集等因素影響。以上研究說明了影響AQI的因素眾多。倘若單從氣象環(huán)境復(fù)雜多變這一特點來說,AQI預(yù)測的準(zhǔn)確度難以保證。并且對于隨機序列來說,找出自變量與因變量之間的關(guān)系關(guān)山難越。因此近些年AQI預(yù)測成為了許多學(xué)科較受關(guān)注的研究課題。目前,針對AQI預(yù)測的方法主要有統(tǒng)計預(yù)測方法[6]、機器學(xué)習(xí)的預(yù)測方法[7]以及組合模型的預(yù)測方法[8]。

1 統(tǒng)計預(yù)測方法

統(tǒng)計預(yù)測方法在AQI預(yù)測中的應(yīng)用可以追溯到20世紀(jì)60年代,當(dāng)時美國政府開始研究如何減少城市空氣污染。為了更好地監(jiān)測和預(yù)測空氣污染情況,政府開始建立空氣質(zhì)量監(jiān)測站,收集氣象數(shù)據(jù)和AQI等信息。隨著數(shù)據(jù)收集的不斷增加,政府開始嘗試使用統(tǒng)計預(yù)測方法來預(yù)測AQI。

統(tǒng)計預(yù)測方法需在建模前對歷史數(shù)據(jù)進行預(yù)處理,然后在處理好的數(shù)據(jù)中選擇部分訓(xùn)練數(shù)據(jù)構(gòu)建數(shù)學(xué)模型,最后使用構(gòu)建好的數(shù)學(xué)模型預(yù)測未來值。在預(yù)測的過程中,需要對模型不斷地迭代更新,使其能夠適應(yīng)數(shù)據(jù)與環(huán)境的變換。統(tǒng)計預(yù)測方法是一種常用的AQI預(yù)測方法。目前,常見的統(tǒng)計預(yù)測方法包括多元線性回歸模型(MLR)、時間序列模型、灰色模型(GM)等。

MLR是根據(jù)實際問題確定自變量與因變量,并利用樣本建立多元線性回歸模型。其因可以處理多個輸入變量的優(yōu)點而被廣泛使用,但該模型的穩(wěn)定性受輸入變量之間相關(guān)性影響,并且過多或過少的自變量都會導(dǎo)致模型預(yù)測能力下降。因此該模型通常用于模型比對,如Lei等[9]通過實驗將MLR與機器學(xué)習(xí)的預(yù)測方法進行比較后發(fā)現(xiàn),MLR的預(yù)測性能不及機器學(xué)習(xí)預(yù)測方法;石曉文等[10]為了證實長短期記憶網(wǎng)絡(luò)(LSTM)的預(yù)測精度,將LSTM與支持向量機(SVM)、MLR、反向傳播神經(jīng)網(wǎng)絡(luò)(BP)進行實驗比對,得出LSTM在AQI≤200時預(yù)測精度高。

時間序列模型包括自回歸模型(AR)、移動平均模型(MA)、自回歸移動平均模型(ARMA)以及差分自回歸移動模型(ARIMA)等預(yù)測模型,這些模型都是用于時間序列數(shù)據(jù)預(yù)測。ARIMA是基于ARMA模型衍生出來的,ARMA由AR與MA模型組成。其中AR描述了當(dāng)前值與前幾個值之間的關(guān)系,MA描述了當(dāng)前值與前幾個隨機誤差之間的關(guān)系。ARMA通常包含兩個參數(shù):自回歸項的數(shù)量p與移動平均項的數(shù)量q,在建模過程中這兩個參數(shù)需要根據(jù)歷史數(shù)據(jù)進行選擇。因ARMA對非平穩(wěn)時間序列的擬合效果較差,故ARIMA在ARMA模型的基礎(chǔ)上引入了差分的概念,使其可以處理非平穩(wěn)序列?；舅枷胧窍葘υ紩r間序列進行階數(shù)為d的差分將信號與噪聲分離,此時非平穩(wěn)序列變?yōu)榱似椒€(wěn)序列,然后再把無噪聲的信號作為數(shù)據(jù)源應(yīng)用ARMA模型來獲得預(yù)測值。Slini等[11]使用ARIMA模型預(yù)測AQI,但該模型的缺點是對于高階差分的計算較為復(fù)雜、僅挖掘時間序列自身的數(shù)據(jù)特征,在實際應(yīng)用中較為簡單。孫朝云等[12]通過對比多個模型(含MLR與ARIMA)的決定系數(shù)(R)值得出LSTM、門控循環(huán)單元(GRU)和雙向長短期記憶網(wǎng)絡(luò)(Bi-LSTM) 性能較好。

灰色模型主要包括GM(1,1)模型和GM(2,1)模型。GM(1,1)模型是通過對原始數(shù)據(jù)進行累加后得到累加數(shù)列,然后通過建立一階線性微分方程對累加數(shù)列進行預(yù)測。GM(2,1)模型在GM(1,1)模型的基礎(chǔ)上增加了一階差分序列,以便獲得更為準(zhǔn)確的預(yù)測結(jié)果?；疑Ｐ蛯π颖緮?shù)據(jù)集具有較好的預(yù)測效果,計算速度快、精度高且不受噪聲影響。但是,當(dāng)數(shù)據(jù)呈周期性、季節(jié)性時,該模型可能出現(xiàn)較大誤差,需要進行適當(dāng)?shù)念A(yù)處理。并且,灰色模型的預(yù)測精度受參數(shù)影響較大,即若參數(shù)選擇不當(dāng),模型的預(yù)測能力將會受到影響。因此有學(xué)者在灰色模型的基礎(chǔ)上進行了優(yōu)化研究,如Xiong 等[13]在灰色模型的基礎(chǔ)上進行了擴展,改進了參數(shù)估計法,獲得了令人滿意的預(yù)測結(jié)果。此外,經(jīng)研究發(fā)現(xiàn)數(shù)據(jù)質(zhì)量也會直接影響該模型的預(yù)測結(jié)果,且因其依賴歷史數(shù)據(jù),導(dǎo)致其在中長期預(yù)測中的效果不佳[14]。

總的來說,統(tǒng)計預(yù)測方法是基于統(tǒng)計數(shù)據(jù)和歷史數(shù)據(jù)來構(gòu)建數(shù)學(xué)模型,利用構(gòu)建好的模型預(yù)測未來的AQI值。統(tǒng)計預(yù)測模型雖然可以高效地分析并處理歷史數(shù)據(jù),以便精確地預(yù)測未來值,但該方法的局限性在于它不能很好地處理非線性關(guān)系和復(fù)雜的影響因素。并且,統(tǒng)計預(yù)測方法一般不考慮空氣污染物濃度,而是只通過分析空氣質(zhì)量的歷史數(shù)據(jù)來預(yù)測AQI,僅根據(jù)局部數(shù)據(jù)規(guī)律建立模型,并不適用于中長期預(yù)測。

2 機器學(xué)習(xí)預(yù)測方法

機器學(xué)習(xí)預(yù)測方法一般需要較大的數(shù)據(jù)量,因其擁有著強大的計算能力,且在不同領(lǐng)域應(yīng)用廣泛。目前,BP、SVM、隨機森林(RF)、極限學(xué)習(xí)機(ELM)、LSTM等機器學(xué)習(xí)預(yù)測模型在AQI預(yù)測研究中備受學(xué)者青睞。

BP神經(jīng)網(wǎng)絡(luò)通過反向傳播算法來訓(xùn)練模型,不斷調(diào)整網(wǎng)絡(luò)中的權(quán)重和偏置,以適應(yīng)新的數(shù)據(jù)從而使預(yù)測結(jié)果更加準(zhǔn)確[15]。該網(wǎng)絡(luò)可以模擬非線性數(shù)據(jù),適合用于復(fù)雜的預(yù)測問題。閆妍等[16]使用BP神經(jīng)網(wǎng)絡(luò)預(yù)測AQI,得出該模型雖然速度慢,但適應(yīng)性較好,得到的結(jié)果比對還算成功。此外,單一的BP神經(jīng)網(wǎng)絡(luò)在搜索較大解空間時,容易陷入局部最優(yōu),影響預(yù)測結(jié)果。

SVM是一種常用的分類回歸算法,對于小樣本的數(shù)據(jù)集具有較好的泛化能力。此外,SVM具有較好的非線性建模能力。其核心思想是通過計算找到一個最佳分離超平面對新數(shù)據(jù)進行分類[17],當(dāng)數(shù)據(jù)不是線性可分時,通過核函數(shù)將數(shù)據(jù)映射到高維的空間中進行處理。在實際應(yīng)用中,SVM需根據(jù)具體情況選擇合適的核函數(shù)與調(diào)節(jié)參數(shù)來優(yōu)化分類效果。王春梅等[18]將BP與SVM用于預(yù)測西安的AQI,通過實驗對比得到BP與SVM在預(yù)測上雖有較大誤差,但預(yù)測結(jié)果與實際值的總體走勢相差不多,并且BP的均方誤差值略大于SVM的均方誤差值。但當(dāng)數(shù)據(jù)集規(guī)模較大時SVM訓(xùn)練速度慢,且其對噪聲及異常值比較敏感,在訓(xùn)練前需進行預(yù)處理。

RF是決策樹(DTs)與裝袋(bagging)算法的結(jié)合,由多個相互獨立的RF組成,各個RF之間隨機抽取特征進行訓(xùn)練,最終綜合各個RF的輸出結(jié)果進行預(yù)測。楊思琪等[19]研究運用了RF及回歸算法構(gòu)建了AQI預(yù)測模型,并進行了分類模型的對比。該研究表明RF泛化性與穩(wěn)健性較好,且在預(yù)測AQI時有較高的預(yù)測精度。但當(dāng)數(shù)據(jù)集中噪聲較大時,使用RF會導(dǎo)致模型產(chǎn)生過擬合現(xiàn)象。

ELM輸入層接收數(shù)據(jù)特征,隱藏層中每個神經(jīng)元隨機生成權(quán)重,采用sigmoid等激活函數(shù)將輸入特征映射到到非線性空間,最后將結(jié)果傳輸至輸出層進行線性回歸得到預(yù)測結(jié)果[20]。該模型訓(xùn)練速度快并且可避免局部最優(yōu)問題[21]。李澤婧等[22]所做實驗中ELM預(yù)測AQI的MSE達到了6.923。但由于ELM的隱藏層節(jié)點數(shù)、權(quán)值和閾值是隨機生成的,從而出現(xiàn)ELM學(xué)習(xí)能力下降的問題[23]。

LSTM是一種循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN),該模型通過一系列LSTM單元處理輸入數(shù)據(jù),每個LSTM單元包括輸入門、遺忘門以及輸出門。其中,輸入門決定當(dāng)前時間步輸入的數(shù)據(jù)對當(dāng)前狀態(tài)的影響,遺忘門決定了上一個時間步的狀態(tài)對當(dāng)前狀態(tài)的影響,輸出門控制當(dāng)前狀態(tài)的輸出。LSTM從時間序列數(shù)據(jù)中學(xué)習(xí)長期依賴關(guān)系,有效地捕獲長序列之間的時間相關(guān)性。張春露等[24]利用LSTM模型預(yù)測空氣質(zhì)量指數(shù),雖然其預(yù)測的精準(zhǔn)度不錯。但將LSTM應(yīng)用到AQI預(yù)測中存在訓(xùn)練時間長、容易過擬合、受限于數(shù)據(jù)質(zhì)量以及需要大量數(shù)據(jù)等缺點。

綜上可知,雖然機器學(xué)習(xí)在AQI預(yù)測研究中廣泛應(yīng)用,但目前對于機器學(xué)習(xí)算法與空氣質(zhì)量問題結(jié)合所需關(guān)注和解決的問題主要有:(1)如何在保證模型精度的同時減少計算復(fù)雜度;(2)如何避免過擬合;(3)如何避免局部最優(yōu);(4)如何選擇合理算法以提高預(yù)測準(zhǔn)確性和預(yù)測精度保持較高水平;(5)如何提高數(shù)據(jù)質(zhì)量等。

3 組合預(yù)測方法

組合預(yù)測方法是指將兩個及兩個以上的模型通過取長補短組合在一起后使預(yù)測效果提高。組合預(yù)測模型不僅具有適應(yīng)未來預(yù)測環(huán)境變化的能力,還彌補了單一模型的預(yù)測精度與穩(wěn)定性差的問題[25]。目前,組合預(yù)測方法主要包含基于模型融合的組合預(yù)測方法、基于數(shù)據(jù)融合的組合預(yù)測方法以及基于參數(shù)優(yōu)化的組合預(yù)測方法。

基于模型融合的組合預(yù)測方法是指將多個不同的預(yù)測模型組合進行預(yù)測,基本思想是綜合多個模型的預(yù)測結(jié)果,以不同權(quán)重進行組合來提高預(yù)測準(zhǔn)確性和魯棒性,例如, Song等[26]為了提高AQI預(yù)測的準(zhǔn)確性,將徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(RBFNN)、集合經(jīng)驗?zāi)B(tài)分解(EEMD)以及ARIMA三種模型通過權(quán)重分配進行組合來預(yù)測AQI,結(jié)果表明組合模型可以提高AQI預(yù)測的精度,有效減少了預(yù)測誤差。

基于數(shù)據(jù)融合的組合預(yù)測方法是通過將多源數(shù)據(jù)進行融合預(yù)測,正如前面所說,影響AQI的因素有很多,因此需要輸入數(shù)據(jù)中納入對其影響大的影響因子。例如,楊豐玉等[27]在預(yù)測AQI時將風(fēng)力這一影響因子融合進來,實驗表明,融合風(fēng)力影響因子與傳統(tǒng)K鄰近(KNN)相比AQI的預(yù)測準(zhǔn)確率顯著提升。

基于參數(shù)優(yōu)化的組合預(yù)測方法是通過優(yōu)化算法來選取模型中的最優(yōu)超參數(shù),例如,為了解決BP收斂速度慢,容易陷入局部極小值的問題,Xu等[28]通過對慣性權(quán)重的優(yōu)化來改善粒子群優(yōu)化算法(PSO)的探索能力,再將優(yōu)化后PSO與BP組合,以優(yōu)化BP的權(quán)值和閾值,最終通過對比實驗表明作者提出的改進PSO-BP算法具有較高的預(yù)測精度。

綜上可知,不同的AQI預(yù)測方法都有各自的獨特之處。組合預(yù)測方法的優(yōu)勢在于其靈活性,它通過將多個不同的預(yù)測模型、數(shù)據(jù)源和特征選擇算法結(jié)合起來,提高AQI預(yù)測的準(zhǔn)確性和魯棒性。為了更好地比較不同方法的特點,本文提供以下表格,總結(jié)了統(tǒng)計預(yù)測方法、機器學(xué)習(xí)預(yù)測方法和組合預(yù)測方法之間的比較,如表1所示。

表1 不同方法比較

4 展望與總結(jié)

近些年,隨著環(huán)境科學(xué)領(lǐng)域?qū)諝赓|(zhì)量問題的日益關(guān)注,AQI預(yù)測成為了該領(lǐng)域的一個重要研究焦點。本文采用綜述的方式,系統(tǒng)闡述并分析了AQI預(yù)測方法的演進歷程以及不同預(yù)測方法的特點和適用場景。盡管目前AQI預(yù)測研究已經(jīng)取得了一定的研究成果,但依然存在以下挑戰(zhàn)和問題:

(1)復(fù)雜的影響因素?？諝赓|(zhì)量波動往往受到多種復(fù)雜因素的共同影響,包括氣象條件、人類活動、地理位置等因素。舉例來說,氣象因素如溫度、濕度、風(fēng)速等會對顆粒物的濃度和擴散產(chǎn)生影響,溫度和濕度的變化可能會導(dǎo)致顆粒物的濃度波動,風(fēng)速則會影響顆粒物的擴散和輸送[29],這些效應(yīng)常常呈現(xiàn)為非線性關(guān)系,使得建立精準(zhǔn)的AQI預(yù)測模型變得更具挑戰(zhàn)。因此在進行AQI預(yù)測時必須考慮這些因素對AQI的直接和間接影響,研究人員須深入挖掘這些復(fù)雜的空間關(guān)聯(lián),特別是考慮到城市內(nèi)不同區(qū)域之間的空氣質(zhì)量差異。

(2)數(shù)據(jù)質(zhì)量問題。在一些發(fā)展中國家和地區(qū),AQI監(jiān)測系統(tǒng)不完善,監(jiān)測站點不足,會導(dǎo)致數(shù)據(jù)缺失、噪聲干擾嚴(yán)重等問題,這會影響預(yù)測模型的準(zhǔn)確性和魯棒性。因此,除了需要加強監(jiān)測設(shè)備的更新和維護、提高數(shù)據(jù)質(zhì)量之外。可以考慮采用合適的數(shù)據(jù)處理方法和技術(shù),以減少數(shù)據(jù)缺失和噪聲干擾的影響。舉例來說,采用數(shù)據(jù)插補的方法填充數(shù)據(jù)缺失的部分,以維持?jǐn)?shù)據(jù)的完整性,例如,線性插值、均值替換、回歸替換等方法根據(jù)歷史數(shù)據(jù)來估算缺失的數(shù)值以確保檢測數(shù)據(jù)的連續(xù)性;采用噪聲濾波技術(shù)減少噪聲對數(shù)據(jù)的干擾,例如,使用平均濾波、中值濾波、小波變換等技術(shù)平滑數(shù)據(jù),去除異常值,以提高數(shù)據(jù)質(zhì)量。

(3)時空特性建模挑戰(zhàn)?？諝赓|(zhì)量受季節(jié)、時間和地理位置等因素的影響,在時空維度上存在明顯的變化?？紤]到時間因素,過去的AQI數(shù)據(jù)對未來的空氣質(zhì)量具有直接影響,因此,需要建立模型來捕捉歷史數(shù)據(jù)對AQI預(yù)測的影響。這可能需要動態(tài)更新權(quán)重以反映不同時間點的重要性,以確保模型能夠充分考慮歷史數(shù)據(jù)的影響。此外,空間依賴性在AQI預(yù)測中也至關(guān)重要,因為不同地理位置之間的氣象條件和人類活動水平差異較大,直接影響空氣質(zhì)量。因此,建立空間關(guān)聯(lián)模型來考慮不同區(qū)域之間的空氣質(zhì)量差異是一項重要挑戰(zhàn)?？臻g插值、地理信息系統(tǒng)(GIS)等方法可以用于建立這種空間依賴性的模型,從而更準(zhǔn)確地預(yù)測不同地區(qū)的AQI。

總的來說,解決這些挑戰(zhàn)需要繼續(xù)深入研究和發(fā)展AQI預(yù)測模型,結(jié)合先進的數(shù)據(jù)處理技術(shù)、時空建模方法以及監(jiān)測設(shè)備的改進,以提高AQI預(yù)測的準(zhǔn)確性和可靠性。這將有助于更有效地管理和改善空氣質(zhì)量,保護環(huán)境和人類健康。