基于空气(普通包装)和真空包装条件下鲈鱼片脂肪氧化指标建立鲈鱼货架期预测模型

鞠 健,胡佳慧,熊光权,乔 宇,*,廖 李,汪 兰,吴文锦,丁安子,石 柳

(1.湖北省农业科学院农产品加工与核农技术研究所,湖北省农业科技创新中心农产品加工分中心,湖北武汉 430064; 2.湖北工业大学生物工程与食品学院,湖北武汉 430068)

鲈鱼(Lateolabraxjaponicus)又名花鲈，其肉质鲜嫩,营养价值丰富,体内富含蛋白质、VA、VB、钙、镁、锌、硒等营养元素,具有补肝肾、益脾胃的功效,深受广大消费者喜爱[1]。特别是鲈鱼中所含的二十二碳六烯酸(DHA)和二十碳五烯酸(EPA),是脑黄金和深海鱼油的主要构成元素,对降低血胆固醇的含量和预防血脂异常等具有一定的作用,十分有益于人体健康[2]。然而,与高含量的不饱和脂肪酸相比,鲈鱼中具有抗氧化活力的维生素E含量较低,因此鲈鱼在贮藏过程中很容易与空气中的氧发生自由基连锁反应,加之鲈鱼体内自身所携带的酶类及在贮藏过程中微生物生长繁殖产生的酶等都会加剧鲈鱼脂肪的氧化,导致一些酸类,醛类和胺类物质的产生,促发腐败,致使鲈鱼因肌肉组织结构和风味劣化而失去营养价值[3]。

水产品的品质在贮藏期间会随着时间的延长而不断变化,如何通过一种有效快速的方法来评估水产品在贮藏期间品质和鲜度的变化并能够较为准确的预测其剩余货架期显得尤为重要。大量研究表明,食品在贮藏期间品质的变化可以通过化学动力学模型得到很好的反应[47]。周果等[8]以挥发性盐基氮(TVBN值)、鱼品鲜度(K值)及菌落总数(TVC)为指标,研究了梭子蟹在冷藏、冰藏和冰温贮藏期间的变化,通过零级化学动力学模型与Arrhenius方程建立了适合梭子蟹的货架期预测模型,结果表明该预测模型具有较高的拟合度,预测值与实测值的相对误差在10%以内;郭全友等[9]根据大黄鱼在不同贮藏温度下的腐败速率与温度的相关性开发出了3种货架期预测模型(Exponential、Schoolfield和Squareroot),经过对模型的验证后表明其误差分别为0%～13.8%、0.9%～9.8%和0.2%～22.1%。但目前关于货架期预测模型的相关研究均是在普通包装即空气条件下进行的,而对真空包装条件下水产品品质变化动力学和相关预测模型的研究还尚未见报道,且有关鲈鱼货架期预测模型的研究尚属空白。因此,本文以鲈鱼为原料,研究其在微冻(2 ℃)、冰温(0 ℃)和冷藏(4 ℃)过程中AV,POV和TBA值的变化规律,并运用化学动力学模型建立基于鲈鱼脂肪氧化指标变化的货架期模型,用以预测鲈鱼鱼片在贮藏期间的脂肪氧化情况及剩余货架期,为提高鲈鱼片在贮藏过程中的品质提供理论依据。

1 材料与方法

1.1 材料与仪器

鲈鱼 规格(800±50) g,活体,购于武商量贩农科城店;尼龙包装袋(PA) 市售;硫酸、盐酸、氯化钠、乙醇、1,1,3,3四乙氧基丙烷、2硫代巴比妥酸、无水碳酸钠、碳酸氢钠、酒石酸钾钠、氯化镁、高氯酸、氢氧化钠 分析纯(AR),国药集团化学试剂有限公司;溴甲酚绿、甲基红、次甲基蓝 连云港市鑫源化工股份有限公司;平板计数营养琼脂 北京奥博星生物技术有限责任公司;高氯酸 北京双旋微生物培养制品厂;硫代巴比妥酸、没食子酸 天津市光复精细化工研究所;正辛醇 天津市大茂化学试剂厂;茶多酚(纯度95%) 上海梦荷生物科技有限公司;牛血清白蛋白 杭州吴天生物技术服务有限公司。

