梁鵬騰，李繼清

(華北電力大學(xué) 可再生能源學(xué)院，北京 102206)

0 引 言

三峽工程是當(dāng)今世界上最大的綜合性水利樞紐工程，具有防洪、發(fā)電、航運、供水和節(jié)能減排等巨大的綜合效益。但是水庫的運行調(diào)度改變了上下游的天然水文情勢，下泄流量均一化，洪水脈沖漲落過程不明顯、次數(shù)偏少[1]，嚴重影響了長江主要經(jīng)濟魚類“四大家魚”在每年5-7月產(chǎn)卵期的繁殖活動。據(jù)《長江三峽工程生態(tài)與環(huán)境監(jiān)測公報》數(shù)據(jù)表明，三峽庫區(qū)四大家魚捕撈量由1998年的2 299 t下降至2013年的352 t，年捕撈量最少的2009年僅有95 t，壩下監(jiān)利斷面監(jiān)測到的四大家魚魚苗徑流量也呈嚴重的下降趨勢。盡管洪水脈沖與生物生命周期的相關(guān)關(guān)系機理還沒有完全被人們所掌握,但許多研究資料已顯示了這種相關(guān)性，如王俊娜等[2]研究發(fā)現(xiàn)某些洪水脈沖指標(biāo)對漁業(yè)產(chǎn)量具有正相關(guān)的促進作用；三峽水庫也曾在2011-2014年先后開展了5次試驗性生態(tài)調(diào)度工作，人工創(chuàng)造出適合四大家魚產(chǎn)卵繁殖的洪水脈沖過程，有效監(jiān)測到產(chǎn)卵活動。因此，為探究三峽水庫運行后洪水脈沖過程在四大家魚產(chǎn)卵期的變化情況，本文根據(jù)洪水脈沖相關(guān)理論提出若干個水文指標(biāo)并進行顯著性檢驗，基于檢驗結(jié)果建立水庫多目標(biāo)生態(tài)調(diào)度模型，論述延長洪水脈沖歷時的可行性，并對開展此類生態(tài)調(diào)度對水庫發(fā)電效益造成的影響進行分析，以供水庫調(diào)度決策人員參考。

1 相關(guān)理論和方法

1.1 洪水脈沖理論

洪水脈沖理論(Flood pulse concept，F(xiàn)PC)是Junk等人[3]基于在亞馬遜河和密西西比河的長期觀測和數(shù)據(jù)積累，于1989年提出的河流生態(tài)理論，該理論認為每年的洪水淹沒是河流～河漫灘區(qū)系統(tǒng)生物生存、生產(chǎn)力和交互作用的主要驅(qū)動力，河流徑流的脈沖式變化是河流洪泛區(qū)生物區(qū)系最主要的控制因子[4]。洪水脈沖帶來的水位漲落會引發(fā)不同生物的特定行為，例如鳥類遷徙、魚類的產(chǎn)卵繁殖和洄游、無脊椎動物的繁殖和遷徙等。特定的河流攜帶著既成的生物生命節(jié)律信息，洪水期間的脈沖過程將這種信息更加豐富和強烈地傳達給生物，助其完成生命活動。魚類是一種依靠水文情勢的豐枯變化來完成產(chǎn)卵、孵化和生長的生物。在巴西的Pantanal河，多種魚類已適應(yīng)在洪水脈沖高峰期進行產(chǎn)卵[5]；長江鱘則在汛期開始時感受到脈沖信息流的傳遞，進入支流活動。洪水脈沖理論在生態(tài)系統(tǒng)的改善方面也有廣泛應(yīng)用：美國基西米河針對歷史洪水脈沖特點，應(yīng)用河流～洪泛灘區(qū)的生態(tài)系統(tǒng)管理，經(jīng)過近20年的實踐和摸索，該地區(qū)的生態(tài)系統(tǒng)得到改善，生物豐富度提升明顯[6]；Bruce R等[7]對亞馬遜河洪泛區(qū)湖泊進行研究，發(fā)現(xiàn)洪水脈沖改變了洪泛區(qū)湖泊的顆粒態(tài)和可溶態(tài)氮磷的濃度；Robertson等[8]提出重建更多的天然春季洪水脈沖來維持洪泛區(qū)林木產(chǎn)量和濕地價值。本文所涉及的與三峽水庫下泄洪水脈沖特點相關(guān)的水文指標(biāo)主要有流量、歷時、變化率、發(fā)生頻率等方面，將在2.2節(jié)中提出并進行詳細研究。

1.2 Mann-Whitney U兩個獨立樣本檢驗

