一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法技术

本发明专利技术公开了一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，包括：通过卫星观测以及地面观测获得总悬浮体积谱分布、波长为500nm的细模态颗粒物与总悬浮颗粒物光学厚度的比值（FMF）、粗模态颗粒物的体积浓度；通过上述观测数据建立粗模态颗粒物体积消光比模型，在该过程中，以颗粒物半径1.0μm作为细模态颗粒物和粗模态颗粒物的分界，对粗模态颗粒物的体积浓度谱分布进行积分，得到粗模态颗粒物的体积浓度；根据观测值计算得到粗模态颗粒物体积消光比；根据观测值统计得到粗模态颗粒物体积消光比关于FMF的统计关系；将得到的粗模态颗粒物体积消光比与细模态颗粒物体积消光比线性加权得到总悬浮颗粒物体积消光比；计算总悬浮颗粒物的质量浓度。

A remote sensing physical estimation method for the mass concentration of total suspended particles near the ground

【技术实现步骤摘要】
一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法
本专利技术涉及卫星遥感
，尤其是一种近地面的大气总悬浮颗粒物质量浓度遥感物理估算方法。
技术介绍
大气悬浮颗粒物也称为大气总悬浮，包含悬浮在大气中的各种固体和液体颗粒物，粒径范围在0.001-100μm之间。随着这经济增长，人口膨胀和城市化导致空气污染物大量排放，大气颗粒物污染已经成为最严重的大气污染物之一，可以带来严重的健康风险，气候变化，建筑物腐蚀和能见度降低。目前主要的监测手段包括地面采样和卫星遥感。但是，布设成本高，监测范围小，难以达到空间均匀分布，无法提供大范围空间连续覆盖的监测结果。因此，对全球人口的暴露以及由远距离地区的运输引起的高污染事件的研究很难通过现有的地面监测网络来实现。然而，卫星遥感技术具有大面积同步观测的特点，极大的弥补了地基监测的不足。卫星遥感观测在于自然存在状态下对颗粒物进行观测，不破坏颗粒物的原有状态，避免了采样测量时加热、干燥等过程对颗粒物原有性状的改变。因此，发展卫星遥感近地面颗粒物遥感估算技术对大气颗粒物污染的治理和监测具有重要意义。目前，用于卫星遥感估算颗粒物质量浓度的方法主要包括数学统计法，化学传输模式法和物理模型方法。数学统计方法得到的结果往往具有较高的精度，然而，其物理意义不明确，并且依赖于地面观测数据。化学传输模式方法，物理机制复杂，无法在小时级别上获得理想的结果，并且消耗大量的运算资源。物理模型方法具有简明的物理机制，不依赖于地面在位测量和不消耗计算资源的优势，被广泛的以用于近地面颗粒物质量浓度的估算。而现有技术中的方法只能得到细粒子总悬浮质量浓度，无法提供总悬浮颗粒物质量浓度，这主要是由于目前不存在覆盖所有气溶胶颗粒物粒径范围的测量仪器。对自然状态下卫星观测得到整层大气的颗粒物消光进行细粒子、高度和湿度订正，利用细粒子光学厚度比计算细颗粒物总悬浮光学厚度。假设颗粒物垂直分布满足负指数分布，采用总悬浮层高计算近地面细粒子总悬浮光学厚度。采用吸湿增长因子获得近地面干燥细颗粒物的光学厚度。基于地基遥感观测网提供的历史资料，分析不同类型细模态总悬浮颗粒物的消光与体积关系的变化规律，寻找对细颗粒物的消光与体积转换特性进行描述的遥感参数。再根据密度体积乘积，得到近地面干燥的细粒子总悬浮质量浓度。
技术实现思路
针对现有技术存在的问题，本专利技术的目的在于提供一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，通过建立粗模态颗粒物体积消光比模型来测量近地面大气总悬浮颗粒物质量浓度。为实现上述目的，本专利技术的技术方案如下：一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，包括：通过卫星观测以及地面观测获得总悬浮体积谱分布、波长为500nm的细模态颗粒物与总悬浮颗粒物光学厚度的比值(FMF)、粗模态颗粒物的体积浓度；通过上述观测数据建立粗模态颗粒物体积消光比模型；在建立粗模态颗粒物体积消光比模型的过程中，以颗粒物半径1.0μm作为细模态颗粒物和粗模态颗粒物的分界，对粗模态颗粒物的体积浓度谱分布进行积分，得到粗模态颗粒物的体积浓度；根据观测值计算得到粗模态颗粒物体积消光比；根据观测值统计得到粗模态颗粒物体积消光比关于FMF的统计关系；将得到的粗模态颗粒物体积消光比与细模态颗粒物体积消光比线性加权得到总悬浮颗粒物体积消光比；计算总悬浮颗粒物的质量浓度。进一步，所述总悬浮颗粒物体积消光比通过以下公式进行计算：其中，VE表示总悬浮颗粒物体积消光比，VEf和VEc分别是细模态和粗模态颗粒物体积消光比，Vcol表示总悬浮颗粒物柱体积浓度，AOD表示总悬浮颗粒物光学厚度，FMF表示细模态颗粒物与总悬浮颗粒物光学厚度的比值，Vf,col和Vc,col表示细模态和粗模态颗粒物的体积浓度，VEf是关于FMF的二次多项式。进一步，通过所述粗模态颗粒物体积消光比与1-FMF的倒数进行线性统计回归得到所述粗模态颗粒物体积消光比关于FMF的统计关系：其中，FMF小于1.0。进一步，在计算总悬浮颗粒物质量浓度过程中，首先通过总悬浮颗粒物体积消光比得到大气中总悬浮颗粒物体积浓度，进而获得近地面总悬浮颗粒物体积浓度，接着利用吸湿增长因子对近地面总悬浮颗粒物体积浓度进行干燥，再通过乘以密度得到近地面总悬浮颗粒物质量浓度：其中，PM表示总悬浮颗粒物的质量浓度，VE表示总悬浮颗粒物体积消光比，ρ表示干燥总悬浮颗粒物密度，这里采用常数1.8g/cm3，PBLH表示行星边界层高度，f(RH)表示吸湿增长因子。进一步，采用北京站点2015年1-2月观测PM10与计算PM进行比对，本专利技术计算PM与观测值变化趋势的一致性相关系数R是0.78。进一步，所述观测数据来自地基AerosolRoboticNetwork(AERONET)7个观测站点10年的观测数据，包括城市/工业型总悬浮主导的北京和GSFC站点，生物质燃烧型总悬浮主导的CuiabaMiranda和Mongu站点，沙尘总悬浮主导的SolarVillage站点，海洋总悬浮主导的AscensionIsland和Lanai站点。进一步，所述北京站点和GSFC站点分别代表光吸收型和非吸收型城市/工业总悬浮颗粒物，所述AscensionIsland和Lanai站点分别代表污染型和清洁型海洋总悬浮颗粒物。本专利技术中一种近地面的大气总悬浮颗粒物质量浓度遥感物理估算方法，通过建立粗模态颗粒物体积消光比模型，通过卫星和/或地面测量的方式，将长时间尺度的数据通过统计分析的方法建立一种通过细模态颗粒物和粗模态颗粒物的体积消光比计算近地面大气总悬浮颗粒物质量浓度的方法，解决了现有技术中卫星遥感观测只能提供细粒子气溶胶质量浓度，无法提供总悬浮颗粒物质量浓度的问题。同时，现有估算技术往往基于传统的数学统计方法，只考虑大气总悬浮的光学厚度与近地面颗粒物质量浓度的线性回归关系或者非线性回归关系，往往都忽略了两个变量之间的物理关系。其统计关系往往与季节和地理位置有关，也受到极端天气的影响，导致现有技术无法完全移植到卫星使用，需要依赖卫星观测和地面测量的数学统计关系。然而，本专利技术不仅可以将适用于全球不同地区、不同季节的、不同气溶胶类型的体积消光比的模型参数固定下来，还具有完全移植到卫星上的能力。同时，还考虑了大气总悬浮颗粒的光学厚度与近地面颗粒物质量浓度的物理关系(颗粒物的体积消光比，垂直分布，吸湿增长等)，可以实现总悬浮颗粒物的体积浓度与消光系数的转换，对抽象的气溶胶进行体积可视化和称重计算。附图说明图1为本专利技术流程图；图2为本专利技术中粗模态颗粒物体积数谱图；图3为本专利技术中粗模态颗粒物体积散点图。具体实施方式为了使本领域的技术人员更好地理解本专利技术的方案，下面结合本专利技术示例中的附图对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的示例仅仅是本专利技术的一部分示例，而不是全部的示例。基于本专利技术的中示例，本领域的普本文档来自技高网...

