任小虎，朱 詠，鄭 洋，王啟優(yōu)

(甘肅省水文站，蘭州 730030)

1 研究背景

近年來(lái)，我國(guó)生態(tài)文明建設(shè)理念、理論不斷形成、建立和發(fā)展，已出臺(tái)許多明確規(guī)定，促進(jìn)生態(tài)文明建設(shè)加速發(fā)展。社會(huì)對(duì)水電站工程建設(shè)影響河道徑流過(guò)程和水文情勢(shì)變化的關(guān)注度也在提高。目前，關(guān)于水電站工程建設(shè)對(duì)下游徑流影響的研究有很多。顧穎等[1]通過(guò)對(duì)瀾滄江梯級(jí)電站的研究，發(fā)現(xiàn)梯級(jí)電站的興建對(duì)出境處徑流有調(diào)節(jié)作用，且這種調(diào)節(jié)作用在小灣水電站和糯扎渡水電站建成后，將更加明顯；郭文獻(xiàn)等[2]以梯級(jí)水庫(kù)下游宜昌水文站流量資料為研究對(duì)象，得出三峽水庫(kù)蓄水后在較大程度上改變了河流水文情勢(shì)，其改變必將影響長(zhǎng)江中下游河流水生生物生態(tài)安全；黃維東等[3-4]分析了大通河流域水能水資源開(kāi)發(fā)利用現(xiàn)狀及其對(duì)河流水文過(guò)程與生態(tài)環(huán)境的影響，結(jié)果表明梯級(jí)水電站群無(wú)序蓄放水使洪水過(guò)程由天然的平穩(wěn)狀態(tài)轉(zhuǎn)變?yōu)槿斯じ深A(yù)的劇烈變化狀態(tài)，上下游洪峰不對(duì)應(yīng)，對(duì)下游地區(qū)的防洪安全產(chǎn)生極大威脅；賈鳳[5]對(duì)豐滿電站、白山、紅石徑流分段進(jìn)行還原分析，結(jié)果表明，豐滿水電站的運(yùn)行改變了河流天然狀態(tài)年際徑流量和徑流量分布，而且人工控制對(duì)下游徑流量年內(nèi)分配具有很大的影響；紀(jì)道斌等[6]研究我國(guó)梯級(jí)水庫(kù)對(duì)水環(huán)境的累積影響，得出亟需建立統(tǒng)一的評(píng)價(jià)標(biāo)準(zhǔn)，系統(tǒng)地認(rèn)識(shí)梯級(jí)水庫(kù)建設(shè)的環(huán)境效應(yīng)。以上文獻(xiàn)，采用不同方法從不同角度研究梯級(jí)水電站對(duì)水生態(tài)或水文情勢(shì)的影響。

本文基于黑河鶯落峽水文站徑流、洪水過(guò)程及上游流域內(nèi)降水等資料，運(yùn)用雙累積曲線法、MK檢驗(yàn)法、基尼系數(shù)法等方法，重點(diǎn)分析梯級(jí)水電站對(duì)黑河上游出山口鶯落峽水文站徑流年內(nèi)、年際變化及洪水過(guò)程的影響。

2 資料與研究方法

2.1 研究區(qū)概況

黑河是我國(guó)第二大內(nèi)陸河，位于西北干旱地區(qū)，發(fā)源于祁連山北麓，由36條大小河流組成，主要由野牛溝和八寶河兩大河源匯入，與干流一起形成“兩源一干”的整體格局，流經(jīng)青海、甘肅、內(nèi)蒙古三省(區(qū))后注入尾閭東居延海、西居延海。黑河干流出山口控制站鶯落峽水文站為上中游分界點(diǎn)，斷面多年平均流量49.3 m3/s，年徑流總量15.96億m3，近來(lái)多年徑流量呈增加趨勢(shì)；上游主河道長(zhǎng)約313 km，流域面積約10 009 km2，約占總流域面積的16.1%，是黑河的產(chǎn)流區(qū)。流域內(nèi)地形山高谷深，河床陡峻，氣候陰濕，植被較好，多年平均降水量350 mm左右，經(jīng)濟(jì)活動(dòng)以牧業(yè)為主[7]。梯級(jí)水電開(kāi)發(fā)河段系黑河上游大峽谷段，長(zhǎng)95 km，首尾落差近1 000 m，平均比降0.91%，水力資源理論蘊(yùn)藏量453 MW，單位河長(zhǎng)出力4.77 MW/km，具備良好的水能開(kāi)發(fā)條件[8]。

