Hydrus D简明使用手册

合集下载

Hydrus D简明使用手册

Hydrus D简明使用手册

用HYDRUS-1D模拟剖面变饱和度地下水流(简明手册)王旭升中国地质大学(北京)目录1. 如何获取HYDRUS-1D (2)2. 版权声明 (2)3. 参考资料 (3)4. HYDRUS-1D的WINDOWS界面 (4)5. 设计模型 (6)6. 使用HYDRUS-1D创建模型 (7)7. 输入模型控制信息 (7)8. 水流模型——迭代计算参数 (10)9. 水流模型——土壤水力特性模型 (11)10. 水流模型——土壤水分特征曲线 (11)11. 水流模型——边界条件 (12)12. 水流模型——定水头或通量边界设置 (13)13. 根系吸水——吸水模型 (14)14. 根系吸水——水分胁迫参数 (15)15. 输入可变边界条件的信息 (16)16. 编辑土壤剖面——使用图形界面 (18)17. 编辑土壤剖面——使用表格 (21)18. 运行模型 (22)19. 察看结果 (22)20. 输出结果 (23)HYDRUS-1D是一个共享专业软件,用于模拟一维变饱和度地下水流、根系吸水、溶质运移和热运移。

本手册只介绍应用HYDRUS1D模拟垂向剖面水流和根系吸水的操作方法。

1. 如何获取HYDRUS-1DHYDRUS-1D由位于欧盟捷克的PC-Progress工程软件开发公司发行,用户可以登录该公司首页: 。

为了下载HYDRUS-1D,应先注册成为用户,然后下载Hydrus-1D的安装文件:H1D_4_14.exe。

这个文件对应目前HYDRUS-1D的最高版本。

2. 版权声明HYDRUS-1D的作者为:(1) J. Simunek, Department of Environmental Sciences, University of California Riverside, Riverside, California, USA.(2) M. Sejna, PC Progress, Prague, Czech Republic.(3) M.Th. van Genuchten, Department of Mechanical Engineering, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil.感谢他们提供了一个如此精美而又免费使用的专业软件,帮助我们从事有关的科学和教育工作。

Hydrus-2D中文程序说明

Hydrus-2D中文程序说明

HYDRUS_2D中文程序说明李道西HYDRUS_2D是一个可用来模拟地下滴灌土壤水流及溶质二维运动的有限元计算机模型。

该模型的水流状态为二维或轴对称三维等温饱和-非饱和达西水流,忽略空气对土壤水流运动的影响,水流控制方程采用修改过的Richards方程,即嵌入汇源项以考虑作物根系吸水。

程序可以灵活处理各类水流边界,包括定水头和变水头边界、给定流量边界、渗水边界、自由排水边界、大气边界以及排水沟等。

水流区域本身可以是不规则水流边界,甚至还可以由各向异性的非均质土壤组成。

通过对水流区域进行不规则三角形网格剖分,控制方程采用伽辽金线状有限元法进行求解。

无论饱和或非饱和条件,对时间的离散均采用隐式差分。

采用迭代法将离散化后的非线性控制方程组线性化。

对于非饱和土壤水力特性,HYDRUS_2D采用VG模型进行描述,嵌入了Scott [1983] 、Kool 和 Parker [1987]经验模型中的假定:吸湿(脱湿)扫描线与主吸湿(脱湿)曲线成比例变化。

并运用一个比例程序,将用户定义的水力传导曲线与参考土壤相比较,通过线性比例变换,在给定的土壤剖面近似水力传导变量。

第一节程序的基本模块HYDRUS_2D程序模块可以顺序嵌套调用,由以下七个基本模块组成:HYDRUS2D:主程序,定义系统的整个计算机环境。

它控制整个程序的运行过程,根据需要调用相应的子程序模块。

程序执行前,首先需选定模拟选项,包括水流、溶质运移、热运移或是否考虑根系吸水等;然后给定时空单位、土壤水力参数以及用来模拟的边界条件。

程序执行后,可输出一系列土壤水力特性曲线、设定观测点处随时间变化的含水率或负压水头曲线,以及沿边界的实际或累积水通量。

输出文件还可提供质量平衡信息和逆向最优结果。

