发明名称 一种电厂SCR脱硝系统催化剂活性在线监测方法
摘要 本发明属于电厂选择性催化还原法进行烟气尾部脱硝技术领域,特别涉及一种电厂SCR脱硝系统催化剂活性在线监测方法。本发明方法主要包括:一、获取实验室进行相同规格新鲜催化剂的活性检测,得到新鲜催化剂的活性a<sub>0</sub>=1,其能达到的最佳脱硝效率为η<sub>0</sub>;二、在机组实际运行中,对SCR脱硝系统进行现场性能试验,获取SCR实际运行数据及性能参数;三、利用本发明提出的SCR催化剂活性修正方法,计算当前催化剂的真实活性。本发明的催化剂活性在线监测方法旨在提供一种SCR脱硝系统催化剂活性在线修正计算方法,修正因实验室取样代表性问题和反应条件过于理想化而造成计算得到的活性偏离真实值,以提高催化剂活性检测的准确性和实时性。
申请公布号 CN103592407A 申请公布日期 2014.02.19
申请号 CN201310549775.X 申请日期 2013.11.07
申请人 华北电力大学 发明人 付忠广;齐敏芳;曹宏芳;陈颖;张越
分类号 G01N31/10(2006.01)I 主分类号 G01N31/10(2006.01)I
代理机构 北京麟保德和知识产权代理事务所(普通合伙) 11428 代理人 韩建功
主权项 1.一种电厂SCR脱硝系统催化剂活性在线监测方法,其特征在于:具体步骤如下:(1)分析确定SCR脱硝系统催化剂结构参数和SCR脱硝系统运行参数;(2)根据所选择的催化剂,作出实验室催化剂活性检测试验报告,新鲜催化剂的活性即相对活性a<sub>0</sub>=1,对应的新鲜催化剂可达到的最高脱硝效率为η<sub>0</sub>,脱硝效率:<img file="FDA0000409665790000011.GIF" wi="640" he="84" />催化剂面速度:AV=Q/(V*Ap),催化剂活性系数:K=-AV*ln(1-η),催化剂活性:a=K/K<sub>0</sub>;其中:η为脱硝效率;<img file="FDA0000409665790000012.GIF" wi="138" he="80" />为SCR入口NO<sub>x</sub>浓度,单位为mg/Nm<sup>3</sup>;<img file="FDA0000409665790000013.GIF" wi="157" he="83" />为SCR出口NO<sub>x</sub>浓度,单位为mg/Nm<sup>3</sup>;AV为催化剂的面速度,单位为Nm/h;Q为烟气流量,单位为Nm<sup>3</sup>/h,V为催化剂的总体积,单位为m<sup>3</sup>;Ap为催化剂的比表面积,单位为m<sup>2</sup>/m<sup>3</sup>;K为实际催化剂活性系数,单位为Nm/h;K<sub>0</sub>为新鲜催化剂的活性系数,单位为Nm/h;a为催化剂活性;(3)引入活性修正因子k对催化剂活性进行修正,以得到催化剂活性的真实值,其中k≥0;(4)首次应用该方法时,对电厂SCR脱硝系统催化剂活性起始值的确定方法如下:根据电厂SCR脱硝系统投运时间及催化剂活性取值范围,选用如下方式之一进行催化剂活性起始值取值:①催化剂活性起始值取1;②催化剂活性起始值取[0,1)之间的值;③参考最近一次催化剂活性检测报告,以检测报告中的催化剂活性作为该方法中催化剂活性的起始值;(5)获取t时刻的理论催化剂活性a<sub>t</sub>:t时刻的理论催化剂活性计算:<img file="FDA0000409665790000021.GIF" wi="558" he="142" />在此忽略催化剂比表面积的变化引起的影响,整理得ln(1-η<sub>t</sub>)=a<sub>t</sub>*ln(1-η<sub>0</sub>);k<sub>0</sub>=ln(1-η<sub>0</sub>)为常数,且为负值,则有ln(1-η<sub>t</sub>)=a<sub>t</sub>*k<sub>0</sub>;引入活性修正因子k对a<sub>t</sub>进行修正,a<sub>t,best</sub>为t时刻催化剂的最佳活性,则t时刻的理论催化剂活性为a<sub>t</sub>=k[j]*a<sub>t,best</sub>;设置k的初始值k[0]=1.0;假设t-1时刻的催化剂真实活性为a<sub>t-1</sub>,则a<sub>t,best</sub>=a<sub>t-1</sub>,即将上一时刻计算得到的催化剂活性真实值作为当前时刻催化剂的最佳活性;整理得ln(1-η<sub>t</sub>)=a<sub>t,best</sub>*k[j]*k<sub>0</sub>=a<sub>t</sub>*k<sub>0</sub>;t时刻的理论催化剂脱硝效率为η<sub>t</sub>=1-exp(a<sub>t</sub>*k<sub>0</sub>);理论SCR出口NO<sub>x</sub>浓度:<maths num="0001"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><msub><mi>NO</mi><mi>x</mi></msub><mo>,</mo><mi>cal</mi></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mi>N</mi><msub><mi>O</mi><mi>x</mi></msub><mo>,</mo><mi>in</mi></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&eta;</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mrow><mi>N</mi><msub><mi>O</mi><mi>x</mi></msub><mo>,</mo><mi>in</mi></mrow></msub><mo>*</mo><mi>exp</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>t</mi></msub><mo>*</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,K<sub>t</sub>为t时刻实际催化剂活性系数,单位为Nm/h;K<sub>0</sub>为新鲜催化剂的活性系数,单位为Nm/h;AV<sub>t</sub>为t时刻的催化剂面速度单位为Nm/h;AV<sub>0</sub>为新鲜催化剂面速度单位为Nm/h;η<sub>t</sub>为t时刻的理论催化剂脱硝效率;a<sub>t</sub>为t时刻的理论催化剂活性;a<sub>t,best</sub>为t时刻催化剂的最佳活性;j为活性修正因子迭代次数;a<sub>t-1</sub>为t-1时刻的催化剂真实活性;<img file="FDA0000409665790000023.GIF" wi="150" he="81" />理论SCR出口NO<sub>x</sub>浓度,单位为mg/Nm<sup>3</sup>;(6)计算SCR反应器实际出口NO<sub>x</sub>浓度<img file="FDA0000409665790000024.GIF" wi="159" he="79" />与根据催化剂活性计算得到的理论SCR出口NO<sub>x</sub>浓度<img file="FDA0000409665790000027.GIF" wi="137" he="67" />偏差的绝对误差相对值<img file="FDA0000409665790000025.GIF" wi="479" he="177" />若ε小于设定的误差因子ε<sub>0</sub>,则保存并输出当前的活性修正因子k[n]以及催化剂真实活性<img file="FDA0000409665790000028.GIF" wi="46" he="67" />,迭代计算结束;若ε大于设定的误差因子ε<sub>0</sub>,则继续进行步骤(7);其中,<img file="FDA0000409665790000026.GIF" wi="147" he="80" />的单位为mg/Nm<sup>3</sup>,n为终止迭代次数;(7)调整活性修正因子:在第j次迭代结束后,ε大于设定的误差因子ε<sub>0</sub>,根据SCR反应器实际出口NO<sub>x</sub>浓度与SCR出口NO<sub>x</sub>浓度理论值的偏差大小,调整活性修正因子;若实际值大于理论值,则减小活性修正因子;若实际值小于理论值,则增大活性修正因子;(8)重复(5)(6)(7)进行迭代计算,直至SCR反应器实际出口NO<sub>x</sub>浓度<img file="FDA0000409665790000031.GIF" wi="148" he="77" />与根据催化剂活性计算得到的理论SCR出口NO<sub>x</sub>浓度<img file="FDA0000409665790000032.GIF" wi="152" he="79" />偏差的绝对误差相对值ε小于设定的误差因子ε<sub>0</sub>。
地址 102206 北京市昌平区回龙观朱辛庄2号