目前，黑河干流自上游至出山口鶯落峽水文站斷面，已建成運(yùn)行的水電站依次為寶瓶河、三道灣、二龍山、大孤山、小孤山、龍首二級(jí)、龍匯、龍首一級(jí)8座水電站；正在建設(shè)的黃藏寺水利樞紐工程(黃藏寺水庫(kù)正常蓄水位2 628 m，最大壩高123 m，總庫(kù)容4.03億m3，水庫(kù)死水位2 580 m，屬大(2)型)位于東西岔交匯后下游11 km處，控制黑河上游來(lái)水量的65%～84%，形成黑河干流梯級(jí)水電開(kāi)發(fā)結(jié)構(gòu)[9]，詳見(jiàn)表1和圖1。梯級(jí)水電站工程的建設(shè)和運(yùn)行，對(duì)出山口徑流量及洪水過(guò)程產(chǎn)生一定影響。

表1 梯級(jí)水電站基本信息

2.2 資料來(lái)源

本研究采用黑河出山口鶯落峽水文站月、年徑流量、歷年最大洪峰流量、洪水過(guò)程等系列資料，共76 a(1945—2020年)，資料來(lái)源于水文局監(jiān)測(cè)資料整編刊印成果(年鑒)，資料數(shù)據(jù)可靠。降水量數(shù)據(jù)資料采用祁連、札馬什克、鶯落峽3個(gè)水文站及14個(gè)氣象站系列監(jiān)測(cè)資料，共17站、系列65 a(1956—2020年)，共計(jì)1 105站·年的數(shù)據(jù)資料，采用泰森多邊形法計(jì)算流域內(nèi)面平均降水量。水文站和主要雨量站分布位置見(jiàn)圖1。

圖1 黑河上游流域水系和水電站分布示意圖Fig.1 Distribution of water system and hydropower stations in upstream Heihe River Basin

2.3 研究方法

2.3.1 雙累積曲線法

雙累積曲線是檢驗(yàn)具有因果關(guān)系的兩個(gè)變量間關(guān)系一致性及其變化的常用方法，可用于檢驗(yàn)水文氣象要素的一致性、趨勢(shì)性。雙累積曲線是在直角坐標(biāo)系中繪制同期內(nèi)一個(gè)變量的連續(xù)累積值與另一個(gè)變量連續(xù)累積值的關(guān)系線。

梯級(jí)水電站工程主要影響和改變河道徑流演變過(guò)程，對(duì)流域下墊面條件和產(chǎn)匯流影響不大。在區(qū)域徑流量變化歸因分析中，降水-徑流雙累積曲線可定量衡量氣候變化對(duì)徑流量的影響程度，降水-徑流雙累積曲線表現(xiàn)為線性關(guān)系，即徑流量隨降水量同步變化，說(shuō)明歷史時(shí)期降水量是影響徑流量的主要因素[10-11]，同時(shí)也說(shuō)明下墊面的變化對(duì)徑流變化影響不大。

2.3.2 MK檢驗(yàn)法

MK檢驗(yàn)法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)法，可用于檢驗(yàn)時(shí)間序列的突變點(diǎn)，已被世界氣象組織推薦并廣泛使用，其優(yōu)點(diǎn)是不需要樣本遵從一定的分布、不受少數(shù)異常值的干擾、計(jì)算簡(jiǎn)便、檢驗(yàn)范圍寬、定量化程度高[12]，適用水文、氣象等非正態(tài)分布的數(shù)據(jù)。對(duì)于具有n個(gè)樣本量的時(shí)間序列xi，構(gòu)造一個(gè)秩序列，即

(1)

其中，

(2)

式中sk是第i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)。

定義統(tǒng)計(jì)變量如下：

(3)

E(sk)=n(n+1)/4 ；

(4)

Var(sk)=n(n-1)(2n+5)/72 。

(5)

式中：E(sk)、Var(sk)分別是累計(jì)數(shù)sk的均值和方差；UFk為標(biāo)準(zhǔn)正態(tài)分布統(tǒng)計(jì)量，給定顯著性水平α，若|UFk|>Uα/2，則表明序列存在明顯的趨勢(shì)變化。將時(shí)間序列xi按逆序排列，再按照上述方法計(jì)算，同時(shí)使：

(6)

式中UBk為反序列統(tǒng)計(jì)量。

