茶油掺伪的数学回归模型研究

1、茶油掺伪的数学回归模型研究彭思敏1林谦2 吴卫国1（食品科学与生物技术湖南省重点试验室1；湖南农业大学食品科技学院，长沙 410128）（中国农业科学院麻类研究所2，长沙 410205）摘 要：本研究以纯茶油中掺伪不同比例大豆油、花生油和菜籽油中任意两种植物油脂来建立茶油掺伪模型收稿日期：2017-02-19作者简介：彭思敏，女，1987年出生，硕士，粮食、油脂与植物蛋白工程通信作者：吴卫国，男，1968年出生，教授，粮食、油脂与植物蛋白工程，实验通过气质联用技术测得各样品脂肪酸的含量，并通过相关性、显著性分析来筛选出合适的标志性脂肪酸建立茶油掺伪的多元回归方程模型。结果表明：某些脂肪酸含量在

2、茶油掺伪模型中变化显著，且与茶油掺伪模型间具有高的相关性，选取该部分脂肪酸建立的相应回归模型方程的相关系数（R2值）均较高，说明本实验获得的回归方程模型真实、可靠。研究结果可以为茶油的产品质量监控和掺假检测提供理论依据。关键词：茶油；掺伪；相关性分析；显著性分析；标志性脂肪酸；回归预测中图分类号：TS225.1 文献标识码：A 文章编号：20170219-1修改稿Mathematical Regression Model of Camellia Oil AdulterationPeng Simin1 Lin Qian2 Wu Weiguo1(Food Science and Biotechno

3、logy of Key Laboratory of Hunan Province1；College of Food Science and Technology of Hunan Agricultural University,Changsha 410128)（Institute of Bast Fiber Crops,Chinese Academy of Agricultural science，Changsha 410205）Abstract：Camellia oil adulteration model was established by using any two kinds of

4、vegetable oils (soybean oil, peanut oil and rapeseed oil) with different proportions in pure camellia oil. The content of fatty acids in each sample was measured by gas chromatography-mass spectrometry (GC-MS). Correlation and significance analysiswere used to screen out the appropriate symbolic fat

5、ty acid to establish the multiple regression equation model of camellia oil adulteration. The results showed that some fatty acids in camellia oil adulteration model changes significantly, and has a high correlation with the camellia oil adulteration model.At the same time, the correlation coefficie

6、nt (R2vaule) of the regression equationswere high, which indicated that the regression equation is true and reliable. The results could provide theoretical basis for quality control and adulterate detection of camellia oil.Key Words：camellia oil;adulterate；correlation analysis; HYPERLINK app:ds:sign

7、ificance%20analysis significant analysis; appropriate symbolic fatty acid； HYPERLINK app:ds:regression regression-forecast HYPERLINK app:ds:model 茶油因其脂肪酸组成符合人体营养需求，同时富含多种活性物质，食用价值高，也因此在价格上较普通食用油更高，又因近几年气候影响导致茶籽减产，其价格更一度超过进口橄榄油。所以一些不良商家将某些低价油脂（如菜籽油、大豆油等）混入茶油中掺假售卖。 目前，用于油脂掺伪检测的方法有理化特性分析法1,2，色谱、光谱分析法，化

8、学计量法3、气质联用法4-7等。理化特性分析法仅适用于特定油脂的掺伪检测，不具有普遍适用性，且操作复杂，对人员操作技能具有一定要求。化学计量法是通过分析油脂量化的特征反应进行鉴定的，此法专一性较强，但局限于鉴别某一种特定油脂类型的掺假，且需要通过油脂样品的大数据库进行验证，该法仅适用于初步定性分析油脂的掺伪现象；气质联用技术则通过获取油脂脂肪酸种类及含量组成，并结合特征脂肪酸进行定性定量分析，不过就实际而言，特征脂肪酸的选择具有一定局限性，以菜籽油为例，随着油料品种的改良，低芥酸菜油已逐步占据市场，这使得芥酸不能作为菜籽油的特征脂肪酸用于菜籽油的鉴别分析中。不同种类的油脂所含脂肪酸种类及含量具

