噴灌與畦灌對拉薩河谷農田小氣候和蒸散量的影響分析

張紫森，湯鵬程，徐 冰，李仙岳，楊 波，王國帥，李澤坤

（1.內蒙古農業大學水利與土木建筑工程學院，呼和浩特 010018；2.中國水利水電科學研究院內蒙古陰山北麓草原生態水文國家野外科學觀測研究站，北京 100038；3.中國水利水電科學研究院牧區水利科學研究所，呼和浩特 010020）

0 引 言

拉薩河谷位于青藏高原“一江兩河”流域，具有強輻射、高寒缺氧、干旱少雨的特點，其水利設施相對落后，田間灌溉多以大水漫灌和畦灌為主[1]。有研究表明拉薩河谷平均年潛在蒸散量在1 190 mm 以上，且正在以8.21 mm/10a 的速度增加[2]，蒸散量的過高不利于作物生長甚至可能導致干旱減產。噴灌通過改變SPAC 系統之間的水熱交換過程，影響農田小氣候并減少蒸散量，使作物生長環境得到改善的同時達到節水灌溉的目的。

適宜的農田小氣候有利于作物高產、穩產，反之則會抑制作物生長[3]。蒸散量是表面能量平衡過程的重要組成部分，能客觀反映作物生理生長狀況，對于決定灌水需求，制定灌溉制度，開展農業水資源管理具有重要指導意義[4]。近年來國內外學者對噴灌對農田小氣候和蒸散的影響進行了大量的研究。劉海軍指出，噴灌的氣溫、溫度梯度、VPD 和蒸騰強度均小于地面灌溉，而作物光合速率、產量和水分利用效率均大于地面灌溉，這種改變有利于作物生長[5-7]。王慶改[8]、王迪[9]、Jose Cavero[10]在噴灌影響田間蒸散的機理研究中指出，噴灌通過影響作物冠層內外水汽交換過程進而抑制田間蒸散，同時證明了這種冠層截留水量在作物生長過程中的有效性。程光遠[11]、閆澤宇[12]結合噴灌條件下大豆不同生育期的需耗水規律，采用模型模擬制定灌溉制度。由此可見，內地從噴灌影響蒸散機理到采用模型模擬優化灌溉制度均較為完善，但由于拉薩河谷高效節水灌溉建設起步較晚[13]，噴灌相關基礎性研究較少，迫切需要加強。

本文以拉薩市林周縣為典型研究區，對比分析兩種灌溉方式蒸散量與農田小氣候的關聯性，明晰噴灌影響蒸散量的原因，并估算噴灌對農田小氣候和蒸散量的改變量，為發展拉薩河谷噴灌技術和制定噴灌條件下灌溉制度提供理論依據。

1 材料與方法

1.1 試驗區概況

試驗區位于西藏自治區拉薩市林周縣縣畜牧場內，地理坐標緯度東經91°11′，北緯29°54′。林周縣地處拉薩河谷農區，平均海拔3 650 m，屬溫帶高原季風氣候。年平均氣溫7.4 ℃，晝夜溫差較大，年日照時數在3 500 h以上，太陽輻射強，平均年潛在蒸散量在1 190 mm 以上。試驗區土壤類型為沙壤土，容重為1.4 g/cm3，計劃濕潤層最深30 cm。

1.2 試驗設計

拉薩河谷燕麥的種植時間一般在5-9月，本試驗開展時間為2021年7-8月，此時燕麥正處于拔節期～抽穗期，燕麥品種為青海444，播量為225 kg/hm2，播種方式為撒播，底肥為復合肥，底肥量180 kg/hm2。本試驗沿道路南北兩側設置噴灌與畦灌兩種灌水處理，噴灌處理采用PY-5022 搖臂式噴頭，直徑為7 mm，噴灑半徑為6～12 m，工作水量1.96～2.49 m3/h。灌溉過程中采用兩組設備對農田小氣候和蒸散量進行監測，每組設備由1 臺氣象站和1 臺蒸滲儀構成，每組設備中雨量筒和微型蒸滲儀的位置相距緊密，具體布置如圖1所示。試驗結束后分別對兩種灌溉方式的農田小氣候與蒸散量進行灰色關聯度分析，并以畦灌各項實測數據為自變量對噴灌各項實測數據進行回歸分析。

