HydrusD簡(jiǎn)明使用手冊(cè)

1、用HYDRUS-1 D模擬剖面變飽和度地下水流（簡(jiǎn)明手冊(cè)）王旭升中國(guó)地質(zhì)大學(xué)（北京）目錄1如何獲取HYDRUS-1D 錯(cuò)誤!未指定書簽。2 版權(quán)聲明.錯(cuò)誤！未指定書簽。3 參考資料.錯(cuò)誤！未指定書簽。4 . HYDRUS-1D的WINDOWS界面錯(cuò)誤!未指定書簽。5 設(shè)計(jì)模型錯(cuò)誤!未指定書簽。6 使用 HYDRUS-1D 創(chuàng)建模型錯(cuò)誤!未指定書簽。7 輸入模型控制信息錯(cuò)誤!未指定書簽。8 水流模型一一迭代計(jì)算參數(shù)錯(cuò)誤!未指定書簽。9 水流模型一一土壤水力特性模型錯(cuò)誤!未指定書簽。10 .水流模型一一土壤水分特征曲線錯(cuò)誤!未指定書簽。11 水流模型一一邊界條件錯(cuò)誤!未指定書簽。12 .水流模型一

2、一定水頭或通量邊界設(shè)置錯(cuò)誤!未指定書簽。13 . 根系吸水吸水模型 錯(cuò)誤!未指定書簽。14 .根系吸水一水分脅迫參數(shù) 錯(cuò)誤!未指定書簽。15 .輸入可變邊界條件的信息錯(cuò)誤!未指定書簽。16 .編輯土壤剖面一一使用圖形界面 錯(cuò)誤!未指定書簽。17 編輯土壤剖面一一使用表格錯(cuò)誤!未指定書簽。18 運(yùn)行模型錯(cuò)誤!未指定書簽。19 .察看結(jié)果錯(cuò)誤!未指定書簽。20 輸出結(jié)果錯(cuò)誤!未指定書簽。HYDRUS-1 D是一個(gè)共享專業(yè)軟件，用于模擬一維變飽和度地下水流、根系吸 水、溶質(zhì)運(yùn)移和熱運(yùn)移。本手冊(cè)只介紹應(yīng)用HYDRUS1D模擬垂向剖面水流和根系吸 水的操作方法。1 如何獲取HYDRUS-1DHYDRUS

3、-1 D由位于歐盟捷克的PC-Progress工程軟件開發(fā)公司發(fā)行，用戶可 以 登錄該公司首頁(yè)：。為了下載HYDRUS-1D，應(yīng)先注冊(cè)成為用戶，然后下載Hydrus- 1D的安裝文件：H1D_4_14.exe。這個(gè)文件對(duì)應(yīng)目前HYDRUS-1D的最高 版本。2 .版權(quán)聲明HYDRUS-1D的作者為：(1) J.Sim un ek,Departme ntofE nviro nmen talScie nces,U niversityofCaliforniaRiversid e,Riverside,California,USA.(2) M.Sej na,PCProgress,Prague, Czec

4、h Republic.(3)感謝他們提供了一個(gè)如此精美而又免費(fèi)使用的專業(yè)軟件，幫助我們從事有關(guān)的 科學(xué)和教育工作。當(dāng)你運(yùn)行H1D_4_14.exe解壓文件后，會(huì)在您的電腦中產(chǎn)生一個(gè)安裝目錄，其中 包含Setup.exe可執(zhí)行文件。運(yùn)行這個(gè)文件即可安裝HYDRUS-1 D軟件。當(dāng)您安裝HYDRUS-1D時(shí)，象安裝其它軟件一樣，會(huì)出現(xiàn)一個(gè)許可協(xié)議，從中可 知本共享軟件也受到美國(guó)法規(guī)的保護(hù)。3參考資料HYDRUS-1D安裝之后、在軟件運(yùn)行目錄下有HYDRS-1DManual.pdf文件。從這 個(gè) 文件您可以了解到HYDRUS-1 D的一些技術(shù)細(xì)節(jié)，如水流、溶質(zhì)運(yùn)移、熱流的方程、 一些處理專門問題的模

5、型、輸入輸出文件等等。有一個(gè)Examples目錄，包含大量的模擬算例可供參考。用戶還可以參考以下文獻(xiàn)：7 * *這些文獻(xiàn)都可以從下載。4 . HYDRUS_1D 的 WINDOWS 界面運(yùn)行HYDRUS-1D，可以看到一個(gè) Windows的界面如下:圖1所有的前后處理在界面中一目了然，左邊是前處理工具，右邊是后處理工具。其中前 處理的各項(xiàng)功能如下圖所示。¥DRIJ： i U) ROO1UPTKFile View Pre t-rcH：e®sijiE CalculaticiLid aflai!綱 sLauij i 材 ROOTUI TKPre-proeesqin 八Maio P

6、togcseesGeDEBtry Informal ion八；Z Time IrsffiTMti anFriat InfbnBntiQnH :怦 t Flow It erat ion CritAriA模擬內(nèi)容選項(xiàng)幾何形狀參數(shù)及剖面方式時(shí)間信息輸出方水流一一迭代求解控制參數(shù)水流 一一土壤水分特征模型水流一一 土壤水分特征曲線參數(shù)水流一一 邊界條件溶質(zhì)運(yùn)移一般信息溶質(zhì)運(yùn)移運(yùn)移參數(shù)溶質(zhì)運(yùn)移一一化學(xué)反應(yīng)參數(shù)溶質(zhì)運(yùn)移邊界條件根系吸水一一模型根系吸水一一水分脅迫函數(shù)可隨 時(shí)間變化的邊界條件土壤剖面一 一圖形界面土壤剖面一一數(shù)據(jù)列隔水底. I Wattr Plow - Ssi I Adraulic Prqp

7、Efty Ifodf m Water Plow -Sail Efc-draulic Pnar 氏 b與 i I W-sti-r* Flow * Boundary ConrditionEImeI -1 J* TrujiLpurt - GArArail lnAsrnfttim23 Jalutc Ireflisrpart - Traiisporr f sruTctcrw，、53 JolwCc Tronepart - F a tLoD ParaEctsrs IBM Solute Tranfig>art 一 BDmdary Conditions $ WiTET - NodelsFLi t (Chu

8、t 久 Water Siresti Re JUG 1 i caiF(r) V-3r: ahlr Boundary Condstlonp vJ Soil RraHle -Graphical Editw 刪、Soil Prgfi 1 9 SiMn ryFvr Help, pirua F圖25 .設(shè)計(jì)模型在使用HYDRUS-1D之前，您需要對(duì)飽和-非飽和水流模擬的基本原理有所了解，并設(shè) 計(jì)出自己想做的模型，準(zhǔn)備好數(shù)據(jù)。一個(gè)剖面水流模型通常包含以下幾個(gè)要素：(1) 土壤剖面從地面算起的深度，準(zhǔn)備模擬那個(gè)時(shí)間段的水分變化(2) 土壤分幾層，每層土壤的滲透性參數(shù)和水分特征曲線是怎樣的(3)根系是怎么分布的

9、。(4)是否已經(jīng)確定地面降雨入滲、蒸發(fā)蒸騰的信息，特別是它們隨時(shí)間的變化(5)是否已經(jīng)確定剖面底部的狀態(tài)屬于哪種類型的邊界條件。F面是一個(gè)參考模型的設(shè)計(jì)圖6,使用HYDRUS-1D倉(cāng)IJ建模型打開HYDRUS-1D軟件，選擇“ File/new”菜單，新建一個(gè)模型。在name 一欄中輸 入本模型的名稱“test”，更改模型存放的目錄。圖4需要注意的是，HYDRUS-1D模型本身在計(jì)算機(jī)中就表現(xiàn)為一系列的輸入輸出文 件，它們存放在與模型名稱一致的目錄中。本例中，軟件會(huì)自動(dòng)創(chuàng)建一個(gè)名稱 為“ test”的目錄，而"CATOOLSHYDR1 D'Projects”中除了 test目

10、錄之外，還有一個(gè) test.hid文件。這是一個(gè)模型項(xiàng)目（project）文件，告訴軟件下次到哪里去尋找模型。 模型創(chuàng)建之后，會(huì)顯示前處理和后處理窗口（圖5）。由于是新模型，還沒有任何模擬結(jié)果，所以后處理窗口是空白的。圖57輸入模型控制信息首先,在前處理窗口雙擊Main Processes在彈出的對(duì)話框中輸入模型的描述:atestmodel然后在Simulate一欄中選中RootWaterllptake,表示想處理根系吸水 問 題。電擊0K之后，前處理窗口將增加處理根系吸水的工具條。圖6下一步，是輸入模型的幾何信息和土層劃分信息。在前處理窗口雙擊 Geometryl nformatio n ,

11、在彈出的對(duì)話框中輸入如圖7所示的數(shù)據(jù)。模型有2個(gè)土層+長(zhǎng)度單位是cm土壤剖面的深+度是300 cm圖7接下來(lái)輸入時(shí)間信息，在前處理窗口雙擊Timeinformation，會(huì)彈出一個(gè)對(duì)話框（圖 8）。時(shí)間單位+上邊界是隨時(shí)間 變化的，每天一組 數(shù)據(jù)哄30組數(shù)據(jù)模擬30d內(nèi)的變化+時(shí)間步長(zhǎng)信息能螂雌24小時(shí)內(nèi)的圖8這個(gè)對(duì)話框中提供了一些靈活的選項(xiàng)來(lái)處理上邊界條件的變化，下面簡(jiǎn)要加以說(shuō)明： (1)蒸騰量的每日周期變化HYDRUS-1D可以使用一個(gè)經(jīng)驗(yàn)公式來(lái)處理每天24小時(shí)潛在蒸騰量的變化，設(shè)某天的 潛在蒸騰量為Tp (例如用Pemman公式獲取的，cm/d)，則其中Tp是瞬時(shí)潛在蒸騰量，t為時(shí)間。模

12、型假設(shè)早上6點(diǎn)之前以及晚上18點(diǎn)-24點(diǎn)的 蒸騰量總和只占全天蒸騰量的1%。注意本例中蒸騰量的單位是cm/do(2)降水量的周期變化如果在你的模型中降水量是周期性變化的，HYDRUS-1D也可以用一個(gè)公式來(lái) 處 理其中P是周期t內(nèi)的平均降雨量。(3)使用氣象數(shù)據(jù)也可以在HYDRUS-1D中輸入氣象數(shù)據(jù)，它將自動(dòng)利用這些數(shù)據(jù)計(jì)算潛在蒸 散量ETp?？梢赃x擇FAO組織推薦的Penman-Monteith公式，也可以選擇Hargreaves 公式。這些公式需要輻射、氣溫、濕度之類的氣象數(shù)據(jù)?？梢源_定需要輸出哪 些時(shí)間點(diǎn)的結(jié)果模型的另一個(gè)控制信息是對(duì)模擬結(jié)果的輸出如何進(jìn)行設(shè)置。在前處理窗口雙擊 Pri

13、ntin formation工具條，彈出一個(gè)對(duì)話框。本例中確定輸出30組模擬結(jié)果，每天輸 出1組。FTIffiT on5Inftriritiica£網(wǎng)嚇ElQutpatP?3 a.1 FI UJ BT (xni L«i-Ia£ rMP： fcr Bhtcr 瑞L：的 S« £.ccP Hit 犁 id* FiTT輸出30組結(jié)*果，每天1組圖98水流模型一一迭代計(jì)算參數(shù)HYDRUS-1D是采用迭代法來(lái)處理非線性Richards方程的。在前處理窗口雙擊 WaterFlow-lterationCHteria工具條，彈出一個(gè)設(shè)置迭代參數(shù)的對(duì)話框(圖10

14、)。迭代控 制參數(shù)的設(shè)置具有高度的專業(yè)技術(shù)性，除非特別了解，一般可以使用默認(rèn)值。如果模擬 結(jié)果出現(xiàn)不收斂的情況，需要對(duì)最大迭代次數(shù)、迭代精度等參數(shù)進(jìn)行調(diào)整，但是在缺乏 經(jīng)驗(yàn)的情況下很難操作。on CriLc r i n最多迭代次數(shù)含水量迭代精度壓力水頭迭代精度增大步長(zhǎng)迭代次數(shù)信號(hào)縮小步長(zhǎng)的迭代次數(shù)信號(hào)步長(zhǎng)增大比例步長(zhǎng)縮小比例最小吸力間距最大吸力間距ioVuinuFi l»b «ar 41|0. XIIB.«Let CcdjLIEiil Tvler-uufieErassiizrc： jtao.4 T<Jor ui=<ii!7liuI3Lover lto!

15、.電 Arc.1 l cm=Vjjptf-r Ipt. nal Tleriiti,Xi i不 <>1 3Lwa Tiint H'jlupljTulGFU T Ippa* TmAi SiapInt arAial Eni &irj *ct luhtamk1 e-flO#!lAivti Lint t *f ilie 7ensi ki IntervilOROPppei- ZJn 11 *E tht I eas： ot Interred.Jacx&r用于生成水分特征曲線的散點(diǎn)值圖10HYDRUS-1D采用自動(dòng)控制時(shí)間步長(zhǎng)的方法來(lái)處理迭代的收斂性。對(duì)于每個(gè)時(shí) 步，如果迭

16、代次數(shù)太多，就縮小時(shí)間步長(zhǎng)；如果沒經(jīng)過幾次迭代就達(dá)到收斂精度，則 適當(dāng)增大時(shí)間步長(zhǎng)。9 .水流模型一一土壤水力特性模型水分特征曲線是非飽和土壤的重要物理性質(zhì)，HYDRUS.1D提供了幾種方法來(lái)處理與之有關(guān)的參數(shù)。在前處理窗口雙擊WaterFlow-SoilHydraulicProperties IM條， 彈出一個(gè)設(shè)置水力特性模型的對(duì)話框（圖11 ）。單孔介質(zhì)模型=雙重介質(zhì)模型雙重孔隙度雙重滲透性吸濕和疏干滯后過程模型Wydrml'i e S 0,A1 # F-sraus. by Hc<dal tr M sK Yil則誨吃e U ALj Fa vbt： Gi-T.-z.rLh -1

17、T1i Br4iaki.-Caf ay 11)OmJ SZ .it y.lfQ.tal. J hJI ta ybkla I.V lyndulL 廣m 迎.即 卬& J 1" VSife fWldbttA , -31$-D'H e Jfaii1： ? (in h bil.t-i.rFdSTd 奇 審物二燈 c. mi tf sii . Oiiil-E JiffltLrjr 區(qū)晦 比廠中“ BH !也11自前羽b««>lw ht« rsaJy f«r aK|i«ri«M<>JrlXf "

18、;Igjf DXjj.-!； irith 4 AHJxA (G-ijrlcii. ood v«ii <»DUjelii4M 1 ±33)U3 叫 1 Bill V/1.Ifcra-tcre-E 5 斤 In Lvklvfl-BS.Ii i.町臼菲：>£ ill k- diL4 id- ilWArx lwbl3 x a. r L*a.i i.-« 3i =£”/ un £ ncainEt lvs iy 電 I tL，n Lwal Lvh. 叱ilhrpk tiS ： |1 口卜 Lefeftialr J)主要在溶質(zhì)

