流域产流与汇流计算.docx
《流域产流与汇流计算.docx》由会员分享，可在线阅读，更多相关《流域产流与汇流计算.docx（50页珍藏版）》请在冰豆网上搜索。
流域产流与汇流计算
流域产流与汇流计算（总47
页）
本页仅作为文档页封面，使用时可以删除
Thisdocumentisforreferenceonly-rar21year.March
第四章流域产流与汇流计算
第一节概述
根据笫二章的论述，山降雨形成流域出口断面径流的过程是非常复杂的，为了进行定量阐述，将这一过程概化为产流和汇流两个阶段进行讨论。
实际上，在流域降雨径流形成过程中，产流和汇流过程儿乎是同时发生的，在这里提到的所谓产流阶段和汇流阶段，并不是时间顺序含义上的询后两个阶段，仅仅是对流域径流形成过程的概化，以便根据产流和汇流的特性，采用不同的原理和方法分别进行计算。
产流阶段是指降雨经植物截留、填洼、下渗的损失过程。
降雨扣除这些损失后，剩余的部分称为净雨，净雨在数量上等于它所形成的径流量，净雨量的计算称为产流计算。
由流域降雨量推求径流量，必须具备流域产流方案。
产流方案是对流域降雨径流之间关系的定量描述，可以是数学方程也可以是图表形式。
产流方案的制定需充分利用实测的流域降雨、蒸发和径流资料，根据流域的产流模式，分析建立流域降雨径流之间的定量关系。
汇流阶段是指净雨沿地面和地下汇入河网，并经河网汇集形成流域出口断面流量的过程。
山净雨推求流域出口断面流量过程称为汇流讣算。
流域汇流过程乂可以分为两个阶段，由净雨经地面或地下汇入河网的过程称为坡面汇流；进入河网的水流自上游向下游运动，经流域出口断面流出的过程称为河网汇流。
山净雨推求流域出口流量过程，必须具备流域汇流方案。
流域汇流方案是根据流域净雨讣算流域出口断面流量过程，应根据流域雨量、流量及下垫面特征等资料条件及计算要求制定。
就径流的来源而论，流域出口断面的流量过程是山地面径流、壤中流、浅层地下径流和深层地下径流组成的，这四类径流的汇流特性是有差别的。
在常规的汇流计算中，为了计算简便，常将径流概化为直接径流和地下径流两种水
源。
地面径流和壤中流在坡面汇流过程中经常相互交换，且相对于河网汇流，坡面汇流速度较快，儿乎是直接进入河网，故可以合并考虑，称为直接径流，但在很多情况仍称为地面径流。
浅层地下径流和深层地下径流合称为地下径流，其特点是坡面汇流速度较慢，常持续数十天乃至数年之久。
U前，在一些描述降雨径流的流域水文模型中，为了更确切地反映流域径流形成的过程，采用了三水源或四水源进行模拟讣算。
第二节流域降雨径流要素的计算
、流域降雨量
（一）流域平均雨量计算
实测雨量只代表雨量站所在地的点雨量，分析流域降雨径流关系需要考虑全流域平均雨量。
一个流域一般会有若干个雨量站，由各站的点雨量可以推求流域平均降雨量，常用的方法有算术平均法、垂直平分法和等雨量线法三种。
算术平均法：
当流域内雨量站分布较均匀且地形起伏变化不大时，可根据各站同时段观测的降雨量用算术平均法推求流域平均降雨量
（4-1）
式中P——流域某时段平均降雨量，mm：
P<流域内笫i个雨量站同时段降雨量，mm；
n——流域内雨量站点数。
垂直平分法：
也称为泰森多边形法，适用于地形起伏变化不大的流域。
这一方法假定流域内各处的雨量可山与之距离最近站点的雨量代表，如图4-1所示。
具体做法是先用直线连接相邻雨量站，构成门-2个三角形（最好是锐角三角形），再作每个三角形各边的垂直平分线，将流域划分成门个多边形，每一多边形内均含有一个雨量站，按多边形面积为权重推求流域平均降雨量
_1”
P=（牛2）
屮r-1
式中力一一第i个雨量站所在多边形的面积，km2；
F流域面积，km2o
等雨量线法：
当流域内雨量站分布较密时，可根据各站同时段雨量绘制等
雨量线（见图4-2）,然后推算流域平均降雨量
（牛3）
式中fj——相邻两条等雨量线间的面积，km2；
Pj——相应面积齐上的平均雨深，一般采用相邻两条等雨量线的平均值，
mm：
m——分块面积数。
（二）雨量过程线
降雨强度过程线：
降雨强度随时间的变化过程线称为降雨强度过程线，通常以时段平均雨强为纵坐标，降雨时程为横坐标的柱状图表示，如图4-3。
如果以时段雨量为纵坐标，则称为雨量过程线，也称为雨量直方图。
累积雨量过程线：
自降雨开始起至各时刻降雨量的累积值随时间的变化过程线，称为累积雨量过程线，如图4・4。
图4-3雨量过程线图4-4累积雨量过程线
山降雨强度过程线转换成累积雨量的公式为：
匕=£（4-4）
式中Pj——至第j时段末的累计雨量，mm;
ik第k时段的降雨强度，mm/h;
△t时段长度，h
反之，根据累积雨量推求时段降雨强度的公式为：
（4-5）
P）-S
A/
二、径流量
流域出口流量过程线除本次降雨形成的径流以外，往往还包括前期降雨径流中尚未退完的水量，在计算本次径流时，应把这部分水量从流量过程线中分割出去。
此外，III于不同水源成分的水流运动规律是不相同的，需对流量过程线中的不同水源进行划分，以便进行汇流计算。
（1）流量过程线的分割
流域蓄水量的消退过程线称为退水曲线，不同次降雨形成的流量过程线的
分割常采用退水曲线。
取多次实测洪水过程的退水部分，绘在透明纸上，然后
沿时间轴平移，使它们的尾部重合，形成一簇退水线，作光滑的下包线，就是流域地下水退水曲线，如图4・5。
有了退水曲线，就可以将各次降雨所形成流量过程线分割，如图4・6,得出对应于本次降雨所形成的流量过程线。
Kh）
图4・5流域退水曲线
流域地下径流退水过程比较稳定且时间较长，地下水退水曲线可以用下式来描述：
0（f）=0（0）严（4-6）
式中Q（t）——r时刻地下水流量；
Q（0）——初始地下水流量；
kg——地下水退水参数，可利用地下水退水曲线来率定。
（2）径流量计算
实测流量过程线割去非本次降雨形成的径流后，可以得出本次降雨形成的
流量过程线。
据此，推求出相应的径流深：
3.6工Q△/
/?
=—
F
（4-7）
式中R径流深，mm；
At——时段长度，h；
Q——第j时段末的流量值，m3/s；
F流域面积，km2o
（三）水源的划分
地面径流和地下径流汇流特性不同，求得次径流总量之后，还需划分地面径流和地下径流。
简便的划分方法是斜线分割法，从流量起涨点到地面径流终止点之间连一直线，直线以上部分为地面径流，直线以下部分为地下径流，如图4-7所示。
地面径流终止点可以用流域退水曲线来确定，使退水曲线的尾部
与流量过程线退水段尾部重合，分离点即为地面径流终止点。
为了避免人为分
析误差，地面径流终止点也可用经验公式确定。
例如，某区域的经验公式为：
（4-8）
式中N——洪峰出现时刻至地面径流终止点的日数:
F——流域面积，km2o
三、土壤含水量
（一）流域土壤含水量的计算
降雨开始时，流域内包气带土壤含水量的大小是影响降雨形成径流过程的一个重要因素，在同等降雨条件下，土壤含水量大则产生的径流量大，反之则小。
流域土壤含水量一般是根据流域询期降雨、蒸发及径流过程，依据水量平衡原理采用递推公式推求：
Wt+i=Wt+Pt-ErRt（4-9）
式中Wt一一笫r时段初始时刻土壤含水量，mm;
Pt——第r时段降雨量，mm;
Et第t时段蒸发量，mm;
Rt第r时段产流量，mm。
流域土壤含水量的上限称为流域蓄水容量必，山于雨量、蒸发量及流量的观测与计算误差，采用公式4-9的讣算出的流域土壤含水量有可能大于必或小于0的情况，这是不合理的，因此还需附加一个限制条件：
OS必必。
采用公式4-9需确定合适的起始时刻及相应土壤含水量。
可以选择前期流域出现大暴雨的次日作为起始日，相应的土壤含水量为必；或选择流域长时间干旱期作为起始日，相应的土壤含水量取为0或较小值；也可以提前较长时间（如15〜30天）作为起始日，假定一个土壤含水量（如取必值的一半）作为初值，经过较长时间计算后，误差会减小到允许的程度。
（二）流域蒸发量
流域蒸发量的大小主要决定于气象要素及土壤湿度，这可以用流域蒸发能力和土壤含水量来表征。
流域蒸发能力是在当日气象条件下流域蒸发量的上
限，一般无法通过观测途径直接获得，可以根据当日水面蒸发观测值通过折算
间接获得：
（4-10）
Em=6E0
式中Em——流域蒸发能力；
Eo——水面蒸发观测值；
6——折算系数。
我国水利部门常用的流域蒸发量计算模式有三种。
1・一层蒸发模式：
假定流域蒸发量与流域土壤含水量成正比
（4-11）
（4-12）
F
E=^W
一层蒸发模式比较简单，但没有考虑土壤水分的垂直分布情况。
当包气带土壤含水量较小，而表层土壤含水量较大时，按一层蒸发模式得岀讣算值偏小，例如，久旱后降了一场小雨，其雨量仅补充了表层土壤含水量，就是这种情况。
2.二层蒸发模式：
将流域蓄水容量必分为上层和下层WS，相应的土壤含水量分别WU和假定降雨量先补充上层土壤含水量，当上层土壤含水
量达WUn后再补充下层土壤含水量；蒸发则先消耗上层土壤含水量，蒸发完了
再消耗下层的土壤含水量,
且上层蒸发EU按流域蒸发能力蒸发,
下层蒸发豆与
EU=
下层土壤含水量成正比，即:
（4-13）
EL=
WL
（E」EU）
（4-14）
流域蒸发量为上下二层蒸发量之和:
（4-15）
E=EU+EL
二层蒸发模式仍存在一个问题，即久旱以后由于下层土壤含水量很小，计
算岀的蒸发量很小，流域土壤含水量难以达到凋萎含水量，不太符合实际情况。
3•三层蒸发模式：
在二层蒸发模式的基础上，确定了一个下层最小蒸发系
数C,上层蒸发仍按公式4-13II-算，下层蒸发按下式计算:
当WL>C（Em时
EL=<
WI
页2）
（4-16）
C（E枫-EU）
当WL（4-17）
EL=WL
（三）前期影响雨量
在很多情况下，采用式4-9推求土壤含水量时，会遭遇径流资料缺乏的问题。
在生产实际中常采用前期影响雨量Pa来替代土壤含水量，计算公式为
Pa.t+1=K（Pa.t+Pt）（牛18）
式牛18的限制条件为PaWm时取Pa=Wm.
在式4J8中，K是与流域蒸发量有关的土壤含水量日消退系数。
如果采用采用一层蒸发模式，对于无雨日：
E
“-E尸（1-—）P“（4J9）
对照无雨日时的公式4-18,即Pa.t+1=KPa.t，可知：
E
K=1一一（牛20）
如果在某一时间段，％取一平均值，则在该时间段的K为常数。
第三节蓄满产流计算
蓄满产流模式
在湿润地区，由于雨量充沛，地下水位较高，包气带较薄，包气带下部含水量经常保持在田间持水量。
在汛期，包气带的缺水量很容易为一次降雨所充满。
因此，当流域发生大雨后，土壤含水量可以达到流域蓄水容量，降雨损失等于流域蓄水容量减去初始土壤含水量，降雨量扣除损失量即为径流量。
这种产流方式称为蓄满产流，方程式表达如下：
R=P-（Wm-Wo）（4-21）
但是，式4-21只适用于包气带各点蓄水容量相同的流域，或用于雨后全流域蓄满的情况。
在实际情况下，流域内各处包气带厚度和性质不同，蓄水容量是有差别的。
因此，在一次降雨过程中，当全流域未蓄满之前，流域部分面积包气带的缺水量已经得到满足并开始产生径流，这称之为部分产流。
随降雨继续，蓄满产流面积逐渐增加，最后达到全流域蓄满产流，称之为全面产流。
在湿润地区，一次洪水的径流深主要是与本次降雨量、降雨开始时的土壤含水量密切相关。
因此，可以根据流域历次降雨量、径流深、雨前土壤含水量，按蓄满产流模式进行分析，建立流域降雨与径流之间的定量关系。
二、降雨径流相关图
（1）降雨径流相关图的编制
根据流域多次实测降雨量P（雨期蒸发量可直接从雨量中扣除）、径流深R、雨前土壤含水量Wo,以Wo为中间变量建立P~Wo~R关系图，即流域降雨径流相关图，见图牛8所示。
当流域降雨量较大时，雨后土壤含水量可以达到流域蓄水容量，故P-Wo-R关系的右上部应是一组等距离的45。
直线，直线方程满足公式4-21o当流域雨前土壤含水量和降雨量较小时，流域部分面积蓄满产流，不满足全流域蓄满产流方程，在P-Wo-R关系线的下部表现为一组向下凹的曲线交汇于坐标轴的0点，见图牛8。
如果点绘在降雨径流相关图上P、R、Wo点据规律不明显，无法绘制出符合上述要求的P-Wo-R关系线，在P、R资料可鼎的前提下，则有可能是Wo的计算结果不合理，需要分析影响Wo计算值的参数。
一般说来，必是一个敬感性不强的参数，而流域蒸散发量对Wo影响比较显著。
因此，关键是对式4J0中的蒸发折算系数3的合理分析和取用，或调整流域蒸发计算模式。
当实测P、R、Wo点据较少时，也可以点绘P+Wq-R相关图，见图4-9所示。
此时，P+Wq^R关系线的上部是满足式4-21的45。
直线，P+Wo-R关系线的下部为向下凹的曲线交汇于坐标轴的0点。
在流域全面产流时，按P-Wo-R关系图或P+Wo-R相关图的查算结果相同；但在流域部分产流时，按P-W0-R关系图的查算结果的精度要高于P+Wo^R相关图。
图4-8严相关图图4-9P+Wo-R相关图
半流域径流资料不充分或分析困难时，可以采用前期影响雨量&代替14/0编制流域降雨径流相关图。
（二）降雨径流相关图的应用
降雨径流相关图、土壤含水量计算模式及相应参数构成了流域产流方案，据此可以进行流域产流计算。
依据产流方案，先由流域前期实测雨量、蒸发、径流资料推求本次雨前土壤含水量Wo,然后山本次降雨的时段雨量过程，查降雨径流相关图上相应于
图4-10由严W（TR相关图査算时段径流深
Wo的关系曲线，便可推求得本次降雨所形成的径流总量及逐时段径流深。
【例44】已知某流域一次降雨的逐时段雨量,见表44的第1、2栏，且计
算得雨前土壤含水量Wo=58mm,请根据Wo^R相关图（图牛10）查算降雨该
次所形成的逐时段径流深。
表4・1由严相关图査算时段径流深
■
J
P;
IP
1R
R；
（At=3h）
（mm）
（mm）
（mm）
（mm）
（1）
（2）
（3）
（4）
（5）
1
50
50
18
18
2
30
80
38
20
3
25
105
63
25
4
25
130
88
25
1.将表4-1第2栏时段降雨量转换为各时段末累积雨量1P,列第3栏;
2•在旷W0^R内插出Wo=58mm的P^R线，见图4-10；
3•由各时段末IP值查图4-10中Wo=58mm的LR线,得各时段末累积径
流深审，见表4-1第4栏；
4.
三.蓄满产流模型
将XR错开时段相减得岀各时段降雨所产生的径流深，见表4-1第5栏。
流域部分产流的现象主要是因为流域各处蓄水容量不同所致。
如果将流域内各点蓄水容量从小到大排列，最大值为"mm,计算大于某一的面积占流域面积的比重血则可绘出必‘〜Q关系曲线，称之为流域蓄水容量曲线，如图411所示。
山于流域蓄水容量在流域内的实际分布是很复杂的，要想用直接测定的办
图441流域蓄水容量曲线
法来建立蓄水容量曲线是困难的。
通常
的做法是通过实测的降雨径流资料来选配线型，间接确定蓄水容量曲线。
多数地区经验表明，流域蓄水容量曲线是一条单增曲线，可用B次抛物线来表示：
z、B
W'
a=l-1（4-22）
W神丿
式中B——反映流域内蓄水容量空间分布不均匀性的参数，取值一般为~。
W'mm流域内最大的点蓄水容量，取值一般为80^150mmo
蓄水容量曲线以下包围的面积（见图4-11）就是流域蓄水容量
%严J1-而竺dW'm=（4-23）
00\tmnJI+D
降雨初始的土壤含水量Wo采用递推公式牛9推求，流域蒸发可以根据不同要求选采用一层、二层或三层蒸发公式计算。
对应于Wo,流域土壤含水量已经达到蓄水容量的面积为血，相应于血的最大点蓄水容量为4见图442
（a）,Wo与A的关系为:
0o\
i-2L,
叭,”
R
）W'
IdW\n=
w1+3
1-1-—
、1+B
mm
将式牛23代入上式，整理后得：
卜宀）
图442流域蓄水容量曲线及部分产流
如果流域降雨量为C当A+P（b）,山P减去降雨损失AW得出产流量
“+P
/?
=—△"=—『（1一讪"為=/>+%一叱“
如果A+P>\A/mm流域为全面产流状态，山式4・22计算产流量。
综上所述，已知流域降雨量P和初始土壤含水量必时，蓄满产流模型的产流计算公式归纳为：
（4-24）
（4-25）
（4-26）
r=“+%-巴,+叱』_（A+P）/W\J+RA+P[P+Wo—巴,A+P>W'inm
\A/mm=（1+B）Wm
产流计算公式、土壤含水量递推公式以及流域蒸发公式，构成蓄满产流模型的主体。
蓄满产流模型不仅可以计算次降雨径流深，还能够进行连续演算，模拟多年径流过程。
因此，可以依据流域多年雨量、蒸发量及流量资料，分析率定模型的参数&必和乩
4.水源划分
按照蓄满产流概念，土壤含水量达到蓄水容量的面积称为产流面积，只有这部分面积上的降雨才能产生径流，其中的一部分按稳定下渗率下渗，形成地下径流，超过稳定下渗率的部分为地面径流。
由图牛12（b）可知，根据流域时段降雨量P及所产生的净雨量儿可以得岀产流面积比
（4-27）
a=h/P
根据稳定下渗率斤和产流面积比e就可以将各时段净雨h划分为地面净
雨忌和地下净雨Rg两部分:
h>c^fah（4-28）
（4-29）
hs=h■hg
流域稳定下渗率分析可以根据雨后流域蓄满的降雨径流资料分析推求。
首先按本章第二节介绍的方法得出P、R及Rg，雨前土壤含水量Wo=P-R,然后进行产流计算得出径流过程，根据式4-27和式4-28釆用试算法就可以求得稳定下渗率农。
【例牛2】已知某流域一次降雨过程及地下径流总量Rg=,并已经推求出时
段净雨量，见表牛2笫1~3栏，试推求稳定下渗率九。
表4・2稳定下渗率九的推求
t
P
（mm）
h
（mm）
a
A=h
A=h
afc't
（mm）
hg
（mm）
a/cAt
（mm）
hg
（mm）
（1）
（2）
（3）
（4）
（5）
（6）
（7）
（8）
2-8
8-14
14-20
20-2
2-8
8-14
1
计算步骤：
1.按式4-27计算产流面积a=h/P,见表4-2第4栏；
2.假定启，计算ot/血，见表4・2第5栏;
3.依据第5栏数值，按式4-28计算各时段為，见第6栏，求和得地下径流总量为，计算值显著大于实际值；
4.经分析后重新假定于尸，计算曲,见表4-2第7栏；
5.依据第7栏数值，按式4-28计算各时段加，求和得地下径流总量为，见表4-2第8栏，计算值与实际值相近。
最终，计算出本次洪水的农訪。
第四节超渗产流计算
一、超渗产流模式
在干旱和半干旱地区，降雨量小，地下水埋藏很深，包气带可达儿十米甚至上白米，降雨过程中下渗的水量不易使整个包气带达到田间持水量，一般不产生地下径流，且只有当降雨强度大于下渗强度时才产生地面径流，这种产流方式称为超渗产流。
在超渗产流地区，影响产流过程的关键是土壤下渗率的变化规律，这可用下渗能力曲线来表达。
如第二章所述，下渗能力曲线是从土壤完全干燥开始，在充分供水条件下的土壤下渗能力过程，见图4-13所示。
土壤下渗过程大体可分为初渗、不稳定下渗和稳定下渗三个阶段。
在初渗阶段，下渗水分主要在土壤分子力的作用下被土壤吸收，加之包气带表层土壤比较疏松，下渗率很大；随着下渗水量增加，进入不稳定下渗阶段，下渗水分主要受毛管力和重力的作用，下渗率随着土壤含水量的增加而减少；随着下渗水量的锋面向土壤下层延伸，土壤密度变大，下渗率随之递减并趋于稳定，也称为稳定下渗率。
与蓄满产流相比，超渗产流的影响因素更为复杂，对讣算资料的要求较高，产流计算成果的精度也相对较差。
因此，必须对干旱地区下渗特性及主要影响要素进行深入分析，充分利用各种资料条件，制定合理的超渗产流计算方案。
图4-13下渗能力曲线图4-14下渗曲线法
二、下渗曲线法
按照超渗产流模式，判别降雨是否产流的标准是雨强i是否超过下渗强度fo因此，用实测的雨强过程/〜r扣除实际下渗过程f〜r,就可得产流量过程/?
〜r,如图4-14中阴影部分。
这种产流计算方法称为下渗曲线法。
在实际降雨径流过程中，流域初始土壤含水量一般不等于0,降雨强度并非持续大于下渗强度，不能直接采用流域下渗能力曲线推求各时段的实际下渗率。
如果将下渗能力曲线转换为下渗能力与土壤含水量的关系曲线，就可以通过土壤含水量推求各时段下渗强度了。
在第二章已经提到，流域下渗能力曲线常用霍顿下渗公式来表达，即：
/（0=（厶一£比"+£（4-30）
根据霍顿下渗公式可以推求累积下渗量曲线
）=j=fct+扣。
一£）一寺（/0-£）严（4-31）
0PP
F（t）为累积下渗量，这部分水量完全被包气带土壤吸收，也就是r时刻流域的土壤含水量，因此有
呵）"+討"）-討rw
连立求解式牛30和式4-32,消除时间变量r,可以得出下渗强度与土壤含水量的关系曲线广W，见图4-15of-W反映了土壤含水量变化对下渗强度的影响。
根据雨前土壤含水量Wo,就可以山降雨过程采用广W关系曲线逐时段进行产流计算，步骤如下：
1.
（4-32）
图445广W关系曲线
从降雨第一时段起，由时段初始土壤含水量M查广W曲线，得到相应的下渗率人，如果时段不长，可以近似代表时段平均下渗率。
2.根据人及时段雨强亦按超渗产流模式计算净雨量加，讣算公式为：
（4-33）
i>f
h={
[0i3.根据水量平衡公式，计算下时段初始土壤含水量：
（4-34）
Wk+1=Wk+Pk-hk
4.重复步骤1~3就可以山降雨过程计算出逐时段的产流量。
采用下渗曲线法进行产流计算时，应该注意到降雨强度时空分布的不均匀性对产流的影响，且流域不同地点的下渗特点也是存在差别的。
因此，为了提高计算精度，降雨时段长度不宜大，常以分钟计，流域应按雨量站分布状况划分为较小的单元区域进行产流计算。
三、初损后损法
采用下渗曲线法进行产流讣算，必须知道计算区域的下渗能力曲线，这需要很多径流资料或实地试验才能获得，在实际工作中往往难以实现。
初损后损法是下渗曲线法的一种简化方法，它把实际的下渗过程简化为初损和后损两个阶段。
产流以前的总损失水量称为初损，以流域平均水深表示；后损主要是流域产流以后的下渗损失，以平均下渗率表示。
一次降雨所形成的径流深可用下式表示：
R=P-lo-JIr-Po（4-35）
式中P一一次降雨量，mm；
Io初损，mm；
7—一平均后渗率，mm/h：
tR
升级会员 