激光雷达方程的三种解

1.3.3激光雷达的基本原理激光雷达最基本的工作原理与普通雷达类似，即由发射系统发送一个信号，与目标作用产生的返回信号被接收系统收集并处理，以获得所需信息。不同的是激光雷达的发射信号为激光束，与普通无线电雷达发送的毫米波相比，波长要短得多。无线电雷达由于波长长，无法探测小型或微粒型目标，而用于激光雷达系统的激光波长一般在微米量级，因而能用于探测极细小的微粒和分子。图1-3激光雷达原理Figl-3Principleoflidar激光雷达的基本原理如图1-3所示，由激光器发出一束波长为λ0，宽度为tP的脉冲,经准直扩束后垂直射入大气，光脉冲在通过大气时受到散射和衰减，其后向散射光被接收面积为Ar的望远镜系统接收，高度z处的后向散射信号功率可用雷达方程表示[36]：zTzzAr2tczYPzP22p0，，，（1.1）式中，λ为接收到的散射信号的波长，P0为发射的激光脉冲的峰值功率，Y(z)为发送器与接收器光路的几何重叠系数，c为光速，z，为大气体积后向散射系数，zT，为大气的透过率，由朗伯定律可知:’’’，，z0dzz-expzT（1.2）式中，z，为大气的消光系数。从理论上，z，为大气数密度N(z)与散射截面dd的乘积，即ddzNz，(1.3)如前所述，大气中与激光雷达脉冲相互作用产生的后向散射信号成分包括了大气气体分子和气溶胶粒子，由于分子尺寸小，所产生的散射光相对较弱，瑞利散射截面与激发波长的四次方成反比，大气气溶胶粒子对激光的散射光为米散射。对某一激光雷达的特定波长，分子散射直接随着大气分子的浓度而发生变化，但气溶胶散射却很复杂，取决于粒径分布以及气溶胶粒子的折射率。这些气溶胶粒子随地域、时间变化明显，所以无法对其准确估算与预测。而分子散射却可以相对较准确地估计出，只要通过标准大气信息或观测点上空的大气温度、压力分布数据即可获得。因此，实际大气体积后向散射系数和大气消光系数就包含了两个部分:分子散射部分和气溶胶散射部分，即:zzzma(1.4)zmzza（1.5）下标a表示气溶胶，m表示大气分子。接收信号如图1-4所示，在0-A段，发射光束还没进入接收视场，接收不到后向散射信号，探测到的信号强度为0。在A-B段，发射光束逐渐进入视场，发射器与接收器的重叠系数逐渐增大，接收到的后向散射信号也逐渐增强，至B点时，重叠系数Y(z)=1接收到的信号最强。此后尽管视场内继续充满发射光束，但由于光速在大气中传输时按-2z规律发散，因此回波信号也将按此规律减小，如图中B-C段，直至遇到气溶胶密集的介质。在遇到气溶胶密集的介质时，回波信号增强，在某一距离点(图中D点)处达到最大，随后又会逐渐衰减至系统无法探测为止。图1-4激光雷达接收的信号Fig.l一4Lidarreceivedsignal1.3.4激光雷达方程的解激光雷达方程提供了激光回波信号与被探测物的光学性质之间的函数关系，因此可以根据激光雷达探测到的回波信号，通过求解激光雷达方程，获得有关大气性质的信息。但从前面激光雷达方程可知，除了激光雷达系统的光电探测器接收到的大气回波信号功率，以及激光器参数和系统常数以外，方程中还有四个未知量:分子和气溶胶粒子后向散射系数zm，za以及它们的消光系数zm，za。由于分子散射较稳定，其散射系数zm和消光系数zm可以通过标准大气模型或观测点上空的大气温度、压力分布数据获得，因而雷达方程中就剩下气溶胶的两个系数za和za为未知。为了求解这个方程，人们想出许多办法，激光雷达方程的求解方法主要有:斜率法、Klett法和Fernald法。a.斜率法在均匀的大气层中，整个后向散射系数z和消光系数z在该范围内是常数，在这种情况下，z可以用简单的斜率法求得。将雷达方程((l.1)简化如下:zdzzezzCzP02-2(1.6)将P(z)乘以2z，再取自然对数得:zdzzzCzzPzD022lnln(1.7)对D(z)求导得:zdzzdzdzzdD21(1.8)如果大气均匀，则z和z均为常数，0dzzd,z由下式给出:dzzdDz21(1.9)即对曲线D(z)进行最小二乘法拟合，求出曲线的斜率即可求出均匀大气的消光系数。斜率法由于要求均匀大气作为条件，实际是很难做到的，因为z和z在某一高度上通常并不是常数，因而这种方法很难达到很高的精度，但它适合于水平探测。b.Klett法为了精确求解雷达方程，必须知道z与z之间的关系。Klett提出，对于弹性散射，假定两者之间满足如下关系:kB(1.10)这里，B和k是与激光雷达波长和气溶胶的折射率以及粒径分布有关的系数。根据经验值，k的取值范围为:3.167.0k，一般情况下，取k=1,将其代入((1.8)式，得:zdzzdzdzzdD21(1.11)取远方某一边界值高度zc，对应的边界值为czD,cz,czKlett求得了该方程的稳定解为:ccczzdzzDzzDzDz2(1.12)ccczzdzzDBzzDzDz2(1.13)上两式中，由于分母中积分项的符号为正，随着积分数值的增加，边界值和噪声等的影响会减少，因此方程的解很容易保持稳定。这种解法的优点在于只要在边界值高度附近大气透过率czT较小，对于较粗略估算出的cz也可反演出相对较精确的z。一般情况下，cz必须以较精确的方法估算出，才能得到精确的z值，通常边界值的灵敏度会随大气的浑浊度和边界值的取值位置改变而变化很大。c.Fernald法当激光在大气中传输时，要受到大气分子和气溶胶粒子的共同作用，因而雷达方程中的z和z实际应包括两部分，即分子散射部分和气溶胶散射部分。Klett法求解得出的是大气总的消光系数，由于米散射信号的强度与辐射波长的一(1~2)次方成比例，瑞利散射信号的强度与辐射波长的一4次方成比例，在波长较长或有云或气溶胶浓度较大的情况下，大气回波信号中米散射信号占据主要成分，而瑞利散射信号相对很弱，可以忽略，在这种只需考虑单一成分的情况下，使用Klett法求解最有效。然而，在波长较短，对流层中气溶胶浓度不大，或观测高层卷云和气溶胶时，瑞利散射信号的影响就不能忽略了。Fernald法在雷达方程中将分子散射和气溶胶散射分开来考虑，即:zzzam(1.14)zzzam(1.15)下标a表示气溶胶，m表示大气分子，代入激光雷达方程:zamampdzzzzzzAtczYPzP02002exp2(1.16)对于气溶胶引起的米散射，其消光系数与后向散射系数的比设为zzSaa1;对于分子引起的瑞利散射，其消光系数与后向散射系数之比满足如下关系:SzzSmm382(1.17)代入雷达方程式((1.16)，并将方程两边同乘以2z，得到:ZZmpzampdzzSSdzzSzActzYPdzzSzSzActzYPzPz001210001200222exp22exp2(1.18)dzdzzSdSActzYPdzzSzActzYPdzzSSzPzzpZpZm01100010001222exp2122exp22exp(1.19)将上式两边取积分，得:12exp202exp22exp20100010000''1221ZpZpZZmdzzSActzYPdzdzzdzzSdActzYPdzdzzSSzPzS(1.20)再将上式两边取对数，得：dzzzSActzYPdzzSActzYPdzdzzSSzPzSActpZpZZp022ln2expln2ln2exp22zYPln100010000''122100(1.21)上式两边对Z求导，得:zdzdzzSSzzPSActzYPzSSzzPzZmpZ02exp22dz2exp0''1221000122(1.22)选取边界点高度cz，假设已知cz处的后向散射系数cz和消光系数cz，代入上式，并设2zzPzD，20ccczzPzD，得到:ZZdzZZdzzSSzDSzzDdzzSSzDCCmcZm''12100122exp22expz(1.23)则气溶胶的后向散射系数为:ZZdzZZdzzSSzDSzzDdzzSSzDzCCmcZmma''12100122exp22expz(1.24)根据气溶胶消光系数与后向散射系数的关系，即可求出气溶胶消光系数:(1)边界点cz以下高度的气溶胶消光系数为(后向积分):(1.25)czCmcmcacCmmazdzZzdzzSSzDzSSzzDZzdzzSSzDzSSz''21210212112exp212exp(2)边界点zc以上高度的气溶胶消光系数为(前向积分):(1.26)zczcmcmcaccmmazdzzdzzSSzDzSSzzDzzdzzSSzDzSSz''21210212112exp212exp若雷达信号中气溶胶信号可以忽略，上两式的边界值(分母第一项)则由瑞利散射信号和大气密度决定。对于瑞利散射信号，S2是个常数;对于米散射信号，S1与气溶胶或云的折射率、粒径分布和形状等有关，取值范围一般在10sr~100sr之间。大气分子的消光系数可以根据美国标准大气模型确定。