CN113669740B - 一种预测燃烧室内局部最高温度的方法 - Google Patents

一种预测燃烧室内局部最高温度的方法 Download PDF

Info

Publication number
CN113669740B
CN113669740B CN202110968277.3A CN202110968277A CN113669740B CN 113669740 B CN113669740 B CN 113669740B CN 202110968277 A CN202110968277 A CN 202110968277A CN 113669740 B CN113669740 B CN 113669740B
Authority
CN
China
Prior art keywords
combustion chamber
fuel
highest temperature
prediction function
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.)
Active
Application number
CN202110968277.3A
Other languages
English (en)
Other versions
CN113669740A (zh
Inventor
李清海
杨潇潇
周会
张衍国
丛堃林
张时语
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN202110968277.3A priority Critical patent/CN113669740B/zh
Publication of CN113669740A publication Critical patent/CN113669740A/zh
Application granted granted Critical
Publication of CN113669740B publication Critical patent/CN113669740B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23GCREMATION FURNACES; CONSUMING WASTE PRODUCTS BY COMBUSTION
    • F23G5/00Incineration of waste; Incinerator constructions; Details, accessories or control therefor
    • F23G5/50Control or safety arrangements

Abstract

本发明公开了一种预测燃烧室内局部最高温度的方法。首先获取基础特性参数,包括燃烧室的定性容积V和截面积F、过量空气系数α、二次风率β、燃料供给量D、燃料热值等,然后通过热电偶测量获得相应的温度分布,并根据温度分布获取局部最高温度Tmax。基于机理模型和数据分析建立燃烧室局部最高温度预测函数形式Tmax=F(x1、x2、x3、α、β),并对预测函数F进行标定,获得具有具体形式的预测函数,并用于预测燃烧室内局部最高温度。本发明具有预测准确率高、便于集成到DCS或者PLC系统,使用方便、友好等优点。

Description

一种预测燃烧室内局部最高温度的方法
技术领域
本发明涉及预测燃烧温度的方法,尤其涉及一种预测燃烧室和料层内局部最高温度的方法,特别适合于多种垃圾混烧、生物质燃烧、煤燃烧等领域。
背景技术
垃圾焚烧已成为中国生活垃圾处理的主流工艺,同时,中国工业有机固废产量逐年增加,这些有机固废种类多、处置难、潜在环境污染严重,亟待解决。如果充分发挥垃圾焚烧设施处置能力和优势,利用垃圾焚烧设备协同处置工业有机固废,则既可有效利用焚烧设施产能,又能处置工业有机固废,不仅具有显著的经济、环境和社会效益,也是实现“无废城市”的重要手段。
垃圾分类后组分、产量及热值变化明显,特别是热值增加较多。任何固废在燃烧室内燃烧,都要经历干燥、热解、固定碳燃烧燃烧、挥发分燃烧、燃尽等几个阶段,这是共同的,但燃烧室内的燃烧过程组织不同,燃烧过程也会略有差异,比如流化床燃烧和炉排燃烧过程就有显著的不同,其干燥、热解、燃烧过程不同也导致温度分布不同。另外,分类垃圾和有机固废自身在挥发分、固定碳、水分、热值方面等存在较大差异,因而燃烧过程必然不同。典型的工业有机固废热值较高,在协同焚烧时会引起燃烧室内燃烧组织的差别,导致炉膛或炉排局部温度快速飞升或者升高幅度超限,严重时可烧毁炉排。
归纳而言,燃烧室内局部最高温度变化可定义为“飞温”,这种飞温包含两个基本特征:1)升温速率较快,可能超过5℃/s,当飞温发生时操作人员很难及时采取动作抑制升温速率;
2)升温幅度过高,可能超过50%,升温幅度越高产生的危害性越大。“飞温”可能与以下因素均有关联:1)几何结构因素,包括燃烧室结构、烟气流动导向方式、烟气冷却形式、炉拱布置、余热锅炉结构等的影响;2)混合燃料的总体热值、燃料中个别组分的热值、混合程度;
3)配风情况,即过量空气系数、二次风率等因素;4)负荷情况,比如容积热负荷、截面热负荷、炉排热负荷等;5)其他一些偶然因素,比如局部漏风、炉排局部卡涩、偏烧等。
为了主动采取措施提前避免燃烧室内(包括炉膛或者炉排内)局部高温的产生,需要一种在线预测燃烧室局部最高温度的方法。
发明内容
本发明提供一种在线预测燃烧室局部最高温度的方法,旨在根据燃烧室几何因素、燃料热值、配风情况、负荷情况以及其他因素预测燃烧室内的局部最高温度。
本发明是通过以下技术方案实现的:
一种预测燃烧室内局部最高温度的方法,燃烧室内设置有若干热电偶,包括:
获取燃烧室几何特征参数,包括燃烧室的定性容积V和定性截面积F;
获取燃料热值特征参数,包括燃料热值Qmix和Qmax;
获取配风参数,包括燃烧室过量空气系数α、二次风率β;
获取负荷参数,包括燃料供给量D;
获取燃烧室局部最高温度测量参数:通过所述热电偶获得相应的温度分布,并根据温度分布获取局部最高温度测量值Tmax;
对获取的参数进行无量纲化,得到最小负荷指数x1=D*Qmix/V,最大负荷指数x2=D*Qmax/F,燃料不均匀指数x3=Qmax/Qmix;
建立燃烧室局部最高温度预测函数形式Tmax=F(x1、x2、x3、α、β),并对预测函数F进行标定,获得具有具体形式的预测函数,并用于预测燃烧室内局部最高温度。
进一步的,所述预测函数F设定为指数函数的组合Tmax=(x1)a(x2)b(x3)cαdβe,其中a、b、c、d、e为待定常数,通过标定获得。
上述技术方案中,标定过程通过获取不同t时刻对应的参数,采用多元线性回归方法进行;所述参数包括配风参数、负荷参数。
上述技术方案中,所述燃料为单一燃料时,Qmin=Qmax;所述燃料为混合燃料时,随机选取不多于5种燃料组分测试其热值,得到Qmin和Qmax。
本发明的有益效果包括:1)基于机理模型和数据分析,综合多种因素预测燃烧室内局部最高温度,预测准确率高;2)通过运行参数提前预判燃烧室内的局部最高温度,可以降低炉膛内密相区结焦或者炉排内高温烧毁炉排片的风险;3)模型便于集成到DCS或者PLC系统,使用方便、友好。
具体实施方式
下面对本发明的具体实施方式及工作过程作进一步的说明。
燃烧室内的温度场、局部最高温度一般与以下因素均有关联:1)几何因素,包括燃烧室结构、容积、面积等;2)混合燃料的总体热值、燃料中个别组分的热值、混合程度;3)配风情况,即过量空气系数、二次风率等因素;4)负荷情况,比如容积热负荷、截面热负荷、炉排热负荷等。这些因素是决定因变量燃烧时内局部最高温度的自变量,可以通过测量、计算等显式获得,而其他如烟气流动导向方式、烟气冷却形式、燃烧方式(流化床、炉排等)、炉拱布置(顺流、逆流、混合流)、余热锅炉结构等的影响可隐式隐含在温度预测函数中。通过测量若干组自变量、因变量的关系,则可标定出温度预测函数。当这些自变量不发生变化时,正常情况不会发生“飞温”。其他一些偶然因素,比如局部漏风、局部卡涩等导致的飞温可归结为偶然现象,不作为本发明预测范围。
为了获得燃烧室局部飞温预测,需要在该区域(即待测区域内)设置能够获得标定所需温度的热电偶。根据后续温度曲线标定需要,通常设置至少3个测温点,测温点内插入热电偶,热电偶可以采用B型、S型、K型等。
获取燃烧室几何特征参数,包括燃烧室的定性容积V和截面积F。几何特征参数为静态参数,不随燃料、工况变化而变化。
获取燃料热值特征参数,包括燃料热值Qmix和Qmax。当燃料为单一燃料时,Qmin=Qmax。当燃料为混合燃料时,少于5种燃料组分混合的,测试每种燃料的热值,以得到Qmin和Qmax;而多于5种时,随机选取不多于5种燃料组分测试其热值,得到Qmin和Qmax。需要解释的是,如城市生活垃圾之类的燃料视为混合燃料,其中含有多种燃料组分,诸如纸张、纺织品、厨余等,按行业相应的规范取样检测。当换燃料时需要重新检测。
获取配风参数,包括燃烧室过量空气系数α、二次风率β。获取负荷参数,包括燃料供给量D。配风参数和负荷参数可以通过燃烧室的既有仪表测量系统DCS或者PLC读取和计算。
获取燃烧室局部最高温度参数:通过热电偶获得相应的温度分布,并根据温度分布获取测量的局部最高温度测量值Tmax。该温度是根据热电偶所测温度进行拟合所得,故视为测量值。
对获取的参数进行无量纲化,得到最小负荷指数x1=D*Qmix/V,最大负荷指数x2=D*Qmax/F,燃料不均匀指数x3=Qmax/Qmix。
建立燃烧室局部最高温度预测函数形式Tmax=F(x1、x2、x3、α、β),在不同t时刻下获得若干组最高温度测量值Tmax及其对应的x1、x2、x3、α、β,即可采用多元线性回归方法拟合,获得具有具体形式的预测函数F(),这个过程即视为对预测函数F进行标定。
进一步的,可以将预测函数F设定为指数函数的组合Tmax=(x1)a(x2)b(x3)cαdβe,通过标定获得待定常数a、b、c、d、e。
获得了具有具体形式的预测函数Tmax=F()即可用于预测燃烧室内任意工况下的局部最高温度。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种预测燃烧室内局部最高温度的方法,燃烧室内设置有若干热电偶,其特征在于,包括:
获取燃烧室几何特征参数,包括燃烧室的定性容积V和定性截面积F;
获取燃料热值特征参数,包括燃料热值Qmix和Qmax;
获取配风参数,包括燃烧室过量空气系数α、二次风率β;
获取负荷参数,包括燃料供给量D;
获取燃烧室局部最高温度测量参数:通过所述热电偶获得相应的温度分布,并根据温度分布获取局部最高温度测量值Tmax;
对获取的参数进行无量纲化,得到最小负荷指数x1=D*Qmix/V,最大负荷指数x2=D*Qmax/F,燃料不均匀指数x3=Qmax/Qmix;
建立燃烧室局部最高温度预测函数形式Tmax=F(x1、x2、x3、α、β),并对预测函数F进行标定,获得具有具体形式的预测函数,并用于预测燃烧室内局部最高温度。
2.根据权利要求1所述的方法,其特征在于,所述预测函数F设定为指数函数的组合Tmax=(x1)a(x2)b(x3)cαdβe,其中a、b、c、d、e为待定常数,通过标定获得。
3.根据权利要求1或2所述的方法,其特征在于,标定过程通过获取不同t时刻对应的参数,采用多元线性回归方法进行;所述参数包括配风参数、负荷参数。
4.根据权利要求1所述的方法,其特征在于,所述燃料为单一燃料时,Qmin=Qmax;所述燃料为混合燃料时,随机选取不多于5种燃料组分测试其热值,得到Qmin和Qmax。
CN202110968277.3A 2021-08-23 2021-08-23 一种预测燃烧室内局部最高温度的方法 Active CN113669740B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110968277.3A CN113669740B (zh) 2021-08-23 2021-08-23 一种预测燃烧室内局部最高温度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110968277.3A CN113669740B (zh) 2021-08-23 2021-08-23 一种预测燃烧室内局部最高温度的方法

