基于细颗粒物与臭氧的空气污染源区确定方法

摘要

本公开提供了一种基于细颗粒物与臭氧的空气污染源区确定方法，包括：获取空气污染源区确定任务；执行空气污染源区确定任务，以获取研究区域后向轨迹数据；将获取的后向轨迹数据网格化；导入研究区域的PM2.5和O3的逐时观测值；根据PM2.5对应的观测值确定PM2.5对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；根据O3对应的观测值确定O3对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；根据PM2.5和O3分别对应的PSCFij以及PM2.5和O3分别对应的CWTij分别确定满足预设条件的多个区域；确定多个区域中的公共区域；将公共区域确定为空气污染源区。

权利要求

1.一种基于细颗粒物与臭氧的空气污染源区确定方法，包括：获取空气污染源区确定任务；

执行所述空气污染源区确定任务，以获取研究区域后向轨迹数据；

将获取的所述后向轨迹数据网格化；

导入所述研究区域的PM2.5和O3的逐时观测值；

根据所述PM2.5对应的观测值确定所述PM2.5对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；

根据所述O3对应的观测值确定所述O3对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；其中，i和j分别代表经度、纬度；

根据所述PM2.5和所述O3分别对应的潜在源贡献因子PSCFij以及所述PM2.5和所述O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域；

确定所述多个区域中的公共区域；

将所述公共区域确定为空气污染源区；

其中，所述满足预设条件的多个区域包括：

所述PM2.5对应的潜在源贡献因子PSCFij大于第一预设阈值的第一区域；

所述O3对应的潜在源贡献因子PSCFij大于第二预设阈值的第二区域；

所述PM2.5对应的浓度权重轨迹CWTij大于第三预设阈值的第三区域；

所述O3对应的浓度权重轨迹CWTij大于第四预设阈值的第四区域；

其中，所述确定所述多个区域中的公共区域包括：确定所述第一区域、所述第二区域、所述第三区域和所述第四区域的第一公共区域，其中，所述第一公共区域作为协同控制污染外来源区；

确定所述第一区域和所述第三区域的第二公共区域，其中，所述第二公共区域作为PM2.5污染外来源区；

确定所述第二区域和所述第四区域的第三公共区域，其中，所述第三公共区域作为O3污染外来源区。

2.根据权利要求1所述的方法，

其中，所述第一预设阈值和所述第二预设阈值均大于等于0.5。

3.根据权利要求1所述的方法，

3 3

其中，所述第三预设阈值大于等于35μg/m；所述第四预设阈值大于等于160μg/m。

4.根据权利要求1所述的方法，还包括：

根据所述协同控制污染外来源区包含的城市的FNR值，确定所述协同控制污染外来源区的PM2、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区；

根据所述O3污染外来源区包含的城市的FNR值，确定所述O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区；

其中，所述FNR值表示HCHO的浓度与NO2/NOv污染物浓度的比值。

5.根据权利要求4所述的方法，所述根据所述协同控制污染外来源区包含的城市的FNR值，确定所述协同控制污染外来源区的PM2.5、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区包括：所述协同控制污染外来源区包含的城市的FNR值≤2.3的情况下，确定为PM2.5、VOCs协同控制区；FNR值≥4.2的情况下，确定为PM2.5、NOx协同控制区；2.3＜FNR＜4.2的情况下，确定为PM2.5、NOx、VOCs协同控制区；

所述根据所述O3污染外来源区包含的城市的FNR值，确定所述O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区包括：所述O3污染外来源区包含的城市的FNR值≤2.3的情况下，确定为VOCs主控区；FNR值≥

4.2的情况下，确定为NOx主控区；2.3＜FNR＜4.2的情况下，确定为NOx、VOCs协同控制区。

6.根据权利要求1所述的方法，在所述确定满足预设条件的多个区域之前还包括：引入权重因子W(nij)，将权重因子W(nij)分别与所述潜在源贡献因子PSCFij和所述浓度权重轨迹CWTij相乘，得到修正潜在源贡献因子WPSCFij和修正浓度权重轨迹WCWTij；其中，nij表示落在网格(i，j)中的端点数量；

