5种鲈形目淡水鱼肌肉脂肪酸及磷脂组成的研究

韩迎雪, 林婉玲, 杨少玲, 李来好, 黄卉, 杨贤庆, 王锦旭, 吴燕燕, 翟红蕾, 郝淑贤

韩迎雪, 林婉玲, 杨少玲, 李来好, 黄卉, 杨贤庆, 王锦旭, 吴燕燕, 翟红蕾, 郝淑贤. 5种鲈形目淡水鱼肌肉脂肪酸及磷脂组成的研究[J]. 南方水产科学, 2019, 15(1): 85-92. DOI: 10.12131/20180184
引用本文: 韩迎雪, 林婉玲, 杨少玲, 李来好, 黄卉, 杨贤庆, 王锦旭, 吴燕燕, 翟红蕾, 郝淑贤. 5种鲈形目淡水鱼肌肉脂肪酸及磷脂组成的研究[J]. 南方水产科学, 2019, 15(1): 85-92. DOI: 10.12131/20180184
HAN Yingxue, LIN Wanling, YANG Shaoling, LI Laihao, HUANG Hui, YANG Xianqing, WANG Jinxu, WU Yanyan, ZHAI Honglei, HAO Shuxian. Composition of fatty acids and phospholipids in five Perciformes freshwater fish muscle[J]. South China Fisheries Science, 2019, 15(1): 85-92. DOI: 10.12131/20180184
Citation: HAN Yingxue, LIN Wanling, YANG Shaoling, LI Laihao, HUANG Hui, YANG Xianqing, WANG Jinxu, WU Yanyan, ZHAI Honglei, HAO Shuxian. Composition of fatty acids and phospholipids in five Perciformes freshwater fish muscle[J]. South China Fisheries Science, 2019, 15(1): 85-92. DOI: 10.12131/20180184

5种鲈形目淡水鱼肌肉脂肪酸及磷脂组成的研究

基金项目: 国家重点研发计划专项 (2016YFD0400201-6);现代农业产业技术体系建设专项资金 (CARS-46);“扬帆计划”引进创新创业团队专项资助 (2015YT02H109)
详细信息
    作者简介:

    韩迎雪(1993—),女,硕士研究生,研究方向为水产品加工。E-mail: 15737935907@163.com

    通讯作者:

    李来好(1963 — ),男,博士,研究员,从事水产品精深加工与质量安全研究。E-mail: laihaoli@163.com

  • 中图分类号: TS 225.2