Publications (2)

Publication Number Publication Date
CN113669740A CN113669740A (zh) 2021-11-19
CN113669740B true CN113669740B (zh) 2022-08-12

Family

ID=78545101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110968277.3A Active CN113669740B (zh) 2021-08-23 2021-08-23 一种预测燃烧室内局部最高温度的方法

Country Status (1)

Country Link
CN (1) CN113669740B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI824487B (zh) * 2022-04-15 2023-12-01 中國鋼鐵股份有限公司 溫控決策指引方法及系統

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009139280A (ja) * 2007-12-07 2009-06-25 Takenaka Komuten Co Ltd 温度予測装置、方法及びプログラム
CN103267771A (zh) * 2013-05-09 2013-08-28 西安交通大学 一种基于试验的临界着火温度确定方法
CN106897537A (zh) * 2017-03-14 2017-06-27 清华大学 含三维或曲面外形结构体的温度场与热流同时重构方法
CN113094935A (zh) * 2021-01-06 2021-07-09 武汉科技大学 一种用于“v”形坡隧道内顶板最高温度预测的方法及应用

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102629104B (zh) * 2011-12-01 2013-10-02 燕山大学 一种水泥回转窑煅烧预测控制系统及方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009139280A (ja) * 2007-12-07 2009-06-25 Takenaka Komuten Co Ltd 温度予測装置、方法及びプログラム
CN103267771A (zh) * 2013-05-09 2013-08-28 西安交通大学 一种基于试验的临界着火温度确定方法
CN106897537A (zh) * 2017-03-14 2017-06-27 清华大学 含三维或曲面外形结构体的温度场与热流同时重构方法
CN113094935A (zh) * 2021-01-06 2021-07-09 武汉科技大学 一种用于“v”形坡隧道内顶板最高温度预测的方法及应用

