含仿生草明渠水流流速结构试验研究

付长静; 王锦国; 赵天龙; 吕毅; Chang-jing FU; Jin-guo WANG; Tian-long ZHAO; Yi LÜ

中国农村水利水电 ›› 2021 ›› (10) : 55-65.

水环境与水生态

含仿生草明渠水流流速结构试验研究

作者信息 +

Experimental Research on the Flow Structure in the Open Channel with Bionic Grass

Author information +

稿件信息 +

摘要

仿生草防护技术是通过降低水下管道附近水流的流速，促进泥沙在管道附近进行沉降淤积以控制冲刷发展的防护技术。目前，对含有仿生草明渠水流流速的研究主要通过测定草前、中、后的流速变化规律，分析仿生草防护段泥砂淤积变化，极少对含仿生草明渠流的全流速场进行有效的测量。利用PIV测试技术对水流作用下仿生草防护段后的流速场进行了测量，分析了仿生草防护段后明渠水流沿流向垂直平面内时均流速分布规律以及草高及铺设间距对流速分布的影响，并深入探讨了仿生草防护段后明渠流脉动流速概率密度分布及脉动流速空间相关性。

Abstract

The problem of suspension treatment of oil and gas pipelines has been highly concerned by engineering construction units and researchers. Research shows that the bionic grass can effectively reduce the flow rate， promote sediment deposition， and control the development of the pipeline suspension area. The velocity distribution of open channel flow with bionic grass is very complex.The height and spacing of bionic grass will affect the river flow velocity distribution. At present， the flow velocity in the open channels containing bionic grass is mainly studied by measuring the velocity variation in front， middle and back of bionic grass， but few effective measurements are made for the full velocity field. This paper describes the use of modern means of advanced test by using standard particle image velocimetry （PIV） measurements with bionic grass along the water channel to the vertical plane of the distribution of velocity field. The probability density distribution and spatial correlation of pulsating velocity in the open channel after bionic grass section is further discussed.

关键词

仿生草 / 时均流速 / 脉动流速 / 空间相关性

Key words

bionic grass / time-averaged velocity / fluctuating velocity / spatial correlation

基金

重庆市科委基础研究与前沿探索项目(cstc2018jcyjAX0084)

水利部土石坝破坏机理与防控技术重点实验室开放基金项目(YK319006)

港口航道泥沙工程交通行业重点实验室开放基金项目(Yk220001-6)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

付长静 , 王锦国 , 赵天龙 , 吕毅. 含仿生草明渠水流流速结构试验研究[J].中国农村水利水电, 2021(10): 55-65

Chang-jing FU , Jin-guo WANG , Tian-long ZHAO , Yi LÜ. Experimental Research on the Flow Structure in the Open Channel with Bionic Grass[J].China Rural Water and Hydropower, 2021(10): 55-65

0　引言

仿生草是基于海洋仿生学原理，采用耐海水浸泡、抗长期冲刷的新型高分子材料加工的^［1］。国内外研究表明，仿生草不仅可以有效降低水流流速，促进泥沙淤积，而且不会产生次生环境污染^［2，3］，主要用于水下悬空管道的防护。目前关于仿生草的物理模型试验研究成果较少，主要讨论仿生草前中后特定位置的流速变化规律以及仿生草防护段泥沙淤积变化情况^［4］，极少对含仿生草防护段后明渠流的全场流速进行有效的测量研究，而含仿生草流道水流流速分布是进一步研究水流中泥砂以及物质的输运规律的基础。

