基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法制造技术

本发明专利技术公开了基于MAX‑DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，基于朗伯‑比尔定律，使用QDOAS计算紫外波段和可见波段的O4的差分斜柱浓度及二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等痕量气体的差分斜柱浓度；利用不同波段的O4 DSCD反演对应波段的气溶胶廓线；提取反演的气溶胶光学厚度与太阳光度计实测的气溶胶光学厚度做对比验证及提取痕量气体廓线底层浓度小时平均后的结果与国控站点地面对应气体浓度做对比验证。本发明专利技术解决地基被动遥感观测垂直廓线反演算法的难题，实现对大气气溶胶及污染气体时空变化特征进行快速、连续、高效监测，满足环境污染研究和监测的需要。

Vertical Profile Inversion of Aerosols and Trace Polluted Gases Based on MAX-DOAS

The invention discloses a vertical profile inversion algorithm for aerosols and trace contaminated gases based on MAX DOAS, and calculates the differential inclined column concentration of O4 in ultraviolet and visible bands and the differential inclined column concentration of trace gases such as nitrogen dioxide, sulfur dioxide, formaldehyde, nitrite and glyoxal using QDOAS based on Lambert Beer's law. The aerosol profiles in corresponding bands are retrieved by O4 DSCD in different bands, and the retrieved aerosol optical thickness is compared with the measured aerosol optical thickness by solar photometer, and the hourly average of the bottom concentration of the extracted trace gas profiles is compared with the corresponding gas concentration on the ground of national control stations. The invention solves the difficult problem of inversion algorithm of vertical profiles for passive ground-based remote sensing observation, realizes rapid, continuous and efficient monitoring of space-time variation characteristics of atmospheric aerosols and polluted gases, and meets the needs of environmental pollution research and monitoring.

