HYDRUS1D在土壤污染物迁移模拟中的应用与实践

1. 项目背景与核心价值

在环境影响评价工作中，土壤污染物垂直入渗模拟是评估污染物迁移扩散风险的关键环节。传统方法多采用经验公式或简化模型，难以准确反映复杂地质条件下的污染物运移规律。HYDRUS1D作为国际公认的土壤水动力和溶质运移模拟软件，其基于Richards方程和对流-弥散方程的数值解法，能够精确模拟污染物在非饱和带中的一维垂向迁移过程。

我在多个石化、矿山类环评项目中实践发现，采用HYDRUS1D进行预测具有三大优势：一是能够考虑土壤分层异质性，通过设置不同土层的hydraulic parameters（如van Genuchten参数）实现精准建模；二是支持多种边界条件设置，可模拟降雨入渗、蒸发、抽水等实际工况；三是输出结果可视化程度高，可直接生成浓度时空分布图用于报告编制。

2. 模型构建全流程解析

2.1 基础数据准备要点

完整的模拟需要四类基础数据：

1. 场地地质数据：包括钻孔柱状图、土壤分层厚度、容重、孔隙度等物理参数。建议至少收集3个代表性钻孔数据，通过加权平均确定各土层参数。
2. 污染物特性数据：重点关注分配系数(Kd)、降解半衰期、分子扩散系数等。对于重金属需注意氧化还原状态对吸附行为的影响。
3. 气象水文数据：至少需要10年日降雨量序列（用于确定上边界条件），建议从当地气象站获取并校验数据完整性。
4. 初始条件：通过现场采样确定污染物在表层土壤中的初始浓度分布，采样深度建议达到预测深度的1.5倍。

关键提示：Kd值对模拟结果影响显著，实验室测得的数据往往与现场条件存在差异，建议通过敏感性分析确定合理取值区间。

2.2 参数化建模实操步骤

1. 软件界面导航：

   • 安装后首次运行需设置单位制（建议采用cm-day单位系）
   • 主界面依次为"Project"-"Geometry"-"Flow"-"Transport"-"Results"

2. 几何模型建立：

   python复制# 示例：定义1.5m厚度的三层土壤结构
   Geometry Editor:
   - Layer 1 (0-30cm): 砂质壤土 
   - Layer 2 (30-100cm): 黏土
   - Layer 3 (100-1