水文计算软件:霍顿下渗曲线计算器

 一、软件简介

霍顿下渗曲线计算器是一款基于霍顿下渗公式（Horton Infiltration Equation）开发的专业水文计算工具。本软件可帮助水文学者、水利工程师及相关专业学生快速计算土壤下渗能力随时间的变化过程，并直观绘制下渗率曲线与累积下渗量曲线。

适用场景：

- 水文分析与产流计算教学演示

- 城市排水工程设计中的下渗量估算

- 土壤入渗实验数据对比参考

- 水文学课程实验与作业辅助

---

二、理论基础

2.1 霍顿下渗公式

软件采用经典的霍顿经验公式：

下渗率公式：

f(t) = fc + (f₀ - fc) × e^(-k×t)

累积下渗量公式：

F(t) = fc × t + (f₀ - fc)/k × (1 - e^(-k×t))

2.2 参数说明

| 参数符号 | 名称 | 单位 | 物理意义 |

|:---:|:---|:---:|:---|

| f₀ | 初始下渗率 | mm/h | 降雨初期、土壤完全干燥时的最大下渗能力 |

| fc | 稳定下渗率 | mm/h | 土壤饱和后趋于稳定的下渗率 |

| k | 衰减系数 | 1/h | 表征下渗率衰减速度，与土壤质地密切相关 |

| t | 时间 | h | 从降雨开始起算的时间 |

 2.3 重要水文概念提醒

请注意：

本软件计算的是 “充分供水条件下的下渗能力曲线”。在真实降雨过程中，若降雨强度小于土壤下渗能力，实际下渗率将等于降雨强度而非计算值。本软件适用于估算最大可能下渗量，如需进行实际产流计算，请结合降雨过程线使用。

---

三、界面布局说明

软件主界面采用分区布局，自上而下、自左而右依次为：

┌────────────────────────────────────────────────────────────┐

│                   霍顿下渗曲线计算器                        │

├───────────────┬────────────────────────────────────────────┤

│               │                                            │

│   参数设置     │                                            │

│   （左上）     │              曲线图显示区                   │

│               │              （右侧主区域）                  │

│   计算结果     │                                            │

│   （左中）     │                                            │

│               │                                            │

│   公式说明     │                                            │

│   （左下）     │                                            │

├───────────────┴────────────────────────────────────────────┤

│                      使用说明栏（底部）                      │

└────────────────────────────────────────────────────────────┘

| 区域 | 位置 | 功能 |

|:---|:---|:---|

| 参数设置区 | 左侧上部 | 输入四个模型参数，含参数提示说明 |

| 功能按钮区 | 左侧中部 | 执行计算与绘图命令 |

| 计算结果区 | 左侧中下部 | 显示典型时刻的计算数值表格 |

| 公式说明区 | 左侧下部 | 显示霍顿公式的数学表达式 |

| 曲线图显示区 | 右侧主区域 | 展示生成的下渗曲线图 |

| 使用说明栏 | 底部通栏 | 简要操作提示 |

---

四、操作指南

4.1 参数输入

1. 在“参数设置”区域依次输入四个参数值

2. 每个输入框下方均有灰色小字提示该参数的物理含义

3. 软件预设的默认参数为：

   - 初始下渗率 f₀：**50.0 mm/h**

   - 稳定下渗率 fc：**8.0 mm/h**

   - 衰减系数 k：**0.5 1/h**

   - 计算时长：**10.0 h**

4.2 典型时刻数值计算

点击 【计算典型时刻数值】按钮，软件将计算以下典型时刻的数值：

- t = 0.1 h（6 分钟）

- t = 0.5 h（30 分钟）

- t = 1.0 h

- t = 2.0 h

- t = 5.0 h

- t = t_max（用户设定的总时长）

计算结果将以表格形式显示在“计算结果”区域，包含：

- 下渗率 f(t) （单位：mm/h）

- 累积下渗量 F(t) （单位：mm）

4.3 绘制下渗率曲线

点击 【绘制下渗率曲线】 按钮，右侧图形区将显示：

- 蓝色实线：下渗率随时间变化的完整曲线

- 红色虚线：稳定下渗率 fc 的参考线

- 绿色虚线：初始下渗率 f₀ 的参考线

图表标题将同步显示当前使用的参数值，便于截图存档。

4.4 绘制累积下渗量曲线

点击 【绘制累积下渗量曲线】 按钮，右侧图形区将显示：

- 绿色实线：累积下渗量随时间的增长曲线

该曲线直观反映了一定时间内土壤可吸纳的总水量。

4.5 图形交互操作

软件集成了 Matplotlib 标准交互功能：

| 操作 | 功能 |

|:---|:---|

| 鼠标左键拖拽 | 平移视图 |

| 鼠标滚轮 | 缩放视图 |

| 右键点击 | 弹出菜单，可保存图片或调整显示范围 |

| 底部工具栏 | 复位视图、切换缩放模式等 |

五、参数取值参考

5.1 不同土壤类型的典型参数值

| 土壤类型 | f₀ (mm/h) | fc (mm/h) | k (1/h) |

|:---|:---:|:---:|:---:|

| 粗砂土 | 200 ~ 250 | 20 ~ 30 | 3.0 ~ 5.0 |

| 细砂土 | 100 ~ 150 | 10 ~ 15 | 1.5 ~ 3.0 |

| 粉砂壤土 | 60 ~ 100 | 5 ~ 10 | 0.8 ~ 1.5 |

| 粘壤土 | 30 ~ 60 | 3 ~ 6 | 0.3 ~ 0.8 |

| 重粘土 | 10 ~ 30 | 1 ~ 3 | 0.1 ~ 0.3 |

注：以上数值为经验参考范围，实际工程中建议通过双环入渗实验测定。

5.2 参数验证规则

软件内置以下参数校验逻辑：

- `f₀ > fc`：初始下渗率应大于稳定下渗率，若违反将给出警告提示

- `k > 0`：衰减系数必须为正值

- `t_max > 0`：计算时长必须为正值

六、应用示例

示例：粉砂壤土下渗能力分析

问题： 某流域土壤为粉砂壤土，经实测 f₀ = 80 mm/h，fc = 6 mm/h，k = 1.0 1/h。求 2 小时后的下渗率及前 2 小时的累积下渗量。

操作步骤：

1. 输入参数：f₀ = 80，fc = 6，k = 1.0，t_max = 2.0

2. 点击“计算典型时刻数值”

3. 读取 t = 2.0 对应的结果：

   - 下渗率 ≈ 16.0 mm/h

   - 累积量 ≈ 78.0 mm

七、常见问题解答

Q1：为什么计算出的下渗率比实际观测值大？

 A：本软件计算的是“充分供水条件”下的下渗能力。实际降雨中若雨强较小，土壤未必处于饱和状态，实际下渗率将小于计算值。

Q2：参数 k 应该如何测定？

 A：通常通过双环入渗实验获取实测 t-f 数据序列，再利用最小二乘法拟合得到。本软件暂不支持自动拟合功能，需用户手动输入已知 k 值。

Q3：软件输出的累积量单位是 mm，如何换算为 m³？

 A：累积下渗量 F（mm）表示单位面积上的水深。换算为体积水量时，乘以流域面积即可：V(m³) = F(mm) / 1000 × A(m²)。

Q4：能否导出计算数据为 Excel 文件？

 A：当前版本不支持直接导出。可截图保存图形，或复制“计算结果”区域的文本数据粘贴至 Excel。

八、免责声明

本软件为教学与工程辅助工具，计算结果仅供参考。在实际工程设计和水文预报中，建议结合现场实测数据和专业判断综合使用。使用者需自行承担因参数选择不当或模型适用性偏差导致的计算误差风险。