根据所述PM2.5和所述O3分别对应的潜在源贡献因子PSCFij以及所述PM2.5和所述O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域包括：根据所述PM2.5和所述O3分别对应的修正潜在源贡献因子WPSCFij以及所述PM2.5和所述O3分别对应的修正浓度权重轨迹WCWTij分别确定满足预设条件的多个区域。

7.根据权利要求6所述的方法，

所述权重因子W(nij)包括：

8.根据权利要求1～7中任一所述的方法，

所述潜在源贡献因子PSCFij的计算方法包括：其中，nij表示落在网格(i，j)中的端点数量；

在计算PM2.5对应的潜在源贡献因子PSCFij时，mij表示为网格(i，j)中，超过所述第一预设阈值的端点数量；

在计算O3对应的潜在源贡献因子PSCFij时，mix表示为网格(i，j)中，超过所述第二预设阈值的端点数量；

网格(i，j)表示经度i和纬度j围成的网格；

所述浓度权重轨迹CWTij的计算方法包括：

其中，CWTij为网格(i，j)的平均污染权重浓度；k为后向轨迹；N为经过网格(i，j)的轨迹总数；Ck为轨迹k经过网格(i，j)时对应的污染物浓度，αijk为轨迹k在网格(i，j)停留的时间。

9.根据权利要求8所述的方法，所述将获取的所述后向轨迹数据网格化包括：将所述后向轨迹数据网格化为分辨率为1°×1°的网格。

说明书

基于细颗粒物与臭氧的空气污染源区确定方法

技术领域

[0001] 本公开涉及空气质量监测技术领域，尤其涉及一种基于细颗粒物与臭氧的空气污染源区确定方法。

背景技术

[0002] 由于空气污染的严峻形势，国内外研究者针对大气污染物的浓度分布、时空变化、区域输送及形成机理等方面开展了大量的研究。各地也积极开展空气污染来源解析、空气质量预报业务和能力建设，以对大气污染进行来源解析并分析识别污染贡献重点区，为污染控制提供依据。

[0003] 目前，虽然细颗粒物(PM2.5)浓度持续下降，但是臭氧(O3)浓度呈快速上升和蔓延态势。研究表明，虽然PM2.5和O3属于不同类型的污染物，但二者之间可以通过多种途径相互作用，以PM2.5防控为重心的减排措施所导致的污染源排放结构的变化，可能是臭氧污染持续上升的重要原因。因此，基于PM2.5与O3的污染传输关系，来准确地确定空气污染来源的主要源区，不仅能为日常生活和生产活动提供指导，还能针对性地进行污染防控。

发明内容

[0004] 针对上述技术问题，本公开提供了一种基于细颗粒物与臭氧的空气污染源区确定方法，以至少部分地解决上述技术问题中的至少之一。

[0005] 为了解决上述技术问题，本公开提供了一种基于细颗粒物与臭氧的空气污染源区确定方法，包括：

[0006] 获取空气污染源区确定任务；

[0007] 执行所述空气污染源区确定任务，以获取研究区域后向轨迹数据；

[0008] 将获取的所述后向轨迹数据网格化；

[0009] 导入所述研究区域的PM2.5和O3的逐时观测值；

[0010] 根据所述PM2.5对应的观测值确定所述PM2.5对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；

[0011] 根据所述O3对应的观测值确定所述O3对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；其中，i和j分别代表经度、纬度；

[0012] 根据所述PM2.5和所述O3分别对应的潜在源贡献因子PSCFij以及所述PM2.5和所述O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域；

[0013] 确定所述多个区域中的公共区域；

[0014] 将所述公共区域确定为空气污染源区。

[0015] 根据本公开的实施例，所述满足预设条件的多个区域包括：

[0016] 所述PM2.5对应的潜在源贡献因子PSCFij大于第一预设阈值的第一区域；

[0017] 所述O3对应的潜在源贡献因子PSCFij大于第二预设阈值的第二区域；

[0018] 所述PM2.5对应的浓度权重轨迹CWTij大于第三预设阈值的第三区域；

[0019] 所述O3对应的浓度权重轨迹CWTij大于第四预设阈值的第四区域；

