HydrusD簡明使用手冊

1、DocumentserialnumberNL89WT-NY98YT-NC8CB-NNUUT-NUT108用HYDRUS-1D模擬剖面變飽和度地下水流（簡明手冊）王旭升中國地質(zhì)大學(xué)（北京）HYDRUS-1D是一個共享專業(yè)軟件，用于模擬一維變飽和度地下水流、根系吸水、溶質(zhì)運移和熱運移。本手冊只介紹應(yīng)用HYDRUS1D模擬垂向剖面水流和根系吸水的操作方法。1. 如何獲取HYDRUS-1DHYDRUS-1D由位于歐盟捷克的PC-Progress工程軟件開發(fā)公司發(fā)行，用戶可以登錄該公司首頁：。為了下載HYDRUS-1D,應(yīng)先注冊成為用戶，然后下載Hydrus-1D的安裝文件：。這個文件對應(yīng)目前HYDRU

2、S-1D的最高版本。2. 版權(quán)聲明HYDRUS-1D的作者為：(1) J.Simunek,DepartmentofEnvironmentaiSciences,UniversityofCaliforniaRiverside,Riverside,California,USA.(2) M.Sejna,PCProgress,Prague,CzechRepublic.(3) .vanGenuchten,DepartmentofMechanicalEngineering,FederalUniversityofRiodeJaneiro,RiodeJaneiro,Brazil.感謝他們提供了一個如此精美而乂免

3、費使用的專業(yè)軟件，幫助我們從事有關(guān)的科學(xué)和教育工作。當(dāng)你運行解壓文件后，會在您的電腦中產(chǎn)生一個安裝目錄，其中包含可執(zhí)行文件。運行這個文件即可安裝HYDRUS-1D軟件。當(dāng)您安裝HYDRUS-1D時，象安裝其它軟件一樣，會出現(xiàn)一個許可協(xié)議，從中可知本共享軟件也受到美國法規(guī)的保護。3. 參考資料HYDRUS-1D安裝之后，在軟件運行目錄下有HYDRS-1D文件。從這個文件您可以了解到HYDRUS-1D的一些技術(shù)細(xì)節(jié)，如水流、溶質(zhì)運移、熱流的方程、一些處理專門問題的模型、輸入輸出文件等等。有一個Examples目錄，包含大量的模擬算例可供參考。用戶還可以參考以下文獻：?im?nektJ,MTh.va

4、nGenuchten,and?ejna,DevelopmentandapplicationsoftheHYDRUSandSTANMODsoftwarepackages,andrelatedcodes,VadoseZoneJournal,doi:,SpecialIssueMVadoseZoneModeling",7(2),587-600,200&Jacques,D,J?im?nek,DMallants,andMTh.vanGenuchten,Modelingcoupledhydrologicalandchemicalprocesses:Long-termuraniumtrans

