一种面向城市臭氧超标污染遥感溯源方法技术

本发明专利技术涉及大气污染溯源技术领域，具体地说是一种面向臭氧超标污染遥感溯源技术方法，通过融合地面观测数据、后向轨迹、卫星遥感和电子地图POI数据构建了臭氧污染溯源技术方法及系统，准确定位污染源，列出污染清单，不仅能够解析臭氧污染成因，还能精准锁定污染来源，指导臭氧污染源精准管控工作，从而弥补传统的污染溯源排查方法中过于依赖污染源清单数据、覆盖面不全以及污染确定性不足导致无法准确锁定污染企业的局限性。锁定污染企业的局限性。锁定污染企业的局限性。

【技术实现步骤摘要】
一种面向城市臭氧超标污染遥感溯源方法


[0001]本专利技术涉及大气环境质量污染溯源
，具体地说是一种面向城市臭氧超标污染遥感溯源方法。

技术介绍

[0002]臭氧超标不像颗粒物超标那样明显影响大气能见度，不易察觉，但高浓度臭氧仍可对人体健康、作物生长造成危害，已成为现阶段夏季主要大气污染物，日益受到社会公众关注。近年来，我国多个城市O3浓度超标，臭氧污染形势不容乐观，污染防治任务仍然较重，特别是在夏季O3污染高发时期，成为影响城市环境空气质量的主要污染因子。O3污染成因复杂，主要是由于人类活动排放了大量的NOx和VOCs，和天然源一起在太阳辐射、温度、湿度等环境中通过复杂的光化学反应生成O3。O3与前体物复杂的非线性关系给O3污染溯源带来了较大困难，目前VOCs尚未列入国家环境空气质量监测体系，主要通过建设超级站点对城市主要VOCs进行监测，并根据VOCs种类特征和O3生成潜势进行溯源。但是，由于超级站点建设成本高昂，有限的监测站点无法对大范围O3污染及来源进行全覆盖监测，同时由于O3生命周期较长，存在城市间O3及其前体物污染传输，因此仅依靠单个城市超级站进行O3污染溯源效果有限。
[0003]卫星遥感作为一种新兴技术，超光谱卫星遥感技术可以实现对O3前体物甲醛(HCHO)和二氧化氮(NO2)，为O3污染溯源监测提供了重要手段，它具有宏观、动态、客观、准确等特点，在信息获取方面与传统的地面监测手段相比，在空间范围上具有其独特的优势，通过超光谱和多光谱监测技术连续获取大范围区域的O3前体物浓度分布以及土地利用分布时空变化，因此可以有效地支撑区域O3污染溯源。
[0004]目前，国内外主流大气环境卫星遥感数据(如OMI、TropOMI等)的O3前体物监测能够实现最高每天1次、空间分辨率最高达到3.5公里；土地利用则最高能达到优于1m的空间分辨率，高分一号、高分二号、高分六号、资源三号等多星组网观测能达到2m分辨率卫星数据每月覆盖1次的监测频次，可以实现针对O3污染源的目标识别和无组织排放精细化解译；相关的天气再分析资料也能达到24小时内的逐小时、公里级别，可以满足O3污染溯源的需要。
[0005]目前，主要基于地面超级监测站点实时测量开展O3污染溯源分析，但是溯源覆盖面小且准确性较低，一方面超级站点的测量设备建设和运行维护成本较高，一般城市难以建设多个超级站点进行O3污染源溯源，因此只是基于有限站点数据很难实现大范围区域的全面溯源监测，无法准确定位到具体区域甚至具体企业；另一方面，为弥补地面监测站点覆盖不足的缺陷，部分研究引入超光谱卫星遥感资料或者大气化学模式模拟资料筛选出几公里甚至几十公里的O3污染高值区，为臭氧污染溯源提供目标参考区域，实际上对于环境管理人员来说，几公里的范围企业量大面广，逐一核查的工作量也相对较大，同时O3污染形成机理较为复杂，除了城市本地的排放还需要考虑城市间的O3污染区域传输影响。因此，传统的溯源方法中监测覆盖面不足或者监测数据分辨率较低等会导致无法对O3污染进行精准
溯源。
[0006]因此，为了解决上述问题，本申请提出了一种面向城市臭氧超标污染遥感溯源方法，通过融合多尺度卫星遥感、后向轨迹和潜在污染源贡献分析方法构建了O3污染遥感溯源技术方法，来综合准确反映O3污染来源，不仅能够反映区域的O3污染前体物时空分布变化和区域污染传输，还能精细化反映污染源的无组织排放和潜在污染风险，解析O3污染来源，指导O3污染源精准管控工作，从而弥补传统的预测方法中监测覆盖面不足或者监测数据分辨率较低等会导致无法对O3污染进行精准溯源的局限性。

技术实现思路

[0007]本专利技术的目的是填补现有技术的空白，提供了一种面向城市臭氧超标污染遥感溯源方法，通过融合地面观测数据、后向轨迹、卫星遥感和电子地图POI数据构建了臭氧污染溯源技术方法及系统，准确定位污染源，列出污染清单，不仅能够解析臭氧污染成因，还能精准锁定污染来源，指导臭氧污染源精准管控工作，从而弥补传统的污染溯源排查方法中过于依赖污染源清单数据、覆盖面不全以及污染确定性不足导致无法准确锁定污染企业的局限性。
[0008]为了达到上述目的，本专利技术提供一种面向城市臭氧超标污染遥感溯源方法，该方法包括步骤：
[0009]S1.根据城市环境空气质量监测站点的O3浓度小时监测数据，利用潜在源贡献因子(PSCF)方法和浓度权重轨迹分析(CWT)方法计算加权潜在源贡献因子(WPSCF)和权重浓度，提取5km潜在污染源网格单元；
[0010]S2.根据卫星遥感监测的区域甲醛(HCHO)和二氧化氮(NO2)柱浓度数据，重构合成高分辨率O3生成前体物分布结果，通过上下文方法，结合潜在污染源网格单元提取1km污染排放热点网格；
[0011]S3.根据热点网格圈定矢量区域范围，利用关键字搜索方法从电子地图中获取污染排放热点网格单位内的相关企业POI；
[0012]S4.利用高分辨率影像对的企业POI进行解译核查，筛选出污染企业清单。
[0013]S1包括：
[0014]S1.1按月收集城市O3浓度小时监测数据和全球数据同化预报系统(GDAS)发布的气象数据；
[0015]S1.2以城市建成区环境空气质量监测站点的几何中心作为起始点，采用500m高度作为大气边界层的平均流场，通过HYSPLIT模型对每天的气团移动24h后向轨迹进行模拟；
[0016]S1.3将评估城市及周边区域划分为5km
×
5km的网格，结合城市O3小时浓度以160μg/m3作为污染轨迹判断标准，即经过某网格的气团轨迹达到城市时所对应的O3浓度超过160μg/m3则判定该轨迹为污染轨迹，反之则判定为清洁轨迹；
[0017]S1.4逐网格计算潜在源加权贡献因子(WPSCF)，计算方式如下：
[0018][0019]其中，WPSCF
ij
为第i行j列网格的潜在源加权贡献因子，p
ij
为经过第i行j列网格的污染轨迹数量，n
ij
为经过第i行j列网格的所有轨迹数量。
[0020]S1.5根据潜在源加权贡献因子(WPSCF),筛选出潜在源加权贡献因子超过0.6的5km
×
5km网格单位作为潜在污染源网格。
[0021]S2包括：
[0022]S2.1所述卫星遥感资料包括通过Sentinel
‑
5P卫星TropOMI传感器探测数据反演获取的O3生成前体物甲醛(HCHO)和二氧化氮(NO2)柱浓度逐日产品数据，并从产品中提取每日5.5km
×
3.5km格点信息，包括格点经度、纬度、HCHO柱浓度、NO2柱浓度以及相应的qa值。
[0023]S2.1按月将逐日的格点数据汇总合并，并根据qa值筛选出qa大于0.5的经度、纬度、HCHO柱浓度和NO2柱浓度月度监测数据集。根据月度监测数据，通过克里金插值方法将甲醛(HCHO)和二氧化氮(NO2)月度数据集采样成城市及本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种面向城市臭氧超标污染遥感溯源方法，其特征在于，该方法包括步骤：S1.根据城市环境空气质量监测站点的O3浓度小时监测数据，利用潜在源贡献因子(PSCF)方法和浓度权重轨迹分析(CWT)方法计算加权潜在源贡献因子(WPSCF)和权重浓度，提取5km潜在污染源网格单元；S2.根据卫星遥感监测的区域甲醛(HCHO)和二氧化氮(NO2)柱浓度数据，重构合成高分辨率O3生成前体物分布结果，通过上下文方法，结合潜在污染源网格单元提取1km污染排放热点网格；S3.根据热点网格圈定矢量区域范围，利用关键字搜索方法从电子地图中获取污染排放热点网格单位内的相关企业POI；S4.利用高分辨率影像对的企业POI进行解译核查，筛选出污染企业清单。2.根据权利要求1所述的一种面向城市臭氧超标污染遥感溯源方法，其特征在于，所述S1包括：S1.1按月收集城市O3浓度小时监测数据和全球数据同化预报系统(GDAS)发布的气象数据；S1.2以城市建成区环境空气质量监测站点的几何中心作为起始点，采用500m高度作为大气边界层的平均流场，通过HYSPLIT模型对每天的气团移动24h后向轨迹进行模拟；S1.3将评估城市及周边区域划分为5km
×
5km的网格，结合城市O3小时浓度以160μg/m3作为污染轨迹判断标准，即经过某网格的气团轨迹达到城市时所对应的O3浓度超过160μg/m3则判定该轨迹为污染轨迹，反之则判定为清洁轨迹；S1.4逐网格计算潜在源加权贡献因子(WPSCF)，计算方式如下：其中，WPSCF
ij
为第i行j列网格的潜在源加权贡献因子，p
ij
为经过第i行j列网格的污染轨迹数量，n
ij
为经过第i行j列网格的所有轨迹数量；S1.5根据潜在源加权贡献因子(WPSCF)，筛选出潜在源加权贡献因子超过0.6的5km*5km网格单位作为潜在污染源网格。3.根据权利要求1所述的一种面向城市臭氧超标污染遥感溯源方法，其特征在于，所述S2包括：S2.1所述卫星遥感资料包括通过Sentinel
‑
5P卫星TropOMI传感器探测数据反演获取的O3生成前体物甲醛(HCHO)和二氧化氮(NO2)柱浓度逐日产品数据，并从产品中提取每日5.5km
×
3.5km格点信息，包括格点经度、纬度、HCHO柱浓度、NO2柱浓度以及相应的qa值；S2.1按月将逐日的格点数据汇总合并，并根据qa值筛选出qa大于0.5的经度、纬度、
HCHO柱浓度和NO2柱浓度月度监测数据集；根据月度监测数据，通过克里金插值方法将甲醛(HCHO)和二氧化氮(NO2)月度数据集采样成城市及周边地区高空间分辨率的浓度分布数据，分别形成1km*1km分辨率的甲醛(HCHO)...

【专利技术属性】
技术研发人员：陈辉，赵爱梅，张丽娟，杨艺，翁国庆，马鹏飞，周春艳，陈琳涵，张连华，代亚贞，周伟，赵少华，王中挺，
申请(专利权)人：生态环境部卫星环境应用中心，
类型：发明
国别省市：