中国激光, 2019, 46 (8): 0810001, 网络出版: 2019-08-13   

LEO-LEO红外激光掩星CO2浓度测量技术研究 下载: 1024次

LEO-LEO Infrared Laser Occultation Technique to Measure Atmospheric Carbon Dioxide Concentration
作者单位
1 中国科学院上海光学精密机械研究所空间激光信息传输与探测技术重点实验室, 上海 201800
2 中国科学院大学材料与光电研究中心, 北京 100049
摘要
利用基于近地轨道(LEO)卫星组网的红外激光掩星(LIO)技术,可对地球大气温室气体垂直廓线进行主动探测,为全球温室气体浓度廓线的测量提供新手段。介绍了基于LIO技术的大气温室气体浓度廓线的测量原理,针对最主要的温室气体CO2建立了LIO信号链路模型,通过仿真分析了工作波数对CO2探测精度的影响,并基于误差最小的仿真结果对CO2探测波数进行优化选择。采用所选波数对LIO技术的探测性能进行分析,最终得到可用于CO2浓度廓线探测的波数为4771.6215 cm -1和4772.0240 cm -1。在5~35 km高度的有效探测范围内,实现了0.6~1.4 km的垂直分辨率,CO2浓度探测的相对随机误差小于0.8%,最小相对随机误差(0.229%)出现在10 km处。研究结果为星载LIO大气探测系统原理样机的设计提供了重要参考。
Abstract
The infrared laser occultation (LIO) technique based on the low earth orbit satellite network can perform an active detection for vertical profiles of greenhouse gases in the earth's atmospheric, thereby providing a new method for global greenhouse gas concentration profile measurement. In this paper, the principle of the concentration profile measurement of the LIO technique is introduced. The influence of the operation wavelength on the detection accuracy of the most important greenhouse gas CO2 is simulated by establishing the LIO signal link model. Subsequently, the detection wavelengths for CO2 are optimized based on the simulation results with the minimum error. Additionally, the detection performance of the LIO technique is analyzed using the selected wavelengths. Finally, the wavenumbers of 4771.6215 cm -1 and 4772.0240 cm -1 available for the CO2 concentration profile detection are selected. The selected wavelengths can achieve a vertical resolution of 0.6-1.4 km in the altitude range of 5-35 km. The relative random error of profile detection is lower than 0.8%. The minimum relative random error (0.229%) appears at 10 km. The research results can provide an important reference for the design of a prototype of spaceborne LIO atmospheric detection system.

1 引言

自工业革命以来,人类活动加剧,随之而来的酸雨、温室效应等环境问题逐渐显露并日趋严峻。研究表明,以CO2为主的温室气体是温室效应的主要推力,大气中CO2浓度的变化对气候的影响备受气候、环境等自然学科的重视。而高时空分辨率的CO2浓度数据是气候科学研究的重要需求,也是研究大气全球变化、人类活动对大气环境影响的重要参数,对于提高数值天气预报的精度以及研究日地能量传输过程具有重要意义。

星载积分路径差分吸收(IPDA)激光雷达具有探测精度高、测量区域广且不受气溶胶干扰等优点[1],刘继桥课题组[2-3]针对CO2 IPDA激光雷达开展了广泛研究,目前已研制出原理样机,其对CO2体积混合比柱浓度的理论测量误差可达到10-6。但是IPDA激光雷达获取的是整层大气CO2柱浓度的平均值,无法提供大气CO2的垂直变化数据。为了兼顾全球CO2浓度的测量精度和垂直分辨能力,可以考虑将IPDA激光雷达技术应用到新的大气探测方式中。通过卫星组网构成掩星星座,对地球大气进行临边探测,是获得大气参量廓线的有效手段。目前利用无线电信号进行探测的全球导航卫星系统无线电掩星技术GNSS RO(Global Navigation Satellite System Radio Occultation,GRO),已实现对大气热力学参量廓线的测量,其数据产品已被广泛应用于数值天气预报和全球气候监测[4]。鉴于GRO技术在大气热力学参量廓线测量中的成功运用,欧洲的ACCURATE(Atmospheric Climate and Chemistry in the UTLS Region And Climate Trends Explorer)[5]提出了在近地轨道(LEO)卫星构成的掩星星座上搭载红外激光发射机和接收机的设想,采用主动临边探测的方式进行大气成分和风速的测量,即激光差分吸收与掩星探测相结合的红外激光掩星(LIO)技术。Clifford等[6]也提出了类似的方案,即将搭载有激光发射机的卫星与多个空间反射镜组网进行临边探测,并针对对流层顶到平流层底高度(UTLS, 5~35 km)区域水汽的探测进行仿真研究。

