衡邵干旱走廊典型灌區水資源優化調度研究

劉元沛　盛　東　山紅翠　鄒　亮　徐幸儀

（湖南省水利水電科學研究所　　長沙市　410007）

衡邵干旱走廊典型灌區水資源優化調度研究

劉元沛盛東山紅翠鄒亮徐幸儀

（湖南省水利水電科學研究所長沙市410007）

年內干旱是干旱研究的重要方向，對農業灌溉具有非常重要的意義。由于地區不同，降水不同，農業種植結構不同，導致各個地區的年內干旱情勢也各不相同。選取南方豐水地區的衡邵干旱走廊作為研究對象，通過對豐水年、平水年和枯水年進行水庫群優化調度研究，揭示衡邵干旱走廊的干旱成因、干旱時間和短缺水量。結果發現：衡邵干旱走廊干旱成因是農業干旱，未來可能會出現社會經濟干旱；干旱時間為每年7～9月，符合該走廊夏秋連旱的特點；水量短缺以龍溪橋水庫最為顯著，豐水年約為1960萬m3，平水年約為2383萬m3，枯水年約為4302萬m3。因此，衡邵干旱走廊應積極發展高效節水灌溉、開展渠道防滲、適當調整農業種植結構和增加水源工程以應對未來的社會經濟干旱。

年內干旱衡邵干旱走廊干旱成因

近年來，我國已在干旱研究方面取得諸多成果。但是，當前干旱研究以年際為主，年內研究相對較少。干旱研究方法主要有基于標準化降水指數（SPI）［1］、馬爾科夫鏈模型（Markov）［2］、自回歸移動平均模型（ARIMA）［3］、帕默爾干旱指數（PDI/PDSI）［4］、相對濕潤度指數（CMI）［5］、降水距平百分率（Pa）［6］、土壤相對濕度干旱指數［7］、綜合氣象干旱指數（CI）［8］、水文模型［9］等方法，通過干旱指數對干旱歷史進行評估，研究干旱的成因和規律，預測未來干旱情勢。年內干旱則受地區的用水情況影響較大，這一情況在南方豐水地區顯得尤為明顯，再加上最嚴格水資源管理制度對用水總量的限制，年內干旱情勢顯得愈發嚴峻。

1 研究區概況

衡邵干旱走廊是湖南省四大干旱區之一，地處湖南省中部，雪峰山脈以南，五嶺山脈以北，是以湘江、資江流域的分水嶺為中軸線分布的廣大丘陵地區，主要包括湘江中上游、資江中上游的丘陵地區，位于北緯25°40′～27°31′，東經110°25′～113°16′之間。

本次研究選擇衡邵干旱走廊內衡南縣龍溪橋灌區為研究對象。龍溪橋灌區多年平均降雨量為1263mm，較全省多年平均降雨量少187mm。灌區包括龍溪橋灌片、唐灣灌片、敏東灌片和大宜沖灌片共4個部分。龍溪橋灌區內有龍溪橋水庫、唐灣水庫、敏東水庫和大宜沖水庫，有引水工程自歐陽海灌區、唐灣水庫、敏東水庫和大宜沖水庫引水至龍溪橋水庫，具體如表1所示。

表1　龍溪橋灌區各水庫基本特征表

2 理論和方法

本文通過灌溉期（4～10月）的降水量數據劃分豐水年、平水年和枯水年；根據灌溉定額和灌溉水有效利用系數核算灌溉需水量；根據月降雨～月徑流關系曲線計算產水模數，再根據集雨面積計算產水量；最后建立基于線性規劃的水庫群優化調度模型，核算缺水量，進而確定年內干旱期情勢。

2.1劃分豐平枯

龍溪橋灌區農業種植以水稻為主，早稻種植自4月份開始7月份結束，晚稻種植自7月份開始10月份結束。龍溪橋灌區季節性干旱嚴重，尤其是在灌溉時期（4～10月）。本文采用1966～2008年衡南縣雨量站日降雨量數據為依據，在此基礎上計算全年降雨量數據和灌溉期降雨量數據，多年月平均降水量如圖1所示。

圖1　龍溪橋水庫多年月平均降水量

根據灌溉期降雨量與多年平均灌溉期降雨量之間的差值采用水文年法將年降雨量情況分為枯水年、平水年、偏豐水年和豐水年。距平百分率P=（某年灌溉期降雨量-多年平均灌溉期降雨量）/多年平均灌溉期降雨量×100%，以10%為劃分標準，具體如表2所示。

表2　豐平枯劃分標

2.2確定典型年

根據統計資料，龍溪橋灌區農業種植以單、雙季稻為主，種植比例為0.15∶0.85。灌溉季節為每年的4～10月。龍溪橋灌區按照0.48的灌溉水有效利用系數核算年均毛需水量為8785萬m3，遠大于龍溪橋灌區的可用水量。同時，龍溪橋灌區多年平均降雨量為1263mm遠小于全省平均1450mm。因此，采用灌溉期的豐平枯年份作為水庫調度的年份以研究不同情況下的水庫調度方案。按照龍溪橋灌區現有資料和豐平枯年份的平均情況，選擇1987年作為灌溉期為枯水年的調度年份，1982年作為灌溉期為豐水年的調度年份，1983年作為灌溉期為平水年的調度年份。經皮爾遜-III型曲線計算得到：P1982=7.17%，P1983= 49.76%，P1983=83.54%，如圖2所示。

圖2　皮爾遜-ⅠⅠ型曲線擬合圖

2.3灌溉需水量計算

根據湖南省地方標準《用水定額》DB43/T-2014的90%供水保證率的灌溉定額和龍溪橋灌區灌水情況確定灌溉綜合定額標準如表3所示，其中旱地灌水定額采取已建工程的的具體情況確定為200m3/畝。

表3　綜合定額計算表m3/畝

按照全區0.48的灌溉水有效利用系數核算龍溪橋灌片的灌溉毛需水量，結果如表4所示。同樣的方法計算唐灣灌片、敏東灌片和大宜沖灌片的灌溉毛需水量。

表4　灌溉需水量計算表m3/畝，萬m3

2.4產水量核算

利用龍溪橋水庫現有資料，查移用井頭江月降雨-月徑流深關系曲線，確定龍溪橋灌區的產水模數。龍溪橋水庫集雨面積43.3km2，3處引水工程集雨面積25.65km2，考慮引水工程本身灌溉及入庫水量損失，故按面積的70%計算入庫水量，具體產水量如表5所示。

表5　1982、1983、1987年產水量計算表萬m3

2.5水庫群優化調度

水庫調度采用Lingo軟件編程實現基于線性規劃的水庫群聯合調度。水庫群聯合調度是對龍溪橋水庫、大宜沖水庫、敏東水庫、唐灣水庫開展聯合調度，考慮大宜沖水庫、敏東水庫和唐灣水庫的可用水量引入龍溪橋水庫進行調度。水庫群聯合調度主要遵循以下原則：

（1）防洪約束：6～7月份是汛期，因此水庫水位在6～7月應低于防洪限制水位：

V（t）＜Vf（1）

式中V（t）——t月份的庫容；

Vf——對于水庫的防洪限制水位對于的庫容。

（2）水量平衡約束：當月的可用水量=前月的可用水量+當月產水量+當月引水量-當月灌溉用水量，考慮到引水渠道破壞嚴重，滲漏損失率取50%，式（2）為龍溪橋灌溉渠水量平衡約束，式（3）為引水的水量平衡約束，式（3）為唐灣水庫的水量平衡約束（敏東水庫和大宜沖水庫亦如此）。

Vl（t）=Vl（t-1）+Vlc（t）+Vy（t）-Vlg（t）（2）

Vy（t）=0.5Vty（t）+0.5Vmy（t）+0.5Vdy（t）+Voy（t）（3）

Vt（t）=Vt（t-1）+Vtc（t）-Vtg（t）-Vty（t）（4）

式中

Vl（t）——龍溪橋水庫t月份可用水量；

Vl（t-1）——龍溪橋水庫t-1月份可用水量；

Vlc（t）——龍溪橋水庫集水區t月份產水量；

Vy（t）——龍溪橋水庫t月份引水量；

Vlg（t）——龍溪橋水庫t月份灌溉可用水量；

Vty（t）——t月份從唐灣水庫的引水量；

Vmy（t）——t月份從敏東水庫的引水量；

Vdy（t）——t月份從大宜沖水庫的引水量。

（3）灌溉約束：以龍溪橋水庫為例，當月的灌溉可用水量與當月歐陽海引水量之和小于等于當月灌溉需水量，如式（5）所示；當月自龍溪橋水庫的灌溉需水量＜當月龍溪橋水庫可用水量，對于大宜沖水庫、敏東水庫、唐灣水庫則不含歐陽海引水量，如式（5）、式（6）所示。

Vlg（t）+Voy（t）≤Vlx（t）（5）

Vlx（t）≤Vl（t）（6）

式中Vlg（t）——t月份龍溪橋灌片的灌溉可用水量；

Voy（t）——t月份龍溪橋灌片的自歐陽海灌片的引水量；

Vlx（t）——t月份龍溪橋灌片毛需水量；

Vl（t）——t月份龍溪橋水庫的可用水量。

（4）引水約束：龍溪橋水庫從唐灣水庫、敏東水庫和大宜沖水庫引水的前提是唐灣水庫、敏東水庫和大宜沖水庫的水量能夠滿足自身的需求。以唐灣水庫為例，敏東水庫和大宜沖水庫類似，如式（7）所示。

式中Vtg（t）——唐灣水庫t月份可用水量；

Vtx（t）——唐灣灌片t月份灌溉毛需水量；Vty（t）——唐灣水庫t月份引水量；

Vt（t-1）——唐灣水庫t-1月份引水量；

Vtg（t）——唐灣水庫t月份灌溉可用水量；Vt（t）——唐灣水庫t月份可用水量。

（5）目標函數：以龍溪橋灌片、唐灣灌片、敏東灌片和大宜沖灌片的灌溉可用水量最大為準。

式中Vlg（t）——t月份龍溪橋水庫的可用灌溉水量；

Vtg（t）——t月份唐灣水庫的可用灌溉水量；

Vmg（t）——t月份敏東水庫的可用灌溉水量；

Vdg（t）——t月份大宜沖水庫的可用灌溉水量；

Voy（t）——t月份自歐陽海灌區的引水量。

表6　典型年灌溉滿足情況　　萬m3

3 結果分析

根據水庫群優化調度的結果，對4個灌片的灌溉可用水量與灌溉需水量的比值作為灌溉滿足度來分析缺水情況，缺水情況如表6所示。從表6中可以看出不論是枯水年、平水年還是豐水年，龍溪橋灌片均有大額的水量差，灌溉無法得到保證；唐灣水庫無論豐水年、平水年還是枯水年灌片灌溉均能得到保證，而敏東水庫和大宜沖水庫在豐水年和平水年灌溉渠灌溉能夠得到保證，枯水年灌溉無法得到保證。

從4～10月份調度結果可以看出龍溪橋灌區干旱期為7月、8月和9月。從1968～2008的43年降雨量數據可以得到月平均降雨量為113mm，而7月、8月的月平均降雨量為107mm，118mm，和43年的月平均降雨量基本持平，9月的月平均降雨量為58mm小于7月、8月和9月。造成7月、8月和9月干旱的原因是7月份開始種植晚稻，需水量，總需水量大約占整個灌溉季節灌溉需水量的60%，如圖3、圖4、圖5、圖6所示。

考慮到未來灌區基礎設施不斷完善，灌溉水有效利用系數將進一步提高，因此對灌溉水有效利用系數提高后的情況進一步優化調度，具體如表7所示，結果表明灌溉水有效每提高0.01，水量滿足情況可提高約1.25%，但是水量仍然有很大的缺口不能滿足灌溉需要。為了滿足灌溉需水，一方面要根據灌溉水量缺口，新建或者擴建水源工程、調水工程；另一方面，要嘗試調整農業種植結構，可以大幅度節省農業用水，例如玉米的灌溉定額約為水稻灌溉定額的1/3。

圖3　溪橋水庫4～10月灌溉情況

圖4　唐灣水庫4～10月灌溉情況

圖5　敏東水庫4～10月灌溉情況

圖6　大宜沖水庫4～10月灌溉情況

表7　提高灌溉水有效利用系數后滿足度

4 結　論

通過對衡邵干旱走廊龍溪橋灌區的基于典型年的水庫群優化調度分析可以得到以下結論：衡邵干旱走廊的干旱期為7月、8月和9月；衡邵干旱走廊地區屬于相對豐水地區，降水量大，造成缺水的主要原因是蓄水工程不足和當月降水量不能滿足當月灌溉需水量造成的；衡邵干旱走廊的干旱成因屬于人類活動造成的農業干旱不屬于氣象干旱和水文干旱，未來隨著經濟的發展可能出現社會經濟干旱。對于衡邵干旱走廊而言未來隨著經濟的發展，可能會出現工業用水擠占農業用水情況，因此，應發展高效節水灌溉、適當調整農業種植結構減少對水量的需要；同時增加水源工程，提升蓄水能力，以應對未來的社會經濟干旱。

［1］楊斌斌，孫萬光.改進POA算法在流域防洪優化調度中的應用［J］.水電能源科學，2010，（12）：36-38+115.

［2］林潔，夏軍，佘敦先，張利平.基于馬爾科夫鏈模型的湖北省干旱短期預測［J］.水電能源科學，2015，（04）：6-9+51.

［3］許文寧，王鵬新，韓萍，嚴泰來，張樹譽.Kappa系數在干旱預測模型精度評價中的應用——以關中平原的干旱預測為例［J］.自然災害學報，2011，（06）：81-86.

［4］劉庚山，郭安紅，安順清，劉巍巍.帕默爾干旱指標及其應用研究進展［J］.自然災害學報，2004，13（4）：21-27.

［5］王明田，王翔，黃晚華，張玉芳，馬均.基于相對濕潤度指數的西南地區季節性干旱時空分布特征［J］.農業工程學報，2013，28（19）：85-92.

［6］王富強，王雷.基于降水距平百分率的河南省干旱特征分析［J］.中國農村水利水電，2014，（12）：84-88.

［7］王素萍，張存杰，宋連春，李耀輝，馮建英，王勁松.多尺度氣象干旱與土壤相對濕度的關系研究 ［J］.冰川凍土，2013，35（4）：865-873.

［8］鄒旭愷，任國玉，張強.基于綜合氣象干旱指數的中國干旱變化趨勢研究［J］.氣候與環境研究，2010，15（04）：371-378.

［9］熊見紅.長沙市農業干旱規律分析及旱情預報模型探討［J］.湖南水利水電，2003，（4）：29-31.

劉元沛（1989-），男，碩士，主要從事水資源管理、3S技術等相關工作，E-mail：530964293@qq.com；盛東（1979-），男，博士，高級工程師，主要從事水資源管理工作，E-mail：shengdong1979@163.com。

（2016-06-17）