神經(jīng)網(wǎng)絡(luò)模型在降雨中長(zhǎng)期預(yù)報(bào)中的應(yīng)用比較

霍麗峰

(晉中市水文水資源勘測(cè)分局，山西 晉中 030600)

我國(guó)目前國(guó)內(nèi)人均水資源僅為世界人均水平的不到30%，全國(guó)用水緊張是近幾年制約國(guó)內(nèi)經(jīng)濟(jì)發(fā)展的重要因素之一[1-2]。降雨是補(bǔ)給地表水資源的主要途徑之一，掌握降雨的一般規(guī)律，對(duì)國(guó)內(nèi)水資源管理有著重要的意義[3]。在水文預(yù)報(bào)中，將預(yù)見(jiàn)期為3～15d的降雨預(yù)報(bào)稱為中期降雨預(yù)報(bào)，而預(yù)見(jiàn)期在15d以上的降雨預(yù)報(bào)為長(zhǎng)期降雨預(yù)報(bào)，目前國(guó)內(nèi)對(duì)中長(zhǎng)期降雨預(yù)報(bào)的研究較滯后，成熟度不高[4]。由于中長(zhǎng)期降雨預(yù)報(bào)過(guò)程中涉及到的因素較多，人文、地理等因素均可能影響降雨預(yù)報(bào)的最終結(jié)果，因此找出合理的方法進(jìn)行降雨預(yù)報(bào)，對(duì)國(guó)內(nèi)水資源評(píng)價(jià)及防災(zāi)防洪措施制定有著突出的指導(dǎo)意義[5]。

目前神經(jīng)網(wǎng)絡(luò)模型由于其計(jì)算速度快，精度高，已逐漸應(yīng)用到了降雨預(yù)報(bào)領(lǐng)域中[6]。

梁國(guó)華等[7]基于BP神經(jīng)網(wǎng)絡(luò)模型對(duì)旬降雨進(jìn)行了預(yù)報(bào)，指出BP模型模擬精度較高，但對(duì)5d以上的降雨預(yù)報(bào)模擬精度較差；劉樂(lè)等[8]基于PCA法改進(jìn)了BP神經(jīng)網(wǎng)絡(luò)模型，改進(jìn)后的模型測(cè)試預(yù)報(bào)率明顯提高，較原模型提高了50%以上；丁晶等[9]基于人工神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)月徑流量，結(jié)果表明，人工神經(jīng)網(wǎng)絡(luò)模型用作過(guò)渡期徑流預(yù)報(bào)可行，且效果優(yōu)于多元回歸方法。

RBF神經(jīng)網(wǎng)絡(luò)模型與BP神經(jīng)網(wǎng)絡(luò)模型是應(yīng)用較廣泛的模型，但2種模型均存在一定的問(wèn)題。BP網(wǎng)絡(luò)模型收斂速度慢，容易陷入局部極小值，造成計(jì)算結(jié)果的不準(zhǔn)確，RBF神經(jīng)網(wǎng)絡(luò)模型則需要較大的數(shù)據(jù)量進(jìn)行訓(xùn)練，否則其精度較低[10-11]，同時(shí)目前對(duì)降雨中長(zhǎng)期預(yù)報(bào)的模型研究較少。本文基于遺傳算法，分別對(duì)RBF神經(jīng)網(wǎng)絡(luò)模型和BP神經(jīng)網(wǎng)絡(luò)模型進(jìn)行優(yōu)化，用于預(yù)報(bào)中長(zhǎng)期降雨，同時(shí)可比較4種模型的預(yù)報(bào)精度，得出最優(yōu)模型。

1 方法與評(píng)價(jià)指標(biāo)

1.1 遺傳算法優(yōu)化神經(jīng)網(wǎng)絡(luò)模型

將遺傳算法原理應(yīng)用到BP神經(jīng)網(wǎng)絡(luò)和RBF神經(jīng)網(wǎng)絡(luò)模型中，可分為3步：

