基于多目標(biāo)調(diào)節(jié)的東莊水庫(kù)運(yùn)行可靠性研究

孫軍平,許建建,張飛儒

(陜西省水利電力勘測(cè)設(shè)計(jì)研究院,陜西 西安 710001)

為了發(fā)揮水庫(kù)綜合效益,對(duì)以防洪、發(fā)電、供水、航運(yùn)、生態(tài)等綜合利用運(yùn)行的水庫(kù),實(shí)施多目標(biāo)調(diào)度逐漸替代傳統(tǒng)單目標(biāo)調(diào)度模式、制定具有可靠性的調(diào)度規(guī)則,已成為實(shí)現(xiàn)水庫(kù)高效調(diào)度運(yùn)行、優(yōu)化水資源配置的重要手段[1]。國(guó)內(nèi)外研究學(xué)者以黃河小浪底梯級(jí)水庫(kù)、溪洛渡-向家壩梯級(jí)水庫(kù)等為研究對(duì)象,構(gòu)建了考慮供水-發(fā)電-生態(tài)效益的多目標(biāo)調(diào)度模型[2-6]。朱金峰等[7]采用混合整數(shù)多目標(biāo)優(yōu)化模型,研究了黑河中游生活、工業(yè)、農(nóng)業(yè)、生態(tài)用水的調(diào)配方案;丁偉等[5]使用兩階段多目標(biāo)水庫(kù)優(yōu)化調(diào)度模型,分析了決策者對(duì)水源保護(hù)的偏好和水庫(kù)入庫(kù)預(yù)測(cè)的不確定性對(duì)防洪和節(jié)水之間權(quán)衡影響。調(diào)度規(guī)則是水庫(kù)優(yōu)化調(diào)度成果的重要表現(xiàn)形式,通常包括調(diào)度圖或者調(diào)度函數(shù),調(diào)度圖因其直觀、簡(jiǎn)明的優(yōu)點(diǎn)在實(shí)際調(diào)度中得到了廣泛應(yīng)用[8],但當(dāng)決策因子數(shù)目增多時(shí),調(diào)度圖表征難度較大,而調(diào)度函數(shù)通過(guò)函數(shù)關(guān)系表征水庫(kù)狀態(tài)量與決策變量之間關(guān)系,最初常利用具有結(jié)構(gòu)簡(jiǎn)單、直觀特點(diǎn)的回歸方程確定,后隨著機(jī)器學(xué)習(xí)的發(fā)展,極限學(xué)習(xí)機(jī)[9]、支持向量機(jī)[10]、神經(jīng)網(wǎng)絡(luò)[11]、隨機(jī)森林[12]等方法為水庫(kù)調(diào)度函數(shù)擬合提供了新思路。

陜西省東莊水庫(kù)是以防洪減淤為主,兼顧供水、發(fā)電、改善生態(tài)的綜合利用水庫(kù),也是關(guān)中渭北地區(qū)水資源配置網(wǎng)絡(luò)中的關(guān)鍵節(jié)點(diǎn)性水庫(kù),為最大化發(fā)揮綜合效益,亟需開發(fā)一套系統(tǒng)的調(diào)度規(guī)則,構(gòu)建調(diào)水調(diào)沙情景下考慮供水-發(fā)電-生態(tài)需求的東莊水庫(kù)多目標(biāo)調(diào)度模型,基于隨機(jī)森林算法提取東莊水庫(kù)調(diào)度規(guī)則。研究成果可豐富東莊水庫(kù)調(diào)度體系,為未來(lái)東莊水庫(kù)運(yùn)行提供參考。

1 工程概況

東莊水利樞紐工程位于涇河干流最后一個(gè)峽谷段出口(張家口水文站)以上29 km,左岸為陜西省淳化縣王家山林場(chǎng),右岸為陜西省禮泉縣叱干鎮(zhèn),壩址距涇河入渭河口約87 km,距西安市約90 km,是黃河水沙調(diào)控體系的重要支流水庫(kù),東莊水庫(kù)位置及供水系統(tǒng)示意圖如圖1所示。

圖1 東莊水庫(kù)位置及供水系統(tǒng)示意圖

