2018/11/1
9
（2）超蓄产流模型的结构 a）点模型 以含气层缺水量为控制条件，就流域中某点而言：
蓄满前: P E WW2 WW1 蓄满后: P E R 式中: P : 时段降雨量 E : 时段蒸散发量 R : 时段产流量
(6 - 1)
WW1 , WW2 : 时段初末的土壤含水量
21.11
30.51 24.27 84.62
0.0427
0.0632 0.0495 0.2091
0.50
0 7.46 17.61
0.50
0 2.22 5.84
0
0 5.23 11.76
29.57
23.68 76.57 113.20
23
24
20.27
-2.79
138.85
156.00
0.4844
1.0000
c）流域产流计算 P－E>0时，产流，否则不产流 ，产流时：
P E A WWMM 时 : R P E ( WM W ) P E A WWMM 时 : R P - E - (WM - W) - WM1 - (P - E A)/WWMM

1 B

13
2018/11/1
10
式(6 1)中 , R RS RG , 即 RG FC RS R - RG P - E - FC FC : 时段稳定下渗量
b）流域蓄水容量曲线（超蓄产流模型的核心）
(6 - 2)
WWM：流域蓄水容量 WWMM:流域最大蓄水容量
WM：流域平均蓄水容量
I
E
2018/11/1
18
2、用试算法求fc f RS i R i i f c t i F
fi RS RS i R i f c t i 1 1 1 F f R 又 i F PE 得: fc
n n n
R
1 n
n
i
RS fc
R
1 n
n
i
RS (忽略雨期蒸散) (8 - 11)
EU EL ED
WUM
WLM
C
出流过程
KE XE
径流 R
径流
R
2018/11/1
8
二、二水源新安江模型的微结构 （一）用超蓄产流（即“蓄满产流”）模型计算总径流 R、地表径流RS 及地下径流RG （1）超蓄产流模型概念 超蓄产流模型是目前我国湿润地区的主要产流模型。 “蓄满”，指含气层的土壤含水量达到田间持水量，而非土壤完全 饱和； “超蓄产流”指土壤达到田间持水量以前不产流，所有降雨都被土 壤吸收，成为薄膜水和张力水；而在土壤达到田间持水量以后，所 有降雨（除去同期蒸发）都产流。这时土壤的下渗能力为稳定下渗 率，稳定下渗量FC补充地下水，形成地下径流，而超渗的部分则形 成地表径流。 与“超渗产流”模型的区别： “超蓄产流”模型先计算R，在分成RS、RG； “超渗产流”模型先计算RS、RG，再合成R。
EL C ( EP EU ) ，ED 0
WL C WLM
若
且 WL C ( EP EU ) ，则
EL WL ， ED C ( EP EU ) WL
2018/11/1
24
If w(1) + p(i) > ep(i) Then e(1) = ep(i) e(2) = 0 e(3) = 0 Else e(1) = w(1) + p(i) e(2) = (ep(i) - e(1)) * w(2) / wm(2) If w(2) <= c * wm(2) Then e(2) = c * (ep(i) - e(1)) e(3) = 0 If w(2) >= c * (ep(i) - e(1)) Then e(2) = c * (ep(i) - e(1)) e(3) = 0 Else e(2) = w(2) e(3) = c * (ep(i) - e(1) - e(2)) End If End If End If w(1) = w(1) + p(i) - r - e(1) w(2) = w(2) - e(2) w(3) = w(3) - e(3) If w(1) > wm(1) Then w(2) = w(1) - wm(1) + w(2) w(1) = wm(1) If w(2) > wm(2) Then w(3) = w(3) + w(2) - wm(2) w(2) = wm(2) End If End If
(6 - 7)
产流计算特点：雨强对产量无影响，产流量取决于P－E与W。
2018/11/1
模型参数：WM与B WM：流域干燥时的缺水量，代表 流域干旱情况，气候因素； B：蓄水容量在流域上的分布不均 匀性，B＝0时分布均匀，愈大愈不均匀， 决定于地形、地质条件。
d）地面、地下径流的划分（分水源） 产流面积变化，则： P E FC时: R RG FC(f/F) FC P-E RS R - RG (6 - 8) P - E FC时: RS 0, RG R
(6 - 3)
大量资料表明，WWM～f/F有如下关系：
f 1 (1 F f 或 1 (1 F
2018/11/1
WWM B ) WWMM WWM B ) WWMM
(6 - 4)
12
则：
WM

1
0
WWMd (f / F)
WWMM 1 B
(6 - 5)
对纵坐标积分 :
A f WWM W (1 )dWWM (1 )dWWM 0 0 F WWMM 1 W 1 B A WWMM 1 - (1 ) (6 - 6) WM A
2018/11/1 21
W=WU+WL+WD
E=EU+EL+ED
模型中流域蓄水容量WM和流域蓄水量W都是上层, 下层和深层之和, 即: WM WUM WLM WDM W WU WL WD WUM, WLM, WDM - - - -上层, 下层, 深层土壤蓄水容量 WU, WL, WD - - - -上层, 下层, 深层土壤蓄水量 当WU EM , 则 当WU EM , 则 EU EM EL 0 EU WU (6 - 14) (6 - 13)
2018/11/1 14
例6-1: 超蓄产流模型产流量计算示例 WM＝120mm，B＝0.3，FC＝18mm/d
年.月.日
1978.7.17 18
P－E
4.22
A
16.89
f/F
0.0338
R
0.16
RG
0.16
RS
0
W
16.61 20.67
19
20 21 22
9.40
-5.98 60.35 54.24
单元流域面积要适中，使得在每块面积上降雨比较均匀，并有一定数 目的雨量站；(泰森多边形) 其次，尽可能是单元流域与自然流域相一致；若流域中有大中型水库， 则水库以上的集水面积即可作为一个单元流域。
2018/11/1
4
Rainfall Averaging Methods Thiessen Polygons
2018/11/1(6 - )(6 - 10)16
（二）稳定下渗率fc的推求
1、求一场洪水的RS、R、RG （1）据上图求RS （2）根据图求R （3）求RG=R-RS （4）fc=RG/T T为净雨时间
2018/11/1 17
Q(m3/s) B
N
本次降雨形成的径流过程
H
C 直接径流 B’ A C’ 地下径流 F G D’ t(h) D
WL EL (EM - EU) WLM EU, EL - - - -上层, 下层的时段蒸发量。
(6 - 15)
（ 6－ 14）、（ 6－ 15）式表明，上层有足够水分是，EU 等于蒸散发能力，下层EL ＝0； 当上层水分不足时，把上层水分全蒸掉，下层蒸发量EL 等于上层剩余蒸发能力 （EM －EU ）与下层含水比WL / WLM 之乘积。
2018/11/1
22
（ 6－ 15）中，只用到 否则，取
EL C（EM －EU ） EL ＝C（EM －EU ） EL ＝WL
但是，若EL C（EM －EU ） ,同时WL C（EM －EU ） , 则取 此时 这时才发生深层蒸发。 ED 深层的时段蒸发量 C 与深层蒸散发有关的系数 ED ＝C（EM －EU ）－EL
13.47
0
11.96
0
1.51
0
120.00
117.21
2018/11/1
15
尤其半湿润地区需要考虑
IMP:不透水面积参数（新安江模型新增参数），流域不透水面 积占总面积的百分比，增加后，需修改（6－5），（6－ 8）式，其它都不变。
WWMM (6 5) WM (1 IMP ) 1 B (6 - 8) P - E FC时 R RG FC IMP P - E RS R RG
马斯京根法或滞 后演算法
KE XE 或 L
7
2018/11/1
2、二水源新安江模型
一、二水源新安江模型基本结构(状态变量（模型参数）)
蒸散发E 降雨P 蒸发皿蒸发EI 透水面积 土壤湿度 W 上层 WU 下层 WL 深层 WD WM B FC 地下径流 RG KKG 地下径 流过程 单元流域 不透水面积 IMP 地面径流 RS UH 地面径 流过程
2018/11/1
11
利用流域蓄水容量曲线计算产流量（右图）：