-
Mie散射理论是各向同性的均匀球体对平面电磁波散射的精确解,研究沙尘粒子的散射特性时,可将沙粒视为球形粒子利用Mie理论进行分析计算。根据Mie散射理论,单个粒子对入射波的消光效率因子
${Q_{\rm e}}$ 、散射效率因子${Q_{\rm s}}$ 和吸收效率因子${Q_{\rm a}}$ 分别为:$${Q_{\rm e}} = \frac{2}{{{x^2}}}\sum\limits_{n = 1}^\infty {(2n + 1)[{Re} ({a_n} + {b_n})]} $$ (1) $${Q_{\rm s}} = \frac{2}{{{x^2}}}\sum\limits_{n = 1}^\infty {(2n + 1)({{\left| {{a_n}} \right|}^2} + {{\left| {{b_n}} \right|}^2})} $$ (2) $${Q_{\rm a}} = {Q_{\rm e}} - {Q_{\rm s}}$$ (3) 式中:
$x = 2{\text{π}}r/\lambda$ ,$r$ 为粒子半径,λ为波长;an和bn为Mie散射系数。效率因子的计算涉及到沙粒的复折射率,复折射率是描述粒子光学特性的一个重要参数,与粒子的构成、尺寸和波长有关。复折射率由实部和虚部两部分组成,其实部表示粒子对入射光散射能力的大小,值越大反映散射能力越大;虚部表示粒子对入射光吸收能力的大小,绝对值越大反映吸收能力越大,文中采用参考文献[11]提供的沙尘粒子在红外波段(0.86~20 μm)的复折射率数值,其实部和虚部如图1所示。
根据Mie散射理论和沙粒的复折射率,分别在近红外,中波红外和长波红外选取了几个典型波长(λ=1.06、4.0、10.6、20.0 μm),利用MTLAB编程计算得到了散射、吸收和消光三个效率因子与粒径的变化关系,结果如图2所示。可以看出,对于1.06 和4.0 μm的红外波而言,沙粒的散射作用比较显著,尤其是小粒径时(D<10.0 μm),吸收作用很弱,消光主要由散射所致,随着沙尘粒径的增大,波长与粒径的差值也增大,散射作用逐渐减弱,吸收作用增强,二者趋于比较接近的值。对于10.6 μm长波红外,粒径小于2 μm时,吸收和散射都较弱,粒径在3.5~17.5 μm之间,散射强于吸收,而后,随着沙尘粒径的增大,二者趋于比较接近的值。对于20 μm长波红外,粒径小于5 μm时,吸收和散射较弱,但吸收略大于散射,粒径在9 ~25 μm范围时,散射效率因子出现第一峰值,粒径大于50 μm之后,散射和吸收趋于比较接近的值。散射效率因子的第一峰值随波长增大逐渐减弱。总之,对于小尺寸的沙尘粒子,在近、中波红外,消光效应主要是散射作用的结果,对于大尺寸的沙尘粒子,红外消光是吸收和散射共同作用的结果。
-
每种沙尘天气同时都含有多种尺度的粒子,沙尘暴对电磁辐射的作用影响与沙粒尺度分布密切相关。浮尘和扬沙天气的沙尘粒子尺寸较小,沙尘暴的沙粒尺寸较大,其粒径一般在0.01~300 μm范围内,研究表明:可采用对数正态分布函数来描述沙尘粒子的尺度分布,其表达式为[12]:
$$N(r) = {N_0}p(r)$$ (4) 式中:N0为粒子数密度;
$p(r) = N(r)/{N_0}$ 为粒子尺度分布密度函数;$r$ 为沙粒半径。$$p(r) = \frac{1}{{2r\sigma \sqrt {2\text{π} } }}\exp \left\{ { - \frac{{{{[\ln (2r) - {m_0}]}^2}}}{{2{\sigma ^2}}}} \right\}$$ (5) 式中:m0和σ分别为ln(2r)的均值和标准偏差。沙源地不同,测量设备规格不同,测得沙尘暴粒子的尺寸分布也不同,文中选取董庆生和牛生杰分别在塔克拉玛干沙漠、贺兰山沙区等地测量的七种沙粒尺寸分布,其相应参数的统计结果如表1所示[13-14]。
图3是根据公式(5)和表1参数绘出的七种沙尘粒子的尺寸分布函数,图中横坐标D为沙粒直径。可以看出,1~6模态的沙粒尺寸比7模态的沙粒尺寸大许多,其中1黄土沙、2毛乌素两种模态的沙粒尺寸主要分布在5~150 μm范围之间,3腾格里、4塔克拉玛干和5甘肃三种模态的沙粒尺寸分布比较相近,粒径主要分布在30~210 μm范围之间,6海岸模态沙粒尺寸相对最大,粒径主要分布在80~250 μm范围之间;7贺兰山区模态沙粒尺寸较小,粒径主要分布在0.01~6 μm范围之间.
表 1 沙粒尺寸分布参数
Table 1. Sand size distribution parameters
Number Sand source m0 σ 1 Loess –3.08 0.491 2 Maowusu –2.96 0.380 3 Tengger –2.31 0.296 4 Taklimakan –2.26 0.276 5 Gansu –2.19 0.279 6 Coast –1.83 0.166 7 Helan mountain –7.48 0.995 沙尘暴是具有多种尺寸分布的多相粒子群,总的消光是所有尺寸粒子消光共同贡献的结果,忽略多次散射对消光的影响,总的消光系数为:
$$\mu = \int\limits_0^\infty {{\text{π}}{r^2}{Q_{\rm e}}(r){N_0}p(r)} {\rm d}r$$ (6) 式中:粒子数密度N0是一个很难测量的物理量,国内外学者通常借助于光学能见度Vb来表示,二者的关系为[15]:
$${N_0} = \dfrac{{15}}{{4.343 \cdot {V_b}\int\limits_0^\infty {2\text{π} {r^2}p(r)} {\rm d}r}}$$ (7) 除消光系数外,单次散射反照率ω和不对成因子g也是描述粒子散射和消光特性的重要物理量,对于多相粒子群,ω和g的计算式分别为:
$$\omega = \frac{{\int\limits_0^\infty {\text{π} {r^2}{Q_{\rm s}}(r){N_0}p(r)} {\rm d}r}}{{\int\limits_0^\infty {\text{π} {r^2}{Q_{\rm e}}(r){N_0}p(r)} {\rm d}r}}$$ (8) $$g = \frac{{\int\limits_0^\infty {\text{π} {r^2}{Q_{\rm s}}(r)g(r){N_0}p(r)} {\rm d}r}}{{\int\limits_0^\infty {\text{π} {r^2}{Q_{\rm s}}(r){N_0}p(r)} {\rm d}r}}$$ (9) 根据以上各式,在能见度Vb=1 km时,计算得7种模态沙尘暴在0.86~20 μm红外波段的3个消光参量随波长的变化情况如图4和图5所示。
图 4 1~6模态沙尘暴消光参量与波长的关系
Figure 4. Relationship between extinction parameters and wavelength of 1-6 mode sandstorm
图 5 模态7沙尘暴消光参量与波长的关系
Figure 5. Relationship between extinction parameters and wavelength of mode 7 sandstorm
根据图4所示的6种大颗粒沙尘暴的消光特性,可以把6种模态的沙尘暴分成3组,第一组1和2模态由于粒径分布接近,3个消光参量随红外波长的变化基本一致,第二组是3、4和5模态,第三组是6模态。图4(a)表明第一组的两种沙尘暴其反照率在0.86~5 μm波段大于另两组的值,在长波红外波段,反照率又小于另两组的值,但在整个0.86~20 μm波段,反照率ω的值都在0.5~0.6之间,说明散射略强于吸收。图4(b)表明第一组的不对称因子小于另两组,但在整个0.86~20 μm波段内,反照率g的值都大于0.9,表明沙尘暴的前向散射较强。图4(c)表明第一组的消光系数明显大于另外两组的值,第三组沙尘暴的颗粒尺寸相对最大,消光系数相对最小,主要因为沙粒尺寸与该红外波段的波长差异最大的缘故。总之,在0.86~20 μm红外波段内,可见度相同时,这6种大颗粒模态沙尘暴的消光系数除在个别波长有小的震荡外,大体是随红外波长的增大而增大。
图5是小颗粒模态7沙尘暴消光参量随红外波长的变化情况。可以看出:在红外波段波长小于5 μm时,其反照率和不对称因子的值都比较大,反映沙尘暴在近、中波红外前向散射比较强;随红外波长增大,反照率震荡减小,不对称因子递减,说明该沙尘暴散射随波长增大而逐渐减弱且趋于各向同性;在红外波长小于2 μm时,消光系数明显大于其他波段的值,即该模态沙尘暴在近红外的消光能力较强,在中、长红外波段的消光能力较弱,且在13~20 μm波段范围内,消光对波长不敏感。
-
沙尘暴对红外辐射信号强度的影响最终归结为其对信号的衰减程度如何,为此,下面将根据前面消光参量的计算,分析计算衰减率与红外波长、沙尘暴模态、可见度的关系。
-
根据Lamber-Beer定律,当强度为I0的光波通过厚度为L的随机介质后,透射光波强度I为:
$$I = {I_O}\exp ( - \mu \cdot L)$$ (10) 式中:μ为消光系数。衰减率常用分贝(dB/km)表示,定义如下:
$$A = \dfrac{{10}}{L}\log \dfrac{I}{{{I_0}}}$$ (11) 根据公式(10),可得衰减率与消光系数的关系为:
$$ A = 4.343\mu $$ (12) 公式(12)表明:消光系数μ的大小就决定了衰减率的大小。但这是以单次散射为前提的,即粒子之间彼此独立地散射光波。在沙尘暴天气条件下,大气中的沙尘粒子比较浓密,能见度较低,这时,单次散射很难真实地反映沙尘暴天气的实际散射衰减情况,因此,就需要考虑多次散射的影响。
-
红外辐射在沙尘暴中的传输不可避免存在多重散射情况,利用Monte Carlo[16-17]法计算电磁辐射在沙尘暴中多次散射现象具有独特的优点。Monte Carlo方法的基本思想是把电磁辐射看成由很多光子组成的光子束,电磁辐射的传输问题就转化为光子的传输问题。每一个光子在沙尘暴中的传递过程,都可以由计算机产生随机数来确定,通过对光子的散射方向、行进路径进行抽样,实现对光子的跟踪,光子权重随着在沙尘暴中的散射次数而减小,直到光子权重小于阈值、或光子离开沙尘暴、或到达接受界面为止。通过对大量的光子进行跟踪,最后对到达接受界面的光子权重进行统计,就可以获得稳定的透过率及衰减率估值。
设光波沿z方向垂直射入到厚度为H的沙尘暴中,光子在相邻两次碰撞间随机运动步长由平均自由程的累计概率分布决定,即:
$$L = - \frac{1}{{{\mu _{}}}}\ln \xi $$ (13) 式中:ξ为(0,1)上均匀分布的随机数;μ为消光系数。
光子在某点与沙粒碰撞后,新的运动方向由Henyey和Greenstein提出的HG散射相函数来确定,根据HG散射相函数可得散射角θ的抽样值为:
$$\theta = {\cos ^{ - 1}}\left\{ \frac{1}{{2g}}[{(1 + g)^2} - {(\frac{{1 - {g^2}}}{{1 - g + 2g\xi }})^2}]\right\} $$ (14) 式中:g为不对称因子。光子散射方位角φ可认为在(0, 2π)内均匀分布,其抽样值为:
$$\phi = 2\text{π} \xi $$ (15) 设(μx, μy, μz)为光子与沙粒碰撞前的方向余弦,碰撞后光子新的方向余弦可通过坐标变换得到:
$$\begin{split} {{\mu '}_x} = & \dfrac{{\sin \theta }}{{\sqrt {1 - \mu _z^2} }}({\mu _x}{\mu _z}\cos \phi - {\mu _y}\sin \phi ) + {\mu _x}\cos \theta \\ {{\mu '}_y} = & \dfrac{{\sin \theta }}{{\sqrt {1 - \mu _z^2} }}({\mu _x}{\mu _z}\cos \phi + {\mu _x}\sin \phi ) + {\mu _y}\cos \theta \\ {{\mu '}_z} = & - \sin \theta \cos \phi \sqrt {1 - \mu _z^2} + {\mu _x}\cos \theta \end{split} $$ (16) 根据加权Monte Carlo的思想,光子与粒子碰撞后,它的部分能量就会被吸收,其权重W会减小:
$${W_i} = {W_{i - 1}} \cdot \omega $$ (17) 式中:ω为单次散射反照率;Wi-1为第i次碰撞前的权重;Wi为碰撞后的权重,设初始权重W0=1。光子在沙尘暴中行进步长L后到达下一碰撞点,其坐标变为:
$$ x' = x + {\mu _x}L,\;y' = y + {\mu _y}L,\;z' = z + {\mu _z}L $$ (18) 如果光子所在位置z'<0或z'>H,则表示光子从沙尘暴中逃逸出去,上述跟踪过程结束;如果光子的权因子小于设定的阈值(文中取10−6),则表示光子被吸收,跟踪过程也结束。如果光子既没被吸收,也没逃逸出去,重复以上步骤,直至光子到达接受界面。
跟踪N个光子(取N=106),平均透过率的估值为:
$$T = \frac{1}{N}\sum\limits_{n = 1}^N {{W_n}} $$ (19) 式中:Wn为第n个出射光子的权重,则得光波在沙尘暴中传播单位距离的衰减率计算式为:
$$A = - \frac{{10}}{H}\log T$$ (20) -
根据以上单次散射和多重散射的计算方法,选取模态1和模态7两种不同尺度分布的沙尘暴,以典型波长1.06 和10.6 μm为例,基于Matlab编程,首先计算比较了两种散射相应的衰减率与可见度的关系,结果如图6所示。可以看出,每种模态沙尘暴相应的两种散射衰减率计算结果存在明显差异,基于多重散射的衰减率小于同一条件下的单次散射衰减率,差异随可见度的增大而减小;模态7沙尘暴对于1.06 μm近红外的两种散射衰减率计算结果差异最大,对于10.6 μm长波红外衰减率计算结果差异相对较小,主要缘于模态7沙尘暴在近红外的反照率很高(图5(a)),则多重散射现象显著,按照单次散射计算误差大。因此,后面的分析计算只采用3.2节反映多重散射的Monte Carlo方法。
图 6 单次散射和多重散射计算结果比较
Figure 6. Comparison of calculation results based on single scattering and multiple scattering
图7是在可见度一定时(V=1 km),基于多重散射的Monte Carlo方法计算得到的7种模态沙尘暴红外衰减率与波长的关系。显然,可见度和波长相同时,大颗粒1~6模态沙尘暴比第7种小颗粒模态沙尘暴的红外衰减强度大,尤其是在12.5 ~20 μm范围的长波红外,第7种模态沙尘暴的衰减率远小于前6种模态沙尘暴的衰减率。图7(a)表明,6种大颗粒模态沙尘暴的衰减率除有小幅度震荡外,基本呈现的是随着波长的增大而增大的趋势。波长处于0.86~5.8 μm范围时,1和2模态沙尘暴小于其余4种模态沙尘暴的衰减率;在8~20 μm长波红外,1和2模态沙尘暴明显大于其余4种模态沙尘暴的衰减率。
图 7 7种模态沙尘暴衰减率与波长的关系
Figure 7. Relationship between attenuation rate and wavelength of 7 kinds of modes
从图7(b)可以看出,第7种小颗粒模态沙尘暴的衰减率随波长变化有明显的起伏,在0.86 μm近红外附近,衰减率的值相对较大,但随波长增大至2.5 μm时,衰减率急剧下降至最小,而后在2.5~6.5 μm范围内衰减率对波长不很敏感,在7.9~9.0 μm范围内衰减率随波长增大又快速增大,在10~12.5 μm范围内衰减率随波长增大又快速减小,而后在13~20 μm范围内衰减率对波长不敏感,变化幅度不大。
比较图4(c)、图5(c)和图7,不难看出,在可见度相同时,同种模态沙尘暴的消光系数和衰减率随波长的变化情况基本趋势大致一致,但并不完全相同,主要原因是消光系数的计算是仅仅考虑单次散射的情况,而图7各模态衰减率的计算是考虑了多重散射得到的结果,沙尘暴中的沙粒浓度大,造成的多重散射现象比较明显,Monte Carlo方法计算时不仅要利用消光系数μ,还要涉及反照率ω和不对称因子g,因此,图7的衰减率比消光系数更能全面反映沙尘暴对红外辐射的影响。
图8是对于几个典型红外波长,基于Monte Carlo方法计算得到的7种模态沙尘暴衰减率随可见度的变化情况。结果表明,6种大颗粒模态的沙尘暴对1.06 μm的近红外衰减率基本相同,对于10.6 μm和20 μm的长波红外而言,1~2模态的沙尘暴比其余几种模态的沙尘暴衰减程度更强一些,各种模态沙尘暴衰减率随可见度的增大而减小。相同条件下,7模态沙尘暴的衰减率相对最小。
Analysis of extinction characteristics of sandstorm to infrared radiation(Invited)
-
摘要: 基于Mie散射理论和Monte Carlo方法,分析了国内七种模态沙尘暴对0.86~20 μm波段红外辐射的消光和衰减特性。研究结果表明:对于小尺寸的沙尘粒子,在近、中波红外,消光效应主要是散射作用的结果;对于大尺寸的沙尘粒子,红外消光是吸收和散射共同作用的结果。比较了沙尘暴单次散射和多重散射产生的衰减率差异,多重散射的衰减率小于同一条件下的单次散射衰减率,差异随可见度的增大而减小,基于Monte Carlo方法计算的多重散射衰减率比消光系数更能全面反映沙尘暴对红外辐射强度衰减的影响。六种大颗粒模态沙尘暴比小颗粒模态沙尘暴的红外衰减强度大,衰减率基本随着波长的增大而增大;小颗粒模态沙尘暴的衰减率随波长的变化有明显的起伏,在7.9~12.5 μm范围内有峰值,在13~20 μm范围内衰减率对波长不敏感。
-
关键词:
- 红外辐射 /
- 沙尘暴 /
- 消光系数 /
- Monte Carlo /
- 衰减率
Abstract: The extinction and attenuation characteristics of seven modes of sandstorm in China for 0.86 -20 μm band infrared radiation were analyzed based on Mie scattering theory and Monte Carlo method. The results show that for the small size dust particles, the extinction effect is mainly the result of scattering in the near and middle infrared, and for the large size dust particles, the infrared extinction is the result of absorption and scattering. The attenuation rate difference of single scattering and multiple scattering of sandstorm are compared. The attenuation rate of multiple scattering is less than that of single scattering under the same condition, and the difference decreases with the increase of visibility. The multiple scattering attenuation rate calculated based on Monte Carlo method is more comprehensive than the extinction coefficient to reflect the influence of sandstorm on infrared radiation intensity attenuation. The infrared attenuation intensity of six kinds of large particle mode sandstorm is greater than that of small particle mode sandstorm, and the attenuation rate increases with the increase of wavelength. The attenuation rate of the small particle mode duststorm fluctuates obviously with the change of wavelength, with a peak value in the range of 7.9–12.5 μm, and is not sensitive to the wavelength in the range of 13–20 μm.-
Key words:
- infrared radiation /
- sandstorm /
- extinction coefficient /
- Monte Carlo /
- attenuation rate
-
表 1 沙粒尺寸分布参数
Table 1. Sand size distribution parameters
Number Sand source m0 σ 1 Loess –3.08 0.491 2 Maowusu –2.96 0.380 3 Tengger –2.31 0.296 4 Taklimakan –2.26 0.276 5 Gansu –2.19 0.279 6 Coast –1.83 0.166 7 Helan mountain –7.48 0.995 -
[1] Ansari A J, Evans B G. Microwave propagation in sand and dust storms [J]. IEE Proc F-Commun Radar & Signal Process, 1982, 129(5): 315−322. [2] Gang Hong, Ping Yang, Fuzhong Weng. Microwave scattering properties of sand particles: Application to the simulation of microwave radiances over sandstorms [J]. Journal of Quantitative Spectroscopy & Radiative Transfer, 2008, 109: 684−702. [3] Dong Qunfeng, Li Yingle, Xu Jiadong, et al. Effect of sand and dust storms on microwave propagation [J]. IEEE Transactions on Antennas & Propagation, 2013, 61(2): 910−916. [4] 周旺, 周东方, 侯德亭, 等. 微波传输中沙尘衰减的计算与仿真[J]. 强激光与粒子束, 2005, 17(8): 1259−1262. Zhou Wang, Zhou Dongfang, Hou Deting, et al. Calculation and simulation of sand and dust attenuation in microwave propagation [J]. High Power Laser and Particle Beams, 2005, 17(8): 1259−1262. (in Chinese) [5] 李曙光, 刘晓东, 侯蓝田. 沙尘暴对低层大气红外辐射的吸收和衰减[J]. 电波科学学报, 2003, 18(1): 43−47. doi: 10.3969/j.issn.1005-0388.2003.01.009 Li Shuguang, Liu Xiaodong, Hou Lantian. The absorption and extinction of infrared radiation in lower atmosphere by sand storm [J]. Chinese Journal of Radio Science, 2003, 18(1): 43−47. (in Chinese) doi: 10.3969/j.issn.1005-0388.2003.01.009 [6] 吴振森, 由金光, 杨瑞科. 激光在沙尘暴中的衰减特性研究[J]. 中国激光, 2004, 31(9): 1075−1080. doi: 10.3321/j.issn:0258-7025.2004.09.012 Wu Zhensen, You Jinguang, Yang Ruike. Study on laser ttenuation character in sand and dust storms [J]. Chinese Journal of Lasers, 2004, 31(9): 1075−1080. (in Chinese) doi: 10.3321/j.issn:0258-7025.2004.09.012 [7] 李学彬, 徐青山, 魏合理. 1次沙尘暴天气的消光特性研究[J]. 激光技术, 2008, 32(6): 566−567. Li Xuebin, Xu qingshan, Wei Heli. Extinction character of one sand and dust blowing [J]. Laser Technology, 2008, 32(6): 566−567. (in Chinese) [8] 董群锋, 李应乐, 许家栋, 等. 沙尘的红外衰减特性研究[J]. 装备环境工程, 2012, 9(2): 4−6. doi: 10.3969/j.issn.1672-9242.2012.02.002 Dong Qunfeng, Li Yingle, Xu Jiadong, et al. Study on infrared attenuation characteristics of sand and dust [J]. Equipment Environmental Engineering, 2012, 9(2): 4−6. (in Chinese) doi: 10.3969/j.issn.1672-9242.2012.02.002 [9] 杨瑞科, 朱传帅, 刘科祥. 红外波在沙尘暴中的后向散射增强[J]. 红外与激光工程, 2017, 46(1): 0104006. doi: 10.3788/IRLA201746.0104006 Yang Ruike, Zhu Chuanshuai, Liu Kexiang. Backscattering enhancement for infrared wave in sand and dust storm [J]. Infrared and Laser Engineering, 2017, 46(1): 0104006. (in Chinese) doi: 10.3788/IRLA201746.0104006 [10] 杨玉峰, 李挺, 李建勋, 等. 沙尘在FSO常用红外波段的散射特性研究[J]. 红外与激光工程, 2017, 46(6): 0604004. doi: 10.3788/IRLA201746.0604004 Yang Yufeng, Li Ting, Li Jianxun, et al. Study on scattering characteristics of sand and dust in the infrared waveband commonly used in FSO [J]. Infrared and Laser Engineering, 2017, 46(6): 0604004. (in Chinese) doi: 10.3788/IRLA201746.0604004 [11] 石广玉. 大气辐射学[M]. 北京: 科学出版社, 2007: 367-368. Shi Guangyu. Atmospheric Radiation Science[M]. Beijing: Science Press, 2007: 367-368 (in Chinese). [12] Ahmed A S. Role of particle-size distributions on millimetre-wave propagation in sand/dust storms [J]. Microwaves, Antennas and Propagation, 1987, 134(1): 55−59. doi: 10.1049/ip-h-2.1987.0011 [13] 董庆生. 我国典型沙区中沙尘的物理特性[J]. 电波科学学报, 1997, 12(1): 15−21. Dong Qingsheng. Physical characteristics of the sand and dust in different deserts of China [J]. Chinese Journal of Radio Science, 1997, 12(1): 15−21. (in Chinese) [14] 牛生杰, 孙继明, 陈跃, 等. 贺兰山地区春季沙尘气溶胶质量浓度的观测分析[J]. 高原气象, 2001, 20(1): 82−87. doi: 10.3321/j.issn:1000-0534.2001.01.014 Niu Shengjie, Sun Jimin, Chen Yue, et al. Observation and analysis of mass concentration of dust and sand aerosol in spring in Helanshan Area [J]. Plateau Meteorology, 2001, 20(1): 82−87. (in Chinese) doi: 10.3321/j.issn:1000-0534.2001.01.014 [15] 杨瑞科, 苏振玲, 刘科祥. 沙尘暴多重散射对毫米波衰减影响研究[J]. 电波科学学报, 2008, 23(3): 530−533. Yang Ruike, Su Zhenling, Liu Kexiang. Influence of multiple scattering on attenuation of millimeter wave propagation in sand and dust storm [J]. Chinese Journal of Radio Science, 2008, 23(3): 530−533. (in Chinese) [16] 王红霞, 竹有章, 田涛, 等. 激光在不同类型气溶胶中传输特性研究[J]. 物理学报, 2013, 62(2): 024214. doi: 10.7498/aps.62.024214 Wang Hongxia, Zhu Youzhang Tian Tao, et al. The characteristics of laser transmission in different types of aerosol [J]. Acta Phys Sin, 2013, 62(2): 024214. (in Chinese) doi: 10.7498/aps.62.024214 [17] Seyoung M, Dongyun K, Eunji S. Monte Carlo study of coherent diffuse photon transport in a homogeneous turbid medium: a degree-of-coherence based approach [J]. Appl Opt, 2008, 47(3): 336−345. doi: 10.1364/AO.47.000336