工程開發(fā)任務(wù)為以防洪減淤為主,兼顧供水、發(fā)電和改善生態(tài)等綜合利用。水庫(kù)正常蓄水位789 m,死水位756 m,汛期限制水位780 m,設(shè)計(jì)洪水位799 m,校核洪水位803 m。水庫(kù)總庫(kù)容32.76億m3,防洪庫(kù)容4.30億m3,調(diào)水調(diào)沙庫(kù)容3.27億m3,攔沙庫(kù)容20.53億m3,調(diào)節(jié)庫(kù)容5.78億m3。水電站裝機(jī)容量110 MW,多年平均發(fā)電量2.85億kW·h。結(jié)合堤防工程建設(shè)可將涇河下游防洪標(biāo)準(zhǔn)提高到20年一遇,顯著減輕渭河下游河道淤積,可提高涇惠渠灌區(qū)灌溉保證率,為三原縣城、西咸新區(qū)、富平縣城及工業(yè)園區(qū)等城鎮(zhèn)生活和工業(yè)供水,并承擔(dān)陜西電網(wǎng)發(fā)電任務(wù)。

2 多目標(biāo)調(diào)度模型構(gòu)建與求解

重點(diǎn)考慮供水-發(fā)電-生態(tài)需求,按照“模型構(gòu)建—優(yōu)化求解—規(guī)則提取”的研究思路展開,首先分別選取供水量最大、發(fā)電量最大和生態(tài)供水量最大為目標(biāo)函數(shù)構(gòu)建東莊水庫(kù)多目標(biāo)調(diào)度模型,其次通過(guò)加權(quán)法和遺傳算法求解獲取調(diào)度過(guò)程,最后采用隨機(jī)森林算法提取調(diào)度函數(shù)。

2.1 多目標(biāo)調(diào)度模型構(gòu)建

(1) 基本調(diào)度原則

東莊水庫(kù)調(diào)度過(guò)程中,需要兼顧協(xié)調(diào)東莊水庫(kù)生產(chǎn)生活及灌溉用水、生態(tài)用水、發(fā)電用水的矛盾,主要圍繞以下基本原則構(gòu)建調(diào)度模型:

a) 生態(tài)基流保障

為確保水庫(kù)工程下游河道內(nèi)需水,需優(yōu)先保障生態(tài)基流5.33 m3/s。

b) 供水順序

在滿足生態(tài)基流的前提下,分灌溉期(灌溉期主要集中在每年7月—8月中旬、11月—12月以及3月—4月)與非灌溉期進(jìn)行供水。灌溉期優(yōu)先保證農(nóng)業(yè)供水,城鎮(zhèn)生活及工業(yè)用水次之;非灌溉期,優(yōu)先供給城鎮(zhèn)生活及工業(yè)用水,農(nóng)業(yè)用水次之。

c) 發(fā)電模式

發(fā)電采取“以水定電”模式,即利用中線灌溉流量發(fā)電。

d) 調(diào)水調(diào)沙

除考慮供水、生態(tài)以及發(fā)電外,在正常運(yùn)用期還考慮調(diào)水調(diào)沙問(wèn)題。調(diào)沙期(7月—9月),當(dāng)水庫(kù)入庫(kù)流量大于等于400 m3/s,水庫(kù)提前敞泄,開啟排沙底孔泄流,減少水庫(kù)淤積;當(dāng)水庫(kù)入庫(kù)流量小于400 m3/s,正常供水部分泄放,非調(diào)沙期正常供水。

(2) 目標(biāo)函數(shù)

根據(jù)基本調(diào)度原則和系統(tǒng)需求,制定兼顧生態(tài)、發(fā)電和供水多目標(biāo)優(yōu)化調(diào)度模型為:

F=Max{F1,F2,F3}

(1)

式中:F、F1、F2、F3分別代表系統(tǒng)優(yōu)化總目標(biāo)、生態(tài)供水、發(fā)電和供水目標(biāo)值。

其中生態(tài)供水、發(fā)電量和供水量目標(biāo)函數(shù)為:

(2)

(3)

(4)

式中:q生態(tài)、q發(fā)電、q供水分別為生態(tài)供水流量、發(fā)電流量和供水流量;k為發(fā)電系數(shù);h(i)為第i時(shí)段的發(fā)電水頭;T為調(diào)度周期;Δt為調(diào)度步長(zhǎng);i為調(diào)度時(shí)段;j為供水路線編號(hào)。

選用加權(quán)法[13],將上述3個(gè)運(yùn)行目標(biāo)轉(zhuǎn)化成綜合目標(biāo)的調(diào)度模型,形式如下:

F4=Max{w1F1+w2F2+w3F3}

(5)

式中:F4代表采用加權(quán)法表征后的系統(tǒng)優(yōu)化總目標(biāo);w1、w2、w3分別為優(yōu)化模型中賦予生態(tài)、發(fā)電、供水的權(quán)重系數(shù)。