. Project Manager:该模块用来管理已建立的工程数据,包括打开、删除、重命名工程和保存工程的输入输出数据等。

每个工程可能是针对不同的具体问题,Project Manager会自动将每个工程单独建立一个以工程名命名的文件夹保存相应的工程数据。

Hydrus软件使用实例

Hydrus软件使用实例

Hydrus软件使用实例
天山北麓平原区包气带水分运移机理与数值分析HYDRU S—1D 软件使用手册
一、软件介绍
二、安装指南
三、使用指南
四、实例介绍
、软件介绍
HYDRU—1D是国际地下水模型中心公布的,计算包气带水分、盐分运移规
律的软件,用它可以解算在不同边界条件制约下的数学模型。

若将坐标原点选在
地面,取z轴向下为正,则一维饱和一非饱和带水分运移基本方程为:
K (― 1) S
t z z
式中:B为含水率;h为负压水头;S为植物根系吸水量,对裸露区为0。

二、安装指南
第一步:安装界面;
第二步:选择安装目录;
第三步:修改存放名称;
时?::峙却叶, —- ?——-―
HTvaf F/arr*h F*LA*r * k
< fc?*
第四步:安装进度显示条;
第五步:完成安装;
lifmi w: li
HYDRLS-ID
5i^fc . *■*?#
CafTi!ng-CBta
u 也闭ki ?°i ■d'rdn E^fyita ^T^drui
1 d. dii
H - K。

HYDRUS—— D中文说明书

HYDRUS—— D中文说明书

边界条件的选取对模拟的结果的影响也非常关键。一共有如下边界:
上边界 Constant pressure head
定水头
Constant flux
定流量
Atmospheric boundary condition with surface layer
大气边界
Atmospheric boundary condition with surface run off
Hydrus 是美国盐土实验室开发的系列软件,用于计算模拟饱和-非饱和带的 水分运动和溶质运移。系列软件分为 1 维、2 维、3 维三种,分别命名维 Hydrus-1D、 Hydrus-2D、Hydrus-3D。
由于非饱和土壤水主要是 1 维垂向运动的形式,因此 Hydrus-1D 的应用非常 广泛;加上学习 Hydrus-1D 的界面之后,再学习 Hydrus-2D、Hydrus-3D,符合 由易变难的学习规律。所以,本课程将重点放在 Hydrus-1D 的水分运动模拟上。对于 溶质运移及 Hydrus-2D、Hydrus-3D,同学们可在此基础上自学。
加入选择的是定水头或定流量等边界,则仅需要输入水头或流量的值。特别 的,隔水边界是一种流量为 0 的边界。
有些边界带有模型参数,如 Deep Drainage 边界。
10、土壤剖面-图像编辑器
剖面-图像编辑器可将模拟土层的离散划分、土壤质地划分、土壤亚区划分、
土壤初始的含水量(负压)等条件进行编辑,并用图像的形式显现出来。
模型分为两大类,单孔模型和多孔模型。最常用的是单孔模型中的 van Genuchten-Mualem 模型。Hysteresis 表示滞后现象,即吸水和脱水并不是的负压-含水量曲 线不一致。一般选择不考虑滞后。 7、土壤水力参数

HYDRUS-3D指导手册

HYDRUS-3D指导手册

