HydrusD简明使用手册

1、DocumentserialnumberNL89WT-NY98YT-NC8CB-NNUUT-NUT108用HYDRUS-1D模拟剖面变饱和度地下水流（简明手册）王旭升中国地质大学（北京）HYDRUS-1D是一个共享专业软件，用于模拟一维变饱和度地下水流、根系吸水、溶质运移和热运移。本手册只介绍应用HYDRUS1D模拟垂向剖面水流和根系吸水的操作方法。1. 如何获取HYDRUS-1DHYDRUS-1D由位于欧盟捷克的PC-Progress工程软件开发公司发行，用户可以登录该公司首页：。为了下载HYDRUS-1D,应先注册成为用户，然后下载Hydrus-1D的安装文件：。这个文件对应目前HYDRU

2、S-1D的最高版本。2. 版权声明HYDRUS-1D的作者为：(1) J.Simunek,DepartmentofEnvironmentaiSciences,UniversityofCaliforniaRiverside,Riverside,California,USA.(2) M.Sejna,PCProgress,Prague,CzechRepublic.(3) .vanGenuchten,DepartmentofMechanicalEngineering,FederalUniversityofRiodeJaneiro,RiodeJaneiro,Brazil.感谢他们提供了一个如此精美而乂免

3、费使用的专业软件，帮助我们从事有关的科学和教育工作。当你运行解压文件后，会在您的电脑中产生一个安装目录，其中包含可执行文件。运行这个文件即可安装HYDRUS-1D软件。当您安装HYDRUS-1D时，象安装其它软件一样，会出现一个许可协议，从中可知本共享软件也受到美国法规的保护。3. 参考资料HYDRUS-1D安装之后，在软件运行目录下有HYDRS-1D文件。从这个文件您可以了解到HYDRUS-1D的一些技术细节，如水流、溶质运移、热流的方程、一些处理专门问题的模型、输入输出文件等等。有一个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创皂旧|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區旧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/模拟内容选项几何形状参数及剖面方式时间信息输出方式水流迭代求解控制参数水流土壤水分特征模型水流土壤水分特征曲线参数水流边界条件溶质运移-般信息溶质运移-一一运移参数溶质运移-一一化学反应参数溶质运移-一一边界条件根系吸水-一一模型根系吸水-一一水分胁迫函数可随时间变化的边界条件土壤剖而一一图形界面土壤剖而一一数据列表图25. 设计模型在使用HYDRUS-ID之前，您需要对饱和-非饱和水流模拟的基本原理有所了解，并设计出自己想做的模型

12、，准备好数据。一个剖面水流模型通常包含以下几个要素：(1) 土壤剖面从地面算起的深度，准备模拟那个时间段的水分变化。(2) 土壤分儿层，每层土壤的渗透性参数和水分特征曲线是怎样的。(3) 根系是怎么分布的。(4) 是否己经确定地面降雨入渗、蒸发蒸腾的信息，特别是它们随时间的变化。(5) 是否己经确定剖面底部的状态属于哪种类型的边界条件。下面是一个参考模型的设计图：6.使用HYDRUS-1D创建模型打HYDRUS-1D软件，选择”File/new”菜单，新建一个模型。在name栏中输入本模型的名称”test”,更改模型存放的目录。图4需要注意的是，HYDRUS-1D模型本身在计算机中就表现为一系

13、列的输入输出文件，它们存放在与模型名称一致的目录中。本例中，软件会自动创建一个需称为”test”的目录，而”CATOOLSHYDRlDProjectsw中除了test目录之外，还有一个文件。这是一个模型项目(project)文件，告诉软件下次到哪里去寻找模型。、模型创建之后，会显示前处理和后处理窗口(图5)。由于是新模型，还没有任何模拟结果，所以后处理窗口是空白的。图57. 输入模型控制信息首先，在前处理窗口双击MainProcesses,在弹出的对话框中输入模型的描述：atestmodel.然后在Simulate一栏中选中RootWaterUptake,表示想处理根系吸水问题。电击0K之后，

14、前处理窗口将增加处理根系吸水的工具条。图6下一步，是输入模型的儿何信息和土层划分信息。在前处理窗口双击GeometryInformation,在弹出的对话框中输入如图7所示的数据。模型有2个土层长度单位是cm土壤剖而的深V度是300cm图7接下來输入时间信息，在前处理窗口双击TimeInformation,会弹出一个对话框（图8）。TimcInformationlimeVniIsfineDiscretizatioCSecondsCDiriutesCKnursr2=FinalTineIniiialTimeStep|D0C1MinaniunTimeStep»001Mx2fTiWiTime

