土壤水特征曲線

1、 研究生課程論文封面課程名稱 土壤水動(dòng)力學(xué) 教師姓名 研究生姓名 研究生學(xué)號(hào) 研究生專業(yè) 所在院系 類別: 日期: 2012 年 1月 7 日 評(píng) 語對(duì)課程論文的評(píng)語:平時(shí)成績：課程論文成績：總 成 績：評(píng)閱人簽名：注：1、無評(píng)閱人簽名成績無效；2、必須用鋼筆或圓珠筆批閱，用鉛筆閱卷無效；3、如有平時(shí)成績，必須在上面評(píng)分表中標(biāo)出，并計(jì)算入總成績。水分特征曲線測(cè)定實(shí)驗(yàn)報(bào)告1 實(shí)驗(yàn)的目的要求 理解水分特征曲線的含義，掌握水分特征曲線的測(cè)定方法，以及比較不同土壤水分特征曲線的特點(diǎn)。2 實(shí)驗(yàn)的原理土壤水的基質(zhì)勢(shì)（或土壤吸力）與土壤含水量之間的關(guān)系曲線稱為土壤水分特征曲線或土壤持水曲線（soil wat

2、er retention function）。土壤水分特征曲線表示土壤水的能量和數(shù)量之間的關(guān)系，是研究土壤水分的保持和運(yùn)動(dòng)所用到的反映土壤水分基本特性的曲線。各種土壤的水分特征曲線均需由實(shí)驗(yàn)測(cè)定。水分特征曲線儀主要由陶土頭、集氣管、壓力傳導(dǎo)管、水銀測(cè)壓計(jì)（由玻璃管和水銀槽組成）、觀測(cè)板以及樣品容器組成，其結(jié)構(gòu)如圖1所示。圖1 水分特征曲線儀結(jié)構(gòu)圖1.樣品容器；2.陶土頭；3.集氣管；4.壓力傳導(dǎo)管；5.水銀測(cè)壓計(jì)；6.觀測(cè)板；7.水銀槽陶土頭是儀器的傳感部件，由具有均勻微細(xì)孔隙的陶土材料制成，當(dāng)儀器內(nèi)充滿水使陶土頭被水飽和時(shí)，陶土頭管壁就形成張力相當(dāng)大的一層水膜，陶土頭與土壤充分接觸后，土壤水

3、與其內(nèi)部的水體通過陶土頭建立了水力聯(lián)系，在一定的壓差范圍內(nèi)，水分和溶質(zhì)可以通過陶土頭管壁，而氣體則不能通過，即所謂透水不透氣。因此，如果陶土頭內(nèi)外之間存在壓力差，水分就會(huì)發(fā)生運(yùn)動(dòng)，直至內(nèi)外壓力達(dá)到平衡為止。這時(shí)，通過水銀壓力表測(cè)定的負(fù)壓值就是陶土頭所在位置土壤水的基質(zhì)勢(shì)。陶土頭所在位置的壓力水頭（基質(zhì)勢(shì)或負(fù)壓）的計(jì)算公式為：式中h為壓力水頭，hm為壓力表中水銀柱高度（以水銀槽水銀液面為基準(zhǔn)面），hm是水銀槽液面到陶土頭中心位置的垂直距離。3 實(shí)驗(yàn)的步驟1） 準(zhǔn)備無氣水若干、與陶土頭直徑相當(dāng)?shù)你@頭一個(gè)、注射器一個(gè)、刻度尺一把、透水石若干、電子天平一臺(tái)、記錄表格若干及其它輔助器材如燒杯、橡皮、鉛筆

4、等。2） 按圖1組裝好儀器，檢查儀器是否漏氣；用無氣水浸泡陶土頭（無氣水要淹沒整個(gè)陶土頭）并注入無氣水與其相連的連通管中，浸泡時(shí)間要一般需24小時(shí)以上，以除去其表面的氣泡。在確認(rèn)儀器沒有問題以及沒有氣泡存在時(shí)，重新?lián)Q入無氣水（換水后，要確保連通管中無氣泡），用電子天平稱量此時(shí)的儀器重量，并做記錄。3） 裝樣：裝樣時(shí)，保證陶土頭與土樣接觸良好，同時(shí)檢查各部件的接口處是否密封完好。4） 飽和土樣，把盛土樣的儀器放在透水石上，水面最好與透水石持平，時(shí)間一般在24個(gè)小時(shí)以上。稱量飽和后的重量（包括儀器及土樣），并做好記錄。5） 一定時(shí)間后，讀取觀測(cè)板水銀柱凸液面的高度，同時(shí)稱重（儀器及土樣總重）。稱重