Also Published As

Publication number Publication date
CN113669740A (zh) 2021-11-19

Similar Documents

Publication Publication Date Title
Carroll et al. Air staging to reduce emissions from energy crop combustion in small scale applications
CN113669740B (zh) 一种预测燃烧室内局部最高温度的方法
JP5996762B1 (ja) 廃棄物の燃焼制御方法およびこれを適用した燃焼制御装置
US4022569A (en) Calcination of coke
Junga et al. Uncertainty estimation of the efficiency of small-scale boilers
JP6342367B2 (ja) 廃棄物の発熱量推算方法およびこれを用いた廃棄物処理装置
Lundgren et al. Experimental studies of a biomass boiler suitable for small district heating systems
Kraszkiewicz et al. Combustion of plant biomass pellets on the grate of a low power boiler
Feldmeier et al. Categorization of small-scale biomass combustion appliances by characteristic numbers
EP2385321A2 (en) A method for regulating the combustion process in solid fuel central heating boilers
US6288528B1 (en) Method and system for evaluating a condition of a combustion vessel
Bourgonje et al. The new method to characterize the gas emissions during torrefaction real-time
KR101053292B1 (ko) 가열로에서의 가연 가스 연소 제어 방법
Mack et al. Influence of pellet length, content of fines, and moisture content on emission behavior of wood pellets in a residential pellet stove and pellet boiler
Panwar et al. Experimental investigation of open core downdraft biomass gasifier for food processing industry
CN110598185A (zh) 一种焦炉炼焦单耗影响因素排序方法及系统
Kajda-Szcześniak et al. Characteristics of the combustion process of woodwork waste in the installation of thermal treatment of municipal solid waste (tpok)
Holubcik et al. The impact of bark content in wood pellets on emission production during combustion in small heat source
Petrocelli et al. A note on calculation of efficiency and emissions from wood and wood pellet stoves
Skopec et al. Specific emissions from biomass combustion
Kraszkiewicz The combustion of wood biomass in low power coal-fired boilers
Naď et al. Thermal load non-uniformity estimation for superheater tube bundle damage evaluation
Väätäinen et al. The effects of air staging and combustion air control on black carbon and other particulate and gaseous emissions from a sauna stove
EP3214370B1 (de) Verfahren und vorrichtung zur verbrennung von festen organischen brennstoffen
Horák et al. Real measurement of carbon monoxide, total suspended particulate, and thermal efficiency in modern biomass household boilers

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