留言板

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

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

环境相对湿度对光学镜面散射特性的影响

陈壮 田博宇 缪麟 孙年春 张彬

陈壮, 田博宇, 缪麟, 孙年春, 张彬. 环境相对湿度对光学镜面散射特性的影响[J]. 红外与激光工程, 2021, 50(6): 20210041. doi: 10.3788/IRLA20210041
引用本文: 陈壮, 田博宇, 缪麟, 孙年春, 张彬. 环境相对湿度对光学镜面散射特性的影响[J]. 红外与激光工程, 2021, 50(6): 20210041. doi: 10.3788/IRLA20210041
Chen Zhuang, Tian Boyu, Miao Lin, Sun Nianchun, Zhang Bin. Effect of ambient humidity on scattering characteristics of optical mirror surface[J]. Infrared and Laser Engineering, 2021, 50(6): 20210041. doi: 10.3788/IRLA20210041
Citation: Chen Zhuang, Tian Boyu, Miao Lin, Sun Nianchun, Zhang Bin. Effect of ambient humidity on scattering characteristics of optical mirror surface[J]. Infrared and Laser Engineering, 2021, 50(6): 20210041. doi: 10.3788/IRLA20210041

环境相对湿度对光学镜面散射特性的影响

doi: 10.3788/IRLA20210041
基金项目: 四川省科技计划项目(2018JY0553);中国科学院自适应光学重点实验室基金(LAOF1801)
详细信息
    作者简介:

    陈壮,男,硕士生,主要从事光学系统优化与杂光分析方面的研究

    张彬,女,教授,博士生导师,博士,主要从事激光与光通信、光电子技术等方面研究

  • 中图分类号: O436.2

