CN109948869A - 基于有序聚类离散化的脱硫系统so2出口浓度预测方法 - Google Patents

基于有序聚类离散化的脱硫系统so2出口浓度预测方法 Download PDF

Info

Publication number
CN109948869A
CN109948869A CN201910323588.7A CN201910323588A CN109948869A CN 109948869 A CN109948869 A CN 109948869A CN 201910323588 A CN201910323588 A CN 201910323588A CN 109948869 A CN109948869 A CN 109948869A
Authority
CN
China
Prior art keywords
sample
section
exit concentration
discretization
stable state
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
Application number
CN201910323588.7A
Other languages
English (en)
Other versions
CN109948869B (zh
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.)
Southeast University
Datang Environment Industry Group Co Ltd
Original Assignee
Southeast University
Datang Environment Industry 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 Southeast University, Datang Environment Industry Group Co Ltd filed Critical Southeast University
Priority to CN201910323588.7A priority Critical patent/CN109948869B/zh
Publication of CN109948869A publication Critical patent/CN109948869A/zh
Application granted granted Critical
Publication of CN109948869B publication Critical patent/CN109948869B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明公开了一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,包括从SIS历史数据库中采集脱硫系统相关参数的历史运行数据;以负荷、SO2进口浓度为条件对历史数据进行稳态判断筛选得到原始稳态样本;将原始稳态样本的SO2出口浓度进行升序排列得到一维有序样本,然后利用有序聚类算法对其进行离散化,得到k个出口浓度区间,并统计得到各个区间的样本数量;样本数量最少的区间内的样本数量为m,对除样本数量最少的区间之外的所有区间进行抽样,使得各个出口浓度区间内的样本数量最接近m,得到SO2出口浓度均匀分布的稳态样本集;基于均匀分布的稳态样本集,采用BP神经网络算法进行模型训练,建立SO2出口浓度预测模型。

Description

