一种基于气象观测数据的可吸入颗粒物浓度估算方法技术

本发明专利技术公开了一种基于大气能见度与相对湿度观测数据估算近地面可吸入颗粒物(PM10)浓度的方法，包括：S1：基于水平能见度观测数据计算近地面气溶胶散射系数；S2：基于近地面气溶胶消光系数与PM10质量浓度计算气溶胶平均质量散射效率；S3：根据季节、背景环境及气象信息，建立气溶胶平均质量散射效率与相对湿度的相关模型；S4：利用S3所述相关模型，基于大气能见度与相对湿度观测数据估算PM10浓度。本发明专利技术通过基于大气能见度、相对湿度等气象数据估算近地面颗粒物污染水平，能有效补充颗粒物地面监测站点的不足，拓展颗粒物浓度监测的空间范围。

【技术实现步骤摘要】
一种基于气象观测数据的可吸入颗粒物浓度估算方法
本专利技术涉及大气环境遥感监测
，特别涉及一种基于地面气象观测数据的近地面可吸入颗粒物(PM10)浓度估算方法。
技术介绍
作为一种主要的大气污染物，气溶胶对公共健康的影响已经是公认的事实，直接威胁着人类生存与可持续发展。其中，可吸入颗粒物PM10(空气动力学直径小于10μm的颗粒物)可达人体的支气管区，粒径小于5μm的颗粒物可达到肺泡区，部分更小的粒子甚至可以通过毛细血管进入人体血液循环系统，对心脏及心血管造成较大危害。据《2006年中国环境状况公报》统计，在我国监测的557个城市中，43.4％的城市大气质量没有达标，颗粒物为主要污染物。准确获取PM10的时空分布、来源及传输路径是衡量其污染影响，制定颗粒物防治政策的重要保障。当前PM10监测主要依靠地面环境站点获取，然而由于观测设备一般昂贵且维护复杂，环境监测站点分布不均且数量有限。而地面气象站点相对于环境监测站点而言，数量更大且空间分布较为均匀，如果基于经验方法由气象站点观测的数据估算周边PM10浓度，将有效补充环境监测站点的数据不足。目前，大气能见度和相对湿度已成为绝大部分气象站点的基本观测数据，数据质量较为可靠。更重要的是，气象数据往往比环境监测数据获取的时间更长，可用于某地区颗粒物污染状况的历史数据分析。
技术实现思路
本专利技术要解决的技术问题是基于地面气象站点观测的大气能见度及相对湿度数据，估算近地面可吸入颗粒物PM10浓度，从而进一步拓展颗粒物监测的空间覆盖能力。为解决上述技术问题，本专利技术提供了一种基于地面气象观测数据的近地面可吸入颗粒物(PM10)浓度估算方法，包括以下步骤：S1、选择建模站点：选择处于同类环境中的、用于建模的地面气象站点(具有能见度及相对湿度观测数据)和地面环境监测站点(具有PM10质量浓度观测数据)；S2、计算散射系数：基于所述地面气象站点观测的大气能见度，依据已知经验关系，计算所述地面气象站点处近地面气溶胶散射系数(aerosolscatteringcoefficient)；S3、计算平均质量散射效率：基于步骤S2所述近地面气溶胶散射系数与步骤S1中所述地面环境监测站点观测的PM10质量浓度数据，计算近地面气溶胶平均质量散射效率(aerosolaveragemassscatteringefficiency)；S4、拟合相关模型：根据目标时间段不同气溶胶类型，基于步骤S3所述气溶胶平均质量散射效率与S1中所述地面气象站点获取的相对湿度，按时间进行匹配，选取非线性经验函数拟合得到二者的相关模型；根据Mie散射理论，假定气溶胶的化学组分与粒径分布在一定时间内变化不大，则气溶胶平均质量散射效率可近似看作相对湿度的函数；S5、估算PM10浓度：利用步骤S4所述相关模型，基于与S1所述用于建模的地面气象站点处于同类环境的另一地面气象站点所观测的大气能见度与相对湿度数据，按不同气溶胶类型估算所述另一地面气象站点近地面PM10质量浓度。估算的过程是：S51、基于所述另一地面气象站点观测的大气能见度计算所述另一地面气象站点处的近地面气溶胶散射系数；S52、根据气溶胶来源类型选取目标时间段所述相关模型，由所述另一地面气象站点处观测的相对湿度计算出所述另一地面气象站点处近地面气溶胶平均质量散射效率；S53、由所述另一地面气象站点处近地面气溶胶散射系数与所述另一地面气象站点处气溶胶平均质量散射效率，计算所述另一地面气象站点处PM10质量浓度。优选地，步骤S4所述不同气溶胶类型，具体包括：根据地面站点所在区域的气象条件与气溶胶来源，将当地目标时间段气溶胶分为无机成分主导型和有机成分主导型两种类型。步骤S1和步骤S5所述同类环境，指城市或农村中设定的距离范围、或污染程度处于相同的量化范围。优选地，在所述同类环境中，所述用于建模的地面气象站点和地面环境监测站点距离A小于所述用于建模的地面气象站点和所述另一地面气象站点之间的距离B。建模时PM10观测数据与能见度、相对湿度的观测最好是在同一地点，或者二者相距越近越好；而在应用这一模型估算另一地面气象站点的PM10时，该气象站点只要与前述建模的站点环境背景相同且距离不要太远即可本专利技术通过基于大气能见度、相对湿度等气象数据估算近地面颗粒物污染水平，能有效补充颗粒物地面监测站点的不足，拓展颗粒物浓度监测的空间范围。附图说明图1是按照本专利技术一种基于地面气象观测数据的近地面可吸入颗粒物(PM10)浓度估算方法的流程图。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。图1为按照本专利技术一种实施方式的基于气象观测数据的可吸入颗粒物(PM10)估算方法的流程图，包括以下步骤：步骤S1、选择建模站点：选择用于建模的地面气象站点(具有能见度及相对湿度观测数据)和环境监测站点(具有PM10观测数据)，要求二者处于同类环境中。例如，满足：1)直线距离较近(根据环境背景特点设定距离范围A，如城市或高污染背景地区不超过5km，城郊或低污染背景地区不超过20km)，2)二者背景环境相同(周边环境或背景污染程度处于相同的量化范围)，并具有一定代表性；步骤S2、计算散射系数：基于地面气象站点观测的大气能见度，依据世界气象组织(WMO)提出的经验关系，计算近地面的气溶胶散射系数(aerosolscatteringcoefficient)；步骤S2具体包括：S21、采用对比度阈值为0.05，由观测的大气能见度计算大气总消光系数，如下式：σtotal(λ)＝3.912/VIS公式1其中，σtotal(λ)为大气总消光系数，λ为消光系数对应的波长，本说明中都取值为550nm，因此后文中λ不再标出；VIS为大气能见度观测值(单位km)。S22、基于经验公式和波长计算大气分子消光系数，如下式：其中，n为大气折射率，在海平面处取值为n-1＝293*10-6；N为分子数密度，在海平面出取值N＝266*1019cm-3；S23、从大气总消光系数中分别减去大气分子消光系数贡献，获得气溶胶消光系数σa，如下式：σa＝σtotal-σm公式3S24、基于经验的气溶胶单次散射反照率(singlescatteringalbedo)由气溶胶消光系数计算气溶胶散射系数，σa,scat＝σa·Ω公式4其中σa,scat为近地面气溶胶散射系数，Ω为气溶胶单次散射反照率，北京地区一般设置为0.85。步骤S3、计算平均质量消光效率：基于步骤S2获取的近地面气溶胶散射系数与S1中地面环境监测站点观测的PM10质量浓度数据，计算定义为二者比值的近地面气溶胶平均质量散射效率(aerosolaveragemassscatteringefficiency)；根据Mie散射理论，假定气溶胶的化学组分与粒径分布在一定时间内变化不大，则气溶胶平均质量散射效率可近似看作相对湿度的函数；步骤S3中，基于S2获取的气溶胶散射系数与S1所述地面环境监测站点的PM10质量浓度，定义气溶胶平均质量散射效率Ea,scat：Ea,scat＝σa,scat/CPM10公式5其中CPM10为PM10质量浓度。根据Mie散射理论，假定气溶胶的化学组分与粒径分布在一定时间内变本文档来自技高网...

