激光引信近程探测海面回波特性研究

白明健,杨 丽,依 研

(沈阳理工大学 装备工程学院, 沈阳 110159)

0 引言

脉冲激光引信因其方向性好、抗光电干扰能力强、探测精度高和良好的距离截止特性等优点,被广泛应用在导弹领域[1-2]。海战场环境下,搭载激光引信的导弹在打击海面目标时,激光引信发射的激光束会有一部分照射在海面上发生散射效应,导致引信产生虚警现象,使引信发生早炸,严重影响导弹的毁伤效能[3]。目前,国内外众多学者在研究激光束的海面散射时,主要研究蓝绿波段激光的海面散射性能,对近红外波段激光研究很少。因此研究近红外波段激光的海面散射特性有非常重要的意义。

路明等[4]提出一种激光引信双波束探测抗海杂波技术,对发射和接收光学系统着重进行分析,最后分析了该引信系统的抗海杂波性能。Zhang等[5]基于线性叠加法利用JONSWAP谱的谱特点完成了海浪的模拟。张莹珞等[6]利用Mie理论和蒙特卡洛法分析了海水中激光的传输特性。Sahoo等[7]基于蒙特卡洛法研究了气-海界面、气泡和粒子散射对水下光场分布的影响。Li等[8]为提高激光引信的探测能力,研究了激光引信的回波信号检测方法和回波功率计算方法。史铭森[9]在研究海面背景下目标的激光束散射时,利用海面模型完成了不同海面背景下散射回波的研究。韩复兴等[10]给出海浪谱信息的动态边界条件,并在此基础上构建起伏海面速度模型。针对上述关键问题,本研究中以脉冲激光引信为背景,通过建模和仿真研究激光引信的海面回波特性。该方法可为后续研究激光引信抗海杂波干扰提供有力支撑。

1 基于海浪谱的海面建模

1.1 线性叠加法海面建模

JONSWAP谱海浪相对于海面的瞬时波高为[11]:

(1)

式(1)中:Aij为振幅;ωi=2πfi为角频率;fi为频率;ki=2π/λi为波数;λi为波长;θj为波在xy水平面的传播方向,代表方向角;εij为初相角;t表示时间。

由国家标准规定可知,风力等级和波高以及风速的关系[12]。具体可见表1所示。

表1 蒲福风力等级

选取风力等级分别为2级、5级和8级,规定风区长度10 km,方向角30°、45°和60°进行海面建模。45°方向角下典型风速对应的JONSWAP谱海面模型如图1所示。相同风速下典型方向角的JONSWAP谱海面模拟如图2所示。

图1 典型风速的JONSWAP谱海面模拟

图2 典型方向角的JONSWAP谱海面模拟

根据图1所示典型风力等级的JONSWAP谱海面模型可以得到表2。

表2 典型风力等级的JONSWAP谱海面波高

根据图2所示典型方向角的JONSWAP谱海面模型可以得到表3。

由于文氏谱和其他海谱的区别,描述海面波高时有所不同。文氏谱海面波高可表示为[13]:

(2)

A=cos[ki·cos(θwind+φj)·x+

ki·sin(θwind+φj)·y+ωi·t+εij]

(3)

表3 典型方向角的JONSWAP谱海面波高

选风力等级为2级、5级和8级,规定风区长度10 km,方向角30°、45°和60°进行海面建模。45°方向角下典型风力等级的文氏谱海面模拟如图3所示。规定风力等级为5级,典型方向角的文氏谱海面模拟如图4所示。根据图3所示典型风速的文氏谱海面模拟可得表4。根据图4所示典型方向角的文氏谱海面模型可得表5。

图3 典型风速的文氏谱海面模拟

图4 典型方向角的文氏谱海面模拟

表4 典型风速的文氏谱海面波高

表5 典型方向角的文氏谱海面波高