HYDRUS (2D/3D)目录摘要HYDRUS图形用户界面介绍1.项目管理和数据管理2.项目几何信息3.水流参数3.1主要过程3.2反向求解3.3时间信息3.4输出信息3.5迭代准则3.6土壤水力模型3.7水流参数3.8神经网络预测3.9在导水率中各向异性3.10溶质迁移3.11溶质迁移参数3.12溶质反应参数3.13溶质迁移参数温度依赖3.14热扩散参数3.15根系吸水模型3.16根系吸水参数3.17根系分布参数3.18时间变量边界条件3.19人工湿地4.迁移域几何4.1边界对象4.1.1点4.1.2线和折线4.1.3弧线和圆4.1.4曲线和样条曲线4.1.5移动、复制、旋转和镜像操作4.1.6附件操作4.2表面4.2.1定义一个二维域的步骤4.2.2几何正确定义几个注意规则4.2.3内部对象4.3打开4.4立体图形4.4.1立方体分为土柱4.4.2立方体分为次层4.4.3在不同厚度矢量上次层不同厚度的各自特征4.4.4定义一个三维层域的步骤4.5厚度4.6辅助对象4.6.1尺寸4.6.2标记4.6.3位图(结构)4.6.4横断面4.6.5网格线4.7对象其它注意事项4.7.1对象编号4.7.2对象间关系4.7.3对象间参考和写一列Index规矩4.8从一文本文件输入几何5.有限元网格5.1有限元网格创建5.2结构的有限元创建5.3非结构的有限元创建5.4有限元网格改良5.5非结构的有限元创建MeshGen2D5.6有限元网格统计5.7有限元网格部分6.域特性、初始与边界条件6.1默认与特性6.2初始条件6.3边界条件6.4域特性7.图形输出7.1结果-图形显示7.1.1显示选项7.1.2编辑Isoband值和彩色光谱7.2结果-其它信息7.2.1转化到ASCII8.图形用户界面部分8.1视图窗口8.1.1屏幕和视图命令8.1.2坐标方格与工作面8.1.3伸缩因子8.1.4渲染模型8.1.5选项和编辑命令8.1.6弹出菜单8.1.7拖拉与放弃8.1.8部分8.2导航栏8.3编辑栏8.4工具栏8.5HYDRUS菜单9.杂项信息9.1程序选项9.2HYDRUS许可证与激活9.2.1请求准则9.2.2安装、移到其它电脑9.3打印选项9.4坐标体系9.5在计算时期DOS窗口9.6可视文件参考文献Introduction to the HYDRUSGraphical User Interface(HYDRUS图形用户界面介绍)过去几十年,模拟在地下水流和污染物迁移的复杂数值模型暴增,包括处理在土壤表面与地下水位之间非饱和或渗流区一和多维水流和迁移过程模型。

Hydrus软件使用实例

Hydrus软件使用实例

天山北麓平原区包气带水分运移机理与数值分析HYDRU S—1D软件使用手册目录一、软件介绍二、安装指南三、使用指南四、实例介绍一、软件介绍HYDRUS —1D 是国际地下水模型中心公布的,计算包气带水分、盐分运移规律的软件,用它可以解算在不同边界条件制约下的数学模型。

若将坐标原点选在地面,取z 轴向下为正,则一维饱和—非饱和带水分运移基本方程为:()(1)h K S t z z θθ∂∂∂⎡⎤=--⎢⎥∂∂∂⎣⎦式中:θ为含水率;h 为负压水头;S 为植物根系吸水量,对裸露区为0。

二、安装指南第一步:安装界面;第二步:选择安装目录;第三步:修改存放名称;第四步:安装进度显示条;第五步:完成安装;三、使用指南(一)初始界面软件开始界面如图所示:(二)数据输入及运行1、选择计算模块2、模型几何结构属性3、时间信息4、选择输出结果5、迭代步长6、土壤水运移方程7、岩性参数8、上、下边界条件确定9、上、下边界条件数据输入10、几何结构确定11、初始条件12、模型运行(三)运行结果1、剖面水头、含水率等计算结果3、参数属性5、水均衡信息四、实例操作流程介绍1、实例介绍对奎屯原位试验场包气带裸露区非均质水分运移规律进行模拟。

剖面岩性结构从上到下共7层;包括:壤土、粉壤土、粉质黏壤土以及砂质壤土4种岩性类型;水位埋深最大不超过9m;模型运行期为4月1号至11月20号,共234天。

选择万方程作为模型计算的方程;上边界为变化的负压边界,下边界为变化的水头边界;初始条件以4月1号观测的负压为准输入模型进行计算;输出每月最后一天的结果信息。