首先，確定神經(jīng)網(wǎng)絡(luò)計(jì)算結(jié)構(gòu)，確定模型計(jì)算長(zhǎng)度。本次降雨預(yù)報(bào)模型中，神經(jīng)網(wǎng)絡(luò)計(jì)算結(jié)構(gòu)指的是模型降雨實(shí)測(cè)數(shù)據(jù)，確定的實(shí)測(cè)數(shù)據(jù)回歸變化趨勢(shì)，系統(tǒng)自動(dòng)形成回歸模型結(jié)構(gòu)。

其次，每個(gè)個(gè)體通過(guò)自適應(yīng)函數(shù)計(jì)算個(gè)體適應(yīng)度值，找出最優(yōu)解?；谶z傳算法中的交叉、變異處理，優(yōu)化2個(gè)神經(jīng)網(wǎng)絡(luò)的權(quán)重賦值計(jì)算，通過(guò)模型計(jì)算權(quán)重，與輸入樣本的個(gè)體適應(yīng)度相乘，得出輸出結(jié)果，通過(guò)公式(1)計(jì)算每個(gè)樣本的個(gè)體適應(yīng)度：

(1)

式中，Y—表樣本的個(gè)體適應(yīng)度值；a—系數(shù)；Ei、Fi—第i個(gè)節(jié)點(diǎn)的期望輸出和預(yù)測(cè)輸出，n—樣本數(shù)量。

通過(guò)交叉、變異等處理，選擇合適的個(gè)體適應(yīng)度，其中每個(gè)個(gè)體的選擇概率Pi可用下式計(jì)算，最終得出輸出層結(jié)果[12]：

(2)

最后是神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的建立，通過(guò)計(jì)算每個(gè)指標(biāo)樣本的個(gè)體適應(yīng)度，選出最優(yōu)解，通過(guò)模型自動(dòng)計(jì)算每個(gè)樣本與實(shí)測(cè)值的誤差，驗(yàn)證預(yù)測(cè)值是否滿足要求，若滿足則輸出最終結(jié)果。通過(guò)上述步驟，對(duì)中長(zhǎng)期降雨進(jìn)行預(yù)報(bào)。本文用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)(GA-BP)模型和RBF神經(jīng)網(wǎng)絡(luò)模型(GA-RBF)，比較不同模型的計(jì)算精度。

1.2 評(píng)價(jià)指標(biāo)體系

本文模型模擬結(jié)果評(píng)價(jià)指標(biāo)體系采用相對(duì)均方根誤差(RMSE)、相對(duì)誤差(RE)、模型效率系數(shù)(Ens)和決定系數(shù)(R2)綜合分析精度，具體公式如下：

(3)

(4)

(5)

(6)

1.3 GPI指數(shù)確定

由于評(píng)價(jià)指標(biāo)過(guò)多，單一的評(píng)價(jià)指標(biāo)對(duì)許多具有不同統(tǒng)計(jì)指標(biāo)的模型進(jìn)行排序或比較是很困難的，因此引入全球績(jī)效指標(biāo)(GPI)對(duì)模型模擬結(jié)果進(jìn)行綜合評(píng)價(jià)。在計(jì)算時(shí)需對(duì)所有指標(biāo)的值都進(jìn)行了縮放，以防止任何特定統(tǒng)計(jì)指標(biāo)出現(xiàn)顯著影響，具體過(guò)程見(jiàn)文獻(xiàn)[13]。

GPI指數(shù)的計(jì)算公式如下：

(7)

式中，αj—常數(shù)，RMSE和RE取1，Ens和R2取-1；gj—不同指標(biāo)的縮放值的中位數(shù)；yij—不同指標(biāo)的尺度值。

2 結(jié)果與分析

2.1 中期降雨預(yù)報(bào)精度分析

表1為不同模型計(jì)算的中期降雨的精度對(duì)比。由表1可以看出，BP模型的計(jì)算精度總體低于RBF模型，用遺傳算法改進(jìn)后的GA-RBF模型計(jì)算精度同樣高于GA-BP模型。BP模型的RMSE為0.684mm/d，而相對(duì)誤差RE達(dá)到了35.78%，綜合GPI指數(shù)來(lái)看，GPI為-0.54，其值最小，表明BP模型的精度最低；RBF模型RMSE為0.501mm/d，RE為22.17%，GPI為-0.17，較BP模型精度明顯提高；改進(jìn)后的BP模型和RBF模型GPI分別達(dá)到了0.16和0.47，表明在中期降雨預(yù)報(bào)中，RBF模型的計(jì)算精度高于BP模型，GA-RBF模型計(jì)算精度高于GA-BP模型，GA-RBF模型計(jì)算精度最高，可作為降雨中期預(yù)報(bào)的模型使用。