5、時(shí)注意不要碰撞儀器，以免漏氣或陶土頭與樣品接觸不良造成實(shí)驗(yàn)失敗。開始測(cè)量時(shí)，要密切注意水銀柱讀數(shù)的變化，不同土樣其變化差異較大；應(yīng)持續(xù)觀察，在有較大變化時(shí)測(cè)量一次；在掌握其變化規(guī)律后，視土面蒸發(fā)變化情況定時(shí)測(cè)量；每天可測(cè)一至兩次。隨著時(shí)間的延長，壓力傳導(dǎo)管以及集氣管中可能存在少量氣泡，在氣泡沒有連通形成斷點(diǎn)時(shí)，不必?fù)Q無氣水；在氣泡連通形成斷點(diǎn)后，必須更換儀器中的水。6） 在測(cè)量后期（土壤負(fù)壓在760cm水柱高度左右），如果更換無氣水后，在下次測(cè)量時(shí)，又形成很多斷點(diǎn)，這時(shí)應(yīng)該停止測(cè)量，即此時(shí)的負(fù)壓已超出了該儀器的測(cè)量范圍。停止測(cè)量后，取儀器中的土樣，放在烘箱中，105oC烘24小時(shí)，稱量儀器與干

6、土樣的重量，確定土樣干容重。7） 計(jì)算土壤的體積含水量（重量含水量*干容重）與相應(yīng)的土壤負(fù)壓值，根據(jù)van Genuchten公式，利用statistic非線性程序包或RETC軟件進(jìn)行參數(shù)擬合，即可獲得此時(shí)的水分特征曲線。式中：為飽和含水量；為凋萎含水量；h代表負(fù)壓；代表體積含水量；、m、n為待定系數(shù)，m=1-1/n，由土壤的性質(zhì)確定。4 實(shí)驗(yàn)數(shù)據(jù)記錄于2011年11月26日開始本次實(shí)驗(yàn)；11月26日至12月12日每天記錄并觀測(cè)板水銀柱凸液面的高度，同時(shí)稱重，數(shù)據(jù)記錄如表所示；2012年1月5日，完成實(shí)驗(yàn)土樣的烘干稱重，以及實(shí)驗(yàn)數(shù)據(jù)的初步整理工作。表1 實(shí)驗(yàn)記錄表儀器+飽水土重（g）5235其

7、它儀器重（除環(huán)刀（g)2103.3環(huán)刀+飽水土重（g)3131.7環(huán)刀+干重（g）1931.1環(huán)刀重（g)1477.6環(huán)刀直徑（cm）14.75環(huán)刀高（cm）5表2 水分特征曲線測(cè)量記錄表水銀柱讀數(shù)（cm）土壤負(fù)壓（cm）儀器+水+土樣重（g）土壤含水量（g）土壤重量含水量（%）土壤體積含水量（%）2.2-25.925216.1361.80.2841 0.4233 2.5-29.75210.2355.90.2795 0.4164 2.78-33.2285196.1341.80.2684 0.3999 3.1-37.265192.7338.40.2657 0.3960 4.3-52.385174

8、.33200.2513 0.3744 4.65-56.795168.5314.20.2467 0.3676 7.3-90.185147292.70.2299 0.3425 7.8-96.485142.5288.20.2263 0.3372 10-124.25129.7275.40.2163 0.3222 10.4-129.245126.4272.10.2137 0.3184 13.7-170.825114.8260.50.2046 0.3048 14.05-175.235113.5259.20.2035 0.3033 17.7-221.225103.8249.50.1959 0.2919 18

9、.32-229.0325102.2247.90.1947 0.2901 19.6-245.165098.8244.50.1920 0.2861 22.9-286.745092.8238.50.1873 0.2791 23.2-290.525091.9237.60.1866 0.2780 24.2-303.125089.9235.60.1850 0.2757 24.4-305.645089.1234.80.1844 0.2747 25.9-324.545087.1232.80.1828 0.2724 32.4-406.445077.1222.80.1750 0.2607 37-464.45071

10、216.70.1702 0.2536 39.5-495.95069.8215.50.1692 0.2522 48.2-605.525059.7205.40.1613 0.2403 56-703.85046.7192.40.1511 0.2251 66.6-837.365036.7182.40.1432 0.2134 71.5-899.15019.2164.90.1295 0.1929 73.5-924.34997.5143.20.1125 0.1676 185 實(shí)驗(yàn)數(shù)據(jù)處理將表2的土壤吸力和體積含水量輸入到RETC軟件中，參數(shù)設(shè)置如表3所示： 表3 RETC軟件參數(shù)設(shè)置設(shè)置項(xiàng)設(shè)置內(nèi)容Type

11、of Problem選擇Retention data onlyScale UnitsLength Units選擇cm ,Time Units選擇DaysType of Retention/Conductivity Model每個(gè)模型都進(jìn)行一次模擬，選取R值最大的一個(gè)進(jìn)行模擬預(yù)測(cè)Soil Hydraulic ParametersFitted參數(shù)勾選Qr QS Alpha n;Soil Catalog for Initial 選擇siltRetentionPressure對(duì)應(yīng)土壤吸力，Theta對(duì)應(yīng)體積含水量，注意體積含水量不能填百分制的，否則擬合溢出，無結(jié)果。Run RETC選擇不同的土壤水分特

12、征曲線函數(shù)進(jìn)行參數(shù)擬合，各特征曲線對(duì)應(yīng)的相關(guān)系數(shù)R和參數(shù)擬合結(jié)果如表4所示：表4 不同土壤特征曲線函數(shù)對(duì)應(yīng)的R值和參數(shù)擬合結(jié)果Retention Curve ModelConductivityRrsn van Genuchten 1980，variable m and n Mualem0.975435590.05423 0.458830.015591.00500van Genuchten 1980，variable m and nBurdine0.826068320.147730.43655 0.016292.00500van Genuchten 1980，m = 1 - 1/nMualem0

13、.975543430.453110.024961.24858van Genuchten 1980， m = 1 - 2/nBurdine0.972072260.430860.027982.22046Brooks and Corey 1964Mualem0.965753770.408900.02569 0.20543Brooks and Corey 1964Burdine0.965753770.408900.025690.20543Kosugi 1996 (log normal distribution model) Mualem0.977453690.8049934.95205.04147Du

14、al-porosity model Mualem0.974513250.468560.103611.143656 研究成果及結(jié)論分析表4，本應(yīng)該選擇R值最大的log normal distribution model ，但是為了方便計(jì)算，并且van Genuchten模型與log normal distribution模型R相差不大，因此采用van Genuchten模型進(jìn)行擬合，土壤水分特征曲線如圖2所示：圖2 土壤水分特征擬合曲線參數(shù)擬合結(jié)果：飽和含水量s=0.45311，凋萎含水量r=0；=0.02496；n=1.24858；m=1-1/n=1-1/1.24858=0.1991，故負(fù)壓h

15、與體積含水量的關(guān)系曲線方程為：由圖2可看出，實(shí)測(cè)值與擬合曲線吻合程度很高，我們認(rèn)為模型的選擇和參數(shù)的擬合結(jié)果是合理的。通過此模型和相關(guān)參數(shù)，可以用來衡量該地區(qū)土壤水分的能量水平，估計(jì)作物生長的有效性，決定是否采取補(bǔ)水措施，推算土壤中孔隙分布和比水容量，為研究土壤釋水性能指標(biāo)和諧調(diào)水氣提供理論依據(jù)。（實(shí)驗(yàn)二）應(yīng)用HYDRUS-1D模型-評(píng)價(jià)欒城地區(qū)淺層土壤水分運(yùn)移1 引言土壤水是四水轉(zhuǎn)化的中樞,無論是灌溉水、還是天然降水,都要轉(zhuǎn)化為土壤水后才能被作物根系吸收。人們對(duì)土壤水的研究已有很長的歷史,但把土壤水作為一種資源加以研究的時(shí)間并不長。自20世紀(jì)70年代蘇聯(lián)地理學(xué)家李沃維奇首次使用“土壤水資源”

16、這一術(shù)語至今,越來越多的理論與實(shí)踐證明,土壤水不但可以計(jì)其量,而且可以辨其質(zhì)。因此,從科學(xué)的角度上,應(yīng)該把它當(dāng)作一種自然資源來加以研究。然而,目前大多數(shù)評(píng)價(jià)土壤水資源的方法是時(shí)段水量平衡法。此種方法誤差較大,特別是對(duì)于灌區(qū)土壤水。而數(shù)值方法具有常規(guī)方法不能比擬的優(yōu)點(diǎn),它可以正確地反映出農(nóng)田土壤水的動(dòng)態(tài)特征?；谝陨?筆者利用欒城地區(qū)的相關(guān)資料,采用數(shù)值模擬的方法對(duì)進(jìn)行了土壤水分運(yùn)移的計(jì)算與評(píng)價(jià)。2 模型的建立2.1 HYDRUS-1D模型的水分運(yùn)動(dòng)的基本方程選用一維飽和非飽和模型來模擬該土壤剖面上的水分變化規(guī)律。其數(shù)學(xué)模型為：式中：土壤含水率(體積%)；K土壤水導(dǎo)水率(cm/d)，在飽和土壤中

