复杂海洋环境下天空光偏振特性模拟研究 下载: 1081次
1 引言
雾是由悬浮在近地面空气中缓慢沉降的水滴或冰晶质点组成的一种气溶胶系统,海雾是在特定的海洋环境和天气条件下形成并维持的气溶胶[1]。大气海雾介质对天空光的影响主要体现在微粒的散射作用上,微粒与非偏振的太阳光相互作用产生偏振光,使天空光具有稳定的偏振分布模式[2-3]。对大气海洋环境下天空偏振分布规律的研究在天文遥感、环境监测、激光通信等领域具有重要的意义[4-5]。
浑浊介质中光传播的研究以求解辐射传输方程为中心,问题的本质是小粒子的光散射问题[6-7]。在雾天辐射传输的研究中,关于矢量辐射传输方程的解法主要有球谐函数法[8]、倍加累加法[9]、逐次散射法[10]以及蒙特卡罗方法[11-12],基于以上研究方法,许多研究者做出了相应的研究。文献[ 13]针对雾天环境,采用蒙特卡罗方法研究辐射雾和平流雾中圆偏光的保偏性能,为雾天环境不同波段下偏振探测中偏振器件的选择做出理论指导;文献[ 14]采用蒙特卡罗方法重点研究波段及探测距离变化下雾天环境中偏振传输对探测范围的影响,为偏振探测中相应探测模式的选择提供了依据。但以上方法都局限在对水平方向雾天环境中的偏振传输,对于垂直方向经过大气、海雾多层海洋环境的偏振分布的情况并未见报道。针对垂直方向的多层介质,Evans等[15]提出的基于矢量辐射传输方程的RT3方法是计算非均匀多层介质传输特性时使用最广泛的方法之一,该方法基于直观物理过程采用倍加累加方法求解两层间的辐射和透射性质,进而求得整个全天空偏振分布模式。文献[ 16]采用该方法对简化的双层大气光波的Stokes参量进行求解,分别模拟晴天、多云及阴天情况下全天空的偏振信息。文献[ 17]采用基于倍加累加法的辐射传输模型研究可见光波段至近红外波段下多种天气条件的天空偏振模型,并由实测验证,仿真结果在80%的区域内与实测结果一致。但以上方法都局限于模拟陆地环境下天空光的偏振分布模式。
基于此,在现有研究的基础上,为了进一步研究复杂海洋环境下海面对天空偏振探测及通信时天空光偏振特性对测试结果的影响,采用基于矢量辐射传输方程的RT3方法仿真从大气顶层到海面下行辐射的传输过程,将目前研究中简化的双层大气结构扩展为更为复杂的三层结构,分别以晴朗天气下大气-海雾双层结构以及阴天环境下大气-云层-海雾三层结构作为研究对象,根据云层、海雾粒子尺寸分布特点,由Mie散射方法计算该层粒子的分布特性,以变化规律较明显的太阳子午线位置为基准,研究在不同云层、海雾浓度下观测高度角随天空下行辐射偏振度(DOP)的变化,掌握复杂海洋环境下天空光偏振分布模式的变化规律,为海洋环境下偏振探测的研究提供了参考。
2 基于倍加累加方法的天空偏振建模
2.1 基于RT3的大气偏振建模
辐射传输方程可描述为单色辐射量与大气的相互作用过程,包含了大气与电磁辐射的吸收、散射和热辐射等,解决大气介质中的光波或辐射传播问题就是求解满足边界条件的辐射传输方程,各种辐射传输问题可采用不同的处理方法。
RT3模型采用的是倍加累加法求解矢量辐射传输方程,两层间相互作用的基本原理是:若两个介质层的反射和透射矩阵已知,则组合层的辐射和透射矩阵可通过两层间的连续反射过程得到。层间作用过程如
图 1. 介质层间相互作用示意图
Fig. 1. Schematic diagram of the interaction principle between the medium layers
基于RT3的大气偏振建模的基本流程如
2.2 单次散射传输过程
由相位矩阵F到散射Mueller矩阵M的偏振传输过程可表示为
式中:i1和i2分别表示入射光线所在散射平面与子午面的夹角和出射光线所在散射平面与子午面的夹角;Θ为散射角;L为偏振旋转矩阵,其中
图 2. 基于RT3的大气偏振建模的基本流程
Fig. 2. Basic procedure of atmospheric polarization mode based on RT3
对于旋转对称粒子,散射相位矩阵中16个元素中只有6个独立元素,表示为
式中,若粒子为球形粒子,则可应用Mie散射求解散射相位矩阵中各元素a1~a6。a1~a6分别由勒让德多项式及广义勒让德多项式展开,表示为
其中,β、δ、γ和ε分别为散射相位矩阵中各元素的勒让德多项式的系数,l表示勒让德多项式的级数,Pl(cos Θ)和
当采用基于RT3的偏振建模方法时,需将上述Mie散射方法求得的散射相位矩阵中各元素全部展开为勒让德多项式,得到
其中,j=1,2,…,6表示散射相位矩阵中各元素,α为勒让德多项式系数。
3 大气-海雾复杂环境的多层粒子分布特性
3.1 大气-海雾复杂环境的分层
采用按大气温度场的垂直结构进行分层,将大气从地面起依次分为对流层、平流层、中层、热层以及外大气层。基本上所有的天气现象都发生在对流层内,其含有地球大气总质量的80%,以及几乎所有的水汽,是产生偏振现象的主要层。在大气偏振模式的研究中,根据散射粒子的不同将对流层进行划分,主要分为Rayleigh散射层和Mie散射层。其中,5~15 km的大气层水汽含量较少,气体粒子半径小,主要由大气分子组成,可用Rayleigh散射进行模拟;在0~5 km高度,除作为主体的气体成分之外,大气中还存在着各种液态和固态,以及液态和固态混合的粒子,其中最常见的是液态水构成的云和雾粒子,雾和各种高低层水云都是接近球形的水滴,当地表环境为海面时,海面平流运动的空气与海面之间的显热和潜热的交换作用,形成平流雾,该雾粒子直径较大,散射过程复杂,主要由Mie散射进行计算。
根据实际海面上海雾天气的情况对大气-海雾环境进行划分时需考虑的情况极其复杂,为了在简化计算的同时,突出海洋环境作用下最下层海雾层与中间云层对大气偏振模式作用的区别,研究晴朗天气无云条件下和有云条件下的海雾环境大气偏振模式,参考美国标准大气分别将大气简化为两层和三层的模式进行计算,如
图 3. 大气-海雾环境分层情况。 (a)晴朗天气无云情况;(b)有云情况
Fig. 3. Stratification of the atmosphere and advection fog environment. (a) Clear sky without cloud; (b) cloudy sky
3.2 海雾和水云层粒子分布特性
海雾和水云的外观和几何特征在空间和时间上的复杂可变性,使得海雾和水云层中球形水滴的尺寸大小不一。实际应用中,海雾和水云在单位体积内的粒径谱分布n
其中,N为粒子每立方厘米中总的数量密度,a、b、α'、γ'均为正实数,a是归一化常数,用于确保尺寸分布的积分在所有粒径产生的数量密度N之内,α'和γ'用于描述粒子尺寸分布的斜率,rmod为模半径,b=α'/rγ'mod·γ',A=Na。
修正的Gamma分布适用于描述具有单峰分布的粒子谱分布,
表 1. 云层和海雾层的修正Gamma分布参数
Table 1. Parameters of modified Gamma distribution for cloud and advection fog layer
|
图 4. 云层和海雾层天气情况下归一化粒子尺寸分布情况。 (a)云层; (b)海雾层
Fig. 4. Normalized particle size distribution for cloud layer and advection fog layer. (a) Cloud; (b) advection fog
根据
4 大气-海雾偏振仿真结果与分析
4.1 大气环境仿真结果及验证
由于目前未见对大气海雾分层环境偏振模式的研究,基于经典的Rayleigh散射理论,针对晴朗天空,假设大气粒子是均匀分布的大气分子,将整个大气层看作为一层进行处理,满足Rayleigh散射条件。将基于RT3的仿真建模方法得到的大气散射天空偏振分布与基于Rayleigh散射的大气理论模型[4]得到的分布结果进行对比,如
从
图 5. 大气散射天空偏振分布。 (a)基于RT3的仿真方法;(b)基于Rayleigh散射的大气理论模型
Fig. 5. Sky polarization distribution of atmosphere scattering. (a) Simulation method based on RT3; (b) atmospheric theory model based on Rayleigh scattering
图 6. 子午线上偏振度与观测高度角的关系
Fig. 6. Relationship between DOP and observed altitude angle on the meridian of the sun
4.2 大气-海雾环境仿真结果
4.2.1 晴朗无云天气下天空偏振特性的研究
在晴朗无云天气条件下,根据3.1节中对大气的分层情况,将其分为大气-海雾两层进行研究,如
图 7. 浓海雾和中海雾环境下偏振度与观测高度角的变化关系
Fig. 7. Relationship between DOP and observed altitude angle under heavy and moderate advection fog
4.2.2 有云影响下天空偏振特性研究
在对大气海雾环境天空偏振特性的研究中,海洋复杂环境多伴有云层的影响,如
为了进一步研究中海雾和浓海雾对偏振特性的影响,在相同云层下对比中海雾和浓海雾所得的结果,如
图 8. 海洋性积云和海洋性层云影响下海雾环境偏振度随观测高度角的变化。(a)海洋性积云和海洋性层云影响下浓海雾;(b)海洋性积云和海洋性层云影响下中海雾
Fig. 8. Relationship between DOP and observed altitude angle under advection fog with maritime cumulus and stratus. (a) Heavy advection fog with maritime cumulus and stratus; (b) moderate advection fog with maritime cumulus and stratus
图 9. 浓海雾和中海雾在云层环境下偏振度随观测高度角的变化。 (a)海洋性积云影响下浓海雾和中海雾;(b)海洋性层云影响下浓海雾和中海雾
Fig. 9. Relationship between DOP and observed altitude angle under heavy and moderate advection fog with cloud layer. (a) Heavy and moderate advection fog with maritime cumulus;(b) heavy and moderate advection fog with maritime stratus
5 结论
从地面对天空进行偏振探测时,天空光的偏振分布对探测结果存在着不可忽视的影响。针对传统天空光偏振分布的研究不适用于复杂海洋环境的问题,对海雾环境进行研究,将传统简化的两层大气结构扩展为三层,由Mie散射算法对云层、海雾层粒子分布特性进行描述,采用基于矢量传输方程的RT3方法仿真天空光经过大气、云层、海雾层的下行辐射过程,研究不同性质云层和不同海雾浓度环境下太阳子午线位置处偏振度与观测高度角的变化关系。研究结果表明:1) 采用基于RT3的仿真建模方法得到的大气散射天空偏振分布与基于Rayleigh散射的大气理论模型得到的分布趋势完全一致。2) 在太阳与观测高度角间的散射角为90°时,可得到最大偏振度值;相反地,在太阳位置处,可得最小偏振度值;在大气Arago中性点处,可得到偏振度的谷值。3) 在波段480 nm下,随着云层和海雾粒子尺寸的增大,以及光学厚度的增加,偏振度逐渐增大。以上研究表明在海面对天空进行偏振观测时,太阳高度角、观测高度角及天气环境等条件的选择至关重要。
在实际的应用中,可用海面上搭载光谱仪与偏振片组合的方式,通过调节仪器的观测位置,测量不同观测方位下天空中某一点的下行辐射偏振度值,最终用仿真结果加以验证。但由于海洋环境复杂多变,需通过大量的测试实验对天空光偏振模式的仿真结果进行修正。本研究对典型海洋环境下天空偏振特性进行了研究,给出了天空偏振特性的变化趋势,为下一步复杂海洋环境下天空光偏振模式的测试提供了理论支撑。
[1] 王家成, 乔延利, 杨世植, 等. 海洋大气气溶胶光学模型参数的相关性研究[J]. 光学学报, 2011, 31(11): 1101005.
[2] 崔岩, 张西光, 周鑫昌, 等. 气溶胶对天空光偏振分布的影响[J]. 光学学报, 2019, 39(6): 0601001.
[3] 吴良海, 张骏, 范之国, 等. 多次散射因素影响下天空偏振光模式的解析模型[J]. 物理学报, 2014, 63(11): 114201.
Wu L H, Zhang J, Fan Z G, et al. An analytical mo del for skylight polarization pattern with multiple scattering[J]. Acta Physica Sinica, 2014, 63(11): 114201.
[4] 崔岩, 赵金勇, 关乐, 等. 黄海海域天空光偏振分布仿真与测试[J]. 光学学报, 2017, 37(10): 1001004.
[5] Zhao H J, Ji Z, Zhang Y, et al. Mid-infrared imaging system based on polarizers for detecting marine targets covered in sun glint[J]. Optics Express, 2016, 24(15): 16396-16409.
[6] Chorfi H, Ayadi K, Gader R, et al. Morphological characterization of particles by the intensity and polarization of the scattered radiation[J]. Optik, 2018, 154: 251-257.
[7] Jin Z, Charlock T P, Rutledge K, et al. Analytical solution of radiative transfer in the coupled atmosphere-ocean system with a rough surface[J]. Applied Optics, 2006, 45(28): 7443-7455.
[8] Tapimo R. Kamdem H T T, Yemele D. Discrete spherical harmonics method for radiative transfer in scalar planar inhomogeneous atmosphere[J]. Journal of the Optical Society of America A, 2018, 35(7): 1081-1090.
[9] Oelkrug D, Brun M, Hubner P, et al. Penetration of light into multiple scattering media: model calculations and reflectance experiments. Part II: the radial transfer[J]. Applied Spectroscopy, 2013, 67(4): 385-395.
[10] 李姣, 段民征, 覃军. 基于逐次散射法的海洋——大气矢量辐射传输模式[J]. 遥感技术与应用, 2014, 29(2): 181-188, 211.
Li J, Duan M Z, Qin J. A coupled successive order scattering vector radiative transfer model for ocean-atmosphere system[J]. Remote Sensing Technology and Application, 2014, 29(2): 181-188, 211.
[11] Zhang S, Zhan J T, Fu Q, et al. Propagation of linear and circular polarization in a settling smoke environment: theory and experiment[J]. Applied Optics, 2019, 58(17): 4687-4694.
[12] 张肃, 战俊彤, 付强, 等. 非均匀烟雾环境对偏振光传输特性的影响[J]. 光学学报, 2018, 38(7): 0729002.
[13] Wright J B, Kemme S A, et al. Superior signal persistence of circularly polarized light in polydisperse, real-world fog environments[J]. Applied Optics, 2018, 57(19): 5464-5473.
[14] Zeng X W, Chu J K, Cao W D, et al. Visible-IR transmission enhancement through fog using circularly polarized light[J]. Applied Optics, 2018, 57(23): 6817-6822.
[15] Evans K F, Stephens G L. A new polarized atmospheric radiative transfer model[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 1991, 46(5): 413-423.
[16] 王威, 褚金奎, 崔岩, 等. 基于矢量辐射传输的大气偏振建模[J]. 中国激光, 2013, 40(5): 0513001.
[17] 张颖, 张熠, 赵慧洁. 多种天气条件下的天空光偏振模型[J]. 红外与毫米波学报, 2017, 36(4): 453-459.
[18] 饶瑞中. 现代大气光学[M]. 北京: 科学出版社, 2015: 152- 155.
Rao RZ. Modern atmospheric optics [M]. 3rd ed. , Beijing: Science Press, 2015: 152- 155.
[19] Hess M, Koepke P, Schult I. Optical properties of aerosols and clouds: the software package OPAC[J]. Bulletin of the American Meteorological Society, 1998, 79(5): 831-844.
[20] 石广玉. 大气辐射学[M]. 北京: 科学出版社, 2007: 389- 390.
Shi GY. Atmospheric radiation[M]. Beijing: Science Press, 2007: 389- 390.
Article Outline
张肃, 战俊彤, 付强, 姜会林. 复杂海洋环境下天空光偏振特性模拟研究[J]. 光学学报, 2020, 40(22): 2201001. Su Zhang, Juntong Zhan, Qiang Fu, Huilin Jiang. Simulation Research on Sky Polarization Characteristics Under Complicated Marine Environment[J]. Acta Optica Sinica, 2020, 40(22): 2201001.