UV3802分光光度计 上海尤尼科仪器有限公司;BS210 型电子天平 德国 Sartorius Instruments 有限公司;CR00 色差计 柯尼卡美能达株式会社;DGX9143B 电热恒温鼓风干燥箱 上海雷磁仪器生产厂;SPX250BZ 生化培养箱 上海博讯实业有限公司医疗设备厂;SWCJ净化工作台 苏州净化设备有限公司;FG2便携式 pH计 梅特勒托利多仪器(上海)有限公司;101B2型均质机 上海一恒科技有限公司;KjeltecTM 8400型凯氏定氮仪 德国FOSS公司;HH6数显恒温水浴锅 国华电器有限公司;真空包装机 诸城市昊坤装备科技有限公司;超纯水系统 上海博饮环保科技有限公司;手压式薄膜封口机 永康市喜尚工有限公司。

1.2 实验方法

1.2.1 原料预处理 将购回的鲈鱼宰杀后用无菌水洗净,去皮,分别取背部和腹部鱼肉切成5 cm×5 cm×0.5 cm的鱼片。在无菌操作台上将切好的鲈鱼片装入已经灭菌的保鲜袋中(灭菌条件121 ℃,15 min),将鱼片样本随机分成2个处理组:a. 空气包装组,即将鱼片直接采用封口机进行封口包装;b. 真空包装组,即采用真空包装机进行抽真空处理并分别置于2、0和4 ℃冰箱中贮藏,每隔3 d 取样进行检测,每份样品准备三份,用作平行实验。

1.2.2 氧化指标的测定

1.2.2.1 酸价(AV)和过氧化值(POV)的测定 称取5 g 1.2.1中处理的样品,分别按照国标GB/T5009.372003中酸价滴定的方法和过氧化值滴定的方法进行测定。其中单位用(mg KOH·g-1)来表示，过氧化值单位为meq·kg-1[10]。

1.2.2.2 TBA值的测定 参照Salih等[11]的方法。分别取背部和腹部鱼肉置于含有25 mL超纯水的烧杯中混均,加入25 mL的5%三氯乙酸(TCA)溶液,混匀,室温下静置30 min,过滤;用5% TCA定容至50 mL,将上清液与TBA溶液(0.02 mol/L)按1∶1混匀。置于80 ℃的恒温水浴中加热40 min,冷却,在532 nm下测定吸光度。TBA值用丙二醛(MDA)的质量分数表示,单位为mg MDA/kg。

式中:C为丙二醛质量,mg;m为样品质量,g。

1.2.3 鲈鱼感官鉴定 感官评定主要是依据鱼肉的色泽、气味、组织形态和组织弹性指标进行评定[13]。评定人员由6名经过培训的评价员组成,定期对鲈鱼片进行评定。具体评定标准如表1所示。实验采用加权评分法,各指标权重设置为:色泽20%、气味30%、组织形态30%、组织弹性20%。各特性的平均分乘以其权重即为该特性分值,各特性分值之和为感官评定分。

表1 鲈鱼感官评分标准Table 1 Standard for sensory evaluation of weevers

1.3 数据处理

采用EXCEL 2007软件和SPSS软件进行数据处理和方差分析,Origin 8.5软件作图,并选择合适的货架期模型进行非线性拟合。

2 结果与分析

2.1 不同贮藏温度下鲈鱼背、腹部鱼片脂肪氧化指标的变化