【技术实现步骤摘要】
基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法
本专利技术涉及到环境科学、地基遥感
，特别涉及基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法。
技术介绍
随着环境问题的日益加重和国际化趋势的日益明显，解决环境问题的必要性和迫切性也越来越广泛地被世界各国所认识。地基遥感相对于传统的地面采样分析技术具有其快速、大信息量、实时性强，无需观测人员值守的特点在环境监测领域具有无可比拟的优势。而目前我国在区域大气污染观测方面开展的工作，主体是以地面采样为主，辅以探空气球等对区域大气的污染状况进行观测，缺少对大气垂直方向污染状况的实时观测。同时，我国又迫切需要获得国内重点区域在近地面及垂直方向上的污染状况。地基遥感网络数据既可以实现对地面及对流层柱总量气溶胶和痕量气体的观测，又可以实现对气溶胶及各种污染气体的实时垂直观测。鉴于目前我国在地基遥感方面主要以对流层柱总量的反演为主，缺乏合理有效的垂直廓线反演算法。
技术实现思路
专利技术的目的在于提供基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，以解决地基被动遥感观测垂直廓线反演算法的难题，拓展地基MAX-DOAS的应用范围，实现对大气气溶胶及污染气体(二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等)时空变化特征进行快速、连续、高效监测，满足我国环境污染研究和环境监测的需要，以解决上述
技术介绍
中提出的问题。为实现上述目的，本专利技术提供如下技术方案：基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，包括如下步骤：步骤一：基于朗伯-比尔定律，使用QDOAS计算紫外波段和可见波段O4的差分斜柱浓度及二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等痕量气体的差分斜柱浓度，朗伯-比尔定律：步骤二：利用不同波段的O4DSCD反演对应波段的气溶胶廓线，基于最优估计算法结合前向模型SCIATRAN反演出气溶胶的垂直廓线，迭代反演算法公式：反演过程中实测结果与前向模型模拟结果拟合程度用最小价值函数来评估，函数式：步骤三：将步骤二中反演的气溶胶廓线作为输入量结合实测的痕量气体差分斜柱浓度，反演痕量气体廓线，反演的迭代过程与步骤二中气溶胶的反演相同；步骤四：提取反演的气溶胶光学厚度与太阳光度计实测的气溶胶光学厚度做对比验证，提取痕量气体廓线最底层的浓度值做小时平均后与国控站点的实测地面气体浓度做对比验证，验证过程包括数据绝对值的对比及相关性的验证。优选的，步骤二和步骤三中，基于遥感实测的气溶胶和痕量气体差分斜柱浓度反演及实时获得气溶胶和痕量气体的垂直廓线，作用于分析气溶胶及痕量气体在垂直方向上的变化规律。优选的，步骤一的公式中，I(λ,L)和I0(λ)分别为测量和初始太阳散射光谱，σj(λ)为不同气体的吸收截面，cj(l)为气体浓度，εM(λ,l)和εR(λ,l)分别为米散射和瑞利散射的消光系数，L为光程长度。优选的，步骤一测得不同气体的差分斜柱浓度为：优选的，步骤二的函数式中，x代表反演的状态量，y为实测量，F(x)为前向模型的模拟值，K为权重函数，χ2(x)代表最小价值函数。与现有技术相比，本专利技术的有益效果是：1、利用地基MAX-DOAS获得的大气吸收光谱信息，结合实测的气体成分二氧化氮、二氧化硫、甲醛、亚硝酸及乙二醛等的吸收截面，基于差分吸收光谱解析方法反演大气污染气体(主要为二氧化氮、二氧化硫、甲醛、亚硝酸及乙二醛等)的对流层柱浓度分布结果及基于最优估计算法和SCIATRAN辐射传输模型反演的气溶胶及大气污染气体的垂直分布廓线，进行立体遥感测量，可以同时反演获取不同仪器观测方位角的垂直分布廓线。2、获取区域范围大气痕量成份数据，为我国污染气体的总量监测和特性研究提供关键的数据支持，以解决地基被动遥感观测垂直廓线反演算法的难题，拓展地基MAX-DOAS的应用范围，实现对大气气溶胶及污染气体(二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等)时空变化特征进行快速、连续、高效监测，在进行大气气溶胶分布立体监测的同时，可以获取大气气溶胶的光学性质，满足我国环境污染研究和环境监测的需要。3、建立适合大气污染气体垂直廓线反演的方法，计算便于程序化，实现无人值守的自动反演测量。附图说明图1为本专利技术的实现流程图；图2为本专利技术的MAX-DOAS反演的气溶胶廓线图；图3为本专利技术的MAX-DOAS反演的痕量气体廓线图；图4为本专利技术的MAX-DOAS反演的气溶胶光学厚度与太阳光度计实测的气溶胶光学厚度的对比图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。请参阅图1，基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，包括如下步骤：第一步：光谱波段选择，差分吸收光谱解析方法是通过测量各种气体成分的特征吸收波段来进行定性和定量反演的，不同大气成分具有不同的吸收波段，地基MAX-DOAS拥有紫外波段和可见波段两个光谱仪，能够探测到O4、二氧化氮、二氧化硫、甲醛、亚硝酸及乙二醛等痕量气体。表1为大气中主要痕量气体成分的特征光谱吸收反演波段范围。反演基于朗伯-比尔定律(如下式)，使用QDOAS对测得的太阳散射光谱进行反演，计算紫外波段(360nm)和可见波段(477.1nm)O4的差分斜柱浓度及二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等痕量气体的差分斜柱浓度。同时反演中要对仪器的暗电流影响及光谱偏移、非线性的影响进行去除和校正。然后根据拟合的残差及相对误差进行判定拟合的质量：式中，I(λ,L)和I0(λ)分别为测量和初始太阳散射光谱，σj(λ)为不同气体的吸收截面，cj(l)为气体浓度，εM(λ,l)和εR(λ,l)分别为米散射和瑞利散射的消光系数，L为光程长度，测得不同气体的差分斜柱浓度S(L)(ΔSCD)为：表1MAX-DOAS气体反演波段气体名称MAX-DOAS反演波段(nm)O4338-370，425-490NO2338-370，411-445，425-490SO2305-317.5，312-326HCHO324.5-359，336.5-359HONO335-373CHOCHO431-461步骤二：利用第一步中得到O4的差分斜柱浓度反演气溶胶的垂直廓线，反演基于最优估计算法(OEM)结合前向模型SCIATRAN。这个过程中需要预设先验廓线及先验廓线的误差。算法中规定了4种先验廓线的形状(E指数衰减型、线型、高斯型及玻尔兹曼型)，先验廓线的误差规定为不同高度先验浓度值的100％。温度、压力廓线可以使用对应维度美国标准大气或者使用WRF-Chem化学模型模拟得到的结果。气溶胶的预设光学特性参数(单次散射反照率和不对称因子)可以从国际开源网站AERONET上下载得到。反演的迭代过程如下：反演过程中实测结果与前向模型模拟结果拟合程度用最小价值函数来评估，函数式如下：式中，x代表反演的状态量，y为实测量，F(x)为前向模型的模拟值，K为权重函数，χ2(x)代表最小价值函数。图2中显示的是反演得到的气溶胶廓线的实例及与实际廓线的对比情况。步骤本文档来自技高网...

【技术保护点】
1.基于MAX‑DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，其特征在于：包括如下步骤：步骤一：基于朗伯‑比尔定律，使用QDOAS计算紫外波段和可见波段O4的差分斜柱浓度及二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等痕量气体的差分斜柱浓度，朗伯‑比尔定律：

【技术特征摘要】
1.基于MAX-DOAS对气溶胶及痕量污染气体的垂直廓线反演算法，其特征在于：包括如下步骤：步骤一：基于朗伯-比尔定律，使用QDOAS计算紫外波段和可见波段O4的差分斜柱浓度及二氧化氮、二氧化硫、甲醛、亚硝酸、乙二醛等痕量气体的差分斜柱浓度，朗伯-比尔定律：步骤二：利用不同波段的O4DSCD反演对应波段的气溶胶廓线，基于最优估计算法结合前向模型SCIATRAN反演出气溶胶的垂直廓线，迭代反演算法公式：反演过程中实测结果与前向模型模拟结果拟合程度用最小价值函数来评估，函数式：步骤三：将步骤二中反演的气溶胶廓线作为输入量结合实测的痕量气体差分斜柱浓度，反演痕量气体廓线，反演的迭代过程与步骤二中气溶胶的反演相同；步骤四：提取反演的气溶胶光学厚度与太阳光度计实测的气溶胶光学厚度做对比验证，提取痕量气体廓线最底层的浓度值做小时平均后与国控站点的实测地面气体浓度做对比验证，验证过程包括数据绝对值的对比及相关性的验证。2.根据权利要求1所述的基于MAX-DOA...

【专利技术属性】
技术研发人员：刘诚，
申请(专利权)人：刘诚，
类型：发明
国别省市：安徽,34