Yang和Choi^［5］将明渠淹没植被化为两层，认为植被层内平均流速均匀分布、植被上层符合对数律分布。Nikora^［6］认为淹没植被水流流速主要由四部分线性叠加构成，①近壁区域植被底层的流速为均匀分布；②植被顶部附近混合层流速为双曲正切分布；③自由水层边界层为对数分布；④植被影响下的尾流函数。区域的划分会根据植被类型、试验布置、水流流态等不同而发生变化。槐文信等^［7，8］开展了柔性淹没植被明渠水流试验研究，发现植被上层区域纵向流速呈对数分布，植被区域纵向流速呈S型分布。吴福生等^［9，10］采用PVC薄片模拟柔性植被，用PIV测量不同位置的流速场，分析发现植被冠层以下区域，流速明显减小，随着流量的增加，植被对纵向流速的影响范围减少，而植物冠层处的流速梯度随流量和密度的增加而增大，且越靠近植物，植物冠层处的速度梯度越大。王忖等^［11］考虑了植物的生态作用，通过水槽试验，对沉水植物（水蕨）和挺水植物（菖蒲）段紊流特性进行测量研究。结果表明，含沉水植物和挺水植物明渠水流的垂向流速分布各自遵循不同的规律。

本文开展水槽试验，用粒子图像流速仪（PIV）对仿生草防护段后流速结构进行了测量研究，分析了不同工况下仿生草防护段后明渠水流的时均流速分布规律、脉动流速概率密度分布规律以及脉动流速空间结构，定量地阐述了仿生草的存在对水流流速结构的影响规律。

1 试验概况

本文物理模型试验在高精度变坡水槽中开展，该水槽长12.6 m，宽0.25 m，高0.25 m。为减小水槽边壁接头对水流结构的影响，水槽侧面和底面均由3.6 m长玻璃板组成，玻璃安装误差小于±0.2 mm，水槽结构变形小于±0.3 mm。为方便PIV进光以及照相机拍照，玻璃底板与边壁玻璃之间采用斜口黏接。水槽入口放置整流格栅，出口设置合页式尾门。试验中，在距仿生草垫铺设结束断面约0.1 m处测定水流瞬时流速，测定距离为0.08 m，如图1所示。

图1 水槽试验示意图

Fig.1 Schematic diagram of flume test

Full size|PPT slide

试验水流为均匀流，流量采用电磁流量计测量，采样频率为10 Hz，测量精度为0.4%；水深采用超声水位计测量。流场采集窗口中部距离水槽入水口约8 m，满足紊流充分发展所需进口段长度要求，距离水槽出口段约4 m，可充分消除尾门对水流的扰动。试验采用二维高频PIV系统沿水槽纵向中轴面进行流场采集，该系统主要由高频CMOS相机、8W半导体连续激光和PIV流场计算软件组成。CMOS相机最高像素为256×1 920，满幅最高帧频为800 Hz。激光束经棱镜转变为45°角的片光，厚度约为1 mm，片光源可从水槽玻璃底板进光，可穿透床面。试验采用的示踪粒子是直径为8~12 μm、密度为1.03 g/mm³的空心玻璃球。PIV流场计算采用WIDIM多重网格迭代图像变形算法，最小诊断窗口为16×16像素。试验以低频独立采样方式采集流场，两次瞬时流场平均时间间隔为1 s，同一个瞬时流场对应2帧图像的时间间隔为1.25 ms。各工况样本容量均为 5 000 次（10 000帧流场图像），样本容量满足研究要求^［12］。试验采用正态模型，模型中的流速根据重力相似准则确定，试验比例尺为1/20。仿生草垫试样中采用的仿生草采用T25型工程样品，仿生草垫实际尺寸为5 m×5 m，草垫中仿生草各排布置间距为1.67 m，草高为1.5 m。试验中仿生草的铺设均以实际工程存在的条件为参照，同时为了分析草高和铺设间距对水流流速结构的影响，本试验设置了不同草高和不同铺设间距的样本，并且每一个工况又分成3组平行样，试验工况参数如表1所示，仿生草布置如图2所示。

表1 试验工况参数

Tab.1 Test conditions and related parameters