考慮到水庫(kù)任務(wù)重要程度,本文選用生態(tài)、發(fā)電和供水的權(quán)重系數(shù)組合為0.3、0.2和0.5。

(2) 約束條件

東莊水庫(kù)調(diào)度模型主要受水庫(kù)的水量平衡條件、水庫(kù)水位上下限、庫(kù)容上下限、死庫(kù)容、機(jī)組容量、需水過(guò)程、變量非負(fù)約束等條件約束。

2.2 調(diào)度模型求解方法

運(yùn)用遺傳算法[14]求解,將目標(biāo)函數(shù)對(duì)應(yīng)的模型指標(biāo)設(shè)置為適應(yīng)度函數(shù),如生態(tài)供水量、發(fā)電量、供水量,每個(gè)個(gè)體所攜帶的解都能計(jì)算出一個(gè)相應(yīng)的適應(yīng)度值,而每一代種群最佳適應(yīng)度值就是用于定量評(píng)價(jià)種群進(jìn)化效果的指標(biāo)。若下一代的最佳適應(yīng)度值超過(guò)歷時(shí)最佳適應(yīng)度值,則說(shuō)明種群進(jìn)化成功,替換歷史最佳適應(yīng)度值后再進(jìn)行交叉變異進(jìn)化;反之,下一代最佳適應(yīng)度值不及歷史最佳適應(yīng)度值,則不進(jìn)行替換,直接進(jìn)行交叉變異進(jìn)化。

以出庫(kù)流量為決策變量,染色體長(zhǎng)度為660,種群規(guī)模為1 000,交叉概率和變異概率分別為0.8和0.1,迭代次數(shù)為5 000次。

2.3 基于隨機(jī)森林算法的調(diào)度規(guī)則提取

2.3.1 隨機(jī)森林算法

隨機(jī)森林算法[15]是一種常用的機(jī)器學(xué)習(xí)算法,是由Breiman在2001年提出的一種集成學(xué)習(xí)方法。它的基本原理是將多個(gè)決策樹合成一個(gè)強(qiáng)大的分類器或回歸器。該算法不僅對(duì)于異常值不敏感,而且在處理高維數(shù)據(jù)時(shí)表現(xiàn)良好,具有強(qiáng)大的數(shù)據(jù)特征挖掘能力,同時(shí)可以處理分類及回歸問(wèn)題。該算法實(shí)現(xiàn)步驟如下:

(1) 隨機(jī)抽取訓(xùn)練樣本,通過(guò)Bootstrap法隨機(jī)地有放回地進(jìn)行抽樣,選擇一部分樣本,構(gòu)成一個(gè)訓(xùn)練集。

(2) 隨機(jī)篩選一部分特征變量,構(gòu)成一個(gè)候選特征集合,候選特征集合的大小是隨機(jī)森林模型的一個(gè)重要參數(shù)。

(3) 決策樹的構(gòu)建,用選定的特征變量和訓(xùn)練集來(lái)訓(xùn)練一個(gè)決策樹模型。

(4) 重復(fù)上述過(guò)程,可以生成多棵決策樹,最終組成隨機(jī)森林。

(5) 針對(duì)分類問(wèn)題,采用投票法來(lái)決定最終的分類結(jié)果,即讓每個(gè)決策樹投票,選擇得票最多的類別作為最終結(jié)果;對(duì)于回歸問(wèn)題,采用平均法來(lái)確定最終預(yù)測(cè)值。

選用優(yōu)化調(diào)度模型的輸出結(jié)果的數(shù)據(jù)樣本,利用隨機(jī)森林?jǐn)M合數(shù)據(jù)樣本之間的關(guān)系,擬合調(diào)度函數(shù),作為東莊水庫(kù)的調(diào)度規(guī)則。

2.3.2 相關(guān)性分析

采用Pearson相關(guān)系數(shù)分析方法。Pearson相關(guān)系數(shù)[16]分析方法是最常用的線性相關(guān)系數(shù)方法,主要用于表征線性相關(guān)性,最適用于相差不大的線性、連續(xù)且符合正態(tài)分布的數(shù)據(jù)。

相關(guān)系數(shù)r的絕對(duì)值越接近于1,說(shuō)明這兩個(gè)變量的相關(guān)程度越高,即這兩個(gè)變量越相似。如果r為正數(shù),說(shuō)明兩個(gè)因子之間的關(guān)系為正相關(guān),反之為負(fù)相關(guān)。

