一种面向城市臭氧超标污染遥感溯源方法与流程

1.本发明涉及大气环境质量污染溯源技术领域，具体地说是一种面向城市臭氧超标污染遥感溯源方法。


背景技术：

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


技术实现要素：

7.本发明的目的是填补现有技术的空白，提供了一种面向城市臭氧超标污染遥感溯源方法，通过融合地面观测数据、后向轨迹、卫星遥感和电子地图poi数据构建了臭氧污染溯源技术方法及系统，准确定位污染源，列出污染清单，不仅能够解析臭氧污染成因，还能精准锁定污染来源，指导臭氧污染源精准管控工作，从而弥补传统的污染溯源排查方法中过于依赖污染源清单数据、覆盖面不全以及污染确定性不足导致无法准确锁定污染企业的局限性。
8.为了达到上述目的，本发明提供一种面向城市臭氧超标污染遥感溯源方法，该方法包括步骤：
9.s1.根据城市环境空气质量监测站点的o3浓度小时监测数据，利用潜在源贡献因子(pscf)方法和浓度权重轨迹分析(cwt)方法计算加权潜在源贡献因子(wpscf)和权重浓度，提取5km潜在污染源网格单元；
10.s2.根据卫星遥感监测的区域甲醛(hcho)和二氧化氮(no2)柱浓度数据，重构合成高分辨率o3生成前体物分布结果，通过上下文方法，结合潜在污染源网格单元提取1km污染排放热点网格；
11.s3.根据热点网格圈定矢量区域范围，利用关键字搜索方法从电子地图中获取污染排放热点网格单位内的相关企业poi；
12.s4.利用高分辨率影像对的企业poi进行解译核查，筛选出污染企业清单。
13.s1包括：
14.s1.1按月收集城市o3浓度小时监测数据和全球数据同化预报系统(gdas)发布的气象数据；
15.s1.2以城市建成区环境空气质量监测站点的几何中心作为起始点，采用500m高度作为大气边界层的平均流场，通过hysplit模型对每天的气团移动24h后向轨迹进行模拟；
16.s1.3将评估城市及周边区域划分为5km
×
5km的网格，结合城市o3小时浓度以160μg/m3作为污染轨迹判断标准，即经过某网格的气团轨迹达到城市时所对应的o3浓度超过160μg/m3则判定该轨迹为污染轨迹，反之则判定为清洁轨迹；
17.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)月度数据集采样成城市及周边地区高空间分辨率的浓度分布数据，分别形成1km*1km分辨率的甲醛(hcho)和二氧化氮(no2)月度分布数据。
[0024]
s2.3根据步骤s1.5提取的潜在污染源网格，提取hcho和no2浓度数据，对潜在污染源网格内的每一个1km
×
1km网格hcho和no2浓度，通过比较网格浓度与背景浓度差异判定是否为热点网格，如满足下列条件(1)则判为no2热点网格，满足条件(2)则判为vocs热点网格：
[0025][0026][0027]
其中，νo2为1km网格no2柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格平均no2柱浓度，单位为10
13
molec/cm2；为5km潜在污染热点网格no2柱浓度的平均绝对偏差，单位为10
13
molec/cm2；hcho为1km网格hcho柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格平均hcho柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格hcho柱浓度的平均绝对偏差，单位为10
13
molec/cm2。
[0028]
s2.4分别以no2和vocs热点网格中心为原点，以800m为半径制作缓冲区，并合并形成no2热点区域和vocs热点区域。
[0029]
s3包括：
[0030]
s3.1所述电子地图包括高德地图、百度地图，根据多边形搜索方法，在no2热点区域内，以电厂、钢铁、冶金、铸造、焦化、锅炉、物流园区等关键字搜索公司企业类型的poi，收
集poi的地理坐标、名称等信息；
[0031]
s3.2所述电子地图包括高德地图、百度地图，根据多边形搜索方法，在vocs热点区域内，以4s、汽修、喷涂、涂料、装饰、家具、溶剂、油墨、印刷、加油站、石化、化工关键字搜索公司企业类型的poi，收集poi的地理坐标、名称信息；
[0032]
s3.3汇总热点区域搜索的poi信息，采用坐标转换方法将poi的火星坐标系(gcj02)转化为wgs84坐标系，形成疑似污染企业清单。
[0033]
s4包括：
[0034]
s4.1所述卫星遥感数据包括高分一号、高分二号、高分六号优于2m的多光谱卫星遥感数据，收集近一个月的多光谱卫星影像，并进行几何校正、图像融合等预处理；
[0035]
s4.2对疑似污染企业清单高分辨影像进行解译，剔除没有生产车间或者工厂面积小于100m2以及位于商场的企业poi干扰源；
[0036]
s4.3利用高分辨率卫星影像对剔除干扰poi后的污染企业无组织排放点等情况进行解译识别，无组织排放包括物料未密闭存储、露天作业等识别，存在10处以上无组织排放点判为高风险污染企业，5-10处无组织排放点判为中风险污染企业，5处以下判为低风险污染企业。
[0037]
本发明同现有技术相比，通过融合地面观测数据、后向轨迹、卫星遥感和电子地图poi数据构建了臭氧污染溯源技术方法及系统，准确定位污染源，列出污染清单，不仅能够解析臭氧污染成因，还能精准锁定污染来源，指导臭氧污染源精准管控工作，从而弥补传统的污染溯源排查方法中过于依赖污染源清单数据、覆盖面不全以及污染确定性不足导致无法准确锁定污染企业的局限性。
附图说明
[0038]
图1为本发明一种面向城市臭氧超标污染遥感溯源流程图
具体实施方式
[0039]
现结合附图对本发明做进一步描述。
[0040]
请参阅图1，一种面向城市臭氧超标污染遥感溯源方法，其特征在于，该方法包括步骤：
[0041]
s1.根据城市环境空气质量监测站点的o3浓度小时监测数据，利用潜在源贡献因子(pscf)方法和浓度权重轨迹分析(cwt)方法计算加权潜在源贡献因子(wpscf)和权重浓度，提取5km潜在污染源网格单元；
[0042]
s2.根据卫星遥感监测的区域甲醛(hcho)和二氧化氮(no2)柱浓度数据，重构合成高分辨率o3生成前体物分布结果，通过上下文方法，结合潜在污染源网格单元提取1km污染排放热点网格；
[0043]
s3.根据热点网格圈定矢量区域范围，利用关键字搜索方法从电子地图中获取污染排放热点网格单位内的相关企业poi；
[0044]
s4.利用高分辨率影像对的企业poi进行解译核查，筛选出污染企业清单；
[0045]
在本实施例中，步骤s1进一步包括：
[0046]
s1.1按月收集城市o3浓度小时监测数据和全球数据同化预报系统(gdas)发布的
气象数据；
[0047]
s1.2以城市建成区环境空气质量监测站点的几何中心作为起始点，采用500m高度作为大气边界层的平均流场，通过hysplit模型对每天的气团移动24h后向轨迹进行模拟；
[0048]
s1.3将评估城市及周边区域划分为5km
×
5km的网格，结合城市o3小时浓度以160μg/m3作为污染轨迹判断标准，即经过某网格的气团轨迹达到城市时所对应的o3浓度超过160μg/m3则判定该轨迹为污染轨迹，反之则判定为清洁轨迹；
[0049]
s1.4逐网格计算潜在源加权贡献因子(wpscf)，计算方式如下：
[0050][0051]
其中，wpscf
ij
为第i行j列网格的潜在源加权贡献因子，p
ij
为经过第i行j列网格的污染轨迹数量，n
ij
为经过第i行j列网格的所有轨迹数量。
[0052]
s1.5根据潜在源加权贡献因子(wpscf),筛选出潜在源加权贡献因子超过0.6的5km*5km网格单位作为潜在污染源网格。
[0053]
在本实施例中，步骤s2进一步包括：
[0054]
s2.1所述卫星遥感资料包括通过sentinel-5p卫星tropomi传感器探测数据反演获取的o3生成前体物甲醛(hcho)和二氧化氮(no2)柱浓度逐日产品数据，并从产品中提取每日5.5km
×
3.5km格点信息，包括格点经度、纬度、hcho柱浓度、no2柱浓度以及相应的qa值。
[0055]
s2.1按月将逐日的格点数据汇总合并，并根据qa值筛选出qa大于0.5的经度、纬度、hcho柱浓度和no2柱浓度月度监测数据集。根据月度监测数据，通过克里金插值方法将甲醛(hcho)和二氧化氮(no2)月度数据集采样成城市及周边地区高空间分辨率的浓度分布数据，分别形成1km*1km分辨率的甲醛(hcho)和二氧化氮(no2)月度分布数据。
[0056]
s2.3根据步骤s1.5提取的潜在污染源网格，提取hcho和no2浓度数据，对潜在污染源网格内的每一个1km
×
1km网格hcho和no2浓度，通过比较网格浓度与背景浓度差异判定是否为热点网格，如满足下列条件(1)则判为no2热点网格，满足条件(2)则判为vocs热点网格：
[0057][0058][0059]
其中，νo2为1km网格no2柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格平均no2柱浓度，单位为10
13
molec/cm2；为5km潜在污染热点网格no2柱浓度的平均绝对偏差，单位为10
13
molec/cm2；hcho为1km网格hcho柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格平均hcho柱浓度，单位为10
13
molec/cm2；为5km潜在污染源网格
hcho柱浓度的平均绝对偏差，单位为10
13
molec/cm2。
[0060]
s2.4分别以no2和vocs热点网格中心为原点，以800m为半径制作缓冲区，并合并形成no2热点区域和vocs热点区域。
[0061]
在本实施例中，步骤s3进一步包括：
[0062]
s3.1所述电子地图包括高德地图、百度地图等，根据多边形搜索方法，在no2热点区域内，以电厂、钢铁、冶金、铸造、焦化、锅炉、物流园区等关键字搜索公司企业类型的poi，收集poi的地理坐标、名称等信息；
[0063]
s3.2所述电子地图包括高德地图、百度地图等，根据多边形搜索方法，在vocs热点区域内，以4s、汽修、喷涂、涂料、装饰、家具、溶剂、油墨、印刷、加油站、石化、化工等关键字搜索公司企业类型的poi，收集poi的地理坐标、名称等信息；
[0064]
s3.3汇总热点区域搜索的poi信息，采用坐标转换方法将poi的火星坐标系(gcj02)转化为wgs84坐标系，形成疑似污染企业清单。
[0065]
在本实施例中，步骤s4进一步包括：
[0066]
s4.1所述卫星遥感数据包括高分一号、高分二号、高分六号等优于2m的多光谱卫星遥感数据，收集近一个月的多光谱卫星影像，并进行几何校正、图像融合等预处理；
[0067]
s4.2对疑似污染企业清单高分辨影像进行解译，剔除没有生产车间或者工厂面积小于100m2以及位于商场的企业poi干扰源；
[0068]
s4.3利用高分辨率卫星影像对剔除干扰poi后的污染企业无组织排放点等情况进行解译识别，无组织排放包括物料未密闭存储、露天作业等识别，存在10处以上无组织排放点判为高风险污染企业，5-10处无组织排放点判为中风险污染企业，5处以下判为低风险污染企业。
[0069]
需要说明的是，基于卫星遥感的o3超标污染溯源方法是一种综合溯源方法,它是一种反映区域o3前体物排放等大气污染源时空分布特征、区域传输影响和污染源无组织排放等污染溯源的量化指标，反映了o3污染来源成因及污染排放企业空间分布特征。因此，采用本方法客服了传统溯源证据不足或者覆盖不全面的劣势，使得通过基于卫星遥感的o3污染溯源更加精准全面。
[0070]
以上仅是本发明的优选实施方式，只是用于帮助理解本技术的方法及其核心思想，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
[0071]
本发明从整体上解决了现有技术中只是基于有限站点监测数据的方向性溯源，没有考虑到o3生成前体物污染排放源动态变化、区域污染传输的问题，通过融合地面观测数据、后向轨迹、卫星遥感和电子地图poi数据构建了臭氧污染溯源技术方法及系统，准确定位污染源，列出污染清单，不仅能够解析臭氧污染成因，还能精准锁定污染来源，指导臭氧污染源精准管控工作，从而弥补传统的污染溯源排查方法中过于依赖污染源清单数据、覆盖面不全以及污染确定性不足导致无法准确锁定污染企业的局限性。