2、软件操作流程:第一步:软件可以进行计算的模块包括水分运移模块、溶质运移模块、热传导模块、植物根系吸收水分模块以及植物根系生长模块等,本次模拟只对包气带水分运移规律进行模拟,故选择水流模块,操作如下图所示:第二步:选定计算长度单位m,剖面岩性类型4种,从上到下共分为7层,计算深度为9m,操作如下图所示:第二步:选定计算时间单位Days,计算时间从4月1号至11月20号,共234天,边界条件随时间变化记录的最小间隔为1天,故随时间变化的边界条件记录数为234天,操作如下图所示:第三步:要求输出每月最后一天的计算结果,统计共有9天的计算结果,操作如下图所示:第四步:迭代误差限的选择以模型默认的为主,如果软件在运行的过程中出现计算不收敛的情况,可以通过适当调整误差限来达到计算收敛的目的,操作如下图所示:第五步:土壤水力传导模块选择土壤水计算中普遍认可的且通用的wan方程作为本次计算的拟合方程,计算中不考虑水分运移滞后现象,操作如下图所示:第六步:包气带4种不同类型土壤岩性参数,可以通过输入实验室土壤颗粒分析结果获得,操作如下图所示:第七步:水流运移模块上边界条件据实际情况选择随时间变化的水头边界,下边界选择随时间变化的水头边界,以变化的地下水位确定,初始条件选择压力水头来确定,操作如下图所示:第八步:输入随时间变化的上、下边界压力水头值,操作如下图所示:第九步:剖面不同岩性层位位置的确定,点击下一步即可。

(完整版)Hydrus-1D中文说明

(完整版)Hydrus-1D中文说明

天山北麓平原区包气带水分运移机理与数值分析HYDRU S—1D软件使用手册目录一、软件介绍二、安装指南三、使用指南四、实例介绍一、软件介绍HYDRUS —1D 是国际地下水模型中心公布的,计算包气带水分、盐分运移规律的软件,用它可以解算在不同边界条件制约下的数学模型。

若将坐标原点选在地面,取z 轴向下为正,则一维饱和—非饱和带水分运移基本方程为:()(1)h K S t z z θθ∂∂∂⎡⎤=--⎢⎥∂∂∂⎣⎦式中:θ为含水率;h 为负压水头;S 为植物根系吸水量,对裸露区为0。

二、安装指南第一步:安装界面;第二步:选择安装目录;第三步:修改存放名称;第四步:安装进度显示条;第五步:完成安装;三、使用指南(一)初始界面软件开始界面如图所示:(二)数据输入及运行1、选择计算模块2、模型几何结构属性3、时间信息4、选择输出结果5、迭代步长6、土壤水运移方程7、岩性参数8、上、下边界条件确定9、上、下边界条件数据输入10、几何结构确定11、初始条件12、模型运行(三)运行结果1、剖面水头、含水率等计算结果3、参数属性5、水均衡信息四、实例操作流程介绍1、实例介绍对奎屯原位试验场包气带裸露区非均质水分运移规律进行模拟。

剖面岩性结构从上到下共7层;包括:壤土、粉壤土、粉质黏壤土以及砂质壤土4种岩性类型;水位埋深最大不超过9m;模型运行期为4月1号至11月20号,共234天。

选择万方程作为模型计算的方程;上边界为变化的负压边界,下边界为变化的水头边界;初始条件以4月1号观测的负压为准输入模型进行计算;输出每月最后一天的结果信息。

2、软件操作流程:第一步:软件可以进行计算的模块包括水分运移模块、溶质运移模块、热传导模块、植物根系吸收水分模块以及植物根系生长模块等,本次模拟只对包气带水分运移规律进行模拟,故选择水流模块,操作如下图所示:第二步:选定计算长度单位m,剖面岩性类型4种,从上到下共分为7层,计算深度为9m,操作如下图所示:第二步:选定计算时间单位Days,计算时间从4月1号至11月20号,共234天,边界条件随时间变化记录的最小间隔为1天,故随时间变化的边界条件记录数为234天,操作如下图所示:第三步:要求输出每月最后一天的计算结果,统计共有9天的计算结果,操作如下图所示:第四步:迭代误差限的选择以模型默认的为主,如果软件在运行的过程中出现计算不收敛的情况,可以通过适当调整误差限来达到计算收敛的目的,操作如下图所示:第五步:土壤水力传导模块选择土壤水计算中普遍认可的且通用的wan方程作为本次计算的拟合方程,计算中不考虑水分运移滞后现象,操作如下图所示:第六步:包气带4种不同类型土壤岩性参数,可以通过输入实验室土壤颗粒分析结果获得,操作如下图所示:第七步:水流运移模块上边界条件据实际情况选择随时间变化的水头边界,下边界选择随时间变化的水头边界,以变化的地下水位确定,初始条件选择压力水头来确定,操作如下图所示:第八步:输入随时间变化的上、下边界压力水头值,操作如下图所示:第九步:剖面不同岩性层位位置的确定,点击下一步即可。

