留言板

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

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

基于曲率/斜率混合传感的镜面视宁度检测方法

安其昌 吴小霞 张景旭 李洪文 陈璐

安其昌, 吴小霞, 张景旭, 李洪文, 陈璐. 基于曲率/斜率混合传感的镜面视宁度检测方法[J]. 红外与激光工程, 2021, 50(7): 20200419. doi: 10.3788/IRLA20200419
引用本文: 安其昌, 吴小霞, 张景旭, 李洪文, 陈璐. 基于曲率/斜率混合传感的镜面视宁度检测方法[J]. 红外与激光工程, 2021, 50(7): 20200419. doi: 10.3788/IRLA20200419
An Qichang, Wu Xiaoxia, Zhang Jingxu, Li Hongwen, Chen Lu. Detection method of mirror seeing based on curvature/ slope hybrid sensing[J]. Infrared and Laser Engineering, 2021, 50(7): 20200419. doi: 10.3788/IRLA20200419
Citation: An Qichang, Wu Xiaoxia, Zhang Jingxu, Li Hongwen, Chen Lu. Detection method of mirror seeing based on curvature/ slope hybrid sensing[J]. Infrared and Laser Engineering, 2021, 50(7): 20200419. doi: 10.3788/IRLA20200419

基于曲率/斜率混合传感的镜面视宁度检测方法

doi: 10.3788/IRLA20200419
基金项目: 中国科学院青年创新促进会(2020221);国家自然科学基金(62005279);中国科学院装备研制项目(YJKYYQ20200057);中国科学院大学2020年度大学生创新实践训练计划(2019005770)
详细信息
    作者简介:

    安其昌,男,助理研究员,博士,主要从事大口径光机系统检测装调方面的研究

    通讯作者: 陈璐,女,助理研究员,博士,主要从事自适应光学系统波前检测方面的研究。
  • 中图分类号: TH751

