1、 1 无无实测实测径流资料地区径流资料地区编制编制施工期施工期水文水文预报方案预报方案探讨探讨 (孙恒兵陈端丹 江力 计强) (十堰市水文水资源勘测局 44200) 摘要摘要:在无实测径流地区编制施工期水文预报方案,面临预报参数的选用与分析,笔者通过对白沙河水电站施工期预报方案的编制, 在作业实践中进行检验, 满足了白沙河水电站等大(2)型工程建设需求。 关键词关键词: 水文预报参数借用 施工期水文预报编制 1 引言引言 白沙河水电站属等大(2)型水利水电工程,位于湖北省竹溪县兵营乡,水库坝址处于堵河西支泗河一级支流泉河上,泉河全长 84.6km,流域面积 901.0km2,坝址控制集水面积
2、812.61km2,占全流域面积的 90.2%。 泉河流域洪水由暴雨形成,洪水的时空分布特性与暴雨相应,径流有明显的季节性。汛期为 510 月,分为夏汛和秋汛,夏汛洪水多发生在 6 月中下旬至 7 月中旬,秋汛洪水则多发生在 9 月中旬至 10 月上旬。夏汛洪水以单峰居多,一次过程约为 13 天,具有陡涨陡落的特点;秋汛洪水多为复式峰,一次洪水过程约为 35 天。 2 水文资料水文资料情况情况 流域内无水文站,缺乏水位、流量、泥沙等水文资料,也缺乏历史洪水调查资料。邻近流域有汇湾河汇湾水文站、竹溪河延坝水文站、泗河新洲水文站以及官渡河老码头水文站。 汇湾河汇湾水文站(现迁建为鄂坪水文站)流域面
3、积 1650 km2,资料年限为 19862006年,2001 年后受白果坪水电站蓄水影响,流量失去代表性。 竹溪河延坝水文站流域面积 598 km2,与泉河流域面积最为接近,但流域内坡降平缓,受上游水利工程以及人类活动影响较大, 人为因素对流域的产汇流条件影响较大, 水文系列缺乏一致性。 泗河新洲水文站流域面积 4660 km2, 比泉河流域的流域面积大 5 倍多, 属于大河控制站,泉河流域属小河站,两者的产、汇流规律不同,其流域特性与泉河流域存在较大差异,水文资料不宜采用。 官渡河老码头水文站流域面积 1854 km2,资料年限为 19862010 年,共计 25 年水文资料。 老码头水文
4、站代表的官渡河流域其自然地理和水文气象条件与泉河流域较相似, 两个流 2 域被同一分水线隔开,流域内植被良好,森林覆盖率达 70%以上;干流两岸支流密布,河道弯曲,为典型的河曲地貌;流域径流由降水产生,其洪水与暴雨在时空分布上基本相近,洪水过程以单峰为主,具有陡涨陡落的特点;且水文站以上的流域内均无大、中型蓄水工程,人为因素对流域的产汇流条件影响较小, 水文系列具有较好的一致性。 故采用老码头水文站为参证站来分析泉河流域的产汇流预报方案。 流域地理特征对照表见表 1 流域多年平均降水量对照表见表 2。 表 1 流域地理特征对照表 水系 河名 河长(km) 流域面积 (km2) 水系形态 河道坡
5、度() 流域平均高程(海拔米) 河流弯曲系数 河网密度(km/ km2) 堵河 泉河 84.6 901 树枝状 7.0 1094 1.5 0.3 堵河 官渡河 127 2885 树枝状 4.6 1349 1.7 0.3 表 2 流域多年平均降水量对照表 1 2 3 4 5 6 7 8 9 10 11 12 年总量 泉河 161 27.6 54.7 87.4 134.6 158.4 184.5 184.2 116.6 95.0 41.8 19.9 1120.8 官渡河 15.7 28.5 51.8 89.8 138.3 161.0 165.6 153.0 113.2 91.4 43.3 18.5
6、 1070.0 3 方案编制方案编制 白沙河施工期水文预报方案产流采用降雨径流相关图法, 汇流采用综合单位线法及谢尔曼单位线法进行编制。 由于泉河缺乏实测的径流资料, 编制方案时根据蓄满产流法的原理, 借用相邻官渡河流域产流模型参数,汇流采用泉河地理参数分析和相邻官渡河流域官渡河单位线经过 S-过程线转换。 3.1 产流产流参数选用与分析参数选用与分析 流域平均降水量的计算流域平均降水量的计算:根据坝址以上流域平均降水量计算站点,采用面积加权法,求得整个流域平均降水量。 月 份 河 名 niiinnFfPFfPfPfPP12211. 3 次洪径流量 R 的计算:采用流域退水曲线分析,选择孤立单
7、式洪水过程,点绘 QT 过程线图,做退水曲线簇,取外包线,作为综合退水曲线,分割地下径流,计算实测次洪径流量 R。 前期影响降水量前期影响降水量 PaPa 的计算的计算:前期影响雨量 Pa,计算出第一场洪水之前的 Pa,以后在整个汛期内采用连算,公式 tt , aa1t , aPPKP,但必须控制 mtaIP1, 。用上式计算也可取连续大暴雨之后的aP等于 Im,由此向后逐日推算。 按经验公式计算:t 日有雨 Pa,t+1=k(Pa,t+Pt) t 日无雨 Pa,t+1=kPa,t Pa,t某日 t 的前期影响雨量(mm) ; Pa,t+1次日(t+1)的前期影响雨量(mm) ; K与土壤蒸发
8、能力有关的日折减系数。 根据已确定的 K 值,用日雨量资料逐日计算 Pa值。 流域最大损失量流域最大损失量 ImIm 的确定的确定:流域的最大损失量由实测雨洪资料分析而得。方法是在历年资料中选择久旱不雨后, 一次降雨量较大能满足流域最大损失值且全流域产流的资料, 计算该次降雨的流域平均雨量 P 及所产生的径流量 R,因为久旱,Pa0,则: Im=P-R-E雨 E雨雨间蒸发量。 本流域最大损失量 Im 值,经多次计算分析取用 80mm。 消退系数消退系数 k k 值的确定值的确定:流域消退系数按下式计算,E 为流域蒸发能力,可用 E601 观测器观测的水面蒸发值作为近似值。 降降 雨 径 流 相
9、 关 图 的 建 立雨 径 流 相 关 图 的 建 立 : 采用老码头站 19862010 年 76 次洪水资料,经过降雨径流关系分析,确定每一场洪水的 P、Pa、R,建立老码头水文站以上流域降雨径流三变数关系图。 4 产流汇流方案率定与检验产流汇流方案率定与检验 汇流参数选用:汇流参数选用: 在编制汇流参数时采用老码头水文站代表的官渡河流域, 其单位线的分析用综合单位线法与经验单位线法,即:一是根据泉河流域地理参数推求综合单位线;二是借用老码头水文站经验单位线, 由 S-过程线转换为泉河流域经验单位线, 单位净雨为 10mm,单位时段t=1h,流域小单位时间取 1 小时比较符合实际。 产流方
10、案预报作业产流方案预报作业成果检验成果检验:用 1986 年2010 年 25 年实测降水量成果,选取场次降Im1EK 4 水 76 次进行预报检验, 预报项目的精度评定用合格率表示, 降雨径流关系合格率为 84.2。 汇流方案检验与评价汇流方案检验与评价:选取老码头水文站 19862010 年实测大、中、小洪水进行洪水预报检验,预报项目的精度评定用合格率表示,洪峰流量合格率 QR=94.4%,确定性系数DC=0.97,按表 5-1 标准,评价为甲等(QR85.0%,DC0.90) 。 作业预报评价:作业预报评价: 5 水文水文预报方案预报方案及需要注意及需要注意的问题及应对方法的问题及应对方
11、法 白沙河施工期水文预报方案编制在作业预报实践中, 我们发现在用 中国洪水预报系统预报时,它最适宜的状况是流域平均降水,均达到蓄满产流的理想情况下,有比较满意的预报结果,当出现下列情况时: 1、流域暴雨中心在上游或下游,预报的的洪峰与实测洪峰相比,可能出现推后或提前的情况,有时误差比较大。因此在预报实践中要考虑分析暴雨中心点位置,合理实施人工实时校正,一般通过多次预报实践,在预报中采用逐步、连续预报作业手段,一般取得满意的效果。 2、受上游水库调蓄产生的洪水过程与自然状态下的洪水过程相比,前者扭曲变形,洪峰提前或推后,洪水过程延长,对预报结果影响较大。采用对策是及时准确地收集上游水库发电、排沙
12、、溢洪等信息,对预报过程进行人工干预,拟合最佳预报过程,从而提高预报时效性和准确性。 3、由于施工进度不停在变化,围堰挡水转为坝体挡水,随着坝体不断浇筑,坝高增高,对来水产生的调蓄作用越来越明显,从而影响预报精度。采用对策是预报时根据施工进度,对预报成果实时修正。 4、由于泉河流域无水文站,缺乏历史实测水文资料,汇流部分没有实测的历史洪水过程来检验汇流成果,因此在预报实践中,要经过实测洪水资料实时检验汇流参数,使预报成果更贴近实际。 6 结论 通过对白沙河水电站施工期水文预报的编制,对今后在十堰地区无径流流域编制水文预报积累了经验,对于在作业预报中出现问题提供参考依据。 6、参考文献 1、 白沙河施工期水文预报方案编制 2、 水文情报预报规范 (GB/T 22482-2008) 5 3 相关水文预报书籍