2.3.3 調(diào)度規(guī)則擬合效果評(píng)價(jià)

選用平均絕對(duì)誤差MAE、平均偏差MBE和確定性系數(shù)R2作為評(píng)價(jià)指標(biāo)量化調(diào)度規(guī)則擬合效果。平均絕對(duì)誤差衡量了預(yù)測(cè)值與真實(shí)值之間平均誤差幅度,均方根誤差用來(lái)模擬預(yù)測(cè)值與真實(shí)值之間的偏差,前三者均為越小越好;確定性系數(shù)衡量了模型擬合成果的優(yōu)劣程度,其值越大,精度越高。

3 案例應(yīng)用與結(jié)果分析

將收集到張家山站1956年—2010年月尺度徑流資料作為東莊水庫(kù)的入庫(kù)徑流,擬定水庫(kù)起調(diào)水位為死水位756 m,結(jié)果分析從優(yōu)化模型調(diào)度結(jié)果和規(guī)則提取結(jié)果兩個(gè)方面展開。

3.1 多目標(biāo)調(diào)度結(jié)果分析

基于遺傳算法獲取的調(diào)度結(jié)果從供水、發(fā)電和生態(tài)三個(gè)角度依次展開。

3.1.1 供水過(guò)程

(1) 水量關(guān)鍵指標(biāo)

經(jīng)驗(yàn)證,模型調(diào)度結(jié)果滿足水量平衡方程,結(jié)果正確,其中調(diào)度過(guò)程的多年平均值統(tǒng)計(jì)結(jié)果如表1所示、棄水變化過(guò)程如圖2所示。

表1 東莊水庫(kù)調(diào)度關(guān)鍵指標(biāo)統(tǒng)計(jì)結(jié)果 單位:億m3

圖2 東莊水庫(kù)月平均棄水量統(tǒng)計(jì)圖

如圖2所示,東莊水庫(kù)棄水主要發(fā)生在夏季,夏季棄水占比為57.5%,由于汛期來(lái)水量較大,且水庫(kù)在滿足正常供水任務(wù)以外,無(wú)法滿足汛限水位的約束,因此產(chǎn)生大量棄水,多年平均棄水率約為16.94%。建議可以在汛期增大機(jī)組出力,提高棄水利用效率。

(2) 三線供水過(guò)程

如圖3所示,中線、北線和南線多年平均缺水量分別約為0.33億m3、0.21億m3、0.07億m3,中線缺水量最大,但是考慮中線需水基數(shù)最大,計(jì)算得各線的缺水率分別約為8.87%、12.48%、15.04%。長(zhǎng)系列的供水過(guò)程中,各線最大年缺水量分別約為1.87億m3(1981年)、1.11億m3(1956年)和0.33億m3(1956年),由于1956年來(lái)水較少,北線和南線出現(xiàn)了年最大缺水,而中線最大年缺水量出現(xiàn)在1981年,是因?yàn)?981年年內(nèi)進(jìn)行了調(diào)水調(diào)沙,大量的水參與了水庫(kù)的沖沙減淤,剩余水量無(wú)法滿足供水需求。若參與調(diào)水調(diào)沙的水量過(guò)多,造成大量缺水,建議考慮此時(shí)調(diào)水調(diào)沙是否具有較高價(jià)值,應(yīng)選擇性地進(jìn)行調(diào)水調(diào)沙。

圖3 東莊水庫(kù)多年平均供水量與缺水量統(tǒng)計(jì)圖

(3) 調(diào)沙與供水關(guān)系

東莊水庫(kù)1956年—2010年調(diào)沙水量與缺水量對(duì)比情況如圖4所示,紅色表示存在調(diào)水調(diào)沙期的年份的缺水量,藍(lán)色表示無(wú)調(diào)水調(diào)沙年份的缺水量。東莊水庫(kù)存在調(diào)沙期的年份,缺水量相對(duì)較大。圖4調(diào)度數(shù)據(jù)顯示:(a)調(diào)沙水量為0的年份,即不調(diào)沙年份,缺水量較小,多年平均缺水量為0.31億m3;(b)存在調(diào)沙期的年份,多年平均缺水量為1.16億m3,調(diào)沙年份缺水量約為不調(diào)沙年份缺水量的4倍。同時(shí),由于調(diào)沙期水庫(kù)水位消落較大,水位回升較慢,此時(shí)若天然來(lái)水相對(duì)較少,還會(huì)影響后續(xù)年份的供水,產(chǎn)生連續(xù)供水不足年份,不利于正常生產(chǎn)生活。