2.1.1 不同温度下AV值的变化趋势 AV值是衡量油脂变质程度的重要指标,也是国际通用的衡量食品腐败程度的一项重要指标[1415]。如图1所示,鲈鱼背部鱼片的初始AV值为1.98 mg KOH·g-1,腹部鱼片的初始AV值为1.46 mg KOH·g-1,三个温度下,鲈鱼背腹部鱼片在整个贮藏期间的AV值均呈显著上升趋势(p<0.05)。在2 ℃贮藏条件下,空气和真空包装处理组中的鲈鱼背腹部鱼片的AV值在贮藏9 d之前没有显著性的差异(p>0.05),从第9 d后空气包装组中的AV值显著高于真空包装组(p<0.05),第12 d时空气背部、空气腹部、真空背部和真空腹部鱼片的AV值已分别达到了4.56、4.77、3.93和4.07 mg KOH·g-1。到贮藏末期第21 d时鲈鱼背、腹部鱼片的AV值分别达到了空气背部5.61 mg KOH·g-1、空气腹部5.96 mg KOH·g-1、真空背部5.08 mg KOH·g-1和真空腹部5.19 mg KOH·g-1,与初始值相比分别提高了255.5%、308.2%、156.6%、183.3%。0和4 ℃贮藏条件下的鲈鱼背、腹部鱼片AV值的变化趋势与2 ℃时的相似,但在整个贮藏期间空气包装组与真空包装组存在显著差异(p<0.05)。在0 ℃贮藏到第18 d时,,空气组鲈鱼背、腹部鱼片的AV值分别为5.47和5.88 mg KOH·g-1,真空组为4.66 和5.44 mg KOH·g-1;在4 ℃贮藏到第12 d时,空气组鲈鱼背、腹部鱼片的AV值分别为5.03和5.24 mg KOH·g-1,真空组为4.42 和4.64 mg KOH·g-1。由此可见,在相同的包装方式,相同贮藏时间下,贮藏温度越高鱼片AV值上升也越快。这说明鱼片脂肪氧化速率随温度升高而加快。然而,在整个贮藏期间真空包装组中鲈鱼背腹部鱼片的AV值始终低于空气包装组。这可能是由于真空处理致使鱼肉中脂肪氧化酶的活力下降,在一定程度上也抑制了一些好氧腐败微生物的生长繁殖。

2.1.2 不同温度下POV值的变化趋势 POV值是衡量油脂氧化酸败程度的又一项重要指标。通常情况下,POV值越低,食品的酸败程度就越低[1617]。如图2所示,鲈鱼背腹部鱼片在2 和0 ℃贮藏期间POV值的变化趋势与AV值的变化趋势相似。4 ℃ 贮藏温度下鲈鱼背、腹部鱼片POV值增加最快,且在贮藏末期9～12 d内增加最为显著,腹部鱼片到第12 d时空气组和真空组腹部鱼片已由初始值0.89 meq·kg-1分别增加到5.93和4.53 meq·kg-1,增加了566.3%和408.9%。而背部鱼片初始值为1.55 meq·kg-1,到第12 d时空气组和真空组分别增加到了5.37和4.10 meq·kg-1,分别增加了246.5%和164.5%。2 ℃温度下贮藏鲈鱼背、腹部鱼片POV值的变化整体也呈现上升趋势,到第21 d时空气组背、腹部鱼片分别为4.68、5.09 meq·kg-1、真空组背、腹部鱼片分别为4.49 和4.18 meq·kg-1。0 ℃ 贮藏温度下鲈鱼背、腹部鱼片POV值的升高趋势则介于其他两个温度之间。在整个贮藏期间鲈鱼背、腹部鱼片POV值不断升高可能是因为鱼片中所含有的不饱和脂肪酸在贮藏期间受到鱼体自身所携带的酶类以及外在环境因素如空气、光照等影响发生了酶促氧化和光氧化等一系列氧化反应,从而产生了氢过氧化物。最终这些物质又经过反复的相互分解,相互作用,继而产生醛类、酮类、酸类、羟甲基和环氧化物等低分子物质,致使鲈鱼片发生氧化变质。

图2 2,0和4 ℃ 贮藏温度下鲈鱼片POV值的变化Fig.2 POV value of weever fillets stored at 2,0 and 4 ℃

2.1.3 不同温度下TBA值的变化趋势 TBA值也被广泛的用于评价鱼肉类制品脂肪的氧化程度,脂肪酸发生氧化的降解产物丙二醛能够与TBA 反应生成稳定的红色复合物[1820]。由图3可知,在整个贮藏期间TBA值与AV和POV值的变化趋势均相似,与贮藏时间呈正相关,这表明TBA、AV和POV值均能够较好地表征鲈鱼片在贮藏期间品质的劣化。在4 ℃ 贮藏温度下鲈鱼背、腹部鱼片的TBA含量在前9 d时增加缓慢,随后出现显著(p<0.05)增加,到贮藏末期第12 d时,分别达到了空气背部1.89 mg MDA ·kg-1、空气腹部2.01 mg MDA ·kg-1、真空背部1.26 mg MDA ·kg-1和真空腹部1.87 mg MDA ·kg-1,比初始值分别增长了687.5%、528.1%、425%和484.4%。在2 ℃下贮藏的鲈鱼背、腹部鱼片的TBA值变化相对较为缓慢,到第21 d时分别为1.42、1.61、1.18和1.34 mg MDA ·kg-1。通常情况下鱼肉中的TBA含量超过1 mg MDA ·kg-1时会开始出现令人难以接受的腥臭味。在4 ℃条件下贮藏的鲈鱼背、腹部鱼片到贮藏第12 d时已经产生了令人难以接受的腥臭味;在0 ℃贮藏鱼片氧化程度比在4 ℃条件下贮藏的鱼片氧化稍缓,但空气组的鱼片第15 d时TBA含量也约为1 mg MDA ·kg-1。因此,贮藏温度越高,油脂氧化速率越快。另外,李婷婷等[21]已用TBA值来预测大西洋鲑的货架期模型并且得到了比较可靠的预测结果。因此采用TBA 值建立的关于鲈鱼背、腹部鱼片的货架期预测模型具有可行性。

图3 2,0和4 ℃贮藏温度下鲈鱼片TBA值的变化Fig.3 TBA value of weever fillets stored at 2,0 and 4 ℃

2.1.4 不同温度下感官评价 鲈鱼背、腹部鱼片的感官评分值在整个贮藏期间不断下降且温度越高下降越快(图4)。在2 ℃下在第15 d时空气包装中背部和腹部鱼片的感官评分值分别为5.4和5.3,而真空包装中背、腹部鱼片在第18 d时感官评分值为5.6和5.5。在0 ℃ 温度下贮藏的普通包装和真空包装中背、腹部鱼片的感官评分值在第15 d时分别为4.9、4.6、5.2和5.1。在4 ℃下空气包装组背部和腹部鱼片的感官评分值在第9 d时分别为4.4、4.0,真空包装组背部和腹部的感官评分值分别为5.3和4.9。从三个不同贮藏温度下的鲈鱼背、腹部鱼片的感官评分值来看,鲈鱼背部鱼片的感官评分值高于腹部鱼片,这与其它三个脂肪氧化指标的变化规律相一致,这表明鲈鱼在贮藏期间其腹部鱼片的脂肪氧化速率快于背部鱼片。

图4 2,0和4 ℃贮藏温度下鲈鱼片感官评价Fig.4 Sensory evaluation of salmon fillets stored at 2,0 and 4 ℃

2.2 货架期预测模型的建立

由一级化学反应动力学方程式(1)改写为(2)如下:

式(1)

A=A0×ekt

式(2)

式中,A、A0、t、k和n分别代表品质因子数、初始值、时间、反应速率和级数。

据相关文献介绍,由微生物的生长繁殖及食品氧化所引起的腐败变质一般遵循一级化学反应动力学[2223],而鲈鱼肉在贮藏期间的腐败变质主要是由微生物和脂肪氧化共同作用导致。且李婷婷等[21]在做大西洋鲑在冷藏期间脂肪氧化的货架期预测模型时也证明了大西洋鲑在冷藏期间脂肪氧化指标的变化符合一级化学反应动力学模型,因此取n=1。然后用一级化学反应动力学模型对以上三个脂肪氧化指标进行拟合,得k值(表2)。

鲈鱼在不同包装方式下的决定系数R2如表2所述。R2决定了脂肪氧化指标与一级化学反应动力学方程的拟合程度,其值越大拟合精度越高。在2 ℃(271K)、0 ℃(273K)和4 ℃(277K)贮藏条件下各包装方式中的鲈鱼片R2≥0.82。其中TBA和POV值的拟合程度较高R2>0.90。由此可见,该动力学方程能够较好地反映鲈鱼在贮藏期间的氧化。