17、，其值與滲透系數(shù)相同；S(z,t)t時(shí)刻z深度處耗水速率，取該處根系吸水率(cm3/d)；h壓力水頭值(cm)；水流方向與垂直方向上的夾角。2.1土壤水分特征方程式中：s飽和含水率；r殘余含水率；，m，n經(jīng)驗(yàn)參數(shù)；Se有效含水率；Ks飽和滲透系數(shù)；Kr相對(duì)滲透系數(shù)；h土壤基質(zhì)勢(shì)。3 研究區(qū)概況欒城縣系石家莊市東南郊區(qū)縣。位 于 北 緯 37°47_36_38°01´， 東經(jīng)114°29´114°47´40_之間，地處燕山東麓的洪積沖積平原，總面積360km2。該縣西靠太行山, 東部是遼闊的華北大平原, 是淖沱河、槐沙河沖、洪

18、積扇的一部分, 地形從西北緩緩向東南傾斜, 地面高程在4566之間, 地面坡度在1/5001/1000之間。該縣屬于暖溫帶半濕潤半干旱季風(fēng)氣候區(qū)，受大陸性季風(fēng)氣候的影響具有干寒同季、雨熱同期、干濕季節(jié)明顯、四季比較分明，春季干燥多風(fēng), 夏季炎熱多雨, 秋季溫和涼爽, 冬季寒冷寡照、少雨雪的氣候特征。年平均氣溫12.2，無霜期平均約200 天，年平均降水為536.7mm，降水主要集中在69 月份，為480.9mm。栗城縣多年平均氣溫12.2, 10的積溫4251, 無霜期191天, 年日照時(shí)數(shù)2544小時(shí), 相對(duì)濕度65%, 多年平均年蒸發(fā)量1644.5mm, 年平均風(fēng)速2.6m/s, 多年平均

19、降雨量一年483.5mm , (19492000)?？h轄5鎮(zhèn)3鄉(xiāng)194村，人口35萬，面積379平方公里；耕地47萬畝，均為沃土良田，適宜耕作；水質(zhì)優(yōu)良，節(jié)水農(nóng)業(yè)正在普及。欒城縣地勢(shì)自西北向東南傾斜，地面高程 4565m, 高差約 20m。從區(qū)域性地貌特征來看全縣可劃分為 3 個(gè)不同的地貌單元，即滹沱河洪積沖積扇中部的西南邊緣部分，槐沙河洪積沖積扇的北部以及兩河之間的扇間碟形交接洼地。由于歷史上歷次洪水泛濫，河流改道而形成一系列古河道、崗地和洼地，再加上許多人工坑塘，以至造成微地形變化較大，大平小不平。同時(shí)洪積沖積物的黃土性母質(zhì)的分異，以及長期以來人為耕作活動(dòng)的影響，土壤的區(qū)域性差異明顯，土壤

20、基層單元類型比較復(fù)雜，并各具相應(yīng)的農(nóng)業(yè)生產(chǎn)特性。 根據(jù)實(shí)測(cè)的欒城地區(qū)的水位動(dòng)態(tài)變化數(shù)據(jù)及土壤初始含水量數(shù)據(jù)繪制下列曲線：圖1 水位動(dòng)態(tài)變化曲線 圖2 降雨量及蒸發(fā)量變化曲線表1 土壤巖性及含水量調(diào)查如下表：土壤深度(cm)土壤質(zhì)地 rsnKsl0-20壤土0.0598 0.4243 0.0134 1.5740 350.520-40壤土0.0694 0.4300 0.0140 1.3820 420.540-80壤土0.0600 0.4300 0.0315 1.1790 200.580-120輕粘壤土0.2294 0.3834 0.0329 1.1500 400.5120-260粘壤土0.0890

21、 0.4111 0.0151 1.1470 150.5260-340輕粘壤土0.1200 0.3836 0.0086 1.1530 100.54 概念模型本研究模擬過程中對(duì)HYDRUS一1D模型進(jìn)行了一些簡(jiǎn)化，通過數(shù)值模擬考察欒城地區(qū)包氣帶土壤水分的運(yùn)移情況。對(duì)欒城地區(qū)土壤水試驗(yàn)場(chǎng)包氣帶裸露區(qū)非均質(zhì)水分運(yùn)移規(guī)律進(jìn)行模擬。剖面巖性結(jié)構(gòu)從上到下共6層；包括：壤土、輕粘壤土及粘壤土3種巖性類型；水位埋深最大不超過33m；模型運(yùn)行期為2003年1月1號(hào)至12月31號(hào)，共365天。上邊界為變化的表層的大氣邊界，下邊界為變化的自由排水邊界；初始條件以1月2號(hào)觀測(cè)的土壤初始含水量為準(zhǔn)輸入模型進(jìn)行計(jì)算；輸出每

22、月最后一天的結(jié)果信息。本研究模擬過程中對(duì)HYDRUS1D模型進(jìn)行了一些簡(jiǎn)化，通過數(shù)值模擬考察欒城地區(qū)包氣帶土壤水分的運(yùn)移情況。對(duì)欒城地區(qū)土壤水試驗(yàn)場(chǎng)包氣帶裸露區(qū)非均質(zhì)水分運(yùn)移規(guī)律進(jìn)行模擬。剖面巖性結(jié)構(gòu)從上到下共6層；包括：壤土、輕粘壤土以及粘壤土亞粘土3種巖性類型；地下水位埋深在3133 m波動(dòng)；模型運(yùn)行期為2003年1月1號(hào)至12月31號(hào)，共365天。由于地下水位埋深比較大，故可以將研究區(qū)邊界概化為：上邊界為表層的大氣邊界，下邊界為自由排水邊界。圖3 剖面巖性分布圖、頻譜5 數(shù)學(xué)模型及其求解HYDRUS1D模型由美國鹽土實(shí)驗(yàn)室開發(fā)，用于模擬非飽和介質(zhì)中一維水分、熱和溶質(zhì)運(yùn)移過程，是模擬多孔介

23、質(zhì)中飽和非飽和水分流動(dòng)和溶質(zhì)運(yùn)移的一維模型，該模型可以考慮多種源匯項(xiàng)影響及邊界條件，其水分和溶質(zhì)控制方程應(yīng)用迦遼金線性有限元法求解。本文只介紹一維水分運(yùn)移過程。模型中水分運(yùn)動(dòng)基本方程: 式中:為土壤體積含水率;h為水頭;K為水力傳導(dǎo)系數(shù);S為根系吸水項(xiàng); e。r為有效土壤含水率; e，為殘余土壤含水率; s為飽和土壤含水率;Ks為滲透系數(shù);n，m和a均為經(jīng)驗(yàn)參數(shù)。6 數(shù)值模擬模型6.1時(shí)間離散和空間離散模擬時(shí)間離散，初始時(shí)間為0d，最終時(shí)間為365d，初始步長為1d，最小時(shí)間步長為0.01d，最大時(shí)間步長為5d?？臻g離散，模型深度取地下340cm，土壤共分為6種類型，按1cm 等間隔剖分成34

24、0個(gè)單元，并每層設(shè)置一個(gè)觀測(cè)點(diǎn)。6.2初始條件(t=0)和邊界條件(t>0)(1) 初始條件：根據(jù)現(xiàn)有資料，將實(shí)驗(yàn)?zāi)M中的初始條件選取為初始剖面體積含水率，即 =。下表為不同深度的含水量：表2 不同土壤深度對(duì)應(yīng)的初始含水量深度/cm含水量深度/cm含水量100.2697 1200.1997 200.2932 1400.2042 300.2035 1600.1673 400.1811 1800.1614 500.1751 2000.3468 600.1843 2200.3516 700.2160 2400.3696 800.2224 2600.4091 900.1821 3000.3576

25、 1000.1799 3400.3971 圖4 初始含水量深度曲線 (2) 邊界條件：根據(jù)現(xiàn)有資料中的地下水水位埋深，可知地下水水位一般在3133 m波動(dòng)，故可將邊界條件概化為：上邊界為表層的大氣邊界，下邊界為自由排水邊界。 上邊界（地表）條件:上邊界設(shè)為無徑流蒸發(fā)入滲通量邊界, 即降水和灌溉后, 水分很快下滲到下層土壤并存在蒸發(fā)，數(shù)學(xué)模型為：或Etp t0，z=0或 或Etp t0，z=0式中：為降雨或灌水強(qiáng)度，Etp潛在蒸發(fā)量。根據(jù)欒城地區(qū)植被覆蓋的情況，Etp=0.6*，為潛在騰發(fā)量3。作物潛在藤發(fā)量的計(jì)算方法有很多, 常用的有空氣動(dòng)力學(xué)、能量平衡法、彭曼公式法和經(jīng)驗(yàn)公式法。一般的模型中

