大伙房流域降雨径流模型

1、.大伙房流域降雨径流模型大伙房流域降雨径流预报模型又简称“dhf”模型，该模型于1973 年由辽宁省大伙房水库管理局刘爱杰、王本德等人提出，至今已使用30 余年，为水库洪水调度做出了很大贡献。“ dhf”模型是适用于我国湿润地区的超渗产流模型，目前已在辽宁省多个水库的水情自动测报系统中使用，效果较好。建立在“ dhf”模型基础上的降雨径流预报方案，在大伙房流域经过调试和精度验证后进行使用，在使用中平均精度令人满意。尤其在“ 957”特大洪水调度中，发挥了显著作用，准确预报出了第一非常溢洪道溢流时间和水库最高库水位，为省防制定调度决策提供了科学依据，使水库工程发挥了强大的调蓄作用，最大限度地配合

2、了下游抢险，共减免下游直接经济损失74.89 亿元。“ dhf”模型由两部分组成，一是八参数超渗产流计算模型，引用双层入渗曲线进行扣损计算，并以抛物线描述表层蓄水量和下层渗率的分布状况；二是八参数变强度、变速度的经验单位线汇流计算模型，参数随降雨分布而变，采用“前期影响净雨”描述汇流速度的变化。这是一个集总的概念模型，模型的参数多半在满足其物理意义的前提下确定，只有6 个需要优选法选定或试错法确定。1 大伙房模型概化流程流域下垫面分为表层、下层和地下水蓄存三部分，计算流程如图10-1 所示。2 大伙房模型产流计算产流模型将下垫面分为表层，下层和深层三部分。表层土壤中的张力水蓄量与植物截流、填洼

3、储存合称表层蓄水量sa ，其极值为表层蓄水容量 s0 ；下层土壤中的张力水蓄量称为下层蓄水量u a ，其极值为下层蓄水容量 u 0 ；地下水储水层的蓄水量以va 表示，其极值为地下水库蓄水容量v0 。;.ed=kc*er/c雨强 p蒸散发雨间蒸发pe =p-ed不透水 面积上pc=pe -y 0直 接 径流蒸 散 发sa/s0y =g*p表 层 存储+sa地 面 壤 中r=pc-s0rc=r-yu流蒸 散 发ua/u0 下层存储rl =rc(u0-ua)蒸 散 发va/v0地下水存储地下径流y l=kw*r l深层地下水存储总出流图 10-1大伙房产流模型概化流程图其中：s0- 表层蓄水容量；

4、u0- 下层蓄水容量；d0- 下层蓄水容量（ u0）与地下水库蓄水容量（v0）之和；g- 不透水面积占全流域面积之比值；kc- 流域蒸散发能力与大水体蒸发量的比值；a- 表层蓄水容量曲线形状参数；b- 下渗率抛物线分配曲线形状参数；k- 下层下渗曲线曲率；2.1不透水面积上的径流计算根据上面的流程图，不透水面积上直接产生径流，通过下式来计算：ycgpepepedpkc er / c（1）（2）式中 p 为时段降雨量，kc 是流域蒸散发能力与er 大水体蒸发量的比值，c 为日降雨观测的时段数，g 为流域不透水面积占总流域面积的比值。;.2.2地面壤中流的计算地面壤中流通过下式计算yu rf（3）

5、2.2.1 表层扣损 r的计算当smpcas0时，r pcsa s0s0 (1 (smpc ) /(as0 ) a（4）当 smpcas0 时， rpc (s0sa )（ 5）式中 sa 是表层流域平均蓄水量，sm是相应于 sa 的点蓄水量； a 为曲线形状参数。sm和 sa 由下式计算得到：smas0 1 (1 sa / s0 ) 1/ a（6）sasm / a 1 (1 sm / sm )（7）式中 sm 为表层点最大蓄水量，当smsm 时，则 sa s0 ，并有smas0（8）其中 pc 通过下式来计算：pcp edyc（9）2.2.2下渗计算时段平均总下渗率f与da的关系 :fz1(

6、d0da )(10)时段平均下层下渗率fu与da的关系 :f uz2(u0ua )(11)其中 z1 1 e k2 vtu0 / d0 ; z21e k2 vt， k1 和 k2 的关系如下 :k1k2u0/ d0(12)k1 和 k2 分别是反映总下渗和下层下渗曲线曲率系数；da 是 ua 的函数，二者的关系式如下 :dad01 (1 ua / u0 )u0 / d0(13)d0 是下层蓄水容量u0 与地下水库蓄水容量v0 之和。2.2.3径流深计算公式;.总径流 y的计算公式：z2un rby r z2 (uau0 ) z2u01z2 bu0(14)当 rz2unz2 bu0 时，得y r

7、 z2 (u0ua ) r fu(15)地面壤中流yu 的计算公式：当 z1dnrz1bd0 时，yur z1 d0 (1 ua / u0 )u0 / d0z1d1 (1z1 dnr ) bz1bd0(16)当z1dnrz1bd0 时，yur z1 d0(1 ua / u0 )u0 / d0(17)d n 和 un 通过下列两式计算dnbd0 1 (1 ua / u0 )u0 /(bd 0 )(18)unbu01 (1 ua / u0 )1/ b(19)地下径流 yl 的计算：ylkw rl(20)其中 kw 为 yl 与 rl 之比值，在干旱地区rl 全部渗入深层，不形成地下径流，则kw0

8、；在湿润地区 rl 全部形成地下径流流入河道，所以kw1。 rl 为地下水库下渗强度，为总径流与壤中流之差，即：rly yu(21)2.3大伙房模型蒸散发计算（ 1）当降雨量超过流域蒸发能力时，即净渗雨强pc0 ，产生下渗，土壤蓄水量增加，此时：表层蓄水量计算公式sa(t 1)s0 1 (1 (smt pct ) / as0 )a(22)式中 t 为单位日，以8:00 为日分界；sm 以式（ 6）计算当 smtpct as0 时，上式可写成;.sa (t 1)s0(23)依照流程图的下渗程序，逐时段计算式，sa (i 1)saipciri(24)式中 i 为时段序；pc 、 r 分别用（ 9）

9、、（ 4）、（ 5）式进行计算。下层土壤蓄水计算公式据流程图，逐时段计算ua 的公式ua (i 1) uai riyi(25)若逐日计算，可用下式计算若 z2 untrtz2 bu0 时ua(t 1)uatz2(u0uat )z2u0 (1 z2untrt )bz2 bu0(26)若 z2untrtz2 bu0 时ua (t1)uatz2 (u0 uat )(27)式中的un与z2分别用式（19）及z21e k2 vt求出。（ 2）当日流域平均降雨量小于流域蒸发能力，即ped 时，土壤蓄水量将减少进行蒸散发计算表层蒸散发公式假定： a. 表层蒸发等于流域平均蒸散发能力；b. 每次下渗后的初次蒸

10、散发皆从地表开始。eu 为表层可蒸发量，计算公式为：es (1eb ,t 1 )a(1eb ,t ) au ,t0as0as0(28)式中：ec, teu, tptt1eb ,t 1ecneaeb,teb ,t 1ec,t(29)因此， t1日表层蓄水量计算式为sa ,t 1sa ,teu, t(30);.下层蓄水量 ua 计算公式假定： a. 表层蓄水量蒸发尽后，该点的下层开始蒸发；b. 下层蒸发量 el 与下层蓄水量 ua 成比例关系。因此，el(eceu )ua/ u0(31)当 ( s eu )0 时，则el( ecsa )ua / u0(32)这样， t1日下层土壤蓄水量以下式计算u

11、a,t 1ua,tel(33)2.4大伙房模型汇流计算大伙房汇流模型是根据多年预报经验，总结出的8参数变强度、变汇流速度的经验单位线，其线型公式为cctk 3 e aa( ?t / n?tm )ddsintntmntm式中t汇流曲线，是tm 与 t 的函数；tm汇流曲线底宽；t 汇流时间（时刻）；k3 汇流曲线的比例系数；n 地下径流与地面壤中流汇流曲线底宽的比例系数；aa、 dd、 cc反映汇流曲线形状的参数；地表径流与地下径流具有不同的汇流曲线形状，可用形状参数区分；不同的降雨分布地表径流也具有不同的汇流曲线，其参数也不同。tm 是时段净雨（径流）与前期影响净雨（已产生的径流）的函数，模型

12、中假定：tmlb ( yayt )k 0ya, t 1k ( ya ,tyt ,t )式中lb 特征河长lbl / b0 ；l 流域最大河长；;.k 0 拟合线的斜率；ya 前期影响净雨量，表征河槽蓄水状况，制约汇流速度等；k 前期净雨影响程度衰减系数，反映已产生的径流在河槽中消退速度；2.5大伙房模型产流参数优选大伙房模型产流参数的最优化问题是多目标的有约束的非线性规划问题。目标函数取为n次洪水产流的拟合精度和n次洪水计算产流与实际产流误差总和的线性加权，各取0.5 的权重进行求解。将问题变为：max w1* acu 1 w 2 * tytc s.t. xi minxixi max式中acu

13、1-n次洪水产流的拟和精度tytc-n次洪水计算产流与实际产流误差总和acu 1n合格 / ntytc(totaltytc) / total其中 n为洪水的次数，total为 n次洪水的实际净雨的总和，tytc为 n次洪水的预报净雨的总和，n 合格为拟合的合格洪水的次数。若某次洪水满足以下条件之一即为合格：计算产流与实际产流的绝对误差小于3mm；计算产流与实际产流的相对误差小于20%且绝对误差小于20mm；计算时判断合格的绝对误差和相对误差的标准可以根据具体情况来设定；本程序采用可变容差法求解问题，把多个约束求极小值问题变为一个单约束求极小值问题。设所求解的非线性规划问题为：minf ( x

14、)xenhi ( x )0(i1 , 2,.,m )s .t .0(jm1 ,., p )gi ( x )本算法将上述两问题化为单约束问题进行求解:minf ( x )xe ns . t .( k)t( x )0式中 t(x) 是约束破坏的估计量，表达式为：mp1/ 2t (x)hi 2 ( x)u j g j 2 (x)i 1j m 1其中， uj 为 heaviside 算子。( k)是第 k步搜索中给出的关于可行性的可变容差准则值，表达式为：;.(k )min(k 1) , m1x(k )x(k )r1ir 2(0)2(m 1)t其中t-初始多面体的大小 ;m-等式约束数目 ;r=n-m

15、 ， f ( x) 的自由度 ;xi ( k )n- e 多面体中第 k次搜索时的第 i 个顶点值 ;xr 2(k )xh-多面体（除）外的形心顶点 ;由于各参数的初值相差较大，而且变化范围也各不相同，所以将参数归一化，即初值都取为1.0 。根据物理意义将e6、 e7、 e8、 e9的变化范围取为 0.5-10.0，其它参数取为 0.5-5.0 。2.6大伙房模型汇流参数优选大伙房汇流模型需要优选的参数有b0、 k0、 k、n、 dd、cc，其中 dd和cc根据降雨的分布分为上游、中游和下游三套参数，所以需要优选的参数一共有10个。 k3和 aa可以通过优选出的参数带入下列两式计算求得：ntm

16、cct )k 3( ntm0 ) /e aa( t / ntm )ddsin(0ntmaacc / dd ( t0 ) /(ntm ) dd1 tan ( t0 ) / ntm大伙房汇流参数优选的基础是多次历史洪水，其最优化问题是多目标的有约束的非线性规划问题。目标函数取为 n次洪水过程的拟合精度、n次洪水计算洪峰与实际洪峰误差总和以及n次洪水的计算峰现时间和实际的峰现时间误差总和的线性加权，各取0.4 、 0.3 、 0.3 的权重进行求解。将问题变为：minw1qe 1w2 qe 2w3 qe 3s.t .xi minxixi max式中qe1-n次洪水汇流的拟和精度;qe2-n次洪水计算

17、峰值与实际峰值误差总和;qe3-n次洪水计算峰现时间与实际峰现时间误差总和;本程序采用可变容差法求解问题，把多个约束求极小值问题变为一个单约束求极小值问题. 设所求解的非线性规划问题为:minf ( x )xnehi ( x )0(i1, 2,., m )s .t .0(jm 1,., p )gi ( x )本算法将上述两问题化为下式进行求解:;.minf ( x ) xe ns.t .( k )t ( x )0式中 t(x) 是约束破坏的估计量，表达式为:m1/ 2pt (x)hi 2 ( x)u j g j 2 (x)i 1j m 1其中， uj 为 heaviside 算子。( k)是第

18、 k步搜索中给出的关于可行性的可变容差准则值，表达式为：( k )min( k1) , m1xi ( k ) xr 2 ( k )(0)r12( m1)t其中t-初始多面体的大小;m-等式约束数目 ;r=n-m ， f ( x) 的自由度 ;xi ( k ) -en 多面体中第 k次搜索时的第 i 个顶点值 ;(k )xr 2- 多面体（除xh ）外的形心顶点 ;2.7大伙房产流模型参数为提高产流预报精度，按降雨开始时的起始流量大小进行分类。分类预报参数如表10-1 表 10-2 所示。表 10-1降雨起始流量100m3/s 的大伙房产流模型参数gabks0u0d00.054.1291.1810.2845100130续表 10-1流域蒸发能力 ed月份6789日雨量 5mm3.74.34.33.75mm 日雨量 1.91.91.91.930mm日雨量 30mm1.01.01.01