基于有序聚类离散化的脱硫系统SO2出口浓度预测方法
技术领域
本发明涉及SO2出口浓度预测方法,尤其涉及一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法。
背景技术
当前电力行业烟气脱硫技术中以石灰石-石膏湿法脱硫工艺为主,用含石灰石粉的浆液为吸收剂吸收烟气中的SO2等酸性气体。脱硫系统耗电率在0.6-0.7%之间,此外还需要消耗大量的石灰石和工艺水,脱硫系统的优化运行对于电厂节能减排具有重要的意义。《中华人民共和国环境保护税法》自2018年1月1日起正式施行,环保税政策有了很大变化。对于脱硫系统而言,在完成区域内污染物排放标准的前提下,可以通过降低污染物排放浓度来达到节约环保税额的目的,脱硫系统在实际运行过程中需要密切关注出口SO2浓度的变化。控制SO2出口浓度不仅有助于实现脱离系统的达标排放,也有利于提升脱硫系统的经济性,因此SO2出口浓度的精确预测显得尤为重要。
近年来,许多学者对脱硫系统优化运行建立了多种数学模型,徐刚等人在研究脱硫系统特性的基础上通过调整运行方式对脱硫系统进行了优化,顾慧等人基于数据挖掘方法建立了脱硫系统目标工况库模型,但目前的研究都未考虑SO2出口浓度对于脱硫系统经济运行的影响。
由于深度调峰以及煤质变化的影响,脱硫系统的历史数据具有动态性、不均匀性等特征,这给脱硫系统建模以及SO2出口浓度预测带来了很大困难。
因此,亟待解决上述技术难题。
发明内容
发明目的:本发明的目的是提供一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,对数据进行离散化和抽样可保证样本分布的均匀性,为后续建模提供均匀样本,从而提升脱硫系统SO2出口浓度预测的准确性。
技术方案:为实现以上目的,本发明公开了一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,包括如下步骤:
(1)、从SIS历史数据库中采集脱硫系统相关参数的历史运行数据,取数间隔为1min;其中相关参数包括负荷、SO2进口浓度、浆液密度、pH值、液气比和SO2出口浓度;
(2)、以负荷、SO2进口浓度为条件对历史数据进行稳态判断,筛选得到原始稳态样本;
(3)、将原始稳态样本的SO2出口浓度进行升序排列得到一维有序样本,然后利用有序聚类算法对其进行离散化,得到k个出口浓度区间,并统计得到各个区间的样本数量;
(4)、样本数量最少的区间内的样本数量为m,对除样本数量最少的区间之外的所有区间进行抽样,使得各个出口浓度区间内的样本数量最接近m,得到SO2出口浓度均匀分布的稳态样本集;
(5)、基于SO2出口浓度均匀分布的稳态样本集,采用BP神经网络算法进行模型训练,建立SO2出口浓度预测模型。
其中,所述步骤(2)中稳态判断的筛选条件为:10min内机组负荷变化不超过10MW、SO2进口浓度变化不超过50mg/m3
优选的,所述步骤(3)中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其进行离散化的具体步骤包括如下:
(3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算公式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(3.2)、计算分类损失函数:
其中b(n.k)表示将n个样本分为k类的一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法为:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(3.4)、计算最小分类损失变化delta(k),当delta(k)开始小于10%时对应的k为最佳分类数kopt,delta(k)的计算方法如下:
delta(k)=[C(n,k)-C(n,k+1)]/C(n,k)×100%
(3.5)、根据kopt值确定离散化区间,具体方法为:
第kopt个区间
第(kopt-1)个区间
以此类推可以得到所有离散区间{Gk,Gk-1,···,G1}。
再者,所述步骤(4)中对于任一出口浓度区间内的样本集{Si}(1≤i≤n,n>m),从其中抽取得到最接近m个样本的具体抽取方法包括如下步骤:
(4.1)、设参数离散化区间数为k1(k1≥2,k1为整数),并初始化为2;
(4.2)、将样本集{Si}中的负荷升序排列得到一维有序样本{Pi}(1≤i≤n},将样本集{Si}中的SO2进口浓度升序排列得到一维有序样本{ci}(1≤i≤n};
(4.3)、分别采用有序聚类算法将{Pi}、{ci}离散化为k1个区间,得到k1×k1个二维单元格,将{Si}的样本分别映射到各个单元格,并统计各个单元格的样本数量;
其中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其离散化为k1个离散化区间的具体步骤包括如下:
(4.3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(4.3.2)、计算分类损失函数:
其中,b(n.k)表示将n个样本分为k类的某一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(4.3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法如下:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(4.3.4)、由离散化区间数k1求离散化区间,具体方法为:
第k1个区间
第(k1-1)区间
以此类推可以得到所有离散区间{Gk,Gk-1,···G1};
(4.4)、对于样本数为1的单元格,直接采集其样本;对于样本数大于1的单元格,计算单元格内各个样本的局部密度,采集局部密度最大的样本,若存在多个密度同时最大的样本则取第一个样本,其中局部密度采用高斯核函数进行计算:
其中:dij为样本i、j之间的欧式距离,dc为截断距离;
(4.5)、统计抽取得到的样本数量总和,若最接近m则退出,若不满足则增大k1并转至步骤4.3循环,直至抽取到最接近m个数量的数据样本。
进一步,所述步骤(5)中SO2出口浓度预测模型的建立方法包括如下步骤为:
(5.1)、对模型的输入参数、输出参数进行归一化,数据为X={x1,x2,···,xi,···,xn},则归一化后的数据X’为:
其中min(X)为采集数据X的最小值,max(X)为采集数据X的最大值;
(5.2)、设置模型参数:输入层大小为5,隐层节点数为6,输出层大小为1,激活函数用sigmoid;
(5.3)、训练神经网络模型:随机值初始化权重矩阵,偏差向量初始化为0,利用梯度下降法进行参数更新,最大训练次数设置为3000,训练要求精度设为0.001。
有益效果:与现有技术相比,本发明具有以下显著优点:
(1)、本发明首先从历史数据库中筛选得到稳态数据样本,再利用有序聚类算法对数据进行离散化,能够自动确定离散化区间,使得各区间内部之间的差异最小,无需人工参与;
(2)、针对各个区间样本数量不一致问题,本发明基于有序聚类算法对样本数量多的区间进行了抽取,使得各个区间的样本数量基本一致,从而得到出口浓度均匀分布的稳态样本集;该抽样方法可发现样本密集单元,删除冗余样本,在保证建模数据完整性的情况下可获取期望的抽样数量;
(3)、本发明采用BP神经网络算法建立了SO2出口浓度预测模型,将SO2出口浓度均匀分布的稳态样本集作为模型训练数据,可对各个区间的SO2出口浓度进行准确的预测。
附图说明
图1为本发明流程示意图;
图2为本发明实施例中最小分类损失与出口浓度分类数的趋势图;
图3为本发明实施例中样本抽取数量总和与离散区间数k1的关系示意图;
图4为本发明实施例中预测效果对比图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示,一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于,包括如下步骤:
(1)、从SIS历史数据库中采集脱硫系统相关参数的历史运行数据,取数间隔为1min;其中相关参数包括负荷、SO2进口浓度、浆液密度、pH值、液气比和SO2出口浓度;
(2)、以负荷、SO2进口浓度为条件对历史数据进行稳态判断,筛选得到原始稳态样本;其中稳态判断的筛选条件为:10min内机组负荷变化不超过10MW、SO2进口浓度变化不超过50mg/m3
(3)、将原始稳态样本的SO2出口浓度进行升序排列得到一维有序样本,然后利用有序聚类算法对其进行离散化,得到k个出口浓度区间,并统计得到各个区间的样本数量;
其中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其进行离散化的具体步骤包括如下:
(3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算公式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(3.2)、计算分类损失函数:
其中b(n.k)表示将n个样本分为k类的一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法为:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(3.4)、计算最小分类损失变化delta(k),当delta(k)开始小于10%时对应的k为最佳分类数kopt,delta(k)的计算方法如下:
delta(k)=[C(n,k)-C(n,k+1)]/C(n,k)×100%
(3.5)、根据kopt值确定离散化区间,具体方法为:
第kopt个区间
第(kopt-1)个区间
以此类推可以得到所有离散区间{Gk,Gk-1,···,G1};
(4)样本数量最少的区间内的样本数量为m,对除样本数量最少的区间之外的所有区间进行抽样,使得各个出口浓度区间内的样本数量最接近m个,得到SO2出口浓度均匀分布的稳态样本集;
对于任一出口浓度区间内的样本集{Si}(1≤i≤n,n>m),从其中抽取得到最接近m个样本的具体抽取方法包括如下步骤:
(4.1)、设参数离散化区间数为k1(k1≥2,k1为整数),并初始化为2;
(4.2)、将样本集{Si}中的负荷升序排列得到一维有序样本{Pi}(1≤i≤n}、将样本集{Si}中的SO2进口浓度升序排列得到一维有序样本{ci}(1≤i≤n};
(4.3)、分别采用有序聚类算法将{Pi}、{ci}离散化为k1个区间,得到k1×k1个二维单元格,将{Si}的样本分别映射到各个单元格,并统计各个单元格的样本数量;
其中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其离散化为k1个离散化区间的具体步骤包括如下:
(4.3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(4.3.2)、计算分类损失函数:
其中,b(n.k)表示将n个样本分为k类的某一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(4.3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法如下:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(4.3.4)、由离散化区间数k1求离散化区间,具体方法为:
第k1个区间
第(k1-1)区间
以此类推可以得到所有离散区间{Gk,Gk-1,···G1};
(4.4)、对于样本数为1的单元格,直接采集其样本;对于样本数大于1的单元格,计算单元格内各个样本的局部密度,采集局部密度最大的样本,若存在多个密度同时最大的样本则取第一个样本,其中局部密度采用高斯核函数进行计算:
其中:dij为样本i、j之间的欧式距离,dc为截断距离;
(4.5)、统计抽取得到的样本数量总和,若最接近m则退出,若不满足则增大k1并转至步骤4.3循环,直至抽取到最接近m个数量的数据样本;
(5)、基于SO2出口浓度均匀分布的稳态样本集,采用BP神经网络算法进行模型训练,建立SO2出口浓度预测模型;
其中SO2出口浓度预测模型的建立方法包括如下步骤为:
(5.1)、对模型的输入参数、输出参数进行归一化,数据为X={x1,x2,···,xi,···,xn},则归一化后的数据X’为:
其中min(X)为采集数据X的最小值,max(X)为采集数据X的最大值;
(5.2)、设置模型参数:输入层大小为5,隐层节点数为6,输出层大小为1,激活函数用sigmoid;
(5.3)、训练神经网络模型:随机值初始化权重矩阵,偏差向量初始化为0,利用梯度下降法进行参数更新,最大训练次数设置为3000,训练要求精度设为0.001。
实施例
以某660MW燃煤机组脱硫系统为例,采集SIS历史数据库中2017年10月1日到2018年10月1日的历史数据,采集间隔为1min,采集参数主要包括负荷、SO2进口浓度、浆液密度、pH值、液气比、SO2出口浓度等。以负荷、SO2进口浓度为边界条件进行稳态筛选,筛选条件为:10min内负荷变化不超过10MW,SO2进口浓度变化不超过50mg/m3。稳态筛选后得到原始稳态样本2077条。
对原始稳态样本的SO2出口浓度进行升序排列,得到其最小浓度为4.0mg/m3,最大浓度为32.9mg/m3,采用有序聚类算法对其进行离散化,如图2所示得到原始稳态样本的最小分类损失随出口浓度分类数k的趋势图,当k=5时,最小分类损失变化开始小于10%,因此确定最佳分类数为5,将出口浓度分为[4.0,8.7]、[8.8,13.0]、[13.0,17.5]、[17.6,22.9]、[23,32.9]五个区间,具体如表1所示,并统计得到各个区间的样本数量及所占的百分比。
表1原始稳态样本各出口浓度区间的分布情况
从表1可以看出原始稳态样本中各个SO2出口浓度区间的样本数量占比差距非常大,[23,32.9]区间样本占比为8.1%,而[13.0,17.5]区间的样本占比高达34%,原始稳态样本的SO2出口浓度分布极不均匀。SO2出口浓度在[13.0,17.5]、[17.6,22.9]之间的样本占比较大,多于23mg/m3的样本非常少,SO2出口浓度样本分布不均匀对于建立SO2出口浓度预测模型非常不利。
以[13.0,17.5]区间为例具体说明发明的抽样方法,该区间的样本总数为706条,需从中抽取样本数量最接近168条,如此[13.0,17.5]与[23,32.9]区间的样本数量才会接近,数据的均匀性才会较好。具体取样方法如下:先令离散化区间数k1=2,利用有序聚类算法分别对升序排列后的负荷、SO2进口浓度进行离散化,得到二维单元格,再将706个样本映射到单元格中,计算每个单元格中各样本的局部密度,抽取密度最大的样本,统计抽取到的样本数量总和,然后增加k1直至抽样数量总和最接近168。[13.0,17.5]区间抽样数量总和与离散化区间数k1的关系如图3所示,当k1=16时抽样数量总和为164,与168最接近。因而最终确定的离散化区间数k1=16,抽取得到样本数量164条。
对[4.0,8.7]、[8.8,13.0]、[17.6,22.9]等三个区间采用同样的抽样方法进行抽样,直接采集[23,32.9]区间的样本,最终得到的均匀分布样本在各个SO2出口浓度区间的分布情况如表2所示,共抽取样本847条,抽取后各个区间的样本数量基本一致,样本数量所占百分比基本持平,数据分布呈现很好的均匀性。本发明对样本多的SO2出口浓度区间进行了抽样,目的是让各个出口浓度区间的样本数量大致相等,保证样本分布的均匀性,为后续建模提供均匀样本。
表2均匀分布样本各出口浓度区间的分布情况
最后基于抽取得到的均匀分布样本集,利用BP神经网络算法建立了SO2出口浓度预测模型。本文选取了150组实际运行数据对模型进行了测试,[4.0,8.7]、[8.8,13.0]、[13.0,17.5]、[17.6,22.9]、[23,32.9]五个区间各30组测试样本。经核算,本发明所采用的预测方法的平均预测偏差为2.5mg/m3,在各个SO2出口浓度区间都有较好的预测效果。为方便对比,本发明还用2077条原始稳态样本作为训练数据建立了SO2出口浓度预测模型,其平均预测偏差为3.1mg/m3。图4给出了两种训练样本的预测效果对比,本发明的预测方法得到的均匀分布稳态样本建立的预测模型不仅平均预测偏差更低,相比原始稳态样本建立的预测模型在[23,32.9]、[4.0,8.7]两个原始样稳态样本数量较少的区间预测效果也得到了很大提升,能够满足各个区间的预测需求。
综上,发明所公开的一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法能够在各个SO2出口浓度区间都取得很好的预测精度。

Claims (5)

1.一种基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于,包括如下步骤:
(1)、从SIS历史数据库中采集脱硫系统相关参数的历史运行数据,取数间隔为1min;其中相关参数包括负荷、SO2进口浓度、浆液密度、pH值、液气比和SO2出口浓度;
(2)、以负荷、SO2进口浓度为条件对历史数据进行稳态判断,筛选得到原始稳态样本;
(3)、将原始稳态样本的SO2出口浓度进行升序排列得到一维有序样本,然后利用有序聚类算法对其进行离散化,得到k个出口浓度区间,并统计得到各个区间的样本数量;
(4)、样本数量最少的区间内的样本数量为m,对除样本数量最少的区间之外的所有区间进行抽样,使得各个出口浓度区间内的样本数量最接近m,得到SO2出口浓度均匀分布的稳态样本集;
(5)、基于SO2出口浓度均匀分布的稳态样本集,采用BP神经网络算法进行模型训练,建立SO2出口浓度预测模型。
2.根据权利要求1所述的基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于:所述步骤(2)中稳态判断的筛选条件为:10min内机组负荷变化不超过10MW、SO2进口浓度变化不超过50mg/m3
3.根据权利要求1所述的基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于:所述步骤(3)中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其进行离散化的具体步骤包括如下:
(3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算公式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(3.2)、计算分类损失函数:
其中b(n.k)表示将n个样本分为k类的一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法为:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(3.4)、计算最小分类损失变化delta(k),当delta(k)开始小于10%时对应的k为最佳分类数kopt,delta(k)的计算方法如下:
delta(k)=[C(n,k)-C(n,k+1)]/C(n,k)×100%
(3.5)、根据kopt值确定离散化区间,具体方法为:
第kopt个区间
第(kopt-1)个区间
以此类推可以得到所有离散区间{Gk,Gk-1,···,G1}。
4.根据权利要求1所述的基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于:所述步骤(4)中对于任一出口浓度区间内的样本集{Si}(1≤i≤n,n>m),从其中抽取得到最接近m个样本的具体抽取方法包括如下步骤:
(4.1)、设参数离散化区间数为k1(k1≥2,k1为整数),并初始化为2;
(4.2)、将样本集{Si}中的负荷升序排列得到一维有序样本{Pi}(1≤i≤n},将样本集{Si}中的SO2进口浓度升序排列得到一维有序样本{ci}(1≤i≤n};
(4.3)、分别采用有序聚类算法将{Pi}、{ci}离散化为k1个区间,得到k1×k1个二维单元格,将{Si}的样本分别映射到各个单元格,并统计各个单元格的样本数量;
其中对于一维有序样本{x1,x2,···xn},利用有序聚类算法对其离散化为k1个离散化区间的具体步骤包括如下:
(4.3.1)、计算离差平方和矩阵Dn×n离差平方和D(i,j)的计算式为:
其中为该类的均值向量,Dn×n的其余元素置为0;
(4.3.2)、计算分类损失函数:
其中,b(n.k)表示将n个样本分为k类的某一种分法,1=j1<j2<···jk<n=jk+1-1为分类点;
(4.3.3)、构建最小分类损失矩阵Cn×n及分类标记矩阵Jn×n,具体方法如下:
J(l,k)=jlk
式中,3≤l≤n,k≤j≤n;p(l,k)表示将l个样本分为k类的最佳分法;jlk表示p(l,k)中第k类的起始样本序号;矩阵Cn×n和矩阵Jn×n中其余元素均置为0;
(4.3.4)、由离散化区间数k1求离散化区间,具体方法为:
第k1个区间
第(k1-1)区间
以此类推可以得到所有离散区间{Gk,Gk-1,···G1};
(4.4)、对于样本数为1的单元格,直接采集其样本;对于样本数大于1的单元格,计算单元格内各个样本的局部密度,采集局部密度最大的样本,若存在多个密度同时最大的样本则取第一个样本,其中局部密度采用高斯核函数进行计算:
其中:dij为样本i、j之间的欧式距离,dc为截断距离;
(4.5)、统计抽取得到的样本数量总和,若最接近m则退出,若不满足则增大k1并转至步骤4.3循环,直至抽取到最接近m个数量的数据样本。
5.根据权利要求1所述的基于有序聚类离散化的脱硫系统SO2出口浓度预测方法,其特征在于:所述步骤(5)中SO2出口浓度预测模型的建立方法包括如下步骤为:
(5.1)、对模型的输入参数、输出参数进行归一化,数据为X={x1,x2,···,xi,···,xn},则归一化后的数据X’为:
其中min(X)为采集数据X的最小值,max(X)为采集数据X的最大值;
(5.2)、设置模型参数:输入层大小为5,隐层节点数为6,输出层大小为1,激活函数用sigmoid;
(5.3)、训练神经网络模型:随机值初始化权重矩阵,偏差向量初始化为0,利用梯度下降法进行参数更新,最大训练次数设置为3000,训练要求精度设为0.001。
CN201910323588.7A 2019-04-22 2019-04-22 基于有序聚类离散化的脱硫系统so2出口浓度预测方法 Active CN109948869B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910323588.7A CN109948869B (zh) 2019-04-22 2019-04-22 基于有序聚类离散化的脱硫系统so2出口浓度预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910323588.7A CN109948869B (zh) 2019-04-22 2019-04-22 基于有序聚类离散化的脱硫系统so2出口浓度预测方法

Publications (2)

Publication Number Publication Date
CN109948869A true CN109948869A (zh) 2019-06-28
CN109948869B CN109948869B (zh) 2020-11-13

Family

ID=67014532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910323588.7A Active CN109948869B (zh) 2019-04-22 2019-04-22 基于有序聚类离散化的脱硫系统so2出口浓度预测方法

Country Status (1)

Country Link
CN (1) CN109948869B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044668A (zh) * 2019-12-26 2020-04-21 中国华电科工集团有限公司 一种scr脱硝系统催化剂寿命在线预测方法
CN112329817A (zh) * 2020-10-12 2021-02-05 鄂尔多斯应用技术学院 一种基于神经网络的采煤机稳态样本分析方法及装置
CN112365065A (zh) * 2020-11-16 2021-02-12 大唐环境产业集团股份有限公司 一种wfgd自适应在线优化调度方法
CN113144844A (zh) * 2021-03-22 2021-07-23 国家能源集团国源电力有限公司 一种脱硫剂流量控制方法、装置及煤燃烧系统
CN113361971A (zh) * 2021-07-13 2021-09-07 浙江菲达环保科技股份有限公司 石膏脱水脱硫废水联控方法及系统
CN113468813A (zh) * 2021-07-07 2021-10-01 大唐环境产业集团股份有限公司 一种脱硫系统入口so2浓度预测方法、装置、电子设备
CN113887785A (zh) * 2021-09-13 2022-01-04 江苏智臻能源科技有限公司 基于离散化采样和神经网络的台区理论线损预测方法
CN114177747A (zh) * 2021-12-02 2022-03-15 昆岳互联环境技术(江苏)有限公司 一种基于机器学习算法的烟气脱硫二氧化硫浓度预测方法
CN114995145A (zh) * 2022-06-13 2022-09-02 华能国际电力股份有限公司德州电厂 一种基于大数据的脱硫岛智能控制方法及其控制系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573050A (zh) * 2015-01-20 2015-04-29 安徽科力信息产业有限责任公司 一种基于Canopy聚类和BIRCH层次聚类的连续属性离散化方法
CN105488586A (zh) * 2015-11-23 2016-04-13 广东电网有限责任公司电力科学研究院 基于入口so2浓度的煤粉炉超低排放系统的fgd方式确定方法
CN105912822A (zh) * 2016-05-10 2016-08-31 东南大学 海水脱硫装置脱硫效率的神经网络模拟方法
CN106251001A (zh) * 2016-07-18 2016-12-21 南京工程学院 一种基于改进模糊聚类算法的光伏功率预测方法
CN106845731A (zh) * 2017-02-20 2017-06-13 重庆邮电大学 一种基于多模型融合的潜在换机用户发现方法
CN109508832A (zh) * 2018-11-22 2019-03-22 李东峰 基于变量压缩bp神经网络的电厂so2排放软测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104573050A (zh) * 2015-01-20 2015-04-29 安徽科力信息产业有限责任公司 一种基于Canopy聚类和BIRCH层次聚类的连续属性离散化方法
CN105488586A (zh) * 2015-11-23 2016-04-13 广东电网有限责任公司电力科学研究院 基于入口so2浓度的煤粉炉超低排放系统的fgd方式确定方法
CN105912822A (zh) * 2016-05-10 2016-08-31 东南大学 海水脱硫装置脱硫效率的神经网络模拟方法
CN106251001A (zh) * 2016-07-18 2016-12-21 南京工程学院 一种基于改进模糊聚类算法的光伏功率预测方法
CN106845731A (zh) * 2017-02-20 2017-06-13 重庆邮电大学 一种基于多模型融合的潜在换机用户发现方法
CN109508832A (zh) * 2018-11-22 2019-03-22 李东峰 基于变量压缩bp神经网络的电厂so2排放软测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙栓柱 等: "样本优选BP神经网络在脱硫效率预测中的应用", 《PROCESS AUTOMATION INSTRUMENTATION》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044668A (zh) * 2019-12-26 2020-04-21 中国华电科工集团有限公司 一种scr脱硝系统催化剂寿命在线预测方法
CN111044668B (zh) * 2019-12-26 2022-12-02 中国华电科工集团有限公司 一种scr脱硝系统催化剂寿命在线预测方法
CN112329817A (zh) * 2020-10-12 2021-02-05 鄂尔多斯应用技术学院 一种基于神经网络的采煤机稳态样本分析方法及装置
CN112329817B (zh) * 2020-10-12 2024-03-15 鄂尔多斯应用技术学院 一种基于神经网络的采煤机稳态样本分析方法及装置
CN112365065A (zh) * 2020-11-16 2021-02-12 大唐环境产业集团股份有限公司 一种wfgd自适应在线优化调度方法
CN112365065B (zh) * 2020-11-16 2024-07-02 大唐环境产业集团股份有限公司 一种wfgd自适应在线优化调度方法
CN113144844A (zh) * 2021-03-22 2021-07-23 国家能源集团国源电力有限公司 一种脱硫剂流量控制方法、装置及煤燃烧系统
CN113468813A (zh) * 2021-07-07 2021-10-01 大唐环境产业集团股份有限公司 一种脱硫系统入口so2浓度预测方法、装置、电子设备
CN113468813B (zh) * 2021-07-07 2022-04-15 大唐环境产业集团股份有限公司 一种脱硫系统入口so2浓度预测方法、装置、电子设备
CN113361971B (zh) * 2021-07-13 2023-04-07 浙江菲达环保科技股份有限公司 石膏脱水脱硫废水联控方法及系统
CN113361971A (zh) * 2021-07-13 2021-09-07 浙江菲达环保科技股份有限公司 石膏脱水脱硫废水联控方法及系统
CN113887785A (zh) * 2021-09-13 2022-01-04 江苏智臻能源科技有限公司 基于离散化采样和神经网络的台区理论线损预测方法
CN114177747A (zh) * 2021-12-02 2022-03-15 昆岳互联环境技术(江苏)有限公司 一种基于机器学习算法的烟气脱硫二氧化硫浓度预测方法
CN114995145A (zh) * 2022-06-13 2022-09-02 华能国际电力股份有限公司德州电厂 一种基于大数据的脱硫岛智能控制方法及其控制系统

Also Published As

Publication number Publication date
CN109948869B (zh) 2020-11-13

Similar Documents

Publication Publication Date Title
CN109948869A (zh) 基于有序聚类离散化的脱硫系统so2出口浓度预测方法
CN110263988B (zh) 一种基于电厂脱硫系统的数据运行优化方法
CN101498459B (zh) 一种锅炉燃烧优化的建模方法
CN108051035B (zh) 基于门控循环单元的神经网络模型的管网漏损识别方法
CN100561509C (zh) 一种改进的混合遗传算法优化水质模型参数的设计方法
CN112215461B (zh) 长距离引调水工程输水调度状态分析及预警系统
Wu et al. A DEA-based improvement of China's green development from the perspective of resource reallocation
CN112149314B (zh) 一种基于虚拟库容修正的多沙水库库容冲淤模拟方法
Wang et al. Evolution of coordination degree of eco-economic system and early-warning in the Yangtze River Delta
CN109343367A (zh) 一种基于神经网络预测控制烟气脱硫的方法
CN114077970B (zh) 一种基于城市形态的碳排放相关因子的确定方法及装置
CN112922582B (zh) 基于高斯过程回归的气井井口油嘴气体流量分析预测方法
CN102402726A (zh) 一种基于区域负荷分析的大规模配电网电力电量预测方法
CN111080157B (zh) 一种梯级水电站排磷量的调度方法和系统
CN112465239B (zh) 基于改进pso-fcm算法的脱硫系统运行优化方法
CN104636834B (zh) 一种改进的联合概率规划模型系统优化方法
CN106169101A (zh) 一种基于结构分解的能源需求预测方法
CN109241665A (zh) 一种叠梁门水温改善效果的分析方法
CN106600055A (zh) 一种基于自激励门限自回归模型的风速预测方法
CN102184450A (zh) 一种多孔介质燃烧器燃烧优化的建模方法
CN110348600A (zh) 一种水库群防洪效益最大的库容分配方法
CN103218516A (zh) 一种聚类再回归的机场噪声单步预测方法
CN114492617B (zh) 一种基于聚类的管网分区及跨区水量调配的方法
CN106383970B (zh) 一种基于泵装置的轴流泵必需汽蚀余量的预测方法
CN106126770A (zh) 基于s形模型的钢铁需求预测方法

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