一种双波长海洋激光雷达探测的波长择优方法与流程

本发明属于激光雷达技术领域，尤其是涉及一种双波长海洋激光雷达探测叶绿素和cdom吸收系数的波长择优方法。

背景技术：

海洋覆盖了地球的71％以上，在全球生态系统和气候变化中扮演着非常重要的角色。所以对海洋的探测显得尤为重要。目前有以下几种方法来检测海洋内部。原位方法如透射仪和散射仪等，可以准确地获取海洋信息，但效率受到限制。基于飞机或卫星的海洋彩色遥感(例如“海洋宽视野传感器”)能够长期有效地收集全球数据。然而，垂直分布状况探测能力不足和对自然光的依赖性限制了其应用。声学被广泛用于对海水进行剖析，但是由于空气-水界面的大量损失，它们只能在水下平台上工作。激光雷达受平台和自然光的限制较少，是一种对海水进行剖面探测的有效方法，具备高时空分辨率、昼夜连续观测、全球尺度测量等优势。到目前为止，海洋激光雷达已能够成功应用于生物探测、浮游植物层分布研究、浅海地形测绘、海洋内波探测等相关研究。

目前，激光雷达中广泛使用单波长来探测海洋。此外，几乎所有的工作波长都是532nm，这是nd：yag激光器的二次谐波波长。例如，正交偏振云气溶胶激光雷达是一种简单的弹性后向散射激光雷达，发射波长为532和1064nm，但是对于1064nm，其在水中的衰减太大，无法提供有用的水下信息。532nm的高光谱分辨率激光雷达(hsrl)是一种更先进的激光雷达，它可以独立分离衰减和反向散射。当海水为仅包含水分子和浮游植物的i类水时，单波长激光雷达可以使用现有的生物光学模型获得叶绿素a。然而，如果海水中还存在诸如有色溶解有机物(cdom)之类的复杂物质，则仅通过使用一个波长就无法准确获得叶绿素a，而多波长激光雷达是解决这一问题的方法之一。

gray等人开发了一种多波长激光雷达，证明了激光雷达波长与水类型之间的关系，但激光雷达仍处于实验室使用阶段，并未用于检测叶绿素a。hoge提出了双波长激光雷达的想法来解决这个问题，但是他没有提供关于波长选择的建议。hostetler等提议使用355nm和532nm的波长，但是，选择355nm和532nm是因为它们是nd：yag激光器的二次谐波和三次谐波波长，但并未证明它们是最佳波长。

技术实现要素：

本发明提出了一种双波长海洋激光雷达探测的波长择优方法，对双波长激光雷达的两个波长进行择优，可获得最佳的叶绿素和cdom吸收系数反演精度。

一种双波长海洋激光雷达探测的波长择优方法，包括以下步骤：

(1)基于双波长海洋高光谱分辨率激光雷达进行探测，利用激光雷达参量参数中的180度体积散射函数和激光雷达衰减系数，得到双波长下的总吸收系数a(λ1)和a(λ2)；

(2)利用反演算法得到叶绿素的反演参数模型c和有色溶解有机物cdom的反演参数模型ag(λ1)，其中，c是叶绿素浓度，ag(λ1)是cdom在第一波长下的吸收系数；

(3)根据相对误差传递公式和反演参数模型得出叶绿素的相对误差模型εc和有色溶解有机物cdom的相对误差模型

(4)根据波长和误差的关系得到双波长任意时的择优方案；再通过改变参数设置得到激光雷达的第一波长λ1取定为532nm时，第二波长λ2对反演误差的影响，从而得到单波长固定时的择优方案。

本发明中，所述的双波长海洋高光谱分辨率激光雷达包括激光发射系统、接收系统、锁频系统、滤波系统、数据采集及处理系统；其中，激光发射系统中使用两种不同波长的激光进行探测，根据不同激光波长的回波信号进行后续的波长择优计算。

步骤(1)中，双波长海洋高光谱分辨率激光雷达同时测量两个波长的水体总吸收系数，水体总吸收系数通过激光雷达衰减系数和后向散射系数得出，其中，后向散射系数通过转换因子从180度体积散射函数中获取；双波长下的总吸收系数a(λ1)和a(λ2)计算表达式为：

