CN106709662B - 一种电力设备运行工况划分方法 - Google Patents
一种电力设备运行工况划分方法 Download PDFInfo
- Publication number
- CN106709662B CN106709662B CN201611263826.2A CN201611263826A CN106709662B CN 106709662 B CN106709662 B CN 106709662B CN 201611263826 A CN201611263826 A CN 201611263826A CN 106709662 B CN106709662 B CN 106709662B
- Authority
- CN
- China
- Prior art keywords
- data
- clustering
- optimal
- value
- algorithm
- 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
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 30
- 238000012549 training Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 3
- 239000013598 vector Substances 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000007418 data mining Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 239000006185 dispersion Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 230000001680 brushing effect Effects 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000010298 pulverizing process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013024 troubleshooting Methods 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/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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/30—Information retrieval; Database structures therefor; File system structures therefor of unstructured textual data
- G06F16/35—Clustering; Classification
-
- 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/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种电力设备运行工况划分方法,包括构建均值聚类模型的训练数据;对训练数据进行预处理,删除无效数据,并且进行归一化处理;采用K‑means++算法确定初始聚类中心;确定聚类划分类别的上下限;循环实现多种类别的聚类过程;对聚类结果进行评价分析,确定最优的聚类类别数K值;选取最优K值模型作为最优工况划分结果等步骤,提高了效率和准确率。
Description
技术领域
本发明涉及火电厂运行工况划分领域,具体涉及一种电力设备运行工况划分方法。
背景技术
我国燃煤火电厂中,锅炉汽水系统应用广泛,而高温过热器作为汽水系统的重要部件,在整个锅炉正常运行中起着重要作用。而在实际火电厂中由于高温过热器在运行中时常会出现受热面结焦超温、安装不合理等问题,长期定点吹刷,蒸汽品质长期不合格等等故障,严重威胁着发电厂的安全和经济运行。目前,国内大部分火电厂的汽水系统高温过热器的正常运行处于依靠专家或运行人员的经验来判断故障的发生,所以在当前大力提倡环保和节约能源的大趋势下,确保高温过热器对实时数据的检测过程中保持与历史相同工况下数据的对比,进而辅助专业人员通过判据实时监测运行状况,提前预测故障发生尤为重要。其中,最优工况的划分是实现这个过程的最重要的前提条件。
当前对于电厂工况划分的方法主要有两大类:基于专家经验的单指标工况划分,基于有监督的数据挖掘方法。
现有技术中已经公开的方式有:《火电厂球磨机制粉系统运行工况划分方法》,CN103301926B;《配电网线路工况分类识别方法》,CN105277852A;《基于支持向量机的风电机组运行工况划分方法》,太阳能学报,2010年第31卷第9期;《风电机组振动监测中运行工况划分方法研究》,华北电力大学,2010年;《基于小波特征的配电网线路运行工况识别》,电网与清洁能源,2010年10月;《基于数据挖掘工况分析技术的研究》,电网与清洁能源,2010年10月;《基于数据挖掘的电站运行优化理论研究与应用》,华北电力大学,2006年;《基于工况划分的火电机组运行优化规则提取》,2009年全国博士生学术会议-电站自动化信息化。
基于专家经验的单指标工况划分方法是根据某一具体的实际应用,电力专家通过自己的专业知识和多年的电厂实践经验,通过单一指标对工况进行划分。这一方法需要较高的专业知识做基础,并且在电厂中需要工作多年,积累一定量的实际经验,才能对电厂内各个系统设备的运行状况有一定的深入理解。然而,电厂发电过程本来就是一个复杂系统,单纯靠人的经验是无法精确判定设备运行状况防止故障发生的。因为,每个人的知识和经验有限,并且现场的运行人员的经验更有限,特别是一些工作年限较短的人员,即使是有丰富的经验的专家也不可能完全解决电厂运行的所有问题。并且从专家的角度来看,只能从单一指标出发,选取有代表性的某一指标进行工况划分,然后对于某些复杂过程,影响正常工况划分的指标不仅仅只有一个,而对于高维的情况只依靠专家的经验很难得到一个最优的工况划分结果。
基于有监督的数据挖掘方法是通过挖掘数据内部的规律,达到划分最优工况的目的。一般对于某些电厂的应用场景,正常的运行工况有多种,并且专业上已经给出每一类工况什么特征,把各个工况的数据作为训练样本集,通过有监督的挖掘学习,得到分类器模型,再去对测试数据进行判别测试。它相对于依靠专家经验的方法,有了质的飞跃。但对于复杂系统的需要专家花费大量的时间对各类进行标记,时间成本和人力成本较高。不过,整体上来说,它可以解决多指标的工况划分问题,不再只依靠专业人员的知识,并且可以处理高维的问题,并可以对测试数据判定是在历史运行过程中的哪一个工况。具体的算法有:支持向量机、神经网络、判别分析、贝叶斯分类器、PCA等。
基于专家经验的方法是建立在专工大量的专业知识和丰富的工作经验的基础上来进行工况划分的。然而,这在现实电厂的实际运行中是不现实的。因为,运行人员的经验不足,专业领域知识有限是很普遍的现象,所以实际的可行性较差。专业人员给出的划分标准是基于单指标来划分的,然而,这在实际电厂应用中是有很大缺陷的,因为有些复杂系统不仅仅依靠某一个指标来工况划分,对于多指标或高维的情况,各个指标之间的关联关系较复杂,仅依靠专业人员自身知识和经验解决较困难。
基于有监督的数据挖掘方法虽然跨越了维度的限制,但是还是脱离不了专业知识的人为主观因素的影响,人为的主观因素产生的误差是不可避免的。但是,如果某一个设备的运行状况专业人员给不出一个工况划分标准,从专业角度给不了指导实际生产的经验,那么,有监督的数据挖掘方法使用将意义不大;并且若是专业人员能够给出相对应的类别标记,对于大量的数据来说,花费的人力、物力也是非常大的,况且人为的操作本身就有它的不可避免的误差。一个模型的准确,源于给定的训练样本的准确性。在电厂的实际应用中,人所能理解的系统的运行状况是有限的,有很大的一部分系统运行状况远远超出了人的理解,专业更给不出有依据的知识。这也是当前急需要解决的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种电力设备运行工况划分方法及评价方法。它是基于无监督学习,仅对数据进行分析就可以得到最优工况结果的方法,采用聚类的思想直接对研究对象进行建模并对建模结果进行评价最终得到最优模型。因此,它可以不用依赖专业人员的知识和丰富的实践经验,从数据角度出发,挖掘出数据内部规律,找到最优工况。并且它可以对多指标划分工况,更符合实际应用要求,最终使得工况划分的准确率会更高。
本发明提供了一种电力设备运行工况划分方法,依次包括如下步骤:
步骤1:从电厂数据库中读取相关设备足够时间的历史运行状态数据,利用设定的筛选规则从历史运行状态数据中筛选出运行状态良好的健康设备数据作为构建均值聚类模型的训练数据;
步骤2:对训练数据进行预处理,删除无效数据,并且进行归一化处理;
步骤3:采用K-means++算法确定初始聚类中心;
步骤4:确定聚类划分类别的上下限;
步骤5:循环实现多种类别的聚类过程;
步骤6:对聚类结果进行评价分析,确定最优的聚类类别数K值;
步骤7:选取最优K值模型作为最优工况划分结果。
进一步地,步骤1还包括假若从一个具有N个观察指标的设备数据库截取M个时刻的设备数据作为训练数据,则在m时刻的各个观察指标数据xmi能够看作为一个N维的向量,其中1≤i≤N,表示为:
X(tm)=[xm1,xm2,xm3,…,xmN]T;
则该训练数据文件应当为M*N的矩阵形式。具体形式如下:
进一步地,步骤2具体步骤为将包含无效数据的观测向量进行删除,采用的公式为:
其中,datau、datal为前后边界数据,momentl、momentu为前后边界时间,datan、momentn为NaN数据及其时间;
对各个参数的测量值根据各自的极值进行归一化处理,采用如下公式对各个参数的测量值进行预处理,将各个参数的测量值全部映射到[0,1]区间内:
其中max(data)和min(data)分别表示对应测量值的最大值和最小值,data表示对应各个参数的测量值。
进一步地,步骤3中确定初始聚类中心具体步骤为:
(1)从输入的数据点集合中随机选择一个点作为第一个聚类中心;
(2)对于数据集中的每一个点x,计算它与最近聚类中心的距离D(X);
(3)选择一个新的数据点作为新的聚类中心,选择的原则是:D(X)较大的点,被选取作为聚类中心的概率较大;
(4)重复(2)和(3)直到K个聚类中心被选出来。
进一步地,步骤5包括如下步骤:
(1)选取k-means++算法得到K个初始聚类中心:Z1(1),Z2(1),...,ZK(1);
(2)按最小距离原则将其余样品分配到K个聚类中心中的某一个,即:
min{||X-Zi(k)||,i=1,2,…,K}=|X-Zj(k)||=Dj(k),则X∈Sj(k);
(3)计算各个聚类中心的新向量值:Zj(k+1)j=1,2,…,K;
Nj为第j类的样本数;
(4)如果Zj(k+1)≠Zj(k),j=1,2,…,K,则回到步骤(2),将模式样本逐个重新分类,重复迭代计算;
(5)如果Zj(k+1)=Zj(k),j=1,2,…,K,则算法收敛,计算完毕。
本发明的电力设备运行工况划分方法,可以实现:
1.本发明为无监督学习方法,能够根据数据特征自动对数据进行分析建模,不需要专业知识和经验,降低了运行人员的门槛。
2.本发明是做进一步研究运行优化的基础条件。在实际电厂应用中,每一时刻的数据对应的最优工况是不一样的,春冬和夏秋对应的外界影响因素也会有很大差距,因此,再对实时数据进行进一步分析时,首先一点先判断出其所在的最优工况是至关重要的。
3.本发明是对多指标和单指标的分析具有通用性。这更加符合实际应用的要求,对于某一设备影响工况的因素有多个时,而这时由于设备的复杂运行特性,使得专业人员没有足够的经验和知识来确定这些因素时,本发明采用的方法将解决这一棘手问题。
4.本发明为故障预警提供前提基础。电厂设备故障的预警是当前电厂最重要的一方面,而本发明通过寻求最优工况,为进一步进行预警提供先决条件,并且使得预警结果更准确,弥补预警方法的固有缺陷。
5.本发明为故障判据的标准值的确定提供一个标准。因为,从专业人员角度来说,判定当前实时值的运行情况要通过历史同工况下的运行情况来确定一个标准值,那么对于最优工况的确定将至关重要。
6.本发明是提高火电厂经济效益的重要手段。使得火电更节能更经济是当前最重要的任务之一,而最优工况的研究为大力提高电厂节能经济提供了一个更加优化的环境,使得各种方法和技术能够更加符合其使用环境,将对解决各类问题更加有意义。
附图说明
图1为建模及模型评价分析过程的步骤流程图
图2为均值聚类的流程图
图3为SIL指标最优类别流程图
图4为CH指标最优类别流程图
具体实施方式
下面详细说明本发明的具体实施,有必要在此指出的是,以下实施只是用于本发明的进一步说明,不能理解为对本发明保护范围的限制,该领域技术熟练人员根据上述本发明内容对本发明做出的一些非本质的改进和调整,仍然属于本发明的保护范围。
图1为本发明建立模型的流程图,整个建模及模型评价过程主要包括以下步骤:
步骤1:获取训练数据。
本发明选取的训练数据为一单设备的多指标历史健康数据,大致操作过程为:以选定的指标为准从电厂数据库中读取该相关设备足够时间的历史运行状态数据;然后利用设定的筛选规则从所有的历史数据中筛选出运行状态良好的健康设备数据作为构建均值聚类模型的训练数据。
假若从一个具有N个观察指标的设备数据库截取M个时刻的设备数据作为训练数据,则在m时刻的各个观察指标数据可看作一个N维的向量,可表示为:
X(tm)=[xm1,xm2,xm3,…,xmN]T (1)
则该训练数据文件应当为M*N的矩阵形式。具体形式如下:
步骤2:数据预处理。
数据预处理需要对训练数据采取删除无效数据,归一化等预处理措施。由于数据获取过程中可能存在的问题,最初获得的原始训练数据可能存在空数据等无效数据,需要将包含无效数据的观测向量进行删除。所采用的公式为:
其中,datau、datal为前后边界数据,momentl、momentu为前后边界时间,datan、momentn为NaN数据及其时间。
另外,由于设备模型相关参数的量纲不同,且不同参数数据绝对值相差很大,为保证使用非线性算子正确衡量不同观测向量之间的距离,需要对各个参数的测量值根据各自的极值进行归一化处理。可采用如(3)式所示的线性归一化方式,对各参数数据进行预处理,将各参数值全部映射到[0,1]区间内:
步骤3:确定初始聚类中心。
本发明确定初始聚类中心的算法采用K-means++算法进行优化。K-means++算法选择初始聚类中心的基本思想是:初始聚类中心之间的相互距离尽可能的远。
具体步骤如下:
(1)从输入的数据点集合中随机选择一个点作为第一个聚类中心;
(2)对于数据集中的每一个点x,计算它与最近聚类中心的距离D(X);
(3)选择一个新的数据点作为新的聚类中心,选择的原则是:D(X)较大的点,被选取作为聚类中心的概率较大;
(4)重复(2)和(3)直到K个聚类中心被选出来;
从上面算法描述上可以看到,关键是(3)中将D(X)上的点确定为被选择的概率上,本发明采用以下方法:
(1)先从我们的数据库随机挑选K个随机点当“种子点”。
(2)对于每个点,我们都计算其和最近的一个“种子点”的距离D(X)并保存在一个数组里,然后把这些距离加起来得到Sum(D(x))。
(3)然后,再取一个随机值,用权重的方式来取计算下一个“种子点”。这个算法的实现是,先取一个能落在Sum(D(x))中的随机值Random,然后用Random-=D(x),直到其D(x)<=0,此时的点就是下一个“种子点”。
(4)重复(2)和(3)直到K个聚类中心被选出来。
步骤4:确定聚类划分类别的上下限。
要想得到划分类别的上下限,从以下两方面考虑。
从专业角度,专工根据自身经验以及电力的实际需求给出一个大体的划分标准,譬如,对温度相关的设备,冬天和夏天差距就会较大。从专业角度,来说划分工况为一个,相对来说作用不大,当然也要根据实际应用需要,确定划分出的最优工况个数是符合实际应用需求的。
从算法人员角度,要通过对算法的研究,制定一个分类类别标准。比如,上限不可能超过样本总个数。但是,一般数据量较大时,我们不可能把所有的样本个数都试一遍,这样一来工作量太大。我们可以去掉重复的数据或近似重复的数据集,使得整个数据集样本大大减小。采用先使用系统聚类法得到一个相对粗略的最优值,确定出划分的上限。再用均值聚类从2至最大值进行循环,采用评价函数,选取最优工况划分的类别数。
步骤5:实现多种类别的聚类过程。
k-means算法也称C-均值算法,是根据函数准则进行分类的聚类算法,基于使聚类准则函数最小化。k-means算法的主要思想是先在需要分类的数据中寻找K组数据作为初始聚类中心,然后计算其他数据距离这三个聚类中心的距离,将数据归入与其距离最近的聚类中心,之后再对这K个聚类的数据计算均值,作为新的聚类中心,继续以上步骤,直到新的聚类中心与上一次的聚类中心值相等时结束算法。
准则函数:聚类集中每一样本点到该类中心的距离平方和。对于第j个聚类集,准则函数定义为:
Sj:第j个聚类集(域),聚类中心为Zj;
Nj:第j个聚类集Sj中所包含的样本个数。
对所有K个模式类有:
k-means算法的聚类准则:聚类中心的选择应使准则函数J极小,即使Jj的值极小。即有:
上式表明,Sj类的聚类中心应选为该类样本的均值。
具体描述:
(1)选取k-means++算法得到K个初始聚类中心:Z1(1),Z2(1),...,ZK(1);
(2)按最小距离原则将其余样品分配到K个聚类中心中的某一个,即:
min{||X-Zi(k)||,i=1,2,…,K}=||X-Zj(k)||=Dj(k) (9)
则X∈Sj(k)。
(3)计算各个聚类中心的新向量值:Zj(k+1)j=1,2,…,K
Nj:第j类的样本数。
(4)如果Zj(k+1)≠Zj(k)j=1,2,…,K,则回到步骤(2),将模式样本逐个重新分类,重复迭代计算。
(5)如果Zj(k+1)=Zj(k)j=1,2,…,K,则算法收敛,计算完毕。
步骤6:对聚类结果进行分析,确定最优K值。
聚类分析的目的就是确定一个最优的聚类类别数K。那么,聚类结果的评价至关重要。
在聚类有效性研究中,评价聚类算法得到聚类结果的方法主要由三类。
外部标准:用事先判定的聚类结构来评价;它用来度量聚类算法发现的聚类结构与某种外部结构的匹配程度。它使用了不在数据集中出现的信息,也称监督的度量,即其簇的个数及每个数据项的正确的分类均为已知。主要的方法有:熵,纯度,F-度量等。
内部标准:用参与聚类的样本来评价聚类结果,比如采用各个簇的误差平方和。如果处理的数据集结构未知,聚类结果的评价就只能依赖数据集自身的特征和量值。在这种情况下,聚类分析的地量追求两个目标:紧密度和分离度。此外,还要考虑单个簇的大小,以达到均衡较好的解。上述几种方法大多是基于统计测试,计算量比较大。主要的方法有:Calinski-Harabasz(CH)指标、Davies-Bouldin(DB)指标、Weighted inter-intra(Wint)指标、Krzanowski-Lai(KL)指标、Hartigan(Hart)指标、In-Group Proportion(IGP)指标等。
本发明中,采用了Calinski-Harabasz(CH)指标和'silhouette'(Sil)指标。
相对标准:用同一种算法的不同结果来评价聚类效果,即通过与其他结果的比较来判断聚类的优劣。相对度量是在确定聚类算法的基础上,采用预定义的评价标准,针对算法不同的参数设置进行算法测试,最终选择最优的参数设置和聚类模式。主要的方法有:簇内和簇间距离的线性组合,簇内和簇间距离的非线性组合,模糊聚类度量等。
实施例
本实施例以某火力发电厂1#机组的高温过热器为对象,一次风机是锅炉汽水系统中重要的设备,它结构复杂,难以建立精确数学机理模型,且易发多发故障人工排查难度大,符合本发明适应场景。通过本实施例的详细阐述,进一步说明本发明的实施过程。
本发明实施例对某电厂高温过热器设备的工况寻优的实施步骤如下:
步骤1:选取除停机以外的高温过热器设备的相关指标共33个。包括高过入(出)口温度(℃),高过壁温(℃),实发功率(MW),出口压力(kPa)等,因此该设备的观测向量为33维的向量:
从13年8月到14年7月半年的历史数据中去除异常部分以获取训练数据,共355438条数据(去掉空数据和停机数据)。
步骤2:对训练数据进行数据预处理工作。
对挑选出的训练数据进行删除无效数据和归一化等预处理措施。按照(1)式进行线性归一化方式,将各参数值全部映射到0到1的区间内。
步骤3:确定初始聚类中心。
利用k-means++算法得到的初始聚类中心,并且选取前1000组数据,以K=2为例,部分指标结果如下表所示:
步骤4:确定聚类划分类别的上下限。
对于高温过热器指标一般都是以温度为主,从专业角度给出了可以按照季节分工况;并且每个月份的情况可能会不一样,因此,上限最高为12个。如果划分一个工况,将对实际应用没有较大意义。因此,确定最优工况范围为[2,12]。
步骤5:实现多种类别的聚类过程。
对多类别进行聚类得到聚类结果,选取K=2,数据量为100000为例,部分结果展示如下:
步骤6:对聚类结果进行分析,确定最优K值。
由于本实施例实际应用中专工没有已知专业的知识提供,因此,适用于内部标准。本发明确定K值的过程。从以下三个角度来考虑:
(1)K值上下限:由步骤4以确定为[2,12]。
(2)数据集分别选取不同数据集进行测试,确定最优的K。按照每分钟取一条历史数据的频率取数,数据集分别选取:100,1000,10000,100000条数据。
(3)不同的评价标准;Calinski-Harabasz(CH)指标和'silhouette'(Sil)指标。
Calinski-Harabasz(CH)指标
CH指标是基于全部样本的类内离差矩阵与类间离差矩阵的测度,其最大值对应的类数作为最优的聚类个数。
其中,n表示聚类的数目,k表示当前的类,trSB(k)表示类间离差矩阵的迹,trSW(k)表示类内离差矩阵的迹。
'silhouette'(Sil)指标
设a(t)为聚类Cj中的样本t和类内所有其他样本的平均不相似度或距离,d(t,Ci)为样本t到另一个类Cj的所有仰恩的平均不相似度或距离,则b(t)=min{d(t,Ci)},i=1,2,…,k,i≠j。Sil指标计算每一个样本与同一聚类中样本的不相似度以及与其他聚类中样本的不相似度,其每个样本t的计算公式如下:
一般以一个数据集的所有样本的平均Sil值来评价聚类结果的质量,Sil指标越大表示聚类质量越好,其最大值对应的类别数作为最优的聚类个数。
通过图3、图4所示,数据量为1000时,两个指标的最优类别个数最高值为2,进而,得到最优工况K值最优结果2。通过不同数量级和不同指标的多类别对比试验,如下表所示:
通过多种评价标准从多维度分析,数量级100数据量较小。100000条大约两个月数据相对较合理,从专业应用角度分析,数据量也足够代表正常运行的状态。所以最终选取K=2作为最优的工况类别数。
步骤7:选取最优K值模型作为最优工况划分结果。
本发明最终选取的最优工况的划分类别为2,较符合专业按月划分工况的要求;各个工况的最优聚类中心的部分结果如步骤5中表格所示。
尽管为了说明的目的,已描述了本发明的示例性实施方式,但是本领域的技术人员将理解,不脱离所附权利要求中公开的发明的范围和精神的情况下,可以在形式和细节上进行各种修改、添加和替换等的改变,而所有这些改变都应属于本发明所附权利要求的保护范围,并且本发明要求保护的产品各个部门和方法中的各个步骤,可以以任意组合的形式组合在一起。因此,对本发明中所公开的实施方式的描述并非为了限制本发明的范围,而是用于描述本发明。相应地,本发明的范围不受以上实施方式的限制,而是由权利要求或其等同物进行限定。
Claims (1)
1.一种电力设备运行工况划分方法,其特征在于,依次包括如下步骤:
步骤1:从电厂数据库中读取相关设备足够时间的历史运行状态数据,利用设定的筛选规则从历史运行状态数据中筛选出运行状态良好的健康设备数据作为构建均值聚类模型的训练数据;
步骤2:对训练数据进行预处理,删除无效数据,并且进行归一化处理;
步骤3:采用K-means++算法确定初始聚类中心;
步骤4:确定聚类划分类别的上下限;
步骤5:循环实现多种类别的聚类过程;
步骤6:对聚类结果进行评价分析,确定最优的聚类类别数K值;
步骤7:选取最优K值模型作为最优工况划分结果;
步骤1还包括假若从一个具有N个观察指标的设备数据库截取M个时刻的设备数据作为训练数据,则在m时刻的各个观察指标数据xmi能够看作为一个N维的向量,其中1≤i≤N,表示为:
X(tm)=[xm1,xm2,xm3,…,xmN]T;
则该训练数据文件应当为M*N的矩阵形式, 具体形式如下:
步骤2具体步骤为将包含无效数据的观测向量进行删除,采用的公式为:
其中,datau、datal为前后边界数据,momentl、momentu为前后边界时间,datan、momentn为NaN数据及其时间;
对各个参数的测量值根据各自的极值进行归一化处理,采用如下公式对各个参数的测量值进行预处理,将各个参数的测量值全部映射到[0,1]区间内:
其中max(data)和min(data)分别表示对应测量值的最大值和最小值,data表示对应各个参数的测量值;
步骤3中确定初始聚类中心具体步骤为:
(1)从输入的数据点集合中随机选择一个点作为第一个聚类中心;
(2)对于数据集中的每一个点x,计算它与最近聚类中心的距离D(X);
(3)选择一个新的数据点作为新的聚类中心,选择的原则是:D(X)较大的点,被选取作为聚类中心的概率较大;
(4)重复(2)和(3)直到K个聚类中心被选出来;
步骤5包括如下步骤:
(1)选取k-means++算法得到K个初始聚类中心:Z1(1),Z2(1),...,ZK(1);
(2)按最小距离原则将其余样品分配到K个聚类中心中的某一个,即:
min{||X-Zi(k)||,i=1,2,…,K}=||X-Zj(k)||=Dj(k),则X∈Sj(k);
(3)计算各个聚类中心的新向量值:Zj(k+1)j=1,2,…,K;
Nj为第j类的样本数;
(4)如果Zj(k+1)≠Zj(k),j=1,2,…,K,则回到步骤(2),将模式样本逐个重新分类,重复迭代计算;
(5)如果Zj(k+1)=Zj(k),j=1,2,…,K,则算法收敛,计算完毕。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611263826.2A CN106709662B (zh) | 2016-12-30 | 2016-12-30 | 一种电力设备运行工况划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611263826.2A CN106709662B (zh) | 2016-12-30 | 2016-12-30 | 一种电力设备运行工况划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106709662A CN106709662A (zh) | 2017-05-24 |
CN106709662B true CN106709662B (zh) | 2021-07-02 |
Family
ID=58905608
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611263826.2A Active CN106709662B (zh) | 2016-12-30 | 2016-12-30 | 一种电力设备运行工况划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106709662B (zh) |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052707B (zh) * | 2017-11-28 | 2021-09-14 | 中国船舶工业系统工程研究院 | 一种基于聚类分析的船舶航行工况划分方法 |
CN107992449B (zh) * | 2017-12-05 | 2021-04-30 | 北京工业大学 | 一种基于低秩表示的地铁异常流量检测方法 |
CN108061877B (zh) * | 2017-12-14 | 2019-10-18 | 电子科技大学 | 一种基于角度信息的无源多站多目标测向交叉定位方法 |
CN108304567B (zh) * | 2018-02-09 | 2020-02-21 | 杭州安脉盛智能技术有限公司 | 高压变压器工况模式识别与数据分类方法及系统 |
CN108573285A (zh) * | 2018-04-25 | 2018-09-25 | 河海大学常州校区 | 一种面向疏浚能耗的工艺参量聚类分析方法 |
CN108596489A (zh) * | 2018-04-25 | 2018-09-28 | 河海大学常州校区 | 一种面向疏浚产量的工艺参量聚类分析方法 |
CN108595381B (zh) * | 2018-04-27 | 2022-03-22 | 厦门尚为科技股份有限公司 | 健康状态评价方法、装置及可读存储介质 |
CN109145988A (zh) * | 2018-08-22 | 2019-01-04 | 广东电网有限责任公司 | 脱硝系统的目标工况的确定方法、装置、设备及存储介质 |
CN109344488B (zh) * | 2018-09-25 | 2023-11-24 | 云南电网有限责任公司电力科学研究院 | 一种火电机组回热加热器性能状态的证据表征构建方法 |
CN109298225B (zh) * | 2018-09-29 | 2020-10-09 | 国网四川省电力公司电力科学研究院 | 一种电压量测数据异常状态自动识别模型系统及方法 |
CN109685348A (zh) * | 2018-12-14 | 2019-04-26 | 广东电网有限责任公司 | 一种电力设备典型运行模式确定方法 |
CN110133488B (zh) * | 2019-04-09 | 2021-10-08 | 上海电力学院 | 基于最优等级数的开关柜健康状态评价方法及装置 |
CN110033141B (zh) * | 2019-04-22 | 2022-10-04 | 大唐环境产业集团股份有限公司 | 一种脱硫系统运行工况数据库的建立方法 |
CN110990495B (zh) * | 2019-12-13 | 2023-05-02 | 华润电力技术研究院有限公司 | 数据处理方法、数据处理系统及计算机存储介质 |
CN111612033A (zh) * | 2020-04-15 | 2020-09-01 | 广东电网有限责任公司 | 基于引力搜索与密度峰值聚类的配电变压器故障诊断方法 |
CN111553580A (zh) * | 2020-04-22 | 2020-08-18 | 浙江农林大学 | 基于k均值和支持向量机的护林员巡护情况的综合分析方法 |
CN111598470B (zh) * | 2020-05-20 | 2023-03-24 | 贵州电网有限责任公司 | 一种配网物资市场价格采集及监测预警方法及系统 |
CN111914942A (zh) * | 2020-08-12 | 2020-11-10 | 烟台海颐软件股份有限公司 | 一种多表合一用能异常分析方法 |
CN112181003A (zh) * | 2020-08-24 | 2021-01-05 | 中国神华能源股份有限公司国华电力分公司 | 一种机组压力控制的方法和电子设备 |
CN112084095B (zh) * | 2020-09-18 | 2021-12-21 | 广域铭岛数字科技有限公司 | 基于区块链的能源网联监控方法、系统及存储介质 |
CN112306414A (zh) * | 2020-10-30 | 2021-02-02 | 北京泽石科技有限公司 | 数据处理方法及装置 |
CN113065766B (zh) * | 2021-04-01 | 2024-05-14 | 中核核电运行管理有限公司 | 一种基于历史数据挖掘分析的汽机运行工况寻优方法 |
CN113339207B (zh) * | 2021-06-29 | 2023-05-26 | 浙江理工大学 | 一种风电机组主传动系统的运行工况划分方法及系统 |
CN113569950B (zh) * | 2021-07-28 | 2024-05-28 | 大唐环境产业集团股份有限公司 | 电站设备故障监测模型生成方法、系统及装置 |
CN113781758A (zh) * | 2021-09-07 | 2021-12-10 | 浙江大学 | 面向高端燃煤发电装备的变量协同动态报警阈值优化方法 |
CN114565209A (zh) * | 2021-12-28 | 2022-05-31 | 安徽省安泰科技股份有限公司 | 一种基于聚类的流程工业能耗状态评估方法 |
CN115079654B (zh) * | 2022-07-01 | 2023-05-23 | 西安热工研究院有限公司 | 火电厂控制过程分析方法、装置、设备及介质 |
CN116304835B (zh) * | 2023-03-31 | 2023-08-29 | 北京博华信智科技股份有限公司 | 基于ai的动设备工况监测管理方法、系统、设备及介质 |
CN116894744A (zh) * | 2023-07-21 | 2023-10-17 | 国网湖北省电力有限公司随州供电公司 | 一种基于改进k均值聚类算法的电网用户数据分析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103154845A (zh) * | 2010-07-16 | 2013-06-12 | 纽约市哥伦比亚大学托管会 | 电网的机器学习 |
CN104102773A (zh) * | 2014-07-05 | 2014-10-15 | 山东鲁能软件技术有限公司 | 一种设备故障预警及状态监测方法 |
WO2015066421A1 (en) * | 2013-11-01 | 2015-05-07 | H. Lee Moffitt Cancer Center And Research Institute, Inc. | Integrated virtual patient framework |
WO2015191312A1 (en) * | 2014-06-09 | 2015-12-17 | Ultragenyx Pharmaceutical Inc. | The effective and efficient control of serum phosphate for optimal bone formation |
-
2016
- 2016-12-30 CN CN201611263826.2A patent/CN106709662B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103154845A (zh) * | 2010-07-16 | 2013-06-12 | 纽约市哥伦比亚大学托管会 | 电网的机器学习 |
WO2015066421A1 (en) * | 2013-11-01 | 2015-05-07 | H. Lee Moffitt Cancer Center And Research Institute, Inc. | Integrated virtual patient framework |
WO2015191312A1 (en) * | 2014-06-09 | 2015-12-17 | Ultragenyx Pharmaceutical Inc. | The effective and efficient control of serum phosphate for optimal bone formation |
CN104102773A (zh) * | 2014-07-05 | 2014-10-15 | 山东鲁能软件技术有限公司 | 一种设备故障预警及状态监测方法 |
Non-Patent Citations (4)
Title |
---|
kmeans聚类理论篇;bourneli(李伯韬);《博客园》;20140404;正文K值选取部分 * |
K-Means聚类算法原理;刘建平Pinard;《博客园》;20161212;正文第3节 * |
第二章(K均值算法);Heavenmoonqxy;《百度文库》;20130120;全文 * |
聚类分析中类数估计方法的实验比较;王开军等;《计算机工程》;20080505;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106709662A (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106709662B (zh) | 一种电力设备运行工况划分方法 | |
WO2022110557A1 (zh) | 一种台区户变关系异常诊断方法及装置 | |
CN109598451B (zh) | 一种基于pca与lstm神经网络的非侵入式负荷辨识方法 | |
CN110458230A (zh) | 一种基于多判据融合的配变用采数据异常甄别方法 | |
CN106779200A (zh) | 基于在历史数据中进行相似搜索的风电机组状态预测方法 | |
CN111754026B (zh) | 光伏电站群功率预测方法、装置、计算机设备及存储介质 | |
CN109670676A (zh) | 基于支持向量数据描述的配网台区风险预警方法及系统 | |
CN111878320B (zh) | 风力发电机组的监控方法、系统及计算机可读存储介质 | |
CN109086793A (zh) | 一种风力发电机的异常识别方法 | |
CN106663086A (zh) | 用于核回归模型的集体的设备和方法 | |
CN109272258B (zh) | 基于K-means聚类的区域风光发电资源评估方法 | |
CN109344518A (zh) | 一种基站热管理系统的故障诊断方法 | |
CN111126499A (zh) | 一种基于二次聚类的用电行为模式分类方法 | |
CN112213687A (zh) | 基于伪异常点辨识的关口电能表数据异常检测方法及系统 | |
CN115163424A (zh) | 基于神经网络的风电机组齿轮箱油温故障检测方法和系统 | |
CN111461565A (zh) | 一种电力调控下的电源侧发电性能评估方法 | |
CN114443338A (zh) | 面向稀疏负样本的异常检测方法、模型构建方法及装置 | |
CN116975670A (zh) | 一种短期调荷场景下线损特征分类处理方法 | |
CN111880090B (zh) | 一种面向百万千瓦超超临界机组的在线故障检测方法 | |
BR112021013153A2 (pt) | Método, sistema e produto de programa de computador para avaliação de consumo de energia em ambientes industriais | |
CN109522657A (zh) | 一种基于相关性网络和svdd的燃气轮机异常检测方法 | |
CN113989073B (zh) | 一种基于大数据挖掘的光伏高占比配网电压时空多维评估方法 | |
CN114046816B (zh) | 基于轻量级梯度提升决策树的传感器信号故障诊断方法 | |
Qiao et al. | Study on K-means method based on Data-Mining | |
Zhou et al. | Study on Optimization of Data-Driven Anomaly Detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 250101 5th floor, block B, Yinhe building, 2008 Xinluo street, high tech Zone, Jinan City, Shandong Province Patentee after: Shandong luruan Digital Technology Co.,Ltd. Address before: 250101 5th floor, block B, Yinhe building, 2008 Xinluo street, high tech Zone, Jinan City, Shandong Province Patentee before: SHANDONG LUNENG SOFTWARE TECHNOLOGY Co.,Ltd. |
|
CP01 | Change in the name or title of a patent holder |