LIO技术结合了激光差分吸收与掩星探测的优点,主动式激光差分吸收光谱技术可达到较高的探测精度,同时LEO-LEO掩星探测方式可以实现对全球高覆盖率的测量,获得垂直分辨能力(约1 km)的探测数据。本文主要研究了利用LIO技术测量全球CO2浓度廓线的方法。建立针对UTLS区域CO2浓度廓线的探测仿真模型,通过计算分析优化选择激光的工作波数,利用选择的工作波数对LIO技术的探测性能进行仿真分析。

2 基本原理

图1为LIO方法的几何结构图,αarθΦ分别为激光路径弯曲角、碰撞参数、卫星到地心的距离、两颗卫星与地心连线的夹角以及激光传输方向与地心连线之间的夹角,Ve分别表示卫星速度矢量和激光传输方向矢量,下标T、R分别表示发射机和接收机。如图1所示,两颗分别搭载红外激光发射机(Tx)和接收机(Rx)的LEO卫星构成掩星星座。在一定几何位置处,发射机发射两束波长相近的激光,波长分别位于CO2吸收谱的强吸收(on-line)和弱吸收(off-line)位置,激光信号穿过地球大气层被接收机接收,然后在两颗卫星相对运动过程中形成下降(或上升)掩星事件,实现对大气层从上到下(或从下到上)近似垂直的扫描测量,获得剖面信息。利用光学积分路径差分吸收原理和Abel积分变换实现CO2浓度廓线的反演。

图 1. LIO几何结构示意图

Fig. 1. Schematic of LIO geometric structure

下载图片 查看所有图片

2.1 激光路径的确定

图1所示,在掩星事件中,激光信号穿过大气层时由于大气折射率的不均匀分布,将发生α角度的弯曲,定义碰撞参数a和切点半径r0。在大气折射率局部球对称的近似条件下对大气分层进行光线追踪[7],αar0满足

a=μ(r0)r0,(1)

α(a)=2ar0rtopdlnμdrdr[(μr)2-a2]12,(2)

式中:μ(r0)为r0位置处的折射率;rtop为积分上限,由于高层大气十分稀薄,对积分的贡献可以忽略,在这里将rtop取为60 km。

在发生掩星事件时,两颗LEO卫星的位置坐标rTrR与弯曲角α存在特殊的几何关系,即

α(a)=θ-arccosarT-arccosarR(3)

结合(1)~(3)式可实时对αar0进行求解。

2.2 浓度廓线的反演

on-line和off-line激光在穿过大气层时被CO2分子吸收,传输路径上的差分吸收光学厚度(DAOD)表示为

DAOD=ρCO2(σon-σoff)ds,(4)

式中:s为路径积分的积分变量,表示长度;ρCO2为CO2的分子数密度;σonσoff分别为CO2在on-line和off-line波数处的吸收截面。结合(2)式,切点高度为r0时传输路径上的DAOD表示为

DAOD(r0)=2r0rtopΔkμr(μ2r2-μ02r02)12dr,(5)

式中:Δk=ρCO2(σon-σoff);μ为折射率。

对(5)式应用Abel积分变化可以得到

Δk(r0)=-1πa'(r0)r0rtopdDAODdrμr(μ2r2-μ02r02)12dr,(6)

根据Beer-Lambert定律,DAOD又可以用发射光强和透射光强表示为

DAOD=lnP0(λon)P(λoff)P0(λoff)P(λon),(7)

式中:P0为激光信号的发射功率;λon为吸收通道的工作波长;P为接收到的激光信号功率;λoff为参考通道的工作波长。通过(7)式求得DAOD的测量值,将其代入(6)式中,再结合理想气体状态方程,可以得到CO2体积分数廓线为

χCO2(r0)=Δk(r0)kBT(σon-σoff)p·106,(8)

式中:χCO2为CO2的体积分数;kB为玻尔兹曼常数;T为温度;p为大气压强。

2.3 链路衰减模型

激光信号从发射机到接收机整个传输过程中的损耗主要由三部分组成:自由空间传播损耗、大气损耗和接收机光学系统损耗。由激光雷达方程可知,接收机接收到的激光信号功率为

PR=PTdRDTRφ2exp(-τ)η,(9)

式中:PRPT分别为接收功率和发射功率;dRDTRφ分别为接收望远镜直径、激光传输距离和激光发散角;τ为大气衰减对应的光学厚度;η为接收机光学系统的效率。

大气衰减主要包含CO2吸收、气溶胶吸收、其他无关气体(H2O、CO2、O3、N2O、CO、CH4、O2等)吸收、大气分子瑞利散射。瑞利散射产生的等效吸收截面为

σR=32π3(μ-1)23λ4Nair,(10)

式中:λ为波长;Nair为空气的分子数密度。

计算过程中的大气温度、湿度、压强采用美国1976标准大气模型[8]获得,红外波段的折射率采用Edlén方程[9]获得。利用HITRAN 2016数据库拟合分子的吸收截面,结合(5)式对光学厚度进行计算。

3 激光波数的优化选择

3.1 on-line波数选择

目前,学术领域针对CO2浓度的主动探测主要集中在1.5 μm和2 μm的红外波段。在临边探测方式下,考虑到激光信号穿过大气层的路径较长,在1.5 μm波段处,空气分子的瑞利散射和气溶胶米散射对信号的影响很大,而2 μm波段处受到的太阳辐射和大气热辐射的影响较小,因此拟将2 μm作为CO2浓度廓线探测的优选波段[5]

off-line处的光学厚度很小,DAOD主要由on-line光学厚度决定。如果所选择的on-line处的CO2吸收过强,则当掩星切点高度较低时,on-line信号就会因穿过大气的路程较长而具有较大的光学厚度,从而被强烈吸收,难以测量;相反,如果on-line处的吸收过弱,则当掩星切点高度较大时,on-line信号穿过大气的路程较小,并且由于高空中的CO2含量减少,导致光学厚度太小,到达接收机的on-line和off-line信号的差别就比较微弱,不利于减小测量误差。所以,应合理选择on-line波数,使UTLS高度范围(5~35 km)具有较高的探测精度。在切点高度r0处的相对随机误差(RRE)为[6]

RRE=ΔDAODDAODgAbel=1DAODN1RSNRgAbel,(11)

式中:ΔDAOD为DAOD的测量误差;N为单位垂直分辨率上的探测次数;RSNR为接收信号的信噪比;gAbel为对测量数据进行Abel变换时引入的误差放大因子[10]

利用2.3节中所述链路衰减模型,分别对不同on-line波数下系统的相对随机误差进行仿真计算,得到了5 km和35 km高度处相对随机误差与波数的关系曲线,如图2所示。

图 2. 5 km和35 km高度处相对随机误差随on-line波数的变化

Fig. 2. RRE as a function of on-line wavenumber at altitudes of 5 km and 35 km

下载图片 查看所有图片

图2可知,与5 km高度相比,on-line波数对35 km高度的相对随机误差影响更大,这是因为高海拔处CO2吸收谱线的线宽较窄,当所选用的on-line波数偏离吸收谱线中央时,CO2光学厚度显著减小,导致差分效果减弱,最终使相对随机误差增大。波数对5 km高度处相对随机误差的影响比较弱,故在选择on-line波数时,优先考虑使35 km高度相对随机误差最小的波数。最终选定on-line波数为4771.6215 cm-1,5 km和35 km高度处CO2浓度测量的相对误差分别为0.331%和0.768%。

3.2 off-line波数选择

为了达到良好的差分效果,off-line波数应选在CO2吸收尽可能弱的位置。为了消除大气宽带效应的影响,同时降低光源研制的技术难度,off-line与on-line之间的波数间隔不能过大。除此之外,还应考虑水汽对on-line和off-line的影响,在选择off-line波数时,应尽可能使on-line和off-line处水汽的光学厚度相近。由于大气中的水汽含量随海拔高度的升高而迅速减小,对CO2测量的干扰主要在低空,仿真计算切点高度为5 km时,激光信号传输路径上CO2和H2O在4771 cm-1附近的吸收光谱如图3所示。

图 3. 5 km高度处CO2和H2O的吸收光谱

Fig. 3. Absorption spectra of CO2 and H2O at altitude of 5 km

下载图片 查看所有图片

由仿真结果可知,在选定的on-line波数处,H2O的光学厚度为0.013,将产生约1.5%的浓度测量误差,影响较大。可以将off-line选在与on-line H2O光学厚度接近的4772.0240 cm-1波数处,以减小水汽的干扰。考虑到在不同高度处温度、压强的差异,水汽在不同波数处的吸收谱线将有不同程度的展宽,故5~35 km各高度上水汽对CO2差分光学厚度的干扰有所差异。水汽对CO2差分光学厚度测量结果造成的误差可表示

