随着石油工业的发展,人类活动造成的突发性石油泄露事故时常发生,如泉州市泉港区东港石化公司码头发生裂解碳九外泄事故[1]。一旦发生溢油事故,会造成生态环境如水体、土壤、大气等的污染,污染物进入到生态循环中,对生物圈造成严重的破坏,同时人类食用被污染过的水产品等,会有中毒、致癌等风险[2-5],因此对水体、水产品、土壤等进行检测显得尤为重要,如多环芳烃[6]、裂解碳九等的检测可以有效预防和评估石油的污染。
裂解碳九是乙烯装置的副产物,是由裂解石脑油经抽提分离出C5馏分、C6~C8馏分(经加氢生产BTX)后剩余的馏分[7],约占乙烯总产量的10%~20%[8-9]。碳九芳烃是裂解碳九的重要组成部分,它是多组分的混合物[10],其馏分组分复杂,组成多达数10种以上[11-14],其中苯系物组成较多,主要成分有双茂、茚、苯乙烯、二甲苯、甲苯、乙苯、萘、三甲苯、双环戊二烯、异丙苯、正丙苯等[15-16]。为了更好地分辨碳九芳烃,需要选择一种合适的检测手段,这对于石油工业的发展具有重要意义。目前国内对碳九芳烃检测的标准和相关研究非常少,而国外对碳九芳烃检测的研究迄今未见报道, SH/T 1793—2015《工业用裂解碳九组成的测定 气相色谱法》[17]对裂解碳九的组成分析做了规范,肖洁等[18]采用HS-GC/MS法对裂解碳九进行分析检测;孙桂芳等[15]采用GC-MS与气相色谱-傅立叶变换红外光谱联用技术法联合测定C9烃馏分组成。本文通过三重四级杆气质联用法(GC-MS/MS)来测定水产品中碳九芳烃的含量,该法不仅具有操作简便、有机物分离效果良好、灵敏度高、重现性和准确度较好、有机溶剂消耗少等优点,而且能够分辨出气相色谱法和气相质谱法所不能分辨的同分异构体[19-22],使更多种类的碳九芳烃的定性定量成为可能。该法的建立为地标乃至国标的实施提供了理论参考依据,填补了食品检测、环境检测、石油工业等行业中碳九芳烃检测市场的空白。
GCMS-TQ8040 NX 三重四级杆气质联用仪,日本岛津公司;GM200 型刀式研磨机,德国莱驰公司;HFJ-25 型匀浆机,天津市恒奥发展有限公司;RGL-200T 型高速冷冻离心机,上海卢湘仪离心机仪器有限公司;MTV-100 型涡旋混合器,东莞市谱标实验器材科技有限公司;JJ224BF 型电子天平,常州市双杰测试仪器厂;FYL-Y5-100L 型恒温冰箱,思锦科技有限公司;PS-100A 型超声清洗机,东莞市谱标实验器材科技有限公司。
双环戊二烯(97.9%)、α-甲基苯乙烯(99.0%)、苯乙烯(99.9%)、萘(99.8%)、茚(99.1%)、1,3,5-三甲苯(99.8%)、异丙苯(99.9%)、正丙苯(99.7%)、邻二甲苯(99.4%)、间二甲苯(99.2%)、乙苯(99.5%)、对二甲苯(99.0%),北京坛墨质检科技有限公司;重蒸正己烷、重蒸丙酮、乙酸乙酯(均为色谱纯);乙二胺-N-丙基硅烷(primary secondary amine,PSA)(填料级);中性氧化铝、无水硫酸镁(均为分析纯)。
用正己烷配制12种碳九芳烃的混合标准储备液,质量浓度为100 μg/mL,并置于4 ℃冰箱避光保存。
用乙酸乙酯将碳九混合标准储备液稀释成质量浓度为10 μg/mL的标准中间溶液,置于4 ℃冰箱避光保存。
碳九标准工作溶液:吸取一定量的标准中间溶液,用乙酸乙酯稀释成适当浓度的标准工作溶液(质量浓度为0.01~1.0 μg/mL)。
将水产品(市场购置)可食部分用研磨机充分捣碎制成试样,称取5.0 g试样(精确至 0.01 g)于50 mL具塞离心管中,加入10 mL乙酸乙酯,均质2 min,-20 ℃下冷冻6 h,于4 000 r/min冷冻离心5 min,取上清液移入50 mL具塞离心管中,加入0.1 g 乙二胺-N-丙基硅烷(primary secondary amine,PSA)、1 g中性氧化铝、1 g无水硫酸镁,涡旋3 min,于4 000 r/min离心5 min,取上清液过0.22 μm微孔滤膜,供GC-MS/MS测定。
1.4.1 气相色谱条件
色谱柱:SH-Rtx-Wax柱(30 m×0.25 mm×0.25 μm);升温程序:初始温度40 ℃,保持3.5 min;然后以5 ℃/min的速率升温至60 ℃,保持3 min;再以20 ℃/min的速率升温至180 ℃,保持2.5 min;最后以20 ℃/min的速率升温至250 ℃,保持1 min。进样口温度220 ℃,色谱柱温度40 ℃,进样量1.0 μL,进样时间1.00 min,进样方式:不分流进样。
1.4.2 质谱条件
离子化方式:电子轰击离子源(electron impact ion source,EI);离子源温度230 ℃;接口温度220 ℃;检测器电压:相对于调谐电压0.2 kV;扫描模式:质谱多反应监测(multiple reaction monitoring,MRM);溶剂延迟时间3 min。12种碳九芳烃优化后的保留时间和MRM质谱参数见表1。
表1 12种碳九芳烃的保留时间和MRM质谱参数
Table 1 Retention time and MRM mass parameters of 12 kinds of C9 aromatics
序号化合物保留时间/min离子对(m/z)碰撞电压/eV91.00>65.00∗131乙苯6.663106.00>91.1013106.00>65.0028106.00>91.10∗132对二甲苯6.81091.00>65.1018106.00>65.1028106.00>91.00∗133间二甲苯6.97291.00>65.1018106.00>65.0028120.00>105.10∗84异丙苯7.888105.00>77.1018105.00>79.108106.00>91.10∗135邻二甲苯8.12091.00>65.1013106.00>65.102891.00>65.00∗186正丙苯8.981120.00>91.0013120.00>65.1028120.00>105.10∗1371, 3, 5-三甲苯10.256105.00>77.1018105.00>79.108132.00>117.00∗88双环戊二烯10.387132.00>66.00866.00>51.0023104.00>77.90∗139苯乙烯10.80878.00>52.0013104.00>51.0028103.00>77.10∗1310α-甲基苯乙烯12.586118.00>91.1023118.00>77.0028116.00>114.20∗1811茚14.540116.00>89.0028116.00>65.1028128.00>102.10∗1812萘16.715128.00>78.1023128.00>76.0028
注:“*”,定量离子对
首先确定气相色谱和质谱的各个参数,然后创建Q3 SCAN全扫描方法,用Q3 SCAN全扫描方式测定12种碳九芳烃的混合标准溶液,得到其总离子流图(total ion chromatography,TIC)。根据TIC,选择m/z比和相对丰度较高的一级碎片为不同碳九芳烃的前体离子,并得到其保留时间。进一步根据已确定的前体离子,在3~45 V范围内间隔5 V设定不同的碰撞电压,进行产物离子扫描模式测定,通过MRM-Optimization-Tool表自动优化碰撞电压和离子对,最后利用Smart Database数据库建立12种碳九芳烃的MRM方法。12种碳九芳烃的TIC见图1,其质谱参数见表1。
图1 十二种碳九芳烃的总离子流图
Fig.1 TIC of 12 kinds of C9 aromatics
对乙酸乙酯、正己烷、丙酮和正己烷∶丙酮(体积比7∶3)混合溶液共4种有机溶剂提取水产品中目标物的能力进行试验,如图2所示,在加标质量分数为0.2 mg/kg下,采用乙酸乙酯为提取溶剂的加标回收率在85%~110%,正己烷和正己烷∶丙酮(体积比7∶3)混合溶液回收率基本>120%,丙酮回收率基本<70%。这4种提取溶剂均能溶解碳九芳烃,但正己烷易挥发,使目标物的浓度增大,回收率变高,丙酮极性大,能更多地提取试样中的蛋白质等杂质,而对碳九芳烃的溶解度较差,导致其提取效果不如乙酸乙酯和正己烷,回收率较低,正己烷∶丙酮(体积比7∶3)混合溶液中正己烷的影响占主导,其回收率偏高,而乙酸乙酯相对不易挥发,且提取杂质少,其回收率满足碳九芳烃提取溶剂的要求。综上分析,选用乙酸乙酯作为最终的提取溶剂。
a-基质:鱿鱼;b-基质:姑鱼;c-基质:对虾
图2 不同提取溶剂对12种碳九芳烃回收率的影响
Fig.2 Effects of different extraction solvents on the recovery of 12 kinds of C9 aromatics
分别采用超声和均质这2种提取方式,经过不同的提取时间,分析其对12种碳九芳烃回收率的影响,如图3所示。在加标质量分数为0.2 mg/kg下,超声5、15和30 min的回收率分别为68%~76%、73%~80%和81%~87%,均质2 min的回收率为82%~89%,回收率均较高,可见在一定范围内,随着超声时间的增加,其回收率逐渐增加,均质2 min和超声30 min的回收率基本相同,但所用时间明显更少。均质提取可以将试样充分打碎,起到降低试样尺度和提高目标物分布均匀性的作用,从而将目标物从试样中提取出来,该提取方法操作简单,效率高,可以充分地提取目标物。因此,选择均质提取作为本试验的提取方式,提取时间为2 min。
图3 不同提取方式和提取时间对12种碳九芳烃 回收率的影响
Fig.3 Effects of different extraction methods and extraction time on the recovery of 12 kinds of C9 aromatics
为了研究不同净化吸附剂用量对碳九芳烃回收率的影响,设定加入0.1 g PSA、1 g无水硫酸镁、1 g中性氧化铝,在维持另外2种吸附剂用量不变的条件下,分别调整PSA、无水硫酸镁和中性氧化铝的用量,得各吸附剂用量与12种碳九芳烃回收率的关系柱状图,如图4所示。在一定范围内,12种碳九芳烃回收率与中性氧化铝、无水硫酸镁、PSA的用量分别呈负相关、正相关、正相关,但回收率均在85%~120%,且各回收率之间相差不大。综合考虑目标物质的回收率、净化效果和耗材成本,折中选择加入的净化吸附剂用量为:0.1 g PSA、1 g无水硫酸镁、1 g中性氧化铝。
a-中性氧化铝;b-无水硫酸镁;c-PSA
图4 不同净化吸附剂用量对12种碳九芳烃回收率的影响
Fig.4 Effects of different dosage of purification adsorbent on the recovery of 12 kinds of C9 aromatics
取大黄鱼空白试样5.0 g,配制成质量浓度分别为0.01、0.05、0.1、0.5和1 μg/mL的基质标准工作溶液,按照“1.3”和“1.4”所示方法进行测定,以质量浓度为横坐标(x)、峰面积为纵坐标(y)绘制标准工作曲线。12种碳九芳烃在质量浓度为0.01~1.0 μg/mL均呈良好的线性关系,相关系数均>0.999,再以3倍和10倍信噪比(S/N)为基准分别得到仪器检出限为0.124~2.502 μg/L、定量限为0.413~8.340 μg/L,再计算出方法检出限为0.248~5.004 μg/kg,如表2所示。
为验证该方法的精密度,选取3种空白水产品:大黄鱼、对虾、鱿鱼,分别加入相对于曲线最低浓度1、5、10和50倍浓度(即0.02,0.1,0.2和1 mg/kg)的4种质量分数的碳九芳烃混合标准溶液进行添加回收试验,按照“1.3”和“1.4”进行测定,每一浓度水平测定6次(n=6),计算得到这3种水产品中12种碳九芳烃的平均回收率和相对标准偏差(relative standard deviation,RSD),如表3所示,平均回收率在75%~112%,RSD<6.7%。
表2 12种碳九芳烃的线性方程、相关系数(R2)、仪器检出限、定量限和方法检出限
Table 2 Linear equations, correlation coefficients(R2), limits of instrument detection, limits of quantitation and limits of method detection of 12 kinds of C9 aromatics
序号化合物线性方程相关系数(R2)仪器检出限/(μg·L-1)定量限/(μg·L-1)方法检出限/(μg·kg-1)1乙苯y=1.2×106x+4.4×1020.999 95 0.590 1.967 1.180 2对二甲苯y=1.3×106x+4.3×1030.999 84 0.429 1.430 0.858 3间二甲苯y=1.2×106x-2.6×1030.999 89 0.492 1.640 0.984 4异丙苯y=1.3×106x-2.4×1030.999 97 0.486 1.620 0.972 5邻二甲苯y=1.2×106x-1.5×1030.999 97 0.423 1.410 0.846 6正丙苯y=1.1×106x-1.7×1030.999 96 0.368 1.227 0.736 71, 3, 5-三甲苯y=9.9×105x-2.6×1030.999 88 0.727 2.425 1.455 8双环戊二烯y=2.2×105x-6.3×1020.999 90 2.502 8.340 5.004 9苯乙烯y=2.3×105x+7.3×1020.999 84 1.400 4.667 2.800 10α-甲基苯乙烯y=6.0×105x-2.3×1030.999 84 1.026 3.420 2.052 11茚y=1.6×106x-5.2×1030.999 89 0.309 1.030 0.618 12萘y=9.6×105x-2.1×1030.999 94 0.124 0.413 0.248
表3 水产品中12种碳九芳烃的平均回收率和RSD
Table 3 Average recoveries and RSD of 12 kinds of C9 aromatics in aquatic products
序号化合物试样加标质量分数/(mg·kg-1)回收率/%RSD/%(n=6)大黄鱼0.02, 0.1, 0.2, 184~921.5~4.91乙苯对虾0.02, 0.1, 0.2, 175~920.8~5.1鱿鱼0.02, 0.1, 0.2, 179~942.8~4.9大黄鱼0.02, 0.1, 0.2, 193~1041.3~3.92对二甲苯对虾0.02, 0.1, 0.2, 192~1011.3~5.8鱿鱼0.02, 0.1, 0.2, 194~1061.9~3.9大黄鱼0.02, 0.1, 0.2, 181~991.9~4.63间二甲苯对虾0.02, 0.1, 0.2, 176~961.5~6.7鱿鱼0.02, 0.1, 0.2, 181~990.9~5.5大黄鱼0.02, 0.1, 0.2, 190~970.8~3.74异丙苯对虾0.02, 0.1, 0.2, 191~1021.7~3.4鱿鱼0.02, 0.1, 0.2, 191~1041.5~3.9大黄鱼0.02, 0.1, 0.2, 189~1010.4~3.55邻二甲苯对虾0.02, 0.1, 0.2, 187~1010.8~4.2鱿鱼0.02, 0.1, 0.2, 184~1031.1~4.4大黄鱼0.02, 0.1, 0.2, 181~940.6~4.66正丙苯对虾0.02, 0.1, 0.2, 183~950.5~4.5鱿鱼0.02, 0.1, 0.2, 180~960.2~3.5大黄鱼0.02, 0.1, 0.2, 185~960.5~2.071, 3, 5-三甲苯对虾0.02, 0.1, 0.2, 192~951.0~6.7鱿鱼0.02, 0.1, 0.2, 191~931.7~3.9大黄鱼0.02, 0.1, 0.2, 190~1050.6~4.58双环戊二烯对虾0.02, 0.1, 0.2, 192~1101.0~2.2鱿鱼0.02, 0.1, 0.2, 194~1121.0~3.0大黄鱼0.02, 0.1, 0.2, 185~1011.0~4.39苯乙烯对虾0.02, 0.1, 0.2, 188~1041.1~6.4鱿鱼0.02, 0.1, 0.2, 181~1000.4~4.5大黄鱼0.02, 0.1, 0.2, 182~922.2~5.110α-甲基苯乙烯对虾0.02, 0.1, 0.2, 183~951.1~4.8鱿鱼0.02, 0.1, 0.2, 187~1101.2~2.5大黄鱼0.02, 0.1, 0.2, 181~1030.5~4.611茚对虾0.02, 0.1, 0.2, 181~1071.4~3.4鱿鱼0.02, 0.1, 0.2, 182~1080.3~2.9大黄鱼0.02, 0.1, 0.2, 184~1060.4~5.012萘对虾0.02, 0.1, 0.2, 195~1051.3~3.7鱿鱼0.02, 0.1, 0.2, 186~1031.3~4.4
采用“1.3”和“1.4”的方法检测了海鲈鱼、白鲳鱼、鲍鱼、对虾、罗非鱼、龙舌鱼、鳗鱼等市场上常见的10种水产品试样的碳九芳烃,结果显示只有龙舌鱼检测出碳九芳烃,具体结果见表4,检测数值为间二甲苯70.510 μg/kg、对二甲苯29.410 μg/kg、苯乙烯631.980 μg/kg、邻二甲苯32.650 μg/kg和萘47.990 μg/kg,造成苯乙烯含量较高的原因可能是龙舌鱼食用了含有较多苯乙烯污染物的食物,并不断在生物体中累积。由此可见,市场上的大部分水产品未被碳九污染,但仍有一部分受到一定程度的污染,其污染源可能来自于炼油厂泄露和油轮漏油,说明保护环境依然任重道远。
表4 龙舌鱼中12种碳九芳烃的检测结果
Table 4 Detection results of 12 kinds of C9 aromatics in sole fish
序号化合物含量/(μg·kg-1)1乙苯N.D2对二甲苯29.4103间二甲苯70.5104异丙苯N.D5邻二甲苯32.6506正丙苯N.D71, 3, 5-三甲苯N.D8双环戊二烯N.D9苯乙烯631.98010α-甲基苯乙烯N.D11茚N.D12萘47.990
注:“N.D”,未检出
本文建立了水产品中12种碳九芳烃的三重四级杆气质联用(GC-MS/MS)分析测定方法,该方法的线性、检出限、回收率、RSD等技术指标均满足检测方法的要求,能够准确地对水产品中的碳九芳烃进行定性定量检测,方法的建立为地标乃至国标的实施提供理论参考依据,填补了食品检测、环境检测、石油工业等行业中碳九芳烃检测市场的空白。
[1] 张宏, 周伯煌.泉港碳九泄漏事件引起的生态环境危机[J].中国环境管理干部学院学报,2019,29(1):38-41.
ZHANG H,ZHOU B H.Ecological environment crisis caused by carbon 9 leakage event in Quangang[J].Journal of the Environmental Management College of China-EMCC,2019,29(1):38-41.
[2] 朱林超, 王东光,顾佳鹏.海上石油的泄漏与防治[J].山东化工,2018,47(23):88-89.
ZHU L C,WANG D G,GU J P.Leakage and treatment methods of offshore oil[J].Shandong Chemical Industry,2018,47(23):88-89.
[3] 余春浩. 石油泄漏对土壤的污染[J].当代化工,2019,48(10):2 385-2 387;2 423.
YU C H.Review of soil pollution in petrochemical industry[J].Contemporary Chemical Industry,2019,48(10):2 385-2 387;2 423.
[4] 余文琴, 李根容,童兰艳,等.食品中苯系物的成因及检测研究进展[J].食品研究与开发,2020,41(9):207-212.
YU W Q,LI G R,TONG L Y,et al.Research progress in the formation and detection of benzene series in food[J].Food Research and Development,2020,41(9):207-212.
[5] 樊孝俊, 邓嘉辉,刘盈智,等.固定污染源中苯系物的便携式气质联用检测方法研究[J].中国环境监测,2015,31(3):139-143.
FAN X J,DENG J H,LIU Y Z,et al.Determination of benzene and its analogies in stationary source emission by portable gas chromatography-mass spectrometry[J].Environmental Monitoring in China,2015,31(3):139-143.
[6] 乌日娜, 金芬,杨莉莉,等.食用油脂中多环芳烃污染来源及控制技术研究展望[J].食品与发酵工业,2015,41(9):225-229;246.
WU R N,JIN F,YANG L L,et al.Review on the pollution sources and control techniques of polycyclic aromatic hydrocarbons (PAHs) in edible oil[J].Food and Fermentation Industries,2015,41(9):225-229;246.
[7] CHOI S,OH S H,KIM Y S,et al.APU SM technology for the production of BTX and LPG from pyrolysis gasoline using metal promoted zeolite catalyst[J].Catalysis Surveys from Asia,2006,10(2):110-116.
[8] 王廷海, 向永生,刘瑞刚,等.裂解碳九加氢工业化试验[J].化工进展,2012,31(5):1 169-1 172.
WANG T H,XIANG Y S,LIU R G,et al.Industrial study on hydrogenation process of cracking C9[J].Chemical Industry and Engineering Progress,2012,31(5):1 169-1 172.
[9] 葛昌平, 李桂明,张力擎,等.碳九关键组分双环戊二烯色谱检测方法应用探析[J].石化技术,2018,25(10):29-30.
GE C P,LI G M,ZHANG L Q,et al.Application of chromatographic determination of dicyclopentadiene,a key component of C9[J].Petrochemical Industry Technology,2018,25(10):29-30.
[10] 徐蓓, 顾正桂,顾黄宗,等.裂解碳九芳烃中苯乙烯的抽提方案设计及优化研究[J].现代化工,2014,34(9):118-121.
XU B,GU Z G,GU H Z,et al.Extractive process design and optimization research on refining styrene from cracked C9 Arene[J].Modern Chemical Industry,2014,34(9):118-121.
[11] 宗旭. 国内乙烯裂解碳九综合利用进展现状[J].化学工程与装备,2013(5):155-157.
ZONG X.Current status of domestic comprehensive utilization of ethylene cracking C9[J].Chemical Engineering & Equipment,2013(5):155-157.
[12] 陈远新. 萃取精馏与连续精馏结合从加氢裂解C9中提取三甲苯馏分的研究[D].南京:南京师范大学,2011.
CHEN Y X.Study on extractive distillation combined with continuous distillation to extract trimethylbenzene fraction from hydrocracking C9[D].Nanjing:Nanjing Normal University,2011.
[13] 杜郢,朱博雅.碳九芳烃的综合利用[J].化学与粘合,2004,26(6):350-353.
DU Y, ZHU B Y. The synthetic utilization of C9 aromatic compound[J]. Chemistry and Adhesion, 2004, 26(6):350-353.
[14] 张敏.裂解碳九综合利用[D].北京:北京化工大学,2013.
ZHANG M. Research on comprehensive utilization of pyrolysis Carbon9[D]. Beijing: Beijing University of Chemical Technology,2013.
[15] 孙桂芳, 刘国文.GC-MS与GC-FTIR联合测定C9烃馏分组成[J].分析测试学报,2004,23(Z1):289-292.
SUN G F,LIU G W.Analysis of the C9 fraction by GC-MS and GC-FTIR[J].Journal of Instrumental Analysis,2004,23(Z1):289-292.
[16] 郑新梅. 乙烯副产品碳九的GC/MS分析[J].河南化工,2001,18(8):33-34.
ZHENG X M.GC/MS analysis of mixed C9 of ethylene by-production[J].Henan Chemical Industry,2001,18(8):33-34.
[17] 中华人民共和国工业和信息化部. SH/T 1793—2015 工业用裂解碳九组成的测定气相色谱法[S].北京:中国石化出版社,2015.
Ministry of Industry and Information Technology of the People’s Republic of China.SH/T 1793—2015 Craking C9 for industrial use-determination of components-gas chromatographic method[S].Beijing:China Petrochemical Press,2015.
[18] 肖洁, 张珰妮,张天闻.海洋污染事故中裂解碳九主要成分的监测分析方法[J].渔业研究,2020,42(2):138-145.
XIAO J,ZHANG D N,ZHANG T W.Analysis methods for detecting cracking C9 in sudden marine pollution accidents[J].Journal of Fisheries Research,2020,42(2):138-145.
[19] 魏富梅. 水中挥发性有机物检测技术的进展[J].中国医药指南,2018,16(32):11-12.
WEI F M.Process in detection technology of volatile organic compounds in water[J].Guide of China Medicine,2018,16(32):11-12.
[20] 房贤文, 石仁德,魏巍.生活饮用水中挥发性有机物的几种检测方法[J].化工管理,2019 (20):45-46.
FANG X W,SHI R D,WEI W.Several detection methods of volatile organic compounds in drinking water[J].Chemical Enterprise Management,2019 (20):45-46.
[21] ANGELI I, CASATI S, RAVELLI A, et al. A novel single-step GC-MS/MS method for cannabinoids and 11-OH-THC metabolite analysis in hair[J]. Journal of Pharmaceutical and Biomedical Analysis, 2018, 155:1-6.
[22] 刘笑笑, 张菁菁,李晨曦,等.低温微萃取-GCMS/MS联用技术检测食用油中20种持久性污染物[J].食品与发酵工业,2020,46(5):292-298.
LIU X X,ZHANG J J,LI C X,et al.Determination of 20 persistent pollutants in edible oil by freeze-degreasing-GCMS/MS[J].Food and Fermentation Industries,2020,46(5):292-298.