陳立剛

(新疆額爾齊斯河流域開(kāi)發(fā)工程建設(shè)管理局，新疆 烏魯木齊 830000)

1 研究方法

水庫(kù)運(yùn)行政策在水庫(kù)運(yùn)行中非常重要，因?yàn)樗畮?kù)運(yùn)行對(duì)社會(huì)和經(jīng)濟(jì)的影響巨大。水庫(kù)運(yùn)營(yíng)商和規(guī)劃師需要規(guī)劃一個(gè)可用于確定決策的策略。水庫(kù)的運(yùn)行政策通常是根據(jù)以前的氣象和水文數(shù)據(jù)制定的。在某些情況下，運(yùn)營(yíng)政策(也稱為運(yùn)營(yíng)規(guī)則)通常用于擬建水庫(kù)早期的規(guī)劃層面。在水庫(kù)運(yùn)行中，決策是需要明智實(shí)施的重要程序之一，以平衡水的需求和供應(yīng)，實(shí)現(xiàn)最佳的社會(huì)、經(jīng)濟(jì)和環(huán)境效益，水庫(kù)泄水的早期決策問(wèn)題通常發(fā)生在不可預(yù)測(cè)的天氣條件下。因此，開(kāi)發(fā)了一些決策程序，以確定水庫(kù)系統(tǒng)中的最佳運(yùn)行規(guī)則至關(guān)重要[1-3]。

本文以新疆地區(qū)的希尼爾水庫(kù)為研究對(duì)象，希尼爾水庫(kù)是新疆地區(qū)最大的多用途水庫(kù)之一，位于新疆庫(kù)爾勒市境內(nèi)，水庫(kù)隸屬開(kāi)孔河流域，是一座注入式反調(diào)節(jié)的平原水庫(kù)，從孔雀河第一分水樞紐引水并經(jīng)庫(kù)塔干渠總干渠輸水，尾閭接孔雀河阿恰樞紐。希尼爾水庫(kù)與其他用途(供水和娛樂(lè))一起用作洪水緩解。在本研究中，使用了2010—2019年的3600個(gè)日水庫(kù)水位(WL)數(shù)據(jù)。數(shù)據(jù)經(jīng)過(guò)預(yù)處理并標(biāo)準(zhǔn)化為-1～1的范圍。水庫(kù)水位的變化，即水位的上升或下降，影響水位的階段性。表1根據(jù)專家分類和標(biāo)稱值顯示了希尼爾水庫(kù)大壩的水位。

表1 水位表示和標(biāo)稱值

該模型采用帶偏差、學(xué)習(xí)率和動(dòng)量的標(biāo)準(zhǔn)反向傳播神經(jīng)網(wǎng)絡(luò)。模型采用基于水庫(kù)水位階段變化的時(shí)間模式。以水庫(kù)水位和水位階段的變化代替水庫(kù)實(shí)際水位作為輸入模式。水庫(kù)水位的這些變化見(jiàn)式(1)：

ΔWLt=WLt-WLt-1

(1)

式中：ΔWLt為t時(shí)水庫(kù)水位的變化；WLt為t時(shí)水庫(kù)水位；WLt-1為t-1時(shí)水庫(kù)水位。

每個(gè)數(shù)據(jù)集由N個(gè)輸入列和1個(gè)輸出列組成。N等于表示不同時(shí)間t、t-1、t-2、…、t-w的水庫(kù)水位的窗口大小，其中t表示時(shí)間，w表示窗口大小。根據(jù)式(2)對(duì)輸入進(jìn)行歸一化，以將值x變換為適合范圍[C，D]。

(2)

式中：C為新的最小值；D的范圍是[-1，1]。輸出是t時(shí)水庫(kù)水位的變化，其中，“1”表示有變化，“-1”表示無(wú)變化。

使用滑動(dòng)窗口技術(shù)保存水庫(kù)水位數(shù)據(jù)的時(shí)間信息，這個(gè)過(guò)程稱為分段過(guò)程。在本研究中，形成了六個(gè)數(shù)據(jù)集，每個(gè)數(shù)據(jù)集表示不同的窗口大小。每個(gè)窗口大小表示未來(lái)的持續(xù)時(shí)間。例如，窗口大小為3表示未來(lái)3 d的時(shí)間。數(shù)據(jù)示例如表2所示，表3顯示了為每個(gè)數(shù)據(jù)集提取的實(shí)例數(shù)。

表2 窗口大小3的數(shù)據(jù)示例

表3 實(shí)例數(shù)