15、Step|1XnitialTime='上边界是随时间变化的，每天yI7组数据，共30<Timtt-Vori.abloBcurtiryCor»8it:cns7ime-VeriallsBoundaryCo&dilionsNinbci-o£T:ncVwiobieBozdzyRec&r4s(e.g.?模拟30d内的变化A时间步长信息组数据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小时内的变化图8这个对话框中提供了一些灵活的选项來处理上边界条件的变化，下面简要加以说明：(1) 蒸腾量的每日周期变化HYDRUS-1D可以使用一个经验

17、公式來处理每天24小时潜在蒸腾量的变化，设某天的潜在蒸腾量为(例如用Pemman公式获取的，cm/d)，则坊(0=0.247；7(/)=2.7578111'2冗1冗'、lday2；/<0.264d./>0.736d/c(0.264cl.0.736d)(2.73)其中7；(广)是瞬时潜在蒸腾量，广为时间。模型假设早上6点之前以及晚上18点-24点的蒸腾量总和只占全天蒸腾量的1%。注意本例中蒸腾量的单位是cm/d。(2) 降水量的周期变化如果在你的模型中降水量是周期性变化的，HYDRUS-1D也可以用一个公式來处理P(z)=P1+cos,Az丿丿(2.74)其中戸是周期

18、?t内的平均降雨量。(3) 使用气象数据也可以在HYDRUS-1D中输入气象数据，它将自动利用这些数据计算潜在蒸散量ETpo可以选择FA0组织推荐的Penman-Monteith公式，也可以选择Hargreaves公式。这些公式需要辐射、气温、湿度之类的气象数据。模型的另一个控制信息是对模拟结果的输出如何进行设置。在前处理窗口双击PrintInformationIE具条，弹出一个对话框。本例中确定输出30组模拟结果，每天输出1组。输出30组结V果，每天1组VPrintTnfonnationXP7-L<iv61TnfonitiitionOK£verynlinesteps：1Co&

19、amp;ccLPrint冬沪laTigIntwrvdrErevioxjs.j"Z厂7Scr6wxiOutputHwt.一PrintFluxes(insteadofl«ip)forObservationNodes*KtIEnteratEnd9II.PrintTinasUiiribwofPrintTinas可以确定需要输出哪些时间点的结果SelectFrintIarits.8. 水流模型一一迭代计算参数HYDRUS-1D是采用迭代法來处理非线性Richards方程的。在前处理窗口双击WaterFlow-IterationCriteriaE具条，弹出一个设置迭代参数的对话框（图1

20、0）o迭代控制参数的设置具有高度的专业技术性，除非特别了解，一般可以使用默认值。如果模拟结果出现不收敛的情况，需要对最大迭代次数、迭代精度等参数进行调整，但是在缺乏经验的情况下很难操作。IteratiomCriteria最多迭代次数含水虽迭代精度压力水头迭代秸度IterationCritcriaPressureHex!5a>:inuTiNanberofIterationsContentTolerance增大步长迭代次数信号缩小步长的迭代次数信号J步长增大比例I步长缩小比例1最小吸力间距lT最大吸力间距qlimeStepControl3OptimoLLt«rtionUpperOp

21、tiffialIterationRangeTim。St&pHuitiplication0.7VppcrTimeSlepflulliplicationFactorIntarpolticnLimitofthalnsionIntervalIOCCOVpperLimitofth«TensionInterv<J用于生成水分特征曲线的散点值图10HYDRUS-ID采用自动控制时间步长的方法來处理迭代的收敛性。对于每个时步，如果迭代次数太多，就缩小时间步长；如果没经过儿次迭代就达到收敛精度，则适当增大时间步长。9. 水流模型一土壤水力特性模型水分特征曲线是非饱和土壤的重要物理性质，H

22、YDRUS-1D提供了儿种方法來处理与之有关的参数。在前处理窗口双击WaterFlow-SoilHydraulicProperties工具条，弹出一个设置水力特性模型的对话框（图11）。Soi1HydraulicModel单孔介质模型HoldSinaiaPoroxityModal*&van-V*Q41en厂Mth.ir-EnurV<J,nc&f-2cnJlodiEj«dvanGenuchtenCBrooks-Corwj*K&zudQcfi-xiomAl)双重介质模型双重孔隙度双重渗透性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主要在溶质运移问题中使用Myiteresis吸湿和疏十滞后过程模型e3?0hystorasisJlysteresjsmretenlioacurveCH'starwKKr«tAntioacurveandconductjvtty'fctrczds5nrctantiacutvqopwpinc,BobIhird)CIrntLilly空yi班.cor"C*H:tidilyerve图11在一般情况下，选择单孔介质模型，并选择用vanGenuchten-Mua1em公式处理土壤的水力特性就可以了。如果还要模拟溶质运移

