□聶義泊 □楊 杰

（1河南省水利勘測設(shè)計研究有限公司；2河南黃河建設(shè)工程有限公司）

0 引言

干旱缺水問題已經(jīng)成為危機世界糧食安全，人類健康和自然生態(tài)系統(tǒng)的重大問題。據(jù)估計，按目前人口增長速度，到2025年世界人口將增加到78.50億，糧食產(chǎn)量必須翻一番才能滿足需求，為提高糧食產(chǎn)量，農(nóng)業(yè)用水將增加14%，即需要20000億m3的灌溉用水。如何才能提供足夠的糧食生產(chǎn)用水已成為全球范圍內(nèi)所面臨的最嚴峻挑戰(zhàn)。因此，無論是發(fā)達國家還是發(fā)展中國家都把發(fā)展農(nóng)業(yè)節(jié)水作為解決這一挑戰(zhàn)的首要途徑。

水稻是我國的主要糧食作物之一，全國水稻產(chǎn)量約占糧食作物總產(chǎn)量的40%，同時，它也是主要耗水大戶，占農(nóng)業(yè)用水量的65%以上（茆智，2002；茆智，崔遠來，2003；彭世彰，2011）。我國水稻節(jié)水灌溉技術(shù)自20世紀90年代以來得到了大面積的推廣應(yīng)用（彭世彰，徐俊增，2011）。水稻節(jié)水灌溉的特點是通過設(shè)置較低的灌水上限、下限，使得水稻在生長期間某些時間田面無水層或土壤含水量低于其飽和含水量，總的灌溉用水量明顯低于常規(guī)淹灌（彭世彰等，2005；2011）。

對水稻生產(chǎn)來說，研究節(jié)水灌溉條件下土壤水分運動規(guī)律，對稻田土壤水分狀況進行調(diào)控，使之在非充分灌溉與充分灌溉之間合理轉(zhuǎn)換，進一步提高作物水分利用效率，是一個值得探究的科學(xué)問題，其成果可以直接用來確定灌水下限、上限等灌溉要素，是對節(jié)水灌溉理論的重大補充。

通過建立稻田土壤水分運動數(shù)學(xué)模型，分析土壤水分運動規(guī)律是進行水分調(diào)控的前提，是完善節(jié)水灌溉技術(shù)的理論基礎(chǔ)。為此，需要解決兩個問題：一是，土壤水分運動基本方程的選取；二是，作物根系吸水模型的建立。

1 土壤水分運動基本方程

從土壤水動力學(xué)的角度來看，現(xiàn)行的各種水稻節(jié)水灌溉模式，無論其灌水上限、下限有何不同，灌水時間如何確定，其土壤水分運動均可以用統(tǒng)一的數(shù)學(xué)模型來描述。土壤水分運動基本方程的選取是關(guān)鍵。

國內(nèi)外學(xué)者在進行農(nóng)田灌溉排水，地下水補給等問題的研究中提出了具有不同特點和用途的土壤水分運動方程。常見的有Kostiakov模型，Kostiakov-Leiws模型，Philip模型，Horton模型，Green-Ampt模型等。但是，這些土壤水分入滲模型均未考慮下墊面、土地利用等因素的影響，導(dǎo)致利用起來誤差較大，適用條件和范圍、主要參數(shù)的確定、理論的近似解等方面均存在問題。在實際應(yīng)用時，需要對其進行補充和完善。

以達西定律為基礎(chǔ)的Richards方程在研究作物生長條件下的土壤水分運動中被廣泛應(yīng)用。Constable等人（1981）利用該方程研究了半濕潤環(huán)境中的作物灌溉問題。Belmans等人（1983）建立了作物生長情況下的土壤水分模擬模型。Kiefer EM（1993）提出了一個模擬非均質(zhì)土不飽和滲流的隨機概念模型。雷志棟，楊詩秀等人（1985）在我國較早地對非飽和均質(zhì)土壤的水分運動進行了多種條件下的數(shù)值模擬研究。申雙和，周英等人（1992，1993）開展了農(nóng)田土壤水分預(yù)測預(yù)報的工作。高煥文等人（1996）建立了保護性耕地土壤水分模型。

Richards方程常見形式為：

以含水量為因變量，可寫為：

式中S(z,t)為該方程的源匯項，即作物根系吸水模型。Richards方程計算和應(yīng)用前，需要確定土壤水力學(xué)特性參數(shù)。在研究有作物生長影響的土壤水分運動時，還需要確定作物根系吸水模型。

2 作物根系吸水模型

通過研究作物根系吸水機理，建立根系吸水模型，對根系吸水過程進行定量描述，一直是土壤物理學(xué)和農(nóng)業(yè)水利學(xué)等學(xué)科研究的熱點。

Gardner(1960)首先開展了這方面的工作，提出了單根吸水的理論模型，假定單根為無限長、半徑均勻并具有均勻吸水特性的圓柱體，植物根系和土壤水吸力分布都是均勻的，它研究流入典型單根的徑向流。Hainsworth，Aylmore經(jīng)過研究發(fā)現(xiàn)，Gardner假定與事實不符（1989）。Aura提出了單根非均勻吸水模型（1996），該模型是否與田間情況一致還無法驗證。Molz于1970年與Remson提出了吸水速率與蒸騰速率、有效根密度和土壤水分擴散率之積成正比的根系吸水模型；1976年，又將根系吸水在土壤中水分向根表面的流動和水分在根組織內(nèi)的流動聯(lián)系起來，從而考慮了根的水力特性，提出了根系吸水的土壤一維水流運動模型。Nimah和Hanks（1973）考慮溶質(zhì)的影響和植物根異管傳導(dǎo)水分的內(nèi)摩阻力，修正了Molz-Remson模型。Hillel（1976）等人在Gardner等人研究的基礎(chǔ)上，考慮到土壤及根系對水流的阻力，利用根系吸水速率與土壤和根系水勢之差關(guān)系建立起來的模型。Herklrath（1977）在模型的基礎(chǔ)上加進飽和度作為根土接觸的校正因子，建立吸水模型。Feddes（1978）認為根系吸水率是水勢差函數(shù)，首先將根系吸水與作物蒸騰聯(lián)系起來，將蒸騰量在根系層土壤剖面上的土水勢建立的根系吸水模型。VanGenuchten在1984年對α(h)進一步修正，提出了非線性的定義。Chandra和Amaresh（1996）改進了前人模型的不足之處，建立了作物蒸騰量參與的非線性根系吸水模型。

