CN115494039A - 基于激光诱导荧光的氦原子密度反演方法 - Google Patents
基于激光诱导荧光的氦原子密度反演方法 Download PDFInfo
- Publication number
- CN115494039A CN115494039A CN202211138153.3A CN202211138153A CN115494039A CN 115494039 A CN115494039 A CN 115494039A CN 202211138153 A CN202211138153 A CN 202211138153A CN 115494039 A CN115494039 A CN 115494039A
- Authority
- CN
- China
- Prior art keywords
- energy level
- laser
- density
- level
- atomic density
- 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
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 title claims abstract description 67
- 229910052734 helium Inorganic materials 0.000 title claims abstract description 55
- 239000001307 helium Substances 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000001499 laser induced fluorescence spectroscopy Methods 0.000 title claims abstract description 29
- 238000001514 detection method Methods 0.000 claims abstract description 38
- 238000012545 processing Methods 0.000 claims abstract description 31
- 230000005283 ground state Effects 0.000 claims abstract description 21
- 230000006698 induction Effects 0.000 claims abstract description 17
- 230000007935 neutral effect Effects 0.000 claims abstract description 17
- 230000007704 transition Effects 0.000 claims description 43
- 230000003287 optical effect Effects 0.000 claims description 13
- 238000002834 transmittance Methods 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 5
- 239000004065 semiconductor Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000009191 jumping Effects 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 description 9
- 238000000862 absorption spectrum Methods 0.000 description 5
- 239000003245 coal Substances 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 3
- 230000008030 elimination Effects 0.000 description 3
- 238000003379 elimination reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004927 fusion Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001739 density measurement Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002536 laser-induced breakdown spectroscopy Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000011425 standardization method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction 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/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/6402—Atomic fluorescence; Laser induced fluorescence
- G01N21/6404—Atomic fluorescence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Optics & Photonics (AREA)
- Probability & Statistics with Applications (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明提供基于激光诱导荧光的氦原子密度反演方法包括:设定氦原子单态能级,处理第一单位时间探测光强;求取c能级的原子密度;考虑激光诱导,设中性氦吸收预设频率的激光已进行激光诱导,处理单位时间探测到的荧光光强;获取探测光强差值与单位时间探测光强的比值;利用激光器对氦原子进行激光诱导,将氦原子激发至预设能级,并处理得到能级c原子密度;根据玻尔兹曼分布处理得到能级a原子密度;计算爱因斯坦系数;计算局部热平衡下的基态密度;计算各能级原子密度及总密度。本发明解决了无法通过荧光信号反演氦密度及设备使用成本高的技术问题。
Description
技术领域
本发明涉及通过激光诱导荧光进行密度反演的技术领域,具体涉及基于激光诱导荧光的氦原子密度反演方法。
背景技术
激光诱导荧光技术是指利用激光激发物质发射荧光,通过采集激光脉冲与介质相互作用产生的荧光信号进行判别的一种检测技术。LIF检测系统主要包括激光器、检测光路、光电探测器、信号处理模块几部分。LIF检测作为高灵敏度的检测技术近年来得到了快速发展,其浓度检验精度高,同时具有检测时间短、可实现在线检测等优点,在大气污染监测等领域得到广泛应用。
关于激光诱导荧光测量氦原子密度,目前主要应用于核聚变领域。国外,国际联合的ITER(International Thermonuclear Experimental Reactor,即国际热核聚变实验堆)和德国的AXDEX Upgrade托卡马克装置,通过激光诱导荧光进行了氦原子密度测量或者进行了方案设计,采取的光谱方案主要有两种:入射激光波长388.9nm,探测荧光587.6nm;入射激光波长667.8nm,探测荧光667.8nm。但文献中都未给出具体的通过荧光信号反演氦密度的方法,而国内关于激光诱导荧光测量氦原子密度尚无相关文献。公开号为CN110146410A的现有发明专利文献《基于差分吸收法的原子密度及布居数的测量装置及方法》包括:采用两个脉冲光源,一个作为信号光,一个作为光学泵浦光。两束光合束后沿两条光路传播,其中一路放置原子池,作为信号光路,一路不放置原子池,作为吸收谱线的校准光路。两路光经探测器转化为电信号后相减,能够去除吸收谱线中由于光路及光源不稳定性带来的影响,提高测量精确度。该现有专利文献中利用吸收谱线的每个精细吸收峰分别确定处于各个能级的原子密度,利用完整吸收谱中的每个单独的吸收谱线与能级跃迁的对应关系,进行数据截取和积分操作,可以计算出原子每个能态的布居数,但无法计算得到特定热平衡条件等场景下处于各个能级的原子密度及总密度。公开号为CN101949852A的现有发明专利文献《一种基于光谱标准化的煤质在线检测方法》以经过光谱标准化后的原子特征谱线的强度为纵坐标,以原子特征谱线对应的元素的浓度作为横坐标,建立定标模型。对于未知成分的煤炭进行检测时,经过光谱标准化,根据定标模型即可得到待测元素的浓度。该方法把等离子体温度和电子密度,以及同一元素的离子原子密度比和玻尔兹曼定律考虑进来,从而达到补偿实验参数的波动,提高测量精度的目的。该现有专利文献通过光谱标准化方法进行LIBS定量化分析,以对煤质进行全元素分析,其从激光器出射出的激光经过聚焦透镜2聚焦后作用在定标煤样品3表面,在聚焦点产生等离子体,对于定标煤样中各元素的原子特征谱线,利用元素的离子原子密度比和玻尔兹曼定律,对原子特征谱线的原始强度进行标准化,可知,该现有技术仅考虑了考虑玻尔兹曼定律和离子密度原子密度比的变化,无法解决同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题,且该现有技术用于煤质元素测定,无法用于氦密度的反演。
综上,现有技术存在无法通过荧光信号反演氦密度且设备使用成本高的技术问题。
发明内容
本发明所要解决的技术问题在于如何解决现有技术无法通过荧光信号反演氦密度且设备使用成本高的技术问题。
本发明是采用以下技术方案解决上述技术问题的:基于激光诱导荧光的氦原子密度反演方法包括:
S1、设定氦原子单态能级数据,其中,氦原子单态能级数据包括:能级a、能级b、能级c以及能级d;
S2、设中性氦从能级c跃迁至能级d,获取c至d能级跃迁发射荧光频率以及c至d能级跃迁爱因斯坦系数;
S3、以预置逻辑处理c至d能级跃迁发射荧光频率、c至d能级跃迁爱因斯坦系数及单位时间探测光强,据以得到无激光诱导c能级原子密度;
S4、基于激光诱导,设中性氦吸收预设频率的激光,以从所述能级a跃迁至所述能级b,在向低能级跃迁时,经过能级c向能级d跃迁,将荧光频率保持为c至d级跃迁发射荧光频率,以通过预置逻辑处理得到探测光强差值;
S5、获取探测光强差值与单位时间探测光强的比值,据以获取第一激光诱导原子密度比例关系;
S6、使用预置能量的激光器,将能级a上的原子全部激发至能级b,设能级a的原子密度为na0,能级b到能级c跃迁的爱因斯坦系数为Abc,据以通过下述逻辑处理激光诱导b能级原子密度的增量nb,以得到能级c原子密度:
nc=nbAbc=na0Abc;
S7、设n1为局部热平衡下的基态(1s2)粒子密度,Te为等离子体温度,g为简并度,E为能级能量,根据玻尔兹曼分布,以下述逻辑处理得到能级a原子密度:
S8、利用下述逻辑处理能级c原子密度、能级a原子密度、第一激光诱导原子密度比例关系以及c至d能级跃迁爱因斯坦系数:
从能级表中读取Ea、Ec,通过测量得到I,I0,Te,据以计算得到Abc;
S9、根据无激光诱导能级c原子密度nc0,将I-I0和I0相比,以得到第二激光诱导原子密度比例关系:
以根据下述逻辑处理得到局部热平衡基态密度:
S10、根据玻尔兹曼分布,以下述逻辑处理各能级能量,以得到各能级原子密度:
将各能级原子密度求和得到总密度。
本发明针对同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题,通过激光诱导荧光计算爱因斯坦系数,进而给出局部热平衡条件下四能级密度反演方法,结合氦原子能级数据表,可得到氦原子各能级密度和总密度,且易于推广,可应用到其他原子四能级反演问题中,对于氦原子其他四能级系统或其它元素的四能级系统同样适用。
在更具体的技术方案中,步骤S2中,以下述逻辑处理c至d能级跃迁发射荧光频率以及c至d能级跃迁爱因斯坦系数:
其中,nc0为c能级的密度,V0为探测体积,Ω为探测立体角,k为光学系统透射率。
在更具体的技术方案中,光学系统包括:收光镜筒和光电倍增管,其中,收光镜筒包括501.6nm滤光片,滤光片为窄带宽、高透射率滤光片。
在更具体的技术方案中,步骤S3中,采用下述逻辑处理c至d能级跃迁发射荧光频率、c至d能级跃迁爱因斯坦系数及单位时间探测光强,以得到无激光诱导c能级原子密度:
在更具体的技术方案中,步骤S4中,以下述逻辑处理获取探测光强差值:
其中,nc为由于激光诱导c能级原子密度的增量。
本发明通过激光诱导荧光,解决了同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题。
在更具体的技术方案中,步骤S5中,以下述逻辑处理得到获取第一激光诱导原子密度比例关系:
在更具体的技术方案中,步骤S1到S3中排除激光入射,计算出当前c能级的原子密度,据以荧光诱导计算基态密度。
在更具体的技术方案中,步骤S6中采用667.8nm的半导体激光器,以将能级a上的原子全部激发至能级b。
本发明的入射激光可通过半导体激光器产生,相对染料激光器造价低。
在更具体的技术方案中,步骤S5到S8中,通过将探测光强差值I-I0与第一单位时间探测光强I0相比,以消去未知的局部热平衡基态密度n1。
本发明通过把I-I0和I0相比,将暂时未知的n1消去,根据能级能量和可测量量I,I0,Te,计算出Abc,解决反演中的主要未知量问题,此方法也可应用到其他原子四能级反演问题中。
在更具体的技术方案中,基于激光诱导荧光的氦原子密度反演系统包括:
能级数据模块,用以设定氦原子单态能级数据,其中,氦原子单态能级数据包括:能级a、能级b、能级c以及能级d;
被动光强获取模块,设中性氦从能级c跃迁至能级d,获取c至d能级跃迁发射荧光频率以及c至d能级跃迁爱因斯坦系数,探测光强获取模块与能级数据模块连接;
无激光诱导原子密度获取模块,用于以预置逻辑处理c至d级跃迁发射荧光频率、c至d级跃迁爱因斯坦系数及单位时间探测光强,据以得到无激光诱导能级c原子密度,能级c原子密度获取模块与被动光强获取模块连接;
探测光强差值模块,用以基于激光诱导,设中性氦吸收预设频率的激光,以从能级a跃迁至能级b,在向低能级跃迁时,经过能级c向能级d跃迁,将荧光频率保持为c至d能级跃迁发射荧光频率,以通过预置逻辑处理得到探测光强差值,探测光强差值模块与能级数据模块连接;
密度比例模块,用以获取探测光强差值与单位时间探测光强的比值,据以获取第一激光诱导原子密度比例关系,密度比例模块与探测光强差值模块及探测光强获取模块连接;
能级c原子密度模块,用以使用预置能量的激光器,将能级a上的原子全部激发至能级b,设能级a的原子密度为na0,能级b到能级c跃迁的爱因斯坦系数为Abc,据以通过下述逻辑处理激光诱导b能级原子密度的增量nb,以得到能级c原子密度,能级c原子密度模块与能级数据模块连接:
nc=nbAbc=na0Abc;
能级a原子密度模块,用以设n1为局部热平衡下的基态(1s2)粒子密度,Te为等离子体温度,g为简并度,E为能级能量,根据玻尔兹曼分布,以下述逻辑处理得到能级a原子密度,能级a原子密度模块与能级数据模块连接:
爱因斯坦系数获取模块,用以利用下述逻辑处理c能级原子密度、能级a原子密度、第一激光诱导原子密度比例关系以及c至d能级跃迁爱因斯坦系数:
从能级表中读取Ea、Ec,通过测量得到I,I0,Te,据以计算得到Abc,爱因斯坦系数获取模块与能级c原子密度模块、能级a原子密度模块及密度比例模块连接;
基态密度模块,用以根据无激光诱导能级c原子密度nc0,将I-I0和I0相比,以得到第二激光诱导原子密度比例关系:
以根据下述逻辑处理得到局部热平衡基态密度,基态密度模块与能级c原子密度模块及探测光强差值模块连接:
原子密度结果获取模块,用以根据玻尔兹曼分布,以下述逻辑处理各能级能量,以得到各能级原子密度:
将各能级原子密度求和得到总密度,原子密度结果获取模块与基态密度模块连接。
本发明相比现有技术具有以下优点:本发明针对同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题,通过激光诱导荧光计算爱因斯坦系数,进而给出局部热平衡条件下四能级密度反演方法,结合氦原子能级数据表,可得到氦原子各能级密度和总密度,且易于推广,可应用到其他原子四能级反演问题中,对于氦原子其他四能级系统或其它元素的四能级系统同样适用。
本发明通过激光诱导荧光,解决了同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题。
本发明的入射激光可通过半导体激光器产生,相对染料激光器造价低。
本发明通过把I-I0和I0相比,将暂时未知的n1消去,根据能级能量和可测量量I,I0,Te,计算出Abc,解决反演中的主要未知量问题,此方法也可应用到其他原子四能级反演问题中。本发明解决了现有技术中无法通过荧光信号反演氦密度及设备使用成本高的技术问题。
附图说明
图1为本发明实施例1的基于激光诱导荧光的氦原子密度反演方法利用的氦原子能级示意图;
图2为本发明实施例1的基于激光诱导荧光的氦原子密度反演方法流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1所示,单态能级1s2p用符号a表示,能级1s3d用符号b表示,能级1s3p用符号c表示,能级1s2s用符号d表示。通过667.8nm激光将氦原子从1s2p能级激发到1s3d能级,探测从能级1s3p到能级1s2s辐射的501.6nm荧光。
如图2所示,步骤S101:设定氦原子单态能级,设中性氦从能级c跃迁至能级d,获取c至d级跃迁发射荧光频率以及c至d级跃迁爱因斯坦系数,处理单位时间探测光强;在本实施例中,氦原子反演中用到的单态各能级1s2p、1s3d、1s3p和1s2s分别用符号a、b、c和d表示,设中性氦从c能级向d能级跃迁时,发射的荧光频率νcd,c能级到d能级跃迁的爱因斯坦系数为Acd,无激光入射时,设单位时间探测到的荧光光强为I0,则:
其中,nc0为c能级的密度。探测体积V0、探测立体角Ω、光学系统透射率k根据实验装置的实际条件代入。光学系统包括收光镜筒和光电倍增管,收光镜筒中应使用窄带宽、高透射率的501.6nm滤光片;
步骤S102:求取c能级的原子密度;在本实施例中,c能级的原子密度可表示为:
在本实施例中,所述步骤S101到S102中无激光入射,先计算出此时的c能级的原子密度,有荧光诱导计算基态密度时需应用。
步骤S103:考虑激光诱导,设中性氦吸收预设频率的激光,以从能级a跃迁至能级b,在向低能级跃迁时,经过能级c向能级d跃迁,将荧光频率保持为c至d级跃迁发射荧光频率,处理得到单位时间探测到的荧光光强;在本实施例中,考虑激光诱导,设中性氦吸收频率为ν的激光后,从a能级跃迁到b能级,向低能级跃迁时,经过c能级向d能级跃迁,荧光频率仍是νcd。设单位时间探测到的荧光光强为I,则:
nc为由于激光诱导c能级原子密度的增量,若系统中使用锁相放大器,则探测光强无需再减去I0;
步骤S104:获取探测光强差值与单位时间探测光强的比值,在本实施例中,把I-I0和I0相比得:
步骤S105:利用激光器对氦原子进行激光诱导,将氦原子激发至预设能级b,并处理得到能级c原子密度;在本实施例中,使用能量足够的667.8nm激光器,a能级上的原子全部激发至b能级。设a能级的原子密度为na0,b能级到c能级跃迁的爱因斯坦系数为Abc,由于激光诱导b能级原子密度的增量为nb,则
nc=nbAbc=na0Abc
激光器应配合波长计,以便在激光器波长不稳时及时调整;在本实施例中,步骤S105中要求,必须使用能量足够的激光器,将a能级上的原子全部激发至b能级。
步骤S106:根据玻尔兹曼分布处理得到第一能级原子密度;在本实施例中,设n1为局部热平衡下的基态(1s2)粒子密度,Te为等离子体温度,g为简并度(由于使用的是单态光谱,g为1),E为能级能量,根据玻尔兹曼分布:
步骤S107:以下述逻辑计算爱因斯坦系数:
从能级表中读取Ea、Ec,通过测量得到I,I0,Te,可计算出Abc,解决反演中的主要未知量问题,此方法也可应用到其他原子四能级反演问题中;
步骤S108:计算局部热平衡下的基态密度,在本实施例中,再利用无激光诱导时计算出的nc0,重新把I-I0和I0相比:
可得局部热平衡下的基态密度。
步骤S109:计算各能级原子密度及总密度,在本实施例中,根据玻尔兹曼分布和各能级能量,可得各能级原子密度为
最后将各能级原子密度求和得到总密度。
其中,单态各能级能量为
通过波长为587.6nm的激光将氦原子从三态1s2p能级激发到1s3d能级,探测从能级1s3p到能级1s2s辐射的388.9nm荧光时,所需三态各能级能量为
综上,本发明针对同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题,通过激光诱导荧光计算爱因斯坦系数,进而给出局部热平衡条件下四能级密度反演方法,结合氦原子能级数据表,可得到氦原子各能级密度和总密度,且易于推广,可应用到其他原子四能级反演问题中,对于氦原子其他四能级系统或其它元素的四能级系统同样适用。
本发明通过激光诱导荧光,解决了同一主量子数、相邻角量子数的能级间爱因斯坦系数未知问题。
本发明的入射激光可通过半导体激光器产生,相对染料激光器造价低。
本发明通过把I-I0和I0相比,将暂时未知的n1消去,根据能级能量和可测量量I,I0,Te,计算出Abc,解决反演中的主要未知量问题,此方法也可应用到其他原子四能级反演问题中。本发明解决了现有技术中无法通过荧光信号反演氦密度及设备使用成本高的技术问题。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.基于激光诱导荧光的氦原子密度反演方法,其特征在于,所述方法包括:
S1、设定氦原子单态能级数据,其中,所述氦原子单态能级数据包括:能级a、能级b、能级c以及能级d;
S2、设中性氦从所述能级c跃迁至所述能级d,获取c至d能级跃迁发射荧光频率以及c至d能级跃迁爱因斯坦系数;
S3、以预置逻辑处理所述c至d能级跃迁发射荧光频率、所述c至d能级跃迁爱因斯坦系数及所述单位时间探测光强,据以得到无激光诱导c能级原子密度;
S4、基于激光诱导,设所述中性氦吸收预设频率的激光,以从所述能级a跃迁至所述能级b,在向低能级跃迁时,经过所述能级c向所述能级d跃迁,将荧光频率保持为所述c至d能级跃迁发射荧光频率,以通过预置逻辑处理得到探测光强差值;
S5、获取所述探测光强差值与所述单位时间探测光强的比值,据以获取第一激光诱导原子密度比例关系;
S6、使用预置能量的激光器,将所述能级a上的原子全部激发至所述能级b,设所述能级a的原子密度为na0,所述能级b到所述能级c跃迁的爱因斯坦系数为Abc,据以通过下述逻辑处理激光诱导b能级原子密度的增量nb,以得到能级c原子密度:
nc=nbAbc=na0Abc;
S7、设n1为局部热平衡下的基态(1s2)粒子密度,Te为等离子体温度,g为简并度,E为能级能量,根据玻尔兹曼分布,以下述逻辑处理得到能级a原子密度:
S8、利用下述逻辑处理所述能级c原子密度、所述能级a原子密度、所述第一激光诱导原子密度比例关系以及所述c至d能级跃迁爱因斯坦系数:
从能级表中读取Ea、Ec,通过测量得到I,I0,Te,据以计算得到Abc;
S9、根据所述无激光诱导能级c原子密度nc0,将I-I0和I0相比,以得到第二激光诱导原子密度比例关系:
以根据下述逻辑处理得到局部热平衡基态密度:
S10、根据玻尔兹曼分布,以下述逻辑处理各能级能量,以得到各能级原子密度:
将各能级原子密度求和得到总密度。
3.根据权利要求2所述的基于激光诱导荧光的氦原子密度反演方法,其特征在于,所述光学系统包括:收光镜筒和光电倍增管,其中,所述收光镜筒包括501.6nm滤光片,所述滤光片为窄带宽、高透射率滤光片。
7.根据权利要求1所述的基于激光诱导荧光的氦原子密度反演方法,其特征在于:所述步骤S1到S3中排除激光入射,计算出当前c能级的原子密度,据以荧光诱导计算基态密度。
8.根据权利要求1所述的基于激光诱导荧光的氦原子密度反演方法,其特征在于:所述步骤S6中采用667.8nm的半导体激光器,以将所述能级a上的原子全部激发至所述能级b。
9.根据权利要求1所述的基于激光诱导荧光的氦原子密度反演方法,其特征在于:所述步骤S5到S8中,通过将所述探测光强差值I-I0与所述所述第一单位时间探测光强I0相比,以消去未知的所述局部热平衡基态密度n1。
10.基于激光诱导荧光的氦原子密度反演系统,其特征在于,所述系统包括:
能级数据模块,用以设定氦原子单态能级数据,其中,所述氦原子单态能级数据包括:能级a、能级b、能级c以及能级d;
探测光强获取模块,设中性氦从所述能级c跃迁至所述能级d,获取c至d能级跃迁发射荧光频率以及c至d能级跃迁爱因斯坦系数,所述探测光强获取模块与所述能级数据模块连接;
无激光诱导原子密度获取模块,用于以预置逻辑处理所述c至d能级跃迁发射荧光频率、所述c至d能级跃迁爱因斯坦系数及所述单位时间探测光强,据以得到无激光诱导能级c原子密度,所述能级c原子密度获取模块与所述探测光强获取模块连接;
探测光强差值模块,用以基于激光诱导,设所述中性氦吸收预设频率的激光,以从所述能级a跃迁至所述能级b,在向低能级跃迁时,经过所述能级c向所述能级d跃迁,将荧光频率保持为所述c至d能级跃迁发射荧光频率,以通过预置逻辑处理得到探测光强差值,所述探测光强差值模块与所述能级数据模块连接;
密度比例模块,用以获取所述探测光强差值与所述第一单位时间探测光强的比值,据以获取第一激光诱导原子密度比例关系,所述密度比例模块与所述探测光强差值模块及所述探测光强获取模块连接;
能级c原子密度模块,用以使用预置能量的激光器,将所述能级a上的原子全部激发至所述能级b,设所述能级a的原子密度为na0,所述能级b到所述能级c跃迁的爱因斯坦系数为Abc,据以通过下述逻辑处理激光诱导b能级原子密度的增量nb,以得到能级c原子密度,所述能级c原子密度模块与所述能级数据模块连接:
nc=nbAbc=na0Abc;
能级a原子密度模块,用以设n1为局部热平衡下的基态(1s2)粒子密度,Te为等离子体温度,g为简并度,E为能级能量,根据玻尔兹曼分布,以下述逻辑处理得到能级a原子密度,所述能级a原子密度模块与所述能级数据模块连接:
爱因斯坦系数获取模块,用以利用下述逻辑处理所述第三能级原子密度、所述第一能级原子密度、所述第一激光诱导原子密度比例关系以及所述c至d能级跃迁爱因斯坦系数:
从能级表中读取Ea、Ec,通过测量得到I,I0,Te,据以计算得到Abc,所述爱因斯坦系数获取模块与所述能级c原子密度模块、所述能级a原子密度模块及所述密度比例模块连接;
基态密度模块,用以根据所述无激光诱导能级c原子密度nc0,将I-I0和I0相比,以得到第二激光诱导原子密度比例关系:
以根据下述逻辑处理得到局部热平衡基态密度,所述基态密度模块与所述能级c原子密度模块及所述探测光强差值模块连接:
原子密度结果获取模块,用以根据玻尔兹曼分布,以下述逻辑处理各能级能量,以得到各能级原子密度:
将各能级原子密度求和得到总密度,所述原子密度结果获取模块与所述基态密度模块连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211138153.3A CN115494039B (zh) | 2022-09-19 | 2022-09-19 | 基于激光诱导荧光的氦原子密度反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211138153.3A CN115494039B (zh) | 2022-09-19 | 2022-09-19 | 基于激光诱导荧光的氦原子密度反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115494039A true CN115494039A (zh) | 2022-12-20 |
CN115494039B CN115494039B (zh) | 2024-05-14 |
Family
ID=84470879
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211138153.3A Active CN115494039B (zh) | 2022-09-19 | 2022-09-19 | 基于激光诱导荧光的氦原子密度反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115494039B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5657340A (en) * | 1996-04-19 | 1997-08-12 | The Aerospace Corporation | Rubidium atomic clock with fluorescence optical pumping and method using same |
CN101949852A (zh) * | 2010-07-30 | 2011-01-19 | 清华大学 | 一种基于光谱标准化的煤质在线检测方法 |
CN105092542A (zh) * | 2014-05-07 | 2015-11-25 | 中国科学院大连化学物理研究所 | 脉冲纳米氦液滴光谱仪 |
CN105571579A (zh) * | 2014-10-24 | 2016-05-11 | 霍尼韦尔国际公司 | 具有多普勒致宽增益介质的芯片上环形激光陀螺仪 |
CN110837109A (zh) * | 2019-10-22 | 2020-02-25 | 山西大学 | 原子激发态光谱获得方法和超精细能级测量方法与装置 |
-
2022
- 2022-09-19 CN CN202211138153.3A patent/CN115494039B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5657340A (en) * | 1996-04-19 | 1997-08-12 | The Aerospace Corporation | Rubidium atomic clock with fluorescence optical pumping and method using same |
CN101949852A (zh) * | 2010-07-30 | 2011-01-19 | 清华大学 | 一种基于光谱标准化的煤质在线检测方法 |
CN105092542A (zh) * | 2014-05-07 | 2015-11-25 | 中国科学院大连化学物理研究所 | 脉冲纳米氦液滴光谱仪 |
CN105571579A (zh) * | 2014-10-24 | 2016-05-11 | 霍尼韦尔国际公司 | 具有多普勒致宽增益介质的芯片上环形激光陀螺仪 |
CN110837109A (zh) * | 2019-10-22 | 2020-02-25 | 山西大学 | 原子激发态光谱获得方法和超精细能级测量方法与装置 |
Non-Patent Citations (1)
Title |
---|
宋根宗;耿平;张多;杨玉玲;: "He―Ne放电管荧光光谱增益的实验测试与研究", 大学物理实验, no. 01, 25 March 2007 (2007-03-25) * |
Also Published As
Publication number | Publication date |
---|---|
CN115494039B (zh) | 2024-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Carlstrom et al. | Design and operation of the multipulse Thomson scattering diagnostic on DIII‐D | |
Dagenais et al. | Investigation of two-time correlations in photon emissions from a single atom | |
Alfier et al. | New Thomson scattering diagnostic on RFX-mod | |
Podesta et al. | The NSTX fast-ion D-alpha diagnostic | |
Rafac et al. | Measurement of the ratio of the cesium D-line transition strengths | |
Seraydarian et al. | Multichordal charge exchange recombination spectroscopy on the Doublet III tokamak | |
CN106018878A (zh) | 等离子体速度测量方法及系统 | |
CN105043930A (zh) | 微结构碱金属气室碱金属蒸汽原子密度的检测装置和方法 | |
Kurskiev et al. | Digital filter polychromator for Thomson scattering applications | |
Rivero-Rodriguez et al. | Upgrade and absolute calibration of the JET scintillator-based fast-ion loss detector | |
Milder et al. | Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering | |
CN115494039B (zh) | 基于激光诱导荧光的氦原子密度反演方法 | |
Morgan et al. | Post-irradiation examination of optical components for advanced fission reactor instrumentation | |
JP2011128078A (ja) | ガス中の煤塵濃度計測方法及び燃焼設備の運転方法 | |
JPS59176649A (ja) | 粒子分析装置 | |
Cima et al. | Correlation radiometry of electron cyclotron radiation in TEXT‐U | |
Martin | Reciprocity between emission and absorption for rare earth ions in glass | |
Tsang et al. | Optical beam profile monitor and residual gas fluorescence at the relativistic heavy ion collider polarized hydrogen jet | |
CN110426653B (zh) | 一种测量光抽运率的方法 | |
Dimock et al. | A compact Thomson scattering system | |
Gong et al. | Development of vertical edge Thomson scattering diagnostic on HL-2M Tokamak | |
Wiesinger et al. | Trap-integrated fluorescence detection with silicon photomultipliers for sympathetic laser cooling in a cryogenic Penning trap | |
McLean et al. | Quantification of radiating species in the DIII-D divertor in the transition to detachment using extreme ultraviolet spectroscopy | |
Tournianski et al. | Use of one dimensional D α camera to measure edge electron density gradients | |
Allen et al. | Measurement of sodium groundstate density profile in a metal-halide lamp using laser absorption spectroscopy |
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 |