5、portfollowingmineralphosphorusfertilization,VadoseZoneJournal,doi:,SpecialIssue”VadoseZoneModelingM,7(2),698-711,200&?im?nektJand(Th.vanGenuchten,ModelingnonequilibriumflowandtransportwithHYDRUS,VadoseZoneJournal,doi:,SpecialIssueMVadoseZoneModelingv,7(2),782-797,2008.這些文獻都可以從下載。4. HYDRUS-1D的WIN

6、DOWS界面運行HYDRUS-ID,可以看到一個Windows的界面如下:HYDRUSIDROOTUPTKEileViw£re-processingCalcvlatiooEesuits(kitionsWindwHelpD|2陶舊I創(chuàng)皂舊|D|<|麗fROOTUPTKPre-jirocessingiXflainProcessesp/GeometryInforimtioni/TimeInforwationeZPrint二Watci:UlWater:二Watci:UlWater:Solutu曲1Solutel£SJSolute磴1Solute前處理工具Information

7、Flow-Itera-tionCriteriaFlow-SoilHydraulicPropertyflodelFlow-SoilHydraulicParanglersFlow-BoundaryConditions:Transport-GeneralInfornia-tion!Transport一TransportParanetersTransport一ReactionParanglers:Transport一BoundaryConditions璧RootWaterIntake-Kadels冬RootWaterUptake-WaterStressReactionF®VariableBo

8、undaryConditions皈SoilProfile-GraphicalEditorI切ISoilProfile-SurnraryPost-processing00«Q00TO00TO«0模擬計算后處理工具ObservationPointsProfi1eInfoTHiation脳turFlow-BoundaryFluxesandHeadsSoluteTransport-ActualandCurwjlariveBoundaryFluxesSoilMjrdraulicPropertiesRunTinsInformationflassBalancenfoixiatiemFot

9、:Help,pressFl圖1所有的前后處理在界面中一目了然，左邊是前處理工具，右邊是后處理工具。其中前處理的各項功能如下圖所示。HYDRUSIDROOTUPTKEileViw£re-processingCalculationEesultdZI區(qū)舊I皂舊Iml麗IROOTUPTKPre-processingiXflainProcessesp/GeometryInforimtioni/TimeInforwationeZPrint二Watci:UlWater:二Watci:UlWater:Solutu曲1Solutel£SJSolute磴1SoluteInformationFl

10、ow-Itera-tionCriteriaFlow-SoilHydraulicPropertyflodelFlow-SoilHydraulicParanglersFlow-BoundaryConditions:Transport-GeneralInfornia-tion!Transport一TransportParanetersTransport一ReactionParanglers:Transport一BoundaryConditions璧RootWaterIntake-Kadels冬RootWaterUptake-WaterStressReactionF®VariableBoun

11、daryConditions皈SoilProfile-GraphicalEditorIWlSoilProfile-SunnaryFot:Help,pressFl1/模擬內(nèi)容選項幾何形狀參數(shù)及剖面方式時間信息輸出方式水流迭代求解控制參數(shù)水流土壤水分特征模型水流土壤水分特征曲線參數(shù)水流邊界條件溶質(zhì)運移-般信息溶質(zhì)運移-一一運移參數(shù)溶質(zhì)運移-一一化學(xué)反應(yīng)參數(shù)溶質(zhì)運移-一一邊界條件根系吸水-一一模型根系吸水-一一水分脅迫函數(shù)可隨時間變化的邊界條件土壤剖而一一圖形界面土壤剖而一一數(shù)據(jù)列表圖25. 設(shè)計模型在使用HYDRUS-ID之前，您需要對飽和-非飽和水流模擬的基本原理有所了解，并設(shè)計出自己想做的模型

12、，準(zhǔn)備好數(shù)據(jù)。一個剖面水流模型通常包含以下幾個要素：(1) 土壤剖面從地面算起的深度，準(zhǔn)備模擬那個時間段的水分變化。(2) 土壤分兒層，每層土壤的滲透性參數(shù)和水分特征曲線是怎樣的。(3) 根系是怎么分布的。(4) 是否己經(jīng)確定地面降雨入滲、蒸發(fā)蒸騰的信息，特別是它們隨時間的變化。(5) 是否己經(jīng)確定剖面底部的狀態(tài)屬于哪種類型的邊界條件。下面是一個參考模型的設(shè)計圖：6.使用HYDRUS-1D創(chuàng)建模型打HYDRUS-1D軟件，選擇”File/new”菜單，新建一個模型。在name欄中輸入本模型的名稱”test”,更改模型存放的目錄。圖4需要注意的是，HYDRUS-1D模型本身在計算機中就表現(xiàn)為一系

13、列的輸入輸出文件，它們存放在與模型名稱一致的目錄中。本例中，軟件會自動創(chuàng)建一個需稱為”test”的目錄，而”CATOOLSHYDRlDProjectsw中除了test目錄之外，還有一個文件。這是一個模型項目(project)文件，告訴軟件下次到哪里去尋找模型。、模型創(chuàng)建之后，會顯示前處理和后處理窗口(圖5)。由于是新模型，還沒有任何模擬結(jié)果，所以后處理窗口是空白的。圖57. 輸入模型控制信息首先，在前處理窗口雙擊MainProcesses,在彈出的對話框中輸入模型的描述：atestmodel.然后在Simulate一欄中選中RootWaterUptake,表示想處理根系吸水問題。電擊0K之后，

14、前處理窗口將增加處理根系吸水的工具條。圖6下一步，是輸入模型的兒何信息和土層劃分信息。在前處理窗口雙擊GeometryInformation,在彈出的對話框中輸入如圖7所示的數(shù)據(jù)。模型有2個土層長度單位是cm土壤剖而的深V度是300cm圖7接下來輸入時間信息，在前處理窗口雙擊TimeInformation,會彈出一個對話框（圖8）。TimcInformationlimeVniIsfineDiscretizatioCSecondsCDiriutesCKnursr2=FinalTineIniiialTimeStep|D0C1MinaniunTimeStep»001Mx2fTiWiTime

15、Step|1XnitialTime='上邊界是隨時間變化的，每天yI7組數(shù)據(jù)，共30<Timtt-Vori.abloBcurtiryCor»8it:cns7ime-VeriallsBoundaryCo&dilionsNinbci-o£T:ncVwiobieBozdzyRec&r4s(e.g.?模擬30d內(nèi)的變化A時間步長信息組數(shù)據(jù)30*jD'ilyVQijr.igm<£&生匕?Hep巴衽iggcr.w.cd3y二>廠SinusoidalVariaxionsof?recipitatioriGeneratedb

16、yHYDXVSAorolog;icalDitaIBtaorlogicdUatQJUNwberofl1eteorclog：LcelXecoris(ez.tKaciaticTi)?er/Tion-Fi1onthcitXEquationCXargreavesFcrjrileCZnergyEalanceSoundaryCone廠laiJyVr)q-1nn>cfHotcoIlwtQHrr】wTl矽Gcnorwt"byfllR門動處理蒸24小時內(nèi)的變化圖8這個對話框中提供了一些靈活的選項來處理上邊界條件的變化，下面簡要加以說明：(1) 蒸騰量的每日周期變化HYDRUS-1D可以使用一個經(jīng)驗

17、公式來處理每天24小時潛在蒸騰量的變化，設(shè)某天的潛在蒸騰量為(例如用Pemman公式獲取的，cm/d)，則坊(0=0.247；7(/)=2.7578111'2冗1冗'、lday2；/<0.264d./>0.736d/c(0.264cl.0.736d)(2.73)其中7；(廣)是瞬時潛在蒸騰量，廣為時間。模型假設(shè)早上6點之前以及晚上18點-24點的蒸騰量總和只占全天蒸騰量的1%。注意本例中蒸騰量的單位是cm/d。(2) 降水量的周期變化如果在你的模型中降水量是周期性變化的，HYDRUS-1D也可以用一個公式來處理P(z)=P1+cos,Az丿丿(2.74)其中戸是周期

18、?t內(nèi)的平均降雨量。(3) 使用氣象數(shù)據(jù)也可以在HYDRUS-1D中輸入氣象數(shù)據(jù)，它將自動利用這些數(shù)據(jù)計算潛在蒸散量ETpo可以選擇FA0組織推薦的Penman-Monteith公式，也可以選擇Hargreaves公式。這些公式需要輻射、氣溫、濕度之類的氣象數(shù)據(jù)。模型的另一個控制信息是對模擬結(jié)果的輸出如何進行設(shè)置。在前處理窗口雙擊PrintInformationIE具條，彈出一個對話框。本例中確定輸出30組模擬結(jié)果，每天輸出1組。輸出30組結(jié)V果，每天1組VPrintTnfonnationXP7-L<iv61TnfonitiitionOK£verynlinesteps：1Co&

19、amp;ccLPrint冬滬laTigIntwrvdrErevioxjs.j"Z廠7Scr6wxiOutputHwt.一PrintFluxes(insteadofl«ip)forObservationNodes*KtIEnteratEnd9II.PrintTinasUiiribwofPrintTinas可以確定需要輸出哪些時間點的結(jié)果SelectFrintIarits.8. 水流模型一一迭代計算參數(shù)HYDRUS-1D是采用迭代法來處理非線性Richards方程的。在前處理窗口雙擊WaterFlow-IterationCriteriaE具條，彈出一個設(shè)置迭代參數(shù)的對話框（圖1

20、0）o迭代控制參數(shù)的設(shè)置具有高度的專業(yè)技術(shù)性，除非特別了解，一般可以使用默認(rèn)值。如果模擬結(jié)果出現(xiàn)不收斂的情況，需要對最大迭代次數(shù)、迭代精度等參數(shù)進行調(diào)整，但是在缺乏經(jīng)驗的情況下很難操作。IteratiomCriteria最多迭代次數(shù)含水雖迭代精度壓力水頭迭代秸度IterationCritcriaPressureHex!5a>:inuTiNanberofIterationsContentTolerance增大步長迭代次數(shù)信號縮小步長的迭代次數(shù)信號J步長增大比例I步長縮小比例1最小吸力間距l(xiāng)T最大吸力間距qlimeStepControl3OptimoLLt«rtionUpperOp

21、tiffialIterationRangeTim。St&pHuitiplication0.7VppcrTimeSlepflulliplicationFactorIntarpolticnLimitofthalnsionIntervalIOCCOVpperLimitofth«TensionInterv<J用于生成水分特征曲線的散點值圖10HYDRUS-ID采用自動控制時間步長的方法來處理迭代的收斂性。對于每個時步，如果迭代次數(shù)太多，就縮小時間步長；如果沒經(jīng)過兒次迭代就達到收斂精度，則適當(dāng)增大時間步長。9. 水流模型一土壤水力特性模型水分特征曲線是非飽和土壤的重要物理性質(zhì)，H

22、YDRUS-1D提供了兒種方法來處理與之有關(guān)的參數(shù)。在前處理窗口雙擊WaterFlow-SoilHydraulicProperties工具條，彈出一個設(shè)置水力特性模型的對話框（圖11）。Soi1HydraulicModel單孔介質(zhì)模型HoldSinaiaPoroxityModal*&van-V*Q41en廠Mth.ir-EnurV<J,nc&f-2cnJlodiEj«dvanGenuchtenCBrooks-Corwj*K&zudQcfi-xiomAl)雙重介質(zhì)模型雙重孔隙度雙重滲透性PnAl-P-:r0?tty/Ikwl-PerntdbiL?iyVod

23、ls'Duil-porsityOixxiQy,da<dvz<kotw：itari-HuaIob)Dual-poroiity(nobjle-iwobilez*atercnastransfeCChiftl-porosity(Aobilft-irnobil«zhwadrm“trftnctar)=Model:bolcvyrccmcftdadorilyEoy4xp4rio&<:ad*11pn?】t7'J：IJirfnl*Vr二(!'IJi."Dtt41-pam<ilil7(Ocrltcgd7&aGo&'x

24、chtc&.1993)CLook-upTablx主要在溶質(zhì)運移問題中使用Myiteresis吸濕和疏十滯后過程模型e3?0hystorasisJlysteresjsmretenlioacurveCH'starwKKr«tAntioacurveandconductjvtty'fctrczds5nrctantiacutvqopwpinc,BobIhird)CIrntLilly空yi班.cor"C*H:tidilyerve圖11在一般情況下，選擇單孔介質(zhì)模型，并選擇用vanGenuchten-Mua1em公式處理土壤的水力特性就可以了。如果還要模擬溶質(zhì)運移

25、，可能需要考慮雙重介質(zhì)模型。雙重介質(zhì)在同一個點有兩個孔隙度或兩個滲透率，相當(dāng)于兩種介質(zhì)的混雜。雙重介質(zhì)模型能夠模擬這兩種“介質(zhì)”之間的水分和鹽分交換。10. 水流模型土壤水分特征曲線在前處理窗口雙擊WaterFlow-SoilHydraulicParametersX具條，彈出一個設(shè)置水分特征曲線參數(shù)的對話框（圖12）o本例中選擇vanGenuchten-Mua1em公式處理水分特征曲線，務(wù)+躋0/<07/>0(2.28)(2.29)其中222,/均為控制因子。HYDRUS-1D軟件中提供了一組土壤經(jīng)驗參數(shù)庫，可供用戶參考。本例中兩層土壤的參數(shù)直接從數(shù)據(jù)庫中調(diào)出：第1層對應(yīng)Sandy

26、loam,第二層對應(yīng)sando圖12在輸入?yún)?shù)時，請注意參數(shù)的單位。11. 水流模型一一邊界條件在前處理窗口雙擊WaterFlow-BoundaryConditionsX具條，彈出一個設(shè)置邊界條件的對話框(圖13)o圖13上邊界條件有6種類型，下邊界條件有8種類型。邊界類型的確定需要考慮實際條件，在本算例中，上邊界選擇大氣邊界條件，在降雨量很大時地表可以產(chǎn)生積水。植被蒸騰量和土壤蒸發(fā)量分開處理，HYDRUS-1D推薦使用一個經(jīng)驗公式來把潛在蒸散量分割為蒸騰潛力和土壤蒸發(fā)潛力：(2.75)Tp=ETp0-)=ETpSCFEp=ETp/g二FTpQ-SCF)其中ETp為潛在蒸散量(可以使用Penm

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

28、窗口雙擊RootWaterUptake-ModelsIE具條，彈出一個處理根系吸水模型的對話框（圖14）o水分脅迫模型鹽分脅迫模型根系吸鹽模型補償吸水域值+（臨界濕潤度）圖14HYDRUS-1D使用水分脅迫和鹽分脅迫模型處理根系的吸水。對于水分脅迫模型，計算公式為Ta-Js（九,.x）dx：二坊Ja（九（2.19）二Lr其中7；是潛在蒸騰量（cm/d）,7；是實際蒸騰量（cm/d）,S（0是吸水強度函數(shù)（cm/,注意x坐標(biāo)實際表示深度），？（力）是水分脅迫函數(shù)，力為土壤壓力水頭（cm）,b（x）是根系吸水分配（密度）函數(shù)，厶為根系層的深度。水分脅迫函數(shù)有2種經(jīng)驗表示方法，即Feddes模型和S

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

30、機制受到抑制，發(fā)生整體的水分脅迫，根系吸水總量將低于潛在蒸騰量，并正比于濕潤度。如果不考慮這種補償吸水機制，可以令2=lo14. 根系吸水一一水分脅迫參數(shù)在前處理窗口雙擊RootWaterUptake-WaterStressReduction工具條,彈出一個處理水分脅迫參數(shù)的對話框（圖15）o本算例中直接從數(shù)據(jù)庫中調(diào)入Wheat的經(jīng)驗值。脅迫函數(shù)r2Hr2L兩個臨界吸水強度數(shù)據(jù)庫吸水強度(cm/d)P3P2LP2HPoptPO壓力水頭(pressurehead,cm)圖15根系水分脅迫Feddes模型參數(shù)15. 輸入可變邊界條件的信息在前處理窗口雙擊VariableBoundaryCondit

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

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

33、需要仔細(xì)處理這些問題。XTimeVariableBoundaryConditions圖16表1時間(d)降水量(cm/d)ETp(cm/d)hCritA(cm):時間(d)降水量(cm/d)ETp(cm/d)hCritA(cm)LAI101000001601000002010000001000003010000018010000040100000-01000005010000020010000060100000210100000010000022010000081000002301000009010000024100000100100000251000001101000002610000012

34、0100000-10000013310000028010000014010000029010000015100000300100000數(shù)據(jù)可以先在Excel中準(zhǔn)備好，如表1。這些數(shù)據(jù)可以拷貝到圖16的電子表格中。這些數(shù)據(jù)顯示的降水量、蒸散潛力和葉面積指數(shù)變化特征如圖16右圖所示。在第17日由于莊稼收割，葉面及指數(shù)大幅度下降。16-編輯土壤剖面使用圖形界面在前處理窗口雙擊SoilProfile-GraphicalEditorL具條，程序?qū)棾鲆粋€處理土壤剖面的軟件(圖17)o這個圖形軟件的使用比較簡單，我們需要注意的是在Conditions菜單下面有很多子菜單，包含處理各種問題的功能菜單。首先要

