皇甫川流域产流产沙数学模型及水沙变化原因分析(王向东谢树南陈海迟)

热度64票 浏览76次 【共0条评论】【我要评论 时间:1999年9月01日 14:42

皇甫川流域产流产沙数学模型及水沙变化原因分析水利论文x7p8PX#P P

王向东谢树南陈海迟
,_H)I0E!{)H%V0(中国水利水电科学研究院;清华大学 水电系;内蒙古水利学校)

w{2Jv#^@R0

摘 要  该模型是从泥沙运动力学的基本理论出发,结合水文学、气象学、土壤学和地理地貌学的原理建立起来的。并着重解决了在大流域计算中流域的地理、地貌的差异性和降雨在时空分布上的不均匀性问题。利用前人采用的遥感技术资料和实地调查资料,将皇甫川流域分为三大流域、14个小流域、510个计算小单元。在径流系数的选取上考虑了前期降雨的影响。经过对1966年至1969年59组944车雨、径流和输沙资料的验证,计算精度良好。应用该模型对皇甫川流域1970年至1989年次暴雨产流产沙进行了计算,分析了皇甫川流域水沙变化原因。结果表明数学模型的计算成果基本合理。水利论文 tj;\1qhy,D0B

关键词  皇甫川流域;暴雨;数学模型;水沙变化

L[7Wa:d0

1 皇甫川流域基本概况水利论文 RR6K Ij

  皇甫川流域位于黄河中游龙门一河口区间的上段,是黄河粗泥沙产区及泥沙来源的主要地区,年平均向黄河输送约0.5亿吨的泥沙。主流长125km,流域面积3240km。皇甫川流域水系主要由干流纳林川和支流长川组成。按不同地表物质分布和侵蚀差异可将流域划分为三个类型区: (1)黄土丘陵沟壑区。主要分布于流域的东部和西南部,例如十里长川以东和中间部位,面积为1627km2,沟壑密度为5~9km/km2,本地除部分梁峁和缓坡地为耕地外,多为天然草场,植被度为20%左右。本区黄土较厚,约20~30m之间,是陕晋黄土的边缘,呈现较典型的黄土梁峁和黄土沟谷地貌。该区以水蚀为主,水蚀、风蚀和重力侵蚀交错,(2)沙土丘陵沟壑区。主要分布于纳林川中下游以东到长川以西地区和库布其沙漠边缘。面积为243km2,平均沟壑密度为4.2km/km2,该区地势较平缓,地形较完整,沟道较浅。表层原为黄土,由于长期开荒及过度放牧的结果,天然植被严重破坏,经风蚀,引起沙化。该区的主要植被有沙蒿、柠条和沙棘等,覆盖度为15%左右。该区水蚀较轻,风蚀为主要侵蚀方式。(3)砒砂岩丘陵沟壑区。主要分布于流域的西北部,例如干察板沟,圪秋沟,虎石沟和尔架麻沟。面积为1379km2,沟壑密度平均为7.42km/km2,该区水土流失极其严重,地形切割十分破碎,坡陡沟深,本区覆盖度很低,基岩大面积外露。该区侵蚀以水蚀为主复合重力侵蚀。
6o,}B8q8zF O4D0  流域内降雨的年际和年内变率大,年内降雨集中于6~9月,占全年的80%以上。且常以暴雨出现,历时短,强度大,一次降雨往往集中在几分钟至几十分种分内。对流域面雨量的分析计算,多年平均雨量为379mm,1959年年雨量高达685mm,1965年只有192mm,两者相差三倍多。水利论文&JDt'xt!e U
  皇甫川出口站皇甫站多年平均径流量为1.7亿m3,输沙量为5676万吨,多年平均含沙量320kg/m3。径流和输沙量在年内分配极不均匀,沙量尤甚,主要集中于汛期(6~9月)。汛期沙量占全年沙量的98.9%,汛期水量占全年水量的82.1%。流域历年最大洪峰流量高达1.16万m3/s,洪水均为暴雨形成。皇甫川经常出现高含沙洪水,最大含沙量为1270km/m3。洪水平均含沙量达422km/m3。悬沙粒径大于0.05mm,占54.5%,大于0.05mm粒径的沙量多年平均为3093万吨,皇甫站多年平均输沙模数为1.78万t/km2·a。水利论文:GH-IiS r"^(L

2 次暴雨产流产沙计算方法研究及模型建立

6e\7VX | U4w4t0

  研究流域产流产沙规律,需要遵循成因分析的途径,在泥沙运动力学理论的基础上,结合水文学、气象学、土壤学和地质地貌学的基本原理,研究包括降雨、地形、土壤、植被、工程等诸多因素的影响和流域产流产沙的物理图形及定量关系。只有在此基础上建立流域次暴雨产流产沙的数学模型,才能实际地反映流域的产流产沙过程和水沙变化。流域产流产沙过程非常复杂,其影响因素很多,有气候因子(降雨量、降雨强度、蒸发等)、下垫面情况(土壤性质、地形、植被等)和人类活动影响等。
i9V]g+c9hp8rz k:eY0  研究未经治理情况下的流域产流产沙规律,需要突出降雨和产流产沙之间的矛盾。作者把皇甫川流域根据具体情况分为三个计算时段:1966至1969年作为第一计算时段,这四年期间流域基本未经治理或水土保持措施尚未发挥作用,作为公式推导和验证;1970至1979年为第二计算时段;1980至1989年为第三计算时段。第二和第三计算时段作为公式应用。根据“黄河皇甫川流域土壤侵蚀系统模型和治理模式”一文中已有的资料,将皇甫川流域划分为三个小流域,即沙圪堵流域、长滩流域和皇甫区间(即沙圪堵流域和长滩流域以下区域);在此基础上,又划分为14个子流域;每个子流域再划分为若干个类型不相同的单元,共510个单元,而每个单元作为模型的基本计算单元。水利论文)j E0`-LbW)Q
2.1 建立次暴雨产流模型
D%a;{h;Wt(r){"j0
  对本流域1963~1973年85场降雨径流分析研究,取Ie=0.2mm/min和Pe=6.0mm作为产流产沙的临界值,即认为降雨开始和结束时雨强小于0.2mm/min和雨量小于6.0mm时不产流产沙,忽略这两部分降雨,而取中间的一段时间作为产流历时,取产流历时内的降雨量和降雨强度分别作为该场暴雨的有效降雨量和有效降雨强度(以下简称降雨量和降雨强度)。水利论文0UB5tgMmV

  将流域采用泰森法划分成由若干个典型雨量站控制的区域,共分12个区域。求出各个流域内的各个区域的面积及面积权重αi。按照面积权重αi计算出次暴雨的加权降雨量P和加权雨强I。根据本流域1966~1969年59组次暴雨径流资料分析得知,径流系数不单是某个单因子的函数,而是降雨量(P)、降雨强度(I)和前期影响雨量(Pa)三因子综合的函数。
0g ?` Y!y fp.?0  图1是径流系数与有效降雨量、降雨强度和前期影响雨量的相关分析图。有两种变化趋势,经研究,当I<10mm/hr时,相关系数为0.89,标准差为0.32;当I>10mm/hr时,相关系数为0.96,标准差为0.26。因此在建立模型时,采用了两种指标进行了分别统计

