留言板

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

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

舰船尾流激光探测跟踪方法与试验

宗思光 张鑫 曹静 梁善永 李斌

宗思光, 张鑫, 曹静, 梁善永, 李斌. 舰船尾流激光探测跟踪方法与试验[J]. 红外与激光工程, 2023, 52(3): 20220507. doi: 10.3788/IRLA20220507
引用本文: 宗思光, 张鑫, 曹静, 梁善永, 李斌. 舰船尾流激光探测跟踪方法与试验[J]. 红外与激光工程, 2023, 52(3): 20220507. doi: 10.3788/IRLA20220507
Zong Siguang, Zhang Xin, Cao Jing, Liang Shanyong, Li Bin. Method and experiment of laser detection and tracking of ship wake[J]. Infrared and Laser Engineering, 2023, 52(3): 20220507. doi: 10.3788/IRLA20220507
Citation: Zong Siguang, Zhang Xin, Cao Jing, Liang Shanyong, Li Bin. Method and experiment of laser detection and tracking of ship wake[J]. Infrared and Laser Engineering, 2023, 52(3): 20220507. doi: 10.3788/IRLA20220507

舰船尾流激光探测跟踪方法与试验

doi: 10.3788/IRLA20220507
基金项目: 国防科研基金(2019-JCJQ-JJ-056, 2020-JCJQ-20-99-00-02)
详细信息
    作者简介:

    宗思光,男,副教授,博士,主要从事水下激光探测、激光通信、激光致声等方面的研究

  • 中图分类号: TN249

Method and experiment of laser detection and tracking of ship wake