以三峽水庫為研究對象，對其運行前后下泄洪水過程的各水文指標(biāo)做顯著性檢驗，以判斷是否產(chǎn)生顯著性變化。水文指標(biāo)序列樣本所屬的總體分布類型往往是未知的，因此，若想判斷這種情況下兩個獨立樣本是否來自相同分布的總體，Mann-Whitney U非參數(shù)檢驗是處理此類問題的有效方法。Mann-Whitney U檢驗又稱秩和檢驗，是比較沒有配對的兩個獨立樣本的非參數(shù)檢驗，它不依賴于總體的分布形式，不考慮被研究對象為何種分布及分布類型是否已知[9]，以中位數(shù)為測度，檢驗獨立非成對樣本總體之間是否存在顯著差異，具體檢驗步驟如下：①從總體A、B中分別獨立抽取樣本nA和nB，將(nA+nB)個觀察值從小到大編序；②分別計算兩個樣本的秩和TA和TB；③建立假設(shè)。H0：兩總體A和B相對次數(shù)分布相同；H1：兩總體相對次數(shù)分布不相同；④計算檢驗統(tǒng)計量，如式(1)，式(2)。⑤設(shè)定顯著性水平，查U表得到臨界值U臨；⑥取統(tǒng)計量中較小值與臨界值進行比較，做出判斷：U≤U臨，拒絕H0，U>U臨，接受H0。

(2)

2 實例研究

2.1 水庫概況

三峽水庫是長江中上游干流的大型水庫，大壩為千年一遇設(shè)計，萬年一遇校核，正常蓄水位175 m，汛期防洪限制水位145 m，枯水期最低消落水位155 m，防洪庫容221.5 億m3，保證出力500 萬kW，裝機容量2 250 萬kW。為滿足防洪需要，水庫一般在每年5月下旬開始降低庫水位，在6月10日左右將水位降至145 m并維持低水位運行直至九月底開始蓄水。長江中游主要控制水文站宜昌站位于三峽-葛洲壩梯級水庫下游，葛洲壩屬于徑流式調(diào)節(jié)電站，對下游徑流特性影響極為有限，因此本文僅探討改善三峽水庫的調(diào)度過程對四大家魚產(chǎn)卵的促進作用。

2.2 三峽水庫水文指標(biāo)樣本檢驗

以宜昌站作為研究站點，獲取該站點1877-2013年的日流量資料，并將這137年劃分為兩個階段，即三峽運行前的1877-2003年和運行后的2004-2013年。本文重點研究水庫運行前后洪水過程在四大家魚產(chǎn)卵期的5-7月的變化情況，因此以每年的這三個月為研究時段，提出具體水文指標(biāo)如下：①5、6、7月的日平均流量；②每次洪水平均漲水歷時、平均落水歷時；③日漲水率、日落水率范圍分別在0～2 000、2 000～4 000和4 000～6 000 m3/(s·d)的漲、落水天數(shù)； ④最大累積漲水流量；⑤最長連續(xù)漲水歷時；⑥次均累積漲水流量。對上述指標(biāo)進行統(tǒng)計，生成每年的樣本數(shù)據(jù)，并使用SPSS統(tǒng)計分析軟件對水庫運行前后兩個階段的獨立樣本進行Mann-Whitney U非參數(shù)檢驗。結(jié)果見表1，平均漲水歷時、平均落水歷時兩項指標(biāo)的漸近顯著性(雙尾)取值均小于0.05，說明這兩項指標(biāo)的樣本值在三峽水庫運行前后不服從同一分布，已發(fā)生顯著性變化，統(tǒng)計值相比于運行前有所減小。其他指標(biāo)的顯著性取值均大于0.05，并未發(fā)生顯著性變化。上述檢驗結(jié)果說明三峽水庫的運行調(diào)度雖然在洪水脈沖的峰、量等方面并未對河流水文情勢造成太大影響，但洪水脈沖漲落歷時由于受水庫防洪要求、調(diào)峰需求等而有下降，造成四大家魚對洪水變化的感知不夠充分，阻礙了其產(chǎn)卵、繁殖和其他生命活動。

2.3 水庫優(yōu)化調(diào)度建模

本文基于表1的檢驗結(jié)果，建立兼顧生態(tài)的三峽水庫優(yōu)化調(diào)度模型，由于洪水過程無非漲水和落水兩種情況，平均漲水歷時的增加勢必也會帶來平均落水歷時的增加，因此目標(biāo)函數(shù)中只將漲水歷時代入計算，將發(fā)電量最大和平均漲水歷時最大耦合作為目標(biāo)函數(shù)。計算時間范圍從2004-2013年，其中每年5-7月為日調(diào)度，其他月份為月調(diào)度。

表1 水文指標(biāo)檢驗結(jié)果表Tab.1 The results of hydrological index test

(1)目標(biāo)函數(shù)：

(3)

其中：

(4)

式中：T為時段數(shù)；Em為時段m的發(fā)電量；K為出力系數(shù)；qm為時段i平均發(fā)電流量；Hm為時段i平均發(fā)電水頭；Δt為時段長；Ei等同理；w1、w2分別為目標(biāo)函數(shù)中總發(fā)電量和平均漲水天數(shù)的權(quán)重；常數(shù)4 000用于消除兩個目標(biāo)間的數(shù)量級差別，由試算得知；CY為該年5-7月的漲水次數(shù)；avg()為平均數(shù)函數(shù)；Durj表示第Y年第j場漲水的歷時。當(dāng)面臨時段為5-7月份時，ξ1=0，ξ2=1，其他月份ξ1=1，ξ2=0。

(2)約束條件：

水位約束：

Zi,min≤Zi≤Zi,max

(5)

發(fā)電流量約束：

qmin≤qi≤qmax

(6)

下泄流量約束：

Qi,min≤Qi≤Qi,max

(7)

生態(tài)約束：

Qst≤Qi

(8)

出力約束：

Nbz≤Ni≤Nyx

(9)

水量平衡約束：

Vi=Vi-1+(Qi,rk-qi-qi,qi)Δt

(10)

日水位最大變幅約束：

-ZBF≤Zi-Zi-1≤ZBF

(11)

式中：Zi,min，Zi，Zi,max分別為水庫在時段i允許最低水位，當(dāng)前水位和允許最高水位；qmin，qi，qmax分別為水輪機最小發(fā)電流量、當(dāng)前時段發(fā)電流量、最大發(fā)電流量；Qi,min，Qi，Qi,max分別為時段i允許最小下泄流量、當(dāng)前下泄流量、允許最大下泄流量；Qst為河流生態(tài)需水量，采用逐月90%保證率下的流量；Nbz，Ni，Nyx分別為水電站保證出力、當(dāng)前時段出力、水電站預(yù)想出力；qi,qi為棄水流量；ZBF為單日庫水位最大變幅。

(3)洪水脈沖歷時控制方法：統(tǒng)計宜昌站1877-2003年不同漲落歷時的歷史洪水次數(shù)和1877-2013年歷史洪水中不同日漲、落水率的天數(shù)，分別見表2、表3。分別將漲水歷時、落水歷時、漲水率、落水率的頻率累加，形成4個頻率區(qū)間。在模型中使程序生成一個介于0～1的隨機小數(shù)，根據(jù)此隨機數(shù)的大小落在漲水歷時的某個頻率區(qū)間內(nèi)，則決定當(dāng)前一次的漲水歷時；在此階段的漲水過程中，再每日生成一個隨機數(shù)，根據(jù)隨機數(shù)大小落在漲水率的某個區(qū)間內(nèi)，確定當(dāng)日的漲水率。當(dāng)前一次漲水完成后，即進入落水階段，同理生成每次的落水歷時和每日的落水率。一次落水完成后又回到漲水階段，如此循環(huán)，直至完成5-7月的日調(diào)度。需要說明的是，洪水脈沖歷時約束雖然以隨機方法確定洪水漲落歷時及漲落水率，但由于受模型中其他約束條件的共同限制作用(如水位約束、單日庫水位最大變幅約束)，且程序計算過程中的解均要符合所有約束條件，否則棄之重新生成狀態(tài)變量(時段末水位值)，因此模型實際上考慮了上游來水條件、水庫防洪、下游生態(tài)等要求，不會出現(xiàn)下泄流量不合理的跳躍，對下游防護對象造成威脅。

表2 不同漲、落水歷時的歷史洪水次數(shù)統(tǒng)計表Tab.2 Counts of different rise and fall duration of historical floods

表3 歷史洪水日漲、落水率天數(shù)統(tǒng)計表Tab.3 The statistics of the numbers of days of rise andfall rate of historical floods