25、，可能需要考虑双重介质模型。双重介质在同一个点有两个孔隙度或两个渗透率，相当于两种介质的混杂。双重介质模型能够模拟这两种“介质”之间的水分和盐分交换。10. 水流模型土壤水分特征曲线在前处理窗口双击WaterFlow-SoilHydraulicParametersX具条，弹出一个设置水分特征曲线参数的对话框（图12）o本例中选择vanGenuchten-Mua1em公式处理水分特征曲线，务+跻0/<07/>0(2.28)(2.29)其中222,/均为控制因子。HYDRUS-1D软件中提供了一组土壤经验参数库，可供用户参考。本例中两层土壤的参数直接从数据库中调出：第1层对应Sandy

26、loam,第二层对应sando图12在输入参数时，请注意参数的单位。11. 水流模型一一边界条件在前处理窗口双击WaterFlow-BoundaryConditionsX具条，弹出一个设置边界条件的对话框(图13)o图13上边界条件有6种类型，下边界条件有8种类型。边界类型的确定需要考虑实际条件，在本算例中，上边界选择大气边界条件，在降雨量很大时地表可以产生积水。植被蒸腾量和土壤蒸发量分开处理，HYDRUS-1D推荐使用一个经验公式来把潜在蒸散量分割为蒸腾潜力和土壤蒸发潜力：(2.75)Tp=ETp0-)=ETpSCFEp=ETp/g二FTpQ-SCF)其中ETp为潜在蒸散量(可以使用Penm

27、an-Monteith公式处理气象数据得到，cm/d),Tp为潜在蒸腾量(cm/d),Ep为土壤潜在蒸发量(cm/d),LAI是叶面积指数，R为消光系数，取决于太阳角度、植被类型及叶片空间分布特征。SCF是一个中间参数，即土壤覆盖度(Soilcoverfraction)o在阔叶植被发育的情况下，消光系数的经验值为扫。12. 水流模型一一定水头或通量边界设置如果边界条件中包含定水头或定通量的边界，则在前处理窗口双击WaterFlow-ConstantBC工具条，弹出一个设置边界数据的对话框。本算例模型中，下边界为定流量边界，实际上就是隔水边界，因此直接输入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是吸水强度函数（cm/,注意x坐标实际表示深度），？（力）是水分胁迫函数，力为土壤压力水头（cm）,b（x）是根系吸水分配（密度）函数，厶为根系层的深度。水分胁迫函数有2种经验表示方法，即Feddes模型和S

29、-Shape模型。Feddes模型是一个梯形函数，只需要知道力值。而S-Shape模型把水分胁迫和叶片气孔的压力水头联系起來，需要知道气孔压力水头的数值力”。土壤的湿润度可以表示为co=TaITlt=ja（h,x）b（x）dxLR但是如果直接用这种方法来计算实际蒸腾量有一定的问题。植被其实可以调节不同深度的水分胁迫响应特征；某个深度土壤干燥吸不上水，植被可以加大在比较湿润的土层的吸水量，以补偿不足。这种现象称为补偿吸水。为了模拟根系补偿吸水，HYDRUS-1D提供了一种简化的模型，即如果湿润度高于某个临界值（62）,植被根系可以通过补偿机制充分吸水达到潜在蒸腾量。如果湿润度低于这个临界值，补偿

30、机制受到抑制，发生整体的水分胁迫，根系吸水总量将低于潜在蒸腾量，并正比于湿润度。如果不考虑这种补偿吸水机制，可以令2=lo14. 根系吸水一一水分胁迫参数在前处理窗口双击RootWaterUptake-WaterStressReduction工具条,弹出一个处理水分胁迫参数的对话框（图15）o本算例中直接从数据库中调入Wheat的经验值。胁迫函数r2Hr2L两个临界吸水强度数据库吸水强度(cm/d)P3P2LP2HPoptPO压力水头(pressurehead,cm)图15根系水分胁迫Feddes模型参数15. 输入可变边界条件的信息在前处理窗口双击VariableBoundaryCondit

31、ionsJ2具条，弹出一个处理时间序列数据的对话框(图16)o本算例中，在步骤(11)中已经把地面处理大气边界，同时乂选择使用消光系数法划分植被蒸腾和土面蒸发，因此需要输入每天的降水、潜在蒸散量、叶面积指数等数据。还有一个需要输入的数据是最小压力水头值，即地面土壤达到最干燥状态时的压力水头。从理论上讲，当土壤十分干燥时，吸力很大，而液态孔隙水的压强很小，与空气湿度保持平衡关系，因此有H,.=exp一组響财(2.72)寸比)其中九为最小压力水头，忆为空气绝对湿度，RT/Mg为空气的摩尔气体常数。空气湿度虽然可以通过气象数据得到，但这里公式需要的是近地面的空气湿度。一般情况下，取饱和水汽湿度是可取

32、的，因为2cm深度以下土壤空气的湿度往往都是饱和的，只不过随温度发生变化。因此，可以根据近地面气温的变化來推算地表土壤的空气湿度(饱和水汽湿度),再换算成压力水头。HYDRUS-1D中需要输入的是最小压力水头的绝对值，缺省值为hCritA=|=106cm=10'm这个数值只会对土壤蒸发起作用。HYDRUS-1D建议：hCritA所对应的土壤含水量应该至少比残余含水量大，在模拟根系吸水的情况下，人还应该低于图15中的P3o否则(九P3),当根系吸水的临界值压力水头(P3)和地面蒸发的最小压力水头(九)满足时，会导致回流(inflow)现象，这是不合理的。除非存在特别干燥的情况，模型一般不

33、需要仔细处理这些问题。XTimeVariableBoundaryConditions图16表1时间(d)降水量(cm/d)ETp(cm/d)hCritA(cm):时间(d)降水量(cm/d)ETp(cm/d)hCritA(cm)LAI101000001601000002010000001000003010000018010000040100000-01000005010000020010000060100000210100000010000022010000081000002301000009010000024100000100100000251000001101000002610000012

34、0100000-10000013310000028010000014010000029010000015100000300100000数据可以先在Excel中准备好，如表1。这些数据可以拷贝到图16的电子表格中。这些数据显示的降水量、蒸散潜力和叶面积指数变化特征如图16右图所示。在第17日由于庄稼收割，叶面及指数大幅度下降。16-编辑土壤剖面使用图形界面在前处理窗口双击SoilProfile-GraphicalEditorL具条，程序将弹出一个处理土壤剖面的软件(图17)o这个图形软件的使用比较简单，我们需要注意的是在Conditions菜单下面有很多子菜单，包含处理各种问题的功能菜单。首先要

35、做的事情，是确定把土壤剖面离散化为多少个节点。本算例土壤模型深度为3m,我们希望节点间距达到1cm,因此需要301个节点。选择菜单Conditions/ProfileDiscretization,在下拉工具条中把Number修改为301。ConditionsProfileDiscretizationMaterialDistributionRootDistributionScalingfactorInitialConditionsConditions菜单剖分节点土壤岩性分层根系分布尺度因子初始条件子区分布观察点ElHydrus-ID-Profile|IlnllX丿土堆剖而分成31个节点，则节点间

