激光雷达测气溶胶消光系数实验研究.doc

1、激光雷达测气溶胶消光系数实验研究摘要：大气气溶胶影响着天气和气候的变化，通过用激光雷达对水平大气中的气溶胶进行连续观测，得到大气气溶胶浓度的高度分布数据，用 Klett 法反演和斜率法得到了气溶胶消光系数数值并利用 MATLAB 程序用计算机对所得实验数据快速方便地直接得出出测量结果和图示。关键词： 气溶胶；激光雷达；探测；Klett 反演算法；斜率法；消光系数；MATLAB0 前言本文为米散射激光雷达测大气气溶胶消光系数的研究报告，包含激光雷达的原理、斜率法等内容。大气气溶胶是指悬浮在大气中的固体和（或）液体微粒，气溶胶的范围很广，如地面的扬尘、烟粒、微生物、植物的孢子和花粉以及水和冰组成的

2、云雾滴、冰晶和雨雪等粒子，都属于气溶胶。除了自然界自然产生的气溶胶粒子以外，人为排放到大气内的污染气体等也会最终形成为气溶胶粒子。气溶胶粒子会直接和间接地影响着气候，通过吸收和散射太阳辐射以及地球的长波辐射而影响着地球和大气的辐射收支，凝结成核参与云的形成，从而影响到天气和气候的变化。 1另外，大气气溶胶的浓度变化还会直接影响到人们的健康和生存环境 2。激光雷达是一种主动的遥感探测工具，已有多年的历史，已广泛应用于激光大气传输、全球气候预测、气溶胶辐射效应及大气环境等研究领域 3。激光雷达探空方法是向大气中发射激光束（微米波） ，利用大气中的气溶胶或大气分子为媒介，进行大气遥感探测。由于激光的

3、波长较短且脉冲的宽度窄，故可进行全大气层内的高精度及高时空分辨率探测。数十年来，激光技术的不断发展为激光雷达大气气溶胶探测提供了所需要的光源。另一方面，信号探测和数据采集极其控制技术的发展使激光雷达在大气气溶胶的探测高度、空间分辨率、时间上的连续监测和测量精度等方面具有全面的优势，是其他探测手段不能比拟的 4。1 研究目的及背景大气中，尘埃、烟雾、云团等气溶胶粒子对大气的化学过程、辐射平衡、气候变化乃至人们的日常生活都有着非常重要的影响。因此，对大气气溶胶粒子的光学特性的探测研究一直是大气科学、气象探测和环境保护的一项重要任务。近年来，中国经济的飞速发展已受到全世界的关注。然而，这种快速的经济

4、增长也伴随着社会体系的变革，高度的工业化和城市化造成许多气溶胶粒子和温室气体被排放到大气，带来了一系列的环境问题，对可持续发展有着严重的负面影响，同时对人们的日常生活和身体健康存在着严重的威胁。如何获取环境变化的第一手资料，准确地提供大气物性及其变化趋势，是当前环境测量领域的一项迫切任务。激光雷达作为一种主动遥感探测工具，已广发用于大气气溶胶辐射效应、大气环境等研究领域 5。在其探测高度、空间分辨率、时间上的连续监测和测量精度等方面具有的独到优势，是其他探测手段不能比的。 利用激光雷达探测城区对流层气溶胶工作已在我国北京、合肥、苏州、西藏那曲等部分地区开展，然而环境监测的数据信息还远远不够，对

5、我国的大气环境状况的评估需要大量的气溶胶监测数据的支撑。如西部地区特殊的地理位置和气候条件及进来的高度的经济发展状况，造成大气中的气溶胶粒子含量较高，空气污染较为严重。对西部地区上空开展气溶胶时空分布探测及其输送特性研究，一直是气象、环保部门所关注的研究课题。特别是起源于我国北方地区频繁发生的沙尘暴事件已引起国内国际的关注，沙尘暴已成为一个重要的地球环境问题。但地面探测仪如太阳光度计、粒子取样计等仪器还无法实现对气溶胶时空分布的剖面极其输送特性的实时探测，也很难实现长期观测。而利用米散射激光雷达实际探测沙尘气候，判断沙尘的相对浓度分布，为利用激光雷达技术进一步研究沙尘的发生及输送规律等奠定基础

