引黃畦灌土壤水分變化特征及作物耗水規律研究

孫耀民，王 勇，周宇偉，胡雅琪，趙永安

0 引 言

目前，全世界有90%的灌溉面積采用的是不同形式的地面灌溉方式。畦灌是通過田埂將耕地分割成長條的田塊，灌溉水流在重力作用下向畦尾推進并不斷滲入土壤的灌溉方式，田間工程簡單，投資小而被廣泛使用[1-3]。灌溉的目的是為了滿足作物正常生長需要而使土壤具有一定的含水率，一方面灌溉既要控制水分較好的分布在一定的土層范圍，另一方面也要使得水分最大限度地被作物吸收，減少灌溉水的深層滲漏。研究表明，畦田規格和灌水技術參數對灌水質量有重要的影響[4]。李益農等[5]以數值模擬方法求解田間水流運動方程，分析并討論了田面平整精度與畦灌系統性能間的定量關系。結果表明，田面平整精度改善下的灌溉效率可提高34%，灌水均勻度可提高28%。洪瑜等[6]研究了不同畦灌方式下冬小麥主要生育期的耗水規律，結果表明灌漿期小麥對水分需求最大，同時提出了適宜當地入畦單寬流量。本文通過開展冬小麥田間試驗，監測田間土壤含水率，對引黃畦灌條件下沿畦長方向的土壤剖面水分變化與冬小麥耗水規律進行了研究。

1 材料與方法

1.1 研究區概況

本田間試驗在山西省運城市尊村引黃灌區進行，尊村引黃工程設計灌溉面積11.066 7×104hm2，5個縣農業工業從中受益，是山西省重要的糧食、經濟作物種植基地。試驗區地勢平坦，土壤質地優良，光照充足，是典型的暖溫帶大陸性季風氣候，年內降水分布不均勻，多年平均降水量525 mm；日照時長為2 039.5 h，年均氣溫為13.3℃，無霜期212 d。運城市年平均水資源利用總量為106 795×104m3，其中地表水資源年均利用量26 177×104m3，提引黃河水16 956×104m3，提引黃河水約占地表水的65%；地下水年均利用量為80 618×104m3，占總用水量的75.5%，開發利用程度在全省最高。試驗區域糧食經濟作物種類繁多，其中夏玉米(6月初至10月初)和冬小麥(10月初至來年6月初)連作方式是最主要糧食作物種植模式。

1.2 測點布置及測定方法

以冬小麥作為研究對象，對田間土壤水分分布和作物耗水規律進行研究。使用TDR時域反射法對土壤水分進行測定，沿畦長方向進行測點布置，每隔50 m布置一個測點，對每個測點深100 cm的土層每10 cm進行一次含水率測定。試驗前對距測點50 cm取土，通過烘干法進行土壤含水率測量，測量結果用于TDR反射法的標定。土壤含水率監測時間為每周三與周日；若進行田間灌溉或有降雨時，在灌水降雨前后進行含水率監測；降雨量、氣溫、風速等氣象數據通過試驗區的小型氣象站獲取。

1.3 計算方法

有效灌溉或降水量的計算如下[4-5]：

P0=F-D

(1)

(2)

式中：P0為有效灌水量或降雨量，m；F為灌水量；D為土壤深層滲漏，mm；θ后為灌水后或降雨后的土壤含水率，%；θ前為灌水前或降雨前的土壤含水率，%；Si為所測土層厚度，mm。

日均作物騰發計算[6]：

(3)

式中：I為灌溉量，mm；P為降水量，mm；D為深層滲漏量，mm；θ1、θ2為在時間1和時間2時的體積含水率；Si為對應土層厚度，mm；%：t為兩個連續測量的時間間隔，d；n為層數。

參考作物騰發量使用Penman-Monteith公式計算[7]：

(4)

式中：ET0為參考作物騰發量，mm；Δ為飽和水汽壓與溫度關系曲線的斜率,kPa/℃；Rn為植被表面凈輻射量，MJ/m2d；G為土壤熱通量，MJ/m2d；γ為溫度計常數，kPa/℃；T為空氣平均溫度，℃；u2為在地面以上2 m高處的風速，m/s；es為空氣飽和水汽壓，kPa；ed為空氣實際水汽壓，kPa。

作物系數可以通過單作物系數方法和雙作物系數方法計算[8,12]。在本文中根據單作物系數方法，利用水平衡方法式(3)計算的ETc和Penman-Monteith公式計算的ET0計算作物系數，如下：

(5)

式中：ETc和ET0為在不同生長階段的作物騰發和參考作物騰發。

式(5)中作物系數Kc包括土壤和作物蒸發的影響。

2 結果與分析