影响微生物生长繁殖的最重要的因素是温度[24]。其中k是关于温度的常数,通常情况下采用Arrhenius方程来表示k与温度的变化关系:

式(3)

式中,k0:频率因子;EA:活化能(J/mol);T:温度(K);R:常数,8.3144 J/(mol·K);k0和 EA是与相关实验物质本性有关的经验常数。

2.3 空气包装下的鲈鱼脂肪氧化指标的非线性拟合

表2中,鲈鱼背、腹部各指标的决定系数R2均>0.820,其中TBA值的整体决定系数R2>0.900,POV值的拟合精度较高。由此可见一级化学反应动力学方程能够较好地描述冷藏期间鲈鱼背、腹部鱼片的脂肪氧化。

表2 空气包装中的非线性拟合参数Table 2 Nonlinear fitting parameters of weever on fat oxidation indicators

k0及EA值的确定:用Arrhenius方程将各氧化指标的k值在不同贮藏温度下拟合即得k0值和EA值:

如表3所示,背、腹部鱼片各脂肪氧化指标的R2均>0.900。R2越大表明相关性越高。相关性越高则越能较好地反映在不同的贮藏温度下鲈鱼背、腹部鱼片反应速率常数随贮藏温度变化的关系。由上表还得到关于鲈鱼背部AV和TBA值的活化能小于腹部,POV值的活化能则大于腹部。这可能是因为AV值和TBA值在背部中的变化更敏感,反之POV值在腹部中更加敏感。

表3 k值对不同贮藏温度的拟合结果Table 3 The fitting results of reaction rate constant k for temperature

由图5可见三个脂肪氧化指标的k值整体与温度呈正相关,这表明它们均可很好地与拟合曲线相拟合。

图5 三种氧化值指标的k值与Arrhenius方程拟合Fig.5 Oxidation indicators k values fitting with Arrhenius equation

由(2)、(3)公式推导得出:

F(A)=k×t=k0×e-E/RT×t=lnAlnA0

式(4)

经上述公式推导可得货架期预测模型公式:

式(5)

式中：A0:初始脂肪氧化指标的数值;A:贮藏 t时间后的脂肪氧化指标的数值;EA:活化能(J/mol);T:绝对温度(K);R:气体常数。

根据鲈鱼背、腹部鱼片在3个贮藏温度下所得到的k0值与EA值带入式(5),可以计算出空气包装中鲈鱼背、腹部鱼片在贮藏过程中AV值、POV值和TBA值的货架期预测模型:

关于AV值的货架期预测模型:

关于POV值的货架期预测模型:

关于TBA值的货架期预测模型:

其中,A0:初始腹部TBA值,A:储藏 t时间后的腹部TBA值。

鲈鱼背、腹部鱼片在贮藏温度下货架期预测模型的验证与评价。贮藏温度下鲈鱼背、腹部鱼片的感官评价以5分作为可接受的最低限值。当冷藏鲈鱼背、腹部鱼片TBA含量超过1 mg MDA ·kg-1时,用所建立的鲈鱼货架期预测模型的预测值与实验值对货架期模型进行验证(表4)。

表4 空气包装条件下鲈鱼货架期模型验证结果Table 4 The shelflife model validation results of weever on fillets in the air package

由表4可以看出,在三个不同的贮藏温度下鲈鱼片关于三个脂肪氧化指标的相对误差均在±10.40%内。这表明预测值与实验值的结果误差较小,在允许的范围内。其中鲈鱼背部鱼片TBA的预测值与实验值的平均相对误差较小为0.83%,其次是背部鱼片POV值平均相对误差为1.09%;而腹部鱼片AV值的预测值与实验值平均相对误差最大为7.95%。这表明关于背部鱼片的TBA值的货架期模型更为准确。

2.4 真空包装下的鲈鱼脂肪氧化指标的非线性拟合

表5中,鲈鱼背、腹部各脂肪氧化指标的决定系数R2均大于0.820,其中TBA 值的决定系数R2最高,均大于0.968,POV值的拟合精度较高。由此可见一级化学反应动力学方程能够较好地描述冷藏期间鲈鱼背、腹部鱼片的脂肪氧化情况。

