區(qū)域灌溉需水量時空變化特征及其影響因素研究——以魯西北地區(qū)為例

鄭潤橋，費良軍，介飛龍，劉利華，劉 念，閻思宇

（西安理工大學(xué) 西北旱區(qū)生態(tài)水利國家重點實驗室，西安 710048）

區(qū)域灌溉需水量時空變化特征及其影響因素研究——以魯西北地區(qū)為例

鄭潤橋，費良軍*，介飛龍，劉利華，劉 念，閻思宇

（西安理工大學(xué) 西北旱區(qū)生態(tài)水利國家重點實驗室，西安 710048）

【目的】探尋魯西北地區(qū)灌溉需水量時空分布特征及影響因素?！痉椒ā窟x取冬小麥為研究對象，利用Penman-Monteith公式對魯西北地區(qū)灌溉需水量進(jìn)行計算，應(yīng)用Mann-Kendall突變檢驗法分析其變化趨勢，得出魯西北地區(qū)灌溉需水量的時間分布規(guī)律，使用反距離權(quán)重插值法得出灌溉需水量空間分布規(guī)律，并利用Pearson相關(guān)分析法對各氣象因素對灌溉需水量的影響程度進(jìn)行分析比較?！窘Y(jié)果】①魯西北地區(qū)的年平均灌溉需水量為339 mm，年際變化呈上升趨勢；②灌溉需水量空間變化呈帶環(huán)狀分布，各生育期內(nèi)規(guī)律不相同；③相對濕度和日照時數(shù)是影響魯西北地區(qū)灌溉需水量較為重要的因素，分別呈極強負(fù)相關(guān)與較強正相關(guān)?！窘Y(jié)論】魯西北地區(qū)冬小麥灌溉需水量逐年緩慢上升，由東北到西南遞增；各站點降水量在播種期、返青期和拔節(jié)期內(nèi)均不滿足需水要求，需要人工灌溉補給。

灌溉需水量；時空變化；影響因素

0 引 言

【研究意義】山東作為我國農(nóng)業(yè)大省，人均水資源占有量僅為334 m3，低于全國平均水平，水資源短缺情況嚴(yán)重[1]。除農(nóng)業(yè)用水量大和效率低等原因外，因山東省屬暖溫帶季風(fēng)氣候，降水集中在夏季，水資源時空分布不均勻也是造成水資源短缺的重要原因[2]。因此能否對農(nóng)作物進(jìn)行高效灌溉和合理配置農(nóng)業(yè)用水規(guī)模，對灌溉需水量的時空變化規(guī)律及影響因素進(jìn)行分析是關(guān)鍵所在。

【研究進(jìn)展】目前，對灌溉需水量的研究主要包括3個方面：①單考慮蒸散量的時空分布規(guī)律。如環(huán)海軍等[3]在計算冬小麥各生育階段需水量時采用了Penman-Monteith模型和單作物系數(shù)法，同時計算了有效降水量、水分盈虧指數(shù)，并對冬小麥在不同生育期的各種灌溉指標(biāo)包括需水量、有效降水量等指數(shù)的時空變化等進(jìn)行了計算。王景雷等[4]建立了基于DEM的山東省冬小麥作物需水量空間分布的數(shù)學(xué)模型，同時也涉及ArcGIS的空間分析功能，對其蒸散量的空間關(guān)系進(jìn)行了分析。②從灌溉需水量的計算確定及其時空關(guān)系特征入手，尋其特征規(guī)律。如劉鈺等[5]選取國內(nèi)近300個氣象站點30 a的氣象數(shù)據(jù)，對我國近30余種主要經(jīng)濟作物灌溉需水量進(jìn)行計算，得出我國不同地區(qū)主要經(jīng)濟作物的灌溉需求指數(shù)。丁浩[6]對陜西寶雞峽灌區(qū)的氣候條件因素、冬小麥和夏玉米等主要經(jīng)濟作物情況和寶雞峽灌區(qū)的灌溉農(nóng)業(yè)現(xiàn)狀進(jìn)行整體研究，得出影響寶雞峽灌區(qū)農(nóng)業(yè)需水的主要因素和時空變化規(guī)律。③氣象因素對灌溉需水量的影響研究，如宋妮等[7]以河南地區(qū)18個氣象站的氣象資料為基礎(chǔ)，對河南1961—2012年的主要經(jīng)濟作物需水量與不同氣象影響因素進(jìn)行了分析比較，得出結(jié)論為：影響河南地區(qū)主要經(jīng)濟作物的氣象因素為日照時間與最高氣溫這2個因素。張山清等[8]在計算參考作物蒸散量時選取的氣象數(shù)據(jù)為1959—2007年的新疆吐魯番地區(qū)的基本氣象數(shù)據(jù)，并采用M-K突變檢驗和線性回歸分析相結(jié)合的方式研究了不同氣象條件對作物需水量的影響，研究表明影響作物需水量的主要因素為風(fēng)速、有效降雨量和平均氣溫?！厩腥朦c】結(jié)合文獻(xiàn)資料可知，對于魯西北地區(qū)灌溉需水量的相關(guān)研究較少且多為單方面的計算分析，制約了當(dāng)?shù)氐默F(xiàn)代節(jié)水農(nóng)業(yè)發(fā)展。【擬解決的關(guān)鍵問題】為此，通過灌溉需水量的計算和影響因素分析，從年際變化尺度、生育期內(nèi)月際尺度下開展山東省西北地區(qū)灌溉需水量時間、空間變化規(guī)律和影響因素的綜合研究。

