大伙房流域降雨径流模型(共10页)

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

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

3、分。表层土壤中的张力水蓄量与植物截流、填洼储存合称表层蓄水量，其极值为表层蓄水容量；下层土壤中的张力水蓄量称为下层蓄水量，其极值为下层蓄水容量；地下水储水层的蓄水量以表示，其极值为地下水库蓄水容量。图10-1 大伙房产(fngchn)流模型概化流程图其中(qzhng)：S0-表层蓄水(x shu)容量；U0-下层蓄水容量；D0-下层蓄水容量（U0）与地下水库蓄水容量（V0）之和；g-不透水面积占全流域面积之比值；Kc-流域蒸散发能力与大水体蒸发量的比值；A-表层蓄水容量曲线形状参数；B-下渗率抛物线分配曲线形状参数；K-下层下渗曲线曲率；2.1 不透水面积上的径流计算 根据上面的流程图，不透水

4、面积上直接产生径流，通过下式来计算： （1） （2）式中P为时段降雨量，KC是流域蒸散发能力与ER大水体蒸发量的比值，C为日降雨观测(gunc)的时段数，g为流域不透水面积占总流域面积的比值。2.2 地面(dmin)壤中流的计算地面壤中流通过(tnggu)下式计算 （3）2.2.1 表层扣损R的计算当时， （4）当时， （5）式中Sa是表层流域平均蓄水量，Sm是相应于Sa的点蓄水量；a为曲线形状参数。Sm和Sa由下式计算得到： （6） （7）式中为表层点最大蓄水量，当时，则，并有 （8）其中通过下式来计算： （9）2.2.2 下渗计算时段平均总下渗率与的关系: (10)时段平均下层下渗率与的关

5、系: (11)其中，和的关系如下: (12) 和分别是反映总下渗和下层下渗曲线曲率系数；是的函数，二者的关系式如下: (13)是下层蓄水(x shu)容量与地下水库(d xi shu k)蓄水容量之和。2.2.3 径流(jngli)深计算公式总径流y的计算公式： (14)当时，得 (15)地面壤中流的计算公式：当时， (16)当时， (17)和通过下列两式计算 (18) (19)地下径流的计算： (20)其中为与之比值，在干旱地区全部渗入深层，不形成地下径流，则；在湿润地区全部形成地下径流流入河道，所以。为地下水库下渗强度，为总径流与壤中流之差，即： (21)2.3 大伙房模型蒸散发计算（1）

6、当降雨量超过流域蒸发能力时，即净渗雨强，产生下渗，土壤蓄水量增加，此时：表层蓄水量计算公式 (22)式中t为单位(dnwi)日，以8:00为日分界；以式（6）计算(j sun)当时，上式可写成 (23)依照流程图的下渗程序(chngx)，逐时段计算式， (24)式中i为时段序；、分别用（9）、（4）、（5）式进行计算。下层土壤蓄水计算公式 据流程图，逐时段计算的公式 (25)若逐日计算，可用下式计算若时 (26)若时 (27)式中的与分别用式（19）及求出。（2）当日流域平均降雨量小于流域蒸发能力，即时，土壤蓄水量将减少进行蒸散发计算表层蒸散发公式假定：a.表层蒸发等于流域平均蒸散发能力；b.

7、每次下渗后的初次蒸散发皆从地表开始。为表层可蒸发量，计算公式为： (28)式中： (29)因此(ync)，日表层(biocng)蓄水量计算式为 (30)下层(xicng)蓄水量计算公式假定：a.表层蓄水量蒸发尽后，该点的下层开始蒸发；b.下层蒸发量与下层蓄水量成比例关系。因此， (31)当时，则 (32)这样，日下层土壤蓄水量以下式计算 (33)2.4 大伙房模型汇流计算大伙房汇流模型是根据多年预报经验，总结出的8参数变强度、变汇流速度的经验单位线，其线型公式为式中 汇流曲线，是与t的函数； 汇流曲线底宽；t 汇流时间（时刻）；K3 汇流曲线的比例系数；n 地下径流与地面壤中流汇流曲线底宽的比