以上模型采用狼群算法進行求解，狼群算法是一種新型的群體智能優(yōu)化算法，它模仿野生狼群捕食行為及其獵物分配方式，將尋優(yōu)過程分成3個步驟：①狼群在捕獵過程中，首先派出少數(shù)精壯的狼在各自活動范圍內(nèi)進行游獵搜索；②發(fā)現(xiàn)獵物的狼通過嚎叫召喚同伴向之奔襲、包圍獵物；③通過優(yōu)勝劣汰原則分配食物。狼群算法具有較強的魯棒性，能夠有效求解水庫優(yōu)化調(diào)度問題。為驗證其尋優(yōu)能力，本文以狼群算法求解三峽水庫優(yōu)化調(diào)度問題，引用2004-2013年三峽月平均入庫流量數(shù)據(jù)，以月為調(diào)度時段計算10年最大發(fā)電量，并以動態(tài)規(guī)劃法優(yōu)化結(jié)果作為比較，設(shè)置水位離散點數(shù)為500，兩種算法各運行10次，由表4的計算結(jié)果可知，狼群算法的優(yōu)化結(jié)果要稍優(yōu)于傳統(tǒng)的動態(tài)規(guī)劃法，證明了該智能算法的全局尋優(yōu)能力。

表4 兩種算法計算結(jié)果比較Tab.4 Comparision between two algorithms

3 模型結(jié)果分析

3.1 發(fā)電量優(yōu)化結(jié)果分析

根據(jù)目標(biāo)函數(shù)中不同的權(quán)重分配得到的2004-2013年三峽電站總發(fā)電量和洪水歷時優(yōu)化解如表5所示?？梢姡瑢θ龒{水庫下泄洪水的漲落水歷時進行約束時，電站的發(fā)電效益下降，且洪水歷時目標(biāo)所占權(quán)重越大，發(fā)電量下降越多。主要原因是在四大家魚產(chǎn)卵密集的5-7月份，為了讓洪水漲落歷時延長到水庫運行之前的天然狀態(tài)，水庫需提前加大下泄流量并且加大出力，造成發(fā)電流量有時會超過三峽機組的最大過流能力，引起水庫棄水，發(fā)電效益無法達到以發(fā)電量最大作為單一目標(biāo)時的優(yōu)化調(diào)度水平。但總體來看，優(yōu)化洪水漲落歷時對發(fā)電量的影響并不大，在四大家魚產(chǎn)卵期的5-7月完全優(yōu)化發(fā)電效益和完全優(yōu)化洪水脈沖歷時的發(fā)電量之差為62.6 億kWh，減少率約為0.7%。而在對兩個目標(biāo)均有側(cè)重時差值會更小，這一結(jié)果說明生態(tài)調(diào)度無需犧牲太多經(jīng)濟效益即可有效提高洪水脈沖的漲落歷時。

表5 不同權(quán)重下的多目標(biāo)優(yōu)化結(jié)果表Tab.5 The optimal results of multi-objects withdifferent weight distributions

3.2 洪水漲落歷時優(yōu)化結(jié)果分析

由表5可知，控制漲水歷時的權(quán)重w2越大，平均漲落歷時也越大。繪制不同權(quán)重分配下的歷時-發(fā)電量關(guān)系如圖1所示，圖中數(shù)字代表w1，直線分別表示三峽水庫下泄洪水、宜昌站徑流實際平均漲落歷時，可見優(yōu)化模型有效地重新分配了下泄流量，提升了洪水脈沖的持續(xù)時間。洪水漲落歷時對四大家魚的影響，除了文獻[2]中的結(jié)論外，李翀、彭靜等[11]也通過分析1964-2001年的野外監(jiān)測四大家魚魚苗發(fā)江量與宜昌站每年月的洪水過程關(guān)系，認為每年的5-6月總漲水日數(shù)以及平均每次漲水過程日數(shù)是決定家魚苗發(fā)江量多寡的一個重要的環(huán)境因子。四大家魚只在漲水期進行產(chǎn)卵活動，漲水歷時的延長有助于增加四大家魚的產(chǎn)卵時間和產(chǎn)卵量。另外，其魚卵屬于“半漂浮性”卵[12,13]，喜好多變的洪水脈沖，一旦沉入水底過久會導(dǎo)致魚卵死亡，漲水歷時的延長使魚卵更不容易沉底，孵化率提升。研究發(fā)現(xiàn)，5-7月總漲水日數(shù)達到22 d左右，每次漲水3～5 d，漲水水位上漲達3～5 m，日上漲率達0.3 m，有助于四大家魚出現(xiàn)大規(guī)模的產(chǎn)卵活動。緩慢、延長的落水過程對四大家魚的繁殖也有重要的作用，它給水體帶來充足的營養(yǎng)供魚類生長，并促進幼魚在洄游過程中獲得性腺成熟。實際上，在水庫調(diào)度中考慮延長下泄洪水的漲落歷時需結(jié)合水文預(yù)報技術(shù)、當(dāng)前來水條件以及整個四大家魚產(chǎn)卵期的洪水脈沖次數(shù)和歷時：如預(yù)報入庫徑流已然是一個歷時較長的漲水過程，則下泄流量只需保持或稍微延長該漲水過程即可；如預(yù)報入庫流量歷時和漲水幅度較小，而前期已有一個理想的洪水脈沖過程，則無需在當(dāng)前再進行生態(tài)調(diào)度。

圖1 不同權(quán)重下的漲落歷時優(yōu)化結(jié)果示意圖Fig.1 The optimal results of flood rise and fall duration with different weights

4 結(jié) 語

根據(jù)宜昌站徑流資料對三峽水庫下泄洪水的若干個水文指標(biāo)進行統(tǒng)計和非參數(shù)檢驗，發(fā)現(xiàn)洪水漲落歷時較水庫運行前顯著減小?；诖私⑺畮於嗄繕?biāo)生態(tài)調(diào)度模型，運用狼群算法求解，結(jié)果顯示不同的目標(biāo)分配權(quán)重對發(fā)電量影響不大，但洪水脈沖歷時可隨其權(quán)重的增加有明顯增加。本文的調(diào)度方法為現(xiàn)階段的水庫生態(tài)調(diào)度提供了一種新的思路，可用于促進四大家魚產(chǎn)卵的實踐探索。漁業(yè)資源的養(yǎng)護與修復(fù)是一項系統(tǒng)工程，三峽下游漁業(yè)資源下降的原因復(fù)雜多樣，流域中其他因素的影響也不可忽視，例如：漁民的過度捕撈行為、魚類產(chǎn)卵場受到破壞、河流連通性降低、水質(zhì)下降等，需要更多的專家學(xué)者關(guān)注并采取系統(tǒng)性措施，來拯救瀕危的長江漁業(yè)。

