大气湍流中光传播的数值模拟

大气湍流中光传播的数值模拟



马保科郭立新 吴振森

1,2 11

，

（1.西安电子科技大学，陕西西安 710071 2.西安工程大学，陕西西安 710048 ）

摘 要 光在大气湍流中传播时，受大气分子、气溶胶等粒子的相互作用，将发生光束扩展、漂移和

相干性退化等大气湍流效应，这些因素严重影响了光波的远场特性。文章从大气湍流中光传播的理论

研究入手，分析了如何构造较为合理的大气湍流相位屏。进而采用McGlamery算法，对Kolmogorov

谱下的大气湍流随机相位屏进行了数值模拟，并分析了光波从发射机经湍流大气传播到达接收机时的

远场变化特性。研究表明，大气湍流的存在对光的远场传播质量造成很大的影响，研究结果也为大气

湍流中与光传播相关的工程应用及自适应光学技术的完善提供了参考。

关键词 大气湍流；McGlamery算法；相位屏模拟； 大气结构常数；

中图分类号 TP391 文献标识码 A

1 引言

大气湍流是一个相当复杂的随机媒质系统，虽然物理学界对湍流的研究已经历了相

当漫长的历史，但因涉及的因素千头万绪，其间的相互作用和关系也错综复杂，人们对

其物理本质至今未能做到较为清楚的认识。因此，光在大气湍流中传播问题的研究仍存

在理论和实验上的挑战[1,2]。通常，当光在湍流大气中传播时，光束截面内包含着许多

的大气漩涡，这些漩涡各自对照射到它的那一部分光束形成衍射作用，可导致光束的强

度和相位随机变化，进而表现出光束扩展，大气闪烁和相位起伏等大气湍流效应，从而

严重降低了接收机的接收效率。目前，突破大气湍流的影响仍是光在随机介质中传播所

要解决的关键问题[3]。早在20世纪中期，苏联的Obukhov便采用Rytov平缓微扰法由

实验反演湍流特征。在闪烁的饱和现象被发现之后，物理学界又将Markov近似引入求

解光场的统计矩，研究大气湍流下的光场特性[1]。然而，在中等起伏条件下，目前仍没

有找到很好的解析处理方法。由于数值模拟能够从光的传播过程出发，较为清楚地反映

出所涉及问题的物理本质，因而成为研究湍流效应的主要方法[4]。本文采用McGlamery

算法[5]，对Kolmogorov谱下的大气随机相位屏进行了数值模拟，进而结合

Huygens-Fresnel原理，模拟了在有无大气湍流的情况下，接收机处光场的变化特性。

2 大气湍流中光的传播

在折射率为的随机媒质中，一束波长，波数为（）的单色波的电场

n



k

k2



E

由Maxwell波动方程来描述

[1,4]

Ekn(r)E2(Elnn)0

222

(1)

收稿日期: 2010年3月14日 收到修改稿日期：2010年 月 日

基金项目：教育部科技重点项目（105164）

通信作者: 马保科（1972），男，在读博士，副教授，主要从事随机介质中波传播方面的研究。Email: baokema2006@

1

r1n()n(r)0

其中，，（，且）为湍流折射率，(1)式左端最后一项

n(r)1n(r)

1

11

反映了偏振特性，当波长远小于湍流的内尺度（）时，此项可以忽略不计，有



l

0

Ekn(r)E0

222

(2)

因而，电场的任一分量波动方程为

E

Ekn(r)E0

222

(3)

如果介质的非均匀尺度远大于波长，可认为只存在前向小角度散射而无后向散射。对如

图1沿方向的传播进行旁轴近似处理，并将光场写为，得

z

Euexp(ikz)

uu

2

u2ikk(n1)u0

2

222

(4)

zz

22222

其中，为横向算符。(4)式还可进一步化简为

xy





(ikikH)(ikikH)uikH,u0







(5)

zzz



222

其中，算子，。如果折射率不随传播距离而变化（或变化较小）

H1k(n1)



