石羊河流域生態治理目標下節水灌溉技術分區

王軍德 金彥兆 程玉菲

摘要：針對目前實施的石羊河流域綜合治理工程提出的蔡旗斷面下泄水量目標，以石羊河流域中游涼州區為研究區域，基于GIS技術，選擇土壤類型、土地利用類型、坡度和高程等因子，利用生態景觀分析中的生態適宜性原理確定節水灌溉技術的適宜性分區；基于最小阻力模型，以灌區渠系為源區，利用距離函數確定節水灌溉的阻力表面特征，基于自然分類法將涼州區耕地劃分為高、中、低三種強度的節水灌溉區及常規灌溉區。

關鍵詞：節水灌溉技術；生態目標；生態適宜性；最小累計阻力模型；空間布局；石羊河流域

中圖分類號：5274 文獻標識碼：A 文獻標志碼：A doi：10.3969/j.issn.1000-1379.2018.09.035

近年來，隨著我國水利事業的快速發展，節水灌溉技術得到了快速發展。據《2015年全國水利發展統計公報》[1]，到2015年底，我國工程節水灌溉面積已達3106萬hm2，其中噴灌、微灌面積901.20萬hm2、低壓管灌面積891.20萬hm2。節水灌溉技術的發展一方面有效緩解了水資源供需矛盾，增加了農業灌溉面積和農業水資源效率，為促進地區社會經濟發展和生態環境的協調起到了積極作用；另一方面節水灌溉技術的發展對干旱半干旱區的綠洲內部生態特別是人工生態產生了一定程度的負面影響，大量農田路網的防護林因灌區渠系襯砌、降低滲漏而死亡，從而導致綠洲農田生態系統遭受干熱風等自然災害的影響，出現作物耗水量增大甚至死亡的現象，綠洲內部生態穩定性遭到一定程度的破壞[2-5]。如何科學地結合節水灌溉技術特點在灌區內部合理布置各項節水灌溉技術，確定節水灌溉技術布置的空間分布方式和強度，從而確保綠洲內部生態穩定、支撐農業發展，是節水灌溉技術空間布局需要考慮的重要內容，也是節水灌溉技術優化研究的熱點問題之一。

近年來，專家學者對節水灌溉技術空間布局開展了大量研究，但研究方向主要集中于綜合考慮地區水土資源和農作物需水特征進行節水灌溉的空間優化布局研究[6-8]，忽略了區域綠洲內部生態穩定問題，對節水灌溉技術的水源支撐條件未作考慮，因此有必要探索一種在強調區域農業水源供給與生態環境良性發展條件下的節水灌溉技術空間優化方法。位于我國西北干旱內陸區的石羊河流域從2007年開始實施重點流域治理項目，核心就是推廣高效節水灌溉技術，從而在流域中上游減少農業用水，實現流域治理中控制蔡旗斷面有充足的水資源，從而保障流域下游農業、生態和地下水的穩定性恢復。本研究以該流域中游的涼州區為例，以構建生態目標良性恢復和農業節水優化布局為出發點，利用生態適宜性評價原理和最小累計阻力模型進行基于生態目標的節水灌溉空間分區，為進一步研究節水灌溉優化提供技術支撐。

1 研究區域

涼州區位于甘肅省河西走廊東段，地處石羊河流域中游，總土地面積4873.85km2。涼州區屬溫帶大陸性干旱氣候區，具有太陽輻射強、干旱少雨、日照充足、蒸發量大、溫差大的特點。全年日照時數為2967.8h，日照百分比為67%；日均氣溫7月最高，為29℃，1月最低，為-14.9℃。多年平均降水量160mm，年蒸發能力2020mm。流經區內的西營河、金塔河、雜木河、黃羊河是全區的主要地表水源，多年（1993-2003年）平均地表年徑流量為6.80億m3。區內地下水年取水量為5.25億m3，其中90%為河渠入滲量、基巖裂隙水入滲量，近l0 a來地下水年均可利用量為2.82億m3。全區地表水、地下水可利用總量多年平均為9.62億m3。

2 數據來源與研究方法

2.1 數據來源

研究區的數字地形資料來源于USGS/NASASRTM數據，土壤和土地利用數據分別來自于中國科學院南京土壤研究所和中國科學院西部環境數據中心。研究區的灌區、水系、河渠和節水灌溉等資料來源于涼州區水務局。

2.2 研究方法

2.2.1 節水灌溉適宜性評價原理

節水灌溉適宜性是指節水灌溉技術及其應用所依存的水源、耕地、地形等條件的適宜與否及適宜程度。這種適宜性是基于自然過程而不致引起節水灌溉技術發展受阻，按照土地適宜性評價結果構建的整體格局，是符合區域水資源條件、土壤條件和氣候條件的，能夠滿足技術的地區適宜性要求。

基于生態目標的節水灌溉空間優化布局，其基本思路是根據干旱區流域上中下游的水力聯系及水資源與生態之間的關系，節水灌溉技術的布局應設置不同的節水程度，從而減小節水灌溉對區域地下水和地表水轉換關系的影響，根據高程、土壤類型、坡度、土地利用類型等生態因子，確立維護節水灌溉和生態良性發展雙贏條件下的節水灌溉布局原則，利用生態適宜性原理確定節水灌溉技術的空間布局。

2.2.2 最小累計阻力模型

最小累計阻力模型是1992年由Knaapen提出的，在生態景觀規劃和布局研究中應用廣泛，后經俞孔堅等[9-12]的修改，可用下式表示：

MCR=fmin∑DijRi（1）

依據原公式參數定義，結合本研究目的，將式（1）中參數定義如下MCR為節水灌溉最小累計阻力值；Dij為從節水灌溉技術源地j到應用單元i的空間距離；Ri為應用單元i對節水灌溉技術擴張的阻力系數；∑表示節水灌溉技術應用單元i與節水灌溉技術源地j之間穿越所有單元的距離和阻力的合計；min表示被評價的節水灌溉技術應用張力對于不同源地取累計阻力最小值；f為單調遞增函數，表示最小累計阻力與節水灌溉技術應用的正相關關系。

最小累計阻力值反映了節水灌溉技術應用的潛在可能性及趨勢，通過最小累計阻力的大小可判斷應用單元與源地的連通性和相似性。

本研究將渠道作為節水灌溉的源地，主要考慮以下兩方面：一是石羊河流域中游武威盆地和下游的民勤盆地地下水單元相對獨立，流域中下游的水力聯系主要為地表水，流域內部灌區主要采用渠道進行水量輸送；二是渠道將地表水和地下水輸送到田間，耗費了大量的資金技術成本，更應該采用節水灌溉技術提升水資源利用效率。因此，本研究借助最小累計阻力模型模擬基于生態目標的節水灌溉技術從渠道向周圍空間擴展的過程，最小累計阻力小的地區適宜開展節水灌溉，而最小累計阻力大的地區相對而言不適宜開展節水灌溉。

3 基于生態目標的節水灌概技術空間優化

3.1 節水灌溉技術生態適宜性評價

3.1.1 評價因子分級

各種灌溉技術都有一定的適用范圍，在選擇時主要應考慮作物、地形、土壤和水源等條件[13]。畦灌、溝灌、淹灌和漫灌等地面灌溉方法適用于地形平坦、水源充足和中等滲透性土壤的區域，而噴灌適用于復雜地形、水量較少和各種滲透性土壤區域，滲灌、滴灌和微噴灌等高效節水灌溉適用于水量缺乏、地形平坦的區域，除滲灌僅適用于弱滲透性土壤外，滴灌和微噴灌適用于各種滲透性土壤。因此，本研究選擇土地利用類型、坡度和土壤類型作為節水灌溉適宜性的評價因子?？紤]到西北干旱區獨特的水文循環特征，流域中下游的自然生態對維護人工綠洲的穩定具有舉足輕重的作用，流域中下游區域應采用傳統灌溉而非強化節水技術措施，從而維護天然生態，因此選擇高程作為評價基于生態目標的節水灌溉技術適應性的重要因子。

3.1.2 適宜性評價原則、模型及過程

基于上述評價因子，以實現石羊河流域中游蔡旗斷面下泄水量生態目標下節水灌溉優化布局為基本目標，適宜性評價的基本原則如下：①由于高海拔地區降水資源豐富，其天然生態靠降水得以維護，節水灌溉技術的發展對當地自然生態的影響程度較小，因此設定區域內高程越高，發展節水灌溉的優先級越高；②受西北干旱區風力較大的限制，噴灌技術的推廣較為緩慢，在坡度因子中僅考慮地面灌溉和高效節水灌溉的方式，因此設定坡度越小，其發展節水灌溉的優先級越高；③受資料的限制（無研究區的種植結構資料），本研究僅將土地利用類型作為節水灌溉發展的評價因子，耕地、林地和草地具有較高的發展節水灌溉的優先級，這三類面積比重越大，優先級越高；④土壤滲透性較大的土壤類型具有較高的發展節水灌溉技術的優先級，同時灌區內部發展節水灌溉的優先級高于灌區外部。

根據篩選的節水灌溉技術發展影響因素，構建以高程、坡度、土地利用類型和土壤類型為因子的節水灌溉技術適宜性評價模型為

Yi=Hi（Pi+Ti+Si）（i=1，2，…，n）（2）式中：Yi為第i個評價單元的適宜性等級分值，即阻力值；Hi為第i個評價單元的高程因子優先級；Pi、Ti、Si分別為第i個評價單元的坡度、土地利用類型、土壤類型的優先級；n為評價單元的個數。

參考其他專家學者的評價過程[11，14]，基于生態目標的節水灌溉適宜性評價的基本過程如下：基于節水灌溉技術評價的基本原則，對高程、坡度、土地利用類型和土壤類型因子進行分級，確定各個因子的優先級（見表1至表4）；利用ArcGIS空間分析功能，根據節水灌溉技術適宜性評價模型基本結構，進行節水灌溉技術阻力表面分析，即適宜性評價分析，結果見圖1。