2.1 冬小麥沿畦長土壤水分動態變化規律

圖1為所布置監測點處100 cm土層土壤體積含水率隨時間變化圖。隨著土層深度增加，表層土壤含水率具有明顯的高低變化，并且高值與低值的區域范圍較小，中層土壤含水率的變化范圍擴大，而底層土壤含水率在較大范圍內變化不明顯。0～10、10～20 cm處的土層在沿畦長方向上隨時間變化具有很好的一致性，同時也反映不同時間的灌水量變化情況；20～40 cm處土層體積含水率隨時間變化逐漸不明顯；而40 cm以下的土層基本不隨時間變化，也即不隨灌水量而變化，分析其原因耕作層(表土層)由于經常翻耕土質疏松，土壤滲透率好；30cm土層處受農機耕犁壓實形成犁底層，通氣透水性不良，影響耕作層與心土層的物質與能量交換。所以從此處開始，土壤體積含水率隨時間變化不明顯，40cm以下的土層土壤體積含水率與土層深度和灌水量的關系不明顯，從畦首到畦尾各觀測點各同層土壤的體積含水率趨勢相同，基本穩定 在0.25。

2.2 冬小麥生育期內騰發特征

圖2為沿畦長方向的冬小麥平均蒸騰量隨時間變化圖。冬小麥在返青期到收獲期，平均日騰發量在0.75～5.69 mm之間。

圖2 冬小麥各監測點平均騰發量隨時間分布圖

從圖2中可知，2月17日至4月8日，各監測點的日蒸騰量變化不大，基本在2.27 mm左右，此時段為冬季與初春時節，氣溫較低，風速較大，而冬小麥處于出苗-拔節期，用于生長的騰發量不大。隨著小麥進入快速生長階段，騰發量逐漸增大，加之此階段內氣溫逐漸升高，所以從4月8日開始小麥蒸騰量快速增加，期間5月10日左右，當地降雨量大導致蒸發量較小，小麥騰發量有所降低。5月20日以后，由于小麥進入成熟收獲期，用于作物本身的騰發量不斷減少導致此階段內小麥騰發量快速降低。另外，從圖2在所布置的6個觀測點的蒸騰量隨時間的變化趨勢具有一致性，通過計算同一時間不同監測點之間騰發量的變異系數，在全部28個時間點上，弱變異有8個，中等變異有20個，見圖3。同時比較冬小麥騰發量與變異系數隨時間變化的關系圖可知，騰發量與變異系數存在相反的變化趨勢，騰發量越大其相對應的變異系數就越小，反之騰發量小的其變異系數就越大。

圖3 冬小麥沿畦長各點變異系數隨時間分布圖

表1列出了冬小麥后4個生育階段騰發量ETc、參考作物騰發量ET0以及作物系數Kc。

表1 冬小麥返青-收獲期騰發量、參考作物騰發量以及作物系數

由表1可以看出，拔節-抽穗期是冬小麥騰發量最大的生育期，平均為108.42 mm；灌漿-收獲期冬小麥騰發量最小，平均為70.71 mm。通過氣象數據并根據Penman-Monteith公式計算的參考作物騰發量，得到后4個生育期的作物系數，冬小麥在拔節-抽穗期和抽穗-灌漿期的生長系數大于1，為1.17和1.14，說明在這兩個生育期內小麥的騰發量大于參考作物騰發量，返青-拔節期的騰發量基本與參考作物騰發量相等；4個生育期總的騰發量為333.47 mm，平均作物系數為1.05，實際騰發量略大于參考作物騰發量。

冬小麥從返青期到收獲期，灌溉水量以及降雨總共為369.21 mm。第一次灌溉田間深層滲漏為21.21 mm，占灌水量的22.09%；第二次田間灌溉深層滲漏為26.02 mm，占灌水量的27.1%，說明在灌水及降雨量大的情況下極易發生深層滲漏。

3 結 論

1) 冬小麥日平均騰發量隨生育期的不同而變化，在全生育期初期變化較小，隨著小麥的快速生長逐漸增加，至成熟收獲期騰發量逐漸減小，平均總騰發量為350.17 mm。

2) 在所布置的6個觀測點處的騰發量在時間軸上的變化有很好的一致性。同時比較冬小麥騰發量與變異系數隨時間的變化，二者存在相反的變化趨勢，騰發量越大其相對應的變異系數就越小；反之騰發量小的其變異系數就越大。

3) 本文采用田間試驗的方法，對引黃灌區冬小麥糧田沿畦長方向土壤剖面水分變化與作物騰發等內容進行了研究。但本試驗所研究的只是各個生育期，并沒有對全育期進行研究分析，所以所得結論還有待進一步的試驗驗證。