青海省大峽灌區典型地塊作物耗水系數研究

呂文星　周鴻文　王永峰等

摘要：采用引排差法和最大蒸發量法對青海省大峽灌區典型地塊耗水系數進行計算研究。在典型地塊設置了2個斗渠引水口和6個退水口，并對其春灌期流量進行監測。同時在典型地塊內設置5眼地下水監測井，用于計算地下退水量。在綜合考慮地表退水和地下退水基礎上，采用引排差法計算耗水系數；同時利用最大蒸發量法計算耗水系數，對引排差法進行驗證。結果表明，采用引排差法計算得到的大峽灌區典型地塊春灌期耗水系數為0.48，采用最大蒸發量法，計算得到耗水系數為0.49，兩種方法計算得出的耗水系數非常接近，說明在試驗區耗水系數計算上，最大蒸發量法可以對引排差法進行較好的驗證。

關鍵詞：耗水系數；引排差法；最大蒸發量法；彭曼公式；青海省大峽灌區

中圖分類號：S157.1，X592 文獻標識碼：A 文章編號：0439-8114（2015）19-4692-06

DOI：10.14088/j.cnki.issn0439-8114.2015.19.009

Abstract： Diversion and drainage method and the biggest evaporation method were used respectively to calculate crop water consumption coefficient of typical plot of Daxia irrigation area in Qinghai province. Two lateral canal irrigating gates and six waterspouts were set up in the typical plot and water flow was monitored during the spring irrigating period， and five groundwater wells were also set up to monitor groundwater. In consideration of the surface water and underground water drainage， diversion and drainage method was used respectively to calculate crop water consumption coefficient. The biggest evaporation method was also used to calculate crop water consumption coefficient， to verify the diversion and drainage method. The results showed that crop water consumption coefficient calculated by diversion and drainage method and the biggest evaporation method were 0.48 and 0.49 respectively. The results were so closely， which indicated that the biggest evaporation method could lead better validation to diversion and drainage method for crop water consumption coefficient calculation in the scope of the study area.

Key words： water consumption coefficient； diversion and drainage method； the biggest evaporation method； Penman-Monteith formula； Daxia irrigation area in Qinghai province

水資源作為重要的戰略資源，是經濟社會發展的基礎支撐，也是現代農業建設的主要約束條件。青海省水資源豐富，多年平均水資源總量為629.3億m3，但水資源總量的93%流向中下游地區，利用率低；區域降水稀少，時空分配不均，黃河一級支流湟水流域水資源總量只占全省的3.5%，卻集中了全省52.3%的耕地面積，平均水資源量7 245 m3/hm2，為全國平均值的33.6%，水資源開發利用程度較低和水資源供需矛盾突出成為制約湟水流域農業發展最主要的因素[1]。灌溉方式粗放、灌區農業用水效率偏低、水利基礎設施建設滯后成為影響農業發展的薄弱環節。

目前研究主要從4個方面進行耗水量的核定[2]：①黃河干流入境年徑流量減出境年徑流量加區間匯入量，稱河段差法；②直接從黃河取水的渠道統計年引水量之和，減各排水溝年排水量，再減灌區以潛流形式直接排入黃河的排水量，稱引排差法；③用彭曼公式計算灌區蒸發蒸騰量，換算灌區耗水量，稱最大蒸發量法；④灌區作物生長期耗水來自灌溉、有效降水和作物根系吸收地下水量，3 項和稱作物總需水量法，計算耗用黃河水量。

本研究采用引排差法和最大蒸發量法進行耗水系數計算。其中引排差法是通過對所有引水干渠全面監測和對主要排水溝（大約占70%～90%的排水量）進行駐測、巡測，對于漏測部分利用鄰近排水溝排水模數進行插補計算，從而能計算出控制區域內的引水量、耗水量和排水量。最大蒸發量法采用彭曼公式計算作物蒸騰蒸發量（也就是植株蒸騰量和棵間蒸發量的總和），再換算成灌區耗水量[3]。FAO Penman-Monteith方法被認為是干旱地區和濕潤地區計算效果最好的一種方法[4]，被認為是參考作物蒸發蒸騰量（ET0）計算的標準化方法和ET0實測資料缺乏地區評價其他方法的標準[5-7]，目前已被列入了中國現行規范[8]，在中國針對不同地區開展的研究方法的評價中，基本上都是采用FAO56-PM方法的計算結果作為評價其他方法的標準[9-11]。在此背景下，對青海省湟水流域大峽灌區典型地塊采用引排差法和最大蒸發量法開展耗水系數研究，對加強區域水資源管理，區域水權轉換、提高水資源利用效率都具有十分重要的意義。