表5 真空包装的鲈鱼背、腹部鱼片脂肪氧化指标非线性拟合参数Table 5 Nonlinear fitting parameters of weever on fat oxidation indicators

将2 ℃(271 K)、0 ℃(273 K)和4 ℃(277 K)贮藏条件下得到的鲈鱼鱼片的AV值、POV值和TBA值的反应速率常数k用Arrhenius方程在不同贮藏温度下拟合,可得到 k0及 EA的数值如下:

表6中,鲈鱼背、腹部各脂肪氧化指标中的相关系数均较大(R2≥0.830),其中,背部 AV 值的相关系数最大为0.999;腹部AV值和POV值的相关系数都是0.993。从表6中还可以看出,背部鱼片AV、POV、TBA指标的活化能EA均小于腹部鱼片。这可能是因为在真空包装条件下这三个指标的变化在背部鱼肉中更容易实现的缘故。

表6 反应速率常数k对不同温度的拟合结果Table 6 The fitting results of reaction rate constant k for temperature

AV值、POV值和TBA值的反应速率常数k值随温度的变化规律与空气包装条件下的变化规律相同。其中POV值的拟合精度最高,表明与Arrhenius 方程的拟合效果最好。这三个实验指标的化学反应速率常数k的数据点均可以很好的与拟合曲线相拟合。这表明用Arrhenius方程可以很好的表达在不同的贮藏温度下鲈鱼片AV值、POV值和TBA值的反应速率常数随贮藏温度变化的关系(图6)。

图6 三种氧化指标的反应速率常数k与Arrhenius方程拟合Fig.6 Oxidation indicators reaction rate constant k fitting with Arrhenius equation

根据鲈鱼背、腹部鱼片在3个贮藏温度下所得到的k0值与EA值带入式(5),可以计算出真空包装中鲈鱼背、腹部鱼片在贮藏过程中AV值、POV值和TBA值的货架期预测模型:

关于AV值的货架期预测模型:

关于POV值的货架期预测模型:

关于TBA值的货架期预测模型:

其中,A0:初始腹部TBA值,A:储藏 t时间后的腹部TBA值。

鲈鱼背、腹部鱼片在贮藏温度下货架期预测模型的验证与评价。贮藏温度下鲈鱼背、腹部鱼片的感官评价以5分作为可接受的最低限值。当冷藏鲈鱼背、腹部鱼片TBA值超过1 mg MDA ·kg-1时,用所建立的鲈鱼货架期预测模型的预测值与实验值对货架期模型进行验证(表7)。

由表7可以看出在真空包装条件下所建立的关于三个指标的实际值与预测值的相对误差在±10.06%内。这表明在真空包装条件下所建立的货架期预测模型可以很好地预测贮藏期间鲈鱼片的脂肪氧化变化和货架期。其中鲈鱼背部鱼片AV值的平均相对误差较小为2.3%,其次是背部POV值为3.23%;而腹部TBA值的平均相对误差最大为7.97%。这表明在真空包装条件下所建立的关于鲈鱼背部的AV值和POV值的货架期模型较为准确。

表7 真空包装下鲈鱼货架期模型验证结果Table 7 The shelflife model validation results of weever on fillets in the vacuum package

3 结论

在三个不同的贮藏温度下(2、0和4 ℃)将鲈鱼片的脂肪氧化指标 AV值、POV值和TBA值与一级化学反应动力学模型拟合发现在空气和真空包装方式下关于AV值和POV值有很好的拟合效果,TBA值拟合效果较好。所以选择将AV值、POV值和TBA值作为货架期评判指标建立预测模型。通过对模型验证后得到预测值与实验值相对误差均在±10.60%以内,且真空包装组所建立的三个值的预测模型平均相对误差较小,但普通包装组所得的货架期模型平均相对误差也在允许的范围内。因此,根据鲈鱼在贮藏期间的脂肪氧化指标所建立的预测模型能很好地预测冷藏温度下鲈鱼背、腹部鱼片的货架期,可为预测和控制鲈鱼片在贮藏期间的品质提供现实参考。

[1]邓锦峰,王安利,周初霞,等.鲈鱼的营养研究进展[J].饲料工业,2006(10):5960.

