CN106548032B - 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 - Google Patents
基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 Download PDFInfo
- Publication number
- CN106548032B CN106548032B CN201610981223.XA CN201610981223A CN106548032B CN 106548032 B CN106548032 B CN 106548032B CN 201610981223 A CN201610981223 A CN 201610981223A CN 106548032 B CN106548032 B CN 106548032B
- Authority
- CN
- China
- Prior art keywords
- boiler
- flue gas
- flue
- concentration
- grid
- 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
- 238000000034 method Methods 0.000 title claims abstract description 61
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 title claims abstract description 43
- 239000003546 flue gas Substances 0.000 title claims abstract description 43
- 238000001514 detection method Methods 0.000 title claims abstract description 16
- 238000002485 combustion reaction Methods 0.000 claims abstract description 22
- 238000012546 transfer Methods 0.000 claims abstract description 12
- 239000012530 fluid Substances 0.000 claims abstract description 6
- 239000002245 particle Substances 0.000 claims description 33
- 239000003245 coal Substances 0.000 claims description 18
- 229910052760 oxygen Inorganic materials 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 239000007787 solid Substances 0.000 claims description 11
- 239000002817 coal dust Substances 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 10
- 230000005855 radiation Effects 0.000 claims description 10
- 238000009792 diffusion process Methods 0.000 claims description 9
- 239000007789 gas Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 6
- 208000011580 syndromic disease Diseases 0.000 claims description 6
- 239000000571 coke Substances 0.000 claims description 4
- 238000010304 firing Methods 0.000 claims description 4
- 239000000779 smoke Substances 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000002950 deficient Effects 0.000 claims description 3
- 239000002594 sorbent Substances 0.000 claims description 3
- 235000019504 cigarettes Nutrition 0.000 claims description 2
- 239000003517 fume Substances 0.000 claims description 2
- 230000003796 beauty Effects 0.000 claims 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 8
- 239000001301 oxygen Substances 0.000 description 8
- 238000005259 measurement Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000000843 powder Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 238000010025 steaming Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Regulation And Control Of Combustion (AREA)
Abstract
本发明涉及一种基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，该方法通过对锅炉省煤器出口烟气CO、O2浓度和烟气压力的网格式多点在线检测，采用计算流体力学的方法建立炉膛的流动、燃烧和传热模型，以追溯烟气CO和O2在炉膛、烟道中的运动轨迹，进而重建烟气CO和O2浓度场，为锅炉燃烧调整提供有效可靠的关键依据。
Description
技术领域
本发明涉及一种基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法。
背景技术
社会的发展和进步，火电厂面临越来越大的压力。电力需求不断增加，发电品质要求日益苛刻，传统的电压、频率已不是唯一的考核指标，环境因素逐渐纳入其中，形成了新的衡量标准；煤质波动以及大型机组日益频繁的参加调频调峰，全工况变负荷运行以及长期中低负荷运行，都对机组运行水平提出了更高的要求。目前，发电厂越来越多的面临来自电网、环境和自身的压力，节能、降耗、减排的呼声高涨，燃烧优化技术日益得到重视。锅炉燃烧效率的高低，进入炉膛的燃料能否充分燃烧，将直接影响锅炉的整体效率与经济性。因此通过燃烧控制实现锅炉最高效地燃烧是提高锅炉整体效率的有效途径。
电站煤粉锅炉燃烧控制设计基本思路是以合适的过量空气系数为原则，实现锅炉燃烧最优控制即是以保证合适的过量空气系数。锅炉的燃烧量变化时，相应地改变送风量，以保证合适的过剩空气系数，减少锅炉未完全燃烧损失和排烟热损失，使锅炉运行于最佳燃烧工况。单纯测量O2并不能反映炉内煤粉和空气混合状况的好坏，即使O2足够，若混合不好等原因，也会使不完全燃烧损失增大。锅炉实际运行中，只会对省煤器出口以后的烟道气的O2进行有限点的在线检测，样本点非常有限。对于大容量锅炉而言，烟气成分分布不均容易使氧量检测值不具备代表性。另外，氧量表本身由于漏风、零位漂移等原因也不能反映烟气中的真实氧量值。基于上述原因，锅炉实际运行过程中，运行人员只能以较大过量空气系数控制燃烧，非常不经济，致使锅炉运行容易发生炉膛高温腐蚀、火焰偏斜、燃烧滞后、结焦结渣、燃烧效率偏低等运行困扰。如何寻找一个完善的检测手段来对运行中的锅炉进行调节以解决锅炉运行存在的问题并提高锅炉的实际运行效率是十分必要的。
已有的相关实现方案：烟气含氧量软测量技术：该技术选用蒸汽流量、磨煤机给煤量、总风量、总一次风压、送风机电流、引风机电流、风箱与炉膛差压、给水及排烟温度等11个参数作为输入信号，将烟气含氧量作为输出参数，应用前向神经网络建立烟气含氧量软测量模型。在各种不同负荷工况下，测取多组实际运行数据，抽取1～2组数据作为模型检测使用，其余作为训练样本。此种方法一定程度上解决了氧量测量滞后、误差大的问题。
烟气含氧量软测量技术——软测量结果的准确性依赖于边界试验的完整性和全面性，往往由于入炉煤煤质的波动、工况的变化导致误差增大。在燃烧调整过程中，软测量技术的应用很难大规模的有效推广。
发明内容
为解决上述技术问题，本发明的目的是提供一种基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，采用计算流体力学的方法建立炉膛的流动、燃烧和传热模型，以追溯烟气CO和O2在炉膛、烟道中的运动轨迹，进而重建烟气CO和O2浓度场，为锅炉燃烧调整提供有效可靠的关键依据。
本发明基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，包括：获取电站锅炉的结构、操作数据，对锅炉省煤器出口烟气CO、O2浓度和烟气压力的网格式多点在线检测，采用计算流体力学的方法建立炉膛的流动、燃烧和传热模型，以追溯烟气CO和O2在炉膛、烟道中的运动轨迹，重建烟气CO和O2浓度场。
进一步地，所述方法具体包括：
获取电站锅炉的结构、操作数据，所述的电站锅炉的结构包括：
所述的电站锅炉的机构、操作数据包括电站锅炉冷灰斗、炉膛、折焰角、燃烧器喷口、省煤器的详细尺寸结构；所述的操作数据至少包括：每台锅炉的总风量、总煤量，每个燃烧器的燃煤量、风量、一次风温、二次风温、一次风率、二次风率；
基于烟道网格多点布置取样方法，从烟道中取烟气的代表性样品，即探头单元的安装在预定的位置通过探头单元获取烟道中的烟气样品；
建立炉膛及烟道的计算网格，其中全锅炉的计算域选至锅炉省煤器出口处，计算域化分为四个部分，分为冷灰斗区域I、燃烧器区域II、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV；冷灰斗I、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV采用结构化网格，燃烧器区域I I 采用非结构化网格；
各区域采用三维稳态计算，采用流动、传热和燃烧的非耦合求解，计算中采用有限容积法对微分控制方程和输运方程离散为代数方程组进行数值求解。以各典型工况下的CO、O2和烟气压力的实际测量值作为输入条件，燃煤量、总风量、一次风温、二次风温、一次风率、二次风率和漏风率等作为边界条件多次迭代直至满足收敛，追溯烟气CO、O2在炉膛至省煤器出口中的运动轨迹，重建锅炉燃烧过程中CO、O2等三维分布场。
进一步地，以美国机械工程师协会(ASME)标准的规定对于烟气取样：测点的数量及布置取决于烟道的尺寸和形状，取样点间隔应不超过915mm，且所用点的总数应不少于4点；对不同容量的锅炉，在离锅炉炉膛出口最近的尾部烟道竖井上选址布置，安装的取样点个数严格按照ASME标准执行，点数选择如下所示：机组容量300MW，安装点数8～16；机组容量600MW，安装点数16～20；机组容量1000MW，安装点数20～24。
进一步地，煤粉燃烧模型采用扩散/反应联合控制模型，主要包括两项，分别为扩散项和传热项。扩散项为：
其中q为焦炭消耗速率，A指前因子，E为表观活化能，Pox，s为颗粒表面处的O2分压，nox为反应阶数。传热项为：
其中mc为颗粒质量，Tg，Tp，Tr分别为气相、颗粒和背景辐射温度，cp为固体比热，fh为被固体吸收的反应热份额，Hreac为气固反应热，εp为颗粒黑度，σ为斯蒂芬-玻尔兹曼常数。
进一步地，煤粉颗粒的运动轨迹，采用颗粒随机轨道模型,根据射流的数目决定煤粉跟踪颗粒数；假定射流数目为20,其含义是把之前跟踪的颗粒分成 20个,每个颗粒占有的流量就是原来的1/20,然后对这20个颗粒赋予不同的速度脉动值,求得各自的颗粒轨迹。
进一步地，在对辐射传热模拟时选用以热流法为基础的P1辐射模型。与离散坐标法相比，P1辐射模型将辐射散射考虑在内。
进一步地，选用标准k-ε双方程湍流模型，非预混燃烧守恒标量的PDF模型以及FLUENT软件提供的污染物生成模型；对于速度与压力的耦合问题，采用 SIMPLE的算法；对壁面处理采用了壁面函数法，；同时为了保持计算的稳定性，采用欠松弛的方法。
借由上述方案，本发明至少具有以下优点：
使用基于CFD的方法，高效省时，可以利用现有很多模型库进行优化处理。
直观反应烟气CO、O2在炉膛中的运动轨迹和分布规律，可以清晰掌握炉内风粉混合状况，有利于有针对性的单独调整单个燃烧器的风粉配比。
通过网格式多点检测技术可有效提升浓度场的准确性和有效性。
同时测量的网格式多点压力构造的实际速度场，可以对基于CFD得出的速度场进行有效校验。
为运行人员调整燃烧提供了有效、可靠、全面的关键依据。
上述说明仅是本实用新型技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
图1是本发明的锅炉结构示意图；
图2是本发明炉膛及烟道结构示意图；
图3是本发明炉膛及烟道网格划分结构示意图；
图4是本发明燃烧器区域网格化结构示意图；
图5是本发明独立求解方法流程图；
图6是本发明耦合求解方法流程图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
本实施例基于锅炉烟道网格检测的烟气CO和O2浓度场的重建方法，所述方法包括：
[1]确定电站锅炉的结构、操作数据，包括电站锅炉冷灰斗、炉膛、折焰角、燃烧器喷口、省煤器的详细尺寸结构，每台锅炉的总风量、总煤量，每个燃烧器的燃煤量、风量、一次风温、二次风温、一次风率、二次风率等。
[2]烟道网格多点布置取样方法。从烟道中取烟气的代表性样品，即探头单元的安装位置选择是保证分析结果准确度的关键。对于烟气取样代表性的问题， ASME(美国机械工程师协会)标准的规定：“测点的数量及布置取决于烟道的尺寸和形状，取样点间隔应不超过3英尺(915mm)，且所用点的总数应不少于4 点。”对不同容量的锅炉，在离锅炉炉膛出口最近的尾部烟道竖井上选址布置。安装的取样点个数严格按照ASME标准执行，一般点数选择如下表所示：
| 机组容量 | 300MW | 600MW | 1000MW |
| 安装点数 | 8～16 | 16～20 | 20～24 |
[3]建立炉膛及烟道的计算网格。全锅炉的计算域选至锅炉省煤器出口处，计算域化分为四个部分，分为冷灰斗区域I、燃烧器区域II、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV。冷灰斗I、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV采用结构化网格，燃烧器区域I I 采用非结构化网格，有效简化整个计算域网格数目，有效提升了系统的收敛速度。锅炉燃烧器喷口均不进行简化，严格按照锅炉实际运行中各喷口进风方向进行数值模拟。
[4]煤的物理结构和化学组分都十分复杂，煤的燃烧过程也会经历诸多复杂的过程。各部分之间又没有明显界限，煤粉可能会同时经历两个或多个过程。所以对煤粉的燃烧过程进行合理的描述并不容易，对各阶段选取恰当的模型十分重要。在数值计算时，建立挥发分热解模型和焦炭燃烧模型以给出能量平衡方程和组份平衡方程中的源项，使控制方程组得以封闭。煤中挥发分的析出过程也称为煤的热解，是一个十分复杂的过程，涉及物理汽化、化学键断裂、不稳定中间产物的产生及演变等诸多过程。合理的挥发分热解模型应该能够比较完整的描述煤粉在加热阶段，从中析出的产物种类、数量及其析出速率。
煤粉燃烧模型：采用扩散/反应联合控制模型，主要包括两项，分别为扩散项和传热项。扩散项为：
其中q为焦炭消耗速率，A指前因子，E为表观活化能，Pox，s为颗粒表面处的O2分压，nox为反应阶数。传热项为：
其中mc为颗粒质量，Tg，Tp，Tr分别为气相、颗粒和背景辐射温度，cp为固体比热，fh为被固体吸收的反应热份额，Hreac为气固反应热，εp为颗粒黑度，σ为斯蒂芬-玻尔兹曼常数。
煤粉炉内的燃烧过程一般都会涉及到典型的气固两相流动。本发明的气固两相流模型选用随机轨道模型。随机轨道模型以气相湍流的统观模型为基础，以考虑气相湍流基本结构特征；由颗粒相的瞬态动量方程随机给定脉动速度，以考虑湍流扩散对颗粒的影响；用蒙特卡洛法计算颗粒的运动轨道及其沿程的复杂变化。这一方法的做法基本都是在求解气态流体时均场的基础上对颗粒相的瞬态动量方程求解。随机轨道模型在研究伴随有蒸发、挥发及异相化学反应等复杂经历的颗粒在流场中的变化时有很好的效果，是目前在两相湍流模拟和固体颗粒燃烧模拟中被最广泛采用的模型。
煤粉颗粒的运动轨迹：采用颗粒随机轨道模型,根据射流的数目决定煤粉跟踪颗粒数。假定射流数目为20,其含义是把之前跟踪的颗粒分成20个,每个颗粒占有的流量就是原来的1/20,然后对这20个颗粒赋予不同的速度脉动值,求得各自的颗粒轨迹。
本发明在对辐射传热模拟时选用以热流法为基础的P1辐射模型。与离散坐标法相比，P1辐射模型将辐射散射考虑在内，在模拟较复杂几何结构中的辐射问题时有一定优势。这种方法相对比较成熟，已被国内外许多学者用于炉内辐射的模拟。
本发明选用标准k-ε双方程湍流模型，非预混燃烧守恒标量的PDF模型以及FLUENT软件提供的污染物生成模型。对于速度与压力的耦合问题，采用SIMPLE 的算法；对壁面处理采用了壁面函数法，保证了数值计算的适应性；同时为了保持计算的稳定性，采用欠松弛的方法。
[5]各区域采用三维稳态计算，采用流动、传热和燃烧的非耦合求解，计算中采用有限容积法对微分控制方程和输运方程离散为代数方程组进行数值求解。以各典型工况下的CO、O2和烟气压力的实际测量值作为输入条件，燃煤量、总风量、一次风温、二次风温、一次风率、二次风率和漏风率等作为边界条件多次迭代直至满足收敛，追溯烟气CO、O2在炉膛至省煤器出口中的运动轨迹，重建锅炉燃烧过程中CO、O2等三维分布场。
求解方法：具体的求解方法可采用下面图5、图6方法中的一种。通过计算每个计算网格中的各项参数，获得锅炉计算域中的浓度场和速度场。
此重建方法具有较好的适用性以及推广性，可用于各种型式的锅炉，对炉内各燃烧区域的烟气CO和O2的浓度场重建是一种行之有效的方法。
本发明结构化网格：指网格区域内所有的内部点都具有相同的毗邻单元。它可以很容易地实现区域的边界拟合，适于流体和表面应力集中等方面的计算。它的主要优点是网格生成的速度快、质量好，数据结构简单。它的最典型的缺点是适用的范围比较窄。尤其随着近几年的计算机和数值方法的快速发展，人们对求解区域的复杂性的要求越来越高，在这种情况下，结构化网格生成技术就显得力不从心了。
非结构化网格：指网格区域内的内部点不具有相同的毗邻单元。即与网格剖分区域内的不同内点相连的网格数目不同。从定义上可以看出，结构化网格和非结构化网格有相互重叠的部分，即非结构化网格中可能会包含结构化网格的部分。非结构化网格可以解决任意形状和任意连通区域的网格剖分。
以上所述仅是本发明的优选实施方式，并不用于限制本发明，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变型，这些改进和变型也应视为本发明的保护范围。
Claims (6)
1.一种基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，包括：获取电站锅炉的结构、操作数据，对锅炉省煤器出口烟气CO、O2浓度和烟气压力进行网格式多点在线检测，采用计算流体力学的方法建立炉膛内煤粉的流动、燃烧和传热模型，以追溯烟气CO和O2在炉膛、烟道中的运动轨迹，重建烟气CO和O2浓度场；
其中，获取电站锅炉的结构、操作数据，所述的电站锅炉的结构包括：
电站锅炉冷灰斗、炉膛、折焰角、燃烧器喷口、省煤器的详细尺寸结构；所述的操作数据至少包括：每台锅炉的总风量、总煤量，每个燃烧器的燃煤量、风量、一次风温、二次风温、一次风率、二次风率；
基于烟道网格多点布置取样方法，从烟道中取烟气的代表性样品，即探头单元安装在预定的位置，通过探头单元获取烟道中的烟气样品；
建立炉膛及烟道的计算网格，其中全锅炉的计算域选至锅炉省煤器出口处，计算域划分为四个部分，分为冷灰斗区域I、燃烧器区域II、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV；冷灰斗区域I、燃烧器以上至炉膛出口区域III、炉膛出口至省煤器出口区域IV采用结构化网格，燃烧器区域II采用非结构化网格；
各区域采用三维稳态计算，采用流动、传热和燃烧模型的非耦合求解，计算中采用有限容积法对微分控制方程和输运方程离散为代数方程组进行数值求解，以各典型工况下的CO、O2和烟气压力的实际测量值作为输入条件，燃煤量、总风量、一次风温、二次风温、一次风率、二次风率和漏风率作为边界条件多次迭代直至满足收敛，追溯烟气CO、O2在炉膛至省煤器出口中的运动轨迹，重建锅炉燃烧过程中CO、O2三维分布场。
2.根据权利要求1所述的基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，以美国机械工程师协会标准的规定对于烟气取样：测点的数量及布置取决于烟道的尺寸和形状，取样点间隔应不超过915mm，且所用点的总数应不少于4点；对不同容量的锅炉，在离锅炉炉膛出口最近的尾部烟道竖井上选址布置，安装的取样点个数严格按照美国机械工程师协会标准执行，点数选择如下所示：机组容量300MW，安装点数8～16；机组容量600MW，安装点数16～20；机组容量1000MW，安装点数20～24。
3.根据权利要求2所述的基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，煤粉燃烧模型采用扩散/反应联合控制模型，包括两项，分别为扩散项和传热项，扩散项为：
其中q为焦炭消耗速率，A指前因子，E为表观活化能，Pox，s为颗粒表面处的O2分压，nox为反应阶数，传热项为：
其中mc为颗粒质量，Tg，Tp，Tr分别为气相、颗粒和背景辐射温度，cp为固体比热，fh为被固体吸收的反应热份额，Hreac为气固反应热，εp为颗粒黑度，σ为斯蒂芬-玻尔兹曼常数。
4.根据权利要求3所述的基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，煤粉颗粒的运动轨迹，采用颗粒随机轨道模型，根据射流的数目决定煤粉跟踪颗粒数；假定射流数目为20，其含义是把之前跟踪的颗粒分成20个，每个颗粒占有的流量就是原来的1/20，然后对这20个颗粒赋予不同的速度脉动值，求得各自的颗粒轨迹。
5.根据权利要求4所述的基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，在对辐射传热模拟时选用以热流法为基础的P1辐射模型，与离散坐标法相比，P1辐射模型将辐射散射考虑在内。
6.根据权利要求5所述的基于锅炉烟道网格检测烟气CO和O2浓度的分布场重建方法，其特征在于，选用标准k-ε双方程湍流模型，非预混燃烧守恒标量的PDF模型以及FLUENT软件提供的污染物生成模型；对于速度与压力的耦合问题，采用SIMPLE的算法；对壁面处理采用了壁面函数法同时为了保持计算的稳定性，采用欠松弛的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610981223.XA CN106548032B (zh) | 2016-11-08 | 2016-11-08 | 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610981223.XA CN106548032B (zh) | 2016-11-08 | 2016-11-08 | 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106548032A CN106548032A (zh) | 2017-03-29 |
| CN106548032B true CN106548032B (zh) | 2019-04-05 |
Family
ID=58394735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610981223.XA Active CN106548032B (zh) | 2016-11-08 | 2016-11-08 | 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106548032B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107191918B (zh) * | 2017-06-14 | 2019-07-19 | 华北电力科学研究院有限责任公司 | 锅炉炉底漏风量确定方法及装置 |
| CN108254518B (zh) * | 2017-12-30 | 2021-02-09 | 聚光科技(杭州)股份有限公司 | 密集分布式水质监测系统及方法 |
| CN111312337A (zh) * | 2020-01-17 | 2020-06-19 | 中国地质大学（武汉） | 基于热重实验的单峰可燃物热解反应机理模型确定方法 |
| CN113406282B (zh) * | 2021-06-18 | 2022-08-30 | 国网安徽省电力有限公司电力科学研究院 | 电站锅炉尾部烟道氧量场标定方法 |
| CN117383518B (zh) * | 2023-12-11 | 2024-03-08 | 福建德尔科技股份有限公司 | 一种基于温度控制的硫磺与氟气反应系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103216846A (zh) * | 2013-04-24 | 2013-07-24 | 广东电网公司电力科学研究院 | 煤粉锅炉改变燃煤种类后的烟气浓度调整方法及系统 |
| CN103235842A (zh) * | 2013-03-29 | 2013-08-07 | 广东电网公司电力科学研究院 | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 |
| CN103955599A (zh) * | 2014-03-28 | 2014-07-30 | 广东电网公司电力科学研究院 | 锅炉改造后不同周界风量下燃烧特性的获取方法 |
| CN105891414A (zh) * | 2014-12-26 | 2016-08-24 | 浙江浙能技术研究院有限公司 | 基于网格布置法的无采样泵烟气浓度测量方法及测量装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7615200B2 (en) * | 2000-12-01 | 2009-11-10 | Fuel Tech, Inc. | Selective catalytic reduction of NOx enabled by urea decomposition in heat-exchanger bypass |
-
2016
- 2016-11-08 CN CN201610981223.XA patent/CN106548032B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103235842A (zh) * | 2013-03-29 | 2013-08-07 | 广东电网公司电力科学研究院 | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 |
| CN103216846A (zh) * | 2013-04-24 | 2013-07-24 | 广东电网公司电力科学研究院 | 煤粉锅炉改变燃煤种类后的烟气浓度调整方法及系统 |
| CN103955599A (zh) * | 2014-03-28 | 2014-07-30 | 广东电网公司电力科学研究院 | 锅炉改造后不同周界风量下燃烧特性的获取方法 |
| CN105891414A (zh) * | 2014-12-26 | 2016-08-24 | 浙江浙能技术研究院有限公司 | 基于网格布置法的无采样泵烟气浓度测量方法及测量装置 |
Non-Patent Citations (1)
| Title |
|---|
| 300MW电站锅炉富氧燃烧的数值模拟与分析;王建强;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20081115(第11期);论文摘要、第21-26、31、38-45页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106548032A (zh) | 2017-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106548032B (zh) | 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法 | |
| Guo et al. | Experimental and numerical investigations on oxy-coal combustion in a 35 MW large pilot boiler | |
| Echi et al. | CFD simulation and optimization of industrial boiler | |
| CN103968371B (zh) | 电力燃烧锅炉及基于数值模拟技术的分离燃尽风调节方法 | |
| Yin et al. | Characterizing and modeling of an 88 MW grate-fired boiler burning wheat straw: Experience and lessons | |
| Askarova et al. | Investigation of turbulence characteristics of burning process of the solid fuel in BKZ 420 combustion chamber | |
| Zhang et al. | Analysis and neural network prediction of combustion stability for industrial gases | |
| CN105930585A (zh) | 基于CFD的Shell气化炉流场和温度场的模拟方法 | |
| CN103148506A (zh) | 前后对冲旋流燃烧煤粉锅炉的二次风配风方法及系统 | |
| CN110263467A (zh) | 一种基于场协同原理的四角切圆锅炉数值分析方法 | |
| Guo et al. | Experimental and numerical investigations on heat transfer characteristics of a 35MW oxy-fuel combustion boiler | |
| CN103955599B (zh) | 锅炉改造后不同周界风量下燃烧特性的获取方法 | |
| Xu et al. | Simulation of the gas temperature deviation in large-scale tangential coal fired utility boilers | |
| CN106838979B (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 | |
| CN103148508A (zh) | 前后对冲旋流燃烧煤粉锅炉的二次风配风方法及系统 | |
| Bianco et al. | CFD analysis and design optimization of an air manifold for a biomass boiler | |
| Cui et al. | Scale-up prediction of supercritical CO2 circulating fluidized bed boiler based on adaptive PSO-SVM | |
| CN103324991B (zh) | 煤粉锅炉NOx排放预测的方法与系统 | |
| Holkar et al. | CFD Anlysis of Pulverised-Coal Combustion of Burner Used In Furnace with Different Radiation Models | |
| Ferrín et al. | Distribution of the coal flow in the mill-duct system of the As Pontes Power Plant using CFD modeling | |
| Hart et al. | Aerodynamics of an isolated slot-burner from a tangentially-fired boiler | |
| Shang et al. | Simulation of gas-particle turbulent combustion in a pulverized coal-fired swirl combustor | |
| CN106021916B (zh) | 一种适用于超超临界锅炉NOx排放量分析的计算方法 | |
| Ariwibowo et al. | Investigation on combustion characteristics of nonpremixed combustion in an 85 MWe wall firing boiler under various air-fuel ratio using CFD |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190703 Address after: 065000 No. 1 Building, No. 20 Qingxiang Road, Guangyang Economic Development Zone, Langfang City, Hebei Province Patentee after: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Address before: 065000 No. 10 Changxiang Road, Guangyang Economic Development Zone, Langfang City, Hebei Province Patentee before: WATSON MEASUREMENT & CONTROL TECHNOLOGY (HEBEI) Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Reconstruction method of distribution field based on grid detection of flue gas CO and O2concentration in boiler flue Effective date of registration: 20221223 Granted publication date: 20190405 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2022980028822 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20231212 Granted publication date: 20190405 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2022980028822 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Reconstruction method of distribution field for detecting CO and O2concentrations in flue gas based on boiler flue grid detection Effective date of registration: 20231215 Granted publication date: 20190405 Pledgee: China Everbright Bank Limited by Share Ltd. Langfang branch Pledgor: WATSON ENERGY TECHNOLOGY (LANGFANG) Co.,Ltd. Registration number: Y2023110000529 |