-
在激光相干探测中,目标反射光受大气湍流和目标散射特性的双重影响。大气湍流和目标散射特性是影响接收光强特性和空间相干特性的重要因素。而光场的相关长度(横向相干长度)反应的是激光的空间相干性,开展大气湍流传输条件下的相干长度变化的研究有助于理解激光传输与接收光场的空间相干特性变化。
已有大量文献研究了目标的散射特性和大气湍流中的空间相干特性[1-4]。1963年,Goodman提出了服从随机高斯分布的目标散射模型[5]。2003年,Korotkova等人采用惠更斯-菲涅尔原理和Rytovy原理,并结合光学参量法得到了高斯光束的随机粗糙目标的散射光互相干函数,并研究了双基激光雷达系统中随机电磁光束经目标散射后的平均光强和偏振特性[4, 6]。目前对目标散射特性的分析主要针对回波光强和偏振特性,针对目标散射特性对空间相干特性影响的研究较少。
Richlin等人运用交叉谱密度函数研究了部分相干光在大气湍流中的光强闪烁和孔径平滑因子[7- 8]。Wu等人研究比较了GSM (Gaussian Schell-Mode)光束在真空和湍流大气中的横向相干长度等传输特性,得到了“与完全相干光相比,部分相干光受大气湍流影响更小”的结论[9]。Friberg研究了部分相干光在自由空间中的光束束宽和波前曲率半径特性,并分析了空间相干特性[10]。王华和李成强课题组研究了高斯-谢尔模光束在大气湍流中的空间相干长度特性[11]。
然而,以上有关空间相干特性的分析主要是基于单程大气湍流传输开展的。对于激光相干探测中目标反射参数对接收光强和相干长度的影响仍需进一步研究和分析。本文利用广义惠更斯-菲涅尔原理结合随机高斯分布的目标散射理论,推导得到了大气湍流中激光探测目标回波光信号的交叉谱密度函数解析表达式,进一步得到激光探测回波信号光的光强分布和相干长度的表达式,并进行了不同光源参数和不同湍流条件下的数值模拟分析。
-
在发射平面
${\textit{z}} = 0$ 处,部分相干高斯光束的交叉谱密度函数如下[12-13]:$$\begin{split} {W_0}\left( {{{r}_1},{{r}_2},0} \right) =& \left\langle {{U_0}\left( {{{r}_1},0} \right)U_0^ * \left( {{{r}_2},0} \right)} \right\rangle \\ =& {I_0}\exp \left[ \!{ - \frac{1}{{w_0^2}}\left( {{r}_1^2 + {r}_2^2} \right)}\! \right] \exp \left[\! { - \frac{1}{{2\sigma _0^2}}{{\left( {{{r}_1} - {{r}_2}} \right)}^2}}\! \right] \end{split} ,$$ (1) 其中,
${{r}_1}$ 和${{r}_2}$ 表示垂直于发射光传输方向的横截面矢量,${U_0}\left( {{{r}_1},0} \right)$ 和${U_0^*}\left( {{{{r}}_2},0} \right)$ 表示发射平面处的随机光场,${I_0}$ 表示光强,${w_0}$ 表示发射光束束腰半径,${\sigma _0}$ 表示光源部分相干长度。根据惠更斯-菲涅尔衍射原理,部分相干光经过大气湍流远距离传输后,到达目标处的交叉谱密度函数可定义为[14]:
$$\begin{split} & W\left( {{{\rho }_1},{{\rho }_2},{\textit{z}}} \right) \\ & = \frac{1}{{{\lambda ^2}{{\textit{z}}^2}}} {\iint {{W_0}\left( {{{r}_1},{{r}_2},0} \right){X_1}{Y_1}{\rm{d}}{{r}}_1^2{\rm{d}}{{r}}_2^2}}, \end{split}$$ (2) 其中:
$${X_1} = \left\langle {\exp \left[ {\phi \left( {{{r}_1},{{\rho }_1}} \right) + {\phi }\left( {{{r}_2},{{\rho }_2}} \right)} \right]} \right\rangle ,$$ $${Y_1} = \exp \left\{ {\frac{{jk}}{{2{\textit{z}}}}\left[ {{{\left( {{{\rho }_1} - {{r}_1}} \right)}^2} + {{\left( {{{\rho }_2} - {{r}_2}} \right)}^2}} \right]} \right\}\;.$$ 式中,
${{\rho }_1}$ 和${{\rho }_2}$ 表示目标处垂直于入射光的横截面矢量,z表示传输距离,$\phi \left( {{{r}_1},{{\rho }_1}} \right)$ 和$\phi \left( {{{r}_2},{{\rho }_2}} \right)$ 表示随机相位变化,部分相干光的复相位互相关函数${X_1}$ 可近似为[11]:$$\begin{split} {X_1} \approx \exp \left\{ { - M\left[ {{\rho }_c^2 + {{\rho }_c}{r_c} + {{r}}_c^2} \right]} \right\} \end{split}, $$ (3) 其中,
${{\rho }_c} = {{\rho }_1} - {{\rho }_2}$ ,${{r}_c} = {{r}_1} - {{r}_2}$ ,M 表示与大气湍流相干长度相关的参数,$M = 0.33{\textit{z}}{{\text π}^2}{k^2}\displaystyle\int_0^\infty$ ${{\kappa ^3}{{\varPhi }_n}\left( \kappa \right)} {\rm d}\kappa$ ,${\varPhi _n}\left( \kappa \right) = 0.033C_n^2{\kappa ^{{{ - 11} / 3}}}{{\rm{e}} ^{{{ - {\kappa ^2}} / {\kappa _m^2}}}}$ ,表示大气折射率功率函数,$\kappa \geqslant 1/{l_0}$ ,${\kappa _m} = 5.92/{l_0}$ ,${l_0}$ 表示大气湍流的内尺度,$C_n^2$ 表示大气折射率结构常数,其衡量的是湍流强度大小,该值越大,表示湍流强度越强。$$\begin{split} & {{r}_c} = \frac{{{{{r}}_1} + {{r}_2}}}{2}\quad ,\quad {{r}_d} = {{r}_1} - {{r}_2} ,\\ & {{\rho }_c} = \frac{{{{\rho }_1} + {{\rho }_2}}}{2},\quad {{\rho }_d} = {{\rho }_1} - {{\rho }_2} . \end{split} $$ (4) 经过积分运算,可以得到部分相干光经过大气湍流传输后,目标处光场的交叉谱密度函数[14]:
$$\begin{split} W\left( {{{\rho }_1},{{\rho }_2},{\textit{z}}} \right) =\; & \frac{{{I_0}}}{{L({\textit{z}})}}\exp \left( { - {X_2}{{\rho }}_d^2} \right) \cdot\\ & \exp \left[ { - \frac{{{\rho }_c^2}}{{2w_0^2L\left({\textit{z}} \right)}}} \right]\exp \left[ { - \frac{{ik{{\rho }_c}{{\rho }_d}}}{{R\left( {\textit{z}} \right)}}} \right], \end{split}$$ (5) 其中,w0为束腰半径;
$${X_2} = \frac{1}{{2\delta _0^2L\left({\textit{z}}\right)}} + M\left( {1 + \frac{2}{{L\left( {\textit{z}}\right)}}} \right) - \frac{{2{M^2}{{\textit{z}}^2}}}{{{k^2}w_0^2L\left( z \right)}};$$ $$R\left({\textit{z}} \right) = \frac{{{\textit{z}}{k^2}w_0^2L\left({\textit{z}} \right)}}{{{k^2}w_0^2L\left({\textit{z}}\right) + M{{\textit{z}}^2} - {k^2}w_0^2}};$$ $$L\left({\textit{z}} \right) = 1 + \frac{{{{\textit{z}}^2}}}{{{k^2}w_0^2\delta _0^2}} + \frac{{2M{{\textit{z}}^2}}}{{{k^2}w_0^2}};$$ $$\frac{1}{{{\delta _0}}} = \frac{1}{{4w_0^2}} + \frac{1}{{\sigma _0^2}}.$$ 当
$\,{{\rho }_1} =\, {{\rho }_2} =\, {\rho}$ 时,$\,{{\rho }_d} = 0$ ,$\,{{\rho} _c} = {\rho}$ ,由交叉谱密度函数公式可以得到目标处的入射光平均光强为:$$\begin{split} & \left\langle {I\left( {{\rho} ,{\textit{z}}} \right)} \right\rangle = W\left( {{\rho} ,{\textit{z}}} \right) \\ & = \frac{{{I_0}}}{{L\left({\textit{z}}\right)}}\exp \left[ { - \frac{{2{{\rho} ^2}}}{{w_0^2L\left({\textit{z}}\right)}}} \right]. \end{split}$$ (6) -
根据Goodman提出的服从高斯随机分布的粗糙目标表面散射模型[5],目标反射光场与入射光场的关系可表示为:
$${U_R}\left( {{\rho}} \right) = F{\alpha _T}{U_{\textit{z}}}\left( {{\rho}} \right)\exp \left[ { - i{\phi _T}\left( {{\rho}} \right)} \right],$$ (7) 其中,
$F$ 表示目标表面的反射率,${U_{\textit{z}}}\left( {{\rho}}\right)$ 表示目标处的入射光场,${\alpha _T} = \left[ {1 + {4 / {\left( {{k^2}l_c^2} \right)}}} \right]$ ,表示目标表面粗糙度的函数,${l_c} = {\sigma _H}/{\delta _\phi }$ ,${\sigma _H}$ 表示随机起伏表面相关长度,$\delta _\phi ^2 = {\left[ {k\left( {1 + \cos \beta } \right)} \right]^2}\delta _H^2$ ,表示粗糙表面引起的相位起伏方差,$\delta _H^2$ 表示粗糙面高度起伏方差,$\beta $ 是入射角。${\phi _T}\left( {{\rho}} \right) = k\left( {1 + \cos \beta } \right)h\left( {{\rho}} \right)$ ,表示目标反射引起的光场随机相位因子,$h\left({{\rho}} \right)$ 表示目标表面任意一点相对于参考平面的高度函数。对于服从高斯分布的随机粗糙目标,当入射光垂直照射目标时,目标反射光场的互相关函数为:$$\begin{split} & \left\langle {{U_R}\left( {{{\rho} _1}} \right)U_R^ * \left( {{{\rho} _2}} \right)} \right\rangle \\ & = {F^2}\alpha _T^2{X_3}\exp \left[ { - \frac{{\delta _\phi ^2}}{{\sigma _H^2}}{{\left( {{{\rho} _1} - {{\rho} _2}} \right)}^2}} \right] \end{split}, $$ (8) 其中,
${X_3} = \left\langle {{U_{\textit{z}}}\left( {{{\rho} _1}} \right)U_{\textit{z}}^ * \left( {{{\rho} _2}} \right)} \right\rangle$ 。当目标为漫反射物体时${l_c} \to 0$ ,目标反射光场的互相干函数可表示为:$$\begin{split} & \left\langle {{U_R}\left( {{{\rho} _1}} \right)U_R^ * \left( {{{\rho} _2}} \right)} \right\rangle \\ & = \frac{{4{\text π} {F^2}}}{{{k^2}}}\left\langle {{U_{\textit{z}}}\left( {{{\rho} _1}} \right)U_{\textit{z}}^ * \left( {{{\rho} _2}} \right)} \right\rangle \delta \left( {{{\rho} _1} - {{\rho} _2}} \right) \end{split}. $$ (9) 结合探测光在目标处的入射光强公式,可得到探测光经过粗糙目标漫反射后的平均光强为:
$${I_T}\left({\textit{z}}\right) = \frac{{4{\text π} {I_0}{F^2}}}{{{k^2}L\left({\textit{z}}\right)}}\exp \left[ { - \frac{{2r_T^2}}{{w_0^2L\left({\textit{z}}\right)}}} \right],$$ (10) 其中,
${r_{\rm{T}}}$ 表示目标偏离光轴中心的距离,当${r_{\rm{T}}} = 0$ 时,式(10)表示目标刚好位于光轴中心位置。 -
设反射光的束腰半径为
${w_T} = \alpha {w_0}$ ,$\alpha $ 表示束腰半径变化参数:当目标为漫反射时$\alpha > 1$ ;当目标为均匀的镜面反射时$\alpha \cong 1$ 。可以近似得到目标反射光场函数和空间相干度分布函数:$${I_T}\left( {{{\rho} _1},{{\rho} _2},0} \right) = {I_T}\left({\textit{z}}\right)\exp \left[ { - \frac{{{\rho} _1^2 + {\rho} _2^2}}{{w_T^2}}} \right],$$ (11) $$\mu \left( {{{\rho} _1},{{\rho} _2},0} \right) = \exp \left[ { - \frac{{{{\left( {{{\rho} _1} - {{\rho} _2}} \right)}^2}}}{{2{\sigma ^2}\left( {\rho} \right)}}} \right],$$ (12) 式中,
${I_T}\left({\textit{z}}\right)$ 表示反射光振幅,$\sigma \left({ \rho} \right)$ 表示反射光的相干长度。再次利用惠更斯-菲涅尔原理,可以得到探测系统接收平面上光场空间频域的交叉谱密度函数:$$\begin{split} & W\left( {{{s}_1},{{s}_2},{\textit{z}}} \right) = \frac{{{k^2}}}{{4{{\text π}^2}{{\textit{z}}^2}}} .\\ & {\iint {{X_4}{Y_4}{W_T}\left( {{{\rho} _1},{{\rho} _2},0} \right){\rm{d}}{{s}}_1^2{\rm{d}}{{s}}_2^2}} , \end{split}$$ (13) 式中:
$${X_4} = \left\langle {\exp \left[ {\phi \left( {{{\rho} _1},{{s}_1}} \right) + {\phi ^ * }\left( {{{\rho} _2},{{s}_2}} \right)} \right]} \right\rangle ;$$ $${Y_4} = \exp \left\{ {\frac{{jk}}{{2{\textit{z}}}}\left[ {{{\left( {{{s}_1} - {{\rho} _1}} \right)}^2} + {{\left( {{{s}_2} - {{\rho} _2}} \right)}^2}} \right]} \right\}.$$ ${{{s}}_1}$ 和${{{s}}_2}$ 表示接受光场的横截面矢量。$$\begin{split} & {W_T}\left( {{{\rho} _1},{{\rho} _2},0} \right) = {I_T}\left({\textit{z}}\right) \times \\ & \exp \left[ { - \frac{{{\rho} _1^2 + {{\rho}} _2^2}}{{w_T^2}}} \right]\exp \left[ { - \frac{{{{\left( {{{\rho} _1} - {{\rho} _2}} \right)}^2}}}{{2{\sigma ^2}\left( {\rho} \right)}}} \right]\text{.} \end{split} $$ (14) 目标反射光束与入射光束传输经历的大气湍流环境相同,所以,假设反射光的复相位互相关函数与入射光的复相位互相关函数一样。利用相似的积分运算过程,得到接收平面上光场的交叉谱密度函数:
$$\begin{split} & W\left( {{{s}_1},{{s}_2},{\textit{z}}} \right) = \frac{{{I_T}\left({\textit{z}}\right)}}{{{L_2}({\textit{z}})}}\exp \left\{ { - {X_5}{{\left( {{{s}_1} - {{s}_2}} \right)}^2}} \right\} .\\ & \exp \left[ { - \frac{{{{\left( {{{s}_1} + {{s}_2}} \right)}^2}}}{{2w_T^2{L_2}\left({\textit{z}}\right)}}} \right]\exp \left[ { - \frac{{ik\left( {{s}_1^2 - {{s}}_2^2} \right)}}{{2{R_2}\left({\textit{z}}\right)}}} \right] , \end{split}$$ (15) 式中
$${X_5} = \frac{1}{{2\delta _1^2{L_2}\left({\textit{z}}\right)}} + M\left( {1 + \frac{2}{{{L_2}\left({\textit{z}}\right)}}} \right) - \frac{{2{M^2}{{\textit{z}}^2}}}{{{k^2}w_T^2{L_2}\left({\textit{z}}\right)}};$$ $${L_2}\left({\textit{z}}\right) = 1 + \frac{{{{\textit{z}}^2}}}{{{k^2}w_T^2\delta _1^2}} + \frac{{2M{{\textit{z}}^2}}}{{{k^2}w_T^2}}\text{;}$$ $$\frac{1}{{{\delta _1}}} = \frac{1}{{4w_T^2}} + \frac{1}{{{\sigma ^2}\left( {\rho} \right)}};$$ $$\sigma \left( {{\rho}} \right)\; = \frac{1}{{\dfrac{1}{{2\sigma _0^2L\left({\textit{z}}\right)}} + M\left( {1 + \dfrac{2}{{L\left({\textit{z}}\right)}}} \right) - \dfrac{{2{M^2}{{\textit{z}}^2}}}{{{k^2}w_0^2L\left({\textit{z}}\right)}}}};$$ $${R_2}\left({\textit{z}}\right) = \frac{{{\textit{z}}{k^2}w_T^2{L_2}\left({\textit{z}}\right)}}{{{k^2}w_T^2{L_2}\left({\textit{z}}\right) + M{{\textit{z}}^2} - {k^2}w_T^2}}.$$ 由公式(15)可以得到接收平面上的平均接收光强函数:
$$\begin{split} \left\langle {I\left( {{s},{\textit{z}}} \right)} \right\rangle =& W\left( {{s},{\textit{z}}} \right) \\ =& \frac{{{I_T}\left({\textit{z}}\right)}}{{{L_2}({\textit{z}})}}\exp \left[ { - \frac{{{{s}^2}}}{{2w_T^2{L_2}\left({\textit{z}}\right)}}} \right], \end{split}$$ (16) 由公式(16)可得,本文得到的目标反射光强公式与向宁静[15-16]推导得到的结果有所不同,本文得到目标反射光强分布公式与光源参数和湍流相关。为了分析湍流大气传输对目标反射光相干特性的影响,本文利用复相干度来分析接收光场的空间相干特性[9, 17]:
$$\mu \left( {{{s}_1},{{s}_2},\omega } \right) = \frac{{W\left( {{{s}_1},{{s}_2},\omega } \right)}}{{\sqrt {W\left( {{{s}_1},{{s}_1},\omega } \right)} \sqrt {W\left( {{{s}_2},{{s}_2},\omega } \right)} }}.$$ (17) 考虑实际工程应用,研究目标反射光场任意一点与轴上一点之间的空间相干度,令
${{s}_1} = {s}$ ,${{s}_2} = 0$ ,将反射光功率谱密度函数带入式(16),整理得:$$\mu \left( {{s},\omega } \right) = \exp \left( { - {X_6}{{s}^2}} \right)\exp \left[ { - \frac{{ik}}{{2{R_2}\left({\textit{z}}\right)}}{{s}^2}} \right],$$ (18) 式中
$${X_6} = \frac{1}{{2\delta _1^2{L_2}\left({\textit{z}}\right)}} + M\left( {1 + \frac{2}{{{L_2}\left({\textit{z}}\right)}}} \right) - \frac{{2{M^2}{{\textit{z}}^2}}}{{{k^2}w_T^2{L_2}\left({\textit{z}}\right)}}.$$ 由式(18)可得,目标反射光的空间相干特性与探测光的光源参数、自由空间衍射效应及大气湍流效应有关,但与目标的平均反射率无关。
根据光束束宽的定义,可以近似定义光束相干长度,反射光束相干长度的平方定义为:
$${\sigma ^2}\left( {{s},{\textit{z}}} \right) = \frac{{\displaystyle{\iint_{ - \infty }^\infty {{{s}^2}\exp \left( { - {X_6}{{s}^2}} \right){\rm{d}^2}{s}} } }}{{\displaystyle{\iint_{ - \infty }^\infty {\exp \left( { - {X_6}{{s}^2}} \right){\rm{d}^2}{s}} } }}.$$ (19) 将空间相干度的幅值带入式(19),经过复杂积分运算可得到目标反射光相干长度的平方:
$${\sigma ^2}\left( {{s},{\textit{z}}} \right) = \frac{1}{{{X_6}}}.$$ (20) 由公式(20)可以看出,目标反射光的相干长度与光源参数、湍流强度、传输距离和目标反射特性相关。本文得到的目标反射光相干长度公式与Friberg和李成强等人得到的单程大气传输的相干长度公式形式上一致。
-
根据理论分析得到的式(16)和式(20),本文对不同光源参数、不同目标反射参数和不同湍流强度下的光强分布和相干长度进行数值仿真,数值仿真中的参数设置如表1所示。数值模拟过程中的激光线宽在几十兆赫兹到上百兆赫兹,与实际应用中的激光线宽相符。
表 1 数值模拟参数
Table 1. Parameters of numerical simulation
系统参数 符号 数值 激光波长 $\lambda $ 532 nm、1 064 nm、1 550 nm 束腰半径 ${w_0}$ 3 cm 湍流内尺度 ${l_0}$ 1~5 cm 束腰半径与相干长度比 ${{{w_0}}/ {{\sigma _0}}}$ 10−4~0.1 大气结构常数 $C_n^2$ $0 : {10^{ - 15} }{\rm{m}^{ { { - 2} / 3} } }$ 目标反射率 F 0~1 首先,分析波长分别为532 、1 064 和1 550 nm时,光束在大气湍流传输中的归一化光强随传输距离变化的情况,目标反射率为0.3,
$C_n^2 =1.7 \times {10^{ - 17}}$ 和$C_n^2 = 1.2 \times {10^{ - 16}}$ 时的归一化光强变化曲线如图1(a)和1(b)所示。由图1可得,光源波长越长归一化光强值越大,作用距离越远,而且随着湍流强度的增大,归一化光强值下降速度加快。图 1 不同波长的归一化光强随传输距离的变化曲线
Figure 1. Curves of normalized intensity varying with transmission distance with different wavelengths
当波长为532 nm,目标反射率
$F = 0.3$ ,z = 1 000 m,反射光与光轴的距离${r_T} = 0$ 时,在不同湍流强度、不同束腰半径与相干长度比的情况下,对接收平面上的接收光强进行数值分析,采用$C_n^2 = 0$ ,${{{w_0}} / {{\sigma _0}}} = 0.4 \times {10^{ - 3}}$ 时的中心光强对接收平均光强进行归一化处理,接收平面上的归一化光强分布模拟结果如图2所示。由图2(a)和2(b)可以看出,在大气湍流传输中,不同光斑尺度与相干长度比值${{{w_0}} / {{\sigma _0}}}$ 下的归一化光强间的差别较小,同时,接收平面上的光斑展宽所受影响也较小。比较图2(a)和2(b)可知,在传输距离和光源参数固定的条件下,随着湍流强度的增加,接收平面上的光束半径展宽明显,同时平均接收光强下降明显,湍流强度越大,光束展宽和光强下降越显著。从上述结果可以看出,激光线宽对接收光强和光斑展宽的影响较小。图 2 不同湍流强度和
${{{{w_0}} / {{\sigma _0}}}}$ 时,接收平面上的归一化光强分布Figure 2. Normalized intensity distributions at the receiving plane with different
$C_n^2$ and${{{{w_0}} / {{\sigma _0}}}}$ 在发射光源参数确定的情况下,不同湍流强度,不同反射光束腰半径与光源束腰半径比下,中心接收光强随探测距离的变化如图3所示。由图3(a)可看出,在弱湍流近距离传输时,中心光强曲线下降缓慢,随着传输距离的增加,中心光强变化曲线呈线性下降趋势。随着湍流强度的增大,如图3(b)、3(c)、3(d)所示,光强曲线下降速度加快,当下降到一定程度后保持稳定。这是因为在弱湍流区,光束传输过程中受湍流影响较小,而受光源束宽和相干长度的影响较大。当湍流强度增大时,光束受大气湍流影响增大,受光源参数的影响相对变小,大气湍流对光强损耗的影响占主导地位。
图 3 不同湍流强度和
${{{w_T}} / {{w_0}}}$ 时,中心接收光强随传输距离z的变化曲线Figure 3. Center average received light intensity varying with transmission distance with different
$C_n^2 $ and${{{w_T}} / {{w_0}}}$ 图4(a)和4(b)是光源波长分别为532 、1 064 和1 550 nm时,不同湍流强度下接收光场的相干长度变化曲线。由数值模拟结果可以看出,其他参数一致的情况下,激光波长越长,相干长度值越大。比较图4(a)和4(b)可知,随着湍流强度的增大,3种波长接收光场的相干长度值均明显下降,在强湍流区接收光场的相干长度变化曲线近似缓慢地线性增加。
图 4 不同湍流强度下3种波长接收光场的相干长度变化曲线
Figure 4. Coherent length curves of the received light field varying with the wavelength and
$C_n^2 $ 图5表示大气湍流中不同光源束腰参数对接收光场空间相干长度的影响,仿真波长为532 nm。由图5(a)可得,在弱湍流条件下,接收光场的相干长度曲线随着传输距离的增大迅速增加,且光源束腰参数越小相干长度曲线增加越快。图5(b)表示大气湍流中不同光源束腰半径接收光场的相干长度变化曲线。由图5(b)可得,在强湍流传输过程中,随着传输距离的增大相干长度曲线先增加而后缓慢下降,然后再增加,即出现了相干长度随传输距离变化的峰值,且光源束腰半径越小峰值现象越明显,光源束腰值越大相干长度曲线增加越缓慢。
图 5 不同光源束宽和湍流强度时,目标反射光的相干长度变化曲线
Figure 5. Coherent length curves of target reflected light varying with the beam width and turbulence intensity
图6表示大气湍流中光源参数一定时,反射光束腰半径变化对接收光相干长度的影响,仿真波长为532 nm。由图6(a)可以看出,随着传输距离的增加接收光的相干长度近似呈线性递增,且在相同传输距离下,目标反射光的束腰半径越大反射光的相干长度越小。图6(b)为大气湍流中的相干长度变化曲线。由图6(b)可以看出,接收光相干长度在近距离传输时先呈线性增加状态,且反射光束腰半径越大相干长度值越小,这与自由空间中的特性相似,随着在大气湍流中传输距离的增大,相干长度缓慢下降,反射光不同束腰参数的相干长度值最终汇聚在一起。
图 6 不同湍流强度下,目标反射光与光源不同束腰比时,目标反射光的空间相干长度变化曲线
Figure 6. Spatial coherence length curves of the reflected light varying with beam waist ratio under different
$C_n^2 $ 产生以上结果的物理解释如下:在弱湍流大气传输中,大气湍流影响较弱,大气湍流对光束展宽的影响可以忽略,接收光场的相干长度主要受光源参数和目标反射参数的影响。光束在大气湍流中传输时,受大气湍流的影响,光束波前相位随机增加,湍流强度越大,湍流大气对波前的影响越大。光束在大气湍流中传输时,随着传输距离的增加,大气湍流线性叠加。当传输距离较近时,光束受大气湍流的影响较小,光源参数和目标反射参数对相干长度的影响占主导地位。随着传输距离的增加,大气湍流对光束的影响呈现线性叠加,此时,大气湍流对接收光场相干长度的影响占主导地位,而光源参数和目标反射参数的影响较弱。因此随着传输距离的增加,不同光源参数和反射光参数的相干长度曲线会聚一起。
-
本文利用广义惠更斯-菲涅尔原理并结合目标散射理论模型,以部分相干光为例,得到了激光探测接收光场的交叉谱密度函数表达式,并得到了激光探测的光强函数和空间相干长度表达式。利用光强函数和相干长度表达式,分析了光源参数和大气湍流对光强分布和相干长度的影响。结果表明:接收光强和相干长度受光源参数、大气湍流和目标反射光参数综合影响,在弱湍流大气传输中,光强分布和相干长度主要受光源参数的影响;在强湍流大气传输中,大气湍流对接收光强的影响占主导地位,光源参数和目标反射光参数的影响减弱,接收光强下降,同时相干长度下降;目标反射光半径对接收光强分布和相干长度的变化影响较大。以上目标探测中接收光强和相干长度的分析可以用于激光相干探测领域,对分析激光相干探测性能具有重要意义。
Light intensity and spatial coherence characteristics of laser coherent detection in a turbulent atmosphere
-
摘要: 本文利用广义惠更斯-菲涅尔原理结合Goodman目标散射理论,推导出了激光探测中目标反射光的交叉谱密度函数解析式,并进一步得到了目标反射光的光强分布和空间相干长度表达式。利用得到的表达式在湍流大气条件下,分析了不同光源参数和目标反射光参数对目标反射光光强分布和相干长度的影响。研究结果表明:光源相干长度对归一化反射光强影响较小;光源束腰半径和反射光斑半径值越大接收光的相干长度值越小;随着传输距离的增加,相干长度增加越来越缓慢;在弱湍流大气传输过程中,光源参数对接收光的影响占主导作用,光源束腰半径越大接收光的光强和相干长度值越小;在强湍流大气传输过程中,大气湍流对反射光的影响起主导作用。Abstract: In this paper, the cross-spectral density function of target reflected light in laser detection is obtained by using generalized Huygens-Fresnel principle and Goodman target scattering theory. On the basis of above, the expression of intensity distribution and spatial coherence length of target reflected light is derived. The influence of different light source and target reflected light parameters on the intensity distribution and coherent length of the target reflected light is simulated by using the expressions obtained in this paper under turbulent atmospheric conditions. The results show that the coherence length of the light source has little effect on the normalized light intensity distribution; the coherence length of the received light is smaller with a larger beam waist radius and reflected light radius, and the coherence length increases at a slower rate as the transmission distance increases. In the process of weak turbulent atmospheric transmission, the influence of light source parameters on the received light is much stronger. The larger the beam waist radius, the smaller the received light intensity and coherence length value. During strong turbulent atmospheric transmission, the influence of atmospheric turbulence on the received light is dominant.
-
Key words:
- laser detection /
- atmospheric turbulence /
- intensity distribution /
- coherence length
-
表 1 数值模拟参数
Table 1. Parameters of numerical simulation
系统参数 符号 数值 激光波长 $\lambda $ 532 nm、1 064 nm、1 550 nm 束腰半径 ${w_0}$ 3 cm 湍流内尺度 ${l_0}$ 1~5 cm 束腰半径与相干长度比 ${{{w_0}}/ {{\sigma _0}}}$ 10−4~0.1 大气结构常数 $C_n^2$ $0 : {10^{ - 15} }{\rm{m}^{ { { - 2} / 3} } }$ 目标反射率 F 0~1 -
[1] LU W, LIU L R, SUN J F, et al. Change in degree of coherence of partially coherent electromagnetic beams propagating through atmospheric turbulence[J]. Optics Communications, 2007, 271(1): 1-8. doi: 10.1016/j.optcom.2006.09.058 [2] 周鑫, 姜鹏, 孙剑峰, 等. 基于点目标大气闪烁的目标回波分布研究[J]. 红外与激光工程,2017,46(S1):74-81. ZHOU X, JIANG P, SUN J F, et al. Investigation on the distribution of target echo based on point target atmospheric scintillation[J]. Infrared and Laser Engineering, 2017, 46(S1): 74-81. (in Chinese) [3] COLLETT E, WOLF E. Beams generated by Gaussian quasi-homogeneous sources[J]. Optics Communications, 1980, 32(1): 27-31. doi: 10.1016/0030-4018(80)90307-7 [4] KOROTKOVA O, ANDREWS L C, PHILLIPS R L. Lidar model for a rough-surface target: method of partial coherence[J]. Proceedings of SPIE, 2004, 5237: 49-60. doi: 10.1117/12.515086 [5] GOODMAN J W. Statistical properties of laser speckle patterns[M]. DAINTY J C. Laser Speckle and Related Phenomena. Berlin, Heidelberg: Springer, 1975. [6] KOROTKOVA O, CAI Y, WATSON E. Stochastic electromagnetic beams for LIDAR systems operating through turbulent atmosphere[J]. Applied Physics B, 2009, 94(4): 681-690. doi: 10.1007/s00340-009-3404-4 [7] RICKLIN J C, DAVIDSON F M. Atmospheric turbulence effects on a partially coherent Gaussian beam: implications for free-space laser communication[J]. Journal of the Optical Society of America A, 2002, 19(9): 1794-1802. doi: 10.1364/JOSAA.19.001794 [8] RICKLIN J C, DAVIDSON F M. Atmospheric optical communication with a Gaussian Schell beam[J]. Journal of the Optical Society of America A, 2003, 20(5): 856-866. doi: 10.1364/JOSAA.20.000856 [9] WU J, BOARDMAN A D. Coherence length of a Gaussian-Schell beam and atmospheric turbulence[J]. Journal of Modern Optics, 1991, 38(7): 1355-1363. doi: 10.1080/09500349114551521 [10] FRIBERG A T, SUDOL R J. Propagation parameters of Gaussian Schell-model beams[J]. Optics Communications, 1982, 41(6): 383-387. doi: 10.1016/0030-4018(82)90161-4 [11] LI M N, TAN L Y, MA J, et al. Performance analysis of a free-space laser communication system with a Gaussian Schell model[J]. Journal of Modern Optics, 2015, 62(19): 1608-1615. doi: 10.1080/09500340.2015.1054907 [12] 柯熙政, 王婉婷. 部分相干光在大气湍流中的光束扩展及角扩展[J]. 红外与激光工程,2015,44(9):2726-2733. doi: 10.3969/j.issn.1007-2276.2015.09.032 KE X ZH, WANG W T. Expansion and angular spread of partially coherent beam propagating in atmospheric turbulence[J]. Infrared and Laser Engineering, 2015, 44(9): 2726-2733. (in Chinese) doi: 10.3969/j.issn.1007-2276.2015.09.032 [13] WU ZH S, LI Y Q. Scattering of a partially coherent Gaussian-Schell beam from a diffuse target in slant atmospheric turbulence[J]. Journal of the Optical Society of America A, 2011, 28(7): 1531-1539. doi: 10.1364/JOSAA.28.001531 [14] 李成强, 张合勇, 王挺峰, 等. 高斯-谢尔模光束在大气湍流中传输的相干特性研究[J]. 物理学报,2013,62(22):224203. doi: 10.7498/aps.62.224203 LI CH Q, ZHANG H Y, WANG T F, et al. Investigation on coherence characteristics of Gauss-Schell model beam propagating in atmospheric turbulence[J]. Acta Physica Sinica, 2013, 62(22): 224203. (in Chinese) doi: 10.7498/aps.62.224203 [15] LI M N, TAN L Y, MA J, et al. Statistical distribution of the optical intensity obtained using a Gaussian Schell model for space-to-ground link laser communications[J]. Journal of Modern Optics, 2016, 63(10): 921-931. doi: 10.1080/09500340.2015.1111452 [16] 向宁静, 王明军, 王太荣. 部分相干高斯-谢尔光束在大气湍流中的平均强度与展宽[J]. 激光杂志,2012,33(5):4-6. doi: 10.3969/j.issn.0253-2743.2012.05.003 XIANG N J, WANG M J, WANG T R. Average intensity and spreading of a partially coherent Gaussian Schell-model beam propagation through atmospheric turbulence[J]. Laser Journal, 2012, 33(5): 4-6. (in Chinese) doi: 10.3969/j.issn.0253-2743.2012.05.003 [17] 埃米尔·沃尔夫. 光的相干与偏振理论导论[M]. 蒲继雄, 译. 北京: 北京大学出版社, 2014. WOLF E. Introduction to the Theory of Coherence and Polarization of light[M]. PU J X, trans. Beijing: Beijing University Press, 2014. (in Chinese) -