Effect of ambient humidity on scattering characteristics of optical mirror surface

  • 摘要: 在成像光学系统中,环境相对湿度会影响光学镜面污染粒子的散射特性,进而对成像光学系统性能造成影响。根据吸湿性,将污染粒子分为亲水性粒子和疏水性粒子。利用镜面污染粒子的光散射模型,分析了亲水性粒子散射效率因子随环境湿度的变化规律。在此基础上,结合亲水性和疏水性污染粒子的统计模型,分析了多粒子污染下,不同种类粒子占比和镜面洁净度等参数对光学镜面散射特性的影响。结果表明,污染粒子造成光学镜面的双向散射分布函数(BRDF)和全积分散射值(TIS)均随着环境相对湿度的增加而增大。并且,亲水性粒子比例越高,光学镜面散射损耗越强。随着光学镜面洁净度等级的减小,环境相对湿度对镜面散射损耗的影响有所减小。
  • 图  1  不同表面洁净度镜面上的粒子分布

    Figure  1.  Particle distribution with different surface cleanliness in the mirror surface

    图  2  (NH4)2SO4粒子散射增长因子随环境相对湿度变化的对比曲线

    Figure  2.  Comparison curve of (NH4)2SO4 particle scattering growth factor with environmental relative humidity

    图  3  不同粒径污染粒子散射效率因子Qsca随湿度变化

    Figure  3.  Scattering efficiency factor Qsca of polluted particles of different particle size varies with humidity

    图  4  不同条件下BRDF随入射角度的变化 (λ=0.55 μm)

    Figure  4.  Variation of BRDF with incident angle under different conditions (λ=0.55 μm)

    表  1  不同种类粒子占比时的TIS (λ=0.55 μm G1=300)

    Table  1.   TIS for surfaces with different proportions of hydrophilic particles (λ=0.55 μm G1=300)

    Humidity00.30.60.9
    TIS(A:30%)7.01×10−57.14×10−57.54×10−51.85×10−4
    TIS(A:40%)7.01×10−57.18×10−57.72×10−52.23×10−4
    TIS(A:50%)7.01×10−57.23×10−57.90×10−52.61×10−4
    TIS(A:60%)7.01×10−57.27×10−58.08×10−52.99×10−4
    TIS(A:70%)7.01×10−57.32×10−58.25×10−53.37×10−4
    TIS(A:80%)7.01×10−57.37×10−58.43×10−53.76×10−4
    下载: 导出CSV

    表  2  不同镜面洁净度等级时的TIS (λ=0.55 μm)

    Table  2.   TIS for different mirror cleanliness levels (λ=0.55 μm)

    Mirror class300500700
    TIS(XRH=0)7.01×10−58.10×10−44.56×10−3
    TIS(XRH=0.3)7.27×10−58.63×10−44.73×10−3
    TIS(XRH=0.6)8.08×10−51.02×10−35.26×10−3
    TIS(XRH=0.9)2.99×10−45.23×10−31.95×10−2
    下载: 导出CSV

    表  3  不同角度入射时的TIS

    Table  3.   TIS at different angles of incidence (λ=0.55 μm G1=300)

    Humidity00.30.60.9
    TIS(0°)7.02×10−57.62×10−58.97×10−53.71×10−4
    TIS(5°)7.02×10−57.61×10−58.97×10−53.71×10−4
    TIS(10°)7.02×10−57.61×10−58.96×10−53.71×10−4
    TIS(15°)7.01×10−57.61×10−58.96×10−53.71×10−4
    TIS(30°)6.99×10−57.59×10−58.93×10−53.70×10−4
    TIS(50°)6.95×10−57.54×10−58.89×10−53.70×10−4
    下载: 导出CSV
  • [1] 张斌. 南京市大气降尘特征及源解析[D]. 南京: 南京大学, 2013.

    Zhang Bin. Characteristics and source apportionment of atmospheric dustfall in Nanjing city[D]. Nanjing: Nanjing University, 2013. (in Chinese)
    [2] Huang Cong, You Xinghai, Zhang Bin. Influence of surface cleanliness of optical element on its surface scattering characteristics [J]. Infrared and Laser Engineering, 2019, 48(1): 0120002. (in Chinese)
    [3] Zhang Zhuo, Wang Weihe, Wang Houmao, et al. Effects of humidity on absorbing aerosol index [J]. Journal of Remote Sensing, 2019, 23(6): 1177-1185. (in Chinese)
    [4] Ma Jing, Zhou Funan, Zhou Yanping, et al. Optical surface pollution detection device for dual-use star and earth [J]. Optical and Precision Engineering, 2016, 24(8): 1878-1883. (in Chinese) doi:  10.3788/OPE.20162408.1878
    [5] Spyak P R, Wolfe W L. Scatter from particulate-contaminated mirrors. part 2: theory and experiment for dust and λ=0.6328 μm [J]. Optical Engineering, 1992, 31(8): 1757-1763. doi:  10.1117/12.58709
    [6] Tang I N. Chemical and size effects of hydroscopic aerosols on light scattering coefficients [J]. Journal of Geopgysical Research: Atmospheres, 1996, 101(D14): 19245-19250. doi:  10.1029/96JD03003
    [7] Zhao Jiajia, Gu Fang, Zhang Jiahong, et al. Hygroscopic growth model and scattering characteristics of single-particle aerosols [J]. Acta Optical Sinica, 2020, 40(5): 0501001. (in Chinese)
    [8] Zhang Xuehai, Dai Congming, Zhang Xin, et al. Effect of relative humigity and particle shapr on the scattering properties of sea saly aerosols [J]. Infrared and Laser Engineering, 2019, 48(8): 0809002. (in Chinese)
    [9] Zhang Zhicha, Ni Changjian, Zhang Chengyu, et al. Relationship between particle size hygroscopic growth and scattering hygroscopic growth [J]. China Enviromental Science, 2020, 40(12): 5198-5204. (in Chinese) doi:  10.3969/j.issn.1000-6923.2020.12.011
    [10] Mao Qianjun, Wang Li. Scattering hygroscopic growth characteristics of typical haze particles and clusters [J]. Acta Optica Sinica, 2020, 40(24): 2429001. (in Chinese)
    [11] 崔芬萍. 南京地区大气气溶胶光学特性的观测研究[D]. 南京: 南京信息工程大学, 2017.

    Cui Fenping. Observational study on optical properties of aerosol in Nanjing[D]. Nanjing: Nanjing University of Information Science and Technology, 2017. (in Chinese)
    [12] 倪刘建. 南京市不同功能区大气降尘性质及其对土壤影响的研究[D]. 南京: 南京农业大学, 2006.

    Ni Liujian. A dissertation submitted to Nanjing agricultural university for the degree of master of environmental engineering[D]. Nanjing: Nanjing Agricultural University, 2006. (in Chinese)
    [13] Wang Jiacheng, Qiao Yanli, Yang Shizhi, et al. Study on the relation of ocean atmospheric aerosol optical model parameters [J]. Acta Optica Sinica, 2011, 31(11): 1101005. (in Chinese)
    [14] Ma Wangjiehui, Liu Yanlei, Chen Zhiying, et al. Near-infrared BRDF of material surfaces at varying temperatures [J]. Chinese Optics, 2020, 13(5): 1115-1123. (in Chinese) doi:  10.37188/CO.2019-0256
    [15] 王瑞雪. 随机表面粗糙度对热辐射特性影响的研究[D]. 北京: 北京科技大学, 2009.

    Wang Ruixue. A study on the heat radiant characteristics influenced by the random surface roughness[D]. Beijing: Beijing Science & Technology University, 2009. (in Chinese)
  • [1] 王建峰, 姜晓军, 李陶然, 张晓明, 葛亮.  雷达定标卫星光学散射特性研究 . 红外与激光工程, 2020, 49(1): 0105003-0105003(8). doi: 10.3788/IRLA202049.0105003
    [2] 战俊彤, 张肃, 付强, 段锦, 李英超.  非球型粒子对激光偏振特性的影响 . 红外与激光工程, 2020, 49(11): 20200150-1-20200150-6. doi: 10.3788/IRLA20200150
    [3] 王翔, 张科鹏, 陈壮, 张彬, 陈坚, 赵建华.  光学镜面污染颗粒吸收特性的影响分析 . 红外与激光工程, 2020, 49(4): 0414004-0414004-7. doi: 10.3788/IRLA202049.0414004
    [4] 黄聪, 游兴海, 张彬.  光学元件表面洁净度对其表面散射特性的影响 . 红外与激光工程, 2019, 48(1): 120002-0120002(7). doi: 10.3788/IRLA201948.0120002
    [5] 王菲菲, 李学彬, 郑显明, 张文忠, 罗涛, 朱文越, 成巍, 邓志武.  相对湿度和风速对海洋大气气溶胶粒子谱的影响 . 红外与激光工程, 2019, 48(S1): 83-88. doi: 10.3788/IRLA201948.S117003
    [6] 张学海, 戴聪明, 张鑫, 魏合理, 朱希娟, 马静.  相对湿度和粒子形态对海盐气溶胶粒子散射特性的影响 . 红外与激光工程, 2019, 48(8): 809002-0809002(8). doi: 10.3788/IRLA201948.0809002
    [7] 徐强, 王东琴, 王旭, 吴振森.  应用T矩阵法对大气灰霾简单非球形粒子散射特性的计算与分析 . 红外与激光工程, 2017, 46(11): 1117003-1117003(7). doi: 10.3788/IRLA201746.1117003
    [8] 杨玉峰, 李挺, 李建勋, 王昭雷.  沙尘在FSO常用红外波段的散射特性研究 . 红外与激光工程, 2017, 46(6): 604004-0604004(6). doi: 10.3788/IRLA201746.0604004
    [9] 张肃, 付强, 段锦, 战俊彤.  湿度对偏振光传输特性影响的研究 . 红外与激光工程, 2016, 45(5): 511001-0511001(6). doi: 10.3788/IRLA201645.0511001
    [10] 潘永强, 陈佳.  光学薄膜减散射特性研究 . 红外与激光工程, 2016, 45(1): 118007-0118007(5). doi: 10.3788/IRLA201645.0118007
    [11] 熊兴隆, 刘雯箐, 符超, 蒋立辉.  考虑多次散射的卷云几何特征和光学特性反演方法 . 红外与激光工程, 2015, 44(9): 2734-2741.
    [12] 张学海, 魏合理.  基于用几何光学和米散射法的球形粒子前向散射特性计算研究 . 红外与激光工程, 2015, 44(5): 1485-1490.
    [13] 李忠洋, 谭联, 袁源, 邴丕彬, 袁胜, 徐德刚, 姚建铨.  太赫兹波光学参量效应放大特性的理论研究 . 红外与激光工程, 2014, 43(8): 2650-2655.
    [14] 孙腾飞, 张骏, 吕海兵, 袁晓东, 曹增辉, 郑田甜.  光学镜面污染对激光传输特性的影响 . 红外与激光工程, 2014, 43(5): 1444-1448.
    [15] 于双双, 杜吉, 史宣.  激光粒度仪光学系统设计方法 . 红外与激光工程, 2014, 43(6): 1735-1739.
    [16] 李浩, 陈晓颖, 单陈华, 唐丽萍.  雾和气溶胶的光学特性对前向散射和总散射能见度仪性能的影响 . 红外与激光工程, 2013, 42(6): 1568-1574.
    [17] 韩意, 孙华燕.  空间目标光学散射特性研究进展 . 红外与激光工程, 2013, 42(3): 758-766.
    [18] 汪之国, 龙兴武, 王飞, 谭中奇.  环形腔中背向散射特性的理论研究 . 红外与激光工程, 2013, 42(2): 339-343.
    [19] 王莲芬, 赵选科, 左翔, 王金金, 孙红辉.  高斯激光束TEM00模散射信号模拟与分析 . 红外与激光工程, 2013, 42(11): 2940-2943.
    [20] 黄朝军, 吴振森, 刘亚锋.  1.06 μm激光气溶胶凝聚粒子散射特性 . 红外与激光工程, 2013, 42(9): 2353-2357.
  • 加载中
