海面舰船尾迹SAR图像仿真

粗糙海面上舰船Kelvin尾迹的电磁散射仿真韦尹煜;吴振森;路越【摘要】基于点源扰动方法建立了舰船的Kelvin尾迹模型,模拟了不同船速下的Kelvin尾迹.实际中,舰船尾迹除了受到舰船自身参数影响外,还会受到风驱海浪的调制.为了建立更为符合实际的Kelvin尾迹模型,借助PM(Pierson-Moscowitz)海谱生成粗糙海面,将其与Kelvin尾迹模型进行线性叠加后,得到了粗糙海面上舰船尾迹.在该复合模型基础上,利用物理光学法获得了不同入射情况下的单站、双站的雷达散射截面,实现了海面舰船尾迹电磁散射的模拟仿真.这一过程的实现,有利于更好地模拟舰船尾迹在合成孔径雷达(Synthetic Aperture Radar,SAR)中的成像仿真,同时为SAR成像中舰船尾迹以及舰船自身的识别奠定了基础.【期刊名称】《电波科学学报》【年(卷),期】2016(031)003【总页数】5页(P438-442)【关键词】Kelvin尾迹;PM谱;物理光学法;雷达散射截面(RCS)【作者】韦尹煜;吴振森;路越【作者单位】西安电子科技大学物理与光电工程学院,西安710071;西安电子科技大学物理与光电工程学院,西安710071;西安电子科技大学信息感知技术协同创新中心,西安710071;西安电子科技大学物理与光电工程学院,西安710071【正文语种】中文【中图分类】O451DOI 10.13443/j.cjors.2015071002对于水面舰船目标,传统探测方法都是以目标本身作为探测客体,而探测的信号源就是来自水中航行舰船的声反射性质、电磁性质等. 但是,随着各种先进消声技术、声对抗技术、消磁技术等隐身技术的发展和应用,使传统探测手段的工作能力被严重削弱.然而,在舰船航行过程中,其产生的尾迹具有尺度大、存在时间较长的特点,且无论采用什么隐身技术也很难被消除掉,完全满足被用作探测信号源的要求. 同时,随着合成孔径雷达(Synthetic Aperture Radar,SAR)的兴起,可以用星载或机载SAR图像中的船尾迹来识别运动舰船,并可通过尾迹形态得到舰船的吃水深度、航速等重要信息,因此研究舰船尾迹在海上目标探测中具有重要作用.近几十年来,国内外学者对舰船尾迹,尤其是对Kelvin尾迹进行了大量研究,Tuck等人[1]于1971年在Kelvin驻相法的基础上计算了船行波.2000年,Fontaine等[2]在考虑了舰船形状和大小的情况下,求解了在无粘流体中的船行波.2010年,孙荣庆等[3]对运动舰船的Kelvin尾迹和湍流尾迹的电磁散射建模作了相关研究.2013年,巩彪等[4]对Kelvin尾迹作了SAR图像仿真.2015年,Zilman等人[5]研究了粗糙海面上Kelvin尾迹在SAR成像中的识别问题．本文在不同环境条件下,对上述模型进行了模拟验证,建立了更为符合实际的Kelvin 尾迹模型;同时,利用PM海谱生成粗糙海面,并结合物理光学(Physical Optics,PO)法,实现了粗糙海面上舰船尾迹的电磁散射的模拟仿真.1.1 Kelvin尾迹建模舰船尾迹主要由分歧波和横断波组成,二者相干涉形成尖头波,由于尖头波的波长短,因而每个波前不能够被独立地分辨开来而表现为亮线,形成开尔文臂.开尔文臂的尾迹角约为±19.5°,如图1、图2所示.点源扰动模型是模拟Kelvin尾迹的基础,其原理是:将舰船目标看作一个以速度v移动的椭球体,对椭球体上每一个点按照航迹模式进行计算,再将每个点的航迹的近似值相叠加,从而得到整个Kelvin尾迹.舰船Kelvin尾迹的自由波模形如下:图3为不同船速下的Kelvin尾迹．Kelvin尾迹轮廓形如扇形，横断波沿舰船的航行方向传播，具有明显的周期特性．船速较小时，Kelvin波长较短，相应的尾迹高度起伏较小，随着船速的增加，波长与振幅明显增加，这与实际情况相符合．1.2 Kelvin尾迹与海面的线性叠加模型高频利用PM海谱[7-8]生成海面模型,风速u19.5=3.5 m/s,风向φw=0°,如图4(a)所示.同时,建立小型舰船尾迹模型,如图4(b)所示,相关参数为d=1.5 m,l=5m,b=2.5 m,v=2 m/s.仿真区域Lx=Ly=19.2 m,取样点数为M=N=512,面元尺寸为0.037 5 m×0.037 5 m.将海面模型与Kelvin尾迹模型进行线性叠加,得到海洋中的舰船尾迹模型,如图4(c)所示.之后,将利用此模型进行电磁散射计算.PO法是求解Helmholtz积分方程的近似方法,在求解电磁散射问题中应用广泛,其基本思想是:当入射波的波长λ和粗糙表面曲率半径ρ满足条件ρ>10λ时,粗糙面可视为由许多小平面拼接而成,电磁波在粗糙面的边缘和尖点绕射、面元间的多次散射等作用可以忽略.入射波不能直接照射的区域是暗区,入射波可以直接照射的区域是亮区,计算每一个亮区散射场,叠加后得到总散射场.粗糙面表面面元电磁流为: 利用物理光学法对1.2节中所建立的复合模型进行电磁散射计算,得到其单、双站雷达散射截面(Radar Cross-Section,RCS).3.1 单站电磁散射单站散射系数的仿真结果如图5所示.可以看到: 镜像区域中,来自介质海面的平均后向RCS明显大于带有舰船Kelvin尾迹的后向RCS; 随着入射角的增大,来自单纯海面的后向RCS明显小于叠加有Kelvin尾迹海面的RCS.这是由于尾迹的存在增大了海面起伏,改变了海面的粗糙度,使散射效果增强.3.2 双站电磁散射图6给出了入射角分别为0°和30°情况下的双站RCS.仿真结果表明：在镜像区域,单纯海面的RCS明显大于带有舰船Kelvin尾迹的RCS. 而在远离了镜像区域,单纯海面的RCS明显小于叠加有舰船Kelvin尾迹的RCS.原因是:尾迹的存在改变了海面粗糙度,使得在镜像区域相应的散射减弱,在其他区域的散射增强.本文利用点源扰动模型建立了舰船的Kelvin尾迹模型,模拟了不同船速下的舰船尾迹,发现当船速越大时,尾迹波长越长,波高也越大.同时,借助PM海谱生成粗糙海面,将其与Kelvin尾迹模型进行线性叠加得到了海面上的舰船尾迹.最后,针对舰船尾迹对粗糙海面电磁散射影这一问题进行了模拟仿真,利用PO获得了不同入射情况下的单、双站RCS,发现由于尾迹的存在增大了海面起伏,改变了海面的粗糙度,从而在镜向方向上的电磁散射减弱,非镜向方向上的散射效果增强.韦尹煜 (1990-),男,甘肃人,西安电子科技大学博士研究生,研究方向为目标与环境复合电磁散射及并行计算方法.吴振森 (1946-),男,湖北人,西安电子科技大学教授,博士生导师,近年来主要从事复杂介质、非均匀介质中的电磁波/光波的传播与散射、目标激光散射特性和电磁散射等方面的研究.【相关文献】[1] TUCK E. O. On line distribution of Kelvin sources[J]. Journal of ship research, 1971,12(2): 105-112.[2] FONTAINE E, FALTINSEN O M, COINTE R. New insight into the generation of ship bow waves[J]. Journal of fluid mechanics, 2000, 421(1): 15-38.[3] 孙荣庆, 罗根, 张民，等. 运动舰船的开尔文尾迹和湍流尾迹的电磁散射模型[C]//2010全国电磁散射与逆散射学术年会. 2010.[4] 巩彪, 黄韦艮, 陈鹏. Kelvin尾迹SAR图像仿真研究[J]. 海洋通报, 2013, 32(2): 208-213. GONG B, HUANG W G, CHEN P. Study on the SAR imaging simulation of kelvin ship wake[J]. Marine science bulletin, 2013, 32(2): 208-213. (in Chinese)[5] ZILMAN G, ZAPOLSKI A, MAROM M. On detectability of a ship's Kelvin wake in simulated SAR images of rough sea surface[J]. IEEE transactions on geoscience & remote sensing, 2015, 53(2):609-619.[6] OUMANSOUR K, WANG Y, SAILLARD J. Multifrequency SAR observation of a ship wake[J]. IEE proceedings-radar, sonar and navigation, 1996, 143(4): 275-280.[7] 李钦, 管荣生. 双尺度粗糙面的电磁波散射问题[J]. 电波科学学报, 1990, 5(3): 59-66.LI Q, GUAN R S. Wave scattering from rough surfaces with two scale[J]. Chinese journal of radio science, 1990, 5(3): 59-66. (in Chinese)[8] PIERSON W J, LIONEL M. A proposed spectral form for fully developed wind seas based on the similarity theory of S[J]. Journal of geophysical research, 1964, 69: 5181-5190.。