Detection method of mirror seeing based on curvature/ slope hybrid sensing

  • 摘要: 随着光学元件的口径增加,无论是在加工检测还是在站址观测过程中,镜面视宁度的影响已经越来越难以忽略。在此情况下,研究可以定量测量镜面视宁度的方法与设备就显得十分重要。基于曲率/斜率混合传感策略,结合冷冻湍流假设实现镜面视宁度评价在时-空域上的转化,实现对镜面视宁度大动态范围的检测。首先,对曲率/斜率混合传感的基本方法进行了理论推导,并以标准化点源敏感性作为评价指标进一步分析了检测过程中的镜面视宁度变化。实验结果显示,在热扰动流场较为均匀时,镜面视宁度的影响较流场反复变化时小,由于湍流不稳定所引入点源敏感性(PSSn)为0.971 8。实验结果证明基于曲率/斜率混合传感的方法可以定量分析大口径光学元件镜面视宁度,对于后续开展的大口径系统设计优化与检测加工具有很好的指导意义。
  • 图  1  镜面视宁度测量过程图

    Figure  1.  Process map of mirror seeing measurement

    图  2  基于斜率的波前检测

    Figure  2.  Slope-based wavefront detection

    图  3  曲率传感实验示意图

    Figure  3.  Schematic diagram of curvature sensing experiment

    图  4  曲率传感实验现场图

    Figure  4.  Site map of curvature sensing experiment

    图  5  曲率传感实验现场采集图

    Figure  5.  Result map of curvature sensing experiment

    图  6  曲率传感实验结果图

    Figure  6.  Curvature sensing experimental results

    图  7  Zernike多项式系数变化

    Figure  7.  Variation of coefficient of Zernike polynomial

    图  8  镜面视宁度动态测量

    Figure  8.  Dynamic measurement of mirror seeing

  • [1] Andrei Tokovinin. Measurement of seeing and the atmospheric time constant by differential scintillations [J]. Appl Opt, 2002: 41(6): 957-964.
    [2] Tokovinin A. From differential image motion to seeing [J]. Publications of the Astronomical Society of the Pacific, 2002, 114: 1156-1166.
    [3] Banyal R K, Ravindra B, Chatterjee S. Opto-thermal analysis of a lightweighted mirror for solar telescope [J]. Optics Express, 2013, 21(6): 7065-7081.
    [4] Zhu Ran, Gu Bozhong, Xu jieqian, et al. Thermal control technology of primary mirror of 2.5 m class solar telescope [J]. Acta Optica Sinica, 2018, 38(11): 140-149. (in Chinese)
    [5] Yang Fei, An Qichang, Zhang Jing, et al. Seeing metrology of large aperture mirror of telescope [J]. Optics and Precision Engineering, 2017, 25(10): 2572-2579. (in Chinese) doi:  10.3788/OPE.20172510.2572
    [6] Yang Fei, An Qichang. Evaluation of specular luminance based on standardized point source sensitivity [J]. Optics and Precision Engineering, 2016, 24(5): 979-985. (in Chinese) doi:  10.3788/OPE.20162405.0979
    [7] Xia Yongbo, Xie Yongjun, Wang Peng, et al. Evaluation on mirror seeing for AIMS solar telescope [C]//Third International Conference on Photonics and Optical Engineering, SPIE, 2019, 11052: 110520O.
    [8] Zhang Enpeng, Cui Xiangqun, Li Guoping, et al. Preliminary numerical simulation of mirror seeing for the Chinese Future Giant Telescope [J]. Research in Astronomy and Astrophysics, 2016, 16(6): 123-128.
    [9] Vogiatzis K, Segurson A, Angeli G. Estimating theeffect of wind loading on extremely large telescope performance usingcomputational fluid dynamics [C]//Proceedings of SPIE - The International Society for Optical Engineering, 2004, 5497: 311-320.
    [10] Masanori I, Takeshi N, Yasuo T, et al. Evaluation of seeing on a 62-cm mirror [J]. Astronomical Society of the Pacific, 1991, 103(665): 712-722.
    [11] Bridgeland M T, Jenkins C R. Measurements of mirror seeing in the laboratory and at the telescope [J]. New Astronomy Reviews, 1998, 42(6-8): 435-440.
    [12] Beckers J M, Williams J T. Performance Of The Multiple Mirror Telescope (MMT) III. Seeing Experiments With The MMT [C]//Astronomy Conferences. International Society for Optics and Photonics, SPIE, 1982, 332: 16-23.
    [13] Tallis M, Bailey V P, Macintosh B, et al. Effects of mirror seeing on high-contrast adaptive optics instruments [J]. Journal of Astronomical Telescopes, Instruments, and Systems, 2020, 6(1): 1.
    [14] Lowne C M. An investigation of the effects of mirror temperature upon telescope seeing [J]. Monthly Notices of the Royal Astronomical Society, 1979, 188(2): 249-259. doi:  10.1093/mnras/188.2.249
    [15] Racine R, Salmon D, Cowley J D, et al. Mirror, dome, and natural seeing at CFHT[J]. Astronomical Society of the Pacific, 1991, 103(667): 1020-1032.
    [16] Roddier F. Curvature sensing and compensation: a new concept in adaptive optics [J]. Applied Optics, 1988, 27(2): 1223-1225.
    [17] An Qichang, Zhang Jingxu, Yang Fei, et al. Sensitivity analysis of standardized point source for large-scale synthetic aperture telescope [J]. Chinese Optics, 2019, 12(3): 567-574. (in Chinese) doi:  10.3788/co.20191203.0567
    [18] Yang Fei, Zhang Xuejun, Zhao Hongchao, et al. Relay optical function and pre-construction results of a Giant Steerable Science Mirror for a thirty meter telescope [J]. Optics Express, 2019, 27(10): 13991.
    [19] An Qichang, Zhang Jingxu, Yang Fei, et al. Normalized point source sensitivity analysis in GSSM prototype [J]. Chinese Optics Letters, 2017, 15(11): 50-53.
    [20] Chen Lisi, Hu Zhongwen. A study on the application of normalized point source sensitivity in wide field optical spectrometer of the thirty meter telescope [J]. Chinese Astronomy and Astrophysics, 2017, 41(4): 588-600.
    [21] Angeli G Z, Xin B, Claver C, et al. An integrated modeling framework for the Large Synoptic Survey Telescope (LSST) [C]//Proceedings of SPIE, 2016, 9911: 991118.
    [22] Angeli G Z, Bernstein R, Walls B, et al. Systems engineering for the Giant Magellan Telescope [C]//Proceedings of SPIE, 2018, 10705: 107050I.
    [23] Xin B, Claver C, Liang M, et al. Curvature wavefront sensing for the large synoptic survey telescope [J]. Applied Optics, 2015, 54(30): 9045-9054. doi:  10.1364/AO.54.009045
  • [1] 王焱, 杨永兴, 李金鹏, 毕勇, 朱庆生.  大型太阳望远镜镜面视宁度的实验研究 . 红外与激光工程, 2024, 53(1): 20230412-1-20230412-11. doi: 10.3788/IRLA20230412
    [2] 安其昌, 吴小霞, 李洪文.  基于曲率传感的稀疏孔径望远镜共焦检测方法研究 . 红外与激光工程, 2023, 52(10): 20230050-1-20230050-9. doi: 10.3788/IRLA20230050
    [3] 朱嘉康, 安其昌, 杨飞.  大口径望远镜镜面视宁度检测方法综述 . 红外与激光工程, 2023, 52(2): 20220488-1-20220488-18. doi: 10.3788/IRLA20220488
    [4] 刘欣, 杨忠明, 刘兆军, 王继红.  大口径快反镜面形测试系统设计 . 红外与激光工程, 2022, 51(5): 20210514-1-20210514-9. doi: 10.3788/IRLA20210514
    [5] 丁磊, 余烈, 谢勤岚, 朱兰艳.  光纤马赫曾德尔干涉曲率传感器的膝关节曲率检测系统研究 . 红外与激光工程, 2021, 50(12): 20210195-1-20210195-9. doi: 10.3788/IRLA20210195
    [6] 安其昌, 吴小霞, 张景旭, 李洪文, 王亮.  大口径主动光学巡天望远镜大动态范围曲率传感 . 红外与激光工程, 2021, 50(10): 20210224-1-20210224-9. doi: 10.3788/IRLA20210224
    [7] 张禹, 杨忠明, 刘兆军, 王继红.  大口径多光谱通道波前测量系统设计 . 红外与激光工程, 2020, 49(8): 20190559-1-20190559-7. doi: 10.3788/IRLA20190559
    [8] 唐境, 张景旭, 安其昌, 李洪文.  大口径巡天望远镜校正镜弹性体支撑 . 红外与激光工程, 2020, 49(S1): 20200124-20200124. doi: 10.3788/IRLA20200124
    [9] 黄宇飞, 白绍竣, 高冀, 吕争, 徐嘉.  大口径空间光学反射镜面形动力学响应分析 . 红外与激光工程, 2019, 48(11): 1114001-1114001(6). doi: 10.3788/IRLA201948.1114001
    [10] 赵宪宇, 薛栋林, 程强.  哈特曼原理子口径斜率扫描检测及误差研究 . 红外与激光工程, 2019, 48(8): 813003-0813003(8). doi: 10.3788/IRLA201948.0813003
    [11] 杨慧珍, 王斌, 刘瑞明, 马良.  模型式无波前探测自适应光学系统抗噪能力分析 . 红外与激光工程, 2017, 46(8): 817002-0817002(6). doi: 10.3788/IRLA201746.0817002
    [12] 朱进一, 谢永军.  采用衍射主镜的大口径激光雷达接收光学系统 . 红外与激光工程, 2017, 46(5): 518001-0518001(8). doi: 10.3788/IRLA201746.0518001
    [13] 王晓鹏, 刘克, 李艳秋.  基于子孔径波前振幅调制的斜率和曲率信号提取算法 . 红外与激光工程, 2017, 46(12): 1217003-1217003(7). doi: 10.3788/IRLA201746.1217003
    [14] 杨慧珍, 刘荣, 刘强.  基于变形镜本征模的模型式无波前探测自适应光学系统 . 红外与激光工程, 2015, 44(12): 3639-3644.
    [15] 李宏壮, 赵勇志, 马鑫雪, 于树海, 殷丽梅.  大口径折反式三组元连续变焦距光学系统设计 . 红外与激光工程, 2015, 44(10): 3037-3042.
    [16] 安其昌, 张景旭, 杨飞, 乔兵.  基于斜率的TMT 三镜面形检测方法 . 红外与激光工程, 2015, 44(6): 1884-1889.
    [17] 马鑫雪, 王建立, 王斌.  利用相位恢复波前传感技术检测球面镜面形 . 红外与激光工程, 2014, 43(10): 3428-3433.
    [18] 杨飞, 刘国军, 安其昌.  结构函数在大口径光学系统评价中的应用 . 红外与激光工程, 2014, 43(11): 3832-3836.
    [19] 程强, 薛栋林, 张学军.  相位差异的三反光学系统波前传感技术 . 红外与激光工程, 2013, 42(6): 1601-1606.
    [20] 明名, 杨飞, 赵金宇, 张丽敏, 吴小霞.  折反式大口径、大视场、宽光谱光学系统 . 红外与激光工程, 2012, 41(1): 149-154.
  • 加载中
