搜索

x

留言板

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

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

大气湍流对接收光场时间相干特性的影响

艾则孜姑丽·阿不都克热木 陶志炜 刘世韦 李艳玲 饶瑞中 任益充

引用本文:
Citation:

大气湍流对接收光场时间相干特性的影响

艾则孜姑丽·阿不都克热木, 陶志炜, 刘世韦, 李艳玲, 饶瑞中, 任益充

Influence of atmospheric turbulence on temporal coherence characteristics of received optical field

Abdikirim Azizigul, Tao Zhi-Wei, Liu Shi-Wei, Li Yan-Ling, Rao Rui-Zhong, Ren Yi-Chong
PDF
HTML
导出引用
计量
  • 文章访问数:  1100
  • PDF下载量:  60
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-06-17
  • 修回日期:  2022-09-01
  • 上网日期:  2022-11-17
  • 刊出日期:  2022-12-05

大气湍流对接收光场时间相干特性的影响

  • 1. 中国科学院合肥物质科学研究院, 安徽光学精密机械研究所, 中国科学院大气光学重点实验室, 合肥 230031
  • 2. 中国科学技术大学研究生院科学岛分院, 合肥 230026
  • 3. 中国科学技术大学环境科学与光电技术学院, 合肥 230026
  • 4. 先进激光技术安徽省实验室, 合肥 230037
  • 通信作者: 任益充, rych@aiofm.ac.cn
    基金项目: 安徽省自然科学基金青年项目(批准号: 1908085QA37)、国家自然科学基金青年科学基金(批准号: 11904369)和脉冲功率激光技术国家重点实验室主任基金(批准号: 2019ZR07)资助的课题.

摘要: 激光在大气传输过程中其相干特性受湍流影响而下降, 进一步影响相干探测过程的效率和性能. 本文定义大气相干时间描述激光经大气传输后接收光场的起伏速度, 大气相干时间与相干过程持续时间的相对大小决定了相干探测过程的性能. 本文基于多层动态相位屏技术仿真激光大气传输过程, 并系统研究大气参数、收发参数、波长等对大气相干时间的影响规律. 研究指出大气相干时间与波长、接收口径、大气相干长度呈正相关, 与风速呈反比; 通过改进光学设计、光束整形、激光选频等方法能够有效减少湍流扰动的影响, 提高接收光场的稳定性, 进而改善湍流对相干探测过程的不良影响. 大气相干时间是衡量湍流对相干探测影响的重要参数, 本研究可为评估大气对相干探测过程的影响提供有力参考.