式中，aw(λ)是纯水各波长吸收系数，λ1为激光雷达的第一波长，λ2为激光雷达的第一波长；a(λ)是叶绿素a的比吸收系数；c是叶绿素a的浓度；e(λ)是叶绿素a的指数吸收系数；s为光谱吸收斜率，单位为nm^-1；λ0是参考波长；ag(λ0)是cdom在参考波长下的吸收系数。

步骤(2)中，叶绿素的反演参数模型c和有色溶解有机物cdom的反演参数模型ag(λ1)通过以下隐式函数的形式表示：

步骤(3)中，叶绿素的参数相对误差模型εc计算过程为：

采用经典的误差传播定律来估计c的相对误差，假设a(λ1)、a(λ2)的相对误差分别表示为δ1a(λ1)、δ2a(λ2)，c的相对误差表示为：

其中：

综合上述公式，得到：

其中，mg＝ag(λ1)/ag(λ2)、mc＝aph(λ1)/aph(λ2)、me＝e(λ1)/e(λ2)。

有色溶解有机物cdom的参数相对误差模型计算过程为：

采用经典的误差传播定律来估计ag(λ1)的相对误差：

其中，

综合上述公式，得到：

步骤(4)中，双波长任意时的择优方案确定方法具体如下：对双波长同时进行随机后作出双波长与反演误差的关系图，在图中找到反演误差最小的点，其对应的双波长即为最优的双波长。

单波长固定时的择优方案确定方法具体如下：

取定激光雷达的第一波长λ1为532nm，在取定初始ag(λ1)、初始叶绿素浓度、光谱吸收斜率s值的情况下，作第二波长与反演误差的关系图，反演误差最小时对应的第二波长即为最优波长。

与现有技术相比，本发明具有以下有益效果：

本发明采用的是双波长的海洋激光雷达，建立参数相对误差的评价方法，择优选取海洋激光雷达的两个波长，提高叶绿素和cdom吸收系数的反演精度。克服了传统单波长海洋激光雷达探测存在cdom等复杂物质的水体时反演的不精确，以及传统双波长海洋激光雷达没有建立参数相对误差传递模型的不足，通过建立参数相对误差的评价方法，对双波长激光雷达的两个波长进行择优，可获得最佳的叶绿素和cdom吸收系数反演精度。

附图说明

图1为本发明方法的流程示意图；

图2为纯水、cdom、浮游植物在不同波长下的吸收系数占比；

图3为本发明实施例中蒙特卡洛法反演和模型反演的比对图；

图4为本发明实施例中任意波长选择对误差的影响；

图5为本发明实施例中λ1＝532时，λ2对误差的影响；

图6为本发明实施例中不同c、ag(532)、s设置时反演的相对误差；

图7为本发明实施例中测量误差对反演误差的影响；

图8为本发明实施例中黄海收集的原位数据-吸收系数图；

图9为本发明实施例中不同站点数据反演与实测结果比较图。

具体实施方式

下面结合附图和实施例对本发明做进一步详细描述，需要指出的是，以下所述实施例旨在便于对本发明的理解，而对其不起任何限定作用。

如图1所示，一种双波长海洋激光雷达探测的波长择优方法，包括以下步骤：

步骤1，基于双波长海洋高光谱分辨率激光雷达，利用激光雷达参量参数——180度体积散射函数和激光雷达衰减系数，得到双波长下的总吸收系数a(λ1)和a(λ2)；

步骤2，利用反演算法得到叶绿素和cdom的反演参数模型——c和ag(λ1)；

步骤3，根据相对误差传递公式得出叶绿素和cdom的参数相对误差模型——εc和

步骤4，根据波长和误差的关系得到双波长任意时的择优方案，再通过改变参数设置得到λ1＝532时，λ2对误差的影响，从而得到单波长固定时的择优方案。

下面将详细介绍各步骤的方法和流程。

步骤1，基于双波长海洋高光谱分辨率激光雷达，所测水体的总吸收系数a(λ)具体通过如下方法获取：

双波长海洋高光谱分辨率激光雷达主要包括激光发射系统、接收系统、锁频系统、滤波系统、数据采集及处理系统，其中激光发射系统中使用两种不同波长的激光进行探测，根据不同激光波长的回波信号进行后续的计算。

激光雷达信号包含两个海水光学特性，分别是180度体积散射函数和激光雷达衰减系数。可以使用多种方法从原始信号中检索两个属性。对于弹性激光雷达，可以使用fernald方法和摄动方法，对于hsrl，两个属性可以独立获取。通过转换因子可以从180度体积的散射函数中获取后向散射系数。

gordon指出，使用kullenberg相位函数可将激光雷达衰减系数视为扩散衰减系数，而扩散衰减系数可表示为吸收系数和反向散射系数的函数。因此，吸收系数可以从激光雷达测得的漫反射衰减系数和反向散射系数中得出。此外，如果海水散射相位函数为两项henyey-greenstein形式，则激光雷达衰减系数可以直接表示为吸收系数和反向散射系数的函数。然后，可以从激光雷达的衰减系数和反向散射系数中得出吸收系数。无论如何，根据现有理论，可以使用海洋激光雷达来获取吸收系数。如果激光雷达具有双波长，则可以同时测量两个波长的吸收系数。

水体的主要成分包括水、浮游植物和cdom等。多成分的总体吸收系数是各成分吸收系数之和。因此水体总吸收系数表达式如下：

a(λ)＝aw(λ)+aph(λ)+ag(λ),(1)

其中，aw(λ)是纯水各波长吸收系数，仅与波长相关，并已被测量。aph(λ)是浮游植物吸收系数，主要由叶绿素a贡献，可由以下公式估算：

aph(λ)＝a(λ)×c^e(λ),(2)

其中，a(λ)是叶绿素a的比吸收系数，单位为m^-1；c是叶绿素a的浓度，以毫克每立方米表示。ag(λ)是cdom的吸收系数，在200-800nm范围内呈指数衰减，可用以下公式估算：

ag(λ)＝ag(λ0)*exp[-s(λ-λ0)],(3)

其中，s为光谱吸收斜率，单位为nm^-1；λ0是参考波长，ag(λ0)是cdom在参考波长下的吸收系数。

图2显示了海水各成分的吸收系数，随波长的变化趋势是不同的。纯水的吸收系数在550nm以下很小，但在550至700nm之间迅速增加。对于浮游植物的吸收光谱，在440nm处有一个主要峰，在680nm附近有一个次要峰；主峰约为次峰的峰值的两倍。另外，浮游植物在310nm处具有大的吸收峰。cdom吸收随着波长的增加呈指数下降。cdom和浮游植物的吸收趋势在360-440nm和550-680nm之间相反，但在440-550nm和310-360nm之间相似。相同的趋势不利于区分这两个组成部分。

步骤2，叶绿素和cdom的反演参数模型具体通过如下方法获取：

假设激光雷达采用双波长λ1和λ2探测，可得到下式：

叶绿素和cdom的关键参数分别为c和ag(λ1)，可以用以下隐式函数的形式表示：

cdom的浓度可以用ag(λ0)表示。因此，只需要将测量的两个吸收系数代入，并输入已知的叶绿素吸收特性和cdom吸收特性，就可以得到所测水域叶绿素和cdom的浓度。

步骤3，叶绿素和cdom的参数相对误差模型具体通过如下方法获取：

实际测量中受仪器、环境等因素影响，激光雷达测量的两个吸收系数存在误差。a(λ)误差的存在，使评估c和ag(λ1)的相对误差成为必要。相对误差越小，说明结果量的精确度受测量误差的影响越小，说明该方法的可行性越高。采用经典的误差传播定律来估计c的相对误差，假设a(λ1)、a(λ2)的相对误差分别表示为δ1a(λ1)、δ2a(λ2)，c相对误差可表示为：

其中:

综合公式(6)、(7)得到：

其中，mg＝ag(λ1)/ag(λ2)、mc＝aph(λ1)/aph(λ2)、me＝e(λ1)/e(λ2)。

同样的，ag(λ1)相对误差是：

其中，

综合公式(9)、(10)得到：

步骤4，对双波长激光雷达进行波长择优具体通过如下方法进行：首先验证模型的正确性，并从双波长任意和单波长固定两种情况提供波长择优的方案，同时利用反演模型对从黄海采集到的真实水体吸收数据进行反演，验证择优结果。

图3展示了蒙特卡洛反演和模型反演的结果，从图3中的(a)和(c)可以看到mc的大多数结果都局限于预测的3-σ误差范围，同时图3中的(b)和(d)误差理论值和mc的均方根基本吻合，这表明理论模型与mc方法之间具有很好的一致性。

图4展示了双波长任意时的误差结果。令测量误差等于0.1，同时设置相关参数为表1中的基准值，作双波长与误差的关系图。可以看到cdom反演的相对误差总体小于c的相对反演误差。当λ1>550nm或λ1<400nm时，c的相对误差都大于100％，相比之下cdom的误差可以通过改变λ2得到显著的降低。另外这种误差和物质浓度的关系很大，若模拟情况下的叶绿素浓度高于设置值，则叶绿素的误差将进一步减小。另外，叶绿素反演和cdom反演的误差有很多相似之处：当400nm<λ1<600nm和λ2<430nm时，误差相对较小；当λ1和λ2的值非常接近时，误差值很大。叶绿素误差最小在λ1＝488nm和λ2＝300nm时得到为51.01％，cdom误差最小在λ1＝500nm和λ2＝300nm时得到为14.58％。

表1

鉴于目前532nm单波长激光器发展得较为成熟，同时已有投入应用的532激光雷达，而本文的目的在于探讨在现有的激光雷达上若增加一个波长通道，波长设置在哪里更优的问题。因此在单波长反演中将波长设置为532nm；双波长反演中，将在λ1固定为532nm的基础上，探究不同λ2设置对反演误差的影响。

同样令测量误差等于0.1，设置除波长外的参数为表1中的基准值。图5展示了当λ1＝532nm时，λ2对反演误差的影响。当λ2在510～570nm之间时，出现几个误差峰值；当λ2>420nm时，误差都比较大。而在小于420nm的波段内，误差相对比较小，并且随着波长的减小，进一步减小。因此为了获得较小的反演误差，应将λ2设置在小于420nm的短波长范围内。

为了说明所得规律的普适性，通过改变初始ag(532)、初始叶绿素浓度c和s值，得到在300～420nm的λ2与误差的关系，图6展示了在不同参数设置情况下λ2和误差的关系，得到了误差随着波长减小而减小，极小值出现在358nm左右的位置的结果。

图7显示了仪器测量误差对反演结果的影响。显然，测量误差越小，反演误差就越小，并且测量误差的值不会影响λ2与误差之间的关系：反演误差通常随λ2的增加而增加，在358nm附近存在极小值。

图8展示了从黄海收集到的原位吸收数据。现场设备通过船后甲板上的收缩孔沉入海水中，以收集水柱的固有光学特性。处理完初始数据后，获得了不同深度的吸收系数，其中s1和s2代表不同的收集位置。s1采集点位于黄海的冷水区，深度很大。s2采集点位于坪山岛附近，其深度较浅。

图9所显示的是真实水体模拟的结果。结果表明，单波长反演获得的叶绿素浓度与实际值相差很大，而双波长反演获得的结果与实际值基本吻合。同时，在两种双波长反演的情况下，当λ2设置为较小的412nm值时，反演曲线与实测曲线更吻合，验证了之前的立论推导结果。同时，在双波长反演中，当λ1固定在532nm时，λ2越小，反演结果越准确。

以上所述的实施例对本发明的技术方案和有益效果进行了详细说明，应理解的是以上所述仅为本发明的具体实施例，并不用于限制本发明，凡在本发明的原则范围内所做的任何修改、补充和等同替换，均应包含在本发明的保护范围之内。