发明名称 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法
摘要 本发明针对现有炉膛出口烟温软测量方法的缺陷,公开了一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法,此方法在烟温软测量中采用水冷壁热有效系数反映炉膛实时结渣情况,并在计算考虑了其实时可变性;计算输出部分主要包括炉膛和部分对流受热面联合传热计算,以及水冷壁热有效系数和炉膛出口烟气温度的输出,在不增加锅炉测点的基础上,采集入炉煤质参数、结构参数和锅炉实时运行参数,通过联合炉膛及对流受热面的传热平衡计算,并考虑炉膛内实时结渣性,从而建立炉膛出口烟温实时推算软测量模型,能够输出炉膛出口烟温的实时值,并同时提供炉膛内部实时结渣情况,可以作为运行人员进行炉膛吹灰以及高温受热面减温水操作的直观数据参考。
申请公布号 CN105276563A 申请公布日期 2016.01.27
申请号 CN201510764823.6 申请日期 2015.11.11
申请人 东南大学 发明人 黄亚继;徐力刚;陈波
分类号 F22B37/38(2006.01)I 主分类号 F22B37/38(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 柏尚春
主权项 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法,其特征在于此方法在烟温软测量中采用水冷壁热有效系数反映炉膛实时结渣情况,并在计算考虑了其实时可变性;计算输出部分主要包括炉膛和部分对流受热面联合传热计算,以及水冷壁热有效系数和炉膛出口烟气温度的输出,具体有以下步骤:步骤1:计算燃料带入炉内的有效热量、理论燃烧温度,辐射吸收减弱系数和火焰综合黑度,为炉膛辐射传热计算作准备:(1)理论冷风焓:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mrow><mi>c</mi><mi>k</mi></mrow><mn>0</mn></msubsup><mo>=</mo><msub><mi>g</mi><mn>1</mn></msub><msubsup><mi>I</mi><mrow><mi>k</mi><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><msubsup><mi>I</mi><mrow><mi>k</mi><mn>2</mn></mrow><mo>&prime;</mo></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000011.GIF" wi="382" he="79" /></maths>理论热风焓:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>I</mi><mrow><mi>c</mi><mi>k</mi></mrow><mn>0</mn></msubsup><mo>=</mo><msub><mi>g</mi><mn>1</mn></msub><msubsup><mi>I</mi><mrow><mi>k</mi><mn>1</mn></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><msubsup><mi>I</mi><mrow><mi>k</mi><mn>2</mn></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000012.GIF" wi="356" he="79" /></maths>(此处的空气焓均由测点空气温度根据空气焓温表查取);进入炉内空气热量:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>&alpha;</mi><mi>f</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>-</mo><msub><mi>&Delta;&alpha;</mi><mi>f</mi></msub><mo>-</mo><msub><mi>&Delta;&alpha;</mi><mrow><mi>p</mi><mi>c</mi><mi>s</mi></mrow></msub><mo>)</mo></mrow><msubsup><mi>I</mi><mrow><mi>h</mi><mi>k</mi></mrow><mn>0</mn></msubsup><mo>+</mo><mrow><mo>(</mo><msub><mi>&Delta;&alpha;</mi><mi>f</mi></msub><mo>+</mo><msub><mi>&Delta;&alpha;</mi><mrow><mi>p</mi><mi>c</mi><mi>s</mi></mrow></msub><mo>)</mo></mrow><msubsup><mi>I</mi><mrow><mi>c</mi><mi>k</mi></mrow><mn>0</mn></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000013.GIF" wi="934" he="84" /></maths>燃料带入炉内有效热量:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>f</mi><mrow><mi>e</mi><mi>f</mi></mrow></msubsup><mo>=</mo><msub><mi>Q</mi><mi>r</mi></msub><mfrac><mrow><mn>100</mn><mo>-</mo><msub><mi>q</mi><mn>3</mn></msub><mo>-</mo><msub><mi>q</mi><mn>4</mn></msub><mo>-</mo><msub><mi>q</mi><mn>6</mn></msub></mrow><mrow><mn>100</mn><mo>-</mo><msub><mi>q</mi><mn>4</mn></msub></mrow></mfrac><mo>+</mo><msub><mi>Q</mi><mi>k</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000014.GIF" wi="654" he="141" /></maths>(2)炉内有效热量<img file="FDA0000843776270000015.GIF" wi="77" he="71" />即为理论燃烧温度T<sub>th</sub>对应焓值,在获得<img file="FDA0000843776270000016.GIF" wi="74" he="71" />后通过烟气焓温表利用插值法查取T<sub>th</sub>;(3)烟气中灰分颗粒的质量浓度:<img file="FDA0000843776270000017.GIF" wi="326" he="148" />烟气中焦炭颗粒的容积浓度:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mrow><mi>c</mi><mi>o</mi><mi>k</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>5.5</mn><msub><mi>C</mi><mrow><mi>a</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mrow><mn>10</mn><mo>+</mo><msub><mi>q</mi><mn>4</mn></msub></mrow><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mn>100</mn><mo>+</mo><msub><mi>V</mi><mrow><mi>d</mi><mi>a</mi><mi>f</mi></mrow></msub><mo>)</mo><msub><mi>V</mi><mi>y</mi></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>t</mi></msub><mo>-</mo><msub><mi>h</mi><mrow><mi>u</mi><mi>n</mi></mrow></msub></mrow><msub><mi>h</mi><mi>f</mi></msub></mfrac></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000018.GIF" wi="749" he="165" /></maths>三原子气体减弱系数:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>k</mi><mi>g</mi></msub><mi>r</mi><mo>=</mo><mrow><mo>&lsqb;</mo><mrow><mfrac><mrow><mn>0.78</mn><mo>+</mo><mn>1.6</mn><msub><mi>r</mi><mrow><msub><mi>H</mi><mn>2</mn></msub><mi>O</mi></mrow></msub></mrow><msqrt><mrow><mi>r</mi><mi>S</mi></mrow></msqrt></mfrac><mo>-</mo><mn>0.1</mn></mrow><mo>&rsqb;</mo></mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mn>0.37</mn><mfrac><msubsup><mi>T</mi><mi>f</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mn>1000</mn></mfrac></mrow><mo>)</mo></mrow><mi>r</mi><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000019.GIF" wi="910" he="158" /></maths>灰分减弱系数:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>k</mi><mrow><mi>a</mi><mi>s</mi><mi>h</mi></mrow></msub><msub><mi>&mu;</mi><mrow><mi>a</mi><mi>s</mi><mi>h</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>5330</mn><msub><mi>&mu;</mi><mrow><mi>a</mi><mi>s</mi><mi>h</mi><mo>,</mo><mi>m</mi></mrow></msub></mrow><msup><mrow><mo>(</mo><msubsup><mi>T</mi><mi>f</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><msub><mi>d</mi><mrow><mi>a</mi><mi>s</mi><mi>h</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></mfrac><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><mfrac><mn>0.65</mn><mrow><mn>1</mn><mo>+</mo><mn>0.0177</mn><mo>/</mo><msup><mrow><mo>(</mo><msub><mi>&mu;</mi><mrow><mi>a</mi><mi>s</mi><mi>h</mi><mo>,</mo><mi>m</mi></mrow></msub><mi>S</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA00008437762700000110.GIF" wi="966" he="165" /></maths>焦炭颗粒减弱系数:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>k</mi><mrow><mi>c</mi><mi>o</mi><mi>k</mi></mrow></msub><msub><mi>&mu;</mi><mrow><mi>c</mi><mi>o</mi><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>10</mn><msub><mi>&mu;</mi><mrow><mi>c</mi><mi>o</mi><mi>k</mi><mo>,</mo><mi>v</mi></mrow></msub></mrow><msup><mrow><mo>(</mo><msubsup><mi>T</mi><mi>f</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><msub><mi>d</mi><mrow><mi>c</mi><mi>o</mi><mi>k</mi></mrow></msub><mo>)</mo></mrow><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></mfrac><mo>;</mo></mrow>]]></math><img file="FDA00008437762700000111.GIF" wi="438" he="143" /></maths>辐射吸收减弱系数:k<sub>a</sub>=k<sub>g</sub>r+k<sub>ash</sub>μ<sub>ash</sub>+k<sub>cok</sub>μ<sub>cok</sub>;(4)炉膛实际火焰黑度:<img file="FDA00008437762700000115.GIF" wi="278" he="77" />炉膛火焰综合黑度:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>&epsiv;</mi><mi>f</mi></msub><mrow><mn>0.32</mn><msub><mi>k</mi><mi>a</mi></msub><msub><mi>R&epsiv;</mi><mi>f</mi></msub><mo>+</mo><mn>1</mn></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA00008437762700000112.GIF" wi="429" he="127" /></maths>其中,<img file="FDA00008437762700000113.GIF" wi="146" he="71" />分别为理论的热空气焓和冷空气焓,kJ/kg;I<sub>k1</sub>′,I<sub>k1</sub>″分别为一次风进出口空气焓,kJ/kg;I<sub>k2</sub>′,I<sub>k2</sub>″分别为二次风进出口空气焓,KJ/kg;g<sub>1</sub>,g<sub>2</sub>分别为一、二次风流量占总空气流量的份额;Q<sub>k</sub>为随单位质量燃料带入炉内的空气含漏风的热量,kJ/kg。Q<sub>r</sub>为单位质量燃料带入炉内的热量,通常等于燃料收到基低位发热量,kJ/kg;q<sub>3</sub>为化学未完全燃烧热损失,%;q<sub>4</sub>为机械未完全燃烧热损失,%;q<sub>6</sub>为其他热损失,%;<img file="FDA00008437762700000114.GIF" wi="74" he="70" />为单位质量燃料带入炉内的有效热,KJ/kg;T<sub>th</sub>为理论燃烧温度,K;k<sub>g</sub>r为三原子气体减弱系数,m<sup>‑1</sup>;k<sub>ash</sub>μ<sub>ash</sub>为灰分颗粒减弱系数,m<sup>‑1</sup>;k<sub>cok</sub>μ<sub>cok</sub>为焦炭颗粒减弱系数,m<sup>‑1</sup>;<img file="FDA00008437762700000210.GIF" wi="89" he="61" />为水蒸汽占烟气容积份额;r为水蒸汽和二氧化物之和占烟气容积份额;μ<sub>ash,m</sub>为烟气中灰分颗粒的质量浓度,kg/kg;μ<sub>cok,v</sub>为烟气中焦炭颗粒的容积浓度,g/Nm<sup>3</sup>;d<sub>ash</sub>为灰分颗粒的平均粒径,μm根据煤种查取;d<sub>cok</sub>为焦炭颗粒的平均粒径,;C<sub>ar</sub>为收到基碳元素,%;A<sub>ar</sub>为收到基灰分,%;α<sub>fa</sub>为飞灰系数;G<sub>y</sub>为单位燃料燃烧产生烟气质量,kg/kg;V<sub>y</sub>为单位燃料燃烧产生烟气体积,m<sup>3</sup>/kg;q<sub>4</sub>为机械未完全燃烧损失,%;V<sub>daf</sub>为干燥无灰基挥发份,%;h<sub>t</sub>为最上排燃烧器布置高度,m;h<sub>un</sub>为最下排燃烧器布置高度,m;ε<sub>f</sub>为炉膛实际火焰黑度;ε<sub>syn</sub>为考虑了火焰辐射强度因介质吸收而减弱的火焰综合黑度;S为炉内辐射层有效厚度,m;R为与炉膛截面等面积圆形的半径,m;步骤2:假设炉膛出口烟温T<sub>f</sub>″进行炉膛传热计算,求出水冷壁热有效系数ψ:(1)炉膛假想黑度:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msubsup><mi>&epsiv;</mi><mi>f</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msubsup><mo>=</mo><mfrac><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub><mrow><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub></mrow><mo>)</mo></mrow><mi>&psi;</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000021.GIF" wi="485" he="140" /></maths>(2)炉膛辐射传热关系式:<img file="FDA0000843776270000022.GIF" wi="726" he="167" />(3)联立本步骤中的关系式(1)和(2),则水冷壁热有效系数:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>&psi;</mi><mo>=</mo><mfrac><mrow><msub><mi>Z&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub></mrow><mrow><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&epsiv;</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mi>Z</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000023.GIF" wi="446" he="151" /></maths>式中<img file="FDA0000843776270000024.GIF" wi="564" he="158" />为求解方便所设参数,无实际意义;其中,<img file="FDA0000843776270000025.GIF" wi="72" he="63" />为炉膛黑度即假想黑度;ψ为水冷壁热有效系数;x<sub>m</sub>为炉膛火焰最高温度位置的相对高度,近似取等于燃烧器布置的相对高度;σ<sub>0</sub>为玻尔兹曼常数,通常取5.67×10<sup>‑11</sup>kW/(m<sup>2</sup>·K<sup>4</sup>);B<sub>j</sub>为计算燃烧量,kg/s;<img file="FDA0000843776270000026.GIF" wi="44" he="53" />为保热系数;H<sub>f</sub>为水冷壁的吸热表面积,m<sup>2</sup>;<img file="FDA0000843776270000027.GIF" wi="108" he="86" />为炉内烟气在理论燃烧温度至炉膛出口温度区间内的平均热容,KJ/(kg·K);步骤3:计算炉膛出口直接辐射量,根据烟温测点列出部分涉及计算的对流受热面热平衡方程,进行炉膛出口烟温的逆烟气流程推算:(1)炉膛出口直接辐射:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>f</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&psi;</mi><mrow><mi>p</mi><mo>,</mo><mi>f</mi></mrow></msub><msubsup><mi>&epsiv;</mi><mi>f</mi><mrow><mi>s</mi><mi>y</mi><mi>n</mi></mrow></msubsup><msub><mi>&sigma;</mi><mn>0</mn></msub><msup><mrow><mo>(</mo><msubsup><mi>T</mi><mi>f</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>)</mo></mrow><mn>4</mn></msup><msub><mi>F</mi><mi>f</mi></msub></mrow><msub><mi>B</mi><mi>j</mi></msub></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000843776270000028.GIF" wi="485" he="157" /></maths>(2)通常电站锅炉烟道中低温对流受热面的烟温是布置测点的,但是靠近炉膛出口的部分高温对流受热面未布置烟温测点,从最靠近炉膛出口的可测得烟温处,逆烟气流动方向列出涉及计算的高温对流受热面的热平衡方程;工质侧:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><mfrac><msub><mi>D</mi><mi>i</mi></msub><msub><mi>B</mi><mi>j</mi></msub></mfrac><mrow><mo>(</mo><mrow><msubsup><mi>h</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>-</mo><msubsup><mi>h</mi><mi>i</mi><mo>&prime;</mo></msubsup></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>Q</mi><mrow><mi>d</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>x</mi><mi>i</mi></msub><msub><mi>Q</mi><mi>f</mi></msub></mrow>]]></math><img file="FDA0000843776270000029.GIF" wi="468" he="142" /></maths>烟气侧:<img file="FDA0000843776270000031.GIF" wi="710" he="71" />上式中ΔQ<sub>i</sub>为变负荷变化速率过快下对流受热面总蓄热量:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msub><mi>&Delta;Q</mi><mi>i</mi></msub><mo>=</mo><msub><mi>&Delta;Q</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>&Delta;Q</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>M</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub><msub><mi>c</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>M</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><msub><mi>c</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mfrac><mrow><mo>&part;</mo><msub><mi>t</mi><mrow><mi>g</mi><mi>i</mi></mrow></msub></mrow><mrow><mo>&part;</mo><mi>&tau;</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000843776270000032.GIF" wi="845" he="127" /></maths>此蓄热量只在负荷变化大于2%/min~3%/min时考虑,稳定负荷或负荷变化较慢时总蓄热量为零;通过联立部分未知烟温的对流受热面的工质及烟气侧方程,则可推算出炉膛出口烟气焓I<sub>f</sub>″,通过烟气焓温表利用插值法求出另一个炉膛出口烟气温度(T<sub>f</sub>″)<sub>2</sub>;其中,下标i代表涉及计算的对流受热面;Q<sub>f</sub>为炉膛出口直接辐射量,KJ/kg;ψ<sub>p,f</sub>为炉膛出口对流受热面获取炉膛直接辐射热有效系数,其实就是炉膛出口处的局部热有效系数,通常用炉膛水冷壁整体热有效系数ψ乘以一个0.8~0.9系数;F<sub>f</sub>为炉膛出口烟窗截面积,m<sup>2</sup>;ΔQ<sub>i</sub>为对流受热面总蓄热量,KJ/kg;ΔQ<sub>gi</sub>、ΔQ<sub>ji</sub>分别为对流受热面工质和金属侧蓄热量,KJ/kg;M<sub>gi</sub>、M<sub>ji</sub>分别为对流受热面工质和金属质量,Kg;c<sub>gi</sub>、c<sub>ji</sub>分别为对流受热面工质和金属比热,kJ/kg·℃;t<sub>gi</sub>为对流受热面工质温度,℃;x<sub>i</sub>为对流受热面接收炉膛出口直接辐射的份额,可根据锅炉设计书得知,纯对流受热面x<sub>i</sub>取0;D<sub>i</sub>为对流受热面工质流量,kg/s;B<sub>j</sub>为计算燃烧量,kg/s;h<sub>i</sub>′、h<sub>i</sub>″吩别为对流受热面进出口工质焓,kJ/kg;Q<sub>di</sub>为对流受热面对流吸热量,KJ/kg;<img file="FDA0000843776270000033.GIF" wi="43" he="60" />为保热系数;I<sub>f</sub>″为炉膛出口烟气焓,KJ/kg;I<sub>i</sub>为最靠近炉膛出口的可测得烟气焓,kJ/kg;Δα<sub>i</sub>为受热面漏风系数;I<sub>k</sub>为漏入空气焓,kJ/kg;步骤4:根据校核式|T<sub>f</sub>″‑(T<sub>f</sub>″)<sub>2</sub>|<1校核假设炉膛出口烟温T<sub>f</sub>″,如满足校核式则输出炉膛出口烟温T<sub>f</sub>″及水冷壁热有效系数ψ;如不满足校核式则假设错误重复步骤2至步骤4。
地址 210096 江苏省南京市四牌楼2号