一种基于大气污染浓度变化率的突出值站点识别方法技术

本发明专利技术公开了一种基于大气污染浓度变化率的突出值站点识别方法，属于大气污染防治技术领域；包括：基于各站点长时间序列历史浓度数据，计算逐时刻的站点浓度平均值；基于逐时刻平均值，依据各污染物浓度标准，划分逐时刻污染等级；逐站点计算研究时刻相较上一时刻的浓度变化率；将逐时刻的浓度变化率与划分的逐时刻污染等级相匹配，得到不同污染等级的时刻变化率数组，分别绘制概率分布图；依据概率分布图，分别确定变化率上限值，作为不同污染等级的突出值识别标准，即：标准变化率限值；基于标准限值，实时动态识别突出值站点与时刻。本发明专利技术基于浓度变化率实时动态识别需要重点管控的站点和管控时刻，可为面向监测站点的污染精细化管控提供技术支撑。精细化管控提供技术支撑。精细化管控提供技术支撑。

【技术实现步骤摘要】
一种基于大气污染浓度变化率的突出值站点识别方法


[0001]本专利技术涉及大气污染防治
，具体涉及一种基于大气污染浓度变化率的突出值站点识别方法。

技术介绍

[0002]大气污染对人类健康、社会经济发展有重要负面影响，大气污染治理具有重要意义。随着我国经济快速发展和工业化、城镇化进程不断加速，城市污染问题凸显，面向城市监测站点的污染精细化管控日渐成为重要的污染防治手段。而要实现站点精细化管控，就必须明确需要重点管控的站点和突出值产生的时刻，识别站点的特异性，进而实现快速精准管控。因此，需要建立一套实时动态识别突出值站点和对应突出值时刻的方法。
[0003]小尺度站点级别的气象条件在大范围内差别不大，而不同监测站点周边的污染源数量和种类却大不相同，对监测站点污染物浓度贡献也不同；因此，站点监测浓度的差异与人为排放（周边源）密切相关，对站点突出值的识别就需要综合考虑站点周边源影响的差异性。此外，在面向管理需求时，实时、动态、准确的识别突出值站点以及突出值时刻是当前站点精细化管控急需解决的难题。
[0004]目前常用的突出值识别方法主要以某一浓度阈值（如：国家二级标准等）或站点浓度平均值为标准，超过阈值或平均值即识别为突出值站点和对应的突出值时段，将站点浓度作为判别标准，周边源数量多的站点浓度通常较高，稳定排放源（如：机动车源早晚高峰）会在固定时刻浓度贡献较高，得到的突出值站点和时刻并未识别出站点的特异性，并且识别的突出值数量众多，人力物力消耗巨大，无法实现高效管控。

技术实现思路

