三峡水库蓄水后重庆段水流状态分析.docxVIP

三峡水库蓄水后重庆段水流状态分析 三峡在防洪、蓄水和交通运输等方面发挥着显著的综合作用，对水库的生态和环境产生了重大影响。水库水位将同时建成至最终水平，并根据各阶段的蓄水要求逐步增加。2003年，水库水位将逐渐增加至135m。该项目开始充分发挥蓄水和通风的作用。2006年10月，水库开始运营初期的156m蓄水位。经过几年的早期蓄水，水库水位下降到最终正常蓄水位的175m。国内外专家对峡谷自身水流问题进行了详细的分析研究，并取得了丰富的研究成果。然而，由于多种复杂因素的影响，人们对水库储存175m后的水流状态研究较多，而对水库储存的平均距离水流预测不足。基于圣维南方程建立的三维水流模型，作者计算了库夏河段及其自然状态（水库前65.8m水位）的整体水流特征。然后，应用模型预测了库夏河段（水库前156m水位）和正式运营（水库前175m水位）后重庆段整体水流的特征。例如，整个河流速度的变化是整个过程的结果，横截面的变化规律是控制排放点和排放方法的重要参考。为了分析不同蓄水位在重庆段水流状态和生态环境中的影响，分析了不同蓄水位的影响。 1 数学模型 1.1 市域出口河流 一维流场(计算区域为重庆境内,长江干流寸滩断面开始到重庆市域出口约500余km长江干流江段.计算范围包括了长寿、涪陵、丰都、忠县、万州、云阳、奉节、巫山等城区江段.并且在计算范围内有乌江、大宁河、御临河等31条一级支流汇入. 1.2 水库公民水位公民n 分别计算了平水期常水位(坝前65.8 m),水库二期蓄水水位(坝前156 m)以及正常蓄水后11、12月高水位(坝前175 m)3个水位的水流参数. 1.3 起点流量和正常蓄积水 一维计算的起点为寸滩断面.常水位(坝前65.8 m)起点流量采用寸滩断面平水期多年90%保证率流量,二期蓄水高水位(坝前156 m)和正常蓄水高水位(坝前175 m)起点流量采用寸滩断面多年11～12月平均流量. 1.4 断面宽度的计算 一维水流计算模型以长江干流沿流向长度为空间坐标,在一定距离和支流汇入口处分段.段内流量、流速、水力坡降等参数保持一致.根据河床数字地形图提供的断面划分精度,计算江段长度1.0～2.0 km,计算范围共划分计算江段363段.计算江段起点断面和终点断面的水面宽度确定后,计算江段内各点的水面宽度在起点断面和终点断面的水面宽度之间线性变化(如图1). 1.5 1层速度计算方法 1.5.1 u3000过水面积及国家应力q的计算 模型基本方程采用一维圣维南方程组的基本方程: 连续方程:?A?t+?Q?x=q.(1)?A?t+?Q?x=q.(1) 动量方程: ?Q?t+??t(α?Q2A)+gA?Ζ?x+gASf=?Q?t+??t(α?Q2A)+gA?Z?x+gASf= q·V分, (2) Κ=CA√R, (3) Sf=Q2Κ2=Q2C2A2R=Q2LC2A3. (4) 式(1)-(4)中A为过水断面的过水面积;q为旁侧入流,若为出流,q0;B为水面宽;Q为流量;Z水位;g为重力加速度;Α为动量校正系数;L为湿周;Sf为摩阻比降;Sf用流量模数K和曼宁公式(或谢才系数C)计算. 在旁侧入流q的动量沿干流流向的分量相对较小,常忽略不计. 1.5.2 数值化河道断面结果 式(1)和式(2)全面考虑了计算江段的旁侧入流流量(包括集中入流量和单位长度入流量)和支流流量在沿干流流向方向上的动量分量.但这给方程组的求解造成了困难.因此,适当简化方程组以降低下一步计算的难度.文中采用如下基本假设:计算江段在一定距离处和支流汇合处分段,计算江段内无支流汇入,为恒定流;支流流量沿干流流向方向动量可忽略不计,动量校正系数α取为1.0.则方程组简化为: 连续方程:Q=A·V, (5) 动量方程:Q2=??x(1A)+gA?Ζ?x+gASf=0. (6) 式(5)中V为断面平均流速,其余符号意义同上. 将连续方程带入动量方程,且把式(3)、式(4)代入得 (1-gQgV3B)?V?x-gQgV3B?B?x+gV(V3LC2Q+Ι)=0, (7) 采用隐式差分格式求解式(7).采用Bi=√BiBi+1带入式得: Vi=Vi-1+[gQ1Vi=1BiBi+1(Bi+1-Bi)Δxi-gn2V73i=1L43iQ43i-