黑龙江高寒山区冰坝最高水位预报模型研究.docxVIP

黑龙江高寒山区冰坝最高水位预报模型研究 1 冰凌洪水的形成 黑龙江省上游的冰坝通常发生在河流的交汇处，属于寒性和山区河流。这条河很窄，深度和宽度不同。它很大，有很大的落差，气候寒冷，地面冰雪，土壤硬化。关闭的时间为170天，有些土壤为300天。大多数流域内都存在永久性冷冻土，这对冷冻土具有明显的水文效应和良好的产水条件。每当冬春气候异常, 河槽蓄水 (冰) 量偏大, 积雪偏厚, 春季气温急剧回升并伴有降雨年份, 都会在上游形成强大冰凌洪水, 在河道弯曲、束狭等输冰不畅江段或冰盖前缘形成冰凌卡塞和冰坝。黑龙江高寒山区性河流冰坝, 是在倒开江的条件下, 由融雪、融冰和降雨径流形成的冰凌洪水。冰坝形成过程是强大的水流动力, 促成冰盖破裂聚集、堆积而形成的叠加组合过程。冰坝程度决定于河槽前期蓄水 (冰) 量、冬季积雪和融冻期降水量以及气温回升过程决定的倒开江程度和冰盖强度。 冰坝凌汛多为单一洪水过程, 洪水总量可表示为: W总=W冰+W水+W雪+W雨(1) 式中:W总为冰坝凌汛洪水总水量;W冰为河槽结冰水量;W水为河槽蓄水量;W雪为融雪水量;W雨为融冻时降雨径流量。 黑龙江上游属山区河流, 落差大, 河槽蓄水量少, 可忽略不计, 则凌汛水量和洪峰高度决定于融雪、融冰和降雨径流量。洪峰附加壅水高度与倒开江程度和冰盖强度有关。当冰凌洪水的峰量关系单一时, 冰坝最高水位可近似表示为: 式中:Hm为冰坝最高水位, 以变幅计, m;H0为融雪、降雨、融冰径流混合过程的正常水位, m;R雪、R雨、R冰分别为融雪、降雨和融冰水量, m3;△H为冰凌叠加和冰盖强度影响下的附加壅高水位, m;h为临界冰厚, m;准为冰盖抗弯强度系数, kg/cm2。 式 (2) 基本反映了冰坝凌汛形成的水量平衡和热量平衡关系, △H还间接反映了倒开江条件下, 冰凌洪水的叠加组合形式和程度。 2 计算冰坝洪水的计算 2.1 融冻期的蒸发能力 融雪和降雨径流, 多为开江期的一次混合径流过程, 在黑龙江上游, 为冰坝凌汛主要水量来源。分析寒区冻土水文模型, 水量平衡方程式采用下式: 式中:R为融冻期降雨径流量, mm;P为冬季降水量, mm;E为蒸散发 量 , mm;Wm为流域最 大蓄水容 量 , mm;W0为起始蓄水量, mm。 在融冻时, 融雪和降雨径流的产流方式和径流量受到冻土融冻深度及土壤含水量影响, 容易达到田间持水量, 产流方式符合蓄满产流理论。但是受冻土特性和土壤水分动态规律的影响, 产流方式和参数的定量随冻土融冻深度有所变化。在融冻期, 产流层土壤可分为融冻 (包气带) 层和冻结 (包冰带) 层。融雪和降雨径流和参数的定量与无冻期雨洪的计算存在区别。融冻期雨雪径流量的计算可表示为: 式中:PE为有效降水 (雪) 量;W0上、W0下分别为上下层计算的起始蓄水量。将流域最大蓄水容量Wm分为上下两层, 从冻土融冻日起, 随气温增加和冻土融冻深度变化, 上层从0~Wm变化, 下层则从Wm~0变化 , 分别计算W0上和W0下。令上层按蒸发能力蒸发, 下层按冻土消融速度变化, 则上下层起始蓄水量为: 式中:W0, t和W0, t+1为时段始末蓄水量, mm;KD t, a为下层蓄水量消退系数和指数, 经分析冻土深度变化近似于直线, a可取0.8。 运用以上公式计算, 需要考虑冻土融冻期上下层最大蓄水容量的变化过程, 蒸发能力随冻土过程和下垫面情况不同而变化。对常用的观测水面蒸发仪器进行器皿修正和高程修正外, 还要作冻土蒸发修正。冻土之上有积雪存在时, 蒸发发生在雪面, 而冻土温度在冰点以下, 蒸发受到抑制, 冻土蒸发修正系数C从0变化到1。 根据上述方法, 分析计算了大小兴安岭地区的科后、山口、北安、清溪, 西沟等水文站的融雪和降雨径流关系, 收到了良好的效果。可将综合的R=f (P+W0) 关系用于计算黑龙江上游冰坝凌汛水量。 2.2 河网蓄冰量估算方法 黑龙江上游无水利工程调节, 冬季河流最小流量不足20m3/s, 河槽蓄水量主要以蓄冰为主。最大蓄冰量, 用断面法估算可表示为: 式中:hm为最大冰厚, m;B为冰面宽度, m;L为河道长度, m。分河段计算, 并考虑起始冰面不同宽度, 式 (6) 可表示为: 式中:B为开始封冻时上下断面平均宽度 , m;BC为最大冰厚时, 上下断面冰面平均宽度, m; hm为上下断面平均最大冰厚m;△L为两断面的间距, m;n为计算河段数。 河网蓄冰量的另一种估算方法是采用无冻期断面退水曲线与冻结期退水曲线水量之差来估算。 以上估算的河槽蓄冰量与冰坝水位高度具有直接关系, 但最大冰厚多出现在3月下旬前后, 离4月下旬开江日期冰厚还存在升温减薄过程。如考虑临界冰厚和冰盖强度, 还需进行冰盖厚度和强度的计算。 2.3 冰流速及冰厚减薄