一种基于线性回归的余热锅炉蒸发器变工况建模方法技术

本发明专利技术公开了一种基于线性回归的余热锅炉蒸发器变工况建模方法，包括：步骤1，利用传统的余热锅炉蒸发器变工况热力计算模型梳理影响余热锅炉蒸发器变工况热力性能的边界参数；步骤2，确定的边界参数的变化范围内，进行余热锅炉蒸发器变工况热力计算，得到一系列计算结果；步骤3，得到多元线性回归算法的特征和多元线性回归的预测结果；步骤4，将步骤3中确定的多元线性回归算法的特征与预测结果作为多元线性回归算法的数据集，并将数据集分作训练集与测试集，开展多元线性回归算法的训练，并基于训练结果构建余热锅炉蒸发器变工况热力计算的简化计算模型。本发明专利技术用于余热锅炉的变工况热力性能分析、在线性能监测等。在线性能监测等。在线性能监测等。

【技术实现步骤摘要】
一种基于线性回归的余热锅炉蒸发器变工况建模方法


[0001]本专利技术属于燃气
‑
蒸汽联合循环机组领域，具体涉及一种基于线性回归的余热锅炉蒸发器变工况建模方法。

技术介绍

[0002]余热锅炉是燃气
‑
蒸汽联合循环机组中一个重要的换热设备，它回收燃气轮机的排气余热，并借以加热汽轮机系统的给水，使其产生高温高压的水蒸汽，进而送至汽轮机中做功，产生一部分机械功，不仅能增加联合循环机组的功率，而且能提高燃料的化学能与机械能之间的转换效率。
[0003]燃气
‑
蒸汽联合循环机组负荷、环境气象条件、燃料组分等因素的变化，将导致燃气轮机常常在偏离设计工况下运行，因此进入余热锅炉的烟气温度、流量以及烟气组分都将经常变化；此外，环境条件的变化，还将影响汽轮机排汽背压的变化，进而影响汽轮机、余热锅炉的热力性能。此外对于热电联供的联合循环机组，其供热量的变化，也会直接影响余热锅炉的热力性能。以上因素都导致余热锅炉时常处于变工况运行状态。因此研究余热锅炉在变工况下的热力性能计算方法，对于预测、分析余热锅炉甚至联合循环机组的变工况下的热力性能是至关重要的。
[0004]余热锅炉由省煤器、蒸发器和过热器组成。因此余热锅炉在变工况下的热力计算也由省煤器、蒸发器和过热器的变工况热力计算联立组成。在传统的余热锅炉蒸发器变工况热力计算中，烟气侧、水侧换热系数需要经过复杂、繁琐的迭代计算，导致传统的蒸发器变工况热力计算方法并不适用于快速的热力性能分析场景，如余热锅炉的在线性能监测等，因此亟需提出一种可用于快速、简便以及准确的建立余热锅炉蒸发器变工况热力计算模型的方法。

技术实现思路

