van Genuchten 土壤水分特征模型（设计校核应用） · van Genuchten Soil Water Retention Model — Design and Verification

01

模型定义

描述土壤基质势与含水率关系的连续函数模型，配合 Mualem 假设可推求非饱和导水率，是土壤水数值模拟的标配。

02

公式

θ(h) = θr + (θs − θr)·[1 + (α|h|)ⁿ]^(−m)，m = 1 − 1/n

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

03

输入与输出

项	内容
输入	设计校核场景：残余含水率 θr、饱和含水率 θs、曲线参数 α（1/m）、n（-）采用设计工况与历史极值，输出用于方案比选与安全校核。
输出	体积含水率 θ 随基质势 h 的曲线；结合 Mualem 得非饱和导水率 K(θ)

04

参考依据

van Genuchten (1980)；RETC 软件用户手册

05

成熟度与边界

成熟度

Established

领域

土壤 · Soil

适用边界参数常用压力膜法实测拟合；砂土 n 大（3~6），黏土 n 接近 1.1；缺乏实测时可用 ROSETTA 神经网络估计。设计校核关注参数敏感性分析与极值工况覆盖，结果用于设计与验收文档。

06

机理与过程链

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

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

07

核心假设与成立条件

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

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

08

误差来源与失效模式

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

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

09

输入数据条件与精度要求

本模型最小可用数据集：设计校核场景：残余含水率 θr、饱和含水率 θs、曲线参数 α（1/m）、n（-）采用设计工况与历史极值，输出用于方案比选与安全校核。

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

输入

设计校核场景：残余含水率 θr、饱和含水率 θs、曲线参数 α（1/m）、n（-）采用设计工况与历史极值，输出用于方案比选与安全校核。

输出

体积含水率 θ 随基质势 h 的曲线；结合 Mualem 得非饱和导水率 K(θ)

10

部署、标定与运行要点

围绕van Genuchten 土壤水分特征模型（设计校核应用）的运行前检查清单。

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