则(5)式中忽略交换算子项，得到如下的广义抛物型方程

u

ik(H1)u

(6)

z

通常，算子不同的近似可得不同的抛粅型方程，而算子最简单的近似即为其Taylor

HH

展式

222

H1（k(n1)）/2



(7)

当湍流场的折射率满足时，(7)式化简为

n(r)1n(r)

1

H12kn

22

1

(8)

将(8)式代入(6)式得

ui

uiknu

2

1

(9)

z2k

或

uuu

22

2ik2knu0

2

1

(10)

22

xyz

(9)式即为忽略后向散射而具有抛粅型近似的波动方程。如果仅考虑介质折射率起伏对场

的作用，则它的右边只需保留与折射率有关的第二项，这时其解反映了在光的传播方向

上由于积分光学路径所导致的相位调制，即

z

u(r,z)u(r,z)exp(ikn(r,)d)

'





'

1







z

(11)

u(r,z)exp(iS)

'

(11)式中，如果折射率起伏引起的相位变化足够小，则可将真空传播和介质相位

S

调制看作是相互独立并同时完成的两个过程。这样，如图1可将连续随机介质分割为一

系列厚度为的平行平板，位于平板前的光场根据(9)式的真空解传播至平板的后面，

z

然后被该平板引起的相位调制；这个场再经同样的真空传播和相位调制传播至下一个平

板，依次形成最终的光场。也即将光在湍流介质中的传播等效为光在真空中通过一系列

2

薄的相位屏后的传播。

图1 随机介质中光传播的相位屏模型

Fig .1 The pha screen model in the atmosphere turbulence

这样，通过第个相位屏后的光场由(9)式得

i



i

z

i

u(r,z)expdzexp(iS(r,z))u(r,z)

iii1





2

(12)



2k

z

i1

其中，，由于的随机特性，不可能利用解析的方法

S(r,z)kn(r,z)dz

i1



S(r,z)

i

z

i1

z

i

获得上式的解，而只能利用数值方法。

3 大气湍流随机相位屏的数值模拟

湍流介质中光传播的特殊性就在于介质折射率是随机起伏的，数值模拟的一个关键

也在于如何构造合理的相位屏，来正确反映介质折射率的变化特性，通常所构造的相位

屏必须满足：

[4,5]

（1）、相位屏所代表的平板厚度应足够小，以确保其对场的振幅没有明显影响而

z

只影响相位，即

2

(13)

z/



n

2

其中，为湍流折射率起伏均方差。



n

（2）、光场的变化特性与相位屏的构造方法应无关，即平板厚度应大于湍流介质非均

匀元尺度，即

zL

0

(14)

这里，为湍流外尺度。

L

0

（3）要使光在平板内的传播满足几何光学近似，则Fresnel尺度应小于湍流内尺

l

F

度，即

l

0

2

zl/

0



(15)

目前，大气湍流随机相位屏的模拟方法较多。这里主要利用McGlamery算法来模

[6]

拟大气湍流相位屏，算法步骤如下：

1）生成二维具有高斯分布的复随机数矩阵。



k

3

2）根据大气折射率功率谱函数，生成二维的功率谱密度函数矩阵 。

3) 将功率谱密度函数矩阵求算术根再乘以复随机数矩阵，得到一个相位均匀分布在



k



,



，振幅受功率谱密度函数调制的复随机数矩阵，即



Y()K()()



(16)

4）对矩阵进行逆Fourier变换得到其空间域形式，即



Y()





Y(l)IFFTY()



Y()exp(j2l)dl







0

L

(17)

5）将矩阵分解为实部和虚部。这样，每一部分均独立代表一种随机的大气湍流相



Y(l)

位屏



ScreenReY(l)ScreenImY(l)

12

， (18)

在模拟相位屏时，相关参数的选择为：网格间距或，屏间

[4,5]

xl/3

0

xL/



2

距。 ，网格长度或（为发射机孔径尺寸）

zl/

0



NxL/



0

Nx2D

D

6

10

S

t

r

u

c

t

u

