非饱和土壤导水率计算要求!
θ:土壤体积含水量;ψ:土壤水势(基质势,负压);K:非饱和导水率
逻辑
土壤水势‑含水量关系就是土壤水分特征曲线 SWCC,非饱和导水率是从水分特征曲线推导出来的,不是简单公式把 θ、ψ 直接代入算 K。
只给一个点的 θ 和 ψ,缺少模型参数,算不出 K;
有完整 SWCC 参数(van‑Genuchten、Brooks‑Corey),给任意 θ 或 ψ,就可以计算对应非饱和导水率。
两个常用模型


现场实测的替代方法
如果你手上只有传感器连续观测:土壤湿度 θ、土壤基质势 ψ,没有实验室测 Ks:
不能直接算出导水率;
可以做瞬时剖面法(instantaneous profile method),利用一段时间内 θ、ψ 动态变化,结合达西定律,反演算法计算原位非饱和导水率,不需要预先 VG 参数。
瞬时剖面法要求多层传感器,不同深度 θ、ψ 时间序列,还要考虑重力势。
简单总结
1. 单点 θ+ψ ⇒ 算不出土壤导水率
2. 多组 θ‑ψ 配对数据 ⇒ 可拟合水分特征曲线,得到相对导水率\(K/Ks)
3. 多组 θ‑ψ + 饱和导水率 Ks ⇒ 得到真实非饱和导水率 K
4. 多层时序 θ、ψ 原位观测,可用瞬时剖面法反演原位 K,不需要预先实验室参数。