RSE=DAOD(H2ODAODCO2)(12)

对所选定的on-line(4771.6215 cm-1)和off-line(4772.0240 cm-1),仿真计算5~35 km不同掩星切点高度处水汽干扰产生的相对系统误差(RSE)如图4所示。结果显示,对于on-line和off-line优选波数,水汽对CO2差分光学厚度测量的干扰主要发生在5~10 km范围,且产生的误差保持在0.15%以内。可见,优选波数较大程度地减弱了水汽干扰造成的影响。

图 4. 水汽干扰产生的相对系统误差随高度的变化

Fig. 4. RSE caused by H2O as a function of altitude

下载图片 查看所有图片

4 探测性能分析

仿真中所采用的激光掩星系统参数和LEO-LEO轨道参数设置分别如表1表2所示。

表 1. 系统仿真参数

Table 1. System simulation parameters

Instrument key parameterValue
Transmitter
Pulse energy@2.09 μm /mJ1.5
Pulse length τ /ms1.5
Laser beam divergence θ /mrad0.5
Receiver
Telescope diameter D /m0.36
System optical efficiency η0.40
Detector(PIN,G12183-203k)
Nominal gain M1
Noise factor F1
Responsivity@2.09 μm R /(A·W-1)1.1
Dark current Id /nA85
Amplifier bandwidth B /kHz3

查看所有表

表 2. LEO-LEO轨道参数

Table 2. LEO-LEO orbital parameters

Orbit parameterLEO(Tx)LEO(Rx)
Apogee altitude /km500650
Perigee altitude /km500650
Inclination /(°)9090
Argument of perigee /(°)9090
Right ascension ofascending node /(°)1800
Mean anomaly /(°)270270

查看所有表

利用卫星工具包STK生成卫星位置坐标,结合(1)~(3)式仿真计算了一次掩星事件中,激光信号扫描高度与时间的关系以及不同高度下的扫描速度,结果如图5所示。

图 5. 掩星事件中激光信号的扫描特点。(a)切点高度随时间的变化; (b)不同高度下的垂直扫描速度

Fig. 5. Scanning characteristics of laser signals in occultation events. (a) Tangent point altitude as a function of time; (b) vertical scanning speed at different altitudes

下载图片 查看所有图片

图5可知,激光信号对5~35 km高度扫描持续时间约15 s,扫描速度为1.2~2.8 km/s,对于50 Hz的脉冲发射频率,25次脉冲累计次数,将获得0.6~1.4 km的垂直探测分辨率。对所选定的CO2 on-line(4771.6215 cm-1)和off-line(4772.0240 cm-1)进行波数探测,图6给出了不同高度处CO2浓度测量的随机误差的仿真计算结果。计算结果表明,随机误差仿真结果随海拔高度上升,呈先减后增的变化趋势,这是DAOD随海拔高度上升由过大变化到过小的缘故,与3.1节中的分析一致。随机误差的最小值(0.229%)出现在约10 km高度处,并且在5~35 km区域内均小于0.8%,若取CO2的体积分数为4×10-4,则35 km高度处对应的最大绝对误差为3.2×10-6,垂直分辨率为1.4 km,10 km高度附近的最小绝对误差约为9.16×10-7,垂直分辨率为0.76 km,显示出较高的测量精度和剖面分辨能力。

图 6. 5~35 km不同高度处测量得到的CO2浓度的随机误差

Fig. 6. RRE of CO2 concentration profile at different altitudes from 5 km to 35 km

下载图片 查看所有图片

图7(a)~(c)分别给出了在一个月、一季度、一年时间内掩星事件在全球分布的情况。图中正三角、倒三角分别表示上升掩星事件和下降掩星事件。由图7可知,一个月、一个季度、一年时间内发生的掩星事件数量分别为1797、5511和21866,平均60次/d,且均匀地分布在地球各个经纬度地区,水平测量分辨率为185~267 km。可见,采用LIO技术进行大气CO2浓度测量具有较高的全球覆盖率。

图 7. 掩星事件全球分布图。(a)一个月; (b)一季度; (c)一年

Fig. 7. Global distributions of occultation events. (a) One month; (b) one quarter; (c) one year

下载图片 查看所有图片

5 结论

