兽药在防治动物疾病、提高生产效率、改善畜产品质量等方面起着十分重要的作用,但部分养殖人员在养殖过程滥用兽药,甚至使用违禁药物、假劣药物,造成了严重的兽药残留问题[1-5]。动物性食物兽药残留在人体内不断的蓄积会引起人体的多种急慢性反应,如莱克多巴胺(ractopamine,RAC)和氯霉素(chloramphenicol,CAP)对人体有毒副作用,盐酸克仑特罗(clenbuterol,CL)有致癌风险,氟苯尼考(florfenicol,FF)具有一定的免疫毒性和胚胎毒性[6-10]。动物源性食品中兽药的残留不仅直接危害着人类生命健康,同时也极大危害着畜牧业的发展和生态环境,已引起国际有关组织及许多国家、地区的高度重视。
现阶段,根据我国第176号农业部公告,RAC禁止在饲料和动物饮用水中使用;根据我国第235号农业部公告,CL及CAP在动物性食品中不得检出,FF在猪肉组织中残留量不得超过300 μg/kg。研究表明,RAC、CL、CAP和FF在猪、牛等动物体内易残留、消除缓慢,主要经尿液排出体外[11-14]。通过检测尿液中排出的药物浓度,可对屠宰前动物体内药物的残留情况进行监测,对有效监控食品安全有重要意义。
目前,国内外报道的关于RAC、CL、CAP、FF及其代谢物氟苯尼考胺(florfenicol amine,FFA)的传统检测方法主要包括气相色谱法、高效液相色谱(high performance liquid chromatography,HPLC)法、超高效液相色谱法、液相色谱-串联质谱法、液相色谱-四极杆-飞行时间串联质谱法等。这些方法准确度高、稳定性好,被广泛应用于兽药残留检测[15-22]。但上述方法多数需要配合大仪器,且操作繁琐、检测耗时较长、前处理复杂,不适用于快速检测分析。常用快速检测方法主要有酶联免疫分析(enzyme-linked immunosorbent assay,ELISA)法、化学发光免疫分析法、胶体金免疫层析(colloidal gold immunochromatography assay,GICA)法等[23-27],这些方法或需专业人员操作、或只能检测单一指标,无法实现兽药多残留的简便快速检测[28-29]。免疫蛋白芯片测定方法相对于以上各种方法具有低成本、高通量、灵敏度高等优点[30-33]。本实验利用免疫蛋白芯片技术建立了同时检测猪尿中RAC、CL、CAP、FF及其代谢物FFA的方法,同时基于智能手机开发了与免疫芯片匹配的检测应用程序(application,APP),适用于现场大量筛选分析。
驴抗山羊血清、HRP标记羊抗鼠血清、HRP羊抗兔血清 广州优抗多生物技术有限公司;RAC单克隆抗体、RAC抗原、CL单克隆抗体、CL抗原、CAP单克隆抗体、CAP抗原、FF单克隆抗体、FF抗原 广东省食品质量安全重点实验室;RAC(GR,100 mg/瓶)、CL(GR,100 mg/瓶)、CAP(GR,250 mg/瓶)、FF(GR,250 mg/瓶)、FFA(GR,100 mg) 德国Dr. Ehrenstorfer公司;紫色沉淀型单组分3,3’,5,5’-四甲基联苯胺(3,3’,5,5’-tetramethylbenzidine,TMB)显色液湖州英创生物科技有限公司;P液:含10%蔗糖、0.4%鱼明胶、0.5%酪蛋白的磷酸盐缓冲液;T液:含8%蔗糖、0.4%鱼明胶、0.5%酪蛋白的三羟甲基氨基甲烷盐酸缓冲液;PBST:含0.1% Tween-20的磷酸缓冲液。
免疫芯片反应盒 广东省食品质量安全重点实验室;硝酸纤维素膜(NC膜) 美国Millipore公司;M227fdw多功能一体机 中国惠普有限公司。
1.3.1 免疫芯片分析方法基本步骤
本免疫芯片测定兽药残留的基础是目标检测物与抗原竞争抗体的反应。免疫芯片上固定有药物抗原,检测时在反应格加入药物的单克隆抗体和药物标准溶液或待测样品,待抗原抗体反应后洗涤,再加入HRP标记二抗,反应后洗涤并加入显色液,最后通过显色强度进行药物定量分析。基本步骤如下:
点样包被:用碳酸盐缓冲液(carbonate buffer solution,CBS)稀释驴抗山羊血清作为质控点样液备用,稀释RAC抗原、CL抗原、CAP抗原、FF抗原作为检测点样液备用。将孔径为0.22 μm的NC膜裁剪成合适大小,粘贴于PVC板上后置于反应盒中,吸取点样液按照0.5 cm×0.5 cm的间距点样,每个样品点样量为1 μL。将点样芯片置于37 ℃恒温水浴锅中孵育1 h使蛋白质固定在芯片上。
封闭:包被抗原或抗体后的NC膜上尚有未被抗原或抗体填充的空隙,需对NC膜进行封闭处理,防止非特异性吸附。将上述芯片用PBST洗涤2 次,双蒸水洗涤1 次,反应格中加入500 μL含1%牛血清蛋白的封闭液,于37 ℃封闭30 min。
一抗孵育:将上述芯片用PBST洗涤2 次,双蒸水洗涤1 次。反应格中加入T液稀释的单克隆抗体及标准品溶液或待测样品液,置于摇床上反应15 min。
二抗孵育:将上述芯片用PBST、双蒸水各洗涤2 次。反应格中加入200 μL用T液稀释的HRP标记羊抗鼠、HRP标记羊抗兔混合液,置于摇床上反应15 min。
显色反应:将上述芯片用PBST、双蒸水各洗涤2 次。反应格中加入200 μL TMB显色液,反应3 min,随后将芯片取出置于洁净滤纸,风干后用扫描仪扫描,运用Photoshop软件获取显色斑点灰度值进行数据处理分析。
1.3.2 标准曲线与智能手机检测系统的建立
以确定的最佳工作质量浓度、最佳反应条件制备免疫芯片,以6 个不同质量浓度的药物标准溶液进行检测,每个质量浓度重复3 次,使用广东省食品质量安全重点实验室开发的智能手机APP“天眼测”获取芯片显色斑点灰度值与空白值(CBS点样)的差值ΔI,以药物质量浓度(ng/mL)为横坐标,灰度值差值ΔI作为纵坐标绘制标准曲线。在检测软件中添加兽药残留检测项目,导入RAC、CL、CAP、FF兽药残留检测标准曲线。检测时配合实验室开发的图片采集暗箱装置,运行检测APP识别芯片显色斑点,软件即可显示兽药残留浓度并做出阴阳性判断。
1.3.3 基于智能手机检测兽药多残留免疫芯片分析方法的验证评价
选取经HPLC确证未被RAC、CL、CAP及FF污染的20 份猪尿样品,离心取上清液后向其中分别添加3 种质量浓度的RAC、CL、CAP、FF/FFA药物标准品,采用基于智能手机免疫芯片的分析方法进行检测,计算检测限(limit of detection,LOD),分析其添加回收率与批内、批间变异系数(coefficient of variation,CV)。
1.3.4 与其他快速检测方法比较实验
采集30 份猪尿样品,经前处理后向猪尿中添加不同质量浓度的RAC、CL、CAP、FF/FFA标准品,采用基于智能手机的免疫芯片的分析方法进行检测,同时与市售胶体金免疫层析检测卡、ELISA试剂盒检测进行比较。
NC膜固定蛋白质的能力与其孔径大小相关,对于低分子质量的蛋白,NC膜孔径越小,膜结合越牢固;点样稀释液、抗体稀释液的选择影响抗体、抗原在液体中的分散均匀性、显色效果,从而影响芯片检测灵敏度;抗原包被时间影响蛋白质固定效果及有效性;显色反应时间影响显色效果,反应时间太短则显色不明显,反应时间太长易产生背景干扰。
选用孔径大小分别为0.22、0.45 μm的NC膜作为固定蛋白质的载体制备检测单种兽药残留的免疫蛋白芯片,选用CBS和PBS作为点样稀释液稀释驴抗羊血清与抗原,选用P液、T液、PBS作为抗体稀释液稀释单克隆抗体,选择0.5、1、2、8 h作为包被时间固定蛋白质,选择1、3、5、10、20 min作为显色时间进行显色反应,观察芯片显色效果与检测灵敏度。结果显示,当使用0.22 μm的NC膜作为蛋白固定载体、CBS作为点样稀释液、T液作为抗体稀释液时,显色反应后芯片背景干扰较少、点样斑点形状较规则、显色较均匀、灵敏度较高;当包被时间为1 h(37 ℃)时芯片对蛋白固定量最大,斑点显色最深;显色反应时间3 min时斑点显色明显,背景干扰较少。
抗原、抗体的工作质量浓度影响芯片显色效果与检测灵敏度,随着点样抗原与抗体质量浓度的增高,芯片斑点显色越深,灰度值越高,药物与抗原竞争抗体的反应越不明显,灵敏度越低。使用棋盘法优化本实验中抗原、抗体的最佳工作质量浓度,各抗原、抗体稀释倍数如表1所示。
表1 各抗原、抗体稀释倍数的确定
Table 1 Determination of optimal dilution factors of coating antigens and monoclonal antibodies
images/BZ_290_1276_1491_2256_1538.pngRAC 1∶250 1∶500 1∶1 000 1∶2 000 1∶4 000 1∶15 000 1∶30 000 1∶6 000 1∶12 000 CL 1∶250 1∶500 1∶1 000 1∶2 000 1∶4 000 1∶10 000 1∶20 000 1∶40 000 1∶80 000 CAP 1∶625 1∶1 250 1∶2 500 1∶5 000 1∶10 000 1∶150 0001∶300 000 1∶60 000 1∶120 000 FF 1∶25 1∶50 1∶100 1∶200 1∶400 1∶50 000 1∶100 0001∶200 000 1∶400 000
在综合检测斑点显色强度、竞争反应灵敏度、原料用量的基础上,根据实验结果最终确定RAC抗原、CL抗原、CAP抗原、FF抗原最佳稀释倍数为1∶1 000、1∶500、1∶2 500、1∶100,RAC抗体、CL抗体、CAP抗体、FF抗体最佳稀释倍数分别为1∶30 000、1∶400 000、1∶300 000、1∶200 000。在最优反应条件下,免疫芯片对RAC、CL、CAP、FF的检测灵敏度分别为2、2、0.5、50 ng/mL。免疫芯片检测不同质量浓度RAC的结果如图1所示(CL、CAP、FF检测结果图类似不再列出)。
图1 免疫芯片分析RAC残留的检测灵敏度
Fig. 1 Sensitivity of immunochip assay for RAC
当4 种检测项目集成于同一免疫芯片时,4 种检测对象的抗原-抗体-药物之间特异性影响检测结果的准确性。根据间接竞争的原理,对抗原-抗体-药物之间的特异性进行评价。以确定的最优条件进行实验,将免疫芯片裁剪成2.5 cm×0.5 cm,吸取质控、检测点样液按照5×1(驴抗山羊血清、RAC抗原、CL抗原、CAP抗原、FF抗原)阵式点样,其余步骤如1.3.1节,其中一抗孵育时1~4号点样芯片所在反应格中分别加入RAC单克隆抗体、CL单克隆抗体、CAP单克隆抗体、FF单克隆抗体和缓冲溶液,5~8号点样芯片所在反应格中加入4 种药物单克隆抗体的混合液及RAC、CL、CAP、FF单种药物标准溶液。
从图2a可知,使用包被了驴抗羊血清及多种药物抗原的免疫芯片进行实验,当一抗孵育只加入一种单克隆抗体及缓冲液时,只有质控区及该单克隆抗体对应的抗原包被区显现斑点,其余抗原包被区不显现斑点,说明该实验中一种单克隆抗体只与驴抗羊血清及其所对应的抗原产生特异性结合,不与其余3 种药物抗原发生明显特异性结合。从图2b可看出,当一抗孵育只加入混合抗体及单种药物标准品时,只有该药物对应检测区不显现斑点,说明一种药物只特异性抑制对应药物抗原与其单克隆抗体的结合,对其余3 种药物抗原与其对应抗体的结合无明显抑制。由此可知,各抗体对其抗原的特异性较好,同时检测4 种兽药残留时各对象的检测无明显干扰。
图2 抗体特异性分析(a)和各分析物交叉反应(b)结果
Fig. 2 Evaluation of antibody specificity (a) and cross-reactivities between analytes (b)
以药物质量浓度(ng/mL)为横坐标、灰度值差值ΔI作为纵坐标绘制标准曲线,曲线方程及免疫芯片检测4 种兽药残留的分析参数如表2所示,该方法对RAC、CL、CAP、FF的检出限(20 份空白样品平均值与3 倍标准差的和)分别为0.146、0.137、0.035、3.7 3 n g/m L,线性范围(I C 2 0~I C 8 0)分别为0.35 3 ~1.0 17、0.309~0.897、0.082~0.249、9.424~35.594 ng/mL,相关系数大于0.99。
表2 免疫芯片检测4 种兽药残留的分析参数
Table 2 Figures of merit of the immunochip assays for the detection of four veterinary drugs
药物 标准曲线 R2 线性范围/(ng/mL)IC50/(ng/mL)LOD/(ng/mL)芯片灵敏度/(ng/mL)images/BZ_291_1299_601_2279_685.png1+( )0.560 CL y=images/BZ_291_1432_798_1601_828.png3.296 -0.095 1+( )0.527 0.998 0.309~0.897 0.527 0.137 2 CAP y=images/BZ_291_1426_940_1590_971.png2.22 -0.095 7 1+( )0.527 0.991 0.082~0.249 0.143 0.035 0.5images/BZ_291_1299_1076_2279_1168.png2.544
使用本方法测出FFA的IC50,计算交叉反应率(交叉反应率/%=IC50(FF)/IC50(FFA)×100)。结果显示IC50(FFA)为22.465 ng/mL,交叉反应率为81.5%。说明该方法可同时检测FF与FFA。
使用20 份猪尿样品进行药物添加检测,回收率与CV如表3所示,添加回收率在85.69%~110.66%之间,批内CV不高于12.53%,批间CV小于15%。由此可见,本方法具有良好的回收率,可同时检测猪尿样品中RAC、CL、CAP、FF/FFA的残留量。
表3 猪尿中各药物回收率和CV值
Table 3Recoveries and inter-batch and intra-batch coefficients of variation of four veterinary drugs spiked into pig urine samples
images/BZ_291_1299_2073_2277_2167.pngRAC 20 0.25 105.45 6.84 99.23 8.65 0.5 110.66 4.61 106.12 10.05 1 90.06 6.13 88.51 7.64 CL 20 0.25 104.5 4.75 103.23 10.90 0.5 106.45 3.41 104.53 12.32 1 94.5 8.33 105.45 14.59 CAP 20 0.062 5 102.9 12.53 106.26 10.70 0.125 88.6 4.59 95.42 12.89 0.25 97.65 11.63 91.166 6.64 FF 20 6.25 98.4 5.41 94.67 7.70 12.5 103.4 3.02 97.54 7.42 25 109.7 8.39 102.78 5.10 FFA 20 6.25 85.69 7.75 86.04 8.09 12.5 87.69 8.08 91.22 9.69 25 86.38 8.43 91.35 12.31
使用基于智能手机的免疫芯片法、GICA、ELISA法检测猪尿中RAC、CL、CAP、FF/FFA残留量,检测结果如表4所示,免疫芯片法与市售胶体金免疫层析法检测结果的符合率为99.2%、与ELISA法检测结果的符合率为98.3%,说明本方法与常用快速检测方法相比有较高的符合性,可用于现场快速检测猪尿中RAC、CL、CAP、FF/FFA的残留量。
表4 免疫芯片法与GICA法、ELISA试剂盒法检测结果对比
Table 4 Comparison of results of detection using immunochip assay,GICA and ELISA
images/BZ_292_201_856_1180_950.pngRAC 20 10 20 10 20 10 CL 17 13 18 12 18 12 CAP 15 15 15 15 14 16 FF/FFA 22 8 22 8 21 9
本实验研制了基于智能手机检测多种兽药的免疫芯片,能同时检测猪尿中RAC、CL、CAP、FF及其代谢物FFA。该方法对RAC、CL、CAP、FF的检出限分别为0.146、0.137、0.035、3.73 ng/mL,线性范围分别为0.353~1.017、0.309~0.897、0.082~0.249、9.424~35.594 ng/mL。猪尿中药物添加回收率为85.69%~110.66%,批内CV不高于12.53%,批间CV小于15%,与GICA和ELISA两种传统快速检测方法相比检测符合率大于98%,是一种方便快速、高通量、低成本检测多种兽药残留的方法,适用于现场大量筛选分析。
[1] 矫秀燕. 我国动物性食品中兽药残留的现状、问题及对策[J]. 中国畜牧兽医文摘, 2017, 33(1): 28.
[2] 黄升达, 袁成进, 严小东. 浅谈动物产品药物残留的原因[J]. 安徽农学通报, 2018, 24(9): 139-140.
[3] 李研东, 霍惠玲, 李肖莉. 动物源性食品中兽药残留检测现状及思考[J]. 中国畜收兽医文摘, 2017, 33(4): 4-5.
[4] MERSMANN H J. Overview of the effects of beta-adrenergic receptor agonists on animal growth including mechanisms of action[J]. Journal of Animal Science, 1998, 1(76): 160. DOI:10.2527/1998.761160x.
[5] GUO J J, TSAI K, KELTON CHRISTINA M L. Risk of serious asthma exacerbations associated with long-acting beta agonists among patients with asthma: a retrospective cohort study[J]. Annals of Allergy, Asthma & Immunology, 2011, 3(106): 214-222. DOI:10.1016/j.anai.2010.12.008.
[6] 罗奕铭. 盐酸克伦特罗与莱克多巴胺联检速测卡的研制[J]. 现代食品科技, 2013, 29(1): 197-200.
[7] 吴军明. 浅析动物性食品中兽药残留对人体健康的危害及解决措施[J]. 畜禽业, 2018, 29(8): 65-66.
[8] 何资颖, 宋小云, 王琤韡. 莱克多巴胺的危害及防范措施[J]. 江西饲料, 2017(5): 15-17. DOI:10.3969/j.issn.1008-6137.2017.05.004.
[9] 陈秀杰, 刘佳林. 氟苯尼考的理化性质及毒性[J]. 养殖技术顾问,2011(1): 173. DOI:10.3969/j.issn.1673-1921.2011.01.173.
[10] 李洪权, 柏玉冰, 宋忠祥, 等. 动物源性食品中氟苯尼考残留水平和分布情况分析[J]. 食品安全质量检测学报, 2019, 10(2): 426-433.DOI:10.3969/j.issn.2095-0381.2019.02.024.
[11] 田甜. 猪尿盐酸克伦特罗快速检测卡优选及应用[D]. 雅安: 四川农业大学, 2014.
[12] 李然, 林泽佳, 杨金易, 等. 酶联免疫法检测动物组织及尿液中氟苯尼考与甲砜霉素的残留[J]. 分析化学, 2018, 46(8): 1321-1328.DOI:10.11895/j.issn.0253-3820.181084.
[13] 汤超华, 张军民, 李励军. 莱克多巴胺在反刍动物体内代谢及残留规律研究进展[C]. 北京中国毒理学会兽医毒理学与饲料毒理学学术讨论会暨兽医毒理专业委员会第4次全国代表大会, 2012.
[14] 郭平. 供港生猪禁限用兽药残留高通量检测方法的建立和应用[D].北京: 中国农业大学, 2017.
[15] GONZÁLEZ A, JESSICA A, VÍCTOR C. Estrogens determination in wastewater samples by automatic in-syringe dispersive liquidliquid microextraction prior silylation and gas chromatography[J].Journal of Chromatography A, 2015, 413: 1-8. DOI:10.1016/j.chroma.2015.08.031.
[16] YANG Y Y, HUI L, ZHANG J X. Analysis of infant formula for steroid hormones by gas chromatography-tandem mass spectrometry using microwave-assisted extraction and gel permeation chromatography clean up[J]. Food Analytical Methods, 2014, 4(7): 798-805.DOI:10.1007/s12161-013-9683-8.
[17] 伍华雯. 高效液相色谱-串联质谱法同时测定供港生猪尿液中29 种限用兽药残留[J]. 分析测试学报, 2016, 35(10): 1261-1266.DOI:10.3969/j.issn.1004-4957.2016.10.007.
[18] 邵会, 冷凯良, 周明莹, 等. 水产品中氯霉素、甲砜霉素、氟苯尼考、氟苯尼考胺多残留的同时测定: GC/MS法[J]. 渔业科学进展,2015, 36(3): 137-141. DOI:10.11758/yykxjz.20150321.
[19] 姚佳, 王昕, 张建新, 等. 免疫亲和柱-高效液相色谱法测定牛奶中氯霉素和玉米赤霉醇及其类似物[J]. 食品科学, 2014, 35(18): 124-127.DOI:10.7506/spkx1002-6630-201418025.
[20] 谭策, 赵琪, 甄珍, 等. 动物源性食品中β-受体兴奋剂类多残留检测方法研究[J]. 现代食品, 2018(19): 141-145. DOI:10.16736/j.cnki.cn41-1434/ts.2018.19.042.
[21] 郝杰, 姜洁, 余建龙, 等. 固相萃取-超高效液相色谱-串联质谱法同时测定动物源性食品中多种兽药残留[J]. 食品科学, 2017, 38(12):266-272. DOI:10.7506/spkx1002-6630-201712041.
[22] 张静, 高玉时, 周倩, 等. 高效液相色谱-电喷雾电离串联质谱法测定鸡肉、鸡蛋中氯霉素、甲砜霉素和氟苯尼考药物残留量[J].食品安全质量检测学报, 2017, 8(6): 2211-2219. DOI:10.3969/j.issn.2095-0381.2017.06.043.
[23] REN M L, XU H Y, HUANG X L. Immunochromatographic assay for ultrasensitive detection of aflatoxin B1 in maize by highly luminescent quantum dot beads[J]. ACS Applied Materials & Interfaces, 2014,16(6): 14215-14222. DOI:10.1021/am503517s.
[24] WANG L Y, YAO M W, FANG X Y, et al. Gold magnetic nanoparticles-based chemiluminescent immunoassay for detection of chloramphenicol in milk[J]. MATEC Web of Conferences, 2017, 88:2007. DOI:10.1051/matecconf/20178802007.
[25] JANNE L, LAMMINMÄKI U, LÖVGREN T. Multiresidue detection of fluoroquinolones: specificity engineering of a recombinant antibody with oligonucleotide-directed mutagenesis[J]. Journal of Agricultural and Food Chemistry, 2013, 49(61): 11981-11985. DOI:10.1021/jf403715n.
[26] 李超辉. 胶体金免疫层析试纸条定量检测猪尿以及猪肉中的克伦特罗残留[D]. 南昌: 南昌大学, 2014.
[27] WEN K, NÖLKE G, SCHILLBERG S. Improved fluoroquinolone detection in ELISA through engineering of a broad-specific single-chain variable fragment binding simultaneously to 20 fluoroquinolones[J]. Analytical and Bioanalytical Chemistry, 2012,403: 2771-2783. DOI:10.1007/s00216-012-6062-z.
[28] 邓家珞, 陆利霞, 熊晓辉, 等. 兽药残留快速检测方法比较分析[J]. 生物加工过程, 2018, 16(2): 42-48. DOI:10.3969/j.issn.1672-3678.2018.02.006.
[29] 朱晓艾, 陈雅君, 陈雪香, 等. 高效液相色谱法检测猪肉中盐酸克伦特罗的前处理方法[J]. 食品与机械, 2015, 31(2): 106-109.
[30] 左鹏. 蛋白芯片法快速测定食品中氯霉素和磺胺二甲嘧啶残留[J]. 食品科学, 2007, 28(2): 254-258. DOI:10.3321/j.issn:1002-6630.2007.02.065
[31] 孙来玉, 尹莉, 姚婷, 等. 瘦肉精克伦特罗残留检测免疫芯片技术的研究[J]. 药物分析杂志, 2014, 34(5): 859-864.
[32] 袁利鹏, 卢倚群, 孙远明, 等. 盐酸克伦特罗抗体纯化及直接竞争ELISA方法研究[J]. 食品科学, 2013, 34(14): 227-231. DOI:10.7506/spkx1002-6630-201314046.
[33] 王玲, 管笛, 周欣, 等. 磁性化学发光酶免疫法检测猪肉中的氯霉素[J].食品科学, 2017, 38(10): 305-309. DOI:10.7506/spkx1002-6630-201710049.
Simultaneous Detection of Four Veterinary Drugs Using Immunochip Assay Based on Smart Phone
WU Ying, FAN Congcong, SU Xiaona, et al. Simultaneous detection of four veterinary drugs using immunochip assay based on smart phone[J]. Food Science, 2021, 42(2): 278-282. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20190927-334. http://www.spkx.net.cn