[0020] 其中，所述确定所述多个区域中的公共区域包括：

[0021] 确定所述第一区域、所述第二区域、所述第三区域和所述第四区域的第一公共区域，其中，所述第一公共区域作为协同控制污染外来源区；

[0022] 确定所述第一区域和所述第三区域的第二公共区域，其中，所述第二公共区域作为PM2.5污染外来源区；

[0023] 确定所述第二区域和所述第四区域的第三公共区域，其中，所述第三公共区域作为O3污染外来源区。

[0024] 根据本公开的实施例，所述第一预设阈值和所述第二预设阈值均大于等于0.5。

[0025] 根据本公开的实施例，所述第三预设阈值大于等于35μg/m3；所述第四预设阈值大3
于等于160μg/m。

[0026] 根据本公开的实施例，上述方法还包括：

[0027] 根据所述协同控制污染外来源区包含的城市的FNR值，确定所述协同控制污染外来源区的PM2、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区；

[0028] 根据所述O3污染外来源区包含的城市的FNR值，确定所述O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区；

[0029] 其中，所述FNR值表示HCHO的浓度与NO2/NOy等污染物浓度的比值。

[0030] 根据本公开的实施例，所述根据所述协同控制污染外来源区包含的城市的FNR值，确定所述协同控制污染外来源区的PM2.5、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区包括：

[0031] 所述协同控制污染外来源区包含的城市的FNR值≤2.3的情况下，确定为PM2.5、VOCs协同控制区；FNR值≥4.2的情况下，确定为PM2.5、NOx协同控制区；2.3

[0032] 所述根据所述O3污染外来源区包含的城市的FNR值，确定所述O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区包括：

[0033] 所述O3污染外来源区包含的城市的FNR值≤2.3的情况下，确定为VOCs主控区；FNR值≥4.2的情况下，确定为NOx主控区；2.3

[0034] 根据本公开的实施例，在所述确定满足预设条件的多个区域之前还包括：

[0035] 引入权重因子W(nij)，将权重因子W(nij)分别与所述潜在源贡献因子PSCFij和所述浓度权重轨迹CWTij相乘，得到修正潜在源贡献因子WPSCFij和修正浓度权重轨迹WCWTij；其中，nij表示落在网格(i，j)中的端点数量；

[0036] 根据所述PM2.5和所述O3分别对应的潜在源贡献因子PSCFij、所述PM2.5和所述O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域包括：

[0037] 根据所述PM2.5和所述O3分别对应的修正潜在源贡献因子WPSCFij、所述PM2.5和所述O3分别对应的修正浓度权重轨迹WCWTij分别确定满足预设条件的多个区域。

[0038] 根据本公开的实施例，所述权重因子W(nij)包括：

[0039]

[0040] 根据本公开的实施例，所述潜在源贡献因子PSCFij的计算方法包括：

[0041]

[0042] 其中，nij表示落在网格(i，j)中的端点数量；

[0043] 在计算PM2.5对应的潜在源贡献因子PSCFij时，mij表示为网格(i，j)中，超过所述第一预设阈值的端点数量；

[0044] 在计算O3对应的潜在源贡献因子PSCFij时，mij表示为网格(i，j)中，超过所述第二预设阈值的端点数量；

[0045] 网格(i，j)表示经度i和纬度j围成的网格；

[0046] 所述浓度权重轨迹CWTij的计算方法包括：

[0047]

[0048] 其中，CWTij为网格(i，j)的平均污染权重浓度；k为后向轨迹；N为经过网格(i，j)的轨迹总数；Ck为轨迹k经过网格(i，j)时对应的污染物浓度，dijk为轨迹k在网格(i，j)停留的时间。

[0049] 根据本公开的实施例，所述将获取的所述后向轨迹数据网格化包括：将所述后向轨迹数据网格化为分辨率为1°×1°的网格。

