CN116128260A - 基于数据样本的重点企业环境分析方法 - Google Patents
基于数据样本的重点企业环境分析方法 Download PDFInfo
- Publication number
- CN116128260A CN116128260A CN202310409806.5A CN202310409806A CN116128260A CN 116128260 A CN116128260 A CN 116128260A CN 202310409806 A CN202310409806 A CN 202310409806A CN 116128260 A CN116128260 A CN 116128260A
- Authority
- CN
- China
- Prior art keywords
- monitoring point
- data
- target monitoring
- moment
- degree
- 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.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 20
- 238000012544 monitoring process Methods 0.000 claims abstract description 307
- 230000005856 abnormality Effects 0.000 claims abstract description 106
- 230000002159 abnormal effect Effects 0.000 claims abstract description 41
- 238000012216 screening Methods 0.000 claims abstract description 9
- 238000000034 method Methods 0.000 claims description 36
- 230000008859 change Effects 0.000 claims description 32
- 238000009826 distribution Methods 0.000 claims description 5
- 238000003891 environmental analysis Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract 1
- 239000003344 environmental pollutant Substances 0.000 description 30
- 231100000719 pollutant Toxicity 0.000 description 30
- 230000005540 biological transmission Effects 0.000 description 14
- 239000007789 gas Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 7
- 239000000356 contaminant Substances 0.000 description 6
- 230000007613 environmental effect Effects 0.000 description 6
- 230000035772 mutation Effects 0.000 description 5
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 4
- 239000002912 waste gas Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- 101150096839 Fcmr gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 229910002091 carbon monoxide Inorganic materials 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 150000002222 fluorine compounds Chemical class 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Images
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
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06311—Scheduling, planning or task assignment for a person or group
- G06Q10/063114—Status monitoring or status determination for a person or group
-
- 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
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
- G06Q30/00—Commerce
- G06Q30/018—Certifying business or products
-
- 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)
- Economics (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Accounting & Taxation (AREA)
- Finance (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明涉及数据处理领域,具体涉及基于数据样本的重点企业环境分析方法,包括:获得目标监测点的异常程度,进而获得目标监测点的第二异常程度。根据参考监测点的第二异常程度、以及目标监测点与参考监测点的监测数据的差异获得目标监测点的第三异常程度。根据目标监测点的第三异常程度、第一数量获得目标监测点的第四异常程度。根据目标监测点的第二异常程度和第四异常程度获得目标监测点的最终异常程度,根据目标监测点的最终异常程度进行监测数据筛选和环境分析。本发明用数据分析方式对采集的监测数据进行异常分析处理,得到更加接近实际的样本数据,使得对污染气体的超标检测更加准确。
Description
技术领域
本发明涉及数据处理领域,具体涉及基于数据样本的重点企业环境分析方法。
背景技术
随着科学技术的发展,各种工业企业得到迅速的发展,同时各种工厂生产过程对环境的影响显而易见,其中包括化工工厂的污染物排放超标事件时有发生,所排放的各种废水、废气破坏着生态环境的平衡,甚至对周围居民的健康产生加大的威胁,所以在实际工厂需要进行污染物的监测,其中由于废气传播速度快,可控性差,所以需要重点监测。
在进行污染物监测时,为了提高检测准确性,需要在采集的数据中筛选出能够反应污染物实际数据的数据样本,一般影响监测数据准确性的因素主要为监测点数据采集异常,主要表现为数据突变,现有操作一般直接利用数据分布差异进行异常识别,但是在实际监测过程中,风向影响着污染物的传播,表现为监测点数据与风向的相关性,并且在风向变化时同样会引起监测点数据的变化,使得数据分布差异不能直接反应数据的异常,即所进行的数据异常分析不准确,从而影响数据样本质量。
发明内容
本发明提供基于数据样本的重点企业环境分析方法,以解决现有的问题。
本发明的基于数据样本的重点企业环境分析方法采用如下技术方案:
本发明一个实施例提供了基于数据样本的重点企业环境分析方法,该方法包括以下步骤:
获取每个时刻每个监测点污染气体浓度的监测数据;将每个时刻下的任意一个监测点记为目标监测点;
根据目标监测点在相邻时刻的监测数据之间的差异获得目标监测点的异常程度;
根据目标监测点的异常程度、目标监测点与相邻时刻风向与风速的变化值获得每个目标监测点的第二异常程度;
将与目标监测点相对于工厂方向与风向之间的夹角最接近的监测点记为目标监测点的参考监测点;根据每个时刻下在参考监测点的第二异常程度、目标监测点相对于工厂方向与风向夹角与参考监测点相对于工厂方向与风向夹角的差异、目标检测点相对于工厂方向与参考监测点相对于工厂方向的夹角以及目标监测点与参考监测点的监测数据的差异获得目标监测点的第三异常程度;
获取每个时刻下所有监测点相对于工厂方向与风向夹角,记为夹角集合S;获取S中与目标监测点相对于工厂方向与风向夹角差异小于预设阈值的夹角所对应的监测点,记为集合S1;获取每个时刻的相邻时刻下所有监测点相对于工厂方向与风向,记为夹角集合S2;获取S2中与目标监测点相对于工厂方向与风向夹角差异小于预设阈值的夹角所对应的监测点,记为集合S3;将S3与S1中相同监测点的数量记为目标监测点在相邻时刻的第一数量;
根据目标监测点的第三异常程度、每个时刻与相邻时刻的时间间隔、目标监测点与相邻时刻的风向变化角度、每个时刻的风速以及第一数量获得目标监测点的第四异常程度;
根据目标监测点的第二异常程度和第四异常程度获得目标监测点的最终异常程度,根据目标监测点的最终异常程度进行监测数据筛选和环境分析。
进一步地,所述目标监测点的异常程度的获取方法包括:
将每个时刻在预设时间窗口内的所有监测数据作为每个时刻的周围最相邻数据,目标监测点在每个时刻的监测数据与每个周围最相邻数据分别计算差异,记为第一差异,所有第一差异的均值记为目标监测点的异常程度。
进一步地,所述目标监测点的第二异常程度的获取方法为:
获取每个时刻工厂到监测点的方向与风向之间的夹角和风速大小;得到每个时刻与相邻时刻风向与风速的变化值,将用风向和风速变化值相乘之后取倒数,求取每个时刻的所有相邻时刻下所述倒数的第一均值;最后将第一均值与风速大小、风向与监测点到工厂方向夹角的倒数以及目标监测点的异常程度的乘积记为目标监测点的第二异常程度。
进一步地,所述目标监测点的第三异常程度的获取方法为:
对于第i时刻的第j个监测点,即为目标监测点,目标监测点的第三异常程度公式为:
其中,表示目标监测点的监测数据,表示目标监测点相对于工厂方向与风向之间的夹角最接近的最接近的第u个监测点数据,即参考监测点;表示第i时刻第u个监测点相对于工厂方向与风向的夹角和目标监测点相对于工厂方向与风向夹角之间的差异;表示第i时刻第u个监测点相对于工厂方向与目标监测点相对于工厂方向之间的夹角;表示参考监测点的个数;表示第u个监测点数据在风向影响下的第二异常程度;为目标监测点的第三异常程度。
进一步地,所述目标监测点的第四异常程度的获取方法为:
目标监测点的第四异常程度公式为:
其中,表示对应第u个监测点在第i时刻与之前相邻的第v时刻风向的变化角度;表示第v时刻与第i时刻的时间间隔;表示第一数量;表示在第v时刻下的风速大小;表示参考监测点的个数;为预设值;表示为目标监测点的第三异常程度;表示目标监测点的第四异常程度。
进一步地,所述目标监测点的最终异常程度的获取方法为:
由目标监测点的第二异常程度和第四异常程度相乘得到目标监测点的最终异常程度。
进一步地,所述进行监测数据筛选方法包括为:
用正态分布的原则进行设置阈值,所述阈值设置为,其中,分别表示所有监测点最终异常程度的均值和标准差;我们将大于阈值的监测数据进行删除,将小于的监测数据保留,最后筛选得到用于环境分析的监测数据。
本发明的技术方案的有益效果是:
(1)根据风向对污染物传播的影响,反应单个监测点与多个监测点数据差异关系对数据异常表现程度,提高监测数据突出表现反应数据异常的准确性,即获得更加准确的数据异常程度。
(2)考虑风向变化引起的污染物传播变化对监测点数据的影响,避免污染物相对于监测点的传播变化所引起监测数据变化被误识别为异常数据,提高所获得异常程度的可靠性。
(3)通过根据监测点数据采集异常表现进行数据样本筛选,获得反应污染物实际分布的检测数据样本,提高环境监测分析的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于数据样本的重点企业环境分析方法的步骤流程图。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种基于数据样本的重点企业环境分析方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的基于数据样本的重点企业环境分析方法的具体方案。
请参阅图1,其示出了本发明一个实施例提供的基于数据样本的重点企业环境分析方法的步骤流程图,该方法包括以下步骤:
步骤S001:获取每个时刻每个监测点的监测数据。
在工厂废气监测过程中,因为废气的传播范围不可控,所以需要在工厂周围多个位置获取原始废气浓度数据。所以我们首先需要在工厂周围安装多个监测点,在此监测点数据我们采用多种污染气体检测仪。在气体污染中有多种成分,比如:二氧化硫、氮氧化物、氟化物、一氧化碳、二氧化碳等。在检测污染气体浓度时每种污染气体有对应的化学试剂进行检测,此实施例我们以二氧化硫浓度检测进行分析。用监测点进行对污染气体浓度数据的监测获取。且我们设定检测点之间的间距为50米。此时多个监测点同时进行监测,获得监测数据传输到监测数据处理中心,数据处理中心进行对原始数据进行分析。后续步骤通过对原始数据分析,确定数据在采集过程中的异常数据,即获得反应真实污染物的数据样本。后续步骤所述数据异常即为采集过程异常,即不具备样本分析的数据。
在后面分析中我们还需要考虑风的影响,所以此处我们还需安装风速传感器,用来获取风速的大小与风向的角度情况。
步骤S002:获得目标监测点的异常程度,进而获得目标监测点的第二异常程度。根据参考监测点的第二异常程度、以及目标监测点与参考监测点的监测数据的差异获得目标监测点的第三异常程度。根据目标监测点的第三异常程度、第一数量获得目标监测点的第四异常程度。根据目标监测点的第二异常程度和第四异常程度获得目标监测点的最终异常程度。
对于监测数据在数据采集过程中的异常主要表现为数据突变,即数据相对于相邻数据具有较大差异,数据差异性越大,其异常程度越大。风速与风向影响监测数据之间的差异关系,从而影响数据异常关系。所以本发明通过分析风向对污染物传播的影响,反应监测数据之间的差异对异常的表现关系,从而进行准确的异常分析。
需要说明的是,工厂大废气污染物在工厂内部产生向四周扩散,一般废气在工厂运行过程中连续产生,所以其传播过程具有一定持续性,此时对于监测点所采集的监测数据具有一定连续性。此时由于采集设备的***或者偶然误差导致监测异常数据主要表现为与相邻数据的突变,此时获得异常程度为:
其中,表示第i时刻的监测数据大小,表示周围最相邻的第t个数据,周围最相邻数据是指以第i时刻为中心的且长度为n的时间窗口内的所有数据,本实施例以n=7为例进行叙述。表示二者差异,其值越大,相对突变的可能性越大,表示与周围相邻的n个数据的平均差异,其值越大,相对于周围数据突变的可能性越大,即其局部异常程度越大。即表示第i个数据的局部异常程度。
实际大气污染物传播具有多变性,其传播方向以及污染物浓度容易受风向与风速的影响,并且随着风向以及风速的变化使得对污染物传播的影响发生变化,所以为了准确判断监测点数据的准确性,需要综合判断风向以及风速对其影响程度。
需要说明的是,一般污染物受风影响主要表现为污染物随着风向传播,污染物随着风向传播程度越大,污染物浓度越稳定,即监测数据变化差异越小。在风向发生变化时,污染物的传播方向发生对应的变化,对于监测点的污染物浓度监测数据即发生变化,此时对于同一监测点监测数据在时间序列上发生变化,即可能导致监测数据发生变化,从而影响监测数据异常识别。
监测数据异常受风向的影响主要表现为:工厂到监测点的方向与风向之间的夹角越小,监测点数据越稳定;风向在当前时间段变化越小,监测点数据越稳定,数据差异受传播影响的可能性越小,数据差异所表现的采集异常可信度越高。
进一步需要说明的是,在用监测器采集气体污染物浓度时会有误差,一个误差来自于采集数据时出现的误差,另一个误差来自于外界因素的干扰造成的误差。现在我们将工厂到监测点的方向与风向之间的夹角进行分析,我们知道,当夹角越小,外界因素风向对它的的影响较小,这时采集的样本更大的误差来自于采集误差,此时对于上述样本异常程度的可信度越高,而样本中出现采集数据突变的异常点时,则样本的气体的污染检测就越能检测出异常;而当夹角越大,外界因素风向对它的影响较大,这是采集的样本更大的误差来自于外界环境的干扰,而在采集时出现的误差占很少,此时对于上述样本异常程度的可信度就相对较低,样本检测出的异常可能更多不是来自于样本气体浓度的异常。
具体获得数据第二异常程度为:
其中,表示第i时刻工厂到监测点的方向与风向之间的夹角,其值越小,监测点数据越稳定,数据差异所表现的采集异常的可信度越高,表示第i时刻风速,其值越大,风向对的影响越大;,分别表示第i个时刻与之前相邻的第v时刻风向与风速的变化值,其值越大,风向导致监测数据变化程度越大,此时数据变化反应的异常的可信度越低,即所反映的异常程度越小;则表示第i时刻与之前相邻的个时刻风向与风速变化对的影响,其值越大,数据异常程度越大,其中设置为3。表示第i个时刻监测点数据受风向与风速影响下的异常程度,即第二异常程度;其中常量0.1是为了避免分母为零,本实施例对该常量的具体取值不进行具体限定,其他实施例中可以设置其他值以保证分母不为零。
至此,上述获得了第i时刻下任意监测点的异常程度。
上述分析获得单个监测点受风向影响下数据异常程度,此时在实际监测中,多个监测点与工厂的相对方向不同,导致风向对不同位置监测点数据影响程度不同,并且在风向的影响下不同位置检测点数据之间的关系反应污染物传播关系,不同监测点数据具有基于风向的相关性关系,所以不同位置检测数据之间的关系同样能够反应数据的异常。
需要说明的是,在同一时刻多个监测点同时受风向影响,此时因为不同监测点方向与风向的关系存在差异,所以不同监测点监测数据之间的相关关系具有明显的差异,具体表现为:监测点方向与风向夹角越一致,对应的监测数据越一致,所以监测点数据之间的差异越能够反映对应数据的异常程度,即第三异常程度。
此时,根据监测点相对于工厂方向与风向夹角之间的差异,反应不同数据的异常程度。需要说明的是,夹角之间的差异是指两个夹角之间差值的绝对值;第三异常程度具体计算方式为:
其中,表示第j个监测点第i时刻的监测数据,表示第i时刻的第j个监测点相对于工厂方向与风向之间的夹角最接近的最接近的第u个监测点数据,表示两个监测点数据差值,反应相对差异,表示第i时刻第u个监测点相对于工厂方向与风向的夹角和第j个监测点相对于工厂方向与风向夹角之间的差异,其值越小,当前两个位置污染物传播越一致,即对应的监测数据越一致,此时反映在相同的第i时刻,第j个监测点相对于第u个监测点的异常,且越大,异常程度越大,表示第i时刻第u个监测点相对于工厂方向与第j个监测点相对于工厂方向之间的夹角,其值越小,污染物传播相对与两监测点的一致性越高,反应异常的程度越大。表示第u个监测点数据在风向影响下的局部异常程度,其值越大,可信度越高,所反映的数据异常可信度越高。表示相对于污染物传播方向最一致的个监测点数据的异常程度,表示的监测点数量,其中本实施例中是预设阈值,设置度。
在实际中风向可能随时发生变化,此时改变污染物相对于每一监测点的传播方向,所以每一监测点监测数据将发生变化,并且监测点之间的差异关系发生变化,从而影响利用监测点数据关系进行异常判断的准确性。
需要说明的是,风向在发生变化时,所有监测点与风向的方向关系发生变化,风向发生变化后,在一段时间内污染物的传播是相对混乱的,导致监测点数据发生相对应的变化,不同监测点数据变化存在差异,最终导致监测点数据之间的差异无法直接反应数据的异常程度,即为第四异常程度。所以风速变化越大,对监测点之间差异所反映的数据异常程度的影响越大。
第四异常程度具体表示为:
其中,表示上述对应第u个监测点在第i时刻与之前相邻的第v时刻风向的变化角度,其值越大,第v时刻的风向改变对第i时刻污染物传播影响越大,所以监测点之间数据差异对数据异常表现程度越小。表示第v时刻与第i时刻的时间间隔,其值越大,表示第v时刻风向变化对第i时刻监测影响越小。表示在相同的第v时刻下,所有监测点相对于工厂方向与风向夹角和第j个监测点相对于工厂方向与风向夹角作比较,选出两个夹角之间的差异小于所对应的监测点,记为第一集合,在相同的第i时刻下,所有监测点相对于工厂方向与风向夹角和第j个监测点相对于工厂方向与风向夹角作比较,选出两个夹角之间的差异小于所对应的监测点,记为第二集合,最后得到第一集合与第二集合中相同监测点数量为;其值越大,第v时刻风向变化相对于第i时刻差异越小,所以第i时刻污染物传播越稳定。表示在第v时刻下的风速大小,其值越大,对监测点的监测数据影响也越大。表示第i时刻之前个时刻对于第u个监测点风向变化对的影响程度,其值越大,对数据异常的表现程度越高。因为实际风向变化时可能存在局部差异,所以对于不同监测点分别进行风向变化分析,获得,即表示个监测点在个时刻风向变化所反映对数据异常的表现程度。表示监测点之间的数据关系所表现的数据异常程度。
进一步说明,在上述分析中,考虑风向变化对多个监测点数据的影响,避免风向在不同监测点的不同表现对监测点之间关系的影响,提高监测点之间数据差异对异常的表现效果,即便于进行更加可靠的异常分析。
由上述分析,根据风向对污染物传播的影响,确定监测点监测数据局部表现以及多个监测点之间关系所反映的数据异常程度。此时综合上述异常关系,确定监测点数据最终异常程度为:
其中,表示第j个监测点第i时刻数据局部关系所反映的异常程度,其值越大,最终异常程度越大;表示第i时刻的第j个监测点与其他监测点之间的数据关系所反映的异常程度。因为和同时受风向的影响,且受的影响,所以二者具有一致相关性,所以采用二者相乘表示最终数据异常程度。
步骤S003:根据目标监测点的最终异常程度进行监测数据筛选和环境分析。
根据上述过程,确定所有监测点数据异常程度,异常程度越大,越不能反应污染物的含量,即无法进行环境监测。所以此时需要根据数据异常程度在所有数据中获得能够反映环境数据的监测样本数据。具体直接删除异常程度大于阈值的数据,剩余数据作为数据样本进行环境监测分析。
我们用正态分布的3原则进行设置阈值,所述阈值在本实施例中设置为,其中分别表示所有的均值和标准差。
最后根据上述过程,利用监测数据的异常程度确定当前环境监测数据中的样本数据,此时根据所获得的数据样本进行环境监测分析,及时获取污染超标,并对超标做出及时处理,专利CN107478778B公开了根据数据样本进行环境监测分析监测的方法,本实施例不再赘述具体过程。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.基于数据样本的重点企业环境分析方法,其特征在于,该方法包括以下步骤:
获取每个时刻每个监测点污染气体浓度的监测数据;将每个时刻下的任意一个监测点记为目标监测点;
根据目标监测点在相邻时刻的监测数据之间的差异获得目标监测点的异常程度;
根据目标监测点的异常程度、目标监测点与相邻时刻风向与风速的变化值获得每个目标监测点的第二异常程度;
将与目标监测点相对于工厂方向与风向之间的夹角最接近的监测点记为目标监测点的参考监测点;根据每个时刻下在参考监测点的第二异常程度、目标监测点相对于工厂方向与风向夹角与参考监测点相对于工厂方向与风向夹角的差异、目标检测点相对于工厂方向与参考监测点相对于工厂方向的夹角以及目标监测点与参考监测点的监测数据的差异获得目标监测点的第三异常程度;
获取每个时刻下所有监测点相对于工厂方向与风向夹角,记为夹角集合S;获取S中与目标监测点相对于工厂方向与风向夹角差异小于预设阈值的夹角所对应的监测点,记为集合S1;获取每个时刻的相邻时刻下所有监测点相对于工厂方向与风向,记为夹角集合S2;获取S2中与目标监测点相对于工厂方向与风向夹角差异小于预设阈值的夹角所对应的监测点,记为集合S3;将S3与S1中相同监测点的数量记为目标监测点在相邻时刻的第一数量;
根据目标监测点的第三异常程度、每个时刻与相邻时刻的时间间隔、目标监测点与相邻时刻的风向变化角度、每个时刻的风速以及第一数量获得目标监测点的第四异常程度;
根据目标监测点的第二异常程度和第四异常程度获得目标监测点的最终异常程度,根据目标监测点的最终异常程度进行监测数据筛选和环境分析。
2.根据权利要求1所述基于数据样本的重点企业环境分析方法,其特征在于,所述目标监测点的异常程度的获取方法包括:
将每个时刻在预设时间窗口内的所有监测数据作为每个时刻的周围最相邻数据,目标监测点在每个时刻的监测数据与每个周围最相邻数据分别计算差异,记为第一差异,所有第一差异的均值记为目标监测点的异常程度。
3.根据权利要求1所述基于数据样本的重点企业环境分析方法,其特征在于,所述目标监测点的第二异常程度的获取方法为:
获取每个时刻工厂到监测点的方向与风向之间的夹角和风速大小;得到每个时刻与相邻时刻风向与风速的变化值,将用风向和风速变化值相乘之后取倒数,求取每个时刻的所有相邻时刻下所述倒数的第一均值;最后将第一均值与风速大小、风向与监测点到工厂方向夹角的倒数以及目标监测点的异常程度的乘积记为目标监测点的第二异常程度。
4.根据权利要求1所述基于数据样本的重点企业环境分析方法,其特征在于,所述目标监测点的第三异常程度的获取方法为:
对于第i时刻的第j个监测点,即为目标监测点,目标监测点的第三异常程度公式为:
其中,表示目标监测点的监测数据,表示目标监测点相对于工厂方向与风向之间的夹角最接近的最接近的第u个监测点数据,即参考监测点;表示第i时刻第u个监测点相对于工厂方向与风向的夹角和目标监测点相对于工厂方向与风向夹角之间的差异;表示第i时刻第u个监测点相对于工厂方向与目标监测点相对于工厂方向之间的夹角;表示参考监测点的个数; 表示第u个监测点数据在风向影响下的第二异常程度;为目标监测点的第三异常程度。
5.根据权利要求4所述基于数据样本的重点企业环境分析方法,其特征在于,所述目标监测点的第四异常程度的获取方法为:
目标监测点的第四异常程度公式为:
其中,表示对应第u个监测点在第i时刻与之前相邻的第v时刻风向的变化角度; 表示第v时刻与第i时刻的时间间隔; 表示第一数量; 表示在第v时刻下的风速大小;表示参考监测点的个数;为预设值;表示为目标监测点的第三异常程度;表示目标监测点的第四异常程度。
6.根据权利要求1所述基于数据样本的重点企业环境分析方法,其特征在于, 所述目标监测点的最终异常程度的获取方法为:
由目标监测点的第二异常程度和第四异常程度相乘得到目标监测点的最终异常程度。
7.根据权利要求1所述基于数据样本的重点企业环境分析方法,其特征在于,所述进行监测数据筛选方法包括为:
用正态分布的原则进行设置阈值,所述阈值设置为,其中,分别表示所有监测点最终异常程度的均值和标准差;我们将大于阈值的监测数据进行删除,将小于的监测数据保留,最后筛选得到用于环境分析的监测数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310409806.5A CN116128260B (zh) | 2023-04-18 | 2023-04-18 | 基于数据样本的重点企业环境分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310409806.5A CN116128260B (zh) | 2023-04-18 | 2023-04-18 | 基于数据样本的重点企业环境分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116128260A true CN116128260A (zh) | 2023-05-16 |
CN116128260B CN116128260B (zh) | 2023-07-11 |
Family
ID=86312158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310409806.5A Active CN116128260B (zh) | 2023-04-18 | 2023-04-18 | 基于数据样本的重点企业环境分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116128260B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116400639A (zh) * | 2023-06-08 | 2023-07-07 | 佛山市星际云数字科技有限公司 | 一种plc采集数据智能清洗方法及*** |
CN116576553A (zh) * | 2023-07-11 | 2023-08-11 | 韦德电子有限公司 | 一种空气调节机数据优化采集方法及*** |
CN116662767A (zh) * | 2023-08-01 | 2023-08-29 | 陕西中科绿能能源研究院有限公司 | 基于多传感器的蒸发冷却机组***温度数据智能采集方法 |
CN117146384A (zh) * | 2023-10-30 | 2023-12-01 | 深圳市宇隆伟业科技有限公司 | 基于数据驱动的空调运行状态监测方法 |
CN117310101A (zh) * | 2023-10-31 | 2023-12-29 | 湖北凌创环保有限公司 | 一种大气污染物浓度监测*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650346A (zh) * | 2020-07-14 | 2020-09-11 | 中科三清科技有限公司 | 大气污染监测数据的自动审核方法、装置及电子设备 |
US20210293651A1 (en) * | 2018-09-11 | 2021-09-23 | Hefei Institute For Public Safety Research, Tsinghua University | Method and device for detecting fault of monitoring device |
CN113777223A (zh) * | 2021-08-12 | 2021-12-10 | 北京金水永利科技有限公司 | 一种大气污染物溯源方法及*** |
CN114329245A (zh) * | 2021-12-29 | 2022-04-12 | 神彩科技股份有限公司 | 一种空气污染的溯源方法、装置、服务器及存储介质 |
CN115861011A (zh) * | 2023-02-15 | 2023-03-28 | 山东优嘉环境科技有限公司 | 一种基于多源数据融合的智慧城市优化管理方法和*** |
CN115860590A (zh) * | 2023-03-02 | 2023-03-28 | 广东慧航天唯科技有限公司 | 一种企业排放污染数据的智能分析预警方法及*** |
CN115876258A (zh) * | 2023-02-14 | 2023-03-31 | 山东大佳机械有限公司 | 基于多源数据的畜禽养殖环境异常监测及报警*** |
-
2023
- 2023-04-18 CN CN202310409806.5A patent/CN116128260B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210293651A1 (en) * | 2018-09-11 | 2021-09-23 | Hefei Institute For Public Safety Research, Tsinghua University | Method and device for detecting fault of monitoring device |
CN111650346A (zh) * | 2020-07-14 | 2020-09-11 | 中科三清科技有限公司 | 大气污染监测数据的自动审核方法、装置及电子设备 |
CN113777223A (zh) * | 2021-08-12 | 2021-12-10 | 北京金水永利科技有限公司 | 一种大气污染物溯源方法及*** |
CN114329245A (zh) * | 2021-12-29 | 2022-04-12 | 神彩科技股份有限公司 | 一种空气污染的溯源方法、装置、服务器及存储介质 |
CN115876258A (zh) * | 2023-02-14 | 2023-03-31 | 山东大佳机械有限公司 | 基于多源数据的畜禽养殖环境异常监测及报警*** |
CN115861011A (zh) * | 2023-02-15 | 2023-03-28 | 山东优嘉环境科技有限公司 | 一种基于多源数据融合的智慧城市优化管理方法和*** |
CN115860590A (zh) * | 2023-03-02 | 2023-03-28 | 广东慧航天唯科技有限公司 | 一种企业排放污染数据的智能分析预警方法及*** |
Non-Patent Citations (4)
Title |
---|
YANG, XIAOQIN等: "Monitoring water color anomaly of lakes based on an integrated method using Ladsat-8 OLI images", 《INTERNATIONAL JOURNAL OF DIGITAL EARTH》, vol. 15, no. 1 * |
李茂洁;: "对环境监测数据中异常数据的处理方法研究", 居舍, no. 16 * |
诸葛燕;方海泉;: "基于大数据的水量监测***有效性分析", 浙江海洋大学学报(自然科学版), no. 06 * |
赵晓庆: "污染源自动监测动态管控***的研究与实现", 《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》, no. 12 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116400639A (zh) * | 2023-06-08 | 2023-07-07 | 佛山市星际云数字科技有限公司 | 一种plc采集数据智能清洗方法及*** |
CN116400639B (zh) * | 2023-06-08 | 2023-08-11 | 佛山市星际云数字科技有限公司 | 一种plc采集数据智能清洗方法及*** |
CN116576553A (zh) * | 2023-07-11 | 2023-08-11 | 韦德电子有限公司 | 一种空气调节机数据优化采集方法及*** |
CN116576553B (zh) * | 2023-07-11 | 2023-09-22 | 韦德电子有限公司 | 一种空气调节机数据优化采集方法及*** |
CN116662767A (zh) * | 2023-08-01 | 2023-08-29 | 陕西中科绿能能源研究院有限公司 | 基于多传感器的蒸发冷却机组***温度数据智能采集方法 |
CN116662767B (zh) * | 2023-08-01 | 2023-10-13 | 陕西中科绿能能源研究院有限公司 | 基于多传感器的蒸发冷却机组***温度数据智能采集方法 |
CN117146384A (zh) * | 2023-10-30 | 2023-12-01 | 深圳市宇隆伟业科技有限公司 | 基于数据驱动的空调运行状态监测方法 |
CN117146384B (zh) * | 2023-10-30 | 2024-01-30 | 深圳市宇隆伟业科技有限公司 | 基于数据驱动的空调运行状态监测方法 |
CN117310101A (zh) * | 2023-10-31 | 2023-12-29 | 湖北凌创环保有限公司 | 一种大气污染物浓度监测*** |
Also Published As
Publication number | Publication date |
---|---|
CN116128260B (zh) | 2023-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116128260B (zh) | 基于数据样本的重点企业环境分析方法 | |
CN111898691A (zh) | 一种河流突发水污染预警溯源方法、***、终端及介质 | |
CN109446189A (zh) | 一种工业参数离群点检测***及方法 | |
US20210216877A1 (en) | Ultra-high sensitive target signal detection method based on noise analysis using deep learning based anomaly detection and system using the same | |
CN112799898B (zh) | 基于分布式故障检测的互联***故障节点定位方法及*** | |
CN111289697B (zh) | Cems运行状态监控与数据作假识别方法及*** | |
AU2020102518A4 (en) | A method of gas sensor calibration based on linear optimization | |
CN107764976B (zh) | 土壤氮素快速诊断方法及在线监测*** | |
CN111414571A (zh) | 一种大气污染物监测方法 | |
CN117195135B (zh) | 一种水污染异常溯源检测方法及*** | |
CN116985183B (zh) | 一种近红外光谱分析仪质量监测管理方法及*** | |
CN114217025B (zh) | 评估空气质量浓度预测中气象数据对其影响的分析方法 | |
CN115659195A (zh) | 一种大气污染在线识别方法 | |
CN116359285A (zh) | 一种基于大数据的油气浓度智能检测***及方法 | |
JPH06117600A (ja) | ガス、蒸気等の漏洩地点および漏洩量推定システム | |
CN111272950B (zh) | 一种智能污染物追踪*** | |
CN115656063A (zh) | 基于高光谱数据的目标污染物含量分析方法及装置 | |
CN114814092A (zh) | 基于bp神经网络的ip指标测量方法 | |
CN113804595A (zh) | 一种多参数空气质量监测*** | |
CN114004138A (zh) | 基于大数据人工智能的建筑监测方法、***及存储介质 | |
CN114626162B (zh) | 一种接触球轴承损耗程度定量识别方法 | |
CN117607019B (zh) | 电力金具表面智能检测方法和检测*** | |
CN117336195B (zh) | 一种基于雷达图法的入侵检测模型综合性能评估方法 | |
Amato et al. | Decision trees in time series reconstruction problems | |
CN113705547B (zh) | 环境弄虚作假行为识别动态管控方法、装置 |
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 |