第50卷 第2期 海 洋 与 湖 沼Vol.50, No.2 2019年3月OCEANOLOGIA ET LIMNOLOGIA SINICAMar., 2019* 国家重点研发计划项目, 2016YFC1401002号; 国家自然科学基金项目, 41606201号, 41620104003号, 41576173号; 国家自然科学基金青年项目, 41706202号; 江苏省高校自然科学研究面上项目，17KJB170017号。

端和阳, 硕士研究生, E-mail: duanheyang@① 通信作者: 何宜军, 博士生导师, 教授, E-mail: yjhe@收稿日期: 2018-09-27, 收修改稿日期:2019-01-17一种改进的基于逆Omega-K 算法的海面场景SAR 原始数据仿真方法*端和阳 刘保昌 何宜军① 李秀仲(南京信息工程大学海洋科学学院 南京 210044)摘要 合成孔径雷达(SAR)海面场景原始数据仿真是研究海洋动力参数(表面波浪、风矢量和洋流)的有效工具。

目前海面场景原始数据仿真方法已经基于逆Omega-K 算法实现了海洋运动参数的空间变化。

但是目前仅仅讨论了正侧视情况下的海面场景仿真, 应用范围有限, 同时没有考虑Stokes 漂流以及Bragg 相速度的影响, 而这两者都是存在于真实海面的。

