CN101004453B - 一种气象与大气环境参数的测定方法 - Google Patents
一种气象与大气环境参数的测定方法 Download PDFInfo
- Publication number
- CN101004453B CN101004453B CN200610105194A CN200610105194A CN101004453B CN 101004453 B CN101004453 B CN 101004453B CN 200610105194 A CN200610105194 A CN 200610105194A CN 200610105194 A CN200610105194 A CN 200610105194A CN 101004453 B CN101004453 B CN 101004453B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- atmospheric
- light
- passage
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开的一种气象与大气环境参数的测定方法,利用激光雷达测定系统将大气的后向散射光接收下来,经分光成为该激光雷达所需的各种光信号,再经光电检测后,根据设定的程序进行数据的分析处理,从而实现大气温度、水蒸气密度垂直分布的气象参数测量,以及大气气溶胶消光系数及散射系数、气溶胶光学厚度、大气能见度和非球形粒子后向散射光的退偏振比垂直分布的大气环境参数的测量。
Description
技术领域
本发明属于气象与环境监测技术领域,涉及一种气象与大气环境参数的测定方法,特别涉及利用紫外波长激光雷达系统进行气象和大气环境观测,来测定其参数的方法。
背景技术
在大气环境研究领域中,大气温度、水蒸气密度和大气成分等是进行大气研究的重要参数,比如:大气温度、水蒸气密度或相对湿度,以及大气气溶胶消光系数及散射系数、气溶胶光学厚度及大气能见度、非球形粒子后向散射光的退偏振比等。
在大气参数遥感探测中,大气温度,湿度及气溶胶是相互依存性很高的三个大气参数,比如:影响降雨。降水的相对湿度分布的导出,需要同时测量大气温度数值,而大气温度分布又影响气溶胶的聚集、输送,气溶胶又影响云的形成等等。虽然各参数的单独数据反演技术目前已比较成熟,但利用一个激光雷达系统同时测量数个参数的实用技术还在研究中。本文提出一种利用数个不同性质的大气回波的测量信号,来同时获得温度,水汽密度(相对湿度),气溶胶等气象及大气环境参数的反演技术。
发明内容
本发明的目的是提供一种气象与大气环境参数的测定方法,来进行大气温度、水蒸气密度(相对湿度)垂直分布等气象参数测量,以及大气气溶胶消光系数及散射系数、气溶胶光学厚度、大气能见度和非球形粒子后向散射光的退偏振比的垂直分布等大气环境参数的测量,而且能同时得到气象和大气环境的参数。
本发明所采用的技术方案是,一种气象与大气环境参数的测定方法,将大气中激光脉冲与大气分子、粒子相互作用后产生的后向散射光接收下来,经分光成为该激光雷达所需的各种光信号,再经光电转换后送入计算机,根据设定的程序进行数据的分析处理,从而实现大气温度、水蒸气密度垂直分布的气象参数测量,以及大气气溶胶消光系数及散射系数、气溶胶光学厚度、大气能见度和非球形粒子后向散射光的退偏振比垂直分布的大气环境参数的测量,具体包括以下步骤:
a.向大气中发射脉冲激光束;
b.接收上述脉冲激光与大气中的分子和粒子相互作用后产生的后向散射光;
c.对上述大气后向散射光进行分光,将后向散射光中的转动拉曼谱线、水汽分子的振动拉曼谱线与米-瑞利散射谱线和太阳背景光进行分离,分别得到2个转动拉曼信号,用于温度测量,得到的水蒸气振动拉曼信号,用于水蒸气密度测量,得到的米-瑞利散射信号,用于大气气溶胶光学特性的测量,得到的大气后向散射光信号的偏振特性,即与激光发射频率相同的散射信号的平行与垂直偏振分量,用于大气中非球形气溶胶形态的测量;
d.将上述分离后的各种散射光信号变为电信号接收;
e.将上述处理后的各种信号按照设定的程序分别进行处理,最终得到大气温度、水蒸气密度的垂直分布,以及大气气溶胶消光系数的垂直分布和大气气溶胶光学厚度和非球形气溶胶粒子的后向散射光的退偏振比.
本发明的气象及大气环境参数反演方法主要利用转动拉曼信号反演大气温度;利用水汽振动拉曼信号结合其中一个转动拉曼信号来消除大气透过率的影响,从而实现水蒸气密度的探测;探测到的温度以及水汽密度可以反演相对湿度的垂直分布;利用米-瑞利散射信号高精度反演大气消光系数,从而求得大气气溶胶消光系数、后向散射系数、气溶胶光学厚度、大气能见度等;利用探测到的2个偏振方向垂直的回波信号,反演非球形粒子后向散射光的退偏振比的垂直分布,从而判明非球形矿物质气溶胶及酸性球形城市气溶胶等大气环境参数的分布。所以本发明的有益效果是实现了一套系统同时进行气象参数和大气环境观测的数据反演。
附图说明
图1是实现本发明方法的激光雷达系统结构原理图;
图2是实现本发明的分光与光电检测步骤的一装置实施例;
图3是一种分光系统中边缘反射镜位置的确定示意图;
图4为本发明的数据分析处理总流程;
图5为大气温度的数据分析与处理流程;
图6为气溶胶的数据分析与处理流程;
图7为非球形粒子的数据分析与处理流程;
图8为水蒸气的数据分析与处理流程;
图9为用本发明方法求得的大气温度分布实施例;
图10为用本发明方法求得的水汽密度分布实施例。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
激光遥感的基础是光辐射与大气中的原子、分子以及气溶胶粒子之间相互作用所产生的各种物理过程。各种散射机理中,米散射(Mie scattering)是一种散射谱的中心波长与入射激光波长相同,散射谱的谱宽近似于入射激光谱宽的弹性散射,它是由粒径相当或大于激光波长的气溶胶粒子引起的散射;瑞利散射(Rayleigh scattering)也是一种中心波长与入射激光波长相同,谱宽依存大气温度变化的弹性散射,它是由散射体粒径比激光波长小的分子或原子引起的散射现象,主要用于大气温度、大气分子密度等参数的测量;拉曼散射(Raman scattering)可分为转动拉曼和振动拉曼散射,是一种由大气分子或原子引起的非弹性散射,散射谱分布于入射激光谱线的两侧其散射截面是各种散射机理中较小的一种,需要高效率的分光及检测系统,但由于其特殊的散射机理,很适合用来探测大气温度、水蒸气密度和大气成分。
本发明对气象和大气环境参数的测定方法就是采用拉曼散射机理,利用大气分子(主要是N2,O2)的纯转动拉曼散射强度用于测量大气的温度,利用大气中的水汽分子产生的振动拉曼散射强度探测大气水汽的密度分布,同时结合大气分子的转动拉曼散射信号,反演得到在气象/气候学研究中具有重要作用的相对湿度的分布。
本发明的气象与大气环境参数的测定方法,利用激光雷达测定系统将大气的后向散射光接收下来,经分光成为该激光雷达所需的各种光信号,再经光电转换后送入计算机,根据设定的程序进行数据的分析处理,从而实现大气温度、水蒸气密度垂直分布的气象参数测量,以及大气气溶胶消光系数及后向散射系数、气溶胶光学厚度、大气能见度和非球形粒子后向散射光的退偏振比垂直分布的大气环境参数的测量,具体包括以下步骤:
a.向大气中发射脉冲激光束;
b.接收上述脉冲激光与大气中的分子和粒子相互作用后产生的后向散射光;
c.对上述大气后向散射光进行分光,将后向散射光中的转动拉曼谱线、水汽分子的振动拉曼谱线与米-瑞利散射谱线和太阳背景光进行分离,分别得到2个转动拉曼信号,用于温度测量,得到的水蒸气振动拉曼信号,用于水蒸气密度测量,得到的米-瑞利散射信号,用于大气气溶胶光学特性的测量,得到的大气后向散射光信号的偏振特性,即与激光发射频率相同的散射信号的平行与垂直偏振分量,用于大气中非球形气溶胶形态的测量;
d.将上述分离后的各种散射光信号变为电信号接收;
e.将上述处理后的各种信号按照设定的程序分别进行处理,最终得到大气温度、水蒸气密度的垂直分布,以及大气气溶胶消光系数的垂直分布和大气气溶胶光学厚度和非球形气溶胶粒子的后向散射光的退偏振比。
激光雷达作为一种主动遥感探测工具已广泛用于激光大气传输、全球气候预测、气溶胶辐射效应及大气环境等研究领域,下面通过一种实现本发明方法的激光雷达系统对上述方案作详细说明。
如图1所示,激光雷达系统由发射系统1、接收系统2、分光与光电检测系统3和数据处理系统5构成。发射系统1包括Nd:YAG激光器6和接收Nd:YAG激光器6发出的脉冲激光束并进行扩束和准直的扩束器4,激光器具有能输出其三次谐波分量(λ0=355nm)并发射出去的光学部件和结构,Nd:YAG激光器6发出波长为355nm的脉冲激光束,经扩束器4扩束后,垂直射向大气;激光与大气中的分子和粒子相互作用后产生的后向散射光由望远镜接收系统2接收下来,经光纤导入分光与光电检测系统3中的分光器,分光器将接收系统2接收到的大气后向散射光分光成为该激光雷达所需的各种光信号,并由光电检测器检测,再送入数据处理系统5进行处理;数据处理系统5包括高速高精度数据采集卡、工控计算机、测试系统软件以及为得到测量结果所进行的数据分析与计算,最终求得大气的温度数值、水蒸气密度、大气气溶胶光学特性参数和非球形粒子后向散射光的退偏振比。
图2是实现本发明的分光与光电检测步骤的一装置实施例。分光与光电检测系统3由分光器和光电检测器组成,包括一个高光谱分辨率光栅Grating、一个边缘反射镜Edge_mirror、5个窄带干涉滤光片IF_1、IF_2、IF_3、IF_4、IF_5、十个聚焦透镜L_1-L_10、一个分光镜BS、一个偏振分光镜PBS以及光电检测部件PMT_1、PMT_2、PMT_3、PMT_4、PMT_5和PMT_6。大气的后向散射光由望远镜接收,经由光纤、聚光镜L_1照射到高光谱分辨率光栅Grating上,利用光栅的衍射效应,将雷达接收到的大气后向散射信号以及太阳光谱从空间上分离开,配合边缘反射镜Edge_mirror,将接收到的大气后向散射光谱中的纯转动拉曼光谱的反斯托克斯谱线与Mie-Rayleigh散射光谱以及太阳光谱从空间分开,让Mie-Rayleigh信号透过边缘反射镜,而转动拉曼信号由边缘反射镜反射。
反射的拉曼信号经透镜L_3入射到窄带干涉滤光片IF_1上,透过窄带干涉滤光片IF_1的信号再透过窄带干涉滤光片IF_2(IF_1与IF_2的中心波长相同,为λ1(如λ1=353.9nm),入射角相同,方向不同),经透镜L_4聚焦后由光电倍增管PMT_1接收,此为通道1;由窄带干涉滤光片IF_1反射的信号,射向窄带干涉滤光片IF_3,透过窄带干涉滤光片IF_3的信号再透过窄带干涉滤光片IF_4(IF_3与IF_4的中心波长相同,为λ2(如λ2=353.1nm),入射角相同,入射方向不同),经透镜L_5聚焦后由光电倍增管PMT_2接收,此为通道2;用这2个波长的窄带干涉滤光片,分离了中心波长为λ1和λ2的转动喇曼散射信号,同时对噪声信号信号进行2次高精度滤除,以满足白天及高密度气溶胶空间的温度测量要求.通道1和通道2的信号用于测量大气温度.
接收的大气回波信号在经由光栅Grating衍射后,还可分出波长为407.5nm的振动拉曼信号,经反射镜M反射后,射向窄带干涉滤光片IF_5(IF_5的中心波长为407.5nm),经透镜L_6聚焦后由光电倍增管PMT_3接收,此为通道3,用于测量水蒸气。
透过边缘反射镜Edge_mirror的Mie-Rayleigh散射信号用于测量气溶胶特性参数,这部分光束经透镜L_7准直后,射向分光镜BS,透过分光镜BS的光束经透镜L_8聚焦后,由光电倍增管PMT_4接收,此为通道4,用于测量气溶胶光学特性参数。由分光镜BS反射的光束经偏振分光镜PBS又分为2部分:由PBS反射的光束经透镜L_9聚焦后,由光电倍增管PMT_5接收,此为通道5;透过偏振分光镜PBS的光束经透镜L_10聚焦后,由光电倍增管PMT_6接收,此为通道6。通道5和通道6用于测量气溶胶粒子的偏振度。
窄带干涉滤光片IF_1和IF_2的中心波长相同,窄带干涉滤光片IF_3和IF_4的中心波长相同,并且两个波长应同时设在激光波长的同一侧,其中心波长位置要根据纯转动拉曼散射光谱函数作相应的理论计算,将窄带干涉滤光片IF_1和IF_2的中心波长设在转动拉曼信号随温度变化率为最小处(λ1),而将窄带干涉滤光片IF_3和IF_4的中心波长设在转动拉曼信号随温度变化率为最大处(λ2),其光束入射角在0~10度之间。图3提供了一种分光系统中边缘反射镜位置的确定示意图;利用高光谱分辨率光栅Grating对系统接收到的大气散射信号进行分光,光栅的谱分辨率为6pm,衍射光聚焦长度为300mm,取1级衍射级次。由光栅方程计算得到,经聚光镜后波长为λ1和λ2的衍射光与波长为λ0=354.67nm的衍射光的空间直线距离分别为0.66mm和1.37mm。可见,利用光栅衍射,可以将Raman波长λ1和λ2的散射信号和Mie-Rayleigh散射信号从空间上和光谱上分离开来,并同时也将日光光谱分开。
散射信号光谱可简化为高斯模型,先取光斑直径d0=0.30mm,温度T=300K,经光栅分光后,转动量子数分别为6和14的转动Raman散射信号强度与Mie-Rayleigh散射信号强度沿x轴的分布如图4所示。由图3可见,将边缘反射镜Edge_mirror的遮光位置设在0.43mm处,就可以使Mie-Rayleigh散射信号以及大部分太阳背景光透过,而将拉曼散射信号反射向窄带滤光片。
根据图3所示的光谱位置,将边缘反射镜Edge_mirror的遮光位置设在0.43mm处,对各散射信号从边缘反射镜Edge_mirror的遮光位置开始沿x轴进行积分,取J=6得到进入通道1的拉曼散射截面强度为1.199×10-34,取J=14得到进入通道2的拉曼散射截面强度为0.62×10-34,进入各通道内的剩余米散射截面强度和瑞利散射截面强度分别为0.3×10-34和1.128×10-34。改变光斑直径d0=0.25mm,Mie-Rayleigh散射信号的光谱分布也示于图中,积分后,得到进入各通道内剩余的瑞利散射截面强度为1.166×10-35,信号强度又降低了一个数量级。
由以上分析可见,利用光栅分光,并借助边缘反射镜,对Mie-Rayleigh散射信号的抑制率可以达到4个数量级,并且同时能将日光光谱分离并抑制。相对于只使用窄带干涉滤光片分光,这种分光方法光谱分辨率要高的多,更益于滤除Mie-Rayleigh散射信号及太阳背景光,同时实现了对拉曼信号的高效率获得。
散射信号经光栅分光以及边缘反射镜后,再分别通过窄带干涉滤光片IF_1、IF_2和IF_3、IF_4,又使Mie-Rayleigh散射信号被抑制3个数量级以上。至此,整个光路系统对Mie-Rayleigh散射信号的抑制率达到7个数量级以上。
将大气回波信号中的各个散射信号光谱成分分离、检测出来后,输入到数据处理系统5,利用其相互关系进行分析求解。
数据处理系统5包括多通道同步高速A/D采集卡,工控计算机,以及为求得大气温度、水蒸气密度(相对湿度)以及大气气溶胶消光系数、气溶胶光学厚度和非球形粒子后向散射光的退偏振比的垂直分布等参数所列的数据反演方法及其对应的应用软件。
分析处理程序包括根据光电检测部件PMT_1和PMT_2检测得到的拉曼散射信号,按拉曼散射激光雷达方程求出探测得到的两个信号的强度以及两个信号的强度比,结合系统温度探测的灵敏度及系统校正参数,求得大气的温度数值。通过利用气温探测的灵敏度及检测得到的拉曼散射信号的信噪比,继而求出系统的气温探测误差。
分析处理程序还包括根据光电检测部件PMT_3检测得到的水蒸气振动拉曼散射强度信号,利用水蒸气拉曼激光雷达方程和转动拉曼激光雷达方程,消除大气消光系数的影响,经反演求得水蒸气密度。
分析处理程序还包括根据光电检测部件PMT_4检测得到的米-瑞利散射信号,按米散射激光雷达方程,经反演求得大气气溶胶光学特性参数。
分析处理程序还包括根据光电检测部件PMT_5和PMT_6检测得到的与激光发射频率相同的散射信号的平行与垂直偏振分量信号,按偏振激光雷达方程,经反演求得非球形粒子后向散射光的退偏振比。
分析处理总程序如图5所示,首先将各通道进入的数据进行平滑处理和Y(z)、z2校正,接着分别进行大气温度数据分析、水蒸气数据分析、气溶胶数据分析、非球形粒子的数据分析。
具体来说,可以通过以下各式反演得到大气温度、水蒸气密度的垂直分布,以及大气气溶胶消光系数的垂直分布和大气气溶胶光学厚度和非球形气溶胶粒子的后向散射光的退偏振比。
首先根据下面的激光雷达方程,计算出进入雷达系统的各散射信号的功率P(z):
其中,K为光学系统效率,E0为激光脉冲能量,τ为激光脉冲间隔时间,Ar为望远镜受光面积,Y(z)为发送与接收器的光路重叠系数,z为探测高度,β(z)为高度z处的后向散射系数(是大气数密度N(z)和散射截面强度σ的函数),α(z)为高度z处的大气消光系数。方程中,除了大气消光系数α(z)和后向散射系数β(z)以外的参数都是系统提供的已知量。
(1)按图6所述的流程,进行大气温度的计算
由通道1和通道2检出的2个拉曼散射信号强度,分别由以下2个激光雷达方程表示:
其中,T为大气温度,J1和J2分别是两个通道探测的拉曼散射信号的转动量子数,σb_1(J1,T)和σb_2(J2,T)分别是温度为T时转动量子数J1的转动拉曼散射信号的散射截面强度和转动量子数J2的转动拉曼散射信号的散射截面强度,并且J1和J2为已知。
则由(2)式与(3)式比,得到通道1与通道2的散射信号强度比H(T,z)为:
经对理论上的σb_1(J1,T)和σb_2(J2,T)求解,并做曲线拟合,可以得到下列关系:
其中,A、B为系统常数,可以通过无线电探空仪对雷达系统标定得到。
则大气温度数值可以由(5)式导出:
(2)根据图7所示的流程进行气溶胶光学特性的计算
由通道4检出的米-瑞利散射信号强度P4(z),由以下激光雷达方程表示:
式中,β(z)=βm(z)+βa(z),α(z)=αm(z)+αa(z),
而βm(z)和βa(z)分别为大气分子和气溶胶粒子的后向散射系数,αm(z)和αa(z)分别为大气分子和气溶胶粒子的消光系数。
设激光雷达回波信号的距离平方修正函数为:
X(z)=P(z)·z2 (8)
如果事先已知某一高度Zc处大气总的消光系数α(Zc)或后向散射系数β(Zc),则Zc处以下各高度上的气溶胶粒子消光系数αa(z)或后向散射系数βa(z)分别为:
式中,αm(z)和βm(z)分别为大气分子的消光系数和后向散射系数,S为激光雷达比,S=αa/βa,它的取值与不同气溶胶的光学特性有关,为一常数。
上面(9)式和(10)式中,高度Z处的大气分子的消光系数αm(z)或后向散射系数βm(z)可以通过美国标准大气模型获得,如果测量高度较高,可以选取一段近乎不含气溶胶粒子的清洁大气层,在这段高度范围内,对激光雷达距离校正对数回波信号进行最小二乘拟合,该回归曲线斜率的一半即为大气分子消光系数,求出大气分子的消光系数高度分布模型,这种取值方法更加切合当时的大气状况。在测量高度较高时,可以选取近乎不含气溶胶粒子的清洁大气层所在高度作为Zc,则边界值α(Zc)或β(Zc)中就只含有大气分子的成分;同样,若测量高度不够高,可以选取一段较为均匀的大气层,对该高度范围内激光雷达距离校正对数回波信号进行最小二乘拟合,回归曲线斜率的一半即为高度Zc处的大气消光系数或后向散射系数的边界值α(Zc)或β(Zc)。
(3)根据图8的流程进行非球形粒子后向散射光的退偏振比计算
非球形粒子后向散射光的退偏振比可以用来大气非球形粒子的形态。
由通道5和通道6检出的2个大气回波信号的平行分量和垂直分量的散射信号强度,分别由以下2个激光雷达方程表示:
下标P和S分别表示与发射激光偏振方向平行和垂直的两个方向。
退偏振比δ(z)定义为
对于随机取向的卷云和大气气溶胶粒子,αP(z)=αS(z),则上式写为:
(4)根据图9所示的流程进行水蒸气密度的计算
由通道3检出的水蒸气振动拉曼散射信号强度PH(z),由以下激光雷达方程表示:
其中,βH(z)为高度z处的水蒸气后向散射系数。利用通道1的激光雷达方程(2)与该水蒸气通道3的激光雷达方程(15)之比,可以得到下式:
其中,N(z)可以由大气标准模型得到,σb_1(z)为已知,βH(z)用下式表示:
βH(z)=NH(z)·[dσH(π)/dΩ] (17)
其中,NH(z)是水蒸气密度,dσH(π)/dΩ是水蒸气拉曼散射的差分截面。因此,水蒸气密度可以由式(16)和式(17)导出:
NH(z)=N(z)·σb_1(z)/SH(z)·[dσH(π)/dΩ]。 (18)
利用Tetens经验公式,可以计算出大气温度t状态下的水面饱和水汽压e(t);
e(t)=610.8×exp[17.27t/(t+237.3] (19)
式中,e(t)的单位:Pa,温度t的单位:℃。
利用饱和水汽压e(t),可以计算出饱和水汽密度SH(t),其单位为g/m3,
SH(t)=2.17×e(t)/(t+273.15). (20)
利用式(15)和(17),可以求得相对湿度UH(z)(单位:%)的高度分布:
UH(z)=NH(z)×100/SH[t(z)]. (21)
式中t(z)为转动拉曼方式测得的大气温度高度分布。
实施例
图9表示大气温度分布的一个测量实例。它通过利用由通道1和通道2检出的2个拉曼散射信号强度,结合公式(2)-(6)算法,求得的一个大气温度高度分布结果。
图10表示相对湿度的分布实例,它通过利用通道1(或通道2)检测出的一个转动拉曼散射信号和通道3检测出的水汽的振动拉曼散射信号,结合公式(2),(15)-(18)算出水汽密度分布,然后利用求得的大气温度(t)分布,结合公式(19)-(21)即可求得的一个大气相对的湿度分布。
Claims (2)
1.一种气象与大气环境参数的测定方法,将大气中激光脉冲与大气分子、粒子相互作用后产生的后向散射光接收下来,经分光成为激光雷达所需的各种光信号,再经光电检测后,根据设定的程序进行数据的分析处理,从而实现大气温度、水蒸气密度垂直分布的气象参数测量,以及大气气溶胶消光系数及散射系数、气溶胶光学厚度、大气能见度和非球形粒子后向散射光的退偏振比垂直分布的大气环境参数的测量,其特征在于,具体包括以下步骤:
a.向大气中发射脉冲激光束;
b.接收上述脉冲激光与大气中的分子和粒子相互作用后产生的后向散射光;
c.对上述大气后向散射光进行分光,将后向散射光中的转动拉曼谱线、水汽分子的振动拉曼谱线、米-瑞利散射谱线、太阳背景光相互进行分离,分别得到2个转动拉曼信号,用于温度测量,得到的水蒸气振动拉曼信号,用于水蒸气密度测量,得到的米-瑞利散射信号,用于大气气溶胶光学特性的测量,得到的大气后向散射光信号的偏振特性,即与激光发射频率相同的散射信号的平行与垂直偏振分量,用于大气中非球形气溶胶形态的测量;
将大气后向散射光先通过一高光谱分辨率光栅Grating后再经透镜聚光,使不同波长的水蒸气振动拉曼散射光、转动拉曼散射光、米-瑞利散射光以及太阳背景光光谱信号在透镜焦平面上按波长呈有序空间排列,通过设定的边缘反射镜Edge_mirror让米-瑞利散射光以及大部分太阳背景光通过,而将转动拉曼散射光反射,使反射的转动拉曼散射光的一部分透过两个窄带干涉滤光片IF_1和IF_2后,以中心波长为λ1的光被有选择性地透过,并由光电检测部件PMT_1检测,作为通道1;反射的转动拉曼散射光的另一部分经窄带干涉滤光片IF_1表面反射后透过窄带干涉滤光片IF_3和IF_4,以中心波长为λ2的光被有选择性地透过,并由光电检测部件PMT_2检测,作为通道2;使光栅Grating衍射分离出的水蒸气振动拉曼散射信号透过窄带干涉滤光片IF_5,并由光电检测部件PMT_3检测,作为通道3;使透过边缘反射镜Edge_mirror的米-瑞利散射信号通过分光镜BS被分为两路:透过的一路由光电检测部件PMT_4检测,作为通道4,使反射的一路再通过偏振分光镜PBS又被分为两路,一路由光电检测部件PMT_5检测,作为通道5,另一路由光电检测部件PMT_6检测,作为通道6。
d.将上述分离后的各种散射光信号变为电信号接收;
e.将上述处理后的各种信号按照设定的程序分别进行处理,最终得到大气温度、水蒸气密度的垂直分布,以及大气气溶胶消光系数的垂直分布和大气气溶胶光学厚度和非球形气溶胶粒子的后向散射光的退偏振比。
2.按照权利要求1所述的测定方法,其特征在于,所述步骤e中将处理后的各种信号按照设定的程序分别进行处理,是按照以下方式进行:
首先根据下面的激光雷达方程,计算出处理后的各散射信号的功率P(z):
其中,K为光学系统效率,E0为激光脉冲能量,τ为激光脉冲间隔时间,Ar为望远镜受光面积,Y(z)为发送与接收器的光路重叠系数,z为探测高度,β(z)为高度z处的后向散射系数,α(z)为高度z处的大气消光系数,光学系统效率K,激光脉冲能量EO,激光脉冲间隔时间τ,望远镜受光面积Ar,发送与接收器的光路重叠系数Y(z),探测高度z根据系统预先输入;
再分别进行:
1)大气温度的计算
根据分光检测得到的2个转动拉曼信号,分别由以下2个激光雷达方程表示:
其中,T为大气温度,J1和J2分别是两个拉曼散射信号的转动量子数,σb_1(J1,T)是温度为T、转动量子数为J1的转动拉曼散射信号的散射截面强度,σb_2(J2,T)是温度为T、转动量子数为J2的转动拉曼散射信号的散射截面强度,N(z)是高度z处的大气分子数密度。
则由(2)式与(3)式比,得到通道1与通道2的散射信号强度比H(T,z)为:
经对σb_1(J1,T)和σb_2(J2,T)的求解,并做曲线拟合,得到下列关系:
其中,A、B为系统常数,通过无线电探空仪对雷达系统标定得到,根据下式
得出大气温度数值;
2)气溶胶光学特性的计算
由通道4检出的米-瑞利散射信号强度P4(z),由以下激光雷达方程表示:
式中,β(z)=βm(z)+βa(z),α(z)=αm(z)+αa(z),
而βm(z)和βa(z)分别为大气分子和气溶胶粒子的后向散射系数,αm(z)和αa(z)分别为大气分子和气溶胶粒子的消光系数,
设激光雷达回波信号的距离平方修正函数为:
X(z)=P(z)·z2 (8)
根据已知某一高度Zc处大气总的消光系数α(Zc)或后向散射系数β(Zc),则Zc处以下各高度上的气溶胶粒子消光系数αa(z)或后向散射系数βa(z)分别为:
式中,αm(z)和βm(z)分别为大气分子的消光系数和后向散射系数,S为激光雷达比,S=αa/βa;
3)非球形粒子后向散射光的退偏振比计算
由通道5和通道6检出的2个大气回波信号的平行分量和垂直分量的散射信号强度,分别由以下2个激光雷达方程表示:
下标P和S分别表示与发射激光偏振方向平行和垂直的两个方向,
退偏振比δ(z)定义为
对于随机取向的卷云和大气气溶胶粒子,αP(z)=αS(z),则上式写为:
4)水蒸气密度的计算
由通道3检出的水蒸气振动拉曼散射信号强度PH(z),由以下激光雷达方程表示:
其中,βH(z)为高度z处的水蒸气后向散射系数,
利用通道1的激光雷达方程(2)与该水蒸气通道3的激光雷达方程(15)之比,可以得到下式:
其中,N(z)可以由大气标准模型得到,σb_1(z)为已知,βH(z)用下式表示:
βH(z)=NH(z)·[dσH(π)/dΩ] (17)
其中,NH(z)是水蒸气密度,dσH(π)/dΩ是水蒸气拉曼散射的差分截面,水蒸气密度由式(16)和式(17)导出:
NH(z)=N(z)·σb_1(z)/SH(z)·[dσH(π)/dΩ]。 (18)
利用Tetens经验公式,计算出大气温度t状态下的水面饱和水汽压e(t);
e(t)=610.8×exp[17.27t/(t+237.3] (19)
式中,e(t)的单位:Pa,温度t的单位:℃;
利用饱和水汽压e(t),可以计算出饱和水汽密度SH(t),其单位为g/m3,
SH(t)=2.17×e(t)/(t+273.15) (20)
利用式(15)和(17),可以求得相对湿度UH(z)(单位:%)的高度分布:
UH(z)=NH(z)×100/SH[t(z)] (21)
式中t(z)为转动拉曼方式测得的大气温度高度分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200610105194A CN101004453B (zh) | 2006-12-20 | 2006-12-20 | 一种气象与大气环境参数的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200610105194A CN101004453B (zh) | 2006-12-20 | 2006-12-20 | 一种气象与大气环境参数的测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101004453A CN101004453A (zh) | 2007-07-25 |
CN101004453B true CN101004453B (zh) | 2010-05-12 |
Family
ID=38703729
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200610105194A Expired - Fee Related CN101004453B (zh) | 2006-12-20 | 2006-12-20 | 一种气象与大气环境参数的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101004453B (zh) |
Families Citing this family (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329173B (zh) * | 2008-07-07 | 2010-06-02 | 武汉大学 | 一种浑浊水体大气校正方法 |
CN101672768B (zh) * | 2008-09-11 | 2011-08-10 | 中国海洋大学 | 海上浓雾条件下大气水平能见度场的获取方法 |
CN101923162B (zh) * | 2009-06-09 | 2012-08-22 | 中国科学院安徽光学精密机械研究所 | 拉曼激光雷达标定装置及其标定方法 |
CN101915757A (zh) * | 2010-08-10 | 2010-12-15 | 北京路源光科技有限公司 | 简易型拉曼远距离爆炸物的快速检测装置 |
CN101963665B (zh) * | 2010-08-23 | 2013-11-06 | 西安理工大学 | 激光雷达几何重叠因子自动调整方法 |
CN102200508A (zh) * | 2011-03-23 | 2011-09-28 | 南京信息工程大学 | 消除背景光和热噪声对大气消光系数测量精度影响的方法 |
CN102539383B (zh) * | 2012-01-09 | 2013-07-10 | 北京大学 | 散射式能见度仪中能见度测量方法 |
CN103852449A (zh) * | 2012-12-05 | 2014-06-11 | 中国辐射防护研究院 | 一种高灵敏度与宽测量范围气溶胶浓度的光学检测方法 |
CN103115872B (zh) * | 2012-12-18 | 2015-02-25 | 中国人民解放军63655部队 | 一种多波长大气消光系数高度分布数据反演方法 |
CN103234877B (zh) * | 2013-04-10 | 2015-03-18 | 中国科学院合肥物质科学研究院 | 一种大气颗粒物粒径谱时空分布激光雷达数据反演方法 |
CN103344611B (zh) * | 2013-07-16 | 2015-07-01 | 中国人民解放军陆军军官学院 | 基于ccd成像技术的侧向激光雷达测量气溶胶参数的方法 |
CN103675791B (zh) * | 2013-12-05 | 2015-12-30 | 北京师范大学 | 基于数值分布均衡化的米散射激光雷达云识别方法 |
CN103942439B (zh) * | 2014-04-24 | 2017-05-17 | 中国科学院遥感与数字地球研究所 | 一种基于气象观测数据的可吸入颗粒物浓度估算方法 |
CN104007088B (zh) * | 2014-06-16 | 2016-03-30 | 中国人民解放军陆军军官学院 | 后向散射激光雷达几何因子的测量方法 |
CN104034700B (zh) * | 2014-06-21 | 2016-09-21 | 中国科学院合肥物质科学研究院 | 一种大气传输激光透过率的测量方法 |
CN104198439A (zh) * | 2014-09-23 | 2014-12-10 | 南昌航空大学 | 一种检测高纯气体中气溶胶的方法及装置 |
CN104359810B (zh) * | 2014-10-30 | 2017-10-27 | 儒亚科技(北京)有限公司 | 一种远程无线pm2.5浓度检测装置 |
CN104655576B (zh) * | 2015-02-09 | 2017-06-13 | 上海安允科技有限公司 | 基于后向散射回光的自由状态气体参数远程测量方法 |
CN104767483B (zh) * | 2015-03-31 | 2017-08-25 | 国家电网公司 | 一种光热特性耦合光伏电池温度预测方法 |
CN104808194A (zh) * | 2015-05-07 | 2015-07-29 | 中国科学院合肥物质科学研究院 | 一种纯转动拉曼激光雷达系统常数的标定方法 |
CN104880711B (zh) * | 2015-05-18 | 2016-08-17 | 中国科学院合肥物质科学研究院 | 单波长四拉曼激光雷达探测系统及探测方法 |
CN105115862A (zh) * | 2015-07-02 | 2015-12-02 | 南京信息工程大学 | 一种云粒子探测方法及探测器 |
CN105204014A (zh) * | 2015-09-08 | 2015-12-30 | 中国科学院合肥物质科学研究院 | 一种双波段纯转动拉曼激光雷达系统常数的自标定方法 |
CN105487083B (zh) * | 2015-12-01 | 2017-11-21 | 中国科学院上海技术物理研究所 | 一种多光谱影像气溶胶像元气溶胶厚度的检测与提取方法 |
CN105676231B (zh) * | 2016-01-25 | 2018-04-27 | 西安理工大学 | 一种基于转动拉曼激光雷达反演大气温度的方法 |
CN105891064B (zh) * | 2016-04-05 | 2018-12-28 | 山东大学 | 一种非球形气溶胶粒子混合率的检测方法及其装置 |
CN105911534B (zh) * | 2016-04-12 | 2018-11-20 | 中国科学院上海光学精密机械研究所 | 透过率随接收视场角变化的光学滤波装置 |
CN105675466B (zh) * | 2016-04-12 | 2018-09-11 | 北京理工大学 | 一种大气相对湿度协同颗粒物消光的在线监测装置 |
CN106226842B (zh) * | 2016-07-13 | 2018-11-13 | 成都信息工程大学 | 一种城市下垫面气溶胶对雷电过程影响机理的测试方法 |
CN106249252B (zh) * | 2016-07-22 | 2019-02-01 | 西安理工大学 | 探测云中过冷水的机载近红外激光雷达系统及反演方法 |
CN106054158A (zh) * | 2016-08-09 | 2016-10-26 | 北方民族大学 | 一种探测拉曼激光雷达光路系统 |
CN106443709A (zh) * | 2016-09-14 | 2017-02-22 | 中国科学技术大学 | 基于时分复用的单探测器全光纤偏振激光雷达 |
CN106569226B (zh) * | 2016-10-28 | 2019-02-05 | 兰州大学 | 一种有云情况下利用激光雷达资料反演边界层高度的方法 |
CN106814371B (zh) * | 2017-01-20 | 2020-06-09 | 武汉大学 | 一种测量大气温度和水汽以及气溶胶的激光雷达系统 |
CN107505630B (zh) * | 2017-08-10 | 2020-11-10 | 成都信息工程大学 | 基于气溶胶特性和成云机制提高天气预报准确性的方法 |
CN107688187B (zh) * | 2017-08-16 | 2021-01-08 | 南京红露麟激光雷达科技有限公司 | 基于空间波长编码的目标探测激光雷达 |
CN107796742B (zh) * | 2017-09-20 | 2020-05-15 | 中国计量科学研究院 | 一种校准大气组分浓度探测激光雷达的装置 |
CN107991282B (zh) * | 2017-11-30 | 2020-05-26 | 青岛大学 | 一种利用卫星分析大气Ring效应的方法及系统 |
CN108957474B (zh) * | 2018-06-12 | 2022-07-15 | 西安理工大学 | 用于检测粒子形态的全偏振激光雷达系统及其检测方法 |
WO2019101247A2 (zh) * | 2018-07-20 | 2019-05-31 | 大舜激光(黄山)科技有限公司 | 激光海面能见度监测仪以及探测海雾的方法 |
CN109633654B (zh) * | 2018-12-04 | 2020-07-31 | 上海无线电设备研究所 | 一种用于太赫兹雷达的卷云微物理参数计算方法 |
CN109782019B (zh) * | 2019-03-15 | 2020-04-24 | 中国科学技术大学 | 大气污染物二维运动速度测量方法及装置 |
CN110006848B (zh) * | 2019-05-05 | 2022-02-01 | 北京是卓科技有限公司 | 一种获取气溶胶消光系数的方法和装置 |
CN110006849B (zh) * | 2019-05-05 | 2021-08-31 | 北京是卓科技有限公司 | 一种获取气溶胶消光系数的方法和装置 |
CN110850392B (zh) * | 2019-10-31 | 2023-06-30 | 无锡中科光电技术有限公司 | 激光雷达退偏增益比的校正方法 |
CN110850443B (zh) * | 2019-11-11 | 2023-03-03 | 西北工业大学 | 激光雷达温湿度数据阶梯分析处理方法 |
CN111736176B (zh) * | 2020-06-23 | 2024-05-31 | 中国人民解放军陆军炮兵防空兵学院 | 一种基于激光雷达的大气环境数据的共享方法 |
CN113378377B (zh) * | 2021-06-08 | 2024-04-05 | 西安理工大学 | 一种建立随时空变化大气密度模型的方法 |
CN114049755B (zh) * | 2021-12-13 | 2024-07-23 | 夏若凡 | 一种环境测试系统和测试方法 |
CN114441505B (zh) * | 2022-03-17 | 2023-08-18 | 中国工程物理研究院机械制造工艺研究所 | 一种用于拉曼探头的水汽原位标定装置、标定方法及应用 |
CN114994711B (zh) * | 2022-08-05 | 2022-10-04 | 南京信息工程大学 | 一种基于菲索干涉仪的激光雷达 |
CN118393509B (zh) * | 2024-06-28 | 2024-08-20 | 中电建(洛阳)绿色建筑科技有限公司 | 一种复杂环境下的大气校正方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1541339A (zh) * | 2002-02-21 | 2004-10-27 | Ӣ�뾫����ʽ���� | 气象观测激光雷达系统 |
-
2006
- 2006-12-20 CN CN200610105194A patent/CN101004453B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1541339A (zh) * | 2002-02-21 | 2004-10-27 | Ӣ�뾫����ʽ���� | 气象观测激光雷达系统 |
Non-Patent Citations (7)
Title |
---|
Dengxin Hua, Cong Dai, Kobayashi T.Daytime Temperature and Humidity ProfileMeasurementwithUV Rayleigh-Mie Raman Lidar.Lasers and Electro-Optics,2005. CLEO/Pacific Rim 2005..2005,1263-1264. |
Dengxin Hua,Cong Dai,Kobayashi T.Daytime Temperature and Humidity ProfileMeasurementwithUV Rayleigh-Mie Raman Lidar.Lasers and Electro-Optics,2005. CLEO/Pacific Rim 2005.2005,1263-1264. * |
吴永华,胡欢陵,胡顺星,周军,岳古明,戚福弟,李琛.瑞利-拉曼散射激光雷达探测大气温度分布.中国激光31 7.2004,31(7),851-856. |
吴永华,胡欢陵,胡顺星,周军,岳古明,戚福弟,李琛.瑞利-拉曼散射激光雷达探测大气温度分布.中国激光31 7.2004,31(7),851-856. * |
杨洋,关荣华.用于大气环境监测的YAG激光雷达技术的现状与动向.红外与激光工程30 6.2001,30(6),466-469. |
杨洋,关荣华.用于大气环境监测的YAG激光雷达技术的现状与动向.红外与激光工程30 6.2001,30(6),466-469. * |
王建宇,洪光烈.激光主动遥感技术及其应用.激光与红外36.2006,36742-748. * |
Also Published As
Publication number | Publication date |
---|---|
CN101004453A (zh) | 2007-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101004453B (zh) | 一种气象与大气环境参数的测定方法 | |
CN100543495C (zh) | 气象与大气环境观测拉曼散射激光雷达系统 | |
CN105738916B (zh) | 高光谱偏振大气探测激光雷达系统及控制方法 | |
CN101833089B (zh) | 多普勒测风激光雷达灵敏度标定系统及其标定方法 | |
US4123160A (en) | Method and apparatus for remotely measuring sub-surface water temperatures | |
Varotsos | Airborne measurements of aerosol, ozone, and solar ultraviolet irradiance in the troposphere | |
CN108303706B (zh) | 一种气溶胶光学参数探测方法及高光谱激光雷达探测系统 | |
CN110045391B (zh) | 一种用于气溶胶尺度谱测量的高光谱激光雷达系统 | |
CN110967704A (zh) | 多波长测量大气二氧化碳浓度及气溶胶垂直廓线的激光雷达系统装置 | |
CN207882443U (zh) | 一种高光谱激光雷达探测系统 | |
Titos et al. | Spatiotemporal evolution of a severe winter dust event in the western Mediterranean: Aerosol optical and physical properties | |
CN102944503A (zh) | 基于太阳光度计和激光雷达的pm2.5质量浓度值自动反演算法 | |
CN103344964B (zh) | 一种遥感区分灰霾类别的激光雷达装置 | |
CN110441792B (zh) | 同时测风测温的瑞利散射激光雷达系统及相关校准方法 | |
CN206638826U (zh) | 一种用于液态水云测量的多视场激光雷达探测系统 | |
WO2011066808A1 (zh) | 实时定标高光谱分辨率激光雷达装置 | |
KR101012279B1 (ko) | 라이다를 이용한 황사 판독방법 | |
CN102879835A (zh) | 一种激光降水天气现象的测量方法和激光降水天气现象仪 | |
CN107064957A (zh) | 一种用于液态水云测量的多视场激光雷达探测系统及方法 | |
Jayaraman et al. | Spatial variations in aerosol characteristics and regional radiative forcing over India: Measurements and modeling of 2004 road campaign experiment | |
CN103399363A (zh) | 基于光衰减和散射原理的天气现象在线观测装置及方法 | |
CN106093915B (zh) | 一种新型拉曼测温激光雷达的分光系统 | |
CN204086139U (zh) | 用于区分雾和霾的双波长激光探测系统 | |
CN106483531B (zh) | 大气拉曼-瑞利散射测温激光雷达及反演方法 | |
Yufeng et al. | A UV multifunctional Raman lidar system for the observation and analysis of atmospheric temperature, humidity, aerosols and their conveying characteristics over Xi'an |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100512 Termination date: 20121220 |