楊 連 海

(甘肅省甘州區(qū)水務(wù)局， 甘肅 張掖 734000)

甘州區(qū)位于甘肅省河西走廊中部、古“絲綢之路”南北兩線和“居延古道”交匯點上，位于東經(jīng)100°6′～100°52′，北緯30°32′～39°24′，最高海拔2 230 m，最低海拔1 410 m，平均海拔1 474 m，屬典型的溫帶大陸性氣候，具有寒冷、干燥、晝夜溫差大、日照時間長、降水稀少變率大、蒸發(fā)大等明顯的水平分帶差異[1]。甘州區(qū)年平均氣溫為7.8 ℃，1月平均氣溫最低，為-9.1 ℃,7月平均氣溫最高，為22.3 ℃，全年極端最低氣溫為-28.2 ℃，極端最高氣溫為39.8 ℃；多年平均降水量為132.6 mm，降水年際變化大，年內(nèi)分配不均；多年平均蒸發(fā)量為1 796.7 mm，年均蒸發(fā)量是降水量的14倍，是全國重點缺水區(qū)域，也是中國西部地區(qū)自然生態(tài)環(huán)境比較脆弱的地區(qū)。2016年甘州區(qū)實際可供水量為7.79 億m3，本地區(qū)的可利用水資源100%利用，農(nóng)業(yè)用水比例85.3%，用水效率低。現(xiàn)狀年全區(qū)人口51.29 萬人，耕地面積約8.92 萬hm2。

1 需水量預(yù)測

1.1 預(yù)測方法

總需水量包括生活需水、牲畜需水、工業(yè)需水、農(nóng)業(yè)需水和生態(tài)需水。區(qū)域需水量預(yù)測模型為：

WZ=WP+WB+WI+WA+WG

(1)

式中：WZ為總需水量；WP為生活需水量；WB為牲畜需水量；WI為工業(yè)需水量；WA為農(nóng)業(yè)需水量；WG為生態(tài)需水量。

(1)生活需水預(yù)測。 生活需水量包括城鎮(zhèn)生活需水和農(nóng)村生活需水。城鎮(zhèn)生活需水分為城鎮(zhèn)居民生活需水和城鎮(zhèn)公共需水；農(nóng)村生活需水包括農(nóng)村居民生活需水和牲畜需水[2]。生活需水在一定范圍內(nèi)，其增長速度是比較有規(guī)律的，因而可以用綜合分析定額方法推求未來需水量。此方法考慮的因素是用水人口需水定額。用水人口以計劃部門預(yù)測數(shù)為基礎(chǔ)。需水定額以現(xiàn)狀用水調(diào)查數(shù)據(jù)為基礎(chǔ)。分析歷年變化情況，考慮不同水平年居民生水平的改善及提高程度，擬定其相應(yīng)的用水定額：

(2)

式中：Pk為k地區(qū)人口數(shù)量，萬人；QPk為k地區(qū)生活綜合用水定額，L/(人·d)。

(2)牲畜需水預(yù)測。牲畜需水量包括大牲畜需水量和小牲畜需水量。牲畜數(shù)量以計劃部門預(yù)測數(shù)為基礎(chǔ)，需水定額以現(xiàn)狀用水調(diào)查數(shù)據(jù)為基礎(chǔ)。分析歷年變化情況，考慮不同水平年牲畜數(shù)量的變化，擬定其相應(yīng)的用水定額：

(3)

式中：Nk為k地區(qū)牲畜的數(shù)量，萬頭；QBk為k地區(qū)牲畜的用水定額，L/(頭·d)。

(3)工業(yè)需水預(yù)測。工業(yè)需水量包括火電工業(yè)需水，國有及規(guī)模以上工業(yè)需水及一般工業(yè)需水3類。工業(yè)需水主要利用萬元增加值需水量法進行預(yù)測：

(4)

式中：QIk為k地區(qū)萬元工業(yè)增加值用水量，m3/萬元；Ek為k地區(qū)工業(yè)增加值，萬元。