HYDRUSDD学习手册

HYDRUSDD学习手册

H Y D R U S(2D/3D)天津城建大学模拟在变饱和介质中水、热、复合溶质的二和三维运动软件包用户手册版本1.022007.05目录摘要HYDRUS图形用户界面介绍1.项目管理和数据管理2.项目几何信息3.水流参数3.1 主要过程3.2 反向求解3.3 时间信息3.4 输出信息3.5 迭代准则3.6 土壤水力模型3.7 水流参数3.8 神经网络预测3.9 在导水率中各向异性3.10 溶质迁移3.11 溶质迁移参数3.12 溶质反应参数3.13 溶质迁移参数温度依赖3.14 热扩散参数3.15 根系吸水模型3.16 根系吸水参数3.17 根系分布参数3.18 时间变量边界条件3.19 人工湿地4.迁移域几何4.1 边界对象4.1.1 点4.1.2 线和折线4.1.3 弧线和圆4.1.4 曲线和样条曲线4.1.5 移动、复制、旋转和镜像操作4.1.6 附件操作4.2 表面4.2.1 定义一个二维域的步骤4.2.2 几何正确定义几个注意规则4.2.3 内部对象4.3 打开4.4 立体图形4.4.1 立方体分为土柱4.4.2 立方体分为次层4.4.3 在不同厚度矢量上次层不同厚度的各自特征 4.4.4 定义一个三维层域的步骤4.5 厚度4.6 辅助对象4.6.1 尺寸4.6.2 标记4.6.3 位图(结构)4.6.4 横断面4.6.5 网格线4.7 对象其它注意事项4.7.1 对象编号4.7.2 对象间关系4.7.3 对象间参考和写一列Index规矩4.8 从一文本文件输入几何5. 有限元网格5.1 有限元网格创建5.2 结构的有限元创建5.3 非结构的有限元创建5.4 有限元网格改良5.5非结构的有限元创建MeshGen2D5.6 有限元网格统计5.7 有限元网格部分6.域特性、初始与边界条件6.1 默认与特性6.2 初始条件6.3 边界条件6.4 域特性7.图形输出7.1 结果-图形显示7.1.1 显示选项7.1.2 编辑Isoband值和彩色光谱7.2 结果-其它信息7.2.1 转化到ASCII8.图形用户界面部分8.1 视图窗口8.1.1 屏幕和视图命令8.1.2 坐标方格与工作面8.1.3 伸缩因子8.1.4 渲染模型8.1.5 选项和编辑命令8.1.6 弹出菜单8.1.7 拖拉与放弃8.1.8 部分8.2 导航栏8.3 编辑栏8.4 工具栏8.5 HYDRUS菜单9.杂项信息9.1程序选项9.2 HYDRUS许可证与激活9.2.1 请求准则9.2.2 安装、移到其它电脑9.3打印选项9.4 坐标体系9.5 在计算时期DOS窗口9.6 可视文件参考文献Introduction to the HYDRUSGraphical User Interface(HYDRUS图形用户界面介绍)过去几十年,模拟在地下水流和污染物迁移的复杂数值模型暴增,包括处理在土壤表面与地下水位之间非饱和或渗流区一和多维水流和迁移过程模型。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

用HYDRUS-1D模拟剖面变饱和度地下水流(简明手册)王旭升中国地质大学(北京)目录1. 如何获取HYDRUS-1D (2)2. 版权声明 (2)3. 参考资料 (3)4. HYDRUS-1D的WINDOWS界面 (4)5. 设计模型 (6)6. 使用HYDRUS-1D创建模型 (7)7. 输入模型控制信息 (7)8. 水流模型——迭代计算参数 (10)9. 水流模型——土壤水力特性模型 (11)10. 水流模型——土壤水分特征曲线 (11)11. 水流模型——边界条件 (12)12. 水流模型——定水头或通量边界设置 (13)13. 根系吸水——吸水模型 (14)14. 根系吸水——水分胁迫参数 (15)15. 输入可变边界条件的信息 (16)16. 编辑土壤剖面——使用图形界面 (18)17. 编辑土壤剖面——使用表格 (21)18. 运行模型 (22)19. 察看结果 (22)20. 输出结果 (23)HYDRUS-1D是一个共享专业软件,用于模拟一维变饱和度地下水流、根系吸水、溶质运移和热运移。