9、有相似性也具有差异性，不同油脂之间的掺伪行为会导致其脂肪酸组成发生相应的变化。本研究借此原理，分析各脂肪酸与油脂含量之间的相关性，相关性系数值越大，则表明该脂肪酸越重要，在茶油掺伪现象鉴别中具有极强影响力。同时，通过不同掺伪油脂含量对各脂肪酸含量的差异性分析，可辨别出在茶油掺伪现象中变化明显的脂肪酸，差异性越强，越具有代表性8-10，从而综合筛选出对茶油掺伪检测具有较大影响的脂肪酸种类。在此基础上，深究不同茶油掺伪模型中各脂肪酸含量变化间的差异性，建立了茶油掺伪油脂含量和具有代表性的脂肪酸之间的数学回归模型11-15，在检测茶油掺伪现象时，通过检测得出相应标志性脂肪酸的含量，依次代入方程组即可

10、获得相应掺伪油脂含量的相关值，在一定程度上实现了茶油掺伪检测技术的科学化，为茶油掺伪鉴别的标准建立提供了有效的理论基础。1 材料与方法1.1 材料与试剂茶油：市售；大豆、花生、菜籽样品：湖南农业大学实验基地；正己烷、甲醇：色谱纯；无水硫酸钠、H2SO4、KOH：分析纯。1.2 设备与仪器SHZ-82水浴恒温振荡器；HH-8数显恒温水浴锅；GCMS-2010气质联用仪（SHIMADZU）。1.3 实验方法1.3.1 样品提取索氏抽提法进行纯油样品提取，具体操作流程参照参考文献16报道进行。1.3.2 茶油掺伪模型设计 根据二因子四水平的正交实验设计原理，取任意两种纯油样品（菜籽油、大豆油、花生油

11、）按5%、15%、25%、50%的梯度加入到茶油油脂中，模拟茶油油脂掺假（见表1）。表1 植物油脂掺假模型设计基础油茶油茶油茶油用于掺假的油脂花生油，大豆油花生油，菜籽油大豆油，菜籽油1.3.3 脂肪酸甲酯化处理运用酸碱结合法进行处理，具体操作流程参照参考文献16报道进行。1.3.4 气质联用分析 GC-MS分析条件为：HB-88毛细管柱（100.0 m 0.32 mm 0.25 m）；载气为He（纯度99.999%)；进样口温度：250 ；柱温：初始温度120 ，保持1 min、10 /min升温至175 ，保持10 min、5 /min升温至210 ，保持5 min、5 /min升温至23

12、0，保持15min；柱流速：1.04 ml/min，分流比：10:0；离子源温度：200 ；接口温度：220 ；质核比扫描范围：50700 m/z；进样量：1 l。1.3.5 脂肪酸含量分析 通过面积归一化法测得样品油脂中各脂肪酸的相对含量。1.3.6 相关性及显著性分析以样品中棕榈酸、肉豆蔻酸、油酸、亚油酸、硬脂酸、亚麻酸、花生酸、花生烯酸和芥酸17九类脂肪酸甲酯的相对含量以及掺伪模型中茶油含量为变量，采用SPSS 17.0 数据统计软件进行双变量相关分析，根据相关性大小判断出随着茶油含量变化而具有较高程度关联的脂肪酸种类。同时结合单因素方差分析（ANOVA），若组间差异显著，则采用Dunc

13、an氏法进行多重比较，显著水平为0.05，极显著水平为0.01。实验结果以平均值标准差表示。1.3.7 回归模型建立 以筛选出的具有显著相关性及差异性的脂肪酸的含量为Y值，茶油掺伪模型中掺入的油脂含量为X值，经SAS逐步回归分析建立回归方程组。2 结果与分析2.1 茶油掺伪菜籽油和大豆油模型的回归方程建立2.1.1 掺伪模型中不同茶油含量与脂肪酸含量的相关性分析表2为茶油掺伪模型中不同的茶油含量水平与其中9种脂肪酸含量间的相关性分析结果，从表2可知，各茶油掺伪模型中不同的茶油含量水平除与其中顺式油酸的含量变化呈极显著正相关（R=0.872，P0.01）外，其他多呈极显著负相关。而其中又与亚麻酸

14、的含量变化呈现的负相关关系最强（R=-0.978，P0.01），其次是花生酸（R=-0.943，P0.01）和顺式亚油酸（R=-0.845，P0.01）。说明在各茶油掺伪模型中茶油的含量越高，则顺式油酸含量越高，而亚麻酸、花生酸和顺式亚油酸含量则越低。 表2 茶油掺伪模型中各脂肪酸含量与茶油含量间的相关系数分析百分含量茶油肉豆蔻酸棕榈酸硬脂酸顺式油酸顺式亚油酸花生酸亚麻酸花生烯酸芥酸茶油1肉豆蔻酸-0.0031棕榈酸-0.003-0.0351硬脂酸-0.692*0.0250.430*1顺式油酸0.872*-0.025-0.159-0.911*1顺式亚油酸-0.845*0.0220.2440.9

15、39*-0.982*1花生酸-0.943*-0.123-0.441*0.560*-0.733*0.698*1亚麻酸-0.978*0.015-0.1050.820*-0.941*0.933*0.893*1花生烯酸-0.603*-0.044-0.861*-0.016-0.1650.1190.740*0.456*1芥酸0.2090.596*0.353*0.0210.028-0.038-0.269-0.152-0.405*1注：表中“*”表示显著相关（P0.05），“*”表示极显著相关（P0.01），下同。2.1.2 茶油掺伪模型中各脂肪酸含量的显著性分析表3 掺假菜籽油和大豆油的茶油中各脂肪酸含量分

16、析掺假5%菜油（实验I组）Exp.1掺假15%菜油（实验II组）Exp.2掺假25%菜油（实验III组）Exp.3掺假50%菜油（实验IV组）Exp.4茶油掺假5%大豆油中脂肪酸百分含量棕榈酸甲酯8.80.06Aa8.740.04Aa9.090.10ABa9.500.26Bb硬脂酸甲酯2.490.02Aa2.520.01Aa3.250.03Bb3.410.05Bb顺式油酸甲酯74.90.57Dd69.00.40Cc62.80.74Bb49.81.35Aa顺式亚油酸甲酯11.350.08Aa16.290.09Bb19.90.23Cc29.20.79Dd花生酸甲酯0.090.02Aa0.130.0

17、1Bb0.230.00Cc0.350.00Dd亚麻酸甲酯1.000.00Aa1.890.05Bb2.970.01Cc4.870.13Dd花生烯酸甲酯0.76 0.00Bb0.750.00Bb0.750.01Bb0.590.01Aa芥酸甲酯0.000.02Aa0.120.01Bb0.180.00Cc0.300.00Dd茶油掺假15%大豆油中脂肪酸百分含量棕榈酸甲酯8.100.14Aa8.350.14Aab8.650.21ABb9.280.19Bc硬脂酸甲酯2.040.04Aa2.350.04Bb2.790.07Cc3.610.07Dd顺式油酸甲酯74.171.37Cd69.371.26Cc62.

18、691.55Bb48.170.97Aa顺式亚油酸甲酯11.680.22Aa15.290.28Bb19.520.48Cc29.840.60Dd花生酸甲酯0.140.00Aa0.190.00Bb0.260.01Cc0.350.00Dd亚麻酸甲酯1.510.03Aa2.140.05Bb3.220.08Cc5.730.11Dd花生烯酸甲酯0.83 0.01Bc0.720.01Aa0.800.02ABbc0.770.01ABb芥酸甲酯0.240.11Aa0.330.01Bb0.380.01Cc0.710.01Dd茶油掺假25%大豆油中脂肪酸百分含量棕榈酸甲酯8.060.21ABab7.720.14Aa8

19、.550.25ABbc8.810.24Bc硬脂酸甲酯2.490.06Aa2.580.04Aa3.000.09Bb3.640.10Cc顺式油酸甲酯70.01.85Cc66.241.25BCc59.381.73Bb46.81.35Aa顺式亚油酸甲酯13.480.35Aa16.980.34Bb21.080.61Cc31.250.90Dd花生酸甲酯0.320.02Bb0.230.01Aa0.370.00BCc0.400.00Cc亚麻酸甲酯2.230.05Aa3.00.05Bb3.940.11Cc6.010.17Dd花生烯酸甲酯1.040.021.00.021.010.020.980.02芥酸甲酯0.4

20、90.01Aa0.640.01Bb0.570.01Cc0.000.00Dd茶油掺假50%大豆油中脂肪酸百分含量棕榈酸甲酯6.870.26a7.000.21ab7.330.27ab7.630.27c硬脂酸甲酯2.350.09Aa2.670.07ABb2.900.10Bb3.650.04Cc顺式油酸甲酯67.162.65Cc62.451.97BCc55.982.10Bb42.761.56Aa顺式亚油酸甲酯15.30.61Aa19.410.64Bb23.70.86Cc33.521.21Dd花生酸甲酯0.430.01Aa0.440.01Aa0.540.02Bb0.700.02Cc亚麻酸甲酯3.690.

21、12Aa4.450.11Ab5.520.20Bc7.880.28Cd花生烯酸甲酯1.380.041.370.051.320.041.310.04注：同行数据肩标不同小写字母表示差异显著（P0.05），不同大写字母表示差异极显著（P0.01），反之则不显著，下同。 从表3可知：各个掺伪阶段，各实验组间顺式亚油酸甲酯、亚麻酸甲酯含量差异均明显（P0.01）。当掺伪大豆油含量固定情况下，其含量均以实验IV组最高，其次为III、II、I组。当掺伪菜籽油含量固定情况下，顺式亚油酸甲酯、亚麻酸甲酯含量均随大豆油含量的增加呈递增规律；茶油在掺伪15%以下大豆油时，顺式油酸甲酯、花生酸甲酯浓度随着掺伪菜籽油浓

22、度的增加在各实验组间分别呈递减、递增变化，但各组间两两比较，差异均极显著（P0.01）；芥酸甲酯含量极低，在茶油掺伪25%以下大豆油时，其含量随着掺伪菜籽油浓度的增加而增加，且各实验组间均差异极显著（P0.01）。2.1.3 标志性脂肪酸的确定及回归方程的建立综上，在茶油掺伪大豆油和菜籽油掺伪模型中，同时具有显著相关性及差异性的脂肪酸中，芥酸甲酯在茶油掺伪5%大豆油和菜籽油的情况下存在不易被检出的可能性；顺式油酸甲酯、花生酸甲酯仅在茶油掺伪15%以下大豆油时具有检测性。从而优化选出顺式亚油酸甲酯、亚麻酸甲酯2种脂肪酸作为标志性脂肪酸。同时以此2种脂肪酸的含量为Y值，以茶油掺伪模型中掺入的菜籽油

23、含量为X值，经过SAS分析法建立了以下回归方程组：Y（顺式亚油酸甲酯）=8.8112+0.0925X1+0.3998 X2（R2=0.9916，P0.01）。Y（亚麻酸甲酯）=0.1524+0.0611X1+0.0908X2（R2=0.9926，P0.01）2.2 茶油掺伪大豆油和花生油模型的回归方程建立2.2.1 掺伪模型中不同茶油含量与脂肪酸含量的相关性分析表4为茶油掺伪模型中不同的茶油含量水平与其中9种脂肪酸含量间的相关性分析结果，从表4可知，各茶油掺伪模型中不同的茶油含量水平除与其中顺式油酸的含量变化呈极显著正相关（R=0.985，P0.01）外，其他多呈极显著负相关。而其中又与顺式亚

24、油酸的含量变化呈现的负相关关系最强（R=-0.991，P0.01），其次是棕榈酸（R=-0.944，P0.01）和硬脂酸（R=-0.897，P0.01）。说明在各茶油掺伪模型中茶油的含量越高，则顺式油酸含量越高，而顺式亚油酸、棕榈酸和硬脂酸含量则越低。 表4 茶油掺伪模型中各脂肪酸含量与茶油含量间的相关系数分析百分含量茶油肉豆蔻酸棕榈酸硬脂酸顺式油酸顺式亚油酸花生酸亚麻酸花生烯酸茶油1肉豆蔻酸-0.1251棕榈酸-0.944*0.0931硬脂酸-0.897*0.1080.931*1顺式油酸0.985*-0.104-0.961*-0.949*1顺式亚油酸-0.991*0.1220.960*0.9

25、44*-0.998*1花生酸-0.818*0.1210.687*0.494*-0.730*0.744*1亚麻酸-0.712*0.0760.782*0.940*-0.805*0.793*0.1891花生烯酸-0.363*0.0920.193-0.050-0.2350.2570.801*-0.355*12.2.2 茶油掺伪模型中各脂肪酸含量的显著性分析表5 掺假大豆油和花生油的茶油中各脂肪酸含量分析掺假5%豆油（实验I组）Exp.1掺假15%豆油（实验II组）Exp.2掺假25%豆油（实验III组）Exp.3掺假50%豆油（实验IV组）Exp.4茶油掺假5%花生油中脂肪酸百分含量棕榈酸甲酯9.32

