CN113376142A - 一种基于prlibs技术的蔬菜富硒元素检测方法 - Google Patents

一种基于prlibs技术的蔬菜富硒元素检测方法 Download PDF

Info

Publication number
CN113376142A
CN113376142A CN202110688582.7A CN202110688582A CN113376142A CN 113376142 A CN113376142 A CN 113376142A CN 202110688582 A CN202110688582 A CN 202110688582A CN 113376142 A CN113376142 A CN 113376142A
Authority
CN
China
Prior art keywords
selenium
plasma
energy
incidence
vector
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.)
Granted
Application number
CN202110688582.7A
Other languages
English (en)
Other versions
CN113376142B (zh
Inventor
徐将
徐媛
姚明印
王晓
刘斌
黄杨生
曾敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangxi Agricultural University
Original Assignee
Jiangxi Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jiangxi Agricultural University filed Critical Jiangxi Agricultural University
Priority to CN202110688582.7A priority Critical patent/CN113376142B/zh
Publication of CN113376142A publication Critical patent/CN113376142A/zh
Application granted granted Critical
Publication of CN113376142B publication Critical patent/CN113376142B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/71Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
    • G01N21/718Laser microanalysis, i.e. with formation of sample plasma

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Plasma & Fusion (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明提供一种基于PRLIBS技术的蔬菜富硒元素检测方法,其特征在于,包括:利用光学透镜组合装置将脉冲激光器发射的低能量激光聚焦于富硒蔬菜表面,得到硒元素等离子体;通过偏振分束器和光谱仪采集硒元素等离子体的偏振光谱图;计算富硒蔬菜表面硒元素等离子体的的光强反射率R和透射率T;在偏振光谱图中提取硒元素峰值光强,根据反射透射模型计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji ;结合菲涅尔方程,根据横波信号I (s) ji 判断计算蔬菜中的富硒元素I se 。本发明的检测方法利用了等离子体横波较活跃的特征,可以在测量时得到不同角度下的反射率和透射率,从而实现蔬菜富硒元素的稳定检测,并且可以降低对激光能量的需求。

Description

一种基于PRLIBS技术的蔬菜富硒元素检测方法
技术领域
本发明属蔬菜富硒检测方法技术领域,具体涉及一种基于PRLIBS技术的蔬菜富硒元素检测方法。
背景技术
硒元素是人体生命活动中许多抗氧化酶的必要组成部分,被誉为“生命的奇效元素”和“抗癌之王”,适量补硒有助于增强机体免疫力、延缓衰老和抗癌防癌。
在硒元素检测技术的研究中,常规化学分析方法具有稳定性好、精准度高的优点,但是设备昂贵且不易搬动。作为一种绿色无污染的测量手段,光学方法可以对富硒蔬菜进行现场快速测量。其中,激光诱导击穿光谱(laser-induced breakdown spectroscopy,LIBS)技术主要依靠脉冲激光激发和灼烧物质产生等离子体,用光谱仪检测并获得元素含量信息。
偏振态是光子受电磁场影响产生的,在等离子体非平衡态时存在各向异性电子速度分布使得硒元素的轫致辐射、复合辐射都具有一定的线偏性。在等离子体复合过程中,处于局部热平衡状态的硒原子和电子的单位时间偏转粒子数不同,导致了辣椒中硒元素复合辐射的偏振化。因此,偏振分辨激光诱导击穿光谱(polarization-resolved laser-induced breakdown spectroscopy,PRLIBS)技术借助等离子体中的偏振信息,有助于提高检测数据的信噪比和信背比,增强硒含量光谱信息的稳定性。
发明内容
本发明旨在提供一种基于PRLIBS技术的蔬菜富硒元素检测方法,以实现蔬菜富硒元素的稳定检测,并且可以降低对激光能量的需求。
为实现上述目的,本发明采用以下技术方案:
本发明提供一种基于PRLIBS技术的蔬菜富硒元素检测方法,包括:
步骤S1:利用光学透镜组合装置将脉冲激光器发射的低能量激光聚焦于富硒蔬菜表面,得到硒元素等离子体;
步骤S2:通过偏振分束器和光谱仪采集步骤S1所述硒元素等离子体的偏振光谱图;
步骤S3:计算步骤S1所述硒元素等离子体的光强反射率R和透射率T
步骤S4:在步骤S2所述偏振光谱图中提取硒元素峰值光强,根据步骤S3所述光强反射率R和透射率T计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
步骤S5:结合菲涅尔方程,根据步骤S4所述横波信号I (s) ji 判断计算蔬菜中的富硒元素I se
优选地,步骤S3所述计算富硒蔬菜表面硒元素等离子体的光强反射率R,包括:
步骤S311:按照坡印廷矢量方程,计算得到步骤S1所述硒元素等离子体的光强系S
步骤S312:根据步骤S311所述光强系S计算得到信号波每秒入射到界面单位面积上的能量J
步骤S313:利用斯涅耳定律,根据步骤S312所述能量J计算得到反射波每秒从界面单位面积带走的能量J (r)
步骤S314:把步骤S313所述能量J (r) 分解成平行于入射面的分量J (r) // 和垂直于入射面的分量J (r)
步骤S315:结合菲涅尔方程,根据步骤S314所述两个分量J (r) // J (r) ,计算得到富硒蔬菜表面硒元素等离子体的光强反射率R
优选地,步骤S3包括富硒蔬菜表面硒元素等离子体的光强反射率R的计算模型为:
Figure 112036DEST_PATH_IMAGE001
式中,J (r) 为反射波每秒从界面单位面积带走的能量,J (r) // J (r) 矢量分解成平行于入射面的分量,J (r) J (r) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,R // R矢量分解成平行于入射面的分量,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角。
优选地,步骤S3所述计算富硒蔬菜表面等离子体的光强透射率T,包括:
步骤S321:利用斯涅耳定律,根据步骤S312所述能量J计算得到折射波每秒从界面单位面积带走的能量J (t)
步骤S322:把步骤S321所述的能量J (t) 分解成平行于入射面的分量J (t) //和垂直于入射面的分量J (t)
步骤S323:结合菲涅尔方程,根据S322所述分量J (t) //和分量J (t) ,计算得到富硒蔬菜表面硒元素等离子体的光强折射率T
优选地,步骤S3包括富硒蔬菜表面硒元素等离子体的光强折射率T的计算模型为:
Figure 105400DEST_PATH_IMAGE002
式中,J (t) 为折射波每秒从界面单位面积带走的能量,J (t) //J (t) 矢量分解成平行于入射面的分量,J (t) J (t) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,T // T矢量分解成平行与入射面的分量,T T矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角。
优选地,步骤S4所述计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji ,包括:
步骤S41:根据玻尔兹曼分布公式,得到激发态能级的粒子分布数与中性原子或者该元素离子的总浓度关系式;
步骤S42:利用S41所述总浓度关系式,得到激光等离子体中硒元素对应两个能级跃迁的原子线强度I λ
步骤S43:根据步骤S42所述原子线强度I λ ,综合光谱接收效率和探头接收角度,得到实测的信号强度I ji
步骤S44:结合步骤S315所述光强反射率R与S43所述信号强度I ji ,推导得到蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
优选地,步骤S4包括的蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji 计算模型为:
Figure 328571DEST_PATH_IMAGE003
式中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角,β为入射光垂直方向和入射面垂直方向之间的夹角。
优选地,步骤S5所述判断计算蔬菜中的富硒元素I se ,包括:
步骤S51:根据步骤S44所述横波信号I (s) ji ,结合布鲁斯特角度和方位角的实测值,得到蔬菜中的富硒元素I se
优选地,步骤S5包括的蔬菜中的富硒元素I se 计算公式为:
Figure 65583DEST_PATH_IMAGE004
其中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量。
优选地,所述基于PRLIBS技术的蔬菜富硒元素检测方法,还包括步骤S6:移动富硒蔬菜表面的检测位置,重复步骤S2至步骤S5,判断计算出富硒蔬菜表面不同位置的富硒元素I se
优选地,步骤S6包括:通过判断计算出富硒蔬菜表面不同位置的富硒元素I se ,得到蔬菜表面的硒元素分布状态。
本发明可实现以下有益效果:
本发明的一种基于PRLIBS技术的蔬菜富硒元素检测方法仅仅通过对传统激光诱导击穿光谱设备的改装,结合菲涅尔方程和散射光的偏振特性分解出了横波信号,用于识别蔬菜中的富硒元素;通过反射率和透射率模型,构建横波光强信息的特征参量,用来分析硒元素等离子体光强产生机理与空间光信号偏振信息分解方式之间的关系;通过提高硒元素检测数据的信噪比和信背比,从而实现了蔬菜硒元素检测的稳定性;此外,本发明通过特征优化降低了激光强度要求,可以通过测量检测等离子体所在位置,得到蔬菜表面的硒元素分布状态。
附图说明
图1是一种典型的硒元素检测光谱原理图。
图2是本发明所采用的偏振分辨激光诱导击穿光谱系统的结构示意图。
图3是本发明实施例的光谱传播路径示意图。
图4是本发明实施例的蔬菜富硒元素检测方法的技术路线图。
图5是本发明实施例所测量得富硒苦瓜中硒元素的检测光谱图。
图6是本发明实施例所测量得富硒辣椒中硒元素的检测光谱图。
图中:1-脉冲激光器,2-数据采集仪,201-穿孔镜,202-1#聚焦透镜,203-蔬菜表面,204-偏振分束器,205-1#反射镜,206-2#反射镜,207-2#聚焦透镜,208-3#聚焦透镜,209-1#光纤探头,210-2#光纤探头,3-数字脉冲延迟发生器,4-光谱仪,5-处理器。
具体实施方式
下面结合附图1-6,给出本发明的较佳实施例,并予以详细地描述,使能更好地理解本发明的功能、特点。
偏振分辨激光诱导击穿光谱(PRLIBS)是指通过超短脉冲激光聚焦样品表面形成等离子体,利用偏振分束器与多通道光谱仪对发光等离子体发射光谱进行采集分析,从而确定物质的成分及含量。偏振是电磁波矢量振动的空间分布对于传播方向失去了对称性的现象。因偏振特征参量值的改变与电磁波空间传播路径及碰撞媒介晶体排列结构等因素相关,且较传统检测方式具有安全高效、无需接触和无损的优势,作为一种新技术,在农产品检测领域被广泛应用。鉴于硒元素光谱信号的特点,可以利用偏振分辨激光诱导击穿光谱系统对蔬菜中硒元素进行监测,所以通过测量检测等离子体所在位置,得到蔬菜表面的硒元素分布状态。
本发明提供了一种基于PRLIBS技术的蔬菜富硒元素检测方法,包括以下步骤:
步骤S1:利用光学透镜组合装置将脉冲激光器发射的低能量激光聚焦于富硒蔬菜表面,得到硒元素等离子体;在激光作用下,蔬菜中的硒元素会以电离度较高的气体形式出现,气体在大气压力下放电,原子和分子密度较大,电子的自由行程较短,电子和重粒子之间频繁碰撞,将获得的动能较快地传递给重粒子,再以光谱的形式表示硒元素的具体含量;
步骤S2:通过偏振分束器和光谱仪采集步骤S1所述硒元素等离子体的偏振光谱图;
步骤S3:计算步骤S1所述硒元素等离子体的光强反射率R和透射率T
步骤S4:在步骤S2所述偏振光谱图中提取硒元素峰值光强,根据步骤S3所述光强反射率R和透射率T计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
步骤S5:结合菲涅尔方程,根据步骤S4所述横波信号I (s) ji 判断计算蔬菜中的富硒元素I se
优选地,步骤S3所述计算富硒蔬菜表面硒元素等离子体的光强反射率R,包括:
步骤S311:按照坡印廷矢量方程,计算得到步骤S1所述硒元素等离子体的光强系S
步骤S312:根据步骤S311所述光强系S计算得到信号波每秒入射到界面单位面积上的能量J
步骤S313:利用斯涅耳定律,根据步骤S312所述能量J计算得到反射波每秒从界面单位面积带走的能量J (r)
步骤S314:把步骤S313所述能量J (r) 分解成平行于入射面的分量J (r) // 和垂直于入射面的分量J (r)
步骤S315:结合菲涅尔方程,根据步骤S314所述两个分量J (r) // J (r) ,计算得到富硒蔬菜表面硒元素等离子体的光强反射率R
优选地,光强系S的计算模型为:
Figure 588968DEST_PATH_IMAGE005
式中,c为光速,ε为介电常数,E为电矢量值,n为复折射率。
优选地,能量J的计算模型为:
Figure 18813DEST_PATH_IMAGE006
式中,c为光速,θ i 为入射角度,n 1为第一介质复折射率,A 2为单位入射面积。
优选地,能量J (r) 的计算模型为:
Figure 712968DEST_PATH_IMAGE007
式中,c为光速,θ i 为入射角度,n 1为第一介质复折射率,R 2为单位反射面积。
优选地,分量J (r) // J (r) 的计算模型为:
Figure 988091DEST_PATH_IMAGE008
式中,c为光速,θ i 为入射角度,n 1为第一介质复折射率,R // 2为单位平行反射面积,R 2为单位垂直反射面积。
优选地,步骤S3包括富硒蔬菜表面硒元素等离子体的光强反射率R的计算模型为:
Figure 631562DEST_PATH_IMAGE001
式中,J (r) 为反射波每秒从界面单位面积带走的能量,J (r) // J (r) 矢量分解成平行于入射面的分量,J (r) J (r) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,R // R矢量分解成平行于入射面的分量,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角。
优选地,步骤S3所述计算富硒蔬菜表面等离子体的光强透射率T,包括:
步骤S321:利用斯涅耳定律,根据步骤S312所述能量J计算得到折射波每秒从界面单位面积带走的能量J (t)
步骤S322:把步骤S321所述的能量J (t) 分解成平行于入射面的分量J (t) //和垂直于入射面的分量J (t)
步骤S323:结合菲涅尔方程,根据S322所述分量J (t) //和分量J (t) ,计算得到富硒蔬菜表面硒元素等离子体的光强折射率T
优选地,能量J (t) 的计算模型为:
Figure 169991DEST_PATH_IMAGE009
式中,c为光速,θ t 为折射角度,n 2为入射介质复折射率,T 2为单位折射面积。
优选地,分量J (t) //J (t) 的计算模型为:
Figure 898913DEST_PATH_IMAGE010
式中,c为光速,θ t 为折射角度,n 2为入射介质复折射率,T // 2为单位平行折射面积,T 2为单位垂直折射面积。
优选地,步骤S3包括富硒蔬菜表面硒元素等离子体的光强折射率T的计算模型为:
Figure 977727DEST_PATH_IMAGE002
式中,J (t) 为折射波每秒从界面单位面积带走的能量,J (t) //J (t) 矢量分解成平行于入射面的分量,J (t) J (t) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,T // T矢量分解成平行与入射面的分量,T T矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角。
优选地,步骤S4所述计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji ,包括:
步骤S41:根据玻尔兹曼分布公式,得到激发态能级的粒子分布数与中性原子或者该元素离子的总浓度关系式;
步骤S42:利用S41所述总浓度关系式,得到激光等离子体中硒元素对应两个能级跃迁的原子线强度I λ
步骤S43:根据步骤S42所述原子线强度I λ ,综合光谱接收效率和探头接收角度,得到实测的信号强度I ji
步骤S44:结合步骤S315所述光强反射率R与S43所述信号强度I ji ,推导得到蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
优选地,离子的总浓度关系式为:
Figure 944546DEST_PATH_IMAGE011
式中,q为子的配分函数,N j 为能级分布数,NN j 的求和,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度。
优选地,原子线强度I λ 的计算模型为:
Figure 716193DEST_PATH_IMAGE012
式中,λ为跃迁波长,N s 为发射原子的数密度,A ji 为该线的跃迁几率,U s (T)为等离子体温度下该类粒子的配分函数,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度。
优选地,信号强度I ji 的计算模型为:
Figure 666832DEST_PATH_IMAGE013
式中,F是实验参数,包括仪器的接收效率和其他,在测量过程中需稳定各种实验条件,保证F值是常数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数。
优选地,步骤S4包括的蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji 计算模型为:
Figure 1867DEST_PATH_IMAGE003
式中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角,β为入射光垂直方向和入射面垂直方向之间的夹角。
优选地,步骤S5所述判断计算蔬菜中的富硒元素I se ,包括:
步骤S51:根据步骤S44所述横波信号I (s) ji ,结合布鲁斯特角度和方位角的实测值,得到蔬菜中的富硒元素I se
步骤S5包括的蔬菜中的富硒元素I se 计算公式为:
Figure 354351DEST_PATH_IMAGE004
其中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量。
优选地,本发明还包括步骤S6:移动富硒蔬菜表面的检测位置,重复步骤S2至步骤S5,判断计算出富硒蔬菜表面不同位置的富硒元素I se
优选地,步骤S6还可以包括:通过判断计算出富硒蔬菜表面不同位置的富硒元素I se ,得到蔬菜表面的硒元素分布状态。
优选地,本发明的PRLIBS技术的系统包括:脉冲激光器1、数据采集仪2、数字脉冲延迟发生器3、光谱仪4、处理器5;其中脉冲激光器1的中心波长为1064nm,能量分别为10mJ、36mJ、56mJ和78mJ;沿着光路的走向,数据采集仪2中依次分布着穿孔镜201,1#聚焦透镜202,蔬菜表面203,偏振分束器204,1#反射镜205,2#反射镜206,2#聚焦透镜207,3#聚焦透镜208,1#光纤探头209,2#光纤探头210;数字脉冲延迟发生器3的采集门宽度为5ns,采集延时步长为0.5μs;光谱仪4的波长范围为200-1050nm,分辨率为0.06-0.13nm;1#聚焦透镜202、2#聚焦透镜207、3#聚焦透镜208的焦距都为100mm。
优选地,其工作原理如下:使用时,为了得到有效的样品数据,对蔬菜表面203的每个监测点进行100次光谱信号采集;数字脉冲延迟发生器3控制脉冲激光器1以脉冲频率激发出激光光束,通过斜45°穿孔镜201后垂直照射在1#聚焦透镜202上,然后聚焦在蔬菜表面203的测量点处,产生高能量等离子体;等离子体跃迁光强信号经过斜135°穿孔镜201后,以56°入射角穿过偏振分束器204,得到纵波和横波,通过1#反射镜205、2#反射镜206与2#聚焦透镜207、3#聚焦透镜208组合光路后聚焦至1#光纤探头209、2#光纤探头210上进行信号采集;最后利用处理器5对采集信号进行优化处理。
本发明实施例的技术路线图如图4所示,对蔬菜中硒元素等离子体的相关物理特性进行归纳总结,构建偏振分光路径;基于路径方程,依托麦克斯韦方程和斯涅耳定律,建立硒元素的偏振分辨激光诱导击穿光谱特征函数;明确等离子体辐射偏振特性,结合特征函数,组成硒元素光学作用方程组;根据扰动分析、数值仿真、实验验证等手段,从蔬菜中硒元素等离子体的能量传输方程组中提炼特征峰信号;结合连续谱降噪处理、系统误差分析,形成硒元素偏振分辨激光诱导击穿光谱监测方法,通过测量检测等离子体所在位置,从而得到蔬菜表面的硒元素分布状态。
本发明实施例的实验结果如图5-图6所示,根据步骤S6得到富硒苦瓜和辣椒中硒元素的监测光谱图,如图5-图6所反映的,对富硒苦瓜和辣椒的表面经过多次击打后,得到的等离子体光谱信号中可以观察到明显的硒元素峰值信号,且可以判断得到硒元素的含量要比钙元素的含量少。
以上所述的,仅为本发明的较佳实施例,并非用以限定本发明的范围,本发明的上述实施例还可以做出各种变化,例如本发明不仅限于以上操作步骤,在不背离本发明原则的范围内可以作出适当调整。即凡是依据本发明申请的权力要求书内容所作的简单、等效变化与修饰,皆落入本发明专利的权利要求保护范围;本发明未详尽描述的均为常规技术内容。

Claims (3)

1.一种基于PRLIBS技术的蔬菜富硒元素检测方法,其特征在于包括:
步骤S1:利用光学透镜组合装置将脉冲激光器发射的低能量激光聚焦于富硒蔬菜表面,得到硒元素等离子体;
步骤S2:通过偏振分束器和光谱仪采集步骤S1所述硒元素等离子体的偏振光谱图;
步骤S3:计算步骤S1所述硒元素等离子体的光强反射率R和透射率T
步骤S4:在步骤S2所述偏振光谱图中提取硒元素峰值光强,根据步骤S3所述光强反射率R和透射率T计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
步骤S5:结合菲涅尔方程,根据步骤S4所述横波信号I (s) ji 判断计算蔬菜中的富硒元素I se
步骤S3所述计算富硒蔬菜表面硒元素等离子体的光强反射率R,包括:
步骤S311:按照坡印廷矢量方程,计算得到步骤S1所述硒元素等离子体的光强系S
步骤S312:根据步骤S311所述光强系S计算得到信号波每秒入射到界面单位面积上的能量J
步骤S313:利用斯涅耳定律,根据步骤S312所述能量J计算得到反射波每秒从界面单位面积带走的能量J (r)
步骤S314:把步骤S313所述能量J (r) 分解成平行于入射面的分量J (r) // 和垂直于入射面的分量J (r)
步骤S315:结合菲涅尔方程,根据步骤S314所述两个分量J (r) // J (r) ,计算得到富硒蔬菜表面硒元素等离子体的光强反射率R
步骤S3包括富硒蔬菜表面硒元素等离子体的光强反射率R的计算模型为:
Figure 667549DEST_PATH_IMAGE001
式中,J (r) 为反射波每秒从界面单位面积带走的能量,J (r) // J (r) 矢量分解成平行于入射面的分量,J (r) J (r) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,R // R矢量分解成平行于入射面的分量,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角;
步骤S3所述计算富硒蔬菜表面等离子体的光强透射率T,包括:
步骤S321:利用斯涅耳定律,根据步骤S312所述能量J计算得到折射波每秒从界面单位面积带走的能量J (t)
步骤S322:把步骤S321所述的能量J (t) 分解成平行于入射面的分量J (t) //和垂直于入射面的分量J (t)
步骤S323:结合菲涅尔方程,根据S322所述分量J (t) //和分量J (t) ,计算得到富硒蔬菜表面硒元素等离子体的光强折射率T
步骤S3包括富硒蔬菜表面硒元素等离子体的光强折射率T的计算模型为:
Figure 490011DEST_PATH_IMAGE002
式中,J (t) 为折射波每秒从界面单位面积带走的能量,J (t) //J (t) 矢量分解成平行于入射面的分量,J (t) J (t) 矢量分解成垂直于入射面的分量,J (i) 为入射波的能量,J (i) // J (i) 矢量分解成平行于入射面的分量,J (i) J (i) 矢量分解成垂直于入射面的分量,T // T矢量分解成平行与入射面的分量,T T矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角;
步骤S4所述计算出蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji ,包括:
步骤S41:根据玻尔兹曼分布公式,得到激发态能级的粒子分布数与中性原子或者该元素离子的总浓度关系式;
步骤S42:利用S41所述总浓度关系式,得到激光等离子体中硒元素对应两个能级跃迁的原子线强度I λ
步骤S43:根据步骤S42所述原子线强度I λ ,综合光谱接收效率和探头接收角度,得到实测的信号强度I ji
步骤S44:结合步骤S315所述光强反射率R与S43所述信号强度I ji ,推导得到蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji
步骤S4包括的蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji 计算模型为:
Figure 209574DEST_PATH_IMAGE003
式中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量,α i 为入射波电矢量与入射面的夹角,β为入射光垂直方向和入射面垂直方向之间的夹角;
步骤S5所述判断计算蔬菜中的富硒元素I se ,包括:
步骤S51:根据步骤S44所述横波信号I (s) ji ,结合布鲁斯特角度和方位角的实测值,得到蔬菜中的富硒元素I se
步骤S5包括的蔬菜中的富硒元素I se 计算公式为:
Figure 408474DEST_PATH_IMAGE004
其中,F是实验修正系数,C s 是等离子体中粒子的含量,A ji 是从j能级跃迁到i能级的跃迁几率,g j j能级的统计权重,ε j j能级的激发能,k为玻尔兹曼常数,T为电离温度,U s (T)为等离子体温度下该类粒子的配分函数,R R矢量分解成垂直于入射面的分量。
2.根据权利要求1所述的基于PRLIBS技术的蔬菜富硒元素检测方法,其特征在于,还包括步骤S6:移动富硒蔬菜表面的检测位置,重复步骤S2至步骤S5,判断计算出富硒蔬菜表面不同位置的富硒元素I se
3.根据权利要求2所述的基于PRLIBS技术的蔬菜富硒元素检测方法,其特征在于,所述步骤S6包括:通过判断计算出富硒蔬菜表面不同位置的富硒元素I se ,得到蔬菜表面的硒元素分布状态。
CN202110688582.7A 2021-06-21 2021-06-21 一种基于prlibs技术的蔬菜富硒元素检测方法 Active CN113376142B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110688582.7A CN113376142B (zh) 2021-06-21 2021-06-21 一种基于prlibs技术的蔬菜富硒元素检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110688582.7A CN113376142B (zh) 2021-06-21 2021-06-21 一种基于prlibs技术的蔬菜富硒元素检测方法

Publications (2)

Publication Number Publication Date
CN113376142A true CN113376142A (zh) 2021-09-10
CN113376142B CN113376142B (zh) 2023-12-22

Family

ID=77578213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110688582.7A Active CN113376142B (zh) 2021-06-21 2021-06-21 一种基于prlibs技术的蔬菜富硒元素检测方法

Country Status (1)

Country Link
CN (1) CN113376142B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2938066A1 (fr) * 2008-11-06 2010-05-07 Centre Nat Rech Scient Systeme et procede d'analyse quantitative de la composition elementaire de la matiere par spectroscopie du plasma induit par laser (libs)
CN204008464U (zh) * 2014-07-31 2014-12-10 邹平县好生镇卫生院 激光诱导击穿光谱测定中药材中微量元素硒含量的装置
WO2015104049A1 (en) * 2014-01-08 2015-07-16 Universidad Publica De Navarra Quantitative analysis method for analyzing the elemental composition of materials by means of libs technique
CN111398400A (zh) * 2020-03-31 2020-07-10 广西壮族自治区地质矿产测试研究中心 一种碱熔电感耦合等离子体质谱法测定Se、Te的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2938066A1 (fr) * 2008-11-06 2010-05-07 Centre Nat Rech Scient Systeme et procede d'analyse quantitative de la composition elementaire de la matiere par spectroscopie du plasma induit par laser (libs)
WO2015104049A1 (en) * 2014-01-08 2015-07-16 Universidad Publica De Navarra Quantitative analysis method for analyzing the elemental composition of materials by means of libs technique
CN204008464U (zh) * 2014-07-31 2014-12-10 邹平县好生镇卫生院 激光诱导击穿光谱测定中药材中微量元素硒含量的装置
CN111398400A (zh) * 2020-03-31 2020-07-10 广西壮族自治区地质矿产测试研究中心 一种碱熔电感耦合等离子体质谱法测定Se、Te的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王彩虹;黄林;杨晖;胡慧琴;陈添兵;刘木华;姚明印;: "蔬菜中Ca含量的LIBS快速检测研究", 江西农业大学学报, no. 03 *

Also Published As

Publication number Publication date
CN113376142B (zh) 2023-12-22

Similar Documents

Publication Publication Date Title
Faure et al. Observation of laser-pulse shortening in nonlinear plasma waves
Amiranoff et al. Electron acceleration in Nd-laser plasma beat-wave experiments
Tosa et al. High-order harmonic generation by chirped and self-guided femtosecond laser pulses. I. Spatial and spectral analysis
Giulietti et al. Production of ultracollimated bunches of multi-MeV electrons by 35 fs laser pulses propagating in exploding-foil plasmas
Goldman et al. Saturation of stimulated backscattered radiation in laser plasmas
CN103411903B (zh) 光丝阵列增强THz辐射方法
Young et al. Propagation of subpicosecond laser pulses through a fully ionized plasma
Watts et al. Measurements of relativistic self-phase-modulation in plasma
CN102798625B (zh) 粉状物质元素含量的激光诱导击穿光谱测量方法
Odagiri et al. Doubly excited states of molecular hydrogen as studied by coincident electron-energy-loss spectroscopy
Ke et al. Influence of ambient pressure on spatial–temporal evolution of local thermodynamic equilibrium for laser-induced plasma
Bowden et al. A Thomson scattering diagnostic system for measurement of electron properties of processing plasmas
Nagao et al. Nanosecond time-resolved Raman spectroscopy on phase transition of polytetrafluoroethylene under laser-driven shock compression
Burlakov et al. A differential-absorption lidar for ozone sensing in the upper atmosphere-lower stratosphere
CN113376142B (zh) 一种基于prlibs技术的蔬菜富硒元素检测方法
Alexander et al. Observation of recollision-based high-harmonic generation in liquid isopropanol and the role of electron scattering
Matsuoka et al. Focus optimization of relativistic self-focusing for anomalous laser penetration into overdense plasmas (super-penetration)
CN101975866B (zh) 电子波包扩散速度的测量装置和测量方法
CN103543131A (zh) 一种基于双脉冲和空间限制作用提高元素测量精度的方法
Apatin et al. Reactions induced in (CF 3 I) n clusters by femtosecond UV laser pulses
CN107505307B (zh) 全包围式椭圆球面镜光路激光诱导击穿光谱仪系统
Liu et al. Large focal depth of THz imaging system based on quasi-Bessel beams
Dangor et al. Generation of uniform plasmas for beat wave experiments
Bukreev et al. Excimer-laser-based lidar system for stratospheric and tropospheric ozone measurements
Zhang et al. Terahertz Liquid Photonics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210910

Assignee: Shangyou Guixin Agricultural Development Co.,Ltd.

Assignor: JIANGXI AGRICULTURAL University

Contract record no.: X2024980008559

Denomination of invention: A method for detecting selenium rich elements in vegetables based on PRLIBS technology

Granted publication date: 20231222

License type: Common License

Record date: 20240701

EE01 Entry into force of recordation of patent licensing contract