留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

不同湿度环境下可见光波段激光偏振特性研究

战俊彤 张肃 付强 段锦 李英超 姜会林

战俊彤, 张肃, 付强, 段锦, 李英超, 姜会林. 不同湿度环境下可见光波段激光偏振特性研究[J]. 红外与激光工程, 2020, 49(9): 20200057. doi: 10.3788/IRLA20200057
引用本文: 战俊彤, 张肃, 付强, 段锦, 李英超, 姜会林. 不同湿度环境下可见光波段激光偏振特性研究[J]. 红外与激光工程, 2020, 49(9): 20200057. doi: 10.3788/IRLA20200057
Zhan Juntong, Zhang Su, Fu Qiang, Duan Jin, Li Yingchao, Jiang Huilin. Laser polarization characteristics of visible light band in different humidity environments[J]. Infrared and Laser Engineering, 2020, 49(9): 20200057. doi: 10.3788/IRLA20200057
Citation: Zhan Juntong, Zhang Su, Fu Qiang, Duan Jin, Li Yingchao, Jiang Huilin. Laser polarization characteristics of visible light band in different humidity environments[J]. Infrared and Laser Engineering, 2020, 49(9): 20200057. doi: 10.3788/IRLA20200057

不同湿度环境下可见光波段激光偏振特性研究

doi: 10.3788/IRLA20200057
基金项目: 国家自然科学基金(61705017,61905025,61890960);吉林省科学技术厅项目(20190103156JH);吉林省教育厅项目(JJKH20181140KJ,JJHK20181089KJ)
详细信息
    作者简介:

    战俊彤(1987-),女,讲师,博士,主要从事偏振特性分析方面的研究。Email:zhanjuntong@cust.edu.cn

    通讯作者: 张肃(1985-),女,讲师,博士,主要从事偏振成像探测方面的研究。Email:susiezhang21@126.com
  • 中图分类号: O436.3