序号	工况	水深d/m	流量Q/（L·s^-1）	流速v/（m·s^-1）	Re	Fr	谢才系数 C	运动黏滞系数υ/（10^-6m²·s^-1）	草高H_g /m
1	B25S	0.109	6	0.18	8 324.7	0.171	624.065	1.27	0.075
2	B25M	0.112	5	0.22	10 268.4	0.217	622.568
3	B25L （H25-1.5、D252）	0.090	6	0.27	11 040.4	0.285	611.579
4	H25-1	0.090	6	0.27	11 040.4	0.285	611.579	1.27	0.050
5	H25-1.2	0.090	6	0.27	11 040.4	0.285	611.579	1.27	0.060
6	D251	0.090	6	0.27	11 040.4	0.285	611.579	1.27	0.075
7	D253	0.090	6	0.27	11 040.4	0.285	611.579	1.27	0.075

图2 仿生草布置图

Fig.2 Schematic diagram of bionic grass layout

Full size|PPT slide

2 试验结果分析

2.1 仿生草防护段后纵向时均流速分布

仿生草防护段后水流纵向时均分布，如图4所示。仿生草的存在增加了水流的阻力，水流流经仿生草后改变了原有流速分布：①渠底附近水流流速明显减小，部分区域出现负值，产生漩涡，并且随着断面平均流速的增加，渠底附近流速全部为负值，出现大范围回流；②仿生草倒伏位置（见表2）偏下（约2 cm）附近靠近仿生草的水流纵向流速明显增加，说明当水流流过仿生草后，过流断面突然增大，导致水面急剧降落，发生水跌；③仿生草倒伏位置偏下（约2 cm）附近远离仿生草的位置水流纵向流速逐渐减小；④仿生草倒伏位置以上水流流速较小，且多为负值，这是由于水流流经仿生草时，倒伏后的草面成为新的床面，使水流发生局部变形而分离，产生旋转运动的副流随着水流向后移动。

图3 不同工况条件下仿生草倒伏

Fig.3 Lodging values of bionic grass at different working condition

Full size|PPT slide

图4 不同工况仿生草防护段后纵向时均流速分布云图

Fig.4 The distribution of the time-averaged longitudinal velocity in different working conditions

Full size|PPT slide

表2 不同工况下仿生草倒伏情况

Tab.2 The lodging of bionic grass under different working conditions

序号	工况	原始草高/cm	倒伏后草高/cm	序号	工况	原始草高/cm	倒伏后草高/cm
1	B25S	7.5	7.0	5	H25-1.2	6.0	5.1
2	B25M	7.5	6.5	6	H25-1.5 （D252）	7.5	5.5
3	B25L	7.5	5.5	7	D251	7.5	5.5
4	H25-1	5.0	4.6	8	D253	7.5	5.4

为分析不同工况下水流纵向流速垂线分布，取流速变化最大断面x=2 cm。由图5可见，当断面平均流速较小时，仿生草段后渠底2 cm范围内水流纵向流速呈线性，且大小接近于0值，在渠底2 cm以上，纵向流速垂线分布呈S型；随着断面平均流速的增大，仿生草段后渠底2 cm范围内水流纵向流速则沿水深呈抛物线分布，在渠底2 cm以上，纵向流速垂线分布近似S型，垂线上纵向时均流速的最大值分别在y=5.4 cm、y=4.9 cm、y=3 cm处取得，最大值分别为0.235、0.30、0.24 m/s，可见垂线上纵向时均流速的最大值位置随着断面平均流速增大而逐渐下移，这与仿生草的倒伏有关。同时，在相同水流条件下，不同草高和铺设间距的仿生草段后纵向流速垂线分布各不相同，尤其是铺设间距对仿生草段后水流纵向流速影响较大。

图5 仿生草防护段段后纵向时均流速垂线分布（x=2 cm）

Fig.5 The vertical distribution of the time-averaged longitudinal velocity in different working conditions

Full size|PPT slide

2.2 仿生草段后脉动流速概率密度分布