r

e

F

u

n

c

t

i

o

n

D

f

(

r

)

Simulant value

Theoretical value

10

5

10

4

10

3

101010

-101

Distance r (m)

图2 模拟生成的大气湍流相位屏 图3 相位结构函数的理论与模拟值比较

(Fig.2 The simulant atmosphere turbulent pha screen,Fig.3 Compare the simulant and the theoretical value

of the pha structure function )

图2是利用McGlamery算法模拟生成的Kolmogorov谱下的大气湍流相位屏，由于

通常单层相位屏上相位分布的重复性较高，这里我们采用了多屏叠加的思路来尽量减小

相位的重复性。模拟屏上的相位结构函数定义为

[7]

D(r)(r)(r)







2



[Screen(r)Screen(r)]



r1

n

2

(19)

n

其中，同时，与Kolmogorov表示相位屏上处的相位，为屏上的网格数。

Screen(r)

r

n

大气湍流功率谱相对应的相位结构函数理论值为

[2]

D(r)6.88(rr)



0

5/3

(20)

4

通常，大气湍流的统计特性可以用相位结构函数来描述，因此，可以将相位结构函

数作为验证模拟的相位屏正确与否的判断标准。图3比较了相位结构函数的理论值与模

拟值。其中，湍流的内、外尺度分别为米和米，大气相干长度。

l0.01mL1m

00

r8mm

0

可见，在垂直于传播方向的平面内横向距离米（即）的范围内，相位结构

r2

rL5

0

函数的理论值与Kolmogorov谱下的模拟值吻合较好，但当时，两者相差较大。

rL5

0

因此，在湍流谱及上述参数选定的情况下，为了达到模拟大气湍流的准确性和避免大尺

度湍流起伏所带来的误差，相位屏的尺寸最好选为大小。

L5L5

00

4 大气湍流下传播光场的数值模拟

图4 发射机平面上的模拟光斑

Fig.4 The stimulant facula on the transmitter plane

结合上述Kolmogorov谱下模拟生成的相位屏以及利用Huygens-Fresnel原理，下面

就光在大气湍流中水平传播到达接收机处的光场进行数值模拟。这里，相关参数的选择

为：发射光波为的平面波,发射孔径为的方形孔径，传播距离约



0.6328m

D20cm

10, 湍流的内、外尺度分别为1和10，传播路径上分布40个相位屏, 垂直于传

km

cm

m

播方向的平面被分成了100×100个小网格, 网格宽度选择为湍流的内尺度大小。

5

图5 无湍流时接收机平面上的模拟光斑图 6 无湍流时接收机平面上光场的振幅模拟

(Fig.5 The stimulant facula on the receiver plane(no turbulence), Fig.6 The stimulant amplitude on the

receiver plane(no turbulence))

图7 有湍流时接收机平面上的模拟光斑 图8 有湍流时接收机平面上光场的振幅模拟

(Fig.7 The stimulant facula on the receiver plane(have turbulence), Fig.8 The stimulant amplitude on

the receiver plane(have turbulence))

图5-图8分别模拟给出了光波从处的矩形发射孔径经大气湍流传播到达

z0

zL

处的接收机时，光场的分布变化。为了能做到比较研究，图4也给出了发射机平面上的光

斑分布。从图5，图6可以看出，当传播路径上无大气湍流存在时，接收机平面上的光斑

形状和发射光斑较为相似，能量也是集中分布在光轴上，光斑较为明亮清晰，不同的只

是由于衍射作用的存在，致使光斑在垂直于光轴的平面内发生了一定的衍射扩展，有衍

射旁瓣存在，且旁瓣呈十字形对称分布，离光轴越远亮度越小。其场强的空间分布表现

为准正态的分布形式，光轴上主峰能量分布集中，锐度较大，由于衍射的作用使得能量

的分布呈现为多峰分布。图7，图8是在光传播路径上引入了Kolmogorov大气湍流后的场

分布。模拟表明：接收机平面上的光场同无大气湍流时的图5，图6相比，图像的光学品

