Hydrus1D简明使用手册

1、用HYDRUS-1D模拟剖面变饱和度地下水流（简明手册）王旭升中国地质人学（北京）目录1 .如何获取 HYDRUS-1D22 .明23 .参考资料24 . HYDRUS-1D 的 WINDOWS 界面25 .设计模型36,使用HYDRUS-1D创建模型47 .输入模型控制信息58 水流模型一迭代计算参数69水流模型土壤水力特性模型710水流模型一土壤水分特征曲线711 .水流模型边界条件812 .水流模型一定水头或通量边界设置913 .根系吸水-吸水模型914 .根系吸水一水分胁迫参数1015 .输入可变边界条件的信息1016 .编辑土壤剖面一使用图形界面1117 .编辑土壤剖面使用表格141

2、8 .运行模型1419 .察看结果1520 .输出结果15HYDRUS-1D是一个共享专业软件，用于模拟一维变饱和度地下水流、根系吸水、溶质运移和热 运移。本手册只介绍应用HYDRUS1D模拟垂向剖面水流和根系吸水的操作方法。1 .如彳腺取HYDRUS-1DHYDRUS-1D由位于欧盟捷克的PC-Pogess工程软件开发公司发行，用户可以登录该公司首 页：o为了下载HYDRUS-1D,应先注册成为用户，然后下载Hydus-1D的安装文件：H1D_44. exe。这 个文件对应目前HYDRUS-1D的最高版本。2 .版权声明HT0RUS-1D的作者为：(1) J. Simunek, Depart

3、ment of Environmental Sciences, University of California Riverside, Riverside, California, USA(2) M Sejnaf PC Progress, Prague, Czech Republic (3) M Th. van Genuch ten, Depar to f Mechanical Engineering, Federal Universi ty of Rio de Janeiro, Rio de Janeiro, Brazil 感谢他们提供了一个如此精美而又免费使用的专业软件，帮助我们从事有关的

4、科学和教育工 作。当你运行HID_4_14.exe解压文件后，会在您的电脑中产生一个安装目录，其中包含Setup, exe可 执行文件。运行这个文件即可安装HYDRUS-1D软件。当您安装HYDRUS-1D时，象安装其它软件一样，会出现一个许可协议，从中可知本共享软件也 受到美国法规的保护。3 .参考资料HTORUS-1D安装之后，在软件运行目录卞有HYDRS-IDManual.pdf文件。从这个文件您可以了解到 HYDRUS-1D的一些技术细节，如水流、溶质运移、热流的方程、一些处理专门问题的模型、输入输 出文件等等。有一个Examples目录，包含大量的模拟算例可供参考。用户还可以参考以下

5、文献： Simunek, J Jh. van Genuchtenr and Sejna, Development and applications of the HYDRUS and STANMOD software packages, and related codes, Vadose Zone Journal, doi:10 - 2136/VZJ2007 - 0077, Special Issue n Vadose Zone Modeling", 7(2), 587-600, 200& Jacques, D ,J - Simunek, D Mallants, and M

6、Th. van Genuchten, Modeling coupled hydrological and chemical processes: Long-term uranium transport following mineral phosphorus fertilization, Vadose Zone Journal> doi:10. 2136/VZJ2007, 0084, Special Issue “ Vadose Zone Modeling ,7(2), 698-711, 2008. Simunek, J and Th. van Genuchten, Modeling n

7、onequilibrium flow and transport with HYDRUS, Vadose Zone Journal> doi:10.2136/VZJ2007. 0074, Special Issue " Vadose Zone Modeling", 7(2), 782-797, 2008.这些文献都可以从下载。4 . HYDRUS-1D 的 WINDOWS 界面运行HYDRUS-1D,可以看到一个Windows的界面如卞: HYDRUS ID ROOTUPTK刘丽IViw Ere-prAcessing Calculation Eesuits ”川：j、

8、：二陶旧 i 创皂旧|rn| ROOTUPTKPre-processing iX flain ProcessesS/ Ceoroetry Inforimtion 0/Time Inforwation前处理工具Pogt-X»rocessirObservation Points模拟计算'一台处理工具Print:二 1Water : CZBopyWitoartSh driteriaWater:Flow - Soil Hydraulic Property flodelSolutu ®wl- SolLtd£SUlic PazanA-tersSolute 磴 1 So

9、lBtew - Boundary Conditions:Transport - General Infornm-tion! Transport Transport Paraneters Transport - Reaction PazanA-ters :Transport Boundary Conditions 璧 Root Water Intake Kodels；令 Root Water Uptake - Water Stress Reaction F® Variable Boundary Conditions 皈Soil Profile - Graphical EditorI 切

10、 I Soil Profile - SurnraryPot: Help, press FlW Profi 1 e InforiTia-tion Flow - Boundary Fluxes and Heads Solute Transport - Actual and Cunfjlarive Boundary Fluxes Soil«Q Mjrdraulic PropertiesRun Tins Information flass Balance nfoixiatiem00TO图1所有的前后处理在界面中一目了然，左边是前处理工具，右边是后处理工具。其中前处理的 各项功能如下图所示。

11、HYDRUS ID ROOTUPTKEile Viw Ere-processing Calculation EesultDie陶旧I创皂旧|m|刘丽I ROOTUPTKPre-processingiX flain ProcessesS/ Ceoroetry Inforimtion 0/Time InforwationPrint ElolbjnWtitai:二 1Water:匚二FjoWa跳应tien CriteriaWater:Flow - Soil Hydraulic Property flodelSolutu 磅的raulic PazanA-tersSolute 磴 1 Sol6tew -

12、 Boundary Conditions :Transport - General Infornm-tion ! Transport Transport Paraneters Transport - Reaction PazanA-ters :Transport Boundary Conditions 璧 Root Water Intake Kodels；令 Root Water Uptake - Water Stress Reaction F® Variable Boundary Conditions 皈Soil Profile - Graphical Editor (包 1 So

13、il Profile - Stnnnary几何形状参数及剖而方式时间信息输出方式水流迭代求解控制参数水流土壤水分特征模型水流 土壤水分特征曲线参数水流边界条件溶质运移-般信息溶质运移-运移参数溶质运移-化学反应参数 溶质运移-边界条件 根系吸水一一模型根系吸水一一水分胁迫函数可随时间变化的边界条件土壤剖而图形界面土壤剖而一一数据列表Pot: Help, press Fl5 .设计模型在使用HYDRUS-1D之前，您需要对饱和-非饱和水流模拟的基本原理有所了解，并设计出自己 想做的模型，准备好数据。一个剖面水流模型通常包含以下几个要素：Q)土壤剖面从地面算起的深度，准备模拟那个时间段的水分变化。

14、土壤分几层，每层土壤的渗透性参数和水分特征曲线是怎样的。根系是怎么分布的。 是否已经确定地面降雨入渗、蒸发蒸腾的信息，特别是它们随时间的变化。是否已经确定剖面底部的状态属于哪种类型的边界条件。降水量F面是一个参考模型的设计图:根系层1.细砂壤2.中砂-初始潜水面隔水底板6.使用HYDRUS-1D创建模型打开HYDRUS-1D软件，选择 File/new”菜单,新建一个模型。在name一栏中输入本模型的名称 test",更改模型存放的目录。图4需要注意的是，HYDRUS-1D模型本身在计算机中就表现为一系列的输入输出文件,它们存放 在与模型名称一致的目录中。本例中，软件会自动创建一个名

15、称为 test”的目录，而 CATOOLSHYDRIDProjectsn中除了 test目录之外，还有一个test, hid文件。这是一个模型项目（project） 文件，告诉软件下次到哪里去寻找模型。模型创建之后，会显示前处理和后处理窗1 （图5 ）。由于是新模型，还没有任何模拟结果， 所以后处理窗II是空白的。7.输入模型控制信息首先，在前处理窗II双击Main Processes,在弹出的对话框中输入模型的描述：a test model.然后在 Simulate 一栏中选中Root Water Uptake表示想处理根系吸水问题。电击0K之后，前处理窗II将增加处理根系吸水的工具条。Ma

16、in Processesng :|a test nodelSinulfttaInverse SolutionI - Vapor Row| XycroLogy Solute Trasport A General SoluteIcn Chemistryr KR (PHIEEQC)Heat Iransrcrt7例 5erVpt ctke :Root Srovih032 Iranzpon下一步，是输入模型的几何信息和土层划分信息。在前处理窗1 I双击Geometry Information,在 弹出的对话框中输入如图7所示的数据。模型有2个土层V长度单位是cm土壤剖而的深V度呈300 cm接卜来输入时

