張曉斌，李抗彬，郝改瑞，張小鵬

基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算研究

張曉斌1，李抗彬2*，郝改瑞3，張小鵬1

(1.運城學(xué)院，山西 運城 044000；2.西安蘭特水電測控技術(shù)有限責(zé)任公司，西安 710043；3.西安理工大學(xué)，西安 710048)

【】克服傳統(tǒng)經(jīng)驗折減系數(shù)法在計算新安江模型初始土壤蓄水量方面的缺點，并提高新安江模型在濕潤半濕潤地區(qū)的應(yīng)用效果。結(jié)合流域初始土壤蓄水量的影響因素和神經(jīng)網(wǎng)絡(luò)模型特點，提出構(gòu)建基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法。在3種輸入因子組合方式下，當BP神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點大于11時，模擬訓(xùn)練期模型應(yīng)用效果達到項目精度評價指標的甲等水平，預(yù)測檢驗期的9個樣本，均有6個以上樣本檢驗合格；當BP神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點數(shù)從4個變化到21個時，模型評價指標納什效率系數(shù)從0.51變到0.97、均方根誤差從11.77降到2.74；與采用傳統(tǒng)經(jīng)驗折減系數(shù)法計算新安江模型初始土壤蓄水量相比，采用BP神經(jīng)網(wǎng)絡(luò)模型應(yīng)用效果明顯占優(yōu)，且能克服經(jīng)驗折減系數(shù)法計算土壤初始蓄水量需要選擇流域一場暴雨或久旱未雨才能開始計算和計算過程數(shù)據(jù)不能中斷的缺點。在濕潤半濕潤地區(qū)采用BP神經(jīng)網(wǎng)絡(luò)模型計算新安江模型初始土壤蓄水量具有可行性和適用性；當神經(jīng)網(wǎng)絡(luò)輸入因子和隱含層節(jié)點數(shù)選擇合理時，模型模擬和預(yù)測精度較高。

初始土壤蓄水量；BP神經(jīng)網(wǎng)絡(luò)；新安江模型；徑流模擬

0 引言

隨著全球氣候變化劇烈活動，洪水災(zāi)害頻繁發(fā)生[1]。針對降水引起的洪水過程，全球水文學(xué)者提出不同的水文模型進行流域的降水徑流過程模擬預(yù)報。劉金濤等[2]指出水文過程模擬除水文模型本身外，模型輸入的數(shù)據(jù)質(zhì)量是影響模擬精度的關(guān)鍵因素。根據(jù)流域降雨形成徑流過程的驅(qū)動因素分析可知，影響流域降雨徑流過程的因素主要包括：氣候變化因素、流域土壤類型、植被覆蓋度、地形和地貌等[3]。如果水文模型在流域降雨徑流過程模擬中輸入的這些影響因素數(shù)據(jù)具有足夠精確度，那么模型對水文過程模擬通常會有比較好的精度[4-6]。【研究意義】對于在我國濕潤半濕潤地區(qū)應(yīng)用比較廣泛和成功的新安江模型來說，模型應(yīng)用過程中存在相同的問題。在模型參數(shù)和流域降雨量數(shù)據(jù)確定情況下，初始土壤蓄水量是影響新安江模型精度的最重要因素，而初始土壤需水量是為反映降雨前流域干旱程度而提出的概念性狀態(tài)量，并非真實的土壤含水率，因而無法通過測量獲得，但其影響因素與流域土壤含水率基本相同[2]。因此通過合理有效的方法對新安江模型初始土壤蓄水量進行確定，對新安江模型應(yīng)用和提高模擬精度具有重要的意義?！狙芯窟M展】目前新安江模型應(yīng)用過程中，流域初始土壤蓄水量主要采用新安江模型預(yù)熱法或采用經(jīng)驗折減系數(shù)法進行賦值，其他計算方法較少。但這2種方法均存在一定的局限性，如：新安江模型預(yù)熱法需要較長的預(yù)熱期；經(jīng)驗折減系數(shù)法計算時首先要選擇流域中一場大雨或久旱未雨之后才能開始計算，且計算過程數(shù)據(jù)不能中斷。而針對流域真實土壤含水率預(yù)測計算方法較多，根據(jù)許秀英等[7]統(tǒng)計，主要方法包括：經(jīng)驗法、水量平衡法、時間序列法、神經(jīng)網(wǎng)絡(luò)法、遙感監(jiān)測法等，不同的方法在進行土壤含水率預(yù)測計算時各有優(yōu)缺點。任罡[8]在汾河灌區(qū)土壤墑情預(yù)報研究中指出，經(jīng)驗法建模公式簡單，但精度較低，系數(shù)易受數(shù)據(jù)時空分布影響；侯瓊等[9]在對內(nèi)蒙古典型草原土壤水分動態(tài)研究中表明，水量平衡法計算簡單，但所需測量和計算因素較多，預(yù)測精度受各分量測定方法和技術(shù)影響明顯；白冬妹等[10]采用時間序列法進行對黃土丘陵半干旱區(qū)檸條林地的土壤含水率進行預(yù)測，結(jié)果表明時間序列法在短期內(nèi)預(yù)測效果良好，但在中長期土壤水分預(yù)測時有一定的局限性；郭慶春等[11]采用神經(jīng)網(wǎng)絡(luò)模型建立駐馬店地區(qū)深層土壤含水率預(yù)測模型，結(jié)果表明預(yù)測模型精度較高，但模型輸入因子對預(yù)測精度有影響；劉影等[12]在土壤含水率高光譜遙感定量反演研究進展中指出遙感監(jiān)測法需要積累大量的觀測數(shù)據(jù)、工作量大，模型不穩(wěn)定，適用于大范圍土壤水分預(yù)測，預(yù)測精度受地面實測資料的多少和分布影響較大?！厩腥朦c】綜上所述，以往的研究多集中在流域真實土壤含水率預(yù)測方法上，而針對新安江模型初始土壤蓄水量計算方法研究較少，且已有的方法存在一定的局限性?！緮M解決的關(guān)鍵問題】因此本文結(jié)合流域真實土壤含水率預(yù)測計算方法的特點，以及新安江模型初始土壤蓄水量的影響因素，提出基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法，期望為新安江模型在濕潤半濕潤地區(qū)應(yīng)用過程中初始土壤蓄水量計算提供新的參考方法。

