本发明涉及一种煤质波动状态在线软测量方法，具体讲是一种基于磨煤机实时运行状态监测数据，运用热平衡原理，在线计算原煤水分的方法。
背景技术：
：目前，对于原煤水分等煤质状态特性指标，主要依靠电厂运行人员每8小时、每天定期抽样、化验原煤得到，该方法分析结果存在两方面不足：首先，由于原煤抽样过程是随机过程，样品的代表性不一定能够保证；其次，由于原煤化验过程每天定期开展，最小周期也需要8小时，无法实时在线分析煤质成分，用每天定期化验值代替8小时平均值，存在一定的偏差。本发明基于磨煤机实时运行状态监测数据，运用热平衡原理，在线计算原煤水分的方法，能够实现煤质状态的实时在线评估，较传统化验结果更具有代表性及可靠性。技术实现要素：本发明的目的是基于磨煤机实时运行状态监测数据，运用热平衡原理，在线计算原煤水分的方法，能够实现磨煤机干燥出力的在线评估，有利于强化运行人员及调度人员调节的准确性及针对性。为了实现上述的目的，本发明是采取以下的技术方案来实现的：一种煤质波动状态在线软测量方法，其特征是，通过在线迭代计算原煤水分Mar，初始时原煤水分Mar＝0，包括以下步骤：(1)计算干燥剂流量：其中Qv代表干燥剂量，Bm代表给煤量，Qs代表密封风量；(2)计算干燥剂物理热：qag1＝cag1×t1×g1，其中t1代表磨煤机进口一次风温，Cag1代表入口干燥剂比热；(3)计算密封风物理热：其中Cs代表密封风比热，ts代表密封风温度；(4)测量磨煤机碾磨热qmac；(5)计算输入总热量：qin＝qag1+qs+qmac；(6)计算蒸发水分消耗的热量qev=ΔM×(2500+cH2O′′×t2-CH2O×trc);]]>其中，t2代表磨煤机出口风粉温度，代表水比热，trc代表进入系统原煤温度，Mpc代表煤粉水分；(7)计算乏气干燥剂带出热量：qag2=[g1+QsBm]×cag2×t2]]>其中Cag2代表出口乏气比热，t2代表磨煤机出口混合物的温度；(8)计算加热燃料消耗的热量：qf=100-Mar100[0.0034×(t2+trc)×0.5+0.8796+cH2O×Mpc100-Mpc]×(t2-trc)]]>(9)计算设备散热损失：q5＝0.02×qin(10)输出总热量：qout＝qev+qag2+qf+q5(11)根据热平衡计算得到蒸发原煤中水分消耗的热量：q′ev＝qin-qag2-qf-q5(12)根据热平衡计算得到单位原煤蒸发掉的水量ΔM′=qev′/(2500+cH2O′′×t2-CH2O×trc)]]>(13)根据热平衡计算得到原煤水分Mar＝(100-Mpc)*ΔM'+Mpc；(14)平衡校验计算若|qin-qout|>0.1，则使用第(13)步计算得到的Mar从第(1)步重新计算，直到满足条件为止。本发明所达到的有益效果：本发明的方法，基于磨煤机实时运行状态监测数据，运用热平衡原理，在线计算原煤水分的方法，能够实现磨煤机干燥出力的在线评估，有利于强化运行人员及调度人员调节的准确性及针对性。具体实施方式下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。一、常量数据表二、水分计算流程初始原煤全水，即原煤水分Mar＝0(1)干燥剂流量(kg/kg)其中Qv代表干燥剂量(磨煤机进口风量)，t/h；Bm代表给煤量，t/h。(2)干燥剂物理热(kj/kg)qag1＝cag1×t1×g1，其中t1代表磨煤机进口一次风温；(3)密封风物理热(kj/kg)qs=QsBmcs×ts]]>(4)磨煤机碾磨热(kj/kg)qmac＝12(5)输入总热量(kj/kg)qin＝qag1+qs+qmac(6)蒸发水分消耗的热量(kj/kg)qev=ΔM×(2500+cH2O′′×t2-CH2O×trc)]]>ΔM=Mar-Mpc100-Mpc]]>其中t2代表磨煤机出口风粉温度(7)乏气干燥剂带出热量(kj/kg)qag2=[g1+QsBm]×cag2×t2]]>(8)加热燃料消耗的热量(kj/kg)qf=100-Mar100[0.0034×(t2+trc)×0.5+0.8796+cH2O×Mpc100-Mpc]×(t2-trc)]]>(9)设备散热损失(kj/kg)q5＝0.02×qin(10)输出总热量(kj/kg)qout＝qev+qag2+qf+q5(11)热平衡计算得到蒸发原煤中水分消耗的热量(kj/kg)q′ev＝qin-qag2-qf-q5(12)热平衡计算得到单位原煤蒸发掉的水量(％)ΔM′=qev′/(2500+cH2O′′×t2-CH2O×trc)]]>(13)热平衡计算得到原煤水分(％)Mar＝(100-Mpc)*ΔM'+Mpc(14)平衡校验计算若|qin-qout|>0.1，则使用第(13)步计算得到的Mar从第(1)步至第(14)步重新计算，直到满足条件即|qin-qout|≤0.1为止。实施例1下面以某1000MW等级机组为例，介绍一下该算法的具体应用过程：该机组实测数据如下：过热蒸汽流量Dgrt/h2887过热器出口蒸汽压力Pgr"Mpa26.07过热器出口蒸汽温度Tgr"℃605过热器出口蒸汽焓hgr"kJ/kg3498.77再热蒸汽流量Dzrt/h2339再热器出口蒸汽压力Pzr"Mpa4.68再热器出口蒸汽温度Tzr"℃603再热器出口蒸汽焓hzr"kJ/kg3676.35再热器进口蒸汽压力Pzr′Mpa4.93再热器进口蒸汽温度Tzr′℃351再热器进口蒸汽焓hzr′kJ/kg3073.67锅炉给水压力PgsMpa29.51锅炉给水温度Tgs℃300锅炉给水焓hgskJ/kg1328.87总有效利用热QkJ/kg7674169.82燃料消耗量Bt/h293.48磨煤机台数台6单台磨煤机出力Bmt/h64.56型号RP983基本出力B0t/h65.3电机输入功率PmaxkW468单位电耗kWh/t7.16运用该算法计算结果为：以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。当前第1页1 2 3