26、0.01Aa9.440.12Aab9.920.36ABb10.570.06Bc硬脂酸甲酯2.420.03Aa2.640.03ABa2.980.14Bb3.730.07Cc顺式油酸甲酯73.20.32Dd66.20.11Cc61.80.89Bb46.40.36Aa顺式亚油酸甲酯12.70.09Aa16.10.05Bb20.70.31Cc30.40.01Dd花生酸甲酯0.050.07a0.140.01ab0.140.04ab0.260.03b亚麻酸甲酯0.080.03Aa1.580.01Bb2.630.03Cc4.930.14Dd花生烯酸甲酯0.65 0.03Bb0.590.03ABab0.590

27、.01ABab0.530.00Aa茶油掺假15%花生油中脂肪酸百分含量棕榈酸甲酯9.190.09Aa9.710.09Bb9.810.07Bb10.80.05Cc硬脂酸甲酯2.380.09Aa2.670.06Bb2.950.00Cc3.830.02Dd顺式油酸甲酯70.40.53Dd63.61.37Cc58.80.10Bb43.70.28Aa顺式亚油酸甲酯14.70.02Aa18.90.28Bb22.80.07Cc32.30.20Dd花生酸甲酯0.240.02Aa0.300.04ABab0.330.03ABbc0.410.00Bc亚麻酸甲酯0.730.07Aa1.570.11Bb2.370.07

28、Cc4.710.09Dd花生烯酸甲酯0.70 0.070.720.070.740.080.580.02茶油掺假25%花生油中脂肪酸百分含量棕榈酸甲酯9.80.12Aa9.980.12Aab10.20.17Ab11.10.19Bc硬脂酸甲酯2.550.02Aa2.850.02Bb3.160.07Cc3.990.04Dd顺式油酸甲酯63.11.2Dd56.60.36Cc53.50.10Bb38.50.28Aa顺式亚油酸甲酯19.20.06Aa23.40.08Bb26.60.10Cc36.00.14Dd花生酸甲酯0.470.00a0.510.02ab0.530.02ab0.570.02b亚麻酸甲酯0

29、.640.00Aa1.60.04Bb2.50.05Cc4.880.00Dd花生烯酸甲酯1.60.83Bc1.60.81Bc1.60.74ABb1.60.67Aa茶油掺假50%花生油中脂肪酸百分含量肉豆蔻酸甲酯0.050.01a0.020.00ab0.050.01ab0.060.00b棕榈酸甲酯10.20.02a10.40.17a10.80.12ab11.20.50b硬脂酸甲酯2.770.08Aa3.180.00Bb3.430.00Bc4.330.05Cd顺式油酸甲酯54.60.81Dd47.40.09Cc43.10.62Bb28.20.09Aa顺式亚油酸甲酯25.80.16Aa30.20.10

30、Bb33.50.11Cc43.90.26Dd花生酸甲酯0.810.010.880.090.900.050.920.00亚麻酸甲酯0.660.07Aa1.640.01Bb2.450.05Cc5.190.00Dd花生烯酸甲酯0.930.06ab0.960.00b0.950.03b0.830.01a 从表5可知：各个掺伪阶段，各实验组间顺式油酸甲酯、顺式亚油酸甲酯、亚麻酸甲酯含量差异均明显（P0.05），且均明显低于实验III、IV组（P0.05），茶油掺伪15%以上花生油时，硬脂酸甲酯浓度随着掺伪大豆油浓度的增加在各实验组间递增，且各组间均差异极显著（P0.05）。2.2.3 标志性脂肪酸的确定及

31、回归方程的建立综上，在茶油掺伪花生油和大豆油掺伪模型中，亚麻酸甲酯虽具有极显著差异性，但因其相关性系数为-0.712，数值较低，相关性不强，不作为选择项。硬脂酸甲酯在茶油中掺伪小含量的花生油和大豆油情况下不易分辨掺伪变化的差异性，但在绝大多数情况下均具有极显著差异性。从而优化选出顺式油酸甲酯、顺式亚油酸甲酯、硬脂酸甲酯3种脂肪酸作为标志性脂肪酸。同时以此3种脂肪酸的含量为Y值，以茶油掺伪模型中掺入的菜籽油含量为X值，经过SAS分析法建立了以下回归方程组： Y（顺式油酸甲酯)=77.9941-0.4237X1-0.5713X2（R2=0.9965，P0.01）Y（顺式亚油酸)=9.0384+0.

32、3035X1+0.3896X2（R2=0.9914，P0.01）Y（硬脂酸)=2.0873+0.0114X1+0.0320X2（R2=0.9855，P0.01）2.3 茶油掺伪菜籽油和花生油模型的回归方程建立2.3.1 掺伪模型中不同茶油含量与脂肪酸含量的相关性分析表6为茶油掺伪模型中不同的茶油含量水平与其中9种脂肪酸含量间的相关性分析结果，从表6可知，各茶油掺伪模型中不同的茶油含量水平除与其中顺式油酸的含量变化呈极显著正相关（R=0.902，P0.01）外，其他多呈极显著负相关。而其中又与花生酸的含量变化呈现的负相关关系最强（R=-0.935，P0.01），其次是花生烯酸（R=-0.929，

33、P0.01）和顺式亚油酸（R=-0.871，P0.01）。说明在各茶油掺伪模型中茶油的含量越高，则顺式油酸含量越高，而花生酸、花生烯酸和顺式亚油酸含量则越低。表6 茶油掺伪模型中各脂肪酸含量与茶油含量间的相关系数分析百分含量茶油肉豆蔻酸棕榈酸硬脂酸顺式油酸顺式亚油酸花生酸亚麻酸花生烯酸芥酸茶油1肉豆蔻酸0.2081棕榈酸0.2900.1641硬脂酸-0.734*0.0300.393*1顺式油酸0.902*0.146-0.038-0.821*1顺式亚油酸-0.871*-0.1080.2120.947*-0.951*1花生酸-0.935*-0.1300.0540.901*-0.958*0.982*

34、1亚麻酸-0.685*-0.244-0.853*0.096-0.3360.2530.399*1花生烯酸-0.929*-0.222-0.523*0.551*-0.702*0.673*0.772*0.867*1芥酸甲酯-0.738*-0.287-0.802*0.172-0.409*0.3310.475*0.983*0.894*12.3.2 茶油掺伪模型中各脂肪酸含量的显著性分析 表7 掺假菜籽油和花生油的茶油中各脂肪酸含量分析掺假5%花生（实验I组）Exp.1掺假15%花生（实验II组）Exp.2掺假25%花生（实验III组）Exp.3掺假50%花生（实验IV组）Exp.4茶油掺假5%菜籽油中脂肪

35、酸百分含量肉豆蔻酸甲酯0.050.00Bb0.040.00Aa0.060.00Cc0.050.04BCc棕榈酸甲酯9.260.09a9.420.12ab9.690.26ab10.000.32b硬脂酸甲酯2.260.02Aa2.340.03Aab2.530.06ABb2.730.07Bc顺式油酸甲酯74.930.79Cc70.941.30Cc66.41.89Bb57.81.36Aa顺式亚油酸甲酯11.160.13Aa13.990.25Bb17.10.51Cc24.30.74Dd花生酸甲酯0.190.03Aa0.370.01Bb0.510.01Cc0.920.03Dd亚麻酸甲酯0.560.00Bb

36、0.530.01ABab0.490.03Aa0.510.01ABa花生烯酸甲酯0.760.01Aa0.880.01Bb0.900.02Bb1.130.03Cc芥酸甲酯0.060.03Aa0.150.01Bb0.100.03Cc0.130.14Dd茶油掺假15%菜籽油中脂肪酸百分含量棕榈酸甲酯8.690.23b8.680.26b9.260.25ab9.690.32a硬脂酸甲酯2.260.06Aa2.470.07ABab2.480.07ABb2.750.09Ac顺式油酸甲酯72.41.9Bc69.122.1Bbc65.31.8Bb56.011.9Aa顺式亚油酸甲酯12.060.32Aa14.70.