1 材料與方法

1.1 研究區(qū)概況

山東省內(nèi)位于黃河以北的區(qū)域稱之為魯西北地區(qū)，地區(qū)面積大約3萬km2。魯西北地區(qū)北連河北，西臨河南，地貌單一，幾乎全部為坦蕩的黃泛平原，屬華北大平原的一部分[9]。降水主要受到大氣環(huán)流、季風(fēng)氣候和地形環(huán)境的影響，年際變化劇烈，年內(nèi)分配極不均勻，1965—2015年平均降水量為526 mm，且集中在夏季，6—8月的降水量占全年降水量的70%，多春旱夏澇。山東省年平均蒸發(fā)量為879 mm，蒸發(fā)量較大[10]。魯西北地區(qū)主要作物以冬小麥、夏玉米、大豆、棉花為主，其中冬小麥種植分布較廣，種植面積最大，占作物種植面積的54%[11]。

1.2 灌溉需水量計算

1.2.1 參考作物蒸發(fā)蒸騰量（ET0）

采用Penman-Monteith公式[12]計算標(biāo)準(zhǔn)參考作物的潛在蒸散量。影響因素包括最高溫度、最低溫度、平均溫度、日照時間、風(fēng)速、降水量等。計算式為：

式中：ET0為參考作物蒸散量（mm/d）；Rn為冠層表面凈輻射（MJ/（m2·d））；G為土壤熱通量（MJ/（m2·d））；Δ為飽和水汽壓與溫度曲線斜率（kPa/℃）；γ為濕度計常數(shù)（kPa/℃）；T為平均氣溫（℃）；U2為2 m高處的風(fēng)速（m/s）；ex為空氣飽和水汽壓（kPa）；ea為空氣實際水汽壓（kPa）。

1.2.2 作物需水量（ETc）

一般采取田間測定進(jìn)行作物需水量的確定，同時也可進(jìn)行理論計算得出[13]。本文采用作物系數(shù)法計算作物需水量，具體計算式為：

式中：Kc為作物系數(shù)；ETc為作物需水量（mm/d）；ET0為參考作物蒸散量（mm/d）。

1.2.3 有效降水量（Pe）

使用簡化法計算有效降水量[14-15]，其計算式為：

式中：P為日降水量（mm/d）；Pe為有效降水量（mm/d）；ETc為作物需水量（mm/d）。

1.2.4 灌溉需水量（IS）

研究區(qū)內(nèi)大部分地區(qū)的地下水埋深較大，平均達(dá)8.3 m，部分站點達(dá)20 m以上[16-17]，忽略地下水補給因素，灌溉需水量可由作物需水量減去有效降水量算出[18]，其計算見式（4）；而濟南地區(qū)為“泉城”，其區(qū)域內(nèi)地下水埋深較淺，地下水補給條件充沛。計算其灌溉量時，將泉域補給灌溉區(qū)的灌溉需水量簡化為0，全濟南地區(qū)的灌溉需水量以式（4）計算得出后，乘以權(quán)重系數(shù)（非泉域灌溉區(qū)面積占總灌溉面積的比值），即得最終的灌溉需水量。

式中：ETc為作物需水量（mm/d）；N為生育期日數(shù)（d）；Pe為有效降水量（mm/d）；IS為灌溉需水量（mm）。

1.3 M-K突變檢驗法

Mann-Kendall檢驗可以用來判斷序列中是否存在氣候的突變，并且可以確定其突變的時間，即突變年份[19]。假定Xn為一時間序列數(shù)據(jù)（X1，X2，…，Xn），其不同時間點的數(shù)據(jù)是相互獨立、隨機的樣本數(shù)據(jù)；假設(shè)Hi是雙邊檢驗，對于任意的k，j≤n，且k≠j，Xi和Xj的分布是不相同的，檢驗統(tǒng)計量S計算式[20]：

其中，

S服從標(biāo)準(zhǔn)正態(tài)分布，方差Vαr(S)=n(n-1)(2n+5)/18。當(dāng)ngt;10時，通過式（7）計算統(tǒng)計量Z值：

對于統(tǒng)計量Z，當(dāng)Zgt;0時表示當(dāng)前時間序列呈增加趨勢；Zlt;0時表示當(dāng)前時間序列呈減少趨勢。統(tǒng)計量|Z|≥1.28、|Z|≥1.64、|Z|≥2.32時，表示當(dāng)前時間序列分別通過了置信度為90%、95%、99%的顯著性檢驗。

Mann-Kendall檢驗還可以用于時間序列的突變節(jié)點檢測，通過構(gòu)造秩序列：

其中：

定義統(tǒng)計變量：

式中：E(Sk)=k(k-1)/4；Var(Sk)=k(k-1)(2k+5)/72；UFk為標(biāo)準(zhǔn)正態(tài)分布，對于給定顯著性水平α，若｜UFk｜＞Uα/2則表明序列存在明顯的趨勢變化，將時間序列x按逆序排列，再按照式（10）計算，同時使：

對于統(tǒng)計序列UFk和UBk進(jìn)行進(jìn)一步分析時間序列的變化趨勢以及其突變的時間節(jié)點。若UFkgt;0，則表明序列呈增加趨勢；UFklt;0則表示呈減少趨勢；當(dāng)其達(dá)到設(shè)定閾值時，表示當(dāng)前時間序列具有顯著性的變化趨勢。如果UFk和UBk的序列統(tǒng)計曲線存在交點，并且交點位于設(shè)定閾值的范圍之內(nèi)時，那么此時交點出現(xiàn)的時間即對應(yīng)時間序列突變時間點[21]。

