CN102879772B - 一种差分吸收大气co2激光雷达大气压力增宽修正方法 - Google Patents
一种差分吸收大气co2激光雷达大气压力增宽修正方法 Download PDFInfo
- Publication number
- CN102879772B CN102879772B CN201210351492.XA CN201210351492A CN102879772B CN 102879772 B CN102879772 B CN 102879772B CN 201210351492 A CN201210351492 A CN 201210351492A CN 102879772 B CN102879772 B CN 102879772B
- Authority
- CN
- China
- Prior art keywords
- broadening
- atmospheric
- laser radar
- air
- gas
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种差分吸收大气CO2激光雷达大气压力增宽修正方法,其特点是:首先采用对波长敏感的纠正公式进行洛伦兹半峰全宽FWHM的计算,兼顾增宽系数、温度相关指数的独立性,采用直接带值和迭代求解两种思路解决高空CO2混合比动值问题,并将差分吸收激光雷达方程浓度解——CO2分子浓度与体积比χ进行转化。采用本发明后,能够最大程度上摒弃经典算法在取值上的不确定性,避免新误差的带入,同时根据实验场景需求,具体采用带值求解或迭代的方法具有灵活性,在差分吸收激光雷达测量大气CO2领域具有很好的应用前景和极大的发展潜力。
Description
技术领域
本发明涉及大气光学领域,尤其涉及一种差分吸收大气CO2激光雷达大气压力增宽修正方法。
背景技术
CO2作为温室气体主要成分对全球气候有重要的影响作用,其浓度的变化,会对全球气候产生重要影响。通过对CO2浓度的实时监测,可以自主、适时地获得全面、可靠的区域CO2源汇动态信息,应对日益增长的碳排放压力问题。差分吸收激光雷达在探测CO2方面具有巨大优势,是观测其垂直廓线分布的理想探测设备,能够提供连续的时间和空间观测,获得较高的时空分辨率。对于差分吸收激光雷达,可用信号表现为被探测大气物质对激光束能量的吸收。同时为了排除其他因素的影响,在差分吸收激光雷达中通常采用两束差别很小的波长作为发射激光束。其中一束激光的波长选在CO2吸收峰的中心,使其受到最大吸收,称为λon,另一束激光的波长在吸收谷处,使其受到的吸收尽可能小,称为λoff。通过测量随高度增加的大气CO2吸收的变化来反演CO2浓度分布。
然而,在利用差分吸收激光雷达探测CO2浓度时,大气压力增宽对其有不可忽视的重要影响,大气压力增宽的修正对于CO2激光雷达的浓度反演是极其重要的。
高度R处厚度为ΔR=R2-R1的大气层中(R即位于大气层较高的高度R2和大气层较低的高度R1之间的中间高度),差分吸收激光雷达探测的CO2浓度方程可以表现为如下形式:
式中Ng代表CO2密度的高度分布,σg(λ)代表被探测组分的吸收截面积,是波长λ的敏感函数,其中σg(λon)是波长λon的敏感函数,σg(λoff)是波长λoff的敏感函数;Ps(,)代表回波信号强度,Ps(λoff,R2)表示采用波长λoff探测高度R2处得到的回波信号强度,Ps(λon,R1)表示采用波长λon探测高度R1处得到的回波信号强度,Ps(λon,R2)表示采用波长λon探测高度R2处 得到的回波信号强度,Ps(λoff,R1)表示采用波长λoff探测高度R1处得到的回波信号强度。由上式可见,被探测组分的密度与两通道回波信号之比和两波长处的吸收截面积之差有关。σg(λ)的值与洛伦兹半峰全宽(FWHM)相关,而大气压强以及温度的变化会直接导致FWHM的变化,称为压力增宽,其谱线形状由洛伦兹廓线给出。随着高度的增加,压强同温度存在较大的波动,导致洛伦兹廓线有所变化,必须对其进行修正,即大气压力增宽修正,以期获取更高的精度。
目前国际上所采用的经典方法是大气辐射的有关理论,求解FWHM,解决压力增宽问题,有如下计算方法:
式中γL(P,T)代表洛伦兹半峰全宽,P是压强数据,T是温度数据,γ0是标准气压P0(1.013KPa)和标准温度T0(273K)时的半高宽,n为温度相关指数。
其局限性在于取值上的不确定性。对于CO2分子,γ0随谱线变化明显,国际上通用取值为0.07cm-1;温度相关指数n在0.5-1的范围内变动,取0.5时称为经典值。这种取值上的模糊性,给反演结果带来了巨大的误差。
本技术领域尚未有相关解决方案出现。
发明内容
本发明的目的是提供一种较为准确的大气压力增宽的修正方法。
为达到上述目的,本发明采用如下的技术方案:
一种差分吸收大气CO2激光雷达大气压力增宽修正方法,其特征在于,包括以下步骤:步骤1,采用以下公式计算洛伦兹半峰全宽,
式中,γL(P,T,λ)代表洛伦兹半峰全宽,P为压强数据,λ为波长,T是温度数据,T0是标准温度,bL(air,λ)为精确测定的大气压力增宽系数,bL(CO2,λ)是CO2气体的自增宽系数,χ是CO2气体同空气的体积比,na(λ)、 值分别为大气和CO2气体温度相关指数;
步骤2,判断是否满足迭代结束条件,是则结束流程,否则进入步骤3;
步骤3,根据步骤1所得洛伦兹半峰全宽,由差分吸收激光雷达方程得到CO2浓度分布,该 结果直接以分子浓度C显示,通过以下公式转化得到CO2气体的浓度与体积比χ
式中,C为每立方厘米空气中CO2分子个数;
返回步骤1,根据计算出的新的CO2气体的浓度与体积比χ重新计算洛伦兹半峰全宽。
而且,步骤2中迭代结束条件采用预设的迭代次数。
而且,步骤1中,根据CO2气体同空气的体积比χ的预设初值,采用公式(3)计算洛伦兹半峰全宽,预设的迭代次数设为1。
而且,步骤1中,CO2气体同空气的体积比χ的预设初值为0.04%。
本申请提出了一种压力增宽的新修正方法,在最大程度上实现了谱线数据同标准状况下FWHM以及温度相关指数n的对应,根据实验波长的选取(本发明采用1.57μm波长)采用对应的半高宽γ0和温度相关指数n,更为科学、合理,反演结果更为精确。同时可以解决了其他波长(谱线)数据面临的相同问题。本发明具有以下优点和积极效果:
1)摒除了经典纠正在取值上的不确定性,避免了新的误差的带入;
2)明确了标准状况下的FWHM同波长的关系;
3)根据实际场景需求,分别采用带值求解和迭代的方法具有灵活性。带值求解运算速度快,可以实现实验反演结果的实时显示;迭代的方法能够大幅度提高反演精度。
4)新纠正算法参数的测量是在广谱范围内进行的,对于其他(如2μm)激光雷达具有参考价值;
具体实施方式
以下通过实施例详细说明本发明技术方案。
实施例的一种差分吸收大气CO2激光雷达大气压力增宽修正方法可采用计算机软件计算实现自动运行,包括以下步骤:
步骤1,采用以下公式计算洛伦兹半峰全宽,
式中,γL(P,T,λ)代表洛伦兹半峰全宽,P为压强数据,λ为波长,T是温度数据,T0是标准温度,bL(air,λ)为精确测定的大气压力增宽系数,bL(CO2,λ)是CO2气体的自增宽系数, χ是CO2气体同空气的体积比,na(λ)、 值分别为大气和CO2气体温度相关指数。
本发明提出的差分吸收大气CO2激光雷达大气压力增宽修正方法,首先明确了采用公式(3)进行洛伦兹半峰全宽的计算。该公式需要随高度分布的温度、压强廓线,一共涉及8个参数,其中5个参数bL(air,λ)、bL(CO2,λ)、χ、na(λ)、 是新增的、具体的,明确了FWHM对波长的敏感性,更加复杂、精确。纠正公式中的增宽系数bL(air)、bL(CO2)及温度相关指数na、 同激光发射波长λ密切相关,一一对应。
根据CO2的吸收谱线,并考虑到相应的水汽干扰应该最小情况,本发明建议采用1.57μm波长,即差分吸收1.57μm大气CO2激光雷达选择1572.018nm的波长λon作为波长,1572.150nm作为波长λoff。根据λon或λoff相应的参数用公式(3)进行计算,可以得到λon或λoff的吸收截面积精度。
实际在计算时,由于波长λoff的吸收截面积同λon相比可忽略,往往只关注λon的相关参数值。选择λon作为λ,增宽系数bL(air,λ)、bL(CO2,λ)及温度相关指数na(λ)、 采用对应的值,采用公式(3)进行计算,结果即可作为本步骤所求洛伦兹半峰全宽。
实施例参考自然界CO2浓度的垂直分布,将CO2气体同空气的体积比χ的初值预设为0.04%。
步骤2,判断是否满足迭代结束条件,是则结束流程,否则进入步骤3。
步骤2中迭代结束条件可以采用预设的迭代次数。一般可以设为3次或4次。流程设计时,可以在第一次执行步骤1时,令当前迭代次数的初值为1,执行到步骤3时当前迭代次数=当前迭代次数+1。
步骤3,通过以下公式转化得到CO2气体的浓度与体积比χ
式中,C为每立方厘米空气中CO2分子个数;
然后返回步骤1,根据计算出的新的CO2气体的浓度与体积比χ重新计算洛伦兹半峰全宽。
本发明涉及不同高度大气中CO2分子浓度与体积比χ的转化,同时考虑到温度、压强的影响。
洛伦兹半峰全宽FWHM直接影响到差分吸收激光雷达求解CO2浓度的关键部分——吸收截面积。根据步骤1所得洛伦兹半峰全宽,由差分吸收激光雷达方程得到CO2浓度分布,差分吸收激光雷达方程即公式(1)。该结果直接以分子浓度C显示。
公式(3)是在实验室条件下进行实验并验证的,是向容器中冲入固定体积比的CO2,得到的是定值理论;在激光雷达领域,CO2的浓度随高度明显变化,这是动值问题。对此,以该定值理论为基础,本发明提出了两种解决思路:一是直接带值求解FWHM,二是采用迭代的方法求解。
具体来说,第一种思路即认为χ为0.04%(参考自然界CO2浓度的垂直分布)且为计算FWHM时的体积比真值,直接带入公式计算,可以迅速得到计算结果,同时对精度的影响不大。其优点在于,进行观测实验时可以快速、实时显示CO2浓度的垂直分布情况,便于直接查看CO2浓度结果,同时不会将错误值重复引入计算,对于激光雷达***调试工作具有很高的参考价值。步骤1中根据CO2气体同空气的体积比χ的预设初值,采用公式(3)计算的洛伦兹半峰全宽即作为本发明的结果,无需再执行后续步骤。具体实施时,可以通过将预设的迭代次数设为1实现。
第二种思路是采用迭代的方式,设定CO2体积比初值0.04%作为迭代变量,经过一次迭代可以反演出CO2浓度分布,该结果直接以分子浓度显示,通过公式(4)转化为χ,此时,χ是一条与高度有关的廓线,与温度、压强值相对应,观察公式(3)可以发现,FWHM同温度、压强、χ密切相关,三个量均是高度的函数,χ廓线可以直接应用到第二次迭代中,经过3~4次迭代可以满足精度。迭代后的反演结果精度明显优于第一种思路,但是对于数据处理***性能要求较高,计算时间延迟较大。
在实验初期,激光雷达回波信号并不理想,经常出现体积比大大超过或小于正常大气水平的情况,采用第一种方案,可快速、实时地查看结果的正确性,便于***调节;实验后期,***稳定后,在各种硬件、软件问题已经克服的前提下,误差二次引用的可能性大大降低,对于数据的精确处理必须运用迭代的方法,来提高精度。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种差分吸收大气CO2激光雷达大气压力增宽修正方法,其特征在于,包括以下步骤:
步骤1,根据CO2气体同空气的体积比χ的预设初值,采用以下公式计算洛伦兹半峰全宽,
式中,γL(P,T,λ)代表洛伦兹半峰全宽,P为压强数据,λ为波长,T是温度数据,T0是标准温度,bL(air,λ)为精确测定的大气压力增宽系数,bL(CO2,λ)是CO2气体的自增宽系数,χ是CO2气体同空气的体积比,na(λ)、值分别为大气和CO2气体温度相关指数;
步骤2,判断是否满足迭代结束条件,是则结束流程,否则进入步骤3;
步骤3,根据步骤1所得洛伦兹半峰全宽,由差分吸收激光雷达方程得到CO2浓度分布,该结果直接以分子浓度C显示,通过以下公式转化得到CO2气体同空气的体积比χ
式中,C为每立方厘米空气中CO2分子个数,P0为标准气压;
返回步骤1,根据计算出的新的CO2气体同空气的体积比χ重新计算洛伦兹半峰全宽。
2.根据权利要求1所述差分吸收大气CO2激光雷达大气压力增宽修正方法,其特征在于:步骤2中迭代结束条件采用预设的迭代次数。
3.根据权利要求2所述差分吸收大气CO2激光雷达大气压力增宽修正方法,其特征在于:步骤1中,选择1572.018nm的波长λon作为波长λ,采用公式(1)计算洛伦兹半峰全宽。
4.根据权利要求1或2或3所述差分吸收大气CO2激光雷达大气压力增宽修正方法,其特征在于:步骤1中,CO2气体同空气的体积比χ的预设初值为0.04%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210351492.XA CN102879772B (zh) | 2012-09-20 | 2012-09-20 | 一种差分吸收大气co2激光雷达大气压力增宽修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210351492.XA CN102879772B (zh) | 2012-09-20 | 2012-09-20 | 一种差分吸收大气co2激光雷达大气压力增宽修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102879772A CN102879772A (zh) | 2013-01-16 |
CN102879772B true CN102879772B (zh) | 2014-04-16 |
Family
ID=47481167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210351492.XA Active CN102879772B (zh) | 2012-09-20 | 2012-09-20 | 一种差分吸收大气co2激光雷达大气压力增宽修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102879772B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106526614B (zh) * | 2016-10-25 | 2019-04-30 | 中国科学院合肥物质科学研究院 | 一种优化激光雷达探测大气成分谱线分析的方法 |
CN107505632B (zh) * | 2017-08-18 | 2020-04-28 | 中国科学院遥感与数字地球研究所 | 一种温压廓线与切高联合反演方法 |
CN108426856B (zh) * | 2018-03-26 | 2019-09-10 | 武汉大学 | 评估激光雷达卫星测量大气co2浓度性能的综合分析*** |
CN109655843A (zh) * | 2019-01-16 | 2019-04-19 | 武汉大学 | 探测二氧化碳浓度廓线的脉冲红外差分吸收激光雷达*** |
CN111736172B (zh) * | 2020-08-24 | 2020-11-20 | 中国人民解放军国防科技大学 | 一种基于大气扰动相干激光探测的空中目标探测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6727492B1 (en) * | 2000-01-13 | 2004-04-27 | Regents Of The University Of Colorado | Cavity ringdown spectroscopy system using differential heterodyne detection |
US20110150035A1 (en) * | 2009-12-17 | 2011-06-23 | Hanson Ronald K | Non-intrusive method for sensing gas temperature and species concentration in gaseous environments |
-
2012
- 2012-09-20 CN CN201210351492.XA patent/CN102879772B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102879772A (zh) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879772B (zh) | 一种差分吸收大气co2激光雷达大气压力增宽修正方法 | |
Zhang et al. | Modified ensemble Kalman filter for nuclear accident atmospheric dispersion: Prediction improved and source estimated | |
CN103868836B (zh) | 一种同时测量大气颗粒物后向散射系数和臭氧浓度廓线的方法 | |
CN104990886A (zh) | 一种傅里叶变换红外光谱仪波数漂移的校正方法 | |
CN102680020A (zh) | 一种基于波长调制光谱技术的气体参数在线测量方法 | |
CN103175759A (zh) | 基于多种地基遥感技术获取城市气溶胶复折射指数的方法 | |
CN104062307B (zh) | 一种车载伽玛能谱大气氡影响修正方法 | |
CN106052556A (zh) | 一种基于激光追踪仪多站位测量的三坐标测量机空域坐标修正方法 | |
CN102628946B (zh) | 大气二氧化硫和臭氧廓线拉曼-瑞利/米多功能激光雷达测量装置及探测方法 | |
CN104614333B (zh) | 一种大气水汽分子吸收系数廓线测量方法 | |
US9759597B2 (en) | Methods for calibrating a fugitive emission rate measurement | |
EP3226035A1 (en) | Measurement method for air kerma conventional true value | |
US10302610B2 (en) | Generation source analyzing device and generation source analyzing method | |
CN104406715A (zh) | 一种遥感估算地表感热/潜热通量的精度评价方法及*** | |
CN106052717A (zh) | 一种利用精密弹道的电波折射修正效果标定方法 | |
Vardag et al. | Comparisons of continuous atmospheric CH 4, CO 2 and N 2 O measurements–results from a travelling instrument campaign at Mace Head | |
CN112131685A (zh) | 一种联合循环机组整体热力性能试验不确定度评估方法 | |
CN103558182B (zh) | 一种用于***体在线分析仪确定气体浓度的方法 | |
US11193917B2 (en) | Composition analysis apparatus and composition analysis method | |
CN103743679A (zh) | 一种加快温室气体垂直柱浓度反演速度的方法 | |
CN103344252B (zh) | 一种航空高光谱成像***定位误差分析方法 | |
CN104570035A (zh) | 一种放射性气体核素β射线自吸收校正方法 | |
CN106569252A (zh) | 一种校正伽玛总量测量型活性炭测氡仪一致性的方法 | |
CN102628947A (zh) | 大气臭氧廓线拉曼-瑞利/米激光雷达测量装置及探测方法 | |
CN103106332A (zh) | 一种测量不确定度的分析方法 |
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 |