图(8)
计量
  • 文章访问数:  241
  • HTML全文浏览量:  95
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-17
  • 修回日期:  2021-01-07
  • 刊出日期:  2021-07-25

基于曲率/斜率混合传感的镜面视宁度检测方法

doi: 10.3788/IRLA20200419
    作者简介:

    安其昌,男,助理研究员,博士,主要从事大口径光机系统检测装调方面的研究

    通讯作者: 陈璐,女,助理研究员,博士,主要从事自适应光学系统波前检测方面的研究。
基金项目:  中国科学院青年创新促进会(2020221);国家自然科学基金(62005279);中国科学院装备研制项目(YJKYYQ20200057);中国科学院大学2020年度大学生创新实践训练计划(2019005770)
  • 中图分类号: TH751

摘要: 随着光学元件的口径增加,无论是在加工检测还是在站址观测过程中,镜面视宁度的影响已经越来越难以忽略。在此情况下,研究可以定量测量镜面视宁度的方法与设备就显得十分重要。基于曲率/斜率混合传感策略,结合冷冻湍流假设实现镜面视宁度评价在时-空域上的转化,实现对镜面视宁度大动态范围的检测。首先,对曲率/斜率混合传感的基本方法进行了理论推导,并以标准化点源敏感性作为评价指标进一步分析了检测过程中的镜面视宁度变化。实验结果显示,在热扰动流场较为均匀时,镜面视宁度的影响较流场反复变化时小,由于湍流不稳定所引入点源敏感性(PSSn)为0.971 8。实验结果证明基于曲率/斜率混合传感的方法可以定量分析大口径光学元件镜面视宁度,对于后续开展的大口径系统设计优化与检测加工具有很好的指导意义。

