温州珊溪水库施工期洪水预报模型及运行(1)程

摘要：论述了应用姜湾径径流模型编制大型水利建设工程——珊溪水利施工期洪水预报模型的结构原理、模型参数、计算分析和模型的运行，包括建模以来作业预报的过程、精度和应用效果。

关键词：洪水预报模型 运行 珊溪水库 浙江

珊溪水库是目前浙江省在建的最大水库工程，建成后总库容18.24亿m³，正常蓄水位142m（黄海基面），相应库容12.24亿m³，防洪库容2.12亿m³，系多年调节水库，装机容量200mw，年发电量3.55亿kw·h。

水库坝址位于飞云江上游文成县珊溪镇附近（图1），坝址以上河长92.4km，控制流域面积1529k m²，占飞云江流域面积的47%。

坝址以上流域地形以高山为主，河道蜿蜒曲折，水流湍急。流域内植被良好。

坝址以上地处亚热带季风区，气候温和，雨量充沛，属暖温带多雨气候。坝址以上流域多年平均降量量1876.9mm，年降水量在1280~2458mm之间。降水量年内分配不均匀，主要集中在5~7月的梅雨期和8~9月的台风期，台风暴雨是造成本流域洪水灾害的主要天气因素。

2  珊溪水库施工期洪水预报模型

珊溪水库施工期洪水预报模型采用浙江省水文勘测局研制的姜湾径流模型，该模型的径流系四水源组成（地表径流；壤中回归水径流；壤中径流；基岩风化层径流），而各种水源认为是产生在介质特性发生急剧变化的交界面上，其模型结构和流程如图2所示。

图2   姜湾径流模型流程图

r－产流量（mm）；α－产流面积系数，即α=r/（p-e）；f₁－淋溶层稳定下渗率（mm/△t），一般取2.5~10.0mm/h；f₃－沉积层稳定下渗率（mm/△t），一般取1.0~1.5mm/h；r_s－地表径流；r_g－地下径流，包括r_i和r_ig(mm)；x₁－淋深层垂向调蓄系数；x₂－沉积层垂向调蓄系数；r₁－壤中总径流；r_ig－基岩风化层径流；r_gg－母质层总入流；r_il－壤中回归径流；r_i2－壤中径流；k_i2－淋溶层容积(mm)；k_s、k_il、k₁₂－坡面淋溶层侧向汇流系数；kk_g－基岩风化径流消退系数；q_s、q_il、q_i2、q_g、q分别为四种水流在出口断面的流量和总出流量。

2.1  预报流域的划分

根据飞云流域的特性和施工期的特殊性，把珊溪坝址以上流域划分为二个单元，即百丈口水文站以上流域，其控制面积为866 km²，流域内有三插溪、司前、仙居雨量站和百丈口水文站，其中三插溪水库在1998年7月建成，其上267.5 km²的来水已受控制调节；百丈口至珊溪区间流域面积为663 k m²，内设有外墙、黄坦、西坑雨量站和珊溪水位站。

2.2  珊溪水库预报模型参数的率定

2.2.1  珊溪坝址以上流域面平均雨量计算

报汛站网的布设原则是充分利用原有的报汛站点，以最经济的站点数达到能够掌握所需水情的变化，满足水库水文预报的精度和将来运行调度的需要为目标。经过对流域内原有的雨量站优化计算分析，最后拟定分布比较均匀的8个报汛站，其中二处报汛站为新增。按单元用算术平均法计算面雨量，单元流域的根据站（表1）。

表1        单元流域面雨量计算根据站表

2.2.2  产流量计算模型参数的率定

采用日模型进行率定，应用浙江省水文勘测局的“多层筛选网格计算法”软件，对年径流、次径流进行同步模拟计算——根据地区、流域的特性、设置各参数的合理取值范围，然后进行第一层面的网格计算，淘汰不合理的参数组；第二层面则对次径流以《水文情报预报规范》的规定为标准进行网格计算，从中选出年径流深误差在许可误差范围内，次径流深合格率最高的参数组群。再从中选出次径流深合格率最高，同时年径流深误差最小的参数作为最佳参数组。

2.2.3  汇流计算模型参数的率定

首先根据洪水退水过程线自下部至上部的拟合程度判断其参数的合理性。由于退水末段陡与缓主要决定于x、f₃和kk_g这三个参数，而kk_g可用退水过程线直接计算，如果各次洪水退水末段拟合好了，就可认为x和f₃求得的数值也是合理的。其次是向上综合壤中流参数k_f2和k_i2，可从两方面判断壤中流参数的合理性：一方面判断流量过程线腰段的拟合程度；另一方面判断壤中流的峰现时间（壤中流的峰现时间应在地表径流峰值出现之后，在基岩风化层径流峰值之前）。最后，对于大暴雨造成的特大洪水，关键是拟合其上部的峰断和峰现时间，即壤中回归水径流和地表径流的汇流参数f1、k_s和k_il的率定。各径流拟合、叠加过程线参见图3；珊溪水库施工期洪水预报模型参数见表2。

表2      珊溪水库施工期洪水预报模型参数表

中国大学排名

共2页: 1