1 研究區概況

試驗布設在青海省黃河支流湟水河流域大峽灌區，地處東經102°15′～102°24′和北緯36°13′～36°30′，海拔范圍在1 650～4 400 m。大峽灌區位于湟水左岸的高店鎮河灘寨村，水源引自湟水，下游有引勝溝等湟水一級支流作為補充水源。灌區貫穿于湟水左岸的高店、雨潤、共和、碾伯、高廟5個鄉鎮的43個行政村和單位。

渠首引水樞紐位于樂都縣高店鎮河灘寨村，干渠渠道全長57 km（于2005年立項維修27 km）。干渠有各類建筑物298座，其中渡槽36座，長1 950 m，隧洞50座，總長19 200 m，倒虹吸1座，長384 m，退水17座，涵洞17座。其他建筑物168座（完好117座），其中斗門137處。農渠退水口多達198處。由于水污染加重和徑流量年內分配不均勻，每年4～6月供需矛盾極為突出，嚴重影響灌溉。

根據地形和海拔高度，該區域地貌類型分為：河谷平原川水區、黃土淺山丘陵區和石質高山腦山區3種。大峽灌區位于河谷平原川水區，該區沿湟水干流及其一級支流呈帶狀分布，由河灘和1～5級階地坡洪積扇組成，土體構型較好，質地松，是全縣的主要產糧區。

灌區屬半干旱的高原大陸性氣候類型，寒冷和干旱是區內主要氣候特點。年平均氣溫4.5～7.5 ℃，無霜期130～150 d，區域地形復雜，海拔高差大，各地降水量不盡一致，山區一般大于川區，腦山大于淺山，川水地區年降雨量為320～340 mm，蒸發量川區大于山區，川區年蒸發量達843 mm，最大凍土深度為86 cm。大峽灌區土壤主要包括灰鈣土和栗鈣土兩種，成土母質有沖積物、洪積物和次生黃土等，土質松散，質地均一，耕性好，結構呈團粒狀或粒狀。在大峽灌區引退水條件較好處設置典型地塊進行試驗觀測（圖1），典型地塊地處N36°29′16.4″和E102°13′34.8″，海拔1 950 m，平均坡度3°，種植作物為大蒜，土壤類型為灰鈣土。

2 材料與方法

2.1 試驗設計

2.1.1 引退水監測斷面選取 在大峽灌區典型地塊開展引退水量監測，典型地塊灌溉面積19.33 hm2，進水口斷面2個，退水口斷面6個，引、退水口設置見圖1和表1。

2.1.2 引退水測流方法及頻次 根據典型地塊主要作物大蒜的生長周期和耗水規律、灌溉制度、氣象條件等因素綜合分析，在3月10-25日春灌期進行監測，共16 d。按照規范要求，3月19、20日對水尺零高進行測量。3月21日之前，每日測流2次，21日通過調節閘門增加測次，完成了斗渠水位流量關系率定，采用水位-流量關系曲線法推求引水量。典型地塊水位采用駐測的方式觀測。每日按2段制觀測水位，并根據引水口斗門開啟變化情況隨時增加斷面的水位觀測次數。同時加強與灌區區段管理人員聯系，及時了解閘門開啟情況。大峽灌區典型地塊退水斷面由于退水渠沒有正規渠道，斷面不規整，且退水時間不固定，退水量小，設立水尺觀測水位難度大，故不監測水位。大峽灌區典型地塊引水口斷面采用率定水位流量關系曲線法推流。典型地塊斗渠水位采用駐測方式監測，每日于9：00、19：00觀測兩次；流量采用間測方式進行監測，每月不少于1次，測流方式為流速斷面法（3條垂線），測速歷時大于100 s。灌區典型地塊退水口斷面采用流量過程線法推流。退水口流量采用巡測方式監測，有退水時隨時監測，退水口斷面流量較小時測流方式為直角三角型量水堰測流。根據文獻[12]，對于自由流直角三角堰，流量計算公式為：

Q=1.343H2.24（1）

式中，H=0.06～0.65 m。薄壁堰厚度1.5 mm，堰頂高0.5 m，頂寬0.5 mm。

2.1.3 地下水觀測井布設 灌區農田土壤性質、透水性能、地下水位埋深以及灌溉定額等因素對田間灌溉滲漏量產生綜合影響。在典型灌區水文地質條件下，因灌溉水滲漏致使地下水位變動，含水層中的重力水體積的變化在疊加降雨入滲因素后，可近似作為地下水補給量，亦是灌溉水滲漏回歸河道的水量。

