光子学报, 2020, 49 (8): 0801001, 网络出版: 2020-11-27   

非球形灰霾的紫外脉冲回波特性 下载: 655次

Characteristics of Ultraviolet Pulse Echo on Non-spherical Haze
作者单位
1 西安理工大学 自动化与信息工程学院, 西安 710048
2 陕西省智能协同网络军民共建重点实验室, 西安 710000
摘要
为了区分不同物理特性的非球形灰霾粒子,基于T矩阵理论和蒙特卡洛方法,建立了“日盲”紫外光后向散射探测灰霾模型,仿真了非球形灰霾条件下的紫外光后向散射过程,并分析了不同宽度的紫外脉冲后向散射回波的峰值功率和半高全宽等特征.研究结果表明,当灰霾浓度在500 μg/m3以内时,紫外脉冲后向散射回波峰值功率随着灰霾浓度的增大而增大,且不同浓度灰霾的脉冲回波峰值功率之间呈线性关系,脉冲回波半高全宽随着灰霾浓度的增大而减小.对于椭球形和圆柱形灰霾粒子,在相同浓度条件下,当粒子的形变量越小时,脉冲回波的峰值功率和半高全宽越大;对于切比雪夫形灰霾粒子,在相同浓度条件下,当形变参数和波纹参数越大时,脉冲回波的峰值功率和半高全宽越大.本文研究结果可为紫外光探测灰霾浓度以及区分非球形灰霾粒子提供依据.
Abstract
In order to distinguish the non-spherical haze particles with different physical characteristics, the haze detection model of "solar-blind" ultraviolet backscatter is established to simulate the process of UV backscatter under the condition of non-spherical haze, based on T-matrix theory and Monte Carlo method. The peak power and the full width at half maximum of UV pulse with different widths are analyzed. The simulation results show that when the haze concentration is less than 500 μg/m3, the peak power of backscatter echo of UV pulse increases with the increase of haze concentration, and there is a linear relationship between the peak power of pulse echo of different concentrations of haze, and the full width at half maximum of pulse echo decrease with the increase of haze concentration. For ellipsoidal and cylindrical haze particles, under the same concentration, the smaller the shape variable of particles is, the larger the peak power and the full width at half maximum of pulse echo are. For chebyshev haze particles, under the same concentration, the larger the deformation parameters and ripple parameters are, the larger the peak power and the full width at half maximum of pulse echo are. The results of this study can provide a basis for the detection of haze concentration and the discrimination of non-spherical haze particles by ultraviolet light.

0 引言

近年来,大气灰霾现象被认为是部分城市的灾害性天气,原因是:1)灰霾影响司机和行人视线,容易引起交通事故;2)灰霾易附带有毒、有害物质,容易导致疾病发生[1];3)灰霾会造成大气光学特性的变化[2],对所在地区的无线光通信产生不可忽视的影响.因此,对灰霾进行有效而准确的监测是十分有必要的,监测结果可以为治理灰霾污染以及提高无线光通信质量提供参考.

目前环境监测部门对大气气溶胶的测量方法主要有:重量法、β射线吸收法、天平振荡法和光散射法[3]等.重量法、β射线吸收法和天平振荡法都需要对气溶胶进行取样,且取样过程和实际操作都很麻烦[4],不能满足实时测量的需求.光散射测量方法是利用气溶胶粒子对光的散射作用,通过分析散射光的特性来反演被测粒子的浓度和形态等特征[5].在气溶胶粒子的光散射领域,国内外许多学者做过了大量研究.2017年,左晨泽[6]等利用13种常见的气溶胶分析了粒子光散射法测量中,复折射率对测量结果的影响.同年,HU S [7]等建立了一种基于时域多分辨方法的非球形气溶胶散射模型.结果表明,该模型计算得到的相函数的相对误差较小.2018年,MIROSHNICHENKO A E[8]等根据球形粒子散射结果推广得出,对于任意有限散射体,其最大部分吸收截面不能超过三维均匀球面和二维圆柱的相应值.2018年,HE S [9]等利用偏振光蒙特卡洛模型模拟了红外光在海面气溶胶粒子中传播的散射过程,生成了不同环境条件和不同探测参数条件下的海面偏振红外图形,用于指导海洋环境探测.2019年,陈柯蓓[10]等利用华北地区的激光雷达数据,对比分析了不同灰霾天气下的大气悬浮颗粒物浓度的垂直分布.结果表明,大气悬浮颗粒物浓度的垂直分布随天气的变化而剧烈变化.2019年,李素文[11]等根据所建立的反演模型获得了粒径在0.1~1.25 μm范围内的气溶胶数密度谱分布,该方法可为今后大气气溶胶探测提供参考.目前,气溶胶粒子的光散射测量主要集中在可见光波段,可见光散射测量方法容易受背景光的影响,需要在暗盒中进行.由于大气臭氧对200~280 nm波段的紫外光具有很强的吸收作用,该波段紫外光又被称为“日盲”紫外光,它具有抗干扰能力强、隐秘性高的特点[12-13],因此,我们将“日盲”紫外光用于灰霾探测领域,可以弥补可见光测量方法容易受背景光干扰的不足.由于实际大气中的灰霾粒子多是非球形的[14],在利用“日盲”紫外光探测灰霾时,需要分析不同物理参数的灰霾粒子的紫外光散射特性差异,并根据此差异来区分非球形灰霾粒子.

由于266 nm的紫外激光器具有操作简单、高效率、光束质量好和体积小等优点,该波长的紫外激光器常被应用于紫外探测和无线紫外光通信领域[15].本文使用266 nm的“日盲”紫外光作为探测光源,通过分析灰霾粒子的紫外光后向散射回波特性来探测不同浓度及形态的灰霾.紫外光探测灰霾属于非取样法,可以满足全天候实时测量的需求,实时测量不仅能够保证对灰霾进行实时监控,而且还可以为分析无线光通信信道的实时变化提供参考.