明渠紊流的脉动流速随时空变化特征可采用概率密度分布表示，为了对紊流运动的基本特征进行预测，本文对脉动流速的概率密度P进行研究，计算出各个工况脉动流速出现的概率，绘制了纵向和垂向脉动流速概率分布曲线，图中“-”号表示方向。图6为不同工况纵向脉动流速概率密度分布。可见，仿生草防护段后明渠水流纵向脉动流速概率密度均呈正态分布；随着断面平均流速增大，纵向脉动流速概率密度减小，正态分布峰度系数减小，峰形走势越平坦；在相同水流条件下，随着仿生草高度的增加，纵向脉动流速概率密度增大，峰形走势越陡峭，而随着仿生草铺设间距的减小，纵向脉动流速概率密度先增大后减小，峰形走势趋于平坦。

图6 不同工况纵向脉动流速概率分布

Fig.6 The probability distribution of longitudinal pulsation velocity in different working conditions

Full size|PPT slide

不同工况垂向脉动流速概率分布如图7所示。可见，仿生草防护段后明渠水流垂向脉动流速概率密度呈正态分布，但与纵向脉动流速概率密度分布相比，整体峰形走势分布更为陡峭，正态分布的峰度更大，峰顶更为尖锐；当断面平均流速较小时，垂向脉动流速概率密度达到62.5%，这说明当来流速度较低时，仿生草段后垂向脉动流速波动较小，而随着断面平均流速的增大，垂向脉动流速概率密度降低，峰度系数减小，整体峰形走势趋于平坦；在相同水流条件下，随着草高的增加，垂向脉动流速概率面密度分布的峰度逐渐增加，随着仿生草铺设间距的减小，垂向脉动概率密度虽然先增大后减小，但垂向脉动流速概率分布曲线的对称性越好，整体峰形势走势更平坦。

图7 不同工况垂向脉动流速概率分布

Fig.7 The probability distribution of vertical pulsation velocity in different working conditions

Full size|PPT slide

2.3 仿生草段后脉动流速空间相关性

通常脉动流速的空间相关结构可采用相关系数定量描述，相关系数越大，流速空间相关性越好，反之则越差。脉动流速的空间相关性包括自相关性和互相关性2类，

C_{u u}

、

C_{v v}

为自相关系数，

C_{u u} = \frac{\frac{1}{n} \sum_{i = 1}^{n} u_{i} (x_{0}, y_{0}) u_{i} (x_{1}, y_{1})}{u^{'} (x_{0}, y_{0}) u^{'} (x_{1}, y_{1})}

，

C_{v v} = \frac{\frac{1}{n} \sum_{i = 1}^{n} v_{i} (x_{0}, y_{0}) v_{i} (x_{1}, y_{1})}{v^{'} (x_{0}, y_{0}) v^{'} (x_{1}, y_{1})}

；C_uv 、C_vu 为互相关系数，

C_{u v} = \frac{\frac{1}{n} \sum_{i = 1}^{n} u_{i} (x_{0}, y_{0}) v_{i} (x_{1}, y_{1})}{u^{'} (x_{0}, y_{0}) v^{'} (x_{1}, y_{1})}

，

C_{v u} = \frac{\frac{1}{n} \sum_{i = 1}^{n} v_{i} (x_{0}, y_{0}) u_{i} (x_{1}, y_{1})}{v^{'} (x_{0}, y_{0}) u^{'} (x_{1}, y_{1})}

；（x ₀，y ₀）和（x ₁，y ₁）为空间任意两点的坐标；n为样本容量、u、v 分别为纵向、垂向脉动流速；

