文章信息
- 蔡振辉
- CAI Zhenhui
- QuEChERS结合GC-MS/MS检测食用菌中20种有机磷农药的基质效应
- Matrix effect of 20 kinds of organophosphorus pesticides in edible fungi by GC-MS/MS combined with QuEChERS
- 亚热带农业研究, 2022, 18(2): 136-141
- Subtropical Agriculture Research, 2022, 18(2): 136-141.
- DOI: 10.13321/j.cnki.subtrop.agric.res.2022.02.011
-
文章历史
- 收稿日期: 2022-02-22
有机磷农药具有高效、广谱、降解快、价格低等特点,被广泛应用于农业生产[1],但其残留问题也日益凸显,严重影响生态环境,甚至危害人体健康。目前,世界各国通过制定农产品中有机磷农药的最大残留限量加以监督,我国颁布的《食品安全国家标准食品中农药最大残留限量:GB 2763—2021》[2]也规定了农产品中有机磷农药的最大残留限量。当前,对有机磷农药残留的检测主要有气相色谱法[3]、气相色谱—串联质谱法(gas chromatography-tandem mass spectrometry, GC-MS/MS)[4]和液相色谱—串联质谱法[5]等。其中,GC-MS/MS相比气相色谱—质谱联用法(GC-MS)优势明显,可通过监测多个离子,有效排除假阳性干扰,提高分析的选择性和灵敏度[6-7]。
基质效应指样品中其他成分对待测物测定造成干扰,导致测定时无法对目标组分进行准确定量[8]。基质匹配校准法是目前农药残留检测领域降低基质效应的常用方法,但现实中待测样品种类繁多,难以对不同样品一一配制基质标准曲线[9]。因此,寻找代表性的样品基质并配制标准曲线进行定量分析,成为提高检测效率及准确性的有效手段[10-12]。QuEChERS法是一种方便快捷、经济实用、安全可靠的样品前处理技术,常用于农产品农药的检测[13]。目前,有关食用菌中有机磷农药的基质效应报道不多,且涉及的食用菌种类较少[14-16]。因此,本研究采用QuEChERS结合GC-MS/MS法探讨20种有机磷农药在10种常见食用菌中的基质效应,以期为食用菌中有机磷农药的检测提供参考。
1 材料与方法 1.1 材料与试剂 1.1.1 材料试材分别为香菇(Lentinula edodes)、秀珍菇(Pleurotus geesteranus)、蟹味菇(Hypsizygus marmoreus)、杏鲍菇(Pleurotus eryngii)、鹿茸菇(Lyophyllum decastes)、金针菇(Flammulina velutipes)、平菇(Pleurotus ostreatus)、双孢菇(Agaricus bisporus)、黑木耳(Auricularia heimuer)和银耳(Tremella fuciformis)等10种常见食用菌,均为市售。试验前对所有试材进行检测,均未发现待测农药。
1.1.2 试剂20种有机磷农药标准样品由农业部环境保护科研监测所提供,质量浓度均为100 μg·mL-1。丙酮、乙腈、正己烷均为色谱纯,上海安谱公司;5982-5056分散固相萃取试剂盒,美国Agilent公司。
1.2 仪器与设备GCMS-TQ8040三重四极杆气相色谱质谱联用仪,日本岛津制作所;3-18KS高速冷冻离心机,德国Sigma公司;EVAP-111氮吹仪,美国Organomation公司;FP3010料理机,德国Braun公司;移液枪(10~100 0 μL),德国VITILAB公司;PL602E电子天平,METTLER TOLEDO科技有限公司。
1.3 试验方法 1.3.1 样品前处理参照文献[17]中有关农药残留检测的方法进行样品前处理。分别将10种食用菌样品充分搅碎、混匀,准确称取10.0 g样品于50 mL离心管,加入20.0 mL乙腈、5.0 g NaCl,充分震荡混匀后超声提取15 min,5 000 r·min-1离心3 min。吸取8.0 mL上清液置于分散固相萃取净化管,振荡1 min后5 000 r·min-1离心3 min,吸取1.0 mL上清液于10.0 mL试管,40 ℃氮吹至近干,加入1.0 mL丙酮与正己烷混合液复溶,收集净化后的空白基质液备用。
1.3.2 GC-MS/MS条件(1) 色谱条件。色谱柱:Aglient DB-5 MS (30 m×0.25 mm,0.25 μm);梯度升温设置:起始温度50 ℃,保持1 min,以25 ℃·min-1升温至125 ℃,再以10 ℃·min-1升温至290 ℃。(2)质谱条件。EI离子源温度:200 ℃,接口温度:250 ℃,溶剂延迟时间:3.0 min;多反应监测(multiple reaction monitoring, MRM)检测。
1.3.3 标准曲线的绘制分别使用10种食用菌基质液和丙酮正己烷(体积比为1∶1)将标准物质稀释为25、50、100、200、300、400 μg·L-1的标准曲线工作液,上机进样,以外标法绘制基质标准曲线和溶剂标准曲线。
1.3.4 基质效应的计算基质效应(ME)的计算公式为:ME/%=
通过GC-MS/MS对20种农药的混合标准溶液(1 μg·mL-1)进行全扫描,扫描范围设置为50~500 amu。依据相对丰度、质荷比等因素确定1~3个离子,并将其作为母离子进行产物离子扫描,获得最佳子离子和优化电压,确定最优MRM扫描质谱条件。经过质谱优化后获得20种农药组分的母离子、子离子和碰撞能量(表 1)。
| 农药 | 目标离子对 | 碰撞能 | 参考离子对1 | 碰撞能 | 参考离子对2 | 碰撞能 | |||||
| m·z-1 | eV | m·z-1 | eV | m·z-1 | eV | ||||||
| 辛硫磷 | 130.0/103.0 | 10 | 103.0/76.0 | 15 | 130.0/76.0 | 25 | |||||
| 敌敌畏 | 185.0/93.0 | 14 | 185.0/109.0 | 14 | 185.0/63.0 | 22 | |||||
| 灭线磷 | 200.0/158.0 | 6 | 200.0/114.0 | 14 | 200.0/97.0 | 24 | |||||
| 甲拌磷 | 260.0/75.0 | 8 | 260.0/231.0 | 4 | 260.0/47.0 | 4 | |||||
| 特丁硫磷 | 231.0/174.9 | 14 | 231.0/128.9 | 26 | 231.0/202.9 | 8 | |||||
| 二嗪磷 | 304.1/179.1 | 10 | 304.1/162.1 | 8 | 304.1/137.1 | 26 | |||||
| 氯唑磷 | 257.0/162.0 | 8 | 257.0/119.0 | 18 | 257.0/146.0 | 18 | |||||
| 甲基对硫磷 | 263.0/109.0 | 14 | 263.0/136.0 | 8 | 263.0/246.0 | 6 | |||||
| 杀螟硫磷 | 277.0/260.0 | 6 | 277.0/109.1 | 14 | 277.0/228.0 | 14 | |||||
| 甲拌磷亚砜 | 153.0/97.0 | 12 | 199.0/171.1 | 6 | 153.0/125.1 | 9 | |||||
| 马拉硫磷 | 173.1/99.0 | 14 | 173.1/127.0 | 6 | 173.1/145.0 | 6 | |||||
| 甲拌磷砜 | 153.0/97.0 | 12 | 153.0/125.0 | 6 | 199.0/143.0 | 15 | |||||
| 毒死蜱 | 313.9/257.9 | 14 | 313.9/285.9 | 8 | 313.9/193.9 | 28 | |||||
| 对硫磷 | 291.1/109.0 | 14 | 291.1/137.0 | 6 | 291.1/81.0 | 24 | |||||
| 水胺硫磷 | 289.1/136.0 | 14 | 289.1/113.0 | 6 | 289.1/108.0 | 28 | |||||
| 甲基异柳磷 | 199.0/121.0 | 14 | 241.1/121.1 | 22 | - | - | |||||
| 丙溴磷 | 336.9/266.9 | 14 | 336.9/308.9 | 6 | 336.9/294.9 | 10 | |||||
| 三唑磷 | 257.0/162.0 | 8 | 257.0/134.0 | 22 | 257.0/119.0 | 26 | |||||
| 亚胺硫磷 | 160.0/133.0 | 14 | 160.0/77.0 | 24 | 160.0/105.0 | 18 | |||||
| 伏杀硫磷 | 182.0/111.0 | 14 | 182.0/138.0 | 8 | 182.0/102.0 | 14 |
分别将25、50、100、200、300、400 μg·L-1供试农药的标准溶液进行GC-MS/MS检测,可获得各标准溶液在设定条件下的总离子流图(total ion chromatogram, TIC)(图 1)。对分离的色谱峰进行积分,并以质量浓度为横坐标、定量离子对的峰面积为纵坐标,通过外标法绘制标准曲线。结果显示,20种农药在25~400 μg·L-1范围内具有良好的线性关系,相关系数(R2)范围为0.998~0.999。
|
A.辛硫磷;B.敌敌畏;C.灭线磷;D.甲拌磷;E.特丁硫磷;F.二嗪磷;G.氯唑磷;H.甲基对硫磷;I.杀螟硫磷;J.马拉硫磷;K.毒死蜱;L.甲拌磷砜;M.甲拌磷亚砜;N.对硫磷;O.水胺硫磷;P.甲基异柳磷;Q.丙溴磷;R.三唑磷;S.亚胺硫磷;T.伏杀硫磷。 图 1 GC-MS/MS中20种供试农药标准溶液(100 μg·L-1)的总离子流图 Figure 1 Total ion chromatogram of 20 pesticide standard solutions (100 μg·L-1) by GC-MS/MS |
20种农药在10种食用菌中的基质效应见表 2。由表 2可知,供试农药总体上表现为基质增强效应,仅特丁硫磷、氯唑磷、甲拌磷亚砜、甲拌磷砜、甲基异柳磷在个别食用菌中表现为基质抑制效应。其中,甲基对硫磷在杏鲍菇、银耳、平菇中表现出强基质效应,基质效应范围在51.54%~62.39%之间;杀螟硫磷在杏鲍菇、银耳中表现强基质效应,分别为55.64%、50.81%;伏杀硫磷在鹿茸菇中表现为强基质效应(50.92%);亚胺硫磷在10种食用菌中均表现出强基质效应,范围在78.20%~209.12%之间。除亚胺硫磷外,其他农药对10种食用菌基质效应的相对标准偏差(relative standard deviation,RSD)均<10%,表明不同食用菌对供试有机磷农药具有相似的基质效应。
| 农药 | 基质效应/% | RSD | ||||||||||
| 香菇 | 秀珍菇 | 蟹味菇 | 杏鲍菇 | 鹿茸菇 | 金针菇 | 黑木耳 | 银耳 | 平菇 | 双孢菇 | % | ||
| 辛硫磷 | 16.35 | 32.69 | 29.36 | 37.38 | 31.75 | 32.05 | 15.70 | 29.04 | 33.99 | 24.01 | 5.71 | |
| 敌敌畏 | 7.51 | 16.03 | 9.86 | 14.60 | 11.68 | 5.91 | 5.80 | 13.97 | 16.75 | 11.04 | 3.62 | |
| 灭线磷 | 9.87 | 7.59 | 1.27 | 14.65 | 7.03 | 3.66 | 4.49 | 12.40 | 15.95 | 8.10 | 4.40 | |
| 甲拌磷 | 8.30 | 16.19 | 16.13 | 27.94 | 10.57 | 7.01 | 11.72 | 6.95 | 16.41 | 6.89 | 5.87 | |
| 特丁硫磷 | -2.67 | -1.43 | -1.65 | 11.05 | 6.45 | -0.05 | 5.14 | 11.22 | 16.00 | 13.26 | 6.55 | |
| 二嗪磷 | 10.01 | 2.51 | 8.94 | 16.98 | 11.27 | 3.71 | 4.80 | 9.06 | 13.64 | 9.28 | 4.09 | |
| 氯唑磷 | -0.20 | -4.23 | 3.60 | 9.39 | -0.01 | -4.23 | 1.53 | 3.79 | 13.55 | 4.69 | 5.44 | |
| 甲基对硫磷 | 21.03 | 46.08 | 39.59 | 62.39 | 48.94 | 41.83 | 37.39 | 55.50 | 51.54 | 34.76 | 8.15 | |
| 杀螟硫磷 | 23.48 | 46.81 | 36.62 | 55.64 | 48.01 | 39.29 | 38.62 | 50.81 | 47.15 | 36.32 | 6.51 | |
| 甲拌磷亚砜 | -4.42 | 11.11 | 6.88 | 14.56 | 10.30 | 1.79 | 7.60 | 15.30 | 18.16 | 6.29 | 6.17 | |
| 马拉硫磷 | 0.11 | 18.11 | 10.38 | 22.45 | 15.41 | 10.07 | 10.56 | 19.11 | 23.30 | 17.82 | 6.13 | |
| 甲拌磷砜 | -3.78 | 6.24 | 4.95 | 12.44 | 7.29 | 0.95 | 8.46 | 14.70 | 15.24 | 4.07 | 5.61 | |
| 毒死蜱 | 5.72 | 14.55 | 10.30 | 20.28 | 15.42 | 4.10 | 9.74 | 15.94 | 23.48 | 4.94 | 5.87 | |
| 对硫磷 | 13.20 | 26.96 | 17.91 | 34.36 | 31.06 | 13.61 | 20.16 | 26.40 | 33.18 | 14.77 | 6.69 | |
| 水胺硫磷 | 9.68 | 20.14 | 6.87 | 26.15 | 20.46 | 8.84 | 13.72 | 21.90 | 26.47 | 10.09 | 6.39 | |
| 甲基异柳磷 | 2.46 | 4.15 | 2.88 | 10.21 | 2.65 | -1.76 | 0.67 | 3.54 | 13.30 | 6.55 | 4.27 | |
| 丙溴磷 | 18.96 | 28.68 | 13.74 | 34.80 | 20.15 | 12.98 | 13.03 | 25.82 | 37.70 | 28.31 | 7.31 | |
| 三唑磷 | 30.87 | 39.08 | 24.46 | 47.09 | 36.24 | 15.37 | 25.72 | 35.76 | 49.13 | 41.03 | 7.83 | |
| 亚胺硫磷 | 78.20 | 209.12 | 177.20 | 154.64 | 194.58 | 127.12 | 107.38 | 135.92 | 111.61 | 135.79 | 16.87 | |
| 伏杀硫磷 | 29.37 | 45.67 | 43.89 | 46.27 | 50.92 | 36.70 | 33.85 | 36.14 | 44.75 | 32.57 | 5.11 | |
由表 2还可知,除亚胺硫磷、甲基对硫磷外其他农药的基质效应RSD均<8%。通过考察多种食用菌中亚胺硫磷的基质效应均值与单种食用菌基质效应的关系,同时以RSD值较大的甲基对硫磷为参考筛选代表基质。亚胺硫磷、甲基对硫磷在10种食用菌基质中的平均基质效应值分别为143.16%、43.91%,比对表 2中各种食用菌的基质效应,选择双孢菇作为代表基质。
2.4 以双孢菇为代表基质时有机磷农药的基质效应比较如表 3所示,以双孢菇作为代表基质配制标准曲线,有效校正辛硫磷、对硫磷、甲基对硫磷、杀螟硫磷、三唑磷、亚胺硫磷、伏杀硫磷等农药的强基质效应,除甲基对硫磷、亚胺硫磷在个别食用菌中表现为中等基质效应外,其他农药均表现为弱基质效应。以溶剂和双孢菇为参比,比较其他供试基质的强、中、弱基质效应所占比例(图 2)。其中,强基质效应由8.0%下降为0,中等基质效应由36.5%降低为2.0%,弱基质效应由55.5%提升为98.0%。由此可知,利用GC-MS/MS检测食用菌中有机磷农药时,将双孢菇作为代表基质配置标准曲线能够有效提高定量分析的准确性。
| 农药 | 基质效应/% | |||||||||
| 香菇 | 秀珍菇 | 蟹味菇 | 杏鲍菇 | 鹿茸菇 | 金针菇 | 黑木耳 | 银耳 | 平菇 | 双孢菇 | |
| 辛硫磷 | -6.17 | 7.01 | 4.32 | 10.78 | 6.24 | 6.49 | -6.70 | 4.06 | 8.06 | 0.00 |
| 敌敌畏 | -3.18 | 4.49 | -1.06 | 3.20 | 0.58 | -4.62 | -4.72 | 2.63 | 5.14 | 0.00 |
| 灭线磷 | 1.64 | -0.47 | -6.32 | 6.05 | -0.99 | -4.11 | -3.34 | 3.98 | 7.26 | 0.00 |
| 甲拌磷 | 1.32 | 8.70 | 8.65 | 19.69 | 3.44 | 0.11 | 4.52 | 0.05 | 8.91 | 0.00 |
| 特丁硫磷 | -14.06 | -12.97 | -13.17 | -1.95 | -6.01 | -11.75 | -7.17 | -1.80 | 2.42 | 0.00 |
| 二嗪磷 | 0.67 | -6.20 | -0.31 | 7.05 | 1.82 | -5.09 | -4.10 | -0.20 | 3.99 | 0.00 |
| 氯唑磷 | -4.67 | -8.52 | -1.04 | 4.49 | -4.49 | -8.52 | -3.02 | -0.86 | 8.47 | 0.00 |
| 甲基对硫磷 | -10.19 | 8.40 | 3.58 | 20.50 | 10.52 | 5.24 | 1.95 | 15.39 | 12.45 | 0.00 |
| 杀螟硫磷 | -9.42 | 7.70 | 0.22 | 14.17 | 8.58 | 2.18 | 1.69 | 10.63 | 7.94 | 0.00 |
| 甲拌磷亚砜 | -10.07 | 4.53 | 0.55 | 7.78 | 3.77 | -4.24 | 1.23 | 8.48 | 11.16 | 0.00 |
| 马拉硫磷 | -15.03 | 0.25 | -6.31 | 3.93 | -2.04 | -6.58 | -6.16 | 1.10 | 4.65 | 0.00 |
| 甲拌磷砜 | -7.54 | 2.08 | 0.84 | 8.04 | 3.09 | -3.01 | 4.21 | 10.21 | 10.73 | 0.00 |
| 毒死蜱 | 0.75 | 9.16 | 5.11 | 14.62 | 9.99 | -0.80 | 4.58 | 10.49 | 17.67 | 0.00 |
| 对硫磷 | -1.36 | 10.63 | 2.74 | 17.07 | 14.20 | -1.01 | 4.70 | 10.14 | 16.05 | 0.00 |
| 水胺硫磷 | -0.37 | 9.13 | -2.93 | 14.59 | 9.42 | -1.14 | 3.30 | 10.73 | 14.88 | 0.00 |
| 甲基异柳磷 | -3.84 | -2.25 | -3.45 | 3.43 | -3.66 | -7.80 | -5.52 | -2.83 | 6.34 | 0.00 |
| 丙溴磷 | -7.29 | 0.29 | -11.36 | 5.05 | -6.36 | -11.95 | -11.91 | -1.94 | 7.32 | 0.00 |
| 三唑磷 | -7.20 | -1.39 | -11.75 | 4.30 | -3.40 | -18.20 | -10.86 | -3.74 | 5.74 | 0.00 |
| 亚胺硫磷 | -24.42 | 31.10 | 17.56 | 8.00 | 24.93 | -3.68 | -12.05 | 0.06 | -10.25 | 0.00 |
| 伏杀硫磷 | -2.41 | 9.88 | 8.54 | 10.33 | 13.84 | 3.12 | 0.96 | 2.69 | 9.18 | 0.00 |
|
A.以溶剂为参比;B.以双孢菇为参比。 图 2 9种食用菌不同基质效应占比 Figure 2 Proportions of different matrix effects of 9 edible fungi |
基质效应是农药定量检测中常见的问题,受多种因素影响。有机磷农药具有P=O或P=S等极性基团,这些基团极易被进样口活性位点吸附,从而加强基质效应[18]。此外,农药化合物的空间结构以及基质中糖类、油脂、蛋白质等含量的不同也会影响基质效应[15]。本研究采用QuEChERS结合GC-MS/MS法探讨20种有机磷农药在10种常见食用菌中的基质效应,结果显示,不同食用菌基质对有机磷农药的基质效应相似,且使用双孢菇作为代表基质配制标准曲线进行定量分析,既保证准确率又可提高效率。
| [1] | 刘欣红, 王艳辉, 吴迪, 等. 蔬菜中13种有机磷农药基质效应研究[J]. 农药科学与管理, 2019, 40(7): 28–35. DOI: 10.3969/j.issn.1002-5480.2019.07.008 |
| [2] | 国家卫生健康委员会, 中华人民共和国农业农村部, 国家市场监督管理总局. 食品安全国家标准食品中农药最大残留限量: GB 2763—2021[S]. 北京: 中国标准出版社, 2021. |
| [3] | WANG L G, LIANG Y C, JIANG X. Analysis of eight organophosphorus pesticide residues in fresh vegetables retailed in agricultural product markets of Nanjing, China[J]. Bulletin of Environmental Contamination and Toxicology, 2008, 81(4): 377–382. DOI: 10.1007/s00128-008-9498-2 |
| [4] | HADIAN Z, SAMIRA S, YAZDANPANAH H. Pesticide residues analysis in Iranian fruits and vegetables by gas chromatography-mass spectrometry[J]. Iranian Journal of Pharmaceutical Research, 2019, 18(1): 275–285. |
| [5] | HAN C, HU B Z, CHEN S B, et al. Determination of Xinjunan pesticide residue in foodstuffs of plant origin by a modified QuEChERS method and ultra performance liquid chromatography-tandem mass spectrometry[J]. LWT, 2021, 151: 112101. DOI: 10.1016/j.lwt.2021.112101 |
| [6] | 钟丹丹, 柳春辉, 窦筱艳. 固体废物中挥发性及半挥发性有机污染物前处理和GC/GC-MS测定方法的研究进展[J]. 环境与发展, 2018, 30(8): 64–68. |
| [7] | 蔡振辉, 张春柳, 许亦峰, 等. QuEChERS结合GC-MS/MS检测铁皮石斛中5种农药残留[J]. 亚热带农业研究, 2021, 17(1): 66–71. |
| [8] | 陈鹏. 气质联用和液质联用中基质效应研究进展[J]. 天津药学, 2020, 32(4): 67–72. DOI: 10.3969/j.issn.1006-5687.2020.04.023 |
| [9] | 黄松, 宋安华, 陈穗, 等. QuEChERS-气相色谱-串联质谱法测定油麦菜中农药多残留的基质效应[J]. 食品安全质量检测学报, 2019, 10(12): 3863–3868. DOI: 10.3969/j.issn.2095-0381.2019.12.038 |
| [10] | 许炳雯, 葛祥武, 刘雪红, 等. 蔬菜水果中有机磷农药的基质效应分析及检测方法选择[J]. 食品安全质量检测学报, 2020, 11(23): 8784–8792. |
| [11] | 吴姻如, 汤旭东, 孙超斌, 等. 基质效应对瓜类蔬菜中有机磷农药残留测定的影响研究[J]. 农产品质量与安全, 2016(4): 62–64. DOI: 10.3969/j.issn.1674-8255.2016.04.019 |
| [12] | ZHAO J, PU J J, WU X P, et al. Evaluation of the matrix effect of pH value and sugar content on the analysis of pesticides in tropical fruits by UPLC-MS/MS[J]. Microchemical Journal, 2021, 168: 106375. DOI: 10.1016/j.microc.2021.106375 |
| [13] | 戴唯, 李巧, 朱明, 等. QuEChERS-同位素内标-高效液相色谱-串联质谱法测定动物源性食品中植物生长调节剂类农药残留[J]. 色谱, 2021, 39(11): 1213–1221. |
| [14] | 夏珍珍, 彭西甜, 夏虹, 等. 基于GC-MS/MS的蔬菜多农药残留基质效应研究[J]. 湖北农业科学, 2016, 55(24): 6578–6581. |
| [15] | 戚文华, 李莺, 王毅红. 气相色谱-串联质谱法测定鲜香菇中25种有机磷类农药残留[J]. 中国瓜菜, 2019, 32(2): 18–22. DOI: 10.3969/j.issn.1673-2871.2019.02.005 |
| [16] | 易盛国, 侯雪, 韩梅, 等. 气相色谱-串联质谱法检测蔬菜农药残留基质效应与基质分类的研究[J]. 西南农业学报, 2012, 25(2): 537–543. DOI: 10.3969/j.issn.1001-4829.2012.02.038 |
| [17] | 国家卫生健康委员会, 农业农村部, 国家市场监督管理总局. 食品安全国家标准植物源性食品中208种农药及其代谢物残留量的测定气相色谱-质谱联用法: GB 23200.113—2018[S]. 北京: 中国标准出版社, 2018. |
| [18] | 方海仙, 耿慧春, 梅文泉, 等. 气相色谱法测定农产品中农药残留的基质效应研究[J]. 食品安全质量检测学报, 2018, 9(14): 3770–3779. DOI: 10.3969/j.issn.2095-0381.2018.14.037 |
2022,Vol. 18