图(4) / 表(3)
计量
  • 文章访问数:  27
  • HTML全文浏览量:  9
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-22
  • 修回日期:  2021-03-22
  • 刊出日期:  2021-06-30

环境相对湿度对光学镜面散射特性的影响

doi: 10.3788/IRLA20210041
    作者简介:

    陈壮,男,硕士生,主要从事光学系统优化与杂光分析方面的研究

    张彬,女,教授,博士生导师,博士,主要从事激光与光通信、光电子技术等方面研究

基金项目:  四川省科技计划项目(2018JY0553);中国科学院自适应光学重点实验室基金(LAOF1801)
  • 中图分类号: O436.2

摘要: 在成像光学系统中,环境相对湿度会影响光学镜面污染粒子的散射特性,进而对成像光学系统性能造成影响。根据吸湿性,将污染粒子分为亲水性粒子和疏水性粒子。利用镜面污染粒子的光散射模型,分析了亲水性粒子散射效率因子随环境湿度的变化规律。在此基础上,结合亲水性和疏水性污染粒子的统计模型,分析了多粒子污染下,不同种类粒子占比和镜面洁净度等参数对光学镜面散射特性的影响。结果表明,污染粒子造成光学镜面的双向散射分布函数(BRDF)和全积分散射值(TIS)均随着环境相对湿度的增加而增大。并且,亲水性粒子比例越高,光学镜面散射损耗越强。随着光学镜面洁净度等级的减小,环境相对湿度对镜面散射损耗的影响有所减小。

