CN107976686A - 一种多视场角海洋激光雷达及其视场角择优方法 - Google Patents
一种多视场角海洋激光雷达及其视场角择优方法 Download PDFInfo
- Publication number
- CN107976686A CN107976686A CN201711013714.6A CN201711013714A CN107976686A CN 107976686 A CN107976686 A CN 107976686A CN 201711013714 A CN201711013714 A CN 201711013714A CN 107976686 A CN107976686 A CN 107976686A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- field angle
- mfrac
- water body
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种多视场角海洋激光雷达及其视场角择优方法。本发明采用工作在蓝绿光波段的多视场角海洋激光雷达，建立参数灵敏度评价方法，择优选取多视场角激光雷达系统的视场角，在水体等强散射体的多次散射效应下，仍能够高精度地遥测水体参数的垂直分布特性。克服了传统的单视场角激光雷达在海洋水体的多次散射效应下模型失配的不足，通过建立参数灵敏度评价方法，对多视场角激光雷达系统的视场角选取进行择优，可获得最佳的水体参数反演精度。
Description
技术领域
本发明属于激光雷达技术领域，特别是涉及一种多视场角海洋激光雷达及其视场角择优方法。
背景技术
海洋覆盖了地球表面71％的面积，但人类已探索的海洋仅占海洋面积的5％。海洋提供了丰富的生物、矿产、空间资源和海洋能源，同时在全球气候调节、地球生物化学循环等方面发挥着重要作用；作为海洋食物链的起点，包括藻类等浮游植物在内的海洋植物主要生活在海水表层，通过光合作用产生供所有陆地生物呼吸的一半氧气，同时海洋吸收了陆地二氧化碳排放总量的四分之一。随着经济的日益发展，陆源污染物排放、局部海域富营养化、过度捕捞等海洋环境问题日益突出；全球气候变暖等气候问题亦会造成海平面上升、海洋生态系统失衡等后果。因此，亟需发展高效的、多样化的海洋探测手段，帮助人类更好地认识海洋、开发海洋，同时保护海洋环境。
原位测量是目前最主要的探测方式之一，通过衰减计、散射计等设备直接接触探测水体，具有较高的探测精度，但探测效率较低。被动遥感技术能够提供全球尺度上的水色观测，但却受限于探测深度和光照条件。声学探测具有较高的探测深度和探测精度，但受限于水气界面的高反射率，无法工作于飞机、卫星等平台。海洋激光雷达是一种主动式光学遥感设备，向海面发射的激光能够穿透海气界面，获得上层水体参数的垂直分布特性，具有高时空分辨率、全天时工作、可覆盖全球的优点。海洋激光雷达是一种探测海洋上层水体十分有效的手段，可应用于全球海洋初级生产力估算、浅海测深与水下地貌绘制、潜艇/鱼雷等军事目标探测、海底底质分类等领域。
海洋激光雷达通过向海面发射的激光并接收海水的后向散射信号，可从回波信号中提取出海洋的水体参数信息。激光在水中传输时具有较强的多次散射效应，造成了激光光束在时间、空间上的扩展，增加了回波信号的复杂性，传统激光雷达采用的单次散射模型难以精确地反映实际情况，不再适用于高精度的海洋激光雷达反演过程。已有的理论研究表明，较小的海洋激光雷达接收光路视场角有利于反演获得光束衰减系数c，较大的海洋激光雷达接收光路视场角有利于反演漫衰减系数Kd，通过可变视场角激光雷达可研究多次散射贡献量。本发明提出了可校正多次散射效应的多视场角激光雷达，并通过建立参数灵敏度评价方法，择优选取视场角数量和大小，对于提高海洋水体参数的遥测精度至关重要。
发明内容
本发明的目的是针对现有技术的不足，提供一种多视场角海洋激光雷达及其视场角择优方法。
本发明解决其技术问题所采用的技术方案如下：
一种多视场角海洋激光雷达，包括激光发射系统、多视场信号接收系统、数据采集与处理系统；激光发射系统中发射光路的光轴与多视场信号接收系统中接收光路的光轴相互平行且相邻，数据采集与处理系统置于后方并通过同轴电缆线与多视场信号接收系统相连；激光发射系统用于生成准直激光束，多视场信号接收系统用于接收不同视场角范围内的激光雷达回波信号并将其转换为电信号，数据采集与处理系统用于对多视场信号接收系统中生成的电信号进行数据采集与处理；
激光发射系统包括脉冲激光器和扩束镜；脉冲激光器和扩束镜与水平面垂直设置，脉冲激光器的光发射口与扩束镜的入光口相对设置，扩束镜的出光口与水平面平行设置；
多视场信号接收系统包括三个结构相同的具有不同视场角的信号接收通道；小视场信号接收通道包括小视场望远镜、第一光阑、第一窄带滤光片、第一会聚透镜、第一光电探测器；中视场信号接收通道包括中视场望远镜、第二光阑、第二窄带滤光片、第二会聚透镜、第二光电探测器；大视场信号接收通道包括大视场望远镜、第三光阑、第三窄带滤光片、第三会聚透镜、第三光电探测器；小视场信号接收通道、中视场信号接收通道、大视场信号接收通道垂直于水平面与脉冲激光器、扩束镜平行设置；各信号接收通道中的望远镜入光口与水平面平行设置，望远镜出光口与光阑中心相对设置；光阑上方设置有窄带滤光片，窄带滤光片上方设置有会聚透镜，会聚透镜上方设置有光电探测器；
数据采集与处理系统包括高速数据采集卡和计算机；高速数据采集卡与多视场信号接收系统中的三个光电探测器通过同轴电缆线连接，高速数据采集卡与计算机通过电连接。
一种多视场角海洋激光雷达的视场角择优方法，具体包括如下：
该视场角择优方法根据拟反演的水体光学参数确定多视场角海洋激光雷达系统的视场角的数量和大小；进入海水中的激光，一部分被吸收，改变光的传输能量，另一部分被散射，改变光的传播方向，吸收系数a和散射系数b描述了上述海水对激光的吸收与散射效应，而前向散射展宽系数m描述了在海水多次散射效应下光场重新分布的特征，以上三个参数代表了最为重要的水体光学参数；通过建立激光雷达回波方程，并在拟定的水体分层方法下，分析这三个参数的灵敏度特征，基于经验值与灵敏度特性分析，分别确定反演各个参数最佳的视场角；最终确定接收激光雷达回波信号采用信号接收通道方案。
所述的视场角的数量和大小具体通过如下步骤获取：
步骤1.设接收光路视场角为γr时，深度为z处的水体产生的激光雷达回波信号可表示为式(1)：
其中，C0为系统常数，n为折射率、H为激光雷达系统距水面的高度，τ(z)为单次散射条件下的光学厚度，为探测器有效区域的傅里叶变换，v为空间频率，ξ代表水深的积分变量，P(z,p)为前向散射相函数的Hankel变换，采用俄罗斯科学家Dolin提出的散射相函数模型，即为式(2)：
其中，m为前向散射展宽系数，p为散射角的角频率；
步骤2.将水体以Δzi的厚度分成N层，其中i＝1,2…N；每层水体内部的参数均匀分布，即吸收系数a、散射系数b、前向散射展宽系数m在每层水体内为定值ai、bi、mi；假设水体具有等厚分层的特征，则第i层水体的厚度可表示为式(3)：
其中，Δz为等厚分层水体的层厚，zN＝z表示深度为z处的水体；则步骤1所述的激光雷达回波信号方程可简化为式(4)：
其中，B(v)是与前向散射展宽系数mi有关的项，可表示为式(5)：
其中z0＝0表示水面位置，zi-1和zi分别表示第i层水体的起始深度和终止深度，arsinh是反双曲正弦函数；
步骤3.定义由深度为z处的单位水体产生的激光雷达回波信号恢复出第i层水体参数的灵敏度为水体参数X包括吸收系数a、散射系数b、前向散射展宽系数m，则有吸收系数灵敏度为式(6)：
散射系数灵敏度为式(7)：
前向散射展宽系数灵敏度为式(8)：
步骤4.将激光雷达回波信号方程式(4)代入式(6)、(7)、(8)并化简可得到式(9)：
其中：
由式(9-a)可知，吸收系数灵敏度仅与吸收系数a、水体层次厚度Δzi呈简单的线性关系，与视场角大小无关；不同的水体参数对视场角变化具有不同的响应特性；
步骤5.通过设置三个具有不同视场角的接收通道，获得吸收系数a、散射系数b、前向散射展宽系数m三个水体参数的最佳反演。
所述的三个具有不同视场角的接收通道，其大小分别是200mrad、14mrad和10mrad。
脉冲激光器工作在蓝绿光波段，脉冲宽度不大于10ns，单脉冲能量不小于5mJ。扩束镜采用抗强激光扩束镜。窄带滤光片采用带宽不大于10nm的窄带滤光片。高速数据采集卡采样率不低于400MSPS，量化位数不低于12位。
本发明有益效果如下：
本发明采用工作在蓝绿光波段的多视场角海洋激光雷达，建立参数灵敏度评价方法，择优选取多视场角激光雷达系统的视场角，在水体等强散射体的多次散射效应下，仍能够高精度地遥测水体参数的垂直分布特性。克服了传统的单视场角激光雷达在海洋水体的多次散射效应下模型失配的不足，通过建立参数灵敏度评价方法，对多视场角激光雷达系统的视场角选取进行择优，可获得最佳的水体参数反演精度。
附图说明
图1为本发明的整体结构模块图；
图2为本发明实施例中多视场角海洋激光雷达结构示意图；
图3为本发明视场角择优方法流程图；
图4本发明实施例中水体分层示意图；
图5本发明实施例中散射系数和前向散射展宽系数的灵敏度情况。
具体实施方式
下面结合实施例和附图对本发明作进一步说明。
如图1和2所示，一种多视场角海洋激光雷达，包括激光发射系统1、多视场信号接收系统2、数据采集与处理系统3。激光发射系统1中发射光路的光轴与多视场信号接收系统2中接收光路的光轴相互平行且相邻，数据采集与处理系统3置于后方并通过同轴电缆线与多视场信号接收系统2相连。激光发射系统1用于生成准直激光束，多视场信号接收系统2用于接收不同视场角范围内的激光雷达回波信号并将其转换为电信号，数据采集与处理系统3用于对多视场信号接收系统2中生成的电信号进行数据采集与处理。
激光发射系统1包括脉冲激光器1-1和扩束镜1-2；脉冲激光器1-1和扩束镜1-2与水平面垂直设置，脉冲激光器1-1的光发射口与扩束镜1-2的入光口相对设置，扩束镜1-2的出光口与水平面平行设置。
多视场信号接收系统2包括三个结构相同的具有不同视场角的信号接收通道；小视场信号接收通道2-1包括小视场望远镜2-1-1、第一光阑2-1-2、第一窄带滤光片2-1-3、第一会聚透镜2-1-4、第一光电探测器2-1-5；中视场信号接收通道2-2包括中视场望远镜2-2-1、第二光阑2-2-1、第二窄带滤光片2-2-3、第二会聚透镜2-2-4、第二光电探测器2-2-5；大视场信号接收通道2-3包括大视场望远镜2-3-1、第三光阑2-3-2、第三窄带滤光片2-3-3、第三会聚透镜2-3-4、第三光电探测器2-3-5；小视场信号接收通道2-1、中视场信号接收通道2-2、大视场信号接收通道2-3垂直于水平面与脉冲激光器1-1、扩束镜1-2平行设置。各信号接收通道中的望远镜入光口与水平面平行设置，望远镜出光口与光阑中心相对设置；光阑上方设置有窄带滤光片，窄带滤光片上方设置有会聚透镜，会聚透镜上方设置有光电探测器。
数据采集与处理系统3包括高速数据采集卡3-1和计算机3-2。高速数据采集卡3-1与多视场信号接收系统2中的三个光电探测器通过同轴电缆线连接，高速数据采集卡3-1与计算机3-2通过电连接。
基于上述多视场角海洋激光雷达，该多视场角海洋激光雷达的视场角择优方法，具体如下：
视场角择优方法根据拟反演的水体光学参数确定多视场角海洋激光雷达系统的视场角的数量和大小；如图3所示，进入海水中的激光，一部分被吸收，改变光的传输能量，另一部分被散射，改变光的传播方向，吸收系数a和散射系数b描述了上述海水对激光的吸收与散射效应，而前向散射展宽系数m描述了在海水多次散射效应下光场重新分布的特征，以上三个参数代表了最为重要的水体光学参数；通过建立激光雷达回波方程，并在拟定的水体分层方法下，分析这三个参数的灵敏度特征，基于经验值与灵敏度特性分析，分别确定反演各个参数最佳的视场角；最终确定接收激光雷达回波信号采用信号接收通道方案。通过以下具体步骤得到视场角的数量和大小：
步骤1.设接收光路视场角为γr时，深度为z处的水体产生的激光雷达回波信号可表示为式(1)：
其中，C0为系统常数，n为折射率、H为激光雷达系统距水面的高度，τ(z)为单次散射条件下的光学厚度，为探测器有效区域的傅里叶变换，v为空间频率，ξ代表水深的积分变量，P(z,p)为前向散射相函数的Hankel变换，采用俄罗斯科学家Dolin提出的散射相函数模型，即为式(2)：
其中，m为前向散射展宽系数，p为散射角的角频率。
步骤2.如图4所示，将水体以Δzi的厚度分成N层，其中i＝1,2…N；每层水体内部的参数均匀分布，即吸收系数a、散射系数b、前向散射展宽系数m在每层水体内为定值ai、bi、mi；假设水体具有等厚分层的特征，则第i层水体的厚度可表示为式(3)：
其中，Δz为等厚分层水体的层厚，zN＝z表示深度为z处的水体；则步骤1所述的激光雷达回波信号方程可简化为式(4)：
其中，B(v)是与前向散射展宽系数mi有关的项，可表示为式(5)：
其中z0＝0表示水面位置，zi-1和zi分别表示第i层水体的起始深度和终止深度，arsinh是反双曲正弦函数；
步骤3.定义由深度为z处的单位水体产生的激光雷达回波信号恢复出第i层水体参数的灵敏度为实施例选取的水体参数X包括吸收系数a、散射系数b、前向散射展宽系数m，则有吸收系数灵敏度为式(6)：
散射系数灵敏度为式(7)：
前向散射展宽系数灵敏度为式(8)：
步骤4.将激光雷达回波信号方程式(4)代入式(6)、(7)、(8)并化简可得到式(9)：
其中：
由式(9-a)可知，吸收系数灵敏度仅与吸收系数a、水体层次厚度Δzi呈简单的线性关系，与视场角大小无关；根据式(9-b)、(9-c)，在由水深z＝6m处的信号反演层次厚度Δzi＝2m的第3层次的条件下进行仿真，得到图5，分别表明了散射系数灵敏度和前向散射展宽系数灵敏度与视场角大小的关系；由图5(a)、图5(b)可知，散射系数灵敏度前向散射展宽系数灵敏度呈现出小视场条件下变化缓慢、大视场条件下快速下降的特性，前者在10mrad处出现快速下降转折点，后者在14mrad处出现了快速下降转折点。可见，不同的水体参数对视场角变化具有不同的响应特性。
步骤5.吸收系数a、散射系数b、前向散射展宽系数m三个水体参数是影响激光在水中传输特性的主要因素，对于建立激光雷达模型至关重要，但这三个参数无法由单一视场角的激光雷达回波信号同时实现准确反演；引入多视场角激光雷达系统，在不同视场角得到不同参数的最佳反演；根据步骤4中所述参数灵敏度的特征，将大视场信号接收通道的视场角设置为200mrad，此时散射系数灵敏度前向散射展宽系数灵敏度均处于5‰量级，对激光雷达回波信号影响较小，同时吸收系数灵敏度为定值，在该视场角下可进行吸收系数a的最佳反演；根据图5(a)、图5(b)中信号的转折点，同时兼顾激光雷达回波信号的信噪比，将中视场信号接收通道、小视场信号接收通道的视场角分别设置为14mrad、10mrad，并将上述已反演得到的吸收系数a代入，即可分别得到前向散射展宽系数m、散射系数b的最佳反演。
由以上步骤可知，为获得吸收系数a、散射系数b、前向散射展宽系数m三个水体参数的最佳反演，应设置三个具有不同视场角的接收通道，其大小分别是200mrad、14mrad和10mrad。
脉冲激光器工作在蓝绿光波段，脉冲宽度不大于10ns，单脉冲能量不小于5mJ。本实施例中脉冲激光器1-1采用脉宽宽度不大于3ns的超窄脉宽脉冲激光器，例如奥地利Montfort Laser公司的半导体泵浦固体激光器，中心波长为531.7nm，脉冲宽度3ns，单脉冲能量10mJ，重复频率10Hz；
扩束镜1-2采用抗强激光扩束镜，例如美国Thorlabs公司的BE03-532型号扩束镜，3倍扩束；
望远镜、光阑、会聚透镜根据所需的视场角大小分别确定，选配合适的参数以实现不同视场角下的最优设计；
窄带滤光片采用带宽不大于10nm的窄带滤光片，例如武汉优光公司生产的带通滤光片，带宽小于10nm，透过率大于70％；
光电探测器必须有高响应速度和高灵敏度，例如日本滨松公司的H10721-110型号光电倍增管，上升时间0.5ns；
高速数据采集卡3-1采样率不低于400MSPS，量化位数不低于12位，例如德国Spectrum公司的M4i-4451型号数据采集卡。
本发明的工作过程如下：
脉冲激光器发射的窄脉宽激光束经扩束镜后得到准直激光束，发散角减小，进入到被探测水体中；窄脉宽激光束在水体中传输时受到水分子、颗粒物等强散射介质的多次散射，多次散射信号中的后向散射部分形成了激光雷达回波信号。信号接收系统接收上述激光雷达回波信号；三路具有不同大小视场角的接收望远镜收集所设视场角范围内的激光雷达回波信号；通过调节三个光阑实现所需的视场角，其大小与分布由视场角择优方法确定；三个窄带滤光片用于滤除回波信号中混杂的以太阳光为主的背景辐射；三个会聚透镜将光信号会聚于后方三个光电探测器的接受面上；三个光电探测器分别将三路通道的光信号转换为电信号，并由高速数据采集卡数字化后传输至计算机；计算机将采集到的激光雷达回波信号作进一步的反演处理，从而得到被探测水体的参数。
Claims (8)
1.一种多视场角海洋激光雷达，其特征在于包括激光发射系统(1)、多视场信号接收系统(2)、数据采集与处理系统(3)；激光发射系统(1)中发射光路的光轴与多视场信号接收系统(2)中接收光路的光轴相互平行且相邻，数据采集与处理系统(3)置于后方并通过同轴电缆线与多视场信号接收系统(2)相连；激光发射系统(1)用于生成准直激光束，多视场信号接收系统(2)用于接收不同视场角范围内的激光雷达回波信号并将其转换为电信号，数据采集与处理系统(3)用于对多视场信号接收系统(2)中生成的电信号进行数据采集与处理；
激光发射系统(1)包括脉冲激光器(1-1)和扩束镜(1-2)；脉冲激光器(1-1)和扩束镜(1-2)与水平面垂直设置，脉冲激光器(1-1)的光发射口与扩束镜(1-2)的入光口相对设置，扩束镜(1-2)的出光口与水平面平行设置；
多视场信号接收系统(2)包括三个结构相同的具有不同视场角的信号接收通道；小视场信号接收通道(2-1)包括小视场望远镜(2-1-1)、第一光阑(2-1-2)、第一窄带滤光片(2-1-3)、第一会聚透镜(2-1-4)、第一光电探测器(2-1-5)；中视场信号接收通道(2-2)包括中视场望远镜(2-2-1)、第二光阑(2-2-1)、第二窄带滤光片(2-2-3)、第二会聚透镜(2-2-4)、第二光电探测器(2-2-5)；大视场信号接收通道(2-3)包括大视场望远镜(2-3-1)、第三光阑(2-3-2)、第三窄带滤光片(2-3-3)、第三会聚透镜(2-3-4)、第三光电探测器(2-3-5)；小视场信号接收通道(2-1)、中视场信号接收通道(2-2)、大视场信号接收通道(2-3)垂直于水平面与脉冲激光器(1-1)、扩束镜(1-2)平行设置；各信号接收通道中的望远镜入光口与水平面平行设置，望远镜出光口与光阑中心相对设置；光阑上方设置有窄带滤光片，窄带滤光片上方设置有会聚透镜，会聚透镜上方设置有光电探测器；
数据采集与处理系统(3)包括高速数据采集卡(3-1)和计算机(3-2)；高速数据采集卡(3-1)与多视场信号接收系统(2)中的三个光电探测器通过同轴电缆线连接，高速数据采集卡(3-1)与计算机(3-2)通过电连接。
2.根据权利要求1所述的多视场角激光雷达的视场角择优方法，其特征在于具体包括如下：
该视场角择优方法根据拟反演的水体光学参数确定多视场角激光雷达系统的视场角的数量和大小；进入海水中的激光，一部分被吸收，改变光的传输能量，另一部分被散射，改变光的传播方向，吸收系数a和散射系数b描述了上述海水对激光的吸收与散射效应，而前向散射展宽系数m描述了在海水多次散射效应下光场重新分布的特征，以上三个参数代表了最为重要的水体光学参数；通过建立激光雷达回波方程，并在拟定的水体分层方法下，分析这三个参数的灵敏度特征，基于经验值与灵敏度特性分析，分别确定反演各个参数最佳的视场角；最终确定接收激光雷达回波信号采用信号接收通道方案。
3.根据权利要求2所述的多视场角激光雷达的视场角择优方法，其特征在于视场角的数量和大小具体通过如下步骤获取：
步骤1.设接收光路视场角为γr时，深度为z处的水体产生的激光雷达回波信号可表示为式(1)：
其中，C0为系统常数，n为折射率、H为激光雷达系统距水面的高度，τ(z)为单次散射条件下的光学厚度，为探测器有效区域的傅里叶变换，v为空间频率，ξ代表水深的积分变量，P(z,p)为前向散射相函数的Hankel变换，采用俄罗斯科学家Dolin提出的散射相函数模型，即为式(2)：
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>m</mi>
<msqrt>
<mrow>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>p</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，m为前向散射展宽系数，p为散射角的角频率；
步骤2.将水体以Δzi的厚度分成N层，其中i＝1,2…N；每层水体内部的参数均匀分布，即吸收系数a、散射系数b、前向散射展宽系数m在每层水体内为定值ai、bi、mi；假设水体具有等厚分层的特征，则第i层水体的厚度可表示为式(3)：
<mrow>
<msub>
<mi>&Delta;z</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
<mo>,</mo>
<mi>i</mi>
<mo><</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>z</mi>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，Δz为等厚分层水体的层厚，zN＝z表示深度为z处的水体，；则步骤1所述的激光雷达回波信号方程可简化为式(4)：
其中，B(v)是与前向散射展宽系数mi有关的项，可表示为式(5)：
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>v</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mi>r</mi>
<mi>sinh</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<mi>a</mi>
<mi>r</mi>
<mi>sinh</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中z0＝0表示水面位置，zi-1和zi分别表示第i层水体的起始深度和终止深度，arsinh()是反双曲正弦函数；
步骤3.定义由深度为z处的单位水体产生的激光雷达回波信号恢复出第i层水体参数的灵敏度为选取的水体参数X包括吸收系数a、散射系数b、前向散射展宽系数m，则有吸收系数灵敏度为式(6)：
<mrow>
<msubsup>
<mi>C</mi>
<mi>i</mi>
<mi>a</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
散射系数灵敏度为式(7)：
<mrow>
<msubsup>
<mi>C</mi>
<mi>i</mi>
<mi>b</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
前向散射展宽系数灵敏度为式(8)：
<mrow>
<msubsup>
<mi>C</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤4.将激光雷达回波信号方程式(4)代入式(6)、(7)、(8)并化简可得到式(9)：
其中：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>v</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mi>r</mi>
<mi>sinh</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<mi>a</mi>
<mi>r</mi>
<mi>sinh</mi>
<mo>(</mo>
<mfrac>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>m</mi>
<mi>i</mi>
</msub>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>m</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mi>m</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
由式(9-a)可知，吸收系数灵敏度仅与吸收系数a、水体层次厚度Δzi呈简单的线性关系，与视场角大小无关；不同的水体参数对视场角变化具有不同的响应特性；
步骤5.通过设置三个具有不同视场角的接收通道，获得吸收系数a、散射系数b、前向散射展宽系数m三个水体参数的最佳反演。
4.根据权利要求3所述的多视场角激光雷达的视场角择优方法，其特征在于所述的三个具有不同视场角的接收通道，其大小分别是200mrad、14mrad和10mrad。
5.根据权利要求1所述的多视场角激光雷达的视场角择优方法，其特征在于脉冲激光器工作在蓝绿光波段，脉冲宽度不大于10ns，单脉冲能量不小于5mJ。
6.根据权利要求1所述的多视场角激光雷达的视场角择优方法，其特征在于扩束镜(1-2)采用抗强激光扩束镜。
7.根据权利要求1所述的多视场角激光雷达的视场角择优方法，其特征在于窄带滤光片采用带宽不大于10nm的窄带滤光片。
8.根据权利要求1所述的多视场角激光雷达的视场角择优方法，其特征在于高速数据采集卡(3-1)采样率不低于400MSPS，量化位数不低于12位。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711013714.6A CN107976686B (zh) | 2017-10-26 | 2017-10-26 | 一种多视场角海洋激光雷达及其视场角择优方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711013714.6A CN107976686B (zh) | 2017-10-26 | 2017-10-26 | 一种多视场角海洋激光雷达及其视场角择优方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107976686A true CN107976686A (zh) | 2018-05-01 |
| CN107976686B CN107976686B (zh) | 2020-12-08 |
Family
ID=62012793
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711013714.6A Active CN107976686B (zh) | 2017-10-26 | 2017-10-26 | 一种多视场角海洋激光雷达及其视场角择优方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107976686B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110865389A (zh) * | 2019-10-29 | 2020-03-06 | 浙江大学 | 一种海洋激光雷达系统响应优化处理方法 |
| CN113655495A (zh) * | 2021-08-31 | 2021-11-16 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
| CN113740873A (zh) * | 2021-08-31 | 2021-12-03 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
| CN114325656A (zh) * | 2021-12-30 | 2022-04-12 | 浙江大学 | 一种探测水体生物光学特性剖面的激光雷达及方法 |
| CN115266509A (zh) * | 2022-09-26 | 2022-11-01 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于激光雷达的水下垂向悬浮物浓度探测方法及系统 |
| CN118294928A (zh) * | 2024-06-06 | 2024-07-05 | 深圳市速腾聚创科技有限公司 | 激光雷达及可移动设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101071171A (zh) * | 2007-06-06 | 2007-11-14 | 中国科学院安徽光学精密机械研究所 | 双波长双视场米散射激光雷达的结构及其探测方法 |
| CN102854514A (zh) * | 2012-09-21 | 2013-01-02 | 武汉大学 | 一种近场同轴双视场米散射大气激光雷达 |
| CN204649962U (zh) * | 2015-05-19 | 2015-09-16 | 浙江大学 | 一种基于望远镜阵列的大气遥感激光雷达光学接收装置 |
-
2017
- 2017-10-26 CN CN201711013714.6A patent/CN107976686B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101071171A (zh) * | 2007-06-06 | 2007-11-14 | 中国科学院安徽光学精密机械研究所 | 双波长双视场米散射激光雷达的结构及其探测方法 |
| CN102854514A (zh) * | 2012-09-21 | 2013-01-02 | 武汉大学 | 一种近场同轴双视场米散射大气激光雷达 |
| CN204649962U (zh) * | 2015-05-19 | 2015-09-16 | 浙江大学 | 一种基于望远镜阵列的大气遥感激光雷达光学接收装置 |
Non-Patent Citations (2)
| Title |
|---|
| 李俊 等: "探测武汉上空大气气溶胶的双视场激光雷达", 《光学学报》 * |
| 李志刚 等: "一种激光雷达测量海洋光学参数的新方法", 《应用光学》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110865389A (zh) * | 2019-10-29 | 2020-03-06 | 浙江大学 | 一种海洋激光雷达系统响应优化处理方法 |
| CN113655495A (zh) * | 2021-08-31 | 2021-11-16 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
| CN113740873A (zh) * | 2021-08-31 | 2021-12-03 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
| CN113740873B (zh) * | 2021-08-31 | 2024-05-28 | 自然资源部第二海洋研究所 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
| CN113655495B (zh) * | 2021-08-31 | 2024-06-11 | 自然资源部第二海洋研究所 | 一种基于水体漫衰减系数光谱依赖性的星载海洋激光雷达最优波段评估方法 |
| CN114325656A (zh) * | 2021-12-30 | 2022-04-12 | 浙江大学 | 一种探测水体生物光学特性剖面的激光雷达及方法 |
| CN114325656B (zh) * | 2021-12-30 | 2024-05-24 | 浙江大学 | 一种探测水体生物光学特性剖面的激光雷达及方法 |
| CN115266509A (zh) * | 2022-09-26 | 2022-11-01 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于激光雷达的水下垂向悬浮物浓度探测方法及系统 |
| CN115266509B (zh) * | 2022-09-26 | 2023-02-24 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于激光雷达的水下垂向悬浮物浓度探测方法及系统 |
| CN118294928A (zh) * | 2024-06-06 | 2024-07-05 | 深圳市速腾聚创科技有限公司 | 激光雷达及可移动设备 |
| CN118294928B (zh) * | 2024-06-06 | 2024-09-10 | 深圳市速腾聚创科技有限公司 | 激光雷达及可移动设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107976686B (zh) | 2020-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107976686A (zh) | 一种多视场角海洋激光雷达及其视场角择优方法 | |
| CN109298410B (zh) | 一种海洋溢油探测激光雷达 | |
| CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
| Guenther et al. | Meeting the accuracy challenge in airborne lidar bathymetry | |
| Zhou et al. | An optical system for suppression of laser echo energy from the water surface on single-band bathymetric LiDAR | |
| Tulldahl et al. | Classification of aquatic macrovegetation and substrates with airborne lidar | |
| CN106019306A (zh) | 一种基于计算鬼成像原理实现水下目标探测的装置 | |
| CN109061668B (zh) | 一种多视场偏振激光雷达探测系统及探测冰云的方法 | |
| CN105372641A (zh) | 激光雷达水体回波信号仿真方法 | |
| Steinvall et al. | Experimental evaluation of an airborne depth-sounding lidar | |
| Tulldahl et al. | Simulation of sea surface wave influence on small target detection with airborne laser depth sounding | |
| CN115797760B (zh) | 主被动融合的水质立体遥感反演方法、系统及存储介质 | |
| Chen et al. | Lidar remote sensing for vertical distribution of seawater optical properties and chlorophyll-a from the East China Sea to the South China Sea | |
| CN115266509B (zh) | 一种基于激光雷达的水下垂向悬浮物浓度探测方法及系统 | |
| Shrestha et al. | Shallow bathymetric mapping via multistop single photoelectron sensitivity laser ranging | |
| Xu et al. | Design and validation of a shipborne multiple-field-of-view lidar for upper ocean remote sensing | |
| Lv et al. | Analysis of wave fluctuation on underwater acoustic communication based USV | |
| CN101464522B (zh) | 基于激光波前变化的水中动态目标尾迹探测系统及方法 | |
| CN109632258A (zh) | 一种基于矢量传感器的收发分离的海洋内波声学检测方法 | |
| Bailly et al. | Airborne LiDAR methods applied to riverine environments | |
| CN101533105B (zh) | 基于激光光斑漂移的水中动态目标尾迹探测系统及方法 | |
| Yu et al. | High-precision 3D imaging of underwater coaxial scanning photon counting Lidar based on spatiotemporal correlation | |
| Guenther et al. | Laser applications for near-shore nautical charting | |
| Dzieciuch et al. | Propagation of sound through a spicy ocean, the SOFAR overture | |
| CN111780726A (zh) | 海洋激光测深方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liu Dong Inventor after: Xu Peituo Inventor after: Zhou Yudi Inventor before: Liu Dong Inventor before: Xu Peituo |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |