Hydrus1D简明使用手册.docx
《Hydrus1D简明使用手册.docx》由会员分享,可在线阅读,更多相关《Hydrus1D简明使用手册.docx(27页珍藏版)》请在冰豆网上搜索。
Hydrus1D简明使用手册
用HYDRUS-1D模拟剖面变饱和度地下水流(简明手册)
王旭升
中国地质人学(北京)
1.如何获取HYDRUS-1D2
2.明2
3.参考资料2
4.HYDRUS-1D的WINDOWS界面2
5.设计模型3
6.使用HYDRUS-1D创建模型4
7.输入模型控制信息5
&水流模型一迭代计算参数6
9.水流模型——土壤水力特性模型7
10・水流模型一土壤水分特征曲线7
11.水流模型边界条件8
12.水流模型一定水头或通量边界设置9
13.根系吸水吸水模型9
14.根系吸水一水分胁迫参数10
15.输入可变边界条件的信息10
16.编辑土壤剖面一使用图形界面11
17.编辑土壤剖面——使用表格14
18.运行模型14
19.察看结果15
20.输出结果15
HYDRUS-1D是一个共享专业软件,用于模拟一维变饱和度地下水流、根系吸水、溶质运移和热运移。
本手册只介绍应用HYDRUS1D模拟垂向剖面水流和根系吸水的操作方法。
1.如何获取HYDRUS-1D
HYDRUS-1D由位于欧盟捷克的PC-Progress工程软件开发公司发行,用户可以登录该公司首页:
。
为了下载HYDRUS-1D,应先注册成为用户,然后下载Hydrus-1D的安装文件:
H1D_4_14.exeo这个文件对应目前HYDRUS-1D的最高版本。
2.版权声明
HTORUS-1D的作者为:
(1)J.Simunek,DepartmentofEnvironmentaiSciences,UniversityofCaliforniaRiverside,Riverside,California,USA・
(2)M・Sejna,PCProgress,Prague,CzechRepublic・
(3)M・Th.vanGenuchten,DepartofMechanicalEngineering,FederalUniversity
ofRiodeJaneiro,RiodeJaneiro,Brazi1・
感谢他们提供了一个如此精美而又免费使用的专业软件,帮助我们从事有关的科学和教育工作。
当你运行HlD_4_14.exe解压文件后,会在您的电脑中产生一个安装目录,其中包含Setup,exe可执行文件。
运行这个文件即可安装HYDRUS-1D软件。
当您安装HYDRUS-1D时,象安装其它软件一样,会出现一个许可协议,从中可知本共享软件也受到美国法规的保护。
3.参考资料
HTORUS-1D安装之后,在软件运行目录卞有HYDRS-lDManual.pdf文件。
从这个文件您可以了解到HYDRUS-1D的一些技术细节,如水流、溶质运移、热流的方程、一些处理专门问题的模型、输入输出文件等等。
有一个Examples目录,包含大量的模拟算例可供参考。
用户还可以参考以下文献:
•Simunek,J・,Th.vanGenuchten,andSejna,DevelopmentandapplicationsoftheHYDRUSandSTANMODsoftwarepackages,andrelatedcodes,VadoseZoneJournal,doi:
10・2136/VZJ2007・0077,SpecialIssue"VadoseZoneModeling",7
(2),587-600,200&
•Jacques,D・,J・Simunek,D・Mallants,andM・Th.vanGenuchten,Modelingcoupledhydrologicalandchemicalprocesses:
Long-termuraniumtransportfollowingmineralphosphorusfertilization,VadoseZoneJournal>doi:
10.2136/VZJ2007.0084,SpecialIssue”VadoseZoneModelingM,7
(2),698-711,2008.
•Simunek,J・andTh.vanGenuchten,ModelingnonequilibriumflowandtransportwithHYDRUS,VadoseZoneJournal>doi:
10.2136/VZJ2007.0074,SpecialIssue”VadoseZoneModeling",7
(2),782-797,2008.
这些文献都可以从下载。
4.HYDRUS-1D的WINDOWS界面
运行HYDRUS-1D,可以看到一个Windows的界面如卞:
□HYDRUSIDROOTUPTK
Pot:
Help,pressFl
图1
所有的前后处理在界面中一目了然,左边是前处理工具,右边是后处理工具。
其中前处理的各项功能如下图所示。
图2
5.设计模型
在使用HYDRUS-1D之前,您需要对饱和-非饱和水流模拟的基本原理有所了解,并设计出自己想做的模型,准备好数据。
一个剖面水流模型通常包含以下几个要素:
(1)土壤剖面从地面算起的深度,准备模拟那个时间段的水分变化。
(2)土壤分几层,每层土壤的渗透性参数和水分特征曲线是怎样的。
(3)根系是怎么分布的。
(4)是否已经确定地面降雨入渗、蒸发蒸腾的信息,特别是它们随时间的变化。
(5)是否已经确定剖面底部的状态属于哪种类型的边界条件。
下面是一个参考模型的设计图:
1.细砂壤
2.中砂
隔水底板
6.使用HYDRUS-1D创建模型
打开HYDRUS-1D软件,选择”File/new”菜单,新建一个模型。
在name—栏中输入本模型的名称”test",更改模型存放的目录。
图4
需要注意的是,HYDRUS-1D模型本身在计算机中就表现为一系列的输入输出文件,它们存放在与模型名称一致的目录中。
本例中,软件会自动创建一个名称为”test”的目录,而C\ATOOLS\HYDRlD\Projectsn中除了test目录之外,还有一个test.hid文件。
这是一个模型项目(project)文件,告诉软件下次到哪里去寻找模型。
模型创建之后,会显示前处理和后处理窗「1(图5)o由于是新模型,还没有任何模拟结果,
所以后处理窗II是空白的。
7.输入模型控制信息
首先,在前处理窗II双击MainProcesses,在弹出的对话框中输入模型的描述:
atestmodel.然后在Simulate一栏中选中RootWaterUptake表示想处理根系吸水问题。
电击0K之后,前处理窗II将增加处理根系吸水的工具条。
InverseSolution
ng:
|atestnodel
Sinulftta
7WatarFlov
I—VaporFlow
|Sx.urt-XycroLogySoluteTrasport
AGeneralSolute
「IcnChemistryrKFI(PHIEEQC)
—HeatIransrcrt
7鈿oterVptctke:
—RootSrovih
032Iranzport
下一步,是输入模型的几何信息和土层划分信息。
在前处理窗1丨双击GeometryInformation,在弹出的对话框中输入如图7所示的数据。
图7
接卜来输入时间信息,在前处理窗II双击TimeInformation,会弹出一个对话框(图8)。
图8
这个对话框中提供了一些灵活的选项来处理上边界条件的变化,卜•面简要加以说明:
(1)蒸腾量的每日周期变化
HTORUS-1D可以使用一个经验公式来处理每天24小时潜在蒸腾量的变化,设某天的潜在蒸腾量为打(例如用Pemman公式获取的,cm/d),则
其中%&)是瞬时潜在蒸腾量,F为时间。
模型假设早上6点之前以及晚上18点-24点的蒸腾量总和只占全天蒸腾量的1%。
注意本例中蒸腾量的单位是cm/do
(2)降水量的周期变化
如果在你的模型中降水量是周期性变化的,HYDRUS-1D也可以用一个公式来处理
其中戸是周期At内的平均降雨量。
(3)使用气象数据
也口I以在HTORUS-1D中输入气象数据,它将自动利用这些数据计算潜在蒸散量ETp。
可以选择FAO组织推荐的Penman-Monteith公式,也可以选择Hargreaves公式。
这些公式需要辐射、气温.湿度之类的气彖数据。
模型的另一个控制信息是对模拟结果的输出如何进行设置。
在前处理窗II双击PrintInformation工具条,弹出一个对话框。
本例中确定输出30组模拟结果,每天输出1组。
8.水流模型——迭代计算参数
HTORUS-1D是采用迭代法来处理非线性Richards方程的。
在前处理窗II双击WaterFlo^IterationCriteria工具条,弹出一个设置迭代参数的对话框(图10)。
迭代控制参数的设置具有高度的专业技术性,除非特别了解,一般可以使用默认值。
如果模拟结果出现不收敛
的情况,需要对最人迭代次数、迭代精度等参数进行调整,但是在缺乏经验的情况下很难操作。
IterationCriteria
图10
HYDRUS-1D采用自动控制时间步长的方法来处理迭代的收敛性。
对于每个时步,如果迭代次数太多,就缩小时间步长;如果没经过几次迭代就达到收敛精度,则适当增人时间步长。
9.水流模型一土壤水力特性模型
水分特征曲线是非饱和土壤的重要物理性质,HYDRUS-1D提供了几种方法来处理与之有关的参数。
在前处理窗II双击WaterFlowSoilHydraulicProperties工具条,弹出一个设置水力特性模型的对话框(图11)。
;'LoclmpTables
吸湿和疏卡滞后
过程模型
冷3khystere总
C}|yst«ir*K)sgratantioacurw
Qliyst^rasisio.r^touU^acurwowiggdactRity
Cfcttresdsdo.retentioo.curve(o.opwpiiuLBobL»Lird)rTn.tidily也yEglot"
CK:
tidily"讥werjrve
图11
在一般情况下,选择单孔介质模型,并选择用vanGenuchten-Mualem公式处理土壤的水力特性就可以了。
如果还要模拟溶质运移,可能需要考虑双重介质模型。
双重介质在同一个点有两个孔隙度或两个渗透率,相当于两种介质的混杂。
双重介质模型能够模拟这两种“介质”之间的水分和盐分交换。
10.水流模型——土壤水分特征曲线
在前处理窗II双击WaterFlowSoilHydraulicParameters工具条,弹出一个设置水
分特征曲线参数的对话框(图12)o
本例中选择vanGenuchten~Mua1em公式处理水分特征曲线,
(2.28)
//<0
7?
>0
(2.29)
其中Gn,/均为控制因子。
HYDRUS-1D软件中提供了一组土壤经验参数库,可供用户参考。
本例中两层土壤的参数直接从数据库中调出:
第1层对应Sandyloam,第二层对应sand。
图12
在输入参数时,请注意参数的单位。
11.水流模型一一边界条件
在前处理‘窗II双击WaterFlow-BoundaryConditions工具条,弹出一个设置边界条件的对话框(图⑶。
图13
上边界条件有6种类型,下边界条件有8种类型。
边界类型的确定需要考虑实际条件,在本算例中,上边界选择人气边界条件,在降雨量很人时地表可以产生积水。
植被蒸腾量和土壤蒸发量分开处理,HTORUS-1D推荐使用一个经验公式来把潜在蒸散量分割为蒸腾潜力和土壤蒸发潜力:
T=ETA\-£肛3\已丁SCF
PP'7P(2.75)
Ep=二£7^(1—SCF)
其中ETp为潜在蒸散量(口J以使用Penman-Monteith公式处理气象数据得到,cm/d),Tp为潜在蒸腾量(cm/d),Ep为土壤潜在蒸发量(cm/d),LAI是叶面积指数,R为消光系数,取决于太阳角度、植被类型及叶片空间分布特征。
SCF是一个中间参数,即土壤覆盖度(Soilcoverfraction)o在阔叶植被发育的情况下,消光系数的经验值为卍0・5-0・75。
12.水流模型——定水头或通量边界设置
如果边界条件中包含定水头或定通量的边界,则在前处理窗II双击WaterFlowConstantBC工具条,弹出一个设置边界数据的对话框。
本算例模型中,下边界为定流量边界,实际上就是隔水边界,因此直接输入0即可。
13.根系吸水——吸水模型
在前处理窗II双击RootWaterUptake-Models工具条,弹出一个处理根系吸水模型的对话框(图14)。
水分胁迫模型
盐分胁迫模型
根系吸盐模型
补偿吸水域值+
(临界湿润度)
图14
HYDRUS-1D使用水分胁迫和盐分胁迫模型处理根系的吸水。
对于水分胁迫模型,计算公式为
Ta-.|s(以胪)/\:
(力九(2.19)
LLr
其中7;是潜在蒸腾量(cm/d),7;是实际蒸腾量(cm/d),SCO是吸水强度函数(cm/(cm.d),注意用坐标实际表示深度),从力)是水分胁迫函数,力为土壤压力水头(cm),方(£是根系吸水分配(密度)函数,厶为根系层的深度。
水分胁迫函数有2种经验表示方法,即Feddes模型和S-Shape模型。
Feddes模型是一个梯形函数,只需要知道力值。
而S-Shape模型把水分胁迫和叶片气孔的压力水头联系起来,需要知道气孔压力水头的数值加。
土壤的湿润度可以表示为
co=T(l/Tp=Ja(/?
x)b{x)dx(2.22)
LR
但是如果直接用这种方法来计算实际蒸腾量有一定的问题。
植被其实可以调节不同深度的水分胁迫响应特征;某个深度土壤干燥吸不上水,植被可以加人在比较湿润的土层的吸水量,以补偿不足。
这种现彖称为补偿吸水。
为了模拟根系补偿吸水,HYDRUS-1D提供了一种简化的模型,即如果湿润度高于某个临界值(。
>纵),植被根系可以通过补偿机制充分吸水达到潜在蒸腾量。
如呆湿润度低于这个临界值,补偿机制受到抑制,发生整体的水分胁迫,根系吸水总量将低于潜在蒸腾量,并正比于湿润度。
如果不考虑这种补偿吸水机制,可以令型二1。
14.根系吸水一水分胁迫参数
在前处理窗II双击RootWaterUptake-WaterStressReduction工具条,弹出一个处理水分胁迫参数的对话框(图15)。
本算例中直接从数据库中调入Whst的经验值。
图15根系水分胁迫Feddes模型参数
15.输入可变边界条件的信息
在前处理窗II双击VariableBoundaryConditionsT具条,弹出一个处理时间序列数据的对话框(图16)。
本算例中,在步骤(11)中已经把地面处理大气边界,同时又选择使用消光系数法划分植被蒸腾和土面蒸发,因此需要输入每天的降水、潜在蒸散量、叶面积指数等数据。
还有一个需要输入的数据是最小压力水头值,即地面土壤达到最干燥状态时的压力水头。
从理论上讲,当土壤十分干燥时,吸力很人,而液态孔隙水的压强很小,与空气湿度保持平衡关系,因此有
(2.72)
其中九为最小压力水头,忆为空气绝对湿度,RT/Mg为空气的摩尔气体常数。
空气湿度虽然可以通过气彖数据得到,但这里公式需要的是近地面的空气湿度。
一般情况下,取饱和水汽湿度是可取的,因为2cm深度以下土壤空气的湿度往往都是饱和的,只不过随温度发生变化。
因此,可以根据近地面气温的变化来推算地表土壤的空气湿度(饱和水汽湿度),再换算成压力水头。
HYDRUS-1D中需要输入的是最小压力水头的绝对值,缺省值为
hCritA=|^|=106cm=10,m
这个数值只会对土壤蒸发起作用。
HYDRUS-1D建议:
hCritA所对应的土壤含水量应该至少比残余含水量大0.005,在模拟根系吸水的情况下,九还应该低于图15中的P3。
否则(QP3),当根系吸水的临界值压力水头(P3)和地面蒸发的最小压力水头(Q满足时,会导致回流(inflow)现彖,这是不合理的。
除非存在特别干燥的情况,模型一般不需要仔细处理这些问
题。
图16
时间
(d)
降水fit
(cm/d)
ETp
(cm./d)
hCritA
(cm)
:
时间
(d)
降水量
(cm/d)
ETp
(czn/d)
hCritA
(cm)
LAI
1
0
0.5
100000
21
0
0.4
100000
2.1
2
0
0.6
100000
21
0
0.6
100000
0.3
3
0
0.1
100000
2.2
18
0
0.7
100000
0.3
1
0
0.5
100000
2.2
19
0
0.7
100000
0.3
5
0
0.5
100000
2.2
20
0
0.7
100000
0.3
6
0
0.5
100000
2.2
21
0
0.7
100000
0.3
7
0
0.3
100000
2.3
22
0
0.6
100000
0.3
8
0.7
0.1
100000
2.3
23
0
0.5
100000
0.3
9
0
0.1
100000
2.3
24
0.3
0.2
100000
0.3
10
0
0.5
100000
2.3
25
2.3
0.1
100000
0.3
11
0
0.6
100000
2.2
26
1.1
0.1
100000
0.3
12
0
0.1
100000
2.2
0.5
0.1
100000
0.3
13
3
0.1
100000
2.2
28
0
0.4
100000
0.3
11
0
0.5
100000
21
29
0
0.7
100000
0.3
15
0.5
0.2
100000
21
30
0
0.7
100000
0.3
数据可以先在Excel中准备好,如表lo这些数据可以拷贝到图16的电子表格中。
这些数据显示的降水量、蒸散潜力和叶面积指数变化特征如图16右图所示。
在第17口由于庄稼收割,叶面及指数大幅度下降。
16.编辑土壤剖面一一使用图形界面
在前处理窗I1双击SoilProfile-GraphicalEditorT具条,程序将弹出一个处理土壤剖面的软件(图17)。
这个图形软件的使用比较简单,我们需要注意的是在Conditions菜单
下面有很多子菜单,包含处理各种问题的功能菜单。
首先要做的事情,是确定把土壤剖面离散化为多少个节点。
本算例土壤模型深度为3m,我们希望节点间距达到1cm,因此需要301个节点。
选择菜单Conditions/ProfileDiscretization,在卜拉工具条中把Number修改为301。
Conditions
ProfileDiscretizationMaterialDistributionRootDistributionScalingfactorInitialConditionsSubregions
Conditions菜单剖分节点土壤岩性分层根系分布尺度因子初始条件子区分布观察点
图17
接下来,确定土层的分布,本模型有2个土层。
缺省的土层编号为index二1,就是图12中的细砂壤土层。
选择菜单Conditions/MaterialDistribution,在下拉工具条中使用Editcondition,把下部土壤层设置index=2,这个土层编号为2,实际上就是图12中的中砂层。
图18
再就是确定根系随深度的分布,根系分布函数b(x)是一个很特殊的函数,它满足以下条件
£心皿=1=工佔
w=i
其中LR为根系层厚度,x是深度;Az是节点间距,人是每个间距中的根系分布函数值,M是根系层占节点数。
本算例中根系层的厚度为1m,假设根系分布函数为线性,并有
b(x)=2x9Ml
容易证明上式满足积分为1的条件。
设置方法为,选择菜单Conditions/RootDistribution,在在下拉工具条中使用Editcondition,划定根系层范围(100cm),把顶部数值设置为0,底部数值设置为2,让程序自动进行线性差值。
OK
图19
卜一步处理初始条件,选择菜单Conditions/InitialCondition。
本算例中初始地卜zK位高于底板100cm,假设土壤剖面初始状态是静力平衡态,则模型底部的压力水头为+100cm,地面的压力水头为-200cm。
使用EditCondition!
具条,分别设置顶部和底部的压力水头,并让程序自动插值形成初始条件。
最后,选择菜单Conditions/ObservationPoints,使用Insert工具条添加若干观察点,有必要在靠近模型底部的位置加一个观察点以判断地下水位的变化。
关闭图形程序退回主程序。
17.编辑土壤剖面——使用表格
在前处理窗II双击SoilProfile-Summary工具条,程序将弹出一个表格(图20)。
在这个表格中可以进一步修改土壤剖面数据。
18.运行模型
我们已经把所有需要的数据都输入了模型中,下面就可以运行模拟计算程序了。
选择菜'Y1-Calculation/ExecuteHydrus,在弹出的对话框中选择0K,则会进入模拟计算模块。
软件将调用HYDRUS-1D的核心程序进行计算,并输出有关的信息,见图21。
图21
一般情况下,如果迭代计算的收敛性较好,模拟时期又不是很长,将很快计算完毕。
但是,一旦出现难以收敛的情况,则可能等待很长时间才有结果,或者非正常中断。
本算例模型在很短的时间内就能够完成计算。
19.察看结果
模拟运算完毕之后,模拟结果可以通过后处理窗II中的各个工具条来察看,见图22。
[*•£*»!
<|Eri>t|
1€!
•>•1
ntBC
耕KootWaterUptake一Models卷RootWaterUptake-WaterStre・・F®VariableBoundaryConditions
<—JUL—|>
ForHelp,pressFl
图22
模型模拟的结果包扌