26、采用彭曼公式進(jìn)行計(jì)算其公式如下： 溫度飽和水汽壓關(guān)系曲線在T處的切線斜率，kPa-1；T平均氣溫，Rn凈輻射，MJ/m2d；ed實(shí)際水汽壓，kpa；Ea飽和水汽壓，kpa；G土壤熱通量，MJ/m2d；濕度計(jì)常數(shù)，mbar-1；潛熱，MJkg-1。而在此模型中根據(jù)實(shí)測(cè)資料可得該區(qū)的蒸發(fā)量。 下邊界條件：地下水埋深邊界：地下水埋深較大，全年基本在31-33m間波動(dòng)，而研究深度為340cm，距離地下水水位很大，所以將下邊界設(shè)為自由排水。6.3土壤水分運(yùn)動(dòng)參數(shù)的確定根據(jù)大量土壤水分特征曲線變化特征，人們先后提出了不同類型和特征的描述土壤水分特征曲線的經(jīng)驗(yàn)公式。常用的描述土壤水分特征曲線的經(jīng)驗(yàn)公式主要有

27、:Gardner（1970)模型、Brooks一Corey(1964)模型、Campber(1974)模型、 van Genuehten(1980)模型、Gardner-Russ(1988)模型、Hutson一cass(1987)分段函數(shù)模型等。本文水分特征曲線采用van Genuchten方程:模型中采用了相對(duì)飽和度來描述土壤含水量狀況，即考慮土壤滯留含水量的影響。van Genuchten模型具有連續(xù)性，適用土壤質(zhì)地范圍比較寬，應(yīng)用土壤含水量范圍較廣，同時(shí)它可以使飽和土壤吸力為0，符合吸濕過程中土壤吸力變化特點(diǎn)。非飽和導(dǎo)水率k(h)由Jackson公式得到:式中:s為土壤飽和含水量;ks為

28、土壤飽和導(dǎo)水率。7 模型識(shí)別參數(shù)的初值可以通過試驗(yàn)室測(cè)定結(jié)果來給出，但最終參數(shù)的確定需要通過對(duì)模型的識(shí)別與驗(yàn)證才能確定。同時(shí)對(duì)參數(shù)的敏感性分析可進(jìn)一步得出土壤數(shù)值模型的準(zhǔn)確性與否。該模型土壤巖性以及土壤土壤水分運(yùn)移參數(shù)分析結(jié)果如下表所示： 表3 土壤巖性及水分運(yùn)移參數(shù)土壤深度(cm)土壤質(zhì)地 rsnKsl0-20壤土0.0598 0.4243 0.0134 1.5740 350.520-40壤土0.0694 0.4300 0.0140 1.3820 420.540-80壤土0.0600 0.4300 0.0315 1.1790 200.580-120輕粘壤土0.2294 0.3834 0.03

29、29 1.1500 400.5120-260粘壤土0.0890 0.4111 0.0151 1.1470 150.5260-340輕粘壤土0.1200 0.3836 0.0086 1.1530 100.58 模擬結(jié)果分析由于所研究的土壤層分為六層，因此在每層設(shè)置了一個(gè)觀測(cè)點(diǎn)，分別位于：10、35、70、100、220、310cm處，以下是模擬后各點(diǎn)的觀測(cè)情況：8.1 不同埋深點(diǎn)包氣帶剖面含水量隨時(shí)間的變化曲線，如下圖所示：h=35cmh=10cmh=310cmh=220cmh=100cmh=70cm8.2 不同時(shí)刻剖面上含水量隨深度變化曲線(t0<t1<t2<t3)8.3 上下邊界水分通量分析1、上邊界土壤水分流量2、 下邊界土壤水分流量3、土壤累計(jì)入滲量和累計(jì)潛在蒸發(fā)量8.4 水均衡信息（水流的輸入及輸出）表3 選取均衡時(shí)段內(nèi)各均衡項(xiàng)數(shù)值Time（day）In-flowL/TTop