随着物质生活水平的不断提高，人们对优质水产品的需求与日俱增，陆基与近海养殖存在的环境污染和养殖规模等问题日益凸显[1]，发展以游弋式养殖工船为代表的深远海养殖已成为中国海洋产业发展的新方向之一[2]。自20世纪70年代末，“未来海洋农牧场”的初步设想、大型养殖工船的系统研究[3-4]及深远海养殖平台的总体技术研究[5]逐步开展，养殖工船鱼舱环境的适渔性问题被提出[6]。养殖工船在深远海水域开展水产养殖作业时，受风浪流作用往往会引起船舶的横摇和垂荡等运动，从而不可避免地会发生鱼舱晃荡现象；鱼舱作为鱼类的生活场所，其大幅晃荡可引起舱内流速急剧变化，使养殖鱼类失去游泳能力[7-8]，并对其摄食和生长等造成不利影响[9]，直接影响养殖工船的经济效益。因此，研究晃荡鱼舱内流速特性对优化鱼舱结构和发展工船养殖具有重要意义。
相较于陆基养殖的进出水推流，养殖工船鱼舱水环境受进出水推流与鱼舱晃荡共同作用。目前，关于鱼池养殖环境流速的研究主要集中在陆基养殖[10-12]。Labatut等[13]研究发现，排水速率(池底中心排水)对养殖池内流速大小和均匀性无显著影响，进水喷嘴射流对养殖池内流速大小有显著性影响，但对流速均匀性无显著性影响，流速大小几乎随喷嘴直径和射流速度呈线性变化。然而，Gorle等[12]则认为，采用带有径向分量的进水射流可以改善养殖池内流速的均匀性。面向深远海的养殖工船，外界环境激励下的鱼舱晃荡和进水射流均是其不可忽略的影响因素，其鱼舱内流速特性是鱼舱适渔性的重要判据。在船舶与海洋工程领域，目前有关液舱晃荡问题的研究多集中在液舱晃荡下舱壁载荷和船体运动等安全性问题上[14-16]，而有关液舱晃荡下舱内流速的统计特性鲜有报道。近期，Guo等[17]和Cui等[18]针对养殖工船鱼舱适渔性问题开展了鱼舱水环境适渔性的数值研究，采用概率统计方法分析了鱼舱内流速分布和滞留时间等特性。Li等[19]提出了鱼舱内适渔性保证率，并通过其研究分析了鱼舱水流环境的适渔性。目前，养殖工船鱼舱内流速特性的试验研究尚未见报道。
本研究中，采用实船试验对3 000 t级养殖工船横摇运动下鱼舱内的流速特性进行了研究，重点分析了不同试验工况下横摇运动特征参数对舱内流速的影响，以期为船载鱼舱流速适渔性问题的研究提供实证试验数据。
1 材料与方法
1.1 试验系统
试验在3 000 t级养殖工船“国信101”号(以下简称“工船”)的鱼舱内进行，船全长为93.0 m，型宽为16.6 m，型深为6.6 m，设计吃水为4.6 m。
工船鱼舱长为8.8 m，宽为7.8 m，高为5.2 m，舱内共布置了16个进水口，1个出水口；进水口分别位于A1、A2、A3、A4四处，进水方向如图1(a)所示，出水口位于舱中央底部位置。A1～A4四处，每个位置垂向布置4个进水口，相邻两进水口间的距离为0.7 m，最下端的进水口距舱底1.5 m(距液面2.7 m)，结构尺寸如图1(b)所示。工船全船布置3个结构尺寸相同的鱼舱，分别为T1、T2、T3，其中T3为试验舱(图1(c))。试验过程中，在145 m3/h的总进水流量下，舱内的液位高度维持在4.2 m左右；T3舱内约有1万尾大黄鱼(Larimichthys crocea)，养殖密度约为17.3 kg/m3。
A，A1～A4为进水口；B1～B4为进水方向；C为出水口；T1～T3为鱼舱；P为AHRS的安装位置；单位均为mm。
A and A1-A4 are the water inlets，and B1-B4 are the water inlet directions；C is outlet；T1-T3 are fish tanks；P is the location where AHRS is installed；all units are mm.
图1 工船鱼舱的布置及结构
Fig.1 Arrangement and structure of fish tanks on the vessel
工船试验历时近1个月，工船锚泊于青岛市长门岩锚地(36°15′071″N，120°54′729″E)；锚泊状态下，受风浪流作用工船运动以横摇和垂荡为主，垂荡运动偶尔发生，而横摇运动是常态。
1.2 试验方法
1.2.1 船体横摇角度和流速数据的采集 本文主要研究不同海况下工船船体横摇运动对鱼舱内水体流速特性的影响，其中，船体横摇运动幅度通过横摇角度进行量化，而在不同横摇角度下鱼舱内的水体流动状态通过流速进行描述，故船体横摇角度和流速是本试验需采集的主要数据。船体横摇角度数据通过数字量航姿参考系统(型号AHR730，中国无锡极锐科技有限公司)采集，采样频率为16 Hz；流速数据通过声学多普勒流速剖面仪(型号DCPS 5400，挪威安德拉数据仪器股份有限公司)采集，每2 min输出一组数据，包括水平流速和垂向流速，输出的流速为采集时间内的平均流速。
试验过程中，数字量航姿参考系统(以下简称为AHRS)固定于船的甲板上(图1(c)中P点位置)，实时采集船体横摇角度数据。声学多普勒流速剖面仪(以下简称“ADCP”)倒置悬挂于鱼舱上方，对不同水深处的流速数据进行采集，ADCP的定位安装如图2(a)所示。
L1～L4 依次为 ADCP、定位绳、悬挂安装绳和辅助安装绳；0～12为流速采集测点。
L1-L4 are ADCP，positioning rope，suspension installation rope，and auxiliary installation rope in sequence；0-12 are the measuring points of the flow velocity.
图2 鱼舱流速的采集
Fig.2 Measurement of flow velocity in a fish tank
ADCP对鱼舱内流速的采集存在盲区，液面下0～2 m内流速无法准确采集；同时，根据ADCP的工作原理[20]，应在待测区域内均匀布置测点，以完成该区域内的流速采集。因T3舱属中心对称结构，故为了避免流速采集测点的重复布置，提高试验效率，仅在舱内1/4区域均匀布置测点；如图2(b)所示，鱼舱内共布置了13个流速采集测点，即0～12测点，其中，0测点位于鱼舱中央位置，9～12测点位于舱进水方向上。为便于后续处理，将各测点至舱中央处(即0测点位置)的直线距离记为距舱中心距离，则0～12测点距舱中心的距离依次为0、3.9、2.6、1.3、0.4、4.3、3.1、2.2、1.8、4.6、4.1、3.4、3.2 m。
1.2.2 横摇角度数据的处理 AHRS所采集的船体横摇角度数据是时历数据，其主要特征有横摇角度幅值和横摇周期。横摇角度幅值可直接从时历数据中获取，但横摇周期的获取需对时历数据做一定的处理，该过程包括分帧、加窗和时频域转换。本研究中，横摇角度数据的处理通过MATLAB R2021a完成。
在不同外界激励的作用下，工船船体横摇周期常是变化的，但其仍具有短时平稳的特点，故可对船体横摇角度数据做分帧和加窗处理，截取一小段时历数据，提取其短时特性；同时，为减少帧与帧之间的不连续性，取帧移为帧长的一半。本研究中，根据横摇角度数据的特点，选取帧长为32 s，帧移为16 s。对船体横摇角度数据的处理主要是为了从其中获取横摇周期，故应选择频率分辨率较高的窗函数，即矩形窗。
在分帧和加窗处理后，仍需通过一定的变换，才能将船体横摇角度的时历数据转换到频域，提取其频率特性，从而计算出相应的横摇周期。采用快速傅立叶变换(fast fourier transform，FFT)完成数据处理[21-22]，计算公式为
(1)
式中：x(n)为FFT变换前的时域数据；X(k)为FFT变换后的频域数据；N为时域内的采样点数，应为2的整数次幂；k为频域内相应的点数，k=0，1，2，…，N-1；e为自然常数，e≈2.72；j为复数中虚部的单位；π为圆周率，π≈3.14。
1.2.3 工况分组 通过对船体横摇角度数据的处理可知，在锚泊状态下，工船横摇运动的周期为4～8 s，且多集中在5.3～6.4 s，工船横摇运动的角度幅值为0°～7.3°。因此，本研究中将横摇周期控制在5.3～6.4 s，并根据船体横摇角度幅值对试验工况进行分组，以分析不同试验工况下工船鱼舱内的流速分布，试验工况分组如表1所示。
表1 试验工况分组
Tab.1 Grouping of experimental conditions
工况condition横摇角度幅值rolling angle工况说明conditions descriptionC11°以内within 1°小幅横摇slight rollingC23°左右about 3°中等幅度横摇moderate rollingC35°以上over 5°大幅横摇large rolling
1.2.4 流速数据的处理 基于鱼舱内0～12测点处所测的流速数据，依次绘制C1、C2和C3 3种试验工况下的流速云图，包括水平流速云图和垂向流速云图；基于舱进水方向上的测点流速数据(9～12测点)，绘制3种试验工况下的舱进水方向上水平流速云图。同时，为分析流速随距舱中心距离的变化，基于0～12测点处所测的流速数据，绘制箱线图，在距舱中心距离维度上分析C1、C2和C3 3种试验工况下舱内流速的变化；为分析舱内流速随水深的变化，在各试验工况下，取相同水深所有流速数据的平均值，作为该水深下的流速，并绘制折线图，在水深维度上分析C1、C2和C3 3种试验工况下舱内流速的变化。
1.2.5 数据的归一化处理 为便于直观地比较各试验工况下舱内流速变化的幅度，需对流速数据进行归一化处理[23]。相较于其他工况，C1试验工况的横摇角度幅值在1°以内，更接近于无横摇状态，因此，以C1试验工况下测得的最大流速数据为基础对其他流速数据进行归一化处理。计算公式为
uN=u/u0，
(2)
vN=v/v0。
(3)
式中：u0为C1试验工况下的最大水平流速(cm/s)；u为归一化前的水平流速(cm/s)；uN为归一化后的水平流速数据；v0为C1试验工况下的最大垂向流速(cm/s)；v为归一化前的垂向流速(cm/s)；vN为归一化后的垂向流速数据。流速数据的处理通过MATLAB R2021a软件完成。
2 结果分析
2.1 横摇对鱼舱水平流速的影响
2.1.1 水平流速分布 在C1、C2和C3 3种试验工况下，鱼舱内水平流速为0.5～28.2 cm/s，3种工况下鱼舱内最大水平流速依次为18.1、18.1、28.2 cm/s(表2)。在C1试验工况下，舱内水平流速从舱壁至舱中央呈递减趋势，最大水平流速出现在2.8 m水深附近，且舱中央存在低速区(图3(a))。因鱼舱最下端进水口处的水深为2.7 m(图1(b))，故在C1试验工况下舱内最大水平流速出现的位置可能与进水射流有关，这一点有待进一步研究证实。
表2 不同试验工况下的鱼舱流速
Tab.2 Flow velocity in a fish tank under different conditions
方向direction流速 velocity/(cm·s-1) C1工况C1 conditionC2工况C2 conditionC3工况C3 condition水平horizontal 0.5～18.10.5～18.10.6～28.2舱进水 inflow2.1～18.10.5～17.81.6～14.4垂直vertical -4.3～9.6-4.8～8.7-2.5～8.3
图3 3种试验工况下鱼舱内水平流速分布
Fig.3 Distribution of horizontal flow velocity in the fish tank under three experimental conditions
对比C1和C2试验工况发现，舱内水平流速的最大值变化不大，而最大水平流速的位置发生了变化，随着横摇角度幅值由1°增大至3°，舱内最大水平流速的位置有向舱中部偏移的趋势；对比C2和C3试验工况发现，随着横摇角度幅值由3°增大至5°，舱内水平流速出现了大幅增加，且最大水平流速的位置已偏移至舱中部(图3(a))。由于鱼舱的进水射流在靠近舱壁区域(图1(a))，故C1试验工况下舱内最大水平流速的位置发生变化可能与其有关。对比3种工况发现，随着横摇角度幅值的逐渐增大，横摇运动对舱内水平流速的影响逐步增大，在大幅横摇运动(C工况)下，舱内水体的动能与势能交替转化，舱内流速大幅增加，在舱中部区域也出现了较高流速，故从整体上看，会发现舱内最大水平流速的位置向舱中部偏移。但舱中央位置处始终存在低速区，且随着横摇角度幅值的增大，低速区的范围略有缩小。
另外，观察舱内水体的流动和鱼的游动行为发现，在C1和C2试验工况下，舱内水体绕着鱼舱的垂向中心轴旋转流动，流动方向呈逆时针方向(从浅水区向深水区看)，而鱼成群游动，游动方向与水流方向相反。相较于C2试验工况，在C3试验工况下观察到舱内有部分鱼的游动略显混乱，游动时不再保持队形(图4)。
T—鱼舱；B—水流方向；Fg—成群游动的鱼；Fo—离群游动的鱼。
T—fish tank；B—direction of the water flow；Fg—fish swimming in groups；Fo—fish swimming out of the group.
图4 3种试验工况下鱼舱水深处鱼的游动状态
Fig.4 Swimming state of fish at the deep area of the fish tank under three experimental conditions
2.1.2 水平流速随距舱中心距离的变化 对舱内的水平流速数据进行归一化，即取u0=18.1 cm/s。对比C1和C2试验工况发现，随着横摇角度幅值的增大，舱内水平流速整体有所增大，且最大水平流速的位置有向舱中部偏移的趋势；在C2试验工况下，除4.1 m测点外，3.1 m和2.2 m测点处均出现了较大的水平流速；C3试验工况下，舱内水平流速相对于其他工况出现了大幅的增加，最大水平流速的增幅达56.2%(较u0)，且最大水平流速的位置已偏移至舱中部(图3(b))。
2.1.3 水平流速随水深的变化 在各试验工况下，取相同水深的所有水平流速数据的平均值，作为该水深下的水平流速，C1、C2和C3 3种试验工况下所有深水处的最大水平流速依次为6.8、8.1、10.4 cm/s。对舱内的水平流速数据进行归一化，即取u0=6.8 cm/s。
对比C1和C2试验工况发现，随着横摇角度幅值的增大，深水区的水平流速变化不大，而浅水区的水平流速出现了一定幅度的增大，C2试验工况下最大水平流速的增幅为19.7%(较u0)；相对于其他工况，C3试验工况下舱内水平流速均出现了大幅的增加，最大水平流速的增幅达52.6%(较u0)(图5)。此外还发现，在不同试验工况下，舱内浅水区的水平流速均明显高于深水区。
图5 3种试验工况下水平流速随水深的变化
Fig.5 Variation in horizontal flow velocity with water depth under three experimental conditions
2.1.4 舱进水方向上的水平流速分布 3种试验工况下，舱进水方向上的水平流速整体为0.5～18.1 cm/s，各工况下最大水平流速依次为18.1、17.8、14.4 cm/s(表2)。对比3种试验工况发现，舱内水平流速分布(图3(a))和进水方向上的水平流速分布(图6)可知，在C1试验工况下，鱼舱内较大的水平流速均出现在进水方向上；随着横摇角度幅值的增大，晃荡与射流共同作用使得舱进水方向上的水平流速整体呈减小，最大水平流速的位置略向深水区偏移。
图6 舱进水方向上水平流速云图
Fig.6 Horizontal flow velocity contour along the inflow direction
2.2 横摇对鱼舱垂向流速的影响
2.2.1 垂向流速分布 在C1、C2和C3 3种试验工况下，舱内垂向流速的范围依次为-4.3～9.6、-4.8～8.7、-2.5～8.3 cm/s(表2)，其中，正值表示流速方向向下，即从浅水区至深水区，负值则相反。在C1试验工况下，舱中央区域的垂向流速整体呈向上趋势，而靠近舱壁区域的垂向流速整体呈向下，这可能与舱内形成的“二次流”有关。此外还发现，靠近舱底部区域的垂向流速均向下，这应与舱底中央处的出口排水有关。对比3种试验工况发现，随着横摇角度幅值的增大，舱壁位置处向下的垂向流速区域逐步扩大，而舱中央位置处向上的垂向流速区域逐步缩小；在C3试验工况下，舱中央区域内向上的最大垂向流速出现了大幅的降低，从4.8 cm/s降至2.5 cm/s(图7(a))。
图7 3种试验工况下鱼舱内垂向流速分布
Fig.7 Distribution of vertical flow velocity in fish tank under three experimental conditions
2.2.2 垂向流速随距舱中心距离的变化 对舱内的垂向流速数据进行归一化，即取v0=9.6 cm/s。对比3种试验工况发现，舱内垂向流速的变化主要出现在近舱壁区域和舱中央区域；随着横摇角度幅值的增大，近舱壁区域的垂向流速出现了一定幅度的增大，而舱中央区域内向上的垂向流速明显减小，尤其在C3试验工况下(图7(b))。
2.2.3 垂向流速随水深的变化 在各试验工况下，取相同水深的所有垂向流速数据的平均值，作为该水深下的垂向流速，C1、C2和C3 3种试验工况下最大的垂向流速依次为5.2、5.2、5.9 cm/s。对舱内的垂向流速数据进行归一化，即取v0=5.2 cm/s。在3种试验工况下，垂向流速均为正值，即流速方向为向下；随着水深的增加，垂向流速呈先减小后增大的趋势，但整体上垂向流速呈增大趋势。由于鱼舱最下端的进水口在2.7 m水深处，出水口在舱底中央处(图1)，舱内水体从进水口流入，并逐渐从出水口排出，因此，对于2.7～4.0 m水深区域内的水体，在进水射流和出口排水的作用下，其垂向流速随着水深的增加而增大；但对于2.2～2.7 m水深区域内的水体，其在进水口上方且远离出水口，故该区域内垂向流速变化不明显，尤其在C3试验工况下；对比3种试验工况发现，随着横摇角度幅值的增大，垂向流速增大，C3试验工况下最大垂向流速的增幅为12.1%(较v0)(图8)。
图8 3种试验工况下垂向流速随水深的变化
Fig.8 Variation in vertical flow velocity with water depth under three experimental conditions
3 讨论
3.1 不同试验工况下的适渔性
衡量工船鱼舱是否能为养殖对象提供适宜的养殖环境，流速是一个重要的物理环境指标。流速引起的强制运动可以改善养殖鱼类的肌肉张力，提高其耐力，从而提高鱼的存活率[24]，但过高的流速对鱼类也是不利的[7]。养殖鱼类的适养流速与其体长有关，对于1龄以上的大黄鱼(体长大于27 cm)，流速在50 cm/s以下时对其的不利影响较小[19，25]。本研究中通过工船现场实测，在小幅横摇(横摇角度幅值在1°以内)、中等幅度横摇(横摇角度幅值在3°左右)和大幅横摇(横摇角度幅值在5°以上)工况下采集舱内流速数据，发现3种试验工况下鱼舱采集区域内的流速(水平流速和垂向流速的合速度)整体为0～28.4 cm/s，小于50 cm/s；同时还发现，在小幅和中等幅度横摇工况下舱内鱼成群逆水流方向游动，在大幅横摇工况下舱内有部分鱼的游动略显混乱，但总的来说，3种试验工况下舱内鱼的游动是呈队列的。因此，对于锚泊状态下的工船，当横摇角度幅值在5°以内时鱼舱流速适合1 龄以上大黄鱼(体长大于27 cm)的养殖；当横摇角度幅值大于5°时鱼舱流速是否会对大黄鱼生长产生影响还需进一步试验研究。
3.2 不同试验工况下的流速分布
与陆基养殖不同，养殖工船鱼舱内流速分布受进出水推流与鱼舱晃荡共同作用。本研究在锚泊状态下的养殖工船上进行，重点分析不同试验工况下横摇运动的角度幅值对鱼舱内流速的影响，结果发现，在小幅横摇(横摇角度幅值在1°内)工况下，鱼舱内水体呈逆时针方向(从浅水区向深水区看)旋转流动，舱内水平流速从舱壁至舱中央呈递减趋势，舱中央存在低速区，这一规律与Labatut等[10，13，26]针对进出水推流下MCR跑道池(mixed-cell raceway)内流速分布的研究结果一致；综合对比不同试验工况下舱内水平流速分布(图3(a))和舱进水方向上的水平流速分布(图6)可知，小幅横摇(横摇角度幅值在1°内)工况下可仅考虑进出水推流，忽略船体横摇运动对舱内水平流速的影响，但在大幅横摇(横摇角度幅值在5°以上)工况下，相对于进出水推流，船体横摇运动对舱内水平流速的影响更显著。同时，相关研究[12-13]也表明，出口排水(池底中心排水)对养殖池内流速大小和均匀性的影响均不显著，而进水喷嘴射流的影响是主要的。本研究中，对比小幅、中幅和大幅横摇工况下舱内水平流速分布还发现，横摇周期一定时，随着横摇角度幅值的增大，舱内水平流速增大，最大水平流速的位置向舱中部偏移，流速大小随横摇角度幅值变化的规律在Guo等[17]和崔铭超等[27]的研究中也被发现。
本研究中的试验鱼舱是类圆柱体结构，进水方向为其切线方向(图1(a))，故在舱内会形成逆时针旋转(从浅水区至深水区看)的主流及叠加于主流之上的二次流[28]，从而会在舱中央区域形成向上的垂向流速，在近舱壁区域形成向下的垂向流速，这与在小幅横摇工况下舱内垂向流速分布相吻合(图7(a))。对比3种试验工况下舱内的垂向流速分布还发现，随着船体横摇角度幅值的增大，舱内垂向流速分布有所变化，尤其在大幅横摇工况下，舱中央位置处向上的垂向流速区域大幅缩小，向上的最大垂向流速大幅降低。这一变化会促进舱中央区域内生物固体的沉降，更有利于鱼舱的集排污。
3.3 本试验结果与相关数值仿真的对比分析
在鱼舱晃荡问题的研究中，数值仿真和实船试验是两种常用的方法，但因部分技术问题尚未解决，目前无法在舱内流场的数值仿真中考虑鱼的存在[17，27，29]。Guo等[17]和崔铭超等[27]针对横摇运动下养殖工船鱼舱水环境适渔性问题进行的数值研究发现，在进出水推流工况下舱内流速为0～12.0 cm/s，这与本研究中的小幅横摇(横摇角度幅值在1°内)工况下测得的舱内流速(0.5～18.1 cm/s)相当；同时，Guo等[17]发现在横摇工况(横摇周期10 s，横摇角度幅值在5°以上)下舱内流速为0～44.0 cm/s，这与本研究中在大幅横摇工况(横摇周期5.3～6.4 s，横摇角度幅值在5°以上)下所测得的舱内流速(0.6～28.2 cm/s)相差较大。Guo等[17]数值仿真研究所得的舱内流速是在恒定激励下的稳定结果，而本研究中的现场试验环境多变，无法获得某一特定激励下的稳定结果，只能获取某一时间段下的平均流速，故在相同横摇工况下所获得的舱内流速应较数值仿真下所得的流速低；此外，本研究中的实船试验是在工船养殖现状下进行的，鱼舱内有鱼的存在，其对舱内流速的影响应予以考虑。根据Masaló等[11]和Plew等[23]的研究结果，鱼的存在会降低舱内的流速，并且随着放养密度的增加，水槽内流速降低的幅度逐步增大；在35.6、79.4 kg/m3放养密度下，鱼池内的平均流速依次降低了15%和57%。本研究中，鱼舱内鱼的放养密度约为17.3 kg/m3，这会进一步地降低工船鱼舱内的流速。横摇角度幅值和横摇周期均是横摇运动下养殖工船鱼舱内流速的重要影响因素，因此，在相关的数值仿真研究[17，27，30]中常重点考虑，但在本研究中，无法在控制横摇角度的前提下分析横摇周期对舱内流速的影响。实际海况下船体横摇运动的周期与船舶固有属性有关，无法因外界激励的作用而发生较大的变化，故本研究中船体横摇周期变化范围较小，多为5.3～6.4 s。
4 结论
1)小幅横摇工况(横摇角度幅值在1°以内)下可仅考虑进水射流，忽略船体横摇运动对舱内水平流速的影响；但在大幅横摇工况(横摇角度幅值在5°以上)下，相对于进水射流，船体横摇运动对舱内水平流速的影响更显著。横摇周期一定时，随着横摇角度幅值的增大，舱内水平流速增大，最大水平流速的位置向舱中部偏移。
2)小幅横摇工况(横摇角度幅值在1°以内)下，舱中央区域的垂向流速方向整体呈向上，近舱壁区域的垂向流速方向整体呈向下，这可能与舱内形成的二次流有关。横摇周期一定时，随着横摇角度幅值的增大，舱内垂向流速分布有所变化；在大幅横摇工况(横摇角度幅值在5°以上)下，舱中央位置处向上的垂向流速区域大幅缩小，向上的最大垂向流速大幅降低，从4.8 cm/s降至2.5 cm/s。这一变化会促进舱中央区域内生物固体的沉降，更有利于集排污。
3)对于锚泊状态下的工船，当横摇角度幅值在5°以内时鱼舱流速适合1 龄以上大黄鱼(体长大于27 cm)的养殖。
[1] 唐启升,丁晓明,刘世禄,等.我国水产养殖业绿色、可持续发展战略与任务[J].中国渔业经济,2014(1):6-14.
TANG Q S,DING X M,LIU S L,et al.Strategy and task for green and sustainable development of Chinese aquaculture[J].Chinese Fisheries Economics,2014(1):6-14.(in Chinese)
[2] 余东华.“十四五”期间我国未来产业的培育与发展研究[J].天津社会科学,2020,11(3):12-22.
YU D H.The cultivation and development of China’s future industries during the 14th five-year-plan[J].Tianjin Social Sciences,2020,11(3):12-22.(in Chinese)
[3] 徐皓,江涛.我国离岸养殖工程发展策略[J].渔业现代化,2012,39(4):1-7.
XU H,JIANG T.Development strategy of offshore aquaculture engineering in China[J].Fishery Modernization,2012,39(4):1-7.(in Chinese)
[4] 徐皓,谌志新,蔡计强,等.我国深远海养殖工程装备发展研究[J].渔业现代化,2016,43(3):1-6.
XU H,CHEN Z X,CAI J Q,et al.Research on the development of deep sea aquaculture engineering equipment in China[J].Fishery Modernization,2016,43(3):1-6.(in Chinese)
[5] 蔡计强,张宇雷,李建宇,等.10万吨级深远海养殖平台总体技术研究[J].船舶工程,2017,39(sup 1):198-203.
CAI J Q,ZHANG Y L,LI J Y,et al.General technology research of 100 thousand ton deep sea aquaculture platform[J].Ship Engineering,2017,39(sup 1):198-203.(in Chinese)
[6] 崔铭超,金娇辉,黄温赟.养殖工船系统构建与总体技术探讨[J].渔业现代化,2019,46(2):61-66.
CUI M C,JIN J H,HUANG W Y.Discussion about system construction and general technology of aquaculture platform[J].Fishery Modernization,2019,46(2):61-66.(in Chinese)
[7] 郑金秀,韩德举,胡望斌,等.与鱼道设计相关的鱼类游泳行为研究[J].水生态学杂志,2010,31(5):104-110.
ZHENG J X,HAN D J,HU W B,et al.Fish swimming performance related to fishway design[J].Journal of Hydroecology,2010,31(5):104-110.(in Chinese)
[8] 赵希坤,韩桢锷.鱼类克服流速能力的试验[J].水产学报,1980,4(1):31-37.
ZHAO X K,HAN Z E.Experiments on the current overcoming ability of some freshwater fishes[J].Journal of Fisheries of China,1980,4(1):31-37.(in Chinese)
[9] 张宇雷,管崇武.船载摇摆胁迫对斑石鲷血液生化指标的影响研究[J].中国农学通报,2017,33(29):145-149.
ZHANG Y L,GUAN C W.Ship rocking affecting serum biochemical indexes of spotted knifejaw(Oplegnathus punctatus)[J].Chinese Agricultural Science Bulletin,2017,33(29):145-149.(in Chinese)
[10] LABATUT R A,EBELING J M,BHASKARAN R,et al.Modeling hydrodynamics and path/residence time of aquaculture-like particles in a mixed-cell raceway (MCR) using 3D computational fluid dynamics (CFD)[J].Aquacultural Engineering,2015,67:39-52.
[11] MASAL
I,OCA J.Influence of fish swimming on the flow pattern of circular tanks[J].Aquacultural Engineering,2016,74:84-95.
[12] GORLE J M R,TERJESEN B F,SUMMERFELT S T.Hydrodynamics of Atlantic salmon culture tank:effect of inlet nozzle angle on the velocity field[J].Computers and Electronics in Agriculture,2019,158:79-91.
[13] LABATUT R A,EBELING J M,BHASKARAN R,et al.Effects of inlet and outlet flow characteristics on mixed-cell raceway (MCR) hydrodynamics[J].Aquacultural Engineering,2007,37(2):158-170.
[14] SOUTO-IGLESIAS A,BULIAN G,BOTIA-VERA E.A set of canonical problems in sloshing.part 2:influence of tank width on impact pressure statistics in regular forced angular motion[J].Ocean Engineering,2015,105:136-159.
[15] 肖凯隆,陈作钢.多液舱晃荡与养殖工船时域耦合运动的数值模拟[J].中国舰船研究,2020,15(1):136-144.
XIAO K L,CHEN Z G.Numerical simulation of aquaculture ship motions coupled with tanks sloshing in time domain[J].Chinese Journal of Ship Research,2020,15(1):136-144.(in Chinese)
[16] ZHUANG Y,WAN D C.Numerical study on ship motion fully coupled with LNG tank sloshing in CFD method[J].International Journal of Computational Methods,2019,16(6): 1840022.
[17] GUO X Y,LI Z S,CUI M C,et al.Numerical investigation on flow characteristics of water in the fish tank on a force-rolling aquaculture platform[J].Ocean Engineering,2020,217:107936.
[18] CUI M C,LI Z S,ZHANG C L,et al.Statistical investigation into the flow field of closed aquaculture tanks aboard a platform under periodic oscillation[J].Ocean Engineering,2022,248:110677.
[19] LI Z S,GUO X Y,CUI M C.Numerical investigation of flow characteristics in a rearing tank aboard an aquaculture vessel[J].Aquacultural Engineering,2022,98:102272.
[20] BRUMLEY B H,CABRERA R G,DEINES K L,et al.Performance of a broad-band acoustic Doppler Current profiler[J].IEEE Journal of Oceanic Engineering,1991,16(4):402-407.
[21] AYINALA M,PARHI K K.FFT architectures for real-valued signals based on Radix-23 and Radix-24 algorithms[J].IEEE Transactions on Circuits and Systems I:Regular Papers,2013,60(9):2422-2430.
[22] 杨丽娟,张白桦,叶旭桢.快速傅里叶变换FFT及其应用[J].光电工程,2004,31(sup 1):1-3,7.
YANG L J,ZHANG B H,YE X Z.Fast Fourier transform and its applications[J].Opto-Electronic Engineering,2004,31(sup 1):1-3,7.(in Chinese)
[23] PLEW D R,KLEBERT P,ROSTEN T W,et al.Changes to flow and turbulence caused by different concentrations of fish in a circular tank[J].Journal of Hydraulic Research,2015,53(3):364-383.
[24] BURROWS R E,CHENOWETH H H.The rectangular circulating rearing pond[J].The Progressive Fish-Culturist,1970,32(2):67-80.
[25] 吴飞飞,王萍,桂福坤,等.大黄鱼续航时间和临界游泳速度的初步研究[J].渔业现代化,2014,41(4):29-33.
WU F F,WANG P,GUI F K,et al.Preliminary study on sustained swimming time and critical swimming speed of Pseudosciaena crocea[J].Fishery Modernization,2014,41(4):29-33.(in Chinese)
[26] LABATUT R A,EBELING J M,BHASKARAN R,et al.Exploring flow discharge strategies of a mixed-cell raceway (MCR) using 2-D computational fluid dynamics (CFD)[J].Aquacultural Engineering,2015,66:68-77.
[27] 崔铭超,王靖,郭晓宇.横摇运动下养殖工船水环境流场特性数值分析[J].中国造船,2020,61(3):204-215.
CUI M C,WANG J,GUO X Y.Numerical analysis of flow field characteristics in a tank of aquaculture vessel under rolling motion[J].Shipbuilding of China,2020,61(3):204-215.(in Chinese)
[28] 湛含辉.二次流现象及其初步研究[J].株洲工学院学报,2001,15(3):27-29.
ZHAN H H.Research and phenomenon of secondary flow[J].Journal of Zhuzhou Institute of Technology,2001,15(3):27-29.(in Chinese)
[29] 方舟,徐红昌,戴锦阳.横摇运动下养殖工船多液舱流场特性数值分析[J].船海工程,2021,50(6):41-45,50.
FANG Z,XU H C,DAI J Y.Numerical analysis of flow field characteristics in muti-tanks of aquaculture vessel under rolling motion[J].Ship & Ocean Engineering,2021,50(6):41-45,50.(in Chinese)
[30] 徐博,姜胜超,刘浩,等.基于STAR-CCM+的液舱晃荡数值研究[C]//第十九届中国海洋(岸)工程学术讨论会论文集.重庆:中国海洋工程学会,2019:178-183.
XU B,JIANG S C,LIU H,et al.Numerical investigation on tank sloshing based on STAR-CCM+[C]//The 19th China Ocean (Shore) Engineering Symposium.Chongqing:China Marine Engineering Society,2019:178-183.(in Chinese)