CN106093964B - 一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 - Google Patents
一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 Download PDFInfo
- Publication number
- CN106093964B CN106093964B CN201610676685.0A CN201610676685A CN106093964B CN 106093964 B CN106093964 B CN 106093964B CN 201610676685 A CN201610676685 A CN 201610676685A CN 106093964 B CN106093964 B CN 106093964B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- scattering signal
- mfrac
- air
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
-
- 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
Abstract
本发明提供了一种测量大气风场和温度的钠荧光多普勒激光雷达及方法,本发明通过将钠荧光多普勒激光雷达接收到的大气散射信号的一部分分离出来,透过碘分子吸收光谱仪,利用碘分子在589nm的吸收谱线检测大气瑞利散射信号的光强,进而测量得到大气瑞利散射信号所对应的大气风速和温度信息,使钠荧光多普勒激光雷达的探测高度范围由大气中间层顶区域扩展至平流层、对流层,扩展了钠荧光多普勒激光雷达的大气探测范围。
Description
技术领域
本发明涉及大气探测技术领域,具体涉及一种测量大气风场和温度的钠荧光多普勒激光雷达及方法。
背景技术
在利用激光进行大气风场和温度的测量中,589nm激光已经被应用到钠荧光多普勒激光雷达中。钠荧光多普勒激光雷达发射的窄线宽589nm激光使中间层顶区域(约75-105km)的钠原子产生共振荧光,通过接收并分析共振荧光的多普勒频移和展宽获得大气风场和温度。
1990年美国科罗拉多州立大学和伊利诺伊大学合作建立了国际首台钠荧光多普勒激光雷达,进行中间层顶区域(约75-105km)大气温度的探测。2002年美国科罗拉多州立大学的钠荧光多普勒激光雷达升级为两个发射光束,实现了大气风场与温度的同时观测。2011年,中国科学院空间科学与应用研究中心研制完成了中国首台钠荧光多普勒激光雷达,实现了对中间层顶区域大气风场和温度的测量。随后,2012年,中国科学技术大学也完成了钠荧光多普勒激光雷达的研制。钠荧光多普勒激光雷达已成为探测中间层顶区域大气风场和温度的重要手段,对中高层大气的科学研究具有重要的意义。
现有的钠荧光多普勒激光雷达探测得到的信号主要包括大气瑞利散射信号和钠荧光散射信号。大气瑞利散射信号主要来自于对流层与平流层区域;钠荧光散射信号主要来自于中间层顶区域。目前,主要使用钠荧光散射信号进行中间层顶区域的大气风场和温度测量,而对流层、平流层的大气散射信号没有得到利用。2009年,科罗拉多大学开展了利用对流层、平流层大气散射信号探测大气风场的研究工作,使用一钠原子磁光滤波器检测589nm激光散射信号的多普勒频移和展宽。由于钠原子磁光滤波器的结构比较复杂,其稳定性受钠原子泡温度、磁场、光路等多种因素制约,从而限制了其推广使用。
发明内容
本发明的目的在于,为了更好地解决利用钠荧光多普勒激光雷达瑞利散射信号测量大气风场和温度的技术问题,使钠荧光多普勒激光雷达的探测高度范围由大气的中间层扩展至平流层和对流层,本发明提供一种利用589nm激光大气瑞利散射信号测量大气风场和温度的钠荧光多普勒激光雷达及方法。该方法不仅容易实现,而且比较稳定可靠。
为实现上述目的,本发明提供的一种测量大气风场和温度的钠荧光多普勒激光雷达,包括:激光发射器、大气散射信号接收器、钠荧光散射信号测量器、大气瑞利散射信号测量器、信号采集器;
所述的激光发射器用于向大气发射不同频率的589nm激光,所述的大气散射信号接收器接收由大气中间层顶区域产生的钠荧光散射信号,以及由对流层、平流层产生的大气瑞利散射信号,并将其输出端分两路分别连接钠荧光散射信号测量器和大气瑞利散射信号测量器,所述的钠荧光散射信号测量器接收不同频率的钠荧光散射信号,并通过分析获得中间层顶区域的温度值和风速值;所述的大气瑞利散射信号测量器接收不同频率的大气瑞利散射信号,并测量各频率的大气瑞利散射信号透过该大气瑞利散射信号测量器中设有的碘分子吸收光谱仪的光强后,计算获得对流层、平流层的温度值和风速值,所述的信号采集器用于采集并判读钠荧光散射信号测量器和大气瑞利散射信号测量器输出的电信号,提取光强数据。
作为上述技术方案的进一步改进,所述的大气瑞利散射信号测量器包括:准直镜、背景光抑制器、碘分子吸收光谱仪、耦合镜和光电探测器;所述的准直镜通过光纤与大气散射信号接收器的输出端连接,并依次连接背景光抑制器、碘分子吸收光谱仪、耦合镜和光电探测器;所述的大气瑞利散射信号依次经准直镜准直操作、背景光抑制器滤除天空背景光及透过碘分子吸收光谱仪后,由耦合镜将大气瑞利散射信号汇聚到光电探测器,并通过该光电探测器测量透过碘分子吸收光谱仪后的信号光强。
作为上述技术方案的进一步改进,所述的大气散射信号接收器包括两个信号接收通道,所述的两个信号接收通道以两个不同的方向接收钠荧光散射信号和大气瑞利散射信号,且每一个信号接收通道的输出端均分两路连接钠荧光散射信号测量器和大气瑞利散射信号测量器,使得钠荧光散射信号测量器以两路接收来自两个不同方向上的钠荧光散射信号,用于测量大气中间层顶区域的温度值和水平风场,使得大气瑞利散射信号测量器以两路接收来自两个不同方向上的大气瑞利散射信号,用于测量对流层、平流层的温度值和水平风场。
本发明还提供了一种测量大气风场和温度的方法,该方法包括:
步骤1)利用激光发射器依次向大气发射三个不同频率的589nm激光;
步骤2)通过大气散射信号接收器依次接收由大气中间层顶区域产生的三个不同频率的钠荧光散射信号后,执行步骤3),通过大气散射信号接收器依次接收由对流层、平流层产生的三个不同频率的大气瑞利散射信号后,执行步骤4);
步骤3)利用钠荧光散射信号测量器测量得到的三个不同频率的钠荧光散射信号所对应的光强N0H(V,T)、NMH(V,T)和NPH(V,T),定义激光雷达风速测量比和激光雷达温度测量比分别为:
则有测量得到的大气中间层顶区域的风速值和温度值分别为:
其中,SWH和STH分别表示风速测量灵敏度和温度测量灵敏度;
所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
步骤4)利用光电探测器测量得到透过碘分子吸收光谱仪的三个不同频率的大气瑞利散射信号所对应的光强N0L(V,T)、NML(V,T)和NPL(V,T),定义激光雷达风速测量比和激光雷达温度测量比分别为:
则有测量得到的大气对流层、平流层的风速值和温度值分别为:
其中,SWL和STL分别表示风速测量灵敏度和温度测量灵敏度;
所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
作为上述技术方案的进一步改进,所述激光发射器向大气发射三个589nm激光的频率分别位于钠原子D2a饱和荧光光谱处f0、f0+630MHz和f0-630MHz。
本发明的一种测量大气风场和温度的钠荧光多普勒激光雷达及方法优点在于:
本发明通过将钠荧光多普勒激光雷达接收到的大气散射信号的一部分分离出来,透过碘分子吸收光谱仪,利用碘分子在589nm的吸收谱线检测大气瑞利散射信号的光强,进而测量得到大气瑞利散射信号所对应的大气风速和温度信息,使钠荧光多普勒激光雷达的探测高度范围由大气中间层顶区域扩展至平流层、对流层,扩展了钠荧光多普勒激光雷达的大气探测范围。
附图说明
图1为碘分子吸收谱线与大气瑞利散射信号透过率曲线的关系视图。
图2为本发明实施例中的一种测量大气风场和温度的装置结构示意图。
图3为本发明实施例中的大气瑞利散射信号测量器结构示意图。
附图标记
1、碘分子吸收谱线
2、碘分子吸收光谱仪对200K、0.267大气压状态下的589nm大气瑞利散射光谱的透过曲线
3、碘分子吸收光谱仪对230K、0.267大气压状态下的589nm大气瑞利散射光谱的透过曲线
4、碘分子吸收光谱仪对260K、0.267大气压状态下的589nm大气瑞利散射光谱的透过曲线
5、钠荧光多普勒激光雷达发射激光频率F0
6、钠荧光多普勒激光雷达发射激光频率FM
7、钠荧光多普勒激光雷达发射激光频率FP
8、激光发射器 9、大气散射信号接收器
10、钠荧光散射信号测量器 11、信号采集器
12、大气瑞利散射信号测量器 13、光纤
14、准直镜 15、背景光抑制器
16、碘分子吸收光谱仪 17、耦合镜
18、光电探测器
具体实施方式
下面结合附图和实施例对本发明所述的一种测量大气风场和温度的钠荧光多普勒激光雷达及方法进行详细说明。
研究发现,在钠荧光多普勒激光雷达的激光发射589nm波长处,碘分子存在吸收光谱线,可以用于检测大气散射信号的多普勒频移和展宽。将一定量的碘分子封装在透明的泡体内,经过加热等制作成碘分子吸收光谱仪。由于多普勒展宽和压力展宽等因素,碘分子吸收光谱线具有一定的宽度,如图1中示出的碘分子吸收谱线1形成“U”字形状。它对于不同频率和宽度的大气散射信号的透过率也不同,因此可以通过测量大气散射信号对于碘分子吸收光谱仪的透过率而获得大气散射信号的多普勒频移和展宽,进而得到大气的风场和温度。
具体做法是:钠荧光多普勒激光雷达接收系统接收到的大气散射信号被分为两部分,一部分保持钠荧光多普勒激光雷达原来的光路结构,不经过碘分子吸收光谱仪,另一部分经过碘分子吸收光谱仪。这两部分信号分别用于测量中间层顶区域和对流层、平流层区域的大气风速和温度。其中大气散射信号包括大气瑞利散射信号和钠荧光散射信号,发射激光为短脉冲激光,不同时刻接收到的散射信号即为不同高度的大气瑞利散射信号或钠荧光散射信号。
如图2所示,本发明提供的一种测量大气风场和温度的钠荧光多普勒激光雷达,包括:激光发射器8、大气散射信号接收器9、钠荧光散射信号测量器10、大气瑞利散射信号测量器12、信号采集器11;
所述的激光发射器8用于向大气发射不同频率的589nm激光,所述的大气散射信号接收器9接收由大气中间层顶区域产生的钠荧光散射信号,以及由对流层、平流层产生的大气瑞利散射信号,并将其收集到的大气散射信号分为两部分,一部分送给钠荧光散射信号测量器10,用于探测中间层大气风场和温度,另一部分送给大气瑞利散射信号测量器12,用于测量平流层、对流层大气风场和温度。在本实施例中,所述的大气散射信号接收器9包括两个信号接收通道,所述的两个信号接收通道以两个不同的方向接收钠荧光散射信号和大气瑞利散射信号,且每一个信号接收通道的输出端均分两路连接钠荧光散射信号测量器10和大气瑞利散射信号测量器12,使得钠荧光散射信号测量器10以两路接收来自两个不同方向上的钠荧光散射信号,用于测量大气中间层顶区域的温度值和水平风场,使得大气瑞利散射信号测量器12以两路接收来自两个不同方向上的大气瑞利散射信号,用于测量对流层、平流层的温度值和水平风场。
所述的钠荧光散射信号测量器10接收各方向上不同频率的钠荧光散射信号,并通过分析获得中间层顶区域的温度值和风速值,所述的大气瑞利散射信号测量器12接收各方向上不同频率的大气瑞利散射信号,并测量各频率的大气瑞利散射信号透过碘分子吸收光谱仪16的光强后,计算获得对流层、平流层的温度值和风速值。所述的信号采集器11用于采集并判读钠荧光散射信号测量器10和大气瑞利散射信号测量器12输出的电信号,由信号采集器11将电信号转化为数字信号,从而提取光强数据。
基于上述结构的钠荧光多普勒激光雷达,如图3所示,在本实施例中,所述的大气瑞利散射信号测量器包括:光纤13、准直镜14、背景光抑制器15、碘分子吸收光谱仪16、耦合镜17和光电探测器18。所述的准直镜14通过光纤13与大气散射信号接收器9的输出端连接,并依次连接背景光抑制器15、碘分子吸收光谱仪16、耦合镜17和光电探测器18;光纤13输出的光信号经准直镜14准直,然后经过背景光抑制器15滤除天空背景光,随后经过碘分子吸收光谱仪16进行多普勒信号频移和展宽的检测,多普勒频移会导致光信号中心频率改变,展宽则是指光信号的频率分布展宽,不同中心频率以及不同宽度的光信号经过碘分子吸收光谱仪之后的透过率不同,而大气瑞利散射光经过碘分子吸收光谱仪后其透射光强就包涵了大气瑞利散射光的频率信息。经过碘分子吸收光谱仪后由耦合镜17将光信号汇聚到光电探测器18的光阴极,并通过该光电探测器18测量透过碘分子吸收光谱仪16后的信号光强。此时可采用光电探测器18内置的光电倍增管将光信号转换为电信号。所述的碘分子吸收光谱仪为密封有碘分子的玻璃泡,并进行恒温控制。
不同温度、风场状态下的大气瑞利散射光信号经过碘分子吸收光谱仪的透过率也不同。碘分子吸收光谱仪的透过率曲线如图中1所示。不同温度下的大气瑞利散射光谱经过碘分子吸收光谱仪的透过率曲线如图中示出的曲线2、3、4所示。
利用上述原理,本发明还提供了一种测量大气风场和温度的方法,包括:
步骤1)利用激光发射器8依次向大气发射三个不同频率的589nm激光;
光信号中心频率位于透过率曲线高的地方则碘分子对其吸收较弱,信号较强;反之信号较弱。图1中示出的曲线5、6、7是钠荧光多普勒激光雷达所发射的频率分别为F0、FM、FP位置的激光。其中,F0位于钠原子D2a饱和荧光光谱处,FM=F0-630MHz,FP=F0+630MHz。三个不同频率的激光交替发射,在时间上是分离的。
步骤2)通过大气散射信号接收器9依次接收由大气中间层顶区域产生的三个不同频率的钠荧光散射信号后,执行步骤3),通过大气散射信号接收器9依次接收由对流层、平流层产生的三个不同频率的大气瑞利散射信号后,执行步骤4);
步骤3)在上述三个发射激光频率F0、FM、FP工作状态下,利用钠荧光散射信号测量器10测量得到的三个不同频率钠荧光散射信号的光强可以表示为:N0H(V,T)、NMH(V,T)、NPH(V,T)。它们与大气的风速V和温度T有关。
定义激光雷达风速测量比和激光雷达温度测量比分别为:
风速和温度的测量灵敏度分别为SWH、STH,那么,测量得到大气中间层顶区域的风速值和温度值分别为:
其中,所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
步骤4)在上述三个激光频率工作状态下,利用光电探测器18测量得到经过碘分子吸收光谱仪后的大气瑞利散射信号的光强可以表示为:N0L(V,T)、NML(V,T)、NPL(V,T)。它们也与大气的风速V和温度T有关。
定义激光雷达风速测量比和激光雷达温度测量比分别为:
风速和温度的测量灵敏度分别为SWL、STL,那么,测量得到大气对流层、平流层的风速值和温度值分别为:
其中,所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
由于N0L(V,T)、NML(V,T)、NPL(V,T)同时与风速V和温度T有关,在进行反演时可以采用迭代算法。先给定一个大气模式下的温度,进行反演大气风速,然后代入大气风速,再进行大气温度反演,如此进行若干次循环,即可得到较为准确的大气风速V和温度T。
使用上述钠荧光多普勒激光雷达,能够同时测量中间层、平流层、对流层大气风场和温度,大大扩展了目前的探测范围。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种测量大气风场和温度的钠荧光多普勒激光雷达,其特征在于,包括:激光发射器(8)、大气散射信号接收器(9)、钠荧光散射信号测量器(10)、大气瑞利散射信号测量器(12)、信号采集器(11);
所述的激光发射器(8)用于向大气发射不同频率的589nm激光,所述的大气散射信号接收器(9)接收由大气中间层顶区域产生的钠荧光散射信号,以及由对流层、平流层产生的大气瑞利散射信号,并将其输出端分两路分别连接钠荧光散射信号测量器(10)和大气瑞利散射信号测量器(12),所述的钠荧光散射信号测量器(10)接收不同频率的钠荧光散射信号,并通过分析获得中间层顶区域的温度值和风速值;所述的大气瑞利散射信号测量器(12)接收不同频率的大气瑞利散射信号,并测量各频率的大气瑞利散射信号透过该大气瑞利散射信号测量器(12)中设有的碘分子吸收光谱仪(16)的光强后,计算获得对流层、平流层的温度值和风速值,所述的信号采集器(11)用于采集并判读钠荧光散射信号测量器(10)和大气瑞利散射信号测量器(12)输出的电信号,提取光强数据。
2.根据权利要求1所述的测量大气风场和温度的钠荧光多普勒激光雷达,其特征在于,所述的大气瑞利散射信号测量器(12)包括:准直镜(14)、背景光抑制器(15)、碘分子吸收光谱仪(16)、耦合镜(17)和光电探测器(18);所述的准直镜(14)通过光纤(13)与大气散射信号接收器(9)的输出端连接,并依次连接背景光抑制器(15)、碘分子吸收光谱仪(16)、耦合镜(17)和光电探测器(18);所述的大气瑞利散射信号依次经准直镜(14)准直操作、背景光抑制器(15)滤除天空背景光及透过碘分子吸收光谱仪(16)后,由耦合镜(17)将大气瑞利散射信号汇聚到光电探测器(18),并通过该光电探测器(18)测量透过碘分子吸收光谱仪(16)后的信号光强。
3.根据权利要求1所述的测量大气风场和温度的钠荧光多普勒激光雷达,其特征在于,所述的大气散射信号接收器(9)包括两个信号接收通道,所述的两个信号接收通道以两个不同的方向接收钠荧光散射信号和大气瑞利散射信号,且每一个信号接收通道的输出端均分两路连接钠荧光散射信号测量器(10)和大气瑞利散射信号测量器(12),使得钠荧光散射信号测量器(10)以两路接收来自两个不同方向上的钠荧光散射信号,用于测量大气中间层顶区域的温度值和水平风场,使得大气瑞利散射信号测量器(12)以两路接收来自两个不同方向上的大气瑞利散射信号,用于测量对流层、平流层的温度值和水平风场。
4.一种测量大气风场和温度的方法,其特征在于,包括:
步骤1)利用激光发射器(8)依次向大气发射三个不同频率的589nm激光;
所述激光发射器(8)向大气发射三个589nm激光的频率分别位于钠原子D2a饱和荧光光谱处f0、f0+630MHz和f0-630MHz;
步骤2)通过大气散射信号接收器(9)依次接收由大气中间层顶区域产生的三个不同频率的钠荧光散射信号后,执行步骤3),通过大气散射信号接收器(9)依次接收由对流层、平流层产生的三个不同频率的大气瑞利散射信号后,执行步骤4);
步骤3)利用钠荧光散射信号测量器(10)测量得到的三个不同频率的钠荧光散射信号所对应的光强N0H(V,T)、NMH(V,T)和NPH(V,T),定义激光雷达风速测量比和激光雷达温度测量比分别为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
则有测量得到的大气中间层顶区域的风速值和温度值分别为:
<mrow>
<mi>V</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>T</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>H</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中,SWH和STH分别表示风速测量灵敏度和温度测量灵敏度;
所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>V</mi>
</mrow>
</mfrac>
</mrow>
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
步骤4)利用光电探测器(18)测量得到透过碘分子吸收光谱仪(16)的三个不同频率的大气瑞利散射信号所对应的光强N0L(V,T)、NML(V,T)和NPL(V,T),定义激光雷达风速测量比和激光雷达温度测量比分别为:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
则有测量得到的大气对流层、平流层的风速值和温度值分别为:
<mrow>
<mi>V</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>T</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mn>0</mn>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>M</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>P</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>,</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
其中,SWL和STL分别表示风速测量灵敏度和温度测量灵敏度;
所述的风速测量灵敏度为风速变化1m/s所引起的风速测量比的变化量,其表示为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>W</mi>
<mi>L</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>V</mi>
</mrow>
</mfrac>
</mrow>
所述的温度测量灵敏度为温度变化1K引起的温度测量比的变化量,其表示为:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>T</mi>
<mi>L</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610676685.0A CN106093964B (zh) | 2016-08-16 | 2016-08-16 | 一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610676685.0A CN106093964B (zh) | 2016-08-16 | 2016-08-16 | 一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093964A CN106093964A (zh) | 2016-11-09 |
CN106093964B true CN106093964B (zh) | 2018-05-04 |
Family
ID=58070343
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610676685.0A Active CN106093964B (zh) | 2016-08-16 | 2016-08-16 | 一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093964B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108917105A (zh) * | 2018-06-13 | 2018-11-30 | 苏州若依玫信息技术有限公司 | 一种基于瑞利散射的智能家居换气系统 |
CN108535739B (zh) * | 2018-06-13 | 2023-11-28 | 中国科学技术大学 | 一种全固态连续波钠测温测风激光雷达 |
CN108828624B (zh) * | 2018-06-29 | 2020-11-24 | 中科院合肥技术创新工程院 | 一种高灵敏度和高线性度兼顾的大气瑞利回波光信号接收方法和装置 |
CN110703278A (zh) * | 2019-11-05 | 2020-01-17 | 中国科学院武汉物理与数学研究所 | 一种钠层层析观测激光雷达及观测方法 |
CN112596075B (zh) * | 2020-11-26 | 2022-03-25 | 兰州大学 | 一种多激发波长光谱仪式荧光激光雷达系统 |
CN113358160B (zh) * | 2021-07-08 | 2022-05-17 | 天津大学 | 一种大气数据测量方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201021940Y (zh) * | 2007-02-13 | 2008-02-13 | 中国科学院武汉物理与数学研究所 | 一种综合型多功能中高层大气探测激光雷达 |
CN101819275A (zh) * | 2010-04-20 | 2010-09-01 | 中国海洋大学 | 测量多气象参数的多普勒激光雷达装置 |
CN102749627A (zh) * | 2012-06-02 | 2012-10-24 | 中国科学院武汉物理与数学研究所 | 一种全高程大气风温密探测激光雷达 |
CN104345319A (zh) * | 2014-10-12 | 2015-02-11 | 中国海洋大学 | 一种非相干多普勒激光雷达风速的误差校正方法 |
-
2016
- 2016-08-16 CN CN201610676685.0A patent/CN106093964B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201021940Y (zh) * | 2007-02-13 | 2008-02-13 | 中国科学院武汉物理与数学研究所 | 一种综合型多功能中高层大气探测激光雷达 |
CN101819275A (zh) * | 2010-04-20 | 2010-09-01 | 中国海洋大学 | 测量多气象参数的多普勒激光雷达装置 |
CN102749627A (zh) * | 2012-06-02 | 2012-10-24 | 中国科学院武汉物理与数学研究所 | 一种全高程大气风温密探测激光雷达 |
CN104345319A (zh) * | 2014-10-12 | 2015-02-11 | 中国海洋大学 | 一种非相干多普勒激光雷达风速的误差校正方法 |
Non-Patent Citations (6)
Title |
---|
532nm测风激光雷达长时间稳频系统;郭文杰 等;《红外与激光工程》;20160531;第45卷(第S1期);1-5 * |
Development of a mobile doppler lidar system for wind and temperature measurements at 30-70 km;Zhaoai Yan 等;《Journal of Quantitative Spectroscopy & Radiative Transfer》;20160504;52-59 * |
Simultaneous three-frequency Na lidar measurements of radial wind and temperature in the mesopause region;C.Y.She 等;《GEOPHYSICAL RESEARCH LETTERS》;19940805;第21卷(第17期);1771-1774 * |
利用瑞利激光雷达观测北京地区上平流层地形重力波活动;郭文杰 等;《地球物理学报》;20151031;第58卷(第10期);3481-3486 * |
星载钠荧光多普勒激光雷达性能分析;闫召爱 等;《雷达学报》;20150228;第4卷(第1期);99-106 * |
钠多普勒激光雷达回波光子数仿真及大气参数反演;徐丽 等;《地球物理学报》;20100731;第53卷(第7期);1520-1528 * |
Also Published As
Publication number | Publication date |
---|---|
CN106093964A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093964B (zh) | 一种测量大气风场和温度的钠荧光多普勒激光雷达及方法 | |
CN104237161B (zh) | 一种煤炭自燃指标气体多组分实时在线远程监测装置及方法 | |
CN102636459B (zh) | 一种前向散射与透射结合的能见度测量仪及其测量方法 | |
CN105334519B (zh) | 基于三通道f-p标准具的多大气参数同时高精度探测激光雷达系统 | |
Hansen et al. | First detection of a noctilucent cloud by lidar | |
CN101477196B (zh) | 振动拉曼激光雷达散射光处理系统及其处理方法 | |
CN105182365B (zh) | 一种基于双波长单接收通道的转动拉曼测温激光雷达 | |
CN108303706A (zh) | 一种气溶胶光学参数探测方法及高光谱激光雷达探测系统 | |
Hoareau et al. | A Raman lidar at La Reunion (20.8 S, 55.5 E) for monitoring water vapour and cirrus distributions in the subtropical upper troposphere: preliminary analyses and description of a future system | |
CN106526611A (zh) | 一种利用激光雷达系统测量海洋固有光学参数的方法 | |
CN106443709A (zh) | 基于时分复用的单探测器全光纤偏振激光雷达 | |
Zhang et al. | Sensitivity of estimated total canopy SIF emission to remotely sensed LAI and BRDF products | |
CN201340455Y (zh) | 振动拉曼激光雷达散射光处理系统 | |
CN107389560A (zh) | 多波段全光纤高光谱分辨率整层大气透过率同时测量装置及测量方法 | |
CN103234942B (zh) | 利用天空散射光测量大气水平能见度的方法及装置 | |
CN107064957A (zh) | 一种用于液态水云测量的多视场激光雷达探测系统及方法 | |
Yan et al. | Modeling the radiation regime of a discontinuous canopy based on the stochastic radiative transport theory: Modification, evaluation and validation | |
Poggio et al. | Scintillometer wind measurements over complex terrain | |
Liu et al. | Micrometeorological methods to determine evapotranspiration | |
CN105352861A (zh) | 浮游植物粒级结构分析方法 | |
CN103149180B (zh) | 一种土壤光谱反射率和电导率检测方法 | |
CN102692191B (zh) | 基于光纤传感测量公路路面表层水膜厚度的方法 | |
CN106769895A (zh) | 一种标定测量整层大气光谱透过率的方法 | |
CN106483531A (zh) | 大气拉曼‑瑞利散射测温激光雷达及反演方法 | |
He et al. | Modis aerosol optical thickness product algorithm verification and analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |