基于二维水流模型的南四湖调水工程水力计算.docxVIP

基于二维水流模型的南四湖调水工程水力计算 1 湖内输水流场计算 水生百合河从长江下游扬州附近提取长江。自江苏省汉江市以来，汉江市已从韩庄老运河和不老河改为山东省南四湖下游湖泊，通过二级水库站将水转移至下游湖泊，然后通过芦沟、邓楼、柳河航道监察桥、东平湖，从东平湖分为两条街道：山东省、山东省和天津市。东线输水渠道的90%是利用京杭大运河及与其平行的河道和湖泊,南四湖则是其重要的输水通道和调蓄湖泊之一［１］。南四湖位于山东省西南部济宁市,是南阳湖、独山湖、昭阳湖和微山湖四个由西北向东南相连而成湖泊的总称,流域面积31 700km２。南四湖内总面积1 266 km２,南北狭长126 km,东西宽5~25km,岸线总长311km。在湖中部最窄的湖腰建有二级闸坝水利枢纽工程,将其分成上级湖和下级湖。二级闸坝以北为上级湖,湖内面积606km２,二级闸坝以南为下级湖,湖内面积660km２。南四湖天然水流是由北向南经上级湖泄水入下级湖再入中运河,而南水北调东线南四湖输水是由下级湖南端进入、经二级坝泵站、由上级湖北端入梁济运河［２］。 目前,对于南四湖输水能力的水力计算一般仍采用一维明渠恒定非均匀渐变流的逐段试算法,其计算结果无法反映湖内水流的实际流动情况,水面曲线的计算误差难以估计,给南四湖输水航道开挖疏浚工程设计带来一定的困难［３－５］。关于浅水湖泊风生流的研究较多,马生伟等［６］采用迎风有限元方法模拟了太湖的风生流;张萍峰等［７］利用二维水动力学模型研究了滇池的风生流;武周虎等［８］研究了调水期南四湖上、下级湖在南风、东南风和北西风情景下的风生流流场;骆辉煌等［９］基于南四湖100m精度的DEM地形数据,采用平面二维水动力学模型研究了南四湖调水情景下不同风向的水动力特性,由于所采用的地形数据精度不能有效地反映南四湖航道和错综复杂的地形特征,忽视了湖内主航道的主要输水功能,其结果与实际情况有较大偏差。总之,平面二维水动力学模型的研究均未涉及采用浅水动力学二维数学模型研究湖泊输水航道开挖疏浚、水面曲线和特征水位等工程计算问题。 本文基于SMS(Surface Water Modeling System)地表水模拟系统创建适合南四湖深度平均平面二维水流数学模型体系,根据南水北调东线一期工程南四湖设计方案,进行南四湖输水航道开挖疏浚方案的优化模拟与计算,确定最佳开挖位置、长度和底宽,给出满足设计输水流量时的水面曲线和特征水位,开展湖内输水流场分析,以增强南四湖输水工程设计的可靠性、合理性和科学性。该研究成果已作为南四湖湖内开挖疏浚工程设计建设的科学依据和技术支持,对其它湖泊、水库类似工程的水力计算具有借鉴意义。 2 数学模型 2.1 南四湖洪水正常设计方案的阵风下水位分布 南四湖湖盆呈浅平形,湖水不深,一般情况下,除京杭大运河湖内航道外,大部分湖区水深不及1 m,最深处在微山岛以南,水深3m左右。基于此,采用深度平均的平面二维数学模型,相应的浅水动力学基本控制方程包括两个运动方程和一个连续方程,控制方程组为［１０］: 式中:h为水深;u、v分别为x、y方向的流速分量;x、y为正交坐标;ρ为水的密度;g为重力加速度;E为水平涡黏性系数张量;zｂ为湖底高程;n为糙率系数;ω 为地球自转角速度;φ为当地纬度;ζ为风应力系数;W为当地风速;Ψ 为风向与x方向的逆时针夹角。 南四湖调水正常设计方案的上、下级湖出湖口水位均为其相应的死水位,特殊干旱设计方案的上、下级湖出湖口水位则低于相应的死水位0.30 m。湖内地形情况复杂,水生植物种类繁多,各级航道、深槽、芦苇、湖草、莲藕、鱼池、湖田、河沟、庄台建筑物等纵横交错情况遍布湖区。湖内水体面积占湖区面积的45.53%以下［１１］,非航道区水深一般很浅或沼泽化,多为湿地植物保护区。受各种植物、芦苇荡分布的阻隔作用,非航道区糙率增大,风生大环流难以形成主导作用。文献在南四湖调水正常设计方案平均风速3.3m/s的条件下,模拟了NW、N、S、SE风向和无风情况下的上、下级湖调水流场和水面线,结果表明有风比无风情况下的上级湖控制点水位相差范围为-0.02~0.02m,下级湖控制点水位相差范围为0~0.02m,其水位差的相对影响误差均小于5.0%。因此,对南四湖调水设计方案来讲,风力的影响可以忽略;文献的研究同样说明了这一点,因此本文不考虑风力作用。 采用伽辽金有限元法求解上述方程组,借助于美国研制的SMS地表水模拟软件支持系统,以南四湖湖区万分之一地形图(2001年CAD版,标有实测苇草分布,1985年国家黄海高程基准,下同)和卫星遥感影像资料［１１］为基础,根据南四湖上、下级湖湖底地形、地势和岸线变化等边界特征,进行结点/网格划分、输入结点高程和单元介质类型、进行网格质量检查、SMS程序调试与运行、参数选取设置与模型