hydrus1d使用说明_Hydrus-1D中文说明.doc

Hydrus-1D中文说明.doc

天山北麓平原区包气带水分运移机理与数值分析HYDRUS1D 软件使用手册目 录一、软件介绍二、安装指南三、使用指南四、实例介绍1一、软件介绍HYDRUS1D 是国际地下水模型中心公布的，计算包气带水分、盐分运移规律的软件，用它可以解算在不同边界条件制约下的数学模型。若将坐标原点选在地面，取 z 轴向下为正，则一维饱和非饱和带水分运移基本方程为 1hKStzz式中 为含水率；h 为负压水头；S 为植物根系吸水量，对裸露区为0。二、安装指南第一步安装界面；2第二步选择安装目录；第三步修改存放名称；3第四步安装进度显示条；4第五步完成安装；三、使用指南(一)初始界面软件开始界面如图所示(二)数据输入及运行51、选择计算模块2、模型几何结构属性土壤材质数量土层分为几层垂直向下轴土壤深度3、时间信息64、选择输出结果5、迭代步长76、土壤水运移方程7、岩性参数88、上、下边界条件确定99、上、下边界条件数据输入10、几何结构确定1011、初始条件12、模型运行11(三)运行结果1、剖面水头、含水率等计算结果122、边界流量信息3、参数属性134、运行时间信息5、水均衡信息14四、实例操作流程介绍1、实例介绍对奎屯原位试验场包气带裸露区非均质水分运移规律进行模拟。剖面岩性结构从上到下共 7 层；包括壤土、粉壤土、粉质黏壤土以及砂质壤土 4 种岩性类型；水位埋深最大不超过 9m；模型运行期为 4 月 1 号至 11 月 20 号，共234 天。选择万方程作为模型计算的方程；上边界为变化的负压边界，下边界为变化的水头边界；初始条件以 4 月 1 号观测的负压为准输入模型进行计算；输出每月最后一天的结果信息。2、软件操作流程第一步软件可以进行计算的模块包括水分运移模块、溶质运移模块、热传导模块、植物根系吸收水分模块以及植物根系生长模块等，本次模拟只对包气带水分运移规律进行模拟，故选择水流模块，操作如下图所示15第二步选定计算长度单位 m，剖面岩性类型 4 种，从上到下共分为 7 层，计算深度为 9m，操作如下图所示第二步选定计算时间单位 Days，计算时间从 4 月 1 号至 11 月 20 号，共234 天，边界条件随时间变化记录的最小间隔为 1 天，故随时间变化的边界条件记录数为 234 天，操作如下图所示16第三步要求输出每月最后一天的计算结果，统计共有 9 天的计算结果，操作如下图所示第四步迭代误差限的选择以模型默认的为主，如果软件在运行的过程中出现计算不收敛的情况，可以通过适当调整误差限来达到计算收敛的目的，操作如下图所示17第五步土壤水力传导模块选择土壤水计算中普遍认可的且通用的 wan 方程作为本次计算的拟合方程，计算中不考虑水分运移滞后现象，操作如下图所示第六步包气带 4 种不同类型土壤岩性参数，可以通过输入实验室土壤颗粒分析结果获得，操作如下图所示18第七步水流运移模块上边界条件据实际情况选择随时间变化的水头边界，下边界选择随时间变化的水头边界，以变化的地下水位确定，初始条件选择压力水头来确定，操作如下图所示第八步输入随时间变化的上、下边界压力水头值，操作如下图所示19第九步剖面不同岩性层位位置的确定，点击下一步即可。第十步根据奎屯原位试验场剖面开挖实际情况，设置包气带不同岩性层，设置结果如下图所示20第十一步包气带剖面按间隔 5cm 进行剖分计算，观测数据以线性插值的方式获得，如下图所示第十二步软件运行的准备工作设置就绪，点击 OK 软件进行运行计算，如下图所示21第十三步下图为软件运行完毕对话框，从下图中可以看出不同时间的计算结果。第十四步不同埋深点包气带剖面压力水头随时间的变化曲线，如下图所22示第十五步不同埋深点包气带剖面含水率随时间的变化曲线，如下图所示23第十六步包气带剖面压力水头变化曲线，如下图所示第十七步包气带剖面含水率变化曲线，如下图所示24第十八步模型上边界流量随时间的变化曲线，如下图所示第十九步剖面含水率与压力水头关系曲线，如下图所示25第二十步模型运行时间和时间步长变化关系曲线，如下图所示第二十一步模型计算结果水均衡信息表，其中包括上、下边界通量的变化值，如下图所示26注意事项1、参数的初值可以通过试验室测定结果来给出，但最终参数的确定需要通过对模型的识别与验证才能确定。2、单一包气带岩性结构且不同埋深的情况可以通过改变下边界条件来实现。3、上边界的给出可以按照大气边界处理，也可以按负压的方式给出。4、初始条件的确定关键在于地表值得确定，常常需要反复计算，才可以得到较为满意的效果。