本手册只介绍应用HYDRUS1D模拟垂向剖面水流和根系吸水的操作方法。

1. 如何获取HYDRUS-1DHYDRUS-1D由位于欧盟捷克的PC-Progress工程软件开发公司发行,用户可以登录该公司首页: 。

为了下载HYDRUS-1D,应先注册成为用户,然后下载Hydrus-1D的安装文件:H1D_4_14.exe。

这个文件对应目前HYDRUS-1D的最高版本。

2. 版权声明HYDRUS-1D的作者为:(1) J. Simunek, Department of Environmental Sciences, University of California Riverside, Riverside, California, USA.(2) M. Sejna, PC Progress, Prague, Czech Republic.(3) M.Th. van Genuchten, Department of Mechanical Engineering, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil.感谢他们提供了一个如此精美而又免费使用的专业软件,帮助我们从事有关的科学和教育工作。

当你运行H1D_4_14.exe解压文件后,会在您的电脑中产生一个安装目录,其中包含Setup.exe可执行文件。

运行这个文件即可安装HYDRUS-1D软件。

当您安装HYDRUS-1D时,象安装其它软件一样,会出现一个许可协议,从中可知本共享软件也受到美国法规的保护。

3. 参考资料HYDRUS-1D安装之后,在软件运行目录下有HYDRS-1D Manual.pdf文件。

从这个文件您可以了解到HYDRUS-1D的一些技术细节,如水流、溶质运移、热流的方程、一些处理专门问题的模型、输入输出文件等等。

有一个Examples目录,包含大量的模拟算例可供参考。

用户还可以参考以下文献:?im?nek, J., M. Th. van Genuchten, and M. ?ejna, Development and applications of the HYDRUS and STANMOD software packages, and related codes, Vadose Zone Journal, doi:10.2136/VZJ2007.0077, Special Issue ”Vadose Zon e Modeling”, 7(2), 587-600, 2008.Jacques, D., J. ?im?nek, D. Mallants, and M. 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.?im?nek, J. and M. Th. van Genuchten, Modeling nonequilibrium 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的界面如下:图 1所有的前后处理在界面中一目了然,左边是前处理工具,右边是后处理工具。

其中前处理的各项功能如下图所示。

模拟内容选项几何形状参数及剖面方式时间信息输出方式水流——迭代求解控制参数图 25. 设计模型在使用HYDRUS-1D之前,您需要对饱和-非饱和水流模拟的基本原理有所了解,并设计出自己想做的模型,准备好数据。

一个剖面水流模型通常包含以下几个要素:(1)土壤剖面从地面算起的深度,准备模拟那个时间段的水分变化。

(2)土壤分几层,每层土壤的渗透性参数和水分特征曲线是怎样的。

(3)根系是怎么分布的。

(4)是否已经确定地面降雨入渗、蒸发蒸腾的信息,特别是它们随时间的变化。

(5)是否已经确定剖面底部的状态属于哪种类型的边界条件。

下面是一个参考模型的设计图:图 36. 使用HYDRUS-1D创建模型打开HYDRUS-1D软件,选择”File/new”菜单,新建一个模型。

在name一栏中输入本模型的名称”test”,更改模型存放的目录。

图 4需要注意的是,HYDRUS-1D模型本身在计算机中就表现为一系列的输入输出文件,它们存放在与模型名称一致的目录中。

本例中,软件会自动创建一个名称为”test”的目录,而”C\ATOOLS\HYDR1D\Projects”中除了test目录之外,还有一个test.h1d文件。

这是一个模型项目(project)文件,告诉软件下次到哪里去寻找模型。

模型创建之后,会显示前处理和后处理窗口(图5)。

由于是新模型,还没有任何模拟结果,所以后处理窗口是空白的。

