A-麦角甾醇；B-胆甾醇；C-菜油甾醇； D-豆甾醇；E-β-谷甾醇
图1 5种植物甾醇的吸收光谱图(扫描波长190～640 nm)
Fig.1 Adsorption chromatography of 5 phytosterols
红糖产品(红糖，黑糖，赤砂糖)因其保留了甘蔗原有的风味和营养物质(蛋白质，矿物质，维生素，酚类物质，有机酸，植物甾醇等非糖成分)，符合天然健康的饮食观念[1-2]，享有健脾养胃、温中补气、缓解精通的温补美誉，深受消费者喜爱。国内外学者研究发现，其保健功能物质基础可能与红糖产品中的非糖成分有密切关系[3-9]。甘蔗甾醇含量高，抗炎、降血脂、缓解非酒精性脂肪肝化能力强[10-13]，红糖作为甘蔗全汁浓缩物，甘蔗甾醇在传统提取浓缩工艺中自然进入红糖产品，但目前对传统工艺红糖产品的甾醇类成分研究报道较少。为探讨传统红糖的保健功能及机制，了解红糖产品中甾醇的含量及其在加工过程中的变化,本文建立了一种高效液相色谱法(high performance liquid chromatography，HPLC)测定红糖产品中5种植物甾醇含量的方法，并对红糖产品中3种常见的甾醇含量进行了分析测定,以期为红糖产品的深入开发提供理论依据。
高效液相色谱仪(Agilent 1290)，安捷伦科技(中国)有限公司；旋转蒸发仪(IKA RV8)、旋涡混合器(VORTEX 3)，艾卡(广州)仪器设备有限公司；恒温水浴锅(DKZ-2B)，上海一恒科技有限公司；医用离心机(CL5R)，湖南湘仪实验室仪器开发有限公司；Milli-Q超纯水系统，默克密理博公司。
甲醇、乙腈(色谱纯)，西格玛奥德里奇(上海)贸易有限公司；石油醚、二氯甲烷、乙醚、乙酸乙酯、KOH，无水乙醇(分析纯)，广州化学试剂有限公司。
色谱柱 Waters ACQUITY UPLC BEH C18(2.1 mm×100 mm，1.7 μm)；Waters ACQUITY UPLC BEH C8(2.1 mm×100 mm，1.7 μm)；ZORBAX SB C8 RRHD(2.1 mm×50 mm，1.8 μm)
标准物质：麦角甾醇、胆甾醇、豆甾醇、β谷甾醇和菜油甾醇(纯度≥95%)，上海安谱实验科技股份有限公司。
96批次食糖样品均来自2017～2019年广东、广西、云南3个不同产地，其中41个红糖样品，11个黑糖样品，44个赤砂糖样品。在分析之前，所有食糖样品都储存在塑料罐中，阴凉避光保存。
1.2.1 色谱条件：
Waters ACQUITY UPLC BEH C8(1.7 μm，2.1 mm×100 mm)液相色谱柱；检测波长205 nm(β-谷甾醇，菜油甾醇，胆甾醇和豆甾醇最大吸收波长)，282 nm(麦角甾醇最大吸收波长)；流速0.4 mL/min；柱温35 ℃；进样量5 μL。流动相：A相为水，B相为乙腈；采用梯度洗脱，0～9.5 min 80%～90% B相。
1.2.2 供试品溶液的制备
取约3.0 g样品置于锥形瓶中，加入无水乙醇溶液旋涡提取30 min，以体积比1∶1加入50%(质量分数)NaOH溶液，摇匀，在80 ℃水浴中皂化60 min，冷却后加入水10 mL，并移入分液漏斗中，用15 mL二氯甲烷分3次萃取，合并萃取液，旋干，用甲醇2 mL复溶，用0.22 μm微孔滤膜过滤，用于HPLC进样。
1.2.3 标准溶液的配制
分别称取适量麦角甾醇、胆甾醇、豆甾醇、β谷甾醇和菜油甾醇，置于10 mL容量瓶中，用甲醇溶解并定容至刻度，分别制成约1.0 mg/mL的储备液，于4 ℃下避光储存备用。取适量部分上述标准储备液用甲醇配制成质量浓度为0.1～100 ng/mL系列标准工作液。
2.1.1 检测波长选择
在190～400 nm下分别对5种甾醇标准溶液进行紫外扫描，发现麦角甾醇的最大吸收波长为282 nm，胆甾醇，豆甾醇，菜油甾醇和β-谷甾醇在205 nm处有最强吸收(见图1)，因此麦角甾醇的检测波长选择282 nm，胆甾醇，豆甾醇，菜油甾醇和β-谷甾醇的检测波长选择205 nm。
A-麦角甾醇；B-胆甾醇；C-菜油甾醇； D-豆甾醇；E-β-谷甾醇
图1 5种植物甾醇的吸收光谱图(扫描波长190～640 nm)
Fig.1 Adsorption chromatography of 5 phytosterols
2.1.2 色谱柱选择
尝试了安捷伦ZORBAX SB C8 RRHD(2.1 mm×50 mm，1.8 μm)，Waters ACQUITY UPLC BEH C8(2.1 mm×100 mm，1.7 μm)，Waters ACQUITY UPLC BEH C18(2.1 mm×100 mm，1.7 μm)三款色谱柱。发现C18对4种成分的分离效果好，但是运行比较长，在超高效色谱系统需要28 min才能完全出峰。不同牌子的C8色谱柱对于分离情况也有区别，相同流动相，相同梯度下ZORBAX SB C8 RRHD柱不能分离菜油甾醇和豆甾醇，Waters ACQUITY UPLC BEH C8则能实现8 min内5种成分完全分离。
2.1.3 流动相选择
流动相：甲醇截止波长为205 nm;乙腈截止波长为190 nm，而胆甾醇，豆甾醇，菜油甾醇和β-谷甾醇的最大吸收波长为205 nm，在甲醇水系统下，基线漂移严重。故选择乙腈水作为流动相。尝试了乙腈-水100∶0、98∶2、95∶5、90∶10、80∶20(体积比)等多种等度洗脱，均无法实现菜油甾醇和豆甾醇的分离。梯度洗脱：0～9.5 min乙腈80%～90%，则能5种成分完全基线分离；流速，柱温对峰位和组分分离影响不大。在此条件下菜籽甾醇、豆甾醇、菜油甾醇、β-谷甾醇分离情况良好，满足定量分析要求，见图2、图3。
A-205 nm；B-282 nm；1-麦角甾醇;2-胆甾醇;3-菜油甾醇;4-豆甾醇;5-β-谷甾醇
图2 标准色谱图
Fig.2 Chromatogram of standard
A-205 nm；B-282 nm；1-麦角甾醇;2-胆甾醇;3-菜油甾醇;4-豆甾醇;5-β-谷甾醇
图3 样品色谱图
Fig.3 Chromatogram of sample
2.1.4 提取方法的优化
文献报道的提取方法大多数使用回流皂化，溶剂萃取再浓缩，耗费时间比较长，操作复杂[10-13]。也有针对样品基质比较简单的，仅需超声或涡旋提取，皂化即得[14-20]。因为糖品的基质较简单，所以本文也采用直接提取皂化法。前期在文献总结的基础上运用响应面法对超声提取时间、料液比、皂化时间和萃取次数进行实验探究和优化，得到5种甾醇的提取工艺条件为：料液比为1∶10，提取时间为19.58 min，皂化时间1 h，萃取次数3次，提取总量预测为20.012 μg。对预测的最佳工艺条件进行了验证并结合实际操作，将验证条件修正为料液比为1∶10，提取时间为20 min，皂化时间1 h，萃取次数3次。
2.2.1 线性范围
精密吸取混合对照品储备液一定体积，按倍数关系稀释成不同质量浓度的标准溶液进行测定，采用优化后的前处理方法以及色谱条件，以5种甾醇色谱峰的峰面积平均值(Y)为纵坐标，对应的标准溶液浓度(X)为横坐标，绘制标准曲线；以信噪比S/N=3确定仪器检出限。结果如表1所示，5种植物甾醇在各自浓度范围内线性关系良好，相关系数均不小于0.999 7；麦角甾醇检出限均为0.023 mg/kg，胆甾醇为0.14 mg/kg，菜油甾醇为0.18 mg/kg，豆甾醇为0.11 mg/kg，β-谷甾醇为0.17 mg/kg，方法灵敏度较高，可满足糖品中甾醇检测需要。
表1 方法的线性范围、回归方程及检出限
Table 1 Linear range，linear equations and limits of detection of the method
成分线性方程R2线性范围/(μg·mL-1)检出限/(mg·kg-1)定量限/(mg·kg-1)精密度(相对标准偏差)/%麦角甾醇Y=129.57X-20.620.999 900.1～100.00.0230.0681.6胆甾醇 Y=75.49X+5.070.999 870.1～100.00.140.433.9菜油甾醇Y=58.21X-10.710.999 830.1～100.00.180.541.3豆甾醇 Y=84.51X+36.860.999 720.1～100.00.110.340.20β-谷甾醇Y=60.57X+0.990.999 900.1～100.00.170.520.29
2.2.2 精密度试验
取供试品溶液按照1.2.1项下色谱条件连续进样6次，结果显示，麦角甾醇、胆甾醇、菜油甾醇、豆甾醇、β-谷甾醇峰面积的精密度(相对标准偏差)分别为1.6%、3.9%、1.3%、0.20%、0.29%表明仪器的精密度良好。
2.2.3 加标回收率
精确称取红糖样品3.0 g，按重复性结果计算得到5个甾醇成分的含量，分别加入3个水平的对照品，在最佳提取条件及色谱条件下，每一浓度平行3次测定红糖中植物甾醇的回收率，结果见表2。由表2可知，红糖中麦角甾醇的平均回收率为74.4%～95.8%，胆甾醇的平均回收率为82.3%～91.4%，菜油甾醇的平均回收率为84.6%～102.5%，豆甾醇的回收率为88.0%～92.1%，β-谷甾醇的回收率为75.4%～92.6%。
表2 方法的回收率(n=6)
Table 2 Recovery of the method(n=6)
加标量/(mg·kg-1)测定均值/(mg·kg-1)平均回收率/%精密度(相对标准偏差)/%0.200.1679.62.5 麦角甾醇1.000.75 74.4 0.7 2.001.8995.83.610.008.7782.33.1 胆甾醇50.0042.28 87.71.9 100.0091.3991.4 1.11010.25 102.5 2.7 菜油甾醇5037.3384.61.1 10091.1891.20.8109.1591.52.3 豆甾醇5044.088.00.2 10092.0892.10.1107.5475.40.2 β-谷甾醇5044.99 90.00.3 10092.60 92.60.1
2.3.1 样品含量测定
取96批不同产地和级别的红糖，黑糖和赤砂糖样品按优化后的最佳提取条件和最优色谱条件进行含量测定，按2.4.1中的线性方程计算麦角甾醇、胆甾醇、菜油甾醇、豆甾醇、β-谷甾醇的含量，结果见表3。由表3可知，所测96个样品中大部分含有5种植物甾醇。5种植物甾醇以豆甾醇和β-谷甾醇含量较高，可达244.41 mg/kg、67.38 mg/kg，麦角甾醇含量最低。
3种糖品中平均总甾醇含量差异明显，分布为红糖>黑糖>赤砂糖。GB/T 35886—2018[21]食糖分类中提到，红糖是糖汁清净处理后直接炼煮不经分蜜制成。而赤砂糖则是经过了分蜜制成，是工业化生产白砂糖的副产品。目前还没有黑糖的国家标准，据DB 4413/T 10—2019[22]和QB/T 4567—2013[23]描述其生产工艺与红糖相似，只是炼煮的颜色差异。因此总甾醇含量差异可能与3种糖品的生产工艺或产地有关，但还需进一步分析研究。
从糖品的等级上看，无论红糖的优级、一级和二级还是赤砂糖的一、二级，在总甾醇含量上并没有多大区分。
表3 实际样品的测定结果 单位：mg/kg
Table 3 Results of sample determination for phytosterols content
样品麦角甾醇胆甾醇菜油甾醇豆甾醇β-谷甾醇总量等级产地样品麦角甾醇胆甾醇菜油甾醇豆甾醇β-谷甾醇总量等级产地S10.445.077.4313.5012.3538.80优级红糖广西崇左S20.22ND4.095.604.8514.76优级红糖广西崇左S30.243.2016.1421.4729.2970.34一级红糖云南镇康S40.28ND11.8216.0219.5847.70一级红糖云南昌宁S50.580.4221.3531.3835.2989.01一级红糖广东英德S616.341.536.349.2410.4043.85一级红糖广西柳城S70.608.489.3513.3517.7649.54一级红糖广东博罗S80.593.8820.0625.3232.9682.82一级红糖广西马山S90.621.9722.6034.2339.8399.26一级红糖广东博罗S100.429.8918.58244.4132.58305.89一级红糖云南昌宁S110.821.5018.0125.5325.1270.98一级红糖广西金秀S120.831.4820.3627.0933.1282.88一级红糖广西马山S130.481.1921.7827.1933.0783.70一级红糖广西马山S1416.231.676.639.1512.8846.56一级红糖广西柳城S150.850.1919.7726.2929.1176.22一级红糖广东英德S160.945.339.8517.0018.0151.12一级红糖广东江门S170.390.8219.7725.6431.2677.86一级红糖云南昌宁S180.153.321.060.552.697.77一级红糖广西农垦S190.120.900.840.862.645.36一级红糖云南镇康S200.300.412.753.123.5810.16一级红糖广西凤山S210.225.475.488.199.5428.90一级红糖广东博罗S220.16ND4.626.237.5218.52一级红糖云南昌宁S230.171.054.336.286.2518.08一级红糖广西农垦S240.121.085.577.449.2523.45一级红糖云南红河S250.28/6.298.909.5425.01一级红糖广东江门S260.291.3311.3416.2719.2148.45一级红糖广东博罗S270.13/4.215.717.2617.31一级红糖云南昌宁S280.351.6010.2916.0916.2344.57一级红糖广西柳城S290.17ND9.5415.4922.6247.82一级红糖云南镇康S300.28ND11.4518.5923.2253.53一级红糖广西凤山S31ND0.731.060.401.423.61一级红糖广西农垦S320.260.1910.8017.2516.3144.80二级红糖云南昆明S330.120.291.190.491.653.74二级红糖广西雒容S340.51ND13.2218.8619.4452.02二级红糖广东江门S35ND0.310.530.381.642.86二级红糖广西柳州S36ND1.590.590.141.293.61二级红糖广西雒容S370.161.373.004.594.6313.76二级红糖广东博罗S380.990.7421.0625.9738.3087.06二级红糖广西马山S390.704.8518.5126.2334.3884.67二级红糖广东佛山S400.85ND21.3228.3640.3290.86二级红糖广东佛山S410.43ND31.8942.9461.45136.72二级红糖云南镇康平均1.161.7511.0920.7718.9753.75S420.624.9420.4129.5639.4594.98一级黑糖广东博罗S431.02ND39.3159.5267.38167.24一级黑糖广东博罗S440.350.977.6114.1113.9336.97一级黑糖广西凤山S450.811.655.048.648.6524.78一级黑糖广西凤山S460.248.702.202.578.9622.68一级黑糖广西雒容S470.082.181.861.293.749.15一级黑糖广西雒容S480.99ND37.0656.8056.80151.66一级黑糖广东博罗S490.904.8016.7624.1424.1470.73一级黑糖广东博罗S500.442.378.5414.3614.3640.08一级黑糖广西柳城S510.351.518.2411.9211.9233.94一级黑糖广西凤山S520.096.99ND12.0112.0131.10一级黑糖广西雒容平均0.543.1013.3721.3623.7662.12S530.10/0.511.041.042.68一级赤砂糖广西环江S54ND2.96NDNDND2.96一级赤砂糖福建漳州S550.091.550.570.13/2.34一级赤砂糖广西凤山S56ND1.030.600.980.983.59一级赤砂糖广西香山S570.43ND13.6123.8023.8061.64一级赤砂糖广西农垦S580.27/7.8016.4016.4040.87一级赤砂糖广西崇左S59ND1.45NDNDND1.45一级赤砂糖广西上思S60ND1.001.061.401.404.85一级赤砂糖广西贵港S610.140.552.454.164.1611.45一级赤砂糖广西廖平S620.130.230.791.251.253.67一级赤砂糖广西北海S630.141.34ND5.125.1211.72一级赤砂糖广西上思S64ND0.390.530.190.191.31一级赤砂糖广西伶俐S650.1018.32ND0.11/18.54一级赤砂糖广西贵港S660.403.099.2116.8716.8746.43一级赤砂糖广西崇左S67ND0.29NDNDND0.29一级赤砂糖广西明阳S680.291.817.0413.4713.4736.08一级赤砂糖广西崇左S69ND3.192.561.301.308.34一级赤砂糖广西柳城S701.600.900.550.370.373.80一级赤砂糖广西凤山S710.25ND3.666.156.1516.22一级赤砂糖广西浦北S720.08ND1.231.761.764.83一级赤砂糖广西融水S730.243.081.100.16/4.58一级赤砂糖广西香山S74ND3.23NDNDND3.23一级赤砂糖广西明阳S750.090.950.921.441.444.85一级赤砂糖广西柳城S760.07NDNDNDND0.07一级赤砂糖广西上思S770.140.331.191.301.304.25一级赤砂糖广东遂溪S78ND/NDNDNDND一级赤砂糖福建漳州S790.250.511.320.990.994.04一级赤砂糖广西灵山S80ND8.72NDNDND8.72一级赤砂糖广西农垦S810.167.242.112.232.2313.96一级赤砂糖广西廖平S820.08ND0.560.470.471.59一级赤砂糖广西柳城S830.251.236.268.688.6825.09一级赤砂糖广西崇左S840.072.430.870.420.424.22一级赤砂糖广西北海S850.120.70NDNDND0.81一级赤砂糖广西农垦S86ND0.63NDNDND0.63一级赤砂糖广西农垦S870.160.183.024.274.2711.90二级赤砂糖广西上林S88ND0.71NDNDND0.71二级赤砂糖广西融安S890.220.582.211.671.676.35二级赤砂糖广西上林S900.141.662.091.791.797.47二级赤砂糖广西上林S910.08/0.730.190.191.18二级赤砂糖广西雒容S920.21ND3.215.715.7114.83二级赤砂糖广西柳江S930.262.234.396.316.3119.51二级赤砂糖广西柳江S940.09NDNDNDND0.09二级赤砂糖广西雒容S950.10ND2.862.842.848.64二级赤砂糖广西六塘S960.103.99NDNDND4.10二级赤砂糖广西融安平均0.161.741.933.023.019.86
注：“/”表示低于检出限；ND表示未检出
2.3.2 聚类分析
综合不同产地和级别的糖品中麦角甾醇、胆甾醇、菜油甾醇、豆甾醇、β-谷甾醇的含量和甾醇总含量这6个指标，借助SPSS软件并采用系统聚类法中的Ward 联结和平方Euclidean距离系数对红糖，黑糖，赤砂糖样品进行聚类分析，结果如图4～图6所示。
图4 红糖聚类分析树状图
Fig.4 Brown sugar cluster analysis dendrogram
图5 黑糖聚类分析树状图
Fig.5 Dark brown sugarcluster analysis dendrogram
图6 赤砂糖聚类分析树状图
Fig.6 Brown granulated suga rcluster analysis dendrogram
由图4可知，41个红糖样品可分为四大类。第一类为2、18～25、27、31、33、35～37号样品，第二类为1、4、6、7、14、16、26、28、29、30、32、34号样品，第三类为3、5、8、9、11～13、15、17、38～41号样品；10号样品与其他样品间距较大，分一类。
由图5可知，11个黑糖样品可分为三大类。第一类为44～47、50～52号样品，第二类为42、49号样品，第三类为43、48号样品；
由图6可知，44个赤砂糖样品可分为五大类。第一类为53～56、59～60、62、64、67、70、72～79、82、84～86、88、91、94、96号样品，第二类为71、83、92～93号样品，第三类为61、63、65、69、80、81、87、89、90、95号样品；第四类为58、66、68号样品；57号样品与其他样品间距较大，分为一类。
本文通过对3种糖品的前处理和色谱条件的优化，建立了超高效液相色谱法同时测定红糖中麦角甾醇、胆甾醇、菜油甾醇、豆甾醇、β-谷甾醇这5种植物甾醇的方法。对3种糖品96个样本中5种植物甾醇的含量进行检测，结果发现其中5种植物甾醇以豆甾醇和β-谷甾醇含量较高，麦角甾醇含量最低，不同糖品中总甾醇成分含量差异较大，分布为红糖>黑糖>赤砂糖。植物甾醇作为具有一定药用价值的成分，研究其在红糖产品中的含量可为其产品质量标准的提高提供参考，为其合理开发利用奠定基础。
[1] 黄苏婷. 红糖香气成分鉴别及其溯源初探[D].南宁：广西大学,2018.
[2] 苏佳廷,李凯,杭方学,等.原生态红糖[J].广西糖业,2014(2):37-40.
[3] YOSHIMOTO M, KURATA R, FUJII M, et al. In vitro and in vivo anticarcinogenesis of sugar cane vinegar[J]. Acta Horticulturae, 2008，765: 17-22.
[4] ASIKIN Y, CHINEN T, TAKARA K, et al. Determination of long-chain alcohol and aldehyde contents in the non-centrifuged cane sugar kokuto[J]. Food Science and Technology Research, 2008, 14(6): 583-588.
[5] 甘蔗脂肪醇、红曲提取物和虾青素对家兔的抗动脉粥样硬化作用[J]. 国外医药(植物药分册), 2006, 21(4): 172-173.
[6] SINGH N, KUMAR D, LAL K, et al. Adverse health effects due to arsenic exposure: Modification by dietary supplementation of jaggery in mice[J]. Toxicology and Applied Pharmacology, 2010, 242(3): 247-255.
[7] SUMIYOSHI M,HAYASHI T,KIMURA Y.Effects of the nonsugar fraction of brown sugar on chronic ultraviolet B irradiation-induced photoaging in melanin-possessing hairless mice[J].Natural medicines，2009(2):130-136.
[8] 邱荣其, 赵振刚, 丰永红, 等. 甘蔗抗氧化物质的提取与分离[J]. 中国甜菜糖业, 2007(3): 1-4.
[9] TAKARA K, OTSUKA K, WADA K, et al. 1,1-diphenyl-2-picrylhydrazyl radical scavenging activity and tyrosinase inhibitory effects of constituents of sugarcane molasses[J]. Bioscience Biotechnology &Biochemistry, 2007, 71(1):183-191.
[10] 冯思敏. 甘蔗中植物化学素组成及其中甾醇功能活性和相关机理的研究[D].杭州：浙江大学,2017.
[11] 冯思敏. 直接酸碱皂化法处理应用于甘蔗中自由及结合态植物甾醇含量测定的研究[C]. 中国食品科学技术学会第十一届年会论文摘要集.北京：中国食品科学技术学会,2014:402-403.
[12] 武晓云. 甘蔗糖厂滤泥中植物甾醇的提取及豆甾醇和β-谷甾醇的分离[D].南宁：广西大学,2007.
[13] 武晓云,郭海蓉,陈赶林,等.从甘蔗滤泥中提取植物甾醇的研究[J].食品研究与开发,2007(2):64-68.
[14] 项雷文,陈国泰,翁佳敏,等.HPLC-DAD结合交替三线性分解二阶校正法测定荔枝果肉中游离植物甾醇[J].食品科学,2015,36(24):172-176.
[15] 冯文欢,吴正章,张鹏,等.HPLC-UV法测定植物甾醇中3种甾醇单体的含量[J].食品工业科技,2019,40(5):244-248.
[16] 林宏琳.HPLC法测定植物甾醇软胶囊中植物甾醇的含量[J].营养学报,2018,40(4):398-402.
[17] 赵海燕,詹未,李兵,等.超高效液相色谱法同时测定食用植物油中的4种植物甾醇[J].中国食品卫生杂志,2016,28(2):223-226.
[18] 张冬阳. 薏米中β-谷甾醇的提取、分离、纯化及其功能性的研究[D].锦州:锦州医科大学,2018.
[19] 江海.HPLC测定食用菜籽油中β-谷甾醇的含量[J].汉中师范学院学报(自然科学),2004(6):46-49.
[20] 刘晓飞,宋洁,王薇,等.发芽糙米植物甾醇的提取优化及抗氧化性研究[J].哈尔滨商业大学学报(自然科学版),2019,35(1):44-48.
[21] 国家标准化管理委员会.GB/T 35886—2018 食糖分类[S].北京:中国标准出版社, 2018.
[22] DB4413/T 10—2019.惠州市市场监督管理局地理标志产品 观音阁红糖(观音阁黑糖)[S].
[23] 国家工业和信息化部. QB/T 4567—2013黑糖[S]. 北京：中国轻工业出版社，2013.
[24] 葛志伟,罗自生.食品中植物甾醇的色谱分析方法研究进展[J].食品安全质量检测学报,2017,8(3):773-778.