17、间信息，在前处理窗II双击Time Information,会弹出一个对话框（图8 ）。时间单位V上边界是随时间变 化的,每天一组数 据，共30组数据自动处理蒸瞌品在每天24小时内的模拟30 d内的A时间步长变化信息这个对话框中提供了一些灵活的选项来处理上边界条件的变化，卜面简要加以说明：(1)蒸腾量的每日周期变化HT0RUS-1D可以使用一个经验公式来处理每天24小时潜在蒸腾量的变化,设某天的潜在蒸腾量为打(例如用Pemman公式获取的,cm/d ),则坊("二0 24 坊/ <0.264d,/>0.736d(2 加乃(2 73)7八(/) = 2,757 ; sin

18、/w(0 - 264d.O - 736d).lday 2 )其中& )是瞬时潜在蒸腾量,F为时间。模型假设早上6点之前以及晚上18点-24点的蒸腾量总 和只占全天蒸腾量的1%。注意本例中蒸腾量的单位是cm/d。(2 )降水量的周期变化如果在你的模型中降水量是周期性变化的，HYDRUS-1D也可以用一个公式来处理玖 7) =P1 + COS(2.74)7Options7 1-Lrvel Infomat ioox可以确定需要输出哪输出30组结V 果，每天1组steps : Priie cal orvScreenQin .putPrint Fluxes 一 一linztaad ofT

19、1;np) for ObzcrvatknIf rt Pnta r at Urvi,"Uunbar of Print Tinas 总S«L<ct Print T:rks.line IntendC »CQLFrcvf OH .些时间点的结果J其中尹是周期At内的平均降雨量。(3 )使用气象数据也口 I以在HTORUS-1D中输入气象数据，它将自动利用这些数据计算潜在蒸散量ETpe可以选 择FAO组织推荐的Penman-Monteith公式，也可以选择Hargreaves公式。这些公式需要辐射、气温.湿 度之类的气彖数据。模型的另一个控制信息是对模拟结果的输出如何

20、进行设置。在前处理窗II双击Print Information 工具条，弹出一个对话框。本例中确定输出30组模拟结果，每天输出1组。Print Inforniation8 .水流模型迭代计算参数HTORUS-1D是采用迭代法来处理非线性Richards方程的。在前处理窗11双击Water FloA Iteration Criteria工具条，弹出一个设置迭代参数的对话框(图10 )。迭代控制参数的设置具有高度的专业技 术性，除非特别了解，一般可以使用默认值。如果模拟结果出现不收敛的情况，需要对最人迭代次数、迭代精度等参数进行调整，但是在缺乏经验的情况下很难操作。最小吸力间距L J 最大吸力间距

21、XIxitcmaL latepolatkn tjlq：)* CC6 LftvAr Limit of the lension Interv<l(1CCCOVpper Limit of th« Tension Iftterv<J用于生成水分特征曲线的散点值在前处理窗11双击Water Flow Soil Hydraulic Parameters工具条,弹出一个设置水单孔介质模型双重介质模型双重孔隙度Hyduulfc HoddSK>CLC Forosity Vod?ls 金 vux 也V*Q<il(nFfVnh Air-EntrA Xlu®offvoa G

22、ea'khte&f grooks-Corey(KosuAQog-ncmU.)IhKU-rvrosi ty/Duwl-P a/-(76d)ilf ty VoddzfDaal-porosityGhun "dnzd Genucht*n - Hmlei) , Dzl-poyogily Sotilirffnobi 达皿3c,gU DrvLL-porAsrty "'obSle-iTGobil he&d riAstransier) =Models be "arerec emend* 8 only for experienced 厂 -pcy ,

23、c。81伊宜。037的皿？3(11 VMl 6ewfcht?nx 1333)主要在溶质运移问题中使用图10HYDRUS-1D采用自动控制时间步长的方法来处理迭代的收敛性。对于每个时步，如果迭代次 数太多，就缩小时间步长；如果没经过几次迭代就达到收敛精度，则适当增人时间步长。9 .水流模型一土壤水力特性模型水分特征曲线是非饱和土壤的重要物理性质，HYDRUS-1D提供了几种方法来处理与之有关的 参数。在前处理窗II双击Water Flow Soil Hydraulic Properties工具条，弹出一个设置水力特性模型的对 话框（图11）。Soil Hydraulic M;,LocImpTab

24、les吸湿和疏卡滞后过程模型冷 3k hystere S|yst«ir* K)s g ratantioa curwQ lfystArasis io. rAtouUAa curv owf ggdactRityC fcttresdsdo. retentioo. curve (0.0 pwpiiuL Bob L»Lird) r . :id：.y 也Eg lotCK :tidily"tflwe rjrve图11在一般情况下,选择单孔介质模型,并选择用van Genuchten-Mualem公式处理土壤的水力特性 就可以了。如果还要模拟溶质运移，可能需要考虑双重介质模型。双