35、做的事情，是確定把土壤剖面離散化為多少個節(jié)點。本算例土壤模型深度為3m,我們希望節(jié)點間距達到1cm,因此需要301個節(jié)點。選擇菜單Conditions/ProfileDiscretization,在下拉工具條中把Number修改為301。ConditionsProfileDiscretizationMaterialDistributionRootDistributionScalingfactorInitialConditionsConditions菜單剖分節(jié)點土壤巖性分層根系分布尺度因子初始條件子區(qū)分布觀察點ElHydrus-ID-Profile|IlnllX丿土堆剖而分成31個節(jié)點，則節(jié)點間

36、距為10圖17接下來，確定土層的分布，本模型有2個土層。缺省的土層編號為index=l,就是圖12中的細(xì)砂壤土層。選擇菜單Conditions/MaterialDistribution，在下拉工具條中使用Editcondition,把下部土壤層設(shè)置ConditionsMaterialDistribution1.點擊工具條2在剖面底部點擊鼠標(biāo)左鍵再放開，從下往上移動亂標(biāo)，選定土層范圍，然后再點擊鼠標(biāo)左鍵，修改土z坐標(biāo)高度設(shè)置2層土壤巖性編號index二2,這個土層編號為2,實際上就是圖12中的中砂層。圖18再就是確定根系隨深度的分布，根系分布函數(shù)方(勸是一個很特殊的函數(shù)，它滿足以下條件lob(x

37、)dx=l=工仇/=!其中LR為根系層厚度，x是深度；？z是節(jié)點間距，3是每個間距中的根系分布函數(shù)值，M是根系層占節(jié)點數(shù)。本算例中根系層的厚度為1m,假設(shè)根系分布函數(shù)為線性，并有b(x)=2x,y?l容易證明上式滿足積分為1的條件。設(shè)置方法為，選擇菜單Conditions/RootDistribution,在在下拉工具條中使用Editcondition,劃定根系層范圍(100cm),把頂部數(shù)值設(shè)置為0,底部數(shù)值設(shè)置為2,讓程序自動進行線性差值。ConditionsRootDistributiontiersLdit71護Ct-ticn?EpId使用工具條*RootslineardisTopval

38、ueDottonvalueOK口&區(qū)w9，根系分布函數(shù)9IZJConditionSpecifCancel圖19下一步處理初始條件，選擇菜單Conditions/InitialCondition。本算例中初始地下水位高于底板100cm,假設(shè)土壤剖面初始狀態(tài)是靜力平衡態(tài)，則模型底部的壓力水頭為+100cm,地面的壓力水頭為?200cm。使用EditCondition工具條，分別設(shè)置頂部和底部的壓力水頭，并讓程序自動插值形成初始條件。最后，選擇菜單Conditions/ObservationPoints,使用InsertX具條添加若干觀察點，有必要在靠近模型底部的位置加一個觀察點以判斷地下水

39、位的變化。關(guān)閉圖形程序退回主程序。17.編輯土壤剖面使用表格在前處理窗口雙擊SoilProfile-SummaryX具條，程序?qū)棾鲆粋€表格(圖20)o在這個表格中可以進一步修改土壤剖面數(shù)據(jù)。SoilProfileuamiryglnes7»%s5ex8*-cu8.4溢總總黑-.。圖2018.運行模型我們已經(jīng)把所有需要的數(shù)據(jù)都輸入了模型中，下面就可以運行模擬計算程序7o選擇菜單Calculation/ExecuteHydrus,在彈出的對話框中選擇OK,則會進入模擬計算模塊。軟件將調(diào)用HYDRUS-1D的核心程序進行計算，并輸出有關(guān)的信息,見圖21。Hydrus-IDCalculati