19、運(yùn)移問題中使用圖11在一般情況下，選擇單孔介質(zhì)模型，并選擇用vanGenuchten-Mualem公式處理土 壤的水力特性就可以了。如果還要模擬溶質(zhì)運(yùn)移，可能需要考慮雙重介質(zhì)模型。雙重介質(zhì)在同一個(gè)點(diǎn)有兩 個(gè)孔隙度或兩個(gè)滲透率，相當(dāng)于兩種介質(zhì)的混雜。雙重介質(zhì)模型能夠模擬這兩種“介 質(zhì)”之間的水分和鹽分交換。10 .水流模型一一土壤水分特征曲線在前處理窗口雙擊WaterFlow-SoilHydraulicParameters工具條、彈出一個(gè)設(shè)置 水 分特征曲線參數(shù)的對(duì)話框（圖12）。本例中選擇vanGenuchten-Mualem公式處理水分特征曲線，其中,nJ均為控制因 子。HYDRUS-1D軟

20、件中提供了一組土壤經(jīng)驗(yàn)參數(shù)庫(kù)，可供用戶參考。本例中兩層土 壤的參數(shù)直接從數(shù)據(jù)庫(kù)中調(diào)出：第1層對(duì)應(yīng)Sandyloam,第二層對(duì)應(yīng)sand。細(xì)砂壤土中砂Water Flow Parameters圖12在輸入?yún)?shù)時(shí)，請(qǐng)注意參數(shù)的單位。11 .水流模型一一邊界條件在前處理窗口雙擊WaterFlowBoundaryConditions工具條，彈出一個(gè)設(shè)置邊界 條件的對(duì)話框(圖13)。Water Flow Boundary Condi地面邊界類型定壓力水頭定水分通量大氣邊界，可積水大氣邊界，產(chǎn)流變水頭變水頭和通量下端邊界類型定壓力水頭定水分通量變地下水位變流量自由下滲排水滲出面水平排水7pp«r