English Abstract

    • 激光的相干性在相干探测、相干成像、量子探测等领域中扮演着重要角色[1-7]. 激光在大气传输过程中, 不可避免地受到环境的影响发生退相干使其相干性下降, 致使相干探测、相干成像、量子探测等性能大幅下降. 在真实大气环境中, 湍流是致使激光相干性退化的重要因素, 湍流引起的光场退相干主要表现在时间和空间两方面. 空间上, 湍流使激光产生光斑漂移、光束扩展, 并引起激光波阵面的畸变, 激光从单模光场退化为多模光场, 致使其相干效率下降. 时间上, 湍流大气以kHz的频率随机起伏, 光强和相位随时间呈随机涨落. 关于湍流对激光相干特性影响的研究主要集中在空间部分, 如湍流对部分相干环状光束、部分相干厄密-高斯光束和斜程传输准单色高斯-谢尔光束等空间相干性[8-12]的影响, 以及接收光场与本振光场的相干效率等[13]. 时间相干特性对相干探测以及逆合成孔径激光雷达(ISAL)等相干过程至关重要, 然而目前湍流对激光光场时间相干特性影响的研究相对匮乏.

      与强度探测相比, 相干探测能够实现弱信号检测、提高探测精度/距离等, 但其对激光的相位稳定性往往有更高要求. 相干探测是持续过程, 在探测过程中湍流引起的激光光场随机起伏非常关键, 如探测时间较短, 在探测过程中湍流符合“冻结”假设, 接收光场保持稳定, 则探测结果受湍流影响较小; 反之若探测时间较长, 在探测过程中湍流的随机起伏引起接收光场的起伏, 致使接收端的相位随之涨落, 则探测结果受湍流影响较大. 以ISAL为例, ISAL向目标发射多个脉冲(啁啾)信号, 并对接收信号进行外差探测以获取目标信息. 在此过程中, 接收光场的空间相干性退化使外差效率低、信噪比下降, 此时可通过加大激光发射功率、增大接收口径、自适应光学校正等方法解决; 接收光场的时间相干性(非激光器的相干时间)退化使脉冲持续时间内接收光场的相位出现随机起伏, 使得ISAL在整个成像过程中无法保持相位稳定性, 致使ISAL成像质量差甚至无法相干合成. 在实际工程中, 湍流引起的接收光场相位的随机起伏是ISAL的技术难点, 这也正是逆合成孔径微波雷达(ISAR)早已成熟应用而ISAL还在研究中的重要原因. 因此研究大气湍流对接收光场时间相干特性的影响规律, 对于评估自由空间中相干探测精度、ISAL成像质量等有着重要参考价值.

      综上所述, 湍流影响激光接收光场时间相干性的研究在评估相干探测性能、ISAL成像质量等领域有着重要价值, 该研究能够为相干探测方案、ISAL成像方案中关键参数选取、性能优化提供有力参考. 本文基于多层动态相位屏仿真技术, 针对高斯光束的大气传输过程开展仿真, 定义大气相干时间用于描述接收光场的时间相干性, 根据仿真得到的接收光场分析大气参数、发射/接收参数以及激光参数等对大气相干时间的影响.

    • 激光在湍流大气中的传输相当于在真空中放置一系列薄的随机相位屏, 光场衍射至相位屏并被相位调制, 此衍射和相位调制过程连续进行直至传输到接收端. 这种用多层相位屏模型来代替连续随机介质的方法是激光传输仿真中最为常见的随机相位屏方法[14]. 其生成方法主要有傅里叶方法[15,16]、分形方法[17]、泽尼克多项式法[18], 添加次谐波的傅里叶方法是目前最常见的相位屏生成方法. 1999年,Assemat等[19]提出基于协方差的多层动态相位屏方法具有高精度、长曝光, 可以仿真湍流参数随时间变化等诸多特点, 本文选择该方法生成相位屏.

      图1所示, 多层动态相位屏方法基于初始随机相位屏的最后几列Z (一般取最后两列即可, 以符号⊙表示) 结合湍流功率谱密度函数生成新列X (以符号 表示), 将新列X添加进相位屏并删除左侧旧列, 即可得到新的相位屏; 在行方向进行同样操作即可生成新行, 在仿真过程中每次生成的新行和新列数量需要根据风速、风向、帧间隔、格点间距等参数确定. 新生成的列X和旧列Z满足关系式[17,20]:

      图  1  多层动态相位屏原理示意图

      Figure 1.  Schematic of infinitely long phase screen principle.

      $ \boldsymbol{X}=\boldsymbol{A}\boldsymbol{Z}+\boldsymbol{B}\boldsymbol{\beta } \text{, } $

      其中, ${\boldsymbol{\beta}} $是均值为0、方差为1的高斯随机向量[21-25]. 矩阵A, B可以从XZ矢量的协方差计算得到(关于A, B的推导以及更多的维度信息, 详见参考文献[20]). 通过多层动态相位屏方法仿真出激光在大气湍流中传输后的接收光场, 为了定量描述光场相干性变化, 下面以大气湍流对外差探测的影响为例定义其相干性.

    • 光场的相干性是较为宽泛的概念, 在不同方向或领域往往有着不同的物理意义. 如量子力学中相干性通常与密度算符非对角项有关, 激光器的相干性则由相干时间或光束质量描述等. 本文主要针对大气湍流对ISAL成像、相干探测等相干过程的影响, 这类相干过程的典型场景是激光雷达向目标发射脉冲信号并接收从目标返回的信号与本地信号进行拍频的外差信号. 对于受湍流影响的信号, 需从相干探测的角度分析大气湍流对接收光场时间相干性的影响.

      频率分别为$ {\omega _{\text{L}}} $$ {\omega _{\text{S}}} $的本振电场$ {E_{\text{L}}} $和信号电场$ {E_{\text{S}}} $的表达式为

      $ {E_{\text{L}}} = {E_{\text{L}}}\left( {x,y} \right){{\text{e}}^{ - {\text{i}}{\omega _{\text{L}}}t}} \text{, } \tag{2a}$

      $ {E_{\text{S}}} = {E_{\text{S}}}\left( {x,y} \right){{\text{e}}^{ - {\text{i}}{\omega _{\text{S}}}t}} . \tag{2b} $

      本振电场$ {E}_{\rm{L}} $和信号电场$ {E}_{\rm{S}} $中不仅包含振幅信息, 还包含$ t=0 $时的相位信息. 电场$ {E}_{\rm{L}} $$ {E}_{\rm{S}} $相干叠加后的光强为

      $ \begin{split}{I_{{\text{tot}}}}\left( {x,y,t} \right) =\;& {\left| {{E_{\text{L}}}\left( {x,y} \right)} \right|^2} + {\left| {{E_{\text{S}}}\left( {x,y} \right)} \right|^2}\\ &+ \big[ {E_{\text{L}}}\left( {x,y} \right)E_{\text{S}}^*\left( {x,y} \right){{\text{e}}^{ - {\text{i}}\left( {{\omega _{\text{L}}} - {\omega _{\text{S}}}} \right)t}} \\ &+ E_{\text{L}}^*\left( {x,y} \right){E_{\text{S}}}\left( {x,y} \right){{\text{e}}^{ - {\text{i}}\left( {{\omega _{\text{L}}} - {\omega _{\text{S}}}} \right)t}} \big] . \end{split} $

      由于探测器测得的信号为一定空间ΔS与一定时间Δt内场强叠加, 因此得到

      $ I_{{\text{tot}}}^{\text{d}} = \displaystyle \int_0^{\Delta t} {\displaystyle \int_{\Delta S} {{I_{{\text{tot}}}}\left( {x,y,t} \right){\text{d}}x{\text{d}}y{\text{d}}t} } . $

      将(3)式代入(4)式后可得

      $ \begin{split}{I}_{\text{tot}}^{\text{d}}=\;&{\displaystyle {\int }_{0}^{\Delta t}{\displaystyle {\int }_{\Delta S}{\left|{E}_{\text{L}}\left(x,y\right)\right|}^{2}+{\left|{E}_{\text{S}}\left(x,y\right)\right|}^{2}\text{d}x\text{d}y\text{d}t}}\\ &+\displaystyle\int_{0}^{\Delta t}\displaystyle\int _{\Delta S}\big[{E}_{\text{L}}\left(x,y\right){E}_{\text{S}}^{*}\left(x,y\right){\text{e}}^{-\text{i}\left({\omega }_{\text{L}}-{\omega }_{\text{S}}\right)t}\\ &+{E}_{\text{L}}^{*}\left(x,y\right){E}_{\text{S}}\left(x,y\right){\text{e}}^{-\text{i}\left({\omega }_{\text{L}}-{\omega }_{\text{S}}\right)t}\big]\text{d}x\text{d}y\text{d}t\text{ }\\ &={I}_{\text{L}}^{\text{d}}+{I}_{\text{S}}^{\text{d}}+{I}_{\text{c}}^{\text{d}}\text{, }\\[-10pt]\end{split} $

      其中, $ {I}_{\rm{L}} $$ {I}_{\rm{S}} $分别为电场$ {E}_{\rm{L}}\left(x, y\right) $$ {E}_{\rm{S}}\left(x, y\right) $的直流分量, $ {I}_{\rm{c}} $$ {E}_{\rm{L}} $, $ {E}_{\rm{S}} $的干涉项或相干项; $\displaystyle\int_0^{\Delta t} {{{\text{e}}^{ - {\text{i}}\left( {{\omega _{\text{L}}} - {\omega _{\text{S}}}} \right)t}}{\text{d}}t} $$\displaystyle\int_{\Delta S} {{E_{\text{L}}}\left( {x, y} \right)} E_{\text{S}}^ * \left( {x, y} \right){\text{d}}x{\text{d}}y$分别为时间部分和空间部分的积分结果, 代表信号光场与本振光场在时间和空间上的相干特性.

      将相干对比度衡量电场$ {E}_{\rm{L}} $$ {E}_{\rm{S}} $的相干性定义为R, 即

      $ R = \dfrac{{\left| {I_{\text{c}}^{\text{d}}} \right|}}{{I_{\text{L}}^{\text{d}} + I_{\text{S}}^{\text{d}}}} = \dfrac{\displaystyle{2\left| {\int_{\Delta S} {{E_{\text{L}}}\left( {x,y} \right)E_{\text{S}}^*\left( {x,y} \right){\text{d}}x{\text{d}}y} } \right|{{\rm{Re}}} \int_0^{\Delta t} {{{\text{e}}^{ - {\text{i}}({\omega _{\text{L}}} - {\omega _{\text{S}}})t + {\phi _{\text{S}}}\left( t \right)}}{\text{d}}t} }}{\displaystyle{\Delta t\int_{\Delta S} {{{\left| {{E_{\text{L}}}\left( {x,y} \right)} \right|}^2} + {{\left| {{E_{\text{S}}}\left( {x,y} \right)} \right|}^2}} {\text{d}}x{\text{d}}y}} . $

      由(6)式可知, 分子部分由$ {E}_{\rm{L}}\left(x, y\right) $$ {E}_{\rm{S}}\left(x, y\right) $在时间和空间上的积分决定. 当$ {E}_{\rm{L}}\left(x, y\right) $$ {E}_{\rm{S}}\left(x, y\right) $横模不匹配时, $ {E}_{\rm{L}}\left(x, y\right){E}_{\rm{S}}^{\rm{*}}\left(x, y\right) $在空间上积分值则趋向于0, 此时意味着二者空间上不相干. 湍流在信号光上附加的相位${\phi }_{\rm{S}}\left(t\right)$随时间涨落, 意味着积分时间$ \Delta t $较长时$ R $趋于0, 干涉项被时间涨落平均从而表现出不相干. 对于信号探测需要持续一定时间的成像过程而言, 电场$ {E}_{\rm{S}}\left(x, y\right) $的时间涨落会导致从目标探测回来的接收信号相位产生波动, 从而无法在脉冲持续时间内保持接收信号的相位稳定, 为准确衡量相干过程的时间相干性, 我们定义大气相干时间$ {\tau }_{\rm{c}} $用于描述$ {E}_{\rm{S}}\left(x, y\right) $的变化情况.

    • 接收光场受大气湍流的影响呈现随机起伏, 其起伏速度对相干过程的影响起决定性作用. 为衡量受大气湍流影响的接收光场的稳定性, 基于电场的自相关函数[26,27]定义大气相干时间$ {\tau }_{\rm{c}} $. 接收光场为$ {E}_{\rm{r}}(x, y, t) $, 其自相关函数$ M\left(t\right) $

      $ M(t) = \left\langle\displaystyle\iint {E}_{\rm{r}}(x,y,{t}_{0}){E}_{\rm{r}}^{\rm{*}}(x,y,{t}_{0} + t) {{\rm{d}}}x{{\rm{d}}}y\right\rangle_{{t}_{0}} , $

      其中, $ t $为相邻接收光场的时间间隔, $ {t}_{0} $为帧间隔, ${\left\langle{\cdot }\right\rangle}_{{t}_{0}}$表示在间隔$ t $时刻对$ {t}_{0} $, $ {t}_{0}+t $不同的两帧接收光场进行求和取系综平均. 根据随机场的统计理论, $ M\left(t\right) $应为单调递减函数, 采用$ a{{\text{e}}^{ - \tfrac{t}{{{\tau _{\text{c}}}}}}} + c $对自相关函数$ M\left(t\right) $进行拟合, 即$ M\left( t \right) = a{{\text{e}}^{ - \tfrac{t}{{{\tau _{\text{c}}}}}}} + c $(其中拟合参数ac为常数), 拟合所得$ {\tau }_{\rm{c}} $与接收光场的时间相干特性有关, 接收光场受湍流扰动产生随机起伏, 其越剧烈则接收光场变化越快时间相干性越弱且$ {\tau }_{\rm{c}} $越小, $ {\tau }_{\rm{c}} $可用来评价接收光场的稳定性或时间相干特性.

    • 分析大气相干时间的影响因素之前, 首先需要指出, 本文所定义的大气相干时间$ {\tau }_{\rm{c}} $与Greenwood相干时间$ {\tau }_{0} $存在相关性, 但二者并不完全一致: Greenwood相干时间${\tau }_{0}=0.314 {{r}_{0}}/{{V}_{\rm{w}\rm{i}\rm{n}\rm{d}}}$描述湍流大气本身时间起伏, 主要由湍流强度r0和风速$ \nu $决定[28]; 本文定义的大气相干时间$ {\tau }_{\rm{c}} $则由接收光场的时间起伏特性所决定, 不仅受湍流和风速的影响, 还受波长、发射接收参数的影响. 光束参数和激光波长不变的情况下, 接收光场的时间起伏特性与湍流本身时间起伏应具有一致性, 即$ {\tau }_{\rm{c}} $$ {\tau }_{0} $原则上应呈线性关系. 在1.1 km传输距离取不同大气相干长度、风速数值仿真得到$ {\tau }_{\rm{c}} $随大气相干长度r0的变化规律, 如图2所示, 大气相干时间与大气相干长度呈线性关系, 且在一定的大气相干长度下, 随着风速的增大大气相干时间减小. 通过对图2线性拟合可以得出表征接收光场时间特性的大气相干时间${\tau }_{\rm{c}}= 0.45 {{r}_{0}}/{\nu }$. 可以看出$ {\tau }_{\rm{c}} $, $ {\tau }_{0} $均与${{r}_{0}}/{\nu }$呈线性关系, 且二者系数在一般情况下差别较小, 但$ {\tau }_{\rm{c}} $不仅与大气湍流参数相关, 还受激光波长和光束参数的影响, 体现在$ {\tau }_{\rm{c}} $${{r}_{0}}/{\nu }$的系数由激光波长$ \lambda $, 束腰半径$ {w}_{0} $和接收口径$ D $等参数共同决定. 其次, 为定量研究大气湍流与时间相干性之间的变化关系, 下面逐步分析大气参数、发射和接收系统光学参数等对大气相干时间的影响.

      图  2  $\lambda = 532\;{\text{ nm}}$, $L = 1100\;{\text{ m}}$, ${w_0} = 0.2\;{\text{ m}}$, 大气相干时间在不同风速下随大气相干长度变化

      Figure 2.  Atmospheric coherent time varies with atmospheric coherence length in different wind speeds at $\lambda = 532\;{\text{ nm}}$, $L = 1100\;{\text{ m}}$, ${w_0} = 0.2\;{\text{ m}}$.

    • 风速$\nu$和大气相干长度$ {r}_{0} $是决定湍流起伏频率的重要因素, 大气相干长度$ {r}_{0} $描述整个传输路径上湍流扰动的综合强度; 横向风速$\nu$则决定大气湍流的横向“漂移”速度, 二者是影响接收光场大气相干时间$ {\tau }_{\rm{c}} $关键参数. 图3为大气相干时间随大气相干长度与风速的变化. 结果表明, 大气相干时间$ {\tau }_{\rm{c}} $始终随大气相干长度$ {r}_{0} $增加而增加; 大气相干时间$ {\tau }_{\rm{c}} $与风速$\nu$始终呈反比关系, 和Greenwood时间与风速的关系一致.

      图  3  大气相干时间随大气相干长度与风速的变化 (a) L = 1100 m, ${w_0} = 0.1\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$; (b) $L = 1100\;{\text{ m}}$, ${w_0} = 0.1\;{\text{ m}}$, ${r_0} = 0.1\;{\text{ m}}$; (c) $\lambda = 532\;{\text{ nm}}$, ${w_0} = 0.1\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$; (d) $\lambda = 532\;{\text{ nm}}$, ${w_0} = 0.1\;{\text{ m}}$, ${r_0} = 0.1\;{\text{ m}}$; (e) $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, $\nu = 10\;{\text{ }}{{\text{m}} \mathord{\left/ {\vphantom {{\text{m}} {\text{s}}}} \right. } {\text{s}}}$; (f) $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, ${r_0} = 0.1\;{\text{ m}}$

      Figure 3.  Atmospheric coherent time varies with atmospheric coherence length and wind speeds: (a) $L = 1100\;{\text{ m}}$, ${w_0} = 0.1\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$; (b) $L = 1100\;{\text{ m}}$, ${w_0} = 0.1\;{\text{ m}}$, ${r_0} = 0.1\;{\text{ m}}$; (c) $\lambda = 532\;{\text{ nm}}$, ${w_0} = 0.1\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$; (d) $\lambda = 532\;{\text{ nm}}$, ${w_0} = 0.1\;{\text{ m}}$, ${r_0} = 0.1\;{\text{ m}}$; (e) $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, $\nu = 10\;{\text{ }}{{\text{m}} \mathord{\left/ {\vphantom {{\text{m}} {\text{s}}}} \right. } {\text{s}}}$; (f) $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, ${r_0} = 0.1\;{\text{ m}}$.

      图3(a), (c), (e) 可以看出, 不同激光波长$ \lambda $, 传输距离$ L $与束腰半径$ {w}_{0} $下, 大气相干时间$ {\tau }_{\rm{c}} $与大气相干长度$ {r}_{0} $的关系亦有所不同. 图3(a)显示波长$ \lambda $愈长则相干时间$ {\tau }_{\rm{c}} $愈大; 是因为相同湍流强度引起的光程差起伏一致, 波长较长时湍流引起的相位起伏越小, 此时接收光场受湍流扰动也相对较小, 故$ \lambda $越长$ {\tau }_{\rm{c}} $越大, 这也是ISAR相比ISAL几乎不受湍流影响的原因(微波波长远大于激光). 图3(c)表明$ {\tau }_{\rm{c}} $仅与$ {r}_{0} $有关, 与$ L $无关; 这是因为当路径上湍流强度相同时, ${r}_{0}={\big[0.423{C}_{n}^{2}{\left(2\text{π}/\lambda \right)}^{2}L\big]}^{-3/5}$代表整个传输路径上大气湍流效应的强弱, 而$ {C}_{n}^{2} $描述大气湍流本身的折射率起伏情况. 当$ {r}_{0} $不变时, 增大$ L $, $ {C}_{n}^{2} $变小, 使整体$ {r}_{0} $不改变. 因此$ {r}_{0} $不变的情况下改变$ L $$ {\tau }_{\rm{c}} $几乎无影响. 图3(e)则显示$ {\tau }_{\rm{c}} $$ {r}_{0} $的关系随$ {w}_{0} $的增大呈先增大后减小趋势, 即选择合适的$ {w}_{0} $能有效减小湍流影响保持接收光场的时间相干特性. 图3 (b), (d), (f)以横向风速$\nu$为横坐标, 分析不同参数情况下$\nu$对相干时间的影响. 左右两侧图对比可知, 激光波长$ \lambda $, 传输距离$ L $和束腰半径$ {w}_{0} $对大气相干时间的影响规律保持一致, 且相干时间$ {\tau }_{\rm{c}} $正比于相干长度$ {r}_{0} $, 反比于横向风速$\nu$, 分析可知${\tau }_{\rm{c}}\propto {r}_{0}/\nu$, 具体的比例系数受其他参数影响.

    • 大气相干时间$ {\tau }_{\rm{c}} $的定义基于接收光场的时间起伏特性, 故$ {\tau }_{\rm{c}} $不仅受大气参数的影响, 还受激光发射/接收及波长等参数的影响. 3.1节已指出激光波长对大气相干时间的影响, 下面将进一步研究发射参数对大气相干时间的影响.

      图4(a)(c)分别描绘不同$ {r}_{0} $、不同传输距离$ L $及不同接收口径$ D $的情况下, 激光束的束腰半径$ {w}_{0} $对大气相干时间$ {\tau }_{\rm{c}} $的影响规律. 由图4(a)(c)可知, 相干时间$ {\tau }_{\rm{c}} $$ {w}_{0} $的增大而出现先增大后减小的趋势: 当束腰$ {w}_{0} $较小时, 激光发散角较大, 致使接收端只接收到部分接收光场, 故$ {\tau }_{\rm{c}} $较小; 反之当束腰$ {w}_{0} $过大时, 虽激光发散角较小, 但过大束腰半径仍导致接收光场直径较大, 使得$ {\tau }_{\rm{c}} $较小; 图4(c), (d)中显示最优束腰半径w0max随接收口径$ D $的增大而增大, 选择合适的发射和接收参数能够有效降低湍流的退相干效应, 综合来看针对米量级的接收口径、10 km量级的传输距离、中等强度湍流而言, 最优束腰半径在5—15 cm之间.

      图  4  大气相干时间随束腰半径变化  (a) $L = 1100\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$, $\lambda = 532\;{\text{ nm}}$, $D = 0.5\;{\text{ m}}$; (b) $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, $D = 0.5\;{\text{ m}}$; (c) $L = 1100\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$; (d) $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$

      Figure 4.  Atmospheric coherent time varies with beam waist radius: (a) $L = 1100\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$, $\lambda = 532\;{\text{ nm}}$, $D = 0.5\;{\text{ m}}$; (b) $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$, $D = 0.5\;{\text{ m}}$; (c) $L = 1100\;{\text{ m}}$, $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$; (d) $\nu = 10 \;{\rm{m/s}}$, ${r_0} = 0.1\;{\text{ m}}$, $L = 1100\;{\text{ m}}$, $\lambda = 532\;{\text{ nm}}$

    • 接收口径D是影响大气相干时间$ {\tau }_{\rm{c}} $的重要参数,大气相干时间$ {\tau }_{\rm{c}} $随接收口径D变化如图5所示, 可以看出$ {\tau }_{\rm{c}} $D的增加而增加, 当口径D较大时相干时间$ {\tau }_{\rm{c}} $达到饱和不再随$ D $的增加而增加. 通过分析图5(a)可知, 若采用湍流冻结假设当D更大时, 接收端相邻两帧光斑将拥有更多重叠部分, 此时相邻两帧间的电场相干性更强使得$ {\tau }_{\rm{c}} $较大. 而光斑尺寸与接收口径大小相当则光束全接收, 此时$ {\tau }_{\rm{c}} $D呈先增大后不变的变化趋势. 图5(b)表明, 不同$ {w}_{0} $下, D$ {\tau }_{\rm{c}} $的影响均与图5(a)中保持一致. 因此, 采用大口径接收端有助于提升接收信号的相位稳定性和时间相干性, 这也是在ISAL中大口径望远镜能够获得比小口径望远镜更高分辨率图像的原因.

      图  5  $\nu = 10 \;{\rm{m/s}} , L=1100\;\text{ m }, \lambda =532\;\text{ nm}$, 大气相干时间随接收口径的变化 (a) ${w_0} = 0.1\;{\text{ m}}$; (b) ${r_0} = 0.1\;{\text{ m }}$

      Figure 5.  Atmospheric coherent time varies with receiving aperture at $\nu = 10 \;{\rm{m/s}} , L=1100\;\text{ m }, \lambda =532\;\text{ nm}$: (a) ${w_0} = 0.1\;{\text{ m}}$; (b) ${r_0} = 0.1\;{\text{ m }}$

    • 大气湍流引起的接收信号随时间的相位起伏, 直接影响相干探测信号的稳定性, 这种不稳定性为ISAL信号处理带来相位误差, 导致成像性能降低. 而受湍流影响后的接收信号相位变化的快慢体现在其大气相干时间上, 保证成像时间在大气相干时间内能够有效抑制产生较大的相位误差, 有助于提高相干探测性能以及ISAL成像质量.

      本文采用多层动态相位屏方法仿真高斯光束在大气湍流中的传输过程, 并定义大气相干时间表征接收光场的时间起伏特性, 进一步研究大气参数、发射/接收参数、激光波长等对大气相干时间的影响. 研究结论包括以下三方面: 1)大气相干时间与大气相干长度呈正相关, 与风速呈反比关系, 和Greenwood时间与风速的关系一致; 2)大气相干时间与激光波长呈正相关; 在大气相干长度不变的情况下(即路径上综合湍流强度不变), 大气相干时间不随传输距离发生变化; 3)大气相干时间随束腰半径的增大出现先增大后减小的趋势, 最优束腰半径随接收口径的增大而增大, 选择合适的发射和接收参数能够有效降低湍流的退相干效应; 大气相干时间与接收口径呈正相关, 当接收口径较大时大气相干时间达到饱和不再随接收口径的增加而增加, 故适当采用大口径接收望远镜有助于提高接收信号的时间相干性. 本文的研究结果为评估大气湍流对相干探测过程的影响提供有力参考, 并在ISAL成像过程中为优化参数提供可靠的参考价值.本文仅考虑影响大气相干时间的湍流和光学参数, 并未考虑如何消除湍流引起的接收相位误差, 对于实际情况的研究将在以后的工作中具体展开.

参考文献 (28)

目录

    /

    返回文章
    返回