趙曉棟
(山西萬水水利勘測設計咨詢有限公司,山西 太原 030024)
農田灌溉受到降水時空分布不均的影響較大,當前國內外學者在農田灌溉空間尺度效應方面的研究成果頗豐,大部分學者主要從作物蒸散發角度以及利用旬降雨尺度進行灌區灌溉制度合理性評價,但對于設計灌溉保證率受時間尺度效應的影響研究較少。設計灌溉保證率是灌區灌溉制度確定的基礎,也是體現灌區灌溉水利用效率的主要指標,其取值大小及合理性直接影響灌區灌溉設施規模及灌溉面積。根據《灌溉與排水工程設計規范》,灌區可根據具體灌水方式、作物類型查表選取設計灌溉保證率,以便灌區灌溉工程設計工作的展開。但是該規范并未將時間尺度對設計灌溉保證率的影響考慮進去,所給出的設計灌溉保證率取值的準確性和適用性值得商榷。
以山西省隰縣2018年果水配套工程為例,工程建設項目由8個灌溉片區組成,分別為桃坡果水配套工程、下崖底果水配套工程、員家莊果水配套工程、瓦窯坡果水配套工程、水堤果水配套工程、下司徒果水配套工程、竹干果水配套工程、馬家果水配套工程。具體實施地點為桃坡、下崖底、員家莊、瓦窯坡、竹干、下司徒、水堤、馬家等8個自然村,因地制宜建設節水灌溉工程。本研究所需氣象數據主要取自中國氣象數據網1960年~2020年的逐日數據。
根據研究目的,將隰縣2018年果水配套工程區氣象數據分為降水數據以及根據修正彭曼公式所計算出的參考作物需水量數據兩個部分[1],其中工程區參考作物需水量主要根據以下公式確定:

式中:ET0為作物需水量,mm/d;Rn為土壤靜輻射量,MJ/(m2·d);G為土壤熱通量,MJ/m2·d;T為日平均氣溫,℃;2為風速,m/s;es為飽和水汽壓,kPa;ea為實際水汽壓,kPa;為飽和水汽壓與溫度關系曲線的斜率值,kPa/℃;γ為干濕溫度計常數,kPa/℃。
即項目區用水量在多年間能得到保證的概率,通常按照正常供水年份在總年份中的占比進行量化,可采用的經驗頻率計算公式如下:

式中:P為設計灌溉保證率,%;m為設計灌溉用水量所對應的供水年份數,a;n為總年份數,a。
根據該區需水過程,并綜合考慮灌區內水資源、降水等情況,本研究中設計灌溉標準主要以灌溉保證率法所確定的值為準,設計灌溉保證率取75%,設計基準年取2017 年。
項目區灌溉制度的有關參數主要根據當地實際經驗確定。隰縣2018年果水配套工程區代表性作物為梨樹和蘋果樹,兩種作物的種植比例分別為35%和65%。在設計灌溉保證率下,根據隰縣2018年果水配套項目區作物種植狀況及結構調整規劃,應主要采用低壓管道輸水的管灌方式。
項目區梨樹和蘋果樹兩種典型作物的灌溉應主要根據《農田低壓管道輸水灌溉工程技術規范》(GB /T 20203-2006),同時結合項目區實際情況,灌水定額設計參數具體見表1。

表1 作物灌水定額設計參數表
按照所確定出的灌水定額,并結合項目區現狀灌溉規律,同時結合相關試驗結果,進行項目區管灌方式下具體灌溉制度的確定,見表2。

表2 作物灌溉制度表
通過P-Ⅲ曲線進行山西省隰縣2018年果水配套工程降雨量數據擬合,設計灌溉保證率分別取P=25%、P=50%、P=75%及P=90%的典型年進行時間尺度降雨量計算。
2.1.1 典型年降雨量分析
對1960 年~2020 年項目區年降雨量進行頻率排列計算,并將實際降雨頻率與設計灌溉保證率之差不超出±5%范圍的情況確定為典型年,具體確定結果見表3,通過對表中結果的分析可以看出,當設計灌溉保證率相同時,不同典型年所對應的降雨量取值之間的差異并不大。