21、 Bouiidlury Condi ti onS' " Ccihl IulL廠 CtmtltFit Fl-iscb Atrnicspheric BZ Witt Surfice Laver Atinitspheri c BT witB Swfni.ce Run. IffC VtriiKle Prssswe_ov*r RonridliHry C on Ji 七 口口MIQPrsEsurA Kead 金 Fl UE廠 JuriaKle Pr*ss,ure K«&dCFL電Tai li. al. CcfficLi £ i vxl* Id Pressuric

22、 m血如In.初始條件使用壓力水頭使用含水量L Fre« Jrai nafieDeep AraiaaAe（斗町 4爐 Fa.ca ; hHflrizsntftlc VC大氣邊界I* Input PET ud LA基于ETp和LAI消光系數(shù)最大積水厚度圖13上邊界條件有6種類型，下邊界條件有8種類型。邊界類型的確定需要考慮實(shí)際 條件，在本算例中，上邊界選擇大氣邊界條件，在降雨量很大時(shí)地表可以產(chǎn)生積水。 植被蒸騰量和土壤蒸發(fā)量分開處理，HYDRUS-1D推薦使用一個(gè)經(jīng)驗(yàn)公式 來(lái)把潛在蒸 散量分割為蒸騰潛力和土壤蒸發(fā)潛力：其中ETp為潛在蒸散量(可以使用Penman- MontHth公式處

23、理氣象數(shù)據(jù)得到，cm/d)，Tp為潛在蒸騰量(cm/d)，Ep為土壤潛在蒸 發(fā)量(cm/d)，LAI是葉面積指數(shù)，k為消光系數(shù)，取決于太陽(yáng)角度、植被類型及葉片空 間分布特征。SCF是一個(gè)中間參數(shù)，即土壤覆蓋度(Soilcoverfraction)。在闊葉植被發(fā)育的情況下，消光系數(shù)的經(jīng)驗(yàn)值為 k=0.5-0.75。 12.水流模型一一定水頭或通量邊界設(shè)置如果邊界條件中包含定水頭或定通量的邊界，則在前處理窗口雙擊WaterFlow-Co nsta ntBC工具條，彈出一個(gè)設(shè)置邊界數(shù)據(jù)的對(duì)話框。本算例模型中，下邊界為定流量 邊界，實(shí)際上就是隔水邊界，因此直接輸入0即13.根系吸水一一吸水模型在前處理

24、窗口雙擊RootWaterllptake-Models工具條，彈出一個(gè)處理根系吸水 模 型的對(duì)話框(圖14)。補(bǔ)償吸水域值 (臨界濕潤(rùn)度)水分脅迫模型鹽分脅迫模型根系吸鹽模型=0圖14HYDRUS-1D使用水分脅迫和鹽分脅迫模型處理根系的吸水。對(duì)于水分脅迫模 型，計(jì)算公式為其中Tp是潛在蒸騰量(cm/d), Ta是實(shí)際蒸騰量(cm/d)，S(x)是吸水強(qiáng)度函數(shù) (cm/(cm.d),注意x坐標(biāo)實(shí)際表示深度)，(h)是水分脅迫函數(shù)，h為土壤壓力水頭(cm)， b(x)是根系吸水分配(密度)函數(shù)，lr為根系層的深度。水分脅迫函數(shù)有2種經(jīng)驗(yàn) 表示 方法，即Feddes模型和S-Shape模型。Fed

25、des模型是一個(gè)梯形函數(shù)，只需要 知道h 值。而S-Shape模型把水分脅迫和葉片氣孔的壓力水頭聯(lián)系起來(lái)，需要知道氣孔壓力 水頭的數(shù)值ho土壤的濕潤(rùn)度可以表示為Ta/Tp (h,x)b(x)dx(2.22) LR但是如果直接用這種方法來(lái)計(jì)算實(shí)際蒸騰量有一定的問題。植被其實(shí)可以調(diào)節(jié)不同深 度的水分脅迫響應(yīng)特征；某個(gè)深度土壤干燥吸不上水，植被可以加大在比較 濕潤(rùn)的土 層的吸水量，以補(bǔ)償不足。這種現(xiàn)象稱為補(bǔ)償吸水。為了模擬根系補(bǔ)償吸水，HYDRUS-1D提供了一種簡(jiǎn)化的模型，即如果濕潤(rùn)度高 于某個(gè)臨界值(>c)，植被根系可以通過補(bǔ)償機(jī)制充分吸水達(dá)到潛在蒸騰量。如果濕潤(rùn)度 低于這個(gè)臨界值，補(bǔ)償機(jī)

