劉 慧，楊澤川，許鳳冉，張學(xué)君，白音包力皋，穆祥鵬，胡宏昌

（1.中國(guó)水利水電科學(xué)研究院 水力學(xué)研究所，北京 100038；2.水利部國(guó)際經(jīng)濟(jì)技術(shù)合作交流中心，北京 100038；3.中國(guó)水利水電科學(xué)研究院 防洪抗旱減災(zāi)研究所，北京 100038；4.清華大學(xué) 水利系，北京 100084）

1 研究背景

干旱是影響最大的自然災(zāi)害之一，往往對(duì)農(nóng)業(yè)、漁業(yè)、供水、環(huán)境等方面造成巨大影響[1-4]。瀾滄江-湄公河流域是連接?xùn)|南亞六國(guó)的重要紐帶，該地區(qū)以農(nóng)業(yè)和漁業(yè)為主，對(duì)降水和湄公河依賴度高[5]。近年該流域在干季（12月至次年5月）干旱頻發(fā)（如2004/2005、2009/2010、2012/2013、2015/2016、2019/2020年等）[6]，給各國(guó)造成了巨大損失[7-10]。

和地球上其他大型流域相比，瀾滄江-湄公河的水文過(guò)程更接近天然過(guò)程[11]。但隨著經(jīng)濟(jì)社會(huì)發(fā)展，該流域的水利工程建設(shè)正在迅速增加。瀾滄江干流“兩庫(kù)七級(jí)”梯級(jí)水庫(kù)在2010年至2015年間建成，其中第一座多年調(diào)節(jié)水庫(kù)——小灣水庫(kù)于2010年汛期開始發(fā)揮調(diào)蓄作用[8]。湄公河干流的第一個(gè)水電站——沙耶武里水電站于2019年建成，湄公河支流截止2009年共有26座大型和中型水庫(kù)，但規(guī)劃和在建大中型水庫(kù)達(dá)99座[5]。

關(guān)于水利工程對(duì)瀾滄江-湄公河水文過(guò)程的影響，已有一些研究基于歷史數(shù)據(jù)或未來(lái)氣候變化情景數(shù)據(jù)進(jìn)行定量評(píng)估，然而由于瀾滄江-湄公河地區(qū)自然地理?xiàng)l件復(fù)雜、觀測(cè)資料和水利工程調(diào)度資料短缺，所采用的數(shù)據(jù)、模型和假設(shè)不同，雖然對(duì)于變化趨勢(shì)（削峰補(bǔ)枯）的結(jié)論一致，但結(jié)果卻有較大的差異[8，12-16]。如湄公河委員會(huì)（MRC）和中國(guó)水利部研究[8]認(rèn)為，受水利工程影響2010-2015年間瀾滄江徑流量平均每年減少4.9%，而氣候變化使年徑流量平均減少20%。HAN等[16]的研究認(rèn)為人類活動(dòng)對(duì)瀾滄江2008—2014年的徑流變化貢獻(xiàn)率約為95%。另外，氣候變化的不確定性也進(jìn)一步增加了水利工程影響評(píng)估的不確定性。Lauri等[17]的研究表明，采用不同的GCM數(shù)據(jù)驅(qū)動(dòng)模型可能導(dǎo)致不同的預(yù)估時(shí)段（2032—2042）和基準(zhǔn)期（1982—1992）的徑流變化，雨季為-11%～+15%，干季為-10%～+13%。

