水利工程論文-岳城水庫洪水預報人工神經網絡模型實現摘要：應用visualbasic6.0編程技術，實現了人工神經網絡bp算法的程序化，并建立了岳城水庫洪水過程預報的反向傳播神經網絡模型，經檢驗，洪水預測精度較好，結果令人滿意，為岳城水庫的入庫洪水預報及調度工作提供新的思路和依據。關鍵詞：人工神經網絡BP算法VB程序洪水預報1前言人工神經網絡（ANN）又稱連接機制模型（ConnectionModel）或并行分布處理模型（ParallalDistributedModel）。作為人工智能的研究方法，目前已廣泛應用于自然科學的各個領域，應用計算機程序來模擬這種特殊的數學模型并應用于實際流域的洪水預報研究中，無疑是一種新的嘗試和有益探索。岳城水庫是海河流域南運河水系漳河上的一座大型控制性工程，入庫洪水突發性強，水猛多沙，為確保下游河北、河南、山東、天津廣大平原地區和京廣鐵路的安全，對水庫入庫洪水進行精確預報，及時采取預泄和分洪措施顯得極其重要，因此，用人工神經網絡模型模擬預報水庫的入庫洪水過程，有重要參考和借鑒意義。2BP網絡的構建人工神經網絡是一個高度復雜的非線性動力學系統，它有大量的簡單處理單元（神經元）廣泛連接而成，他對人腦的功能作了某種簡化、抽象和模擬，具有很強的非線性映射能力，其中對多層前向神經網絡bp模型的研究相對成熟，應用最為廣泛，其模型結果如圖：結構中，輸入層、隱層和輸出層神經元的個數根據具體情況設定，其中隱層層數不一，不失一般性對輸出層中只含有一個神經元的三層前向神經網絡分析如下：假設輸入層中有個神經元，隱層中有個神經元，輸出層神經元的輸出，即整個網絡的輸出為Y，網絡中輸入層的輸入分別為,，則隱層神經元的輸入分別是：(i=1,2,m)(2.1)在上式中，為隱層神經元i與輸入層神經元j的連接權，為隱層神經元的閾值，選擇函數作為隱層神經元的激發函數,則隱層神經元的輸出為：(i=1,2,m)(2.2)輸出神經元的激發函數取為線性函數，輸出層神經元的輸出及整個網絡的輸出為：(2.3)其中，Vi為輸出層神經元與隱層神經元i的連接權。定義由、組成的向量為網絡的連接權向量(ij，i，i)。設有學習樣本(，；)(=1,2,p；p為樣本數)。對某樣本(，；)在給出網絡向量后，可以通過公式(1.1)(1.3)計算出網絡的輸出值，對于樣本定義網絡的輸出誤差為：(2.4)定義誤差函數為：(2.5)(ij，i，i)隨機給出，計算式(2.5)定義的誤差值較大，網絡計算精度不高，在確定網絡結構后，通過調整(ij，i，i)的值，以逐步降低誤差，以提高網絡的計算精度，下面給出根據誤差信息調整(ij，i，i)的具體計算過程。在反向傳播算法中，是沿著誤差函數隨(ij，i，i)變化的負梯度方向對進行休整。設的修正值為：(2.6)式中：為第n次迭代計算時連接權的修正值；為前一次迭代計算時計算所得的連接權修正值；為學習率，取01間的數；為動量因子，一般取接近1的數。將式(1.4)和(1.5)代入式(1.6)中，有(2.7)定義=（,），則：(2.8)(2.9)(2.10)采用迭代式對修正計算，得到新的連接權向量。對于所有的學習樣本均按照樣本排列順序進行上述計算過程，然后固定的值，對于p個樣本分別進行正向計算，從而求出學習樣本的能量函數值：(2.11)這樣結束了一個輪次的迭代過程，當滿足某一精度要求時，就停止迭代計算，所得(ij，i，i)即為最終模型參數，否則就要進行新一輪的計算。3BP算法的VB程序實現因程序代碼太多，不再給出。網絡學習程序界面如下圖2：4洪水預報網絡模型構建41資料收集岳城水庫的入庫水文站為觀臺水文站，該站上游有清漳河匡門口水文站和濁漳河天橋段水文站，距觀臺分別為66km和64km。上游匡門口、天橋段與下游觀臺的區間流域面積為1488km2，見流域水系圖3。資料采用年鑒1962、1976、1977、1988年四次洪水和相應年份的區間時段降雨量共118組調查數據作為模型的學習訓練樣本，另取1971年和1982年兩次大洪水作為模型的檢驗數據。42預報模型構建網絡模型采用輸出層中有一個神經元的三層前向人工神經網絡，洪水預報模型的輸出節點為岳城水庫的入庫站觀臺水文站的時刻的流量，即網絡。考慮河道洪水演進時間和區間流域的產匯流時間，分別取清漳河匡門口站和天橋斷以及流域平均降雨量、作為模型的輸入節點值；隱層神經元節點數和輸入層節點數相同取為4。模型參數優化：計算中，學習率越大，學習速度會越快，但是過大時會引起振蕩效應；動量因子取得過大可能導致發散，過小則收斂速度太慢。據有關文獻介紹，取，算法收斂速度較快。本次計算取，；網絡中的初始值取(0.10.1)之間的隨機數（由VB程序產生）。樣本的歸一化處理：為了有效利用型函數的特性，以保證網絡神經元的非線性作用，對于數值型的學習樣本要進行歸一化處理。對樣本(p=1,2,p)定義，歸一化處理計算就是按照公式：(4.1)將樣本轉化為01之間的數據。對于網絡的輸出數據還應進行還原計算恢復實際值，公式為：(4.2)使用VB程序對網絡模型進行訓練學習，經102135次學習后，網絡輸出能量函數值為3.210-3，此時得到模型最優參數如表1。表1模型參數表序號11.625608-2.361247-3.2047142.846384-2.718568-1.185164-1.02973620.138017-0.041076-0.6707810.8445030.936162-2.7469741.54636231.770369-6.0486571.1342710.7454760.58331-3.0635082.5033320.250.94-0.8183131.471144-1.4842650.875319-1.995911-2.6673082.0759745-1.9730730.3