质很差，由于大气湍流的存在而导致光斑主瓣在扩张的同时变得较为模糊，亮度下降，

光斑出现了漂移，且与光轴不再表现为严格的十字形对称分布；同时，主瓣的周围出现

了许多大小不一，亮度不同且无规律的旁瓣分布，这是由于大气湍流的存在而导致光斑

出现了破碎。与图6相比，图8中场的主能量仍为准正态分布但幅度减小，能量重心偏移，

主能量峰亦集中了最多的能量，在主峰周围，能量又呈现出强弱不一的多峰空间分布，

但同图6相比多峰的锐度均变小，可见其与单纯受衍射作用的光场截然不同。

[8-13,15,16]

5结束语

当光在湍流大气中传播时，受大气分子、气溶胶等粒子的相互作用而导致发射光束

的光学特性发生改变，从而严重制约了与光在大气湍流中传播有关的许多工程应用的发

展。目前，在中等和强起伏区, 理论和实验都没有能够很好地解决上述问题，这为模拟

研究提供了一定的必要性。文章采用McGlamery算法，对Kolmogorov谱下的大气湍流

6

随机相位屏进行了数值模拟，进而在模拟相位屏的基础上，分析了从方形孔径发射的平

面波经湍流大气传播到达接收机时的光场变化。研究表明，大气湍流的存在严重影响了

光波的远场聚焦特性，对光的传播质量造成很大的影响。模拟分析为大气湍流情况下的

光通信以及自适应光学技术的发展提供了参考。

参考文献：

[1] Ishimaru, Akira, Wave Propagation and Scattering in Random Media[M], Volume 2, Academic Press, New York, 1978.

[2] Andrews, L C., Lar Beam Propagation Through Random Media[M], The Society of Photo-Optical Instrumentation

Engineers, Washington, 1998.

[3] 饶瑞中，激光大气闪烁统计特征的研究进展[J] ,光电子技术与信息，2000，13(5):12-17.(Rao Ruizhong, Rearch progress

on Lar Atmospheric scintillation Statistical Characteristics, Opto-Electronic technique and massage , 2000，13(5):12-17.)

[4] 饶瑞中，光在湍流大气中的传播[M]，合肥，安徽科学技术出版社，2005:24-112.(Rao ruizhong, Light propagation in the

turbulent atmosphere[M],Hefei, Anhui Science and Technology Press. 2005:24-112.)

[5] Benjamin ery. Computer simulation studies of compensation of turbulence degraded images[J]. SPIE，1976，74:

225-233.

[6] Roddier N. Atmospheric wavefront simulation using Zernike polynomials[J]. Opt. Eng., 1990,29:1174-1179.

[7] 王立瑾，李强，魏宏刚等，大气湍流随机相位屏的数值模拟和验证[J]，光电工程，2007,34(3):1-4.(WANG Li-jin, LI Qiang,

WEI Hong-gang, Numerical simulation and validation of pha screen distorted by atmospheric turbulence[J] ,

Opto-Electronic Engineering, 2007,34(3): 1-4.)

[8] Jakobson H, Simulation of time ries of atmospherically distorted waves fronts[J], Appl. Opt.,1996, 35: 1561-1564.

[9] Coles W A, Filice J P, Frehlich R G ,Simulation of wave propagation in three-dimensions media[J]. Appl. Opt.,1995, 34:

2089-2100.

[10] Flatte S M, Martin J M, Intensity images and statistics from numerical simulation of wave propagation in 3-D random

media[J] , Appl. Opt.,1988, 27(11): 2111-2124.

[11] Frehlich R. Simulation of lar beam in a turbulent atmosphere[J], Appl. Opt.,2000,39: 393-397.

[12] Knepp D L. Multiple pha-screen calculation of the temporal behavior of stochastic waves [J], Proc. IEEE,1983,71:722-737

