李 匡,朱成涛,胡宇丰,徐海卿
(1.中国水利水电科学研究院 北京中水科水电科技开发有限公司,北京 100038;2.二滩水电科技开发有限公司,四川 成都 610021)
水电站施工期水位预报主要为施工断面导进口、导出口的水位预报,对施工安全有着重要的作用。对于受雍水影响较小的施工断面,可以利用水位流量关系曲线进行水位预报(根据预报流量查水位流量关系曲线得到预报水位)。施工期利用水位流量关系曲线进行水位预报时存在两个难点:(1)施工断面的水位流量关系曲线不容易获得。一般来说,施工断面的水位资料可以通过水位计测量得到,而施工断面流量不易取得,导致在率定水位流量关系曲线时存在困难。有学者提出可通过水面比降法或水力学方法计算得到施工断面的水位流量关系曲线,但是受各种因素影响,往往存在较大的误差[1-3];(2)由于电站施工,经常会对施工断面河道造成一定的影响,导致天然河道发生改变,例如:块石滚落河道导致河道底部升高等。此时,如果直接查水位流量关系曲线,得到的水位预报值误差较大。
针对以上两个问题,提出了解决方法:(1)水电站的设计代表站一般距离施工断面较近,因此可以将设计代表站的流量值作为施工断面的流量。而设计代表站一般为条件较好的水文站,水文站流量容易获得;(2)由于预报发布时间的水位和流量是实测得到的,预报时间的流量可以通过流量预报得到,可以利用已知条件,利用水位系数法计算得到预报水位。水位系数法不直接利用水位流量关系曲线进行预报,克服了河道变化对水位预报造成的影响。
在水位流量关系上任取一点q,我们称q点的切线斜率为水位系数。
如图1所示,水位流量关系曲线上任意一点q处切线的斜率为
如果知道水位流量关系曲线的函数表达式,则可以求出任意一点切线的斜率。在实际应用中,如果水位流量关系曲线函数表达式难以获得,也可将水位流量关系曲线横坐标Q按照相同的Δq进行离散,近似计算出各点的水位系数。
图1 水位流量关系曲线
图2 流量水位系数曲线
将流量和水位系数组合,则可得到流量水位系数曲线,如图2所示。可以通过曲线拟合得到其函数表达式k=f(Q),也可以将其制作成离散的水位系数表。
利用水位系数法进行预报时,根据预报发布时间的实测流量qt查出水位系数kt,预报时间的预报流量qc查出水位系数kc,利用下式即可计算出预报时间的预报水位。
式中:zc为预报时刻的预报水位;zt为预报发布时刻的实测水位;k为水位系数,取kt和kc的平均值。
如图3所示:ΔZt为采用kt进行计算得到的水位增幅;ΔZc为采用kc进行计算得到的水位增幅;ΔZ为采用kt和kc的平均值进行计算得到的水位增幅。在k取kt、kc、kt和 kc的平均值时 3 种情况下,当 k取kt和kc的平均值时,预报值最接近实际值,产生的误差最小,因此实际计算中,k取kt和kc的平均值。
图3 水位系数法计算示意
雅砻江是金沙江左岸的最大支流,干流上目前有两河口、锦屏一级、锦屏二级、官地等电站正在施工,需要对各电站各施工断面的水位进行预报。两河口电站的设计代表站为雅江水文站,锦屏一级和锦屏二级电站的设计代表站为锦屏水文站,官地电站的设计代表站为打罗水文站。通过设计代表站测得流量,结合各施工断面测得的水位,利用水位系数法进行水位预报,其步骤如下。
(1)利用已有的施工断面水位、设计代表站流量资料建立施工断面水位流量关系曲线。
(2)计算施工断面水位系数,得到水位系数计算公式或水位系数表。
(3)对施工断面进行流量预报。
(4)利用水位系数法公式进行水位预报。
利用2009年各水位站水位资料及其对应的设计代表站的流量资料建立的水位系数计算公式为:
两河口导进口:
两河口导出口:
锦屏一级导进口:
锦屏一级导出口:
锦屏一级景峰桥:
锦屏二级猫猫滩:
锦屏二级大水沟:
官地导进口:
官地导出口:
式中:k为水位系数;Q为流量。
利用以上步骤分别对各电站各施工断面进行水位预报,统计2010,2011两年汛期内的预报精度。表1为施工断面的流量预报精度表,表2为施工断面水位预报成果表。
分析表1、表2可以得出以下结论:
表1 断面流量预报精度
(1)水位系数法预报精度较高,锦屏一级导进口、锦屏二级猫猫滩、锦屏二级大水沟、官地导进口12h预见期平均绝对误差在0.15m左右;锦屏一级导出口、锦屏一级景峰桥、官地导出口的12h预见期平均绝对误差在0.1m左右。两河口断面水位预报精度较低,12h预见期平均绝对误差在0.25m左右,最大绝对误差为0.3m。各断面12h预见期平均相对误差均在0.01%左右,可以忽略不计。
(2)在利用水位系数法进行预报时,决定水位预报精度的因素有3个:流量预报精度、水位系数准确度和实测水位准确度。其中水位系数可以通过多年的资料进行率定可以达到较高的准确度;实测水位也可以通过多次测量或者平滑处理达到较高的准确度。而流量预报的精度对水位预报影响最大,流量预报精度高的断面水位预报精度也较高,流量预报精度低的断面水位预报精度也较低。
(3)有的断面最大绝对正误差和最大绝对负误差较大,例如两河口导进口、锦屏一级导进口、官地导进口,其最大绝对正误差和最大绝对负误差均比相对应导出口的误差大。这主要是由于水位测量误差造成的,导进口进水口由于壅水影响,受回水波动影响较大,造成测量出现较大的误差。
(4)随着预见期的增长,水位预报的精度也随之降低,这是由于流量预报精度的降低造成的。
利用施工断面的实测水位,设计代表站的实测流量建立施工区断面水位流量关系,再采用水位系数法进行施工区断面水位预报,克服了施工对施工区断面造成的影响,方法简单易操作,取得了较好的预报精度,值得推广使用。在使用过程中,应该注意以下几点:(1)每年率定一次水位系数,减小由于水位系数造成的误差;(2)提高流量预报的精度,因为流量预报精度是水位预报精度的前提条件;(3)减小实测资料的测量误差,提高预报精度;(4)水位系数法适用于无围堰水库或者围堰水库无调节性能的施工断面,对于有围堰水库调节作用的施工断面,可以采用调洪演算法进行水位预报[4-7]。
表2 施工断面水位预报精度
[1]林三益.水文预报(第二版)[M].北京:中国水利水电出版社,2001.
[2]詹道江,叶守泽.工程水文学[M].北京:中国水利水电出版社,2000.
[3]范可旭,伍宁.拟定坝址水位流量关系的水力学方法[J].人民长江,1999,30(3):30-31.
[4]熊金和,雷辉光.构皮滩水电站施工期洪水预报综述[J].贵州水力发电,2010,24(2):4-7.
[5]邓柯,仇海兵,李艳.皂市水库施工洪水预报模型[J].人民长江,2005,36(7):43-44.
[6]刘文俊.街面水电站库区施工期洪水预报研究[J].中国农村水利水电,2005(9):38-41.
[7]邹少龙.棉花滩水电站施工期洪水预报[J].水力发电,2001(7):17-18.