t0901.gif (2286 bytes)水利论文Z;u"f h9QP*pv|

图1 f和P、I、Pa的关系水利论文%{)k}2}6A

当I<10mm/hr时水利论文0E.@Aw-S,G'a

f=0.0006017 P0.91I1.44Pa0.14

!N6R.I/dDO9cLu0

(1)水利论文hC2@ wX'V+I

当I>10mm/hr时

|4lP5S0x BV3p0

f=0.002693 P0.38I0.90

,o `N7Q,RQlk0

(2)

2h:Y&e&kp;E-`.P0

  公式(2)中没有考虑前期雨量的影响,这是因为对于高强度降雨,前期雨量的影响较小。
tb&`,IX.fo0
2.2 产沙模型的建立
bG'kS%^P't|m0  关于产沙公式,到目前为止已有许多成果,但仍存在许多问题,国外在这方面的发展较国内迅速,但直接用于我国黄土高原产沙计算是不合适的。国内产沙公式大多是从研究黄土高原上的实验数据或观测资料回归分析得到的经验公式,对于研究产沙的物理过程方面主要是在90年代,仍带有深刻经验烙印。为此,皇甫川流域产沙模型在谢树楠等的产沙模型基础上,结合水文、气象、土壤和地质地貌学的基本原理,着重解决了流域的不均匀性和K值率定等问题建立的。水利论文:ia9|;`%O~1ud
2.2.1 基本假定水利论文Q2RtRWy}W R
  (1)暴雨产生的径流为坡面一维流动;(2)压强按静水压强公布;(3)流动中的动量系数为常量;(4)泥沙不考虑粘性;(5)沟道泥沙输移比为1。
-L%`8P*a.qs1Hk0
2.2.2 公式推导水利论文G2|s1bufoDlU?9];j
  (1)挟沙力公式 根据爱因斯坦挟沙力的概念,水力参数Θ和输沙参数Φ的关系可以写为水利论文3fR ~:\#K k_

903.gif (955 bytes)

j%SWd.t0

(3)

i6RII&`+i N1N{0

水利论文K%l2z@9n(n'D `P

904.gif (1462 bytes)

(4)水利论文/X3C#f-m qYZO%P W

905.gif (1201 bytes)

7i!\lmN)P0

(5)

N4YsZn}A0

将(3)式和(4)式代入(5)式,整理合并为

Nh6}NM0

906.gif (1052 bytes)水利论文_0Dta2cq

(6)水利论文6W6sX$X1_F?

式中 β1=αg-β+1/2/{[(γs-γ)/γ]g}β+1/2D-β+3/2;qs为输沙率;γs,γ为分别为泥沙和水的比重;g为重力加速度;τ0为泥沙的中数粒径;α为常数;β为指数。水利论文^^e)_slA&]7m
  由式(6)可见,坡面输沙率主要与作用在坡面的水流剪切力成正比。
:b|2@A+y~~$J+A7k.Gx0  (2)泥沙输移的连续方程水利论文-j?ggO~%h

dqs/dx=Ps水利论文K8L}^F4fW%|E0i

(7)水利论文JiS0z Ey

式中Ps单位面积上的起冲率;x为在流动方向上沿坡面的距离。水利论文gR#N#T{VD
  (3)水流的连续方程水利论文-LgKE TOW

q=q0+q*x水利论文j a:|s1iM-u%t

(8)水利论文.sPTl!C a

式中 q为单宽流量;q0为斜坡上段流入的单宽流量; q*为单位时间、单位面积斜坡上由降雨产生的径流量。水利论文/E ?GxL_Sg.?

  (4)水流运动方程 对于明渠,具有侧向入流的水流运动方程

:R(B1R5_pa AHi,fND0

dy/dx=[s0-sf-(2β2q*q)/gy2+q*Ucosφ/gy]/(cosθ-β2q2/gy2)水利论文0w3V"U9p-q

(9)

r"yu Or}0

式中 y为水深;s0,sf为坡面坡度和能坡;β2为动量因子;U为进入主流中的侧向水流的流速(即雨点的终极速度);φ为x方向与速度U的夹角;θ为x方向与水平方向之间的夹角
/A }E)~+V U5Vu0  (5)泥沙起冲率与边界剪切力 泥沙的起冲率可以表达为边界剪切力的函数。即水利论文h!s zxt!r

Ps=aτb0水利论文o,Z!| Ukv4v5|M

(10)

(X!OI^X0p;Wg0

式中 a为描述坡面土壤侵蚀特性的系数;b为指数。水利论文qj.@1i$yF(^8F_
将(6)和(7)联立求导,整理合并得

%lY;x!D bev4m SD0
Ps=(ββ1β)(dτ0/dx)τβ-1

(11)

%_#~v1biSG0az0

与(10)式比较得出水利论文]N_)Q)S G.f

a =(ββ1β)(dτ0/dx)水利论文%\5bh%b ^5OU

(12)

v#wI f.c4qWo0

b=β-1

4|)O$udxl v0

(13)

