长江上游水利枢纽下泄流量对兰叙段航道水位影响研究.docVIP

长江上游水利枢纽下泄流量对兰叙段航道水位影响研究 　　摘 要：该文根据近年来该河段的实测资料对其进行河流演变分析，以及对其典型河段进行数学模型的建立和分析，以此来分析向家坝水电站的排洪下泄对该河段的流速及水位的影响。 　　关键字：排洪下泄 数学模型 流速 　　中图分类号：O 319.56 文献标志码：A 文章编号：1000?C6915(2008)01?C0001?C03 　　 　　 　　一、概述 　　向家坝水电站位于四川省宜宾市和云南省水富县交界处的金沙江干流下游，是金沙江梯级开发规划中的最末一级电站，该工程是以发电为主，是同时具有防洪、拦沙及改善航运、兼顾灌溉等综合作用的特大型水利枢纽工程。但向家坝水电站由于其电站日调节，大坝泄洪以及电站事故等工况引起的非恒定流将改变下游河道的天然水流条件，为了研究非恒定流对下游河段通航水流条件的影响，通过本项目的研究，探究支流金沙江、岷江水利枢纽下泄流量对本河段航道水位演变规律以及典型碍航水道特征，为科学的制定航道维护计划，充分利用水利枢纽蓄放水对航道水位调节作用，提供一定的技术保障。 　　二、河床演变分析 　　所研究的河段上起向家坝，下至泸州沱江口，河段内滩险众多，其浅滩成因主要与具体的河型、河势有关。总体来讲，研究河段可主要分为两大段，一段为向家坝至宜宾河段，该河段段为金沙江，地形地貌上看属深切河谷，河道狭窄，河段内宽谷与峡谷相间，纵向上也是深潭与浅脊相间，该河段内浅滩主要位于放宽段内，其特点是河道内礁石、石盘较多，枯水期水流散乱，航行条件差。第二段是宜宾至泸州段，该段地貌上属低山丘陵带，河道在此段内汇入众多支流，主要有岷江、赤水河等，该段内浅滩主要位于河道放宽段，且其成因多受局部地形影响或受支流入汇影响形成，从总体上看，研究河段河岸、河床多由卵石、基岩组成，抗冲性较强，河势是基本稳定的，但由于研究河段地形复杂，水流流速较大，在局部河段受河型、河势及水沙条件影响较为复杂，但年际冲淤变化可基本维持稳定。 　　三、数学模型的建立 　　（一）一维水沙模型 　　一维水流数学模型常用来模拟天然河道时间和空间尺度较大的水流的运动。在数值格式的研究中, 近期具有代表性的研究成果是借鉴计算空气动力学中的优秀算法, 获得的精度高、守恒性好的多种计算格式。 　　1.一维水流的控制方程[1] 　　明渠或天然河流被考虑作为一维非恒定流动时，基于数学模型为圣维南方程组，以流量Q(x, t) 和水位Z(x, t) 为未知变量，它的方程组为： 　　连续方程： （1） 　　动量方程： （2） 　　其中A为过水面积；Q为流量；Z为水位；Sf为河床摩降；g为重力加速度， ， 　　（其中：n为Manning 糙率系数；R 为水力半径。） 　　则上述方程可写为向量形式： （3） 　　2.一维方程组的求解 　　采用清华大学的钟德钰[2]最近提出来的一种新的基于有限体积法的TVD的数值格式求解对式(3)采用采用一阶显格式有限体积法TVD离散,即： 　　 　　但当河道断面沿程变化剧烈，水位变化幅度较大，河床冲淤幅度大的时候，为避免出现动量方程计算流量时常出现锯齿状振荡现象，本文模型对谭维炎[3]等的方法进行了改进，即将流量与水位、含沙量交错布置,如下图1所示： 　　 　　 　　 　　以此动量守恒方程中水面坡度一项自然离散为： 　　 ， 式中A可用(Aj+Aj+1)/2计算。 　　3.边界条件 　　边界条件一般可以分为三种形式给出：1.已知水位随时间的变化，2.已知流量随时间的变化，3.已知流量随水位的变化。由于采用交错网格布置，使得对于进出口的边界布置非常简单，在上游直接根据时间流量关系来进行代入第一个流量。下游的水位根据计算额流量，然后根据流量水位关系直接插值得到下游边界的水位。 　　（二）二维水流模型[4.5] 　　1.二维水流的控制方程及求解： 　　 　　以上各式中，t为时间，x,y为空间坐标，z为水位， h为水深，h=z-zb，zb为河床高程，υ,ν分别为沿x,y方向的垂线平均流速，g为重力加速度，n为糙率系数， νt为湍流动能扩散系数。 　　利用有限元单元法来计算控制方程的数学解，基本求解过程为： 　　①将控制方程分别在时间和空间上进行离散。其中时间离散采用差分法，空间离散采用有限元法。 　　②通过伽辽金加权余量法将控制方程从偏微分方程转化为代数方程组； 　　根据边界条件和初始条件数值求解代数方程组； 　　2.二维边界条件 　　在二维的计算河段的上下边界处，根据一维计算的水位和流量进行插值到二维的网格中，在二维的模型中边界网格直接用已知的数据来得到。 　　四、计算范围及网格的布置 　　由于从向家坝至泸州的全河流的资料不足，主要的地形资料