1. 概述
汉江属于长江最长支流,发源于秦岭南麓,干流流经陕西、湖北2省,于武汉市汇入长江。汉江干流通常分为三段:丹江口以上为上游,丹江口至碾盘山为中游,碾盘山以下为下游。丹江口水库位于湖北省丹江口市,丹江口水库至王甫洲河段河道现状设计防洪标准为1964年同大洪水,堤防加固设计水位按照1964年实测洪水位确定,丹江口水库建成至今河道已发生了明显的冲刷变化,在现状河道下发生1964年同大洪水河道水位较1964年变化程度有待进一步的研究。目前,国内外研究学者针对汉江中下游河道冲淤等方面进行了深入研究,并取得了一定的成果。林云发 [1] 对汉江中游近期冲刷状况研究表明丹江口水库建成后河道呈现自上游至下游冲刷变化,丹江口–襄阳河段基本达到冲淤平衡;郭海晋 [2] 等对汉江雅口枢纽坝下河道冲淤和水面线变化研究表明雅口枢纽坝建成后,坝下河道表现为近坝段微冲,其余河段淤积。其研究结果表明,汉江中下游河段水利工程建设均会对坝下河道产生冲刷,但随着时间变化会逐渐达到冲淤平衡状态。分析河段位于丹江口–襄阳河段,经过五十多年的河道冲刷变化,河道从冲刷状态逐渐达到冲淤平衡状态。
针对河道的冲淤变化对河道水位影响方面研究,国内外目前仅进行了定性分析,深入研究的相对较少。国内余明辉 [3] 等对长江中下游洪水位河床冲淤关系研究表明,长江中下游河道淤积和河道形态变化是造成洪水位抬高和高洪水位持续时间加长的主要原因;郭海晋 [2] 等对汉江雅口枢纽坝下河道冲淤对水面线变化研究表明雅口枢纽坝下河道水位随着河道的冲刷下切而降低,随着河床淤积而抬升;宋涛 [4] 论文表明随着河道整治开挖能够明显降低河道水位。对丹江口–襄阳河段河道冲刷后水位变化目前尚无明确的研究,本文结合已有研究方法,对丹江口水库至王甫洲河段河道冲刷条件下中高水洪水位变化进行分析,明确在河道冲刷情况下中高水的水位变化规律,为沿线的岸线保护、水利工程设计以及河道管理提供参考。
1.1. 河段介绍
丹江口水库枢纽位于湖北省丹江口市,是亚洲第一大人工淡水湖,也是国家南水北调中线工程水源地。丹江口水库一期工程始建于1958年,1973年底建成,坝顶高程162 m,正常蓄水位157 m,相应水库面积745 km2,库容174.5亿m3。二期工程自2005年9月开始进行,2011年7月大坝坝顶从162 m提高至176.6 m,正常蓄水位170.0 m,相应库容增加到290.5亿m3。是一座以防洪、供水为主,结合发电、航运等综合利用的大型水利工程。丹江口水库建成后,下游河道来沙量仅占建库前的0.45% [5]。
王甫洲水利枢纽位于丹江口水库枢纽下游30 km,属于湖北省老河口市,王甫洲水利枢纽1995年开始新建,于1999年并网发电,坝址设计洪水(P = 2%) 18,070 m3/s,校核洪水(P = 0.667%) 22,000 m3/s;丹江口大坝加高后,王甫洲水利枢纽提高至100年一遇设计洪水,设计洪水位88.11 m,正常蓄水位86.23 m,调节库容0.28亿m3。王甫洲水利枢纽来水在11,100 m3/s以上时,全部闸门开启并按泄水闸泄洪能力敞泄。
本次以丹江口水库坝址下游王家营水位站为起始断面,王甫洲水利枢纽库区老河口水位站为终点断面进行分析,河段全长26.65 km,河段内现在无其他大的水电工程,本次根据实测1967~2019年实测断面进行分析,河段断面分布见图1。
分析河段从1967年至今已建有3座桥梁,在DM1~DM6段、DM15~DM18段建有河道整治工程,岸边修建有堤防,其对河道水位变化影响有限。其他河段受人类活动影响很小,河道变化主要影响因素为河道冲刷。
Figure 1. Network of hydrographic stations and analysis of section distribution
图1. 河段水文站网及分析断面分布图
1.2. 河段冲淤变化
自丹江口水库建库以来河段基本呈现冲刷状态,合计冲刷量达到7850万m3 [6],从1967年~1978年冲刷量2761万m3、1978年~1987年冲刷量1401万m3、1987年~2005年冲刷量869万m3、2005年~2012年冲刷量2776万m3、2012年~2019年冲刷量2804万m3,河段冲刷进入冲淤平衡阶段,此段主河槽粗化已基本完成,卵石大量出露,2005年~2019年冲刷量加大主要为人为采挖影响 [1]。河道纵向深泓线整体呈现冲刷状态,冲刷最大深度达到7 m,河段冲刷总体呈逐年减弱的趋势。本次将分析年份1967年、1978年、1987年、2005年、2012年、2019年实测断面深泓变化进行了对比,其变化见图2。
Figure 2. Variation chart of the channel line in different years
图2. 河段不同年份深泓线变化图
2. 计算方法
本次计算采取HEC-RAS (Hydrologic Engineering Center-River Analysis System)进行计算 [7],该软件是由美国陆军工程兵团水文工程中心开发的水面线计算软件包,适用于河道稳定和非稳定流一维水力计算,可进行各种涉水建筑物的水面线分析计算 [8]。
本次主要计算固定洪峰情况下的沿线水面线,因此采取河道稳定流进行一维计算。计算原理为伯努利能量方程并考虑流速水头损失,采用试算法求解 [9]。
计算公式为:
(1)
式中:Z1、Z2为断面1、2的水位(m);V1、V2为断面1、2的流速(m/s);α1、α2为断面1、2的动能校正系数;Δhf、Δhj为沿程水头损失与局部水头损失。
3. 计算边界条件
1) 河道断面选择
根据沿线水利工程修建时间以及发生洪水的统计,本次从1967年至2019年选择1967年、1978年、1987年、2005年、2012年、2019年等6年的实测大断面作为断面进行分析,整个河段布设18个断面,断面分布图见图1。
2) 水位流量选择
河段起点建有王家营水位站、终点建有老河口水位站,区间建有黄家港水文站,河段之间无大的支流汇入。因此,进口边界流量选择黄家港水文站实测流量成果,出口采用对应时段老河口水文站实测成果作为边界。按照洪峰的大小分别选择1964年洪水、1974年洪水、1983年洪水、2005年洪水和2011年洪水,洪峰流量涵盖实测中高水成果。边界成果见表1。
Table 1. The boundary conditions of water level and flow discharge calculation
表1. 计算水位流量边界条件
3) 河道糙率选择
河道糙率是河流阻力的一个综合性参数,其值的大小与河流的水深、比降、河床质组成、断面几何形状、水流平面形态等诸多水力因素有关,为此,张秉文 [10] 通过对天然河道糙率影响因素的定性分析,给出了单式断面、复式断面、大小河道及滩地漫流情况下的经验糙率取值;马玉红 [11]、方神光 [12] 等针对天然河道糙率影响因素及特点,提出了用曼宁公式反求糙率对河段及资料选样的要求,明确了确定糙率的基本原则;腾凯 [13] 针对目前率定河道糙率在实测洪水资料利用、控制断面水位流量关系曲线建立、复核洪水选用等方面存在的问题提出了相应对策及改进措施;张红武 [14] 通过对冲积河流糙率的计算提出了新的计算方法;王万站 [15]、张世安提出影响糙率的主要因素;王克仁 [16] 等对天然河道糙率的影响因素与分析方法进行了讨论。本次分别选择其推荐的常用曼宁公式和张红武提出的新的计算方法分别计算,最终选择合适的参数。
曼宁公式计算公式如下:
(2)
(3)
(4)
式中:Q为流量;A为过水断面面积;R为水力半径;Sf水力坡度;n为糙率系数;χ为湿周。
通过对王家营水位站、黄家港水文站、老河口水位站实测资料和比降分析计算,其从1954年至今实测资料反算河道综合糙率见图3。
参考张红武文献提出的公式 [14] 如下:
(5)
式中:n为河道糙率,Fr为傅汝德数。按照上式计算得到不同量级河道糙率成果见图4。
通过图3和图4对比可以看出,随着流量的增加河道糙率均呈现减小的趋势,这与基本规律是相符的。河道糙率在10,000 m3/s以下时同一量级洪水糙率均波动变化较大,这主要原因是河道的冲刷加剧了河段的河床粗化,低水过水断面主要为主河槽,同时由于2005年以下王甫洲水利枢纽建成后对于低水的顶托影响。但中高水时河道糙率基本稳定,受河道冲刷影响相对较小。
Figure 3. Calculation of flow roughness diagram using Manning formula
图3. 曼宁公式计算流量糙率关系图
Figure 4. Calculation of flow roughness diagram using Zhang Hongwu’s empirical formula
图4. 张红武经验公式计算流量糙率关系图
对图3和图4对比分析,张红武经验公式计算出的成果在10,000 m3/s以上流量时河床糙率基本位于0.02以下,曼宁公式计算成果主要在0.018~0.04之间,高水时期曼宁公式计算成果偏大。通过对两个糙率带入计算与实测值对比分析,曼宁公式计算成果与实际实测成果更为接近。因此,本次采取曼宁公式计算的河道综合糙率。
通过两种方法计算对比,其均表现为10,000 m3/s以上流量时,同一量级洪水时,河道综合糙率基本稳定,不同量级洪水河道糙率不一致且呈现随流量增大而减小趋势。范北林 [17] 等在南水北调中线工程调水后对汉江中下游河势的影响中河道糙率选择同一流量对应不同组水位的方法,求得不同流下不同组水位时的糙率,本次参照这种方法在河道糙率同一量级洪水选择同一个值进行分析计算,成果见表2。
4. 计算结果
按照上述边界条件和断面进行分析计算,计算得到同一量级洪水条件下,不同年份断面水位变化成果。通过计算对比分析,其变化规律基本一致,本次仅列出1964年洪水成果,成果见表3。
Table 2. Analysis of comprehensive roughness n value
表2. 分析工况综合糙率取值表
Table 3. Results of water surface lines in different flood years in 1964
表3. 1964年洪水不同年份断面水面线成果
5. 对比分析
通过对1964年、1974年、1983年、2005年、2011年洪水不同年份断面进行对比分析,其规律基本一致。本次以其发生的最大洪水1964年洪水为例进行分析。
当遭遇1964年洪水时,河道水位不同年份断面均存在一定的降低,其不同年份断面水面线成果见图5。
由图5可知,1964年洪水在1967年断面水位最高,在1978年断面时水位出现了明显的下降,最大下降幅度达到0.79 m,平均下降幅度0.61 m,该段时间内河道的冲刷明显,说明河道的冲刷对水位的降低产生了显著
Figure 5. Results of water surface lines at different sections of the 1964 flood
图5. 1964年洪水不同断面水面线成果图
的作用。从1978年断面至1987年断面,河道水位变化幅度最大为0.57 m,平均下降幅度0.26 m,水位变化幅度出现了明显的减小,这与该段时间冲刷减弱是吻合的。从1987年到2005年断面,水位最大变幅0.41 m,平均下降幅度0.28 m,水位变化幅度基本稳定,这与河道基本趋于稳定是吻合的,其最大水位变化逐步减小。
从2005年至2016年河道冲刷量达到4979万m3,河段河道冲刷已基本达到平衡,主要原因为人工采砂造成 [3]。在2016年~2019年,针对河道采砂进行了整治限制,河道冲刷量仅为601万m3,这与河道基本达到冲淤平衡是相符的。
从2005年到2012年,水位最大变幅0.52 m,平均下降幅度0.33 m,这与该段时间河道采砂影响有关,该段时间采砂量较水库建后时期增加一倍,这是水位变幅有所增长的主要原因。2012年至2019年断面,水位最大变化幅度0.45 m,平均下降幅度0.19 m,水位最大变幅主要是2012年~2016年人工采砂影响造成的。
综上,在丹江口水库蓄水后,下游河道的冲刷对水位降低产生了明显的作用,初期冲刷强度大,对整个河道的过洪产生明显影响,对河道的行洪有利;随着时间推移,下游河道的演变逐渐趋于向平衡方向发展,断面冲淤变化对水位影响逐步缩小,预计当下游河道冲淤达到平衡时,河道水位将逐渐趋于稳定。在2005年~2012年期间,河道已达到冲淤平衡,但是人为采砂影响,河道的冲刷量接近建库初期,因此水位仍有一定幅度的变化。该结论与郭海晋 [2] 对汉江雅口枢纽坝下河道冲淤及水面线变化数值模拟结论雅口枢纽坝下河道水面线受河床冲刷下切而降低是基本一致的。
6. 结论
1) 河道综合糙率受河道流量影响显著,在中高水时表现的更为明显。
2) 分析河段河道受丹江口水利工程影响,河段整体呈现冲刷状态,在建库初期冲刷现象十分明显,河道的冲刷对同频率洪水水位有明显的降低作用,运行初期平均降幅达到0.61 m,近几年平均降幅为0.19 m,随着时间的推移,影响逐步趋于稳定。
3) 对河段防洪标准1964年同大洪水,现状断面情况下水位较1964年实际洪水位有明显的降低,平均降幅达到1.6 m,采取1964年洪水水位作为河道设计洪水水位可以抵御超1964年同大洪水。河道中高洪水位均有不同程度的降低。
4) 目前河段已基本达到冲淤平衡,在对采砂后河段逐步达到新的平衡后,河道水位将趋于稳定。
5) 在实际工程中,当河道尚未达到冲淤平衡时,采取冲刷时段断面进行水面线计算将会造成计算水位偏低,不利于工程的安全。
6) 国内对河道冲淤变化后水位变化研究相对较少,本次仅在宋涛 [4]、苏程佳 [18] 相关文章中找到相关参考。本文分析丹江口水库建成后对下游影响,考虑2005年~2016年河道受人为采砂影响严重,河道水位是否受采砂影响有待进一步的验证。
7) 随着丹江口水库的加高后运行,丹江口水库对河道泥沙的拦蓄将进一步加大,河道是否会发生新的冲刷有待监测,河道水位将会随着河道的冲刷变化而发生新的改变。
参考文献