圖1 田間試驗布置圖（單位：m）Fig.1 Field trial layout

噴灌與畦灌同時進行，本試驗為與當地農牧民的灌水時間保持一致，灌水時間大部分在日間11∶00-19∶00選取且在午間13∶30-14∶30 停止灌溉，具體灌水時間見表1。噴灌試驗通過水表監測水量，為了消除噴灌過程中蒸發漂移損失引起的誤差，試驗中以雨量筒中水量為蒸滲儀中實際進入水量；在畦灌條件下，通過水表監測的水量與灌溉面積確定單位面積灌水量，根據微型蒸滲儀面積折算出蒸滲儀的灌水量，用量筒將水一次性灌入蒸滲儀中。兩組處理的其他農業措施均保持一致。

表1 噴灌與畦灌試驗灌水時間表Tab.1 Irrigation schedule for sprinkler and border irrigation tests

1.3 測定指標

試驗觀測內容主要包括氣象指標、蒸散量和噴灌灌溉量。

（1）氣象指標。采用兩臺NHQXZ601 全天候氣象站，進行農田小氣候觀測，觀測內容包括：大氣溫度、相對濕度、氣壓、2 m 處風速、總輻射、反輻射、凈輻射。每隔10 min 記錄一次，試驗結束后對30 min 內的大氣溫度、相對濕度、氣壓、2 m處風速求平均值，對總輻射、反輻射、凈輻射求和。

（2）蒸散量。采用兩臺LYS80 微型蒸滲儀每2 min 監測一次桶內土柱質量，得到土壤貯水量的變化量。蒸滲儀尺寸（長×寬×深）為0.8 m×0.8 m×0.8 m。

（3）噴灌灌溉量。為了消除噴灌過程中蒸發漂移損失帶來的水量平衡法中灌溉量的誤差，實際灌溉量采用氣象站中雨量筒收集到的水量為微型蒸滲儀中實際進入水量，雨量筒尺寸（長×寬×高）為0.2 m×0.2 m×0.4 m，噴灌水量如圖2所示。

圖2 微型蒸滲儀實際進入水量Fig.2 The actual amount of water entered by the micro lysimeter

1.4 計算方法

1.4.1 基于微型蒸滲儀的蒸散量計算方法

本文中實際蒸散量的計算采用水量平衡法進行，具體公式如下：

式中：ET為作物實際蒸散量；P為階段降雨量；I為階段灌溉量；△SWS為階段土壤貯水量變化；Q為階段地下水補給量和滲漏量。由于試驗需要，上述指標均折算為mm。

1.4.2 灰色關聯度的計算方法與分析方法

本試驗以蒸散量為參考序列，氣象指標為比較序列，首先求出各序列的差序列與各差序列最大值、最小值，然后計算關聯系數與關聯度，最后將關聯度從大到小排序。關聯度越大說明比較序列與參考序列的貢獻度越大，關聯系數與關聯度的計算公式如下，詳細計算方法不做贅述[14-17]。

1.4.3 數據處理

采用Excel 2007 對數據進行整理數據和灰色關聯度分析，采用Spss 23對數據進行回歸分析并建立回歸模型。

2 結果與分析

2.1 噴灌與畦灌對農田小氣候和蒸散量的影響分析

2.1.1 噴灌與畦灌對大氣溫度和相對濕度的影響

