[发明专利]一种用于测量土壤垂向饱和渗透系数的方法

权利要求书

1.一种用于测量土壤垂向饱和渗透系数的方法，其特征在于该方法包括以下步骤：

(1)将单环(1)缓慢均匀插入土壤，记录单环(1)插入土壤的深度L；

(2)将土壤含水量测量仪(3)插入单环(1)内的土壤中，测量记录土壤初始含水量；

(3)将自记水位计(2)放入单环(1)内，往单环(1)内注水，记录单环(1)水位变化，即时间序列{t_i}及单环内水深变化序列{H_i}；

(4)在单环(1)中水入渗结束时，用土壤含水量测量仪(3)测定单环(1)内土壤含水量；

(5)根据单环(1)变水头入渗试验中水深变化过程的计算原理，将单环(1)内水流入渗分为两个阶段，单环(1)中的一维入渗及出单环(1)后的三维入渗；单环(1)中一维入渗采用一维变水头Green-Ampt入渗模型计算：

t=ΔθK(1-Δθ)2(H0-hΔθ×(1-Δθ)+(H0+C)×ln(H0+CH0+C+(H0-h)×(1-Δθ)/Δθ)),]]>其中H₀-h≤LΔθ(1)

式中：Δθ为土壤缺水量，H₀为单环初始水深，C为土壤吸力，h为单环内水深，t为入渗时间，K为土壤垂向饱和渗透系数，L为单环插入土壤深度；

出单环(1)后的三维入渗包含两部分：重力引起的下渗土壤水以及单环水压力与土壤吸力引起的土壤水扩散，重力入渗水量等于土壤垂向饱和渗透系数乘入渗时间，土壤水扩散采用Nestingen提出的球缺扩散模型计算，此时段单环水位变化计算公式为：

h=H0-K(t-t0)-LΔθ-2R3+3LR2-2r03-3Lr033r12Δθ,]]>其中H₀-h＞LΔθ  (2)

式中：LΔθ为一维入渗水量，K(t-t₀)为重力引起的扩散水量，其中t₀为一维土壤水入渗结束时间，计算公式为：

t0=ΔθK(1-Δθ)2(L×(1-Δθ)+(H0+C)×ln(H0+CH0+C+L×(1-Δθ)))---(3)]]>

为单环水压力与土壤吸力引起的土壤水扩散，其中r₁为单环半径，R为土壤水球缺扩散半径，r₀为初始球缺半径，R可通过求解如下微分方程获得：

dRdt=H0+C-K(t-t0)-LΔθ-2R3+3LR2-2r03-3Lr023r12Δθ(π28r0(L+r0)Lln(RR+Lr0+Lr0)+L)2R(R+L)ΔθKr02R|t=t0=r0---(4)]]>

式中C为土壤吸力，计算公式为：

C=(π28r0(L+r0)Lln(RR+Lr0+Lr0)+L)(|dhKdt|-1)-h---(5)]]>

(6)整理试验数据，根据实测时间序列{t_i}及单环内水深变化序列{H_i}计算出入渗末期单环水位入渗率dH/dt|_end；然后假设土壤垂向饱和渗透系数K，再联立求解下列方程组(6)、(7)和(8)，得出土壤吸力C：

C=(π28r0(L+r0)Lln(RendRend+Lr0+Lr0)+L)(|dHKdt|end-1)-Hend---(6)]]>

Hend=H0-K(tend-t0)-LΔθ-2Rend3+3LRend2-2r03-3Lr023r12Δθ---(7)]]>

t0=ΔθK(1-Δθ)2(L×(1-Δθ)+(H0+C)×ln(H0+CH0+C+L×(1-Δθ)))---(8)]]>

式中：H_end为入渗末期单环水位，t_end为入渗时间，R_end为入渗末期球缺扩散半径；

然后再将土壤垂向饱和渗透系数K、土壤吸力C及实测时间序列{t_i}代入式(1)求出一维入渗时单环内水深；再代入式(4)采用Adams-Bashforth-Moulton变阶法求解该微分方程，获得球缺扩散半径R，然后再将R代入(2)求出三维入渗时单环内水深；此时，求得的单环一维与三维入渗状态下单环内水深系列即为单环变水头入渗模型中实测时间序列{t_i}对应的单环内水深变化序列{h_i}，比较其与实测内水深序列{H_i}的差异调整K重新计算{h_i}，直到{h_i}与{H_i}的误差最小，此时的K即为土壤垂向饱和渗透系数。