[0005]针对现有技术中存在的不足之处，本专利技术提供一种基于大气污染浓度变化率的突出值站点识别方法，其综合考虑站点周边源排放的差异性，对突出值站点和对应的突出值时段进行实时动态识别，进而为采取有针对性的管控措施提供技术支撑，这对有效改善空气质量意义重大。
[0006]本专利技术公开了一种基于大气污染浓度变化率的突出值站点识别方法，包括：步骤1、基于各站点各污染物的长时间序列历史浓度数据，计算研究时段的逐时刻站点浓度平均值；步骤2、基于逐时刻站点浓度平均值，依据各污染物浓度标准，划分逐时刻污染等级；步骤3、逐站点计算当前时刻相较上一时刻的浓度变化率，将得到的逐时刻浓度变化率与划分的逐时刻污染等级相匹配，得到不同污染等级的时刻变化率数组，并对不同污染等级的时刻变化率数组进行概率拟合，绘制不同污染等级的变化率概率分布图；步骤4、依据不同污染等级的变化率概率分布图，分别确定其浓度变化率上限值，即标准变化率限值；
步骤5、基于不同污染等级的标准变化率限值，实时动态识别突出值站点与对应的突出值时刻；步骤6、依照不同地区不同污染物的时间分布特征对研究时刻进行时段划分，再基于得到的历史突出值站点及其出现突出值的时刻，对突出值站点和划分时段的突出值进行统计分析，得到重点管控站点以及重点管控时段。
[0007]作为本专利技术的进一步改进，在所述步骤1中，不同污染物研究时段内每t时刻的站点浓度平均值的计算公式为：式中，为t时刻的站点浓度平均值；为i站点t时刻的污染物浓度；t为时刻，其时间尺度包括分钟、小时、日和月中的一种；I为站点个数。
[0008]作为本专利技术的进一步改进，所述步骤2，具体包括：基于逐时刻站点浓度平均值，依照污染物相关浓度标准（具体标准依据实际情况而定），将污染物的污染等级划分为多个等级，并得到逐时刻污染等级，统称为l
t
。
[0009]作为本专利技术的进一步改进，所述步骤3，具体包括：步骤31、计算各站点污染物每t时刻相较于t
‑
1时刻的浓度变化率：式中，为i站点t时刻的污染物浓度变化率；t为时刻，其时间尺度包括分钟、小时、日和月中的一种；、分别为t时刻、t
‑
1时刻的污染物浓度；步骤32、将逐时刻的浓度变化率与划分的逐时刻污染等级l
t
相匹配，得到不同污染等级对应的M组时刻变化率数组；步骤33、基于不同污染等级对应的M组时刻变化率数组，分别进行概率拟合，计算得到不同污染等级对应的M组时刻变化率数组未重复值的概率f后，分别绘制不同污染等级的变化率直方图。
[0010]作为本专利技术的进一步改进，在步骤33中，拟合方法为：分别计算M组数据去除重复值后的概率：式中，f为l等级数组未重复值的概率；
为l等级数组的未重复值；为l等级数组中未重复值的频数；为l等级数组的浓度变化率总样本数；l为污染等级。
[0011]作为本专利技术的进一步改进，所述步骤4具体包括：基于不同污染等级的变化率概率分布图，计算不同污染等级的平均值μ和标准偏差σ，得到站点不同污染等级的浓度变化率上限值Q
l
，将Q
l
作为不同污染等级的突出值识别标准，即标准变化率限值；式中，Q
l
为l等级的标准变化率限值；为l等级的平均值和标准差；系数1,2,3的选择依据实际需求而定。
[0012]作为本专利技术的进一步改进，所述步骤5，具体包括：步骤51、确定站点研究时刻的污染等级，计算站点研究时刻的变化率，与对应污染等级的标准变化率限值进行对比；步骤52、超出标准变化率限值，则判定为突出值站点，并识别突出值站点与对应的突出值时刻。
[0013]作为本专利技术的进一步改进，所述步骤6，具体包括：步骤61、依照不同地区不同污染物的时间分布特征对研究时刻进行时段划分，划分为时段1、时段2、
……
等多个时段；步骤62、基于得到的研究时段的历史突出值站点以及其对应突出值的时刻，对突出值站点和划分的各时段内的突出值进行概率统计，出现突出值频率最高的站点为重点管控站点，出现突出值频率最高的时段为重点管控时段。
[0014]与现有技术相比，本专利技术的有益效果为：本专利技术通过对站点不同污染等级的污染物浓度变化率进行概率拟合确定标准变化率限值，综合考虑不同站点、不同污染物、不同污染等级周边源排放的差异性导致的浓度变化率差异以及管理需求的时效性，实时动态识别突出值站点和突出值时段并统计历史数据得到突出值出现的高频站点和高频时段，得到的突出值相对有效，研究成果可为管理部门提供重要的参考依据，为站点的靶向性管控提供科学支撑，对精细化管控效率的提升有着重要作用。
附图说明
[0015]图1为本专利技术公开的基于大气污染浓度变化率的突出值站点识别方法的流程图；图2 临沂市国控站点分布图；图3 站点不同污染等级概率分布图；图4 站点1的PM2.5突出值识别结果；
图5 站点2的PM2.5突出值识别结果。
具体实施方式
[0016]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0017]下面结合附图对本专利技术做进一步的详细描述：如本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于大气污染浓度变化率的突出值站点识别方法，其特征在于，包括：步骤1、基于各站点各污染物的长时间序列历史浓度数据，计算研究时段的逐时刻站点浓度平均值；步骤2、基于逐时刻站点浓度平均值，依据各污染物浓度标准，划分逐时刻污染等级；步骤3、逐站点计算当前时刻相较上一时刻的浓度变化率，将得到的逐时刻浓度变化率与划分的逐时刻污染等级相匹配，得到不同污染等级的时刻变化率数组，并对不同污染等级的时刻变化率数组进行概率拟合，绘制不同污染等级的变化率概率分布图；步骤4、依据不同污染等级的变化率概率分布图，分别确定其浓度变化率上限值，即标准变化率限值；步骤5、基于不同污染等级的标准变化率限值，实时动态识别突出值站点与对应的突出值时刻；步骤6、依照不同地区不同污染物的时间分布特征对研究时刻进行时段划分，再基于得到的历史突出值站点及其出现突出值的时刻，对突出值站点和划分时段的突出值进行统计分析，得到重点管控站点以及重点管控时段。2.如权利要求1所述的基于大气污染浓度变化率的突出值站点识别方法，其特征在于，在所述步骤1中，不同污染物研究时段内每t时刻的站点浓度平均值的计算公式为：式中，为t时刻的站点浓度平均值；为i站点t时刻的污染物浓度；t为时刻，其时间尺度包括分钟、小时、日和月中的一种；I为站点个数。3.如权利要求1所述的基于大气污染浓度变化率的突出值站点识别方法，其特征在于，所述步骤2，具体包括：基于逐时刻站点浓度平均值，依据各污染物浓度标准，将污染物的污染等级划分为多个等级，并得到逐时刻污染等级，统称为l
t
。4.如权利要求1所述的基于大气污染浓度变化率的突出值站点识别方法，其特征在于，所述步骤3，具体包括：步骤31、计算各站点污染物每t时刻相较于t
‑
1时刻的浓度变化率：式中，为i站点t时刻的污染物浓度变化率；t为时刻，其时间尺度包括分钟、小时、日和月中的一种；、分别为t时刻、t
‑
1时刻的污染物浓度；
步骤32、将逐时...

【专利技术属性】
技术研发人员：郎建垒，杨欣欣，周颖，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：