发明名称 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法
摘要 一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,属于颗粒光学特性测量技术领域。解决了现有的球形颗粒光谱复折射率测量方法仍然存在测量过程复杂和精确度低的问题。通过建立球形颗粒系频域反射信号和频域透射信号测量的正问题和逆问题求解模型,采用粒度分析仪测量得到颗粒系的粒径分布情况,最后基于这些信号结合逆问题求解技术获得颗粒的光谱复折射率在已知颗粒其他物性参数的前提下,结合频域辐射传输模型,运用Mie理论结合微粒群优化算法反演获得球形颗粒光谱复折射率的方法。本发明适用于测量颗粒的光谱复折射率。
申请公布号 CN103471968A 申请公布日期 2013.12.25
申请号 CN201310467043.6 申请日期 2013.10.09
申请人 哈尔滨工业大学 发明人 齐宏;贺振宗;任亚涛;孙双成;阮立明;谈和平
分类号 G01N15/00(2006.01)I 主分类号 G01N15/00(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 岳泉清
主权项 1.一种利用单频调制激光辐照技术测量球形颗粒光谱复折射率的方法,其特征在于,该方法的具体步骤为:步骤一、将待测颗粒装在有机玻璃样本容器中,将待测颗粒溶于溶液中,并搅拌均匀,使颗粒系处于悬浮流动状态;步骤二、利用脉冲宽度为t<sub>p</sub>的波长为λ的角频率为ω的单频调制激光沿着与样本容器表面法线成θ<sub>c</sub>角的方向入射到样本容器左侧表面;其中,10<sup>-12</sup>≤t<sub>p</sub>≤10<sup>-9</sup>s,0&lt;θ<sub>c</sub>&lt;π/2,400MHz&lt;ω&lt;800MHz,0.3μm&lt;λ&lt;2.5μm;步骤三、采用光电接收器在样本颗粒容器的左侧表面测量频域半球反射信号,在样本颗粒容器的右侧表面测量频域半球透射信号,并记录样本颗粒的左侧边界处的时频域半球反射信号<img file="FDA0000392892870000014.GIF" wi="172" he="62" />和右侧边界处的频域半球透射信号<img file="FDA0000392892870000015.GIF" wi="203" he="69" />步骤四、利用粒径分析仪测量样本颗粒系的粒径,并根据粒径的分布规律获得粒径分布函数P(D);步骤五、利用待测颗粒光谱复折射率方程m(λ)=n(λ)+ik(λ),通过待测颗粒的折射率范围对该待测物质的折射率n进行取值,通过待测颗粒的吸收因子的取值范围对该待测颗粒的吸收因子k进行取值,获得该种待测颗粒单个颗粒的光谱复折射率m(λ);其中i为虚数单位;步骤六、通过现有Mie理论,利用步骤五获得的该种待测颗粒的光谱复折射率m(λ),获得单个颗粒的光谱吸收截面C<sub>abs,pred</sub>(D,λ)和光谱散射截面C<sub>sca,pred</sub>(D,λ);步骤七、利用步骤六获得的单个颗粒的光谱吸收截面C<sub>abs,pred</sub>(D,λ)和光谱散射截面C<sub>sca,pred</sub>(D,λ),获得颗粒系的吸收系数κ<sub>a</sub>(λ)和散射系数κ<sub>s</sub>(λ);步骤八、利用步骤七获得颗粒系的吸收系数κ<sub>a</sub>(λ)和散射系数κ<sub>s</sub>(λ),通过对辐射传输方程求解,获得计算域内的辐射强度场I<sub>c</sub>(ω,λ);步骤九、通过公式:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>R</mi><mo>^</mo></mover><mrow><mi>pred</mi><mo>,</mo><mi>&lambda;</mi></mrow></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>I</mi><mn>0</mn></msub></mfrac><msubsup><mo>&Integral;</mo><mrow><mi>&pi;</mi><mo>/</mo><mn>2</mn></mrow><mi>&pi;</mi></msubsup><msub><mover><mi>I</mi><mo>^</mo></mover><mrow><mi>d</mi><mo>,</mo><mi>&lambda;</mi></mrow></msub><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>cos</mi><mi></mi><mi>&theta;</mi><mi>sin</mi><mi>&theta;d&theta;</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mi>pred</mi><mo>,</mo><mi>&lambda;</mi></mrow></msub><mrow><mo>(</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>I</mi><mn>0</mn></msub></mfrac><mo>[</mo><mn>2</mn><mi>&pi;</mi><msubsup><mo>&Integral;</mo><mn>0</mn><mrow><mi>&pi;</mi><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mover><mi>I</mi><mo>^</mo></mover><mrow><mi>d</mi><mo>,</mo><mi>&lambda;</mi></mrow></msub><mrow><mo>(</mo><mi>L</mi><mo>,</mo><mi>&theta;</mi><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mi>cos</mi><mi></mi><mi>&theta;</mi><mi>sin</mi><mi>&theta;d&theta;</mi><mo>+</mo><msub><mover><mi>I</mi><mo>^</mo></mover><mrow><mi>c</mi><mo>,</mo><mi>&lambda;</mi></mrow></msub><mrow><mo>(</mo><mi>L</mi><msub><mrow><mo>,</mo><mi>&theta;</mi></mrow><mi>c</mi></msub><mo>,</mo><mi>&omega;</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>获得左侧边界的频域半球反射信号的估计值<img file="FDA0000392892870000013.GIF" wi="188" he="79" />和右侧边界的频域半球透射信号的估计值<img file="FDA0000392892870000021.GIF" wi="206" he="79" />式中I<sub>0</sub>是频域激光的辐射强度峰值;<img file="FDA0000392892870000022.GIF" wi="218" he="77" />为波长为λ角频率为ω的激光,在θ方向上x=0处的左侧边界上散射光的辐射强度,θ为辐射方向角;<img file="FDA0000392892870000023.GIF" wi="227" he="76" />为波长为λ角频率为ω的激光,θ方向上x=L处的右侧边界上散射光的辐射强度,θ为辐射方向角;<img file="FDA0000392892870000024.GIF" wi="240" he="79" />为波长为λ角频率为ω的激光,频域激光沿着入射方向θ<sub>c</sub>衰减到样本右侧壁面时的辐射强度,θ<sub>c</sub>为时频域激光入射方向角,L为样本厚度;步骤十、利用步骤三中光电接收器测量获得的左侧边界处的频域半球反射信号<img file="FDA0000392892870000025.GIF" wi="182" he="79" />和右侧边界处的频域半球透射信号<img file="FDA0000392892870000026.GIF" wi="204" he="79" />与步骤九中获得左侧边界的频域半球反射信号的估计值<img file="FDA0000392892870000027.GIF" wi="195" he="79" />和右侧边界的频域半球透射信号的估计值<img file="FDA0000392892870000028.GIF" wi="207" he="79" />获得目标函数F<sub>obj</sub>;步骤十一、判断步骤十中的获得目标函数F<sub>obj</sub>是否小于设定阈值ξ,若是,则将步骤五中对待测物质的折射率n的取值和对待测颗粒的吸收因子k取值,作为结果,完成利用单频调制激光辐照技术对球形颗粒光谱复折射率的测量,否则返回步骤五。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号