考虑天气影响的火电机组空冷系统性能预测方法
Performance Prediction Method for Air Cooling System of Thermal Power Unit Considering Weather Effect
收稿日期: 2023-09-06 修回日期: 2023-10-26
基金资助: |
|
Received: 2023-09-06 Revised: 2023-10-26
作者简介 About authors





目的 直接空冷机组是一部分缺水地区常用的火力发电形式,由于其采用空气作为冷却介质,因此运行受到诸多限制。为解决直接空冷机组受环境影响大和耗煤量高的问题,对空冷岛换热性能进行预测研究。 方法 基于河北省某超临界2×600 MW机组的历史运行数据,利用MATLAB软件计算历史空冷岛性能,将历史数据作为训练集和测试集,通过长短期记忆(long short-term memory,LSTM)神经网络机器学习算法,实现对未来一段时间内的空冷岛性能预测。在不改变模型参数的条件下,通过去除各项特征的方式确定特征重要性排名,基于此确定最佳特征选择策略,进一步优化模型。考虑到空冷岛性能受天气影响大,为提升模型在特殊天气下的预测能力,将原数据集结合历史天气预报数据,编写考虑天气因素的预测程序,以预测空冷岛未来一段时间性能,并对预测结果进行可视化分析。 结果 所采用的预测模型预测准确度明显高于传统自回归移动平均模型(autoregressive integrated moving average model,ARIMA),对未来1 h以内的直接空冷机组换热性能预测拟合优度均在0.90以上。 结论 模型所采用的数据特征及算法可以为直接空冷机组的稳定运行提供数据支撑,为智慧电厂的建设提供技术基础。
关键词:
Objectives Direct air-cooled unit is a common equipment of thermal power generation in some water-deficient areas. The operation is subject to many restrictions because it uses air as its cooling medium. Heat transfer performance of air-cooled island was studied to solve these problems that direct air-cooled units are greatly affected by the environment and have high coal consumption. Methods Based on history-data of a supercritical 2×600 MW unit in Hebei Province, the performance of its air-cooled island was calculated with MATLAB software, this study considered the acquired data as the training set and the test set,which were used to predict future performance in virtue of long short-term memory (LSTM) neural network machine learning algorithm. Under the condition that the model parameters were not changed, the feature importance ranking was determined by removing all features, based on which the best feature selection strategy was determined to further optimize the model. Considering the great impact from the weather, a prediction procedure, taking into account weather factors, was written to improve the accuracy of predicting air-cooled island performance, by combining the original data set with historical weather data. Accordingly prediction results were subjected to visualization and analyzation. Results The prediction accuracy of the adopted prediction model is significantly higher than that of the traditional autoregressive integrated moving average model (ARIMA), and the goodness of fit of the direct air-cooled unit heat transfer performance prediction within the next hour is above 0.90. Conclusions The data characteristics and algorithms used in the model can provide data support for the stable operation of the direct air-cooled unit and provide a technical basis for the construction of intelligent power plants.
Keywords:
本文引用格式
董建宁, 安吉振, 陈衡, 潘佩媛, 徐钢, 王修彦.
DONG Jianning, AN Jizhen, CHEN Heng, PAN Peiyuan, XU Gang, WANG Xiuyan.
0 引言
随着中国经济飞速发展,用电需求随之攀升,高功率火力发电机组在中国华北、东北及西北等地区广泛投入使用。而以上地区普遍“富煤贫水”。直接空冷机组凭借其优越的节水能力以及占地面积小等优势,成为上述地区火力发电的常规冷却手段[1]。直接空冷机组采用空气作为汽轮机排汽的冷却介质,在空冷岛上对高温蒸汽进行冷却。相较于水冷,其性能受环境影响更大。较差的运行环境会直接影响空冷岛性能,存在着夏季温度高、背压高、负荷受限以及冬季低温冷凝管束冻结等问题[2]。通过气象站天气预报信息能实现未来天气的提前预知,如果能在考虑未来天气影响的情况下对空冷岛性能进行预测,则可以据此有针对性地调整运行策略来减小极端天气对机组运行的影响。
数据是机器学习算法的“燃料”,用于训练模型并使其能够做出预测和决策。机器学习通过分析数据中的模式和关系来学习,从而构建出能够处理新数据的模型。而实际运行的火电机组会产生大量的运行数据,将其与气象数据相结合,二者作为典型的时序数据,可通过LSTM算法进行分析,实现对空冷岛实时性能的预测。同时考虑到直接空冷机组受环境影响大,结合历史天气预报数据扩充数据集,本文着重分析天气因素对空冷岛性能的影响,可以提升直接空冷机组在不同气象条件下的应对能力,为机组运行提供数据依据。
1 长短期记忆神经网络
图1
图1所示结构公式表示如下:
式中:σ为激活函数,将输出量控制在[0,1];xt为输入量;ht-1为上一神经元的输出量,在这一层作为输入量输入; Wf、 Wi、 bf、bi、bc和 b0为参数矩阵,隐藏层通过这些参数来对数据进行运算和舍弃;
其中,
2 数据处理和模型建立
2.1 数据来源
研究数据来源于河北省某超临界2×600 MW直接空冷机组,包括大气压力、环境温度、风速、风向、排汽压力、排汽温度、主汽压力、主汽温度、风机转速、主汽流量、给水流量、凝汽器真空、凝结水流量、凝结水温度14项运行数据,其空冷岛由56台风机单元组成,按照8排7列布置,总换热面积为13 739 m2。利用机组的监控信息系统提取2018年6月30日至2018年9月26日的一部分运行数据,时间步长选定为5 min,并结合某超临界2×600 MW直接空冷机组所在城市气象站对应时间的历史天气预报数据,包括温度、风速、风向、气压在内的共计25 633条数据。
2.2 数据预处理
实测数据一般存在噪声、缺失值和离群值等异常情况[14]。在建立预测模型前,需要对实测数据进行预处理。首先填补缺失值,对于运行数据以及气象数据缺失值,由于其为连续时序数据,可以利用缺失位置的上下平均值来进行填补。
与此同时,对数据范围进行筛查,去除负荷小于250 MW的数据值,对数据进行统计后发现离群值处在8月25日至8月31日的停机阶段,为保证时序数据的时间步长尽可能均匀,选择剔除这一时间段的数据,以保证对数据的计算和分析是准确、有意义的,最终数据共23 743条。
2.3 空冷岛性能计算
利用处理后的运行实时数据,估算空冷岛实时换热性能。在MATLAB软件上使用相关热力学计算函数,利用排汽压力与排汽干度计算汽轮机末级焓值;利用凝结水温度以及凝汽器真空计算凝结水焓值,再进行关于传热过程的热力学计算:
式中:fn为凝结水流量,kg/s;fs为排汽流量,kg/s;hs为汽轮机末级焓值,kJ;hn为凝结水焓值,kJ;Qn为蒸汽凝气放热量,kJ。
根据风机相似定律[15],计算变工况下的空气流量:
式中:n为实际风机转速,Hz;n0为额定风机转速,Hz;
根据能量守恒定律,计算空气出口温度:
式中:ta1为环境温度,即进口空气温度,℃;Cp为空气的定压比热容。
利用空气和水蒸气进出口温度计算对数平均换热温差:
式中:T1为排汽温度,℃;T2为凝结水温度,℃。
对表征空冷岛换热性能的传热系数K进行计算:
式中:A0为空冷岛总换热面积,m2。
图2
2.4 数据集构造
选定预测数据为未来60 min内的空冷岛实时性能,学习时长为4 h,利用python计算机语言对数据进行维度处理和时序对应。为保证数据的有效性,需要对数据进行
式中:Xnorm为归一化后的数据;X、Xmean、Xstd分别为原数据值、平均值和标准差。
2.5 模型建立
图3
表1 初次建模预测值与实际值拟合优度
Tab. 1
时间/ min | 5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 | 60 |
---|---|---|---|---|---|---|---|---|---|---|---|---|
R2 | 0.95 | 0.95 | 0.94 | 0.93 | 0.92 | 0.91 | 0.89 | 0.87 | 0.85 | 0.83 | 0.81 | 0.78 |
2.6 特征重要性分析
对于LSTM来说,在Tensorflow平台无法直接查看特征重要性[20],而灰色关联度分析法对于LSTM模型并不准确,其计算方法难以代表长序列历史数据,这就需要定义特征重要性[21-22]。对于机器学习模型来说,不同特征会不同程度地影响模型表现,因此,可以通过人为地在某项特征中添加噪声或去除某项特征,在不改变原模型参数的情况下重新预测,并用同一种误差评价指标来对预测效果进行评价。如果误差波动幅度大,说明此项特征较为重要;如果误差波动幅度小,说明此项特征相对不重要。将归一化处理后的测试集某项特征数值设为平均值,以消除该特征的影响,重新进行预测,误差评价指标设为R2。对每一项特征分别进行上述处理,多次预测后得到R2变化见表2。
表2 去除各项特征后拟合优度对比
Tab. 2
参数 | 时间/min | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 | 60 | |
原测试集 | 0.95 | 0.95 | 0.94 | 0.93 | 0.92 | 0.91 | 0.89 | 0.87 | 0.85 | 0.83 | 0.81 | 0.78 |
去大气压 | 0.97 | 0.97 | 0.96 | 0.96 | 0.95 | 0.93 | 0.92 | 0.90 | 0.89 | 0.87 | 0.85 | 0.82 |
去温度 | 0.94 | 0.93 | 0.91 | 0.90 | 0.89 | 0.88 | 0.87 | 0.85 | 0.85 | 0.83 | 0.82 | 0.80 |
去风速 | 0.94 | 0.94 | 0.92 | 0.92 | 0.91 | 0.89 | 0.87 | 0.85 | 0.83 | 0.81 | 0.78 | 0.76 |
去风向 | 0.96 | 0.95 | 0.94 | 0.94 | 0.92 | 0.91 | 0.89 | 0.87 | 0.86 | 0.83 | 0.81 | 0.78 |
去排汽压力 | 0.90 | 0.90 | 0.87 | 0.86 | 0.85 | 0.82 | 0.78 | 0.76 | 0.73 | 0.70 | 0.66 | 0.63 |
去排汽温度 | 0.94 | 0.93 | 0.91 | 0.90 | 0.88 | 0.86 | 0.83 | 0.81 | 0.77 | 0.75 | 0.72 | 0.68 |
去主汽压力 | 0.94 | 0.94 | 0.93 | 0.93 | 0.92 | 0.91 | 0.89 | 0.88 | 0.86 | 0.84 | 0.82 | 0.80 |
去主汽温度 | 0.96 | 0.95 | 0.94 | 0.94 | 0.92 | 0.91 | 0.89 | 0.88 | 0.86 | 0.83 | 0.81 | 0.79 |
去风机转速 | 0.93 | 0.93 | 0.92 | 0.92 | 0.91 | 0.89 | 0.87 | 0.86 | 0.85 | 0.83 | 0.81 | 0.79 |
去主汽流量 | 0.96 | 0.96 | 0.95 | 0.94 | 0.93 | 0.92 | 0.91 | 0.89 | 0.88 | 0.85 | 0.83 | 0.81 |
去给水流量 | 0.96 | 0.96 | 0.94 | 0.94 | 0.93 | 0.91 | 0.89 | 0.88 | 0.86 | 0.84 | 0.82 | 0.79 |
去凝器真空 | 0.87 | 0.88 | 0.87 | 0.87 | 0.86 | 0.85 | 0.85 | 0.83 | 0.83 | 0.80 | 0.78 | 0.77 |
去凝结水流量 | 0.94 | 0.94 | 0.93 | 0.92 | 0.91 | 0.90 | 0.88 | 0.87 | 0.85 | 0.83 | 0.81 | 0.79 |
去凝结水温度 | 0.94 | 0.93 | 0.92 | 0.92 | 0.91 | 0.90 | 0.88 | 0.86 | 0.85 | 0.83 | 0.80 | 0.78 |
计算每次训练R2的变化幅度,根据其大小绘制特征重要性排名如图4所示。其中,大气压、主汽流量、给水流量、主汽温度、主汽压力、风向等特征会降低预测准确度,在新预测模型将进行去除。由于与空冷岛换热性能相关的环境特征为风速和环境温度,因此可以结合相关的天气预报数据来对模型准确度进行提升。
图4
3 模型优化与对比分析
3.1 结合天气预报数据优化模型
重新构造数据集,根据特征重要性排名结果,在不改变其他参数的情况下对原无效特征进行去除,重新构造模型,新模型拟合优度见表3。
表3 去除无效特征后预测值R2
Tab. 3
时间/min | 5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 | 60 |
---|---|---|---|---|---|---|---|---|---|---|---|---|
R2 | 0.98 | 0.98 | 0.98 | 0.97 | 0.97 | 0.96 | 0.95 | 0.94 | 0.93 | 0.91 | 0.90 | 0.89 |
相较原模型,新模型准确度明显提升,模型训练时间缩短。根据特征重要性排名,有针对地结合历史天气预报数据中温度和风速2项重要特征,重新构造数据。将未来4 h气象数据上移,与实时运行数据对应,使天气预报数据相较于运行数据具有超前性,模型在实际应用中能结合未来天气预报对空冷岛性能进行预测。重新创建模型并训练,带入测试集得到预测值,计算拟合优度见表4。
表4 结合天气预报数据后预测值拟合优度
Tab. 4
时间/min | 5 | 10 | 15 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 | 60 |
---|---|---|---|---|---|---|---|---|---|---|---|---|
R2 | 0.99 | 0.98 | 0.98 | 0.98 | 0.98 | 0.97 | 0.96 | 0.95 | 0.95 | 0.93 | 0.92 | 0.91 |
图5
图5
空冷岛预测性能与实际性能对比
Fig. 5
Comparison between predicted performance and actual performance of air-cooled island
3.2 LSTM模型与ARIMA模型对比
为体现LSTM在时序数据预测处理方面的优越性,利用相同的训练数据创建整合自回归移动平均模型(autoregressive integrated moving average model,ARIMA)来进行预测效果对比。ARIMA是一种数据驱动的时间序列预测模型[23]。取2018年9月20日时间序列和空冷岛换热系数作为训练数据,进行ARIMA建模。
首先对数据进行平稳性校验,采用单位根校验法(augmented dickey-fuller,ADF),先假设时序数据非稳定。统计学根据显著性检验方法所得到的假定值p-value,一般p-value<0.05为显著,p-value<0.01为非常显著。计算所选取的时序数据假定值得到p-value=0.390 121>0.05,T检验的检验统计量T-statistic=-1.780 628>临界值critical value(1%)=-3.455 754。根据统计学指标,无法拒绝原假设,数据非稳定,无法直接绘制自相关(auto correlation function,ACF)图和偏自相关(partial auto correlation function,PACF)图。由于空冷岛总体换热系数大于0,因此可以对数据进行对数变换和一阶差分,以提升数据稳定性。处理后p-value=0.011 384<0.05,T-statistic=-3.487 738<critical value(1%)=-3.455 853,经过处理的数据满足稳定条件,绘制ACF图和PACF图[24]如图6所示。
图6
自相关图有拖尾现象,偏自相关图为1阶截尾,因此选定自回归项数p=1,滑动平均项数q=0,差分阶数d=1。创建ARIMA模型,输入复原后的空冷岛性能数据,对空冷岛性能进行预测,计算预测误差指标函数,与同一时间的LSTM预测结果进行可视化对比,如图7所示。
图7
图7
LSTM预测效果与ARIMA预测效果对比
Fig. 7
Comparison of prediction effect between LSTM and ARIMA
由图7可知,尽管ARIMA模型所需数据量相较LSTM模型少,但使用LSTM的拟合优度明显好于ARIMA,LSTM模型误差评价指标低于ARIMA误差评价指标。因此,LSTM算法在空冷岛性能预测方面相较传统ARIMA模型具有优越性。
4 结论
基于河北省某电厂历史实时运行数据对空冷岛换热性能进行计算,基于LSTM长短期记忆神经网络建立空冷岛性能预测模型,利用历史数据对模型进行训练,使模型具有一定的预测未来空冷岛性能的能力。通过对特征重要性进行分析,得到增删各项特征的最优策略,进一步优化模型性能。具体结论如下:
1)通过对特征重要性计算,在数值层面上证明空冷岛换热性能除了受一部分运行数据的影响外,同时还受环境温度、风速等因素影响。
2)去除与空冷岛性能不相关特征能明显提升模型训练时间,同时有效提升模型预测准确度。有针对地结合天气预报温度和风速数据,进一步提升了模型准确度,模型对未来1 h空冷岛性能预测值与实际值的拟合优度均大于0.9,拟合效果很好。
3)通过对比ARIMA模型与LSTM模型预测结果,发现LSTM在空冷岛实时性能预测方面具有优越性,预测效果明显优于ARIMA模型。
4)所建空冷岛换热性能预测模型能够实现对空冷岛低性能情况进行预警,进而指导机组运行中风机风量、背压等参数的调节,以提升机组的热效率,提高经济性。
5)所建立的模型着重考虑了未来天气对于空冷性能的影响,结合未来气象数据,使直接空冷机组能够针对不同气象的条件做出运行调整,同时提高运行安全性,为“智慧电厂”的建设和普及提供了技术基础。
参考文献
发电厂直接空冷技术研究现状及展望
[J].
Research status and prospect of direct air cooling technology in power plants
[J].
基于空冷岛性能在线监测的风机优化调节
[J].
Optimization and regulation of fans based on online monitoring of air-cooled island
[J].
基于大数据与人工智能的燃煤电站直接空冷系统状态评估和运行优化
[D].
State evaluation and operation optimization of direct air-cooling system in coal-fired power plant based on big data and artificial intelligence
[D].
基于模型预测控制的源荷储低碳经济调度
[J].
Low-carbon economic dispatching of source,load and storage based on model predictive control
[J].
面向新型电力系统智能化提效的多源异构数据融合技术研究
[J].
Research on multi-source heterogeneous data integration technology for smart efficiency improvement in novel power system
[J].
基于Attention-LSTM与多模型集成的短期负荷预测方法
[J].
Short-term load forecasting method based on Attention-LSTM and multi-model integration
[J].
基于模型预测的VSC-MTDC协调控制策略
[J].
Coordinate control strategy of VSC-MTDC based on model predictive control
[J].
短期电力负荷预测模型的建立与应用
[J].
Establishment and application of short-term power load forecasting
[J].
基于深度学习LSTM网络的短期电力负荷预测方法
[J].
Short-term electrical load forecasting based on deep learning LSTM networks
[J].
基于长短期记忆神经网络的火电厂NO x 排放预测模型
[J].
Prediction model for NOx emissions from thermal power plants based on long-short-term memory neural network
[J].
Long short-term memory
[J].
Speech recognition with deep recurrent neural networks
[C]//
一种改进组合神经网络的超短期风速预测方法研究
[J].
Research on an ultra-short-term wind speed prediction method based on improved combined neural networks
[J].
基于优化ELM算法的电站凝汽器故障诊断
[D].
Fault diagnosis for power plant condenser based on optimized ELM algorithm
[D].
图像降噪的自适应高斯平滑滤波器
[J].
Adaptive Gaussian smoothing filter for image denoising
[J].
CaO carbonation kinetics determined using micro-fluidized bed thermogravimetric analysis
[J].
谷歌TensorFlow机器学习框架及应用
[J].
Application of Google TensorFlow machine learning framework
[J].
卷积神经网络研究综述
[J].
Review of convolutional neural network
[J].
联合循环发电站电力输出预测
[J].
Electrical power output prediction of combined cycle power station
[J].
卷积神经网络特征重要性分析及增强特征选择模型
[J].
Convolution neural network feature importance analysis and feature selection enhanced model
[J].
碳中和背景下多通道特征组合超短期风电功率预测
[J].
A multi-channel feature combination model for Ultra-short-term wind power prediction under carbon neutral background
[J].
Time series analysis forecasting and control-rev.ed
.[J].
/
〈 |
|
〉 |
