留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

湍流变化对多孔径光学系统成像特性的影响

刘婷 唐善发 刘何伟 钱俊宏 张蓉竹

刘婷, 唐善发, 刘何伟, 钱俊宏, 张蓉竹. 湍流变化对多孔径光学系统成像特性的影响[J]. 红外与激光工程, 2021, 50(12): 20210189. doi: 10.3788/IRLA20210189
引用本文: 刘婷, 唐善发, 刘何伟, 钱俊宏, 张蓉竹. 湍流变化对多孔径光学系统成像特性的影响[J]. 红外与激光工程, 2021, 50(12): 20210189. doi: 10.3788/IRLA20210189
Liu Ting, Tang Shanfa, Liu Hewei, Qian Junhong, Zhang Rongzhu. Influence of turbulence on imaging characteristics of multi-aperture optical system[J]. Infrared and Laser Engineering, 2021, 50(12): 20210189. doi: 10.3788/IRLA20210189
Citation: Liu Ting, Tang Shanfa, Liu Hewei, Qian Junhong, Zhang Rongzhu. Influence of turbulence on imaging characteristics of multi-aperture optical system[J]. Infrared and Laser Engineering, 2021, 50(12): 20210189. doi: 10.3788/IRLA20210189

湍流变化对多孔径光学系统成像特性的影响

doi: 10.3788/IRLA20210189
基金项目: 四川省重大科学仪器设备专项(2019ZDZX0038)
详细信息
    作者简介:

    刘婷,女,硕士生,主要从事光学设计与成像理论方面的研究

  • 中图分类号: TH751

