基于CROPWAT模型的灌溉需水分析

張志堅 錢俊

摘要：為研究西非基礎資料匱乏地區的灌溉需水預測，提出了基于CROPWAT模型的灌溉規劃設計思路與方法，從作物需水量、降雨模型、土壤參數、種植規劃以及灌溉流量設計等方面展開了分析，提出了項目灌溉設計流量。結果表明：①對于資料匱乏的國際工程項目，利用CROPWAT開展分析工作具有較充分的理論基礎和較高的計算精度，具有一定的實用性;②作物種植計劃的變化對于灌溉毛需水量峰值（GIWRpeak）的影響較大，而灌溉系數、工作時間等參數對于設計流量的影響較大，在實際設計過程中，應根據項目實際情況調整優化。

關鍵詞：灌溉需水;CROPWAT模型;干旱地區;西非

中圖法分類號：S274 文獻標志碼：A DOI：10.15974/j.cnki.slsdkb.2022.02.011

文章編號：1006 - 0081（2022）02 - 0063 - 05

0 引 言

研究項目所在地區北部為熱帶沙漠氣候，從20世紀70年代開始的干旱導致該地區的自然資源遭到嚴重破壞，農業和畜牧業發展受到限制。水資源時空分布不均衡性導致該地區干旱缺水問題已持續多年，采取工程措施進行合理的水資源調配是解決區域性水資源短缺，有效緩解區域干旱缺水問題，改善區域農業生產條件的重要舉措。灌溉需水量、灌溉設計流量等是確定引水工程規模的重要指標，但該地區經濟相對落后，缺乏一般灌溉設計[1-2]所需的各種觀測數據及實驗資料。CROPWAT模型目前比較廣泛應用于各類作物的需水預測[3-4]、單個作物灌溉制度分析[5-7]及區域蒸發量的分析研究等，鑒于中國已有比較完善的灌溉定額、灌溉制度等，直接采用CROPWAT模型進行工程規劃設計的情況比較少。因此，本文提出利用基于FAO的CROPWAT模型開展灌溉需水預測與規劃設計的思路。本文結合西非地區的工程實際，按照作物需水量、作物系數、降雨、土壤、種植規劃等步驟分析凈灌溉需水量（NIWR），利用灌溉系數計算灌溉毛需水量峰值（GIWRpeak），依據灌溉面積、灌溉工作時間等最終確定灌溉設計流量，可為類似地區引水工程的規劃設計研究提供思路。

1 研究區域概況

1.1 基本概況

項目位于非洲西部地區，屬于沙漠草原氣候，年平均氣溫29 ℃，年均降雨量約314 mm。項目所在國為全球最不發達國家之一，農業是其主要支柱產業之一，該國超過50%的人從事農業。擬建工程從河道引水，利用渠道、湖泊以及疏浚河道作為輸水通道，并在輸水通道末端修建泵站和渠道為規劃灌區進行灌溉，工程疏挖河道約87 km，新建主干渠約150 km。

1.2 灌區現狀及存在問題

項目所在地區為平原地形，平均高程約25 m，平均坡度約0.1‰，高程逐漸向西和西北方向降低。平原受樹狀河網（現狀無水）切割，形成面積多樣的扇形條帶，區域現狀大部分為灌木草原和裸地。地區約90%的人口從事農業和畜牧業，農業主要種植花生、小米等。由于該地區處于撒哈拉沙漠邊緣地帶，屬于沙漠草原氣候，降雨量少、農田產量低。此外土壤質量較差、有機肥料施用少、農藥使用量低、農田設備設施不足，以及耕種技能缺乏也是制約農業發展的主要因素。按照聯合國糧農組織標準，目前該區域糧食無法自給。根據研究，因水資源缺乏，該區域現狀農田較少，呈零星斑塊狀分布。

2 研究方法與參數

研究項目地處非洲西部的干旱地區，缺乏各種基礎水文氣象及灌溉實驗數據等，綜合考慮各種因素，擬采用FAO（世界糧農組織）開發的CROPWAT模型進行需水預測分析，主要研究方法與工作步驟如圖1所示。

2.1 作物需水量

采用FAO-56號文[8]中單作物系數方法計算作物需水量，計算步驟如下：①計算參考作物蒸發蒸騰量ETo;②確定作物的生長發育階段及各發育階段的長度，選定作物在此階段的作物系數Kc（可參考FAO-56），并根據區域的氣象條件調整;③ETo乘以Kc得到ETcrop。

2.1.1 參考作物蒸發蒸騰量ETo

CROPWAT模型是在Penman-Monteith方程的基礎上建立的，該方程以能量平衡和水汽擴散理論為基礎，考慮了作物的生理特征和空氣動力學參數的變化，具有較充分的理論基礎和較高的計算精度，經過多年研究改進，FAO-56號文推薦使用的公式形式為

式中：ETo 為 參考作物蒸發蒸騰量，mm/d;[Δ]為飽和水汽壓-溫度曲線斜率，KPa/℃;Rn 為冠層凈輻射，MJ/（m2 ·d）;G為土壤熱通量，MJ/（m2 ·d）;T為日平均氣溫，℃;u2 為 2 m高度處的風速，m/s;es為飽和水汽壓，kPa;en為實際水汽壓，kPa。

CROPWAT模型計算需要的參數數據包括溫度、濕度、日照時間、蒸發、風，上述氣象數據依據FAO提供的CLIMWAT模型導出，再輸入CROPWAT模型中，項目所在區域的ETo計算結果如表1所示。

2.1.2 作物系數

作物系數Kc反映作物和參考作物之間需水量的差異，體現了作物本身的生物學特性，主要取決于作物冠層的生長發育。作物系數可用一個系數來綜合反映，也可用兩個系數分別來描述蒸發和蒸騰的影響，即所謂的單作物系數和雙作物系數，一般設計可采用單作物系數，FAO建議將作物的生育期劃分為4個階段：初期、發育期、中期和后期，本項目主要依據FAO-56號文提供的作物系數，并適當調整。

2.1.3 作物需水量

作物需水量（ETcrop，mm/d）定義為作物在整個蒸散過程中所需的水的深度，可用作物系數（Kc）乘以參考作物蒸發蒸騰量（ETo）得出，具體公式如下：

式中： ETo 為參考蒸發蒸騰量，mm;Kc 為計劃耕作作物的作物系數，取自FAO-56號文。

將作物系數輸入CROPWAT模型，即可計算出作物的需水量ETcrop，項目典型作物（玉米）的作物需水量（ETcrop）計算結果如表2所示。

2.2 降雨參數

國外項目一般采用P=80%頻率的降雨進行灌溉設計，依據當地1962～2016年降雨統計資料，按經驗頻率公式，采用頻率曲線分析P=80%的年均降雨量，并利用式（3）計算設計枯水年份的月均降雨量，作為程序輸入數據。經分析，項目所在地區年平均降雨量為314.8mm，P=80%的年降雨量為234.4 mm。設計枯水年（P=80%）降雨計算結果見圖2。

式中：Piav 為月均降雨量，mm;Pidry 為設計枯水年月均降雨量，mm; Pav 為年平均降雨量，mm;Pdry 為設計枯水年平均降雨量，mm，P=80%。

有效降雨量（effective rainfall），即考慮了地表徑流和深層滲透的降雨損失后可供作物有效利用的降雨量。CROPWAT模型中提供了多種有效降雨量分析方法，本項目有效降雨采用FAO的AGLW公式進行計算，經計算有效年降雨量為102.6 mm，模型輸出計算結果如表3所示。

式中：Pe為有效降雨量，mm;Ptot 為平均月降雨量，mm。

2.3 土壤參數

2.3.1 土壤情況

項目所在河谷地區為沙漠草原氣候環境，土壤中松散層厚度超過0.6 m，且不存在土質差、卵石多、板結、含鹽、風蝕或水蝕敏感、較大坡度等限制因素，評估認為該區土壤適合種植及灌溉。對河谷兩側各20 km的帶狀區域進行研究，發現大部分的研究區域都是潛在的適合耕種和灌溉的土壤，初步勘測了區域適合、一般、不適合發展農業灌溉的土地范圍，河谷原河道兩岸區域大部分土地都適合耕種和灌溉。

2.3.2 土壤參數

CROPWAT模型土壤數據參數分為水稻和非水稻作物兩種，本項目主要作物為非水稻作物，土壤數據模塊參數包括，總可用水量 （TAW）、最大滲透率（maximum infiltration rate）、最大生根深（maximum rooting depth）、初始土壤水分（initial soil moisture depletion），依據項目地區土壤情況選定。

2.4 種植規劃

收集地區已有作物及相關研究資料，初步擬定項目規劃灌區的主要種植作物為玉米、土豆、甘蔗、番茄、洋蔥等，主要作物種植計劃如表4所示。

3 結果分析與討論

依據上述氣象、降雨、土壤以及種植規劃等參數，CROPWAT模型可以計算出項目凈灌溉需水量（NIWR），利用模型計算結果進一步分析項目灌溉設計流量。

