第47卷 第9期 2017年9月 激光与红外 LASER & INFRARED Vo1.47,No.9 September,2017 文章编号:1001—5078(2017)09—1076-06 ・激光应用技术・ 相位屏法研究湍流大气对干涉图样的影响 单聪淼,孙华燕,赵延仲,李 帅 (装备学院,北京101416) 摘要:基于马赫一泽德光路发射的相干光为模型,利用功率谱反演法生成随机相位屏的方法 模拟湍流大气,分析了相干高斯光在湍流大气中的传输过程并建立了传输模型,推导得到了接 收屏上的干涉光强。数值仿真分析了湍流大气对干涉图样分布以及条纹间距和对比度等特征 参数的影响。结果表明,随着湍流强度增大,干涉图样变得模糊甚至破碎,条纹间距增大,对比 度下降。此研究结果对于激光主动干涉探测技术的应用具有一定参考意义。 关键词:随机相位屏;马赫一泽德干涉原理;干涉图样;条纹间距;对比度 中图分类号:TN249 文献标识码:A DOI:10.3969/j.issn.1001-5078.2017.09.004 Influence of atmosphere turbulence on interference pattern by phase screen method SHAN Cong-miao,SUN Hua—yan,ZHAO Yan-zhong,LI Shuai (Academy of Equipment Command and Technology,Beijing 101416,China) Abstract:Taking the coherent laser based on Mach—Zehnder light path as the model,random phase screen is generated by the power spectrum method of inversion to simulate the atmospheric turbulence.The transmission of coherent Gauss beams in atmosphere turbulent was analyzed,and the transmission model was established,and the interference light in— tensity on the receiving screen is derived.The influences of atmosphere turbulence on the distribution of interference pattern and the characteristic parameters such as fringe spacing and contrast are analyzed by numerical simulation.The results show that with the increase of turbulence intensity,the interference pattern becomes blurred or even broken,the fringe spacing increases and the contrast ratio decreases.The results of this study have a certain reference value for the application of laser active interference detection technology. Key words:random phase screen;Mach—Zehnder interferometer principle;interference pattern;fringe spacing;contrast ratio 1 引 言 对组合光束反射特性的差异,能够从复杂漫反射背 景中快速识别出光学目标¨ 。然而,湍流大气是影 激光主动干涉探测技术具有抗背景光干扰性 强、回波信息丰富的特点,是激光主动探测技术在未 来的一个重要发展方向。因此,课题组提出了一种 响激光主动干涉探测作用距离的最主要因素,湍流 大气对在其中传输激光的闪烁、漂移、展宽等效应都 会影响激光传输相干性,进而影响目标处的干涉场 分布,因此如何消除湍流大气对激光主动干涉探测 利用相干合成线阵高斯光束扫描识别漫反射背景中 的光学目标的新方法,利用光学目标和漫反射目标 基金项目:国家自然科学基金项目(No.61302183)资助。 作者简介:单聪淼(1988一),女,博士,主要从事激光主动探测,光学信息获取与处理的研究。E—mail:dandanscmiao @163.com 收稿日期:2017-01—19;修订日期:2017-02—19 激光与红外No.9 2017 单聪淼等相位屏法研究湍流大气对干涉图样的影响 1077 的不利影响,可以获得更多的目标特征信息,是值得 研究的课题。 咖( Yi)=∑ ∑ h(x )。  ̄/ (,c , )exp[i(Kx +Kyy )] (1) 式中,h(,c ,,c )为均值为0,方差为1的Hermitian 复高斯随机矩阵; ( ,K )为符合Von Karmon谱 的功率谱函数; 、,c 分别为 方向和Y方向的波矢 分量。其中: 光波在湍流介质中传播的研究方法可以概括为 两类:理论研究的解析方法 和数值模拟研究方 法_6.9 J。解析方法在弱起伏和强起伏时可以通过引 入某种微扰近似获得较好的结果,但对于中等强度 起伏的情况,至今依然没有很好的解析处理方法。 因此,文中采用数值仿真的方法模拟湍流介质中的 相干高斯光束的传输过程,以基于马赫一泽德光路 (,c ,,c )=0.033C (,c +,c:+,J )-1 ・ exp{[一(K:+,c:) f0/5.92】 ) (2) 发射的相干光为模型,利用功率谱反演法生成随机 相位屏的方法模拟湍流大气,分析了相干高斯光在 湍流大气中的传输过程并建立了传输模型,推导得 到了接收屏上的干涉光强。数值仿真分析了湍流大 气对干涉图样分布以及条纹间距和对比度等特征参 数的影响。此研究结果对于激光主动干涉探测技术 的应用具有一定参考意义。 2功率谱反演法生成随机相位屏 利用数值仿真的方法模拟相干高斯光束在湍流 大气中的传输过程,其主要思想是利用数值方法生 成随机相位屏模拟湍流大气,目的是生成对相干光 场分布产生影响的随机相位起伏。现在广泛使用的 功率谱反演法是一种频域间接模拟,这种方法根据 大气湍流的功率谱密度函数得到相位空间随机场, 然后进行傅里叶变换得N--维随机相位的空间分 布,通常也称为快速傅里叶变换(F兀1)法 “ 。如 图1所示,假设激光器发出的光束经过Ⅳ个相位屏 到达距离为 的接收屏上,则屏间距离 = Ⅳ。 图1 光波经过相位屏的传输过程 Fig.1 The transmission of light wave through phase screen 在湍流谱密度函数的选择上,文中采用同时 考虑了湍流内尺度和外尺度、更符合大气折射率 起伏特性的Von Karmon谱模型。具体过程为:在 频率域生成一均值为0、方差为1的复数高斯随机 矩阵,用符合Von Karmon大气相位扰动的功率谱 (,c)对其进行滤波,再进行逆傅里叶变换得到 湍流大气的扰动相位 ( ,Y ): 式中,c 为大气折射率结构常数,当传输路径为水 平传输时,c 可视为常数。 将上式离散化,设 =max,Y =nAy,其中△ 、 △ 为空间域中 方向和Y方向的取样间隔,Ax= / ,ny=R /Ny,m、n为整数,其中R ×R 为相 位屏尺寸, 、Ⅳv分别为 方向和Y方向的采样点个 数; =m AK ,K =n △,c ,其中△,c 、AK 为波数域 中 方向和Y方向的取样间隔,m 、n 为整数,△ = 27#(NxAx),△,( =2订/(NyAy)。则可以得到上式 的离散化结果为: N.Z2一】N/2—1 6(m,n)=’ ∑ ∑一N,/2“ 一Ny/2 h(m ,n )・ xp ( + )】 Nx/2。‘1 N?/2。。1 =, =一N ,2n。=一N ,2 … ,{ (m,,n,) 丽‘ exp i2竹( m + )】 N 2—1 Ny/2—1 =m, =一N√ /2n, :一N /2 圩(m ,n,)exp ( i、 + )i、v 】 NJ2—1 ⅣV2—1 =m , :N 2、n =一N ,2 ,,{, , (m,,n )exp ( ni、v ,)])J ‘ exp[i2'rr(m. )] =FFvr{FFT[日(m ,n )]l ,)J , (3) 式中: 日(m ,n ):h(m ,n ) ̄/ (m , ) (4) 上式即对日(m ,n )的二维FFr变换。 3基于马赫一泽德光路的相干合成高斯光束 由图2可知,一束基模高斯光束经马赫一泽德 干涉光发射装置的分束镜 分为z。和z 两束光, 则两束光在分束镜 上的复振幅分布可以分别 表示为: 1078 激光与红外 第47卷 )= toO唧[_ exp(一ikAS)exp(一ikysin2a) Ut2]・ (5) 2 4相干光经过相位屏的传输过程 相干光经过相位屏传输最后到达接收屏的过程 如图3所示。相干光从初始位置开始在自由空间中 ( 。)= 4-exp[-X 2 "y]z()(6) 传输△ 到达第一个相位屏,经过相位屏时,相位屏 对相干光场的相位扰动作用叠加在相干光场上,再 经过 的传输到达下一个相位屏,此过程依次重 复,直到到达接收屏为止。 “1 ItI U Uf 式中: △Ls=f一4 ̄-csina+ S1nZ + C0S Ot+ HN UN F 一2f (7) _●● ● ( )=too√ + 吉) (8) 其中,△.s为两束光在 上的光程差,Ol为反射镜 的偏转角;z为镜间距离;DG=C,h为 和 的厚 度;to。为高斯光束的束腰半径,即z=0处的 的值。 Fig.2 Mach—Zehnder intefferometer light emitting device principle diagram 则两束光在激光发射端的相干光场可以表 示为: =Uq+UI= 2oto{exp[一 ]・ exp(_i s)ex i in2 )+exp[. ]) (9) 利用式(3)中的方法对上式进行离散化,则相 干光场的离散化表示为: Uo(m, )= 』I r一 p卜——— :± —一二 :].J。 exp(一ikAS)exp(一iknAysin2a)+ 唧[_ 每 ]) (10) Fig.3 Coherent light passing through the phase screen 则由菲涅尔衍射积分公式,可以得到高斯光束 在自由空间传输 到达第一个相位屏的复振幅分 布为: 一 exp(i )exp(i )× ( exp(i )exp[_i2盯( + '7 )I da: …) 对上式进行离散化处理得: )一 exp(i )eXp(i )× N./2-1 NJ2-1 . m—=N 2n √ :一N√2{、 ( xp[i 厶 exp【一i2 c 会 m + — n ,] 一 exp(i )eXp(i )・ Nx/2—1 N√2一l m =一N 2n =一N /2 ∑ ∑s(,n ,n )・ exp-i2,rr( m + n )] =一 exp(i )exp(i ) FFT{FFFES(m',n,)]I ,)l , (12) 式中: S(m ,n )= (m ,n )・ exp【i ] (13) 079 没 i个卡Hf 1I乏示的相f 为 ( , ),激_)匕束 卡_II 入射 fIl川{射 的复振f】{ffi分别记为 ( , 造成接收 的某些点处的光 差 小阿 满址}:涉的条什, 此T涉 佯的肜j 也就发,卜 十H 的变化 , )f『I,, (Y,, ) !J!『J仃: f“ 、 )=“ ( )t x i(f) ( ,, )1 (14) 一0 04 ×、J l 』 进f 离f投化处理i导: “ (7”,,?)=“ (,,?,II)exp【id,,(,7},/1)1 (15) 目—0 02 …欠利川 尔衍射积分公 , 光束 过第i 0 /『、fltf,'rlI ̄ 他输距离 到达第( +I)个相位 ̄i"-l'l"J入射 [f1i_j 的 振幅为: (,, ”,11 ) i 一 r\l )(i )× z ’”一 l,f ,,z k(ma.r) , /! l『iL 一 p 耵( + : ] =一志\『l(i )( i) , Ⅲ 2 + )] 一 i ll(i心 .1- ” ) FF] ;FI一…1[7’(Ⅲ,¨)j【 }} . .(16) 川,,『)_f (『7 m ] (17) 复进行 述过弹,肖刮划达接收屏终L卜没 到达接收J砰 的州I 光场的复振幅分布为 ,则 接收 f 卡H=F场的比 分布为: 1、=(/、 (1 8) 5 数值仿真  ̄ff.tl2值为A=1064 I1I11,, =100 I11,/、’=10, =I,/ =l0Ⅲ, X R =().1『Il×0.1 Ill,,v、= =、1000,, =I I11,/f】=1 lllll1, ()=0.002 111,O/= 0.001 ’.,2=1.5.f =0.01 Il1.『]=0.005 in 、 4为小同湍流强度F,接收胼』二的 涉 样 分伽 }{I图 叮知,随荷湍流的增强,干涉 样 的条纹分 由清晰变得馍糊, I,湍流强度增大到 ・定 ,1 涉 样的整体形念 经破碎,只有儿 条形念 舰删的条纹 山于湍流强度增人,比束 在 过湍流…时的随机起伏也变得剧烈,这就使 得 l 光之问的 程差发,卜小规则的变化. O 02 0.04 .004 .0.02 0 0 02 0.04 x/m (a) 。10’ l1] ? 一O 04 0.02 0 0 O.02 0.O4 ∞ b ( __ O恤“ 、 一—0 04 O 、 —0.02 g 0 002 0 04 -0 04 -0(/2 U 【J112 0 U肆 X/In ㈦(’ m 【 1 4 t: 【 】ft分n J Fig.4 l11 P Ⅲ path … I ihutim ̄ 5为接收 的卡臼对j匕曲 分打i情况 为 便 下lfⅡ}I‘算 涉 样的条纹 1条纹对比瞍. 图 『Il标 光强分 ,j波峰和波 的化嚣,了1:i 求每一 个波峰币¨波谷所刈’心的卡门对光 的值 [1I 还l,JI以 ,卡II对光强分 i曲线与} 涉比慢的-:fli::分 是 十u刈-J 的、 1080 激光与红外 第47卷 躺 靛 如 (a)c暑=10’ m,2 罂 y/m (b)c暑=10mm 蠢 罂 y/m (c) =10‘ om 图5相对光强 Fig.5 Relative intensity 条纹间距是描述干涉图样特征的一个重要参 数,表示相邻两条亮条纹或暗条纹之间的距离。当 不考虑相干光在湍流大气环境中传输时,可以通过 物理光学中对光程差的推导方法,进而推导出接收 屏上干涉图样的条纹间距。但由于湍流大气对在其 中传输的光束具有闪烁、漂移和光束展宽等影响效 应,因此在研究湍流大气中的相干光传输时,已经不 能通过简单的推导给出条纹间距的表达式。这里, 我们从条纹间距的定义出发,通过计算两个相邻的 波峰或者波谷之间的距离得到经过湍流大气后的相 干光形成的干涉图样的条纹间距。 图6为相对光强曲线中波峰和波谷问的峰值点 距离。由图可知,由于湍流大气的影响,干涉图样已 经不再是等问距的,且随着湍流的增强,条纹间距的 变化的随机性越来越明显,对比图6(a)和(b),可 以发现条纹间距有增大的趋势,这是由于湍流大气 对光束的影响在宏观体现上有整体展宽的趋势。当 湍流增大到一定程度时,条纹分布已经失去规律性, 此时的干涉图样已经不再适合主动干涉探测。 O.06 哕符 曰0.059 波峰 掘丑 ‘0058 \ 广一・ 。\ \ / \ / \ 2 4 6 ll肇I值点 (a) =10一 m一 r 0.06 娅 / / ≥ 、 / V 鹫0 059 / \ \ 如盘 磐0.058 { /_/ j 卜、坝峰 0.057 0.056 2 4 6 峰傅点 (b)C ̄=lO n3 爱 量 电; i妥 磐 图6峰值点间距 Fig.6 Peak spacing 条纹对比度是描述干涉图样特征的另外一个重 激光与红外No.9 2017 单聪淼等相位屏法研究湍流大气对干涉图样的影响 1081 要参数,用来表征干涉场中某一点附近的条纹的清 晰程度。在主动干涉探测时,必须保证传输到目标 处的干涉图样具有满足一定清晰度的条纹,以确保 相干探测光经目标反射原路返回到接收端探测器上 时能够获得目标时间化的光强变化信息。条纹对比 度的定义为: K= (19) 式中, 为该点附近条纹强度的极大值;lm为该点 附近条纹强度的极小值。上式表明,条纹对比度是 在0—1之间变化的,其值大小与条纹亮暗差和背景 光强有关。 图7为对应的波峰的峰值光强与其相邻的波谷 之间计算得到的条纹对比度。由图可知,对比度整 体呈现先增大后减小的趋势,这是由于干涉图样的 光强分布整体呈现类高斯的分布,因此,中间波峰点 的对比度会较两边的大。另外,随着湍流强度的增 大,对比度呈现减小趋势。由于湍流大气对光束的 闪烁和漂移效应加剧,使得光束中各点在接收屏上 的落点的随机性变强,光分布的能量变得分散,进而 使得体现光强集中性的对比度减小。 峰值点 图7干涉图样条纹对比度 Fig.7 Fringe contrast of interference pattern 6结论 本文采用数值仿真的方法模拟湍流介质中的相 干高斯光束的传输过程,以基于马赫一泽德光路发 射的相干光为模型,利用功率谱反演法生成随机相 位屏的方法模拟湍流大气,分析了相干高斯光在湍 流大气中的传输过程并建立了传输模型,推导得到 了接收屏上的干涉光强。数值仿真分析了湍流大气 对干涉图样分布以及条纹间距和对比度等特征参数 的影响。此研究结果对于激光主动干涉探测技术的 应用具有一定参考意义。 参考文献: [1] Zhao Yanzhong,Sun Huayan,Shan Congmiao,et a1.A new identification method aimed at optical targets using an active interference laser beam『J].IEEE Photonics Tech— nology Letters,2014,26(10):1019—1022. 『2] Alan M Whitman.Mark J Beran.Two.scale solution for at. mospherie scintillation from a point source[J].J.Opt. Soc.Am.A,1993,5(5):735—737. [3] Beran M J,Whitman A M.Scintillation index calculation using an altitude dependent structure constant[J].App1. Opt,1988,27(11),2178—2182. [4] JWR Brown.rough moment of a wave propagating in a random medium[J].J.Opt.Soc.Am.A,1972,62(8): 966—971. [5] Andrew Forbes.Laser beam propagation in non—linearly absorbing media[J].SPIE,2006,6290:629003. [6]YU Jiping,QI Wenzong,GUO Chunfeng,et a1.Numerical simulation of laser propagation through atmospheric turbu— lence[J].Laser&Infrared,2008,38(6):523—527.(in Chinese) 于继平,齐文宗,郭春凤,等.激光大气传输特性的数 值模拟[J].激光与红外,2008,38(6):523—527. [7]LI Bo,WANG Tingfeng,WANG Dinan,et a1.Simulation of laser beam propagation through turbulence[J].Chinese Optics,2012,5(3):289—295.(in Chinese) 李波,王挺峰,王弟男,等,激光大气传输湍流扰动仿 真技术[J].中国光学,2012,5(3):289—295. [8]XU Jiawen,SUN Huayan,WU Weiwei,et a1.Research of laser propagation simulation in turbulent atmosphere[J]. Equipment Manufaetring Technology.2009,4:43—45. (in Chinese) 许嘉纹,孙华燕,吴伟伟,等.湍流大气中的激光传输 仿真研究[J].装备制造技术,2009,4:43—45. [9]XU Guangyong,wu Jian,YANG Chunping,et 1a.Simula— tion and optical scintillation research of Gaussian beam in atmosphere turbulence[J].Laser Technology,2008,32 (5):548—550.(in Chinese) 徐光勇,吴健,杨春平,等.高斯光束在大气湍流中的 数值模拟和光强起伏[J].激光技术,2008,32(5): 548—550. [10]WANG Lijin,LI Qiang,WEI Honggang,et a1.Numerical simulation and validation of phase screen distorted by at・ mospheric turbulence[J].Opto—Electronic Engineering, 2007,34(3):1—5.(in Chinese) 王立瑾,李强,魏宏刚,等.大气湍流随机相位屏的数 值模拟和验证[J].光电工程,2007,34(3):1—5. [11]WANG Long,SHEN Xueju,ZHANG Weian,et a1.Study on propagation of flattened Gaussian beams through turbu— lent atmosphere utilizing multiple phase—screens[J].IJa- ser&Infrared,2012,42(8):852—856.(in Chinese) 王龙,沈学举,张维安,等.相位屏法研究平顶高斯光 束的大气传输特性[J].激光与红外,2012,42(8): 852—856.