技术领域
[0001] 本
发明涉及煅烧过程优化方法和系统，尤其涉及煅烧过程
煤耗优化方法和系统。
背景技术
[0002] 新型干法
水泥窑炉煅烧过程主要包括
分解炉、
回转窑和篦冷机三个装置。
水泥窑炉煅烧过程主要流程如下。粉磨后的
生料从悬浮预热器顶部进入预热器，生料由于重
力作用下落，热
风进行充分热交换达到脱水预热的目的。预热后的生料进入分解炉，约90％的
碳酸盐在分解炉进行煅烧分解；分解后的生料由分解炉进入回转窑进行固相、液相反应，最终形成
水泥熟料。高温熟料进入篦冷机进行冷却、
破碎，同时换热后的热风可用于回转窑二次风、分解炉三次风以及煤磨用风、余热发电，提高
能量利用效率。
[0003] 水泥窑炉煅烧过程具有高能耗、高污染的特点，国家层面高度重视节能降耗，因此，水泥窑炉煅烧过程迫切需要有效的优化方法，起到节能降耗的作用，提高企业利润。
发明内容
[0004] 以下给出一个或多个方面的简要概述以提供对这些方面的基本理解。此概述不是所有构想到的方面的详尽综览，并且既非旨在指认出所有方面的关键性或决定性要素亦非试图界定任何或所有方面的范围。其唯一的目的是要以简化形式给出一个或多个方面的一些概念以为稍后给出的更加详细的描述之序。
[0005] 本发明的目的在于提供一种煅烧过程煤耗优化方法。该方法适用于包括但不限于水泥窑炉煅烧过程节能降耗。
[0006] 本发明第一方面提供一种煅烧过程煤耗优化方法，所述煅烧过程包括分解炉、回转窑和篦冷机装置，所述方法包括：
[0007] (1)采集物料参数和系统参数，
[0008] (2)利用
质量守恒、能量守恒和/或最小二乘法对分解炉、回转窑和篦冷机装置构建数学模型，对回转窑喂煤量、分解炉喂煤量、单位时间进入分解炉生料量进行预处理，[0009] (3)以低煤耗为优化目标，控制回路设定值作为决策变量，利用约束条件构建优化模型，所述控制回路设定值包括分解炉出口
温度、回转窑主
电机电流和篦冷机篦下压力；
[0010] (4)根据约束条件获得控制回路设定值。
[0011] 在一个或多个实施方案中，物料参数包括分解炉喂煤量CoalCal、进料量FRaw、窑尾烟室温度TTail、三次风温度TThird、回转窑喂煤量CoalKiln。
[0012] 在一个或多个实施方案中，系统参数包括高温风机转速SFan。
[0013] 在一个或多个实施方案中，步骤(2)中分解炉的模型结构如下：
[0014] TCal＝f(CoalCal,FRaw,TTail,TThird,SFan)
[0015] 其中模型输出TCal为分解炉出口温度，模型输入CoalCal为分解炉喂煤量，FRaw为进料量，TTail为窑尾烟室温度，TThird为三次风温度，SFan为高温风机转速。
[0016] 在一个或多个实施方案中，分解炉模型基于分解炉内的质量守恒和能量守恒构建。
[0017] 在一个或多个实施方案中，分解炉内包含生料分解反应和煤燃烧反应。
[0018] 在一个或多个实施方案中，生料分解反应是碳酸
钙分解反应。
[0019] 在一个或多个实施方案中，生料分解反应的吸热QR_Caco3为
[0020] QR_Caco3＝HRaw*RRaw
[0021] 其中HRaw为每千克碳酸钙分
解吸热量，RRaw为碳酸钙反应速率，与分解炉下料管温度和分解炉内二
氧化碳分压相关。
[0022] 在一个或多个实施方案中，碳酸钙反应速率RRaw为
[0023]
[0024] 其中为ACaco3反应界面的表面积，MCaco3为碳酸钙分子量，PCo2为分解炉中CO2分压，Pe为生料温度为TCal时CO2平衡分解压力，其中ks为
[0025]
[0026] 其中Pe为
[0027]
[0028] 在一个或多个实施方案中，分解炉中CO2分压与风机转速SFan相关。
[0029] 在一个或多个实施方案中，煤燃烧反应的放热QR_Coal为
[0030] QR_Coal＝HCoal*RCoal*ηCoal
[0031] 其中HCoal为
煤粉的
燃烧热，RCoal为煤反应速率，ηCoal为分解炉内煤粉燃尽率。
[0032] 在一个或多个实施方案中，煤反应速率RCoal与分解炉温度、分解炉内O2含量相关。
[0033] 在一个或多个实施方案中，煤反应速率RCoal为
[0034]
[0035] 其中，kc为反应速率常数，ACoal为煤燃烧反应表面积，PO2为氧气分压，氧气分压与窑尾高温风机转速SFan关联。
[0036] 在一个或多个实施方案中，分解炉内质量守恒如下：
[0037] MThird+MTail+MClear+MRaw+MCoal＝MGout+MSout+MError
[0038] 其中MThird为三次风流量，MTail为窑尾烟室风流量，MClear为入炉一次风流量，MRaw入炉物料流量，MCoal入炉煤粉流量。MGout和MSout分别为分解炉出口气体流量和入窑物料流量，MError为质量守恒补偿项，基本为0。
[0039] 在一个或多个实施方案中，步骤(2)中分解炉的模型包括如下能量守恒模型：
[0040] QThird+QTail+QClear+QRaw+QCoal+QR_Coal-QR_Caco3＝QGout+QSout+QError[0041] 其中QThird为三次风
显热，QTail为二次风显热，QClear为入炉一次风显热，QRaw为入炉物料显热，QCoal为入炉煤粉显热，QR_Coal为煤燃烧放热，QR_Caco3为碳酸钙分解吸热，QGout为分解炉出口气体显热，QSout为入窑物料显热，QError为补偿项。
[0042] 在一个或多个实施方案中，QSout＝Cs*MSout*TCal，QError＝Cs*MError*TCal。
[0043] 在一个或多个实施方案中，求解能量守恒模型，即可获得TCal。
[0044] 在一个或多个实施方案中，步骤(2)中回转窑模型结构如下：
[0045] IKiln＝f(TCal,CoalKiln,TThird)
[0046] 其中回转窑主电机电流Ikiln为模型输出，TCal为分解炉出口温度，TThird为三次风温度，CoalKiln为回转窑喂煤量，Coalcal为分解炉喂煤量。
[0047] 在一个或多个实施方案中，采用最小二乘法构建如下回转窑模型：
[0048] IKiln＝A1*TCal+A2*CoalKiln+A3*TThird+A4
[0049] 其中A1,A2,A3,A4为模型参数，参数采用最小二乘法进行辨识。相比于现有模型回转窑机理模型中气固反应反应动力学参数难确定，回转窑内部无检测情况，本发明采用基于输入输出数据的数据驱动建模方法，只需要输入输出数据即可获得精确优化模型。
[0050] 在一个或多个实施方案中，步骤(2)中篦冷机模型结构如下：
[0051] TThird＝f(PCool,IKiln,SFan)
[0052] 其中三次风温度TThird为模型输出，PCool为篦冷机篦下压力，IKiln为回转窑主电机电流，SFan为高温风机转速。
[0053] 在一个或多个实施方案中，篦冷机主体部分为一段篦冷机。
[0054] 在一个或多个实施方案中，一段篦冷机的气体热量守恒模型如下：
[0055] QG,In+QC,G＝QG,Out，
[0056] 且，一段篦冷机的固体能量守恒模型如下：
[0057] QC,In＝QC,G+QC,Out，
[0058] 其中QG,In和QG,Out分别为进入一段篦冷机气体和离开一段篦冷机气体热量，QC,In和QC,out分别为进入一段篦冷机熟料和离开一段篦冷机一段熟料的热量，其中QC,G为熟料和气体间换热。
[0059] 在一个或多个实施方案中，进入一段篦冷机气体热量为
[0060] QG,In＝CG,In*MG,In*TG,In
[0061] MG,In为篦冷机底层风机吹入篦冷机气体质量流量。
[0062] 在一个或多个实施方案中，离开一段篦冷机气体热量为
[0063] QG,Out＝CThird*MThird*TThird+CTail*MTail*TTail
[0064] 离开一段篦冷机气体热量是二次风和三次风热量，与风机转速SFan关联。
[0065] 在一个或多个实施方案中，
[0066] QC,In＝CC,In*MC*TC,In
[0067] 其中TC,In为进入篦冷机熟料温度。
[0068] 在一个或多个实施方案中，进入篦冷机熟料温度为
[0069]
[0070] 其中TBase,IBase为基准温度值和基准电流值。
[0071] 在一个或多个实施方案中，
[0072] QC,Out＝CC,Out*MC*TC,Out
[0073] 其中TC,Out为离开篦冷机一段熟料温度。
[0074] 在一个或多个实施方案中，熟料与气体之间换热QC,G为
[0075] QC,G＝K*a*(Tc-Tg)*VCooler
[0076] 其中K为气体间的综合换热系数，a为单位体积床层的有效表面积，(Tc-Tg)为气固之间平均温差，VCooler为一段篦冷机熟料的体积，如下所示
[0077]
[0078]
[0079] (Tc-Tg)＝(TC,In+TC,out-TG,In-TG,out)
[0080] VCooler＝LCooler*WCooler*HCooler
[0081] 其中h为气体
对流换热系数，λc为熟料导热系数，Φ为颗粒形状校正系数，x为颗粒的透热深度。其中ε为
空隙率，dp为物料尺寸，TC,In,TC,Out,TG,In,TG,out分别为进入一段篦冷机物料温度，离开一段篦冷机物料温度，进入一段篦冷机气体温度，离开一段篦冷机气体温度，LCooler,WCooler,HCooler分别为一段篦冷机长度、宽度和高度。
[0082] 在一个或多个实施方案中，采用篦冷机篦下压力PCool表征料层高度HCooler，优选为[0083] HCooler＝B1*PCool2+B2*PCool+B3
[0084] 其中B1,B2,B3为模型参数。
[0085] 在一个或多个实施方案中，步骤(3)中优化目标为单产熟料煤耗最小
[0086] min Energy＝(CoalCal+CoalKiln)/Clinker
[0087] 其中Clinker为熟料产量。
[0088] 在一个或多个实施方案中，步骤(3)中决策变量底层控制回路设定值。优选地，步骤(3)中决策变量包括分解炉出口温度，回转窑主电机电流，篦冷机篦下压力。
[0089] 在一个或多个实施方案中，步骤(3)中约束条件包括关键质量指标游离
氧化钙含量Fcao约束，分解炉预热器出口氧气含量O2约束，入窑生料分解率Rate约束，和/或决策变量TCal、IKiln及PCool约束。
[0090] 在一个或多个实施方案中，游离氧化钙含量Fcao模型结构如下所示：
[0091] Fcao＝f(TCal,IKiln,PCool)
[0092] 游离氧化钙模型采用辅助变量为分解炉出口温度TCal，回转窑主电机电流IKiln和篦冷机篦下压力PCool。
[0093] 在一个或多个实施方案中，分解炉预热器出口氧气含量O2模型结构如下所示：
[0094] O2＝f(CoalCal,SFan)。
[0095] 在一个或多个实施方案中，分解炉预热器出口O2含量为：
[0096] O2＝O2Third+O2Tail-O2CalCoal
[0097] 其中O2Third主要成分是空气，O2Tail成分与回转窑喂煤量CoalKiln有关，O2CalCoal成分与分解炉喂煤量CoalCal有关。
[0098] 在一个或多个实施方案中，
[0099]
[0100]
[0101] 其中kCarbon为煤粉含碳量，MO2为O2摩尔质量，MC为
焦炭摩尔质量。
[0102] 在一个或多个实施方案中，入窑生料分解率Rate模型结构如下所示：
[0103] Rate＝f(TCal,FRaw)
[0104] 在一个或多个实施方案中，入窑生料分解率Rate为
[0105] Rate＝RRaw/FRaw
[0106] 其中RRaw为生料分解速率，FRaw为单位时间进入分解炉生料。
[0107] 决策变量TCal、IKiln及PCool约束
[0108] 根据本发明的新型干法水泥窑炉煅烧过程全流程优化方法的一
实施例，步骤(4)中使用模式搜索
算法求解。
[0109] 在一个或多个实施方案中，决策变量约束选自以下一个或多个：分解炉出口温度：700℃≤T_Cal≤1000℃，回转窑主电机电流：600A≤I_Kiln≤1000A，篦冷机篦下压力：2Kpa≤P_Cool≤6Kpa。
[0110] 本发明还公开了一种煅烧过程参数调节方法，所述煅烧过程包括使用分解炉、回转窑和篦冷机装置，所述方法包括：
[0111] (1)采集物料参数和系统参数，
[0112] (2)采用根据
权利要求1-6中任一项所述的方法构建的模型对煅烧过程的控制回路设定值进行优化，
[0113] (3)根据优化的控制回路设定值调节煅烧过程参数。
[0114] 在一个或多个实施方案中，物料参数包括分解炉喂煤量、进料量、窑尾烟室温度、三次风温度、回转窑喂煤量，系统参数包括高温风机转速S_Fan。
[0115] 本发明还公开了一种煅烧过程煤耗优化系统，包括：
[0116]
数据采集模
块，采集物料参数和系统参数，
[0117] 模型构建模块，利用质量守恒、能量守恒和/或最小二乘法对分解炉、回转窑和篦冷机装置构建数学模型，对回转窑喂煤量、分解炉喂煤量、单位时间进入分解炉生料量进行预处理，
[0118] 优化模型构建模块，以低煤耗为优化目标，控制回路设定值作为决策变量，利用约束条件构建优化模型，所述控制回路设定值包括分解炉出口温度、回转窑主电机电流和篦冷机篦下压力，
[0119]
数据处理模块，根据约束条件获得控制回路设定值。
[0120] 在一个或多个实施方案中，物料参数包括分解炉喂煤量CoalCal、进料量FRaw、窑尾烟室温度TTail、三次风温度TThird、回转窑喂煤量CoalKiln。在一个或多个实施方案中，系统参数包括高温风机转速SFan。
[0121] 本发明还公开了一种煅烧过程煤耗优化系统，包括计算机以及运行于计算机上的一
计算机程序，计算机程序在计算机上运行如上所述的方法。
[0122] 本发明还公开了一种存储计算机程序的计算机可读存储介质，其特征在于，存储介质上所存储的计算机程序运行后执行如上所述的方法。
[0123] 本发明对比
现有技术有如下的有益效果：本发明的煅烧过程煤耗优化方法能够通过煅烧全流程过程机理和过程数据建立优化模型，构建能耗最低优化模型，采用
优化算法获得优化设定值。
[0124] 本发明的主要优点在于：
[0125] 1、相比于传统子单元优化，本发明综合考虑了分解炉、回转窑和篦冷机三个单元全流程优化，即窑炉煅烧过程全流程优化；针对传统回转窑机理模型参数难以获得的问题，本发明采用基于数据驱动建模方法获得有效的模型。
[0126] 2、相比于设计优化，本发明无需改变任何装置结构，只需
修改操作参数即可达到节能降耗的目的，具有成本低、收益好的特点。
[0127] 总的来说，本发明有效
结合水泥窑炉煅烧过程工艺和过程数据，构建优化模型充分考虑各类约束，获得底层回路优化设定值。
附图说明
[0128] 在结合以下附图阅读本公开的实施例的详细描述之后，能够更好地理解本发明的上述特征和优点。在附图中，各组件不一定是按比例绘制，并且具有类似的相关特性或特征的组件可能具有相同或相近的附图标记。
[0129] 图1示出了本发明煅烧过程煤耗优化方法的
流程图。
[0130] 图2示出了本发明一实例中水泥窑炉煅烧过程流程图及优化变量。
[0131] 图3示出了分解炉出口温度模型值与实际值对比的示意图。
[0132] 图4示出了回转窑主电机电流模型值与实际值对比的示意图。
[0133] 图5示出了三次风温度模型值与实际值对比的示意图。
具体实施方式
[0134] 以下结合附图和具体实施例对本发明作详细描述。注意，以下结合附图和具体实施例描述的诸方面仅是示例性的，而不应被理解为对本发明的保护范围进行任何限制。
[0135] 图1示出了本发明本发明煅烧过程优化方法的一个实施方案的流程。该实施方案具体以水泥窖炉煅烧过程为例进行说明，但是本发明并不局限于水泥窖炉煅烧，而是能扩展到任何与之类似的场景中。图2示出了水泥窑炉煅烧流程，主要包括预热器、分解炉、回转窑和篦冷机，同时示出了本发明的优化决策变量：分解炉出口温度TCal，回转窑主电机电流IKiln，篦冷机篦下压力PCool。分解炉
燃料用量的比例对整个分解系统的热力分布有着重要影响，而分解炉的燃料用量又与分解炉
温度控制有关；主电机电流是反映回转窑烧成状况的关键变量；篦下压力反映篦冷机料层厚度，影响篦冷机换热效果。实施方案数据全部来源于实际现场。
[0136] 本发明中，煤耗表示煅烧过程中所消耗的煤与产物量的比值。在本文中，煤耗包括单产熟料煤耗。在一个或多个实施方案中，煤耗是分解炉用煤量和回转窑用煤量与熟料产量的比值。
[0137] 步骤S1：建立全流程优化模型
[0138] 本步骤是对煅烧过程建模的过程，包括针对分解炉、回转窑和篦冷机，利用质量守恒、能量守恒结合最小二乘法等分别得到模型。
[0139] (1)分解炉模型：
[0140] 1)分解炉机理模型结构如下：
[0141] TCal＝f(CoalCal,FRaw,TTail,TThird,SFan)
[0142] 其中模型输出TCal为分解炉出口温度，模型输入CoalCal为分解炉喂煤量，FRaw为进料量，TTail为窑尾烟室温度，TThird为三次风温度，SFan为高温风机转速。分解炉出口温度模型基于分解炉质量守恒和能量守恒，主要包含过程反应为a)生料分解反应和b)煤燃烧反应两部分。
[0143] a)生料分解反应主要是生料中的碳酸钙反应，碳酸钙反应吸热QR_Caco3[0144] QR_Caco3＝HRaw*RRaw
[0145] 其中HRaw为每千克碳酸钙分解吸热量，在一个或多个实施方案中，HRaw＝1660kj/kg；RRaw为碳酸钙反应速率。碳酸钙反应速率RRaw，主要与分解炉下料管温度和分解炉内二氧化碳分压相关。
[0146]
[0147] 其中ACaco3为反应界面的表面积，与颗粒细度d关系为ACaco3＝π*d2。颗粒细度可为10-150um，优选为50-80um。在一个或多个实施方案中，颗粒细度d＝70um，MCaco3＝100g/mol。
PCo2为分 解炉中CO2分压 ，CO2分压与高温风机转速SFan成比 例关系
其中SFan_base＝600rmp,PCo2_base＝0.36Kpa。Pe为生料温
度为TCal时CO2平衡分解压力。
[0148] 其中ks为
[0149]
[0150] 其中Er可为1.6*105-2*105kj/mol，优选1.8*105-1.9*105kj/mol。在一个或多个实施方案中，Er＝1.85*105kj/mol，常数R＝8.314J/(mol·K)。
[0151] 其中Pe为
[0152]
[0153] b)煤燃烧放热QR_Coal
[0154] QR_Coal＝HCoal*RCoal*ηCoal
[0155] 其中HCoal为煤粉的燃烧热，RCoal为煤燃烧反应速率，ηCoal是分解炉内煤粉燃尽率，为80％-95％，优选85-92％。在一个或多个实施方案中，ηCoal＝90％。煤燃烧反应速率RCoal与分解炉温度、分解炉内O2含量相关。
[0156]
[0157] 其中，kc为反应速率常数。ACoal为煤燃烧反应表面积，与颗粒细度d关系为ACoal＝π*d2。PO2为氧气分压，氧气分压与窑尾高温风机转速SFan成比例关系，其中 Po2_base＝0.5Kpa。Ec可为1.2*104-1.4*104kj/mol，优选1.3*104-
1.35*104kj/mol在一个或多个实施方案中，kc＝1.32*105mol/(m2·s)，EC＝6.81*104kj/mol。
[0158] 2)分解炉内质量守恒如下：
[0159] MThird+MTail+MClear+MRaw+MCoal＝MGout+MSout+MError
[0160] 其中MThird为三次风流量，MTail为窑尾烟室风流量，MClear为入炉一次风流量，MRaw入炉物料流量，MCoal入炉煤粉流量。MGout和MSout分别为分解炉出口气体流量和入窑物料流量，MError为质量守恒补偿项，基本为0。
[0161] 3)分解炉内能量守恒如下：
[0162] QThird+QTail+QClear+QRaw+QCoal+QR_Coal-QR_Caco3＝QGout+QSout+QError[0163] 其中QThird为三次风显热，QTail为二次风显热，QClear为入炉一次风显热，QRaw为入炉物料显热，QCoal为入炉煤粉显热，QR_Coal为煤燃烧放热，QR_Caco3为碳酸钙分解吸热，QGout为分解炉出口气体显热，QSout为入窑物料显热，QError为补偿项。
[0164] 其中
[0165] QSout＝Cs*Msout*TCal
[0166] QError＝CS*MError*TCal
[0167] 据此求解能量守恒方程，即可获得分解炉模型输出TCal，模型输出值和实际值对比结果如图3所示。可以看出，模型输出值和实际值非常接近，表明分解炉模型准确性高，符合窑炉实际情况。
[0168] (2)回转窑模型
[0169] 在一个或多个实施方案中，回转窑模型结构如下：
[0170] IKiln＝f(TCal,CoalKiln,TThird)
[0171] 其中回转窑主电机电流Ikiln为模型输出，TCal为分解炉出口温度，TThird为三次风温度，CoalKiln为回转窑喂煤量，Coalcal为分解炉喂煤量。
[0172] 回转窑过程涉及气固、气液和固液反应，具有反应机理复杂，过程噪声大等特点。在一个或多个实施方案中，采用最小二乘建模方法确定回转窑模型参数。回转窑模型如下：
[0173] IKiln＝A1*TCal+A2*CoalKiln+A3*TThird+A4
[0174] 参数采用最小二乘方式进行辨识。在一个或多个实施方案中，A1＝8.0188,A2＝40,A3＝0.8311,A4＝-7882，模型输出值和实际值对比结果如图4所示。可以看出，模型输出值和实际值非常接近，表明回转窑模型准确性高，符合窑炉实际情况。相比于现有模型回转窑机理模型中气固反应反应动力学参数难确定，回转窑内部无检测情况，本发明采用基于输入输出数据的数据驱动建模方法，只需要输入输出数据即可获得精确优化模型。
[0175] (3)篦冷机模型
[0176] TThird＝f(PCool,IKiln,SFan)
[0177] 其中三次风温度TThird为模型输出，PCool为篦冷机篦下压力，IKiln为回转窑主电机电流，SFan为高温风机转速。篦冷机单元不涉及反应，满足质量守恒。
[0178] 篦冷机主体部分为一段篦冷机急冷换热，对于一段篦冷机，气体热量守恒如下：
[0179] QG,In+QC,G＝QG,Out
[0180] 固体能量守恒如下：
[0181] QC,In＝QC,G+QC,Out
[0182] 其中QG,In和QG,Out分别为进入一段篦冷机气体和离开一段篦冷机气体热量；QC,In和QC,Out分别为进入一段篦冷机熟料和离开一段篦冷机一段熟料的热量；其中QC,G为熟料和气体间换热。它们分别通过下述模型表示。
[0183] QG,In＝CG,In*MG,In*TG,In
[0184] 其中MG,In为篦冷机底层风机吹入篦冷机气体质量流量；
[0185] QG,out＝CThird*MThird*TThird+CTail*MTail*TTail
[0186] 其中主要是二次风和三次风热量，二次风、三次风流量与高温风机转速SFan直接关联；
[0187] QC,In＝CC,In*MC*TC,In
[0188] 其中TC,In为进入篦冷机熟料温度，因为现场无实时检测，为表征温度，将进入篦冷机熟料温度与回转窑主电机电流关联，主电机电流越大，进入篦冷机熟料温度越高。
[0189]
[0190] 其中TBase是基准温度值，为1300-1500℃，优选1300-1500℃。IBase是基准电流值，为200-1200A，优选500-900A。在一个或多个实施方案中，TBase＝1400℃,IBase＝790A。
[0191] QC,Out＝CC,Out*MC*TC,Out
[0192] TC,Out为离开篦冷机一段熟料温度，温度范围大约为600～700℃。
[0193] 熟料与气体之间换热QC,G是篦冷机模型重点，QC,G表征如下所示：
[0194] QC,G＝K*a*(Tc-Tg)*VCooler
[0195] 其中K为气体间的综合换热系数，a为单位体积床层的有效表面积，(Tc-Tg)为气固之间平均温差，VCooler为一段篦冷机熟料的体积。
[0196]
[0197]
[0198] (Tc-Tg)＝(TC,In+TC,Out-TG,In-TG,out)
[0199] VCooler＝LCooler*WCooler*HCooler
[0200] 其中气体对流换热系数h为45-60w/(m2·K)，优选50-55w/(m2·K)；熟料导热系数λc为1-3w/(m·K)，优选1.5-2.5w/(m·K)；颗粒形状校正系数Φ为0.2-0.4，优选0.25-0.3，颗粒的透热深度x为0.003-0.01，优选0.005-0.0055；空隙率ε为0.4-0.6，优选0.4-0.55，物料尺寸dp为0.008-0.02m，优选0.01-0.02m。TC,In,TC,Out,TG,In,TG,out分别为进入一段篦冷机物料温度，离开一段篦冷机物料温度，进入一段篦冷机气体温度，离开一段篦冷机气体温2
度。在一个或多个实施方案中，h＝54w/(m·K)，λc＝2w/(m·K)，Φ＝1/4，x＝0.005m，ε＝
0.5，dp＝0.01m。在一个或多个实施方案中，一段篦冷机长度LCooler＝50m、宽度,WCooler＝9m。
由于高度HCooler无法实时监测，采用篦冷机篦下压力PCool表征料层高度，关系为非线性关系：
[0201] HCooler＝B1*PCool2+B2*PCool+B3
[0202] 本领域周知，可以根据历史记录数据获取篦下压力与料层厚度的关系。在一个或多个实施方案中，模型参数B1＝-0.1006,B2＝1.1051,B3＝-2.2134。
[0203] 据此，求解能量守恒，即可获得三次风温度TThird，模型输出值和实际值对比结果如图4所示。可以看出，模型输出值和实际值非常接近，表明篦冷机模型准确性高，符合窑炉实际情况。
[0204] 步骤S2：构建优化模型
[0205] 本步骤根据目标、变量和条件构建优化模型。
[0206] (1)优化目标：单产熟料煤耗最小
[0207] min Energy＝(CoalCal+CoalKiln)/Clinker
[0208] 其中Clinker为熟料产量，分解炉用煤量CoalCal、回转窑用煤量CoalKiln如步骤S1所示。
[0209] (2)决策变量：对应于底层回路设定值，三个决策变量为分解炉出口温度TCal，回转窑主电机电流IKiln，篦冷机篦下压力PCool。
[0210] (3)约束条件：包括但不限于关键质量指标游离氧化钙含量Fcao约束，分解炉预热器出口氧气含量O2约束，入窑生料分解率Rate约束，以及决策变量TCal、IKiln及PCool上下限约束。在一个或多个实施方案中，上述参数的模型结构和约束条件如下所示。
[0211] 其中Fcao模型结构如下所示：
[0212] 0.6≤Fcao＝f(TCal,IKiln,PCool)≤1.2
[0213] 游离氧化钙模型采用辅助变量为分解炉出口温度TCal，回转窑主电机电流IKiln和篦冷机篦下压力PCool。对游离氧化钙模型采用输入输出数据进行参数辨识得到如下模型：
[0214] Fcao＝6.0161-0.002*TCal-0.003*IKiln-0.2*PCool
[0215] 其中O2模型结构如下所示：
[0216] O2＝f(CoalCal,SFan)≥0.6
[0217] 预热器出口O2含量如下：
[0218] O2＝O2Third+O2Tail-O2CalCoal
[0219] 其中三次风氧气O2Third主要成分是空气，氧气成分约为21％，窑尾烟室氧气O2Tail成分与回转窑喂煤量CoalKiln有关，分解炉氧气O2CalCoal成分与分解炉喂煤量CoalCal有关。
[0220]
[0221]
[0222] 其中kCarbon为煤粉含碳量，O2摩尔质量MO2＝32g/mol，焦炭摩尔质量MC＝12g/mol。
[0223] 其中Rate模型结构如下所示：
[0224] Rate＝f(TCal,FRaw)≥0.9
[0225] 分解率模型如下：
[0226] Rate＝RRaw/FRaw
[0227] 其中RRaw为生料分解速率，FRaw为单位时间进入分解炉生料。在一个实施方式中，FRaw＝220kg/s。
[0228] 在一个或多个实施方案中，决策变量约束上下限为：分解炉出口温度：860℃≤TCal≤880℃，回转窑主电机电流：650A≤IKiln≤850A，篦冷机篦下压力：4.2Kpa≤PCool≤5Kpa。优选地，决策变量约束上下限为：分解炉出口温度：862℃≤TCal≤875℃，回转窑主电机电流：700A≤IKiln≤820A，篦冷机篦下压力：4.4Kpa≤PCool≤4.8Kpa。在一个或多个实施方案中，决策变量约束上下限为：分解炉出口温度：865℃≤TCal≤875℃，回转窑主电机电流：740A≤IKiln≤800A，篦冷机篦下压力：4.45Kpa≤PCool≤4.55Kpa
[0229] 步骤S3：求解优化模型
[0230] 本步骤中，对步骤S2中的优化模型求解，可以获得底层回路优化设定值分解炉出口温度TCal，回转窑主电机电流IKiln，篦冷机篦下压力PCool。通过修改设定值达到节能降耗的目的。
[0231] 在一个或多个实施方案中，采用模式搜索算法求解优化模型。模式搜索算法求解过程如下。
[0232] 从初始点开始，包括两种类型搜索：探测移动和模式移动，其中探测移动沿着n个坐标轴方向进行，用以确定新的基点和有利于函数值下降的方向，模式移动沿相邻两个基点连线方向进行，目的使得函数值下降更快。算法步骤如下：
[0233] (1)给定初始点x1∈Rn，初始步长δ，
加速因子α≥1，缩减率0＜β＜1，
精度ε＞0。令y1＝x1，k＝1，j＝1。
[0234] (2)轴向搜索：
[0235] 如果f(yj+δej)＜f(yj)，令yj+1＝yj+δej，转步骤(3)；
[0236] 如果f(yj-δej)＜f(yj)，令yj+1＝yj-δej，转步骤(3)；
[0237] 否则，令yj+1＝yj。
[0238] (3)如果j＜n，则令j：＝j+1，转步骤(2)。如果f(yn+1)＜f(xk)，转步骤(4)。否则，转步骤(5)。
[0239] (4)模式搜索：令xk+1＝yn+1，y1＝xk+1+α(xk+1-xk)。令k：＝k+1，j＝1，转步骤(2)。
[0240] (5)如果δ≤ε，则停止，得到点xk；
[0241] 否则，令δ:＝βδ,y1＝xk，xk+1＝xk。令k：＝k+1，j＝1，转步骤(2)。
[0242] 在一个或多个实施方案中，模式搜索算法初始步长δ设置为1，加速因子α设置为1.2，缩减率β设置为0.6。精度设置为0.01。对步骤S2中优化模型进行求解。
[0243] 求解结果表1所示。
[0244] 表1优化结果对比
[0245]优化变量 优化前 优化后
分解炉出口温度 871℃ 870.4℃
回转窑主电机电流 782A 742A
篦冷机篦下压力 4.51Kpa 4.55Kpa
优化目标 优化前 优化后
单产熟料煤耗 107.2kg/t 103.9kg/t
[0246] 从表1中可以看出，通过全流程过程优化，单吨熟料煤耗可以降低约3.3kg/t。与现有的基于专家知识的优化系统相比，本发明能够精确表征过程运行状态，在满足优化目标和约束条件的前提下，得到操作决策变量，指导系统优化。与传统设计优化相比，在不增加企业
硬件设施的前提下，通过精确的过程建模与优化求解，有效起到节能降耗的目的，提高企业利润。
[0247] 此外，本发明还公开了煅烧过程煤耗优化系统，包括计算机以及运行于计算机上的计算机程序，该计算机程序在该计算机上运行如前述实施例的煅烧过程煤耗优化方法。方法的具体步骤不再赘述。
[0248] 此外，本发明还揭示了存储计算机程序的计算机可读存储介质，其特征在于，存储介质上所存储的计算机程序运行后执行如前述实施例的煅烧过程煤耗方法。方法的具体步骤不再赘述。
[0249] 尽管为使解释简单化将上述方法图示并描述为一系列动作，但是应理解并领会，这些方法不受动作的次序所限，因为根据一个或多个实施例，一些动作可按不同次序发生和/或与来自本文中图示和描述或本文中未图示和描述但本领域技术人员可以理解的其他动作并发地发生。
[0250] 本领域技术人员将进一步领会，结合本文中所公开的实施例来描述的各种解说性逻辑板块、模块、
电路、和算法步骤可实现为
电子硬件、计算机
软件、或这两者的组合。为清楚地解说硬件与软件的这一可互换性，各种解说性组件、框、模块、电路、和步骤在上面是以其功能性的形式作一般化描述的。此类功能性是被实现为硬件还是软件取决于具体应用和施加于整体系统的设计约束。技术人员对于每种特定应用可用不同的方式来实现所描述的功能性，但这样的实现决策不应被解读成导致脱离了本发明的范围。
[0251] 结合本文所公开的实施例描述的各种解说性逻辑板块、模块、和电路可用通用处理器、数字
信号处理器(DSP)、
专用集成电路(ASIC)、现场可编程
门阵列(FPGA)或其它
可编程逻辑器件、分立的门或晶体管逻辑、分立的硬件组件、或其设计成执行本文所描述功能的任何组合来实现或执行。通用处理器可以是
微处理器，但在替换方案中，该处理器可以是任何常规的处理器、
控制器、
微控制器、或状态机。处理器还可以被实现为计算设备的组合，例如DSP与微处理器的组合、多个微处理器、与DSP核心协作的一个或多个微处理器、或任何其他此类配置。
[0252] 结合本文中公开的实施例描述的方法或算法的步骤可直接在硬件中、在由处理器执行的
软件模块中、或在这两者的组合中体现。软件模块可驻留在RAM
存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、
硬盘、可移动盘、CD-ROM、或本领域中所知的任何其他形式的存储介质中。示例性存储介质耦合到处理器以使得该处理器能从/向该存储介质读取和写入信息。在替换方案中，存储介质可以被整合到处理器。处理器和存储介质可驻留在ASIC中。ASIC可驻留在用户终端中。在替换方案中，处理器和存储介质可作为分立组件驻留在用户终端中。
[0253] 在一个或多个示例性实施例中，所描述的功能可在硬件、软件、
固件或其任何组合中实现。如果在软件中实现为计算机程序产品，则各功能可以作为一条或更多条指令或代码存储在计算机可读介质上或藉其进行传送。计算机可读介质包括计算机存储介质和通信介质两者，其包括促成计算机程序从一地向另一地转移的任何介质。存储介质可以是能被计算机
访问的任何可用介质。作为示例而非限定，这样的计算机可读介质可包括RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁存储设备、或能被用来携带或存储指令或数据结构形式的合意程序代码且能被计算机访问的任何其它介质。任何连接也被正当地称为计算机可读介质。例如，如果软件是使用同轴
电缆、光纤电缆、双绞线、数字订户线(DSL)、或诸如红外、无线电、以及
微波之类的无线技术从web
网站、
服务器、或其它远程源传送而来，则该同轴电缆、光纤电缆、双绞线、DSL、或诸如红外、无线电、以及微波之类的无线技术就被包括在介质的定义之中。如本文中所使用的盘(disk)和碟(disc)包括压缩碟(CD)、激光碟、光碟、数字多用碟(DVD)、
软盘和蓝光碟，其中盘(disk)往往以磁的方式再现数据，而碟(disc)用激光以光学方式再现数据。上述的组合也应被包括在计算机可读介质的范围内。
[0254] 提供对本公开的先前描述是为使得本领域任何技术人员皆能够制作或使用本公开。对本公开的各种修改对本领域技术人员来说都将是显而易见的，且本文中所定义的普适原理可被应用到其他变体而不会脱离本公开的精神或范围。由此，本公开并非旨在被限定于本文中所描述的示例和设计，而是应被授予与本文中所公开的原理和新颖性特征相一致的最广范围。
