CN112989719A - 一种燃气锅炉机组动态特性建模方法 - Google Patents
一种燃气锅炉机组动态特性建模方法 Download PDFInfo
- Publication number
- CN112989719A CN112989719A CN202110233993.7A CN202110233993A CN112989719A CN 112989719 A CN112989719 A CN 112989719A CN 202110233993 A CN202110233993 A CN 202110233993A CN 112989719 A CN112989719 A CN 112989719A
- Authority
- CN
- China
- Prior art keywords
- model
- steam
- water
- heat
- drum
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012546 transfer Methods 0.000 claims abstract description 23
- 230000009471 action Effects 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 62
- 239000007789 gas Substances 0.000 claims description 30
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 21
- 239000003546 flue gas Substances 0.000 claims description 21
- 238000004088 simulation Methods 0.000 claims description 17
- 230000001174 ascending effect Effects 0.000 claims description 15
- 229920006395 saturated elastomer Polymers 0.000 claims description 8
- 238000013178 mathematical model Methods 0.000 claims description 6
- CDBYLPFSWZWCQE-UHFFFAOYSA-L Sodium Carbonate Chemical compound [Na+].[Na+].[O-]C([O-])=O CDBYLPFSWZWCQE-UHFFFAOYSA-L 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 4
- 239000000779 smoke Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 239000002918 waste heat Substances 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims description 2
- 238000012795 verification Methods 0.000 abstract description 6
- 238000010438 heat treatment Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000000446 fuel Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000002737 fuel gas Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009835 boiling Methods 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000008400 supply water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
本发明涉及综合能源系统设备建模领域,本发明公开一种燃气锅炉机组动态特性建模方法,包括以下步骤:S1、分析燃气锅炉系统的组成部分,建立相应的模型框架;S2、考虑模型的准确性,建立模型部件的传热、导热方程;S3、根据模型部件的动态特性建立动态微分方程模型;S4、根据模型部件的具体作用建立相应的精细化模型;S5、将建立模型在apros系统中进行仿真并通过实际数据验证。本发明给出了燃气锅炉的精细化模型框架,并给出了相应的精细化模型分块,继而在此模型框架的基础上,建立了考虑省煤器、过热器、蒸发器、汽包的燃气锅炉精细化模型方程,并进行了相应的模型验证,即利用apros仿真软件进行模型搭建并与实际案例进行验证比较。
Description
技术领域
本发明涉及综合能源系统设备建模领域,具体的是一种燃气锅炉机组动态特性建模方法。
背景技术
随着经济发展,用户的用电量飞速增长,对电网、发电侧造成巨大的压力。而长时间以来,我国的用户侧新能源又与电力系统规划工作相分离,用户侧的分布式能源对配电网的安全运行产生较大影响。随着综合能源概念的提出,构建能源互联网、解决化石燃料逐渐枯竭所造成的环境问题、使多能源开放互联、实现能源梯级利用成为研究热点。
用户侧分布式能源系统是一种新型的能源系统,其集合新型技术和设备为一体,可实现对输出能量梯级利用,达到更高的能源利用率。目前的分布式能源系统种类繁多,不仅包括以燃气轮机或内燃机为核心的冷热电联产系统,还包括太阳能、风能、生物能等可再生能源综合利用系统,以及由具有极高效率的新型燃料电池与超级电容组成的能源综合利用系统。
在能源互联网框架下的用户侧分布式能源的调度优化,促进能源发展科学规划,促进各类能源系统协调运行,助力用户侧能源实现综合能效最优至关重要。因此,本申请提供一种基于用户侧分布式能源系统的燃气锅炉精细化建模方法。
发明内容
为解决上述背景技术中提到的不足,本发明的目的在于提供一种燃气锅炉机组动态特性建模方法,本发明给出了燃气锅炉的精细化模型框架,并给出了相应的精细化模型分块,继而在此模型框架的基础上,建立了考虑省煤器、过热器、蒸发器、汽包的燃气锅炉精细化模型方程,并进行了相应的模型验证,即利用apros仿真软件进行模型搭建并与实际案例进行验证比较。
本发明的目的可以通过以下技术方案实现:
一种燃气锅炉机组动态特性建模方法,包括以下步骤:
S1、分析燃气锅炉系统的组成部分,建立相应的模型框架;
S2、考虑模型的准确性,建立模型部件的传热、导热方程;
S3、根据模型部件的动态特性建立动态微分方程模型;
S4、根据模型部件的具体作用建立相应的精细化模型;
S5、将建立模型在apros系统中进行仿真并通过实际数据验证。
进一步优选地,所述步骤S1中燃气锅炉系统包含两台燃气热水锅炉,所述燃气热水锅炉包括省煤器、过热器、换热器和汽包。
进一步优选地,所述步骤S2中模型部件为省煤器和过热器,省煤器用于回收所排烟的余热,将锅炉给水加热成汽包压力下的饱和水的受热面,过热器将蒸汽从饱和温度进一步加热至过热温度,所述传热、导热方程具体为:
烟气侧导热微分方程:
式(1)中:T为温度,K;γ为比热容比;y为空间体积,m;Vg为烟气体积,m3;Qgm为烟气侧传热量,W;ug为烟气比热力学能,kJ/kg;λg为烟气导热系数,W/(m·K);cv为烟气侧比定容热容,kJ/(kg·K);ρg为烟气密度,kg/m3;
汽水侧导热微分方程:
式(2)中,Qfm为汽水侧传热量,W;cp为汽水侧比定压热容,kJ/(kg·K);f为汽水;Vf为汽水体积,m3;uf为烟气比热力学能,kJ/kg;λf为汽水侧导热系数,W/(m·K);cp为汽水侧比定容热容,kJ/(kg·K);ρf为汽水密度,kg/m3;
烟气侧换热方程:
Qgm=Aohgm(Tg-Tm) (3),
汽水侧换热方程:
Qfm=Aihfm(Tf-Tm) (4),
式(3)、(4)中,hgm为烟气侧对流换热系数,W/(m2·K);hfm为汽水侧对流换热系数,W/(m2·K);Ai为管内总表面积,m2;A0为管外总表面积,m2。
进一步优选地,所述步骤S3中模型部件为蒸发器,蒸发器中包含下降管,上升管和集总管,所述动态微分方程模型包括:
1)下降管数学模型
动量方程:
式(5)中,Ldc为下降管管长,m;Ddc为下降管管径,m;Adc是下降管面积m2;fdc为下降管阻力系数;mdc为下降管流量,kg/s;vdc为下降管比容,m3/kg;ξ为入口损失系数;pwh为上升管压强,MPa;pdr为汽包压强,MPa;
2)上升管数学模型
动量方程:
式(6)中:Lr为上升管管长,m;Dr为上升管管径,m;Ar是上升管面积m2;fr为上升管阻力系数;mr为上升管流量,kg/s;vr为下降管比容,m3/kg;ξ为入口损失系数;pwh为上升管压强,MPa;pdr为汽包压强,MPa;
能量方程:
式(7)中:ur为上升管比热力学能,kJ/kg;Qr为上升管换热量,W;Mr为上升管流体总质量,kg;hr为上升管对流换热系数,W/(m2·K);hdc为下降管对流换热系数,W/(m2·K);
质量守恒方程:
进一步优选地,所述步骤S4中模型部件为汽包,所述汽包的精细化模型具体建模方法如下:
质量平衡:
式(9)中:Mdr为汽包流体总质量,kg;mec为省煤器流量,kg/s;ms为蒸汽流量,下同;
能量平衡:
式(10)中:udr为汽包比热力学能,kJ/kg;udc为下降管热力学能,kJ/kg;Qmf为壳侧对汽包放热量(W);hec为省煤器对流换热系数,W/(m2·K);hs为蒸汽对流换热系数,W/(m2·K);
汽包水体积:
VW=[(1-xdr)/ρm]Mdr (12),
式(11)、(12)中,Vw为汽包水实际体积,m3;Vdr为汽包水总体积,m3;Ldr为汽包长度;ρw为汽包水密度,kg/m3;rdr是汽包半径,m;Y为汽包液位,m;xdr为汽包干度;α为汽包剖面圆心到汽包水平面与圆剖面相交的点的半径和通过圆心的垂直半径形成的夹角;
汽包水体积关于时间微分变换后得到:
进一步优选地,所述步骤S5具体为:将建立模型在apros系统中进行仿真得到仿真模型数据,再根据历史数据最小二乘-分段线性连续拟合得到实际数据,然后将实际数据与仿真模型数据进行比较,考察仿真模拟精度。
进一步优选地,所述最小二乘-分段线性连续拟合包括以下步骤:
a、首先给定历史样本数据集(xi,yi),i=1、2、3、…、N;并统计xi的最小xmin和最大值xmax;
b、根据线性分段数,将xi等分为L个区间,得到各子区间为:[X0,X1],[X1,X2],…,[XL-1,XL];其中X0=xmin,XL=xmax,Xl=xmin+(xmax-xmin)/L×l;
c、将历史样本集按xi的大小划分到各子区间:若xi≥Xl且xi≤Xl+1,则样本(xi,yi)属于第l分段;
d、将各区间内的xi进行单位化,得到1-(xi-Xl)/(Xl-Xl-1),(xi-Xl)/(Xl-Xl-1),分别记为x′l,i和x″l,i;
e、利用最小二乘技术,求得Y0,Y1,…,YL,使各分段总体拟合误差最小,即:
f、根据端点(X0,Y0),(X1,Y1),,…,(XL,YL)计算各分段斜率ρi。
本发明的有益效果:
本发明给出了燃气锅炉的精细化模型框架,并给出了相应的精细化模型分块,继而在此模型框架的基础上,建立了考虑省煤器、过热器、蒸发器、汽包的燃气锅炉精细化模型方程,并进行了相应的模型验证,即利用apros仿真软件进行模型搭建并与实际案例进行验证比较。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明燃气锅炉机组动态特性建模方法的流程示意图;
图2是本发明燃气锅炉模型总览;
图3是本发明燃气锅炉控制系统;
图4是本发明燃气锅炉动态特性监测结果;
图5是本发明燃气锅炉模型验证结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“开孔”、“上”、“下”、“厚度”、“顶”、“中”、“长度”、“内”、“四周”等指示方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的组件或元件必须具有特定的方位,以特定的方位构造和操作,因此不能理解为对本发明的限制。
一种燃气锅炉机组动态特性建模方法,具体包括以下步骤:
S1、分析燃气锅炉系统的组成部分,建立相应的模型框架:
燃气锅炉主要原理是水通过进水口进入锅炉,经过锅炉加热后达到标准水后通过循环泵送入室内散热器,通过辐射和对流换热进行供暖。通过散热器的水重新回到锅炉里进行加热,然后重新进入散热器,如此循环往复地进行。用户还可以根据供热范围的大小,选择合适的循环水泵。燃气锅炉燃料在炉内燃烧放出来的热量,加热锅内的水,使其汽化成蒸汽。水在锅炉中不断被炉里气体燃料燃烧释放出来的能量加热,温度升高并产生带压蒸汽,由于水的沸点会随压力的升高而升高,因此水蒸气在里面的膨胀受到限制而产生压力形成热动力作为一种能源广泛使用。
热水锅炉燃烧天然气产生高温燃气,通过换热器与供热回水进行热量交换,加热冷却回水,进而对用户供热供暖。系统中包含有两台燃气热水锅炉。模型中燃气热水锅炉包含有省煤器、过热器、换热器、汽包等模型装置。
S2、考虑模型的准确性,建立省煤器、过热器的传热导热方程:
省煤器安装于锅炉尾部烟道下部用于回收所排烟的余热,将锅炉给水加热成汽包压力下的饱和水的受热面,由于省煤器吸收高温烟气的热量,降低烟气的排烟温度,给循环增加一个回热过程。提高吸热平均温度。从而增加循环效率。过热器将蒸汽从饱和温度进一步加热至过热温度。各部件的传热、导热方程具体建模如下所示:
烟气侧导热微分方程
汽水侧导热微分方程
金属侧导热微分方程
烟气侧换热方程
Qgm=Aohgm(Tg-Tm)
汽水侧换热方程
Qfm=Aihfm(Tf-Tm)
S3、,根据模型的动态特性建立蒸发器的动态微分方程模型:
蒸发器中包含下降管,上升管,集总管组成,下降管中的水一般看作饱和水,可以简化为汽包压力下的饱和水。
1)下降管数学模型
动量方程
2)上升管数学模型
动量方程
能量方程
质量守恒方程
S4、根据汽包的具体作用建立相应的精细化模型:
汽包是自然循环锅炉的关键部件.它的工作好坏直接关系到水循环的安全和输出蒸汽的品质,在锅炉中用以进行汽水分离和蒸汽净化,组成水循环回路并蓄存锅水的筒形压力容器。主要作用为接纳省煤器来水,进行汽水分离和向循环回路供水,向过热器输送饱和蒸汽。汽包中存有一定水量,具有一定的热量及工质的储蓄,在工况变动时可减缓汽压变化速度,当给水与负荷短时间不协调时起一定的缓冲作用。汽包的具体建模方法如下所示
质量平衡
能量平衡
汽包水体积
VW=[(1-xdr)/ρm]Mdr
汽包水体积关于时间微分变换后得到:
S5、将建立模型在apros系统中进行仿真并通过实际数据验证:
作为一个可以对自动化以及能源系统进行动态建模的仿真软件,apros可以有效的描述动力系统的动态过程,Apros仿真支撑系统由仿真引擎和图形用户界面组成。仿真引擎包括各种求解器和模型库,而图形用户界面提供了方便的图形组态和运行方式。Apros支持在建厂过程的所有阶段进行动态仿真,从而避免不必要的数据传递以及仿真模型的重新组态。
在实际的锅炉建模精度验证的过程中,一般直接给定的是设备的输入和输出历史数据。然后根据历史数据,拟合得到其函数关系式,再利用分段线性化技术获得各分段的线性关系。
利用分段连续-最小二乘技术得到的函数(粉色虚线)相比传统方法得到的函数(绿色实线)更为符合实际数据分布。
最小二乘-分段线性连续拟合包括以下步骤:
a、首先给定历史样本数据集(xi,yi),i=1、2、3、…、N;并统计xi的最小xmin和最大值xmax;
b、根据线性分段数,将xi等分为L个区间,得到各子区间为:[X0,X1],[X1,X2],…,[XL-1,XL];其中X0=xmin,XL=xmax,Xl=xmin+(xmax-xmin)/L×l;
c、将历史样本集按xi的大小划分到各子区间:若xi≥Xl且xi≤Xl+1,则样本(xi,yi)属于第l分段;
d、将各区间内的xi进行单位化,得到1-(xi-Xl)/(Xl-Xl-1),(xi-Xl)/(Xl-Xl-1),分别记为x′l,i和x″l,i;
e、利用最小二乘技术,求得Y0,Y1,…,YL,使各分段总体拟合误差最小,即:
f、根据端点(X0,Y0),(X1,Y1),,…,(XL,YL)计算各分段斜率ρi。
而后,将实际数据与仿真模型数据进行比较,考察其仿真模拟精度。
本发明实例仿真均采用江苏省南京市某工业园区内某智能楼宇楼宇典型锅炉工况数据为优化对象,该数据是由项目组中的同学处理后得到的。采用apros软件进行详细建模仿真,具体结果如附图5所示。
在本项目中,选取某天实际运行数据作为边界条件,对所搭建仿真模型进行验证,从图中可以看出,仿真模型运行数据与实际运行数据相对误差在6%以内,符合工程实际要求。
在本说明书的描述中,参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。
Claims (7)
1.一种燃气锅炉机组动态特性建模方法,其特征在于,包括以下步骤:
S1、分析燃气锅炉系统的组成部分,建立相应的模型框架;
S2、考虑模型的准确性,建立模型部件的传热、导热方程;
S3、根据模型部件的动态特性建立动态微分方程模型;
S4、根据模型部件的具体作用建立相应的精细化模型;
S5、将建立模型在apros系统中进行仿真并通过实际数据验证。
2.根据权利要求1所述的燃气锅炉机组动态特性建模方法,其特征在于,所述步骤S1中燃气锅炉系统包含两台燃气热水锅炉,所述燃气热水锅炉包括省煤器、过热器、换热器和汽包。
3.根据权利要求1所述的燃气锅炉机组动态特性建模方法,其特征在于,所述步骤S2中模型部件为省煤器和过热器,省煤器用于回收所排烟的余热,将锅炉给水加热成汽包压力下的饱和水的受热面,过热器将蒸汽从饱和温度进一步加热至过热温度,所述传热、导热方程具体为:
烟气侧导热微分方程:
式(1)中:T为温度,K;γ为比热容比;y为空间体积,m;Vg为烟气体积,m3;Qgm为烟气侧传热量,W;ug为烟气比热力学能,kJ/kg;λg为烟气导热系数,W/(m·K);cv为烟气侧比定容热容,kJ/(kg·K);ρg为烟气密度,kg/m3;
汽水侧导热微分方程:
式(2)中,Qfm为汽水侧传热量,W;cp为汽水侧比定压热容,kJ/(kg·K);f为汽水;Vf为汽水体积,m3;uf为烟气比热力学能,kJ/kg;λf为汽水侧导热系数,W/(m·K);cp为汽水侧比定容热容,kJ/(kg·K);ρf为汽水密度,kg/m3;
烟气侧换热方程:
Qgm=Aohgm(Tg-Tm) (3),
汽水侧换热方程:
Qfm=Aihfm(Tf-Tm) (4),
式(3)、(4)中,hgm为烟气侧对流换热系数,W/(m2·K);hfm为汽水侧对流换热系数,W/(m2·K);Ai为管内总表面积,m2;A0为管外总表面积,m2。
4.根据权利要求1所述的燃气锅炉机组动态特性建模方法,其特征在于,所述步骤S3中模型部件为蒸发器,蒸发器中包含下降管,上升管和集总管,所述动态微分方程模型包括:
1)下降管数学模型
动量方程:
式(5)中,Ldc为下降管管长,m;Ddc为下降管管径,m;Adc是下降管面积m2;fdc为下降管阻力系数;mdc为下降管流量,kg/s;vdc为下降管比容,m3/kg;ξ为入口损失系数;pwh为上升管压强,MPa;pdr为汽包压强,MPa;
2)上升管数学模型
动量方程:
式(6)中:Lr为上升管管长,m;Dr为上升管管径,m;Ar是上升管面积m2;fr为上升管阻力系数;mr为上升管流量,kg/s;vr为下降管比容,m3/kg;ξ为入口损失系数;pwh为上升管压强,MPa;pdr为汽包压强,MPa;
能量方程:
式(7)中:ur为上升管比热力学能,kJ/kg;Qr为上升管换热量,W;Mr为上升管流体总质量,kg;hr为上升管对流换热系数,W/(m2·K);hdc为下降管对流换热系数,W/(m2·K);
质量守恒方程:
5.根据权利要求1所述的燃气锅炉机组动态特性建模方法,其特征在于,所述步骤S4中模型部件为汽包,所述汽包的精细化模型具体建模方法如下:
质量平衡:
式(9)中:Mdr为汽包流体总质量,kg;mec为省煤器流量,kg/s;ms为蒸汽流量,下同;
能量平衡:
式(10)中:udr为汽包比热力学能,kJ/kg;udc为下降管热力学能,kJ/kg;Qmf为壳侧对汽包放热量(W);hec为省煤器对流换热系数,W/(m2·K);hs为蒸汽对流换热系数,W/(m2·K);
汽包水体积:
VW=[(1-xdr)/ρm]Mdr (12),
式(11)、(12)中,Vw为汽包水实际体积,m3;Vdr为汽包水总体积,m3;Ldr为汽包长度;ρw为汽包水密度,kg/m3;rdr是汽包半径,m;Y为汽包液位,m;xdr为汽包干度;α为汽包剖面圆心到汽包水平面与圆剖面相交的点的半径和通过圆心的垂直半径形成的夹角;
汽包水体积关于时间微分变换后得到:
6.根据权利要求1所述的燃气锅炉机组动态特性建模方法,其特征在于,所述步骤S5具体为:将建立模型在apros系统中进行仿真得到仿真模型数据,再根据历史数据最小二乘-分段线性连续拟合得到实际数据,然后将实际数据与仿真模型数据进行比较,考察仿真模拟精度。
7.根据权利要求6所述的燃气锅炉机组动态特性建模方法,其特征在于,所述最小二乘-分段线性连续拟合包括以下步骤:
a、首先给定历史样本数据集(xi,yi),i=1、2、3、…、N;并统计xi的最小xmin和最大值xmax;
b、根据线性分段数,将xi等分为L个区间,得到各子区间为:[X0,X1],[X1,X2],…,[XL-1,XL];其中X0=xmin,XL=xmax,Xl=xmin+(xmax-xmin)/L×l;
c、将历史样本集按xi的大小划分到各子区间:若xi≥Xl且xi≤Xl+1,则样本(xi,yi)属于第l分段;
d、将各区间内的xi进行单位化,得到1-(xi-Xl)/(Xl-Xl-1),(xi-Xl)/(Xl-Xl-1),分别记为x′l,i和x″l,i;
e、利用最小二乘技术,求得Y0,Y1,…,YL,使各分段总体拟合误差最小,即:
f、根据端点(X0,Y0),(X1,Y1),,…,(XL,YL)计算各分段斜率ρi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110233993.7A CN112989719B (zh) | 2021-03-03 | 2021-03-03 | 一种燃气锅炉机组动态特性建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110233993.7A CN112989719B (zh) | 2021-03-03 | 2021-03-03 | 一种燃气锅炉机组动态特性建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112989719A true CN112989719A (zh) | 2021-06-18 |
CN112989719B CN112989719B (zh) | 2023-12-22 |
Family
ID=76352190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110233993.7A Active CN112989719B (zh) | 2021-03-03 | 2021-03-03 | 一种燃气锅炉机组动态特性建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112989719B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114294637A (zh) * | 2022-01-04 | 2022-04-08 | 华润电力技术研究院有限公司 | 一种基于机器学习的低温省煤器状态监测系统及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106326534A (zh) * | 2016-08-15 | 2017-01-11 | 上海交通大学 | 变工况亚临界火电机组的锅炉‑汽机控制模型构建方法 |
CN107292006A (zh) * | 2017-06-08 | 2017-10-24 | 西安交通大学 | 一种超临界锅炉的启动计算方法及系统 |
CN107657104A (zh) * | 2017-09-20 | 2018-02-02 | 浙江浙能台州第二发电有限责任公司 | 基于在线支持向量机的锅炉燃烧系统动态建模方法 |
CN109977560A (zh) * | 2019-03-29 | 2019-07-05 | 江苏方天电力技术有限公司 | 一种余热锅炉动态运行特性建模方法 |
CN111523248A (zh) * | 2020-05-12 | 2020-08-11 | 国电新能源技术研究院有限公司 | 燃煤电站动态机理模型建模方法 |
-
2021
- 2021-03-03 CN CN202110233993.7A patent/CN112989719B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106326534A (zh) * | 2016-08-15 | 2017-01-11 | 上海交通大学 | 变工况亚临界火电机组的锅炉‑汽机控制模型构建方法 |
CN107292006A (zh) * | 2017-06-08 | 2017-10-24 | 西安交通大学 | 一种超临界锅炉的启动计算方法及系统 |
CN107657104A (zh) * | 2017-09-20 | 2018-02-02 | 浙江浙能台州第二发电有限责任公司 | 基于在线支持向量机的锅炉燃烧系统动态建模方法 |
CN109977560A (zh) * | 2019-03-29 | 2019-07-05 | 江苏方天电力技术有限公司 | 一种余热锅炉动态运行特性建模方法 |
CN111523248A (zh) * | 2020-05-12 | 2020-08-11 | 国电新能源技术研究院有限公司 | 燃煤电站动态机理模型建模方法 |
Non-Patent Citations (2)
Title |
---|
崔凝;王兵树;高建强;贡献;: "大容量余热锅炉动态模型的研究与应用", 中国电机工程学报, no. 19 * |
王亚欧;陶谦;陈波;管诗骈;张恩先;: "自然循环双压燃气轮机余热锅炉动态特性分析", 中国电力, no. 02 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114294637A (zh) * | 2022-01-04 | 2022-04-08 | 华润电力技术研究院有限公司 | 一种基于机器学习的低温省煤器状态监测系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112989719B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Modeling and dynamic simulation of a steam generation system for a parabolic trough solar power plant | |
CN106998079A (zh) | 一种热电联合优化调度模型的建模方法 | |
CN103868389A (zh) | 一种独立熔盐蓄热电站 | |
CN112228853B (zh) | 多孔介质传蓄热装置、传蓄热发电系统及储能电站 | |
CN206073441U (zh) | 强制循环燃油燃气水管热水锅炉 | |
CN111680417B (zh) | 一种适用于钠水直流蒸汽发生器的热工设计方法 | |
Dzierwa et al. | Transient CFD simulation of charging hot water tank | |
CN112989719B (zh) | 一种燃气锅炉机组动态特性建模方法 | |
CN203772087U (zh) | 一种独立熔盐蓄热电站 | |
CN86105222A (zh) | 利用排烟余热的低压节能复合系统 | |
CN103925820B (zh) | 一种具有中间蓄热的多级闪蒸系统 | |
CN107906498A (zh) | 超临界二氧化碳循环流化床燃煤锅炉及其驱动的发电系统 | |
CN209840064U (zh) | 一种交叉混流过热器 | |
Xu et al. | The Badaling 1MW parabolic trough solar thermal power pilot plant | |
Cao et al. | Thermal system for comprehensive utilization of boiler and steam turbine energy of coal-fired power units based on the S-CO2 cycle | |
CN113051752B (zh) | 一种与煤电耦合的高压空气储能系统最优热源的确定方法 | |
Pang et al. | Dynamic characteristics of a solar‐aided coal‐fired power generation system under direct normal irradiance (DNI) feedforward control | |
CN112446097B (zh) | 一种蒸汽发生器体积和负荷的多目标优化方法 | |
Rączka et al. | Methods of thermal calculations for a condensing waste-heat exchanger | |
CN212430800U (zh) | 一种用于大功率光热电站的蒸汽发生装置 | |
CN210921355U (zh) | 一种蒸汽锅炉余热回收装置 | |
CN106288371A (zh) | 强制循环燃油燃气水管热水锅炉 | |
CN206695662U (zh) | 一种基于相变储热单元的模块化梯级储热装置 | |
CN107508370B (zh) | 电锅炉供电和智能监控系统及其监控方法 | |
Vuillerme et al. | Experimental results of 2018 test campaign of the ALSOLEN 450 prototype for direct steam generation and thermal storage |
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 |