CN103969212B - 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 - Google Patents
利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 Download PDFInfo
- Publication number
- CN103969212B CN103969212B CN201410156255.7A CN201410156255A CN103969212B CN 103969212 B CN103969212 B CN 103969212B CN 201410156255 A CN201410156255 A CN 201410156255A CN 103969212 B CN103969212 B CN 103969212B
- Authority
- CN
- China
- Prior art keywords
- sample
- frequency range
- quantitative
- rice
- grain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000001514 detection method Methods 0.000 title claims abstract description 39
- 238000005516 engineering process Methods 0.000 title claims abstract description 25
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 title claims abstract 10
- 239000000447 pesticide residue Substances 0.000 title abstract description 12
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 34
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000001157 Fourier transform infrared spectrum Methods 0.000 claims abstract description 17
- 230000006835 compression Effects 0.000 claims abstract description 17
- 238000007906 compression Methods 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 239000004698 Polyethylene Substances 0.000 claims abstract description 9
- -1 polyethylene Polymers 0.000 claims abstract description 9
- 229920000573 polyethylene Polymers 0.000 claims abstract description 9
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 8
- 241000209094 Oryza Species 0.000 claims description 25
- 235000007164 Oryza sativa Nutrition 0.000 claims description 25
- 235000009566 rice Nutrition 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 23
- CVXBEEMKQHEXEN-UHFFFAOYSA-N carbaryl Chemical compound C1=CC=C2C(OC(=O)NC)=CC=CC2=C1 CVXBEEMKQHEXEN-UHFFFAOYSA-N 0.000 claims description 17
- 230000001419 dependent effect Effects 0.000 claims description 15
- 239000013598 vector Substances 0.000 claims description 12
- 238000002790 cross-validation Methods 0.000 claims description 10
- 238000012417 linear regression Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000227 grinding Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 239000000575 pesticide Substances 0.000 abstract description 22
- 235000013339 cereals Nutrition 0.000 description 25
- 239000000843 powder Substances 0.000 description 12
- YCIHPQHVWDULOY-FMZCEJRJSA-N (4s,4as,5as,6s,12ar)-4-(dimethylamino)-1,6,10,11,12a-pentahydroxy-6-methyl-3,12-dioxo-4,4a,5,5a-tetrahydrotetracene-2-carboxamide;hydrochloride Chemical compound Cl.C1=CC=C2[C@](O)(C)[C@H]3C[C@H]4[C@H](N(C)C)C(=O)C(C(N)=O)=C(O)[C@@]4(O)C(=O)C3=C(O)C2=C1O YCIHPQHVWDULOY-FMZCEJRJSA-N 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 6
- 230000003115 biocidal effect Effects 0.000 description 4
- 238000002329 infrared spectrum Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 239000000470 constituent Substances 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 235000013305 food Nutrition 0.000 description 3
- 229920001903 high density polyethylene Polymers 0.000 description 3
- 239000004700 high-density polyethylene Substances 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- KXDHJXZQYSOELW-UHFFFAOYSA-M Carbamate Chemical group NC([O-])=O KXDHJXZQYSOELW-UHFFFAOYSA-M 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000004570 mortar (masonry) Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- VYZAMTAEIAYCRO-UHFFFAOYSA-N Chromium Chemical compound [Cr] VYZAMTAEIAYCRO-UHFFFAOYSA-N 0.000 description 1
- 231100000703 Maximum Residue Limit Toxicity 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 208000031320 Teratogenesis Diseases 0.000 description 1
- 239000004098 Tetracycline Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 150000001336 alkenes Chemical class 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000004817 gas chromatography Methods 0.000 description 1
- 238000004128 high performance liquid chromatography Methods 0.000 description 1
- 238000003859 hyphenated technique Methods 0.000 description 1
- 238000002649 immunization Methods 0.000 description 1
- 230000003053 immunization Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002207 metabolite Substances 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 description 1
- 229960002180 tetracycline Drugs 0.000 description 1
- 229930101283 tetracycline Natural products 0.000 description 1
- 235000019364 tetracycline Nutrition 0.000 description 1
- 150000003522 tetracyclines Chemical class 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 235000012431 wafers Nutrition 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,通过将待测粮食样品研磨后压片得到的待测粮食压片样品,直接采用太赫兹频段FTIR光谱***对其进行测试,得到0‑20.4THz频段的吸收光谱,同时将目标农药与聚乙烯的混合压片样品采用相同的方式测试,以确定所述目标农药的特征波段,并将所述待测粮食压片样品和所述混合压片样品在所述特征波段下的吸收光谱随机划分为训练集样本吸收光谱和验证集样本吸收光谱,用偏最小二乘回归方法建立定量分析模型,获得各所述待测粮食样品的定量检测值,本发明所述方法能够真实、有效地实现对粮食中农药残留进行快速准确的定量检测,所述定量分析模型的预测相关系数(Rv)的平均值高达0.9927。
Description
技术领域
本发明涉及一种利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,属于农药检测技术领域。
背景技术
随着现代农业的快速发展,农业及农产品对农药的需求和依赖日益增长。加之一些从业人员缺乏农药知识,长期大量使用农药甚至滥用,使得农药对环境和人类健康造成了极大的影响和危害。对于农作物而言,农药施用到农作物之后,农作物内往往会有少量残留,长时间摄食残留农药会严重影响人体健康,亦会对身体和生态环境造成重大危害。特别是近年来,食物中农药残留超标而引起的致畸致残及中毒事件等也越来越受到关注,农药痕迹残留及其代谢物在土壤、水及农产品中均有发现。因而,在大力倡导科学使用农药的前提下,如何快速、准确的检测因农药引起的食品、环境的污染等问题就显得极为重要和迫切。
世界各国对农产品及食品中农药最大残留量均有相应检测方法和标准,现有的农药检测方法主要有气相色谱法、高效液相色谱法、气相色谱与质谱联用技术、液相色谱与质谱联用技术、免疫法、生物传感器法等,但上述方法均存在着样品前处理复杂,检测时间长,检测成本高,对检测人员要求较高,无法在线检测等不足。
远红外光谱区是指波长在0.025-1mm,频率在0.3-12THz之间的电磁波。在红外区域中是开拓最晚的一个波段,长期以来农残的红外光谱检测主要采取近红外或中红外光谱仪。近年来,随着科学技术的发展,傅里叶变换红外光谱(FTIR)技术已变为主流的红外光谱技术,分束器的使用已将光谱范围覆盖到远红外区域。该技术通过傅里叶去卷积处理将复杂的样品IR重叠谱峰有效分开,降低了峰宽,提高了分辨率;计算功能的显著增强减少了干扰以及杂质造成的信号漂移,从而提高了灵敏度、信噪比,扩大了波长范围,缩短扫描时间。太赫兹辐射(也称“THz辐射”)是指频率在0.1THz-10THz,波长在0.03-3mm之间的电磁波,其波段位于微波和红外线之间,是宏观电子学向微观光子学过渡的区域,在电磁波频谱中占有很特殊的位置。许多极性大分子在振动能级间的跃迁正好处于太赫兹频率范围,因此,生物分子的太赫兹光谱可以反映由分子内或分子间集体振动和晶格振动引起的低频振动膜的本征特性。太赫兹电磁波具有较低的光子能量,在进行样品检测时,不会产生有害的光致电离,是一种有效的无损探测方法。太赫兹光谱技术可作为近红外光谱技术的有益补充用于药品检测领域。
中国专利文献CN103472032A公开了一种利用太赫兹时域光谱技术检测盐酸四环素的方法,其包括如下步骤:(1)将盐酸四环素粉末与高密度聚乙烯粉末以不同比例混合研磨,用压片机一一压成圆盘形薄片,得到含有不同质量百分比的盐酸四环素压片;(2)在0.1-3.5THz频段范围内,用太赫兹时域光谱***一一采集不同质量百分比的盐酸四环素压片的太赫兹时域光谱;(3)以氮气的时域波形作为参考信号,以盐酸四环素压片的时域波形作为样品信号,分别进行傅里叶变换,得到两种信号的频域分布,利用公式得到压片样本的吸收系数和折射率;(4)根据各个压片样本的质量百分比及其对应的吸收谱建立校正模型,进行校正模型的验证与评价。上述方法虽然能够实现对食品中抗生素盐酸四环素进行定性和定量检测,但是在样品制备过程中,是将待检测的抗生素盐酸四环素和高密度聚乙烯粉末混合压片样品,利用聚乙烯为“背景”检测抗生素的含量,因而该方法最终检测的是存在于聚乙烯中抗生素的含量的相对数值,这就导致与实际的检测样本相差甚远,影响了检测样本的真实性;此外,引入新的无关物质聚乙烯,不可避免将会对样本检测的准确度造成影响,同时还浪费了药品。
发明内容
本发明所要解决的技术问题在于提供一种利用太赫兹频段FTIR技术结合化学计量学直接对粮食中农药残留进行定量检测的方法。
为解决上述技术问题,本发明是通过以下技术方案实现的:
本发明提供一种利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,其包括如下步骤:
(1)取待测粮食样品研磨后压片,即得待测粮食压片样品;
(2)取欲检测的目标农药样品,按照质量比1:4的比例与聚乙烯进行混合,并研磨后压片,即得含有目标农药与聚乙烯的混合压片样品;
(3)利用太赫兹频段FTIR光谱***对所述待测粮食压片样品和所述混合压片样品进行逐一测试,得到各样品在0-20.4THz频段内的吸收光谱,并根据所述混合压片样品的吸收光谱,选定吸收光谱的重现性较好的波段区间作为特征波段,并将所述待测粮食压片样品和所述混合压片样品在所述特征波段下的吸收光谱随机划分为训练集样本吸收光谱和验证集样本吸收光谱;
(4)利用偏最小二乘回归方法建立所述训练集样本吸收光谱和所述验证集样本吸收光谱的定量分析模型,获得各所述待测粮食样品的定量检测值。
所述农药为氨基甲酸酯类农药。
所述的氨基甲酸酯类农药为西维因。
所述步骤(3)中,所述特征波段为1.8-6.3THz。
所述步骤(3)中,所述FTIR的测试条件是:温度为25℃,波数范围为0-680cm-1,分辨率为4cm-1,扫描次数为200次。
所述步骤(3)中,所述训练集样本和所述验证集样本吸收光谱的划分利用留N交叉验证法进行;
所述留N交叉验证法的具体计算过程如下:
(a)假设数据集有n个样本,将样本随机排序;
(b)将该数据集划分成相等大小的m组(m=n/N);
(c)以第一组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(d)以第二组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(e)以此类推,共进行m次,直至所有组被抽出一次并进行预测为止。
所述步骤(3)中,每个样品重复测量3次。
所述待测粮食压片样品和混合压片样品均为直径13mm、厚度0.15mm的圆薄片。
所述步骤(4)中,所述偏最小二乘回归方法的具体计算原理如下:
偏最小二乘方法是建立在自变量X和因变量Y矩阵基础上的模型,通过建立自变量的潜变量关于因变量的潜变量的线性回归模型,进而反应自变量与因变量之间的关系;
X=TPT+E=∑tapa t;
Y=UQT+F=∑uaqa t;
式中:T为X的得分矩阵,P为X的载荷矩阵,E为X的残差矩阵,ta为得分向量,pa为相应的载荷向量,U为Y的得分矩阵,Q为Y的载荷矩阵,F为Y的残差矩阵,ua为得分向量,qa为相应的载荷向量;
偏最小二乘回归分别在X和Y中提取各自的潜变量,它们分别为自变量与因变量的线形组合,同时应满足以下条件:
(a)两组潜变量分别最大程度地承载自变量和因变量的变异信息;
(b)二者之间的协方差最大化。
所述偏最小二乘回归是一个以迭代法逐步提取成分的方法,在迭代计算中互相利用对方的信息,每一次迭代不断根据X、Y的剩余信息调整ta、ua进行第二轮的成分提取,直到残余矩阵中的元素绝对值近似为零,算法停止。由此得到的系数能更好的反映X和Y的关系。
本发明的上述技术方案相比现有技术具有以下优点:
(1)本发明所述的利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,通过将待测粮食样品研磨后压片得到的待测粮食压片样品,直接采用太赫兹频段FTIR光谱***对其进行测试,得到0-20.4THz频段的吸收光谱,同时将目标农药与聚乙烯的混合压片样品采用相同的方式测试,以确定所述目标农药的特征波段,并将所述待测粮食压片样品和所述混合压片样品在所述特征波段下的吸收光谱随机划分为训练集样本吸收光谱和验证集样本吸收光谱,用偏最小二乘回归方法建立定量分析模型,获得各所述待测粮食样品的定量检测值,较之现有技术中已公开的利用太赫兹光谱检测农药残留的方式而言,本发明所述方法用于检测的样品中不需掺入其他任何物质,样品制备简单,不需进行任何预处理,能够真实、有效地实现对粮食中农药残留进行快速准确的定量检测,所述定量分析模型的预测相关系数(Rv)的平均值高达0.9927,即所述方法的检测准确率高达99%以上。
(2)本发明所述的利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,在所述步骤(3)中采用留N交叉验证法对所述特征波段进行划分得到所述训练集样本和所述验证集样本吸收光谱,所述留N交叉验证法是检验模型稳健性的一种方法,在计算过程中,首先将样本随机排序,之后每次从数据集中抽出N个样本,用剩余的样本建模并预测被抽出的N个样本,保证每个样本用于且仅用于1次预测,保证了不同质量分数的压片样本在训练集和验证集中以同样的比例出现,从而实现对所建模型预测能力的无偏评价,使鉴定模型更可靠,分析结果更具有统计学意义。
(3)本发明所述的利用太赫兹频段FTIR技术定量检测粮食中农药残留的方法,在所述步骤(4)中,对所述训练集样本吸收光谱和所述验证集样本吸收光谱采用偏最小二乘回归方法建立定量分析模型,并采用校正均方根误差(RMSEC)、预测均方根误差(RMSEP)、预测相关系数(Rv)作为模型性能评判的依据,RMSEC、RMSEP越小,Rv越大,模型越好。
附图说明
为了使本发明的内容更容易被清楚的理解,下面结合附图,对本发明作进一步详细的说明,其中,
图1是本发明实施例1所述待测粮食样品在1.8-6.3THz特征波段区间的吸收光谱图;
图2是本发明所述预测均方误差与主成分数的关系图;
图3是本发明所述模型预测值与实验值之间的关系图。
具体实施方式
本实施例提供一种利用太赫兹频段FTIR技术结合化学计量学定量检测大米样品中西维因含量,其中,所述“大米样品”采用本身不含有任何农药残留的大米(吉林市家福米业有限公司生产)添加已知量西维因(北京北纳创联生物技术研究院提供)后制备得到,并将所述“大米样品”作为盲样采用本发明方法进行检测,具体步骤如下:
(1)将所述大米放入粉碎机中粉碎,过筛(100目),放入烘箱中烘干,得到大米粉末;将大米粉末与西维因粉末按照不同质量比例混合,转移至玛瑙研钵中进一步研磨得到所述“大米样品”细粉,其中,所述“大米样品”中西维因质量百分含量依次为0%,1.0%,1.5%,2.0%,2.5%,3.0%,3.5%,4.0%,4.5%,5.0%,5.5%,6.0%,6.5%,7.0%,7.5%,8.0%,8.5%,9.0%,9.5%,10.0%,11.0%,12.0%,13.0%,14.0%,16.0%,18.0%,20.0%;
称取上述每种所述“大米样品”细粉约24mg,置于Specac公司的模具中,用压片机在7t的压力下保持3-4min,从而将所述“大米样品”细粉压制成直径13mm、厚度0.15mm的圆薄片,即得含有不同质量百分比西维因的待测粮食压片样品,所述待测粮食样品两表面平行、表面光滑且没有裂缝;
(2)将西维因粉末与聚乙烯粉末按照质量比1:4混合,转移至玛瑙研钵中进一步研磨,采用与步骤(1)相同的方法压制得到所述混合压片样品;
(3)用太赫兹频段FTIR光谱***(Bruker公司,VERTEX/80V)对所述待测粮食压片样品和所述混合压片样品(共28个压片)进行逐一测试,每个压片重复测量3次得到0-20.4THz频段的吸收光谱,取平均值;根据所述混合压片样品的测试结果,选取1.8-6.3THz(即60-210cm-1)作为特征波段,在所述特征波段区间内吸收光谱的重现性较好;其中,所述FTIR的测试条件是:温度为25℃,波数范围为0-680cm-1,分辨率为4cm-1,扫描次数为200次;
如图1所示为其中6个不同西维因质量百分含量的所述待测粮食压片样品在1.8-6.3THz特征波段区间的吸收光谱图,6个不同西维因质量百分含量分别为1.5%、3.0%、6.0%、9.0%、11.0%、16.0%;
对于1.8-6.3THz特征波段的吸收光谱,采用留N交叉验证法将其划分为训练集样本吸收光谱和验证集样本吸收光谱,这里N值取7;于是将28个样本随机分为4组,每次均有21个样本作为训练集,7个样本作为验证集,重复4次,保证每个样本仅作为1次验证集,且每个样本作为训练集3次;
其中,所述留N交叉验证法,具体计算过程的说明如下:
(a)假设数据集有n个样本,将样本随机排序;
(b)将该数据集划分成相等大小的m组(m=n/N);
(c)以第一组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(d)以第二组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(e)以此类推,共进行m次,直至所有组被抽出一次并进行预测为止;
(4)对所述训练集样本和所述验证集样本光谱采用偏最小二乘回归方法建立定量分析模型,其具体原理描述如下:
偏最小二乘方法是建立在自变量X和因变量Y矩阵基础上的模型,通过建立自变量的潜变量关于因变量的潜变量的线性回归模型,进而反应自变量与因变量之间的关系;
X=TPT+E=∑tapa t;
Y=UQT+F=∑uaqa t;
式中:T为X的得分矩阵,P为X的载荷矩阵,E为X的残差矩阵,ta为得分向量,pa为相应的载荷向量,U为Y的得分矩阵,Q为Y的载荷矩阵,F为Y的残差矩阵,ua为得分向量,qa为相应的载荷向量;
偏最小二乘回归分别在X和Y中提取各自的潜变量,它们分别为自变量与因变量的线形组合,同时应满足以下条件:
(a)两组潜变量分别最大程度地承载自变量和因变量的变异信息;
(b)二者之间的协方差最大化;
偏最小二乘回归是一个以迭代法逐步提取成分的方法。在迭代计算中互相利用对方的信息,每一次迭代不断根据X、Y的剩余信息调整ta、ua进行第二轮的成分提取,直到残余矩阵中的元素绝对值近似为零,算法停止。由此得到的系数能更好的反映X和Y的关系。
如图2所示为预测方根误差与主成分数的关系图;根据图2选择主成分数为2,用留N交叉验证法对所建模型进行验证,将校正均方根误差(RMSEC)、预测均方根误差(RMSEV)、预测相关系数(Rv)作为模型性能评判的依据,结果如表1所示。
如图3所示为所述模型预测值与实验值之间的关系图,说明模型的预测值与实验值很接近,预测相关系数可达到0.9874。
表1-留N交叉验证法结合偏最小二乘回归建立分析模型的性能结果
从表中数据可以看出,上述分析模型的RMSEC、RMSEP值比较小,Rv较大,从而说明本发明方法建立的模型可靠可行,可以用来对大米样品中西维因进行定量检测。
采用上述描述的方法,不仅可以对大米中的西维因进行定量检测,也可以推广到其它种类的粮食样品和其它农药的检测。在其它种类粮食的样品的检测中,亦需将样品制备成粉末压片。所检测的农药应当在太赫兹频段内具有特征吸收。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (7)
1.一种利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,包括如下步骤:
(1)取待测大米样品研磨后压片,即得待测大米压片样品;
(2)取欲检测的目标西维因样品,按照质量比1:4的比例与聚乙烯进行混合,并研磨后压片,即得含有目标西维因与聚乙烯的混合压片样品;
(3)利用太赫兹频段FTIR光谱***对所述待测大米压片样品和所述混合压片样品进行逐一测试,得到各样品在0-20.4THz频段内的吸收光谱,并根据所述混合压片样品的吸收光谱,选定吸收光谱的重现性较好的波段区间作为特征波段,并将所述待测大米压片样品和所述混合压片样品在所述特征波段下的吸收光谱随机划分为训练集样本吸收光谱和验证集样本吸收光谱;
(4)利用偏最小二乘回归方法建立所述训练集样本吸收光谱和所述验证集样本吸收光谱的定量分析模型,获得各所述待测大米样品的定量检测值。
2.根据权利要求1所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述步骤(3)中,所述特征波段为1.8-6.3THz。
3.根据权利要求1或2所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述步骤(3)中,所述FTIR的测试条件是:温度为25℃,波数范围为0-680cm-1,分辨率为4cm-1,扫描次数为200次。
4.根据权利要求3所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述步骤(3)中,所述训练集样本和所述验证集样本吸收光谱的划分利用留N交叉验证法进行;
所述留N交叉验证法的具体计算过程如下:
(a)假设数据集有n个样本,将样本随机排序;
(b)将该数据集划分成相等大小的m组,m=n/N;
(c)以第一组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(d)以第二组为验证集,剩余m-1组为训练集,即用n-N个训练集样本建立分析模型,并预测被抽出来的N个样本;
(e)以此类推,共进行m次,直至所有组被抽出一次并进行预测为止。
5.根据权利要求1、2或4所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述步骤(3)中,每个样品重复测量3次。
6.根据权利要求5所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述待测大米压片样品和混合压片样品均为直径13mm、厚度0.15mm的圆薄片。
7.根据权利要求1、2、4或6所述的利用太赫兹频段FTIR技术定量检测大米中西维因残留的方法,其特征在于,所述步骤(4)中,所述偏最小二乘回归方法的具体计算原理如下:
偏最小二乘方法是建立在自变量X和因变量Y矩阵基础上的模型,通过建立自变量的潜变量关于因变量的潜变量的线性回归模型,进而反应自变量与因变量之间的关系;
X=TPT+E=∑tapa t;
Y=UQT+F=∑uaqa t;
式中:T为X的得分矩阵,P为X的载荷矩阵,E为X的残差矩阵,ta为得分向量,pa为相应的载荷向量,U为Y的得分矩阵,Q为Y的载荷矩阵,F为Y的残差矩阵,ua为得分向量,qa为相应的载荷向量;
偏最小二乘回归分别在X和Y中提取各自的潜变量,它们分别为自变量与因变量的线形组合,同时应满足以下条件:
(a)两组潜变量分别最大程度地承载自变量和因变量的变异信息;
(b)二者之间的协方差最大化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410156255.7A CN103969212B (zh) | 2014-04-17 | 2014-04-17 | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410156255.7A CN103969212B (zh) | 2014-04-17 | 2014-04-17 | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103969212A CN103969212A (zh) | 2014-08-06 |
CN103969212B true CN103969212B (zh) | 2016-08-17 |
Family
ID=51238978
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410156255.7A Expired - Fee Related CN103969212B (zh) | 2014-04-17 | 2014-04-17 | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103969212B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104237157A (zh) * | 2014-09-28 | 2014-12-24 | 首都师范大学 | 应用太赫兹时域光谱技术检测粮食中氨基酸含量的方法 |
CN104330384B (zh) * | 2014-11-14 | 2017-10-20 | 首都师范大学 | 应用太赫兹频域光谱技术检测粮食中氨基酸含量的方法 |
CN104406923A (zh) * | 2014-11-14 | 2015-03-11 | 首都师范大学 | 一种利用THz-TDS技术定量检测粮食中氨基酸含量的方法 |
CN108195791A (zh) * | 2017-12-18 | 2018-06-22 | 武汉中科志康生物科技有限公司 | 一种基于太核兹技术检测农药残留和非法添加剂的方法 |
CN108645811A (zh) * | 2018-04-18 | 2018-10-12 | 上海理工大学 | 一种采用太赫兹技术检测中草药三七的方法 |
CN113655023B (zh) * | 2021-07-15 | 2023-09-05 | 上海理工大学 | 一种检测低浓度莠去津快速检测的方法 |
CN114112979A (zh) * | 2021-12-02 | 2022-03-01 | 南京林业大学 | 一种基于太赫兹光谱的农产品中农药残留定量检测方法 |
CN117929347B (zh) * | 2023-12-14 | 2024-06-07 | 珠海康龙源生物科技有限公司 | 一种中药红花的农残快速检测方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590135A (zh) * | 2012-03-02 | 2012-07-18 | 中国计量学院 | 一种基于最小二乘支持向量机的除草剂鉴别方法 |
CN102930158A (zh) * | 2012-10-31 | 2013-02-13 | 哈尔滨工业大学 | 基于偏最小二乘的变量选择方法 |
CN103487399A (zh) * | 2013-09-24 | 2014-01-01 | 浙江大学 | 一种定量检测盐酸四环素溶液的太赫兹时域光谱方法 |
-
2014
- 2014-04-17 CN CN201410156255.7A patent/CN103969212B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590135A (zh) * | 2012-03-02 | 2012-07-18 | 中国计量学院 | 一种基于最小二乘支持向量机的除草剂鉴别方法 |
CN102930158A (zh) * | 2012-10-31 | 2013-02-13 | 哈尔滨工业大学 | 基于偏最小二乘的变量选择方法 |
CN103487399A (zh) * | 2013-09-24 | 2014-01-01 | 浙江大学 | 一种定量检测盐酸四环素溶液的太赫兹时域光谱方法 |
Non-Patent Citations (5)
Title |
---|
Qualitative and Quantitative Detection of Pesticides With Terahertz Time-Domain Spectroscopy;Yuefang Hua and Hongjian Zhang;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;20100731;第58卷(第7期);2064-2070 * |
Qualitative and quantitative identification of nitrofen in terahertz region;Wang Qiang,Ma Ye-Hao;《Chemometrics and Intelligent Laboratory Systems》;20130604;第127卷;43-48 * |
基于太赫兹时域光谱技术的农药定性和定量分析;花月芳;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20110715(第7期);第15-16页,2.3.2.1 固体样品制备,第29-31页,3.2.3.2 偏最小二乘PLS,第53页,4.2.1.3理论模拟与太赫兹光谱解析,第65页,4.2.4 单成分农药混合物光谱分析,第83页,4.3.1.2 吡虫啉与聚乙烯、面粉混合物的定量分析,表4.4 * |
基于太赫兹时域光谱技术的农药残留检测方法;曹丙花等;《红外与毫米波学报》;20081231;第27卷(第6期);429-432 * |
菊酯农药的太赫兹时域光谱定性和定量检测;郝国徽等;《光谱学与光谱分析》;20120531;第32卷(第5期);1184-1188 * |
Also Published As
Publication number | Publication date |
---|---|
CN103969212A (zh) | 2014-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103969212B (zh) | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 | |
CN103969214B (zh) | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 | |
Nakajima et al. | An overview of and issues with sky radiometer technology and SKYNET | |
CN103364362B (zh) | 一种利用THz-TDS结合化学计量学鉴定中草药的方法 | |
CN104181122A (zh) | 利用太赫兹时域光谱技术检测粮食中农药含量的方法 | |
CN104215599A (zh) | 一种利用THz-TDS技术定量检测粮食中农药残留的方法 | |
CN103196838B (zh) | 一种海岸河口富营养化高光谱遥感监测方法 | |
CN103091274B (zh) | 近红外漫反射光谱测定注射用丹参多酚酸水分含量的方法 | |
CN105486655A (zh) | 基于红外光谱智能鉴定模型的土壤有机质快速检测方法 | |
CN104749132A (zh) | 一种测定面粉中偶氮甲酰胺的含量的方法 | |
CN104406923A (zh) | 一种利用THz-TDS技术定量检测粮食中氨基酸含量的方法 | |
CN104849233A (zh) | 一种检测谷物新陈度的方法及装置 | |
Hu et al. | Optimization of quantitative detection model for benzoic acid in wheat flour based on CARS variable selection and THz spectroscopy | |
CN103487395A (zh) | 一种重楼药材多指标快速检测方法 | |
CN105784628A (zh) | 一种利用中红外光谱检测土壤有机质化学组成的方法 | |
CN105758819A (zh) | 一种利用近红外光谱检测土壤的有机组分的方法 | |
CN105136737A (zh) | 基于近红外光谱的馒头中马铃薯全粉含量快速测定的方法 | |
CN107064057A (zh) | 一种低含水率油田污泥含油量的快速测试方法 | |
CN105486663B (zh) | 一种利用近红外光谱检测土壤的稳定碳同位素比值的方法 | |
CN107505179A (zh) | 一种土壤预处理与养分近红外光谱检测方法 | |
CN104297202B (zh) | 利用THz‑TDS频域谱定量检测粮食中农药残留的方法 | |
CN106990056A (zh) | 一种土壤全氮光谱估算模型校正样本集构建方法 | |
CN106126879B (zh) | 一种基于稀疏表示技术的土壤近红外光谱分析预测方法 | |
CN114169165A (zh) | 一种三波段植被指数估算镉胁迫下水稻叶绿素的模型方法 | |
CN103105369B (zh) | 液态物质光谱基线校正定量分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 |
|
CF01 | Termination of patent right due to non-payment of annual fee |