本发明专利技术涉及一种使用纯转动拉曼激光雷达测量大气温度的方法，属于测量技术领域。本发明专利技术方法首先获取纯转动拉曼激光雷达的原始高低阶回波信号，并得到去除背景的高低阶通道回波信号；然后判断是否出现近地面光子数饱和现象，如果出现，则通过采用选取不同倍增系统的电子鉴别水平μ值，来实现对近地面光子数饱和的修正，这样就能够从光子计数非线性的深层次原因，即脉冲高度分布和电子鉴别水平，来分析出现近地面饱和现象的原因，利用修正后的高低阶通道回波信号进行温度反演过程，提高纯转动拉曼激光雷达测量大气温度廓线的精度。
【技术实现步骤摘要】
本专利技术涉及，属于测量
。
技术介绍
当前对大气温度探测的激光雷达方法主要有瑞利散射密度法、差分吸收光谱法、振动拉曼光谱法和纯转动拉曼光谱法等。瑞利散射密度法是基于瑞利散射强度与大气分子密度之间的依存关系来反演大气分子密度，其适用条件是洁净的空气，认为大气透过率函数近似为1，低层大气中较高浓度的气溶胶使得其适用范围为30公里以上的高层大气。差分吸收光谱法是基于氧气的吸收光谱，要求高稳定性的两束激光不同的频率，通过探测氧气吸收光谱内和吸收谱外的信号来反演大气分子溶度，其要求激光器参数较高而忽视了大气中各种散射的影响。振动拉曼光谱法是基于振动拉曼散射强度与温度的依存关系，其信号较弱，信噪比较低，对后续的接收设备有很高的要求。同时低层大气气溶胶分布的剧烈变化会给大气透过率校正带来影响，也导致其不能很好的测量大气温度。转动拉曼光谱法是基于转动拉曼散射强度与温度之间的关系，通过检测散射光中心频率一定偏移处的回波信号来反演大气温度。回波信号为光子量级，使得此种方法适用于低层大气。气溶胶和分子消光项除了影响回波信号的强度外，对测量精度基本没有影响，能够很好的测量低层大气温度。 1972年，Cooney提出了利用转动拉曼激光雷达测量大气温度，其基本原理是基于转动拉曼谱线强度与温度之间的关系来进行大气温度廓线的反演。在对于对流层中的温度进行反演的激光雷达中，转动拉曼激光雷达已经成为可靠性最高的。同时，在测量对流层温度廓线时，其仍然有较高的空间分辨率。现今，测量对流层内大气温度廓线的光源波长主要是532nm和355nm波长，前者主要是由于其高的激光能量和接收效率而被接受，而后者主要是因为转动拉曼后向散射系数与波长的四次方成反比而被大家所推崇，同时紫外波段的光源对于眼睛的伤害更小。 对于大气温度廓线的反演过程主要是有两种方法，即一次项式拟合和二次项式拟合，在拟合的过程中，都会采用一定高度区间的回波数据与无线电探空仪测得的探空温度曲线或者是已经定标的卤素灯来完成定标过程，求解出相应的参数来完成温度廓线的反演过程。在修正技术中，主要是从几何因子和反演温度与探空温度之间差异来实现。对于近地面的强回波信号在经过倍增系统后，考虑到脉冲高度分布和电子鉴别水平的选择这两个因素，会出现光子计数器计数非线性的情况(饱和现象)。德国的RAMSES雷达系统考虑采用的是在高度层进行分层采集的形式来避免饱和现象的产生，对近地面层和高层采用不同的采集系统，同时采用不同的衰减片，这样在保证信号信噪比的同时，能够有效的降低饱和现象的发生率，但该方法的缺陷是对近地面的大气温度测量不准确。另外，美国NASA的GSFC和CNRS两个雷达系统采用的是饱和修正的方式，GSFC系统依靠的是系统的时间分辨率，而CNRS系统依靠的是探测系统的计数速率的最大值，而这两项参数在实际的实验中很难能够准确的得到，因此对近地面大气温度的测量误差较大。
技术实现思路
有鉴于此，本专利技术提供了一种基于纯转动拉曼激光雷达测量温度时近地面光子数饱和的修正方法，本方法能够在近地面出现光子数饱和的情况下，通过修正高阶回波信号廓线和低阶回波信号廓线，提高纯转动拉曼激光雷达测量大气温度廓线的精度。 本专利技术的目的是通过下述技术方案实现的。 ，其操作步骤为: 步骤1:获取纯转动拉曼激光雷达的原始高阶回波信号(用符号P1表示)和原始低阶回波信号(用符号P2表示)； 步骤2:将纯转动拉曼激光雷达获得的原始高阶回波信号P1的后5公里回波信号的平均值作为高阶通道的背景噪声(用符号Ptll表示)；将纯转动拉曼激光雷达获得的原始低阶回波信号P2的后5公里回波信号的平均值作为低阶通道的背景噪声(用符号Ptl2表示)； 步骤3:用原始高阶回波信号P1减去高阶通道的背景噪声Ptll，得到去除背景的高阶通道回波信号(用符号Pltl表示)，即Pltl = P1-Poi ;用原始低阶回波信号P2减去低阶通道的背景噪声Ptl2，得到去除背景的低阶通道回波信号(用符号P2tl表示)，即P2tl = P2-P02 ； 步骤4:判断是否出现近地面光子数饱和现象，如果出现，则进行步骤5的操作。否贝U，进行步骤6的操作。 所述判断是否出现近地面光子数饱和现象的方法为: 步骤4.1:求解出I公里以下高度层内高低阶回波信号比值的最小值以及对应的高度；所述高低阶回波信号比值是高阶回波信号和低阶回波信号的比值。 步骤4.2:计算得到3公里处高低阶回波信号比值； 步骤4.3:绘制3公里以内高低阶回波信号比值与高度曲线图；所述曲线图位于二维直角坐标系上，横坐标轴为高低阶回波信号比值，纵坐标轴为高度值。在所述曲线图上，找到步骤4.1得到的I公里以下高度层内高低阶回波信号比值的最小值对应的点以及步骤 4.2得到的3公里处高低阶回波信号比值对应的点，并进行连线，将该连线作为基准线。若在3公里高度层以下，高低阶回波信号比值廓线的一部分位于基准线的右侧，则判断出现近地面光子数饱和现象。 步骤5:考虑系统几何因子对于系统采集回波信号的影响，对于1.5-12公里内，去除背景的高阶通道回波信号Pltl和去除背景的低阶通道回波信号P2tl进行修正，具体过程为: 步骤5.1:得到修正高度范围内高阶通道的最大回波光子数(用符号Nmax表示)和低阶通道的最大回波光子数(用符号N, _表示)。 步骤5.2:通过公式⑴修正预设高度范围内高阶通道的回波信号，得到修正后的去除背景高阶通道回波信号(用符号P, K1表示)。 本文档来自技高网...
【技术保护点】
一种针对纯转动拉曼激光雷达测量温度的饱和修正方法，其特征在于：其操作步骤为：步骤1：获取纯转动拉曼激光雷达的原始高阶回波信号P1和原始低阶回波信号P2；步骤2：将纯转动拉曼激光雷达获得的原始高阶回波信号P1的后5公里回波信号的平均值作为高阶通道的背景噪声P01；将纯转动拉曼激光雷达获得的原始低阶回波信号P2的后5公里回波信号的平均值作为低阶通道的背景噪声P02；步骤3：用原始高阶回波信号P1减去高阶通道的背景噪声P01，得到去除背景的高阶通道回波信号P10，即P10＝P1‑P01；用原始低阶回波信号P2减去低阶通道的背景噪声P02，得到去除背景的低阶通道回波信号P20，即P20＝P2‑P02；步骤4：判断是否出现近地面光子数饱和现象，如果出现，则进行步骤5的操作；否则，进行步骤6的操作；步骤5：考虑系统几何因子对于系统采集回波信号的影响，对于1.5‑12公里内，去除背景的高阶通道回波信号P10和去除背景的低阶通道回波信号P20进行修正，具体过程为：步骤5.1：得到修正高度范围内高阶通道的最大回波光子数Nmax和低阶通道的最大回波光子数N′max；步骤5.2：通过公式(1)修正预设高度范围内高阶通道的回波信号，得到修正后的去除背景高阶通道回波信号P′10；P10′=(1-μ)P10exp[-P10Nmax]+(μ-μ22!)P102Nmaxexp[-P10Nmax]---(1)]]>其中，μ为电子鉴别水平参数，μ∈[0,1]；exp[·]表示以自然常数e为底的指数函数；步骤5.3：通过公式(2)修正高度范围内低阶通道的回波信号，得到修正后的去除背景低阶通道回波信号P′20；P20′=(1-μ)P20exp[-P20Nmax′]+(μ-μ22!)P202Nmax′exp[-P20Nmax′]---(2)]]>步骤6：利用小波分析方法对修正后的去除背景高阶通道的回波信号P′10进行滤波操作后，得到去噪后的高阶通道回波信号P″10，利用小波分析方法对修正后的去除背景低阶通道的回波信号P′20进行滤波操作后，得到去噪后的低阶通道回波信号P″20步骤7：选取3至9公里范围内的高度层作为定标高度层，将定标高度层内无线电探空仪测得的温度廓线与定标高度层内的去噪后的高低阶回波信号比值曲线，通过公式(3)在定标高度层内进行二次项式拟合；A[ln(P10′′P20′′)]2+B[ln(P10′′P20′′)]+C=1Tcli---(3)]]>其中，A、B和C为3个定标系数值，可通过公式(3)拟合得到；Tcli为定标高度层内无线电探空仪测得的温度；步骤8：利用步骤7得到的A、B和C定标系数值，通过公式(4)对K公里以下高度内进行温度廓线的反演计算，K∈[9,12]；T=1A[ln(P10′′P20′′)]2+B[ln(P10′′P20′′)]+C---(4)]]>其中，T为0至K公里范围内的反演温度；经过上述步骤的操作，即可得到纯转动拉曼激光雷达测量温度经过修正后的温度廓线。...
【技术特征摘要】
1.一种针对纯转动拉曼激光雷达测量温度的饱和修正方法，其特征在于:其操作步骤为: 步骤1:获取纯转动拉曼激光雷达的原始高阶回波信号P1和原始低阶回波信号P2 ；步骤2:将纯转动拉曼激光雷达获得的原始高阶回波信号P1的后5公里回波信号的平均值作为高阶通道的背景噪声Ptll ;将纯转动拉曼激光雷达获得的原始低阶回波信号P2的后5公里回波信号的平均值作为低阶通道的背景噪声Ptl2 ； 步骤3:用原始高阶回波信号P1减去高阶通道的背景噪声Ptll,得到去除背景的高阶通道回波信号Pltl,即P1O = P1-Poi ;用原始低阶回波信号P2减去低阶通道的背景噪声Ptl2，得到去除背景的低阶通道回波信号P2Q，即P2tl = P2-P02 ； 步骤4:判断是否出现近地面光子数饱和现象，如果出现，则进行步骤5的操作；否则，进行步骤6的操作； 步骤5:考虑系统几何因子对于系统采集回波信号的影响，对于1.5-12公里内，去除背景的高阶通道回波信号Pltl和去除背景的低阶通道回波信号P2tl进行修正，具体过程为:步骤5.1:得到修正高度范围内高阶通道的最大回波光子数Nmax和低阶通道的最大回波光子数N' max; 步骤5.2:通过公式(I)修正预设高度范围内高阶通道的回波信号，得到修正后的去除背景高阶通道回波信号P' 10;其中，μ为电子鉴别水平参数，μ e [O, I] ;exp[.]表示以自然常数e为底的指数函数； 步骤5.3:通过公式(2)修正高度范围内低阶通道的回波信号，得到修正后的去除背景低阶通道回波信号P' 20；步骤6:利用小波分析方法对修正后的去除背景高阶通道的回波信号P' 1(1进行滤波操作后，得到去噪...
【专利技术属性】
技术研发人员：张寅超，王中昆，陈思颖，陈和，郭磐，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。