对比基于线性叠加法创建的2种海浪谱模型,JONSWAP谱海面模型海面波高误差为10%,文氏谱海面模型海面波高误差为6%,文氏谱海面模型更加注重海浪的低频部分,当方向角发生变化时,海浪受方向角的影响发生变化的程度小于JONSWAP谱建模时方向角改变所引起的变化。对照蒲福风力等级表,基于文氏谱的海浪建模,海面误差更小,海浪更加稳定。真实性更好,实时性更好,海面细节更好。具体可见表6所示。

表6 2种海浪谱建模对照表

1.2 线性滤波法海面建模

线性滤波法就是选择合适的白噪声,经过傅里叶变换对白噪声进行处理,将处理后的白噪声与选定的海浪谱进行卷积,然后再对其进行快速傅里叶逆变换,最后就可通过滤波器得到所需的海面模型。

海面波高表达式为[14]:

(4)

式(4)中:ξr、ξi是随机数生成的余弦和正弦值;S(K)为随机粗糙面功率谱密度,kx=2πm/Lx,ky=2πm/Ly。Lx和Ly为在海面坐标系下,X和Y方向的长度;m和n为离散间隔数。

根据线性滤波法海面建模方式,选取风力等级为5级,规定风区长度10 km,方向角为45°完成海面建模,网格数分别为128×128、256×256和512×512。根据以上设定的参量完成海面建模。不同网格划分的JONSWAP谱海面模型如图5所示。不同网格划分的文氏谱海面模型如图6所示。

图5 不同网格数的JONSWAP谱海面模拟

由图5不同网格数的JONSWAP谱海面模型可得表7。由图6不同网格数的文氏谱海面模型绘制表8。

表7 不同网格数的JONSWAP谱海面波高

表8 不同网格数的文氏谱海面波高

对比2种基于线性滤波法的海面建模方式,采用JONSWAP谱进行海面建模时,海面波高更加尖锐,相比文氏谱建模,计算耗时长,实时性差。

2 激光束海面散射模拟

2.1 JONSWAP谱激光束模拟

选取基于线性滤波法的JONSWAP谱海面模型。用几何光学的方法加入光束进行反射光的模拟,图7为JONSWAP谱模型下激光入射角分别为15°、45°和75°时入射光与反射光的状态。图8是45°入射角下风速为3、10、20 m/s的反射情况。其中,蓝色光线为入射光线,黑色光线为反射光线。

由图7可知,当风力等级为5级时,随着入射角的增大,海面反射趋近于镜面反射,反射光线与入射光线以海面激光入射点呈圆锥形。当入射角增加到75°时,入射激光趋近于掠入射,反射光线偏离镜面反射方向。当激光束趋于掠入射状态时,由于海浪之间会出现相互遮挡,反射激光的反射率会减少。由图8所示,当入射角一定时,随着风速的增大,粗糙海面对散射光的分布范围逐渐变宽,并且以海面激光入射点所呈现的圆锥形会随着风速的增大而增大。

2.2 文氏谱激光束模拟

选取基于线性滤波法的文氏谱模型。进行激光束的模拟。图9为文氏谱模型下激光入射角分别为15°、45°和75°时入射光与反射光的状态。图10是文氏谱模型下45°入射角风速为3、10、20 m/s的反射情况。

由图9可知,对比基于JONSWAP谱的相同风力等级不同入射角的激光束模拟,基于文氏谱的海面激光反射模拟,当风力等级为5级时,随着入射角的增大,海面激光束反射趋近于镜面反射,当入射角增加到75°时,入射光线趋近于掠入射,反射光线未偏离镜面反射方向。

由图10所示,当入射角一定时,随着风速的增大,粗糙海面对散射光的分布范围未发生变化,以海面激光入射点所呈现的圆锥形会随着风速的增大而增大。

图7 不同入射角的JONSWAP谱激光反射模拟

图8 不同风速JONSWAP谱激光反射模拟

图9 不同入射角下文氏谱海面激光反射模拟

