水-草-土耦合边坡流场特性研究_陈彬.pdf

第 卷 第 期 年 月人 民 长 江 ，收稿日期：作者简介：陈 彬，男，高级工程师，博士，研究方向为水土保持与边坡开挖。：通信作者：姚锡伟，男，副教授，博士，研究方向为水土保持与边坡开挖。：文章编号：（）引用本文：陈彬，姚锡伟，毛朝勋 水 草 土耦合边坡流场特性研究 人民长江，（）：水 草 土耦合边坡流场特性研究陈彬，姚 锡 伟，毛 朝 勋（内江职业技术学院，四川 内江）摘要：植被的存在会改变天然边坡的流场特性，其对水流产生的阻力较为复杂。将含植被边坡划分为水层、草层和土层，用 方程描述水层的流体运动，用 空隙理论描述草层和土层中的流体运动，在切向和法向的速度和应力保持连续性的条件下，推导出含植被边坡立面二维情形下沿坡度方向的流速沿水深分布的解析解。研究表明：当草层孔隙度较大时，地面下水流流动不可忽略；随着草径的增大，流速逐渐变小，最大剪应力由地面过渡到草顶面。研究成果可为土壤侵蚀研究和水土流失防治提供参考。关 键 词：流场特性；植被边坡；流速分布；孔隙度；剪应力；土壤侵蚀中图法分类号：文献标志码：引 言随着绿色生态理念的深入人心，明渠或边坡中的植被流引起了大量研究者的关注。当地表水流经植被区域时，植被按水流条件主要分为淹没和非淹没两种类型。宋滢汀等通过玻璃棒模拟刚性植被，采用三维激光测速仪对含植被明渠的水流流速和紊动强度进行测量，得到了含刚性植被明渠水流特性。陈甜甜等将植被水流分为树干层、树冠层和自由水层，在 涡流黏度理论和经典对数律流速公式的基础上，推求出淹没植被垂向流速分布解析解，并设计相应试验验证解析解的合理性。韩宇等通过随机森林、神经网络、支持向量回归 和 神经网络等机器学习算法研究了明渠出露刚性植被流速分布，发现模拟结果与试验结果吻合良好。任姗等利用人造草模拟柔性植被，通过 技术研究了柔性淹没植被的垂向流速、应力和紊流强度随曼宁系数的变化特性，发现应力的峰值出现在植被顶部，并与植被密度呈正相关。此外，一些研究者也通过数值模拟的方式研究含植被明渠或边坡中水流特性。采用大尺度格子玻尔兹曼方法分别模拟了刚性和柔性植被下明渠垂向流速分布。李雅琪等和丁雪采用 模拟含刚性植被的明渠流动特性。由上述分析可知，对于含植被渠道内水流流动特性的研究大多是基于物理试验和数值模拟，并且研究流体对象都是地表上的水流，而忽视了地表下的水流。本文提出一个含植被边坡的流场计算模型，模型主要分为 层：地表上有水层和草层；地表下为土层。联立 方程和 空隙理论，推导出含刚性植被的垂向流速分布解析解。理论分析图 描绘了二维表面流在恒定情况下经过植被区的场景。在图中，很明显有 层区域：水层、草层和土壤层。液体是假设为黏性的、不可压缩的均匀流，草和土层假定为均质的多孔介质，而在孔隙中运动的水视为 第 期 陈 彬，等：水 草 土耦合边坡流场特性研究图 含植被边坡示意 多孔介质流。所有的物理变量和参数在一个代表截面上取平均基本体积值（），从而在多孔介质内流动的流体可以用偏微分方程来描述。此外，地表流和地下流假设流动已达到稳定状态。当土层达到饱和状态时，流体的垂直速度可以忽略不计。由于流动长度比流动深度大得多，意味着切向导数比法向导数小，所以切向导数可以被忽略。基本控制方程基于上述假设，如果水层较浅，大部分冲击动量将转移到下垫层土壤中。如果忽略雨滴向水层的动量传递，同时忽略降雨在水流表层产生的水波，则 方程可以简化为 方向：，（）方向：（）式中：是水层流体在 方向（切向）上的速度，；是水层（即层）的流体压力，；和 分别为流体的密度和动力黏滞系数，（）；是重力加速度，；为坡脚，；（）为水层厚度，；是草的高度，；（）为水深，；、为图 所示坐标。和 简化了 的孔隙弹性理论，以模拟植被覆盖的倾斜地面上的低雷诺数流动，尺度多孔介质内部孔隙水的动量方程可以表示为 方向：（）方向：（）式中：为孔隙水在 层 方向（切向）上的速度，；为 层孔隙介质内部孔隙水压力，；为水的动力黏滞系数，（）；为 层孔隙度；为 层的本征渗透率，；下标 表示草层（）；下标 表示土层（）；为土层厚度，。边界条件在无风力、无表面水波、恒定流和使用表面压力的假设下，本研究含有以下 种边界条件：（）在自由液面处（）。流体应力在法向（方向）上保持连续：（）流体应力在切向（方向）上保持连续：（）（）在水层和草层的交界面处（）。流速在切向（方向）上保持连续：（）流体应力在切向（方向）上保持连续：（）流体应力在法向（方向）上保持连续：（）（）在草层和土壤层的交界面处（）。流速在切向（方向）上保持连续：（）流体应力在切向（方向）上保持连续：（）流体应力在法向（方向）上保持连续：（）（）在土壤层底部（）。由于无滑移条件，切向（方向）上的流速为零：（）（）在边界层上部（）。（）（）在边界层下部（，是沿坡地测量的坡度长度）。边界层下部水头零深度梯度条件：（）解析解通过求解式（），（）和（），然后将结果代入边界条件式（），（）和（），压力分布解析解如下：()，（）()，（）()，（）人 民 长 江 年式中：和 分别是水层和草层的高度；为坡脚。然后将压力分布代入式（），（）和（），可得到流速分布解析解表达式：()，（），（）()，（）式中：是待求系数。将式（）（）代入式（）（）和（），可求出 的表达式（，）。()()（）()()()()()()()()()()()()()()()()()()()（）()()()()()()()()()（）()()（）（）()()()()()（）()()()()()()()()()（）()()()()()()()()()（）分析与讨论下文中取水的黏度 ，水的密度 ，重力加速度 ，坡长 。粗砂的孔隙度和渗透系数分别为 ，。草层密度是由草层孔隙度决定的。植物直径（）与孔隙度（）之间的关系式及本征渗透率（）与孔隙率的关系式为 （）()（）为了验证该解的有效性，笔者发现推导的解析解在无草地面上的流速分布可以退化为 等提出的解的形式（见图）。换句话说，当坡面无植被时，等的结果是本文的一个特例。图 中的其他曲线描述了不同草密度下的速度分布。当草密度增加时，草对流体的阻力会增加，水位上升，流速降低，但草层和土壤层流速分布变得更加均匀。草的基本参数如下：，土层厚度 。由于自然界中草的直径大小不一，且不同的直径会对坡面流体产生不同的阻 第 期 陈 彬，等：水 草 土耦合边坡流场特性研究图 不同孔隙度下淹没草模型速度分布（）（）力，所以图 中分别选取 为，和 进行讨论。剪切应力可由 得到。如图 所示，其中 和 分别为 时的最大速度和最大剪应力。从图（）中，可以发现垂直速度分布中存在拐点，这与一般认知一致。从图（）中可知，最大剪应力存在于水草层界面或地面上，其取决于草的分布密度。此外，参考，流体通过柱状纤维多孔介质时单位长度的阻力可表示为 （）（）式中：是平均流速。水流通过草层的拖曳力变化见图。图（）为各位置的平均速度分布，为坡比 时的最大平均速度。平均速度随着斜率的增加而增加，这与传统的曼宁公式的结论相同。由于式（）中通过草层的阻力与平均速度成正比，因此阻力的无量纲分布与平均速度分布相似，如图（）所示，其中，时，为最大阻力。值得注意的是，随着坡度的增加，阻力分布的偏移量也在增加。土层厚度对模型流速分布的影响见图，可知水土交界面附近有轻微的差异，总体上土层厚度对模型整体流速分布影响较小。图 草直径对模型结果的影响 图 坡度变化对模型结果的影响 草的直径对阻力的影响如图 所示。流经草层的水流所受到的阻力随草径的增大而减小，而随草层孔 人 民 长 江 年隙度的增大而增大。这是因为相同情况下，草的直径越大孔隙就会越小。平均速度会随着孔隙度的减少而降低，从而阻力会降低。图 土体厚度对模型速度分布的影响（）（）图 草的直径变化对阻力的影响 结 论（）本文研究了水层、草层和土壤层的切向流速分布。采用 方程描述流体在水层的运动，通过 空隙理论描述流体在草层和土壤层中的运动。在速度和应力连续条件下推导出流速切向分布解析解。（）通过与前人无草条件下的流速分布进行对比，验证了模型的有效性。研究发现，土层厚度变化对模型流速分布影响不大，流速分布主要受坡度的影响。（）最大剪应力存在于水 草层交界面或草 土层界面上。随着草径的增大，最大剪应力主要出现在水 草层交界面处，这将有利于防治土壤侵蚀。参考文献：槐文信，钟娅，杨中华 明渠漂浮植被水流内部能量损失和传递规律研究 水利学报，（）：，米云彤，王丹，蔡暾，等 植被群河道水流紊动特性研究 中国农村水利水电，（）：焦军丽，王丹，李文奇，等 植被群密度对河道水流紊动特性影响的试验研究 水力发电学报，（）：，：宋滢汀，景何仿，王维红，等 含多重刚性植被明渠水流特性实验研究 水力发电，（）：陈甜甜，赵明登，张颖 淹没树状植被水流速度分布规律研究武汉大学学报（工学版），（）：，韩宇，刘家备，穆同，等 基于机器学习的含植被水流流速分布研究 中国农村水利水电，（）：任姗，冯民权 含柔性植被明渠水流水力特性的试验研究 水资源与水工程学报，（）：，江雅赛，张景新 有限植被群水流尾流场特征的数值模拟研究 水动力学研究与进展（辑），（）：丁锐，李坤芳，黄尔，等 植被作用下明渠交汇水流特性数值模拟研究 中国农村水利水电，（）：，周睿，程永光，吴家阳，等 用浸没边界 格子 方法模拟双层刚性植被明渠水流特性 水动力学研究与进展（辑），（）：，（）：李雅琪，刘善均，刘超，等 关于刚性植被尾流区两种数值模拟方法对比 水电能源科学，（）：丁雪 含刚性沉水植物明渠水流结构试验和数值模拟研究扬州：扬州大学，：，：，（）：，（）：，；（）：，第 期 陈 彬，等：水 草 土耦合边坡流场特性研究 ，：（编辑：郑 毅），（，）：，：；（上接第 页），（，（），；，）：，：；