3.1 凈灌溉需水量

ETcrop 和有效降雨量（Pe） 的差額為凈灌溉需水量（NIWR），計算見式（6）所示，通過CROPWAT模型可自動輸出計算結果，如表5所示。

式中：ETcrop為作物需水量，mm/d;Pe 為有效降雨量，mm。

3.2 毛灌溉需水量

毛灌溉需水量（GIWR）即滿足作物的蒸發蒸騰需求，以及不可避免的損失水量和滲漏水量，計算公式如下：

式中：Ep 為項目灌溉系數，%。

Ep主要考慮灌溉方法、滲流損失、渠道和相關結構的損失以及取水口與田間的損失，本項目灌溉系數取為60%。

3.3 灌溉設計流量

依據CROPWAT程序計算結果，分析各月中的峰值灌溉需水量作為灌溉設計流量（m3/s），計算公式如下：

式中：A 為灌溉面積，hm2;GIWRpeak 為 灌溉毛需水量峰值，mm/d;h 為每天的工作小時數（項目取16 h）;v為每周的工作日（項目取6 d），若用月天數，則v為每月工作日。

結合模型計算結果，項目灌溉凈需水量峰值為2.6 mm/d，據此估算灌溉毛需水量峰值GIWRpeak為4.33 mm/d，項目灌溉面積24 348 hm2，分析項目設計流量Q設=21.4 m3/s。

3.4 討 論

CROPWAT模型有效降雨公式有4種，有效降雨計算結果為102.6～207.7 mm，對于作物需水量有一定影響，本文選擇偏保守的公式;由于項目灌溉高峰月在旱季，有效降雨量的變化對于最終設計流量影響不大。種植計劃的變化特別是作物的起始種植時間的變化對灌溉毛需水量峰值有一定影響，需要依據項目實際情況調整。灌溉設計流量對于灌溉系數、工作時間的變化較為敏感，若灌溉系數在0.55～0.65變化，則相應設計流量變化幅度約為15%;工作小時數若在16～24 h間變化，相應設計流量變化幅度約為33%，規劃設計過程應合理選取各類參數。

4 結 論

（1） 國際項目工程灌溉規劃設計面臨各種困難，特別是非洲等國別農業發展需求迫切，但難以獲取各種農業灌溉的試驗資料、土壤以及作物灌溉制度等資料，為工程規劃設計帶來了諸多實際困難。

（2） CROPWAT是FAO開發的軟件，可用于作物需水量預測、灌溉工程設計以及灌溉制度分析等工作，對于在資料缺乏地區開展規劃設計工作具有重要意義，工作過程中發現很多國外咨詢機構也經常采用該模型開展規劃設計工作。

（3） CROPWAT程序中的氣象、作物系數、降雨資料以及土壤資料數據是重要的輸入數據，應盡可能地依據當地實際情況獲取或者調整;種植計劃變化對于灌溉毛需水量峰值有一定影響，而灌溉設計流量對于灌溉系數、工作時間的變化較為敏感，在實際設計過程中，應根據項目實際情況調整優化。

參考文獻：

[1] 于濤.定額法計算灌溉需水量[J].現代農業科技，2013（14）：202-202.

[2] 索麗生，劉寧. 水工設計手冊第二冊[M]. 北京：中國水利水電出版社，2014.

[3] 段琪彩，張雷，黃英，等.? CROPWAT模型在滇中南部灌水量模擬中的應用研究[J].節水灌溉，2020（6）：81-86.

[4] 趙福民，王新華，萬大紅，等. 基于CROPWAT的昆明主要農作物虛擬水分析[J]. 安徽農業科學，2017，45（25）：218-220.

[5] 徐冰，湯鵬程，李奇，等.? 基于CROPWAT模型的拉薩地區燕麥優化灌溉制度研究[J].干旱地區農業研究，2015（6）：35-39，183.

[6] 李紅芳，王會肖，楊雅雪，等. 基于WOFOST與CROPWAT模型的松嫩平原春玉米灌溉制度研究[J]. 北京師范大學學報（自然科學版），2020，56（3）：424-435.

[7] 王忠波，王欣亮，聶堂哲. 基于CROPWAT模型水稻需水量及灌溉制度研究[J]. 東北農業大學學報，2019，50（7）：87-96.

[8] RICHARD G， ALLEN，LUIS S. Crop evapotranspiration? guidelines for computing crop water requirements—FAO irrigation and drainage paper 56[M]. Rome：FAO，1998.

（編輯：江 文）