感潮河道河口水闸闸外水位震动现象的初步探讨.pdfVIP

大泖港等三大支流以及太浦河，主要为太湖流域杭嘉湖地区的运东片，区域总面积约为2000km2。 2．2模型基本方程 Venant方程： 描述水流在明渠中运动的一维非恒定流基本方程是Saint 连续方程：挲+掣：留 dX dt U Q=。 动舫程：鲁+掣+gA警+可gn2 式中：t，工为时、空变量；Z，Q，U分别是各断面的水位、流量和流速；A，B，R分别是各断面的 过水断面面积、水面宽和水力半径；口滩地系数；q为单位河长的均匀旁侧入流(包括降雨产汇流)； n为河道综合糙率系数；g为重力加速度。 由于SaintVenant方程组是一个二元一阶拟线性双曲型偏微分方程组，求得精确解非常困难，模 型采用四点隐式非线性直接差分法离散求数值解。其差分格式如下： 溯．M…j+1+M?+1一M二1一M? af 2△f ‘ 丝：堡!丝墨二丝竺!±!!二皇丛丝矗!二丝!! 孤 血 一==一 M=aM／+1+(1一O)M／ 式中：目是权重系数(O．5s秒51．0)：M为任何变量，如水位、流量等；下标代表断面号；上标代表计 算时刻。时段初的_『时刻为已知，时段末的一1时刻为待求。 2．3过闸流量计算 在闸坝泄流时，如果闸门开启度较大可能属于堰流，如果闸门开启度较小则可能属于孔流。一般 采用下列条件(宽顶堰式)判别堰流和孔流： 二≥0．65时，为堰流 月 二0．65时，为孔流 且 式中，e为闸门开启高度，日为闸上水头。 在堰流情况下，过闸流量可按水力学宽顶堰公式计算： 自由出流：Q=聊即玩3儿 淹没出流：Q=缎√29(乙一乙) 顶的水深(m)；Zo为闸底高程(m)；Zu为闸上游水位(m)；Zd为闸下游水位(m)。 在孔流情况下，过闸流量计算公式： Q=cr／．eB＼／29Ho 式中，盯为淹没系数，∥为孔流流量系数。 2．4模型率定与验证 水文资料和1997年11号台风暴潮水位资料等资料进行了模型率定与验证计算。率定计算成果表明水 位、流量过程与实测情况基本吻台良好(见罔1)。 蛳 _嚣 ，＼ f、 』I n ，、 ；：㈧㈧＼刊 瑚V 。 ” 、』 ＼』 ，产，上_ ：鲥～：= l二=Ⅳ，『Lql～ 圉1 I 997年1{号台风犟定计算黄浦公园站水位和吴凇断面流量过程 3 闸外水位震动现象及影响因素分析 对于挡潮闸方案，由于闸门门型的选择具有较太的不 确定性，本文分析计算考虑最为常见的直升川]型方案。 根据有关研宄闸址方案可能有多个方寨．奉文主要阻距离 河口约58km的闸址方案为例进行分析计算，考虑该闸址- 处河道断面宽度、通航宽度和高度等因素．计算中闸孔净 宽按500m考虑。计算时间步长取1分钟。 在黄浦江吴淞口出现1000年一通设计潮位情况下， 关闸历时取30丹钟，考虑整个关闸过程，计算成果表明． 瞳2不同关闸水位情I兄下闸外计算水位过程 闸门关闭后闸外水位有明显的壅商．并且出现明显 的水位振荡现象。不同关闸水位情况下闸外汁算水位过程见圈2 由图4可见．闸外水位振荡与关闸水位密切相关，闸门 起共水位越高(即相对于高潮位关闸时间越晚)．闸外堆高 水位越高．水他振荡的幅度也越大。可见挡潮闸应根据预报 情况尽早关闸为宣。 为了分析边界条件对闸外水位振荡的影响，在黄浦汀曼 凇口外连接一段长扛河段进行分忻计算，长匝河段宽度远大 于黄浦江，可以看作是一个天然“调蓄水库”，边界条件外移 后完全可以认为边界对波浪的影响很小了。司算结果表明闸 外水位振荡仍然存在f见圉3)，并且可以发现圈5与圈2 囤3边界条件外移后计算的闸外水位过程 蓬本致。口f见闸外水位的振荡上要不是