□

[1] 陳 進, 李清清. 三峽水庫試驗性運行期生態(tài)調(diào)度效果評價[J]. 長江科學(xué)院院報, 2015,(4):1-6.

[2] 王俊娜, 馮順新, 駱輝煌. 近50年長江中游洪水脈沖與漁業(yè)產(chǎn)量的關(guān)系演變[J]. 四川大學(xué)學(xué)報:工程科學(xué)版, 2012,(S2):198-205.

[3] Junk W J. The Flood Pulse Concept in River-Floodplain Systems[J]. Canadian Journal of Fisheries & Aquatic Sciences, 1989,106:110-127.

[4] 盧曉寧, 鄧 偉, 張樹清. 洪水脈沖理論及其應(yīng)用[J]. 生態(tài)學(xué)雜志, 2007,26(2):269-277.

[5] Wantzen KM, Machado FA.Seasonal isotopic changes in fish of the pantanal wetland[J].Brazil.Aquatic Sciences,2002,64:239-251.

[6] Anderson D H. Interim Hydrologic Responses to Phase I of the Kissimmee River Restoration Project, Florida[J]. Restoration Ecology, 2014,22(3):353-366.

[7] Bruce R. Forsberg, Devol A H, Richey J E, et al. Factors controlling nutrient concentrations in Amazon floodplain lakes 1[J]. Limnology & Oceanography, 1988,33(1):41-56.

[8] Robertson A I, Bacon P, Heagney G. The responses of floodplain primary production to flood frequency and timing[J]. Journal of Applied Ecology, 2001,38(1):126-136.

[9] 陳永剛, 丁麗霞, 葛宏立,等. 基于Mann-Whitney非參數(shù)檢驗和SVM的竹類高光譜識別[J]. 光譜學(xué)與光譜分析, 2011,31(11):3 010-3 013.

[10] 王建群, 賈洋洋, 肖慶元. 狼群算法在水電站水庫優(yōu)化調(diào)度中的應(yīng)用[J]. 水利水電科技進展, 2015,(3):1-4.

[11] 李 翀, 彭 靜, 廖文根. 長江中游四大家魚發(fā)江生態(tài)水文因子分析及生態(tài)水文目標(biāo)確定[J]. 中國水利水電科學(xué)研究院學(xué)報, 2006,4(3):170-176.

[12] 張曉敏, 黃道明, 謝文星,等. 漢江中下游“四大家魚”自然繁殖的生態(tài)水文特征[J]. 水生態(tài)學(xué)雜志, 2009,2(2):126-129.

[13] 陳永柏, 廖文根, 彭期冬,等. 四大家魚產(chǎn)卵水文水動力特性研究綜述[J]. 水生態(tài)學(xué)雜志, 2009,2(2):130-133.