Funds: The National Defense Foundation of China (2019-JCJQ-JJ-056, 2020-JCJQ-20-99-00-02)
  • 摘要: 舰船尾流激光探测跟踪是水下航行器对舰船进行探测、识别、跟踪的新手段。论文基于舰船尾流分布特性、气泡目标特性,采用蒙特卡洛仿真方法,实现了多尺度、宽数密度、大厚度舰船尾流气泡群的后向散射回波信号特性仿真,得到了水下航行器载激光探测系统在搜索、跟踪阶段信号的变化趋势,以及不同目标舰船的激光后向回波信号变化强度,可有效模拟激光探测系统对舰船尾流目标特性的真实跟踪状态。对于大型船只,当激光探测系统位于尾流之下时,航行器距舰船目标越近,尾流气泡激光回波越强,脉冲宽度展宽幅度越大;当激光探测系统位于尾流之中时,航行器距舰船目标越近,尾流气泡激光回波越弱,脉冲宽度变窄幅度越大。探测系统位于尾流之下时与探测系统位于尾流之中时,信号变化相反。小型船只信号变化趋势基本与大型船只保持一致,但尾流激光探测回波强度变低。开展了湖泊环境下船舶尾流激光探测跟踪试验,当探测系统在尾流之下时,大型船只尾流激光回波信号信噪比高,小型船只尾流激光难以检测。探测系统位于尾流之中时,大小船只尾流激光探测系统都可实现有效探测。论文可为舰船尾流探测实际工程应用提供支撑。
  • 图  1  舰船尾流的分布特性。(a) 俯视图;(b) 横向剖面图;(c) 侧向剖面图

    Figure  1.  Distribution characteristics of ship wake. (a) Top view; (b) Transverse section; (c) Lateral profile

    图  2  光子模拟运动程序图

    Figure  2.  Program diagram of photon simulated motion

    图  3  探测可能存在的形式。 (a) 探测位于尾流之中;(b) 探测位于尾流之下

    Figure  3.  Possible forms of detection. (a) Detection in wake; (b) Detection under wake

    图  4  信号变化趋势。(a) 尾流下; (b) 尾流中

    Figure  4.  Signal change trend. (a) Under the wake; (b) In wake

    图  5  信号变化趋势。 (a)尾流下; (b) 尾流中

    Figure  5.  Signal change trend. (a) Under the wake; (b) In wake

    图  6  信号变化趋势。(a)尾流下; (b) 尾流中

    Figure  6.  Signal change trend. (a) Under the wake; (b) In wake

    图  7  信号变化趋势。 (a)尾流下; (b) 尾流中

    Figure  7.  Signal change trend. (a) Under the wake; (b) In wake

    图  8  试验系统及试验过程。(a)探测系统;(b)试验船;(c)探测系统;(d)湖泊测试跟踪过程;(e)跟踪过程示意图

    Figure  8.  Experimental system and process. (a) Detection system; (b) Experimental ship; (c) Detection system; (d) Lake test tracking pro-cess; (e) Schematic diagram of tracking process

    图  9  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

    Figure  9.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    图  10  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号;(b)平滑处理后的信号

    Figure  10.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    图  11  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

    Figure  11.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    图  12  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

    Figure  12.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    表  1  典型舰船尾流的分布[11]

    Table  1.   Wake distribution of common ships[11]

    Destroyer
    model
    Speed/knThickness
    at 360 m/m
    Width/mDraft
    Rathborn10-126.5±1.1--
    Hopewell107.7±1.2125.8
    Evian134.4±1.17.52.9
    下载: 导出CSV

    表  2  实测舰船尾流中气泡分布密度[16]

    Table  2.   Measured bubble distribution density in ship wake[16]

    Bubble
    radius/μm
    1 min/m−33 min/m−35 min/m−3
    10704.9×1024.1×101-
    4001.79×1045.6×1037.8×102
    1604.6×1051.98×1059.3×104
    805.5×1062.61×1061.45×106
    10-10702.83×1081.32×1087.16×107
    下载: 导出CSV

    表  3  大型船只仿真条件设置

    Table  3.   Simulation conditions for large vessel

    Simulation
    conditions
    Zone 1Zone 2Zone 3Zone 4Zone 5Zone 6Zone 7Zone 8
    Movement
    mode
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    Bubble size/μm20-8080-140140-200200-260260-320320-380380-440440-500
    Bubble thickness/m2-33-44-55-66-77-88-99
    Bubble density/m−3105-106 106-3×1063×106-6×1066×106-107 107-108108-2×1082×108-3×1083×108
    下载: 导出CSV

    表  4  小型船只仿真条件设置

    Table  4.   Small vessel simulation condition settings

    Simulation
    conditions
    Zone 1Zone 2Zone 3Zone 4
    Movement
    mode
    In/under
    wake
    In/under
    wake
    In/under
    wake
    In/under
    wake
    Bubble size/μm20-6060-100100-300300-500
    Bubble
    thickness/m
    0-11-22-33-4
    Bubble
    density/m−3
    104-105 105-106 106-107 107
    下载: 导出CSV
  • [1] Wang Huili, Qi Yi, Liu Huanying. Boundary detection method of infrared image of ship wake [J]. Infrared and Laser Engineering, 2013, 42(2): 524-527. (in Chinese)
    [2] Wang Ping, Du Yongcheng, Yang Li, et al. Numerical and experimental study on the buoyancy and diffusion law of submarine thermal wake based on overlapping grid technology and VOF model [J]. Infrared and Laser Engineering, 2019, 48(4): 0404002. (in Chinese) doi:  10.3788/irla201948.0404002
    [3] Jiang Xin, Chen Wuxiong, Nie Haitao, et al. Real-time ship target detection based on aerial remote sensing images [J]. Optics and Precision Engineering, 2020, 28(10): 2360-2369. (in Chinese) doi:  10.37188/OPE.20202810.2360
    [4] Tang Meng, Zhang Yu. Study on polarization characteristics of wake microbubbles detected by laser [J]. Infrared and Laser Engineering, 2020, 49(1): 0105006. (in Chinese) doi:  10.3788/IRLA202049.0105006
    [5] Li Yiyue. Turbulence spectrum separation of wind lidar using independent component analysis [J]. Optics and Precision Engineering, 2020, 28(5): 1029-1037. (in Chinese)
    [6] Zhang Hui, Gao Xiaocheng. Analysis of optical characteristics of ship wake and research on wake detection technology [J]. Ship Science and Technology, 2021, 43(24): 52-54. (in Chinese)
    [7] Gao Kexin, Jin Liangan, Yuan Zhijiang, et al. Study on bubble number density attenuation model of ship bubble wake [J]. China Test, 2019, 45(8): 61-66. (in Chinese)
    [8] Jin Liangan, Yan Xuefei, Wang Yong, et al. Study on bubble coalescence in ship wake based on large eddy simulation method and modified bubble equilibrium equation [J]. Science, Technology and Engineering, 2014, 14(24): 141-145. (in Chinese)
    [9] Chen P, Li X N, Zhang G. Rapid detection to long ship wake in synthetic aperture radar satellite imagery [J]. Journal of Oceanology and Limnology, 2019, 37(5): 1523-1532. (in Chinese) doi:  10.1007/s00343-019-8221-y
    [10] Li Hai, Jia Hongguang, Chen Zaibin. Analysis and experiment on aerodynamic characteristics of coaxial rotor system [J]. Optics and Precision Engineering, 2021, 29(9): 2140-2148. (in Chinese) doi:  10.37188/OPE.20212909.2140
    [11] Zhang Qun, Wang Yingmin. Multi bubble model and finite element analysis in wake [J]. Torpedo Technology, 2014, 22(4): 316-320. (in Chinese)
    [12] Zhang Yinbo, Li Sining, Jiang Peng, et al. PCA feature extraction and elastic BP neural network for underwater bubble recognition [J]. Infrared and Laser Engineering, 2021, 50(6): 20200352. (in Chinese) doi:  10.3788/IRLA20200352
    [13] Gao Jiang, Zhang Jingyuan Yang Li. Research status of ship bubble wake characteristics [J]. Ship Science and Technology, 2008(4): 27-32. (in Chinese)
    [14] Tian Hengdou, Jin Liangan, Wang Yong, et al. Study on bubble distribution in ship wake considering single bubble motion characteristics [J]. Acta Armamentarii, 2011, 32(9): 1126-1130. (in Chinese)
    [15] Guo C Y, Wu T C, Zhang Q, et al. Numerical simulation and experimental research on wake field of ships under off-design conditions [J]. China Ocean Engineering, 2016, 30(5): 821-834. (in Chinese) doi:  10.1007/s13344-016-0053-3
    [16] Zhang X X, Wu Z S, Su X. Influence of breaking waves and wake bubbles on surface-ship wake scattering at low grazing angles [J]. Chinese Physics Letters, 2018, 35(7): 074101. (in Chinese) doi:  10.1088/0256-307X/35/7/074101
    [17] Wang B L, Guo X Y, He C. Numerical simulations of wake signatures around high-speed ships [J]. Journal of Hydrodynamics, 2014, 26(6): 986-989. (in Chinese) doi:  10.1016/S1001-6058(14)60109-8
    [18] Li Xiaobo, Li Wei, Sun Tianyu, et al. Analysis and design of onboard interface of Tianwen-1 high-resolution camera [J]. Optics and Precision Engineering, 2022, 30(2): 227-236. (in Chinese)
    [19] Feng Z F, Wang P P, Yang X, et al. Evaluation of light penetration of LED phototherapy apparatus in skin [J]. Optics and Precision Engineering, 2022, 30(10): 1139-1150. (in Chinese) doi:  10.37188/OPE.20223010.1139
    [20] Cheng Yan, Yu Xuelian, Qian Weixian, et al. Ship wake extraction and detection from infrared remote sensing images [J]. Infrared and Laser Engineering, 2022, 51(2): 20210844. (in Chinese) doi:  10.3788/IRLA20210844
  • [1] 蒋衍, 朱精果, 刘汝卿, 李锋, 姜成昊, 孟柘, 王宇.  脉冲激光多回波峰值检测电路设计 . 红外与激光工程, 2023, 52(5): 20220625-1-20220625-8. doi: 10.3788/IRLA20220625
    [2] 张鑫, 宗思光, 李斌, 余扬.  浑浊水体激光后向散射特性仿真与实验研究 . 红外与激光工程, 2023, 52(1): 20220280-1-20220280-7. doi: 10.3788/IRLA20220280
    [3] 苗澍茁, 安宁, 高健, 温冠宇, 宋清丽, 董雪, 马磊, 范存波.  SLR系统地靶数值仿真及数据处理 . 红外与激光工程, 2021, 50(9): 20200402-1-20200402-9. doi: 10.3788/IRLA20200402
    [4] 崔晓宇, 陶雨婷, 刘群, 徐沛拓, 刘志鹏, 王晓彬, 陈扬, 周雨迪, 刘东.  采用半解析蒙特卡洛技术模拟星载海洋激光雷达回波信号的软件 . 红外与激光工程, 2020, 49(2): 0203009-0203009. doi: 10.3788/IRLA202049.0203009
    [5] 陈鹏, 赵继广, 宋一铄, 王燊.  气溶胶环境下FMCW与脉冲激光探测性能对比 . 红外与激光工程, 2020, 49(6): 20190399-1-20190399-10. doi: 10.3788/IRLA20190399
    [6] 战俊彤, 张肃, 付强, 段锦, 李英超, 姜会林.  不同湿度环境下可见光波段激光偏振特性研究 . 红外与激光工程, 2020, 49(9): 20200057-1-20200057-7. doi: 10.3788/IRLA20200057
    [7] 王云鹏, 胡以华, 雷武虎, 郭力仁.  典型旋翼形状参数微多普勒激光探测计算方法 . 红外与激光工程, 2018, 47(9): 906003-0906003(9). doi: 10.3788/IRLA201847.0906003
    [8] 马俊凯, 罗海波, 常铮, 惠斌, 周晓丹, 侯德飞.  基于可变形模型的目标跟踪算法 . 红外与激光工程, 2017, 46(9): 928001-0928001(9). doi: 10.3788/IRLA201746.0928001
    [9] 黄爱萍, 张莹珞, 陶林伟.  蒙特卡洛仿真的水下激光通信信道特性 . 红外与激光工程, 2017, 46(4): 422004-0422004(6). doi: 10.3788/IRLA201746.0422004
    [10] 邢运龙, 李艾华, 崔智高, 方浩.  改进核相关滤波的运动目标跟踪算法 . 红外与激光工程, 2016, 45(S1): 214-221. doi: 10.3788/IRLA201645.S126004
    [11] 邵春艳, 丁庆海, 罗海波, 李玉莲.  采用高维数据聚类的目标跟踪 . 红外与激光工程, 2016, 45(4): 428002-0428002(10). doi: 10.3788/IRLA201645.0428002
    [12] 崔雄文, 吴钦章, 蒋平, 周进.  子空间模型下的仿射不变目标跟踪 . 红外与激光工程, 2015, 44(2): 769-774.
    [13] 高昕, 李希宇, 冯灵洁, 唐嘉.  强度相干成像中符合计数方法的优化与仿真 . 红外与激光工程, 2015, 44(11): 3454-3462.
    [14] 郑勇辉, 孙华燕, 赵延仲, 张令军.  基于Zemax的探测激光气动光学畸变快速仿真 . 红外与激光工程, 2015, 44(S1): 80-85.
    [15] 黎志华, 李新国.  基于OpenCV的红外弱小运动目标检测与跟踪 . 红外与激光工程, 2013, 42(9): 2561-2565.
    [16] 彭晨, 陈钱, 钱惟贤, 徐富元.  复杂地面场景下的红外运动目标跟踪 . 红外与激光工程, 2013, 42(6): 1410-1414.
    [17] 查冰婷, 张合, 张祥金, 许敬.  水下单光束脉冲激光方位识别系统角度参数优化设计 . 红外与激光工程, 2013, 42(4): 895-899.
    [18] 秦刚, 杨郁, 张建生.  船舰远程尾流散射光偏振特性的蒙特卡洛模拟 . 红外与激光工程, 2013, 42(7): 1730-1736.
    [19] 万俊, 张晓晖, 饶炯辉, 梁瑞涛.  基于独立成分分析的舰船气泡尾流后向散射光信号处理 . 红外与激光工程, 2013, 42(1): 244-250.
    [20] 王慧丽, 齐异, 刘焕英.  舰船尾流红外图像边界检测方法 . 红外与激光工程, 2013, 42(2): 524-527.
  • 加载中
