基于环形临边探测中层大气臭氧浓度廓线非线性反演方法

1.本发明涉及中层大气气体成分遥感反演领域，特别涉及基于环形紫外临边探测散射辐射反演大气臭氧浓度廓线的遥感反演方法。


背景技术：

2.臭氧是大气中的关键成分之一，平流层臭氧层含有约90％的大气臭氧分子，通过阻止太阳中的有害紫外辐射来保护地球上的生命；尽管中间层大气臭氧含量很少(在50公里以上臭氧含量占大气臭氧总含量的0.3％),但是它在辐射过程中却起着重要的作用。臭氧有三个主要的吸收带：哈特莱hartley吸收带(200～300nm)、哈金斯huggins吸收带(300～360nm)和查普斯chappuis弱吸收带(430～750nm)。臭氧将吸收的太阳紫外辐射转换为热能，影响大气温度结构，从而影响大气环流。因此，研究臭氧垂直分布及变化对地球气候、区域环境以及生物圈都有着非常重要的意义。
3.臭氧测量的方法和技术手段很多，从地面的无线电探空仪、激光雷达、微波辐射计、傅里叶红外光谱仪、乌姆克尔(umkehr)测量、差分吸收光谱仪doas等等到卫星上的各种载荷，如用于获得臭氧总量的toms(mcpeters等，1996)、omi(bak等，2015)、gome系列(van der a等，1998；miles等，2015)，用于获得臭氧廓线的sbuv系列、sage系列、osiris、sciamachy、omps等。还有高空红外微波临边探测仪器mls(waters等,1993)、mipas(fischer等,2008)、tes(beer，2006)等等，在此就不一一赘述。地基仪器探测的高度有限，卫星载荷才能探测高层大气。卫星观测一般分为天底和临边两种观测方式。天底观测具有较好的水平分辨率，利用天底观测大气发射光谱也能获得臭氧廓线信息，但其中的遥感反演方法是基于先验臭氧廓线，一般在假设廓线形状不变的基础上反演，而且垂直分辨率不高。只有临边探测方式直接分层探测，才能提供较高垂直分辨率的臭氧廓线。临边探测臭氧因其具有较高的垂直分辨率和信息量大的优点，是近年来卫星探测臭氧的主要手段。遥感反演方法的提出都是基于具体的遥感探测仪器，每一个星载临边探测仪器正常运行以后，都会有针对该仪器性能和特点的遥感反演方法。现在主流的遥感反演方法是最优化方法，传统上通过高斯牛顿迭代算法来实现，但高斯牛顿迭代算法一般应用于非线性关系较弱的反演。而对于临边散射辐射和大气成分廓线之间非常复杂的强非线性关系，需要考虑强非线性的反演算法。
4.我国第一个真正意义上的空间实验室天宫二号于2016年9月15日发射，2019年7月19日离轨结束。在距地面约400公里的轨道上运行1036天，获得了大量的科学观测数据。其中临边探测仪器组(多波段前向临边光谱仪和环形紫外临边成像仪)是其中的重要组成部分，特别是环形紫外临边成像仪作为世界首创的仪器，它由三个波段(中心波长分别为265、295、360nm)组成，像元分辨率为3km，采用超大视场(360
°
环形)同时获得二维(高度和空间)观测数据，成像平面上的高度和空间位置都不相同，因此相比其它卫星载荷获得的一维高度观测数据复杂，处理困难，而从二维观测数据直接反演大气臭氧浓度廓线，得到臭氧浓度的三维分布，还没有人做过这方面的工作，数据处理方法和反演算法都值得探索研究。因为
360nm通道波段相对于265nm和295nm来说，臭氧的吸收较弱，敏感高度偏低，因此中层中部大气臭氧浓度廓线的遥感反演只考虑265nm和295nm两个通道波段。


技术实现要素：

5.为解决上述技术问题，本发明提出了基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，本方法能在同一时刻得到臭氧廓线的三维分布，对于研究气候变化和大气环境都具有重要的意义。
6.为实现上述目的，本发明提供了基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，包括：
7.对通道波段的二维临边辐射观测数据进行预处理；
8.计算所述通道波段的权重函数，获取敏感高度范围和反演高度范围；
9.基于所述敏感高度范围和反演高度范围，反演所述通道波段大气臭氧浓度廓线；
10.基于所述通道波段的反演结果，获取中层大气(35-75km)大气臭氧浓度廓线；
11.基于所述大气臭氧浓度廓线，获取大气臭氧浓度环形三维分布。
12.可选的，所述对通道波段二维临边辐射观测数据进行预处理包括：
13.对所述二维临边辐射观测数据进行筛选，获取一维临边辐射数据；
14.对所述一维临边辐射数据进行优化，采用五点平滑法获取新临边辐射数据，其中，所述五点平滑法为所述临边辐射数据中的临边辐射值使用其前后共五点平均值替代。
15.可选的，对所述二维临边辐射观测数据进行筛选包括：
16.将所述二维临边辐射观测数据按照1度间隔视方位角读取，然后按照切点高度排序，得到一维临边辐射数据；
17.将所述一维临边辐射数据，按照方位角进行排列，获取预设数量的一维临边辐射数据。
18.可选的，对所述临边辐射数据进行优化包括：
19.计算所述一维临边辐射数据的切点均值，并计算所述临边辐射数据与切点均值的相对误差；
20.剔除所述相对误差大于预设阈值的所述一维临边辐射数据。
21.可选的，反演所述通道波段大气臭氧浓度廓线包括：采用最优估计方法和levenberg-marquardt非线性反演算法，反演所述通道波段大气臭氧浓度廓线。
22.可选的，采用所述最优估计化方法包括：
23.获取大气臭氧浓度廓线；
24.基于所述大气臭氧浓度廓线，计算模式临边辐射值，并获取通道反演模式向量；
25.基于所述模式通道反演模式向量和所述通道反演观测向量，计算代价函数；
26.通过计算所述代价函数，判断所述代价函数的结果是否满足要求，其中，满足则输出大气臭氧浓度廓线，不满足则将所述大气臭氧浓度廓线输入所述levenberg-marquardt非线性反演算法。
27.可选的，所述levenberg-marquardt非线性反演算法包括：
28.将所述最优估计化方法中所述代价函数的结果不满足要求的所述大气臭氧浓度廓线输入到所述levenberg-marquardt非线性反演算法，引入非线性算子并重新计算所述
大气臭氧浓度廓线；
29.基于引入的所述非线性算子，重新计算所述代价函数；
30.通过重新计算的所述代价函数，判断所述代价函数的结果是否满足要求；
31.其中，满足则输出所述大气臭氧浓度廓线，并将所述大气臭氧浓度廓线输入到所述最优估计化方法中，重新计算所述大气臭氧浓度廓线；不满足则重新计算代价函数，直至满足要求为止。
32.可选的，所述代价函数为；
[0033][0034]
其中，j为代价函数，y为通道反演观测向量，f(x)为通道反演模式向量，x表示大气臭氧浓度廓线一维向量，j为切点高度编号，n为最大切点高度编号，yj为通道反演观测向量y中第j个切点高度的值，f(xj)为通道反演模式向量f(x)中第j个切点高度的值，xj为第j个切点高度大气臭氧浓度值。
[0035]
可选的，所述大气臭氧廓线采用迭代公式计算，
[0036]
所述迭代公式为：
[0037][0038]
其中，x
i+1
为第i+1次迭代计算的大气臭氧浓度廓线一维向量，xi为第i次迭代计算的大气臭氧浓度廓线一维向量，γ为levenberg-marquardt非线性算子，sa为先验臭氧廓线协方差矩阵，ki为第i次迭代计算的权重函数矩阵，为ki的转置矩阵，s
ε
为观测噪声协方差矩阵，y为通道反演观测向量，f(xi)为第i次迭代计算的通道反演模式向量，xa为初始假设的先验大气臭氧浓度廓线一维向量。
[0039]
与现有技术相比，本发明具有如下优点和技术效果：
[0040]
本发明基于天宫二号上搭载的环形紫外临边成像仪二维探测数据，经过环形临边辐射数据的预处理、采用最优化方法和levenberg-marquardt非线性反演算法得到两个通道波段的大气臭氧浓度廓线、综合两个通道波段的反演结果得到35-75km的大气臭氧浓度廓线、推广到整个探测范围内的点得到同时观测的大气臭氧浓度的环形三维分布。相比于现有方法，只能得到某一个时刻的臭氧廓线，这是第一次能在同一时刻得到臭氧廓线的三维分布。对于研究气候变化和大气环境都具有重要的意义。也为未来紫外-可见-红外临边成像仪的研制及应用研究提供了技术支持。
附图说明
[0041]
构成本技术的一部分的附图用来提供对本技术的进一步理解，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。在附图中：
[0042]
图1为本发明实施例的环形临边探测中层中部大气臭氧浓度廓线非线性反演方法的流程示意图。
具体实施方式
[0043]
需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本技术。
[0044]
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
[0045]
实施例
[0046]
本实施例提出了基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，如图1所示，包括：
[0047]
1.原始观测数据的预处理，即从二维观测数据处理出一维高度数据，如图1所示，包括以下几个步骤：
[0048]
1.1分别针对两个通道波段(265nm和295nm)二维临边辐射观测数据，按照1度间隔视方位角读取数据，一般有15组数据，每组数据包含几百个切点高度，按照切点高度从小到大排列，得到15条一维临边辐射数据。然后按照方位角大小排列，考虑到边缘数据受仪器噪声影响较大，只选择中间的11条临边辐射数据。
[0049]
1.2对上述11条临边辐射数据，首先排除异常点，即沿切点高度方向，上下两点均值的相对误差大于50％，则排除该点。如：
[0050]
对于第j个切点高度的临边辐射rj，分别计算其上下两点均值和相对误差σj：
[0051][0052][0053]
其中，为临边辐射数据的第j个切点均值，r
j-1
为第j-1个切点高度的临边辐射值，r
j+1
为第j+1个切点高度的临边辐射值，σj为第j个切点高度临边辐射数据的相对误差，rj为第j个切点高度的临边辐射值。
[0054]
如果σj》50％，则删除第j个切点高度的临边辐射值，然后重新排列临边辐射数据。
[0055]
其次对上述有n个切点高度(1,2，
…
，n)的临边辐射值采用五点平滑方法，第j个切点高度的临边辐射值用其前后五点平均值取代，即
[0056][0057]
针对中层大气选择切点高度在[30,80]km范围内平滑后的一维随切点高度变化的临边辐射数据用于反演。
[0058]
2.由辐射传输模式sciatran分别计算每个通道波段权重函数，根据权重函数最大值所在的高度为该通道的敏感高度，权重函数值中大于(20％
×
最大值)的高度范围为反演高度范围，确定敏感高度和反演高度范围，得到大气臭氧在265nm通道波段的敏感高度约为60km，可反演45-75km高度范围的廓线，而在295nm通道波段的敏感高度约为50km，可反演35-65km高度范围的廓线。
[0059]
3.采用最优化估计方法和levenberg-marquardt非线性反演算法分别反演单通道
波段大气臭氧浓度廓线，包括以下几个步骤：
[0060]
3.1定义单通道反演向量y＝[y1,y2,
…
,yn]：
[0061][0062]
其中，j＝1,2
…
n代表第j个切点高度，n为可反演的最大切点高度，rj为第j个切点高度的临边辐射值，r0为参考高度的临边辐射值，yj为反演向量y中第j个切点高度的值。
[0063]
3.2采用最小平方函数作为代价函数j，最优化方法就是最小化代价函数，即使观测得到的反演向量与模式计算的临边辐射向量的差值平方最小。
[0064][0065]
其中，j为代价函数，y为通道反演观测向量，f(x)为通道反演模式向量，x表示大气臭氧浓度廓线一维向量，j为切点高度编号，n为最大切点高度编号，yj为通道反演观测向量y中第j个切点高度的值，f(xj)为通道反演模式向量f(x)中第j个切点高度的值，xj为第j个切点高度大气臭氧浓度值。
[0066]
3.3采用levenberg-marquardt非线性反演算法计算大气臭氧浓度廓线，迭代公式为：
[0067][0068]
其中，x
i+1
为第i+1次迭代计算的大气臭氧浓度廓线一维向量，xi为第i次迭代计算的大气臭氧浓度廓线一维向量，γ为levenberg-marquardt非线性算子，sa为先验臭氧廓线协方差矩阵，ki为第i次迭代计算的权重函数矩阵，为ki的转置矩阵，s
ε
为观测噪声协方差矩阵，y为通道反演观测向量，f(xi)为第i次迭代计算的通道反演模式向量，xa为初始假设的先验大气臭氧浓度廓线一维向量。
[0069]
因为引入了非线性算子，该算法相比于传统的高斯牛顿迭代算法能更好地描述临边散射辐射和大气臭氧浓度之间的非线性关系。γ非线性算子是一个变化的量，先给定一个小值γ0＝10-5
，计算初始代价函数j0，然后采用线性迭代方法更新γ值，计算公式如下：
[0070]
γ
k+1
＝γk×
ratio_factor
[0071]
其中，γ
k+1
表示第k+1次迭代，γk表示第k次迭代，ratio_factor表示比例因子，是固定值，可取10。每次更新γ值，重新计算代价函数，如果更新后的代价函数比更新前小，则继续更新，否则停止。
[0072]
4.综合两个通道波段的反演结果，即相同的反演高度求平均值，最终得到中层大气(35-75km)大气臭氧浓度廓线。
[0073]
5.对环形紫外临边成像仪一次成像探测范围内的数据均采用上述方法处理，就可以同时得到地面多个点的大气臭氧浓度廓线，形成大气臭氧浓度的环形三维分布。
[0074]
以上，仅为本技术较佳的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应该以权利要求的保护范围为准。

技术特征：
1.基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，包括：对中层大气通道波段的二维临边辐射观测数据进行预处理；计算所述通道波段的权重函数，获取敏感高度范围和反演高度范围；基于所述敏感高度范围和所述反演高度范围，反演所述通道波段的大气臭氧浓度廓线；基于所述通道波段的反演结果，获取中层大气臭氧浓度廓线；基于所述大气臭氧浓度廓线，获取大气臭氧浓度环形三维分布。2.根据权利要求1所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，对通道波段的所述二维临边辐射观测数据进行预处理包括：对所述二维临边辐射观测数据进行筛选，获取一维临边辐射数据；对所述一维临边辐射数据进行优化，采用五点平滑法获取新临边辐射数据，其中，所述五点平滑法为所述临边辐射数据中的临边辐射值使用其前后共五点平均值替代。3.根据权利要求2所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，对所述二维临边辐射观测数据进行筛选包括：将所述二维临边辐射观测数据按照1度间隔视方位角读取，然后按照切点高度排序，得到一维临边辐射数据；将所述一维临边辐射数据，按照方位角进行排列，获取预设数量的所述一维临边辐射数据。4.根据权利要求2所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，对所述一维临边辐射数据进行优化包括：计算所述一维临边辐射数据的上下两切点均值，并计算所述临边辐射数据与切点均值的相对误差；剔除所述相对误差大于预设阈值的所述一维临边辐射数据。5.根据权利要求1所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，反演所述通道波段大气臭氧浓度廓线包括：采用最优估计方法和levenberg-marquardt非线性反演算法，反演所述通道波段大气臭氧浓度廓线。6.根据权利要求5所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，采用所述最优估计化方法包括：获取大气臭氧浓度廓线；基于所述大气臭氧浓度廓线，计算模式临边辐射值，并获取通道反演模式向量；基于所述通道反演模式向量和所述通道反演观测向量，计算代价函数；通过计算所述代价函数，判断所述代价函数的结果是否满足要求，其中，满足则输出大气臭氧浓度廓线，不满足则将所述大气臭氧浓度廓线输入所述levenberg-marquardt非线性反演算法。7.根据权利要求6所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，所述levenberg-marquardt非线性反演算法包括：将所述最优估计化方法中所述代价函数的结果不满足要求的所述大气臭氧浓度廓线输入到所述levenberg-marquardt非线性反演算法，引入非线性算子并重新计算所述大气臭氧浓度廓线；
基于引入的所述非线性算子，重新计算所述代价函数；通过重新计算的所述代价函数，判断所述代价函数的结果是否满足要求；其中，满足则输出所述大气臭氧浓度廓线，并将所述大气臭氧浓度廓线输入到所述最优估计化方法中，重新计算所述大气臭氧浓度廓线；不满足则重新计算代价函数，直至满足要求为止。8.根据权利要求6所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，所述代价函数为；其中，j为代价函数，y为通道反演观测向量，f(x)为通道反演模式向量，x表示大气臭氧浓度廓线一维向量，j为切点高度编号，n为最大切点高度编号，y
j
为通道反演观测向量y中第j个切点高度的值，f(x
j
)为通道反演模式向量f(x)中第j个切点高度的值，x
j
为第j个切点高度大气臭氧浓度值。9.根据权利要求7所述的基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，其特征在于，所述大气臭氧浓度廓线采用迭代公式计算，所述迭代公式为：其中，x
i+1
为第i+1次迭代计算的大气臭氧浓度廓线一维向量，x
i
为第i次迭代计算的大气臭氧浓度廓线一维向量，γ为levenberg-marquardt非线性算子，s
a
为先验臭氧廓线协方差矩阵，k
i
为第i次迭代计算的权重函数矩阵，为k
i
的转置矩阵，s
ε
为观测噪声协方差矩阵，y为通道反演观测向量，f(x
i
)为第i次迭代计算的通道反演模式向量，x
a
为初始假设的先验大气臭氧浓度廓线一维向量。

技术总结
本发明公开了基于环形临边探测中层大气臭氧浓度廓线非线性反演方法，包括：对通道波段的观测数据进行预处理；计算所述通道波段的权重函数，获取敏感高度范围和反演高度范围；基于所述敏感高度范围和反演高度范围，反演所述通道波段大气臭氧浓度廓线；基于所述通道波段的反演结果，获取中层大气(35-75km)大气臭氧浓度廓线；基于所述大气臭氧浓度廓线，获取大气臭氧浓度环形三维分布。本发明能在同一时刻得到臭氧廓线的三维分布。对于研究气候变化和大气环境都具有重要的意义。和大气环境都具有重要的意义。和大气环境都具有重要的意义。


技术研发人员：宗雪梅 章卫星 段民征 吕达仁
受保护的技术使用者：中国科学院大气物理研究所
技术研发日：2023.04.14
技术公布日：2023/7/13