混合产流模型和新安江模型在半干旱地区的应用比较.docVIP

混合产流模型和新安江模型在半干旱地区的应用比较 路芳，任立良，袁飞 河海大学水文水资源学院，南京(210098) E-mail： 摘 要：为了比较混合产流模型和新安江模型在半干旱地区的应用效果，选定辽宁老哈河支流坤头河流域的小城子水文站的集水流域作为研究区域，分别采用这两个模型进行模拟计算，并对模拟的结果进行分析。 关键词：混合产流模型，新安江模型，半干旱地区 1. 引言 水文模型的研究起始于20世纪60年代，经历了从黑箱模型到集中式模型再到分布式模型的发展过程。20世纪80年代中期以来，随着计算机技术、地理信息系统和遥感技术的出现，分布式水文模型产生了巨大的发展，例如TOPMODEL模型，SHE模型，SWMM模型等[1]。 如着眼于影响次降雨—径流关系的因素，则自然界的产流只有两种基本的模式：蓄满产流模式（次降雨与总径流不受雨强的影响）和超渗产流模式（次降雨与总径流的关系受雨强的影响）[2]。所谓的蓄满产流是指在土湿满足田间持水量以前不产流，所有的降雨都被土壤吸收；而在土湿满足田间持水量以后，所有的降雨（减去同期的蒸散发）都产流。而超渗产流是指在土湿未达到田间持水量以前因雨强大于渗强而产流[3]。湿润地区的产流方式是蓄满产流，所以蓄满产流模型适用于湿润地区；干旱地区的产流方式是超渗产流，所以超渗产流模型适用于干旱地区。但实际上，在许多地区，既非干旱也不湿润，常称为半干旱半湿润地区。在这类地区的流域，有些洪水是蓄满产流，有些是超渗产流。同一场洪水，前期可能是超渗产流，到后期又变为蓄满产流；在较干旱地区，一个以超渗产流为主的流域，遇上长期连绵的低强度降雨的洪水，也可是蓄满产流；类似地，在较湿润地区，以蓄满产流为主的流域，久旱后遇雨强特大的暴雨，也会有超渗产流发生等等。这就表示混合产流模型研究的必要性。 下面就把混合产流模型和新安江模型分别应用到半干旱地区，分析比较结果。 2. 模型简介 胡春歧研究的混合产流模型的结构如图2-1所示。该模型通过土壤蓄水容量曲线与下 式的有机结合。分为α与1-α两部分，在α部产流，面积上，量后才产流。强很小，- 1 - 满而产流[4]。 新安江模型是一个完整的降雨径流流域模型，用于湿润半湿润地区。模型的产流部分就是蓄满产流模图 型，蒸散发计算采用三层模型，流域汇流中地面汇流采用经验单位线，地下径流汇流采用线性水库，河道汇流采用马斯京法[3]。其结构见图2-2。 图2-2 新安江模型结构示意图 3. 模型应用 为了检验混合产流模型的适用性，并与新安江模型比较，采用辽宁老哈河支流坤头河流域1967年到1976年十年的实测资料计算分析。该领域的出口处是小城子水文站，集水面积有803平方公里，位置是北纬41°30′－41°54′，东经118°24′－119°12′，含有2-7个雨量站，1个蒸发站，1个水文站。 本研究应用数字高程流域水系模型（Digital Elevation Drainage Network Model，缩写为DEDNM），利用栅格型高程数据自动生成流域水系、河网及其拓扑关系。栅格DEM数据来自美国国家地球物理数据中心提供的全球陆地一公里基础高程GLOBE（Global Land One-kilometer Base Elevation）数据，其空间分辨率为30秒。之后，利用DEM结果把流域资料分配到空间上，以此作为模型的输入。 小城子流域混合产流模型和新安江模型模拟的参数数值及结果见下表。 表3-1 小城子流域混合产流模型模拟参数数值表 参数 数值 表中：kc—蒸散发系数，Qrgo—初始地下径流流量，b—张力水蓄水容量曲线指数，bx—下渗能力分布曲线指数，fo—流域平均最大下渗能力，fc—流域平均稳定下渗能力，k—下 - 2 - 渗能力衰减系数，kg—地下水日消退系数，imp—不透水面积占全流域面积百分数，c—深层蒸散发系数，x—子河段的马斯京根演算系数。 表3-2 小城子流域新安江模型模拟参数数值表 参数 数值 表中：kc—蒸散发系数， Fr—初始产流面积，s—初始流域自由水水深，Qrsso—初始壤中流流量，Qrgo—初始地下径流流量，sm—自由水蓄水库容量，b—张力水蓄水容量曲线指数，ex—自由水蓄水容量曲线指数，imp—不透水面积占全流域面积百分数，kss—自由水蓄水库壤中流日出流系数，kg—自由水蓄水库地下水日出流系数，kkss—壤中流日消退系数，kkg—地下径流日消退系数，c—深层蒸散发系数，x—子河段的马斯京根演算系数。 表3-3 小城子流域混合产流模型模拟结果 年份 确定性系数（%） 年径流误差（%） 表3-4 小城子流域新安江模型模拟结果 年份 确定性系数（%） 年径流误差（%） 总的 总的 140120流量（m3/s）1008060402001