8、例系数；AA、DD、CC反映汇流曲线形状的参数；地表径流与地下径流具有不同的汇流曲线形状，可用形状参数区分；不同的降雨分布地表径流也具有不同的汇流曲线，其参数也不同。是时段净雨（径流）与前期影响净雨（已产生的径流）的函数(hnsh)，模型中假定：式中 特征(tzhng)河长 ；流域(liy)最大河长；拟合线的斜率；前期影响净雨量，表征河槽蓄水状况，制约汇流速度等；前期净雨影响程度衰减系数，反映已产生的径流在河槽中消退速度；2.5 大伙房模型产流参数优选大伙房模型产流参数的最优化问题是多目标的有约束的非线性规划问题。目标函数取为N次洪水产流的拟合精度和N次洪水计算产流与实际产流误差总和的线性加权

9、，各取0.5的权重进行求解。将问题变为：式中 ACU1-N次洪水产流的拟和精度 TYTC-N次洪水计算产流与实际产流误差总和 其中N为洪水的次数，TOTAL为N次洪水的实际净雨的总和，TYTC为N次洪水的预报净雨的总和， N合格为拟合的合格洪水的次数。若某次洪水满足以下条件之一即为合格：计算产流与实际产流的绝对误差小于3mm；计算产流与实际产流的相对误差小于20%且绝对误差小于20mm；计算时判断合格的绝对误差和相对误差的标准可以根据具体情况来设定；本程序采用可变容差法求解问题，把多个约束求极小值问题变为一个单约束求极小值问题。设所求解的非线性规划问题为：本算法将上述两问题化为单约束问题进行(

10、jnxng)求解:式中T(x)是约束(yush)破坏的估计量，表达式为：其中(qzhng)，Uj为Heaviside算子。是第K步搜索中给出的关于可行性的可变容差准则值，表达式为：其中 t-初始多面体的大小; m-等式约束数目;r=n-m，的自由度;-多面体中第K次搜索时的第i个顶点值;-多面体（除）外的形心顶点;由于各参数的初值相差较大，而且变化范围也各不相同，所以将参数归一化，即初值都取为1.0。根据物理意义将E6、E7、E8、E9的变化范围取为0.5-10.0，其它参数取为0.5-5.0。2.6 大伙房模型汇流参数优选大伙房汇流模型需要优选的参数有B0、K0、K、n、DD、CC，其中DD

11、和CC根据降雨的分布分为上游、中游和下游三套参数，所以需要优选的参数一共有10个。K3和AA可以通过优选出的参数带入下列两式计算求得：大伙房汇流参数优选的基础是多次历史洪水，其最优化问题是多目标的有约束的非线性规划问题。目标函数取为N次洪水过程的拟合精度、N次洪水计算洪峰与实际洪峰误差总和以及N次洪水的计算峰现时间和实际的峰现时间误差总和的线性加权，各取0.4、0.3、0.3的权重进行求解。将问题变为：式中 QE1-N次洪水汇流(huli)的拟和精度; QE2-N次洪水(hngshu)计算峰值与实际峰值误差总和; QE3-N次洪水计算峰现时间与实际(shj)峰现时间误差总和;本程序采用可变容差

12、法求解问题，把多个约束求极小值问题变为一个单约束求极小值问题.设所求解的非线性规划问题为:本算法将上述两问题化为下式进行求解:式中T(x)是约束破坏的估计量，表达式为:其中，Uj为Heaviside算子。是第K步搜索中给出的关于可行性的可变容差准则值，表达式为：其中 t-初始多面体的大小; m-等式约束数目;r=n-m，的自由度;-多面体中第K次搜索时的第i个顶点值;-多面体（除）外的形心顶点;2.7 大伙房产流模型参数为提高产流预报精度，按降雨开始时的起始流量大小进行分类。分类预报参数如表10-1表10-2所示。表10-1 降雨起始流量100m3/s的大伙房产流模型参数gABKS0U0D00

13、.054.1291.1810.2845100130续表10-1 流域(liy)蒸发能力Ed月份6789日雨量 5mm3.74.34.33.75mm日雨量 30mm1.91.91.91.9日雨量30mm1.01.01.01.0表10-2 降雨起始(q sh)流量100m3/s的大伙房产流模型(mxng)参数gABKS0U0D00.054.1291.3710.46980100140续表10-2 流域蒸发能力Ed月份6789日雨量 5mm35.15.135mm日雨量 30mm2.32.32.32.3日雨量30mm1.41.41.41.42.8 大伙房汇流模型参数表10-3 大伙房模型汇流参数表KaLNB0K0CclDdl0.8661202.50.8250.430.350.77续表10-3 大伙房模型汇流参数表降雨中心CcDdCoe上游0.82.80.19中游1.01.2