CN113111588A - 一种燃气轮机nox排放浓度预测方法及装置 - Google Patents

一种燃气轮机nox排放浓度预测方法及装置 Download PDF

Info

Publication number
CN113111588A
CN113111588A CN202110442134.9A CN202110442134A CN113111588A CN 113111588 A CN113111588 A CN 113111588A CN 202110442134 A CN202110442134 A CN 202110442134A CN 113111588 A CN113111588 A CN 113111588A
Authority
CN
China
Prior art keywords
feature
input
data
sample
data set
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
CN202110442134.9A
Other languages
English (en)
Other versions
CN113111588B (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.)
Datang Boiler Pressure Vessel Examination Center Co Ltd
East China Electric Power Test Institute of China Datang Corp Science and Technology Research Institute Co Ltd
Original Assignee
Datang Boiler Pressure Vessel Examination Center Co Ltd
East China Electric Power Test Institute of China Datang Corp Science and Technology Research Institute 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 Datang Boiler Pressure Vessel Examination Center Co Ltd, East China Electric Power Test Institute of China Datang Corp Science and Technology Research Institute Co Ltd filed Critical Datang Boiler Pressure Vessel Examination Center Co Ltd
Priority to CN202110442134.9A priority Critical patent/CN113111588B/zh
Publication of CN113111588A publication Critical patent/CN113111588A/zh
Application granted granted Critical
Publication of CN113111588B publication Critical patent/CN113111588B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种燃气轮机NOX排放浓度预测方法及装置,所述方法包括:建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;对向量集合进行预处理,得到预处理后的数据集D*;分别采用偏互信息PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure DDA0003035413660000011
利用非线性多特征选择组合预测模型进行NOX排放浓度预测;本发明的优点在于:预测精度较高。

Description

一种燃气轮机NOX排放浓度预测方法及装置
技术领域
本发明涉及NOX浓度预测领域,更具体涉及一种燃气轮机NOX排放浓度预测方法及装置。
背景技术
煤炭发电以低廉的经济优势占据着我国能源结构的主要部分,煤炭发电过程中会释放大量NOX污染物,随着严格的环境法规的推进,加大了煤电环保工作的难度。与煤炭发电相比,燃气蒸汽联合循环机组以天然气为介质燃料,发电产物清洁,NOX排放浓度仅为煤炭电厂的1/10。
燃机NOX排放浓度能够反映燃机的燃烧健康水平,是建立燃烧调整模型的重要变量。异常的排放特性影响机组燃烧效率,燃烧脉动不稳定,触发负荷超驰机制。准确的预测NOx可以对异常工况进行预警,消除安全隐患。
构造准确预测模型的核心是针对研究问题找到适合的建模方法。在已有研究中,建模方法以传统物理分析方法、统计学方法和数据驱动的先进智能方法三种。其中,曹军等发表文献《F级燃气-蒸汽联合循环机组高精度全范围仿真***的开发》,利用APROS软件对燃机热力过程进行质量、动量分析,遵循能量守恒法则对F级燃气蒸汽机组全范围过程进行实时动态仿真。然而,该机理方法在建模精度方面有提升空间。厂级数据监控***(Supervisory Information System,SIS)作为智慧电厂控制的子***实现了机组过程控制中大量参数数据状态的监测记录与实时访问,为先进智能建模控制策略提供了保障。其中,李景轩等发表《燃气轮机机理-数据混合建模方法研究》,设计了一种以智能算法作为机理模型误差补偿的混合模型控制器方法,并对不同组合方式的设计进行了基于分布式控制*** (Distributed Control System,DCS)数据的验证实验,将预测精度提高。云世豪等发表《基于RS-RBF的燃气轮机控制***传感器故障诊断研究》,将传感器故障征兆属性通过改进的等频离散方法处理构建粗糙集,并进行下一步RBF网络的建立减少传感器故障的误判率。已有研究在电站参数建模方面取得了成功,然而以上研究均为浅层机器学习方法,不能够捕捉隐藏在数据底层的深层次有用信息,导致预测精度不够高。
发明内容
本发明所要解决的技术问题在于现有技术NOX浓度预测方法预测精度不够高的问题。
本发明通过以下技术手段实现解决上述技术问题的:一种燃气轮机 NOX排放浓度预测方法,所述方法包括:
步骤一:TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;
步骤二:对向量集合进行预处理,得到预处理后的数据集D*
步骤三:分别采用偏互信息PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure BDA0003035413640000021
步骤四:基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000031
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。
本发明分别通过PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择并排序得到了新的样本集
Figure BDA0003035413640000032
基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000033
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测,预测过程不是单一算法预测结果,而且通过DNN算法对多个预测结果拟合得到较为准确的预测结果,预测精度较高。
进一步地,所述步骤一包括:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure BDA0003035413640000034
N 为特征个数,
Figure BDA0003035413640000035
特征中含有m个数值,记为m个样本;输出特征 Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure BDA0003035413640000036
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure BDA0003035413640000037
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure BDA0003035413640000038
Figure BDA0003035413640000039
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
进一步地,所述步骤二包括:采用DBSCAN异常点检测方法与Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
更进一步地,所述步骤三包括:
步骤3-1:数据集D*作为初始数据集,通过公式
Figure BDA0003035413640000041
计算初始数据集中特征
Figure BDA0003035413640000042
Figure BDA0003035413640000043
的互信息
Figure BDA0003035413640000044
并选出最大互信息时对应的特征移入初始集合S中;其中, f(·)为基于m组样本的估计密度函数;
步骤3-2:通过公式
Figure BDA0003035413640000045
计算数据集D*中剩余特征
Figure BDA0003035413640000046
中的每一个特征剔除S信息影响后输入残差V,通过公式
Figure BDA0003035413640000047
计算
Figure BDA0003035413640000048
中剔除S信息影响后的输出残差U;其中,E(·) 表示条件期望;
步骤3-3:计算I(V,U),并找出互信息最大时的变量
Figure BDA0003035413640000049
步骤3-4:将
Figure BDA00030354136400000410
放入集合S中更新集合S得到更新后的集合S’,将更新后的集合S’作为初始集合S,将变量集
Figure BDA00030354136400000411
作为初始数据集,返回执行步骤3-2-步骤3-4,直到变量集C为空为止;
步骤3-5:根据所有变量在计算最大互信息时所移入集合S中的变量顺序,得到所有输入变量重要信息排序
Figure BDA00030354136400000412
更进一步地,所述步骤三还包括:
步骤3-6:通过公式
Figure BDA00030354136400000413
计算集合D*中某一个特征的第 i个样本的Gini系数,其中,
Figure BDA00030354136400000414
为集合D*中某一个特征的第i个样本划分为左集合后的方差,
Figure BDA00030354136400000415
为集合D*中某一个特征的第i个样本划分为右集合后的方差;
步骤3-7:循环步骤3-6计算集合D*中所有特征,得到每个特征的Gini 系数以及其对应的切分点;
步骤3-8:选择Gini系数最小的特征以及对应的切分点为最优特征和最优切分点,根据最小的特征以及对应的切分点将该属性下的样本划分为二部分子集D1与子集D2;
步骤3-9:将子集D1作为集合D*返回执行步骤3-6至步骤3-8,对子集D1继续划分,将子集D2作为集合D*返回执行步骤3-6至步骤3-8,对子集D2继续划分,直到每个特征不再有除该特征本身以外的子集,所有特征划分完成,生成CART树结构;
步骤3-10:通过公式
Figure BDA0003035413640000051
计算特征重要性信息;其中,N是特征总数,m_t是当前节点特征数,Gini表示当前特征的基尼系数,m_t_R为当前节点***右孩子数,GiniR对应当前节点***的右孩子基尼系数,m_t_L为当前节点***左孩子数,GiniL对应当前节点***的左孩子基尼系数;
步骤3-11:将获得的所有特征重要性信息从大到小排序,排序结果对应的特征的顺序作为输入变量重要信息排序
Figure BDA0003035413640000052
更进一步地,所述步骤三还包括:
步骤3-12:调用matlab中LASSO特征选择变量工具包中的LASSO(X,Y) 函数,计算得到D*中的输入特征
Figure BDA0003035413640000053
对应的系数矩阵H;
步骤3-13:对系数矩阵H从左到右分析,记录每一行最先出现0元素时对应的列,将所有的列中最先出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000061
的最后一位,将所有的列中最后出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000062
的第一位,按照此排序规则对D*中的所有的输入特征进行排序得到最终的输入变量重要信息排序
Figure BDA0003035413640000063
步骤3-14:将输入变量重要信息排序
Figure BDA0003035413640000064
输入变量重要信息排序
Figure BDA0003035413640000065
输入变量重要信息排序
Figure BDA0003035413640000066
中前10个特征分别与输出特征构成新的样本集
Figure BDA0003035413640000067
更进一步地,所述步骤四包括:
步骤4-1:将新的样本集
Figure BDA0003035413640000068
分别输入深度置信网络DBN进行建模预测,分别得到NOx浓度预测结果值y1、y2、y3;
步骤4-2:将DNN网络的输入神经元个数设置为3,输出神经元个数设置为1,构建DNN模型;
步骤4-3:将NOx浓度预测结果值y1、y2、y3作为DNN模型的建模输入变量,实际测量值y作为DNN模型的建模输出变量,对NOx浓度预测结果进行误差修正,得到最终预测结果。
本发明还提供一种燃气轮机NOX排放浓度预测装置,所述装置包括:
数据集获取模块,用于TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;
预处理模块,用于对向量集合进行预处理,得到预处理后的数据集D*
特征选择模块,用于分别采用偏互信息PMI、决策树CART、套索回归 LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure BDA0003035413640000071
预测模块,用于基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000072
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。
进一步地,所述数据集获取模块还用于:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure BDA0003035413640000073
N 为特征个数,
Figure BDA0003035413640000074
特征中含有m个数值,记为m个样本;输出特征 Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure BDA0003035413640000075
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure BDA0003035413640000076
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure BDA0003035413640000077
Figure BDA0003035413640000078
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
进一步地,所述预处理模块还用于:采用DBSCAN异常点检测方法与 Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
更进一步地,所述特征选择模块还用于:
步骤3-1:数据集D*作为初始数据集,通过公式
Figure BDA0003035413640000079
计算初始数据集中特征
Figure BDA00030354136400000710
Figure BDA00030354136400000711
的互信息
Figure BDA0003035413640000081
并选出最大互信息时对应的特征移入初始集合S中;其中, f(·)为基于m组样本的估计密度函数;
步骤3-2:通过公式
Figure BDA0003035413640000082
计算数据集D*中剩余特征
Figure BDA0003035413640000083
中的每一个特征剔除S信息影响后输入残差V,通过公式
Figure BDA0003035413640000084
计算
Figure BDA0003035413640000085
中剔除S信息影响后的输出残差U;其中,E(·) 表示条件期望;
步骤3-3:计算I(V,U),并找出互信息最大时的变量
Figure BDA0003035413640000086
步骤3-4:将
Figure BDA0003035413640000087
放入集合S中更新集合S得到更新后的集合S’,将更新后的集合S’作为初始集合S,将变量集
Figure BDA0003035413640000088
作为初始数据集,返回执行步骤3-2-步骤3-4,直到变量集C为空为止;
步骤3-5:根据所有变量在计算最大互信息时所移入集合S中的变量顺序,得到所有输入变量重要信息排序
Figure BDA0003035413640000089
更进一步地,所述特征选择模块还用于:
步骤3-6:通过公式
Figure BDA00030354136400000810
计算集合D*中某一个特征的第 i个样本的Gini系数,其中,
Figure BDA00030354136400000811
为集合D*中某一个特征的第i个样本划分为左集合后的方差,
Figure BDA00030354136400000812
为集合D*中某一个特征的第i个样本划分为右集合后的方差;
步骤3-7:循环步骤3-6计算集合D*中所有特征,得到每个特征的Gini 系数以及其对应的切分点;
步骤3-8:选择Gini系数最小的特征以及对应的切分点为最优特征和最优切分点,根据最小的特征以及对应的切分点将该属性下的样本划分为二部分子集D1与子集D2;
步骤3-9:将子集D1作为集合D*返回执行步骤3-6至步骤3-8,对子集D1继续划分,将子集D2作为集合D*返回执行步骤3-6至步骤3-8,对子集D2继续划分,直到每个特征不再有除该特征本身以外的子集,所有特征划分完成,生成CART树结构;
步骤3-10:通过公式
Figure BDA0003035413640000091
计算特征重要性信息;其中,N是特征总数,m_t是当前节点特征数,Gini表示当前特征的基尼系数,m_t_R为当前节点***右孩子数,GiniR对应当前节点***的右孩子基尼系数,m_t_L为当前节点***左孩子数,GiniL对应当前节点***的左孩子基尼系数;
步骤3-11:将获得的所有特征重要性信息从大到小排序,排序结果对应的特征的顺序作为输入变量重要信息排序
Figure BDA0003035413640000092
更进一步地,所述特征选择模块还用于:
步骤3-12:调用matlab中LASSO特征选择变量工具包中的LASSO(X,Y) 函数,计算得到D*中的输入特征
Figure BDA0003035413640000093
对应的系数矩阵H;
步骤3-13:对系数矩阵H从左到右分析,记录每一行最先出现0元素时对应的列,将所有的列中最先出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000094
的最后一位,将所有的列中最后出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000095
的第一位,按照此排序规则对D*中的所有的输入特征进行排序得到最终的输入变量重要信息排序
Figure BDA0003035413640000096
步骤3-14:将输入变量重要信息排序
Figure BDA0003035413640000097
输入变量重要信息排序
Figure BDA0003035413640000098
输入变量重要信息排序
Figure BDA0003035413640000099
中前10个特征分别与输出特征构成新的样本集
Figure BDA0003035413640000101
更进一步地,所述预测模块还用于:
步骤4-1:将新的样本集
Figure BDA0003035413640000102
分别输入深度置信网络DBN进行建模预测,分别得到NOx浓度预测结果值y1、y2、y3;
步骤4-2:将DNN网络的输入神经元个数设置为3,输出神经元个数设置为1,构建DNN模型;
步骤4-3:将NOx浓度预测结果值y1、y2、y3作为DNN模型的建模输入变量,实际测量值y作为DNN模型的建模输出变量,对NOx浓度预测结果进行误差修正,得到最终预测结果。
本发明的优点在于:本发明分别通过PMI、决策树CART、套索回归 LASSO对预处理后的数据集D*进行特征选择并排序得到了新的样本集
Figure BDA0003035413640000103
基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000104
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测,预测过程不是单一算法预测结果,而且通过DNN算法对多个预测结果拟合得到较为准确的预测结果,预测精度较高。。
附图说明
图1为本发明实施例所提供的一种燃气轮机NOX排放浓度预测方法的流程图;
图2为本发明实施例所提供的一种燃气轮机NOX排放浓度预测方法的预测结果与PMI、CART、LASSO三种特征选择算法进行比较的示意图;
图3为采用相同LASSO特征选择算法不同建模模型的NOx排放预测结果的相对误差箱型图;
图4为本发明实施例所提供的一种燃气轮机NOX排放浓度预测方法与 3种特征选择算法经过非线性组合前后的NOx排放预测值拟合实际测量值的散点图对比,其中(a)、(b)、(c)、(d)分别表示的是PMI、CART、LASSO、本发明的预测方法的NOx排放预测值拟合实际测量值的散点图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1所示,一种燃气轮机NOX排放浓度预测方法,所述方法包括:
步骤S1:TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;具体过程为:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure BDA0003035413640000111
N 为特征个数,
Figure BDA0003035413640000112
特征中含有m个数值,记为m个样本;输出特征 Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure BDA0003035413640000113
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure BDA0003035413640000121
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure BDA0003035413640000122
Figure BDA0003035413640000123
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
步骤S2:对向量集合进行预处理,得到预处理后的数据集D*;具体过程为:采用DBSCAN异常点检测方法与Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
步骤S3:分别采用偏互信息PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure BDA0003035413640000124
具体过程为:
步骤3-1:数据集D*作为初始数据集,通过公式
Figure BDA0003035413640000125
计算初始数据集中特征
Figure BDA0003035413640000126
Figure BDA0003035413640000127
的互信息
Figure BDA0003035413640000128
并选出最大互信息时对应的特征移入初始集合S中;其中, f(·)为基于m组样本的估计密度函数;
步骤3-2:通过公式
Figure BDA0003035413640000129
计算数据集D*中剩余特征
Figure BDA00030354136400001210
中的每一个特征剔除S信息影响后输入残差V,通过公式
Figure BDA00030354136400001211
计算
Figure BDA00030354136400001212
中剔除S信息影响后的输出残差U;其中,E(·) 表示条件期望;
步骤3-3:计算I(V,U),并找出互信息最大时的变量
Figure BDA0003035413640000131
步骤3-4:通过公式
Figure BDA0003035413640000132
计算AIC值,其中,u为已选变量计算的U回归残差;p为已选变量个数;AIC值用于表示
Figure BDA0003035413640000133
对应 U的剩余信息,将
Figure BDA0003035413640000134
放入集合S中更新集合S得到更新后的集合S’,将更新后的集合S’作为初始集合S,将变量集
Figure BDA0003035413640000135
作为初始数据集,返回执行步骤3-2-步骤3-4,直到变量集C为空为止,得到所有变量对应的AIC 值,变量对应的AIC值能够反映该变量对模型复杂度的影响情况,对于模型每加入一个变量这个模型复杂度越高了那么就可以少选几个变量特征,相反,如果对于模型每加入一个变量这个模型复杂度越低了那么就可以多选几个变量特征,本实施例中选10个特征;
步骤3-5:根据所有变量在计算最大互信息时所移入集合S中的变量顺序,得到所有输入变量重要信息排序
Figure BDA0003035413640000136
步骤3-6:通过公式
Figure BDA0003035413640000137
计算集合D*中某一个特征的第i个样本的Gini系数,其中,
Figure BDA0003035413640000138
为集合D*中某一个特征的第i个样本划分为左集合后的方差,
Figure BDA0003035413640000139
为集合D*中某一个特征的第i个样本划分为右集合后的方差;
步骤3-7:循环步骤3-6计算集合D*中所有特征,得到每个特征的Gini 系数以及其对应的切分点;
步骤3-8:选择Gini系数最小的特征以及对应的切分点为最优特征和最优切分点,根据最小的特征以及对应的切分点将该属性下的样本划分为二部分子集D1与子集D2;
步骤3-9:将子集D1作为集合D*返回执行步骤3-6至步骤3-8,对子集D1继续划分,将子集D2作为集合D*返回执行步骤3-6至步骤3-8,对子集D2继续划分,直到每个特征不再有除该特征本身以外的子集,所有特征划分完成,生成CART树结构;
步骤3-10:通过公式
Figure BDA0003035413640000141
计算特征重要性信息;其中,N是特征总数,m_t是当前节点特征数,Gini表示当前特征的基尼系数,m_t_R为当前节点***右孩子数,GiniR对应当前节点***的右孩子基尼系数,m_t_L为当前节点***左孩子数,GiniL对应当前节点***的左孩子基尼系数;
步骤3-11:将获得的所有特征重要性信息从大到小排序,排序结果对应的特征的顺序作为输入变量重要信息排序
Figure BDA0003035413640000142
步骤3-12:调用matlab中LASSO特征选择变量工具包中的LASSO(X,Y) 函数,计算得到D*中的输入特征
Figure BDA0003035413640000143
对应的系数矩阵H;
步骤3-13:对系数矩阵H从左到右分析,记录每一行最先出现0元素时对应的列,将所有的列中最先出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000144
的最后一位,将所有的列中最后出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000145
的第一位,按照此排序规则对D*中的所有的输入特征进行排序得到最终的输入变量重要信息排序
Figure BDA0003035413640000146
步骤3-14:将输入变量重要信息排序
Figure BDA0003035413640000147
输入变量重要信息排序
Figure BDA0003035413640000148
输入变量重要信息排序
Figure BDA0003035413640000149
中前10个特征分别与输出特征构成新的样本集
Figure BDA00030354136400001410
步骤S4:基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA00030354136400001411
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。具体过程为:
步骤4-1:将新的样本集
Figure BDA0003035413640000151
分别输入深度置信网络DBN进行建模预测,分别得到NOx浓度预测结果值y1、y2、y3;
步骤4-2:将DNN网络的输入神经元个数设置为3,输出神经元个数设置为1,构建DNN模型;
步骤4-3:将NOx浓度预测结果值y1、y2、y3作为DNN模型的建模输入变量,实际测量值y作为DNN模型的建模输出变量,对NOx浓度预测结果进行误差修正,得到最终预测结果。
以下分析本发明的实验结果:
本研究数据来自某热电电厂1#美国GE公司生产的,型号为PG9371FB 的燃机TCS数据采集***,采集了与NOx排放浓度相关的参数特征30个,采样频率为1min,共252组样本。
为了验证MFSA算法(本发明提供的NOx排放浓度预测方法)性能,与PMI、CART、LASSO三种特征选择算法进行比较,为不失一般性,均采样DBN网络作为建模模型。图2为各方法在验证集下的预测值折线图,从图中可以看出三种特征选择算法与经过非线性组合预测的特性选择算法均对实际NOx排放有很高的预测能力,但本发明提供的组合方法更能真实反应排放量的变化趋势。
根据评价指标对各方法进行进一步的比较分析,如表1所示。从表1可以看出MFSA算法的三种指标值均比其他算法有所提高,采用MFSA算法比次优的PMI算法的MAPE、RMSE、MAE、R2精度分别提高了10%、4.76%、11.76%、0.1%,说明MFSA算法达到了提高特征选择精度的要求,算法适用。
表1不同算法评价指标结果
Figure BDA0003035413640000161
图3所示为采用相同LASSO特征选择算法不同建模模型的NOx排放预测结果的相对误差箱型图,其中对比建模模型分别为BP算法与SVR算法。从线型图可以明显看出,深度网络DBN的预测能力比其他传统浅层机器学习算法BP与SVR更好,预测值相对误差上下四分位数更接近0附近。
图4中(a)、(b)、(c)、(d)分别表示的是3种特征选择算法以及本发明的 MFSA算法经过非线性组合前后的NOx排放预测值拟合实际测量值的散点图,其中黑色星号代表拟合分布,黑色实对角线为理想分布线,R2代表测量值拟合预测值的程度,值越大,模型预测精度越高。从图4可以看出图(d) 的R2值最大并且预测值均匀的集中在理想曲线附近,实验结果说明所提算法能提高模型的预测精度。
通过以上技术方案,本发明分别通过PMI、决策树CART、套索回归 LASSO对预处理后的数据集D*进行特征选择并排序得到了新的样本集
Figure BDA0003035413640000162
基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000171
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测,预测过程不是单一算法预测结果,而且通过DNN算法对多个预测结果拟合得到较为准确的预测结果,预测精度较高。
实施例2
与本发明实施例1相对应的,本发明实施例2还提供一种燃气轮机NOX排放浓度预测装置,所述装置包括:
数据集获取模块,用于TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;
预处理模块,用于对向量集合进行预处理,得到预处理后的数据集D*
特征选择模块,用于分别采用偏互信息PMI、决策树CART、套索回归 LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure BDA0003035413640000172
预测模块,用于基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure BDA0003035413640000173
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。
具体的,所述数据集获取模块还用于:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure BDA0003035413640000174
N 为特征个数,
Figure BDA0003035413640000181
特征中含有m个数值,记为m个样本;输出特征 Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure BDA0003035413640000182
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure BDA0003035413640000183
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure BDA0003035413640000184
Figure BDA0003035413640000185
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
具体的,所述预处理模块还用于:采用DBSCAN异常点检测方法与 Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
更具体的,所述特征选择模块还用于:
步骤3-1:数据集D*作为初始数据集,通过公式
Figure BDA0003035413640000186
计算初始数据集中特征
Figure BDA0003035413640000187
Figure BDA0003035413640000188
的互信息
Figure BDA0003035413640000189
并选出最大互信息时对应的特征移入初始集合S中;其中, f(·)为基于m组样本的估计密度函数;
步骤3-2:通过公式
Figure BDA00030354136400001810
计算数据集D*中剩余特征
Figure BDA00030354136400001811
中的每一个特征剔除S信息影响后输入残差V,通过公式
Figure BDA00030354136400001812
计算
Figure BDA00030354136400001813
中剔除S信息影响后的输出残差U;其中,E(·) 表示条件期望;
步骤3-3:计算I(V,U),并找出互信息最大时的变量
Figure BDA00030354136400001814
步骤3-4:将
Figure BDA00030354136400001815
放入集合S中更新集合S得到更新后的集合S’,将更新后的集合S’作为初始集合S,将变量集
Figure BDA0003035413640000191
作为初始数据集,返回执行步骤3-2-步骤3-4,直到变量集C为空为止;
步骤3-5:根据所有变量在计算最大互信息时所移入集合S中的变量顺序,得到所有输入变量重要信息排序
Figure BDA0003035413640000192
更具体的,所述特征选择模块还用于:
步骤3-6:通过公式
Figure BDA0003035413640000193
计算集合D*中某一个特征的第 i个样本的Gini系数,其中,
Figure BDA0003035413640000194
为集合D*中某一个特征的第i个样本划分为左集合后的方差,
Figure BDA0003035413640000195
为集合D*中某一个特征的第i个样本划分为右集合后的方差;
步骤3-7:循环步骤3-6计算集合D*中所有特征,得到每个特征的Gini 系数以及其对应的切分点;
步骤3-8:选择Gini系数最小的特征以及对应的切分点为最优特征和最优切分点,根据最小的特征以及对应的切分点将该属性下的样本划分为二部分子集D1与子集D2;
步骤3-9:将子集D1作为集合D*返回执行步骤3-6至步骤3-8,对子集D1继续划分,将子集D2作为集合D*返回执行步骤3-6至步骤3-8,对子集D2继续划分,直到每个特征不再有除该特征本身以外的子集,所有特征划分完成,生成CART树结构;
步骤3-10:通过公式
Figure BDA0003035413640000196
计算特征重要性信息;其中,N是特征总数,m_t是当前节点特征数,Gini表示当前特征的基尼系数,m_t_R为当前节点***右孩子数,GiniR对应当前节点***的右孩子基尼系数,m_t_L为当前节点***左孩子数,GiniL对应当前节点***的左孩子基尼系数;
步骤3-11:将获得的所有特征重要性信息从大到小排序,排序结果对应的特征的顺序作为输入变量重要信息排序
Figure BDA0003035413640000201
更具体的,所述特征选择模块还用于:
步骤3-12:调用matlab中LASSO特征选择变量工具包中的LASSO(X,Y) 函数,计算得到D*中的输入特征
Figure BDA0003035413640000202
对应的系数矩阵H;
步骤3-13:对系数矩阵H从左到右分析,记录每一行最先出现0元素时对应的列,将所有的列中最先出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000203
的最后一位,将所有的列中最后出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure BDA0003035413640000204
的第一位,按照此排序规则对D*中的所有的输入特征进行排序得到最终的输入变量重要信息排序
Figure BDA0003035413640000205
步骤3-14:将输入变量重要信息排序
Figure BDA0003035413640000206
输入变量重要信息排序
Figure BDA0003035413640000207
输入变量重要信息排序
Figure BDA0003035413640000208
中前10个特征分别与输出特征构成新的样本集
Figure BDA0003035413640000209
更具体的,所述预测模块还用于:
步骤4-1:将新的样本集
Figure BDA00030354136400002010
分别输入深度置信网络DBN进行建模预测,分别得到NOx浓度预测结果值y1、y2、y3;
步骤4-2:将DNN网络的输入神经元个数设置为3,输出神经元个数设置为1,构建DNN模型;
步骤4-3:将NOx浓度预测结果值y1、y2、y3作为DNN模型的建模输入变量,实际测量值y作为DNN模型的建模输出变量,对NOx浓度预测结果进行误差修正,得到最终预测结果。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种燃气轮机NOX排放浓度预测方法,其特征在于,所述方法包括:
步骤一:TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;
步骤二:对向量集合进行预处理,得到预处理后的数据集D*
步骤三:分别采用偏互信息PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure FDA0003035413630000011
步骤四:基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure FDA0003035413630000012
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。
2.根据权利要求1所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤一包括:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure FDA0003035413630000013
N为特征个数,
Figure FDA0003035413630000014
特征中含有m个数值,记为m个样本;输出特征Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure FDA0003035413630000015
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure FDA0003035413630000016
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure FDA0003035413630000017
Figure FDA0003035413630000018
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
3.根据权利要求1所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤二包括:采用DBSCAN异常点检测方法与Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
4.根据权利要求3所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤三包括:
步骤3-1:数据集D*作为初始数据集,通过公式
Figure FDA0003035413630000021
计算初始数据集中特征
Figure FDA0003035413630000022
Figure FDA0003035413630000023
的互信息
Figure FDA0003035413630000024
并选出最大互信息时对应的特征移入初始集合S中;其中,f(·)为基于m组样本的估计密度函数;
步骤3-2:通过公式
Figure FDA0003035413630000025
计算数据集D*中剩余特征
Figure FDA0003035413630000026
中的每一个特征剔除S信息影响后输入残差V,通过公式
Figure FDA0003035413630000027
计算
Figure FDA0003035413630000028
中剔除S信息影响后的输出残差U;其中,E(·)表示条件期望;
步骤3-3:计算I(V,U),并找出互信息最大时的变量
Figure FDA0003035413630000029
步骤3-4:将
Figure FDA00030354136300000210
放入集合S中更新集合S得到更新后的集合S’,将更新后的集合S’作为初始集合S,将变量集
Figure FDA00030354136300000211
作为初始数据集,返回执行步骤3-2-步骤3-4,直到变量集C为空为止;
步骤3-5:根据所有变量在计算最大互信息时所移入集合S中的变量顺序,得到所有输入变量重要信息排序
Figure FDA00030354136300000212
5.根据权利要求4所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤三还包括:
步骤3-6:通过公式
Figure FDA0003035413630000031
计算集合D*中某一个特征的第i个样本的Gini系数,其中,
Figure FDA0003035413630000032
为集合D*中某一个特征的第i个样本划分为左集合后的方差,
Figure FDA0003035413630000033
为集合D*中某一个特征的第i个样本划分为右集合后的方差;
步骤3-7:循环步骤3-6计算集合D*中所有特征,得到每个特征的Gini系数以及其对应的切分点;
步骤3-8:选择Gini系数最小的特征以及对应的切分点为最优特征和最优切分点,根据最小的特征以及对应的切分点将该属性下的样本划分为二部分子集D1与子集D2;
步骤3-9:将子集D1作为集合D*返回执行步骤3-6至步骤3-8,对子集D1继续划分,将子集D2作为集合D*返回执行步骤3-6至步骤3-8,对子集D2继续划分,直到每个特征不再有除该特征本身以外的子集,所有特征划分完成,生成CART树结构;
步骤3-10:通过公式
Figure FDA0003035413630000034
计算特征重要性信息;其中,N是特征总数,m_t是当前节点特征数,Gini表示当前特征的基尼系数,m_t_R为当前节点***右孩子数,GiniR对应当前节点***的右孩子基尼系数,m_t_L为当前节点***左孩子数,GiniL对应当前节点***的左孩子基尼系数;
步骤3-11:将获得的所有特征重要性信息从大到小排序,排序结果对应的特征的顺序作为输入变量重要信息排序
Figure FDA00030354136300000410
6.根据权利要求5所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤三还包括:
步骤3-12:调用matlab中LASSO特征选择变量工具包中的LASSO(X,Y)函数,计算得到D*中的输入特征
Figure FDA0003035413630000041
对应的系数矩阵H;
步骤3-13:对系数矩阵H从左到右分析,记录每一行最先出现0元素时对应的列,将所有的列中最先出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure FDA0003035413630000042
的最后一位,将所有的列中最后出现0元素的列对应的输入特征放在输入变量重要信息排序
Figure FDA0003035413630000043
的第一位,按照此排序规则对D*中的所有的输入特征进行排序得到最终的输入变量重要信息排序
Figure FDA0003035413630000044
步骤3-14:将输入变量重要信息排序
Figure FDA0003035413630000045
输入变量重要信息排序
Figure FDA0003035413630000046
输入变量重要信息排序
Figure FDA0003035413630000047
中前10个特征分别与输出特征构成新的样本集
Figure FDA0003035413630000048
7.根据权利要求6所述的一种燃气轮机NOX排放浓度预测方法,其特征在于,所述步骤四包括:
步骤4-1:将新的样本集
Figure FDA0003035413630000049
分别输入深度置信网络DBN进行建模预测,分别得到NOx浓度预测结果值y1、y2、y3;
步骤4-2:将DNN网络的输入神经元个数设置为3,输出神经元个数设置为1,构建DNN模型;
步骤4-3:将NOx浓度预测结果值y1、y2、y3作为DNN模型的建模输入变量,实际测量值y作为DNN模型的建模输出变量,对NOx浓度预测结果进行误差修正,得到最终预测结果。
8.一种燃气轮机NOX排放浓度预测装置,其特征在于,所述装置包括:
数据集获取模块,用于TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA并将样本数据集DATA转换为向量形式,得到向量集合;
预处理模块,用于对向量集合进行预处理,得到预处理后的数据集D*
特征选择模块,用于分别采用偏互信息PMI、决策树CART、套索回归LASSO对预处理后的数据集D*进行特征选择,并根据变量排序原则得到各算法进行特征选择后的重要信息排序,选择每个重要信息排序的前10个特征分别构成新的样本集
Figure FDA0003035413630000051
预测模块,用于基于深度置信网络DBN对PMI、CART、LASSO算法下构成的样本集
Figure FDA0003035413630000052
分别进行建模预测,并采用DNN算法构建非线性多特征选择组合预测模型,利用非线性多特征选择组合预测模型进行NOX排放浓度预测。
9.根据权利要求8所述的一种燃气轮机NOX排放浓度预测装置,其特征在于,所述数据集获取模块还用于:
从TCS数据采集***中获取与NOX排放浓度相关的燃烧状态参数与控制参数原始数据,建立样本数据集DATA=X∪Y,其中,
Figure FDA0003035413630000053
N为特征个数,
Figure FDA0003035413630000054
特征中含有m个数值,记为m个样本;输出特征Y={y1,y2,…,ym},输出特征Y中的每一个数值对应
Figure FDA0003035413630000055
中每个样本的NOX排放浓度值,构造输入输出矩阵
Figure FDA0003035413630000056
将输入输出矩阵D改写为由N+1个列向量构成的向量集合
Figure FDA0003035413630000057
Figure FDA0003035413630000058
表示向量集合中第一个元素,其实际是输入输出矩阵D中第一列向量,对应于样本数据集的第一个特征。
10.根据权利要求8所述的一种燃气轮机NOX排放浓度预测装置,其特征在于,所述预处理模块还用于:采用DBSCAN异常点检测方法与Savitzky-Golay平滑滤波相结合的方法对向量集合进行数据预处理,并将处理后数据映射到[0,1]区间,得到预处理后的数据集D*
CN202110442134.9A 2021-04-23 2021-04-23 一种燃气轮机nox排放浓度预测方法及装置 Active CN113111588B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110442134.9A CN113111588B (zh) 2021-04-23 2021-04-23 一种燃气轮机nox排放浓度预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110442134.9A CN113111588B (zh) 2021-04-23 2021-04-23 一种燃气轮机nox排放浓度预测方法及装置

