方形人工鱼礁单体流场效应的PIV 试验研究

2014-02-15付东伟陈勇陈衍顺栾曙光

大连海洋大学学报 2014年1期

付东伟，陈勇，陈衍顺，栾曙光

(1.大连海洋大学海洋与土木工程学院，辽宁大连116023;2.大连海洋大学辽宁省海洋牧场工程技术研究中心，辽宁大连116023;3.福建省水产设计院，福建福州350003)

人工鱼礁(以下简称鱼礁)具有良好的环境功能，不仅能够修复海洋生态环境和优化渔业结构，而且对于净化固体废弃物和发展休闲渔业也有良好作用。随着海洋牧场建设的深入展开，对鱼礁流场效应的研究取得了突破性进展。李珺等［1］采用大涡模拟(LES)紊流模式，对单体方形鱼礁的流场进行三维模拟试验，得到礁体不同结构的通透系数与礁体周围流场变化的数学模型;郑延璇等［2］利用计算流体力学CFD 软件，采用RNG k-ε湍流模型，研究不同来流流速(0.3、0.5、0.8 m/s)下星体形鱼礁模型的流场效应，分析了礁体投放方式与流场效应的关系;李晓磊等［3］利用CFD 软件，对定常流作用下方形鱼礁模型进行数值模拟，研究了礁体周围的三维涡结构。但有关人工鱼礁开口比和来流流速对流场效应影响的研究目前尚未见报道。本研究中，作者采用粒子图像测速(PIV)二维流场测速技术对单体鱼礁进行定量物理模型试验，再应用双因素方差分析法讨论鱼礁开口比和来流流速对流场效应的影响程度，以期为鱼礁开口比的选取和鱼礁的合理布局提供参考。

1 材料与方法

1.1 试验装置

物理模型试验在大连理工大学海岸及近海工程国家重点实验室进行。试验用循环水槽长、宽、高分别为22.0、0.45、0.60 m。试验通过ADV 流速仪在试验区前端控制所需流速，待流速稳定后，采用激光片光源照射模型礁体与水流方向垂直的中垂面流场，用CCD 相机连续记录图像场，得到一个时间段的12 位灰度图像，然后采用美国TSI 公司的PIV 测量系统，对鱼礁模型进行二维片流场的流动测量和流动显示，并由粒子成像系统对返回计算机的数据进行分析，得到片流场的流速矢量图。PIV 测量系统试验装置示意图见图1。试验水流动力源由变频器(0 ～50 Hz)提供，当水槽中水深为0.5 m 时可以提供0 ～250 mm/s 的水流流速。根据重力相似准则，试验中几何相似的线性长度比λ(原型的线性长度与模型的线性长度的比值)为20［4］。本研究中选取原型所在位置的水流流速为0.8 m/s，对应模型所在位置的水流流速为180 mm/s，且变频器可以满足水流流速的需求。示踪粒子根据试验环境要求选用聚氯乙烯(PVC)粉。

1.2 模型类型和试验工况

开口比是礁体迎流面孔洞在横向中垂面的投影面积与迎流面投影面积的比值。本次试验模拟实际鱼礁尺寸为边长2 m 的立方块体，开口比分别为0、0.04、0.16、0.25。试验模型材料为有机玻璃，根据重力相似准则，模型边长为100 mm 的立方块体，并且中心分别开有边长为0、20、40、50 mm的孔洞，4种单体鱼礁模型示意图如图2所示。

图1 PIV 测量系统试验装置示意图Fig.1 Experimental equipment in the PIV system

图2 开口比为0、0.04、0.16、0.25 的单体鱼礁模型Fig.2 Models of artificial reefs of cut－opening ratio of 0，0.04，0.16，and 0.25

2 结果与讨论

2.1 流场效应评价指标

鱼礁投置于海底后产生的上升流可以促进上下层海水的交换，扩大基础饵料和营养盐的分布水平，易于诱集各层鱼类形成渔场;背涡流可以形成流体力学阴影区域［5］，提供鱼类休憩和躲避强水流的场所。本研究中流场效应的评价指标选取平均上升流流速和背涡流面积，如果流场中平均上升流流速较大和背涡流面积较大，则判定鱼礁流场效应较显著。流场矢量数据提取于连续时间段内鱼礁纵向中轴流场矢量图中流场效应最明显的时刻(图3)。平均上升流流速为礁体迎流面一端的水流分离点前一定范围内垂直流速的平均值(9 组数据);背涡流面积为礁体后侧产生最大涡旋区域的面积。

2.2 流场效应分析