1 材料與方法

1.1 研究流域概況

石頭河水庫位于陜西省寶雞市眉縣，水庫壩址以上控制流域面積為673 km2。流域?qū)儆谂瘻貛О霛駶櫞箨懶詺夂?，植被覆蓋度高，流域降水主要集中在7—9月，多年平均降水量在816 mm以上，多年平均氣溫在12 ℃左右[13]。

1.2 方法介紹

1.2.1 BP神經(jīng)網(wǎng)絡(luò)

神經(jīng)網(wǎng)絡(luò)是由神經(jīng)元之間通過相互連接而組成[14-15]，通常神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)是由輸入層、隱含層和輸出層組成。在神經(jīng)網(wǎng)絡(luò)應(yīng)用中，BP神經(jīng)網(wǎng)絡(luò)最常用，其是一種多層前饋型網(wǎng)絡(luò)[16-17]，其算法核心思想為梯度下降法，是在誤差傳遞過程中采用梯度搜索技術(shù)，以使網(wǎng)絡(luò)的實際輸出值和期望輸出值之間的誤差平方和為最小[18-19]。BP神經(jīng)網(wǎng)絡(luò)算法具體原理可參見文獻[20]。

1.2.2 新安江模型

新安江模型是趙人俊[21]提出的一種集總式概念性模型，主要適用于濕潤半濕潤地區(qū)的降雨徑流過程模擬預(yù)報。模型以蓄滿產(chǎn)流理論為主，指出在流域中包氣帶土壤含水率在沒有達到土壤田間持水率之前，流域的降水將全部被土壤吸收用來補充土壤含水率，不會形成產(chǎn)流；當流域包氣帶土壤含水率達到田間持水率時，降水將全部形成產(chǎn)流[22]。新安江模型詳細介紹參見文獻[23]。

1.2.3 經(jīng)驗折減系數(shù)法[23]

在新安江模型進行場次洪水計算時，降雨起始時間對應(yīng)的初始土壤蓄水量通常采用前期影響雨量代替。前期影響雨量常采用經(jīng)驗折減系數(shù)法進行計算，計算式為：

式中：Pa為時刻前期影響雨量（mm）；P為時段內(nèi)降雨量（mm）；為土壤含水率折減系數(shù)；Pa△為+Δ時刻前期影響雨量（mm）。

1.3 基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算模型構(gòu)建

新安江模型在流域降雨徑流過程模擬預(yù)報中，模型輸入的降雨開始前流域初始土壤蓄水量受多種因素影響[24]，這些因素之間相互作用，構(gòu)成流域內(nèi)復(fù)雜的非線性水循環(huán)系統(tǒng)。對于這種復(fù)雜的非線性系統(tǒng)模擬，神經(jīng)網(wǎng)絡(luò)方法具有很強的優(yōu)勢，因此，可采用BP神經(jīng)網(wǎng)絡(luò)構(gòu)建計算新安江模型初始土壤蓄水量的方法。

1.3.1 模型輸入和輸出的確定

在影響流域初始土壤蓄水量眾多因素中，流域地形、土壤類型隨時間變化較小，而流域植被覆蓋度和氣候影響因素均隨時間實時發(fā)生變化，因此構(gòu)建基于BP神經(jīng)網(wǎng)絡(luò)初始土壤蓄水量計算模型時，可將反映流域植被覆蓋度和氣候變化的因素作為模型的輸入因子[25-26]。由文獻[2]可知，在流域次洪模擬過程中，流域初始土壤蓄水量與流域降水前的15~22 d內(nèi)氣候變化關(guān)系最為密切，因此模型輸入因子可采用洪水場次降雨開始前的15~22 d內(nèi)的數(shù)據(jù)統(tǒng)計獲得。

1.3.2 隱含層節(jié)點數(shù)的確定