[0005]本专利技术所要解决的技术问题是：针对传统的余热锅炉蒸发器变工况热力计算模型中，烟气侧和水侧换热系数需要经过复杂的迭代计算，导致余热锅炉蒸发器变工况热力计算过程较为繁琐的问题，提供了一种基于线性回归的余热锅炉蒸发器变工况建模方法。
[0006]本专利技术采用如下技术方案来实现的：
[0007]一种基于线性回归的余热锅炉蒸发器变工况建模方法，包括以下步骤：
[0008]步骤1，利用传统的余热锅炉蒸发器变工况热力计算模型梳理影响余热锅炉蒸发器变工况热力性能的边界参数；
[0009]步骤2，采用传统的余热锅炉蒸发器变工况热力计算模型，并在步骤1确定的边界参数的变化范围内，进行余热锅炉蒸发器变工况热力计算，得到一系列计算结果；
[0010]步骤3，取余热锅炉蒸发器设计工况下的热力计算结果作为基准工况，将步骤2余热锅炉蒸发器变工况热力计算时所取边界参数相对于基准工况边界参数的变化量并进行对数变换后作为多元线性回归算法的特征，步骤2余热锅炉蒸发器变工况热力计算结果中
的总体传热系数相对于基准工况传热系数的变化量并进行对数变换后作为多元线性回归的预测结果；
[0011]步骤4，将步骤3中确定的多元线性回归算法的特征与预测结果作为多元线性回归算法的数据集，并将数据集分作训练集与测试集，开展多元线性回归算法的训练，并基于训练结果构建余热锅炉蒸发器变工况热力计算的简化计算模型。
[0012]本专利技术进一步的改进在于，步骤1和2中，传统的余热锅炉蒸发器变工况热力计算模型为：
[0013]根据余热锅炉蒸发器变工况下的水侧参数以及蒸发器管束直径、管束数量几何参数求解得蒸发器水侧传热系数：
[0014][0015]式中：h
w
为变工况蒸发器中水侧传热系数，W/m2.℃；c
pw
为蒸发器水侧比热，J/kg.℃；μ
w
为蒸发器水侧黏度，kg/m.s；k
w
为蒸发器水侧热导率，W/m.℃；di为蒸发器管束内径，mm。
[0016]本专利技术进一步的改进在于，步骤1和2中，余热锅炉蒸发器变工况下的烟气侧传热系数计算公式为：
[0017][0018]式中：h
g
为变工况蒸发器中烟气侧传热系数，W/m2.℃；C1、C3、C5为系数；d
o
为蒸发器管束外径，mm；h为翅片高度，mm；T
g
为进出口烟气温度平均值，℃；Tf为翅片处温度，℃；W
g
为烟气流量，kg/s；c
pg
为烟气侧比热，J/kg.℃；μ
g
为蒸发器水侧黏度，kg/m.s；k
g
为蒸发器水侧热导率，W/m.℃。
[0019]本专利技术进一步的改进在于，步骤1和2中，余热锅炉蒸发器变工况下的总体传热系数计算公式为：
[0020][0021]式中：U
evp
为变工况蒸发器中总体传热系数，W/m2.℃；A
t
为单位翅片长度下翅片的外表面面积，m2/m；A
i
为单位翅片长度下翅片的内表面面积，m2/m；ff
i
为翅片管内的污垢系数，m2℃/W；A
w
为单位翅片长度下翅片的内、外表面积平均值，m2/m；K
m
为翅片壁面导热系数，W/m2.℃；ff
o
为翅片管外的污垢系数，m2K/W；η为翅片的效率。
[0022]本专利技术进一步的改进在于，步骤1和2中，基于余热锅炉蒸发器变工况下的总体换热系数计算得蒸发器的对数平均温差为：
[0023][0024]式中：LMTD为蒸发器的对数平均温差；e
w2
为蒸发器进口水焓，kJ/kg；e
w1
为蒸发器出口水焓，kJ/kg；A
evp
为蒸发器换热面积，m2。
[0025]本专利技术进一步的改进在于，步骤2中，蒸发器的对数平均温差还可通过余热锅炉蒸
发器的进、出口烟气侧、水侧参数进行热力计算得到：
[0026][0027]式中：T
g1
为蒸发器进口烟气温度，℃；T
g2
为蒸发器出口烟气温度，℃；T
w1
为蒸发器进口水温度，℃。
[0028]本专利技术进一步的改进在于，步骤3中，计算公式如下：
[0029]将余热锅炉设计工况下的计算结果作为基准工况，即：
[0030]U
b
＝f(W
g
，b，cp
g，b
，v
g，b
，km
g，b
，W
w，b
，cp
w，b
，v
w，b
，km
w，b
)
[0031]式中：U
b
为基准工况下蒸发器总体传热系数，W/K；W
g，b
为基准工况下烟气侧流量，kg/s；cp
g，b
为基准工况下烟气侧比热，J/kg/K；v
g，b
为基准工况下烟气侧粘度，Pa.s；km
g，b
为基准工况下烟气侧导热系数，W/m/K；W
w，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，包括以下步骤：步骤1，利用传统的余热锅炉蒸发器变工况热力计算模型梳理影响余热锅炉蒸发器变工况热力性能的边界参数；步骤2，采用传统的余热锅炉蒸发器变工况热力计算模型，并在步骤1确定的边界参数的变化范围内，进行余热锅炉蒸发器变工况热力计算，得到一系列计算结果；步骤3，取余热锅炉蒸发器设计工况下的热力计算结果作为基准工况，将步骤2余热锅炉蒸发器变工况热力计算时所取边界参数相对于基准工况边界参数的变化量并进行对数变换后作为多元线性回归算法的特征，步骤2余热锅炉蒸发器变工况热力计算结果中的总体传热系数相对于基准工况传热系数的变化量并进行对数变换后作为多元线性回归的预测结果；步骤4，将步骤3中确定的多元线性回归算法的特征与预测结果作为多元线性回归算法的数据集，并将数据集分作训练集与测试集，开展多元线性回归算法的训练，并基于训练结果构建余热锅炉蒸发器变工况热力计算的简化计算模型。2.根据权利要求1所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤1和2中，传统的余热锅炉蒸发器变工况热力计算模型为：根据余热锅炉蒸发器变工况下的水侧参数以及蒸发器管束直径、管束数量几何参数求解得蒸发器水侧传热系数：式中：h
w
为变工况蒸发器中水侧传热系数，W/m2.℃；c
pw
为蒸发器水侧比热，J/kg.℃；μ
w
为蒸发器水侧黏度，kg/m.s；k
w
为蒸发器水侧热导率，W/m.℃；d
i
为蒸发器管束内径，mm。3.根据权利要求2所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤1和2中，余热锅炉蒸发器变工况下的烟气侧传热系数计算公式为：式中：h
g
为变工况蒸发器中烟气侧传热系数，W/m2.℃；C1、C3、C5为系数；d
o
为蒸发器管束外径，mm；h为翅片高度，mm；T
g
为进出口烟气温度平均值，℃；T
f
为翅片处温度，℃；W
g
为烟气流量，kg/s；c
pg
为烟气侧比热，J/kg.℃；μ
g
为蒸发器水侧黏度，kg/m.s；k
g
为蒸发器水侧热导率，W/m.℃。4.根据权利要求3所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤1和2中，余热锅炉蒸发器变工况下的总体传热系数计算公式为：式中：U
evp
为变工况蒸发器中总体传热系数，W/m2.℃；A
t
为单位翅片长度下翅片的外表面面积，m2/m；A
i
为单位翅片长度下翅片的内表面面积，m2/m；ff
i
为翅片管内的污垢系数，m2℃/W；A
w
为单位翅片长度下翅片的内、外表面积平均值，m2/m；K
m
为翅片壁面导热系数，W/m2.℃；ff
o
为翅片管外的污垢系数，m2K/W；η为翅片的效率。
5.根据权利要求4所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤1和2中，基于余热锅炉蒸发器变工况下的总体换热系数计算得蒸发器的对数平均温差为：式中：LMTD为蒸发器的对数平均温差；e
w2
为蒸发器进口水焓，kJ/kg；e
w1
为蒸发器出口水焓，kJ/kg；A
evp
为蒸发器换热面积，m2。6.根据权利要求5所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤1和2中，蒸发器的对数平均温差还可通过余热锅炉蒸发器的进、出口烟气侧、水侧参数进行热力计算得到：式中：T
g1
为蒸发器进口烟气温度，℃；T
g2
为蒸发器出口烟气温度，℃；T
w1
为蒸发器进口水温度，℃。7.根据权利要求6所述的一种基于线性回归的余热锅炉蒸发器变工况建模方法，其特征在于，步骤3中，计算公式如下：将余热锅炉设计工况下的计算结果作为基准工况，即：U
b
＝f(W
g,b
,cp
g,b
,v
g,b
,km
g,b
,W
w,b
,cp
w,b
,v
w,b
,km
w,b
)式中：...

【专利技术属性】
技术研发人员：夏林，肖俊峰，胡孟起，高松，连小龙，李晓丰，
申请(专利权)人：西安热工研究院有限公司，
类型：发明
国别省市：