1 紫外光探测灰霾粒子的理论基础

1.1 灰霾粒子的物理特性

1.1.1 粒子尺度参数

灰霾粒子的紫外光散射特性受粒子尺度参数的影响,粒子尺度参数定义为αD/λD为灰霾粒子直径,λ为入射光的波长.灰霾是由多个粒径不同、形态各异的粒子组成,一般条件下,灰霾粒子的粒径和浓度如表 1所示.

表 1. 灰霾的半径和浓度[16]

Table 1. Radius and density of haze [16]

TypeRadius/μmConcentration/cm3
Aitken10-3~10-210-4~102
Haze10-2~110~103
Droplet1~1010~100

查看所有表

1.1.2 粒子复折射率

粒子复折射率是表征粒子散射和吸收的一个光学常数,其实部和虚部分别反映了粒子对光的散射和吸收能力.对于灰霾粒子而言,由于其组成成分多为硫酸盐、有机碳等物质,其实部变化不大.

1.1.3 粒子尺度谱分布

大气中的灰霾并不是单一组成成分的,常见的有沙尘型灰霾粒子、水溶型灰霾粒子和黑碳型灰霾粒子.因此,本文采用三对数正态分布来模拟灰霾粒子的分布情况,三对数正态分布表达式n(r)为

式中,Ci表示归一化常数,σi是标准偏差,$r_{\rm{c}}^i$是不同类型灰霾粒子的平均半径,i表示灰霾的第i种成分.表 2给出了三种灰霾粒子的尺度谱分布[17].

表 2. 灰霾的尺度谱分布

Table 2. Scale spectrum distribution of haze

TypeSize distribution
WASO(σ=0.81,rc=0.03 μm)
BC(σ=0.69,rc=0.011 μm)
DUST(σ=0.92,rc=0.47 μm)

查看所有表

1.2 非球形灰霾粒子的紫外光散射理论

一般将非球形灰霾粒子等效成椭球形、圆柱形和切比雪夫形这几种常见的粒子.图 1(a)是椭球形灰霾粒子示意图,其中,ba分别为垂直半轴和水平半轴,a/b的大小决定椭球灰霾粒子的形状,大于1为扁椭球形灰霾粒子,小于1是长椭球形灰霾粒子.图 1(b)是圆柱形灰霾粒子示意图,圆柱形灰霾粒子的形状由底面圆直径与长度之比D/L决定,D/L小于1是长圆柱形灰霾粒子,等于1是紧实圆柱形灰霾粒子,大于1是扁圆柱形灰霾粒子.图 1(c)是切比雪夫形灰霾粒子示意图,切比雪夫形灰霾粒子是一种直径遵循n阶切比雪夫多项式的非球形灰霾粒子,其形状由形变参数ξ和波纹参数n决定.

图 1. 非球形粒子示意图

Fig. 1. Schematic diagram of non-spherical particle

下载图片 查看所有图片

本文使用T矩阵方法来处理非球形灰霾粒子的紫外光散射问题.T矩阵中的矩阵元仅与粒子的物理特性和几何结构有关,与入射场和散射场的传播方向无关.入射场Einc和散射场Esca展开式[18]分别为

式中,k1=2π/λ是散射环境介质中的波数,r0为能包围非球形粒子的最小外接球的半径,Mmn(kr)、Nmn(kr)是Henkel函数$h_n^{(1)}(kr)$的矢量球谐函数,RgMmn(kr)和RgNmn(kr)是Bessel函数的矢量球谐函数.入射场的展开系数amnbmn的表示式为

同理,散射场展开系数pmnqmn的表达式为

入射场与散射场的展开系数满足线性关系,用T矩阵表示为