根據(jù)水循環(huán)不同環(huán)節(jié)的水分虧缺，干旱分為氣象干旱、水文干旱和農(nóng)業(yè)干旱。當(dāng)前瀾滄江-湄公河的干旱研究仍多集中于氣象或農(nóng)業(yè)干旱分析，如基于模型或再分析降水?dāng)?shù)據(jù)重建歷史干旱事件、基于有限的地面站觀測(cè)分析干旱的時(shí)空特征[18-21]。Zhang等[18]基于VIC模擬瀾滄江-湄公河土壤水分，并建立對(duì)應(yīng)的干旱指標(biāo)，較好的重建了2004/2005和2009/2010年干旱事件，并能夠?qū)υ摰貐^(qū)的農(nóng)業(yè)干旱進(jìn)行實(shí)時(shí)評(píng)估。Tang和Cao[19]基于瀾滄江流域的8個(gè)氣象站數(shù)據(jù)，計(jì)算SPI指標(biāo)，發(fā)現(xiàn)瀾滄江下游區(qū)域發(fā)生洪旱的頻率更高，東南亞季風(fēng)和青藏高原季風(fēng)是影響該地區(qū)洪水和干旱的關(guān)鍵因素。Guo等[20]的研究表明，湄公河南部和北部較中部更容易發(fā)生干旱，發(fā)生在三角洲地區(qū)的干旱往往具有時(shí)間長(zhǎng)強(qiáng)度高的特點(diǎn)。周婷等[21]利用1977—2010年的38個(gè)氣象站降水資料，采用SPI分析瀾滄江-湄公河流域的氣象干旱時(shí)空分布特征以及2009/2010年干旱。目前對(duì)流域水文干旱及其與氣象干旱之間關(guān)系的相關(guān)研究較少。

本文基于氣象干旱指標(biāo)和水文干旱指標(biāo)，分析該流域干旱演變趨勢(shì)；選擇2010年前后的兩場(chǎng)典型流域氣象干旱（2009/2010和2012/2013）作為研究對(duì)象，分析干旱特征，以及氣象干旱和水文干旱之間的聯(lián)系；分析瀾滄江梯級(jí)水庫(kù)的補(bǔ)水作用。相關(guān)工作有利于科學(xué)認(rèn)識(shí)近來(lái)發(fā)生的全流域干旱，而且有利于流域內(nèi)國(guó)家從全流域視角調(diào)整合作戰(zhàn)略，采取合理的措施共同應(yīng)對(duì)旱情，減少災(zāi)害損失，保障流域經(jīng)濟(jì)社會(huì)發(fā)展。

2 研究區(qū)域概況

瀾滄江-湄公河流域及其干流主要水文站位置如圖1所示，水文站信息如表1所示。景洪水電站是瀾滄江最末一級(jí)水電站，因此選擇位于景洪水電站下游的景洪水文站作為瀾滄江流域的代表水文站開展相關(guān)分析?；诳色@取的數(shù)據(jù)，在湄公河干流選擇了7個(gè)水文站開展分析，并分別選取清盛（Chiang Saen）、穆達(dá)漢（Mukdahan）和上?。⊿tung Treng）站作為湄公河干流上游、中游和下游的代表站。圖1中不同顏色的區(qū)域代表對(duì)應(yīng)水文站的集水面積/區(qū)間匯水面積。

瀾滄江-湄公河流域在2012年10月至2013年4月期間遭遇了氣象干旱災(zāi)害，干旱的范圍和強(qiáng)度與2009/2010年干季發(fā)生的干旱相當(dāng)。而不同的是，中國(guó)小灣水庫(kù)在2010年汛期開始具備調(diào)蓄能力，并在2012年11月至2013年4月的干季期間對(duì)下游進(jìn)行了“補(bǔ)水”。由于兩場(chǎng)干旱發(fā)生在臨近年份，本文認(rèn)為其自然地理和社會(huì)經(jīng)濟(jì)變化可以忽略，即流域內(nèi)用水等因素基本不變，氣象干旱和水文干旱關(guān)系的演變主要是氣候因素和水利工程引起的。

圖1 水文站及其集水區(qū)域

表1 瀾滄江-湄公河干流主要水文站信息

3 數(shù)據(jù)和方法

3.1 數(shù)據(jù)降水?dāng)?shù)據(jù)采用全球陸面數(shù)據(jù)同化產(chǎn)品（Global Land Data Assimilation System，GLDAS）。GLDAS數(shù)據(jù)是基于全球觀測(cè)數(shù)據(jù)以及四個(gè)陸面過(guò)程模型模擬結(jié)果，采用數(shù)據(jù)同化方法得到覆蓋全球的高精度（2.5°到1km）的數(shù)據(jù)集。研究中采用的GLDAS降水序列長(zhǎng)度為1948—2014年，空間分辨率為0.25°×0.25°，范圍為瀾滄江-湄公河流域。

水文數(shù)據(jù)包括清盛、瑯勃拉邦、廊開、那空帕儂、穆達(dá)漢、巴色和上丁等湄公河干流的主要水文站點(diǎn)的1985—2016年逐日流量數(shù)據(jù)。

