-
根据米氏散射理论,消光效率因子Qe为[9]:
$ \begin{array}{l} {Q_{\rm{e}}}(\alpha , n) = \frac{2}{{{\alpha ^2}}}{\rm{Re}}\sum\limits_{m = 1}^\infty {} \left( {2m + 1} \right)\left( {{a_m} + {b_m}} \right) = \\ \frac{4}{{{\alpha ^2}}}{\rm{Re}}\left[ {{s_1}\left( 0 \right)} \right] \end{array} $
(1) 式中, n为折射率,且n=n1-n2i,是粒子的复折射率;$ \alpha = \frac{{2{\rm{ \mathsf{ π} }}r}}{\lambda }$为粒子的尺度参量;r为粒子半径;am和bm为Mie散射系数[10],表达式如下:
$ \left\{ \begin{array}{l} {a_m} = \frac{{{\mathit{\Psi} _m}\prime \left( {n\alpha } \right){\mathit{\Psi} _m}\left( \alpha \right) - n{\mathit{\Psi} _m}\left( {n\alpha } \right){\mathit{\Psi} _m}\prime \left( \alpha \right)}}{{{\mathit{\Psi} _m}\prime \left( {n\alpha } \right){Z_m}\left( \alpha \right) - n{\mathit{\Psi} _m}\left( {n\alpha } \right){Z_m}\prime \left( \alpha \right)}}\\ {b_m} = \frac{{n{\mathit{\Psi} _m}\prime \left( {n\alpha } \right){\mathit{\Psi} _m}\left( \alpha \right) - {\mathit{\Psi} _m}\left( {n\alpha } \right){\mathit{\Psi} _m}\prime \left( \alpha \right)}}{{n{\mathit{\Psi} _m}\prime \left( {n\alpha } \right){Z_m}\left( \alpha \right) - {\mathit{\Psi} _m}\left( {n\alpha } \right){{\rm Z}_m}\prime \left( \alpha \right)}} \end{array} \right. $
(2) 式中,
$ \left\{ \begin{array}{l} {\mathit{\Psi} _m}\left( \alpha \right) = \frac{{\sqrt {{\rm{ \mathsf{ π} }}\alpha } }}{2}{{\rm{J}}_{m + 1/2}}\left( \alpha \right)\\ {{\rm Z}_m}(\alpha ) = \frac{{\sqrt {{\rm{ \mathsf{ π} }}\alpha } }}{2}{{\rm{H}}_{m + 1/2}}(\alpha ) \end{array} \right. $
(3) 式中, ${{\rm{J}}_{m + 1/2}}(\alpha ), {{\rm{H}}_{m + 1/2}}(\alpha ) $分别表示半整阶的贝塞耳函数和第2类汉克函数[11]。Ψm′(α)和Zm′(α)分别为Ψm(α)和Zm(α)的导数。
1μm, 2μm和3μm波长正常温度下的折射率为1.322+3×10-6i,1.297+1.101×10-3i和1.352+2.721×10-1i。根据消光效率因子的计算公式[12],对雨滴直径在0μm~200μm之间,近红外辐射分别为1μm, 2μm和3μm波长的消光效率因子计算如图 1所示。
(1) 近红外辐射波段,1μm波长的辐射在雨滴粒子直径最小的位置(D=2.089μm),消光效率因子达到最大值3.973,其次是2μm波长的辐射在雨滴粒子直径D=2.089μm时, 消光效率因子取得最大值3.911,最后是3μm波长的辐射在雨滴直径D=5.360μm时,消光效率因子取得最大值4.036。
(2) α=πD/λ为粒子的尺度参量,计算最消光效率因子最大值的位置与近红外辐射波长的比值可得,1μm波长尺度参量α1=6.5651;2μm波长尺度参量α2=6.4241;3μm波长尺度参量α3=5.6135。由此可以看出,尺度参量相对比较接近,通过(2)式也可以解释此现象,即粒子的尺度参量和折射率共同决定散射系数,进而决定雨滴粒子的消光效率因子。尺度参量α1, α2, α3比较接近却并不相同,原因在于近红外波长不同时,折射率也不相同,通过消光效率最大值以及曲线形状的不同也恰恰反映出折射率的影响。
(3) 消光效率因子的曲线出现了一系列的极大值和极小值。极值变化幅度较大,是由于散射光和入射光干涉导致的干涉结构;极值变化较小且具有一定周期性, 是由于粒子“谐振”而导致的[12]。
(4) 观察发现:3μm波长震荡的局部极大值大于2μm波长的局部极大值,2μm波长震荡的局部极大值大于1μm波长的局部极大值;同时,3μm波长震荡的局部极小值小于2μm波长的局部极小值,2μm波长震荡的局部极小值小于1μm波长的局部极小值[13]。即3μm波长的消光效率曲线的包络线包含了2μm波长的消光效率曲线的包络线,2μm波长的消光效率曲线的包络线包含了1μm波长的消光效率曲线的包络线。
-
LAWS和PARSONS在美国华盛顿地区对不同类型降雨的雨滴尺寸分布进行了广泛的测量[14],对相同降雨量的雨滴尺寸分布进行了平均,得到了离散性的雨滴尺寸分布的数据表,即不同降雨量的平均雨滴尺寸分布Laws-Parsons分布数据表[15]。根据Laws-Parsons分布数据表作图,如图 2所示,为不同降雨量条件下的雨滴直径与所占百分比的关系曲线。
有学者首先对Laws-Parsons分布数据表中的雨滴直径间隔为0.25mm的区间进行了平均,得到雨滴直径的平均值,雨滴直径平均值、所占体积百分比以及降雨量构成3维坐标,然后将所占百分比看成是雨滴直径平均值和降雨量的二元函数进行了拟合,得到了雨滴体积分布的公式[16]:
$ \begin{array}{l} f\left( {D, R} \right) = \left( {33.44 \pm 1.8} \right){R^{ - 1.28}}{D^{5.93}} \times \\ {\rm{exp}}( - 0.538{R^{ - 0.186}}D)({{\rm{m}}^{ - 3}}\cdot{\rm{m}}{{\rm{m}}^{ - 1}}) \end{array} $
(4) 式中,D为雨滴直径(mm),R为降雨量(mm/h)。由(4)式可以看出,其拟合的雨滴谱分布形式非常复杂,实践经验证明,这个雨滴谱模型在不同地理区域内并没有非常好的适用性。因此本文中提出了建立雨滴谱分布数据库的新方法,即根据雨滴谱分布数据,考虑降雨量固定情况下,雨滴直径平均值和所占体积百分比构成的一元函数关系式。利用MATLAB软件绘制降雨量R分别为0.25mm/h, 1.25mm/h, 2.5mm/h, 12.5mm/h, 25mm/h和50mm/h折线图,如图 2所示。
由观察发现,不同降雨量条件下雨滴谱分布具有正态分布的3个特点:(1)集中性。曲线的高峰位于正中央;(2)对称性。曲线左右对称,曲线两端永远不与横轴相交;(3)均匀变动性。曲线左右两侧逐渐均匀下降。因此文中选择正态分布模型拟合雨滴谱分布。雨滴谱分布的拟合,首先将不同降雨量条件下,雨滴直径和所占百分比数据导入工作区MATLAB的workspace作为变量(x,y),接着在曲线拟合工具箱中分别选择自变量x和因变量y,然后在拟合函数类型中选择Gaussian分布即正态高斯分布,得到降雨量分别为0.25mm/h, 1.25mm/h, 2.5mm/h, 12.5mm/h, 25mm/h和50mm/h时的拟合雨滴谱分布图形,如图 3所示。
得到的雨滴谱分布函数N(D)(单位为m-3·mm-1)的公式见下:
$ \begin{array}{l} N(D) = 27.44{\rm{exp}}\left[ { - {{\left( {\frac{{D - 0.9349}}{{0.5059}}} \right)}^2}} \right], \\ (R = 0.25{\rm{mm}}/{\rm{h}}) \end{array} $
(5) $ \begin{array}{l} N(D) = 20.66{\rm{exp}}\left[ { - {{\left( {\frac{{D - 1.267}}{{0.5968}}} \right)}^2}} \right], \\ (R = 1.25{\rm{mm}}/{\rm{h}}) \end{array} $
(6) $\begin{array}{l} N(D) = 17.6{\rm{exp}}\left[ { - {{\left( {\frac{{D - 1.41}}{{0.7677}}} \right)}^2}} \right], \\ (R = 2.5{\rm{mm}}/{\rm{h}}) \end{array} $
(7) $ \begin{array}{l} N(D) = 13.3{\rm{exp}}\left[ { - {{\left( {\frac{{D - 1.904}}{{1.025}}} \right)}^2}} \right], \\ (R = 12.5{\rm{mm}}/{\rm{h}}) \end{array} $
(8) $ \begin{array}{l} N(D) = 11.81{\rm{exp}}\left[ { - {{\left( {\frac{{D - 2.156}}{{1.159}}} \right)}^2}} \right], \\ (R = 25{\rm{mm}}/{\rm{h}}){\rm{ }} \end{array} $
(9) $ \begin{array}{l} N(D) = 10.43{\rm{exp}}\left[ { - {{\left( {\frac{{D - 2.449}}{{1.319}}} \right)}^2}} \right], \\ (R = 50{\rm{mm}}/{\rm{h}}) \end{array} $
(10) 通过拟合的雨滴谱分布函数N(D)可以发现,降雨量R越大,雨滴直径D的平均值越大,雨滴直径D的方差也越大,图形上表现为:雨滴谱分布曲线向右侧移动,并且曲线形状更加扁平。
-
近红外辐射在雨中进行传播的过程中,并不是只和单个雨滴发生作用,而是与多个雨滴或者雨滴群发生作用。因此在求解计算时,必须得到单位体积内雨滴粒子群对近红外辐射的衰减系数β(dB/km),β反映了近红外辐射在降雨中传输时衰减大小,计算公式如下式所示:
$ \begin{array}{l} \beta = \int_0^\infty {} {\sigma _{\rm{e}}}N(D){\rm{d}}D = \int_0^{6.5} {} {\rm{ \mathsf{ π} }}{\left( {\frac{D}{2}} \right)^2} \times \\ {Q_{\rm{e}}}(\alpha , n)N(D){\rm{d}}D \end{array} $
(11) 式中, N(D)是雨滴谱分布,σe是雨滴的消光截面。根据图 1中消光效率因子与雨滴直径的关系,由于雨滴粒子相对于近红外波段较大[17],尺度参量很大,因而消光效率因子近似为2[18]。
将研究者DAVID拟合的雨滴谱分布公式(4)式和雨滴谱分布数据库中的分布公式(5)式~(10)式代入(11)式,即可以计算出降雨量分别为0.25mm/h,1.25mm/h,2.5mm/h,12.5mm/h,25mm/h和50mm/h时,近红外辐射在雨中的衰减系数如图 4所示。
观察图 4中DAVID和本文中拟合的雨滴谱分布数据库计算得到的衰减系数可以发现,本文中计算的衰减系数略大于DAVID模型计算的衰减系数,曲线递增趋势相同,随着降雨量的增加,在降雨量相对较小时,随着降雨量的增加,衰减系数会以一个较大的斜率迅速增加;而当降雨量增大到一定强度时,随着降雨量的继续增加,衰减系数仍然会继续增加,但斜率相对较小,曲线近似为指数模型。从衰减系数的数值上来看,随着降雨量的增加,近红外辐射在雨中的衰减也随之增大[19],具体体现降雨量为50mm/h时的衰减是降雨量约为0.25mm/h衰减系数的6.8倍,衰减更为明显。
在实际应用过程中,衰减系数与降雨强度的函数关系在工程中具有非常重要的应用。因而基于图 4中降雨量和利用本文中拟合出雨滴谱分布模型计算的衰减系数,分别采用了拟合曲线1和拟合曲线2两种曲线模型进行拟合,并进行了方差分析。拟合曲线1和曲线2的表达式分别见下:
$ y = {A_1}{\rm{exp}}\left( { - \frac{x}{{{t_1}}}} \right) + {y_0} $
(12) $ y = {A_1}{\rm{exp}}\left( { - \frac{x}{{{t_1}}}} \right) + {A_2}{\rm{exp}}\left( { - \frac{x}{{{t_2}}}} \right) + {y_0} $
(13) 式中,t1,t2, A1, A2和y0是拟合曲线的参量。
原曲线、(12)式拟合曲线和(13)式拟合曲线如图 5所示。可以看出, 两种曲线的拟合效果均比较好,进一步进行方差分析可知,拟合曲线1的方差和相关系数为123.2和0.9840, 拟合曲线2的方差和相关系数为2.869和0.9996,拟合曲线2的相关系数大于拟合曲线1的相关系数,拟合曲线2的方差小于拟合曲线1的方差,因此拟合曲线2的效果较好。
计算得到,衰减系数和降雨强度的函数关系式采用拟合曲线2的结果,其中A1=-251.8,t1=48.92,y0=353.6, A2=-70.85,t2=2.848。由此,衰减系数β和降雨量R拟合的函数关系见下:
$ \begin{array}{l} \beta = - 251.8{\rm{exp}}\left( { - \frac{R}{{48.92}}} \right) - 70.85 \times \\ {\rm{exp}}\left( { - \frac{R}{{2.848}}} \right) + 353.6{\rm{ }}({\rm{dB}}/{\rm{km}}) \end{array} $
(14)
基于雨滴谱分布数据库的近红外辐射衰减研究
Study on near infrared radiation attenuation based on distribution database of raindrops
-
摘要: 为了研究降雨对近红外辐射的衰减,定量评估降雨对近红外设备的影响,在雨滴谱分布二元函数模型建立的基础上,提出了一种雨滴谱分布数据库的建立方法,对LAWS和PARSONS在美国华盛顿地区测量的雨滴谱分布数据进行拟合,建立雨滴谱分布数据库。基于雨滴谱分布数据库,结合米散射理论和衰减系数公式,得到了1μm,2μm和3μm波长近红外辐射的消光效率因子与雨滴直径的关系曲线,并计算了在降雨量分别为0.25mm/h,1.25mm/h,2.5mm/h,12.5mm/h,25mm/h和50mm/h条件下的近红外辐射衰减系数。结果表明,不同降雨量条件下,对于雨滴谱分布的函数关系式,正态分布具有更好的拟合效果;消光效率曲线的包络线随波长的增大而增大;拟合衰减系数和降雨强度的函数关系发现,指数函数具有很好的拟合效果。该计算结果对研究降雨对近红外辐射的衰减具有重要意义。Abstract: In order to study the attenuation of near infrared radiation by rainfall and quantitatively estimate the effects of rainfall on near infrared devices, on the basis of the establishment of two-element function model of raindrop distribution, a method for establishing the database of raindrop distribution was proposed. After the fitting of raindrop distribution data measured by LAWS and PARSONS in Washington region of the United States, raindrop distribution database was established. And then, combining with Mie scattering theory and attenuation coefficient formula, the relationship curves between extinction efficiency factor and raindrop diameter under 1μm, 2μm and 3μm wavelength near infrared radiation were obtained. Under the condition that rainfall was respectively 0.25mm/h, 1.25mm/h, 2.5mm/h, 12.5mm/h, 25mm/h and 50mm/h, the attenuation coefficients of near infrared radiation were calculated. The results show that, under different rainfall conditions, normal distribution has better fitting effect on function relation of raindrop spectral distribution. The envelope of extinction efficiency curve increases with the increasing of wavelength. After the fitting of function relation between attenuation coefficient and rainfall intensity, it is found that, the exponential function has the good fitting effect. The calculated results are important for studying the attenuation of near infrared radiation by rainfall.
-