40、on:test6O?2"32229526SW72“X26力2QR271630S07刖212325S66666662222222213E«98-G.S4E«0i036E,*6i-0.i3E-9*««*S3?8.1胸14E*胸-9.54E電1g.36E«61-«.23E-9S»*w631?1陽14E«IMO.S4EH»1a.36IQle.23EdSi»«*S31Q.1IM14E*l>a-a.54E«01B.36E«01-a.23E-BS»

41、71;*w*w*-53ai.1»Ml-1E«»a-H.54EH91B.36E«fllHI.23E-BS*«M*w*-S2»7.IBM13E*BH-O.54F*ei036EtHY23E騎12E*86-0.G4E«eie.36E*61-«.?9E-«*MM<*M-S27316012E-0.54t»«l8.36E-61-«.23E-«-*-5242186TintHUItCuRvropSvTopEvRootSvBotnTophRocnPro”仍62752314&

42、;33?5%鎳陽U陌7172K7678tl絲竝也僦呎z2t.2?.2t.2-?.2?.ee?i29.mci2222232732)和於414345祁357342“666666G6C6622222222222contiiKi®13E*W-013E胸-013E»l>a-M13E*»a-W11E*WJuiee*B8-eiiE»ee-0liE.aaqliE*O3-«HE*»a-w-0.S4E«ei-$24854E«91S.36E*61-«.23E-Www*w-$24254E©8.36E*61-0.2

