摘要:夾巖水利樞紐的建設(shè)將對(duì)水庫(kù)上下游的生態(tài)環(huán)境產(chǎn)生一定影響,如何在保證區(qū)域防洪、灌溉、供水等多功能需求的前提下,最大程度上滿足河道生態(tài)流量要求,將成為水庫(kù)建成后在管理上的主要問題.根據(jù)夾巖水利樞紐的工程設(shè)計(jì)成果,以壩址處典型年的月均入庫(kù)量為邊界條件,綜合考慮防洪、灌溉、供水等約束條件,采用NSGA-Ⅱ算法,以?shī)A巖水利樞紐的發(fā)電量、年內(nèi)生態(tài)流量滿足率為目標(biāo),進(jìn)行了多目標(biāo)生態(tài)調(diào)度研究.結(jié)果表明:豐水年年發(fā)電量為2.791億kW·h,年內(nèi)生態(tài)流量滿足率為100%;平水年年發(fā)電量為2.116億kW·h,年內(nèi)生態(tài)流量滿足率為96.14%;枯水年年發(fā)電量為1.278億kW·h,年內(nèi)生態(tài)流量滿足率為92.27%,即通過對(duì)夾巖水利樞紐開展水庫(kù)優(yōu)化調(diào)度,可在滿足片區(qū)防洪安全、灌溉和生產(chǎn)生活用水等要求的情況下,最大程度提高不同典型年內(nèi)的生態(tài)流量滿足率、保證水電站的發(fā)電量,研究成果可為水庫(kù)調(diào)度管理部門提供決策指導(dǎo).

　　關(guān)鍵詞:夾巖水利樞紐;多目標(biāo)優(yōu)化;生態(tài)調(diào)度;NSGA-Ⅱ

　　夾巖水利樞紐工程是貴州省規(guī)模最大的水利工程,該工程的功能以供水和灌溉為主,兼顧發(fā)電、生態(tài)等綜合效益,建成后將對(duì)區(qū)域的經(jīng)濟(jì)建設(shè)提供強(qiáng)大助力.但隨著水庫(kù)的建設(shè),河流水文、水動(dòng)力條件及河流生態(tài)系統(tǒng)也會(huì)發(fā)生一定的變化[1-2],如何能在保證各種生態(tài)要素健康穩(wěn)定的條件下,最大程度上發(fā)揮工程的綜合效益是主管部門當(dāng)下及以后值得長(zhǎng)期探究的問題.

　　水庫(kù)的建設(shè)對(duì)河流生態(tài)系統(tǒng)會(huì)產(chǎn)生一定的影響已成為學(xué)者們的共識(shí)[3-5],水庫(kù)的建設(shè)會(huì)直接改變河流的水文情勢(shì)和水動(dòng)力條件,進(jìn)而降低河流的連通性,部分魚類洄游路線和產(chǎn)卵場(chǎng)發(fā)生變化,嚴(yán)重會(huì)導(dǎo)致部分物種的滅絕[6-8];水庫(kù)建設(shè)對(duì)河流水溫、水質(zhì)也會(huì)有較大的影響,水庫(kù)大壩的建設(shè)會(huì)改變水體的垂向水溫結(jié)構(gòu),且水體下泄可能會(huì)導(dǎo)致水體中的溶解氣體過飽和,引起魚類的氣泡病等[9-12].

　　針對(duì)水庫(kù)建成后的生態(tài)影響進(jìn)行生態(tài)補(bǔ)償,國(guó)內(nèi)外諸多學(xué)者開展了相關(guān)研究.Schluter[13]等最早于1971年提出了通過開展水庫(kù)調(diào)度維持河流多樣性的命題;隨后20世紀(jì)90年代,美國(guó)的田納西河流域管理局(TennesseeValleyAuthority,TVA)開展了保護(hù)水庫(kù)下游魚類及大型無(wú)脊椎動(dòng)物的生態(tài)調(diào)度實(shí)踐[14];Babel等[15]在已有實(shí)踐成果的基礎(chǔ)上,提出了以總經(jīng)濟(jì)效益最大為目標(biāo),最大程度上滿足不同部門需水的水庫(kù)優(yōu)化調(diào)度模型.國(guó)內(nèi)關(guān)于水庫(kù)生態(tài)調(diào)度的研究起步稍晚.2002年,黃河萬(wàn)家寨、三門峽、小浪底水庫(kù)開展了泥沙生態(tài)調(diào)度研究,通過人工異重流排沙,進(jìn)而減少了水庫(kù)淤積,加快了黃河口造陸過程[16-17];2003年,郭生練等[18-19]研究了丹江口水庫(kù)枯水期調(diào)度的問題;2009年,粟曉玲等[20]以保障流域生態(tài)需水、用水等要素為目標(biāo)。

　　建立了石羊河流域水資源多目標(biāo)優(yōu)化調(diào)度模型;2011年,盧有麟等[21]針對(duì)三峽梯級(jí)水庫(kù)的發(fā)電量和生態(tài)需水量等進(jìn)行了研究;2016年,劉德富等[22]以三峽水庫(kù)支流水華頻發(fā)為目標(biāo),提出了“潮汐式”生態(tài)調(diào)度方法;2018年,張連鵬[23]等為保護(hù)額爾齊斯河的魚類棲息地、河谷林草等,開展了面向生態(tài)的水庫(kù)群中長(zhǎng)期調(diào)度研究.總的來(lái)說,當(dāng)前有關(guān)水庫(kù)生態(tài)調(diào)度的研究較為豐富,但主要基于已建成運(yùn)行的水庫(kù)或已存在明顯生態(tài)影響的水庫(kù),缺少對(duì)未建、待建水庫(kù)的生態(tài)調(diào)度研究.為此,本文主要基于夾巖水利樞紐的設(shè)計(jì)參數(shù)、工程位置水文情勢(shì)和生態(tài)需求等要素,采用NSGA-Ⅱ優(yōu)化算法[24-26],對(duì)水庫(kù)建成后開展生態(tài)調(diào)度的可能性及生態(tài)調(diào)度方案進(jìn)行探討,為水庫(kù)建成后實(shí)際生態(tài)調(diào)度決策提供指導(dǎo).

　　1研究區(qū)域概況

　　夾巖水利樞紐的水庫(kù)壩址位于畢節(jié)市六沖河中游段,壩址以上流域面積4312km2,多年平均徑流量19億m3,多年平均流量為59.7m3/s,其中水庫(kù)校核洪水位為1326.01m,正常蓄水位為1323.0m,總庫(kù)容為13.23億m3,壩后電站裝機(jī)容量90MW,多年平均發(fā)電量2.198億kW·h.

　　工程的主要功能為供水和灌溉,同時(shí)兼顧發(fā)電、生態(tài)等綜合效益.水庫(kù)下游的生活、灌溉用水均由水庫(kù)供水,下游河道無(wú)航運(yùn)功能,因此下泄流量主要保證下游生態(tài)景觀用水.根據(jù)已有研究成果[1-2],工程以壩址90%保證率最小月均流量值作為壩址斷面處的生態(tài)流量成果,為12.1m3/s.根據(jù)項(xiàng)目初設(shè)成果,水庫(kù)設(shè)計(jì)調(diào)度規(guī)則為:1)優(yōu)先滿足下游河道基流,再保障供水、灌溉取水,最后多余水量經(jīng)過壩后電站機(jī)組發(fā)電產(chǎn)能后下放;2)當(dāng)天然來(lái)水流量不滿足最小下泄生態(tài)流量時(shí),按天然實(shí)際來(lái)水流量下泄.

　　2模型建立

　　夾巖水利樞紐工程是對(duì)貴州水資源開發(fā)的關(guān)鍵性工程,根據(jù)設(shè)計(jì)調(diào)度規(guī)則,在滿足安全運(yùn)行的條件下優(yōu)先滿足下游河道基流,再滿足供水、灌溉用水要求,最后利用可發(fā)電水量保證發(fā)電效益最大化,據(jù)此工程的生態(tài)調(diào)度數(shù)學(xué)模型見式(1)~(8).

　　3求解方法