25、重介质在同一个点有两个孔隙度或 两个渗透率，相当于两种介质的混杂。双重介质模型能够模拟这两种介质之间的水分和盐分 交换。10.水流模型,土壤水分特征曲线分特征曲线参数的对话框（图12 ） o本例中选择van Genuchten-Mua 1 em公式处理水分特征曲线，/<0（2.28）7?>0(2.29)其中G n, /均为控制因子。HYDRUS-1D软件中提供了一组土壤经验参数库，可供用户参考。本例中两层土壤的参数直接从数据库中调出：第1层对应Sandy loam,第二层对应sand。细砂壤土中砂图12在输入参数时，请注意参数的单位。11 .水流模型一一边界条件在前处理窗II双击W

26、ater Flow- Boundary Conditions工具条，弹出一个设置边界条件的对话框（图。Water Flow Boundary Condi 地面边界类型定压力水头定水分通量大气边界，可积水大气边界，产流变水头下端边界类型定压力水头定水分通虽变地下水位变流虽I勺山下渗排水渗出而水平排水Upper Boundary Condition1 Constant Pressure HeadConstAnt AumQ Atmospheric BC vith Surface Layer1 Atmospheric BC vith Surface Run Off' Variable Pres

27、sure HeadLower Boundary Condition 'Constant Pressure Head GCon 泄ant Flux1 Variable Pressure Headc Variable Flux1 Free Drainage' Deep Drainage* Saapaga; I' HorLzontai DrainsInitial ConditionS In Pressm-e Heads初始条件使用压力水头使用含水量rAtmospkeric BC大气边界基于ETp和LAI消光系数最大积水厚度图13上边界条件有6种类型，下边界条件有8种类型。边界

28、类型的确定需要考虑实际条件，在本 算例中，上边界选择人气边界条件，在降雨量很人时地表可以产生积水。植被蒸腾量和土壤蒸发 量分开处理，HT0RUS-1D推荐使用一个经验公式来把潜在蒸散量分割为蒸腾潜力和土壤蒸发潜力：P P'T=ETA | -£肛 3 BTSCF7P(2.75)Ep =二 £7A(1 SCF)图14Ta 卜似胪)八：LR(力九(2.19)其中ETp为潜在蒸散量(口 J以使用Penman-Monteith公式处理气象数据得到，cm/d), Tp为潜 在蒸腾量 (cm/d), Ep为土壤潜在蒸发量(cm/d), LAI是叶面积指数，R为消光系数，取决于太阳

29、角度、植被类型及 叶片空间分布特征。SCF是一个中间参数，即土壤覆盖度(Soil cover fraction)。在阔叶植被发育的情况下， 消光系数的经验值为0 . 5-0 . 75。12 .水流模型定水头或通量边界设置如果边界条件中包含定水头或定通量的边界，则在前处理窗I I双击Water Flow Constant BC工具 条，弹出一个设置边界数据的对话框。本算例模型中，下边界为定流量边界，实际上就是隔水边界，因此直接输入0即可。13 .根系吸水吸水模型在前处理窗I I双击Root Water Uptake- Models工具条,弹出一个处理根系吸水模型的 对话框(图 14)。水分胁迫模

