留言板

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

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

近180°水中悬浮颗粒物体积散射函数测量

陈都 刘秉义 杨倩 唐军武 吴松华

陈都, 刘秉义, 杨倩, 唐军武, 吴松华. 近180°水中悬浮颗粒物体积散射函数测量[J]. 红外与激光工程, 2021, 50(6): 20211029. doi: 10.3788/IRLA20211029
引用本文: 陈都, 刘秉义, 杨倩, 唐军武, 吴松华. 近180°水中悬浮颗粒物体积散射函数测量[J]. 红外与激光工程, 2021, 50(6): 20211029. doi: 10.3788/IRLA20211029
Chen Du, Liu Bingyi, Yang Qian, Tang Junwu, Wu Songhua. Measurement of volume scattering function of suspended particles in water approaching to 180°[J]. Infrared and Laser Engineering, 2021, 50(6): 20211029. doi: 10.3788/IRLA20211029
Citation: Chen Du, Liu Bingyi, Yang Qian, Tang Junwu, Wu Songhua. Measurement of volume scattering function of suspended particles in water approaching to 180°[J]. Infrared and Laser Engineering, 2021, 50(6): 20211029. doi: 10.3788/IRLA20211029

近180°水中悬浮颗粒物体积散射函数测量

doi: 10.3788/IRLA20211029
基金项目: 国家重点研发计划(2016YFC1400905, 2016YFC1400904)
详细信息
    作者简介:

    陈都,男,硕士生,主要从事海洋激光探测与水体散射测量方面的研究

  • 中图分类号: O430

Measurement of volume scattering function of suspended particles in water approaching to 180°

  • 摘要:

    针对目前国内外现有的体积散射函数测量系统在后向小角度散射测量上的局限性,提出了基于离轴反射式光路的近180°水体体积散射函数测量方法并研发了实验室测量系统。系统采用离轴抛物面反射镜,将后向小角度散射光和入射激光分离,减小了系统后向小角度散射的测量盲区,而且能够获取全方位角的后向小角度散射光信号。选取聚苯乙烯标准粒子用于测量系统定标检验,结果表明,定标后的测量系统能够完成在173°~179.4°范围内水中悬浮颗粒物体积散射函数的测量,角度分辨率为0.01°。经对比分析,体积散射函数测量值与米散射理论值具有很好的一致性,验证了系统测量近180°水体体积散射函数的准确性和可行性。

  • 图  1  体积散射函数定义图[2]

    Figure  1.  Geometry for defining the volume scattering function[2]

    图  2  球形粒子米散射示意图

    Figure  2.  Schematic diagram of scattering by a sphere

    图  3  近180°水体体积散射函数测量系统原理图(a)和实物模型(b)

    Figure  3.  Schematic diagram (a) and physical model (b) of the system for measuring volume scattering function of water approaching to 180°

    图  4  离轴抛物面反射镜截面图

    Figure  4.  Sectional view of off-axis parabolic mirror

    图  5  测量系统CMOS相机像素坐标与颗粒物散射角度的对应关系

    Figure  5.  Correspondence between CMOS camera pixel coordinates and particle scattering angle of measurement system

    图  6  亚微米级(a)和微米级(b)聚苯乙烯标准粒子散射相函数最小值归一化曲线

    Figure  6.  Minimum normalized scattering phase function of Submicron (a) and Micron (b) polystyrene standard particles

    图  7  1 μm (a)和10 μm (b)粒径标准粒子体积散射函数测量值与理论值(方位角Φ=90°)

    Figure  7.  Measured value and theoretical value of volume scattering function of 1 μm (a) and 10 μm (b) size standard particle (azimuth Φ=90°)

    图  8  1 μm (a)和10 μm (b)粒径标准粒子体积散射函数测量值与理论值(方位角Φ=0°)

    Figure  8.  Measured value and theoretical value of volume scattering function of 1 μm (a) and 10 μm (b) size standard particle (azimuth Φ=0°)

    表  1  文中使用的聚苯乙烯标准粒子规格

    Table  1.   Specification of polystyrene standard particles used in this paper

    MaterialμD/μmPS σD/μmCV<5%nN0/mL
    Polystyrene0.09651%0.00191.97%1.5981.82×1013
    Polystyrene0.9771%0.02542.6%1.5981.82×1010
    Polystyrene10.081%0.242.38%1.5981.82×107
    下载: 导出CSV

    表  2  不同浓度聚苯乙烯标准粒子悬浮溶液散射系数

    Table  2.   Scattering coefficient of polystyrene standard particle suspension solution with different concentrations

    ConcentrationD=0.1 μmD=1 μmD=10 μm
    Dilute 1 μL0.0634/m1.67/m0.1366/m
    Dilute 5 μL0.3154/m8.345/m0.685/m
    下载: 导出CSV
  • [1] Seelye Martin. An Introduction to Ocean Remote Sensing[M]. UK: Cambridge University Press, 2013.
    [2] Mobley C D. Light and Water: Radiative Transfer in Natural Waters[M]. USA: Academic Press, 1994.
    [3] Kong Xiaojuan, Liu Bingyi, Yang Qian, et al. Simulation of water optical property measurement with shipborne lidar [J]. Infrared and Laser Engineering, 2020, 49(2): 0205010. (in Chinese) doi:  10.3788/IRLA202049.0205010
    [4] Berthon J F, Shybanov E, Lee E G, et al. Measurements and modeling of the volume scattering function in the coastal northern Adriatic Sea [J]. Applied Optics, 2007, 46(22): 5189-5203. doi:  10.1364/AO.46.005189
    [5] Chami M, Shybanov E B, Khomenko G A, et al. Spectral variation of the volume scattering function measured over the full range of scattering angles in a coastal environment [J]. Applied Optics, 2006, 45(15): 3605-3619. doi:  10.1364/AO.45.003605
    [6] Petzold T J. Volume scattering functions for selected ocean waters[R]. San Diego: Scripps Institute of Oceanography, 1972: 72–78.
    [7] Kullenberg G. Light-scattering observations in frontal zones [J]. Journal of Geophysical Research, 1978, 83(C9): 4683-4690. doi:  10.1029/JC083iC09p04683
    [8] Babin M, Stramski D, Reynolds R A, et al. Determination of the volume scattering function of aqueous particle suspensions with a laboratory multi-angle light scattering instrument [J]. Applied Optics, 2012, 51(17): 3853-3873. doi:  10.1364/AO.51.003853
    [9] Shao B, Jaffe J S, Chachisvilis M, et al. Angular resolved light scattering for discriminating among marine picoplankton: modeling and experimental measurements [J]. Optics Express, 2007, 14(25): 12473-12484.
    [10] Tan H, Doerffer R, Oishi T, et al. A new approach to measure the volume scattering function [J]. Optics Express, 2013, 21(16): 18697-18711. doi:  10.1364/OE.21.018697
    [11] Chami M, Thirouard A, Harmel T. POLVSM (Polarized Volume Scattering Meter) instrument: An innovative device to measure the directional and polarized scattering properties of hydrosols [J]. Optics Express, 2014, 22(21): 26403-26428. doi:  10.1364/OE.22.026403
    [12] Lee M E, Lewis M R. A new method for the measurement of the optical volume scattering function in the Upper Ocean [J]. Journal of Atmospheric & Oceanic Technology, 2003, 20(4): 563-571.
    [13] Sullivan J M, Twardowski M S. Angular shape of the oceanic particulate volume scattering function in the backward direction [J]. Applied Optics, 2009, 48(35): 6811-6819. doi:  10.1364/AO.48.006811
    [14] Slade W H, Agrawal Y C, Mikkelsen O A. Comparison of measured and theoretical scattering and polarization properties of narrow size range irregular sediment particles[C]//OCEANS – San Diego. IEEE, 2013: 1-6.
    [15] Li Cai, Cao Wenxi, Jing Yu, et al. An instrument for in situ measuring the volume scattering function of water: design, calibration and primary experiments [J]. Sensors, 2012, 12(4): 4514-4533. doi:  10.3390/s120404514
    [16] Wang W, Yang K, Li W, et al. Measuring the three-dimensional volume scattering functions of microsphere suspension: Design and laboratory experiments [J]. IEEE Photonics Journal, 2019, 11(5): 1-17.
    [17] Guo Yilu, Tao Bangyi, Huang Haiqing, et al. Wide angle volume scattering function measurement methods for particles in water [J]. Infrared and Laser Engineering, 2020, 49(2): 0203011. (in Chinese) doi:  10.3788/IRLA202049.0203011
    [18] Bohren C F, Huffman D R. Absorption and Scattering of Light by Small Particles[M]. Weinheim: Wiley, 1983: 83-129.
  • [1] 李路, 邢昆明, 赵明, 邓迁, 王邦新, 庄鹏, 施云.  探测气溶胶-水汽的拉曼-米散射激光雷达系统 . 红外与激光工程, 2023, 52(4): 20220484-1-20220484-11. doi: 10.3788/IRLA20220484
    [2] 熊志航, 麦浩基, 黄庄钒, 黎经腾, 孙培韬, 王嘉霖, 谢永韬, 何梓熙, 曾亚光, 王宏剑, 郭志明, 廖然, 马辉.  基于偏振光散射与荧光测量的水中悬浮颗粒现场快速分类仪 . 红外与激光工程, 2023, 52(9): 20230030-1-20230030-13. doi: 10.3788/IRLA20230030
    [3] 沈法华, 徐菁苑, 范安冬, 谢晨波, 王邦新, 杨亮亮, 周慧, 徐华.  基于双法布里-珀罗干涉仪多纵模米散射多普勒激光雷达技术 . 红外与激光工程, 2023, 52(7): 20220762-1-20220762-10. doi: 10.3788/IRLA20220762
    [4] 陶雨婷, 赵泓楷, 周雨迪, 卓文淇, 刘群, 崔晓宇, 刘斌, 乐成峰, 潘孙强, 刘崇, 刘东.  CALIOP反演海洋颗粒物后向散射系数方法概述 . 红外与激光工程, 2021, 50(6): 20211037-1-20211037-9. doi: 10.3788/IRLA20211037
    [5] 陈慧敏, 马超, 齐斌, 郭鹏宇, 杨尚贤, 高丽娟, 霍健.  脉冲激光引信烟雾后向散射特性研究 . 红外与激光工程, 2020, 49(4): 0403005-0403005-7. doi: 10.3788/IRLA202049.0403005
    [6] 郭乙陆, 陶邦一, 黄海清, 吴超钒, 宋宏, 冷建兴.  水中颗粒物大角度范围体散射函数测量方法 . 红外与激光工程, 2020, 49(2): 0203011-0203011. doi: 10.3788/IRLA202049.0203011
    [7] 陈鹏, 赵继广, 杜小平, 宋一铄.  基于粒子群优化的近似散射相函数拟合方法 . 红外与激光工程, 2019, 48(12): 1203005-1203005(7). doi: 10.3788/IRLA201948.1203005
    [8] 滕曼, 庄鹏, 张站业, 李路, 姚雅伟.  大气气溶胶污染监测中应用的新型全天时户外型拉曼-米散射激光雷达系统 . 红外与激光工程, 2019, 48(7): 706001-0706001(7). doi: 10.3788/IRLA201948.0706001
    [9] 单良, 徐良, 洪波, 曹丽霞, 王道档, 郭天太, 孔明.  小角前向散射偏振比法颗粒粒度分布反演 . 红外与激光工程, 2019, 48(1): 117001-0117001(8). doi: 10.3788/IRLA201948.0117001
    [10] 王雪敏, 申晋, 徐敏, 黄钰, 高明亮, 刘伟, 王雅静.  多角度动态光散射角度误差对权重估计的影响 . 红外与激光工程, 2018, 47(10): 1017004-1017004(7). doi: 10.3788/IRLA201847.1017004
    [11] 王凤杰, 陈慧敏, 马超, 龙胤宇.  云雾后向散射激光回波特性研究 . 红外与激光工程, 2018, 47(5): 506002-0506002(5). doi: 10.3788/IRLA201847.0506002
    [12] 樊长坤, 李琦, 赵永蓬, 陈德应.  不同粗糙度抛物面的2.52 THz后向散射测量 . 红外与激光工程, 2017, 46(11): 1125005-1125005(6). doi: 10.3788/IRLA201746.1125005
    [13] 石振华, 林冠宇, 王淑荣, 于磊, 曹佃生.  基于激光散射理论的微小颗粒测量的数值分析 . 红外与激光工程, 2015, 44(7): 2189-2194.
    [14] 张学海, 魏合理.  基于用几何光学和米散射法的球形粒子前向散射特性计算研究 . 红外与激光工程, 2015, 44(5): 1485-1490.
    [15] 陶会荣, 张福民, 曲兴华.  无合作目标测量中目标表面后向散射特性的实验研究 . 红外与激光工程, 2014, 43(S1): 95-100.
    [16] 王安祥, 吴振森.  光散射模型中遮蔽函数的参数反演 . 红外与激光工程, 2014, 43(1): 332-337.
    [17] 薛向尧, 高云国.  利用光束后向散射提高激光跟踪发射精度方法 . 红外与激光工程, 2013, 42(8): 2003-2007.
    [18] 万俊, 张晓晖, 饶炯辉, 梁瑞涛.  基于独立成分分析的舰船气泡尾流后向散射光信号处理 . 红外与激光工程, 2013, 42(1): 244-250.
    [19] 崔骥, 付欢, 刘烨, 王清华, 李振华, 贺安之.  典型粗糙表面激光后向散射特性实验 . 红外与激光工程, 2013, 42(5): 1208-1211.
    [20] 李琛, 舒志峰, 夏海云, 孙东松, 周广超, 汪建业.  利用Fabry-Perot标准具的米散射测风激光雷达接收机 . 红外与激光工程, 2012, 41(1): 107-111.
  • 加载中