3.2 研究方法從氣象、水文角度出發(fā)，計(jì)算SPI和SRI干旱指標(biāo)，作為氣象干旱和水文干旱的監(jiān)測(cè)診斷指標(biāo)，并利用Mann-Kendall趨勢(shì)檢驗(yàn)方法進(jìn)行干旱趨勢(shì)分析。進(jìn)一步選取大型水利工程建設(shè)前后的兩場(chǎng)干旱（2009/2010年和2012/2013年）進(jìn)行干旱特征對(duì)比分析。

SPI是氣象干旱評(píng)估的常用指標(biāo)[21-22]，被推薦應(yīng)用于所有國(guó)家的氣象干旱指標(biāo)，其定義和計(jì)算過(guò)程參見國(guó)際氣象組織（WMO）相關(guān)文件[23]。表2為國(guó)家標(biāo)準(zhǔn)《氣象干旱等級(jí)》（GB/T 20481-2017）對(duì)SPI值的干旱等級(jí)劃分。SRI（或SSI）是參考SPI計(jì)算原理提出的標(biāo)準(zhǔn)化徑流指數(shù)，該指數(shù)基于長(zhǎng)序列月徑流計(jì)算，用于月尺度以上的水文干旱評(píng)估[24-26]，其干旱等級(jí)劃分同SPI。

當(dāng)SPI的統(tǒng)計(jì)尺度達(dá)到6個(gè)月時(shí)，其信息往往可以和河道水文變化情況聯(lián)系起來(lái)[23]。Gibson等[24]用6個(gè)月尺度的SPI和對(duì)應(yīng)尺度的徑流指標(biāo)開展流域干旱分析。因此，本文分析湄公河干流主要水文站集水面積上的氣象干旱和水文干旱聯(lián)系時(shí)，采用6個(gè)月尺度的SPI（SPI6）和6個(gè)月尺度的SRI（SRI6）。在展示氣象干旱年際波動(dòng)時(shí)，采用12個(gè)月尺度的SPI（SPI12），SPI12更能反映大范圍長(zhǎng)歷時(shí)的嚴(yán)重干旱，若非連續(xù)12個(gè)月時(shí)段內(nèi)顯著濕潤(rùn)或干燥，該指標(biāo)會(huì)接近0。

Mann-Kendall檢驗(yàn)法是降水、干旱趨勢(shì)檢測(cè)的常用方法，該方法可判斷變化趨勢(shì)是否顯著，以及是否發(fā)生突變，并確定突變發(fā)生的時(shí)間[27-30]。具體計(jì)算過(guò)程可參考相關(guān)文獻(xiàn)，本文中采用的顯著性水平為0.05。

表2 基于SPI的干旱等級(jí)劃分

4 結(jié)果及分析

4.1 年際變化特征瀾滄江-湄公河流域的SPI12序列如圖2所示?？傮w而言瀾滄江-湄公河流域的SPI無(wú)顯著變化趨勢(shì)，呈洪旱交替變化，基于SPI12指標(biāo)分析得到，1948—2014年間瀾滄江-湄公河流域達(dá)到中旱等級(jí)的年份包括：1968、1992、1971、1958、1993、2010、1967和1955年（按照SPI遞增順序）。本研究的結(jié)果和Guo等[20]的研究結(jié)論基本一致，其中1992—1994年是1981—2014年間干旱最嚴(yán)重的時(shí)期，SPI12小于-1。圖3展示了SPI序列的Mann-Kendall趨勢(shì)檢驗(yàn)結(jié)果為顯著的網(wǎng)格及其變化趨勢(shì)，從圖3中可知流域內(nèi)局部有SPI增大（降水增加）的趨勢(shì)，包括瀾滄江源頭、瑯勃拉邦周邊和巴色至上丁區(qū)間，這些區(qū)域在干季、雨季和全年尺度上的趨勢(shì)相似。

瀾滄江-湄公河干流主要水文站集水面積上的SRI趨勢(shì)分析結(jié)果如表3所示。從表3可知，瀾滄江-湄公河全流域年尺度SRI變化不顯著。干季SRI（徑流）總體呈增加趨勢(shì)；雨季各區(qū)間變化趨勢(shì)不同，下游上丁站SRI略有下降趨勢(shì)。

