滩地刚性植物消浪作用特性研究.pdf

年 月水利学报 第 卷第 期文章编号：（）收稿日期：；网络首发日期：网络首发地址：?基金项目：广东省水利科技创新项目（）；广东省科技计划（）作者简介：张之琳（），高级工程师，主要从事水动力、水环境研究。：通信作者：黄本胜（），教授级高级工程师，主要从事水力学及河流动力学研究。：滩地刚性植物消浪作用特性研究张之琳，黄本胜，刘达，邱静，（广东省水利水电科学研究院，广东 广州 ；河口水利技术国家地方联合工程实验室，广东 广州 ；南方海洋科学与工程广东省实验室（珠海），广东 珠海 ）摘要：滩地植物有利于保障沿海地区的水安全，量化成熟期植物消浪作用是提升生态海岸韧性的基础。本文基于一般试验规律和两种波高消减模型，分析提出了新的滩地植物消浪系数计算方法，结合波高条件可计算林带消浪后波高，相关表达式由大量试验数据进行了验证，最后对比了由本文方法与规范方法计算的林带后波高。结果表明，阻尼因数及消浪系数能直观体现植物消浪作用特性，波浪传播方向上植物面积占横截面比例是影响植物消浪机理的关键无量纲参数。植物非淹没时，本文方法得到的林带后波高与规范方法接近；植物淹没时，本文方法能更准确描述植物消浪作用。相较于规范方法，本文方法考虑了植物淹没情况，直观揭示了林带密度、宽度、高度，植物平均直径，水深等参数的影响，适用于更多工况。基于植物消浪作用的模型可快速研判相关参数变化对植物消浪效果的影响，为生态海堤滩地刚性植物规划设计及数模研究提供科学和技术支撑。关键词：植物消浪；波高消减模型；消浪系数；拖曳力系数；物理模型试验中图分类号：文献标识码：?研究背景红树林等滩地植物的消浪作用实质上是波浪与植物相互作用导致波浪能量的耗散，而林带与水体相互作用复杂，量化成熟期滩地植物的消浪作用是提升生态海岸韧性、合理利用滩涂资源的基础 。目前针对植物消浪作用的研究方式主要有原型观测 ，、物模试验 及数模计算 ，其中，植物消浪系数表征波高在林带前后的变化，是量化植物消浪效果的关键系数。原型观测方面，由于观测地点及其地形地貌、植物类型及其种植结构的差异，不同学者观测到的植物消浪作用虽有所差异，但已知盐沼湿地的消浪效果远优于潮泥滩地、林带内波高呈指数衰减等规律 ，。物模研究方面，章家昌 根据试验结果得到植物消浪系数与来波波长、林带宽度、林木主干和枝叶遮蔽系数的定量表达式。虽然后续研究中纳入了波浪类型 、初始波高 、波周期 、波陡 ，、波流耦合作用 、植物淹没度 、植物垂直密度分布 、杆群刚度 、植物形态 等影响因素，但植物消浪系数的计算方法仍是以文献 为主。数模研究中，植物作用一般通过摩擦项体现 ，模型中往往以植物对波浪水流的阻力研究作为机理研究的基础，其中，求解拖曳力系数是量化植被在水中阻力的关键。在以往研究中 ，基于波浪通过林带后波高消减这一规律，根据波高消减模型率定得到阻尼因数，再通过阻尼因数与拖曳力系数及其他易于测量的参数（如滩地水深、平均植物密度）之间的联系，得到拖曳力系数的值，最后通过非线性回归得到拖曳力系数与雷诺数、邱卡数等参数之间的关系式以供数模计算。然而，拖曳力系数受不同试验工况（如植物密度）影响与重要参数呈不同的非线性关系 ，。本文在文献 基础上，基于“林带越宽、植物越密、水深越浅，植物消浪作用越强”的一般试验规律 ，结合 组物模试验结果，分析提出了新的植物消浪系数计算方法，通过 组试验数据开展了验证，并与规范方法结果进行了对比，揭示了影响植物消浪机理的关键无量纲参数，可根据地区消浪要求，结合植物直径、高度和水深、波要素条件，设计合理的植物种植宽度、密度，实现滩涂资源的有效利用。此外，得到的植物拖曳力系数计算方法可为相关数模研究提供基础。理论基础 计算植物消浪系数的规范方法常以波高的变化表征波浪能量的变化，以植物消浪系数量化滩地植物对波浪的消减作用。海堤工程设计规范（?）基于文献 提出防浪林消浪系数（）可按下列公式计算：?()?()?（）（）珚?（）珚?（）则林带消浪后的波高（）为：（）（）式中：为波长，；为林带宽度，；为入射波高，；为林木枝叶遮蔽系数；为林木主干遮蔽系数；为林木主干的平均直径，；为林木整体（包括主干和枝叶在内）的平均直径，；珚为林木成正方形排列的密度，株?；即 珚为林木成品字形交错排列时的密度。其中，式（）的适用范围为 ，。两种波高消减模型波高在种树带中的衰减规律主要可由两种波高消减模型描述。一种是 提出的林带沿程波高（）符合互反函数的消减模型：（）?（）?（）（?）（）式中：为无量纲缩尺相对波高，（）即为消浪系数；为阻尼因数，；（）为缩尺阻尼因数；为与林带前端的距离，；为与林带前端的缩尺距离。另一种是 等 根据深度淹没的人工巨藻（）试验并假设（）（），得到通过林带的沿程波高呈指数消减：（）?()()（）式中：为指数阻尼因数，；（）为缩尺指数阻尼因数，根据已有试验数据其取值范围为 ，）。等 结合两种模型，得到缩尺阻尼因数与缩尺指数阻尼因数的关系式：?（）（）当植物深淹没消浪效果有限时 ，值接近于 ，值大于但接近于 值，此时 。传统拖曳力系数求解方法 、等 、等 提出了阻尼因数或指数阻尼因数与拖曳力系数（）的关系式：?（）（）（）（）（）（）（）（）（）（）()（）式中：为植物直径，；为单位面积植物数量，株?；为滩地水深，；为波数，；为植物在水中的高度，若植物为非淹没工况其值等于滩地水深 ；为植物平均直径，即波浪传播方向上单位高度的面积。因此，拖曳力系数的率定方法可归纳如下：结合 提出的消减模型式（）及阻尼因数与拖曳力系数的关系式（）；结合式（）及 等 建立的关系式（）；结合 等 提出的消减模型式（）及指数阻尼因数与拖曳力系数的关系式（）；根据式（）率定指数阻尼因数，基于关系式（）得到阻尼因数，再结合经典关系式（）。根据文献 分析，方法和仅适用于植物淹没工况，方法和能在植物淹没和非淹没工况下表现良好。将率定得到的拖曳力系数与雷诺数、邱卡数等重要参数关联得到非线性关系式 ，并结合下式是数模中计算水流通过植物杆群的阻力（）的常用方法：（）式中：为流体密度，?；为阻力面积，；为流速，?。阻尼因数与拖曳力系数植物消浪系数越大代表植物消减波浪能量的作用越强，但式（）较复杂且尚未能直观显示影响植物消浪作用的关键因素。根据波高消减模型，若想消浪系数越大，则林带需要更大的（指数）阻尼因数。然而，阻尼因数历来是通过率定得到，现有成果对其仍了解甚少。而且，式（）和式（）显示，植物平均直径、植物密度、入射波高越小，水深、阻尼因数越大，则拖曳力系数越大。其中，“植物直径、植物密度越小，拖曳力系数越大”与以往对植物消浪效果研究的认识不相符，可能因为式中的阻尼因数本身又与植物直径、植物密度等参数有关，拖曳力系数的大小受多因子共同影响。因此，进一步研究影响阻尼因数的关键因素及定量表达式，有助于增加对植物消浪系数和拖曳力系数的认识。植物消浪作用特性分析图 波浪传播方向上植物面积占横截面比例示意图 阻尼因数法计算植物消浪系数为研究植物消浪作用，本 文 收 集 了 等 、等 、等 、等 等文献中共计 组在水平滩地上开展的消浪试验数据。各试验工况不同，例如模拟植物有的为淹没状态，有的为非淹没状态。首先通过收集的 组数据研究缩尺阻尼因数计算表达式。已有研究表明，种树带情况对波浪的影响巨大，且林带密度（）越大，植物消浪效果越好 ；植物淹没度（?）大于等于 时，消浪作用强于其小于 时 ，。结合文献 归纳的一般试验规律，分析得到可以表达林带情况的无量纲参数：波浪传播方向上植物面积占横截面比例?，其中，单位植物占比：?（为 方向距离，）；植物数量：。因此，把缩尺阻尼因数看作与具有物理意义的无量纲因子?有关，则该因数可用下式表示：（?）（）植物淹没状态是影响植物消浪作用的关键因素之一 ，。因此，区分淹没和非淹没工况，基于不同模拟植物的 组数据进行规划求解，得：?()（）式中对于植物非淹没工况 且?，对于淹没工况 。即缩尺阻尼因数在非淹没和淹没工况的表达式分别如式（）和式（）。不管林带内的 株植物如何排列，以 和 方向的株间距一致计，根据大量试验结果，单位 方向距离、单位高度内植物面积占横截面比例槡的范围为（，）。此外，对于特定工况（以 、等为参数），若想要植物淹没时与非淹没时有相同的缩尺阻尼因数 值，?需要成为非淹没时的 倍，说明在量化植物消浪作用时区分植物淹没情况是必要的：（）（非淹没工况）（）?()（淹没工况）（）阻尼因数法计算植物消浪系数 的表达式为：?（）（）所有组次植物消浪系数 法的计算值与实测值对比见图 。其中，实测值为根据实测沿程波高和波高消减模型所得的结果；计算值为每组工况根据实测参数由表达式计算得的结果；预测值为根据表达式的理论计算结果。图中不同符号表示不同学者的研究，实心符号表示植物非淹没工况，半填充符号表示植物淹没工况（下同）。结果显示，来自不同实验室的试验数据能较好地被综合分析。组工况得到的实测值与计算值结果接近，线性拟合结果斜率为 ，修订判定系数达 。计算值几乎在预测值 范围内，考虑试验误差，上述消浪系数计算方法是可行的。进一步地，根据式（）和（）可以快速预测相关参数变化对植物消浪效果的影响。例如，根据试验结果，非淹没林带宽 、密度 株?时植物消浪系数值为 ；根据式（），此时缩尺阻尼因数值为 。当密度 株?时，增加一倍，根据式（），为原来的 倍，即此时 值为 ，对应 值为 ；与试验结果“林带宽、密度 株?时的波高消减率约为 ”一致。指数阻尼因数法计算植物消浪系数类似的，对缩尺指数阻尼因数规划求解得：（）（非淹没工况）（）?()（淹没工况）（）指数阻尼因数法计算植物消浪系数 的表达式为：（）（）所有组次植物消浪系数 法的计算值与实测值对比见图 。结果显示，组工况得到的实测值与计算值线性关系良好，拟合结果修订判定系数为 ，且计算值基本在预测值 范围内。结果说明式（）和（）可以较好地描述缩尺指数阻尼因数，与式（）搭配可以有效计算植物消浪系数。同理，根据式（）至（）可以预测参数变化后植物消浪效果的变化情况。例如，根据试验结果，非淹没林带宽 、密度 株?时，根据式（），值为 。增加一倍时，根据式（），为原来的 倍，即此时 值为 ，对应 值为 ；与试验结果接近。此外，若考虑更多影响因素，在式（）至（）中增加其他无量纲数可能得到更好的拟合结果，但本文方法抓住了影响植物消浪作用的最主要因素且关系式相对简单，后续可加入例如植物排列方式、波浪要素等因素开展研究。图 消浪系数 法计算值与实测值对比图 消浪系数 法计算值与实测值对比 两种阻尼因数的关系不同研究中，根据实测参数由式（）得的缩尺阻尼因数计算值与由式（）（）得的缩尺指数阻尼因数计算值的关系见图 。结果显示，值小于约 时，与 计算值接近 。而且，对于淹没和非淹没工况，通过式（）可以很好地拟合 与 的关系，佐证了阻尼因数与指数阻尼因数的计算方法以及式（）的合理性。图 组工况 与 计算值对比 根据消浪系数计算林带后波高至此，通过入射波高和消浪系数计算表达式可以得到林带消浪后波高。根据 法和 法的计算值与实测值的对比分别见图和图 。结果显示，两种方法都可以良好地预测林带后波高；两者计算值与实测值线性关系的修订判定系数都达到约 ，斜率接近于 。考虑到试验误差及不同工况，结果令人满意。两种方法的计算值对比结果显示（图 略），线 性 关 系 斜 率 ，修 订 判 定 系 数 ，进一步说明两种波高消减模型都可以很好地描述波浪在林带拖曳作用下波高的消减情况。图 阻尼因数法计算值与实测值对比图 指数阻尼因数法计算值与实测值对比 根据阻尼因数计算拖曳力系数如果将阻尼因数由式（）带入式（）和（），可得到拖曳力系数的计算表达式：（）（）（）（）（）()（）（）（）当植物非淹没时，两式可简化为：（）（）（）（）（）()（）（）（）结果显示，拖曳力系数 与 正相关，与、等参数负相关。然而，根据试验结果，植物消浪作用与上述参数并无明显直接关系，可能因为拖曳力系数本身是一个数值解，物理意义并不明确，植物消浪作用的物理意义主要由阻尼因数及消浪系数体现。拖曳力系数根据式（）（）计算值与根据经典 等 公式率定值（第 节中方法）的对比关系分别见图 、图 。结果显示，计算值与率定值接近，绝大多数情况下，计算值在预测值 倍范围内。考虑到试验误差、多种工况以及拖曳力系数的性质，计算结果基本可行，上述拖曳力系数计算方法可酌情应用于数模研究。物模试验概况 本次试验设计本次物理模型试验在广东省水利水电科学研究院河口水利技术国家地方联合工程实验室的波浪水槽中进行，水槽长 、宽 、高 。图 拖曳力系数式（）计算值与率定值对比图 拖曳力系数式（）计算值与率定值对比本次水力模型水平滩地工况的总体设计及模型布局、仪器设备、模型树与文献 一致，但模拟植物材料、种树带宽度、波浪参数等设计不同。此外，基于自然规律下滩地坡度一般为 ，在水平滩地模型的基础上设计了坡度为 缓坡滩地模型试验（见图 ），以研究根据水平滩地试验数据所得的植物消浪计算方法在缓坡滩地上的适用性。试验模拟植物除了模型树 还采用了直径 材料的刚性柱。根据文献 ，中国大陆地区最大面积的连片红树林生长区内引种的无瓣海桑林平均高约 ；本地的白骨壤、木榄、红海榄、秋茄、桐花林高度分别约 、；保护区内小树密度 株?、大树密度 株?。因此，本试验设计原型高 和 的模型树，即模型树、刚性柱高度采用 、，植物密度采用 株?、株?。模型中，种树带宽度采用 、，滩地水深采用 、。缓坡滩地前缘水深为 时，植物为非淹没；前缘水深为 时，种树带计为淹没状态。缓坡滩地时计算滩地水深（）考虑两种：植物区水深的平均值；缓坡滩地前缘水深。经试验验证，两种考虑的结果几乎一致，为简单化处理，计算滩地水深采用缓坡滩地前缘水深。表 中罗列了不同影响因素组合下的 种工况，对应都测量了无植物时的空白试验，为 处初始波高；处的波高为种树带入射波高。历史数据归集为更充分研究植物消浪作用特性，本文还收集了谭超等 、陈晖 文献中不同工况共 组试验数据。谭超等 在长 、宽 、高 的波浪水槽中，采用 种不同刚度但外形相同的圆柱形硅胶棒开展了植物杆群消浪试验；模拟植物直径 、高度 ，林带密度 株?、宽度 ；滩地水深 ，规则波波高 、波周期 。陈晖 试验水槽长 、宽 、高 ，植被区长 ，模拟植物 管直径 、高度 ；滩地水深采用 、，植物密度采用 株?、株?；入射波高采用 、，波周期采用 、；共开展 组规则波实验。组试验资料的纳入有利于验证本文所提方法的适用性。图 缓坡滩地工况试验设置（滩地前缘水深 时，单位：）表 刚性柱和模型树试验工况工况?（株?）?模拟植物?、刚性柱?、刚性柱?刚性柱?、刚性柱?、模型树?、刚性柱?、模型树?、刚性柱?刚性柱?、刚性柱?、模型树?、刚性柱?、模型树?刚性柱?、模型树?、刚性柱?刚性柱?、刚性柱?刚性柱?、刚性柱?刚性柱?、刚性柱?刚性柱?、刚性柱?、刚性柱?、刚性柱?、刚性柱?、刚性柱?、刚性柱?、刚性柱；、模型树?、刚性柱?、刚性柱?、模型树?、模型树?、模型树?、模型树?、模型树注：工况 采用不规则波，其他工况采用规则波。工况 为缓坡滩地试验，其他工况为水平滩地试验。结果与讨论 植物消浪系数计算方法验证用本次试验的 组和收集的 组数据验证消浪系数计算方法。不同工况消浪系数根据阻尼因数法和指数阻尼因数法的计算值与实测值对比分别见图 和图 。结果显示，消浪系数实测值较小的工况，其计算值可能略高于实测值；消浪效果较好时，计算值接近于实测值。无论模型树或刚性柱，规则波或不规则波，植物淹没或非淹没，各工况消浪系数计算值基本在预测值 范围内。异常值为波周期 工况，可能跟试验测量波高精准度有关 。值得注意的是，消浪系数计算方法尚未考虑不同刚度杆群的影响。然而，根据文献 ，植物抗弯弹性模量值不小于 时可达到预期消浪效果。根据木材弹性模量资料，所有常见树木都满足该刚度，即计算方法适用于常见防浪树木。图 验证工况消浪系数 法计算值与实测值对比图 验证工况消浪系数 法计算值与实测值对比 林带后波高计算方法验证用 组试验数据验证林带消浪后波高计算方法。不同工况 根据法和 法的计算值与实测值对比分别见图 和图 。结果显示，所有工况两种方法计算值与实测值线性关系的修订判定系数分别高达 、；计算值接近于实测值，说明 计算方法效果好，能有效计算植物消浪作用。波高较大时，计算值可能偏小，即计算方法略高估了消浪效果较差时（例如刚性柱 较小）滩地植物的消浪作用。此外，基于水平滩地工况的计算方法略高估缓坡滩地工况下的试验效果，可能由于波浪波高较大时在水深较浅时易破碎，此时设计波长等参数与实际操作不符导致计算结果出现偏差。实际中可通过增加安全系数的方法进行修正。图 验证工况 法计算值与实测值对比图 验证工况 法计算值与实测值对比 与规范计算方法对比针对模型树试验，包括文献 的 组工况，水平和缓坡滩地工况下，规则波和不规则波作用下，根据规范方法的计算值与实测值对比见图 ，根据规范方法和本文方法（以阻尼因数法为例）所得的计算值对比见图 。结果显示，本文方法计算得到的林带后波高与规范方法接近，尤其对于水平滩地植物非淹没工况，进一步说明了本文方法的可靠性。规范方法仅考虑种树带宽度、波长、植物遮蔽系数、植物密度，未考虑植物高度、水深等垂向要素，在植物淹没时可能导致预测的消浪系数偏高从而预测的林带后波高偏低（图 ）。本文方法考虑了植物淹没情况，有利于准确描述更多工况下的植物消浪作用。此外，两种方法理论上都未考虑滩地形态的影响，由计算值低于实测值可推测缓坡滩地时植物的消浪作用比同等条件下水平滩地偏弱，实际中可通过增加安全系数的方法稍作修正。图 规范方法计算值与实测值对比图 规范方法与本文方法计算值对比 结论通过海岸防浪林消减波浪能量是一种基于自然的解决方案。为了解波高消减规律，文献 提供了两种可靠的消减模型：互反函数和指数函数。本文通过 组试验资料，基于波高消减模型和一般规律分析提出了影响植物消浪作用的关键无量纲参数，可用于求解植物消浪系数、林带后波浪波高、植物拖曳力系数。相关植物消浪作用计算方法又由本次试验的 组和收集的 组数据进行了验证，最后对比了由本文方法与规范方法的计算结果。主要结论如下：（）波高消减模型中的缩尺（指数）阻尼因数与波浪传播方向上植物面积占横截面比例?这个无量纲数密切正相关，揭示了植物消浪机理，可快速研判相关参数变化对植物消浪效果的影响。植物非淹没时，缩尺阻尼因数 （），缩尺指数阻尼因数 （）。和 基本符合?（），这增加了对消减模型的约束，有助于更准确研判植物消浪作用。相应的植物消浪系数、林带后波高、拖曳力系数结果也说明了 和 定量表达式的有效性。而且，相较于拖曳力系数，阻尼因数及消浪系数更能直观体现植物对波浪的消减效果。（）组试验数据良好地验证了消浪系数和林带后波高的计算方法。波高较大时，计算方法可能略高估消浪效果较差时滩地植物的消浪作用；滩地缓坡形态也可能对结果产生一定影响，实际中可考虑相关系数保证水安全。（）与规范方法相比，本文方法考虑了植物淹没程度的影响，拓展了消浪机理计算方法的适用范围，直观揭示了种树带密度、宽度、高度、植物平均直径、水深等对植物消浪作用的影响，可根据消浪要求合理设计植物种植宽度、密度以实现滩涂资源的有效利用，为指导植物防护工程设计提供依据。参考文献：黄本胜，赖冠文海堤外滩地种树效果及对行洪影响 人民珠江，（）：？，：，：，（）：，：，?，章家昌防波林的消波性能 水利学报，（）：吉红香植物消浪护岸试验研究 乌鲁木齐：新疆农业大学，刘达，黄本胜，邱静，等海岸破碎带植物消浪数值模拟研究 海洋工程，（）：，：，：龚尚鹏，陈杰，蒋昌波，等孤立波作用下植物带消浪特性对比分析 中国海洋大学学报（自然科学版），（）：，：，：，：谭超，黄本胜，刘达，等不同刚度植物杆群对规则波传播及紊动特性影响研究 海洋工程，（）：，：高学平，吕建璋，孙博闻，等含植物河道等效床面阻力试验研究 水利学报，（）：，陈杰，何飞，蒋昌波，等规则波作用下刚性植物拖曳力系数实验研究 水利学报，（）：，：，：，?，?：，（）：，：，（）：，陈清华中国红树林生物多样性调查（广东卷）青岛：中国海洋大学出版社，陈晖基于物理模型实验的波流耦合作用下植物消浪机理研究 广州：中山大学，（下转第 页），（，；，）：，：，；，：；（责任编辑：李娜）（上接第 页），（，；，；（），）：，；，：；（责任编辑：李福田）