繪制UFk和UBk的變化過(guò)程，若UFk>0，表明序列呈上升趨勢(shì)；UFk<0，表明呈下降趨勢(shì)，當(dāng)它們超過(guò)臨界置信水平直線時(shí)(檢驗(yàn)置信水平α=0.05時(shí)，置信水平線為±1.96)，表明上升或下降趨勢(shì)顯著。如果UFk和UBk兩條曲線出現(xiàn)交點(diǎn)，且交點(diǎn)在臨界線之間，那么交點(diǎn)對(duì)應(yīng)的時(shí)刻就是突變開(kāi)始的時(shí)間。

2.3.3 基尼系數(shù)法

基尼系數(shù)是20世紀(jì)初意大利經(jīng)濟(jì)學(xué)家基尼根據(jù)洛倫茲曲線提出的判斷國(guó)民收入分配平等程度的指標(biāo)，其實(shí)質(zhì)是對(duì)分布均勻度的量化分析，可將其應(yīng)用到其他與均勻度分析相關(guān)的各方面[13-16]。鑒于基尼系數(shù)在不均勻度量化中的優(yōu)勢(shì)，本文按照基尼系數(shù)的求值方法，參考文獻(xiàn)[13]河流水文要素年內(nèi)分布均勻度模型，計(jì)算鶯落峽水文站1945—2020年每年月徑流量年內(nèi)分配均勻度。為便于理解，本研究分配均勻度的量化值等于1減去經(jīng)濟(jì)學(xué)中所求出的基尼系數(shù)GI，該值越大，則表示年內(nèi)分配均勻度越大，即越均勻；反之則均勻度越小，越不均勻。分配均勻度的量化值公式為

uf=1-GI=SB/(SA+SB) 。

(7)

式中：uf為月徑流量年內(nèi)分配均勻度；SA為月徑流量年內(nèi)分配洛倫茲曲線上半部面積；SB為月徑流量年內(nèi)分配洛倫茲曲線下半部面積(見(jiàn)圖2)。

圖2 月徑流量年內(nèi)分配洛倫茲曲線Fig.2 Lorentz curve of annual distribution of monthly runoff

2.3.4 洪水漲落頻次及變化率

在一場(chǎng)連續(xù)的洪水過(guò)程中，洪水漲落變化的頻繁程度可用流量的漲落頻次N表示。根據(jù)文獻(xiàn)[8]，一次洪峰過(guò)程即流量從上漲到回落，漲落變化頻次N為1，如果有k次洪峰過(guò)程，漲落頻次N=2k-1。單位時(shí)間內(nèi)洪水流量的變化量表示洪水的漲落率，其表達(dá)式為

(8)

式中：q為漲落率；Qi、Qi-1分別為時(shí)間ti、ti-1對(duì)應(yīng)的流量；ΔQ為流量差；Δt為時(shí)間差。

3 結(jié)果與分析

3.1 對(duì)徑流的影響分析

3.1.1 年徑流量變化分析

繪制鶯落峽站1945—2020年年徑流量過(guò)程線和距平累積曲線(圖3)，從圖3可以得出：年徑流量呈線性增加趨勢(shì)，線性回歸系數(shù)T檢驗(yàn)的P值為7.81×10-7，小于顯著性水平0.01，故增加趨勢(shì)顯著。按照曲線出現(xiàn)連續(xù)5 a以上變化趨勢(shì)一致的原則，可將鶯落峽歷年徑流量變化劃分為3個(gè)階段：1980年以前為枯水期，1981—2001年為平水期，2002年以后為豐水期。各階段徑流量統(tǒng)計(jì)特征見(jiàn)表2。由表2可知，豐水期的徑流量多年平均值相對(duì)平水期和枯水期分別增加了19.4%和30.8%；變差系數(shù)為0.10，相對(duì)平水期和枯水期分別減小41.2%和33.3%。

圖3 年徑流量過(guò)程線和距平累積變化Fig.3 Annual runoff hydrograph and cumulative change of annual runoff anomaly

表2 黑河干流上游徑流量階段性特征

采用MK檢驗(yàn)法分析鶯落峽站1945—2020年歷年徑流量變化趨勢(shì)和突變點(diǎn)，MK趨勢(shì)變化見(jiàn)圖4。由圖4可得2005年發(fā)生突變，并于2008年突破置信水平線表現(xiàn)為顯著增加趨勢(shì)。

圖4 年徑流量系列MK統(tǒng)計(jì)量變化過(guò)程Fig.4 Changes of MK statistics of annual runoff series

