CN116151022A - 一种基于热平衡计算的水泥回转窑温度实时估计方法 - Google Patents
一种基于热平衡计算的水泥回转窑温度实时估计方法 Download PDFInfo
- Publication number
- CN116151022A CN116151022A CN202310211711.2A CN202310211711A CN116151022A CN 116151022 A CN116151022 A CN 116151022A CN 202310211711 A CN202310211711 A CN 202310211711A CN 116151022 A CN116151022 A CN 116151022A
- Authority
- CN
- China
- Prior art keywords
- kiln
- heat
- temperature
- rotary kiln
- cement
- 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
- 239000004568 cement Substances 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000004364 calculation method Methods 0.000 title claims abstract description 37
- 239000000463 material Substances 0.000 claims abstract description 129
- 238000006243 chemical reaction Methods 0.000 claims abstract description 44
- 230000008569 process Effects 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 5
- 239000003245 coal Substances 0.000 claims description 42
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 32
- 239000003546 flue gas Substances 0.000 claims description 32
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 claims description 22
- 239000000779 smoke Substances 0.000 claims description 21
- 230000005855 radiation Effects 0.000 claims description 18
- 229910000019 calcium carbonate Inorganic materials 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 11
- 230000007246 mechanism Effects 0.000 claims description 10
- 239000002994 raw material Substances 0.000 claims description 7
- 238000002485 combustion reaction Methods 0.000 claims description 6
- 239000007789 gas Substances 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 235000010216 calcium carbonate Nutrition 0.000 claims description 5
- 230000004927 fusion Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 229910018072 Al 2 O 3 Inorganic materials 0.000 claims description 4
- 229910004298 SiO 2 Inorganic materials 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 claims description 3
- AGWMJKGGLUJAPB-UHFFFAOYSA-N aluminum;dicalcium;iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Ca+2].[Ca+2].[Fe+3] AGWMJKGGLUJAPB-UHFFFAOYSA-N 0.000 claims description 3
- 229910052918 calcium silicate Inorganic materials 0.000 claims description 3
- 235000012241 calcium silicate Nutrition 0.000 claims description 3
- JHLNERQLKQQLRZ-UHFFFAOYSA-N calcium silicate Chemical compound [Ca+2].[Ca+2].[O-][Si]([O-])([O-])[O-] JHLNERQLKQQLRZ-UHFFFAOYSA-N 0.000 claims description 3
- HOOWDPSAHIOHCC-UHFFFAOYSA-N dialuminum tricalcium oxygen(2-) Chemical compound [O--].[O--].[O--].[O--].[O--].[O--].[Al+3].[Al+3].[Ca++].[Ca++].[Ca++] HOOWDPSAHIOHCC-UHFFFAOYSA-N 0.000 claims description 3
- BCAARMUWIRURQS-UHFFFAOYSA-N dicalcium;oxocalcium;silicate Chemical compound [Ca+2].[Ca+2].[Ca]=O.[O-][Si]([O-])([O-])[O-] BCAARMUWIRURQS-UHFFFAOYSA-N 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 239000007788 liquid Substances 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 235000019976 tricalcium silicate Nutrition 0.000 claims description 3
- 229910021534 tricalcium silicate Inorganic materials 0.000 claims description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 2
- 238000002844 melting Methods 0.000 claims description 2
- 230000008018 melting Effects 0.000 claims description 2
- 239000011707 mineral Substances 0.000 claims description 2
- 235000010755 mineral Nutrition 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 9
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000014509 gene expression Effects 0.000 description 8
- 239000000292 calcium oxide Substances 0.000 description 7
- ODINCKMPIJJUCX-UHFFFAOYSA-N calcium oxide Inorganic materials [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 description 7
- 238000010801 machine learning Methods 0.000 description 4
- 238000001354 calcination Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 3
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 2
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 235000012054 meals Nutrition 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- BRPQOXSCLDDYGP-UHFFFAOYSA-N calcium oxide Chemical compound [O-2].[Ca+2] BRPQOXSCLDDYGP-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000376 reactant Substances 0.000 description 1
- 235000012239 silicon dioxide Nutrition 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Curing Cements, Concrete, And Artificial Stone (AREA)
Abstract
本发明公开了一种基于热平衡计算的水泥回转窑温度实时估计方法,包括如下6个步骤:步骤1:获取计算所需参数数据;步骤2:构建零维模型;步骤3:设置物料温度初始值,初步求解物料化学反应热;步骤4:通过水泥回转窑内部的物料平衡、热平衡初步求解回转窑内温度;步骤5:利用水泥回转窑输入输出的热平衡对求解效果进行评估;步骤6:使用步骤4得到的温度参数求解物料化学反应热,然后重新求解回转窑温度,迭代上述过程,输出回转窑温度的计算结果。本发明实现了水泥回转窑运行数据实时输入、温度参数实时输出的分钟级、秒级实时温度估计,显著提高了水泥回转窑温度模型的实用性,有助于生产企业对回转窑温度进行实时的监控和调整。
Description
技术领域
本发明属于水泥生产过程建模、监控和仿真领域,特别涉及一种基于热平衡计算的水泥回转窑温度实时估计方法。该方法主要适用于分解炉、回转窑、篦冷机等设备参与的干式水泥生产工艺中水泥回转窑温度的实时估计。
背景技术
在水泥生产过程中,水泥回转窑温度的监测是行业普遍重视的一大问题,其原因有三:其一,企业常常采用水泥熟料的游离钙(以下简称“f-CaO”)含量作为水泥产品质量的评价指标,f-CaO表征了水泥生料煅烧过程中氧化钙(以下简称“CaO”)与二氧化硅(以下简称“SiO2”)、氧化铝(以下简称“Al2O3”)和氧化铁(以下简称“Fe2O3”)等成分结合的完全程度,对水泥的强度和安定性影响较大,而水泥熟料的f-CaO含量与水泥回转窑的煅烧温度高度相关;其二,烧成系统乃至整个水泥生产过程的煤耗与回转窑的煅烧温度直接挂钩,水泥回转窑的温度监测可以为水泥生产的能耗情况提供参考;其三,水泥回转窑温度的异常可能会引发红窑、掉窑砖等问题,对设备、人员造成损坏和伤害,因此水泥回转窑的温度监测对于水泥回转窑的安全运行起到关键作用。
回转窑温度的影响因素众多,包括生料喂料量、生料成分、回转窑喂煤量、二次风温度、二次风量等等,甚至回转窑窑速的变化都会对回转窑温度造成一定的影响。由于回转窑筒体为转动部件,无法安装热电偶等接触式测温器件,且回转窑内部高温、高灰的环境也不利于红外测温仪等非接触式测温装置的部署,因此生产过程中大多采用二次风温度、窑尾烟气温度等容易测量的运行参数间接反映回转窑内的温度,或者通过机理建模、机器学习建模等方法,通过已知的运行参数对回转窑的温度进行推理和建模。
其中,回转窑温度的机器学习建模方法主要依靠传统机器学习算法或神经网络算法从标签数据中学习、总结、归纳运行参数和回转窑温度之间的关系,然而机器学习模型存在的局部最优值问题以及神经网络训练的不稳定问题会对模型的可靠性产生很大的影响。
用于回转窑温度建模的机理模型根据维度划分主要分为多维模型和一维模型两类,多维模型依赖Fluent、COMSOL等CFD软件对回转窑内部空间进行网格划分,并在所设定的复杂边界条件下开展大量的迭代计算,虽然其建模结果可以良好地反映窑内温度和气体流动的分布状况,但是其建模操作复杂,计算时间长,且存在计算结果不收敛的风险,从生产应用的角度来看,其实用性和可靠性并不强。
一维模型忽略了回转窑内物料和能量在除轴向以外其他方向上的运动和传递,而只考虑轴向上物质和能量的变化情况,相比多维模型在计算速度和灵活性方面有了较为明显的提升,然而在实际应用中,一维模型也暴露出了两个不可忽视的问题:
(1)现有一维模型对于窑内煤粉的运动方式、火焰的形态以及在各个长度处生料的热状态等因素都进行了理想化假设,然而在实际运行条件下,进入回转窑的煤粉风呈现出紊流、旋流的形态,煤粉的运动方式、燃烧方式和一维模型的描述之间存在明显的出入,且由于热惯性的原因,在回转窑内某一特定长度处,生料的实际状态与模型推理结果之间也会产生较大的差异,而不稳定的窑况会进一步放大上述差异,致使一维模型对于回转窑温度分布的描述精度欠佳,甚至出现在模型验证过程中回转窑输入、输出能量明显不守恒的现象,最终导致模型无法正确描述窑内的温度状况;(2)尽管相比多维模型,一维模型复杂度有所下降,但在实际应用的过程中,构建一维温度场的过程中依然存在较大的计算量,并不能满足工业生产中回转窑温度分钟级、秒级实时监测的需求。
发明内容
为了克服现有回转窑温度建模技术中存在的不足,特别是模型复杂度、实时性方面的问题,本发明的目的在于提供一种基于热平衡计算的水泥回转窑温度实时估计方法,在实时估计水泥回转窑温度的同时,也取得较高的估计精度。方法具体包括了如下6个步骤:
步骤1:获取计算所需的水泥回转窑静态尺寸参数、动态运行参数数据;
步骤2:综合水泥回转窑内部结构状态、传热机理、物化反应,构建零维模型;
步骤3:设置物料温度初始值,初步求解物料化学反应热;
步骤4:通过水泥回转窑内部的物料平衡、热平衡初步求解回转窑内温度;
步骤5:利用水泥回转窑输入输出的热平衡对求解效果进行评估并调整求解方法;
步骤6:使用步骤4得到的回转窑内温度求解物料化学反应热,然后使用求解得到的物料化学反应热,综合回转窑内部的物料平衡、热平衡重新求解回转窑温度,迭代上述过程直到相邻两次温度计算结果的差异为可接受的范围,输出回转窑温度的计算结果。
进一步地,所述的步骤1中的水泥回转窑静态尺寸参数主要包括窑长、窑外径、窑内径、窑倾角、物料休止角等回转窑的结构尺寸参数,动态运行参数主要包括生料CaO、SiO2、Al2O3、Fe2O3含量化验值、燃煤低位发热量化验值、分解炉出口温度、二次风温度、窑尾烟室温度、喂料量、喂煤量、回转窑进风量。
进一步地,所述的步骤2中水泥回转窑内部结构状态的表示方法如下:
其中,h为物料料床高度,Ag为窑内的气体横截面积,R为窑内半径,β为窑倾角,γ为窑内物料的休止角,φv为窑内物料的体积流量,n为窑转速,θ为料床角,其计算公式为:
θ=2arcos((R-h)/R)。
进一步地,所述的步骤2中水泥回转窑传热机理的表示方法如下:
其中,分别表示烟气对物料的热对流、烟气对内壁的热对流、烟气对物料的热辐射、烟气对内壁的热辐射、内壁对物料的热传导、内壁对物料的热辐射、内壁对窑壳的热传导、窑壳对外界的热对流以及窑壳对外界的热辐射,Tg、Ts、Tw、Tsh、Text分别表示烟气、物料、内壁、窑壳和外界的热力学温度,hgs、hgw、hnv分别表示烟气与物料、烟气与内壁、窑壳与外界之间的换热系数,kws、kcomp分别表示内壁与物料、内壁与窑壳之间的导热系数,∈g、∈s、∈w、∈sh分别表示烟气、物料、内壁、窑壳的发射率,αg表示烟气的吸收率,σ为斯蒂芬-玻尔兹曼常数,Pgs、Pgw、Pcws、Prws、Prsw、Psh分别表示烟气与物料、烟气与内壁、内壁与物料传导、内壁到物料辐射、物料到内壁辐射、窑壳与外界之间的换热周长;
进一步地,所述的步骤2中水泥回转窑物化反应的表示方法如下:
其中,Emelt表示物料在熔融过程中吸收的热量,其中ms为喂料量,emelt为物料熔融潜热,其值为416kJ/kg,TS、TL分别为物料完全为固态、完全转换为液态时的温度,分别为1273K和2200K:Echem表示物料中的矿物成分在高温下发生化学反应所产生的热量,即物料化学反应热,其中,和分别为物料中碳酸钙(CaCO3)、硅酸二钙(C2S)、硅酸三钙(C3S)、铝酸三钙(C3A)、铁铝酸四钙(C4AF)的单位摩尔反应焓值,rCacO3、rC2s、rC3s、rC3A、rC4AF分别为物料中CaCO3、C2S、C3S、C3A和C4AF的单位千克反应摩尔数。
进一步地,所述的步骤4中水泥回转窑内部物料平衡、热平衡的表示方法如下:
ma=Agρgug-ηcoalmcoal;
其中,ma为水泥回转窑的进风量,Cpa、Cps分别为水泥回转窑进风和物料的比热容,Ta、Ts、Tprec分别为二次风、物料、分解炉出口的热力学温度,mcoal为水泥回转窑的喂煤量,ηcoal为煤粉燃烧率,qnet为燃煤的低位发热量,ρg为烟气的密度,ug为烟气的流速。
进一步地,所述的步骤5中水泥回转窑输入输出热平衡的表示方法如下:
maCpaTa+ηcoalmcoalqnet+CpsTprecms=CpgTtail(ma+ηcoalmcoal)+CpsTsms;
其中,Cpg为烟气比热容,Ttail为窑尾烟室温度,等式左边为水泥回转窑的热量输入,主要包括二次风携带的热量、煤粉燃烧放热以及分解炉输入水泥回转窑的物料携带的热量;等式右边为水泥回转窑的热量输出,主要包括窑尾烟气中携带的热量以及输出的熟料中携带的热量。由于在温度建模的过程中窑尾烟室温度未知,故等式左右两边计算得到的热量一般不会完全相等。定义热平衡误差如下:
模型计算所得的输入输出热量之间的误差可以作为模型准确度的一个评价指标。一般而言,5%以下的热平衡误差是可以接受的,如热平衡误差过高,则需调整模型参数,或者重新寻找其他的求解起点对水泥回转窑的温度进行求解。
本发明的有益效果:
本发明放弃了水泥回转窑多维、一维的传统机理建模方法,将水泥回转窑看作一个整体构建零维模型,然后使用热平衡求解水泥回转窑温度,并采用回转窑输入输出的热平衡误差评估单次求解的效果,最后通过对水泥回转窑温度和化学反应热的迭代计算确保求解结果的准确度,规避了水泥回转窑多维建模存在的操作复杂、计算时间长的问题,同时弥补了回转窑一维建模精度欠佳、实时性差的缺陷,实现了水泥回转窑运行数据实时输入、温度参数实时输出的分钟级、秒级实时温度估计,显著提高了水泥回转窑温度模型的实用性,有助于生产企业对回转窑温度进行实时的监控和调整。
附图说明
图1是本发明所展示的一种基于热平衡计算的水泥回转窑温度实时估计方法的流程图。
图2是本发明所展示的一种基于热平衡计算的水泥回转窑温度实时估计方法的水泥回转窑热平衡示意图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。需要注意的是,本发明所涉及到的表达式中的各个符号仅为简化和表达清楚服务,可以在不同的例子中具有不同的表述方法,因而本发明的关键在于水泥回转窑温度求解的出发点和求解思路;除此之外,本发明省略了对公知方法,如气体成分密度、比热容的计算方法,回转窑内部传热尺寸的计算方法,换热系数的计算方法以及物料化学反应热的计算方法的描述,以避免本发明受到不必要的限制。
本发明的目的在于提出一种简单、有效的基于热平衡计算的水泥回转窑温度实时估计方法,总体来说,本发明的实施方法流程图如图1所示,具体包括了如下6个步骤:
步骤1:获取计算所需的水泥回转窑静态尺寸参数、动态运行参数数据;
步骤2:综合水泥回转窑内部结构状态、传热机理、物化反应,构建零维模型;
步骤3:设置物料温度初始值,初步求解物料化学反应热;
步骤4:通过水泥回转窑内部的物料平衡、热平衡初步求解回转窑内温度;
步骤5:利用水泥回转窑输入输出的热平衡对求解效果进行评估并调整求解方法;
步骤6:使用步骤4得到的回转窑内温度求解物料化学反应热,然后使用求解得到的物料化学反应热,综合回转窑内部的物料平衡、热平衡重新求解回转窑温度,迭代上述过程直到相邻两次温度计算结果的差异为可接受的范围,输出回转窑温度的计算结果。
步骤1为实施本发明所述方法的预备工作,其目的在于为步骤2中所构建的零维模型提供必要的输入。其中,水泥回转窑静态尺寸参数主要包括窑长、窑外径、窑内径、窑倾角、物料休止角等回转窑的结构尺寸参数,动态运行参数主要包括生料CaO、SiO2、Al2O3、Fe2O3含量化验值、燃煤低位发热量化验值、分解炉出口温度、二次风温度、窑尾烟室温度、喂料量、喂煤量、回转窑进风量等可以反映回转窑运行状态的、实时更新的参数,动态运行参数的数据一般通过工业DCS的数据库获取。
步骤2中,首先通过获取到的数据计算出物料料床高度:
其中,h为物料料床高度,R为窑内半径,β为窑倾角,γ为窑内物料的休止角,φv为窑内物料的体积流量,n为窑转速。计算获得的料床高度可以用于料床角θ的计算:
θ=2arcos((R-h)/R);
计算所得的料床角可以进一步算得水泥回转窑内气体的横截面积Ag:
回转窑内部结构状态的计算主要是为后续模型的构建服务。更具体地,通过上述计算式得到的h和θ主要用于后续传热机理建模过程中传热尺寸的计算,而Ag则主要用于水泥回转窑内气体质量流量的计算。
回转窑内部结构状态计算完毕后,对水泥回转窑内的传热机理进行表达:
其中,分别表示烟气对物料的热对流、烟气对内壁的热对流、烟气对物料的热辐射、烟气对内壁的热辐射、内壁对物料的热传导、内壁对物料的热辐射、内壁对窑壳的热传导、窑壳对外界的热对流以及窑壳对外界的热辐射,Tg、Ts、Tw、Tsh、Text分别表示烟气、物料、内壁、窑壳和外界的热力学温度,hgs、hgw、hnv分别表示烟气与物料、烟气与内壁、窑壳与外界之间的换热系数,kws、kcomp分别表示内壁与物料、内壁与窑壳之间的导热系数,∈g、∈s、∈w、∈sh分别表示烟气、物料、内壁、窑壳的发射率,αg表示烟气的吸收率,σ为斯蒂芬-玻尔兹曼常数,Pgs、Pgw、Pcws、Prws、Prsw、Psh分别表示烟气与物料、烟气与内壁、内壁与物料(传导)、内壁到物料(辐射)、物料到内壁(辐射)、窑壳与外界之间的传热周长,这些周长的计算本质上是涉及h、θ、R和水泥回转窑窑长的简单几何计算,故其计算方法不再一一赘述。向上述所有传热表达式代入步骤1中采集到的参数数据和步骤2中计算得到的内部结构状态数据,最终将所有表达式化简为只与温度有关的计算式。
考虑到物料在水泥回转窑内煅烧的过程中会伴随着物化反应的发生,因此还需对水泥回转窑内物料的物化反应热进行表达,更具体地,是对物料的熔融热和化学反应热进行表达。
其中,物料熔融热的表达式为:
其中,ms为喂料量,emelt为物料熔融潜热,其值为416kJ/kg,TS、TL分别为物料完全为固态、完全转换为液态时的温度,分别为1273K和2200K。由表达式可知,在喂料量ms己知的条件下,Emelt只与物料温度Ts有关。物料化学反应热的表达式如下:
其中,和分别为物料中碳酸钙(CaCO3)、硅酸二钙(C2S)、硅酸三钙(C3S)、铝酸三钙(C3A)、铁铝酸四钙(C4AF)的单位摩尔反应焓值,rCaCO3、rC2s、rC3S、rC3A、rC4AF分别为物料中CaCO3、C2S、C3S、C3A和C4AF的单位千克反应摩尔数。该方法通过联立各个反应物的阿伦尼乌斯方程,结合化学反应过程中的元素质量守恒联立方程组,最终求得rcacO3、rC2S、rc3s、rC3A和rC4AF的数值。在生料化验成分己知的情况下,可以认为物料的化学反应热Echem只与物料温度有关。
步骤3设置了一个初始的物料温度Ts大小Ts_init,并初步计算物料的化学反应热Echem_init,其目的在于解决求解过程中物料温度、化学反应热之间相互作用的问题,即物料温度决定化学反应热,化学反应热的变更又会影响物料温度。因此,将化学反应热和回转窑温度一并求解的结果往往并不精确甚至不能收敛,需要对化学反应热和回转窑温度分开求解,并将化学反应热求解结果作为常值用于步骤4中回转窑温度的求解过程,并对上述过程在步骤6中反复迭代计算,方可使物料温度和该温度条件下的化学反应热完全对应。
步骤4基于步骤2中的水泥回转窑建模结果,结合如图2所示的水泥回转窑内部物料平衡、热平衡规律建立如下的五元多阶方程组:
ma=Agρgug-ηcoalmcoal;
其中,ma为水泥回转窑的进风量,Cpa、Cps分别为水泥回转窑进风和物料的比热容,Ta、Ts、Tprec分别为二次风、物料、分解炉出口的热力学温度,mcoal为水泥回转窑的喂煤量,ηcoal为煤粉燃烧率,qnet为燃煤的低位发热量,ρg为烟气的密度,ug为烟气的流速。
在上述方程组中,未知量为ug、Tg、Ts、Tw和Tsh,对应方程组中的五个等式。对上述未知量选取合适的求解起点并通过方程组求解,可以初步获得一组水泥回转窑的温度参数。
步骤5针对步骤4的计算结果,通过水泥回转窑的输入输出热平衡误差评估温度建模的准确度。如图2所示,水泥回转窑输入输出热平衡的表示方法如下:
maCpaTa+ηcoalmcoalqnet+CpsTprecms=CpgTtail(ma+ηcoalmcoal)+CpsTsms;
其中,Cpg为烟气比热容,Ttail为窑尾烟室温度,等式左边为水泥回转窑的热量输入,主要包括二次风携带的热量、煤粉燃烧放热以及分解炉输入水泥回转窑的物料携带的热量;等式右边为水泥回转窑的热量输出,主要包括窑尾烟气中携带的热量以及输出的熟料中携带的热量。接着,定义热平衡误差如下:
由于窑尾烟室温度不参与水泥回转窑温度的建模,因此这一热平衡误差可以较好地反映模型对于水泥回转窑各部分温度估计的精确与否。一般而言,5%以下的热平衡误差是可以接受的,如热平衡误差过高,则需调整模型参数,或者重新寻找其他的求解起点对水泥回转窑的温度进行求解。
步骤6为物料化学反应热、水泥回转窑温度反复迭代计算的过程,其流程如下:
(1)通过步骤4中解得的物料温度Ts重新求解物料化学反应热Echem;
(2)将(1)中解得的物料化学反应热Echem通过步骤4中的水泥回转窑内部物料平衡、热平衡方程组重新计算水泥回转窑温度;
(3)将(2)中解得的水泥回转窑温度与上一次迭代过程中的计算结果进行对比,若温度差异为可接受的范围,则输出温度求解结果,否则,重复(1)、(2)、(3)三个环节,若求解Ts和Echem的方程组本身收敛,则每经过一次迭代计算,相邻两次求解得到的水泥回转窑温度差异会不断减小,直至满足输出要求。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明的保护范围应以所附权利要求为准。
Claims (6)
1.一种基于热平衡计算的水泥回转窑温度实时估计方法,其特征在于,包括如下6个步骤:
步骤1:获取计算所需的水泥回转窑静态尺寸参数、动态运行参数数据;
步骤2:综合水泥回转窑内部结构状态、传热机理、物化反应,构建零维模型;
步骤3:设置物料温度初始值,初步求解物料化学反应热;
步骤4:通过水泥回转窑内部的物料平衡、热平衡初步求解回转窑内温度;
步骤5:利用水泥回转窑输入输出的热平衡对求解效果进行评估并调整求解方法;5%以下的热平衡误差予以接受,热平衡误差过高则调整模型参数,或者重新寻找其他的求解起点对水泥回转窑的温度进行求解;
步骤6:使用步骤4得到的回转窑内温度求解物料化学反应热,然后使用求解得到的物料化学反应热,综合回转窑内部的物料平衡、热平衡重新求解回转窑温度,迭代上述过程直到相邻两次温度计算结果的差异为可接受的范围,输出回转窑温度的计算结果。
2.如权利要求1所述的方法,其特征在于,步骤1中,所述的水泥回转窑静态尺寸参数包括窑长、窑外径、窑内径、窑倾角、物料休止角结构尺寸参数,所述的动态运行参数为反映回转窑运行状态的、实时更新的参数,包括生料CaO、SiO2、Al2O3、Fe2O3含量化验值、燃煤低位发热量化验值、分解炉出口温度、二次风温度、窑尾烟室温度、喂料量、喂煤量、回转窑进风量。
3.如权利要求1所述的方法,其特征在于,步骤2、3中:
所述的水泥回转窑内部结构状态的表示方法如下:
其中,h为物料料床高度,Ag为窑内的气体横截面积,R为窑内半径,β为窑倾角,γ为窑内物料的休止角,φv为窑内物料的体积流量,n为窑转速,θ为料床角,其计算公式为:
θ=2arcos((R-h)/R);
所述的水泥回转窑传热机理的表示方法如下:
其中,分别表示烟气对物料的热对流、烟气对内壁的热对流、烟气对物料的热辐射、烟气对内壁的热辐射、内壁对物料的热传导、内壁对物料的热辐射、内壁对窑壳的热传导、窑壳对外界的热对流以及窑壳对外界的热辐射,Tg、Ts、Tw、Tsh、Text分别表示烟气、物料、内壁、窑壳和外界的热力学温度,hgs、hgw、hnv分别表示烟气与物料、烟气与内壁、窑壳与外界之间的换热系数,kws、kcomp分别表示内壁与物料、内壁与窑壳之间的导热系数,∈g、∈s、∈w、∈sh分别表示烟气、物料、内壁、窑壳的发射率,αg表示烟气的吸收率,σ为斯蒂芬-玻尔兹曼常数,gs、Pgw、Pcws、Prws、Prsw、Psh分别表示烟气与物料、烟气与内壁、内壁与物料传导、内壁到物料辐射、物料到内壁辐射、窑壳与外界之间的换热周长;
所述的水泥回转窑物化反应的表示方法如下:
6.如权利要求1所述的方法,其特征在于,步骤6中,迭代的流程如下:
(1)通过步骤4中解得的物料温度Ts重新求解物料化学反应热Echem;
(2)将(1)中解得的物料化学反应热Echem通过步骤4中的水泥回转窑内部物料平衡、热平衡方程组重新计算水泥回转窑温度;
(3)将(2)中解得的水泥回转窑温度与上一次迭代过程中的计算结果进行对比,若温度差异为可接受的范围,则输出温度求解结果,否则,重复(1)、(2)、(3)三个环节,若求解Ts和Echem的方程组本身收敛,则每经过一次迭代计算,相邻两次求解得到的水泥回转窑温度差异会不断减小,直至满足输出要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310211711.2A CN116151022B (zh) | 2023-03-07 | 2023-03-07 | 一种基于热平衡计算的水泥回转窑温度实时估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310211711.2A CN116151022B (zh) | 2023-03-07 | 2023-03-07 | 一种基于热平衡计算的水泥回转窑温度实时估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116151022A true CN116151022A (zh) | 2023-05-23 |
CN116151022B CN116151022B (zh) | 2024-04-23 |
Family
ID=86338999
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310211711.2A Active CN116151022B (zh) | 2023-03-07 | 2023-03-07 | 一种基于热平衡计算的水泥回转窑温度实时估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116151022B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570244A (zh) * | 2016-10-25 | 2017-04-19 | 浙江邦业科技股份有限公司 | 一种预测水泥回转窑熟料质量的一维仿真方法 |
CN111174569A (zh) * | 2020-01-16 | 2020-05-19 | 武汉科技大学 | 一种在线预测回转窑内煅烧段烟气温度的方法及系统 |
CN112947094A (zh) * | 2021-04-07 | 2021-06-11 | 宁波大学 | 一种水泥回转窑温控pid参数自调整方法 |
CN113139275A (zh) * | 2021-03-22 | 2021-07-20 | 浙江大学 | 一种基于多层矿焦比分布模型的高炉炉喉温度估计方法 |
-
2023
- 2023-03-07 CN CN202310211711.2A patent/CN116151022B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106570244A (zh) * | 2016-10-25 | 2017-04-19 | 浙江邦业科技股份有限公司 | 一种预测水泥回转窑熟料质量的一维仿真方法 |
WO2018076403A1 (zh) * | 2016-10-25 | 2018-05-03 | 浙江邦业科技股份有限公司 | 一种预测水泥回转窑熟料质量的一维仿真方法 |
CN111174569A (zh) * | 2020-01-16 | 2020-05-19 | 武汉科技大学 | 一种在线预测回转窑内煅烧段烟气温度的方法及系统 |
CN113139275A (zh) * | 2021-03-22 | 2021-07-20 | 浙江大学 | 一种基于多层矿焦比分布模型的高炉炉喉温度估计方法 |
CN112947094A (zh) * | 2021-04-07 | 2021-06-11 | 宁波大学 | 一种水泥回转窑温控pid参数自调整方法 |
Non-Patent Citations (3)
Title |
---|
S. MUJUMDAR等: "T=零维,没有找到可用的对比文件", JOURNALCHEMICAL ENGINEERING, vol. 84, no. 3, 31 March 2006 (2006-03-31), pages 165 - 177 * |
楼波;马晓茜: "高温空气回转窑内燃料着火的突变特性", 华南理工大学学报(自然科学版)., no. 6, 15 June 2008 (2008-06-15), pages 118 - 123 * |
郑坤灿;温治;刘训良;楼国锋;张欣;张衍国;武文斐;: "高温散料气-固换热过程通用数学模型的研究", 冶金能源, vol. 29, no. 02, 31 March 2010 (2010-03-31), pages 27 - 30 * |
Also Published As
Publication number | Publication date |
---|---|
CN116151022B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018076403A1 (zh) | 一种预测水泥回转窑熟料质量的一维仿真方法 | |
CN111665809B (zh) | 一种水泥回转窑的分段机理建模方法及系统 | |
CN111174569B (zh) | 一种在线预测回转窑内煅烧段烟气温度的方法及系统 | |
Mujumdar et al. | Modeling of rotary cement kilns: applications to reduction in energy consumption | |
Wu et al. | Improving the sustainability of cement clinker calcination process by assessing the heat loss through kiln shell and its influencing factors: A case study in China | |
CN110981240B (zh) | 煅烧过程优化方法和系统 | |
CN104498706B (zh) | 一种链篦机-回转窑-环冷机三位一体的球团生产优化方法 | |
Yin et al. | Design requirements and performance optimization of waste heat recovery systems for rotary kilns | |
CN103939941A (zh) | 一种融入了不可逆热力学的锅炉燃烧优化方法 | |
CN116151022B (zh) | 一种基于热平衡计算的水泥回转窑温度实时估计方法 | |
CN114283895B (zh) | 一种基于热力学机理的氧化锌回转挥发窑高温反应区温度场监测方法 | |
CN109657345B (zh) | 回转式空预器直接漏风流量的实时估计方法 | |
Stadnyk et al. | Energy saving thermal systems on the mobile platform of the mini-bakery | |
Ma et al. | Energy and thermodynamic analysis of a typical cement production system based on experimental and simulated investigations | |
Mian et al. | Energy Integration in the cement industry | |
CN102690072A (zh) | 一种水泥生产过程不完全燃烧碳化物在线监测控制方法 | |
CN113868952A (zh) | 一种回转窑内温度场及结圈厚度分布在线检测方法 | |
JP2009084620A (ja) | 熱風炉の燃焼制御方法及びその燃焼制御装置 | |
CN112270145B (zh) | 一种基于多相流模型的结圈分析方法及系统 | |
CN109135778A (zh) | 一种外热式回转炉干馏装置的反应控制系统及方法 | |
Aldina et al. | Implementation of a mathematical modelling of a rotary cement kilns | |
Fu et al. | Research on soft measurement method of temperature field in cement rotary kiln | |
CN104832914B (zh) | 一种环壁保护燃烧炉及其传热传质数值模拟方法 | |
Mian et al. | 23 European Symposium on Computer Aided Process Engineering: Energy Integration in the cement industry | |
Hu et al. | Process fault detection of cement calcination system based on median filtered DPCA |
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 |