基于WEAP模型的葉爾羌河流域需水量模擬※

(新疆塔里木河流域管理局， 新疆 庫爾勒 841000)

基于WEAP模型的葉爾羌河流域需水量模擬※

周海鷹

(新疆塔里木河流域管理局， 新疆 庫爾勒 841000)

本文以葉爾羌河流域為研究對象，基于WEAP模型對流域現狀需水量進行了模擬。結果表明：干旱年，研究區總需水量為60.71億m3，缺水10.597億m3； 重度干旱年，總需水量為56.559億m3，缺水14.749億m3； 極端干旱年，總需水量為55.164億m3，缺水16.143億m3。造成缺水的主要原因是：河川徑流年內年際分配不均，缺乏山區控制性水利工程，平原水庫蒸發，滲漏損失較大及河道沿程損失大。研究結果為該流域水資源的可持續利用提供了科學依據。

水資源配置； WEAP模型； 需水量計算； 葉爾羌河流域

1 引 言

需水量的模擬與預測是掌握區域需水規律及特征的重要手段[1]，在水資源規劃、水資源系統模擬、供水風險分析等方面廣泛應用[2]。許多發達國家從20世紀60年代就開始重視對國民經濟各部門需水量的預測，我國從20世紀80年代末開始對需水量預測進行了深入研究[3]。目前常用的預測方法主要有趨勢分析法、時間序列法、定額法以及回歸分析法[4]等，為了減少不確定性，一些非線性理論，如灰色系統理論、混沌理論、人工神經網絡、小波理論等[5-6]被引入到需水量模擬中。

本文以WEAP(water evaluation and planning system) 理論為基礎，以葉爾羌河流域為研究對象，對流域現狀需水量進行模擬，以期實現流域水資源的優化調配。

2 WEAP模型介紹

WEAP模型是由瑞典斯德哥爾摩環境研究所開發的將供水、水質和生態系統保護背景之下考慮水資源開發的綜合模型[7-8]。模型以月為時間步長計算系統中“節點”和“連接”的水量和污染物的質量平衡。近年來，WEAP系統已廣泛用于流域未來水資源供需平衡評估和氣候、決策驅動的水資源管理情景分析等。

2.1 系統構成要素

WEAP模型基于水資源系統節點對水資源利用情況進行規劃和管理。它主要由12個要素組成，分別為：需求點、集水盆地、水庫、地下水、河流、分流、輸送連接、回流、流量要求、流量測站、過水河流發電、廢水處理廠，通過對它們的組合連接，實現對水資源系統的模擬[9]。

2.2 模型運算流程

2.2.1 系統定義

主要包括以下內容[10]：

a.設置時間跨度。主要包括時間范圍、時間步長和水文年起始。

b.設置空間界限。主要包括研究區域所在的地理范圍。

c.建立水資源系統概化節點網絡圖。

2.2.2 現狀基準設置

現狀基準是對現狀基準年的水資源系統數據和運行情況的準確描述，主要包括研究的第一年逐月的供給和需求數據的詳細說明。模型通過現狀基準年的詳細數據，預設模型參數。

2.2.3 情景預案設計

情景預案設計是WEAP模型的核心，WEAP模型可以生成并比較預案，以評價不同情景預案的需水、成本和環境影響。如：人口增長和經濟發展模式發生變化，需水量會發生多大變化？若大量開發利用地下水，水環境和生態環境會發生怎么樣變化？實施高效的節水灌溉技術后，生產成本、產量與需水量間關系會發生什么變化？產業結構調整和農作物構成發生變化時，土地利用方式及生態環境會如何變化？

2.3 運算法則

WEAP模型采用線性規劃方法求解能最大化滿足需求點的最優解，以月為間隔計算系統中每個節點和連接的水和污染物質量平衡，受需求優先順序、供給擇優順序、質量平衡和其他限制約束。

WEAP模型均以月為時間間隔(步長)運算，從現狀基準年的第一個月計算到預案最后一年的最后一個月，每個月都獨立于前面的一個月。因此，在一個月中，所有進入系統的水(如源頭來水、地下水補給、河道徑流)或者儲存在潛水層或水庫中、或者在月末之前離開系統(例如：河流末端的出流、需求點消耗、水庫或河段蒸發、輸送及回流損失)。由于時間尺度相對較長(月)，所有流量被假定瞬時出現。因此，一個需求點可以從河流取水，消耗一部分，將其余部分返回到廢水處理廠加以處理或返回到河流。該回流可為下游需求點在同一個月所用。

每個月計算遵循以下步驟：

a.需求點和流量要求的年需求和月供給要求。

b.集水盆地徑流和下滲，假定沒有灌溉入流。

c.系統中每個節點和連接的水的入流和出流，這包括計算從供水水源取水以滿足需求和分派水庫存水。該步驟由線性規劃(LP)求解，試圖優化需求點和河道內流量要求滿足度，受需求優先順序、供給擇優順序、質量平衡和其他限制。

d.需求點產生的污染物，污染物的量和處理，受體水體負荷，河流中濃度。

e.水力發電。

f.資本和運行成本及收入。

2.3.1 連接規則

在WEAP模型中，需求點按需求優先順序和供給擇優順序配水。

2.3.2 質量平衡約束

質量平衡方程是WEAP模型中水的月收支計算的基礎，質量平衡按式(1)計算：

∑QZI=∑QZO-∑QZX

(1)

式中QZI——月總入流,m3/s；

QZO——月總出流,m3/s；

QZX——月總消耗,m3/s。

2.3.3 滿足度

滿足度是指需求被滿足的百分比，是為每個需求點生成的一個新的線性規劃變量。WEAP模型的目標是最大化所有需求點的滿足度。在沒有足夠水量滿足優先順序相同的所有需求時，WEAP模型以其需求的相同百分比滿足所有需求。計算式為

∑QI=QXSβ

(2)

式中QI——入流量，m3/s；

QXS——需水量，m3/s；

β——滿足度，%。

3 葉爾羌河流域的WEAP模型構建

根據葉爾羌河流域內水資源調配的實際，基于流域內的水源條件和引用水關系，結合流域行政分區，將水資源系統中的主要要素、水源以及用戶進行概化，將研究區劃分為葉城、澤普、莎車、麥蓋提、巴楚、岳普湖和前海7個計算子分區，分區水源主要是河道水源、水庫水源和地下水源?？紤]到地表供水實際情況，將部分分區分開概化，得到葉爾羌河流域水資源系統概化圖，如下圖所示。

葉爾羌河流域灌區WEAP模型概化圖

葉爾羌河流域灌區共概化23個水庫、57個需求點、7個地下水和4條河流。

a.水庫：分別是保拉水庫、宗郎一二水庫(合并)、白來克其亞水庫、蘇蓋提水庫、桑水庫、米吉東水庫和墩巴克水庫(合并)、吉仁力瑪水庫、汗克爾水庫 、塔合其水庫、色里勿衣水庫、艾力西湖水庫、東方紅上下庫、邦克爾水庫、衛星水庫、紅海子水庫、草龍水庫、蘇庫恰克水庫、小海子水庫、永安壩水庫和依干其水庫。

b.需求點：總計57個需求點，分別是葉城灌區1、葉城灌區2、澤普灌區、莎西灌區、莎東灌區1、莎東灌區2、麥西灌區、麥東灌區1、麥東灌區2、麥東灌區3、岳普湖灌區、巴楚北灌區、巴楚南灌區、前進灌區、小海子灌區、葉城城鄉、葉城工業、葉城牲畜1、葉城牲畜2、葉城漁業1、葉城漁業2、澤普城鄉、澤普工業、澤普牲畜、澤普漁業、莎車城鄉、莎車工業、莎西牲畜、莎西漁業、莎東牲畜1、莎東漁業1、莎東牲畜2、莎東漁業2、麥蓋提城鄉、麥蓋提工業、麥西牲畜、麥西漁業、麥東牲畜1、麥東漁業1、麥東牲畜3、麥東漁業3、岳普湖城鄉、岳普湖工業、岳普湖牲畜、岳普湖漁業、巴楚城鄉、巴楚工業、巴楚北牲畜、巴楚北漁業、巴楚南牲畜、巴楚南漁業、農三師、前海工業、前進牲畜、前進漁業、小海子牲畜和小海子漁業。

c.地下水：7個地下水分別是葉城地下水、澤普地下水、莎車地下水、麥蓋提地下水、岳普湖地下水、巴楚地下水和前海地下水(各地下水節點用來儲存對應灌區內的水庫、河道、渠系滲漏水量)。

d.河流：4條河流分別是葉爾羌河、提孜那甫河、柯克亞河和烏魯克河。

用水關系根據水資源供給現狀進行概化。模型的輸入包括需水量、入流量和有關參數，其中需水量包括農業、城鄉、工業、牲畜、漁業需水量；入流量包括河流來水量；有關參數包括水庫庫容、河道損失、供給優先順序等。

3.1 河流來水量

河流來水量指的是在不同干旱等級下進入研究區的水量，見表1。

表1 葉爾羌河流域不同干旱等級下的月河道來水量

3.2 需求點概化

葉爾羌河流域分為葉城灌區、澤普灌區、莎車灌區、麥蓋提灌區、麥東灌區1、岳普湖灌區、巴楚灌區、前海灌區。

各灌區包含若干需求點，為方便模型運行和考慮實際情況，將莎車灌區，麥蓋提灌區，巴楚灌區和前海灌區分別分成莎西灌區、莎東灌區1、莎東灌區2，麥西灌區、麥東灌區1、麥東灌區2、麥東灌區3，巴楚北灌區、巴楚南灌區，前進灌區、小海子灌區。

4 結果及分析

4.1 基準年模擬結果

WEAP模型運行結果見表2。

由表2可知，現狀年(2012年)葉爾羌河流域總需水量為71.31億m3，通過模型對不同干旱等級水資源供需計算得：干旱年，灌區全年總供水量為60.71億m3，缺水10.597億m3，開采地下水10.22億m3(含泉水1.605億m3)；重度干旱年，灌區全年總供水量為56.559億m3，缺水14.749億m3，開采地下水10.22億m3(含泉水1.605億m3)；在極端干旱年，灌區全年總供水量為55.164億m3，缺水16.143億m3，開采地下水10.22億m3(含泉水1.605億m3)。

表2 不同干旱等級下葉爾羌河流域供需平衡模擬成果

4.2 結果分析

從水資源供需分析結果可知，研究區現狀年干旱情況下缺水較為嚴重，不同干旱等級下的缺水量分別為10.597億m3、14.749億m3和16.143億m3，分別占總需水量的14.86%、20.68%和22.64%。缺水主要集中在春季3—6月，為農業灌溉季節性缺水。

通過分析，其主要原因如下：

a.葉爾羌河徑流年內分配不均，干旱年來水總量極小，現狀又缺乏山區控制性水庫，干流和支流來水得不到有效調控，汛期水量得不到充分利用，造成春季農業灌溉嚴重缺水。

b.平原水庫蒸發、滲漏損失較大。研究區平原水庫多年水庫損失為2.3億m3，約占水庫需水總量的30%左右，成為流域內缺水的一個重要原因。

c.河道損失水量很大。在不同干旱等級下，葉爾羌河和提孜那甫河的天然來水總量分別為63.365億m3、56.638億m3和53.855億m3，但在扣除河道損失后，灌區內的供水量僅有51.758億m3、46.135億m3和44.364億m3，沿程損失水量分別達11.607億m3、10.503億m3和9.491億m3，全年河道損失水量占河道來水量的18%。

5 結 論

本文基于WEAP模型對葉爾羌河流域現狀水量配置進行了模擬，結果表明，不同干旱等級下研究區均存在缺水現象。缺水主要集中在春季3—6月，為農業灌溉季節性缺水。造成缺水的主要原因是河川徑流年內、年際分配不勻，缺乏山區控制性水利工程，平原水庫蒸發、滲漏損失較大及河道沿程損失大。

[1] 彭飛,劉兵.基于生態的新疆干旱灌區水資源的優化配置[J].石河子大學學報(自然科學版),2017(1):291-297.

[2] 姜志嬌,楊軍耀,任興華.基于“三條紅線”及SE-DEA模型的水資源優化配置[J].節水灌溉,2016(11):81-84.

[3] 李承紅,姜卉芳,何英.基于WRMM模型的水資源配置及方案優選研究[J].水資源與水工程學報,2016(3):32-38.

[4] 熊雪珍,何新玥,陳星.基于改進TOPSIS法的水資源配置方案評價[J].水資源保護,2016(2):14-20.

[5] 吳丹,王士東.基于需求導向的城市水資源優化配置模型[J].干旱區資源與環境,2016(2):31-37.

[6] 孟凡臣.基于MIKEBASIN模型的再生水灌區水資源配置研究[J].灌溉排水學報,2014(6):10-13.

[7] 王帥兵.基于WEAP模型的董志塬水資源合理利用研究[J].干旱地區農業研究,2015(4):238-245.

[8] 王鵬全,方三朝.基于WEAP模型的石羊河水資源安全分析[J].中國農村水利水電,2014(9):92-94.

[9] 徐慧,張運超.基于WEAP的塔里木河干流灌溉滿足度模擬[J].干旱區研究,2013(4):582-587.

[10] 宋瑞勇,劉青勇.基于WEAP模型的龍口市水資源優化配置研究[J].中國農村水利水電,2011(4):30-33.

WaterdemandsimulationofYarkantRiverbasinbasedonWEAPmodel※

ZHOU Haiying

(XinjiangTarimRiverBasinAdministration,Korla841000,China)

In the paper, Yarkant River basin is adopted as a research object for simulating river basin current (2012) water demand on the basis of WEAP model. The results show that total water demand in the study area is 6.071 billion m3during drought year, and water deficiency is 1.0597 billion m3. Total water demand was 5.6559 billion m3in the severe drought year, and water deficiency is 1.4749 billion m3. Total water demand is 5.516.4 billion m3in extreme drought year, and water deficiency is 1.6143 billion m3. The main reason for the water shortage is shown as follows: uneven interannual distribution in the river runoff year, no controllable water conservancy project in mountainous area, evaporation of plain reservoir, great leakage loss and large loss along the river. The research results provide scientific basis for sustainable utilization of water resources in the river basin.

water resources allocation; WEAP model; calculation of water demand; Yarkant River basin

10.16616/j.cnki.10-1326/TV.2017.12.004

水利部公益行業科研專項經費(201501059)資助

TV214

A

2096-0131(2017)12-0012-05