表1 中期降雨預(yù)報(bào)計(jì)算精度對(duì)比

2.2 長(zhǎng)期降雨預(yù)報(bào)精度分析

圖1為不同模型長(zhǎng)期降雨預(yù)報(bào)趨勢(shì)分析對(duì)比。由圖1可以看出，不同模型模擬值與實(shí)測(cè)值的變化趨勢(shì)基本一致，均呈現(xiàn)開(kāi)口向下的拋物線型式，降雨主要集中在了每年的3—10月的主要作物生長(zhǎng)期，占全年降雨的68.74%左右。改進(jìn)后的GA-BP模型和GA-RBF模型更接近實(shí)測(cè)曲線，BP模型曲線在第20個(gè)月的時(shí)候誤差最高，達(dá)到了27.38%，RBF模型在第20個(gè)月和第30個(gè)月的誤差較高，分別達(dá)到了26.37%和21.74%。圖中可以清晰看出，GA-BP模型和GA-RBF模型的計(jì)算精度明顯較高。

表2為不同模型計(jì)算的長(zhǎng)期降雨的精度對(duì)比。由表2可以看出，不同模型在長(zhǎng)期降雨預(yù)報(bào)中的精度要明顯高于中期降雨預(yù)報(bào)，BP模型、RBF模型、GA-BP模型和GA-RBF模型GPI均有了明顯提高。在長(zhǎng)期降雨預(yù)報(bào)中，GA-BP模型和GA-RBF模型的計(jì)算精度基本一致，GA-RBF模型的計(jì)算精度略高，GPI達(dá)到了0.63，而B(niǎo)P模型和RBF模型的計(jì)算精度較低，GPI分別為0.27和0.35。表明在長(zhǎng)期降雨預(yù)報(bào)中，GA-RBF模型的計(jì)算精度最高，可作為降雨長(zhǎng)期預(yù)報(bào)的模型使用。

圖1 長(zhǎng)期降雨預(yù)報(bào)趨勢(shì)分析

計(jì)算模型RMSE/(mm/d)RE/%EnsR2GPI指數(shù)BP模型0.48517.480.7930.7870.27GA-BP模型0.2988.670.9710.9780.51RBF模型0.40115.730.8710.8380.35GA-RBF模型0.1975.770.9760.9690.63

3 結(jié)語(yǔ)

本文研究了降雨中長(zhǎng)期預(yù)報(bào)的標(biāo)準(zhǔn)模型，計(jì)算結(jié)果顯示通過(guò)遺傳算法優(yōu)化后的RBF神經(jīng)網(wǎng)絡(luò)模型在降雨中長(zhǎng)期預(yù)報(bào)中表現(xiàn)出了較高的精度，其全球績(jī)效指標(biāo)GPI指數(shù)最高，可作為中長(zhǎng)期降雨預(yù)報(bào)中的標(biāo)準(zhǔn)模型使用。遺傳算法通過(guò)模擬達(dá)爾文生物進(jìn)化論的自然選擇和生物進(jìn)化過(guò)程的計(jì)算模型，更方便了神經(jīng)網(wǎng)絡(luò)模型選取最優(yōu)解，同時(shí)貝葉斯概率原理可通過(guò)后驗(yàn)分布的計(jì)算結(jié)果迭代算出模型先驗(yàn)分布的計(jì)算結(jié)果，在今后的研究中可將貝葉斯原理應(yīng)用于神經(jīng)網(wǎng)絡(luò)模型中，綜合比較降雨中長(zhǎng)期降雨預(yù)報(bào)計(jì)算精度，進(jìn)一步得出標(biāo)準(zhǔn)模型。