隱含層節(jié)點數(shù)過多或過少不僅關(guān)系神經(jīng)網(wǎng)絡(luò)性能好壞，同時也影響神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的復(fù)雜度，但目前神經(jīng)網(wǎng)絡(luò)模型隱含層節(jié)點數(shù)確定并無較好辦法或統(tǒng)一的計算規(guī)則[27-28]，只能根據(jù)經(jīng)驗確定BP神經(jīng)網(wǎng)絡(luò)模型隱含層節(jié)點數(shù)變化范圍在/2~21之間（其中為輸入層因子個數(shù)），然后通過試算選擇合理的網(wǎng)絡(luò)隱含層節(jié)點數(shù)。

1.3.3 激活函數(shù)的確定

神經(jīng)元激活函數(shù)主要是用來提高神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的非線性表達能力。最常見的是將S型函數(shù)和線性函數(shù)相互組合，這樣既能兼顧神經(jīng)網(wǎng)絡(luò)的非線性映射能力，又可以提高神經(jīng)網(wǎng)絡(luò)的計算性能[29]。

1.3.4 權(quán)重、閾值調(diào)整

BP神經(jīng)網(wǎng)絡(luò)在網(wǎng)絡(luò)權(quán)重和閾值調(diào)整過程中，采用不同的訓(xùn)練函數(shù)會影響網(wǎng)絡(luò)的計算耗時和收斂速度。網(wǎng)絡(luò)訓(xùn)練時選擇哪一個訓(xùn)練函數(shù)并無統(tǒng)一要求，因此可根據(jù)具體情況選擇合適的訓(xùn)練函數(shù)[30-31]。

1.4 BP神經(jīng)網(wǎng)絡(luò)初始土壤蓄水量計算方法與新安江模型耦合

通過對降雨前流域初始土壤蓄水量影響因素統(tǒng)計，構(gòu)建基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法，然后將計算的初始土壤蓄水量作為新安江模型的輸入量代入到新安江模型中，可將基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法與新安江模型耦合起來，用于流域的洪水過程降雨徑流過程模擬和預(yù)報。

1.5 研究流域數(shù)據(jù)收集及整理

根據(jù)構(gòu)建的基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算模型與新安江模型耦合應(yīng)用要求，收集研究流域內(nèi)遙感、水文以及氣象數(shù)據(jù)包括：1998—2018年不同月份流域歸一化植被指數(shù)（）[32]空間分布數(shù)據(jù)（數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心：http://www.resdc.cn）、1991—2018年水庫入庫流量數(shù)據(jù)和流域雨量站（高碼頭、杜家莊、桃川、牛心石）降水量數(shù)據(jù)以及1991—2018年流域周圍氣象站（寶雞縣、太白縣、眉縣）的氣象數(shù)據(jù)。流域各站點分布及流域位置如圖1所示。

圖1 研究流域位置及站點分布

對研究流域收集的數(shù)據(jù)進行整理和統(tǒng)計，步驟如下：

1）對雨量站數(shù)據(jù)采用泰森多邊形法計算流域面平均雨量數(shù)據(jù)序列；對氣象站數(shù)據(jù)采用算術(shù)平均數(shù)法計算流域各氣象因子數(shù)據(jù)序列；根據(jù)收集的流域歸一化植被指數(shù)空間分布數(shù)據(jù)，統(tǒng)計流域不同年份各月歸一化植被指數(shù)平均值，并計算流域多年各月歸一化植被指數(shù)平均值數(shù)據(jù)序列，然后通過樣條插值獲得流域年內(nèi)不同日期多年歸一化植被指數(shù)平均值數(shù)據(jù)序列。

2）根據(jù)計算的流域面平均雨量及收集的水庫入流量數(shù)據(jù)序列，整理流域的洪水場次數(shù)據(jù)序列。

3）根據(jù)整理的洪水場次數(shù)據(jù)序列、流域面平均雨量數(shù)據(jù)序列、流域各氣象因子數(shù)據(jù)序列以及流域多年歸一化植被指數(shù)平均值數(shù)據(jù)序列，提取各場洪水降雨起始時間前21天內(nèi)流域面平均雨量數(shù)據(jù)序列、流域各氣象因子數(shù)據(jù)序列以及多年平均歸一化植被指數(shù)平均值數(shù)據(jù)序列。

4）采用經(jīng)驗折減系數(shù)法進行流域新安江模型參數(shù)率定，在新安江模型參數(shù)確定后反推得到流域各洪水場次降雨起始時間對應(yīng)的最優(yōu)初始土壤蓄水量數(shù)據(jù)（反推計算時，使每一場洪水場次降雨起始時間對應(yīng)的初始土壤蓄水量在代入新安江模型后計算的徑流深與實測徑流深誤差小于±0.1 mm），將該數(shù)據(jù)作為建立BP神經(jīng)網(wǎng)絡(luò)初始土壤蓄水量計算模型的輸出樣本。

5）根據(jù)提取的流域各洪水場次前21天內(nèi)面平均雨量數(shù)據(jù)序列、流域各氣象因子數(shù)據(jù)序列以及多年歸一化植被指數(shù)平均值數(shù)據(jù)序列，統(tǒng)計得到建立BP神經(jīng)網(wǎng)絡(luò)初始土壤蓄水量計算模型的輸入樣本。

根據(jù)以上步驟，整理統(tǒng)計獲得流域年內(nèi)不同日期多年歸一化植被指數(shù)平均值數(shù)據(jù)序列如圖2所示；在新安江模型參數(shù)確定后，根據(jù)各洪水場次實測徑流深反推得到各洪水場次降雨起始時間對應(yīng)的流域初始土壤蓄水量如圖3所示。

圖2 流域年內(nèi)不同日期多年平均歸一化植被指數(shù)平均值

圖3 各洪水場次降雨起始時間對應(yīng)初始土壤蓄水量

由圖3可知，研究流域不同洪水場次起始時間對應(yīng)的初始土壤蓄水量分布在35~100 mm之間，數(shù)值比較離散，具有代表性。

對于BP神經(jīng)網(wǎng)絡(luò)初始土壤蓄水量計算模型的輸入樣本，通過流域不同洪水場次起始時間前21天水文氣象因素數(shù)據(jù)序列以及多年日平均數(shù)據(jù)序列統(tǒng)計得到。由于統(tǒng)計時段劃分不同，輸入因子的組成也不相同，統(tǒng)計時段劃分過小則模型因子較多，統(tǒng)計時段劃分過大則模型因子較少，因此在實際應(yīng)用過程中可根據(jù)資料收集情況選擇不同的統(tǒng)計時段進行因子統(tǒng)計。本文采用以下3種因子組合方式進行輸入樣本統(tǒng)計，第一種：降雨前第1~3、第4~6、第7~9、第10~12、第13~15、第16~18、第19~21天累積降雨量、降雨前20天平均氣溫、平均濕度、平均氣壓、平均日照時間、平均風(fēng)速、平均；第二種：降雨前第1~4、第5~8、第9~12、第13~16、第17~20天累積降雨量、其他因子（平均濕度、平均氣壓、平均日照時間、平均風(fēng)速、平均）統(tǒng)計同第一種；第三種：降雨前第1~5、第6~10、第11~15、第16~20天累積降雨量、其他因子（平均氣溫、平均濕度、平均氣壓、平均日照時間、平均風(fēng)速、平均）統(tǒng)計同第一種。

1.6 方法評價指標

為驗證基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法的效果，選擇相對誤差、納什效率系數(shù)、均方根誤差作為模型應(yīng)用好壞的評價指標[33]。

相對誤差計算式為：

納什效率系數(shù)計算式為：

均方根誤差計算式為：

2 模型應(yīng)用結(jié)果及對比分析

2.1 模型應(yīng)用結(jié)果分析

采用研究流域147場洪水過程數(shù)據(jù)、反推獲得的流域初始土壤蓄水量數(shù)據(jù)以及通過統(tǒng)計獲得的神經(jīng)網(wǎng)絡(luò)3種輸入因子組合數(shù)據(jù)構(gòu)成基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法與新安江模型耦合應(yīng)用的模擬訓(xùn)練期和預(yù)測檢驗樣本數(shù)據(jù)。模型訓(xùn)練期樣本數(shù)據(jù)來自1991—2016年，預(yù)測檢驗期樣本數(shù)據(jù)來自2017—2018年。

采用MATLAB軟件編寫相應(yīng)模型程序，并以各洪水場次計算徑流深與實測徑流深相對誤差小于實測徑流深的20%作為樣本合格的評價指標（實測徑流深的20%大于20 mm時，以20 mm為誤差允許上限），通過訓(xùn)練計算得到在不同輸入因子和隱含層節(jié)點下模擬訓(xùn)練期合格樣本數(shù)和模型評價指標變化過程線圖和預(yù)測檢驗期合格樣本數(shù)變化過程線圖，如圖4、圖5所示。

圖4 模擬訓(xùn)練期樣本合格數(shù)及模型評價指標變化過程線

由圖4可知，在模擬訓(xùn)練期，基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法在隱含層節(jié)點數(shù)小于18時，隨著隱含層節(jié)點數(shù)的增加，不同輸入因子下模型耦合應(yīng)用得到合格樣本數(shù)呈逐漸增加趨勢，當隱含層節(jié)點數(shù)超過18個以后，訓(xùn)練合格樣本數(shù)和納什效率系數(shù)增加趨勢以及均方根誤差下降趨勢明顯變緩。當隱含層節(jié)點數(shù)從4個變化到21個，模擬訓(xùn)練樣本合格率從54.3%變化到92.8%、納什效率系數(shù)從0.51變到0.97、均方根誤差從11.77降到2.74。其中當隱含層節(jié)點數(shù)大于11時，基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法與新安江模型耦合應(yīng)用樣本合格率達到水情預(yù)報項目精度評價等級的甲等水平，效果良好，表明該方法具有可行性。

圖5 預(yù)測檢驗期樣本合格數(shù)變化過程線

由圖5可知，在預(yù)測檢驗期，不同輸入因子和隱含層節(jié)點下合格樣本數(shù)均在6場以上，達到水情預(yù)報項目精度評價等級的丙等水平以上。當隱含層節(jié)點數(shù)小于18時，第一種與第二種因子組合在檢驗期的樣本合格率要好于第三種因子組合，預(yù)測效果達到水情預(yù)報項目精度評價等級的乙等水平以上。

新安江模型在濕潤半濕潤地區(qū)場次洪水過程模擬預(yù)報時，采用基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法計算洪水場次降雨開始時新安江模型輸入的初始土壤蓄水量具有可行性和適用性。

2.2 與經(jīng)驗折減系數(shù)法應(yīng)用結(jié)果對比分析

在研究流域，新安江模型參數(shù)是通過參數(shù)率定獲得。在參數(shù)率定過程中各洪水場次降雨開始時間對應(yīng)的初始土壤蓄水量是采用經(jīng)驗折算系數(shù)法計算得到，因此模型參數(shù)率定過程中目標函數(shù)最優(yōu)所對應(yīng)的洪水場次合格數(shù)即為經(jīng)驗折減系數(shù)法計算的模擬期合格樣本數(shù)。預(yù)測檢驗期合格樣本數(shù)可通過新安江模型參數(shù)帶入計算統(tǒng)計獲得。通過計算統(tǒng)計，采用經(jīng)驗折減系數(shù)法計算新安江模型初始土壤蓄水量，并代入到新安江模型得到模擬期合格樣本數(shù)為104場，預(yù)測檢驗期合格樣本數(shù)為5場。與采用BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法相比，當BP神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點數(shù)大于11時，3種輸入因子組合方式下模擬訓(xùn)練期和預(yù)測檢驗期結(jié)果均要優(yōu)于傳統(tǒng)的經(jīng)驗折減系數(shù)法。這也表明基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法具有可行性和適用性。

3 討論

3.1 基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法與傳統(tǒng)方法計算特點對比

基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法是通過統(tǒng)計各洪水場次起始時間之前15~22天的水文、氣象等要素構(gòu)建模型輸入因子，并采用網(wǎng)絡(luò)訓(xùn)練構(gòu)建模型，建模和應(yīng)用過程相對簡單，輸入因子和隱含層節(jié)點數(shù)選擇合理時模型的精度較高，與郭慶春等[11]、李興旺等[34]采用神經(jīng)網(wǎng)絡(luò)法進行流域真實土壤含水率預(yù)測計算時結(jié)論相同；而經(jīng)驗折減系數(shù)法首先需要根據(jù)流域氣象資料確定經(jīng)驗折減系數(shù)，其次，需要選擇起算值（通常為久旱之后或暴雨之后），最后在起算值確定后連續(xù)往后計算[35]；新安江模型預(yù)熱法通過模型連續(xù)計算來調(diào)整初始土壤蓄水量，計算過程由于初始土壤蓄水量初值誤差較大，會導(dǎo)致開始計算的前幾場洪水過程計算值與實測值誤差較大，因此需要較長的預(yù)熱期才能應(yīng)用。由3種方法的計算特點可知，經(jīng)驗折減系數(shù)法和新安江模型預(yù)熱法由于需要連續(xù)計算，所以對數(shù)據(jù)的連續(xù)性要求要高于基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法。

3.2 基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法不足和未來研究方向

通過驗證可知，基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法在濕潤半濕潤地區(qū)具有可行性和適用性，但是由于建模過程是通過網(wǎng)絡(luò)訓(xùn)練完成的，因此模型訓(xùn)練樣本數(shù)據(jù)的好壞和模型輸入因子的合理性會影響模型應(yīng)用效果。文中雖然通過3種因子組合方式構(gòu)建了基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算法，且模擬和預(yù)測效果要優(yōu)于經(jīng)驗折減系數(shù)法，但是否還有其他更有效的因子組合方式，需要進一步對比研究。例如：增加流域真實土壤含水率衛(wèi)星產(chǎn)品、模型模擬產(chǎn)品等作為模型輸入因子或刪減一些對計算影響較小的模型輸入因子等，從而優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)，選擇更加合理的模型輸入因子，達到提高模型的模擬預(yù)報精度目標。

4 結(jié)論

1）在研究流域新安江模型參數(shù)確定情況下，采用基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算方法進行洪水場次降雨開始時初始土壤蓄水量計算完全可行，且具有較好的適用性。

2）當BP神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點數(shù)選擇合理時，該方法在模擬期和檢驗期計算結(jié)果要明顯優(yōu)于傳統(tǒng)經(jīng)驗折減系數(shù)法，且能夠克服經(jīng)驗折減系數(shù)法的缺點。

[1] 曾維華, 程聲通. 環(huán)境災(zāi)害學(xué)引論[M]. 北京: 中國環(huán)境科學(xué)出版社, 2000.

ZENG Weihua, CHENG Shengtong. Introduction to Environmental Disaster Science [M]. Beijing: China Environment Science Press, 2000.

[2] 劉金濤, 張佳寶. 前期土壤含水量對水文模擬不確定性影響分析[J]. 冰川凍土, 2006, 28(4): 519-525.

LIU Jintao, ZHANG Jiabao. Effect of antecedent soil water content on the uncertainty of hydrological simulation[J]. Journal of Glaciology and Geocryology, 2006, 28(4): 519-525.