[0050] 根据本公开的实施例，通过将后向轨迹数据网格化，并导入研究区PM2.5和O3的逐时观测值，然后根据PM2.5和O3对应的逐时观测值分别确定PM2.5和O3分别对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；基于PM2.5和O3分别对应的潜在源贡献因子PSCFij以及PM2.5和O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域，从而确定多个区域的公共区域，作为空气污染源区。因此，根据PM2.5和O3的污染传输关系，将城市按照污染的复合程度进行分类，可以为地方在选择符合各地污染特点的大气污染防治措施时提供判断依据，同时有助于研究人员在进行复合污染的研究时，在目标城市的选择上可以贴近实际需求。

[0051] 另外，采用本公开提供的空气污染源区的确定方法，操作简便、准确率较高，能够实现对区域的大气污染来源解析、对识别污染贡献重点区的论证，进而实现对污染贡献重点区的识别，为研究大气污染治理分区管控提供有力支撑。

附图说明

[0052] 图1示意性示出了本公开实施例中基于PM2.5与O3的空气污染源区确定方法的流程图；

[0053] 图2示意性示出了本公开另一实施例中基于细颗粒物与臭氧的空气污染源区确定方法的流程图；

[0054] 图3示意性示出了基于FNR值确定协同控制污染外来源区类型的方法的流程图；

[0055] 图4示意性示出了基于FNR值确定O3污染外来源区类型的方法的流程图。

具体实施方式

[0056] 为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开作进一步的详细说明。

[0057] PM2.5和O3相互作用机制存在较多争议，对流层O3和二次PM2.5的产生均来自NOx和VOCs的光化学氧化过程，但转化机制的相关研究还较为薄弱。目前，联防联控机制在解决区域性大气污染问题中，发挥重要作用。尤其在急需解决PM2.5与O3污染，而化学机制又不明的情况下，可以充分发挥联防联控机制的作用，梳理城市间污染物传输关系，厘清城市的污染类型，按照污染的复合程度进行分类，有利于联防联控政策的制定，以及指导地方进行符合各地污染特点的大气污染防治措施，同时有助于研究人员选择性的进行研究。

[0058] 拉格朗日混合单粒子轨道模型，结合潜在源贡献因子法(Potential Source Contribution Function，PSCF)，浓度权重轨迹分析法(Concentration‑Weighted Trajectory，CWT)由于操作简便，可在不依赖排放清单的前提下，进行大气污染外来源区的研究。

[0059] 图1示意性示出了本公开实施例中基于细颗粒物与臭氧的空气污染源区确定方法的流程图。

[0060] 如图1，根据本公开的实施例，提供一种基于细颗粒物与臭氧的空气污染源区确定方法，包括操作S101～S109。

[0061] 在操作S101中，获取空气污染源区确定任务。

[0062] 在操作S102中，执行空气污染源区确定任务，以获取研究区域后向轨迹数据。

[0063] 根据本公开的实施例，获取研究区域的后向轨迹数据例如可以包括：采用MeteoInfo软件中的Trajectory Calculation功能(用于计算轨迹)进行后向轨迹计算，其中，MeteoInfo是一个气象数据显示、分析软件。根据本公开的实施例，可以选择研究区域坐标作为后向轨迹起点，模拟研究区域的PM2.5和O3在距离地面预设高度(如500米高度)，预设时间范围(如24h)内的后向轨迹，每次模拟预设时长(如1个月)的后向轨迹，例如，可以每天逐时模拟后向轨迹，直到模拟获得预设时长的后向轨迹。根据研究所需时间段(如季节、年份)将模拟的后向轨迹合并形成后向轨迹数据。

[0064] 后向轨迹是空气微团在大气中移动的实际路径，气团轨迹由其移动速度确定。气团轨迹可由风向风速资料进行计算，亦可通过随动“标记粒子”的观测而得到。气团轨迹可用于了解污染物平均搬运路径和烟团中心轴线的确定。

[0065] 在操作S103中，将获取的后向轨迹数据网格化。

[0066] 在操作S104中，导入研究区域的PM2.5和O3的逐时观测值。

[0067] 在操作S105中，根据PM2.5对应的观测值确定PM2.5对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij。

[0068] 在操作S106中，根据O3对应的观测值确定O3对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；其中，i和j分别代表经度、纬度。

[0069] 在操作S107中，根据PM2.5和O3分别对应的潜在源贡献因子PSCFij以及PM2.5和O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域。

