CN107918368A - 钢铁企业煤气产生量与消耗量的动态预测方法及设备 - Google Patents
钢铁企业煤气产生量与消耗量的动态预测方法及设备 Download PDFInfo
- Publication number
- CN107918368A CN107918368A CN201711210058.9A CN201711210058A CN107918368A CN 107918368 A CN107918368 A CN 107918368A CN 201711210058 A CN201711210058 A CN 201711210058A CN 107918368 A CN107918368 A CN 107918368A
- Authority
- CN
- China
- Prior art keywords
- mtd
- msub
- mtr
- production
- gas
- 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
Links
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 36
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 25
- 239000010959 steel Substances 0.000 title claims abstract description 25
- 229910052742 iron Inorganic materials 0.000 title claims abstract description 18
- 239000003034 coal gas Substances 0.000 title abstract description 16
- 238000004519 manufacturing process Methods 0.000 claims abstract description 105
- 238000012549 training Methods 0.000 claims abstract description 53
- 230000008439 repair process Effects 0.000 claims abstract description 6
- 238000012423 maintenance Methods 0.000 claims description 15
- 238000010276 construction Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000035484 reaction time Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 2
- 238000009434 installation Methods 0.000 claims description 2
- 239000003245 coal Substances 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 10
- 239000000571 coke Substances 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 239000006227 byproduct Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 239000011541 reaction mixture Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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/41865—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 job scheduling, process planning, material flow
-
- 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/32252—Scheduling production, machining, job shop
-
- 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]
-
- 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/80—Management or planning
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种钢铁企业煤气产生量与消耗量的动态预测方法及设备，方法包括：获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应煤气产耗历史数据周期的生产计划和检修计划；根据生产计划和检修计划，确定所有设备生产工况；针对煤气产耗历史数据中的每一个设备的生产工况，获取每一设备的工况点，并根据所有设备的工况点对煤气量历史数据进行分类；依据分类后的各工况点的历史数据训练LSSVM模型；并采用训练的最优预测模型对煤气产耗历史数据的部分数据进行预测，得到输出的预测结果，输出的预测结果包括：未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。上述方法能够提高动态预测的煤气量的准确性。
Description
技术领域
本发明属于钢铁企业煤气预测技术，尤其涉及一种钢铁企业煤气产生量与消耗量的动态预测方法及设备。
背景技术
在我国，钢铁工业是国民经济的基础性支柱产业，是资源、能源密集型产业，也是节能减排工作的重点。钢铁企业高炉煤气、焦炉煤气和转炉煤气是企业重要的二次能源，占企业总能耗的30％左右。钢铁企业煤气系统的运行规律与设备、原材料和产品种类等密切相关，同时也受到生产计划、检修计划和设备运行工况等因素影响。煤气预测是煤气系统的重要组成部分，根据生产计划、检修计划等工况信息，对煤气的产生量和消耗量进行预测，以便进行后续的煤气资源分配与优化调度，为煤气系统的安全平稳运行，降低放散，提高煤气利用率提供依据。因此，煤气预测成为钢铁企业节能减排的重要环节之一。
钢铁企业副产煤气由高炉、焦炉、转炉生产过程中产生的，广泛应用于钢铁生产各个流程，富余的煤气还会被用来产生电力、蒸汽或其他副产品。煤气产生量是指高炉、焦炉、转炉产生的煤气量；煤气消耗量是指钢铁企业主要生产工序消耗的煤气量。实际生产中，由于钢铁生产的波动性，煤气的供需呈现动态变化，供需不平衡时会出现煤气放散或影响生产的情况。例如高炉检修结束后的复风阶段，煤气发生量较少而煤气消耗量正常，需要煤气柜供给煤气，如果此时气柜内煤气不足将会影响生产，导致严重后果；正常生产工况下如果煤气柜煤气较多，遇到煤气产生量较大的情况将不能起到平衡作用，致使煤气放散。为了合理调度煤气，调整煤气柜柜位，我们需要提前对煤气产生量和消耗量进行预测。如果能够提前预测煤气供需状况，将大幅度提高煤气的使用效率，减少不必要的放散，减少污染物的排放。
目前煤气动态预测方法大多是针对某一特定设备在特定工况条件下的动态预测，由于钢铁企业实际生产过程受到生产计划、检修计划等人为因素和炉况等不确定因素的影响，设备工况变动频繁，特定工况条件下的动态预测在设备工况变化时难以准确预测，导致煤气预测值准确性不高。因此提供一种针对多工况条件下钢铁企业煤气产生量与消耗量的动态预测方法显得极其重要。
发明内容
针对现有技术中的问题，本发明提供一种钢铁企业煤气产生量与消耗量的动态预测方法及设备。
第一方面，本发明提供一种钢铁企业煤气产生量与消耗量的动态预测方法，包括：
步骤A：获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应所述煤气产耗历史数据周期的生产计划和检修计划；
步骤B：根据所述生产计划和检修计划，确定所有设备生产工况；
步骤C：针对煤气产耗历史数据中的每一个设备的生产工况，获取每一设备的工况点，并根据所有设备的工况点对煤气量历史数据进行分类；
依据分类后的各工况点的历史数据训练LSSVM模型；并采用训练后的最优预测模型对煤气产耗历史数据的部分数据进行预测，得到输出的预测结果，所述输出的预测结果包括：未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。
可选地，所述步骤C包括：
C1、以集合Vi，j表示第i个设备第j个工况的煤气量历史数据集；
C2、根据集合Vi，j，通过小波分析方法和向量空间重构方法构造LSSVM模型的训练数据集；
C3、利用所述训练数据集，并通过启发式参数寻优方式对所述LSSVM模型进行训练，获得训练后的最优预测模型；
C4、依据待预测的设备、及待预测设备的工况和所述煤气产耗历史数据，构造待处理的预测数据集；
C5、采用最优预测模型处理所述待处理的预测数据集，获得输出的预测结果。
可选地，步骤C2包括：
依据企业煤气调度需求确定初始样本采样周期t、向量空间变换样本维数m以及模型中其它参数；
获取煤气产生量和消耗量V0，并基于小波分析方法，提取低频信号提取低频信号Di(n)；
对所述煤气产耗历史数据中各用户的煤气产耗历史数据进行归一化处理，使数据线性映射到0-1之间；
对归一化处理后的煤气产耗历史数据进行向量空间重构，将线性排列的数据构造为最优预测模型所需要的训练数据集，构造过程可表示为：
其中，输入数据中i时刻的煤气产耗历史数据记为xi，输出数据中的i时刻预测结果记为Yi，每一行中的一组xi和一组Yi分别为模型的一个输入样本和一个输出样本，训练数据集构造方法记为
利用输入输出样本xt-m+1 xt-m+2 … xt|Yt+1 Yt+2 … Yt+i对模型进行训练，以及应用启发式参数寻优，对LSSVM预测模型迭代求逆，由公式计算预测值
其中，I为单位矩阵，Ω为满足的Mercer定理核函数；
根据预测值的判断标准/约束条件，直至将所述训练数据集中的所有样本训练完，得到训练数据集。
可选地，所述待处理的预测数据集中各样本的训练方式和所述训练数据集中各样本的训练方式是相同的。
可选地，所述方法还包括：
将预测结果上传到企业的煤气调度服务器上，煤气调度系统根据一侧结果进行调度。
可选地，所述C2中的根据预测值的判断标准，获得训练数据集，包括：
在预测下一个样本之前需要计算第j个预测值误差平方：
其中，式中，为第j个预测值，Yj为第j个预测样本；
进而得出训练的修正值βj＝ηjej
将修正值加入下一时刻的样本学习预测。
可选地，所述C2中获取训练数据集时，当新的数据集(xi，yi)，i∈(t，t+N)进入训练数据集时，排除之前N个数据，重新得到输入样本。
可选地，所述输出的预测结果是反归一化之后的结果。
可选地，高炉设备的工况点包括：
正常生产、计划年修、计划月修、临时休风、热风炉减风1500km3/h、热风炉减风1000km3/h和热风炉减风1000km3/h以下。
第二方面，本发明还提供一种钢铁企业煤气产生量与消耗量的动态预测设备，包括：
处理器，用于存储所述处理器的可执行指令的存储器；
其中所述处理器被配置为：
获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应所述煤气产耗历史数据周期的生产计划和检修计划；
根据所述生产计划和检修计划，确定所有设备生产工况；
针对煤气产耗历史数据中的每一个设备的生产工况，获取每一设备的工况点，并根据所有设备的工况点对煤气量历史数据进行分类；
依据分类后的各工况点的历史数据训练LSSVM模型；并采用训练后的最优预测模型对煤气产耗历史数据的部分数据进行预测，得到输出的预测结果，所述输出的预测结果包括：未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。
本发明具有的有益效果如下：
本发明的钢铁企业煤气产生量与消耗量的动态预测方法及设备，通过获取每一设备的工况点，并采用工况点对煤气量历史数据进行分类，依据分类的工况点历史数据训练LSSVM模型，进而采用训练后的LSSVM模型(即最优预测模型)对构建的预测数据集进行预测，获得预测结果，预测结果相对于现有技术，预测准确性更高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。
图1为本发明一实施例提供的钢铁企业煤气产生量与消耗量的动态预测方法的流程图；
图2为本发明实施例中高炉工况点的分类示意图；
图3为本发明实施例中改进的LSSVM预测模型流程图；
图4为本发明实施例中小波分析示意图。
具体实施方式
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
在以下的描述中，将描述本发明的多个不同的方面，然而，对于本领域内的普通技术人员而言，可以仅仅利用本发明的一些或者全部结构或者流程来实施本发明。为了解释的明确性而言，阐述了特定的数目、配置和顺序，但是很明显，在没有这些特定细节的情况下也可以实施本发明。在其它情况下，为了不混淆本发明，对于一些众所周知的特征将不再进行详细阐述。
如图1所示，图1示出了本发明一实施例提供的钢铁企业煤气产生量与消耗量的动态预测方法的流程图，本实施例的动态预测方法包括下述步骤：
步骤S1：获取训练所需历史数据，通过从数据库中获得各设备煤气量历史数据，以及对应时期的生产计划和检修计划；通过生产检修计划确定设备生产工况。
具体地，上述步骤S1包括下述子步骤：
步骤101：通过企业调度应用服务器从综合数据集成平台服务器获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应时期的生产计划和检修计划；
步骤102：通过生产计划和检修计划确定设备生产工况信息。
步骤S2：按照设备生产工况，设立工况点，进行工况点划分，对设备生产工况具体细分；根据设备及工况对煤气量历史数据进行分类，以集合Vi，j表示第i个设备第j个工况煤气量历史数据集。
步骤S2可包括下述子步骤：
步骤201：设立工况点，按照设备生产工况进行工况点划分，对设备生产工况具体细分。
举例来说，如图2所示，高炉运行工况点划分为高炉基准工况、高炉年修、高炉月修、临时休风、热风炉减风1500km3/h、热风炉减风1000km3/h等。
步骤202：根据步骤102中的设备及生产工况信息对煤气量历史数据进行分类，以集合Vi，j表示第i个设备第j个工况煤气量历史数据集；
步骤S3：获取步骤202的煤气量历史数据集Vi，j，通过小波分析方法和向量空间重构方法构造模型的训练数据集；利用训练数据集对改进的LSSVM预测模型训练，通过启发式参数寻优，得到适应该设备该工况的最优预测模型，将第i个设备第j个工况的最优预测模型记为Li，j，相应的训练数据集构造方法记为对待预测的设备及工况依据提取所需历史数据，构造预测数据集并利用模型Li，j进行预测。
相应的训练数据集构造方法是指归一化方法和向量空间重构，需要记录的是构造过程中的一些参数。针对不同的设备和工况需要不同的参数，而对同一设备和工况的训练和预测必须使用相同的模型参数。
对待预测的设备及工况依据训练数据集构造方法提取所需历史数据，可举例说明：
对某一用户使用前一小时的历史数据预测未来一小时的数据，那么在预测8：00到9：00的数据时就需要7：00到8：00之间的数据。这一关系是可以调整的，可以由向量空间重构方法中的参数确定。
另外，如图3所示，图3示出了改进的LSSVM预测模型的预测流程。需要说明的是，图3中满足模型收敛条件，可理解为则本次模型训练成功，否则要继续训练。图3中显示的属于LSSVM模型训练的一部分，本实施例不对其限定，仅为参考说明。
具体地，步骤S3包括下述子步骤：
步骤301：通过企业调度应用服务器从综合数据集成平台服务器获取未来调度周期内的生产计划和检修计划，得到未来调度周期内各用户的生产工况；
步骤302：通过企业调度应用服务器从综合数据集成平台服务器获取对应生产工况的煤气产耗历史数据Vi，j；
步骤303：根据企业煤气调度的需求确定初始样本采样周期t、向量空间变换样本维数m以及模型中其它参数。
步骤304：获取煤气产生量和消耗量利用小波分析方法，提取低频“信号”Di(n)。
在步骤304中是对一个设备工况的所有数据进行处理的。
从图4(a)至图4(c)可以看出，小波分析方法的小波滤波器分离了趋势数据和波动数据，滤波器分解函数为
式(1)中，为尺度函数，cj，k为尺度系数；Ψm，k(t)为小波函数，dm，k为小波函数。
步骤305：对各用户的煤气产耗历史数据进行归一化处理，使数据线性映射到0-1之间，其转换函数表示为：
式(2)中，x表示原始煤气产耗历史数据，x*表示归一化处理后的煤气产耗历史数据，min表示该用户煤气产耗历史数据中的最小值，max表示该用户煤气产耗历史数据中的最大值。
本实施例中，归一化是对各用户的煤气产耗历史数据分组进行的，每一个用户进行一次归一化。
步骤306：对归一化处理后的煤气产耗历史数据进行向量空间重构，将线性排列的数据构造为改进的LSSVM预测模型所需要的输入输出数据，也即训练数据集，输入数据中i时刻的煤气产耗历史数据记为xi，输出数据中的i时刻预测结果记为Yi，构造多输入多输出数据的过程可表示为：
其中每一行中的一组xi和一组Yi分别为模型的一个输入样本和一个输出样本，相应的训练数据集构造方法记为
步骤307：使用步骤306中训练数据集的第一组(即第一行)输入输出样本xt-m+1xt-m+2 … xt|Yt+1 Yt+2 … Yt+i对模型进行训练。xt的下标t为该第一组输入输出样本对应的当前时刻。
进而应用启发式参数寻优，对LSSVM预测模型迭代求逆，计算、b，由公式计算，
公式(3)中I为单位矩阵，Ω为满足的Mercer定理核函数。
由公式(3)求出预测值：
本实施例中选取的核函数为高斯径向机(RBF)核函数，其表达式为
步骤308：计算第j个预测值误差平方；
式(5)中，为第j个预测值，Yj为第j个预测样本。将全部i个预测误差的集合{et+1，et+2，…，et+i}记为E(即下述步骤309中的预测误差)。
得出修正值βj＝ηjej，修正值是LSSVM预测模型中的一个参数。
加入下一时刻的样本学习预测，这里的下一时刻的样本即第t+i+1个样本或者n个时段之后也就是第t+i+n个样本。
该步骤308可理解为预测值的判断标准/约束条件。预测误差E可理解为步骤308中计算的误差的集合{et+1，et+2，…，et+i}。
步骤309：转至步骤307，将步骤308中的预测误差E加入LSSVM模型优化目标的约束条件中，使用步骤306中训练数据集的下一组输入输出样本xt-m+2 xt-m+3 … xt+1|Yt+2 Yt+3… Yt+i+1进行训练，直至训练数据集中的最后一组样本。
输入样本是模型的预测依据，例如在预测8：00到9：00的数据时就需要7：00到8：00之间的数据，这个7：00到8：00之间的数据就是输入样本。当需要进行新的预测时，就需要新的对应这段需要预测时间的历史数据。这里写的这段时间是(t,t+N)。
步骤310：根据待预测时间段(即未来调度周期)确定预测所需的历史数据时间段，对该段历史数据按照步骤305和步骤306的方法做出相同的处理，得到预测数据集；
步骤311：将预测数据集代入改进的LSSVM预测模型中，得到输出数据，对输出数据进行反归一化得到预测数据，即未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量；反归一化转换函数为：
y＝y*×(max-min)+min (6)
式(6)中，y表示未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量，y*表示改进的LSSVM预测模型的输出数据，min表示该用户煤气产耗历史数据中的最小值，max表示该用户煤气产耗历史数据中的最大值，其中min和max与步骤305中的值相等。
步骤S4：将预测结果上传到企业的煤气调度服务器上，煤气调度系统依据预测结果进行调度，企业能源管理人员依据煤气预测结果与煤气调度方案进行调度。
本发明综合考虑了钢铁企业设备工况不同对煤气产生和消耗的影响，建立了基于设备多工况运行的煤气产生量与消耗量动态预测方法来解决本发明所提到的复杂工况条件下的煤气预测问题；
采用本发明提出的钢铁企业煤气产生量与消耗量的动态预测方法，可以有效的预测未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。
最后应说明的是：以上所述的各实施例仅用于说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或全部技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种钢铁企业煤气产生量与消耗量的动态预测方法，其特征在于，包括：
步骤A：获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应所述煤气产耗历史数据周期的生产计划和检修计划；
步骤B：根据所述生产计划和检修计划，确定所有设备生产工况；
步骤C：针对煤气产耗历史数据中的每一个设备的生产工况，获取每一设备的工况点，并根据所有设备的工况点对煤气量历史数据进行分类；
依据分类后的各工况点的历史数据训练LSSVM模型；并采用训练后的最优预测模型对煤气产耗历史数据的部分数据进行预测，得到输出的预测结果，所述输出的预测结果包括：未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。
2.根据权利要求1所述的方法，其特征在于，所述步骤C包括：
C1、以集合Vi，j表示第i个设备第j个工况的煤气量历史数据集；
C2、根据集合Vi，j，通过小波分析方法和向量空间重构方法构造LSSVM模型的训练数据集；
C3、利用所述训练数据集，并通过启发式参数寻优方式对所述LSSVM模型进行训练，获得训练后的最优预测模型；
C4、依据待预测的设备、及待预测设备的工况和所述煤气产耗历史数据，构造待处理的预测数据集；
C5、采用最优预测模型处理所述待处理的预测数据集，获得输出的预测结果。
3.根据权利要求2所述的方法，其特征在于，步骤C2包括：
依据企业煤气调度需求确定初始样本采样周期t、向量空间变换样本维数m以及模型中其它参数；
获取煤气产生量和消耗量V0，并基于小波分析方法，提取低频信号提取低频信号Di(n)；
对所述煤气产耗历史数据中各用户的煤气产耗历史数据进行归一化处理，使数据线性映射到0-1之间；
对归一化处理后的煤气产耗历史数据进行向量空间重构，将线性排列的数据构造为最优预测模型所需要的训练数据集，构造过程可表示为：
<mrow>
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>4</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>5</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&RightArrow;</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>4</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>5</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>6</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mn>7</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>7</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>6</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>5</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>4</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>x</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
<mo>|</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mn>6</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>7</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>8</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mn>7</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>8</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>9</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mn>8</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>9</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mn>10</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Y</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Y</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mrow>
其中，输入数据中i时刻的煤气产耗历史数据记为xi，输出数据中的i时刻预测结果记为Yi，每一行中的一组xi和一组Yi分别为模型的一个输入样本和一个输出样本，训练数据集构造方法记为
利用输入输出样本xt-m+1 xt-m+2 … xt|Yt+1 Yt+2 … Yt+i对模型进行训练，以及应用启发式参数寻优，对LSSVM预测模型迭代求逆，由公式计算预测值
其中，I为单位矩阵，Ω为满足的Mercer定理核函数；
根据预测值的判断标准/约束条件，直至将所述训练数据集中的所有样本训练完，得到训练数据集。
4.根据权利要求3所述的方法，其特征在于，所述待处理的预测数据集中各样本的训练方式和所述训练数据集中各样本的训练方式是相同的。
5.根据权利要求1所述的方法，其特征在于，所述方法还包括：
将预测结果上传到企业的煤气调度服务器上，煤气调度系统根据一侧结果进行调度。
6.根据权利要求3所述的方法，其特征在于，所述C2中的根据预测值的判断标准，获得训练数据集，包括：
在预测下一个样本之前需要计算第j个预测值误差平方：
<mrow>
<msub>
<mi>e</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>Y</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mi>i</mi>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
其中，式中，为第j个预测值，Yj为第j个预测样本；
进而得出训练的修正值βj＝ηjej；
将修正值加入下一时刻的样本学习预测。
7.根据权利要求6所述的方法，其特征在于，所述C2中获取训练数据集时，当新的数据集(xi，yi)，i∈(t，t+N)进入训练数据集时，排除之前N个数据，重新得到输入样本。
8.根据权利要求7所述的方法，其特征在于，所述输出的预测结果是反归一化之后的结果。
9.根据权利要求1至8任一所述的方法，其特征在于，高炉设备的工况点包括：
正常生产、计划年修、计划月修、临时休风、热风炉减风1500km3/h、热风炉减风1000km3/h和热风炉减风1000km3/h以下。
10.一种钢铁企业煤气产生量与消耗量的动态预测设备，其特征在于，包括：
处理器，用于存储所述处理器的可执行指令的存储器；
其中所述处理器被配置为：
获取未来调度周期内的生产工况下煤气产耗历史数据，以及对应所述煤气产耗历史数据周期的生产计划和检修计划；
根据所述生产计划和检修计划，确定所有设备生产工况；
针对煤气产耗历史数据中的每一个设备的生产工况，获取每一设备的工况点，并根据所有设备的工况点对煤气量历史数据进行分类；
依据分类后的各工况点的历史数据训练LSSVM模型；并采用训练后的最优预测模型对煤气产耗历史数据的部分数据进行预测，得到输出的预测结果，所述输出的预测结果包括：未来调度周期内煤气产生装置的煤气发生量和各生产用户的煤气消耗量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711210058.9A CN107918368B (zh) | 2017-11-28 | 2017-11-28 | 钢铁企业煤气产生量与消耗量的动态预测方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711210058.9A CN107918368B (zh) | 2017-11-28 | 2017-11-28 | 钢铁企业煤气产生量与消耗量的动态预测方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107918368A true CN107918368A (zh) | 2018-04-17 |
| CN107918368B CN107918368B (zh) | 2019-07-26 |
Family
ID=61898054
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711210058.9A Active CN107918368B (zh) | 2017-11-28 | 2017-11-28 | 钢铁企业煤气产生量与消耗量的动态预测方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107918368B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109086949A (zh) * | 2018-09-20 | 2018-12-25 | 鞍钢集团自动化有限公司 | 基于煤气成份变化的高炉煤气发生量及其热值预测方法 |
| CN109948858A (zh) * | 2019-03-25 | 2019-06-28 | 东北大学 | 基于数据驱动的工况识别方法 |
| CN110428053A (zh) * | 2019-08-08 | 2019-11-08 | 东北大学 | 一种蒸汽产耗量的动态预测方法 |
| CN112966861A (zh) * | 2021-02-24 | 2021-06-15 | 浙江中控技术股份有限公司 | 一种基于案例匹配的产耗平衡测算方法 |
| CN113177369A (zh) * | 2021-06-15 | 2021-07-27 | 中冶赛迪技术研究中心有限公司 | 一种能源调度评估方法及系统 |
| CN116596231A (zh) * | 2023-05-09 | 2023-08-15 | 重庆赛迪奇智人工智能科技有限公司 | 一种铁水调度计划修正方法、装置、设备及存储介质 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070136120A1 (en) * | 2005-11-28 | 2007-06-14 | Tohru Watanabe | System and method for providing service |
| CN101109952A (zh) * | 2007-08-23 | 2008-01-23 | 济南钢铁集团总公司 | 基于柜位预测的钢铁企业煤气动态平衡实时控制方法 |
| CN102147273A (zh) * | 2010-01-29 | 2011-08-10 | 大连理工大学 | 一种基于数据的冶金企业高炉煤气动态预测方法 |
| US20140012498A1 (en) * | 2012-03-19 | 2014-01-09 | Thomas F. Gustafson | Informed Traveler Program and Application |
| CN103793767A (zh) * | 2014-02-26 | 2014-05-14 | 大连理工大学 | 基于炼钢节奏估计的冶金企业转炉煤气发生量长期预测方法 |
| CN104102212A (zh) * | 2014-07-29 | 2014-10-15 | 浙江中控软件技术有限公司 | 一种钢铁企业中煤气和蒸汽系统的调度方法、设备和系统 |
| CN104238483A (zh) * | 2013-12-23 | 2014-12-24 | 东北大学 | 基于生产工况的高炉煤气产生量与消耗量的标定方法 |
| CN105303243A (zh) * | 2014-06-09 | 2016-02-03 | 鞍钢股份有限公司 | 一种高炉煤气预测平衡调度系统及方法 |
-
2017
- 2017-11-28 CN CN201711210058.9A patent/CN107918368B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070136120A1 (en) * | 2005-11-28 | 2007-06-14 | Tohru Watanabe | System and method for providing service |
| CN101109952A (zh) * | 2007-08-23 | 2008-01-23 | 济南钢铁集团总公司 | 基于柜位预测的钢铁企业煤气动态平衡实时控制方法 |
| CN102147273A (zh) * | 2010-01-29 | 2011-08-10 | 大连理工大学 | 一种基于数据的冶金企业高炉煤气动态预测方法 |
| US20140012498A1 (en) * | 2012-03-19 | 2014-01-09 | Thomas F. Gustafson | Informed Traveler Program and Application |
| CN104238483A (zh) * | 2013-12-23 | 2014-12-24 | 东北大学 | 基于生产工况的高炉煤气产生量与消耗量的标定方法 |
| CN103793767A (zh) * | 2014-02-26 | 2014-05-14 | 大连理工大学 | 基于炼钢节奏估计的冶金企业转炉煤气发生量长期预测方法 |
| CN105303243A (zh) * | 2014-06-09 | 2016-02-03 | 鞍钢股份有限公司 | 一种高炉煤气预测平衡调度系统及方法 |
| CN104102212A (zh) * | 2014-07-29 | 2014-10-15 | 浙江中控软件技术有限公司 | 一种钢铁企业中煤气和蒸汽系统的调度方法、设备和系统 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109086949A (zh) * | 2018-09-20 | 2018-12-25 | 鞍钢集团自动化有限公司 | 基于煤气成份变化的高炉煤气发生量及其热值预测方法 |
| CN109086949B (zh) * | 2018-09-20 | 2021-09-21 | 鞍钢集团自动化有限公司 | 基于煤气成份变化的高炉煤气发生量及其热值预测方法 |
| CN109948858A (zh) * | 2019-03-25 | 2019-06-28 | 东北大学 | 基于数据驱动的工况识别方法 |
| CN110428053A (zh) * | 2019-08-08 | 2019-11-08 | 东北大学 | 一种蒸汽产耗量的动态预测方法 |
| CN112966861A (zh) * | 2021-02-24 | 2021-06-15 | 浙江中控技术股份有限公司 | 一种基于案例匹配的产耗平衡测算方法 |
| CN113177369A (zh) * | 2021-06-15 | 2021-07-27 | 中冶赛迪技术研究中心有限公司 | 一种能源调度评估方法及系统 |
| CN113177369B (zh) * | 2021-06-15 | 2024-03-01 | 中冶赛迪技术研究中心有限公司 | 一种能源调度评估方法及系统 |
| CN116596231A (zh) * | 2023-05-09 | 2023-08-15 | 重庆赛迪奇智人工智能科技有限公司 | 一种铁水调度计划修正方法、装置、设备及存储介质 |
| CN116596231B (zh) * | 2023-05-09 | 2024-02-06 | 重庆赛迪奇智人工智能科技有限公司 | 一种铁水调度计划修正方法、装置、设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107918368B (zh) | 2019-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107918368B (zh) | 钢铁企业煤气产生量与消耗量的动态预测方法及设备 | |
| CN111353656B (zh) | 一种基于生产计划的钢铁企业氧气负荷预测方法 | |
| Abrell et al. | How effective is carbon pricing?—A machine learning approach to policy evaluation | |
| CN108009667A (zh) | 一种能源需求总量及结构预测系统 | |
| CN104181900B (zh) | 一种多能源介质分层动态调控方法 | |
| US20200219027A1 (en) | Knowledge transfer-based modeling method for blast furnace gas scheduling systems | |
| Xie et al. | Robust stochastic configuration network multi-output modeling of molten iron quality in blast furnace ironmaking | |
| Canyurt et al. | Three different applications of genetic algorithm (GA) search techniques on oil demand estimation | |
| CN106779384B (zh) | 一种基于信息粒度最优分配的钢铁工业高炉煤气长期区间预测方法 | |
| CN109559154A (zh) | 基于Johnson-Copula模型的电力需求总量相关因素分析系统及方法 | |
| CN110880044B (zh) | 一种基于马尔科夫链的负荷预测方法 | |
| Matino et al. | Application of echo state neural networks to forecast blast furnace gas production: pave the way to off-gas optimized management | |
| CN112541625B (zh) | 一种自适应多工况钢铁二次能源发生量动态预测方法 | |
| CN107706938A (zh) | 一种基于分位点回归的风电功率波动区间分析方法 | |
| Li et al. | Influential Factor Analysis and Projection of Industrial CO 2 Emissions in China Based on Extreme Learning Machine Improved by Genetic Algorithm. | |
| CN117113022A (zh) | 一种区域碳排放计算方法、系统、存储介质及设备 | |
| CN116127345A (zh) | 基于深度聚类生成对抗网络的转炉炼钢工艺模式设计方法 | |
| Jin et al. | A causal model-based scheduling approach for coke oven gas system in steel industry | |
| CN114759579A (zh) | 一种基于数据驱动的电网有功优化控制系统、方法和介质 | |
| CN114859718A (zh) | 电氢氨综合能源系统的控制方法、控制系统和计算机可读介质 | |
| CN114201912A (zh) | 一种成本驱动的可再生能源发电功率区间预测方法 | |
| CN112036713A (zh) | 一种基于主成分分析的煤改电综合效益评价方法和系统 | |
| CN105528681B (zh) | 一种基于隐树模型的冶金企业副产能源系统实时调整方法 | |
| CN116050615B (zh) | 一种基于超混沌储备池集成计算的钢铁企业副产煤气系统流量预测方法 | |
| CN114912335B (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 |