CN112882062B - 天基co2通量激光探测装置 - Google Patents
天基co2通量激光探测装置 Download PDFInfo
- Publication number
- CN112882062B CN112882062B CN202110053184.8A CN202110053184A CN112882062B CN 112882062 B CN112882062 B CN 112882062B CN 202110053184 A CN202110053184 A CN 202110053184A CN 112882062 B CN112882062 B CN 112882062B
- Authority
- CN
- China
- Prior art keywords
- laser
- flux
- frequency
- resolving
- wind field
- 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
- 230000004907 flux Effects 0.000 title claims abstract description 63
- 238000001514 detection method Methods 0.000 title claims abstract description 42
- 230000003287 optical effect Effects 0.000 claims abstract description 63
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000005259 measurement Methods 0.000 claims description 36
- 238000010521 absorption reaction Methods 0.000 claims description 35
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 22
- 238000000034 method Methods 0.000 description 21
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 18
- 229910002092 carbon dioxide Inorganic materials 0.000 description 9
- 230000001427 coherent effect Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 6
- 229910052799 carbon Inorganic materials 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004177 carbon cycle Methods 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
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
-
- 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/006—Theoretical aspects
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种天基CO2通量激光探测装置,包括激光模块、收发模块以及解算模块;激光模块包括主激光器、第一单频激光器、第二单频激光器及光开关;光开关用于控制交替发出第一单频激光器、第二单频激光器的输出激光;收发模块包括光学收发装置及伺服系统,光学收发装置将输出激光发射至测量区域,以及接收测量区域的回波信号;解算模块用于根据回波信号进行信号处理解算得到测量区域的风场廓线信息及CO2浓度廓线信息,并计算得到CO2通量数据。本发明的天基CO2通量激光探测装置设置于天基平台可以高精度、高可靠性地测量地表的CO2收支量。
Description
技术领域
本发明涉及光电子技术领域,具体而言,涉及一种天基CO2通量激光探测装置。
背景技术
大气中CO2等温室气体浓度不断升高引起的全球变暖、极地冰盖消融、海平面上升、生态系统的物种组成改变等全球变化已成为人类最为关注的环境问题,研究地球系统的碳循环机理、全球陆地生态系统碳收支及其环境变换响应已成为当前一系列大型国际研究计划的共同核心问题。在全球碳循环和碳收支研究中发现了“碳失踪”问题,并且推测这些“失踪的碳”可能被蓄积于陆地生态系统之中。因此。如何准确地预测地球陆地生态系统与大气间净CO2交换量、分析其时空分布格局,便成为地球物理、生物学家们一致致力于解决的重大科学问题。
目前观测全球生态系统的碳收支的方法是涡度相关法。它是通过快速测定大气中CO2浓度与风速的协方差来计算湍流通量的一种方法,它是基于大气湍流理论和数据统计分析相结合的一种计算。其中涡度相关法经过长期的理论发展和技术改进,已经实现了对森林、草地等局部生态系统CO2通量的测量。随着涡度相关通量观测在全球范围内的广泛开展,一个覆盖全球陆地生态系统的国际碳通量观测研究网络应运而生。但通量观测站点均处于远离市区的野外,经常会由于恶劣天气、系统故障或者其他外界干扰而造成通量观测仪器的失灵,并由此使得通量观测数据在数量和质量上达不到研究的要求。同时全面反演全球的碳分布将需要大量的CO2探测站,因此实现难度大,成本高。除此之外,多光谱衍射光栅也是测量CO2浓度分布的常用方法,它是在可见光和近红外谱段,利用分子吸收谱线探测二氧化碳浓度,但该测量方法无法做到全天候、全天时工作,且仅能测量CO2浓度,无法得到风场信息,因此也就无法反演CO2通量。
实现全球覆盖的CO2浓度测量的高效手段是利用卫星遥感。例如,基于星载差分光谱探测原理的被动遥感方案,是当前较为主流的一种获取全球CO2浓度分布的方式,可以在可见光和近红外谱段,利用分子吸收谱线探测二氧化碳浓度,但该测量方法无法做到全天候、全天时工作,且仅能实现大尺度范围的CO2浓度测量,进而获得CO2通量数据产品,但该产品已经是载荷4-5级数据产品,不仅分辨率极低,而且数据的可靠性大幅降低。目前已经规划的天基差分吸收激光雷达,可以提供更高分辨率的浓度数据,但其任务仍是以柱浓度为探测目标,不具备通量探测能力。
因此,现有基于卫星遥感的CO2通量探测方式存在分辨率及可靠性低的问题。
发明内容
本发明解决的是现有基于卫星遥感的CO2通量探测方式存在分辨率及可靠性低的问题。
为解决上述问题,本发明提供一种天基CO2通量激光探测装置,所述装置包括激光模块、收发模块以及解算模块;所述激光模块包括主激光器、第一单频激光器、第二单频激光器及光开关;所述主激光器的输出波长为 CO2吸收峰波长,所述第一单频激光器的输出波长由所述CO2吸收峰波长进行天基平台对应的多普勒频移补偿得到,所述第二单频激光器的输出波长远离所述CO2吸收峰波长;所述光开关用于控制交替发出所述第一单频激光器、所述第二单频激光器的输出激光;所述收发模块包括光学收发装置及伺服系统,所述光学收发装置用于将所述输出激光发射至测量区域,以及接收所述测量区域的回波信号;所述伺服系统用于控制所述光学收发装置的出射角度;所述解算模块用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息及CO2浓度廓线信息,并根据所述风场廓线信息及所述CO2浓度廓线信息计算得到CO2通量数据。
可选地,所述激光模块还包括:脉冲激光器及激光放大器;所述光开关、所述脉冲激光器、所述激光放大器、所述光学收发装置沿光路依次设置;所述脉冲激光器用于将所述光开关输出的连续激光转换为脉冲激光输出;所述激光放大器用于将所述脉冲激光放大至探测所需功率。
可选地,所述CO2吸收峰波长λpeak,所述天基平台对应的多普勒频移对应波长为Δλ1,所述第二单频激光器的输出波长与所述CO2吸收峰波长的差值为Δλ2;根据跟随锁定技术将所述第一单频激光器的输出波长λon锁定至λpeak+Δλ1,以及将所述第二单频激光器的输出波长λoff锁定至λpeak+Δλ2。
可选地,所述伺服系统用于控制所述光学收发装置将所述输出激光前向发射至所述测量区域,以及将所述输出激光后向发射至所述测量区域;所述前向发射的角度与所述后向发射的角度相同;所述光学收发装置用于接收前向发射至所述测量区域的回波信号,以及接收后向发射至所述测量区域的回波信号。
可选地,所述解算模块包括风场廓线解算单元、CO2浓度廓线解算单元及CO2通量解算单元;所述风场廓线解算单元用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息;所述CO2浓度廓线解算单元用于根据所述回波信号解算得到所述测量区域的CO2浓度廓线信息;或者,所述解算模块包括信号处理单元及CO2通量解算单元;所述信号处理单元用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息及CO2浓度廓线信息;所述CO2通量解算单元用于根据所述风场廓线信息及所述CO2浓度廓线信息计算得到CO2通量数据。
可选地,所述解算模块还包括分光装置;所述分光装置用于将所述测量区域的回波信号分成两路并进行光电转换,将两路电信号分别输入所述风场廓线解算单元及所述CO2浓度廓线解算单元。
可选地,所述第一单频激光器、所述第二单频激光器通过分光光路与所述风场廓线解算单元连接,用于将激光输出至所述风场廓线解算单元作为本振光信号;所述风场廓线解算单元,用于根据所述本振光信号及所述测量区域对应回波信号的外差相干进行风场廓线的反演,得到风场廓线信息;所述测量区域对应回波信号包括所述第一单频激光器的回波信号及所述第二单频激光器的回波信号。
可选地,所述CO2浓度廓线解算单元用于测量所述回波信号的功率,以及根据以下浓度反演公式解算CO2浓度廓线:
其中,Δσ=σon-σoff,所述第一单频激光器的输出激光的波长为λon、吸收系数为σon,所述第二单频激光器的输出激光的波长为λoff、吸收系数为σoff,距离R处的回波功率值为P(R,λON)和P(R,λOFF),距离R+ΔR处的回波功率P(R+ΔR,λON)和P(R+ΔR,λOFF)。
可选地,所述CO2浓度廓线解算单元包括低通滤波器,所述低通滤波器用于过滤所述回波信号中的高频噪声。
本发明实施例提供的天基CO2通量激光探测装置,设置于天基平台,激光模块输出波长处于CO2吸收峰和偏离吸收峰的两种单频激光,通过光学收发装置将激光发射至测量区域以及接收回波信号得到测量结果,解算模块接收到测量结果后可以进行风场廓线及CO2浓度廓线双重测量,测量结果可以直接用于反演测量区域的CO2通量,可以高精度、高可靠性地测量地表的CO2收支量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例的一种天基CO2通量激光探测装置的结构示意图;
图2为本发明实施例的天基水平风场的前向和后向组合观测示意图;
图3为本发明实施例的另一种天基CO2通量激光探测装置的结构示意图。
附图标记说明:
101-主激光器;102-第一单频激光器;103-第二单频激光器;104-光开关;105-脉冲激光器;106-激光放大器;201-光学收发装置;202-伺服系统; 301-风场廓线解算单元;302-CO2浓度廓线解算单元;303-CO2通量解算单元;304-分光装置;305-信号处理单元;306-相干探测单元。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施例做详细的说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例提供了一种天基CO2通量激光探测装置,包括激光模块、收发模块以及解算模块。图1为本发明实施例的一种天基CO2通量激光探测装置的结构示意图。
(1)激光模块包括主激光器101、第一单频激光器102、第二单频激光器103及光开关104。可选地,激光模块可以输出2μm单频脉冲激光。
主激光器101的输出波长为CO2吸收峰波长λpeak,第一单频激光器102 的输出波长由CO2吸收峰波长进行天基平台对应的多普勒频移补偿得到,第二单频激光器103的输出波长远离CO2吸收峰波长;CO2吸收峰波长天基平台对应的多普勒频移对应波长为Δλ1,第二单频激光器103的输出波长与CO2吸收峰波长的差值为Δλ2;根据跟随锁定技术将第一单频激光器102 的输出波长λon锁定至λpeak+Δλ1,以及将第二单频激光器103的输出波长λoff锁定至λpeak+Δλ2。
如图1所示,波长λpeak的主激光器101(采用CO2吸收池作为参考源将频率锁定至CO2吸收峰处),波长λpeak+Δλ1 (Δf1=|cΔλ1/(λpeak-Δλ1)2|)的第一单频激光器102,以主激光器101 作为基准,采用跟随锁定技术将第一单频激光器102的频率锁定至λpeak+Δλ1,同时可以调节外部参考源的频率实现Δλ1的动态调节。Δλ1(Δf1) 用于补偿天基平台运动引起的多普勒频移。
波长λpeak+Δλ2(Δf2=|cΔλ2/(λpeak+Δλ2)2|)的第二单频激光器103,同样以主激光器101作为基准,采用跟随锁定技术将第二单频激光器103 的频率锁定至λpeak+Δλ2,同时可以调节外部参考源的频率实现Δλ2的动态调节。Δλ2(Δf2)的作用是调节激光器的频率使之远离CO2吸收峰。同时假设v为天基平台在激光传播方向的分速度,那么可以通过以下公式计算由于天基平台运动产生多普勒频移:
Δf′1=2v/(λpeak+Δλ1),Δf′2=2v/(λpeak+Δλ2)。
然后,根据跟随锁定技术设定频率偏移Δλ1(或Δλ2),计算确定相对于静止第一单频激光器102的出射激光波长λon和第二单频激光器103出射的激光波长λoff。
光开关104用于控制交替发出第一单频激光器102、第二单频激光器 103的输出激光。第一单频激光器102、第二单频激光器103输出的单频连续激光注入光开关104中,光开关104可以实现不同波长选通,从而使得波长为λpeak+Δλ1和λpeak+Δλ2的单频连续激光交替依次发出。
如图1所示,在激光模块中还包括:脉冲激光器105及激光放大器106。光开关104、脉冲激光器105、激光放大器106、光学收发装置201沿光路依次设置。脉冲激光器105用于将光开关输出的连续激光转换为脉冲激光输出;激光放大器106用于将脉冲激光放大至探测所需功率。
(2)收发模块包括光学收发装置201及伺服系统202。
光学收发装置201将输出激光发射至测量区域,以及接收该测量区域的回波信号;伺服系统202用于控制上述光学收发装置201的出射角度。
伺服系统202用于控制光学收发装置201将输出激光前向发射至测量区域,以及将输出激光后向发射至测量区域。上述前向发射的角度与后向发射的角度相同。光学收发装置201可以接收前向发射至测量区域的回波信号,以及接收后向发射至测量区域的回波信号。
伺服系统202可以控制光学收发装置201偏转。为了实现水平风场廓线的反演,伺服系统202通过控制发光学装置201实现前向和后向两个方向的探测,利用两次回波信号反演水平风速。参见图2所示的天基水平风场的前向和后向组合观测示意图,天基平台沿箭头方向运行,A点为前向激光出射位置,经大气传输至F点后散射,并在B接收;D点为后向激光出射位置,经大气传输至F点后散射,并在E点接收。
伺服系统202还与光开关104连接,用于控制光开关104发出的激光,即控制光开关104允许第一单频激光器102的输出激光通过,或者允许第二单频激光器103的输出激光通过。
(3)本实施例提供了解算模块的两种可行方案。作为第一种可行方案,解算模块包括风场廓线解算单元301、CO2浓度廓线解算单元302及CO2通量解算单元303,作为第二种可行方案,解算模块包括信号处理单元305 及CO2通量解算单元303。
在上述第一种可行方案中,风场廓线解算单元301采用相干探测解算风场的多普勒频移,反演测量区域的风场廓线信息。CO2浓度廓线解算单元 302采用直接探测解算回波信号功率,反演测量区域的CO2浓度廓线信息。 CO2通量解算单元303用于根据风场廓线信息及CO2浓度廓线信息计算得到CO2通量数据。
进一步,解算模块还可以包括分光装置304;该分光装置304用于将测量区域的回波信号分成两路,将两路回波光信号分别输入风场廓线解算单元301及CO2浓度廓线解算单元302。需要说明的是,第一单频激光器102、第二单频激光器103通过分光光路与风场廓线解算单元301连接,用于将激光输出至风场廓线解算单元301作为本振光信号。
具体地,上述风场廓线解算单元301,用于根据本振光信号及测量区域对应回波信号的外差相干进行风场廓线的反演,得到风场廓线信息;测量区域对应回波信号包括第一单频激光器的回波信号及第二单频激光器的回波信号。光学收发装置201接收回波信号连接至分光装置304,将回波光信号分为两路并进行光电转换;两路信号分别通过风场廓线解算单元301及 CO2浓度廓线解算单元302进行信号处理,得到测量区域的风场廓线和CO2浓度廓线信息。根据前向和后向的径向风场值以及前向、后向在地球坐标系下的位置,解算得到水平风场廓线,并传输至CO2通量解算单元303,根据水平风场廓线和CO2浓度廓线得到CO2通量。
在上述第二种可行方案中,信号处理单元305采用相干探测方式解算测量区域的多普勒频移和强度信息,以反演风场廓线信息及CO2浓度廓线信息;CO2通量解算单元303用于根据风场廓线信息及CO2浓度廓线信息计算得到CO2通量数据。参见图3所示的另一种天基CO2通量激光探测装置的结构示意图,与图1中示出的天基CO2通量激光探测装置的结构不同,上述解算模块包括信号处理单元305、相干探测单元306及CO2通量解算单元303。天基CO2通量激光探测装置的除解算模块外的其他组成部分与上述第一种可行方案中相同,在此不再赘述。
光学收发装置201接收回波信号连接至相干探测单元306,信号处理单元305用于根据回波信号进行信号处理解算。需要说明的是,信号处理单元305可以实现上述风场廓线解算单元301及CO2浓度廓线解算单元302 的全部功能。
本发明实施例提供的天基CO2通量激光探测装置,设置于天基平台,激光模块输出波长处于CO2吸收峰和偏离吸收峰的两种单频激光,通过光学收发装置将激光发射至测量区域以及接收回波信号得到测量结果,解算模块接收到测量结果后可以进行风场廓线及CO2浓度廓线双重测量,测量结果可以直接用于反演测量区域的CO2通量,可以高精度、高可靠性地测量地表的CO2收支量。
可选地,以解算模块的第一种可行方案为例,风场廓线解算单元301、 CO2浓度廓线解算单元302的具体功能分述如下:
风场廓线解算单元301采用激光相干探测方式,利用回波光信号与本振光的外差相干实现风场廓线的反演。
具体流程为:回波光信号与对应的本振光信号在探测器端面相干产生射频信号,该信号含有激光的多普勒频移,此时λon和λoff均参与风场廓线解算。CO2浓度廓线解算单元包括低通滤波器,低通滤波器用于过滤回波信号中的高频噪声。
由于激光出射后远距离目标的回波功率极低,该区域的信号被认为完全是噪声,将其作为噪声基底来对消近距离回波的噪声;然后依次通过数据累加提高回波载噪比、带通滤波器滤除带宽以外的噪声、射频放大器提高回波信号强度。
风场廓线测量是通过短时傅里叶变换实现风场廓线的解算,具体为设定距离分辨率,调节短时傅里叶变换的重叠因子得到随时间变化的多普勒频移,根据时间-多普勒关系换算得到风场廓线。例如,回波射频信号依次通过带通滤波器滤除噪声、射频放大器提高回波信号强度,最后设定风场廓线的距离分辨率,调节短时傅里叶变换的重叠因子得到随时间变化的多普勒频移,根据时间-多普勒关系换算得到风场廓线信息v(R)。
CO2浓度廓线解算单元302采用直接探测方式,通过连续依次测量回波信号λon和λoff的功率来反演CO2浓度。
具体流程为:回波光信号进入单光子探测器,得到回波信号功率,经过数据解析分辨ON/OFF回波信号,然后依次通过数据累加降低功率抖动。由于激光出射后远距离目标的回波功率极低,该区域的信号被认为完全是噪声,将其作为噪声基底来对消近距离回波的噪声。
回波功率是通过检测外差信号的强度来实现的,由于大气湍流和散斑噪声的影响,信号功率难以精确测量。
回波信号的载噪比CNR=<Pi>/<Pn>(其中,<Pi>是回波信号的平均值,<Pn>是回波信号的噪声)与信噪比SNR=<Pi>/σ(<Pi>)满足以下公式:
其中,Mp是累计次数,δtR为测量时间,Tc为脉冲宽度,通常δtR≤Tc。为了提高信号的强度解调精度,回波功率抖动需要尽可能减低,即SNR越高,功率测量精度越高。从上述公式可以发现,当载噪比较大时,回波信噪比与成正比,与成正比。因此提高累计次数能显著提高回波信噪比,降低强度抖动。当CO2浓度数据的刷新率固定时,提高脉冲的重复频率是提高测量精度的主要方法。
然后,CO2浓度廓线解算单元测量回波信号的功率,以及根据以下浓度反演公式解算CO2浓度廓线:
其中,Δσ=σon-σoff,第一单频激光器的输出激光的波长为λon、吸收系数为σon,第二单频激光器的输出激光的波长为λoff、吸收系数为σoff,距离R处的回波功率值为P(R,λON)和P(R,λOFF),距离R+ΔR处的回波功率 P(R+ΔR,λON)和P(R+ΔR,λOFF)。由于仅测量功率值,因此通过低通滤波器滤去高频噪声。
联合分析CO2浓度廓线和风场廓线,获得CO2通量数据。例如,在近地面处,根据边界层理论,根据CO2浓度廓线和风场廓线迭代计算可得到大气边界层CO2通量数据。
需要说明的是,信号处理单元305可以实现风场廓线解算单元301及 CO2浓度廓线解算单元302的全部功能。
本发明实施例提供的天基CO2通量激光探测装置,利用单频脉冲激光的外差相干探测原理实现径向风场测量,两个波长分别处于CO2吸收峰和偏离吸收峰的直接差分吸收探测原理实现CO2浓度测量,具备风场廓线和 CO2浓度廓线分布的双重测量功能。为了反演水平风速廓线,提出了前向和后向组合观测的方案,测量结果能直接用于反演测量区域的CO2通量,进而高精度地测量地表的CO2收支量,构建全球CO2循环测量体系。
当然,本领域技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程度来指令控制装置来完成,所述的程序可存储于一计算机可读取的存储介质中,所述程序在执行时可包括如上述各方法实施例的流程,其中所述的存储介质可为存储器、磁盘、光盘等。
在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (9)
1.一种天基CO2通量激光探测装置,其特征在于,所述装置包括激光模块、收发模块以及解算模块;
所述激光模块包括主激光器、第一单频激光器、第二单频激光器及光开关;所述主激光器的输出波长为CO2吸收峰波长,所述第一单频激光器的输出波长由所述CO2吸收峰波长进行天基平台对应的多普勒频移补偿得到,所述第二单频激光器的输出波长远离所述CO2吸收峰波长;所述光开关用于控制交替发出所述第一单频激光器、所述第二单频激光器的输出激光;
所述收发模块包括光学收发装置及伺服系统,所述光学收发装置用于将所述输出激光发射至测量区域,以及接收所述测量区域的回波信号;所述伺服系统用于控制所述光学收发装置的出射角度;
所述解算模块用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息及CO2浓度廓线信息,并根据所述风场廓线信息及所述CO2浓度廓线信息计算得到CO2通量数据。
2.根据权利要求1所述的天基CO2通量激光探测装置,其特征在于,所述激光模块还包括:脉冲激光器及激光放大器;
所述光开关、所述脉冲激光器、所述激光放大器、所述光学收发装置沿光路依次设置;
所述脉冲激光器用于将所述光开关输出的连续激光转换为脉冲激光输出;
所述激光放大器用于将所述脉冲激光放大至探测所需功率。
3.根据权利要求1所述的天基CO2通量激光探测装置,其特征在于,所述CO2吸收峰波长λpeak,所述天基平台对应的多普勒频移对应波长为Δλ1,所述第二单频激光器的输出波长与所述CO2吸收峰波长的差值为Δλ2;
根据跟随锁定技术将所述第一单频激光器的输出波长λON 锁定至λpeak+Δλ1,以及将所述第二单频激光器的输出波长λOFF 锁定至λpeak+Δλ2。
4.根据权利要求1-3任一项所述的天基CO2通量激光探测装置,其特征在于,所述伺服系统用于控制所述光学收发装置将所述输出激光前向发射至所述测量区域,以及将所述输出激光后向发射至所述测量区域;所述前向发射的角度与所述后向发射的角度相同;
所述光学收发装置用于接收前向发射至所述测量区域的回波信号,以及接收后向发射至所述测量区域的回波信号。
5.根据权利要求1所述的天基CO2通量激光探测装置,其特征在于,所述解算模块包括风场廓线解算单元、CO2浓度廓线解算单元及CO2通量解算单元;所述风场廓线解算单元用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息;所述CO2浓度廓线解算单元用于根据所述回波信号解算得到所述测量区域的CO2浓度廓线信息;或者,
所述解算模块包括信号处理单元及CO2通量解算单元;所述信号处理单元用于根据所述回波信号进行信号处理解算得到所述测量区域的风场廓线信息及CO2浓度廓线信息;
所述CO2通量解算单元用于根据所述风场廓线信息及所述CO2浓度廓线信息计算得到CO2通量数据。
6.根据权利要求5所述的天基CO2通量激光探测装置,其特征在于,所述解算模块还包括分光装置;
所述分光装置用于将所述测量区域的回波信号分成两路并进行光电转换,将两路电信号分别输入所述风场廓线解算单元及所述CO2浓度廓线解算单元。
7.根据权利要求6所述的天基CO2通量激光探测装置,其特征在于,所述第一单频激光器、所述第二单频激光器通过分光光路与所述风场廓线解算单元连接,用于将激光输出至所述风场廓线解算单元作为本振光信号;
所述风场廓线解算单元,用于根据所述本振光信号及所述测量区域对应回波信号的外差相干进行风场廓线的反演,得到风场廓线信息;所述测量区域对应回波信号包括所述第一单频激光器的回波信号及所述第二单频激光器的回波信号。
9.根据权利要求5-7任一项所述的天基CO2通量激光探测装置,其特征在于,
所述CO2浓度廓线解算单元包括低通滤波器,所述低通滤波器用于过滤所述回波信号中的高频噪声。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110053184.8A CN112882062B (zh) | 2021-01-15 | 2021-01-15 | 天基co2通量激光探测装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110053184.8A CN112882062B (zh) | 2021-01-15 | 2021-01-15 | 天基co2通量激光探测装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112882062A CN112882062A (zh) | 2021-06-01 |
CN112882062B true CN112882062B (zh) | 2022-12-02 |
Family
ID=76047967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110053184.8A Active CN112882062B (zh) | 2021-01-15 | 2021-01-15 | 天基co2通量激光探测装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112882062B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113433570B (zh) * | 2021-07-15 | 2022-11-01 | 北京遥测技术研究所 | 一种大气二氧化碳浓度探测差分吸收激光雷达系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101493406A (zh) * | 2008-11-04 | 2009-07-29 | 中国人民解放军电子工程学院 | 基于差分吸收的大气二氧化碳远距离激光相干探测装置 |
CN104035102A (zh) * | 2014-06-12 | 2014-09-10 | 中国科学院上海技术物理研究所 | 一种探测大气二氧化碳浓度的激光雷达系统及方法 |
CN107102309A (zh) * | 2017-04-28 | 2017-08-29 | 北京怡孚和融科技有限公司 | 风场空间分布转化为风场时间方向的方法及气溶胶时空分布与风场时空分布的叠加方法 |
CN109459761A (zh) * | 2018-12-20 | 2019-03-12 | 南京牧镭激光科技有限公司 | 一种激光雷达 |
CN109655843A (zh) * | 2019-01-16 | 2019-04-19 | 武汉大学 | 探测二氧化碳浓度廓线的脉冲红外差分吸收激光雷达系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012083267A (ja) * | 2010-10-13 | 2012-04-26 | Japan Aerospace Exploration Agency | マルチライダーシステム |
-
2021
- 2021-01-15 CN CN202110053184.8A patent/CN112882062B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101493406A (zh) * | 2008-11-04 | 2009-07-29 | 中国人民解放军电子工程学院 | 基于差分吸收的大气二氧化碳远距离激光相干探测装置 |
CN104035102A (zh) * | 2014-06-12 | 2014-09-10 | 中国科学院上海技术物理研究所 | 一种探测大气二氧化碳浓度的激光雷达系统及方法 |
CN107102309A (zh) * | 2017-04-28 | 2017-08-29 | 北京怡孚和融科技有限公司 | 风场空间分布转化为风场时间方向的方法及气溶胶时空分布与风场时空分布的叠加方法 |
CN109459761A (zh) * | 2018-12-20 | 2019-03-12 | 南京牧镭激光科技有限公司 | 一种激光雷达 |
CN109655843A (zh) * | 2019-01-16 | 2019-04-19 | 武汉大学 | 探测二氧化碳浓度廓线的脉冲红外差分吸收激光雷达系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112882062A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5847817A (en) | Method for extending range and sensitivity of a fiber optic micro-doppler ladar system and apparatus therefor | |
CN110441792B (zh) | 同时测风测温的瑞利散射激光雷达系统及相关校准方法 | |
CN109143263B (zh) | 一种混合型测风激光雷达 | |
US5889490A (en) | Method and apparatus for improved ranging | |
US8077294B1 (en) | Optical autocovariance lidar | |
US11125879B2 (en) | Method for processing a signal arising from coherent lidar and associated lidar system | |
CN110749872A (zh) | 一种相干差分吸收激光雷达与一种检测气体浓度的方法 | |
Onori et al. | Coherent interferometric dual-frequency laser radar for precise range/Doppler measurement | |
US10295672B2 (en) | System and method to measure an atmospheric thermodynamic profile with a compact, all-fiber and eye-safe Lidar | |
CN107688187B (zh) | 基于空间波长编码的目标探测激光雷达 | |
US20210382164A1 (en) | Multi-tone continuous wave detection and ranging | |
CN114019525B (zh) | 一种基于光梳的高精度激光光谱测距方法 | |
CN115290599B (zh) | 一种测量温室气体浓度的激光雷达系统 | |
Torun et al. | Realization of multitone continuous wave LiDAR | |
CN112882062B (zh) | 天基co2通量激光探测装置 | |
CN111721968B (zh) | 一种基于双光梳系统测定气体流速的方法 | |
Sun et al. | Comparison of IPDA lidar receiver sensitivity for coherent detection and for direct detection using sine-wave and pulsed modulation | |
Liu et al. | Differential discrimination technique for incoherent Doppler lidar to measure atmospheric wind and backscatter ratio | |
CN111208084A (zh) | 一种基于相干探测方法的光纤气体浓度遥感探测装置和方法 | |
CN115267828A (zh) | 基于差分相关脉冲的激光雷达测风方法及系统 | |
CN115290175B (zh) | 一种海水声速测量装置、方法以及海洋测距系统 | |
US4867564A (en) | Apparatus for and method of remotely sensing sub-surface water temperatures | |
Singh et al. | Evaluation of 2-μm pulsed integrated path differential absorption lidar for carbon dioxide measurement—technology developments, measurements, and path to space | |
CN114167440B (zh) | 一种基于相位噪声补偿的相干调频连续波测距方法及系统 | |
Totems et al. | Advanced signal processing methods for pulsed laser vibrometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |