CN103472033B - 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 - Google Patents
一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 Download PDFInfo
- Publication number
- CN103472033B CN103472033B CN201310467760.9A CN201310467760A CN103472033B CN 103472033 B CN103472033 B CN 103472033B CN 201310467760 A CN201310467760 A CN 201310467760A CN 103472033 B CN103472033 B CN 103472033B
- Authority
- CN
- China
- Prior art keywords
- theta
- lambda
- particle
- refractive index
- laser
- 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.)
- Active
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,属于颗粒光学特性测量技术领域,解决了现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。本发明通过建立球形颗粒系反射信号和透射信号测量的正问题和逆问题求解模型,采用连续激光,测量得到颗粒系的半球反射信号和半球透射信号,采用现有Mie理论模型,能很精确的反应出颗粒的电磁散射特性。本发明适用于测量球形颗粒光谱复折射率。
Description
技术领域
本发明属于颗粒光学特性测量技术领域。
背景技术
在自然界及工业生产中,很多物质及产品呈现颗粒状态,据统计,工业中有50%以上的产品与中间产物呈颗粒状。颗粒的吸收散射以及发射等性质在许多工程及环境系统内扮演着重要的角色,因此对颗粒的光谱复折射率的求解就显得尤为重要。颗粒的光谱复折射率属于颗粒的固有物性参数,与颗粒的组成成分、温度水平及表面的状况等因素有关。然而,颗粒的光谱复折射率一般不能通过实验直接测量得到,通常都采用某些实验测得的参数结合相关的反演理论模型进行计算。
常见的颗粒的光谱复折射率的实验研究方法有反射法、透射法和散射法。其中,反射法是通过测量粒子压片的反射率,结合相应的理论模型,求颗粒的光谱复折射率的反问题研究方法。该方法得到主要是块状的物质的光谱复折射率,由于实验所制压片的表面并不光滑,不服从镜反射原理,因此使得实验测得结果有很大偏差。透射法是通过测量颗粒压片的透射率,并结合K-K关系式和相应的理论模型,反演求解颗粒的光谱复折射率的研究方法。该方法可以使得颗粒保持自然状态,并且实验所需设备简单、精度高、操作方便和适应范围广等优点,但是K-K关系式需要考虑全波长范围的积分,而实验中的激光器通常都是有限波长范围,这就增加了实验测量和理论计算的难度。散射法是通过测量激光照射到颗粒系上后,不同方向的散射光的强度,然后结合相应的理论模型求解颗粒的光谱复折射率的反问题研究方法。该方法能很好的保持颗粒的自然状态,但是确定某一波长下的光谱复折射率至少需要测三个不相关方向的散射量,这就增加实验测量过程的难度。
综上所述现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。
发明内容
本发明为了解决现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。提出了一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法。
本发明所述一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用连续激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,0<θc<π/2;
步骤三、采用光电接收器在样本颗粒容器的左侧表面测量半球反射信号,在样本颗粒容器的右侧表面测量半球透射信号,并记录样本颗粒的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数P(D);
步骤五、利用待测颗粒光谱复折射率方程m(λ)=n(λ)+ik(λ),通过待测颗粒的折射率范围对该待测物质的折射率n进行取值,通过待测颗粒的吸收因子的取值范围对该待测颗粒的吸收因子k进行取值,获得该种待测颗粒的光谱复折射率m(λ);其中i为虚数单位;
步骤六、通过现有Mie理论,利用步骤五获得的该种待测颗粒的光谱复折射率m(λ),获得单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ);
步骤七、利用步骤六获得的单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ),获得颗粒系的吸收系数κa(λ)和散射系数κs(λ);
步骤八、利用步骤七获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),通过对辐射传输方程,获得计算域内的辐射强度场Iλ(L,θ);
步骤九、通过公式:
获得左侧边界的反射信号的估计值Rpred,λ和右侧边界的透射信号的估计值τpred,λ;
式中,I0,λ是波长为λ的连续激光的辐射强度;Iλ(0,θ)是波长为λ的激光在θ方向上x=0处的散射光的辐射强度,θ为辐射方向角;
Iλ(L,θ)是波长为λ的激光在θ方向上x=L处的散射光的辐射强度,θ为辐射方向角;
Ic,λ(L,θc)为波长为λ的连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,L为样本厚度;
步骤十、利用步骤三中光电接收器测量获得的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ与步骤九中获得左侧边界的半球反射信号的估计值Rpred,λ和右侧边界的半球透射信号的估计值τpred,λ,获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值,和对待测颗粒的吸收因子k取值,作为结果,完成利用连续激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
本发明通过建立球形颗粒系反射信号和透射信号测量的正问题和逆问题求解模型,采用连续激光,测量得到颗粒系的半球反射信号和半球透射信号,采用Mie理论模型,精确的反应出颗粒的电磁散射特性;采用微粒群优化算法,使测量精确度提高了10%,同时简化了测量过程。
附图说明
图1为连续激光辐照下球形颗粒系反射信号和透射信号测量示意图;
图2为本发明所述一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法流程图。
具体实施方式
具体实施方式一、结合图1、图2说明本实施方式,本实施方式所述一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用连续激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,0<θc<π/2;
步骤三、采用光电接收器在样本颗粒容器的左侧表面测量半球反射信号,在样本颗粒容器的右侧表面测量半球透射信号,并记录样本颗粒的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数P(D);
步骤五、利用待测颗粒光谱复折射率方程m(λ)=n(λ)+ik(λ),通过待测颗粒的折射率范围对该待测物质的折射率n进行取值,通过待测颗粒的吸收因子的取值范围对该待测颗粒的吸收因子k进行取值,获得该种待测颗粒的光谱复折射率m(λ);其中i为虚数单位;
步骤六、通过现有Mie理论,利用步骤五获得的该种待测颗粒的光谱复折射率m(λ),获得单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ);
步骤七、利用步骤六获得的单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ),获得颗粒系的吸收系数κa(λ)和散射系数κs(λ);
步骤八、利用步骤七获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),通过对辐射传输方程,获得计算域内的辐射强度场Iλ(L,θ);
步骤九、通过公式:
获得左侧边界的反射信号的估计值Rpred,λ和右侧边界的透射信号的估计值τpred,λ;
式中,I0,λ是波长为λ的连续激光的辐射强度;Iλ(0,θ)是波长为λ的激光在θ方向上x=0处的散射光的辐射强度,θ为辐射方向角;
Iλ(L,θ)是波长为λ的激光在θ方向上x=L处的散射光的辐射强度,θ为辐射方向角;
Ic,λ(L,θc)为波长为λ的连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,L为样本厚度;
步骤十、利用步骤三中光电接收器测量获得的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ,与步骤九中获得左侧边界的半球反射信号的估计值Rpred,λ和右侧边界的半球透射信号的估计值τpred,λ,获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值,和对待测颗粒的吸收因子k取值,作为结果,完成利用连续激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
本实施方式首先设计悬浮颗粒系的稳态辐射传输物理模型,然后建立相应的数学模型和求解方法,通过测量得到颗粒系的半球反射信号和半球透射信号,利用逆问题理论模型的重建出球形颗粒光谱复折射率。
具体实施方式二、本实施方式是对具体实施方式一所述的一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤七中获得颗粒系的吸收系数κa(λ)和散射系数κs(λ)为,通过公式:
计算获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),式中,D为所测颗粒的粒径Dmax和Dmin分别为所测颗粒的最大粒径和最小粒径,N0所测颗粒的总个数;
具体实施方式三、本实施方式是对具体实施方式一所述的一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤八获得计算域内的辐射强度场I(x,θ,t,λ)的方法为:利用瞬态辐射传输方程:
实现,式中Iλ(x,θ)是波长为λ的激光在θ方向x处的辐射强度,x为待求辐射场中位置,θ待求辐射方向,κa(λ)为粒子系的吸收系数,κs(λ)为粒子系的散射系数;Iλ(x,θ′)是波长为λ的激光在θ′方向x处的辐射强度;θ′为入射方向,Φ(θ′,θ)是从θ′方向入射并从θ方向散射出去的散射相函数,Iλ(0,θ)是波长为λ的激光在θ方向x=0处的辐射强度,Iλ(L,θ)是波长为λ的激光在θ方向x=L处的辐射强度;
c是介质中的光速;
ρ0是由环境进入粒子系时的反射率;
ρ1为由粒子系进入环境时的反射率;
λ是激光在球形颗粒系中的波长;
I0,λ(θc)是波长为λ沿着θc角度入射的连续激光的辐射强度。
具体实施方式四、本实施方式是对具体实施方式一所述的一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤十中获得目标函数Fobj通过公式:
实现。
Claims (3)
1.一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用连续激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,0<θc<π/2;
步骤三、采用光电接收器在样本颗粒容器的左侧表面测量半球反射信号,在样本颗粒容器的右侧表面测量半球透射信号,并记录样本颗粒的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数P(D);
步骤五、利用待测颗粒光谱复折射率方程m(λ)=n(λ)+ik(λ),通过待测颗粒的折射率范围对该待测物质的折射率n进行取值,通过待测颗粒的吸收因子的取值范围对该待测颗粒的吸收因子k进行取值,获得该待测颗粒的光谱复折射率m(λ);其中i为虚数单位;
步骤六、通过现有Mie理论,利用步骤五获得的该种待测颗粒的光谱复折射率m(λ),获得单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ);
步骤七、利用步骤六获得的单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ),获得颗粒系的吸收系数κa(λ)和散射系数κs(λ);
步骤八、利用步骤七获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),通过瞬态辐射传输方程,获得计算域内的辐射强度场Iλ(L,θ);
所述获得计算域内的辐射强度场的方法为利用瞬态辐射传输方程:
实现,式中Iλ(x,θ)是波长为λ的激光在θ方向x处的辐射强度,x为待求辐射场中位置,θ待求辐射方向,κa(λ)为粒子系的吸收系数,κs(λ)为粒子系的散射系数;Iλ(x,θ′)是波长为λ的激光在θ′方向x处的辐射强度;θ′为入射方向,Φ(θ′,θ)是从θ′方向入射并从θ方向散射出去的散射相函数,Iλ(0,θ)是波长为λ的激光在θ方向x=0处的辐射强度,Iλ(L,θ)是波长为λ的激光在θ方向x=L处的辐射强度;
c是介质中的光速;
ρ0是由环境进入粒子系时的反射率;
ρ1为由粒子系进入环境时的反射率;
I0,λ(θc)是波长为λ沿着θc角度入射的连续激光的辐射强度;
步骤九、通过公式:
获得左侧边界的反射信号的估计值Rpred,λ和右侧边界的透射信号的估计值τpred,λ;
式中,I0,λ是波长为λ的连续激光的辐射强度;Iλ(0,θ)是波长为λ的激光在θ方向上x=0处的散射光的辐射强度,θ为辐射方向角;
Iλ(L,θ)是波长为λ的激光在θ方向上x=L处的散射光的辐射强度,θ为辐射方向角;
Ic,λ(L,θc)为波长为λ的连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,L为样本厚度;
步骤十、利用步骤三中采用光电接收器测量获得的左侧边界处的半球反射信号Rexp,λ和右侧边界处的半球透射信号τexp,λ与步骤九中获得左侧边界的半球反射信号的估计值Rpred,λ和右侧边界的半球透射信号的估计值τpred,λ,获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值,和对待测颗粒的吸收因子k取值,作为结果,完成利用连续激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
2.根据权利要求1所述的一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,步骤七中获得颗粒系的吸收系数κa(λ)和散射系数κs(λ)为,通过公式:
计算获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),式中,D为所测颗粒的粒径Dmax和Dmin分别为所测颗粒的最大粒径和最小粒径,N0所测颗粒的总个数。
3.根据权利要求1所述的一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,步骤十中获得目标函数Fobj通过公式:
实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310467760.9A CN103472033B (zh) | 2013-10-09 | 2013-10-09 | 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310467760.9A CN103472033B (zh) | 2013-10-09 | 2013-10-09 | 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103472033A CN103472033A (zh) | 2013-12-25 |
CN103472033B true CN103472033B (zh) | 2015-12-02 |
Family
ID=49796973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310467760.9A Active CN103472033B (zh) | 2013-10-09 | 2013-10-09 | 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103472033B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104634705B (zh) * | 2015-03-11 | 2017-06-23 | 哈尔滨工业大学 | 一种基于连续激光的球形颗粒光谱复折射率与颗粒系粒径分布的获得方法 |
CN105203437B (zh) * | 2015-09-21 | 2017-11-17 | 哈尔滨工业大学 | 基于前向散射多角度测量的球形颗粒光学常数与颗粒系粒径分布的同时重构方法 |
CN106644852B (zh) * | 2016-10-17 | 2019-03-29 | 哈尔滨工业大学 | 基于超短脉冲激光辐照同时获取球形颗粒光学常数与粒径分布的测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003936A (zh) * | 2010-09-14 | 2011-04-06 | 浙江大学 | 同时测量液滴位置、粒径和复折射率的方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7106442B2 (en) * | 2003-04-29 | 2006-09-12 | Silcott David B | Multi-spectral optical method and system for detecting and classifying biological and non-biological particles |
US20120116739A1 (en) * | 2010-11-08 | 2012-05-10 | The Board Of Trustees Of The University Of Illinois | Infrared Microspectroscopy for Intact Fibers |
-
2013
- 2013-10-09 CN CN201310467760.9A patent/CN103472033B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102003936A (zh) * | 2010-09-14 | 2011-04-06 | 浙江大学 | 同时测量液滴位置、粒径和复折射率的方法和装置 |
Non-Patent Citations (3)
Title |
---|
Inverse Radiation Problem for Determination of Optical Constants of Fly-Ash Particles;L. M. Ruan et al.;《International Journal of Thermophysics》;20070831;第28卷(第4期);1322-1341 * |
哈尔滨地区气溶胶粒子光学常数的实验研究;林莉等;《节能技术》;20130731;第31卷(第4期);340-344 * |
悬浊液光谱透射法测量烟尘粒子的复折射率;邢键 等;《光谱学与光谱分析》;20101231;第30卷(第12期);3371-3374 * |
Also Published As
Publication number | Publication date |
---|---|
CN103472033A (zh) | 2013-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487356B (zh) | 一种基于短脉冲激光透反射信号的球形颗粒光谱复折射率测量方法 | |
CN103499521B (zh) | 纳米颗粒关键几何特征量的测量方法 | |
CN104634705B (zh) | 一种基于连续激光的球形颗粒光谱复折射率与颗粒系粒径分布的获得方法 | |
CN109061668B (zh) | 一种多视场偏振激光雷达探测系统及探测冰云的方法 | |
WO2021088376A1 (zh) | 一种利用散射光的偏振差异测量颗粒折射率的方法及系统 | |
CN104515722A (zh) | 一种颗粒粒度及浓度光散射测量方法 | |
CN105928850B (zh) | 一种光散射法颗粒粒径分布反演的索套回归方法 | |
CN103398964A (zh) | 一种基于腔增强技术的气体探测方法 | |
CN105988124A (zh) | 基于吸收消光比标定廓线的气溶胶污染物分布探测方法 | |
CN101980000A (zh) | 浑浊介质微粒运动特性的完整、高分辨测试方法 | |
CN103472033B (zh) | 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 | |
CN103592103A (zh) | 基于激光消光法的小通道液固两相流参数测量装置及方法 | |
CN109883902A (zh) | 基于日盲紫外光圆偏振的雾霾粒子检测装置及其检测方法 | |
CN103499393A (zh) | 光谱的测量方法 | |
CN103471968B (zh) | 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 | |
Zhang et al. | Simultaneous identification of optical constants and PSD of spherical particles by multi-wavelength scattering–transmittance measurement | |
Wang et al. | A wide dynamic range and high resolution all-fiber-optic turbidity measurement system based on single photon detection technique | |
CN105203437B (zh) | 基于前向散射多角度测量的球形颗粒光学常数与颗粒系粒径分布的同时重构方法 | |
CN105092444A (zh) | 纳米颗粒浓度及几何特征量联合分布的测量方法 | |
CN203163700U (zh) | 一种测量部分相干高斯光束波前相位半径的装置 | |
Bao et al. | Study of a laser echo in an inhomogeneous dust environment with a continuous field Monte Carlo radiative transfer model | |
GB2460305A (en) | Determining the particle size distribution of a suspension from the effective index of refraction of the suspension | |
Patil et al. | Refractometric fiber optic sensor for detecting salinity of water | |
CN115266509B (zh) | 一种基于激光雷达的水下垂向悬浮物浓度探测方法及系统 | |
CN107271336B (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 |