Composition of fatty acids and phospholipids in five Perciformes freshwater fish muscle

  • 摘要:

    为研究鲈形目淡水鱼在不同科和种间脂肪酸及磷脂组成的差异,该研究以高体革䱨 (Scortum barcoo)、赤鲈(Perca fluviatilis)、梭鲈(Lucioperca lucioperca)、大口黑鲈(Micropterus salmoides)和罗非鱼(Oreochromis mossambicus)为对象,采用气相色谱-质谱联用仪(GC-MS)、高效液相色谱-蒸发光散射检测(HPLC-ELSD)对以上5种鱼肌肉中的脂肪酸与磷脂进行了分析。结果表明,5种淡水鱼肌肉中共检测出21种脂肪酸,共有的脂肪酸有8种,分别为C14:0,C16:0,C18:0,C16:1,C18:1,C18:2,C20:5和C22:6。多不饱和脂肪酸(PUFA)含量丰富,PUFA含量除高体革䱨为26.81%之外,其他均大于35%。5种淡水鱼肌肉中共检出6种磷脂,分别为磷脂酰乙醇胺(PE)、鞘磷脂(SM)、磷脂酰肌醇(PI)、溶血磷脂酰胆碱(LPC)、磷脂酰丝氨酸(PS)、磷脂酰胆碱(PC)。5种鱼肌肉磷脂中PC含量最高,占据磷脂总量的45.84%~62.55%,PE和PI次之,PS和SM含量相对较少,仅占总磷脂的2.19%~5.27%和2.06%~3.65%。

    Abstract:

    In order to study the difference of fatty acid and phospholipid composition in different families and species of Perciformes freshwater fish, we analyzed the fatty acids and phospholipids in muscles of Scortum barcoo, Perca fluviatilis, Lucioperca lucioperca, Micropterus salmoides and Oreochromis mossambicus by gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC-ELSD). Twenty-one kinds of fatty acids had been detected in those fish. Eight kinds of common fatty acids were detected (C14:0, C16:0, C18:0, C16:1, C18:1, C18:2, C20:5 and C22:6). Polyunsaturated fatty acids (PUFA) were rich and the content was higher than 35% except for S. barcoo (26.81%). Six kinds of phosphatidylcholine were detected in those fish, including phosphatidylethanolamine (PE), sphingomyelin (SM), phosphatidylinositol (PI), lysophosphatidylcholine (LPC), phosphatidylcholine (PS) and phosphatidylcholine (PC). The PC contents in those fish muscles phospholipids were the highest, accounting for 45.84%−62.55% of total phospholipid, followed by PE and PI contents. The PS and SM contents were relatively low, accounting for only 2.19%−5.27% and 2.06%−3.65% of total phospholipids.

  • 中国是世界第一水产养殖大国。为解决近岸养殖面积不足、病害频发等问题,海水养殖产业正在积极向深远海拓展[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) 的适养性水平,以期为超大型养殖平台养殖舱结构设计优化提供理论指导。

    本研究将养殖工船养殖舱内的流体视为三维不可压缩流体,该流体的控制方程为连续性方程、纳维-斯托克斯(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} $ 由湍动能及湍流生成项计算。

    本研究基于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  模型试验验证
    Figure  1.  Model testing and validation

    图1-b所示,模拟结果与试验结果的变化规律几乎相同,整体吻合性较好,误差基本保持在10%以内,其中最大误差处是由于进水射流的影响造成,精度满足本次仿真的要求,因此该模型可以用于工船船载舱养的数值仿真计算。

    本研究的养殖工船,船长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。

    图  2  30万吨级养殖工船养殖舱的布置与结构
    Figure  2.  Arrangement and structure of aquaculture tank of 300 000-ton aquaculture vessel

    本研究采用控制变量的方法探究日换水次数、进水口数量和进水口布设位置对舱内流场的影响,试验工况如表1所示。

    表  1  仿真试验分组
    Table  1.  Grouping of simulation tests
    工况组别
    Working condition
    group
    日换水次数
    Number of water change
    per day
    进水口数量
    Number of
    water inlets
    A1 16 16
    A2 14
    A3 12
    A4 10
    B1 16 16
    B2 12
    B3 8
    B3 4
    C1 16 8
    C2
    C3
    下载: 导出CSV 
    | 显示表格

    养殖舱工作时,进水口与出水口同时打开,在保证初始液位高度不变的条件下,完成1次水体更新即为1次换水,A类工况的进水口布设位置和进水方向如图3-a所示。B类工况的进水口布设位置和进水方向与A类工况相同,并均匀的分布在4个切角处。在研究进水口布设位置时,受MCR跑道池的启发,将4个进水口按照不同的进水方向布置在养殖池的内侧形成C2和C3工况,进水口位置和进水方向如图3-b、3-c所示。

    图  3  养殖舱进水口位置和进水方向示意图
    Figure  3.  Diagram of inlet location and water flow direction of aquaculture tank

    定义养殖舱内水面为自由液面,压力设为标准大气压,舱壁和舱底定义为无滑移壁边界,湍流动能和湍流耗散率设置为二阶迎风格式。基于广义最小余量法(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、10046215481613。在距舱底z=1.2 m的平面上选取7个监测点,流速如图4-b所示。由图可见,粗糙网格的仿真结果与精细网格的仿真结果有一定的差距,而0.3 m的中等网格与0.2 m的精细网格在精度上比较接近,综合考虑计算精度和计算时长,本研究的模拟仿真均采用0.3 m的网格尺寸。

    图  4  网格无关性验证
    Figure  4.  Mesh independence verification

    将基于有限体积法输出的速度数据作为集合{φ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)

    使用高斯拟合的方法计算标准差ω和均值μ,以便更加直观地比对各个流场之间的差异。

    计算时长统一设为7 000 s,采用CFD软件监测计算网格内流体的平均动能和平均湍动能,发现6 000 s后计算趋于稳定状态。图5展示了不同深度截面在第7 000秒时的速度分布。通过观察可知:在进水推流的作用下,高流速区域集中在靠近舱壁处,低流速区域则主要分布在舱体中心以及切角处,此种流场特征与Zhang等 [25]、崔铭超等[14]的研究结果基本一致。对于C2工况,在进水射流的相互影响下难以形成稳定的流场,流速分布比较紊乱,水体旋转流动的特性并未显现;C3工况则以底部2对出水口为中心形成了较为稳定的对称流场。通过对比分析图5中A类和B类工况在同一高度的流速分布,发现减少日换水次数或增加进水口数量均会导致中心低流速区域的扩大。

    图  5  各类工况俯视剖面速度分布云图
    注:z表示到舱底的距离。
    Figure  5.  Contour maps of velocity magnitudes of top-view profile under different working conditions
    Note: z represents the distance to the bottom of breeding tank.

    为了更准确地揭示舱内流速的分布规律,对整个舱内水体的流速进行了定量统计分析。流速的概率密度函数和累积分布函数如图6所示。图6-a显示舱内流速与日换水次数之间存在明显的负相关性。随着日换水次数的减少,流速均值μ下降了47.3%,标准差ω下降了37.6%,表明流速分布趋向集中。另一方面,随着进水口数量的减少,流速均值上升了30.7%,标准差上升了66.3% (图6-b)。C1和C3工况的流速分布大致相同,而C2工况较低流速区域的占比更高(图6-c)。

    图  6  养殖舱内水体流速数理统计分布
    Figure  6.  Mathematical and statistical distribution of water velocity in aquaculture tank

    由于养殖舱进水口设置在舱壁切角处,其附近水体受到进水推流的影响较大,因此会在舱壁处产生较大的流速,随着输入养殖舱的能量在抵抗水体质点之间的相对运动产生损耗时,远离舱壁的流速逐渐降低,并且切角处的水流撞击也会有能量损失,进而在切角处产生低速区。随着日换水次数的减少,单位时间内输入养殖舱的能量降低,舱内水流难以维持较高的流速,造成低流速区域逐渐扩大;在进水流量一定的条件下,减少进水口数量,使进水流速大幅提升,舱内整体流速有所提高。但是相对于改变日换水次数,进水口数量的减少对流速的影响并不显著。在改变进水口布设位置时,对于C2工况,进水射流之间的相互影响造成了能量的损耗,整个养殖舱处在低流速水平;对于C3工况,在进水射流的影响下左右两侧舱体形成反方向转动的水流,而整体则呈“∞”型流动,高流速区域出现在两侧舱壁和舱体中心。

    流体均匀性是评价养殖舱内水体环境的重要指标,均匀的水流有利于鱼类的生长。为了量化这一指标,本研究引入流动均匀性指数,其计算公式为[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工况的整体流动均匀性有明显提升。

    图  7  不同工况下流动均匀性指数
    Figure  7.  Flow uniformity index under different working conditions

    日换水次数的减少导致养殖舱内水流旋转速度变慢,低流速区域扩散,破坏了流速的分布规律,进而降低了舱内流体的均匀性水平。在进水流量一定的条件下,减少进水口的数量使得进水流速增加,对舱内水体的作用面积减少,高流速区域的集中加大了舱内水体的流速差,导致流速分布紊乱,进而降低了流体的均匀性。与减少日换水次数相比,减少进水口的数量对舱内流体均匀性的影响更为显著。进水射流之间的对冲作用打破了原有的旋转运动,使得流速分布变得无规律,显著降低了舱内流体的均匀性。

    为评价不同进水方式下养殖舱内的能量利用率,引入养殖池阻力系数 $ {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]的研究结果相似。而当进水口数量减少时,平均流速与阻力系数则表现出正相关性,阻力系数随着平均流速的增加而逐渐增大。

    图  8  平均流速与养殖舱阻力系数
    Figure  8.  Average flow velocity and resistance coefficient Ct

    随着日换水次数的减少,单位时间内输入养殖舱的总能量降低,而射流水体与池壁、池底等的摩擦撞击导致的能量损失占比相对增加,使得维持水体质点间相对运动的能量消耗比率下降,进而导致养殖池内水体的整体流速降低,养殖池系统的阻力系数相应升高。但在进水流量一定时,减少进水口数量,会使输入舱内的能量更集中,导致舱内水体的均匀性降低,阻力系数上升。在C2工况下,由于进水口的布设位置不当,进水射流之间的相互影响导致舱内的能量利用率最低。

    流速大小是评估流场是否适合鱼类生长繁殖的重要指标,适宜的流速可以增强鱼类的游动能力、改善肌肉张力和提高存活率[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工况变化幅度并不大。

    图  9  不同工况下适养流速区域占比
    Figure  9.  Proportion of velocity suitable for aquaculture under different working conditions

    由于压力、光照以及生物学特性等原因导致大黄鱼大部分时间生活在养殖舱的底部,所以底部流场的状态对其生长具有重要影响。若底层水流出现不稳定和不可预测性,鱼类在这种水体环境中难以及时调整游泳姿势,产生应激反应,导致不必要的能量消耗。同时,良好的水体混合性能也可使溶解氧分布均匀,进而避免鱼群因过度聚集产生损伤。为评估养殖舱底部流场速度分布的均匀性,引入速度分布均匀性系数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%。相比之下,日换水次数的变动对底层流速分布的均匀性影响较小。

    图  10  不同工况下流速分布均匀性系数
    Figure  10.  Velocity distribution uniformity coefficient under different working conditions

    日换水次数的减少导致养殖舱内的整体流速大幅下降,适合大黄鱼生长的流速区域占比也随之减小,但其对底部流速分布的均匀性影响较小。进水口数量的减少造成底部流速分布的均匀性大幅下降,但对大黄鱼生长的流速区域占比影响较小。综合分析可得:适当提升日换水次数,可以扩大养殖舱内大黄鱼的适养范围;同时适当增加进水口的数量,可以改善底部流场环境;在改变进水口的布设位置时,及时避免进水射流之间的对冲影响可有效提升大黄鱼的适养水平。对于30万吨养殖渔船养殖舱,在16次日换水的条件下,将16个进水口均匀地布置在矩形养殖舱的切角处,最有利于大黄鱼的养殖。

    本研究通过对不同进水方式的矩形养殖舱水动力特性分析,探究了流场的变化规律,并对1龄大黄鱼的适渔性水平进行评估,结论如下:

    1) 改变日换水次数对养殖舱内的流速变化幅度影响较大。随着日换水次数由16次减少到10次,养殖舱内的平均流速下降了46.7%,并且低流速区域的扩散致使养殖舱内适渔流速区域占比下降了44.13%。因此,可通过适当增加日换水次数,提高舱内流速以及大黄鱼的适养性水平。

    2) 改变进水口数量对养殖舱内的流动均匀性影响较大。随着进水口的数量由16个减少到4个,舱内的流体均匀性指数下降了18.4%,高流速区域集中,导致养殖舱内流体的流动均匀性明显下降,并且对舱底流场产生不利影响。因此,可通过适当增加进水口数量,提升养殖舱内的流场均匀性,从而改善大黄鱼的生长环境。

    3) 改变进水口布设位置时应避免进水射流的相互影响,造成舱内流场处于低流速、低均匀性的水平。

    本研究在无鱼条件下对养殖舱内流场进行数值仿真计算,存在局限性。事实上,鱼类的游动可以显著提升水体的混合特性,并且在大黄鱼养殖过程中,温度、光照等因素对其生长、繁殖起着至关重要的作用。因此,针对有鱼群存在的养殖舱进行探究将是下一阶段的工作目标。

  • 图  1   7种磷脂混合标准品的HPLC-ELSD图谱

    Figure  1.   HPLC-ELSD profile of seven phospholipid standards

    图  2   5种鲈形目淡水鱼肌肉磷脂的HPLC-ELSD图谱

    a. 梭鲈;b. 罗非鱼;c. 赤鲈;d. 高体革䱨;e. 大口黑鲈

    Figure  2.   HPLC-ELSD profile of phospholipid in five Perciformes freshwater fish muscle

    a. L. lucioperca;b. O. mossambicus;c. P. fluviatilis;d. S. barcoo;e. M. salmonides

    表  1   升温程序表

    Table  1   Heating program

    阶段
    stage
    速率/℃·min–1
    rate
    终点温度/℃
    final temperature
    保持时间/min
    retention time
    0 110 4
    1 10 160 1
    2 4 210 3
    3 4 240 8
    下载: 导出CSV

    表  2   梯度洗脱表

    Table  2   Gradient elution

    t/min 正己烷 (含0.04%三乙胺)/%
    n-hexane (contain 0.04% three ethylamine)
    异丙醇/%
    isopropanol
    13%乙酸溶液/%
    13% acetic acid solution
    0 40 57 3
    8 40 50 10
    15 40 50 10
    15.1 40 57 3
    24 40 57 3
    下载: 导出CSV

    表  3   5 种淡水鱼肌肉的脂质含量

    Table  3   Contents of lipid in muscles of five freshwater fish g·kg–1

    鱼种
    fish species
    脂质含量
    lipid content
    大口黑鲈 M. salmonides 88.27±5.42b
    高体革䱨 S. barcoo 143.84±16.14c
    赤鲈 P. fluviatilis 61.27±8.08a
    梭鲈 L. lucioperca 53.11±3.09a
    罗非鱼 O. mossambicus 61.18±0.20a
     注:同列不同小写字母表示差异显著 (P<0.05)
     Note: Different lowercase superscripts in the same column indicate significant difference (P< 0.05).
    下载: 导出CSV

    表  4   5 种淡水鱼肌肉脂肪酸组成及相对含量

    Table  4   Composition and relative contents of fatty acids in muscles of five freshwater fish %

    脂肪酸
    fatty acid
    罗非鱼
    O. mossambicus
    大口黑鲈
    M. salmonides
    高体革䱨
    S. barcoo
    赤鲈
    P. fluviatilis
    梭鲈
    L. lucioperca
    月桂酸 (C12:0) lauric acid 0.19±0.04 0.07±0.03
    肉豆蔻酸 (C14:0) myristic acid 2.60±0.11b 3.18±0.13c 0.70±0.05a 0.83±0.02a 0.53±0.28a
    十五烷酸 (C15:0) pentadecanoic acid 0.78±0.07 0.62±0.08
    棕榈酸 (C16:0) palmitic acid 17.19±0.49a 18.17±0.12a 25.56±1.92b 25.47±0.17b 25.92±4.34b
    十七烷酸 (C17:0) heptadecanoic acid 1.35±0.17c 0.81±0.04b 0.18±0.04a
    硬脂酸 (C18:0) stearic acid 8.32±0.23bc 9.09±0.05c 7.81±0.65b 5.59±0.23a 8.12±0.28b
    总饱和脂肪酸 ∑SFA 30.43±1.12 31.88±0.32 34.06±1.22 32.13±0.41 34.57±4.34
    十四烯酸 (C14:1) tetradecenoic acid 0.11±0.01
    棕榈油酸 (C16:1) Palmitoleic acid 6.73±0.35d 7.66±0.22e 2.33±0.18b 4.01±0.27c 1.14±0.16a
    十七烯酸 (C17:1) heptadecenoic acid 0.33±0.00 0.76±0.13
    油酸 (C18:1) oleic acid 25.45±0.65c 17.42±0.04a 36.80±1.80d 23.72±0.27bc 21.74±2.07b
    二十烯酸 (C20:1) eicosenoic acid 1.61±0.05 1.46±0.03
    总单不饱和脂肪酸 ∑MUFA 34.22±1.05c 27.3±0.43b 39.13±1.99d 27.73±0.54b 22.87±1.91a
    亚油酸 (C18:2) linoleic acid 13.67±0.42c 2.51±0.00a 11.94±0.54b 10.94±0.36b 19.77±0.60d
    亚麻酸 (C18:3) calendic acid 0.36±0.03 0.53±0.35 0.40±0.07 0.60±0.54
    二十碳二烯酸 (C20:2) eicosadienoic acid 0.97±0.07
    二十碳三烯酸 (C20:3) eicosatrienoic acid 1.44±0.08 2.11±0.57
    二十碳四烯酸 (C20:4) eicosatetraenoic acid 4.58±3.59 6.50±0.07 1.96±0.42 1.37±0.04
    二十碳五烯酸 (C20:5) eicosapentaenoic acid (EPA) 5.29±2.46bc 7.27±0.05c 2.00±0.23a 2.80±0.10ab 2.05±0.25a
    二十一碳四烯酸 (C21:4) heneicosatetraenoic acid 3.20±1.95 1.45±0.02
    二十一碳五烯酸 (C21:5) heneicosapentaenoic acid 0.83±0.02
    二十二碳六烯酸 (C22:6) docosahexenoic acid (DHA) 6.19±0.33a 22.73±0.62b 8.28±2.63a 23.81±1.07b 20.15±2.61b
    总多不饱和脂肪酸 ∑PUFA 35.35±2.17b 40.82±0.74bc 26.81±3.20a 40.14±0.94bc 42.56±2.43c
     注:−. 表示未检测到;同行不同小写字母表示差异显著 (P<0.05);后表同此
     Note: −. not detected; different lowercase superscripts in the same row indicate significant difference (P<0.05); the same case in the following tables.
    下载: 导出CSV

    表  5   磷脂标准曲线方程及 R2

    Table  5   Equations and R2 values of calibration curves for phospholipid standards

    磷脂种类
    phospholipid species
    方程
    equation
    相关系数
    R2
    磷脂酰乙醇胺 PE y=59 679.1x−5 828.79 0.997
    溶血磷脂酰乙醇胺 LPE y=283 030x+4 490.18 0.991
    磷脂酰肌醇 PI y=45 988.6x−923.945 0.996
    磷脂酰丝氨酸 PS y=174 250x+2 475.06 0.994
    磷脂酰胆碱 PC y=422 685x+1 779.12 0.994
    鞘磷脂 SM y=246 526x−3 738.84 0.997
    溶血磷脂酰胆碱 LPC y=561 390x−9 140.99 0.998
    下载: 导出CSV

    表  6   5 种淡水鱼磷脂组成及相对含量

    Table  6   Composition and relative content of phospholipid in five freshwater fish %

    磷脂
    phospholipid
    大口黑鲈
    M. salmonides
    梭鲈
    L. lucioperca
    罗非鱼
    O. mossambicus
    赤鲈
    P. fluviatilis
    高体革䱨
    S. barcoo
    磷脂酰乙醇胺 PE 20.60±3.08a 28.39±1.55b 21.40±0.81a 30.27±0.80b 28.23±0.19b
    磷脂酰肌醇 PI 11.43±2.13ab 9.16±0.91a 11.56±0.35ab 13.41±0.11b 17.39±0.37c
    磷脂酰丝氨酸 PS 3.26±0.65b 3.45±0.21b 5.27±0.06c 2.19±0.15a 4.91±0.03c
    磷脂酰胆碱 PC 62.55±6.34c 56.94±3.18bc 57.03±0.50bc 49.71±2.01ab 45.84±0.78a
    鞘磷脂 SM 2.17±0.48a 2.06±0.52a 2.12±0.18a 2.85±0.78ab 3.65±0.58b
    溶血磷脂酰胆碱 LPC 2.63±0.90 1.59±0.16
    下载: 导出CSV
  • [1] 王立巧. 天然磷脂在食品及药品加工等方面的应用[J]. 河北化工, 2011, 34(2): 24-25, 45. doi: 10.3969/j.issn.1003-5059.2011.02.010
    [2]

    HAN X. Neurolipidomics: challenges and developments[J]. Front Biosci, 2007, 12(4): 2601-2615.

    [3]

    SUD M, FAHY E, COTTER D, et al. LMSD: LIPID MAPS structure database[J]. Nucleic Acids Res, 2007, 35(Database issue): D527-D532.

    [4] 郇延军, 周光宏, 徐幸莲. 脂类物质在火腿风味形成中的作用[J]. 食品科学, 2004, 25(1): 186-190. doi: 10.3321/j.issn:1002-6630.2004.01.050
    [5] 徐为民, 王道营, 诸永志, 等. 肌内磷脂在腌腊肉制品风味形成中的作用机制研究进展[J]. 江苏农业学报, 2010, 26(1): 204-207. doi: 10.3969/j.issn.1000-4440.2010.01.038
    [6]

    KIDD P M. (PhosphatidylSerine), vital lipid for memory, mood and stress[J]. Total Health, 2008, 30(2): 28.

    [7]

    BREWER J W. Phospholipids: “Greasing the wheels” of humoral immunity[J]. Biochimica et Biophysica Acta-Molecular and Cell Biology of Lipids, 2013, 1831(3, SI): 642-651. doi: 10.1016/j.bbalip.2012.09.018

    [8] 陈文娟. 大黄鱼鱼卵磷脂的提取、鉴定及抗氧化活性研究[D]. 福州: 福建农林大学, 2010: 40-44.
    [9] 陆丽丽, 陈舜胜. 鲤、鲫卵的脂质和脂肪酸特性[J]. 食品科学, 2007(12): 100-104.
    [10] 程新伟, 梁鹏, 涂晓玲, 等. 养殖大黄鱼各部位磷脂组分及其脂肪酸组成分析[J]. 食品工业科技, 2017, 38(4): 49-53.
    [11]

    LIU Z Y, ZHOU D Y, WU Z X, et al. Extraction and detailed characterization of phospholipid-enriched oils from six species of edible clams[J]. Food Chem, 2018, 239: 1175-1181. doi: 10.1016/j.foodchem.2017.07.035

    [12] 莫建光, 卢安根, 徐慧, 等. 食品中磷脂酰胆碱检测技术的研究进展[J]. 食品工业科技, 2010(1): 409-411.
    [13]

    FOLCH J, LEES M, SLOANE STANLEY G H. A simple method for the isolation and purification of total lipides from animal tissues[J]. J Biol Chem, 1957, 226(1): 497-509.

    [14] 蔡秋杏, 吴燕燕, 李来好, 等. 厦门白姑鱼腌制加工过程中的脂肪酸变化分析[J]. 食品科学, 2015, 36(12): 76-81. doi: 10.7506/spkx1002-6630-201512014
    [15] 曹松敏, 吴燕燕, 李来好, 等. 蓝圆鲹传统腌干过程中内源脂肪酶和脂质降解氧化的变化分析[J]. 食品科学, 2017, 38(7): 36-42.
    [16] 吴燕燕, 刘法佳, 李来好, 等. GC-MS检测咸鱼中N-亚硝胺的条件优化[J]. 南方水产科学, 2012, 8(4): 16-22. doi: 10.3969/j.issn.2095-0780.2012.04.003
    [17] 刘邦辉, 王广军, 郁二蒙, 等. 投喂蚕豆和普通配合饲料草鱼肌肉营养成分比较分析及营养评价[J]. 南方水产科学, 2011, 7(6): 58-65. doi: 10.3969/j.issn.2095-0780.2011.06.010
    [18] 丁丽丽, 李来好, 吴燕燕. 咸鱼风味的研究进展[C]//“食品加工与安全”学术研讨会暨2010年广东省食品学会年会论文集. 湛江, 2010: 148-151.
    [19] 徐迪迪. 蛋黄磷脂不同方法提取与检测分析的研究[D]. 武汉: 武汉工程大学, 2010: 15-20.
    [20] 卢航, 里慧, 赵景华, 等. 利用高效液相色谱-蒸发光检测器法分析鮰鱼脑中磷脂组成[J]. 食品科学, 2015, 36(24): 177-180. doi: 10.7506/spkx1002-6630-201524032
    [21] 朱瑶, 陈慧民, 卢航, 等. 刺参性腺中6种磷脂含量的HPLC-ELSD分析[J]. 大连海洋大学学报, 2015, 30(4): 422-425.
    [22]

    GRAEVE M, JANSSEN D. Improved separation and quantification of neutral and polar lipid classes by HPLC-ELSD using a monolithic silica phase: application to exceptional marine lipids[J]. J Chromatogr B, 2009, 877(20): 1815-1819.

    [23]

    KIEŁBOWICZ G, MICEK P, WAWRZEŃCZYK C. A new liquid chromatography method with charge aerosol detector (CAD) for the determination of phospholipid classes[J]. Talanta, 2013, 105(105): 28-33.

    [24] 郭辉, 胡晨, 钱俊青, 等. 鸦胆子油HPLC-ELSD指纹图谱研究[J]. 中成药, 2012, 34(7): 1311-1315. doi: 10.3969/j.issn.1001-1528.2012.07.027
    [25] 陈文娟, 陈建福. 大黄鱼鱼卵脂质含量及脂肪酸组成分析[J]. 肉类工业, 2012(4): 21-23. doi: 10.3969/j.issn.1008-5467.2012.04.010
    [26] 刁全平, 郭华, 吕琳琳, 等. 银鲳鱼脂肪酸的气相色谱-质谱分析[J]. 鞍山师范学院学报, 2015, 32(6): 51-53. doi: 10.3969/j.issn.1008-2441.2015.06.013
    [27] 宋长虹, 唐生, 郝克非, 等. 中国居民日常食物中脂肪酸含量的分析[J]. 食品与机械, 2014, 30(05): 61-63.
    [28] 张雪琰, 牟志春, 高建国, 等. 4种海水鱼肉中脂肪酸组成分析及营养评价[J]. 食品研究与开发, 2013(23): 111-113. doi: 10.3969/j.issn.1005-6521.2013.23.030
    [29]

    SWANSON D, BLOCK R, MOUSA S A. Omega-3 fatty acids EPA and DHA: health benefits throughout life[J]. Adv Nutr, 2012, 3(1): 1-7. doi: 10.3945/an.111.000893

    [30] 王琦, 薛长湖, 徐杰, 等. HPLC-ELSD法测定几种水产品卵及生殖腺中磷脂的组成[J]. 分析仪器, 2012(5): 18-22. doi: 10.3969/j.issn.1001-232X.2012.05.005
    [31]

    YAGUCHI T, NAGATA T, NISHIZAKI T. Dilinoleoylphosphatidylcholine ameliorates scopolamine-induced impairment of spatial learning and memory by targeting alpha 7 nicotinic ACh receptors[J]. Life Sci, 2009, 84(9/10): 263-266.

    [32]

    SU Y F, MEI L J, LIU B, et al. Effect of milk powder adding phosphatidylcholine from soybean lecithin on the memory of rats[J]. China Dairy Ind, 2010, 38(1): 29-30.

    [33]

    TEAGUE J, SOUBIAS O, PETRACHE H A, et al. Elastic properties of polyunsaturated phosphatidylethanolamines influence rhodopsin function[J]. Faraday Discuss, 2013, 161: 383-395.

    [34] 马方, 杨宜婷, 陈则华. 不同类型n-3多不饱和脂肪酸对心血管疾病的防治作用及其机制研究进展[J]. 中国油脂, 2018, 43(2): 65-69. doi: 10.3969/j.issn.1003-7969.2018.02.015
    [35]

    WAT E, TANDY S, KAPERA E, et al. Dietary phospholipid-rich dairy milk extract reduces hepatomegaly, hepatic steatosis and hyperlipidemia in mice fed a high-fat diet[J]. Atherosclerosis, 2009, 205(1): 144-150. doi: 10.1016/j.atherosclerosis.2008.12.004

    [36] 童铃, 金毅, 徐坤华, 等. 3种鲣鱼背部肌肉的营养成分分析及评价[J]. 南方水产科学, 2014, 10(5): 51-59. doi: 10.3969/j.issn.2095-0780.2014.05.008
  • 期刊类型引用(2)

    1. 张凯,张麟,彭凌风,陈鑫,刘合刚,胡志刚. 基于性别差异的少棘巨蜈蚣蛋白组和转录组联合分析. 时珍国医国药. 2024(04): 892-898 . 百度学术
    2. 韩财安,李安东,周美玉,廖怀生. 小龙虾幼苗培育关键技术. 江西水产科技. 2022(05): 30-31+34 . 百度学术

    其他类型引用(1)

图(2)  /  表(6)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 3
出版历程
  • 收稿日期:  2018-08-15
  • 修回日期:  2018-09-19
  • 录用日期:  2018-10-11
  • 网络出版日期:  2018-12-04
  • 刊出日期:  2019-02-04

目录

/

返回文章
返回