我國學(xué)者對此也進行了大量的探索（許迪，1997；黃冠華等，1995；趙成義等，1999；郭慶榮等，1999）。邵明安（1987）根據(jù)植物根系吸水的物理過程，提出了一個能反映根系吸水機理的宏觀數(shù)學(xué)模型?？到B忠等（1992）從作物根系吸水機理出發(fā)，用動態(tài)模擬的方法估算作物根系吸水速率，分別建立了冬小麥、夏玉米的根系吸水模型。龔道枝，康紹忠，張建華等（2004；2010）系統(tǒng)地研究了果樹的二維根系吸水特性，根據(jù)根長密度資料建立的根系吸水模型?；⒛憽ね埋R爾白（1999）以土壤水動力學(xué)理論為基礎(chǔ)，綜合考慮作物蒸騰，有效重量根密度分布以及土壤含水率的影響建立了二維根系吸水模型。

經(jīng)過眾多專家學(xué)者對根系吸水的多年研究，其吸水機理逐漸清晰，模型建立方法、途徑逐漸明朗，所建模型的適用性更加廣泛。但是，這些研究多數(shù)限于旱作物，在旱田土壤水分環(huán)境中，圍繞旱作物根系吸水進行機理分析與理論建模。由于水稻與旱作物生理特性的差異及田間土壤水分狀況的不同，尤其是在節(jié)水灌溉形成的稻田土壤水分環(huán)境情況下，現(xiàn)有的旱作物根系吸水模型能否直接用來解決水稻根系吸水問題，以及如何建立水稻根系吸水模型，都值得進一步探討。

3 存在的問題

3.1 構(gòu)建準確描述節(jié)水灌溉稻田土壤水分環(huán)境中根系吸水特性的吸水模型

水稻根系吸水模型作為Richards方程的源匯項，與水稻根系吸水能力有關(guān)。根系吸水能力與土壤水分環(huán)境密切相關(guān)，尤其是在節(jié)水灌溉條件下，土壤水分在飽和非飽和之間變化劇烈，干濕交替頻繁，水稻水分生理活動將對此做出相應(yīng)的調(diào)整，改變根系吸水能力。準確的水稻根系吸水模型對于建立正確的土壤水分運動基本方程，進而描述節(jié)水灌溉條件下的土壤水分運動過程至關(guān)重要。

3.2 建立節(jié)水灌溉稻田土壤水分運動數(shù)學(xué)模型

研究節(jié)水灌溉稻田土壤水分運動的邊界條件及初始條件，通過建立土壤水分運動數(shù)學(xué)模型，研究節(jié)水灌溉土壤水分運動規(guī)律，對于夯實節(jié)水灌溉理論具有重要意義，可以為建立相應(yīng)的土壤水分調(diào)控模式奠定堅實的理論基礎(chǔ)。

[1]彭世彰,徐俊增.水稻控制灌溉理論與技術(shù)[M].南京:河海大學(xué)出版社,2011:8-13.

[2]Sun H Y,Shen Y J,Yu Q,Flerchinger G N,Zhang Y Q,Liu C M,Zhang X Y.Effect of precipitation change on water balance and WUE of the winter wheat-summer maize rotation in the Nor th China Plain.Agricul tural Water Management,2010,97:1139-1145.

[3]郭相平,袁靜,郭楓,陳治平.水稻蓄水-控灌技術(shù)初探[J].農(nóng)業(yè)工程學(xué)報,2009,4:70-73.

[4]陳朱葉,郭相平,姚俊琪.水稻蓄水控灌的節(jié)水效應(yīng)[J].河海大學(xué)學(xué)報(自然科學(xué)版),2011,4:426-430.

[5]郭以明,郭相平,樊峻江,張曉柳.蓄水控灌模式對水稻產(chǎn)量和水分生產(chǎn)效率的影響[J].灌溉排水學(xué)報,2010,3:61-63+73.

[6]解河海,馮杰,馮青,等.考慮大孔隙的土壤入滲模型[J].水利水電科技進展,2011,31(5):35-38.

[7]彭振陽,黃介生,伍靖偉,等.基于分層假設(shè)的Green-Ampt模型改進[J].水科學(xué)進展,2012,23(1):59-66.

[8]郭向紅,孫西歡,馬娟娟,等.不同入滲水頭條件下的Green-Ampt模型[J].農(nóng)業(yè)工程學(xué)報,2010,26(3):64-68.

[9]毛麗麗,雷廷武.用修正的Green-Ampt模型確定土壤入滲性能的速算方法[J].農(nóng)業(yè)工程學(xué)報,2010，26(12):53-57.

[10]史曉楠,王全九,巨龍.微咸水入滲條件下Phil ip模型與Green-Ampt模型參數(shù)的對比分析[J].土壤學(xué)報,2007,44(2):360-363.