CN103471968A - 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 - Google Patents
一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 Download PDFInfo
- Publication number
- CN103471968A CN103471968A CN2013104670436A CN201310467043A CN103471968A CN 103471968 A CN103471968 A CN 103471968A CN 2013104670436 A CN2013104670436 A CN 2013104670436A CN 201310467043 A CN201310467043 A CN 201310467043A CN 103471968 A CN103471968 A CN 103471968A
- Authority
- CN
- China
- Prior art keywords
- particle
- lambda
- omega
- theta
- frequency domain
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,属于颗粒光学特性测量技术领域。解决了现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。通过建立球形颗粒系频域反射信号和频域透射信号测量的正问题和逆问题求解模型,采用粒度分析仪测量得到颗粒系的粒径分布情况,最后基于这些信号结合逆问题求解技术获得颗粒的光谱复折射率在已知颗粒其他物性参数的前提下,结合频域辐射传输模型,运用Mie理论结合微粒群优化算法反演获得球形颗粒光谱复折射率的方法。本发明适用于测量颗粒的光谱复折射率。
Description
技术领域
本发明属于颗粒光学特性测量技术领域。
背景技术
在自然界及工业生产中,很多物质及产品呈现颗粒状态,例如黏土、灰尘和沙子;人类的食物也是往往是颗粒状,比如大米,豆子、蔗糖和盐;还有很多加工物,例如煤炭、水泥和化肥等等。颗粒的吸收散射以及发射等性质在许多工程及环境系统内扮演着重要的角色,因此对颗粒的光谱复折射率的求解就显得尤为重要。颗粒的光谱复折射率属于颗粒的基本物性参数,与颗粒的组成成分、温度水平及表面的状况等因素有关。然而,颗粒的光谱复折射率不能通过实验直接测量得到,通常都采用某些实验测得的参数结合相关的反演理论模型进行计算。
常见的颗粒光学特性的实验研究方法有反射法、透射法和散射法。其中,反射法是通过测量颗粒压片的反射率,结合相应的理论模型,求颗粒光谱复折射率的反问题研究方法。该方法得到主要是块状的物质的光谱复折射率,由于实验所制压片的表面并不光滑,不服从镜反射原理,因此使得实验测得结果有很大偏差。透射法是通过测量颗粒压片的透射率,并结合K-K关系式和相应的理论模型,反演求解颗粒的光谱复折射率的研究方法。该方法可以使得颗粒保持自然状态,并且实验所需设备简单、精度高、操作方便和适应范围广等优点,但是K-K关系式需要考虑全波长范围的积分,而实验中的激光器通常都是有限波长范围,这就增加了实验测量和理论计算的难度。散射法是通过测量激光照射到颗粒系上后,不同方向的散射光的强度,然后结合相应的理论模型求解颗粒的光谱复折射率的反问题研究方法。该方法能很好的保持颗粒的自然状态,但是确定某一波长下的光谱复折射率至少需要测三个不相关方向的散射量,这就增加实验测量过程的难度。
综上所述现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。
发明内容
本发明为了解决现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。提出了一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法。
本发明所述一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用脉冲宽度为tp的波长为λ的角频率为ω的单频调制激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,10-12s≤tp≤10-9s,0<θc<π/2,400MHz<ω<800MHz,0.3μm<λ<2.5μm;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数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(λ),通过对辐射传输方程求解,获得计算域内的辐射强度场Ic(ω,λ);
步骤九、通过公式:
式中I0是频域激光的辐射强度峰值;
步骤十、利用步骤三中光电接收器测量获得的左侧边界处的频域半球反射信号和右侧边界处的频域半球透射信号与步骤九中获得左侧边界的频域半球反射信号的估计值和右侧边界的频域半球透射信号的估计值获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值和对待测颗粒的吸收因子k取值,作为结果,完成利用单频调制激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
本发明通过建立球形颗粒系频域反射信号和频域透射信号测量的正问题和逆问题求解模型,采用粒度分析仪测量得到颗粒系的粒径分布情况,最后基于这些信号结合逆问题求解技术获得颗粒的光谱复折射率在已知颗粒其他物性参数的前提下,结合频域辐射传输模型,运用Mie理论结合微粒群优化算法反演获得球形颗粒光谱复折射率的方法,该方法解决了现有测量过程复杂和精确度低的问题,且与现有方法相比精确度提高了10%。
附图说明
图1为具体实施方式一所述一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法测量示意图;
图2为具体实施方式一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法的流程图。
具体实施方式
具体实施方式一、结合图1、图2说明本实施方式,本实施方式所述一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用脉冲宽度为tp的波长为λ的角频率为ω的单频调制激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,10-12s≤tp≤10-9s,0<θc<π/2,400MHz<ω<800MHz,0.3μm<λ<2.5μm;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数P(D);
步骤五、利用待测颗粒光谱复折射率方程m(λ)=n(λ)+ik(λ),通过待测颗粒的折射率范围对该待测物质的折射率n进行取值,通过待测颗粒的吸收因子的取值范围对该待测颗粒的吸收因子k进行取值,获得该种待测颗粒单个颗粒的光谱复折射率m(λ);其中i为虚数单位;(煤灰颗粒的n∈[1.18,1.92],k∈[0.01,1.13])
步骤六、通过现有Mie理论,利用步骤五获得的该种待测颗粒的光谱复折射率m(λ),获得单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ);
步骤七、利用步骤六获得的单个颗粒的光谱吸收截面Cabs,pred(D,λ)和光谱散射截面Csca,pred(D,λ),获得颗粒系的吸收系数κa(λ)和散射系数κs(λ);
步骤八、利用步骤七获得颗粒系的吸收系数κa(λ)和散射系数κs(λ),通过对辐射传输方程求解,获得计算域内的辐射强度场Ic(ω,λ);
步骤九、通过公式:
式中I0是频域激光的辐射强度峰值;
步骤十、利用步骤三中光电接收器测量获得的左侧边界处的频域半球反射信号和右侧边界处的频域半球透射信号与步骤九中获得左侧边界的频域半球反射信号的估计值和右侧边界的频域半球透射信号的估计值获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值和对待测颗粒的吸收因子k取值,作为结果,完成利用单频调制激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
本发明提出了利用单频调制激光辐照技术的颗粒系光谱复折射率的测量方法,采用单频调制激光,相比于时域脉冲激光而言,获取辐射信号在某段时间内的特征值(包括幅值和相位两个信息)能够较容易实现;采用Mie理论模型,该模型是非偏振平面电磁波透射均质球形粒子时得到的Maxwell方程远场解,能很精确的反应出粒子的电磁散射特性;采用微粒群优化算法,该算法求解优化问题时有简单、高效和灵敏度高等优点;该项发明为研究球形颗粒光谱复折射率的测量方法,对航天、国防和工业生产等具有十分重要的意义。
Iin(λ)——为激光的峰值强度;
x——指方形脉冲在介质中传播后沿与入射边界的距离;
c——光速;
δ——狄拉克函数;
β(λ)——衰减系数;
θ,θc——分别为光在颗粒系中传播方向和激光在颗粒系外的入射方向。
具体实施方式三、本实施方式是对具体实施方式一所述的一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤七中获得颗粒系的吸收系数κa(λ)和散射系数κs(λ)为通过公式:
计算获得,式中,D为所测颗粒的粒径,Dmax和Dmin分别为所测颗粒的最大粒径和最小粒径,N0所测颗粒的总个数。
具体实施方式四、本实施方式是对具体实施方式一所述的一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤八中通过对辐射传输方程求解,获得计算域内的辐射强度场的方法为:
利用瞬态辐射传输方程:
c——光速
θ,θ'——分别指光在颗粒系中的传输方向和散射光的入射方向;
Φ(θ',θ)——从θ′方向入射并从θ方向散射出去的散射相函数。
具体实施方式五、本实施方式是对具体实施方式一所述的一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法的进一步说明,步骤十中获得目标函数Fobj通过公式:
实现。
Claims (5)
1.一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,该方法的具体步骤为:
步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;
步骤二、利用脉冲宽度为tp的波长为λ的角频率为ω的单频调制激光沿着与样本容器表面法线成θc角的方向入射到样本容器左侧表面;其中,10-12≤tp≤10-9s,0<θc<π/2,400MHz<ω<800MHz,0.3μm<λ<2.5μm;
步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数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(λ),通过对辐射传输方程求解,获得计算域内的辐射强度场Ic(ω,λ);
步骤九、通过公式:
式中I0是频域激光的辐射强度峰值;
步骤十、利用步骤三中光电接收器测量获得的左侧边界处的频域半球反射信号和右侧边界处的频域半球透射信号与步骤九中获得左侧边界的频域半球反射信号的估计值和右侧边界的频域半球透射信号的估计值获得目标函数Fobj;
步骤十一、判断步骤十中的获得目标函数Fobj是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值和对待测颗粒的吸收因子k取值,作为结果,完成利用单频调制激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
3.根据权利要求1所述的一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,步骤七中获得颗粒系的吸收系数κa(λ)和散射系数κs(λ)为通过公式:
计算获得,式中,D为所测颗粒的粒径,Dmax和Dmin分别为所测颗粒的最大粒径和最小粒径,N0所测颗粒的总个数。
5.根据权利要求1所述的一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,步骤十中获得目标函数Fobj通过公式:
实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310467043.6A CN103471968B (zh) | 2013-10-09 | 2013-10-09 | 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310467043.6A CN103471968B (zh) | 2013-10-09 | 2013-10-09 | 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103471968A true CN103471968A (zh) | 2013-12-25 |
CN103471968B CN103471968B (zh) | 2015-08-19 |
Family
ID=49796909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310467043.6A Active CN103471968B (zh) | 2013-10-09 | 2013-10-09 | 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103471968B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501214A (zh) * | 2016-10-21 | 2017-03-15 | 厦门大学嘉庚学院 | 基于实折射角法的吸收性介质的复折射率测量方法 |
CN106644852A (zh) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学 | 基于超短脉冲激光辐照同时获取球形颗粒光学常数与粒径分布的测量方法 |
CN107144501A (zh) * | 2017-06-13 | 2017-09-08 | 南京航空航天大学 | 利用短脉冲激光测量球形颗粒分形聚集特征参数的方法 |
CN109211837A (zh) * | 2018-08-21 | 2019-01-15 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157282A1 (en) * | 2004-01-21 | 2005-07-21 | Asml Netherlands B.V. | Optical element, lithographic apparatus comprising such optical element and device manufacturing method |
CN101055241A (zh) * | 2007-05-10 | 2007-10-17 | 中国科学院安徽光学精密机械研究所 | 多通道散射反演微颗粒三维形状检测仪及检测方法 |
-
2013
- 2013-10-09 CN CN201310467043.6A patent/CN103471968B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157282A1 (en) * | 2004-01-21 | 2005-07-21 | Asml Netherlands B.V. | Optical element, lithographic apparatus comprising such optical element and device manufacturing method |
CN101055241A (zh) * | 2007-05-10 | 2007-10-17 | 中国科学院安徽光学精密机械研究所 | 多通道散射反演微颗粒三维形状检测仪及检测方法 |
Non-Patent Citations (3)
Title |
---|
林莉等: "《哈尔滨地区气溶胶粒子光学常数的实验研究》", 《节能技术》, vol. 31, no. 4, 31 July 2013 (2013-07-31) * |
邢键等: "《悬浊液光谱透射法测量烟尘粒子的复折射率》", 《光谱学与光谱分析》, vol. 30, no. 12, 31 December 2010 (2010-12-31) * |
阮立明等: "《采用GMM方法分析聚集粒子的辐射特性》", 《哈尔滨工业大学学报》, vol. 41, no. 3, 31 March 2009 (2009-03-31) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106644852A (zh) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学 | 基于超短脉冲激光辐照同时获取球形颗粒光学常数与粒径分布的测量方法 |
CN106501214A (zh) * | 2016-10-21 | 2017-03-15 | 厦门大学嘉庚学院 | 基于实折射角法的吸收性介质的复折射率测量方法 |
CN107144501A (zh) * | 2017-06-13 | 2017-09-08 | 南京航空航天大学 | 利用短脉冲激光测量球形颗粒分形聚集特征参数的方法 |
CN107144501B (zh) * | 2017-06-13 | 2019-07-09 | 南京航空航天大学 | 一种利用短脉冲激光测量球形颗粒分形聚集特征参数的方法 |
CN109211837A (zh) * | 2018-08-21 | 2019-01-15 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
CN109211837B (zh) * | 2018-08-21 | 2020-12-25 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103471968B (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487356B (zh) | 一种基于短脉冲激光透反射信号的球形颗粒光谱复折射率测量方法 | |
Baran | From the single-scattering properties of ice crystals to climate prediction: A way forward | |
CN100451678C (zh) | 高光谱全偏振三维成像集成探测系统 | |
CN103471968B (zh) | 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法 | |
CN104634705B (zh) | 一种基于连续激光的球形颗粒光谱复折射率与颗粒系粒径分布的获得方法 | |
CN109061668B (zh) | 一种多视场偏振激光雷达探测系统及探测冰云的方法 | |
CN104792672A (zh) | Pm2.5质量浓度的高度分布探测装置及其探测方法 | |
CN103398964A (zh) | 一种基于腔增强技术的气体探测方法 | |
CN105988124A (zh) | 基于吸收消光比标定廓线的气溶胶污染物分布探测方法 | |
Gupta et al. | Microwave response of rough surfaces with auto-correlation functions, RMS heights and correlation lengths using active remote sensing | |
CN105203437B (zh) | 基于前向散射多角度测量的球形颗粒光学常数与颗粒系粒径分布的同时重构方法 | |
Bao et al. | Study of a laser echo in an inhomogeneous dust environment with a continuous field Monte Carlo radiative transfer model | |
CN103472033A (zh) | 一种利用连续激光辐照技术测量球形颗粒光谱复折射率的方法 | |
CN107144501B (zh) | 一种利用短脉冲激光测量球形颗粒分形聚集特征参数的方法 | |
CN107271336B (zh) | 基于连续激光多角度散射测量的球形颗粒分形聚集特征参数的反演方法 | |
CN103454244B (zh) | 基于多频调制激光辐照的半透明介质辐射特性测量方法 | |
Cui et al. | The study on THz wave propagation feature in atmosphere | |
Xu et al. | Calculating model for equivalent consumption efficiency in polarization measurement system of oil-spilled on the sea | |
Atemkeng et al. | Radiative transfer technique for retrieving the radiative properties of agricultural soils | |
Guo et al. | Laser backscattering of multi-scaled large particles based on superimposed scattering | |
Yang et al. | A Seven‐Parameter BRDF Model with Double‐Peak Characteristic Suitable for Sandy Soil | |
Golmohammady et al. | Scrutinizing the effects of turbulence on partially coherent Gaussian beam propagating through Caspian Sea | |
CN203772414U (zh) | 一种微型光谱仪 | |
CN107300514B (zh) | 利用单频调制激光辐照技术测量球形颗粒分形聚集特征参数的方法 | |
Li et al. | Progress of airborne lidar of ocean chlorophyll observations including algorithm and instruments |
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 |