激光雷达方程的三种解

《激光雷达方程的三种解》由会员分享，可在线阅读，更多相关《激光雷达方程的三种解（7页珍藏版）》请在金锄头文库上搜索。

1、1.3.3 激光雷达的基本原理激光雷达最基本的工作原理与普通雷达类似，即由发射系统发送一个信号，与目标作用产生的返回信号被接收系统收集并处理，以获得所需信息。不同的是激光雷达的发射信号为激光束，与普通无线电雷达发送的毫米波相比，波长要短得多。无线电雷达由于波长长，无法探测小型或微粒型目标，而用于激光雷达系统的激光波长一般在微米量级，因而能用于探测极细小的微粒和分子。图 1-3 激光雷达原理Figl-3 Principle of lidar激光雷达的基本原理如图 1-3 所示，由激光器发出一束波长为 0，宽度为 tP 的脉冲,经准直扩束后垂直射入大气，光脉冲在通过大气时受到散射和衰减，其后向散射

2、光被接收面积为 Ar 的望远镜系统接收，高度 z 处的后向散射信号功率可用雷达方程表示 36：（1.1）zTAr2tczYP2p0 ， 式中， 为接收到的散射信号的波长，P 0 为发射的激光脉冲的峰值功率，Y(z)为发送器与接收器光路的几何重叠系数，c 为光速， 为大气体积后向散射系数，zT，为大气的透过率，由朗伯定律可知: （1.2）， z0d-expT式中， 为大气的消光系数。从理论上， 为大气数密度 N (z)与散射， z，截面 的乘积，即 (1.3)ddzN，如前所述，大气中与激光雷达脉冲相互作用产生的后向散射信号成分包括了大气气体分子和气溶胶粒子，由于分子尺寸小，所产生的散射光相对较

3、弱，瑞利散射截面与激发波长的四次方成反比，大气气溶胶粒子对激光的散射光为米散射。对某一激光雷达的特定波长，分子散射直接随着大气分子的浓度而发生变化，但气溶胶散射却很复杂，取决于粒径分布以及气溶胶粒子的折射率。这些气溶胶粒子随地域、时间变化明显，所以无法对其准确估算与预测。而分子散射却可以相对较准确地估计出，只要通过标准大气信息或观测点上空的大气温度、压力分布数据即可获得。因此，实际大气体积后向散射系数和大气消光系数就包含了两个部分:分子散射部分和气溶胶散射部分，即:(1.4)zzma（1.5）下标 a 表示气溶胶，m 表示大气分子。接收信号如图 1-4 所示，在 0-A 段，发射光束还没进入接

4、收视场，接收不到后向散射信号，探测到的信号强度为 0。在 A-B 段，发射光束逐渐进入视场，发射器与接收器的重叠系数逐渐增大，接收到的后向散射信号也逐渐增强，至 B 点时，重叠系数 Y(z)=1接收到的信号最强。此后尽管视场内继续充满发射光束，但由于光速在大气中传输时按规律发散，因此回波信号也将按此规律减小，如图中 B-C 段，直至遇到气溶胶密集的-2z介质。在遇到气溶胶密集的介质时，回波信号增强，在某一距离点(图中 D 点)处达到最大，随后又会逐渐衰减至系统无法探测为止。图 1-4 激光雷达接收的信号Fig.l 一 4Lidar received signal1.3. 4 激光雷达方程的解激

5、光雷达方程提供了激光回波信号与被探测物的光学性质之间的函数关系，因此可以根据激光雷达探测到的回波信号，通过求解激光雷达方程，获得有关大气性质的信息。但从前面激光雷达方程可知，除了激光雷达系统的光电探测器接收到的大气回波信号功率，以及激光器参数和系统常数以外，方程中还有四个未知量:分子和气溶胶粒子后向散射系数， 以及它们的消光系数 ， 。由于分子散射较稳定，其散射系数zma zma和消光系数 可以通过标准大气模型或观测点上空的大气温度、压力分布数据zm获得，因而雷达方程中就剩下气溶胶的两个系数 和 为未知。为了求解这个方程，za人们想出许多办法，激光雷达方程的求解方法主要有:斜率法、Klett

6、法和 Fernald 法。a.斜率法在均匀的大气层中，整个后向散射系数 和消光系数 在该范围内是常数，在zz这种情况下， 可以用简单的斜率法求得。将雷达方程(l. 1) 简化如下:z(1.6)zdeCzP02-2将 P(z)乘以 ，再取自然对数得:(1.7)zzD02lnl对 D(z)求导得:(1.8)zdzd1如果大气均匀，则 和 均为常数， , 由下式给出:0dz(1.9)dzDz21即对曲线 D(z)进行最小二乘法拟合，求出曲线的斜率即可求出均匀大气的消光系数。斜率法由于要求均匀大气作为条件，实际是很难做到的，因为 和 在某一高z度上通常并不是常数，因而这种方法很难达到很高的精度，但它适

7、合于水平探测。b. Klett 法 为了精确求解雷达方程，必须知道 与 之间的关系。Klett 提出，对于弹性散z射，假定两者之间满足如下关系:(1.10)kB这里，B 和 k 是与激光雷达波长和气溶胶的折射率以及粒径分布有关的系数。根据经验值，k 的取值范围为: ，一般情况下，取 k=1,将其代入 (1.8)式，得:3.167.0(1.11)zdzdD2取远方某一边界值高度 zc，对应的边界值为 , , Klett 求得了该方程czDccz的稳定解为:(1.12)cczdDz2(1.13)cczBz上两式中，由于分母中积分项的符号为正，随着积分数值的增加，边界值和噪声等的影响会减少，因此方程

8、的解很容易保持稳定。这种解法的优点在于只要在边界值高度附近大气透过率 较小，对于较粗略估算czT出的 也可反演出相对较精确的 。一般情况下， 必须以较精确的方法估算czz出，才能得到精确的 值，通常边界值的灵敏度会随大气的浑浊度和边界值的取值位置z改变而变化很大。c. Fernald 法当激光在大气中传输时，要受到大气分子和气溶胶粒子的共同作用，因而雷达方程中的 和 实际应包括两部分，即分子散射部分和气溶胶散射部分。zKlett 法求解得出的是大气总的消光系数，由于米散射信号的强度与辐射波长的一(1 2)次方成比例，瑞利散射信号的强度与辐射波长的一 4 次方成比例，在波长较长或有云或气溶胶浓度

9、较大的情况下，大气回波信号中米散射信号占据主要成分，而瑞利散射信号相对很弱，可以忽略，在这种只需考虑单一成分的情况下，使用 Klett 法求解最有效。然而，在波长较短，对流层中气溶胶浓度不大，或观测高层卷云和气溶胶时，瑞利散射信号的影响就不能忽略了。Fernald 法在雷达方程中将分子散射和气溶胶散射分开来考虑，即: (1.14)zzam(1.15)下标 a 表示气溶胶，m 表示大气分子，代入激光雷达方程:(1.16) zamamp dzzzAtcYPz 0200 2exp对于气溶胶引起的米散射，其消光系数与后向散射系数的比设为 ;对zSa1于分子引起的瑞利散射，其消光系数与后向散射系数之比满

10、足如下关系:(1.17) SzSm382代入雷达方程式(1. 16 )，并将方程两边同乘以 ，得到:2z(1.18) ZZmp ap dzSdzSzActzYPt001210002 2exp(1.19)dzSSActzYPzztzSpZpZm 0110000122 2exp2ex将上式两边取积分，得:(1.20) 12exp20exp20100 10 221 ZpZZ ZmdzSActzYPzdzzS再将上式两边取对数，得：(1.21)dzSActzYPt dzSzPtpZZ Zp 02lnexplnl 2eln100 010 012100 上式两边对 Z 求导，得:(1.22) z dzSP

11、SActzYPzz Zmp Z02exp2d0110 01选取边界点高度 ，假设已知 处的后向散射系数 和消光系数 ，代入上式，ccczc并设 ， ，得到:2zzD20cz(1.23)ZdzZSzSzCCmc 1210 0exp2则气溶胶的后向散射系数为:(1.24)ZdzZSzDSz zCCmcma 1210 0exp2z 根据气溶胶消光系数与后向散射系数的关系，即可求出气溶胶消光系数:(1)边界点 以下高度的气溶胶消光系数为( 后向积分): (1.25)cz cz Cmcmca Cma dzZzSDSz zdS 212102121 exp(2)边界点 zc 以上高度的气溶胶消光系数为( 前向积分): (1.26)zc zcmcmca cma dzSDSz zdS 212102121 exp若雷达信号中气溶胶信号可以忽略，上两式的边界值(分母第一项) 则由瑞利散射信号和大气密度决定。对于瑞利散射信号，S 2 是个常数;对于米散射信号， S1 与气溶胶或云的折射率、粒径分布和形状等有关，取值范围一般在 10sr100sr 之间。大气分子的消光系数可以根据美国标准大气模型确定。