新安江流域水文模型

第二章新安江流域水文模型

60年代初，河海大学(原华东水利学院)水文系赵人授等开始研究蓄满产流模型，配合一定的汇流计算，将模型应用于水文预报和水文设计。1973年，他们在对新安江水库做人库流量预报的工作中，把他们的经验归纳成一个完整的降雨径流流域模型——新安江模型。模型可用于湿润地区和半湿润地区的湿润季节径流模拟和计算。

最初的新安江模型为两水源模型，只能模拟地表径流和地下径流。80年代初期，模型研制者将萨克拉门托模型与水箱模型中，用线性水库函数划分水源的概念引入新安江模型，提出了三水源新安江模型，模型可以模拟地面径流、壤中流、地下径流。1984至1986年，又提出了四水源新安江模型，可以模拟地面径流、壤中流、快速地下径流和慢速地下径流。三水源新安江模型一般应用效果较好，但模拟地下水丰富地区的日径流过程精度不够理想。在新安江三模型中增加慢速地下水结构就成为四水源新安江模型。

当流域面积较小时，新安江模型采用集总模型，当面积较大时，采用分块模型。分块模型把流域分成许多块单元流域，对每个单元流域做产、汇计算，得到单元流域的出口流量过程。再进行出口以下的河道洪水演算，求得流域出口的流量过程。把每个单元流域的出流过程相加，就求得了流域出口的总出流过程。

划分单元流域的主要目的是处理降雨分布的不均匀性，因此单元流域应当大小适当，使得每块面积上的降雨分布比较均匀．并有一定数目的雨量站。其次尽可能使单元流域与自然流域相一致，以便于分析与处理问题，并便于利用已有的小流域水文资料。如果流域内有大中型水库，则水库以上的集水面积即应作为一个单元流域。因为各单元流域的产汇、流计算方法基本相同，以下只讨论一个单元流域的情况。

2.1新安江两水源模型

1．模型结构和参数

新安江两水源模型的产流子模型采用蓄满产流模型，蒸发计算采用三层蒸发计算模型。利用稳定下渗率FC将径流划分为地面径流和地下径流两种水源。地面径流采用单位线汇流，地下径流采用一次线性水库汇流。模型把流域面积划分为透水面积和不透水面积两部分，不透水面积上的降水在满足蒸发后将直接转化为地面径流。透水面积上将发生下渗，下渗的水量一部分存储于土壤层，后期耗于蒸发；满足了流域土壤蓄水容量后的下渗水量才能转化为径流。

不透水面积用参数IMP表示，它是用流域内不透水面积占全流域面积的百分比表示的。新安江模型的输出是流域出流过程t

Q～和流域蒸散发过程t

E～，输入则为时段降雨量P、蒸发皿观测蒸发量EI。

新安江两水源模型共有9个参数，一条单位线。

K——流域蒸发折算系数，是流域蒸散发能力与蒸发皿蒸发量之比；

C——深层蒸散发系数；

IMP——不透水面积占全流域面积的百分比，％；

W——流域平均蓄水容量(指张力水)，mm；

m

WUM——流域平均上层蓄水容量, mm；

WLM——流域平均下层蓄水容量, mm；

B——蓄水容量曲线指数；

mm/；

FC——稳定下渗率h

KKG——地下水消退系数；

UH——单位线。

模型结构如下图所示：

图0 新安江两水源模型结构示意图

图中方框内标注为状态变量，方框外标注的是模型参数。

2．模型参数的调试

当模型初定后就可上机调试，在调试的过程中，应注意各参数的径流响应特征，以及参数之间相依性的影响，以便调试能有目的地顺利进行。以下简介的调试技术可供参考。

流域蒸发折算系数K的调试，对一个具体流域来说，此参数完全靠优选。调试此参数主要考察年径流模拟误差，当年径流模拟误差达最小时，此值最优，但个别参数的最优并不表示模型总体最优。根据现有经验，年径流模拟误差控制在5％～8％左右就可以了。

W的调试，一般根据实测的降雨径流资料分析或经验选流域平均蓄水容量

m

取后，适当微调即可。经验选定后，此值调试范围不大，比较容易确定。通常在W可选80～150mm左右，江淮一带约在110mm左右,燕山东北东南方湿润地区

m

部地区，选150mm左右。

流域平均上层蓄水容量WUM，一般取10～20mm，流域平均下层蓄水容量WLM一般取60～90mm。

蓄水容量曲线指数的调试，B表示流域蓄水容量分布的不均匀性，当全流域蓄满后，这个参数就不起作用了，因此，应该选取流域没有达到蓄满的那些洪水点据作为调试的依据。此值一般在0.2～0.5之间取值。注意B值与

W值之间有

m

相依性，二者对径流模拟的结果相互有影响。

稳定下渗率FC的调试，FC对水源划分起决定性作用，但至今对FC的研究还不充分，目前可对FC作一些简单处理，如果FC取为一个常数使模拟效果不好时，可考虑FC作为变动参数。考察FC是否合适主要观察地下径流的模拟精度。

地下水消退系数KKG的调试，主要观察洪水退水段的拟合精度。

不透水面积参数IMP对模型拟合精度影响很小，一般取0～0.05。

2.2新安江三水源模型

一、新安江三水源模型结构

新安江三水源模型包括4个计算环节：流域产流计算；径流的划分；蒸散发计算；汇流计算。流域的产流计算和蒸散发计算与新安江二水源模型相同，水源划分则完全不同，因此汇流计算也不相同。新安江三水源模型结构示意图如下

二、新安江三水源模型流域产流计算

新安江模型产流部分的计算是蓄满产流模式， 蓄满产流指在流域包气带土湿满足田间持水量以前不产流，所有的降雨都被土壤吸收；而在土湿达到田间持水量之后，所有的降雨（除去同期的蒸散发）都产流。在产流后，流域包气带土壤的下渗能力为稳定下渗率，下渗的水分成为地下径流和壤中流，超蓄的部分成为地面径流。考虑到流域内各点的蓄水容量并不相同，实际产流时常常是在部分面积上产流，新安江模型引入流域蓄水容量曲线来刻划流域内各点蓄水容量的不均匀性，把流域内各点的蓄水容量概化成如图2.3所示的一条抛物线（也可概化成其它函数形式），其方程为：

b

mm m W W F f ⎪⎪⎭⎫ ⎝⎛--='11 （2-1） 式中： FR F

f =为产流面积； 'm W ——是点蓄水容量,mm ；

f ——是蓄水容量小于等于蓄水量'm

W 的全部点的面积,2Km ； F ——是单元流域面积, 2Km ；

b ——是蓄水容量曲线指数。

mm W ——是流域内点最大蓄水容量,mm ；其与流域平均蓄水容量M W 的关系为: ()M mm W b W ⋅+=1 （2-2） 引入蓄水容量曲线后，在降雨过程中，只有在满足了包气带蓄水容量的面积上才可能产生径流，其余面积上不产生径流。扣除蒸发后的有效降雨量e P 中，未转化为净雨产流量R 的水分通过下渗进入流域包气带土层，补充土层蓄水量W 。引入流域蓄水容量曲线后的产流计算式为：

当mm e W a P <+时为部分面积产流

b mm e M M e W a P W W W P R +⎪⎪⎭⎫ ⎝

⎛+-⋅++-=1101 （2-3）

当mm e W a P ≥+时为全流域面积产流 0W W P R M e +-= （2-4）