红外激光掩星技术是一种重要的全球大气CO2浓度测量手段,具有高浓度测量精度和高全球覆盖率的优点,同时具备垂直分辨能力,可以成为IPDA大气柱浓度测量方法的重要补充。通过对CO2浓度测量随机误差进行理论分析和仿真计算,得到了UTLS区域边界CO2廓线测量随机误差随工作波数的变化,基于UTLS区域整体测量随机误差水平最小和水汽干扰误差最小的原则,选定on-line和off-line波数分别为4771.6215 cm-1和4772.0240 cm-1。探测性能仿真结果显示,所选工作波数对UTLS区域CO2浓度测量的垂直分辨率为0.6~1.4 km,相对测量误差小于0.8%,最高探测精度在10 km高度附近可达到0.229%(对应的体积分数误差为9.16×10-7)。研究结果表明,LIO技术在CO2浓度廓线探测方面具有广阔应用前景,对开展相关设计具有重要的参考价值。

参考文献

[1] 胡顺星, 陈亚峰, 刘秋武, 等. 差分吸收激光雷达系统探测背景大气SO2和NO2[J]. 中国激光, 2018, 45(9): 0911009.

    Hu S X, Chen Y F, Liu Q W, et al. Differential absorption lidar system for background atmospheric SO2 and NO2 measurements[J]. Chinese Journal of Lasers, 2018, 45(9): 0911009.

[2] 谢杨易, 刘继桥, 姜佳欣, 等. 使CO2浓度测量误差减小的星载激光雷达波长优化[J]. 红外与激光工程, 2014, 43(1): 88-93.

    Xie Y Y, Liu J Q, Jiang J X, et al. Wavelengths optimization to decrease error for a space-borne lidar measuring CO2 concentration[J]. Infrared and Laser Engineering, 2014, 43(1): 88-93.

[3] 穆永吉, 李蕊, 万渊, 等. 星载激光雷达系统杂散光分析与抑制[J]. 中国激光, 2018, 45(5): 0510005.

    Mu Y J, Li R, Wan Y, et al. Stray light analysis and suppression for spaceborne lidar system[J]. Chinese Journal of Lasers, 2018, 45(5): 0510005.

[4] 王树志, 朱光武, 白伟华, 等. 风云三号C星全球导航卫星掩星探测仪首次实现北斗掩星探测[J]. 物理学报, 2015, 64(8): 089301.

    Wang S Z, Zhu G W, Bai W H, et al. For the first time fengyun3 C satellite-global navigation satellite system occultation sounder achieved spaceborne Bei Dou system radio occultation[J]. Acta Physica Sinica, 2015, 64(8): 089301.

[5] ProschekV, KirchengastG, SchweitzerS, et al. ACCURATE: greenhouse gas profiles retrieval from combined IR-laser and microwave occultation measurements[C]∥EGU General Assembly Conference 2010, May 2-7, 2010, Vienna, Austria. [S.l.: s.n.], 2010: 14352

[6] CliffordD, HoffmannA, WeitnauerC, et al. A novel mission concept for upper air water vapour observations: active limb sounding with a constellation of retroreflectors[C]∥American Geophysical Union, Fall Meeting 2011. [S.l.: s.n.], 2011: A21D- 0115.

[7] Fjeldbo G, Kliore A J, Eshleman V R. The neutral atmosphere of Venus as studied with the mariner V radio occultation experiments[J]. The Astronomical Journal, 1971, 76: 123-140.

[8] Anderson GP, Clough SA, Kneizys FX,et al. AFGL atmospheric constituent profiles (0-120 km)[R]. [S.l.: s.n.], 1986.

[9] Edlén B. The refractive index of air[J]. Metrologia, 1966, 2(2): 71-80.

[10] Sofieva VF, KyröläE. Abel integral inversion in occultation measurements[M] ∥Kirchengast G, Foelsche U, Steiner A K. Occultations for probing atmosphere and climate. Berlin, Heidelberg: Springer, 2004: 77- 85.

李文冬, 刘继桥, 朱亚丹, 陈卫标. LEO-LEO红外激光掩星CO2浓度测量技术研究[J]. 中国激光, 2019, 46(8): 0810001. Wendong Li, Jiqiao Liu, Yadan Zhu, Weibiao Chen. LEO-LEO Infrared Laser Occultation Technique to Measure Atmospheric Carbon Dioxide Concentration[J]. Chinese Journal of Lasers, 2019, 46(8): 0810001.

本文已被 8 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!