CN117875576B - 一种基于结构化案例库的城市大气污染分析方法 - Google Patents
一种基于结构化案例库的城市大气污染分析方法 Download PDFInfo
- Publication number
- CN117875576B CN117875576B CN202410283575.2A CN202410283575A CN117875576B CN 117875576 B CN117875576 B CN 117875576B CN 202410283575 A CN202410283575 A CN 202410283575A CN 117875576 B CN117875576 B CN 117875576B
- Authority
- CN
- China
- Prior art keywords
- data
- pollution source
- feature
- pollution
- weather
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 42
- 238000010586 diagram Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000000007 visual effect Effects 0.000 claims abstract description 12
- 238000012800 visualization Methods 0.000 claims abstract description 11
- 239000002893 slag Substances 0.000 claims description 11
- 238000009792 diffusion process Methods 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 8
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000004140 cleaning Methods 0.000 claims description 6
- 239000013618 particulate matter Substances 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 6
- 238000010801 machine learning Methods 0.000 claims description 5
- 230000003993 interaction Effects 0.000 claims description 4
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 claims description 3
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 3
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 claims description 3
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 claims description 3
- 238000003915 air pollution Methods 0.000 abstract description 3
- 239000000126 substance Substances 0.000 description 7
- 238000004590 computer program Methods 0.000 description 5
- 239000003344 environmental pollutant Substances 0.000 description 5
- 231100000719 pollutant Toxicity 0.000 description 5
- 238000011109 contamination Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012731 temporal analysis Methods 0.000 description 2
- 238000000700 time series analysis Methods 0.000 description 2
- 239000000443 aerosol Substances 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 229910052729 chemical element Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013079 data visualisation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012764 semi-quantitative analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于结构化案例库的城市大气污染分析方法,涉及城市大气污染分析技术领域。该方法包括获取污染源数据和气象数据,并进行数据结构化处理;将污染源结构化数据和气象结构化数据进行时序特征可视化处理;将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果;对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系。本发明能够有效的分析出污染源对空气污染指数的贡献度,并且可精确到单独的日期并量化分析当日特征的贡献度及重要性。
Description
技术领域
本发明涉及城市大气污染分析技术领域,具体涉及一种基于结构化案例库的城市大气污染分析方法。
背景技术
城市化的发展过程中伴随着三种主要污染源的产生,分别是交通源、工业源与扬尘源,使得城市大气受到不同程度的污染,给市民的工作和生活带来严重的不便。因此寻找大气污染的原因,采取相应的措施进行治理和预防,不仅能够改善城市的空气质量,降低大气的污染率,还能够保护人们的生存环境以及生命健康。大气污染的源解析技术主要包括源排放清单、扩散模型和受体模型等三类方法。
目前排放清单的定制方法有两种:一种是自上而下的方法,该方法通过能源统计年鉴中能源消耗、人口等信息,粗略估计各污染区排放量。该方法容易实施但精度不高;第二种是自下而上的方法,该方法通过大量实地测试获取,非常耗费人力物力,但排放数据的计算结果相对可靠。扩散模型主要用于处理非稳态下污染物的长距离输送,常用来模拟空气中污染物的排放、扩散、传输、物理化学反应和清除机制等过程。但该类方法的模型种类繁多,很难确定合适的空气质量模型,难以确定合适的模拟时空尺度,数据源的获取难度、技术复杂程度和软硬件平台都限制了该类模型的效果。同时,源排放清单的准确性直接影响了扩散模型模拟结果的准确性,采用不同的源排放清单数据或采用不同的数据处理方法都会导致模型的模拟结果差异。受体模型通过分析颗粒物中化学成分和物理特性来推断污染物来源,实现方法主要有显微镜法和化学法。显微镜法要建立庞大的源数据库,适用于定性或半定量分析形态特征明显的气溶胶定量分析需具备大量单个粒子结果才能代表整个样品。化学法主要有化学元素平衡法(CEB)、因子分析法(FA)、富集因子法、多元分析法、空间系列分析和时间系列分析等。化学质量平衡法很难满足排放源物质成分线性独立的要求,需要测量潜在污染源的化学成分谱,这需要全面的成分谱和准确的机物降解因子,需要大量的人力和物力。富集因子法只能定性给出研究元素污染程度有较大的局限性。时间序列分析法用于间接推定污染源不适合特定的污染源;空间模式分析法利用样品的空间分布与已知污染源的位置和化学组成之间的关系来确定污染源的影响,可间接的推定污染源但不能分析出污染源的有关资料。
发明内容
针对现有技术中的上述不足,本发明提供了一种基于结构化案例库的城市大气污染分析方法。
为了达到上述发明目的,本发明采用的技术方案为:
一种基于结构化案例库的城市大气污染分析方法,包括以下步骤:
获取污染源数据和气象数据,并进行数据结构化处理;
将污染源结构化数据和气象结构化数据进行时序特征可视化处理;
将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果;
对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系。
可选地,所述污染源数据具体包括:
空气质量数据、运渣车污染源数据、交通拥堵指数数据、交通流量数据、工厂用电数据和气象扩散条件数据。
可选地,所述空气质量数据具体包括:
细颗粒物、可吸入颗粒物、一氧化碳、二氧化氮、二氧化硫和臭氧。
可选地,所述气象数据具体包括:
温度、大气压强、湿度、风速和太阳辐射强度。
可选地,所述将污染源结构化数据和气象结构化数据进行时序特征可视化处理具体包括以下步骤:
对污染源结构化数据和气象结构化数据进行数据清洗和数据填充处理;
根据处理后的污染源结构化数据和气象结构化数据按照时序生成柱状图,得到时序特征可视化结果。
可选地,所述将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果,具体包括以下步骤:
将污染源结构化数据和气象结构化数据按同一时间尺度划分数据样本;
采用皮尔逊相关系数计算划分的数据样本中污染源特征与气象特征的相关系数;
根据污染源特征与气象特征的相关系数生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果。
可选地,所述对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系,具体包括以下步骤:
根据污染源结构化数据和气象结构化数据,采用机器学习方法建立污染源特征和气象特征与空气质量数据的映射关系;
根据污染源特征和气象特征与空气质量数据的映射关系,确定污染源特征和气象特征对大气污染的影响关系。
可选地,所述根据污染源特征和气象特征与空气质量数据的映射关系,确定污染源特征和气象特征对大气污染的影响关系具体包括以下步骤:
将污染源结构化数据中的污染源特征与气象结构化数据中的气象特征进行任意组合,建立特征联盟;
按照各个特征数据加入特征联盟的顺序,计算新加入特征联盟的特征数据对特征联盟的边际贡献度;
根据各个特征数据对特征联盟的边际贡献度将空气质量数据分配给各个特征数据,确定各个特征数据对空气质量数据的贡献度以及各种特征组合数据对空气质量数据的组合贡献度。
本发明具有以下有益效果:
本发明首先将城市污染源数据和气象数据结构化,并将污染源数据与气象数据做时序特征可视化展示;之后利用相关性矩阵图可视化污染源特征和气象特征的分布情况;最后进行特征重要性分析,能够有效的分析出污染源对空气污染指数的贡献度,并且可精确到单独的日期并量化分析当日特征的贡献度及重要性,从而为政府与企业的决策与污染源管控提供科学的、合理的依据,实现对环境污染源做到更精细和更有效的管控策略。
附图说明
图1为本发明中一种基于结构化案例库的城市大气污染分析方法的流程示意图;
图2为本发明中相关性矩阵图的示意图;
图3为本发明中个体样本中特征贡献图的示意图;
图4为本发明中全量样本中特征贡献图的示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,本发明实施例提供了一种基于结构化案例库的城市大气污染分析方法,包括以下步骤S1至S4:
S1、获取污染源数据和气象数据,并进行数据结构化处理;
在本发明的一个可选实施例中,本实施例通过步骤S1获取的分析数据包括污染源数据和气象数据。其中污染源数据具体包括:空气质量数据、运渣车污染源数据、交通拥堵指数数据、交通流量数据、工厂用电数据和气象扩散条件数据。
本实施例获取的空气质量数据具体包括:细颗粒物、可吸入颗粒物、一氧化碳、二氧化氮、二氧化硫和臭氧。
本实施例获取运渣车污染源数据的方法为:
首先获取全量的运渣车GPS轨迹点,挖掘运渣车的停驻点,然后根据停驻点信息判断活跃的运渣车,以此为依据对运渣车的车牌信息去重后得到活跃的运渣车数量。
本实施例获取交通拥堵指数数据的方法为:
首先计算早晚高峰的交通拥堵指数,然后计算早晚高峰期间的道路平均拥堵指数,计算道路拥堵指数的90%分位数,读取大于该拥堵指数的道路然后计算城市范围全域的平均道路拥堵指数。
本实施例获取交通流量数据的方法为:
根据城市范围内的交通卡口拍摄记录,根据每小时内的车辆记录信息对车牌号码进行去重得到每小时内的交通流量,对24小时内的交通流量求和得到以天为单位的交通流量数据。
本实施例获取工厂用电数据的方法为:
对工厂用电记录进行数据清洗,删除异常值,然后计算工厂用电的每天的用电量,得到工厂用电数据。
本实施例获取的气象数据具体包括:温度、大气压强、湿度、风速和太阳辐射强度。获取方法为对温度、大气压强、湿度、风速数据计算24小时平均值,对太阳辐射强度计算8小时的滑动平均值,选8小时滑动平均值最大值作为当天的平均值,得到各种气象数据。
本发明获取污染源数据和气象数据及空气质量数据后,剔除数据中的重复值和异常值;然后对数据的缺失值进行线性插补,将所获得的数据按照每日的时间粒度进行聚合得到结构化的数据表,数据样本量大小等于数据收集结束日期减去数据收集开始日期。
S2、将污染源结构化数据和气象结构化数据以柱状图的形式进行时序特征可视化处理;
在本发明的一个可选实施例中,本实施例通过步骤S2将污染源结构化数据和气象结构化数据进行时序特征可视化处理具体包括以下步骤S21至S22:
S21、对污染源结构化数据和气象结构化数据进行数据清洗和数据填充处理;
S22、根据处理后的污染源结构化数据和气象结构化数据按照时序生成柱状图,得到时序特征可视化结果。
具体而言,本实施例先将污染源结构化数据中的运渣车污染源数据、交通拥堵指数数据、交通流量数据、工厂用电数据和气象扩散条件数据以及气象结构化数据中的温度、大气压强、湿度、风速和太阳辐射强度进行数据清洗;然后将运渣车污染源数据、交通拥堵指数数据、交通流量数据、工厂用电数据、温度、大气压强、湿度、风速和太阳辐射强度中的缺失值采用线性插值方法进行填充,将气象扩散条件数据中的缺失值采用均值方法进行填充;最后根据处理后的污染源结构化数据和气象结构化数据按照时序生成柱状图,得到时序特征可视化结果。
S3、将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果;
在本发明的一个可选实施例中,本实施例通过步骤S3将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果,具体包括以下步骤S31至S33:
S31、将污染源结构化数据和气象结构化数据按同一时间尺度划分数据样本;
S32、采用皮尔逊相关系数计算划分的数据样本中污染源特征与气象特征的相关系数;
S33、根据污染源特征与气象特征的相关系数生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果。
具体而言,本实施例经过数据可视化及数据清洗后,先将数据的时间尺度从小时级别缩放到日期级别,并将以日期为时间单位的尺度划分数据样本。
然后采用皮尔逊相关系数计算划分的数据样本中污染源特征与气象特征的相关系数;
最后根据污染源特征与气象特征的相关系数生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果,如图2所示。
S4、对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系。
在本发明的一个可选实施例中,本实施例通过步骤S4对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系,具体包括以下步骤S41至S42:
S41、根据污染源结构化数据和气象结构化数据,采用机器学习方法建立污染源特征和气象特征与空气质量数据的映射关系,所建立的机器学习模型为XGboost树模型f;
S42、根据污染源特征和气象特征与空气质量数据的映射关系,确定污染源特征和气象特征对大气污染的影响关系,具体包括以下步骤:
S421、将污染源结构化数据中的污染源与气象结构化数据中的气象特征进行任意组合,建立个特征组合/>,其中/>表示污染源特征和气象类特征的维度大小;
S422、计算特征的shapley值,并从全量数据样本矩阵X中抽取部分数据样本/>。按照各个特征数据加入特征联盟的顺序,计算新加入特征联盟的特征数据对特征联盟的边际贡献度;根据各个特征数据对特征联盟的边际贡献度将空气质量数据分配给各个特征数据,确定各个特征数据对空气质量数据的贡献度以及各种特征组合数据对空气质量数据的组合贡献度。其中某个特征/>的Shapley值的计算公式如下所示:
其中,表示特征向量集合,/>表示特征向量的维度,/>表示计算得到的特征Shapley值,!表示阶乘。
S423、固定特征子集S中特征的值和样本的特征值,除子集/>和第/>个特征值不变之外其他特征值使用平均值z代替,构造样本/>,有特征/>的情况:。
S424、固定特征子集中特征的值,除子集/>中的特征值保持不变外其他特征值使用平均值z代替,构造样本/>,没有特征/>的情况:。
S425、计算Shapley值的平均值。
具体而言,本实施例为了定量分析城市污染源类型贡献度,首先根据污染源结构化数据和气象结构化数据,采用机器学习方法建立污染源特征和气象特征与空气质量数据的映射关系;
然后根据污染源特征和气象特征与空气质量数据的映射关系,采用合作博弈论方法计算各个特征对空气质量数据的贡献度,从而确定污染源特征和气象特征对大气污染的影响关系。
根据以上计算过程可得到个体样本中特征shapley值,如图3所示,以及全量样本的shapley值的均值,如图4所示。
计算两个特征交互作用下的Shapley值,可用于探索分析两种特征i,j之间的交互作用,其计算公式如下所示:
本发明能够单独分析个例样本中各污染源浓度的贡献度,也能泛化推广至一定时间范围内的城市污染物特征重要性;通过特征重要性分析,既能分析单个污染源特征的贡献值,也能分析多个污染特征组合的特征贡献值,同时能够在结合气象条件对污染物排放进行定性定量的分析。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (5)
1.一种基于结构化案例库的城市大气污染分析方法,其特征在于,包括以下步骤:
获取污染源数据和气象数据,并进行数据结构化处理;所述污染源数据具体包括:空气质量数据、运渣车污染源数据、交通拥堵指数数据、交通流量数据、工厂用电数据;
将污染源结构化数据和气象结构化数据进行时序特征可视化处理;
将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果;
对污染源结构化数据和气象结构化数据进行特征重要性分析,确定污染源特征和气象特征对大气污染的影响关系,具体包括以下步骤:
根据污染源结构化数据和气象结构化数据,采用机器学习方法建立污染源特征和气象特征与空气质量数据的映射关系;
根据污染源特征和气象特征与空气质量数据的映射关系,确定污染源特征和气象特征对大气污染的影响关系,具体包括以下步骤:
将污染源结构化数据中的污染源与气象结构化数据中的气象特征进行任意组合,建立个特征组合/>,其中/>表示污染源特征和气象类特征的维度大小;
计算特征的shapley值,并从全量数据样本矩阵X中抽取部分数据样本/>;按照各个特征数据加入特征联盟的顺序,计算新加入特征联盟的特征数据对特征联盟的边际贡献度;根据各个特征数据对特征联盟的边际贡献度将空气质量数据分配给各个特征数据,确定各个特征数据对空气质量数据的贡献度以及各种特征组合数据对空气质量数据的组合贡献度;其中某个特征/>的Shapley值的计算公式如下所示:
其中,表示特征向量集合,/>表示特征向量的维度,/>表示计算得到的特征Shapley值,!表示阶乘;
固定特征子集S中特征的值和样本的特征值,除子集/>和第/>个特征值不变之外其他特征值使用平均值z代替,构造样本/>,有特征/>的情况:;
固定特征子集中特征的值,除子集/>中的特征值保持不变外其他特征值使用平均值z代替,构造样本/>,没有特征/>的情况:/>;
计算Shapley值的平均值;
计算两个特征交互作用下的Shapley值,用于探索分析两种特征i,j之间的交互作用,其计算公式如下所示:
。
2.根据权利要求1所述的一种基于结构化案例库的城市大气污染分析方法,其特征在于,所述空气质量数据具体包括:
细颗粒物、可吸入颗粒物、一氧化碳、二氧化氮、二氧化硫和臭氧。
3.根据权利要求1所述的一种基于结构化案例库的城市大气污染分析方法,其特征在于,所述气象数据具体包括:
温度、大气压强、湿度、风速、太阳辐射强度和气象扩散条件数据。
4.根据权利要求1所述的一种基于结构化案例库的城市大气污染分析方法,其特征在于,所述将污染源结构化数据和气象结构化数据进行时序特征可视化处理具体包括以下步骤:
对污染源结构化数据和气象结构化数据进行数据清洗和数据填充处理;
根据处理后的污染源结构化数据和气象结构化数据按照时序生成柱状图,得到时序特征可视化结果。
5.根据权利要求1所述的一种基于结构化案例库的城市大气污染分析方法,其特征在于,所述将污染源结构化数据和气象结构化数据按同一时间尺度生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果,具体包括以下步骤:
将污染源结构化数据和气象结构化数据按同一时间尺度划分数据样本;
采用皮尔逊相关系数计算划分的数据样本中污染源特征与气象特征的相关系数;
根据污染源特征与气象特征的相关系数生成相关性矩阵图,得到污染源特征和气象特征的可视化分布结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410283575.2A CN117875576B (zh) | 2024-03-13 | 2024-03-13 | 一种基于结构化案例库的城市大气污染分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410283575.2A CN117875576B (zh) | 2024-03-13 | 2024-03-13 | 一种基于结构化案例库的城市大气污染分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117875576A CN117875576A (zh) | 2024-04-12 |
CN117875576B true CN117875576B (zh) | 2024-05-24 |
Family
ID=90590428
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410283575.2A Active CN117875576B (zh) | 2024-03-13 | 2024-03-13 | 一种基于结构化案例库的城市大气污染分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117875576B (zh) |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2012149045A (ru) * | 2012-11-19 | 2014-05-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Государственный университет по землеустройству" | Способ оперативного контроля атмосферных загрязнений локальных территорий |
CN104881546A (zh) * | 2015-06-01 | 2015-09-02 | 中国科学院上海高等研究院 | 一种提高大气污染模型预测效率的方法 |
KR20170079648A (ko) * | 2015-12-30 | 2017-07-10 | 대한민국(국민안전처 국립재난안전연구원장) | 미래위험 변화예측 분석 시스템 |
CN110427533A (zh) * | 2019-07-25 | 2019-11-08 | 东北师范大学 | 基于时序粒子跟踪的污染传播模式可视分析方法及*** |
CN111538957A (zh) * | 2020-04-21 | 2020-08-14 | 中科三清科技有限公司 | 大气污染物来源贡献度的获取方法、装置、设备及介质 |
CN112418609A (zh) * | 2020-10-30 | 2021-02-26 | 暨南大学 | 一种基于面-网格-点的大气二次污染精准溯源方法 |
CN112613675A (zh) * | 2020-12-29 | 2021-04-06 | 南开大学 | 一种分析污染源与气象因素对不同程度pm2.5污染影响贡献和效应的机器学习模型 |
WO2021208393A1 (zh) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | 一种大气污染物排放清单的反演估算方法 |
CN113804829A (zh) * | 2021-08-20 | 2021-12-17 | 重庆市生态环境监测中心 | 一种大气污染天空地一体化实时监测***及方法 |
CN114359002A (zh) * | 2022-03-21 | 2022-04-15 | 四川国蓝中天环境科技集团有限公司 | 基于移动监测趋势挖掘的大气污染小尺度溯源方法及*** |
CN114356880A (zh) * | 2021-12-15 | 2022-04-15 | 浙江工业大学 | 基于数据驱动的小尺度区域大气污染物快速溯源方法 |
CN115293360A (zh) * | 2022-07-12 | 2022-11-04 | 中国环境科学研究院 | 一种大气污染物浓度分析方法、装置及存储介质 |
CN115542429A (zh) * | 2022-09-20 | 2022-12-30 | 生态环境部环境工程评估中心 | 一种基于XGBoost的臭氧质量预测方法及*** |
CN115757534A (zh) * | 2022-10-20 | 2023-03-07 | 电子科技大学 | 基于潜在源贡献分析的深度学习模型的空气质量预测方法 |
CN116128421A (zh) * | 2022-09-23 | 2023-05-16 | 北京清创美科环境科技有限公司 | 一种耦合污染源排放和解析的大气污染管控方案生成方法 |
CN116662911A (zh) * | 2023-06-27 | 2023-08-29 | 安徽中科蓝壹信息科技有限公司 | 污染源和气象因素对VOCs初始浓度影响贡献评估方法 |
CN116821633A (zh) * | 2023-08-01 | 2023-09-29 | 深圳博沃智慧科技有限公司 | 大气污染成因分析方法、装置、设备及存储介质 |
KR20230144301A (ko) * | 2022-04-07 | 2023-10-16 | 전북대학교산학협력단 | 환경보건 모니터링 시스템 및 그 방법 |
CN117171546A (zh) * | 2023-09-27 | 2023-12-05 | 安徽熵沃智能科技有限公司 | 一种结合污染扩散指数的城市空气质量预测方法 |
CN117332815A (zh) * | 2023-09-22 | 2024-01-02 | 温州大学 | 一种工业园区大气污染的预测方法及预测预警*** |
CN117495637A (zh) * | 2023-11-17 | 2024-02-02 | 华北科技学院(中国煤矿安全技术培训中心) | 基于气象因素和污染传输的大气污染物评估方法及*** |
CN117708617A (zh) * | 2024-02-05 | 2024-03-15 | 四川国蓝中天环境科技集团有限公司 | 基于多源大数据及污染特征时空匹配的大气污染溯源方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20220068053A1 (en) * | 2020-08-25 | 2022-03-03 | ANI Technologies Private Limited | Determination of health status of vehicular systems in vehicles |
CN115293567A (zh) * | 2022-08-03 | 2022-11-04 | 四川省生态环境科学研究院 | 一种大气污染动态管控***及方法 |
-
2024
- 2024-03-13 CN CN202410283575.2A patent/CN117875576B/zh active Active
Patent Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2012149045A (ru) * | 2012-11-19 | 2014-05-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Государственный университет по землеустройству" | Способ оперативного контроля атмосферных загрязнений локальных территорий |
CN104881546A (zh) * | 2015-06-01 | 2015-09-02 | 中国科学院上海高等研究院 | 一种提高大气污染模型预测效率的方法 |
KR20170079648A (ko) * | 2015-12-30 | 2017-07-10 | 대한민국(국민안전처 국립재난안전연구원장) | 미래위험 변화예측 분석 시스템 |
CN110427533A (zh) * | 2019-07-25 | 2019-11-08 | 东北师范大学 | 基于时序粒子跟踪的污染传播模式可视分析方法及*** |
WO2021208393A1 (zh) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | 一种大气污染物排放清单的反演估算方法 |
CN111538957A (zh) * | 2020-04-21 | 2020-08-14 | 中科三清科技有限公司 | 大气污染物来源贡献度的获取方法、装置、设备及介质 |
CN112418609A (zh) * | 2020-10-30 | 2021-02-26 | 暨南大学 | 一种基于面-网格-点的大气二次污染精准溯源方法 |
CN112613675A (zh) * | 2020-12-29 | 2021-04-06 | 南开大学 | 一种分析污染源与气象因素对不同程度pm2.5污染影响贡献和效应的机器学习模型 |
CN113804829A (zh) * | 2021-08-20 | 2021-12-17 | 重庆市生态环境监测中心 | 一种大气污染天空地一体化实时监测***及方法 |
CN114356880A (zh) * | 2021-12-15 | 2022-04-15 | 浙江工业大学 | 基于数据驱动的小尺度区域大气污染物快速溯源方法 |
CN114359002A (zh) * | 2022-03-21 | 2022-04-15 | 四川国蓝中天环境科技集团有限公司 | 基于移动监测趋势挖掘的大气污染小尺度溯源方法及*** |
KR20230144301A (ko) * | 2022-04-07 | 2023-10-16 | 전북대학교산학협력단 | 환경보건 모니터링 시스템 및 그 방법 |
CN115293360A (zh) * | 2022-07-12 | 2022-11-04 | 中国环境科学研究院 | 一种大气污染物浓度分析方法、装置及存储介质 |
CN115542429A (zh) * | 2022-09-20 | 2022-12-30 | 生态环境部环境工程评估中心 | 一种基于XGBoost的臭氧质量预测方法及*** |
CN116128421A (zh) * | 2022-09-23 | 2023-05-16 | 北京清创美科环境科技有限公司 | 一种耦合污染源排放和解析的大气污染管控方案生成方法 |
CN115757534A (zh) * | 2022-10-20 | 2023-03-07 | 电子科技大学 | 基于潜在源贡献分析的深度学习模型的空气质量预测方法 |
CN116662911A (zh) * | 2023-06-27 | 2023-08-29 | 安徽中科蓝壹信息科技有限公司 | 污染源和气象因素对VOCs初始浓度影响贡献评估方法 |
CN116821633A (zh) * | 2023-08-01 | 2023-09-29 | 深圳博沃智慧科技有限公司 | 大气污染成因分析方法、装置、设备及存储介质 |
CN117332815A (zh) * | 2023-09-22 | 2024-01-02 | 温州大学 | 一种工业园区大气污染的预测方法及预测预警*** |
CN117171546A (zh) * | 2023-09-27 | 2023-12-05 | 安徽熵沃智能科技有限公司 | 一种结合污染扩散指数的城市空气质量预测方法 |
CN117495637A (zh) * | 2023-11-17 | 2024-02-02 | 华北科技学院(中国煤矿安全技术培训中心) | 基于气象因素和污染传输的大气污染物评估方法及*** |
CN117708617A (zh) * | 2024-02-05 | 2024-03-15 | 四川国蓝中天环境科技集团有限公司 | 基于多源大数据及污染特征时空匹配的大气污染溯源方法 |
Non-Patent Citations (8)
Title |
---|
High-resolution mapping of regional VOCs using the enhanced space-time extreme gradient boosting machine (XGBoost) in Shanghai;Bingqing Lu et al.;Science of the Total Environment;20230913;全文 * |
Machine learning model to predict vehicle electrification impacts on urban air quality and related human health effects;V. Calatayud et al.;Environmental Research;20230404;全文 * |
Time-sensitivepredictionofNO2concentrationin Chinausinganensemblemachinelearningmodel frommulti-source data;ChenliangTao et al.;Journal of Environmental Sciences;20230226;全文 * |
北京城市大气环境污染机理与调控原理;徐祥德;丁国安;卞林根;;应用气象学报;20061225(第06期);全文 * |
城市大气污染的来源特征分析;阎静;;资源节约与环保;20151025(第10期);全文 * |
基于多源数据的城市大气污染物传输溯源方法研究;李皓维;CNKI硕士电子期刊;20230601;全文 * |
新排放标准下的硫回收装置技术改造方案;韩科;郭勇;;煤炭加工与综合利用;20170425(第04期);全文 * |
泰州市PM2.5与常规空气污染物及气象因素的相关性分析;王莉;贺玉川;于光;黄灵;;中国卫生检验杂志;20181125(第22期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117875576A (zh) | 2024-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110824110B (zh) | 一种基于拉格朗日轨迹模式和化学盒子模式的区域臭氧污染溯源*** | |
CN110489785B (zh) | 一种大气污染物在线源解析方法及*** | |
CN111368401A (zh) | 污染源的溯源方法、装置和存储介质 | |
Qiu et al. | Atmospheric dispersion prediction and source estimation of hazardous gas using artificial neural network, particle swarm optimization and expectation maximization | |
Xing et al. | Quantifying the emission changes and associated air quality impacts during the COVID-19 pandemic on the North China Plain: a response modeling study | |
CN111444633B (zh) | 一种针对大气污染过程的定量解析方法和*** | |
Qiao et al. | Modeling dry and wet deposition of sulfate, nitrate, and ammonium ions in Jiuzhaigou National Nature Reserve, China using a source-oriented CMAQ model: Part I. Base case model results | |
Zhou et al. | A new statistical approach for establishing high-resolution emission inventory of primary gaseous air pollutants | |
CN111798928B (zh) | 大气颗粒物污染源解析方法和装置 | |
Kormi et al. | Estimation of landfill methane emissions using stochastic search methods | |
Foley et al. | Dynamic evaluation of CMAQ part I: Separating the effects of changing emissions and changing meteorology on ozone levels between 2002 and 2005 in the eastern US | |
CN114661849B (zh) | 污染溯源方法及装置 | |
Cai et al. | Impacts of winter heating on the atmospheric pollution of northern China’s prefectural cities: evidence from a regression discontinuity design | |
CN110348746B (zh) | 基于单个污染源的空气质量影响评估方法和装置 | |
CN111753426B (zh) | 颗粒物污染的源解析方法和装置 | |
CN105512485B (zh) | 一种估算细颗粒物及其前体物环境容量的新方法 | |
Ying et al. | Modeling air quality during the California Regional PM10/PM2. 5 Air Quality Study (CRPAQS) using the UCD/CIT source-oriented air quality model–Part I. Base case model results | |
CN114858976A (zh) | 一种工业园区大气质量智能分析方法及*** | |
CN112967764B (zh) | 多技术耦合的污染物源解析方法、装置 | |
Volta et al. | GAMES, a comprehensive gas aerosol modelling evaluation system | |
CN112991132A (zh) | Pm2.5和臭氧协同管控区域识别方法及装置 | |
Barbulescu et al. | Modeling the carbon monoxide dissipation in Timisoara, Romania | |
CN115629159A (zh) | 一种基于多源数据的臭氧及前体物溯源方法及装置 | |
Xing et al. | Mimicking atmospheric photochemical modeling with a deep neural network | |
Hassani et al. | Citizen-operated mobile low-cost sensors for urban PM2. 5 monitoring: field calibration, uncertainty estimation, and application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |