CN103926833B - 一种选择性催化还原脱硝装置入口参数确定方法 - Google Patents
一种选择性催化还原脱硝装置入口参数确定方法 Download PDFInfo
- Publication number
- CN103926833B CN103926833B CN201410164248.1A CN201410164248A CN103926833B CN 103926833 B CN103926833 B CN 103926833B CN 201410164248 A CN201410164248 A CN 201410164248A CN 103926833 B CN103926833 B CN 103926833B
- Authority
- CN
- China
- Prior art keywords
- economizer exit
- denitration device
- solid phase
- distribution
- concentration distribution
- 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.)
- Active
Links
Abstract
本发明是一种选择性催化还原脱硝装置入口参数确定方法,该方法以锅炉炉膛、水平烟道和尾部烟道为计算对象,借助Fluent软件进行全系统燃烧、传热和流动过程的数值模拟,获得不同工况下省煤器出口即SCR脱硝装置入口的烟气参数,包括速度分布、温度分布以及烟气成分的浓度分布;针对锅炉炉膛、水平烟道和尾部烟道(至省煤器出口)全系统,进行燃烧、传热和流动过程的数值模拟,并考虑固相颗粒运动,获得不同工况条件下省煤器出口即SCR脱硝装置入口的固相颗粒参数,解决以往因入口烟气参数设置不合理而造成的SCR脱硝装置设计计算与工程实际偏差过大、不考虑固相颗粒运动而造成催化剂磨损积灰等问题,能够有效提高SCR反应器脱硝效率、延长催化剂寿命。
Description
技术领域
本发明涉及一种选择性催化还原脱硝装置入口参数确定方法,属于燃煤锅炉大气污染物控制技术领域。
背景技术
选择性催化还原(SCR)技术因其脱硝效率高、二次无污染等显著优点,成为燃煤电厂烟气脱硝技术的首选。采用CFD技术对SCR脱硝装置的实际烟气流动和氨喷射过程进行数值模拟,已成为SCR工程设计的重要内容。
目前,常规的流场模拟方法都采用锅炉运行参数和局部测试参数作为计算初始条件。在进行边界条件的设置时,通常采用理想的烟气参数,即均匀的烟气速度、温度和成分浓度分布,如通过加权平均测点温度表示SCR脱硝装置入口温度、通过测量单点浓度表示SCR脱硝装置入口浓度、通过烟气流量的面平均计算表示SCR脱硝装置入口速度,然而,已有的研究表明,不管是烟气的温度分布、浓度分布,还是经尾部受热面整流作用后的速度分布,均存在较大的不均匀性,因此,如果采用均匀性入口边界条件,势必带来与工程实际偏差较大的模拟结果,以此为依据设计的导流构件显然无法满足流场均匀性要求,进而导致目前运行SCR脱硝装置的脱硝效率和氨逃逸率普遍不达标。
而且,几乎所有的流场模拟方法仅仅考虑了气相的流动,并以计算获得的流动特征参数(速度相对标准偏差≤15%、烟气入射角≤±10°、温度分布最大绝对偏差≤±10℃、氨氮比相对标准偏差≤5%)作为流场均匀性的评价指标,而忽略了固相的作用,但在实际过程中,固相颗粒自身在重力场中具有运动特性,在气相场的作用下又会产生颗粒碰撞、团聚、解聚和分离等现象,这些运动又会反作用于气相并影响其流动,尤其是对于燃用高灰份煤种的锅炉,因此,忽略固相作用的评价指标会导致不合理的流场设计,这也是目前运行SCR脱硝装置催化剂磨损和积灰堵塞的一个重要原因。
通过以上分析可知,常规的SCR流场模拟方法存在以均匀入流代替非均匀入流、仅考虑气相运动而忽略固相运动等问题,为克服这些问题,本发明提出了一种选择性催化还原脱硝装置入口参数确定方法,旨在提高SCR脱硝装置的脱硝效率、降低SCR脱硝装置的氨逃逸率、改善SCR脱硝装置的催化剂积灰和磨损。
发明内容
发明目的:本发明针对以上技术问题,提供了一种选择性催化还原脱硝装置入口参数确定方法,用于解决以往SCR脱硝装置流场数值模拟时因入口烟气参数设置不合理而造成的SCR反应器设计计算与工程实际偏差过大、未考虑固相颗粒运动而造成催化剂磨损积灰等问题。
技术方案:为解决上述技术问题,本发明提供一种选择性催化还原脱硝装置入口参数确定方法,该方法以锅炉炉膛、水平烟道和尾部烟道(至省煤器出口)为计算对象,借助Fluent软件进行全系统燃烧、传热和流动过程的数值模拟,并考虑固相颗粒运动,获得不同工况(主要指煤种、负荷变化)下省煤器出口即SCR脱硝装置入口的烟气参数,包括气固两相的速度分布、温度分布以及烟气成分的浓度分布;
所述的方法步骤为:
1)根据实际运行锅炉设备的结构和尺寸,采用Gambit软件建立包含锅炉炉膛、水平烟道和尾部烟道的全系统三维物理模型;
2)采用Fluent软件,以实际运行设备的入口条件为物理模型入口边界条件,采用可实现k-ε湍流模型和离散相模型,数值计算不同工况下煤粉在炉内的燃烧过程,烟气在锅炉炉膛、水平烟道和尾部烟道的传热、流动过程,获得不同工况下省煤器出口处即SCR脱硝装置入口烟气成分(NOx、O2和未燃尽碳)的温度、速度和浓度分布;
3)采用Fluent软件,以除尘器中采样得到的飞灰物性参数作为入口边界条件,数值计算取样飞灰的流动过程,获得不同工况下省煤器出口处即SCR脱硝装置入口飞灰在烟道截面上的粒径分布;
4)结合省煤器出口截面步骤2)中未燃尽碳的浓度分布和步骤3)中飞灰的粒径分布,并比对现场测试数据,得到省煤器出口处即SCR脱硝装置入口不同粒径固相颗粒的浓度分布。
省煤器出口处不同粒径固相颗粒的浓度分布通过以下步骤获得:
1)建立锅炉炉膛、水平烟道和尾部烟道的三维物理模型,分别设置燃烧器入口和省煤器出口为入口边界和出口边界,并将出口边界划分为若干子区域;
2)采用Fluent软件数值计算不同工况下煤粉在炉内的燃烧过程,烟气在炉膛和尾部烟道的传热、流动过程,获得不同工况下省煤器出口截面各子区域未燃尽碳的浓度分布;
3)对除尘器灰斗的飞灰进行取样,并进行粒径分析,获得不同工况下对应的取样飞灰粒径分布;
4)在步骤1)中所述物理模型的基础上,以取样飞灰的物性参数作为固相入口边界条件,选择气固两相流模型,数值计算取样飞灰在锅炉炉膛、水平烟道和尾部烟道的流动过程,获得不同工况下省煤器出口截面各子区域飞灰颗粒的粒径分布;
5)整理不同工况下省煤器出口截面各子区域的未燃尽碳的浓度分布和飞灰颗粒的粒径分布的数据,借助概率统计、加权平均、最小二乘法等数据处理方法,计算得到省煤器出口截面不同粒径固相颗粒的浓度分布;
6)将步骤4)中的计算结果与省煤器出口取样测试数据进行比对,引入修正因子,获得符合实际情况的省煤器出口截面不同粒径固相颗粒的浓度分布。
有益效果:与常规的SCR脱硝装置入口参数确定方法相比较,本发明具有如下的特色及优点:
1、传统SCR流场模拟方法采用均匀性的入口边界条件,而本发明通过对锅炉炉膛、水平烟道和尾部烟道(至省煤器出口)全系统进行数值计算获得更符合实际的SCR脱硝装置入口条件,克服了传统方法由于通过加权平均测点温度表示SCR脱硝装置入口温度、通过测量单点浓度表示SCR脱硝装置入口浓度、通过烟气流量的面平均计算表示SCR脱硝装置入口速度而引起模拟结果与工程实际偏差较大的缺陷,以此为依据设计的导流构件更能满足流场均匀性要求,从而提高SCR脱硝装置的脱硝效率,降低氨逃逸率。
2、传统SCR流场模拟方法仅仅考虑气相的流动,而本发明通过对锅炉炉膛至省煤器全系统燃烧、传热和流动过程以及冷态飞灰流动过程进行数值计算获得SCR脱硝装置入口固相颗粒的流动特性,为基于气固两相流动的SCR流场优化设计奠定基础,避免传统方法由于忽略固相颗粒重力场下自身运动特性和固相颗粒与气相耦合作用下的运动特性而引起流场设计不合理的缺点,改善了催化剂积灰和磨损状况,对于燃用高灰份煤种的锅炉其改善效果尤为明显。
附图说明
图1为本发明方法中的锅炉炉膛至省煤器出口的全系统示意图,其中有:旋流燃烧器1、锅炉炉膛2、过热器和再热器3、水平烟道4、尾部烟道5、省煤器6和省煤器出口7。
具体实施方式
所述的方法包括以下步骤:
1)根据实际运行锅炉设备的结构和尺寸,采用Gambit软件建立包含锅炉炉膛、水平烟道和尾部烟道(至省煤器出口)的全系统三维物理模型,设置燃烧器入口和省煤器出口为入口边界和出口边界,并将出口边界划分为若干子区域;
2)采用Fluent软件,以实际运行设备的入口条件为物理模型入口边界条件,采用可实现k-ε湍流模型和离散相模型计算不同工况下煤粉在炉内的燃烧过程,烟气在锅炉炉膛、水平烟道和尾部烟道的传热、流动过程,获得不同工况下省煤器出口即SCR脱硝装置入口烟气成分(NOx、O2和未燃尽碳)的温度、速度和浓度分布;
3)对除尘器灰斗的飞灰进行取样,并进行粒径分析,获得不同工况下对应的取样飞灰粒径分布;
4)采用Fluent软件,在步骤1)中所述物理模型的基础上,以取样飞灰的物性参数为入口固相边界条件,选择气固两相流模型,数值计算取样飞灰在炉膛、水平烟道和尾部烟道(至省煤器出口)的流动过程,获得不同工况下省煤器出口截面各子区域飞灰颗粒的粒径分布;
5)整理中不同工况下步骤2)中省煤器出口截面上各子区域的未燃尽碳浓度分布和步骤4)飞灰粒径分布的数据,借助概率统计、加权平均、最小二乘法等数据处理方法,计算得到省煤器出口截面不同粒径固相颗粒的浓度分布;
6)并将步骤5)中的计算结果与省煤器出口取样测试数据进行比对,引入修正因子,得到符合实际情况的省煤器出口截面不同粒径固相颗粒的浓度分布。
7)由步骤2)中得到的省煤器出口处气相参数和步骤6)得到的省煤器出口处固相参数,最终确定SCR脱硝装置的入口参数。
以下参照附图详细说明本发明实施。本实施案例在以本技术方法为前提下进行实施,给出了详细的实施方式和操作过程,但本发明的保护范围不仅限于下述的实施例。
本实施案例针对某电厂660MW机组锅炉进行锅炉炉膛至省煤器全系统燃烧、传热和流动过程的数值模拟,以确定SCR脱硝装置入口参数。该锅炉旋流燃烧器1为对冲布置,在额定负荷工况下,总风量为2292t/h,燃煤消耗308t/h。
本实施案例SCR脱硝装置入口参数确定方法的步骤如下:
1)根据实际运行设备简化结构建立包含锅炉炉膛2、水平烟道4和尾部烟道5全系统的三维物理模型,分析给定的条件参数(风率比、煤粉细度、煤质分析等)计算物理模型的边界条件为数值模拟做准备;
2)采用Fluent软件计算全系统的燃烧、传热和流动过程,其中,气相湍流采用RNGk-ε双方程模型,颗粒相采用随机轨道模型,煤粉燃烧采用非预混燃烧模型,辐射传热采用以热流法为基础的P1模型,过热器和再热器3均采用多孔介质模型,NOx的生成采用燃烧后处理方法;
3)计算不同负荷下煤粉在炉内的燃烧过程,烟气在锅炉炉膛2、水平烟道4和尾部烟道5的流动、传热过程,获得不同负荷省煤器出口7处烟气的温度、速度和成分(NOx、O2和未燃尽碳)的浓度分布;
4)对不同负荷下尾部除尘器灰斗的飞灰进行取样,并进行粒径分析,获得不同工况对应的取样飞灰粒径分布;
5)根据取样飞灰粒径分布,对锅炉炉膛2、水平烟道4和尾部烟道5全系统进行气固两相流动的数值计算,模拟不同负荷下煤粉燃烧后飞灰在锅炉炉膛2、水平烟道4和尾部烟道5的流动过程,获得不同负荷下省煤器出口7处飞灰在的粒径分布;
6)根据步骤3)中未燃尽碳在烟道截面的浓度分布和步骤5)中飞灰在省煤器出口7截面的粒径分布,借助概率统计、加权平均、最小二乘法等数据处理方法,计算出省煤器出口7处不同粒径固相颗粒的浓度分布;
7)并将步骤5)中的计算结果与省煤器出口取样测试数据进行比对,引入修正因子,得到符合实际情况的省煤器出口截面不同粒径固相颗粒的浓度分布。
8)将步骤3)中省煤器出口7处烟气的温度、速度和成分(NOx、O2和未燃尽碳)的浓度分布参数和步骤6)中省煤器出口7处不同粒径固相颗粒的浓度分布参数作为符合实际的SCR脱硝装置的入口气固两相边界条件。
Claims (1)
1.一种选择性催化还原脱硝装置入口参数确定方法,其特征在于,该方法以锅炉炉膛、水平烟道和尾部烟道为计算对象,借助Fluent软件进行全系统燃烧、传热和流动过程的数值模拟,并考虑固相颗粒运动,获得不同工况下省煤器出口即SCR脱硝装置入口的烟气参数,包括气固两相的速度分布、温度分布以及烟气成分的浓度分布;
所述的方法步骤为:
1)根据实际运行锅炉设备的结构和尺寸,采用Gambit软件建立包含锅炉炉膛、水平烟道和尾部烟道的全系统三维物理模型;
2)采用Fluent软件,以实际运行设备的入口条件为物理模型入口边界条件,数值计算不同工况下煤粉在炉内的燃烧过程,烟气在炉膛和尾部烟道的传热、流动过程,获得不同工况下省煤器出口处即SCR脱硝装置入口烟气成分的温度、速度和浓度分布;
3)采用Fluent软件,以除尘器中采样得到的飞灰物性参数作为入口边界条件,数值计算取样飞灰的流动过程,获得不同工况下省煤器出口处即SCR脱硝装置入口飞灰在烟道截面上的粒径分布;
4)结合省煤器出口截面步骤2)中未燃尽碳的浓度分布和步骤3)中飞灰的粒径分布,并比对现场测试数据,得到省煤器出口处即SCR脱硝装置入口不同粒径固相颗粒的浓度分布;
其中,
所述的省煤器出口处不同粒径固相颗粒的浓度分布,通过以下步骤获得:
4.1)建立锅炉炉膛、水平烟道和尾部烟道的三维物理模型,分别设置燃烧器入口和省煤器出口为入口边界和出口边界,并将出口边界划分为若干子区域;
4.2)采用Fluent软件数值计算不同工况下煤粉在炉内的燃烧过程,烟气在炉膛和尾部烟道的传热、流动过程,获得不同工况下省煤器出口截面各子区域未燃尽碳的浓度分布;
4.3)对除尘器灰斗的飞灰进行取样,并进行粒径分析,获得不同工况下对应的取样飞灰粒径分布;
4.4)在步骤4.1)中所述物理模型的基础上,以取样飞灰的物性参数作为固相入口边界条件,数值计算飞灰在锅炉炉膛、水平烟道和尾部烟道的流动过程,获得不同工况下省煤器出口截面各子区域飞灰颗粒的粒径分布;
4.5)整理不同工况下省煤器出口截面各子区域的未燃尽碳的浓度分布和飞灰颗粒的粒径分布的数据,借助概率统计、加权平均、最小二乘法数据处理方法,计算得到省煤器出口截面不同粒径固相颗粒的浓度分布;
4.6)将步骤4.4)中的计算结果与省煤器出口取样测试数据进行比对,引入修正因子,得到符合实际情况的省煤器出口截面不同粒径固相颗粒的浓度分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164248.1A CN103926833B (zh) | 2014-04-22 | 2014-04-22 | 一种选择性催化还原脱硝装置入口参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410164248.1A CN103926833B (zh) | 2014-04-22 | 2014-04-22 | 一种选择性催化还原脱硝装置入口参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103926833A CN103926833A (zh) | 2014-07-16 |
CN103926833B true CN103926833B (zh) | 2016-06-08 |
Family
ID=51145093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410164248.1A Active CN103926833B (zh) | 2014-04-22 | 2014-04-22 | 一种选择性催化还原脱硝装置入口参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103926833B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104607042B (zh) * | 2014-12-26 | 2017-11-17 | 东南大学 | 一种基于约束预测控制的scr脱硝系统及方法 |
CN105441982A (zh) * | 2016-01-14 | 2016-03-30 | 江苏立典机床科技有限公司 | 一种基于数值模拟的电解铝烟气净化系统 |
CN107149875B (zh) * | 2017-05-15 | 2023-04-11 | 大唐环境产业集团股份有限公司 | 用于scr冷态模化试验的非均匀性喷射装置及方法 |
CN107764605A (zh) * | 2017-09-07 | 2018-03-06 | 大唐东北电力试验研究所有限公司 | 粉尘采样载体的优选、储存和转运方法及滤筒存放装置 |
CN109289517A (zh) * | 2018-09-13 | 2019-02-01 | 大唐东北电力试验研究院有限公司 | 选择性催化还原脱硝装置催化剂防磨方法 |
CN111437721B (zh) * | 2020-04-03 | 2023-01-17 | 大唐珲春发电厂 | 一种选择性催化还原脱硝装置混合器多工况性能优化方法 |
CN111760445A (zh) * | 2020-06-28 | 2020-10-13 | 淮阴工学院 | 一种脱硫塔流场模拟方法 |
CN112528547A (zh) * | 2020-11-20 | 2021-03-19 | 西安西热锅炉环保工程有限公司 | 基于气固两相流数值研究解决锅炉尾部烟道磨损的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1699951A (zh) * | 2005-05-12 | 2005-11-23 | 西安建筑科技大学 | 除尘器气流分布数字实验方法 |
CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
CN102073273A (zh) * | 2010-05-13 | 2011-05-25 | 安徽兴泉喷射液化装备有限责任公司 | 基于Optimus软件与Fluent软件相结合对蒸汽喷射液化装置优化的方法 |
CN103077269A (zh) * | 2012-12-28 | 2013-05-01 | 厦门龙净环保技术有限公司 | 一种电袋复合除尘器气流均布优化设计方法 |
CN103218504A (zh) * | 2013-05-15 | 2013-07-24 | 西安近代化学研究所 | 固体推进剂配方设计用羽流场数据结构网格化处理方法 |
EP2623966A1 (en) * | 2010-09-28 | 2013-08-07 | Mitsubishi Heavy Industries, Ltd. | Gas analysis device, mercury removal system, gas analysis method, and method of removing mercury from exhaust gas |
-
2014
- 2014-04-22 CN CN201410164248.1A patent/CN103926833B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1699951A (zh) * | 2005-05-12 | 2005-11-23 | 西安建筑科技大学 | 除尘器气流分布数字实验方法 |
CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
CN102073273A (zh) * | 2010-05-13 | 2011-05-25 | 安徽兴泉喷射液化装备有限责任公司 | 基于Optimus软件与Fluent软件相结合对蒸汽喷射液化装置优化的方法 |
EP2623966A1 (en) * | 2010-09-28 | 2013-08-07 | Mitsubishi Heavy Industries, Ltd. | Gas analysis device, mercury removal system, gas analysis method, and method of removing mercury from exhaust gas |
CN103077269A (zh) * | 2012-12-28 | 2013-05-01 | 厦门龙净环保技术有限公司 | 一种电袋复合除尘器气流均布优化设计方法 |
CN103218504A (zh) * | 2013-05-15 | 2013-07-24 | 西安近代化学研究所 | 固体推进剂配方设计用羽流场数据结构网格化处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103926833A (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103926833B (zh) | 一种选择性催化还原脱硝装置入口参数确定方法 | |
CN103927420B (zh) | 选择性催化还原脱硝装置的气固两相流场与喷氨优化方法 | |
Choi et al. | Process design and optimization of MEA-based CO2 capture processes for non-power industries | |
Kang et al. | Study of the performance, simplification and characteristics of SNCR de-NOx in large-scale cyclone separator | |
CN103235842B (zh) | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 | |
CN103969073B (zh) | 烟风模拟换热装置 | |
Liu et al. | Optimization and numerical simulation of the flow characteristics in SCR system | |
CN103148508B (zh) | 前后对冲旋流燃烧煤粉锅炉的二次风配风方法及系统 | |
CN103968413B (zh) | 锅炉改造后不同负荷下燃烧特性的获取方法 | |
CN114239430A (zh) | 一种基于数值模拟的炉膛出口NOx预测方法及系统 | |
Zhou et al. | Numerical investigation of gas-particle flow in the primary air pipe of a low NOx swirl burner–The DEM-CFD method | |
CN106548032A (zh) | 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 | |
Deng et al. | Effects of flue gas recirculation on combustion and heat flux distribution in 660 MW double-reheat tower-type boiler | |
CN203848311U (zh) | 火电厂烟道高温段烟尘控制装置 | |
CN106918464B (zh) | 600mw超超临界前后对冲锅炉运行诊断方法及装置 | |
CN116128136A (zh) | 一种基于LSO-CatBoost的燃煤电厂锅炉NOX排放预测方法 | |
Du et al. | Influence of the parallel oil‐secondary air and F‐layer secondary air distribution on the flow, combustion, and NOx generation characteristics of FW down‐fired boilers retrofitted with a stable combustion technology | |
Liu et al. | Review of simulation research on pulverized coal combustion in industrial boilers | |
Pryiomov et al. | Increasing the energy efficiency of cyclone dust collectors | |
Zhang et al. | Integrated Gas and Solid Model to Determine the Flow Performance in the SCR Reactor | |
Ran et al. | Numerical study on optimum designing of the air distribution structure of a new cyclone combustor | |
CN103968412A (zh) | 锅炉改造后不同ccofa风与sofa风比例下燃烧特性的获取方法 | |
Guihao et al. | Dynamic NO x emission prediction based on composite models adapt to different operating conditions of coal-fired utility boilers | |
CN204213918U (zh) | 一种多股流燃尽风喷嘴装置 | |
CN112944384B (zh) | 一种超临界二氧化碳锅炉空气预热器模块与脱硝装置耦合系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |