付瑩瑩 宋李俊 唐永晟 張月月
重慶理工大學 重慶 400054
由于極端天氣的出現(xiàn)越漸頻繁,全球變暖越來越引起大家的關注。對于全球變暖有很多種解釋,其中最主要的還是溫室氣體排放和溫室效應不斷積累所致[1]溫室效應的產(chǎn)生則是由于問題氣體能夠很好的吸收地球產(chǎn)生的長波輻射,而對太陽輻射不具有吸收性,使得大氣與地面間的溫度不斷升高而產(chǎn)生的。文獻指出,溫室氣體影響了地氣系統(tǒng)的能量發(fā)射和吸收平衡[2],導致全球變暖。
從氣候角度研究全球溫度變化需要全球范圍長時間的觀測積累,但是由于過去這方面的時空數(shù)據(jù)并不完整,給統(tǒng)計計算帶來極大困難。不僅如此,海洋吸收熱量對全球氣候變化的影響很大[5][6],如年代際太平洋震蕩[3]、厄爾尼諾現(xiàn)象[4]、拉尼娜現(xiàn)象[5]等。這些因素使得研究全球溫度變化更加困難。
目前對于氣候預測的研究已經(jīng)很廣泛了,對于用來預測天氣的方法也更是多種多樣。其中最主要的還是以線性回歸函數(shù)和均值回復的方法進行建模[7]。MasakoIkefuji等根據(jù)氣候變化的動態(tài)提出了一種基于回歸數(shù)值方法來求解的,具有潛在的重尾不確定性和通用效用函數(shù)的動態(tài)有限水平經(jīng)濟-氣候模型,但該模型具有一定的不確定性,并且研究者的偏好對其影響較大[8];Wakjira T.Dibaba等[9]發(fā)現(xiàn)降水比溫度更好地再現(xiàn)了季節(jié)異常的年際變化,同時,發(fā)現(xiàn)使用多個模型組合分析降雨量對氣候的影響的效果更好;J.HN Palma等[10]將模擬的森林生長與觀測和模擬的氣候數(shù)據(jù)進行比較發(fā)現(xiàn)輻射和溫度的變化對于氣候的影響是巨大的。
因此本文對比各種因素對氣候變化的影響,考慮地球輻射(長波輻射和短波輻射)、降雨量和全球云量變化對全球氣候多尺度時間變化的影響。利用小波分析對氣象要素的歷史數(shù)據(jù)進行處理,然后結(jié)合NAR人工神經(jīng)網(wǎng)絡對全球氣溫變化趨勢進行預測,以了解全球氣溫變化規(guī)律及氣候突變性,更好地服務于短期氣候預測。
2.1 小波分析 小波分析是一種對函數(shù)或信號的時間軸和頻率軸上不同尺度和不同周期的演變進行多尺度細化分析的方法,這種方法可以很好的補償傅立葉變換的缺陷。氣候變化總是具有多個時間尺度,因此,小波變換非常適合用于氣候變化的研究,并已在大氣科學領域中取得了一系列成果,是研究大氣中多尺度變換的有力工具。本文利用小波分析的方法,對全球1994-2001年的氣候變化數(shù)據(jù)進行分析。
本文小波分析采用Mallat算法,將要素的歷史氣候數(shù)據(jù)通過一定尺度分解近似部分和細節(jié)分量;然后通過對分解后的部分進行平滑處理和重構(gòu);最后用小波強制消噪法去除分解后的數(shù)據(jù)中的噪聲,僅保留余下的重構(gòu)值。其數(shù)學表達式為:
逼近系數(shù)和小波系數(shù)[11]可以通過公式求得:
其中:α為伸縮尺度因子,決定小波的寬度;b為平移因子,反映小波位移量,其中α,b∈R且α≠0,α可以任意取值。根據(jù)文獻[12],α的取值為1-96。
2.2 人工神經(jīng)網(wǎng)絡 人工神經(jīng)網(wǎng)絡近年來在全球氣候研究和天氣預測中應用廣泛,特別是對于非線性問題的解決優(yōu)于其他傳統(tǒng)的方法。在進行小波分析之后,采用NAR人工神經(jīng)網(wǎng)絡對歷史氣候數(shù)據(jù)的時間序列進行預測。其原理為:把消噪后的歷史氣候數(shù)據(jù)作為輸入項輸入到輸入層,通過隱藏神經(jīng)元的加權(quán)形成傳遞函數(shù)輸入;網(wǎng)絡輸出歷史氣候數(shù)據(jù)的擬合值,并將輸出的擬合值通過反饋連接反饋到前面輸入層形成遞歸調(diào)用[11];同時將誤差和擬合值反饋到網(wǎng)絡中進行調(diào)整,得到最后神經(jīng)網(wǎng)絡。滯后結(jié)構(gòu)采用再次調(diào)用歷史氣候數(shù)據(jù)為輸入,求出預定步長內(nèi)的氣候預測值。則通過氣候要素預測時間序列的數(shù)學模型的表達式如下:
式中:Ft為未來的預測值;Ft-u為氣候要素u階滯后變量;p是模型的階數(shù);l為步長;q是隱藏神經(jīng)元個數(shù);αio是隱藏神經(jīng)元的偏差權(quán)值;αio為u階滯后的變量連接輸出的神經(jīng)元的權(quán)值;bo神經(jīng) 元輸入時的偏差權(quán)值;bi是隱藏層與輸出層神經(jīng)元之間的連接權(quán)值。
本文采用db4小波對全球1994年1月至2001年4月的太陽輻射、云層數(shù)據(jù)和降雨量分解為3個層次的多尺度小波分析,以得到數(shù)據(jù)的對比如圖2.1。從圖中可以看出,去噪后的數(shù)據(jù)(藍色曲線)完整的保留了原始數(shù)據(jù)的數(shù)據(jù)特征,僅僅去除了噪聲數(shù)據(jù),對后續(xù)的分析不會產(chǎn)生影響。
圖2.1 小波去噪結(jié)果對比圖
表2.1 數(shù)據(jù)集分類
然后在matlab上進行NAR人工神經(jīng)網(wǎng)絡模型的計算。具體求解如下:
首先,對NAR神經(jīng)網(wǎng)絡進行輸入訓練。即將輸入的數(shù) 據(jù)按上表自動分成training set(訓練集)、validation set(檢驗集)及test set(測試集)三部分,其中訓練集為訓練樣本數(shù)據(jù);檢驗集為驗證樣本數(shù)據(jù),測試集為測試樣本數(shù)據(jù),三個數(shù)據(jù)集是不重疊的。在訓練時,用訓練集訓練,每訓練一次,系統(tǒng)自動會將檢驗集中的樣本數(shù)據(jù)輸入神經(jīng)網(wǎng)絡進行驗證,在檢驗集輸入后會得出一個誤差。在對網(wǎng)絡的訓練過程中,網(wǎng)絡的輸出會逐步逼近目標輸出,從而減小目標誤差,但是在這一訓練過程中,會利用確認樣本來檢驗網(wǎng)絡對確認樣本的擬合能力,即網(wǎng)絡在的輸出誤差不在減小(甚至增大)時,訓練終止。每次訓練的誤差通過ms e(均方誤差)計算。輸出三個影響因素的訓練誤差圖如圖2.2。
圖2.2 從左到右依次為輻射、降雨量、云層數(shù)據(jù)的訓練集誤差
從圖中可以看出,太陽輻射數(shù)據(jù)的最小誤差為0.3632;降雨量數(shù)據(jù)的誤差為0.0080007;云層數(shù)據(jù)的誤差為0.22305。云層數(shù)據(jù)和太陽輻射值的誤差相對較大。圖中三個影響因素的誤差曲線都呈現(xiàn)下降趨勢達到最小誤差之間,三個數(shù)據(jù)集間具有合理的的相關性。圖2.3的橫坐標是年份因子,縱坐標為尺度因子;每一個點代表一個數(shù)據(jù)。從圖中截取的一段數(shù)據(jù)的變化趨勢可以看出,三個訓練集大部分在相關限以內(nèi),除了個別數(shù)據(jù)的相關性較差。而數(shù)據(jù)的波動變化的正負值是相對的。分別對各個因素的各個訓練集進行相關性分析,r是相關系數(shù),絕對值始終小于1,越接近+-1相關程度越大。圖2.4分別對各個影響因素對氣候溫度的變化的相關性進行了分析,根據(jù)輻射和降雨量對溫度的相關性的擬合可以看出,訓練集的數(shù)據(jù)都集中在對角線附近,檢驗集的分布相似,測試集的數(shù)據(jù)分布較散,但是各個數(shù)據(jù)集的相關系數(shù)接近于1,說明相關性較好,說明模式是可靠的。
圖2.4 輻射(左)和降雨量(右)的相關性分析
圖2.3 從左到右為輻射、降雨量、云層數(shù)據(jù)的局部自相關
使用上述模型對2001年1月至2026年12月全球的溫度進行預測,其結(jié)果如圖2.5:其中藍色曲線為1994-2001年的全球溫度曲線,紅色曲線為對2001年以后25年的全球氣溫的預測。從圖上可以看出,預測的溫度曲線于實際的溫度曲線的變化趨勢大致相同,溫度是呈規(guī)律性的變化。同時對于2001年的溫度實際值和預測值是相吻合的,說明模型對于溫度的預測結(jié)果是可靠的。
圖2.5 模型預測的25年溫度(從上到下依次溫度輻射、降雨量和云層)
通過分析每個因素對于氣候變化的影響我們發(fā)現(xiàn),太陽輻射對于全年全球溫度變化的影響較小,但是影響的幅度較大;降雨量對于問的的影響是比較宏觀的,在較短的時間內(nèi)變化不大,但是對于長期的影響具有周期性的特征;全球云量變化對于氣候的預測是最為詳細的,能夠很好的體現(xiàn)氣候的變化波動,對于異常氣候的變化也能進行一定的預測。
本文采用Mallat小波分析和NAR人工神經(jīng)網(wǎng)絡對全球的氣候進行多尺度分析,并對今后25年全球的氣候變化趨勢進行了預測,通過結(jié)果分析我們可以知道:太陽輻射、云量變化和降雨量對于氣候變化預測是相關的,可以直接根據(jù)三個變量的變化建立預測模型;通過對25年的全球溫度進行預測發(fā)現(xiàn),在未來25年內(nèi),全球的溫度波動較小,可能是預測的年限較短的原因;通過對三個變量對氣候預測的結(jié)果我們發(fā)現(xiàn),全球云量變化對氣候的影響是最劇烈的。