天然河道一维数模中综合糙率确定方法.docVIP

天然河道一维数模中综合糙率确定方法 　　摘 要：天然河道断面形态沿程变化比较大，水力计算时需要考虑局部阻力和沿程阻力双重影响。根据能量守恒原理，提出了天然河道一维水力计算时糙率系数的确定方法，并通过实验对该方法进行了验证，取得了预期的结果。该方法减少了率定糙率时的工作量，对缺乏历史观测资料的地区同样适用。 　　关键词：综合糙率；天然河道；一维数值模拟；水头损失 　　中图分类号：TV12 文献标识码：A 　　前言 　　综合糙率系数是一个反映河床边壁粗糙程度、边界形状以及其他因素的对水流影响的一个系数，是进行水力计算、河流模拟研究的重要参数，它的正确与否直接关系着数值模型的精度。天然河流，特别是山区型河流，往往有许多弯道、断面形状突变段，水流结构沿程变化比较大，此时如果仅考虑由沿程阻力造成的能量损失，就会产生较大的误差。本文通过研究一维水力计算中糙率系数的确定方法，导出了利用能量方程来确定综合糙率的计算方法，并进行了验证。 　　1 概述 　　确定河流综合糙率系数的方法大致归纳起来可分为2种： 　　1.1 利用实测资料对糙率进行率定 　　此种方法根据实测的流量、水位湿周等水力要素，利用最优化原理，定出使计算所得的相关水力要素值与实测值最接近的糙率值，这种方法确定的糙率能够综合反映河流的沿程和局部阻力特性，因此应用较广。但它是利用实测资料反演得到糙率，不可避免地受到观测误差的影响，当研究河段较多或缺乏完整的观测数据时，参数取值具有很大的不确定性。 　　1.2 利用河床的壁面特征及状况确定的经验值 　　此种情况确定的糙率值，仅反映由于切应力的作用水流产生的阻力，当计算河段的距离较短，局部能量损失远大于沿程能量损失时，由该方法确定的糙率值不能真实反映河流的阻力特征。 　　2 糙率计算方法 　　天然河道的断面形状沿程变化比较明显，河流结构沿程变化比较大，进行一维水力计算时，可以将研究河段分为多个断面形状近似不变的子河段，在每个子河段内综合考虑沿程阻力和局部阻力的影响，确定出能正确反映该河段阻力特性的综合糙率。 　　假设研究子河段长为L，第i个子河段的长度为Li，在该子河段上下游运用能量方程： 　　（1） 　　式中，、分别为上下游水位；、为动能修正系数；、分别为上下游断面流速。式中沿程水头损失可以表示为： 　　（2） 　　式中，n为河道糙率，设该子河段总水头损失为，则： 　　（3） 　　将该子河段内的局部水头损失全部计入沿程水头损失中，将表示为： 　　（4） 　　式中反映了在沿程阻力和局部阻力的综合作用下河道的阻力特性。整理后得到综合糙率为 　　（5） 　　3 算例 　　为了测试本文提出的确定综合糙率的方法是否可行以及研究将它运用到水力数值模拟时可能产生的影响，本文选择了一条具有3处局部能量损失的明渠河流作为测试的算例，首先根据式（5）计算综合糙率，并运用该综合糙率进行一维水力计算，然后运用伯努利方程推算水面线，并将2者的计算结果进行比较。计算时，为了保证伯努利方程能正确运用，假设水流为稳态。 　　此次所选的计算明渠为梯形断面，底宽b为50m，边坡系数m为1.5，河底纵坡I为2×10-5，糙率n为0.022，入流流量为300m3/s，在2000m处有两处断面急变段，在4000m处有一急弯，明渠出口处渠底高程为0m，由于出口处距离急变流段足够远，经验证此处水流流态为均匀流，该断面处的正常水深为7.40m。为了更好地表示出局部阻力对水流的影响效果，此次数值模拟进行了考虑局部水头损失和不计局部水头损失2种情况下明渠水面线的计算。本次计算中，试验渠道被31个断面分为30个子河段，断面间距为200m。计算结果见图1，由图1可以看出，在同样的水力要素下，考虑局部水头损失和不计局部水头损失对河道的水位计算有较大的影响。 　　图2 调整后数值模拟水位曲线 　　为了与数值模拟方法对比，运用伯努利方程推算水面线时，断面间距也为200m，同样分为考虑局部水头损失和不计局部水头损失2种情况，计算结果见图1。2种计算方法的结果在不计局部水头损失时水面曲线基本一致，而在考虑局部水头损失时，数值模拟计算的水位曲线中局部阻力引起的能量损失均匀分配在计算河段内（该情况不影响所有断面的水位值）。为了解决以上问题，更精确的模拟出发生局部水头损失的断面附近水流的变化情况，将断面间距缩小为20m重新进行一维数值计算，调整后的模拟结果见图2，由图2可以看出，当断面间距比较小时，可以明确观察出发生局部水头损失处的水位变化情况，此结果与伯努利方程推算的水位沿程变化曲线一致。因此，本文提出的确定综合糙率的方法不但可以对一般河段进行糙率估计，也可以通过在敏感河段增加断面的方法对复杂河段进行精细模拟。 　　4 结语 　　对于存在弯道、