30、型盐分胁迫模型补偿吸水域值+(临界湿润度)根系吸盐模型HYDRUS-1D使用水分胁迫和盐分胁迫模型处理根系的吸水。对于水分胁迫模型，计算公式为其中7;是潜在蒸腾量(cm/d), 7 ;是实际蒸腾量(cm/d), SCO是吸水强度函数(cm/(cm. d),注意用坐标实际 表示深度)，从力)是水分胁迫函数，力为土壤压力水头(cm),方(£是根系吸水分配(密度)函数，公为根 系层的深度。水分胁迫函数有2种经验表示方法，即Feddes模型和S-Shape模型。Feddes模型是一个 梯形函数，只需要知道力值。而S-Shape模型把水分胁迫和叶片气孔的压力水头联系起来，需要知 道气孔压力水头

31、的数值加。土壤的湿润度可以表示为co = Ta/Tp 二 J a （ /?, x/ bx dx（ 2. 22 ）LR但是如果直接用这种方法来计算实际蒸腾量有一定的问题。植被其实可以调节不同深度的水分胁 迫响应特征；某个深度土壤干燥吸不上水，植被可以加人在比较湿润的土层的吸水量，以补偿不 足。这种现彖称为补偿吸水。为了模拟根系补偿吸水，HYDRUS-1D提供了一种简化的模型，即如果湿润度高于某个临界值 （。纵），植被根系可以通过补偿机制充分吸水达到潜在蒸腾量。如呆湿润度低于这个临界 值，补偿机制受到抑制，发生整体的水分胁迫，根系吸水总量将低于潜在蒸腾量，并正比于湿润 度。如果不考虑这种补偿吸水机

32、制，可以令型二1。14 .根系吸水一水分胁迫参数胁迫 函数痔P3数据库吸水强度（cm/d）&Hr32L在前处理窗11双击Root Water Uptake- Water Stress Reduction工具条,弹出一个处理水分胁迫参数 的对话框（图15 ）。本算例中直接从数据库中调入Whst的经验值。-1!_L.P曲 P2LP2H PoplpPopt PO压力水头（晔溶q0卜般4娥-'两个临界吸水强度图15根系水分胁迫Feddes模型参数15 .输入可变边界条件的信息在前处理窗II双击Variable Boundary Conditions T具条，弹出一个处理时间序列数据的对

33、话框（图 16 )e本算例中，在步骤（11）中已经把地面处理大气边界，同时又选择使用消光系数法划分植被蒸 腾和土面蒸发，因此需要输入每天的降水、潜在蒸散量、叶面积指数等数据。还有一个需要输入的数据是最小压力水头值，即地面土壤达到最干燥状态时的压力水头。从理 论上讲，当土壤十分干燥时，吸力很人，而液态孔隙水的压强很小，与空气湿度保持平衡关系，因 此有(2.72)其中九为最小压力水头，忆为空气绝对湿度，RT/Mg为空气的摩尔气体常数。空气湿度虽然可以通 过气彖数据得到，但这里公式需要的是近地面的空气湿度。一般情况下，取饱和水汽湿度是可取 的，因为2 cm深度以下土壤空气的湿度往往都是饱和的，只不过

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

35、riable Boundary Conditions(=J 年/KfiUcmd)LAIETpcmd )05-It_!t_I_,*-Al»1»-AlI-_I_r_r_r，' 户 声口| 414 Line gl.S LiajTime (dj图16时间(d)降水fit(cm/d)ETp(cm./ d)hCritA(cm)时间(d)降水量(cm/d)ETp(czn/d)hCritA(cm)LAI100.51000002100.41000002. 1200.61000002100.61000000.3300. 11000002.21800.71000000.3100.5100

36、0002.21900.71000000.3500.51000002.22000.71000000.3600.51000002.22100.71000000.3700.31000002.32200.61000000.380.70. 11000002.32300.51000000.3900. 11000002.3240.30.21000000.31000.51000002325230.11000000.31100.61000002.2261.10.11000000.31200. 11000002.20.50.11000000.31330.11000002.22800.41000000.31100.

37、5100000212900.71000000.3150.50.2100000213000.71000000.3数据可以先在Excel中准备好，如表Io这些数据可以拷贝到图16的电子表格中。这些数据显示的降水量、蒸散潜力和叶面积指数变化特征如图16右图所示。在第17 口由于庄稼收割，叶面及指数大幅度下降。16 .编辑土壤剖面一一使用图形界面这个图形软件的使用比较简单，我们需要注意的是在Conditions菜单有很多子菜单，包含处理各种问题的功能菜单。首先要做的事情，是确定把土壤剖面离散化为多少个节点。本算例土壤模型深度为3 m,我们希望节点间距达到1 cm,因此需要301个节点。选择菜单Cond