將每個(gè)數(shù)據(jù)集隨機(jī)分為三個(gè)數(shù)據(jù)集：訓(xùn)練集(80%)、驗(yàn)證集(10%)和測(cè)試集(10%)。訓(xùn)練集用于神經(jīng)網(wǎng)絡(luò)的訓(xùn)練階段，而驗(yàn)證集用于在訓(xùn)練期間驗(yàn)證神經(jīng)網(wǎng)絡(luò)的性能。測(cè)試集用于在訓(xùn)練完成后測(cè)試神經(jīng)網(wǎng)絡(luò)的性能。在這項(xiàng)研究中，開(kāi)發(fā)了六個(gè)神經(jīng)網(wǎng)絡(luò)模型，每個(gè)神經(jīng)網(wǎng)絡(luò)模型用一個(gè)數(shù)據(jù)集進(jìn)行訓(xùn)練。每個(gè)模型都采用隱含層、學(xué)習(xí)率和學(xué)習(xí)動(dòng)量的不同組合進(jìn)行訓(xùn)練。當(dāng)驗(yàn)證錯(cuò)誤在幾個(gè)時(shí)期內(nèi)持續(xù)出現(xiàn)時(shí)，會(huì)執(zhí)行提前停止。神經(jīng)網(wǎng)絡(luò)訓(xùn)練程序的目的是獲得最佳的結(jié)果組合。

2 結(jié)果與討論

神經(jīng)網(wǎng)絡(luò)訓(xùn)練、驗(yàn)證和測(cè)試的結(jié)果如表4所示。總體而言，訓(xùn)練的最低誤差為0.136 307，驗(yàn)證和測(cè)試的最低誤差為0。訓(xùn)練的最佳結(jié)果為92.62%，驗(yàn)證和訓(xùn)練的最佳結(jié)果均為100%。從訓(xùn)練、驗(yàn)證和測(cè)試中獲得的最高和最低結(jié)果之間存在微小差異。因此，差異表明神經(jīng)網(wǎng)絡(luò)已經(jīng)很好地學(xué)習(xí)了數(shù)據(jù)。在此基礎(chǔ)上，選擇數(shù)據(jù)集1作為水庫(kù)水位階段變化預(yù)測(cè)模型的最佳數(shù)據(jù)集。訓(xùn)練結(jié)果為92.5%，驗(yàn)證和測(cè)試結(jié)果均為100%，誤差分別為0.150 001、0和0。數(shù)據(jù)集1由窗口大小為2的50個(gè)實(shí)例組成。

表4 訓(xùn)練、驗(yàn)證和測(cè)試結(jié)果

訓(xùn)練階段獲得的網(wǎng)絡(luò)參數(shù)值如表5所示。對(duì)于數(shù)據(jù)集1，獲得的最佳結(jié)果是學(xué)習(xí)率為0.7，動(dòng)量為0.7，實(shí)現(xiàn)的最佳網(wǎng)絡(luò)架構(gòu)是4-17-1。本研究的發(fā)現(xiàn)表明，神經(jīng)網(wǎng)絡(luò)架構(gòu)4-17-1在訓(xùn)練(92.5%)、驗(yàn)證(100%)和測(cè)試(100%)期間產(chǎn)生了可接受的性能，見(jiàn)圖1。研究結(jié)果還表明，未來(lái)的最佳持續(xù)時(shí)間為2 d。建議對(duì)水庫(kù)水位變化和水庫(kù)水位階段進(jìn)行為期2 d的觀察。

圖1 訓(xùn)練、驗(yàn)證和測(cè)試的結(jié)果

表5 神經(jīng)網(wǎng)絡(luò)參數(shù)

3 結(jié) 論

水庫(kù)水位已成為水庫(kù)泄水決策的重要指標(biāo)之一，其中,最主要的信息是水庫(kù)水位的變化。基于神經(jīng)網(wǎng)絡(luò)模型對(duì)水庫(kù)的水位進(jìn)行了預(yù)測(cè)，其中,實(shí)現(xiàn)的最佳網(wǎng)絡(luò)架構(gòu)是4-17-1，預(yù)測(cè)未來(lái)的最佳持續(xù)時(shí)間為2 d。根據(jù)未來(lái)和當(dāng)前水庫(kù)水位的信息做出早期放水決策，早期放水對(duì)于為上游流入預(yù)留空間至關(guān)重要，此外，下游河流的容量將由釋放的水量控制。因此，可以避免下游河流因水庫(kù)大量放水而發(fā)生洪水的風(fēng)險(xiǎn)。水庫(kù)水位的預(yù)測(cè)考慮了水位、降雨量和流量等歷史數(shù)據(jù)的可用性。日降雨量和水位數(shù)據(jù)已用于多用途水庫(kù)預(yù)測(cè)模型。因此，這項(xiàng)研究有助于預(yù)測(cè)水庫(kù)水位的變化和水庫(kù)水位的分期，從而有助于防汛抗旱。