繪制鶯落峽水文站斷面以上1956—2020年面降雨和年徑流量的雙累積曲線，見(jiàn)圖5。根據(jù)年徑流MK檢驗(yàn)結(jié)果，以2005年為轉(zhuǎn)折點(diǎn)，繪制1956—2004年和2005—2020年線性趨勢(shì)線，線性回歸系數(shù)T檢驗(yàn)的P值分別為3.44×10-9、8.22×10-8，均小于顯著性水平0.01，故線性關(guān)系顯著，且前后斜率變化不大，說(shuō)明流域內(nèi)降水是徑流量增加的主要因素之一。由此也可得出，鶯落峽站斷面以上黑河干流梯級(jí)水電站庫(kù)容調(diào)節(jié)能力不大，對(duì)徑流量的影響不大，徑流量的增加主要受降水或其他氣候因素的影響。

圖5 徑流與降水量的雙累積曲線Fig.5 Double cumulative curve of runoff and precipitation

3.1.2 年內(nèi)分配變化分析

黑河干流梯級(jí)水電站的運(yùn)行改變了徑流年內(nèi)分配特征。根據(jù)《黑河干流2019—2020年度水量調(diào)度工作總結(jié)》，2020年鶯落峽斷面全年來(lái)水呈兩頭豐、中間枯的特征。目前徑流量年內(nèi)分配不均勻度的定量指標(biāo)多為年內(nèi)豐、枯水期或季節(jié)性等某些特定時(shí)段內(nèi)徑流量占年徑流量的比例，難以全面反映年徑流年內(nèi)分配的不均勻度。本文采用基尼系數(shù)法計(jì)算鶯落峽1945—2020年歷年月徑流量年內(nèi)分配均勻度，歷年分配均勻度變化見(jiàn)圖6。由圖6可得UF線性回歸系數(shù)T檢驗(yàn)的P值0.003 4，小于顯著性水平0.01，故增加趨勢(shì)顯著。

圖6 年內(nèi)分配均勻度變化過(guò)程Fig.6 Change of distribution uniformity in a year

采用MK趨勢(shì)檢驗(yàn)法對(duì)歷年月徑流量年內(nèi)分配均勻度趨勢(shì)及突變點(diǎn)檢驗(yàn)分析，MK趨勢(shì)變化過(guò)程見(jiàn)圖7。由圖7可得，在2001年年內(nèi)分配均勻度的統(tǒng)計(jì)量UFk值開(kāi)始上升，到2009年突破置信水平線達(dá)到最大值，2009—2016年在置信水平線附近變動(dòng)，2017年再次突破置信水平線并發(fā)生突變，2017—2020年年內(nèi)分配均勻度呈顯著增加趨勢(shì)。

圖7 年內(nèi)分配均勻度系列MK統(tǒng)計(jì)量變化過(guò)程Fig.7 Change of MK statistics of distribution uniformity series

3.2 對(duì)洪水過(guò)程的影響分析

3.2.1 年最大洪峰流量變化分析

繪制1948—2020年歷年最大洪峰流量變化過(guò)程(圖8)，線性回歸系數(shù)T檢驗(yàn)的P值為0.87，大于顯著性水平0.01，故增加趨勢(shì)不顯著。采有MK檢驗(yàn)最大洪峰流量的變化趨勢(shì)和突變情況，結(jié)果見(jiàn)圖9。由圖9可知，最大洪峰流量在2002年發(fā)生轉(zhuǎn)折，這與龍首一級(jí)水電站建成2001年7月建成投入運(yùn)行密切相關(guān)(該電站距鶯落峽水文站斷面最近，對(duì)最大洪峰流量的影響亦最大，其調(diào)度運(yùn)行方案直接影響鶯落峽水文站徑流過(guò)程的平穩(wěn)、漲落變化等特征)，UFk值在2011年開(kāi)始為正，呈增長(zhǎng)趨勢(shì)，但未突破置信水平，故增長(zhǎng)趨勢(shì)不顯著。轉(zhuǎn)折點(diǎn)前最大洪峰流量均值為510 m3/s，轉(zhuǎn)折點(diǎn)后均值為512 m3/s，增加0.39%。

圖8 年最大洪峰流量過(guò)程線Fig.8 Hydrograph of annual maximum peak discharge over years

圖9 年最大洪峰流量MK統(tǒng)計(jì)量變化過(guò)程Fig.9 MK test results of annual maximum peak discharge over years

采用P-Ⅲ型曲線適線法，對(duì)鶯落峽站歷年最大洪峰流量序列進(jìn)行頻率分析，該站最大洪峰流量為1 280 m3/s，變差系數(shù)CV為0.41，偏差系數(shù)CS與變差系數(shù)CV的比值CS/CV為3.50，不同頻率和重現(xiàn)期的最大洪峰流量見(jiàn)表3。

表3 歷年最大洪峰流量頻率分析結(jié)果 Table 3 Frequency analysis results of annual maximum peak discharge