English Abstract

    • 为了对更加深远的宇宙进行探索,进一步围绕宇宙起源、暗物质暗能量等科学目标开展相关研究,需要加大望远镜的口径,以获得更高的分辨率与灵敏度。传统的望远镜口径较小,其光学元件尺度小于大气相干长度,因此很少考虑由于大气湍流影响所引起像质的损失,尤其是在加工检测过程中的气流影响往往容易被忽略。但是,随着光学系统口径的增加,镜面附近湍流影响已经越来越不可忽略。单纯依靠多次测量平均所获得的测量结果,不仅在效率上难以保证,同时,由于测量环境难以满足严格的各向同性,平均的结果可能存在一定的统计学偏置,不仅会影响光学表面散射特性,严重时甚至会造成光学加工不收敛。不仅如此,对于下一代地基大口径望远镜,其口径将会超过30 m,以美国的30 m望远镜为例,为了充分利用望远镜的集光能力,在光学平台上配备了诸多科学终端,系统所收集的光线在中继至对应终端的过程中,主镜镜面附近的大气会对光线发生偏折,进而影响光线中继的质量。

      目前,国内外诸多学者已针对高层大气与圆顶内部湍流所引入光程差展开了诸多相关研究[1-3],而由于镜面局部热不均匀性所引入“镜面视宁度”却缺少相关的分析手段,不仅如此,与镜面视宁度类似,系统终端附近也会产生类似的相位扰动。

      与高层大气所引入的波前误差相比,镜面视宁度受边界条件影响较大,无法很好地通过统计学模型进行描述[4~6]。镜面视宁度检测主要可以分为直接检测、间接估计以及数值仿真三个技术流派。具体来说:在计算仿真方面,Xia等[7]针对太阳望远镜观测阶段的镜面湍流情况,采用流仿真软件对镜面视宁度进行了分析,从数值计算的角度验证了强迫对流对镜面视宁度的显著提升作用。基于同样流体力学的原理,Zhang等[8]也利用仿真软件对镜面视宁度对成像探测的影响进行了分析。但是,使用流体仿真分析法[9]对于模型的精确程度与计算机算力均有着较高的要求。因此,为了获得与实际边界更接近的结果,多使用在线实际检测。Iye.M等针对620 mm的反射镜面,采用哈特曼波前传感器测量了镜面附近的湍流情况与周围空气温差的关系[10~12]。Melisa Tallis等[13]利用GIMINI south的科学终端实现了镜面视宁度的测量,分析了镜面视宁度对高度比度的系外行星观测的影响。除使用成像终端进行测量外,镜面视宁度也可以通过温度数据进行间接估计[11],即通过光学元件背部附加的温度传感器建立温差数据与镜面视宁度之间的函数关系,进而对镜面视宁度进行估计。

      对于镜面视宁度而言,其测量的难点与重点均在于与实际检测边界以及温度场之间的交互与耦合,因此,针对镜面视宁度,采用直接测量是最为可行的方案。目前,镜面视宁度直接检测多采用干涉仪[14]或焦面探测部件[15],但是以上两种技术路线均需要占用光学元件光轴方向的空间。在需要在线检测的场景,反复法如拆装迭代将会严重影响元件的加工周期。

      针对以上情况,文中提出了一种基于曲率/斜率混合传感的镜面视宁度检测方法,基于“冷冻湍流”假设(大气湍流为固定分布的“相位屏”按照一定速度匀速前进),利用湍流的定向运动,通过孔径分时复用换取相当的空间分辨能力(可类比非快照系统)。基于测量所得到的斜率与曲率信号,不仅可以对镜面视宁度进行评价,而且可以结合横向风速对测量光路上的大气湍流情况进行重建。

    • 在实际应用中,使用泽尼克(Zernike)多项式等圆域连续基底,不仅会出现非圆域正交性退化的问题,同时在离散取样时,其表征精度会出现相应的下降。在此选取傅里叶级数对波前进行表达,同时,由于傅里叶变换可以使用快速算法,其计算速度也可以大幅提高。

      设最高的拟合阶数为N,将波前Φ(x, y)表示为离散傅里叶级数,如公式(1)所示:

      $$\varPhi (x,\;y) = \sum\limits_{n = 1}^N {\sum\limits_{m = 1}^N {\tilde \varPhi \left( {n,m} \right)} } {{\rm{e}}^{j\frac{{nx + my}}{N}}}$$ (1)

      式中:n、m为整数;xy为离散空间坐标;$\tilde \varPhi $(n, m)为离散傅里叶级数系数。将公式(1)变形可得公式(2):

      $$\tilde \varPhi (m,\;n) = \sum\limits_{x = 1}^N {\sum\limits_{y = 1}^N {\varPhi \left( {x,\;y} \right)} } {{\rm{e}}^{ - j\frac{{nx + my}}{N}}}$$ (2)

      将公式(1)两端进行梯度算子:

      $${\nabla ^2}\varPhi (x,\;y) = \sum\limits_{n = 1}^N {\sum\limits_{m = 1}^N {\frac{{{n^2} + {m^2}}}{{{N^2}}}\tilde \varPhi \left( {n,\;m} \right)} } {{\rm{e}}^{j\frac{{nx + my}}{N}}}$$ (3)

      曲率传感器由Roddier于1988年提出[16],其基本原理是光瞳处波前局部的曲率变化,所对应的焦内像与焦外像的光强分布会发生对应的变化。根据近场电磁波的传输方程可以解算出波前信息,曲率传感因具有非干涉(无需参考光)、结构简单、环境适应性好、解算稳定(无需迭代)、孔径遮拦影响小以及动态范围大(无需相位解缠)等诸多独特优势,具有非常好的应用价值。

      根据曲率传感基本原理,波前相位的曲率与光强沿光轴方向的差分成正比,如公式(4)所示:

      $$\frac{{{I_1}(x,\;y) - {I_2}(x,\;y)}}{{{I_0} \cdot 2\Delta {\textit{z}}}} \approx - {\nabla ^2}\varPhi (x,\;y)$$ (4)

      式中:I1(x, y)、I2(x, y)为焦前焦后的能量分布;I0为光强总值;Δz为离焦量。结合公式(3)~(4)可得:

      $$\begin{array}{l} \dfrac{{{I_1}(x,\;y) - {I_2}(x,\;y)}}{{{I_0} \cdot 2\Delta {\textit{z}}}} = - \displaystyle\sum\limits_{n = 1}^N {\sum\limits_{m = 1}^N {\dfrac{{{n^2} + {m^2}}}{{{N^2}}}\tilde \varPhi \left( {n,m} \right)} } {{\rm{e}}^{j\frac{{nx + my}}{N}}} \end{array} $$ (5)

      对公式(5)两端进行变形可得公式(6):

      $$ \begin{split} & \displaystyle\sum\limits_{x = 1}^N {\sum\limits_{y = 1}^N {\dfrac{{{I_1}(x,\;y) - {I_2}(x,\;y)}}{{{I_0} \cdot 2\Delta {\textit{z}}}}{{\rm{e}}^{ - j\frac{{nx + my}}{N}}}} } = \\ & -\displaystyle\sum\limits_{x = 1}^N {\sum\limits_{y = 1}^N {\sum\limits_{n = 1}^N {\sum\limits_{m = 1}^N {\frac{{{n^2} + {m^2}}}{{{N^2}}}\tilde \varPhi \left( {n,\;m} \right)} } {\rm{e}}^{j\frac{{nx + my}}{N}}} {{\rm{e}}^{ - j\frac{{nx + my}}{N}}}} \end{split} $$ (6)

      进而可得公式(7),通过公式(1)可得波前Φ(x, y):

      $$\tilde \varPhi \left( {n,\;m} \right){\rm{ = }}\sum\limits_{x = 1}^N {\sum\limits_{y = 1}^N {\frac{{{I_2}(x,\;y) - {I_1}(x,\;y)}}{{\left( {{n^2} + {m^2}} \right){I_0} \cdot 2\Delta {\textit{z}}}}{{\rm{e}}^{ - j\frac{{nx + my}}{N}}}} } $$ (7)

      类似地,可以建立$\tilde \varPhi $(n, m)与斜率S(x, y)的关系,如公式(8)所示:

      $$\tilde \varPhi \left( {n,\;m} \right){\rm{ = }}\sum\limits_{x = 1}^N {\sum\limits_{y = 1}^N {\frac{{S(x,\;y)}}{{n + m}}{{\rm{e}}^{ - j\frac{{nx + my}}{N}}}} } $$ (8)
    • 大口径光学元件系统性能评价和误差分配是一个复杂的系统工程,其涉及面广、环节众多、与实际应用需求联系极为紧密。误差分析、分配准则的研究本质上是建立大口径望远镜科学目标与望远镜实际性能指标之间的关系。不仅如此,性能评价和误差分配实际上是相互关联的,即误差分配本身是一个平衡过程:在既定的性能评价标准下,需要综合考虑系统误差的敏感性和实际加工、生产能力以及成本等因素,进而实现误差的分配。标准化点源敏感性(Point Source Sensitivity,PSSn)是点扩散函数在成像区域的平均[17-20]。PSSn的定义与计算方法如公式(9)所示:

      $$PSSn = \dfrac{{\displaystyle\iint\limits_{A(x,\;y)} {{{\left| {PS{F_e}} \right|}^2}}{{\left| {PS{F_{t + a}}} \right|}^2}}}{{\displaystyle\iint\limits_{A(x,\;y)} {{{\left| {PS{F_{t + a}}} \right|}^2}}}} = \dfrac{{\displaystyle\iint\limits_{A(fx,\;fy)} {{{\left| {OT{F_e}} \right|}^2}}{{\left| {OT{F_{t + a}}} \right|}^2}}}{{\displaystyle\iint\limits_{A\left( {x,\;y} \right)} {{{\left| {OT{F_{t + a}}} \right|}^2}}}} $$ (9)

      式中:PSFePSFt+a为望远镜误差的点扩散函数以及理想望远镜在视宁影响下的点扩散函数;OTFeOTFt+a为望远镜误差的光学传递函数以及理想望远镜在视宁影响下的光学传递函数;A(x, y)与A(fx, fy)为空间域与频域的积分区域。根据公式(1),可选择合适的方式计算PSSn。

      Seo在2009年提出PSSn之后,TMT全频段的误差分配均使用PSSn进行表征。Angeli在2016年对LSST进行了系统级建模,利用PSSn分析了各个频段的误差[21]。Angeli在2018年利用PSSn对GMT自然视宁与地表自适应光学两种工作模式行了误差分配,引入进入大于1的PSSn,用于模拟自适应光学中频误差校正环节[22]

    • 对流方向模拟实际使用过程中的不同方向采样,并通过对自准直仪示数的直接分析定性地获得湍流强度,指导系统进行初步的温控与湍流调整,同时利用公式(5)也可以对镜面视宁进行初步计算,如图1所示。

      图  1  镜面视宁度测量过程图

      Figure 1.  Process map of mirror seeing measurement

      为了兼顾测量动态范围以及测量精度,在此采用量级测量的策略,即采用动态范围较大的斜率测量对强湍流情况进行估计,之后采用曲率传感测策略对镜面视宁度进行更加精细与准确的测量。

      在进行斜率探测后,通过自准直仪后端光路的调焦将系统进行一定量的离焦(两台自准直仪的离焦方向需要相反),并通过波前沿光轴方向的微分获得波前的估计。获得数据后,基于PSSn这一频域指标对结果进行表征与分析,以期更全面地理解镜面视宁度作用的内在物理机制与统计学特性,最终平衡镜面视宁度、光学误差以及加工成本之间的关系。由于斜率测量采用大像元尺寸(提高读出速度与读出噪声),当其测量值接近其测角精度时,说明热控系统已经抑制部分湍流,在此情况下可切换曲率传感。

      具体来说,采用系统测量分为两级,分别是大范围的斜率测量以及小量程高精度的曲率测量。由于斜率探测仅需要几个像元,动态范围较大,因此作为初始阶段的镜面视宁度检测手段。

      在镜面视宁度精确测量中,曲率传感因具有非干涉(无需参考光)、结构简单、环境适应性好、解算稳定(无需迭代)、孔径遮拦影响小以及动态范围大(无需相位解缠)等诸多独特优势,具有非常好的应用价值。但高精度的离焦器件系统限制了其应用的场景,在此采用分体球面镜实现固定的离焦量,并最终实现单次曝光测量的小型化、静态化的镜面视宁度测量仪。

    • 在实际的检测过程中,使用极坐标的便利程度远远超过了直角坐标,在此研究了基于极坐标的镜面视宁度检测方法,通过额定160 W的设备,配合风扇进行镜面视宁模拟,其检测结果如图2所示。

      图  2  基于斜率的波前检测

      Figure 2.  Slope-based wavefront detection

      为了对镜面视宁度进行更为精确的解算,在斜率测量的基础上需要利用曲率传感进行测量。曲率传感器所对应的焦内像与焦外像的光强分布会发生对应的变化[23]。在此,采用通过固定部件(错位球面镜)分视场单次曝光获取离焦像的物理架构,模拟实际过程中的曲率传感过程,并根据近场电磁波的传输方程解算出波前信息。采用的球面镜F数为10,在此忽略由于光束数值孔径所引入的误差,曲率传感实验示意图如图3所示,曲率传感实验现场如图4~5所示,曲率传感实验结果如图6所示。实验的目的在于分析望远镜镜面视宁的测量方法以及气流运动对成像效果与波前探测的关系,面向实际的系统,可以将该实验装置与实际光路进行耦合(可采用共路或垂直光路安放),对瞬时的局部视宁度进行测量,在与相机的采样保持同步的基础上,可以挑选视宁度影响较小的短曝光图像,并在后续的计算过程中针对筛选过的图像进行平均等后续处理,可从探测机理上提升主动光学系统精度溯源的能力。在此情况下,错位型曲率传感器单次曝光的优势更加凸显。

      图  3  曲率传感实验示意图

      Figure 3.  Schematic diagram of curvature sensing experiment

      图  4  曲率传感实验现场图

      Figure 4.  Site map of curvature sensing experiment

      图  5  曲率传感实验现场采集图

      Figure 5.  Result map of curvature sensing experiment

      图  6  曲率传感实验结果图

      Figure 6.  Curvature sensing experimental results

      波前传感所对应的各个Zernike多项式系数如图7所示。将视宁度的影响限制在低阶部分,一方面考虑由于镜面局部热不平衡带来的影响在空间频率上较由于多层大气所引入的大气视宁度低,另一方面,考虑到主动校正系统主要针对低阶像差的因素。可见,由于强迫对流的影响,各个空间频率的像差均发生了改变,但是较低阶的影响要更加明显。冷冻湍流假设认为在采样率足够的情况下,大气湍流为固定分布的“相位屏”按照一定速度匀速前进,在固定口径中,通过足够长时间的积分可获得接近高斯分布的相位分布,基于该假设在大多数情况下可实现大气扰动的估计,如果进行更高精度的波前传感,一般来说需要更长的积分时间。但是积分时间是涉及系统带宽、控制残差以及观测效率等多个系统关键参数的物理量,难以无限增加。该设备可在不依靠大量数据的前提下体现湍流本身的统计学特性,通过假设检验对所采数据的波前进行分析,从探测机理上提高面对非平稳过程的鲁棒性。

      图  7  Zernike多项式系数变化

      Figure 7.  Variation of coefficient of Zernike polynomial

      在静态测量的基础上进行动态测量,镜面视宁度动态测量如图8所示。考虑对镜面视宁度进行动态采集,并在50 s前对湍流源进行了方向上的微动,使气流更加均匀,在50 s后仅进行静态的强迫对流,系统的镜面视宁度进一步下降。之前的研究充分肯定了强迫对流方法对镜面视宁度的提升效果,在此,针对强迫对流的稳定性与镜面视宁度之间的关系进行了研究。可见,均匀且各项同性的流场对于波前检测过程的影响较小,由于湍流不稳定所引入PSSn为0.9718,对望远镜观测跟踪中的镜面视宁度控制具有非常好的指导作用。

      图  8  镜面视宁度动态测量

      Figure 8.  Dynamic measurement of mirror seeing

    • 近年来,大口径望远镜的发展几乎主要集中于克服大气所带来的影响。镜面视宁度是大口径光学元件在制造、检测、装调过程中都十分重要的影响因素。

      因此,文中提出基于曲率/斜率混合传感的镜面视宁度检测方法,结合标准化点源敏感性,不仅可以在加工检测的过程中对镜面视宁度静态分布与动态变化进行有效的测量、评价,而且可以在横向风边界条件下对路线上的大气湍流的空间分布情况进行反演,实现非充分曝光下的大气湍流影响解耦,即在大气扰动没有被完全平均的情况下实现高精度波前传感。

      不仅如此,针对主动光学系统波前传感系统局部的大气扰动,也可以在不附加其他硬件的情况下,通过该方法进行气流影响解耦。与传统的利用热控系统进行控制的方法相比,成本更低,适用更强,对于物理空间难以排布主动热控系统的米级口径望远镜具有十分良好的应用场景。

参考文献 (23)

目录

    /

    返回文章
    返回