一种基于功率值计算的汽轮机流量曲线辨识及优化方法技术领域
本发明涉及汽轮机发电技术领域，尤其涉及一种基于功率值计算的汽轮机流量曲线
辨识及优化方法。
背景技术
汽轮机是现代火力发电厂中应用最广泛的原动机，其调门阀组的流量特性曲线是根
据原始设计而定，然而许多在役的老机组，尤其是经过DEH改造、汽机通流部分改造或
者经过大修的机组，经常发生调门阀组流量特性设定曲线与实际曲线不一致的情况。汽
轮机调门阀组作为DEH系统的主要执行机构，其流量特性偏差过大会导致汽轮机节流损
失加大、一次调频的响应负荷不足或者过大、AGC响应变慢、阀门切换负荷波动等问题，
最终影响机组的安全稳定运行。
在汽轮机机组的运行过程中，DCS系统每日存储大量历史数据，这些数据来源于运
行实际，是对机组运行特性的全面、客观反映，不可改变，具备绝对真实性。如何正确
运用历史数据，提取机组运行特性，使其更好地为运行实际服务，也是摆在工程技术人
员面前的一个课题。
发明内容
本发明的目的在于提供一种基于功率值计算的汽轮机流量曲线辨识及优化方法，能
够辨识并优化汽轮机调门阀组在特定主汽压力下的实际流量特性曲线，使调门阀组流量
特性设定曲线与实际曲线一致，保证汽轮机机组的安全稳定运行。
为实现上述目的，本发明采用的技术方案是：
一种基于功率值计算的汽轮机流量曲线辨识及优化方法，包括以下步骤：A、在汽轮机
机组设置若干个采样测点，并设定采样周期和采样时长，采集汽轮机机组的实时工况数
据，分析每个采样测点的实时工况数据包含的参数类型是否完整，如果参数类型有缺失
则重新采集，如果参数类型无缺失则进入步骤B；
B、在同一采样周期内，若一个或多个采样测点的实时工况数据中参数大小为零，或者
超出受其物理意义约束的数值范围，则清除该采样周期内所有采样测点的实时工况数据，
得到基本采样数据，然后进入步骤C；
C、从基本采样数据中选取输出功率值、主汽压力值和总阀位值均连续稳定若干个采样
周期的实时工况数据，并对实时工况数据中各个类型的参数进行均值化处理，得到稳定
采样数据，然后进入步骤D；
D、根据得到的稳定采样数据，采用聚类分析方法中的K-中心点算法挖掘出表征汽轮机
机组运行特性的目标数据集，得到采样工况数据，然后进入步骤E；
E、通过对汽轮机机组的输出功率方程式求全微分，计算在总阀位值不变的情况下，将
机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量，然后根据采样工况
下机组的输出功率值计算修正到特定主汽压力下时机组的输出功率值，从而辨识出特定
主汽压力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲线，
然后进入步骤F；
F、提取特定主汽压力值为额定主汽压力的汽轮机调门阀组流量特性曲线，根据机组负
荷调整的快速性和稳定性对机组流量特性曲线的要求，分析额定主汽压力下汽轮机调门
阀组流量特性曲线的线性度，然后采用分段线性方法对线性度不合理的区段进行优化，
得到优化后的机组总阀位值与输出功率值以及各阀阀位值之间的关系。
所述的步骤E中汽轮机机组的输出功率方程式为：
P i = DΔh t m a c η i 3.6 ]]>
式中Pi、D分别为采样工况下机组的输出功率标幺值和机组进汽量标幺值，ηi为全机相
对内效率，为机组理想比焓降，通过对汽轮机机组的输出功率方程式求全微分，
得到机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量ΔPi
ΔP i = Δh t m a c η i 3.6 ∂ D ∂ p 0 Δp 0 + Dη i 3.6 ∂ Δh t m a c ∂ p 0 Δp 0 + DΔh t m a c 3.6 ∂ η i ∂ p 0 Δp 0 ]]>
式中Δp0＝p'0-p0，p'0为修正后特定主汽压力，p0为采样工况的主汽压力，当汽轮机机
组的主汽压力变化范围在8MPa以内时，全机相对内效率ηi视为不变，即偏导数
汽轮机机组的总阀位值不变时，凝汽式机组或调节级为临界工况的机组进
气量与主汽压力成正比，即偏导数将水蒸汽视为理想气体，对于非再热整机
机组有
Δh t m a c = κ κ - 1 RT 0 [ 1 - ( p z p 0 ) κ - 1 κ ] ]]>
∂ Δh t m a c ∂ p 0 = RT 0 ( p z p 0 ) κ - 1 κ 1 p 0 = p 0 υ 0 ( p z p 0 ) κ - 1 κ ]]>
式中R为理想气体参数，k为水蒸气绝热系数，T0、v0分别为采样工况的温度和比体积，
pz为汽轮机排汽压力，将三个偏导数代入全微分方程，得到
ΔP i P i = [ 1 + p 0 υ 0 Δh t m a c ( p z p 0 ) κ - 1 κ ] Δp 0 p 0 ]]>
故将机组的采样工况修正到特定主汽压力p'0下时，机组的输出功率值P'i为
P i ′ = P i ( 1 + [ 1 + p 0 υ 0 Δh t m a c ( p z p 0 ) κ - 1 κ ] Δp 0 p 0 ) ]]>
所述的步骤F中提取辨识得到额定主汽压力下的汽轮机调门阀组流量特性曲线，进行高
阶次多项式拟合和线性化分析，然后采用分段线性方法对汽轮机调门阀组流量特性曲线
线性度不合理的区段进行优化时，根据汽轮机实际运行中机组通流能力与机组输出负荷
的唯一确定关系，固定机组输出功率值与各阀阀位值的对应关系，然后对汽轮机调门阀
组流量特性曲线线性度不合理的区段进行线性优化处理。
所述的步骤D中将稳定采样数据按照数值距离分配到相应的簇，目标函数采用平方误差
准则，采用K-中心点算法挖掘出包含150个以上采样工况的目标数据集。
所述的步骤C中输出功率值、主汽压力值和总阀位值连续稳定的采样周期为5个以上。
本发明通过采集并挖掘汽轮机机组运行的历史数据，利用基于功率值计算的方法辨
识出汽轮机调门阀组在特定主汽压力下的实际流量特性曲线，并采用分段线性方法对实
际流量特性曲线线性度不合理的区段进行优化，调整机组总阀位值与输出功率值以及各
阀阀位值之间的关系，使优化后的汽轮机调门阀组流量特性曲线与设定曲线一致，改善
机组运行调节品质。
附图说明
图1为本发明的流程示意图；
图2为采用本发明方法辨识得到的汽轮机调门阀组流量特性曲线；
图3为分段线性优化后的汽轮机调门阀组流量特性曲线；
图4为分段线性优化后的总阀位值与各阀阀位值对应关系。
具体实施方式
如图1所示，本发明所述的基于功率值计算的汽轮机流量曲线辨识及优化方法，包
括以下步骤：
A、在汽轮机机组设置若干个采样测点，并设定采样周期和采样时长，采集汽轮机
机组的实时工况数据，分析每个采样测点的实时工况数据包含的参数类型是否完整，如
果参数类型有缺失则重新采集，如果参数类型无缺失则进入步骤B。
B、针对所有采样测点，若在同一采样周期内，一个或多个采样测点的实时工况数
据中有参数大小为零，即出现空缺工况数据，或参数大小超出受其物理意义约束的数值
范围，即出现异常工况数据时，剔除该采样周期内所有采样测点的实时工况数据，得到
基本采样数据，分析基本采样数据中各个类型参数的最大值和最小值，以确定各个类型
参数的值域范围，然后进入步骤C。
C、根据汽轮机机组的控制逻辑，从基本采样数据中选取输出功率值、主汽压力值
和总阀位值均连续稳定5个以上采样周期的实时工况数据，亦即输出功率值、主汽压力
值和总阀位值在5个以上连续采样周期内的变化偏差均小于各自的规定精度，然后对实
时工况数据中各个类型的参数进行均值化处理，得到稳定采样数据，然后进入步骤D。
D、根据得到的稳定采样数据，采用聚类分析方法中的K-中心点算法，将稳定采样
数据按照数值距离分配到相应的簇，目标函数采用平方误差准则，挖掘出包含150个以
上采样工况的目标数据集，得到能够表征汽轮机机组运行特性的采样工况数据，然后进
入步骤E。
E、通过对汽轮机机组的输出功率方程式求全微分，计算在总阀位值不变的情况下，
将机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量，然后根据采样工
况下机组的输出功率值计算修正到特定主汽压力下时机组的输出功率值，从而辨识出特
定主汽压力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲
线。
汽轮机机组的输出功率方程式为：
P i = DΔh t m a c η i 3.6 ]]>
式中Pi、D分别为采样工况下机组的输出功率标幺值和机组进汽量标幺值，ηi为全机相
对内效率，为机组理想比焓降。通过对汽轮机机组的输出功率方程式求全微分，
得到机组的采样工况修正到特定主汽压力下时机组输出功率值的改变量ΔPi
ΔP i = Δh t mac η i 3.6 ∂ D ∂ p 0 Δp 0 + Dη i 3.6 ∂ Δh t mac ∂ p 0 Δp 0 + DΔh t mac 3.6 ∂ η i ∂ p 0 Δp 0 ]]>
式中Δp0＝p'0-p0，p'0为修正后特定主汽压力，p0为采样工况的主汽压力。
当汽轮机机组的主汽压力变化范围在8MPa以内时，全机相对内效率ηi视为不变，
即偏导数汽轮机机组的总阀位值不变时，凝汽式机组或调节级为临界工况
的机组进气量与主汽压力成正比，即偏导数将水蒸汽视为理想气体，对于非
再热整机机组有
Δh t m a c = κ κ - 1 RT 0 [ 1 - ( p z p 0 ) κ - 1 κ ] ]]>
∂ Δh t m a c ∂ p 0 = RT 0 ( p z p 0 ) κ - 1 κ 1 p 0 = p 0 υ 0 ( p z p 0 ) κ - 1 κ ]]>
式中R为理想气体参数，k为水蒸气绝热系数，Y0、v0分别为采样工况的温度和比体积，
pz为汽轮机排汽压力，将三个偏导数代入全微分方程，得到
ΔP i P i = [ 1 + p 0 υ 0 Δh t m a c ( p z p 0 ) κ - 1 κ ] Δp 0 p 0 ]]>
故将机组的采样工况修正到特定主汽压力p'0下时，机组的输出功率值P'i为
P i ′ = P i ( 1 + [ 1 + p 0 υ 0 Δh t m a c ( p z p 0 ) κ - 1 κ ] Δp 0 p 0 ) ]]>
由此得到特定主汽压力下不同总阀位值对应的输出功率值，从而辨识出特定主汽压
力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲线，然后
进入步骤F。
F、提取特定主汽压力值为额定主汽压力的汽轮机调门阀组流量特性曲线，进行高
阶次多项式拟合和线性化分析，然后根据汽轮机机组调节的快速性和稳定性要求，以及
参与调频时总阀位值叠加相应的调频分量，机组能够贡献相应的调频负荷的要求，分析
额定主汽压力下汽轮机调门阀组流量特性曲线的线性度，然后采用分段线性方法对线性
度不合理的区段进行优化。
根据汽轮机实际运行中机组通流能力与机组输出负荷的唯一确定关系，固定机组输
出功率值与各阀阀位值的对应关系，然后线性化汽轮机调门阀组流量特性曲线线性度不
合理区段总阀位值与输出功率值的关系，得到优化后的机组总阀位值与输出功率值以及
各阀阀位值之间的关系，即优化后的分配系数和阀位管理函数。
下面结合具体实施例说明本发明所述的基于功率值计算的汽轮机流量曲线辨识及
优化方法。在荥阳某电厂超临界、单轴、三缸四排汽、一次中间再热、双背压凝汽式汽
轮机机组设置13个采样测点，采集到该机组2013年07月9日至2013年7月13日的
120个数据文件，每个数据文件的采样时间为1小时，采样周期为5秒，剔除空缺及异
常工况数据，整合所有数据文件后得到86400个基本采样数据，从中选出输出功率值、
主汽压力值和总阀位值均连续稳定5个采样周期的实时工况数据，并对各个类型的参数
进行均值化处理，得到稳定采样数据，采用聚类分析方法中的K-中心点算法从稳定采样
数据中挖掘出包含150个采样工况的目标数据集，得到能够表征汽轮机机组运行特性的
采样工况数据。通过对汽轮机机组的输出功率方程式求全微分，计算在总阀位值不变的
情况下，将机组的采样工况修正到额定主汽压力24.2MPa下时机组输出功率值的改变量，
然后根据采样工况下机组的输出功率值计算修正到额定主汽压力24.2MPa下时机组的
输出功率值，从而辨识出额定主汽压力24.2MPa下机组总阀位值与输出功率值之间的关
系，图2所示即为额定主汽压力24.2MPa下的汽轮机调门阀组流量特性曲线。
根据汽轮机机组调节的快速性和稳定性要求，以及参与调频时总阀位值叠加相应的
调频分量，机组能够贡献相应的调频负荷的要求，分析图2所示的汽轮机调门阀组流量
特性曲线，可知总阀位值在80.4％以下的区段出现太过分散的块状，曲线整体线性度较
好，但是斜率过小，负荷调节响应缓慢；总阀位值在80.4～82.35％、89.2～93.31％之间
的区段明显呈非线性，顺序阀对应前后阀重叠度设置不合理，负荷调节时易引起阀门晃
振、影响一次调频响应；总阀位值在82.35～89.2％之间的区段线性度较好，但斜率过大，
机组调整可以满足快速性的要求，但稳定性不好，容易引起负荷抖动；总阀位值93.31％
以上的区段线性度较好，能满足机组调节的快速性和稳定性要求。
采用分段线性方法对线性度不合理的区段进行优化，在给定的机组总阀位值运行行
程内，以总阀位值89％为分段点，根据汽轮机实际运行中机组通流能力与机组输出负荷
的唯一确定关系，固定机组输出功率值与各阀阀位值CV1/至CV4的对应关系，将汽轮
机调门阀组流量特性曲线分为两段进行线性优化，图3和图4分别为优化后的汽轮机调
门阀组流量特性曲线以及机组总阀位值与各阀阀位值CV1/至CV4的对应关系，能够较
好地满足机组调节的快速性和稳定性要求。
本发明有效利用汽轮机机组在实际运行过程中存储的大量历史数据，基于功率值计
算得到将机组的采样工况修正到特定主汽压力下的输出功率值，从而辨识出特定主汽压
力下机组总阀位值与输出功率值之间的关系，即为汽轮机调门阀组流量特性曲线，最后
通过分段线性方法优化曲线线性度不合理的区段，得到优化后的机组总阀位值与输出功
率值以及各阀阀位值之间的关系，使得汽轮机调门阀组流量特性设定曲线与实际曲线一
致，满足机组调节的快速性和稳定性要求。