CN113504765A - 一种电厂循环流化床锅炉操作优化方法 - Google Patents
一种电厂循环流化床锅炉操作优化方法 Download PDFInfo
- Publication number
- CN113504765A CN113504765A CN202110679824.6A CN202110679824A CN113504765A CN 113504765 A CN113504765 A CN 113504765A CN 202110679824 A CN202110679824 A CN 202110679824A CN 113504765 A CN113504765 A CN 113504765A
- Authority
- CN
- China
- Prior art keywords
- data
- boiler
- coal
- main steam
- power plant
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000003245 coal Substances 0.000 claims abstract description 72
- 238000005457 optimization Methods 0.000 claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 238000013528 artificial neural network Methods 0.000 claims abstract description 11
- 238000010219 correlation analysis Methods 0.000 claims abstract description 11
- 230000002159 abnormal effect Effects 0.000 claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 8
- 238000005065 mining Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 54
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 15
- 239000001301 oxygen Substances 0.000 claims description 15
- 229910052760 oxygen Inorganic materials 0.000 claims description 15
- 239000007789 gas Substances 0.000 claims description 12
- 230000000717 retained effect Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 230000007246 mechanism Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 238000010977 unit operation Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000009466 transformation Effects 0.000 abstract description 2
- 238000007621 cluster analysis Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000013213 extrapolation Methods 0.000 description 7
- 230000002776 aggregation Effects 0.000 description 6
- 238000004220 aggregation Methods 0.000 description 6
- 238000005265 energy consumption Methods 0.000 description 5
- 238000007418 data mining Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- RHZUVFJBSILHOK-UHFFFAOYSA-N anthracen-1-ylmethanolate Chemical compound C1=CC=C2C=C3C(C[O-])=CC=CC3=CC2=C1 RHZUVFJBSILHOK-UHFFFAOYSA-N 0.000 description 1
- 239000003830 anthracite Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002224 dissection Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 239000003925 fat Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32339—Object oriented modeling, design, analysis, implementation, simulation language
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Automation & Control Theory (AREA)
- Fluidized-Bed Combustion And Resonant Combustion (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
本发明公开了一种电厂循环流化床锅炉操作优化方法，包括：进行实时数据采集；通过剔除异常值、稳态分析及相关性分析进行数据预处理；使用神经网络算法建立锅炉数字模型；运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型；得到锅炉吨标煤产汽量最大的运行操作参数；经过关联规则计算后得到锅炉吨标煤产汽量最大的区间；运用多维聚类算法对区间进行聚类分析，得到聚类中心，建立优化操作模型；求出锅炉效率最高时所对应的各操作参数，并与当前运行参数对比，制定需要调整参数及范围；可为全厂锅炉的操作优化提供指导，产生明显的经济效益；无硬件改造，不影响电厂正常运行。
Description
技术领域
本发明属于智能工业领域，具体涉及一种电厂循环流化床锅炉操作优化方法。
背景技术
火电厂锅炉由于技术设计偏差、工程建设偏差和长期运行造成设备损坏、效率降低、以及运行管理的原因，导致在实际运行中偏离最佳运行点，从而导致效率偏低，能耗偏高。目前，尽管火电厂节能减排已经受到广泛的重视，但是在效果方面却始终不够显著。由于火力发电是一个基于朗肯循环的闭式系统，单一参数的改变会影响到很多参数，这种影响既有正向影响，又有负向影响。因此，对于综合权衡参数的大小对火电站效率的提高和能耗的降低具有十分重要的意义。能够在安全稳定的前提下，可靠的提供重要参数的优化值，为运行人员在实际运行中提供决策支持是当务之急和大势所趋。
随着信息、通讯、计算机、智能控制、变频技术的发展，为火力发电厂的高效、节约运作、科学管理，以及过程优化提供了前所未有的手段，进而促进火力发电厂的科学管理和自动化水平的提高。同时，计算机技术为火电厂保存大量的数据，它们能够反映电厂历史的运行情况。如果能够对这些数据进行分析处理，就能够针对电厂进行合理的“解剖”，从而能完全了解电厂实际性能和特性。并且能够根据历史优异的运行经验为现在运行操作提供参照。
当前能处理如此众多繁杂的数据，必须依赖大数据技术。大数据技术，是指从各种各样类型的数据中，快速获得有价值信息的能力。适用于大数据的技术，包括大规模并行处理(MPP)数据库，数据挖掘电网，分布式文件系统，分布式数据库，云计算平台，互联网，和可扩展的存储系统。有人把数据比喻为蕴藏能量的煤矿。煤炭按照性质有焦煤、无烟煤、肥煤、贫煤等分类，而露天煤矿、深山煤矿的挖掘成本又不一样。与此类似，大数据并不在“大”，而在于“有用”。价值含量、挖掘成本比数量更为重要。对于很多行业而言，如何利用这些大规模数据是赢得竞争的关键。对于火电厂而言，大数据技术通过分析和挖掘电厂运行数据，能够以能耗为指标，通过分析电厂实际运行相关参数与能耗的关系，分析并获得在特定工况下最优的操作建议，为运行人员提供决策支持。
本锅炉燃烧深度优化系统以DCS的数据以及日常其他数据为数据源，以大数据技术为核心，通过神经网络、关联规则等算法寻找出特定工况下最优运行参数，为电厂实际运行提供决策支持，从而有效的降低电厂能耗，为电厂带来可观的经济效益。
发明内容
为了实现上述目的，本发明涉及：一种电厂循环流化床锅炉操作优化方法，包括如下步骤：
步骤一.基于OPC服务进行实时数据采集；
步骤二.通过剔除异常值、稳态分析及相关性分析进行数据预处理；
步骤三.使用神经网络算法建立锅炉数字模型；
步骤四.运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型；得到锅炉吨标煤产汽量最大的运行操作参数；经过关联规则计算后得到锅炉吨标煤产汽量最大的区间；
步骤五.运用多维聚类算法对区间进行聚类分析，得到聚类中心，以此建立优化操作模型；
步骤六.求出锅炉效率最高时所对应的各操作参数，并与当前运行参数对比，制定需要调整参数及范围，给予运行人员操作指导。
进一步的，所述步骤一基于OPC服务进行实时数据采集的方法为：组建OPC数据采集网络以采集锅炉操作数据。建立专用数据库系统用于保存采集到的锅炉操作数据，为建立锅炉的精确模型做准备。
进一步的，所述步骤二中剔除异常值是指将数据中的NaN空数据和数值为0及违反常理的离群值删除。
进一步的，所述违反常理的判断标准为：开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中不可能为负值或者不可能为正值的数据；
开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中超出正常范围的数据；
停机时，开度、一次风量、二次风量、给水流量、氧量、真空、排烟温度、主蒸汽流量、循环水泵电流、磨煤机电流、负荷、给煤流量数据的非零显示；
空气温度高于100℃的数据。
进一步的，所述稳态分析的方法为：
以“主汽压力(MPa)”用滑动窗口法，对数据按照时间顺序进行稳态状态的判别，得到稳定无异常的运行数据；火电厂机组的运行状态经常发生变化，为了准确获得机组各变量间的关系，需要从采样数据中选择机组稳态运行数据；机组主汽压力是反映机组运行工况是否稳定的重要因素，因此在本方案中采用滑动窗口法对“主汽压力(MPa)”进行稳态区间检测；计算窗口内数据的波动情况，若波动情况较大，则认为窗口内数据处于非稳态，不保留数据；否则数据处于稳态，保留数据，通过对稳态数据进行关联分析，对负荷进行合理的切割，通过聚类分析，寻找出各类负荷下最优运行工况，从而为未来运行操作提供支持；
计算公式：
式中：δ表示稳态判定指标，t表示滑动窗口内采样时间；N表示滑动窗口的宽度，
表示从t到t+N-1之间实时的主蒸汽压力的均值；xt表示第t个工况的实时的主蒸汽压力值；λ为主汽压力量程范围的2％；
上式成立时候,则保留数据,不成立时，则不保留数据。
进一步的，所述步骤二中相关性分析的方法为：对给水流量、给水温度、总风量、一次风量、二次风量、氧量、引风机电流、送风机电流、一次风机电流、循环水进口温度、循环水泵电流、真空泵电流、磨煤机电流、磨煤机转速与煤耗进行相关性分析，将各变量与煤耗间的相关性系数进行排序，保留相关性强的前N个变量，其余变量全部删除。
进一步的，N的取值范围为14≥N≥11。
进一步的，所述步骤三中使用神经网络算法建立锅炉数字模型的方法为：采用BP神经网络算法，负荷、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压等17个变量，以预测进煤总热量与测试时间段内实际实时进煤总热量的方差和小于设定值；以预测数据与实际数据平均值差值小于设置值；以实际实时进煤总热量的差值最大值小于设定值。通过最优模型寻找法，寻找出拟合准确比率和模型总体偏移程度小的最优训练模型。
进一步的，所述锅炉数字模型包括两层：全厂锅炉数字模型和单锅炉数字模型；
1)全厂锅炉数字模型基于各给煤煤质和总负荷下整体锅炉效率数字模型，该模型最大作用是定量分析负荷和给煤煤质对全厂效率的影响，从而指导配煤和负荷分配；
2)单锅炉数字模型基于每台锅炉运行数据建立的反映单锅炉效率的数字模型，与单锅炉机理模型协调配合，最大程度再现锅炉运行状态，从而为优化模型做铺垫。
进一步的，所述步骤四中运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型的方法为；获取历史数据，完成数据的整理以及离散化后，根据Apriori法获取与机组高效率、低煤耗点相关的关联规则。将得到的关联规则进行整理形成关联规则库。进行优化时首先获取机组运行的实时数据，根据负荷及环境温度从关联规则库当中搜索符合负荷及环境温度的全部关联规则。然后在众多检索出的规则进行选取，规则选取时考虑机组的运行规程，工艺指标，同时对当前运行工况的综合性能与关联规则库中选取的工况综合性能进行评估，假若当前运行工况下的综合性能比关联规程库中的综合性能优异，则继续保持当前工况控制参数值，最终选取可控参数最优值对实际操作进行指导；
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
(1)本发明的电厂循环流化床锅炉操作优化方法，通过开发锅炉效率的机理模型和数学模型，建立针对全厂锅炉深度优化的优化模型并求解，可为全厂锅炉的操作优化提供指导，产生明显的经济效益。
(2)本发明的电厂循环流化床锅炉操作优化方法，无硬件改造，不影响电厂正常运行。
附图说明
图1为本发明较佳实施例的锅炉燃烧优化系统技术原理结构图；
图2为本发明较佳实施例的聚类分析流程图；
图3为本发明较佳实施例聚类分析对锅炉数据离散化处理示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
请参考图1-图2，一种电厂循环流化床锅炉操作优化方法，包括如下步骤：
步骤一.基于OPC服务进行实时数据采集；
步骤二.通过剔除异常值、稳态分析及相关性分析进行数据预处理；
步骤三.使用神经网络算法建立锅炉数字模型；
步骤四.运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型；得到锅炉吨标煤产汽量最大的运行操作参数；经过关联规则计算后得到锅炉吨标煤产汽量最大的区间；
步骤五.运用多维聚类算法对区间进行聚类分析，得到聚类中心，以此建立优化操作模型；
步骤六.求出锅炉效率最高时所对应的各操作参数，并与当前运行参数对比，制定需要调整参数及范围，给予运行人员操作指导。
所述步骤一基于OPC服务进行实时数据采集的方法为：组建OPC数据采集网络以采集锅炉操作数据。建立专用数据库系统用于保存采集到的锅炉操作数据，为建立锅炉的精确模型做准备。
所述步骤二中剔除异常值是指将数据中的NaN空数据和数值为0及违反常理的离群值删除。
所述违反常理的判断标准为：开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中不可能为负值或者不可能为正值的数据；
开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中超出正常范围的数据；
停机时，开度、一次风量、二次风量、给水流量、氧量、真空、排烟温度、主蒸汽流量、循环水泵电流、磨煤机电流、负荷、给煤流量数据的非零显示；
空气温度高于100℃的数据。
进一步的，所述稳态分析的方法为：
以“主汽压力(MPa)”用滑动窗口法，对数据按照时间顺序进行稳态状态的判别，得到稳定无异常的运行数据；火电厂机组的运行状态经常发生变化，为了准确获得机组各变量间的关系，需要从采样数据中选择机组稳态运行数据；机组主汽压力是反映机组运行工况是否稳定的重要因素，因此在本方案中采用滑动窗口法对“主汽压力(MPa)”进行稳态区间检测；计算窗口内数据的波动情况，若波动情况较大，则认为窗口内数据处于非稳态，不保留数据；否则数据处于稳态，保留数据，通过对稳态数据进行关联分析，对负荷进行合理的切割，通过聚类分析，寻找出各类负荷下最优运行工况，从而为未来运行操作提供支持；
计算公式：
式中：δ表示稳态判定指标，t表示滑动窗口内采样时间；N表示滑动窗口的宽度，
表示从t到t+N-1之间实时的主蒸汽压力的均值；xt表示第t个工况的实时的主蒸汽压力值；λ为主汽压力量程范围的2％；
上式成立时候,则保留数据,不成立时，则不保留数据。
所述步骤二中相关性分析的方法为：对给水流量、给水温度、总风量、一次风量、二次风量、氧量、引风机电流、送风机电流、一次风机电流、循环水进口温度、循环水泵电流、真空泵电流、磨煤机电流、磨煤机转速与煤耗进行相关性分析，将各变量与煤耗间的相关性系数进行排序，保留相关性强的前N个变量，其余变量全部删除。
N的取值范围为14≥N≥11。
所述步骤三中使用神经网络算法建立锅炉数字模型的方法为：采用BP神经网络算法，负荷、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压等17个变量，以预测进煤总热量与测试时间段内实际实时进煤总热量的方差和小于设定值；以预测数据与实际数据平均值差值小于设置值；以实际实时进煤总热量的差值最大值小于设定值。通过最优模型寻找法，寻找出拟合准确比率和模型总体偏移程度小的最优训练模型。
进一步的，所述步骤三中锅炉数字模型包括两层：全厂锅炉数字模型和单锅炉数字模型；
1)全厂锅炉数字模型基于各给煤煤质和总负荷下整体锅炉效率数字模型，该模型最大作用是定量分析负荷和给煤煤质对全厂效率的影响，从而指导配煤和负荷分配；
2)单锅炉数字模型基于每台锅炉运行数据建立的反映单锅炉效率的数字模型，与单锅炉机理模型协调配合，最大程度再现锅炉运行状态，从而为优化模型做铺垫。
所述步骤四中运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型的方法为；获取历史数据，完成数据的整理以及离散化后，根据Apriori法获取与机组高效率、低煤耗点相关的关联规则。将得到的关联规则进行整理形成关联规则库。进行优化时首先获取机组运行的实时数据，根据负荷及环境温度从关联规则库当中搜索符合负荷及环境温度的全部关联规则。然后在众多检索出的规则进行选取，规则选取时考虑机组的运行规程，工艺指标，同时对当前运行工况的综合性能与关联规则库中选取的工况综合性能进行评估，假若当前运行工况下的综合性能比关联规程库中的综合性能优异，则继续保持当前工况控制参数值，最终选取可控参数最优值对实际操作进行指导；
外推优化值
考虑机组在某些工况下可能从未达到最优运行，因此除从历史工况数据中选择的最优关联规则可控参数最优值外，可能存在更佳的可控参数值使机组运行效率更高、煤耗更低。
基于往坏的方向走将更坏，往好的方向走将更好的原理，通过以差类关联规则平均值为支点往好类关联规则优化点方向外推，控制步长控制外推范围，可以外推得到比历史数据中更好地可控参数最优值，如果外推可控参数最优值超过可操作范围，改变外推歩长，再次进行外推，最终得到不超过可操作范围的外推可控参数最优值。
以关联规则为支点往好关联规则优化点方向外推，得到比历史数据中更好地外推优化点，如果外推优化点操作变量值超过可操作范围，可以改变预测歩长，再次进行外推，得到操作变量值不超过可操作范围的外推优化点。
外推优化：
考虑机组在某些工况下可能从未达到最优运行，因此除从历史工况数据中选择的最优关联规则可控参数最优值外，可能存在更佳的可控参数值使机组运行效率更高、煤耗更低。
基于往坏的方向走将更坏，往好的方向走将更好的原理，通过以差类关联规则平均值为支点往好类关联规则优化点方向外推，控制步长控制外推范围，可以外推得到比历史数据中更好地可控参数最优值，如果外推可控参数最优值超过可操作范围，改变外推歩长，再次进行外推，最终得到不超过可操作范围的外推可控参数最优值。
以关联规则为支点往好关联规则优化点方向外推，得到比历史数据中更好地外推优化点，如果外推优化点操作变量值超过可操作范围，可以改变预测歩长，再次进行外推，得到操作变量值不超过可操作范围的外推优化点。
请参考图3，聚类分析是一种将研究对象分为相对同质的群组的统计分析技术。电厂中的参数大多是连续型数据，需要将连续型数据转换为离散型数据才能用关联规则进行挖掘。在本方案中聚类分析有三个作用：作用一：通过聚类分析查看各变量的聚集程度以及波动情况，然后根据各变量的聚集程度以及波动情况确定各变量类别数量。
作用二：通过聚类分析将电厂连续型数据划分成有意义的组，转换为离散型数据，用于关联规则数据挖掘。
作用三：通过聚类分析将电厂数据转换为离散型数据，提高了数据挖掘的查询和处理速度。数据离散化是指将连续的数据进行分段，使其变为一段段离散化的区间。分段的原则有基于等距离、等频率或优化的方法。
通过聚类分析，将发电煤耗分为三大类，其余变量分为两大类
首先通过聚类分析查看各变量的聚集程度以及波动情况，然后根据各变量的聚集程度以及波动情况确定各变量类别数量，多分配变量波动情况大且聚集程度比较分散的变量类别数量，适当减少变量波动情况小且聚集程度集中的变量类别数量。
每一类的标准为各变量的数据到此聚类中心的距离最小，将数据对象划分到距离最近的类中，以“距离”作为分类的依据，距离越小，数据越“相似”。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种电厂循环流化床锅炉操作优化方法，其特征在于，包括如下步骤：
步骤一.基于OPC服务进行实时数据采集；
步骤二.通过剔除异常值、稳态分析及相关性分析进行数据预处理；
步骤三.使用神经网络算法建立锅炉数字模型；
步骤四.运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型；得到锅炉吨标煤产汽量最大的运行操作参数；经过关联规则计算后得到锅炉吨标煤产汽量最大的区间；
步骤五.运用多维聚类算法对区间进行聚类分析，得到聚类中心，以此建立优化操作模型；
步骤六.求出锅炉效率最高时所对应的各操作参数，并与当前运行参数对比，制定需要调整参数及范围，给予运行人员操作指导。
2.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤一基于OPC服务进行实时数据采集的方法为：组建OPC数据采集网络以采集锅炉操作数据；建立专用数据库系统用于保存采集到的锅炉操作数据，为建立锅炉的精确模型做准备。
3.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤二中剔除异常值是指将数据中的NaN空数据和数值为0及违反常理的离群值删除。
4.根据权利要求3所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述违反常理的判断标准为：开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中不可能为负值或者不可能为正值的数据；
开度、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压、负荷、给煤流量中超出正常范围的数据；
停机时，开度、一次风量、二次风量、给水流量、氧量、真空、排烟温度、主蒸汽流量、循环水泵电流、磨煤机电流、负荷、给煤流量数据的非零显示；
空气温度高于100℃的数据。
5.根据权利要求3所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述稳态分析的方法为：
以“主汽压力(MPa)”用滑动窗口法，对数据按照时间顺序进行稳态状态的判别，得到稳定无异常的运行数据；火电厂机组的运行状态经常发生变化，为了准确获得机组各变量间的关系，需要从采样数据中选择机组稳态运行数据；机组主汽压力是反映机组运行工况是否稳定的重要因素，因此采用滑动窗口法对“主汽压力(MPa)”进行稳态区间检测；计算窗口内数据的波动情况，若波动情况较大，则认为窗口内数据处于非稳态，不保留数据；否则数据处于稳态，保留数据，通过对稳态数据进行关联分析，对负荷进行切割，通过聚类分析，寻找出各类负荷下最优运行工况，从而为未来运行操作提供支持；
计算公式：
式中：δ表示稳态判定指标，t表示滑动窗口内采样时间；N表示滑动窗口的宽度，
表示从t到t+N-1之间实时的主蒸汽压力的均值；xt表示第t个工况的实时的主蒸汽压力值；λ为主汽压力量程范围的2％；
上式成立时候,则保留数据,不成立时，则不保留数据。
6.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤二中相关性分析的方法为：对给水流量、给水温度、总风量、一次风量、二次风量、氧量、引风机电流、送风机电流、一次风机电流、循环水进口温度、循环水泵电流、真空泵电流、磨煤机电流、磨煤机转速与煤耗进行相关性分析，将各变量与煤耗间的相关性系数进行排序，保留相关性强的前N个变量，其余变量全部删除。
7.根据权利要求6所述的电厂循环流化床锅炉操作优化方法，其特征在于，N的取值范围为14≥N≥11。
8.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤三中使用神经网络算法建立锅炉数字模型的方法为：采用BP神经网络算法，负荷、一次风量、二次风量、给水流量、给水温度、空预器入口温度、氧量、真空、排烟温度、主蒸汽压力、主蒸汽温度、主蒸汽流量、循环水泵电流、循环水入口温度、磨煤机电流、炉膛负压17个变量，以预测进煤总热量与测试时间段内实际实时进煤总热量的方差和小于设定值；以预测数据与实际数据平均值差值小于设置值；以实际实时进煤总热量的差值最大值小于设定值；通过最优模型寻找法，寻找出拟合准确比率和模型总体偏移程度小的最优训练模型。
9.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤三中锅炉数字模型包括两层：全厂锅炉数字模型和单锅炉数字模型；
其中，全厂锅炉数字模型基于各给煤煤质和总负荷下整体锅炉效率数字模型，该模型最大作用是定量分析负荷和给煤煤质对全厂效率的影响，从而指导配煤和负荷分配；
单锅炉数字模型基于每台锅炉运行数据建立的反映单锅炉效率的数字模型，与单锅炉机理模型协调配合，最大程度再现锅炉运行状态，从而为优化模型做铺垫。
10.根据权利要求1所述的电厂循环流化床锅炉操作优化方法，其特征在于，所述步骤四中运用关联规则对锅炉运行数据进行挖掘，建立操作优化模型的方法为；获取历史数据，完成数据的整理以及离散化后，根据Apriori法获取与机组高效率、低煤耗点相关的关联规则；将得到的关联规则进行整理形成关联规则库；进行优化时首先获取机组运行的实时数据，根据负荷及环境温度从关联规则库当中搜索符合负荷及环境温度的全部关联规则；然后在众多检索出的规则进行选取，规则选取时考虑机组的运行规程，工艺指标，同时对当前运行工况的综合性能与关联规则库中选取的工况综合性能进行评估，假若当前运行工况下的综合性能比关联规程库中的综合性能优异，则继续保持当前工况控制参数值，最终选取可控参数最优值对实际操作进行指导。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110679824.6A CN113504765A (zh) | 2021-06-18 | 2021-06-18 | 一种电厂循环流化床锅炉操作优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110679824.6A CN113504765A (zh) | 2021-06-18 | 2021-06-18 | 一种电厂循环流化床锅炉操作优化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113504765A true CN113504765A (zh) | 2021-10-15 |
Family
ID=78010522
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110679824.6A Pending CN113504765A (zh) | 2021-06-18 | 2021-06-18 | 一种电厂循环流化床锅炉操作优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113504765A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113790766A (zh) * | 2021-11-15 | 2021-12-14 | 江苏未来智慧信息科技有限公司 | 电厂智慧侦测系统及方法 |
| CN114413249A (zh) * | 2022-03-29 | 2022-04-29 | 朗坤智慧科技股份有限公司 | 一种电站锅炉效率对标优化的数据分析方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106094744A (zh) * | 2016-06-04 | 2016-11-09 | 上海大学 | 基于关联规则挖掘的火电厂主运行参数目标值的确定方法 |
| CN111091251A (zh) * | 2020-03-24 | 2020-05-01 | 汉谷云智(武汉)科技有限公司 | 一种基于大数据技术的锅炉操作优化方法及系统 |
| CN111158239A (zh) * | 2020-01-09 | 2020-05-15 | 上海电力大学 | 基于关联规则算法与神经网络的制粉系统性能优化方法 |
| CN111401652A (zh) * | 2020-03-24 | 2020-07-10 | 汉谷云智(武汉)科技有限公司 | 一种基于co在线检测的锅炉优化方法及系统 |
-
2021
- 2021-06-18 CN CN202110679824.6A patent/CN113504765A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106094744A (zh) * | 2016-06-04 | 2016-11-09 | 上海大学 | 基于关联规则挖掘的火电厂主运行参数目标值的确定方法 |
| CN111158239A (zh) * | 2020-01-09 | 2020-05-15 | 上海电力大学 | 基于关联规则算法与神经网络的制粉系统性能优化方法 |
| CN111091251A (zh) * | 2020-03-24 | 2020-05-01 | 汉谷云智(武汉)科技有限公司 | 一种基于大数据技术的锅炉操作优化方法及系统 |
| CN111401652A (zh) * | 2020-03-24 | 2020-07-10 | 汉谷云智(武汉)科技有限公司 | 一种基于co在线检测的锅炉优化方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 王晓峰: "基于数据挖掘的电站锅炉运行优化研究", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113790766A (zh) * | 2021-11-15 | 2021-12-14 | 江苏未来智慧信息科技有限公司 | 电厂智慧侦测系统及方法 |
| CN114413249A (zh) * | 2022-03-29 | 2022-04-29 | 朗坤智慧科技股份有限公司 | 一种电站锅炉效率对标优化的数据分析方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108954375B (zh) | 锅炉节煤控制方法 | |
| CN104534507B (zh) | 一种锅炉燃烧优化控制方法 | |
| CN112066355B (zh) | 基于数据驱动的余热锅炉阀门自适应调节方法 | |
| CN113504765A (zh) | 一种电厂循环流化床锅炉操作优化方法 | |
| CN101893877A (zh) | 一种发电厂基于能耗分析的优化运行方法及其系统 | |
| CN112861379B (zh) | 基于稀疏大数据挖掘的火电机组汽轮机优化方法及系统 | |
| CN101187804A (zh) | 基于数据挖掘的火电机组运行优化规则提取方法 | |
| CN109872012A (zh) | 基于工况划分的火电厂运行多目标优化的确定方法 | |
| CN111539845B (zh) | 一种基于用电模式隶属评分的企业环保管控响应研判方法 | |
| CN104807039A (zh) | 一种锅炉燃烧优化的变量降维建模方法 | |
| CN112801388B (zh) | 一种基于非线性时间序列算法的电力负荷预测方法及系统 | |
| CN112149714A (zh) | 基于数据挖掘确定煤电机组能效特征指标基准值的方法 | |
| CN111210170A (zh) | 基于90%用电分布特征指标的环保管控监测及评价方法 | |
| CN111199304B (zh) | 基于数据驱动融合策略的多目标燃烧优化方法 | |
| CN111401652A (zh) | 一种基于co在线检测的锅炉优化方法及系统 | |
| Zhang et al. | Research on a TOPSIS energy efficiency evaluation system for crude oil gathering and transportation systems based on a GA-BP neural network | |
| CN112488380B (zh) | 基于相似度动态模型的机组稳态工况匹配方法及系统 | |
| CN114970766A (zh) | 一种基于线性拟合的电站机组运行参数基准值获取方法 | |
| CN114139619A (zh) | 基于改进K-means算法的锅炉燃烧优化控制方法及装置 | |
| CN114066212A (zh) | 一种基于历史工况的机组升降负荷工况寻优方法 | |
| CN113222216A (zh) | 冷热电负荷预测方法、装置及系统 | |
| CN111091251A (zh) | 一种基于大数据技术的锅炉操作优化方法及系统 | |
| Zheng et al. | Combustion process modeling based on deep sparse least squares support vector regression | |
| CN111651938B (zh) | 基于热力计算与大数据的变煤质机组出力预测方法 | |
| Zhang et al. | Markov hot blast stove prediction based on hybrid intelligent algorithm |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211015 |
|
| RJ01 | Rejection of invention patent application after publication |