　　當(dāng)前對(duì)水庫(kù)優(yōu)化調(diào)度計(jì)算的方法較多,主要有動(dòng)態(tài)規(guī)劃(DP)[27-29]、遺傳算法(GA)[30-31]、粒子群算法(PSO)[32-33]、神經(jīng)網(wǎng)絡(luò)算法(ANN)[34-35]等,根據(jù)算法的特性可將諸多算法分為兩大類,一種是以動(dòng)態(tài)規(guī)劃為代表的線性求解法,該類方法能求得數(shù)學(xué)模型的最優(yōu)解,但計(jì)算效率較低;另一類是以遺傳算法為代表的智能優(yōu)化算法,此類方法計(jì)算效率快,但無(wú)法保證計(jì)算結(jié)果為最優(yōu)解.由于本研究對(duì)夾巖水利樞紐的生態(tài)調(diào)度目的是對(duì)未建成的水利樞紐,就施行生態(tài)調(diào)度方案可行性的探討.考慮到求解的生態(tài)調(diào)度方案缺少實(shí)際運(yùn)行情況的對(duì)比分析,為此,在保證優(yōu)化調(diào)度計(jì)算效率的前提下,提出可行的生態(tài)調(diào)度方案為目標(biāo).

　　本次研究主要采用成熟的智能優(yōu)化算法對(duì)模型進(jìn)行求解,未對(duì)計(jì)算方法做出較大改進(jìn).通過對(duì)比現(xiàn)行的求解方法,確定了采用NSGA-Ⅱ算法對(duì)模型進(jìn)行求解.NSGA-Ⅱ算法是由Deb等[36-37]于2002年,通過改進(jìn)非支配排序遺傳算法(NSGA)提出的一種求解多目標(biāo)優(yōu)化計(jì)算問題的方法,該方法主要采用非支配排序、Pareto最優(yōu)等,具有良好的搜索速度和最優(yōu)序列遺傳性.根據(jù)夾巖水利樞紐生態(tài)調(diào)度模型的目標(biāo)函數(shù)和約束條件,結(jié)合NSGA-Ⅱ的計(jì)算流程,確定模型計(jì)算步驟,主要包含內(nèi)容如下:

　　1)通過隨機(jī)創(chuàng)造一個(gè)父代種群P0,通過指定一個(gè)適應(yīng)度值,采用非支配排序、交叉變異等方法,得到后代種群Q0,大小為N.2)以此類推,通過t次循環(huán)計(jì)算后,將子代Qt與父代Pt合并組成Rt,種群大小為2N,結(jié)合適應(yīng)度值和非支配排序,得到系列非支配集F1,并計(jì)算擁擠度.3)由于Rt中包含子代和父代樣本,經(jīng)過非支配排序以后的非支配集F1中所包含的樣本是最優(yōu)良的子代,也是新的父代P(t+1).4)如果F1的樣本數(shù)小于N,則向P(t+1)中補(bǔ)充非支配集F2中的較優(yōu)樣本,直到P(t+1)樣本數(shù)量為N,然后通過交叉變異得到新的子代Q(t+1),模型進(jìn)化計(jì)算流程示意.如此逐步計(jì)算,得到最優(yōu)的子代種群,進(jìn)而得到計(jì)算成果中的最優(yōu)解.

　　4計(jì)算結(jié)果與分析

　　根據(jù)夾巖水利樞紐的工程參數(shù)、調(diào)度規(guī)則、生態(tài)流量成果等,分別對(duì)不同典型水文年的來(lái)水過程進(jìn)行了優(yōu)化調(diào)度計(jì)算.參考夾巖水利樞紐項(xiàng)目設(shè)計(jì)方案中對(duì)夾巖壩址處1957—2012共55年的歷年逐月平均徑流量分析成果,確定選取1968—1969、2002—2003、1993—1994分別作為豐、平、枯各典型年的計(jì)算輸入條件。以各典型年的年徑流過程作為輸入條件,計(jì)算不同典型年來(lái)水條件下的生態(tài)優(yōu)化調(diào)度方案。

　　以?shī)A巖水利樞紐發(fā)電量和年內(nèi)生態(tài)流量滿足率最大為生態(tài)調(diào)度目標(biāo),在考慮防洪、供水、灌溉等用水約束條件下,分別對(duì)豐、平、枯不同來(lái)水條件下的調(diào)度過程進(jìn)行了計(jì)算.由計(jì)算結(jié)果可知:

