格林-安普特入渗模型（实时预报应用） · Green-Ampt Infiltration Model — Real-Time Forecasting

01

模型定义

基于毛管理论的土壤入渗模型，用湿润锋推进描述入渗率随时间衰减的过程，适合短历时强降雨产流入渗分析。

02

公式

f(t) = Ks·(1 + (ψf·Δθ)/F(t)) F 为累计入渗量，f 为入渗率

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

03

输入与输出

项	内容
输入	实时预报场景：饱和导水率 Ks、湿润锋吸力 ψf、初始与饱和含水率差 Δθ 接入实时监测数据流，预报时效与更新频率优先。
输出	入渗率 f(t)、累计入渗量 F(t)

04

参考依据

Green & Ampt (1911)；Mein & Larson (1973) 简化

05

成熟度与边界

成熟度

Established

领域

土壤 · Soil

适用边界假定湿润锋为陡界面；非均质土壤与地表结皮时需修正；参数可用双环入渗试验率定。实时预报需配套数据质量过滤与滚动更新机制，预报结果按时效分级发布。

06

机理与过程链

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

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

07

核心假设与成立条件

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

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

08

误差来源与失效模式

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

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

09

输入数据条件与精度要求

本模型最小可用数据集：实时预报场景：饱和导水率 Ks、湿润锋吸力 ψf、初始与饱和含水率差 Δθ 接入实时监测数据流，预报时效与更新频率优先。

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

输入

实时预报场景：饱和导水率 Ks、湿润锋吸力 ψf、初始与饱和含水率差 Δθ 接入实时监测数据流，预报时效与更新频率优先。

输出

入渗率 f(t)、累计入渗量 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 方程物理模型	过程完整、可预测	参数多，优先流难描述
经验/统计模型	简单、快速	外推弱，极端条件失效
遥感土壤水分	面尺度覆盖	表层为主，精度受植被与地表粗糙度影响