1.4 空間插值法

基于相似性原理，反距離加權(quán)插值法(IDW)以插值點同樣本點之間的距離為權(quán)重進(jìn)行加權(quán)平均。采樣點離插值點越近，分配給它的權(quán)重就越大[22]。采用反距離加權(quán)法對魯西北灌區(qū)冬小麥需水量進(jìn)行插值，得到魯西北灌區(qū)冬小麥需水量的空間分布圖。

1.5 數(shù)據(jù)來源

本研究選取魯西北地區(qū)5個代表氣象站點1961—2015年的逐日氣象資料，所使用的氣象數(shù)據(jù)均由中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供。收集的氣象資料包括：2 m處風(fēng)速、日照時間、日降雨量、相對濕度、最低氣溫、平均氣溫和最高氣溫等。

冬小麥生育期的劃分依據(jù)國家氣象局編制的《農(nóng)業(yè)氣象觀測規(guī)范》[23]進(jìn)行，冬小麥作物系數(shù)基于《北方地區(qū)主要作物灌溉用水定額》[24]確定，共劃定5個生育期，播種期（1013—1215）內(nèi)作物系數(shù)Kc為0.60；越冬期（1216—0226）內(nèi)作物系數(shù)Kc為0.40；返青期（0227—0411）內(nèi)作物系數(shù)Kc為1.17；成熟期（0610）作物系數(shù)Kc為0.40。

2 結(jié)果與分析

2.1 灌溉需水量的計算

利用式（1）—式（4），得出魯西北地區(qū)惠民、朝城、濟南3個氣象站點1986—2015年冬小麥年際灌溉需水量。由于德州、羊角溝站部分年份數(shù)據(jù)存在連續(xù)缺測，因此選取德州站點1986—2000年、羊角溝站點1986—2008年的氣象數(shù)據(jù)進(jìn)行計算，計算結(jié)果如圖1、表1所示。在選取的年份時段內(nèi)，德州地區(qū)冬小麥年際灌溉需水量在313～355 mm之間，平均值為336 mm，灌溉需水量最大的年份為1997年，最小的年份為1991年；惠民地區(qū)冬小麥年際灌溉需水量在315～378 mm之間，平均值為339 mm，灌溉需?水量最大的年份為2011年，最小的年份為1997年；羊角溝地區(qū)冬小麥年際灌溉需水量在324～390 mm之間，平均值為345 mm，灌溉需水量最大的年份為2007年，最小的年份為1999年；朝城地區(qū)冬小麥年際灌溉需水量在277～319 mm之間，平均值為303 mm，灌溉需水量最大的年份為1993年，最小的年份為2011年；濟南地區(qū)冬小麥年際灌溉需水量在356～387 mm之間，平均值為371 mm，灌溉需水量最大的年份為2001年，最小的年份為1997年；魯西北地區(qū)多年平均灌溉需水量最小地區(qū)為朝城，多年平均灌溉需水量最大地區(qū)為濟南。結(jié)合魯西北地區(qū)各氣象站的地理分布可以得出，中部多年平均灌溉需水量最大，東北部多年平均灌溉需水量適中，西南部多年平均灌溉需水量最小。

圖1 灌溉需水量年際變化Fig.1 Annual variation of irrigation water demand

表1 冬小麥多年平均灌溉需水量Table 1 Perennial average irrigation water requirement of winter wheat mm

2.2 灌溉需水量的時間變化規(guī)律研究

2.2.1 灌溉需水量的年際變化規(guī)律

使用Mann-Kendall突變檢驗法和線性傾向估計法對各氣象站點的年際灌溉需水量進(jìn)行突變檢驗和趨勢變化分析，最終得出魯西北灌溉區(qū)各氣象站點的冬小麥灌溉需水量的年際變化趨勢情況。5個站點灌溉需水量M-K檢驗如圖2所示，各站點年際灌溉需水量線性傾向估計值如表2所示。由圖2（a）可知，德州地區(qū)冬小麥年灌溉需水量變化趨勢在1994、1996年發(fā)生2次突變，結(jié)合線性傾向估計可知，德州地區(qū)1986—1994年為上升趨勢，1994—1996年內(nèi)為下降趨勢，1996年后又轉(zhuǎn)為上升趨勢，根據(jù)顯著性檢驗知全時間序列為不顯著變化；由圖2（b）可知，惠民地區(qū)2011年是唯一突變年，結(jié)合線性傾向估計可知，1986—2011年為下降趨勢，其中1996—1998年趨勢變化顯著，2011年發(fā)生突變后2011—2014年轉(zhuǎn)為不顯著上升趨勢；由圖2（c）可知，羊角溝地區(qū)全序列顯著性檢驗均未能通過0.05顯著性檢驗，因此全序列變化趨勢不顯著，而存在2005年為突變年份，結(jié)合線性傾向估計可知，1986—2005年為不顯著下降趨勢，2005—2008年趨勢變?yōu)轱@著性上升趨勢；由圖2（d）可知，朝城站點灌溉需水量的M-K檢驗圖中UFk與UBk曲線交點較多，不易進(jìn)行判別，因此對數(shù)據(jù)進(jìn)行滑動檢驗后發(fā)現(xiàn)：2005年為突變年份，1986—2005年為上升趨勢，其中1989—1998年內(nèi)變化趨勢顯著，2005年發(fā)生突變后，2005—2014年趨勢變?yōu)椴伙@著下降趨勢；由圖2（e）得出，濟南地區(qū)1986—2014年內(nèi)灌溉需水量時間變化趨勢均未發(fā)生突變，全序列為不顯著上升趨勢。

表2 各站點年際灌溉需水量線性傾向估計Table 2 Annual linear trend estimation of irrigation water demand for each site

注*表示通過α=0.05顯著性檢驗。

圖2 各站點年際灌溉需水量M-K檢驗Fig.2 Mann-Kendall test of inter-annual irrigation water requirement of each station

2.2.2 灌溉需水量的月際變化規(guī)律

將魯西北地區(qū)5個氣象點的冬小麥生育期內(nèi)多年平均月灌溉需水量進(jìn)行計算，詳情見圖3所示，各氣象站點生育期內(nèi)冬小麥多年平均灌溉需水量各月百分比如表3所示。總體上看，魯西北地區(qū)冬小麥在生育期內(nèi)灌溉需水量月際分布并不均勻，表現(xiàn)為從10月開始到次年6月結(jié)束，各站點灌溉水量的變化趨勢基本類似：先小幅度增加后減少，再大幅度增加和大幅度減少。具體各個站點情況見表3，德州地區(qū)冬小麥灌溉需水量最大的時期為拔節(jié)期，灌溉需水量模比系數(shù)為0.40，灌溉需水量最小的時期為越冬期，灌溉需水量模比系數(shù)為0.10；惠民地區(qū)冬小麥生育期內(nèi)灌溉需水量最大的時期為拔節(jié)期，模比系數(shù)為0.44，最小時期為越冬期，模比系數(shù)為0.08；羊角溝地區(qū)冬小麥生育期內(nèi)灌溉需水量最大的時期為拔節(jié)期，模比系數(shù)為0.41，最小時期為成熟期，模比系數(shù)為0.07；朝城地區(qū)冬小麥生育期內(nèi)灌溉需水量最大的時期為拔節(jié)期，模比系數(shù)為0.41，最小時期為越冬期，模比系數(shù)為0.08；濟南地區(qū)冬小麥生育期內(nèi)灌溉需水量最大的時期為拔節(jié)期，模比系數(shù)為0.39，最小時期為成熟期，模比系數(shù)為0.08。

對比5個氣象站點的灌溉需水量月際變化規(guī)律，并分析其原因：前半年的小幅度變化是因為此時處于冬小麥的越冬時期，冬小麥停止生長，需水量和灌溉量都很小；3月，魯西北地區(qū)氣溫回升，冬小麥度過越冬期進(jìn)入返青期，開始正常生長，因此作物需水量增大，灌溉需水量也開始大幅度上升；4、5月時，冬小麥的生長進(jìn)入拔節(jié)期，此時大量需要水分，是冬小麥作物需水量最大的時期，并且伴隨氣溫持續(xù)升高，蒸騰作用增大，因此該時期的灌溉用水量達(dá)到峰值，6月10日前后，魯西北地區(qū)冬小麥已經(jīng)進(jìn)入成熟期，不再進(jìn)行灌溉，因此灌溉需水量開始減小。

圖3 魯西北地區(qū)冬小麥灌溉需水量月際變化Fig.3 Monthly variation of irrigation water requirement of winter wheat in northwest Shandong province

表3 各站點各生育期月灌溉需水量及模比系數(shù)Table 3 Irrigation water requirement and modulus ratio coefficient of each station in each growth period

2.3 灌溉需水量的空間變化規(guī)律研究

由于站點的經(jīng)緯度、降雨、地形地貌、氣象、水文等條件不同，因此冬小麥的灌溉需水量在空間分布上也存在一定差異。使用ArcGIS運用反距離權(quán)重插值法對魯西北地區(qū)不同氣象站點的冬小麥多年平均灌溉需水量進(jìn)行空間插值，最終繪制出魯西北地區(qū)冬小麥灌溉需水量空間分布圖。如圖4所示，魯西北地區(qū)冬小麥灌溉需水量成帶環(huán)狀分布，朝城地區(qū)多年平均灌溉需水量最小，在303～318 mm范圍內(nèi)上下波動，處于第一級；其次為德州、惠民、羊角溝地區(qū)，多年平均灌溉需水量在318～334 mm范圍內(nèi)上下波動，處于第二級；濟南地區(qū)的多年平均灌溉需水量最大，其數(shù)值波動范圍為356～371 mm，處于第五級。總體灌溉需水量呈“南多北少，西多東少”，“由西南至東北逐漸遞增”的分布規(guī)律。

圖4 魯西北地區(qū)冬小麥多年平均灌溉需水量空間分布Fig.4 Spatial distribution of perennial average irrigation water requirement of winter wheat in northwest Shandong

以4個主要生育期為研究尺度，對冬小麥不同生育期的灌溉需水量進(jìn)行反距離權(quán)重插值，并繪制各生育期內(nèi)灌溉水量空間分布圖5。由圖5得知，各生育期內(nèi)冬小麥灌溉需水量空間分布有所差異。播種期內(nèi)，朝城地區(qū)冬小麥的灌溉需水量最小，在35～40 mm上下波動，處于第一等級，羊角溝地區(qū)冬小麥播種期灌溉需水量為40～44 mm，處于第二等級，惠民地區(qū)灌溉需水量則在43～45 mm范圍內(nèi)上下波動，位于第三等級，德州地區(qū)灌溉需水量則在45～48 mm范圍內(nèi)上下波動，位于第四等級，濟南地區(qū)灌溉需水量最大，范圍為48～51 mm，位于第五等級；越冬期內(nèi)，惠民、朝城地區(qū)的灌溉需水量較小，為26～30 mm，處于第一等級，羊角溝地區(qū)灌溉需水量位于第三等級，為32～34 mm，德州地區(qū)灌溉需水量位于第四等級，為34～37 mm，濟南地區(qū)灌溉需水量最大，位于第五等級，范圍為37～40 mm；返青期內(nèi)，位于第一等級的地區(qū)為德州和朝城地區(qū)，其灌溉需水量較小，為88～94 mm，惠民地區(qū)次之，灌溉需水量為94～96 mm，位于第二等級，羊角溝地區(qū)灌溉需水量則在98～102 mm范圍內(nèi)上下波動，位于第四等級，濟南地區(qū)最大，為102～106 mm，位于第五等級；拔節(jié)期內(nèi)，朝城地區(qū)灌溉需水量位于第一等級，為124～128 mm，德州地區(qū)灌溉需水量為133～137 mm，位于第三等級，羊角溝、惠民、濟南地區(qū)灌溉需水量為141～148 mm，處于第五等級；成熟期內(nèi)，朝城、羊角溝地區(qū)灌溉需水量較小，為24～27 mm，位于第一等級，濟南地區(qū)為28～29 mm，位于第四等級，德州和惠民地區(qū)灌溉需水量最大，范圍為29～31 mm，處于第五等級。

對比5個氣象站在不同生育期內(nèi)灌溉需水量空間分布情況發(fā)現(xiàn)：返青期與拔節(jié)期的魯西北冬小麥灌?溉需水量空間分布規(guī)律變化較明顯，結(jié)合站點的多年氣象數(shù)據(jù)得知，返青期與拔節(jié)期正處春季，此時的灌溉需水量最大，受不同年份冬季溫、濕度的影響，不同站點在春季的起始雨熱條件也不同，若遇旱冬，返青期田間起始濕度較低，加之降水不足，在春季的冬小麥灌溉量則會增加。因此推測造成不同生育期內(nèi)灌溉需水量空間分布規(guī)律不同的主要因素為各站點春冬季節(jié)的溫濕雨熱條件不同。

圖5 各生育期冬小麥灌溉需水量空間分布Fig.5 Spatial distribution of irrigation water requirement of winter wheat at different growth stages

2.4 灌溉需水量的影響因素研究

對魯西北地區(qū)冬小麥生育期內(nèi)逐年灌溉需水量與對應(yīng)生育期內(nèi)站點的年平均風(fēng)速、年降水量、年日照時間、年相對濕度、年平均氣溫、年最高氣溫、年最低氣溫之間進(jìn)行Pearson相關(guān)性分析。冬小麥灌溉需水量與各個氣象因素相關(guān)系數(shù)和顯著性見表4。

通過計算可知，各站點灌溉需水量與各氣象因素的相關(guān)性存在差異。其大小關(guān)系為：德州：相對濕度＞日照時間＞風(fēng)速＞降雨量＞最高溫度＞最低溫度＞平均溫度；惠民：相對濕度＞日照時間＞平均溫度＞風(fēng)速＞降雨量＞最低溫度＞最高溫度；羊角溝：相對濕度＞日照時間＞風(fēng)速＞降雨量＞最高溫度＞最低溫度＞平均溫度；朝城：相對濕度＞日照時間＞風(fēng)速＞最低溫度＞降雨量＞平均溫度＞最高溫度；濟南：相對濕度＞日照時間＞降雨量＞風(fēng)速＞最低溫度＞最高溫度＞平均溫度。綜合來看，影響灌溉需水量的主要氣象因素為日照時間和相對濕度，其中相對濕度和冬小麥灌溉需水量呈極強負(fù)相關(guān)關(guān)系，即相對濕度越大，灌溉需水量越小，而日照時間與冬小麥灌溉需水量呈較強正相關(guān)，即日照時間越長，冬小麥的灌溉需水量越大。

表4 冬小麥灌溉需水量與各氣象因素的相關(guān)系數(shù)和顯著性Table 4 Correlation coefficient and significance between irrigation water requirement of winter wheat and meteorological factors

注*表示通過α=0.05顯著性檢驗；**表示通過α=0.01顯著性檢驗。

3 討 論

本文得出的魯西北地區(qū)平均冬小麥作物需水量（517 mm）基本符合劉鈺等[25]華北區(qū)冬小麥生育期內(nèi)作物需水量的合理區(qū)間400～500 mm，本文計算出的魯西北地區(qū)平均灌溉需水量為339 mm，符合其給出的華北地區(qū)凈灌溉需水量范圍150～350 mm；再與任鴻瑞等[26]在魯西北的冬小麥耗水量試驗結(jié)果比較，其得到魯西北地區(qū)冬小麥10 a內(nèi)平均耗水量為450 mm，降水占比35%，推得灌溉量為292 mm。綜合以上文獻(xiàn)的相似與差異分析得出：①有效降水量的大小與所選取的計算時段長度有關(guān)，簡化計算時選取旬為計算時段長度。為保障結(jié)果精度，本文采用逐日數(shù)據(jù)計算，因此所得灌溉需水量較大。②由于時間局限性，任鴻瑞等[26]在魯西北研究區(qū)做耗水試驗時，選取了1994—2002年時段長度，以禹城實驗站為測量點測得魯西北冬小麥的耗水量，本文選取研究區(qū)內(nèi)5個站點近30 a的氣象數(shù)據(jù)進(jìn)行推算灌溉蓄水量，因此結(jié)果在精度與代表性上有所不同。

