土壤水分特征曲线（参数率定应用） · Soil Water Retention Curve — Parameter Calibration

01

模型定义

描述土壤含水率与基质势（吸力）关系的曲线，常用 van Genuchten 模型。

02

公式

van Genuchten：θ(h) = θ_r + (θ_s − θ_r)/[1 + |αh|^n]^m m = 1 − 1/n

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

03

输入与输出

项	内容
输入	参数率定场景：参数 α、n、θ_r（残余含水率）、θ_s（饱和含水率）使用多站多年观测序列率定模型参数，输出率定后的参数集与误差统计。
输出	任意吸力下的含水率，或反之

04

参考依据

van Genuchten (1980)；Mualem (1976)

05

成熟度与边界

成熟度

Established

领域

土壤 · Soil

适用边界参数随土壤质地差异巨大，必须由实测（压力膜/离心法）率定；存在滞后效应。参数率定需划分率定与验证期，检验参数外推能力，避免过拟合单场过程。

06

机理与过程链

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

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

07

核心假设与成立条件

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

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

08

误差来源与失效模式

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

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

09

输入数据条件与精度要求

本模型最小可用数据集：参数率定场景：参数 α、n、θ_r（残余含水率）、θ_s（饱和含水率）使用多站多年观测序列率定模型参数，输出率定后的参数集与误差统计。

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

输入

参数率定场景：参数 α、n、θ_r（残余含水率）、θ_s（饱和含水率）使用多站多年观测序列率定模型参数，输出率定后的参数集与误差统计。

输出

任意吸力下的含水率，或反之

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