Influence of turbulence on imaging characteristics of multi-aperture optical system

  • 摘要: 为了研究大气湍流变化对多孔径光学系统成像质量的影响,针对Golay3结构建立了一个理论分析模型,推导了湍流影响下望远系统的点扩散函数(point spread function, PSF)的表达式。分别针对近似圆形和近似长条状分布的湍流结构对成像特性的影响进行了讨论,具体构建了两种湍流影响下的光场相位计算模型,并计算比较了两种湍流的调制传递函数(modulation transfer function, MTF)。结果表明,光路中存在湍流将导致系统成像质量下降。在近似圆形湍流中,湍流强度越小,系统MTF影响越小。不同强度湍流影响的MTF在归一化空间频率0.16、0.45和0.69处均下降0.05左右。在近似长条状湍流中,对流湍流风速越小,系统MTF影响越小。不同风速的湍流影响的MTF在空间频率0.16、0.42和0.69处分别下降0.25、0.09和0.05左右。比较两种湍流表明,近似长条状湍流对成像系统MTF影响更为明显。
  • 图  1  Golay3系统阵列结构

    Figure  1.  Golay3 system array structure

    图  2  MTF分布((a)无误差;(b)活塞误差π;(c)倾斜误差0.012 μrad)和面物体成像((d)无误差;(e)活塞误差π;(f)倾斜误差0.012 μrad)

    Figure  2.  MTF distribution( (a) Error free; (b) Piston error π; (c) Tilt error 0.012 μrad) and surface object imaging ((d) Error free; (e) Piston error π; (f) Tilt error 0.012 μrad)

    图  3  (a)近似圆形湍流在视场中的示意图及(b)湍流区域相位分布

    Figure  3.  (a) Schematic diagram of similar circular turbulence in the field of view and (b) phase distribution in turbulent region

    图  4  (a)近似长条状湍流在视场中的示意图;(b)湍流气体密度分布及(c)相位分布

    Figure  4.  (a) Schematic diagram of similar elongated turbulence in the field of view; (b) Turbulent gas density distribution and (c) phase distribution

    图  5  近似圆形湍流对系统的MTF影响(相位差(a) π/2 (3π/2);(b) π)及面物体成像(相位差(c) π/2 (3π/2);(d) π;(e)不同相位差MTF曲线)

    Figure  5.  Influence of similar circular turbulence on MTF of system ((a) π/2 (3π/2); (b) π) and plane object imaging phase difference ((c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences)

    图  6  不同高度湍流影响的MTF ((a) 8 km; (b) 10 km;(c) 15 km; (d) 20 km)

    Figure  6.  MTF affected by turbulence at different heights ((a) 8 km; (b) 10 km; (c) 15 km; (d) 20 km)

    图  7  不同强度湍流的MTF ((a) 10−8 m3−α ;(b) 10−13 m3−α;(c) 10−15 m3−α;)及(d) 不同湍流强度的MTF曲线

    Figure  7.  MTF of different intensities of turbulence ((a) 10−8 m3−α; (b) 10−13 m3−α; (c) 10−15 m3−α) and (d) MTF curves of different turbulence intensities

    图  8  近似长条状湍流对系统MTF影响(相位差(a) π/2 (3π/2);(b) π面物体成像相位差;(c) π/2 (3π/2);(d) π;(e)不同相位差的MTF曲线

    Figure  8.  Influence of similar elongated turbulence on MTF of system ((a) π/2 (3π/2); (b) π plane object imaging phase difference; (c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences

    图  9  不同风速影响的MTF ((a) 10 m/s;(b) 20 m/s;(c) 30 m/s)及(d)不同风速影响的MTF曲线

    Figure  9.  MTF affected by different wind speeds ((a) 10 m/s; (b) 20 m/s; (c) 30 m/s) and (d) MTF curves of different wind speed influence

  • [1] Golay M J E. Point arrays having compact, nonredundant autocorrelations [J]. JOSA, 1971, 61(2): 272-273. doi:  10.1364/JOSA.61.000272
    [2] Lampkin C M, Flint G W, MacFarlane M J. Optical and mechanical design aspects of a four-telescope array for combined imaging [J]. Opt Eng, 1988, 27(9): 749-754.
    [3] Christopher R, De Hainaut K P, Hentz L D, et al. Design of a wide field of view phased array telescope [J]. Opt Eng, 1988, 27(9): 736-739.
    [4] Flores Jorge L, Strojnik Marija, Páez Gonzalo, et al. Effects of misalignment errors on the optical transfer functions of synthetic aperture telescopes [J]. Applied Optics, 2004, 43(32): 5926-5932. doi:  10.1364/AO.43.005926
    [5] Kendrick R L, Aubrun J N, Bell R, et al. Wide-field fizeau imaging telescope: Experimental results [J]. Appl Opt, 2006, 45(18): 4235-4240. doi:  10.1364/AO.45.004235
    [6] Ni M, Benson L, Camp J, et al. Autonomous tip/tilt alignment and phasing of adistributed aperture imaging testbed [J]. Opt Express, 2010, 18(12): 13051-13056. doi:  10.1364/OE.18.013051
    [7] 吴泉英. 稀疏孔径光学系统成像研究[D]. 苏州大学, 2006.

    Wu Quanying. Research on imaging of sparse aperture optical system[D]. Suzhou: Soochow University, 2006. (in Chinese)
    [8] Zhang Wei, Deng Jian, Long Funian, et al. Research on image quality evaluation of sub-aperture synthetic optical imaging system [J]. Acta Optica Sinica, 2008, 28(4): 687-691. (in Chinese) doi:  10.3321/j.issn:0253-2239.2008.04.014
    [9] Sun Jiangqin, Qian Lin, Wu Quanying. Modulation transfer function and imaging properties of axisymmetric synthetic aperture optical system [J]. Applied Optics, 2008, 29(4): 548-552, 564. (in Chinese)
    [10] Xie Zongliang, Ma Haotong, Li Bo. Experimental demonstration of enhanced resolution of a Golay3 sparse-aperture telescope [J]. Chinese Optics Letters, 2017, 15(4): 34-37.
    [11] Fiete R D, Tantalo T A, Calus J R, et al. Image quality of sparse-aperture designs for remote sensing [J]. Opt Eng, 2002, 41(8): 1957-1969. doi:  10.1117/1.1490555
    [12] Tang Shanfa, Yan Jiawei, Yong Kangle, et al. Propagation characteristics of vortex beams in anisotropic atmospheric turbulence [J]. Opt Soc Am, 2020, 37(1): 133-137. doi:  10.1364/JOSAB.37.000133
    [13] Zhang Huimin, Li Xinyang. Research on numerical simulation method of atmospheric turbulence distortion phase screen [J]. Optoelectronic Engineering, 2006, 33(1): 14-19. (in Chinese) doi:  10.3969/j.issn.1003-501X.2006.01.004
    [14] 张滔. 大气湍流中涡旋光束轨道角动量谱的扩散及修正的研究[D]. 电子科技大学, 2017.

    Zhang Tao. Diffusion and correction of orbital angular momentum spectrum of vortex beams in atmospheric turbulence[D]. Chengdu: University of Electronic Science and Technology of China, 2017. (in Chinese)
    [15] Victor Klimenko, Jimmy Krozel. Clear-air turbulence impact modeling based on flight route analysis [C]//AIAA Guidance, Navigation, and Control Conference, 2011, 6: 5199-5211.
    [16] 范源丹. 基于垂直载荷因子的晴空湍流检测方法研究[D]. 中国民航大学, 2019.

    Fan Yuandan. Study on the clear-air turbulence detection method based on vertical load factor[D]. Tianjin: Civil Aviation University of China, 2019. (in Chinese)
    [17] 曹原. 喷流冷却条件下高速飞行器光学窗口气动特性研究[D]. 哈尔滨工业大学, 2020.

    Cao Yuan. Aerodynamic characteristics of high speed aircraft optical window under jet cooling condition[D]. Harbin: Harbin Institute of Technology, 2020. (in Chinese)
  • [1] 吴绍华, 黄攀, 赵劲松, 赵跃进, 郑丽和, 董汝昆.  三温区梯度化学气相沉积ZnS制备 . 红外与激光工程, 2021, 50(10): 20210004-1-20210004-6. doi: 10.3788/IRLA20210004
    [2] 张爱武, 张希珍, 赵江华.  斜模超分辨率成像最佳角度 . 红外与激光工程, 2019, 48(8): 826001-0826001(7). doi: 10.3788/IRLA201948.0826001
    [3] 马昊军, 王国林, 陈德江, 张军, 刘丽萍, 罗杰.  热环境条件下红外窗口气动光学传输效应实验研究 . 红外与激光工程, 2017, 46(9): 904001-0904001(7). doi: 10.3788/IRLA201746.0904001
    [4] 袁影, 王晓蕊, 吴雄雄, 穆江浩, 张艳.  多孔径压缩编码超分辨率大视场成像方法 . 红外与激光工程, 2017, 46(8): 824001-0824001(7). doi: 10.3788/IRLA201746.0824001
    [5] 丁浩林, 易仕和, 付佳, 朱杨柱, 何霖.  超声速湍流边界层气动光学效应的实验研究 . 红外与激光工程, 2016, 45(10): 1018007-1018007(7). doi: 10.3788/IRLA201645.1018007
    [6] 施展, 樊祥, 程正东, 朱斌, 陈熠.  关联成像的点扩散函数分析法 . 红外与激光工程, 2016, 45(11): 1124001-1124001(6). doi: 10.3788/IRLA201645.1124001
    [7] 郭小虎, 赵跃进, 吴益剑.  波前编码技术在同轴三反系统的应用及其分析 . 红外与激光工程, 2015, 44(7): 2075-2079.
    [8] 张发强, 樊祥, 朱斌, 程正东, 方义强.  折衍混合长波红外光学系统消热差设计 . 红外与激光工程, 2015, 44(4): 1158-1163.
    [9] 高震宇, 方伟, 张浩, 杨振岭, 王玉鹏.  小型宽光谱太阳光谱仪光学设计 . 红外与激光工程, 2015, 44(2): 590-595.
    [10] 胡玲, 王霞, 延波, 李帅帅.  水下距离选通成像系统调制传递函数模型分析 . 红外与激光工程, 2015, 44(11): 3262-3269.
    [11] 张傲, 汪清, 杨敬钰, 孙懿.  编码掩模红外成像的建模与性能分析 . 红外与激光工程, 2015, 44(10): 2891-2899.
    [12] 贺小军, 张贵祥, 郑亮亮, 曲宏松.  高分辨相机传函测试误差分析 . 红外与激光工程, 2014, 43(S1): 89-94.
    [13] 武星星, 刘金国, 孔德柱, 张立平, 龙科慧.  基于地球椭球的离轴式双线阵相机像移补偿 . 红外与激光工程, 2014, 43(3): 838-844.
    [14] 范哲源, 高立民, 张志, 陈卫宁, 杨洪涛, 张建, 武力, 曹剑中.  中波红外三视场变焦光学系统设计 . 红外与激光工程, 2014, 43(2): 523-527.
    [15] 陈宇, 宋玉龙, 霍富荣.  用于导弹逼近告警的“日盲”紫外光学系统设计 . 红外与激光工程, 2014, 43(9): 2964-2969.
    [16] 孙鸣捷, 于康龙, 孟照魁.  像素感光面形状对基于微扫描的红外超分辨成像的影响 . 红外与激光工程, 2014, 43(10): 3183-3187.
    [17] 任秉文, 金光, 张元, 钟兴, 孔林.  畸变对TDI成像相机的像移影响研究 . 红外与激光工程, 2014, 43(12): 3951-3957.
    [18] 刘可辉, 王晓蕊, 张卫国.  红外点目标成像信噪比建模分析 . 红外与激光工程, 2014, 43(7): 2143-2147.
    [19] 孟伟, 金龙旭, 李国宁, 傅瑶.  调制传递函数在遥感图像复原中的应用 . 红外与激光工程, 2014, 43(5): 1690-1696.
    [20] 葛琪, 王可东, 张弘, 李桂斌, 邸超.  长曝光大气湍流退化图像点扩散函数估计 . 红外与激光工程, 2014, 43(4): 1327-1331.
  • 加载中
图(9)
计量
  • 文章访问数:  38
  • HTML全文浏览量:  24
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-19
  • 修回日期:  2021-05-20
  • 刊出日期:  2021-12-31

湍流变化对多孔径光学系统成像特性的影响

doi: 10.3788/IRLA20210189
    作者简介:

    刘婷,女,硕士生,主要从事光学设计与成像理论方面的研究

基金项目:  四川省重大科学仪器设备专项(2019ZDZX0038)
  • 中图分类号: TH751

摘要: 为了研究大气湍流变化对多孔径光学系统成像质量的影响,针对Golay3结构建立了一个理论分析模型,推导了湍流影响下望远系统的点扩散函数(point spread function, PSF)的表达式。分别针对近似圆形和近似长条状分布的湍流结构对成像特性的影响进行了讨论,具体构建了两种湍流影响下的光场相位计算模型,并计算比较了两种湍流的调制传递函数(modulation transfer function, MTF)。结果表明,光路中存在湍流将导致系统成像质量下降。在近似圆形湍流中,湍流强度越小,系统MTF影响越小。不同强度湍流影响的MTF在归一化空间频率0.16、0.45和0.69处均下降0.05左右。在近似长条状湍流中,对流湍流风速越小,系统MTF影响越小。不同风速的湍流影响的MTF在空间频率0.16、0.42和0.69处分别下降0.25、0.09和0.05左右。比较两种湍流表明,近似长条状湍流对成像系统MTF影响更为明显。

English Abstract

    • 为了提高望远系统的空间分辨率,增大光学孔径是一个有效的手段。然而保证光学系统的成像质量需要高质量的光学元件,大口径光学元件的加工难度大、周期长成为望远系统口径进一步扩大的瓶颈。采用离散多孔径结构,借助多孔径技术在扩大总视场的同时采用多个小口径元件控制成本,是提高系统分辨率的一个有效途径。

      为了获得清晰的成像结果,多孔径望远系统需要一个合理的结构设计,以排除结构误差对成像质量的影响。Golay于1971年提出了Golay6和Golay9结构的多孔径阵列结构[1],奠定了多孔径成像系统的基本研究模型。1988年,美国空军实验室建造了多用途多望远镜试验台(Multipurpose Multiple Telescope Testbed, MMTT)[2],在此基础上,C. R. De Hainaut 等人对多孔径系统共相误差检测和校正进行了初步探讨[3]。2004年,J. L. Flores等人具体研究了多孔径望远系统像差对光学系统MTF的影响[4]。2006年,Lockhedd-Martin公司研制了STAR-9多孔径系统[5],为了解决该系统的共相问题,针对活塞误差提出了利用MTF次峰变化结合点扩散函数对称性的检测方法[6]。2006年,苏州大学吴泉英等人研究了多孔径阵列结构对系统成像性能的影响,并提出了一种新型稀疏孔径结构 [7]。2008年,哈尔滨工业大学张伟等人利用基于三子镜合成光学成像系统干涉仪实测参量拟合的泽尼克多项式研究了像差对多孔径成像系统的影响[8-9]。2017年,中国科学院光电技术研究所谢宗良等人搭建了三孔成像系统实验平台[10],基于该实验平台分析了共相误差对多孔径成像系统成像性能的影响,并成功利用SPGD算法校正了该系统的共相误差。

      目前国内外对影响多孔径光学系统成像质量的研究大多针对孔径结构及孔径间的共相误差分析。在实际应用中,探测光路上的大气湍流变化会带来相位不确定性,从而影响光学系统的成像质量。对多孔径系统而言,每个子孔径都会因为湍流产生相位畸变导致像质降低,同时孔径间的像面叠加使得湍流的影响远比单孔径成像系统复杂。为了研究多孔径结构下湍流对光学系统的影响,文中以Golay3结构为研究对象,针对大气湍流的影响进行了分析。

    • Golay3阵列结构如图1所示。

      图  1  Golay3系统阵列结构

      Figure 1.  Golay3 system array structure

      Golay3光学系统的点扩散函数(point spread function, PSF)为:

      $$ \begin{split} PSF( {x,y} ){\rm{ = }}&PS{F_{{\rm{one}}}}( {x,y} ) \cdot \Bigg\{ {3 + } \Bigg. 2{\rm{cos}}\Bigg[ \frac{{2\pi }}{{\lambda f}}\Bigg( \frac{\sqrt 3 ( {D - d} )}{2}x + \\ & \frac{{D - d}}{2}y \Bigg) \Bigg] + {2{\rm{cos}}\Bigg[ {\frac{{2\pi }}{{\lambda f}}\Bigg( { \pm \frac{{\sqrt 3 ( {D - d} )}}{{\rm{4}}}x + \frac{{{\rm{3}}( {D - d} )}}{{\rm{4}}}y} \Bigg)} \Bigg]} \Bigg\} \\ \end{split} $$ (1)

      式中:$PS{F_{{\rm{one}}}}\left( {x,y} \right)$为单孔径PSF;$d$为子孔径直径;$D$为外接圆直径;$\lambda $为入射光波长;$f$为系统焦距。

      Golay3光学系统的调制传递函数(modulation transfer function, MTF) [11]为:

      $$ \begin{split} MTF\left( {{f_x},{f_y}} \right){\rm{ = }}&MT{F_{{\rm{one}}}}\left( {{f_x},{f_y}} \right) * \left[ {\delta \left( {{f_x},{f_y}} \right) + } \right. \\ &\dfrac{1}{3}\delta \left( {{f_x} \pm \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{2\lambda f}},{f_y} \pm \dfrac{{D - d}}{{2\lambda f}}} \right) + \\ &\dfrac{1}{3}\delta \left( {{f_x} \pm \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{4\lambda f}},{f_y} \pm \dfrac{{3\left( {D - d} \right)}}{{4\lambda f}}} \right) + \\ &\left. {\dfrac{1}{3}\delta \left( {{f_x} \mp \dfrac{{\sqrt 3 \left( {D - d} \right)}}{{4\lambda f}},{f_y} \pm \dfrac{{3\left( {D - d} \right)}}{{4\lambda f}}} \right)} \right] \\ \end{split} $$ (2)

      式中:$MT{F_{{\rm{one}}}}\left( {{f_x},{f_y}} \right)$表示单孔径MTF;$({f_x},{f_y})$表示空间频率。

      由公式(2)可知,当子孔径间存在相位误差时,孔径间的相对位置关系不再对称,导致系统的MTF发生改变,成像质量下降。

      子孔径间的相位误差对光学系统成像质量影响非常大,在实验调试阶段必须对其进行严格校正。图2给出了Golay3结构中子孔径3存在活塞误差或倾斜误差时的MTF分布图以及相应的面物体成像。

      图  2  MTF分布((a)无误差;(b)活塞误差π;(c)倾斜误差0.012 μrad)和面物体成像((d)无误差;(e)活塞误差π;(f)倾斜误差0.012 μrad)

      Figure 2.  MTF distribution( (a) Error free; (b) Piston error π; (c) Tilt error 0.012 μrad) and surface object imaging ((d) Error free; (e) Piston error π; (f) Tilt error 0.012 μrad)

      图2可以看出,无误差的多孔径MTF分布呈中心对称分布,中心能量最大,周围对称分布六个旁瓣,这与子孔径的排列方式有关。当最大活塞误差为π时,系统MTF偏移量最大。倾斜误差不改变MTF中心位置,但使得MTF分布地更加杂乱无章。多孔径子光束间存在活塞误差和倾斜误差都会导致成像结果模糊。所以,使用多孔径光学系统观测前必须对子孔径进行调整,以排除系统自身的活塞误差和倾斜误差。

    • 湍流是大气的一种运动形式,是一种复杂的现象,在湍流大气中,折射率随时间、空间随机变化。大多数情况下,湍流会充满整个视场空间,但由于大气环境复杂,湍流的尺度变化很大,因此还会存在仅影响视场部分区域的小尺度湍流。局部小尺度湍流虽然没有充满整个视场,但是仍会对目标观察产生干扰。文中重点针对这种外尺度较小的湍流情况进行分析。

      实际观测中若多孔径望远系统全视场角为0.2°,观测距离为10 km目标,则子孔径可观测到的视场范围约为350 m。当成像光路中出现近似为团状湍流,湍流外尺度小于子孔径的视场时,则望远系统视场中有一部分近似圆形的区域入射光相位值与其他部分不同发生了畸变。为了分析近似圆形湍流对成像系统的影响,对其进行了简化模型构建,如图3(a)中阴影部分所示。

      图  3  (a)近似圆形湍流在视场中的示意图及(b)湍流区域相位分布

      Figure 3.  (a) Schematic diagram of similar circular turbulence in the field of view and (b) phase distribution in turbulent region

      近似圆形湍流引起的相位差由以下方法计算。根据参考文献[12]使用各向异性大气湍流功率谱${\varPhi _n}({\kappa _x},{\kappa _y})$

      $${\varPhi _n}({\kappa _x},{\kappa _y}) = \dfrac{{A(\alpha )\overline C_n^2{\mu _x}{\mu _y}}}{{{{(\mu _x^2\kappa _x^2 + \mu _y^2\kappa _y^2)}^{{\alpha / 2}}}}}$$ (3)

      式中:${\mu _x}$${\mu _y}$分别表示$x$$y$方向各向异性分子;$({\kappa _x},{\kappa _y})$表示空间频率;$\bar C_n^2$表示折射率广义结构常数,单位为${\rm{m} ^{3 - \alpha }}$,是表示湍流强度的一个重要参数;$A(\alpha )$表示保持其指数结构函数与其功率谱一致性的常数。

      $$A(\alpha ) = \frac{1}{{4{\pi ^2}}}\varGamma (\alpha - 1)\cos ({{\alpha \pi } / 2})\quad (3 < \alpha < 4)$$ (4)

      式中:$\varGamma (x)$表示$\gamma $函数;$\alpha $表示与高度有关的光谱指数。

      通过对大气湍流的功率谱进行滤波,再做逆傅里叶变换可以计算出湍流影响的相位值[13-14]分布如图3(b)所示。

      随着湍流随机运动,如果空间存在一个分布方向正好垂直于视场的近似长条状湍流,例如对流湍流中山波湍流[15-16],望远系统视场中有一部分近似长条状的区域入射光相位值与其他部分不同,发生了畸变,对其进行简化模型构建,如图4(a)阴影部分所示。

      图  4  (a)近似长条状湍流在视场中的示意图;(b)湍流气体密度分布及(c)相位分布

      Figure 4.  (a) Schematic diagram of similar elongated turbulence in the field of view; (b) Turbulent gas density distribution and (c) phase distribution

      文中近似长条状湍流引起的相位差由以下方法计算。对流湍流的分析模型需要结合空气的流体动力学情况来建立,利用Fluent软件计算强风通过山脉产生的山波湍流,得到湍流区域的密度分布,如图4(b)所示,再利用Gladstone-Dale关系:

      $$n = 1 + K\rho $$ (5)

      式中:$n$为折射率;$\;\rho $为密度;$K$为Gladstone-Dale系数[17]

      $$K = 2.23 \times {10^{ - 4}}\left( {1 + \frac{{7.52 \times {{10}^{ - 3}}}}{{{\lambda ^2}}}} \right)$$ (6)

      式中:$\lambda $为入射光波长。

      计算出湍流区域的折射率场分布,进而可求出光场经过近似长条状湍流的相位分布,如图4(c)所示。

    • 近似圆形湍流影响下的系统PSF可表示为:

      $$ \begin{split} PSF\left( {x,y} \right) =& \left| {\sum\limits_{n = 1}^3 {\left\{ {\left( {\frac{{\pi {{{d}}^2}}}{{4\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi {{d}}{r_n}} / {\lambda f}}} \right)}}{{{{\pi {{d}}{r_n}} / {\lambda f}}}}} \right]} \right.} } \right. + \\ &{\left. {\left( {\frac{{\pi {l^2}}}{{4\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi l{r_{ln}}} / {\lambda f}}} \right)}}{{{{\pi l{r_{ln}}} / {\lambda f}}}}} \right] \cdot \left( {{{\rm{e}}^{i\varphi }} - 1} \right)} \right|^2} \\ \end{split} $$ (7)

      式中:${J_1}$表示一阶贝塞尔函数;${r_n} = \sqrt {x_n^2 + y_n^2} $$({x_n},{y_n})$表示子孔径中心坐标;$d$表示子孔径直径;${r_{ln}} = \sqrt {x_{ln}^2 + y_{ln}^2} $$({x_{{l_n}}},{y_l}_{_n})$表示近似圆形湍流中心坐标;$\varphi $表示湍流影响附加相位差。后文中湍流影响导致的附加相位差都用相位差表示。公式(7)等号右边绝对值内第一、二项分别表示没有湍流影响和有湍流影响的振幅扩散函数,振幅扩散函数模平方为PSF。

      湍流变化引起的大气折射率起伏导致光波波前畸变,改变望远镜系统视场内相位均匀性,造成光学图像的模糊,影响光学系统PSF,进而恶化系统MTF。

      系统的MTF为:

      $$MTF\left( {{f_x},{f_y}} \right) = \int {\int_{ - \infty }^\infty {PSF\left( {x,y} \right)} } {{\rm{e}}^{ - i2\pi \left( {{f_x}x + {f_y}y} \right)}}{\rm{d}}x{\rm{d}}y$$ (8)

      通过改变近似圆形湍流在视场中的不同位置,计算得出湍流位于子孔径边缘对系统成像性能影响最大。并且经过计算发现,当湍流区域直径小于d/6时,系统MTF分布基本不变,即尺度较小的湍流对系统成像影响不大。

      经过计算发现,湍流在视场中的位置和尺度大小会影响系统MTF,但其主要的影响规律相同。如图3(a)所示,近似圆形湍流直径为d/3,湍流尺度约为100 m,取湍流层的传输距离为512 m,分析湍流变化对系统成像的影响。

      图5给出了近似圆形湍流区域相位变化$0 \sim 2\pi $的MTF分布以及相应的面物体成像和不同相位差的MTF曲线。

      图  5  近似圆形湍流对系统的MTF影响(相位差(a) π/2 (3π/2);(b) π)及面物体成像(相位差(c) π/2 (3π/2);(d) π;(e)不同相位差MTF曲线)

      Figure 5.  Influence of similar circular turbulence on MTF of system ((a) π/2 (3π/2); (b) π) and plane object imaging phase difference ((c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences)

      图5可知,当近似圆形湍流影响改变,系统MTF中心位置并没有发生移动,随着相位差由0增加到π,MTF畸变程度越来越严重,中心由圆形变成椭圆,椭圆两尖端逐渐分裂开,周围的六个旁瓣范围逐渐缩小。从图5(e)定量分析,以相位差为$0\pi (2\pi )$的MTF曲线为基准,相位差为${\pi / 2}({{3\pi } / 2})$$\pi $的MTF分别在归一化空间频率0.16处下降了0.06和0.19,在频率0.45处下降了0.04和0.09,在频率0.69下降了0.034和0.08。

      另外,当湍流强度为$\overline C_n^2 = {10^{ - 8}}\;{{\rm{m}}^{3 - \alpha }}$,距离入瞳距离分别为h=8、10、15、20 km时,由参考文献[12]可知当湍流强度一定时,湍流导致的相位差值与湍流层高度有关,所以在近似圆形湍流区域构建不同高度的随机相位屏,相位分布的值放入512×512的矩阵中,矩阵中每点的相位差值可表示对应空间位置的相位起伏,以此来研究近似圆形湍流影响下的MTF。便于分析计算出对应随机相位屏的相位均值分别为70.82 π、8.09 π、0.61 π、0.027 π。

      图6可知,当高度为8、10、15 km时,MTF的畸变比较明显,中央类似一个椭圆形,周围六个旁瓣位置并没有发生明显变化。当高度为20 km时,圆形湍流对MTF分布的影响不大。

      图  6  不同高度湍流影响的MTF ((a) 8 km; (b) 10 km;(c) 15 km; (d) 20 km)

      Figure 6.  MTF affected by turbulence at different heights ((a) 8 km; (b) 10 km; (c) 15 km; (d) 20 km)

      当湍流高度为8 km时改变湍流强度,令$ {\overline{C}}_{n}^{2}= {\rm{10}}^{\rm{-8}}\rm{ }{\text{、}} $$ {\rm{ 10}}^{\rm{-13}}\rm{ }{\text{、}}{\rm{ 10}}^{\rm{-15}}\;{\rm{m}}^{3-\alpha }$,在近似圆形湍流区域构建不同湍流强度的随机相位屏,并计算出对应随机相位屏的相位均值分别为70.82 π、1.11 π、0.14 π。

      图7(d)所示MTF曲线为过MTF分布中心且平行于${f_x}$方向的数据,结合随机湍流相位屏均值,可知同一高度湍流强度越小对系统的MTF影响越小。以理想情况的MTF曲线为基准,湍流强度最小的MTF曲线和理想情况相似,而另外两条强度大的MTF曲线,在归一化空间频率0.16、0.45和0.69处都有所下降。强度为${\rm{1}}{{\rm{0}}^{{\rm{ - 13}}}}\;{{\rm{m}}^{3 - \alpha }}$${\rm{1}}{{\rm{0}}^{{\rm{ - 8}}}}\;{{\rm{m}}^{{\rm{3 - }}\alpha }}$MTF曲线分别在频率0.16处下降0.06和0.09,0.45处下降0.04和0.05,0.69处下降0.03和0.04。

      图  7  不同强度湍流的MTF ((a) 10−8 m3−α ;(b) 10−13 m3−α;(c) 10−15 m3−α;)及(d) 不同湍流强度的MTF曲线

      Figure 7.  MTF of different intensities of turbulence ((a) 10−8 m3−α; (b) 10−13 m3−α; (c) 10−15 m3−α) and (d) MTF curves of different turbulence intensities

    • 近似长条状湍流影响下的系统PSF可表示为:

      $$ \begin{split} PSF\left( {x,y} \right) =& \left| {\sum\limits_{n = 1}^3 {\left\{ {\left( {\frac{{\pi {d^2}}}{{\lambda f}}} \right)\left[ {\frac{{2{J_1}\left( {{{\pi d{r_n}} / {\lambda f}}} \right)}}{{{{\pi d{r_n}} / {\lambda f}}}}} \right]} \right.} } \right. + b\sqrt {{d^2} - {b^2}} {\rm{sin}}c\cdot \\ & {\left. {\left( {\frac{{x - {x_{ln}}}}{{\lambda f}}b,\frac{{y - {y_{ln}}}}{{\lambda f}}\sqrt {{d^2} - {b^2}} } \right) \cdot \left( {{{\rm{e}}^{i\varphi }} - 1} \right)} \right|^2} \\ \end{split} $$ (9)

      式中:$({x_{{l_n}}},{y_l}_{_n})$表示近似长条状湍流中心坐标;$b$表示近似长条状湍流宽度。公式(9)等号右边绝对值内第一、二项分别表示没有湍流影响和有湍流影响的振幅扩散函数。

      通过改变近似长条状湍流在视场中的位置,计算得出位于子孔径中心的湍流,能更好地突出对光学系统性能的影响。并且计算发现,湍流区域宽度越窄,对系统MTF分布变化影响越大。

      图4(a)所示,近似长条状湍流宽度为${d / {10}}$,取湍流层的传输距离为512 m,分析湍流变化对系统成像的影响。图8给出了近似长条状湍流区域相位变化$0 \sim 2\pi $的MTF分布以及相应的面物体成像和不同相位差的MTF曲线。

      图  8  近似长条状湍流对系统MTF影响(相位差(a) π/2 (3π/2);(b) π面物体成像相位差;(c) π/2 (3π/2);(d) π;(e)不同相位差的MTF曲线

      Figure 8.  Influence of similar elongated turbulence on MTF of system ((a) π/2 (3π/2); (b) π plane object imaging phase difference; (c) π/2 (3π/2); (d) π; (e) MTF curves with different phase differences

      图8可知,当近似长条状湍流引起相位差逐渐改变时,MTF的中心区域没有移动,但中心区域范围减小,中心亮斑和旁瓣产生线条状的畸变,除中心处外,线条处的能量最大。图8(c)(d)为面物体成像,相比近似圆形湍流,近似长条状湍流对系统的成像质量影响更加严重。从图8(e)定量分析,以相位差$0\pi $$ (2\pi )$为基准,相位差${\pi / 2}({{3\pi } / 2})$的MTF在归一化空间频率0.17处下降0.537,0.39处(MTF为0)下降0.28,0.69处下降0.144;相位差π的MTF畸变程度非常大,MTF值在频率0.078、0.21、0.34、0.48、0.63和0.75处下降为0,同时MTF除了在中央零频有极大值,分别在频率0.14、0.27、0.42、0.54、0.68和0.81处有极大值,分别为0.21、0.40、0.14、0.33、0.028和0.044。由此可知,近似长条状湍流对系统的MTF影响非常明显。

      另外,利用Fluent软件计算光波长550 nm、风速分别为10、20、30 m/s的强风通过山脉产生的山波湍流,再通过计算得到近似长条状湍流区域的相位差分布,以此来研究近似长条状湍流对系统MTF的影响。便于分析计算出不同风速湍流对应相位差的相位均值分别为9.39 π、9.54 π、9.63 π。图9给出不同风速影响下的MTF分布和MTF曲线。

      图  9  不同风速影响的MTF ((a) 10 m/s;(b) 20 m/s;(c) 30 m/s)及(d)不同风速影响的MTF曲线

      Figure 9.  MTF affected by different wind speeds ((a) 10 m/s; (b) 20 m/s; (c) 30 m/s) and (d) MTF curves of different wind speed influence

      图9结合Fluent仿真,对流湍流的密度分布在两种气流的交界处变化梯度非常大,进而不同风速下近似长条状湍流,除了交界处相位差有明显差异,其余湍流域相位差相差不大,因此MTF只在较小空间频率处有明显的差异。由图9(d)MTF曲线定量分析,以理想情况的MTF为基准,风速为10、20、30 m/s的MTF曲线在归一化空间频率0.16处分别下降0.23、0.24、0.25,在频率0.42处均下降约0.089,在频率0.69处均下降0.046。同时三条MTF曲线除在中央零频有极大值,分别在频率0.25和0.50处出现极大值,分别为0.37和0.32。

      图3图4的近似湍流示意图分析,成像系统实际上是一个低通滤波器,根据空域与频域的对应关系可知,通光孔径不同空间位置将对应不同的空间频率范围。在近似圆形和近似长条状湍流面积占成像系统视场比相差不大情况下,近似圆形湍流只对目标图像某一频段产生影响,近似长条状湍流由于贯穿了整个孔径,因此对目标图像低、中、高频整个频段产生影响。同时由图5图8分析,附加相位差相同时,近似长条状湍流影响的MTF变化更加明显。由图7图9MTF曲线分析,除去多孔径系统附加相位差2 π周期效应,近似圆形湍流强度为${\rm{1}}{{\rm{0}}^{{\rm{ - 8}}}}\; {{\rm{m}}^{{\rm{3 - }}\alpha }}$的有效均值相位差为0.82 π,近似长条状湍流风速为20 m/s的有效相位差为0.54 π,虽然近似圆形湍流的有效相位差较大,但是近似长条状湍流对MTF的影响更显著。

    • 为了分析较小尺度湍流对多孔径光学系统的影响,对多孔径望远系统设计使用提供参考,建立了一个Golay3结构的湍流成像分析模型,分别推导了近似圆形湍流和近似长条状湍流影响下的光学系统PSF,并构建了光场相位计算模型,计算得到了两种稳态湍流影响下的成像系统MTF分布。当光路中存在湍流将导致光学系统成像质量下降。近似圆形湍流中,计算了不同高度和不同强度的湍流带来的影响,当高度大于或等于20 km时,湍流引起的相位差非常小,对系统MTF影响不大。此外,同一高度湍流强度越小,产生的相位差越小,对系统的MTF影响越小。对近似长条状湍流,计算了不同风速湍流带来的影响,当风速变化时,对MTF的影响并未发生明显改变。近似圆形湍流对成像目标某一频段产生影响,近似长条状湍流对成像目标低、中、高频整个频段产生影响。在较强湍流下,近似圆形湍流影响的MTF在归一化空间频率0.16、0.45、0.69处分别下降0.08、0.05、0.03左右。在较大风速下,近似长条状湍流影响的MTF在频率0.16、0.42、0.69处分别下降0.25、0.09、0.05左右。相比而言,后者对系统成像影响非常大,使得MTF分布产生明显条状畸变。

参考文献 (17)

目录

    /

    返回文章
    返回