【技术保护点】
一种基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，包括以下步骤： 选择处于同类环境中的、用于建模的地面气象站点和地面环境监测站点，所述用于建模的地面气象站点具有能见度及相对湿度观测数据，所述地面环境监测站点具有PM10质量浓度观测数据； 基于所述地面气象站点观测的大气能见度，依据已知经验关系，计算所述地面气象站点处近地面气溶胶散射系数； 基于所述近地面气溶胶散射系数与所述地面环境监测站点观测的PM10质量浓度数据，计算近地面气溶胶平均质量散射效率； 根据目标时间段不同气溶胶类型，基于所述气溶胶平均质量散射效率与所述地面气象站点获取的相对湿度，按时间进行匹配，选取非线性经验函数拟合得到二者的相关模型； 利用所述相关模型，基于与所述用于建模的地面气象站点处于同类环境的另一地面气象站点所观测的大气能见度与相对湿度数据，按不同气溶胶类型估算所述另一地面气象站点近地面PM10质量浓度。

【技术特征摘要】
1.一种基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，包括以下步骤：选择处于同类环境中的、用于建模的地面气象站点和地面环境监测站点，所述用于建模的地面气象站点具有能见度及相对湿度观测数据，所述地面环境监测站点具有PM10质量浓度观测数据；基于所述地面气象站点观测的大气能见度，依据已知经验关系，计算所述地面气象站点处近地面气溶胶散射系数；基于所述近地面气溶胶散射系数与所述地面环境监测站点观测的PM10质量浓度数据，计算近地面气溶胶平均质量散射效率；根据目标时间段不同气溶胶类型，基于所述气溶胶平均质量散射效率与所述地面气象站点获取的相对湿度，按时间进行匹配，选取非线性经验函数拟合得到二者的相关模型；利用所述相关模型，基于与所述用于建模的地面气象站点处于同类环境的另一地面气象站点所观测的大气能见度与相对湿度数据，按不同气溶胶类型估算所述另一地面气象站点近地面PM10质量浓度。2.如权利要求1所述基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，所述同类环境，指城市或农村中设定的距离范围、或污染程度处于相同的量化范围。3.如权利要求1所述基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，在所述同类环境中，所述用于建模的地面气象站点和地面环境检测站点距离A小于所述用于建模的地面气象站点和所述另一地面气象站点之间的距离B。4.如权利要求1所述基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，所述不同气溶胶类型包括无机成分主导型和有机成分主导型两种类型。5.如权利要求1所述基于气象观测数据的可吸入颗粒物浓度估算方法，其特征在于，所述估算另一地面气象站点近地面PM10质量浓度的过程是：S51、基于所述另一地面气象站点观测的大气能见度计算所述另一地面气象站点处的近地面气溶胶散射系数；S52、根据气溶胶来源类型选取目标时间段所述相关模型，由所述另一地面气象站点处观测的相对湿度计算出所述另一地面气象站点处近地面气溶胶平均质量散射效率；S53、由所述另一地面气象站点处近地面气溶胶散射系数与所述另一地面气象站点处气溶胶平均质量散射效率，计算所述另一地面...

【专利技术属性】
技术研发人员：王子峰，陈良富，陶金花，陶明辉，苏林，王中挺，邹明敏，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京;11