[0070] 在操作S108中，确定多个区域中的公共区域。

[0071] 在操作S109中，将公共区域确定为空气污染源区。

[0072] 根据本公开的实施例，通过将后向轨迹数据网格化，并导入研究区PM2.5和O3的逐时观测值，然后根据PM2.5和O3对应的逐时观测值分别确定PM2.5和O3分别对应的潜在源贡献因子PSCFij和浓度权重轨迹CWTij；基于PM2.5和O3分别对应的潜在源贡献因子PSCFij以及PM2.5和O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域，从而确定多个区域的公共区域，作为空气污染源区。因此，本公开提供的方法根据PM2.5和O3的污染传输关系，按照污染的复合程度进行分类，可以为地方在选择符合各地污染特点的大气污染防治措施时提供判断依据，同时有助于研究人员在进行复合污染的研究时，在目标城市的选择上可以贴近实际需求。

[0073] 另外，采用本公开提供的空气污染源区的确定方法，操作简便、准确率较高，能够实现对区域的大气污染来源解析、对识别污染贡献重点区的论证，进而实现对污染贡献重点区的识别，为研究大气污染治理分区管控提供有力支撑。

[0074] 根据本公开的实施例，逐时观测值来自空气质量实时发布平台。

[0075] 根据本公开的实施例，满足预设条件的多个区域包括：

[0076] PM2.5对应的潜在源贡献因子PSCFij大于第一预设阈值的第一区域；

[0077] O3对应的潜在源贡献因子PSCFij大于第二预设阈值的第二区域；

[0078] PM2.5对应的浓度权重轨迹CWTij大于第三预设阈值的第三区域；

[0079] O3对应的浓度权重轨迹CWTij大于第四预设阈值的第四区域；

[0080] 根据本公开的实施例，确定多个区域中的公共区域包括：

[0081] 确定第一区域、第二区域、第三区域和第四区域的第一公共区域，其中，第一公共区域作为协同控制污染外来源区；

[0082] 确定第一区域和第三区域的第二公共区域，其中，第二公共区域作为PM2.5污染外来源区；

[0083] 确定第二区域和第四区域的第三公共区域，其中，第三公共区域作为O3污染外来源区。

[0084] 根据本公开的实施例，第一预设阈值和第二预设阈值均大于等于0.5。

[0085] 根据本公开的实施例，第一预设阈值和第二预设阈值的预设范围可以根据实际应用进行设置，本公开不对其进行限定。例如，第一预设阈值可以为0.5～1.5，具体地可以为0.5、0.6、0.7等；第二预设阈值可以为0.5～2.5，具体地可以为0.5、0.6、0.7、1.5等。

[0086] 根据本公开的实施例，第一预设阈值和第二预设阈值可以相同，也可以不同。

[0087] 根据本公开的实施例，第三预设阈值大于等于35μg/m3；第四预设阈值大于等于3
160μg/m。

[0088] 根据本公开的实施例，第三预设阈值和第四预设阈值的预设范围可以根据实际应3
用进行设置，本公开不对其进行限定。例如，第三预设阈值可以为25～50μg/m ，具体地可以
3 3 3 3
为35μg/m 、38μg/m、40μg/m 等；第四预设阈值可以为160～250μg/m ，具体地可以为160μg/
3 3 3
m、170μg/m、180μg/m等。

[0089] 图2示意性示出了本公开另一实施例中基于细颗粒物与臭氧的空气污染源区确定方法的流程图。

[0090] 如图2，根据本公开的实施例，提供一种基于细颗粒物与臭氧的空气污染源区确定方法，包括操作S201～S210。

[0091] 操作S201～S206与图1中的操作S101～S106相同，在此不再赘述。

[0092] 在操作S207，筛选PM2.5对应的PSCFij大于第一预设阈值的第一区域；筛选O3对应的PSCFij大于第二预设阈值的第二区域；筛选PM2.5对应的CWTij大于第三预设阈值的第三区域；筛选O3对应的CWTij大于第四预设阈值的第四区域。

[0093] 在操作S208，确定第一区域、第二区域、第三区域和第四区域的第一公共区域作为协同控制污染外来源区。