通常情况下为了反演得到海面流场的二维速度矢量, 雷达需要从两个不同的方位方向观察海面的同一个区域, 因此这就需要考虑大斜视的雷达波束, 同时Stokes 漂流和Bragg 相速度是SAR 海表面流场观测不容忽视的两种运动。

本文在不改变原有正侧视逆Omega-K 算法的情况下, 通过增加重新计算零方位时刻的斜视波束中心位置坐标, 并据此确定SAR 原始数据在多普勒域的位置来将其扩展到大斜侧视逆Omega-K 算法, 并通过时域Stokes 漂流公式到频域内离散化Stokes 漂流公式的推导来加入Stokes 漂流, 以及根据Bragg 散射机制加入了Bragg 相速度。

59中国设备工程Engineer ing hina C P l ant中国设备工程 2017.04（下）1 SAR 图像舰船尾迹检测的意义为了保卫海域及岛屿的主权，监测是否有其他国家的舰船侵入我国海域进行军事活动是十分重要以及必要的。

雷达（R a d i o D e t e c t i n g a n d R a n g i n g , RADAR）是一种基本的无线电探测装置。

它可实现对目标的探测和定位，在军事领域起着十分重要的作用。

随着科技的发展与进步，合成孔径雷达（Synthetic Aperture Radar, SAR）技术也日益成熟，其可获取更多的信息。

SAR 是利用一个小天线沿着长线阵的轨迹等速移动并辐射相参信号, 把在不同位置接收的回波进行相干处理, 从而获得较高分辨率的成像雷达。

其具有全天时、全天候的观测能力。

SAR 由于其全天时全天候的特性可以广泛应用在海洋遥感中。

研究人员现利用海洋 SAR 图像进行舰船检测。

为了在短时间内有效地检测舰船及舰船尾迹，关于自动检测算法的研究十分重要且有意义。

一般来说，检测舰船及舰船尾迹都十分有必要；但有时研究人员们对于检测舰船尾迹更感兴趣，因为通过检测尾迹从而来发现舰船目标有如下优点。

（1）在SAR 图像中，舰船尾迹结构通常比舰船本身大许多，而且更为明显，可以更容易进行检测和识别；（2）舰船尾迹可以用于估计运动舰船的实际位置、航速和航向。

比如说湍流尾迹的方向可以帮助人们确定舰船运动的方向；开尔文臂的波长可以用于估计舰船的速度等；（3）舰船尾迹可以用于检测神秘舰船。

所谓的神秘舰船，比如说潜艇等。

对于神秘舰船，由于它在水下运动，在SAR 图像上并不可见。

但在某些情况下，潜艇等神秘舰船可以在海面上产生尾迹图像。

因此，尾迹检测技术可以用于神秘舰船的研究，对于国防有重大意义。

早在20世纪80年代科学家们已经发现SAR 图像上运动舰船存在四种尾迹表现形式：开尔文（Kelvin）尾迹、湍流尾迹、窄V 型尾迹与内波尾迹，其尾迹结构如图1所示。