本文對魯西北地區(qū)冬小麥的灌溉需水量進(jìn)行時空規(guī)律研究中，指出冬小麥灌溉需水量呈現(xiàn)不顯著上升趨勢的結(jié)論與劉曉英等[27]在華北地區(qū)的研究結(jié)果一致。不同的是，在時間尺度上，本文利用M-K突變檢驗確定其突變點，再以突變點為分段點，對一致性不同的時間序列進(jìn)行了線性傾向估計，得出各站點冬小麥趨勢不同的年際變化規(guī)律及趨勢顯著程度，同時也得出生育期月尺度下灌溉需水量存在變化：灌溉所需水資源無法長期滿足灌溉要求時，最先應(yīng)確保冬小麥的拔節(jié)期的水量。在10—11月，山東冬小麥開始進(jìn)入播種期，灌溉需水量維持在40 mm即可。12—2月，冬小麥進(jìn)入越冬期，自身對生長所需水量要求很低，基本上不需要進(jìn)行灌溉。3—5月，冬小麥灌溉需水量變化趨勢為快速上升，所以灌溉量也應(yīng)及時跟隨加大，增加到135 mm。到6月夏季中旬，魯西北地區(qū)的冬小麥進(jìn)入收獲階段，基本上不需要灌溉；在空間尺度下，本文得出冬小麥的灌溉需水量呈“西南多、東北少，由西南至東北逐漸遞增”的分布規(guī)律，從朝城、德州，到惠民、羊角溝，再到濟南，灌溉水量逐漸增多，再將研究尺度調(diào)整至各生育期內(nèi)灌溉量的空間分布規(guī)律，本文發(fā)現(xiàn)在不同生育期內(nèi)，魯西北地區(qū)各站點灌溉需水量的空間變化趨勢也存在差異。

本文認(rèn)為在日照時間、風(fēng)速、相對濕度、最高氣溫、最低氣溫、平均氣溫、降雨量這7個氣象因素中，影響魯西北地區(qū)冬小麥灌溉需水量的最顯著因素為日照時間和相對濕度，相對濕度越大，冬小麥的灌溉需水量越??；日照時間越長，冬小麥的灌溉需水量越大。結(jié)論與馬婉棣[28]的結(jié)果一致。另外，魯西北地區(qū)各個灌溉分區(qū)內(nèi)各個氣象因素的影響率排名也存在一定差異。具體原因為不同灌溉分區(qū)的經(jīng)緯度及海拔高程的不同造成了雨熱條件不同，從而不同程度地影響了灌溉需水量。

4 結(jié) 論

1）魯西北地區(qū)冬小麥灌溉需水量逐年上升，西南區(qū)域高于東北區(qū)域。

2）播種、返青期和拔節(jié)期內(nèi)降水無法滿足魯西北地區(qū)冬小麥的需水要求，不同生育期內(nèi)灌溉需水量空間分布存在差異，農(nóng)業(yè)水資源配置時可作以參考。

3）溫濕條件是魯西北冬小麥灌溉需水量的主要影響因素，日照時間與灌溉需水量呈較強正相關(guān)；相對濕度與灌溉需水量呈極強負(fù)相關(guān)。

[1] 李明輝.山東糧食生產(chǎn)水資源配置及優(yōu)化策略研究[D].泰安: 山東農(nóng)業(yè)大學(xué),2019.

LI Minghui.Study on Water Resources Allocation and Optimization Strategies for Grain Production in Shandong[D].Tai’an: Shandong Agricultural University,2019.

[2] 孟林,劉新建,鄔定榮,等.華北平原夏玉米主要生育期對氣候變化的響應(yīng)[J].中國農(nóng)業(yè)氣象,2015,36(4): 375-382.

MENG Lin,LIU Xinjian,WU Dingrong,et al.Responses of Summer Maize Main Phenology to Climate Change in the North China Plain[J].Chinese Journa of Agrometeorology,2015,36(4): 375-382.

[3] 環(huán)海軍,楊再強,劉巖,等.魯中地區(qū)冬小麥水分盈虧及灌溉需水量的時空變化特征[J].干旱氣象,2016,34(5): 866-872.

HUAN Haijun,YANG Zaiqiang,LIU Yan,et al.Spatial and temporal variation of crop water surplus or deficit and irrigation water requirement for winter wheat in middle Shandong Province[J].Journal of Arid Meteorology,2016,34(5): 866-872.

[4] 王景雷,孫景生,周祖昊,等.基于DEM的作物需水量估算方法[J].灌溉排水學(xué)報,2005,24(5): 36-38.

WANG Jinglei,SUN Jingsheng,ZHOU Zuhao,et al.A method for evaluation crop water requirements based on digital elevation model[J].Journal of Irrigation and Drainage,2005,24(5): 36-38.

[5] 劉鈺,L S Pereira.對FAO推薦的作物系數(shù)計算方法的驗證[J].農(nóng)業(yè)工程學(xué)報,2000(5): 26-30.