37、45Ab17.870.51Bc25.750.87Cd花生酸甲酯0.230.00Aa0.420.01Bb0.610.01Cc0.960.03Dd亚麻酸甲酯1.260.03Bb1.140.03ABa1.120.03ABa1.000.03Aa花生烯酸甲酯0.90.02Ba1.10.03ABa1.10.03ABa1.180.04Ab芥酸甲酯0.160.00Aa0.170.01ABb0.230.00Bc0.180.00Cd茶油掺假25%菜籽油中脂肪酸百分含量棕榈酸甲酯8.140.24a8.430.26ab8.730.31ab9.150.39b硬脂酸甲酯2.270.07Aa2.430.08ABab2.54

38、0.09ABbc2.781.10Bc顺式油酸甲酯71.921.2Bc67.660.36Bbc63.690.10ABb54.320.28Aa顺式亚油酸甲酯12.300.36Aa15.460.48ABb18.70.60Bc26.491.14Cd花生酸甲酯0.290.00Aa0.510.01Bb0.600.02Bc1.060.04Cd亚麻酸甲酯1.610.041.750.041.700.051.620.06花生烯酸甲酯1.040.02Aa1.20.02ABb1.190.04ABb1.310.05Bb茶油掺假50%菜籽油中脂肪酸百分含量肉豆蔻酸甲酯0.050.01a0.040.00b0.040.01a

39、b0.040.00b棕榈酸甲酯7.090.2a7.530.28ab7.740.30ab8.250.33b硬脂酸甲酯2.310.08Aa2.440.09ABa2.50.10ABab2.80.12Bb顺式油酸甲酯67.742.25Bc63.892.2Bbc59.852.38ABc50.322.34Aa顺式亚油酸甲酯15.080.16Aa17.840.10ABb20.980.11Bc28.860.26Cd花生酸甲酯0.550.01Aa0.700.09Ab0.900.05Bc1.20.00Cd亚麻酸甲酯3.060.093.020.113.080.123.030.14花生烯酸甲酯1.340.04Aa1.

40、410.05Aab1.530.07ABb1.730.08Bc芥酸甲酯0.520.02a0.480.01ab0.540.02ab0.550.02b从表7可知：各个掺伪阶段，各实验组间顺式亚油酸甲酯、花生酸甲酯含量差异均明显（P0.01），均以实验IV组最高，其次为III、II、I组；在茶油掺假15%以下花生油基础上掺假5%、15%、25%、50%的菜籽油时，芥酸甲酯含量在实验各组两两相比较，差异均极显著（P0.05）。花生烯酸甲酯、亚麻酸甲酯其含量变化无规律性，各组间也无明显差异性。2.3.3 标志性脂肪酸的确定及回归方程的建立综上，在茶油掺伪花生油和菜籽油掺伪模型中，同时具有显著相关性及差异性

41、的脂肪酸中，芥酸甲酯在茶油掺伪15%以下大豆油情况下掺伪菜籽油时，具有检测性。顺式油酸甲酯、花生烯酸甲酯、亚麻酸甲酯因其不具有显著差异性而不列为脂肪酸选择项。从而优化选出顺式亚油酸甲酯、花生酸甲酯2种脂肪酸作为标志性脂肪酸。同时以此2种脂肪酸的含量为Y值，以茶油掺伪模型中掺入的菜籽油含量为X值，经过SAS分析法建立了以下回归方程组：Y（顺式亚油酸)= 8.9184+0.0884X1+0.3065X2（R2=0.9904，P0.01）Y（花生酸) =0.0738+ 0.0075X1+0.0159X2（R2=0.9875，P0.01）3结论本研究通过筛选出具有显著相关性和差异性的脂肪酸作为标志性脂

42、肪酸，并以脂肪酸的含量为Y值，以茶油掺伪模型中掺入的植物油含量为X值，经过SAS分析法建立的回归方程组，相关系数R2值均较高，说明回归方程组关系真实且可靠。在茶油掺伪大豆油和菜籽油掺伪模型中，优化选出顺式亚油酸甲酯、亚麻酸甲酯2种标志性脂肪酸，建立了以下回归方程组：Y（顺式亚油酸甲酯）=8.8112+0.0925X1+0.3998 X2（R2=0.9916，P0.01）。Y（亚麻酸甲酯）=0.1524+0.0611X1+0.0908X2（R2=0.9926，P0.01）在茶油掺伪花生油和大豆油掺伪模型中，优化选出顺式油酸甲酯、顺式亚油酸甲酯、硬脂酸甲酯3种标志性脂肪酸，建立了以下回归方程组：Y