圖3是2021年7-8月噴灌與畦灌平均溫度和相對濕度的對比圖。結合表1，溫度在11∶00～19∶00 這段時間相較于內地偏低，這是由于高海拔因素所導致的。溫度整體變化相對穩定，變化區間集中在13～24°C之間，其中噴灌變化區間為13～20°C，畦灌變化區間為13～24°C，噴灌對比畦灌可以使溫度的波動區間縮短，對溫度峰值的削弱較為明顯。造成這種現象的主要原因為，水滴蒸發吸熱帶走了由于強輻射產生的部分熱量，使溫度降低。夜間灌溉條件下，噴灌與畦灌對溫度均明顯降低且無顯著變化。濕度方面，噴灌與畦灌的整體趨勢相對一致，但由于強輻射引起的水滴蒸發造成空氣中的水分含量提升，所以畦灌條件下變化區間為40%～78%；噴灌條件下為59%～85%，噴灌條件下平均相對濕度比畦灌約高13%。

圖3 噴灌與畦灌對大氣溫度和相對濕度變化關系Fig.3 Relationship between sprinkler irrigation and border irrigation on atmospheric temperature and relative humidity

2.1.2 噴灌與畦灌對總輻射、反輻射和凈輻射的影響

噴灌與畦灌條件下總輻射、凈輻射與反輻射的變化情況如圖4所示。結合表1在11∶00～19∶00 這段時間內，噴灌與畦灌對總輻射的改變量并不明顯，在圖中的表現為不同符號的重疊。凈輻射方面，兩種灌溉方式對凈輻射的改變并不顯著但畦灌的凈輻射略高于噴灌的凈輻射。反輻射方面，畦灌的反輻射要高于噴灌的反輻射值，分析此現象產生的原因是：水滴均勻地落在土壤表面緩慢提高了土壤含水率，改變了下墊面原有的地表覆蓋類型，導致地表輻射長波的能力和地表反射輻射通量均減少[18]。圖中第4、6、9、13、23 次試驗輻射量均接近零且沒有明顯變化，可能是由于氣象因子中云量的增加而引起地面太陽輻射的減少[19]，進而導致反射率的變化，影響反輻射和凈輻射。第19 次試驗由于在夜間進行所以輻射通量均為零且總輻射、凈輻射、反輻射均為負值。

圖4 噴灌與畦灌對總輻射、反輻射和凈輻射變化關系Fig.4 The relationship between sprinkler irrigation and border irrigation on total radiation,anti-radiation and net radiation

2.1.3 噴灌與畦灌對蒸散量的影響

噴灌與畦灌燕麥蒸散量的變化關系如圖5所示。噴灌條件下蒸散量的變化區間集中在0.01～0.65 mm 之間；畦灌條件下蒸散量的變化區間集中在0.05～0.95 mm 之間。第20 和21 次試驗結果偏大，分析原因由于此時風速較大為4 m/s，使燕麥莖稈和葉片擺動造成了試驗誤差。

圖5 噴灌與畦灌對蒸散量變化關系Fig.5 Relationship between sprinkler irrigation and border irrigation on evapotranspiration

2.2 噴灌與畦灌農田小氣候與蒸散量的灰色關聯度分析

本試驗以蒸散量為參考序列，氣象指標作為比較序列，求出關聯度并進行排序，關聯度如表2所示。畦灌條件下蒸散量與氣象指標的關聯度大小排序為：大氣溫度＞氣壓＞風速＞相對濕度＞總輻射＞反輻射＞凈輻射；噴灌條件下蒸散量與氣象指標的關聯度大小排序為：氣壓＞大氣溫度＞反輻射＞相對濕度＞總輻射＞凈輻射＞風速。排序越靠前，對蒸散量的影響越大，反之對蒸散量的影響越小。畦灌條件下排序靠前的是大氣溫度、氣壓、風速、相對濕度；噴灌條件下排序靠前的是氣壓、大氣溫度、反輻射、相對濕度。綜合分析噴灌與畦灌氣象指標的關聯度排序，大氣溫度、氣壓和相對濕度為影響噴灌與畦灌蒸散量的主要指標；風速、總輻射、反輻射、凈輻射為影響噴灌與畦灌蒸散量的次要指標。

表2 噴灌與畦灌蒸散量與氣象因素的關聯度Tab.2 Correlation between evapotranspiration and meteorological factors of sprinkler irrigation and border irrigation

