CN111339499B - 一种基于时间轨迹相似度的污水处理厂运行状态评估方法 - Google Patents
一种基于时间轨迹相似度的污水处理厂运行状态评估方法 Download PDFInfo
- Publication number
- CN111339499B CN111339499B CN202010245116.7A CN202010245116A CN111339499B CN 111339499 B CN111339499 B CN 111339499B CN 202010245116 A CN202010245116 A CN 202010245116A CN 111339499 B CN111339499 B CN 111339499B
- Authority
- CN
- China
- Prior art keywords
- plant
- water
- sewage
- similarity
- effluent
- 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
- 239000010865 sewage Substances 0.000 title claims abstract description 115
- 238000011156 evaluation Methods 0.000 title claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 140
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000002159 abnormal effect Effects 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims abstract description 7
- 238000012544 monitoring process Methods 0.000 claims abstract description 7
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 34
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 12
- 229910052757 nitrogen Inorganic materials 0.000 claims description 6
- 238000005315 distribution function Methods 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 description 8
- 238000004065 wastewater treatment Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000003911 water pollution Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- 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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Analysis (AREA)
- Marketing (AREA)
- Pure & Applied Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Primary Health Care (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Activated Sludge Processes (AREA)
Abstract
本发明公开了一种基于时间轨迹相似度的污水处理厂运行状态评估方法,包括:建立多个污水厂实时在线监测出水水质数据集,保留修正数据集,并按照不同污水厂不同出水水质指标进行标记,选择所有污水厂和单个污水厂以特定时间周期对水质数据的数据集进行分段;利用厂群总体水质数据数据集和单个水厂水质数据数据集求出相似度,得到所有水厂相似度的数据集;利用分位数将得到所有水厂水质相似度的数据集进行分类,设置上下侧分位数,根据水质相似度的污水厂判定其运行状态。本发明根据相似度值的大小与污水厂实际运行来表征污水厂的行为特性。该方法能够快速识别具有异常风险的污水厂,为污水厂监管评估提供一种有效可行的方法。
Description
技术领域
本发明属于污水处理领域,涉及一种通过对区域中厂群总体污水厂和单独污水厂出水水质数据的相关性分析,从而表征出污水厂特性及异常状态识别的方法。
背景技术
近年以来,环境水污染问题愈发得到人们的重视,目前,通常通过有效的污水厂的在线水质指标排放在线数据方法来评估治理环境水污染问题的实施效果。然而,通过对这些在线数据进行评估,从而评估所有污水处理厂治理环境水污染问题的有效性,其难度和成本是特别大的。这也是因为一些污水处理厂存在不法行为,比如:以临时达标为目的,对监测数据进行作假。这使得环保部门对每个污水处理厂进行合理评估的难度加大,这种造假行为大大的影响了对环保部门的监管能力。对此,我们提出了一种通过对一定区域内多个污水厂出水水质数据的数学分析方法,来为当前相关的环保部门提供一个可行有效的监管方法和绩效评估方法,并表征出污水厂的行为特征。
发明内容
本发明的目的在于提供一种可监管污水厂排放行为和对污水厂进行绩效评估的方法。我们根据区域内污水厂厂群总体数据以特定时间周期的水质数据来计算相似度,通过相似度计算出的上下侧分位数,识别出具有数据异常或运行情况良好的污水厂,同时根据相似度的大小来表征出每个污水厂的运行状况。
本发明的目的是通过下述技术方案来实现的。
一种基于时间轨迹相似度的污水处理厂运行状态评估方法,包括下述步骤:
(1)获取区域内多个污水厂实时在线监测的出水水质数据,建立数据集;
(2)剔除数据集中所包含的缺失值和异常值,保留修正数据集,并将修正数据集按照不同污水厂不同出水水质指标进行标记如下:
Tij表示第i个污水厂第j个出水水质指标的数据集;i=1,2,3...;j=1,2,3...;
Tj表示所有污水厂第j个出水指标的总体出水水质数据集;
(3)选择所有污水厂中第j个出水指标的厂群总体出水水质数据集Tj,以特定时间周期对水质数据的数据集进行分段;
(4)将第i个污水厂第j个出水水质指标的数据集Tij单独以特定时间周期对水质数据的数据集进行分段,获取到单个污水厂的特定时间周期水质数据的数据集;
(5)利用厂群总体水质数据数据集Tj和单个水厂水质数据数据集Tij求出厂群与单个水厂之间的相似度,得到所有水厂相似度的数据集;
(6)利用分位数将得到所有水厂水质相似度的数据集进行分类,设置上侧分位数a和下侧分位数b,在分位数为0-a%水质相似度的污水厂,表明其运行状态差;在分位数b%-1水质相似度的污水厂,表明其运行状态良好。
进一步,所述步骤1)中,出水水质数据包括BOD、氨氮、pH值、COD 和总氮。
进一步,所述步骤2)中,缺失值为数据在收集过程中导致部分数据为空值;所述异常值为监测设备器件出现异常,产生过大的数据漂移,导致数据变化异常大或为负值。
进一步,所述步骤3)、4)中,特定时间周期可为以节气、月份或星期为节点的时间,以节气为时间周期的,将时间段通过24个节气分为24个时间段。
进一步,所述步骤5)中,利用厂群总体与单个水厂修正后的包括BOD、氨氮、pH值、COD和总氮水质数据,求出特定时间周期的特征值,所述特征值为平均值或中位数。
进一步,所述步骤5)中,利用特征值求出厂群总体与单个水厂之间的水质相似度,求相似度的计算通过求出总体水厂与每个水厂之间的最大互信息系数 MIC值得到:
式中,I[X;Y]为X和Y的联合概率,p(X)为X的概率密度分布,p(Y) 为Y的概率密度分布,p(X,Y)为X、Y的联合概率密度分布函数,|X|为X 方向共被分成的段数;|Y|为Y方向共被分成的段数;B为分区网格的大小的限制值,|X||Y|<B为X方向段数乘Y方向段数的值小于网格大小限制值时间的情况,min(|X|,|Y|)为X方向和Y方向被分成的段数中较小的段数,MIC为所求得的最大互信息系数。
进一步,所述步骤6)中,设置上侧分位数a为10%,下侧分位数b为90%。
进一步,所述步骤6)中,根据相似度的大小与污水厂实际运行相结合,来表征污水厂的行为特性如下:
对于出水NH,通过分位数将水厂的相似度分为三个范围0-10%、10%-90%和90%-100%,根据相似度大小可以表征如下污水厂的行为特性:
(1)在分位数为90%-1的污水厂,即水质相似度越高的水厂,单个水厂的运行情况与厂群总体的运行情况相似,水厂出水氨氮规律符合厂群总体出水氨氮规律,表明其运行状态良好;
(2)在分位数为0-10%的污水厂,即水质相似度低的水厂,单个水厂的运行情况与厂群总体的运行情况不相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态差;
(3)在分位数为10%-90%的污水厂,即水质相似度趋近于中间值的水厂,单个水厂运行情况与厂群总体的运行情况部分相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态中等。
本发明具有以下优点:
1)成本低,迅速便捷。
该方法只需要收集污水厂的出水水质数据,通过该方法分析,可对污水厂的运行状态进行监管评估,减少了传统中的耗时长等问题,相关人员不用亲自到现场即完成对污水厂的监管,减少了中间的成本。
2)具有普适性与公平性。
普适性为在样本量足够的情况下,能够捕获不同区域水质数据的关联,而不限定于特定的水质类型。所谓公平性,是指在样本量足够大时能为不同的污水厂相似的相关关系给出相近的系数。
本发明根据相似度值的大小与污水厂实际运行来表征污水厂的行为特性,该方法能够快速识别具有异常风险的污水厂,为污水厂监管评估提供一种有效可行的方法。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的不当限定,在附图中:
图1为本区域总体污水厂以出水氨氮为水质指标求出的与单个污水厂的最大互信息系数(MIC值)。
图2(a)为MIC值最高值(0.79)的0号污水厂,以二十四节气分段的氨氮指标中位数的散点图。
图2(b)为MIC值最低值(0.05)的91号污水厂,以二十四节气分段的氨氮指标中位数的散点图。
图2(c)为2018年报道过有异常的且MIC值为(0.46)的51号污水厂,以二十四节气分段的氨氮指标中位数的散点图。
图2(d)为MIC值为中位数(0.31)的1号污水厂,以二十四节气分段的氨氮指标中位数的散点图。
具体实施方式
下面通过附图及实施例对本发明做进一步的说明,在此本发明的示意性实施例以及说明一解释本发明,但并不作为对本发明的限定。
本发明通过发现一种基于时间轨迹相似度的污水处理厂运行状态评估方法,进而选择合适的管理运行,包括以下步骤:
(1)获取区域内多个污水厂的实时在线监测的出水水质数据,出水水质数据包括:BOD、氨氮、pH值、COD和总氮等多个指标数据,建立数据集。
(2)剔除数据集中所包含的缺失值和异常值,数据在收集过程中,会由于人为或者设备原因,导致部分数据为空值,我们称为缺失值,缺失排除在外,不计入数据样本中;同时监测设备器件会出现异常,产生过大的数据漂移,导致数据变化异常大或者值为负值,对于这一部分数据,我们称为异常值,为避免数据分析结果偏差过大,我们根据仪器的所测范围选择正常范围内的数据。保留修正数据集,并将修正数据集按不同污水厂不同出水水质指标进行标记,如下:Tij(i=1,2,3...j=1,2,3...),表示第i个污水厂第j个出水水质指标的数据集; Tj(j=1,2,3...)表示所有污水厂第j个出水指标的总体出水水质数据集。
(3)选择所有污水厂中第j个出水指标的厂群总体出水水质数据集Tj,以特定时间周期对水质数据的数据集进行分段。
(4)将第i个污水厂第j个出水水质指标的数据集Tij单独以特定时间周期对水质数据的数据集进行分段,获取到每个污水厂的特定时间周期水质数据的数据集。
(5)利用厂群总体水质数据数据集Tj和单个水厂水质数据数据集Tij求出厂群与单个水厂之间的相似度,得到所有水厂相似度的数据集;利用厂群总体与单个水厂修正后的包括BOD、氨氮、pH值、COD和总氮水质数据,求出特定时间周期的特征值,特征值为平均值或中位数。
其中,利用特征值求出厂群总体与单个水厂之间的水质相似度,求相似度的计算通过求出总体水厂与每个水厂之间的最大互信息系数MIC值得到:
式中,I[X;Y]为X和Y的联合概率,p(X)为X的概率密度分布,p(Y) 为Y的概率密度分布,p(X,Y)为X、Y的联合概率密度分布函数,|X|为X 方向共被分成的段数;|Y|为Y方向共被分成的段数;B为分区网格的大小的限制值,|X||Y|<B为X方向段数乘Y方向段数的值小于网格大小限制值时间的情况,min(|X|,|Y|)为X方向和Y方向被分成的段数中较小的段数,MIC为所求得的最大互信息系数。
(6)利用分位数将得到所有水厂水质相似度的数据集进行分类,设置上侧分位数a为10%,下侧分位数b为90%,并且利用MIC值的数值大小与污水厂实际运行相结合,来表征污水厂的行为特性。
根据相似度的大小与污水厂实际运行相结合,来表征污水厂的行为特性如下:
对于出水NH,通过分位数将水厂的相似度分为三个范围0-10%、10%-90%和90%-100%,根据相似度大小可以表征如下污水厂的行为特性:
(1)在分位数为90%-1的污水厂,即水质相似度越高的水厂,单个水厂的运行情况与厂群总体的运行情况相似,水厂出水氨氮规律符合厂群总体出水氨氮规律,表明其运行状态良好;
(2)在分位数为0-10%的污水厂,即水质相似度低的水厂,单个水厂的运行情况与厂群总体的运行情况不相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态差;
(3)在分位数为10%-90%的污水厂,即水质相似度趋近于中间值的水厂,单个水厂运行情况与厂群总体的运行情况部分相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态中等。
下面通过具体实施例进一步说明本发明效果。
1)通过西安绿标水环境科技有限公司提供的WBM400型污水处理智慧运行工作站,收集本地区92个污水厂2019年实时在线监测的出水数据;在此以出水氨氮为数据样本进行说明。
2)对出水氨氮数据进行预处理。
数据在收集过程中,导致部分数据为空值的缺失值不计入数据样本中;同时对仪器本身出现过大的数据漂移导致的数据变化异常大或者值为负值的异常值,进行剔除。因此为避免数据分析结果偏差过大,我们根据仪器的所测范围,选择在范围0-200mg/L的出水氨氮数据。
假定,Ti1(i=1,2,3…92)表示第i个污水厂出水氨氮的数据集,T1为所有污水厂的出水氨氮数据集;
(3)特性时间周期选用节气,求出T1每个节气时间段内的出水水质数据的中位数,总计24个中位数;
(4)求出Ti1(i=1,2,3…92)92个污水厂每个节气时间段内的出水水质数据的中位数,每个水厂24个中位数;
(5)相关度算法选用选最打互信息数,求出总体水厂与每个水厂之间的最大互信息系数(MIC值):
(6)利用分位数将污水厂进行分类,分位数0%-10%、10%-90%和90%-100%三类,分位数在0%-10%的为识别出的异常水厂包括0号水厂;分位数在90%-10 0%的为识别出来运行情况良好的水厂包括0、31、33、36、37、44、47、48、5 7、75号水厂,其余水厂为分位数在10%-90%的水厂。
图1为本区域92个污水厂以出水氨氮为水质指标求出的厂群与单个污水厂的最大互信息系数(MIC值),其中17号和32号为MIC值0的无效水厂,因为其数据有一定时间段的缺失,MIC值最大为0号水厂0.79,MIC值最小为91 号水厂0.05,MIC值中位数为1号水厂0.31。
图2(a)为MIC值最高值(0.79)的0号污水厂,单个水厂的运行情况与总体水厂的运行情况相似,水厂出水水质规律符合总体水厂出水水质规律,数据造假的嫌疑低,原因可能为污水处理技术水平高,也可能其污水厂地区处于符合总体污水厂类型的典型地区,且运行管理水平较好。
图2(b)为MIC值最低值(0.05)的91号污水厂,单个水厂的运行情况与总体水厂的运行情况不相似,水厂出水水质规律不符合总体水厂出水水质规律,或者存在做假风险,将数据较大值去除,也可能其污水厂地区处于不符合总体污水厂类型的严寒等地区。
图2(c)为2018年报道过有异常的且MIC值为(0.46)的51号污水厂,虽然此水厂MIC值不低,但是也说明它的运行情况与总体污水厂有部分不同,此污水厂运行可能存在一定时间段内与通常污水厂运行不同的情况,有一定不规律的水质波动。
图2(d)为MIC值为中位数(0.31)的1号污水厂,单个水厂运行情况与总体水厂的运行情况部分相似部分不相似,其原因较为复杂,可能是由部分节气时间段数据差异造成,或者是在某一段时间范围内存在工艺改造。
由此可以看出通过厂群与单个污水厂的相关度,能够表征污水厂的行为特征,识别出具有异常风险的污水厂,为相关环保部门对污水厂的监管和对污水处理厂的评估提供了可行有效的办法。
本发明并不局限于上述实施例,在本发明公开的技术方案的基础上,本领域的技术人员根据所公开的技术内容,不需要创造性的劳动就可以对其中的一些技术特征作出一些替换和变形,这些替换和变形均在本发明的保护范围内。
Claims (8)
1.一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,包括下述步骤:
(1)获取区域内多个污水厂实时在线监测的出水水质数据,建立数据集;
(2)剔除数据集中所包含的缺失值和异常值,保留修正数据集,并将修正数据集按照不同污水厂不同出水水质指标进行标记如下:
Tij表示第i个污水厂第j个出水水质指标的数据集;i=1,2,3...;j=1,2,3...;
Tj表示所有污水厂第j个出水指标的总体出水水质数据集;
(3)选择所有污水厂中第j个出水指标的厂群总体出水水质数据集Tj,以特定时间周期对水质数据的数据集进行分段;
(4)将第i个污水厂第j个出水水质指标的数据集Tij单独以特定时间周期对水质数据的数据集进行分段,获取到单个污水厂的特定时间周期水质数据的数据集;
(5)利用厂群总体水质数据数据集Tj和单个水厂水质数据数据集Tij求出厂群与单个水厂之间的相似度,得到所有水厂相似度的数据集;
(6)利用分位数将得到所有水厂水质相似度的数据集进行分类,设置上侧分位数a和下侧分位数b,在分位数为0-a%水质相似度的污水厂,表明其运行状态差;在分位数b%-1水质相似度的污水厂,表明其运行状态良好。
2.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤1)中,出水水质数据包括BOD、氨氮、pH值、COD和总氮。
3.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤2)中,缺失值为数据在收集过程中导致部分数据为空值;所述异常值为监测设备器件出现异常,产生过大的数据漂移,导致数据变化异常大或为负值。
4.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤3)、4)中,特定时间周期可为以节气、月份或星期为节点的时间,以节气为时间周期的,将时间段通过24个节气分为24个时间段。
5.根据权利要求4所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤5)中,利用厂群总体与单个水厂修正后的包括BOD、氨氮、pH值、COD和总氮水质数据,求出特定时间周期的特征值,所述特征值为平均值或中位数。
6.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤5)中,利用特征值求出厂群总体与单个水厂之间的水质相似度,求相似度的计算通过求出总体水厂与每个水厂之间的最大互信息系数MIC值得到:
式中,I[X;Y]为X和Y的联合概率,p(X)为X的概率密度分布,p(Y)为Y的概率密度分布,p(X,Y)为X、Y的联合概率密度分布函数,|X|为X方向共被分成的段数;|Y|为Y方向共被分成的段数;B为分区网格的大小的限制值,|X||Y|<B为X方向段数乘Y方向段数的值小于网格大小限制值时间的情况,min(|X|,|Y|)为X方向和Y方向被分成的段数中较小的段数,MIC为所求得的最大互信息系数。
7.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤6)中,设置上侧分位数a为10%,下侧分位数b为90%。
8.根据权利要求1所述的一种基于时间轨迹相似度的污水处理厂运行状态评估方法,其特征在于,所述步骤6)中,根据相似度的大小与污水厂实际运行相结合,来表征污水厂的行为特性如下:
对于出水NH,通过分位数将水厂的相似度分为三个范围0-10%、10%-90%和90%-100%,根据相似度大小可以表征如下污水厂的行为特性:
(1)在分位数为90%-1的污水厂,即水质相似度越高的水厂,单个水厂的运行情况与厂群总体的运行情况相似,水厂出水氨氮规律符合厂群总体出水氨氮规律,表明其运行状态良好;
(2)在分位数为0-10%的污水厂,即水质相似度低的水厂,单个水厂的运行情况与厂群总体的运行情况不相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态差;
(3)在分位数为10%-90%的污水厂,即水质相似度趋近于中间值的水厂,单个水厂运行情况与厂群总体的运行情况部分相似,水厂出水氨氮规律不符合总体水厂出水氨氮规律,表明其运行状态中等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010245116.7A CN111339499B (zh) | 2020-03-31 | 2020-03-31 | 一种基于时间轨迹相似度的污水处理厂运行状态评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010245116.7A CN111339499B (zh) | 2020-03-31 | 2020-03-31 | 一种基于时间轨迹相似度的污水处理厂运行状态评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111339499A CN111339499A (zh) | 2020-06-26 |
CN111339499B true CN111339499B (zh) | 2023-04-07 |
Family
ID=71186351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010245116.7A Active CN111339499B (zh) | 2020-03-31 | 2020-03-31 | 一种基于时间轨迹相似度的污水处理厂运行状态评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111339499B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112037106B (zh) * | 2020-08-07 | 2023-12-15 | 汉威科技集团股份有限公司 | 一种基于特征互相关性和概率密度的数据异常分析方法 |
CN113988359B (zh) * | 2021-09-08 | 2024-08-13 | 中南大学 | 基于非对称Laplace分布的风功率预测方法和系统 |
CN114662981B (zh) * | 2022-04-15 | 2023-01-03 | 广东柯内特环境科技有限公司 | 基于大数据应用的污染源企业监管方法 |
CN117349777B (zh) * | 2023-12-04 | 2024-02-23 | 安徽新宇环保科技股份有限公司 | 一种水环境在线监测数据真伪智能识别系统及方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011145846A (ja) * | 2010-01-14 | 2011-07-28 | Hitachi Ltd | 異常検知方法、異常検知システム、及び異常検知プログラム |
KR20150111672A (ko) * | 2014-03-26 | 2015-10-06 | 부산대학교 산학협력단 | 동일한 방류수계를 공유하는 다수의 하수처리장의 통합운영관리장치 및 그 방법 |
-
2020
- 2020-03-31 CN CN202010245116.7A patent/CN111339499B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111339499A (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111339499B (zh) | 一种基于时间轨迹相似度的污水处理厂运行状态评估方法 | |
CN112906738B (zh) | 一种水质检测及处理方法 | |
CN111428985B (zh) | 区域水生态承载力的评估方法 | |
CN117786584B (zh) | 基于大数据分析的畜牧业水源污染监测预警方法及系统 | |
CN115099677B (zh) | 尾矿库安全生产风险分级预警方法 | |
CN117952308A (zh) | 一种区域可持续发展动态监测与评估的方法及系统 | |
CN118469157B (zh) | 基于大数据分析的河流综合管理系统 | |
CN112101969A (zh) | 一种基于时序滑动窗口离散系数的环保数据造假检测方法 | |
CN117498958A (zh) | 一种基于大数据的通信传输设备数据监管方法及系统 | |
CN115907461A (zh) | 一种基于机理推导方程的电力工程方法 | |
CN117274827B (zh) | 一种智慧环境环保远程实时监测预警方法及系统 | |
CN110807174B (zh) | 一种基于统计分布的污水厂厂群出水分析及异常识别方法 | |
CN117436688A (zh) | 一种基于电力营销数据的电费回收风险预警方法及系统 | |
CN112634113A (zh) | 一种基于动态滑窗的污染废气关联分析方法 | |
CN104519511B (zh) | 一种检测通信网络小区场景突变的方法和装置 | |
CN116523174A (zh) | 基于自组织神经网络的污水处理厂工况评估方法及系统 | |
CN113516850B (zh) | 基于空间句法分析的管道沿线交通流量数据采集方法 | |
CN112309129B (zh) | 一种基于卡口车牌识别的车辆出行特征分类的od矩阵形成方法 | |
CN115146720A (zh) | 基于机器学习的河湖水环境数据融合与样本标注方法、系统及计算机设备 | |
CN109741091A (zh) | 基于基础负荷消减策略的用户负荷分类方法 | |
CN111178772A (zh) | 一种基于犹豫模糊集的灾害风险评估方法和系统 | |
CN113378903B (zh) | 一种广流域水体环境数据集成方法及系统 | |
CN116319262B (zh) | 一种电能计量装置通信故障数据检测模型 | |
CN118350285B (zh) | 一种基于分区分类阈值的河流生态基流目标确定方法 | |
CN113988571A (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 |