基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法技术

本发明专利技术公开了一种利用激光雷达的量数据计算整层气溶胶光学厚度的计算方法，可有效解决白天多云情况和夜晚气溶胶光学参数无法获得的难题。一般利用激光雷达计算气溶胶光学厚度的方法是对消光系数积分得到的，但激光雷达的探测高度有限，特别是微脉冲激光雷达，一般白天探测3‑5公里的高度，夜晚探测5‑10公里的高度，远小于整层大气的高度。为了解决该问题，本发明专利技术提出了一种基于激光雷达垂直与水平测量数据计算整层大气气溶胶光学厚度的方法。激光雷达水平测量出近地面气溶胶消光系数，垂直测量数据计算出气溶胶标高，然后解算出整层大气气溶胶光学厚度。该方法为白天多云天气条件下与夜晚探测整层大气气溶胶光学厚度提供了新的技术手段。

【技术实现步骤摘要】

本专利技术涉及激光雷达大气探测，属于大气光学领域，具体是一种基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法。技术背景大气气溶胶光学厚度(AOD)是表征大气浑浊度和气溶胶含量的一个重要物理量。目前，整层大气气溶胶光学参数主要是在晴朗无云的天气条件下利用太阳光度计测量，作为测量大气光学厚度精度最高的工具得到了广泛的应用，但其无法获得太阳被云遮挡时刻、夜间以及特定高度处气溶胶光学参数。激光雷达作为一种主动遥感工具已有几十年的历史，随着激光雷达技术的发展，已经成为目前对大气、海洋和陆地进行高精度遥感探测的有效手段，广泛用于激光大气传输、全球气候监测、气溶胶辐射效应及大气环境等研究领域。激光雷达可探测气溶胶的垂直分布，打破了以太阳为光源的被动遥感只能探测整层气溶胶特性的局限。激光雷达以其高垂直空间分辨率，高时间分辨率，时间上的连续监测和测量精度等方面的优势在大气气溶胶监测中有重要作用，是其它探测手段不可比拟的。一般计算大气气溶胶光学厚度是通过测量大气后向米散射信号反演得到气溶胶消光系数的高度分布，再通过对气溶胶消光系数的积分得到气溶胶的光学厚度，但激光雷达的探测高度有限，特别是微脉冲激光雷达，一般白天可以探测3-5公里的高度，夜晚可以探测5-10公里的高度，远远小于整层大气的高度。鉴于上述原因本专利技术提出了一种利用激光雷达测量数据计算整层大气气溶胶光学厚度的新方法，即利用气溶胶垂直消光系数廓线计算出气溶胶标高，再结合地面消光系数计算整层大气气溶胶光学厚度。
技术实现思路
本专利技术提供了一种利用激光雷达水平和垂直观测数据反演得到的整层大气气溶胶光学厚度的新方法。针对激光雷达测量得到的数据，根据Fernald方法反演得到消光系数垂直廓线进而得到标高，再结合斜率法反演得到近地面消光系数计算得到气溶胶光学厚度。本专利技术是通过以下技术方案来实现：一种基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法，其特征在于，由激光雷达分别垂直和水平测量大气气溶胶，激光雷达垂直和水平方向测量时间一致，根据探测到的数据进行反演得到大气气溶胶垂直消光系数廓线和大气气溶胶近地面消光系数。将大气气溶胶垂直消光系数廓线按高度分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，并根据这四种类型大气气溶胶垂直消光系数廓线计算得到气溶胶标高，将气溶胶标高结合大气气溶胶近地面消光系数计算整层大气气溶胶光学厚度，具体包括如下步骤：(1)基于激光雷达测量数据获取大气气溶胶近地面消光系数：基于激光雷达分别垂直和水平方向测量时的回波光子数，利用Fernald方法反演激光雷达的回波数据，得到大气气溶胶消光系数的高度分布廓线，利用斜率法反演得到近地面气溶胶的消光系数；(2)大气气溶胶标高的反演：假设大气气溶胶浓度随高度指数衰减时，且气溶胶消光截面不随高度发生变化的情况下，可以得到：σa(z,λ)≈σa(0,λ)exp(-z/Ha)(1)将(1)式对高度积分可得气溶胶光学厚度的计算公式，τa(λ)≈σa(0,λ)Ha(2)其中：σa(z,λ)为高度z处波长为λ时的气溶胶消光系数，σa(0,λ)为波长为λ时近地面的消光系数，τa(λ)是波长为λ时整层大气气溶胶光学厚度，Ha是大气气溶胶标高。一般来讲将步骤(1)中得到的气溶胶消光廓线按照(1)式进行指数拟合，可以得到气溶胶标高，但是通过对气溶胶消光廓线的统计分析得到，在晴朗无云天气条件下，并不是所有的气溶胶消光廓线都随高度指数衰减，因此本专利技术将气溶胶垂直廓线分为指数衰减型、均匀型、近地面高值型和边界层高值型四种类型4种类型，进行分段拟合得到气溶胶标高；(3)大气气溶胶光学厚度的反演：将步骤(2)拟合得到的气溶胶标高和步骤(1)反演得到的近地面消光系数带入(2)式的气溶胶光学厚度的计算公式，计算得到大气气溶胶光学厚度：τa(λ)＝σa(0,λ)Ha。与现有技术相比，本专利技术具有以下有益的技术效果：本专利技术提供了一种利用激光雷达测量数据计算整层大气气溶胶光学厚度的新方法，利用气溶胶标高和近地面消光系数计算得到光学厚度，避免了因激光雷达有效探测高度远小于整层大气高度，直接对气溶胶消光系数积分带来的误差。本专利技术提供了一种利用激光雷达测量数据计算整层大气气溶胶光学厚度的新方法，将标高分为四种类型，对不同时刻气溶胶消光系数廓线按照不同的类型进行分段拟合计算得到气溶胶标高，避免了将气溶胶标高在长时间内看成常数带来的误差。本专利技术原理：使用激光雷达测量得到的回波光子数根据Fernald方法和斜率法反演得到气溶胶垂直消光廓线和近地面气溶胶的消光系数。在假设气溶胶浓度随高度指数衰减且气溶胶消光截面不随高度发生变化的情况下，利用激光雷达的消光系数廓线拟合得到气溶胶标高，结合近地面消光系数计算得到整层大气气溶胶光学厚度。附图说明图1为气溶胶标高随时间的变化图；图2为气溶胶光学厚度随时间的变化图。具体实施方式以下结合具体实施例，对本专利技术进行详细说明。一种基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法，其特征在于，由激光雷达分别垂直和水平测量大气气溶胶，激光雷达垂直和水平方向测量时间一致，根据探测到的数据进行反演得到大气气溶胶垂直消光系数廓线和大气气溶胶近地面消光系数，将大气气溶胶垂直消光系数廓线分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，并根据这四种类型大气气溶胶垂直消光系数廓线计算得到气溶胶标高，将气溶胶标高结合大气气溶胶近地面消光系数计算整层大气气溶胶光学厚度，具体包括如下步骤：(1)基于激光雷达测量数据获取大气气溶胶消光系数：气溶胶消光系数的反演一般都是建立在求解激光雷达方程的基础之上，比较成功地算法有Collis的斜率算法、Klett法和Fernald法。米散射激光雷达探测气溶胶的原理是通过雷达发射机向大气中发射脉冲激光束，再由望远镜接收由大气中分子、粒子产生的后向散射光，通过激光雷达方程求解出气溶胶消光系数等参数。激光雷达方程可表示为P(z)＝Cz-2β(z)T2(z)其中P(Z)是在距离Z处的雷达的回波强度，C是雷达系统常数，包括雷达发射强度、接收系统的衰减、几何校正等。β(z)表示大气的后向散射系数，表示大气的透过率，σ(z)是大气的消光系数。本专利技术水平方向消光系数是由斜率法反演得到的，具体步骤为先对激光雷达方程两边同时取对数得到Ln(P(z)z2)=Ln[Cβ(z)]-2∫0zσ(z)dz]]>假设大气均匀，则全程的β(z)和σ(z)不随Z变化，则上式可表示为Ln[P(z)z2]＝Ln(Cβ)-2σz该式是一个直线方程，以距离校正信号的对数为纵坐标，距离Z为横坐标，用作图法或最小二乘法，可获得一直线，其直线斜率的-1/2就是σ值，由此可得水平消光系数。对于气溶胶垂直消光廓线利用Fernald方法反演。气溶胶垂直消光系数的后向积分方程可表示为σ1(z)=X(z)exp[2(S1S2-1)∫zzcσ2(z)dz]X(zc)σ1(zc)+S1S2σ2(zc)+2∫zzcX(z)exp[(S1S2-1)∫Z本文档来自技高网...

【技术保护点】
一种基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法，其特征在于：由激光雷达分别垂直和水平测量大气气溶胶，根据探测到的数据进行反演得到大气气溶胶垂直消光系数廓线和大气气溶胶近地面消光系数，将大气气溶胶垂直消光系数廓线分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，并根据这四种类型大气气溶胶垂直消光系数廓线计算得到气溶胶标高，将计算得到的气溶胶标高结合大气气溶胶近地面消光系数计算整层大气气溶胶光学厚度，具体包括如下步骤：(1)基于激光雷达测量数据获取大气气溶胶近地面消光系数：基于激光雷达分别垂直和水平方向测量时的回波光子数数据，利用Fernald方法反演激光雷达的回波数据，得到大气气溶胶消光系数的高度分布廓线，利用斜率法反演得到近地面气溶胶的消光系数；(2)大气气溶胶标高的反演：将测量得到的气溶胶垂直方向的消光系数廓线分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，拟合得到大气气溶胶标高；(3)大气气溶胶光学厚度的反演：将步骤(2)拟合得到的气溶胶标高和步骤(1)反演得到的近地面消光系数带入下面的气溶胶光学厚度的计算公式，计算得到大气气溶胶光学厚度：τa(λ)＝σa(0,λ)Ha其中：τa(λ)为波长为λ时整层大气气溶胶光学厚度；σa(0,λ)为波长为λ时大气气溶胶近地面消光系数；Ha是大气气溶胶标高。...

【技术特征摘要】
1.一种基于激光雷达测量数据计算整层大气气溶胶光学厚度的方法，其特征在于：由激光雷达分别垂直和水平测量大气气溶胶，根据探测到的数据进行反演得到大气气溶胶垂直消光系数廓线和大气气溶胶近地面消光系数，将大气气溶胶垂直消光系数廓线分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，并根据这四种类型大气气溶胶垂直消光系数廓线计算得到气溶胶标高，将计算得到的气溶胶标高结合大气气溶胶近地面消光系数计算整层大气气溶胶光学厚度，具体包括如下步骤：(1)基于激光雷达测量数据获取大气气溶胶近地面消光系数：基于激光雷达分别垂直和水平方向测量时的回波光子数数据，利用Fernald方法反演激光雷达的回波数据，得到大气气溶胶消光系数的高度分布廓线，利用斜率法反演得到近地面气溶胶的消光系数；(2)大气气溶胶标高的反演：将测量得到的气溶胶垂直方向的消光系数廓线分成指数衰减型、均匀型、近地面高值型和边界层高值型四种类型，拟合得到大气气溶胶标高；(3)大气气溶...

【专利技术属性】
技术研发人员：陈莎莎，徐青山，徐赤东，
申请(专利权)人：中国科学院合肥物质科学研究院，
类型：发明
国别省市：安徽;34