(4)農(nóng)業(yè)需水預(yù)測。農(nóng)業(yè)灌溉需水預(yù)測采用定額法。農(nóng)業(yè)灌溉需水預(yù)測分鄉(xiāng)鎮(zhèn)進行，各鄉(xiāng)鎮(zhèn)需水量之和即為全區(qū)農(nóng)業(yè)灌溉需水量。灌溉需水量預(yù)測涉及3個關(guān)鍵指標(biāo)：各類作物的凈灌溉定額、灌溉水利用系數(shù)和灌溉面積。定額法的計算公式為：

(5)

式中：QAk為k地區(qū)凈灌溉定額，m3/hm2；Ak為k地區(qū)灌溉面積，hm2；ηk為k地區(qū)灌溉水利用系數(shù)。

(5)生態(tài)需水預(yù)測。生態(tài)環(huán)境需水與農(nóng)業(yè)需水類似，仍采用定額法，用不同種類植物及相應(yīng)定額計算。計算公式：

(6)

1.2 預(yù)測結(jié)果

以2016年為現(xiàn)狀年，預(yù)測年為2025年和2030年。需水量預(yù)測結(jié)果如表1所示。

由表1可知，2025年甘州區(qū)總需水量為75 351.70 萬m3，其中生活需水量為1 879.56 萬m3，牲畜需水量為1 639.21萬m3，工業(yè)需水量為1 480 萬m3，農(nóng)業(yè)灌溉需水量為65 909 萬m3，生態(tài)需水量為4 444 萬m3。2030年甘州區(qū)總需水量為78 408.27 萬m3，其中生活需水量為2 114.31 萬m3，牲畜需水量為2 091.96 萬m3，工業(yè)需水量為1 379 萬m3，農(nóng)業(yè)灌溉需水量為68 188 萬m3，生態(tài)需水量為4 635 萬m3。

表1 甘州區(qū)2025、2030年預(yù)測需水量 萬m3Tab.1 Forecasted water demand summary Tables of 2025 and 2030 in Ganzhou District

1.3 供需平衡分析

根據(jù)甘肅省水資源管理“三條紅線”要求，實際分配給甘州區(qū)2016年、2025年、2030年可供水總量分別為7.79、6.92、7.02 億m3；2016年、2025年、2030年甘州區(qū)總需水量分別為8.42、7.54和7.84 億m3，供水明顯不能滿足區(qū)域需水要求，分別缺水0.63、0.62、0.82 億m3，缺水率分別達到7.48%、8.22%、10.46%。甘州區(qū)發(fā)展道路嚴(yán)重受阻，社會經(jīng)濟發(fā)展將受到嚴(yán)重制約，在沒有外來水源的情況下，水資源緊缺的局面短期內(nèi)將難以緩解，供需矛盾將難以消除?？偭靠刂浦笜?biāo)按“以供定需”的原則確定[3]：當(dāng)需水預(yù)測量大于區(qū)域可利用水資源總量時，即過度開發(fā)，采用區(qū)域可利用水資源總量作為分配總量；當(dāng)需水預(yù)測量小于區(qū)域可利用水資源總量時，即預(yù)留部分水量，采用區(qū)域預(yù)測需水總量作為分配總量；當(dāng)需水預(yù)測量等于區(qū)域可利用水資源總量時，即飽和狀態(tài)下，通過多次供需分析采用區(qū)域預(yù)測需水總量或可利用總量作為分配總量。甘州區(qū)預(yù)測的需水量均大于實際可利用的水資源總量，所以按照“以供定需”的原則，采用可供用水資源總量作為分配總量。

2 水資源合理配置

2.1 模型構(gòu)建