主要指標方面大氣溫度、相對濕度代表的熱量因素直接影響作物蒸散過程，而氣壓代表的空氣中有效氣體濃度（O2、CO2）因素直接影響作物氣孔導度，氣孔導度的大小與光合作用和呼吸作用有著密切聯系[20]，間接影響作物的蒸散過程，拉薩河谷的氣壓約為內地的50%～60%，所以相對于內地氣體濃度因素影響蒸散量不能忽視。次要指標方面在風速和反輻射上產生了較為明顯的差異，畦灌處理中風速與蒸散量的關聯度為0.729，反輻射與蒸散量的關聯度為0.682；噴灌處理中風速與蒸散量的關聯度為0.638，反輻射與蒸散量的關聯度為0.714，但噴灌和畦灌均不能直接影響風速，這就在一定程度上提高了反輻射對蒸散量的貢獻度，降低了風速對蒸散量的貢獻度。

2.3 噴灌與畦灌對農田小氣候和蒸散量的回歸分析

以畦灌條件下的氣象指標和蒸散量為自變量對噴灌條件下的氣象指標和蒸散量進行回歸分析，同時結合各氣象指標的最大值和最小值進行量化計算，回歸分析結果見表3。由表3可以看出，噴灌的各項指標與畦灌均為正相關，但噴灌對氣壓、風速、總輻射基本不產生影響，對凈輻射減少為7%，只能略微減少凈輻射量，對反輻射的減少顯著為37%，為0.01～0.32 MJ/（m2·h）；對蒸散量的減少顯著為29%，具體改變量為0.03～0.38 mm/h；溫度和濕度的減少顯著，但由于涉及常數項所以不便于量化影響。為便于描述噴灌對比畦灌在溫度和濕度的改變量，對表3中溫度和濕度模型設置固定截距為0 分別優化，優化結果見表4。由此可見，噴灌對溫度降低17%，在2～4°C變化；對濕度提升28%，在11%～21%變化。

表3 噴灌與畦灌農田小氣候和蒸散量的回歸模型Tab.3 Regression model of microclimate and evapotranspiration in sprinkler irrigation and border rrigation

表4 噴灌與畦灌優化溫度和濕度模型Tab.4 Optimal temperature and humidity models for sprinkler and f border irrigation

3 結 論

通過對噴灌與畦灌實測氣象指標與蒸散量進行灰色關聯度分析，大氣溫度、氣壓、相對濕度的排序均在前列，所以大氣溫度、氣壓、相對濕度為影響噴灌與畦灌蒸散量的主要指標；風速、總輻射、反輻射、凈輻射為次要指標。主要指標揭示出大氣溫度、相對濕度代表的熱量因素和氣壓代表的氣體濃度（O2、CO2）因素，對拉薩河谷的蒸散量有很大影響。

通過實測噴灌與畦灌數據進行回歸分析，主要指標方面噴灌使溫度降低2～4°C，相對濕度提升11%～21%，對氣壓不產生影響；次要指標方面噴灌使反輻射約減少0.01～0.32 MJ/（m2·h），對其他氣象指標影響較小，綜合主要指標和次要指標分析，日間噴灌對比畦灌使蒸散量減少0.03～0.38 mm/h。主要指標方面大氣溫度越高則空氣中水分含量越低，所以蒸散量升高，反之，蒸散量降低，而相對濕度越高則空氣中的水分越趨于飽和導致蒸散過程被抑制，所以使蒸散量降低；次要指標方面反輻射的減少會直接影響下墊面環境，導致主要指標的溫度降低、濕度升高，進而降低蒸散量。對于高海拔、高蒸散的拉薩河谷，噴灌可以減少蒸散量，有利于作物生長的同時降低灌溉水量，進而指導制定適宜的噴灌灌溉制度。

本文在結合農田小氣候、蒸散量與作物指標共同分析方面有所欠缺。結合作物分析有利于從作物生理角度揭示農田小氣候和蒸散量的改變對作物生長機理的影響，需要后期的研究中加以補充和重視。