xWpc{~/N\4AN0

又∵τ0=γysf水利论文`6hO&B$|

7z;y,[W }K@0
914.gif (1092 bytes)

(14)

F]Wu9ao0

其中 γ为水流容量;914-1.gif (875 bytes)为平均摩阻流速。水利论文!n*km.M)]1tq|l

将(14)式代入(12)式得水利论文joY#^ ct ^W

915.gif (1156 bytes)

RFG8g/c["~0

(15)水利论文7HO eGME&TY g

  (6)径流深的求解 一般说来,对于坡面流,水深很小,可近似认为是均匀流,则(9)式可简化为

Q+Gb!kw R0

sf=s0-(2β2q*q)/gy2+q*Ucosφ/gy水利论文4`e.^%?1Zj

  φ为x方向与速度U的夹角,可以考虑降雨与径流垂直,即cosφ=0,则上式可进一步简化为

K4f9C5B#^ a1y W!G0

sf=s0-(2β2q*q)/gy2

-W$b.M j W`:bo/Hi'[0

则径流深为

,b W*sj-~ Y+|s*[0

y=[(2β2q*q)/g(s0-sf)]1/2

B-Y.|J3Z_!{0

(16)水利论文n9[c3L_$}e

对(16)式微分,得水利论文[H(f|+`Be9J

dy/dx=1/2q*3/2[(2β2)/gq(s0-sf)]1/2

N&S LtF0

(17)水利论文 K$PLw'm;u$n.CF R

将(17)式代入(15)式,得水利论文,iU u$j5[{

a =ρββ1sf/2γβ[(2β2)/gq(s0-s)]1/2q*3/2

![i7B;a/QgD0

(18)水利论文0]+iX'l\

  (7)地表降雨下的边界剪切应力 明渠的阻力可表示为

P`6KQ.{o0

918-1.gif (1411 bytes)水利论文ed SqFf k6U

  在坡面上的水流由于降雨的加入,其阻力系数F已不同于明渠阻力,此时降雨强度对决定阻力系数起重要作用。为解决这个问题,70年代初先是南朝鲜的Y.N.Yoon,后是李日明和沈学汶,都先后做过室内的水槽实验。试验表明:当雷诺数Re<2000时,雨强I越大,F就越大;当Re>2000时,雨强I对F不起作用。
4B | Y_3DS6J0  将阻力系数分为两部分:一部分是没有降雨时的阻力系数F0,另一部分是由降雨所增加的阻力系数Fr,认为Fr在光滑和粗糙面上不变。当Re在2000~12600时水利论文.c:c4iCu9A#LWo

F=F0+Fr=(K2/Re1/4)+(0.012/Re1/4

0aB$C V`)s2P&\"t0

(19)水利论文*tG)\.HS7yM-k?2O

由于黄河中游各流域坡度陡,故在计算阻力时可以按(19)式计算。水利论文Hb a]1wo9Pj

水利论文{sl+HpDs9Hu5^

919-1.gif (1925 bytes)

∴

e|C$J QJ*qP%u5i0

920.gif (1427 bytes)

3jm'ZQi;A0

(20)

E s+G;E!R\S0

               ∵y=τ0/γs0

U?/_5_c0[0

水利论文?Q%U?*a

921.gif (1980 bytes)

(21)

J*j I#K&_va3r4O`0

  (8)有关各参数的确定α,β值的确定 谢树楠等分析了团山沟和段川14个径流试验(1966~1967年)的645组资料,建立了水力参数Θ和输沙参数Φ的关系水利论文0F,Mn"^W ?^tN)[y+J

Φ=1.81Θ2.16

6Oa wwNXG ^0

(22)

(iD[pv _ @!p1k#T0

将(22)与(3)式比较,得

-O { m6H6A0

α=1.81; β=2.16

Dsnhl V0

  α,β2值的确定 Yoon把变流速的运动方程写成

TH7yT {"z Y0

923.gif (1735 bytes)

{1ST%d$txp0

(23)水利论文!a2hlw;~%E

  Yoon采用实验资料对等式右边各项的影响程度进行了分析,发现s0对sf的影响最为显著,第二项及第三项对sf的影响几乎相同,而且均相当于s0的1/10,最后一项可以忽略不计,据此可得

DQ4sO1Y w)XnD-Q0
sf= s0--(1/10)s0-(1/10)s0=4/5s0
∴s=s0/sf=1.25

  另据沈学汶及李日明的研究,当雷诺数大于2000时,s=1.05。根据以上研究,Komura建议采用s=1.20,由降雨产生的径流动量因子β2值,可从Yoon的实验结果确定。Yoon等人研究了降雨强度,雷诺数和渠道坡度对动量因子的影响,认为β2值随雷诺数和雨强的增大而减小,但坡度影响不显著。据他们的研究成果,取β2=1.1。水利论文:M:?(T:\2pFWN
  其它参数值的确定水利论文,kixoV M }pA
  据日本学者的研究,认为K2=0.6为宜。水利论文9o VB d8@E

v=1.0×10-6m2/s
γs/γ=2.65

b=β-1=1.16

Dl!q:UH3KM X%A0

  (9)产沙公式 单位面积、单位时间的产沙量Ws可由下式表示水利论文~F'Z.f3[t1Ri$w

924.gif (1310 bytes)

9X^3}ecmP7}$V0

(24)

[ze.k'wR8u0

将(8)、(10)、(13)、(18)和(21)式代入(24)式,积分并整理得

ZH,zEN)Va1@"{0

925.gif (1987 bytes)

+?,TWze0

(25)

oA*O.kpT(~0

若忽略坡顶来流q0,则(25)式可简化为水利论文x'u-Z6D P,A7J9|C1O!R

Ws=CHq*(7b+12)/12L(7b-6)/12S0(4b+3)/6水利论文x]T C;XDF(x)S

(26)水利论文7?L N2c;r`

其中

d.ra@,r{6~j-_M0
926-1.gif (2076 bytes)

把以上各参数代入(26)式,整理得水利论文 H$l Y[9` d

Ws=0.8237q*1.68L0.18S01.27D50-0.66

pOT(Xn+]0

(27)水利论文[,~@%IY,s b,k+n

  暴雨产沙不仅仅与上式右端各因子(即径流量、坡长、坡降和中数粒径等)有关,而且与土壤性质(SB)、土壤侵蚀特性(CE)和裸露率(CA)有关。土壤性质、土壤侵蚀特性和裸露率的不同导致产沙量大小差异。据研究,Ws与SB、CE和CA成正比关系。因此在式(27)中右边乘以SB、CE和CA值,就组成了次暴雨产沙公式。
'zw D$S0S Z;M6U,e(m0式中 q*表示由降雨产生的单位面积径流量,单位是m3/m2·s,实际为单位时间降雨产生的径流深,可表示为:q*=fI水利论文%sb5D5R6K
故次暴雨产沙公式为水利论文}y$wM.s8GJ

Ws=0.82SBCACE(fI)1.68L0.18S01.27D50-0.66(28)

W9Y K%Nf0M0
Wsb=0.82SBbCACE(fI)1.68L0.18S0b1.27D50-0.66 (坡面产沙)
Wsg=0.82SBgCACE(fI)1.68L0.18S0g1.27D50-0.66(沟道产沙)

次暴雨总产沙量为水利论文i0z _c)P:yr

929.gif (1272 bytes)水利论文F,Y5vl+\n

(29)水利论文Wa/],`B` |

上面各式中各因子的含义及单位:Wsi为次暴雨产沙量(104t);Ws为次暴雨的总产沙量(104t);Wsp为坡面产沙量(104t);Wsg为沟道产沙量(104t);SB为土壤种类所占比例(%);SBp为某一土壤类型坡面所占比例;SBg为某一土壤类型沟道所占比例;CA为地表裸露率;CAp为坡面地表裸露率;CAg为沟道地表裸露率;CE为侵蚀因子;CEp为坡面侵蚀因子;CEg为沟道侵蚀因子;D50为土壤的中数粒径(mm);f为径流系数;I为降雨强度(mm/hr);L为坡面长度(m);Lp为单元坡面长度(m);Lg为单元沟道长度(m);S0为坡面比降;S0p为单元坡面平均比降;S0g为单元沟道平均比降;A为单元面积(km2);Δti为计算时间步长(hr)。

X#lQo]'q#EM0

3 模型中各参数的率定水利论文/~ VB/`1O8{.C2KI%v

3.1 径流系数f
9xK8C.~#N B!V.w0
  本模型采用的径流系数f是根据本流域降雨、径流资料,按有效降雨量、有效降雨强度和前期影响雨量进行回归分析所得到的公式:f=0.0006017P0.91I1.44Pa0.14, (I<10mm/hr);f=0.0026930P0.38I0.90, (I>10mm/hr)。水利论文}WX\@@-W WC
3.2 有效降雨强度I水利论文-Fk W(V2X"}(R
  流域产沙不仅取决于降雨量,而且更主要取决于降雨强度及分布。因此,对于产沙的降雨强度不能笼统地用次暴雨的有效降雨量除以有效降雨历时。水利论文#@&D_\+kH8o
  采用谢树楠等研究出的降雨强度分布模拟公式:I(t+Δt)={[(t+Δt)0.366-t0.366]/Δt}PT-0.366
0y9Njg$Lr6E Ac9i0
3.3 有关流域下垫面因子的确定水利论文Bh,w9C kM~)I2FzX
  下垫面因子的确定依据地形图和实际调查资料。
YN?r#L(a6~0  地形图由中国测绘资料信息中心提供,由中国人民解放军总参谋部测绘局1971年11月航测、1977年5月调绘,1979年第一版。比例尺1∶50000,等高距10m。水利论文,R5Loj7Gt P
  流域总面积A及单元面积Ak: 采用工具是求积仪,并结合米格法确定面积。应用求积仪时,顺时针量测一圈,再逆时针量测一圈,取其平均值,最后平差,即为面积。
)I!Rip `6XG7z_0  主沟长度L和单元坡长Lk: 主沟长度L是从流域出口沿主沟道一直延伸至接近分水岭。坡长是坡面的平均长度,对于单元坡长Lk来讲,Lk是定值。
2y.g`2?O1LPoD0  单元沟坡坡度S: 对于流域内不同级别的沟道其沟坡平均坡度是不同的,即使在同一沟道内,上、中、下游各部位的沟坡坡度也有差别。计算产沙时,若采用平均坡度计算,势必带来极大的误差。因此,坡度要采用单元沟坡坡度,这样坡度是均匀的。
^jX O,w"C0  裸露率CA: 裸露率CA就是单元内无植被的面积和整个单元面积的比值,由实测资料知单元的植被覆盖率CB,则CA=1-CB水利论文:| zK-\7L)];K7v
  土壤中数粒径D50:D50由该流域表土的级配曲线率定,不同的土壤类型其D50不同.
.Q:Q(P J2K0  侵蚀因子CE: 侵蚀因子CE是反映流域侵蚀程度的综合系数,它主要决定于流域坡面土的缔固程度、粘土比、分散比及沟壑密度等因素。
wh*E'W~IN0  本流域采用1970年以前(无水保工程前)的多年平均侵蚀模数Ws0作为确定各个单元的侵蚀因子的根据,即:CE=KWs0水利论文5@j eZ%]Y!eyj
  研究皇甫川三个流域59组944场暴雨产沙的K值变化规律,分析结果表明,流域的K值系数为0.00175。

{I)Y'I0H$Dn8@9I0

4 模型验证

FiSD7?;J P`0

  为了验证产流产沙模型的适用性,采用1966~1969年的59组944车雨资料进行验证。三个流域逐次计算出每次暴雨的产流产沙量,再与实测值比较。 沙圪堵流域产流产沙量的计算值与实测值比较,相关系数分别为0.898和0.995;长滩流域和皇甫区间的计算值与实测值比较,相关系数分别为0.752和0.931;皇甫川流域的计算值与实测值比较,相关系数分别为0.967和0.962 (见图2和3)。这三个流域的共同特征是总径流量或总产沙量计算值和实测值相差很小,相关性较好(只有长滩流域和皇甫区间差一些,这是由于其实测值是全流域的径流量减去沙圪堵流域的径流量得出的原因),说明模型精度较高,可以作为本流域产流模型应用。(注:图中Qc表示计算产流量值,Qm表示实测径流量值;Wsc表示计算产沙量值,Wsm表示实测输沙量值。

4N*NO"a2I(PS0

5 产流产沙模型应用及水沙变化原因分析水利论文9Z6f l$^_Z/}

  为了分析流域的水沙变化原因,采用1970~1989年共213组3408场暴雨,逐单元、逐次、分时段进行次暴雨产流产沙计算,然后求其总和,即某次降雨的产流产沙量。年产流产沙量由次暴雨累加,产沙量的计算按坡面、沟道产沙及砒砂岩、黄土和沙土等三种土壤类型分别计算。在分析水沙变化原因时,按三大区即沙圪堵站以上流域、长滩站以上流域及皇甫区间和全流域分别统计。

{S1~^ZXt!K K0

t0902.gif (3789 bytes)

!FD6jS6n"A$v0

t0903.gif (4302 bytes)水利论文P^jm:`:tg'C3f

图2 全流域次暴雨产流量验证结果

b$Y'J)p%fB~@o0

图3 全流域次暴雨产沙量验证结果水利论文D s x(O3v%G6Bc Fy

5.1 应用产流模型分析水量变化原因水利论文R&w,J3b^:us
  沙圪堵流域在进入70年代以后:(1)降雨影响水量变化:70年代径流量增加了21857.1万m3,年均2185.7万m3,占39.09%;80年代径流量减少了5335.9万m3,年均533.6万m3,占9.54%;1970~1989年径流量增加了8260.6万m3,年均826.1万m3,占14.77%。(2)人类活动和其他因素影响水量变化:70年代径流量减少了184.8万m3,年均18.5万m3,占0。33%;80年代径流量减少了3179.9万m3,年均318.0万m3,占5.69%;1970~1989年径流量减少了1682.3万m3,年均168.2万m3,占3.01%。详见表1。说明该流域70年代和80年代人类活动和其他因素影响(主要是水保工程作用)水量减少甚微。
Z4^ ee'\$A0  长滩流域和皇甫区间在进入70年代以后:(1)降雨影响水量变化:70年代径流量增加了5474.5万m3,占7.45%;80年代径流量减少了13603.6万m3,占18.51%;1970~1989年径流量减少了4064.5万m3,占5.53% 。(2)人类活动和其他因素影响水量变化:70年代径流量减少12037.8万m3,占16.38%;80年代径流量减少8658.0万m3,占11.78%;1970~1989年径流量减少了10347.9万m3,占14.08%。详见表1。说明该流域人类活动和其他因素影响(主要是水保工程)的减水作用较显著。
9Nd"[X0NP0  因此,整个流域进入70年代以后,有一定的减水效益:(1)降雨影响水量变化:70年代径流量增加了27331.6万m3,年均2733.2万m3,占21.12%;80年代径流量减少了18939。5万m3,年均1894万m3,占14.64%;1970~1989年径流量增加了4196.1万m3,年均419.6万m3,占3.24%。(2)人类活动和其他因素影响水沙量变化:70年代径流量减少了12222.6万m3,年均1222.3万m3,占9.45%;80年代径流量减少了11837.9万m3,年均1183.8万m3,占9.15%;1970~1989年径流量减少了11880.3万m3,年均1188.0万m3,占9.18%(详见表1)。说明皇甫川流域进入70年代以后,人类活动和其他因素影响(主要是水保工程)水量有一定的减水作用,但不明显。沙圪堵流域水保工程较少,减水作用很小;长滩和皇甫区间由于有较多的水土保持工程,使得减水效益较高。水利论文 jaj7K4Uu[
5.2 应用产沙模型分析沙量变化原因
8Y/V%xO lQ)h5D05.2.1 年降雨产沙量的计算方法水利论文7Sz'rf W
  次暴雨的产沙量总和基本上代表了年产沙量。沙圪堵流域汛期实测输沙量占年输沙量的94.2%;皇甫川流域汛期实测输沙量占年输沙量的98.9%。故采用由次暴雨产沙累积求年产沙量的办法水利论文+DiR,H%Y

表1 皇甫川流域水沙变化水量分析成果表(单位:104m3)

-HG]:D%J*F"n0

年雨量计算实测降雨影响其他因素影响
(mm)产流量径流量差值占基准年%差值占基准年%
(1) (2)(3)(4)(5)(6)(7)(8)

沙 圪 堵 流 域水利论文-\S2\Lvy

60年代389.855911.155911.1*
70年代384.777768.277583.4-21857.1-39.09184.80.33
80年代304.350575.247395.35335.99.543179.95.69
70、80年代344.564171.762489.4-8260.6-14.771682.33.01

长 滩 流 域 和 皇 甫 区 间水利论文*g ~ruA)g;O6{(F!ok

60年代405.073482.773482.7*   
70年代377.478957.266919.4-5477.5-7.4512037.816.38
80年代363.959879.151221.113603.618.518658.011.78
70、80年代370.769418.259070.34064.55.5310347.914.08

全 流 域

.Vx3_G[-Ij)j,nSM0
60年代419.8129393.8129393.8*
70年代378.8156725.4144502.8-27331.6-21.1212222.69.45
80年代336.3110454.398616.418939.514.6411837.99.15
70、80年代357.6133589.9121709.6-4196.1-3.2411880.39.18

  注:(5)=(4*)-(3),(6)=(5)/(4*),(7)=(3)-(4*),(8)=(7)/(4*)

Z:zN,D |2T0

5.2.2 结果分析水利论文nd2S&A9X$b
  70年代,沙圪堵流域坡面产沙7773万吨,占总产沙量的29.4%,沟道产沙18683万吨,占总产沙量70.6%;而砒砂岩类产沙11383万吨,占总产沙量的43%,黄土类产沙13774万吨,占总产沙量的52.1%,沙土类产沙1297万吨,占4.9%(参见表2)。70年代实测总沙量比60年代增沙13969万吨,其中降雨因素增加6708万吨,占48.0%,而人类活动和其他因素增沙7262万吨,占52.0%(参见表3)。80年代,沙圪堵流域坡面产沙4483万吨,占总产沙量的28%,沟道产沙11527万吨,占72%;而砒砂岩类产沙6767万吨,占总产沙量的42%,黄土类产沙8333万吨,占总产沙量的52%,沙土类产沙910万吨,约占6%。80年代比60年代增沙390万吨,其中降雨因素减沙3738万吨, 人类活动和其他因素增沙4128万吨。80年代沙圪堵流域增沙主要是由于人类活动和其他因素引起的。说明沙圪堵流域增沙的趋势应引起足够的重视。水利论文k6_bkf J/T&m+u"SM
  长滩流域及皇甫区间70年代该流域坡面产沙5028万吨,占总产沙量的15.1%, 沟道产沙28349万吨,占84.9%,而砒砂岩类产沙10432万吨,占31.3%,黄土类产沙20924万吨,占62.7%,沙土类产沙2020万吨,占6.1%(参见表2)。70年代比60年代减沙1224万吨,其中降雨因素增沙3779万吨,而水保工程减沙5002万吨,说明水保工程的作用是显著的(参见表3)。80年代该流域坡面产沙4088万吨,占总产沙量的14.6%,沟道产沙24003万吨,占85.4%。砒砂岩类产沙8784万吨,占31.3%,黄土类产沙17559万吨,占62.5%,沙土类产沙1747万吨,占6.2%。80年代比60年代减沙8350万吨,其中降雨减沙1507万吨,而水保工程减沙6843万吨。说明80年代水保工程的减沙作用比70年代更好。

n"tf%e!R)o!\c9M)}H0

表2 皇甫川不同土壤类型坡面及沟道产沙量统计略表(单位:万吨)水利论文K%w,o7x*fy F1m


年代

沟 道 产 沙水利论文+cCeW!dh3G8X%q hD


坡 面 产 沙水利论文gT2X4J8{,h H/n


总计水利论文grc q c(r'd b6{


砒砂岩黄土类沙土类砒砂岩黄土类沙土类计算值实测值计—实

沙 圪 堵 流 域
1970~19798110.79480.61092.43273.04294.6205.126456.533718.4-7261.9
1980~19894876.25869.3781.41890.62464.1129.016010.720138.7-4128.0

长 滩 流 域 和 皇 甫 区 间水利论文)IdY`0U(I

1970~19798646.917959.81742.71785.52965.2277.133377.128374.95002.3
1980~19897346.915130.41526.21438.42429.0221.228092.121249.16843.0

全 流 域水利论文#Usx@8G8A8g5Y

1970~197916757.627440.42835.15058.57259.8482.259833.662093.3-2259.6
1980~198912223.120999.72307.63329.04893.1350.244102.841387.82715.0


L~Cw(uY0表3 皇甫川流域水沙变化水量分析成果表(单位:万吨)

|'`#^ qbd+a.Ba?0

时段水利论文8b]q|2r5gA!l

年雨量计算实测降雨影响其他因素影响
(mm)产流量输沙量差值占基准年%差值占基准年%
(1) (2)(3)(4)(5)(6)(7)(8)

沙 圪 堵 流 域水利论文(gR[K;|

1966~1969389.819748.919748.9*
1970~1979384.726456.633718.5-6707.7-33.96-7261.8-36.77
1980~1989304.316010.820138.73738.118.93-4128.0-20.90
1970~1989344.542467.453857.2-1484.8-7.52-5694.9-28.84

长 滩 流 域 和 皇 甫 区 间水利论文qZ [[{ E5Om%N(_g

1966~1969405.029598.729598.7*   
1970~1979377.433377.228374.9-3778.5-12.775002.216.90
1980~1989363.928092.121249.01506.65.096843.023.12
1970~1989370.761469.349623.9-1136.0-3.845922.620.01

全 流 域

6l.Z[F/d6X'IDmN0
1966~1969419.849347.649347.6*
1970~1979378.859833.562093.3-10485.9 TAG: 数学模型 陈海迟
顶:11 踩:8
【已经有45人表态】
10票
极差
3票
很差
4票
较差
7票
稍差
8票
稍好
4票
较好
4票
很好
5票
极好
下一篇:大凌河白石水库淤积分析(金宝琛,王立强,刘宇聪)
上一篇:边坡形状对崩塌的影响(徐永年,匡尚富,李文武,王力)
查看全部回复【已有0位网友发表了看法】