Laser polarization characteristics of visible light band in different humidity environments

  • 摘要: 雾霾天气严重干扰了可见光成像效果,利用可见光偏振特性可以有效提升探测效率。雾霾环境会受到气溶胶颗粒湿度的影响,湿度是雾霾环境重要物理参数。为了获得可见光在雾霾环境的偏振特性规律,分析了环境湿度对偏振特性的影响。在非偏振光气溶胶单粒子散射特性基础上,采用改进蒙特卡罗方法建立了偏振传输模型,对不同湿度水雾环境下可见波段偏振光传输特性进行研究,重点分析水雾环境湿度改变对不同可见光波段偏振光偏振特性的影响情况并建立了接近真实水雾环境,通过室内实验对偏振模型进行验证,对不同湿度环境下的 450、532 、671 nm 线偏振光与圆偏振光的偏振度与偏振态改变进行了比较与分析,仿真模型置信度>60%。研究结果表明:偏振光的偏振度随水雾环境湿度是呈下降趋势。随着波长的增大,偏振度变化趋于平缓,出射偏振度随着波长的增加而变大,当激光波长为 450、532、671 nm 时,偏振度下降点的湿度值分别为 50%、70%、90%;圆偏振光入射,圆偏振光的旋性对于偏振度没有影响,且高于线偏振度值。对于水雾这种容易受湿度影响较大的环境而言,在可见光波段,应该尽量选择较长波长的偏振光进行传输探测,因为湿度对较短波长的影响比较长波长的影响要大,所以,在湿度较大的环境中,较长波长的圆偏振光是偏振保持特性最好的。在湿度较大的环境中,应尽量选取波长较长的偏振光成像,以达到较好的成像效果。
  • 图  1  参考平面的定义

    Figure  1.  Definition of reference plane

    图  2  波长450 nm偏振度与湿度关系图

    Figure  2.  Relationship between the polarization degree and humidity of 450 nm polarized light

    图  3  波长532 nm偏振度与湿度关系图

    Figure  3.  Relationship between the polarization degree and humidity of 532 nm polarized light

    图  4  波长671 nm偏振度与湿度关系图

    Figure  4.  Relationship between the polarization degree and humidity of 671 nm polarized light

    图  5  系统组成图

    Figure  5.  System composition diagram

    图  6  系统实物图

    Figure  6.  Real system

    图  7  波长450 nm偏振度与湿度关系实验图

    Figure  7.  Relationship between the polarization degree and humidity of 450 nm polarized light experiment

    图  8  波长532 nm偏振度与湿度关系实验图

    Figure  8.  Relationship between the polarization degree and humidity of 532 nm polarized light experiment

    图  9  波长671 nm偏振度与湿度关系实验图

    Figure  9.  Relationship between the polarization degree and humidity of 671 nm polarized light experiment

    表  1  水雾环境粒子湿度与折射率的关系

    Table  1.   Relationship between humidity and refractive index of particles in water mist environment

    Relative
    humidity
    Particle
    radius /μm
    Index of
    real part
    Index of
    imaginary part
    30% 0.548 1.369 1.0298E-4
    40% 0.570 1.365 1.3907E-4
    50% 0.597 1.361 1.7635E-4
    60% 0.632 1.356 2.1518E-4
    70% 0.681 1.352 2.5608E-4
    80% 0.755 1.347 3.0002E-4
    90% 0.902 1.341 3.4926E-4
    下载: 导出CSV
  • [1] Pilinis C, Pandis S N, Seinfeld J H. Sensitivity of direct climate forcing by atmospheric aerosols to aerosol size and composition [J]. Journal of Geophysical Research, 1995, 100: 18739-18754. doi:  10.1029/95JD02119
    [2] Boesche Eyk, Stammes Piet, Bennartz Ralf. Aerosol influence on polarization and intensity in near-infrared O<sub>2</sub> and CO<sub>2</sub> absorption bands observed from space [J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2009, 110: 223-239.
    [3] 杨利红, 柯熙政, 马冬冬. 偏振激光在大气传输中的退偏研究[J]. 光电工程, 2008, 35(11): 62-67. doi:  10.3969/j.issn.1003-501X.2008.11.013

    Yang Lihong, Ke Xizheng, Ma Dongdong. Depolarization characteristics of the polarized laser in atmosphere [J]. Opto-Electronic Engineering, 2008, 35(11): 62-67. (in Chinese) doi:  10.3969/j.issn.1003-501X.2008.11.013
    [4] 蔡嘉, 高隽, 范之国, 等. 湿度对气溶胶单粒子偏振特性的影响研究[J]. 传感器与微系统, 2013, 32(9): 35-38. doi:  10.3969/j.issn.1000-9787.2013.09.010

    Cai Jia, Gao Jun, Fan Zhiguo, et al. Research on influence of humidity on polarization properties of aerosol single particles [J]. Transducer and Microsystem Technologies, 2013, 32(9): 35-38. (in Chinese) doi:  10.3969/j.issn.1000-9787.2013.09.010
    [5] 郝增周, 龚芳, 潘德炉, 等. 沙尘气溶胶粒子群的散射和偏振特性[J]. 光学学报, 2012, 32(1): 0101002. doi:  10.3788/AOS201232.0101002

    Hao Zengzhou, Gong Fang, Pan Delu, et al. Scattering and polarization characteristics of dust aerosol particles [J]. Acta Optica Sinica, 2012, 32(1): 0101002. (in Chinese) doi:  10.3788/AOS201232.0101002
    [6] 战俊彤, 付强, 张肃, 等. 前向散射角度对偏振成像对比度影响的研究[J]. 光学学报, 2015, 35(7): 0711007.

    Zhan Juntong, Fu Qiang, Zhang Su, et al. Research on the effects of forward scattering angle to polarization imaging contrast [J]. Acta Optica Sinica, 2015, 35(7): 0711007. (in Chinese)
    [7] 杨蔚, 顾国华, 陈钱, 等. 红外偏振图像的目标检测方法[J]. 红外与激光工程, 2014, 43(8): 2746-2751. doi:  10.3969/j.issn.1007-2276.2014.08.059

    Yang Wei, Gu Guohua, Chen Qian, et al. Method of target detection for infrared polarization image [J]. Infrared and Laser Engineering, 2014, 43(8): 2746-2751. (in Chinese) doi:  10.3969/j.issn.1007-2276.2014.08.059
    [8] 付跃刚, 朱启凡, 张玉慧, 等. 仿生虾蛄眼偏振特性结构研究[J]. 红外与激光工程, 2019, 48(8): 0827001. doi:  10.3788/IRLA201948.0827001

    Fu Yuegang, Zhu Qifan, Zhang Yuhui., et al. Research on the polarization structures inspired by mantis shrimp eyes [J]. Infrared and Laser Engineering, 2019, 48(8): 0827001. (in Chinese) doi:  10.3788/IRLA201948.0827001
  • [1] 张鑫, 宗思光, 李斌, 余扬.  浑浊水体激光后向散射特性仿真与实验研究 . 红外与激光工程, 2023, 52(1): 20220280-1-20220280-7. doi: 10.3788/IRLA20220280
    [2] 胡玮娜, 吕勇, 耿蕊, 李宇海, 牛春晖.  光电探测器表面损伤状态偏振成像式探测系统 . 红外与激光工程, 2022, 51(6): 20210629-1-20210629-9. doi: 10.3788/IRLA20210629
    [3] 陈世杰, 牛春晖, 李晓英, 吕勇.  CCD损伤进程中光学成像系统猫眼回波特性研究 . 红外与激光工程, 2021, 50(9): 20200425-1-20200425-6. doi: 10.3788/IRLA20200425
    [4] 战俊彤, 张肃, 付强, 段锦, 李英超.  非球型粒子对激光偏振特性的影响 . 红外与激光工程, 2020, 49(11): 20200150-1-20200150-6. doi: 10.3788/IRLA20200150
    [5] 王辉, 王进, 李校博, 胡浩丰, 刘铁根.  一种基于圆偏光的偏振去雾成像优化方法 . 红外与激光工程, 2019, 48(11): 1126001-1126001(5). doi: 10.3788/IRLA201948.1126001
    [6] 王稼禹, 李英超, 史浩东, 江伦, 王超, 刘壮, 李冠霖.  折反式望远系统全视场全口径偏振特性研究 . 红外与激光工程, 2019, 48(3): 318004-0318004(8). doi: 10.3788/IRLA201948.0318004
    [7] 谢绍禹, 叶茂, 赵毅强, 吴环宝, 贾晓东.  偏振探测激光引信目标判别技术 . 红外与激光工程, 2018, 47(S1): 5-9. doi: 10.3788/IRLA201847.S106002
    [8] 李金金, 孙晓兵, 康晴, 李树, 殷玉龙.  偏振光谱仪偏振探测精度分析 . 红外与激光工程, 2018, 47(1): 123002-0123002(7). doi: 10.3788/IRLA201847.0123002
    [9] 李小路, 李昀晔, 谢鑫浩, 徐立军.  基于叶片含水量的激光偏振成像模型研究 . 红外与激光工程, 2017, 46(11): 1106004-1106004(6). doi: 10.3788/IRLA201746.1106004
    [10] 肖文健, 马东玺, 陈志斌, 张勇, 肖程, 秦梦泽.  大尺寸空间角测量系统光轴指向不确定度评定 . 红外与激光工程, 2016, 45(11): 1118004-1118004(7). doi: 10.3788/IRLA201645.1118004
    [11] 林芬芳, 张东彦, 王秀, 吴太夏, 陈新福.  基于偏振光谱的叶片尺度下玉米与杂草识别研究 . 红外与激光工程, 2016, 45(12): 1223001-1223001(10). doi: 10.3788/IRLA201645.1223001
    [12] 张肃, 付强, 段锦, 战俊彤.  湿度对偏振光传输特性影响的研究 . 红外与激光工程, 2016, 45(5): 511001-0511001(6). doi: 10.3788/IRLA201645.0511001
    [13] 林勇, 徐智勇, 汪井源, 宋超, 王荣, 耿常锁.  雾环境下非视距散射光通信最佳链路分析 . 红外与激光工程, 2015, 44(2): 705-710.
    [14] 于洵, 朱磊, 姜旭, 武继安, 李建强.  基于微透镜阵列偏振探测器的噪声性能研究 . 红外与激光工程, 2015, 44(S1): 189-195.
    [15] 颛孙晓博, 武文远, 黄雁华, 龚艳春, 吴成国, 李兆兆.  基于MB模型的简化偏振BRDF模型建立与仿真 . 红外与激光工程, 2015, 44(3): 1098-1102.
    [16] 杨雨迎, 崔占忠, 张万君.  基于双高斯拟合的目标激光偏振特性实验研究 . 红外与激光工程, 2015, 44(2): 449-454.
    [17] 陶会荣, 张福民, 曲兴华.  无合作目标测量中目标表面后向散射特性的实验研究 . 红外与激光工程, 2014, 43(S1): 95-100.
    [18] 王霞, 夏润秋, 金伟其, 刘敬, 梁建安.  红外偏振成像探测技术进展 . 红外与激光工程, 2014, 43(10): 3175-3182.
    [19] 崔骥, 付欢, 刘烨, 王清华, 李振华, 贺安之.  典型粗糙表面激光后向散射特性实验 . 红外与激光工程, 2013, 42(5): 1208-1211.
    [20] 秦刚, 杨郁, 张建生.  船舰远程尾流散射光偏振特性的蒙特卡洛模拟 . 红外与激光工程, 2013, 42(7): 1730-1736.
  • 加载中
图(9) / 表(1)
计量
  • 文章访问数:  1243
  • HTML全文浏览量:  821
  • PDF下载量:  80
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-25
  • 修回日期:  2020-03-18
  • 网络出版日期:  2020-05-25
  • 刊出日期:  2020-09-22

不同湿度环境下可见光波段激光偏振特性研究

doi: 10.3788/IRLA20200057
    作者简介:

    战俊彤(1987-),女,讲师,博士,主要从事偏振特性分析方面的研究。Email:zhanjuntong@cust.edu.cn

    通讯作者: 张肃(1985-),女,讲师,博士,主要从事偏振成像探测方面的研究。Email:susiezhang21@126.com
基金项目:  国家自然科学基金(61705017,61905025,61890960);吉林省科学技术厅项目(20190103156JH);吉林省教育厅项目(JJKH20181140KJ,JJHK20181089KJ)
  • 中图分类号: O436.3

摘要: 雾霾天气严重干扰了可见光成像效果,利用可见光偏振特性可以有效提升探测效率。雾霾环境会受到气溶胶颗粒湿度的影响,湿度是雾霾环境重要物理参数。为了获得可见光在雾霾环境的偏振特性规律,分析了环境湿度对偏振特性的影响。在非偏振光气溶胶单粒子散射特性基础上,采用改进蒙特卡罗方法建立了偏振传输模型,对不同湿度水雾环境下可见波段偏振光传输特性进行研究,重点分析水雾环境湿度改变对不同可见光波段偏振光偏振特性的影响情况并建立了接近真实水雾环境,通过室内实验对偏振模型进行验证,对不同湿度环境下的 450、532 、671 nm 线偏振光与圆偏振光的偏振度与偏振态改变进行了比较与分析,仿真模型置信度>60%。研究结果表明:偏振光的偏振度随水雾环境湿度是呈下降趋势。随着波长的增大,偏振度变化趋于平缓,出射偏振度随着波长的增加而变大,当激光波长为 450、532、671 nm 时,偏振度下降点的湿度值分别为 50%、70%、90%;圆偏振光入射,圆偏振光的旋性对于偏振度没有影响,且高于线偏振度值。对于水雾这种容易受湿度影响较大的环境而言,在可见光波段,应该尽量选择较长波长的偏振光进行传输探测,因为湿度对较短波长的影响比较长波长的影响要大,所以,在湿度较大的环境中,较长波长的圆偏振光是偏振保持特性最好的。在湿度较大的环境中,应尽量选取波长较长的偏振光成像,以达到较好的成像效果。

English Abstract

    • 大量细颗粒物的排放容易引起的雾霾天气,可见光在这种混浊介质中传输时由于散射,会产生偏振特性,能够降低不良天气状况对可见光传输的影响。因此需要开展可见光偏振传输特性的研究。雾霾环境会受到气溶胶颗粒湿度的影响,湿度是雾霾环境物理参数之一,湿度的变化会导致凝聚现象的发生,从而改变了粒子的物理特性。湿度与颗粒半径,颗粒折射率有很大的关系,半径与折射率的变化最终会影响散射光的偏振特性[1]

      关于偏振光传输特性,国外学者Sebastian将蒙特卡罗方法应用到均一稳定的大气环境偏振光传输[2],对散射特性进行了分析。西安理工大学杨利红分析了10.6 μm波段激光在浓雾(20~30 μm)与薄雾(1~10 μm)环境下偏振角变化关系,偏振激光在弱雾中传输时,偏振角变化相比于浓雾中更为规则,在浓雾中退偏角度更大[3]。但上述研究蒙特卡洛计算方法中直接通过改变粒径大小来控制湿度,并未包含对湿度与偏振特性的数学表达式。合肥工业大学蔡嘉对于0.5 μm波段非偏振光在气溶胶单粒子传输中的散射偏振特性进行了仿真研究,发现在散射角120°~150°时气溶胶单粒子偏振度随湿度变化明显,130°~140°时,随湿度增加而增加[4]。国家海洋局郝增周根据Mie散射理论对0.55 μm波段非偏振光在沙尘粒子群环境的偏振散射特性进行了仿真研究,发现沙尘粒子单次反照率与偏振特性都随着湿度的增加而增大,后向散射高于前向散射,后向散射受湿度影响较大[5]。由此可见,环境湿度影响着气溶胶粒子的粒径和分布变化,从而影响在粒子中传输的偏振光特性,因此深入研究环境湿度与偏振特性的关系对于偏振探测效率的提升十分关键[6-8]

      针对以上问题,文中在非偏振光气溶胶单粒子散射特性基础上,采用蒙特卡罗方法建立了偏振光粒子群传输特性模型,对不同湿度水雾环境下可见波段偏振光传输特性进行研究,重点分析水雾环境湿度改变对不同可见光波段偏振光偏振特性的影响情况。并建立了接近真实水雾环境的半实物仿真系统,通过室内实验对偏振光传输特性模型进行验证,该方法既适合变化湿度环境又适合于特性湿度环境的偏振光特性研究。

    • 湿度与颗粒半径,颗粒折射率有很大的关系,半径与折射率的变化最终会影响散射光的偏振特性,干洁环境中颗粒半径与湿度关系如下:

      $$\frac{{{r_h}}}{{{r_0}}} = {\left( {1 - f} \right)^{ - 1/u}}$$ (1)

      式中:${r_h}$为受湿度影响变化的粒子半径;${r_0}$为干洁环境中的颗粒半径;$f$为湿度;u为常数,对于雾天环境而言,气溶胶吸湿性强,u=3.9;${r_0}$为初始半径。

      湿度导致了粒径的变化,而粒径变化又会导致复折射率变化:

      $${m_{re}} = {m_{rw}} + \left( {{m_{r0}} - {m_{rw}}} \right){\left[ {\frac{{{r_h}}}{{{r_0}}}} \right]^{ - 3}}$$ (2)
      $$\frac{{{m_{ie}}}}{{m_{re}^2 + 2}} = \dfrac{{{m_{iw}}}}{{m_{rw}^2}} + \left( {\frac{{{m_{i0}}}}{{m_{r0}^2 + 2}} - \frac{{{m_{iw}}}}{{m_{rw}^2}}} \right){\left[ {\frac{{{r_h}}}{{{r_0}}}} \right]^{ - 3}}$$ (3)

      ${m_{re}}$为复折射实部;${m_{ie}}$为虚部。rie,0和w分别表示实部、虚部吸湿后气溶胶,干气溶胶粒子和水。湿度影响下的气溶胶折射率为:

      $${m_e} = {m_{re}} + {m_{ie}}$$ (4)

      水雾环境中粒子尺度较大,通常都在0.1~10 μm之间,水雾粒子等效半径0.5 μm,折射率为1.33+1.96E-4i,通过公式(1)、(4)对不同湿度的粒子物理参数进行了计算,如表1所示。

      表 1  水雾环境粒子湿度与折射率的关系

      Table 1.  Relationship between humidity and refractive index of particles in water mist environment

      Relative
      humidity
      Particle
      radius /μm
      Index of
      real part
      Index of
      imaginary part
      30% 0.548 1.369 1.0298E-4
      40% 0.570 1.365 1.3907E-4
      50% 0.597 1.361 1.7635E-4
      60% 0.632 1.356 2.1518E-4
      70% 0.681 1.352 2.5608E-4
      80% 0.755 1.347 3.0002E-4
      90% 0.902 1.341 3.4926E-4

      表1可以看出:水雾粒子折射率实部随着湿度的增加而变小,水雾环境折射率实部值基本不变。由于水雾粒子吸湿性强,所以粒子折射率改变值小。

    • 为了得到湿度与偏振光传输特性的关系,定义参考平面如图1所示。首先,确定光子传播方向,夹角θ为散射角,x轴与z轴和散射平面分别构成了散射前后的参考平面,散射前后的斯托克斯分量由垂直和平行于这两个平面的EE分量定义,每次散射发生时,都必须对斯托克斯分量进行调整以使参考平面为新的参考面。

      图  1  参考平面的定义

      Figure 1.  Definition of reference plane

      参考平面定义以后,设定激光波长、粒子直径、粒子复折射率等,散射系数${\mu _s}$,吸收系数${\mu _a}$,大气的消光系数${\mu _t} = {\mu _s} + {\mu _a}$,不同偏振态的光的斯托克斯参量初始值等。光子沿z轴正方向入射,初始位置${u_0}$(0,0,0),初始方向余弦${D_0}$(0,0,1),偏振参考平面为X轴与Z轴构成的平面,定义光子初始斯托克斯参量,并且先将散射后各斯托克斯参量都设为0。

      粒子传输的过程的步长就是光子的自由程, 其中$\xi $为关于0.5对称分布的随机数,消光系数与光子自由程${S_1}$关系为:

      $${S_1} = \ln \left( \xi \right){\mu _t}$$ (5)

      一旦光子传输的自由程确定后,光子将发生移动,当前位置$\left( {x,y,z} \right)$和传输方向$\left( {{\mu _x},{\mu _y},{\mu _z}} \right)$确定光子到达的下一个散射点的坐标$\left( {x',y',z'} \right)$为:

      $$\left\{ {\begin{aligned} & {x' = x + {\mu _x}{S_1}} \\ & {y' = y + {\mu _y}{S_1}} \\ & {z' = z + {\mu _z}{S_1}} \end{aligned}} \right.$$ (6)

      并且根据坐标位置,判断光子是否射出参考平面所确定的边界,如果射出参考平面确定的边界,则对光子进行斯托克斯参量计算。反之,必须先对斯托克斯分量进行调整以使参考平面为新的子午面,再重新进行粒子半径选择。

      光子与粒子碰撞后,散射角$\alpha $和方位角$\beta $抽样由联合概率密度函数(PDF)得到,该函数与入射光斯托克斯分量[I0Q0U0V0]T的关系为:

      $$\rho \left( {\alpha ,\beta } \right) = {m_{11}}\left( \alpha \right) + {m_{12}}\left( \alpha \right)\left[ {{Q_0}\cos \left( {2\beta } \right) + {U_0}\sin \left( {2\beta } \right)} \right]/{I_0}$$ (7)

      ${m_{11}}\left( \alpha \right)$${m_{12}}\left( \alpha \right)$分别为球形粒子Mueller矩阵$M\left( \alpha \right)$中相应元素。

      $$M\left( \alpha \right) = \left[ {\begin{array}{*{20}{c}} {{m_{11}}\left( \alpha \right)}&{{m_{12}}\left( \alpha \right)}&0&0 \\ {{m_{12}}\left( \alpha \right)}&{{m_{11}}\left( \alpha \right)}&0&0 \\ 0&0&{{m_{33}}\left( \alpha \right)}&{{m_{34}}\left( \alpha \right)} \\ 0&0&{ - {m_{34}}\left( \alpha \right)}&{{m_{33}}\left( \alpha \right)} \end{array}} \right]$$ (8)

      ${m_{11}}\left( \alpha \right)$${m_{12}}\left( \alpha \right)$${m_{33}}\left( \alpha \right)$${m_{34}}\left( \alpha \right)$与散射幅度值${S_1}$,${S_2}$有如下关系:

      $$\begin{split} & {m_{11}}\left( \alpha \right) = \frac{1}{2}\left( {{{\left| {{S_1}} \right|}^2} + {{\left| {{S_2}} \right|}^2}} \right),\;\;{{m_{12}}\left( \alpha \right) = \frac{1}{2}\left( {{{\left| {{S_1}} \right|}^2} - {{\left| {{S_2}} \right|}^2}} \right)} \\ & {m_{33}}\left( \alpha \right) = \frac{1}{2}\left( {{S_1}S_2^ * - S_{_1}^ * {S_2}} \right),\;\;{{m_{44}}\left( \alpha \right) = \frac{i}{2}\left( {{S_1}S_2^ * - S_2^ * {S_1}} \right)} \end{split}$$ (9)

      由于散射幅度值${S_1}$,${S_2}$表达式为:

      $$\left\{ {\begin{aligned} & {{S_1}(\theta ) = \sum\limits_{n = 1}^\infty {\frac{{2 n + 1}}{{n(n + 1)}}} \left[ {{a_n}\frac{{P_n^1(\cos \theta )}}{{\sin \theta }} + {b_n}\frac{{P_n^1(\cos \theta )}}{{d\theta }}} \right]}\\ & {{S_2}(\theta ) = \sum\limits_{n = 1}^\infty {\frac{{2 n + 1}}{{n(n + 1)}}} \left[ {{a_n}\frac{{P_n^1(\cos \theta )}}{{d\theta }} + {b_n}\frac{{P_n^1(\cos \theta )}}{{\sin \theta }}} \right]} \end{aligned}} \right.$$ (10)

      式中:anbn为Mie散射系数,它们与球形颗粒的大小,折射率有很大关系,其表达式为:

      $$\left\{ {\begin{aligned} & {{a_n} = \frac{{\psi _n^\prime (mx){\psi _n}(x) - m{\psi _n}(mx)\psi _n^\prime (x)}}{{\psi _n^\prime (mx){\xi _n}(x) - m{\psi _n}(mx)\xi _n^\prime (x)}}}\\ & {{b_n} = \frac{{m\psi _n^\prime (mx){\psi _n}(x) - {\psi _n}(mx)\psi _n^\prime (x)}}{{m\psi _n^\prime (mx){\xi _n}(x) - {\psi _n}(mx)\xi _n^\prime (x)}}} \end{aligned}} \right.$$ (11)

      式中:${\psi _n}\left( x \right)$${\xi _n}\left( x \right)$为Riccati-Bessel函数,可用第一类Bessel球函数${J_{n + 1/2}}\left( x \right)$和半整数阶第二类Hankel函数$H_{n + 1/2}^{\left( 2 \right)}$表示:

      $$\left\{ {\begin{aligned} & {{\psi _n}(x) = \sqrt {\pi x/2} {J_{n + 1/2}}(x)}\\ & {{\xi _n}(x) = \sqrt {\pi x/2} H_{n + 1/2}^{(2)}(x)} \end{aligned}} \right.$$ (12)

      所以,相应的散射强度函数为:

      $${i_1}(\theta ) = {\left| {{S_1}(\theta )} \right|^2} = {\left| {\sum\limits_{n = 1}^\infty {\frac{{2 n + 1}}{{n(n + 1)}}} \left( {{a_n}{\pi _n} + {b_n}{\tau _n}} \right)} \right|^2}$$ (13)
      $${i_2}(\theta ) = {\left| {{S_2}(\theta )} \right|^2} = {\left| {\sum\limits_{n = 1}^\infty {\frac{{2 n + 1}}{{n(n + 1)}}} \left( {{a_n}{\pi _n} + {b_n}{\tau _n}} \right)} \right|^2}$$ (14)

      光子经过n次散射后的能量权重变为:

      $${W_n} = {W_{n - 1}} \cdot {\mu _s}/\left( {{\mu _s} + {\mu _a}} \right)$$ (15)

      当光子能量权重低于某一阈值(一般为10−4)或飞离边界时,光子传输终止。光子从多分散系边界射出时,其斯托克斯矢量必须经过最后一次旋转以保证参考平面与探测器所在的平面相同,旋转的角度为:

      $$\omega = \pm {\rm{arctan}}\left( {{\mu _{{y}}}/{\mu _{{x}}}} \right)$$ (16)

      式中:反射模式取正号,透射取负号。由于散射路径的不同,光子到达探测器的时间不同,对于偏振分量形如 [I(t),Q(t),U(t),V(t)]T的光束,时域偏振度定义为:

      $${\mathop{\rm DOP}\nolimits} (t) = \frac{{\sqrt {{Q^2}(t) + {U^2}(t) + {V^2}(t)} }}{{I(t)}}$$ (17)

      光束经散射后总的偏振度定义为:DOP=(Q2+U2+V2)1/2/I,其中,IQUV表示不同时刻到达探测器的光子偏振分量的累加值。

      将公式(1)、(2)、(3)代入公式(11),其中f为湿度,得到湿度影响下的Mie散射系数为:

      $$\left\{ {\begin{aligned} & {a{{\left( f \right)}_n} = \frac{{\psi _n'(m(f)x(f)){\psi _n}(x(f)) - m(f){\psi _n}(m(f)x(f))\psi _n'(x(f))}}{{\psi _n'(m(f)x(f)){\xi _n}(x(f)) - m(f){\psi _n}(m(f)x(f))\xi _n'(x(f))}}}\\ & {b{{(f)}_n} = \frac{{m(f)\psi _n'(m(f)x(f)){\psi _n}(x(f)) - {\psi _n}(m(f)x(f))\psi _n'(x(f))}}{{m(f)\psi _n'(m(f)x(f)){\xi _n}(x(f)) - {\psi _n}(m(f)x(f))\xi _n'(x(f))}}} \end{aligned}} \right.$$ (18)

      再将公式(18)代入公式(10),可以得到湿度影响下的散射振幅表达式:

      $$\left\{ {\begin{array}{*{20}{l}} {{S_1}(\theta ) = \displaystyle\sum_{{n} = 1}^\infty {\dfrac{{2{n} + 1}}{{{n}({n} + 1)}}} \left[ {a{{(f)}_{n}}\dfrac{{P_n^1(\cos \theta )}}{{\sin \theta }} + b{{({\rm{f}})}_{n}}\dfrac{{P_n^1(\cos \theta )}}{{{\rm{d}}\theta }}} \right]}\\ {{S_2}(\theta ) = \displaystyle\sum_{{n} = 1}^\infty {\dfrac{{2{n} + 1}}{{{n}({n} + 1)}}} \left[ {a{{(f)}_{n}}\dfrac{{P_n^1(\cos \theta )}}{{d\theta }} + b{{({\rm{f}})}_{n}}\dfrac{{P_n^1(\cos \theta )}}{{\sin \theta }}} \right]} \end{array}} \right.$$ (19)

      再将公式(19)代入公式(9),可以计算出Mueller矩阵中相应元素,再乘以入射偏振光,最终通过公式(17)即可得到改进的雾霾粒子湿度与偏振传输特性的关系模型。

    • 采用蒙特卡洛模拟光子散射过程,计算了可见光波段的3种波长(450、532、671 nm) 6种入射偏振态的偏振光(0°线偏振光、90°线偏振光、45°线偏振光、135°线偏振光、左旋圆偏振光、右旋圆偏振光)经过水雾粒子传输后的偏振光特性,计算曲线如图2~4所示。

      图  2  波长450 nm偏振度与湿度关系图

      Figure 2.  Relationship between the polarization degree and humidity of 450 nm polarized light

      图  3  波长532 nm偏振度与湿度关系图

      Figure 3.  Relationship between the polarization degree and humidity of 532 nm polarized light

      图  4  波长671 nm偏振度与湿度关系图

      Figure 4.  Relationship between the polarization degree and humidity of 671 nm polarized light

      图2~4中可以看出:偏振度随水雾环境湿度是呈下降趋势,随着波长的增大,偏振度变化趋于平缓,出射偏振度变大。当激光波长分为450、532、671 nm时,偏振度下降点时的湿度为50%、70% 、90%。圆偏振光经过水雾粒子散射后偏振度改变情况一致。圆偏振光受湿度影响情况比线偏振光要小,在图中表现为比线偏振光变化平缓。这是因为对于水雾环境中的粒子而言,短波情况下前向散射大于后向散射,所以,出射偏振光的偏振度要大于波长较长的偏振光。当波长增加到671 nm时,相比于其他两种波长变化趋于平缓,这是由于波长的增加,使得$x = 2\pi /\lambda $变小,即粒子尺度参数变小,使得Muller矩阵中${m_{11}}$${m_{12}}$${m_{13}}$${m_{14}}$变小,这样会导致散射光强度的变小,偏振度改变不大,也就显示为震荡不明显。而且由于水雾环境的强吸湿性,使得粒子尺度明显增大,导致了粒子的散射作用增强,故而随着湿度的增加水雾环境粒子的偏振度不断下降。

    • 室内测试搭建实验系统如图5所示,系统实物图如图6所示。发射端组成为激光器、衰减片、偏振片、波片、滤光片组成。激光器输出功率为50 mW的固体激光器,根据实验波长来选择相应激光器,偏振片透过波段为400~700 nm,各波长透过率可以达到80%以上。接收端分成两路一路由波片、偏振片、光功率计组成,一路由偏振态测量仪组成。

      图  5  系统组成图

      Figure 5.  System composition diagram

      图  6  系统实物图

      Figure 6.  Real system

      接收端偏振角度为0°、60°、120°,采集的强度信息通过下式计算出相应的偏振信息,并与偏振态测量仪数据做比对。

      $$\left\{\begin{aligned} & {{{I}} = \dfrac{1}{2}\left( {{{{I}}^\prime }\left( {60,45,\dfrac{\pi }{2}} \right) + {{{I}}^\prime }\left( {120,45,\dfrac{\pi }{2}} \right) - 3{{{I}}^\prime }\left( {0,45,\dfrac{\pi }{2}} \right)} \right)}\\ & Q = 2{I^\prime }(0,0,0) - \dfrac{1}{2}{I^\prime }(60,45,0) - \dfrac{1}{2}{I^\prime }\left( {120,45,\dfrac{\pi }{2}} \right) \\ & \quad\quad + \dfrac{3}{2}{I^\prime }\left( {0,45,\dfrac{\pi }{2}} \right)\\ & {U = \dfrac{{2\sqrt 3 }}{3}\left( {{{{I}}^\prime }\left( {60,45,\dfrac{\pi }{2}} \right) - {{{I}}^\prime }\left( {120,45,\dfrac{\pi }{2}} \right)} \right)}\\ & {{{V}} = \dfrac{1}{2}\left( {{{{I}}^\prime }\left( {60,45,\dfrac{\pi }{2}} \right) + {{{I}}^\prime }\left( {120,45,\dfrac{\pi }{2}} \right) - {{{I}}^\prime }\left( {0,45,\dfrac{\pi }{2}} \right)} \right)} \end{aligned}\right.$$ (20)

      该测量方法增加了圆偏振分量,实现全偏振度测量,弥补了以往仅针对线偏振分量进行偏振特性测量造成的误差。

    • 通过对450、532、671 nm的6种偏振态在不同湿度的水雾环境中实验数据分析可知:随着水雾环境湿度的增加,偏振度是不断降低的,而且圆偏振光与线偏振光在低湿度30%−40%情况下较为接近,随着湿度的增加,圆偏振光的偏振度明显高于线偏振光。

      图7~图9中每一种偏振态都是随着波长的增大而增大,由此可见对于可见光波段水雾环境,较长波长偏振特性保持良好,即在湿度较大的环境中,应尽量选取波长较长的偏振光成像,以达到较好的成像效果。

      图  7  波长450 nm偏振度与湿度关系实验图

      Figure 7.  Relationship between the polarization degree and humidity of 450 nm polarized light experiment

      图  8  波长532 nm偏振度与湿度关系实验图

      Figure 8.  Relationship between the polarization degree and humidity of 532 nm polarized light experiment

      图  9  波长671 nm偏振度与湿度关系实验图

      Figure 9.  Relationship between the polarization degree and humidity of 671 nm polarized light experiment

      为了验证仿真与实验的一致性,依据公式(21)计算了仿真与实测的偏振度的置信度。

      $$M = 1 - \left( {\sum\limits_1^n {\frac{{\left| {R - {R_m}} \right|}}{{{R_m}}}} } \right)/n \times 100\% $$ (21)

      式中:RRm分别为同一条件下的仿真和实测得到的偏振度值,由于对7种湿度进行了仿真与实测实验,所以n=7,分别计算了6种偏振态的偏振度置信度,450 nm的0°、90°、45°、135°、左旋、右旋偏振光的置信度分别为63.9%、64.1%、64.2%、64.1%、75.4%、75.5%;532 nm的置信度分别为73.1%、73.2%、73.1%、73%、82.6%、82.5%。671 nm的置信度分别为74.0%、74.1%、74.4%、74.1%、85.4%、85.3%。仿真模型置信度均大于60%,结果可信。

    • 采用改进蒙特卡洛的仿真方法计算了不同湿度环境对可见光偏振特性的影响情况,并与通过室内实验分析偏振光传输特性影响。在450 、532 、671 nm三个波段条件下,偏振光的偏振度随水雾环境湿度是呈下降趋势,但是随着波长的增大,偏振度变化趋于平缓,出射偏振度随着波长的增加而变大。圆偏振光入射,圆偏振光的旋性对于偏振度没有影响,且高于线偏振度值。因此,对于水雾这种容易受湿度影响较大的环境而言,在可见光波段,应该尽量选择较长波长的偏振光进行传输探测,因为湿度对较短波长的影响比较长波长的影响要大,所以,在湿度较大的环境中,较长波长的圆偏振光是偏振保持特性最好的。在湿度较大的环境中,应尽量选取波长较长的偏振光成像,以达到较好的成像效果。接下来需要开展多波段条件下的不同粒子偏振传输特性研究工作,丰富偏振特性规律,为雾天的偏振成像探测提供更有利的理论依据。

参考文献 (8)

目录

    /

    返回文章
    返回