水土保持措施下的土壤入渗研究及次暴雨地表产流计算方法(王秀英, 曹文洪)

热度96票 浏览53次 【共0条评论】【我要评论 时间:1999年11月01日 14:46

水土保持措施下的土壤入渗研究及次暴雨地表产流计算方法

(is9a[@0

王秀英曹文洪
a9Xu(FHB1A?v!N r)z0(中国水利水电科学研究院
)水利论文Rd"QsSIon

摘   要从土壤水分下渗能力曲线出发,从理论和实验分析两方面说明水土保持措施下土壤入渗性能的改良作用。考虑到流域实际的下渗能力,建立了将下渗能力曲线、不稳定降雨过程和前期土壤含水量有机联系起来的流域场(时段)暴雨地表径流计算模式。水利论文#@}|X aFD2F+g

关键词下渗曲线;入渗率;超渗产流水利论文PD6aH~uU

1前言

#o efu)yF B1u8J2HD,w0

  因为降水是沿着地表和地下两种主主要原因是由于大暴雨下高度集中的地表径流及其冲刷侵蚀挟带的泥沙而引起的。鉴于水土保持措施可拦截降水,增加土壤入渗,减少坡面地表径流量而增加地下径流量,改变汇流特性,从而削减洪峰,降低径流侵蚀力,起到减水减沙作用,所以,近年来许多专家指出:应充分重视土壤入渗的研究,主张“拦蓄降水、就地入渗”作为治理江河灾害的战略性措施。因此,研究土壤入渗规律,进而推导降雨地表径流计算模式,对于土壤侵蚀和洪水预报,各种水土保持措施,特别是各种措施的配置方案和设计及其效益评价都具有极为重要的意义。水利论文e0Eeis

2土壤水分下渗能力曲线水利论文mgk'H(Wr9OYg

2.1下渗能力曲线

A$P;k5BQ0

  干旱和半干旱地区以超渗产流为主。下渗过程是一个消退过程。从水动力平衡角度分析,在只考虑垂向一维水流的情况下,水流在非饱和土壤中的运动的基本表达式为水利论文`1a"Z#p$b3b9l.J_@nwI

990618e1.gif (1161 bytes)

YU:n(~#P\0

(1)

xb] d-Ti(}I0

式中v为渗透速度;k为渗透系数;φ为水势;z为垂向距离。

0A5Vvr1b7@r'eZ0

根据质量守恒原理有

\E.j0V jfV7PEl0

990618e2.gif (1069 bytes)

,[9c nMU(u1m0

(2)水利论文b!y"T#G3q3^:Q4g4FfA4T

式中θ为包气带含水量,其余符号同前。联解式(1)和式(2)

d N6Hw7@v }0

990618e3.gif (1438 bytes)

w^@E9S9KjJq0

(3)

%d$iGH;SB Ez,d z~X0

式中D为扩散系数,其余符号同前。此即非饱和下渗的微分方程。水利论文-W*Y7O,Z&{y*_uf%?*o*[

  菲利浦在上述方程的基础上,考虑均质土壤,起始含水量均匀分布及充分供水条件下,经过代换和微分展开,取展开式的前两项求得如下近似计算公式

;Pq.g)|f|0

f(t)=(1/2)St-1/2+A

:LU j] L(RD0

(4)水利论文+B?%q4U2m }

式中f为下渗率;S为吸水系数;A相当于稳定下渗率fc水利论文BY Qx;K8Ju![

  若取扩散系数和渗透系数均为常数,通过分离变量解法,解式(3)和式(1)可得水利论文.@0Hc!O'f'T$wv

f(t)=fc+(f0-fc)e-βt

_nmji9`x0

(5)水利论文2tU!]O:i"GS

式中β为常数,下渗曲线的递减参数(随土质而变);f0为初始下渗率(也称最大下渗率);fc为稳定下渗率;t为时间;e为自然对数底数。

i's3qQDYK3A0

3水土保持措施下的土壤入渗能力水利论文Q6v{5t4hwq/a9y5|$J

  总的说来,水土保持措施的水文效应主要表现在涵养水源,保持水土,改良土壤,调节径流,改造局部地区水文循环等。而凡是构成对土壤物理性质影响的措施因子,都对土壤的入渗能力有重要影响。

4O;N4|~+R["m_0

  在土壤物理性质中影响土壤入渗率的主要因素是土壤的孔隙状况,土壤质地,土壤结构的稳定性,地表是否有结皮,土内是否有间层等。一般说来,地表土壤的孔隙度越高,土壤结构的通透性越好,初期入渗率也越高;剖面通体的孔隙度越高土壤的稳渗率也越高。如果土壤表面被压实或有地表结皮产生,地表层将会被导水率极低的结皮盖住,极大地降低土壤入渗率。

l/Je0s.y"[4sR0

6p40.gif (4607 bytes)

~^ P x#D%[z0

图1 火山灰地带不同地被的下渗率曲线
P;\i-w5A%I0(i=400mm/h雨强下)
水利论文t~)V8L+hRKK

3.1水土保持林草措施的土壤入渗改良作用

+o0una|"DH0

  林草措施对土壤性质起到一定的改良作用,并且随着林龄的增加,草场的发育,这种作用逐渐增强,反映在下渗曲线的上移和曲线形状发生变化(见图1)。这主要是因为植物根系的生长发育过程及枯枝落叶层的形成分解过程,就是对土壤物化性质不断改善的过程,其有机质和腐殖质含量、根系活动使土壤孔隙度发达,整个剖面上的渗透能力比较高;同时,促进生物作用来改良土壤。

zg,afU0
林地枯枝落叶层还可减少雨滴的击溅,防止地表在雨滴打击作用下颗粒分散引起的结皮,同时对泥沙具有一定的过滤作用,防止了土壤孔隙的堵塞。

B0~MQSFbkM0  对晋西黄土残壑区实验林草地及对比裸地的实测资料进行分析,得土壤入渗方程各参数(见表1)。从表1可知,不同林分及草种对土壤入渗性能的影响是不同的。水利论文(bI)n+N t"k

1晋西黄土残壑区各实验地类土壤入渗性能

V |1i8N5?0

霍顿公式参数沙棘林地油松林地草地裸地

fc(mm/min)水利论文9o4zEL.r)gi

0.230.210.130.09

f0(mm/min)

.R+t$Kc:{2r i3b0
0.710.580.480.41

β(1/min)水利论文7f yp'xa'w d8L8v

0.05930.05250.05080.0602

3.2水土保持梯田及耕作措施对土壤入渗的影响水利论文S;c6};x3S6S

  梯田对于改变地形、改良土壤有很大的作用。坡地修成水平梯田后,能够很好地拦截天然降水,增加土壤的蓄水供水能力,利于肥力的保持和作物的生长。根据四川盆地中部丘陵区50个梯田土样及23个坡地土样的分析结果:土壤有机质梯田比坡地高145.8%。而有机质的存在增加了土粒间的胶结而促进土壤团粒的形成,改良了土壤结构。

c^*aN%[NR0  以多次洪水的平均误差最小为原则,直接用下渗曲线公式来优选参数,分析子洲径流试验站内两个小流域:(0.18km2)和(0.133km2)的雨、洪资料,得出如图2所示的两小流域下渗曲线。这两个小流域的地貌植被情况是相仿的,但从图2可知,墨矾沟的下渗能力明显高于团山沟,究其原因是由于它们的梯田面积有所不同。墨矾沟的水平梯田的面积占相当的比例,而团山沟的梯田面积很少。

Vg.}l+_#r%E0

6p41.gif (3609 bytes)水利论文/R!h,Pi5U }
水利论文9U K0_lK5DXK+vV9d
2团山沟、墨矾沟下渗曲线