[2]Pan S, Wu S. Effect of chitooligosaccharides on the denaturation of weever myofibrillar protein during frozen storage[J]. International Journal of Biological Macromolecules,2014(64):549552.

[3]赵宏强,吴金鑫,张苑怡,等.超高压处理对冷藏鲈鱼片品质及组织结构变化的影响[J].高压物理学报,2017(4):494504.

[4]余丹丹,安康鱼调味烤鱼片加工工艺研究与货架期预测[D].浙江:浙江海洋大学,2016.

[5]刘卓然,翁佩芳,敖雨燕,等.动力学模型预测泥鳅(Misgurnusanguillicaudatus)半干制品的货架期[J].海洋与湖藻,2014(2):382386.

[6]吴雪丽, 刘红英, 韩冬娇. 扇贝贮藏货架期预测模型的建立与评价[J].食品科学,2014(22):315319.

[7]刘文丽, 邢少华, 车长远, 等. 变温条件下大黄鱼品质变化动力学模型的建立与评价[J].食品科学,2015(24):286289.

[8]周果, 崔燕, 杨文鸽, 等. 冰温贮藏对梭子蟹品质影响及其货架期模型的建立[J].核农学报,2017(4):719727.

[9]郭全友, 王锡昌, 杨宪时,等. 不同贮藏温度下养殖大黄鱼货架期预测模型的构建[J].农业工程学报,2012(10):267273.

[10]GB/T 5009. 372003. 食用植物油卫生标准的分析方法[S].北京：中国标准出版社，2013.

[11]Salih A M, Smith D M, Price J F, et al. Modified extraction 2thiobarbituric acid method for measuring lipid oxidation in poultry[J]. Poultry Science, 1987(66):14831488.

[12]Hui G, Liu W, Feng H,et al.Effects of chitosan combined with nisin treatment on storage quality of large yellow croaker(Pseudosciaenacrocea)[J]. Food Chemistry,2016(203):276282.

[13]杨文鸽. 大黄鱼冰藏期间 ATP 关联物含量变化及其鲜度评价[J].农业工程学报,2007(06):217222.

[14]Rao Y, Xiang B, Zhou X, et al. Quantitativeand qualitative determination of acid value of peanutoil using nearinfrared spectrometry[J]. Journal of Food Engineering,2009(93):249256.

[15]Cebi N., Yilmaz M T, Sagdic O, et al.Prediction of peroxide value in omega3 rich microalgae oil by ATRFTIR spectroscopy combined with chemometrics[J].Food Chemistry,2017(225):188196.

[16]王悦齐,吴燕燕,李来好,等.抗氧化乳酸菌对发酵腌干带鱼脂肪氧化的影响及其主成分分析[R].中国水产学会会议论文集,2016:400401.

[17]Giménez B, GómezGuillén M C, PérezMateos M, et al.Evaluation of lipid oxidation in horse mackerel patties covered with borage containing film during frozen storage[J].Food Chemistry,2011(124):13931403.

[18]鞠健, 汪超, 廖李,等.低剂量辐照对冷藏期间鲈鱼品质的影响[J].食品科技,2016(7):157162.

[19]Salih A M, Smith D M, Price J F, et al. Modified extraction 2thiobarbituric acid method for measuring lipid oxidation in poultry[J]. Poultry Science, 1987(66):14831488.

[20]Ojagh S M, Rezaei M, Razavi S H,et al. Effect of chitosan coatings enriched with cinnamon oil on the quality of refrigerated rainbow trout[J].Food Chemistry,2010(120):193198.

[21]李婷婷,丁婷,胡文忠,等.冷藏大西洋鲑背腹部鱼片脂肪氧化变化规律及货架期模型的建立[J].水产学报,2015(2):265274.

[22]田秋实,谢晶, 励建荣. 白鲢鱼糜制品货架期模型的建立[J].食品工业科技,2009(4):7073.

[23]Labuza T P. Application of chemical kinetics todeterioration of food[J].Journal of Chemical Education,1984(61):348351.

[24]陈睿, 徐幸莲, 周光宏.真空包装鸡肉早餐肠货架期预测模型的建立[J].食品科学,2014(4):209213.