圖4 東莊水庫(kù)1956年—2010年缺水量統(tǒng)計(jì)圖

因此,東莊水庫(kù)調(diào)水調(diào)沙期對(duì)水庫(kù)正常供水影響十分顯著,側(cè)面反映出調(diào)沙水量與供水量之間的矛盾較為突出。

3.1.2 發(fā)電量

東莊水電站的發(fā)電量主要與天然來(lái)水過(guò)程和水庫(kù)是否調(diào)沙有關(guān)。模型所得多年平均發(fā)電量為 2.50 億kW·h,略低于水庫(kù)設(shè)計(jì)多年平均發(fā)電量2.58 億kW·h。如圖5所示,若當(dāng)年參與調(diào)沙的水量多,則用于發(fā)電的水量相對(duì)較少,從而導(dǎo)致發(fā)電量降低,影響水庫(kù)的發(fā)電效益。調(diào)沙水量與發(fā)電量之間存在競(jìng)爭(zhēng)關(guān)系。不調(diào)沙年份,多年平均發(fā)電水量為7.79億m3;調(diào)沙年份,多年平均發(fā)電水量為5.31億m3;不調(diào)沙年份多年平均發(fā)電水量約為調(diào)沙年份的1.47倍。

圖5 東莊水庫(kù)1956年—2010年調(diào)沙水量與發(fā)電水量對(duì)比圖

統(tǒng)計(jì)長(zhǎng)系列結(jié)果發(fā)現(xiàn)1958、1995等年份發(fā)電量急劇下降,導(dǎo)致此問(wèn)題發(fā)生的原因是:(1)在調(diào)水調(diào)沙期間,水庫(kù)泄水,導(dǎo)致水量無(wú)法儲(chǔ)存,水位庫(kù)容降低,影響后續(xù)發(fā)電;(2)調(diào)沙期水中含沙量過(guò)高,下泄水流不能通過(guò)水輪機(jī),無(wú)法用于發(fā)電;(3)在不調(diào)沙的年份,天然水量相對(duì)較少,水庫(kù)需要?jiǎng)佑门d利庫(kù)容進(jìn)行灌溉、生活和生態(tài)供水,導(dǎo)致水庫(kù)水位消落較大,影響發(fā)電效率。

3.1.3 生態(tài)供水

根據(jù)東莊水庫(kù)調(diào)度原則,優(yōu)先滿足生態(tài)基流,模型多年平均生態(tài)供水量為1.66億m3,生態(tài)供水保證率98%,基本保證生態(tài)基流。

3.2 基于Pearson相關(guān)性分析的調(diào)度函數(shù)決策因子篩選

隱隨機(jī)優(yōu)化調(diào)度的核心是建立起時(shí)段輸入因子與決策變量之間的映射關(guān)系[12]。因此,調(diào)度規(guī)則提取的關(guān)鍵是選取恰當(dāng)?shù)臎Q策變量和輸入因子。研究東莊水庫(kù)的調(diào)度問(wèn)題,需要預(yù)測(cè)未來(lái)出庫(kù)水量,故選取t+1時(shí)段的出庫(kù)流量Q作為決策因子;為了分析輸入因子對(duì)調(diào)度決策的影響,考慮了t時(shí)段入庫(kù)流量q_int、t+1時(shí)段入庫(kù)流量q_int+1、t+1時(shí)段總需水量q_xst+1、t時(shí)段初水位Zt初、t時(shí)段末水位Zt末、t時(shí)段末庫(kù)容Vt末作為初選輸入因子集。調(diào)水調(diào)沙的結(jié)果受基本模型原則影響,出庫(kù)流量波動(dòng)很大,因此選擇不調(diào)沙的年份做調(diào)度規(guī)則提取樣本,累計(jì)41年。利用隨機(jī)森林方法建立東莊水庫(kù)調(diào)度函數(shù),其形式為:

Q=f(q_int,q_int+1,q_xst+1,Zt初,Zt末,Vt末)

(6)

對(duì)初選輸入因子集和決策變量進(jìn)行Pearson相關(guān)性分析,計(jì)算二者的相關(guān)項(xiàng)系數(shù)并加以分析,計(jì)算結(jié)果如表2所示。

表2 Pearson相關(guān)性系數(shù)計(jì)算表

