发明名称 一种SCR脱硝催化剂最佳更换时间的实时计算方法
摘要 一种SCR脱硝催化剂最佳更换时间的实时计算方法,确定能够反映催化剂状态的SCR脱硝系统参数,依据边界条件获取各负荷工况下参数的实时运行数据,并进行基于拉依达准则(3σ准则)的数据滤波处理,根据数据滤波处理后的有效数据绘制SCR脱硝系统参数时间变化曲线,以此为依据计算脱硝催化剂的理论更换时间,并从经济性的角度实时计算并预测脱硝催化剂的最佳更换时间。
申请公布号 CN105844083A 申请公布日期 2016.08.10
申请号 CN201610153795.9 申请日期 2016.03.17
申请人 西安西热电站信息技术有限公司 发明人 王大鹏;王智微;马鹏波;杜保华
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 刘国智
主权项 一种SCR脱硝催化剂最佳更换时间的实时计算方法,其特征在于:包括以下步骤:1)确定反映催化剂状态的SCR脱硝系统在线参数即脱硝效率E<sub>SCR</sub>、喷氨流量S<sub>NH3</sub>、出口氨浓度或氨逃逸率PPM、NO<sub>X</sub>排放浓度C<sub>NOx‑Out</sub>,确定SCR反应器温度T、烟气流量F<sub>Gas</sub>和SCR系统NO<sub>X</sub>入口浓度C<sub>NOx‑In</sub>为SCR系统稳态运行条件参数;2)确定方案的边界条件机组运行状态良好,锅炉燃烧方式、燃烧温度、过量空气系数、烟气在炉内停留时间稳定,且入炉煤质稳定;3)获取各负荷工况下参数的实时运行数据并进行数据过滤处理3‑1)获取机组设计负荷P<sub>i</sub>分别处于低负荷、中负荷和高负荷各要求负荷工况下机组稳态运行时四个SCR脱硝系统参数E<sub>SCR</sub>,S<sub>NH3</sub>,PPM,C<sub>NOx‑Out</sub>的实时运行数据;其中低负荷、中负荷和高负荷分别定义为50%P<sub>e</sub>、75%P<sub>e</sub>和100%P<sub>e</sub>;3‑2)对3‑1)获取的实时运行数据进行基于拉依达准则即3σ准则的数据滤波处理并剔除奇异数据;4)绘制脱硝系统参数时间变化曲线以12小时为时间单位,分别计算过去12小时内所有参数经过数据滤波处理得到的n‑m个有效数据的算术平均值,并用这些算术平均值分别绘制低负荷、中负荷和高负荷工况下四个SCR脱硝系统参数随时间的变化曲线<img file="FDA0000944035210000011.GIF" wi="893" he="87" />5)通过对变化曲线进行数据拟合得到脱硝催化剂的理论更换时间;5‑1)脱硝效率<img file="FDA0000944035210000012.GIF" wi="170" he="87" />降低到设计要求最低脱硝效率<img file="FDA0000944035210000013.GIF" wi="114" he="86" />时,脱硝催化剂达到理论更换时间<img file="FDA0000944035210000014.GIF" wi="165" he="95" />如式(4)所示;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>E</mi><mrow><mi>S</mi><mi>C</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>t</mi><mo>=</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><msub><mi>E</mi><mrow><mi>S</mi><mi>C</mi><mi>R</mi></mrow></msub></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup></mrow></msub><mo>=</mo><msubsup><mi>E</mi><mrow><mi>S</mi><mi>C</mi><mi>R</mi><mn>0</mn></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000021.GIF" wi="1510" he="156" /></maths>5‑2)NO<sub>X</sub>排放浓度<img file="FDA0000944035210000022.GIF" wi="211" he="87" />升高至排放要求极限值<img file="FDA0000944035210000023.GIF" wi="163" he="87" />时,脱硝催化剂达到理论更换时间<img file="FDA0000944035210000024.GIF" wi="217" he="95" />如式(5)所示;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>C</mi><mrow><mi>N</mi><mi>O</mi><mi>x</mi><mo>-</mo><mi>O</mi><mi>u</mi><mi>t</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>t</mi><mo>=</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><msub><mi>C</mi><mrow><mi>N</mi><mi>O</mi><mi>x</mi><mo>-</mo><mi>O</mi><mi>u</mi><mi>t</mi></mrow></msub></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup></mrow></msub><mo>=</mo><msubsup><mi>C</mi><mrow><mi>N</mi><mi>O</mi><mi>x</mi><mo>-</mo><mi>O</mi><mi>u</mi><mi>t</mi><mn>0</mn></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000025.GIF" wi="1510" he="141" /></maths>5‑3)出口氨浓度或氨逃逸率PPM<sup>Pi</sup>(t)升高至排放要求极限值<img file="FDA0000944035210000026.GIF" wi="179" he="86" />时,脱硝催化剂达到理论更换时间<img file="FDA0000944035210000027.GIF" wi="165" he="87" />如式(6)所示;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>PPM</mi><mrow><mi>P</mi><mi>i</mi></mrow></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mo>|</mo><mrow><mi>t</mi><mo>=</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><mi>P</mi><mi>P</mi><mi>M</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup></mrow></msub><mo>=</mo><msubsup><mi>PPM</mi><mn>0</mn><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000028.GIF" wi="861" he="133" /></maths>在Pi负荷工况下,<img file="FDA0000944035210000029.GIF" wi="453" he="87" />和PPM<sup>Pi</sup>(t)三条曲线各自出现转折点,分别对应式(4)、式(5)和式(6)分别推算的理论更换时间<img file="FDA00009440352100000210.GIF" wi="432" he="87" />和<img file="FDA00009440352100000211.GIF" wi="197" he="79" />通常这三个时间不尽相同,因此,定义Pi负荷工况下,脱硝催化剂的理论更换时间<img file="FDA00009440352100000212.GIF" wi="93" he="79" />如式(7)所示:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>=</mo><mi>M</mi><mi>i</mi><mi>n</mi><mo>{</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><msub><mi>E</mi><mrow><mi>S</mi><mi>C</mi><mi>R</mi></mrow></msub></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>,</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><msub><mi>C</mi><mrow><mi>N</mi><mi>O</mi><mi>x</mi><mo>-</mo><mi>O</mi><mi>u</mi><mi>t</mi></mrow></msub></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>,</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi><mo>-</mo><mi>P</mi><mi>P</mi><mi>M</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009440352100000213.GIF" wi="1078" he="95" /></maths>6)计算脱硝催化剂的最佳更换时间在Pi负荷工况下,如果达到理论更换时间<img file="FDA00009440352100000214.GIF" wi="93" he="77" />而不更换脱硝催化剂,通过增加的喷氨消耗量与液氨成本单价计算增加的喷氨成本,并计入相关脱硝设备增加的电量,计算得到增加的运行消耗成本,如式(8)所示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>C</mi><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>p</mi><mrow><mi>N</mi><mi>H</mi><mn>3</mn></mrow></msub><mo>&CenterDot;</mo><msubsup><mo>&Integral;</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><msup><mi>t</mi><mrow><mi>P</mi><mi>i</mi></mrow></msup></msubsup><mrow><mo>(</mo><msubsup><mi>S</mi><mrow><mi>N</mi><mi>H</mi><mn>3</mn></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msubsup><mi>S</mi><mrow><mi>N</mi><mi>H</mi><mn>3</mn></mrow><mrow><mo>&prime;</mo><mi>P</mi><mi>i</mi></mrow></msubsup><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mi>d</mi><mi>t</mi><mo>+</mo><msub><mi>p</mi><mrow><mi>K</mi><mi>W</mi><mi>H</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mo>&Integral;</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><msup><mi>t</mi><mrow><mi>P</mi><mi>i</mi></mrow></msup></msubsup><mi>K</mi><mi>W</mi><mi>H</mi><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00009440352100000215.GIF" wi="1605" he="135" /></maths>式中,p<sub>NH3</sub>为液氨成本单价,更换脱硝催化剂后的喷氨流量<img file="FDA00009440352100000216.GIF" wi="206" he="86" />p<sub>KWH</sub>为上网电价,KWH为相关脱硝设备每12小时增加平均电耗;脱硝催化剂更换成本R<sub>C</sub>(t)由催化剂购置成本现值B<sub>C</sub>、安装费用现值I<sub>C</sub>及废催化剂处理费现值D<sub>C</sub>通过复利法计算得到,如式(9)所示,其中y为计息年化利率,t为计息时间,12小时为一个时间单位;<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>R</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>B</mi><mi>C</mi></msub><mo>+</mo><msub><mi>I</mi><mi>C</mi></msub><mo>+</mo><msub><mi>D</mi><mi>C</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>y</mi><mo>)</mo></mrow><mfrac><mrow><mn>12</mn><mi>t</mi></mrow><mn>365</mn></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000031.GIF" wi="966" he="119" /></maths>当额外的运行消耗成本<img file="FDA0000944035210000032.GIF" wi="150" he="79" />与脱硝催化剂更换成本R<sub>C</sub>(t)在延长期内的经济收益相等时,如式(10)所示,可实时计算并预测不同Pi负荷工况下脱硝催化剂的最佳更换时间<img file="FDA0000944035210000033.GIF" wi="102" he="79" /><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>C</mi><mi>C</mi><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>y</mi><mo>&CenterDot;</mo><mfrac><mrow><mo>(</mo><msubsup><mi>t</mi><mrow><mi>O</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>-</mo><msubsup><mi>t</mi><mrow><mi>T</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>)</mo></mrow><mn>365</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000034.GIF" wi="910" he="151" /></maths>或将不同Pi负荷工况对应的额外增加的运行消耗成本<img file="FDA0000944035210000035.GIF" wi="149" he="78" />与脱硝催化剂更换成本R<sub>C</sub>(t)绘制与时间相关的曲线,不同Pi负荷工况对应不同最佳更换时间<img file="FDA0000944035210000036.GIF" wi="100" he="78" />其他边界条件一致的前提下,负荷Pi越高,时间<img file="FDA0000944035210000037.GIF" wi="74" he="79" />越短;根据式式(10)计算的脱销催化剂在Pi低负荷、中负荷和高负荷工况对应不同最佳更换时间<img file="FDA0000944035210000038.GIF" wi="100" he="78" />通过不同Pi负荷工况在机组运行过程中所占比例权重α<sup>Pi</sup>,计算脱硝催化剂的最佳更换时间t<sub>OR</sub>,如式(11)所示<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>t</mi><mrow><mi>O</mi><mi>R</mi></mrow></msub><mo>=</mo><munder><mo>&Sigma;</mo><mi>i</mi></munder><msup><mi>&alpha;</mi><mrow><mi>P</mi><mi>i</mi></mrow></msup><mo>&CenterDot;</mo><msubsup><mi>t</mi><mrow><mi>O</mi><mi>R</mi></mrow><mrow><mi>P</mi><mi>i</mi></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000944035210000039.GIF" wi="622" he="119" /></maths>式中,α<sup>Pi</sup>的取值,根据历史负荷情况考虑,也可根据年度或月度调度计划考虑。
地址 710054 陕西省西安市雁翔路99号博源科技广场A座