CN113376142A - 一种基于prlibs技术的蔬菜富硒元素检测方法 - Google Patents
一种基于prlibs技术的蔬菜富硒元素检测方法 Download PDFInfo
- 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
Links
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 title claims abstract description 152
- 229910052711 selenium Inorganic materials 0.000 title claims abstract description 151
- 239000011669 selenium Substances 0.000 title claims abstract description 151
- 235000013311 vegetables Nutrition 0.000 title claims abstract description 62
- 238000001514 detection method Methods 0.000 title claims abstract description 18
- 238000005516 engineering process Methods 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 28
- 230000010287 polarization Effects 0.000 claims abstract description 27
- 238000002310 reflectometry Methods 0.000 claims abstract description 20
- 238000002834 transmittance Methods 0.000 claims abstract description 9
- 230000003287 optical effect Effects 0.000 claims abstract description 7
- 239000002245 particle Substances 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 15
- 230000007704 transition Effects 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 12
- 239000000523 sample Substances 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 10
- 230000005284 excitation Effects 0.000 claims description 9
- 238000005192 partition Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 150000002500 ions Chemical class 0.000 claims description 4
- 102220530114 Linker for activation of T-cells family member 2_S43I_mutation Human genes 0.000 claims description 3
- 102220574555 Neuronal acetylcholine receptor subunit alpha-2_S41I_mutation Human genes 0.000 claims description 3
- 102220635540 Vacuolar protein sorting-associated protein 33A_S44I_mutation Human genes 0.000 claims description 3
- 230000005281 excited state Effects 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 102200143559 rs11550103 Human genes 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract 1
- 229940091258 selenium supplement Drugs 0.000 description 106
- 210000002381 plasma Anatomy 0.000 description 49
- 239000013307 optical fiber Substances 0.000 description 8
- 238000002536 laser-induced breakdown spectroscopy Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 235000002568 Capsicum frutescens Nutrition 0.000 description 3
- 240000008574 Capsicum frutescens Species 0.000 description 3
- 235000009811 Momordica charantia Nutrition 0.000 description 3
- 244000078912 Trichosanthes cucumerina Species 0.000 description 3
- 235000008322 Trichosanthes cucumerina Nutrition 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 238000005215 recombination Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000005461 Bremsstrahlung Effects 0.000 description 1
- BHPQYMZQTOCNFJ-UHFFFAOYSA-N Calcium cation Chemical compound [Ca+2] BHPQYMZQTOCNFJ-UHFFFAOYSA-N 0.000 description 1
- 235000002566 Capsicum Nutrition 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 239000006002 Pepper Substances 0.000 description 1
- 235000016761 Piper aduncum Nutrition 0.000 description 1
- 235000017804 Piper guineense Nutrition 0.000 description 1
- 244000203593 Piper nigrum Species 0.000 description 1
- 235000008184 Piper nigrum Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000001093 anti-cancer Effects 0.000 description 1
- 239000003963 antioxidant agent Substances 0.000 description 1
- 230000003078 antioxidant effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000009614 chemical analysis method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 125000003748 selenium group Chemical group *[Se]* 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/71—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
- G01N21/718—Laser 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技术的蔬菜富硒元素检测方法。
背景技术
硒元素是人体生命活动中许多抗氧化酶的必要组成部分，被誉为“生命的奇效元素”和“抗癌之王”，适量补硒有助于增强机体免疫力、延缓衰老和抗癌防癌。
在硒元素检测技术的研究中，常规化学分析方法具有稳定性好、精准度高的优点，但是设备昂贵且不易搬动。作为一种绿色无污染的测量手段，光学方法可以对富硒蔬菜进行现场快速测量。其中，激光诱导击穿光谱（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的计算模型为：
式中，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的计算模型为：
式中，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 计算模型为：
式中，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 计算公式为：
其中，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的计算模型为：
式中，c为光速，ε为介电常数，E为电矢量值，n为复折射率。
优选地，能量J的计算模型为：
式中，c为光速，θ i 为入射角度，n 1为第一介质复折射率，A 2为单位入射面积。
优选地，能量J (r) 的计算模型为：
式中，c为光速，θ i 为入射角度，n 1为第一介质复折射率，R 2为单位反射面积。
优选地，分量J (r) // 和J (r) ⊥ 的计算模型为：
式中，c为光速，θ i 为入射角度，n 1为第一介质复折射率，R // 2为单位平行反射面积，R ⊥ 2为单位垂直反射面积。
优选地，步骤S3包括富硒蔬菜表面硒元素等离子体的光强反射率R的计算模型为：
式中，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) 的计算模型为：
式中，c为光速，θ t 为折射角度，n 2为入射介质复折射率，T 2为单位折射面积。
优选地，分量J (t) //和J (t) ⊥ 的计算模型为：
式中，c为光速，θ t 为折射角度，n 2为入射介质复折射率，T // 2为单位平行折射面积，T ⊥ 2为单位垂直折射面积。
优选地，步骤S3包括富硒蔬菜表面硒元素等离子体的光强折射率T的计算模型为：
式中，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 。
优选地，离子的总浓度关系式为：
式中，q为子的配分函数，N j 为能级分布数，N是N j 的求和，g j 是j能级的统计权重，ε j 是j能级的激发能，k为玻尔兹曼常数，T为电离温度。
优选地，原子线强度I λ 的计算模型为：
式中，λ为跃迁波长，N s 为发射原子的数密度，A ji 为该线的跃迁几率，U s (T)为等离子体温度下该类粒子的配分函数，g j 是j能级的统计权重，ε j 是j能级的激发能，k为玻尔兹曼常数，T为电离温度。
优选地，信号强度I ji 的计算模型为：
式中，F是实验参数，包括仪器的接收效率和其他，在测量过程中需稳定各种实验条件，保证F值是常数，C s 是等离子体中粒子的含量，A ji 是从j能级跃迁到i能级的跃迁几率，g j 是j能级的统计权重，ε j 是j能级的激发能，k为玻尔兹曼常数，T为电离温度，U s (T)为等离子体温度下该类粒子的配分函数。
优选地，步骤S4包括的蔬菜中硒元素在不同入射角、方位角下的横波信号I (s) ji 计算模型为：
式中，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 计算公式为：
其中，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的计算模型为：
式中，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的计算模型为：
式中，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 计算模型为：
式中，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 计算公式为：
其中，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 ，得到蔬菜表面的硒元素分布状态。
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)
| 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的方法 |
-
2021
- 2021-06-21 CN CN202110688582.7A patent/CN113376142B/zh active Active
Patent Citations (4)
| 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)
| 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 |