水資源優(yōu)化配置是以區(qū)域社會經(jīng)濟發(fā)展與水環(huán)境相協(xié)調(diào)作為優(yōu)化目標(biāo)，在環(huán)境經(jīng)濟學(xué)的基礎(chǔ)上，依據(jù)經(jīng)濟發(fā)展水平條件，對區(qū)域有限水資源進行科學(xué)合理的優(yōu)化配置，實現(xiàn)水資源的合理髙效開發(fā)利用[4]。模型的建立首先要確定目標(biāo)函數(shù)，代表區(qū)域可持續(xù)發(fā)展過程中各方面的主要發(fā)展指標(biāo)，用于衡量水資源優(yōu)化配置的方案優(yōu)劣，其決策變量由經(jīng)濟系統(tǒng)和水資源系統(tǒng)內(nèi)的諸多約束條件組成。對于甘州區(qū)水資源優(yōu)化配置是多目標(biāo)的，以經(jīng)濟效益、環(huán)境效益、社會效益等綜合效益最佳為優(yōu)化配置目標(biāo)[5]。甘州區(qū)水資源在進行優(yōu)化配置中，從環(huán)境生態(tài)的角度來進行水資源配置的綜合考慮，其優(yōu)化目標(biāo)的函數(shù)表達為。

區(qū)域水資源優(yōu)化配置的目標(biāo)是經(jīng)濟、環(huán)境、社會綜合效益最大,并保持系統(tǒng)協(xié)調(diào)持續(xù)發(fā)展。其一般模型形式為：

Z=opt {f1(X),f2(X),f3(X)}

G(X)≤0X≥0

式中：X為決策向量；f1(X)、f2(X)、f3(X)分別為經(jīng)濟效益、環(huán)境效益和社會效益;G(X)為約束條件集,表示水資源承載力、環(huán)境容量、上地資源、其他社會約束和子系統(tǒng)狀態(tài)方程等。

(1)目標(biāo)1:區(qū)域供水凈效益最大。

(7)

(2)目標(biāo)2：區(qū)域重要污染物排放量最小。

(8)

(9)

2.2 約束條件

(3)用戶需水能力約束。

(10)

(4)污染物約束。

(11)

(5)變量非負(fù)約束。

改革開放近40年來，我國得到了快速發(fā)展，但農(nóng)業(yè)農(nóng)村的發(fā)展相對滯后。當(dāng)前處于全面建成小康社會的決勝階段，更為迫切的是要同步推進農(nóng)業(yè)現(xiàn)代化，這就要求必須把農(nóng)業(yè)農(nóng)村放在優(yōu)先發(fā)展的戰(zhàn)略地位，解決城鄉(xiāng)發(fā)展不平衡的問題，縮小城鄉(xiāng)差距，實現(xiàn)與城鎮(zhèn)的同步發(fā)展。

(12)

2.3 模型的求解

粒子群優(yōu)化算法(Particle Swarm Optimization，PSO)是由Kennedy和Eberhart于1995年提出的一種新型的群智能進化算法,它可以靈活方便地處理具有大量等式約束、不等式約束和同時包含連續(xù)變量、離散變量的混合整數(shù)優(yōu)化問題，因此,對于水資源合理配置問題,采用粒子群優(yōu)化算法也是一種可行方案,為水資源合理配置提供了一種很有前景和潛力的新型方法。算法描述如下：

PSO算法一般是采用下面的公式對粒子的位置和速度進行更新：

(13)

(14)

式中：k為迭代代數(shù)；學(xué)習(xí)因子c1、c2是2個正常數(shù),一般取值為[1.5，2.05]；r1、r2是均勻分布于[0，1]的2個隨機數(shù)。

公式(13)主要通過3部分來計算粒子i新的速度：粒子i前一時刻的速度、粒子i當(dāng)前位置與自己最好位置之間的距離、粒子i當(dāng)前位置與群體最好位置之間的距離。粒子i通過公式(14)計算新位置的坐標(biāo)。

由于不同問題對局部搜索和全局搜索能力要求不同,所以在式(13)中引入慣性權(quán)重因子來平衡算法全局搜索和局部搜索,得到新的速度更新公式：

從式(15)中可以看出,慣性權(quán)重ω決定粒子先前速度對當(dāng)前速度的影響程度,具有平衡算法全局搜索和局部搜索的能力。在求解配水優(yōu)化問題時,與遺傳算法相比,粒子群優(yōu)化算法具有較好的收斂性和收斂速度,能以較少的迭代次數(shù)收斂到與遺傳算法用較多迭代次數(shù)才能得到的解甚至更好的解,算法靈活簡單,易于處理不可微、非線性、非凸的優(yōu)化間題。雖然其初衷主要用于求解連續(xù)變量問題，但經(jīng)過簡單處理，PSO算法也能夠很好地求解離散變量以及包含連續(xù)變量和離散變量在內(nèi)的混合整數(shù)規(guī)劃問題。