26、制受到抑制，發(fā)生整體的水分脅迫，根系吸水總量將低于潛 在蒸騰量，并正比于濕潤(rùn)度。如果不考慮這種補(bǔ)償吸水機(jī)制，可以令c=1 o 14.根系吸水 水分脅迫參數(shù)在前處理窗口雙擊 RootWaterllptake-WaterStressReductior，彈出一個(gè)處理水分脅迫參數(shù)的對(duì)話框(圖15)。本算例中直接從數(shù)據(jù)庫(kù)中調(diào)入Wheat的經(jīng) 驗(yàn)值。P3P2LP2HPoptP0壓力水頭(pressurehead,cm)一一兩個(gè)臨界吸水強(qiáng)度-數(shù)據(jù)庫(kù)脅迫 函數(shù)=1圖15根系水分脅迫Feddes模型參數(shù)15 .輸入可變邊界條件的信息在前處理窗口雙擊VariableBoundaryConditions，彈出一個(gè)處

27、理時(shí)間序列數(shù)據(jù)的對(duì)話框(圖16) o本算例中，在步驟(11)中已經(jīng)把地面處理大氣邊界，同時(shí)又選擇使用消光系數(shù) 法 劃分植被蒸騰和土面蒸發(fā)，因此需要輸入每天的降水、潛在蒸散量、葉面積指數(shù)等數(shù) 據(jù)。還有一個(gè)需要輸入的數(shù)據(jù)是最小壓力水頭值，即地面土壤達(dá)到最干燥狀態(tài)時(shí)的壓 力水頭。從理論上講，當(dāng)土壤十分干燥時(shí)，吸力很大，而液態(tài)孔隙水的壓強(qiáng)很小，與 空氣濕度保持平衡關(guān)系，因此有其中hA為最小壓力水頭，H為空氣絕對(duì)濕度，RT/Mg 為空氣的摩爾氣體常數(shù)?？諝鉂穸入m然可以通過氣象數(shù)據(jù)得到，但這里公式需要的是 近地面的空氣濕度。一般情況下，取飽和水汽濕度是可取的，因?yàn)? cm深度以下土 壤空氣的濕度往往都是飽

28、和的，只不過隨溫度發(fā)生變化。因此，可以根據(jù)近地面氣溫 的變化來(lái)推算地 表土壤的空氣濕度(飽和水汽濕度)，再換算成壓力水頭。HYDRUS-1D 中需要輸入的是最小壓力水頭的絕對(duì)值，缺省值為hCritA=|hA|=106 cm=104 m這個(gè)數(shù)值只會(huì)對(duì)土壤蒸發(fā)起作用。HYDRUS-1D建議：hCiitA所對(duì)應(yīng)的土壤含水量應(yīng) 該至少比殘余含水量大6005,在模擬根系吸水的情況下，hA還應(yīng)該低于圖15中的P3。 否則(hA>P3),當(dāng)根系吸水的臨界值壓力水頭(P3)和地面蒸發(fā)的最小壓力水頭(hA)滿足 時(shí)，會(huì)導(dǎo)致回流(inflow)現(xiàn)象，這是不合理的。除非存在特別干燥的情況，模型一般不 需要仔細(xì)

29、處理這些問題。圖16時(shí)間(d)降水量(cm/d)ETp(cm/d)hCritA(cm)LAI時(shí)間(d)降水量(cm/d)ETp(cm/d)hCritA(cm)LAI100.51000002-1I。00.41000002.1200.6100000200.61000000.3300.41000002.21800.71000000.3400.51000002.2-00.71000000.3500.5100000222000.71000000.3600.5100000222100.71000000.3 gondltianB f旦Id ： lc 且尋學(xué)修改節(jié)點(diǎn)數(shù)pFerr 1土壤剖面分成 31個(gè)節(jié)點(diǎn),則

30、 節(jié)點(diǎn)間距為 10 cm700.31000002.3:2200.61000000.380.70.11000002.32300.51000000.3900.41000002.3240.30.21000000.31000.51000002.3252.30.11000000.31100.61000002.2261.10.11000000.31200.41000002.2270.50.11000000.31330.11000002.22800.41000000.31400.510000022900.71000000.3150.50.2100000213000.71000000.3數(shù)據(jù)可以先在Excel

