利用临边遥感数据反演大气臭氧剖面的层析成像方法-山东科威数控机床有限公司

专利名称：利用临边遥感数据反演大气臭氧剖面的层析成像方法
技术领域：
本发明涉及一种大气化学遥感，尤其是利用氧红外波段大气临边遥感数据通过层析成像反演方法与奇氧族光化学模型结合，反演大气臭氧二维剖面的层析成像方法。
背景技术：
臭氧是大气中最重要的温室气体之一。主要分布在大气平流层内的臭氧是地球生命免遭紫外线伤害的“保护伞”；但是，人类生活的对流层中过高的臭氧含量是有害的甚至有毒的大气污染，能够损害人的眼、肺等，以及影响农作物的产量。1929年英国人Dobson开始观测研究大气臭氧，特别是自1979年卫星上的 TOMS(Total Ozone Mapping Spectrometer)禾口 SBUV(Solar BackscatteredUltraviolet) 观测大气臭氧，发现了南极臭氧洞以后，大气臭氧的观测研究进入了空前的快速发展轨道。近年来，研究发现北极大气臭氧在衰减，中国青藏高原的大气臭氧含量也在减少。在新一代对地观测系统中，美国宇航局发射的Terra-Aqua卫星中都搭载了的大气化学探测的遥感器，随后于2004年7月15日又专门发射了大气化学卫星Aura。同时，在欧空局2002年发射的最大的环境卫星Envisat中也搭载了多颗大气化学成分探测的遥感器。当然，所有这些大气化学成分探测中，仍然以大气臭氧等温室气体探测为主。相对传统地面和航空观测技术，卫星遥感大气臭氧能够提供一个全球三维动态的大气臭氧剖面和总量测量。卫星遥感大气化学成分，从卫星观测几何特征上，主要有天底、掩星和临边三种观测方式，分别利用大气的反射、透射、后向散射，以及大气的热辐射来反演大气的化学成分。不同的观测几何特征，决定了卫星遥感大气臭氧的覆盖范围和大气化学成分数据的水平和垂直分辨率。天底观测能够确保卫星的全球观测范围，但只是提供平流层臭氧总量信息和有限垂直分辨率的臭氧剖面，如比较早的TOMS和GOME (Global Ozone Monitoring Experiment)等仪器。掩星观测类似临边观测，充分利用垂直网格内丰富的光谱信息和相对高的探测深度，能够提供大气成分的剖面信息，但其观测范围只限制在50-80纬度范围内，而且还取决于季节的影响，如SAGEIII、Ρ0ΑΜ、G0M0S等仪器。相对而言，临边观测只记录大气散射的太阳辐射，能够提供一个更好的地面覆盖和更高垂直分辨率的大气剖面测量，因此，自2001年以来的多颗大气探测遥感卫星都采用了这种观测技术，如 Odin JlMJl Combined Optical Spectrograph and Infrared Imager System(OSIRIS), Envisat 卫星上的 Scanning Imaging Absorption SpectroMeter forAtmospheric CHartographY(SCIAMACHY), Hitchhiker Junior 上的 the ShuttleOzone Limb Sounding Experiment (SOLSE) ,Aura 卫星上的 High ResolutionDynamics Limb Sounder (HIRDLS)禾口 Microwave Limb Sounder (MLS)。因此，论文将以卫星遥感大气Limb散射反演臭氧剖面为例，研究层析成像技术在臭氧剖面反演中的应用。针对不同的卫星观测结果，近年来开发了许多不同算法来反演大气臭氧剖面。 OSIRIS臭氧剖面反演采用了与SOLSE反演一样的技术，基于正规化的Limb辐射剖面，综合利用辐射传输模型(LIMBTRAN)和最优化估计(OE)方法；SCIAMACHY臭氧剖面利用辐射传输模型(SCIARAYS)和臭氧在可见光波长范围内的Chappuis波段进行反演的；GOME臭氧剖面反演算法也是基于辐射传输模型模拟的辐射和GOME测量的辐射差的最小化方法。这些反演大气臭氧剖面的物理算法都是利用了一些已知的大气臭氧剖面测量数据作为先验知识或者它们的正规化、最优化方法来反演的，同时假设大气化学成分具有均勻的水平结构。尽管这些方法都得到了相对较高的剖面反演精度，即具有较高的垂直分辨率，但是大气臭氧的水平结构信息反演却无能为力，即其水平分辨率却非常低。且基于地基、火箭和卫星平台，利用臭氧对紫外辐射在Hartley波段的吸收特性反演臭氧浓度的应用已有很多。这些直接的观测能得到40-75km高度间的精确的臭氧数据，但通常由于吸收不足反演高度仅限于75km以下。现有的反演方法只能得到臭氧总量信息或垂直分辨率很低的廓线。

发明内容
本发明的目的就在于针对上述现有技术的不足，提供一种利用临边遥感数据反演大气臭氧剖面的层析成像方法。要解决的技术问题是要解决视线与大气网格截距的计算问题；要解决如何由红外辐射原理将气辉体发射率与辐射亮度联系起来组成正向模拟方程，并选择数值算法来对大型稀疏超正定方程进行解算，从而反演大气体发射率；要解决如何由大气网格体发射率计算相应的臭氧数密度，从而实现臭氧剖面的反演。所采用的技术方案是利用临边遥感数据反演大气臭氧剖面的层析成像方法是通过分别计算视线向量与大气网格径向垂直分界和大气层水平分界交点，合并两种类型交点再计算相邻交点距离得到截距；将层析成像方法应用到联系体发射率与辐射亮度方程上，并选用乘迭代重建技术，得到了体发射率的反演迭代模型，对体发射率进行解算；最后利用奇氧族光化学模型将体发射率与臭氧数密度联系起来，选择牛顿迭代法实现臭氧剖面的反演。利用临边遥感数据反演大气臭氧层剖面的层析成像方法，包括如下顺序和步骤a、临边遥感轨道数据的读取与预处理先读取卫星轨道辐射亮度数据、卫星地理坐标及观测时间，并读取辐射亮度对应观测视线的切点的地理坐标、高度、太阳天顶角、相对方位角数据；b、根据读取的卫星地理坐标和观测视线切点的地理坐标，计算卫星和切点的地心三维直角坐标及轨道面内的升交点地理坐标；C、依据卫星升交点地理坐标确定反演剖面沿轨道方向大气角度边界，由切点升交点地理坐标确定反演大气高度范围；d、然后由卫星和切点的轨道面升交点坐标系的地理坐标计算卫星和每个切点的直角坐标；e、大气气辉体发射率反演首先根据反演大气范围内的卫星和切点升交点直角坐标计算切点到卫星的视线向量，然后计算视线向量与反演范围内大气网格截距；f、由氧红外波段大气临边观测辐射亮度实质是体发射率沿视线积分原理出发，建立层析成像反演模型，选择乘代数重建方法进行解算；其中包括初值估算，和迭代计算两个步骤；
权利要求
一种利用临边遥感数据反演大气臭氧剖面的层析成像方法，其特征在于，包括如下顺序和步骤a、临边遥感轨道数据的读取与预处理先读取卫星轨道辐射亮度数据、卫星地理坐标及观测时间，并读取辐射亮度对应观测视线的切点的地理坐标、切高、太阳天顶角、相对方位角数据；b、根据读取的卫星地理坐标和观测视线切点的地理坐标，计算卫星和切点的地心三维直角坐标及轨道面内的升交点地理坐标；c、依据卫星升交点地理坐标确定反演大气剖面沿轨道方向角度边界，由切点升交点地理坐标确定反演大气高度范围；d、然后由卫星和切点的轨道面升交点坐标系的地理坐标计算卫星和每个切点的直角坐标；e、大气气辉体发射率反演首先根据反演大气范围内的卫星和切点升交点直角坐标计算切点到卫星的视线向量，然后计算视线向量与反演范围内大气网格截距；f、由氧红外波段大气临边观测辐射亮度实质是体发射率沿视线积分原理出发，建立层析成像反演模型，选择乘代数重建方法进行解算；其中包括初值估算，和迭代计算两个步骤； <mrow><msubsup> <mi>V</mi> <mi>j</mi> <mrow><mo>(</mo><mi>n</mi><mo>)</mo> </mrow></msubsup><mo>=</mo><msubsup> <mi>V</mi> <mi>j</mi> <mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow></msubsup><munder> <mi>Σ</mi> <mi>i</mi></munder><mrow> <mo>(</mo> <mfrac><msub> <mi>O</mi> <mi>i</mi></msub><msubsup> <mi>O</mi> <mi>iest</mi> <mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow></msubsup> </mfrac> <msub><mi>β</mi><mi>ij</mi> </msub> <mo>)</mo></mrow> </mrow> <mrow><msubsup> <mi>V</mi> <mi>j</mi> <mrow><mo>(</mo><mn>1</mn><mo>)</mo> </mrow></msubsup><mo>=</mo><munder> <mi>Σ</mi> <mi>i</mi></munder><mrow> <mo>(</mo> <mfrac><msub> <mi>O</mi> <mi>i</mi></msub><mrow> <msub><mi>Σ</mi><mi>j</mi> </msub> <msub><mi>L</mi><mi>ij</mi> </msub></mrow> </mfrac> <msub><mi>β</mi><mi>ij</mi> </msub> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>Σ</mi> <mi>j</mi></msub><msub> <mi>L</mi> <mi>ij</mi></msub><msubsup> <mi>V</mi> <mi>j</mi> <mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow></msubsup><mo>=</mo><msubsup> <mi>O</mi> <mi>iest</mi> <mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo> </mrow></msubsup> </mrow> <mrow><msub> <mi>β</mi> <mi>ij</mi></msub><mo>=</mo><mfrac> <msub><mi>L</mi><mi>ij</mi> </msub> <mrow><msub> <mi>Σ</mi> <mi>i</mi></msub><msub> <mi>L</mi> <mi>ij</mi></msub> </mrow></mfrac> </mrow>其中，和分别为第j个大气网格中体发射率初值和经过第n次迭代后的反演值，Oi是第i条视线辐射亮度，O(n 1)iest是第i条视线模拟估计辐射亮度，Lij是第i条视线通过第j个大气网格截距，。βij称为权因子。g、大气臭氧剖面反演首先由比尔定律，用初始大气数据库计算大气网格在不同太阳天顶角下的太阳辐射，由此计算奇氧族光化学反应的光解反应系数；h、由奇氧族光化学反应出发，建立气辉体发射率连续方程； <mrow><mo>[</mo><mi>O</mi><mrow> <mo>(</mo> <mmultiscripts><mi>D</mi><mprescripts/><none/><mn>1</mn> </mmultiscripts> <mo>)</mo></mrow><mo>]</mo><mo>=</mo><mfrac> <mrow><msub> <mi>J</mi> <mi>H</mi></msub><mo>[</mo><msub> <mi>O</mi> <mn>3</mn></msub><mo>]</mo><mo>+</mo><mrow> <mo>(</mo> <msub><mi>J</mi><mi>SRC</mi> </msub> <mo>+</mo> <msub><mi>J</mi><mrow> <mi>Ly</mi> <mo>-</mo> <mi>α</mi></mrow> </msub> <mo>)</mo></mrow><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mo>]</mo> </mrow> <mrow><msub> <mi>A</mi> <mi>D</mi></msub><mo>+</mo><msub> <mi>k</mi> <mi>D</mi></msub><mo>[</mo><msub> <mi>Q</mi> <mi>D</mi></msub><mo>]</mo> </mrow></mfrac> </mrow> <mrow><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mrow> <mo>(</mo> <msup><mi>b</mi><mn>1</mn> </msup> <mi>Σ</mi> <mo>)</mo></mrow><mo>]</mo><mo>=</mo><mfrac> <mrow><msub> <mi>g</mi> <mi>Σ</mi></msub><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mo>]</mo><mo>+</mo><mi>φ</mi><msub> <mi>k</mi> <mrow><mi>D</mi><mo>,</mo><msub> <mi>O</mi> <mn>2</mn></msub> </mrow></msub><mo>[</mo><mi>O</mi><mrow> <mo>(</mo> <mmultiscripts><mi>D</mi><mprescripts/><none/><mn>1</mn> </mmultiscripts> <mo>)</mo></mrow><mo>]</mo><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mo>]</mo> </mrow> <mrow><msub> <mi>A</mi> <mi>Σ</mi></msub><mo>+</mo><msub> <mi>k</mi> <mi>Σ</mi></msub><mo>[</mo><msub> <mi>Q</mi> <mi>Σ</mi></msub><mo>]</mo> </mrow></mfrac> </mrow> <mrow><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mrow> <mo>(</mo> <msup><mi>a</mi><mn>1</mn> </msup> <mi>Δ</mi> <mo>)</mo></mrow><mo>]</mo><mo>=</mo><mfrac> <mrow><msub> <mi>J</mi> <mi>H</mi></msub><mo>[</mo><msub> <mi>O</mi> <mn>3</mn></msub><mo>]</mo><mo>+</mo><msub> <mi>g</mi> <mi>Δ</mi></msub><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mo>]</mo><mo>+</mo><msub> <mi>k</mi> <mi>Σ</mi></msub><mo>[</mo><msub> <mi>Q</mi> <mi>Σ</mi></msub><mo>]</mo><mo>[</mo><msub> <mi>O</mi> <mn>2</mn></msub><mrow> <mo>(</mo> <msup><mi>b</mi><mn>1</mn> </msup> <mi>Σ</mi> <mo>)</mo></mrow><mo>]</mo> </mrow> <mrow><msub> <mi>A</mi> <mi>Δ</mi></msub><mo>+</mo><msub> <mi>k</mi> <mi>Δ</mi></msub><mo>[</mo><msub> <mi>Q</mi> <mi>Δ</mi></msub><mo>]</mo> </mrow></mfrac> </mrow>其中，J，g，k和A分别为光解离系数、光激发系数、碰撞湮灭系数和自发辐射系数。i、结合光解反应系数、碰撞湮灭和自发系数，输入初始大气廓线，对气辉体发射率进行模拟；j、基于层析成像反演的二维体发射率剖面，利用牛顿迭代法进行臭氧二维剖面的反演；k、用迭代计算后的臭氧再进行模拟，然后利用上面的迭代计算式再迭代计算臭氧廓线。FSA00000269377700015.tif,FSA00000269377700016.tif
全文摘要
本发明涉及一种利用临边遥感数据反演大气臭氧剖面的层析成像方法。利用临边遥感数据反演大气臭氧剖面的层析成像方法是通过分别计算视线向量与大气网格径向垂直分界和大气层水平分界交点，合并两种类型交点再计算相邻交点距离得到截距；将层析成像方法应用到联系体发射率与辐射亮度方程上，并选用乘迭代重建技术，得到了体发射率的反演迭代模型，对体发射率进行解算；最后利用奇氧族光化学模型将体发射率与臭氧数密度联系起来，选择牛顿迭代法实现臭氧剖面的反演。
文档编号G01N21/17GK101936881SQ20101028156
公开日2011年1月5日申请日期2010年9月15日优先权日2010年9月15日
发明者孟治国, 杨莹, 汪自军, 路鹏, 陈圣波申请人:吉林大学