38、itions/ Profile Discretization,在卜拉工具条中把Number修改为301。接下来，确定土层的分布,ConditionsProfile Discretization MaterialDistribution Root DistributionScaling factor Initial Conditions SubregionsConditions 菜 单 剖分节点土壤岩性分 层根系分布尺度因子 初始条件子区分布观 率占土壤剖面分I成31 个节点，/则节点间距为 10 cm本模型有2个土层。缺省的土层编 号为index二1,就是图12中的细砂 壤土层。选择菜单Con

39、ditions/ Material Distribution,在下拉工具条中使用Edit condition,把下部土壤层设置 index"这个土层编号为2,实际上就是图12中的中砂层。ConditionsMaterial Distribution1 .点击工具条KHydrus-ID Profile 0A44eC ondi t i ons Edit View Options HelpQuwnU%2 在剖面底部点击鼠标 左键再放开，从下往上 移动鼠标，选定土层范 围，然后再点击鼠标 左 键，修改土层编号HaBsrwxikReady昌 QL:; c E扇聚m用目 t t - . I图18

40、再就是确定根系随深度的分布，根系分布函数b(x)是一个很特殊的函数，它满足以下条件£心皿二1=工彳占其中LR为根系层厚度，x是深度；Az是节点间距，人是每个间距中的根系分布函数值，M是根系 层占节点数。本算例中根系层的厚度为1 m,假设根系分布函数为线性，并有b(x) =2/9 Ml容易证明上式满足积分为1的条件。设置方法为，选择菜单Conditions/ Root Distribution,在在下拉工 具条中使用Edit condition,划定根系层范围(100 cm),把顶部数值设置为0,底部数值设置为2,让程序自动 进行线性差值。ConditionsRoot Distribu

41、tion使用工具条V13 Hydrus-ID - Profile Informnt ionBopIs lineu disCondition Specif. |X|限系分布函 故1 -r«pvd-ucVse top valne for both图19卜一步处理初始条件，选择菜单Conditions/ Initial Conditiono本算例中初始地卜水位高于底板 100 cm,假设土壤剖面初始状态是静力平衡态,则模型底部的压力水头为+ 100 cm,地面的压力水头为- 200 cme使用Edit Condition!具条分别设置顶部和底部的压力水头,并让程序自动插值形成初始条件。最后

42、,选择菜单Conditions/ Observation Points,使用Insert工具条添加若干观察点,有必要在靠近模 型底部的位置加一个观察点以判断地下水位的变化。关闭图形程序退回主程序。17 .编辑土壤剖面使用表格在前处理窗11双击Soil Profile- Summary工具条，程序将弹出一个表格（图20 ）。在这 个表格中可以进一步修改土壤剖面数据。图20xkaiak鲁淳|«?|«)出*rm>> W*，*必1 11d1IX WDIM, »'m w1 -awoaous.<00xzs00 MKK/I工 318 .运行模型我们已经

43、把所有需要的数据都输入了模型中，下IEI就可以运行模拟计算程序了。选择菜W-Calculation/Execute Hydrus,在弹出的对话框中选择0K,则会进入模拟计算模块。软件将调用HYDRUS-1D 的核心程序进行计算，并输出有关的信息，见图21。I 口卜 Hydrus-ID Calculation : test12?39112255126 B3.S4E>O1- 36R<H Y>,23E.：*29.(1-3222£O5a.(4E«»a-u.54P«1P36 丘- 01 T.23E RS-IBM2As 29S22687.54EH9

44、1B360制 a23E-BS'«NIBM226096i4&«e6e.七22EIB1J812t-65X22621e.i3Etee-9.,54E»«19 - B6F - 01 F - 23 f：发能109671*22623a.iZgg A.£4E 炮 Xha29.M.7322625a.i2E*©»-«.54?«i036&wie,23Efl5".',A-5273.IBMTir>frItU1 tCw- ulop9vTop29.707422627a.i2E*®

45、;a-u.54?«1B36E*023EBSW，«w-V*-S2 48.IBM>222639.54EHMO.S6E«fil-e.23E-BS*«IV1*-5242 , e36E、eiivnK7264ftO.nR<M4J.R4P0i气23EV，"sS2?$lec!；2264?6.19E«a8-0.QE由IM2t.7fc142±64S0.BE*i>a-O.S3 仙1IM29.78S332649a.HE»(ia-U.53P4MB36E*01<l23E>£.«/W"

46、52fi7IBMZJ.916B226SCO.ieE«tM H.53EHI1B36E«fila23E-BS'«NIMM*-Siy7. 0 ?6E -IBM2?第26?2$9.nE«e8>e.*七 22ESrtO*IS3263a nDtiA/_A, 53询1092今七；6 22662e.iiE«(M-a036&<ei23Ea5”,s51ELleaIBM4268a.HE*BHL.53E*«1D.36E*ai-a.23E-m»-MlAA-5iiZReal tSe lAoc J O.TAV/i Pre

47、71;oEntAr "comimjo图21一般情况下，如果迭代计算的收敛性较好，模拟时期又不是很长，将很快计算完毕。但是, 一旦出现难以收敛的情况，则可能等待很长时间才有结果，或者非正常中断。本算例模型在很短的时间内就能够完成计算。19 .察看结果梢±1'；云筲寺比夕后 梢出结里于口；函;卜理由TT由的幺小丁目名宓裒看E国99HHYDRTIS-1D | ± Qst LcObservation NodesKAU4C.tUVerticil二HJr?XObservation Nodes: Water Content后处理 观察点结果剖而 曲线图像边界水 分通童水

48、力特性曲线运行迭代信息物质平衡信息Post-processinqObservation PointsProfile Information迎旦f em. .Flo 兄二.Bound 熬坯Fluxes.Wlcul at ion Results Options Window Help(*£*»!< | Eri>t I1 !> 1 nt BCm Koot Water Uptake Models 卷Root Water Uptake - Water Stre - F® VariableBoundary Conditions< -JUL |>F

49、or Help, press Fl图22模型模拟的结果包才舌土壤含水量、压力水头、实际根系吸水（蒸腾）量、实际地表通量、 土壤水分存贮总量等的变化。20 .输出结果如果我们想把模拟结果保存为其它文件格式，用其它软件进行后处理，则可以做进一步的转换。HYDRUS-1D的结果图包含曲线信息，点击右键,会看到有一个Edit Chart Data菜单，用这个菜单 可以通过电子表格输出模拟结果，见图23。这个电子表格的个别数据可以修改后重新显示图像。图 像可以打印。Boundary Water Fluxes and 一KcrtzoAtciLVert ic id|TimaActual Root Water

50、 UptakePrevious RxtDcfeAjl t PrintClose Data Grid Edi tor图23然而，图23中的电子表格不具备输出数据的功能。为此，你可以打开模型目录(文件夹),找到以OUT为后缀的文件，图24。这些文件都可以用记事本打开察看数据。B test向冈文件(E)编群空看(Y)，存放模型的目录收箴®工具帮助(e文件夹1 O C:aToolshydrlDProjuctsVtust文件夹HJ田名称AATNOSPH. IN2 印/EL - OUT BALANCE. OUT I.CHECK. OUT N3DJNF. OUT OBS NODE. OUT PRO

51、FILE. OUT RUN INF. d it-level OUT：EJ转到大 水类型4 KB IM文件3 KB IM文件4 KB OUT文件20 KB OUT 文件PRCi?ILE. DAT£ DESCRIPT. TKT可以用记事本打开的结果文件XXV33 49190 136)UTOUT7SVBKBKBKBKBKBII选定8个对我L45 iW £犹的电筋图24有人量的计算结果保存在T.LEVEL OUT文件中，查阅HYDRUS-1D的用户手册，就知道文件保存了那些结果变量，例如根系吸水强度vROOT等等。附表13.1给出了变量缩写名称的含义。Table 13.1. T_LEVELOUT pressure heads and fluxes on the boundanes and in the root zone.TPne Time. I. at cuiTent Tune-level TJ.rTop Potential surface flux LT1 (infilttation evaporation: /+).rRootPote