图(12) / 表(4)
计量
  • 文章访问数:  95
  • HTML全文浏览量:  12
  • PDF下载量:  47
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-07-21
  • 修回日期:  2022-08-26
  • 网络出版日期:  2023-03-20
  • 刊出日期:  2023-03-25

舰船尾流激光探测跟踪方法与试验

doi: 10.3788/IRLA20220507
    作者简介:

    宗思光,男,副教授,博士,主要从事水下激光探测、激光通信、激光致声等方面的研究

基金项目:  国防科研基金(2019-JCJQ-JJ-056, 2020-JCJQ-20-99-00-02)
  • 中图分类号: TN249

摘要: 舰船尾流激光探测跟踪是水下航行器对舰船进行探测、识别、跟踪的新手段。论文基于舰船尾流分布特性、气泡目标特性,采用蒙特卡洛仿真方法,实现了多尺度、宽数密度、大厚度舰船尾流气泡群的后向散射回波信号特性仿真,得到了水下航行器载激光探测系统在搜索、跟踪阶段信号的变化趋势,以及不同目标舰船的激光后向回波信号变化强度,可有效模拟激光探测系统对舰船尾流目标特性的真实跟踪状态。对于大型船只,当激光探测系统位于尾流之下时,航行器距舰船目标越近,尾流气泡激光回波越强,脉冲宽度展宽幅度越大;当激光探测系统位于尾流之中时,航行器距舰船目标越近,尾流气泡激光回波越弱,脉冲宽度变窄幅度越大。探测系统位于尾流之下时与探测系统位于尾流之中时,信号变化相反。小型船只信号变化趋势基本与大型船只保持一致,但尾流激光探测回波强度变低。开展了湖泊环境下船舶尾流激光探测跟踪试验,当探测系统在尾流之下时,大型船只尾流激光回波信号信噪比高,小型船只尾流激光难以检测。探测系统位于尾流之中时,大小船只尾流激光探测系统都可实现有效探测。论文可为舰船尾流探测实际工程应用提供支撑。