43、3E-«*-52J553E»01B.36E*01-«.23E-BS*iw“5：22753EH910.36E«fll-e.23E-BS*M*w*-S2l8B3EHJ10.36E-eiHI.23E-«52C75U»01936F01F23M»*w52353E仙0.36E«ei-e.23E-a5-*-517653E-9iD.36K*«i-a.23E-OS»M-5iii53EH910.36E«fll-e.23E-0S*«MMM*-Sil2，斑I1陽100IBMIBMiBU1胸1陽IMI

44、BMXBM圖21-般情況下，如果迭代計算的收斂性較好，模擬時期乂不是很長，將很快計算完畢。但是，一旦出現(xiàn)難以收斂的情況，則可能等待很長時間才有結(jié)果，或者非正常中斷。本算例模型在很短的時間內(nèi)就能夠完成計算。9察看結(jié)果模擬運算完畢之后，模擬結(jié)果可以通過后處理窗口中的各個工具條來察看，見圖22。MPHYDRIC"10-rtRtl回岡ObservationNodesfXKmbe.essingCalculationResultsOptionsWindowHelp-(5x10153D25ObservationNodes:WaterContent32oo.】22L1Ert>t|1!»

45、;«1后處理觀察點結(jié)果剖而曲線圖像邊界水分通雖水力特性曲線運行迭代信息物質(zhì)平衡信息ntBC它KootWaterUptake一Models粢RootWaterUptake-WaterStreF®VariableBoundaryConditions4刖.1團ForHelp,pressFlPost-processingOWObservationPoints帆)ProfileInforination、000迎貞號.FI。見.二.$0曲4瓠坯.卩1曲豁.4.J呱SoilHydraulicPropertiesWJRunTiiceInformation000MassBalanceInfo