4.2 典型干旱特征瀾滄江-湄公河干流水文站集水面積2009/2010、2012/2013年干季的SPI6計(jì)算結(jié)果以及對(duì)應(yīng)SRI6結(jié)果如圖4所示。2009/2010年和2012/2013年干季所有水文站對(duì)應(yīng)集水面積的SPI6均小于0（即降水偏少），整個(gè)流域平均SPI6值分別為-0.38和-0.43，即2012/2013年旱情總體比2009/2010年嚴(yán)重。從圖4中可看出2009/2010年干季比2012/2013年干季干旱面積大。但是，2012/2013年出現(xiàn)干旱的地區(qū)其強(qiáng)度達(dá)到了中旱程度，而2009/2010年為輕旱。對(duì)比圖4中2009/2010年和2012/2013年干季的SRI可知，2009/2010年湄公河干流主要水文站均發(fā)生了水文干旱，而2012/2013均未出現(xiàn)水文干旱。

圖2 瀾滄江-湄公河流域12個(gè)月尺度的SPI序列

圖3 具有顯著變化趨勢(shì)的SPI變化趨勢(shì)空間分布（1948—2014）

表3 1986—2016年SRI變化趨勢(shì)

由圖4可知，2009/2010年的SPI6和SRI6表現(xiàn)出較好的一致性，各水文站集水面積的氣象干旱多為輕旱，水文干旱大多為輕旱或中旱；而2012/2013年的SPI6和SRI6表現(xiàn)差異顯著，氣象干旱為無(wú)旱或中旱，水文條件為多年平均狀態(tài)或偏豐狀態(tài)。這說(shuō)明在2012/2013年干季上游實(shí)際來(lái)水接近或大于天然徑流量，本文認(rèn)為湄公河上游段的上述現(xiàn)象主要源于瀾滄江水庫(kù)發(fā)揮的“補(bǔ)枯”作用。而那空帕農(nóng)站以下河段的SRI指標(biāo)增大，可能來(lái)自湄公河中下游產(chǎn)流能力較強(qiáng)的支流的徑流補(bǔ)給或其他水利工程對(duì)徑流的調(diào)節(jié)。正如有關(guān)研究指出，湄公河干流的水文過(guò)程變化不只來(lái)自瀾滄江梯級(jí)水庫(kù)，湄公河支流建設(shè)的水庫(kù)也有影響[14]。如湄公河支流Nam Theun河上的跨流域引水發(fā)電項(xiàng)目Nam Theun 2水電站對(duì)湄公河干流穆達(dá)漢水文站水情的影響，泰國(guó)境內(nèi)Ubol Ratana、Sirindhorn和Lam Pao3個(gè)水庫(kù)對(duì)湄公河干流巴色水文站水情的影響等，這些影響還有待進(jìn)一步研究。

圖4 2009/2010、2012/2013年干季的SPI6空間分布以及各水文站的SRI6

圖5 2009/2010年、2012/2013年和多年平均干季徑流量

4.3 瀾滄江梯級(jí)水庫(kù)的補(bǔ)水作用根據(jù)瀾滄江-湄公河干流水文站點(diǎn)1960—2009年的多年平均月流量[8]和2009/2010年、2012/2013年干季逐月流量資料，計(jì)算得到其干季徑流量，如圖5所示。結(jié)果顯示，2012/2013年干季，景洪水文站干季徑流量由多年平均118.2億m3增大至169.0億m3，增加了50.8億m3，增幅達(dá)43%；2012/2013年干季徑流量比2009/2010年增加67.0億m3，增幅66%。清盛水文站2012/2013年干季徑流量由多年平均177.9億m3增大至231.5億m3，增加了53.6億m3，增幅30%；2012/2013年干季徑流量比2009/2010年增加58.9億m3，增幅34%。景洪站的多年平均（1960—2009）干季徑流量占清盛站的比例為66%，2012/2013年干季增大到73%。2012/2013年瀾滄江干流梯級(jí)電站發(fā)揮的調(diào)豐補(bǔ)枯作用使湄公河干流，尤其是清盛站的旱季徑流量較多年平均增加顯著。