[3] 徐玉英. 土壤含水量計算方法[J]. 東北水利水電, 2005, 23(7): 29-30.

XU Yuying. Calculation method of soil moisture[J]. Water Resources & Hydropower of Northeast China, 2005, 23(7): 29-30.

[4] 牛宏飛, 張鐘莉莉, 孫仕軍, 等. 土壤墑情預(yù)測模型對比[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報, 2018, 23(8): 142-150.

NIU Hongfei, ZHANGZHONG Lili, SUN Shijun, et al. Comparative study on soil moisture content prediction model[J]. Journal of China Agricultural University, 2018, 23(8): 142-150.

[5] 仝兆遠, 張萬昌. 土壤水分遙感監(jiān)測的研究進展[J]. 水土保持通報, 2007, 27(4): 107-113.

TONG Zhaoyuan, ZHANG Wanchang. Progress of soil moisture monitoring by remote sensing[J]. Bulletin of Soil and Water Conservation, 2007, 27(4): 107-113.

[6] 刁萬英, 劉剛, 胡克林. 基于高光譜特征與人工神經(jīng)網(wǎng)絡(luò)模型對土壤含水量估算[J]. 光譜學(xué)與光譜分析, 2017, 37(3): 841-846.

DIAO Wanying, LIU Gang, HU Kelin. Estimation of soil water content based on hyperspectral features and the ANN model[J]. Spectroscopy and Spectral Analysis, 2017, 37(3): 841-846.

[7] 許秀英, 衣淑娟, 黃操軍. 土壤含水量預(yù)報現(xiàn)狀綜述[J]. 農(nóng)機化研究, 2013, 35(7): 11-15.

XU Xiuying, YI Shujuan, HUANG Caojun. Soil moisture content prediciton situation review[J]. Journal of Agricultural Mechanization Research, 2013, 35(7): 11-15.

[8] 任罡. 汾河灌區(qū)土壤墑情預(yù)報經(jīng)驗?zāi)Ｐ蚚J]. 科技情報開發(fā)與經(jīng)濟, 2008, 18(34): 88-90.

REN Gang. The empirical model for the forecast of soil moisture in Fenhe irrigated area[J]. Sci-Tech Information Development & Economy, 2008, 18(34): 88-90.