本研究在典型地塊內挖掘5眼地下水觀測井，開展農田灌溉水下滲及對地下水動態影響試驗研究。地下水位觀測井位置如圖2所示，在湟水左岸邊設立直立式水尺1組共2支，1、2、3、4、5號井距河邊水尺P1距離分別為68.3、68.6、48.8、29.0、29.9 m。為使地下水位及河道水位在同一個高程系統內反映灌溉用水下滲及河道水位的變化情況，大峽灌區典型地塊設有2個水準點，經緯度分別為（E102°13′35.8″，N36°29′11.36″和E102°13′36.0″，N36°29′7.42″）埋深為1.5 m，高程分別為100.000 m和97.763 m。兩水準點相距約124 m。工作中每月對各水準點進行互校，同時校測河道水尺高程及地下水井口高程。

每次灌溉前一天觀測5眼地下水井水位，灌溉后期每日9：00、14：00、19：00觀測3次，地下水水位穩定后停止觀測。每次觀測地下水位時，同步觀測河道水位。地下水位觀測采用PD-26型便攜式激光測距儀結合懸垂式電子感應器人工觀測，激光測距儀測量精度為±2mm，測量范圍為0.2～60.0 m，激光等級為2級，波長為635 mm，工作溫度為-10～50 ℃。依照文獻[13]要求，每次監測地下水位應測量兩次，間隔時間不應少于1 min，當兩次測量數值之差小于0.02 m時，取兩次水位的平均值。當兩次測量偏差超過0.02 m時，應重復測量。在實際觀測中，2次測量偏差在0.005以內時，采用2次平均值，高于規范要求的標準。測量成果當場核查，及時繪制各個地下水井的水位過程線，發現反常及時補測，保證監測資料真實、準確、完整、可靠。灌溉期每日9：00、14：00、19：00觀測3次，水位穩定后每日9時觀測一次。

2.1.4 氣象因子監測 采用安裝在試驗現場的CR1000型自動監測氣象站測定空氣溫度、相對濕度、氣壓、風速、風向、太陽總輻射、降雨量和日照時間，測定時間間隔為每1或2 h測定一次，根據監測結果計算日平均或日累積值。各監測探頭技術指標參數如表2所示。試驗期間研究區氣象條件見表3。2.2 引排差法

2.4 數據處理

采用EXCEL進行數據整理和分析，采用SPSS統計軟件進行相關分析，采用CAD作圖。

3 結果與分析

3.1 基于引排差法耗水系數分析

3.1.1 地表引退水量 大峽灌區典型地塊進水口流量監測情況見表4，典型地塊退水口流量監測情況見表5。由表4可知，引水量中，斗渠1和斗渠2分別測流18次和16次，最大流量分別為0.103 m3/s和0.079 m3/s。由表5可知，退水口共計觀測44次，其中退水口1有水36次、退水口2有水1次；退水口3有水10次；退水口4有水31次；退水口5有水4次；退水口6有水20次。

大峽灌區典型地塊春灌期引退水量見圖3。大峽灌區典型地塊春灌期引水量共計4.916 2萬m3，其中斗渠1引水量3.922 6萬m3，斗渠2引水量0.993 6萬m3。退水量共計2.116 5萬m3。

3.1.2 地下退水量 試驗地下水水位變幅為灌溉前地下水位與灌溉后地下水位上升至最高水位之差，計算中剔除水位異常變化影響。春灌期1號井至5號井水位變幅分別為0.29、0.28、0.29、0.18和0.24 m，平均變幅Δh為0.26 m。典型地塊土壤質地為沙壤土，因此給水度μ參考中細沙下限值，取0.085；典型地塊面積F為19.33 hm2。經計算得到典型地塊春灌期灌溉后地下退水量為4 272.67m3。

3.1.3 引排差法計算耗水系數 綜合考慮地表和地下退水量后，根據公式8采用引排差法可計算得到大峽灌區典型地塊春灌期耗水系數為0.48。

3.2 基于最大蒸發量法耗水系數分析

3.2.1 參考作物蒸發蒸騰量變化規律 參考作物蒸發蒸騰量（ET0）在觀測時段內逐漸增大，達到最大后逐漸降低（圖4）。根據 ET0的變化過程，其變化趨勢可劃分為3個階段。第一階段為J=69～77 d，日值變化幅度較小，大部分ET0值在2～4 mm/d范圍內變化；第二階段為J=77～80 d，ET0較大，4～7 mm/d范圍內變化；第三階段為J=81～84 d，ET0一般為4～5 mm/d，變化幅度相對于第二階段小，比第一階段稍大。 ET0最大值為7.38 mm/d，出現在J=80（3月21日）；最小值為2.26 mm/d，出現在J=71（3月12日）。