[13] 饶瑞中,王世鹏，刘晓春，被湍流大气退化的激光光斑，尺度测量与形变特征描述[J]，光学学报，1998，18（4）:451-456.(Rao

Ruizhong, Wang Shipeng, Liu Xiaochun, Atmospheric Turbulence-Degraded Light Intensity Images: Size Measurement and

Description of Deformation Character istics[J],ACTC OPTICAL SNICA, 1998，18（4）:451-456.)

[14] 胥杰，赵尚弘，侯睿等，强激光大气传输的远场光束亮度仿真研究[J],光学技术，2008，34(2):233-235.(XU Jie, ZHAO

Shanghong, HOU Rui, Simulation rearch on the far field beam irradiance of high energy lar through atmospheric

propagation, OPTICAL,TECHNIQUE, 2008，34(2):233-235.)

[15] 吴晗平，激光光束质量的评价与应用分析[J]，光学精密工程，2000,8(2):130-132.(WU Hanping, Evaluation and applied

analysis of lar beam quality[J],Optical and precision engineering, 2000,8(2):130-132.)

[16] 杨连臣，沈忙作，郭永洪. 天体目标斑点成像的模拟[J]. 光子学报，2000，29(12):1108-1111.)(YANG Lian-chen，SHEN

Mang-zuo，GUO Yong-hong. The Speckle Imaing Simulation of Space Objects[J]. ACTA PHOTONICA SINICA，2000，

29(12):1108-1111.)

Simulation the Light Beam Propagation Through the

Atmospheric Turbulence

Ma Bao-Ke Guo Li-XinWu Zhen-Sen

1,21 1

(1.School of Science , Xidian University, Xi’an, 710071, China. 2. School of Science , XIan

7

Polytechnic University,Xi’an,710048 ,China. )

Abstract: When light propagates in the atmospheric turbulence, due to interaction of the

atmospheric molecules, aerosols and the other particles, it well experience the atmospheric

turbulence effects of the beam spread, drift, and the coherence degradation and so on, the

factors riously affect the far-field characteristics of the light. Bad on the theory of the light

propagates in the atmospheric turbulence, this article analyzed how to construct a more

reasonable pha screen. And by using the McGlamery algorithm, the atmospheric turbulence

pha screen under the Kolmogorov power spectrum was simulated. And then, the far-field

characteristics of the light which propagates from the transmitter reach to the receiver plane

through the atmospheric turbulence was detailedly analyzed. The studies shows that the

prence of the atmospheric turbulence can bring a great affect on the optical quality, and the

rearch results also provided a reference for the projects which relate on the light propagation

in the atmospheric turbulence and the development of adaptive optics technology.

Keywords: Atmosphere turbulence; McGlamery algorithm; pha screen simulation;

atmosphere turbulent structure constant

作者简介：

马保科(1972-)，男，宁夏人，西安工程大学理学院副教授，西安电子科技大学无线

电物理专业在读博士。Email： baokema2006@

郭立新(1968-)，男，出生于陕西省西安市，现为西安电子科技大学理学院教授，博

士生导师，享受政府特贴的专家。Email : lxguo@

吴振森(1946)，男，出生于湖北沙市，西安电子科技大学理学院教授，博士生导师，

享受政府特贴的专家，美国纽约科学院成员，IEEE高级成员。Email :

wuzhs@

文章主要创新点如下：

1). 对于大气湍流下的光传输研究，在中等强度的起伏区, 目前理论上仍未获得较有价值

的结果，因而本文的模拟研究是必要的。

2). 分析推导了在大气湍流折射率起伏的情况下，空间光场的变化特点。

3). 分析了怎样构造较为合理的大气湍流相位屏。

4). 采用McGlamery算法，结合Huygens-Fresnel原理，模拟生成了Kolmogorov谱下的

大气湍流相位屏，并分析了其有效性尺寸。

5). 利用生成的Kolmogorov大气湍流相位屏，模拟比较了在有无大气湍流两种情况下，

光从发射机孔径经湍流大气传播到达接收机时的光场变化。

6). 模拟结果为大气湍流情况下的光学工程应用及自适应光学技术的完善提供了一定的

参考。

8