3.2.2 洪水過(guò)程漲落率變化分析

以水電站(龍首一級(jí))開(kāi)始建設(shè)時(shí)間2001年為分界點(diǎn)，點(diǎn)繪水電站建設(shè)之前的4場(chǎng)洪峰流量較大的洪水過(guò)程線，包括1960年、1974年、1980年、1996年；點(diǎn)繪水電站建設(shè)以后6場(chǎng)洪峰流量較大的洪水過(guò)程線，包括2002年、2008年、2014年、2016年、2018年、2020年，見(jiàn)圖10。由圖10可知，水電站建設(shè)以前，天然洪水過(guò)程呈陡漲緩落特點(diǎn)；水電站建設(shè)以后，受水電站(群)蓄放水影響，呈明顯的陡漲直落特點(diǎn)。

梯級(jí)水電站在調(diào)度運(yùn)行中，依據(jù)河道上游來(lái)水量大小和電網(wǎng)調(diào)度需要，頻繁蓄放水，致使洪水過(guò)程由“陡漲緩落”的形態(tài)轉(zhuǎn)變?yōu)椤岸笣q直落、直漲直落或平直過(guò)程”，平緩變化的天然退水過(guò)程形態(tài)消失。主要表現(xiàn)為洪峰時(shí)間推遲、峰值變小、洪水過(guò)程拉長(zhǎng)。隨著水電站的增多及黃藏寺水庫(kù)的建成運(yùn)行，這種影響更加明顯。

經(jīng)分析，以上10場(chǎng)洪水過(guò)程的最大上漲率、最大回落率、變化頻次及水電站數(shù)量見(jiàn)表4。由表4可知，受梯級(jí)水電站閘門(mén)頻繁升降的影響，洪水過(guò)程的漲落率和變化頻次明顯增加。

4 結(jié)語(yǔ)與建議

4.1 結(jié) 語(yǔ)

黑河干流梯級(jí)水電站水能開(kāi)發(fā)始于2001年，至2012年已建成梯級(jí)電站8座，水電站(群)人為蓄放水對(duì)天然徑流過(guò)程的干擾增大，進(jìn)而對(duì)鶯落峽水文站徑流年內(nèi)分配、年最大洪峰流量及洪水過(guò)程漲落變化產(chǎn)生一定影響，正在興建的黃藏寺水利樞紐工程(黃藏寺水庫(kù)屬大(2)型，系年調(diào)節(jié)水庫(kù))對(duì)這些變化的影響尤為顯著。

黑河干流梯級(jí)水電站庫(kù)容調(diào)節(jié)能力不大，電站聯(lián)合調(diào)度運(yùn)行對(duì)鶯落峽水文站年徑流量變化趨勢(shì)影響不大(在黃藏寺水庫(kù)建成蓄水之前)；徑流年內(nèi)分配均勻度在2001年發(fā)生轉(zhuǎn)折，這與水電站開(kāi)始建設(shè)

圖10 鶯落峽站不同時(shí)段典型洪水過(guò)程線Fig.10 Typical flood hydrograph at Yingluoxia station in different periods

表4 洪水過(guò)程特征變化情況

時(shí)間基本吻合，2017—2020年呈顯著增加趨勢(shì)，亦是8座水電站建成投入聯(lián)合調(diào)度運(yùn)行后，改變了天然徑流過(guò)程，對(duì)徑流年內(nèi)分配產(chǎn)生顯著影響；年最大洪峰流量多年平均值由突變點(diǎn)前510 m3/s增加到突變點(diǎn)后512 m3/s，增加0.39%，增加趨勢(shì)不顯著；受人為頻繁調(diào)節(jié)閘門(mén)的影響，水位呈直上直下臺(tái)階狀變化，漲落率和變化頻次增加，洪水過(guò)程呈明顯陡漲直落形態(tài)。

4.2 建 議

(1)科學(xué)制定黃藏寺水利樞紐工程調(diào)度運(yùn)行方案，科學(xué)調(diào)度，嚴(yán)格管理，有效控制上游來(lái)水過(guò)程，保障中游區(qū)用水需求，保證下游生態(tài)關(guān)鍵期用水量。

(2)合理分配黃藏寺水庫(kù)下泄水量和過(guò)程，減少中游灌區(qū)閉口時(shí)間，既能緩解中下游用水矛盾，又能提升水庫(kù)和下游梯級(jí)電站發(fā)電效益，避免人為造峰頻繁發(fā)生，威脅下游平原(走廊)區(qū)防洪安全。