31、中準(zhǔn)備好，如表1。這些數(shù)據(jù)可以拷貝到圖16的電子表格中。這些數(shù)據(jù)顯示的降水量、蒸散潛力和葉面積指數(shù)變化特征如圖16右圖所示在第17日由于莊稼收割，葉面及指數(shù)大幅度下降。16 .編輯土壤剖面一一使用圖形界面在前處理窗口雙擊SoilProfile-GraphicalEdtor工具條，程序?qū)棾鲆粋€(gè)處理土 壤剖面的軟件（圖17）。這個(gè)圖形軟件的使用比較簡(jiǎn)單，我們需要注意的是在 Conditions菜單下面有很多子菜單，包含處理各種問題的功能菜單。首先要做的事情，是確定把土壤剖面離散化為多少個(gè)節(jié)點(diǎn)。本算例土壤模型深度為3 m,我們希望節(jié)點(diǎn)間距達(dá)到1 cm，因此需要301個(gè)節(jié)點(diǎn)。選擇菜單Conditio

32、ns/ProfileDiscretization，在下拉工具條中把 Number 修改為 301。OHydrus-lD - Profile Conditions ProfileDiscretization MaterialDistribution RootDistribution Scalingfactor Initialconditions Subregions ObservationPointsConditions 菜單剖分節(jié)點(diǎn)土壤巖性分層根系分布尺度因子初始條件子區(qū)分布觀察點(diǎn)圖17接下來(lái)，確定土層的分布，本模型有2個(gè)土層。缺省的土層編號(hào)為in dex=1，就是 圖12中的細(xì)砂壤土層。選擇

33、菜單Conditions/MaterialDistribution，在下拉工具條中 使用Editcondition，把下部土壤層設(shè)置index=2，這個(gè)土層編號(hào)為2，實(shí)際上就是圖 12中的中砂層。 Hydrun-W - Profile 口 n.xConditionsMaterialDistribution?；馗被薸jaf崟(1*幣幣1.點(diǎn)擊工具條2.在剖面底部點(diǎn)擊鼠標(biāo)左 鍵再放開，從下往上移動(dòng) 鼠標(biāo)，選定土層范圍，然 后再點(diǎn)擊鼠標(biāo)左鍵，修改a-修改節(jié)點(diǎn)數(shù)7右下角是鼠標(biāo)所在的Z坐標(biāo)高度H Hydnur-t b - Prof k le Informtlonw C X使用具條圖19圖18再就是確定根

34、系隨深度的分布，根系分布函數(shù)b(x)是一個(gè)很特殊的函數(shù)，它滿足 以下條件其中LR為根系層厚度,x是深度；z是節(jié)點(diǎn)間距，出是每個(gè)間距中的根系分布函數(shù)值， M是根系層占節(jié)點(diǎn)數(shù)。本算例中根系層的厚度為1 m,假設(shè)根系分布函數(shù)為線性，并有b(x) 2x ,x1容易證明上式滿足積分為1的條件。設(shè)置方法為，選擇菜單Conditions/RootDistribution，在在下拉工具條中使用Editcondition，劃定根系層范圍 (100 cm)，把頂部數(shù)值設(shè)置為。，底部數(shù)值設(shè)置為2,讓程序自動(dòng)進(jìn)行線性差值。ConditionsRootDistribution下一步處理初始條件，選擇菜單Conditio

35、ns/lnitialCondition。本算例中初始地下 水位高于底板10。cm,假設(shè)土壤剖面初始狀態(tài)是靜力平衡態(tài)，則模型底部的壓力水頭 為+100 cm，地面的壓力水頭為200 cm。使用EditCondition工具條，分別設(shè)置頂部和 底部的壓力水頭，并讓程序自動(dòng)插值形成初始條件。最后、選擇菜單Conditions/ObservationPoints,使用Insert工具條添加若干觀察 點(diǎn)，有必要在靠近模型底部的位置加一個(gè)觀察點(diǎn)以判斷地下水位的變化。關(guān)閉圖形程 序退回主程序。17 .編輯土壤剖面一一使用表格在前處理窗口雙擊SoilProfile-Summary工具條，程序?qū)棾鲆粋€(gè)表格（圖20）。 在這個(gè)表格中可以進(jìn)一步修改土壤剖面數(shù)據(jù)。圖2