高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

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

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

海洋激光雷达漫射衰减系数反演方法研究

雷子昂 杨颂 沈振民 孙倩 张景豪 郑永超

引用本文:
Citation:

海洋激光雷达漫射衰减系数反演方法研究

    通讯作者: 郑永超, laser0371@163.com
  • 中图分类号: TN958.98;O436

Research on inversion method of diffuse attenuation coefficient of ocean LiDAR

    Corresponding author: ZHENG Yongchao, laser0371@163.com
  • CLC number: TN958.98;O436

  • 摘要: 为了解决漫射衰减系数Kd的反演依赖于原位测量数据的问题,提出了一种基于Klett和Fernald法反演Kd的融合算法。该融合算法以Klett法反演结果为基准,计算出激光雷达比作为Fernald法迭代时的先验信息,再用Fernald法反演出更加精确的水体漫射衰减系数Kd; 并运用数学解析模型作为融合算法的仿真数据源,解决了实测信号数据量不足的问题,最后基于自研的双频激光雷达系统,利用海试数据对融合算法进行了实验验证。结果表明,融合算法反演得到Kd值的均方根误差值为0.74 m-1,精度优于传统的Klett法和Fernald法;所提出的融合算法能够满足近岸浅海水体漫射衰减系数Kd反演的收敛和精度要求,且能在没有原位现场环境先验信息的条件下快速计算激光雷达比。该研究对浅海水域的深度测量以及海洋剖面光学参数的测量是有帮助的。
  • 图 1  融合算法步骤图

    Figure 1.  Flow diagram of fusion algorithm

    图 2  海洋激光雷达系统示意图

    Figure 2.  Schematic of the marine LiDAR system

    图 3  实测和建模的激光雷达回波信号对比图

    Figure 3.  Comparison chart of measured and simulated LiDAR echo signals

    图 4  仿真回波信号

    Figure 4.  Simulated echo signal

    图 5  Klett算法和不同激光雷达比的Fernald法仿真反演结果

    Figure 5.  Inversion results of Klett algorithm and Fernald method with different LiDAR ratios in simulation

    图 6  融合算法反演Kd值仿真结果

    Figure 6.  Diagram of inversion Kd value of fusion algorithm in simulation

    图 7  Klett算法和不同激光雷达比的Fernald法实测反演结果

    Figure 7.  Inversion results of Klett algorithm and Fernald method with different LiDAR ratios in experiment

    图 8  融合算法反演Kd值实测结果

    Figure 8.  Diagram of inversion Kd value of fusion algorithm in experiment

    表 1  激光雷达系统参数

    Table 1.  LiDAR system parameters

    parameter value parameter value
    wavelength 532 nm/1064 nm pulse energy 3 mJ
    pulse width 10 ns repetition frequency 10 Hz
    excessive noise factor 3 electronic bandwidth 300 MHz
    input optical efficiency 0.5 dark current noise 10-8 A
    sample interval 0.2 ns sample rate 5 GHz
    下载: 导出CSV

    表 2  解析模型的数值参数

    Table 2.  Numerical parameters for simulation models

    parameter value parameter value
    λ 532 nm β 0.0015/(m·sr)
    H 200 m Is 0.025 W/(m2·sr)
    Pe 2×107 W nw 1.34
    T0 10 ns G 3
    Tatm2 0.9 Id 10-8 A
    Ar 0.025 m2 Rλ 0.3 A/W
    ne 0.9 Δλ 1 nm
    nr 0.5 γr 0.15
    F 1 B 300 MHz
    e 1.6×10-19 C E0 3 mJ
    下载: 导出CSV
  • [1] 黄田程, 陶邦一, 毛志华, 等. 基于多通道海洋激光雷达的海陆波形分类[J]. 中国激光, 2017, 44(6): 0610002.

    HUANG T Ch, TAO B Y, MAO Zh H, et al. Classification of sea and land waveform based on multi-channel ocean lidar[J]. Chinese Journal of Lasers, 2017, 44(6): 0610002(in Chinese). 
    [2] 余雪敏. 用航空遥感器对海岸带海洋水色的分析[J]. 航天返回与遥感, 1994, 15(2): 50-65.

    YU X M. Analysis of coastal marine water color with aerial remote sensors[J]. Spacecraft Recovery & Remote Sensing, 1994, 15(2): 50-65(in Chinese). 
    [3] 李德仁, 王密. 高分辨率光学卫星测绘技术综述[J]. 航天返回与遥感, 2020, 41(2): 1-11.

    LI D R, WANG M. A review of high resolution optical satellite surveying and mapping technology[J]. Spacecraft Recovery & Remote Sensing, 2020, 41(2): 1-11(in Chinese). 
    [4] 仝迟鸣, 鲍云飞, 黄巧林, 等. 太阳诱导叶绿素荧光卫星遥感技术研究进展[J]. 航天返回与遥感, 2022, 43(2): 45-55.

    TONG Ch M, BAO Y F, HUANG Q L, et al. Research progress of solar induced chlorophyll fluorescence satellite remote sensing techno-logy[J]. Spacecraft Recovery & Remote Sensing, 2022, 43(2): 45-55(in Chinese).
    [5] 刘志鹏, 刘东, 徐沛拓, 等. 海洋激光雷达反演水体光学参数[J]. 遥感学报, 2019, 23(5): 944-951.

    LIU Zh P, LIU D, XU P T, et al. Marine LiDAR inverts the optical parameters of water[J]. Journal of Remote Sensing, 2019, 23(5): 944-951(in Chinese). 
    [6] 周雨迪. 用于水体光学特性探测的海洋激光雷达研究[D]. 杭州: 浙江大学, 2020: 41-58.

    ZHOU Y D. Marine LiDAR research for the detection of optical pro-perties of water[D]. Hangzhou: Zhejiang University, 2020: 41-58 (in Chinese).
    [7] 刘航. 基于激光雷达回波的海洋光学参数反演研究[D]. 杭州: 浙江大学, 2020: 31-63.

    LIU H. Research on inversion of marine optical parameters based on LiDAR echoes[D]. Hangzhou: Zhejiang University, 2020: 31-63 (in Chinese).
    [8] 徐沛拓, 陶雨婷, 刘志鹏, 等. 海洋激光雷达实验与仿真结果的对比[J]. 红外与激光工程, 2020, 49(2): 0203007.

    XU P T, TAO Y T, LIU Zh P, et al. Comparison of oceanic lidar experiments and simulation results[J]. Infrared and Laser Engineering, 2020, 49(2): 0203007(in Chinese). 
    [9]

    FERNALD F G. Analysis of atmospheric LiDAR observations: Some comments[J]. Applied Optics, 1984, 23(5): 652-653. doi: 10.1364/AO.23.000652
    [10]

    FERNALD F G, HERMAN B M, REAGAN J A. Determine of aerosol height distribution by LiDAR[J]. Journal of Applied Meteorology, 1972, 11(3): 482-489. doi: 10.1175/1520-0450(1972)011<0482:DOAHDB>2.0.CO;2
    [11]

    KLETT J D. Stable analytical inversion solution for processing LiDAR returns[J]. Applied Optics, 1981, 20(2): 211-220. doi: 10.1364/AO.20.000211
    [12]

    LIU R, LING Q, ZHANG Q, et al. Detection of chlorophyll a and CDOM absorption coefficient with a dual-wavelength oceanic LiDAR: Wavelength optimization method[J]. Remote Sensing, 2020, 12(18): 3021. doi: 10.3390/rs12183021
    [13] 李凯鹏, 贺岩, 侯春鹤, 等. 双波长海洋激光雷达探测近岸到大洋水体的叶绿素剖面[J]. 中国激光, 2021, 48(20): 2010002.

    LI K P, HE Y, HOU Ch H, et al. Dual-wavelength ocean LiDAR detects chlorophyll profiles from nearshore to oceanic water[J]. Chinese Journal of Lasers, 2021, 48(20): 2010002(in Chinese). 
    [14]

    CHURNSIDE J H. LiDAR detection of plankton in the ocean[C]// IEEE International Geoscience & Remote Sensing Symposium. Barcelona, Spain: IEEE, 2007: 3174-3177.
    [15]

    HALTRIN V I. One-parameter two-term Henyey-Greenstein phase function for light scattering in seawater. Applied Optics, 2002, 41(6): 1022-1028. doi: 10.1364/AO.41.001022
    [16]

    ABDALLAH H, BAGHDADI N, BAILLY J S, et al. Wa-LiD: A new LiDAR simulator for waters[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 744-748. doi: 10.1109/LGRS.2011.2180506
    [17]

    ZHOU Y, CHEN W, CUI X, et al. Validation of the analytical model of oceanic LiDAR Returns: Comparisons with Monte Carlo simulations and experimental results[J]. Remote Sensing, 2019, 11(16): 1870. doi: 10.3390/rs11161870
    [18]

    KOPLLEYICH Y I. FEYGLES V I, SURKOV A I. Mathematical modeling of input signals for oceanographic lidar systems[J]. Proceedings of the SPIE, 2003, 5155: 30-39. doi: 10.1117/12.506980
    [19] 朱峻可, 李丽娟, 林雪竹. 激光雷达测量系统的测量场规划研究[J]. 激光技术, 2021, 45(1): 99-105.

    ZHU J K, LI L J, LIN X Zh. Research on the measurement field planning of LiDAR measurement system[J]. Laser Technology, 2021, 45(1): 99-105(in Chinese). 
    [20] 张熠星. 机载测深激光雷达的海底回波提取技术[D]. 上海: 东华大学, 2018: 32-60.

    ZHANG Y X. Sea bottom echo extraction technology of airborne sounding LiDAR[D]. Shanghai: Donghua University, 2018: 32-60 (in Chinese).
    [21] 薛文佳. 机载激光雷达回波信号建模及处理技术研究[D]. 天津: 天津大学, 2018: 25-42.

    XUE W J. Research on modeling and processing technology of airborne laser radar echo signal[D]. Tianjin: Tianjin University, 2018: 25-42 (in Chinese).
    [22]

    SULLIVAN J M, TWARDOWSKI M S. Angular shape of the oceanic particulate volume scattering function in the backward direction. Applied Optics, 2009, 48(35): 6811-6819. doi: 10.1364/AO.48.006811
    [23] 陈玉宝, 王箫鹏, 步志超, 等. 超大城市试验气溶胶激光雷达标定及结果分析[J]. 激光技术, 2022, 46(4): 435-440.

    CHEN Y B, WANG X P, BU Zh Ch, et al. Calibration and result analysis of aerosol LiDAR in megacity experiment[J]. Laser Technology, 2022, 46(4): 435-440(in Chinese). 
    [24]

    LEE J H, CHURNSIDE J H, MARCHBANKS R D, et al. Oceanographic lidar profiles compared with estimates from in situ optical measurements. Applied Optics, 2013, 52(4): 786-794. doi: 10.1364/AO.52.000786
    [25]

    WANG L H. JACQUES S L, ZHENG L Q. MCML—Monte Carlo modeling of light transport in multi-layered tissue[J]. Computer Methods and Programs in Biomedicine, 1995, 47(2): 131-146. doi: 10.1016/0169-2607(95)01640-F
  • [1] 卢益民罗志祥陈文革杨宗凯 . 海洋激光雷达系统的噪声相空间重建的研究. 激光技术, 1998, 22(4): 199-202.
    [2] 卢益民罗志祥陈文革杨宗凯 . 机载海洋激光雷达信号中噪声的混沌行为初探. 激光技术, 1998, 22(3): 139-143.
    [3] 胥克涛袁扬胜冯霞屈军 . PCFT光束阵列在海洋湍流中的传输特性. 激光技术, 2015, 39(6): 877-884. doi: 10.7510/jgjs.issn.1001-3806.2015.06.031
    [4] 闫顺生李海花 . 利用喇曼-米方法确定激光雷达几何因子. 激光技术, 2013, 37(4): 511-514. doi: 10.7510/jgjs.issn.1001-3806.2013.04.021
    [5] 张艳红卢腾飞刘永欣陈子阳孙顺红 . 非均匀偏振光束在海洋湍流中的光强特性. 激光技术, 2020, 44(3): 310-314. doi: 10.7510/jgjs.issn.1001-3806.2020.03.007
    [6] 陈文革黄铁侠卢益民 . 机载海洋激光雷达发展综述. 激光技术, 1998, 22(3): 147-152.
    [7] 王元庆何炬马泳余寅张念梁琨 . 随机噪声对海洋布里渊激光雷达测量的影响. 激光技术, 2015, 39(1): 6-12. doi: 10.7510/jgjs.issn.1001-3806.2015.01.002
    [8] 马泳冀航梁琨林宏 . 载波调制激光雷达技术在海洋探测中的应用. 激光技术, 2008, 32(4): 346-349.
    [9] 王玉诏 . 一种用于激光雷达信号的自适应分段平滑算法. 激光技术, 2014, 38(5): 599-602. doi: 10.7510/jgjs.issn.1001-3806.2014.05.005
    [10] 戴璨王元庆徐帆 . 基于粒子群算法的3维激光雷达回波分解. 激光技术, 2016, 40(2): 284-287. doi: 10.7510/jgjs.issn.1001-3806.2016.02.028
    [11] 黄敏王玉兰王娜陈涌任鹏周鼎富史晓丁冯力天 . 机载测风激光雷达下视VAD反演及算法仿真. 激光技术, 2012, 36(1): 22-25,41. doi: 10.3969/j.issn.1001-3806.2012.01.007
    [12] 徐文静孙东松舒志峰唐磊董吉辉胡冬冬王国成 . 三通道瑞利散射测风激光雷达风速反演算法. 激光技术, 2011, 35(4): 481-485,491. doi: 10.3969/j.issn.1001-3806.2011.04.011
    [13] 徐强沈思谢修敏吴鹏周强邓光伟王浟宋海智 . 可用于激光雷达的量子光学技术. 激光技术, 2021, 45(1): 44-47. doi: 10.7510/jgjs.issn.1001-3806.2021.01.008
    [14] 魏华江邢达鲁建军沈行燕巫国勇 . 直接和间接法测量人膀胱对激光的总衰减系数. 激光技术, 2005, 29(4): 420-422.
    [15] 彭长华译邹声荣校 . 激光雷达专题文摘. 激光技术, 1985, 9(6): 62-65.
    [16] 王春勇谢俊卞保民李振华 . 提高激光雷达测距精度的最佳信号门限比. 激光技术, 2007, 31(4): 408-411.
    [17] 彭涛王茜石磊姜勇陈涌周鼎富 . 测风激光雷达光学参量对相干效率的影响. 激光技术, 2023, 47(6): 751-756. doi: 10.7510/jgjs.issn.1001-3806.2023.06.004
    [18] 伯广宇钟志庆王邦新刘东 . 喇曼激光雷达探测云及其附近气溶胶光学参量. 激光技术, 2012, 36(5): 597-601. doi: 10.3969/j.issn.1001-3806.2012.05.006
    [19] 蒋佳佳沈法华周慧杨亮亮仇成群徐华 . 测风激光雷达双棱镜2维扫描系统的光学设计. 激光技术, 2023, 47(3): 425-431. doi: 10.7510/jgjs.issn.1001-3806.2023.03.022
    [20] . 国外激光雷达之发展. 激光技术, 1983, 7(5): 1-7.
  • 加载中
图(8) / 表(2)
计量
  • 文章访问数:  169
  • HTML全文浏览量:  103
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-04-27
  • 录用日期:  2023-06-06
  • 刊出日期:  2024-05-25

海洋激光雷达漫射衰减系数反演方法研究

    通讯作者: 郑永超, laser0371@163.com
  • 1. 中国空间技术研究院 北京空间机电研究所,北京 100094,中国
  • 2. 中国空间技术研究院 空间激光信息感知技术核心专业实验室,北京 100094,中国

摘要: 为了解决漫射衰减系数Kd的反演依赖于原位测量数据的问题,提出了一种基于Klett和Fernald法反演Kd的融合算法。该融合算法以Klett法反演结果为基准,计算出激光雷达比作为Fernald法迭代时的先验信息,再用Fernald法反演出更加精确的水体漫射衰减系数Kd; 并运用数学解析模型作为融合算法的仿真数据源,解决了实测信号数据量不足的问题,最后基于自研的双频激光雷达系统,利用海试数据对融合算法进行了实验验证。结果表明,融合算法反演得到Kd值的均方根误差值为0.74 m-1,精度优于传统的Klett法和Fernald法;所提出的融合算法能够满足近岸浅海水体漫射衰减系数Kd反演的收敛和精度要求,且能在没有原位现场环境先验信息的条件下快速计算激光雷达比。该研究对浅海水域的深度测量以及海洋剖面光学参数的测量是有帮助的。

English Abstract

    • 海洋光学参数与海水中的叶绿素、悬浮物、赤潮、营养盐等物质密切相关,海洋光学参数的时空变化对研究海洋碳循环、海洋生态有重要意义[1]。漫射衰减系数Kd是海洋光学参数中常用的一个表观光学量。目前用于测量漫射衰减系数Kd的方法主要包括原位测量和被动光学遥感。原位测量漫射衰减系数Kd虽然有很高的探测精度,但探测效率低,对极浅水域以及暗礁区域等船只无法抵达的区域探测存在盲区。被动光学遥感依靠飞机、卫星等平台,探测效率高,覆盖区域广,但是自然光穿水深度有限且缺乏时间维度信息,只能对30 m以浅的水体进行一个积分探测。近年来,激光雷达(ligth detection and ranging, LiDAR)用于海洋光学参数测量逐渐受到关注,由于其发射的蓝绿波段激光对海洋水体的穿透性高,回波携带时间维度信息能对海洋剖面进行直接探测,且搭载在飞机、卫星等平台,具备高时空分辨率、高效率、无盲区等优点,能够与原位测量和被动光学遥感技术形成相互补充,丰富海洋光学观测的数据层次[2-4]

      海洋激光雷达测量海洋光学参数的理论基础是海洋激光雷达方程。海洋激光雷达方程中的后向体散射系数β与衰减系数α这两个物理量表征光学参数且共同影响激光回波功率值,因而难以通过一个海洋雷达方程同时求解这两个物理量[5-8]。海洋光学参数的反演受水质因素影响显著,常用的方法包括Fernald法和Klett法[9-11]。Fernlad算法广泛适用于各类水质的水体,但其精度依赖于水质的先验信息,不同水质的水体需要选择不同的激光雷达比[9]。LIU等人基于Fernald法对漫射衰减系数Kd进行了反演研究,采用原位测量数据作为最优激光雷达比的确定依据[5]。但是全球海洋的原位探测数据稀疏不均匀,并且实时性不足,对全球大多数海区往往缺乏待测海域的水质先验信息和原位探测数据[10]。LIU等人通过Klett法研究了近岸和大洋水体叶绿素浓度剖面[12-15]。Klett法在浑浊水体条件下Kd值的求解与激光雷达无关,但是对Kd值的求解精度低于Fernald法[11]

      本文作者提出了一种反演水体漫射衰减系数Kd的融合算法,该方法适用于近岸浑浊水体,通过Klett法的Kd计算结果来选择待测水体的激光雷达比,再基于Fernald法对Kd进行精确求解,这种方法不依赖于原位测量数据,能够广泛应用于全球近岸水体的Kd测量,并且反演精度优于单一Klett法的反演结果。为了验证算法的有效性,采用自主研发的双频(近红外、绿光)海洋激光雷达,在某水域开展了近岸浅海水域的光学参数垂直廓线探测试验。本次实验中采用的是岸基平台,水体的水质条件单一,激光雷达的测点密度稀疏,空间离散型大,因而需要大量的仿真数据作为Kd反演算法的验证数据源。使用Walid解析模型产生不同水质条件下的激光雷达回波波形作为仿真数据源,并通过相同水质条件下的激光雷达回波信号对比验证Walid模型数据的有效性。最后,将融合算法用于仿真波形和实测波形的反演验证,结果证明了融合算法具备高精度反演Kd的能力。

    • 海洋激光雷达方程用于描述激光在海水中的辐射传输特性,激光器向海洋发射蓝绿激光脉冲,经大气、气海界面、海水介质,最后被探测器接收。海洋激光雷达方程中表征水体光学特性主要参数是海水后向体散射系数β和激光雷达衰减系数α,把其它常数因子归集为因子A,因此激光雷达方程可以表示为[6]

      $ \begin{array}{l} P(z)=A\left[\beta_{\mathrm{p}}(z)+\beta_{\mathrm{w}}(z)\right] \exp \{-2 \times \\ \left.\int_0^z\left[\alpha_{\mathrm{p}}\left(z^{\prime}\right)+\alpha_{\mathrm{w}}\left(z^{\prime}\right)\right] \mathrm{d} z^{\prime}\right\}(n H+z)^{-2} \end{array} $

      (1)

      式中: $A$为激光雷达系统常数, $A=P_{\mathrm{e}} T_{\mathrm{atm}}{ }^2 A_{\mathrm{r}} \eta_{\mathrm{e}} \eta_{\mathrm{r}} F(1-$ $\left.L_{\mathrm{s}}\right)^2 ; P_{\mathrm{e}}$为激光发射功率; $T_{\mathrm{atm}}{ }^2$为双向大气衰减; $A_{\mathrm{r}}$为接收望远镜面积; $\eta_{\mathrm{e}}$和$\eta_{\mathrm{r}}$为出射光学系数和人射光学系数; $F$为望远镜视场损失因子; $L_{\mathrm{s}}$为海表传输损失; $\beta(z)$为后向体散射函数; $\alpha(z)$为激光雷达衰减函数, 下标$\mathrm{p}$和$\mathrm{w}$表示水体悬浮颗粒物和海水分子; $z$为海水深度; $z^{\prime}$表示积分里对应的$z$值; $n$为海水折射率; $H$为距水面高度。当激光雷达系统的接收视场足够大, 水体的漫射衰减系数$K_{\mathrm{d}}$近似等于激光雷达的衰减系数[12]。因此, 激光雷达衰减系数$\alpha$的测量结果等同于水体漫射衰减系数$K_{\mathrm{d}}$ :

      $ K_{\mathrm{d}} \approx \alpha $

      (2)
    • 海洋水体漫射衰减系数Kd是重要的海洋光学参数,可以用于初级生产力、叶绿素含量、颗粒有机碳含量等的评估[15]。从一个海洋激光雷达方程中求解两个未知光学参数存在困难,因此在反演漫射衰减系数Kd时需要设计合适的算法。目前比较成熟的反演算法为Fernald法和Klett法。

      Fernald法需要对激光雷达比有先验信息,即假设的激光雷达比为常数,具体表示为:

      $ S=\alpha / \beta $

      (3)

      式中: S为所假设的激光雷达比。当假设的激光雷达比选取准确时,反演激光雷达的衰减系数精度更高。但是不同区域水体悬浮物的激光雷达比差异大,对经验的依赖度高, 其迭代格式如下所示[9]

      $ \alpha_{\mathrm{p}}(I-1)+\alpha_{\mathrm{w}}(I-1) S_{\mathrm{p}} / S_{\mathrm{w}}=\frac{X(I-1) \exp [+K(I-1, I)]}{\frac{X(I)}{\alpha_{\mathrm{p}}(I)+\alpha_{\mathrm{w}}(I) S_{\mathrm{p}} / S_{\mathrm{w}}}+\{X(I)+X(I-1) \exp [+K(I-1, I)]\} \Delta Z} $

      (4)

      式中: SpSw分别表示水体悬浮物激光雷达比和海水分子激光雷达比,本文中纯水光学特性为已知量,纯水激光雷达比Sw近似为216 sr; αw(I)近似采用纯水的漫射衰减系数为0.0519 m-1[8]; X(I)和X(I-1)均表示距离修正的回波信号; I为迭代次数; K(I-1, I)=(SwSp) ΔZ; ΔZ为每次迭代距离间隔。针对近岸浅海水域,考虑到当αp(I)$\gg $αw(I)时,因此可以忽略αw(I)Sp/Sw项的影响,则式(4)可以改写成式(5),得到Klett法迭代格式[11]

      $ \alpha(I-1)=\frac{X(I-1)}{\frac{X(I)}{\alpha_{\mathrm{p}}(I)}+[X(I)+X(I-1)] \Delta Z} $

      (5)

      Fernald算法是目前比较成熟的用于反演漫射衰减系数Kd的后向迭代算法,影响其算法精度的主要因素就是所假设的激光雷达比。Klett算法没有分开考虑纯海水和海水悬浮颗粒物,理论上其精度不如Fernald算法,但是在近岸海域,水体悬浮颗粒物的Kd值远大于水分子的Kd值时,求解出的水体悬浮颗粒物的漫射衰减系数Kd的精度接近于Fernald算法,且在反演漫射衰减系数时不依赖于所假设的激光雷达比。因此, 综合利用两种算法的优点,针对近岸浅海水域,本文作者提出了一种基于Fernald和Klett算法的融合算法,算法流程如图 1所示。该融合算法以Klett法反演结果为基准,计算出激光雷达比作为Fernald法迭代时的先验信息,最后用Fernald法反演出更加精确的水体漫射衰减系数Kd。具体来说:(a)输入标定值α0,基于Klett方法求解出一个精度较粗略的激光雷达衰减系数αk; (b)基于Fernald法中设置不同的水体悬浮物激光雷达比Si作为输入条件,计算出不同激光雷达比条件下对应的激光雷达衰减系数αi; (c)将其中与αk最接近的αi记作αf,与其对应的Sf作为最佳激光雷达比; (d)将最优的激光雷达比Sfαf作为初始值代入Fernald法,迭代计算出激光雷达的衰减系数α,求解出水体光学参数Kd

      图  1  融合算法步骤图

      Figure 1.  Flow diagram of fusion algorithm

    • 由于海试数据测点密度稀疏,空间离散型大,为了验证融合算法的有效性,需要大量的仿真数据作为融合算法反演Kd的验证数据源。近年来,学者们发展了不同形式的海洋激光雷达方程来建立激光雷达回波功率与海水深度、海水衰减系数、后向散射系数之间的联系[16-19]。在目前常用的海洋激光雷达回波数学模型中,Walid模型全面考虑了发射激光的脉冲展宽、激光雷达系统参数、海水特性和各种噪声,适合用此模型作为验证融合算法有效性的仿真数据源,下面简要介绍Walid模型的表达形式。

      $ P_{\mathrm{t}}(t)=P_{\mathrm{s}}(t)+P_{\mathrm{c}}(t)+P_{\mathrm{b}}(t)+P_{\mathrm{g}}(t)+P_{\mathrm{n}}(t) $

      (6)

      式中: $P_{\mathrm{t}}(t)$为接收回波功率总功率; $P_{\mathrm{s}}(t)$为海表回波功率; $P_{\mathrm{c}}(t)$为海水回波功率; $P_{\mathrm{b}}(t)$为海底回波功率; $P_{\mathrm{g}}(t)$为背景噪声的回波功率; $P_{\mathrm{n}}(t)$为探测器噪声回波功率; $t$为时间范围。本文中研究的水体光学参数主要影响海洋激光雷达的海水回波, 海水回波的表达形式可以看作是高斯发射光束与水深$z$处海水回波功率的卷积[20-22], 具体形式为:

      $ \begin{array}{c} P_{\mathrm{c}}(t)=w\left(t_{\mathrm{c}}\right) * \\ P_{\mathrm{e}} T_{\mathrm{atm}}{ }^2 A_{\mathrm{r}} \eta_{\mathrm{e}} \eta_{\mathrm{r}} F\left(1-L_{\mathrm{s}}\right)^2 \beta(z) \times \\ \exp \left[-2 \int_0^z \alpha\left(z^{\prime}\right) \mathrm{d} z^{\prime}\right](n H+z)^{-2} \end{array} $

      (7)

      其中定义w(tc)为:

      $ w\left(t_{\mathrm{c}}\right)=\frac{2}{T_0} \sqrt{\frac{\ln 2}{\pi}} \exp \left[-4(\ln 2) \frac{\left(t-t_{\mathrm{c}}\right)^2}{T_0{ }^2}\right] $

      (8)

      式中: tc为高斯光束从海水到探测器的双程时间; T0为脉冲宽度。

      海水中太阳辐射背景噪声回波功率被定义为一个高斯白噪声(均值为0,标准差为1),其与一个冲击响应Pg的卷积[16]:

      $ P_{\mathrm{g}}=I_{\mathrm{s}} A_{\mathrm{r}} T_{\mathrm{atm}}{ }^2\left(1-\gamma_{\mathrm{r}}{ }^2\right) \frac{\theta^4}{4} \Delta \lambda \eta_{\mathrm{e}} $

      (9)

      式中: Is为从海水中反射回来的太阳辐射;Δλ为接收光学滤波带宽;γr为接收遮拦比。

      探测器内部噪声被定义其均值为0,其标准差为σn(t)的随信号水平变化的正态分布[16]

      $ \sigma_{\mathrm{n}}(t)=\sqrt{\frac{2 e B\left[P_{\mathrm{ext}}(t) G+I_{\mathrm{d}}\right]}{R_\lambda}} $

      (10)

      式中: Pext(t)=Pg(t)+Ps(t)+Pc(t)+Pb(t);e为电子电荷;B为探测器电子带宽;G为过量噪声因素;Id为暗电流;Rλ为响应率。

    • 本文作者采用的双频海洋激光雷达系统主要由发射系统、接收系统以及综合控制处理系统三大部分构成,如图 2所示。发射系统采用532 nm/1064 nm的Nd ∶YAG双频脉冲激光器,波长为1064 nm的脉冲激光器发射的激光无法透过海水,用于判断海表的位置。波长为532 nm的脉冲激光器出射能量为3 mJ,脉宽约为10 ns,重频10 Hz,经过发射光路扩束后激光射入水中,经海水的散射和海底的反射后被探测器接收。接收系统的望远镜采用口径为200 mm的R-C(Richey-Chretien)光学系统的望远镜,视场约50 mrad。光电探测器进行距离选通和放大增益调节,由光电倍增管和雪崩二极管组成,用于接收海表通道和海水通道,以满足各个探测通道的不同探测需求。综合控制处理系统具有控制、采集、存储及处理等功能,采集卡的采样频率为5 GHz,采样间隔为0.2 ns。主要参数见表 1

      图  2  海洋激光雷达系统示意图

      Figure 2.  Schematic of the marine LiDAR system

      表 1  激光雷达系统参数

      Table 1.  LiDAR system parameters

      parameter value parameter value
      wavelength 532 nm/1064 nm pulse energy 3 mJ
      pulse width 10 ns repetition frequency 10 Hz
      excessive noise factor 3 electronic bandwidth 300 MHz
      input optical efficiency 0.5 dark current noise 10-8 A
      sample interval 0.2 ns sample rate 5 GHz
    • 解析模型的系统参数主要来自自研的激光雷达系统,海水环境参数如海表损失因子,海水的散射系数由对海域的原位测量给出,具体仿真模型参数见表 2

      表 2  解析模型的数值参数

      Table 2.  Numerical parameters for simulation models

      parameter value parameter value
      λ 532 nm β 0.0015/(m·sr)
      H 200 m Is 0.025 W/(m2·sr)
      Pe 2×107 W nw 1.34
      T0 10 ns G 3
      Tatm2 0.9 Id 10-8 A
      Ar 0.025 m2 Rλ 0.3 A/W
      ne 0.9 Δλ 1 nm
      nr 0.5 γr 0.15
      F 1 B 300 MHz
      e 1.6×10-19 C E0 3 mJ

      为了验证仿真波形的有效性,针对同一水质Kd=0.32 m-1条件下,进行了建模信号与实测信号的吻合性比较,如图 3所示。

      图  3  实测和建模的激光雷达回波信号对比图

      Figure 3.  Comparison chart of measured and simulated LiDAR echo signals

      利用Walid模型[16]仿真了该区域的激光雷达回波信号,如图 3中建模信号蓝色实线所示,海水实测回波信号如图 3中红色散点所示[23-25]。1.5 m以内的实测信号可能受到海表海浪、气泡等因素的影响,有信号失真的情况,将舍去这部分信号。建模信号在1.5 m以下能够很好地复现实测信号,Walid模型对实测水体的激光回波的模拟具有良好效果。因此,Walid建模波形能够作为反演算法的数据源。

    • 采用Walid模型在不同Kd值的情况下进行了融合算法的仿真验证分析, 以设置Kd值为0.2 m-1和0.3 m-1为例,利用融合算法对仿真波形进行了反演Kd值的验证。

      两组仿真波形如图 4所示。图中黑色标记点分别代表海表、海底回波信号位置,仿真波形分别设置的Kd值为0.25 m-1和0.3 m-1,海水深度10 m。迭代过程中受到海表脉冲展宽和海浪因素等影响,取距离海面1 m深度处为收敛值。如图 5所示, 首先利用斜率法求出海底处位置的Kd值,因为海底波形有脉冲展宽影响,因此,以约9.5 m处用斜率法得到的初始值进行后向迭代;然后使用Klett法反演出Kd值的收敛曲线,设置不同水体悬浮物的激光雷达比,得出Fernald法的反演曲线。计算出Kd=0.25 m-1处的最佳激光雷达比为160 sr;Kd=0.3 m-1处的最佳激光雷达比为200 sr。从反演曲线图可知,Kd设置为0.25 m-1处的融合算法第1步收敛值为0.254 m-1Kd设置为0.3 m-1处的融合算法第1步收敛值为0.307 m-1。利用得到收敛的Kd值和最佳激光雷达比作为Fernald法下一步迭代的初始值,最后得出融合算法的反演曲线, 如图 6所示。

      图  4  仿真回波信号

      Figure 4.  Simulated echo signal

      图  5  Klett算法和不同激光雷达比的Fernald法仿真反演结果

      Figure 5.  Inversion results of Klett algorithm and Fernald method with different LiDAR ratios in simulation

      图  6  融合算法反演Kd值仿真结果

      Figure 6.  Diagram of inversion Kd value of fusion algorithm in simulation

      图 6可知,利用融合算法得到的Kd值设置为0.25 m-1处的融合算法收敛值为0.249 m-1,融合算法得到的Kd值设置为0.3 m-1处的融合算法收敛值为0.301 m-1。使用均方根误差(root mean square error, RMSE)来衡量算法反演的精度。当Kd设置为0.25 m-1时,计算得到Fernald法的RMSE为0.105 m-1,Klett法的RMSE为0.103 m-1,融合算法的RMSE为0.098 m-1;当Kd设置为0.3 m-1时,计算得到Fernald法的RMSE为0.975 m-1,Klett法的RMSE为0.984 m-1,融合算法的RMSE为0.957 m-1。仿真结果显示出融合算法在近岸水域反演漫射衰减系数Kd的优势。

    • 2021年10月在某水域开展了海洋激光雷达探测漫射衰减系数的实验。

      实验平台为岸基扫描平台,双频激光雷达设备以向水中发射脉冲激光,并进行小区域扫描。激光经由海面、海水、海底的反射被探测器接收。同时,采用某自由落体式高光谱辐射剖面测量系统采集,获得高精度的向下辐照度,计算Kd值(漫射衰减系数)的公式为Kd(z)=-(dEd/dz)/Ed, Ed为下行辐照度。将得到原位测量数据进行一次拟合,计算得到此海域Kd =0.32 m-1,利用原位测量得到水深为5.5 m。

      实测信号采用某水域的海试数据,其近岸浅水深度约为5.5 m,迭代初始值利用斜率法计算海底处的Kd值,考虑到海底处的脉冲展宽会使海底有一个突起的波峰,因此选用5.2 m处作为迭代初始位置。首先利用Klett算法反演出Kd值曲线,然后分别用不同激光雷达比的Fernald算法得出反演Kd的曲线,从Kd=0.317 m-1处迭代,在1.5 m深度Kd收敛到0.325 m-1左右。图 7中显示了Klett算法最为接近的水体悬浮物激光雷达比为210 sr,以此为依据作为Fernald法最后迭代的初始值。

      图  7  Klett算法和不同激光雷达比的Fernald法实测反演结果

      Figure 7.  Inversion results of Klett algorithm and Fernald method with different LiDAR ratios in experiment

      融合算法最后反演出的Kd收敛值如图 8所示。融合算法反演结果与原位测量拟合数据在1.5 m处收敛值一致,测得此近岸海域最佳激光雷达比为210 sr,海水Kd=0.317 m-1。针对同一实测波形,Klett算法反演Kd值的RMSE为0.85 m-1, Fermald算法反演Kd值的RMSE为0.8 m-1,融合算法反演Kd值的RMSE为0.74 m-1。实验结果进一步验证了所提出的融合算法的有效性,表明了此算法可以很好地反演漫射衰减系数Kd

      图  8  融合算法反演Kd值实测结果

      Figure 8.  Diagram of inversion Kd value of fusion algorithm in experiment

    • 针对近岸浅海水域,提出了一种基于Fernald法和Klett法的融合算法,该融合算法以Klett法反演结果为基准,将计算出的激光雷达比作为Fernald法迭代时的先验信息,最后用Fernald法反演出更加精确的水体漫射衰减系数Kd。利用解析模型仿真了激光雷达回波信号,通过设置与海试环境一致的环境参数和系统参数, 验证了融合算法的可行性,将融合算法应用于仿真波形数据和实测波形数据,通过与原位测量数据比较,利用融合算法反演得到Kd值的RMSE为0.74 m-1,精度优于传统的Klett法和Fernald法。实验结果表明,所提出的融合算法在近岸水体条件下能够快速找到水体悬浮物激光雷达比, 并能准确反演出Kd值。

参考文献 (25)

目录

    /

    返回文章
    返回