對(duì)初選輸入因子進(jìn)行P值檢驗(yàn),檢驗(yàn)六種輸入因子的P值均小于0.01,表明相關(guān)性結(jié)果顯著;六種決策因子的相關(guān)系數(shù)均大于0.3,其中,t+1時(shí)段入庫(kù)流量相關(guān)系數(shù)最大為0.4852,t時(shí)段初水位相關(guān)系數(shù)最小為0.3296,因而都可作為隨機(jī)森林模型的輸入因子。

3.3 基于隨機(jī)森林算法的調(diào)度函數(shù)訓(xùn)練與驗(yàn)證

利用隨機(jī)森林模型對(duì)決策因子進(jìn)行排序,該模型最小葉子數(shù)設(shè)置為5,決策樹樹木設(shè)置為100。本文設(shè)置前300個(gè)數(shù)據(jù)序列樣本做訓(xùn)練集,后192個(gè)數(shù)據(jù)作為驗(yàn)證集。圖6展示了訓(xùn)練集和測(cè)試集預(yù)測(cè)結(jié)果對(duì)比圖,表3統(tǒng)計(jì)了樣本擬合指標(biāo)評(píng)價(jià)結(jié)果。

表3 調(diào)度規(guī)則擬合效果評(píng)價(jià)表

圖6 模型訓(xùn)練和測(cè)試結(jié)果

由圖6可知: 在測(cè)試集預(yù)測(cè)結(jié)果中可以看到,擬合出現(xiàn)了8次較大偏差,這是因?yàn)檎{(diào)度策略是根據(jù)歷史信息做出最優(yōu)決策,而徑流不確定性很大程度上會(huì)影響預(yù)測(cè)的精度;同時(shí)模型在調(diào)度函數(shù)進(jìn)行擬合的過(guò)程中,由于僅僅考慮了對(duì)樣本給出的調(diào)度過(guò)程進(jìn)行研究,而并未結(jié)合實(shí)際調(diào)度運(yùn)行過(guò)程的中的約束來(lái)考慮調(diào)度運(yùn)行,使得最終擬合出現(xiàn)偏差較大的情況。

由表3可知:相比訓(xùn)練集,測(cè)試集預(yù)測(cè)指標(biāo)R2、MAE和MBE指標(biāo)相對(duì)較差,是因?yàn)橛捎跇颖救萘亢湍Ｐ蛷?fù)雜程度不匹配,出現(xiàn)了訓(xùn)練集過(guò)擬合現(xiàn)象,但測(cè)試集確定性系數(shù)R2為0.6,說(shuō)明該模型提取的調(diào)度規(guī)則擬合效果較好。

整體而言,預(yù)測(cè)結(jié)果與模擬過(guò)程總體趨勢(shì)一致,利用隨機(jī)森林提取的調(diào)度規(guī)則可有效預(yù)測(cè)東莊水庫(kù)出庫(kù)徑流過(guò)程。

4 結(jié)論與展望

以東莊水庫(kù)為研究對(duì)象,開展了調(diào)水調(diào)沙情景下考慮供水-發(fā)電-生態(tài)需求的多目標(biāo)調(diào)度研究,主要研究結(jié)論與成果如下:

(1) 東莊水庫(kù)基本可以保障下游河道的生態(tài)基流,生態(tài)供水保證率98%,多年平均生態(tài)供水量為1.66億m3,滿足生態(tài)需求任務(wù);考慮東莊水庫(kù)供水受調(diào)水調(diào)沙的影響,中線、北線和南線多年平均缺水量分別約為0.33億m3、0.21億m3、0.07億m3,各線的缺水率分別約為8.87%、12.48%、15.04%。

(2) 通過(guò)Pearson相關(guān)系數(shù)分析方法以及隨機(jī)森林算法提取了東莊水庫(kù)的調(diào)度規(guī)則。選取了t時(shí)段入庫(kù)流量、t+1時(shí)段入庫(kù)流量、t+1時(shí)段總需水量、t時(shí)段初水位、t時(shí)段末水位、t時(shí)段末庫(kù)容作為隨機(jī)森林模型的輸入因子,選定t+1時(shí)段的出庫(kù)流量作為決策變量,基于隨機(jī)森林理論的調(diào)度函數(shù)可作為調(diào)度規(guī)則指導(dǎo)水庫(kù)運(yùn)行。

(3)通過(guò)缺水頻率分析,解決東莊水庫(kù)供水可持續(xù)問(wèn)題,需增加相應(yīng)調(diào)蓄水庫(kù),在滿足生態(tài)、發(fā)電效益的同時(shí),滿足供水保證率。