CN111178594A - 一种火电机组调峰能力预测方法、装置及系统 - Google Patents

一种火电机组调峰能力预测方法、装置及系统 Download PDF

Info

Publication number
CN111178594A
CN111178594A CN201911274961.0A CN201911274961A CN111178594A CN 111178594 A CN111178594 A CN 111178594A CN 201911274961 A CN201911274961 A CN 201911274961A CN 111178594 A CN111178594 A CN 111178594A
Authority
CN
China
Prior art keywords
thermal power
generating unit
power generating
parameters
clustering
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.)
Pending
Application number
CN201911274961.0A
Other languages
English (en)
Inventor
刘永江
郭瑞君
张叔禹
纪煜
张国斌
韩义
韩俊飞
李晓波
禾志强
周磊
张艳飞
殷建华
张伟
王银河
张成煜
刘文哲
任资龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baotou No3 Thermal Power Plant Of North United Power Co Ltd
Hunan Huayin Energy Technology Co ltd
Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd
Original Assignee
Baotou No3 Thermal Power Plant Of North United Power Co Ltd
Hunan Huayin Energy Technology Co ltd
Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Baotou No3 Thermal Power Plant Of North United Power Co Ltd, Hunan Huayin Energy Technology Co ltd, Inner Mongolia Electric Power Research Institute of Inner Mongolia Power Group Co Ltd filed Critical Baotou No3 Thermal Power Plant Of North United Power Co Ltd
Priority to CN201911274961.0A priority Critical patent/CN111178594A/zh
Publication of CN111178594A publication Critical patent/CN111178594A/zh
Pending legal-status Critical Current

Links

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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G06F18/232—Non-hierarchical techniques
    • G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G06N3/084—Backpropagation, e.g. using gradient descent
    • 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

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Marketing (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Probability & Statistics with Applications (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种火电机组调峰能力预测方法、装置及系统，方法包括以下步骤：从火电机组运行的历史数据中筛选出数据集；根据BIC值确定最优聚类个数；基于K‑means算法对数据集进行一次聚类；基于GMM算法对数据集进行二次聚类；将采集的实时数据在所有类簇中确定相似度最高的类簇，将相似度最高的类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值作为预测结果。本方法分别通过BIC值、K‑means算法确定样本集的最优聚类个数、GMM算法的初始参数，有效提升利用EM算法求解高斯分布参数的准确性，从而提升最终预测的准确性，本发明提供的火电机组调峰能力预测方法、装置及系统能够对火电机组调峰能力进行在线预测，能为电网编制负荷调度策略提供数据基础。

Description

一种火电机组调峰能力预测方法、装置及系统
技术领域
本发明涉及火电厂发电技术领域，特别涉及一种火电机组调峰能力预测方法、装置及系统。
背景技术
火力发电是一个极为复杂的生产过程，在不同煤质和季节下，机组带负荷能力存在显著差异，要实现火电机组调峰能力监测和评估，采用传统方法很难取得较好效果，因此，如何对海量运行历史数据进行挖掘，自动划分机组运行工况，并预测出每个工况下的机组调峰能力是电网安全可靠调度、运行以及智能电网建设的要求，也是目前仍待研究的问题。
发明内容
本发明的目的在于至少解决现有技术中存在的技术问题之一，提供一种火电机组调峰能力预测方法、装置及系统，对火电机组调峰能力进行在线预测，能为电网编制负荷调度策略提供数据基础。
本发明提供了一种火电机组调峰能力预测方法、装置及系统。
根据本发明实施例的一种火电机组调峰能力预测方法，包括以下步骤：
从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，所述训练样本中包含工况参数和寻优参数；
根据BIC值确定所述数据集的最优聚类个数；
基于K-means算法按照所述最优聚类个数对所述数据集进行一次聚类，得到一次聚类结果；
基于GMM算法并以所述最优聚类个数以及所述一次聚类结果作为输入参数对所述数据集进行二次聚类，得到二次聚类结果，并且在所述二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
将采集的实时数据在所有类簇中确定相似度最高的类簇，将所述相似度最高的类簇中寻优参数的最大值、以及所述所有类簇中寻优参数的最小值作为预测结果。
根据本发明的一些实施例，所述工况参数包括：入炉煤质参数、循环水入口温度和磨煤机给煤量，所述入炉煤质参数包括：收到基水分、收到基灰分、干燥无灰基挥发分以及低位发热量。
根据本发明的一些实施例，所述将采集的实时数据在所有类簇中确定相似度最高的类簇之前，还包括步骤：
基于软测量方法对所述入炉煤质参数进行实时预测。
根据本发明的一些实施例，所述软测量方法为三层BP神经网络算法。
根据本发明的一些实施例，所述在所有类簇中确定相似度最高的类簇，进一步包括：
设定最优聚类个数为K，将采集的实时数据与所述二次聚类结果中的K个高斯模型期望值进行相似度计算并选取相似度最高的类簇。
根据本发明实施例的一种火电机组调峰能力预测装置，包括：
样本数据获取组件、最优聚类个数确定组件、多重聚类计算组件、以及预测结果确定组件；
所述样本数据获取组件用于从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，所述训练样本中包含工况参数和寻优参数；
所述最优聚类个数确定组件用于根据BIC值确定所述数据集的最优聚类个数；
所述多重聚类计算组件用于基于K-means算法按照所述最优聚类个数对所述数据集进行一次聚类，得到一次聚类结果；
所述多重聚类计算组件还用于基于GMM算法并以所述最优聚类个数以及所述一次聚类结果作为输入参数对所述数据集进行二次聚类，得到二次聚类结果，并且在所述二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
所述预测结果确定组件用于将采集的实时数据在所有类簇中确定相似度最高的类簇，将所述相似度最高的类簇中寻优参数的最大值、以及所述所有类簇中寻优参数的最小值作为预测结果。
根据本发明的一些实施例，还包括实时预测组件；
所述实时预测组件用于基于软测量方法对所述工况参数中的入炉煤质参数进行实时预测。
根据本发明实施例的一种火电机组调峰能力预测系统，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述所述的一种火电机组调峰能力预测方法。
相对于现有技术，本公开实施例披露的技术方案具有以下优点：
本发明方法分别通过BIC值、以及K-means算法确定训练样本集的最优聚类个数、以及GMM算法的初始参数，能够有效提升GMM算法中利用EM算法求解高斯分布参数的准确性，从而提升最终预测的准确性；本发明提供的一种火电机组调峰能力预测方法、装置及系统能够对火电机组调峰能力进行在线预测，能为电网编制负荷调度策略提供数据基础。
本公开实施例所实现的更多特点和优势，将在具体实施方式或实践中给出。
附图说明
下面结合附图和实施例对本发明进一步地说明；
图1为本发明实施例提供的一种火电机组调峰能力预测方法的流程示意图；
图2为本发明实施例提供的不同聚类个数对应的BIC值的示意图；
图3为本发明实施例提供的火电机机组调峰能力在线预测结果示意图；
图4为本发明实施例提供的低位发热量的实时预测结果示意图；
图5为本发明实施例提供的收到基水分的实时预测结果示意图；
图6为本发明实施例提供的干燥无灰基挥发分的实时预测结果示意图；
图7为本发明实施例提供的收到基灰分的实时预测结果示意图；
图8为本发明实施例提供的一种火电机组调峰能力预测装置的结构示意图；
图9为本发明实施例提供的一种火电机组调峰能力预测方法的流程示意图。
具体实施方式
下面将结合附图，对本公开实施例的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本公开一部分实施例，而不是全部的实施例。基于本公开的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。需要说明的是，在不冲突的情况下，本公开实施例及实施例中的特征可以相互组合。另外，附图的作用在于用图形补充说明书文字部分的描述，使人能够直观地、形象地理解本公开的每个技术特征和整体技术方案，但其不能理解为对本公开保护范围的限制。
在本发明的描述中，若干的含义是一个或者多个，多个的含义是两个以上，大于、小于、超过等理解为不包括本数，以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
参照图1，本发明的一个实施例，提供了一种火电机组调峰能力预测方法，包括以下步骤：
步骤S100、从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，训练样本中包含工况参数和寻优参数；
作为一种可实施方式，获取火电厂机组运行的历史数据，筛选正常运行历史数据得到正常状态数据集X＝{xi∈RD,yi|i＝1,2,…,n}作为训练样本，其中，xi为第i个样本中的工况参数，D为工况参数的个数，yi为第i个样本中寻优参数(机组负荷)，n为数据集中样本的个数；这里的工况参数包括：入炉煤质参数、循环水入口温度、以及磨煤机给煤量，其中，入炉煤质参数包括：收到基水分、收到基灰分、干燥无灰基挥发分、低位发热量；这里的寻优参数即机组负荷。
步骤S200、根据BIC值确定数据集的最优聚类个数；
在本步骤中，根据最佳聚类数的搜索范围，初步确定聚类个数k的范围为：
Figure BDA0002315303250000051
这里，以贝叶斯信息准则(Bayesian Information Criterion，BIC)作为判断聚类个数的评价标准的依据，确定BIC值最大的聚类个数为最优聚类个数K，式如下：
BIC＝2 lk(X,Θ)-Dlog(n) (1)
在式(1)中，lk(X,Θ)表示k个聚类时训练样本的似然函数，D表示工况参数个数，n表示训练样本个数。
步骤S300、根据最优聚类个数，基于K-means算法对数据集进行一次聚类，得到一次聚类结果；
详细步骤如下：
步骤S310、确定正常状态数据集X中工况参数数值范围，在参数范围内任选K个工况参数数据点，作为初始聚类中心点
Figure BDA0002315303250000052
步骤S320、计算数据集X中每一组样本xi到每一聚类中心
Figure BDA0002315303250000053
的欧式距离，选取距离它最近的聚类中心
Figure BDA0002315303250000054
划分到聚类中心
Figure BDA0002315303250000055
所属类簇
Figure BDA0002315303250000056
其中
Figure BDA0002315303250000057
步骤S330、采取求均值的方法重新计算分类后的各聚类中心
Figure BDA0002315303250000058
步骤S340、计算距离函数
Figure BDA0002315303250000059
若E(X,U0)收敛，则输出最终的聚类中心
Figure BDA00023153032500000510
和K个类簇
Figure BDA00023153032500000511
若E(X,U0)不收敛，则跳转至步骤S320；
步骤S350、确定每个类簇的
Figure BDA00023153032500000512
和
Figure BDA00023153032500000513
作为下一步高斯混合模型的初始参数，计算如下；
Figure BDA00023153032500000514
Figure BDA00023153032500000515
Figure BDA00023153032500000516
Figure BDA00023153032500000517
上述式(2)、(3)、以及(5)中
Figure BDA00023153032500000518
和
Figure BDA00023153032500000519
分别表示第j个类簇的权值、期望和协方差(即一次聚类结果)。
步骤S400、根据一次聚类结果以及最优聚类个数，基于GMM算法对数据集进行二次聚类，得到二次聚类结果，并且在二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
本步骤是将上述步骤S200得到的最优聚类个数以及步骤S300得到的聚类结果作为输入参数，根据GMM算法进行聚类，最后确定收敛后的每个聚类中寻优参数的最大值、以及所有聚类中寻优参数的最小值；
详细步骤如下：
首先，根据上述步骤S200中的BIC确定的最优聚类个数K和K-means算法确定的初始输入参数
Figure BDA0002315303250000061
利用式(8)、(9)、(10)、以及(11)不断迭代，直到前后两次参数的差值小于设定好的阈值，得到最终GMM(高斯混合模型)的参数Θ。
其中，GMM的详细描述如下：
概率密度函数可以用高斯混合模型表示为：
Figure BDA0002315303250000062
在式(6)中，K为高斯模型的最优聚类个数，ωk为第k个高斯分量的权重，μk、Σk分别为局部高斯模型的均值和协方差，g(xi|μk,Σk)为第k个高斯分量的多元高斯密度函数，可以表示为：
Figure BDA0002315303250000063
然后，通过期望最大化(EM)算法对GMM方法中的未知参数进行求解，即在p(X；μ,Σ)为最大值条件下(样本点xi已经发生，故可认为p(X；μ,Σ)是样本xi发生的最大概率)，求得μ和Σ。
根据训练样本X＝{xi∈RD,yi|i＝1,2,…,n}、最优聚类个数K和由K-means聚类算法确定的初始值：
Figure BDA0002315303250000064
利用EM算法通过不断重复E-step和M-step来更新参数，以保证训练数据似然度单调增加到一定值。
其中，EM算法的迭代步骤如下：
E-step：
Figure BDA0002315303250000071
式(8)中，p(s)(Ck|xi)为第s次迭代后第i个训练样本xi属于第k个高斯分量的后验概率。
M-step：
Figure BDA0002315303250000072
Figure BDA0002315303250000073
Figure BDA0002315303250000074
在式(9)、(10)、(11)中，
Figure BDA0002315303250000075
和
Figure BDA0002315303250000076
分别为第(s+1)次迭代后第k个高斯分量的期望、协方差和先验概率。
在得到高斯混合模型的数学求解结果后，基于EM算法不断求解迭代可以得到各个模型参数，从而进一步实现了对训练样本X的K个工况划分。
最后，根据越大越优的寻优目标，确定每个类簇中寻优参数的最大值maxk；根据越小越优的寻优目标，确定所有类簇中寻优参数的最小值min。
步骤S500、将采集的实时数据在所有类簇中确定相似度最高的类簇，将相似度最高的类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值作为预测结果。
详细步骤如下：
计算一组实时数据与K个类簇的相似系数sim(μk|xi)，具体式如下：
Figure BDA0002315303250000077
j＝argmax(sim(μk|xi)) (13)
xi∈μj (14)
在式(12)中，xi为现场采集的一组实时数据(即实时工况参数)，分别与高斯混合模型中的K个高斯模型期望值μk进行相似度计算，相似度最高的类簇作为实时数据xi所属类簇j。
该类簇j下的寻优参数的最大值maxj和所有类簇中寻优参数的最小值min作为采集实时数据的预测结果。
本实施例的预测结果如图2与图3所示，在图2中，BIC值最小为最优聚类个数，从图中可以看出聚类个数为20的时候为最优聚类个数(聚类个数即运行工况)；在图3中，预测了20个运行工况带负荷能力的最小值和最大值。
由于在本领域中，通过EM算法求解高斯分布参数对初始参数以及聚类个数的要求是很高的，不同的初始参数以及聚类个数会输出不同的估计结果，本实施例提供的方法通过算法手段让输入参数尽可能贴合实际情况和业务要求；本实施例通过BIC值、以及K-means算法确定训练样本集的最优聚类个数、以及GMM算法的初始参数，能够显著提高EM算法的收敛速度，从而提高EM算法求解高斯分布参数的准确性，从而提升最终预测的准确性；本实施例方法能够对火电机组调峰能力进行在线预测，能为电网编制负荷调度策略提供数据基础。
优选的，在步骤S500之前还包括步骤：基于软测量方法对入炉煤质参数进行实时预测。
入炉煤煤质参数包括：收到收到基水分、收到基灰分、干燥无灰基挥发分、低位发热量，这些参数在本领域中难以直接测量，而且又是火电机组运行工况划分的必须参数，因此需要对这些参数进行预测。
为了方便理解，对以下术语进行解释：
软测量技术(又叫做软仪表技术)，是依据较易测得的辅助过程变量(也称自变量)与难以直接测量的待测主导变量(也称因变量)之间的相关关系，通过数学或者其他方法计算获得待测变量的技术。软测量的核心是寻找辅助过程变量和待测主导变量之间的相关关系并建立相应数学模型。建模的方法非常多，比如回归分析、模糊数学、模式识别、人工神经网络，其中，神经网络具有很好的泛化能力，容错性较高。
为了提高计算结果精确度以及计算效率，本步骤使用BP神经网络算法对入炉煤质参数进行预测，即对收到基水分、收到基灰分、干燥无灰基挥发分和低位发热量进行实时预测。
作为一种可实施方式，本实施例采用3层BP网络结构，选取的输出层参数、输入层参数和隐含层数目如下表1所示：
Figure BDA0002315303250000091
表1
例如，从某厂600MW火电机组运行数据中共选取1873份样本，选取1700份样本进行模型训练，173份样本进行预测，结果如图4至图7所示，图中为实际值与预测值的趋势对比，由图可见，相对误差在5％以内，由此验证了实施例提供的模型、方法的实用性。
在本实施例中，通过BP神经网络对工况参数中的收到基水分、收到基灰分、干燥无灰基挥发分和低位发热量进行在线预测，有效解决了目前入炉煤质化验结果严重滞后的问题，达到了实时预测的要求。
参照图8，本发明的另一个实施例，提供了一种火电机组调峰能力预测装置，包括：
样本数据获取组件、最优聚类个数确定组件、多重聚类计算组件、以及预测结果确定组件；
样本数据获取组件用于从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，训练样本中包含工况参数和寻优参数；
最优聚类个数确定组件用于根据BIC值确定数据集的最优聚类个数；
多重聚类计算组件用于根据最优聚类个数，基于K-means算法对数据集进行一次聚类，得到一次聚类结果；
多重聚类计算组件还用于根据一次聚类结果以及最优聚类个数，基于GMM算法对数据集进行二次聚类，得到二次聚类结果，并且在二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
预测结果确定组件用于将采集的实时数据在所有类簇中确定相似度最高的类簇，将相似度最高的类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值作为预测结果。
优选的，还包括实时预测组件，实时预测组件用于基于软测量方法对工况参数中的入炉煤质参数进行实时预测。
可以理解的是，由于本实施例中的一种火电机组调峰能力预测装置与上述一种火电机组调峰能力预测方法基于相同的发明构思，因此，方法实施例中的相应内容同样适用于本装置实施例，此处不再详述。
本发明的另一个实施例，提供了一种火电机组调峰能力预测系统，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现如上述实施例中的一种火电机组调峰能力预测方法。
可以理解的是，由于本实施例中的一种火电机组调峰能力预测系统与上述一种火电机组调峰能力预测方法、以及一种火电机组调峰能力预测装置基于相同的发明构思，因此，方法实施例、以及装置实施例中的相应内容同样适用于本系统实施例，此处不再详述。
参照图9，本发明的另一个实施例，提供了一种火电机组调峰能力预测方法，包括以下步骤：
步骤1：给定训练样本X；
步骤2：提取工况参数和寻优参数；
步骤3：根据BIC值确定最优聚类个数；
步骤4：采用K-means算法进行聚类，直至收敛，输出聚类结果权值、期望和协方差；
步骤5：采用GMM算法进行进一步聚类，直至收敛；
步骤6：存储每个类簇工况参数的算数平均值、寻优参数的最大值和所有类簇寻优参数最小值；
步骤7：实时采集BP神经网络输入参数；
步骤8：采用BP神经网络算法对入炉煤质参数进行实时预测；
步骤9：实时采集循环水入口温度和磨煤机给煤量；
步骤10：结合入炉煤质参数和循环水入口温度和磨煤机给煤量形成实时工况参数；
步骤11：根据实时工况参数与各类簇工况参数算数平均值计算相似度，从所有类簇中找出相似性最高的类簇；
步骤12：把该类簇下寻优参数的最大值和所有类簇寻优参数最小值作为预测结果。
在本实施例中，步骤1至步骤6的具体过程等同上述方法实施例中的步骤S100-步骤S400，此处不再详述；步骤7和步骤8是基于软测量方法对入炉煤质参数进行实时预测，这里使用BP神经网络算法对入炉煤质参数进行预测，以解决本领域中入炉煤质参数难以直接测量的缺陷；步骤9和步骤10是采集循环水入口温度和磨煤机给煤量并结合进行预测后的入炉煤质参数形成实时输入数据；步骤11和步骤12等同上述方法实施例中的步骤S500，此处不再详述。
本实施例提供的一种火电机组调峰能力预测方法，通过BIC值、以及K-means算法确定训练样本集的最优聚类个数、以及GMM算法的初始参数，能够显著提高EM算法的收敛速度，从而提高EM算法求解高斯分布参数的准确性，从而提升最终预测的准确性；
通过BP神经网络对工况参数中的收到基水分、收到基灰分、干燥无灰基挥发分和低位发热量进行在线预测，有效解决了目前入炉煤质化验结果严重滞后的问题，达到了实时预测的要求；
本实施例对火电机组调峰能力进行在线预测，能为电网编制负荷调度策略提供数据基础。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
通过以上的实施方式的描述，本领域技术人员可以清楚地了解到各实施方式可借助软件加通用硬件平台的方式来实现。本领域技术人员可以理解实现上述实施例方法中的全部或部分流程是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(ReadOnly Memory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。
上面结合附图对本发明实施例作了详细说明，但是本发明不限于上述实施例，在所述技术领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (8)

1.一种火电机组调峰能力预测方法，其特征在于，包括以下步骤：
从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，所述训练样本中包含工况参数和寻优参数；
根据BIC值确定所述数据集的最优聚类个数；
根据所述最优聚类个数，基于K-means算法对所述数据集进行一次聚类，得到一次聚类结果；
根据所述一次聚类结果以及所述最优聚类个数，基于GMM算法对所述数据集进行二次聚类，得到二次聚类结果，并且在所述二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
将采集的实时数据在所有类簇中确定相似度最高的类簇，将所述相似度最高的类簇中寻优参数的最大值、以及所述所有类簇中寻优参数的最小值作为预测结果。
2.根据权利要求1所述的一种火电机组调峰能力预测方法，其特征在于，所述工况参数包括：入炉煤质参数、循环水入口温度和磨煤机给煤量，所述入炉煤质参数包括：收到基水分、收到基灰分、干燥无灰基挥发分以及低位发热量。
3.根据权利要求2所述的一种火电机组调峰能力预测方法，其特征在于，所述将采集的实时数据在所有类簇中确定相似度最高的类簇之前，还包括步骤：
基于软测量方法对所述入炉煤质参数进行实时预测。
4.根据权利要求3所述的一种火电机组调峰能力预测方法，其特征在于：所述软测量方法为三层BP神经网络算法。
5.根据权利要求1所述的一种火电机组调峰能力预测方法，其特征在于，所述在所有类簇中确定相似度最高的类簇，进一步包括：
设定最优聚类个数为K，将采集的实时数据与所述二次聚类结果中的K个高斯模型期望值进行相似度计算并选取相似度最高的类簇。
6.一种火电机组调峰能力预测装置，其特征在于，包括：
样本数据获取组件、最优聚类个数确定组件、多重聚类计算组件、以及预测结果确定组件；
所述样本数据获取组件用于从火电机组正常运行的历史数据中筛选出包含若干个训练样本的数据集，所述训练样本中包含工况参数和寻优参数；
所述最优聚类个数确定组件用于根据BIC值确定所述数据集的最优聚类个数；
所述多重聚类计算组件用于根据所述最优聚类个数，基于K-means算法对所述数据集进行一次聚类，得到一次聚类结果；
所述多重聚类计算组件还用于根据所述一次聚类结果以及所述最优聚类个数，基于GMM算法对所述数据集进行二次聚类，得到二次聚类结果，并且在所述二次聚类结果中确定每个类簇中寻优参数的最大值、以及所有类簇中寻优参数的最小值；
所述预测结果确定组件用于将采集的实时数据在所有类簇中确定相似度最高的类簇，将所述相似度最高的类簇中寻优参数的最大值、以及所述所有类簇中寻优参数的最小值作为预测结果。
7.根据权利要求6所述的一种火电机组调峰能力预测装置，其特征在于，还包括实时预测组件；
所述实时预测组件用于基于软测量方法对所述工况参数中的入炉煤质参数进行实时预测。
8.一种火电机组调峰能力预测系统，包括：存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5中任一所述的一种火电机组调峰能力预测方法。
CN201911274961.0A 2019-12-12 2019-12-12 一种火电机组调峰能力预测方法、装置及系统 Pending CN111178594A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911274961.0A CN111178594A (zh) 2019-12-12 2019-12-12 一种火电机组调峰能力预测方法、装置及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911274961.0A CN111178594A (zh) 2019-12-12 2019-12-12 一种火电机组调峰能力预测方法、装置及系统

Publications (1)

Publication Number Publication Date
CN111178594A true CN111178594A (zh) 2020-05-19

Family

ID=70651959

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911274961.0A Pending CN111178594A (zh) 2019-12-12 2019-12-12 一种火电机组调峰能力预测方法、装置及系统

Country Status (1)

Country Link
CN (1) CN111178594A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967733A (zh) * 2020-07-29 2020-11-20 国网甘肃省电力公司电力科学研究院 一种聚合群电源调峰潜力模糊综合评估方法
CN112181003A (zh) * 2020-08-24 2021-01-05 中国神华能源股份有限公司国华电力分公司 一种机组压力控制的方法和电子设备
CN112555084A (zh) * 2020-12-04 2021-03-26 湖南五凌电力科技有限公司 一种实时状态预警方法、装置、电子设备及存储介质
CN112749508A (zh) * 2020-12-29 2021-05-04 浙江天行健智能科技有限公司 一种基于gmm和bp神经网络的路感模拟方法
CN116681184A (zh) * 2023-07-28 2023-09-01 厦门农芯数字科技有限公司 生猪生长指标的预测方法、装置、设备及可读存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150317589A1 (en) * 2012-11-09 2015-11-05 The Trustees Of Columbia University In The City Of New York Forecasting system using machine learning and ensemble methods
CN108549220A (zh) * 2018-03-29 2018-09-18 广东电网有限责任公司电力调度控制中心 燃煤机组运行状态实时评估方法及其系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150317589A1 (en) * 2012-11-09 2015-11-05 The Trustees Of Columbia University In The City Of New York Forecasting system using machine learning and ensemble methods
CN108549220A (zh) * 2018-03-29 2018-09-18 广东电网有限责任公司电力调度控制中心 燃煤机组运行状态实时评估方法及其系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
巨林仓;李磊;赵强;: "基于遗传神经网络的锅炉入炉煤质软测量研究" *
龙建平;江平;丁伟;: "基于设备健康状态的燃煤机组发电能力评估系统设计" *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967733A (zh) * 2020-07-29 2020-11-20 国网甘肃省电力公司电力科学研究院 一种聚合群电源调峰潜力模糊综合评估方法
CN112181003A (zh) * 2020-08-24 2021-01-05 中国神华能源股份有限公司国华电力分公司 一种机组压力控制的方法和电子设备
CN112555084A (zh) * 2020-12-04 2021-03-26 湖南五凌电力科技有限公司 一种实时状态预警方法、装置、电子设备及存储介质
CN112749508A (zh) * 2020-12-29 2021-05-04 浙江天行健智能科技有限公司 一种基于gmm和bp神经网络的路感模拟方法
CN112749508B (zh) * 2020-12-29 2024-03-05 浙江天行健智能科技有限公司 一种基于gmm和bp神经网络的路感模拟方法
CN116681184A (zh) * 2023-07-28 2023-09-01 厦门农芯数字科技有限公司 生猪生长指标的预测方法、装置、设备及可读存储介质
CN116681184B (zh) * 2023-07-28 2023-10-10 厦门农芯数字科技有限公司 生猪生长指标的预测方法、装置、设备及可读存储介质

Similar Documents

Publication Publication Date Title
CN111178594A (zh) 一种火电机组调峰能力预测方法、装置及系统
CN111148118B (zh) 基于时间序列的流量预测和载波关断方法和系统
CN110417011B (zh) 一种基于互信息与迭代随机森林的在线动态安全评估方法
CN112633317A (zh) 基于注意力机制的cnn-lstm风机故障预测方法及系统
CN111105104A (zh) 一种基于相似日和rbf神经网络的短期电力负荷预测方法
CN113256066A (zh) 基于PCA-XGBoost-IRF的作业车间实时调度方法
CN111523778A (zh) 基于粒子群算法和梯度提升树的电网运行安全评估方法
CN111144552B (zh) 一种粮食品质多指标预测方法及装置
CN111190349A (zh) 船舶机舱设备状态监测及故障诊断方法、系统及介质
CN108510180B (zh) 一种生产设备所处性能区间的计算方法
CN111028100A (zh) 考虑气象因素的精细化短期负荷预测方法、装置及介质
CN112861459B (zh) 全灵敏度对抗重要性抽样产量优化方法和装置
CN110986407A (zh) 一种离心式冷水机组故障诊断方法
CN108830417A (zh) 一种基于arma和回归分析的生活能源消费预测方法及系统
CN115526258A (zh) 基于Spearman相关系数特征提取的电力系统暂稳评估方法
CN115822887A (zh) 风电机组的性能评估与能效诊断方法及系统
CN112186761A (zh) 一种基于概率分布的风电功率场景生成方法及系统
CN115329669A (zh) 基于深度贝叶斯主动学习的电力系统暂态稳定评估方法
CN112990627A (zh) 电能质量评估方法
CN109116300B (zh) 一种基于非充分指纹信息的极限学习定位方法
TAŞABAT et al. Using long-short term memory networks with genetic algorithm to predict engine condition
CN117609818A (zh) 基于聚类与信息熵的电网关联关系发现方法
CN116415732A (zh) 一种基于改进arnn的用户侧电力负荷数据处理方法
CN114970353A (zh) 基于缺失数据填充的mswi过程二噁英排放软测量方法
CN110991743B (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
CB02 Change of applicant information

Address after: 010020 No.21 Xilin South Road, Hohhot, Inner Mongolia Autonomous Region

Applicant after: Inner Mongolia Electric Power Research Institute Branch of Inner Mongolia Electric Power (Group) Co.,Ltd.

Applicant after: Hunan Datang Xianyi Technology Co.,Ltd.

Applicant after: Baotou No.3 thermal power plant of North United Power Co.,Ltd.

Address before: 1001-1061, building 1, chuanggu Industrial Park, 568 queyuan Road, Tianxin District, Changsha, Hunan 410000

Applicant before: Hunan Datang Xianyi Technology Co.,Ltd.

Applicant before: Inner Mongolia Electric Power Research Institute Branch of Inner Mongolia Electric Power (Group) Co.,Ltd.

Applicant before: Baotou No.3 thermal power plant of North United Power Co.,Ltd.

CB02 Change of applicant information
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200519

WD01 Invention patent application deemed withdrawn after publication