霍顿入渗模型（临近预报应用） · Horton Infiltration Model — Nowcasting

01

模型定义

描述入渗率从初渗率随时间指数衰减到稳定入渗率的经验模型，适合长时间连续降雨的产流模拟。

02

公式

f(t) = fc + (f0 − fc)·e^(−kt)

公式中各量的定义与单位见"输入 / 输出"。

03

输入与输出

项	内容
输入	临近预报场景：初始入渗率 f0、稳定入渗率 fc、衰减常数 k 以最近观测与趋势外推为输入，更新频率分钟级，输出未来 0~6 小时短时预报。
输出	任意时刻入渗率 f(t) 与累计入渗量

04

参考依据

Horton (1933)；《土壤水动力学》

05

成熟度与边界

成熟度

Established

领域

土壤 · Soil

适用边界参数经验性强，需试验率定；地表积水后入渗受控于供水率而非模型值；与 Green-Ampt 相比更便于工程应用。临近预报需与雷达/卫星实况融合校正，逐分钟滚动更新，并按强度分级发布预警。

06

机理与过程链

土壤模型以水热耦合为核心：水分运动按 Richards 方程（非饱和导水率随基质势变化），热量按热传导方程并计入相变潜热；冻融期还需冰-水相变与冻结特征曲线。

工程判别落在"含水率—基质势—导水率"三者的一致性：田间持水量（约 −33 kPa）与凋萎点（−1500 kPa）是两条基准线；含水率突变、滞后环异常提示优先流或传感器问题。

07

核心假设与成立条件

假设不成立时，模型偏差通常最先体现在趋势而非绝对值上。

假设	工程含义与成立条件
均质各向同性	土层按层均一，层内变异系数 <20%；耕作层与压实层须分层
局部平衡	水热交换瞬时平衡；快速入渗与优先流条件下不成立
无优先流	基质流主导；有大孔隙/裂隙流证据时需双孔隙模型
滞后可忽略	吸湿/脱湿曲线差异 <10%；强滞后土需滞后模型

08

误差来源与失效模式

出现下列表现时应先复核假设与输入，而不是调整输出。

失效模式	表现与判据
优先流	大孔隙使入渗前锋超前，模型低估深层入渗 2~5 倍
冻融相变	冻融期含水率读数受冰水相变干扰，需温度修正
土壤异质性	容重与质地空间变异使单点模型难以外推
盐分干扰	高盐土壤电导率影响电容式水分传感器读数 5%~15%

09

输入数据条件与精度要求

本模型最小可用数据集：临近预报场景：初始入渗率 f0、稳定入渗率 fc、衰减常数 k 以最近观测与趋势外推为输入，更新频率分钟级，输出未来 0~6 小时短时预报。

数据项	采样与精度要求
含水率	5~30 min，剖面 3~5 层，精度 ±3% VWC
土壤温度	5~30 min，同剖面，精度 ±0.5 ℃
容重与质地	分层实测，用于标定与饱和含水率换算
降雨/灌溉	5 min，用于入渗响应与水量平衡核对

输入

临近预报场景：初始入渗率 f0、稳定入渗率 fc、衰减常数 k 以最近观测与趋势外推为输入，更新频率分钟级，输出未来 0~6 小时短时预报。

输出

任意时刻入渗率 f(t) 与累计入渗量

10

部署、标定与运行要点

围绕霍顿入渗模型（临近预报应用）的运行前检查清单。

环节	要点
采样频率	含水率与温度 5~30 min，剖面 3~5 层（10/20/40/60 cm）
分层标定	按质地与容重分别标定，每层独立
冻融处理	冻融期加温度修正与相变标志位
事件核对	与降雨/灌溉联动核对入渗响应滞后（0.5~4 h）

11

验证指标与参考阈值

阈值为工程常用参考区间，具体项目按合同与规范取值。

指标	参考阈值
含水率 RMSE	0.03~0.05 m³/m³
土壤温度 RMSE	<1.5 ℃
入渗响应滞后误差	<2 h
田间持水量偏差	<15%（相对实验室测定）

12

常见问答

含水率读数为什么雨后突变又回落？
优先流或传感器接触不良皆可能；先核对相邻层响应是否同步，再看容重与孔隙度是否合理。

冻融期数据能用吗？
需温度修正并标注相变期；冰水混合段读数不可直接当液态含水率使用。

单点土壤数据怎么外推到田块？
按质地/容重分区，每区至少 1 点标定；空间变异大时用多点平均并给出变异系数。

14

与相邻方法的取舍

同一目标量常有多种建模路线，取舍取决于数据条件与精度要求。

方法	适用	局限
Richards 方程物理模型	过程完整、可预测	参数多，优先流难描述
经验/统计模型	简单、快速	外推弱，极端条件失效
遥感土壤水分	面尺度覆盖	表层为主，精度受植被与地表粗糙度影响