46、rmation圖22模型模擬的結(jié)果包括土壤含水量、壓力水頭、實際根系吸水（蒸騰）量、實際地表通量、土壤水分存貯總量等的變化。20.輸出結(jié)果如果我們想把模擬結(jié)果保存為其它文件格式，用其它軟件進行后處理，則可以做進一步的轉(zhuǎn)換。HYDRUS-1D的結(jié)果圖包含曲線信息，點擊右鍵，會看到有一個EditChartData菜單，用這個菜單可以通過電子表格輸出模擬結(jié)果，見圖23。這個電子表格的個別數(shù)據(jù)可以修改后重新顯示圖像。圖像可以打印。051015202530TimeIdavsf-DcfeL'iltIPrintFrcviOusBcxtCloseKorirpalalVertical|Time|Act&

47、#39;jud.RootWterVptlscActualRootWaterUptakeBoundaryWaterFluxesandDataGridEdit*Re*心妙0.0Ic八111|艸o.»ro?5B2I002o.»ro?513Z).003o.»ro?54/J.DM3o.»ro?5S/ooceo00082o0OUoB80DU8o9ooiseoEditInsertdeleteGidScej8ca*«屈7圖23然而，圖23中的電子表格不具備輸出數(shù)據(jù)的功能。為此，你可以打開模型目錄（文件夾），找到以O(shè)UT為后綴的文件，圖24。這些文件都可以用記事

48、本打開察看數(shù)據(jù)。'test回岡文件（E）編斡空看邊收藏工具幫助（9存放模型的目錄叟需花文件夾1圍地址（E）|&C:laToolshydrlDProjcctslcst"|IH轉(zhuǎn)到文件夾6DSGMGShy©LJQ.0SkotLaniples»DirectHP!名稱iATKOSFttINSELECTOR,IN、A丄眈LOUT血MCE.OUTI.CfflO.OUTN3D.INF.OUTOBS.WDE.OUTPROFILE.OUTRUNNFOUTTLBWL.OUTfiHYDRUSW.DAT血PROFILEDATEDESCRIPT.TKT大小類型4KBIN文件3 KBIM文件4 KBOUT文件20KBOUT文件丄口nrrr七可以用記事本打開的結(jié)果文件33KBOUT文仁49KBOUT文件190KBOUT文件1KBVCD視頻文件36KBVCD視頻文件1KB文本文檔<><>選定8個對界1.45MBJ猶的電筋亡3Settings圖24有大量的計算結(jié)果保存在文件中，查閱HYDRUS-1D的用戶手冊，就知道文件保存了那些結(jié)果變量，例如根系吸水強度vROOT等等。附表給出了變量縮寫名稱的含義。Table13.1.TLEXEL.OUTpressureheadsandfluxesoiltheboundariesandint