CN103234219A - 煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 - Google Patents
煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 Download PDFInfo
- Publication number
- CN103234219A CN103234219A CN2013101454986A CN201310145498A CN103234219A CN 103234219 A CN103234219 A CN 103234219A CN 2013101454986 A CN2013101454986 A CN 2013101454986A CN 201310145498 A CN201310145498 A CN 201310145498A CN 103234219 A CN103234219 A CN 103234219A
- Authority
- CN
- China
- Prior art keywords
- coal
- boiler
- air quantity
- fired
- pulverized
- 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.)
- Granted
Links
Images
Landscapes
- Regulation And Control Of Combustion (AREA)
Abstract
本发明公开了一种煤粉锅炉改变燃煤种类后的周界风量调整方法,包括:建立所述锅炉的网格化结构模型;建立煤粉燃烧所形成的各个理化过程的数学模型;对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。此外,还公开了一种煤粉锅炉改变燃煤种类后的周界风量调整系统。本发明可提高煤粉锅炉改变燃煤种类后周界风量控制的准确率和安全性,克服了采用人工控制所容易引起的炉内燃烧不稳定和燃烧效率低下的问题,同时大大降低了燃烧过程中产生的污染物。
Description
技术领域
本发明涉及煤粉锅炉技术领域,尤其涉及一种煤粉锅炉改变燃煤种类后的周界风量调整方法和一种煤粉锅炉改变燃煤种类后的周界风量调整系统。
背景技术
随着用电量越来越大,对煤的需求量也越来越多,电煤的生产和运输与用电需求的矛盾越来越突出,已经成为制约我国当前经济发展的瓶颈。煤炭供应持续紧张的同时,煤价也不断上涨,造成发电成本不断增加。为降低发电成本,电厂纷纷开始寻找低价位煤种,这就造成多数电厂很难保证长期燃用单一的设计煤种。另外,煤炭资源分布不平衡,使得北煤南运,西煤东运成为常态。因运力的严重不足等因素,导致很多电厂必须对锅炉进行掺烧甚至改烧。
在这种情况下,电厂已经开始燃用非设计煤种,在锅炉掺烧非设计煤种。但由于客观条件的限制,掺烧的依据基本采用的是煤中各成分简单的算术加和法,该方法对部分煤质的燃烧特性、灰熔融特性接近的煤种是基本可行的,但对煤的燃烧特性或灰特性差异较大时,掺烧存在锅炉运行不安全和不经济的问题。与此同时,长期燃烧、掺烧劣质煤种带来的温室气体和有害气体排放也是引起广泛关注的一大问题。煤炭消费又是导致SO2(二氧化硫)排放的主要原因,80%以上的SO2排放来自于燃煤排放,大量的SO2排放已经产生大面积的酸雨区。能源消费也是造成温室气体排放的主要原因,减排温室气体的压力已经变得十分紧迫。因此,对于长期燃用非设计煤种的电厂,为了增强大型火电厂对煤种的适应性,降低发电成本,提高电厂的经济性和安全性,众多电厂不得不开始对锅炉进行系统改造来满足要求。例如天津大港发电厂对300MW机组进行了改造,改烧烟煤取代原有的贫煤设计煤种;又如云浮发电厂针对来煤煤质变化大所带来的一系列问题,对锅炉系统的燃烧器区域和制粉系统进行了改造。
但对于电厂中广泛使用的四角切圆煤粉锅炉而言,改烧非设计煤种后(例如贫煤锅炉改烧烟煤),在不同的周界风量情况对锅炉燃烧特性的影响方面还缺乏成熟的工程应用,例如锅炉改烧其它燃煤种类后煤质特性变化对锅炉燃尽、结渣和积灰等方面的影响,目前大多是依靠人工经验进行控制,严重影响了锅炉运行的安全性和经济性。
发明内容
基于此,本发明提供了一种煤粉锅炉改变燃煤种类后的周界风量调整方法和一种煤粉锅炉改变燃煤种类后的周界风量调整系统。
一种煤粉锅炉改变燃煤种类后的周界风量调整方法,包括以下步骤:
根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
与一般技术相比,本发明煤粉锅炉改变燃煤种类后的周界风量调整方法建立与煤粉燃烧所形成的各个理化过程相对应的数学模型,采用数值模拟的方法获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系,并根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。本发明在数值模拟过程中耗时短,计算成本低,实现过程中不需电厂现场的调控配合,并且模拟结果具有很好的可视性。本发明可提高煤粉锅炉改变燃煤种类后周界风量控制的准确率和安全性,克服了采用人工控制所容易引起的炉内燃烧不稳定和燃烧效率低下的问题,同时大大降低了燃烧过程中产生的污染物。
一种煤粉锅炉改变燃煤种类后的周界风量调整系统,包括结构模型建立模块、数学模型建立模块、模拟模块和调整模块;
所述结构模型建立模块,用于根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
所述数学模型建立模块,用于根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
所述模拟模块,用于根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
所述调整模块,用于根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
与一般技术相比,本发明煤粉锅炉改变燃煤种类后的周界风量调整系统建立与煤粉燃烧所形成的各个理化过程相对应的数学模型,采用数值模拟的方法获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系,并根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。本发明在数值模拟过程中耗时短,计算成本低,实现过程中不需电厂现场的调控配合,并且模拟结果具有很好的可视性。本发明可提高煤粉锅炉改变燃煤种类后周界风量控制的准确率和安全性,克服了采用人工控制所容易引起的炉内燃烧不稳定和燃烧效率低下的问题,同时大大降低了燃烧过程中产生的污染物。
附图说明
图1为本发明煤粉锅炉改变燃煤种类后的周界风量调整方法的流程示意图;
图2为炉膛结构的网格化示意图;
图3为不同周界风量时煤粉燃尽率的对应关系示意图;
图4为不同周界风量时飞灰含碳量的对应关系示意图;
图5为不同周界风量时NOx排放量的对应关系示意图;
图6为不同周界风量时下炉膛出口温度的对应关系示意图;
图7为不同周界风量时沿炉膛高度的平均温度分布示意图;
图8为不同周界风量时沿炉膛高度的氧量分布示意图;
图9为不同周界风量时沿炉膛高度的CO生成量分布示意图;
图10为不同周界风量时沿炉膛高度的NOx生成量分布示意图;
图11为本发明煤粉锅炉改变燃煤种类后的周界风量调整系统的结构示意图。
具体实施方式
为更进一步阐述本发明所采取的技术手段及取得的效果,下面结合附图及较佳实施例,对本发明的技术方案,进行清楚和完整的描述。
请参阅图1,为本发明煤粉锅炉改变燃煤种类后的周界风量调整方法的流程示意图。本发明煤粉锅炉改变燃煤种类后的周界风量调整方法,包括以下步骤:
S101根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
根据燃烧器出口的流动特性,作为锅炉燃烧器入口边界条件,构建求解区域并进行网格划分,在炉膛出口增加4m的水平和4m的收缩段,以在出口处获得充分发展流。为了获得高质量的计算网格,采用分区网格划分方法,采用结构化六面体网格,为了准确模拟燃烧器出口空气动力场,燃烧器入口区域需要进行网格局部加密,总的网格数为200万。如图2所示,为炉膛结构的网格化示意图。
作为其中一个实施例,在所述建立所述锅炉的网格化结构模型的步骤中,所述设计参数包括蒸发量、主蒸汽压力、主蒸汽温度、再热蒸汽流量、再热蒸汽进口压力、再热蒸汽出口压力、再热蒸汽进口温度、再热蒸汽出口温度、给水温度、排烟温度、热风温度、锅炉效率、计算燃料消耗量、实际燃料消耗量以及燃烧器设计参数。
建立网格化模型时,必要的设计参数需全部包括,这样可确保能够建立有效的模型。另外,设计参数越多,则建立的模型越准确。表1是一个锅炉主要设计参数的实施例。
表1锅炉主要设计参数
作为其中一个实施例,所述燃烧器设计参数包括一次风、二次风、三次风、一次风周界风、三次风周界风和炉膛漏风。
建立网格化模型时,燃烧器设计参数越多越好,这样可确保能够建立有效和准确的模型。表2为一个燃烧器设计参数的实施例。
表2燃烧器设计参数
FLUENT提供了分离式和耦合式两类求解器,而耦合式又分为隐式和显式两种。
分离式求解器(segregated solver)是顺序地、逐一地求解各方程(关于u、v、w、p和T的方程)。也就是先在全部网格上解出一个方程(如u动量方程)后,再解另外一个方程(如v动量方程)。由于控制方程是非线性的,且相互之间是耦合的,因此,在得到收敛解之前,要经过多轮迭代。每一轮迭代由如下步骤组成:
根据当前的解的结果,更新所有流动变量。如果计算刚刚开始,则用初始值来更新;
按顺序分别求解u、v和w动量方程,得到速度场。在计算时,压力和单元界面的质量流量使用当前的已知值;
因上一步骤得到的速度很可能不满足连续方程,因此,用连续方程和线性化的动量方程构造一个Poisson型的压力修正方程,然后求解该压力修正方程,得到压力场与速度场的修正值;
利用新得到的速度场与压力场,求解其他标量(如温度、湍动能和组分等)的控制方程;
对于包含离散相的模拟,当内部存在相间耦合时,根据离散相的轨迹计算结果更新连续相的源项;
检查方程组是否收敛。若不收敛,回到第一步,重复进行。
耦合式求解器(coupled solver)是同时求解连续方程、动量方程、能量方程及组分疏运方程的耦合方程,然后,再逐一地求解湍流等标量方程。由于控制方程是非线性的,且相互之间是耦合的,因此,在得到收敛解之前,要经过2迭代。每一轮迭代由如下步骤组成:
根据当前的解的结果,更新所有流动变量。如果计算刚刚开始,则用初始值来更新;
同时求解连续方程、动量方程、能量方程及组分输运方程的耦合方程组(后两个方程视需要进行求解);
根据需要,逐一地求解湍流、辐射等标量方程。注意在求解之前,方程中用到的有关变量要用前面得到的结果更新;
对于包含离散相的模拟,当内部存在相间耦合时,根据离散相的轨迹计算结果更新连续相的源项;
检查方程组是否收敛。若不收敛,回到第一步,重复进行。
两种求解器都适用于从不可压到高速可压的很大范围的流动,但总的来讲,当计算高速可压流动时,耦合式求解器更具有优势。耦合式求解器常常可以很快收敛,但所需要的内存大约是分离式求解器的1.5到2倍。
此外,分离式求解器中提供的几个物理模型,在耦合式求解器中是没有的。这些物理模型包括:流体体积模型(VOF)、多相混合模型、欧拉混合模型、PDF燃烧模型、预混合燃烧模型、指定质量流量的周期流动模型、周期性热传导模型和壳传导模型等。
由于本发明将运用PDF燃烧模型模拟锅炉燃烧,且计算机硬件配置有限,故选择分离式求解器。
FLUENT中提供了三种算法:SIMPLE算法、SIMPLEC算法和PISO算法。SIMPLEC和PISO算法具有较快的收敛速度,相对于SIMPLE算法的计算时间减少30%~50%。对于瞬态问题,PISO算法具有明显的优势,而对于稳态问题,可选择SIMPLE系列算法。本发明可选用SIMPLE算法。
SIMPLE算法可描述如下:对于给定的压力场(它可以是假定的值,或是上一次迭代计算所得到的结果),求解离散形式的动量方程,得出速度场,因此,必须对给定的压力场加以修正。修正的原则是:与修正后的压力场相对应的速度场能满足这一迭代层次上的连续方程。据此原则,我们把由动量方程的离散形式所规定的压力与速度的关系代入连续方程的离散形式,从而得到压力修正方程,由压力修正方程得出压力修正值。接着,根据修正后的压力场,求得新的速度场。然后检查速度场是否收敛。若不收敛,用修正后的压力值作为给定的压力场,开始下一层次的计算。如此反复,直到获得收敛的解。
S102根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
湍流应力方程模型:
湍流流动的基本控制方程为:
连续方程:
动量方程(Navier-Stokes方程):
其他变量的输运方程:
可见,方程组中多了六个应力(3个正应力和3个切应力),即:
针对这六个应力的处理方式,得到了各种湍流数值模拟的方法:
本发明选用Realizable k-ε模型进行模拟。标准k-ε模型对时均应变率特别大的情形,有可能导致负的正应力。而Realizable k-ε模型主要变化是:
湍流粘度计算公式发生了变化,引入了与旋转和曲率有关的内容:
ε方程发生了很大变化,方程中的产生项不再包含有k方程中的产生项Gk,这样,现在的形式更好地表示了光谱的能量转换。
ε方程中的倒数第二项不具有任何奇异性,即使k值很小或为零,分母也不会为零。这与标准k-ε模型和RNG k-ε模型有很大区别。
Realizable k-ε模型已被有效得应用于各种不同类型的流动模拟,包括旋转均匀剪切流、包含有射流和混合流的自由流动、管道内流动、边界层流动,以及带有分离的流动等。
辐射模型:
FLUENT中可以用5种模型来计算辐射换热问题,这5中模型分别是离散换热辐射模型(DTRM)、P-1辐射模型、Rosseland辐射模型、表面辐射模型(S2S)和离散坐标(DO)辐射模型。
针对本发明研究的问题为锅炉的燃烧,故选用P-1辐射模型。
P-1辐射模型,辐射换热方程是一个计算相对较小的扩散方程,同时模型中包含了散射效应,在燃烧等光学厚度很大的计算问题中,P-1模型的计算效果都比较好。
组分输运与化学反应模型:
FLUENT中提供了四种模拟方法:通用有限速率模型;非预混合燃烧模型;预混合燃烧模型;部分预混合燃烧模型。
模型选择的大致方针如下:
通用有限速度模型主要用于:化学组分混合、输运和反应的问题;壁面或者粒子表面反应的问题(如化学蒸气沉积);
非预混合燃烧模型主要用于:包括湍流扩散火焰的反应系统,这个系统接近化学平衡,其中的氧化物和燃料以两个或者三个流道分别流入所要计算的区域;
预混合燃烧模型主要用于:单一、完全预混合反应物流动;
部分预混合燃烧模型主要用于:区域内具有变化等值比率的预混合火焰的情况;
本发明选择非预混合燃烧模型,用混合分数-概率密度函数模拟气相湍流燃烧。
NOx生成模型:
NOx(氮氧化物)生成类型有三种:热力型、快速型、燃料型。生成机理已于第一章介绍。基于本发明研究对象为煤粉燃烧的数值模拟,故不考虑快速型NOx(其主要由CH燃料燃烧生成),但考虑了NOx的再燃效应。此外,模型中不考虑N2O的生成。
S103根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
作为其中一个实施例,所述锅炉改变后的燃煤种类为烟煤。
目前,烟煤是被广泛采用的一种非设计煤种,锅炉改烧烟煤降低煤粉的消耗率。
作为其中一个实施例,所述锅炉燃烧性能指标包括煤粉燃尽率、飞灰含碳量、氮氧化物排放量和下炉膛出口温度。
在模拟过程中,分析的燃烧性能指标越多,在对锅炉运行进行控制的准确率越高,安全性也越好,对后续的工程实践越有指导意义。
在模拟之前,需要对锅炉设计煤种和实际燃烧煤种进行一个对比分析。表3为一个设计煤常规分析数据的实施例。
表3设计煤常规分析数据
表4为一个实际来煤与设计煤种的特性进行对比分析的实施例。
表4实际来煤与设计煤种的特性比较
请参阅图3,为不同周界风量时煤粉燃尽率的对应关系示意图。
请参阅图4,为不同周界风量时飞灰含碳量的对应关系示意图。
周界风量为119(10.8%)、159(14.5%)和200t/h(18.2%)时,煤粉燃尽率分别为99.67%、99.85%和99.80%,飞灰含碳量分别为0.513%、0.311%和0.390%。标准工况下,即周界风量为159t/h时,煤粉燃尽率最大,飞灰含碳量则最小。当周界风量由119t/h增加到159t/h时,煤粉燃尽率增加了0.18%,比较可观,再当增加到200t/h时,燃尽率在表况下降低了0.05%,影响程度不是十分明显,说明周界风量在标况下,即是159t/h是锅炉煤粉燃烧效率最高。而飞灰含碳量的区别也相当小,说明变周界风对于锅炉整体燃烧效率影响程度不是十分巨大。
请参阅图5,为不同周界风量时NOx排放量的对应关系示意图。
周界风量为119、159和200t/h时,NOx排放量分别为203.42、211.15和211.46mg/m3(6%O2)。几种变周界风工况下的NOx排放量并无显著差别。
请参阅图6,为不同周界风量时下炉膛出口温度的对应关系示意图。
周界风量为119、159和200t/h时,下炉膛出口烟气温度分别为1336、1344和1346K,基本相同。这是因为改变周界风量但是燃烧器区域的总风量是不变的,周界风量的改变引起的波动又由上下层补充的二次风弥补,最终周界风的改变形成的影响相当小。
煤粉燃尽率和飞灰含碳量随周界风量的改变会稍有影响,但变化都很小;NOx排放量随周界风量的改变基本没有变化;下炉膛出口烟气温度随燃烧器出力变化基本没有变化。
作为其中一个实施例,在所述建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型的步骤之后,包括以下步骤:
根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取在各种周界风量情况下炉膛平均温度、氧量、一氧化碳生成量和氮氧化物生成量与炉膛高度之间的对应关系。
获取炉膛平均温度、氧量、一氧化碳生成量和氮氧化物生成量与炉膛高度之间的对应关系有利于在锅炉运行过程中更可靠地进行控制。
请参阅图7,为不同周界风量时沿炉膛高度的平均温度分布示意图。
沿炉膛高度的平均温度分布趋势与变氧量下的分布趋势是一致的,不同周界风量时的分布也基本是一致的。这表明,变周界风量对炉膛平均烟气温度基本没有影响。
请参阅图8,为不同周界风量时沿炉膛高度的氧量分布示意图。
沿炉膛高度的平均氧量分布趋势与变氧量下的分布趋势是一致的,不同周界风量时的分布也基本是一致的。这表明,变周界风量对炉膛平均氧量基本没有影响。
请参阅图9,为不同周界风量时沿炉膛高度的CO生成量分布示意图。
沿炉膛高度的平均CO质量分数分布趋势与变氧量下的分布趋势是一致的,不同周界风量时的平均CO质量分数分布也基本一致,只是在冷灰斗部分由区别,这是由于最下层风量的不同造成煤粉掉入冷灰斗的比例不同,造成小部分煤粉颗粒在冷灰斗底部不完全燃烧,CO含量也有些许不同。
请参阅图10,为不同周界风量时沿炉膛高度的NOx生成量分布示意图。
沿炉膛高度的平均NOx生成量分布趋势与变氧量下的分布趋势是一致的,最终变周界风量对炉膛出口NOx排放量基本没有影响。
S104根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
可采用各种方法对锅炉的周界风量进行控制,例如可采用计算机控制,也可人工进行控制。
作为其中一个实施例,在所述对所述锅炉的周界风量进行调整的步骤之后,包括以下步骤:
在所述锅炉的运行过程中,当锅炉燃烧性能指标超过预设告警指标时,则发出告警信号。
当锅炉燃烧性能指标超过预设告警指标时,进行报警,可更大程度地确保锅炉控制的安全性。
与一般技术相比,本发明煤粉锅炉改变燃煤种类后的周界风量调整方法建立与煤粉燃烧所形成的各个理化过程相对应的数学模型,采用数值模拟的方法获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系,并根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。本发明在数值模拟过程中耗时短,计算成本低,实现过程中不需电厂现场的调控配合,并且模拟结果具有很好的可视性。本发明可提高煤粉锅炉改变燃煤种类后周界风量控制的准确率和安全性,克服了采用人工控制所容易引起的炉内燃烧不稳定和燃烧效率低下的问题,同时大大降低了燃烧过程中产生的污染物。
请参阅图11,为本发明煤粉锅炉改变燃煤种类后的周界风量调整系统的结构示意图。
本发明煤粉锅炉改变燃煤种类后的周界风量调整系统,包括结构模型建立模块201、数学模型建立模块202、模拟模块203和调整模块204;
所述结构模型建立模块201,用于根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
所述数学模型建立模块202,用于根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
所述模拟模块203,用于根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
所述调整模块204,用于根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
作为其中一个实施例,所述锅炉燃烧性能指标包括煤粉燃尽率、飞灰含碳量、氮氧化物排放量和下炉膛出口温度。
在模拟过程中,分析的燃烧性能指标越多,在对锅炉运行进行控制的准确率越高,安全性也越好,对后续的工程实践越有指导意义。
作为其中一个实施例,还包括告警模块;
所述告警模块,用于在所述锅炉的运行过程中,当锅炉燃烧性能指标超过预设告警指标时,发出告警信号。
当锅炉燃烧性能指标超过预设告警指标时,进行报警,可更大程度地确保锅炉控制的安全性。
与一般技术相比,本发明煤粉锅炉改变燃煤种类后的周界风量调整系统建立与煤粉燃烧所形成的各个理化过程相对应的数学模型,采用数值模拟的方法获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系,并根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。本发明在数值模拟过程中耗时短,计算成本低,实现过程中不需电厂现场的调控配合,并且模拟结果具有很好的可视性。本发明可提高煤粉锅炉改变燃煤种类后周界风量控制的准确率和安全性,克服了采用人工控制所容易引起的炉内燃烧不稳定和燃烧效率低下的问题,同时大大降低了燃烧过程中产生的污染物。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,包括以下步骤:
根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
2.根据权利要求1所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,在所述建立所述锅炉的网格化结构模型的步骤中,所述设计参数包括蒸发量、主蒸汽压力、主蒸汽温度、再热蒸汽流量、再热蒸汽进口压力、再热蒸汽出口压力、再热蒸汽进口温度、再热蒸汽出口温度、给水温度、排烟温度、热风温度、锅炉效率、计算燃料消耗量、实际燃料消耗量以及燃烧器设计参数。
3.根据权利要求2所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,所述燃烧器设计参数包括一次风、二次风、三次风、一次风周界风、三次风周界风和炉膛漏风。
4.根据权利要求1所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,所述锅炉改变后的燃煤种类为烟煤。
5.根据权利要求1所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,所述锅炉燃烧性能指标包括煤粉燃尽率、飞灰含碳量、氮氧化物排放量和下炉膛出口温度。
6.根据权利要求1所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,在所述对所述锅炉的周界风量进行调整的步骤之后,包括以下步骤:
在所述锅炉的运行过程中,当锅炉燃烧性能指标超过预设告警指标时,发出告警信号。
7.根据权利要求1所述的煤粉锅炉改变燃煤种类后的周界风量调整方法,其特征在于,在所述建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型的步骤之后,包括以下步骤:
根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取在各种周界风量情况下炉膛平均温度、氧量、一氧化碳生成量和氮氧化物生成量与炉膛高度之间的对应关系。
8.一种煤粉锅炉改变燃煤种类后的周界风量调整系统,其特征在于,包括结构模型建立模块、数学模型建立模块、模拟模块和调整模块;
所述结构模型建立模块,用于根据四角切圆煤粉锅炉的设计参数,按照直流燃烧器、冷灰斗区域、燃烧器区域、燃烧器上方区域和屏式过热器区域的结构,建立所述锅炉的网格化结构模型;
所述数学模型建立模块,用于根据所述网格化结构模型,建立煤粉燃烧所形成的气相湍流流动过程、气固两项流动过程、辐射传热过程和氮氧化物生成过程的数学模型;
所述模拟模块,用于根据建立的所述数学模型,对所述锅炉改变燃煤种类后的煤粉燃烧过程进行模拟,获取所述锅炉的各种周界风量情况与锅炉燃烧性能指标之间的对应关系;
所述调整模块,用于根据获取的所述对应关系,对所述锅炉的周界风量进行调整,使所述锅炉满足预设的燃烧性能指标。
9.根据权利要求8所述的煤粉锅炉改变燃煤种类后的周界风量调整系统,其特征在于,所述锅炉燃烧性能指标包括煤粉燃尽率、飞灰含碳量、氮氧化物排放量和下炉膛出口温度。
10.根据权利要求8所述的煤粉锅炉改变燃煤种类后的周界风量调整系统,其特征在于,还包括告警模块;
所述告警模块,用于在所述锅炉的运行过程中,当锅炉燃烧性能指标超过预设告警指标时,发出告警信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310145498.6A CN103234219B (zh) | 2013-04-24 | 2013-04-24 | 煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310145498.6A CN103234219B (zh) | 2013-04-24 | 2013-04-24 | 煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234219A true CN103234219A (zh) | 2013-08-07 |
CN103234219B CN103234219B (zh) | 2016-02-24 |
Family
ID=48882271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310145498.6A Active CN103234219B (zh) | 2013-04-24 | 2013-04-24 | 煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234219B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103712235A (zh) * | 2013-12-26 | 2014-04-09 | 广东电网公司电力科学研究院 | 一种一次风机的控制方法和装置 |
CN103953942A (zh) * | 2014-05-17 | 2014-07-30 | 国家电网公司 | 一种煤粉锅炉一次风速调平方法 |
CN103955599B (zh) * | 2014-03-28 | 2017-02-01 | 广东电网公司电力科学研究院 | 锅炉改造后不同周界风量下燃烧特性的获取方法 |
CN106500955A (zh) * | 2016-10-20 | 2017-03-15 | 广东电网有限责任公司电力科学研究院 | 电站锅炉煤粉在变负荷下燃烧特性模拟方法及装置 |
CN106765281A (zh) * | 2016-12-30 | 2017-05-31 | 上海外高桥第二发电有限责任公司 | 一种900mw电站锅炉燃烧优化配风调整方法 |
CN110410819A (zh) * | 2019-08-05 | 2019-11-05 | 东南大学 | 燃煤锅炉燃烧动力场快速计算和实时显示系统及计算方法 |
CN110906361A (zh) * | 2019-12-12 | 2020-03-24 | 广州珠江电力有限公司 | 一种基于低氮燃烧的优化燃烧控制方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH094839A (ja) * | 1995-06-16 | 1997-01-10 | Babcock Hitachi Kk | 微粉炭ボイラの制御装置 |
JPH1026343A (ja) * | 1996-07-09 | 1998-01-27 | Ishikawajima Harima Heavy Ind Co Ltd | 石炭焚ボイラのミル残炭パージ時における出炭量算出方法及び装置 |
CN101996260A (zh) * | 2009-08-14 | 2011-03-30 | 华东电力试验研究院有限公司 | 数值模拟预测炉内结渣特性的方法 |
CN102799775A (zh) * | 2012-07-06 | 2012-11-28 | 广东电网公司电力科学研究院 | 利用数值模拟辅助超超临界锅炉屏式过热器管壁温度的计算方法 |
-
2013
- 2013-04-24 CN CN201310145498.6A patent/CN103234219B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH094839A (ja) * | 1995-06-16 | 1997-01-10 | Babcock Hitachi Kk | 微粉炭ボイラの制御装置 |
JPH1026343A (ja) * | 1996-07-09 | 1998-01-27 | Ishikawajima Harima Heavy Ind Co Ltd | 石炭焚ボイラのミル残炭パージ時における出炭量算出方法及び装置 |
CN101996260A (zh) * | 2009-08-14 | 2011-03-30 | 华东电力试验研究院有限公司 | 数值模拟预测炉内结渣特性的方法 |
CN102799775A (zh) * | 2012-07-06 | 2012-11-28 | 广东电网公司电力科学研究院 | 利用数值模拟辅助超超临界锅炉屏式过热器管壁温度的计算方法 |
Non-Patent Citations (1)
Title |
---|
刘丽萍: "四角切圆煤粉炉炉内燃烧及配风的数值模拟", 《中国优秀硕士学位论文全文数据库》, 15 May 2009 (2009-05-15) * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103712235A (zh) * | 2013-12-26 | 2014-04-09 | 广东电网公司电力科学研究院 | 一种一次风机的控制方法和装置 |
CN103712235B (zh) * | 2013-12-26 | 2016-02-24 | 广东电网公司电力科学研究院 | 一种一次风机的控制方法和装置 |
CN103955599B (zh) * | 2014-03-28 | 2017-02-01 | 广东电网公司电力科学研究院 | 锅炉改造后不同周界风量下燃烧特性的获取方法 |
CN103953942A (zh) * | 2014-05-17 | 2014-07-30 | 国家电网公司 | 一种煤粉锅炉一次风速调平方法 |
CN103953942B (zh) * | 2014-05-17 | 2016-08-24 | 国家电网公司 | 一种煤粉锅炉一次风速调平方法 |
CN106500955A (zh) * | 2016-10-20 | 2017-03-15 | 广东电网有限责任公司电力科学研究院 | 电站锅炉煤粉在变负荷下燃烧特性模拟方法及装置 |
CN106765281A (zh) * | 2016-12-30 | 2017-05-31 | 上海外高桥第二发电有限责任公司 | 一种900mw电站锅炉燃烧优化配风调整方法 |
CN106765281B (zh) * | 2016-12-30 | 2018-12-11 | 上海外高桥第二发电有限责任公司 | 一种900mw电站锅炉燃烧优化配风调整方法 |
CN110410819A (zh) * | 2019-08-05 | 2019-11-05 | 东南大学 | 燃煤锅炉燃烧动力场快速计算和实时显示系统及计算方法 |
CN110906361A (zh) * | 2019-12-12 | 2020-03-24 | 广州珠江电力有限公司 | 一种基于低氮燃烧的优化燃烧控制方法及系统 |
CN110906361B (zh) * | 2019-12-12 | 2021-09-24 | 广州珠江电力有限公司 | 一种基于低氮燃烧的优化燃烧控制方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103234219B (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103225819B (zh) | 煤粉锅炉改变燃煤种类后的氧量调整方法及系统 | |
CN103234219B (zh) | 煤粉锅炉改变燃煤种类后的周界风量调整方法及系统 | |
Al-Abbas et al. | CFD modelling of air-fired and oxy-fuel combustion in a large-scale furnace at Loy Yang A brown coal power station | |
CN103364212A (zh) | 煤粉锅炉改烧烟煤变负荷燃烧特性测量方法与系统 | |
CN103235842B (zh) | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 | |
CN103244953B (zh) | 煤粉锅炉改变燃煤种类后的燃烧器出力调整方法及系统 | |
CN103216846B (zh) | 煤粉锅炉改变燃煤种类后的烟气浓度调整方法及系统 | |
Al-Abbas et al. | Numerical modelling of oxy-fuel combustion in a full-scale tangentially-fired pulverised coal boiler | |
Jovanović et al. | Mathematical modelling of swirl oxy-fuel burner flame characteristics | |
CN103324991B (zh) | 煤粉锅炉NOx排放预测的方法与系统 | |
CN103323573A (zh) | 煤粉锅炉改烧烟煤变磨煤机组合燃烧特性测量方法与系统 | |
CN103335327B (zh) | 煤粉锅炉改烧烟煤变燃尽风量燃烧特性测量方法与系统 | |
Perrone et al. | Numerical simulation of MILD (moderate or intense low-oxygen dilution) combustion of coal in a furnace with different coal gun positions | |
Holkar et al. | CFD Anlysis of Pulverised-Coal Combustion of Burner Used In Furnace with Different Radiation Models | |
Jessee et al. | Multi-dimensional analysis of turbulent natural gas flames using detailed chemical kinetics | |
López et al. | A numerical analysis of the effect of atmospheric pressure on the performance of a heating system with a self-recuperative burner | |
Trinchenko et al. | Numerical Research of Nitrogen Oxides Formation for Justification of Modernization of P-49 Nazarovsky State District Power Plant Boiler on the Low-temperature Swirl Technology of Burning | |
Rahman et al. | Ammonia as a Hydrogen Vector: Validated Large Eddy Simulation of Ammonia Co-Firing in a Pilot-Scale Coal Combustor | |
Vega et al. | Geometrical parameter evaluation of a 0.5 MWth bench-scale oxy-combustion burner | |
Sampat | Automatic generation of chemical reactor networks for combustion simulations | |
Xie et al. | Numerical Simulation of Radiation Heat Transfer Characteristics of Oxy-fuel Combustion Furnace | |
Manickam et al. | Enriched oxygen combustion simulation for rotary kiln application | |
Li et al. | Effect of inclined angle of fuel jet on NO x emission in high temperature air combustion | |
Granados et al. | Numerical simulation of oxy-fuel combustion in a cement kiln | |
Nakate et al. | Reactive Turbulent Flow Model Of Anode Baking Furnace To Estimate NOx Through Zeldovich Mechanism |
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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Patentee after: ELECTRIC POWER RESEARCH INSTITUTE, GUANGDONG POWER GRID CO., LTD. Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Patentee before: Electrical Power Research Institute of Guangdong Power Grid Corporation |
|
CP01 | Change in the name or title of a patent holder |