计量
  • 文章访问数:  487
  • HTML全文浏览量:  162
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-05-06
  • 修回日期:  2021-05-21
  • 刊出日期:  2021-06-30

近180°水中悬浮颗粒物体积散射函数测量

doi: 10.3788/IRLA20211029
    作者简介:

    陈都,男,硕士生,主要从事海洋激光探测与水体散射测量方面的研究

基金项目:  国家重点研发计划(2016YFC1400905, 2016YFC1400904)
  • 中图分类号: O430

摘要: 

针对目前国内外现有的体积散射函数测量系统在后向小角度散射测量上的局限性,提出了基于离轴反射式光路的近180°水体体积散射函数测量方法并研发了实验室测量系统。系统采用离轴抛物面反射镜,将后向小角度散射光和入射激光分离,减小了系统后向小角度散射的测量盲区,而且能够获取全方位角的后向小角度散射光信号。选取聚苯乙烯标准粒子用于测量系统定标检验,结果表明,定标后的测量系统能够完成在173°~179.4°范围内水中悬浮颗粒物体积散射函数的测量,角度分辨率为0.01°。经对比分析,体积散射函数测量值与米散射理论值具有很好的一致性,验证了系统测量近180°水体体积散射函数的准确性和可行性。

English Abstract

    • 海洋水体中的环境参数主要包括温、盐度剖面信息、深度分布以及海水组分信息等,近些年来,海洋水色遥感技术一直是监测海洋环境参数重要的海洋探测手段,能够提供水体光学参数、叶绿素浓度、浅海水深等海洋水体的特征参数,这些观测数据能作为在海洋环境监测、气候变化、灾害预警等领域的重要数据支持[1]。近些年来,海洋激光雷达作为新型的主动式海洋水色遥感技术,在剖面探测、夜间观测等方面具有显著优势,而且不受大气校正的影响。无论主动式或被动式海洋水色遥感的研究,都离不开海水固有光学性质的精确测量。

      体积散射函数是表征海洋水体颗粒物散射特性的重要固有光学参数[2],是海洋光学最基础的物理量之一,通过体积散射函数和吸收系数可以得到海水的所有固有光学性质。后向小角度体积散射函数的准确测量,对于海洋辐射传输模型的建立、海洋激光雷达原理印证和水中悬浮颗粒物特征信息的获取具有重要意义[2-5]。体积散射函数在后向小角度范围内存在一定的起伏特征,目前国内外现有的体积散射函数测量系统几乎都是针对广角和前向小角度范围内的散射测量系统[6-17];而且受杂散光的干扰,能覆盖到后向小角度范围的散射测量系统的测量值普遍出现明显上升趋势[17]

      水体后向小角度散射测量的难点在于:(1)海洋水体中大多数悬浮颗粒物的近180°体积散射函数的大小通常比0°附近体积散射函数小六七个数量级,信号十分微弱;(2)在测量系统设计上,激光器和探测器都具有一定的体积,通常由于自身元件遮挡导致后向小角度散射光信号难以测量;(3)激光在入射到样品池侧壁上时,会产生高能量反射光,对后向小角度散射产生很大的噪声干扰。现有测量手段难以满足后向小角度范围内体积散射函数的测量需求。

      文中针对以上问题,提出了一种基于离轴反射式光路的近180°水体体积散射函数测量系统,散射角测量范围为173°~179.4°,角度分辨率达到0.01°。通过测量聚苯乙烯标准粒子悬浮溶液,获取的体积散射函数实验测量值与米散射理论值具有较好的一致性。

    • 体积散射函数就是描述散射光在空间上的角分布情况,定义示意图如图1所示[2]$ {\varPhi }_{i}(\lambda ) $为入射到$ \varDelta{A} $面积上的光功率,$ \varDelta {\varPhi }_{a}(\lambda ) $是在$ \varDelta{r} $距离内吸收的光功率,$ {{\varDelta }^{2}\varPhi }_{s}\left(\varPsi ,\lambda \right) $是散射到$ \varDelta \varOmega $立体角范围内圆环上的光功率,$ \varDelta {\varPhi }_{t}(\lambda ) $是透射光功率。假设水中粒子各向同性且忽略粒子的偏振特性,体积散射函数定义为下面公式,见公式(1):

      $$ \beta \left( {\varPsi ,\lambda } \right) = \mathop {\lim }\limits_{\varDelta r \to 0} \mathop {\lim }\limits_{\varDelta \varOmega \to 0} \frac{{{\varDelta ^2}{\varPhi _s}\left( {\varPsi ,\lambda } \right)}}{{{\varPhi _i}\left( \lambda \right)\varDelta r\varDelta \varOmega }} $$ (1)

      其物理意义表示为单位距离、单位立体角的散射率。光谱辐射强度的定义为单位立体角上的散射功率,因此体积散射函数$ \beta \left(\varPsi ,\lambda \right) $又可以写成公式(2):

      $$ \beta \left(\varPsi ,\lambda \right)=\mathop {{\rm{lim}}}\limits_{\Delta V \to 0} \frac{\varDelta {I}_{s}\left(\varPsi ,\lambda \right)}{{E}_{i}(\lambda )\varDelta V} $$ (2)

      式中:$ \varDelta {I}_{s}\left(\varPsi ,\lambda \right) $$ \Delta\varOmega $立体角内的相对散射强度;$ {E}_{i}(\lambda ) $是入射辐照度;$ \Delta{V} $是散射体体积。公式(2)物理意义表示为单位体积内、单位入射辐照度的散射强度。因此,可以通过测量不同散射角度的光强来获取水体的体积散射函数。

      图  1  体积散射函数定义图[2]

      Figure 1.  Geometry for defining the volume scattering function[2]

    • 海水中光学活性成分包括纯海水、颗粒物和可溶有机物,粒径大小集中于0.1 μm到几十微米不等,散射分布符合米散射基本特征。如图2所示为球形粒子米散射示意图[18],一束波长为$ \lambda $、强度为$ {I}_{0} $的线偏振光入射到O点处一直径为d的均匀介质球形粒子上,光束振幅为E,传播方向沿z轴方向,其电场、磁场的振动方向分别沿x轴和y轴方向;OP为散射光,散射角$ \theta $z轴与散射光之间的夹角,方位角$ \varphi $为散射光再xOy平面上的投影与x轴的夹角,OPP'为散射平面。根据米散射理论,P点处的散射光强度可以表示为:

      $$ {I}_{s}=\frac{{\lambda }^{2}{I}_{0}}{4{\pi }^{2}{r}^{2}}[{\left|{S}_{1}\left(\theta \right)\right|}^{2}{\rm {sin}}^{2}\varphi +{\left|{S}_{2}\left(\theta \right)\right|}^{2}{\rm {cos}}^{2}\varphi ] $$ (3)

      式中:$ r $P点到散射体中心的距离;$ {\left|{S}_{1}\left(\theta \right)\right|}^{2} $$ {\left|{S}_{2}\left(\theta \right)\right|}^{2} $分别表示为平行于和垂直于散射平面的强度分量;$ {S}_{1}\left(\theta \right) $$ {S}_{2}\left(\theta \right) $为散射振幅函数,均可通过相应公式获得。

      根据球形粒子米散射理论,可以在已知入射光偏振方向的条件下,得到球形粒子在空间内任意位置的散射强度,进而获取不同方位角散射平面上的散射相函数及体积散射函数,为利用聚苯乙烯标准粒子对测量系统进行定标检验提供参考依据。

      图  2  球形粒子米散射示意图

      Figure 2.  Schematic diagram of scattering by a sphere

    • 文中研制的近180度水体体积散射函数测量系统的基本原理图(a)和实物模型(b)如图3所示,主要包括激光发射器、离轴抛物面反射镜、散射样品、非球面透镜组、孔径光阑、CMOS偏振敏感相机以及计算机,散射样品位于离轴抛物面反射镜的焦点处。如图4为离轴抛物面反射镜截面图,离轴抛物面反射镜的反射面是母抛物面的截面,此截面不在母抛物面光轴中心。离轴抛物面反射镜的焦点与母抛物面的焦点重合,光轴1为母抛物面的光轴,光轴2是离轴抛物面反射镜的光轴,两光轴相互平行且互相偏离。图中母焦距即为母抛物面的焦距,反射有效焦距为离轴抛物面反射镜的中心点与母抛物面焦点之间的距离,反射有效焦距与离轴抛物面反射镜光轴之间的夹角θ为离轴角,文中选取的离轴抛物面反射镜离轴角为90°。离轴抛物面反射镜背面开有锥形结构的通孔,该通孔穿过反射表面而形成小孔,通孔的中轴线与母抛物面的焦轴重合。入射激光经过离轴抛物面镜背面小孔,穿过反射面,入射到散射样品上,根据抛物面反射镜光学特性,散射样品的近180°后向散射光经离轴抛物面反射镜反射变换为平行光束,再经过非球面透镜组与孔径光阑组成的双侧远心光路准直且调整束宽后由CMOS偏振敏感相机接收散射光信号。根据后向散射角与CMOS像素点位置的对应关系,即可获得水体的各个方位角散射平面的近180°体积散射函数。

      图  3  近180°水体体积散射函数测量系统原理图(a)和实物模型(b)

      Figure 3.  Schematic diagram (a) and physical model (b) of the system for measuring volume scattering function of water approaching to 180°

      图  4  离轴抛物面反射镜截面图

      Figure 4.  Sectional view of off-axis parabolic mirror

      综上所述,根据测量系统几何关系,CMOS相机上成像的不同像素点位置与后向散射角具有一一对应的关系,在定标检验过程中需要根据像素点位置对应的散射角度完成不同角度与光强关系的反演。以抛物面顶点为原点建立笛卡尔三维坐标系,假设在CMOS相机上的任一像素点横纵坐标分别为XY,根据几何模型关系,该像素点坐标XY与水体后向小角度散射角$ \varPsi $的标准方程函数关系见公式(4)。式中:$ {L}_{1} $是透镜1的焦距;$ {L}_{2} $是透镜2的焦距;$ {L}_{PFL} $是离轴抛物面反射镜母焦距的长度;$ {L}_{RFL} $是离轴抛物面镜反射有效焦距的长度,单位均为mm;$ {n}_{w} $是水体折射率。

      在文中测量系统中,颗粒物散射光能够覆盖的相机成像区域是以中心像素为圆心,直径为6 mm的像区,大约覆盖230万个像素点,在体积散射函数反演过程中,需要根据在成像区域内的所有像素点对应的散射角来计算不同角度的散射强度。通过光学仿真软件光线追迹来确定测量系统中的颗粒物散射角度和像素坐标的关系。图5是利用光学仿真软件光线追迹获取的颗粒物散射角度和像素坐标的对应关系,水平中线测量角度范围是172.5°~179.4°,垂直中线测量角度范围是173°~179.4°。由于离轴抛物面反射镜的非球面特点,系统的标准散射角度关于横轴对称,关于纵轴不对称,在172.5°~173.5°范围内,并不能获取在全方位角范围内的散射强度信息。

      图  5  测量系统CMOS相机像素坐标与颗粒物散射角度的对应关系

      Figure 5.  Correspondence between CMOS camera pixel coordinates and particle scattering angle of measurement system

      $$ \varPsi ={\rm {arcsin}}\left({\rm {sin}}\left(180-{\rm {arctan}}\dfrac{\sqrt{{\left(\dfrac{{L}_{1}}{{L}_{2}}Y\right)}^{2}+{\left(\dfrac{1}{4{L}_{PFL}}\left[{\left(-{L}_{RFL}-\dfrac{{L}_{1}}{{L}_{2}}X\right)}^{2}+{\left(\dfrac{{L}_{1}}{{L}_{2}}Y\right)}^{2}\right]-{L}_{PFL}\right)}^{2}}}{\left|-{L}_{RFL}-\dfrac{{L}_{1}}{{L}_{2}}X\right|}\right)/{n}_{w}\right) $$ (4)
    • 在选取样品溶液时,需要考虑海洋水体中颗粒物粒径的真实分布范围,同时选取在系统的散射角度测量范围内颗粒物散射光强度变化趋势具有代表性的聚苯乙烯标准粒子作为测量样品。选取要求主要依据以下两点,一是选取在系统的散射角度测量范围内的散射光强度变化趋势平缓的标准粒子用于测量系统的强度定标;二是选取在系统的散射角度测量范围内的散射光强度变化具有波动起伏趋势的标准粒子作为定标后系统的验证。图6展示了亚微米级(a)、微米级(b)聚苯乙烯标准粒子在测量范围内的散射相函数的最小值归一化曲线,通过利用各自在角度测量范围内的理论散射相函数的最小值做完归一化之后,可以看出不同粒径标准粒子的散射相函数变化趋势的差异。如图6(a)所示,在散射角度测量范围内,亚微米级(a)聚苯乙烯标准粒子的散射相函数呈单调变化趋势,随着粒径的增大,散射相函数最大值与最小值的差异逐渐增大;粒径小于0.5 μm的标准颗粒的散射相函数的幅值变化很小,1 μm粒径标准颗粒的散射相函数的幅值变化最大,最大值与最小值相差约5倍。如图6(b)所示,微米级聚苯乙烯标准粒子在测量角度范围内的散射相函数具有典型的波动起伏特征,其波动的幅值和频率均随粒径增大。

      图  6  亚微米级(a)和微米级(b)聚苯乙烯标准粒子散射相函数最小值归一化曲线

      Figure 6.  Minimum normalized scattering phase function of Submicron (a) and Micron (b) polystyrene standard particles

      综上所述,根据亚微米级和微米级聚苯乙烯标准粒子在角度测量范围内的散射相函数变化特征,文中选用0.1 μm粒径的聚苯乙烯标准粒子用于系统定标,选用1 μm和10 μm粒径的聚苯乙烯标准粒子用于定标后系统测量结果的验证。表1是文中选取的聚苯乙烯标准粒子的规格,μD为实际平均直径,PS为固含量(固相含量1% g/mL),${\sigma }_{\rm D}$为粒径标准偏差,CV为变异系数(coefficient of variation),定义为粒径标准偏差比上实际平均直径,n为标准粒子折射率,$ {N}_{0} $为每mL样品中的微球个数。实际中,标准粒子悬浮溶液中的粒径分布遵循对数正态分布,依据实际平均直径、粒径标准偏差、折射率以及每mL样品中的微球个数便可以根据球形粒子的米散射特征计算出理论的体积散射函数。

      在目前国内外现有的体积散射函数测量系统中,无论是反演体积散射函数还是根据米散射理论计算标准颗粒的散射特性,都有一个基本的假设,就是激光在样品溶液中传输时只考虑单次散射,通常以散射样品的光学厚度τ作为参考依据[8]。文中稀释溶液去离子水的溶液量为25 mL,表2是三种标准粒子在抽取不同样品溶液量稀释后的情况下,根据米散射理论计算出的样品溶液理论散射系数,以不同海洋水体光学性质的数量级作为参考,0.1、1、10 μm粒径聚苯乙烯标准粒子在分别抽取1~5 μL的溶液量时,理论散射系数数量级能与不同的海洋水体光学性质相匹配,且样品溶液的光学厚度τ均小于0.1,则可以认为样品溶液中不会产生多次散射。在标准粒子的样品溶液配制完成后,使标准颗粒微球能够在去离子水中分散均匀。

      表 1  文中使用的聚苯乙烯标准粒子规格

      Table 1.  Specification of polystyrene standard particles used in this paper

      MaterialμD/μmPS σD/μmCV<5%nN0/mL
      Polystyrene0.09651%0.00191.97%1.5981.82×1013
      Polystyrene0.9771%0.02542.6%1.5981.82×1010
      Polystyrene10.081%0.242.38%1.5981.82×107

      表 2  不同浓度聚苯乙烯标准粒子悬浮溶液散射系数

      Table 2.  Scattering coefficient of polystyrene standard particle suspension solution with different concentrations

      ConcentrationD=0.1 μmD=1 μmD=10 μm
      Dilute 1 μL0.0634/m1.67/m0.1366/m
      Dilute 5 μL0.3154/m8.345/m0.685/m
    • 根据米散射理论,当线偏振光入射到已知浓度和粒径分布的球形粒子上时,可以计算出在不同方位角散射平面上的体积散射函数。测量过程中选取0.1 μm标准粒子作为定标粒子,根据在特定方位角散射平面上测量的颗粒物散射强度$ {I}_{\rm{cal}}(\varPsi ,\varPhi ) $和相应的颗粒物体积散射函数$ {\beta }_{\rm{cal}}(\varPsi ,\varPhi ) $的理论值,获取定标系数R

      $$ R = \frac{{I}_{\rm{cal}}(\varPsi ,\varPhi )}{{\beta }_{\rm{cal}}(\varPsi ,\varPhi )} $$ (5)

      式中:$ {I}_{\rm{cal}}(\varPsi ,\varPhi ) $是定标粒子在方位角为$ {\varPhi} $的散射平面上的散射强度;$ {\beta }_{\rm{cal}}(\varPsi ,\varPhi ) $是定标粒子的理论体积散射函数,定义为颗粒物散射相函数与其散射系数的乘积。

      选取1 μm和10 μm标准粒子作为检验粒子,测量其在特定方位角散射平面上的散射强度,利用定标系数R获取其体积散射函数的数值:

      $$ {\beta }_{\rm{val}}(\varPsi ,\varPhi )=\frac{{I}_{\rm{val}}(\varPsi ,\varPhi )}{R}=\frac{{I}_{\rm{val}}(\varPsi ,\varPhi )}{{I}_{\rm{cal}}(\varPsi ,\varPhi )}\times {\beta }_{\rm{cal}}(\varPsi ,\varPhi ) $$ (6)

      式中:$ {I}_{\rm{val}}(\varPsi ,\varPhi ) $是检验粒子在方位角$ \varPhi $散射平面上的散射强度;$ {\rm{\beta }}_{\rm{val}}(\varPsi ,\varPhi ) $是检验粒子在方位角$ \varPhi $散射平面上的体积散射函数测量值。

    • 当入射激光偏振方向为水平方向时,离轴抛物面反射镜的长轴中线和短轴中线所在的散射平面的方位角分别为0°和90°,其散射光强度分别对应于CMOS相机水平中线和垂直中线的灰度值。由于散射平面的不同,在定标过程中,CMOS相机不同方位角上散射强度的定标系数不同。在文中定标检验过程中,选用相机垂直中线所在的90°方位角散射平面作为参考散射平面。相机垂直中线对应离轴抛物面反射镜的短轴中线,由于离轴抛物面反射镜的非轴对称特点,限制了在该散射平面的散射角度测量范围在173°~179.4°。图7(a)(b)的蓝线分别为1 μm和10 μm粒径标准粒子在水平线偏光入射时,利用颗粒物散射强度测量值和定标系数R得到的90°方位角散射平面上的体积散射函数测量值曲线,角度分辨率为0.01°;红线是根据米散射计算的两种标准粒子体积散射函数理论值曲线。

      图  7  1 μm (a)和10 μm (b)粒径标准粒子体积散射函数测量值与理论值(方位角Φ=90°)

      Figure 7.  Measured value and theoretical value of volume scattering function of 1 μm (a) and 10 μm (b) size standard particle (azimuth Φ=90°)

      图7可见,水平线偏振光入射条件下,在173°~177.7°散射角范围内,1 μm粒径标准粒子的体积散射函数的测量值和理论值的数值和变化趋势基本吻合。在177.7°之后,由于1 μm粒径标准粒子的理论散射相函数大小比其他两种标准粒子小一到两个数量级,因此在同样的实验条件下,1 μm标准粒子在系统测量角度范围内的散射强度比其他两种标准粒子小约1个数量级。当散射角大于177.7°时,由于背景杂散光较强,导致体积散射函数测量值信噪比较低。10 μm粒径标准粒子的体积散射函数测量值与理论值曲线在波动起伏和波动周期上具有较好的一致性,且在趋近180°附近的测量值具有较高的信噪比。两种标准粒子在体积散射函数的测量数值与理论数值上具有一些差异,在90°方位角散射平面上,1 μm和10 μm粒径标准粒子体积散射函数测量值与理论值在可信区间内的相对平均偏差为0.34和0.42。造成该误差的主要原因一是两次溶液抽取过程中存在浓度误差;二是悬浮溶液中存在粒子聚集的现象。

    • 当入射激光偏振方向为垂直方向时,离轴抛物面反射镜的长轴中线和短轴中线所在的散射平面的方位角分别为90°和0°,其散射光强度分别对应于CMOS相机水平中线和垂直中线的灰度值。与水平线偏光入射时相同,选取相机垂直中线对应的0°方位角散射平面作为参考散射平面。图8(a)图8(b)的蓝线分别为1 μm和10 μm粒径标准粒子在垂直线偏光入射时,利用颗粒物散射强度测量值和定标系数R得到的0°方位角散射平面上的体积散射函数测量值曲线,角度分辨率为0.01°;红线是根据米散射计算的两种标准粒子体积散射函数理论值曲线。

      图  8  1 μm (a)和10 μm (b)粒径标准粒子体积散射函数测量值与理论值(方位角Φ=0°)

      Figure 8.  Measured value and theoretical value of volume scattering function of 1 μm (a) and 10 μm (b) size standard particle (azimuth Φ=0°)

      图8可见,垂直线偏振光入射条件下,在173°~177.7°散射角范围内,1 μm粒径标准粒子的体积散射函数的测量值的数值和趋势都与理论值基本一致,在177.7°由于信噪比较低,测量值不具有参考意义;10 μm粒径标准粒子的体积散射函数测量值与理论曲线的数值和变化趋势都具有很好的一致性,在趋近180°附近的测量值仍具有较高的信噪比和可信度。在0°方位角散射平面上,1 μm和10 μm粒径标准粒子体积散射函数测量值与理论值在可信区间内的相对平均偏差为0.45和0.35。

    • 文中提出了一种基于离轴光路的近180度水体体积散射函数测量系统。测量系统采用离轴抛物面反射镜和双侧远心光路相结合的方式,有效解决了光路遮挡和各种杂散光的干扰问题,实现了对水中悬浮颗粒物的近180°体积散射函数的测量,散射角测量范围在173°~179.4°,角度分辨率达到0.01°。

      选取在测量角度范围内体积散射函数具有代表性变化趋势的三种粒径尺度的标准粒子完成测量系统定标和检验。通过改变入射激光的线偏方向,实现了1 μm和10 μm粒径标准粒子在90°和0°散射平面上体积散射函数的测量,测量结果与米散射理论值具有较好的一致性,验证了文中测量系统设计原理的可行性和测量结果的准确性。

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回