图10 不同风速下文氏谱激光反射模拟

3 激光引信海面散射特性

3.1 回波参数设置

激光在大气中传输要受到大气的折射、散射、吸收、热晕等影响。这些影响将导致激光束发生衰减[15]。朗伯-比尔定律是光吸收的基本定律,描述了物质对光的吸收程度和吸光介质厚度、吸光物质浓度之间的关系,其表达式为[16]:

(5)

式(5)中:A为吸光度;I0(λ)为入射光强度;It(λ)为透射光强度;T为透过率;ε(λ)为吸光系数;b为介质的厚度;c为吸光介质的浓度。根据朗伯-比尔定律可以得到透过率的值。

选取南海某域海面,当激光波长为905 nm时,基于Cox-Munk模型模拟的菲涅尔系数范围介于0.021～0.043。

根据本文所需,选取适当的海面回波参数,设定引信探测海面的距离为15 m,对其进行200 ms的采集,具体的参数如表9所示。

表9 海面回波参数

3.2 激光束海面回波特性

根据激光回波方程,当激光引信距离海面为Rn时,海面划分的微小面元Sn被激光引信发射的激光束所照射,其反射的回波被光学接收系统所接收,此时,回波能量可用PRn表示为[17]:

(6)

式(6)中:PT为发射功率;ρ为海面面元Sn的菲涅尔反射率;τR为接收光学系统透过率;τT为发射光学系统透过率;τA为2倍海面距离Rn的大气透过率;σn为划分的海面微小面元面积;φn为入射光与海面面元法线的夹角;An为接收光学系统的有效口径;Rn为距离海面的距离。St(Rn)为距离海面Rn处发射光束截面积。

选取基于线性滤波法的JONSWAP谱海面模型,进行激光引信海面探测仿真。图11为激光引信探测海面的回波仿真场景图。

图11 激光引信探测海面示意图

由表9可知,回波仿真的弧矢方向视场角为60°,子午方向视场角为1°,选取入射光与面元法线的夹角分别为0°、15°和30°进行仿真。分别选取风速3、10、20 m/s 3种不同风速的海况进行仿真。设定激光引信与海面的交会速度为 200 m/s,激光引信距离海面15 m。采集200 ms的回波能量。图12为不同风速下激光引信的回波能量随入射角的变化示意图。

由图12可知,当海面风速一定时,激光束入射角的变化,会对海面的反射回波功率产生很大的影响。当入射角为0°时,海面的回波功率最大,随着激光束入射角的逐渐增大,海面会产生镜面反射,海面面元能够被激光束照射的部分越来越小,海面反射回波随入射角的增加连续降低;当激光束入射角一定时,随着风速的增大,海面波动幅度增加,发生大斜率海面的概率增大,能够进入引信接收视场的海面面元面积就会增大,因此接收到的海面反射回波的功率就会增加。

图12 不同风速下激光引信回波示意图

4 结论

1) 对比2种基于线性叠加法生成的海面模型,基于线性叠加法的文氏谱海面模型海面波高误差约为6%,海面真实性、实时性更好,海面呈现的细节更全面。基于线性叠加法的JONSWAP谱海面模型海面波高误差约为10%,建模时更加注重海面的高频部分;2种基于线性滤波法生成的海面模型中,文氏谱生成的海面模型,具有计算耗时短,实时性好的特点。

2) 对比基于线性滤波法的2种不同海浪谱激光束模拟,当风速一定时,入射角由15°增加到75°,海浪之间会出现相互遮挡,反射激光的反射率会减少;当入射角一定时,随着风速的增加,JONSWAP谱海面模型对散射光的分布范围逐渐变宽,而文氏谱海面模型对散射光的分布范围未发生变化。

3) 激光束入射角和风速是影响海面反射回波的重要因素,入射角和风速的改变会影响激光引信光学接收系统对海面面元反射激光的接收,从而影响海面反射回波能量的大小。