　　1)通過實(shí)施生態(tài)調(diào)度,能夠在保證供水、灌溉的條件下,最大程度滿足生態(tài)流量要求,根據(jù)計(jì)算成果,豐水年的年內(nèi)生態(tài)流量滿足率為100%,平水年次之為96.14%,枯水年最低為92.27%.2)開展生態(tài)調(diào)度,能最大程度利用水能,根據(jù)計(jì)算成果,豐水年的最大發(fā)電量為2.791億kW·h,高于設(shè)計(jì)發(fā)電量成果值;平水年最大發(fā)電量為2.116億kW·h,枯水年最大發(fā)電量為1.278億kW·h,均低于設(shè)計(jì)發(fā)電量成果值.

　　3)對(duì)比夾巖水利樞紐的多年平均發(fā)電量與優(yōu)化調(diào)度計(jì)算得到的不同典型年發(fā)電量可知,豐水年水量豐沛,可以增加水庫(kù)的發(fā)電量;平水年、枯水年,水庫(kù)來(lái)水量偏低,且本研究考慮的年內(nèi)生產(chǎn)生活供水量、灌溉供水量的滿足率均為100%,高于項(xiàng)目設(shè)計(jì)的灌溉80%保證率的要求,一定程度上減少了水庫(kù)的發(fā)電量.

　　通過實(shí)施生態(tài)調(diào)度,能夠保證灌溉需水和生產(chǎn)生活用水的需求,也能最大程度滿足生態(tài)流量要求,總體上開展生態(tài)調(diào)度并不會(huì)造成整體經(jīng)濟(jì)效益的虧損.從社會(huì)生態(tài)角度分析,對(duì)夾巖水利樞紐工程展開生態(tài)調(diào)度的研究,致力于改善區(qū)域水環(huán)境,是生態(tài)文明建設(shè)的必要.

　　5結(jié)論

　　1)通過實(shí)施生態(tài)調(diào)度,在保證供水、灌溉的條件下,豐水年的最大發(fā)電量為2.791億kW·h,年內(nèi)生態(tài)流量滿足率為100%;平水年最大發(fā)電量為2.116億kW·h,年內(nèi)生態(tài)流量滿足率為96.14%;枯水年最大發(fā)電量為1.278億kW·h,年內(nèi)生態(tài)流量滿足率為92.27%.

　　2)平水年、枯水年的發(fā)電量低于設(shè)計(jì)多年平均發(fā)電量,年內(nèi)生態(tài)流量滿足率低于100%,主要是因?yàn)榭菟�年上游�?lái)水量較少.3)對(duì)夾巖水利樞紐開展生態(tài)調(diào)度,在保證灌溉需水和生產(chǎn)生活用水的情況下,能提高生態(tài)流量滿足率和發(fā)電量,有利于片區(qū)生態(tài)環(huán)境保護(hù)和社會(huì)經(jīng)濟(jì)發(fā)展.

　　參考文獻(xiàn):

　　[1]鄧偉鑄,徐婉明,劉斌,等.大型水庫(kù)不同取水方式對(duì)下游魚類生態(tài)環(huán)境影響研究:以貴州省夾巖水利樞紐工程為例[J].人民珠江,2019,40(8):57-62.

　　[2]鄧偉鑄,廉浩,等.大型水利樞紐工程最小下泄流量的確定:以貴州夾巖水利樞紐為例[J].人民珠江,2016,37(5):68-71.

　　[3]郭生練,陳炯宏,劉攀,等.水庫(kù)群聯(lián)合優(yōu)化調(diào)度研究進(jìn)展與展望[J].水科學(xué)進(jìn)展,2010,21(4):496-503.

　　[4]董哲仁,孫東亞,趙進(jìn)勇.水庫(kù)多目標(biāo)生態(tài)調(diào)度[J].水利水電技術(shù),2007,38(1):28-32.

　　[5]陳慶偉,劉蘭芬,劉昌明.筑壩對(duì)河流生態(tài)系統(tǒng)的影響及水庫(kù)生態(tài)調(diào)度研究[J].北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,43(5):578-582.

　　[6]王俊娜,李翀,廖文根.三峽-葛洲壩梯級(jí)水庫(kù)調(diào)度對(duì)壩下河流的生態(tài)水文影響[J].水力發(fā)電學(xué)報(bào),2011,30(2):84-90.

　　作者：侯業(yè)楹1,2萬(wàn)曉安1,2姚軍1,2付占魁1,2