[0094] 在操作S209，确定第一区域和第三区域的第二公共区域作为PM2.5污染外来源区。

[0095] 在操作S207，确定第二区域和第四区域的第三公共区域作为O3污染外来源区。

[0096] 根据本公开的实施例，上述空气污染源区确定方法还包括：根据协同控制污染外来源区包含的城市的FNR值，确定协同控制污染外来源区的PM2、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区。

[0097] 根据本公开的实施例，上述空气污染源区确定方法还包括：根据O3污染外来源区包含的城市的FNR值，确定O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区；其中，FNR值表示HCHO的浓度与NO2/NOy等污染物浓度的比值。

[0098] HCHO作为多种VOCs的短寿命氧化产物，可用于指示VOCs的反映速率，其与NO2等污染物浓度的比值FNR被广泛用于判断O3生成敏感性。FNR指示剂法基于卫星观测数据，测量并计算得到目标区域FNR值。根据区间范围判定O3生成的控制类型。

[0099] 根据本公开的实施例，根据协同控制污染外来源区包含的城市的FNR值，确定协同控制污染外来源区的PM2.5、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区包括：协同控制污染外来源区包含的城市的FNR值≤2.3的情况下，确定为PM2.5、VOCs协同控制区；FNR值≥4.2的情况下，确定为PM2.5、NOx协同控制区；2.3

[0100] 根据本公开的实施例，根据O3污染外来源区包含的城市的FNR值，确定O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区包括：O3污染外来源区包含的城市的FNR值≤2.3的情况下，确定为VOCs主控区；FNR值≥4.2的情况下，确定为NOx主控区；2.3

[0101] 图3示意性示出了基于FNR值确定协同控制污染外来源区类型的方法的流程图。

[0102] 如图3，该方法包括操作S301～S305。

[0103] 在操作S301，确定协同控制污染外来源区。

[0104] 在操作S302，根据协同控制污染外来源区包含的城市的FNR值，确定协同控制污染外来源区的PM2.5、VOCs协同控制区，PM2.5、NOx协同控制区或PM2.5、NOx、VOCs协同控制区。

[0105] 在操作S303，协同控制污染外来源区包含的城市的FNR值≤2.3的情况下，确定为PM2.5、VOCs协同控制区。

[0106] 在操作S303，协同控制污染外来源区包含的城市的FNR值≥4.2的情况下，确定为PM2.5、NOx协同控制区。

[0107] 在操作S305，协同控制污染外来源区包含的城市的FNR值在2.3

[0108] 图4示意性示出了基于FNR值确定O3污染外来源区类型的方法的流程图。

[0109] 如图4，该方法包括操作S401～S405。

[0110] 在操作S401，确定O3污染外来源区。

[0111] 在操作S402，根据O3污染外来源区包含的城市的FNR值，确定O3污染外来源区的VOCs主控区、NOx主控区或NOx、VOCs协同控制区。

[0112] 在操作S403，O3污染外来源区包含的城市的FNR值≤2.3的情况下，确定为VOCs主控区。

[0113] 在操作S404，O3污染外来源区包含的城市的FNR值≥4.2的情况下，确定为NOx主控区。

[0114] 在操作S405，O3污染外来源区包含的城市的FNR值在2.3

[0115] 根据本公开的实施例，在确定满足预设条件的多个区域之前还包括：

[0116] 引入权重因子W(nij)，将权重因子W(nij)分别与潜在源贡献因子PSCFij和浓度权重轨迹CWTij相乘，得到修正潜在源贡献因子WPSCFij和修正浓度权重轨迹WCWTij；其中，nij表示落在网格(i，j)中的端点数量；

[0117] 根据PM2.5和O3分别对应的潜在源贡献因子PSCFij以及PM2.5和O3分别对应的浓度权重轨迹CWTij分别确定满足预设条件的多个区域包括：根据PM2.5和O3分别对应的修正潜在源贡献因子WPSCFij以及PM2.5和O3分别对应的修正浓度权重轨迹WCWTij分别确定满足预设条件的多个区域。

[0118] 根据本公开的实施例，权重因子W(nij)包括：

[0119]

[0120] 根据本公开的实施例，潜在源贡献因子PSCFij的计算方法包括：

[0121]

[0122] 其中，nij表示落在网格(i，j)中的端点数量；

[0123] 在计算PM2.5对应的潜在源贡献因子PSCFij时，mij表示为网格(i，j)中，超过第一预设阈值的端点数量；

[0124] 在计算O3对应的潜在源贡献因子PSCFij时，mij表示为网格(i，j)中，超过第二预设阈值的端点数量；

[0125] 网格(i，j)表示经度i和纬度j围成的网格；

[0126] 浓度权重轨迹CWTij的计算方法包括：

[0127]

[0128] 其中，CWTij为网格(i，j)的平均污染权重浓度；k为后向轨迹；N为经过网格(i，j)的轨迹总数；Ck为轨迹k经过网格(i，j)时对应的污染物浓度，αijk为轨迹k在网格(i，j)停留的时间。

[0129] 根据本公开的实施例，将获取的后向轨迹数据网格化包括：将后向轨迹数据网格化为分辨率为1°×1°的网格。

[0130] 为了进一步理解本公开，以具体案例作进一步说明：

[0131] 采用本公开提供的方法分别对城市1的PM2.5和O3污染外来源区进行分析，污染物数据来源于空气质量环境监测平台，时间范围为2017年3月‑2020年2月。

[0132] 其中，时间划分为春季是3～5月、夏季是6～8月、秋季是9～11月、冬季是12月～次年2月。

[0133] 城市1的PM2.5与O3的污染主要外来源区类别，以及O3污染的控制类型，使用不同颜色可视化表示。

[0134] 城市1的PM2.5和O3的季节性变化趋势基本相反，PM2.5的浓度ρ(PM2.5)四季变化趋势3
呈“U”形，夏季PM2.5污染最轻，冬季PM2.5污染最严重，平均值达64μg/m 。O3的浓度ρ(O3)四季变化趋势呈“n”型，夏季O3污染最严重，冬季O3污染最轻。因此，对于城市1，春、秋两季需同时控制PM2.5与O3、夏季重点控制O3、冬季重点控制PM2.5。

[0135] 城市1西南方向部分地区在各个季节对城市1的PM2.5污染有较大贡献，在春季A省南部、B省北部、C省西部对城市1的PM2.5污染贡献较大；在秋季对城市1的PM2.5污染贡献较大的区域集中在A省南部、B省北部；冬季对城市1的PM2.5污染贡献较大的区域明显较大，连片出现A省中部与南部、B省北部、C省西部。城市2、城市3、城市4、城市5、城市6在春、秋、冬季均为PM2.5对应的WPSCF、WCWT高值区所覆盖的城市，上述城市应为城市1的PM2.5联防联控重点关注的城市。

[0136] 城市1的O3污染主要集中在夏季，夏季对城市A的O3污染贡献较大的区域明显最大，连片出现在A省南部，B省北部，C省部分城市；春季O3污染相对较轻，贡献明显的区域在A省东南部分、B省中部与北部、C省西部与北部。城市2、城市4、城市5、城市7、城市8、城市9在春、夏两季均为O3对应的WPSCF、WCWT高值区所覆盖的城市，应为城市1的O3联防联控重点关注城市。引入城市2、城市4、城市5、城市7、城市8、城市9的FNR值，根据FNR值的预设阈值可分为VOCs主控城市，NOx主控城市，VOCs、NOx联防联控城市。

[0137] 另外，存在部分城市对城市1的PM2.5与O3污染均产生较大影响，春季为城市2、城市4、城市5、城市6、城市7、城市10，秋季为城市7。引入城市2、城市4、城市5、城市6、城市7、城市
10的FNR值，根据FNR值的预设阈值划分城市2、城市4、城市5为VOCs主控城市，城市6、城市7、城市10为VOCs、NOx联防联控城市。

[0138] 以上的具体实施例，对本公开的目的、技术方案和有益效果进行了进一步详细说明，应理解的是，以上仅为本公开的具体实施例而已，并不用于限制本公开，凡在本公开的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。