English Abstract

    • 降尘作为大气气溶胶的主要组成成分,是指在空气动力作用下自然沉降到地面的大气悬浮颗粒物。大气降尘的两个主要来源是降尘粒径小于10 μm的细粒物源(大气环境粉尘污染的背景)和短距离悬移搬运的10~50 μm粗粒自然源[1]。考虑到大口径封窗的制造难度以及封窗对红外波段的截止特性,大口径光学系统通常倾向于使用敞开式结构,导致大量元件表面曝露于复杂外界环境中。大气降尘粒子在重力作用下会沉降到光学表面而造成镜面污染,进而影响光学系统的性能[2]

      在镜面散射特性研究中,常采用Mie理论。然而,现有Mie理论多针对干燥粒子进行分析[2],而真实应用环境中的环境湿度会对其有所影响。以南京地区著名的紫金山天文望远镜为例,南京市全年相对湿度维持在一个较高水平,相对湿度范围在40%以上的时间占全年的93.48%,夏季相对湿度最大而冬季相对湿度最小。受环境湿度的影响,降尘粒子的粒径和形状会发生改变[3]。同时,吸附在降尘粒子外的水层也会改变降尘粒子的光学常数,从而影响降尘的散射特性[3]。对于这类暴露在湿度较高环境中的大型光学系统,忽略环境湿度的影响将难以准确评估其光学系统性能。

      光学元件表面散射特性的相关研究已有诸多报道。国内外学者基于理论和实验分析了不同镜面污染颗粒对镜面散射分布的影响[4-5];也有学者对受湿度影响的单粒子散射特性进行了深入研究[6~8]。目前,鲜有报道有关外界环境湿度变化对多污染粒子镜面散射特性造成的影响。文中在考虑相对湿度影响的情况下,给出了Mie理论的修正模型,并以湿度较大的南京地区为例,分析了环境湿度对光学镜面散射特性的影响。

    • 空气中的各种悬浮颗粒附着在镜面上形成表面污染物,这些颗粒污染物可近似视作球形粒子,因而在分析镜面污染粒子光学特性时通常采用Mie理论[2]。大多数镜面散射研究没有考虑环境湿度对污染粒子散射特性的影响,但实际系统中,污染粒子不可避免地会遭受早晚湿度变化、季节降水变化的影响,致使其光学特性会发生变化。为更好地描述污染粒子随环境湿度的变化,笔者基于Mie理论和不同吸湿粒子模型,获得了环境相对湿度影响下的光学镜面散射特性分析模型。

      环境相对湿度(Relative Humidity,XRH)是指单位体积空气内实际所含的水气密度(用d1表示)和同温度下饱和水气密度(用d2表示)的百分比,而另一种计算方法则是实际的空气水气压强(用p1表示)和同温度下饱和水气压强(用p2表示)的百分比[8]

      $${X_{\rm RH}}({\text{%}} )=\frac{{{d_1}}}{{{d_2}}} \times 100{\text{%}} =\frac{{{p_1}}}{{{p_2}}} \times 100{\text{%}} $$ (1)

      随着环境相对湿度的增加,大气中的水蒸气会凝结在污染粒子上,在增加污染粒子粒径尺寸的同时,也逐渐改变其化学成分和复折射率,进而影响粒子的散射效应[7]

      依照粒子吸湿性强弱,通常可分为亲水性和疏水性两类[7]。无机降尘粒子(包括硫酸盐、硝酸盐、铵盐、海盐等)和部分吸湿性粒子属于亲水性污染粒子;而黑碳、沙尘及弱极性和非极性有机物不具有吸湿性,属于疏水性污染粒子[9]。亲水性粒子具有较强的吸湿性,且在高湿度的条件下具有潮解特性,粒子发生潮解时的相对湿度被称为潮解相对湿度(XDRH)。结合已有的实验结果,分两个阶段对亲水性粒子特性进行表征。当XRH<XDRH时,粒子表面吸附水量表示为[7]

      $$\frac{1}{{{X_{\rm RH}}}}=1 + k\frac{{{V_{\rm dry}}}}{{{V_{\rm water}}}}$$ (2)

      XRH=XDRH,粒子发生潮解形成液滴。当时XRH>XDRH,粒子继续吸收水分并被稀释,其粒径与湿度的关系可以表示为[7]

      $$\frac{{{X_{\rm RH}}}}{{\exp \left[ {\dfrac{A}{{{D_{{{\rm dry}}}}G({X_{\rm RH}})}}} \right]}}=\dfrac{{G{{({X_{\rm RH}})}^3} - 1}}{{G{{({X_{\rm RH}})}^3} - (1 - {{t}})}}$$ (3)

      式中:t=0.3为吸湿性参数;Ddry为干燥粒子的直径;A为常量,A=2.09×10−3G(XRH)=D(XRH)/Ddry为湿度增长因子;D(XRH)为随XRH变化的湿粒子直径。

      疏水性粒子的吸湿性较弱,随着环境相对湿度的增加,其表面会吸附一层水膜,疏水性粒子表面的水膜厚度LXRH之间的函数关系为[7]

      $$L={4.34^{1/3}}{\left[ {{\rm ln}\left(\frac{1}{{{X_{\rm RH}}}}\right) + \frac{{1.08}}{{{r_{\rm a}}}}} \right]^{ - 1/3}} \times 0.4$$ (4)

      式中:ra为被吸附粒子的半径。

      Hanel通过大量的实验和理论研究得到等效均匀球形粒子的复折射率me的实部mre、虚部mie、相对湿度XRH和粒径比之间的关系为[8]

      $${m_{\rm re}}={m_{\rm rw}} + ({m_{\rm ro}} - {m_{\rm rw}}) \cdot {[{r_{\rm h}}/{r_{\rm o}}]^{ - 3}}$$ (5)
      $$\frac{{{m_{\rm ie}}}}{{m_{\rm re}^2{\rm{ + }}2}}=\frac{{{m_{\rm iw}}}}{{m_{\rm rw}^2}} + \left(\frac{{{m_{\rm i0}}}}{{m_{\rm ro}^2 + 2}} - \frac{{{m_{\rm iw}}}}{{m_{\rm rw}^2}}\right){\left[\frac{{{r_{\rm h}}}}{{{r_{\rm o}}}}\right]^{ - 3}}$$ (6)

      式中:下标r、i、e、o和w分别表示实部、虚部、吸湿后气溶胶、干燥气溶胶粒子和水。在湿度影响下,气溶胶的复折射率表示为:

      $${m_{\rm e}}={m_{\rm re}} + {m_{\rm ie}}i$$ (7)

      Mie理论的关键参数包括用于计算散射场振幅的Mie系数anbn,以及用于计算内部场振幅的cndn[2]。通过这些系数的组合可以计算粒子的光学特性参数,如散射效率因子、吸收效率因子、不对称因子等。基于公式(8)~(11),可得到考虑环境相对湿度影响下修正后的Mie系数anbncndn分别为:

      $$ {a_{{n}}}{\rm=}\frac{{{m_{\rm{e}}}^2{j_{{n}}}\left( {{m_{\rm e}}x} \right)[xj_{{n}}^{(1)}(x)]' - {j_n}(x)[{m_{\rm e}}x{j_{{n}}}({m_{\rm e}}x)]'}}{{{m_{{\rm e}}}^2{j_{{n}}}({m_{\rm e}}x)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[{m_{\rm e}}x{j_{{n}}}({m_{\rm e}}x)]'}}$$ (8)
      $${b_{{n}}}{\rm=}\frac{{{j_{{n}}}\left( {mx} \right)[x{j_n}(x)]' - {j_n}(x)[{m_{\rm{e}}}x{j_{{n}}}({m_{\rm e}}x)]'}}{{{j_{{n}}}({m_{\rm e}}x)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[{m_{\rm e}}x{j_{{n}}}({m_{\rm e}}x)]'}}$$ (9)
      $${c_{{n}}}{\rm=}\frac{{{j_{{n}}}\left( x \right)[xh_{{n}}^{(1)}(x)]' - h_{{n}}^{(1)}(x)[x{j_{{n}}}({m_{\rm{e}}}x)]'}}{{{j_{{n}}}({m_{\rm e}}x)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[{m_{\rm e}}x{j_{{n}}}({m_{\rm e}}x)]'}}$$ (10)
      $${d_{{n}}}=\frac{{{m_{\rm e}}{j_{{n}}}\left( x \right)[xh_{{n}}^{(1)}(x)]' - {m_{\rm e}}h_{{n}}^{(1)}(x)[x{j_{{n}}}({m_{\rm e}}x)]'}}{{{m_{\rm e}}^2{j_{{n}}}({m_{\rm e}}x)[xh_{{n}}^{(1)}(x)]' - {h_{{n}}}^{(1)}(x)[{m_{\rm e}}x{j_{{n}}}({m_{\rm e}}x)]'}}$$ (11)

      式中:x=ka/2用于表征颗粒的尺寸,为无量纲量,k为波数,a为随湿度变换的球形颗粒直径;me为修正的球形污染粒子相对于周围环境的折射率;z=xz=mxjn(x)和hn(x)分别为Bessel函数和Hankel函数。

      设单色光沿Z轴正向传播,电场振动方向沿X轴,磁场振动方向沿Y轴,散射角θ定义为散射光矢量与Z轴正方向的夹角,方位角φ定义为散射光在XOY平面的投影与X轴的夹角,入射光的强度为I0,散射粒子的直径为d,且污染颗粒满足各向同性,r表示球形粒子到任意观测点P处的直线距离,λ表示波长。于是,自由空间中任意一点P的散射强度可以定义为[2]

      $$I\left( {r,\theta ,\varphi } \right)=\frac{{{\lambda ^2}{I_0}}}{{4{\pi ^2}{r^2}}}\left[ {{i_1}\left( \theta \right){{\sin }^2}\varphi + {i_2}\left( \theta \right){{\cos }^2}\varphi } \right]$$ (12)

      式中:i1(θ)、i2(θ)分别为平行和垂直于散射面的强度分量,称为散射强度函数,它们是振幅函数S1(θ)和S2(θ)模的平方,即

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

      式中:S1S2为散射振幅函数;πnτn为散射角函数。

      值得指出的是,当相对湿度为0时,该修正模型可退化为干燥条件下的Mie散射模型。利用该模型,在参考文献[2]给出的相关参数条件下,其复算结果与文献结果基本一致,从而验证文中所用计算模型和方法的正确性。

    • 在实际应用中,通常使用表面洁净度来度量光学元件的污染程度。根据美国军用标准MIL-STD-209D,镜面的污染程度可用表面洁净度来表征。表面洁净度用一个数值描述,记作G1,表示每平方英尺镜面上大于该数值的颗粒数目刚好为1。根据该标准,镜面上粒子直径与其数量分布近似满足[2]

      $$ \log {N_{\rm a}}=C \cdot \left[ {{{\left( {\log {G_1}} \right)}^2} - {{\left( {\log G} \right)}^2}} \right] $$ (15)

      式中:G为颗粒直径(微米量级且G大于等于1);Na为元件表面每平方英尺上颗粒直径大于等于G的粒子数量;C是归一化常量,对于清洁表面,C的值近似为0.926。

      进一步可推导出每平方毫米的表面上,颗粒直径为G的粒子数目为:

      $$ {F_{\left( G \right)}}=\frac{1}{{{\rm{92\;903}}{\rm{.04}}}}\left[ {{{10}^{C\left( {{{\log }^2}{G_{\rm{l}}} - {{\log }^2}G} \right)}} - {{10}^{C\left( {{{\log }^2}{G_{\rm{l}}} - {{\log }^2}\left( {G + 1} \right)} \right)}}} \right] $$ (16)

      根据公式(16),图1给出了在不同表面洁净度等级下,颗粒尺寸G(μm)与粒子数量F(G)之间的关系。

      图  1  不同表面洁净度镜面上的粒子分布

      Figure 1.  Particle distribution with different surface cleanliness in the mirror surface

      图1可以看出,在同一洁净度等级条件下,镜面上半径越大的污染粒子数量越少。并且,随着镜面洁净度等级的升高,镜面上的污染粒子数增多。

    • 此节利用修正的Mie散射模型分析了湿度变化对单个污染粒子散射效率因子的影响。

      散射效率因子是粒子散射截面与粒子几何截面之比[7],用Qsca表示。对于受湿度影响的包裹性球形粒子,单一污染粒子的散射效率因子为:

      $$ {Q_{{\rm{sca}}}}=\frac{2}{{{x^2}}}\sum\limits_{n=1}^\infty {(2n + 1)} ({\left| {{a_n}} \right|^2} + {\left| {{b_n}} \right|^2}) $$ (17)

      散射吸湿增长因子f(XRH)定义为某一相对湿度下的降尘粒子散射系数与干燥状态下的降尘粒子散射系数的比值[9]

      $${{f(}}{{{X}}_{\rm RH}}{\rm{)=}}\frac{{{Q_{\rm sca}}({X_{\rm RH}})}}{{{Q_{\rm sca}}(X_{\rm dry})}}$$ (18)

      当入射波长为532 nm时,利用修正的Mie模型计算了粒径为0.2 μm的(NH4)2SO4颗粒在不同的XRH时的散射增长因子f(XRH),红点表示的数据参考文献[10],其结果见图2。由图2可知,该模型的计算结果与参考文献[10]结果相符合,说明文中建立的修正Mie散射模型的正确性和有效性。

      仿真时选取湿度较大的南京为典型例,南京地区干燥降尘粒子的复折射率为me=1.368+0.014i[11],利用公式(5)和(6)可分别计算得到折射率实部与虚部随环境相对湿度的变化。由公式(2)~(4)可知污染粒子半径随相对湿度的增大而增大。其中,环境湿度变化对亲水性污染粒子的影响显著,而对疏水性污染粒子影响较小。由于疏水性污染粒子受湿度影响较小,而亲水性污染粒子受环境湿度影响明显,因而笔者重点分析单一亲水性污染粒子散射特性随湿度的变化,此时取XDRH=0.7。图3给出了不同粒径下,亲水性粒子散射效率因子Qsca随湿度的变化曲线。

      图  2  (NH4)2SO4粒子散射增长因子随环境相对湿度变化的对比曲线

      Figure 2.  Comparison curve of (NH4)2SO4 particle scattering growth factor with environmental relative humidity

      图  3  不同粒径污染粒子散射效率因子Qsca随湿度变化

      Figure 3.  Scattering efficiency factor Qsca of polluted particles of different particle size varies with humidity

      图3可以看出,对于粒径相同的污染粒子,散射效率因子随着相对湿度的增加而逐渐增大。在相同的湿度条件下,散射效率因子随着粒径的增大而减小。环境相对湿度在60%以下时,不同粒径条件下的散射系数均增长缓慢,此时,散射效率因子与干燥条件下的散射效率因子相差较小。然而,随着相对湿度的进一步增加,散射效率因子显著变大。导致这一现象的主要原因是亲水性粒子达到潮解相对湿度时发生了潮解所致。McMurry使用前后差分迁移亲水分析仪(Tandem Differential Mobility Analyzer,TDMA)对气溶胶粒子粒径随相对湿度的变化进行了大量研究发现,当环境相对湿度高于60%时,环境湿度的影响变得非常重要[3]。亲水性降尘粒子在吸湿增长过程中会发生潮解,水汽的凝结作用以及固液表面的化学反应和非均向传质过程使得粒径突然增大,致使颗粒物的基本物理状态发生明显的改变。由此可知,相对湿度越大对单个污染粒子造成的散射影响越大,光学损耗加剧,且在相对湿度大于60%的高湿度条件下尤为明显。

    • 近海地区受海洋影响沉积区域上空大气颗粒物以海盐为主,受人类活动和大气环流影响,不同区域降尘粒子成分占比不同[12]。根据近海气溶胶成分和各地降尘成分综合分析,在不同环境条件下,亲水性污染粒子占比大约在30%~80%之间。文中亲疏水粒子统计模型适用于受海洋影响较大的近海地区污染粒子分布。南京市地处中国东部近海地区,受到亚热带季风影响降水充沛。针对南京当地污染粒子成分进行分析[1],文中计算分析了亲水性粒子占比为30%~80%条件下的光学镜面散射情况。

      假设散射到前向方向上并由镜面反射的辐射不会受污染颗粒所影响,来自镜面污染粒子的散射函数(BRDF)可定义[13]为:

      $${\rm{BRDF}}=\frac{{{\lambda ^2}D}}{{4{\pi ^2}}}{i_1}$$ (19)

      式中:D为镜面上的颗粒密度(每平方米上的粒子数);i1为散射强度函数。

      全积分散射TIS (Total Integrated Scatter,TIS)的物理含义是某一散射表面散射到半球空间的辐射通量与入射到该表面的辐射通量的比值,它是BRDF在半球空间的二重积分[2]

      $${\rm TIS}=\int_0^{2\pi } {\int_0^ {\frac{\pi }{2}} {{\rm BSDF}({\theta _{\rm i}},{\phi _{\rm i}};} } {\theta _{\rm s}},{\phi _{\rm s}})\cos ({\theta _{\rm s}})\sin ({\theta _{\rm s}}){\rm{d}}{\theta _{\rm s}}{\rm{d}}{\varphi _{\rm s}}$$ (20)

      式中:θi${\phi _{\rm i}}$为入射光的入射角和方位角;θs${\varphi _{\rm s}}$为散射光的散射角和方位角。

      南京市受人类活动影响,市内风景区、郊区、居民区、道路、工业混合区和钢铁工业区的各类降尘粒子含量各不相同[12]。结合南京当地环境条件,根据不同亲疏水污染粒子占比在镜面分布情况,笔者进一步分析了多粒子分布的光学镜面散射特性随环境湿度的变化。对于不同类型的亲水性粒子和疏水性粒子,利用公式(19)和(20)可以分别计算单一粒径条件下的BRDF和TIS。真实环境下,不同区域光学镜面上亲水性粒子和疏水性粒子比例可能各有不同[9, 14],可以通过线性叠加亲水性粒子和疏水性粒子的BRDF和TIS求得整个镜面的BRDF和TIS,如图4所示。为便于讨论,将亲水性粒子占比表示为A,剩余部分为疏水性粒子占比。为便于计算且不失一般性,做如下简化:(1)所有污染粒子均匀分布;(2)不考虑相干效应和凝结效应[2]

      图  4  不同条件下BRDF随入射角度的变化 (λ=0.55 μm)

      Figure 4.  Variation of BRDF with incident angle under different conditions (λ=0.55 μm)

      表1~表3分别计算了不同亲疏水粒子占比、不同光学镜面洁净度等级和不同入射角度条件下相对应的TIS值。

      表 1  不同种类粒子占比时的TIS (λ=0.55 μm G1=300)

      Table 1.  TIS for surfaces with different proportions of hydrophilic particles (λ=0.55 μm G1=300)

      Humidity00.30.60.9
      TIS(A:30%)7.01×10−57.14×10−57.54×10−51.85×10−4
      TIS(A:40%)7.01×10−57.18×10−57.72×10−52.23×10−4
      TIS(A:50%)7.01×10−57.23×10−57.90×10−52.61×10−4
      TIS(A:60%)7.01×10−57.27×10−58.08×10−52.99×10−4
      TIS(A:70%)7.01×10−57.32×10−58.25×10−53.37×10−4
      TIS(A:80%)7.01×10−57.37×10−58.43×10−53.76×10−4

      表 2  不同镜面洁净度等级时的TIS (λ=0.55 μm)

      Table 2.  TIS for different mirror cleanliness levels (λ=0.55 μm)

      Mirror class300500700
      TIS(XRH=0)7.01×10−58.10×10−44.56×10−3
      TIS(XRH=0.3)7.27×10−58.63×10−44.73×10−3
      TIS(XRH=0.6)8.08×10−51.02×10−35.26×10−3
      TIS(XRH=0.9)2.99×10−45.23×10−31.95×10−2

      表 3  不同角度入射时的TIS

      Table 3.  TIS at different angles of incidence (λ=0.55 μm G1=300)

      Humidity00.30.60.9
      TIS(0°)7.02×10−57.62×10−58.97×10−53.71×10−4
      TIS(5°)7.02×10−57.61×10−58.97×10−53.71×10−4
      TIS(10°)7.02×10−57.61×10−58.96×10−53.71×10−4
      TIS(15°)7.01×10−57.61×10−58.96×10−53.71×10−4
      TIS(30°)6.99×10−57.59×10−58.93×10−53.70×10−4
      TIS(50°)6.95×10−57.54×10−58.89×10−53.70×10−4

      图4(a)可知,不同的光学镜面洁净度对散射造成的影响不一样,镜面洁净度越高其光学镜面散射损耗越大。由图4(b)可知,不同湿度影响的镜面散射差异较大,当相对湿度大于60%时造成散射损耗明显增大。

      分析表1~表3可知,环境湿度越高导致相应的光学镜面TIS值越大。高湿度条件下的TIS相比于干燥环境下的TIS有数倍的增加,从而会对光学系统的 性能造成不良影响。由表1可知,亲疏水性粒子占比不同致使光学镜面散射特性不同,在相同的湿度条件下,亲水性粒子占比越大,光学镜面损耗越大。需要指出的是,无论亲水性粒子占比如何,在环境湿度达到潮解相对湿度时,均会出现明显的散射增强效应,且亲水性粒子占比越大,散射增强越明显。分析表2可知,镜面洁净度等级越高,镜面散射越强烈。由表3可知,在同一环境相对湿度条件下,光学镜面的TIS值随着入射角度的增加而逐渐减小。其结果与参考文献[15]的结论相吻合,说明该模型适应于不同角度入射情形下的镜面散射特性分析。

      对南京市降尘粒子成分进行分析[1],南京市钢铁工业区降尘量远远高于其他地区,主要是因为工业污染区的降尘主要是由于工业燃煤飞灰引起的。因此,南京市钢铁工业区降尘成分中亲水性粒子占比大约为30%,疏水性粒子占比大约为70%。结合对不同环境条件下镜面散射特性分析,南京市钢铁工业区在环境相对湿度为60%时镜面散射的TIS值为7.54×10−5。由于不同地区的污染粒子种类、分布、吸湿特性均有所差异,进而会导致镜面的散射特性各不相同。为了更精确地分析不同地区不同的环境湿度对镜面散射特性的影响,需要对不同地区污染粒子的相关数据进行实地测量。

    • 文中在考虑环境湿度影响的情况下,对Mie理论进行了修正。以湿度较大的南京地区为例,分析了环境相对湿度对单个污染粒子以及多粒子条件下光学镜面散射特性的影响。结果表明,环境湿度的升高会导致污染粒子的粒径增加,进而使其折射率的实部与虚部均随着相对湿度的增加而减小。对于亲水性污染粒子而言,相对湿度的增加会导致其散射效率因子明显增大,散射损耗增强。在相同的湿度条件下,亲水性粒子占比越大,光学镜面散射越大。在不同的镜面洁净度等级条件下,散射损耗随着洁净度等级增大而明显增大。文中研究成果有助于对不同湿度环境中的光学系统性能进行分析与评估。

参考文献 (15)

目录

    /

    返回文章
    返回