其中,QRgQ是2×2的矩阵,即$\mathit{\boldsymbol{Q}} = \left[ \begin{array}{l}{Q^{11}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {Q^{12}}\\{Q^{21}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {Q^{22}}\end{array} \right]$Q中的各元素与MmnNmn及RgMmn、RgNmn的关系分别为

根据求解的T矩阵,散射振幅矩阵的各个矩阵元表达式分别为

式中,${\alpha _{mnm'n'}} = {{\rm{i}}^{n' - n - 1}}{( - 1)^{m + m'}}{\left[ {\frac{{(2n + 1)(2n' + 1)}}{{n(n + 1)n'(n' + 1)}}} \right]^{1/2}}$,且

求解方程组可得到非球形粒子消光截面Cext和散射截面Csca的表达式为

2 紫外光后向散射探测灰霾模型

2.1 灰霾粒子的消光系数和散射相函数

灰霾粒子的散射系数ksca以及吸收系数kabs的表达式为

式中,Csca(r)和Cabs(r)分别表示单个灰霾粒子的散射截面和吸收截面.rminrmax分别表示灰霾粒子尺度谱半径的最小值和最大值,分别取值0.01 μm和1.25 μm.在实际灰霾环境中,空气中大气分子的Rayleigh散射也对紫外光具有散射和吸收作用,因此,总的消光系数ke=ksca+kabs+ksr+kar,总的散射系数ks=ksca+ksr,总的吸收系数ka=kabs+kar.其中,kar是Rayleigh吸收系数,ksr是Rayleigh散射系数.散射相函数表示光子在与散射粒子发生碰撞后在各个方向的散射强度,根据散射粒子的不同,可分为大气分子Rayleigh散射相函数和灰霾粒子的散射相函数.Rayleigh散射相函数见式(19).

式中,β是光子散射夹角,γ是模型参数.

为了更准确地体现灰霾粒子对紫外光的后向散射特性,本文使用改进的RH-G散射相位函数.RH-G散射相位函数不仅能反映球形灰霾粒子的后向散射特征,还能很好的模拟非球形灰霾粒子的后向散射特征[19].服从粒子尺度谱分布的灰霾粒子的散射相函数PRHG(βgn(r))的表达式为

综上所得,紫外光探测灰霾时,总的散射相函数表达式为

2.2 紫外脉冲收发模型

图 2所示,紫外脉冲收发模型是将紫外脉冲的发射端TX和接收端RX固定在同一位置上.图 2中将收发位置设为坐标原点,发射端向灰霾粒子发出紫外脉冲光束,紫外光与灰霾发生散射,Sn为第n次散射点.紫外光波长为266 nm,峰值功率由发端发射的光子数决定,紫外脉冲的光束发散半角和接收端接收孔径等参数在初始化中进行设定.

图 2. 紫外光后向散射探测灰霾模型

Fig. 2. The model of haze detection by UV backscattering

下载图片 查看所有图片

为了更准确地计算光子在灰霾中传输时,每次散射的新位置和散射方向,本文对原有xyz坐标系进行了坐标系变换.如图 2所示,通过绕z轴旋转αt角度,绕y'轴旋转$\frac{\rm \pi }{2} - {\theta _{\rm{t}}}$角度,得到新的坐标系x'y'z',新坐标系的z'轴与接收视场椎体Cr的中心轴一致.坐标系xyz与变换后的新坐标系x'y'z'的关系表达为

图 2所示,某个传输光子的传输方向与x'y'z'坐标系的z'轴所成夹角为θ,投影与x'轴所成夹角为ψ,则光子的传输方向可由(θψ)表示.设ξ(1)ξ(2)是[0, 1]上均匀分布的随机变量,γ是紫外脉冲发端发散半角,则它们的关系表达式分别为

因此,光子在新坐标系x'y'z'中的传输方向余弦表达式为

由式(22)可知,光子在xyz坐标系中的传输方向余弦表达式为

2.3 光子散射路程

紫外光在灰霾中传输时与散射粒子随机发生碰撞,光子两次碰撞之间的路程rt

式中,ξ(t)是在[0, 1]上平均分布的随机变量.由于将发射端设在坐标原点,所以光子的第一次散射点S1xyz坐标系中的坐标可表示为

式中,r0是发射端与第一次散射点S1的距离.光子在S1处发生散射的散射角β1

式中,ξ(s)是在[0, 1]上平均分布的随机变量,P(cosβ)是散射相函数.

图 2可知,再以第一次散射点S1为坐标原点,进行第二次坐标系变换,得到坐标系x''y''z''.根据光子在S1处发生散射时的散射角β1,利用式(29)和(30)可以得到光子下一次散射点S2xyz坐标系中的位置坐标和散射角β2,按照这种推导,依次可以求出光子多次散射后的位置坐标和散射角βn.

2.4 后向接收端脉冲响应

紫外光子需要满足以下3个条件才能最终到达后向接收端:1)最后一次散射点在接收端接收视场角内;2)光子传输方向指向接收面;3)光子必须经过灰霾粒子散射,且传输距离为rn.

图 2所示,最后一次散射点Sn在后向接收端接收视场角内即(ζΦ),光子才有可能到达后向接收端,光子指向后向接收端接收面的概率P1n

式中,A是接收孔径面积,P(cosβn)是光子最后一次散射时的散射相函数.

紫外光子经过灰霾粒子散射,传输距离为rn的概率P2n

所以,紫外光子经过散射后能够被后向接收端捕获的概率Pn

式中,Wn是光子到达后向接收端之前的存活概率,其表示式为

紫外光子在灰霾中传输时,发生散射的概率为ks/ke,将光子第一次散射之前的存活概率设为W0W0=ks/ke),随着散射次数增加,紫外光子存活概率变小.紫外光子在到达后向接收端之前,可能经历N次散射,所以一个紫外光子经过N次散射到达后向接收端的总概率PN的表达式为

因此,紫外光子从发射端发出直至最终被后向接收端捕获,总的散射路程L

紫外光子在灰霾中散射传输所经历的时间tn=L/cc是光速.取均匀时间节点ti,Δt是时间间隔.则发端发出M个光子,每个光子最多经历N次散射,到达后向接收端的脉冲响应h(ti)的表达式为

式中,Pjn是第j个光子经历n次散射到达后向接收端的概率.

3 仿真结果与分析

3.1 非球形灰霾粒子仿真条件

非球形灰霾粒子的尺度谱分布n(r)服从对数正态分布,由非球形灰霾粒子的有效半径reff和有效方差νeff决定,表达式分别为

对于三种非球形灰霾,椭球形灰霾粒子的半轴比a/b分别取1/2、1/3、2/1、2/3、3/1、3/2;圆柱形灰霾粒子的直径与长度比D/L分别取1/2、1/3、2/1、2/3、3/1、3/2;切比雪夫形灰霾粒子的形变参数ξ分别取0.1、0.05和0.02;波纹参数n分别取2、4和6.非球形灰霾粒子的分布也服从三对数正态分布,其粒子尺度谱分布的概率密度函数如式(40)所示.其中,沙尘型、水溶型和黑碳型等三种灰霾粒子所占的比例分别为20%、70%和10%,相对湿度为60%,在紫外光波长为266 nm的条件下,混合形灰霾粒子的复折射率m=1.53+0.005i.在此,我们以其中的沙尘型灰霾为例,有效半径reff取0.47 μm,有效方差veff取0.006 9.通过计算得出的相关散射参数见表 3.

表 3. 灰霾粒子的相关散射参数表

Table 3. Relevant scattering parameters of haze particles list

HazeCext/μm2Csca/μm2Cabs/μm2AlbedoAsymmetry factor
a/b=1/22.144 22.008 40.135 80.936 60.779 1
a/b=1/32.124 71.994 30.130 40.938 50.798 1
a/b=2/12.136 22.001 70.134 50.9370.747 3
a/b=2/32.194 32.052 60.141 70.935 50.779 6
a/b=3/12.065 21.944 30.120 90.941 40.841 2
a/b=3/22.189 42.046 30.143 10.934 60.815 7
D/L=1/22.048 11.925 10.1230.939 70.792 4
D/L=1/32.024 21.897 60.126 60.937 20.821 6
D/L=2/12.034 31.922 20.112 10.944 80.784 1
D/L=2/32.098 31.9710.127 30.939 30.777 1
D/L=3/12.021 11.898 30.122 80.939 10.833 6
D/L=3/22.071 71.944 20.127 50.938 60.774 2
Chebyshev haze2.211 32.068 20.143 10.935 30.773 7
Spherical haze2.225 62.081 50.144 10.935 20.831 3

查看所有表

图 3是不同半轴比a/b的椭球形灰霾粒子的散射相函数随角度的变化图.由图 3可以看出,散射相函数整体上呈先减小后增大的趋势,不同a/b的椭球形灰霾粒子的散射相函数在后向散射角度上区分相对明显.当散射角大于160°时,球形灰霾粒子的散射强度略大于椭球形灰霾粒子的散射强度.图 4是不同直径长度比D/L的圆柱形灰霾粒子的散射相函数随角度的变化图.与椭球形灰霾粒子类似,不同D/L的圆柱形灰霾粒子的散射相函数在后向散射角度上区分也相对明显.当散射角大于130°时,不同D/L的圆柱形灰霾粒子的散射强度呈上升趋势,且整体上小于球形灰霾粒子的散射强度.

图 3. 不同a/b的椭球形灰霾的散射相函数

Fig. 3. Scattering phase function for ellipsoidal

下载图片 查看所有图片

图 4. 不同D/L的圆柱形灰霾的散射相函数

Fig. 4. Scattering phase function for cylindrical

下载图片 查看所有图片

图 5是具有不同形变参数和波纹参数的切比雪夫形灰霾粒子的散射相函数随角度的变化图.从整体上看,散射相函数呈先减小后增大的趋势,且在后向散射角度上区分相对明显.由图 5(a)可以看出,当散射角大于160°时,形变参数越大,切比雪夫形灰霾粒子的散射强度相对越大,且大于球形灰霾粒子的散射强度.由图 5(b)可以看出,当散射角大于150°时,波纹参数越大,切比雪夫形灰霾粒子的散射强度相对越大,且大于球形灰霾粒子的散射强度.

图 5. 不同参数的切比雪夫形灰霾的散射相函数

Fig. 5. Scattering phase function for Chebyshev haze with different parameters

下载图片 查看所有图片

根据紫外光后向散射探测灰霾模型,利用蒙特卡洛方法,可得到灰霾条件下紫外光后向散射冲击响应函数h(ti),即是紫外光在该灰霾系统中的传递函数.因此,将不同宽度的紫外脉冲信号与单位冲击响应函数进行卷积,可以得到该灰霾条件下,到达后向接收端的脉冲回波信号.发射紫外脉冲为

式中, P0是峰值功率,δ是半功率宽度,本文中设定发射脉冲峰值功率不变,脉冲宽度从2 ns依次增加到40 ns,当发射脉冲宽度为40 ns时,脉冲能量为50 mJ.部分仿真参数见表 4.

表 4. 部分仿真参数

Table 4. Partial simulation parameter

ParameterValue
Laser wavelength/nm266
Laser divergence angle/mrad5
Emission optical lens diameter /mm12
Receiving optical lens diameter/mm25
Emission photon number106
Haze concentration/(μg·m-3)100,200,300,400,500

查看所有表

大气灰霾的质量浓度一般在500 μg/m3以内[20],灰霾的单位体积粒子个数N(个/m3)可以转化为质量浓度W(μg/m3).本文中的灰霾质量浓度分别取100 μg/m3、200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3.

3.2 非球形灰霾粒子的仿真结果分析

3.2.1 椭球形灰霾粒子的紫外脉冲回波特性

基于紫外光后向散射探测灰霾模型和T矩阵方法仿真了不同浓度的椭球形灰霾粒子的紫外脉冲回波特性,以半轴比a/b为1/2的椭球形灰霾粒子为例,紫外脉冲回波如图 6所示.

图 6. 不同浓度条件下的椭球形灰霾的紫外脉冲后向散射回波

Fig. 6. Backscattering echo of ellipsoidal haze with different concentration of UV pulse

下载图片 查看所有图片

图 6可以看出,随着灰霾浓度的升高,椭球形灰霾的回波峰值功率也在逐渐增加,且发射脉冲越窄,回波畸变越明显.回波波形畸变随着发射脉冲宽度的增加而逐渐变小,回波峰值功率随着发射脉冲宽度的增加而逐渐增大,波形逐渐近似为高斯波形.

图 7是半轴比a/b为1/2的椭球形灰霾粒子在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.当灰霾浓度在500 μg/m3以内时,紫外脉冲回波峰值功率随着灰霾浓度的增大而增大,对于中低浓度的椭球形灰霾,脉冲回波峰值功率的区分相对比较明显.在发射脉冲宽度小于20 ns时,浓度为400 μg/m3和500 μg/m3的椭球形灰霾的回波峰值功率相差很小.因此,为了提高探测的准确度,我们可以设置较大的发射脉冲宽度,通过分析紫外脉冲回波峰值功率来区分不同浓度的椭球形灰霾.在发射脉冲宽度相同的条件下,不同浓度的椭球形灰霾的回波峰值功率之间近似呈线性关系,拟合函数的近似表达式为

图 7. 椭球形灰霾的脉冲回波峰值功率

Fig. 7. Peak power of pulse echo for ellipsoidal haze

下载图片 查看所有图片

式中,P是浓度为100 μg/m3的椭球形灰霾对应的回波峰值功率,P200P300P400P500是浓度分别为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的椭球形灰霾对应的回波峰值功率.

图 8是半轴比a/b为1/2的椭球形灰霾粒子在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 8可以看出,发射脉冲宽度每增加2 ns,回波半高全宽约增加3 ns.椭球形灰霾的脉冲回波半高全宽随着灰霾浓度的增大而减小,且灰霾浓度越大,回波半高全宽的差距越小.当发射脉冲宽度为40 ns时,浓度为100 μg/m3的椭球形灰霾对应的回波半高全宽比浓度为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的椭球形灰霾对应的回波半高全宽分别大5 ns、9 ns、12 ns和14 ns.因此,当灰霾浓度在500 μg/m3以内时,我们可以通过分析紫外脉冲回波半高全宽的变化规律来区分不同浓度的椭球形灰霾粒子,回波半高全宽越小,灰霾浓度越大.

图 8. 椭球形灰霾的脉冲回波半高全宽

Fig. 8. Full width at half maximum of pulse echo for ellipsoidal haze

下载图片 查看所有图片

图 9是半轴比a/b分别为1/2、1/3、2/1、2/3、3/1和3/2的椭球形灰霾粒子,在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.

图 9. 不同半轴比的椭球形灰霾的脉冲回波峰值功率

Fig. 9. Peak power of pulse echo for ellipsoidal haze with different half axis ratio

下载图片 查看所有图片

图 9可以看出,椭球形灰霾的回波峰值功率整体上随着灰霾浓度的增大而不断变大.对于不同半轴比的椭球形灰霾,当垂直半轴和水平半轴相差越小时,脉冲回波峰值功率就越大.在半轴比互为倒数的情况下,扁椭球形灰霾比长椭球形灰霾的脉冲回波峰值功率要大.在相同浓度条件下,椭球形灰霾的脉冲回波峰值功率均小于球形灰霾的脉冲回波峰值功率,这是因为球形灰霾粒子对紫外光的后向散射能力强于椭球形灰霾粒子对紫外光的后向散射能力.随着灰霾浓度的升高,不同半轴比的椭球形灰霾的脉冲回波峰值功率差距也在增大,在灰霾浓度分别为100 μg/m3、300 μg/m3和500 μg/m3,发射脉冲宽度为40 ns的条件下,半轴比为3/2的扁椭球形灰霾比半轴比为1/3的长椭球形灰霾的回波峰值功率分别大4.01 μW、6.59 μW和7.51 μW.在高浓度灰霾条件下,半轴比越接近1的椭球形灰霾的回波峰值功率和球形灰霾的回波峰值功率的差距越小.

图 10是半轴比a/b分别为1/2、1/3、2/1、2/3、3/1和3/2的椭球形灰霾粒子,在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.

图 10. 不同半轴比的椭球形灰霾的脉冲回波半高全宽

Fig. 10. Full width at half maximum of pulse echo for ellipsoidal haze with different half axis ratio

下载图片 查看所有图片

图 10可以看出,发射脉冲宽度越大,椭球形灰霾的回波半高全宽就越大,两者之间呈线性关系.在整体上,脉冲回波半高全宽随着灰霾浓度的增大而减小.对于不同半轴比的椭球形灰霾,垂直半轴和水平半轴相差越小,回波半高全宽就越大,但均小于球形灰霾粒子的回波半高全宽.在半轴比互为倒数的情况下,扁椭球形灰霾比长椭球形灰霾的脉冲回波半高全宽要大.因此,我们可以结合紫外脉冲回波半高全宽和回波峰值功率的变化规律来区分不同半轴比的椭球形灰霾.

3.2.2 圆柱形灰霾粒子的紫外脉冲回波特性

图 11是基于紫外光后向散射探测灰霾模型和T矩阵方法仿真了D/L为2/1的圆柱形灰霾粒子的紫外脉冲回波波形图.由图 11可以看出,随着不断增加发射脉冲的宽度,圆柱形灰霾的回波波形畸变逐渐变小,波形逐渐近似为高斯波形.

图 11. 不同浓度条件下的圆柱形灰霾的紫外脉冲后向散射回波

Fig. 11. Backscattering echo of cylindrical haze with different concentration of UV pulse

下载图片 查看所有图片

图 12D/L为2/1的圆柱形灰霾粒子在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.圆柱形灰霾的紫外脉冲回波峰值功率随着灰霾浓度的增大而增大,且发射脉冲宽度越大,回波峰值功率也越大,当发射脉冲宽度相对较大时,回波波形也逐渐趋于平缓.对于中低浓度的圆柱形灰霾,根据脉冲回波峰值功率来区分灰霾浓度的效果较好.不同于椭球形灰霾,圆柱形灰霾在高浓度条件下的脉冲回波峰值功率的差距也相对比较明显,当发射脉冲宽度大于18 ns时,根据回波峰值功率可以有效地对圆柱形灰霾浓度进行区分.在发射脉冲宽度相同的条件下,不同浓度的圆柱形灰霾的紫外脉冲回波峰值功率之间近似呈线性关系,拟合函数的近似表达式为

图 12. 圆柱形灰霾的脉冲回波峰值功率

Fig. 12. Peak power of pulse echo for cylindrical haze

下载图片 查看所有图片

式中,P是浓度为100 μg/m3的圆柱形灰霾对应的回波峰值功率,P200P300P400P500是浓度分别为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的圆柱形灰霾对应的回波峰值功率.

图 13D/L为2/1的圆柱形灰霾粒子在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 13可以看出,圆柱形灰霾的回波半高全宽与发射脉冲宽度之间呈线性关系,发射脉冲宽度越大,回波半高全宽越大.圆柱形灰霾的脉冲回波半高全宽随着灰霾浓度的增大而减小,在中低浓度灰霾条件下,脉冲回波半高全宽的差距相对明显,随着灰霾浓度的升高,脉冲回波半高全宽的差距逐渐减小.当发射脉冲宽度为40 ns时,浓度为100 μg/m3的圆柱形灰霾对应的回波半高全宽比浓度为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的圆柱形灰霾对应的回波半高全宽分别大5 ns、10 ns、12 ns和14 ns.不同浓度的圆柱形灰霾的脉冲回波半高全宽的差距比球形灰霾的脉冲回波半高全宽的差距要小,和半轴比a/b为1/2的椭球形灰霾接近.因此,我们可以通过分析紫外脉冲回波半高全宽的变化规律来区分不同浓度的圆柱形灰霾,脉冲回波半高全宽越小,灰霾浓度越大.

图 13. 圆柱形灰霾的脉冲回波半高全宽

Fig. 13. Full width at half maximum of pulse echo for cylindrical haze

下载图片 查看所有图片

图 14是底面圆直径与长度之比D/L分别为1/2、1/3、2/1、2/3、3/1和3/2的圆柱形灰霾粒子,在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系.

图 14. 不同直径长度比的圆柱形灰霾的脉冲回波峰值功率

Fig. 14. Peak power of pulse echo for cylindrical haze with different diameter to length ratio

下载图片 查看所有图片

图 14可以看出,圆柱形灰霾的脉冲回波峰值功率整体上也是随着灰霾浓度的增大而不断变大,发射脉冲宽度越大,回波峰值功率也越大.对于D/L不同的圆柱形灰霾,底面圆直径与长度越相近,脉冲回波峰值功率就越大.另外,不同于椭球形灰霾,在D/L之比互为倒数的情况下,长圆柱形灰霾比扁圆柱形灰霾的脉冲回波峰值功率要大.在相同浓度条件下,球形灰霾粒子的脉冲回波峰值功率都要高于圆柱形灰霾粒子的脉冲回波峰值功率,这是因为球形灰霾粒子对紫外光的后向散射能力强于圆柱形灰霾粒子对紫外光的后向散射能力.随着灰霾浓度的升高,D/L之比互为倒数的圆柱形灰霾的脉冲回波峰值功率的差距在逐渐减小,特别是对于底面圆直径与长度差异较大的圆柱形灰霾,脉冲回波峰值功率的差距更小.在灰霾浓度为500 μg/m3,发射脉冲宽度为40 ns的条件下,D/L之比为1/3和D/L之比为3/1的脉冲回波峰值功率的差距仅为0.19 μW,也就是说,在高浓度灰霾条件下,D/L互为倒数的圆柱形灰霾粒子的紫外光后向散射强度比较相近.

图 15是底面圆直径与长度之比D/L分别为1/2、1/3、2/1、2/3、3/1和3/2的圆柱形灰霾粒子,在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 15可以看出,圆柱形灰霾的回波半高全宽与发射脉冲宽度也具有良好的呈线性关系.从整体上看,灰霾浓度越大,圆柱形灰霾的回波半高全宽越小.在相同浓度条件下,不同D/L的圆柱形灰霾的脉冲回波半高全宽均小于球形灰霾的脉冲回波半高全宽.对于不同D/L的圆柱形灰霾,底面圆直径与长度差异越小,脉冲回波半高全宽越大.在低浓度灰霾条件下,不同D/L的圆柱形灰霾粒子的回波半高全宽区分相对明显.在D/L互为倒数的情况下,长圆柱形灰霾比扁圆柱形灰霾的脉冲回波半高全宽要大,且随着灰霾浓度的升高,脉冲回波半高全宽的差距在减小.在高浓度灰霾条件下,球形灰霾与圆柱形灰霾的回波半高全宽的差距相对大一些.所以,根据上述规律,我们可以结合紫外脉冲回波半高全宽和回波峰值功率的变化来区分球形与圆柱形灰霾,同样,该规律也可以用于区分不同D/L的圆柱形灰霾.

图 15. 不同直径长度比的圆柱形灰霾的脉冲回波半高全宽

Fig. 15. Full width at half maximum of pulse echo for cylindrical haze with different diameter to length ratio

下载图片 查看所有图片

3.2.3 切比雪夫形灰霾粒子的紫外脉冲回波特性

图 16是形变参数ξ=0.1,波纹参数n=4的切比雪夫形灰霾粒子的紫外脉冲后向散射回波波形图.

图 16. 不同浓度条件下的切比雪夫形灰霾的紫外脉冲后向散射回波

Fig. 16. Backscattering echo of Chebyshev haze with different concentration of UV pulse

下载图片 查看所有图片

图 16可以看出,在灰霾浓度较低和发射脉冲较窄的条件下,切比雪夫形灰霾的回波波形畸变较大,随着发射脉冲宽度的增大及灰霾浓度的升高,切比雪夫形灰霾的回波波形畸变逐渐减小,波形逐渐近似为高斯波形.

图 17是形变参数ξ=0.1、波纹参数n=4的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.由图 17可以看出,当发射脉冲宽度增大到30 ns后,回波峰值功率走势趋于平缓.对于浓度低于200 μg/m3的切比雪夫形灰霾,其脉冲回波峰值功率随灰霾浓度的变化而相差较大.因此,通过脉冲回波峰值功率来区分低浓度切比雪夫形灰霾的效果较好.随着灰霾浓度的升高,不同浓度的切比雪夫形灰霾的脉冲回波峰值功率的差距在逐渐减小.对于400 μg/m3和500 μg/m3的高浓度切比雪夫形灰霾,在发射脉冲宽度由小增大的过程中,它们的回波峰值功率的差距一直维持在较小水平.在发射脉冲宽度相同的条件下,不同浓度的切比雪夫形灰霾的紫外脉冲回波峰值功率之间近似呈线性关系,近似拟合函数表达式为

图 17. 切比雪夫形灰霾的脉冲回波峰值功率

Fig. 17. Peak power of pulse echo for Chebyshev haze

下载图片 查看所有图片

式中,P是浓度为100 μg/m3的切比雪夫形灰霾对应的回波峰值功率,P200P300P400P500是浓度分别为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的切比雪夫形灰霾对应的回波峰值功率.由拟合函数式(44)可以看出,切比雪夫形灰霾与椭球形灰霾在不同浓度条件下的线性相关系数较为相似,但其脉冲回波峰值功率却大于椭球形灰霾的脉冲回波峰值功率.

图 18是形变参数ξ=0.1、波纹参数n=4的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 18可以看出,不同浓度条件下的切比雪夫形灰霾的回波半高全宽与发射脉冲宽度之间近似呈线性关系,随着灰霾浓度的升高,脉冲回波半高全宽在减小,脉冲回波半高全宽的差距也在逐渐减小.当发射脉冲宽度为40 ns时,浓度为100 μg/m3的切比雪夫形灰霾对应的脉冲回波半高全宽比浓度为200 μg/m3、300 μg/m3、400 μg/m3和500 μg/m3的切比雪夫形灰霾对应的脉冲回波半高全宽分别大5 ns、9 ns、11 ns和13 ns.因此,我们可以通过分析紫外脉冲回波半高全宽的变化规律来区分不同浓度的切比雪夫形灰霾,脉冲回波半高全宽越小,灰霾浓度越大.

图 18. 切比雪夫形灰霾的脉冲回波半高全宽

Fig. 18. Full width at half maximum of pulse echo for Chebyshev haze

下载图片 查看所有图片

图 19是波纹参数n=4,形变参数ξ分别为0.1、0.05和0.02的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.由图 19可以看出,当ξ大于0.02时,切比雪夫形灰霾比球形灰霾的脉冲回波峰值功率要大,这是因为切比雪夫形灰霾粒子对紫外光的后向散射能力强于球形灰霾粒子对紫外光的后向散射能力,且ξ越大,其后向散射能力越强,脉冲回波峰值功率也就越大.不同形变的切比雪夫形灰霾的脉冲回波峰值功率随着灰霾浓度的增大而不断增大,且灰霾浓度越大,切比雪夫形灰霾与球形灰霾的脉冲回波峰值功率的差距越大.

图 19. 不同形变参数的切比雪夫形灰霾的脉冲回波峰值功率

Fig. 19. Peak power of pulse echo for Chebyshev haze with different deformation parameters

下载图片 查看所有图片

图 20是波纹参数n=4,形变参数ξ分别取0.1、0.05和0.02的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 20可以看出,切比雪夫形灰霾比球形灰霾的脉冲回波半高全宽要大,在整体上,切比雪夫形灰霾的回波半高全宽随着灰霾浓度的增大而减小.对于ξ不同的切比雪夫形灰霾,脉冲回波半高全宽随着ξ的增大而增大.根据上述规律,我们可以结合紫外脉冲回波半高全宽和回波峰值功率来区分ξ取值不同的切比雪夫形灰霾.

图 20. 不同形变参数的切比雪夫形灰霾的脉冲回波半高全宽

Fig. 20. Full width at half maximum of pulse echo for Chebyshev haze with different deformation parameters

下载图片 查看所有图片

图 21是形变参数ξ取0.1,波纹参数n分别取2、4和6的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波峰值功率和发射脉冲宽度的关系图.由图 21可以看出,不同波纹参数的切比雪夫形灰霾的脉冲回波峰值功率随着灰霾浓度的增大而不断增大.不同波纹参数的切比雪夫形灰霾比球形灰霾的脉冲回波峰值功率要大,且n越大,其后向散射能力越强,脉冲回波峰值功率也就越大.

图 21. 不同波纹参数的切比雪夫形灰霾的脉冲回波峰值功率

Fig. 21. Peak power of pulse echo for Chebyshev haze with different ripple parameters

下载图片 查看所有图片

图 22是形变参数ξ取0.1,波纹参数n分别取2、4和6的切比雪夫形灰霾粒子,在不同浓度条件下的紫外脉冲回波半高全宽和发射脉冲宽度的关系图.由图 22可以看出,在形变参数ξ=0.1的条件下,切比雪夫形灰霾比球形灰霾的回波半高全宽要大,不同波纹参数的切比雪夫形灰霾的脉冲回波半高全宽随着灰霾浓度的增大而减小.对于不同波纹参数的切比雪夫形灰霾,当n越大时,其回波半高全宽也越大.根据上述规律,我们可以结合紫外脉冲的回波半高全宽和回波峰值功率来区分n取值不同的切比雪夫形灰霾.

图 22. 不同波纹参数的切比雪夫形灰霾的脉冲回波半高全宽

Fig. 22. Full width at half maximum of pulse echo for Chebyshev haze with different ripple parameters

下载图片 查看所有图片

4 结论

本文根据紫外光后向散射探测灰霾模型和T矩阵方法仿真了不同浓度条件下,不同形状的椭球形、圆柱形及切比雪夫形灰霾粒子的紫外脉冲回波特性.研究结果表明:灰霾的脉冲回波峰值功率随着灰霾浓度的增大而增大,脉冲回波半高全宽随着灰霾浓度的增大而减小.对于椭球形灰霾,当垂直半轴和水平半轴相差越小,灰霾粒子形状越接近球形时,脉冲回波的峰值功率和半高全宽越大;对于圆柱形灰霾,当底面圆直径与长度相差越小时,脉冲回波的峰值功率和半高全宽越大;对于切比雪夫形灰霾粒子,当ξn越大时,脉冲回波的峰值功率和半高全宽越大.因此,我们可以结合紫外脉冲回波的峰值功率和半高全宽来区分不同浓度及不同形状的非球形灰霾粒子.

参考文献

[1] ZENG Xian-gang, 曾 贤刚, 阮 芳芳, RUAN Fang-fang, PENG Yan-yan, 彭 彦彦. 基于空间网格尺度的中国PM2.5污染健康效应空间分布[J]. 中国环境科学, 2019, 39(6): 2624-2632.

[2] GAO Ai-zhen, 高 爱臻, 胡 利军, HU Li-jun, YANG Hao, 杨 豪. 基于宁波地区大气颗粒物监测激光雷达的一次灰霾过程分析[J]. 气象水文海洋仪器, 2019, 36(2): 26-30.

[3] DAN De-zhong, 但 德忠. 环境空气PM_(2.5)监测技术及其可比性研究进展[J]. 中国测试, 2013, 39(2): 1-5.

[4] LING Tao, HARLEY R A. Changes in fine particulate matter measurement methods and ambient concentrations in California[J]. Atmospheric Environment, 2014, 98(12): 676-684.

[5] YE Chao, 叶 超, 孟 睿, MENG Rui, GE Bao-zhen, 葛 宝臻. 基于光散射的粒子测量方法综述[J]. 激光与红外, 2015, 45(4): 343-348.

[6] ZUO Chen-ze, 左 晨泽, 吕 且妮, LV Qie-ni, GE Bao-zhen, 葛 宝臻. 气溶胶颗粒折射率在光学粒径测量中的影响[J]. 光学精密工程, 2017, 25(7): 1777-1782.

[7] HU S, GAO T, LI H. Light scattering computation model for non-spherical aerosol particles based on multi-resolution time-domain scheme:Model development and validation[J]. Optics Express, 2017, 25(2): 1463-1486.

[8] MIROSHNICHENKO A E, TRIBELSKY M I. Ultimate absorption in light scattering by a finite obstacle[J]. Physical Review Letters, 2018, 120(3): 033902.

[9] HE S, WANG X, XIA R. Polarimetric infrared imaging simulation of a synthetic sea surface with Mie scattering[J]. Applied Optics, 2018, 57(7): 150-159.

[10] CHEN Ke-bei, 陈 柯蓓, 何 梦迪, HE Meng-di, ZHANG Qi, 张 琦. 基于激光雷达技术的局域TSP特性分析[J]. 智能计算机与应用, 2019, 9(4): 258-262.

[11] LI Su-wen, 李 素文, 牟 福生, MOU Fu-sheng, 胡 丽莎, HU Li-sa. 双光路DOAS系统获取大气颗粒物粒谱分布方法研究[J]. 光谱学与光谱分析, 2019, 39(11): 3403-3406.

[12] WANG Qiang, 王 强, 何 廷尧, HE Ting-yao, HUA Deng-xin, 华 灯鑫. 不同紫外波段激发生物气溶胶荧光雷达系统性能分析[J]. 光子学报, 2017, 46(12): 133-140.

[13] ZHAO Tai-fei, 赵 太飞, 冯 艳玲, FENG Yan-ling, KE Xi-zheng, 柯 熙政. "日盲"紫外光通信网络中节点覆盖范围研究[J]. 光学学报, 2010, 30(8): 2229-2235.

[14] ZHENG Kai-duan, CHEN Jian, ZHOU Jie. Research on monitoring of a severe haze pollution over Yangtze River Delta based on multi-source remote sensing data[J]. Remote Sensing for Land & Resources, 2018, 30(1): 224-232.

[15] DONG Ye, 董 野. 关于266nm紫外激光器的探究[J]. 通讯世界, 2015, 24(12): 313-314.

[16] 柯熙政.紫外光自组织网络理论[M].北京:科学出版社, 2011:56-59.

[17] WANG Zhong-ting, 王 中挺, 陈 良富, CHEN Liang-fu, LI Qing, 厉 青. 灰霾组分的多角度偏振反射特性模拟[J]. 光学学报, 2015, 35(9): 9-17.

[18] 饶瑞中.现代大气光学[M].北京:科学出版社, 2012:31-35.

[19] CHEN Chen, 程 晨, 徐 青山, XU Qing-shan, ZHU Lin, 朱 琳. 非球形气溶胶粒子散射相函数经验公式[J]. 光谱学与光谱分析, 2019, 39(1): 7-13.

[20] ZHOU C, SHEN X, LIU Z. Simulating aerosol size distribution and mass concentration with simultaneous nucleation, condensation/coagulation, and deposition with the GRAPES-CUACE[J]. Journal of Meteorological, 2018, 32(2): 265-278.

赵太飞, 王世奇, 张健伟, 崔真. 非球形灰霾的紫外脉冲回波特性[J]. 光子学报, 2020, 49(8): 0801001. Tai-fei ZHAO, Shi-qi WANG, Jian-wei ZHANG, Zhen CUI. Characteristics of Ultraviolet Pulse Echo on Non-spherical Haze[J]. ACTA PHOTONICA SINICA, 2020, 49(8): 0801001.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!