3.2 節水灌溉技術空間格局分區

3.2.1 確立源地

節水灌溉技術源地的確立參考生態空間源地的確立方法，西北干旱內陸區農田水利建設特別是渠系建設對農業生產起著十分重要的作用，水源條件較好的地區應優先發展節水灌溉。以研究區現有的農田水利工程為主要參考因素，考慮石羊河流域關井壓田工程及研究區上游均為井灌區的特點，將上游雜木灌區、黃羊灌區、金塔灌區和西營灌區內的渠系作為節水灌溉技術空間格局生態優化的源地。這些源地為灌區農田提供水資源，是灌區農業生產和糧食安全的命脈。

3.2.2 基于生態目標的節水灌溉技術空間格局分區

基于確立的源地和適宜性評價中建立的阻力表面圖（圖1），利用ArcGIS10.1空間分析手段，采用最小累計阻力模型（MCR），生成研究區發展節水灌溉技術最小累計阻力表面圖（圖2）。

以最小累計阻力表面圖為基礎，利用ArcGIS I0.I統計分析工具和自然分類方法確定變化拐點，結合研究區耕地分布圖，利用空間分析功能生成基于生態目標的耕地節水灌溉技術分區，見圖3（圖中白色區域為沙漠、戈壁、居民用地、裸巖裸土、水域等非耕地區）。

4 結論

（1）傳統的節水灌溉分區僅僅從作物需水量和水土資源的優化配置角度加以劃分，沒有考慮區域土壤、坡度和高程等因子對節水灌溉技術適應性的影響，生態適宜性評價和最小累計阻力模型等生態學方法為解決這一問題提供了一種可能。本研究選擇石羊河流域中游的涼州區作為研究區域，基于上述生態學方法對該區節水灌溉技術進行分區，是對該方法的一種嘗試。

（2）在考慮石羊河流域生態目標的前提下，將高程作為西北干旱區內陸河流域節水灌溉技術發展的生態目標的控制因子，這是出于減弱節水灌溉技術對流域中下游地下水和地表水循環的影響角度出發的，有利于改善中下游生態環境質量和保持中下游人工綠洲的穩定。通過中上游節水灌溉技術的發展，節約中上游的地表水，從而保障下游可獲得的地表水資源量，在下游農田灌溉和生態補水過程中補充下游地下水資源，減緩地下水位漏斗擴大的趨勢。

（3）基于生態適宜性評價和最小累計阻力模型將研究區域劃分為4個分區：高強度節水灌溉區、中強度節水灌溉區、低強度節水灌溉區和傳統灌溉區。對于不同強度區域節水灌溉方式的選擇及其節水定量化指標的確定，還有待進一步分析研究。

參考文獻：

[1]中華人民共和國水利部.2015年全國水利發展統計公報[M].北京：中國水利水電出版社，2016：8，

[2]張新民，沈冰，金彥兆.節水灌溉與河西綠洲可持續發展[J].中國農村水利水電，2000（12）：9-11.

[3]王根緒，程國棟，徐中民.中國西北干旱區水資源利用及其生態環境問題[J].自然資源學報，1999，14（2）：109-116.

[4]王根緒，程國棟，沈永平.近50年來河西走廊區域生態環境變化特征與綜合防治對策[J].自然資源學報，2002，17（1）：78-86.

[5]齊學斌，龐鴻賓.節水灌溉的環境效應研究現狀與研究重點[J].農業工程學報，2000，16（4）：37-40.

[6]康紹忠，賀正中，張學.陜西省作物需水量及分區灌溉模式[M].北京：水利電力出版社，1992；239-267.

[7]李英能，黃修橋，吳景社，等.水土資源評價與節水灌溉規劃[M].北京：中國水利水電出版社，1998：86-121.

[8]吳景社，康紹忠，王景雷，等.基于主成分分析和模糊聚類方法的全國節水灌溉分區研究[J].農業工程學報，2004，20（4）：64-68.

[9]俞孔堅.生物保護的景觀生態安全格局[J].生態學報，1999，19（1）：8-15.

[10]陳利頂，傅伯杰，徐建英，等.基于源—匯生態過程的景觀格局識別方法：景觀空間負荷對比指數[J].生態學報，2003，23（11）：2406-2413.

[11]黎曉亞，馬克明，傅伯杰，等.區域生態安全格局：設計原則與方法[J].生態學報，2004，24（5）：1055-1062.

[12]張惠遠，萬軍.GIS支持下的山地景觀生態優化途徑[J].水土保持研究，1999，6（4）：69-75.

[13]郭元裕.農田水利學[M].北京：中國水利水電出版社，2002：56-57.

[14]趙筱青，王海波，楊樹華，等.基于GIS支持下的土地資源空間格局生態優化[J].生態學報，2009，29（9）：4892-4901.