从图3可见:单体鱼礁流场中开口比为0 的鱼礁背涡流区域较大，涡流明显，并且随着来流流速的减小，涡心位置升高，开口比由0 增至0.25 时，礁体中心区域水流流量增大，背涡流区域涡流面积减小，且形成礁体后侧上下两个涡流区域;上升流流速随着来流流速的增加而增大，来流流速较大时礁体周围的流场矢量分布呈现出明显的变化规律，即来流流速越大上升流高度越大，对背涡流影响就越大，当来流流速为67 mm/s 时流场矢量分布规律不明显，这是因为实际物理模型试验中，流速越小，流场受到周围环境因子的影响就越大，如水槽两侧玻璃对水流的黏滞性和PVC 粉的不稳定性等。

2.2.1 上升流流速从图4可见:在相同来流流速工况下，单体鱼礁随着礁体开口比由0 增至0.25 时，平均上升流流速逐渐降低，平均上升流流速与来流流速的比值逐渐降低，分别降低了34.67%、21.35%、15.31%、10.7%，这是因为开口比越小，礁体迎流面受到空隙的分流就越少，礁体上侧的流管效应也越强;在开口比相同的工况下，随着来流流速的降低，平均上升流流速逐渐降低，平均上升流流速与来流流速的比值呈降低趋势，比值分别降低了 24.67%、21.40%、15.45%，即来流流速越大，流场中上升流流速越大。其中，开口比为0 的鱼礁流场中，平均上升流流速与来流流速比值为0.25 ～0.42。刘洪生等［6］的风洞试验结果表明，正方体模型在实际来流流速为0.4、0.6、0.8 m/s 的工况下，平均上升流流速为来流流速的0.3 倍。本试验结果与此结果基本吻合。

图4 平均上升流流速与来流流速比值随开口比的变化Fig.4 The curves of the ratio between average upwelling velocity and incoming flowing velocity with different cut－opening ratios in the flow field

2.2.2 背涡流面积从图5可见:在相同来流流速工况下，单体鱼礁的开口比由0 增至0.25 时，背涡流面积逐渐减小，平均背涡流面积分别为7734.7、4325.4、2578.5、1148.1 mm2，即随着开口比的减小，背涡流面积逐渐增大，这是因为礁体中心贯穿空隙越小，流经礁体中心至礁体后侧的水流流量就越低，也就越有利于形成背涡流，所以选择开口比较小的鱼礁有利于增大礁体后侧背涡流的面积;在开口比相同的工况下，随着来流流速的增大，背涡流面积逐渐增大，平均背涡流面积分别为4925.4、3850.1、3064.5 mm2。

图5 背涡流面积与礁体迎流面面积比值随开口比的变化Fig.5 The curves of ratio between area of reverse vortex and area of incident flow surface with different cut－opening ratios in the flow field

2.3 开口比和来流流速对流场效应影响的双因素方差分析

采用无重复双因素方差分析法［7］对试验数据进行分析，步骤如下:

本研究假设开口比为因素A，来流流速为因素B。通过双因素方差分析研究鱼礁流场开口比和来流流速对平均上升流流速与来流流速比值的影响［8］。从表1和表2可见:当显著性水平α 为0.05 时，因素A和B 均表现出显著性影响(P＜0.05)，表明在本次物理试验中，鱼礁开口比和来流流速对平均上升流流速与来流流速比值均有显著性影响;当显著性水平α 为0.01 时，因素A表现出极显著性影响(P＜0.01)，而因素B 则无极显著性影响(P＞0.01)。根据上述分析可知，礁体的开口比是影响鱼礁流场效应的主要因素，而流场中来流流速对鱼礁流场效应的影响程度次之。

表1 平均上升流流速与来流流速比值的双因素方差分析Tab.1 Variance analysis of ratio between average upwelling velocity and incoming flowing velocity

表2 不同因素的显著性分析Tab.2 Significance analysis between different factors

通过双因素方差分析研究鱼礁流场开口比和来流流速对背涡流面积的影响。从表3和表4可见:当显著性水平α 为0.05 时，因素A和B 均表现出显著性影响(P＜0.05)，表明在本次物理试验中，鱼礁开口比和来流流速对背涡流面积均有显著性影响;当显著性水平α 为0.001 时，因素A表现出极显著性影响(P＜0.001)，而因素B 则无极显著性影响(P＞0.001)。根据上述分析可知，在鱼礁流场效应中，礁体的开口比是影响背涡流面积的主要因素，而流场中来流流速对其影响程度次之。