CN104915567B - 基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 - Google Patents
基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 Download PDFInfo
- Publication number
- CN104915567B CN104915567B CN201510342228.3A CN201510342228A CN104915567B CN 104915567 B CN104915567 B CN 104915567B CN 201510342228 A CN201510342228 A CN 201510342228A CN 104915567 B CN104915567 B CN 104915567B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- coal
- clinker
- heat
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000004568 cement Substances 0.000 title claims abstract description 27
- 238000012544 monitoring process Methods 0.000 title claims abstract description 22
- 238000002386 leaching Methods 0.000 claims abstract description 74
- 239000003245 coal Substances 0.000 claims abstract description 56
- 238000002485 combustion reaction Methods 0.000 claims abstract description 43
- 239000002994 raw material Substances 0.000 claims abstract description 43
- 238000003556 assay Methods 0.000 claims abstract description 9
- 238000007689 inspection Methods 0.000 claims abstract description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 22
- 239000003546 flue gas Substances 0.000 claims description 22
- 238000010304 firing Methods 0.000 claims description 18
- 238000000354 decomposition reaction Methods 0.000 claims description 14
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 12
- 238000013461 design Methods 0.000 claims description 12
- 239000001301 oxygen Substances 0.000 claims description 12
- 229910052760 oxygen Inorganic materials 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 12
- 239000007789 gas Substances 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 10
- 239000000446 fuel Substances 0.000 claims description 8
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 6
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 claims description 6
- 239000002802 bituminous coal Substances 0.000 claims description 6
- 229910052593 corundum Inorganic materials 0.000 claims description 6
- 238000007599 discharging Methods 0.000 claims description 6
- 239000000428 dust Substances 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 6
- 239000002912 waste gas Substances 0.000 claims description 6
- 229910001845 yogo sapphire Inorganic materials 0.000 claims description 6
- 238000001035 drying Methods 0.000 claims description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 3
- 229910052681 coesite Inorganic materials 0.000 claims description 3
- 229910052906 cristobalite Inorganic materials 0.000 claims description 3
- 238000001704 evaporation Methods 0.000 claims description 3
- 230000008020 evaporation Effects 0.000 claims description 3
- JEIPFZHSYJVQDO-UHFFFAOYSA-N iron(III) oxide Inorganic materials O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims description 3
- 239000000377 silicon dioxide Substances 0.000 claims description 3
- 229910052682 stishovite Inorganic materials 0.000 claims description 3
- 229910052905 tridymite Inorganic materials 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 239000000126 substance Substances 0.000 abstract description 7
- 230000002411 adverse Effects 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000017525 heat dissipation Effects 0.000 abstract description 2
- 239000002817 coal dust Substances 0.000 description 3
- RHZUVFJBSILHOK-UHFFFAOYSA-N anthracen-1-ylmethanolate Chemical compound C1=CC=C2C=C3C(C[O-])=CC=CC3=CC2=C1 RHZUVFJBSILHOK-UHFFFAOYSA-N 0.000 description 2
- 239000003830 anthracite Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
Landscapes
- Manufacture And Refinement Of Metals (AREA)
- Air Supply (AREA)
Abstract
本发明公开了一种基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，包括以下步骤：步骤1：采集熟料烧成系统的DCS系统运行参数数据及检化验设备中的生料、熟料化学分析数据；步骤2：判断熟料烧成系统是否处于稳定运行状态，当熟料烧成系统处于稳定运行状态时，执行步骤3；否则，结束计算过程；步骤3：根据熟料烧成系统的燃煤燃烧计算模型及热平衡分析模型，计算燃煤低位热值；步骤4：根据步骤3计算的燃煤低位热值，计算熟料烧成热耗、熟料烧成热效率。本发明克服了燃煤煤质波动给熟料烧成耗热计算带来的不利影响，计算简便，成本低，具有良好的可实施性和可操作性。
Description
技术领域
本发明涉及水泥行业熟料烧成系统烧成热耗、热效率及燃煤热值的实时监测方法，属于软测量领域。
背景技术
水泥熟料烧成系统是新型干法水泥生产的主要部分，其中的能源类型主要为燃煤燃烧释放和熟料生成所消耗的化学能，以及系统各部分耗散的热能，作为一个耗能巨大且复杂的热力学和化学系统，其能耗的高低直接影响和决定整个水泥生产线的能源利用率和经济效益。
目前，熟料烧成系统操作人员主要依靠经验调整回转窑运行工况，当出现燃煤热值波动、生料化学成分波动等情况时，难以及时、准确地调整相关运行参数，导致烧成热耗增加。同时，水泥生产企业对燃煤热值的监测主要是通过定期离线取样化验获得，其与实际进入系统的燃煤之间存在着较大的制样误差及严重的分析时间滞后。
因此，如何在线测算燃煤热值，并在此基础实现熟料烧成热耗、烧成系统热效率的实时监测，指导操作人员及时、准确地调整烧成系统运行工况，保证烧成系统稳定、经济运行，是相关技术人员普遍关心的问题。
中国发明专利201310463399.2提出了一种水泥分解炉入炉煤粉低位热值在线估计方法，通过实时监测分解炉出口温度、分解炉喷煤量、生料下料量、三次风温和三次风量的变化来表征煤粉低位热值的变化，进而计算煤粉低位热值的估计值，但其所建立的煤粉低位热值变化计算模型(详见说明书中式(13))中，并未考虑三次风量及三次风温的变化，导致计算出的热值变化仅取决于分解炉出口温度变化，计算结果不够准确，且该方法仅能计算出热值变化值，如何获取热值初始值，在该专利中并未提及，故该方法具有一定的局限性。
发明内容
本发明的目的在于提供了基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，该方法利用燃煤燃烧计算及熟料烧成系统热平衡计算，实现了燃煤热值与熟料烧成热耗、热效率的在线监测。
为达成上述目的，本发明的技术方案是：基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于包括以下步骤：
步骤1：采集熟料烧成系统的DCS系统运行参数数据及检化验设备中的生料、熟料化学分析数据；
步骤2：判断熟料烧成系统是否处于稳定运行状态，当熟料烧成系统处于稳定运行状态时，执行步骤3；否则，结束计算过程；
步骤3：根据熟料烧成系统的燃煤燃烧计算模型及热平衡分析模型，计算燃煤低位热值；
步骤4：根据步骤3计算的燃煤低位热值，计算熟料烧成热耗、熟料烧成热效率。
作为本发明进一步改进的技术方案，上述判断熟料烧成系统处于稳定运行状态判据为：
1)回转窑主电机电流在其历史数据平均值的±30％范围内波动；
2)窑尾气体温度在1050～1200℃之间波动；
同时满足上述两点判据时，认为熟料烧成系统处于稳定运行状态；否则，认为熟料烧成系统处于不稳定状态。
作为本发明进一步改进的技术方案，上述的熟料烧成系统的热平衡分析基准分别为：温度基准：0℃；质量基准：1kg熟料；平衡范围是从回转窑熟料出口到预热器废气出口。
作为本发明进一步改进的技术方案，所述的熟料烧成系统的燃煤燃烧计算模型包括：理论空气量计算、过量空气系数计算、入窑实际空气量计算、三次空气量计算、系统总漏风量计算以及总烟气量计算。
作为本发明进一步改进的技术方案，所述的熟料烧成系统的燃煤燃烧计算模型具体为：
1)理论空气量计算子模型
Vlk0＝K1mrQar.net/1000 (1)
式中：Vlk0——每千克熟料燃煤燃烧理论空气量，Nm3/kg；
Qar,net——燃煤低位热值，kJ/kg煤，待计算；
K1——由燃煤热值计算理论空气量系数，可依煤种选取，对于烟煤，取0.2620；对于无烟煤，取0.2659；
mr——每千克熟料燃煤消耗量，kg/kg；
Mfr——分解炉瞬时喂煤量，kg/h；
Myr——窑瞬时喂煤量，kg/h；
Ms——生料瞬时喂料量，kg/h；
2)过量空气系数计算子模型
式中：α——过量空气系数；
O2——烟气中的含氧量，％；
——每千克熟料生料分解生成的CO2体积量，Nm3/kg；
ε——碳酸盐分解率，％，窑尾处可取0.05～0.25，分解炉出口及预热器一级筒出口取为1；
ms——每千克熟料的生料消耗量，kg/kg，可取统计值；
CaOs、MgOs——生料中的CaO、MgO百分含量，％；
3)入窑实际空气量计算子模型
Vyk＝αywVlk0kyr (5)
式中：Vyk——每千克熟料实际入窑干空气体积量，Nm3/kg；
αyw——窑尾处的过量空气系数，由窑尾处的烟气含氧量O2yw按(3)式计算；
kyr——入窑燃煤量占总燃煤量的比例；
根据入窑实际干空气量，分别计算入窑一次空气量、二次空气量及漏风量：
Vyk1＝kyk1Vyk (6)
Vyk2＝kyk2Vyk (7)
Vlok1＝klok1Vyk (8)
式中：Vyk1、Vyk2、Vlok1——分别为每千克熟料实际入窑一次空气体积量、二次空气体积量以及漏风量，Nm3/kg；
kyk1、kyk2、klok1——分别为入窑一次空气、二次空气及漏风量的风量比，可按设计值或热平衡试验值选取，且应满足kyk1+kyk2+klok1＝1；
4)三次空气量计算子模型
Vf3k＝V1+V2-V3-V4 (9)
V1＝(αf-1)Vlk0 (10)
V2＝Vlk0kfr (11)
V3＝(αyw-1)Vlk0kyr (12)
V4＝klok4Vlk0kfr (13)
式中：Vf3k——每千克熟料三次空气量，Nm3/kg；
V1、V2、V3、V4——分别为每千克熟料分解炉出口处的过量空气量，分解炉燃煤燃烧所需空气量，窑尾处的过量空气量以及分解炉及窑尾漏入空气量，Nm3/kg；
αf——分解炉出口处的过量空气系数，由窑尾处的烟气含氧量O2f按(3)式计算；
kfr——入分解炉燃煤量占总燃煤量的比例；
klok4——分解炉及窑尾漏风占分解炉燃煤燃烧所需空气量的比例，可按设计值或热平衡试验值选取，一般取为0.02～0.06；
5)系统漏风计算子模型
Vlok＝Vlok1+Vlok2 (14)
Vlok2＝V4+V5 (15)
V5＝klok5Vlk0 (16)
式中：Vlok——每千克熟料系统漏风量，Nm3/kg；
Vlok2——每千克熟料窑尾漏入空气量，Nm3/kg；
V5——每千克熟料预热器漏入空气量，Nm3/kg；
klok5——预热器漏风占理论空气量的比例，可按设计值或热平衡试验值选取，一般取为0.05～0.15；
6)预热器出口总烟气量计算
Vgy0＝K2Vlk0 (18)
Vk＝(αy-1)Vlk0 (22)
式中：Vf——每千克熟料预热器出口总烟气量，Nm3/kg；
Vgy0——每千克熟料燃料燃烧生成的理论干烟气量，Nm3/kg；
——每千克熟料燃料燃烧生成的H2O量，Nm3/kg；
——每千克熟料生料中物理含水释放出的H2O量，Nm3/kg；
——每千克熟料生料中化学含水释放出的H2O量，Nm3/kg；
Vk——每千克熟料烟气中的过量空气量，Nm3/kg；
K2——燃煤燃烧所需的理论干空气量与生成的理论干烟气量之间的折算系数，在0.90～0.98之间取值；
——燃煤收到基H元素含量及水分，％，对于烟煤，可取值为3，按工业分析值选取；
Ws——生料水分，％；
Al2O3 s——生料中的Al2O3百分含量，％；
αy——预热器出口处的过量空气系数，由预热器出口处的烟气含氧量O2y按(3)式计算。
作为本发明进一步改进的技术方案，熟料烧成系统收入总热量Qzs计算如下：
Qzs＝QrR+Qr+Qs+Qyh+Qy1k+Qy2k+Qf3k+Qlok (23)
QrR＝mrQar,net (24)
Qr＝mrcrtr (25)
Qs＝mscsts (26)
Qyh＝myhcyhtyh (27)
Qy1k＝Vy1kcy1kty1k (28)
Qy1k＝Vy2kcy2kty2k (29)
Qf3k＝Vf3kcf3ktf3k (30)
Qlok＝Vlokcloktlok (31)
式中：QrR——每千克熟料燃煤燃烧热，kJ/kg；
Qr——每千克熟料燃煤带入显热，kJ/kg；
Qs——每千克熟料生料带入显热，kJ/kg；
Qyh——每千克熟料入窑回灰显热，kJ/kg；
Qy1k——每千克熟料入窑一次空气带入显热，kJ/kg；
Qy2k——每千克熟料入窑二次空气带入显热，kJ/kg；
Qf3k——每千克熟料入分解炉三次空气带入显热，kJ/kg；
Qlok——每千克熟料系统漏风带入显热，kJ/kg；
cr——燃煤比热，kJ/(kg·℃)，一般取为1.154kJ/(kg·℃)；
cs——生料比热，kJ/(kg·℃)，一般取为0.878kJ/(kg·℃)；
cyh——回灰比热，kJ/(kg·℃)，一般取为0.836kJ/(kg·℃)；
tr——燃煤温度，℃，若无计量，一般取为60℃；
ts——生料温度，℃，若无计量，一般取为50℃；
tyh——回灰温度，℃，若无计量，一般取为50℃；
cy1k、cy2k、cf3k、clok——分别为一次空气比热、二次空气比热、三次空气比热及漏风比热，kJ/(kg·℃)，一般可分别取为1.298、1.403、1.377及1.298；
ty1k、ty2k、tf3k、tlok——分别为一次空气温度、二次空气温度、三次空气温度及漏风温度，℃。
作为本发明进一步改进的技术方案，熟料烧成系统支出总热量Qzs计算如下：
Qzc＝Qsh+Qss+Qf+Qysh+Qfh+Qjb+Qb (32)
Qsh＝17.19Al2O3 sh+27.10MgOsh+32.01CaOsh-21.40SiO2 sh-2.47Fe2O3 sh
(33)
Qf＝Vfcftpy (35)
Qysh＝cshtsh (36)
Qfh＝myhcfhtpy (37)
式中：Qsh——每千克熟料形成热，kJ/kg；
Qss——每千克熟料蒸发生料中的水分耗热量，kJ/kg；
Qf——每千克熟料废气带走热量，kJ/kg；
Qysh——每千克熟料出窑熟料带走热量，kJ/kg；
Qfh——每千克熟料出预热器飞灰带走热量，kJ/kg；
Qjb——每千克熟料机械不完全燃烧热损失，kJ/kg；
Qb——每千克熟料散热损失，kJ/kg，熟料烧成系统稳定运行时，散热损失占总热量的比例基本保持不变，即其中值可根据设计值或试验值选取；
Al2O3 sh、MgOsh、CaOsh、SiO2 sh Fe2O3 sh——熟料中相应成分的质量百分数，％；
qqh——水的汽化潜热，一般取值2380kJ/kg；
cf——预热器出口废气比热，kJ/(kg·℃)，一般取为1.4750kJ/(kg·℃)；
csh——出窑熟料比热，kJ/(kg·℃)，一般取为1.078kJ/(kg·℃)；
cfh——出预热器飞灰比热，kJ/(kg·℃)，一般取为0.895kJ/(kg·℃)；
tpy——预热器出口废气温度，℃；
tsh——熟料出窑温度，℃；
作为本发明进一步改进的技术方案，所述的熟料烧成系统热平衡分析模型为：
Qzs＝Qzc (39)
式中：Qzs——烧成系统收入总热量，kJ/kg；
Qzc——烧成系统支出总热量，kJ/kg；
将烧成系统燃煤燃烧计算模型及各项收入热量与支出热量计算模型(式(1)～(38))带入式(39)，组成关于燃煤热值Qar,net的一元一次方程，进而求解出当前状态下的燃煤热值。
作为本发明进一步改进的技术方案，求解出燃煤热值Qar,net之后，熟料烧成热耗及熟料烧成热效率计算如下：
QrR＝mrQar,net (40)
式中：QrR——单位熟料烧成热耗，kJ/kg；
Qsh——每千克熟料形成热，kJ/kg；
ηs——烧成系统热效率，％。
作为本发明进一步改进的技术方案，检化验设备包括荧光分析仪。
本发明利用相关数据接口，从水泥熟料烧成系统DCS系统读取运行参数数据，从荧光分析仪等检化验系统获取生料、熟料化学分析数据，借助熟料烧成系统的热平衡分析模型，实时计算燃煤低位热值，并在此基础上，实现熟料烧成热耗、热效率的在线监测。本发明克服了燃煤煤质波动给熟料烧成耗热计算带来的不利影响，计算简便，成本低，具有良好的可实施性和可操作性。
本发明的优点在于：
1、作为一种软测量方法，熟料烧成热耗、热效率及燃煤热值的测算过程中所需的参数均可从DCS(或SIS)系统的实时数据库中直接读取，现场不需要额外增加分析或测量仪表灯昂贵的辅助设备，可集成在已有的DCS(或SIS)系统中，成本低。
2、本方法无需参考煤质成分之间的关联关系，只需满足烧成系统的热平衡关系，可对不同煤质的燃煤热值进行测算，适用性广。
附图说明
图1是水泥熟料烧成系统热平衡原理图；
图2是基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法流程图；
图3是实施本发明提供方法的系统结构图。
以下将结合附图，对本发明的技术方案进行详细说明。应强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
具体实施方式
实施例1
本发明以水泥熟料烧成系统热平衡原理为理论基础。图1是熟料烧成系统热平衡原理图，如图1所示，熟料烧成系统的收入总热量应当等于支出总热量。
图2是基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法流程图。
图3是实施本发明提供方法的系统结构图。如图3所示，在实时本发明之前，先要搭建硬件系统。配置一台安全隔离网关，用于采集熟料烧成DCS系统数据，并与DCS系统进行安全隔离；一台IO服务器用于数据采集及实时数据库；一台应用服务器用于运行计算程序；一台WEB服务器用户向企业内网发布计算展示画面及报表；一台工作站用于实时监视。
参见图1、图2和图3，本基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，包括以下步骤：
步骤1：采集熟料烧成系统的DCS系统运行参数数据及检化验设备中的生料、熟料化学分析数据；
步骤2：判断熟料烧成系统是否处于稳定运行状态，当熟料烧成系统处于稳定运行状态时，执行步骤3；否则，结束计算过程；
步骤3：根据熟料烧成系统的燃煤燃烧计算模型及热平衡分析模型，计算燃煤低位热值；
步骤4：根据步骤3计算的燃煤低位热值，计算熟料烧成热耗、熟料烧成热效率。
作为优选方案，上述判断熟料烧成系统处于稳定运行状态判据为：
1)回转窑主电机电流在其历史数据平均值的±30％范围内波动；
2)窑尾气体温度在1050～1200℃之间波动；
同时满足上述两点判据时，认为熟料烧成系统处于稳定运行状态；否则，认为熟料烧成系统处于不稳定状态。
上述的熟料烧成系统的热平衡分析基准分别为：温度基准：0℃；质量基准：1kg熟料；平衡范围是从回转窑熟料出口到预热器废气出口。所述的熟料烧成系统的燃煤燃烧计算模型包括：理论空气量计算、过量空气系数计算、入窑实际空气量计算、三次空气量计算、系统总漏风量计算以及总烟气量计算。10.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于：检化验设备包括荧光分析仪。
所述的熟料烧成系统的燃煤燃烧计算模型具体为：
1)理论空气量计算子模型
Vlk0＝K1mrQar.net/1000 (1)
式中：Vlk0——每千克熟料燃煤燃烧理论空气量，Nm3/kg；
Qar,net——燃煤低位热值，kJ/kg煤，待计算；
K1——由燃煤热值计算理论空气量系数，可依煤种选取，对于烟煤，取0.2620；对于无烟煤，取0.2659；
mr——每千克熟料燃煤消耗量，kg/kg；
Mfr——分解炉瞬时喂煤量，kg/h；
Myr——窑瞬时喂煤量，kg/h；
Ms——生料瞬时喂料量，kg/h；
2)过量空气系数计算子模型
式中：α——过量空气系数；
O2——烟气中的含氧量，％；
——每千克熟料生料分解生成的CO2体积量，Nm3/kg；
ε——碳酸盐分解率，％，窑尾处可取0.05～0.25，分解炉出口及预
热器一级筒出口取为1；
ms——每千克熟料的生料消耗量，kg/kg，可取统计值；
CaOs、MgOs——生料中的CaO、MgO百分含量，％；
3)入窑实际空气量计算子模型
Vyk＝αywVlk0kyr (5)
式中：Vyk——每千克熟料实际入窑干空气体积量，Nm3/kg；
αyw——窑尾处的过量空气系数，由窑尾处的烟气含氧量O2yw按(3)式计算；
kyr——入窑燃煤量占总燃煤量的比例；
根据入窑实际干空气量，分别计算入窑一次空气量、二次空气量及漏风量：
Vyk1＝kyk1Vyk (6)
Vyk2＝kyk2Vyk (7)
Vlok1＝klok1Vyk (8)
式中：Vyk1、Vyk2、Vlok1——分别为每千克熟料实际入窑一次空气体积量、二次空气体积量以及漏风量，Nm3/kg；
kyk1、kyk2、klok1——分别为入窑一次空气、二次空气及漏风量的风量比，可按设计值或热平衡试验值选取，且应满足kyk1+kyk2+klok1＝1；
4)三次空气量计算子模型
Vf3k＝V1+V2-V3-V4 (9)
V1＝(αf-1)Vlk0 (10)
V2＝Vlk0kfr (11)
V3＝(αyw-1)Vlk0kyr (12)
V4＝klok4Vlk0kfr(13)
式中：Vf3k——每千克熟料三次空气量，Nm3/kg；
V1、V2、V3、V4——分别为每千克熟料分解炉出口处的过量空气量，分解炉燃煤燃烧所需空气量，窑尾处的过量空气量以及分解炉及窑尾漏入空气量，Nm3/kg；
αf——分解炉出口处的过量空气系数，由窑尾处的烟气含氧量O2f按(3)式计算；
kfr——入分解炉燃煤量占总燃煤量的比例；
klok4——分解炉及窑尾漏风占分解炉燃煤燃烧所需空气量的比例，可按设计值或热平衡试验值选取，一般取为0.02～0.06；
5)系统漏风计算子模型
Vlok＝Vlok1+Vlok2 (14)
Vlok2＝V4+V5 (15)
V5＝klok5Vlk0 (16)
式中：Vlok——每千克熟料系统漏风量，Nm3/kg；
Vlok2——每千克熟料窑尾漏入空气量，Nm3/kg；
V5——每千克熟料预热器漏入空气量，Nm3/kg；
klok5——预热器漏风占理论空气量的比例，可按设计值或热平衡试验值选取，一般取为0.05～0.15；
6)预热器出口总烟气量计算
Vgy0＝K2Vlk0 (18)
Vk＝(αy-1)Vlk0 (22)
式中：Vf——每千克熟料预热器出口总烟气量，Nm3/kg；
Vgy0——每千克熟料燃料燃烧生成的理论干烟气量，Nm3/kg；
——每千克熟料燃料燃烧生成的H2O量，Nm3/kg；
——每千克熟料生料中物理含水释放出的H2O量，Nm3/kg；
——每千克熟料生料中化学含水释放出的H2O量，Nm3/kg；
Vk——每千克熟料烟气中的过量空气量，Nm3/kg；
K2——燃煤燃烧所需的理论干空气量与生成的理论干烟气量之间的折算系数，在0.90～0.98之间取值；
——燃煤收到基H元素含量及水分，％，对于烟煤，可取值为3，按工业分析值选取；
Ws——生料水分，％；
Al2O3 s——生料中的Al2O3百分含量，％；
αy——预热器出口处的过量空气系数，由预热器出口处的烟气含氧量O2y按(3)式计算。
熟料烧成系统收入总热量Qzs计算如下：
Qzs＝QrR+Qr+Qs+Qyh+Qy1k+Qy2k+Qf3k+Qlok (23)
QrR＝mrQar,net (24)
Qr＝mrcrtr (25)
Qs＝mscsts (26)
Qyh＝myhcyhtyh (27)
Qy1k＝Vy1kcy1kty1k (28)
Qy1k＝Vy2kcy2kty2k (29)
Qf3k＝Vf3kcf3ktf3k (30)
Qlok＝Vlokcloktlok (31)
式中：QrR——每千克熟料燃煤燃烧热，kJ/kg；
Qr——每千克熟料燃煤带入显热，kJ/kg；
Qs——每千克熟料生料带入显热，kJ/kg；
Qyh——每千克熟料入窑回灰显热，kJ/kg；
Qy1k——每千克熟料入窑一次空气带入显热，kJ/kg；
Qy2k——每千克熟料入窑二次空气带入显热，kJ/kg；
Qf3k——每千克熟料入分解炉三次空气带入显热，kJ/kg；
Qlok——每千克熟料系统漏风带入显热，kJ/kg；
cr——燃煤比热，kJ/(kg·℃)，一般取为1.154kJ/(kg·℃)；
cs——生料比热，kJ/(kg·℃)，一般取为0.878kJ/(kg·℃)；
cyh——回灰比热，kJ/(kg·℃)，一般取为0.836kJ/(kg·℃)；
tr——燃煤温度，℃，若无计量，一般取为60℃；
ts——生料温度，℃，若无计量，一般取为50℃；
tyh——回灰温度，℃，若无计量，一般取为50℃；
cy1k、cy2k、cf3k、clok——分别为一次空气比热、二次空气比热、三次空气比热及漏风比热，kJ/(kg·℃)，一般可分别取为1.298、1.403、1.377及1.298；
ty1k、ty2k、tf3k、tlok——分别为一次空气温度、二次空气温度、三次空气温度及漏风温度，℃。
熟料烧成系统支出总热量Qzs计算如下：
Qzc＝Qsh+Qss+Qf+Qysh+Qfh+Qjb+Qb (32)
Qsh＝17.19Al2O3 sh+27.10MgOsh+32.01CaOsh-21.40SiO2 sh-2.47Fe2O3 sh
(33)
Qf＝Vfcftpy (35)
Qysh＝cshtsh (36)
Qfh＝myhcfhtpy (37)
式中：Qsh——每千克熟料形成热，kJ/kg；
Qss——每千克熟料蒸发生料中的水分耗热量，kJ/kg；
Qf——每千克熟料废气带走热量，kJ/kg；
Qysh——每千克熟料出窑熟料带走热量，kJ/kg；
Qfh——每千克熟料出预热器飞灰带走热量，kJ/kg；
Qjb——每千克熟料机械不完全燃烧热损失，kJ/kg；
Qb——每千克熟料散热损失，kJ/kg，熟料烧成系统稳定运行时，散热损失占总热量的比例基本保持不变，即其中值可根据设计值或试验值选取；
Al2O3 sh、MgOsh、CaOsh、SiO2 sh Fe2O3 sh——熟料中相应成分的质量百分数，％；
qqh——水的汽化潜热，一般取值2380kJ/kg；
cf——预热器出口废气比热，kJ/(kg·℃)，一般取为1.4750kJ/(kg·℃)；
csh——出窑熟料比热，kJ/(kg·℃)，一般取为1.078kJ/(kg·℃)；
cfh——出预热器飞灰比热，kJ/(kg·℃)，一般取为0.895kJ/(kg·℃)；
tpy——预热器出口废气温度，℃；
tsh——熟料出窑温度，℃；
Lsh——熟料烧失量，％。
所述的熟料烧成系统热平衡分析模型为：
Qzs＝Qzc (39)
式中：Qzs——烧成系统收入总热量，kJ/kg；
Qzc——烧成系统支出总热量，kJ/kg；
将烧成系统燃煤燃烧计算模型及各项收入热量与支出热量计算模型(式(1)～(38))带入式(39)，组成关于燃煤热值Qar,net的一元一次方程，进而求解出当前状态下的燃煤热值。
求解出燃煤热值Qar,net之后，熟料烧成热耗及熟料烧成热效率计算如下：
QrR＝mrQar,net (40)
式中：QrR——单位熟料烧成热耗，kJ/kg；
Qsh——每千克熟料形成热，kJ/kg；
ηs——烧成系统热效率，％。
Claims (9)
1.基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于包括以下步骤：
步骤1：采集熟料烧成系统的DCS系统运行参数数据及检化验设备中的生料、熟料化学分析数据；
步骤2：判断熟料烧成系统是否处于稳定运行状态，当熟料烧成系统处于稳定运行状态时，执行步骤3；否则，结束计算过程；
步骤3：根据熟料烧成系统的燃煤燃烧计算模型及热平衡分析模型，计算燃煤低位热值；
所述的熟料烧成系统的燃煤燃烧计算模型具体为：
1)理论空气量计算子模型
Vlk0＝K1mrQar.net/1000 (1)
<mrow>
<msub>
<mi>m</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>f</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>y</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mi>s</mi>
</msub>
<mo>/</mo>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：Vlk0——每千克熟料燃煤燃烧理论空气量，Nm3/kg；
Qar,net——燃煤低位热值，kJ/kg煤，待计算；
K1——由燃煤热值计算理论空气量系数，可依煤种选取，对于烟煤，取0.2620；对于无烟煤，取0.2659；
mr——每千克熟料燃煤消耗量，kg/kg；
Mfr——分解炉瞬时喂煤量，kg/h；
Myr——窑瞬时喂煤量，kg/h；
Ms——生料瞬时喂料量，kg/h；
2)过量空气系数计算子模型
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>21</mn>
<mo>+</mo>
<msub>
<mi>O</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&epsiv;V</mi>
<mrow>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<mn>0.02</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>21</mn>
<mo>-</mo>
<msub>
<mi>O</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<mn>1.977</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<mi>CaO</mi>
<mi>s</mi>
</msup>
<mfrac>
<mn>44</mn>
<mn>56</mn>
</mfrac>
<mo>+</mo>
<msup>
<mi>MgO</mi>
<mi>s</mi>
</msup>
<mfrac>
<mn>44</mn>
<mn>40.3</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：α——过量空气系数；
O2——烟气中的含氧量，％；
——每千克熟料生料分解生成的CO2体积量，Nm3/kg；
ε——碳酸盐分解率，％，窑尾处可取0.05～0.25，分解炉出口及预热器一级筒出口取为1；
ms——每千克熟料的生料消耗量，kg/kg，可取统计值；
CaOs、MgOs——生料中的CaO、MgO百分含量，％；
3)入窑实际空气量计算子模型
Vyk＝αywVlk0kyr (5)
式中：Vyk——每千克熟料实际入窑干空气体积量，Nm3/kg；
αyw——窑尾处的过量空气系数，由窑尾处的烟气含氧量O2yw按(3)式计算；
kyr——入窑燃煤量占总燃煤量的比例；
根据入窑实际干空气量，分别计算入窑一次空气量、二次空气量及漏风量：
Vyk1＝kyk1Vyk (6)
Vyk2＝kyk2Vyk (7)
Vlok1＝klok1Vyk (8)
式中：Vyk1、Vyk2、Vlok1——分别为每千克熟料实际入窑一次空气体积量、二次空气体积量以及漏风量，Nm3/kg；
kyk1、kyk2、klok1——分别为入窑一次空气、二次空气及漏风量的风量比，可按设计值或热平衡试验值选取，且应满足kyk1+kyk2+klok1＝1；
4)三次空气量计算子模型
Vf3k＝V1+V2-V3-V4 (9)
V1＝(αf-1)Vlk0 (10)
V2＝Vlk0kfr (11)
V3＝(αyw-1)Vlk0kyr (12)
V4＝klok4Vlk0kfr (13)
式中：Vf3k——每千克熟料三次空气量，Nm3/kg；
V1、V2、V3、V4——分别为每千克熟料分解炉出口处的过量空气量，分解炉燃煤燃烧所需空气量，窑尾处的过量空气量以及分解炉及窑尾漏入空气量，Nm3/kg；
αf——分解炉出口处的过量空气系数，由窑尾处的烟气含氧量O2f按(3)式计算；
kfr——入分解炉燃煤量占总燃煤量的比例；
klok4——分解炉及窑尾漏风占分解炉燃煤燃烧所需空气量的比例，可按设计值或热平衡试验值选取，一般取为0.02～0.06；
5)系统漏风计算子模型
Vlok＝Vlok1+Vlok2 (14)
Vlok2＝V4+V5 (15)
V5＝klok5Vlk0 (16)
式中：Vlok——每千克熟料系统漏风量，Nm3/kg；
Vlok2——每千克熟料窑尾漏入空气量，Nm3/kg；
V5——每千克熟料预热器漏入空气量，Nm3/kg；
klok5——预热器漏风占理论空气量的比例，可按设计值或热平衡试验值选取，一般取为0.05～0.15；
6)预热器出口总烟气量计算
<mrow>
<msub>
<mi>V</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>g</mi>
<mi>y</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mi>O</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>CO</mi>
<mn>2</mn>
</msub>
</mrow>
<mi>s</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
Vgy0＝K2Vlk0 (18)
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mi>O</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<mn>1.24</mn>
<mrow>
<mo>(</mo>
<mn>9</mn>
<msubsup>
<mi>H</mi>
<mrow>
<mi>a</mi>
<mi>r</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>a</mi>
<mi>r</mi>
</mrow>
<mi>r</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mn>100</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<msub>
<mi>W</mi>
<mi>s</mi>
</msub>
</mrow>
<mn>80.4</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mn>0.00353</mn>
<msub>
<mi>m</mi>
<mi>s</mi>
</msub>
<msub>
<mi>Al</mi>
<mn>2</mn>
</msub>
<msup>
<msub>
<mi>O</mi>
<mn>3</mn>
</msub>
<mi>s</mi>
</msup>
</mrow>
<mn>80.4</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
Vk＝(αy-1)Vlk0 (22)
式中：Vf——每千克熟料预热器出口总烟气量，Nm3/kg；
Vgy0——每千克熟料燃料燃烧生成的理论干烟气量，Nm3/kg；
——每千克熟料燃料燃烧生成的H2O量，Nm3/kg；
——每千克熟料生料中物理含水释放出的H2O量，Nm3/kg；
——每千克熟料生料中化学含水释放出的H2O量，Nm3/kg；
Vk——每千克熟料烟气中的过量空气量，Nm3/kg；
K2——燃煤燃烧所需的理论干空气量与生成的理论干烟气量之间的折算系数，在0.90～0.98之间取值；
——燃煤收到基H元素含量及水分，％，对于烟煤，可取值为3，按工业分析值选取；
Ws——生料水分，％；
Al2O3 s——生料中的Al2O3百分含量，％；
αy——预热器出口处的过量空气系数，由预热器出口处的烟气含氧量O2y按(3)式计算；
步骤4：根据步骤3计算的燃煤低位热值，计算熟料烧成热耗、熟料烧成热效率。
2.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于上述判断熟料烧成系统处于稳定运行状态判据为：
1)回转窑主电机电流在其历史数据平均值的±30％范围内波动；
2)窑尾气体温度在1050～1200℃之间波动；
同时满足上述两点判据时，认为熟料烧成系统处于稳定运行状态；否则，认为熟料烧成系统处于不稳定状态。
3.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于上述的熟料烧成系统的热平衡分析基准分别为：温度基准：0℃；质量基准：1kg熟料；平衡范围是从回转窑熟料出口到预热器废气出口。
4.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于所述的熟料烧成系统的燃煤燃烧计算模型包括：理论空气量计算、过量空气系数计算、入窑实际空气量计算、三次空气量计算、系统总漏风量计算以及总烟气量计算。
5.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于熟料烧成系统收入总热量Qzs计算如下：
Qzs＝QrR+Qr+Qs+Qyh+Qy1k+Qy2k+Qf3k+Qlok (23)
QrR＝mrQar,net (24)
Qr＝mrcrtr (25)
Qs＝mscsts (26)
Qyh＝myhcyhtyh (27)
Qy1k＝Vy1kcy1kty1k (28)
Qy1k＝Vy2kcy2kty2k (29)
Qf3k＝Vf3kcf3ktf3k (30)
Qlok＝Vlokcloktlok (31)
式中：QrR——每千克熟料燃煤燃烧热，kJ/kg；
Qr——每千克熟料燃煤带入显热，kJ/kg；
Qs——每千克熟料生料带入显热，kJ/kg；
Qyh——每千克熟料入窑回灰显热，kJ/kg；
Qy1k——每千克熟料入窑一次空气带入显热，kJ/kg；
Qy2k——每千克熟料入窑二次空气带入显热，kJ/kg；
Qf3k——每千克熟料入分解炉三次空气带入显热，kJ/kg；
Qlok——每千克熟料系统漏风带入显热，kJ/kg；
cr——燃煤比热，kJ/(kg·℃)，一般取为1.154kJ/(kg·℃)；
cs——生料比热，kJ/(kg·℃)，一般取为0.878kJ/(kg·℃)；
cyh——回灰比热，kJ/(kg·℃)，一般取为0.836kJ/(kg·℃)；
tr——燃煤温度，℃，若无计量，一般取为60℃；
ts——生料温度，℃，若无计量，一般取为50℃；
tyh——回灰温度，℃，若无计量，一般取为50℃；
cy1k、cy2k、cf3k、clok——分别为一次空气比热、二次空气比热、三次空气比热及漏风比热，kJ/(kg·℃)，一般可分别取为1.298、1.403、1.377及1.298；
ty1k、ty2k、tf3k、tlok——分别为一次空气温度、二次空气温度、三次空气温度及漏风温度，℃。
6.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于熟料烧成系统支出总热量Qzs计算如下：
Qzc＝Qsh+Qss+Qf+Qysh+Qfh+Qjb+Qb (32)
Qsh＝17.19Al2O3 sh+27.10MgOsh+32.01CaOsh-21.40SiO2 sh-2.47Fe2O3 sh
(33)
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>s</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0.804</mn>
<mrow>
<mo>(</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>p</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>w</mi>
<mi>c</mi>
</mrow>
<mi>s</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>q</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>34</mn>
<mo>)</mo>
</mrow>
</mrow>
Qf＝Vfcftpy (35)
Qysh＝cshtsh (36)
Qfh＝myhcfhtpy (37)
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>j</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mn>100</mn>
</mfrac>
<mn>33874</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>38</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：Qsh——每千克熟料形成热，kJ/kg；
Qss——每千克熟料蒸发生料中的水分耗热量，kJ/kg；
Qf——每千克熟料废气带走热量，kJ/kg；
Qysh——每千克熟料出窑熟料带走热量，kJ/kg；
Qfh——每千克熟料出预热器飞灰带走热量，kJ/kg；
Qjb——每千克熟料机械不完全燃烧热损失，kJ/kg；
Qb——每千克熟料散热损失，kJ/kg，熟料烧成系统稳定运行时，散热损失占总热量的比例基本保持不变，即其中值可根据设计值或试验值选取；
Al2O3 sh、MgOsh、CaOsh、SiO2 shFe2O3 sh——熟料中相应成分的质量百分数，％；
qqh——水的汽化潜热，一般取值2380kJ/kg；
cf——预热器出口废气比热，kJ/(kg·℃)，一般取为1.4750kJ/(kg·℃)；
csh——出窑熟料比热，kJ/(kg·℃)，一般取为1.078kJ/(kg·℃)；
cfh——出预热器飞灰比热，kJ/(kg·℃)，一般取为0.895kJ/(kg·℃)；
tpy——预热器出口废气温度，℃；
tsh——熟料出窑温度，℃；
Lsh——熟料烧失量，％。
7.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于所述的熟料烧成系统热平衡分析模型为：
Qzs＝Qzc (39)
式中：Qzs——烧成系统收入总热量，kJ/kg；
Qzc——烧成系统支出总热量，kJ/kg；
将烧成系统燃煤燃烧计算模型及各项收入热量与支出热量计算模型(式(1)～(38))带入式(39)，组成关于燃煤热值Qar,net的一元一次方程，进而求解出当前状态下的燃煤热值。
8.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于，求解出燃煤热值Qar,net之后，熟料烧成热耗及熟料烧成热效率计算如下：
QrR＝mrQar,net (40)
<mrow>
<msub>
<mi>&eta;</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>Q</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<msub>
<mi>Q</mi>
<mrow>
<mi>r</mi>
<mi>R</mi>
</mrow>
</msub>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>41</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：QrR——单位熟料烧成热耗，kJ/kg；
Qsh——每千克熟料形成热，kJ/kg；
ηs——烧成系统热效率，％。
9.根据权利要求1所述的基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法，其特征在于：检化验设备包括荧光分析仪。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510342228.3A CN104915567B (zh) | 2015-06-18 | 2015-06-18 | 基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510342228.3A CN104915567B (zh) | 2015-06-18 | 2015-06-18 | 基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104915567A CN104915567A (zh) | 2015-09-16 |
| CN104915567B true CN104915567B (zh) | 2017-10-27 |
Family
ID=54084628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510342228.3A Expired - Fee Related CN104915567B (zh) | 2015-06-18 | 2015-06-18 | 基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104915567B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109033511B (zh) * | 2018-06-15 | 2019-06-07 | 武汉理工大学 | 一种结合数据驱动与数据挖掘的水泥窑系统热耗分析方法 |
| CN109374046A (zh) * | 2018-10-12 | 2019-02-22 | 安徽海螺建材设计研究院有限责任公司 | 一种水泥熟料线烧成系统在线智能化热工标定的方法 |
| CN110991090B (zh) * | 2019-12-26 | 2023-04-11 | 广东工业大学 | 一种马蹄焰玻璃窑炉数据预处理及热效率计算方法及系统 |
| CN111174569B (zh) * | 2020-01-16 | 2021-10-22 | 武汉科技大学 | 一种在线预测回转窑内煅烧段烟气温度的方法及系统 |
| CN113111511B (zh) * | 2021-04-13 | 2022-10-28 | 厦门信达环能科技有限公司 | 一种基于烟气成分分析的水泥熟料烧成热耗计算方法 |
| CN113606932B (zh) * | 2021-08-13 | 2022-07-01 | 济南大学 | 水泥回转窑热效率在线计算参数的评估及校正方法和装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559382A (zh) * | 2013-10-09 | 2014-02-05 | 济南大学 | 一种水泥分解炉入炉煤粉低位热值在线估计方法 |
| CN104061589A (zh) * | 2014-07-17 | 2014-09-24 | 烟台龙源电力技术股份有限公司 | 适应煤质变化的低氮燃烧控制方法和装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5342355B2 (ja) * | 2009-07-22 | 2013-11-13 | 株式会社神戸製鋼所 | ボイラの灰付着抑制方法及び灰付着抑制装置 |
-
2015
- 2015-06-18 CN CN201510342228.3A patent/CN104915567B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103559382A (zh) * | 2013-10-09 | 2014-02-05 | 济南大学 | 一种水泥分解炉入炉煤粉低位热值在线估计方法 |
| CN104061589A (zh) * | 2014-07-17 | 2014-09-24 | 烟台龙源电力技术股份有限公司 | 适应煤质变化的低氮燃烧控制方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104915567A (zh) | 2015-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104915567B (zh) | 基于热平衡分析的水泥熟料烧成热耗、热效率及燃煤热值实时监测方法 | |
| CN110375289B (zh) | 一种燃煤电站锅炉监测系统及其热效率测试方法 | |
| CN102012968B (zh) | 煤粉锅炉热效率和煤质数据实时监测方法 | |
| CN110866856A (zh) | 一种钢铁企业温室气体排放监测系统及方法 | |
| CN104699937A (zh) | 一种基于烟气测试的锅炉效率自校正计算方法 | |
| CN106323657B (zh) | 一种基于化学平衡的燃煤锅炉能效在线监测方法 | |
| CN110866857A (zh) | 一种水泥厂温室气体排放监测系统及方法 | |
| CN106093336B (zh) | 一种在线确定锅炉燃烧煤质元素成分的方法 | |
| CN103995987A (zh) | 掺烧高炉煤气的煤粉锅炉的热效率测算方法 | |
| CN212061242U (zh) | 一种钢铁企业温室气体排放监测系统 | |
| CN210740365U (zh) | 一种燃煤电站锅炉监测系统 | |
| CN109855892A (zh) | 基于碳平衡的固体不完全燃烧热损失快速检测方法 | |
| CN108982683A (zh) | 一种煤矿乏风通入燃煤锅炉混烧氧化的分析装置及方法 | |
| CN104615895A (zh) | 煤粉与高炉煤气混烧锅炉空气预热器漏风率的测算方法 | |
| CN101306312B (zh) | 一种炭素制品生产过程中的烟气净化处理方法 | |
| CN114091276B (zh) | 生活垃圾焚烧锅炉的锅炉热效率在线计算方法 | |
| CN102690072A (zh) | 一种水泥生产过程不完全燃烧碳化物在线监测控制方法 | |
| CN105823863A (zh) | 一种基于恒温热重分析的煤质在线工业分析测量方法 | |
| CN112495568B (zh) | 一种基于煤的比热变化的磨煤机煤种切换装置及判别方法 | |
| CN109655489A (zh) | 一种燃煤机组标准煤耗的计量方法及装置 | |
| CN109002687B (zh) | 基于钙平衡确定循环流化床锅炉的灰渣份额的方法 | |
| CN112417684B (zh) | 一种燃煤机组煤质成分动态计算方法 | |
| CN116697392A (zh) | 一种煤粉锅炉碳燃尽系数的在线监测方法 | |
| Sergeyev et al. | Researches of air and fuel rate influence on oxygen level in emissions of new type medium power coal boiler | |
| CN208313922U (zh) | 一种超低挥发分碳基燃料变压力燃烧与排放特性试验台系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C53 | Correction of patent for invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yin Jie Inventor after: Liu Jianxiong Inventor after: Niu Honghai Inventor after: Peng Xing Inventor after: Chen Jun Inventor after: Geng Xin Inventor after: Li Bing Inventor before: Yin Jie Inventor before: Niu Honghai Inventor before: Peng Xing Inventor before: Chen Jun Inventor before: Li Bing Inventor before: Geng Xin |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: YIN JIE NIU HONGHAI PENG XING CHEN JUN LI BING GENG XIN TO: YIN JIE LIU JIANXIONG NIU HONGHAI PENG XING CHEN JUN GENG XIN LI BING |
|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171027 Termination date: 20200618 |