LIU Yu,L S Pereira.Validation of FAO methods for estimating crop coefficients[J].Transactions of the CSAE,2000(5): 26-30.

[6] 丁浩.氣候變化條件下寶雞峽灌區(qū)主要作物需水量時空演變及節(jié)水潛力分析[D].楊凌: 西北農(nóng)林科技大學(xué),2017.

DING Hao.Effect of climate change on crop water requirement and water saving potential of Baojixia irrigation district in Northwestern China[D].Yangling: Northwest Aamp;F University,2017.

[7] 宋妮,孫景生,王景雷,等.河南省冬小麥需水量的時空變化及影響因素[J].應(yīng)用生態(tài)學(xué)報,2014,25(6): 1 693-1 700.

SONG Ni,SUN Jingsheng,WANG Jinglei,et al.Temporal and spatial variation of water requirement of winter wheat and its influencing factors in Henan Province[J].Chinese Journal of Applied Ecology,2014,25(6): 1 693-1 700.

[8] 張山清,普宗朝,宋良孌,等.吐魯番地區(qū)氣候變化對參考作物蒸散量的影響[J].中國農(nóng)業(yè)氣象,2009,30(4): 532-537,542.

ZHANG Shanqing,PU Zongchao,SONG Liangluan,et al.Effects of climate change on evapotranspiration of reference crops in Turpan Region[J].Chinese Agrometeorology,2009,30(4): 532-537,542.

[9] 楊學(xué)斌,代玉田.魯西北地區(qū)近50年氣溫和降水量的變化趨勢分析[J].中國農(nóng)業(yè)氣象,2009,30(S1): 25-28.

YANG Xuebin,DAI Yutian.Analysis on Temperature and Precipitation in Northwest of Shandong in Last 50 Years[J].Chinese Journal of Agricultural Meteorology,2009,30(S1): 25-28.

[10]陳一.山東省降水蒸發(fā)時空變化特征及水資源演變影響因素研究[D].濟南: 山東大學(xué),2021.

CHEN Yi.Spatial-temporal variation characteristics of precipitation and evaporation and influential fctors of water resources evolution in Shandong Province[D].Ji’nan: Shandong University,2021.

[11]趙慶慶,姜魯光,李文葉,等.魯西北平原冬小麥種植格局時空變化:2000—2014[J].國土資源遙感,2017,29(2): 173-180.

ZHAO Qingqing,JIANG Luguang,LI Wenye,et al.Spatial-temporal pattern change of winter wheat area in northwest Shandong Provinceduring 2000—2014[J].Remote Sensing for Land and Resources,2017,29(2): 173-180.

[12]杜江濤,張楠,龔珂寧,等.基于單作物系數(shù)法指導(dǎo)南疆棉花膜下滴灌試驗研究[J].中國農(nóng)村水利水電,2021(11): 141-148.

DU Jiangtao,ZHANG Nan,GONG Kening,et al.Experimental research on the irrigation system of mulched-cotton by drip irrigation in south Xinjiang based on the single crop coefficient[J].China Rural Water and Hydropower,2021(11): 141-148.

[13]ALLEN R G,SMITH M,PERRIER A,et al.An update for the definition of reference evapotranspiration[J].ICID Bulletin,1994,43(2):1-34.

[14]SMITH M.Computer program for irrigation planning and management[J].Food amp; Agriculture Organization of the United Nations,1992(7): 20-21.

[15]劉鈺,蔡甲冰,蔡林根,等.黃河下游灌區(qū)農(nóng)田灌溉制度與供需平衡分析[J].水利學(xué)報,2005(6): 701-708.

LIU Yu,CAI Jiabing,CAI Lingen,et al.Analyses of irrigation scheduling and water balance for irrigation districtin lower reaches of Yellow River[J].Journal of Water Conservancy,2005(6): 701-708.

[16]王金曉.山東省魯西北地區(qū)地下水資源評價[D].北京: 中國地質(zhì)大學(xué),2014.

WANG Jinxiao.Evaluation of resources of groundwater in northwest in Shandong Province[D].Beijing: China University of Geosciences,2014.

[17]陳特.魯西北地區(qū)引黃灌溉工程地下水回灌綜合評價研究[D].濟南:濟南大學(xué),2017.

CHEN Te.A thesis submitted to the university of Ji’nan in partial fulfillment of the requirements for the degree of master of engineering[D].Ji’nan: Ji’nan University,2017.

[18]BROUWER C,HEIBLOEM M.Irrigation water management:Irrigation water needs[J].Training Manual,1986(19): 63-68.

[19]冷曼曼,張志強,于洋,等.昕水河流域徑流變化及其對氣候和人類活動的響應(yīng)[J].水土保持學(xué)報,2020,34(3): 113-119,128.

LENG Manman,ZHANG Zhiqiang,YU Yang,et al.Runoff variation of Xinshui River watershed and its response to climate change and human activities[J].Journal of Soil and Water Conservation,2020,34(3):113-119,128.

[20]張林齊,任立良,江善虎,等.1982—2015年渭河流域植被變化特征及氣候因素影響[J].水文,2018,38(2): 66-72.

ZHANG Linqi,REN Liliang,JIANG Shanhu,et al.Vegetation change characteristics in Weihe river basin during 1982—2015 and its related climatic factors[J].Journal of Hydrology,2018,38(2): 66-72.

[21]龐家泰,段金亮,張瑞,等.2000—2019年渭河流域植被覆蓋度時空演變特征及氣候響應(yīng)[J].水土保持研究,2021,28(5): 230-237.