u^{'} = \sqrt[]{\frac{1}{n} \sum_{i = 1}^{n} {u^{2}}_{i}}

；

v^{'} = \sqrt[]{\frac{1}{n} \sum_{i = 1}^{n} {v^{2}}_{i}}

分别为纵向、垂向紊动强度。

本文选取流场采集窗口的几何中心点为（x ₀，y ₀），其周围空间点为（x ₁，y ₁），计算自相关系数和互相关系数，分析含仿生草段后明渠水流脉动流速空间相关结构特征及变化规律。

2.3.1 空间自相关性

图8和图9分别为仿生草防护段后明渠流脉动流速的自相关系数C_uu 、C_vv 云图分布。结果表明，中心点周围C_uu 和C_vv 呈正相关的空间结构；C_uu 、C_vv 值随分析点与采集窗口的几何中心点距离的增大而减小，说明两点距离越大，脉动流速空间自相关越弱；C_uu 和C_vv 云图最内圈大致呈圆形，且圆心附近C_uu 和C_vv _{基本为各向同性；随着平均断面流速的增大，C_uu 和C_vv 云图最内圈更接近圆形，且最内圈的面积先增大后减小；随着仿生草高度的增加，C_uu 和C_vv 云图最内圈的面积同样先增大后减小。}

**图8 不同工况C_uu 云图**

Fig.8 The distribution of C_uu in different working conditions

Full size|PPT slide

**图9 不同工况C_vv 云图**

Fig.9 The distribution of C_vv in different working conditions

Full size|PPT slide

为进一步定量分析C_uu 、C_vv 云图分布特征，定义面积系数ζ=A_r /A，其中A_r 为相关系数r对应等值线所包围的面子，A为分析窗口总面积。显然ζ越大，表明该相关系数等值线所包围的面积占窗口总面积的比例越大。分别采用ζ_uu 、ζ_vv 、ζ_uv 和ζ_vu 表示C_uu 、C_vv 、C_uv 和C_vu 的面积系数。表3和表4为不同工况下ζ_uu 和ζ_vv 值。

**表3 不同工况下的ζ_uu 值**

Tab.3 The value of ζ_uu in different working conditions

相关系数r	B25S	B25M	B25L	H25-1	H25-1.2	H25-1.5	D251	D252	D253
<-0.14	0	0	0	0	0	0	0	0	0
-0.14 ~ -0.10	0	0	0.001	0	0.001	0.001	0.002	0.001	0.012
-0.10 ~ -0.06	0.011	0.008	0.037	0.012	0.027	0.037	0.031	0.037	0.055
-0.06 ~ -0.02	0.080	0.057	0.203	0.145	0.139	0.203	0.127	0.203	0.077
-0.02 ~ 0.02	0.214	0.158	0.352	0.325	0.326	0.352	0.225	0.352	0.170
0.02 ~ 0.06	0.328	0.278	0.255	0.317	0.306	0.255	0.236	0.255	0.303
0.06 ~ 0.1	0.210	0.236	0.113	0.115	0.115	0.113	0.217	0.113	0.195
0.10 ~ 0.14	0.081	0.119	0.029	0.051	0.046	0.029	0.119	0.029	0.104
0.14 ~ 0.18	0.047	0.067	0.004	0.018	0.016	0.004	0.031	0.004	0.053
0.18 ~ 0.22	0.021	0.041	0.002	0.007	0.011	0.002	0.005	0.002	0.017
0.22 ~ 0.26	0.005	0.022	0.001	0.003	0.004	0.001	0.002	0.001	0.004
0.26 ~ 0.30	0.001	0.009	0.001	0.002	0.002	0.001	0.001	0.001	0.003
>0.30	0.001	0.004	0.003	0.006	0.007	0.003	0.003	0.003	0.007

**表4 不同工况下的ζ_vv 值**

Tab.4 The value of ζ_vv in different working conditions

相关系数r	B25S	B25M	B25L	H25-1	H25-1.2	H25-1.5	D251	D252	D253
<-0.14	0	0	0	0	0	0	0	0	0
-0.14～-0.10	0	0	0	0	0	0	0	0	0
-0.10～-0.06	0.011	0.006	0.008	0.006	0.004	0.008	0.018	0.008	0.008
-0.06～-0.02	0.136	0.121	0.174	0.115	0.137	0.174	0.226	0.174	0.184
-0.02～0.02	0.405	0.397	0.506	0.397	0.502	0.506	0.499	0.506	0.490
0.02～0.06	0.329	0.313	0.271	0.316	0.299	0.271	0.222	0.271	0.242
0.06～0.10	0.096	0.117	0.037	0.097	0.046	0.037	0.029	0.037	0.056
0.10～0.14	0.017	0.034	0.003	0.037	0.007	0.003	0.003	0.003	0.012
0.14～0.18	0.003	0.009	0.001	0.019	0.003	0.001	0.001	0.001	0.004
0.18～0.22	0.001	0.002	0	0.006	0.001	0	0	0	0.001
0.22～0.26	0	0	0	0.003	0	0	0	0	0.001
0.26～0.30	0	0	0	0.002	0	0	0	0	0
>0.30	0	0.001	0	0.002	0.001	0	0	0	0.001

由表3和表4可知，ζ_uu 和ζ_vv 随r值的增大一般先增大后降低，r越大（大于0.02时），递减变化越明显，说明r较大时等值线所包围面积的变化更敏感；r约在0.02时，ζ_uu 和ζ_vv 值达到最大值；随着平均断面流速的增加，ζ_uu 和ζ_vv 的最大值均先减小后增大；在相同水流条件下，随着仿生草高度的增加，ζ_uu 和ζ_vv 的最大值逐渐递增，随着仿生草铺设间距的增大，ζ_uu 和ζ_vv 的最大值先增大后减小；

2.3.2 空间互相关性

图10和图11分别为仿生草防护段后明渠流脉动流速的互相关系数C_uv 、C_vu 云图分布。结果表明，C _uv 、C _vu 值随分析点与采集窗口的几何中心点距离的增大而减小，说明两点距离越大，脉动流速空间互相关越弱；C _uv 和C _vu 云图呈点状分布，中心附近C _uv 和C _vu 存在各向异性；随着平均断面流速的增大，C _uv 、C _vu 值先增大后减小；在同样的水流条件下，随着仿生草高度的增加，C _uv 、C _vu 值逐渐减小，随着仿生草铺设间距的增加，C _uv 、C _vu 值逐渐增大。

**图10 不同工况C_uv 云图**

Fig.10 The distribution of C_uv in different working conditions

Full size|PPT slide

**图11 不同工况C_vu 云图**

Fig.11 The distribution of C_vu in different working conditions

Full size|PPT slide

表5和表6为不同工况下ζ_uv 和ζ_vu 值。可见，ζ_uv 和ζ_vu 随r值的增大一般先增大后降低，递增和递减趋势基本一致，且r=0时，ζ_uv 和ζ_vu 值达到最大值；随着平均断面流速的增加，ζ_uv 的最大值逐渐增大，ζ_vu 的最大值先减小后增大；在相同水流条件下，随着仿生草高度的增加，ζ_uv 的最大值先增大后减小，ζ_vu 的最大值逐渐减小，随着仿生草铺设间距的增大，ζ_uv 的最大值先减小后增大，ζ_vu 的最大值先增大后减小。

**表5 不同工况下的ζ_uv 值**

Tab.5 The value of ζ_uv in different working conditions

相关系数r	B25S	B25M	B25L	H25-1	H25-1.2	H25-1.5	D251	D252	D253
<-0.10	0	0	0	0.001	0	0	0.001	0	0
-0.10～-0.08	0.003	0.001	0.002	0.007	0.002	0.002	0.007	0.002	0.001
-0.08～-0.06	0.015	0.010	0.017	0.028	0.015	0.017	0.033	0.017	0.009
-0.06～-0.04	0.062	0.045	0.074	0.066	0.055	0.074	0.096	0.074	0.045
-0.04～-0.02	0.145	0.114	0.162	0.129	0.144	0.162	0.200	0.162	0.145
-0.02～0	0.214	0.188	0.246	0.201	0.241	0.246	0.267	0.246	0.269
0～0.02	0.222	0.225	0.243	0.225	0.259	0.243	0.224	0.243	0.259
0.02～0.04	0.168	0.205	0.164	0.176	0.178	0.164	0.121	0.164	0.164
0.04～0.06	0.103	0.133	0.073	0.104	0.084	0.073	0.041	0.073	0.076
0.06～0.08	0.048	0.056	0.017	0.046	0.020	0.017	0.009	0.017	0.026
0.08～0.10	0.016	0.016	0.002	0.014	0.002	0.002	0.001	0.002	0.006
0.10～0.12	0.004	0.005	0	0.002	0	0	0	0	0
>0.12	0.001	0.001	0	0	0	0	0	0	0

**表6 不同工况下的ζ_vu 值**

Tab.6 The value of ζ_vu in different working conditions

相关系数r	B25S	B25M	B25L	H25-1	H25-1.2	H25-1.5	D251	D252	D253
<-0.10	0.001	0	0.001	0	0	0.001	0.001	0.001	0
-0.01～-0.08	0.006	0.003	0.004	0.001	0.002	0.004	0.005	0.004	0.001
-0.08～-0.06	0.020	0.014	0.026	0.006	0.017	0.026	0.026	0.026	0.008
-0.06～-0.04	0.054	0.049	0.089	0.034	0.072	0.089	0.091	0.089	0.043
-0.04～-0.02	0.118	0.111	0.201	0.118	0.186	0.201	0.197	0.201	0.135
-0.02～0	0.193	0.180	0.281	0.246	0.282	0.281	0.276	0.281	0.250
0～0.02	0.230	0.215	0.240	0.306	0.255	0.240	0.237	0.240	0.279
0.02～0.04	0.191	0.193	0.118	0.204	0.133	0.118	0.121	0.118	0.186
0.04～0.06	0.113	0.132	0.034	0.070	0.044	0.034	0.038	0.034	0.078
0.06～0.08	0.050	0.067	0.006	0.013	0.009	0.006	0.007	0.006	0.017
0.08～0.10	0.018	0.026	0.001	0.002	0.001	0.001	0.001	0.001	0.002
0.10～0.12	0.005	0.008	0	0	0	0	0	0	0
>0.12	0.001	0.002	0	0	0	0	0	0	0

3 结论

本文通过开展室内仿生草物理模型试验，对仿生草防护段后水流纵向时均流速、脉动流速概率密度、流速空间相关性进行了分析，得出以下结论。

（1）纵向流速分布的实测结果表面，在仿生草的影响下，水流的纵向流速垂向分布受很大的影响，不再遵循典型的对数分布规律，存在明显的分区。当断面平均流速较小时，仿生草段后渠底附近水流纵向流速基本呈线性分布，在渠底2 cm以上，纵向流速垂线分布呈S型，随着断面平均流速的增大，仿生草段后渠底附近水流纵向流速沿水深呈抛物线分布，在渠底2 cm以上，纵向流速垂线分布近似S型，并且纵向流速垂向分布的最大值的位置与仿生草的倒伏高度有关。

（2）纵向脉动流速的概率密度分布柱状图呈较好的对称性，分布曲线近似呈正态分布。随着断面平均流速的增加，正态分布的峰态系数减小，整体峰形的走势越趋于平坦。随着仿生草高度的增加，纵向脉动流速概率密度增大，峰形走势越陡峭。随着仿生草铺设间距的减小，纵向脉动流速概率密度先增大后减小，峰形走势趋于平坦。垂向脉动流速概率分布整体峰形的走势比纵向脉动流速概率分布更为陡峭。随着仿生草高度的增加，垂向脉动频率正态分布的峰度逐渐增加；随着仿生草铺设间距的减小，垂向脉动概率分布最大值呈先增大后减小的趋势。

（3）C_uu 和C_vv 云图最内圈大致呈圆形，且圆心附近C_uu 和C_vv _{基本为各向同性，C_uu 与C_vv 的最大相关系数约为0.3，r与ζ_uu 和ζ_vv 呈非线性关系；相关系数约在0.02时，ζ_uu 和ζ_vv 达到最大值，r与0.02差值的绝对值越小，ζ_uu 和ζ_vv 越大；C _uv 和C _vu 云图呈点状分布，中心附近C _uv 和C _vu 存在各向异性，C_uv 与C_vu 的最大相关系数约为0.12，r与ζ_uv 和ζ_vu 呈非线性关系；相关系数约在0时ζ_uv 和ζ_vu 最大，r>0时，r与0差值的绝对值越小，ζ_uv 和ζ_vu 越大。}

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	李俊杰，侯志民，田培胜.巨浪冲蚀威胁下的海底管道仿生草防护技术［J］. 海岸工程， 2017，36（4）：37-43. 本文引用 [1]

2	庄丽华，阎军，范奉鑫，等.茂名30万吨级单点码头输油海管悬空段人工海草试验研究［J］.海洋工程，2010，28（2）：76-81. 本文引用 [1]

3	冯卫兵，汪涛，邓伟. 柔性植物消波特性试验研究［J］.科学技术与工程，2012，12（26）：6 687-6 690. 本文引用 [1]

4	孟强，吴云云.滩海工程中仿生草消浪特性试验研究［J］.水利与建筑工程学报，2014（4）：185-190. 本文引用 [1]

5	YANG W J， CHOI S U ．A two-layer approach for depth limited open channel flows with submerged vegetation ［J］.Journal of Hydraulic Research， 2010，48（4）：466-475. 本文引用 [1]

6	NIKORA N， NIKORA V， O’ DONOGHUE T. Velocity profiles in Vegetated open-channel flows： Combined effects of multiple mechanisms ［J］.Journal of Hydraulic Engineering， 2013，139（10）：1 021-1 032. 本文引用 [1]

7	槐文信，韩杰，曾玉红，等.淹没柔性植被明渠恒定水流水力特性的试验研究［J］.水利学报，2009，40（7）：791-797. 本文引用 [1]

8	槐文信，韩杰，曾玉红，等.基于掺长理论的淹没柔性植被水流流速分布研究［J］.应用数学和力学，2009，30（3）：325-332. 本文引用 [1]

9	吴福生，王文野，姜树海. 含植物河道水动力学研究进展［J］. 水科学进展，2007，18（3）：456-461. 本文引用 [1]

10	吴福生，姜树海，邢磊，等.用PIV测量含淹没刚性植物明渠水流流速场［J］.武汉大学学报（工学版），2009，42（5）：587-591. 本文引用 [1]

11	王忖，王超.含挺水植物和沉水植物水流紊动特性［J］.水科学进展，2010，21（6）：816-822. 本文引用 [1]

12	张建梅，钟亮，刘珺洁.密排粗糙床面明渠水流的流速结构［J］.水电能源科学，2018，36（10）：108-112. 本文引用 [1]

PDF(14092 KB)

623

访问

引用

详细情况

段落导航

摘要
Abstract
关键词
Key words
基金
引用本文
0　引言
1 试验概况
图1 水槽试验示意图
表1 试验工况参数
图2 仿生草布置图
2 试验结果分析
2.1 仿生草防护段后纵向时均流速分布
图3 不同工况条件下仿生草倒伏
图4 不同工况仿生草防护段后纵向时均流速分布云图
表2 不同工况下仿生草倒伏情况
图5 仿生草防护段段后纵向时均流速垂线分布（x=2 cm）
2.2 仿生草段后脉动流速概率密度分布
图6 不同工况纵向脉动流速概率分布
图7 不同工况垂向脉动流速概率分布
2.3 仿生草段后脉动流速空间相关性
2.3.1 空间自相关性
图8 不同工况Cuu 云图
图9 不同工况Cvv 云图
表3 不同工况下的ζuu 值
表4 不同工况下的ζvv 值
2.3.2 空间互相关性
图10 不同工况Cuv 云图
图11 不同工况Cvu 云图
表5 不同工况下的ζuv 值
表6 不同工况下的ζvu 值
3 结论
参考文献

收稿日期	出版日期
2021-01-13	2021-10-15
发布日期
2021-11-09

选择文件类型/文献管理软件名称

选择包含的内容