作物生長田間水分平衡的系統(tǒng)模擬

1、作物生長田間水分平衡的系統(tǒng)模擬 06-02-21 10:33:00 作者：呂 軍編輯：studa9ngns摘要：本文以作物生長模型MACROS(L1D)為基礎，建立了土壤水分運動和作物生長耦合關系綜合模擬系統(tǒng). 著重介紹了SPAC體系水分運動子模型，并對土壤水分運動模擬的結果和各類參數(shù)的靈敏度等進行了分析.最后，應用該模擬系統(tǒng)對南方多雨地區(qū)冬小麥田土壤水分平衡和土壤漬害的發(fā)生過程進行了模擬分析. 關鍵詞 關鍵詞：SPAC水分運動 水分與作物生長耦合模型 水分平衡 漬害 地表排水 近年來數(shù)值模擬方法在農田土壤水分運動研究中被廣泛采用1,4，并逐步從實驗室引向田間5,6，也考慮了農田作物狀況對土壤

2、水分運動的影響2,3,7,9.但是，國內大多數(shù)水分運動模型沒有涉及土壤水分條件對作物生長過程的動態(tài)影響，因而不能很好地反映土壤水分運動與作物生長動態(tài)之間所具有的強烈的反饋和耦合作用，實際上，不良水分條件限制作物生長，而作物生長狀況的改變又對農田水分平衡過程產生重要影響.同時，作物的不同種植方式也將影響農田水分平衡狀況.在作物生長動態(tài)模型MACROS(L1D)8和田間水分平衡模型的基礎上，建立并引進了農田水分平衡與作物生長動態(tài)耦合子模型，本文著重討論水分運動子模型的模擬方法，分析模型模擬結果的可靠性及其對主要參數(shù)的敏感性，并應用該模型研究了多雨地區(qū)冬小麥生長田間土壤水分動態(tài)，討論多雨地區(qū)土壤漬、

3、澇的發(fā)生和消退過程.1土壤水分運動數(shù)學模型以近作物冠層的大氣為上邊界，土體一米深處為下邊界，構成作物生長田間農田水分平衡系統(tǒng).作物生長模型MACROS(L1D)以1d(24h)為時間步長，在作物出苗后開始模擬，并向土壤水分平衡模型提供所需作物參數(shù).而土壤水分運動模型則以可變時間步長(0.1-0.001d)進行模擬.因此，土壤水分狀態(tài)經(jīng)24h循環(huán)計算后才取值與作物發(fā)生相互作用.當不良土壤水分條件發(fā)生時，作物生長受到限制，從而也改變了向土壤水分平衡模型提供的有關作物參數(shù)，這樣就構成了土壤水分平衡與作物生長動耦合關系.模型中，僅考慮作物特性、土壤水分特性和氣象要素，并認為其它條件均滿足作物生長要求.

4、1.1水分液態(tài)流 土壤水分液態(tài)流為有根系吸水條件下的垂直一維流，即/t=-q/Z-Sr,其中q為水流通量，q=-K(h)h/Z+K(h),因此有：C(h)(h)/(t)=()/(Z)K(h)(h)/(Z)-K(h)-Sw, (1)式中 h(cmH2O)為壓力水頭，K(h)為導水率(cm d-1-1), C=d/dh為比水容量，Z為深度座標(cm)，t為時間(d)，Sw為根系吸水函數(shù).田間條件下，土壤水分運動有多種邊界條件，而它們又反復交替變化.本模型有：初始條件： h=hi, Z0, t=0.上邊界條件：(1)表層土壤接近飽和或有極薄積水時，h=0, Z=0, t0. (2)土表向土壤供水，供

5、水強底低于入滲率，K(h)(h)/(Z)-K(h)=-i, Z=0, t0. (3)土表蒸發(fā)而土壤供水充足，K(h)(h)/(Z)-K(h)=E, Z=0, t0. (4)土壤蒸發(fā)而土體向土表供水速率低于蒸發(fā)速率，土表趨于風干狀態(tài)，h=had, Z=0, t0; 其中腳標ad表示風干土.下邊界條件：下邊界條件主要考慮地下水位埋深ZW(m)的影響(1) ZW3.0; h=hi, Z=Zmax, t0. (2)3.0ZW1.0; h=ZW-Zmax, Z=Zmax, t0. (3)ZW1.0; h=0, Z=ZW, t0. 其中，當ZW1.0m時，地下水位埋深即為模擬系統(tǒng)下邊界，流入地下水位中的水

6、亦即為土壤水滲漏損失.土壤中發(fā)生飽和流時， /t=0, 有：q=-Ki(hi+1-hi)/(Zi+1-Zi)i,(h2-h1)+(h3-h2)+(hn+1-hn)=hn+1-h1, (2)i=n, 為土壤層次數(shù)，得到n+1個線性方程求解.對土壤中的非飽和流，考慮到土壤水分運動模型與作物生長動態(tài)模型的綜合模擬計算量較大，故采用基質流勢(Matric flow potential)8,9的概念求解水流方程.定義U為基質流勢U-hK(h)dh,則水流通量方程可改寫為：q=-(dU)/(dZ)+K(h) ,(3)假設兩相鄰土層中心的流量隨深度的分布不變，微薄土層中心之間的水勢變化為線性的，則式(3)可

7、以導出水流通量方程的積分形式為：(4)當水流發(fā)生在非均質土層之間時，式(4)中的積分項由下式修正.(5)作物根系吸水函數(shù)Sw, 先由作物冠層蒸騰量E確定其總量，當土壤供水充足時，各土層中根系吸水量Swi，根據(jù)已確定的根系總量和由實驗測定的根系在土壤剖面上的分布比例i計算，即Swi=iE .(6)不良土壤水分條件發(fā)生時，根據(jù)旱害、漬害程度對式(6)進行修正，修正參數(shù)由實驗確定9.1.2水分汽態(tài)流Ep=ESR+E (7)其中(8)式(8)中，SR為輻射吸收量，為水的汽化潛熱，S為飽和水汽壓隨溫度變化曲線的斜率，P和*分別為空氣干燥力和表觀溫度計常數(shù)項，它們又分別為：(9)式中es和ea分別為飽和水

8、汽壓和水汽壓實測值，Cp為空氣的體積熱容量，為濕度計常數(shù)，而h、b和s分別為邊界層熱交換阻力、邊界層水汽交換阻力和氣孔阻力.式(8)中，ESR和E受作物冠層葉面積狀況的影響，因此作物潛在蒸騰速率Ecp,應由葉面積系數(shù)ALV進行校正8，則式(8)可改寫為：(10)式中ES為土壤蒸發(fā)率，bL和bS分別為作物遮陰修正系數(shù)和土壤濕度修正系數(shù)，即(12)式中m和n分別經(jīng)驗系數(shù)，分別與作物種類和土壤性質有關，ALV為包括死葉和活葉的總葉面積系數(shù)，a和f則分別為風干土壤濕度和田間持水量.考慮到實際土壤蒸發(fā)并不都發(fā)生在一個峰面上，則應由各土層蒸發(fā)貢獻率mi將ES值分配到各土層中去：(13)式中腳標i為土壤層次，d為土層厚度，L為土層深度，K為經(jīng)驗系數(shù).此外，為避免地下側滲流的引入而造成模型過于復雜，地下水埋深按實際觀測值輸入.1.3模型軟件及其輸入和輸出整個模型模擬軟件用CSMP語言編制(連續(xù)系統(tǒng)模擬程序語言).模型所需輸入資料包括作物資料(品種特性和生長特性資料)、氣象資料(日最高和最低