[9] 侯瓊, 王英舜, 楊澤龍, 等. 基于水分平衡原理的內(nèi)蒙古典型草原土壤水動態(tài)模型研究[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2011, 29(5): 197-203.

HOU Qiong, WANG Yingshun, YANG Zelong, et al. Establishing soil water dynamic model for typical steppe based on water balance principle[J]. Agricultural Research in the Arid Areas, 2011, 29(5): 197-203.

[10] 白冬妹, 郭滿才, 郭忠升, 等. 時間序列自回歸模型在土壤水分預(yù)測中的應(yīng)用研究[J]. 中國水土保持, 2014(2): 42-45.

BAI Dongmei, GUO Mancai, GUO Zhongsheng, et al. Application of time series autoregressive model in soil moisture prediction[J]. Soil And Water Conservation In China, 2014(2): 42-45.

[11] 郭慶春, 何振芳. 基于人工神經(jīng)網(wǎng)絡(luò)的土壤含水量預(yù)報模型[J]. 山西農(nóng)業(yè)科學(xué), 2012, 40(8): 892-895.

GUO Qingchun, HE Zhenfang. Forecast model of soil water content based on artificial neural network[J]. Journal of Shanxi Agricultural Sciences, 2012, 40(8): 892-895.

[12] 劉影, 姚艷敏. 土壤含水量高光譜遙感定量反演研究進展[J]. 中國農(nóng)學(xué)通報, 2016, 32(7): 127-134.

LIU Ying, YAO Yanmin. Research progress of soil moisture quantitative inversion by hyperspectral remote sensing[J]. Chinese Agricultural Science Bulletin, 2016, 32(7): 127-134.

[13] 劉拓. 智慧水庫灌區(qū)信息系統(tǒng)建設(shè)技術(shù)研究及應(yīng)用[D]. 西安: 西安理工大學(xué), 2017.

LIU Tuo. Research and application of information system construction in intelligent reservoir irrigation[D]. Xi’an: Xi’an University of Technology, 2017.

[14] ZENG X Y, SHU L, HUANG G M, et al. Triangular fuzzy series forecasting based on grey model and neural network[J]. Applied Mathematical Modelling, 2016, 40(3): 1 717-1 727.

[15] 羅黨, 王澮婷. 灰色神經(jīng)網(wǎng)絡(luò)下的多變量土壤含水量預(yù)測模型[J]. 華北水利水電大學(xué)學(xué)報(自然科學(xué)版), 2017, 38(5): 70-75.

LUO Dang, WANG Huiting. Prediction model of multivariate soil water content constructed by the grey model combined with neural network and three-parameter interval grey number[J]. Journal of North China University of Water Resources and Electric Power, 2017, 38(5): 70-75.

[16] 趙文剛, 馬孝義, 劉曉群, 等. 基于神經(jīng)網(wǎng)絡(luò)算法的廣東省典型代表站點0簡化計算模型研究[J]. 灌溉排水學(xué)報, 2019, 38(5): 91-99.

ZHAO Wengang, MA Xiaoyi, LIU Xiaoqun, et al. Using neural network model to simplify0calculation for representative stations in Guangdong Province[J]. Journal of Irrigation and Drainage, 2019, 38(5): 91-99.

[17] RASTEGAR R, HARIRI A. A step forward in studying the compact genetic algorithm[J]. Evolutionary Computation, 2006, 14(3): 277-289.

[18] 江衍銘, 郝偌楠, 李楠楠, 等. 基于靜態(tài)與動態(tài)神經(jīng)網(wǎng)絡(luò)的運河水位預(yù)報[J]. 天津大學(xué)學(xué)報, 2017, 50(3): 245-254.

CHIANG Yenming, HAO Ruonan, LI Nannan, et al. Water level forecasting by static and dynamic neural networks[J]. Journal of Tianjin University, 2017, 50(3): 245-254.

[19] 王建金, 石朋, 瞿思敏, 等. 改進BP神經(jīng)網(wǎng)絡(luò)算法在中小流域洪水預(yù)報中的應(yīng)用研究[J]. 西安理工大學(xué)學(xué)報, 2016, 32(4): 475-480.

WANG Jianjin, SHI Peng, QU Simin, et al. Application of improved BP neural network algorithm in flood forecasting in the middle and small watershed[J]. Journal of Xi’an University of Technology, 2016, 32(4): 475-480.

[20] 楊章. 基于降維LBP與葉片形狀特征的植物葉片識別方法[D]. 武漢: 湖北工業(yè)大學(xué), 2018.

YANG Zhang. Plant leaves recognition method based on dimension reduction local binary pattern and shape features of leaves[D]. Wuhan: Hubei University of Technology, 2018.

[21] 趙人俊. 流域水文模擬--新安江模型與陜北模型[M]. 北京: 水利電力出版社, 1984.

ZHAO Renjun. Watershed hydrological modeling – Xin’anjiang model and northern Shaanxi model[M] Beijing: China Water Power Press, 1984.

[22] 陳琳, 王健, 宋鵬帥, 等. 降雨對坡耕地地表結(jié)皮土壤水穩(wěn)性團聚體變化研究[J]. 灌溉排水學(xué)報, 2020, 39(1): 98-105.

CHEN Lin, WANG Jian, SONG Pengshuai, et al. Effect of rainfall on water stability aggregates of crust soil on slope surface[J]. Journal of Irrigation and Drainage, 2020, 39(1): 98-105.

[23] 包為民. 水文預(yù)報[M]. 北京: 中國水利水電出版社, 2006.

BAO Weimin. Hydrological Forecast [M]. Beijing: China Water Power Press, 2006.

[24] 劉燕. 氣候變化對干旱區(qū)資源植物開發(fā)模式的影響[D]. 南京: 南京信息工程大學(xué), 2017.

LIU Yan. Impacts of Climate Change on Plant Exploitation Patterns in Arid Areas [D]. Nanjing: Nanjing University of Information Science & Technology, 2017.

[25] 張愉, 吳柏清, 陳輝. 氣候因子變化與植被指數(shù)關(guān)聯(lián)性研究[J]. 內(nèi)江科技, 2012 (6): 35-36.

ZHANG Yu, WU Baiqing, CHENG Hui. Study on the relationship between climate factor change and vegetation index [J]. Neijiang Keji, 2012, (6): 35-36.

[26] KIM S E, SEO I W. Artificial neural network ensemble modeling with exploratory factor analysis for streamflow forecasting[J]. Journal of Hydroinformatics, 2015, 17(4): 614-639.

[27] BAHADIR E. Prediction of prospective mathematics teachers’ academic success in entering graduate education by using back-propagation neural network[J]. Journal of Education and Training Studies, 2016, 4(5): 113-122.

[28] 王嶸冰, 徐紅艷, 李波, 等. BP神經(jīng)網(wǎng)絡(luò)隱含層節(jié)點數(shù)確定方法研究[J]. 計算機技術(shù)與發(fā)展, 2018, 28(4): 31-35.

WANG Rongbing, XU Hongyan, LI Bo, et al. Research on method of determining hidden layer nodes in BP neural network[J]. Computer Technology and Development, 2018, 28(4): 31-35.

[29] 崔麗杰. 基于BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練函數(shù)選取研究[J]. 科技創(chuàng)新導(dǎo)報, 2014(36): 4, 29.

CUI Lijie. Research on the training function selection based on BP neural network [J]. Science and Technology Innovation Herald, 2014 (36): 4, 29.

[30] 胡健偉, 周玉良, 金菊良. BP神經(jīng)網(wǎng)絡(luò)洪水預(yù)報模型在洪水預(yù)報系統(tǒng)中的應(yīng)用[J]. 水文, 2015, 35(1): 20-25.

HU Jianwei, ZHOU Yuliang, JIN Juliang. Flood forecasting model on BP neural networks and its application in flood forecasting systems[J]. Journal of China Hydrology, 2015, 35(1): 20-25.

[31] 劉威, 劉尚, 白潤才, 等. 互學(xué)習(xí)神經(jīng)網(wǎng)絡(luò)訓(xùn)練方法研究[J]. 計算機學(xué)報, 2017, 40(6): 1 291-1 308.

LIU Wei, LIU Shang, BAI RunCai, et al. Research of mutual learning neural network training method[J]. Chinese Journal of Computers, 2017, 40(6): 1 291-1 308.

[32] 程杰, 楊亮彥, 黎雅楠. 2000—2018年陜北地區(qū)時空變化及其對水熱條件的響應(yīng)[J]. 灌溉排水學(xué)報, 2020, 39(5): 111-119.

CHENG Jie, YANG Liangyan, LI Ya’nan. Spatiotemporal variation inand its response to hydrothermal change from 2000 to 2018 in northern Shaanxi Province[J]. Journal of Irrigation and Drainage, 2020, 39(5): 111-119.

[33] 宋曉猛, 孔凡哲. 新安江模型和人工神經(jīng)網(wǎng)絡(luò)的耦合應(yīng)用[J]. 水土保持通報, 2010, 30(6): 135-138, 144.

SONG Xiaomeng, KONG Fanzhe. Application of xinanjiang model coupling with artificial neural networks[J]. Bulletin of Soil and Water Conservation, 2010, 30(6): 135-138, 144.

[34] 李興旺, 馮寶平. 基于BP神經(jīng)網(wǎng)絡(luò)的土壤含水量預(yù)測[J]. 水土保持學(xué)報, 2002, 16(5): 117-119.

LI Xingwang, FENG Baoping. Prediction of soil water content based on BP neural network[J]. Journal of Soil and Water Conservation, 2002, 16(5): 117-119.

[35] 葉守澤, 詹道江. 工程水文學(xué)[M]. 北京: 中國水利水電出版社, 2000.

YE Shouze, ZHAN Daojiang. Engineering Hydrology[M]. Beijing: China Water Power Press, 2000.

Using BP Network to Estimate Initial Soil Water Storage in Xin’anjiang Model

ZHANG Xiaobin1, LI Kangbin2*, HAO Gairui3, ZHANG Xiaopeng1

(1. Yuncheng University, Yuncheng 044000, China; 2. Xi’an Land Water and Electricity Measurement and Control Co.Ltd, Xi’an 710043, China; 3. Xi’an University of Technology, Xi’an 710048, China)

【】Xin’anjiang model is a hydrological model widely used for catchment modelling, but it application needs to know the initial soil moisture storage. Such initial soil moistures were traditionally estimated using the empirical reduction coefficient method which has some shortcomings, and this paper aims to present an alternative method to improve the estimate of this initial soil moisture storage when applying the model to humid and semi humid areas.【】The proposed method is based on the BP neural network and calculates the initial soil moisture storage using some easy-to-measure factors that are thought to affect moisture distribution in soil.【】Using three input factors, when the number of nodes in the hidden layer was more than 11, the accuracy index of the BP network model reached first-class level during the training stage. Of the nine samples used in the test of the model, six met the required criterion. It was also found that when the number of the nodes in the hidden layer varied between 4 and 21, the Nash–Sutcliffe efficiency coefficient in the model evaluation increased from 0.51 to 0.97, with the associated root mean square errors decreasing from 11.77 to 2.74. Compared with the traditional empirical reduction coefficient method, the BP neural network model is superior in resolving the constraints in the former, including that it needs a rainstorm or a long drought to start the calculation and that the calculation needs to be continuous in time.【】The BP neural network model proposed in this paper is feasible to calculate the initial soil water storage when applying the Xin’anjiang model to humid and semi humid areas. It can accurately estimate the initial soil water storage in a catchment if the number of input factors and the nodes in the hidden layers are rationally selected.

initial soil water storage; BP neural network; Xin’anjiang model; runoff simulation

P338.9

A

10.13522/j.cnki.ggps.2020324

1672 - 3317（2021）03 - 0015 - 08

張曉斌, 李抗彬, 郝改瑞, 等. 基于BP神經(jīng)網(wǎng)絡(luò)的新安江模型初始土壤蓄水量計算研究[J]. 灌溉排水學(xué)報, 2021, 40(3): 15-22.

ZHANG Xiaobin, LI Kangbin, HAO Gairui, et al. Using BP Network to Estimate Initial Soil Water Storage in Xin’anjiang Model[J]. Journal of Irrigation and Drainage, 2021, 40(3): 15-22.

2020-06-16

國家自然科學(xué)基金項目（51479215）；山西省水利廳科技項目（TZ2019026）；運城學(xué)院博士科研項目（YQ-2020003）

張曉斌（1977-），男，山西臨猗人。教授，主要從事流域水資源水環(huán)境模擬及污染控制研究。E-mail: 277723402@qq.com

李抗彬（1981-），男，陜西韓城人。高級工程師，主要從事流域水文水資源模擬及水利信息化研究。E-mail: winer891@sina.com

責(zé)任編輯：白芳芳