下载文档原格式
HYDRUS_2D中文程序说明李道西HYDRUS_2D是一个可用来模拟地下滴灌土壤水流及溶质二维运动的有限元计算机模型。
该模型的水流状态为二维或轴对称三维等温饱和－非饱和达西水流，忽略空气对土壤水流运动的影响，水流控制方程采用修改过的Richards方程，即嵌入汇源项以考虑作物根系吸水。
程序可以灵活处理各类水流边界，包括定水头和变水头边界、给定流量边界、渗水边界、自由排水边界、大气边界以及排水沟等。
水流区域本身可以是不规则水流边界，甚至还可以由各向异性的非均质土壤组成。
通过对水流区域进行不规则三角形网格剖分，控制方程采用伽辽金线状有限元法进行求解。
无论饱和或非饱和条件，对时间的离散均采用隐式差分。
采用迭代法将离散化后的非线性控制方程组线性化。
对于非饱和土壤水力特性，HYDRUS_2D采用VG模型进行描述，嵌入了Scott [1983] 、Kool 和 Parker [1987]经验模型中的假定：吸湿（脱湿）扫描线与主吸湿（脱湿）曲线成比例变化。
并运用一个比例程序，将用户定义的水力传导曲线与参考土壤相比较，通过线性比例变换，在给定的土壤剖面近似水力传导变量。
第一节程序的基本模块HYDRUS_2D程序模块可以顺序嵌套调用，由以下七个基本模块组成：HYDRUS2D：主程序，定义系统的整个计算机环境。
它控制整个程序的运行过程，根据需要调用相应的子程序模块。
程序执行前，首先需选定模拟选项，包括水流、溶质运移、热运移或是否考虑根系吸水等；然后给定时空单位、土壤水力参数以及用来模拟的边界条件。
程序执行后，可输出一系列土壤水力特性曲线、设定观测点处随时间变化的含水率或负压水头曲线，以及沿边界的实际或累积水通量。
输出文件还可提供质量平衡信息和逆向最优结果。
. Project Manager：该模块用来管理已建立的工程数据，包括打开、删除、重命名工程和保存工程的输入输出数据等。
每个工程可能是针对不同的具体问题，Project Manager会自动将每个工程单独建立一个以工程名命名的文件夹保存相应的工程数据。
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_2D中文程序说明李道西HYDRUS_2D是一个可用来模拟地下滴灌土壤水流及溶质二维运动的有限元计算机模型。
该模型的水流状态为二维或轴对称三维等温饱和－非饱和达西水流，忽略空气对土壤水流运动的影响，水流控制方程采用修改过的Richards方程，即嵌入汇源项以考虑作物根系吸水。
程序可以灵活处理各类水流边界，包括定水头和变水头边界、给定流量边界、渗水边界、自由排水边界、大气边界以及排水沟等。
水流区域本身可以是不规则水流边界，甚至还可以由各向异性的非均质土壤组成。
通过对水流区域进行不规则三角形网格剖分，控制方程采用伽辽金线状有限元法进行求解。
无论饱和或非饱和条件，对时间的离散均采用隐式差分。
采用迭代法将离散化后的非线性控制方程组线性化。
对于非饱和土壤水力特性，HYDRUS_2D采用VG模型进行描述，嵌入了Scott [1983] 、Kool 和 Parker [1987]经验模型中的假定：吸湿（脱湿）扫描线与主吸湿（脱湿）曲线成比例变化。
并运用一个比例程序，将用户定义的水力传导曲线与参考土壤相比较，通过线性比例变换，在给定的土壤剖面近似水力传导变量。
第一节程序的基本模块HYDRUS_2D程序模块可以顺序嵌套调用，由以下七个基本模块组成：HYDRUS2D：主程序，定义系统的整个计算机环境。
它控制整个程序的运行过程，根据需要调用相应的子程序模块。
程序执行前，首先需选定模拟选项，包括水流、溶质运移、热运移或是否考虑根系吸水等；然后给定时空单位、土壤水力参数以及用来模拟的边界条件。
程序执行后，可输出一系列土壤水力特性曲线、设定观测点处随时间变化的含水率或负压水头曲线，以及沿边界的实际或累积水通量。
输出文件还可提供质量平衡信息和逆向最优结果。
. Project Manager：该模块用来管理已建立的工程数据，包括打开、删除、重命名工程和保存工程的输入输出数据等。
每个工程可能是针对不同的具体问题，Project Manager会自动将每个工程单独建立一个以工程名命名的文件夹保存相应的工程数据。
天山北麓平原区包气带水分运移机理与数值分析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种不同类型土壤岩性参数，可以通过输入实验室土壤颗粒分析结果获得，操作如下图所示：第七步：水流运移模块上边界条件据实际情况选择随时间变化的水头边界，下边界选择随时间变化的水头边界，以变化的地下水位确定，初始条件选择压力水头来确定，操作如下图所示：第八步：输入随时间变化的上、下边界压力水头值，操作如下图所示：第九步：剖面不同岩性层位位置的确定，点击下一步即可。
天山北麓平原区包气带水分运移机理与数值分析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种不同类型土壤岩性参数，可以通过输入实验室土壤颗粒分析结果获得，操作如下图所示：第七步：水流运移模块上边界条件据实际情况选择随时间变化的水头边界，下边界选择随时间变化的水头边界，以变化的地下水位确定，初始条件选择压力水头来确定，操作如下图所示：第八步：输入随时间变化的上、下边界压力水头值，操作如下图所示：第九步：剖面不同岩性层位位置的确定，点击下一步即可。
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图形用户界面介绍)过去几十年，模拟在地下水流和污染物迁移的复杂数值模型暴增，包括处理在土壤表面与地下水位之间非饱和或渗流区一和多维水流和迁移过程模型。
用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-1D模拟剖面变饱和度地下水流（简明手册）目录1. 如何获取HYDRUS-1D (2)2. 版权声明 (2)3. 参考资料 (2)4. HYDRUS-1D的WINDOWS界面 (2)5. 设计模型 (3)6. 使用HYDRUS-1D创建模型 (4)7. 输入模型控制信息 (5)8. 水流模型——迭代计算参数 (6)9. 水流模型——土壤水力特性模型 (7)10. 水流模型——土壤水分特征曲线 (7)11. 水流模型——边界条件 (8)12. 水流模型——定水头或通量边界设置 (9)13. 根系吸水——吸水模型 (9)14. 根系吸水——水分胁迫参数 (10)15. 输入可变边界条件的信息 (10)16. 编辑土壤剖面——使用图形界面 (11)17. 编辑土壤剖面——使用表格 (13)18. 运行模型 (14)19. 察看结果 (14)20. 输出结果 (14)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水流和溶质运移模拟软件介绍HYDRUS是一个运行于Windows系统下的环境模拟软件，主要用于变量饱和多孔介质的水流和溶质运移。
HYDRUS包括用于模拟变量饱和多孔介质下的水、热和多溶质运移的二维和三维有限元计算，包括一个参数优化算法，用于各种土壤的水压和溶质运移参数的逆向估计。
该模型互动的图形界面，可进行数据前处理、结构化和非结构化的有限元网格生成以及结果的图形展示。
HYDRUS一共五个版本，用户可以选择最适合自己版本。
用户可以选择局限于一般功能的二维应用（2D-Standard版本，与之前含有MeshGen-2D的Hydrus-2D功能一致）或者二维和三维应用（如3D-Standard 或3D-Professional）。
用户也可以选择相对简单的（二维直角几何图形—3D-Lite, 与之前不含MeshGen-2D的Hydrus-2D功能一致）或三维的几何立体图形– 3D-Lite）或更复杂的几何图形（用于普通二维几何图形的2D-Standard或在二维基础上以及分层三维的3D-Standard，以及用于普通三维几何图形的3D-Professional）。
用户也可以选择从低版本升级到高版本。
标准计算模型HYDRUS是模拟变量饱和多孔介质下的水、热和多溶质二维和三维运动的有限元计算模型。
HYDRUS数值求解饱和非饱和水流的Richards方程和热传递和溶质运移的对流扩散型方程。
水流方程包含一个下沉期，可导致植物根系吸水。
热传递方程考虑了水流传导和对流运动。
对流扩散的溶质运移方程的管理是一个非常普遍的形式，包括固体和液态非线性非平衡反应的规定以及液体和气体的线性平衡反应。
因此，不管是吸附溶质还是挥发溶质（如杀虫剂）都已经考虑到了。
溶质运移方程还包括了零阶生产的影响、其他溶质的独立一级降解以及一阶衰减和生产反应，以便提供连续一级链中溶质间所需的耦合。
运移模拟也会引起液相对流和扩散、气相扩散，因此次模型在液态和气态条件下可同时模拟溶质运移。
天山北麓平原区包气带水分运移机理与数值分析HYDRU S—1D软件使用手册
目录
一、软件介绍
二、安装指南
三、使用指南
四、实例介绍
一、软件介绍
HYDRUS—1D是国际地下水模型中心公布的，计算包气带水分、盐分运移规律的软件，用它可以解算在不同边界条件制约下的数学模型。
若将坐标原点选在地面，取z轴向下为正，则一维饱和—非饱和带水分运移基本方程为：式中：θ为含水率；h为负压水头；S 为植物根系吸水量，对裸露区为0。
二、安装指南
第一步：安装界面；
第二步：选择安装目录；
第三步：修改存放名称；
第四步：安装进度显示条；
第五步：完成安装；
三、使用指南
（一）初始界面
软件开始界面如图所示：
（二）数据输入及运行1、选择计算模块
2、模型几何结构属性
3、时间信息
4、选择输出结果
5、迭代步长
6、土壤水运移方程
7、岩性参数
8、上、下边界条件确定
9、上、下边界条件数据输入
10、几何结构确定
11、初始条件
12、模型运行
（三）运行结果
1、剖面水头、含水率等计算结果
2、边界流量信息
3、参数属性
4、运行时间信息
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种不同类型土壤岩性参数，可以通过输入实验室土壤颗粒分析结果获得，操作如下图所示：
第七步：水流运移模块上边界条件据实际情况选择随时间变化的水头边界，下边界选择随时间变化的水头边界，以变化的地下水位确定，初始条件选择压力水头来确定，操作如下图所示：
第八步：输入随时间变化的上、下边界压力水头值，操作如下图所示：
第九步：剖面不同岩性层位位置的确定，点击下一步即可。
第十步：根据奎屯原位试验场剖面开挖实际情况，设置包气带不同岩性层，设置结果如下图所示：
第十一步：包气带剖面按间隔5cm进行剖分计算，观测数据以线性插值的方式获得，如下图所示：
第十二步：软件运行的准备工作设置就绪，点击OK软件进行运行计算，如下图所示：
第十三步：下图为软件运行完毕对话框，从下图中可以看出不同时间的计算结果。
第十四步：不同埋深点包气带剖面压力水头随时间的变化曲线，如下图所示：第十五步：不同埋深点包气带剖面含水率随时间的变化曲线，如下图所示：第十六步：包气带剖面压力水头变化曲线，如下图所示：
第十七步：包气带剖面含水率变化曲线，如下图所示：
第十八步：模型上边界流量随时间的变化曲线，如下图所示：
第十九步：剖面含水率与压力水头关系曲线，如下图所示：
第二十步：模型运行时间和时间步长变化关系曲线，如下图所示：
第二十一步：模型计算结果水均衡信息表，其中包括上、下边界通量的变化值，如下图所示：
注意事项：
1、参数的初值可以通过试验室测定结果来给出，但最终参数的确定需要通过对模型的识别与验证才能确定。
2、单一包气带岩性结构且不同埋深的情况可以通过改变下边界条件来实现。
3、上边界的给出可以按照大气边界处理，也可以按负压的方式给出。
4、初始条件的确定关键在于地表值得确定，常常需要反复计算，才可以得到较为满意的效果。