3.2.2 參考作物蒸發蒸騰量影響因子 參考作物蒸發蒸騰量與太陽輻射、空氣平均溫度表現出極顯著的正相關性，與風速表現出顯著的正相關性，而與最大相對濕度表現出顯著負相關性（表6）。除了氣象因子，灌溉和降雨也是 ET0變化的重要影響因子，這與徐俊增等[14]研究結論相同。由圖4可知，灌溉后 ET0一般表現出增大的趨勢，原因在于灌溉一般發生在晴好天氣，輻射較大，氣溫較高，大氣蒸發力較大，因此ET0較大；試驗期間未發生降雨，因此未分析降雨對ET0的影響。

3.2.3 最大蒸發量法計算耗水系數 以春灌期（3月10-25日）典型地塊逐日參考作物蒸發蒸騰量累加值作為耗水量，根據公式30采用最大蒸發量法可計算得到耗水系數為0.49。

4 結論

大峽灌區典型地塊春灌期引水量4.916 2萬m3，退水量2.116 5萬m3。春灌期1號井至5號井地下水水位平均變幅為0.26 m，經計算得到典型地塊春灌期灌溉后地下退水量為272.67 m3。綜合考慮地表和地下退水量后，采用引排差法計算得到大峽灌區典型地塊春灌期耗水系數為0.48。

參考作物蒸發蒸騰量在觀測時段內逐漸增大，達到最大后逐漸降低。最大值為7.38 mm/d，出現在J=80（3月21日）；最小值為2.26 mm/d，出現在J=71（3月12日）。參考作物蒸發蒸騰量與太陽輻射、空氣平均溫度表現出極顯著的正相關性，與風速表現出顯著的正相關性，而與日最大相對濕度表現出顯著負相關性。除了氣象因子，灌溉也是 ET0變化的重要影響因子。以春灌期（3月10-25日）典型地塊逐日參考作物蒸發蒸騰量累加值作為耗水量，采用最大蒸發量法計算得到耗水系數為0.49。

傳統的引排差法是通過灌溉引退水計算耗水系數，為間接法；而最大蒸發量法計算耗水系數為直接法。引排差法對引排水量監測控制全面，基礎資料扎實可靠，可全面反映灌區引用水量情況，對加強灌溉管理、節約用水、灌區內部水量調配均有指導意義，但資料涉及面廣，工作量大，一年總結一次時效較慢。本研究中兩種方法計算得出的耗水系數非常接近，這說明在研究區范圍內耗水系數計算方面，最大蒸發量法可以對引排差法進行較好的驗證。

參考文獻：

[1] 曹進冬.以水為源，促發展惠民生保生態[J].青海科技，2013（2）：54-56.

[2] 王學全，盧 琦，高前兆.內蒙古河套灌區引用黃河水量分析[J].干旱區研究，2005，22（2）：146-151.

[3] 秦大庸，于福亮，裴源生.寧夏引黃灌區耗水量及水均衡模擬[J].資源科學，2003，25（6）：19-24.

[4] ALLEN R G， SMITH M， PERRIER A. An update for the definition of reference evapotranspiration[J]. ICID Bulletin，1994， 43（2）：1-34.

[5] TRAJKORIC S. Temperature-based approaches for estimating reference evapotranspiration[J]. Journal of Irrigation and Drainage Engineering，2005，131（4）：316-325.

[6] PEREIRA A R，PRUITT W O. Adaptation of the thornthwaite scheme for estimating daily reference evapotranspiration[J]. Agricultural Water Management， 2004，66： 251-257.

[7] 劉曉英，李玉中，王慶鎖.幾種基于溫度的參考作物蒸散量計算方法的評價[J].農業工程學報，2006，22（6）：12-18.

[8] SL13-2004，灌溉試驗規范[S].

[9] 高傳昌，朱明鋒，秦海霞.Hargraeves公式計算參考作物騰發量在新疆地區的適用性研究[J].灌溉排水學報，2008，27（3）：60-62.

[10] 胡順軍，潘 渝，康紹忠.Penman-Monteith與Penman修正式計算塔里木盆地參考作物潛在騰發量比較[J].農業工程學報，2005，21（6）：31-37.

[11] 王衛光，彭世彰，羅玉峰.參考作物騰發量的混沌性識別及預測[J].水利學報，2008，39（9）：1030-1036.

[12] SL 537-2011，水工建筑物與堰槽測流規范[S].

[13] SL 183-2005，地下水監測規范[S].

[14] 徐俊增，彭世彰，丁加麗，等.基于蒸滲儀實測數據的日參考作物蒸發騰發量計算方法評價[J].水利學報，2010，41（12）：1497-1505.