【技术保护点】
1.一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，其特征在于，包括：/n通过卫星观测以及地面观测获得总悬浮体积谱分布、波长为500nm的细模态颗粒物与总悬浮颗粒物光学厚度的比值FMF、粗模态颗粒物的体积浓度；/n通过上述观测数据建立粗模态颗粒物体积消光比模型；/n在建立粗模态颗粒物体积消光比模型的过程中，以颗粒物半径1.0μm作为细模态颗粒物和粗模态颗粒物的分界，对粗模态颗粒物的体积浓度谱分布进行积分，得到粗模态颗粒物的体积浓度；/n根据观测值计算得到粗模态颗粒物体积消光比；/n根据观测值统计得到粗模态颗粒物体积消光比关于FMF的统计关系；/n将得到的粗模态颗粒物体积消光比与细模态颗粒物体积消光比线性加权得到总悬浮颗粒物体积消光比；/n计算总悬浮颗粒物的质量浓度。/n

【技术特征摘要】
1.一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，其特征在于，包括：
通过卫星观测以及地面观测获得总悬浮体积谱分布、波长为500nm的细模态颗粒物与总悬浮颗粒物光学厚度的比值FMF、粗模态颗粒物的体积浓度；
通过上述观测数据建立粗模态颗粒物体积消光比模型；
在建立粗模态颗粒物体积消光比模型的过程中，以颗粒物半径1.0μm作为细模态颗粒物和粗模态颗粒物的分界，对粗模态颗粒物的体积浓度谱分布进行积分，得到粗模态颗粒物的体积浓度；
根据观测值计算得到粗模态颗粒物体积消光比；
根据观测值统计得到粗模态颗粒物体积消光比关于FMF的统计关系；
将得到的粗模态颗粒物体积消光比与细模态颗粒物体积消光比线性加权得到总悬浮颗粒物体积消光比；
计算总悬浮颗粒物的质量浓度。


2.根据权利要求1所述的近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，其特征在于，所述总悬浮颗粒物体积消光比通过以下公式进行计算：



其中，VE表示总悬浮颗粒物体积消光比，VEf和VEc分别是细模态和粗模态颗粒物体积消光比，Vcol表示总悬浮颗粒物柱体积浓度，AOD表示总悬浮颗粒物光学厚度，FMF表示细模态颗粒物与总悬浮颗粒物光学厚度的比值，Vf,col和Vc,col表示细模态和粗模态颗粒物的柱体积浓度，VEf是关于FMF的二次多项式。


3.根据权利要求2所述的近地面大气总悬浮颗粒物质量浓度遥感物理估算方法，其特征在于，通过所述粗模态颗粒物体积消光比与1-FMF的倒数进行线性统计回归得到所述粗模态颗粒物体积消光比关于FMF的统计关系：



其中，FMF小于1.0。


...

【专利技术属性】
技术研发人员：李正强，魏瑗瑗，张莹，陈杰，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京;11