CN102252779A - 基于烟气能量平衡的炉膛出口烟温优化测量方法 - Google Patents
基于烟气能量平衡的炉膛出口烟温优化测量方法 Download PDFInfo
- Publication number
- CN102252779A CN102252779A CN201110107881.3A CN201110107881A CN102252779A CN 102252779 A CN102252779 A CN 102252779A CN 201110107881 A CN201110107881 A CN 201110107881A CN 102252779 A CN102252779 A CN 102252779A
- Authority
- CN
- China
- Prior art keywords
- flue gas
- heat
- furnace outlet
- furnace
- gas temperature
- 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
Abstract
一种信息处理技术领域的基于烟气能量平衡的炉膛出口烟温优化测量方法,通过机组中和模型参数相关的实际运行数据,根据炉膛燃烧与辐射能量平衡,过热器烟气侧与工质侧能量平衡,再热器烟气侧与工质侧能量平衡以及尾烟气能量衡算,从烟气侧正反两个方向对炉膛烟温进行了优化测量,并分别测量结果的基础上复合计算最终烟温结果。在实际运行中可根据机组运行状况修正方法参数,使得该发明方法具有较好的通用性与鲁棒性,且对实际现场监测系统要求较低,适用于火力发电系统的运行状况检测与燃烧状况的优化。
Description
技术领域
本发明涉及的是一种信息处理技术领域的方法,具体是一种基于烟气能量平衡的炉膛出口烟温优化测量方法。
背景技术
锅炉作为火电机组的重要组成部分,随着机组容量的增加,其结构越来越复杂,所需检测和控制的参数也越来越多,给其控制和优化带来了新的课题。炉膛作为进入锅炉系统燃料的燃烧空间与锅炉系统辐射传热的主要受热面,其具有体积大、输入输出量多且内部工况复杂特点。炉膛热力计算作为锅炉热力计算中的重要组成部分,其主要结果炉膛出口烟温对整个锅炉的建模、热力计算以及效率评估具有重大意义。
20世纪70年代起,国内一些研究机构利用各种计算语言开发了许多具有一定实用性的热力计算系统,如1998年浙江大学与杭州锅炉厂合作开发的热力计算系统,应用在杭州锅炉厂、太原锅炉厂厂家;2002年浙江大学与嘉兴东方日立锅炉厂联合开发的针对燃机余热锅炉的热力计算系统。因受开发技术、方法、语言及计算模型自身结构复杂和离散度大因素的影响,在软件开发技术方面存在通用性、稳定性、可扩充性、可维护性差及计算精度低问题。由于煤种分布不同,我国锅炉使用分布也呈现一定的地域性,目前的锅炉热力计算系统对煤种的选择存在一定的局限性。
经对现有文献检索发现,《内蒙古农业大学学报》上的“电站锅炉炉膛的热力计算及仿真”在分析锅炉热力计算的现状和热力计算软件的基础上,针对锅炉炉膛的热力计算进行归纳、总结和简化,并提出了相关热力计算程序。该方法具备通用性好、计算简单特点。但该方法由于只采用单向能量平衡算法,缺乏有效地现场数据验证或者是多种方法结果比对,说服力不够。
华中科技大学煤燃烧国家重点实验室周怀春提出一种锅炉及加热炉温度场可视化及锅炉及加热炉温度场可视化及燃烧监控技术,基于热辐射图像处理温度检测原理和技术对大型锅炉炉内二维/三维温度进行检测。该方法具备检测精度高的特点,且能对炉内三维温度场进行全面地检测,但其相关技术因需要加入特殊的热辐射检测仪表,对锅炉炉膛的检测系统提出了较高的要求,暂时无法大范围推广。
发明内容
本发明针对现有技术存在的上述不足,提供一种基于烟气能量平衡的炉膛出口烟温优化测量方法,通过机组中和模型参数相关的实际运行数据,根据炉膛燃烧与辐射能量平衡,过热器烟气侧与工质侧能量平衡,再热器烟气侧与工质侧能量平衡以及尾烟气能量衡算,从烟气侧正反两个方向对炉膛烟温进行了优化测量,并分别测量结果的基础上复合计算最终烟温结果。在实际运行中可根据机组运行状况修正方法参数,使得该发明方法具有较好的通用性与鲁棒性,且对实际现场监测系统要求较低,适用于火力发电系统的运行状况检测与燃烧状况的优化。
本发明是通过以下技术方案实现的,本发明包括以下步骤:
步骤一:设置包括炉膛体积、换热面积以及各级换热器的换热面积、金属质量的锅炉炉膛及烟道特性,通过OPC通讯协议,建立OPC客户端向DCS自带的OPC服务器端发送数据请求,进行实时采样获取锅炉运行数据,对数据进行五点平滑处理。
所述的采样的时间间隔为5秒;
步骤二:在现场控制系统工程师站上建立科学计算平台与人机交互平台,即利用Matlab建立数据处理平台,基于Matlab自带的OPC客户端用于实现数据接收与处理。
步骤三:基于锅炉炉膛的燃烧放热平衡与水冷壁辐射换热平衡,计算出炉膛出口烟温T1;基于锅炉烟道的尾烟气能量、再热器能量平衡以及过热器能量平衡,利用烟气物性参数库,计算出炉膛出口烟温T2,具体步骤包括:
a)利用正向燃煤在炉膛中放热与传热平衡,计算出炉膛出口温度T1
b)利用反向烟气在各级换热器中的产热平衡,利用尾气能量加上各级换热器传递热量计算出炉膛出口烟气能量,并利用烟气物性反推炉膛出口烟温T2
步骤四:根据正向燃烧放热平衡与反向过热器再热器热交换平衡计算得到炉膛出口烟温T1与T2,采用线性叠加的方式,即T=a×T1+b×T2计算出最终炉膛烟温结果,其中系数a取0.3,系数b取0.7。
步骤五:根据实际燃料特性以及炉膛、烟道、换热器特性随着时间推移所发生的变化修正方法执行中炉膛换热系数等各组参数,即对优化测量方法进行参数滚动更新,具体是:
a)每天利用最新煤质在线分析值对烟气物性参数库进行更新;
b)每次大修以后对炉膛换热系数等各种参数进行重新设置。
与现有技术相比,本发明具备以下三个优点:一是无需新增硬件投入,只需利用现有检测系统所获得的数据以及部分火力发电系统的化验数据,避免了新增投资与技术改造;二是根据不同型号以及规模的火力发电系统,只需通过修正部分涉及到锅炉结构的参数,无需对模型结果推倒重来;三是采用了两种不同的烟气侧能量平衡,贯穿了整个烟气通道,覆盖了多种不确定性,并能进行互相校验。本发明的结果更有助于对锅炉整体运行情况的进行判断以及进一步优化。
附图说明
图1为烟气物性曲线图。
图2为炉膛出口烟温最终计算结果与负荷数据对照图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本实施例中的锅炉满足以下条件:该锅炉为负荷300MV、亚临界、自然循环、中间再热、双拱单炉膛、“W”火焰燃烧方式、固态排渣、全钢架悬吊式燃煤汽包炉。
本实施例包括以下步骤:
步骤一:设置锅炉炉膛以及烟道特性,并设定燃煤的成分与低位发热量,特性参数如表1:
表1锅炉物理特性参数
该锅炉DCS控制系统自带OPC服务端,该优化测量方法执行过程中可通过OPC客户端实时读取锅运行数据,采样周期为5秒。计算烟气物理性质,根据现场的实时煤质分析仪提供的数据,计算1Kg燃煤燃烧后的产物的各项特性,并建立烟气物性数据库。
且烟气的焓值,可以根据O2,CO2,H2O,N2气体的体积成分和烟气温度利用各种成分气体的焓值特性计算得到,典型平均成分烟气(rCO2=0.13,rH2O=0.11,rN2=0.76),其比热值如图1所示;工质水的焓值根据IAPWS-IF97公式获得。
步骤二:在工程师站上建立科学计算平台与人机交互界面。
步骤三:根据炉内传热相似理论计算方法可得炉膛出口温度与进入炉膛燃煤量之间的集总参数模型:
其中:Bj为计算燃料消耗量,为实测值;σ0为波尔兹曼常量;F1为炉膛换热壁面积,为2900m2;为考虑炉膛散热损失的保热系数,为0.85;I”1为炉膛出口截面上燃烧产物的焓Q1为基准燃料量送入炉膛内的有效热量;VCp为燃烧产物的平均比热容,查烟气物性参数库可得;Ta为燃烧产物的绝对燃烧温度,为2400℃;T”1为炉膛出口烟气温度;Ψ为水冷壁有效吸热系数为0.45,a1为炉膛黑度,为0.63。
根据后续烟气能量平衡,炉膛出口高温烟气能量满足下式:
高温烟气能量=烟气热交换中减少的能量+烟气尾气能量
基于烟道的绝热性较好这一假设可以认为烟气通过换热器时减少的能量与换热器中工质增加的能量相。即烟气侧能量平衡可改写为:
高温烟气能量=换热器中工质增加的能量+烟气尾气能量
由所研究燃煤锅炉结构可知换热器主要分为过热器和再热器,因而反向推导计算需主要考虑过热器与换热器换热量以及尾烟气能量。
(a)过热器热交换模型
根据过热器中工质质量平衡可得:
ΔM1=M1(k)-M1(k-1)=[ρ1(k)-ρ1(k-1)]V1
式中:k为当前时刻,k-1为前一时刻,Din1为汽包出口饱和蒸汽流量,即过热器系统输入质量,Dout1为高压缸入口过热蒸汽流量,即过热器系统输出质量,M1为过热器系统中工质总质量,V1为过热器容积,p1为过热器系统内工质的平均密度。
根据过热器中工质能量平衡可得:
ΔE1=E1(k)-E1(k-1)=[ρ1(k)H1(k)-ρ1(k)H1(k)]V+cmetal1Mmetal1ΔTmetal1
式中:k为当前时刻,k-1为前一时刻,Qin1为饱和蒸汽能量,即过热器系统输入能量,Qout1为过热蒸汽能量,即过热器系统输出能量,E1为过热器中工质总能量,Hin1和Hout1分别为饱和蒸汽焓值和过热蒸汽焓值,cmetal1,Mmetal1和ΔTmetal1为过热器管道金属壁的比热、质量和温度波动,ΔQgas1为烟气通过过热器传递给工质的能量,即烟气在过热器热交换中减少的能量。
(b)再热器热交换模型
根据再热器中工质能量平衡可得:
Din2=Dmain-Ds1-Ds2
ΔM2=M2(k)-M2(k-1)=[ρ2(k)-ρ2(k-1)]V2
式中:k为当前时刻,k-1为前一时刻,Din2为高压缸出口冷再热蒸汽流量,即再热器系统输入质量,Dout2为中压缸入口热再热蒸汽流量,即再热器系统输出质量,Ds1,Ds2分别为高压缸送入高加1#和高加2#的蒸汽流量,Dmain为主蒸汽流量,与式(5)中Dout1相,M2为再热器系统中工质总质量,V2为再热器容积,ρ2为过热器系统内工质的平均密度。
根据再热器能量平衡可得:
ΔE2=E2(k)-E2(k-1)=[ρ2(k)H2(k)-ρ2(k)H2(k)]V+cmetal2Mmetal2ΔTmetal2
式中:k为当前时刻,k-1为前一时刻,Qin2为冷再热蒸汽能量,即再热器系统输入能量,Qout2为热再热蒸汽能量,即再热器系统输出能量,E2为再热器中工质总能量,Hin2和Hout2分别为冷再热蒸汽焓值和热再热蒸汽焓值,cmetal2,Mmetal2和ΔTmetal2为过热器管道金ΔQgas2为烟气通过再热器传递给工质的能量,即烟气在再热器热交换中减少的能量。
(c)尾烟气能量衡算模型
烟气是燃料与氧气发生化学反应后的产物。烟气的主要成分有氮气,二氧化碳,氧气和水蒸气,还有少量的氢气、甲烷和其他氧化物(如二氧化硫、二氧化氮、一氧化碳)。
尾气能量计算方法采用基于烟气成分比例的方法。计算方法如下:
Qgas3=K(Qpy-Qlf)
式中:K为燃煤固体未完全燃烧修正值,Qpy为排烟能量,Qlf为冷空气能量。
因而,炉膛出口烟气的能量计算方法即为:
Q=Qgas1+Qgas2+Qgas3
基于炉膛出口烟气能力可通过烟气物性参数表反推炉膛出口烟气温度T2。
步骤四:根据基于两种软测量手段所获得炉膛出口烟温T1与T2,采用一定比例线性叠加的方式,即T=a×T1+b×T2计算出最终炉膛烟温结果。
本实施例中,设定a=0.3,b=0.7计算最终的炉膛出口温度T,其值在950~1050℃之间。其一天的测算值如图2所示。
步骤五:在优化测量方法执行过程中,基于化学班每天四班的煤种分析值实时更新煤质状况,并在每次大修后重新设定锅炉物理参数。
Claims (4)
1.一种基于烟气能量平衡的炉膛出口烟温优化测量方法,其特征在于,包括以下步骤:
步骤一:设置包括炉膛体积、换热面积以及各级换热器的换热面积、金属质量的锅炉炉膛及烟道特性,通过OPC通讯协议,建立OPC客户端向DCS自带的OPC服务器端发送数据请求,进行实时采样获取锅炉运行数据,对数据进行五点平滑处理;
步骤二:在现场控制系统工程师站上建立科学计算平台与人机交互平台,即利用Matlab建立数据处理平台,基于Matlab自带的OPC客户端用于实现数据接收与处理;
步骤三:基于锅炉炉膛的燃烧放热平衡与水冷壁辐射换热平衡,计算出炉膛出口烟温T1;基于锅炉烟道的尾烟气能量、再热器能量平衡以及过热器能量平衡,利用烟气物性参数库,计算出炉膛出口烟温T2;
步骤四:根据正向燃烧放热平衡与反向过热器再热器热交换平衡计算得到炉膛出口烟温T1与T2,采用线性叠加的方式,即T=a×T1+b×T2计算出最终炉膛烟温结果,其中系数a取0.3,系数b取0.7;
步骤五:根据实际燃料特性以及炉膛、烟道、换热器特性随着时间推移所发生的变化修正方法执行中炉膛换热系数等各组参数,即对优化测量方法进行参数滚动更新。
2.根据权利要求1所述的基于烟气能量平衡的炉膛出口烟温优化测量方法,其特征是,所述的采样的时间间隔为5秒。
3.根据权利要求1所述的基于烟气能量平衡的炉膛出口烟温优化测量方法,其特征是,所述的步骤三具体包括:
a)利用正向燃煤在炉膛中放热与传热平衡,计算出炉膛出口温度T1
b)利用反向烟气在各级换热器中的产热平衡,利用尾气能量加上各级换热器传递热量计算出炉膛出口烟气能量,并利用烟气物性反推炉膛出口烟温T2
4.根据权利要求1所述的基于烟气能量平衡的炉膛出口烟温优化测量方法,其特征是,所述的步骤五具体包括:
a)每天利用最新煤质在线分析值对烟气物性参数库进行更新;
b)每次大修以后对炉膛换热系数等各种参数进行重新设置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110107881.3A CN102252779B (zh) | 2011-04-28 | 2011-04-28 | 基于烟气能量平衡的炉膛出口烟温优化测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110107881.3A CN102252779B (zh) | 2011-04-28 | 2011-04-28 | 基于烟气能量平衡的炉膛出口烟温优化测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102252779A true CN102252779A (zh) | 2011-11-23 |
CN102252779B CN102252779B (zh) | 2013-07-03 |
Family
ID=44980179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110107881.3A Expired - Fee Related CN102252779B (zh) | 2011-04-28 | 2011-04-28 | 基于烟气能量平衡的炉膛出口烟温优化测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102252779B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799775A (zh) * | 2012-07-06 | 2012-11-28 | 广东电网公司电力科学研究院 | 利用数值模拟辅助超超临界锅炉屏式过热器管壁温度的计算方法 |
CN103728055A (zh) * | 2013-12-17 | 2014-04-16 | 上海交通大学 | 一种火电机组锅炉炉膛出口烟气能量的实时估计方法 |
CN103728054A (zh) * | 2013-12-14 | 2014-04-16 | 云南新立有色金属有限公司 | 一种高钛渣冶炼大型直流密闭电炉热损失的计算方法 |
CN105276563A (zh) * | 2015-11-11 | 2016-01-27 | 东南大学 | 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法 |
CN106352320A (zh) * | 2016-10-31 | 2017-01-25 | 东南大学 | 一种电站锅炉炉膛结渣分区段实时软测量方法 |
CN116976162A (zh) * | 2023-06-16 | 2023-10-31 | 武汉光谷蓝焰新能源股份有限公司 | 基于生物质锅炉有限元分析的脱硝设计仿真方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060064669A1 (en) * | 2004-09-20 | 2006-03-23 | The Mathworks, Inc | Automatic generation of code for component interfaces in models |
CN1858489A (zh) * | 2006-05-18 | 2006-11-08 | 上海交通大学 | 锅炉汽包水位模型辨识的方法 |
CN101799170A (zh) * | 2010-03-18 | 2010-08-11 | 华北电力大学 | 燃煤锅炉燃料发热量实时校正方法 |
CN101825503A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机表面式加热器的出水及疏水温度测算方法 |
-
2011
- 2011-04-28 CN CN201110107881.3A patent/CN102252779B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060064669A1 (en) * | 2004-09-20 | 2006-03-23 | The Mathworks, Inc | Automatic generation of code for component interfaces in models |
CN1858489A (zh) * | 2006-05-18 | 2006-11-08 | 上海交通大学 | 锅炉汽包水位模型辨识的方法 |
CN101799170A (zh) * | 2010-03-18 | 2010-08-11 | 华北电力大学 | 燃煤锅炉燃料发热量实时校正方法 |
CN101825503A (zh) * | 2010-04-16 | 2010-09-08 | 东南大学 | 汽机表面式加热器的出水及疏水温度测算方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799775A (zh) * | 2012-07-06 | 2012-11-28 | 广东电网公司电力科学研究院 | 利用数值模拟辅助超超临界锅炉屏式过热器管壁温度的计算方法 |
CN102799775B (zh) * | 2012-07-06 | 2016-06-08 | 广东电网公司电力科学研究院 | 利用数值模拟辅助获得超超临界锅炉屏式过热器管壁温度的方法 |
CN103728054A (zh) * | 2013-12-14 | 2014-04-16 | 云南新立有色金属有限公司 | 一种高钛渣冶炼大型直流密闭电炉热损失的计算方法 |
CN103728054B (zh) * | 2013-12-14 | 2016-01-13 | 云南新立有色金属有限公司 | 一种高钛渣冶炼大型直流密闭电炉热损失的计算方法 |
CN103728055A (zh) * | 2013-12-17 | 2014-04-16 | 上海交通大学 | 一种火电机组锅炉炉膛出口烟气能量的实时估计方法 |
CN103728055B (zh) * | 2013-12-17 | 2016-03-02 | 上海交通大学 | 一种火电机组锅炉炉膛出口烟气能量的实时估计方法 |
CN105276563A (zh) * | 2015-11-11 | 2016-01-27 | 东南大学 | 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法 |
CN106352320A (zh) * | 2016-10-31 | 2017-01-25 | 东南大学 | 一种电站锅炉炉膛结渣分区段实时软测量方法 |
CN106352320B (zh) * | 2016-10-31 | 2019-01-25 | 东南大学 | 一种电站锅炉炉膛结渣分区段实时软测量方法 |
CN116976162A (zh) * | 2023-06-16 | 2023-10-31 | 武汉光谷蓝焰新能源股份有限公司 | 基于生物质锅炉有限元分析的脱硝设计仿真方法 |
CN116976162B (zh) * | 2023-06-16 | 2024-01-30 | 武汉光谷蓝焰新能源股份有限公司 | 基于生物质锅炉有限元分析的脱硝设计仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102252779B (zh) | 2013-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102252779B (zh) | 基于烟气能量平衡的炉膛出口烟温优化测量方法 | |
Si et al. | Exergy analysis of a 1000 MW double reheat ultra-supercritical power plant | |
CN103726887B (zh) | 一种燃煤机组汽轮机性能在线监测方法 | |
CN103674333B (zh) | 一种燃煤电站入炉煤低位发热量的实时辨识方法 | |
Chaibakhsh et al. | A simulated model for a once-through boiler by parameter adjustment based on genetic algorithms | |
CN102252784B (zh) | 基于燃煤发热量辨识的火力电站锅炉效率在线监测系统 | |
CN101598688B (zh) | 基于煤质在线测量的锅炉灰污监测及吹灰优化方法 | |
Naserabad et al. | Multi-objective optimization of HRSG configurations on the steam power plant repowering specifications | |
CN102494714B (zh) | 电站锅炉效率与燃煤热值、灰分和水分的同步测算方法 | |
Madejski et al. | Calculation methods of steam boiler operation factors under varying operating conditions with the use of computational thermodynamic modeling | |
Sreepradha et al. | Mathematical model for integrated coal fired thermal boiler using physical laws | |
CN103778334A (zh) | 一种燃煤电站锅炉热效率的实时测量方法 | |
CN105276563A (zh) | 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法 | |
CN106326534A (zh) | 变工况亚临界火电机组的锅炉‑汽机控制模型构建方法 | |
Zhao et al. | Coal-fired utility boiler modelling for advanced economical low-NOx combustion controller design | |
Wang et al. | Flexibility and efficiency co-enhancement of thermal power plant by control strategy improvement considering time varying and detailed boiler heat storage characteristics | |
CN103679549A (zh) | 火电机组节能潜力分析方法 | |
Li et al. | Modeling for the performance evaluation of 600 MW supercritical unit operating No. 0 high pressure heater | |
Eboh et al. | Evaluating improvements in a waste-to-energy combined heat and power plant | |
Wang et al. | Mechanism modeling of optimal excess air coefficient for operating in coal fired boiler | |
CN104122291B (zh) | 超超临界火电机组水冷壁向工质传热速率的实时辨识方法 | |
Zhang et al. | A dynamic heat transfer model to estimate the flue gas temperature in the horizontal flue of the coal-fired utility boiler | |
CN110298502A (zh) | 基于能效最优的锅炉最佳氧量计算方法 | |
CN102032956A (zh) | 锅炉水冷壁吸热量的实时测量方法 | |
CN103728055B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130703 Termination date: 20160428 |