Nutrional analysis and quality evaluation of Grateloupia livida
-
摘要: 为评价舌状蜈蚣藻(Grateloupia livida)的营养价值,提高其综合利用率,对其基本营养成分、氨基酸、脂肪酸及矿物元素进行了测定。结果显示,粗蛋白、粗脂肪、灰分和粗纤维的质量分数分别为12.28%、0.18%、13.32%和1.80%;舌状蜈蚣藻含有17种氨基酸,包括8种必需氨基酸,限制性氨基酸为色氨酸;共含有15种脂肪酸,包括具有重要生理功能的亚油酸、二十碳四烯酸、二十碳五烯酸和二十二碳六烯酸;舌状蜈蚣藻富含镁(Mg)、钙(Ca)、铁(Fe)、锌(Zn)和铜(Cu)。结果表明舌状蜈蚣藻具有较好的营养保健价值。Abstract: In order to evaluate the nutritive value of Grateloupia livida and improve its comprehensive utilization rate, we determined the composition of the basic nutrients, amino acids, fatty acids and minerals of G. livida. The results show that the mass fractions of crude protein, crude fat, ash and crude fiber were 12.28%, 0.18%, 13.32% and 1.80%, respectively. Seventeen kinds of amino acids were found which contained eight kinds of essential amino acids, and the limiting amino acid was tryptophan. G. livida contains 15 kinds of fatty acids, such as linoleic acid, arachidonic acid, eicosapentaenoic acid and docosahexaenoic acid which have important physiological functions. The results indicate that G. livida is nutritive and rich in Mg, Ca, Fe, Zn and Cu.
-
Keywords:
- Grateloupia livida /
- nutritional composition /
- nutrition evaluation
-
中国是世界第一水产养殖大国。为解决近岸养殖面积不足、病害频发等问题,海水养殖产业正在积极向深远海拓展[1-2]。在实现深远海养殖的过程中,大型养殖工船是不可或缺的工程装备,从20世纪70年代至今,养殖工船技术经历了快速的发展[3-6]。相对于传统的陆基养殖池,大型养殖工船具有游弋避灾、稳定产出、集约高效、产品优质等优点[7-9]。养殖舱是养殖工船重要的组成部分,其内的流场环境对养殖产量起着至关重要的作用,因此通过优化养殖舱的设计来改善养殖环境、提升鱼类适养水平,逐步成为研究热点[10-11]。
有关养殖工船的船体运动以及养殖舱内水动力特性已有一些研究报道。方舟等[12]对横摇运动下养殖工船多液舱流场特性进行了数值分析,指出将养殖舱沿船中对称布置可有效降低舱内流体的晃荡程度,获得更优的水体环境。苑中排等[13]研究了液舱晃荡与船舶的耦合运动,结果表明液舱晃荡对浮式生产储油船 (Floating production storage and offloading, FPSO) 在风浪下的垂向运动基本无影响。崔铭超等[14]研究发现,方形养殖舱在横摇运动条件下,横摇周期与角度对舱内流速分布有显著影响。Guo等[15]模拟分析了横摇运动条件下养殖舱内流体均匀性、水体停留时间等水力特性,发现流速与横摇角度呈现正相关关系,而水体停留时间受到的影响并不明显。高瑞等[16]通过对横摇运动下适渔性的研究发现,板式隔板可有效地抑制水体晃荡;且在总进水流量一定的条件下,随着进水口数量的增加,舱内流速逐渐下降[17],随着出水口数量的减少,舱内流速上升[18]。以上研究均在养殖工船横摇晃荡的条件下进行。然而,在养殖作业过程中,系泊状态下的养殖工船运动幅度较小,舱内的流场环境主要由进出水口结构、舱体形状等因素决定,与陆基鱼池存在一定的相似性[19]。
Klebert等[20]通过实验与计算流体力学 (Computational Fluid Dynamics, CFD) 数值模拟相结合的方法探究了陆基循环水养殖池的水动力特性。常见的陆基养殖池多为矩形和圆形,研究表明矩形养殖池的长宽比保持在1.43~1.91内,可有效减少死水区的产生[21];相较于矩形池,圆形养殖池在水体均匀性和流动稳定性方面更佳,但空间利用率较低。养殖池内的水体流动主要依赖于进水射流的推动,因此探究进水方式对流场的影响具有重要意义。于林平等[22]研究了不同射流孔数和进水管位置对矩形圆弧角养殖池内流场的影响,发现将进水管布置在圆弧角处有利于改善底部流场环境,且随着射流孔数的增加 (日循环次数增加),出水口处的流速明显增大,这有利于提高集排污效率。将传统的垂直进水方式改为水平进水有利于减少低速漩涡区,增加水流均匀性,进而避免水体环境持续恶化[23-24]。进水射流的方向和角度对舱内流场环境也有重要影响。Zhang等[25]分析了不同进水管数量下养殖池内的水动力特性,结果表明,随着进水管数量的增多,进水射流之间的相互影响增强,形成了复杂的湍流,这降低了流场的均匀性。此外,进水角度对养殖舱内的流场也有显著性影响[26-28]。
与陆基鱼池通过直立柱提供进水相比,工船养殖舱大多将进水口侧置于舱体切角处,而且水深往往更深,二者的流场环境存在差异,因此探究静止条件下养殖舱内的流场变化规律非常必要。本研究以30万吨养殖工船矩形养殖舱为研究对象,其养殖舱体积较大,长宽比接近1.43,并利用CFD技术建立三维数值模型,对不同进水方式下养殖舱内的流场特性进行研究;进一步通过对底部流场的分析评估了大黄鱼 (Pseudosciaena crocea) 的适养性水平,以期为超大型养殖平台养殖舱结构设计优化提供理论指导。
1. 材料与方法
1.1 数值模型
本研究将养殖工船养殖舱内的流体视为三维不可压缩流体,该流体的控制方程为连续性方程、纳维-斯托克斯(N-S)方程。
不可压缩流体的连续性方程为:
$$ \frac{\mathrm{\mathrm{\partial}}u}{\partial x}\mathrm{+}\frac{\partial v}{\partial y}\mathrm{+}\frac{\partial w}{\partial z}\mathrm{=}0 $$ (1) 式中:$ (u,v,w) $ 为坐标轴 $ (x,y,{ {z}}) $ 上的3个速度分量。
纳维-斯托克斯(N-S)方程为:
$$ \frac{{{\mathrm{d}}v}}{{{\mathrm{d}}t}} {\mathrm{=}} \frac{{\partial v}}{{\partial t}} {\mathrm{+}} \left( {v \cdot \nabla } \right)v {\mathrm{=}} {\text{−}} \frac{1}{\rho }\nabla P {\mathrm{+}} \mu {\nabla ^2}v $$ (2) 式中:$ v $ 为速度 (m·s−1);$ t $ 为时间 (s);$ \rho $ 为流体密度 (kg·m−3);P为压强 (Pa);$ \mu $ 为黏度系数 (m2·s−1)。
本研究中,养殖舱壁处的边界层效应和流体旋转过程的湍流现象不可忽略。RNG k-ε湍流模型具有广泛的适用性和较高的准确性,且充分考虑了低雷诺数对湍流结构的影响,提高了复杂涡旋流场的计算精度,特别适用于描述大应变率剪切流、分离流和旋转流量[29]。方程为:
$$ \begin{split} &\begin{gathered} \frac{{\partial {k_T}}}{{\partial t}} {\mathrm{+}} \frac{1}{{{V_F}}}\left( {u{A_x}\frac{{\partial {k_{_T}}}}{{\partial x}} {\mathrm{+}} v{A_y}\frac{{\partial {k_T}}}{{\partial y}} {\mathrm{+}} w{A_{z}}\frac{{\partial {k_T}}}{{\partial {z}}}} \right) {\mathrm{=}} \\ {P_T} {\mathrm{+}} {D_{{k_T}}} {\text{−}} {\varepsilon _T}\end{gathered}\\[-8pt]& \end{split} $$ (3) $$ \begin{split} &\begin{gathered}\frac{{\partial {\varepsilon _T}}}{{\partial t}} {\mathrm{+}} \frac{1}{{{V_F}}}\left( {u{A_x}\frac{{\partial {\varepsilon _T}}}{{\partial x}} {\mathrm{+}} v{A_y}\frac{{\partial {\varepsilon _T}}}{{\partial y}} {\mathrm{+}} w{A_w}\frac{{\partial {\varepsilon _T}}}{{\partial {z}}}} \right) {\mathrm{=}}\\ {c_1}\frac{{{P_T}{\varepsilon _T}}}{{{k_T}}} {\mathrm{+}} {D_{{\varepsilon _T}}} {\text{−}} {c_2}\frac{{\varepsilon _T^2}}{{{k_T}}} \end{gathered}\\[-16pt]& \end{split}$$ (4) 式中:kT为湍流动能;(Ax, Ay, Az)分别为在坐标轴(x, y, z) 方向上流动的分数面积;VF为分数体积;$ {P_T} $ 为湍动能生成项;$ {D_{{k_T}}} $ 为湍动能扩散相;${\varepsilon _T} $ 为湍流耗散率;$ {D_{{\varepsilon _T}}} $ 为湍流耗散扩散相;$ {c_1} $、$ {c_2} $ 为系数,在经典的 k-ε模型中 $ {c_1} $=1.44、$ {c_2} $=1.92,在重整群模型中 $ {c_1} $=1.42,$ {c_2} $ 由湍动能及湍流生成项计算。
1.2 模型试验验证
本研究基于Oca和Masaló[30]和任效忠等[31]的试验结果开展数值模型仿真验证。模型采用圆形养殖池,直径为0.49 m,初始水深为0.06 m,底部排水口的直径为0.01 m,进水口的直径为0.007 m,进水流速为0.47 m·s−1。如图1-a所示,整个模型利用四面体网格划分,在进水管、进水口、出水口的位置进行网格加密,网格的数量级约为55万。在距离池底高度0.03 m处的截面进行流速监测,将数值模拟结果与试验结果平均值进行对比分析。
如图1-b所示,模拟结果与试验结果的变化规律几乎相同,整体吻合性较好,误差基本保持在10%以内,其中最大误差处是由于进水射流的影响造成,精度满足本次仿真的要求,因此该模型可以用于工船船载舱养的数值仿真计算。
1.3 计算模型
本研究的养殖工船,船长342.6 m、宽66 m、深32 m、设计吃水18.5 m (图2-a)。其养殖舱模型如图2-b所示,养殖舱长42 m、宽29.2 m、高19.5 m;初始水位为16 m,水体体积约为
18562 m3。为提高舱内流场的均匀性,全舱16个进水口分为4层布置在壁面拐角处,每层的间隔为3.5 m。在底部等间距设置4个出水口;进水口半径为0.395 m,出水口半径为0.778 m。1.4 数值仿真试验分组
本研究采用控制变量的方法探究日换水次数、进水口数量和进水口布设位置对舱内流场的影响,试验工况如表1所示。
表 1 仿真试验分组Table 1. Grouping of simulation tests工况组别
Working condition
group日换水次数
Number of water change
per day进水口数量
Number of
water inletsA1 16 16 A2 14 A3 12 A4 10 B1 16 16 B2 12 B3 8 B3 4 C1 16 8 C2 C3 养殖舱工作时,进水口与出水口同时打开,在保证初始液位高度不变的条件下,完成1次水体更新即为1次换水,A类工况的进水口布设位置和进水方向如图3-a所示。B类工况的进水口布设位置和进水方向与A类工况相同,并均匀的分布在4个切角处。在研究进水口布设位置时,受MCR跑道池的启发,将4个进水口按照不同的进水方向布置在养殖池的内侧形成C2和C3工况,进水口位置和进水方向如图3-b、3-c所示。
1.5 边界条件和网格划分
定义养殖舱内水面为自由液面,压力设为标准大气压,舱壁和舱底定义为无滑移壁边界,湍流动能和湍流耗散率设置为二阶迎风格式。基于广义最小余量法(Generalized minimal residual method, GMRES) 的速度-压力求解器进行速度和压力的迭代求解,自由液面的模拟采用单相的流体体力分数 (Volume of fluid, VOF) 模型。进出水口分别采用源和汇模拟,时间步长采用自适应时间步长。
本研究采用A1工况进行网格无关性验证,以六面体网格对养殖舱模型进行划分 (图4-a),且由1.2的模型验证可知,进出水口的高流速容易与现实情况产生较大误差,因此需在进出水口处进行局部加密。最小网格尺寸均为0.041 m,最大网格尺寸分别为0.4、0.3、0.2 m,网格数量分别为216 889、
1004621 、5481613 。在距舱底z=1.2 m的平面上选取7个监测点,流速如图4-b所示。由图可见,粗糙网格的仿真结果与精细网格的仿真结果有一定的差距,而0.3 m的中等网格与0.2 m的精细网格在精度上比较接近,综合考虑计算精度和计算时长,本研究的模拟仿真均采用0.3 m的网格尺寸。1.6 数据统计处理
将基于有限体积法输出的速度数据作为集合{φi},并以Δφ的跨度进行分组(本研究取Δφ=0.01 m·s−1),所有网格的体积数据作为集合{Vi}。为更全面地分析流场状态,本研究使用概率密度函数 (Probability density function, PDF) 和累积分布函数 (Cumulative distribution function, CDF) 对流速数据进行统计分析[32]。对于数据 {φi},设定一速度值为θ,则流场中不大于此值的概率表示为:
$$ {F}\left(\varphi \leqslant \theta \right){\mathrm{=}}\frac{{\sum }_{i{\mathrm{=}}1}^{n}{V}_{i}\left({\varphi }_{i} \leqslant \theta \right)}{{\sum }_{i{\mathrm{=}}1}^{n}{V}_{i}} $$ (5) 并且此值的PDF值表示为:
$$ {f}\left(\varphi {\mathrm{=}}\theta \right){\mathrm{=}}{F}\left(\theta {\mathrm{+}}\frac{\Delta \varphi }{2}\right){\text{−}}{F}\left(\theta {\text{−}}\frac{\Delta \varphi }{2}\right) $$ (6) 使用高斯拟合的方法计算标准差ω和均值μ,以便更加直观地比对各个流场之间的差异。
2. 结果与分析
2.1 进水方式对舱内流场的影响
2.1.1 不同进水方式下的流速分布
计算时长统一设为7 000 s,采用CFD软件监测计算网格内流体的平均动能和平均湍动能,发现6 000 s后计算趋于稳定状态。图5展示了不同深度截面在第7 000秒时的速度分布。通过观察可知:在进水推流的作用下,高流速区域集中在靠近舱壁处,低流速区域则主要分布在舱体中心以及切角处,此种流场特征与Zhang等 [25]、崔铭超等[14]的研究结果基本一致。对于C2工况,在进水射流的相互影响下难以形成稳定的流场,流速分布比较紊乱,水体旋转流动的特性并未显现;C3工况则以底部2对出水口为中心形成了较为稳定的对称流场。通过对比分析图5中A类和B类工况在同一高度的流速分布,发现减少日换水次数或增加进水口数量均会导致中心低流速区域的扩大。
为了更准确地揭示舱内流速的分布规律,对整个舱内水体的流速进行了定量统计分析。流速的概率密度函数和累积分布函数如图6所示。图6-a显示舱内流速与日换水次数之间存在明显的负相关性。随着日换水次数的减少,流速均值μ下降了47.3%,标准差ω下降了37.6%,表明流速分布趋向集中。另一方面,随着进水口数量的减少,流速均值上升了30.7%,标准差上升了66.3% (图6-b)。C1和C3工况的流速分布大致相同,而C2工况较低流速区域的占比更高(图6-c)。
由于养殖舱进水口设置在舱壁切角处,其附近水体受到进水推流的影响较大,因此会在舱壁处产生较大的流速,随着输入养殖舱的能量在抵抗水体质点之间的相对运动产生损耗时,远离舱壁的流速逐渐降低,并且切角处的水流撞击也会有能量损失,进而在切角处产生低速区。随着日换水次数的减少,单位时间内输入养殖舱的能量降低,舱内水流难以维持较高的流速,造成低流速区域逐渐扩大;在进水流量一定的条件下,减少进水口数量,使进水流速大幅提升,舱内整体流速有所提高。但是相对于改变日换水次数,进水口数量的减少对流速的影响并不显著。在改变进水口布设位置时,对于C2工况,进水射流之间的相互影响造成了能量的损耗,整个养殖舱处在低流速水平;对于C3工况,在进水射流的影响下左右两侧舱体形成反方向转动的水流,而整体则呈“∞”型流动,高流速区域出现在两侧舱壁和舱体中心。
2.1.2 不同进水方式下的流场均匀性
流体均匀性是评价养殖舱内水体环境的重要指标,均匀的水流有利于鱼类的生长。为了量化这一指标,本研究引入流动均匀性指数,其计算公式为[33]:
$$ \gamma {\mathrm{=}}1{\text{−}}\frac{1}{2S}\sum _{i{\mathrm{=}}1}^{n}\frac{\sqrt{{\left({V}_{i}{\text{−}}{V}_{\mathrm{a}\mathrm{v}\mathrm{g}}\right)}^{2}}}{{V}_{\mathrm{a}\mathrm{v}\mathrm{g}}}{S}_{i} $$ (7) 式中:$ \gamma $ 为流动均匀性指数,取值范围为 [0,1],随着区间内数值的增大,均匀性越好,其中,1表示理想状态下均匀流动,0表示流体只经过1个测试点,这2个条件均为假设条件,在现实中并不存在;$ {V}_{i} $ 和 $ {V}_{\mathrm{a}\mathrm{v}\mathrm{g}} $ 分别为监测点的速度和监测平面上的平均流速;$ n $ 为截面上的计算单元数;$ {S}_{i} $ 为计算单元的面积;S为计算截面的面积。本研究截取7个截面,距离池底的深度分别为2、4、6、8、10、12、14 m。
图7给出了不同工况下监测面上的流动均匀性指数。如图7-a所示,随着日换水次数的减少,养殖舱内的整体流动均匀性水平逐渐下降,平均下降率约为3.745%,最大下降幅度约为4.7%。但日换水次数的改变并未影响流体均匀性的分布规律,即随着距池底深度的增加,均匀性指数逐渐减小,在12 m处小幅回升,14 m处 (第一排进水口)降到最低。如图7-b所示,随着进水口数量的减少,养殖舱内的整体流动均匀性逐渐下降,平均下降率约为5.58%,最大降幅约为18.4%,最大下降点出现在B4工况的8 m处 (进水口)。这与秦康等[17]的研究结果相似。如图7-c所示,C2工况由于受到进水射流的影响,其均匀性指数明显下降,而C3工况的整体流动均匀性有明显提升。
日换水次数的减少导致养殖舱内水流旋转速度变慢,低流速区域扩散,破坏了流速的分布规律,进而降低了舱内流体的均匀性水平。在进水流量一定的条件下,减少进水口的数量使得进水流速增加,对舱内水体的作用面积减少,高流速区域的集中加大了舱内水体的流速差,导致流速分布紊乱,进而降低了流体的均匀性。与减少日换水次数相比,减少进水口的数量对舱内流体均匀性的影响更为显著。进水射流之间的对冲作用打破了原有的旋转运动,使得流速分布变得无规律,显著降低了舱内流体的均匀性。
2.1.3 不同进水方式下的养殖舱阻力
为评价不同进水方式下养殖舱内的能量利用率,引入养殖池阻力系数 $ {C}_{\mathrm{t}} $[34],阻力系数越小表示养殖池中能量利用率越高,计算公式为:
$$ {C}_{\mathrm{t}}{\mathrm{=}}\frac{2Q({V}_{\mathrm{i}\mathrm{n}}{\text{−}}{V}_{\mathrm{a}\mathrm{v}\mathrm{g}})}{A\times {V}_{\mathrm{a}\mathrm{v}\mathrm{g}}^{2}} $$ (8) 式中:$ Q $ 为进水流量 (kg·s−1);$ A $ 为湿表面积 (池底和被水浸湿的侧壁表面积之和,m2);$ {V}_{\mathrm{i}\mathrm{n}} $ 为进水口流速 (m·s−1);$ {V}_{\mathrm{a}\mathrm{v}\mathrm{g}} $ 为养殖舱内平均流速 (m·s−1)。
图8展示了不同进水方式下养殖舱阻力系数与舱内流速平均值的关系。当进水口数量保持不变 (A类与C类工况),养殖舱内的平均流速与阻力系数呈现负相关性;随着平均流速的降低,阻力系数逐渐增大,这与于林平等[22]的研究结果相似。而当进水口数量减少时,平均流速与阻力系数则表现出正相关性,阻力系数随着平均流速的增加而逐渐增大。
随着日换水次数的减少,单位时间内输入养殖舱的总能量降低,而射流水体与池壁、池底等的摩擦撞击导致的能量损失占比相对增加,使得维持水体质点间相对运动的能量消耗比率下降,进而导致养殖池内水体的整体流速降低,养殖池系统的阻力系数相应升高。但在进水流量一定时,减少进水口数量,会使输入舱内的能量更集中,导致舱内水体的均匀性降低,阻力系数上升。在C2工况下,由于进水口的布设位置不当,进水射流之间的相互影响导致舱内的能量利用率最低。
2.2 进水方式对大黄鱼适养性水平的影响
流速大小是评估流场是否适合鱼类生长繁殖的重要指标,适宜的流速可以增强鱼类的游动能力、改善肌肉张力和提高存活率[35]。适合鱼类生长的流速与其体长有关,对于1龄大黄鱼(体长25 cm左右),其所能适应的最大流速不应超过0.5 m·s−1[36-37],当大黄鱼长期处于此流速或高于此流速的流场中,会加速其触底,影响生长发育。但当流场中的流速过小又会造成残饵、鱼粪等固体颗粒物的堆积,影响水体环境,从而导致鱼群病害或死亡[38]。本研究将1龄大黄鱼的适宜流速设为0.05~0.50 m·s−1,引入适渔流速区域占比 (Proportion of velocity suitable for aquaculture, PVA)的概念[12],其计算公式为:
$$ \mathrm{P}\mathrm{V}\mathrm{A}{\mathrm{=}}\frac{{\sum }_{i{\mathrm{=}}1}^{n}{V}_{i}\left({\varphi }_{\mathrm{m}\mathrm{i}\mathrm{n}}\leqslant {\varphi }_{i}\leqslant {\varphi }_{\mathrm{m}\mathrm{a}\mathrm{x}}\right)}{{\sum }_{i{\mathrm{=}}1}^{n}{V}_{i}} $$ (9) 式中:$ {V}_{i} $ 为网格的体积 (m3);$ {\varphi }_{i} $ 为网格的流速 (m·s−1);$ {\varphi }_{\mathrm{m}\mathrm{a}\mathrm{x}} $ 为适渔最大流速 (0.5 m·s−1);$ {\varphi }_{\mathrm{m}\mathrm{i}\mathrm{n}} $ 为适渔最小流速 (0.05 m·s−1);n为计算网格总数。
图9展示了各个工况的适渔流速区域占比。随着日换水次数的减少,养殖舱内适合大黄鱼生长的区域逐渐缩小,最大降幅为44.13%。而进水口数量的改变对适渔流速区域占比的影响并不显著,均保持在较高水平。对于进水口位置的改变,C2工况内的流场状态明显不适合大黄鱼的养殖,C3工况相较于C1工况变化幅度并不大。
由于压力、光照以及生物学特性等原因导致大黄鱼大部分时间生活在养殖舱的底部,所以底部流场的状态对其生长具有重要影响。若底层水流出现不稳定和不可预测性,鱼类在这种水体环境中难以及时调整游泳姿势,产生应激反应,导致不必要的能量消耗。同时,良好的水体混合性能也可使溶解氧分布均匀,进而避免鱼群因过度聚集产生损伤。为评估养殖舱底部流场速度分布的均匀性,引入速度分布均匀性系数DU50,公式为:
$$ {\mathrm{D}\mathrm{U}}_{50}{\mathrm{=}}\frac{{V}_{50}}{{V}_{\mathrm{a}\mathrm{v}\mathrm{g}}}\times 100{\text{%}} $$ (10) $$ {\mathrm{D}\mathrm{U}}_{50(\mathrm{修}\mathrm{正})}{\mathrm{=}}\frac{{\sum }_{i{\mathrm{=}}1}^{n}{\mathrm{D}\mathrm{U}}_{50}\times {h}_{i}}{{\sum }_{i{\mathrm{=}}1}^{n}{h}_{i}} $$ (11) 式中:$ {h}_{i} $ 为监测平面的水深 (m), $ {h}_{1} $=0.7 m、$ {h}_{2} $=1.4 m、$ {h}_{3} $=2.1 m、$ {h}_{4} $=2.8 m;$ {V}_{50} $ 为监测面速度中前50%的平均值;$ {V}_{\mathrm{a}\mathrm{v}\mathrm{g}} $为监测面速度的平均值;DU50对做深度加权平均处理,记为DU50(修正)。
图10展示了不同进水方式下,养殖舱底部的流速分布均匀性系数,随着日换水次数的减少,流速分布均匀性系数逐渐下降,降幅为2.26%;同样随着进水口数量的减少,流速分布均匀性系数也逐渐下降,降幅为6.24%;当改变进水口位置时,最大降幅为21.91%。相比之下,日换水次数的变动对底层流速分布的均匀性影响较小。
日换水次数的减少导致养殖舱内的整体流速大幅下降,适合大黄鱼生长的流速区域占比也随之减小,但其对底部流速分布的均匀性影响较小。进水口数量的减少造成底部流速分布的均匀性大幅下降,但对大黄鱼生长的流速区域占比影响较小。综合分析可得:适当提升日换水次数,可以扩大养殖舱内大黄鱼的适养范围;同时适当增加进水口的数量,可以改善底部流场环境;在改变进水口的布设位置时,及时避免进水射流之间的对冲影响可有效提升大黄鱼的适养水平。对于30万吨养殖渔船养殖舱,在16次日换水的条件下,将16个进水口均匀地布置在矩形养殖舱的切角处,最有利于大黄鱼的养殖。
3. 结论与展望
本研究通过对不同进水方式的矩形养殖舱水动力特性分析,探究了流场的变化规律,并对1龄大黄鱼的适渔性水平进行评估,结论如下:
1) 改变日换水次数对养殖舱内的流速变化幅度影响较大。随着日换水次数由16次减少到10次,养殖舱内的平均流速下降了46.7%,并且低流速区域的扩散致使养殖舱内适渔流速区域占比下降了44.13%。因此,可通过适当增加日换水次数,提高舱内流速以及大黄鱼的适养性水平。
2) 改变进水口数量对养殖舱内的流动均匀性影响较大。随着进水口的数量由16个减少到4个,舱内的流体均匀性指数下降了18.4%,高流速区域集中,导致养殖舱内流体的流动均匀性明显下降,并且对舱底流场产生不利影响。因此,可通过适当增加进水口数量,提升养殖舱内的流场均匀性,从而改善大黄鱼的生长环境。
3) 改变进水口布设位置时应避免进水射流的相互影响,造成舱内流场处于低流速、低均匀性的水平。
本研究在无鱼条件下对养殖舱内流场进行数值仿真计算,存在局限性。事实上,鱼类的游动可以显著提升水体的混合特性,并且在大黄鱼养殖过程中,温度、光照等因素对其生长、繁殖起着至关重要的作用。因此,针对有鱼群存在的养殖舱进行探究将是下一阶段的工作目标。
-
表 1 舌状蜈蚣藻和其他几种藻类基本营养成分
Table 1 Basic nutrition of G. livida and some other algae
% 样品
sample水分
moisture灰分
ash粗纤维
crude fiber粗脂肪
crude fat粗蛋白
crude protein舌状蜈蚣藻 G. livida 8.94 13.32 1.80 0.18 12.28 带形蜈蚣藻[13] G. turuturu 8.92 16.25 6.60 2.18 22.77 海带[14] S. japonica 1.98~2.33 35.3~40.7 6.3~6.9 0.3~0.7 10.2~10.7 皱紫菜[15] P. crispata 10.48~12.95 6.83~8.80 2.14~2.90 0.74~2.12 32.81~40.38 脆江蓠[15] G. chouae 10.2~10.9 6.83~8.45 2.74~3.35 0.46~0.72 10.12~14.13 羊栖菜[16] S. fusiforme 7.38 15.07 10.16 0.74 6.05 表 2 舌状蜈蚣藻与带形蜈蚣藻氨基酸质量分数比较
Table 2 Comparison of amino acid mass fraction of G. livida and G. turuturu
mg·g–1 氨基酸
amino acid舌状蜈蚣藻
G. livida带形蜈蚣藻[13]
G. turuturu天冬氨酸** Asp 13.07 20.09 苏氨酸* Thr 6.37 9.37 丝氨酸** Ser 6.48 10.25 谷氨酸** Glu 12.19 43.01 脯氨酸** Pro 5.93 3.12 甘氨酸** Gly 6.81 10.12 丙氨酸** Ala 7.58 17.27 缬氨酸* Val 7.69 16.98 蛋氨酸* Met 0.82 10.61 异亮氨酸* Ile 5.93 9.83 亮氨酸* Leu 10.21 14.24 酪氨酸 Tyr 3.18 7.06 苯丙氨酸* Phe 7.69 7.26 赖氨酸* Lys 7.03 8.65 组氨酸 His 1.03 2.27 精氨酸 Arg 6.81 28.00 色氨酸* Trp 1.03 2.55 氨基酸总量 TAA 109.85 222.41 必需氨基酸 EAA 46.77 81.21 非必需氨基酸 NEAA 63.08 131.86 呈味氨基酸 DAA 52.05 103.86 必需氨基酸占总氨基酸的比值/% EAA/TAA 42.58 36.52 必需氨基酸与非必需氨基酸的比值/% EAA/NEAA 74.15 61.59 注:*. 必需氨基酸;**. 呈味氨基酸 Note: *. essential amino acid;**. delicious amino acid 表 3 舌状蜈蚣藻蛋白质中必需氨基酸的组成与比较分析
Table 3 Composition and comparative analysis of essential amino acids in protein of G. livida
必需氨基酸
EAAFAO/WHO推荐模式
FAO/WHO evaluation mode全蛋白模式
egg protein evaluation mode质量分数/mg·g–1
mass fraction氨基酸评分
AAS化学评分
CS苏氨酸 Thr 40 51 47.23 118.07 92.61 缬氨酸 Val 50 73 57.00 114.00 78.08 异亮氨酸 Ile 40 66 43.97 109.93 66.63 亮氨酸 Leu 70 88 75.73 108.19 86.06 苯丙氨酸+酪氨酸 Phe+Tyr 60 100 80.62 134.36 80.62 赖氨酸 Lys 55 64 52.12 94.76 81.43 色氨酸 Trp 10 16 7.65 76.54 47.84 表 4 舌状蜈蚣藻和带形蜈蚣藻的氨基酸比值系数 (RC) 及氨基酸比值系数分 (SRC)
Table 4 Ratio coefficient (RC) and score of ratio coefficient (SRC) of amino acids of G. livida and G. turuturu
必需氨基酸
EAA舌状蜈蚣藻
G. livida带形蜈蚣藻
G. turuturu苏氨酸 Thr 1.09 0.98 缬氨酸 Val 1.06 1.42 异亮氨酸 Ile 1.02 1.03 亮氨酸 Leu 1.00 0.85 苯丙氨酸+酪氨酸 Phe+Tyr 1.24 1.00 赖氨酸 Lys 0.88 0.66 色氨酸 Trp 0.71 1.07 氨基酸比值系数分 SRC 83.09 76.84 注:表中数据为氨基酸比值系数 Note: The data in the table are the ratio coefficients of amino acids. 表 5 舌状蜈蚣藻和带形蜈蚣藻的脂肪酸种类及其质量分数
Table 5 Type and mass fraction of fatty acids of G. livida and G. turuturu
% 脂肪酸
fatty acid舌状蜈蚣藻
G. livida带形蜈蚣藻[24]
G. turuturu十一碳酸 C11∶0 0.56 − 十二碳酸 C12∶0 0.13 − 十四碳酸 C14∶0 6.63 5.76 十五碳酸 C15∶0 0.63 3.24 十六碳酸 C16∶0 36.80 51.97 十七碳酸 C17∶0 0.27 0.95 十八碳酸 C18∶0 4.93 4.49 二十一碳酸 C21∶0 − 7.9 二十二碳酸 C22∶0 − 2.23 二十四碳酸 C24∶0 − 0.74 十六碳一烯酸 C16∶1 2.96 − 十八碳一烯酸 C18∶1 13.63 − 二十碳一烯酸 C20∶1 0.22 − 十七碳三烯酸 C17∶3 0.17 − 亚油酸 C18∶2 3.07 − 二十碳二烯酸 C20∶2 − 1.08 二十碳四烯酸 C20∶4 18.97 − 二十碳五烯酸 (EPA) C20∶5 0.40 − 二十二碳六烯酸 (DHA) C22∶6 1.01 − 饱和脂肪酸总量 ∑SFA 49.95 77.28 单不饱和脂肪酸总量 ∑MUFA 16.81 0 多不饱和脂肪酸总量 ∑PUFA 23.61 1.08 注:−. 未检出 Note: −. not detected 表 6 舌状蜈蚣藻和带形蜈蚣藻的矿物元素质量分数
Table 6 Mass fraction of mineral elements of G. livida and G. turuturu
mg·g–1 元素
element舌状蜈蚣藻
G. livida带形蜈蚣藻[24]
G. turuturu钠 Na 32 594.22 − 镁 Mg 9 267.25 − 铝 Al 60.18 − 钙 Ca 582.00 2 728.63 铬 Cr 1.53 1.31 锰 Mn 11.95 15.38 铁 Fe 98.18 102.69 钴 Co 0.13 0.38 镍 Ni 7.43 1.13 铜 Cu 3.27 2.69 锌 Zn 45.36 98.38 砷 As 15.76 6.25 硒 Se 0.04 1.56 镉 Cd 0.63 0.94 铅 Pb 0.28 2.50 注:−. 未测定 Note: −. not determined -
[1] MARCOS M R, EVA G V, YAISEL J, et al. Molecular barcoding confirms the presence of exotic Asian seaweeds (Pachymeniopsis gargiuli and Grateloupia turuturu) in the Cantabrian Sea, Bay of Biscay[J]. Peerj, 2017, 3116: 2-19.
[2] GARCIAJIMENEZ P, MONTEROFERNANDEZ M, ROBAINA R R. Molecular mechanisms underlying Grateloupia imbricata (Rhodophyta) carposporogenesis induced by methyl jasmonate[J]. J Phycol, 2017, 53(6): 1340-1344.
[3] 夏邦美. 中国海藻志· 第二卷· 红藻门· 第三册[M]. 北京: 科学出版社出版, 2004: 118-122. [4] 陈伟洲, 宋志民, 黄中坚. 温度、光照强度对舌状蜈蚣藻早期发育的影响[J]. 南方水产科学, 2013, 9(6): 14-19. [5] TANG L, CHEN Y, JIANG Z, et al. Purification, partial characterization and bioactivity of sulfated polysaccharides from Grateloupia livida[J]. Int J Biol Macromol, 2017, 94(Pt A): 642-652.
[6] NGO D H, KIM S K. Sulfated polysaccharides as bioactive agents from marine algae[J]. Int J Biol Macromol, 2013, 62: 70-75.
[7] MUNIER M, MORANÇAIS M, DUMAY J, et al. One-step purification of R-phycoerythrin from the red edible seaweed Grateloupia turuturu[J]. J Chromatogr B, 2015, 992: 23-29.
[8] LE G C, DUMAY J, DONNAY M C, et al. Ultrasound-assisted extraction of R-phycoerythrin from Grateloupia turuturu with and without enzyme addition[J]. Algal Res, 2015, 12: 522-528.
[9] 李卫东, 吴茜, 蒋丹, 等. 蜈蚣藻属海藻的应用及人工养殖前景展望[J]. 河北渔业, 2012, 11: 48-54. [10] 蔡秋杏, 吴燕燕, 李来好, 等. 厦门白姑鱼腌制加工过程中的脂肪酸变化分析[J]. 食品科学, 2015, 36(12): 76-81. [11] 王哲平, 刘淇, 曹荣, 等. 野生与养殖刺参营养成分的比较分析[J]. 南方水产科学, 2012, 8(2): 64-70. [12] 朱圣陶, 吴坤. 蛋白质营养价值评价——氨基酸比值系数法[J]. 营养学报, 1988, 10(2): 187-190. [13] 付晓婷, 秦子涵, 许加超, 等. 带形蜈蚣藻营养品质的分析和评价[J]. 营养学报, 2011, 33(2): 199-200. [14] 姚海芹, 王飞久, 刘福利, 等. 食用海带品系营养成分分析与评价[J]. 食品科学, 2016, 37(12): 95-98. [15] 陈伟洲, 蔡少佳, 刘婕, 等. 养殖海藻皱紫菜和脆江蓠的主要营养成分分析[J]. 营养学报, 2013, 35(6): 613-615. [16] 张晓梅, 郭芮, 苏红, 等. 羊栖菜营养成分分析与安全性评价[J]. 食品工业科技, 2018, 39(4): 296-300. [17] 何建瑜, 赵荣涛, 刘慧慧. 舟山海域厚壳贻贝软体部分营养成分分析与评价[J]. 南方水产科学, 2012, 8(4): 37-42. [18] 杨贤庆, 杨丽芝, 黄卉, 等. 南海鸢乌贼墨汁营养成分分析与评价[J]. 南方水产科学, 2015, 11(5): 138-142. [19] 邬宇茜. L-苯丙氨酸分子印迹聚合物(膜)的合成及性能研究[D]. 昆明: 云南大学, 2012: 10-11. [20] 崔海峰. 不同品系羊栖菜形态和营养成分的初步分析[D]. 哈尔滨: 东北林业大学, 2013: 21-22. [21] 张春娥, 张惠, 刘楚怡, 等. 亚油酸的研究进展[J]. 粮油加工, 2010, 5: 18-21. [22] 沈雷, 唐俊, 张云霞, 等. 花生四烯酸油微胶囊配方优化及稳定性研究[J]. 中国油脂, 2016, 41(4): 14-18. [23] 李妍, 王静, 李麒龙, 等. EPA与DHA最新研究进展[J]. 农产品加工 (学刊), 2013, 3: 6-13. [24] 芮雯, 岑颖洲, 伍秋明. 蜈蚣藻和带形蜈蚣藻脂肪酸成分和无机元素含量分析[J]. 广东药学院学报, 2010, 26(1): 48-50. [25] 杨庆, 樊亚鸣, 曾源, 等. 海带中有害元素As的研究现状及展望[J]. 食品研究与开发, 2015, 36(12): 127-131, 138. [26] RODRIGUES D, FREITAS A C, PEREIRA L, et al. Chemical composition of red, brown and green macroalgae from Buarcos bay in Central West Coast of Portugal[J]. Food Chem, 2015, 183: 197-207.
[27] DENIS C, MORANÇAIS M, LI M, et al. Study of the chemical composition of edible red macroalgae Grateloupia turuturu from Brittany (France)[J]. Food Chem, 2010, 119(3): 913-917.
-
期刊类型引用(2)
1. 张凯,张麟,彭凌风,陈鑫,刘合刚,胡志刚. 基于性别差异的少棘巨蜈蚣蛋白组和转录组联合分析. 时珍国医国药. 2024(04): 892-898 . 百度学术
2. 韩财安,李安东,周美玉,廖怀生. 小龙虾幼苗培育关键技术. 江西水产科技. 2022(05): 30-31+34 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 3974
- HTML全文浏览量: 1443
- PDF下载量: 34
- 被引次数: 3