本研究采用粒子群算法(PSO)求解數(shù)學(xué)模型，粒子群算法是基于迭代的優(yōu)化算法，期望找到優(yōu)化問題中的一個潛在的解，實際中可以理解為是為了在空間中找到粒子(particle)的位置。具體流程見圖1。

圖1 粒子群算法流程Fig.1 Flow chart of particle swarm optimization

2.4 甘州區(qū)水資源合理配置結(jié)果

針對甘州區(qū)自然環(huán)境特點和經(jīng)濟生活發(fā)展?fàn)顩r確定水資源配置時的優(yōu)先序為：先充分滿足城鄉(xiāng)居民及牲畜等生活用水，保障穩(wěn)定生態(tài)環(huán)境發(fā)展的生態(tài)用水，滿足重點工業(yè)用水及一般工業(yè)的生產(chǎn)用水，公平保障農(nóng)民生存的農(nóng)業(yè)基本用水，協(xié)調(diào)分配各水平年可供用水。經(jīng)過對甘州區(qū)2025年近期水平年和2030年遠(yuǎn)景水平年的社會經(jīng)濟指標(biāo)預(yù)測及用水定額和效率變化趨勢的分析，分別預(yù)測各水平年各區(qū)域和各行業(yè)的用水需求得：2025年甘州區(qū)各行業(yè)總需水量為7.54 億m3，其中農(nóng)業(yè)6.59 億m3、工業(yè)0.15 億m3、生活0.19 億m3、生態(tài)0.44 億m3；2030年張掖各行業(yè)總需水量為7.84 億m3，其中農(nóng)業(yè)用水6.82 億m3、工業(yè)0.14 億m3、生活0.21 億m3、生態(tài)0.46 億m3。

對甘州區(qū)水資源平衡分析后的水量配置如表2所示。甘州區(qū)水資源匱乏，供水量小于需水量。經(jīng)過粒子群優(yōu)化算法的水資源合理配置后，2025、2030年缺水量比配置前分別降低了516.7、1 194.3 萬m3。通過合理配置當(dāng)?shù)厮Y源，已經(jīng)大大提高了當(dāng)?shù)厮Y源利用率，促進當(dāng)?shù)亟?jīng)濟的發(fā)展。

3 結(jié) 論

結(jié)合甘州區(qū)的用水現(xiàn)狀，采用粒子群優(yōu)化算法對水資源進行合理優(yōu)化配置，并對規(guī)劃水平年(2025年和2030年)的水資源優(yōu)化配置結(jié)果進行了分析，對甘州區(qū)水資源合理配置提出了以下建議和對策。

(1)當(dāng)?shù)亟邓纬傻乃Y源時空分布不均，年際變化較大，蒸發(fā)量大?，F(xiàn)有水利工程調(diào)蓄能力不足，汛期洪水資源得不到充分利用。應(yīng)該加強以豐補枯的措施，提髙地表水利用率。

(2)針對甘州區(qū)水資源短缺、農(nóng)業(yè)用水較多的實際情況，開展農(nóng)業(yè)節(jié)水灌溉和高效用水是十分必要的，為其他產(chǎn)業(yè)的發(fā)展提供富余水資源。

(3)在制定資源開發(fā)利用與治理保護總體實施方案時，做到總體目標(biāo)、任務(wù)與措施相協(xié)調(diào)，建設(shè)規(guī)模與發(fā)展機制和生產(chǎn)力發(fā)展水平相協(xié)調(diào)，因地制宜，大中小工程相結(jié)合，努力提高用水效率，合理利用地表水與地下水資源，有效保護水資源，積極治理利用廢污水、雨水等其他水源。

表2 甘州區(qū)2025、2030年水資源合理配置成果 萬m3Tab.2 Water resources rational allocation results of 2025 and 2030 in Ganzhou District