PANG Jiatai,DUAN Jinliang,ZHANG Rui,et al.Characteristics of spatiotemproal evolution and climate response of vegetation cover in the Wei River Basin from 2000 to 2009[J].Research of Soil and Water Conservation,2021,28(5): 230-237.

[22]馬貴仁,屈忠義,王麗萍,等.基于ArcGIS空間插值的河套灌區(qū)土壤水鹽運移規(guī)律與地下水動態(tài)研究[J].水土保持學(xué)報,2021,35(4):208-216.

MA Guiren,QU Zhongyi,WANG Liping,et al.Research on soil water and salt movement and groundwater dynamics in Hetao Irrigation distirct based on ArcGIS spatial interpolation[J].Journal of Soil and Water Conservation,2021,35(4): 208-216.

[23]國家氣象局.農(nóng)業(yè)氣象觀測規(guī)范[M].北京: 氣象出版社,1993.

National Weather Service.Code for agricultural meteorological observation[M].Beijing: Meteorological press,1993.

[24]段愛旺.北方地區(qū)主要農(nóng)作物灌溉用水定額[M].北京: 中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2004.

DUAN Aiwang.Quota of irrigation water for major crops in northern China[M].Beijing: China Agricultural Science and Technology Press,2004.

[25]劉鈺,汪林,倪廣恒,等.中國主要作物灌溉需水量空間分布特征[J].農(nóng)業(yè)工程學(xué)報,2009,25(12): 6-12.

LIU Yu,WANG Lin,NI Guangheng,et al.Spatial distribution characteristics of irrigation water requirement for main crops in China[J].Transactions of the CSAE,2009,25(12): 6-12.

[26]任鴻瑞,羅毅.魯西北平原冬小麥和夏玉米耗水量的實驗研究[J].灌溉排水學(xué)報,2004,23(4): 37-39.

REN Hongrui,LUO Yi.The experimental research on the water consumpion of winter wheat and summer maize in the northwest plain of Shandong Province[J].Journal of Irrigation and Drainage,2004,23(4): 37-39.

[27]劉曉英,林而達(dá).氣候變化對華北地區(qū)主要作物需水量的影響[J].水利學(xué)報,2004(2): 77-82,87.

LIU Xiaoying,LIN Erda.Effects of climate change on water requirement of main crops in North China[J].Journal of Water Conservancy,2004(2): 77-82,87.

[28]馬婉棣.山東省主要農(nóng)作物灌溉需水量時空分布規(guī)律研究[D].泰安:山東農(nóng)業(yè)大學(xué),2019.

MA Wandi.Study on temporal and spatial distribution of irrigation water demand for main crops in Shandong Province[D].Tai’an:Shandong Agricultural University,2019.

Spatiotemporal Variation in Crop Demand for Irrigation Water and Its Determinants: Take Northwest Shandong Province as an Example

ZHENG Runqiao,FEI Liangjun*,JIE Feilong,LIU Lihua,LIU Nian,YAN Siyu
(Institute of Water Resources and Hydro-electric Engineering,Xi’an University of Technology,Xi’an 710048,China)

【Objective】Agricultural production in most regions in northern China relies on irrigation.Understanding spatiotemporal variation of the demand of different crops for water is critical to improving agricultural management.Taking the northwest Shandong province as an example,this paper studies the change in its demand for irrigation.【Method】We took winter wheat as the model plant,and calculated its water demand based on the Penman-Monteith formula.Spatiotemporal variation in the water demand was analyzed using the Mann-Kendall mutation test,the inverse distance weighting and the interpolation method.The factors that affect the water demand was calculated using the Pearson correlation analysis method.【Result】①Annual average irrigation water demand of the winter wheat is 339 mm,but has been increasing.②Spatially,distribution of irrigation water demand is ring-shaped across the studied region,with the shape of the rings varying with growth season.③The relative humidity and sunshine duration are the factors that affect the irrigation water most,with the former negatively correlated while the latter positively correlated with the irrigation water demand,both at significant level.【Conclusion】The irrigation water demand of winter wheat in northwest Shandong province has been increasing.Spatially,it increases from the northeast to the southwest.Precipitation cannot meet the demand of the wheat during its seeding,returning green and elongation stages.Topping up the soil water by irrigation is hence required.

irrigation water requirement; spatiotemporal variation; influencing factors

S274.4

A

10.13522/j.cnki.ggps.2021453

1672-3317（2022）04-0084-09

OSID：

鄭潤橋,費良軍,介飛龍,等.區(qū)域灌溉需水量時空變化特征及其影響因素研究：以魯西北地區(qū)為例[J].灌溉排水學(xué)報,2022,41(4): 84-92.

ZHENG Runqiao,FEI Liangjun,JIE Feilong,et al.Spatiotemporal Variation in Crop Demand for Irrigation Water and Its Determinants: Take Northwest Shandong Province as an Example[J].Journal of Irrigation and Drainage,2022,41(4): 84-92.

2021-09-22

國家自然科學(xué)基金項目（52079105）；陜西省水利科技項目（2020SLKJ-19）

鄭潤橋（1997-），男。碩士研究生，主要從事農(nóng)業(yè)水資源研究。E-mail: zhengrunqiao@163.com

費良軍（1963-），男。教授，博士生導(dǎo)師，主要從事節(jié)水灌溉理論和水資源利用研究。E-mail: feiliangjun2008@163.com

責(zé)任編輯：白芳芳