图 57. 输入模型控制信息首先,在前处理窗口双击Main Processes,在弹出的对话框中输入模型的描述: a test model. 然后在Simulate 一栏中选中Root Water Uptake,表示想处理根系吸水问题。

电击OK之后,前处理窗口将增加处理根系吸水的工具条。

图 6下一步,是输入模型的几何信息和土层划分信息。

在前处理窗口双击Geometry Information ,在弹出的对话框中输入如图7所示的数据。

图 7 接下来输入时间信息,在前处理窗口双击Time Information ,会弹出一个对话框(图8)。

图 8这个对话框中提供了一些灵活的选项来处理上边界条件的变化,下面简要加以说明:(1) 蒸腾量的每日周期变化 时间模拟30 d 上边界是随时间变自动处理蒸腾量在土壤剖面HYDRUS-1D 可以使用一个经验公式来处理每天24小时潜在蒸腾量的变化,设某天的潜在蒸腾量为T(例如用Pemman公式获取的, cm/d),则p其中T p(t)是瞬时潜在蒸腾量,t为时间。

模型假设早上6点之前以及晚上18点-24点的蒸腾量总和只占全天蒸腾量的1%。

注意本例中蒸腾量的单位是 cm/d。

(2) 降水量的周期变化如果在你的模型中降水量是周期性变化的,HYDRUS-1D也可以用一个公式来处理其中P是周期?t内的平均降雨量。

(3) 使用气象数据也可以在HYDRUS-1D中输入气象数据,它将自动利用这些数据计算潜在蒸散量ETp。

可以选择FAO组织推荐的Penman-Monteith公式,也可以选择Hargreaves公式。

这些公式需要辐射、气温、湿度之类的气象数据。

模型的另一个控制信息是对模拟结果的输出如何进行设置。

在前处理窗口双击Print Information工具条,弹出一个对话框。

本例中确定输出30组模拟结果,每天输出1组。

可以确定需要输出30图 98. 水流模型——迭代计算参数HYDRUS-1D是采用迭代法来处理非线性Richards方程的。

在前处理窗口双击Water Flow- Iteration Criteria工具条,弹出一个设置迭代参数的对话框(图10)。

迭代控制参数的设置具有高度的专业技术性,除非特别了解,一般可以使用默认值。

如果模拟结果出现不收敛的情况,需要对最大迭代次数、迭代精度等参数进行调整,但是在缺乏经验的情况下很难操作。

最多迭代次数增大步长迭代次数最小吸力间距用于生成水分图 10HYDRUS-1D 采用自动控制时间步长的方法来处理迭代的收敛性。

对于每个时步,如果迭代次数太多,就缩小时间步长;如果没经过几次迭代就达到收敛精度,则适当增大时间步长。

9. 水流模型——土壤水力特性模型水分特征曲线是非饱和土壤的重要物理性质,HYDRUS-1D提供了几种方法来处理与之有关的参数。

在前处理窗口双击Water Flow- Soil Hydraulic Properties 工具条,弹出一个设置水力特性模型的对话框(图11)。

双重介质模吸湿和疏干图 11在一般情况下,选择单孔介质模型,并选择用van Genuchten-Mualem公式处理土壤的水力特性就可以了。

如果还要模拟溶质运移,可能需要考虑双重介质模型。

双重介质在同一个点有两个孔隙度或两个渗透率,相当于两种介质的混杂。

双重介质模型能够模拟这两种“介质”之间的水分和盐分交换。

10. 水流模型——土壤水分特征曲线在前处理窗口双击Water Flow- Soil Hydraulic Parameters工具条,弹出一个设置水分特征曲线参数的对话框(图12)。

本例中选择van Genuchten-Mualem公式处理水分特征曲线,其中?, n, l均为控制因子。

HYDRUS-1D软件中提供了一组土壤经验参数库,可供用户参考。

本例中两层土壤的参数直接从数据库中调出:第1层对应Sandy loam,第二层对应sand。

图 12在输入参数时,请注意参数的单位。

11. 水流模型——边界条件在前处理窗口双击Water Flow- Boundary Conditions工具条,弹出一个设置边界条件的对话框(图13)。

相关文档
最新文档