English Abstract

    • 舰船在航行的过程中由于螺旋桨的空化作用、海面波浪的破碎以及吃水线部分卷入大量空气,在舰船尾部形成了含有大量气泡的气幕带,即舰船尾流,其拥有与四周水体环境差别极大的光学特性[1]。通过对舰船尾流激光特性的研究,可以进一步判断舰船在海洋中的航行轨迹、速度等特征,进而实现探测系统等水下航行器的精确制导和毁伤打击[2]。由于激光尾流探测具有比较强的军事应用背景,国外很少有相关的报道。20世纪90年代,加拿大海洋科学研究院采用激光后向散射技术对尾流的持续时间及微气泡幕的浓度等特征信息进行测量;Stramski最早进行了试验研究,他利用全息照相方法得到了气泡群尺度的分布,并发现尾流中气泡的激光后向散射强度的近十分之一由气泡群产生[3]。国内对舰船尾流方面的研究起步也较迟,饶炯辉、夏珉等人利用蒙特卡洛模拟方法研究了气泡群的后向光学散射特性,检测到11.5 s以内实船尾流气泡的散射信号;石晟伟、马治国等人基于蒙特卡洛模拟仿真对气泡幕位置、气泡尺寸、气泡数密度进行了模拟仿真[4];王向伟、周田华等人建立了尾流气泡仿真模型,对不同海域以及尾流远区到尾流中心区域气泡的后向散射特性进行了模拟[5]。舰船尾流是一个动态变化的环境,在不同舰船和不同环境下的舰船尾流分布特性及气泡目标特性也不尽相同,要想实现对舰船的精确打击,必须基于舰船尾流分布特性及气泡目标特性进行研究,通过模拟不同舰船尾流情况下的回波信号变化趋势,为舰船尾流在湖泊、海洋等外场跟踪探测中提供理论与仿真支撑。

      文中以舰船尾流激光探测工程应用为背景,分析了舰船尾流分布特性及气泡目标特性变化规律。基于舰船尾流分布特性及气泡目标特性建立仿真环境,采用蒙特卡洛方法,对多尺度、宽数密度、大厚度舰船尾流气泡群进行模拟仿真,通过分析不同情况下舰船尾流后向散射回波信号,可有效模拟探测系统在舰船尾流目标特性下的真实状态,研究了探测系统在搜索、跟踪阶段的信号变化趋势,以及不同目标舰船的回波信号变化强度,当激光探测系统位于尾流之下时,大型船只的气泡回波幅度缓慢上升,气泡脉冲宽度显著展宽,且距舰船越近,气泡回波变化越明显。激光探测系统位于尾流之中时,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。激光探测系统位于尾流之下时与探测系统位于尾流之中时,信号变化相反。开展了湖泊环境下舰船尾流激光跟踪探测试验,验证了探测系统在尾流之下时,气泡回波信号变化为气泡幅度缓慢上升,气泡脉冲宽度显著展宽;探测系统位于尾流之中时,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。这可为舰船尾流探测在实际工程应用中提供支撑。

    • 舰船尾流的分布特性如图1所示,舰船尾流分布特性一般指舰船尾流场厚度[6]、舰船尾流有效检测长度[7]、舰船尾流的宽度[8]等特性。

      图  1  舰船尾流的分布特性。(a) 俯视图;(b) 横向剖面图;(c) 侧向剖面图

      Figure 1.  Distribution characteristics of ship wake. (a) Top view; (b) Transverse section; (c) Lateral profile

      舰船尾流的长度一般为20~50 LL指船体长度),其持续时间可达10~15 min,舰船尾流的有效检测长度与航速及可检测能力有关[9]。其公式可表示为:

      $$ {L}_{T}={C}_{T}\cdot v $$ (1)

      式中:${L}_{T} $为可检测舰船尾流的长度;$ v $为舰船航速;$ {C}_{T} $为常数,代表尾流检测能力。

      已有的研究及实船测试数据表明舰船尾流的宽度和舰船的船型有关。其具体示意图如图1(a)所示,舰船尾流的起始宽度约为舰船宽度的一半,随后以40°~60°扩散角向两边扩散,到达一定距离后,舰船尾流的宽度大约为舰船的2.5倍。此后,舰船尾流以不大于1°的扩散角向两边扩散,直至舰船尾流消失[10]。其公式可以表示为:

      $$ W=\left\{\begin{array}{l}0.5B+2vt\mathrm{tan}\alpha \qquad\qquad\qquad\;\; W\leqslant 2.5B\\ 2.5B+2\mathrm{tan}\beta vt-2B\mathrm{tan}\beta /\mathrm{t}\mathrm{a}\mathrm{n}\alpha \;\;W > 2.5B\end{array}\right. $$ (2)

      式中:W为舰船尾流的宽度;B为舰船的宽度;t为航行时间;$ \alpha $$ \,\beta $为舰船尾流的扩散角。

      舰船尾流的初始厚度与舰船的吃水深度有关。其具体示意图如图1(b)、(c)所示,大型船只为吃水深度的2倍,小型船只为吃水深度的4倍左右,而后又与舰船的航行速度有关,通常为航速越高,舰船尾流越厚,高速航行的舰船可达吃水深度的7倍。尾流达到最大厚度后,气泡主要受自身浮力的影响,舰船尾流的厚度随时间的增加而逐渐降低。典型的舰船尾流分布如表1所示,CSS PAR IIEAU船在航速10 kn时、经过175 s后可达9 m深[11]

      表 1  典型舰船尾流的分布[11]

      Table 1.  Wake distribution of common ships[11]

      Destroyer
      model
      Speed/knThickness
      at 360 m/m
      Width/mDraft
      Rathborn10-126.5±1.1--
      Hopewell107.7±1.2125.8
      Evian134.4±1.17.52.9
    • 尾流气泡目标特性通常指气泡数密度[12]、气泡层厚度[13]、气泡尺度[14]等特征。这里的尾流气泡层厚度特性亦可作为舰船尾流分布特性也可作为气泡目标特性。

      舰船尾流气泡数密度主要由螺旋桨工况和舰船航速所决定。一般螺旋桨转速越快,航速越快,初始数密度越大。舰船尾流气泡数密度分布比较均匀,同一深度呈高斯分布。随气泡半径大致呈线性分布,半径越小的气泡其数密度越大。将舰船尾流分为近程尾流与远程尾流,近程尾流一般为舰船到3倍舰船长度左右距离。近程舰船尾流与远程舰船尾流存在明显的气泡特性差异,近程舰船尾流场中气泡受螺旋桨搅动及舰船行驶过后形成的空穴力因素影响较大,气泡尺度较大,一般在1000 μm以下,主要集中在100~500 μm左右,根据表2实际测量结果,气泡数密度可达到106~3×108 m−3,尺度较大的气泡上浮速度较快,且气泡尺度越大,气泡数密度越低[15]

      表 2  实测舰船尾流中气泡分布密度[16]

      Table 2.  Measured bubble distribution density in ship wake[16]

      Bubble
      radius/μm
      1 min/m−33 min/m−35 min/m−3
      10704.9×1024.1×101-
      4001.79×1045.6×1037.8×102
      1604.6×1051.98×1059.3×104
      805.5×1062.61×1061.45×106
      10-10702.83×1081.32×1087.16×107

      远程舰船尾流场中气泡主要受自身浮力的影响,大尺度气泡上升到海平面消亡,小尺度气泡由于上升缓慢,在海水中存留时间比较长,尺度一般为10~300 μm,其中以40~80 μm居多。根据实际测量,海水背景中的气泡数密度一般为104~105 m−3左右,远程舰船尾流气泡数密度一般可达105~6×106 m−3。随时间的增加,远程舰船尾流数密度逐渐降低到海水背景气泡数密度,其中10~20 μm气泡数密度最大,气泡层厚度一般为2~4 m[17]

      从远程尾流到近程尾流,舰船尾流分布特性及气泡目标特性变化规律为:尾流宽度逐渐减小;气泡群集中尺度逐渐增大;气泡数密度逐渐增大;气泡层厚度逐渐增大,到达最大厚度后开始减小。

    • 蒙特卡洛方法是基于光子运动轨迹的计算机模拟,其基本过程是:光子以一个特殊的方向进入介质,确定发生碰撞时光子运动的距离,若发生散射,则由适当的散射相位函数就能选取散射后新的运动方向。若发生吸收,光子便消亡[18]。这些过程随机重复进行,直到光子被接收面所接收。文中考虑水体介质多次散射的情况,即水中散射粒子很稠密,同时考虑单次、二次及更高次的散射和衰减路径上的衰减[19],能适用于绝大多数海域。

      蒙特卡洛模拟可以分为六个步骤,(1) 发射条件:波长,能量的确定。(2) 运动轨迹(自由程):传输$ l $距离后,发生碰撞。(3) 散射过程:散射角计算。(4) 碰撞后的运动方向。(5) 新的自由程。(6) 终止条件:接收或消亡。蒙特卡洛光子模拟运动程序图如图2所示。

      图  2  光子模拟运动程序图

      Figure 2.  Program diagram of photon simulated motion

    • 不同舰船的舰船尾流目标特性不同,其回波特性也不尽相同,可利用其回波特性识别反演不同舰船,并探究其信号变化规律。文中采用多尺度、宽数密度、大厚度气泡场进行模拟,将舰船分为大型船只和小型船只,大型船只的航速高于小型船只。以远程舰船尾流气泡场作为探测的基本环境,将探测系统寻找舰船的过程分为搜索阶段和跟踪阶段,从远程尾流到近程舰船尾流过程进行模拟,近程舰船尾流模拟状态会存在于远程舰船尾流模拟中。舰船尾流是一个动态变化的环境,舰船尾流会随时间的增加,受湍流、扰动等非均匀因素影响,气泡会向四周扩散[20]。探测系统可能存在于尾流中或者尾流下,进行水下探测时探测系统以固定姿态运动,即与海平面保持固定距离,探测系统可能存在的运动方式如图3所示。

      图  3  探测可能存在的形式。 (a) 探测位于尾流之中;(b) 探测位于尾流之下

      Figure 3.  Possible forms of detection. (a) Detection in wake; (b) Detection under wake

    • 为了仿真模拟大型船只真实情况下的尾流场,仿真条件设置为激光波长532 nm,光子数为106,权值为10−5,接收口径为0.05 m,接收角为5 mrad;水质环境设为纯净海水,衰减系数为0.056 m−1,水体折射率为1.3。大型船只产生的气泡厚度在2~9 m动态变化,气泡尺度主要为20~500 μm动态变化,气泡数密度为105~3×108 m−3动态变化,气泡尺度、数密度、厚度逐渐随时间的增大而增大。仿真条件设置见表3

      表 3  大型船只仿真条件设置

      Table 3.  Simulation conditions for large vessel

      Simulation
      conditions
      Zone 1Zone 2Zone 3Zone 4Zone 5Zone 6Zone 7Zone 8
      Movement
      mode
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      Bubble size/μm20-8080-140140-200200-260260-320320-380380-440440-500
      Bubble thickness/m2-33-44-55-66-77-88-99
      Bubble density/m−3105-106 106-3×1063×106-6×1066×106-107 107-108108-2×1082×108-3×1083×108
    • 假设搜索阶段只在区域1进行,将探测搜索阶段大体分为两种情况:(1) 探测系统切入尾流场到0.5 W处;(2) 探测系统从0.5 W处切出尾流场;两种情况对称,其他情况大体包含在这两种情况中。将探测系统可能存在的位置划分为0、0.5、1、1.5、2 m,则位置划分为位置1~5。

      (1)探测系统位于尾流场之下

      x=3 m,其余条件不变,得到探测系统从水体背景中切入到尾流场0.5 W处的回波信号图,如图4(a)所示,在位置2时,水体回波凹陷,出现气泡回波,出现气泡回波时间为20~30 ns,接收光子数从2360~2795降为33,气泡回波信号逐渐左移,且气泡回波幅度逐渐增加,脉冲宽度逐渐展宽。此时,气泡回波幅度尚未饱和。情况2下的状态与情况1的状态保持对称,其情况2与此相反,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。

      (2)探测系统位于尾流场之中

      现探测系统处于尾流之中,其余条件不变。得到探测系统从水体背景中切入到尾流场0.5 W处的回波信号图,如图4(b)所示,从位置2开始出现变化,回波信号整体先上升后下降,水体初始接收光子数在降低,峰值接收光子数从2700降到2414,气泡回波信号逐渐右移,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄,较水体背景有明显变化。情况2时为气泡回波幅度逐渐增高,脉冲宽度逐渐展宽。出现气泡回波较水体回波降低的原因主要是探测系统位于舰船尾流之中,舰船尾流气泡遮挡了大部分光束,回波信号失去了部分水体的回波信号。

      图  4  信号变化趋势。(a) 尾流下; (b) 尾流中

      Figure 4.  Signal change trend. (a) Under the wake; (b) In wake

      信号强度变化最明显的位置位于0.5 W处,探测系统位于尾流之中与尾流之下探测系统信号变化相反,且探测系统位于尾流之下时,变化更为明显,若信号在一直处于某个阈值之外,则搜索失败,继续保持搜索。

    • 跟踪阶段为搜索阶段经系统检测后探测系统的下一个阶段,其已具备比较明显的变化特征,对从探测系统进入跟踪阶段到摧毁目标舰船进行模拟,将区域划分为8个区域,模拟其信号变化,仿真条件见表3

      (1)探测位于尾流场之下

      x=10 m,其余条件不变,其跟踪阶段回波信号图见图5(a),气泡回波信号逐渐左移、幅度持续凹陷、幅度缓慢增加,脉冲宽度逐渐展宽,接收光子数从1 012~2 721范围降为0,出现气泡信号时间范围为22~85 ns,气泡回波自区域2开始饱和,至区域7时,恢复不到原来水体回波信号,较搜索阶段信号变化明显。

      (2)探测系统位于尾流场之中

      现探测系统在尾流中,其余条件不变。其跟踪阶段回波信号图如图5(b)所示,信号初始接收光子数在降低,在区域2时初始接收光子数为零,大型船只峰值接收光子数从2453降到1683,气泡回波信号整体右移,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄,较搜索阶段有明显变化。

      图  5  信号变化趋势。 (a)尾流下; (b) 尾流中

      Figure 5.  Signal change trend. (a) Under the wake; (b) In wake

      探测系统位于尾流之中与尾流之下探测系统信号变化相反,且探测系统位于尾流之下时,变化更为明显,若信号在一直处于某个阈值之外,则跟踪失败,重新保持搜索。

    • 试验目标设置为小型船只,小型船只产生的气泡厚度在1~4 m之间动态变化,气泡数密度在104~107 m−3之间变化,其余条件不变,仿真条件如表4所示。

      表 4  小型船只仿真条件设置

      Table 4.  Small vessel simulation condition settings

      Simulation
      conditions
      Zone 1Zone 2Zone 3Zone 4
      Movement
      mode
      In/under
      wake
      In/under
      wake
      In/under
      wake
      In/under
      wake
      Bubble size/μm20-6060-100100-300300-500
      Bubble
      thickness/m
      0-11-22-33-4
      Bubble
      density/m−3
      104-105 105-106 106-107 107
    • 将探测系统可能存在的位置划分为0、0.25、0.5、0.75、1 m,则位置划分为位置1~5,仿真环境为在区域1进行,其余条件不变。

      (1)探测系统位于尾流场之下

      x=3 m,得到探测系统从水体背景中切入到尾流场0.5 W处的回波信号图,如图6(a)所示,信号变化趋势与大型船只变化相同,出现气泡时接收光子数从2443降为323,出现气泡时间范围为7~17 ns,信号变化强度不如大型船只,幅值和脉宽变化较小。

      图  6  信号变化趋势。(a)尾流下; (b) 尾流中

      Figure 6.  Signal change trend. (a) Under the wake; (b) In wake

      (2)探测系统位于尾流场之中

      现探测系统处于尾流之中,得到探测系统从水体背景中切入到尾流场0.5 W处的回波信号图,如图6(b)所示,信号变化趋势与大型船只变化相同,峰值接收光子数从2812降为2482,变化强度不如大型船只,幅值和脉宽变化同样较小。

    • 将小型船只跟踪区域划分为四个区域,对四个区域的信号变化进行模拟。仿真条件如表4所示。

      (1)探测系统位于尾流场之下

      x=10 m,其余条件不变,其跟踪阶段回波信号图见图7(a),气泡变化趋势与大型船只相同,但出现气泡回波时间较大型船只晚,在50 ns后才出现,出现气泡回波接收光子数从856~1 362范围降为0,出现气泡时间范围为11~33 ns,变化强度不如大型船只。

      (2)探测系统位于尾流场之中

      现探测系统在尾流中,其余条件不变。其跟踪阶段回波信号图如图7(b)所示,气泡变化趋势与大型船只相同,但幅度和脉宽变化范围不如大型船只,小型船只接收光子数从2510降到了1 960左右,区域4时初始接收光子数为0。

      图  7  信号变化趋势。 (a)尾流下; (b) 尾流中

      Figure 7.  Signal change trend. (a) Under the wake; (b) In wake

      通过仿真分析可得:当激光探测系统位于尾流之下时,大型船只的气泡回波幅度缓慢上升,气泡脉冲宽度显著展宽,且距舰船越近,气泡回波变化越明显。激光探测系统位于尾流之中时,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。激光探测系统位于尾流之下时与探测系统位于尾流之中时,信号变化相反。小型船只信号变化趋势基本与大型船只保持一致,但尾流激光探测回波强度变低。

    • 试验地点选择某湖泊,其环境水质清澈,水质浊度为0.01 NTU,可模拟纯净海水水质。设备采用离轴探测模式,样机选取的激光脉冲能量为(30±1) mJ,重频1~10 Hz可调,波束角为2~10 mrad可调,脉冲宽度为8~12 ns,APD接收增益为0~30 dB可调。将探测系统密封至罐中,由固定装置将其拖入水中,通过水密电缆将回波信息传递回电脑端。目标舰船1设置为运砂船,运砂船上满载砂子,可模拟大型船只;目标舰船2设置为游艇,可模拟小型船只。探测系统设置如图8(a)~(c)所示。探测系统在目标船驶过后,切入尾流中,对舰船尾流进行实时采集,采集时间为3 min,直至逼近目标船,试验过程如图8(d)、(e)所示。

      图  8  试验系统及试验过程。(a)探测系统;(b)试验船;(c)探测系统;(d)湖泊测试跟踪过程;(e)跟踪过程示意图

      Figure 8.  Experimental system and process. (a) Detection system; (b) Experimental ship; (c) Detection system; (d) Lake test tracking pro-cess; (e) Schematic diagram of tracking process

    • 首先,将系统放置于水下3 m位置,此时系统处于舰船尾流之中,水体激光散射背景信号与存在尾流时的激光散射回波信号如图9(a)所示,将左图进行平滑处理,得到处理后的数据图如图9(b)所示,可直观地显示其最大、最小值的变化。如图9(a)、(b)所示,真实舰船尾流激光散射回波信号相对背景水体散射信号变化明显,回波信号幅值逐渐降低,峰值幅度变化0.33 v,脉冲宽度略微降低。将其与仿真进行对比,探测系统切入尾流场后,仿真和试验所得信号幅值都在减小,脉宽都在变窄,将仿真与试验所得信号归一化处理后,两者幅值最大变化分别为0.516和0.598,试验所得信号变化较仿真变化更强。出现幅值降低的原因是激光束向上发射到尾流气泡后,探测系统位于舰船尾流之中,舰船尾流气泡遮挡了大部分光束,回波信号失去了部分水体的回波信号。

      图  9  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

      Figure 9.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

      现将系统放置于水下6 m位置,此时系统处于舰船尾流之下,水体激光散射背景信号与存在尾流时的激光散射回波信号如图10(a)、(b)所示,如图所示,真实舰船尾流激光散射回波信号相对背景水体散射信号变化明显,后向散射信号幅值逐渐增高,信号幅值变化0.42 V,脉冲宽度逐渐增加。将其与仿真进行对比,探测系统切入尾流场后,仿真和试验所得信号幅值一直都在升高,但之后回波信号限幅,幅值都不再增加,回波信号脉宽开始展宽,将仿真与试验所得信号归一化处理后,两者幅值最大变化分别为0.362和0.365,基本保持一致,但试验所得更强一些。

      图  10  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号;(b)平滑处理后的信号

      Figure 10.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    • 将系统放置于水下3 m位置,系统处于舰船尾流之中,水体激光散射背景信号与存在尾流时的激光散射回波信号如图11(a)、(b)所示,如图所示,真实舰船尾流激光散射回波信号相对背景水体散射信号变化不明显,基本拟合在一块。出现这种情况的原因主要是由于舰船尾流较为稀疏,目标回波微弱。

      将系统放置于水下6 m位置,系统处于舰船尾流之下,水体激光散射背景信号与存在尾流时的激光散射回波信号如图12(a)、(b)所示,真实舰船尾流激光散射回波信号相对背景水体散射信号比较明显,峰值变化0.22 V,后向散射信号幅值逐渐降低,脉冲宽度略微降低。

      图  11  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

      Figure 11.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

      通过试验得到:探测系统位于尾流之下时,运砂船尾流激光回波信号相对背景水体背景激光回波信号变化明显,测量点距船舶越近,激光回波信号强度越高,游艇尾流激光探测幅度基本无变化。探测系统位于尾流之中时,大小船只尾流激光回波信号相对背景水体散射信号变化都较明显,激光探测系统都可对尾流实现有效探测。同时将运砂船、游艇尾流激光探测回波对比,运砂船尾流激光探测系统灵敏度高。

      图  12  水体背景激光散射信号与尾流激光散射回波信号对比。(a)原始信号; (b)平滑处理后的信号

      Figure 12.  Water background laser scattering signal and bubble laser scattering echo signal. (a) Original signal; (b) Smoothed signal

    • 文中以舰船尾流激光探测工程应用为背景,分析了舰船尾流分布特性及气泡目标特性的变化规律。基于舰船尾流分布特性及目标特性建立仿真环境,采用蒙特卡洛模拟方法,对多尺度、宽数密度、大厚度舰船尾流气泡群进行模拟仿真,通过对相应情况下舰船尾流后向散射回波信号总结分析,可有效模拟探测系统在舰船尾流目标特性下的真实状态,研究了探测系统在搜索、跟踪阶段的信号变化趋势,以及不同目标舰船的回波信号变化强度。通过分析可得:当激光探测系统位于尾流之下时,大型船只的气泡回波幅度缓慢上升,气泡脉冲宽度显著展宽,且距舰船目标越近,气泡回波变化越明显。激光探测系统位于尾流之中时,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。激光探测系统位于尾流之下时与探测系统位于尾流之中时,信号变化相反。小型船只信号变化趋势基本与大型船只保持一致,但尾流激光探测回波强度变低。搭建外场激光后向散射回波试验系统,验证探测系统在尾流之下时,气泡回波信号变化为气泡幅度缓慢上升,气泡脉冲宽度显著展宽;探测系统位于尾流之中时,气泡回波幅度逐渐降低,脉冲宽度逐渐变窄。这可为舰船尾流探测在实际工程应用中提供支撑。

参考文献 (20)

目录

    /

    返回文章
    返回