36、距为10图17接下來，确定土层的分布，本模型有2个土层。缺省的土层编号为index=l,就是图12中的细砂壤土层。选择菜单Conditions/MaterialDistribution，在下拉工具条中使用Editcondition,把下部土壤层设置ConditionsMaterialDistribution1.点击工具条2在剖面底部点击鼠标左键再放开，从下往上移动亂标，选定土层范围，然后再点击鼠标左键，修改土z坐标高度设置2层土壤岩性编号index二2,这个土层编号为2,实际上就是图12中的中砂层。图18再就是确定根系随深度的分布，根系分布函数方(劝是一个很特殊的函数，它满足以下条件lob(x

37、)dx=l=工仇/=!其中LR为根系层厚度，x是深度；？z是节点间距，3是每个间距中的根系分布函数值，M是根系层占节点数。本算例中根系层的厚度为1m,假设根系分布函数为线性，并有b(x)=2x,y?l容易证明上式满足积分为1的条件。设置方法为，选择菜单Conditions/RootDistribution,在在下拉工具条中使用Editcondition,划定根系层范围(100cm),把顶部数值设置为0,底部数值设置为2,让程序自动进行线性差值。ConditionsRootDistributiontiersLdit71护Ct-ticn?EpId使用工具条*RootslineardisTopval

38、ueDottonvalueOK口&区w9，根系分布函数9IZJConditionSpecifCancel图19下一步处理初始条件，选择菜单Conditions/InitialCondition。本算例中初始地下水位高于底板100cm,假设土壤剖面初始状态是静力平衡态，则模型底部的压力水头为+100cm,地面的压力水头为?200cm。使用EditCondition工具条，分别设置顶部和底部的压力水头，并让程序自动插值形成初始条件。最后，选择菜单Conditions/ObservationPoints,使用InsertX具条添加若干观察点，有必要在靠近模型底部的位置加一个观察点以判断地下水

39、位的变化。关闭图形程序退回主程序。17.编辑土壤剖面使用表格在前处理窗口双击SoilProfile-SummaryX具条，程序将弹出一个表格(图20)o在这个表格中可以进一步修改土壤剖面数据。SoilProfileuamiryglnes7»%s5ex8*-cu8.4溢总总黑-.。图2018.运行模型我们已经把所有需要的数据都输入了模型中，下面就可以运行模拟计算程序7o选择菜单Calculation/ExecuteHydrus,在弹出的对话框中选择OK,则会进入模拟计算模块。软件将调用HYDRUS-1D的核心程序进行计算，并输出有关的信息,见图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-般情况下，如果迭代计算的收敛性较好，模拟时期乂不是很长，将很快计算完毕。但是，一旦出现难以收敛的情况，则可能等待很长时间才有结果，或者非正常中断。本算例模型在很短的时间内就能够完成计算。9察看结果模拟运算完毕之后，模拟结果可以通过后处理窗口中的各个工具条來察看，见图22。MPHYDRIC"10-rtRtl回冈ObservationNodesfXKmbe.essingCalculationResultsOptionsWindowHelp-(5x10153D25ObservationNodes:WaterContent32oo.】22L1Ert>t|1!»

45、;«1后处理观察点结果剖而曲线图像边界水分通虽水力特性曲线运行迭代信息物质平衡信息ntBC它KootWaterUptake一Models粢RootWaterUptake-WaterStreF®VariableBoundaryConditions4刖.1團ForHelp,pressFlPost-processingOWObservationPoints帆)ProfileInforination、000迎貞號.FI。見.二.$0曲4瓠坯.卩1曲豁.4.J呱SoilHydraulicPropertiesWJRunTiiceInformation000MassBalanceInfo

46、rmation图22模型模拟的结果包括土壤含水量、压力水头、实际根系吸水（蒸腾）量、实际地表通量、土壤水分存贮总量等的变化。20.输出结果如果我们想把模拟结果保存为其它文件格式，用其它软件进行后处理，则可以做进一步的转换。HYDRUS-1D的结果图包含曲线信息，点击右键，会看到有一个EditChartData菜单，用这个菜单可以通过电子表格输出模拟结果，见图23。这个电子表格的个别数据可以修改后重新显示图像。图像可以打印。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中的电子表格不具备输出数据的功能。为此，你可以打开模型目录（文件夹），找到以OUT为后缀的文件，图24。这些文件都可以用记事

48、本打开察看数据。'test回冈文件（E）编斡空看边收藏工具帮助（9存放模型的目录叟需花文件夹1圍地址（E）|&C:laToolshydrlDProjcctslcst"|IH转到文件夹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七可以用记事本打开的结果文件33KBOUT文仁49KBOUT文件190KBOUT文件1KBVCD视频文件36KBVCD视频文件1KB文本文档<><>选定8个对界1.45MBJ犹的电筋亡3Settings图24有大量的计算结果保存在文件中，查阅HYDRUS-1D的用户手册，就知道文件保存了那些结果变量，例如根系吸水强度vROOT等等。附表给出了变量缩写名称的含义。Table13.1.TLEXEL.OUTpressureheadsandfluxesoiltheboundariesandint