Publications (2)

Publication Number Publication Date
CN113111588A true CN113111588A (zh) 2021-07-13
CN113111588B CN113111588B (zh) 2023-01-24

Family

ID=76719852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110442134.9A Active CN113111588B (zh) 2021-04-23 2021-04-23 一种燃气轮机nox排放浓度预测方法及装置

Country Status (1)

Country Link
CN (1) CN113111588B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444295A (zh) * 2022-01-21 2022-05-06 南京天洑软件有限公司 一种基于组合时序预测的NOx废气浓度预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108226046A (zh) * 2017-12-26 2018-06-29 中国科学技术大学 一种机动车尾气遥测结果修正方法
CN108647483A (zh) * 2018-04-26 2018-10-12 国家电投集团河南电力有限公司技术信息中心 一种基于模糊树建模方法的scr入口nox浓度的软测量方法
CN110379463A (zh) * 2019-06-05 2019-10-25 山东大学 基于机器学习的海洋藻类成因分析及浓度预测方法和***
WO2019209753A1 (en) * 2018-04-22 2019-10-31 Viome, Inc. Systems and methods for inferring scores for health metrics
CN110532674A (zh) * 2019-08-27 2019-12-03 东北电力大学 一种燃煤电站锅炉炉膛温度测量方法
CN110675920A (zh) * 2019-10-22 2020-01-10 华北电力大学 一种基于MI-LSTM的锅炉NOx预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108226046A (zh) * 2017-12-26 2018-06-29 中国科学技术大学 一种机动车尾气遥测结果修正方法
WO2019209753A1 (en) * 2018-04-22 2019-10-31 Viome, Inc. Systems and methods for inferring scores for health metrics
CN108647483A (zh) * 2018-04-26 2018-10-12 国家电投集团河南电力有限公司技术信息中心 一种基于模糊树建模方法的scr入口nox浓度的软测量方法
CN110379463A (zh) * 2019-06-05 2019-10-25 山东大学 基于机器学习的海洋藻类成因分析及浓度预测方法和***
CN110532674A (zh) * 2019-08-27 2019-12-03 东北电力大学 一种燃煤电站锅炉炉膛温度测量方法
CN110675920A (zh) * 2019-10-22 2020-01-10 华北电力大学 一种基于MI-LSTM的锅炉NOx预测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ZHIHONG ZHANG等: "High-order covariate interacted Lasso for feature selection", 《PATTERN RECOGNITION LETTERS》 *
彭岩等: "PM2.5浓度预测与影响因素分析", 《北京邮电大学学报》 *
曹薇: "《区域技术创新影响因素、网络特征及空间效应研究》", 30 September 2018 *
杨国田等: "基于变量选择的深度置信神经网络锅炉烟气NOX排放预测", 《热力发电》 *
胡忠义: "《基于计算智能技术的电力负荷预测理论及应用》", 31 January 2019 *
黄孝平: "《当代机器深度学习方法与应用研究》", 30 November 2017 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114444295A (zh) * 2022-01-21 2022-05-06 南京天洑软件有限公司 一种基于组合时序预测的NOx废气浓度预测方法
CN114444295B (zh) * 2022-01-21 2023-04-25 南京天洑软件有限公司 一种基于组合时序预测的NOx废气浓度预测方法