U/A.w9bt&RaK9s`s0

  水土保持耕作措施通过外界的机械力量的作用来改善土壤水分-物理性状。耕作可改善土壤的孔隙状况和通气性,利于土壤中微生物活动和有机物的积累;耕作可破坏雨后土壤表层的板结层,增强土壤的渗透能力。根据中科院西北水土保持研究所19571958年的研究资料表明:耕作能显著地提高土壤透水性和土壤的含水量,一般能提高15%30%,随着耕作深度的增加,土壤的含水量也随之提高。水利论文$s.b.q6nj)Ar9m

4下渗曲线法产流计算公式的推导和应用水利论文"FXc!@m$n

  土壤非饱和水流基本微分方程的求解是把下渗的出现作为地表饱和的结果,就是说是在充分供水条件下导出的下渗能力曲线(霍顿和菲利浦公式在充分供水条件下与实际资料配合较好便是很好的例证)。实际的降雨过程并不一定都大于下渗过程,这就使得供水条件不充足,所以在产流计算中就得考虑流域实际下渗曲线了。流域实际下渗曲线受控于下渗能力曲线、降雨强度(降雨过程)和前期土壤含水量。一定的土壤含水量和降雨过程便有一条与之对应的下渗曲线。这就使得利用下渗能力曲线计算流域产流存在一定的困难。水利论文gM:Vw2p7U

  若通过实验已知土壤在最干燥情况下的最大下渗率f0,下面以霍顿公式为例推导出考虑下渗能力曲线、降雨过程和前期土壤含水量有机结合下的次降雨(990618e6a.gif (844 bytes)>fc)产流公式。

+\_b\L t ~kJ0

4.1990618e6a.gif (844 bytes)j>f0水利论文+x~dQlY#~ B

  若次降雨的前期土壤含水量为Pa,则下渗率相对要低,下渗能力曲线需向左平移Δta(见图3),平移后有水利论文tKf,y1H `"O T*w

