明渠漂浮植被水流内部能量损失和传递规律研究.docx

? ? 明渠漂浮植被水流内部能量损失和传递规律研究 ? ? 槐文信，钟 娅，杨中华 （武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 430072） 1 研究背景 漂浮植被在生态河道中很常见，相比于普通的淹没植被和非淹没植被，漂浮植被漂在水面上，并不扎根于河底，植被层和底部河床之间有自由的水流区域［1］。常见的漂浮植被有浮萍科、浮叶类和睡莲水葫芦等［2］。江河湖泊中由于水体污染及富营养化，漂浮植被肆意“疯长”，对水流产生过大阻力，导致河道过水能力减小，造成缓慢流动的明渠水流［3］，而漂浮植被其本身也没有生存天敌，不会受到自然条件的阻碍，会迅速覆盖整个水面，对行洪排涝和某些水生生物繁衍产生了不利影响。但另一方面，漂浮植被为水生生物提供了好的栖息地，生态工程上广泛地应用漂浮植被净化污水［4］和富营养化的水体［5］。对漂浮植被如何适度把握、趋利避害以及化害为利成为关键，因此本文对漂浮植被水流内部能量的研究是非常有必要且有意义的。 对淹没及非淹没植被的研究目前比较广泛，比如研究带有沉水植被、坡面植被的矩形或者复式断面上流速分布和水头损失等特性［6-8］。漂浮植被对水环境的生态影响方面的研究比较成熟，但很少有学者研究明渠漂浮植被水流内部的水动力学特性。近年来一些学者利用二维［9］和三维［10］水动力模型、水流分区及混合长度理论［11］还有数值模拟［12］和实验测量［13］等方法研究漂浮植被水流的流速分布和水流结构等特性［14］。漂浮植被的顶部是自由水面，植被区水流受到的作用力主要有重力、紊动切应力和植被拖曳力，植被以下的流区除了会受到植被的影响之外，还会受到底部河床的制约，无植被的流区受力主要有重力和紊动切应力。在研究淹没植被和非淹没植被的拖曳力时，由于底部河床阻力远小于植被拖曳力，所以在底部植被区主要考虑植被拖曳力［15-16］。但在漂浮植被水流中，植被拖曳力和河床阻力都对水流有较大的影响，则淹没植被及非淹没植被水流的理论及经验方法对漂浮植被并不适用，所以很有必要研究漂浮植被对水流结构的影响。 1945年Bakhmeteff针对明渠均匀流，研究过水断面垂向上的能量分布，得出断面上能量的传递机理，分析了能量的耗散过程，并且给出了断面上能量提供、能量损失和能量传递的表达式，建立了这三者之间的关系［17］。结果表明在水流内部的任一点，当地能量提供等于能量损失和能量传递之和，整个断面能量传递总和为零。基于这个理论，明渠淹没植被水流的断面能量传递及转化也得到了分析，结果显示能量损失主要集中在植被区，能量提供亦等于能量损失和能量传递的总和［18］。本文同样基于Bakhmeteff断面能量平衡理论，力求揭示明渠漂浮植被水流的断面能量传递规律，并利用Plew［13］的水槽实验结果进行验证。 2 理论分析 2.1 水流提供能量 对图1所示的有漂浮植被的明渠恒定均匀流，根据总流能量方程可得任意两个断面之间的水头损失为hf=E1-E2，其中E1、E2为任意两个断面的总水头，因此明渠均匀流的水力坡度S=hf/?l=(E1-E2)/?l，其中?l为两个断面沿水流方向的距离。均匀流中水力坡度S就是单位重量液体在单位流程内水流提供的能量，单位体积液体在单位流程内水流提供的能量为γS，其中γ为水的重度。将γS乘上流速u便得到单位体积液体单位时间内在该位置上提供的能量Wb： W0=γSU即为断面平均的单位液体在单位时间内提供的能量，或称为断面平均的能量提供率，其中U为断面平均流速。 二维恒定均匀紊流中，沿水流方向带有植被的平均动量方程［19］为： 式中：z为垂向坐标；τ=-ρu′w′为时间平均和空间平均的雷诺切应力；ρ为水的密度；u′为纵向脉动流速；w′为垂向脉动流速；F为由植被引起的拖曳力，在植被存在的区域通常被参数化为一个关于流速u、拖曳力系数Cd和植被密度a（单位体积植被的投影面积）的二次函数［20］： 在有植被的区域，植被拖曳力由拖曳力系数、植被密度以及流速决定，然而在无植被的区域拖曳力为零。漂浮植被形态、坐标系及流速分布情况如图1中所示，其中x和z分别为纵向和垂向坐标，H为水深，hc为植被高度，hg为植被底部至河床底部的距离。 图1 有漂浮植被的明渠恒定均匀流示意图 将式（2）代入式（1）可以得到水流能量提供率： 由此可知，能量提供率由该点流速、切应力梯度和拖曳力决定。对Wb求任一点z至水面H的积分，即可得到自z至H流区内的能量提供率的累积值，亦即单位时间内该流区单位水体所提供的总能量。 从式（5）可以看出漂浮植被水流的能量提供累积值与无植被水流情况下是一致的，则表示从水底到水面整个流区的总能量提供。 2.2 水流损失能量 明渠漂浮植被水流中的能量损失大致分为两部分：（1）水体由于切应力梯度的存在会产生变形，水体克服水流阻力损失的能量；（2