6、。2 激光雷达的原理2.1 激光雷达的基本原理及方程激光雷达其最根本的工作原理与普通雷达相类似，即由发射系统发送一个信号，用接收系统采集与目标相作用所产生的返回信号并进行处理。不同之处为，激光雷达发射的信号为激光束，相比于普通无线电雷达所发送的毫米波的波长要短得多，故激光雷达系统的激光波长一般在微米量级，因而只能用于探测极细小的微粒和分子。图 1 激光雷达原理激光雷达的基本原理如图 1 所示，由激光器发出一束波长为 ，宽度为 脉冲0 的经准直扩束后垂直射入大气，光脉冲在通过大气时所受到散射和衰减，其后向散射光被接收面积为 Ar 的望远镜系统接收，高度 z 处的后向散射信号功率可用雷达方程表示6

7、：P( (2.1),)=0()22(,)2(,) 式中， 到的散射信号的波长， 为发射的激光脉冲的峰值功率，Y(z)为发为 接收 0送器与接收器光路的集合重叠系数，c 为光速，为大气的透过率，由朗伯定律可知：(,)为 大气体 积 后向散射系数 ,(,)(2.2)(,)=0(,)式中, 为大气的消光系数。从理论上， 为大气数密度 N(z)与散,z),(z射截面 的乘积，即d/(2.3)dzN)(,(如前所述，大气中与激光雷达脉冲相互作用产生的后向散射信号成分包括了大气气体分子和气溶胶粒子，由于分子尺寸小，所产生的散射光相对较弱，瑞利散射截面与激发波长的四次方成反比，大气气溶胶粒子对激光的散射光为

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

9、射光束逐渐进入视场，发射器与接收器的重叠系数逐渐增大，接收到的后向散射信号也逐渐增强，至 B 点时，重叠系数Y(z)=1，接收到的信号最强。此后尽管视场内继续充满发射光束，但由于光速在大气中传输时按 规律发散，因此回波信号也将按此规律减小，如图 2 中 B-C 段，直至2z遇到气溶胶密集的介质。在遇到气溶胶密集的介质时，回波信号增强，在某一距离点(图中 D 点) 处达到最大，随后又会逐渐衰减至系统无法探测为止。图 2 激光雷达接收的信号2.2 用斜率法求激光雷达方程的解激光雷达方程提供了激光回波信号与被探测物的光学性质之间的函数关系，因此可根据激光雷达探测到的回波信号，通过求解激光雷达方程，获

10、得有关大气性质的信息。但从前面激光雷达方程可知，除了激光雷达系统的光电探测器接收到的大气回波信号功率，以及激光器参数和系统常数以外，方程中还有 4 个未知量：分子和气溶胶粒子后向散射系数 (z)、 (z)以及它们的消光系数 (z)、 (z)。由于分子散射较稳 定，其散射系数 (z)和消光系数 (z)可以通过标准大气模型或观测点上空的大气温 度、压力分布数据获得，因而雷达方程中就剩下气溶胶的两个系数 (z)和 (z)未知。 激光雷达的求解方法主要有：斜率法、Klett 法和 Fernald 法。这里对斜率法求解进行解释，下文会提及 Klett 反演算法。在均匀的大气层中，整个后向散射系数 (z)

11、和消光系数 (z)在该范围内是常数，在 这种情况下， (z)可用简单的斜率法求得。将雷达方程(1.1)简化如下：P(z)=C (2.6)2()20()将 P(z)乘 2,再取自然 对 数得：(2.7)()=()2=()20()对 D(z)求导得：-2 (2.8)()=1()()()如果大气均匀，则 均为常数， /dz=0，()和 () () ()则为 下式：(2.9)()=12()即对曲线 D(z)进行最小二乘法拟合，求出曲线的斜率即可求出均匀大气的消光系数。斜率法由于要求均匀大气作为条件，实际是很难做到的，因为 在某一高()和 ()度上通常并不是常数，因为这种方法很难达到很高的精度，但它适合

12、于水平探测。3 米散射激光雷达实验系统3.1 系统结构示意图米散射激光雷达实验系统构成示意图如图 3，系统采用 Nd：YAG 脉冲激光器作为光源，为了提高白天测量的能力并考虑到入眼安全，选用其三次谐波波长 355nm 作为探测波长。激光束经准直扩束后垂直射入大气，大气的后向散射光由视场角为0.1mrad、直径为 250mm 的望远镜接收，并耦合进入多模光纤。光纤输出信号经准直后入射到一块高光谱分辨率光栅（光栅刻划密度为 2400gr/mm，谱分辨率为 6pm） ，光栅衍射后的大气回波信号经过透镜聚焦，在焦平面上形成光谱分布，利用小孔光阑，分离出主要的大气回波信号（米散射和瑞利散射） ，同时剔除

13、大部分太阳背景及非弹性散射信号，最后米散射和瑞利散射信号由光电倍增管监测。这里，光栅、透镜和光阑组成了一个带宽为 1nm 的分光系统，光电倍增管探测到的散射信号再送入计算机进行数据处理与参数反演。激光雷达系统参数示于图 3。图 3 米散射激光雷达实验系统构成3.2 系统实物图图 4 激光束图 5 实验装置4 气溶胶消光系数的反演算法激光雷达接受的高度 z 处的大气后向散射回波信号能量 P(z)由激光雷达方程决定：P(z)= 02()20()(4.1)式中 激光发射的功率，C 为激光雷达系统常数，z 为探测距离， 和 分别为0为 () ()高度 z 处大气总的后向散射系数和消光系数。在利用激光雷

14、达信号反演大气消光系数的算法中，比较常用的有 Klett 法和Fernald 法，对一般大气环境，常采用 Klett 法反演大气消光系数。Klett 法中假设气溶胶后向散射系数 和消光系数 之间存在如下指数关系： (4.2)=其中，B 和 k 与激光雷达波长以及气溶胶粒子的性质和尺度谱分布有关，k 的取值在0.67-1.3 之间，取 k=1，令 S=1/B，则有 =S,S 称为激光雷达比，这里假设 S=40sr。/设激光雷达回波信号的距离平方修正函数为：D(z)=P(z) (4.3)2若事先已知某一边界值高度 则处 大气 总 的消光系数 ()或后向散射系数 (),处 以下各高度上的气溶胶粒子消

15、光系数 ()或后向散射系数 ()分 别为 ：(4.4)()= ()()()+2()()(4.5)()= ()()()+2()()式中 分别为大气分子消光系数和后向散射系数。() 和 ()上两式中，高度 z 处的大气分子消光系数 和后向散射系数 可以通过美() ()国标准大气模型获得。如果测量高度较高，可选取一段近乎不含气溶胶粒子的清洁大气层，在这段高度范围内，对激光雷达距离校正对数回波信号进行最小二乘拟合，根据 Collis 的斜率法可知，该回归曲线斜率的一半为大气分子的消光系数，从而得到大气分子的消光系数高度分布模型 。这种取值方法更加切合当时的大气状况，取其()中某一高度作为边界值高度 ，

16、此时边界值 中就只含有大气分子的成分。 ()或 ()如果测量高度不够高，就选取一段较为均匀的大气层，对该高度范围内激光雷达距离校正对数回波信号进行最小二乘拟合，回归曲线斜率的一半即为高度 处的边界值。()5 用 MATLAB 对实验数据进行处理5.1 实验数据存储格式1.原始数据文件名为每组数据采集的时间， “年月日时分秒_A.txt”表示，如“2014 年 10 月 14 日 195602_A.txt”.2.数据格式：每个 txt 文件只有一列数据，代表采集信号，时间分辨率为 1 分钟，距离分辨率为 15 米。5.2 MATLAB 代码5.2.1 将 txt 格式的原始数据转换为 xls 格

17、式的表格形式clc;clear;cd E:原始数据 1;data1=dir;for i=1:length(data1)-2data2i=importdata(data1(i+2).name);%读取文件所有的数据time(i)=data1(i+2).datenum;%获取文件保存时间endm,n=size(data2);%得到 data2 的横列大小Ttop=70;%最大探测距离对应的时间点（ us）t=0.1:0.1:Ttop;%设置时间t0=0;%延时h=(t+t0)*150/1000;%将时间转换成高度for j=1:nfor i=1:msign(:,j)=data2i,j;%将 dat

18、a2 的全部数据存入一个数组中 end endfilename=E:数据 1.xls;xlRange1=B1;xlRange2=A2;xlRange3=B2;xlswrite(filename,time,sheet1,xlRange1);xlswrite(filename,h,sheet1,xlRange2);xlswrite(filename,sign,sheet1,xlRange3);5.2.1 求水平大气消光系数clc;clear;%清理内存和屏幕close all% 数据读取filename=E:实验报告附件 3.xls;time=xlsread(filename,Sheet1,B1:

19、K1);data_time=datestr(time,HH:MM);range=xlsread(filename,sheet1,A2:A701);p=xlsread(filename,Sheet1,B2:K701);m,n=size(p);%得到 p 的横列大小hg=300; %背景高度对应的序列数% 扣除背景和距离平方修正for j=1:nfor i=1:mp(i,j)=p(i,j)-mean(p(670:700,j);Pz2(i,j)=p(i,j)*range(i).2;endend% 作垂直方向的平均for i=2:m-1lidar_data_v(i-1,:)=(Pz2(i-1,:)+P

20、z2(i,:)+Pz2(i+1,:)/3;end% 作水平方向的平均for i=2:n-1lidar_data_h(:,i-1)=(lidar_data_v(:,i-1)+lidar_data_v(:,i)+lidar_data_v(:,i+1)/3;%作 3 点平滑endlidar_data_ln=log(lidar_data_h);%取自然对数% 线性拟合h1=40;h2=60;for j=1:n-2 index(:,j)=polyfit(range(h1:h2),lidar_data_ln(h1:h2,j),1);lidar_data_nihe(h1:h2,j)=polyval(inde

21、x(:,j),range(h1:h2); endfigureplot(range(h1:h2),lidar_data_ln(h1:h2,2),*,range(h1:h2),lidar_data_nihe(h1:h2,2),-);% 求水平大气消光系数alpha1_110=-(mean(index(1,:)/2;6 总结通过用激光雷达对水平大气中的气溶胶进行连续观测，得到大气气溶胶浓度的高度分布数据，用 Klett 法反演和斜率法得到了气溶胶消光系数数值并利用 MATLAB 程序用计算机对所得实验数据快速方便地直接得出出测量结果和图示。总体实验流程操作并不复杂困难，易于实现，且具有一定的实用价值

22、。随着该项实验研究的继续深入以及将来观测数据的不断积累，该研究成果将大力有助于大气环境颗粒状气溶胶的发生、传输等时空特性，以及大气辐射和城市热岛等城市气息特征的产生机理、传输规律的研究，并为这些研究工作提供实时的科学数据支撑。参考文献： 1 刘君. 大气温度及气溶胶激光雷达探测技术研究D. 2008. 1.2 盛裴轩, 毛节泰, 李建国等. 大气物理学M.北京：北京大学出版社，2003: 25-28.3 钟志庆，周军，戚福第等. 探测大气气溶胶消光系数的便携式米散射激光雷达J.强激光与粒子束, 2003,15(12).4 周军, 岳古明, 戚福第等. 大气气溶胶光学特性激光雷达探测J. 量子电子学报, 1998, 15(2).5 毛敏娟, 吴永华, 戚福第等. 车载式 1064nm 和 532nm 双波长米散射激光雷达J.强激光与粒子束,2005, 17(5):677-680.6 J. A. Reagan, M. P. Mccormick, J. D. Spinhirne. Lidar sensing of aerosols and clouds in the troposphere and stratosphereJ. Proceedings of the IEEE, 1989, 77(3): 433-447.