6p42.gif (2917 bytes)

0fdx;J0\/f6J|'n0

3990618e6a.gif (844 bytes)j>f0Pa下下渗曲线的变化

/CU I;va b'a0

f1(t)=f(t+Δta)

dP5|W}0

(6)水利论文D'eV]{1j

平移后图中阴影部分的面积在数量上等于前期土壤含水量为Pa,即水利论文 Kev)o8P

990618e7.gif (1226 bytes)水利论文#vlfW&L0^@

(7)水利论文#oDQG7X:M

联立式(5)、式(6)和式(7),得

v/oNj!X*|:Es0

990618e8a.gif (1551 bytes)水利论文 bs4`;in,@9b3P%Y

(8)水利论文)O9yr.c5~2w.]WM

990618e8.gif (1347 bytes)

d%h7T3^]2V/w I0

Δta也可作如下计算:

`6`*X3r3X.\G8@0

  下渗水量部分将形成地下径流,形成地下径流的那部分下渗率即为下渗能力曲线中包含稳定下渗率的项,即为水利论文,d^ zy R3]3r

fg(t)=fc(1-e-βt)水利论文 Yypp;`1y

(9)

5F$TY k-C#C'M0

990618e10.gif (1478 bytes)水利论文 d7Cz|4P mS5u

(10)水利论文a1D z BO&m8V

W在数值上等于在一定f0条件下的雨前土壤缺水量,其值为流域蓄水容量Wm与流域平均土壤含水量之差,所以

T2qJ5vp hb3~0

f0=β(Wm-Pa)

D e.s7uH&d J0

(11)水利论文R4IS4h6e

由式(11)可知,若Pa0,则f0=βWm,达到最大值,即土壤最干燥时的最大初渗率。从图3中可得水利论文K4H$b:`j+rw5x

990618e12.gif (1558 bytes)

(12)

*OK4^v8n0

+HCf%sU!_p"k:_ w0
990618e12b.gif (1165 bytes)

整理上式亦可得式(8)水利论文8J8\.`$f%e0k*X

  次降雨地表产流量为

9YK%|%hI&H!I H0

990618e13a.gif (1284 bytes)

{.aA1q^I0

(13)

ki;VpK8m&Z(u0

990618e13b.gif (1318 bytes)

/e P0V9s&\ vXx0

式中Rl为第l个网格单元的产流量;t(i)为降雨强度的连续函数;tr为降雨历时;P0-tr为次降雨总量。水利论文a1QDq4TB

  时段(Δtj=tj- tj-1)地表产流量为

3Eg.C,w3Mq,A{0

990618e14a.gif (1541 bytes)

@ bc G[1wm/b+_0

(14)

&K(r/Lz~;|(D0

990618e14b.gif (1391 bytes)

P.\ u Ix@+O"k0

  注意若式(8)、式(13)和式(14)中的f0为土壤在最干燥情况下的下渗率,则Pa为次降雨前的土壤含水量;若f0是土壤含水量为Pa0时的初渗率,那么Pa需用Pa-Pa0代替。从式(8)可算得,当Pa=0时,Δta=0,即曲线不作平移。

rl[:ys:p3?0

4.2 f0>990618e6a.gif (844 bytes)j>fc  

+X0d8Bo ~"g/A;Ae0

  在下渗最初阶段,下渗率具有较大的数值,但在降雨初期若雨强很小,小于下渗能力(i<f),很显然f大于i的那一部分下渗能力没能得到满足,需由i>f的部分降雨来补充(见图4,图中ta为产流时刻),下渗曲线需往右移,即水利论文[9q A:FJ

6p43.gif (4526 bytes)

[Y7]&|Cw!a$A0

4   f0>990618e6a.gif (844 bytes)j>fc时下渗曲线的变化

:A AE!W LJ0

f1(t)=f(t-Δta)

.?5N;JG9J{7o-W@s0

(15)

w6W2jm7|Y3m ](le0

4中平移量

T$Y(w jDyGf0

Δta=ta-t0

g\O'E Jb0

(16)水利论文aKW!C1t;a+X+Mj8_

平移后阴影部分1的面积应与阴影部分2的面积相等,即水利论文,l'_i2~"z RO

990618e17.gif (1488 bytes)水利论文6[)yb5WpL-V

(17)水利论文fl"N2f.GK hB

式中P0-ta0ta时刻的降雨量。联立式(5)、式(15)、式(16)和式(17),得水利论文vz nKJ)gyz

990618e18.gif (1338 bytes)

a$c#v?:lLiQn0

(18)水利论文5a$F@9g{kn-Ui5[

990618e18b.gif (1239 bytes)

N:w.pf+r]Y0

整理上式得

?{Nh xi0

Pa+P0-ta-fcta=(1/β)(f0-fc)(1-e-βt0)

&Wwz\Yw1Y0
由图4,次降雨地表径流在产流时刻ta所需的平均临界雨强990618e6a.gif (844 bytes)ta等于下渗率f01,所以可以结合下式计算出产流时刻

990618e19a.gif (1446 bytes)

rE/Vm:k"\1A \y0

整理得

%`.WPa?7Y0

990618e19b.gif (1174 bytes)

7mt9VW6l4P0

(19)

f/af-K+I G,thV;[)vQn0

代入式(18),得

*c{[ Cw!A+v0

990618e20.gif (1287 bytes)

"s(i;n7QO9?0

(20)水利论文(Il6LS bvJ

由式(20)可知,产流时刻不仅与前期土壤含水量Pa有关,而且与平均雨强990618e6a.gif (844 bytes)ta有关(即与产流前的降雨量有关)。可以用试算法找出式(20)中的ta990618e6a.gif (844 bytes)ta

  次降雨地表产流量为水利论文 x q)?+s~P!V

990618e21.gif (1781 bytes)水利论文'E R:B+l!EQ)CN

(21)水利论文)AJt^f-|

  将式(20)中通过试算法找出的ta990618e6a.gif (844 bytes)ta代入式(21)即可求出次降雨地表产流量。时段(Δtj=tj-tj-1)地表产流量的计算则以tjtj-1及相应历时的降雨量分别代入式(21)后相减即得。

@mg6x9Pq0
  根据以上的计算公式,对子洲径流实验站团山沟站1969511日的暴雨进行计算,其下渗与雨强历时曲线见图5,土壤初期入渗率为2.02mm/min(fc=0.42mm/mim,β=0.0538min-1),求得产流时刻在11.8分;各时段产流总和为16.6mm,与实测值15.5mm很接近,误差为7.1%。对1969720日的(前期土壤含水量为14mm)暴雨进行计算,土壤初期入渗率为1.47mm/min,产流总和为0.94mm,实测值0.9mm,误差为4.44%。计算产流量均略大于实测值,这主要是因为用下渗曲线法计算的流域次暴雨产流量是流域的当地产流量,而实测流域出口的产流量是经流域汇流后的流量,而在汇流过程中必有一定的填凹等损失。所以,可以认为计算结果基本上是符合实际的。

6p44.gif (4188 bytes)

c \%ru+`7aHF0

图5 下渗与雨强历时曲线(1969.5.11)水利论文"u9fa$Oay4o7B&v

5结论水利论文}%q#TP+v.LKBR

  1.水土保持措施对土壤水文性质有一定的改良作用,反映在下渗曲线的上移和形状的变化上,从而增加了土壤入渗量,将地表径流更多地转化为地下径流,可作为治理流域,减少洪水和泥沙灾害的战略性措施。

ly)TKavs1cuy8t0

  2.本文在超渗产流计算中考虑到流域实际的下渗能力,即建立了一个将下渗能力曲线、不稳定降雨过程和前期土壤含水量有机联系起来的流域产流计算方法。水利论文,VL8E$U @'J

参 考 文 献水利论文U%_!GT.~1[O

1王礼先主编。水土保持学。中国林业出版社,1995.水利论文J/f5^ d Xs"H,mnh#Lv

2汤立群,陈国祥。流域产沙动力学模型。全国泥沙基本理论研究学术讨论会论文集,1992.水利论文:z ekP Y4d5Vn

3余新晓。土壤动力水文学及其应用。中国林业出版社,1995.水利论文C.N s^9U-J k ~

4赵人俊,王佩兰。霍顿与菲利蒲下渗公式对子洲径流站资料的拟合。人民黄河1982,(1)

AL4w1m2Lr[0K0

 

J6^jhS b0
TAG: 曹文洪 计算方法 王秀英 入渗 产流
顶:11 踩:15
【已经有70人表态】
7票
极差
7票
很差
11票
较差
8票
稍差
9票
稍好
9票
较好
7票
很好
12票
极好
下一篇:分流对高含沙输送影响试验的初步分析(梁志勇,王兆印,匡尚富,徐永年)
上一篇:风沙运动规律的初步研究(陈东 曹文洪 傅玲燕 程义吉)
查看全部回复【已有0位网友发表了看法】