43、（顺式油酸甲酯)=77.9941-0.4237X1-0.5713X2（R2=0.9965，P0.01）Y（顺式亚油酸)=9.0384+0.3035X1+0.3896X2（R2=0.9914，P0.01）Y（硬脂酸)=2.0873+0.0114X1+0.0320X2（R2=0.9855，P0.01）在茶油掺伪花生油和菜籽油掺伪模型中，优化选出顺式亚油酸甲酯、花生酸甲酯2种标志性脂肪酸，建立了以下回归方程组：Y（顺式亚油酸)= 8.9184+0.0884X1+0.3065X2（R2=0.9904，P0.01）Y（花生酸) =0.0738+ 0.0075X1+0.0159X2（R2=0.9875，P

44、0.01）本实验所得回归方程可以把茶油掺伪油脂含量与脂肪酸含量之间的关系量化，在检测茶油掺伪现象时，检测得出相应脂肪酸的含量，代入方程即可获得相应掺伪油脂的相关信息，在一定程度上实现了茶油掺伪检测技术的科学化。参考文献1谢婷婷，黄丽. 山茶油掺假检测技术的研究现状J.科技与创新，2014(21):9-13XIE Tingting，HUANG Li.Research status camellia oil adulteration detection technologyJ.Science and Technology & Innovation,2014(21):9-132支敏，陆鑫，李红霞.

45、食用油品质检测方法的综述J.河北联合大学学报（自然科学版），2014,1(36):109-112ZHI Min，LU Xin，LI Hongxia. Review on detecting methods about the quality of edible oilJ.Journal of Hebei United university natural sclence edition,2014,1(36):109-1123吴卫国，刘真知，彭思敏.基于特征脂肪酸及脂肪酸比值的食用植物油掺假判别J.食品科学，2013,34:270-273WU Weiguo，LIU Zhengzhi，PENG S

46、imin,et al.Adulteration identification of edible vegetable oils based on fatty acid profileJ. Food Science,2013,34:270-2734Miloudi,Hilali,Zoubida,et al. Detection of argan oil adulteration using quantitative campesterol GC-analysisJ. Journal of the American Oil Chemists Society, 2007,(8)4:761-7645严晓

47、丽,徐昕. 气相色谱法鉴别掺假山茶油定性及定量研究J. 食品工程，2011,(2):47-49YAN Xiaoli，XU Xin.A qualitative & quantitative study of distinguishing adulteration in camellia oil by gas chromatography J. Food Engineering,2011,(2):47-496刘飞. 气相色谱法同时测定食品中34种脂肪酸J. 福建分析测定， 2010,19(3):15-18LIU Fei. Simultaneous determination of 34 kinds

48、 of fatty acids in food by gas chromatoraphyJ. Fujian Analysis&Testing, 2010,19(3):15-187何艳华. HYPERLINK :8080/qk/90732X/200809/28312231.html t _blank 应用气相色谱法方法测定花生油掺伪的探讨J. 广西质量监督导报，2008,(9):71-74HE Yanhua. Determination of Adulteration of Peanut Oil by Gas ChromatographyJ.Guangxi Quality Supervision

49、 Guide Periodical，2008,(9):71-748喻凤香. 食用油脂气相色谱指纹图谱分析J.食品研究与开发，2016,2(37):120-123YU Fengxiang. Gas chromatographic fingerprint analysis of edible oils and waste oilsJ. Food Research and Development,2016,2(37):120-1239李艳青，于小玲，宋文静等. 胎儿体重相关因素分析及多元回归方程的建立J.现代生物医学进展，2012,12(25):4896-4898LI Yanqin，YU Xiaol

50、in，SONG Wenjing,et al. The analysis of factors about the weight of fetal and establishment of multiple regression equationJ.Modern Biomedical Progress,2012,12(25):4896-489810蔡冰，陈钢，温露等. 鸦胆子油微囊中油酸与亚油酸含量的气相色谱法测定J.时珍国医国药，2011,5(22):1057-1058CAI Bing，CHEN Gang，WENG Lu,et al.Determination of oleic acid and linoleic acid i