表3 設計灌溉保證率下典型年降雨量
當設計灌溉保證率P=25%時,比較年降雨量較為接近的情況所對應的月降水量值。比較結果表明,在設計灌溉保證率相同時,各典型年實際降水量之間雖然不存在較大差異,但其所對應的月降水量卻相差較大,例如,4月份降雨量最大值和最小值分別為1987 年的227.8 mm和1965 年的42.4 mm,差距達到5 倍之多;9 月份降雨量最大值和最小值分別為2008 年的273.1 mm和2004 年的55.4 mm,差距最大達到5 倍之多。這也表明以年尺度為標準選擇設計灌溉保證率時,必須充分考慮降雨量時空分布不均的問題,否則,工程區典型作物若不能充分吸收和利用區域內降雨資源,便會造成部分雨水資源浪費。
設計灌溉保證率P=25%的情況下比較各典型年4月、7月、8月上、中、下旬的實際降雨量不難看出,4月上旬降雨量最大最小值分別為1974年的18.7 mm和2008 年的0.8 mm;4月中旬降雨量最大最小值分別為1989 年的85.2 mm和2002 年的0.6 mm;4 月下旬降雨量最大最小值分別為1983 年的8.0 mm和1960 年的0 mm。頻率相同情況下各典型年旬降雨量值差距較大,而且極值分布年份也不盡相同。在P=25%的設計灌溉保證率下,7月、8月上、中、下旬降雨量極值變動趨勢規律與4 月基本一致,故省去分析過程。
隰縣2018 年果水配套工程區典型作物生長周期均在1年以內,故在進行典型作物關鍵生長期確定時必須精確至旬或月,為此,該項目區應按照年尺度進行設計灌溉保證率的選定和分析,據此所確定出的保證程度不高。此外,《灌溉與排水工程設計規范》中在選取典型年的過程中必須根據年內分配最不利原則,根據該原則所確定出的典型年中作物生育期內降水量稀少,考慮到作物生長過程中需水性方面的要求,必須增大管道灌溉系統的可供水量[2]。但是,這種情況下,會因此產生一些不利的影響,例如不利于項目區典型作物對降雨資源的充分利用;增大灌溉工程投資,與灌區實際生產需要嚴重脫節等。
2.1.2 月降雨量分析
對1960年~2020年項目區4月、7月、8月降雨量進行頻率排列計算,同時參照頻率不超出±5%的要求進行所對應典型年的選擇。對于設計灌溉保證率相同的情況,主要依據時間尺度選擇典型年,并比較不同時間尺度所對應的4 月、7 月、8月降水量極差。當工程區設計灌溉保證率取P=25%、P=50%、P=75%及P=90%時,4月、7月、8月 降水量極差年尺度大于月尺度,充分表明,根據年尺度進行典型年選取時,月降水量極差比所對應的月尺度下降水量極差大,降雨量分配的不均勻程度十分嚴重。
國內一般根據項目區社會經濟狀況、主要作物類型及灌溉方式等進行設計灌溉保證率的確定,國外農業生產中設計灌溉保證率通常在75%及以上,部分國家甚至達到98%。我國《灌溉與排水工程設計規范》規定的灌區設計灌溉保證率應保持在50%~95%之間,并要求設計灌溉保證率應隨現狀灌溉設施、主要作物類型、水源狀況等進行調整。對于且缺水地區旱作等而言,設計灌溉保證率可以適當降低;對于風水區水稻等作物,設計灌溉保證率應適當提高;灌溉設施較為先進的地區,設計灌溉保證率也應高一些。遺憾的是,現行相關規范中,設計灌溉保證率的確定主要依據降雨量排列頻率,所以,依據降雨量或來水量確定設計灌溉保證率的做法會造成灌區供水時間和需水時間尺度的不同步。
在進行山西省隰縣2018 年果水配套工程區試驗田灌溉制度制定時,典型年的選取仍采用的是相關規范所規定的典型年法,將所確定出的典型年內實際降雨量按照時間尺度由大至小劃分為季尺度、月尺度、旬尺度和日尺度等,并根據不同降雨量尺度確定所對應的灌溉制度。通過水量平衡法確定項目區不同典型年代表性作物灌溉制度,并基于此,結合代表性作物生育期內濕潤層計劃深度具體計算儲水量可能的上下限。當設計灌溉保證率相同情況下,典型年各年份代表性作物需水量差異并不大,只是在灌溉定額與降雨資源利用率方面存在較大差異;且隨時間尺度的減小,項目區代表性作物灌溉定額表現出明顯的下降趨勢,對降雨資源的利用程度呈增大趨勢。這說明項目區內代表性作物灌溉定額和降雨資源利用量之間存在互補關系[3],即降雨資源利用程度越高,灌溉定額越大;反之亦然。
山西省隰縣2018年果水配套工程區不同時間尺度下灌溉定額和降雨量利用情況見表1,其中設計灌溉保證率P=75%時,典型年僅為1989年,故未進行顯著性分析。通過方差分析進行時間尺度影響灌溉定額及降雨資源利用程度的顯著性研究,檢驗水平α=0.05。根據分析可以看出,當設計灌溉保證率P=25%時,時間尺度為季時年灌溉定額顯著較高,且其年降雨量利用率比其余時間尺度低。這表明該項目區灌溉制度制定時,豐水年若采用季時間尺度,則會增大降雨資源的浪費,故不適合以季節作為降雨計算的時間尺度。當設計灌溉保證率P=90%時,隨著時間尺度從季變化到月、旬、日,降雨資源利用程度表現出增大趨勢,年灌溉定額則以較小幅度下降。這表明當降雨計算的時間尺度縮短時,項目區作物對降雨資源的利用程度增大,灌溉水量會得到很大程度的節省。

表4 不同時間尺度下灌溉定額和降雨量利用情況
本研究利用山西省隰縣2018 年果水配套工程區1960 年~2020 年降雨量資料,并按照年降雨量與作物需水所對應的月、旬降雨量進行了設計灌溉保證率時間尺度效應分析。結果顯示,該工程區所處地理位置主要在降水季節性變化大且年內降雨資源分布不均衡的季風氣候區,同時按照月、旬等較小時間尺度進行典型年的選擇,這種處理方式能在一定程度上解決降雨資源時空分布不均的問題,顯著提升項目區內典型作物生育期對降雨資源的利用程度。而且,隨著時間尺度的減小,針對該工程區典型作物灌溉制度,年灌溉定額表現為降低趨勢,但年降雨資源利用率反而上升。隰縣2018年果水配套工程區選擇灌溉設計保證率時主要依據《灌溉與排水工程設計規范》,該規范中雖然對灌區設計灌溉保證率的計算過程以及選擇原則作出了明確規定,但對時間尺度效應缺乏充分考慮。由本研究分析可知,時間尺度的選擇對灌區灌溉水保證程度及節水潛力和效果均存在較大程度的影響,為提升設計灌溉保證率保證程度,并充分利用各類水資源,必須進行合理時間尺度的選擇與論證。