根據(jù)景洪站和清盛站1960—2009年和2010—2015年的流量資料，統(tǒng)計(jì)分析了這兩個(gè)水文站干季徑流量占年徑流量百分比在瀾滄江干流梯級(jí)水庫(kù)運(yùn)行前后的變化情況。結(jié)果顯示，瀾滄江干流梯級(jí)水庫(kù)運(yùn)行后，景洪水文站干季徑流量由天然情況下的118.2億m3（1960—2009年間干季徑流量占年徑流量的21%）增大至177.7億m3（2010—2015年干季徑流量占年徑流量的41%），增加了59.5億m3（占比增加20%）。清盛站干季徑流量占年徑流量的比例由21%增大到33%。可以看出，瀾滄江干流梯級(jí)水庫(kù)不僅在2012/2013年干季發(fā)揮了補(bǔ)水作用，而且在其運(yùn)行后總體上增加了湄公河干流干季徑流量。考慮到干季是湄公河流域農(nóng)業(yè)用水高峰期[31]，瀾滄江干流梯級(jí)水庫(kù)的運(yùn)行對(duì)于緩解湄公河干季的旱情以及保障全流域可持續(xù)發(fā)展方面具有積極作用。

5 結(jié)論和建議

（1）瀾滄江-湄公河流域年降水趨勢(shì)性變化不顯著。1948—2014年間瀾滄江-湄公河流域達(dá)到中旱等級(jí)的年份包括：1968、1992、1971、1958、1993、2010、1967和1955年。SPI變化趨勢(shì)顯著的地區(qū)包括瀾滄江源頭、瑯勃拉邦周邊和巴色至上丁區(qū)間，這些區(qū)域的SPI呈增大趨勢(shì)（降水增多）。

（2）瀾滄江-湄公河年徑流量變化趨勢(shì)不顯著。干季SRI（流量）總體呈增加趨勢(shì)；雨季各區(qū)間變化趨勢(shì)不同，下游上丁站SRI（流量）略有下降趨勢(shì)。

（3）瀾滄江-湄公河流域在2009/2010和2012/2013年干季經(jīng)歷了兩場(chǎng)氣象干旱，其中2012/2013年總體較前者嚴(yán)重。2009/2010年湄公河干流沿程均發(fā)生了水文干旱，而2012/2013年各站均未發(fā)生水文干旱。鑒于兩個(gè)氣象干旱年份臨近，流域內(nèi)其他變化較小，湄公河干流上游段出現(xiàn)上述現(xiàn)象應(yīng)主要源于瀾滄江水庫(kù)發(fā)揮的“補(bǔ)枯”作用。

（4）2012/2013年干季，瀾滄江梯級(jí)水庫(kù)發(fā)揮徑流調(diào)節(jié)作用，景洪水庫(kù)下泄流量比多年平均徑流量增加50.8億m3，增幅達(dá)43%；2012/2013年干季比2009/2010年同期徑流量增加67.0億m3，增幅66%。清盛水文站2012/2013年干季徑流量比多年平均徑流量增加53.6億m3，增幅30%；同時(shí)2012/2013年干季徑流量比2009/2010年增加58.9億m3，增幅34%。即2012/2013年瀾滄江干流梯級(jí)電站發(fā)揮的調(diào)豐補(bǔ)枯作用使湄公河干流，尤其是清盛站的旱季徑流量較多年平均增加顯著。多年數(shù)據(jù)分析結(jié)果表明，瀾滄江干流梯級(jí)水庫(kù)運(yùn)行后總體上增加了湄公河干流干季徑流量，對(duì)于緩解湄公河干季旱情和保障流域可持續(xù)發(fā)展方面有積極作用。

（5）瀾滄江-湄公河流域各國(guó)應(yīng)從全流域視角加強(qiáng)水資源開發(fā)利用協(xié)作，制定科學(xué)的流域規(guī)劃，充分發(fā)揮流域內(nèi)水利工程的抗旱功能。同時(shí)，應(yīng)意識(shí)到瀾滄江水庫(kù)調(diào)度只對(duì)干流水文過(guò)程有影響，無(wú)法緩解支流的大面積旱情，各國(guó)應(yīng)挖掘自身潛力，提高基礎(chǔ)設(shè)施水平，同時(shí)加強(qiáng)水利工程聯(lián)合調(diào)度。