Also Published As

Publication number Publication date
CN113111588B (zh) 2023-01-24

Similar Documents

Publication Publication Date Title
CN110377984B (zh) 一种工业设备剩余有效寿命预测方法、***及电子设备
CN109491816A (zh) 基于知识的故障诊断方法
CN110119854A (zh) 基于代价敏感lstm循环神经网络的稳压器水位预测方法
Lin et al. Machine learning templates for QCD factorization in the search for physics beyond the standard model
CN111144609A (zh) 一种锅炉废气排放预测模型建立方法、预测方法及装置
CN111639783A (zh) 一种基于lstm神经网络的线损预测方法及***
CN109145948A (zh) 一种基于集成学习的注塑机堵料异常检测方法
CN111461463A (zh) 一种基于tcn-bp的短期负荷预测方法、***及设备
CN110987436B (zh) 基于激励机制的轴承故障诊断方法
CN112434848A (zh) 基于深度信念网络的非线性加权组合风电功率预测方法
CN114492191A (zh) 基于dbn-svr的热力站设备剩余寿命评估方法
CN112861436A (zh) 一种发动机排放实时预测方法
CN114169434A (zh) 一种负荷预测方法
CN113449919A (zh) 一种基于特征和趋势感知的用电量预测方法及***
CN116542701A (zh) 一种基于cnn-lstm组合模型的碳价预测方法及***
CN115345297A (zh) 一种基于生成对抗网络的台区样本生成方法及***
CN117556369B (zh) 一种动态生成的残差图卷积神经网络的窃电检测方法及***
CN113111588B (zh) 一种燃气轮机nox排放浓度预测方法及装置
CN110837857A (zh) 工业用电负荷预测方法、***及其存储介质
CN108537581B (zh) 基于gmdh选择性组合的能源消费量时间序列预测方法及装置
CN112801388B (zh) 一种基于非线性时间序列算法的电力负荷预测方法及***
CN114548494A (zh) 一种可视化造价数据预测智能分析***
CN113762591A (zh) 一种基于gru和多核svm对抗学习的短期电量预测方法及***
CN113033898A (zh) 基于k均值聚类与bi-lstm神经网络的电负荷预测方法及***
CN112232570A (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
GR01 Patent grant
GR01 Patent grant